0.   引言

我国西南岩溶山区地质环境脆弱，多形成“上硬下软、上陡下缓”的孕灾地貌，加上人类工程活动强烈，极端降雨频发，是我国高位远程崩滑灾害的高易发区^[1]。这种高位远程滑坡不仅监测预警难、突发性高，并且动力致灾过程复杂，形成了大量造成群死群伤的灾害事件^{[2 − 4]}。

2019年7月23日20时，贵州省水城县鸡场镇坪地村岔沟组发生高位滑坡，大约1.50×10⁶ m³的滑体高速下滑，铲刮坡面的松散堆积物，形成体积约为2.00×10⁶ m³的堆积体，造成了43人遇难，9人失踪，直接经济损失约10 300万元。水城滑坡发生后，很多学者对其启动机理、动力学过程、铲刮作用等进行了分析。Zhang等^[5]基于DAN模拟和滑震信号反演对水城滑坡动力学过程进行了分析；Zhuang等^[6]采用电镜扫描、高密度电法、数值模拟等手段多源协同研究了水城滑坡的高速远程动力学特征；高浩源等^[7]分析了水城滑坡运动过程和冲击铲刮特征，认为铲刮可分为冲击嵌入-剪切推覆-裹挟混合三个阶段；郑光等^[8]阐述了水城滑坡发生的成因机理，并分析了滑坡的动力学过程；李华等^[9]探讨了公路切坡和持续降雨对水城滑坡稳定性的影响，并分析了滑坡的形成机制。刘建强等^[10]基于主成分分析方法，研究了水城滑坡水−岩（土）作用过程与对滑坡稳定性的影响。Guo等^[11]基于地形地貌对水城滑坡动力学过程进行了详细的研究。

水城滑坡后，斜坡残留大量滑体堆积物，在降雨诱发条件下，存在二次滑坡的隐患。然而，当前对水城滑坡的研究主要集中于“7•23”水城滑坡的失稳机理与动力学过程，针对滑后的残留堆积体和边坡的稳定性分析与潜在动力致灾范围的研究还十分欠缺。

本文基于无人机、钻探等滑坡现场调查手段，探明了水城滑坡地质环境条件，阐述了滑坡的基本特征，分析了滑坡成因与破坏模式；基于传递系数法，首次对滑坡后的残留堆积物与附近不稳定斜坡进行了不同降雨条件的稳定性分析，并对二次滑坡动力致灾机理进行了研究分析。

1.   水城滑坡地质环境条件

1.1   地质条件

水城位于云贵高原中部，地形总体西北高、东南低，属于侵蚀、溶蚀地貌。构造位于扬子准地台，断层发育，地质构造复杂，节理裂隙较为发育，风化、卸荷等地质作用明显，对水城滑坡的发育起到综合控制影响作用。

水城区气候温和，雨量充沛，年平均气温11～17 °C，降水主要集中在6—8月，年均降雨940～1450 mm。据六盘水市气象局提供资料，滑坡发生前一周内水城区日平均降雨量较大。7月18日20时—23日20时，水城区鸡场镇累计降雨量鸡场站点141.8 mm，坪地村站点189.1 mm。7月23日15时03分六盘水市自然资源局联合市气象局发布地质灾害气象风险黄色预警。

研究区滑前地形地貌图如图1所示，斜坡顶部高程约为2070 m，坡底高程1250 m；整体地形较为平缓，滑坡区域发育有两条冲沟；斜坡整体呈折线型，存在三级缓坡。

图  1  （a）滑前地形地貌图和 （b）滑前两侧冲沟（Google Earth，2018.11.14）

Figure  1.  (a) Pre-event topography and geomphology map of the Shuicheng landslide and (b) Pre-event gullies on both sides (Google Earth, 2018.11.14)

下载: 全尺寸图片 幻灯片

滑坡主要地层岩性为二叠系峨眉山玄武岩组；区内出露地层有第四系残坡积，下三叠统飞仙关组（T₁f），上二叠统龙潭组（P₃l），上二叠统峨眉山玄武岩组（P₃β）及中二叠统茅口组（P₂m）。

