This technique uses a mixed operator length to approximate spatial derivatives to improve accuracy and stability. The numerical solutions closely resemble the analytical solutions and all waves predicted by Biot’s theory are generated. Application to 2D homogenous model, two-layered model including substantial amounts of CO2 sequestration and 3D poroelastic model has demonstrated the efficacy of the technique with reduced computational time.
E. Anthony and N. Vedanti (2022) Geophysics vol. 87(3), 225-236, https://doi.org/10.1190/geo2020-0240.1