基于波前构建的25维射线走时与振幅计算.pdfVIP

基于波前构建的2．5维射线走时与振幅计算 孙建国何洋 (吉林大学地球探测科学与技术学院地球物理系，长春，130026) 2．5维波场是3维波源在2维介质中产生的场。常规上，在2．5维方法中只考虑在包含场源的垂 直平面内的波场现象。如果采用几何理论模拟波场的传播过程，2．5维模拟方法中的射线路径和传播 时间计算与2维情况下的相应计算完全相同。只不过在计算振幅时要考虑场的3维效应。因此，在缺 少面积性(3维)观测资料的条件下。基于2．5维理论的算法和技术可以考虑场源的3维特点和3维 几何扩散对振幅的影响。 为了实现2．5维快速模拟，我们采用了在上个世纪80年代到90年代之间出现的波前构建法。在 原理上，波前构建法基于传统的射线追踪方程组，但在参数选取上以传播时间代替射线的弧长。因此， 波前构建法的直接输出量是波前面的空间坐标和波前顶上各点处的射线方向。 由于受其高频近似特性的限制，射线理论要求速度在一个波长的范围内变化不大。为了满足这个 要求，需要对速度模型进行光滑处理。在公开的文献中提出的光滑方法主要有两种： (1)改变输入的 结点值。这种方法类似于低通滤波； (2)在输入的结点间采用插值曲线来减小速度梯度，同时保持结 点处的值不变。在本文中，我们采用第一种方法，使用三点平均算子来对速度模型进行光滑处理。 对速度模型进行光滑处理的第二个目的是为了避免在速度或其导数的不连续面上应用反射与折射 定律，加快计算速度。我们的实践表明，适当的光滑可以加快射线追踪的计算速度，但会使计算结果 中含有一定的误差。在极端的光滑处理下，模型的地质结构特征将会发生质的变化。由此而产生的后 果是计算结果将变得与预期结果完全不同。 利用波前构建法进行波场数值模拟得另一个问题是射线插值。由于几何扩散效应的影响，在经过 一定的外推距离之后，射线问的距离将超出傍轴射线理论的有效范围，射线的覆盖率也将会变得越来 越小。因此，要采用适当得方法提高射线的覆盖率。在本文中，我们分别采用了两种射线插值方法， 其所给出的计算精度基本相当。 为了实现2．5雒数值模拟，我们从下列3维动力学射线追踪系统出发， 3 鲁=v2p,(渊驯2 ：： 拿=_V-1罢‘(㈦，2，3) d苫 舐j 这里，x，是位置坐标分量，V是速度矢量，P，是慢摩矢量分量，f是时间。在2．5维条件下，Y方向 的速度不产生变化。因此，加／砂=O，咖J，／dr=0。这说明，Py和dy／dr=v2p，，均为常数。 这样一来，闯题就由3维转化成了2维。相应的射线追踪方程简化为： lIV |l～ V一 p Il |l 出一如出一打 yp 监如盟如 一 V— 加一缸加一七 以(f+Ar)·为了检验上述计算方案的精度，我们与均匀介质中的分析解进行了对比，其误差小于 5％。 · 58·