This paper introduces a comprehensive thermo-hydro-mechanical (THM) modeling framework tailored for high-temperature aquifer thermal energy storage (HT-ATES) systems. Our framework presents a novel dual-assessment approach that simultaneously evaluates thermal performance and geomechanical stability of HT-ATES systems. The framework combines advanced sensitivity analysis with multi-objective optimization to concurrently boost thermal efficiency and maintain geomechanical safety. The model simulates the cyclic injection-extraction process while capturing the interdependent effects of heat transfer, fluid flow, and mechanical stress evolution. A distance-based Generalized Sensitivity Analysis (DGSA) is applied to identify and rank the most critical parameters influencing system performance and stability, particularly in regions such as the cold well and overlying caprock. Furthermore, surrogate models constructed with eXtreme Gradient Boosting (XGBoost) facilitate a computationally efficient Non-dominated Sorting Genetic Algorithm II (NSGA-II) optimization that investigates the trade-offs between enhancing heat production and minimizing failure risks. Validation against high-fidelity simulations reveals that, compared to a benchmark model with a thermal recovery efficiency of approximately 85% and a caprock slip tendency of 34°, the optimized designs achieve around 88% efficiency and reduce the caprock slip tendency to 29°. These quantitative improvements demonstrate that the proposed framework significantly enhances both energy production and geomechanical stability, offering valuable guidance for the design of robust HT-ATES systems under fixed geological conditions.

Geothermal energy extraction through deep mine systems offers the potential to reduce the cost of geothermal systems while meeting the cooling needs of deep mines. However, the injection of cold water into the subsurface triggers strongly coupled thermo-hydro-mechanical (THM) processes that can affect the stability of underground excavations. This study evaluates the impact of geothermal energy extraction on the temperature and stability of a deep mine. By quantifying the sensitivity of the mine temperature and stability to various parameters, we propose a scheme to optimize geothermal energy production, while achieving rapid mine cooling and maintaining stability. We first evaluate the impact of geothermal operations on mine temperature and stability through THM numerical modeling. The simulations show that poro-elastic stress quickly affects mine stability, while thermal stress has a more significant impact on the long-term stability. We then use Distance-based Generalized Sensitivity Analysis (DGSA) to quantify parameter sensitivity. The analysis identifies the distance between the mine system and the geothermal system as the most influential factor. Other important parameters include the injection rate, injection temperature, well spacing, coefficient of thermal expansion, permeability, Young's modulus, and heat capacity. Finally, we propose a DGSA-based optimization framework that accounts for subsurface uncertainty and validate the optimized results. Our results indicate that, with favorable geological conditions, a rational selection of system design parameters can enhance geothermal energy production while ensuring rapid mine cooling and stability. This study provides essential insights for the optimization of deep mine geothermal systems and supports effective decision-making.

Hydraulic fracturing of coalbed methane wells has been widely practised as an effective method to increase drainage efficiency in low-permeability, low-pressure and low-saturated coal seams. To investigate hydraulic fracture performance and associated seismic response in coal, hydraulic fracturing experiments were carried out on two cubic coal blocks containing a host of natural fractures using a true triaxial rock testing machine equipped with loading, injection and acoustic systems. The acoustic system uses transducers with active sources to repetitively generate and receive ultrasonic P/S wave pulses for characterising mechanical properties of the coal blocks and revealing fracture growth. Silicon oil was injected into the middle of coal blocks to create hydraulic fractures under deviatoric stress conditions, and the stress and displacement, borehole pressure and volume, and seismic response were recorded over the injection process. X-ray computed tomography (CT) was conducted before and after the experiments to identify the location and geometry of hydraulic and natural fractures. Results have shown that the fracturing behaviour, the drawdown period of borehole pressure and the intrusion of fracturing fluid are dominated by the complexity and insulation offered by internal natural fracture networks of coal blocks. In addition, seismic spectrograms captured both fracture initiation and its subsequent interaction with natural fractures, which indicates that the induced fracture and fracturing fluid interfere with the propagation of seismic waves and influence ultrasonic seismic characteristics. Seismic velocity tomography of ultrasonic acoustic signals recorded also provided the spatial information of fractures, such as approximate locations of pre-existing fractures and injection-disturbed regions.

Two 0.3 m × 0.3 m × 0.3 m shale blocks, one representing a homogeneous sample while the other representing a naturally fractured sample, are modelled using the lattice based DEM code, XSite. The synthetic rock mass approach (SRM), which assigns the smooth joint contacts (SJM) to the weakness planes, is used to represent the natural fractures in shale block-2. Firstly, the developed models are compared with the findings of previously conducted true-triaxial hydraulic fracturing experiments with acoustic measurements, and their subsequent computed tomography (CT) and seismic velocity tomography results. The 3D model results confirmed the curved shape hydraulic fractures, which propagated perpendicular to the minimum stress directions in both shale blocks. Model results also captured the natural fracture (NF) and hydraulic fracture (HF) interaction, particularly the arrest, the dilation of major NFs, followed by crossing with offset mechanism, in shale block-2. Secondly, the parametric studies are carried out to investigate the role of fluid flow rate (q), and fluid viscosity (µ) on different NF/HF interaction mechanisms. The effects of q and µ are discussed based on the total stimulated area including the tensile and shear microcracks, the pipe apertures, and the pressure evolutions within NFs.