This work presents a mathematical algorithm for modeling the propagation of poroelastic waves. We have shown how the classical Biot equations can be put into Ursin’s form in a plane-layered 3D porous medium. Using this form, we have derived explicit for- mulas that can be used as the basis of an efficient computational algorithm. To validate the algorithm, numerical simulations were performed using both the poroelastic and equivalent elastic models. The results obtained confirmed the proposed algorithm’s reliability, identify- ing the main wave events in both low-frequency and high-frequency regimes in the reservoir and laboratory scales, respectively. We have also illustrated the influence of some physical parameters on the attenuation and dispersion of the slow wave.
Numerical simulation of coseismic electromagnetic fields associated with seismic waves due to finite faulting in porous media
