徐世刚^1,2,,
刘洋^1,2,,
1. 中国石油大学(北京) 油气资源与探测国家重点实验室, 北京 102249
2. 中国石油大学(北京) CNPC物探重点实验室, 北京 102249

基金项目: 国家自然科学基金项目（41474110）和国家重大科技专项课题（2017ZX05018-005）联合资助


详细信息

作者简介: 徐世刚, 男, 1991年生, 中国石油大学(北京)博士研究生, 主要从事各向异性地震波传播与成像等方面的研究工作.E-mail:xushigang2310@163.com

通讯作者: 刘洋, 男, 主要从事地震波正演、成像、反演及多分量地震等方面的研究工作.E-mail:wliuyang@vip.sina.com

中图分类号: P631

收稿日期:2017-07-18
修回日期:2017-11-29
上线日期:2018-07-05

3D acoustic and elastic VTI modeling with optimal finite-difference schemes and hybrid absorbing boundary conditions

XU ShiGang^1,2,,
LIU Yang^1,2,,
1. State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China
2. CNPC Key Laboratory of Geophysical Prospecting, China University of Petroleum, Beijing 102249, China


More Information

Corresponding author: LIU Yang,E-mail:wliuyang@vip.sina.com

MSC: P631

--> Received Date: 18 July 2017
Revised Date: 29 November 2017
Available Online: 05 July 2018


摘要
摘要:传统有限差分系数是通过泰勒级数展开求取的，这样导致所计算的频散曲线在大波数区域会产生较强的数值误差.针对二阶空间偏导数的显式有限差分离散，本文发展了一种新的优化差分系数方法：首先将泰勒级数展开与多点采样方法结合应用于空间频散关系，基于最大范数建立直观有效的优化目标函数，采用Remez算法求解该目标函数，从而获得最优化差分系数.利用优化有限差分方法求解三维垂直对称轴横向各向同性（VTI）介质中的声波和弹性波方程.另外，本文将二维混合吸收边界条件推广到三维VTI介质中，用于吸收人工截断边界反射；基于各向异性特征，合理调整了边界区域的速度值来提高吸收效果.考虑到三维情况下计算效率的问题，本文波场外推过程中采用图形处理器（GPU）取代传统的中央处理器（CPU）.数值精度分析表明，相比较于传统的泰勒级数展开方法，优化有限差分方法在大波数区域对频散误差的压制效果更明显.在三维均匀和修改的Hess VTI模型中的数值模拟实验证明了本文方法具有更高的精度与效率，混合吸收边界条件在三维VTI介质中具有良好的边界吸收效果.
关键词: 数值模拟/
有限差分/
优化方法/
各向异性波动方程/
混合吸收边界条件/
图形处理器

Abstract:The coefficients of the conventional finite-difference (FD) stencil are commonly determined by the Taylor-series expansion method (TEM), which can produce major numerical dispersion errors at large wavenumbers. To solve this problem, we develop a novel optimal explicit FD method (FDM) for the second-order spatial derivative. We first combine the TEM and multiple sampling to the spatial dispersion relation, and then construct an intuitive and effective objective function using the maximum-norm method (MNM) and solve it by the Remez iteration algorithm (RIA). The maximum-norm (MN)-based FDM is adopted to solve the 3D acoustic and elastic vertical transversely isotropic (VTI) wave equations. Moreover, we extend the 2D hybrid absorbing boundary conditions (HABCs) to 3D acoustic and elastic VTI media, and improve the absorbing effectiveness through reasonably adjusting wavefield velocities in the absorption areas based on anisotropy properties. To improve computational efficiency for a 3D case, we utilize graphic processing unit (GPU) instead of traditional central processing unit (CPU) architecture. Numerical accuracy analyses indicate that the optimal scheme has better tolerance to the numerical dispersion at large wavenumbers than the conventional TEM. Numerical modeling experiments for 3D homogeneous and modified Hess VTI models demonstrate that the proposed schemes have greater accuracy and efficiency than the conventional schemes and achieve significant absorbing effects.
Key words:Numerical modeling/
Finite-difference/
Optimized scheme/
Anisotropic wave equation/
Hybrid absorbing boundary condition/
Graphic processing unit



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_14602