刘威1,2,
雍学善1,
王小卫1,
韩令贺1,2,
田彦灿1,2
1. 中国石油勘探开发研究院西北分院, 兰州 730020
2. 中国石油天然气集团有限公司油藏描述重点实验室, 兰州 730020
基金项目: 国家重点研发项目"面向E级计算的能源勘探高性能应用软件系统与示范"(2017YFB0202905)和国家科技重大专项"下古生界-前寒武系地球物理勘探关键技术研究"(2016ZX05004-003)联合资助
详细信息
作者简介: 胡自多, 男, 1971年生, 博士, 高级工程师, 主要从事波动方程数值模拟与成像方法研究.E-mail: huzd@petrochina.com.cn
中图分类号: P631 收稿日期:2020-08-04
修回日期:2021-03-18
上线日期:2021-08-10
Mixed-grid finite-difference method for numerical simulation of 3D wave equation in the time-space domain
HU ZiDuo1,2,,LIU Wei1,2,
YONG XueShan1,
WANG XiaoWei1,
HAN LingHe1,2,
TIAN YanCan1,2
1. Research Institute of Petroleum Exploration and Development-Northwest, PetroChina, Lanzhou 730020, China
2. Key Laboratory of Petroleum Resources of CNPC, Lanzhou 730020, China
MSC: P631
--> Received Date: 04 August 2020
Revised Date: 18 March 2021
Available Online: 10 August 2021
摘要
摘要:常规高阶和时空域高阶有限差分方法广泛应用于三维标量波动方程的数值模拟，这两种差分方法仅利用笛卡尔坐标系中的坐标轴网格点构建三维Laplace差分算子，相应的差分离散波动方程本质上仅具有2阶差分精度，模拟精度低.本文将三维笛卡尔坐标系中非坐标轴网格点分为两类：坐标平面内的非坐标轴网格点和坐标平面外的非坐标轴网格点，系统推导出了两类非坐标轴网格点构建三维Laplace差分算子的方法，进而提出了一种利用坐标轴网格点和非坐标轴网格点共同构建三维Laplace差分算子的混合网格有限差分方法，并利用时空域频散关系和泰勒展开建立差分系数方程，推导出了差分系数的通解.相比常规高阶和时空域高阶差分格式的2阶差分精度，时空域混合网格差分离散波动方程理论上能够达到任意偶数阶差分精度，模拟精度显著提高，同时稳定性更强.频散分析表明：相比常规高阶和时空域高阶差分格式，在计算效率基本相同时，时空域混合网格差分格式能更有效地减小数值频散，减弱数值各向异性，模拟精度更高；在模拟精度基本相当时，混合网格差分格式能采用更大的时间采样间隔，计算效率更高.数值模拟实例进一步验证了混合网格差分格式在提高模拟精度和计算效率方面的先进性，也验证了其普遍适用性.
关键词: 差分格式/
混合网格/
差分系数算法/
数值频散/
三维Laplace差分算子
Abstract:Conventional high-order and time-space domain high-order Finite-Difference (FD) methods have been widely used for 3D wave equation numerical simulation. These two FD schemes adopt the same FD stencil which only uses the axial grid nodes to construct the 3D Laplace difference operator, and the corresponding discrete difference wave equation can only reach 2nd-order accuracy. In this paper, we classify the off-axial grid nodes in the 3D Cartesian coordinate system into two categories: off-axial grid nodes in the coordinate plane and off-axial grid nodes outside the coordinate plane, and systematically derive the methods to construct the 3D Laplace difference operator with these two categories of off-axial grid nodes. Then a new kind of 3D mixed-grid FD scheme is proposed, which adopts a novel FD stencil combining the axial nodes and off-axial nodes to construct the Laplace difference operator and the FD coefficients are calculated based on the time-space domain dispersion relationship and Taylor expansion. The discrete difference wave equation derived from this new mixed-grid FD scheme can reach 4th-order, 6th-order and arbitrary even-order difference accuracy. So it can significantly improve the modeling accuracy comparing to conventional high-order and time-space domain high-order FD schemes, and also has better stability. Dispersion analysis shows comparing to conventional high-order and time-space domain high-order FD schemes, the mixed-grid FD scheme can more effectively suppress the numerical dispersion and weaken the numerical anisotropy to obtain higher modeling accuracy with almost the same computational efficiency, and it can obtain higher computational efficiency by adopting larger time interval with almost the same modeling accuracy. Numerical modeling experiments further verify the superiority of the mixed-grid FD scheme in improving the modeling accuracy and computational efficiency, and also demonstrate its universal applicability.
Key words:Finite-difference scheme/
Mixed-grid/
Difference coefficients algorithm/
Numerical dispersion/
3D Laplace difference operator
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=bf65655f-9b88-4c25-8eb8-05ee45ecfb13
