时间-空间高阶精度矩形交错网格隐式有限差分声波正演模拟 |
| |
引用本文: | 王静,刘洋,周泓宇.时间-空间高阶精度矩形交错网格隐式有限差分声波正演模拟[J].地球物理学报,2023(1):368-382. |
| |
作者姓名: | 王静 刘洋 周泓宇 |
| |
作者单位: | 1. 中国石油大学(北京)油气资源与探测国家重点实验室;2. 中国石油大学(北京)CNPC物探重点实验室;3. 中国石油大学(北京)克拉玛依校区 |
| |
基金项目: | 国家自然科学基金项目(42274147,41874144); |
| |
摘 要: | 有限差分方法因其操作简单、计算消耗低而成为地震勘探领域中最为常用的数值模拟方法之一,然而用离散的显式差分算子数值逼近地震波动方程中的连续导数容易导致数值频散,并且基于正方形网格离散形式的有限差分方法对不同地质模型的适应性较低.针对一阶变密度声波方程的数值模拟,本文发展了一种适用于矩形网格离散形式的时间高阶空间隐式有限差分格式,可以有效压制时间和空间频散,同时灵活的网格剖分增强了其应用的广泛性.基于本文矩形交错网格时间高阶空间隐式有限差分格式的时空域频散关系和变量替换的思想,首先采用泰勒级数展开方法求解不同方向的非轴上时间差分系数及轴上空间差分系数,使本文差分格式可以获得任意偶数阶时间和空间精度.为了进一步提高本文差分格式在更大波数区域的空间模拟精度,我们采用线性优化方法来求取新的轴上空间差分系数用于一阶变密度声波方程的波场迭代求解中.频散、稳定性分析及数值模拟算例表明:相比于传统十字形空间域隐式有限差分格式,本文矩形交错网格时间高阶空间隐式有限差分格式在精度、稳定性和效率方面均具有优势.
|
关 键 词: | 矩形交错网格 隐式有限差分 线性优化算法 时间、空间高阶精度 |
|