滑坡区域内人类工程活动强烈，包括斜坡中下部居民房屋及斜坡上部X244公路的建设；其中，X244公路1984年建成，2016、2018年两次改扩建，未实施防护工程，对周边地质环境破坏较为严重，削弱了滑体的抗剪强度，对滑坡的形成有一定的影响。

1.2   滑坡基本特征

水城滑坡发生前后影响对比见图2。滑坡主滑方向20°～26°，整个滑坡区域呈长条形，平均坡度约24°，滑坡体总长1250 m，宽200～600 m，滑坡体后缘陡壁高程约1665 m，坡脚岔沟处海拔约为1200 m，相对高差465 m，沿主滑方向最大运动距离约为1250 m，滑坡面积约4×10⁴ m²，滑坡体积约2.00×10⁶ m³。

图  2  （a） 滑坡滑后影像图 和（b） 滑坡滑前影像图

Figure  2.  (a) Post-event image of the Shuicheng landslide and (b) Pre-event image of the Shuicheng landslide

下载: 全尺寸图片 幻灯片

滑坡源区后缘为小平台，滑坡发生后在源区后部形成了一处高差50～60 m的玄武岩基岩后壁（图3），壁面呈弧形起伏，为一贯通性较好的长大节理面组成，基岩壁壁面产状：20°～40°∠38°～53°（顺坡向面理）；岩体发育三组主要片理：①40°∠53°（外倾），②20°∠38°（外倾），③184°∠52°（内倾）；滑坡后壁上覆为厚约2 m的残坡积层，下伏为弱风化块状岩体，岩性为杏仁状玄武岩、凝灰岩。基于原公路处钻孔结果，钻探施工深度超过原路面高程近10 m均未见到原公路路基物质成分，故推测滑坡剪出口位于县道以下位置，高程范围为1505～1515 m。

图  3  滑坡源区后缘基本特征

Figure  3.  Basic characteristics of the rear edge of the landslide source area

下载: 全尺寸图片 幻灯片

滑坡左侧残留体较薄，至左侧边界陡坎下部，残留体被水流冲蚀后，底部出露墨绿-绿白色凝灰岩基岩面。滑坡残留体与滑坡左侧边界之间为陡坎，陡坎为土层覆盖，颗粒粒径较细，推测为弯道超高导致的土体薄层覆盖。

滑坡右侧为小山脊，岩土体滑走后在滑坡区右侧和右后侧形成了基岩陡壁（图4），高差30～40 m。由于陡壁突然临空，上部的残积层和强风化岩体二次失稳，形成数个小型垮塌体覆盖在滑坡残留体后部，这些小滑塌宽度30～40 m，前后高差约20 m，体积6000～8000 m³。原位于滑坡区右后侧的一处混凝土民房随滑体滑移后破坏（图3、图4），经测量运动距离约122 m，运动方向为24°。

图  4  滑坡右侧边界特征

Figure  4.  Characteristics of the right boundary of the landslide

下载: 全尺寸图片 幻灯片

基于研究区域的基本孕灾条件、工程活动扰动与滑后侧面临空条件，水城滑坡源区右侧存在不稳定隐患区域，如图4标识所示。不稳定斜坡右侧边界根据地貌特征，以其右侧冲沟为界，如图4黄线标识所示。斜坡坡向为25°～32°，整个斜坡区域呈“上窄下宽”的长条形，斜坡现状地形整体呈起伏波状，斜坡单元后缘区域较平缓，中部有一个平台，前缘为缓坡区（坪寨组居民区），坡度约15°～32°，斜坡体总长1030 m，宽70～305 m，斜坡体后缘高程约1595 m，坡脚坪寨组前缘冲沟处海拔约为1155 m，相对高差440 m，整个斜坡面积约19.4×10⁴ m²，根据现场调查及物探结果显示，斜坡岩土体结构从上至下为第四系残破积含砾粉质黏土、碎石土及上二叠统峨眉山玄武岩组（P₃β）的玄武岩，第四系堆积层及散体-碎裂结构的强风化玄武岩平均厚度约15 m。第四系残破积含砾粉质黏土、碎石土，碎石块粒径一般为3～6 cm，碎石含量约35%，偶含直径大于1.0 m的块石，碎块石原岩成分为玄武岩；在坪寨组居民区前缘坡脚公路边坡处可见堆积层具一定韵律性，如图5所示。图5中，①为碎石土，顶部为根植层；②③为含碎石粉质黏土，分界线处渗水。

