地震波走时有限差分法计算.pdfVIP

第 3s卷 第 l期 地 球 物 理 学 报 vo1．3S，No．1 1992年 1月 ACT^ GEOPHYSICA SINICA Jan．．1992 地震波走时的有限差分法计算 朱金明 王丽燕 (石油地球物理莓I探局，河北嚣州 072750) 摘 要 类似于 Claerbout思想，将波动方程在射线理论申前对偶——镜像方程写为上、下行波方 程的形式．对变换得到的守恒型偏徽分方程 ，用 E—O格式，可计算 出任意速度构造的二维两 格上各点的地震波走时．文中证 明了守恒型偏微分方程的通量为一凸函数，采用单调的守恒 型差分格式，精度符合要求．计算速度比目前所用的各种计算走时的方法都快得多，算法也便 于向量化． 关键词 地震波走时，有限差分，镜像方程． 一 、 引 言 地震渡走时的计算在地球物理学中起着极重要的作用．例如，确定震中位置时，较古 老的方法是，便所有地震台观测到的走时的剩余值最小作为准则；Kirchhoff偏移需要计 算 Green函数，而 Green函数描述了观测点与地下介质 内部点之间的走时．走时可按多 种方法计算 ．对均匀介质 ，走时很易计算；对水平层状介质，波的传播方向将偏离直线而 成 曲线；对纵、横 向变化的介质，走时计算一般通过射线追踪实现．射线追踪法基于以下 原理： 高频近似条件下地震波传播轨迹遵循射线方程 ．但射线追踪法有如下问胚： 首 先 ，速度变化剧烈时，两点间可能出现多条射线路径 ，很易丢失最有意义的那条射线路 径——最小走时路径；其次，对多源点、多接收点作射线追踪 ，需要作大量存储和大量内 插 ，算法常难以实现；第三 ，介质性质剧变处常出现射线相交、射线难以穿越浅层等现象， 内插非常困难． Vidalc提出了一种在二维规则网格上直接计算走时的方法 ．基本原理是：给定某 一 网格点(如源点)的初始走时值，用平面波或球面波前近似外推计算相邻各点的走时． 它避免了走时插值问题 ，但仍有缺点，第一，方法只能从所谓的 “相对最小走时点 开始外 推 ，否则算法将不稳定；第二 ，必须知道每一网格点处的波前 曲率 ，如，在震源附近 ，常常得 用比较笨拙的球面波前近似作外推．这不仅便其计算成本提高 ，也使算法不易向量化，不 易推广． 本文的方法同 Vidale方法一样，可在发散区域、浅层区中计算出地震波走时，且 比 Vidale方法计算速度更快 ，便于向量化． 但方法只限于研究单一方 向传播的波，且计 算 出的走时只取单一值 ，这可能对某些应用有限制 ，但对我们所惑兴趣的几个问题 ，如 本文 1990年 5月 28日收到，1991年 2月 20日收到任改稿 ． 1期 朱金明等 ：地 震波走时的有限差分法计算 87 Kirchho~f型计算 [31、层析反演、提高波动方程有限差分正演计算速度 、和 VSP层速度迭 代反演 ，本方法是非常适宜的． 二、 原 理 为描述方便起见，下按二维进行描述 ，若耍推广到三维在原理上是直接的．方法研究 的是，介质某处设一震源 ，求介质内各网格点的最小射线走时． 1．镜像方程 二维几何射线或二维波动的波前是 由射线的镜像方程 (鲁 +(告 — 控制的 ，其中 、=为坐标轴 ， 为慢度 ，l为时间场函数．镜像方程的解可 以按波前和射 线来解释 ． 注意到镜像方程与声波波动方程 丝 Oxz+。等0=2一、’警0 、(2) 形式上的类同，事实上，(1)式是(2)式的特征线方程 (詈)+(詈)一 (