引言

边坡是一个复杂难解的结构体，这种复杂难解源自于其自身材料特性和坡体内部结构等的不确定性.要想全面地认识、研究边坡在地震作用下的动力稳定性状态，需要对实际工程中的边坡进行总结分类，以便实现不同类型边坡地震稳定性的快速识别、评价和预测，同时为边坡工程的勘察评价、设计施工和运营维护等工作提供依据.

目前，边坡的分类方式主要为单因素分类法和多因素分类法两种.单因素分类法，一般按照边坡的成因、高度、材料、结构特征、使用年限和用途等单因素划分.例如孙玉科和徐义芳 (1959)根据坡体材料的地质特性将边坡分为岩质岩石边坡、半岩质岩石边坡、黏性土质边坡、砂性土质边坡和特殊边坡, 又根据边坡的工程类别将边坡分为路堑边坡、水坝边坡、露天矿边坡和建筑边坡等.多因素分类法的应用实例，例如《建筑边坡工程技术规范》(中华人民共和国住房和城乡建设部，2014) 中综合岩体完整程度、结构面结合程度、产状和自稳能力等4种因素来划分岩体；Deere等 (1967)提出岩石质量指标 (rock quality designation，简写为RQD) 分类法，将岩体质量划分为5个等级，从Ⅰ级到Ⅴ级岩体质量由极好递减到极差；Bieniawski (1974)提出了一种更为全面的综合分类方法，即岩体地质力学分类法 (rock mass rating，简写为RMR).但上述分类方法的对象均为地下工程岩体.

总结前人的研究可知，根据单一因素对边坡进行分类，虽可以更直接地认识边坡，但因影响边坡地震稳定性的因素很多，只考虑单一因素不足以分析其在地震作用下的稳定性，因此需要综合多种因素对边坡进行分类.本文将综合考虑边坡地质和材料力学性质等因素，对地震作用下的土质边坡进行划分，可为地震作用下土质边坡稳定性评价分析奠定基础.

1.   划分条件

对土质边坡进行划分，首先要确定划分指标.前人关于地震滑坡影响因素的研究已有很多.例如Keefer (2000)对1989年Loma Prieta地震引发的滑坡进行了调查分析，认为地震引发的滑坡分布主要与震级、震中距、边坡岩性和坡度等因素有关. 2008年5月12日汶川M_S8.0大地震和2013年芦山M_S7.0地震均引发了大量的滑坡灾害 (黄润秋等，2008；裴向军，黄润秋，2013).许多地震触发的滑坡分布的研究 (黄润秋，李为乐，2008；许冲等，2009；王秀英等, 2009, 2010；李果等，2011；崔鹏等，2013；裴向军，黄润秋，2013) 表明，地形坡度和高程、水系距离、地震因素、地震断裂、坡体结构和岩性等均应作为边坡地震稳定性的主要控制因素.

通过上述研究分析可以看出，地震诱发滑坡的影响因素主要为地形、水文、岩性和地震参数等4个方面，具体表现为坡度、坡高、斜坡结构类型、降水、水系距离、黏聚力、内摩擦角和地震动参数等8个因素.由于这些影响因素均具有一定的不确定性和模糊性 (袁仁茂等，1999；黄润秋，2007)，若将所有影响因素作为划分边坡类型的依据，则划分结果势必杂乱，且操作性和实用性不强，因此须计算上述影响因素的权重，选取权重值较大的主导因素作为划分指标.为了更准确、更方便地评价地震作用下土质边坡的稳定性，本文基于以下原则对边坡进行划分：①综合性.影响地震滑坡的因素有多种，结合地质因素和力学参数，更全面地划分土质边坡类型；②操作性.将影响因素中的主导因素作为划分依据，以实现实际工程运用中的可操作性和高效性；③实用性.在实际工程的基础上进行分类，通过文中得到的土质边坡类型表实现地震作用下边坡稳定性的快速识别、评价和预测.

2.   地震作用下边坡稳定性影响因素的权重计算

2.1   层次分析法

层次分析法 (analytical hierarchy process，简写为AHP)，是美国运筹学家Saaty和Wong (1983)于20世纪70年代提出的一种对复杂的定性问题进行定量描述的简便、灵活而又实用的多目标评价决策方法，该方法把复杂问题中的多个影响因素进行两两比较，而后建立结构模型和数学模型来计算分析各元素的相对重要性次序和影响权重，为无法完全用定量或定性方法解决的问题提供了一种行之有效的策略.近年来层次分析法已被广泛运用于地震滑坡的危险性评估中 (樊晓一等，2004；李萍，2006；陈晓利等, 2008, 2012)，并取得了良好的效果.

本文将利用层次分析法，研究地震作用下和非地震作用下影响土质斜坡稳定性的关键因素及其权重计算等问题.

2.2   各影响因素的权重计算

