基于SBAS-InSAR技术的舟曲县潜在滑坡灾害早期识别及降水相关性分析

近年来,甘肃省舟曲县白龙江两岸古滑坡的发育复活呈现活跃趋势,下游已经开发建设的大量民用设施受其威胁,因此这些滑坡的早期识别成为该地区地质灾害预防需要关注的重要课题。基于Sentinel-1A降轨SLC数据,采用SBAS-InSAR技术,对舟曲县,尤其是白龙江沿岸滑坡灾害隐患进行早期识别,获取舟曲县潜在滑坡灾害点列表,并划分其风险等级,重点分析大蠕滑量的四个典型滑坡的时序形变过程,总结舟曲县滑坡发育与降水之间的相互关系,研究结果表明：舟曲县白龙江干流及支流河谷区西北-东南方的老滑坡体均发生了一定程度的滑动,共识别出高风险等级坡体8个,分别为寺儿沟、立节北山、杭嘎村、锁儿头、三眼峪、南峪、大川乡—两河口及曲告纳乡滑坡,严重威胁15个村镇的居民生命财产安全。其中锁儿头、三眼峪和南峪滑坡的形变最为严重,对其重点分析的结果显示,监测时段内平均形变速率达-80.1421.48mm/y,形变量达-482.01493.31mm,坡体在未来复发的可能性较大,而且,近年来逐渐增加的降水更有可能促使其提前复发。舟曲县典型滑坡的后缘形变与降水量具有一定的相关性,滑坡的发生受降水影响极大,同时,地质构造、地形地貌、地震、地层岩性等也对滑坡的发生具有一定影响。     

斜坡失稳及其启程速度的折迭突变模型

临滑斜坡内聚集了巨大剪切形变能,坡体临滑释放的形变能大于滑带岩体破裂、贯通所需能量的部分,将转变成滑体弹冲剧发的动能.将坡体视为变形体,计算了重力作用下简化坡体模型的剪切形变能,计入滑带厚度,建立斜坡失稳的折迭突变模型,对滑体启动方式进行分析.研究表明,当坡体剪切刚度k小于滑带锁固段岩体抗剪力曲线软化段拐点处刚度|Q'(u1)|,即系统刚度比K=k-Q1(ut)<1时,滑动为剧动式滑坡.给出了滑体弹冲的动能表达式,由此得到的滑体弹冲剧发速度、加速度关系表明,滑体的启程弹冲效应与滑体规模成反比,给出剧动式滑体启程速度的算例.用格林公式证明了临滑坡体剪切形变能释放量,可以用折迭突变模型中的几何图形来表示.当参数变化使系统剐度比K=k-Q1(ut)≥1时,所建立的突变模型可以对缓动式滑坡进行描述.     

基于多时相InSAR技术的黄河上游龙羊峡库区不稳定边坡识别与形变监测

黄河上游龙羊峡水库的建设改变了两岸的水文地质条件,库区水位周期性升降引起的渗透和侵蚀导致库岸坍塌和滑坡局部变形,非常有必要开展不稳定边坡的识别与监测,为库区坡体失稳因素研究以及预警预报提供技术支撑。采用2015年1月至2021年8月升降轨Sentinel-1A SAR影像,基于相干点目标分析(IPTA)方法对龙羊峡库区边坡进行识别,共识别出15处不稳定边坡。其中,查纳滑坡历史上多次发生塌滑现象,现今仍处于活跃状态,对库区人民生命财产与大坝工程安全造成潜在威胁,为此对查纳滑坡开展重点监测。首先,采用DS-InSAR技术获取该滑坡升降轨雷达视线向(LOS)形变; 然后,结合升降轨形变结果进行二维时序形变分解; 最后,采用小波分析定量分析了地表时序形变反演结果与降雨量、库区水位变化的相关性。结果表明:DS-InSAR技术获取的形变场比相干点目标分析方法更加完整; 查纳滑坡东部坡体位移大于中部与西部坡体位移,且垂直向位移大于东西向位移; 查纳滑坡位移时间序列与降雨量序列相关性不高,降雨可能不会显著促进滑坡的位移; 滑坡位移与库区水位均具有很强的年度周期变化,并且滑坡位移时间序列滞后库区水位序列约半年。     

