Pub Date : 2025-11-10DOI: 10.1016/j.ijrmms.2025.106327
Zhenting Sun , Lei Ma , Quan Li , Yaping Deng , Han Qiu , Haichun Ma , Cihai Chen , Yongshuai Yan , Jiazhong Qian
Characterizing fracture networks is critical for groundwater development, geothermal exploitation, hydrocarbon recovery, and geological CO2 sequestration, yet their complex and uncertain spatial distribution poses a persistent challenge. This study proposes an intelligent inversion framework that integrates a 3D-UNet surrogate model, reversible jump Markov Chain Monte Carlo (rjMCMC), and multi-source data fusion for three-dimensional discrete fracture network (DFN) characterization at the field scale. Within this framework, a 3D-UNet model trained on large datasets of fracture configurations, hydraulic head, and electrical potential provides an efficient initial inversion of fracture parameters. Fracture geometries are then extracted with the RANSAC algorithm and iteratively refined via rjMCMC, where the surrogate 3D-UNet replaces conventional forward modeling. This innovation reduces computational costs by an order of magnitude, enabling efficient large-scale inversion. Furthermore, the fusion of electrical potential with hydraulic head data enhances inversion accuracy by about 10 %. Validation demonstrates that the framework reliably reconstructs the spatial distribution of fracture networks, capturing both low-density zones and dominant hydraulic pathways in highly heterogeneous domains. By combining computational efficiency with improved accuracy, this approach offers a practical and scalable solution for field-scale fracture network characterization in a wide range of hydrogeological and engineering applications.
{"title":"Three-dimensional discrete fracture network identification based on deep learning and reversible jump Markov chain Monte Carlo algorithm","authors":"Zhenting Sun , Lei Ma , Quan Li , Yaping Deng , Han Qiu , Haichun Ma , Cihai Chen , Yongshuai Yan , Jiazhong Qian","doi":"10.1016/j.ijrmms.2025.106327","DOIUrl":"10.1016/j.ijrmms.2025.106327","url":null,"abstract":"<div><div>Characterizing fracture networks is critical for groundwater development, geothermal exploitation, hydrocarbon recovery, and geological CO<sub>2</sub> sequestration, yet their complex and uncertain spatial distribution poses a persistent challenge. This study proposes an intelligent inversion framework that integrates a 3D-UNet surrogate model, reversible jump Markov Chain Monte Carlo (rjMCMC), and multi-source data fusion for three-dimensional discrete fracture network (DFN) characterization at the field scale. Within this framework, a 3D-UNet model trained on large datasets of fracture configurations, hydraulic head, and electrical potential provides an efficient initial inversion of fracture parameters. Fracture geometries are then extracted with the RANSAC algorithm and iteratively refined via rjMCMC, where the surrogate 3D-UNet replaces conventional forward modeling. This innovation reduces computational costs by an order of magnitude, enabling efficient large-scale inversion. Furthermore, the fusion of electrical potential with hydraulic head data enhances inversion accuracy by about 10 %. Validation demonstrates that the framework reliably reconstructs the spatial distribution of fracture networks, capturing both low-density zones and dominant hydraulic pathways in highly heterogeneous domains. By combining computational efficiency with improved accuracy, this approach offers a practical and scalable solution for field-scale fracture network characterization in a wide range of hydrogeological and engineering applications.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"197 ","pages":"Article 106327"},"PeriodicalIF":7.5,"publicationDate":"2025-11-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145476166","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-09DOI: 10.1016/j.ijrmms.2025.106338
Kai Guan , Runze Zhu , Ignacio Pérez-Rey , Wancheng Zhu , Xige Liu , Jianyu Zhou
Rock creep and dynamic behaviors are distinct mechanical responses with different strain rates, and their interaction can cause asynchronous deformation in anchoring systems, reducing bolt prestress and increasing time-dependent instability risk. Using the self-developed rock creep-impact testing machine, this study highlights that under combined creep and dynamic loading, unbolted specimens tend to experience delayed failure, whereas bolt-reinforced specimens fail more promptly during impact, indicating improved predictability and stability due to reinforcement. The application of bolt prestress significantly enhances impact resistance by suppressing axial strain increases and damage during dynamic events, thereby extending the time-to-failure and improving overall performance. During creep, bolt strain gradually increases, but impact causes rapid escalation, demonstrating that transient disturbances are more effective in activating bolt reinforcement than slow creep. Repeated dynamic impacts diminish anchoring effectiveness, increasing acoustic emission energy, but higher prestress levels delay weakening and facilitate a transition to more controlled energy dissipation. Prestress initially decline rapidly before stabilizing, with subsequent impacts inducing stepwise reductions and occasional abnormal rebounds that may serve as early-warning signals for potential failure. Prestress relaxation arises from bolt elongation pre-impact and time-dependent damage to the rock mass post-impact, necessitating timely re-tensioning in vibration-prone environments. The progression and failure of cracks are significantly influenced by prestress levels, with higher prestress shifting through-cracking extending along the joint towards both the top and bottom to propagating laterally across the specimen, especially near the tray region, thereby reducing localized damage. Overall, the findings underscore the critical role of prestress management and reinforcement strategies in improving the resilience of anchoring systems under creep stress and dynamic impact conditions, contributing to safer and more durable rock engineering applications.
{"title":"Creep behavior and prestress relaxation mechanism of bolt-reinforced jointed specimen disturbed by dynamic impact","authors":"Kai Guan , Runze Zhu , Ignacio Pérez-Rey , Wancheng Zhu , Xige Liu , Jianyu Zhou","doi":"10.1016/j.ijrmms.2025.106338","DOIUrl":"10.1016/j.ijrmms.2025.106338","url":null,"abstract":"<div><div>Rock creep and dynamic behaviors are distinct mechanical responses with different strain rates, and their interaction can cause asynchronous deformation in anchoring systems, reducing bolt prestress and increasing time-dependent instability risk. Using the self-developed rock creep-impact testing machine, this study highlights that under combined creep and dynamic loading, unbolted specimens tend to experience delayed failure, whereas bolt-reinforced specimens fail more promptly during impact, indicating improved predictability and stability due to reinforcement. The application of bolt prestress significantly enhances impact resistance by suppressing axial strain increases and damage during dynamic events, thereby extending the time-to-failure and improving overall performance. During creep, bolt strain gradually increases, but impact causes rapid escalation, demonstrating that transient disturbances are more effective in activating bolt reinforcement than slow creep. Repeated dynamic impacts diminish anchoring effectiveness, increasing acoustic emission energy, but higher prestress levels delay weakening and facilitate a transition to more controlled energy dissipation. Prestress initially decline rapidly before stabilizing, with subsequent impacts inducing stepwise reductions and occasional abnormal rebounds that may serve as early-warning signals for potential failure. Prestress relaxation arises from bolt elongation pre-impact and time-dependent damage to the rock mass post-impact, necessitating timely re-tensioning in vibration-prone environments. The progression and failure of cracks are significantly influenced by prestress levels, with higher prestress shifting through-cracking extending along the joint towards both the top and bottom to propagating laterally across the specimen, especially near the tray region, thereby reducing localized damage. Overall, the findings underscore the critical role of prestress management and reinforcement strategies in improving the resilience of anchoring systems under creep stress and dynamic impact conditions, contributing to safer and more durable rock engineering applications.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106338"},"PeriodicalIF":7.5,"publicationDate":"2025-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145473193","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-08DOI: 10.1016/j.ijrmms.2025.106328
Licheng Sun , Lixin Wu , Youyou Xu , Guangrui Dong , Tao Zheng , Wenfei Mao
The spatial differences of pressure-stimulated rock potential (PSRP) and its correlation with crack development and distribution were poorly studied, despite their importance in understanding seismic anomalies and monitoring engineering rock mass. In this study, compressive stresses were partly applied to specially designed diorite specimens in shape of composite cuboid-platform, which make the specimens being comprised of loaded cuboid and free platform. The PSRP differences between each end of the free platform and the loaded cuboid were monitored separately, and acoustic emission was monitored simultaneously to observe crack development. Results revealed distinct stage-specific PSRP characteristics: exponential drop during initial loading, rapid rise with local macro failure, pulse appear with instability failure, and large fluctuation rise occur with post-peak loading. Preceding the instability failure phase, the PSRPs variation across different regions exhibited consistent trend with differences in amplitude and timing. At instability failure phase, however, significant spatial differences emerged, with opposite polarity pulses observed. The detected PSRP could be attributed to the fluid electrokinetic effect and the activation of positive holes (P-holes). The PSRP corresponds to crack development, as larger accumulated crack areas activate more P-holes and consequently generate higher PSRP rises. Spatial distribution in macro fracture surfaces further governed the spatial differences of PSRP, where open fractures blocked the upward transmission of P-holes, affecting the response of P-holes in PSRP. The revealed relationship between PSRP and crack development/distribution provides novel insights for monitoring and early warning of engineering rock mass instability and earthquakes.
{"title":"Characteristics and spatial differences of pressure-stimulated rock potential in diorite partly loaded to fracturing","authors":"Licheng Sun , Lixin Wu , Youyou Xu , Guangrui Dong , Tao Zheng , Wenfei Mao","doi":"10.1016/j.ijrmms.2025.106328","DOIUrl":"10.1016/j.ijrmms.2025.106328","url":null,"abstract":"<div><div>The spatial differences of pressure-stimulated rock potential (PSRP) and its correlation with crack development and distribution were poorly studied, despite their importance in understanding seismic anomalies and monitoring engineering rock mass. In this study, compressive stresses were partly applied to specially designed diorite specimens in shape of composite cuboid-platform, which make the specimens being comprised of loaded cuboid and free platform. The PSRP differences between each end of the free platform and the loaded cuboid were monitored separately, and acoustic emission was monitored simultaneously to observe crack development. Results revealed distinct stage-specific PSRP characteristics: exponential drop during initial loading, rapid rise with local macro failure, pulse appear with instability failure, and large fluctuation rise occur with post-peak loading. Preceding the instability failure phase, the PSRPs variation across different regions exhibited consistent trend with differences in amplitude and timing. At instability failure phase, however, significant spatial differences emerged, with opposite polarity pulses observed. The detected PSRP could be attributed to the fluid electrokinetic effect and the activation of positive holes (P-holes). The PSRP corresponds to crack development, as larger accumulated crack areas activate more P-holes and consequently generate higher PSRP rises. Spatial distribution in macro fracture surfaces further governed the spatial differences of PSRP, where open fractures blocked the upward transmission of P-holes, affecting the response of P-holes in PSRP. The revealed relationship between PSRP and crack development/distribution provides novel insights for monitoring and early warning of engineering rock mass instability and earthquakes.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106328"},"PeriodicalIF":7.5,"publicationDate":"2025-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145461551","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-07DOI: 10.1016/j.ijrmms.2025.106314
Huchen Duan , Huaning Wang , Mingjing Jiang , Fei Song , Yongheng Yang
Geomaterials generally exhibit time-dependency and plasticity, meanwhile, significantly influenced by the loading-unloading paths and temperature settings. Furthermore, the mechanical properties heavily depend on temperatures. In this study, an alternative viscoelastic-plastic analytical solution for stress and displacements is developed for spherical cavity closure, innovatively taking into account the time-dependent plastic behaviours of geomaterials, the loading-unloading paths, and the unsteady heat conduction and its effect on the mechanical behaviours.
In order to do that, the Burgers viscoelastic model and plastic slider characterised by the Mohr-Coulomb model are selected to characterise the time-dependency and plasticity of geomaterials, respectively. As a verification and validation step, the developed analytical solutions are compared with numerical predictions and experimental tests. After that, comprehensive parametric analyses are performed and some significant conclusions are achieved: (1) The critical supporting pressure for cavity stability increases linearly with both temperature and elastic modulus of rock. (2) Moderate heating is beneficial for stability, while excessive heating may cause secondary plastic yielding. The cavity reaches its most dangerous state during the early stages of storage. (3) The thermal conductivity and viscosity coefficient primarily affect the transient stress paths and stability of the cavity but do not impact the initial or final stress states. The developed time-dependent plastic analytical solution provides an alternative and efficient tool with a high potential for application to several relevant case studies, such as nuclear waste storage constructed in salt rocks.
{"title":"Time-dependent plastic analytical solutions for spherical cavity considering loading-unloading paths and temperature effects","authors":"Huchen Duan , Huaning Wang , Mingjing Jiang , Fei Song , Yongheng Yang","doi":"10.1016/j.ijrmms.2025.106314","DOIUrl":"10.1016/j.ijrmms.2025.106314","url":null,"abstract":"<div><div>Geomaterials generally exhibit time-dependency and plasticity, meanwhile, significantly influenced by the loading-unloading paths and temperature settings. Furthermore, the mechanical properties heavily depend on temperatures. In this study, an alternative viscoelastic-plastic analytical solution for stress and displacements is developed for spherical cavity closure, innovatively taking into account the time-dependent plastic behaviours of geomaterials, the loading-unloading paths, and the unsteady heat conduction and its effect on the mechanical behaviours.</div><div>In order to do that, the Burgers viscoelastic model and plastic slider characterised by the Mohr-Coulomb model are selected to characterise the time-dependency and plasticity of geomaterials, respectively. As a verification and validation step, the developed analytical solutions are compared with numerical predictions and experimental tests. After that, comprehensive parametric analyses are performed and some significant conclusions are achieved: (1) The critical supporting pressure for cavity stability increases linearly with both temperature and elastic modulus of rock. (2) Moderate heating is beneficial for stability, while excessive heating may cause secondary plastic yielding. The cavity reaches its most dangerous state during the early stages of storage. (3) The thermal conductivity and viscosity coefficient primarily affect the transient stress paths and stability of the cavity but do not impact the initial or final stress states. The developed time-dependent plastic analytical solution provides an alternative and efficient tool with a high potential for application to several relevant case studies, such as nuclear waste storage constructed in salt rocks.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106314"},"PeriodicalIF":7.5,"publicationDate":"2025-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145462347","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-07DOI: 10.1016/j.ijrmms.2025.106324
Jung-Wook Park , Chan-Hee Park , Eui-Seob Park , Olaf Kolditz
This study presents the development, verification, and validation of a novel coupling methodology for simulating coupled hydro-mechanical (HM) processes in fractured geological media. The proposed numerical approach, OGS–3DEC, integrates the finite element method-based OpenGeoSys (OGS) for fluid flow and the discrete element code 3DEC for mechanical deformation. These two simulators are linked through a hierarchical sequential coupling scheme that enables two-way data exchange. This approach modifies the effective stress in the fractured medium through evolving pore pressure while updating permeability and porosity based on stress and deformation. The methodology was verified through benchmark tests including one-dimensional consolidation in a porous matrix, radial fluid flow along a fracture plane, and fluid injection-induced fracture opening and slip. The results demonstrated excellent agreement with analytical solutions. Additionally, the model was validated against field-scale data from a fault reactivation experiment at the Mont Terri Rock Laboratory, with simulated pressures, flow rates, and displacements closely matching observed values. The OGS-3DEC simulator provides a reliable and efficient tool for modeling the HM behavior of fractured rock masses, with applications in geotechnical and subsurface engineering.
{"title":"A novel FEM–DEM coupling methodology for hydro-mechanical modeling in fractured geological media","authors":"Jung-Wook Park , Chan-Hee Park , Eui-Seob Park , Olaf Kolditz","doi":"10.1016/j.ijrmms.2025.106324","DOIUrl":"10.1016/j.ijrmms.2025.106324","url":null,"abstract":"<div><div>This study presents the development, verification, and validation of a novel coupling methodology for simulating coupled hydro-mechanical (HM) processes in fractured geological media. The proposed numerical approach, OGS–3DEC, integrates the finite element method-based OpenGeoSys (OGS) for fluid flow and the discrete element code 3DEC for mechanical deformation. These two simulators are linked through a hierarchical sequential coupling scheme that enables two-way data exchange. This approach modifies the effective stress in the fractured medium through evolving pore pressure while updating permeability and porosity based on stress and deformation. The methodology was verified through benchmark tests including one-dimensional consolidation in a porous matrix, radial fluid flow along a fracture plane, and fluid injection-induced fracture opening and slip. The results demonstrated excellent agreement with analytical solutions. Additionally, the model was validated against field-scale data from a fault reactivation experiment at the Mont Terri Rock Laboratory, with simulated pressures, flow rates, and displacements closely matching observed values. The OGS-3DEC simulator provides a reliable and efficient tool for modeling the HM behavior of fractured rock masses, with applications in geotechnical and subsurface engineering.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106324"},"PeriodicalIF":7.5,"publicationDate":"2025-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145461552","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-06DOI: 10.1016/j.ijrmms.2025.106325
Shun Liang , Hongye Luo , Derek Elsworth , Qiangling Yao , Xuehai Fu , Qiang Wang , Xuehua Li , Weisheng He , Zhi Ma , Guangli Huang , Furong Wang
In inter-stratified coal-petroleum basins, the integrity of oil/gas wells penetrating coal seams is frequently compromised by longwall coal mining-induced disturbances, promoting leakage of highly toxic hydrogen sulfide (H2S) gas from below. This released H2S migrates via mining-induced fractures into underground workings like longwall faces and roadways, posing acute exposure risks for miners. This study investigates H2S leakage, migration, hazard, and mitigation methods using a case study of the Shuangma coal mine (Ordos Basin western margin, China), where upper coal seams are mined above deep oil reservoirs. We developed a novel coupled mechanical-hydraulic-chemical model simulating H2S migration and coal seam alkali injection for sulfur immobilization. Field measurements and simulations reveal that: (1) Leaking H2S primarily adsorbs in coal matrix pores, forming enrichment zones until saturation, then distributes into free and water-soluble states in dynamic equilibrium; (2) H2S pressure increases exponentially near wellbores, with radii of influence after 30 years measuring 238–536 m for wellbore pressures of 0.32–1.52 MPa; (3) Optimal alkali injection parameters for effective H2S mitigation are 10 MPa pressure, 10 m borehole spacing, and 30 h grouting duration. These parameters suppressed H2S concentrations below the safety threshold. The results: (1) elucidate coupled transport-immobilization mechanisms governing H2S behavior in fractured coal-reservoir systems, and (2) provide a validated engineering protocol for abandoned well remediation in inter-stratified coal-hydrocarbon basins. This work advances fundamental understanding and practical solutions for H2S risk management in mining overlying oil/gas reservoirs.
{"title":"Migration mechanisms of leaking hydrogen sulfide in inter-stratified coal-petroleum basins and mitigation with alkali injection","authors":"Shun Liang , Hongye Luo , Derek Elsworth , Qiangling Yao , Xuehai Fu , Qiang Wang , Xuehua Li , Weisheng He , Zhi Ma , Guangli Huang , Furong Wang","doi":"10.1016/j.ijrmms.2025.106325","DOIUrl":"10.1016/j.ijrmms.2025.106325","url":null,"abstract":"<div><div>In inter-stratified coal-petroleum basins, the integrity of oil/gas wells penetrating coal seams is frequently compromised by longwall coal mining-induced disturbances, promoting leakage of highly toxic hydrogen sulfide (H<sub>2</sub>S) gas from below. This released H<sub>2</sub>S migrates via mining-induced fractures into underground workings like longwall faces and roadways, posing acute exposure risks for miners. This study investigates H<sub>2</sub>S leakage, migration, hazard, and mitigation methods using a case study of the Shuangma coal mine (Ordos Basin western margin, China), where upper coal seams are mined above deep oil reservoirs. We developed a novel coupled mechanical-hydraulic-chemical model simulating H<sub>2</sub>S migration and coal seam alkali injection for sulfur immobilization. Field measurements and simulations reveal that: (1) Leaking H<sub>2</sub>S primarily adsorbs in coal matrix pores, forming enrichment zones until saturation, then distributes into free and water-soluble states in dynamic equilibrium; (2) H<sub>2</sub>S pressure increases exponentially near wellbores, with radii of influence after 30 years measuring 238–536 m for wellbore pressures of 0.32–1.52 MPa; (3) Optimal alkali injection parameters for effective H<sub>2</sub>S mitigation are 10 MPa pressure, 10 m borehole spacing, and 30 h grouting duration. These parameters suppressed H<sub>2</sub>S concentrations below the safety threshold. The results: (1) elucidate coupled transport-immobilization mechanisms governing H<sub>2</sub>S behavior in fractured coal-reservoir systems, and (2) provide a validated engineering protocol for abandoned well remediation in inter-stratified coal-hydrocarbon basins. This work advances fundamental understanding and practical solutions for H<sub>2</sub>S risk management in mining overlying oil/gas reservoirs.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106325"},"PeriodicalIF":7.5,"publicationDate":"2025-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145461558","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-05DOI: 10.1016/j.ijrmms.2025.106320
Yeonkyeong Lee, Insun Song, Chan Park
This paper presents an analytical assessment of the CO2-storage efficiency and geomechanical responses to pore pressure buildup. The analytical assessment was based on the results of poromechanical experiments conducted on sandstone and shale samples from a prospective CO2 geological-storage site off the southeastern coast of Korea. Pore compressibility, bulk compressibility, and specific storage coefficients were determined using a custom-designed hydraulic test system under varying pore and confining pressures. Static and dynamic rock properties were further characterized through uniaxial compression tests and ultrasonic velocity measurements. The data from these methods were integrated to evaluate CO2-storage efficiency, rock deformability, and stress changes under different stress conditions. We found pronounced nonlinearity and hysteresis during loading and unloading cycles, owing predominantly to microcrack activity in sandstone and to hydration/dehydration effects in shale. Stress-dependent mechanical changes underscore the need to consider depth-specific conditions and cyclic loading during CO2 injection. Under semi-closed boundary conditions, we estimated a storage-efficiency coefficient of 0.8–1.3 % and approximately 0.06 % vertical deformation for 10 MPa pore pressure buildup at 2,700 m depth. Consistent poromechanical coupling across methods insights for optimizing storage capacity and refining geomechanical models to ensure the long-term storage stability and effectiveness.
{"title":"Analytical assessment of CO2-storage efficiency and geomechanical-coupling behaviors using poromechanical experiments","authors":"Yeonkyeong Lee, Insun Song, Chan Park","doi":"10.1016/j.ijrmms.2025.106320","DOIUrl":"10.1016/j.ijrmms.2025.106320","url":null,"abstract":"<div><div>This paper presents an analytical assessment of the CO<sub>2-</sub>storage efficiency and geomechanical responses to pore pressure buildup. The analytical assessment was based on the results of poromechanical experiments conducted on sandstone and shale samples from a prospective CO<sub>2</sub> geological-storage site off the southeastern coast of Korea. Pore compressibility, bulk compressibility, and specific storage coefficients were determined using a custom-designed hydraulic test system under varying pore and confining pressures. Static and dynamic rock properties were further characterized through uniaxial compression tests and ultrasonic velocity measurements. The data from these methods were integrated to evaluate CO<sub>2</sub>-storage efficiency, rock deformability, and stress changes under different stress conditions. We found pronounced nonlinearity and hysteresis during loading and unloading cycles, owing predominantly to microcrack activity in sandstone and to hydration/dehydration effects in shale. Stress-dependent mechanical changes underscore the need to consider depth-specific conditions and cyclic loading during CO<sub>2</sub> injection. Under semi-closed boundary conditions, we estimated a storage-efficiency coefficient of 0.8–1.3 % and approximately 0.06 % vertical deformation for 10 MPa pore pressure buildup at 2,700 m depth. Consistent poromechanical coupling across methods insights for optimizing storage capacity and refining geomechanical models to ensure the long-term storage stability and effectiveness.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106320"},"PeriodicalIF":7.5,"publicationDate":"2025-11-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145447342","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-04DOI: 10.1016/j.ijrmms.2025.106301
Zheng Teng , Hui Wu , Jize Zhang , Xin Ju , Shengwen Qi
Characterization of discrete fracture networks (DFNs) in the shallow crust is essential for understanding subsurface flow and transport processes and guiding reservoir exploitation such as water/oil/gas/geothermal/mineral recovery and nuclear waste/CO2 storage. However, characterizing the geometry of subsurface DFNs is extremely difficult, due to the inherent complexity of DFNs and the generally spatially sparse, low-resolution geological/geophysical data. Traditional DFN parameterization methods may result in a high-dimensional parameter space, making DFN inversion ill-posed and computationally expensive. In this study, we develop a deep learning-based low-dimensional parameterization method to effectively generate complex DFNs from low-dimensional latent spaces, thus significantly alleviating the ill-posedness and computational burden associated with DFN characterization in a data scarce environment. The Wasserstein generative adversarial network with gradient penalty (WGAN-GP) is used to generate random DFNs from latent spaces. Through both qualitative and quantitative comparisons of fracture characteristics between the generated and training DFNs, we demonstrate the extraordinary capability of the method in generating high-fidelity DFNs from extremely low-dimensional latent spaces. The generated DFNs faithfully honor fracture prior knowledge imposed in training samples, including fracture statistics regarding location, length and orientation as well as fracture existence and connectivity identified from geological/hydrogeological surveys. We also demonstrate the ability of the method in generating DFNs that resemble realistic fracture networks mapped from limestone and glacier outcrops. A synthetic DFN characterization case study illustrates the effectiveness of the proposed method in inversion tasks, showing such an effective low-dimensional and conditional parameterization method is particularly useful to facilitate subsurface DFN characterization.
{"title":"Generating high-fidelity discrete fracture networks from low-dimensional latent spaces using generative adversarial network","authors":"Zheng Teng , Hui Wu , Jize Zhang , Xin Ju , Shengwen Qi","doi":"10.1016/j.ijrmms.2025.106301","DOIUrl":"10.1016/j.ijrmms.2025.106301","url":null,"abstract":"<div><div>Characterization of discrete fracture networks (DFNs) in the shallow crust is essential for understanding subsurface flow and transport processes and guiding reservoir exploitation such as water/oil/gas/geothermal/mineral recovery and nuclear waste/CO<sub>2</sub> storage. However, characterizing the geometry of subsurface DFNs is extremely difficult, due to the inherent complexity of DFNs and the generally spatially sparse, low-resolution geological/geophysical data. Traditional DFN parameterization methods may result in a high-dimensional parameter space, making DFN inversion ill-posed and computationally expensive. In this study, we develop a deep learning-based low-dimensional parameterization method to effectively generate complex DFNs from low-dimensional latent spaces, thus significantly alleviating the ill-posedness and computational burden associated with DFN characterization in a data scarce environment. The Wasserstein generative adversarial network with gradient penalty (WGAN-GP) is used to generate random DFNs from latent spaces. Through both qualitative and quantitative comparisons of fracture characteristics between the generated and training DFNs, we demonstrate the extraordinary capability of the method in generating high-fidelity DFNs from extremely low-dimensional latent spaces. The generated DFNs faithfully honor fracture prior knowledge imposed in training samples, including fracture statistics regarding location, length and orientation as well as fracture existence and connectivity identified from geological/hydrogeological surveys. We also demonstrate the ability of the method in generating DFNs that resemble realistic fracture networks mapped from limestone and glacier outcrops. A synthetic DFN characterization case study illustrates the effectiveness of the proposed method in inversion tasks, showing such an effective low-dimensional and conditional parameterization method is particularly useful to facilitate subsurface DFN characterization.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106301"},"PeriodicalIF":7.5,"publicationDate":"2025-11-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145441365","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-03DOI: 10.1016/j.ijrmms.2025.106323
Xiaotian Wu , Yingchun Li , Tianjiao Li , Guanglei Cui , Chun'an Tang
Hydraulic fracturing is critical to considerably elevate the permeability of the hot dry rock in an enhanced geothermal system. However, the real-time evolution of multi-physical fields over hydraulic fracturing with varying geomechanical properties and injection strategies remains unclear. Here, we established an integrated numerical model that fully couples thermal, hydraulic and mechanical fields with damage (THMD) to simulate hydraulic fracturing in a heterogeneous geothermal formation. We tracked the fracture initiation and propagation via quantitative rock damage detected by the maximum shear/tensile stress criterion, and examined the evolution of temperature, pore pressure and stress. Our numerical simulations revealed that fluid injection first creates uniform diffusion of cooling and elevates the pore pressure around the borehole with increasing minimum principal stress. Subsequently, fluid migration becomes preferentially channeled through newly formed fractures and thus cooling and high-pressure zones are concentrated along these highly permeable pathways. Higher principal stress difference, formation temperature, and injection rate all promote fracture propagation and decrease initiation pressure. Conversely, a greater rock permeability hinders fracture growth and lowers initiation pressure. The quantified sensitivity analysis shows that the principal stress difference primarily impacts the fracture initiation pressure, and the injection rate predominantly dictates the fracture length. The rock heterogeneity complicates the fracture morphology by introducing additional bifurcations but imposes negligible influence on the fracture length and initiation pressure. These insights deliver practical guidance to optimize operational parameters in geothermal reservoir stimulation and fracture propagation modulation.
{"title":"Numerical insights into thermal-hydraulic-mechanical-damage (THMD) fields evolution over hydraulic fracturing in hot dry rock","authors":"Xiaotian Wu , Yingchun Li , Tianjiao Li , Guanglei Cui , Chun'an Tang","doi":"10.1016/j.ijrmms.2025.106323","DOIUrl":"10.1016/j.ijrmms.2025.106323","url":null,"abstract":"<div><div>Hydraulic fracturing is critical to considerably elevate the permeability of the hot dry rock in an enhanced geothermal system. However, the real-time evolution of multi-physical fields over hydraulic fracturing with varying geomechanical properties and injection strategies remains unclear. Here, we established an integrated numerical model that fully couples thermal, hydraulic and mechanical fields with damage (THMD) to simulate hydraulic fracturing in a heterogeneous geothermal formation. We tracked the fracture initiation and propagation via quantitative rock damage detected by the maximum shear/tensile stress criterion, and examined the evolution of temperature, pore pressure and stress. Our numerical simulations revealed that fluid injection first creates uniform diffusion of cooling and elevates the pore pressure around the borehole with increasing minimum principal stress. Subsequently, fluid migration becomes preferentially channeled through newly formed fractures and thus cooling and high-pressure zones are concentrated along these highly permeable pathways. Higher principal stress difference, formation temperature, and injection rate all promote fracture propagation and decrease initiation pressure. Conversely, a greater rock permeability hinders fracture growth and lowers initiation pressure. The quantified sensitivity analysis shows that the principal stress difference primarily impacts the fracture initiation pressure, and the injection rate predominantly dictates the fracture length. The rock heterogeneity complicates the fracture morphology by introducing additional bifurcations but imposes negligible influence on the fracture length and initiation pressure. These insights deliver practical guidance to optimize operational parameters in geothermal reservoir stimulation and fracture propagation modulation.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106323"},"PeriodicalIF":7.5,"publicationDate":"2025-11-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145428796","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Chemical degradation of bonded granular materials sustaining mechanical loads is a critical process governing the long-term stability of natural geological systems and the safety of subsurface energy engineering operations. To investigate the interplays controlling reactive cracking in geomaterials, this study develops a multi-scale model that couples the phase-field fracture mechanics with damage-enhanced reactive diffusion. The novelty lies in a micromechanically-derived degradation law embedded in the phase-field fracture model, where material properties are informed by Discrete Element Method (DEM) simulations of intergranular bond dissolution. Our results show that a higher initial cementation level substantially postpones the triggering of the brittle fracture, as more time is required for the accumulated mass removal to cause the critical amount of degradation in the matrix ahead of the crack tip. A chemical ductilization effect is identified upon the onset of fracturing, where a higher environmental acidity counter-intuitively results in a delay in the initiation of the brittle fracture. This phenomenon is attributed to the acidity-enhanced softening around the crack tip, leading to a less steep increase in the maximum tensile stress perpendicular to the crack growth direction, which reaches the peak at a lower value. Furthermore, a competition between the environmentally induced softening effect that delays the onset of fracturing and a direct degradation of the intrinsic fracture toughness that promotes it, is quantitatively illustrated. The time required for inducing the brittle fracturing in a chemically degrading carbonate-cemented geomaterial shows a linear dependence on the material’s critical energy release rate, under the dynamic interplay between the two counteracting mechanisms. The findings are applicable to a broad category of subsurface engineering concerning geomechanics, providing a fundamental basis for assessing the long-term integrity and maintenance of geostructures subject to reactive environments. The findings could be calibrated through controlled laboratory fracturing experiments equipped with acoustic emission detection.
{"title":"Onset of reactive brittle cracking in sandstones: DEM-informed phase-field modeling","authors":"Fanyu Wu , Alexandre Sac-Morane , Hadrien Rattez , Manolis Veveakis , Manman Hu","doi":"10.1016/j.ijrmms.2025.106319","DOIUrl":"10.1016/j.ijrmms.2025.106319","url":null,"abstract":"<div><div>Chemical degradation of bonded granular materials sustaining mechanical loads is a critical process governing the long-term stability of natural geological systems and the safety of subsurface energy engineering operations. To investigate the interplays controlling reactive cracking in geomaterials, this study develops a multi-scale model that couples the phase-field fracture mechanics with damage-enhanced reactive diffusion. The novelty lies in a micromechanically-derived degradation law embedded in the phase-field fracture model, where material properties are informed by Discrete Element Method (DEM) simulations of intergranular bond dissolution. Our results show that a higher initial cementation level substantially postpones the triggering of the brittle fracture, as more time is required for the accumulated mass removal to cause the critical amount of degradation in the matrix ahead of the crack tip. A chemical ductilization effect is identified upon the onset of fracturing, where a higher environmental acidity counter-intuitively results in a delay in the initiation of the brittle fracture. This phenomenon is attributed to the acidity-enhanced softening around the crack tip, leading to a less steep increase in the maximum tensile stress perpendicular to the crack growth direction, which reaches the peak at a lower value. Furthermore, a competition between the environmentally induced softening effect that delays the onset of fracturing and a direct degradation of the intrinsic fracture toughness that promotes it, is quantitatively illustrated. The time required for inducing the brittle fracturing in a chemically degrading carbonate-cemented geomaterial shows a linear dependence on the material’s critical energy release rate, under the dynamic interplay between the two counteracting mechanisms. The findings are applicable to a broad category of subsurface engineering concerning geomechanics, providing a fundamental basis for assessing the long-term integrity and maintenance of geostructures subject to reactive environments. The findings could be calibrated through controlled laboratory fracturing experiments equipped with acoustic emission detection.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"196 ","pages":"Article 106319"},"PeriodicalIF":7.5,"publicationDate":"2025-11-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145428797","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}