地球物理资料数字处理(第二十五讲).docVIP

第二节 有限差分法波动方程偏移 （Finite Difference Wave Equation Migration） 波动方程偏移是建立在射线理论基础上的绕射扫描偏移具有明显的优越性。因此地震波在地下介质的传播遵循波动方程规律，而射线理论只是近似地描述地震波的传播规律。因此，应用波动方程进行偏移，一方面偏移归位效果好，特别是为复杂的地质条件下反射层的偏移提供了可能性；另一方面在实现反射层偏移归位的过程中保持了反射波的特征。 1979年，J．F．Claerbout提出了用有限差分法直接解简化的波动方程，以实现地震波场的空间归位。这种方法是对标量波动方程进行了差分，从波场延拓的观点获得地震偏移剖面。因为用了差分网格，可以在垂向和横向改变速度参数，适用于一般的非均匀介质。 这种方法是首先对波动方程进行Claerbout坐标变换，将标量波动方程分为上行波和下行波两部分，这样可以避免用标量波动方程做波场延拓所带来的不适定性问题；然后忽略深度方向的二次偏导数项，得到变换后的简化波动方程；最后用有限差分法求解波动方程，使发射层偏移归位到反射界面空间真实位置。 在这一节首先介绍波动方程偏移的成像原理，然后介绍叠后有限差分法波动方程偏移的原理和差分方程的建立及其求解。 波动方程偏移的成像原理（Imaging Principle of Wave Equation Migration ） 无论那种偏移方法，包括射线理论的和波动理论的偏移方法，抽象出来都不过是把地表地震剖面转换成近似的地质剖面的一个过程，它包括两个步骤：波场延拓和成像。 波场的延拓就是已地表地震剖面为边界条件，求取地下各个点的波场值，也称为波场重建。成像就是在得到了地下各个点的波场之后， 从重建波场中检测反映界面信息的波场值的过程。J。K CLAERBLOUT成像原理指出：存在于地壳里的反射界面，是处在下行波到达且刚刚出现上行波那一瞬间。 爆炸反射界面成像原理 爆炸反射界面成像原理是最常用，最简单的一种成像原理。它把地下反射界面想象成旧有爆炸性的物质后爆炸源的形状，位置与反射界面的形状和位置一致。它所产生的波动脉冲波，其强度，极性与界面的反射系数大小和正负一致，因此称地下反射界面为爆炸反射面。假定在t=0时刻，所有的爆炸反射界面同时爆炸，发射上行波到达地面各个观测点，波的传播速度为介质速度的一半，如果利用波动方程将地面测得到的波场记录作反时间方向传播，即向下延拓，则t=0时刻的波场值就正确地描述了地下反射界面的位置，自动实现了偏移成像。 此成像原理使用与对水平叠加地震剖面的偏移处理，因为水平叠加剖面相当于零炮检距剖面，自炮点发出的下行波到达反射点的路径与自该点反射返回地面的上行波的路径一致，由于做半速代换，所以叠后偏移方法就是利用延拓波场的t=0时刻的值成像。 2. 线形下延成像原理 现行下延成像原理是用于有炮检距的的记录，常在叠前偏移时使用。叠前偏移的基本思想就是以地表记录为边界条件。在整条射线上交替的把炮点和检波点向下沿拓，直到炮点和检波器相当于埋到地下同一位置为止。由于这是对地下每一点而言，就好像即埋置了震源又埋置了检波器一样，根据J .F.Claerbout 的原始成像原理，显然这时要取t=0时刻的波场值输入。也适用于零炮检距记录的偏移成像，但是对零炮检距记录使用爆炸反射界面成像原理更单。 3. 波场传播的时间一致性成像原理 波场传播时间的一致性成像原理是弹性波偏移中的一种成像原理。因为弹性波偏移同时要同时处理纵波、横波和转换波，因此它应有自己特殊的成像原理形式。在弹性波偏移中，一般只向下沿拓检波点而炮点不动，即沿拓后的波场好像是炮点仍在地表位置而地下各点都埋置了检波器接收到的波场一样。为了说明这种成像原理，以入射波P波反射纵波PP波和转换波PS为例。如图7-5-a所示，入射纵波传播时间为t反射纵波传播时间为t转换波传播时间为t。只有在反射点Q处才有t=t=t。即这时反射纵波PP和转换波PS在反射点Q上结合在一起 。如图7-5（b） 是弹性波的成像原理 其中 图中ｔ＝ ｔ＝ ｔ＝ 图７－５波场传播时间一致性成像原理示意图 由于沿拓到Ｒ点处时，如果Ｒ不在反射界面上，一般来说，ｔｔｔ，只有当Ｒ反射界面上时，才有ｔ＝ｔ＝ｔ。因此，在弹性波的沿拓波场中，必须取沿拓的纵波和横波（包括转换波）波场的传播时间与入射波传播时间一致的那一时间来构成偏移输出剖面。 二 15有限差分法波动方程的导出 (15 Finite Difference Wave Equation Derivation) Claerbout差分偏移是对水平叠加时间剖面，既自激自收地震记录进行时的，根据爆炸反射界面成像原理，认为水平叠加剖面是位于地下反射界面上的炮点零时刻产生的脉冲波，沿界面的法向以半速度上行到地