Frequency domain spectral element method for modelling poroelastic waves in 3-D anisotropic, heterogeneous and attenuative porous media. Issue 2 (15th July 2021)

• Record Type: Journal Article
• Authors: Zhan, Weichen; Zhuang, Mingwei; Liu, Qi Qiang; Shi, Linlin; Sun, Yuefeng; Liu, Qing Huo
• Is Part Of: Geophysical journal international. Volume 227:Issue 2(2021)
• Page Start: 1339
• Page End: 1353
• Publication Date: 2021-07-15
• Subjects: Numerical modelling -- Computational seismology -- Seismic anisotropy -- Seismic attenuation -- Theoretical seismology -- Wave propagation; Geophysics -- Periodicals; 550
• DOI: 10.1093/gji/ggab269
• Languages: English
• ISSNs: 0956-540X

Abstract

SUMMARY: Simulating poroelastic waves in large-scale 3-D problems having porous media coupled with elastic solids and fluids is computationally challenging for traditional methods. It is well established that the spectral element method (SEM) is more effective than the traditional methods like the finite element method (FEM) when dealing with complex geophysical problems, for its high-order accuracy with exponential convergence. However, at present, little research has been done for SEM in the frequency domain, which will be more efficient than the time-domain SEM for narrowband simulations with multiple sources, material dispersion and attenuation. Herein, we systematically develop a SEM in the frequency domain to simulate coupled poroelastic, elastic and acoustic waves in anisotropic (i.e. porosity, permeability and elastic coefficients with anisotropy), heterogeneous, and lossy media. Furthermore, we completely remove the dimension inconsistency between the displacement field and the pressure in porous media to reduce the condition number of the system matrix by around 16 orders of magnitude while maintaining the symmetry of the system matrix. To solve the multiphysics coupling problems, we apply different coupling conditions to different interface types, and use basis functions to discretize the corresponding governing equations. Numerical examples show that the proposed SEM can obtain higher accuracy with much fewer unknowns compared with the FEM and has the capacity to solve the large-scale real coupling problems.