[发明专利]绕射波成像方法、装置和电子设备

说明书

本发明提供了一种绕射波成像方法、装置和电子设备，其中，该方法包括：获取角度域下的共成像道集；共成像道集包括多个成像道集；每个成像道集对应一个成像点；每个成像道集包括多个采样点；每个采样点包括时间参数；对每个成像道集进行下述处理：按照时间参数对采样点进行划分得到多个时间切片，计算每个时间切片的每一采样点的衰减因子与衰减系数，对每个时间切片的若干采样点及若干采样点对应的衰减系数进行叠加计算，得到每个时间切片的成像结果，将其进行集合，得到每个成像点的成像结果，将每个成像点的成像结果进行集合，得到共成像道集的成像结果。本发明有效地压制了反射波和噪声，增强了绕射波，提高了绕射波成像结果的分辨率。

技术领域

本发明涉及地震勘探技术领域，尤其是涉及一种绕射波成像方法、装置和电子设备。

背景技术

煤炭开采过程中，突水、瓦斯突出、顶板塌陷等地质灾害严重威胁煤矿安全，这些地质灾害常与一部分地质异常体有关，地震波在传播过程中，遇到地质异常体时会形成新震源，以球面波形式继续向前传播，形成绕射波或散射波，绕射波携带有地质异常体的有效信息，但是由于反射波能量较强，绕射波成像结果会受到反射波的影响，因此有效分离绕射波可以有效提高绕射波成像结果。

现有平面波分解方法通过构建平面波解构滤波器对反射波进行压制，增强绕射波来达到绕射波分离的目的，但是，该方法的精度严重依赖于反射波局部倾角值的预测，且分离后的绕射波会包含噪声等干扰，影响绕射波成像结果的分辨率。

发明内容

有鉴于此，本发明的目的在于提供一种绕射波成像方法、装置和电子设备，以压制反射波和噪声，增强绕射波，提高绕射波成像结果的分辨率。

第一方面，本发明实施例提供了一种绕射波成像方法，其中，该方法包括：获取角度域下的共成像道集；共成像道集包括多个成像道集；每个成像道集对应一个成像点，每个成像道集包括多个采样点；采样点包括时间参数；对每个成像道集执行以下处理步骤：根据时间参数对采样点进行时间划分，得到多个时间切片；时间参数相同的若干采样点的集合为一个时间切片；对每个时间切片的每一采样点进行第一计算，得到每一采样点对应的衰减因子；对衰减因子进行第二计算，得到衰减系数；对每个时间切片的若干采样点和若干采样点对应的衰减系数进行叠加计算，得到每个时间切片的成像结果；将成像道集对应相关的多个时间切片的成像结果进行集合，得到成像道集对应的成像点的成像结果；通过对共成像道集的每个成像道集进行以上处理步骤，得到多个成像道集对应的多个成像点的成像结果；将共成像道集的多个成像点的成像结果进行集合，得到共成像道集的成像结果。

结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，采样点还包括倾角参数、方位角参数；成像道集为cig_num(dip,azi,t)；其中，dip为倾角参数，azi为方位角参数，t为时间参数；成像道集包括多个采样点x_(i,j,k)；其中，k为时间参数t的一个取值，(i,j)分布表示时间参数k对应的时间切片的采样点x的分布，i的取值为(1,n)，j的取值为(1,m)，n、m分别为大于0的自然数。

结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第二种可能的实施方式，其中，该方法还包括：利用下述公式计算衰减因子：其中，x_(i,j,k)表示采样点；u表示时间参数k对应的时间切片的所有采样点的均值；σ表示时间参数k对应的时间切片的所有采样点的标准差；dm(x_(i,j,k))表示采样点x_(i,j,k)对应的衰减因子。

结合第一方面的第二种可能的实施方式，本发明实施例提供了第一方面的第三种可能的实施方式，其中，对衰减因子进行第二计算，得到衰减系数的步骤，包括：利用下述公式计算衰减系数：