In the long-term mining of geothermal resources in hot dry rock (HDR), the change of thermal stress and pore pressure will increase fracture conductivity evolution, further improving production performance. The optimization and decision-making of the development scheme based on the impact of damage from fractures have yet to be reported. The damage to fractures is essential in designing and adjusting geothermal resource development schemes, particularly in selecting optimal schemes. Therefore, the production performances of HDR resources under different parameters are analyzed to establish a database. Then, minimizing flow resistance, maximizing net power, and maximizing economic benefits are set as optimization goals. Various injection-mining parameters and fracture characteristics are treated as decision variables. Multi-objective optimization and multi-attribute decision analysis is conducted to obtain optimal schemes. Finally, optimal schemes are evaluated and compared, considering damage and non-damage scenarios. Results show that the NSGA-II algorithm is more suitable for optimizing geothermal development questions. Net power and economic benefits of the optimal scheme considering damage increase by 45.84 % and 21.35 % compared to the control scheme with damage. For the non-damage scenario, the above values increased by 31.55 % and 5.15 %, respectively. Compared to not considering the damage, higher mass flow and well spacing of optimal scheme can be selected for the case when damaged. Moreover, the parametric design of the optimal scheme becomes more conservative as the production cycle increases.

@en

The development of hot dry rocks (HDRs) is of great significance to adjusting energy structure, alleviating energy shortage, reducing pollution, etc. Low-permeability granite is the predominant rock type in deep HDRs, making fractures the primary pathways for fluid circulation and heat extraction. The production of HDRs is significantly influenced by variable fracture conductivity, but current conductivity characterization primarily relies on the elastic deformation of the matrix, neglecting the impact of damage. Accordingly, we propose an experimental method and a supporting apparatus, which is used to unveil the conductivity evolution characteristics resulting from the comprehensive effects of damage and elastic deformation. The experimental results demonstrate that when subjected to confining force squeezing inward, the fracture conductivity experiences varying degrees of decrease compared to its initial state before the experiment. By utilizing the conductivity evolution rate as the evaluation criterion and conducting grey correlation analysis, it has been determined that temperature exerts the most significant influence on the conductivity evolution, followed by injection flow, and lastly, confining pressure. Moreover, rock particle types and production cycles also have different degrees of effect. After considering the comprehensive effects of damage-elastic deformation at the field-scale, the damage has a positive effect on conductivity enhancement. Our study provides a new approach for the characterization of fracture conductivity evolution for deep geothermal projects.

@en

Geothermal is an important renewable energy source, but the high cost of drilling limits its popularization. The oilfield area is rich in geothermal resources and has a large number of high-temperature abandoned wells. It is an economic and effective method to transform the abandoned wells into geothermal wells. At present, the heat extraction research of abandoned wells mostly focuses on single-well closed systems, while the most common in oilfields is the well patterns open system, which involves the flow and heat extraction of oil-water multiphase. More, the main utilization mode of oilfield geothermal is heating, which is an intermittent operation, and the existence of the heat recovery period makes the solution process more complicated compared to continuous production. Therefore, the numerical simulation of oil-water two-phase and the scheme optimization under intermittent operation are the key problems in the current oilfields geothermal research, that is, it should be made clear what the effects of operation parameters are and how to get the optimization design for the intermittent mining of heat-oil cogeneration from abandoned wells. For this reason, the oil-water two-phase heat-flow coupling model and economic evaluation model are established. Then, database is established through parameter sensitivity analysis, and multi-objective optimization research is carried out on the lifetime and economic benefits using Non Sorting Genetic Algorithm II (NSGA-II). The results show that economic benefit and end-point temperature of the optimized scheme are enhanced by 6.75 million US dollars and 7.10 K after 15-year production, respectively. The heat-oil cogeneration performance is superior, so the influence of oil phase cannot be ignored in the heating process of the well patterns system for oilfield geothermal production.

@en

