第9章边坡稳定性分析

第9章边坡稳定性分析
学习指导：本章介绍了边坡的破坏类型，即：岩崩和岩滑；着重介绍了边坡稳定性分析与评价基本方法，包括圆弧法岩坡稳定分析、平面滑动法岩坡稳定分析、双平面滑动岩坡稳定分析、力多边形法岩坡稳定分析及近代理论计算法；介绍了岩坡处理的措施。

重点：1边坡的变形与破坏类型；
2影响边坡稳定性的因素；
3边坡稳定性分析与评价。

9.1 边坡的变形与破坏类型
9.1.1 概述
随着社会进步及经济发展，越来越多地在工程活动中涉及边坡工程问题，通过长期的工程实践，工程地质工作者已对边坡工程形成了比较完善的理论体系，并通过理论对人类工程活动，进行有效地指导。

近年来，随着环境保护意识的增加及国际减轻自然灾害十年来的开展，人类已认识到：边坡诞生不仅仅是其本身的历史发展，而是与人类活动密切相关；人类在进行生产建设的同时，必须顾及到边坡的环境效应，并且把人类的发展置于环境之中，因而相继开展了工程活动与地质环境相互作用研究领域，在这些领域中，边坡作为地质工程的分支之一，一直是人们研究的重点课题之一。

在水电、交通、采矿等诸多的领域，边坡工程都是整体工程不可分割的部分，为保证工程运行安全及节约经费，广大学者对边坡的演化规律、边坡稳定性及滑坡预测预报等进行了广泛研究。

然而，随着人类工程活动的规模扩大及经济建设的急剧发展，边坡工程中普遍出现了高陡边坡稳定性及大型灾害性滑坡预测问题。

在我国，目前的露天采矿的人工边坡已高达300—500m，而水电工程中遇到的天然边坡高度已达500—1000米，其中涉及的工程地质问题极为复杂，特别是
在西南山区，边坡的变形、破坏极为普遍，滑坡灾害已成为一种常见的危害人民生命财产安全及工程正常运营的地质灾害。

因此，广大工程地质和岩石力学工作者对此问题进行了长期不懈的探索研究，取得了很大的进展；从初期的工程地质类比法、历史成因分析法等定性研究发展到极限平衡法、数值分析法等定量分析法，进而发展到系统分析法、可靠度方法灰色系统方法等不确定性方法，同时辅以物理模拟方法，并且诞生了工程地质力学理论、岩(土)体结构控制论等，这些无疑为边坡工程及滑坡预报研究奠定了坚实的基础，为人类工程建设做出了重大贡献。

在工程中常要遇到岩坡稳定的问题，例如在大坝施工过程中，坝肩开挖破坏了自然坡脚，使得岩体内部应力重新分布，常常发生岩坡的不稳定现象。

又如在引水隧洞的进出口部位的边坡、溢洪道开挖的边坡、渠道的边坡以及公路、铁路、采矿工程等等都会遇到岩坡稳定的问题。

如果岩坡由于力过大和强度过低，则它可以处于不稳定的状态，一部分岩体向下或向外坍滑，这一种现象叫做滑坡。

滑坡造成危害很大，为此在施工前，必须做好稳定分析工作。

岩坡不同于一般土质边坡，其特点是岩体结构复杂、断层、节理、裂隙互相切割，块体极不规则，因此岩坡稳定有其独特的性质。

它同岩体的结构、块体密度和强度、边坡坡度、高度、岩坡表面和顶部所受荷载，边坡的渗水性能，地下水位的高低等有关。

岩体内的结构面，尤其是软弱结构面的的存在，常常是岩坡不稳定的主要因素。

大部分岩坡在丧失稳定性时的滑动面可能有三种。

一种是沿着岩体软弱岩层滑动；另一种是沿着岩体中的结构面滑动；此外，当这两种软弱面不存在时，也可能在岩体中滑动，但主要的是前面两种情况较多。

在进行岩坡分析时，应当特别注意结构面和软弱层的影响。

软弱岩层主要是粘土页岩、凝灰岩、泥灰岩、云母片岩、滑石片岩以及含有岩盐或石膏成分的岩层。

