甘肃黑方台“10·5”黄土滑坡启动及运动特征分析

数十年的塬上大面积灌溉导致甘肃黑方台地下水位不断提升,迄今为止已有约200余起黄土滑坡发生,造成了大量伤亡,给当地居民生活和财产安全带来了巨大的影响。针对黑方台黄土滑坡启动迅速、运动速度快、影响范围广的特征,文中以2019年10月5日发生在黑方台党川村附近的一起黄土滑坡为研究对象,通过野外调研和无人机手段,获取该滑坡的发育特征和地形数据。基于Massflow软件再现了该滑坡运动全过程,最终得到堆积厚度及运动速度2个方面特征,得出以下结论:（1）“10·5”黄土滑坡体积约2.7×104 m3,最大滑距355 m,属于小型黄土滑坡;（2）模拟结果显示,滑坡运动持续时间约45 s,最大运动速度为26 m/s,反演结果与实际情况较吻合。运动过程可分为启动加速阶段、稳定加速阶段、减速堆积阶段三个阶段;（3）对滑坡堆积厚度进行分析发现,该滑坡在主滑方向上平均厚度约0.59 m,最大堆积厚度为1.35 m,沿滑坡运动方向,堆积厚度呈现先增大后减小的趋势,堆积区最大堆积厚度为0.82 m。该研究成果可为黑方台黄土滑坡风险防控提供一定的技术支持。     

深圳“12.20”渣土场远程流化滑坡动力过程分析

文章采用DAN3D数值方法对深圳人工堆填体滑坡运动过程进行了模拟研究,探讨了深圳“12.20”滑坡远程动力成灾过程。通过研究得到以下几点结论:（1）滑坡后破坏运动主要分为两个阶段:前一阶段为滑源区内运动,体现了高孔隙水压力下滑剪切;后一阶段为在流通区和堆积区内运动,体现了高饱和度滑体流动（涌动）剪切。（2）饱水渣土滑坡远程流化运动分析中,摩擦模型适合模拟孔隙水压力作用下的滑源区渣土体的失稳下滑运动过程;宾汉姆模型适合模拟非牛顿流体饱和渣土体的流化剪切过程;摩擦-宾汉姆组合模型更适用于该类型滑坡全过程的反演运动分析。（3）深圳滑坡后破坏运动速度变化主要经历了“启动-加速-持速-减速”的运动过程,高含水渣土的固-流转化致使滑坡远程运动,并造成巨大伤亡损失。（4）模拟结果显示:堆积区平均堆积厚度为11 m,堆积范围为0.4 km 2,最大运动速度为30 m/s,最大速度发生于距滑坡后缘620 m处,堆积范围、堆积厚度和运动速度同滑坡实际值基本一致。上述研究思路和方法对城市地质中渣土滑坡灾害的危险区划和渣土场科学选址评估具有一定借鉴意义。     

四川泸定昔格达组滑坡灾害运动过程模拟分析

四川泸定昔格达组以半成岩为主,工程地质特性复杂,在高陡斜坡中常发生浅层蠕滑变形,在强降雨作用下失稳后可转化为泥石流。本文以四川省泸定县海子坪环环村滑坡为例,基于遥感解译、地面调查、数值模拟等方法,对滑坡发育特征、潜在失稳模式和滑坡-泥石流运动过程进行分析。结果表明,环环村滑坡主要发育于昔格达组内,以深度3～5 m的浅层变形为主,整体处于蠕滑变形阶段。滑坡平面上分为强变形区(A区)和弱变形区(B区),体积分别约为5.5×10⁴m³、5.8×10⁴m³,在不同降雨条件下,存在仅有A区下滑和A区牵引B区一起下滑并转化为沟道泥石流2种致灾模式。当仅有A区失稳下滑时,最远运动距离可达1325 m,最大堆积厚度为5.2 m,最大运动速度41.6 m/s,滑坡破坏沟口居民区及道路。当A区和B区同时失稳下滑时,最远运动距离可达1345 m,最大堆积厚度为7.7 m,最大运动速度为44.3 m/s,滑体最远能够冲至河流对岸,形成高约3 m的滑坡坝。研究结果对于深化浅层滑坡-泥石流远程致灾效应的认识和防灾减灾具有一定的...     

贵州省六盘水水城高位远程滑坡流态化运动过程分析

