模糊数学评判和数值模拟相结合的土质边坡稳定性综合评价

王崇敬，张 龙，刘国伟

（1.山西工程技术学院地球科学与工程系，山西 阳泉 045000；2.辽宁工程技术大学矿业学院，辽宁 阜新 123000；3.中国建筑第六工程局有限公司，天津 300450；4.中国地质环境监测院，北京 100081）

0 引言

滑坡是威胁山西地区人民生命和财产安全的主要地质灾害之一，广布的黄土层为滑坡灾害的发生提供了物质基础，尤其是汾河流域多为非自重的湿陷性黄土，具有显著的结构疏松与压缩性高的性质[1-2]，而土质边坡以易发生、危害强、难预测等特点危害性更为突出[3-4]。滑坡现象多为内部因素和外部因素共同作用触发，内部因素包括坡体本身的性质弱化以及较差的地质构造环境等，降雨、工程开挖等外部因素是诱发土质边坡失稳的直接原因[5-6]。

众多学者围绕土质边坡稳定性问题开展了一系列卓有成效的研究，尤其是20 世纪70年代后，有限元法借助计算机开始被广泛应用于边坡稳定分析。然而，在实际应用中，有限元法仍存在一些局限性，如较难确定边坡的初始应力状态以及把握边坡临近破坏时的弹塑性本构关系等。极限平衡法逐渐形成了垂直条分法和滑移线法两个独立的分支[7-8]。我国学者在滑动机理及监测方面也取得了一系列研究成果，彭建兵等[9]研究了黄土高原滑坡灾害形成的动力学机制问题，认为区域构造应力、边坡构造应力及易灾特性是黄土滑坡的“三大元凶”；董震等[1]聚焦山西省黄土地质灾害问题，认为黄土地质灾害的发生受其特有的物质成分和工程性质控制，并具有地域性分布的特点；支泽民等[10]研究发现，边坡稳定性影响因子相互作用会促进滑坡发育；蔡欣育等[11]开展了不同降雨类型对边坡稳定性的影响研究，重点分析了不同降雨类型条件下边坡孔隙水压力的分布特征；侯天顺等[12]也进行了人工降雨作用下的滑坡试验，模拟长期连续降雨边坡稳定系数的变化规律；马蓓青等[13]探讨了降雨条件下边坡裂缝发育规律及其对稳定性的影响；苗青等[14]通过LASSO 算法建立了能够预测滑坡危险区域的模型。

然而，在工程扰动下土质边坡稳定性愈加复杂，很难用一种模式概括其滑动机理与影响因素。调查表明，多数土质边坡的滑坡与降雨有直接或间接的关系。在近年极端天气频发的背景下，探讨工程场地客观地质条件下极端降雨对土质边坡稳定性的影响是十分必要的。本文在山西太原古交市马兰镇营立村拟建水泥厂地质调研基础上，基于模糊数学理论获取该土质边坡稳定的主要影响因素，确定权值并计算稳定系数，结合FLAC3D 软件模拟得到不同降雨强度下边坡稳定系数变化特征，进而得到现有条件下维持该边坡稳定所能承受的极端降雨条件，从而综合判断边坡的稳定状态。

1 工程地质条件

1.1 地形与岩性特征

研究区位于山西省太原市西部山区的古交市马兰镇营立村水泥厂，属于地形起伏较大的山间河谷地貌单元区，拟建场地地形起伏较大，各钻孔孔口标高介于1 182.43～1 194.64 m，相对高差为12.21 m。区内地层发育不完全，据勘察及钻孔揭露的地层情况，勘察深度范围内地层主要由全新统粉土、圆砾与中奥陶统石灰岩组成（图1）。粉土（Qh）呈褐黄色，稍密，稍湿，土质较松散；场地局部出露，可见零星煤屑、云母、氧化铁及植物根系，摇振反应中等，无光泽反应，干强度和韧性低。圆砾（Qhal+pl）为杂色，中密，稍湿，母岩成分主要为灰岩和砂岩，呈椭圆状，磨圆度较好，一般粒径为3～10 cm，最大粒径为18 cm，以粉土、中砂充填其中，含大量漂石，中奥陶统马家沟组石灰岩（O2m）为灰白色，强风化-中风化，岩体节理裂隙较发育，岩芯较完整，呈短柱状或柱状，一般在10.0～30.0 cm，敲击声较清脆，不易击碎，为较硬岩。

图1 拟建场地地质条件Fig.1 Geological conditions of the planned construction site

1.2 水文地质与气候

场地内水文地质条件简单，在勘察深度内未发现稳定地下水位。研究区属于典型的季风气候，当地历年降雨资料显示年平均降雨量为426.1 mm，降雨主要集中在7—8 月，其降雨量约在60～90 mm，但近年极端天气明显增多[14]，多次出现了日降水量超过100 mm 的大暴雨，极端天气频现给区域内土质边坡稳定性带来了挑战，也增加了研究的复杂性和难度。

