The migration of gases from deep to shallow reservoirs can cause damageable events. For instance, some gases can pollute the biosphere or trigger explosions and eruptions. Seismic tomography may be employed to map the accumulation of subsurface bubble‐bearing fluids to help mitigating such hazards. Nevertheless, how gas bubbles modify seismic waves is still unclear. We show that saturated rocks strongly attenuate seismic waves when gas bubbles occupy part of the pore space. Laboratory measurements of elastic wave attenuation at frequencies <100 Hz are modeled with a dynamic gas dissolution theory demonstrating that the observed frequency‐dependent attenuation is caused by wave‐induced‐gas‐exsolution‐dissolution (WIGED). This result is incorporated into a numerical model simulating the propagation of seismic waves in a subsurface domain containing CO2‐gas bubbles. This simulation shows that WIGED can significantly modify the wavefield and illustrates how accounting for this physical mechanism can potentially improve the monitoring and surveying of gas bubble‐bearing fluids in the subsurface.

Using the forced oscillation method, we measure the extensional‐mode attenuation and Young's modulus of a Berea sandstone sample at seismic frequencies (0.5–50 Hz) for varying levels of water saturation (~0–100%) and confining pressures (2–25 MPa). Attenuation is negligible for dry conditions and saturation levels <80%. For saturation levels between ~91% and ~100%, attenuation is significant and frequency dependent in the form of distinct bell‐shaped curves having their maxima between 1 and 20 Hz. Increasing saturation causes an increase of the overall attenuation magnitude and a shift of its peak to lower frequencies. On the other hand, increasing the confining pressure causes a reduction in the attenuation magnitude and a shift of its peak to higher frequencies. For saturation levels above ~98%, the fluid pressure increases with increasing confining pressure. When the fluid pressure is high enough to ensure full water saturation of the sample, attenuation becomes negligible. A second series of comparable experiments reproduces these results satisfactorily. Based on a qualitative analysis of the data, the frequency‐dependent attenuation meets the theoretical predictions of mesoscopic wave‐induced fluid flow (WIFF) in response to a heterogeneous water distribution in the pore space, so‐called patchy saturation. These results show that mesoscopic WIFF can be an important source of seismic attenuation at reservoir conditions.

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.