地震偏移成像基本原理课件.ppt

④ 根据成像原理，对所有地下点（z0）取t=0时的波场值，即可实现三维偏移成像。此时，成像值为 （1.2.33） Kirchhoff积分法的偏移效果见图1-23。 三．有限差分法波动方程偏移 下面讨论使用有限差分法对水平叠加地震剖面的偏移问题。为了把上行波方程表示为空间-时间域的表达式，需要把上行波方程表示为某种近似式。然后在空间-时间域研究其差分方程及求解问题。最后讨论一些计算方法和效果。 1．上行波的空间-时间域方程 为了适应介质速度的空间变化，我们要在空间-时间域中进行偏移成像或地震图的模拟工作。首先就要把上行波方程表示在空间-时间域中，这需要用到某种根式展开。 1）二项式展开 下面我们将用到 这样的二项式展开，在这里我们介绍几种展开式。 （1）Taylor展开 这是一个众所周知的显式展开式，它一般表达为： （1.2.34） 展开条件 。 如果把这种展开式用于微分算子，在不进行辅助处理时将找不到稳定的有限差分方程来解相应的微分方程。因此我们在使用二级近似以上的展开式时不能用这种展开式。 （2）连续分式展开，或称为Pade展开 这个展开式表示为 如下形式（ ）： （1.2.35） 这是一种隐式展开式。其各级展开式如下。 一级展开式： （1.2.36a） 二级展开式： （1.2.36b） 三级展开式： （1.2.36c） 高级展开式可依此类推。 （3）迭代展开 这种隐式展开法，是把前一级的展开结果代入下一级的展开式中。设 则逐次迭代展开式可表示为（ ）： （1.2.37） 用这种展开方法得到的各级展开式如下。 一级展开式： （1.2.38a） 二级展开式： （1.2.38b） 三级展开式： （1.2.38c） 高级近似式可依此类推出来。 从（1.2.36）和（1.2.38）公式组可以看出，后两种展开是等价的。 2）上行波的空间-时间域方程 在第一节已经求出了频率-波数域的上行波方程(1.1.12)式： 用迭代展开法展开上行波方程： （1.2.39） 式中 由（1.2.39）式求出各级近似式如下。 一级近似式： （1.2.40a） 二级近似式： （1.2.40b） 三级近似式： （1.2.40c） 高级近似式可依次类推。现在，我们把（1.2.40a）式转换到空间-时间域，求出一级近似方程。 上行波方程（1.2.40a）式可改写为： （1.2.41） 对（1.2.41）式进行傅里叶反变换： （1.2.42） 根据傅里叶变换的微分性质： （1.2.43a） （1.2.43b） （1.2.43c） 把(1.2.43a)、(1.2.43b)和(1.2.43c)式代入(1.2.42)式，得到： （1.2.44） 上式就是空间-时间域的一级近似的上行波方程，常常被称为 方程。 同理可求出空间-时间域的二级及二级以上近似的上行波方程。经推导，空间-时间域的二级近似的上行波方程为： （1.2.45） 上式常常被称为 波动方程。 3）浮动坐标系中的单程波方程 上行波方程在一定的浮动坐标系中可以简化。 我们对二维波动方程： （1.2.46） 进行如下的坐标转换： 坐标变换前后波场本身是不变的，因此存在： （1.2.47） （1.2.48） 从（1.2.47）和（1.2.48）导出下列导数等式： （1.2.49a） （1.2.49b） （1.2.49c） 将（1.2.49）各式代入（1.2.46）式中得到新坐标系中的波动方程为： （1.2.50） 上式变换到频率-波数域为： （1.2.51） （1.2.51）式可改写为： （1.2.52） 从上式得到下列关系式： 它表示坐标变换前后的算子关系。因此，上行波方程可表示为： （1.2.53） （1.2.53）式的右端项可展开为各级近似式，便得到上行波各级近似方程。 一级近似式 （1.2.54） 二级近似式 （1.2.55） 三级近似式 （1.2.56） 高级近似式可依次类推。由以上各式用前述方法可求出空间-时间域的各级近似方程。下面给出一级和二级近似方程。 一级近似方程 用推导（1.2.44）式那样的方法可以求出一级近似方程为： （1.2.57） 与（1