金沙江上游超大古堰塞湖及其相关问题

历史大型堵江滑坡研究不仅是认识区域构造背景和河流演变规律的重要工作,而且对山区河流工程安全亦有十分重要的意义。以金沙江上游河段的王大龙古滑坡为研究对象,基于野外调查与资料分析,对其成因与演变过程进行了分析。主要结论如下：（1）王大龙滑坡位于金沙江缝合带,源区基岩主要为三叠系中心绒群下段板岩 (T1-2zh1),其次为二叠系嘎金雪山群下段石英砂岩 (Pgj1),滑坡总体上为王大龙断裂与中心绒群板理面切割形成的楔形体,前缘受雄松-苏洼龙活动断裂切割,方量约4.0×108m3;（2）诱发滑坡的内动力为地震,外动力为河流凹岸侵蚀,滑坡时间应为晚更新世大理冰期;（3）堰塞坝长约1700m,宽约3000m,高度超过450m,形态右高左低,右岸高程约2770m,左岸垭口高程约2735m;（4）堰塞湖规模约266×108m3,干流库尾到达叶巴滩水电站坝址的降曲河口,长度约176 km;（5）堰塞坝发生过3次溃决,溃口底面高程分别为2460m、2400m和2358m（现河面高程）;（6）第一次溃决极有可能雄松-苏洼龙断裂错动导致的,时间早于1900 a BP,估算溃口流量21.0~39.8×104m3/s,远远大于长江历史洪水记录。     

金沙江白格滑坡裂缝区失稳概率分析

金沙江上游白格滑坡于2018年10月10日和11月3日连续发生两次滑坡堵江-溃坝事件,给沿岸居民带来巨大灾害。滑坡发生之后,在滑坡体后缘形成3个裂缝区,存在再次滑坡堵江的可能。鉴于白格滑坡岩土体试验数量有限,试验成果不能反映坡体抗剪强度参数的统计特征,本文针对白格11.03滑坡,采用基于极大似然估计的优化算法,对滑坡岩土体的抗剪强度参数进行了反演,并利用反演参数,采用蒙特卡洛方法对滑坡后缘裂缝区进行了可靠度计算和稳定性评价。研究表明：（1）将响应面法与基于极大似然的概率反分析方法结合起来对滑坡体参数进行反演,计算过程简便,计算方法适应性强。（2）根据白格11.03滑坡反演得到的滑坡岩土体抗剪强度后验分布为c~N(13.05kPa, 32),tanφ~N(0.64, 0.042),计算结果表明通过反演可以降低先验参数的不确定性。（3）采用自动搜索最危险滑动面的方法得到了各裂缝区的局部破坏模式。计算得到的最危险滑动区范围与现场调查及监测得到的失稳区范围基本一致,且局部破坏模式下的破坏概率均明显大于整体破坏模式,表明裂缝区破坏以渐进牵引破坏为主,与实际的裂缝发展规律及破坏模式一致。（4）白格滑坡后缘3个裂缝区的潜在失稳区域C1-1、C2-1和C3-1的体积分别为7×105 m3、3.2×106 m3和1.3×106 m3,对应的失稳概率分别为21.70%、33.90%和27.30%,均为极高危险性。因此,有必要对滑坡堵江-溃坝洪水灾害链进行风险评价,并提出合理的处置措施和应急预案。     

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

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,存在再次滑坡与堵江的风险,需要采取合理的工程措施消除隐患。     

金沙江白格堰塞湖溃决过程数值模拟

2018年10月10日和11月3日，在我国四川省与西藏自治区交界处的白格村同一位置连续发生两次滑坡，完全堵塞了金沙江形成堰塞湖。由于“10?10”滑坡形成的堰塞湖水位抬升迅速，堰塞湖于10月12日自然漫顶溃决。“11?03”滑坡堵塞了“10?10”堰塞体溃决形成的流道，形成了更大的堰塞湖，鉴于客观条件允许，采取了开挖泄流槽降低堰塞湖溃决水位的措施，至11月12日，堰塞湖发生漫顶溃决，溃口洪水峰值流量为31000m3/s。由于“11?03”白格堰塞湖溃决案例拥有较为完整的实测资料，为堰塞湖溃决过程的研究提供了宝贵的基础数据。基于堰塞体的溃决机理，建立了可考虑堰塞湖的水动力条件、堰塞体的形态和材料特征的堰塞湖溃决过程数学模型。模型采用宽顶堰公式模拟溃口洪水流量，并根据堰塞湖入湖和溃口流量以及堰塞湖的水位-湖面面积关系曲线确定堰塞湖水位的变化；采用基于水流剪应力和堰塞体材料临界剪应力，并可考虑宽级配堰塞体材料特性的冲蚀公式模拟材料的冲蚀过程；假设溃口在纵向下切和横向展宽过程中坡角保持不变，采用极限平衡法分析溃口在发展过程中发生的边坡失稳；采用按时间步长迭代的数值计算方法模拟堰塞湖溃决时的水土耦合过程。采用建立的模型对“11?03”白格堰塞湖溃决案例进行反演分析后发现，模型计算获得的溃口流量过程、堰塞湖水位变化过程、溃口发展过程与实测值基本吻合。参数敏感性分析表明，冲蚀系数对溃决过程具有重要的影响，残留坝高通过影响下泄库容也对溃决过程产生作用；另外，开挖泄流槽可降低堰塞湖溃决时的库容，从而对溃口流量过程产生影响，是降低灾害损失的有效手段。     