图  5  斜坡底部出露土层特征

Figure  5.  Characteristics of the exposed soil layers at the slope bottom

下载: 全尺寸图片 幻灯片

2.   滑坡成因分析

水城滑坡的变形特征、产生现象并结合环境地质条件综合分析，滑坡发生主要形成因素为地形特征、斜坡岩土体特征、水文地质条件及极端天气、人类工程活动等。

2.1   形成条件

滑坡发生前斜坡整体地形较为平缓，滑坡部位斜坡呈波状起伏，有三级缓坡平台（后缘顶部平台、中上部平台及中下部缓坡平台），中上部平台区域（滑源区）没有系统的排水通道，不利地表水排泄，且在斜坡上呈“陡缓陡”地形，陡坡部位临空条件较好，坡度为35°～50°，临空高度达大于10 m，为滑坡失稳提供了条件。

调查区覆盖层主要成分为斜坡中上部的残坡积粉质黏土及中下部的滑坡堆积层碎石土，滑坡堆积层较厚（图5），下伏地层为上二叠统峨眉山玄武岩组（P₃β）的玄武岩，岩体节理裂隙发育，强风化层呈散体-碎裂结构（图6）。

图  6  滑坡右侧边界处山脊岩体特征（强风化玄武岩）

Figure  6.  Characteristics of the ridge rock at the right boundary of the Shuicheng landslide (highly weathered basalt)

下载: 全尺寸图片 幻灯片

综合现场调查、钻孔揭露及物探成果，上覆土层及破碎强风化玄武岩厚度约15～41 m；含大量碎石，抗剪强度较低，在受扰动的情况下，因自身重力容易下滑。其次，大量地表水易下渗在岩土交界面及强风化层界面上汇聚，长期受地下水浸泡影响，强风化玄武岩和表层黏土层抗剪强度降低，进而造成斜坡岩土体失稳。

地下水的主要类型为第四系孔隙水和基岩（玄武岩）裂隙水，主要补给源为大气降水。根据现场调查来看，滑坡区受大气降雨影响明显。滑坡体后缘汇水面积较大，降雨条件下，后缘大量地表水漫流下渗，致使地下水水位抬高，并向滑坡体前缘斜坡体上排出。高水位的地下水对堆积体产生强大浮托力，不利斜坡稳定。

2.2   诱发因素

降雨是水城滑坡形成的主要诱发因素。首先，雨水浸润后缘岩土体自重增大；其次，降雨导致滑面强风化玄武岩和黏土层软化，抗剪强度降低，后缘岩土体剪切下错，对中部岩土体产生水平推力，推动中部及前缘岩土体滑动，从而促进了滑动面的形成；最后，大量雨水沿覆盖层下渗，富集于滑带处，抬高了地下水水位，对滑体土产生强大的浮托力，且长期浸泡降低了前缘抗滑岩土体的物理力学指标，从而引起坡体失稳。

滑源区公路建设开挖，改变了原有地形地貌，破坏原始斜坡力学平衡，卸荷作用强烈，且公路开挖形成的高陡边坡又为滑坡的形成提供了良好的位移空间。

综上来看，滑坡区处于地质破碎带，覆盖层松散且较厚，岩体风化强烈，山高坡陡，强降雨导致土体和强风化岩体饱和，形成了地下高水位，降低了岩土体的摩擦系数及抗滑力，斜坡上 X244 县道施工坡面清理降低了山体阻滑的能力，在重力和降雨的作用下，导致山体失稳，并转化为高位远程滑动，形成特大型滑坡地质灾害。滑坡体变形、破坏模式分析详见表1。

表  1  滑坡成因分析一览表

Table  1.  Overview of landslide causation analysis of the Shuicheng landslide

