0.   引言

滑带作为滑坡的重要结构要素，其强度变化对滑坡稳定性评价具有控制性影响。国内外有不少中大型水库岸坡因库水波动及滑带强度弱化而引起坡体变形及滑塌失稳，如三门峡水库、凤滩水库、二滩水库及意大利瓦依昂水库等。库水位升降作用导致的滑带强度弱化和产生的渗流驱动力对滑坡稳定性影响不容忽视，而如何量化分析渗流稳定性并提出精确的启滑判据成为库岸滑坡亟待解决的问题。在动水作用下，滑坡体渗流场变化、滑带强度弱化机制较复杂，给库区重大涉水滑坡的失稳预测预报带来了巨大挑战。因此，开展库水升降和滑带弱化作用下岸坡渗流场演变规律及启滑机制等问题的研究迫在眉睫，可为库区涉水滑坡灾害预测预报及工程建设提供可靠的科学依据和理论方法。

近年来，国内外诸多学者对降雨、水位升降及地下水等作用下的滑坡失稳机制开展了一系列的研究。谭淋耘等^[1]分析了滑坡变形在库水位7个不同变化过程下的响应特征和诱发机制。肖捷夫等^[2]制作了大型物理试验模型，分别模拟库水涨落和降雨工况，以及两者联合作用的工况，通过试验获得了坡体位移、孔隙水压力和土压力全时程曲线以及滑坡演化数字影像资料。李跃等^[3]以澜沧江流域大华桥水电站大华滑坡为例，针对其滑带土开展了土-水特征曲线试验，分析了非饱和渗透过程的特性及演化规律。He等 ^[4]以三峡库区某典型滑坡为研究对象，建立滑坡模型，研究了水位波动对滑坡的影响。Tu等 ^[5]采用基于有限元方法的模型，考虑入渗过程和滑坡稳定性，研究了该滑坡在水渗透作用下的活化作用。Cao等^[6]对排土场滑坡的触发启滑机制和动力过程进行分析评价，得知滑体抗剪强度与排土场土体的含水量、孔隙压力和级配关系密切。Chen等^[7]研究了滑坡启动阶段抗剪强度的衰减动力分布机理，推导了边坡的动、静平衡统一表达式，提出了用抗剪强度折减法计算稳定系数的方法。而滑坡稳定性主要受滑带土和软弱面强度的控制^[8]，滑带土强度参数的选取直接影响到边坡的稳定性分析和防护设计。Gao等^[9]提出了一种基于滑动带的土质边坡稳定性分析方法，并通过实例验证了新方法的有效性。Zhu等^[10]通过一系列直剪蠕变试验，研究了红层滑动土在不同竖向荷载和含水率条件下的蠕变特性。Liu等^[11]对原状滑带土的残余强度进行了研究，并进行了残余状态蠕变试验，选取不同的恒蠕变应力与残余强度之比，研究了恒蠕变应力与蠕变速率的关系，试图对滑体稳定进行预测。而Liao等^[12]认为强度参数的弱化规律可以用指数函数模型来描述和表征，并在考虑不同干湿循环次数的情况下，利用不平衡推力法和可靠度理论对浅、深滑动面的长期可靠度进行了定量分析。Zhao等^[13]研究了强降雨条件下不同因素对边坡稳定性的影响程度，结合数理统计方法，采用Geo-studio软件中Seep/W和Slope/W模块，耦合计算了边坡稳定系数。

当前研究大多数探讨分析了降雨及库水位波动下滑坡体的稳定性，对滑坡变形、水位变化及降雨等参数之间的相互影响关系及演变规律有了初步的认识，但对于其定量分析相对较少，对于滑带弱化机理及启滑判据的研究更为少见。鉴于以上，针对库水升降和滑带弱化作用下岸坡启滑机制与渗流稳定性问题，本文以三峡库区凉水井滑坡为研究对象，采用理论分析与数值模拟的研究手段，构建了滑带强度弱化模型、提出了渗流驱动的滑坡启滑判据，应用有限元程序分析了库水不同升降速率对滑坡稳定性的影响，查明了库水升降作用下滑坡渗流场演变规律，进而揭示了滑坡在渗流驱动下的启滑机制，为开展类似滑坡研究和工程防治提供了参考依据。

