公路边坡稳定性分析范文

摘要:在人为扰动、地震、降雨等因素中，降雨是滑坡主要的诱发因素。以玉溪新平大戛高速公路边坡为研究工点，结合地质勘测、钻探、物探、现场调查、室内试验、资料收集等数据，采用无限斜坡模型、瑞典条分法及极限平衡法计算分析降雨条件下斜坡的流向流量及其稳定系数，研究斜坡开挖前后的稳定性及预测水流冲刷严重的区域。结果表明:开挖前在水位高度等于土层厚度的极限饱和工况下，斜坡处于稳定状态;开挖后按1:1．75的坡比，斜坡亦处于稳定状态。所以坡体开挖施工是安全的，但需做好坡面排水及防护工作。

关键词:道路工程;降雨;公路边坡;无限斜坡;极限饱和;GIS

0引言

随着全球气候变暖，降雨耦合人为扰动、地震、改移河道、工程地质等因素，滑坡等自然灾害频发［1］。如1991年9月，云南昭通市盘河乡头寨沟，因持续强降雨的影响，导致山体崩塌，死216人、伤7人、经济遗失近百万元。2011年7月，受强降雨影响，陕西略阳县爆发山体滑坡，致使26400人受灾，27480亩农作物受损，18人死亡，2人失踪，5人受伤，多个乡镇受灾严重［2］。2016年5月，福建三明市泰宁县开善乡，因暴雨诱发山体滑坡，致使35人遇难，1人失联。2017年5月，暴雨影响，斯里兰卡爆发了大规模的泥石流、滑坡，死亡人数达91人，另有100人失踪，约52万人受灾。2017年8月，受台风影响，暴雨频发，致使贵州省纳雍县张家湾镇发生山体滑坡，截止目前已有15人遇难，8人受伤。2017年8月，正好是印度多暴雨季节，在印度北部喜马偕尔邦突发的山体滑坡，导致46人遇难，多人受伤。从以上滑坡实例可知，降雨是诱发滑坡的主要因素。降雨入渗导致坡体强度降低，当抗滑力不足以支撑下滑力时，斜坡发生失稳，给人们的生命、财产安全带来极大的威胁。因此笔者主要考虑降雨因素的影响，结合相应的斜坡稳定计算模型，对斜坡的稳定性进行分析与预测。很多学者对降雨诱发型斜坡稳定性进行了探索与研究，如:均匀岩质边坡稳定性分析的有限元算法［3］;基于动态和整体分析边坡稳定性的强度折减法［4］;人为划分条块的Bishop法［5］;假定滑动面的瑞典条分法［6］;基于GIS及相关模型［7-14］在斜坡稳定性分析中的运用等。然而以上研究都是假设地下水位在滑动面(基岩)以上，但实际情况可能是地下水位在滑动面(基岩)以下，此时雨水下渗时主要考虑动水压力，随着降雨入渗的深入，基岩透水性较弱，基岩以上土层逐渐趋于饱和，此时可能存在一定的静水压力，当降雨持续进行时，降雨转换成地下水来不及排泄，使基岩以上土层达到饱和，若土层内水位高度等于土层厚度，这种情况就是极限饱和工况。如果斜坡在这种极限饱和工况下能保持稳定，则其他状况也能保持稳定。基于此，结合实验工点地质勘测、钻探、物探、现场调查、室内试验、资料收集等获取试验工点的相关数据，研究开挖前极限饱和工况状态下斜坡稳定性及开挖后坡体稳定性，为边坡的开挖施工提供一定的理论依据。

1工点概况

本工点边坡里程为K65+640～K65+940，全长300m，位于戛洒镇达哈村附近。坡脚处为达哈河，该河流为常年流水，河床一般宽约6～8m，最大宽度可达12m，水深0．8～1．5m。区域上属元江水系，河床纵坡平缓，流速慢，勘察期间流速约1．5m/s，流量受大气降雨影响较大。河流对研究坡体的地下水位影响较小。拟建的路线里程地形起伏稍大，斜坡海拔722～802m，相对高差约80m。地处青藏川滇歹字型构造体系，川滇经向构造体系及南岭纬向构造体系的交接地带，应力集中，区域构造现象极为复杂。据有关资料表明:实验工点位于新平-新化褶皱区及水塘断裂附近，受褶皱及断裂带影响岩体较破碎，地层为三叠系上统干海子组(T3g)砾岩、砂岩。上覆发育第四系冲洪积层(Q4dl+el)，层厚2～4m，主要含粉质黏土。地下水位较深，远低于基岩底部。其地下水位受季节性影响较大，地下水为第四系孔隙水类型及基岩裂隙水类型。孔隙水多赋存于松散土体中，孔隙潜水的形式较多，水位受季节影响大;基岩裂隙水赋存于基岩裂隙中，基岩富水性较差，水量不大。该段边坡位于果林中，果林常年滴管，地表水常年冲刷。根据该地区水文气象数据，年平均气温为17．3～18．6℃，最热月平均气温21．3～23．1℃，最冷月平均气温10．6～11．7℃;年平均降雨量为970．7～1301．7mm。

2计算模型