高位远程滑坡是中国西南山区常见的一类灾难性地质灾害,其发生往往伴随有碰撞解体效应,导致滑体碎裂化,转化为碎屑流或泥石流,具有流化运动堆积的特征。2019年7月23日发生于中国贵州省六盘水市水城县的鸡场镇滑坡是典型的高位远程流态化滑坡,滑坡前后缘高差430 m,水平运动距离1340 m,堆积体体积200×104 m3,导致21幢房屋被掩埋,51人遇难。基于野外详细调查和滑前滑后地形对比,采用DAN-W软件对水城滑坡的整个运动堆积过程进行了模拟,结果显示:水城滑坡在滑源区残留堆积体厚度最大为27 m,堆积区最大堆积厚度为15 m,滑坡碎屑流前缘最大运动速度为27 m/s,最大动能为6.57×106 J;滑坡高位剪出,由于势能转化为动能,滑坡快速达到速度峰值,并铲刮地表松散土层;由于强降雨,滑体高速运动使基底孔隙水来不及排出,导致基底摩擦力下降,降低能量损耗,滑体解体促进颗粒流化运动,减少了摩擦,也是滑坡远程运动的重要原因。      

新疆伊犁河谷黄土滑坡冻融失稳机理研究

伊犁河谷地处中国天山山脉西段,属于典型的季节性冻土区,季节性冻融作用对斜坡稳定性影响强烈。本文以皮里青河"3.24"滑坡为例,通过现场详细调查、多期遥感影像动态比对、数值模拟等方法,研究了黄土滑坡冻融失稳机理与滑坡运动特征,得到以下几点认识:(1)滑坡失稳与水密切相关,一方面是地表水侵蚀作用,另一方面是地下水的冻融作用;(2)根据变形特征,将滑坡失稳过程分为坡脚侵蚀、冻结滞水和冻融循环破坏三个阶段,模拟结果表明,地下水水位上升导致孔隙水压力增大是滑坡冻融失稳的主要原因;(3)DAN-W数值模拟软件的Voellmy模型和Frictional模型可较好地模拟滑坡的运动堆积过程,模拟结果显示,滑坡历时12.8 s,最大运动速度17.7 m/s,平均堆积厚度4.9 m,运动距离139 m。本研究为伊犁河谷地区黄土冻融滑坡的早期识别与风险评估提供了重要的技术支撑。     

917灞桥灾难性黄土滑坡形成因素与运动模拟

持续性强降雨是诱发黄土滑坡的主要因素。 2011年9月份的持续强降雨达到50a来最大值, 在陕西省关中地区诱发了多处滑坡, 造成交通中断, 带来了严重灾难。尤其9月17日发生于白鹿塬的灞桥滑坡(917灞桥滑坡), 共造成32人伤亡的灾难性滑坡事件。为了揭示滑坡成因和运动过程, 通过对917灞桥滑坡灾害的调查和分析, 揭示了灾害的特点和诱发因素, 通过滑坡运动模拟还原了其运动过程和成灾范围。(1)9.17灞桥滑坡体相对高度大约90m, 宽170m, 滑动厚度约10m; 滑动方向为60, 滑动距离约150m, 总方量约15104m3; 滑坡共发生3次滑动, 滑动方量分别为9.5104m3、3.5104m3和2104m3, 平均堆积厚度12m左右; (2)917灞桥滑坡诱发因素主要是长历时的强降雨、开挖后的高陡边坡和特殊的黄土结构性。充足的前期降雨和当日的强降雨是滑坡发生的主要激发因素; 开挖后形成的高陡边坡发育一系列裂隙, 促使滑坡的发生; 裂隙的产生为降雨的优势渗流提供了通道, 加速了滑坡的发生。(3)利用LS_RAPID对917灞桥滑坡进行模拟, 根据模拟结果, 滑坡可以分为3个阶段, 起动加速滑动阶段(0～7.5s), 沿x方向速度从0迅速增加到4.9ms-1, 沿y方向速度从0迅速增加到8.4ms-1; 滑坡体在这一阶段共滑动了65m; 减速滑动极端, 速度开始降低(7.5～14s), x、y方向的平均速度分别降到1.5ms-1和2.5ms-1, 滑坡在这一阶段向前运动了50m; 堆积停止阶段(15～28.6s), 滑坡运动速度持续降低, 堆积厚度逐渐变薄, 滑坡体共向前滑动了175m, 滑坡最大堆积厚度为14m, 与野外调查结果基本一致。结果对认识和研究此类滑坡成因机理具有重要的借鉴意义, 也可为今后该地区防灾和减灾提供参考。     

基于DAN-W模型的高速远程滑坡灾变过程分析

高速远程滑坡因其突发性、剧烈性、高速性、远程性等特点对人员和财产造成重大危害,是中国最严重的地质灾害类型之一,研究其灾变过程对中国防灾减灾具有重要的意义。DAN-W是基于拉格朗日有限差分法建立的动力模拟方法,反演模拟碎屑流运动过程,更适合高速远程滑坡的模拟分析。笔者以哀牢山地区芭蕉树滑坡为研究对象,研究芭蕉树滑坡的基本特征,利用DAN-W软件建立芭蕉树滑坡的计算模型,分析滑坡的灾变过程,依据滑坡特征考虑滑坡分区和铲刮效应,将滑动路径分为3段:滑源段、铲刮段、堆积段进行模拟计算;结果显示:滑坡前缘最大速度为26 m/s,后缘最大速度为30 m/s,后缘最大速度发生在剪出口位置,最大滑距为420 m,滑动14 s后滑体基本滑离剪出口,至运动30 s时停止滑动。所得结果与实际情况吻合度高。     