这类岩层遇水浸泡后易软化，强度大大地降低，形成软弱层。

在坚硬的岩层中(如石英岩、砂岩等等)应当查明有无这类软弱夹层存在。

结构面包括沉积作用的层面、假整合面、不整合面；火成岩侵入结构面以及冷缩结构面；变质作用的片理，构造作用的断裂结构面等等。

岩质边坡稳定分析时，应当研究岩体中应力场和各种结构面的组合关系。

岩坡的滑动就是在应力作用下岩体破坏了平衡而沿着某种面(很可能是结构面)产生的。

岩体的应力是由岩体重量、渗透压力、地
质构造应力以及外界因素，如地震惯性力、风力、温度应力等所形成的边坡剪应力，这种剪应力超过结构面的抗剪强度就促使岩体沿着结构面滑动。

有时沿某一结构面滑动，有时沿着多种结构面所组合的滑动面滑动。

通常以后者为多数。

结构面中如夹有粘土或其它泥质充填物，则就成为软弱结构面。

地质构造作用形成的断裂和节理在地壳表层是最多的，这种结构面往往都夹有粘土或泥质充填物，遇水浸泡后，结构面中的软弱充填物就容易软化，强度大大地降低，促使岩坡沿着它发生滑动。

因此，岩坡分析中，对结构面，特别是软弱结构面的类型、性质、组合形式、分布特征以有及由各种软弱面切割后的块体形等进行仔细分析是十分重要的。

9.1.2 岩坡的破坏类型
岩坡的破坏类型从形态上来看可分为岩崩和岩滑两种。

岩崩一般发生在边坡过陡的岩坡中，这时大块的岩体与岩坡分离而向前倾倒，如图9-1(a)所示，或者坡顶岩体因某种原因脱落而在坡脚下堆积，见图9-1(b)、(c)，它经常产生于坡顶裂隙发育的地方。

其起因或由于风化等原因减弱了节理面的凝聚力，或由于雨水进入裂隙产生水压力所致；或者也可能由于气温变化、冻融松动岩石的结果；其它如植物根造成膨胀压力、地震、雷击等都可造成岩崩现象。

岩滑是指一部分岩体沿着岩体较深处某种面的滑动。

岩滑可分为平面滑动、楔形滑动以及旋转滑动。

平面滑动是一部分岩体在重力作用下沿着某一软面(层面、断层、裂隙)的滑动，见图9-2(a)，滑动面的倾角必大于该平面的内摩擦角。

平面滑动不仅滑体克服了底部的阻力，而且也克服了两侧的阻力。

在软岩中(例如页岩)，如底部倾角远陡于内摩擦角，则岩石本身的破坏即可解除侧边约束，从而产生平面滑动。

而在硬岩中，如果不连续面横切坡顶，边坡上岩石两侧分离，则也能发生平面滑动。

楔形滑动是岩体沿两组(或两组以上)的软弱面滑动的现象，见图9-2(b)。

在挖方工程中，如果两个不连续面的交线出露，则楔形岩体失去下部支撑作用而滑动。

法国马尔帕塞坝的崩溃(1959年)就是岩基楔形滑动的结果。

旋转滑动的滑动面通常呈弧形状，见图9-2(c)，这种滑动一般产生于非成层的均质岩体中。

岩坡的滑动过程一般可分为三个阶段。

初期是蠕动变形阶段，这
一阶段中坡面和坡顶出现拉张裂缝并逐渐加长和加宽，滑坡前缘有时出现挤出现象，地下水位发生变化，有时会发出响声。

第二阶段是滑动破坏阶段，此时滑坡后缘迅速下陷，岩体以极大的速度向下滑动，此一阶段往往造成极大的危害。

最后是逐渐稳定阶段，这一阶段中，疏松的滑体逐渐压密，滑体上的草木逐渐生长，地下水渗出由浑变清等。

在进行岩坡稳定性分析时，首先应当查明岩坡可能的滑动类型，然后对不同类型采用相应的分析方法。

严格而言，岩坡滑动大多属空间滑动类型，然后对只有一个平面构成的滑裂面或者滑裂面由多个平面组成而这些面的走向又大致平行者，且沿着走向长度大于坡高时，则也可按平面滑动进行分析，其结果偏于安全方面，在平面分析中，常常把滑动面进行稳定验算。

