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

权利要求书

1.一种绕射波成像方法，其特征在于，包括：

读取待成像的地震数据，并根据地震数据计算倾角场及观察倾角；

根据倾角场及观察倾角构建反极性反稳相滤波器；

基于反极性反稳相滤波器得到边缘绕射波成像模型；

根据所述边缘绕射波成像模型及地震数据计算绕射波成像图像；

所述根据倾角场及观察倾角构建反极性反稳相滤波器的步骤包括：

所述反极性反稳相滤波器表述为：

式中，s为炮点位置，x为成像点位置，r为检波点位置，θ₀(x)为倾角场，θ(s,x,r)为观测倾角。

2.根据权利要求1所述的绕射波成像方法，其特征在于，所述读取待成像的地震数据为U(t,s,r)，该地震数据为采用均值滤波方法对初始地震数据U(t',s',r')进行噪音压制后得到的，式中，t为入射射线与绕射射线的射线双程旅行时间，s为炮点位置，r为检波点位置。

3.根据权利要求1所述的绕射波成像方法，其特征在于，所述倾角场θ₀(x)的计算方法包括：

根据地震数据计算偏移剖面，根据偏移剖面计算倾角场θ₀(x)。

4.根据权利要求1所述的绕射波成像方法，其特征在于，所述观测倾角θ(s,x,r)的计算方法为：

式中，Ps表示炮点位置到成像点位置的矢量，Pr代表检波点位置到成像点位置的矢量，z代表垂直向下的单位矢量。

5.根据权利要求1所述的绕射波成像方法，其特征在于，基于反极性反稳相滤波器w_dp(s,x,r)得到的边缘绕射波成像模型为：

I(x)＝∫∫∫w_dp(s,x,r)U(t,s,r)δ(t-t_d(s,x,r))dtdrds；

式中，δ(t-t_d(s,x,r)为狄拉克函数，U(t,s,r)为地震数据，t为入射射线与绕射射线的射线双程旅行时间，t_d(s,x,r)代表射线自炮点位置s到检波点位置r再到成像点位置x经过的成像双程旅行时间。

6.根据权利要求5所述的绕射波成像方法，其特征在于，所述成像双程旅行时间t_d(s,x,r)的计算方法为：

式中(x_s,y_s),(x₀,y₀),(x_r,y_r)分别代表炮点位置s，成像点位置x，检波点位置r的地表坐标。

7.一种绕射波成像装置，其特征在于，包括：

数据读取模块，用于读取待成像的地震数据，并根据地震数据计算倾角场及观察倾角；

滤波器构建模块，用于根据倾角场及观察倾角构建反极性反稳相滤波器；

绕射波模型模块，用于基于反极性反稳相滤波器得到边缘绕射波成像模型；

绕射波成像模块，用于根据所述边缘绕射波成像模型及地震数据计算绕射波成像图像；

所述滤波器构建模块中的所述反极性反稳相滤波器表述为：

式中，s为炮点位置，x为成像点位置，r为检波点位置，θ₀(x)为倾角场，θ(s,x,r)为观测倾角。

8.一种电子设备，其特征在于，包括处理器和存储器，所述处理器和所述存储器连接；

所述存储器用于存储程序；

所述处理器用于运行存储在所述存储器中的程序，执行如权利要求1-5中任一项所述的方法。

9.一种存储介质，其特征在于，其上存储有计算机程序，所述计算机程序在被计算机运行时执行如权利要求1-5中任一项所述的方法。