四川省汉源县中海村滑坡动力学特征数值分析

通过对中海村滑坡灾害现场的调查研究,阐明了滑坡的致灾因素和失稳机制,基于连续介质力学理论数值模型对滑坡进行了动力参数反演分析。研究结果表明: 连续降雨和第四系昔格达半成岩的特殊岩性特征是中海村滑坡发生的关键因素,滑坡运动持续时间约70 s,滑坡孔隙水压力系数为0.8,最大堆积厚度为8 m,最大运动速度为16.2 m/s; 在整个滑坡运动形态和速度演化过程中,主要堆积区域为1级平台和3级平台,运动速度出现3峰特征,阶梯型地形是出现阶段性速度峰值的主要原因。     

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

高位远程滑坡是中国西南山区常见的一类灾难性地质灾害,其发生往往伴随有碰撞解体效应,导致滑体碎裂化,转化为碎屑流或泥石流,具有流化运动堆积的特征。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;滑坡高位剪出,由于势能转化为动能,滑坡快速达到速度峰值,并铲刮地表松散土层;由于强降雨,滑体高速运动使基底孔隙水来不及排出,导致基底摩擦力下降,降低能量损耗,滑体解体促进颗粒流化运动,减少了摩擦,也是滑坡远程运动的重要原因。      

贵州省岑巩县大榕滑坡特征及形成机制研究

2012年6月29日,岑巩县思旸镇大榕村突发大型滑坡灾害,约310104m3。大榕滑坡为古滑坡堆积区失稳。基于滑坡破坏特征分析和地质原型分析,定性判断滑坡失稳模式为蠕滑-拉裂-牵引式滑移。滑坡启动区在不合理人工填土及强降雨作用下,坡脚蠕滑并发展为失稳,由此导致主滑区坡脚支撑作用明显减弱,主滑区中下部因此滑移失稳,并牵引右侧主滑体中上部坡体逐步失稳。滑带主要位于下伏强风化基岩中。主滑体左侧向西滑移,右侧主体向SW向滑移。基于渗流场-应力场耦合数值分析,再现了滑坡失稳过程及发生机理。大榕滑坡形成机制深入研究对于西部山区类似滑坡分析及识别具有重要的参考价值。     

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

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

金沙江上游沃达滑坡发育特征与堵江危险性分析

金沙江上游沃达滑坡自1985年开始出现变形,现今地表宏观变形迹象明显,存在进一步失稳滑动和堵江的风险。采用遥感解译、地面调查、工程地质钻探和综合监测等方法,分析了沃达滑坡空间结构和复活变形特征,阐明了滑坡潜在复活失稳模式,并采用经验公式计算分析了滑坡堵江危险性。结果表明:沃达滑坡为一特大型滑坡,体积约28.81×106 m3,推测其在晚更新世之前发生过大规模滑动;滑坡堆积体目前整体处于蠕滑变形阶段,局部处于加速变形阶段;复活变形范围主要集中在中前部,且呈现向后渐进变形破坏特征,复活区右侧变形比左侧强烈。滑坡存在浅层和深层两级滑面,平均埋深分别约15.0,25.5 m,相应地可能出现两种潜在失稳模式:滑坡强变形区沿浅层滑带滑动失稳时,形成的堵江堰塞坝高度约87.2 m;滑坡整体沿深层滑带滑动失稳时,形成的堵江堰塞坝高度约129.2 m。沃达滑坡存在形成滑坡-堵江-溃决-洪水链式灾害的危险性,建议进一步加强滑坡监测,针对性开展排水、加固等防治工程。     

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

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

锁固段之间的力学作用模式

