Analysis of elastic wave simulation accuracy with discontinuous Galerkin finite element method based on triangular meshes

被引:0
|
作者
Han D. [1 ,2 ]
Liu W. [1 ,2 ]
Si W. [1 ,2 ]
机构
[1] SINOPEC Geophysical Research Institute, Nanjing
[2] SINOPEC Key Laboratory of Geophysics, Nanjing
关键词
Discontinuous Galerkin finite-element method; Numerical dispersion; Numerical dissipation; Runge-Kutta time scheme; Stability analysis; Triangular mesh;
D O I
10.13810/j.cnki.issn.1000-7210.2021.04.009
中图分类号
学科分类号
摘要
Accuracy analysis is the foundation for numerical simulation of seismic waves. With regard to the numerical stability, dispersion, and dissipation of the discontinuous Galerkin finite element method (DGFEM) based on triangular meshes, a triangular periodic mesh model is constructed, which can be used to study the effects of different triangular elements on simulation accuracy. The theoretical and numerical results show that the stability condition of the DGFEM based on the Runge-Kutta time scheme is related to the shape of triangle elements. The maximum time step for stable mode-ling has a linear relationship with the radius of the inscribed circle of the element, and the equilateral triangle element has the least rigorous stability condition. Meanwhile, the wave field from DGFEM simulation based on the local Lax-Friedrichs flux shows weak dispersion but strong dissipation, and both dispersion and dissipation present directivity in the periodic mesh. In addition, the logarithm of the modeling error has a linear relationship with that of the mesh size. The numerical experiments compare the influence of different mesh shapes on the wave field and verify the theoretical directional difference. The results of this paper can provide a theoretical basis for the triangular mesh division, parameter setting, and selection of numerical flow in DGFEM. © 2021, Editorial Department OIL GEOPHYSICAL PROSPECTING. All right reserved.
引用
收藏
页码:758 / 770
页数:12
相关论文
共 46 条
  • [1] PEI Zhenglin, WANG Shangxu, DI Bangrang, Et al., Study on acquisition methods without geophone array by numeric modeling of 3D seismic wavefield, Oil Geophysical Prospecting, 40, 2, pp. 138-142, (2005)
  • [2] CUI Hongliang, WANG Ruizhen, CHEN Jingguo, Et al., Application of numerical simulation in vibroseis seismic acquisition, Oil Geophysical Prospecting, 52, 2, pp. 209-213, (2017)
  • [3] HE Baoqing, XIE Xiaobi, A numerical investigating on wide-line seismic data acquisition and processing [J], Oil Geophysical Prospecting, 54, 3, pp. 500-511, (2019)
  • [4] MA Guoqi, CAO Danping, YIN Jiaojian, Et al., Numerical simulation of detecting seismic signals in DAS wells, Oil Geophysical Prospecting, 55, 2, pp. 311-320, (2020)
  • [5] WANG Wenfeng, YUE Dali, ZHAO Jiyong, Et al., Research on stratigraphic structure based on seismic forward modeling: A case study of the third member of the Yanchang Formation in Heshui area, Ordos Basin, Oil Geophysical Prospecting, 55, 2, pp. 411-418, (2020)
  • [6] Virieux J, Calandra H, Plessix R E., A review of the spectral, pseudo-spectral, finite-difference and finite-element modelling techniques for geophysical imaging[J], Geophysical Prospecting, 59, 5, pp. 794-813, (2011)
  • [7] Kelly K R, Ward R W, Treitel S, Et al., Synthetic seismograms:A finite-difference approach, Geophysics, 41, 1, pp. 2-27, (1976)
  • [8] Dablain M A., The application of high-order differencing to the scalar wave equation, Geophysics, 51, 1, pp. 54-66, (1986)
  • [9] DONG Liangguo, MA Zaitian, CAO Jingzhong, Et al., A staggered-grid high-order difference method of one-order elastic wave equation, Chinese Journal of Geophysics, 43, 3, pp. 411-419, (2000)
  • [10] Yang D., A nearly analytic discrete method for acoustic and elastic wave equations in anisotropic media[J], Bulletin of the Seismological Society of America, 93, 2, pp. 882-890, (2003)