白格滑坡裂缝区演变过程及其发展趋势分析

白格“10.10”滑坡后,其源区边界外残留有方量较大的裂缝区,“11.3”滑坡就是裂缝区部分失稳的结果。“11.3”滑坡后裂缝区范围进一步扩大,存在进一步失稳,进而滑坡堵江的可能性。为此,2019年作者对滑坡区进行了补充地质调查,并开展了裂缝区深部变形与渗压监测等工作。结果表明：(1)滑坡区岩性属于金沙江缝合带的构造混杂岩,以片麻岩为主,夹杂大理岩、碳质板岩,以及超基性和酸性侵入岩。虽然构造活动扰动强烈,边坡整体为逆向坡。蚀变和风化作用对滑坡的孕育演化作用显著。(2)自“11.3”滑坡发生至2019年6月,滑坡裂缝区范围一直在扩展。C1区扩展约500m,到达白格村;C2区越过山脊,至滑坡边界直线距离约130m;C3区范围相对稳定。从宏观变形来看,“11.3”滑坡诱发的裂缝变形最大,此后发展趋缓。(3)C1-1、C2-1和C3区均已形成剪切带,变形深度C1-1区最小(18m),C2-1区最大(67m),C3-1区居中(22.5m)。平均变形速率C1-1区最大,推测大于20mm/天;C3-1区其次,达到3.33mm/天;C2-1区最小,为0.29mm/天。(4)裂缝区不存在稳定的地下水位,变形与降雨无关。(5)裂缝区目前变形均未收敛,均处于蠕变变形状态。后期需要密切关注C1-1和C2-1区,特别是变形特征与局部基岩解体破坏情况。     

八宿巨型滑坡的发现及其意义（特约稿）

巨型滑坡全球不常见,但是破坏性极强,其识别对于加深区域构造背景的认识,防灾减灾和工程建设均有指导意义。基于遥感解译、无人机地形测绘和现场调查,报道了新发现的八宿巨型滑坡及其残留堆积体特征,初步分析了成因机制与演变过程,估算了滑坡方量。主要结论如下：（1）八宿滑坡发育于怒江缝合带的夏里-八宿裂谷带内,岩性以古生界加玉桥岩群大理岩组和侏罗系马里组变质砂泥岩为主,上部岩层产状为陡倾顺坡向,下部岩层产状为陡倾逆坡向,岸坡结构总体上为上硬下软;（2）滑坡区横跨冷曲的左右两岸,长约7200m,宽约4800m,面积约22.5km2,估算方量约35×108m3,目前残余方量约14×108m3,是国内已经确认的方量最大的滑坡;（3）滑坡极可能属于历史地震诱发的高速堵江滑坡,堰塞坝高度约185m,堆积体碎裂化严重,且在冷曲右岸爬升超过600m。滑坡形成后,堆积区历经了堰塞坝溃决、多拉寺次级滑坡、泥石流堆积与冲刷、表面流水冲刷等改造过程,但是滑坡地貌特征整体保存良好。由于滑坡和泥石流堵江,冷曲在滑坡区先后两次改道;（4）从滑坡与冷曲的演变过程来看,推测滑坡的发生时间应在晚更新世。即八宿滑坡为古滑坡;（5）八宿巨型滑坡的发现说明迄今为止我们对青藏高原地质环境的认知仍然十分有限。川藏铁路等工程建设应加强基础地质与工程地质研究,以及地质灾害风险论证,确保施工与运行期安全。     

湖水位壅高条件下堰塞体土料变形特性DEM-PFV耦合模拟研究