1) 建立递阶层次结构模型.将第一节选取出的坡度、坡高、斜坡结构类型、降水、水系距离、黏聚力、内摩擦角和地震动参数等8个因素作为研究对象，建立如图 1的层次分析结构模型.

图  1  层次分析结构模型

Figure  1.  AHP structure model

下载: 全尺寸图片 幻灯片

2) 构造判断矩阵.建立层次分析结构模型后，需要构造判断矩阵，用以对模型第二层 (指标层) 中各因素对上一层因素的相对重要性进行量化，量化的判断标度及含义如表 1所示，若因素i与因素j的比较值为a_ij，则因素j与因素i的比较值a_ji=1/a_ij.

表  1  判断矩阵的标度及含义

Table  1.  Scale and implication of judgment matrix

标度	含义
1	两因素相比，重要性相同
3	两因素相比，一个因素比另一个因素稍微重要
5	两因素相比，一个因素比另一个因素明显重要
7	两因素相比，一个因素比另一个因素强烈重要
9	两因素相比，一个因素比另一个因素极端重要
2，4，6，8	上述相邻判断标度的中间值

下载: 导出CSV 

| 显示表格

根据目前对地震滑坡影响因素的认识和相关专家的经验，结合表 1，对指标层中8个因素之间的相对重要性进行判断、取值，并将取值结果组成判断矩阵A=(a_ij)_8×8.

为了度量该判断矩阵A的一致性和可信度，需要对判断矩阵进行一致性检验，为此引入两个指标，即一致性指标CI和平均随机一致性指标RI.

一致性指标CI的计算公式为

式中，λ_max为判断矩阵的最大特征根，m为判断矩阵的元素个数，此时m=8.

平均随机一致性指标RI，是通过随机构造500个判断矩阵A₁，A₂，A₃，…，A₅₀₀作为参照对比，分别计算各矩阵对应的一致性指标CI₁，CI₂，CI₃，…，CI₅₀₀，并取其平均数得到平均随机一致性指标：

1—9标度下判断矩阵的平均随机一致性指标RI值如表 2所示 (Saaty，Wong，1983).

表  2  判断矩阵的平均随机一致性指标RI

Table  2.  Average random consistency scale RI of judgment matrix

判断矩阵标度		1	2	3	4	5	6	7	8	9
RI		0	0	0.58	0.90	1.12	1.24	1.32	1.41	1.45

下载: 导出CSV 

| 显示表格

当判断矩阵的阶数大于2时，其一致性指标CI和平均随机一致性指标RI的比值称为随机一致性比率，记为CR=CI/RI.当CR＜0.1时，判断矩阵才具有较好的一致性，即各因素的相对重要性取值的排序是合理的，否则需要重新排列判断矩阵的取值.本文通过多次调整取值并试算，最终获得上述8个因素的合理的判断矩阵为

具体结果列于表 3.

表  3  地震滑坡影响因素判断矩阵

Table  3.  Judgment matrix of seismic landslide influence factors

	黏聚力	降水	斜坡结构类型	坡度	坡高	内摩擦角	水系距离	地震动参数
黏聚力	1	3	5	1/2	4	2	6	1/3
降水	1/3	1	3	1/4	2	1/2	4	1/5
斜坡结构类型	1/5	1/3	1	1/6	1/2	1/4	2	1/7
坡度	2	4	6	1	5	3	7	1/2
坡高	1/4	1/2	2	1/5	1	1/3	3	1/6
内摩擦角	1/2	2	4	1/3	3	1	5	1/4
水系距离	1/6	1/4	1/2	1/7	1/3	1/5	1	1/8
地震动参数	3	5	7	2	6	4	8	1

下载: 导出CSV 

| 显示表格

计算上述判断矩阵的最大特征值，得到λ_max=8.2883，由一致性指标计算公式，即式 (1) 求得CI=(8.2883－8)/(8－1)=0.0412，根据表 2得到平均随机一致性指标RI为1.41，则随机一致性比率CR=0.0412/1.41=0.0292＜0.1，说明该判断矩阵具有较好的一致性，各因素的相对重要性排序是合理的.

之后求出与判断矩阵A的最大特征根λ_max=8.2883对应的特征向量W＝(0.3449，0.1556，0.0717，0.5059，0.1046，0.2323，0.0517，0.7267)，并对其进行归一化处理，得到的向量W′＝(0.1572，0.0709，0.0327，0.2306，0.0477，0.1059，0.0236，0.3313) 即为各影响因素的权重值 (表 4).由表 4可知，对土质边坡稳定性影响最大的因素为地震动参数，之后依次为坡度、黏聚力、内摩擦角、降水、坡高、结构类型和水系距离.

表  4  地震作用下土质边坡稳定性影响因素的权重

Table  4.  Weight of each influence factor of soil slop stability under seismic action

