复杂介质中地震波模拟的无单元法-中国科技论文在线.PDF

中国科技论文在线 石j由地球如理参探 第41卷 第1期 · 综合研究 · 复杂介质中地震波模拟的无单元法 贾晓峰 胡天跃 王润秋 (北京大学地球 与空间科学学 院) (石 油大学 (北京 )资 源与信息学院) 摘 要 贾晓峰，胡天跃 ，王润秋．复杂介质 中地震波模拟的无单元法．石油地球 物理勘探 ，2006。41(1)：43～48 复杂介质 中地 震波的传播 多是通过求取单程或双程波动方 程的数值解进行模拟 的。波动方程的数值 解法 多种 多样 ，本文将工程和岩石力学领域 中广泛应用的一种无 网格算法 ‘ 无单元法应用于地震模拟 。由于无单 元法的节点完全独立存在．因此在计算成本上与有限元相 比优势明显；再采用滑动最小二乘 的拟合方法使得无 单元法具有精度高、j虫立变量解高攻连续等优点。本文在介绍无单元法原理的基础上 ，通过薄膜震动算例表 明 该方法与有限差分法相 比有精度优势：然后将无单元法用于叠前地震模拟。并讨论 了几种吸收边界方法在地震 模拟中与无单元法的结合。模型计茸 的结果表 明。辅以有效的吸收边界，无单元法 能够获得波场完整、走 时精确 的理论地 震图，其精度和稳定性 皆令 ^＼满意 。 关键 词 地震模拟 波动方程 无网格算法 无单 元法 吸 收边界 精度 量，因而计算成本有所降低 。在有限元法中，形函 1 引言 数 的构造是个难题 。一方 面 ，形函数 的插值阶次不 能太高 ，否则很难满足其必备 的正交完备性 ；另一 地震波的传播过程理论上可 以用波动方 程描 方面 ，形函数阶次太低会损害精度 。而无单元法对 述 ，因此人们一直致力于探讨各种各样 的波动方程 形函数的要求则大大降低，仅需满足一个最小二乘 数值解法 以获得满意的正演剖面 。近 2O年来 ，有限 的精度即可 。 差分法和有 限元法是这一研究领域 的热点 。这两种 (2)有 限元法多采用拉格 朗 日插值 ，解是二次不 方法 虽经过 多年 的研 究和发展 已 日臻成熟 ，但依然 连续的；但无单元法则放弃插值的思想，而采用最小 在如有限差分 的稳定性和精度及有 限元的成本等许 二乘拟合的方法，因此所得解是高次连续的。 多方面存在一定 问题 。 由于在解 的连续性上 的独特优势，无单元法 已 从工程力学 中发展起来 的无单元法 (Element— 经被成功地应用于如破裂模拟等力学研究领域[1]， FreeMethod，简称 EFM)一。。可望在地震勘探 中获得 但在地震勘探方面的应用却依然是空 白。本文在验 有效的应用 。无单元法是近年来大量涌现 的无网格 证无单元法对波动 问题的求解效果之后，结合几种 方法 中的一种 。它是针对有 限元法所面临的预 流行的吸收边界方法进一步探讨该算法在地震波高 处理复杂 、独立变量一次导数不连续等 问题应运而 精度模拟方面的应用意义 。 生 的。无单 元 法 采 用 滑 动 最小 二 乘 法 (Moving LeastSquares，简称 MIS)来拟合 波场 函数 ，它 虽 2 基本原理 然继承 了有 限元法 的一些特点．但更具备某些独特 的优势 ：