Geochemical modeling of CO2 storage in deep reservoirs: The Weyburn Project (Canada) case study

Abstract Geological storage is presently one of the most promising options for reducing anthropogenic emissions of CO 2 . Among the several projects investigating the fate of CO 2 stored at depth, the EnCana's CO 2 injection EOR ( Enhancing Oil Recovery ) project at Weyburn (Saskatchewan, Canada) is the most important oil production development that hosts an international monitoring project. In the Weyburn EOR Project CO 2 is used to increase recovery of heavy oil from the Midale Beds, a Mississippian reservoir consisting of shallow marine carbonate, where about 3 billions standard m 3 of supercritical CO 2 have been injected since 2000 with an injection rate of 5000 ton/day. In this work the available dataset (bulk mineralogy of the reservoir, gas-cap composition and selected pre- and post-CO 2 injection water samples) provided by the International Energy Agency Weyburn CO 2 Monitoring & Storage Project has been used in order to: i) reconstruct the pre-injection reservoir chemical composition (including pH and the boundary conditions at 62 °C and 15 MPa); ii) assess the evolution of the reservoir subjected to CO2 injection and predict dissolution/precipitation processes of the Weyburn brines over 100 years after injection; iii) validate the short-term (September 2000–2003) evolution of the in situ reservoir fluids due to the CO2 injection, by comparing the surface analytical data with the composition of the computed depressurized brines. To achieve these goals the PRHEEQC (V2.14) Software Package was used with both modified thermodynamic database and correction for supercritical CO 2 fugacity. The oil–gas–water interaction and the non-ideality of the gas phase (with exception of CO 2 ) were not considered in the numerical simulations. Despite intrinsic limitations and uncertainties of geochemical modeling, the main results can be summarized, as follows: 1) the calculated pre-injection chemical composition of the Midale Beds brine is consistent with the analytical data of the waters collected in 2000 (baseline survey), 2) the main reservoir reactions (CO 2 and carbonate dissolution) take place within the first year of simulation, 3) the temporal evolution of the chemical features of the fluids in the Weyburn reservoir suggests that CO 2 can safely be stored by solubility (as CO 2(aq) ) and mineral trapping (via dawsonite precipitation). The short-term validation performed by calculating chemical composition of the reservoir fluids (corrected for surface conditions) after the simulation of 3 years of CO 2 injection is consistent (error ≤ 5%) with the analytical data of the wellhead water samples collected in 2003, with the exception of Ca and Mg (error > 90%), likely due to complexation effect of carboxilic acid.

[1]  David J. C. Mundy,et al.  DIAGENESIS AND POROSITY DEVELOPMENT OF A SUBCROPPED MISSISSIPPIAN CARBONATE OIL RESERVOIR, AN EXAMPLE FROM THE ALIDA BEDS OF THE PHEASANT RUMP POOL, SOUTHEAST SASKATCHEWAN , 1998 .

[2]  M. Lagache.  New data on the kinetics of the dissolution of alkali feldspars at 200°C in CO2 charged water , 1976 .

[3]  D. Wesolowski,et al.  Aqueous high-temperature solubility studies. II. The solubility of boehmite at 0.03 m ionic strength as a function of temperature and pH as determined by in situ measurements , 2001 .

[4]  D. B. Stewart,et al.  The IEA Weyburn CO2 Monitoring and Storage Project , 2003 .

[5]  E. Busenberg,et al.  The dissolution kinetics of feldspars at 25°C and 1 atm CO2 partial pressure , 1976 .

[6]  B. Mayer,et al.  A model for partitioning gases among brines and hydrocarbons in oil reservoirs: Examples from the IEA-GHG Weyburn CO2 Monitoring and Storage Project, Saskatchewan, Canada , 2006 .

[7]  L. Marini,et al.  Use of reaction path modeling to identify the processes governing the generation of neutral Na-Cl and acidic Na-Cl-SO , 2003 .

[8]  M. Velbel.  Influence of Surface Area, Surface Characteristics, and Solution Composition on Feldspar Weathering Rates , 1987 .

[9]  E. Merino.  Internal consistency of a water analysis and uncertainty of the calculated distribution of aqueous species at 25°C , 1979 .

[10]  Beverly Z. Saylor,et al.  Computer simulation of CO2 trapped through mineral precipitation in the Rose Run Sandstone, Ohio , 2006 .

[11]  Yousif K. Kharaka,et al.  A Compilation of Rate Parameters of Water-Mineral Interaction Kinetics for Application to Geochemical Modeling , 2004 .

[12]  H. Feely,et al.  Origin of Gulf Coast Salt-Dome Sulphur Deposits , 1957 .

[13]  W. Giggenbach.  Geothermal gas equilibria , 1980 .

