A Legendre spectral element method with optimal mass matrix for seismic wave modeling

Liu ShaoLin,Yang DingHui,Meng XueLi,Wang WenShuai,Xu XiWei,Li XiaoFan
DOI: https://doi.org/10.6038/cjg2022Q0145
2022-01-01
Abstract:Numerical solving seismic wave equation of motion plays an important role in investigating wave phenomenon and imaging subsurface structure. The accuracy of a numerical method for forward modeling can eventually determinates the effectiveness of numerical result and the image quality. At present, the spectral element method (SEM) is successfully applied to seismic wave propagation in models with different scales. However, the accuracy of SEM is not fully explored because the integrations involved in the mass matrix of SEM are not accurately approximated by Gauss-Lobatto-Legendre (GLL) quadrature rule. The numerical errors originating from numerical quadrature rule reduce the accuracy of SEM. In this article, we propose a strategy to increase the accuracy of SEM by constructing an optimal mass matrix. The optimal mass matrix is developed in two steps. First, we design a least-square problem in terms of the GLL quadrature weights and the integrations on the diagonal of mass matrix. Second, we adopt the conjugate gradient method to solve the least-square problem and obtain the optimal quadrature weights. The optimal quadrature weights can reduce the errors associated with the mass matrix and eventually increase the accuracy of SEM. The validity of the SEM with optimized mass matrix is demonstrated by numerical dispersion analysis and numerical examples.
What problem does this paper attempt to address?