Multiphase flow dynamics during CO2 disposal into saline aquifers

Injection Of CO2 into saline aquifers is described by mass conservation equations for the three components water, salt (NaCl), and CO2. The equations are discretized using an integral finite difference method, and are solved using methods developed in geothermal and petroleum reservoir engineering. Phase change processes are treated through switching of primary thermodynamic variables. A realistic treatment of PVT (fluid) properties is given which includes salinity and fugacity effects for partitioning of CO2 between gaseous and aqueous phases. Chemical reactions and mechanical stress effects are neglected. Numerical simulations are presented for injection Of CO2 into a brine aquifer, and for loss Of CO2 from storage through discharge along a fault zone. it is found that simulated pressures are much more sensitive to space discretization effects than are phase saturations. CO2 discharge along a fault is a self-enhancing process whose flow rates can increase over time by more than an order of magnitude, suggesting that reliable containment Of CO2 will require multiple barriers.