二维复杂介质中地震波走时和射线的计算方法及其应用 - 中国地震信息网.pdf

维普资讯 第29卷 第3期 地 震 研 究 VOl_29．No．3 2006年 7月 JOURNAL OF SEISMoLoGICAL RESEARCH Ju1．，2006 二维复杂介质 中地震波走时和射线 的计算方法及其应用’ 潘纪顺 张先康 赵成斌 张景超 (1．中国地震局地球物理研究所，北京 100029；2．中国地震局地球物理勘探中心，郑州450002； 3．河南省地质调查院，郑州450007) 摘要：将 Podvin和Lecomte的精确局部格点走时计算方法和Schneider等人的动态规划方法结合起来，可得到一 种快速、有效的有限差分波前计算方法。使用该方法对各种类型地震波的走时和射线的计算进行了讨论 ，并给 出了这种有限差分走时计算在叠前深度偏移中的应用实例。 关键词：有限差分；走时；波前 ；射线 ；偏移；复杂结构 中图分类号：P315．3 l 文献标识码：A 文章编号：1000—0666一(2006)03—0245—06 对比模型中算法失败。Qin等人 (1992)提出扩展 0 引言 波前法，采用一个扩张网络结构，对 Vidale的方 法进行了改进 ，解决了初至波走时的计算问题，但 目前，地震波走时和射线的计算方法主要有两 该方法同样只考虑前向传播问题，不能处理强速度 大类 ：射线追踪法和有限差分法。射线追踪法的基 对比介质。Podvin和Lecomte(1991)提出了一种 本思想是无限高频近似的地震能量遵循一个射线追 方法，基于平面波假设，在计算 中融合了费马原 踪方程表述的轨迹 ，它又分为试射法和弯曲法。通 理，在每个格点考虑前 向和反 向两个方向的透射、 常射线追踪存在三类问题 ：首先，对于强变化速度 衍射和折射，能够处理强速度对 比介质。但是，如 场，连接两点间的射线有许多条，容易错过最小走 果在近源点没有额外的格点时，平面波的假定是不 时；其次，如果射线路径过多、点过多，则计算量 准确的。Schneider等人 (1992)提出一种动态规 很大；第三，即使在光滑介质中，也会存在阴影 划方法，使得计算的效率大大提高。Cao和 优n— 带。因此，射线追踪法对结构简单、数据量较小的 halgh (1994)提出一个角节点算法，改善了计算 地震走时计算是 比较有效的。但随着地震观测技术 的精度。Zhao (1996)利用弯曲波前假设对 Pod— 的快速发展、探测系统的不断完备、地震数据 的大 vin和Lecomte的平面波波前假设进行了改进，使 量积累，对复杂介质结构的探测需求 日益加大，而 得计算的准确度得到了提高。笔者将 P0dviIl和Le— 射线方法在精细复杂结构中的应用又往往失败 。 comte的局部格点走时计算方法和 Schneider等人 Vidale(1988)提出了一种有限差分方法来求 的动态规划方法结合起来，得到一种快速、有效的 解程函方程 ，其基本思想是 以波前替代射线在介质 有限差分波前计算方法，并使用这种方法对各种类 中传播，用中心扩张法模拟地震波场 的向外扩张， 型的地震波走时和射线的计算进行了讨论，给出了 然而该方法只能计算前 向传播 的波。1990年， 这种有限差分走时计算在叠前深度偏移 中的应用 Vidale(1990)将其提出的算法推广到三维状态， 实例 。 首波得到正确处理，阴影带得