A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage

CO₂ injection and storage in deep saline aquifers involves many coupled processes, including multiphase flow, heat and mass transport, rock deformation and mineral precipitation and dissolution. Coupling is especially critical in carbonate aquifers, where minerals will tend to dissolve in response to the dissolution of CO₂ into the brine. The resulting neutralization will drive further dissolution of both CO₂ and calcite. This suggests that large cavities may be formed and that proper simulation may require full coupling of reactive transport and multiphase flow. We show that solving the latter may suffice whenever two requirements are met: (1) all reactions can be assumed to occur in equilibrium and (2) the chemical system can be calculated as a function of the state variables of the multiphase flow model (i.e., liquid and gas pressure, and temperature). We redefine the components of multiphase flow codes (traditionally, water and CO₂), so that they are conservative for all reactions of the chemical system. This requires modifying the traditional constitutive relationships of the multiphase flow codes, but yields the concentrations of all species and all reaction rates by simply performing speciation and mass balance calculations at the end of each time step. We applied this method to the H₂O–CO₂–Na–Cl–CaCO₃ system, so as to model CO₂ injection into a carbonate aquifer containing brine. Results were very similar to those obtained with traditional formulations, which implies that full coupling of reactive transport and multi-phase flow is not really needed for this kind of systems, but the resulting simplifications may make it advisable even for cases where the above requirements are not met. Regarding the behavior of carbonate rocks, we find that porosity development near the injection well is small because of the low solubility of calcite. Moreover, dissolution concentrates at the front of the advancing CO₂ plume because the brine below the plume tends to reach high CO₂ concentrations quite rapidly. We conclude that carbonate dissolution needs not to be feared.