Finite-element numerical simulations of seismic attenuation in finely layered rocks

P-wave conversion to slow diffusion (Biot) modes at mesoscopic (small-scale) inhomogeneities in porous media is believed to be the most important attenuation mechanisms at seismic frequencies. This study considers three periodic thin layers saturated with gas, oil, and water, respectively, a realistic scenario in hydrocarbon reservoirs, and perform finite-element numerical simulations to obtain the wave velocities and quality factors along the direction perpendicular to layering. The results are validated by comparison to the Norris-Cavallini analytical solution, constituting a cross-check for both theory and numerical simulations. The approach is not restricted to partial saturation but also applies to relevant properties in reservoir geophysics, such as porosity and permeability heterogeneities. This paper considers two cases, namely, the same rock skeleton and different fluids, and the same fluid and different dry-rock properties. Unlike the two-layer case (two fluids), the results show two relaxation peaks and the agreement between numerical and analytical solutions is excellent.