各向异性介质中波动方程有限元法模拟及稳定性.pdfVIP

第 卷 第 期 地 球 物 理 学 报 , 年 月 , 各 向异性介质 中波动方程有 限元法 模拟及其稳定性 周 辉 徐世 浙 刘 斌 何樵登 , 青 岛海洋大学地质地球物 理研 究所 , 青 岛 长春地 质学 院地球物理 系 长春 摘 要 利用六 面 和 三 函 , 导 出在 意弹性 向异性介质 中三 维 三 量 动 体单元 线性插值 数 推 任 各 分 波 、 用 , 空 方程 所满足 的有 限元方 程 给 出 时域有 限差 分求解 该有 限元方程 时 时 间 间采样 率和 弹性参数应满足 的稳定性条件 提 出避 免每一 时 间步都求解 一大 型线性代数 方程组 的方法 在开始 波场 , 阵 的逆 的 因子 式 , 以后 时 间 只 需 刚 矩 的 的因子 求 前 先求质量矩 形 每 步 作 度 阵 逆 形 式 与 向量 的乘 积 运算 , 可 大 大提 高有 限元 素法 在微机上 进行 波 动方程 数值 模拟 的运 算速 度 两 , 明 方法 的可 , 与 法 了 速度 的 比 , 限 结合 个数值模拟算例 说 该 行性 并 伪谱 作 计算 较 有 元 法 的计算速 度高于伪谱法 的计算 速度 词 向异 , 动 , 有 限 法 , 定 , 关键 各 性介质 波 方程 元 稳 性条件 数值模拟 引 言 波 动方 程 的数值模拟 在 地 震勘探 中起 着 较 大 的作 用