A coupled thermo-hydro-mechanical modeling of fracture aperture alteration and reservoir deformation during heat extraction from a geothermal reservoir

Abstract Hot water extraction and cold water injection into an underground geothermal reservoir cause mechanical deformation of rock matrix and rock joints/fractures. That leads to alteration of hydraulic transmissivity. To study the evolution of reservoir transmissivity we performed coupled Thermo-Hydro-Mechanical (THM) simulations using a robust code called Finite Element for Heat and Mass Transfer (FEHM) for a 3-D domain with a single fracture connecting the injection and production wells. Rock fracture was modeled as a thin equivalent porous medium. We established dynamic relations between the properties of the equivalent porous medium and fracture aperture. In this paper we discuss the alteration of fracture aperture due to heat extraction. The channeling of flow between injection and production wells by THM effects causes faster temperature drawdown and reduces energy production. The model also predicted fracture opening near injection well and closure at far field locations. We also simulated the aperture alteration for different joint stiffness, thermal expansion coefficients and rock matrix permeabilities. Increase in rock matrix permeability not only causes the leakage of injected water but also increases matrix contraction due to cooling and therefore the aperture growth. Additionally we reported the effect of thermo-poro-elastic deformation on the expansion and contraction of the formation for different reservoir properties. We established that in the early-stages the compaction/expansion of the formation was controlled by pore pressure change but in the late-stage it was controlled by thermal contraction.

[1]  R. Bertani.  Geothermal power generation in the world 2010–2014 update report , 2016 .

[2]  R. J. Pine,et al.  A hydro-thermo-mechanical numerical model for HDR geothermal reservoir evaluation , 1996 .

[3]  Enrique Merino,et al.  Geochemical self-organization I; reaction-transport feedbacks and modeling approach , 1987 .

[4]  Thomas Kohl,et al.  Coupled hydraulic, thermal and mechanical considerations for the simulation of hot dry rock reservoirs , 1995 .

[5]  Ghazal Izadi,et al.  The influence of thermal-hydraulic-mechanical- and chemical effects on the evolution of permeability, seismicity and heat production in geothermal reservoirs , 2015 .

[6]  Karsten Pruess,et al.  Thermal–hydrodynamic–chemical (THC) modeling based on geothermal field data , 2004 .

[7]  Maurizio Vaccaro,et al.  Numerical simulation of geothermal reservoirs for the sustainable design of energy plants: A review , 2014 .

[8]  Nasser Khalili,et al.  A thermo-hydro-mechanical coupled model in local thermal non-equilibrium for fractured HDR reservoir with double porosity , 2012 .

[9]  A. Ghassemi,et al.  A three-dimensional thermo-poroelastic model for fracture response to injection/extraction in enhanc , 2011 .

[10]  Yangsheng Zhao,et al.  THM (Thermo-hydro-mechanical) coupled mathematical model of fractured media and numerical simulation of a 3D enhanced geothermal system at 573 K and buried depth 6000–7000 M , 2015 .

[11]  M. Saar,et al.  Combining geothermal energy capture with geologic carbon dioxide sequestration , 2011 .

[12]  Ahmad Ghassemi,et al.  Changes in fracture aperture and fluid pressure due to thermal stress and silica dissolution/precipitation induced by heat extraction from subsurface rocks , 2007 .

[13]  Christopher C. Pain,et al.  Non-linear regimes of fluid flow in rock fractures , 2004 .

[14]  Mario J. Martinez,et al.  Coupled multiphase flow and geomechanics model for analysis of joint reactivation during CO2 sequestration operations , 2013 .

[15]  Thomas Stephen Lowry,et al.  Modeling of heat extraction from variably fractured porous media in Enhanced Geothermal Systems , 2016 .

[16]  Darius Mottaghy,et al.  Modeling contribution to risk assessment of thermal production power for geothermal reservoirs , 2013 .

[17]  J. Carrera,et al.  Thermal coupling may control mechanical stability of geothermal reservoirs during cold water injection , 2013 .

[18]  N. Barton,et al.  FUNDAMENTALS OF ROCK JOINT DEFORMATION , 1983 .

[19]  Sheik S. Rahman,et al.  Numerical simulation of Fluid-Rock coupling heat transfer in naturally fractured geothermal system , 2011 .

[20]  Heat transport through rough channels , 2003, cond-mat/0311391.

[21]  O. Sarychikhina,et al.  Spatio-temporal evolution of aseismic ground deformation in the Mexicali Valley (Baja California, Mexico) from 1993 to 2010, using differential SAR interferometry , 2015 .

[22]  Lehua Pan,et al.  Fully coupled wellbore-reservoir modeling of geothermal heat extraction using CO2 as the working fluid , 2015 .

[23]  J. Schmittbuhl,et al.  Fracture roughness and thermal exchange: A case study at Soultz-sous-Forêts , 2010 .