金沙江下游永善段隐蔽性滑坡隐患综合遥感识别

隐蔽性滑坡隐患是金沙江下游最普遍的地质灾害发育形式,具有隐蔽性强、突发性强、高位远程运动等特点。近年来,特大山区隐蔽性滑坡灾害案件频繁发生,对人民的生命财产造成了极大威胁。如何突破传统地质灾害调查手段的局限性、滞后性,提前有效识别隐蔽性滑坡隐患并探索其发育特征,对指导中国西南地区防灾减灾、工程规划建设具有重大科学意义。本文选择金沙江下游永善段地质灾害高易发区,利用升降轨时序InSAR–光学遥感综合识别方法,精细识别区域性时序地表形变、隐蔽性滑坡隐患光学遥感信息,通过野外考察,深入探索隐蔽性滑坡隐患发育特征。研究显示：1）通过升轨时序InSAR技术识别隐蔽性滑坡隐患26处,降轨时序InSAR技术识别隐蔽性滑坡隐患28处,光学遥感识别隐蔽性滑坡隐患48处（与升降轨时序InSAR识别结果有10处重合）,合计识别滑坡隐患92处;对识别结果进行100%的野外考察,将升降轨InSAR和光学遥感识别结果划分为完全一致、部分一致、仅有光学遥感识别结果、仅有InSAR识别变形结果4种情况,识别准确率分别为82.86%、80.77%、75.00%和63.64%,总体识别准确率达78.26%,略高于目前国内滑坡隐患识别平均水平,验证了滑坡隐患识别的可靠性和有效性。2）通过对比分析综合遥感识别结果可知,InSAR技术和光学遥感的识别结果与二者的识别方式、影像成像条件、滑坡活动性关系密切,二者不能直接进行互检。3）通过分析滑坡发育特征可知,隐蔽性滑坡隐患发育规律随着地形地貌、地质条件的变化而变化,升降轨InSAR技术和光学遥感识别的隐蔽性滑坡隐患在地貌空间分布、地层岩性均存在一定差别。结果表明,综合遥感识别技术充分利用了升降轨InSAR技术和光学遥感识别方法的互补性,解决了隐蔽性滑坡隐患看不见、看不清、看不准的难题,提高了滑坡识别的准确率。     

金沙江结合带巴塘段滑坡群InSAR探测识别与形变特征

金沙江结合带结构破碎,软弱岩层发育,流域性特大高位地质灾害频繁发生。针对该区域开展大范围滑坡调查与监测研究,对减灾防灾具有重要意义。以金沙江结合带巴塘段为试验区,采用堆叠InSAR技术分别利用升轨、降轨Sentinel-1A卫星数据对该区域滑坡隐患开展了调查研究。在此基础上,以中心绒乡滑坡群为重点研究区,利用多维小基线子集技术获取了区域二维形变速率(水平东西向和垂直向)及二维时间序列结果。通过对4处典型滑坡体的形变时间序列结果进行分析,发现在两年时间段内安里克米滑坡、仁娘村滑坡、贡伙村滑坡1和贡伙村滑坡2水平方向累积位移量分别达到44.3、-26.6、65.3和-77.1 mm,垂直向累积位移量分别达到-30.2、-88.3、-80.9和-56.9 mm,且这4处滑坡呈现缓慢蠕滑变形趋势。通过对贡伙村滑坡2的形变监测二维时间序列与降雨数据分析发现,强降雨对滑坡变形有一定短暂影响。由于滑坡群处于地质条件脆弱地区,构造活动强烈,在强降雨条件下极易导致滑坡失稳,建议对其进行持续监测,同时该研究成果对流域内其他区域的滑坡调查与监测研究具有参考意义。     

