Seismic Wave Attenuation and Dispersion Due to Partial Fluid Saturation: Direct Measurements and Numerical Simulations Based on X-Ray CT

Quantitatively assessing seismic attenuation caused by fluid pressure diffusion (FPD) in partially saturated rocks is challenging because of its sensitivity to the spatial fluid distribution. To address this challenge we performed depressurization experiments to induce the exsolution of carbon dioxide from water in a Berea sandstone sample. In a first set of experiments we used medical X-ray computed tomography (CT) to characterize the fluid distribution. At an equilibrium pressure of approximately 1 MPa and applying a fluid pressure decline rate of approximately 0.6 MPa per minute, we allowed a change in saturation of less than 1%. The gas was heterogeneously distributed along the length of the sample, with most of the gas exsolving near the sample outlet. In a second set of experiments, at the same pressure and temperature, following a very similar exsolution protocol, we measured the frequency dependent attenuation and modulus dispersion between 0.1 and 1,000 Hz using the forced oscillation method. We observed significant attenuation and dispersion in the extensional and bulk deformation modes, however, not in the shear mode. Lastly, we use the fluid distribution derived from the X-ray CT as an input for numerical simulations of FPD to compute the attenuation and modulus dispersion. The numerical solutions are in close agreement with the attenuation and modulus dispersion measured in the laboratory. Our approach allows for accurately relating attenuation and dispersion to the fluid distribution, which can be applied to improving the seismic monitoring of the subsurface.