暴雨情形下含裂隙土质边坡的瞬态稳定性

王晓磊，王旭春，管晓明，张 鹏

(1.青岛理工大学土木工程学院，山东 青岛 266033;2.北京交通大学土木建筑工程学院，北京 100044;3.中国矿业大学力学与建筑工程学院，北京 100083)

降雨诱发边坡失稳占据了滑坡灾害相当的比重，造成了大量的人员伤亡与财产损失。《中国重大地质灾害实例分析》［1］列举了27个国家重大地质灾害，其中有15个就是由于降雨所引发。《中国典型滑坡》［2］一书中实录了90 多例滑坡，90%以上发生在雨中和雨后。降雨诱发边坡失稳在我国广泛存在，造成了大量的人员伤亡与财产损失［3-4］。

以往的研究中，膨胀土、杂填土和裂隙黏土等含裂隙土质边坡的稳定分析大多采用传统的静力平衡法和统计比较判别法等，坡体统一按饱和状态进行计算，并且忽略了裂隙的作用［5-6］，这与现场实际情况差别较大。近年来，人们逐渐认识到土体的非饱和状态以及边坡裂隙对边坡稳定性的影响，并做了大量的研究工作，结果表明考虑裂隙和边坡降雨入渗影响的稳定性分析方法较为合理和实用［7-10］。殷宗泽等［11］提出了一种以条分法为基础，近似反映裂隙影响的膨胀土边坡稳定性分析方法，将膨胀土坡划分成3个不同强度指标的亚层。阙云等［12］对受降雨入渗影响的裂隙性黏土边坡进行了研究，分析了降雨条件下不同裂隙深度和开口宽度等因素对裂隙性黏土边坡渗流场的影响。以上研究对含裂隙土质边坡稳定性分析理论的发展起到了极大的促进作用，为后续研究工作奠定了基础。目前研究中多采用在整个边坡均匀布置裂隙的方式，未考虑边坡裂隙的出现形式，并且对降雨条件下坡体存在裂隙时的应变情况研究较少。

本文采用边坡稳定性分析及治理软件Slide 对暴雨情形下含裂隙土质边坡的稳定性进行分析，确定边坡潜在滑动面的位置。根据分析结果布置促使边坡发生推移式滑坡和牵引式滑坡的裂隙，最后应用MIDAS/GTS 软件模拟分析降雨条件下的坡体应变规律和边坡稳定情况。计算过程中考虑降雨边界条件、体积含水率、渗透系数与孔隙压力之间的关系。

1 暴雨情形下裂隙对非饱和土质边坡稳定性的影响

降雨尤其是暴雨是滑坡活动最重要的触发因素，而边坡裂隙的存在为降雨入渗提供了进入内部的通道，改变了降雨的入渗路径，使得降雨这一危险因素对灾害的触发能力得到充分的发挥，可以归纳为以下几点:

a.裂隙增加了入渗能力。边坡裂隙增大了雨水与坡体的接触面积，在降雨过程中，地表径流的雨水大量涌入裂隙之中，增加了雨水入渗量。随着裂隙的发展，接触范围越来越大，并且在裂隙水头的作用下，雨水的入渗能力大幅度增强，裂隙周围土体的含水量迅速增大，边坡土体基质吸力下降，自重增加。基质吸力的下降改变了非饱和土体的应力状态，使土体产生应变，其外在表现为坡体的宏观位移。而土体在饱和之后，自重的增加导致下滑力变大，同时岩土体黏聚力、内摩擦角及基质吸力减小，导致抗剪强度降低。

b.静水压力作用。强降雨下裂隙中充满积水，对边坡产生静水压力，压力P可分解为平行和垂直于滑动面的分力，在此压力作用下边坡将产生向临空面的位移，促使裂隙的发展。同时平行于滑动面的分力增大了边坡的下滑力，而平行于滑移线的分力使坡体正压力减小，降低边坡局部稳定性。裂隙开裂越深，裂隙积水的静水压力越大，静水压力作用越强。

c.裂隙积水渗流作用。裂隙积水使裂隙处存在水头差，水体在裂隙周围土体孔隙中流动，产生动水压力，沿流线方向的单位体积渗透力为

式中:ρw为水的密度，kg/m3;g为重力加速度，m/s2;i为渗流坡降;L为渗流长度，m;h为水头损失，m。动水压力对土体骨架产生拖曳力，导致土体中的应力与变形发生变化。一旦边坡渗流通道形成，土体在此渗透力的作用下，边坡松散土体中的土颗粒被冲蚀，加剧边坡的不稳定性。

2 实例分析

2.1 计算模型、初始条件和边界条件

2010年雨季期间，由于几次强降雨，山西省安太堡露天煤矿东北角发生滑坡。本文以该煤矿周边陡帮边坡一勘探剖面为研究对象，具体地层尺寸及分布状况见图1。

图1 安太堡露天煤矿剖面(单位:m)