1.   滑坡启滑机制理论分析

1.1   滑带强度弱化模型

滑带作为控制滑坡稳定性的关键区域，其强度在滑坡体中最为薄弱，在三峡库区长时间的低水位疏干和高水位浸泡的干湿循环作用过程中，滑带土体抗剪强度及其他力学性质均有不同程度的弱化效应，而这些参数的弱化将影响滑坡整体稳定性。库水升降条件下水对滑带土体的弱化作用可分为软化作用和劣化作用^[14]，软化作用与渗透压本身无直接关系，而劣化作用可以直接破坏土体结构致使强度劣化，其与渗透压直接相关。弱化过程中存在一个与渗透压（P）、渗透时间（t）有关的弱化变量$D_{P,t}$，其定义为：

式中：$D_{P,t}$——弱化变量；

${\tau }_0$——滑带土体初始抗剪强度/kPa；

∆τ——抗剪强度差，$\mathrm{\Delta }\tau \text{ = }{\tau }_0-{\tau }_{P,t}$；

${\tau }_{P,t}$——渗透压P、渗透时间t下土体抗剪强度/kPa。

根据文献[14]分析，由于软化作用和劣化作用，弱化变量$D_{P,t}$可用软化参数和劣化参数表征，软化参数只与渗透时间t有关，表示为${\chi }_t$。劣化参数与渗透压P、渗透时间t有关，表示为${\omega }_{P,t}$。即：

式中：$D_{P,t}$——弱化变量；

${\chi }_t$——转化参数；

${\omega }_{P,t}$——劣化参数；

${\tau }_{0,t}$——渗透压0 kPa、渗透时间t下土体抗剪强度/kPa；

${\tau }_{P,t}$——渗透压P、渗透时间t下土体抗剪强度/kPa。

结合式（3）（4）可得：

1.2   渗流驱动启滑模型

库水位上涨时，坡体内会形成指向坡内的渗流场，渗流力指向坡内，具有提升坡体阻滑力的作用。库水位下降时，坡体前缘水头迅速降低，使坡体前缘和后缘形成水头差，坡体承受向外的渗流力。随着库水升降作用对滑坡土体的逐渐弱化，坡体前缘土体沉降变形，后缘发展张拉等类型裂缝，在渗流驱动作用下滑带土体发生压剪破坏而启滑，最终导致坡体渐进破坏失稳。

对于渗流驱动的滑坡压剪破坏理论分析模型如图1所示，图中h₁ 、h₂分别为初水位、终水位，$\alpha$为滑坡滑动面与水平的平均夹角，$L$为渗流路径水平投影长度，由Dupuit公式^[15]可得渗流区域岸坡任一截面单宽流量$q_x$为：

图  1  渗流驱动的滑坡压剪破坏分析模型

Figure  1.  Analysis model for seepage-driven landslides under pressure shear failure

下载: 全尺寸图片 幻灯片

由达西定律知：

因此：

可得任一截面水力梯度为：

则任一截面渗流力为：

进而可得考虑库水位降落速度（$V_{\mathrm{s}}$）、降落时间、库水位初始高度（$h_1$）、浸润线方程（$h_x$）以及岸坡滑面方程（$b_x$）的渗流总驱动力为^[16]：

式中：D——渗流总驱动力/kN；

$q_x$——任一截面单宽流量/（m³·s⁻¹）；

$k$——坡体渗透系数；

$h_1$、$h_2$——初水位、终水位/ m；

$L$——渗流路径水平投影长度/ m；

J——任一截面水力梯度；

A——垂直渗流方向的流经面积/m²；

$v_{\mathrm{s}}$——库水位降落速度/（m·d^-1）；

$h_x$——浸润线方程；

$b_x$——滑面关于$x$的函数。

根据传递系数法，在考虑重力、渗流驱动下滑带弱化后坡体压剪启滑判据为：

式中：$F_{\mathrm{s}}$——稳定性系数；

$R_i$——第i条块的抗滑力/（kN·m⁻¹）；

$T_i$——第i条块下滑力/（kN·m⁻¹）；

${\psi }_j$——第i条块剩余下滑力传递至i+1块段时的传递系数（ j=i时）；