2 边坡稳定性模糊数学分析

2.1 边坡稳定性评判模型

边坡稳定性受到内因和外因的双重影响。在现场调研基础上，采用层次分析法确定了影响该水泥厂边坡稳定性的因素及其分级标准（表1），确定了地层岩性、地貌特征、风化作用、降雨及人类活动[15]等 5 项影响因素，构建了边坡稳定性的模糊综合评价模型

表1 稳定性影响因素与量化取值表Table 1 Stability influencing factors and quantified values

根据构建的边坡稳定性模糊综合评价模型[16]，建立评价对象的关联因素集，表示为：

根据边坡实际情况，结合半定量取值法，确定不同因素权重集，表示为：

通过评价对象可能出现的结果建立备择集，表示为：

根据单个影响因素，联立式（2）与（3），确定各离散值对备择集的隶属程度，得出评判矩阵，表示为：

根据式（2）与（4），依据权重集与影响因素评判关系矩阵的乘积可以有效地表达所有因素对判定对象的综合影响指数，表示为：

2.2 边坡稳定影响因素关系矩阵

基于地质资料及调查结果，对各影响因素在不同稳定系数下的影响程度进行评价。基于所构建边坡稳定性评价模型的降雨量、人类活动、地貌特征、地层岩性与风化程度等5 项因素建立关系矩阵，矩阵主对角线上数据表示各影响因素在不同稳定系数下的影响程度。矩阵中各行表示该因素对其它因素在不同稳定系数下的影响程度，设该因素为A。矩阵中各列表示其它因素对其在相同稳定系数下的影响程度，设其为B。不同因素对不同稳定系数的影响具有差异性，因而选用编码方式加以区分。按影响程度由弱到强分别赋值为0，1，2，3，4。其中，Aa为关系矩阵中不同稳定系数的稳定性影响程度，该值越大表示影响越大；Bb为关系矩阵中相同稳定系数的稳定性重要程度，该值越大表示越重要。

基于以上原则建立了关系矩阵，计算公式如下：

式中：Aa——稳定性影响程度；

Bb——稳定性重要程度；

Hi——影响因素的权重（表2）。

表2 营立村水泥厂切坡稳定性影响因素赋值表Table 2 Impact factors of artificial slope cutting stability in Yingli Village quarry

根据地质调查结果及相关地质资料与规范，由式（6）—（8）计算，最终确定各影响因素的权重，降雨为26%、人类活动为22%、地貌特征为17%、地层岩性为19%、风化程度为16%。确定权重集为：

一般土质边坡的稳定系数区间为[1.20,2.00]，区间为0.25 均匀分布，建立备择集表示如下：

由此，计算得出综合评价结果为：

最后，对综合评价结果加权平均求得边坡稳定系数，计算结果如下：

式中：Fs——边坡稳定系数。

综上可知，基于模糊数学求得Fs=1.429>1.15，根据《滑坡防治工程勘察规范》（GBT 32864—2016）[17]第13.3.4 条滑坡稳定状态划分标准，该土质边坡天然条件下处于稳定状态。

3 边坡稳定性FLAC3D 分析

3.1 模型构建与物理力学参数获取

基于模糊数学分析已知降雨是营立村水泥厂土质边坡稳定性的主要影响因素，为进一步量化分析该土质边坡稳定性，使用FLAC3D 软件模拟不同降雨条件下该边坡的稳定性演变特征[18]。模型主要依据边坡的坡高、坡角等几何尺寸建立，由于边坡开挖为阶梯状，z轴分别为8.5 m、13.5 m、25.3 m，x轴长度为42 m，y轴长度为40 m。模型初始边界条件为：流体边界设置边坡表层为渗透面，坡脚和边坡底部为渗漏面；位移边界固定边坡底面，约束模型两侧的y轴位移。此外，为便于分析坡体在不同降雨条件下各层岩土体弹性模量、容量、剪应力、摩擦角等物理力学参数的变化情况[19]，本次模拟在所建模型每级台阶坡脚处在对应台阶内部提取位移监测点与应力监测点（图2）。

岩土体物理力学参数基于地质调查结果及土工试验确定（表3），所建边坡模型采用Mohr-Coulomb 准则和强度折减法求取稳定系数[20]。潜在滑动面主要计算原理如下：

表3 边坡岩土体物理参数指标Table 3 Physical parameter indicators of slope rock and soil

式中：c、c′——边坡土体初始状态和临界状态下的黏聚力/kPa；

φ、φ′——边坡土体初始状态和临界状态下的内摩擦角/（°）；

α——最危险滑动面与最小主应力夹角/（°）；

δ1、δ3——土体单元的第1 和第3 主应力/kPa。

3.2 天然状态下边坡稳定模拟分析