本章从第四节起将分别阐述各种分析方法。

经验证明，许多滑坡的发生都与岩体内的渗水作用有关，这是由于岩体内渗水后岩石强度恶化和应力增加的缘故。

因此，做好岩坡的排水工作是防止滑坡的手段之一。

意大利瓦依昂(Vajont)水库岩坡滑动而造成的事故是闻名于全世界的。

水库的岸坡由分层的石灰岩组成，水库蓄水后在1960年10月就发现上坡附近有主要裂隙，同时直接在沿河的陡坡上曾经发生过一
9.2 影响边坡稳定性的因素
影响边坡稳定性的因素主要有内在因素和外部因素两方面，内在因素包括组成边坡的地貌特征、岩土体的性质、地质构造、岩土体结构、岩体初始应力等。

外部因素包括水的作用、地震、岩体风化程度、工程荷载条件及人为因素。

内在因素对边坡的稳定性起控制作用，外部因素起诱发破坏作用。

1)岩土性质和类型
岩性对边坡的稳定及其边坡的坡高和坡角起重要的控制作用。

坚硬完整的块状或厚层状岩石如花岗岩、石灰岩、砾岩等可以形成数百米的陡坡，如长江三峡峡谷。

而在淤泥或淤泥质软土地段，由于淤泥的塑性流动，几乎难以开挖渠道，边坡随挖随塌，难以成形。

黄土边坡在干旱时，可以直立陡峻，但一经水浸土的强度大减，变形急剧，滑动速度快，规模和动能巨大，破坏力强且有崩塌性。

松散地层边坡的坡度较缓。

不同的岩层组成的边坡，其变形破坏也有所不同，在黄土地区，边坡的变形破坏形式以滑坡为主；在花岗岩、厚层石灰岩、沙岩地区则以崩塌为主；在片岩、板岩、千枚岩地区则往往产生表层挠曲和倾
倒等蠕动变形。

在碎屑岩及松散土层地区，则产生碎屑流或泥石流等。

2)地质构造和岩体结构的影响
在区域构造比较复杂，褶皱比较强烈，新构造运动比较活动的地区，边坡稳定性差。

断层带岩石破碎，风化严重，又是地下水最丰富和活动的地区极易发生滑坡。

岩层或结构的产状对边坡稳定也有很大影响，水平岩层的边坡稳定性较好，但存在陡倾的节理裂隙，则易形成崩塌和剥落。

同向缓倾的岩质边坡(结构面倾向和边坡坡面倾向一致，倾角小于坡角)的稳定性比反向倾斜的差，这种情况最易产生顺层滑坡。

结构面或岩层倾角愈陡，稳定性愈差。

如岩层倾角小于10°～15°的边坡，除沿软弱夹层可能产生塑性流动外，一般是稳定的；大于25°的边坡，通常是不稳定的；倾角在15°～25°的边坡，则根据层面的抗剪强度等因素而定。

同向陡倾层状结构的边坡，一般稳定性较好，但由薄层或软硬岩互层的岩石组成，则可能因蠕变而产生挠曲弯折或倾倒。

反向倾斜层状结构的边坡通常较稳定，但垂直层面或片理面的走向节理发育且顺山坡倾斜，则亦易产生切层滑坡。

3)水的作用
地表水和地下水是影响边坡稳定性的重要因素。

不少滑坡的典型实例都与水的作用有关或者水是滑坡的触发因素，处于水下的透水边坡将承受水的浮托力的作用，而不透水的边坡，将承受静水压力；充水的张开裂隙将承受裂隙水静水压力的作用；地下水的渗流，将对边坡岩土体产生动水压力。

水对边坡岩体还产生软化或泥化作用，使岩土体的抗剪强度大为降低；地表水的冲刷，地下水的溶蚀和潜蚀也直接对边坡产生破坏作用。

不同结构类型的边坡，有其自身特有的水动力模型。

(1) 静水压力
作用于边坡的静水压力主要包括两种情况：其一是当边坡被水库淹没时，库水对边坡面所产生的静水压力；其二是当裂隙岩石边坡的张裂隙充水时，裂隙中的水压力。

