涂正楠，冯 君，邓应进，潘激扬
（西南交通大学土木工程学院，四川 成都 610031）
2019 年6 月至8 月，四川省甘洛县在短时间内接连遭受数次严重暴雨灾害[1-2]。持续强降雨引发了大面积山洪、滑坡、泥石流灾害，比依市村滑坡从2019 年8 月3 日降雨后发生初始变形，房屋出现裂缝，房屋基础局部垮塌，部分房屋院落地面拉裂，村道拉裂破坏。受2020 年8 月30 日强降雨影响，该滑坡于9 月3 日再次发生剧烈滑动，滑坡范围扩大、向坡下扩展，形成一个大的滑坡体。目前该滑坡还处于持续蠕滑变形中，后期在降雨作用下，存在继续扩大变形的可能[3]，共威胁31 户134 人，且滑坡还会破坏坡体上4.67×104m2耕地，共计威胁资产500 万元。
随着近年来数值模拟方法和计算机性能的飞跃，基于有限差分法[4-6]（finite difference method，FDM）、有限单元法[7]（finite element method，FEM）等方法在滑坡研究上取得极大进步，对滑坡全过程运动破坏机理的研究也逐渐从二维平面转为真实三维地形[8-11]。传统的数值模拟方法虽然精度较高，但滑坡是一个存在滑动、平移、转动的复杂过程[12-13]，其过大的变形会导致有限元网格的畸变，使计算结果失真。为了克服这一困难，一些新兴的无网格法：点物质法[14-15]（material point method，MPM）、光滑粒子流体动力学法[16-18]（SPH）、无网格伽辽金法[19]（element free Galerkin method，EFG）等被广泛应用于滑坡的研究中，均得到了较好的效果，为后人提供了理论基础与经验参考。
基于此，本文以比依市村滑坡为例，通过采用现场地质调绘、现场监测并结合数值模拟等方法，详细分析了滑坡的诱发因素，并结合大型商业软件ABAQUS 模拟研究比依市村滑坡的破坏过程，采用真实地形建模方法，引入SPH 方法，对滑坡不同关键位置进行位移、速度监测，以期为比依市村滑坡及该类滑坡后期工程防治方案等提供可靠依据。
甘洛县位于四川盆地西部，居凉山彝族自治州北部，北隔大渡河与汉源县相望，东临峨边县、金口河工农区；
南连美姑、越西两县；
西与石棉、越西接壤。滑坡区位于甘洛县中部，属构造侵蚀中单面山宽谷地貌区，由于受横断山脉的褶皱、隆起与断裂，加上河流急剧的切割作用等因素的影响，因而境内山峦起伏、沟壑纵横，河段与两侧山地高差悬殊，谷壁陡峭，河床狭窄，显现出典型的高山峡谷地貌特征。
比依村滑坡所在斜坡区总体地势东高西低，总体为较陡的斜坡地貌，斜坡总体坡向275°，区内基岩为中侏罗统新村组（J2x）钙质泥岩、粉砂，上覆昔格达组黏土岩及滑坡堆积等。滑坡处于甘洛河右岸陡坡与分水岭山峰间相对较缓的斜坡区，微地貌上为浅缓的槽谷（图1），滑坡区总体为较为开阔、顺直的槽谷、缓坡地貌，平均宽度150 m、长度约700 m，滑坡区面积约9.0×104m2。滑坡后缘右侧以剪切裂缝为界，为陡坡、局部呈陡坎状，坡度25°～40°，滑壁高6～8 m 并导致一段长35 m 村道全毁；
滑坡左侧以剪切裂缝和冲沟为界，为冲沟、右侧为低矮的山脊，形成滑壁高约2～3 m；
滑坡中部发育浅丘，高度10～15 m；
前缘为陡坡，坡度30°～35°。滑坡后缘高程1 291 m、前缘高程1 186 m，相对高差105 m，滑坡区平均坡度20°，滑坡平面形态呈“月牙”状。比依市村滑坡导致大量房屋见图2（a）、图2（d），耕地见图2（c），村道、既有公路（图2（b））遭到损毁，滑坡全貌及局部形态如图3 所示。
图1 滑坡区浅缓槽谷地貌（仰视）Fig.1 Shallow valley morphology trough landform in landslide area(upward view)
图2 比依市村滑坡危害情况Fig.2 Landslide hazard assessment in Biyi Village
图3 比依市村滑坡全貌与局部特征Fig.3 Overview and local features of the landslide in Biyi Village
根据滑坡成因及钻孔揭露滑体情况，滑体平均厚度8.7 m，滑坡体积约78.3×104m3，滑坡规模为中型。滑坡滑体主要由斜坡上的残坡积物组成，表层为薄层耕作土层，以粉质黏土为主，下层土体以黏性土为主，粉土次之，滑床为上新统昔格达组（N2x）黄褐色粉砂岩及紫红色黏土岩，节理裂隙较为发育，为水体向下渗流提供了良好的通道，滑坡工程地质剖面图2-2’如图4所示。
图4 比依市村滑坡2-2′断面图Fig.4 Cross-sectional diagram of the landslide at Biyi Village (Section 2-2′)
比依市村滑坡的发生是其内外因共同作用的结果，内部原因主要有以下四点：①滑坡区纵向上呈阶梯状斜坡地貌，地形临空条件发育，地形临空条件发育为滑坡发育提供了有利的地形条件；
②滑坡为缓倾顺向坡结构，且岩层倾角小于斜坡坡度，为易发生滑坡的斜坡结构类型；
③斜坡区表层分布第四系全新统松散岩类残坡积粉质黏土地层，出露基岩为上新统昔格达组（N2x）粉砂岩，属极软岩类。松散岩类结构松散，属易滑地层，不利条件下易沿基覆界面发生滑动变形。而极软岩类粉砂岩，层理近水平，裂隙较发育，泥质胶结，固结状态差，遇水易软化、崩解，抗风化能力差，遇水后其力学性质指标迅速降低，因此该地层结构组成的边坡易发生滑坡、边坡变形等不良地质现象；
④在水文地质方面，滑坡区由于汇水面积大，且村民生活生产用水无序排放，使坡体局部区域极为富水，另基岩风化带裂隙水也较丰富，滑坡区主要为地下水径流排泄区，地下水活动较为强烈，特别是当地下水沿基覆界面径流，或当基岩裂隙水水位抬升，裂隙水压力增大时，常成为引起滑坡变形的重要因素。由于该区具备上述四个方面滑坡发育的基础地质条件，一旦降雨等不利因素具备，即引发滑坡地质灾害的发生。据调查，2019 年汛期以来，甘洛县新市坝镇持续降雨，降雨量大，持续时间长。在这种长历时、高强度的降雨作用下，孔隙水压力上升，岩土体有效应力降低，岩土体c、φ值急剧减小，抗剪强度严重降低。孔隙水使滑体重度增大，且明显增大了下滑力，从而引发滑坡后缘下错变形。
2.1 SPH 方法简介
光滑粒子动力学方法（smoothed particle hydrodynamics，SPH）是一种拉格朗日无网格粒子法，其不仅避免了传统有限元网格重构和畸变等问题，还可以自然追踪运动界面，因此在模拟大变形仿真方面具有很强的优势。岩土体中常用的本构关系M-C[20-22]（Mohr-Coulomb）、D-P[23]（Drucker-Prager）、Cam-Clay[24-25]、Binghams[26]等均能适用于SPH 法。由于不需要网格划分和小变形假设的特性，其可以不受到网格畸变的影响，每个生成的粒子被单独赋予状态变量（质量、速度、位移、应力、应变等）和材料属性（重度、弹性模量、泊松比、内摩擦角、剪胀角、黏聚力等）。目前SPH 法在高速远程滑坡模拟、灾害防治预测已有部分研究成果[16-18]。
核近似（kernel approximation）与粒子近似（particle approximation）是SPH 的核心方法，通过将连续介质离散为SPH 粒子，分别计算分析粒子受力与运动状态。核近似的一般形式为：
式中：W（x-x′,h）被称为光滑核函数（smoothing kernel function），或简称为核（kernel），基点粒子的变量值可以通过其影响区域内相邻粒子的集合的贡献值和来近似；
h是定义光滑核函数影响区域的光滑长度，其决定了有多少粒子可以影响特定点的插值；
Ω为x的支持域。SPH 方法是基于移动粒子的求积公式，在上式中f(x)近似可以用离散形式定义每个粒子的物理量由相邻粒子基于核函数的插值和确定。
式中：mj与 ρj——粒子j的质量与密度；
N——粒子在x影响域内的粒子总数。
SPH 方法将每个元素点称为基点，所有计算参数和结果信息都记录在基点上。有连续性方程及动量方程：
式中：ρ——基点密度；
t——计算时间；
m——粒子质量；
v——基点速度；
x——基点的坐标；
T——人工黏度系数，其作用为降低计算过程中的非物理振荡；
W——光滑核函数。
2.2 模型建立
通过提取比依市村滑坡等高线、高程点数据建立真实三维滑坡模型如图5，模型中滑床采用有限元网格四面体C3D4 单元，滑体采用SPH 粒子PC3D 单元，共867 184 个SPH 粒子单元。该模型长（x方向）984.4 m，宽（y方向）550.6 m，高（z方向）245.9 m。滑体与滑床均采用摩尔库伦（M-C）本构模型。
图5 三维滑坡模型Fig.5 Three-dimensional landslide model
通过室内试验和参数反演分析相结合的方法确定得到滑坡岩土体物理力学参数，因昔格达组岩体饱水，且属强风化，直接做岩石试验较困难，本次试验将所取的昔格达组粉砂岩和黏土岩按土工试验处理，岩石的参数采用类比获得如表1 所示，滑带通过罚函数的方式实现，其摩擦系数f通过下式计算[27]：
表1 岩土体参数Table 1 Rock and soil mass parameters
式中：φ——滑带土内摩擦角/（°）；
c——滑带土黏聚力/Pa；
L——滑坡长度/m；
W——滑体重力/N；
α——滑体坡角/（°）。
为了探究滑坡运动特征与时空演化，此次模拟在滑坡前缘、滑坡后缘、滑坡左右两侧共设置了4 组监测点，每组监测点三个特征粒子，分别监测滑坡过程中滑体的速度、位移与时间的关系。
2.3 运动过程及分析
利用SPH 方法计算得到的比依市村滑坡数值模拟的全过程如图6 所示。
图6 比依市村滑坡数值模拟结果Fig.6 Numerical simulation results of the village landslide in Biyi City
模拟结果表明：比依市村滑坡全程约5 s，滑坡在1 s 内动能达到最大，如图7 所示，随后滑坡中部与滑坡右侧坡壁开始错动，滑坡右侧后缘开始出现拉裂缝，滑坡中部物质塌落并开始在前缘堆积，最终形成鼓胀区。第2 秒开始，滑坡整体能量开始消散，到第4 秒时动能已消散殆尽，滑坡开始逐渐趋于稳定，并在第5 秒时，最终呈现出与现场勘察报告吻合的三个滑坡区域H1、H2、H3，如图8 所示：在H1 区域，滑坡左侧沉降约为2.8 m，在H2 区域后缘拉裂缝沉降约为6.6 m，在H3 区域前缘鼓胀区处，土体隆起为0.7～2 m。
图7 比依市村滑坡动能-时间曲线Fig.7 Kinetic energy-time curve of the village landslide in Biyi City
图8 比依市村滑坡最终状态竖直位移图Fig.8 Vertical displacement diagram of Biyi Village landslide in final state
滑坡前缘、滑坡后缘、滑坡左侧、滑坡右侧三个特征粒子的位移、速度与时间的关系曲线如图9—12 所示，整体上滑坡各部除滑坡后缘外速度几乎同时在第0.8 秒时达到峰值，随后趋于静止，滑坡后缘粒子峰值速度持续时间最长，滑坡后缘运动时间与速度均大于其余滑坡各部位，表明比依市村滑坡是在滑坡后缘持续推动下造成的。最终滑坡前缘位移峰值约为1.6 m，速度峰值约为2.2 m/s，滑坡后缘位移峰值约为7 m，速度峰值约为3.5 m/s，滑坡中部（左右两侧）位移峰值为1.2～2.23 m，速度峰值为2.4～2.8 m/s。
图9 滑坡前缘特征粒子数据Fig.9 Particle data of the frontal features of the landslide
图10 滑坡后缘特征粒子数据Fig.10 Particle data of the rear features of the landslide
图11 滑坡左侧特征粒子数据Fig.11 Particle data of the left-side features of the landslide
图12 滑坡右侧特征粒子数据Fig.12 Particle data of the right-side features of the landslide
滑坡运动特征参数从后缘到前缘呈现明显地从大到小的分布规律，整个滑坡破坏模式表现为典型的后部推移式破坏，即在H1、H2 的推动加载下导致的滑坡。滑坡总历时5 s，整体水平运动距离约为6 m，SPH 数值模拟结果与实际变形情况基本一致。
（1）比依市村滑坡是由多种因素共同造成的，其中地形临空条件发育，缓倾顺向坡结构，不良地质，长历时、高强度的持续降雨是滑坡的主要诱发因素。
（2）利用SPH 方法得到的滑坡全过程模拟中，对于滑坡运动距离，滑坡运动趋势，滑体前缘、后缘的鼓胀高度与拉裂缝深度，滑体左侧、右侧的错动距离均有较好的模拟效果，证明SPH 方法模拟滑坡具有较好的准确度。
（3）滑坡中部与后缘陡坡处最早开始错动，形成贯通张拉裂缝，随后H1、H2 滑坡继续推动H1 滑坡，在前缘处形成鼓胀区，导致整个滑坡破坏模式表现为推移式破坏，与现场勘查结果基本一致。
（4）目前比依市村滑坡依然处于蠕变阶段，在降雨条件下或其它外界因素干扰作用下，产生进一步向前滑移变形破坏的可能性较大。若滑坡继续滑移，滑坡还可能向后扩展，导致滑坡后缘及后部两侧房屋损坏。因此利用SPH 模拟大变形的优势可以预测滑坡在极端条件下造成的危害，可以为此类灾害的机理研究和工程防治提供新的研究思路。
猜你喜欢 后缘前缘滑坡 滑坡推力隐式解与显式解对比分析——以河北某膨胀土滑坡为例河北地质(2021年1期)2021-07-21一种飞机尾翼前缘除冰套安装方式民用飞机设计与研究(2020年1期)2020-05-21机翼后缘连续变弯度对客机气动特性影响北京航空航天大学学报(2017年8期)2017-12-20浅谈公路滑坡治理北方交通(2016年12期)2017-01-15柔性后缘可变形机翼气动特性分析北京航空航天大学学报(2016年2期)2016-12-01茶南方文学(2016年4期)2016-06-12基于Fluent的滑坡入水过程数值模拟水利科技与经济(2016年6期)2016-04-22“监管滑坡”比“渣土山”滑坡更可怕山东青年(2016年3期)2016-02-28TNF-α和PGP9.5在椎体后缘离断症软骨终板的表达及意义川北医学院学报(2015年5期)2015-12-05深水沉积研究进展及前缘问题中国海上油气(2015年1期)2015-07-01本文来源：http://www.triumph-cn.com/fanwendaquan/gongwenfanwen/2024/1111/128265.html