基于表3 的岩土体力学参数采用FLAC3D 数值模拟得到营立村水泥厂土质边坡在天然状态下边坡位移和剪应力空间分布图（图3）。图3（a）表明，天然状态下边坡的位移集中在边坡顶部区域，最大位移为29.3 cm，自坡顶向坡脚位移量具有减小趋势，直至边坡中下部无位移量。图3（b）表明，边坡内部应力变化最为显著，坡体在自重及岩性影响下剪应力集中分布于石灰岩及其上覆第四系沉积层界面附近，剪切力增量向边坡底部集中，最大剪切力累积在坡脚位置，说明在没有断裂的情况下马家沟组灰岩与其上覆第四系沉积层界面存在明显的渗透系数差，是最主要的结构面。模拟分析结果显示，此状态下边坡Fs=1.426，与模糊数学分析法得出的Fs=1.429 几乎一致，综合判定该边坡天然状态下处于稳定状态。

图3 天然状态下边坡位移和剪应力空间分布图Fig.3 Spatial distribution of slope displacement and shear stress in natural state

3.3 降雨状态下边坡稳定性变化

降雨会引起土体黏聚力和内摩擦角改变，进而影响边坡稳定性[21-23]。为了简化研究，将均匀雨型日降雨量设为基础分析条件。FLAC3D 模拟不同降雨强度条件下该边坡的稳定性，结果显示，随着日降水量增加稳定系数小幅升高后下降，整体呈现单峰型曲线形态（图4）。由图4 可知，较小的降雨强度对边坡渗流场影响有限，基本无法引起地下水位变化，而雨水渗透引起边坡粉土的黏聚力增大超过了内摩擦角减小的效应，因此边坡稳定性出现了小幅升高。随着日降雨量继续增加，土体黏聚力与内摩擦角均为下降趋势，导致稳定系数快速下降，当降雨量增长至75 mm/d 时Fs=1.15，边坡处于稳定状态；当降雨量增长至120 mm/d 时Fs=1.05，边坡处于基本稳定状态；当降雨量增长至150 mm/d 时，稳定系数为1.004，此时边坡已处于极限平衡的欠稳定状态，存在发生滑坡的潜在风险。

图4 日降雨量与稳定系数拟合图Fig.4 Fitted diagram of daily rainfall and safety and stability coefficient

3.4 极端降雨条件下边坡稳定状态分析

150 mm/d 的极端降雨引起了潜在滑动带由坡体表面逐渐向坡体深部转移，直至滑动带沟通了岩土界面，计算模型中出现了一条明显的塑性剪切带，见图5（a）。降雨约为0.5 h 时，雨水快速入渗引起土体自重增大，变形区域主要集中在表土层，局部具有轻微向外隆起的趋势，Y2、Y3 监测点应力随之上升，分别为178.2 kN、132.4 kN；降雨约为3 h 时，因入渗土体中水不能及时排出，孔隙水压力抵消一部分土体压力，内部的应力则呈现减小的趋势，Y2、Y3 监测点应力分别为123.5 kN、53.5 kN；此后，Y3 监测点应力基本不变，Y2 监测点经历快速下降（3～9 h）后进入缓慢下降阶段，约至20 h 时，Y2 监测点应力稳定为59.8 kN。处于坡脚岩土界面附近的Y1 监测点应力随降雨历时整体呈上升趋势，变化趋势大致呈三个阶段，6 h 内为应力较快上升阶段，6～16 h 为应力缓慢上升阶段，16 h 后应力增长达到最大限度，最大应力增量约为78 kN，见图5（b）。随降雨历时增加，边坡顶面浅表层首先产生了拉裂缝，裂缝产生的位移量由顶部表层向内部减小，即相同降雨历时条件下，监测点位移量W4>W3>W2>W1。降雨约18 h 时，各监测点位移量达到最大值，见图5（c），边坡已基本处于极限平衡状态。

图5 极限平衡状态边坡特征Fig.5 Characteristics of slope in limit equilibrium state

4 结论

综上可知，降雨对土质边坡稳定性有着显著的影响，雨水入渗导致土体自重增大、强度减小，使边坡内部发生局部剪切作用。随着降雨强度增大，边坡岩土层剪切面贯通形成完整的滑动面，最终导致土质边坡失稳。

（1）基于地质资料及野外调研基础上，构建了边坡稳定性模糊综合评价模型，结合关系矩阵确定降雨为营立村水泥厂边坡稳定性最主要的影响因素（26%），并求得Fs=1.429>1.15，表明该土质边坡天然条件下处于稳定状态。

（2）使用FLAC3D 模拟分析结果显示天然状态下边坡Fs=1.426；当降雨量75～120 mm/d 时，Fs由1.15 降至1.05；降雨量120～150 mm/d 时，Fs由1.05 降至1.004，此时边坡已处于极限平衡的欠稳定状态，存在发生滑坡的潜在风险。

（3）模糊数学与FLAC3D 评价营立村水泥厂土质边坡天然状态下稳定性结果几乎一致，联合运用两种方法相互佐证能够提高评价准确性，使结果更符合边坡实际情况。