主要 影响 因素	地形 地貌	滑坡体发育于波状起伏地形的斜坡体上，后缘为平台区，汇水面积较大。坡体呈“陡缓陡”交替分布，地面起伏较大，滑源区没有系统的排水通道，不利地表水排泄。斜坡上呈“陡缓陡”地形，陡坡部位临空条件较好，坡度35°～50°，临空高度大于10 m，为滑坡失稳提供了条件
	地层 岩性	根据钻探及物探资料显示，滑坡体主要为第四系坡堆积层及强风化玄武岩，强风化层呈散体-碎裂结构，抗剪强度较低；其次，大量地表水易下渗在强风化层界面上汇聚，长期受地下水浸泡影响，界面处岩土体强度易降低，不利坡体稳定
	水文 条件	根据现场调查来看，滑坡区受大气降雨影响明显。滑坡体后缘汇水面积较大，降雨条件下，后缘大量地表水漫流下渗，致使地下水水位抬高，并向滑坡体前缘斜坡体上排出。高水位的地下水对堆积体产生强大浮托力，不利斜坡稳定
	诱发 外因	降雨是滑坡形成的主要诱发因素，一方面增加了坡体自重，另一方面大量雨水沿覆盖层下渗，富集于滑带处，长期浸泡降低了滑带土的物理力学指标，不利于坡体稳定。同时后缘大面积的汇水抬高地下水水位，向沟道排出，对堆积体产生强大的浮托力，降低坡体稳定； 斜坡上公路开挖形成高边坡，改变了原有地形地貌，破坏原始斜坡力学平衡，且公路开挖形成的高陡边坡又为滑坡的形成提供了良好的位移空间
变形模式		滑坡在强降雨作用下，大量雨水下渗浸泡后，产生强大浮托力，坡体在自重作用下，易沿底部强风化层界面向下部临空侧整体发生滑动变形，为一推移式变形破坏滑坡
破坏模式		综合现场调查、钻孔揭露及物探成果，水城区鸡场镇坪地村岔沟组滑坡物质组成分为三层，自上而下分别为：第四系堆积层，强风化玄武岩及中风化玄武岩层；根据钻探、物探资料及现场变形特征推测，滑坡的破坏模式可能为在堆积层及强风化玄武岩层范围内产生似圆弧型滑动，及在堆积层与基岩界面产生折线滑动

下载: 导出CSV 

| 显示表格

3.   斜坡稳定性与二次滑坡动力致灾分析

3.1   稳定性分析

水城滑坡后，斜坡残留大量滑体堆积物、两条冲沟之间存在未滑区域；沿1−1′剖面，滑体强风化层较厚，拟定强风化层界面为斜坡整体滑动的潜在滑动面，潜在剪出口位于斜坡底部（图7黄色线）；将1−1′剖面原水城滑坡剪出口作为对照组，采用滑前地形，累积滑源区条形地基的混凝土民房作为附加荷载，稳定性分析验算。1−1′剖面对照滑面与主要滑面的稳定性分析模型如图8a—b所示。

图  7  滑后稳定性分析区域

Figure  7.  Post-landslide stability analysis area

下载: 全尺寸图片 幻灯片

图  8  1−1′剖面示意图

注：a为对照滑面稳定性模型；b为主要潜在滑面稳定性模型。

Figure  8.  The longitudinal profile of cross-section 1−1′

下载: 全尺寸图片 幻灯片

其次，滑后地形导致右侧边界附近斜坡临空，卸荷严重；在强降雨条件下，原水城滑坡区域边界右侧外容易形成不稳定斜坡，如图7中2−2′所示。对于右侧不稳定斜坡2−2′剖面（图9）：①右侧斜坡地形地貌、岩土体结构、水文地质等条件基本与已滑动区一致，加之X244公路建设期间工程活动影响，切坡卸荷现象严重，结合钻孔资料，钻探施工深度超过原路面高程近10 m均未见到原公路路基物质成分，故拟定右侧不稳定斜坡主要潜在剪出口位于县道以下位置；②基于右侧斜坡现状基本特征，堆积层及强风化层较厚，拟定强风化层界面长节理面作为次级潜在滑面，剪出口位于右侧斜坡缓坡平台后缘陡坡底部。2−2′剖面主要潜在滑面与次级潜在滑面的稳定性分析模型如图9 a—b所示。

图  9  2−2′剖面示意图

注：a为主要潜在滑面稳定性模型；b为次级潜在滑面稳定性模型。