我国西南地区多为高山峡谷地貌,易发生滑坡堵江事件,形成堰塞湖。堰塞湖水位壅高过程中,堰塞体为常剪应力路径,即剪应力保持不变,孔隙水压力不断增大。但已有的研究主要集中于固结排水剪和固结不排水剪,与堰塞湖水位壅高过程中的实际应力路径有别,因此开展常剪应力路径下堰塞体材料的变形特性研究,是深入分析湖水位壅高过程中堰塞体动态响应的有效途径。鉴于此,本文以2018年10月11日白格堰塞湖为例,以堰塞体的实际高度、漫顶前的最高堰塞湖水位及湖水位壅高过程中的实际应力路径为基础,基于细观尺度的离散元（DEM）—孔隙有限体积法（PFV）流固耦合方法,从敏感性分析的角度出发,开展了不同围压、不同初始应力条件下堰塞体土料的常剪应力剪数值模拟试验,并从材料的应力应变关系（宏观）及内部接触力（微观）的分布规律等角度出发,揭示了湖水位壅高过程中堰塞体不同位置的变形响应及其微观力学机理。研究表明,湖水位壅高条件下堰塞体土料变形特性受到土料位置、强度和湖水位壅高程度的联合影响。处于堰塞体不同位置的土料,围压与初始应力比条件不同,并且在堰塞体漫顶之前所遭遇的最大孔隙水压力也不同,从而导致在堰塞湖水位壅高过程中,不同位置的堰塞体土料呈现出不同的变形特性,一般呈现由里及外变形逐渐增大的规律。在相同围压条件下,靠近堰塞体上游外缘的土料,初始应力比相对较高,且遭遇的最大孔隙水压力也相对较高,从而在堰塞湖水位壅高过程中其应力路径会穿越失稳线,导致颗粒之间的接触力减弱,从而产生较大的变形,且大变形区的厚度与范围受到初始应力比及最高湖水位的限制。堰塞体内部初始应力比相对外缘较少,在湖水位壅高过程中应力状态穿过失稳线的可能性降低,从而变形也相对较小。     

湖水位壅高条件下堰塞体土料变形特性DEM–PFV耦合模拟研究

中国西南地区多为高山峡谷地貌,易发生滑坡堵江事件,形成堰塞湖。堰塞湖水位壅高过程中,堰塞体为常剪应力路径,即剪应力保持不变,孔隙水压力不断增大。但已有的研究主要集中于固结排水剪和固结不排水剪,与堰塞湖水位壅高过程中的实际应力路径有别,因此开展常剪应力路径下堰塞体材料的变形特性研究,是深入分析湖水位壅高过程中堰塞体动态响应的有效途径。鉴于此,本文以2018年10月11日白格堰塞湖为例,以堰塞体的实际高度、漫顶前的最高堰塞湖水位及湖水位壅高过程中的实际应力路径为基础,基于细观尺度的离散元(DEM)–孔隙有限体积法(PFV)流固耦合方法,从敏感性分析的角度出发,开展了不同围压、不同初始应力条件下堰塞体土料的常剪应力剪数值模拟试验,并从材料的应力应变关系(宏观)及内部接触力(微观)的分布规律等角度出发,揭示了湖水位壅高过程中堰塞体不同位置的变形响应及其微观力学机理。研究表明,湖水位壅高条件下堰塞体土料变形特性受到土料位置、强度和湖水位壅高程度的联合影响。处于堰塞体不同位置的土料,围压与初始应力比条件不同,并且在堰塞体漫顶之前所遭遇的最大孔隙水压力也不同,从而导致在堰塞湖水位壅高过程中,不同位置的堰塞体土料呈现出不同的变形特性,一般呈现由里及外变形逐渐增大的规律。在相同围压条件下,靠近堰塞体上游外缘的土料,初始应力比相对较高,且遭遇的最大孔隙水压力也相对较高,从而在堰塞湖水位壅高过程中其应力路径会穿越失稳线,导致颗粒之间的接触力减弱,从而产生较大的变形,且大变形区的厚度与范围受到初始应力比及最高湖水位的限制。堰塞体内部初始应力比相对外缘较少,在湖水位壅高过程中应力状态穿过失稳线的可能性降低,从而变形也相对较小。     

地震堰塞湖排险技术与治理保护

“5·12”汶川特大地震形成了上百处堰塞体高度大于10 m,蓄水量大于1.0×105 m3,集雨面积大于20 km2的堰塞湖。根据各类堰塞湖因对下游威胁程度、地质、地貌、水文条件、堰塞体规模和颗粒组成制定相应的排险方案。极高危和高危堰塞湖具有潜在的溃坝风险,必须及时排险。相对稳定的堰塞湖及其所在河道则应治理保护。结合北川唐家山、都江堰枷担湾等堰塞湖排险实例,提出了地震堰塞湖排险技术,包括险情勘查、险情评估、排险方案、施工技术。同时,对灾后堰塞湖观测、治理与保护利用提出了建议。     

