有限频层析成像方法.docVIP

有限频体波层析成像简介 徐小兵 赵亮 传统射线走时层析成像基于高频射线理论。在地球内部的速度异常变化不大，且异常体的尺寸要远大于地震波的特征波长的情况下，可以根据费马原理确定地震波沿走时最短的路径，此时可以认为到时差主要取决于射线路径上的速度结构。然而，实际的地震波具有有限频率，当地震波经过尺寸较小或相近（相对地震波波长）的速度异常体时会发生绕射现象。随着传播距离的增加，地震波波前面由于非均匀介质导致的超前或滞后特征能逐渐愈合甚至消失（Nolet and Dahlen, 2000；Hung et al., 2001）。地震波的这种有限频率特性表明地震到时差并不主要取决于射线路径上的速度结构。针对上述问题，Dahlen和Hung 等人推导了地震到时在有限频地震波传播波场下的灵敏场函数（或根据其几何形状被简称为“香蕉—甜面圈”理论） （Dahlen et al.,2000；Hung et al., 2000），进而发展了有限频层析成像方法。相对传统射线走时层析成像，该方法由于更真实地反映地震波的传播特征，因此能够更好地重建地球深部结构速度图像。 一、有限频层析成像方法原理 地震波的有限频率特性决定了不同频率的地震波的走时对射线周围三维速度结构具有不同的敏感程度。Dahlen 等（Dahlen et al., 2000）结合体波传播理论和线性Born单次散射近似来模拟通过地球内部任意一点的散射波到达台站的波形扰动，利用不同频率的体波与未经散射的直达体波之间的相互干涉结果，得出有限频率走时变化与中心几何射线附近第一菲涅尔带内的三维速度扰动之间的关系为： ……(1) 式中，是经过波形互相关处理得到的散射波与直达波的走时差；为空间点x处的体波速度扰动量；是该点三维灵敏度核，代表了某个散射点x处的速度异常对走时变化的贡献。从上式可以看出，敏感场K不仅是空间位置x的函数，同时还是频率w的函数，即不同频率的地震波对不同空间范围的速度结构进行了采样。 在射线理论下地震波走时表示为慢度沿射线路径的线积分，而在有限频率灵敏核理论中地震波走时可以更精确地表示为灵敏度核函数在射线周围一定空间内的体积分，这是有限频率理论与射线理论最本质的差别。因此，确定在一定速度模型下的有限频率灵敏核函数是进行有限频层析成像的关键所在。 二、有限频率地震体波灵敏度核函数及其性质 获得有限频率地震体波灵敏度核函数是开展有限频层析成像的基础。 图1. 有限频体波走时在球对称介质下的几何路径示意图（Hung et al 2004） ……（2） 公式（2）表明有限频敏感度核函数是对时间窗内所有单点散射波进行积分求和的形式。由于对所有满足条件的单点散射波进行射线追踪是一个非常耗时的过程，因此在实际成像过程很难实现。针对此情况，Dahlen等（2000）在计算敏感场函数时引入旁轴近似，即认为走时残差只对第一菲涅尔带内的速度扰动比较敏感，从而来有效的减少射线追踪的运算量。采用的旁轴近似的公式为： ..........(3) 图2. 公式2得出的直达P波灵敏度核函数（Hung et al 2000） 图3. 根据旁轴近似（公式3）给出的直达P波灵敏度核函数（Hung et al 2000） 从图2和图3的对比中我们可以看出通过旁轴近似由公式3得出的灵敏场可以很好地对由公式2得出的理论敏感场进行近似。 图5. 灵敏度核函数与频率的关系（Hung et al., 2004） 由高、中、低三个不同频带内的灵敏度核函数（图5）的形态可以看出，地震波的频段越高其灵敏度核越窄且强度越大，最接近于传统的射线；而频带范围越低则灵敏度核越宽，这时的灵敏度核强度也会降低。理论分析表明，对宽频带波形分频带处理，不同频带内的分析结果可以分别约束中心射线路径周围不同范围处的速度结构。 三、有限频理论与射线理论的比较 图6. 冰岛下方射线理论与有限频方法的灵敏场函数分布对比。（a）射线方法的灵敏场；（b）有限频方法的灵敏场。图件来自Hung et al., 2004。 图6展示了射线方法和有限频率方法的灵敏场分布：相对射线方法，有限频率方法对空间的采样更加密集、平滑，且空间范围更大，因此能够充分利用宽频带资料中的信息，减少数据不均匀覆盖的影响，进而能更加充分对速度异常体进行重建。如同很多对比实验所示，有限频