融合多源数据的稀土矿区滑坡危险性定量识别方法

为了提前对离子吸附型稀土矿区潜在的滑坡危险性进行识别,本研究以岭北矿区为例,提出了融合多源数据的稀土矿区滑坡危险性定量识别方法.首先基于时序哨兵一号(Sentinel-1A)数据,采用小基线集(Small Baseline Subset InSAR,SBAS-InSAR)技术获取研究区地表形变;再以数字高程模型(DEM)、Landsat 8陆地成像仪(OLI)等多源遥感数据为辅,提取潜在滑坡点,构建矿区滑坡危险性信息量模型;最后再结合同时期高空间分辨率遥感影像进行潜在滑坡危险性识别.结果表明：矿区整体地表年平均形变速率在-20.28～20.08 mm/a,共提取183个潜在滑坡点;构建的稀土矿区滑坡灾害危险性信息量模型是可行的,总结了研究区诱发滑坡的最佳信息量组合方式;对甲子背和大坑2个典型矿点进行具体分析,发现稀土开采活动会加速滑坡的产生,开采之后即使进行复垦,也容易诱发滑坡.     

汶川县地震诱发崩滑灾害影响因素的敏感性分析

震后遥感影像与调查结果表明,汶川Ms 8.0级地震在汶川县城约4 083 km2的区域内诱发了 5 528余处崩塌、滑坡灾害.利用GIS技术对地震地质灾害的分布与地层岩性、距断裂距离、坡度、高程等因素的关系进行统计分析,结合崩滑的确定性系数方法(CF),对诸影响因子进行敏感性分析,确定该区域内各因子最利于地震滑坡发育的数值区间,结果表明:地震崩塌滑坡主要分布在元古界、震旦系、二叠系硬质岩石分布的地层里,坡度为大于40°的坡体上;活动断裂10 km缓冲区内,未动断裂6 km缓冲区内崩滑密集分布;崩滑发育较有利的高程范围位于1 000～2 500 m;相对高差范围位于600～1 000 m之间,岷江干流对崩滑分布起绝对控制作用.分析得出影响本次崩塌滑坡分布的因子具有地震波传播介质效应和距离效应;由地震产生的崩滑与降雨诱发的崩滑在岩性、坡度的发生区间上有较为明显的区别;并为进行其他地区区域地震崩滑灾害稳定性评价提供了可以借鉴的科学依据.     

汶川地震山地灾害对环境因素的响应机制

为了研究汶川地震山地灾害灾害分布对环境因素的响应机制,通过对汶川地震灾区的调查和遥感解译,详细分析灾害分布特点,用GIS平台对4722个典型灾害点的环境因素进行分析,得出以下结论:1)地震能量释放的长短轴效应造成地震能量主要沿断裂带释放,灾害点沿中央断裂呈条带状分布,在距中央断裂20km范围内集中灾害总数的64.42%,逆冲断裂的上下盘效应使上盘灾害发育程度明显高于下盘和下盘灾害密度经历先降低后升高再降低的变化;2)灾害点在山地地区集中分布在河谷两岸,河谷岸坡对地震的放大作用使坡体上部岩土首先崩解,带动坡体下部的岩土一起运动,形成"山剥皮"特点的崩滑型灾害;3)灾害总数的99.43%集中分布在地震烈度Ⅶ度及其以上地区,地震的场地效应使灾害在易滑岩层中较为发育;4)在区域地势变化与地震断裂破裂一致的方向上有利于地震能量的传播,灾害更容易发生,规模较大的滑坡主滑方向也大多数发生在两者的并集方向.     

支撑拱式锁固型滑坡变形演化特征的物理模型试验