金沙江白格滑坡堰塞坝溃决洪水灾害调查与致灾浅析

2018年10月10日和11月3日,西藏自治区江达县波罗乡白格村金沙江右岸同一位置先后两次发生滑坡堵江事件并形成了巨大的堰塞湖,其堰塞坝在自然泄流和人工开挖泄流槽两种处置方式后溃决。其中,第二次滑坡堰塞坝的溃决洪水给下游西藏、四川和云南3省（自治区）受灾范围内的道路、桥梁、耕地和房屋造成了巨大破坏。为了应对类似的极端、超常规、特大堰塞坝溃决洪水威胁以及相关的基础性研究需要,课题组于2018年12月21日至29日对这次金沙江白格堰塞湖溃坝洪水对下游的受灾情况进行了考察调研。考察以受灾最为严重的巴塘县巴楚河（又称巴曲河）与金沙江的交汇口为起点,直至洪水威胁基本消除的梨园水库库区为终点沿江共计488.6km的受灾河段为主。考察重点为沿岸房屋、道路、桥梁和水利基础设施等受损情况,并对溃坝洪水的最大淹没水位（洪痕）、考察时的河道水位,河道两岸堆积的泥沙及其颗粒级配,桥梁致灾水位等进行了分析,得到了一些有价值的灾情数据与成果,这些成果可为进一步的基础性研究提供一定的数据支撑。     

“10·10”西藏白格滑坡运动特征反演分析

为获取“10·10”白格滑坡运动规律与特征，基于“10·10”白格滑坡发育的3大区（启动区、堆积区、冲击区）和6个运动阶段（主体失稳破坏、牵引区启动、高速临空滑跃、冲击对岸、折返相撞及水砂射流、堆积坝次级滑移）,对滑坡碎屑体堆积状态及其在四川岸的冲击形态进行分析。采用谢徳格尔法及能量转化计算方法，选取滑坡体冲击高度的5个特征点，计算滑坡碎屑体在各运动阶段的运动速度。结果表明：主滑区阻滑区滑体以2.2 m/s的初速度启动，从启动区至剪出口，速度不断增大，5个特征点达到最大速度，分别为H₁ 67.0 m/s、H₂ 73.0 m/s、H₃ 73.7 m/s、H₄ 73.2 m/s、H₅ 68.3 m/s;牵引区滑体到达剪出口时，速度为70.2 m/s;主滑区阻滑区滑体的滑动速度具有从中间向两端递减的态势，其中滑坡主滑方向的中间位置速度最高，达73.7 m/s;滑体整体剧动时所释放的能量E至少为10^10.8J,引起震动相当于4.0～4.7级的表层基岩地震。通过滑坡...     

金沙江上游白格滑坡黏土化蚀变岩的灾变力学行为研究

金沙江构造结合带是典型的板块碰撞作用形成的构造混杂岩带,不仅地质结构复杂,而且发育了工程特性极差的特殊岩土体。在金沙江上游白格滑坡成因机理调查与分析过程中,新发现了巨型滑坡体后缘滑壁发育厚度较大、分布较连续的黏土化蚀变岩,是构成滑动带的重要组分。基于相关样品的原位测试和室内系统的物质组成、耐崩解、环形剪切、动三轴等物理力学试验研究,揭示了黏土化蚀变岩的灾变力学特性及其在白格滑坡形成演化中的作用。主要结论如下：（1）白格滑坡黏土化蚀变岩是金沙江构造结合带中蛇绿岩蚀变的产物,蚀变程度高,蚀变系数介于0.64～0.86之间,黏土矿物成分以伊利石为主、伊/蒙混层矿物次之,遇水极易崩解,崩解指数达50%以上;（2）黏土化蚀变岩中工程性质最差的是蚀变黏土,水和振动荷载作用下强度弱化显著,含水率由10%增至20%后剪切强度降幅达34%,动强度远小于其静强度,并随着振动次数的增加不断降低,是促进斜坡岩体结构面贯通、破坏的重要因素;（3）白格滑坡滑动面的形成受节理结构面与黏土化蚀变岩联合控制,整体上可分为2种模式：黏土化蚀变岩厚层软弱夹层型和充填节理裂隙贯通型;（4）金沙江上游发育的不少巨型高位滑坡与黏土化蚀变岩不良工程特性具有密切的关系,岩体结构与特殊岩性的组合是值得关注的区域性易滑地质结构。相关认识对于青藏高原东缘大江大河岸坡稳定性分析和防灾减灾具有较好的启示意义。     