The chemical reaction in the reservoir causes fracture deformation during the heat extraction of enhanced geothermal systems (EGS), affecting thermal performance. The reaction rate is sensitive to temperature, concentration, and reaction-specific surface area. While previous research mainly focuses on the influence of temperature and concentration on fracture deformation, conversely, the effect of fracture morphology(aperture and tortuosity) is ignored. In this study, the deformation characteristics of rough and flat fractures are compared, and the influences of aperture and tortuosity on fracture deformation are analyzed. According to the influence law, the fracture deformation relationship equation between the aperture deformation rate with tortuosity and aperture is fitted. Results show that the deformation of rough fracture is significantly higher than that of flat fracture, and the variations of fracture aperture increase with the aperture and tortuosity. Furthermore, the influence of tortuosity (the variation of aperture increased by 47.97% when the tortuosity increased from 1 to 1.5) is greater than the aperture (that increased by 4.8% when the aperture increased from 0.2 to 1). The rate of aperture change is a logarithmic function of tortuosity and a power function of the aperture. These results provide significant references for the study of EGS, subsurface karst et al.

@en

The fractures are the main flow and heat transfer channel for fluids in deep high-temperature enhanced geothermal systems (EGS). The deformation of the fracture controlled by reactive flow is a common phenomenon during geothermal development, which might lead to a reduction in the system's thermal performance and operating life. While most previous research focuses on the influence of fracture deformation on system heat extraction performance, conversely the fracture deformation mechanism caused by the reactive flow is ignored. In this paper, a coupled thermal-hydraulic-chemical-deformation (THCD) model is established to investigate the fracture deformation mechanism. The deformation behavior of quartz and anorthite is compared, and the Damkohler number (Da) is adopted to explore the reaction mechanism. Meanwhile, the influence mechanism of concentration and temperature on fracture deformation is also analyzed. Results show that the Da at the fracture surface is 10⁻¹² -10⁻⁸, which means that the deformation of the fracture is controlled by the reaction rate. Compared with the concentration, the influence of temperature on fracture deformation is complex, and there is an inflection point. The inflection point is governed by the mineral reaction kinetics. These results provide significant references for the efficient heat extraction of the EGS.

@en

The multi-physics coupling process during the heat extraction from enhanced geothermal system, encompassing thermo(T)-hydro(H)-mechanical(M)-chemical(C) interactions, plays a pivotal role in changing geothermal reservoir characteristics. However, a comprehensive quantitative assessment of these multi-physics behaviors has been lacking. In this study, a novel approach was proposed to calculate the magnitude of mechanical, chemical, strong mechanical-chemical coupling, and weak mechanical-chemical coupling effects on the variations of reservoir characteristics. In particular, mechanical-chemical coupling effects are quantified for the first time. They are obtained by the fracture aperture difference results across five distinct coupling models (thermo-hydro, thermo-hydro-chemical, thermo-hydro-mechanical, partially-coupled four-field, and fully-coupled four-field models). The findings indicate that mechanical effects lead to an increase in fracture aperture, while chemical effects contribute to its reduction under underbalanced injection conditions. Strong mechanical-chemical coupling effects, exhibiting a negative correlation with chemical effects, conversely result in a diminished fracture aperture. The influences of these effects are investigated from the temporal and spatial perspectives. Temporally, mechanical effects dominate early production while chemical effects become prominent in later stages. Spatially, there mainly exists two zones when stable production: a mechanical-controlled region surrounding injection wells, and a chemical-controlled area distant from the injection wells. Furthermore, sensitivity analysis of injection concentration indicates its alternation changes the reservoir traits and production performance by modifying the magnitudes of chemical and mechanical-chemical coupling effects. This quantification of multi-physics effects offers insights into optimizing injection strategies for better geothermal development. The approach could hold promising potential in other geo-energy scenarios like carbon and hydrogen storage in reservoirs.

@en