青藏高原察达高速远程滑坡运动过程与形成机理

为了了解青藏高原察达高速远程滑坡的运动过程与形成机理,运用遥感测绘、无人机地形测绘和现场勘查资料对滑坡进行分区,对滑坡形成机理进行研究,并利用PFC2D数值模拟对地震工况下滑坡运动过程进行模拟.将察达高速远程滑坡分为源区,流通区和堆积区;数值模拟结果得到滑坡平均运动速度为15~20 m/s,运动时间150 s,最大运动距离为2 800 m.察达滑坡为地震条件下诱发的高速远程滑坡,源区砾岩对上部堆积体后缘铲刮推移,使得上部堆积体产生整体变形,其运动过程可分为崩滑→铲刮→滑移→堆积4个阶段.      

复合型滑坡固液耦合过程数值模拟分析——以无山坪滑坡为例

固液耦合作用是碎屑流向泥石流转化形成复合型滑坡灾害的关键因素, 会导致成灾范围和规模放大, 是防灾减灾领域研究中的热点和难点问题之一。文中采用自主研发的滑坡后破坏数值模拟平台(LPF3D, Landslides post failure 3D), 以2014年9月强降雨诱发的重庆奉节无山坪滑坡为例, 探讨了滑坡在水动力作用下远程成灾的动力过程, 揭示了固液耦合影响机制。研究结果显示: 水动力作用在滑坡运动过程中主要体现为液化和拖曳两种, 两种力学作用的增程效应明显, 往往使得碎屑流转化为泥石流, 导致远程成灾; 基于光滑粒子流体动力学(SPH)方法的两相耦合计算模型, 考虑流体状态方程、固体黏塑性本构方程和相间作用力的共同影响, 基本还原了强降雨条件下重庆奉节无山坪滑坡两相运动过程; 数值计算结果显示无山坪滑坡最大运动速度为34 m/s, 最大堆积厚度为21.5 m, 堆积面积为0.12 km2, 最远运动距离为1300 m, 模拟结果同实际滑坡的堆积形态基本一致。综上认为, 在高位远程滑坡风险调查与预测过程中, 需充分考虑强降雨工况下孔隙水压力和固液相间作用, 基于LPF3D方法的数值模拟为高位远程滑坡的风险定量评估提供了依据。      

高速远程滑坡超前冲击气浪三维动力学分析

为研究沟谷多级折射高速远程滑坡超前冲击气浪的特征,以牛圈沟滑坡三维地形形态为原型,采用CFD软件,引入Voellmy准则定义运动阻力,在反演碎屑流运动全过程的同时,对其冲击气浪的流场特征进行了数值模拟,分析其形成的动力学机制。结果表明：(1) 牛圈沟滑坡从启动到停止历时为119 s,最大速度出现在滑体内部靠近前缘处,峰值分别为14 s的52 m/s、27 s的55 m/s和50 s的49 m/s。(2) 气浪速度的最大值50 s时为38 m/s,其压强的最大值达657 Pa,相当于11级暴风。(3) 当碎屑流前方出现高大的障壁时,来不及扩散的高压气流产生极为明显的压强梯度;当其前方遇跌坎或翻越山脊冲向山坡下方在其速度瞬间增大的同时,强烈压缩前方空气并导致其压强值局部出现陡增。     

马边滑坡运动特征及冲击强度的数值研究

滑坡运动堆积特征及其冲击强度研究对滑坡风险定量评估具有重要意义。通过对四川乐山市马边滑坡基本特征调查,利用支持向量机模型(SVM)和颗粒流方法(PFC),对滑坡岩土体细观强度参数进行反演和标定,结合UAV数据生成滑坡区高精度DEM,在此基础上,重构马边滑坡三维颗粒流数值模型,模拟并研究滑坡的运动堆积和冲击过程。结果表明:马边滑坡运动时长32 s,主滑时间16 s,运动开始5 s后速度达到峰值,为10.2 m/s;滑坡中后部岩土体运动迹线为直线型,中前部运动迹线成扩散状态,最终呈扇形堆积;滑坡在坡脚处的冲击力可达1.5×109 N,并随着坡脚距的增大,冲击力呈现出指数衰减特征。研究结果与滑坡运动过程实际视频解译结果及堆积现状基本一致,相关研究方法为滑坡定量风险评估提供借鉴。      