采用自行设计的滑坡物理模型系统,通过单、双段拉力筋牵引式启动的形式,分别设置了无锁固体和加装锁固体条件下单、双段滑坡试验(共4组),并利用应力-应变数据采集系统和高清摄像机等设备监测滑坡过程中坡体变形的多维表征信息。试验结果表明:在无锁固体条件下,边坡的牵引力-时间曲线峰值最大为39.5 N且快速陡降,坡体破坏形式为整体快速大滑动;在有锁固体条件下,边坡的牵引力-时间曲线峰值最大为55.7 N且维持近400 s,峰值提升41%,坡体的破坏形式转变为慢速蠕滑,可知拱支座间形成的支撑拱对拱后土体具有稳定作用。     

基于GIS的滑坡灾害危险性区划研究——以重庆市万州区为例

为了弥补滑坡灾害危险性区划研究中影响因子和等级划分的不确定性,结合前人研究成果,依据斜坡几何形态、岩性、地质构造、河流侵蚀、土地利用类型、人类工程活动、降水条件等影响因子与研究区实际已发生的滑坡灾害数之间的关系,编制重庆市万州区滑坡灾害危险性评价标准,并基于GIS技术和信息量模型法,计算滑坡评价因子的信息量,就万州区滑坡危险性进行区划,最后基于乡镇行政区对该区滑坡危险性区划进行细化。结果表明:建设用地、坡高为90～200m的地形、1024～1060mm的年降雨量以及侏罗系中统上沙溪庙组岩层等因素对万州区滑坡发生影响较大;根据滑坡灾害危险性评价标准,万州区滑坡灾害被划分为高、中、低、极低等4个危险区;应用信息量模型法得到的万州区滑坡危险性区划与实际情况比较吻合;高危险区和中危险区面积分别为564.4km²和848.6km²,分别占万州区总面积的16.3%和24.5%,主要分布于长江干流及支流两岸的居民相对集中区以及公路干线地段;高危险和中危险乡镇主要分布在万州区经济较为发达的长江干流两岸,尤其是左岸的黄柏乡、太龙镇、天城镇、李河镇等以及万州主城区。     

联合InSAR和质量守恒法估计西藏贡觉地区雄巴滑坡厚度

滑坡厚度是灾害风险评估和防控的关键参数。传统的滑坡厚度估计方法只能确定滑坡体局部稀疏点的深度,无法反映整个坡体深度分布情况,且存在探测成本高、效率低等问题。基于此,提出了一种联合合成孔径雷达干涉测量(InSAR)和质量守恒法估计平移式滑坡厚度分布的方法,并将其应用于西藏贡觉地区雄巴滑坡。使用时序InSAR技术分别计算Sentinel-1A升轨、Sentinel-1A降轨和ALOS-2升轨影像的一维雷达视线向(LOS)形变速度场,联合SAR成像参数与滑坡空间几何关系构建坡面坐标系,解算得到滑坡沿坡面的三维形变场; 在此基础上,联合InSAR三维形变场和质量守恒准则估计滑坡厚度。实验共收集了覆盖雄巴滑坡的140景Sentinel-1A升轨、138景Sentinel-1A降轨和12景ALOS-2升轨数据。结果表明:雄巴滑坡总运动面积约为5.33 km²,当流变参数取值0.5时,滑坡厚度为0～109.57 m,集中分布在9～73 m,滑坡体积约3.12×10⁸ m³,估算结果与现场测量结果一致。本文提出的方法可以获取滑坡连续的厚度分布情况,可为灾害风险评估、分析与预防提供关键性数据支撑。     

库水位升降作用下土质岸坡变形特征实验研究

库水位升降作用下三峡库区土质岸坡坡体吸水、应力集中及滑移变形,形成不同时段和不同空间部位的裂缝体系。通过模型试验的方法,分析了土质岸坡在一个蓄水降水循环周期内裂缝体系的时空演化分期配套规律。试验结果表明：蓄水初期,裂缝主要出现在岸坡前缘水位线附近,且出现频率高、规模逐渐增大,岸坡出现局部坍塌破坏;蓄水中期,岸坡前后缘均有裂缝出现,但出现频率低、规模减小;蓄水后期,前缘裂缝发育基本消失,后缘裂缝继续发育,岸坡以沿滑动面整体蠕滑为主;水位上升过程中暂停蓄水时裂缝发展缓慢。降水阶段,拉张裂缝主要出现在岸坡坡体后缘,初期降水裂缝出现频率较低,发展缓慢;后期降水裂缝出现频率、规模变大,产生下座变形及下错台坎;水位下降阶段岸坡沿滑动面发生整体缓慢蠕滑。     

