基于溃决机理的堰塞湖溃决快速风险评估方法研究.pdf

第 卷 第 期 年 月自 然 灾 害 学 报 .收稿日期:修回日期:基金项目:国家自然科学基金项目()作者简介:王 琳()女副教授博士主要从事溃决机理及灾害链效应研究:.通讯作者:钟启明()男教授级高级工程师博士主要从事溃决机理及灾害链效应研究:.文章编号:()./.基于溃决机理的堰塞湖溃决快速风险评估方法研究王 琳苑鹏飞钟启明胡 亮单熠博薛一峰(.西安理工大学 省部共建西北旱区生态水利国家重点实验室陕西 西安 .南京水利科学研究院江苏 南京.水利部土石坝破坏机理与防控技术重点实验室江苏 南京 .陕西省水利电力勘测设计研究院 陕西 西安)摘 要:堰塞湖发生后极易形成溃决灾害链亟需构建基于有限数据包的溃决快速定量风险评估方法 通过快速获取堰塞湖影响区域的三维地理信息构建溃决致灾的快速定量评估模型 基于堰塞湖坝体颗粒级配组成实现精细化、简单化的稳定性快速评价基于冲蚀特性和崩塌过程溃决机理实现溃口流量变化过程分析与洪水演进过程模拟(内)基于极限学习机网络模型建立风险人口与生命损失的函数实现了生命损失评估预警明晰了溃决致灾的灾害链效应 将其应用于唐家山堰塞湖实例发现:研究方法能够较好地预警溃决灾害链开挖引流槽可降低堰塞湖的溃口峰值流量、最大流速、溃口宽度和溃决库容但无法防止溃决发生若不开挖引流槽溃决库容将达到.溃口峰值流量达到./溃口顶宽增大到.开挖引流槽可使溃口峰值流量减少.溃决库容减少.降低了堰塞湖的溃决风险当提前预警时间超过 后及时疏散下游淹没范围内的居民可使生命损失率降低为 研究方法可实现堰塞湖应急处置时的快速定量风险评估并为其应急处置决策提供技术支撑关键词:堰塞湖稳定性溃决过程生命损失风险评估唐家山中图分类号:文献标识码:(.):.().自 然 灾 害 学 报第 卷 ./.:引言堰塞湖是受地震活动、气候变化、人类活动等影响的一种典型地质灾害 据统计的堰塞湖在 内溃决在 内溃决在 个月内溃决 其一旦溃决将形成溃决致灾灾害链威胁下游沿岸的工程设施及人民群众生命财产安全 目前我国堰塞湖溃决灾害频发例如:年 月易贡堰塞湖形成 月溃决溃决时峰值流量达到 /年最具代表性的唐家山和小岗剑堰塞湖溃决洪峰分别达到、/使接近百万人的生命财产安全受到威胁 年白格堰塞湖更是 次堵塞金沙江上游干流溃决洪峰流量分别达到、/月 日形成的堰塞湖更是在 内就发生溃决 月 日堰塞湖溃决导致下游在建的苏洼龙水电站围堰拆除在建的叶巴滩、拉哇和巴塘等水电站工程延期 堰塞湖应急处置易受到基础资料收集困难、周边环境危险、分析研究时间紧和溃决灾害扩大效应强等因素制约由此导致应急处置难度极大 因此考虑其致灾灾害链效应基于极其有限的数据信息建立堰塞湖溃决风险快速、定量评估方法是目前应急处置亟需解决的难题目前堰塞湖溃决风险评估研究主要分为两类:统计分析和定量分析 统计分析是基于堰塞湖统计信息评估溃决灾害损失划分不同类型灾害风险等级代表性的有水利部行业标准/堰塞湖风险等级划分与应急处置技术规范 该标准规定可根据堰塞湖坝体规模、物质组成和高度运用查表法和模糊数学法将其风险划分为极高、高、中和低 级 崔鹏等和 等基于汶川地震堰塞湖溃决实例根据堰塞湖坝体高度、库容和物质组成 项指标提出风险等级划分建议 等采用模糊数学方法对红石岩堰塞湖的风险等级进行了评估 定量分析模型多基于堰塞湖溃决洪水淹没情况与损失之间的关系构建评估模型 代表性的有石振明等、等构建的综合考虑人员疏散情况、生命损失和溃决参数等因素的定量分析模型以上两类评估方法均可实现溃决风险评估 第一类方法主要是基于经验的定性分析对堰塞湖风险快速评估具有较好的指导作用但其主要考虑堰塞湖坝体地貌学特征未考虑坝体的物质组成结构 而目前不同诱因导致堰塞湖坝体级配特征存在极大差异若未合理考虑其物质组成结构将无法准确评估堰塞湖风险 第二类方法虽然综合考虑了水力学参数和灾害损失但其风险评估过程中并未开展坝体稳定性快速评价无法快速实现稳定性的量化评估也并未考虑到溃决机理对溃决特征值变化的影响机制导致其溃决致灾的灾害链效应尚未明晰本文构建了一套堰塞湖溃决致灾的快速定量评估方法实现快速定量风险评估 本方法基于极其有限的数据信息构建考虑了堰塞湖和坝体相关物理力学指标、下游河道和影响区地形数据、人口分布情况实现了考虑颗粒级配组成的稳定性快速评价、基于冲蚀特性和崩塌过程的溃口流量变化过程分析、洪水演进过程模拟、生命损失评估预警明晰了溃决致灾的灾害链效应 该方法包括五部分:三维地理信息模块、稳定性快速评价模块、溃口流量快速定量分析模块、溃决洪水演进模块和生命损失快速评估模块将该方法应用于拥有完整溃决及演进过程实测资料的唐家山堰塞湖案例一方面验证了本文模型的可行性同时也为堰塞湖的风险评估和应急处置提供技术参考堰塞湖溃决快速风险评估方法堰塞湖溃决快速定量风险评估方法主要包括:信息收集稳定评价溃决洪水分析风险评估损失第 期王 琳等:基于溃决机理的堰塞湖溃决快速风险评估方法研究预警五部分 详细步骤如下:运用 插件快速获取研究区域的三维地形数据建立坝体和河道的三维地理信息模块实现信息快速收集基于逻辑回归算法建立考虑堰塞湖坝体形态特征、物质组成和上游堰塞湖水动力条件的稳定性快速评价模块实现堰塞湖坝体稳定性快速评价基于堰塞湖溃决机理开发基于物理机制并可在 内快速计算溃决洪峰流量、溃口崩塌过程、溃决时间和溃决库容等数据的溃决洪水分析模块实现溃决参数变化过程快速预测将三维地理信息数据和溃决洪水定量分析数据导入 软件构建溃决洪水演进计算模块开展风险评估在地理信息系统()中构建可视化溃决洪水淹没范围图实现风险快速评估预警揭示溃决洪水影响区域内的生命损失分布特征.三维地理信息模块基于中科院地理空间数据云(:/.)获取堰塞湖区域范围内的高程数据运用插件快速获取溃决及演进范围的三维地理信息 主要步骤如下:基于中科院地理空间数据云平台获取堰塞湖具体的经纬度信息及其等高线数据图高程数据精度越高数据量几何倍增加数据处理效率越低因此在实际应用时应根据具体工程实际选择所需数据精度 对于堰塞湖规模、下游河道信息和周边城镇地区等数据使用精度 的高程数据可满足应用需求 对于坝体高度、宽度和坝体体积等数据可在堰塞湖所在经纬度范围内使用精度为 的高程数据 运用 插件在 平台将河道数据、高程数据和其它三维地形信息数据直接导入并叠加实现图像数据和数字数据的转化 迅速获取堰塞湖的几何特征、下游影响区域等关键数据(内)堰塞湖实际应急处置时间非常紧迫例如唐家山堰塞湖在 个月内溃决 月 日形成的白格堰塞湖更是在 内发生溃决 应急处置面临着基础资料收集少、分析研究时间紧的困难在满足精度的要求下快速处理三维地理信息能够为后续稳定性评价溃决洪水分析等一系列工作提供数据支持.稳定性快速评价模块本模块采用的稳定性快速评价方法考虑了堰塞湖坝体级配特征、物质组成以及上游的水动力条件而现有评价方法多未考虑堰塞湖坝体物质组成 该方法可根据应急处置时获取的级配特征(如块石、砾石、粗细颗粒等)实现精细化和简化稳定性快速评价 具体表达式为().().().().()式中:()为堰塞湖坝体稳定性评价指标当()时堰塞湖坝体可判定为稳定/其中 为堰塞湖坝体高度 为堰塞湖坝体宽度即顺河向距离为堰塞湖坝体体积为堰塞湖流域面积为颗粒特征参数的取值方法在文献中具体阐述该稳定性快速评价方法的绝对准确率、保守准确率、判错率分别为:.、.、.同时与国内外同类方法比较(法、法、法和 法)结果表明该方法具有可靠性与优越性 其具体比较过程已在文献中详细说明文中不再进行赘述.溃口流量快速定量分析模块堰塞湖溃决风险评估方法中的溃决洪水定量分析模块多基于历史上的溃决资料未考虑溃决机理 虽然参数模型可利用较短计算时间完成溃决洪峰流量的模拟但无法获取溃决流量过程且未考虑土水耦合作用机制和冲蚀破坏过程 溃决洪水快速定量分析是风险评估和应急处置方案制定的基础 因此构建基于溃决机理考虑坝体材料冲蚀特性和溃口崩塌过程的快速、定量溃决过程计算模型是提高堰塞湖溃决洪水风险评估准确性的必然需求本文针对 世纪以来典型的唐家山、易贡、小岗剑、红石岩和白格等堰塞湖的实测溃决数据进行梳理并在此基础上提出了基于物理机制的溃决洪水分析模型(:/./)免费下载使用如图 所示 为水库水位()为时间()为库容()为溃口宽度()为跌落系数 为溃口底高程()为水深()为流量系数(/)、为双曲线侵蚀模型系数、分别为流速和初始流速(/)为容重(/)为黏聚力()为摩擦角()、为双曲线倾角扩展系数、分别为溃口倾角和初始倾角()此模型采用宽顶堰公式计算溃决流量基于圆筒型冲蚀试验设备()的试验数据提出双曲线形式的冲蚀本构模型在计算溃口崩塌时基于总应力法采用双曲线图表形式快速确定和模拟溃口崩塌过程 当溃口不断被冲深时发生侧向破坏可通自 然 灾 害 学 报第 卷过溃口边坡稳定分析逐级确定临界坡面和滑裂面 溃决分析程序可保证工程师在 内完成目标工况计算在堰塞湖应急处置时可快速应用图 堰塞湖溃决洪水分析模型().().洪水演进模块将 插件中构建的河道断面三维模型导入 软件构建河道模型将其分为左漫滩、中心区域和右漫滩 个区域 确定不同区域的曼宁系数以及演进计算的上下游边界确定溃决及演进区域的中心线、河岸线、行洪区线、区域面积和糙率等数据 溃决演进计算采用开源的 .软件的非恒定流模型计算实现溃决洪水演进的快速准确模拟 基于非恒定流的方程采用差分法进行求解通过划分河道断面、水流路径等模拟河道信息水位流量关系为边界条件时间步长设置为 基于四点隐式差分法进行迭代计算非恒定流模拟是基于水流的连续方程与动量方程展开的其连续方程为()()式中:为液体密度 为流速动量方程为 ()式中:为质量力 为压力 为运动黏滞系数.基于学习机网络模型的生命损失评估模块堰塞湖溃决产生的洪水会对下游地区造成严重损失其中以生命损失最为重要 溃决洪水造成生命损失的致灾环境复杂亟需在堰塞湖应急处置时开展预警计算溃决可能造成的损失及时撤离群众目前的生命损失评估方法多存在参数量化过程不统一影响因素定义模糊评估精度不高等问题 本文方法基于极限学习机网络模型建立风险人口与生命损失间的函数关系考虑了主要影响因素与次要影响因素间第 期王 琳等:基于溃决机理的堰塞湖溃决快速风险评估方法研究的作用关系与影响程度并通过对生命损失影响因素的统一量化结合前文计算的溃决洪水定量分析数据与地理信息数据获取溃决洪水深度、淹没面积和水位变化过程绘制淹没范围图确定溃决洪水到达城镇时间以及淹没水深和淹没范围在堰塞湖灾害所造成的生命损失中风险人口和人口密度等人口统计学变量常作为关键参数参与生命损失模型的计算如王志军等的模型 在实际运用中由于城镇的人口分布通常较为集中在城镇中心区域的人口密度通常较大其他区域人口密度较小对生命损失计算影响较大 在本模型中通过溃口流量与洪水演进模块所计算出的下游淹没水深能够快速预估出下游淹没范围并将该范围作为生命损失模块的计算起始条件 在这个范围内根据人口密度将受灾区域进一步划分为高密度区域、低密度区域不同密度区域中的修正系数、取值不同 在计算过程中根据式()对各区域分别计算求和最终获得淹没范围内的生命损失 为 ()()式中:为生命损失 为主要因素修正系数 为次要因素修正系数 为人口密度参数 取.案例分析.唐家山堰塞湖溃决过程概述 年 月 日受汶川.级地震影响四川通口河右岸发生特大滑坡堵塞河流形成唐家山堰塞湖 其距下游北川县城 绵阳市 下游有将近 万人民还有宝成铁路、输油管道等重要基础设施如图 所示 为减少唐家山堰塞湖溃决带来的巨大灾害应急处置开挖了宽为 、深为 、总长为 的引流槽开挖后的堰塞湖槽底高程为.年 月 日槽内水深和宽度开始逐步增加此时上游来流量为 /堰塞湖水位逐渐抬升 堰塞湖水位至 月 日 达到最大值.从 到 测量的平均流速值在./之间相应的水流剪应力在 之间 堰塞湖溃决时颗粒的启动流速大约为./至 年 月 日溃口达到峰值流量约为 /溃决过程中槽底的平均侵蚀率约为./对应的流速为./剪应力为 月 日溃口流速(约./)及流量基本保持不变表明堰塞湖溃决基本结束此时堰塞湖水位为.整个过程中唐家山堰塞湖的下泄水量约.溃口最终深度为 顶部宽度约为 底部宽度约为 溃口洪水最大流速为./月 日 洪水到达绵阳市未出现较大洪峰排险成功结束 唐家山堰塞湖拥有珍贵的溃决全过程实测资料为堰塞湖溃决研究提供了翔实的基础数据堰塞湖相关参数如表 所示 下游北川(距唐家山坝址.)、通口(距唐家山坝址.)、涪江桥(距唐家山坝址.)等水文站的实测洪水特征值见表 根据前述堰塞湖溃决风险快速定量评估方法将其应用于唐家山堰塞湖风险评估具体模拟方法如下图 唐家山堰塞湖及附近水文站.表 唐家山堰塞湖相关参数 名称参数数值堰塞湖坝体体积/().最高点高程/最低点高程/.顺河向长度(堰塞湖底部)/垂直河向长度/.左侧顺河向长度与高度比.流域面积/堰塞湖最大可能库容/().溃决库容/().原始河床高程/入库流量/(/)溃决洪峰流量/(/)启动流速/(/).自 然 灾 害 学 报第 卷表 唐家山堰塞湖及下游水文站的水文特征值 站名洪峰流量(/)最高水位/洪峰出现时刻唐家山.月 日 北川.月 日 通口.月 日 涪江桥.月 日.三维地理信息模型建立唐家山堰塞湖位于通口河沿通口河向下游途经北川县城、通口镇、含增镇、青莲镇、龙凤镇、石马镇以及绵阳市等城市和乡镇约有 万居民 因此选取的溃决计算及演进计算区域为唐家山堰塞湖所在地沿通口河至绵阳市的矩形区域(、)基于 插件共确定 个断面河道中心线长度为.唐家山堰塞湖底高程约为 计算区域高程急剧下降约 如图 所示另外图 标注了沿河段的 个典型水文站位置 将基于 插件构建的河道模型导入 软件构建溃决洪水演进模型 在唐家山至绵阳市河道选择 个断面作为典型断面其中 号、号、号、号断面分别为唐家山堰塞湖、北川水文站、通口水文站、涪江桥水文站和绵阳市如图 所示 图 通口河段剖面图 图 溃决洪水演进模型构建 .稳定性快速评价唐家山堰塞湖坝体高度为 堰塞湖坝体宽度为 、堰塞湖坝体体积为.堰塞湖流域面积为 其颗粒组成以碎石为主粒径多小于 粉质黏土达到 碎石约为 块石约为 颗粒特征参数 取为.将各参数代入式()获得().堰塞湖坝体稳定性判断结果为“不稳定”亟需开展溃决洪水模拟分析.溃决洪水过程定量分析采用 开展溃决洪水过程定量分析输入参数如表 所示 图 和表 为溃口流量、堰塞湖水位、溃口水流流速、溃口宽度、溃口底高程的实测值和计算值的对比情况 反演分析结果显示 月 日溃口出现峰值流量为./溃决库容达到.溃口顶部宽度为.溃决水流最大流速为./若不开挖宽为 、深为 的引流槽则在坝顶 时堰塞湖将溃决库容将达到.峰值将达到./溃口顶部宽度为.溃决最大流速为./开挖引流槽可使溃口峰值流量减少.洪量减少.反演分析的误差均在 以内因此可以认为反演结果能够验证模型的适用性能够在保证精度需求的前提下进行快速反演 因此运用基于物理机制的溃决洪水分析程序可快速、定量地评价堰塞湖溃决过程第 期王 琳等:基于溃决机理的堰塞湖溃决快速风险评估方法研究表 溃决过程反演分析 输入参数 参数类别参数名称符号输入值堰塞湖特征库容、.、.、初始库水位/.入库流量/(/)初始底高程/溃口初始宽度.水力学参数宽顶堰系数.跌落系数.启动流速/(/).冲蚀参数、.、.岩土力学参数容重/(/).凝聚力/摩擦角/()初始溃口倾角/()注:库容关系曲线由相应 组库水位高程()和库容()、内插决定为相对高程起点图 唐家山堰塞湖溃决计算结果与实测数据对比.表 唐家山堰塞湖溃决参数计算结果对比 溃决参数实测值有引流槽(反演分析)无引流槽误差/溃口峰值流量/(/).峰值流量出现时间.最大流速/(/).最大流速出现时间.溃口最终顶宽/.溃口最终底高程/.溃口下泄库容/().溃决洪水演进分析溃决洪水演进分析输入参数如表 所示 将计算结果与具有实测资料的下游 个水文站资料进行分析自 然 灾 害 学 报第 卷对比验证洪水演进模块的可行性 唐家山堰塞湖下游 个水文站反演计算与实测数据对比结果如图 及表 所示表 洪水演进输入参数表 左岸糙率 右岸糙率 河道糙率 缩放系数扩张系数流量/(/)上游关系曲线下游坡度.非恒定流 曲线图 唐家山堰塞湖下游 个水文站洪峰流量过程计算与实测数据对比.表 唐家山堰塞湖下游 个水文站实测与计算结果对比 水文站名称洪峰流量 实测/(/)反演分析/(/)误差/洪峰到达时间 实测时间反演分析时间误差/最大淹没水深 实测/反演分析/误差/北川水文站.通口水文站.涪江桥水文站.针对北川水文站实测北川水文站溃决洪峰时间为 月 日 溃决流量为 /最高水位为.反演分析的溃决洪峰时间为 溃决流量为 /最高水位为.溃决洪峰的误差为.最高水位的误差为.到达溃决洪峰时间误差为.反演分析与实测数据间相关系数为.针对通口水文站实测通口水文站溃决洪峰时间为 月 日 溃决流量为./最高水位第 期王 琳等:基于溃决机理的堰塞湖溃决快速风险评估方法研究为.反演分析的溃决洪峰时间为溃决流量为./最高水位为.溃决洪峰的误差为.最高水深的误差为.到达溃决洪峰时间误差为.反演分析与预测数据间相关系数为.针对涪江桥水文站实测涪江桥水文站溃决洪峰时间为 月 日 溃决流量为./最高水位为.反演分析的溃决洪峰时间为溃决流量为./最高水位为.溃决洪峰的误差为.最高水深的误差为.到达溃决洪峰时间误差为.反演分析与实测数据间相关系数为.唐家山堰塞湖溃决后到达下游北川、通口、涪江桥水文站的溃决流量、洪峰到达时间、最大水深等溃决参数的计算误差均在 以内且反演数据与实测数据的相关系数均到达.以上在洪峰局部与过流全过程均验证了洪水演进模块的可行性 计算误差的来源主要为地形信息误差、河道糙率数据、断面划分密度、理论计算依据等的影响这是洪水演进研究中不可避免的且随着模拟时间与洪水演进这些缺陷将被不可避免地逐渐放大.生命损失风险评估图 唐家山堰塞湖溃决后最大淹没范围.基于溃决洪水演进分析数据在 中可视化并构建溃决洪水淹没范围图获取溃决后影响区域的淹没范围、淹没水深和淹没历时评估堰塞湖溃决洪水影响区域内的生命损失分布特征 唐家山堰塞湖溃决后最大淹没范围如图 所示为规避面临的溃决风险绵阳市大约有.万人在大坝溃决前 被安全撤离并没有人员丧生 本文主要研究若没有进行及时疏散唐家山堰塞湖溃决后将面临的生命损失采用生命损失模块预测实际溃决过程和相应的预警时间 基于洪水演进计算结果唐家山堰塞湖下游 个城镇和绵阳市均有可能遭受溃决洪水的威胁如表 所示不同预警条件下的生命损失计算结果如表 所示 计算结果表明生命损失与预警时间存在明显的相关性若预警时间仅为.时生命损失率达到.当预警时间大于.时基本上可做到无生命损失 而唐家山堰塞湖提前 预警保证了无人员丧生表 唐家山堰塞湖下游信息 下游城镇距离/河岸高度/人口数量/人北川镇.通口镇.含增镇.青莲镇.龙凤镇.石马镇.绵阳市.表 不同警报时间下生命损失评估 预警时间/生命损失/人北川镇通口镇含增镇青莲镇龙凤镇石马镇绵阳市.结论)本文提出了在堰塞湖突发情况下的快速风险评估方法该方法分为:信息收集稳定评价溃决洪水分析风险评估损失预警共五部分可实现坝体稳定性评价、溃决洪水分析、演进过程快速计算和生命损失分布特征预警的溃决快速定量风险评估 方法考虑了堰塞湖坝体的级配特征实现精细化和简化稳定性快速评价考虑溃决机理建立了基于冲蚀特性和崩塌过程的溃口流量变化过程分析物理模型实现溃决自 然 灾 害 学 报第 卷流量过程快速分析(内)基于极限学习机网络模型考虑风险人口与生命损失间的函数关系实现了生命损失分布特征预警明晰了溃决致灾的灾害链效应)堰塞湖的溃决风险与坝体形态特征、物质组成和堰塞湖水动力条件密切相关基于有限信息合理评价堰塞湖坝体的稳定性是溃决风险评估的先决条件 构建基于溃决机理考虑坝体材料冲蚀特性和溃口崩塌过程的快速、定量溃决过程计算模型获取水位、溃口流量、溃口宽度、流速和溃口底高程等溃决数据的变化过程是堰塞湖风险评估的关键)基于 技术可视化构建溃决洪水淹没范围图迅速获取溃决后影响区域的淹没范围、淹没水深和淹没历时是准确开展风险评估的重要保证 根据溃决洪水的淹没情况充分考虑洪水影响区域内的人口分布特征开展生命损失评估是风险评估和抢险救灾的重要目标)在唐家山堰塞湖应急处置中若不开挖引流槽则在坝顶 时堰塞湖将溃决库容将达到.峰值将达到./溃口宽度为.溃决最大流速为./开挖引流槽可使洪峰流量减少.洪量减少.因此对于库容较大的堰塞湖在条件允许的情况下可通过开挖引流槽降低溃决风险)本文的评估方法基于堰塞湖和坝体的物理力学特性有效结合三维地理信息考虑堰塞湖坝体的稳定性、冲蚀特性和溃口崩塌以及复杂条件下的溃决洪水演进过程可快速、定量评估堰塞湖溃决洪水风险为堰塞湖的应急处置提供技术支撑仍需考虑经济损失和环境影响等类别的风险后果需进一步研究考虑致灾后果的堰塞湖溃决风险评估方法参考文献:李君 张新 陈燕 等.萨雷兹湖东向引水最低成本路径及可行性分析.自然灾害学报 ():.():.()蔡耀军 杨兴国 张利民 等.堰塞湖风险评估快速检测与应急抢险技术和装备研发研究构想与成果展望.工程科学与技术():.():.()刘威 何思明.金沙江沃达潜在滑坡诱发灾害链成灾过程数值模拟.工程科学与技术 ():.():.().:.():.王琳 武钰淼 段庆伟 等.堰塞湖泄流槽开挖措施快速优化研究.水力发电学报 ():.():.().():.:.:.():.():.:.():.():.():.第 期王 琳等:基于溃决机理的堰塞湖溃决快速风险评估方法研究 许强 郑光 李为乐 等.年 月和 月金沙江白格两次滑坡堰塞堵江事件分析研究.工程地质学报 ():.():.()邓建辉 高云建 余志球 等.堰塞金沙江上游的白格滑坡形成机制与过程分析.工程科学与技术 ():.():.()赵洋 李守义.堰塞体逐渐溃决数学模型研究.自然灾害学报 ():.():.().:.():.():.刘宁 张建新 林伟 等.汶川地震唐家山堰塞引流除险工程及溃坝洪水演进过程.中国科学(辑:技术科学)():.(:)():.().:.():.徐凡献.一种滑坡堰塞坝稳定性的快速评估模型.自然灾害学报 ():.():.()/堰塞湖风险等级划分与应急处置技术规范.北京:中国水利水电出版社./.:.()崔鹏 韩用顺 陈晓清.汶川地震堰塞湖分布规律与风险评估.四川大学学报(工程科学版)():.()():.().():.():.石振明 熊永峰 彭铭 等.堰塞湖溃坝快速定量风险评估方法 以 年鲁甸地震形成的红石岩堰塞湖为例.水利学报():.:.():.().():.杨志强 裴向军 裴钻 等.基于有效应力原理的溃坝分析.自然灾害学报 ():.():.()赵安 王婷君.基于过程机理的洪灾生命损失评价模型框架初探.自然灾害学报 ():.():.()单熠博 陈生水 钟启明.堰塞体稳定性快速评价方法研究.岩石力学与工程学报 ():.():.().():.():.(/):.自 然 灾 害 学 报第 卷 .:.:.():.():.:.:():.介玉新 谢琪 周婷.边坡滑动面构造的基函数.水力发电学报 ():.():.().:():.:.:.:.():.:.王志军 顾冲时 刘红彩.基于模糊物元与指数平滑法的溃坝生命损失估算.长江科学院院报 ():.():.().:().():.:():.():.:.():.