[发明专利]一种复杂海底上下缆地震采集数据边界元延拓校正方法

权利要求书

1.一种复杂海底上下缆地震采集数据边界元延拓校正方法，包括以下步骤：

1）构造波动方程基准面校正算子

2）根据复杂海底地形数据构造空间滤波算子filt(x,z)；

3）选择某个复杂海底最高点之上的水平面作为延拓初始面，校正基准面选择位于复杂海底最低点之下，对延拓初始面和复杂海底之间的海水层充填复杂海底速度；

4）自延拓初始面激发初始波场，并开始向下做深度步延拓，在每个延拓的层的高程面上判断是否有新波场加入，有则累加到初始延拓波场中，形成一新的延拓波场，以该波场继续向下延拓，否则以原波场继续向下延拓，利用空间滤波方法来校正上下缆采集的地震数据,在每一步延拓中滤除非真实海底空间上的波场，逐层延拓直到校正基准面完成整个道集的延拓过程。

2.如权利要求1所述的一种复杂海底上下缆地震采集数据边界元延拓校正方法，其特征在于：所述步骤1）中，构造波动方程基准面校正算子包括以下步骤：

①根据边界元波动方程理论，均匀各向同性弹性介质中的频域弹性波方程为：

μ▿2u(x,ω)+(λ+μ)▿▿·u(x,ω)+ρω2u(x,ω)=-f(x,ω)---(1)]]>

式中μ、λ为介质拉梅系数，ρ为介质密度，x为测线方向坐标，ω为圆频率，u(x,ω)为位移场，f(x,ω)为震源矢量；

②利用边界元波动方程得到的介质波场的频散关系式为：

u^(kx,z)=FTX[u(x,z)exp(ik0Δz(n(x)-1))],---(2)]]>

式中FT_x为从x→k_x的正向傅立叶变换，k_x为测线方向坐标x处对应的波数值，k₀为深度z处对应的波数值，n(x)为测线方向坐标x处对应的折射率，z为波场值的深度，Δz为波场延拓步长的深度间隔；

③利用频散关系确定分裂步算子的宽带特性，对应的波场延拓算子为：

P~(x,z+Δz,ω)=IFTx{exp(ik0Δz)FTx[P(x,z,ω]}---(3)]]>

式中为深度z+Δz处对应的波场，[P(x,z,ω]为深度z处对应的波场。