$R_n$——第n条块的抗滑力/（kN·m⁻¹）；

$T_n$——第n条块下滑力/（kN·m⁻¹）；

$W_i$——第i条块的有效重量/（kN·m⁻¹），水位线以上取天然重度，水位线以下取浮重度；

$D_i$——第i条块的动水压力/（kN·m⁻¹）；

${\alpha }_i$——第i条块滑面倾角/（°）；

${\beta }_i$——第i块地下水位线平均倾角/（°）。

${\gamma }_{\mathrm{w}}$——水的重度，取10 kN/m³；

$l_i$——第i条块滑面长度/m；

$h_{\mathrm{w}}$——第i条块在水位面以下、库水位面以上范围内的高度/m；

${\phi }_i$——第i条块滑面内摩擦角标准值/（°）；

$c_i$——第i条块的滑面黏聚力/kPa；

根据式（12），当$F_{\mathrm{s}}<1$时，滑坡体在渗流驱动作用下发生压剪破坏而启滑；当$F_{\mathrm{s}}\text{ = }1$时，滑坡体处于临界状态；当$F_{\mathrm{s}}>1$时，滑坡体处于稳定状态。

2.   凉水井滑坡启滑机制验证

2.1   滑坡概况

凉水井滑坡位于重庆市云阳县水让村8组三峡库区段长江右岸的斜坡地带，滑坡整体平面形态呈“U”形，后部地形呈近似圈椅状，南高北低，中后部地形较陡，前部地形较缓，自然坡度30°～35°，滑坡全貌如图2所示。滑坡前缘高程约100 m，后缘高程约319.5 m，相对高差约221.5 m，平面纵向长度约434 m，横向宽358 m，面积约11.82×10⁴ m²，滑体平均厚度约34.5m，总体积约4.0779×10⁶ m³。

图  2  凉水井滑坡全貌

Figure  2.  The overall view of Liangshuijing landslide

下载: 全尺寸图片 幻灯片

凉水井滑坡滑体为滑坡堆积层，主要由含角砾粉质黏土、人工填土、含碎石、块石黏土以及砂、泥岩块石组成；滑床为中侏罗统沙溪庙组（J₂s）互层砂岩、泥岩；滑动带为滑坡堆积与下伏基岩接触带，滑动方向与现坡向基本一致。自三峡库区175 m试验性蓄水后，受水位上升影响，凉水井滑坡出现了不同程度的地表变形，主要表现为滑坡后缘地表拉张裂缝全部贯通，滑坡中部横向地表出现拉裂缝、中前部出现剪切裂缝以及两侧出现斜裂缝，典型剖面如图3所示。

图  3  凉水井滑坡典型地质剖面图

Figure  3.  Typical geological cross-section of the Liangshuijing landslide

下载: 全尺寸图片 幻灯片

2.2   变形特征分析

图4为凉水井滑坡滑面形态特征，滑带土厚度为3～5 cm，可见明显的擦痕，呈镜面，有光泽，局部有揉搓现象，混杂的少量泥岩碎石有挤压、角砾磨光现象，且后壁光滑，此为滑坡移动的特点。

图  4  凉水井滑坡滑面形态特征

Figure  4.  Morphological characteristics of the sliding surface of the Liangshuijing landslide

下载: 全尺寸图片 幻灯片

根据凉水井滑坡监测结果^[17]，由图5知，地表累积位移随时间推移逐渐增加，累计水平位移变形量平均为52.92 mm，最大已达93.88 mm。该滑坡变形从后部到前部、右侧到左侧、高处到低处存在轮流交接变形的特征。滑坡地表裂缝宽度累积与时间的变化关系如图6所示，图中所示裂缝宽度最大达到了28.35 mm，裂缝宽度累计变化量较大的监测点均分布在滑坡的右侧和后缘，说明滑坡主要变形部位位于滑坡的右侧和后缘，而位于滑坡左侧的监测点变化量较小，表明滑坡的整体变形方向是从滑坡的后缘右侧向滑坡前缘左侧推覆挤压。随着库水位的波动，坡体地表累积位移量逐渐增加，滑移变形较明显。

图  5  地表累积位移与时间的关系

Figure  5.  The relationship between cumulative surface displacement and time

下载: 全尺寸图片 幻灯片