①边坡坡面上的静水压力：当边坡被水淹没，而边坡的表部相对不透水时，坡面上将受一定的静水压力，静水压力的方向与坡面正交。

当边坡的滑动面(软弱结构面)的倾角θ小于坡角α时，则坡面静水压力传到滑动面上的切向分量为抗滑力，对边坡稳定有利。

当α
θ>时，则切向分量为下滑力，则不利于边坡的稳定。

②边坡裂隙静水压力：有张裂隙发育的岩石边坡以及长期干旱的裂隙粘土边坡，如果因降雨或地下水活动使裂隙充水，则裂隙面将承受静水压力(图9-5)。

静水压力的作用方向与裂隙面相垂直，其大小
P w(取与裂隙水水头有关。

对部分充水的高角度裂隙，裂隙静水压力
w
由于裂隙水活动的不规律性，岩体中的地下水位通常不是圆滑的曲线。

在相邻裂隙的地下水位不同时，地下水位高的裂隙较地下水位低的裂隙承受较大的静水压力，这种静水压力的差别，有时是使边坡失稳的原因之一。

由于地下水出口节理裂隙敞开情况不同，也影响裂隙水压力的大小，因而影响边坡的稳定。

如图9-6所示，出口节理张开，地下水位低，裂隙水压力小；出口节理闭合，透水性差，则地下水位高，裂隙水压力大。

如作用在岩块底部滑面上的静水压力，有时可使覆岩块隆胀(静水压力等于上覆岩块重)，而使边坡稳定严重恶化。

(2)浮托力
处于水下的透水边坡，承受浮托力的作用，使坡体的有效重量减轻，这对边坡的稳定不利。

不少水库周围松散堆积层边坡，在水库蓄水时发生变形，浮托力的影响是原因之一。

对处于极限稳定状态，依靠坡脚岩体重量保持暂时稳定的边坡，坡脚被水淹没后，浮托力对边坡稳定的影响就更加显著。

(3)动水压力
动水压力是地下水在流动过程中所施加于岩土体颗粒上的力。

它
是一种体积力，其数值为：
I V D w γ= (9-2)
其中：V 流动水体体积；
w γ 水的块体密度；
I 水力梯度。

动水压力的方向和水流方向平行，在近似计算中，多假定与地下
水面或滑面平行，如果动水压力方向和滑体滑动方向不一致，则应分
解为垂直和平行于滑面的两个分量参与稳定计算。

在边坡稳定的实际
计算中，由于渗流方向不是定值，且水力梯度不易精确确定，一般则
作简化假定，以采用不同的滑体块体密度将动水压力的影响计入。

即
在地下水位以下静水位以上有渗流活动的滑体，计算下滑力时，采用
饱和块体密度；计算抗滑力时，采用浮块体密度。

4) 工程荷载 在水利水电工程中，工程荷载的作用影响边坡的稳定性。

例如，
拱坝坝肩承受的拱端推力、边坡坡顶附近修建大型水工建筑物引起的
坡顶超载、压力隧洞内水压力传递给边坡的裂隙水压力、库水对库岸
的浪击淘涮力、为加固边坡所施加的力，如预应力锚杆时所加的预应
力等都影响边坡的稳定性。

由于工程的运行也可能间接地影响边坡的
稳定，例如由引水隧洞运行中的水锤作用，使隧洞围岩承受超静水荷
载，引起出口边坡开裂变形等。

5）地震作用
地震对边坡稳定性的影响表现为累积和触发(诱发)等两方面效
应。

(1)累积效应
边坡中由地震引起的附加力S 的大小，通常以边坡变形体的重量
W 与地震振动系数k 之积表示(S =kW )。

在一般边坡稳定性计算中，将
地震附加力考虑为水平指向坡外的力。

但实际上应以垂直与水平地震
力的合力的最不利方向为计算依据。

总位移量的大小不仅与震动强度
有关，也与经历的震动次数有关，频繁的小震对斜坡的累进性破坏起
着十分重要的作用，其累积效果使影响范围内岩体结构松动，结构面
强度降低。

