第4章-波动方程法叠前深度偏移3.doc

§4.6 广义屏法波动方程叠前深度偏移 本节基于波的散射理论，首先从波动方程Green函数解出发，借助Born近似等一系列数学手段，导出广义屏（GS）偏移算子。然后在小角度近似条件下，得到相屏（PS）偏移算子（分步Fourier算子）。为了克服广义屏算子在计算中遇到奇点，特采用Taylor展开方法得到扩展的局部Born近似的广义屏算子。接着我们对相屏算子和扩展的局部Born近似的广义屏（ELBF）算子在数值上做了脉冲测试，并用于凹陷模型的叠后深度偏移处理和复杂模型叠前深度偏移处理。 一．概述 三维叠前深度偏移广泛地应用于地下复杂地质体成像。对当前工作量较大的三维成像问题一般采用基于射线追踪或有限差分走时计算的Kirchhoff积分方法，但该方法在处理复杂构造（如盐丘或陡倾地层）时总存在困难（Hu McMechan, 1986）。虽然基于有限差分方法的全方程逆时偏移对存在剧烈横向速度变化的非均匀介质也具有非常好的成像精度，但是，把该方法应用于三维叠前深度偏移需要花费较多的机时和占用巨大的计算机内存容量，这些要求在目前是不太现实的。于是从20世纪80年代后期直至目前，地球物理学家们一直在寻求一种既方便快捷又准确可靠的叠前深度偏移方法。 要研究波动方程偏移问题，搞清楚波的传播问题是非常必要和有帮助的。在80年代中后期，Knepp与de Hoop等人为了研究波在随机介质中的传播问题，在地震波场的叠加原理基础上，基于速度场分解、波的散射理论以及波动方程的Green函数解法来求解Helmholtz方程，提出了一种在频率-波数域和频率-空间域（双域）交替进行的波动方程解法。这可以说是屏方法的雏形。20世纪90年代初，Wu R.S.与de Hoop等在屏的雏形理论基础上，通过一系列近似处理手段，发展成了较实用的广义屏算法。这类算法既可用于研究波（声波或弹性波）的传播问题，又可用于地震波场成像。该类方法认为速度场可分解为层内常速背景和层内变速扰动。对背景场相当于解常速的声波方程，可通过相移实现；对变速扰动，可认为这种非均匀性相当于散射源（二次源），入射波场作用于这些散射源上，由此产生散射波场。一般假设波场延拓层厚度较小，在薄板近似条件下，对散射场计算式的积分核采用不同的近似方法，如Born近似、De Wolf近似或Rytov近似等等，可以得到不同的散射波场延拓式。随后在此基础上，Wu R.S.和Huang L.J.等人把该类方法发展成了广义屏方法、相屏方法以及局部Born近似的屏方法和局部Rytov近似的屏方法。接着他们针对上述方法的一些不足之处做了进一步的扩展和改进，提出了一系列新的屏算子，如扩展的局部Born近似和Rytov近似的屏方法、拟线性局部Born近似的屏方法，等等。复杂模型及实际地震资料试验证明：这些算子能较准确地描述地震波在随机介质中的传播特征（入射、反射、散射和透射等等），并在波动方程偏移成像应用中具有较好的效果。近年来，西方许多石油公司对该类算法产生了浓厚的兴趣，并逐步投向生产。 本节基于这类屏方法的思路，首先简要介绍各种广义屏算子的基本原理。为了验证相屏算子和扩展的局部Born近似的广义屏算子，我们进行了常速和变速介质中的脉冲响应试算，然后把它们应用于深度偏移成像中，进行了凹陷模型爆炸反射记录的叠后深度偏移试验和对Marmousi模型炮记录的叠前深度偏移试验。结果表明：相屏算子在速度场中等程度横向变化的情况下，能较准确地描述地震波场的传播过程，对较复杂的地质构造具有较好的成像效果，且该算子具有稳定、高效的优点。扩展的局部Born近似的广义屏算子是条件稳定的。它一般要求深度延拓步长要小，且步长与频率呈反比变化，所以最好采用变步长延拓方式。最后是结论与讨论部分。 二．广义屏偏移方法的基本原理 声波方程的Helmholtz形式为： （4-135） 式中，为频率-空间域波场，为水平坐标分量。且，为圆频率，为介质速度场。假定为背景速度，它在层内是一常数。从而（4-135）式可写成： （4-136） 若定义， （4-137） 则（4-136）式可表示为： （4-138） 式中，，等号右边为速度扰动引起的散射源。按照地震波场的叠加原理，如果把总的波场认为是入射波场和散射波场的和，则有： （4-139） 对于背景介质中的波场满足方程：