㊀第４９卷第５期
煤炭科学技术
Ｖｏｌ ４９㊀Ｎｏ ５㊀
㊀２０２１年
５月
ＣｏａｌＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ
㊀Ｍａｙ
２０２１㊀
移动扫码阅读
倪向龙，龙建辉，张玲玲．黄土地区高速公路隧道洞口滑坡变形机理分析［Ｊ］．煤炭科学技术，２０２１，４９（５）：２３４－２４１ ｄｏｉ：１０ １３１９９／ｊ ｃｎｋｉ ｃｓｔ ２０２１ ０５ ０３０ＮＩＸｉａｎｇｌｏｎｇ，ＬＯＮＧＪｉａｎｈｕｉ，ＺＨＡＮＧＬｉｎｇｌｉｎｇ．Ａｎａｌａｙｓｉｓｏｎｄｅｆｏｒｍａｔｉｏｎｍｅｃｈａｎｉｓｍｏｆｌａｎｄｓｌｉｄｅａｔｔｕｎｎｅｌｅｎｔｒａｎｃｅｏｆｅｘｐｒｅｓｓｗａｙｉｎｌｏｅｓｓａｒｅａ［Ｊ］．ＣｏａｌＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，２０２１，４９（５）：２３４－２４１ ｄｏｉ：１０ １３１９９／
ｊ ｃｎｋｉ ｃｓｔ ２０２１ ０５ ０３０
黄土地区高速公路隧道洞口滑坡变形机理分析
倪向龙１，２，龙建辉１，张玲玲１
（１．太原理工大学矿业工程学院，山西太原㊀０３００２４；２．中国建筑西南勘察设计研究院有限公司，四川成都㊀６１００５１）
摘㊀要：黄土地区修建高速公路隧道时，时常会遇到由于隧道洞口的施工造成上覆黄土发生蠕滑，使
隧道因滑坡滑动而产生变形㊁开裂等病害，造成隧道无法正常运行㊂基于上述问题，通过收集黄土地区已发生的隧道洞口滑坡资料，选取６组典型案例分析该类型滑坡的诱发因素与形成条件，结合地质结构模型对隧道洞口滑坡的变形发展过程和运动特征进行了理论分析，并利用数值模拟方法研究隧道洞口开挖对上覆黄土的影响规律，确定隧道在滑坡体中的受力模式，结合典型工程实例对该规律进行验证㊂研究结果表明：黄土地区隧道洞口滑坡多发生于黄土塬梁地区，坡脚常有冲沟存在，滑坡体主要为第四系风积黄土与马兰黄土，其中人工扰动与降雨入渗为滑坡主要诱发因素㊂诱发机理为隧道施工坡脚出现邻空面，坡体内部开挖，造成坡体应力重分布，同时坡体在降雨条件下降低了岩土体的力学性质㊂根据诱发原因，滑坡的破坏机理与形成过程分为４个阶段：隧道开挖应力调整㊁蠕滑剪切与拉裂㊁变形加速阶段和坡体失稳㊂有限元数值模拟案例也很好地反映了滑坡体内的变形场㊁应力场及破坏区分布，且隧道与滑坡的空间位置关系控制滑坡的性质和规模，研究结果可为后续类似的隧道洞口滑坡提供实践经验与理论指导㊂
关键词：黄土地区；隧道洞口；滑坡；滑移特征；变形规律
中图分类号：ＴＤ１６３㊀㊀㊀文献标志码：Ａ㊀㊀㊀文章编号：０２５３－２３３６（２０２１）０５－０２３４－０８
Ａｎａｌａｙｓｉｓｏｎｄｅｆｏｒｍａｔｉｏｎｍｅｃｈａｎｉｓｍｏｆｌａｎｄｓｌｉｄｅａｔｔｕｎｎｅｌｅｎｔｒａｎｃｅｏｆ
ｅｘｐｒｅｓｓｗａｙｉｎｌｏｅｓｓａｒｅａ
ＮＩＸｉａｎｇｌｏｎｇ１，２，ＬＯＮＧＪｉａｎｈｕｉ１，ＺＨＡＮＧＬｉｎｇｌｉｎｇ１
（１．ＣｏｌｌｅｇｅｏｆＭｉｎｉｎｇＥｎｇｉｎｅｅｒｉｎｇ，ＴａｉｙｕａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｔａｉｙｕａｎ㊀０３００２４，Ｃｈｉｎａ；２．ＣｈｉｎａＳｏｕｔｈｗｅｓｔＧｅｏｔｅｃｈｎｉｃａｌ
Ｉｖｅｓｔｉｇａｔｉｏｎ＆ＤｅｓｉｇｎＩｎｓｔｉｔｕｔｅＣｏ．，Ｌｔｄ．，Ｃｈｅｎｇｄｕ㊀６１００５１，Ｃｈｉｎａ）
收稿日期：２０２０－１２－１２；责任编辑：曾康生
基金项目：山西省自然科学基金资助项目（２０１８０１Ｄ１２１０３３，２０１７０１Ｄ１２１０１５）
作者简介：倪向龙（１９９４ ），男，内蒙古呼和浩特人，硕士研究生㊂Ｅ－ｍａｉｌ：２５７７５０５６０７＠ｑｑ．ｃｏｍ
Ａｂｓｔｒａｃｔ：Ｗｈｅｎｃｏｎｓｔｒｕｃｔｉｎｇｅｘｐｒｅｓｓｗａｙｔｕｎｎｅｌｓｉｎ
ｌｏｅｓｓａｒｅａ，ｃｒｅｅｐｏｆｏｖｅｒｌｙｉｎｇｌｏｅｓｓｃａｕｓｅｄｂｙｔｕｎｎｅｌｏｐｅｎｉｎｇｃｏｎｓｔｒｕｃｔｉｏｎｏｆｔｅｎｏｃ⁃ｃｕｒｒｅｄ，ｗｈｉｃｈｃａｕｓｅｄｄｅｆｏｒｍａｔｉｏｎａｎｄｃｒａｃｋｉｎｇｏｆｔｕｎｎｅｌｓｄｕｅｔｏｌａｎｄｓｌｉｄｅｓｌｉｄｉｎｇ，ａｎｄｍａｄｅｔｈｅｔｕｎｎｅｌｓｕｎａｂｌｅｔｏｏｐｅｒａｔｅｎｏｒｍａｌｌｙ．Ｂａｓｅｄｏｎｔｈｅａｂｏｖｅｐｒｏｂｌｅｍｓ，ｃｏｌｌｅｃｔｅｄｔｈｅｄａｔａｏｆｔｕｎｎｅｌｅｎｔｒａｎｃｅｌａｎｄｓｌｉｄｅｉｎＬｏｅｓｓａｒｅａ，ｓｅｌｅｃｔｅｄｓｉｘｔｙｐｉｃａｌｃａｓｅｓｔｏａｎａｌｙｚｅｔｈｅｉｎｄｕ⁃ｃｉｎｇｆａｃｔｏｒｓａｎｄｆｏｒｍｉｎｇｃｏｎｄｉｔｉｏｎｓｏｆｔｈｉｓｔｙｐｅｏｆｌａｎｄｓｌｉｄｅ，ｃｏｍｂｉｎｅｄｗｉｔｈｇｅｏｌｏｇｉｃａｌｓｔｒｕｃｔｕｒｅｍｏｄｅｌｔｏａｎａｌｙｚｅｔｈｅｄｅｆｏｒｍａｔｉｏｎｄｅｖｅｌｏｐ⁃ｍｅｎｔｐｒｏｃｅｓｓａｎｄｍｏｖｅｍｅｎｔｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｔｕｎｎｅｌｅｎｔｒａｎｃｅｌａｎｄｓｌｉｄｅ．Ｔｈｅｉｎｆｌｕｅｎｃｅｌａｗｏｆｔｕｎｎｅｌｏｐｅｎｉｎｇｅｘｃａｖａｔｉｏｎｏｎｏｖｅｒｌｙｉｎｇｌｏｅｓｓｗａｓｓｔｕｄｉｅｄｂｙｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｍｅｔｈｏｄ，ａｎｄｔｈｅｓｔｒｅｓｓｍｏｄｅｏｆｔｕｎｎｅｌｉｎｌａｎｄｓｌｉｄｅｗａｓｄｅｔｅｒｍｉｎｅｄ，ｗｈｉｃｈｗａｓｖｅｒｉｆｉｅｄｂｙｔｙｐｉｃａｌｅｎ⁃ｇｉｎｅｅｒｉｎｇｅｘａｍｐｌｅｓ．Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｔｈａｔｔｈｅｌａｎｄｓｌｉｄｅｓａｔｔｈｅｔｕｎｎｅｌｅｎｔｒａｎｃｅｉｎｔｈｅｌｏｅｓｓａｒｅａｍｏｓｔｌｙｏｃｃｕｒｒｅｄｉ
ｎｔｈｅＬｏｅｓｓＰｌａｔｅａｕＢｅａｍａｒｅａａｎｄｔｈｅｇｕｌｌｉｅｓｏｆｔｅｎｅｘｉｓｔｅｄａｔｔｈｅｆｏｏｔｏｆｔｈｅｓｌｏｐｅ．ＬａｎｄｓｌｉｄｅｂｏｄｙｗａｓｍａｉｎｌｙＱｕａｔｅｒｎａｒｙＡｅｏｌｉａｎＬｏｅｓｓａｎｄＭａｌａｎｌｏｅｓｓ．Ａｒｔｉｆｉｃｉａｌｄｉｓｔｕｒｂａｎｃｅａｎｄｒａｉｎｆａｌｌｉｎｆｉｌｔｒａｔｉｏｎｗｅｒｅｔｈｅｍａｉｎｉｎｄｕｃｉｎｇｆａｃｔｏｒｓｏｆｌａｎｄｓｌｉｄｅ．Ｔｈｅｉｎｄｕｃｉｎｇｍｅｃｈａｎｉｓｍｗａｓｔｈａｔｔｈｅｒｅｗａｓａｎｅｍｐｔｙｓｕｒ⁃ｆａｃｅａｔｔｈｅｆｏｏｔｏｆｔｈｅｓｌｏｐｅｉｎｔｕｎｎｅｌｃｏｎｓｔｒｕｃｔｉｏｎ，ａｎｄｔｈｅｅｘｃａｖａｔｉｏｎｉｎｓｉｄｅｔｈｅｓｌｏｐｅｃａｕｓｅｄｔｈｅｒｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆｓｔｒｅｓｓｉｎｔｈｅｓｌｏｐｅｂｏｄｙ．
Ａｔｔｈｅｓａｍｅｔｉｍｅ，ｔｈｅｍｅｃｈａｎｉｃａｌｐｒｏｐｅｒｔｉｅｓｏｆｒｏｃｋａｎｄｓｏｉｌｍａｓｓｗｅｒｅｒｅｄｕｃｅｄｕｎｄｅｒｔｈｅｃｏｎｄｉｔｉｏｎｏｆｒａｉｎｆａｌｌ．Ａｃｃｏｒｄｉｎｇｔｏｔｈｅｉｎｄｕｃｅｄｃａｕｓｅｓ，ｔｈｅｆａｉｌｕｒｅｍｅｃｈａｎｉｓｍａｎｄｆｏｒｍａｔｉｏｎｐｒｏｃｅｓｓｏｆｌａｎｄｓｌｉｄｅｃａｎｂｅｄｉｖｉｄｅｄｉｎｔｏｆｏｕｒｓｔａｇｅｓ：ｓｔｒｅｓｓａｄｊｕｓｔｍｅｎｔｓｔａｇｅｏｆｔｕｎｎｅｌｅｘｃａｖａ⁃
ｔｉｏｎ，ｃｒｅｅｐｓｈｅａｒａｎｄｔｅｎｓｉｏｎｃｒａｃｋｓｔａｇｅ，ｄｅｆｏｒｍａｔ
ｉｏｎａｃｃｅｌｅｒａｔｉｏｎｓｔａｇｅａｎｄｓｌｏｐｅｉｎｓｔａｂｉｌｉｔｙｓｔａｇｅ．Ｔｈｅｃａｓｅｏｆｆｉｎｉｔｅｅｌｅｍｅｎｔｎｕｍｅｒｉｃａｌ
４
３２
倪向龙等：黄土地区高速公路隧道洞口滑坡变形机理分析２０２１年第５期ｓｉｍｕｌａｔｉｏｎｗｅｌｌｒｅｆｌｅｃｔｓｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄ，ｓｔｒｅｓｓｆｉｅｌｄａｎｄｆａｉｌｕｒｅａｒｅａｉｎｌａｎｄｓｌｉｄｅｂｏｄｙ，ａｎｄｔｈｅｓｐａｔｉａｌｒｅｌａｔｉｏｎｓｈｉｐｂｅ⁃ｔｗｅｅｎｔｕｎｎｅｌａｎｄｌａｎｄｓｌｉｄｅｃｏｎｔｒｏｌｓｔｈｅｎａｔｕｒｅａｎｄｓｃａｌｅｏｆｌａｎｄｓｌｉｄｅ．Ｔｈｅｒｅｓｅａｒｃｈｒｅｓｕｌｔｓｃａｎｐｒｏｖｉｄｅｐｒａｃｔｉｃａｌｅｘｐｅｒｉｅｎｃｅａｎｄｔｈｅｏｒｅｔｉｃａｌｇｕｉｄａｎｃｅｆｏｒｓｉｍｉｌａｒｔｕｎｎｅｌｅｎｔｒａｎｃｅｌａｎｄｓｌｉｄｅ．
Ｋｅｙｗｏｒｄｓ：ｌｏｅｓｓａｒｅａ；ｔｕｎｎｅｌｅｎｔｒａｎｃｅ；ｌａｎｄｓｌｉｄｅ；ｍｏｖｅｍｅｎｔｃｈａｒａｃｔｅｒｉｓｔｉｃｓ；ｄｅｆｏｒｍａｔｉｏｎｌａｗｓ
０㊀引㊀㊀言
近年来，黄土地区随着城市建设的逐步加深，高速公路的建设也在向山区不断拓展㊂在山区修建高速公路时常需要在坡体内开挖隧道，而在隧道施工过程中容易发生滑坡等地质灾害，隧道因坡体滑动而产生变形㊁开裂等病害，严重威胁到工程的建设及运营安全［１－３］㊂
针对隧道开挖诱发滑坡的变形机理与治理措施，国内外学者进行了大量的研究㊂吴红刚等［４］根据隧道与滑动面的空间位置关系，将隧道－滑坡体系划分为平行㊁正交和斜交３大体系类型，并通过工程实例验证了３种模式的受力机制㊂王国欣等［５］通过利用数值模拟与实际的监测数据分析了某隧道洞口滑坡的变形机制㊂张伟等［６］在分析隧道洞口滑坡的形成条件的基础上，使用数值模拟的方法探讨了隧道上方围岩的变形分布特征，得到了滑面位置㊂陈羽等［７］采用不平衡传递系数法对隧道洞口滑坡稳定性进行分析，并提出了相应的治理措施㊂上述文献对隧道洞口滑坡的研究大部分都是基于岩质边坡与碎石坡积物堆积的边坡［８－１０］，而对黄土地区的隧道洞口滑坡研究相对较少㊂黄土边坡由于黄土特殊的结构性㊁水敏性等，决定了黄土滑坡的孕灾机制㊁成灾模式㊁成因机理㊁演化过程㊁链生规律等明显有别于岩质边坡［１１］㊂因此黄土地区隧道洞口滑坡相比较岩质边坡的滑移拉裂，更多的表现为蠕滑拉裂，笔者通过总结黄土地区已发生的隧道洞口滑坡的基本特征与形成条件，分析滑坡的形成机理，并通过数值模拟具体分析某一案例滑坡的发育规律进行验证，为后续发生的类似滑坡的防治提供理论支持㊂
１㊀黄土地区隧道滑坡形成机理
１．１㊀隧道洞口滑坡形成条件
对黄土地区隧道洞口发生的滑坡进行数据收集分析其滑坡的形成原因与形成条件（表１），其中主要以工程地质条件㊁水文地质条件与人类工程活动进行分析㊂
１）工程地质条件㊂黄土地区隧道洞口滑坡多发生于黄土塬梁地区，滑坡体主要为第四系风积黄
表１㊀典型案例形成原因分析
Ｔａｂｌｅ１㊀Ａｎａｌｙｓｉｓｏｆｃａｕｓｅｓｏｆｔｙｐｉｃａｌｃａｓｅｓ工程实例形成原因
洪亮营隧道洞
口滑坡［１２］
洪亮营斜坡是由马兰黄土和离石黄土组成的，
由于隧道开挖扰动黄土内部应力分布，在降雨
的影响下，上部黄土饱和度增大，下部形成隔
水层，造成该边坡发生滑动
龙王庙隧道洞
口滑坡
岢临高速沿线为黄土区，土质疏松，在降雨持
续时间较长的条件下土壤饱和易于产流，导致
龙王庙１号隧道岢岚端洞门端站墙在洪水与
土体的推动下所破坏
青银高速雅沟
隧道洞口滑坡
发生崩塌的山体上方土质十分松软，部分土体
悬空，在降雨作用下使上覆土体含水量增加，
自重增加，同时在土体下部为临空面，扰动破
坏原生地质环境平衡，引起斜坡应力重分布，
集中与坡脚部位，造成斜坡滑动破坏
黄延高速隧道
洞口滑坡［１３］
马岔沟隧道施工队在滑坡体顶部安置有简易
储水罐供下方生活用水，滑坡发生前４ｄ出现
供水异常，可能是引水管渗漏，造成滑坡部位土体
饱和，重度增加，力学强度降低，导致发生滑坡
卧牛山隧道洞
口滑坡［１４］
该滑坡为一古滑坡，由于坡体结构为第四系堆
积物与基岩的组合，在裂隙水与地下水的影响
下会对坡体稳定性产生影响
杜家坪隧道洞
口滑坡［１５］
该隧道洞口边坡地层组合为Ｑ３＋Ｑ２＋红黏土，
上覆土体结构疏松，下伏土体为良好的隔水
层，易软化，在降雨等条件下，容易产生滑坡
土与马兰黄土，土体松散，有利于地表水入渗，下伏地层为离石黄土或基岩接触层，是较好的隔水层㊂这种透水性呈上强下弱的二元地层结构易使水在离石黄土和基岩接触层这类透水性差的岩土层顶面阻滞，形成上层滞水，对滑带土软化后降低其强度，有利滑坡的形成㊂同时隧道洞口滑坡发生其围岩类型多为破碎松散岩石或软弱的塑性岩类与碎裂结构岩体及半坚硬的薄层状结构岩体［１６］㊂
２）水文地质条件㊂黄土地区隧道洞口滑坡发生于沟谷斜坡过渡带，坡脚常有冲沟存在，汇集周边山体雨水，为季节性流水，水流量小，强降雨水量增大，旱季干枯㊂地表水是诱发滑坡形成的外在因素，滑坡体后壁在地形㊁地貌上呈明显的陡峭状，加之滑体土质结构松散，有利于地表水入渗，故在雨季特别是长期降雨时，使滑坡体湿度㊁重度加大，而滑床通常为较好的隔水层，致使滑坡下滑力增大，滑带土的
５３２
２０２１年第５期
煤炭科学技术
第４９卷
抗剪强度下降，对滑坡的形成创造了有利的条件［１７］㊂
３）人类工程活动㊂隧道洞口开挖，多采取机械
与爆破的方法进行开挖，其中隧道爆破中所产生的冲击力对土体结构破坏较严重，降低了土体的抗剪强度，对滑坡的形成创造了有利条件㊂隧道洞口边坡岩体的完整性较差，由于隧道开挖产生新的变形
空间，使得岩体内径向方向由原来的紧密压缩状态向松弛状态过渡，在洞室周围形成松动圈，围岩内的应力状态也将因松动圈的应力被释放而重新调整，形成一定的应力分带，其中产生位移的范围大约为地下洞室横断面中最大尺寸的３ ５倍，且隧道洞口距离上覆黄土距离较近，在应力分带的影响下，黄土坡体产生相应的应力调整，为滑坡的发生提供了动力㊂１．２㊀隧道洞口滑坡形成机理
隧道洞口滑坡的形成主要是坡脚隧道开挖，造成上覆黄土发生蠕滑，由于黄土滞后性和长期蠕变效应，使得上方坡体受控于坡脚隧道开挖的阶段㊂因此根据已发生隧道洞口滑坡的变形发展过程和运动特征，隧道洞口滑坡形成机理可分为４个阶段㊂１．２．１㊀隧道开挖应力调整阶段
岩土体在天然条件下均处于一定的初始应力状态㊂隧道开挖前需要先对坡脚切削出便于施工的掌子面，而坡脚开挖引起应力重分布与应力集中等效应（图１ａ），斜坡周围主应力迹线发生明显的偏转，
在临空面附近造成应力集中带，坡脚附近最大主应力显著增高，形成一最大剪应力增高带，通常是斜坡中最容易发生变形的部位㊂同时隧道洞口初期开挖在原有的基础上更加破坏了岩土体原有的应力平衡状态，使得开挖部分以上岩土体失去支撑而发生自身重力转移，斜坡体上部出现张拉应力，坡体出现蠕滑的迹象［１８－１９］（图１ｂ）㊂
bug free１．２．２㊀蠕滑剪切与拉裂阶段
隧道洞口继续开挖，并对洞口进行初期支护，但是开挖的同时坡体内部应力不断调整，其放大效应影响上覆黄土的蠕滑效果，其中卸荷作用使得坡体内裂缝逐渐加深，形成应力集中部位，产生蠕滑剪切带㊂坡顶在张拉应力作用下，出现拉张裂缝（图１ｃ）㊂１．２．３㊀变形加速阶段
坡体后缘拉张裂缝不断向下扩张降低岩土体本身的力学性质，为雨水进入坡体内部提供了入渗通道㊂在降雨的条件下，雨水进入坡体软化滑带，降低滑带的剪切强度，进而加速土体的破坏［２０－２１］（图１ｄ）㊂１．２．４㊀坡体失稳阶段
随着蠕滑变形的加剧，坡体应力不断调整，变形逐渐向前缘积累，坡体前缘开始显现出由后向前的鼓胀㊁隆起等现象；同时在外力及降雨入渗作用下，后缘裂隙逐渐加深加宽，最终拉张裂缝与蠕变区贯通，滑带形成，坡体发生整体滑动（图
１ｅ）㊂
１ 张应力集中区；２ 剪应力集中区；３ 隧道开挖影响范围；４ 隧道洞口初期开挖；５ 坡顶卸载裂缝；
６ 剪切蠕滑带；７ 降雨入涌；８ 滑面贯通；９ 边坡失稳
图１㊀隧道洞口滑坡破坏机理模式
Ｆｉｇ．１㊀Ｆａｉｌｕｒｅｍｅｃｈａｎｉｓｍｍｏｄｅｌｏｆｔｕｎｎｅｌｅｎｔｒａｎｃｅｌａｎｄｓｌｉｄｅ
２㊀典型案例分析
２．１㊀工程概况
某隧道原设计为分离式中隧道，其左洞长８４６ｍ㊁右洞长６５５ｍ㊂根据现场勘察和详勘资料揭示，
传真机怎么发传真该隧道右线出口隧道洞顶主要为Ｑａｌ＋ｐｌ
４黄土状粉土，
洞室围岩以强风化砂岩㊁泥质砂岩为主，夹砂质泥
岩㊁泥岩薄层，局部以互层状产出，属软岩 较软岩，
６
端午节 英文３２
the ball
倪向龙等：黄土地区高速公路隧道洞口滑坡变形机理分析２０２１年第５期
country road
corpus岩体较破碎，岩石节理裂隙较发育㊂地下水以基岩裂隙水为主，水量较贫乏㊂隧道洞口为Ⅴ级围岩浅埋段，长度１２０ｍ㊂隧道右线出口段的现场施工采用台阶法开挖，目前仅开挖了上台阶，并完成了上台阶初期支护，暗洞施工段长度７６．５ｍ，上台阶高约
６．２ｍ，初期支护采用锚喷支护㊂在强降雨的影响下，导致已完成开挖及初支段隧洞内拱顶下沉㊁周边严重收敛，挤压钻孔简易台车，套拱出现贯通裂缝，喷射砼表面剥落开裂，初支拱架变形，临时支撑无法顺利施作，洞外右侧临时住房地表开裂下陷，隧道顶部缓坡地表出现多处开裂，右侧拱脚处水平岩层挤压错位㊁施工安全无法保证（图２）㊂
beefcake
blue box图２㊀隧道洞口滑坡现场
Ｆｉｇ．２㊀Ｏｎ－ｓｉｔｅｉｎｖｅｓｔｉｇａｔｉｏｎｍａｐｏｆｔｕｎｎｅｌｅｎｔｒａｎｃｅｌａｎｄｓｌｉｄｅ
㊀㊀经调查发现，该隧道洞口段有滑坡发生，滑坡体发育于北侧黄土斜坡上，滑坡体后缘发育有落水洞，
拟建线路从滑坡体中部通过，滑坡体黄土沿基岩面产生滑动㊂从地貌上看，地形整体呈上陡 中缓
下陡之势，滑坡后壁高陡，坡度３０ʎ ３６ʎ，高３０ｍ；滑坡前部地形局部较为平缓，坡度约１０ʎ；滑坡前部坡体邻空，坡度较陡，平均坡度３０ʎ㊂受雨水冲刷，
坡面波状起伏，沟壑发育㊂因隧道开挖以及暂时性流水冲蚀前缘有基岩出露，标程＋９８６ ＋９９２ｍ，滑体后部西侧沿滑坡周界有一条小型冲沟，受冲沟冲刷，滑坡后壁西侧出现崩塌，并有基岩出露，滑坡体上长满杂草㊂滑体中部宽度约１５０ｍ，轴长约１３０ｍ，厚度约１７ｍ，体积约３３ˑ１０４ｍ３，为黄土切基岩工程滑坡，主滑方向为７３ʎ，与路线呈３９ʎ相交（图
３）㊂
图３㊀隧道洞口滑坡剖面
Ｆｉｇ．３㊀Ｌａｎｄｓｌｉｄｅｐｒｏｆｉｌｅｏｆｔｕｎｎｅｌｐｏｒｔａｌ
２．２㊀滑坡变形及形成机制分析谷歌英文版>翻译中文
该隧道洞口的洞顶为第四系黄土覆盖，下伏为基岩层，根据现场调查发现，上部黄土有明显的蠕滑特征，下部基岩节理裂隙发育，剪切脆性破坏特征明显，因此根据滑坡的变形发展过程和运动特征，将该
滑坡的形成机理与上述总体滑坡形成机理相匹配，并采用有限元模拟方法对形成机理进行分析㊂
２．２．１㊀有限元计算基本假设
对该隧道洞口滑坡进行数值模拟的前提下，需要做如下假设：①岩土体为理想的弹塑介质；②岩土
７
３２
２０２１年第５期
煤炭科学技术
第４９卷
体的材料属性为各向同性均匀介质；③该边坡所受的力为垂直的地应力场；④降雨对该滑坡的诱发影响仅为坡面渗流，无地下水的影响；⑤隧道开挖后，每个剖面隧道及围岩的受力和变形分析都是平面应变问题；⑥由于隧洞洞口开挖深度较小，初期支护对坡体影响较小，不考虑初期支护的加固作用㊂
２．２．２㊀计算范围与边界条件
根据目前发生的滑坡边界确定建立３维地质模
型的范围，同时根据圣维南原理，隧道开挖后的应力应变在隧道中心点３ ５倍的隧道开挖宽度的范围内有影响，因此模拟的范围为宽２００ｍ，高９０ｍ㊂计算范围如图４所示㊂在坡体底部与边界进行垂直于水平约束［２２］㊂
２．２．３㊀计算参数的选取
计算模型中物理参数的取值根据场地的勘察报告并对照‘公路隧道设计规范“进行取值，物理参数见表２㊂表２㊀数值模拟物理参数取值
Ｔａｂｌｅ２㊀Ｎｕｍｅｒｉｃａｌｌｙｓｉｍｕｌａｔｅｄｏｆｐｈｙｓｉｃａｌｐａｒａｍｅｔｅｒｓ
岩性弹性模量Ｅ／ＭＰａ泊松比υ容重γ／（ｋＮ㊃ｍ３）
黏聚力Ｃ／ｋＰａ内摩擦角φ／（ʎ）渗透系数ｋ／（ｍ㊃ｄ－１）孔隙比ｅ上覆黄土５０．３０２０３０３６０．８６４００．６下伏基岩
４ˑ１０４
０．４５
２１
３００００
６０
０．０２６４
０．４
２．２．４㊀有限元模型的建立
利用Ｍｉｄａｓ软件进行数值模拟，所采用的本构模型为摩尔－库伦准则，模型有５个参数，即控制弹性行为的２个参数：弹性模量（Ｅ）和泊松比（υ）及控制塑性行为的３个参数：有效黏聚力ｃ㊁内摩擦角φ和剪胀角ψ㊂而对于降雨入渗的模拟采用属性赋值的非饱和特性函数，设置节点水头与节点流量控制降雨入渗速度，其中降雨强度根据实际统计为１３０ｍｍ／ｄ㊂
经过有限元网格划分产生２０６８０个应变等参元和４３９４个节点，有限元模型如图４所示
㊂
图４㊀隧道洞口边坡有限元模型
Ｆｉｇ．４㊀Ｆｉｎｉｔｅｅｌｅｍｅｎｔｍｏｄｅｌｏｆｔｕｎｎｅｌｅｎｔｒａｎｃｅｓｌｏｐｅ
２．２．５㊀结果分析
１）隧道洞口未进行施工前，坡体整体处于稳定
状态，总体位移为零，因此在数值模拟的过程中原始状态位移清零㊂但是该边坡坡面松散堆积，由图５中可以看出坡面沿坡脚方向有最大主应力集中，容易产生沿坡面的浅层破坏㊂同时可以看出黄土地区上覆黄土通常工程性质较差，在原始状态下，坡体易受外力作用而发生破坏
㊂图５㊀隧道洞口边坡原始状态下最大主应力集中Ｆｉｇ．５㊀Ｍａｘｉｍｕｍｐｒｉｎｃｉｐａｌｓｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｄｉａｇｒａｍｏｆ
ｔｕｎｎｅｌｅｎｔｒａｎｃｅｓｌｏｐｅｉｎｏｒｉｇｉｎａｌｓｔａｔｅ
２）隧道开始施工后，为了便于开挖隧道，需要
在坡脚处开挖掌子面，坡脚开挖的工作面近乎直立㊂通过数值模拟开挖过程发现，坡脚处应力集中，坡顶有塑性应变集中区，其范围向下延伸㊂坡面的最大位移为２．７ｃｍ，对其进行稳定性分析，边坡内部会产生１条贯通的塑性应变带，但是坡体整体会处于稳定状态㊂该状态下边坡前缘处于锁固状态，后缘处于蠕滑状态，后缘位移大于前缘位移㊂
３）隧道洞口开始开挖后，破坏了坡体的整体的
应力平衡，使岩土体内部发生应力调整㊂坡脚处变形不断向上扩展，其影响范围符合圣维南原理，该处最大位移为１５ｃｍ，且坡脚已经开始出现破坏，同时坡体后缘有效塑性应力比坡脚开挖时的塑性区更加明显［２３］㊂根据塑性区判断，坡体后缘可能有拉张裂缝出现，坡体已经有蠕滑发生的征兆［２４］㊂
４）隧道洞口施工过程中，由于强降雨的发生诱
发上覆黄土蠕滑㊂因此根据实际降雨发生的时长建立降雨入渗模型，分析降雨持续时间对坡体饱和度与孔隙压力的影响，间接分析其对边坡的滑移影响㊂
８
３２