坡体在计算范围内没有地下水位，所以在计算过程中，初始状态不考虑边坡两侧土体渗流情况。Mein等［13］提出:当降雨量小于土体的入渗能力时，雨水全部入渗，当降雨量超过土体的入渗能力时产生地表径流。降雨过程中，由于考虑的是暴雨情况，降雨强度超过土体的入渗能力，边坡裂隙很快被地表径流充满，因此在降雨过程中，始终认为裂隙处的压力水头为裂隙内水的深度，坡面不考虑积水情况。

本文的降雨入渗模型以土体体积含水率作为控制变量，通过在坡面施加节点流量的方式实现降雨，运用土壤水动力学的基本理论，求解给定降雨条件下的边坡土体瞬态体积含水率的数值解，并将降雨计算结果应用于下一步边坡稳定性的计算之中。边坡土体非稳态渗流控制方程为

模型边界条件:底部完全约束，两侧边界水平位移约束，顶部为自由边界。

2.2 岩土计算参数

计算中所用岩土计算参数来自现场的勘测资料，同一地层取试验平均值，具体见表1。研究区域土体为非饱和土，其体积含水率的变化取决于应力状态的变化和土体特性的变化。体积含水率是孔隙水压力的函数，渗透系数又是体积含水率的函数，所以渗透系数也是孔隙水压的函数。渗透系数和压力水头之间的关系用渗透系数曲线表示，体积含水率和孔隙水压力之间的土水特征曲线采用Van Genuchten模型。降雨过程中体积含水率的变化对土体的强度影响采用毕肖普提出的非饱和土的强度公式［9］:

式中:σ'为有效应力;σ为总应力;ua为孔隙气压力;uw为孔隙水压力;χ为反映土中气相影响的系数。

表1 岩土计算参数

由于滑坡下部为岩质边坡，本文的研究对象主要为发生滑坡的上层土质边坡，限于篇幅只给出粉土的渗透系数曲线和土水特征曲线，见图2和图3。

图2 粉土渗透系数和压力水头的关系

图3 粉土体积含水率和孔隙水压力的关系

2.3 数值模拟方案

针对安太堡矿区露天边坡几次不同程度滑坡均出现在暴雨之后的情况，本次模拟暴雨量为120 mm/d，分析不同降雨持续时间对边坡应变和稳定性的影响规律。为了使边坡裂隙位置更符合实际状态，首先采用Slide 软件对边坡进行分析，确定边坡滑移线位置。然后根据模拟结论将裂隙布置如下:①潜在滑体的上部滑移线位置布置裂隙，模拟边坡上部滑移线附近土体进水软化，发生推移式滑坡;②边坡中下部布置裂隙，模拟雨水进入下部土体，下部滑移线区域土体软化，发生牵引式滑坡，详细模拟方案见表2。

表2 数值模拟方案

2.4 计算结果分析

2.4.1 降雨过程中坡体应变规律

研究区域回填土边坡结构松散，孔隙性和透水性较强，使其更易受到降雨的影响，表层土体可以产生较大的塑性变形和显著的流变变形。以下针对1号边坡顶部平台布置裂隙与顶部平台和坡面中部同时布置裂隙两种情况，对暴雨影响下的坡体应变进行分析:

a.1号边坡顶部平台布置裂隙坡体应变规律。以方案6为例，仅在1号边坡顶部平台出现张拉裂隙时，坡体内部剪应变见图4。随着雨水入渗，降雨持续1 d 时，1号边坡整个回填土区域出现较大的剪应变，最大值出现在1号边坡坡脚位置。由于土体自重作用，坡脚位置区域存在应力集中现象。随着雨水的不断入渗，降雨持续2 d 时，1号边坡回填土层和下层粉土接触面附近区域也出现了较大的剪应变，且塑性区不断向下部平台扩展。到降雨持续3 d时，回填土与粉土接触区域的剪应变范围变大，这是由于雨水入渗使非饱和土含水量增加，基质吸力下降，而基质吸力的降低又使得边坡非饱和物料产生应变，同时滑移线的轮廓变得非常明显。

图4 方案6 坡体剪应变

b.1号边坡顶部平台和坡面中部布置裂隙坡体应变规律。以方案8为例，1号边坡顶部平台和坡面中部同时布置裂隙时，坡体剪应变见图5。降雨持续1 d 时，1号边坡整个回填土层发生较大的剪应变，并且1号边坡坡脚以及回填土与粉土接触区域都出现剪应变极值。下方2号边坡回填土与粉土接触区域也形成了连贯的剪应变带。随着降雨的持续，剪应变区域缩小，极值只出现在1号边坡坡脚下方回填土与粉土接触区域。降雨持续3 d 时，随着渗流的进行水分在坡体内部运移，剪应变极值区几乎消失。

图5 方案8 坡体剪应变

对比以上两种裂隙布置情况，后者降雨入渗影响的速度更快，且影响范围较前者大。降雨持续1 d时，后者的塑性区已经扩展到2号边坡上方区域。而前者始终对1号边坡存在较大的影响。降雨持续3 d后，后者的塑性区从2号边坡一直贯通到1号边坡的坡顶。

