节理倾角及组数对铁路隧道围岩稳定性的影响.pdf

收稿日期:20230210基金项目:国家自然科学基金高铁联合基金项目(U1934213)。作者简介:罗胜利(1980),男,2003 年毕业于西南交通大学土木工程专业,工学学士,高级工程师,E-mail:1351256741 。通信作者:张志强(1968),男,1990 年毕业于重庆大学采矿专业,工学博士,教授,E-mail:214317032 。文章编号:16727479(2023)05004806节理倾角及组数对铁路隧道围岩稳定性的影响罗胜利1 伍容兵1 张志强2(1.中铁二院昆明勘察设计研究院有限责任公司,昆明 650200;2.西南交通大学土木工程学院,成都 610031)摘 要:我国西南地区崇山峻岭、地质构造条件复杂多变,岩体内部节理、片理、层理发育。隧道穿越节理发育围岩时,极易引发围岩大变形、掌子面失稳坍塌、钢架变形扭曲、初支掉块和二次衬砌开裂等工程灾害。为了分析节理对隧道围岩稳定性的影响规律,依托玉磨铁路西双版纳隧道,利用 ABAQUS建立计算分析模型,得到不同节理条件下围岩塑性破坏特征。(1)节理对称分布时,节理屈服、围岩塑性应变呈现出对称分布于拱部、两侧拱肩和仰拱两侧区域的特点。(2)节理倾角较陡时,岩体性质是影响主控因素;节理倾角较平缓时,围岩发生沿节理面的剪切滑移破坏,节理是围岩整体发生塑性破坏的主控因素;当节理倾角为 60或 120时,围岩的塑性应变最大,最大塑性应变为 0.197。(3)当节理倾角为 90时,围岩及节理屈服区域主要沿着节理方向垂向分布,且影响范围深入地层中。(4)2 组节理条件下造成围岩塑性破坏的主要原因是节理面的塑性屈服;当节理倾角组合为 60+90时,围岩的塑性应变最大,最大塑性应变为 0.521。关键词:铁路隧道;节理倾角;节理组数;塑性区;剪切滑移;塑性屈服中图分类号:U239.9;U451.2;TU454 文献标识码:ADOI:10.19630/ki.tdkc.202302100001开放科学(资源服务)标识码(OSID):The Influence of Joint Inclination Angle and Set Number on Stability of Railway Tunnel Surrounding RockLUO Shengli1 WU Rongbing1 ZHANG Zhiqiang2(1.Kunming Design and Research Institute,China Railway Eryuan Engineering Group Co.,Ltd.,Kunming 650200,China;2.School of Civil Engineering,Southwest Jiaotong University,Chengdu 610031,China)Abstract:The southwest of China is a mountainous region with complex and variable geological and tectonic conditions.The rock has developed internal joints,lamellae and laminae.When the tunnel crosses the jointed rock,it is very easy to cause engineering disasters such as large deformation of the surrounding rock,destabilization and collapse of the palm face,deformation of the steel frame,falling blocks of the initial support and cracking of the secondary lining.In order to analyze the influence of joints on the stability of tunnel surrounding rock,a study was carried out based on the Xishuangbanna Tunnel of Yuxi-Mohan Railway.ABAQUS was used to establish a computational analysis model to obtain the plastic damage characteristics of the surrounding rock under different jointing conditions.(1)When the joints are distributed symmetrically,the yield of the joints and the plastic strain of the surrounding rock are symmetrically distributed in the areas of the arch,spandrels on both sides,and both sides of the invert;(2)When the dip angle of the joint is steep,the rock mass property is the main controlling factor.When the joint dip angle is relatively gentle,the surrounding 84铁 道 勘 察2023 年第 5 期rock will be damaged by shear slip along the joint surface.In this case,the joints are the main controlling factors for the plastic failure of the surrounding rock as a whole.When the joint dip angle is 60 or 120,the plastic strain of the surrounding rock is the largest,and the maximum plastic strain is 0.197;(3)When the joint dip angle is 90,the surrounding rock and joint yield area are mainly distributed vertically along the joint direction,and the influence range is deep into the stratum;(4)The main reason for the plastic failure of the surrounding rock under the two groups of joint conditions is the plastic yield of the joint surface.When the joint dip angle combination is 60+90,the plastic strain of the surrounding rock is the largest,and the maximum plastic strain is 0.521.Key words:railway tunnel;joint dip angle;joint set number;plastic zone;shear slip;plastic yield引言节理作为岩体中存在的天然软弱结构面,普遍存在于工程地质中。而这些节理面又将完整的岩体切割得极为破碎,使得围岩的结构完整性大幅降低1-4。越来越多的隧道工程表明,岩体内部节理面的存在,是造成隧道工程围岩变形过大、隧道整体塌方失稳的主要原因之一5-9。目前,隧道围岩失稳分析常用方法主要分为室内试验、现场监测、数值模拟等。在室内试验与现场监测方面,张志强等采用正交试验研究裂隙岩体节理对隧道稳定性的影响10;索超峰等通过模型试验研究多组节理情况下隧道施工扰动后围岩的塌落范围变化情况11;李唱唱等利用水压致裂法和三维水压致裂法对地应力进行现场监测与分析,研究深埋高地应力下引水隧洞节理围岩的稳定性问题12;戚伟等采用理论分析、现场监测多种手段结合,对焦家金矿优势节理组不同组合形成的块体稳定性进行分析13。常用的数值模拟方法主要包括有限单元法和离散单元法。其中有限单元法是基于连续介质理论的,把物体划分为有限个单元,节点之间用数学方程联系起来。马天辉等采用有限元模拟节理岩体中隧洞围岩损伤破坏过程,研究节理岩体中隧洞围岩体的破坏机理,分析岩体中节理倾角对隧洞围岩稳定性的影响规律14;袁铁等通过数值模拟方法研究节理岩体隧道开挖后围岩塑性区分布,得出水平节理使拱顶及仰拱周边围岩的塑性区增大的结论15;王昊等建立导洞遍布节理模型,进行参数敏感性分析16。离散单元法是典型的非连续数值计算方法,已大量应用在地下巷道的稳定性分析中。李军等采用 UDEC 对隧道岩体节理进行模拟,得出单组节理下隧道拱顶变形为多组节理下的 1/4 左右17;郑余朝等采用 3DEC 软件对不同节理情况下的隧道安全系数进行分析,评价各因素对隧道稳定性的影响18;袁彬等利用 UDEC 离散元软件分析了在 090倾角范围内的隧道围岩变形破坏特征19;张斌等基于三维离散元计算方法,模拟分析不同倾角下层状隧道围岩的稳定性20。综上所述,很少有学者针对节理围岩发生塑性破坏的主控因素进行深入研究。基于此,以下采用遍布节理本构模型,分析节理倾角、节理组数对围岩稳定性的影响,并探究围岩发生塑性破坏的主控因素,为以后进一步研究节理发育围岩失稳破坏提供参考。1 工程概况玉磨铁路线路位于云南省,由昆玉铁路的玉溪西站起,经峨山、元江、普洱、景洪止于勐腊边境磨憨口岸,是中老铁路的重要组成部分。玉磨铁路线路全长 507 km,设计速度 160 km/h,是“一带一路”建设的标志性工程之一。西双版纳双线铁路隧道位于野象谷站西双版纳站区间,全长 10 680 m,最大埋深 620 m。隧址区属低中山地貌,地面高程 5801 325 m,相对高差 20750 m。受区域构造影响,测段构造极为发育,发育 3 条断层、1 个背斜,走向均为 NW 向。隧道施工现场实际开挖揭示为:级围岩为 0.4%,级围岩为 29%,级围岩为70.6%。DK354+000DK357+750 段以炭质片岩为主,受软弱节理面影响,呈弹塑性体,岩体强度较低。掌子面岩体节理发育且呈碎裂状,并存在许多分离体,掌子面岩体状况见图 1。图 1 玉磨铁路西双版纳隧道掌子面岩体状况94节理倾角及组数对铁路隧道围岩稳定性的影响:罗胜利 伍容兵 张志强根据隧道已开挖段初支表面监控量测情况,西双版纳隧道 2 号斜井工区小里程方向,DK354+131 DK354+075 段拱顶沉降累计最大为 31.2 cm,最大累计收敛 416 mm,DK354+103 左拱腰最大侵限 为56.3 cm。针对大变形段,为了防止变形进一步扩大,进而发生大规模的围岩垮塌事故,现场采取临时套拱和临时斜撑的方案,对变形未稳定段落进行临时支护。在变形稳定后,为了满足隧道限界要求,需要凿除已侵限的围岩和初支,进行换拱施工。2 计算模型为揭示穿越节理发育地层时由于施工开挖引起的隧道周边围岩变形及破坏特征,计算模型见图 2,根据圣维南原理,隧道距离计算模型边界均大于 3.5 倍洞径,计算模型尺寸为 103.8 m(X)30 m(Y)101.6 m(Z)。考虑计算结果的收敛性,隧道计算埋深取100 m,隧道拱顶至模型上表面的高度为 50 m,其余覆土厚度采用施加等效表面荷载形式表示。计算模型共离散成 65 580 个六面体网格,除模型上表面为荷载边界外,其余表面均采用位移边界,左右表面限制 X 方向自由度,前后表面限制 Y 方向自由度,底面限制 X、Y和 Z 方向自由度。为了突显断层破碎带节理围岩在隧道开挖影响下的变形和破坏特征,计算过程采用全断面开挖,且不考虑支护结构,隧道处于毛洞状态,开挖进尺取 1.0 m。图 2 三维整体计算模型根据隧道节理发育区域的围岩资料以及 TB 100032016铁路隧道设计规范,确定计算模型中围岩及节理的参数,具体参数见表 1。在节理发育的断层带中,岩体在节理的“切割”作用下,岩体呈现出破碎、碎裂的状态,围岩越破碎,对隧道施工越不利。对于隧道施工而言,围岩的节理倾角和节理数量不同,隧道开挖后围岩的变形特征和破坏形态将会有较大的差异。现场显示平均节理间距0.05 m,节理之间的间距相比起模型的结构尺寸足够小。遍布节理模型不考虑节理面厚度,且节理贯通整个遍布节理模型。单组节理和两组节理工况见图 3。节理倾角、节理组数设置情况见表 2。表 1 围岩及节理力学参数项目弹性模量/GPa泊松比内摩擦角/()黏聚力/MPa容重/(kNm-3)围岩参数1.00.39260.1119.4节理参数230.02图 3 节理工况示意表 2 模型工况设置节理数量节理组合倾角1 组节理0、20、40、60、80、90、100、120、140、160、1802 组节理0+60、10+60、20+60、30+60、40+60、50+60、70+60、80+60、90+60、100+60、110+60、120+60、150+603 节理倾角对隧道围岩稳定性的影响针对单组节理围岩,开挖完成后,截取计算模型Y=15 m 处截面的节理屈服云图和塑性应变云图,见图 4图 9。图 4 0节理各工况下围塑性破坏特征图 5 20节理各工况下围塑性破坏特征05铁 道 勘 察2023 年第 5 期图 6 40节理各工况下围塑性破坏特征图 7 60节理各工况下围塑性破坏特征图 8 80节理各工况下围塑性破坏特征图 9 90节理各工况下围塑性破坏特征由图 4 可知,当节理倾角为 0,由于模型的节理、模型的几何对称性,节理屈服、围岩塑性应变呈现出左右对称、上下对称的特点,主要分布于拱部两侧拱肩和仰拱两侧区域,此区域受施工扰动明显,在施工影响、重力以及应力重分布等综合作用下,发生了塑性屈服,施工过程中不及时加以控制,易发生较大变形甚至是坍塌事故。由图 5 图 8 可知,当倾角由 0逐渐增大至 90时,随着倾角增大,受隧道开挖扰动围岩及节理塑性屈服分布特征发生改变,不再是对称分布。相较于水平节理,围岩及节理塑性屈服分布形态和发展趋势发生偏转,倾角较小(20、40)时,塑性区主要沿着垂直于节理的方向分布,此时岩体性质是影响塑性区形成和分布主控因素;随着倾角进一步增加(60、80),塑性分布区及发展域包括两部分,一部分是沿着节理方向;另一部分是垂直于节理方向。其中,沿着节理方向的塑性区范围更大,由此说明在此倾角情况下,节理是围岩整体发生塑性破坏的主控因素。由图9 可知,当节理倾角为 90时,由于模型的节理、模型的几何对称性,围岩及节理屈服区域主要分布在隧道两侧,沿着节理方向垂向分布且影响范围深入地层中。节理是影响围岩发生塑性破坏的主控因素。由此可知,对于垂直节理岩体,施工过程中一旦节理屈服且发生显著的剪切滑移,则可能导致拱顶岩体整体垮塌,造成施工灾难,所以隧道穿越90节理围岩施工较为危险。当倾角由 90逐渐增大至 180时,节理倾角越接近 90,围岩倾角越陡,节理倾角越接近 180,围岩倾角越平缓。根据反对称性,节理倾角较陡时,岩体性质是影响塑性区形成和分布主控因素;节理倾角较平缓,围岩发生沿节理面的剪切滑移破坏节理是围岩整体发生塑性破坏的主控因素,与节理倾角由 090时的规律相同。提取单组节理各倾角工况的塑性应变最大值,并绘制曲线,见图 10。由图 10 可知,当节理倾角为 60或 120时,围岩的塑性应变最大,最大塑性应变为0.197,说明节理倾角为 60或 120时,在开挖扰动下,隧道周边围岩发生塑性破坏的趋势最明显,对隧道施工最为不利。该结论与通过隧道洞周围岩变形量分析结论一致。图 10 单组节理各倾角组合工况的塑性应变4 节理组数对隧道围岩稳定性的影响计算完成后,截取 2 组节理工况下计算模型 Y=15 m 处截面的围岩塑性应变及节理屈服云图,见图11图 14。由图 11图 14 可知,2 组节理将围岩进一步“切割”,围岩变得较为破碎,根据节理屈服云图,节理屈服区域主要分布在洞周两侧拱肩至拱脚区域,节理屈服区域不再沿着节理或垂直节理方向深入围岩内部,由此说明,2 组节理条件下节理面塑性屈服与节理方向无明显关系。岩体发生塑性应变的部位与节理屈服部位重合,且塑性应变区域面积和深度均小于节理塑15节理倾角及组数对铁路隧道围岩稳定性的影响:罗胜利 伍容兵 张志强图 11 0+60节理各工况下围塑性破坏特征图 12 30+60节理各工况下围塑性破坏特征图 13 90+60节理各工况下围塑性破坏特征图 14 120+60节理各工况下围塑性破坏特征性屈服。因此,对于存在 2 组节理的地层,围岩整体破碎,开挖后洞周塑性区主要分布在隧道两侧拱肩至墙脚区域,造成围岩塑性破坏的主要原因是节理面的塑性屈服。图 15 2 组节理各倾角组合工况的塑性应变提取 2 组节理各倾角工况的塑性应变最大值,并绘制曲线,见图 15。由图 15 可知,当节理倾角组合为60+90时,围岩的塑性应变最大,最大塑性应变为0.521,说明节理倾角为 60+90时,在开挖扰动下,隧道周边围岩发生塑性破坏的趋势最明显,对隧道施工最为不利。该结论与通过隧道洞周围岩变形量分析结论一致。5 结论采用有限元软件,研究节理倾角和节理组数对隧道围岩塑性破坏特征的影响,得出以下主要结论。(1)节理倾角水平时,节理屈服、围岩塑性应变呈现出左右对称、上下对称的特点,主要分布于拱部两侧拱肩和仰拱两侧区域。(2)节理倾角较陡时,岩体性质是影响塑性区形成和分布主控因素;节理倾角较平缓,围岩发生沿节理面的剪切滑移破坏节理是围岩整体发生塑性破坏的主控因素。当节理倾角为 60或 120时,围岩的塑性应变最大,最大塑性应变为 0.197。(3)当节理倾角为 90时,围岩及节理屈服区域主要分布在隧道两侧,沿着节理方向垂向分布,且影响范围深入地层中,节理是影响围岩发生塑性破坏的主控因素。(4)2 组节理条件下开挖后,洞周塑性区主要分布在隧道两侧拱肩至墙脚区域,造成围岩塑性破坏的主要原因是节理面的塑性屈服。当节理倾角组合为60+90 时,围岩的塑性应变最大,最大塑性应变为 0.521。参考文献1 李婷,姜谙男,张峰瑞,等.节理岩体对盾构开挖稳定性和地层损失率的影响J.工业建筑,2022,52(3):29-37.LI Ting,JIANG Annan,ZHANG Fengrui,et al.Effects of Jointed Rock Masses on Stability of Shield Tunnels under Tunnelling and Ground Loss Rates J.Industrial Construction,2022,52(3):29-37.2 吴斐,温森,孔庆梅.节理岩体力学参数对双护盾 TBM 卡机的影响J.北京交通大学学报,2022,46(1):105-114.WU Fei,WEN Sen,KONG Qingmei.Influence of Mechanical Parameters for Jointed Rock Masses on Double-shield TBM JammingJ.Journal of Beijing Jiaotong University,2022,46(1):105-114.3 杨卫领,郭佳奇,黄猛,等.爆破与开挖卸荷作用下岩溶隧道突水灾变过程研究J.公路,2022,67(10):432-438.YANG Weiling,GuoJia qi,HUANG Meng,et al.Study on Sudden Flood Process of Karst Tunnel under the Action of Blasting and Excavation UnloadingJ.Highway,2022,67(10):432-438.4 周君,周洋,陈林,等.节理倾角对软岩隧道掌子面变形及失稳模式的影响J.路基工程,2022(6):196-200.ZHOU Jun,ZHOU Yang,CHEN Lin,et al.Influence of Joint Dip Angle on Deformation and Instability Mode of Soft Rock Tunnel Face25铁 道 勘 察2023 年第 5 期J.Subgrade Engineering,2022(6):196-200.5 张子新,王帅峰,黄昕.节理岩体中 TBM 掘进施工围岩稳定性及可视化分析J.现代隧道技术,2018,55(2):36-43,64.ZHANG Zixin,WANG Shuaifeng,HUANG Xin.Analysis and 3D Visualization of Surrounding Rock Stability during TBM Tunnelling in Jointed Rock Masses J.Modern Tunnelling Technology,2018,55(2):36-43,64.6 高阳,孙浩凯,刘德军,等.强降雨影响下破碎复理岩地层隧道洞口段失稳机理J.中南大学学报(自然科学版),2019,50(9):2295-2303.GAO Yang,SUN Haokai,LIU Dejun,et al.Collapse Mechanism of Tunnel Portal Sectionin Broken Flysch under Influence of Heavy Rainfall J.Journal of Central South University(Science and Technology),2019,50(9):2295-2303.7 潘文韬,何川,吴枋胤,等.不同大变形等级下层理角度对层状软岩隧道的影响J.土木与环境工程学报(中英文),2023,45(5):94-105.PAN Wentao,HE Chuan,WU Fangyin,et al.Effect of bedding angle of layered soft rock tunnels with different large deformation gradesJ.Journal of Civil and Environmental Engineering,2023,45(5):94-105.8 王亚琼,杨强,潘红伟,等.基于 3DEC 模拟的高地应力水平层状隧道围岩变形破坏特征分析J.现代隧道技术,2022,59(4):127-136,146.WANG Yaqiong,YANG Qiang,PAN Hongwei,et al.Analysis on Deformation and Failure Characteristics of Surrounding Rock in Horizontal Layered Tunnels under High Crustal Stress Based on 3DEC SimulationJ.Modern Tunnelling Technology,2022,59(4):127-136,146.9 王心泉,王智猛,牛犇,等.8 度地震烈度区新民隧道出口处边坡的稳定性J.山东大学学报(工学版),2023,53(3):23-30,40.WANG Xinquan,WANG Zhimeng,NIU Ben,et al.Slope stability at the exit of Xinmin tunnel in 8 degree seismic intensity area J.Journal of Shandong University(Engineering Science),2023,53(3):23-30,40.10 张志强,何本国,李永珑,等.节理岩体隧道稳定性影响因素正交试验及锚杆作用要点研究J.武汉理工大学学报,2013,35(2):90-96.ZHANG Zhiqiang,HE Benguo,LI Yonglong,et al.Study on Orthogonal Test of Influencing Factors of Tunnel in Joined Rock Mass and Mechanical Effect of Bolt Support J.Journal of Wuhan University of Technology,2013,35(2):90-96.11 索超峰,石益东,李军.节理特征对破碎围岩稳定性影响的模型试验J.公路交通科技,2013,30(4):82-87.SUO Chaofeng,SHI Yidong,LI Jun.Model Test of Effect of Joint Characteristics on Stability of Broken Rock Mass J.Journal of Highway and Transportation Research and Development,2013,30(4):82-87.12 李唱唱,侍克斌,姜海波.深埋高地应力引水隧洞节理围岩稳定性研究J.水资源与水工程学报,2020,31(2):219-224.LI Changchang,SHI Kebin,JIANG Haibo.Study on stability of jointed surrounding rock of a deep-lying diversion tunnel with high ground stress J.Journal of Water Resources and Water Engineering,2020,31(2):219-224.13 戚伟,付建新,李腾.考虑节理劣化效应的深部围岩应力演化规律及敏感性分析J.采矿与安全工程学报,2020,37(2):327-337.QI Wei,FU Jianxin,LI Teng.Stress Evolution Law and Sensitivity Analysis of Surrounding Rock in Deep Mine Considering Deterioration Effect Caused by JointJ.Journal of Mining&Safety Engineering,2020,37(2):327-337.14 马天辉,张文东,徐涛.节理岩体中隧洞围岩的损伤破坏机理J.东北大学学报(自然科学版),2013,34(10):1485-1489.MA Tianhui,ZHANG Wendong,XU Tao.Damage and Failure Mechanism of Tunnels in Jointed Rock Mass J.Journal of Northeastern University(Natural Science),2013,34(10):1485-1489.15 袁铁,李畅,申志军,等.节理对泥岩隧道稳定性影响分析J.公路,2020,65(2):288-294.YUAN Tie,LI Chang,SHEN Zhijun,et al.Effect of Joints on Stability of Mudstone TunnelJ.Highway,2020,65(2):288-294.16 王昊,姜谙男,郑帅,等.基于遍布节理模型的围岩稳定性分析J.工业建筑,2022,52(9):193-197.WANG Hao,JIANG Annan,ZHENG Shuai,et al.Stability Analysis of Surrounding Rock Based on the Ubiquitous Joint Model J.Industrial Construction,2022,52(9):193-197.17 李军,吴惠卿,王彦杰,等.节理特征和围岩支护对隧道变形的影响分析J.施工技术,2013,42(3):82-86.LI Jun,WU Huiqing,WANG Yanjie,et al.Analysis of the Influence of Joints Characteristics and Surrounding Rock Supporting on Tunnel DeformationJ.Construction Technology,2013,42(3):82-86.18 郑余朝,张文胜,孙克国,等.基于离散元强度折减法的大跨隧道稳定性研究J.现代隧道技术,2020,57(1):18-25.ZHENG Yuchao,ZHANG Wensheng,SUN Keguo,et al.Study on the Stability of Large-span Tunnel by Distinct Element Based Strength Reduction Method J.Modern Tunnelling Technology,2020,57(1):18-25.19 袁彬,徐凡献,廖欢,等.控制型结构面几何参数对隧道围岩稳定性的控制效应:以郑万高铁向家湾隧道为例J.长江科学院院报,2020,37(4):115-121.YUAN Bin,XU Fanxian,LIAO Huan,et al.Control Effect of Controlling Structural Planes Geometric Parameters on Stability of Tunnels Surrounding RockCase Study on Xiangjiawan Tunnel of Zhengzhou-Wanzhou High-speed Railway J.Journal of Yangtze River Scientific Research Institute,2020,37(4):115-121.20 张斌,李英杰,穆鹏华,等.不同节理倾角下隧道围岩变形特征研究J.地下水,2022,44(6):174-175.ZHANG Bin,LI Yingjie,MU Penghua,et al.Study on Deformation Characteristics of Tunnel Surrounding Rock under Different Joint Dip AnglesJ.Ground Water,2022,44(6):174-175.35节理倾角及组数对铁路隧道围岩稳定性的影响:罗胜利 伍容兵 张志强