影响因素		黏聚力	降水	斜坡结构类型	坡度	坡高	内摩擦角	水系距离	地震动参数
权重值		0.1572	0.0709	0.0327	0.2306	0.0477	0.1059	0.0236	0.3313

下载: 导出CSV 

| 显示表格

根据计算得到的各因素的权重值，本文选取权重值较大的前4个因素，即地震动参数、坡度、黏聚力和内摩擦角作为划分土质边坡类型的指标.

3.   非地震作用下边坡稳定性影响因素的权重计算

在不考虑地震作用的情况下，利用层次分析法计算坡度、坡高、斜坡结构类型、降水、水系距离、黏聚力和内摩擦角等7个影响因素的权重.根据各影响因素相互作用的大小作多次取值调整并试算，最终获得合理的判断矩阵为

具体结果列于表 5.

表  5  诱发滑坡的影响因素判断矩阵

Table  5.  Judgment matrix of induced landslide influence factors

影响因素	黏聚力	降水	坡高	坡度	水系距离	斜坡结构类型	内摩擦角
黏聚力	1	1/3	5	1/2	4	2	3
降水	3	1	7	2	6	4	5
坡高	1/5	1/7	1	1/6	1/2	1/4	1/3
坡度	2	1/2	6	1	5	3	4
水系距离	1/4	1/6	2	1/5	1	1/3	1/2
斜坡结构类型	1/2	1/4	4	1/3	3	1	2
内摩擦角	1/3	1/5	3	1/4	2	1/2	1

下载: 导出CSV 

| 显示表格

计算上述判断矩阵的最大特征值，得到λ_max=7.1955，对于7阶矩阵，m=7，则由一致性指标计算公式求得CI=(7.1955－7)/(7－1)=0.0326，根据表 2得到平均随机一致性指标RI为1.32，则随机一致性比率CR=0.0326/1.32=0.0247＜0.1，说明该判断矩阵具有较好的一致性，各因素的相对重要性排序是合理的.

之后进一步求得判断矩阵A的最大特征根λ_max=7.1955所对应的特征向量W＝(-0.3333，-0.7444，-0.0655，-0.5041，-0.0941，-0.2177，-0.1420)，并对其进行归一化处理，得到向量W′＝(0.1586，0.3543，0.0312，0.2399，0.0448，0.1036，0.0676)，即为各影响因素的权重值 (表 6).

表  6  非地震作用下土质边坡稳定性影响因素的权重

Table  6.  Weight of each influence factor of soil slop stability under non-seismic action

影响因素		黏聚力	降水	坡高	坡度	水系距离	斜坡结构类型	内摩擦角
权重值		0.1586	0.3543	0.0312	0.2399	0.0448	0.1036	0.0676

下载: 导出CSV 

| 显示表格

根据计算得到的各因素权重值，对比分析表 4与表 6可知，边坡稳定性的影响因素在有无地震作用下所得的权重是有所差别的.在非地震作用下，降水对土质边坡稳定性的影响最大，之后依次为坡度、黏聚力、斜坡结构类型、内摩擦角、水系距离和坡高; 与地震作用下各影响因素的权重相比，降水、斜坡结构类型和水系距离的影响权重显著增大.相关统计显示，自20世纪80年代以来发生在我国大陆范围内的滑坡灾害，约一半是由降雨引发的 (黄润秋，2007)，这也很好地印证了本文的结论.

4.   土质边坡类型

4.1   坡体材料性质

土质边坡材料具有多样性，而不同土质的力学性质不同，在地震作用下表现出的稳定性亦不同.根据土体的工程分类 (刘祖典，1997；唐大雄等，1999；孔宪立，石振明，2001；姜德义，王国栋，2003；李萍，2006；《工程地质手册》编委会，2011)，将土质边坡材料的内摩擦角和黏聚力划分为以下4段：

1) 内摩擦角介于15°—24°，黏聚力介于30—65 kPa.典型代表为黏土边坡.以黏性土为主的边坡称为黏土边坡.黏土以细密颗粒为主，含有较多黏粒和亲水性黏土矿物，是修筑堤坝的常用土料.但由于黏土的生成时间及空间环境的差异，其物理力学特性等也有较大的差别，对边坡稳定性的影响也不同.尤其是稠度，如固态或者硬塑态的黏土, 抗剪强度较高，其组成的边坡一般具有较好的稳定性.

2) 内摩擦角介于40°—76°，黏聚力介于20—33 kPa.典型代表为黄土边坡.黄土，顾名思义，是颜色呈黄色的土，天然状态下含水量较少，对湿度变化很敏感.黄土边坡在我国分布较广，甘、陕、晋地区尤为常见，宁夏、河北和河南等部分地区也有分布.在地震和降雨的作用下，黄土边坡极易发生失稳破坏，造成严重的人员伤亡与经济损失.

