Modified symplectic scheme with finite element method for seismic wavefield modeling

SU Bo,LI HuaiLiang,LIU ShaoLin,YANG DingHui

DOI: https://doi.org/10.6038/cjg2019M0538

2019-01-01

Abstract:Finite element method (FEM) with triangular mesh has the properties of flexibility and adaptability in complex medium. However, the traditional finite element method is inefficient in seismic wavefield modeling because this method commonly requires a large amount of computation resources. Here we propose a so-called improved kernel matrix storage (IKMS) strategy to improve the computation efficiency of FEM. The improved FEM does not need to assemble global stiffness matrix and the memory requirement is many times less than that of the conventional FEM. In terms of temporal discretization, the elastic wave equation after FEM discretization is first transformed into a Hamilton system. Then, an additional spatial discretization term is added to the second-order explicit Runge-Kutta-Nystrom (RKN) scheme. Based on Taylor series expansion, we obtain a fourth-order symplectic scheme with all positive symplectic coefficients. Theoretical analysis shows that the temporal-spatial modified numerical scheme is superior over the traditional symplectic FEM in suppressing numerical dispersion and increasing numerical stability. Furthermore, numerical results verify that this method requires less computer resource and achieves higher numerical accuracy.