首页 > 最新文献

International Journal of Rock Mechanics and Mining Sciences最新文献

英文 中文
Three-dimensional discrete fracture network identification based on deep learning and reversible jump Markov chain Monte Carlo algorithm 基于深度学习和可逆跳跃马尔可夫链蒙特卡罗算法的三维离散断裂网络识别
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-10 DOI: 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.
裂缝网络的特征对地下水开发、地热开发、油气开采和地质二氧化碳封存至关重要,但其复杂和不确定的空间分布给人们带来了持续的挑战。本研究提出了一种集成3D-UNet代理模型、可逆跳跃马尔可夫链蒙特卡罗(rjMCMC)和多源数据融合的智能反演框架,用于现场尺度的三维离散裂缝网络(DFN)表征。在此框架内,3D-UNet模型在裂缝配置、水头和电势的大型数据集上进行训练,提供了裂缝参数的有效初始反演。然后使用RANSAC算法提取裂缝几何形状,并通过rjMCMC进行迭代优化,其中代理3D-UNet取代了传统的正演建模。这一创新将计算成本降低了一个数量级,实现了高效的大规模反演。此外,电势与水头数据的融合使反演精度提高了约10%。验证表明,该框架可靠地重建了裂缝网络的空间分布,既捕获了低密度区域,也捕获了高度非均质区域的主要水力路径。该方法结合了计算效率和更高的精度,为广泛的水文地质和工程应用中的现场裂缝网络表征提供了一种实用且可扩展的解决方案。
{"title":"Three-dimensional discrete fracture network identification based on deep learning and reversible jump Markov chain Monte Carlo algorithm","authors":"Zhenting Sun ,&nbsp;Lei Ma ,&nbsp;Quan Li ,&nbsp;Yaping Deng ,&nbsp;Han Qiu ,&nbsp;Haichun Ma ,&nbsp;Cihai Chen ,&nbsp;Yongshuai Yan ,&nbsp;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}
引用次数: 0
Creep behavior and prestress relaxation mechanism of bolt-reinforced jointed specimen disturbed by dynamic impact 动冲击扰动下锚杆加固节理试件蠕变行为及预应力松弛机制
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-09 DOI: 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 ,&nbsp;Runze Zhu ,&nbsp;Ignacio Pérez-Rey ,&nbsp;Wancheng Zhu ,&nbsp;Xige Liu ,&nbsp;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}
引用次数: 0
Characteristics and spatial differences of pressure-stimulated rock potential in diorite partly loaded to fracturing 部分压裂加载闪长岩压激岩势特征及空间差异
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-08 DOI: 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.
压力激岩势(PSRP)的空间差异及其与裂缝发育和分布的相关性,尽管在认识地震异常和监测工程岩体方面具有重要意义,但研究较少。在本研究中,对特殊设计的闪长岩试样施加部分压应力,使其呈长方体-台地复合形状,使试样由加载长方体和自由台地组成。分别监测自由平台与加载长方体两端的PSRP差异,同时监测声发射,观察裂纹发展情况。结果表明,PSRP具有明显的阶段性特征:初始加载时呈指数下降,局部宏观破坏时快速上升,失稳破坏时呈脉冲状,峰后加载时出现大幅波动上升。在失稳失效阶段之前,不同区域PSRPs的变化趋势一致,但在幅度和时间上存在差异。然而,在不稳定失效阶段,出现了显著的空间差异,观察到相反的极性脉冲。检测到的PSRP可归因于流体电动力学效应和正孔(p孔)的激活。PSRP与裂缝发育相对应,裂缝累积面积越大,激活的p孔越多,PSRP上升越高。宏观裂缝面的空间分布进一步控制了PSRP的空间差异,开放裂缝阻断了p孔的向上传播,影响了PSRP中p孔的响应。揭示了PSRP与裂缝发育/分布的关系,为工程岩体失稳和地震监测预警提供了新的思路。
{"title":"Characteristics and spatial differences of pressure-stimulated rock potential in diorite partly loaded to fracturing","authors":"Licheng Sun ,&nbsp;Lixin Wu ,&nbsp;Youyou Xu ,&nbsp;Guangrui Dong ,&nbsp;Tao Zheng ,&nbsp;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}
引用次数: 0
Time-dependent plastic analytical solutions for spherical cavity considering loading-unloading paths and temperature effects 考虑加载-卸载路径和温度影响的球腔时效塑性解析解
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-07 DOI: 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.
岩土材料普遍表现出时间依赖性和可塑性,同时受加载-卸载路径和温度设置的显著影响。此外,机械性能很大程度上取决于温度。在这项研究中,开发了一种用于球形空腔闭合的应力和位移的粘弹塑性替代解析解,创新地考虑了岩土材料的时间依赖塑性行为、加载-卸载路径、非稳态热传导及其对力学行为的影响。为了做到这一点,选择Burgers粘弹性模型和以Mohr-Coulomb模型为特征的塑性滑块来分别表征岩土材料的时间依赖性和可塑性。作为验证和验证步骤,将开发的解析解与数值预测和实验测试进行了比较。在此基础上,进行了全面的参数分析,得到了一些有意义的结论:(1)空腔稳定临界支撑压力随温度和岩石弹性模量的增加而线性增加。(2)适度加热有利于稳定性,过度加热可能导致二次塑性屈服。在储存的早期阶段,蛀洞达到了最危险的状态。(3)导热系数和粘滞系数主要影响腔体的瞬态应力路径和稳定性,而不影响初始和最终应力状态。所开发的随时间变化的塑性分析解决方案提供了一种替代的高效工具,具有很高的应用潜力,可用于几个相关案例研究,例如在盐岩中建造的核废料储存。
{"title":"Time-dependent plastic analytical solutions for spherical cavity considering loading-unloading paths and temperature effects","authors":"Huchen Duan ,&nbsp;Huaning Wang ,&nbsp;Mingjing Jiang ,&nbsp;Fei Song ,&nbsp;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}
引用次数: 0
A novel FEM–DEM coupling methodology for hydro-mechanical modeling in fractured geological media 裂缝性地质介质中流体力学建模的一种新型FEM-DEM耦合方法
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-07 DOI: 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.
本研究提出了一种新的耦合方法,用于模拟裂缝性地质介质中耦合的水-力学过程,并进行了开发、验证和验证。所提出的数值方法OGS - 3DEC集成了基于有限元方法的流体流动OpenGeoSys (OGS)和机械变形离散单元代码3DEC。这两个模拟器通过分层顺序耦合方案连接起来,从而实现双向数据交换。该方法通过孔隙压力的变化来修正裂缝介质中的有效应力,同时根据应力和变形来修正渗透率和孔隙度。该方法通过基准测试进行了验证,包括多孔基质中的一维固结、沿裂缝平面的径向流体流动以及流体注入引起的裂缝张开和滑移。结果与解析解非常吻合。此外,该模型还通过Mont Terri岩石实验室断层再激活实验的现场数据进行了验证,模拟的压力、流速和排量与观测值非常吻合。OGS-3DEC模拟器为模拟裂隙岩体的HM行为提供了可靠、高效的工具,可用于岩土工程和地下工程。
{"title":"A novel FEM–DEM coupling methodology for hydro-mechanical modeling in fractured geological media","authors":"Jung-Wook Park ,&nbsp;Chan-Hee Park ,&nbsp;Eui-Seob Park ,&nbsp;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}
引用次数: 0
Migration mechanisms of leaking hydrogen sulfide in inter-stratified coal-petroleum basins and mitigation with alkali injection 层间油气盆地泄漏硫化氢运移机理及注碱治理
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-06 DOI: 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.
在层间煤-油气盆地中,穿透煤层的油气井的完整性经常受到长壁采煤引起的扰动的破坏,促进了高毒性硫化氢(H2S)气体从下面泄漏。这些释放出的H2S通过开采引起的裂缝进入地下工作,如长壁工作面和巷道,给矿工带来了严重的暴露风险。本研究以鄂尔多斯盆地西缘双马煤矿为例,研究了H2S的泄漏、迁移、危害和缓解方法。建立了模拟H2S运移和煤层注碱固硫的力学-水力-化学耦合模型。现场实测和模拟结果表明:(1)泄漏的H2S主要吸附在煤基质孔隙中,形成富集带直至饱和,然后在动态平衡状态下分布为游离态和水溶性;(2)井筒附近H2S压力呈指数增长,井筒压力为0.32 ~ 1.52 MPa时,30年后影响半径为238 ~ 536 m;(3)有效减缓H2S的最佳注碱参数为压力10 MPa、井距10 m、注浆时间30 h。这些参数将H2S浓度抑制在安全阈值以下。研究结果:(1)阐明了裂缝性煤储层系统中控制H2S行为的耦合运输-固定机制;(2)为层间煤-烃盆地弃井修复提供了一种有效的工程方案。这项工作促进了对油气储层开采中H2S风险管理的基本认识和实际解决方案。
{"title":"Migration mechanisms of leaking hydrogen sulfide in inter-stratified coal-petroleum basins and mitigation with alkali injection","authors":"Shun Liang ,&nbsp;Hongye Luo ,&nbsp;Derek Elsworth ,&nbsp;Qiangling Yao ,&nbsp;Xuehai Fu ,&nbsp;Qiang Wang ,&nbsp;Xuehua Li ,&nbsp;Weisheng He ,&nbsp;Zhi Ma ,&nbsp;Guangli Huang ,&nbsp;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}
引用次数: 0
Analytical assessment of CO2-storage efficiency and geomechanical-coupling behaviors using poromechanical experiments 基于孔隙力学实验的co2封存效率及地质力学耦合行为分析评价
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-05 DOI: 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.
本文分析了孔隙压力累积对co2储存效率和地质力学响应的影响。分析评估是基于对韩国东南海岸一个潜在的二氧化碳地质储存地点的砂岩和页岩样品进行的孔隙力学实验结果。在不同的孔隙压力和围压下,使用定制的水力测试系统来确定孔隙压缩性、体积压缩性和特定存储系数。通过单轴压缩试验和超声测速进一步表征了岩石的静态和动态特性。综合这些方法的数据,评估不同应力条件下的co2储存效率、岩石变形能力和应力变化。我们发现在加载和卸载循环过程中存在明显的非线性和滞后,这主要是由于砂岩中的微裂缝活动和页岩中的水化/脱水效应。应力相关的机械变化强调了在二氧化碳注入过程中需要考虑特定深度条件和循环载荷。在半封闭边界条件下,我们估计在2,700 m深度积累10 MPa孔隙压力时,存储效率系数为0.8 - 1.3%,垂直变形约为0.06%。不同方法间一致的孔隙力学耦合有助于优化储层容量,完善地质力学模型,确保长期储层的稳定性和有效性。
{"title":"Analytical assessment of CO2-storage efficiency and geomechanical-coupling behaviors using poromechanical experiments","authors":"Yeonkyeong Lee,&nbsp;Insun Song,&nbsp;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}
引用次数: 0
Generating high-fidelity discrete fracture networks from low-dimensional latent spaces using generative adversarial network 利用生成对抗网络从低维潜在空间生成高保真离散断裂网络
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-04 DOI: 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.
浅层地壳离散裂缝网络(DFNs)的表征对于理解地下流动和输送过程以及指导储层开发(如水/油/气/地热/矿物回收和核废料/二氧化碳储存)至关重要。然而,由于地下DFNs固有的复杂性以及通常空间稀疏、低分辨率的地质/地球物理数据,表征地下DFNs的几何形状是非常困难的。传统的DFN参数化方法可能导致参数空间高维,使得DFN反演不适定且计算量大。在本研究中,我们开发了一种基于深度学习的低维参数化方法,从低维潜在空间中有效地生成复杂的DFN,从而显著减轻了在数据稀缺环境中与DFN表征相关的不适定性和计算负担。采用Wasserstein梯度惩罚生成对抗网络(WGAN-GP)从潜在空间生成随机dfn。通过对生成和训练DFNs之间的断裂特征进行定性和定量比较,我们证明了该方法在从极低维潜在空间生成高保真DFNs方面的非凡能力。生成的dfn忠实地遵循了训练样本中的裂缝先验知识,包括有关裂缝位置、长度和方向的统计数据,以及从地质/水文地质调查中确定的裂缝存在性和连通性。我们还证明了该方法在生成类似于石灰石和冰川露头绘制的真实裂缝网络的dfn方面的能力。一个综合DFN表征案例研究说明了该方法在反演任务中的有效性,表明这种有效的低维条件参数化方法对于促进地下DFN表征特别有用。
{"title":"Generating high-fidelity discrete fracture networks from low-dimensional latent spaces using generative adversarial network","authors":"Zheng Teng ,&nbsp;Hui Wu ,&nbsp;Jize Zhang ,&nbsp;Xin Ju ,&nbsp;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}
引用次数: 0
Numerical insights into thermal-hydraulic-mechanical-damage (THMD) fields evolution over hydraulic fracturing in hot dry rock 热-水力-机械损伤(THMD)场在干热岩石水力压裂过程中的演化
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-03 DOI: 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.
在增强型地热系统中,水力压裂对于显著提高热干岩的渗透率至关重要。然而,在不同地质力学性质和注入策略的情况下,水力压裂过程中多物理场的实时演化仍不清楚。为了模拟非均质地热地层的水力压裂,建立了热、水力和力学场与损伤完全耦合的综合数值模型(THMD)。通过最大剪切/拉应力准则定量检测岩石损伤,跟踪裂缝的起裂和扩展过程,并研究了温度、孔隙压力和应力的演化过程。我们的数值模拟表明,流体注入首先产生均匀的冷却扩散,并提高井眼周围的孔隙压力,同时增加最小主应力。随后,流体运移优先通过新形成的裂缝,因此冷却区和高压区集中在这些高渗透性通道上。较高的主应力差、地层温度和注入速度均促进裂缝扩展,降低起裂压力。相反,较大的岩石渗透率会阻碍裂缝的扩展,降低起裂压力。定量敏感性分析表明,主应力差主要影响起裂压力,注入速度主要决定裂缝长度。岩石的非均质性通过引入额外的分叉使裂缝形态复杂化,但对裂缝长度和起裂压力的影响可以忽略不计。这些见解为优化地热储层增产和裂缝扩展调节的操作参数提供了实用指导。
{"title":"Numerical insights into thermal-hydraulic-mechanical-damage (THMD) fields evolution over hydraulic fracturing in hot dry rock","authors":"Xiaotian Wu ,&nbsp;Yingchun Li ,&nbsp;Tianjiao Li ,&nbsp;Guanglei Cui ,&nbsp;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}
引用次数: 0
Onset of reactive brittle cracking in sandstones: DEM-informed phase-field modeling 砂岩中反应性脆性开裂的开始:dem通知相场建模
IF 7.5 1区 工程技术 Q1 ENGINEERING, GEOLOGICAL Pub Date : 2025-11-03 DOI: 10.1016/j.ijrmms.2025.106319
Fanyu Wu , Alexandre Sac-Morane , Hadrien Rattez , Manolis Veveakis , Manman Hu
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.
承受机械载荷的粘结颗粒材料的化学降解是控制自然地质系统长期稳定和地下能源工程运行安全的关键过程。为了研究控制岩土材料反应性开裂的相互作用,本研究建立了一个将相场断裂力学与损伤增强反应扩散耦合的多尺度模型。其新颖之处在于相场断裂模型中嵌入了微观力学推导的降解规律,其中材料性能是通过离散元法(DEM)模拟晶间键溶解得到的。我们的研究结果表明,较高的初始胶结水平大大推迟了脆性断裂的触发,因为在裂纹尖端之前,积累的质量需要更多的时间来清除,从而导致基体中临界量的降解。在压裂开始时,发现了化学延展性效应,较高的环境酸度反而会延迟脆性裂缝的形成。这种现象是由于裂纹尖端周围的酸性增强软化,导致垂直于裂纹扩展方向的最大拉应力增加幅度较小,在较低的值处达到峰值。此外,还定量说明了环境诱导的软化效应(延迟压裂发生)与促进压裂发生的固有断裂韧性的直接退化之间的竞争。在两种抵消机制的动态相互作用下,化学降解碳酸盐胶结岩土材料发生脆性破裂所需的时间与材料的临界能量释放率呈线性关系。研究结果适用于地质力学的地下工程,为评估受反应环境影响的土工结构的长期完整性和维护提供了基本依据。这些发现可以通过配备声发射检测的受控实验室压裂实验进行校准。
{"title":"Onset of reactive brittle cracking in sandstones: DEM-informed phase-field modeling","authors":"Fanyu Wu ,&nbsp;Alexandre Sac-Morane ,&nbsp;Hadrien Rattez ,&nbsp;Manolis Veveakis ,&nbsp;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}
引用次数: 0
期刊
International Journal of Rock Mechanics and Mining Sciences
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1