A numerical procedure to model and monitor CO2 sequestration in aquifers

Carbon Dioxide (CO2) sequestration into geologic formations is a means of mitigating greenhouse effect. In this work we present a new numerical simulation technique to model and monitor CO2 sequestration in aquifers. For that purpose we integrate numerical simulators of CO2-brine flow and seismic wave propagation (time-lapse seismics). The simultaneous flow of brine and CO2 is modeled applying the Black-Oil formulation for two phase flow in porous media, which uses the Pressure-Volume-Temperature (PVT) behavior as a simplified thermodynamic model. Seismic wave propagation uses a simulator based on a space-frequency domain formulation of the viscoelastic wave equation. In this formulation, the complex and frequency dependent coefficients represent the attenuation and dispersion effect suffered by seismic waves travelling in fluid-saturated heterogeneous porous formations. The spatial discretization is achieved employing a nonconforming finite element space to represent the displacement vector. Numerical examples of CO2 injection and time-lapse seismics in the Utsira formation at the Sleipner field are analyzed. The Utsira formation is represented using a new petrophysical model that allows a realistic inclusion of shale seals and fractures. The results of the simulations show the capability of the proposed methodology to monitor the spatial distribution of CO2 after injection.