[发明专利]一种起伏地表条件下的叠前深度偏移方法

权利要求书

1.一种起伏地表条件下的叠前深度偏移方法，其特征在于，所述的方法包括如下步骤：

步骤1：以最高高程为基准面z_beg=0，在基准面与起伏地表之间填充近地表速度；

步骤2：假设在最初基准面z_beg=0时初始波场u(x,z_beg)=0，记录波场为u_e(x,z)，将z_beg=z=0时对应的波场加入到初始波场u(x,z_beg)=u_e(x,z)；

步骤3：在单炮偏移孔径范围内将波场u(x,z)划分成多个小波束窗，将波场u(x,z)变换到小波束域，分别将下行波场和上行波场沿深度方向延拓Δz至深度z+Δz，并变回到空间域，延拓后的下行波场为和延拓后的上行波场为

步骤4：对波场下行波场和上行波场进行滤波处理；

步骤5：根据互相关成像条件对所述下行波场和上行波场进行互相关成像，得到深度层z+Δz的成像结果；

步骤6：判断深度z+Δz层是否存在波场u_e(x,z+Δz)，如果存在，将相应波场叠加到所述延拓波场中，生成深度z+Δz层的开始延拓波场u(x,z+Δz)=u_con(x,z+Δz)+u_e(x,z+Δz)；

步骤7：重复所述步骤3至步骤6，逐步延拓并成像直至最大深度层z_max，得到单炮的成像结果，最终的偏移结果为所有单炮成像结果的叠加。

2.根据权利要求1所述的起伏地表条件下的叠前深度偏移方法，其特征在于，所述的步骤3包括：

将每一深度层z处的波场分解为：

u(x,z,w)=ΣmΣnu^(x‾n,ξ‾m;z)bmn(x)---(1)]]>

将（1）式代入频率-空间域波动方程[∂x2+∂y2+w2v2(x,y)]u(x,z,w)=0]]>中得到标量声波方程ΣnΣmu^(x‾n,ξ‾m;z)[∂x2+∂z2+w2v2(x,z)]bmn(x)=0;]]>

将小波束的演化表示为amn(x)=e±iAnΔzbmn(x);]]>

根据波场延拓算子A_n将所述波场u(x,z)延拓至下一深度层z+Δz，生成延拓后的波场u_con(x,z+Δz)；

其中，w是频率；v(x,y)代表层速度；u(x,z,w)是地表(x,z)处接收的地震记录；b_mn(x)代表空间位置波数为的小波束；是一系列关于深度z的系数，a_mn是小波束b_mn经过介质传播后的波束。