第3章Kirchhff积分法叠前深度偏移.docVIP

第三章 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）