Many underground reservoirs for storing water have been constructed in China’s western coal mines to protect water resources. Coal pillars which work as dams are subjected to a long-term soaking environment of concentrated salty water. Deterioration of the coal dam under the attack of the salty solution poses challenges for the long-term stability and serviceability of underground reservoirs. The evolution of the physical and mechanical properties of coal subjected to salty solutions are investigated in this paper. Coal from a western China mine is made to standard cylinder samples. The salty solution is prepared according to chemical tests of water in the mine. The coal samples soaked in the salty solution for different periods are tested by scanning electron microscope, nuclear magnetic resonance, and ultrasonic detector techniques. Further, uniaxial compression tests are carried out on the coal specimens. The evolutions of porosity, mass, microstructures of coal, solution pH values, and stress–strain curves are obtained for different soaking times. Moreover, a damage constitutive model for the coal samples is developed by introducing a chemical-stress coupling damage variable. The result shows that the corrosion effect of salty solution on coal samples becomes stronger with increasing immersion time. The degree of deterioration of the longitudinal wave velocity (v_p) is positively correlated with the immersion time. With the increase in soaking times, the porosity of coal gradually increases. The relative mass firstly displays an increasing trend and then decreases with time. The peak strength and elastic modulus of coal decreases exponentially with soaking times. The developed damage constitutive model can well describe the stress–strain behavior of coal subjected to salty solution under the uniaxial compression.

The instability and failure of coal pillars is one of the important factors leading to the catastrophic consequences of coal mine goaf collapse. Coal mine water has the characteristics of high salinity. Long-term mine water erosion can easily deform the coal pillar structure, eventually leading to instability and damage. This study carried out tests on coal samples soaked in salt solutions with different concentrations, and the nuclear magnetic resonance (NMR) method was used to obtain the dynamic evolution of the pore-fracture structure of coal. On the basis of fractal theory, the changes in fractal dimension of pore structure during the soaking process were discussed. The damage variable based on the pore fractal dimension was defined and the evolution relationship between the damage variable and immersion time was characterized. The findings demonstrated that the porosity change rate has an exponentially increasing relationship with the immersion time, and with the increasement of concentration of salt solution, the porosity change rate also shows increasing trends. The number of seepage pores and total pores increased with the immersion time. While, with the extension of soaking time, the number of adsorption pores first increased and then decreased. The connectivity between pores was enhanced. The relationship between the fractal dimension and the immersion time is linearly decreasing. The damage variable showed an increasing trend with the immersion time. As the concentration of salt solution increased, the damage of coal increased. The research results are of great significance for rationally evaluating the stability of coal pillars and ensuring the safe operation of underground engineering.

Numerical Research of Fluid Flow and Solute Transport in Rough Fractures under Different Normal Stress: The effects of roughness and normal stress on hydraulic properties of fractures are significant during the coupled shear flow test. Knowing the laws of fluid flow and solute transport in fractures is essential to ensure the nature and safety of geological projects. Although many experiments and numerical simulations of coupled shear flow test have been conducted, there is still a lack of research on using the full Navier-Stokes (N-S) equation to solve the real flow characteristics of fluid in three-dimensional rough fractures. The main purpose of this paper is to study the influence of roughness and normal stress on the fluid flow and solute transport through fractures under the constant normal stiffness boundary condition. Based on the corrected successive random addition (SRA) algorithm, fracture surfaces with different roughness expressed by the Hurst coefficient () were generated. By applying a shear displacement of 5 mm, the sheared fracture models with normal stresses of 1 MPa, 3 MPa, and 5 MPa were obtained, respectively. The hydraulic characteristics of three-dimensional fractures were analyzed by solving the full N-S equation. The particle tracking method was employed to obtain the breakthrough curves based on the calculated flow field. The numerical method was verified with experimental results. It has been found that, for the same normal stress, the smaller the fracture value is (i.e., more tough the fracture is), the larger the mechanical aperture is. The ratio of hydraulic aperture to mechanical aperture ( ) decreases with the increasing of normal stress. The smaller the value, the effect of the normal stress on the ratio is more significant. The variation of transmissivity of fractures with the flow rate exhibits similar manner with that of . With the normal stress and value increasing, the mean velocity of particles becomes higher and more particles move to the outlet boundary. The dispersive transport behavior becomes obvious when normal stress is larger.

