波动方程叠前深度偏移原理要点.docVIP

3 波动方程法叠前深度偏移 常用的叠前深度偏移方法包括射线法和波动方程法。射线法主要指基于绕射旅行时计算的Kirchhoff积分法，在绕射旅行时计算方法上可以采用基于函程方程的变速射线追踪法、基于费马原理的二维有限差分法和稳健高效的三维迎风有限差分法；而波动方程叠前深度偏移是复杂介质成像的有效手段，能够解决强横向变速条件下复杂地质体的地震波成像问题。基于共炮集的波动方程叠前深度偏移的基本思路是，首先对每一炮进行单炮偏移成像，然后再把各炮成像结果在对应地下位置上叠加，从而得到整个剖面成像。从计算角度而言，成像过程是很简单的步骤，波场外推算子决定了偏移方法的效率、成像精度及其适应范围。一般要求偏移算子能够适应陡倾角反射的成像及剧烈的横向速度变化，同时具有较高的计算效率。 波动方程叠前深度偏移的基本思路 基于共炮集的波动方程叠前深度偏移的基本思路是，首先对每一炮进行单炮偏移成像，然后再把各炮成像结果在对应地下位置上叠加，从而得到整个成像剖面。对于每一炮，标准的波动方程叠前深度偏移可以分为三步：震源波场的正向延拓、炮集记录波场的反向延拓和应用成像条件求取成像值（Clearbout, 1971）。 为了方便叙述基于共炮集的波动方程叠前深度偏移的基本过程，我们引入基于单程波方程的波场传播算子（Berkhout, 1987），并以频率域二维波场为例加以说明。对震源波场和炮集记录波场做如下定义： （1）：它是炮点处频谱为的点源激发产生的震源波场，有 （4-1） （2）：它是点处激发，排列接收到的记录波场，该波场可以写成： （4-2） 其中，含有一非零道，即在接收点处的记录道，它满足： （4-3） （3）：它表示在深度处的正向延拓波场，如果引入表征波场从地面传播到深度的传播算子，则有： （4-4） （4）：它表示记录波场在深度的反向延拓波场： （4-5） 其中，为记录波场的反向传播算子。因为波场传播算子描述上行波从深度到地面的传播过程，故描述了（向上传播的）记录波场从地面到深度的反向延拓过程（Berkhout, 1987）。和分别称为下行波和上行波的深度外推算子。实际计算过程中，逐层实现上、下行波的波场延拓和求取成像值。层内波场延拓如图4-1所示。 图4-1 叠前深度偏移波场延拓示意图 其中，波场深度外推算子和分别表征下行波从深度到深度的正向延拓过程和上行波从深度到深度的反向延拓过程。 注意到所有波场空间上离散分布在采样间隔为的地震道上，故空间函数可以用长度为、振幅为的加窗函数表示，而且以上公式中的积分可以用离散求和替代。 用传统的偏移公式，我们可以得到对点处震源在处的单一记录道的叠前深度偏移结果具有如下形式： （4-6） 式中表示复共轭，表示取复数的实部。这种频率域的成像公式相当于时间域震源波场正向延拓值与记录波场反向延拓值的互相关。由（4-6）式，可以得到共炮集数据的叠前深度偏移成像公式： （4-7） 如果采用的是有限差分算法进行波场延拓计算，则上式中关于的求和并非是显式的，它一般在对整个单炮记录波场的反向延拓中自动实现。因而，共炮集记录叠前深度偏移公式可表述为： （4-8） 其中，为炮整个记录波场的反向延拓波场值。 从计算角度而言，成像过程是很简单的步骤，波场外推算子的数学形式和计算实现才是地震波偏移成像的核心。目前常用的地震成像的方法主要有：Kirchhhoff积分法偏移技术、有限差分偏移技术、域偏移技术和域偏移技术等。但各种方法都有其适用范围和限制条件：对于起伏地表条件下的叠前深度偏移，Kirchhoff积分法可以灵活地处理起伏的地表条件，但它需要利用射线追踪技术计算偏移成像所需要的旅行时，这样在处理地下复杂构造时，多路径问题经常使射线追踪不能获得准确的旅行时，导致对复杂构造成像的精度较低，效果不理想；有限差分法直接求解函数方程的Kirchhoff积分偏移方法，一般也仅能计算初至旅行时，无法处理在复杂地质体的成像效果。基于波动方程波场外推理论的方法能较好的解决多路径问题，使成像结果更精确。目前所有基于波动方程的波场延拓算子有波动方程有限差分波场延拓算子和