Figure  9.  The longitudinal profile of cross-section 2−2′

下载: 全尺寸图片 幻灯片

对水城滑坡区域基覆面土体进行现场取样（图8钻孔处），取土体的残余剪抗剪强度指标；由滑带土样室内试验成果数理统计得，天然抗剪强度c值为10.2 kPa，φ值为22.6°；饱和抗剪强度c值为9.8 kPa，φ值为22.3°。

对右侧不稳定斜坡区域基覆面土体进行钻孔取样（图9钻孔处），取土体峰值剪抗剪强度指标；由滑带土样室内试验成果数理统计得，天然抗剪强度c值为10.8 kPa，φ值为23.1°；饱和抗剪强度c值为10.1 kPa，φ值为22.5°。

选取①自重+天然水位作为基础工况、②自重+ 1/3饱水作为水城滑坡降雨工况以及③自重+2/3饱水作为极端降雨工况进行分析。根据现场野外调研情况，滑坡后壁上覆为厚约2 m的残破积层，呈浸湿状态。基于李华等^[9]研究，水城滑坡后缘玄武岩壁见大量的水从玄武岩-残坡积土间隙中大量渗出；此外，李华等认为存在3 m 的毛细水上升高度，将地下水位线以下的风化玄武岩视为饱和状态；因此，结合稳定性分析中滑体条块的厚度数据，选取了1/3饱水平均厚度作为水城滑坡降雨条件下稳定性分析的入渗。2/3饱水厚度是对应百年一遇降雨或郑州“7•20”暴雨等极端情况做的饱水厚度假设，郑州“7•20”暴雨短时内高达200 mm/h的雨型达到了超蓄产流条件，土壤基本处于饱和状态；此外，水城滑坡汇水面积较大，因此本文选取了2/3饱水厚度作为极端降雨条件下的土壤饱和厚度。

本文基于传递系数法对斜坡稳定性进行计算，其控制方程为：

式中：$R$——抗滑力；

$T$——下滑力；

$\Psi$——传递系数，下标代表条块序号。

依据《滑坡防治工程勘察规范》评价标准，结果如表2所示。

表  2  稳定性计算结果表

Table  2.  Slope stability calculation results

剖面	滑面	计算工况	稳定系数	评价
1−1′	主要滑面	自重+天然水位	1.174	稳定
		自重+ 1/3饱水	1.073	基本稳定
		自重+2/3饱水	1.035	欠 稳定
	对照滑面	自重+附加荷载+天然水位	1.122	基本稳定
		自重+附加荷载+1/3饱水	1.037	欠 稳定
		自重+附加荷载+2/3饱水	0.996	不 稳定
2−2′	主要滑面	自重+天然水位	1.135	基本稳定
		自重+ 1/3饱水	1.021	欠 稳定
		自重+ 2/3饱水	0.989	不 稳定
	次级滑面	自重+天然水位	1.302	稳定
		自重+1/3饱水	1.105	基本稳定
		自重+ 2/3饱水	1.064	基本稳定

下载: 导出CSV 

| 显示表格

在自重+天然水位工况时：1−1′剖面的主要滑面与对照滑面稳定系数为1.174和1.122，处于稳定与基本稳定状态；2−2′剖面的主要滑面与次级滑面的稳定系数为1.135和1.302，处于稳定与基本稳定状态；

在自重+1/3饱水工况时：1−1′剖面的主要滑面与对照滑面稳定系数为1.073和1.037，处于基本稳定与欠稳定状态；2−2′剖面的主要滑面与次级滑面的稳定系数为1.021和1.105，处于欠稳定与基本稳定状态；

在自重+2/3饱水工况时：1−1′剖面的主要滑面与对照滑面稳定系数为1.035和0.996，处于欠稳定与不稳定状态；2−2′剖面的主要滑面与次级滑面的稳定系数为0.989和1.064，处于不稳定与基本稳定状态。