该工点斜坡的地下水位远低于粉质黏土层，笔者采用无限斜坡模型［8，13-14］结合GIS计算开挖前坡体的稳定系数;使用极限平衡法［6］、瑞典条分法［6］计算开挖后，斜坡按1∶1．75坡比设计的稳定系数。然后根据实验工点的TIN数据，通过GIS的3D表面分析工具、数据管理工具等，得出坡体表面径流的流向流量，预测因雨水汇集而导致冲刷严重的区域。

2．1无限斜坡模型的选用SHALSTAB模型、TRIGRS模型、降雨入渗模型［8］、基于GIS的无限斜坡耦合模型［13-14］，这4种模型都是降雨型无限斜坡稳定性的计算模型，如图2。对比分析4种模型，选取比较符合假设的计算模型。

2．1．1SHALSTAB模型SHALSTAB模型主要考虑土体为饱和状态时，坡体的稳定性。

2．1．2TRIGRS模型TRIGRS模型主要考虑降雨入渗引起孔隙水压力发生改变这一现象，计算坡体的稳定系数。

2．1．3降雨入渗模型显然上述的SHALSTAB、TRIGRS模型不能很好的描述降雨条件下，无限斜坡的稳定性。

2．1．4基于GIS的无限斜坡耦合模型降雨入渗模型没有考虑静水压力，因此笔者在计算时作如下假设:粉质黏土基岩透水性较差，渗透能力可以忽略不计。

2．2相关参数以及饱和度计算公路边坡相关参数，而斜坡坡度θ由GIS的spatialanalysisttools提取;假定极限饱和状态时水位深度Zw与土层厚度一致。

3讨论

3．1降雨入渗-水位变化分析斜坡稳定性的问题，考虑降雨入渗引起的动水压力、静水压力变化。降雨条件下，雨水一部分被植被截留，当地表结皮或者下渗土层饱和后，一部分雨水通过地表径流流走，最后一部分通过植被根系、生命体通道、坡体裂缝下渗，下渗的水流中，一部分渗入土体内成包气带内，一部分转换成地下水，引起地下水位变化，地下水位不同，对土体的作用效果也就不同。地下水位位于基岩以上，降雨时，不仅受到动水压力的作用，而且还会受到静水压力的作用［15］。如果地下水位位于基岩以下，静水压力可以忽略，但随着降雨入渗的不断增加，土体饱和程度发生了明显的变化。考虑基岩渗透性弱，土体内水位上升，可能还存在一定的静水压力。

3．2坡体稳定性分析无限斜坡模型的计算结果是基于钻探、室内试验等土体参数;开挖后的边坡稳定系数计算是基于设计值;流向流量的预测结果是基于坡体表面的TIN数据。参照王学鹏［16］对降雨入渗模型稳定系数值的划分:1～1．25(基本稳定);1．25～1．5(稳定);1．5以上(极稳定);基于GIS的无限斜坡耦合模型稳定系数值都在1．5以上(极稳定)。说明开挖前坡体处于稳定状态，这与实际情况相符合。开挖后，瑞典条分法及极限平衡法按1∶1．75的坡比计算，稳定系数为1．492、1．179，则边坡亦处于稳定。根据计算结果及水流流向从微观和宏观两方面对斜坡稳定性进行分析，微观方面:随着降雨持续，达到了诱发渗流的水头梯度，渗流发生，导致土层内基质吸力降低，含水率升高，土体由非饱和趋向于饱和［17］发展，内摩擦角、黏聚力发生变化，土体抗剪强度降低，坡体下滑力逐步增大;宏观方面:随着坡面雨水汇集，其流向流量如图7，水流对坡体的冲刷作用凸显，在坡体表面易形成沟壑。随着冲刷次数的累积，坡体表面的原状土渐渐流失，逐渐形成沟壑。若沟壑未进一步发展，或降雨量减小直至停止以及坡体自身应力调节作用，坡体又再次达到平衡状态。但如果降雨持续进行，雨水入渗将加剧，引起坡体内部土体抗剪强度降低，再加上外部水流的冲刷，带走坡体表面的原状土，应力平衡终遭到破坏，最终将导致斜坡稳定性降低，直至发生失稳。

3．3坡面防护分析结合以上分析及开挖后的边坡稳定系数，需要注意坡面开挖后的稳定系数值。

4结论

通过降雨入渗模型、基于GIS的无限斜坡耦合模型、瑞典条分法及极限平衡法计算出玉溪新平大戛高速开挖前后边坡的稳定系数，得出以下结论:1)在水位高度等于土层厚度的极限饱和工况下，降雨入渗模型及基于GIS的无限斜坡耦合模型，预测结果为斜坡处于稳定状态，符合现场的实际情况。2)根据1∶1．75的设计坡比，按照瑞典条分法及极限平衡法计算的稳定系数值为:1．429、1．179，则开挖后坡体处于稳定状态。3)根据流向流量图，可以预测坡面的严重冲刷区域。在坡面防护时，水流冲刷严重区域，应注重坡面排水。4)施工中可以采取框格梁、锚杆、挡土墙、抗滑桩等，确保施工安全。5)结合坡体开挖前后稳定系数分析，对公路边坡的稳定性分析与预测提供了新思路。

作者：何玉琼1；李春明1；杨兴华2；黄犀1 单位：1．昆明理工大学，交通工程学院，2．云南建设基础设施投资股份有限公司