(2)触发(诱发)效应
触发效应可有多种表现形式。

在强震区，地震触发的崩塌、滑坡
往往与断裂活动相联系。

高陡的陡倾层状边坡，震动可促进陡倾结构面(裂缝)的扩展，并引起陡立岩层的晃动。

它不仅可引发裂缝中的空隙水压力(尤其是在暴雨期)激增而导致破坏，也可因晃动造成岩层根部岩体破碎而失稳。

碎裂状或碎块状边坡，强烈的震动(包括人工爆破)甚至可使之整体溃散，发展为滑塌式滑坡。

结构疏松的饱和砂土受震液化或敏感粘土受震变形，也可导致上覆土体产生滑坡。

海底斜坡失稳，不少也与地震造成饱水固结土体的液化有关，这也是为什么在十分平缓的海底斜坡中会产生滑坡的重要原因之一。

我国岩质边坡工程实践中，为量化评价爆破的影响，根据经验采取降低计算结构面的抗剪强度的方法实施，f值降低15%～30%，c值降低20%～40%。

理论计算，降低的低值和高值分别相当于地震烈度Ⅷ度和Ⅸ度时造成的影响。

9.3 边坡稳定分析与评价
随着人类工程活动向更深层次发展，在经济建设过程中，遇到了大量的边坡工程，且规模越来越大，其重要程度也越高，有时会影响人类工程活动；并且人们更注重由于边坡失稳造成的地质灾害，故边坡稳定性研究一直是重中之重。

边坡稳定性分析与评价的目的，一是对与工程有关的天然边坡稳定性作出定性和定量评价；二是要为合理地设计人工边坡和边坡变形破坏的防治措施提供依据。

边坡稳定性分析评价的方法主要有：地质分析法(历史成因分析法)、力学计算法、工程地质类比法、过程机制分析法、理论体边坡已有的变形迹象，阐明其形成演变机制。

分析中要特别注意变形模式的转化标志，它往往是失稳的前兆。

边坡稳定性分析方法很多，简要归纳如下。

9.3.1 边坡稳定性分析方法简介
1）定性分析方法
主要是分析影响边坡稳定性的主要因素、失稳的力学机制、变形破坏的可能方式及工程的综合功能等，对边坡的成因及演化历史进行分析，以此评价边坡稳定状况及其可能发展趋势。

该方法的优点是综
合考虑影响边坡稳定性的因素，快速地对边坡的稳定性做出评价和预测。

常用的方法有：
（1）地质分析法（历史成因分析法）
根据边坡的地形地貌形态、地质条件和边坡变形破坏的基本规律，追溯边坡演变的全过程，预测边坡稳定性发展的总趋势及其破坏方式，从而对边坡的稳定性做出评价，对已发生过滑坡的边坡，则判断其能否复活或转化。

（2）工程地质类比法
其实质是把已有的自然边坡或人工边坡的研究设计经验应用到条件相似的新边坡的研究和人工边坡的研究设计中去。

需要对已有边坡进行详细的调查研究，全面分析工程地质因素的相似性和差异性，分析影响边坡变形发展的主导因素的相似性和差异性，同时，还应考虑工程的类别、等级及其对边坡的特定要求等。

它虽然是一种经验方法，但在边坡设计中，特别是在中小型工程的设计中是很通用的方法。

（3）图解法
图解法可以分为两类：
①用一定的曲线和偌谟图来表征边坡有关参数之间的定量关系，由此求出边坡稳定性系数，或已知稳定系数及其它参数（ 、c、r、结构面倾角、坡角、坡高）仅一个未知的情况下，求出稳定坡角或极限坡高。

这是力学计算的简化。

②利用图解求边坡变形破坏的边界条件，分析软弱结构面的组合关系，分析滑体的形态、滑动方向，评价边坡的稳定程度，为力学计算创造条件。

常用的为赤平极射投影分析法及实体比例投影法。

（4）边坡稳定专家系统
工程地质领域最早研制出的专家系统是用于地质勘察的专家系统Propecter，由斯坦福大学于70年代中期完成的。

另外，MIT在80 年代中期研制的测井资料咨询的专家系统也得到成功地应用。

在国内，许多单位正在进行研制，并取得了很多的成果。