对水城滑坡区域（1−1′剖面）及右侧不稳定斜坡（2−2′剖面）稳定性进行的计算分析，由于对照组计算结果与滑坡实际情况基本吻合，验证了原滑坡区域残留堆积物、未滑区域与右侧不稳定斜坡在暴雨工况下的稳定性计算结果可信。在极端降雨工况下，水城滑坡残留堆积物与未滑区域处于欠稳定状态，右侧不稳定斜坡沿主要滑面处于不稳定状态。极端暴雨下，右侧不稳定斜坡首先发生破坏。因此，需要对右侧不稳定斜坡进行动力致灾分析。

3.2   右侧不稳定斜坡动力致灾分析

高速远程滑坡-碎屑流往往呈现出流态化的动力学特征^[12]，滑体运动控制方程的形式与浅水方程类似，其笛卡尔坐标系下的守恒方程表述为：

式中：h——滑体厚度；

u、v——x、y方向的速度；

t——时间；

z_t (x, y)——底坡高程；

τ——摩擦项，其下标代表摩擦方向。

Xing等^[13]基于环剪试验与数值模拟方法，验证了包括摩擦项和湍流项的 Voellmy流变模型适合模拟高速远程滑坡动力学过程。Voellmy流变模型为：

式中：f——摩擦系数；

ξ——湍流系数；

σ_z——正应力，定义方程为：

式中：R——曲率半径；

α——坡度。

根据稳定性分析结果，基于Matlab编程，采用Lax-Wendroff差分格式对方程组（2）—（4）进行离散并数值求解，得到右侧不稳定斜坡破坏后的动力致灾过程（图10）。

图  10  不稳定斜坡动力致灾过程

Figure  10.  Dynamic disaster-causing process of post-failure unstable slopes

下载: 全尺寸图片 幻灯片

如图10所示，当不稳定斜坡（2−2′）沿主要滑面发生破坏时，滑体沿左前方剪出，向左侧水城滑坡源区发生倾倒；5 s后土质滑体逐渐解体，其中大量滑体沿原水城滑坡的滑移路径在重力作用下加速滑动，少量滑体解体后与水城滑坡右侧边界外下滑；20 s后，滑体主体沿水城滑坡右侧冲沟滑移，少量滑体沿水城滑坡右侧边界外沟谷滑动；最终滑体主体停在水城滑坡右侧冲沟前缘，少量滑体沿右侧边界外的沟谷堆积分布。

水城滑坡滑体高速滑移，铲刮裹挟了大量的表层松散土体，为右侧不稳定斜坡二次滑坡创造了优势滑移路径；右侧冲沟内的曲率半径较小，导致滑体在冲沟内受到的支撑力与摩擦阻力较大，造成二次滑坡在冲沟内的减速与堆积。

受控制方程的限制，模拟中将滑床视为不可铲刮的刚体；然而，水城滑坡后，仍有大量滑坡堆积物残留于斜坡上。不稳定斜坡二次滑坡冲击、铲刮等运动过程中，铲刮裹挟松散的残留堆积物，甚至造成残留堆积物的破坏与启动，可能导致滑体大量增积、滑移路径增加和动力致灾范围变大，严重威胁斜坡下部岔沟组及坪寨组居民及斜坡前缘区域内公共设施的生命财产安全。

4.   结论

（1） 水城滑坡处于地质破碎带，覆盖层松散且较厚，岩体风化强烈，山高坡陡，强降雨导致土体和强风化岩体饱和，形成了地下高水位，降低了岩土体的摩擦系数及抗滑力，斜坡上X244县道施工坡面清理降低了山体阻滑的能力，在重力和降雨的作用下，导致山体失稳，并转化为高位远程滑动，形成特大型滑坡地质灾害。

（2） 对水城滑坡区域及右侧不稳定斜坡稳定性进行的计算分析，基于对照组验证了稳定性计算结果可信。在百年一遇暴雨工况下，水城滑坡残留堆积物与未滑区域处于基本稳定状态，右侧不稳定斜坡沿主要滑面处于欠稳定状态。

（3） 水城滑坡后，仍有大量滑坡堆积物残留于斜坡上。极端降雨可能诱发不稳定斜坡二次滑坡运动，滑体铲刮裹挟松散的残留堆积物，可能导致残留堆积物的破坏与启动，造成滑体大量增积、滑移路径增加及动力致灾范围变大，严重威胁斜坡下部岔沟组及坪寨组居民及斜坡前缘区域内公共设施的生命财产安全。