地震波场数值模拟方法.docVIP

第 42卷第 2期 石 　 油 　 物 　 探 V ol. 42,N o. 2 文章编号 :100021441(2003 0220143206 地震波场数值模拟方法 张永刚 (中国石油化工股份有限公司科技发展部 , 北京 100029 摘要 :简要总结了地震波场数值模拟的各种方法的基本原理及其主要特点 , 对最近在该领域出现的一些方法和 研究结果做了简要的阐述 , 并对比了各种方法的优缺点 。 在此基础上提出了运用波动方程数值模拟作为基础 , 结合射线方法辅助识别波场类型 , 用于分析异常波的产生机理和出现特点的基本思想 , 这对复杂条件下的地震 勘探具有指导和借鉴意义 。 关键词 :地震波场 ; 数值模拟 ; 射线追踪 ; 有限元 ; 伪谱法 ; 正演模拟 中图分类号 :P63114+1　　　　　　 文献标识码 :A On numerical simulations of seismic w avefield Zhang Y onggang (Department of Science and T echnology Development , , , Abstract :This paper reviews the principles of of seismic wavefield , and com 2 pares the merits and defects the and results are briefly discussed. The author pro 2 poses to study of abnormal waves based on wave equation numerical simulation supplemented by K ey w ords :seismic wavefield ;numerical simulation ;ray tracing ;finite element ;pseudo 2spectrum ;forward m odeling 　　 地震波场数值模拟是研究复杂地区地震资料 采集 、 处理和解释的有效辅助手段 , 地震波场数值 模拟的主要方法包括 2大类 , 即波动方程法和几何 射线法 。 波动方程数值模拟方法实质上是求解地 震波波动方程 , 因此模拟的地震波场包含了地震波 传播的所有信息 , 但其计算速度相对于几何射线法 要慢 。 几何射线法也就是射线追踪法 , 属于几何地 震学方法 , 由于它将地震波波动理论简化为射线理 论 , 主要考虑的是地震波传播的运动学特征 , 缺少 地震波的动力学信息 , 因此该方法计算速度快 。 因 为波动方程模拟包含了丰富的波动信息 , 为研究地 震波的传播机理和复杂地层的解释提供了更多的 佐证 , 所以波动方程数值模拟方法一直在地震模拟 中占有重要地位 。 1　 几何射线法 几何射线法又称为射线追踪法 , 其主要理论基 础是 , 在高频近似条件下 , 地震波的主能量沿射线 轨迹传播 。 基于这种认识 , 运用惠更斯原理和费马 原理来重建射线路径 , 并利用程函方程来计算射线 的旅行时 。 在旅行时计算中可以应用有限差分等 方法 , 以期获得快速的解 。 随着地震勘探技术的发 展 , 新的射线追踪技术也不断涌现 , 以满足大的数 据处理 (如三维数据 和较高精度要求下对复杂地 质体研究的需要 。这些技术主要研究焦点是如何 精确地划分地质体 , 如何实现旅行时场的快速准确 的计算以及对已有方法的改良 。 射线法的主要优点是概念明确 , 显示直观 , 运 算方便 , 适应性强 ; 其缺陷是应用有一定限制条件 , 计算结果在一定程度上是近似的 , 对于复杂构造进 行两点三维射线追踪往往比较麻烦 。为了计算波 沿射线的旅行时和波沿射线的振幅变化 , 首先都必 须知道波的传播路径 。 所谓射线追踪 , 狭义来说指 的就是根据地震波的传播规律确定地震波在实际 地层中传播的射线路径 。 在地震学中 , 有 2类地震射线追踪问题 :一类 是 1点射线追踪 , 即已知射线初始点位置和初始出 射方向求地震波的传播路径问题 ; 另一类是 2点射 线追踪 , 即已知射线初始点和另一观察点 (接收点 的位置 , 不知射线初始出射方向 , 求 2点之间的射 线路径问题 。用射线理论制作 VSP 模型 , 不论是 收稿日期 作者简介 :张永刚 (1956— , 男 , 高级工程师 , 现从事科技发展规划 和管理工作。 零偏移距或非零偏移距 , 因为震源和接收点不是同 一点 (除零偏移距井口记录道外 , 遇到的都是 2点 射线追踪问题 。 1 有关初值问题的试射法 (Shooting method 。 这种方法根据由震源出发的一束射线