Hot dry rocks (HDRs), as an essential renewable energy source, its development has received widespread attention, especially for heat extraction. The fracture is the main seepage and heat transfer channel of circulating fluid in dense HDR reservoirs, and its conductivity evolution significantly affects the production performance. Most existing studies have focused on the change of fracture conductivity under elastic deformation without considering the additional conductivity induced by rock damage. However, the additional conductivity may have significant implications for rational design and timely adjustment of the production scheme. Therefore, a three-dimensional model at the field-scale is established, and it is used to analyze the effect of additional conductivity on production performance and economic efficiency. To simplify the calculation, the actual forms of damage are equivalent to the macroscopic physical evolution of the matrix. Results show that the rock is mainly tensile failure affected by thermal stress during production. The occurrence of damage will increase the reservoir permeability and porosity, reduce Young's modulus, and then reduce the differential pressure and production temperature, with a maximum reduction of 2.21 MPa and 14.21 °C in the control case, respectively. The effects of injection temperature, Young's modulus, and injection mass flow on the production performance are significant, followed by Poisson's ratio. In contrast, production pressure and fracture initial permeability had less influence. The maximum differential economic benefit of the control case is up to 2.289 million RMB. This research proves the necessity of damage study during the long-term production of HDRs.

@en

The efficient development of hydrothermal resource is an important focus in the global geothermal industry. A novel hydrothermal open-loop geothermal system in a multi-lateral horizontal well is presented to achieve the simultaneous exploitation of multi-layer reservoirs and accelerate the large-scale development of the hydrothermal resource. To investigate the detailed performance of the novel geothermal system, a novel 3D fluid flow and heat transfer numerical model is established based on the local thermal non-equilibrium and dual porosity theories. The flow and temperature fields are analyzed, while the pressure and temperature interferences between the lateral wellbores of the multi-lateral horizontal well system are comprehensively investigated. The influences of the injection flow rate, lateral-wellbore number, spacing between the injection segment and the main wellbore, fracture aperture on the heat extraction are studied. Furthermore, the thermal characteristics of the novel geothermal system and traditional doublet well system are compared for the first time. The investigation indicates that the flow impedance of the novel geothermal system can decrease by 75.38%, and its energy efficiency reaches 8 times that of the doublet well system. The findings can provide a key reference for scientific research and application of the novel and traditional geothermal systems.

@en

There are chemical dissolution and precipitation in geothermal formations with the long-term exploitation of hot dry rocks, which would induce the remarkable alteration of fracture transmissibility and evolution of thermal performance. The objective of this study is to identify the influence of various parameters incorporating chemical reactions to provide suggestions for the operation of the enhanced geothermal system (EGS). A 3D thermal-hydraulic-chemical coupling model is established for a fractured EGS considering silica-water reaction kinetics at fracture surface. Parametric study and sensitivity analysis are innovatively conducted to investigate the effects of key operation and formation parameters, including injection concentration (c_in), injection temperature (T_in), injection rate (Q_in) and initial formation silica reactant content (C_i), on thermal extraction performance. Results indicate that undersaturated injection causes chemical dissolution while oversaturated injection leads to silica precipitation. Higher c_in relieves the thermal breakout (10 °C) but enhances the injection-production pressure difference (13.5 MPa). The selection of c_in is an efficient way to adjust geothermal production, especially for oversaturated injection. The influence degrees of c_in on production temperature (T_out) and pressure difference (Δp) under oversaturated injection are 8 times and 2 times larger than those under undersaturated injection. T_in (<80 °C) exerts little influence on the chemical reaction. The remarkable water storage phenomenon will be induced with the increase of Q_in when precipitation occurs. Higher C_i results in lower T_out and Δp. The sensitivity analysis provides the rank of major parameters to thermal performance: Q_in, T_in, c_in and C_i. The standard sensitivities of Q_in to T_out, Δp and net thermal power N are 0.78, 0.63, and 0.88 respectively. As for the contribution ratio, c_in contributes the most (89%) to the variation of fracture aperture while T_in (<80 °C) contributes the least (nearly 0%). The research will provide a significant reference for the selection of parameters during the development of geothermal energy.