图  6  地表裂缝宽度累积与时间的关系

Figure  6.  The relationship between cumulative surface crack width and time

下载: 全尺寸图片 幻灯片

2.3   滑带弱化分析

按照《三峡库区三期地质灾害防治工程地质勘察技术要求》，滑坡稳定性分级情况如下：当稳定性系数$F_{\mathrm{s}}<1$时，为不稳定状态；$1.00<F_{\mathrm{s}}⩽1.05$时，为欠稳定状态；$1.05<F_{\text{s}}⩽F_{\text{st}}$时，为基本稳定状态；$F_{\text{s}}＞F_{\text{st}}$时，为稳定状态。根据前述稳定性计算方法可得出图3所示剖面的稳定系数$F_{\mathrm{s}}$随库水位的变化，如图7（a）所示，可见坡体稳定性随库水位的下降由基本稳定状态过渡为欠稳定状态。根据勘查所取滑带土进行的室内抗剪试验结果分析，得出滑带抗剪强度与库水位变化关系，如图7（b）所示，凉水井滑坡随库水位周期波动，滑带抗剪强度逐渐弱化，库水位从175 m降到145 m，抗剪强度降低了19%，滑坡在渗流驱动下启动变形，总体呈现蠕滑变形特征。以上分析印证了在库水位升降作用下滑带土体抗剪强度的弱化效应。

图  7  凉水井滑坡稳定系数F_s-库水位（左）变化与滑带抗剪强度-库水位（右）变化

Figure  7.  Stability coefficient F_s - change of reservoir water level (left) and shear strength of sliding zone - change of reservoir water level (right) of Liangshuijing landslide

下载: 全尺寸图片 幻灯片

图8为滑带土抗剪参数反演关系，其根据滑坡体的特征及现状，由已知的稳定状态（取$F_{\mathrm{s}}=1.02$），反算出滑坡的综合抗剪强度参数取值。在c值不变时按线性差值法确定剖面稳定系数对应的φ值。根据反演计算结果，当黏聚力c > 11 kPa、内摩擦角φ=25°时，滑坡处于基本稳定状态；当11 kPa < c < 20 kPa、φ=24°时，滑坡处于欠稳定状态；当11 kPa < c < 19 kPa、φ=23°时，滑坡处于不稳定状态。根据前述渗流驱动力计算方法可得渗流驱动力与库水位降落速度关系如图9所示，可知渗流驱动力随库水位降落速度增大而增大，水流在滑坡土体中形成的流动管道被扩宽，孔隙与孔隙间逐渐被连通，土颗粒随水流流失，孔隙率逐渐增加，呈现渗流驱动力不断增大的现象。而库水位降落速度较小时，由于土体本身存在的滞后效应，导致渗流驱动力较小。

图  8  滑带土抗剪参数反演

Figure  8.  Inversion of shear resistance parameters for the sliding zone soil

下载: 全尺寸图片 幻灯片

图  9  渗流驱动力与库水位降落速度变化关系

Figure  9.  The relationship between seepage driving force and falling rate of reservoir water level

下载: 全尺寸图片 幻灯片

2.4   滑移机制验证

根据现场调查，结合该滑坡地形地貌、滑面形态、物质组成等特征分析，综合确定其形成原因为原岩质顺层老滑坡前缘在江水侵蚀、剥蚀以及河床切割作用下，前缘临空，产生滑移，并堆积于原地面上而形成该新滑坡。由于滑坡中前部及前部滑体主要为砂岩块石，其强度、刚度均较大，在滑动过程中，将原地面以下松散土层及强风化泥岩推移至新滑坡前缘，并在江水长期作用下产生库岸再造，最终形成现状地形地貌。且调查区可见基岩陡崖产生的崩坡积层以及陡崖上方的残坡积层在水、自身重力作用下堆积于新近形成的滑坡表层，最终形成现状滑坡滑体中前部较厚，后部逐渐变薄的趋势。而基于前文分析，滑坡体的变形和滑动取决于坡体的应力分布和强度特性，滑坡在库水位升降作用下，滑带土体抗剪强度及其它力学性质会有不同程度的弱化效应，使得滑坡前缘阻滑段有效抗滑力减小，而库水的长期侵蚀作用又使前缘滑体被剥蚀，最终导致滑坡稳定性降低，当下滑力超过抗滑力时，坡体将产生滑移破坏。