2015年8月12日,陕西省山阳县中村镇烟家沟村发生大规模山体滑坡,滑坡堆积体摧毁并掩埋了18间房屋,造成7人死亡,64人失踪。本文在大量野外地质调查的基础上,结合无人机航拍等技术,详细描述了山阳滑坡基本特征并对其进行工程地质分区,初步揭示滑坡的成因机制和运动过程。利用离散元软件PFC^3D,对山阳滑坡运动全过程进行模拟,再现了滑坡体失稳后,在重力作用下沿视倾向启动、碰撞和堆积全过程。本文取得主要认识如下：(1)山阳滑坡属于典型的顺层岩质斜坡视倾向失稳,滑坡体长期顺层蠕滑,加之表水溶蚀作用,最终导致锁固段剪断发生整体视倾向失稳破坏;(2)滑坡体在运动过程中由于前部稳定山体的阻挡,运动方向发生了两次碰撞偏转,呈现出层层超覆、拆离运动和侧向挤滑等特殊运动堆积特征;(3)依据山阳滑坡的堆积特征与滑体结构分布特征,将山阳滑坡运动破坏过程划分为7个阶段：滑坡体裂缝发展-贯通阶段、锁固段剪断-滑体启动阶段、整体高速滑动阶段、碰撞“急刹车”制动阶段、拆离-侧向挤滑阶段、二次碰撞抛撒流动阶段和断壁区崩滑阶段;(4)数值模拟显示山阳滑坡从滑坡体启滑到最终停积用时仅35 s,滑体最大运动距离500 m,最大运动速度36.2 m·s^-1。     

甘肃黑方台“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。该研究成果可为黑方台黄土滑坡风险防控提供一定的技术支持。     

金沙江白格滑坡运动过程特征数值模拟与危险性预测研究

2018年10月、11月于金沙江川藏交界处江达县波罗乡白格村先后发生两次体积约2400×104 m3和850×104 m3的滑坡,两次滑坡平均运动距离1400 m,堵塞金沙江形成堰塞湖。首次形成的堵江滑坡坝天然溃决,未造成人员伤亡;然而第2次滑坡堵塞第1次滑坡自然溃口,导致堰塞湖库容迅速增加到3.85×108 m3。政府部门立即开展抢险工作,通过人工修建溢洪道的方法成功泄洪,极大程度上降低洪水风险。本文利用PFC3D颗粒流软件模拟两次滑坡的发生、运动、堆积过程,并在反演结果的基础上对白格滑坡滑源区残留潜在不稳定部分未来失稳的运动路径和堆积范围进行预测,对其危险性进行科学评价。结果表明:（1）滑坡在重力作用下失稳,除了受初始势能的影响外,微地貌也是决定滑坡运动路径与距离的关键因素之一;（2）PFC3D颗粒流数值模拟方法适用于类似于白格滑坡这类碎屑流类型的滑坡,两次滑坡反演得到的堆积厚度、堆积范围均与真实结果相近;（3）利用两次事件反演所得参数,可以预测若滑源区潜在不稳定部分同时失稳,则形成约70 m高的滑坡坝,可能再次堵塞金沙江。     

新疆叶城“7·6”滑坡泥石流灾害调查与形成机理研究

2016年7月6日晨0时40分左右,新疆喀什地区叶城县柯克亚乡六村发生特大滑坡泥石流灾害,造成重大人员伤亡和财产损失。本研究采用地面调查、访问和遥感解译方法,从滑坡、泥石流形成的地形地貌、地层岩性和水文气象条件着手,分析了滑坡特征与成因,泥石流形成条件、灾害链过程与致灾机理,并对未来本区滑坡泥石流的发展趋势进行了预测。现场调查表明:（1）位于西昆仑山北坡中高山区的六村与七村毗邻区域,在6日凌晨发生的大暴雨,触发了38处群发性浅层黄土滑坡,属蠕滑-拉裂机制,且部分大滑坡表现为远程滑坡-泥流特征;（2）这些滑坡体堆积在“V”形沟谷中形成堰塞坝,其中六村上游发育15处滑坡和2个滑坡堰塞坝,在持续降水和小型沟谷泥石流的作用下,滑坡坝发生串联式溃决而形成堵溃型泥石流,冲毁六村居民区和道路;（3）同时,该毗邻区域两侧存在超过227处的黄土滑坡变形体,未来在区域气候由暖干转向暖湿条件下,强降雨的极端气候事件会增多,若遭遇强降雨,研究区仍然会爆发滑坡-堰塞坝-溃决泥石流模式的灾害,其规模可能比“7·6”事件还大。建议深化本区浅层黄土滑坡变形机理、临界雨量及泥流运动机理等方面的深入研究,同时加强滑坡泥石流的监测和预测研究。     

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

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