陕西泾阳东风高速远程黄土滑坡运动过程的模拟

陕西泾阳南塬的东风滑坡是一个典型的高速远程黄土滑坡,由于斜坡顶部的黄土台塬长期灌溉所诱发。通过野外调查,采取滑坡后缘斜坡的黄土试样,对斜坡中部的Q2黄土做天然含水率的固结不排水剪试验,对斜坡底部地下水位面附近的Q1黄土做饱和固结不排水剪试验,结果表明非饱和的Q2黄土和饱和的Q1黄土在破坏后都具有显著的应力降,但前者是由结构强度引起,而后者则是由剪切形成的超孔隙水压力引起,由此表明该斜坡一旦破坏,就会有较高的启动速度。对滑坡坡脚以外一级阶地上的砂卵层做环剪慢剪试验,测得其有效残余摩擦角,根据滑坡各部分含水状态给定孔隙水压力参数,利用Sassa. K滑坡运动模型对该滑坡运动过程进行模拟,从模拟结果反映出滑坡运动过程中的覆盖范围和速度的变化,模拟的滑坡运动距离和覆盖范围与实测比较吻合,最高速可达到37m s-1,为高速滑坡。黄土滑坡以高的初速度下滑到坡脚下一级阶地上的饱和砂卵层上,砂卵层上液化是导致其远程运动的主要原因。     

新疆伊宁县喀拉亚尕奇滑坡动力学特征研究

通过对高速远程黄土滑坡动力学特征的研究,提出黄土高速远程滑坡空间预测的模拟方法。以新疆伊宁县喀拉亚尕奇黄土滑坡为例,基于野外地质调查和无人机航拍影像图,结合滑坡研究区的工程地质条件,分析了该滑坡的基本特征和形成条件。研究发现,该滑坡的主要诱发因素是冰雪融水入渗,其孕灾模式主要为四个阶段:后缘拉裂阶段,黄土节理冻胀扩展阶段,融雪入渗失稳阶段,高速下滑阶段。同时利用Rapid模型对滑坡运动全过程进行模拟,计算得到滑坡运动持续时间为26 s,最大运动速度达到22 m/s,堆积体的平均厚度达到5 m等运动特征要素,结果表明Rapid模型可以较好的模拟分析黄土高速远程滑坡动力学效应,为黄土地区类似滑坡的成灾机理和动力学效应分析提供参考。      

基于MatDEM的烟家沟滑坡演化过程数值模拟分析

在对陕西省山阳县烟家沟滑坡基本地质特征调查的基础上,利用无人机获取滑坡高程数据,基于高性能矩阵离散元软件MatDEM建立起烟家沟滑坡运动及堆积特征的三维离散元模型,还原了滑坡启动发生的全过程,最终得到了运动速度、滑移距离以及堆积体厚度3个方面特征,取得以下认识：烟家沟滑坡属于高速滑坡;滑坡运动堆积过程分为加速碰撞堆积阶段、整体滑动堆积阶段和减速堆积阶段;其堆积层厚度分布沿运动方向呈现出先增加后减小的总体变化趋势,且堆积厚度在堆积体后部达到最大;MatDEM矩阵离散元法用来模拟滑坡演化过程具有可行性.     

舟曲江顶崖滑坡的早期判识及风险评估研究

如何提前判识滑坡变形并对其进行早期风险评估已成为地质灾害防治领域的研究热点。文章以舟曲白龙江流域江顶崖堆积层滑坡为反分析案例,进行了滑坡变形早期判识及风险评估综合研究,提出了小基线集雷达干涉（SBAS-InSAR）技术解译分析、地质-力学联合分析、动力过程数值模拟分析三者相结合的滑坡变形早期判识与风险评估全流程分析模式。基于SBAS-InSAR技术解译能够准确地判识江顶崖滑坡的分布范围及早期形变特征,江顶崖滑坡的变形破坏模式为牵引式,滑坡体长度约680 m,宽度约210 m。基于早期识别信息,地质-力学联合分析表明:江顶崖滑坡为典型的老堆积层滑坡,前缘局部变形,破坏模式为牵引式,滑坡体平均厚度约35 m,滑床整体坡度较缓,失稳后运移速度不大。选取符合江顶崖滑坡体滑移摩擦特征的库伦摩擦模型,基于深度积分连续介质方程,分析计算滑坡体的动力学过程,结果表明:滑坡体滑移速度不大,最大值约为2.2 m/s,运动方式表现为推挤白龙江河道,堵江可能性较小,并且江顶崖滑坡体前缘错动完成后,该滑坡体滑移速度从前缘到后缘快速降为0,表现为牵引式运动特征。本次分析结果与实际相符,吻合度较高,采取的综合分析方法及研究模式可用于舟曲白龙江沿岸类似滑坡的早期判识及风险评估。