@en

Geothermal energy is gaining attention as an environmentally friendly and sustainable alternative to fossil fuels. There exist complicated coupling processes among fluid flow, heat transfer, mechanical deformation, and chemical reaction in fractured geothermal reservoirs. Understanding of thermo-hydro-mechanical-chemical (THMC) coupled process is of significance for the efficient development of enhanced geothermal systems. This research aims to quantify the contributions of mechanical and chemical behaviors to reservoir feature variation during geothermal production. Herein, a fully coupled THMC model is developed for the production of an enhanced geothermal system based on fracture aperture alternation. The key improvement lies in the inclusion of mechanical-chemical (MC) coupling where chemical behaviors cause stress to be redistributed. On a geographical and temporal scale, the unique proportion technique is employed to study ratios of mechanical to chemical effects and poroelasticity to thermoelasticity. The findings emphasize the necessity of incorporating the MC coupling when analyzing reservoir characteristic variations. The mechanical effect plays a major role near the injection well and in a short term. The chemical effect is dominated far from the injection well and its time scale of action is long-term. The intensity of thermoelasticity is stronger than that of poroelasticity. Moreover, the proportion analysis reveals that when stable production is achieved, there are primarily thermo-driven and chemical-driven regions. To improve production performance and adjust geothermal reservoir features, injection temperature and solute concentration should be optimized. This paper serves as a valuable reference for both the theoretical multi-physical coupling model and practical geothermal production.

@en

Global warming induced by excessive carbon dioxide emission generated from fossil fuels consumption is becoming increasingly severe. Development and utilization of clean and renewable energy is a significant measure to contribute for the carbon neutrality. Geothermal resource is an important clean and renewable energy, because of its advantages of abundant reserves, exploitation without restrictions from the climate and season. Extracting heat from the high temperature geothermal reservoir requires well drilling, completion and fracturing. However, due to some particular conditions in geothermal development, such as high temperature, erosion, leakage, and hard rocks, traditional petroleum well drilling and completion technologies are not well applicable in the geothermal industry. In this paper, the state-of-the-art geothermal drilling, and completion technologies as well as several fracturing methods, such as foam drilling, air hammer drilling, hydrothermal jet drilling, expansion casing technology, high temperature resistant drilling fluid and cement, chemical stimulation and thermal stimulation are comprehensively reviewed. Besides, heat extraction performances of water and CO₂ as working fluids in the geothermal system are discussed. A few novel heat extraction methods including closed loop, open loop, enhanced geothermal systems are compared and analyzed. This review is expected to provide suggestions for research fields worthy of study for geothermal drilling and exploitation.

@en

Horizontal ground heat exchanger (HGHE) is a low-cost shallow geothermal exploitation method compared with vertical ones, but it needs more land area. Optimizing geometrical configurations and arrangements of ground heat exchangers is desirable to enhance heat exchange per unit land area. However, a comprehensive study on different HGHE configurations and arrangements considering atmosphere-soil-HGHE interaction is still inadequate. To this end, thermal performances of four configurations, including slinky, spiral, multi-tube with serial connection and tiled tubes, are compared based on a 3D numerical model for HGHE considering atmosphere conditions variation and humidity migration. Effects of main geometrical parameters on HGHE performance and economical efficiency are investigated for these four configurations. Finally, different tube arrangements of multi-tube HGHE are designed and their thermal performances are analyzed. Results indicated that longer HGHE tube length has better thermal performance but lower economical efficiency, thus being not advisable for the HGHE design. Among these four configurations, multi-tube with serial connection shows the best thermal performance. Increasing the tube spacing is beneficial for reducing the tube mutual interference and enhance the heat exchange per unit land area, which is recommendable. Results are expected to provide meaningful suggestions for the HGHE design.

@en