天水市大沟滑坡-泥石流运动过程模拟分析

滑坡-泥石流通常由降雨触发,是一种危害极大的灾害链。本文以甘肃天水市大沟村泥岩滑坡-泥石流为例,从滑坡后壁取原状泥岩样进行固结不排水（CU）三轴剪切试验,获取强度参数,利用LS_RAPID软件对该滑坡-泥石流运动过程进行模拟,再现了滑坡从发生到转化为泥石流的过程。试验表明饱和泥岩的摩擦角仅11,小于沟道坡度。模拟表明滑坡下滑至坡脚时滑体平均速度为5.2ms-1,对沟道滑床具有极大的冲击作用,滑床泥岩在滑坡强大的冲击作用下强度进一步降低,有利于滑坡体迅速转化为泥石流,并造成更大的致灾范围,模拟堆积范围和形态与实测基本一致。     

大光包滑坡工程地质研究

大光包滑坡是512汶川MS8.0级特大地震触发的规模最大的滑坡。滑坡位于汶川地震发震断裂上盘,距发震断裂直线距离3.0～4.5km,覆盖面积7.12km2,体积11.59108m3,是我国有史料记载以来规模最大的滑坡,也是目前世界上已知为数不多的几个10109m3以上的超大规模滑坡之一。本文作者自2011以来,在过去工作基础上,对该滑坡开展了1：2000滑坡工程地质测绘,结合物探、坑槽探及浅孔钻探等工作,编制了滑坡工程地质系列图件,从而进一步查明了大光包滑坡的平面和空间形态、滑体结构、滑面位置和堆积特征等,获得了较为完整的滑坡要素定量数据。研究成果还原了大光包滑坡发生时的真实场景;揭示了滑坡堆积地貌特征、植被分布、岩性分布以及堆积体结构特征等,认为大光包滑坡是受强震和特定地形条件、岩体结构条件（层间剪切断层和两组陡裂结构面）控制的楔形体失稳。在此基础上,根据滑坡的要素组成和运动堆积特征,对滑坡进行了分区,分为滑坡断壁区、主滑堆积区和次滑堆积区3个大区,进一步划分为10个小区,分析阐明了各区特征。最后,对大光包滑坡的发生过程和形成机理进行了概要分析,将其概括为5个主要阶段：即强震拉裂阶段、锁固段剪断及楔形体失稳阶段、高速滑动和急刹车运动堆积阶段、拆离滑动阶段、断壁崩滑阶段。研究还发现：强震作用下滑带的进一步碎裂化,以及可能出现的水击作用可能是滑坡骤然启动的主要原因。     

贵州都匀马达岭滑坡碎屑流形成泥石流特征

2006年5月18日,贵州都匀市马达岭发生滑坡,滑坡体位于马达岭沟上游采煤区,部分滑坡体滑入沟内并碰撞解体形成碎屑流堆积在沟道中上游,形成约800m长的淤积段,沟道内固体物质总量约16×104m3。滑坡后采煤区内老窑水迅速释放汇流至沟道启动沟道内物源形成泥石流运动至下游,从流通区下游逐渐淤积并少量冲出沟口,沟道内淤积长350m,堆积区体积约0.7×104m3。马达岭沟泥石流呈滑坡—碎屑流—水作用形成泥石流的连续性过程,具有典型的滑坡碎屑流形成泥石流特征。由于沟道内仍有较多物源,同时沟道内形成两处巨石淤塞形成的堰塞坝,后期强降雨作用下仍可能形成泥石流。     