2.4.2 降雨过程中稳定系数的变化规律

通过降雨渗流计算，将得到的瞬态孔隙水压力作为荷载应用于边坡的稳定性分析中。采用有限元法的强度折减法，计算边坡稳定性。

a.1号边坡顶部平台布置裂隙稳定性计算。在1号边坡顶部布置深度为0.5 m的裂隙，随着降雨的入渗，不同开口宽度的裂隙所对应的稳定系数如图6 所示。可以看出，裂隙的宽度对于安全系数的影响不大，这是因为开口宽度和裂隙充满积水的速度有关，本文设定裂隙瞬间被雨水充满，因而裂隙的开口宽度对边坡稳定性未有较强的影响;裂隙深度不同时，降雨影响下边坡稳定性变化见图7。随着裂隙深度的开展，稳定系数逐步降低，在1～1.5 d之间稳定系数降低速率最大，且在1.5 d 左右边坡最不稳定，随后稳定系数又出现反弹的趋势。这是由于当雨水入渗到一定程度时土体达到饱和则入渗速度降到很低，此时随着坡体内部水体的运移，边坡稳定性出现回升。

图6 边坡顶部布置0.5 m 裂隙时降雨稳定系数

图7 边坡顶部不同裂隙深度时降雨稳定系数

b.1号边坡顶部平台和坡面中部同时布置裂隙时稳定性计算。考虑1号边坡顶部平台和坡面中部同时布置裂隙，当坡面裂隙开裂较浅时，降雨产生的裂隙积水未能影响到重力荷载作用下坡体最危险滑移线的位置，此时随着雨水的入渗，其稳定系数的降低程度相比于只在坡顶布置裂隙时变化不大，如图8 所示;但当裂隙开裂达到能够影响到潜在滑动面位置深度时，边坡潜在滑动面区域土体遇水软化，强度降低，稳定系数急剧下降。由图5(c)可以看出坡体塑性变形向上延伸，坡体产生牵引式滑坡。

图8 边坡顶部和坡面中部布置裂隙时降雨稳定系数

3 结 论

a.降雨过程中边坡塑性区首先在坡脚位置出现，随着降雨的不断进行，塑性区自坡脚位置向上部不断发展，直到形成贯通的塑性区。

b.降雨过程中稳定系数并不是持续降低。当裂隙周围土体完全饱和之后，土体的渗透系数变得很小，此时降雨入渗变得极为缓慢，坡体内部的水体发生了运移，边坡的稳定性又得到了回升。

c.模拟结果表明在持续降雨1.5 d 左右边坡的稳定系数最小。当裂隙开裂较深时，边坡稳定系数降低到1 以下，存在发生滑坡危险。但分析结果表明，滑坡的发生存在于边坡的表层土体。

d.本文提出的边坡裂隙布置方法体现了边坡实际情况下裂隙的产生以及发展，相比以往坡体均布裂隙的布置方法更接近边坡裂隙的真实情况。

［1］钟立勋.中国重大地质灾害实例分析［J］.中国地质灾害与防治学报，1999，10(3):1-7.

［2］殷跃平.中国典型滑坡［M］.北京:中国大地出版社，2007.

［3］张先发，李明华，张小刚.长江上游暴雨与滑坡崩塌关系［J］.地理，1995(3):102-106.

［4］杨文东.降雨型滑坡特征及其稳定稳定分析研究［D］.武汉:武汉理工大学，2006.

［5］廖世文.膨胀土与铁路工程［M］.北京:中国铁道出版社，1984.

［6］刘特洪.工程建设中的膨胀土问题［M］.北京:中国建筑工业出版社，1997.

［7］姚海林，郑少河，李文斌，等.降雨入渗对非饱和膨胀土边坡稳定性影响的参数研究［J］.岩石力学与工程学报，2002，21(7):1034-1039.

［8］姚海林，郑少河，陈守义，等.考虑裂隙及雨水入渗影响的膨胀土边坡稳定性分析［J］.岩土工程学报，2001，23(5):606-609.

［9］陈铁林，邓刚，陈生水，等.裂隙对非饱和土边坡稳定性的影响［J］.岩土工程学报，2006，28(2):210-215.

［10］李焕强，孙红月，孙新民，等.降雨入渗对边坡性状影响的模型实验研究［J］.岩土工程学报，2009，31(4):589-594.

［11］殷宗泽，徐彬.反映裂隙影响的膨胀土边坡稳定性分析［J］.岩土工程学报，2011，33(3):454-459.

［12］阙云，胡昌斌，姚晓琴.降雨入渗对裂隙性黏土边坡稳定性作用机理的分析［J］.福州大学学报，2009，37(3):423-429.

［13］MEIN R G，LARSON C L.Modelling infiltration during a steady rain［J］.Water Resource Research，1973，9(2):384-393.