[发明专利]基于起伏地表的井间地震波动方程叠前深度偏移成像方法

权利要求书

1.一种基于起伏地表的井间地震波动方程叠前深度偏移成像方法，其特 征在于，该方法包括以下步骤：

A、预设基准面，将反射波场进行初始化后得到初始反射波场；

B、对所述初始反射波场使用带误差补偿的频率-空间域有限差分波场延 拓算子计算波场值，所述带误差补偿的频率-空间域有限差分波场延拓算子具 体为：

∂u~∂z=±jωv(x,y,z)[Q+E]u~,]]>

其中：为频率-空间域反射波场；v(x，y，z)为介质速度；ω为频率；j为虚数 单位；±符号分别对应频率-空间域下行反射波场和上行反射波场； Q=1+Σi=1nαiRx1+βiRx+Σi=1nαiRy1+βiRy,]]>其中，Rx=ν2(x,y,z)ω2∂2∂x2,]]>Ry=ν2(x,y,z)ω2∂2∂y2,]]>α_i、 β_i是连分式展开系数；E=1-(vkxω)2-(vkyω)2-[1-α(vkxω)21-β(vkxω)2-α(vkyω)21-β(vkyω)2],]]>k_x、 k_y分别是频率-波数域x和y方向的波数；

C、在某一延拓步长上，如果检测到炮点或检波器，则计算出炮点或者 检波器对应的横向位置的波场值，所述计算出炮点或者检波器对应的横向位 置的波场值具体为：将步骤B计算出的所述波场值加上所述某一延拓步长位 置处的检波器对应的地震道上得到的波场值或者炮点对应的地震道上得到的 波场值；

D、对成像空间范围内所有的炮点或检波器重复步骤C，得到所有炮点 或者检波器对应的横向位置的波场值；

E、对步骤D得到的所有波场值都进行处理后进行叠加成像，具体为：

E1、根据相关成像条件提取每个波场值的成像值，并将每个成像值放置 到对应的横向位置的波场中，得到成像剖面图；

E2、对每个炮点或者检波器重复步骤A、B、C、D及E1直至得到所有 炮点或者检波器的多个成像剖面图；

E3、对所述多个成像剖面图进行叠加成像。