[14]  W. Gunter,et al.  Aquifer disposal of CO2-rich greenhouse gases: Extension of the time scale of experiment for CO2-sequestering reactions by geochemical modelling , 1997 .

[15]  J. Palandri,et al.  Reconstruction of in situ composition of sedimentary formation waters , 2001 .

[16]  Kenneth S. Pitzer,et al.  Thermodynamics of electrolytes. I. Theoretical basis and general equations , 1973 .

[17]  D. Savage,et al.  The effect of organic acids on the dissolution of K-feldspar under conditions relevant to burial diagenesis , 1989, Mineralogical Magazine.

[18]  J. Dandurand,et al.  Experimental study of aluminum-acetate complexing between 60 and 200°C , 1994 .

[19]  S. Hem,et al.  Influence of Acetate, Oxalate, and Citrate Anions on Precipitation of Aluminum Hydroxide , 1983 .

[20]  Karsten Pruess,et al.  Multiphase flow dynamics during CO2 disposal into saline aquifers , 2002 .

[21]  David J. Wesolowski,et al.  Aqueous high-temperature solubility studies. I. The solubility of boehmite as functions of ionic strength (to 5 molal, NaCl), temperature (100–290°C), and pH as determined by in situ measurements , 2001 .

[22]  R. Berner,et al.  Mechanism of pyroxene and amphibole weathering; II, Observations of soil grains , 1982 .

[23]  K. Pruess,et al.  Numerical simulation of CO2 disposal by mineral trapping in deep aquifers , 2004 .

[24]  D. Grandstaff.  Changes in surface area and morphology and the mechanism of forsterite dissolution , 1978 .

[25]  J. Denis,et al.  Une méthode globale d'estimation des températures des réservoirs alimentant les sources thermales. Exemple du Massif Central Français , 1981 .

[26]  G. S. Parks,et al.  Selected values of chemical thermodynamic properties , 1953 .

[27]  Z. Pang,et al.  THEORETICAL CHEMICAL THERMOMETRY ON GEOTHERMAL WATERS: PROBLEMS AND METHODS , 1998 .

[28]  Carl I. Steefel,et al.  Reactive Transport Modeling of Geologic CO{sub 2} Sequestration in Saline Aquifers: The Influence of Intra-Aquifer Shales and the Relative Effectiveness of Structural, Solubility, and Mineral Trapping During Prograde and Retrograde Sequestration , 2001 .

[29]  E. Oelkers,et al.  Experimental study of K-feldspar dissolution rates as a function of chemical affinity at 150°C and pH 9 , 1994 .

[30]  P. Fenter,et al.  Orthoclase dissolution kinetics probed by in situ X-ray reflectivity: effects of temperature, pH, and crystal orientation , 2003 .

[31]  C. Tsang,et al.  A study of caprock hydromechanical changes associated with CO2-injection into a brine formation , 2002 .

[32]  S. Brantley,et al.  Chemical weathering rates of silicate minerals , 1995 .

[33]  Albert Wegelin,et al.  GEOLOGY AND RESERVOIR PROPERTIES OF THE WEYBURN FIELD, SOUTHEASTERN SASKATCHEWAN , 1984 .

[34]  D. L. Parkhurst,et al.  User's guide to PHREEQC (Version 2)-a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations , 1999 .

[35]  Christopher A. Rochelle,et al.  The IEA Weyburn CO2 monitoring and storage project : final report of the European research team , 2005 .

[36]  Barry Freifeld,et al.  Real‐time quadrupole mass spectrometer analysis of gas in borehole fluid samples acquired using the U‐tube sampling methodology , 2006 .

[37]  Y. Kharaka,et al.  Deep Fluids in the Continents: I. Sedimentary Basins , 2003 .

[38]  P. Aagaard,et al.  Thermodynamic and kinetic constraints on reaction rates among minerals and aqueous solutions. II. Rate constants, effective surface area, and the hydrolysis of feldspar , 1984 .

[39]  John H. Weare,et al.  An equation of state for the CH4-CO2-H2O system: I. Pure systems from 0 to 1000°C and 0 to 8000 bar , 1992 .

[40]  Y. Kharaka,et al.  Organic ligand distribution and speciation in sedimentary basin brines, diagenetic fluids and related ore solutions , 1994, Geological Society, London, Special Publications.

[41]  S. Arnórsson,et al.  The chemistry of geothermal waters in Iceland. II. Mineral equilibria and independent variables controlling water compositions , 1983 .

[42]  F. Quattrocchi,et al.  Strontium isotope (87SR/86SR) chemistry in produced oilfield waters: The IEA Weyburn CO2 Monitoring and Storage Project , 2006 .

[43]  G. R. Holdren,et al.  Mechanism of feldspar weathering—II. Observations of feldspars from soils , 1979 .