综上，将其灾变滑移演化机制分为以下四个阶段：

阶段一：岸坡长期受到长江水的冲刷和侵蚀，滑带抗剪强度逐渐弱化，前缘涉水区域逐渐被掏空而产生临空面，前缘发生了局部滑移，造成原始岸坡抗滑力减小。

阶段二：原始岸坡前缘发生局部变形后，中后缘坡体失去前部支撑作用，当滑带抗剪强度弱化至临界强度时，原始岸坡整体启滑，形成老滑坡。

阶段三：在原始岸坡发育为老滑坡后，由于库水位的波动及降雨影响，滑带抗剪强度损伤弱化进一步加强，老滑坡发生多次多级滑移，坡体后缘物质逐渐向中前部迁移堆积，形成了新滑坡的滑体堆积体。

阶段四：滑坡堆积体形成后，库水位的继续波动与降雨影响，新滑坡滑带抗剪强度等力学参数再一次被弱化损伤，滑坡变形加剧，滑坡后缘地表出现拉张裂缝，中前部出现剪切裂缝等，且全部贯通，形成了最终的凉水井滑坡。

凉水井滑坡形成机制及启滑过程如图10所示。

图  10  凉水井滑坡启动及演化机制

Figure  10.  Initiation and evolution mechanism of Liangshuijing landslide

下载: 全尺寸图片 幻灯片

3.   数值分析

3.1   渗流稳定性计算模型

选取如图3所示的滑坡典型剖面作为计算分析剖面，采用Geo-studio有限元软件中SEEP/W模块计算滑坡渗流场，采用SLOPE/W模块计算渗流驱动下的滑坡稳定性。根据凉水井滑坡地质模型，对于饱和与非饱和的渗流问题建立如图11所示的数值计算模型，模型节点数为2061，单元数为1955。将滑面设置为隔水边界，库水位175 m以上地表为定流量边界，以下为变水头边界。模拟库水位下降条件的渗流演变规律时，以175 m水位的稳态渗流作为初始渗流场，模拟库水位上升条件的渗流演变规律则以145 m水位的稳态渗流作为初始渗流场。

图  11  凉水井滑坡渗流分析有限元模型

Figure  11.  Finite element model for seepage analysis of Liangshuijing landslide

下载: 全尺寸图片 幻灯片

3.2   计算参数及工况

凉水井滑坡地勘报告及室内试验获取的滑坡物理力学参数如表1所示，三峡库区每年汛期为降雨量较大的5—9月，库水位波动特征总体是每年1—9月^[18-22]，为探究库水升降作用下滑坡渗流场演变规律，设置两种模拟工况，库水位从175 m以0.5 m/d、0.7 m/d、0.9 m/d、1.1 m/d的速率下降到145 m为工况一，汛期过后从145 m又以0.5 m/d、0.7 m/d、0.9 m/d、1.1 m/d的速率上升到175 m水位为工况二，将30 m的水位差作为库水位涨落值。

表  1  滑坡计算参数一览表

Table  1.  List of calculation parameters for Liangshuijing landslide

岩土体	弹性模量/GPa	泊松比	重度/（kN·m−3）			黏聚力/kPa			内摩擦角/（°）
			天然	饱和		天然	饱和		天然	饱和
滑床	24	0.2	25.1	25.3		5650	5420		47.04	43.80
滑带	3	0.4	23.0	23.8		19.17	14.41		25.03	24.29
滑体	3	0.4	23.0	23.8		21.48	19.29		34.83	27.68

下载: 导出CSV 

| 显示表格

3.3   计算结果分析

图12为库水升降作用下凉水井滑坡地下水渗流场演变规律，库水位上升到175 m后，坡体中前缘部位的总水头均进一步的升高，见图12（A2），在水位降落过程中，后缘有足够的水量补充，由后缘流向前缘，在与库水接触溢流区流速最大，见图12（A1）。由图12（B1）、（B2）知，库水位升降条件下，滑坡体前缘孔隙水压为正值，中后缘为负值，最大值为1675.8 kPa。由图12（C1）、（C2）所示的体积含水率变化知，库水位涨落过程中由于滞后效应，体积含水率较大区域均集中在滑坡体前缘部位。而土体饱和度变化趋势，上升过程呈上凹状，下降过程呈下凹状，饱和到非饱和的过渡区域较小，见图12（D1）、图12（D2）。无论库水位上升还是下降，坡体内的浸润面中部都呈现缓慢上升或下降，越远离库水面，水位影响越小，见图12（E1）、图12（E2）。