专家系统使得一般工程技术人员在解决工程地质问题时能象有经验的专家给出比较正确的判断并做出结论，因此，专家系统的应用为工程地质的发展提供了一
条新思路。

2）定量评价方法
实质是一种半定量的方法，虽然评价结果表现为确定的数值，但最终判定仍依赖人为的判断。

目前，所有定量的计算方法都是基于定性方向之上。

（1）极限平衡法
极限平衡法在工程中应用最为广泛，这个方法以摩尔—库仑抗剪强度理论为基础，将滑坡体划分为若干条块，建立作用在这些条块上的力的平衡方程式，求解安全系数。

这个方法，没有象传统的弹、塑性力学那样引入应力－应变关系来求解本质上为静不定的问题，而是直接对某些多余未知量作假定，使得方程式的数量和未知数的数量相等，因而使问题变得静定可解。

根据边坡破坏的边界条件，应用力学分析的方法，对可能发生的滑动面，在各种荷载作用下进行理论计算和抗滑强度的力学分析。

通过反复计算和分析比较，对可能的滑动面给出稳定性系数。

刚体极限平衡分析方法很多，在处理上，各种条分法还在以下几个方面引入简化条件：（1）对滑裂面的形状作出假定，如假定滑裂面形状为折线、圆弧、对数螺旋线等；（2）放松静力平衡要求，求解过程中仅满足部分力和力矩的平衡要求；（3）对多余未知数的数值和分布形状做假定。

该方法比较直观、简单，对大多数边坡的评价结果比较令人满意。

该方法的关键在于对滑体的范围和滑面的形态进行分析，正确选用的滑面计算参数，正确地分析滑体的各种荷载。

基于该原理的方法很多，如条分法、圆弧法、Bishop法、Janbu 法、不平衡传递系数法等。

目前，刚体极限平衡方法已经从二维发展到目前的三维。

有关边坡稳定三维极限平衡方法，已有众多文献介绍研究成果。

下面就几种常用的方法做介绍。

（2）数值分析方法
主要是利用某种方法求出边坡的应力分布和变形情况，研究岩体中应力和应变的变化过程，求得各点上的局部稳定系数，由此判断边
坡的稳定性。

主要有以下几种：
①有限单元法（FEM）
该方法是目前应用最广泛的数值分析方法。

其解题步骤已经系统化，并形成了很多通用的计算机程序。

其优点是部分地考虑了边坡岩体的非均质、不连续介质特征，考虑了岩体的应力应变特征，因而可以避免将坡体视为刚体、过于简化边界条件的缺点，能够接近实际地从应力应变分析边坡的变形破坏机制，对了解边坡的应力分布及应变位移变化很有利。

其不足之处是：数据准备工作量大，原始数据易出错，不能保证整个区域内某些物理量的连续性；对解决无限性问题、应力集中问题等其精度比较差。

②边界单元法（BEM）
该方法只需对已知区的边界极限离散化，因此具有输入数据少的特点。

由于对边界极限离散，离散化的误差仅来源于边界，区域内的有关物理量是用精确的解析公式计算的，故边界元法的计算精度较高，在处理无限域方面有明显的优势。

其不足之处为：一般边界元法得到的线性方程组的关系矩阵是不对称矩阵，不便应用有限元中成熟的对稀疏对称矩阵的系列解法。

另外，边界元法在处理材料的非线性和严重不均匀的边坡问题方面，远不如有限元法。

③离散元法（DEM）
是由Cundall（1971）首先提出的。

该方法利用中心差分法解析动态松弛求解，为一种显式解法，不需要求解大型矩阵，计算比较简便，其基本特征在于允许各个离散块体发生平动、转动、甚至分离，弥补了有限元法或边界元法的介质连续和小变形的限制。

因此，该方法特别适合块裂介质的大变形及破坏问题的分析。

其缺点是计算时步需要很小，阻尼系数难以确定等。

离散单元法可以直观地反映岩体变化的应力场、位移场及速度场等各个参量的变化，可以模拟边坡失稳的全过程。

④块体理论（BT）
是由Goodman和Shi（1985）提出的，该方法利用拓扑学和群论。