Uniaxial compression tests were conducted on coal samples subjected to different dry–wet cycling treatments to investigate the damage and degradation mechanisms of coal samples under the dry–wet cyclic action of acidic, high-salinity solutions. The damage process of the coal samples was monitored in situ using acoustic emission (AE). The degradation evolution of the mechanical parameters and macroscopic failure modes with the number of cycles was analyzed. Based on the AE ringing parameters, the RA-AF distribution and the AE fractal dimension’s variation characteristics were studied. Additionally, scanning electron microscopy (SEM) was used to observe the microstructure of the coal samples. The results showed that with the increase in the number of dry–wet cycles, both the peak strength and elastic modulus of the coal samples exhibited varying degrees of degradation, and the failure mode gradually shifted from tensile failure to shear failure. AE ringing counts decreased progressively, while the proportion of shear cracks based on the RA-AF classification increased. At the same time, the mean AE fractal dimension of the coal samples increased, and the fractal dimension decreased with an increase in AE ringing counts. The sharp drop in fractal dimensions could serve as an early warning signal for a major failure in the coal samples. Furthermore, under the influence of dry–wet cycling in acidic, high-salinity solutions, defects such as pores and cracks in the microstructure of the coal samples became more pronounced, and the degradation effect continuously intensified.

The hydraulic properties of fractures are greatly affected by the stress. Knowing the fluid flow behavior of fractures is of great importance to underground engineering construction and environmental safety. The main purpose of this paper is to study the fluid flow characteristics of rough fractures under different stress states. First, rough fracture surfaces were generated by using the corrected successive random addition (SRA) algorithm. Then, the sheared fracture models subjected to different stress condition were obtained under the boundary condition of constant normal stiffness (CNS). Finally, the hydraulic characteristics of the three-dimensional rough rock fractures were analyzed by numerically solving the full Navier–Stokes equation. It has been found that (1) the aperture of fractures all obeys the Gaussian distribution. The dilatancy effect is gradually obvious and aperture becomes larger with the increase of shear displacement. (2) When the initial normal stress increases, the contact area of fracture becomes larger and the reverse flow can be observed around the contact area. (3) The relationship between hydraulic gradient and flowrate exhibits nonlinearity which can be described by the Forchheimer’s law. The linear coefficient a and the nonlinear coefficient b gradually decrease with the increase of shear displacement and finally stabilize. The values of a and b are reduced by 1–2 and 1–3 orders of magnitude respectively during the shear. The critical Reynolds number increases with the increase of shear displacement and decrease as the initial normal stress increases.

Closed-loop geothermal systems (CLGSs) avoid groundwater production and offer stable deep heat supply, but their long-term performance hinges on reliable coupling between the wellbore, the near-well interface and the surrounding formation. Using the D22 well in the Xiongan New Area (deep carbonate reservoir), we built a three-domain thermo-hydraulic framework that updates CO₂ properties with temperature and pressure and explicitly accounts for wellbore-formation thermal resistance. Two geometries (U-tube and single-well coaxial) and two working fluids (CO₂ and water) were compared and optimized under field constraints. With the coaxial configuration, CO₂ delivers an average thermal power of 186.3 kW, exceeding that of water by 44.9%, while the fraction of wellbore heat loss drops by 3–5%. Under field-matched conditions, the predicted outlet temperature (76.8 °C) agrees with the measured value (77.2 °C) within 0.52%, confirming the value of field calibration for parameter transferability. Long-term simulations indicate that after 30 years of continuous operation the outlet temperature decline remains <8 °C for CO₂, outperforming water and implying better reservoir utilization and supply stability. Sensitivity and Pareto analyses identify a practical operating window, i.e., flow velocity of 0.9–1.1 m s⁻¹ and depth of 3000–3500 m, favoring the single-well coaxial + CO₂ scheme. These results show how field-calibrated modeling narrows uncertainty and yields implementable guidance on geometry, operating conditions, and wellbore insulation strategy. This study provides quantitative evidence that CO₂-CLGSs in deep carbonate formations can simultaneously increase thermal output and limit long-term decline, supporting near-term engineering deployment.