[44]  Nicolas Spycher,et al.  Calculation of pH and mineral equilibria in hydrothermal waters with application to geothermometry and studies of boiling and dilution , 1984 .

[45]  R. A. Robie,et al.  The entropy and Gibbs free energy of formation of the aluminum ion , 1977 .

[46]  K. Bateman,et al.  The Underground Disposal of Carbon Dioxide: Final Report , 1996 .

[47]  H. L. Miller,et al.  Climate Change 2007: The Physical Science Basis , 2007 .

[48]  David L. Parkhurst,et al.  Revised chemical equilibrium data for major water-mineral reactions and their limitations , 1990 .

[49]  L. Evans,et al.  DISSOLUTION OF FELDSPARS BY LOW‐MOLECULAR-WEIGHT ALIPHATIC AND AROMATIC ACIDS , 1986 .

[50]  Luigi Marini,et al.  Geological Sequestration of Carbon Dioxide: Thermodynamics, Kinetics, and Reaction Path Modeling , 2006 .

[51]  S. Holloway,et al.  The underground disposal of carbon dioxide : summary report , 1996 .

[52]  P. Fenter,et al.  Resolving orthoclase dissolution processes with atomic force microscopy and X-ray reflectivity , 2001 .

[53]  Susan D. Hovorka,et al.  The U-Tube: A Novel System for Acquiring Borehole Fluid Samples from a Deep Geologic CO2 Sequestration Experiment , 2005 .

[54]  C. R. Evans,et al.  Alteration of Crude Oil by Waters and Bacteria--Evidence from Geochemical and Isotope Studies , 1973 .

[55]  Kathryn L. Nagy,et al.  Chemical weathering rate laws and global geochemical cycles , 1994 .

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

[57]  Karsten Pruess,et al.  Reactive geochemical transport simulation to study mineral trapping for CO2 disposal in deep arenaceous formations , 2003 .

[58]  O. Pokrovsky,et al.  Experimental determination of the effect of dissolved CO2 on the dissolution kinetics of Mg and Ca silicates at 25 °C , 2005 .

[59]  William D. Gunter,et al.  Aquifer disposal of CO2-rich gases: Reaction design for added capacity , 1993 .

[60]  G. R. Holdren,et al.  Reaction rate-surface area relationships during the early stages of weathering. II. Data on eight additional feldspars , 1987 .

[61]  Karsten Pruess,et al.  Numerical modeling of injection and mineral trapping of CO2 with H2S and SO2 in a Sandstone Formation , 2007 .

[62]  R. C. Weast.  CRC Handbook of Chemistry and Physics , 1973 .

[63]  T. Pačes.  The kinetics of base cation release due to chemical weathering , 1990 .

[64]  C. Steefel,et al.  A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution rea , 1994 .

[65]  W. Gunter,et al.  Aquifer disposal of acid gases: modelling of water–rock reactions for trapping of acid wastes , 2000 .

[66]  Nicolas Spycher,et al.  Fugacity coefficients of H2, CO2, CH4, H2O and of H2O- CO2-CH4 mixtures: A virial equation treatment for moderate pressures and temperatures applicable to calculations of hydrothermal boiling , 1988 .

[67]  H. Barnes,et al.  Oxidation of pyrite in low temperature acidic solutions: Rate laws and surface textures , 1986 .

[68]  L. M. Walter,et al.  Kinetics of feldspar and quartz dissolution at 70–80°C and near-neutral pH: effects of organic acids and NaCl , 1999 .

[69]  Wolery,et al.  EQ3NR: a computer program for geochemical aqueous speciation-solubility calculations. User`s guide and documentation , 1983 .

[70]  Zhenhao Duan,et al.  An improved model calculating CO2 solubility in pure water and aqueous NaCl solutions from 273 to 533 K and from 0 to 2000 bar , 2003 .

[71]  Zhenhao Duan,et al.  An improved model for the calculation of CO2 solubility in aqueous solutions containing Na+, K+, Ca2+, Mg2+, Cl−, and SO42− , 2006 .

[72]  Feasibility study ( I stage ) of CO 2 geological storage by ECBM techniques in the Sulcis coal province ( SW Sardinia , Italy ) , 2006 .

[73]  D. D. Wagman,et al.  Selected Values of Chemical Thermodynamic Properties. Tables for the First Thirty-Four Elements in the Standard Order of Arrangement. , 1968 .

[74]  W. Giggenbach.  Geothermal solute equilibria. Derivation of Na-K-Mg-Ca geoindicators , 1988 .

[75]  James W. Ball,et al.  User's manual for WATEQ4F, with revised thermodynamic data base and text cases for calculating speciation of major, trace, and redox elements in natural waters , 1991 .

[76]  L. Marini,et al.  Fluid geochemistry of the Acqui Terme-Visone geothermal area (Piemonte, Italy) , 2000 .

[77]  Ernie Perkins,et al.  Monitoring of fluid–rock interaction and CO2 storage through produced fluid sampling at the Weyburn CO2-injection enhanced oil recovery site, Saskatchewan, Canada , 2005 .

[78]  T. J. Wolery,et al.  EQ3NR, a computer program for geochemical aqueous speciation-solubility calculations: Theoretical manual, user`s guide, and related documentation (Version 7.0); Part 3 , 1992 .