[24]  Rajesh J. Pawar,et al.  Modeling caprock bending stresses and their potential for induced seismicity during CO2 injection , 2014 .

[25]  Harihar Rajaram,et al.  Investigation of permeability alteration of fractured limestone reservoir due to geothermal heat extraction using three-dimensional thermo-hydro-chemical (THC) model , 2014 .

[26]  Thomas Kohl,et al.  Coupled thermal–hydraulic–chemical modelling of enhanced geothermal systems , 2005 .

[27]  Jens Birkholzer,et al.  Coupled thermo-hydro-mechanical processes in the near field of a high-level radioactive waste repository in clay formations , 2012 .

[28]  J. Vrugt,et al.  Inverse Modeling of Subsurface Flow and Transport Properties: A Review with New Developments , 2008 .

[29]  M. Saar,et al.  Brine displacement by CO2, energy extraction rates, and lifespan of a CO2-limited CO2-Plume Geothermal (CPG) system with a horizontal production well , 2015 .

[30]  Darius Mottaghy,et al.  The geothermal project Den Haag: 3D numerical models for temperature prediction and reservoir simula , 2011 .

[31]  Ki-Bok Min,et al.  Determination of optimum parameters of doublet system in a horizontally fractured geothermal reservoir , 2012 .

[32]  Harihar Rajaram,et al.  Fracture transmissivity evolution due to silica dissolution/precipitation during geothermal heat extraction , 2015 .

[33]  Michael A. Hicks,et al.  A prototype design model for deep low-enthalpy hydrothermal systems , 2015 .

[34]  Ahmad Ghassemi,et al.  A three‐dimensional integral equation model for calculating poro‐ and thermoelastic stresses induced by cold water injection into a geothermal reservoir , 2009, International Journal for Numerical and Analytical Methods in Geomechanics.

[35]  Seyyed A. Hosseini,et al.  Hydro-thermo-mechanical analysis during injection of cold fluid into a geologic formation , 2015 .

[36]  Daniel Sutter,et al.  Sustainable heat farming: Modeling extraction and recovery in discretely fractured geothermal reservoirs , 2013 .

[37]  Mei Ding,et al.  Modeling solute transport through saturated zone ground water at 10 km scale: example from the Yucca Mountain license application. , 2010, Journal of contaminant hydrology.

[38]  Christian Vogt,et al.  Maximum potential for geothermal power in Germany based on engineered geothermal systems , 2015, Geothermal Energy.

[39]  Joshua Taron,et al.  Thermal–hydrologic–mechanical–chemical processes in the evolution of engineered geothermal reservoirs , 2009 .

[40]  Olaf Kolditz,et al.  A comparative simulation study of coupled THM processes and their effect on fractured rock permeability around nuclear waste repositories , 2009 .

[41]  A. Cheng,et al.  Effects of heat extraction on fracture aperture: A poro–thermoelastic analysis , 2008 .

[42]  J. Hansen,et al.  The impact of the 2009–10 El Niño Modoki on U.S. West Coast beaches , 2011 .

[43]  P. N. Razdan,et al.  Geothermal Energy Resources and its Potential in India , 2008 .

[44]  S. Rahman,et al.  A numerical study on the long term thermo-poroelastic effects of cold water injection into naturally fractured geothermal reservoirs , 2011 .

[45]  H. Viswanathan,et al.  Alteration of fractures by precipitation and dissolution in gradient reaction environments: Computational results and stochastic analysis , 2008 .

[46]  Satish Karra,et al.  A simulator for modeling coupled thermo-hydro-mechanical processes in subsurface geological media , 2014 .

[47]  Ahmad Ghassemi,et al.  A reactive thermo-poroelastic analysis of water injection into an enhanced geothermal reservoir , 2014 .

[48]  Y. Hao,et al.  A geomechanical mechanism that counteracts flow channeling induced by reservoir thermal drawdown , 2013 .

[49]  Roy Baria,et al.  HDR/HWR reservoirs: concepts, understanding and creation , 1999 .

[50]  D. Bruel.  Impact of Induced Thermal Stresses During Circulation Tests in an Engineered Fractured Geothermal Reservoir: Example of the Soultz-Sous-Forêts European Hot Fractured Rock Geothermal Project, Rhine Graben, France , 2002 .

[51]  Sergey V. Samsonov,et al.  Ground deformation in the Taupo Volcanic Zone, New Zealand, observed by ALOS PALSAR interferometry , 2011 .

[52]  Pengcheng Fu,et al.  Thermal drawdown-induced flow channeling in a single fracture in EGS , 2016 .

[53]  H. Viswanathan,et al.  Early‐stage hypogene karstification in a mountain hydrologic system: A coupled thermohydrochemical model incorporating buoyant convection , 2013 .

[54]  Harihar Rajaram,et al.  Influence of aperture variability on dissolutional growth of fissures in Karst Formations , 1998 .