适于复杂地质模型的三维射线跟踪方法.pdfVIP

第 卷 第 期 地 球 物 理 学 报 , 年 月 』 , ’ 适于复杂地质模型的三维射线跟踪方法 杨长春 冷传波 李幼铭 国 地 理 所 , 北京 叨 中 科学院 球物 研究 摘 要 复杂地质体模型 的三 维射线跟踪方法在数值计算时随波 阵面 自动生 成射线跟 踪所 需 的 计算单元 , 无 需对复杂地 质体模型 预先进 行模型离散化 基于 计算 中 自动生 成 的计算单 元 , 进行射线跟踪计算时就无需求解射线路径 与计算单元数值边 界 的交 点, 因而 具有甚 高 的计算效率 本文方法不仅可作为地震 资料叠前成像 的算法基础 , 而 且还 为复杂介质 中波 场高频渐近解 的计算提供一种最佳的计算途径 关键词 射线跟踪 , 地震 资料叠前成像 , 计算单元 自动生成 , 复杂地质模型 引 言 近 , 对于 日益 杂的地质勘查 目 , 杂地质体构造成像研究 已成为热 门课 年来 复 标 复 题 这就要求人们从地震叠前资料 出发 , 利用层析成像方法求取反 映地 下介质速度分布 的速度模型 , 继而 利用基于波动方程 的叠前深度偏移方法来实现复杂构造 的精确成像 除层析成像和偏移成像 的 自身算法外, 射线跟踪 的数值计算精度和效率则直接影 响着介 质成像质量和实用化 为此 , 对地震 波场射线解 的快速有效算法开展 了多方面 的研 究, 并获得令人瞩 目的进展 在射线跟踪方法 中, 一类属于全局算法 , 这类算法在射线跟踪计算过程 中, 同时考 虑空 间所有离散点上 的走时和 射线路径 的计算 , 如程 函方程 的差 分解法 一 , 波前法 ’ 和最短路径法 一 等 另一类属 于局 部算法 , 它在计算过程 中, 只考虑两 点之 间的射线 路径和走时计算 , 即每条射线 的计算独立 于其它射线 , 如程 函方程 的积分解法 等 现有 的这些射线跟踪方法 都需要 在射线跟 踪计算之 前把地质模型 离散化成 网格单 元 , 以实现地质模型 的数值化描述 于是射线跟踪算法则转 向侧重于处理 网格单元 内和 穿过 网格边界时的射线路径和走 时 全局 法大多 用立方 网格 , 地质 型 据 计算 算 使 体 使 模 数 离散化 对于三 维射线跟踪 以及复杂地质模型而 言, 复杂地质体 的速度分布形态需要用 大量 的正方体单元才能精确地 予 以描述 于是 , 射线跟踪计算过程就必须涉及大量 的网格 体单元和 网格体边界 由于全局算法 同时考虑所有空 间离散点上 的走时和射线路径 的计 算, 从而可 部分地缓解 了因使用立方体单元造成 的计算效率过低 问题 , 但其