云南省彝良县麻窝后山滑坡变形破坏机理离心机模拟试验

结合具体滑坡进行大型土工离心机模型试验,研究滑坡在天然状态下的变形破坏特征及危害性,再现了滑坡坡体裂缝形成和破坏后的实际运动过程,并分析探讨了滑坡的变形破坏机制和危害范围。试验结果表明:滑坡的变形破坏模式为牵引式软弱基座型蠕滑-拉裂破坏;滑坡危害性特征显示,滑坡并未发生整体破坏,只是前缘的局部垮塌破坏,滑坡运动最远距离0.41m,还原到实际滑动距离为82m。同时由于土体力学参数的衰减,暂时未滑动坡体具有潜在危害性。     

汉中市宁强县凤凰街滑坡运动特征及冲击强度的数值研究

为了对正在演化的潜在滑坡开展变形、失稳、运动及冲击致灾过程的超前模拟,以凤凰街滑坡为例,基于PFC~(3D)对降雨诱发滑坡失稳后的位移、速度、能量及冲击力等指标进行了定量研究。结果表明,凤凰街滑坡变形开裂从后缘开始,并逐步向前推进,表现为"后推前"特征,属典型的推移式滑坡;滑坡主滑时间约30s,最大平均运动速度3.3m/s,最大位移约40m;在峰值速度下滑坡动能仅占重力势能减小量的17.5%,碰撞和摩擦耗能始终占据主导地位;滑坡的发生将摧毁下方建筑物,对十几户居民的生命财产安全构成重大威胁。该研究成果可给决策者提供科技支撑,为堆积层滑坡的风险预测提供一定的参考,具有重要的防灾减灾意义。     

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

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

鄂西清江隔河岩水库茅坪滑坡蠕滑变形及其稳定性

茅坪滑坡是清江隔河岩库区一个大型的基岩古滑坡体。自1993年4月隔河岩水库蓄水以来,滑坡体变形日趋明显,表现为整体蠕滑的特征。根据该滑坡12年的变形监测资料,分析了滑坡的蠕滑变形特征及其变化趋势,认为茅坪滑坡变形具有累进性破坏的特点,目前处于加速蠕变阶段。稳定计算结果也表明,滑坡体处于一种潜在临界不稳定状态。若遇强降雨、库水位骤降、岩崩或地震等触发因素,茅坪滑坡有可能整体失稳下滑。     

贵州惠水鑫朋滑坡变形破坏机制

以贵州省惠水县宁旺乡鑫朋滑坡为例,依据野外调查成果,分析了滑坡结构及变形特征,结合滑坡所处的地质环境条件讨论了其稳定性影响因素,借助Geostudio分析了坡体后缘崩塌堆积体加载及降雨对滑坡的作用机制。分析结果表明:持续降雨和坡体后缘加载共同作用使第四系地层与基岩接触面产生压应力集中,坡脚处孔隙水压力持续增大,第四系地层因受到雨水浸润稳定性降低,在堆积体牵引下发生浅层蠕滑变形,因此认为鑫朋滑坡变形失稳机制为整体蠕滑-拉裂。     

三峡库区康家嘴滑坡破坏过程及演绎研究

三峡库区康家嘴滑坡是在其老滑坡泥石流的前缘堆积扇体中形成的新滑动变形体。滑体表面为第四纪风化土层,滑带为灰绿色粉质粘土或粘土,在降雨与库区蓄水的周期性浸泡作用下,饱和后蠕滑变形,处于欠稳定-极限稳定状态,亟需整治。本文针对其破坏过程及形态加以详细描述,并对滑坡的阶段性演绎过程加以追溯。在查清滑坡地质条件与发展趋势的基础上,通过经济对比分析,对康家嘴滑坡不同变形阶段提出了具有针对性的防治方案:在新滑坡体(强变形区)采用抗滑桩+排水沟+局部回填压脚、后部变形体采用抗滑桩+排水沟。经对滑体计算,该滑坡体处于稳定状态,达到了防灾减灾的目的。