Study on the potential instability patterns of tunnel type underground caverns for compressed air energy storage
SUN Guanhua1，2，YI Qi1，2，YAO Yuanfeng3，SHANG Haoliang4，JI Wendong4
(1. State Key Laboratory of Geotechnical Mechanics and Engineering，Institute of Rock and Soil Mechanics，Chinese Academy of Sciences，Wuhan，Hubei 430071，China；2. University of Chinese Academy of Sciences，Beijing 100049，China；3. Central Southern China Electric Power Design Institute Co.，Ltd. of China Power Engineering Consulting Group，Wuhan，Hubei 430071，China；4. China Energy Digital Technology Group Co.，Ltd.，Beijing 100044，China)
Abstract：With the rapid development of new energy sources，the compressed air energy storage(CAES) underground caverns，which differ from traditional tunnel conditions，face stability issues in the overlying rock mass. The potential instability modes of these shallow-buried tunnel-type CAES caverns urgently need to be studied. Adopting the Mohr-Coulomb strength criterion and based on the ultimate stress field of the rock mass and considering both tensile and shear failure modes，the problem of potential instability modes in shallow buried high internal pressure tunnel-type lined rock caverns is reduced to the initial value problem of a cluster of ordinary differential equations，and the solution method is given. The reliability of the method is verified by comparing and analyzing the results with the physical model test results. In addition，the effects of burial depth，tunnel diameter，in-situ stress ratio and internal friction angle on the instability pattern are investigated，and the results show that the in-situ stress ratio and internal friction angle have a significant influence. The horizontal stress direction is perpendicular to the axis of the tunnel，which is more conducive to the stability of the overlying rock. On the other hand，the influence of the internal friction angle on the instability pattern depends on the in-situ stress ratio. When the in-situ stress ratio is less than 1，the smaller the internal friction angle is，the closer the potential rupture surface is to the horizontal；when the in-situ stress ratio is equal to 1，the internal friction angle has almost no effect；when the in-situ stress ratio is greater than 1，the larger the internal friction angle is，the closer the potential rupture surface is to the horizontal. In view of this，the design of cavern layout，burial depth and diameter should consider both in-situ stress and internal friction angle. Finally，potential failure surface morphology is parameterized in a simple geometric way，and the correspondence between the geometric parameters and the in-situ stress ratio and internal friction angle is given for engineering reference. In terms of failure modes，when the geostress coefficient is small，the potential failure surface exhibits a tensile-shear composite mode，with the tension region located near the cavern wall and the ground surface. When the geostress coefficient is large，the failure mode is primarily shear.
孙冠华1，2，易 琪1，2，姚院峰3，商浩亮4，纪文栋4. 压缩空气储能电站隧道式地下硐库潜在失稳模式研究[J]. 岩石力学与工程学报, 2024, 43(1): 41-49.
SUN Guanhua1，2，YI Qi1，2，YAO Yuanfeng3，SHANG Haoliang4，JI Wendong4. Study on the potential instability patterns of tunnel type underground caverns for compressed air energy storage. , 2024, 43(1): 41-49.
PERAZZELLI P，ANAGNOSTOU G. Design issues for compressed air energy storage in sealed underground cavities[J]. Journal of Rock Mechanics and Engineering，2016：8(3)：314-328.
[13]
郑 宏，刘德富，罗先启. 基于变形分析的边坡潜在滑面的确定[J]. 岩石力学与工程学报，2004，23(5)：709-716.(ZHENG Hong，LIU Defu，LUO Xianqi. Determination of potential slide line of slopes based on deformation analysis[J]. Chinese Journal of Rock Mechanics and Engineering，2004，23(5)：709-716.(in Chinese))
[1]
李仲奎，马芳平，刘 辉. 压气蓄能电站的地下工程问题及应用前景[J]. 岩石力学与工程学报，2003，22(增1)：2 121-2 126.(LI Zhongkui，MA Fangping，LIU Hui. Underground engineering problems in compressed air energy storage and its developing future[J]. Chinese Journal of Rock Mechanics and Engineering，2003，22(Supp.1)：2 121-2 126.(in Chinese))
[11]
GHALY A，HANNA A. Ultimate pullout resistance of single vertical anchors[J]. Canadian Geotechnical Journal，1994，31(5)：661-672.
[9]
徐英俊，夏才初，周舒威，等. 基于极限分析上限定理的压气储能洞室抗隆起破坏准则[J]. 岩石力学与工程学报，2022，41(10)：1 971-1 980.(XU Yingjun，XIA Caichu，ZHOU Shuwei，et al. Anti-uplift failure criterion of caverns for compressed air energy storage based on the upper bound theorem of limit analysis[J]. Chinese Journal of Rock Mechanics and Engineering，2022，41(10)：1 971-1 980.(in Chinese))
[5]
KIM H M，PARK D，RYU D W，et al. Parametric sensitivity analysis of ground uplift above pressurized underground rock caverns[J]. Engineering Geology，2012；135/136：60-65.
[15]
王 娇，林 姗，罗 柳，等. 张拉-剪切破坏模式下边坡三维临界滑面的Cauchy问题[J]. 岩土力学，2022，43(9)：2 634-2 642. (WANG Jiao，LIN Shan，LUO Liu，et al. Cauchy problem of three-dimensional critical slip surface of slope under tension-shear failure mode[J]. Rock and Soil Mechanics，2022，43(9)：2 634-2 642.(in Chinese))
[10]
LITTLEJOHN G S，BRUCE D A. Rock anchors-state of the art，Part 1：design (1)[J].Ground Engineering，1975，8(3)：25-32.
[2]
徐新桥，杨春和，李银平. 国外压气蓄能发电技术及其在湖北应用的可行性研究[J]. 岩石力学与工程学报，2006，25(增2)：3 987-3 992.(XU Xinqiao，YANG Chunhe，LI Yinping. Review on compressed air energy storage abroad and its feasibility application to Hubei province[J]. Chinese Journal of Rock Mechanics and Engineering，2006，25(Supp.2)：3 987-3 992.(in Chinese))
[4]
蔡晓鸿，蔡勇平. 水工压力隧洞结构应力计算[M]. 北京：中国水利水电出版社，2004：40-48.(CAI Xiaohong，CAI Yongping. Stress calculation of hydraulic pressure tunnel[M]. Beijing：China Water and Power Press，2004：40-48.(in Chinese))
[6]
PERAZZELLI P，ANAGNOSTOU G，AMBERG J. Uplift analysis for CAES tunnels[C]// The 14th International Conference of the International Association for Computer Methods and Advances in Geomechanics. Kyoto：CRC Press；2014：1 691-1 696.
[7]
BRANDSHAUG T，CHRISTIANSON M，DAMJANAC B. Technical review of the lined rock caverns (LRC) concept and design methodology：mechanical response of rock mass[R]. Minnesota：Itasca Consulting Group，Inc.，2001.
[12]
张国祥，刘新华，魏 伟. 二维边坡滑动面及稳定性弹塑性有限元分析[J]. 铁道学报，2003，(2)：79-83.(ZHANG Guoxiang，LIU Xinhua，WEI Wei. Analysis of elastic-plastic finite element for two dimension slope slip surface and stability[J]. Journal of the China Railway Society，2003，(2)：79-83.(in Chinese))
[14]
朱以文，吴春秋，蔡元奇. 基于滑移线场理论的边坡滑裂面确定方法[J]. 岩石力学与工程学报，2005，24(15)：2 609-2 616.(ZHU Yiwen，WU Chunqiu，CAI Yuanqi. Determination of slip surface slope based on theory of slip line field[J]. Chinese Journal of Rock Mechanics and Engineering，2005，24(15)：2 609-2 616.(in Chinese))
[8]
TUNSAKUL J，JONGPRADIST P，SOPARAT P，et al. Analysis of fracture propagation in a rock mass surrounding a tunnel under high internal pressure by the element-free Galerkin method[J]. Computers and Geotechnics，2014，55：78-90.
2024, Vol. 43 