Geothermal reservoir scaling controlled by reactive flow is a common phenomenon during geothermal development, which might lead to the decline of geothermal productivity and abandonment of geothermal projects. It is significant to reveal chemical influence on the geothermal reservoir characteristics and production performance to provide the reference for geothermal sustainable development. In this work, a thermal-hydraulic-chemical coupling model is developed for a triple-well fractured enhanced geothermal system considering the silica dissolution and precipitation at the fracture surface. It is focused on the evolution of fracture aperture and geothermal productivity. The results indicate that undersaturated injection will cause mineral dissolution which leads to the rise of fracture transmissibility and a thermal break. The fracture width is enhanced by 0.093 mm, and the permeability is 8 times the initial value after 20-year production. The oversaturated injection will cause mineral precipitation which leads to the reduction of fracture aperture and the increase of pressure difference. The opening is reduced by 0.041 mm and the permeability goes down by 2 orders of magnitude. The maximum increases of average production temperature and differential pressure are 10 °C and 13.5 MPa respectively in the 20th year from undersaturation to supersaturation. Besides, it is figured out that reaction rate denotes the evolution of fracture aperture. Both of them are distributed as a band at a fracture plane, which is determined by saturation index and temperature together. The former controls the reaction pattern (dissolution/precipitation) while the latter controls the reaction speed. Furthermore, the evolution of opening is uneven around the production well at different exploitation stages. Fracture width is higher on the strong side in the first 8.7 years but is higher on the weak side in later development when undersaturated water (0 mol/kg) is injected. The phenomenon can be employed to judge the production status, guide the descaling, etc. Production analysis indicates that oversaturation delays the heat breakthrough and undersaturation reduces the surface pump input. This study will provide a significant reference for the geothermal sustainable production influenced by chemical reactions in an enhanced geothermal system.

@en

As a kind of clean renewable energy, the production and utilization of geothermal resources can make a great contribution to optimizing the energy structure and energy conservation and emission reduction. The circulating heat extraction process of working fluid will disturb the equilibrium state of physical and chemical fields inside the reservoir, and involve the mutual coupling of heat transfer, flow, stress, and chemical reaction. Revealing the coupling mechanism of flow and heat transfer inside the reservoir during geothermal exploitation can provide important theoretical support for the efficient exploitation of geothermal resources. This paper reviews the research advances of the multi-field coupling model in the reservoir during geothermal production over the past 40 years. The thrust of this paper is on objective analysis and evaluation of the importance of each coupling process and its influence on reservoir heat extraction performance. Finally, we discuss the existing challenges and perspectives to promote the future development of the geothermal reservoir multi-field coupling model. An accurate understanding of the multi-field coupling mechanism, an efficient cross-scale modeling method, as well as the accurate characterization of reservoir fracture morphology, are crucial for the multi-field coupling model of geothermal production.@en

Salt cavern is one of the significant types of underground gas storage. During the construction of the salt cavern gas storage, part of the insoluble particles will remain in the cavity. The insoluble particles gradually deposit in the bottom of the salt cave, which reduces the cavity construction speed and lessens the volume of gas storage. The high-pressure rotating water jet technology is simulated to flush insoluble particles, and the flushing system is mainly composed of surface facilities and downhole flushing tools. The seven high-pressure rotating water jets can flush the insoluble particles at the bottom of the cavity. A mathematical model with the solid-liquid two-phase flow is established and numerical simulated first. Then, this paper analyzes the influence of the injection flowrate, the distance from inlet to outlet, the rotational velocity on the flushing efficiency. The results show that the injection flowrate and the distance from inlet to outlet are the two main controlling factors. Under the conditions of this paper, the injection flowrate is 100 m³/h, the rotational velocity is 300 r/min~400 r/min, the distance from inlet to outlet is between 0.5 m~1.0 m, and the flushing efficiency is higher than that of other range studied. The research is helpful to understand the flushing rule and optimize the flushing parameters during the construction of the gas storage through water jetting technology.

@en

