冉林，马鹏辉，彭建兵，孔嘉旭，张作蓬，杨炬，李泽坤，马哲
（1. 长安大学地质工程与测绘学院，陕西 西安 710054；
2. 西部矿产资源与地质工程教育部重点实验室，陕西 西安 710054；
3. 中国有色金属工业西安勘察设计研究院有限公司，陕西 西安 710054）
我国是世界上最大的黄土分布区，总面积约6.3×105km2，占国土总面积的6.6 %[1]。由于黄土的类卡斯特结构，导致黄土具有强烈的水敏性，天然状态下的黄土具有极高的强度，遇水后，强度急剧降低，易发生软化变形，进而形成黄土地质灾害。据统计，90 %以上的黄土地质灾害发生与水有关[2−3]。甘肃黑方台是我国最著名的水致黄土滑坡研究场地之一，被称为黄土滑坡研究的天然试验场。从1968 年至今，几十年农业漫溉使黑方台地下水位大幅抬升，塬边黄土滑坡频繁发生[4]，引起了国内外学者广泛关注，由此开展了大量的研究。例如：在滑坡类型及分布规律方面，任晓虎[5]基于多种滑坡识别指标对黑方台黄土滑坡发育早期特征识别进行研究。许领等[6]等利用IKONOS 影像与DEM 叠加重现了滑坡场景，对黑方台滑坡进行了系统的分类；
在成因机理方面，张茂省等[7−8]、王家鼎等[9]、潘攀[10]、蔺晓燕[11]对灌溉水诱发黄土滑坡的形成机理进行了较系统研究；
在风险评估方面，彭大雷[12]在研究归纳了各类型灾害的形成条件、致灾因子、成灾模式等的基础上，提出黄土滑坡潜在隐患地质识别方法。薛强等[13]基于多期DEM 数据，利用GIS 技术，通过滑坡体积与速度乘积计算滑坡强度，对黑方台南塬典型滑坡进行了危险性定量评估。然而随着滑坡的频繁发生，当地居民安全问题异常突出，而黑方台黄土滑坡的运动特征方面的研究直接关乎着滑体最终的致灾范围和灾害效应，影响当地居民生命财产。
针对黄土滑坡的运动特征研究，数值模拟是最常用也是最可视化的手段。例如：王思源[14]运用FLAC3D有限差分模拟软件，分析了黑方台焦家滑坡的实际变形破坏过程；
胡炜等[15]在实际测试及工程地质测试的基础上，基于有限差分法的Sassa 运动学模型模拟焦家崖头滑坡运动全过程，并分析滑坡的速度变化规律；
贾俊等[16]以离散元模型对焦家崖头13 号滑坡为例，研究地下水对滑坡作用机制及灌溉型黄土滑坡破坏过程和高速远程运动学特征;赵纪飞等[17]以黑方台焦家滑坡为例，基于野外调查和室内试验，通过Sassa 运动学模型模拟该滑坡泥流运动过程。考虑到黑方台地区的黄土滑坡具有很强的连续性和流动性，而Massflow 基于深度积分的连续介质力学方法，将运动物体视作连续流体，可用来模拟碎屑流、滑坡、泥石流、山洪等地质灾害动力演化全过程，且建模效率高、计算速率快，故采用该软件反演黑方台黄土滑坡较为适用，已在多个案例[12，18−21]得到应用。
文中以黑方台2019 年10 月5 日发生在党川村附近的一起黄土滑坡为研究对象，通过野外科学考察，分析该滑坡的基本特征，并采用Massflow 软件反演该滑坡运动全过程，为黑方台地区滑坡风险评估及防灾减灾提供一定的科学参考。
黑方台位于甘肃省永靖县盐锅峡镇境内，地处黄河与湟水河交汇处上游，距兰州市约40 km（图1），属于黄河Ⅳ级阶地，台塬总面积约14 km2，塬面东西向长约7.7 km，南北向宽约2.5 km。台塬被虎狼沟切割为黑台和方台，东侧为黑台，面积约9 km2，西侧为方台，面积约2 km2。台塬出露地层自上而下为：全新统崩滑堆积物（）；
上更新统马兰黄土（）；
中更新统粉质黏土（）；
中更新统砂卵石层（）；
下白垩统河口群粉砂质泥岩（K1hk）（表1）。
表1 研究区地层岩性Table 1 Stratigraphic lithology of the study area
研究区内黄土滑坡具有群体性分布特征[6]（图2），根据黑方台地区内物质组成和运动特征，将黑方台滑坡分为四大类：（1）黄土—错落型；
（2）黄土—碎屑流型；
（3）黄土—泥流型；
（4）黄土基岩—滑移型。
图2 黑方台滑坡分布图Fig.2 Distribution of landslides in the Heifangtai platform
“10·5”滑坡发生于2019 年10 月5 日，位于甘肃黑方台党川村，据当地居民介绍，该滑坡运动时间极短，不超过1 min。该滑坡为黄土—碎屑流型滑坡，基本没有水的影响，多为崩塌型的黄土内部滑动，发生于塬边裂缝发育和内部存在软弱结构面的Q3黄土中，通常情况下前期次滑坡增加坡体的临空条件，并在坡体后缘形成一定数量的拉张裂缝，为黄土—碎屑流型滑坡的发生创造条件，坡体失稳破坏之后，在运动过程中逐渐破碎解体。“10·5”滑坡存在多期性、隐蔽性、突发性及影响范围广等特点（图3）。
图3 滑坡基本特征Fig.3 The characteristics of the landslide
2.1 地形地貌
滑坡区在地貌单元上属于黄土高原丘陵沟壑区，整体位于黄河Ⅳ级基座阶地，滑坡所处党川段为目前以来滑坡发育最活跃的区域。滑坡区斜坡坡体上陡下缓，滑坡后壁坡度范围为60°～85°，滑坡前缘为较平缓开阔农田，坡度范围20°～25°，平均坡度约50°，这种后缘高陡、前缘开阔的地形特点，为滑坡的发生提供了良好的临空条件。
2.2 滑体物质及坡体结构
“10·5”滑坡的滑体物质成分主要由风成黄土组成，位于滑坡上部，厚度约35 m，岩性为马兰黄土，浅黄色，土体松散，遇水后土体体积含水量迅速增大，导致黄土黏聚力大幅降低，发生软化，坡体稳定性降低，滑动面在黄土层非饱和带，剪出口在斜坡中上部，多在坡脚形成堆积。坡体及周边出露地层主要有马兰黄土、粉质黏土、砂卵石层、粉砂质泥岩[22−24]（图4）。
图4 “10·5”滑坡剖面图Fig.4 The profile of the “10·5” landslide
2.3 滑坡形态与规模
根据野外实地调查，该滑坡整体平面形态呈长椅型，后壁近似圆弧状，后壁宽约110 m，顺坡长约355 m，后缘高程1 710 m，前缘高程1 605 m，相对高差105 m，坡体后缘后退约11 m，滑源区滑体平均厚度30 m，滑体在坡脚形成堆积区，呈扇形展布，堆积区地形较为平缓，滑体极为松散，滑坡前缘堆积区宽度不一，范围在70～140 m（图3），滑坡区域面积为700 m2，体积约2.7×104m3，属于小型黄土滑坡。
2.4 滑坡变形特征
坡体后缘发育数条拉张裂缝，平均宽度5～20 cm，走向平行于坡边，竖向上呈上宽下窄的特点，深度受坡高、坡度等因素控制，分布形态为直线型，呈现出渐近式后退的趋势，有明显错台，其错台最大下错量可达90 cm，平均下错大于20 cm。选取一条具有明显变化的裂缝1，通过几期不同时间段的对比，裂缝横向长度与宽度逐渐增大，在201912 期前，斜坡土体沿裂缝1 贯通发生崩塌滑移，并在坡顶后部发育多条拉张裂缝（图5）。
图5 滑坡后缘错台及局部裂缝变化过程Fig.5 Dislocation of the trailing edge of the landslide and the change of local fractures
2.5 滑坡运动过程特征
据野外考察发现，滑源区土体高度临空，后缘受拉张应力影响，裂缝不断扩张，最终形成滑动面发生剪切破坏（图6），滑源区土体在重力作用下失稳破坏以235°方向开始滑动，坡下地形无约束，滑坡在运动过程中未发生较大的改道行为，滑坡东临山脊，部分滑体受到右侧山脊的碰撞，滑至底部平缓地区时，向西南方向发展。滑坡剪出口位置较高，滑体储有较大势能，滑坡前缘到达斜坡底部之后继续向前运动，穿过底部农田和灌溉渠道直至停止，滑坡前缘最远处距公路60 m（图7）。
图6 滑坡发育过程Fig.6 The developping process of landslide
图7 滑坡运动过程路径Fig.7 Motion path of the landslide
3.1 运动模型建立
采用Massflow 软件进行滑坡特征反演，该软件是基于深度积分的连续介质力学方程而建立的有限差分法数值软件，能够很好地考虑运动过程中地形地貌对其运动的影响。该软件把复杂地质模型简化为上下界面，计算公式如下[18,21]：
式中：ρ1——流动层密度/（kg·m−3）；
τ1——上部流动层对底部的剪力/kPa；
z10——顶部流动层的基底边界/m；
E1——底部稳定层的夹带率；
u1、v1——材料边界速率/（m·s−1）；
h1——流体厚度/m；
g——重力加速度，取9.8 m/s2。
该软件内有Coulomb、Manning、Voellmy 三种计算模型，其中Coulomb 模型适用于滑坡、岩崩等碎屑流（固相）、Manning 模型适用于溃坝、洪水等水力学灾害（液相）、Voellmy 模型适用于泥浆，泥石流等灾害（固液相）。文中采取Coulomb 计算模型，进行初步模拟计算，Coulomb 模型其主要计算方程如式（2），关于软件的基本力学方程，可参考欧阳朝军等的相关研究[18]：
式中：τb——基底摩檫力/kPa；
c、φ——土体黏聚力和内摩擦角/（°）；
g——重力加速度，取9.8 m/s2；
h——滑坡竖向厚度/m；
λ——基底液化系数，取值为0～1。
关于滑坡前后地形数据及影像的获取，文中通过无人机航测获得，将所获DEM 数据导入ArcGIS10.7 中进行处理，得到滑坡物源区范围，将滑坡地形和物源区栅格数据转换为Massflow 能够识别的ASCII 格式。
3.1.1 参数选取
模拟所需参数有土体密度（ρ）、黏聚力（c）、基底液化系数（λ）、岩土体内摩擦角（φ），参考黑方台地区已有室内实验数据及参数[12,25]，将模拟所需非关键参数土体密度与内摩擦角视为定值。对于黏聚力与基底液化参数，这两者主要影响着滑坡的堆积形态和堆积厚度，滑坡的基底液化系数和黏聚力通过控制变量法确定，主要参数取值见表2。
表2 模拟计算参数Table 2 Simulation parameters
通过反演计算得到黏聚力和基底液化系数最佳取值，滑体物质主要是黄土，滑床基本为老的滑坡堆积体，因此滑坡黏聚力可简化为一个定值，黄土浸水后其黏聚力大幅下降[26−27]，故黏聚力取值范围设为0～7 kPa。第一轮反演将基底液化系数视为定值（取0.4），黏聚力在0.5 kPa 的梯度取值条件进行反演试算，模拟结果与滑坡实际范围进行对比，确定黏聚力值；
第二轮反演在黏聚力为5.5 kPa 时，基底液化系数 λ以0.1 的梯度进行反演试算，发现在0.6～0.8 区间结果较好，进而将参数以0.01 的梯度细化，得到更精确的基底液化系数λ。
3.2 反演结果分析
根据表2 中模拟计算参数进行计算，模拟滑坡不同时间段的运动特征及滑坡体的堆积范围和厚度，直至运动停止。该滑坡为多期次滑坡，所获滑坡前后地形数据时间分别为2019 年5 月及2019 年12 月，通过前后DEM数据所获滑源区土体大于“10·5”滑坡实际滑源区，并且期间滑坡上部土体存在多次崩落的可能性，崩落土体因重力作用落于坡面上，改变了滑坡下垫面微地形，使得堆积区右侧地形高于左侧，滑体向下运动过程中向两侧扩展，因滑坡右侧为一山脊，故在滑坡模拟过程中，滑源区土体在重力作用下向下运动时，东侧滑体受到右侧山脊的碰撞，发生一定程度的破碎、折返运动，随即向西南方向发展运动堆积，模拟所得结果与实际堆积范围出现些许偏差。
3.2.1 堆积厚度分析
通过对滑坡各时间段运动过程的模拟，得到了不同时间段的堆积范围（图8），滑坡体初始厚度最大约16 m，滑坡启动后，快速冲向坡脚和斜坡两侧，并最终堆积在滑坡前缘耕地处。图8 中（b）可看出t=5 s 时滑坡已启动，滑源逐渐向坡脚滑动，此时滑坡最大堆积厚度为7.08 m，t=10 s 时滑坡主体部分已基本抵达坡脚，最大堆积厚度为3.36 m，t=15 s 时，滑坡在滑坡前缘耕地处堆积，t>15 s 时滑坡体在前缘耕地处缓慢流动至停止，45 s 左右时滑坡已基本停止运动。
图8 滑坡不同时刻堆积范围Fig.8 The simulation results of landslide accumulation boundary at different times
滑坡运动是一个动态运动过程，为定量分析滑坡堆积厚度分布特征，选取t=45 s 的模拟结果，在滑坡的主滑方向、滑源区与滑坡前缘三个位置提取堆积厚度，与真实堆积情况对比。根据滑坡前后高程数据对比，获得滑坡真实堆积厚度分布（图9）。对于主滑方向的1-1′剖面，滑移距离在约355 m 处基本达到收敛，滑坡运动方向的反演堆积厚度与实际堆积特征相符，堆积区堆积厚度整体趋势为先增大减小，平均堆积厚度逐渐变小，最大堆积厚度约为0.82 m。对于3-3′剖面，因滑坡前缘地形无约束且滑体流动性较强，在滑坡物源体积一定的情况下，滑体向前运动的同时向两侧扩散堆积，主要以散状堆积在坡脚，呈扇形展布，实际平均堆积厚度约0.63 m，模拟平均堆积厚约0.49 m。由2-2′剖面看出，在滑坡滑源区附近，坡体崩塌下滑，故滑源区堆积厚度为负值，2-2′剖面真实下滑最大深度约14.86 m，反演结果与真实下滑深度一致。
图9 滑坡剖面位置及堆积厚度图（t=45 s）Fig.9 Location of typical profile and accumulate thickness (t=45 s)
通过三个剖面堆积厚度的比较，可发现在主滑方向与滑源方向上的堆积厚度变化和实际情况一致，滑坡前缘模拟堆积厚度与实际情况存在误差，这是因为滑坡堆积区的整体厚度均小于1 m，在小范围厚度上误差体现比较明显，同时也因滑坡运动路径受地形和前期堆积体的限制，运动过程中发生了碰撞折返，故在模拟滑坡运动时，滑体到达实际堆积区范围后，继续向西南方向运动堆积，模拟所得堆积区范围略大于实际堆积区范围，使得滑坡前缘处模拟堆积厚度小于实际堆积厚度。
3.2.2 运动速度分析
图10 为该滑坡不同时刻运动速度模拟结果。为获取滑坡不同部位在不同时刻的运动速度，在滑坡主滑方向上布置速度监测点1−5 号，得出1−5 号监测点滑坡运动速度变化图。
图10 不同时刻滑坡运动速度Fig.10 The movement speed of the landslide at different times
滑坡自启动到停止运动总用时45 s，可将其运动堆积过程分为启动加速-稳定加速-减速堆积三个阶段（图11）：0～7.2 s 内为启动加速阶段；
7.2～11.3 s 为稳定加速阶段；
11.3～45 s 为减速堆积阶段，滑坡动能逐渐耗尽，最终停止。其中启动加速阶段与稳定加速阶段时间较短，仅占运动总时间的25.1 %，说明“10·5”滑坡具有突发性，坡体上部迅速崩塌滑动，在极短的时间内速度达到最大值，减速堆积阶段占据滑坡运动主要部分，约74.9 %，说明滑坡具有较强的流动性。
图11 运动速度-时间变化曲线Fig.11 Relationship between the movement speed and time
从5 个速度监测点获得的数据来看，到达每个监测点的最大速度数值相差较小，滑坡各点运动速度均表现为先快速增大后迅速减小。坡下横向地形无约束，能量消散快，故滑坡运动速度变化迅速。对于1 号监测点，坡体上部土体开始崩塌失稳，滑动开始约1 s 运动速度达到最大值16 m/s，因坡体上部存在些许土体掉落，在减速阶段速度未能减至0；
2 号监测点上，滑坡滑源物质经4 s 左右达到最大速度23 m/s，17 s 时速度趋于0；
3号监测点上，滑坡历经7 s 达到最大速度26 m/s，25 s 时速度趋于0；
4 号监测点上，滑坡历经11 s 达到最大速度24.6 m/s，30 s 时速度趋于0；
5 号监测点位于堆积区前缘，滑坡动能已耗散部分，故此处最大速度较低为18.5 m/s，在45 s 时趋于0。结合运动速度-时间变化曲线图可以看出，滑坡从后缘至前缘处逐步达到静止状态，且在滑坡不同部位，滑坡运动速度从加速直至降低为0 的时间间隔都比较小，进一步说明了“10·5”滑坡运动速度快。
文中通过野外考察研究，对黑方台党川村“10·5”滑坡的基本特征及运动特征进行了一定的讨论，并利用Massflow 数值模拟软件反演滑坡运动全过程。得到以下结论：
（1）“10·5”滑坡体积约为2.7×104m3，主滑方向为235°，最大滑距355 m，具有突发性，影响范围广的特点。
（2）根据滑坡运动速度可将“10·5”滑坡运动堆积过程分为启动加速-稳定加速-减速堆积三个阶段，自启动到停止整个过程仅持续45 s，其中滑坡最大运动速度可达26 m/s，堆积区最大堆积厚度为0.82 m。
（3）通过1—5 号速度监测点的数据可看出，“10·5”滑坡运动速度呈先快速增大后迅速减小的趋势，所用时间极短。
（4）经过Massflow 数值模拟获得的滑坡危险范围较大于真实范围，存在些许偏差，结果偏于保守。
猜你喜欢滑体坡体黑方降雨与水位波动作用下的水库边坡稳定性影响研究陕西水利(2022年11期)2022-11-29四川美姑拉马阿觉滑坡复活特征与影响因素分析中国地质灾害与防治学报(2022年4期)2022-08-30采动-裂隙水耦合下含深大裂隙岩溶山体失稳破坏机理煤炭学报(2021年11期)2021-12-09基于FLAC 3D的推移式滑坡变形破坏模式及稳定性分析四川建筑(2020年4期)2020-09-18强降雨作用下滑坡稳定性分析中国房地产业·中旬(2020年12期)2020-06-09露天矿反铲挖掘机处理滑体的方式中国科技纵横(2017年5期)2017-05-12棋规问答棋艺(2017年1期)2017-03-31干砌石护坡灌注混凝土施工方法在水电站大坝施工中的应用科学与财富(2017年2期)2017-03-15滑坡运动中滑体变形解体研究综述兵器装备工程学报(2017年2期)2017-01-01棋盘上的三十六计小溪流(成长校园)(2016年6期)2016-06-22