Modeling 3-D Elastic Wave Propagation in TI Media Using Discontinuous Galerkin Method on Tetrahedral Meshes

Abstract

Transversely isotropic (TI) medium is a widely studied anisotropic solid medium in seismology. The numerical simulation of seismic wave propagation in TI media is an effective tool to analyze the mechanism of seismic waves in complex anisotropic media. In this study, we introduce a double-weighted Runge-Kutta discontinuous Galerkin (RKDG) method for numerically solving wave propagation problems in 3D TI media with surface topography. This method incorporates the discontinuous Galerkin spatial discretization with an explicit double-weighted two-step iteration time discretization. The local Lax-Friedrichs flux is used as the numerical flux in the formulations. This method can solve the first-order velocity-stress seismic wave equations including TI media as well as more general anisotropic media. Due to the large scale of 3D problems, the parallel technology is adopted. Regions with irregular boundaries are discretized into unstructured tetrahedral meshes. Numerical experiments for various anisotropic media including the vertical and tilted TI media are presented. The results demonstrate the effectiveness of the double-weighted RKDG method in wavefield simulations in 3D complicated anisotropic media.