地球物理资料数字处理第二十六讲.docVIP

克希霍夫积分法波动方程偏移 （Kirchoff Integral Wave Equation Migration） 克希霍夫积分偏移建立在波动方程积分解的基础上,这种偏移方法一般不受反射界面倾角的限制,但是偏移归位后的剖面上噪声干扰背景比较强. 本节要点： 克希霍夫积分公式； 克希霍夫积分法偏移； 频率—波数域波动方程偏移； 一 地震偏移的克希霍夫积分公式（Kirchoff Integralformula） 如果介质是均匀的,各向同性和完全弹性的,纵波波动方程为: (7-3-1) 如图7-9所示,S为其内部无源的任一闭曲面, 是内的任意观测点, 是 面上的任意点，已知闭合曲面的波动函数求曲面内部空间任意一点上的波场值的克希霍夫积分为： （7-3-2） 其中为闭合曲面的外法线，r是从点到Q点的距离： 在地震勘探中，我们的工作平面是地表面，此时闭合曲面可理想化为两部分：,表示地表平面，是无穷远处的边界面。如图7-10所表示，当R（把无穷远处的边界面看成是一个半径趋于无穷大的半球面）时，就代表下半空间区域。因此可以认为面上的扰动对点无贡献，因此（7-3-2）式的闭曲面积分仅与地面观测平面有关，这时有： 图7-9 内部无源的任一闭曲面 图7-10 闭合曲面 （7-3-3） 在地平面上，因为外法线方向与坐标方向相反，因而 所以，（7-3-3）式又可以写成： （7-3-4） 上式中含有项，但在偏移问题中，我们无法提供此值，所以为了用克希霍夫积分公式做偏移，必须从上式中消除此项。 为此，取关于地表面的镜象点，点位于下半空间之外，点的波场值为零，在闭曲面的克希霍夫积分公式为： （7-3-5） 其中， 在平面上，有， 因此，（7-3-5）变为： 将与上式与（7-3-4）相加得： （7-3-6） 根据爆炸反射界面成像原理，用代替（7-3-6）式中的 ，并且地表面上的波场值是由界面点发出的，因此还应该用负速度取代中的速度，得到： （7-3-7） 上式的意义就是利用地表时刻的波场值求取点的时刻的延拓波场值，这正体现了延拓与正演的过程相反，是地震波的反传播过程。 根据爆炸反射界面成像原理，应用=0的时刻的延拓波场来构组偏移输出，因此得到克希霍夫积分法偏移的基本公式： （7-3-8） 在地表面上，，因此， 代入 （7-3-8） 得到 ： （7-3-9） 把深度换成时间则： 代入（7-3-9）中 得到 ： = = （7-3-10） 实际上，不可能无限大，只能是某一工区面积 ，或是，工区某一部分面积，假设为，则（7-3-10）就变成： = （7-3-11） 上式就是最终的克希霍夫积分法偏移公式，其中可以在地表水平叠加时间剖面上对各时间道关于求导得到 二 克希霍夫积分法偏移（Kirchoff Integral Wave Equation Migration） 由于地下介质的速度变化情况不同，我们分三种情况来讨论克希霍夫积分法偏移的实现》 常速介质克希霍夫积分法偏移 这里所说的常速介质有两种含义：一是地下介质本身确实是常速的，二是偏移时仅利用速度谱上的读取的速度，因为这时实际上已经把用时间表示的深度以上的地层用一常速介质所代替，这种情况应用比较广泛。 另外，在实际处理中，范围，一般并不是三维面积观测的整个地表范围，而是，取一个滑动的有限范围（二维时为一段，三维时为一块），称为偏移孔径，如取道，此时，（7-3-11）式的离散形式为： （7-3-13） 令 = （7-3-13） 则（7-3-12）可以写成： + （7-3-14） 其中， （7-3-15） 常速介质克希霍夫积分法偏移的特点是地下任意一点的偏移都是利用了地表水平叠加时间剖面做边界值，采用的是非递推的延拓形式。 层状介质克希霍夫积分法偏移 在层状介质中，