汶川地震诱发滑坡方向效应研究

针对在野外调查过程中发现汶川地震灾区的一些特殊区域,滑坡的滑动方向呈现出一些与地震动力直接相关的异常现象,作者对典型研究区进行了详细的遥感解译,统计了不同坡向上滑坡灾害的发育密度和大型滑坡滑向与发震断层的关系,得出汶川地震诱发滑坡具有"背坡面效应"和"断层错动方向效应"."背坡面效应"是指在与发震断裂带近于垂直的沟谷斜坡中,在地震波传播的背坡面一侧的滑坡发育密度明显大于迎坡面一侧."断层错动方向效应"是指滑坡的滑动方向与发震断裂的活动方式和错动方向密切相关.在分析强震诱发滑坡成因机理时,应加强对2种因素和效应的研究:1)地震时断裂的突然错动所产生的直接动力对斜坡失稳的初始驱动作用;2)地震过程中地震波在斜坡山体内传播的折射、反射以及由此引起的坡体内应力的变化.     

滑坡稳定性的束口锁固控制效应研究

为研究滑坡稳定性的束口锁固控制效应,以豫西某典型束口型滑坡为例,建立了该滑坡的三维数值计算模型,利用大型岩土数值分析软件FLAC 3D对该滑坡在天然和饱和工况下的变形和破坏特征进行了计算分析。结果表明:天然状态下,该滑坡以发生在滑体中的局部变形为主;饱和状态下,滑坡滑体呈现出由于滑带破坏而整体变形的破坏特征;受地形特征的影响,滑坡束口处的锁固控制效应明显。对束口锁固控制效应的研究可从束口前缘的阻滑效应和大主应力的拱形分布特征两个角度进行分析。分析发现束口对滑坡具有显著的阻滑效应,表现为滑体前缘和后缘变形小、中部变形大的特征;束口锁固段大主应力呈拱形分布,作用模式为土体内部应力集中的转移与重分布,可为坡体的整体稳定起到控制作用。     

堰塞金沙江上游的白格滑坡形成机制与过程分析

2018年10月10日和11月3日,西藏自治区江达县白格村金沙江右岸先后2次发生滑坡堵江事件,堰塞湖与溃坝洪水给金沙江上游沿岸居民及其生产和生活设施带来巨大灾害。滑坡发生后,作者先后2次赶赴现场,参与灾害调查与救灾工作。基于现场调查,结合相关资料,对滑坡的形成机制与过程进行系统分析。结果表明:1)白格"10·10"滑坡是一个高位、高剪出口、高速非完全楔形体基岩滑坡,方量约107 m~3。2)滑坡地处金沙江缝合带,岩性为元古界熊松群片麻岩组,具有多期、多次变形与变质特点,糜棱岩化和蚀变很严重。3)滑坡按高程划分为3区,即前缘的阻滑区、中部的主滑区和后缘的牵引区,分割高程大致为3 500和3 000 m。主滑区为楔形体,系2组发育良好的结构面切割形成;阻滑区为四面体,由2组发育较差的结构面切割形成;牵引区为完全风化的岩土体夹团块状碎裂岩体。4)滑坡存在2个滑动方向,即主滑区的S80°E方向和阻滑区的N70°E方向,剪出口高程约2 950 m。5)主滑区楔形体重力是滑坡的主要动力来源,滑坡的孕育过程是相对完整的阻滑区岩体在主滑区重力驱动下的渐进破坏过程。6)滑坡过程如下:首先,主滑区和阻滑区启动;其次,失去支撑的牵引区再启动;随后,先启动的滑体高速撞击四川岸,逆坡爬高约95 m,并在两侧形成碎屑冲刷区;然后,先启动的滑体折返,并与后启动的滑体在河面上方相撞,冲击河水形成高速水砂射流,在两岸形成水砂射流冲刷区;而后,堰塞坝下游坡滑动,形成次级滑移区;最后,冲击产生的雨雾降落,完成滑坡坝表面冲刷。7)白格"11·3"滑坡是牵引区的部分岩土体在起阻滑作用的碎裂岩体渐进解体后下滑的结果,方量约3×106 m~3。8)牵引区目前严重变形的方量约5.50×106 m~3,存在再次滑坡与堵江的风险,需要采取合理的工程措施消除隐患。     