Operational parameters optimization is of great significance to improve overall heat extraction performance from the hydrothermal or enhanced geothermal system. Injection flowrate, injection temperature, and production pressure are several of the easily human-controlled parameters to make the best of limited geothermal resources during the planned life. The net heat power and flow impedance are two contradictory production indexes to be optimized for efficient exploitation of long-term geothermal production. In this study, an integrated approach of finite element, multiple regression, non-dominated sorting genetic algorithm, and the technique for order preference by similarity to ideal solution is proposed and applied to the multilateral-well system to realize the optimization of geothermal development. Firstly, parametric cases coupling with thermal and hydraulic models are analyzed. Then, multiple regression is employed to obtain the net heat power and flow impedance functions considering human-controlled operational parameters and reservoir physical properties. Afterward, the multi-objective optimization algorithm is used to gain the Pareto solution set of injection and production parameters. Finally, the technique for order preference by similarity to ideal solution is employed to select the optimal combination of operational parameters for geothermal extraction. It is concluded that water loss and thermal drawdown are necessarily considered in the optimization process. The proposed approach represents global optimization. Operational parameters for the optimal case are (Q_in, T_in, p_out) which equal (49.98 °C, 62.21 kg/s, 27.44 MPa) under the conditions of this study. From the comparison between the base case and optimal case, it is observed that horizontal spread length, 2D swept area and 3D swept volume of the low-temperature scope is reduced by 77.9 m, 10 × 10⁴ m², and 16 × 10⁶ m³. Besides, the thermal breakthrough has been delayed by 1.1 years. The water loss decreases by 36.23%. The optimal case demonstrates a great improvement in production performance. Sustainable exploitation is achieved through multi-objective optimization for operational parameters. The proposed method reflects superiority, efficiency, and intelligence in geothermal development.

@en

The effective exploitation of deep hydrothermal resource is an emphasis in the global geothermal industry. A horizontal-well coaxial closed-loop geothermal system (CCGS) is presented to realize its efficient development by significantly increasing the heat exchange area based on a long horizontal well. However, few researchers focus on the complex convection process in the hydrothermal reservoir, making it difficult to effectively evaluate the thermal characteristics of a geothermal system. Therefore, a 3D transient flow and heat transfer model considering the forced and natural convection in the reservoir is established to achieve a comprehensive CCGS analyzes. Then, the flow and temperature fields in the formation with forced and natural convection are compared. The thermal characteristics of the horizontal-well CCGS with and without convection are discussed. Also, the influences of key factors on the heat production of the horizontal-well CCGS are investigated. The differences between the horizontal-well and vertical-well CCGSs in exploiting the hydrothermal resource are analyzed. The results show that the fluid convection in the reservoir, particularly the forced convection, can greatly enhance heat transfer by affecting heat convection. The temperature field in the reservoir may return to the original state before the next heating season when a high convection velocity occurs, which is beneficial to realize 100% sustainable heat production. The horizontal length, permeability, and reverse circulation should be the primary considerations in improving heat extraction. Compared with the limited performance of vertical-well CCGS, the horizontal-well CCGS can fully extract the heat stored in the hydrothermal reservoir via a long horizontal well. The findings provide a vital reference for the CCGS scientific study and engineering application.

@en

Operational parameter optimization is of great significance to improve overall heat extraction performance from hydrothermal or enhanced geothermal systems. Injection flowrate/temperature and production pressure are relatively easy to control to optimize the exploitation of geothermal resources during the planned reservoir lifetime. The net heat power and flow impedance are two contradictory production indexes for describing the exploitation effect. One indicates energy efficiency from reservoirs, the other represents the mining difficulty or artificial energy input. In this study, a multi-objective optimization procedure is proposed and applied to a synthetic multilateral-well system for 30 years. Firstly, a multilateral-well geothermal model coupled with thermal and hydraulic parameters is established. Then, a multiple regression method is employed to obtain the net heat power and flow impedance functions with injection/production parameters and physical properties of the reservoir. Finally, a multi-objective genetic algorithm is used to gain a Pareto solution set of injection and production parameters. A comparison with the base case indicates the superiority, high efficiency, and intelligence of multi-objective optimization.

@en