时空域波动方程混合网格有限差分数值模拟

• 作者: 刘威 杨哲 王艳香 王述江 胡自多 韩令贺
• 来源期刊: 岩性油气藏
• 学科: 地球科学
• 关键词: 混合网格 有限差分 数值模拟 时空域
• 年，卷（期）: 2018,（2）
• 所属期刊栏目: 技术方法
• 页码范围: 93-109
• 页数: 17页
• 分类号: P631.4
• 字数: 13065字
• 语种: 中文
• DOI: 10.12108/yxyqc.20180211

摘要

传统2M阶有限差分格式(Tradtional 2 Mth-order Finite-Difference Schemes,T2 M-FD)和时空域2 M阶有限差分格式(Time-Space Domain 2 Mth-order Finite-difference,TS2 M-FD)均是目前应用较普遍且具代表性的高精度有限差分方法.T2 M-FD仅基于空间域频散关系求解差分系数,模拟精度较低.TS2 M-FD基于时空域频散关系和平面波理论求解差分系数,模拟精度较高.T2 M-FD和TS2 M-FD的差分格式相同,都是只利用常规直角坐标系中坐标轴上的网格点差分近似波动方程中的Laplace算子,而没能充分利用旋转直角坐标系中距离差分中心点更近的网格点来进一步提高模拟精度.本次研究提出利用常规直角坐标系和旋转直角坐标系中的网格点一起差分近似波动方程中的Laplace算子,并将Laplace算子表示为常规直角坐标系中M个Laplace算子和旋转直角坐标系中N个Laplace算子的加权平均,构建出一种新的混合2 M+N型有限差分格式(M2 M+N-FD).推导出M2 M+N-FD基于时空域频散关系和平面波理论的差分系数计算方法,进行频散及稳定性分析.频散分析表明:与T2 M-FD和TS2 M-FD相比,M2 M+N-FD能更有效地压制数值频散,模拟精度更高.稳定性分析表明:M2 M+N-FD和TS2 M-FD的稳定性基本相当,比T2 M-FD的稳定性强.最后,利用M2 M+N-FD进行均匀介质和层状介质模型的数值模拟试验,并将其推广应用于Marmousi模型的逆时偏移,高精度的数值模拟结果和偏移成像质量证明了M2 M+N-FD的优越性和普遍适用性.

作者信息

序号	姓名	单位	发文数	被引次数	H指数	G指数
1	胡自多		20	127	7.0	11.0
2	王述江		8	109	5.0	8.0
3	刘威		36	49	3.0	6.0
4	杨哲		28	255	7.0	15.0
5	王艳香		12	45	4.0	6.0
6	韩令贺		7	31	3.0	5.0