地下水雍高诱发黄土滑坡离心模型试验研究

为了研究地下水位雍高条件下坡体的变形破坏特征和失稳的机理,以黑方台焦家滑坡为地质原型。通过在离心模型后部安置排水挡板和临时水箱,其能较好得模拟灌溉水通过集中通道快速入渗补给坡体地下水,再现地下水位雍高诱发黄土滑坡变形破坏的全过程。此外,在离心模型试验期间,对边坡的孔隙水压力、土压力和变形演化特征进行实时监测和研究。研究表明:地下水雍高造成底部黄土层软化,潜蚀作用和不均匀沉降致使边坡产生塑性流动,从而导致滑坡的顶部形成不同程度的平行错拉裂缝。造成这类滑坡的后壁形态大多呈现直立状。地下水位雍高造成坡体内土体的局部剪缩变形导致孔隙水压力持续以及急剧增高,最终诱发此类滑坡发生突发性破坏,呈现出高速泥流破坏特征。因此地下水位雍高诱发滑坡的形成机理过程为:牵引式蠕滑-半坡以及后缘拉裂-锁固段累进性滑移剪断-滑带贯通高速溃滑破坏。     

基于地球多传感器网络信息的潜在滑坡判识模型

潜在滑坡的判识是滑坡预测预报中的关键问题,潜在滑坡早期判识能够有效减少灾害的发生。基于滑坡灾害诱发的各种影响因素,利用地球多源时空信息和多传感器网络,遥感监测潜在滑坡体,利用监测设备获取滑坡岩性、坡体结构、地貌形态、活动迹象等关键控制信息,从中选取滑坡灾害诱发的主要控制要素作为判识指标,建立基于多源信息的潜在滑坡多因素判识模型;特别是根据滑坡的成灾规律,分析潜在滑坡孕育过程中地貌形态改变、成灾条件变化与滑坡发生的临滑诱发条件,建立基于不同信息源的滑坡控制因素判识模型。最后,通过汶川地震灾区垮梁子滑坡体实例分析,验证潜在滑坡综合判识模型并进行优化,从而为滑坡灾害的早期判识预测提供依据。     

重大滑坡灾害全过程调控减灾技术及试验验证与潜在应用

青藏高原其独特的地质背景和复杂的地理环境,滑坡灾害极为发育。由于构造活动强烈、河谷深切、斜坡高陡,内外动力耦合作用下诱发的高位崩滑往往具有极高的位能和超远的位移,对其危害范围内的构筑物和人类活动等造成严重破坏。因此,青藏高原重大滑坡灾害防治应从源头区—运动区—堆积区进行全过程调控和综合治理,因势利导地进行总体设计。本文对重大滑坡灾害源区、运动区的处置措施、减灾机理及其应用范围进行了探讨,根据青藏高原滑坡特点,研发了适应高海拔地区的变管径高扬程虹吸排水技术,抵抗滑体大变形的新型抗滑桩结构、以及抗冲击耗能减震棚洞结构。室内虹吸排水试验和混凝土梁试样的四点弯试验结果表明：1)对于不同扬程的试验,4 mm + 6.5 mm的虹吸系统排水能力均强于4 mm的,变管径使得排水能力提升比例为64%~117%;2)随着扬程的升高,6.5 mm虹吸排水系统排水能力下降明显比4 mm + 6.5 mm快,变管径使得随扬程升高排水能力下降比例减小了28%~42%;3)相较于传统抗滑桩,新型钢绞线抗滑桩不仅提高了抗折强度,而且韧性更好,可承载更大的变形而不发生破坏。以金沙水电站山梁子高位崩塌滑坡灾害治理为例,开展了从滑源区到运动区全过程的综合治理设计和应用实施,为金沙水电工程的安全施工与长期运行提供保障。     

