NUMERICAL SIMULATION OF SEISMIC WAVES IN POROUS MEDIA

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 formulas 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, identifying 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.