图  12  库水升降作用下凉水井滑坡地下水渗流场演变

Figure  12.  Evolution of groundwater seepage field of Liangshuijing landslide under the action of reservoir water level rise and fall

下载: 全尺寸图片 幻灯片

图13表示库水位升降速率分别为0.5 m/d、0.7 m/d、0.9 m/d、1.1 m/d的滑坡稳定性，可见稳定系数$F_{\mathrm{s}}$的变化随升降速率的增大而增大，库水上升过程中，滑坡内部地下水位的爬升滞后于库水位，形成内低外高的压力差，稳定性增加。而在库水位下降过程中，形成内高外低的压力差，渗流力成为向外推动坡体的驱动力，滑坡稳定性降低。图14为凉水井滑坡平均位移随降水时间变化曲线，可知随着库水位的下降，坡体地表累积平均位移量逐渐增加，滑移变形越明显。

下载: 全尺寸图片 幻灯片

图  13  库水升降作用下凉水井滑坡稳定性变化

Figure  13.  Stability change of Liangshuijing landslide under the action of reservoir water level rise and fall

下载: 全尺寸图片 幻灯片

图  14  凉水井滑坡平均位移随降水时间变化

Figure  14.  Variation of average displacement of Liangshuijing landslide with precipitation over time

下载: 全尺寸图片 幻灯片

图15所示为降水条件下滑面法向力变化，可知法向力峰值随降落时间的推移逐渐减小，滑面从非饱和向饱和过渡过程中，法向力弱化度逐渐增加，从175 m水位降到145 m，法向力最大降低了38.19%。图16为降水条件下滑面剪力变化，其变化规律与法向力变化基本一致，且剪力最大降低了22.20%。

图  15  降水条件下滑面法向力变化

Figure  15.  Variation of normal force on sliding surface under precipitation conditions

下载: 全尺寸图片 幻灯片

图  16  降水条件下滑面剪力变化

Figure  16.  Variation of shear force on sliding surface under precipitation conditions

下载: 全尺寸图片 幻灯片

图17为库水升降条件下滑坡孔隙水压变化，从非饱和向饱和过渡过程中，降水条件下土体孔隙水压增长幅度随降落时间的推移逐渐减小，根据有效应力原理，该区段有效法向力衰减将随降落时间的推移而逐渐增大，抗剪强度衰减规律亦如此，其变化曲线如图18所示。图中有效法向应力最大落差为168.64 kPa，抗剪强度最大落差为63.45 kPa。

图  17  库水升降条件下滑坡孔隙水压变化

Figure  17.  Changes in pore water pressure of the landslide under reservoir water level rise and fall conditions

下载: 全尺寸图片 幻灯片

图  18  降水条件下滑面有效法应力与抗剪强度变化

Figure  18.  Variation of effective normal stress and shear strength on sliding surface under precipitation conditions

下载: 全尺寸图片 幻灯片

4.   结论

（1）库水升降作用对滑带强度会产生弱化效应，且渗透压力与渗透时间的变化是滑带土强度弱化的关键因素，弱化到临界强度时，在渗流驱动作用下滑带土体会发生压剪破坏而启滑。

（2）库水升降过程中，坡体内孔隙水压力滞后性较明显，水位升降速率会影响坡体地下水响应时程，升降速率越大，渗流驱动力越大，滑坡稳定性变化越快，越趋近于渐进破坏。且库水位不论以多大的速率升降，坡体内浸润面中部都呈现缓慢上升或下降趋势变化。

（3）库水从175 m水位以0.5 m/d的速率降到145 m，凉水井滑坡滑面法向应力最大降低了38.19%，剪应力最大降低了22.20%，有效法向应力最大落差为168.64 kPa，抗剪强度最大落差为63.45 kPa。