冰川密林陡谷区活动性滑坡InSAR精细识别-以澜沧江梅里雪山段为例

西南山区地质灾害频发,近年来,连续多起如茂县滑坡、白格滑坡、水城滑坡等大型滑坡灾害表明,灾难型滑坡事件的发生不仅与该地区易孕育地质灾害的条件相关,无法准确识别与监控滑坡灾害的发育发展情况也是主要因素之一。位于青藏高原东南缘的“三江”地区,地形高陡、地质灾害频发,其中澜沧江流域环境更为恶劣,海拔高、地形陡,夏秋季山高林密、多降雨,春冬季天气寒冷、多降雪,传统的地面式地质灾害调查手段在此效果有限,光学遥感技术也受到天气条件的诸多限制。而基于主动式微波遥感成像的InSAR技术被证明适合进行人类不易到达地区的地质灾害识别工作。本文以澜沧江流域梅里雪山段为研究区,使用PALSAR1、PALSAR2与Sentinel-1多源雷达遥感数据,结合D-InSAR与IPTA-InSAR技术开展冰川密林陡谷区的滑坡灾害精细识别工作。过程中发现长波长的PALSAR1和PALSAR2数据有更好的适应性,短波长的Sentinel-1数据通过数据量多的优势弥补了自身穿透能不足的缺陷；成功解译了研究区内92处活动性滑坡灾害；平面面积超过20万m_²的有26处；时间序列InSAR结果结合地面调查验证表明,澜沧江河谷典型的大型滑坡,梅里水滑坡和亚贡滑坡活动速率快,危险性高；德钦县成周边典型滑坡,丁羊丁滑坡与归巴顶滑坡整体处于蠕变状态,局部变形速率高；同时德钦县城周边泥石流与冰碛物运动较为活跃。     

猴子岩水库色玉滑坡涌浪灾害链CFD-DEM耦合数值模拟

滑坡涌浪是一种常见的地质灾害现象,由于滑坡体与水体之间存在复杂的流固耦合作用,使得传统的单一介质模型无法进行准确求解。为此,介绍一种基于计算流体力学方法（CFD）与离散单元法（DEM）的流固耦合模型CFD–DEM,采用计算流体力学方法（CFD）求解水体流动,采用离散单元法（DEM）模拟散粒体滑坡运动,充分利用不同计算模型的优势,对滑坡及涌浪演进过程进行数值模拟分析。首先,利用该耦合模型对Robbe–Saule开展的颗粒堆积体坍塌–涌浪试验进行了相同工况下的数值计算,从颗粒坍塌运动过程、涌浪高度演化过程等方面进行对比,结果表明模拟结果与试验结果吻合很好,验证了CFD–DEM流固耦合模型的有效性。然后,将该方法应用于四川省猴子岩水库色玉滑坡–涌浪灾害的演进过程分析,重现了该事件滑坡失稳运动、涌浪产生及传播、涌浪爬升、涌浪回流的全过程,计算结果显示：计算得到的电站进水口处涌浪高度与实测数据较为接近;色玉滑坡从失稳运动至静止堆积的持续时间约为20 s,颗粒平均速度最大达到16.12 m/s;滑坡引起的涌浪约在滑坡失稳10 s后传播到对岸,之后开始沿坡面向上爬升,最大爬升高度达到27.32 m。研究表明CFD–DEM流固耦合模型能够很好地应用于模拟山区河谷大规模滑坡涌浪灾害,可为库区防灾减灾提供高效的技术支持。