• CSCD核心中文核心科技核心
  • RCCSE(A+)公路运输高质量期刊T1
  • Ei CompendexScopusWJCI
  • EBSCOPж(AJ)JST
二维码

隧道建设(中英文) ›› 2026, Vol. 46 ›› Issue (7): 1402-1412.DOI: 10.3973/j.issn.2096-4498.2026.07.002

• 研究与探索 • 上一篇    下一篇

基于非弹性断裂相场法的储气硐室围岩裂缝扩展规律研究

夏才初1,刘少华1, 李兴2, 彭奕亮2, 徐英俊1, 徐晨1   

  1. (1. 宁波大学岩石力学研究所 全省岩石力学与地质灾害重点实验室, 浙江 宁波 315211;2. 中国电建集团河南省电力勘测设计院有限公司, 河南 郑州 450007)
  • 出版日期:2026-07-20 发布日期:2026-07-20
  • 作者简介:夏才初(1963—),男,浙江萧山人,1992年毕业于中南大学,采矿工程专业,博士,教授,主要从事隧道与地下工程方面的研究工作。 E-mail: tjxiaccb@126.com。

Crack Propagation Patterns of Surrounding Rock in Gas Storage Caverns Based on Inelastic Phase-Field Method

XIA Caichu1, LIU Shaohua1, LI Xing2, PENG Yiliang2, XU Yingjun1, XU Chen1   

  1. (1. Zhejiang Key Laboratory of Rock Mechanics and Geohazards, Institute of Rock Mechanics, Ningbo University, Ningbo 315211, Zhejiang, China; 2. PowerChina Henan Electric Power Survey & Design Institute Co., Ltd., Zhengzhou 450007, Henan, China)
  • Online:2026-07-20 Published:2026-07-20

摘要: 在压气储能人工开挖的地下储气硐室中,围岩是承担高内压气体的主要载体; 在高内压作用下,围岩裂缝的扩展会威胁地下硐室的密封性和安全性。为研究高内压作用下储气硐室围岩裂缝的扩展规律,在经典的弹性断裂相场模型理论框架的基础上,建立非弹性断裂相场法。引入非弹性应变,以非弹性应变能密度函数作为相场损伤演化的驱动力,驱动裂缝的损伤演化。采用非弹性断裂相场法计算不对称缺口试样在拉应力作用下的裂缝扩展过程,并与真实的试验结果进行对比,验证非弹性断裂相场法在计算弹脆性材料裂缝扩展过程中的正确性。自然界中围岩是非均质的,基于Weibull分布函数引入非均质度参数建立储气硐室围岩的非均质分布模型,采用非弹性断裂相场法计算储气硐室围岩裂缝的分布。计算结果表明: 围岩非均质性会对裂缝扩展规律产生重要影响,同等内压条件下,围岩非均质性越强,越容易产生裂缝; 高内压条件下,非均质性强的围岩会产生更多的弥散性裂缝。因此,在压气储能地下储气硐室选址过程中,应选择围岩性质更为均质的地理环境; 在储气硐室的开挖阶段,应对硐室周围围岩进行灌浆处理,以加强硐室围岩强度和均质性,减少裂缝扩展。

关键词: 压气储能, 储气硐室, 相场法, 裂缝

Abstract: In underground caverns excavated for compressed air energy storage (CAES), the surrounding rock serves as the primary carrier of high-pressure gas. Under high internal pressure, crack propagation in the surrounding rock can compromise the tightness and structural integrity of underground caverns. To investigate crack propagation patterns of the surrounding rock of gas storage caverns under high internal pressure, an inelastic fracture phase-field method is developed based on the classical elastic fracture phase-field model. By introducing inelastic strain, the inelastic strain energy density function is adopted as the driving force for phase-field damage evolution, thereby governing crack damage development. The proposed method is applied to simulate crack propagation in an asymmetric notched specimen under tensile stress, and the simulation results are compared with experimental data to verify the accuracy of the method for calculating crack propagation in elastobrittle materials. As the natural surrounding rock is inherently heterogeneous, a heterogeneity parameter is introduced based on the Weibull distribution function to construct a heterogeneous distribution model for gas storage cavern surrounding rock. The inelastic phase-field method is then employed to calculate the crack distribution in the surrounding rock. The numerical results indicate that surrounding rock heterogeneity markedly influences crack propagation patterns. Under identical internal pressure, surrounding rock with stronger heterogeneity is more prone to crack initiation. Under high internal pressure, highly heterogeneous surrounding rock tends to generate more diffuse cracks. Accordingly, geographical environments with more homogeneous surrounding rock should be prioritized during site selection for underground CAES caverns. During cavern excavation, grouting treatment shall be applied to the surrounding rock to improve its strength and homogeneity and suppress crack propagation.

Key words: compressed air energy storage, gas storage caverns, phase field method, crack propagation