A Novel Numerical Method for Geothermal Reservoirs Embedded with Fracture Networks and Parameter Optimization for Power Generation

Geothermal recovery involves a coupled thermo-hydro-mechanical (THM) process in fractured rocks. A fluid transient equilibrium equation, considering thermal conduction, convection, and heat exchange, is established. The evolution of the reservoir permeability and the variance in the fracture aperture due to a change in the stress field are derived simultaneously. THM coupling is accomplished through iterative hydromechanical and thermo-hydro processes. To overcome the difficulty of geometric discretization, a three-dimensional THM coupler model embedded with discrete fracture networks, using a zero-thickness surface and line elements to simulate fractures and injection/production wells, is established to evaluate the geothermal production. The reliability of the method is verified by a case study. Then, this method is applied to evaluate the influence of the geometric topological characteristics of fracture networks and the fracture aperture on the reservoir temperature evolution and heat extraction effectiveness. The results show that the power generation efficiency and geothermal depletion rate are significantly affected by the injection–production pressure. Injection wells and production wells with pressures higher than the initial fluid pressure in the fractures can be used to significantly increase power generation, but the consumption of geothermal energy and loss of efficiency are significant and rapid. To achieve better benefits for the geothermal recovery system, an optimization algorithm based on simultaneous perturbation stochastic approximation (SPSA) is proposed; it takes the power generation efficiency as the objective function, and the corresponding program is developed using MATLAB to optimize the position and pressure values for each production well. The results show that the heat transfer for the entire EGS reservoir becomes more uniform after optimization, and the heat transfer efficiency is greatly improved.

[1]  I. Sass,et al.  Thermo-hydro-mechanical modeling of an enhanced geothermal system in a fractured reservoir using carbon dioxide as heat transmission fluid- A sensitivity investigation , 2022, Energy.

[2]  Zetian Zhang,et al.  Migration of the Industrial Wastewater in Fractured Rock Masses Based on the Thermal-Hydraulic-Mechanical Coupled Model , 2021, Geofluids.

[3]  M. Singh,et al.  Impact of Well Placement in the Fractured Geothermal Reservoirs Based on Available Discrete Fractured System , 2021, Geosciences.

[4]  G. P. Oliveira,et al.  Poiseuille-Number-Based Kozeny–Carman Model for Computation of Pore Shape Factors on Arbitrary Cross Sections , 2021, Transport in Porous Media.

[5]  M. Meng,et al.  Critical review of wellbore ballooning and breathing literature , 2021 .

[6]  Gangwei Fan,et al.  Impacts of longwall mining speeds on permeability of weakly cemented strata and subsurface watertable: a case study , 2021, Geomatics, Natural Hazards and Risk.

[7]  Tai-lu Li,et al.  Synergetic mechanism of fracture properties and system configuration on techno-economic performance of enhanced geothermal system for power generation during life cycle , 2020 .

[8]  K. Zosseder,et al.  Porosity–permeability relationship derived from Upper Jurassic carbonate rock cores to assess the regional hydraulic matrix properties of the Malm reservoir in the South German Molasse Basin , 2020, Geothermal Energy.

[9]  Dongxiao Zhang,et al.  Coupled thermo-hydro-mechanical analysis of stimulation and production for fractured geothermal reservoirs , 2019, Applied Energy.

[10]  H. M. Nick,et al.  Towards optimisation of geothermal heat recovery: An example from the West Netherlands Basin , 2019, Applied Energy.

[11]  B. Liu,et al.  Numerical manifold method for thermal–hydraulic coupling in fractured enhance geothermal system , 2019, Engineering Analysis with Boundary Elements.

[12]  Sheng-Qi Yang,et al.  A 2D coupled hydro-thermal model for the combined finite-discrete element method , 2019 .

[13]  Zhangxin Chen,et al.  An analysis of stochastic discrete fracture networks on shale gas recovery , 2018, Journal of Petroleum Science and Engineering.

[14]  Guowei Ma,et al.  Evaluation of geothermal development in fractured hot dry rock based on three dimensional unified pipe-network method , 2018 .