3) 内摩擦角小于15°，黏聚力小于20 kPa.典型代表为软土边坡.软土边坡是指坡体主要由淤化、泥变、淤泥实土及其它抗剪强度极低的土组成的边坡.软土边坡主要分布于我国东海、黄海、渤海和南海等沿海地区，内陆地区也有少量分布，这类土体的抗剪强度较低.

4) 内摩擦角介于24°—42°，黏聚力近似为零.典型代表为砂土边坡.砂土边坡是指主要由砂或砂性土组成的边坡.砂土的特性与砂粒的大小和密度有关.一般情况下，砂土的透水性较强.因砂土中黏粒含量很少，其黏聚力也很小，几乎为零.在地震作用下，砂土边坡多易产生液化而导致滑坡.

4.2   边坡坡度

坡度是反映边坡面倾斜程度的一个物理量，坡度对滑坡的发育具有重要的影响 (樊晓一等，2004).边坡坡度决定了坡体有效临空面的大小，而有效临空面又为滑坡的发生提供了必要的空间条件.在地震作用下，滑坡等级的分布规律有所不同，但是大量研究统计显示，地震滑坡多发生在20°—50°的坡度范围内 (孙崇绍，蔡红卫，1997；黄润秋，2007；许冲等，2009).本文参考姜德义和王国栋 (2003)一文，将坡度划分为4个等级, 详见表 7.

表  7  边坡坡度等级类型

Table  7.  Slop gradient type

等级	边坡名称	坡度/°
Ⅰ	缓坡	＜20
Ⅱ	缓陡坡	20—35
Ⅲ	陡坡	35—50
Ⅳ	急坡	＞50

下载: 导出CSV 

| 显示表格

4.3   土质边坡类型划分

由于地震动参数的影响因素较多，使得探讨变得较为复杂且不具操作性.作者此前已对地震动参数的土质边坡响应特性进行过讨论 (张江伟，2015)，研究结果表明，地震动的峰值加速度、频谱和持时对地震边坡稳定性均有不同程度的影响，其中峰值加速度和频谱的影响较为显著.地震动中幅值、频谱和持时的不同组合对边坡稳定性的影响程度也不同.本文考虑的地震作用主要体现在影响因素的选取，比如降水等因素在非地震作用下对边坡稳定性的影响作用显著，但在地震作用下则相反.根据地震作用下边坡影响因素及其权重的研究结果可知，以坡度、黏聚力和内摩擦角这3个指标进行土质边坡分类比较合适，其分类结果见表 8.首先将土质边坡按照土体的内摩擦角和黏聚力划分为4个大类，然后根据坡度等级将土质边坡进一步细分.

表  8  土质边坡分类

Table  8.  Soil slope classification

类型	内摩擦角/°	黏聚力/kPa	坡度				参考材料
			缓坡/°	缓陡坡/°	陡坡/°	急坡/°
Ⅰ	15—24	30—65	＜20	20—35	35—50	＞50	黏土
Ⅱ	40—76	20—33	＜20	20—35	35—50	＞50	黄土
Ⅲ	＜15	＜20	＜20	20—35	35—50	＞50	软土
Ⅳ	24—42	≈0	＜20	20—35	35—50		砂土

下载: 导出CSV 

| 显示表格

5.   讨论与结论

影响土质边坡地震稳定性的因素有很多，在研究时需要对其进行筛选，利用主导因素来研究土质边坡的地震稳定性，这样更具操作性和效率性.本文在对地震作用下土质边坡类型研究分析的过程中，得到以下几点结论：

1) 通过对边坡地震稳定性影响因素的研究分析，认为影响地震诱发滑坡的主要因素涉及地形、水文、岩性和地震参数等4个方面，具体表现为坡度、坡高、斜坡结构类型、降水、水系距离、黏聚力、内摩擦角和地震动参数等8个主要因素.

2) 利用层次分析法对上述8个影响因素建立的层次结构模型进行分析计算，得到各因素在有、无地震作用下土质边坡稳定性的影响权重.结果表明：在地震作用下，影响土坡稳定性的主要因素为地震动参数，之后依次为坡度、黏聚力、内摩擦角、降水、坡高、斜坡结构类型和水系距离；非地震作用下，降水对土质边坡稳定性的影响作用最大，之后依次为坡度、黏聚力、斜坡结构类型、内摩擦角、水系距离和坡高.

3) 以地震作用下影响权重较高的坡度、黏聚力和内摩擦角等3个因素为指标，建立地震作用下土质边坡类型的划分方案，将土质边坡划分为4类，每个类型中又根据坡度的不同划分为缓坡、缓陡坡、陡坡和急坡等4个亚类.

本文的划分结果可为土坡地震稳定性的数值模拟分析及危险性评价提供基础和依据.值得注意的是，对于上述各类型土坡的易滑性有待进一步研究.