[79]  P. Barak,et al.  Solubility and Dissolution Kinetics of Dolomite in Ca–Mg–HCO3/CO3 Solutions at 25°C and 0.1 MPa Carbon Dioxide , 2000 .

[80]  Christophe Kervévan,et al.  Improvement of the Calculation Accuracy of Acid Gas Solubility in Deep Reservoir Brines: Application to the Geological Storage of CO2 , 2005 .

[81]  A. Blum,et al.  Feldspar dissolution kinetics , 1995 .

[82]  H. Heasler,et al.  Organic-Inorganic Interactions and Sandstone Diagenesis , 1989 .

[83]  David R. Cole,et al.  Gas-water-rock interactions in Frio Formation following CO2 injection: Implications for the storage of greenhouse gases in sedimentary basins , 2006 .

[84]  Ernie Perkins,et al.  Geochemical monitoring of fluid-rock interaction and CO2 storage at the Weyburn CO2-injection enhanced oil recovery site, Saskatchewan, Canada , 2004 .

[85]  Z. Aizenshtat,et al.  Mechanisms of sulfate removal from subsurface calcium chloride brines: Heletz-Kokhav oilfields, Israel , 1995 .

[86]  I. Gaus,et al.  Reactive transport modelling of the impact of CO2 injection on the clayey cap rock at Sleipner (North Sea) , 2005 .

[87]  R. A. Robie,et al.  Thermodynamic properties of minerals and related substances at 298.15 K and 1 bar (10[5] pascals) pressure and at higher temperatures , 1995 .

[88]  Joseph N. Moore,et al.  Mineralogical and geochemical consequences of the long-term presence of CO2 in natural reservoirs: An example from the Springerville-St. Johns Field, Arizona, and New Mexico, U.S.A , 2005 .

[89]  Koorosh Asghari,et al.  CO2 storage during EOR process in the Weyburn Oil Pool: Modelling and monitoring results , 2005 .

[90]  L. M. Walter,et al.  Effects of organic acids on the dissolution of orthoclase at 80°C and pH 6 , 1996 .

[91]  K. Brown,et al.  Weyburn CO 2 monitoring and storage project , 2005 .

[92]  Sally M. Benson,et al.  The role of hydrogeological and geochemical trapping in sedimentary basins for secure geological storage of carbon dioxide , 2004, Geological Society, London, Special Publications.

[93]  J. Seewald.  Model for the origin of carboxylic acids in basinal brines , 2001 .

[94]  W. Giggenbach.  Geothermal mineral equilibria , 1981 .

[95]  E. Oelkers,et al.  Are quartz dissolution rates proportional to B.E.T. surface areas , 2001 .

[96]  V. A. Alekseyev,et al.  Change in the dissolution rates of alkali feldspars as a result of secondary mineral precipitation and approach to equilibrium , 1997 .

[97]  Fabrizio Gherardi,et al.  Numerical modeling of self-limiting and self-enhancing caprock alteration induced by CO2 storage in a depleted gas reservoir , 2007 .

[98]  J. B. Rapp,et al.  Short chain aliphatic acid anions in oil field waters and their contribution to the measured alkalinity , 1975 .

[99]  Lawrence M. Anovitz,et al.  Dawsonite synthesis and reevaluation of its thermodynamic properties from solubility measurements: Implications for mineral trapping of CO2 , 2007 .

[100]  B. Metz.  IPCC special report on carbon dioxide capture and storage , 2005 .

[101]  M. Wilson.  CHEMICAL WEATHERING OF SOME PRIMARY ROCK‐FORMING MINERALS , 1975 .

[102]  P. Aharon,et al.  Microbial sulfate reduction rates and sulfur and oxygen isotope fractionations at oil and gas seeps in deepwater Gulf of Mexico , 2000 .

[103]  A. Navarre‐Sitchler,et al.  Basalt weathering across scales , 2007 .

[104]  A. Faaij,et al.  Health, Safety and Environmental Risks of Underground Co2 Storage – Overview of Mechanisms and Current Knowledge , 2006 .

[105]  D. D. Wagman,et al.  Selected Values of Chemical Thermodynamic Properties. Tables for the Lanthanide (Rare Earth) Elements (Elements 62 through 76 in the Standard Order of Arrangement). , 1971 .