[发明专利]起伏地表的地震叠后正演方法

权利要求书

1、起伏地表的地震叠后正演方法，其特征在于采用以下步骤：

(1)采用矩形框对目标地质模型进行整体封装，矩形框的顶面即为输出的基准面，矩形顶面位于起伏地表的最高点之上，起伏地表上方采用任意速度值或密度值进行充填，并采用规则网格对速度模型和密度模型进行网格化，如果只有速度模型，则采用密度与速度的经验公式计算密度模型；

(2)采用斯奈尔定理计算反射系数，并利用雷克子波进行褶积处理，得到反射系数模型，雷克子波的主频根据工区的地层特点进行选择，一般介于10Hz～90Hz之间；

(3)将起伏地表的反射系数赋值为0，消除起伏地表的反射影响，并对反射系数数组进行正傅立叶变换，变换到频率域作为初始延拓波场；

(4)利用单程波动方程的傅立叶有限差分延拓算子沿着深度轴从下往上进行波场延拓，当延拓深度z_i介于输入基准面与输出基准面之间时，在深度z_i+1与z_i之间延拓的波场增量值赋值为零，即保持输入基准面和输出基准面之间的波场延拓增量为零，根据工区的地层特点，参考初始波和终止波的时间，确定波场延拓的频率范围f1和f2，f1一般介于1Hz～5Hz之间，f2一般介于50Hz～100Hz之间，对其它频率的波场值充零处理，

(5)将延拓后的频率域的波场进行反傅立叶变换，拾取时间域波场的实部，得到最后的成像结果；

(6)通过常规的显示软件将成像结果处理为地震剖面图像。

2、根据权利要求1所述的起伏地表的地震叠后正演方法，其特征在于：对目标模型进行矩形封装，矩形顶面位于起伏地表的最高点之上；并通过在输入基准面和输出基准面之间设置一个零波场层来参与正常的波场传播和延拓，即当延拓深度z_i介于输入基准面与输出基准面之间时，在深度z_i+1与z_i之间延拓的波场增量值赋值为零。

3、根据权利要求1所述的起伏地表的地震叠后正演方法，其特征在于：采用斯奈尔定理计算反射系数，并利用雷克子波进行褶积处理，得到反射系数模型，雷克子波的主频根据目标工区的地层特征选定，一般介于10Hz～90Hz之间。

4、根据权利要求1或2所述的起伏地表的地震叠后正演方法，其特征在于：采用频率域的傅立叶有限差分方法算子沿深度轴从下往上进行波场延拓计算。