金沙江结合带巴塘段滑坡群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的形变监测二维时间序列与降雨数据分析发现,强降雨对滑坡变形有一定短暂影响。由于滑坡群处于地质条件脆弱地区,构造活动强烈,在强降雨条件下极易导致滑坡失稳,建议对其进行持续监测,同时该研究成果对流域内其他区域的滑坡调查与监测研究具有参考意义。     

G4216仁沐新高速公路新市互通区巨型滑坡处治方案

G4216线仁寿经沐川至新市高速公路止点发育多处巨型滑坡,且新市镇作为仁沐新高速、宜新高速、新金高速、串新高速等4条高速公路转换的重要结点,有设落地互通的需求;因此给新市互通区的工程防护、处治设计提出了巨大挑战。在室内试验结合工程现状反算的基础上综合确定滑带土力学参数,再通过滑坡剩余推力计算,得到主滑断面最大剩余下滑力为16 720 kN/m,远超支挡结构提供的极限抗滑强度。结合互通区构造物设计、路基标高、滑坡推力等因素确定了分区处治方案:清方减载实现整体稳定;前缘设桩加固次级滑面;后缘设锚杆(索)加固清方形成的高边坡;采用小截面抗滑桩支挡顶部中小型滑坡;结合坡面截排水措施共同防护。     

库水复活型滑坡牵引滑动机理模型试验研究

为深入探索库水复活型滑坡的破坏模式和力学机制,设计研发新型试验装置,装置主体是由若干渗透盒构成的分段式滑面,能够模拟各种几何形态的滑动面。通过向不同区段的渗透盒注水,可模拟滑带土分阶段饱水软化,实现库水位上升条件下相邻滑块间的牵引滑动过程。采用传感元件、数码摄像以及数字图像处理等测试手段,获取滑带土体积含水率、孔隙水压力、坡面水平位移的变化规律以及后缘面演化特征,并探讨滑坡体失稳破坏模式。结果表明:滑带土饱水后抗剪强度大幅度下降是滑坡发生的重要条件,滑动面处的孔压变动是激发老滑坡复活的重要诱因;滑带处孔隙水压力的增大与滑坡位移的增大是同时发生的,滑带土强度衰减与滑体变形具有良好的相关关系;坡面变形区域分为强、弱变形区和牵引区3部分,坡面变形为1～1.5倍的失稳滑带长度;随着失稳滑带长度增加,坡面变形区域变大,对后侧稳定坡体的牵引变形影响变小;后缘破裂面多呈折线型滑面形态,后缘破裂面倾角试验值受失稳滑段位置和滑体厚度影响显著。研究结果为库区滑坡灾变机制的认识和治理提供重要依据。     

滑坡–堵江–涌浪灾害链模拟的DEM–CFD耦合分析方法及其应用

滑坡-堵江-涌浪灾害链的演化过程涉及复杂的滑坡-河流相互作用。为了真实再现这一演进过程,本文提出了一种扩展的离散元法（DEM）与计算流体动力学（CFD）耦合数值方法,通过在局部平均DEM-CFD耦合方法中引入流体体积法（VOF）建立了一套能够描述河流自由水面演化的扩展DEM-CFD耦合数值框架,实现了对滑坡涌浪形成及传播过程的追踪;并提出了一种新的局部孔隙率计算方法即虚拟球模型,来克服网格颗粒临界尺寸比的限制,从而实现复杂地形建模要求;进一步通过多个典型算例验证其准确性和有效性。在此基础上,开展了滑坡-堵江-涌浪灾害链演化模拟,并应用至白格滑坡坝形成过程复演,深入分析了滑坡-河流动态演化过程及耦合作用机理。结果表明：所提出的方法能够准确模拟滑坡-河流相互作用;滑坡-河流-河谷地形的耦合作用驱动了灾害演化,滑坡体驱动了涌浪的传播,而河流则增加了滑坡体的动能耗散,且对其沉积过程在河流流向上产生相反影响;模拟得到的白格滑坡运移路径、沉积形态及涌浪侵蚀面积与实地调查结果吻合度较高,成功再现了2018年10月11日白格滑坡堵江事件。本文所提出的扩展DEM-CFD耦合数值方法为深入研究滑坡-堵江-涌浪灾害链演化机制提供了有力的数值工具,对于灾害预测及防灾减灾策略的制定具有一定的参考价值。