Physics and Seismic Modeling for Monitoring CO2 Storage

We present a new petro-elastical and numerical-simulation methodology to compute synthetic seismograms for reservoirs subject to CO2 sequestration. The petro-elastical equations model the seismic properties of reservoir rocks saturated with CO2, methane, oil and brine. The gas properties are obtained from the van der Waals equation and we take into account the absorption of gas by oil and brine, as a function of the in situ pore pressure and temperature. The dry-rock bulk and shear moduli can be obtained either by calibration from real data or by using rock-physics models based on the Hertz-Mindlin and Hashin-Shtrikman theories. Mesoscopic attenuation due to fluids effects is quantified by using White's model of patchy saturation, and the wet-rock velocities are calculated with Gassmann equations by using an effective fluid modulus to describe the velocities predicted by White's model. The simulations are performed with a poro-viscoelastic modeling code based on Biot's theory, where viscoelasticity is described by generalizing the solid/fluid coupling modulus to a relaxation function. Using the pseudo-spectral method, which allows general material variability, a complete and accurate characterization of the reservoir can be obtained. A simulation, that considers the Utsira sand of the North Sea, illustrates the methodology.