Dissolution and deformation in fractured carbonates caused by flow of CO2-rich brine under reservoir conditions

Abstract Geometrical alteration in fractures, caused by mineral dissolution, changes the contact area between fracture surfaces and affects the mechanical strength of fractures. It is difficult to determine the influence of dissolution on fracture porosity and permeability given the competition between fracture opening due to dissolution and fracture closure caused by mechanical deformation. Therefore, simulating flow in fractured reservoirs during enhanced oil recovery and CO 2 sequestration, where local porosity changes may significantly alter permeability, remains a fundamental challenge. Here, we present results from experiments and numerical simulations that explore the influence of coupled geochemical alteration and mechanical deformation on calcium carbonate fracture geometry. We scan the fracture surfaces, before and after the flow experiments, using high-resolution optical profilometry to map changes in aperture fields. Flow of brine equilibrated with CO 2 at 60 °C and pore pressure of 14 MPa leads to significant dissolution in fractured calcium carbonate cores. The evolution of the dissolution process depends, to first order, on the dimensionless Damkohler number D a . We vary D a in experiments by changing the flow rate through fractured cores and observe a transition in dissolution behavior over the range of flow rates from 0.1 ml/min to 20 ml/min. At large D a (small flow rate), dissolution causes the formation of large-scale channels aligned with the fractures. At small D a (large flow rate), dissolution occurs more uniformly over the fracture surfaces. However, the area of contacting asperities, controlled by the significant mechanical stresses (28 MPa) and fracture surface roughness, constrains the spatial extent of the dissolution. We use the measured aperture fields as input to a single-species, reactive-transport model. Simulated Ca 2+ dissolution agrees with temporal evolution of measured Ca 2+ concentrations in outflow fluid. Comparison of the simulated dissolved aperture fields with measured aperture fields after flow-through shows qualitative agreement. We also use our reactive-transport model to explore the role of changing fluid properties and flow rates on dissolution and the resulting alteration of porosity and permeability beyond our experimental conditions. Our results emphasize the importance of length-scales in the coupled response of chemo-mechanical deformation in fractured reservoirs under varied stress and flow conditions.

[1]  O. Pokrovsky,et al.  Calcite, dolomite and magnesite dissolution kinetics in aqueous solutions at acid to circumneutral pH, 25 to 150 °C and 1 to 55 atm pCO2: New constraints on CO2 sequestration in sedimentary basins , 2009 .

[2]  C. Noiriel,et al.  Hydraulic Properties and Microgeometry Evolution Accompanying Limestone Dissolution by Acidic Water , 2005 .

[3]  David J. Brush,et al.  Fluid flow in synthetic rough‐walled fractures: Navier‐Stokes, Stokes, and local cubic law simulations , 2003 .

[4]  F. Boadu,et al.  Effects of fractures on seismic-wave velocity and attenuation , 1996 .

[5]  E. L. Sjöberg,et al.  A fundamental equation for calcite dissolution kinetics , 1976 .

[6]  Etienne Brosse,et al.  Modelling Fluid-Rock Interaction Induced by the Percolation of CO2-Enriched Solutions in Core Samples: the Role of Reactive Surface Area , 2005 .

[7]  M. Ghiaasiaan.  Convective Heat and Mass Transfer: Index , 2011 .

[8]  S. Vialle,et al.  Laboratory measurements of elastic properties of carbonate rocks during injection of reactive CO2‐saturated water , 2011 .

[9]  D. Elsworth,et al.  A numerical model simulating reactive transport and evolution of fracture permeability , 2006 .

[10]  V. Balakotaiah,et al.  Effect of medium heterogeneities on reactive dissolution of carbonates , 2009 .

[11]  Harihar Rajaram,et al.  Saturated flow in a single fracture: evaluation of the Reynolds Equation in measured aperture fields , 1999 .

[12]  A. Hill,et al.  Numerical and experimental investigation of acid wormholing during acidization of vuggy carbonate rocks , 2010 .

[13]  Pierre M. Adler,et al.  Dissolution and deposition in fractures , 1997 .

[14]  R. Detwiler.  Experimental observations of deformation caused by mineral dissolution in variable-aperture fractures , 2008 .

[15]  Laura J. Pyrak-Nolte,et al.  Single fractures under normal stress: The relation between fracture specific stiffness and fluid flow , 2000 .

[16]  D. Lockner,et al.  Reduction of Permeability in Granite at Elevated Temperatures , 1994, Science.

[17]  Pierre M. Adler,et al.  Fractures and Fracture Networks , 1999 .

[18]  B. Berkowitz,et al.  Measurement and analysis of dissolution patterns in rock fractures , 2002 .

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

[20]  Jishan Liu,et al.  A fully-coupled hydrological-mechanical-chemical model for fracture sealing and preferential opening , 2006 .

[21]  D. Rickard,et al.  Temperature dependence of calcite dissolution kinetics between 1 and 62°C at pH 2.7 to 8.4 in aqueous solutions , 1984 .

[22]  H. S. Fogler,et al.  Pore evolution and channel formation during flow and reaction in porous media , 1988 .

[23]  C.A.J. Appelo,et al.  User's guide to PHREEQC - a computer program for speciation,batch-reaction, one-dimensional transport, and inversegeochemical calculations. , 1999 .

[24]  Fred Kofi Boadu,et al.  Predicting the transport properties of fractured rocks from seismic information: numerical experiments , 2000 .

[25]  E. Oelkers,et al.  An experimental study of calcite and limestone dissolution rates as a function of pH from −1 to 3 and temperature from 25 to 80°C , 1998 .

[26]  R. Detwiler,et al.  Predicting dissolution patterns in variable aperture fractures: Evaluation of an enhanced depth‐averaged computational model , 2007 .

[27]  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 .

[28]  Steffen Mehl,et al.  Comparison of an Algebraic Multigrid Algorithm to Two Iterative Solvers Used for Modeling Ground Water Flow and Transport , 2002, Ground water.

[29]  J. C. Santamarina,et al.  CO2 geological storage — Geotechnical implications , 2011 .

[30]  A. Grader,et al.  Spontaneous switching of permeability changes in a limestone fracture with net dissolution , 2002 .

[31]  Pathegama Gamage Ranjith,et al.  A review of studies on CO2 sequestration and caprock integrity , 2010 .

[32]  E. Burton,et al.  Direct observation of reactive flow in a single fracture , 2001 .

[33]  Rick Chalaturnyk,et al.  IEA GHG Weyburn CO2 monitoring and storage project , 2005 .

[34]  Daccord.  Chemical dissolution of a porous medium by a reactive fluid. , 1987, Physical review letters.

[35]  A. Ladd,et al.  The initial stages of cave formation: Beyond the one-dimensional paradigm , 2010, 1011.4111.

[36]  Michel Quintard,et al.  On the ability of a Darcy-scale model to capture wormhole formation during the dissolution of a porous medium , 2002, Journal of Fluid Mechanics.