第3章Kirchhoff积分法叠前深度偏移

第三章 Kirchhoff积分法叠前深度偏移大家知道，叠前偏移的概念早在70年代中期就提出来了，但由于叠前记录的信噪比较低，偏移的初始模型又很难选准，加之当时的计算机无法承受叠前偏移较大的计算量，直到90年代叠前偏移才开始尝试应用于油气勘探地震数据的精细处理中。常见的叠前深度偏移方法可以分为两类：第一类是基于绕射扫描叠加原理的Kirchhoff积分法，另一类是基于波动方程的偏移方法（如有限差分偏移方法、Fourier偏移方法等）。本章重点讨论Kirchhoff积分法叠前深度偏移。Kirchhoff积分法叠前深度偏移被认为是一种高效实用的叠前深度偏移方法，目前主要完善三维采集和叠前深度偏移软件。积分法具有高偏移角度、无频散、占用资源少和实现效率高的特点，并且积分法能够适应变化的观测系统和起伏的地表，优化的射线追踪法和改进的有限差分法能够在速度场变化的情况下快速准确地计算绕射波和反射波旅行时，从而使积分法能够适应复杂的构造成像。地震偏移成像问题，经过最近十多年的研究与发展，已经基本解决了和正在解决三维偏移，叠前深度偏移和多分量地震偏移等诸问题。但是偏移中有诸多问题尚未解决，例如真振幅偏移问题和各向异性介质中的地震偏移问题。近年来，解决真振幅偏移问题就是偏移地震数据得到真正的振幅和相位信息，从而为岩性解释服务。由于积分法具有许多优点，因此研究Kirchhoff型保幅叠前深度偏移具有很高的理论价值和实用价值。下面就变速射线追踪法计算走时、有限差分法计算走时以及Kirchhoff型常规叠前深度偏移和保幅叠前深度偏移做详细讨论和分析。3.1 变速射线追踪法计算走时Kirchhoff积分法叠前深度偏移已在实际生产中应用了多年，并解决了不少复杂构造的成像问题（Zhu & Lines, 1998）。Kirchhoff积分法的关键是绕射旅行时的计算，目前常用的计算方法是射线追踪法和有限差分法（Schneider, 1992, 1995）。有限差分绕射旅行时计算是基于费马原理，可在直角坐标系或球坐标系实现，具体方法原理将在本章第二节介绍。射线追踪法计算绕射旅行时可分为常速法和变速法，常速法很简单，在此不再赘述；下面主要介绍变速法。考虑到地下介质在纵横向上通常是变速的，为通过射线追踪较精确地生成CMP道集中各道的反射旅行时，下面我们基于Langan（1985）的思想，推导变速介质条件下的射线追踪方程。 由程函方程可推出如下的射线方程 （3-1）其中，是波速，是空间位置，s是与路径长度有关的仿射参数。路径长度l由（3-2）式给出 （3-2）其中， （3-3）是位置处指向射线传播方向的矢量，射线旅行时为 （3-4）若为正实数且时不变，则有 （3-5）由射线方程，经推导得出 （3-6）其中，是炮点(s=0)相对于原点的位置，为s=0处射线方向的单位矢量。 当射线通过恒定速度梯度或恒定慢度梯度介质传播时，可建立方程（3-6）的简单近似解。如果射线通过恒定速度梯度介质传播，速度场可表示为 （3-7）其中，是原点处的速度，为速度梯度。借助于Taylor展开()，由（3-6）、（3-3）和（3-4）分别求得 （3-8） （3-9） （3-10）其中，, 是射线入射点处的速度。方程（3-8）、（3-9）、和（3-10）就是恒定速度梯度介质条件下的射线追踪方程。3.2 有限差分法计算走时一二维有限差分绕射旅行时计算方法在Kirchhoff叠前深度偏移中，绕射走时的计算精度直接影响了成像精度。为此，我们利用Eikonal方程、费马原理和有限差分近似，基于二维介质在规则网格上进行绕射走时的有限差分计算，它无需走时的内插，并且把它扩展到三维是很简单的。1方形网格情况下的绕射走时有限差分计算在高频近似下，二维几何路径及二维波前的传播符合射线追踪的程函方程， （3-11）方程（3-11）建立了走时梯度与速度模型的基本定量关系。坐标轴是x和z，s0是慢度（速度的倒数）。方程（3-11）中的两个微分项能够应用有限差分来近似表示： （3-12a） （3-12b）t0、t1S0t2t3其中t0、t1、t2和t3的含义如图3-1所示，和分别是横向和纵向采样间距，方形网格情况下，。在平面波近似下，已知t0、t1和t2可以求t3的值。即已知四个点中的任意三个都可以把另一点计算出来。 （3-13）为获取待求行或列上第一个点处的旅行时，利用有限差分得出如下计算公式 图3-1 方形网格走时计算示意图 （3-14）其中，分别是内行或内列上旅行时的相对极小值、两侧的旅行时、和待求行或列上的旅行时。2矩形网格情况下的绕射走时有限差分计算为提高绕射走时有限差分计算方法的适应性，考虑具有不同纵、横向采样间距的矩形网格。矩形网格走时计算见图3-2a。基于方程（3-11）和（3-12），经推导可以得出：（3-15）图3-2b 十字型中, O点是震源点的位置, 十字型各点的走时已知, 各顶点的走时可由相邻三个点的走时求出。OABCD图3-2a 矩形网格中水平距离为, 垂直距离为,慢度s0为常数, t0、t1和t2的走时已知, 可以求t3的走时值。S0t3t2 t0、t1经过（3-15）式，“十”字形的四个顶点值都可以求取出来，这样十字形的各个顶点构成一个环，下面要做的是如何把该环进行外推，求取整个模型上的走时。我们把距离前一个环上的极小和极大走时点最近的点定义为该环的相对极小和极大走时点。对于沿z方向边界上的第一个点可由方程（3-11）的非中心有限差分求出，即方程（3-16a）， （3-16a）t3是待求C和D的走时，t0是图3-2b中点O的走时（相对极小值走时），t1和t2、分别为A和B点的走时。同理可求得沿x方向边界上的第一个点的走时为： （3-16b）图3-3 扩展矩阵图 R1、R2、R3和R4分别为源点O到下边、左边、上边和右边的距离。先由源点O进行矩阵扩展直到R1环，如虚线所示；然后基于R1环向左、向上和向右外推直到R2环；再由R2环向上和向右外推至R3环；最后由R3环向上外推至R4环，这样整个模型上的走时就可计算出来。R1R2R3R4O式（3-15）和（3-16）考虑了平面透射波的情况。当计算下一环走时时，在四条边上依次进行（当震源点在地下时，震源下面的点的走时不必要计算）。逐个从相对极小走时点出发，沿边界向前、向后逐点扫描计算走时，直到遇到极大走时点或顶点为止。当所有四个边上的网格点的走时都计算出来后，再计算四个顶点的走时。在扫描过程中，相对极大值被计算了两次，选其中较小的一个作为该点的走时，这等价于几何射线从两边到达同一个点，我们只考虑了初至射线。另外，当（3-15）式右边根号下出现负值时有， （3-17）当（3-16）式方程的右边根号下出现负值时，我们选取： （3-18）当向z方向外推时，当向x方向外推时。式（3-17）和（3-18）考虑了首波出现的情况。若考虑顶点的散射波，则有： （3-19）式（3-19）考虑了散射波传播的情况。依据费马原理，在考虑了平面透射波、首波和散射波情况下的最终初至走时为： （3-20）具体的实现步骤如下：如图3-3所示，从源点到各边的长度分别为R1、R2、R3和R4，以R1R2R31时，其中为波数，上式可简化为（Berkhout, 1983） （3-28）（3-28）式是波场延拓的计算公式，需要的仅是计算偏移成像点的成像值，这时只从外推的记录上取时间的波场值，即为r点的成像振幅值，用公式表示为 （3-29）如果将反射系数引进积分式，（3-29）式可写为 （3-30）R(r)表示反射系数加权后的偏移成像点的波场；A(r, rs)表示从震源到成像点（反射点）的振幅；A(r, rg)表示从成像点到接收点的振幅。绕射走时的射线追踪法或有限差分法见本章的第一节和第二节。2采集和传播效应及其补偿1）传播效应反射波由于在地面和反射层之间的传播而被模糊，偏移就是要去掉传播效应。地震能够搞清的细节是半波长，更小的细节是不能解决的，但这并不意味着小比例尺的介质参数可以忽视。通过一维介质的充分研究已经表明，层间的多次反射可能严重地影响着地震波场的可传播特征，其主要影响是与角度有关的频散现象。但当前的宏观模型对这种影响不予考虑，因此这种影响在偏移中也被忽略了，这会导致成像波形频散和错误的振幅随入射角的变化（AVA）。要进行真振幅偏移，就是要对几何扩散和反射系数随入射角的变化和透射损失等进行补偿。在震源附近，地震子波是脉冲型的。但在远离震源后的传播中子波具有了一定长度的延续时间。研究表明：反射面越深，反射波的优势频率也越低，这种变化主要与波场衰减有关。研究波场衰减为仪器设计和参数设置提供参考、对高保真和高分辨率地震信息的精确接收十分有价值。反射波的波场损失主要包括三方面：发散损失、透过损失和非弹性衰减损失。发散损失是波前扩散的必然结果。无论反射界面或速度梯度存在与否都要发生这种损失，并且它们的存在控制着具体的发散特点。透过损失是只有当部分能量被反射的反射界面存在时才发生的一种损失，界面越多损失越大。因此，视反射界面的分布情况，透过损失可以大于或小于发散损失。而地震波的非弹性吸收机理目前还不是很清楚，较为一致的意见是由于介质的内摩擦和振动做功引起的。这与地震能量转换为热能的吸收过程有密切关系，也和物质的不均匀性引起的散射而使信号失去相干性有关。其中，岩石孔隙和裂隙中的流体对地震波的非弹性衰减的研究于20世纪50年代，由Biot建立的均匀各向异性多孔线性粘弹性介质中波传播的半唯象理论展开。但由于其复杂性，目前对其具体的定量研究还没有开始。总之，波场传播衰减的研究对于指导爆炸和保真接收均有实际意义。2）采集效应对于复杂地区的三维地震数据采集有许多限制，常采用一些不规则观测系统。Kirchhoff偏移在三维叠前深度偏移领域中已经取得了很大的成功，然而对于一些盐丘或者是采集观测系统比较复杂的情况，偏移振幅的精度不够高。这主要是由不规则采集观测系统的采样和波场穿过地下构造复杂地区的扭曲相互作用而产生的。目前，不规则空间采集分布是三维地震勘探中一个常见的问题。Chemingui和Biondi已经说明了不规则采样在偏移成像剖面上会留下印痕，尤其是在3-D多道处理和3-D叠前偏移中更是如此，为此设计了许多使这些假象最小化的方法。为了考虑到不规则空间分布，近来提出了两种不同的方法：一种是局部方法，它是权函数基于等效数据理论（Canning and Gardner, 1998）；另一种是全局方法，介绍的是加权真振幅偏移核的反演（Wu, 1998）。为了计算合适的保幅叠前深度偏移的权函数，Philippe等提出了地面道位置的几何研究。由此把上述两种方法综合成3-D保幅叠前深度偏移方法。在这种方法中，保幅偏移的权函数包含了振幅补偿和采集观测系统补偿两大部分。对于采集观测系统补偿应该考虑道密度和采集（传播）效应。把一特定的密度权函数因子直接包含在偏移核中，这样就考虑了不规则采集的影响，由此应能提高最终的保幅成像结果。注意，其中的加权面是利用对整个观测系统的每一炮检距范围的中心点位置的三角剖分来建立的，且对每一叠前道给出一加权因子。因为处理四维空间（震源和检波点）是比较困难的，因此把工作空间简化为一理想的三维地震观测，其中对每一炮检距范围考虑中点坐标（和）。每一个对目的层有贡献的地震道必须关于地面分布加权，此地面分布可以为每一地震道关于它的相邻道（这样寻找最接近的道）定义一个积分区。在我们的三维工作空间内，在每一炮检距范围，把相邻道定义为相邻中点。3）偏移核内的权因子三维叠前深度偏移是基于渐进近似（射线理论+Born或Kirchhoff近似），格林函数可以在三维不均匀光滑速度场中通过基于动力学射线追踪的波前重建方法求得。另外，即使成像算法设计成保幅的，但由于地震数据的不规则空间采样，振幅假象和相位畸变将会影响成像的质量。因此，不规则采集的三维数据体的偏移成为一个重点研究的课题。下面主要说明一下如何把地震道的不规则采样的影响考虑到叠前深度偏移的核函数内。在一海上三维地震情况下，对于每一深度点x0，三维偏移/反演的最终表达式可概述如下 （3-31）其中，q是一个矢量，它是由与分别来自于炮点和接收点位置(s, r)的射线有关的慢度矢量和ps+pr确定的，A和T是由射线追踪计算的振幅与走时。是地震数据，LSR是地震道。在三角剖分和权函数计算后，能够对每一地震道LSR定义适当的密度，并可以用包含一传播项和一采集项的二个矩阵表示雅可比矩阵偏移核通过矩阵补偿了双程振幅（1/A），并通过矩阵考虑了采集观测系统。在该矩阵中，含有由射线追踪计算的旁轴量和一个加权密度（wmid）。wmid(s, r)来自于三角剖分，它直接考虑了采样的不规则性。因而，在不均匀模型中，基于三维射线追踪的偏移方法具有如下几个优点：计算效率高；能提供有关反射层的定量信息（即反射系数、速度扰动或阻抗反差）和能解释多值波至。另外，三角剖分是纯几何问题，步骤比较简单。以二维为例，（1）输入中点谱（对一给定炮检距范围的中点坐标），输出是一系列三角形（中点为顶点）。因为不想内插地震数据，中点坐标是固定的。每一中点与其它中点通过三角形的边相联接，边的数目取决于该中点与其它中点的相对位置；（2）划分每一三角形面并把部分面分配给最近的中点。如果在保幅叠前深度偏移的核中选用合适的权因子，则可以消除实际采集所产生的印痕。三角剖分和权函数的计算实现了偏移算子的保幅能力。3地震波的散射地震剖面的真振幅分析日益成为地震勘探研究的重点。一个反射层的波阻抗反差和沿该反射层的振幅变化用来定量表征反射层的物理和地质性质。对于勘探问题，关于反射层性质的推论有可能确定油/气/水界面。在地震勘探中广泛接受的观测方式是：一般反射层越深，相对应的反射信号越复杂。一个可能的原因是由于随着深度的增加反射层本身会变得更加复杂；其二，复杂的覆盖层在很大程度上影响地震波从震源到反射层再到检波器的传播。对于第二种情况，即使在反射层本身很简单的情况下，在传播过程中由于散射会造成波的累积复杂性和能量损失。此时，如果采用这样的记录信号做振幅研究，那么将欠估计反射系数，且按照该反射层的内部结构，将对偏移剖面做错误的解释。下面给出了对给定的覆盖层的统计参数，如何定量估计散射损失。根据这些信息，就可以校正由于散射引起的振幅损失，并依次可靠地估计反射系数和AVO趋势。Kirchhoff型真振幅叠前深度偏移成像是生成振幅解释所依据的剖面的最新技术。通常用于Kirchhoff积分中的加权函数仅仅考虑了炮点和检波点之间的几何扩散损失。有些人已尝试考虑吸收和衰减的综合影响，另一些人做了薄层介质影响的研究。更进一步地，人们已经认识到了统计不均匀覆盖层对反射振幅的影响，但是却没有尝试校正相应的振幅损失。下面描述的统计不均匀二维介质对振幅的影响是依据广义的ODoherty Anstey理论，这种方法是把前面所提到的散射损失在Kirchhoff积分中引入附加的权函数。权函数依赖于覆盖层的统计参数、信号的主频和其传播路径的长度。实现过程简单直接，只需少量计算即可。1）随机介质中的散射损失地震波通过岩石传播的基本性质是吸收和衰减。为了描述脉冲传播和定量表征散射衰减，采用了利用Born近似的平均场理论或走时校正的平均场形式的理论方法。平均场理论高估了散射衰减；而走时校正的平均场排除了大波数，因而没有把大规模不均匀性考虑进去。此外，需要正确合理地选取截止波数，这往往仅能通过数值试验来确定。Shapiro and Hubral提出了广义ODoherty Anstey理论，用来在一维随机介质中描述一次波场。基于Rytov and Bourret近似，Muller and Shapiro描述了在二维和三维随机介质中的地震脉冲传播。这一工作用因果性原理进行了延伸，以便于对在随机不均匀弹性固体中传播的瞬态平面波推出Green函数。这一过程可以理解为ODoherty Anstey理论向二维和三维随机介质的一种延伸。对于二维介质中的点源情况，通过类比法，可以得到下列的关于格林函数的近似式 （3-32）其中，为散射衰减系数，（3-33）且表示相位增量 （3-34）在方程（3-32）-（3-34）中，表示波数，为常数背景速度，是传播距离。是扰动谱，它含有介质扰动的二阶统计量，也就是说，在岩石中纵波（或横波）速度的方差和相关长度。表示Heaviside阶跃函数。函数分别表示Fresnel余弦和正弦积分。如果，按照波参数，（3-32）式中的格林函数的有效范围为。其中为波长。注意到：方程（3-32）同时也局限于弱波场扰动范围。传播损失由散射衰减系数方程（3-33）来表征。根据方程（3-33），图3-9描述了品质因数的倒数作为无量刚波数ka的函数用不同的扰动谱作为参数表示的曲线。在图3-9中的品质因数的倒数由进行了归一化处理。 图3-9 2-D随机介质中的ka与1/Q的关系曲线 其中，L/a=50，用不同的扰动谱作为参数，波是从一点震源辐射的2）Kirchhoff型真振幅叠前深度偏移Kirchhoff型叠前深度偏移通常表述为对观测地震波场U的一加权积分，这一积分过程是沿相应的绕射曲线在地面对所有接收点进行的 （3-35）加权函数取决于地下绕射点坐标和接收点坐标，且主要考虑震源和检波点之间的几何扩散损失。为了校正散射效应，必须用上面推出的逆格林函数与观测波场进行褶积。这样，扩展后的偏移公式变为 （3-36）格林函数依赖于传播路径的长度，因而它是震源点、检波点和地下绕射点坐标的函数。这意味着从理论上我们必须计算格林函数，且对每一震源点、检波点和地下绕射点完成褶积。这会导致计算走时的误差增加。下面给出了一些近似，这些近似最终变为实用的公式，且在计算中仅有小的总误差，这样计算总误差的增加可以忽略。第一步近似是在求取格林函数时，不对频率进行积分。而是从观测波场，选择一个主频，忽略相位增量，计算一与时间无关的逆格林函数因子 （3-37）这意味着在Kirchhoff积分中与逆格林函数的褶积变成与这一因子的相乘 （3-38）最花费时间的是散射衰减系数的计算，它直接依赖于传播路径的长度。幸运的是这种依赖性是近似线性的，这样没有必要对每一个地下绕射点计算，而是用预先计算的粗略参考表，通过在偏移过程中的线性插值得到。下面是一个应用散射理论对模型偏移的例子。图3-10是所采用的模型。图3-11的上半部分是对均匀模型的正演模拟炮记录，基于均匀宏观模型，应用散射理论进行克希霍夫叠前偏移的结果；图3-11的下半部分是对统计不均匀模型的正演模拟炮记录，基于均匀背景速度宏观模型，应用散射理论进行克希霍夫叠前偏移的结果。图3-12是图3-11中各偏移剖面沿反射层的振幅。其中黑线是均匀介质偏移结果的反射振幅，红虚线是对不均匀介质未做校正的反射振幅，红实线是对不均匀介质校正后的反射振幅。 图3-10 统计不均匀模型 红框为偏移区域 图3-11 Kirchhoff叠前偏移剖面 上图：炮记录是基于均匀介质，偏移是利用均匀宏观模型下图：炮记录是基于不均匀介质，偏移是利用均匀背景速度宏观模型图3-12 图3-11中的偏移结果沿反射层的振幅变化曲线黑线为均匀模型情况，点红线为校正后的不均匀模型，实红线为未经校正的不均匀模型4真振幅权函数真振幅偏移公式与叠前深度偏移公式相似，但真振幅核项与积分加权项不同。反射点M的成像积分式为 （3-39）其中，M为反射点坐标；为地面点坐标；是从炮点到反射点M的走时；是从反射点M到地面点的走时。真振幅核为 （3-40）式中，是垂直方向的单位矢量；s(SM)和s(MG)为SM段和MG段的慢度；r为反射点M的水平面上的径向坐标。在单一界面情况下 （3-41）式中和分别为入射角和反射角；为速度。表示为Beylkin行列式（Beylkin,1985; Bleistein,1987） （3-42）Beylkin行列式可以通过动力学射线追踪计算出来，真振幅权函数是基于Beylkin行列式的反演积分。从上述可以看出，真振幅偏移主要是补偿地震波在传播过程中的各种损失，使偏移后的振幅能得到保持。尽管理论上的推导是可行的，但离实际计算还存在一段距离。目前，这段距离在逐步缩小。Schleicher等（1993）给出了三维真振幅有限偏移孔径偏移，它是基于绕射叠加原理提出的，所有的绕射叠加都是基于对地下每一绕射点沿惠更斯面完成一加权求和，求和可以用下面的数学积分来表示 （3-43）其中，和分别表示从到M和从M到的旅行时。权函数是采用稳相法对叠加积分的渐进估算得到的，并且在计算时采用两次动力学射线追踪，一次从震源点到成像点，另一次从成像点到接收点。这是因为复权函数的模是借助于射线变换子矩阵来给定，而复权函数的相位因子是由分支射线的焦散数目来确定。权函数的形式如下 （3-44）其中，下标s, g和p分别表示震源点、接收点和成像点，下标sp和gp分别表示从震源点rs和接收点rg到成像点P的传播，是一对震源-接收组合参数化的向量，和是描述震源-接收组合的22矩阵，（或）是一22常矩阵，（）是起始锐角（或出射角），它是地面法线与连接（或）和P点的射线之间的夹角，是震源点的纵波速度。Schleicher方法与其它方法的不同之处在于使用了显式零级射线理论（例如，描述了一次反射波的旅行时和振幅值），因此叠加运算和其性质的推出可以通过几何理论来理解。这一方面克服了Born近似存在的问题；另一方面可以直接在弹性介质中进行保幅偏移。由此可以得到一新的在绕射叠加积分中涉及到的权函数的表达式。并且考虑了沿射线路径存在的焦散问题，因此使存在波场焦散情况下的震源脉冲恢复的问题得以解决。由于三维Kirchhoff型真振幅叠前深度偏移通常采用动力学射线追踪计算真振幅权函数，所以计算量特别大。孙建国等人（1996）说明了利用运动学射线追踪可以很简单地计算射线变换子矩阵和分支射线的焦散数目，并给出了一新的真振幅共炮集权函数的计算公式。该公式如下 （3-45）依据以上理论，对一盐丘模型分别进行了常规Kirchhoff型叠前深度偏移和Kirchhoff型真振幅叠前深度偏移。其中图3-13a 是没有做真振幅处理的Kirchhoff型前深度偏移剖面，图3-13b 是做真振幅处理后的Kirchhoff型叠前深度偏移剖面。从图上我们可以看出，应用真振幅处理后的偏移剖面的分辨率有明显提高，盐丘内部的构造特征更清晰。图3-13 (a) 常规Kirchhoff型叠前深度偏移剖面 (b) Kirchhoff型真振幅叠前深度偏移剖面5薄互层效应偏移成像详细表征了地震波长一半时的结果，很明显，薄层问题不能很好地得到解决。在偏移中忽略了与薄层有关的依赖于角度的反射和相位扭曲，这将导致频散的成像结果和错误的AVA效应。波在层状介质中的传播特性与波长有关，当波长足够长，或地震波的频率足够低时整个介质可以与一个均匀的各向同性的介质等效，当波长与层的厚度可比时，会有散射波出现，这种波长的波显示出频散效应。在地震勘探频带范围，实际地层厚度都较小，即d，薄层反射系数随阻抗差的增加而增加，随层厚度的增大而线性增加，随频率增高而线性增加，并且薄层的反射响应是高频滤波器，透射系数在调谐频率处具有极小值，且与该频率处反射系数的极大值对应。因此在调谐频率以下，薄层的透射响应是低通的。一系列薄层，不管是否是周期性沉积，都存在一个相应的低通透射效应。上述说明薄层的反射响应是高通的，而透射响应是低通的。因此要想做到真振幅偏移必须考虑薄互层的效应。C.P.A.Wapenaar等人（1996）提出了考虑薄层的波动方程法真振幅偏移，它是通过考虑薄层效应的改进型匹配滤波来消除AVA畸变和相位畸变，从而实现没有频散图像的真振幅偏移。由积分法与波动方程法之间的关系，可以在Kirchhoff型真振幅偏移中，通过考虑薄层效应，从而消除波的频散和振幅相位畸变，达到真振幅偏移的目的。6成像孔径的优化成像孔径在应用Kirchhoff型偏移方法中起着关键的作用，成像孔径选取的大小和形状会明显影响计算效率、成像的信噪比和成像的保幅性。传统的全孔径Kirchhoff型偏移是把一时间样点模糊到沿一准椭圆的成像空间上。当偏移一单道信息时，成像孔径是整个成像空间，因此成本高、噪音干扰强。稳相法是处理成像孔径问题的另一种方法，它需要计算振幅加权函数并与叠加道相乘。它的成像孔径是一在地面道域中的概念，且地下成像点处的真振幅通常是研究的目的。它也是一种成像驱动的方法，首先完成一全孔径Kirchhoff偏移，并把所得到的成像结果用来求取加权计算所必需的稳相射线；然后用局部化到稳相道附近的成像孔径进行另一Kirchhoff偏移。波径偏移法也是处理成像孔径问题的一种方法，它也利用了稳相法的特点，是稳相法的一种延伸，波径法与传统的成像驱动方法的不同之处就在于波径法偏移依赖于从接收点到成像空间发射的镜像（反射）射线。很显然这是一种数据驱动的方法，且成像孔径是一在成像空间域的概念。在波径偏移法中，需要对一特定反射波至通过几道数据计算一入射角，然后利用该入射角计算镜像射线，最后沿该射线确定镜像反射点。与成像驱动的方法相比，波径偏移的计算效率更高，因为它只需一次偏移，并且成像孔径由离散Fresnel带组成。但是它也存在缺陷，对强反射很敏感但对弱反射的偏移不足。为了充分利用波径偏移的优点和进一步提高构造细节的成像效果，人们提出了一优化的Kirchhoff型偏移（OKM）方法。它是一种将常规Kirchhoff偏移和波径偏移相结合的联合偏移方法。利用波径偏移来计算一单道的镜像反射点和确定密集充填镜像成像点的成像孔径；然后用Kirchhoff偏移把该道偏移到成像孔径上。与成像孔径由离散Fresnel带组成的波径偏移相比，该联合方法的成像孔径既包含Fresnel带又包含该带周围的小空间，这些小空间大部分是与弱反射层有关的Fresnel带。这样，通过OKM法对强、弱反射层都能很好的成像。相比较而言，这种方法的计算成本要高于波径偏移法，但低于传统的Kirchhoff偏移法。1）稳相解在地下成像点处的叠前Kirchhoff偏移成像公式可以用下式表示 （3-46）其中，是震源子波的振幅谱；和分别代表震源和检波器的位置；是几何扩散因子；是相位函数，其中是观测旅行时，和分别是地震能量从到和从到的计算旅行时；积分范围为一检波孔径。应用稳相原理，且在积分范围内假定一个一阶稳相点，这样公式（3-46）可以近似表示为 （3-47）其中，是的稳相点，。对公式（3-47）的数值实现，需要求取稳相点。可以在成像域求取，但事先必须做全偏移处理。公式（3-47）的求取方法有底-顶射线追踪和顶-底射线追踪。2）一优化成像孔径下面给出了一利用一优化成像孔径的Kirchhoff偏移方法。其算法如下，首先把自动增益控制应用到各个共炮点道集上，然后用下述方法确定的一优化孔径逐道偏移炮道集。优化偏移孔径的实现步骤如下：（1）用大于一预定门槛值的振幅拾取反射波至。（2）利用在几个地震道上的局部倾斜叠加对拾取的波至计算入射角。（3）从检波点向地下介质发射射线。（4）沿该射线路径寻找镜像反射点。（5）确定含有该镜像反射点的一成像孔径。利用抛物回归分析确定该孔径的中心位置，在该位置镜像反射点分布最密，孔径的宽度由剩余标准偏差控制。（6）把该道偏移到优化成像孔径上。3）数值实例依据以上优化成像孔径的Kirchhoff偏移方法对Marmousi模型作了试算，而且与常规的Kirchhoff偏移和波径法Kirchhoff偏移作了对比。图3-14是以上三种方法对单道进行偏移的结果，目的是比较它们的成像效果。从图3-14看到：常规Kirchhoff法有较强的偏移噪声；波径法只突出了强反射；而优化孔径法对强、弱反射都有较好的成像效果。图3-15是应用以上三种方法对Marmousi模型偏移的结果。从图中可以看出，优化孔径的Kirchhoff偏移无论在剖面的分辨率还是在同相轴的连续性等方面都有明显的改进，尤其是对弱反射同相轴的成像效果有较大改善。 (a) 常规Kirchhoff偏移 (b) 波径法Kirchhoff偏移 (c) 优化孔径的Kirchhoff偏移图3-14 对Marmousi模型单道的成像结果(a) 常规Kirchhoff偏移 (b) 波径法Kirchhoff偏移 (c) 优化孔径的Kirchhoff偏移图3-15 对Marmousi模型的成像结果7各向异性真振幅偏移上世纪九十年代以来，对各向异性的偏移问题给予了极大的关注。随着油气勘探对地震成像精度要求的提高，对偏移速度精度的要求也越来越高。速度的各向异性在偏移成像中对反射层位置的确定起着重要的作用。研究表明（Larner, 1993）：如果偏移中忽略横向各向异性问题，那么在该介质中的偏移误差对陡倾角反射来说是很大的。用各向同性的偏移算法处理横向各向同性（TI）介质的数据会产生平反射面的错误的位置，甚至漏掉陡构造。目前已经开始研究各向异性介质的偏移算法。真振幅叠前深度偏移的目的是在偏移处理中保持地震记录数据的振幅和相位信息。所得到的成像结果包含了有关在每一成像点处的反射系数随偏移距变化的信息（AVO），从而可以间接推出有关地下介质的弹性参数和岩性反差的信息。基于真振幅偏移技术对弹性各向异性介质的2-D合成数据进行了试算，结果表明除了能实现反射层归位外，还能提供依赖于入射角的反射系数以及AVO信息。因为旅行时和振幅在Kirchhoff型真振幅偏移中起核心作用，因此它们的计算必须是精确的。可以利用射线追踪通过三角剖分的各向异性模型计算初至的旅行时和振幅波场。图3-16a为各向异性推覆体模型，其中浅颜色区表示各向同性介质，深颜色区表示各向异性介质。各向同性区的P波和S波的速度分别为2.74km/s和1.38km/s；各向异性（横向各向同性）区的快P波和慢P波的速度分别为3.365km/s和2.925km/s，vs/vp的比值为0.5214。对应的弱各向异性参数为，和。图3-16b为各向异性推覆体模型的炮记录，图3-16c为该模型的初至旅行时图；图3-16d为该模型的初至振幅谱；图3-16e为该模型的叠前真振幅偏移成像结果。图3-16a 各向异性推覆体模型 图3-16b 由图(a)模型生成的炮记录图3-16c 图(a)模型的初至旅行时图 图3-16d 图(a)模型的初至振幅谱上面所使用的各向异性