[15]  X. Zhuang,et al.  Cracking elements: A self-propagating Strong Discontinuity embedded Approach for quasi-brittle fracture , 2018 .

[16]  X. Jia,et al.  Prospects of power generation from an enhanced geothermal system by water circulation through two horizontal wells: A case study in the Gonghe Basin, Qinghai Province, China , 2018 .

[17]  M. Dąbrowski,et al.  The Effective Transmissivity of a Plane‐Walled Fracture With Circular Cylindrical Obstacles , 2018 .

[18]  Pardeep Garg,et al.  Techno-economic comparison of solar organic Rankine cycle (ORC) and photovoltaic (PV) systems with energy storage , 2017 .

[19]  Wenjiong Cao,et al.  An analytical method to determine the fluid-rock heat transfer rate in two-equation thermal model for EGS heat reservoir , 2017 .

[20]  Wenbo Huang,et al.  Heat extraction performance of EGS with heterogeneous reservoir: A numerical evaluation , 2017 .

[21]  J. Latham,et al.  The use of discrete fracture networks for modelling coupled geomechanical and hydrological behaviour of fractured rocks , 2017 .

[22]  R. Xu,et al.  Convective heat transfer of supercritical CO2 in a rock fracture for enhanced geothermal systems , 2017 .

[23]  Satish Karra,et al.  Fracture size and transmissivity correlations: Implications for transport simulations in sparse three‐dimensional discrete fracture networks following a truncated power law distribution of fracture size , 2016 .

[24]  Quan Gan,et al.  Production optimization in fractured geothermal reservoirs by coupled discrete fracture network modeling , 2016 .

[25]  Hal Gurgenci,et al.  Annual performance variation of an EGS power plant using an ORC with NDDCT cooling , 2016 .

[26]  B. Faybishenko,et al.  Permeability variations within mining-induced fractured rock mass and its influence on groundwater inrush , 2016, Environmental Earth Sciences.

[27]  Satish Karra,et al.  dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport , 2015, Comput. Geosci..

[28]  C. Tsang,et al.  A study of changes in deep fractured rock permeability due to coupled hydro-mechanical effects , 2015 .

[29]  Tianfu Xu,et al.  An integrated study of fluid–rock interaction in a CO2-based enhanced geothermal system: A case study of Songliao Basin, China , 2015 .

[30]  Rainer Helmig,et al.  “Non-linearities and upscaling in porous media“ Multi-physics modeling of non-isothermal compositional flow on adaptive grids , 2014 .

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

[32]  Jiliang Chen,et al.  Designing multi-well layout for enhanced geothermal system to better exploit hot dry rock geothermal energy , 2015 .

[33]  Jinchao Xu,et al.  Well-posedness and Robust Preconditioners for the Discretized Fluid-Structure Interaction Systems , 2014, 1403.0046.

[34]  Jian Hu,et al.  Numerical simulation of heat production potential from hot dry rock by water circulating through a novel single vertical fracture at Desert Peak geothermal field , 2013 .

[35]  Yu-Chao Zeng,et al.  Numerical simulation of heat production potential from hot dry rock by water circulating through two horizontal wells at Desert Peak geothermal field , 2013 .

[36]  Fangming Jiang,et al.  A novel three-dimensional transient model for subsurface heat exchange in enhanced geothermal systems , 2013 .

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

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

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

[40]  J. Fairley.  Fracture/matrix interaction in a fracture of finite extent , 2010 .

[41]  Mustafa Versan Kok,et al.  Optimization of well placement geothermal reservoirs using artificial intelligence , 2010, Comput. Geosci..

[42]  K. Min,et al.  Chemically- and mechanically-mediated influences on the transport and mechanical characteristics of rock fractures - eScholarship , 2009 .

[43]  Gudmundur S. Bodvarsson,et al.  A modeling approach for analysis of coupled multiphase fluid flow, heat transfer, and deformation in fractured porous rock , 2002 .

[44]  Gudmundur S. Bodvarsson,et al.  Hydraulic conductivity of rock fractures , 1996 .

[45]  Stephen R. Brown,et al.  closure of rock joints , 1986 .