Reactive transport of CO2 in saline aquifers with implicit geomechanical analysis

Geological storage of CO2 in saline aquifers is a promising way to reduce the concentration of the greenhouse gas in the atmosphere. Injection of CO2 will, however, lead to dissolution of minerals in regions of lowered pH and precipitation of minerals from transported ions in regions of higher pH. The geomechanical implications of these changes on the stability of the reservoir are of crucial importance in the evaluation of potential injection reservoirs. The possible injection rate for given over-pressures of the injected CO2 depends on the porosity and permeability of the rock matrix in the vicinity of the injection well. Local fracturing in this region can be a tool for increasing the injection flow rate but a geomechanical analysis will be needed in order to make sure that this fracturing will not affect the geomechanical stability outside this limited region to a significant degree. This paper presents a new rewritten version of RetrasoCodeBright (RCB) which have been extended to simulations of CO2 storage in saline aquifers. An advantage of this code compared to other codes is the implicit geo-mechanical module. The code has been rewritten to account for non-ideal gas through corrections of gas density and gas solubility in all transport terms. Newton–Raphson method used to solve the flow and mechanics in RCB has been improved so as to improve convergence even under high gas injecting pressures. A 2D hydro-chemical-mechanical problem is used to illustrate the modified RCB code.

[1]  Koji Fujiwara,et al.  Method for determining relaxation factor for modified Newton-Raphson method , 1993 .

[2]  Koji Fujiwara,et al.  Improvements of convergence characteristics of Newton-Raphson method for nonlinear magnetic field analysis , 1992 .

[3]  Guglielmo Rubinacci,et al.  Numerical procedures for the solution of nonlinear electromagnetic problems , 1992 .

[4]  J. Morse,et al.  The incorporation of Mg2+ and Sr2+ into calcite overgrowths: influences of growth rate and solution composition , 1983 .

[5]  Koji Fujiwara,et al.  SUMMARY OF RESULTS FOR BENCHMARK PROBLEM 13 (3‐D NONLINEAR MAGNETOSTATIC MODEL) , 1992 .

[6]  R. Haszeldine.  Deep Geological CO2 Storage: Principles Reviewed, and Prospecting for Bio-energy Disposal Sites , 2006 .

[7]  G. Soave.  Equilibrium constants from a modified Redlich-Kwong equation of state , 1972 .

[8]  Jesús Carrera,et al.  RETRASO, a code for modeling reactive transport in saturated and unsaturated porous media , 2004 .

[9]  A. Gens,et al.  Nonisothermal multiphase flow of brine and gas through saline media , 1994 .

[10]  L. Gránásy,et al.  Phase field modeling of CH4 hydrate conversion into CO2 hydrate in the presence of liquid CO2. , 2007, Physical Chemistry, Chemical Physics - PCCP.

[11]  B. Kvamme,et al.  CO2 storage in the Utsira Formation – ATHENA 3D reactive transport simulations , 2005 .

[12]  J. Carrera,et al.  Numerical formulation for a simulator (CODE_BRIGHT) for the coupled analysis of saline media , 1996 .

[13]  T. Buanes,et al.  Storage of CO2 in natural gas hydrate reservoirs and the effect of hydrate as an extra sealing in cold aquifers , 2007 .

[14]  Joseph N. Moore,et al.  Natural CO2 Reservoirs on the Colorado Plateau and Southern Rocky Mountains, USA: A Numerical Model , 2003 .

[15]  A. Lasaga.  Chemical kinetics of water‐rock interactions , 1984 .