求解弹性波动方程的频率域近似解析离散化波场模拟方法

摘要：

为提高频率域弹性波动方程数值求解的计算效率,本文引入近似解析离散化(NAD)方法将其进行数值离散并得到大型线性代数方程组.在详细分析了相应系数矩阵的稀疏分块结构与数学性质之后,本文提出采用不精确旋转分块三角预处理子加速Krylov子空间迭代方法来快速求解该线性方程组,并利用数值试验证实这种方法在弹性波场模拟方面的数值效率.通过与另外两种经典数值方法(常规有限差分方法和交错网格有限差分方法)对多种介质模型进行波场模拟、数值频散分析以及与解析解的波形对比,NAD方法显示了其在压制数值频散和提高计算效率方面的优势以及对复杂介质模型弹性波场数值模拟的有效性.

关键词：频率域弹性波动方程;近似解析离散化;预处理迭代方法;波场模拟;频散分析

所属期刊栏目：64

分类号：P315(大地（岩石界）物理学（固体地球物理学）)

资助基金：中国地震局地壳应力研究所中央级公益性科研院所基本科研业务专项资助项目;国家自然科学青年基金;北京市教委科技一般项目

在线出版日期：2021-08-23（万方平台首次上网日期，不代表论文的发表时间）

页数：共20页

页码：2838-2857