本发明涉及泥石流危险性评价方法领域，特别是一种泥石流实时动态危险性评价方法。
背景技术：
：地震引发的大量崩塌、滑坡等山地灾害可以为泥石流的形成提供丰富的物源，大大降低泥石流的致灾条件，导致其规模、发生频率显著增加。如四川省绵竹市清平乡在汶川地震后的2010年8月13日局地大暴雨的激发下，11条沟同时暴发泥石流，形成了长达3.5km、宽400～500m、平均厚约5m、总方量约600×104m3的堆积区，覆盖面积约达120×104m2，造成379户房屋受损，6000余人遭受不同程度灾害，直接经济损失达6亿元左右。震后泥石流严重阻碍了山区的经济发展和生态建设，对它们进行精准的危险性评价能为区域宏观规划和防灾减灾提供重要决策依据。目前国内外学者提出的泥石流危险性评价方法大致可以分为统计分析方法、人工智能方法、数值模拟及其他方法四类。例如，linetal(2012)基于模糊数学理论建立了泥石流风险评估模型，用以评估台湾地区泥石流发生的可能性。niuetal(2014)运用逐步判别分析等方法评价中国西南山区泥石流的危险性。changetal(2010年)运用神经网络(nn)和遗传算法(ga)预测了台湾东部地区泥石流的危险性。castellietal(2017)运用flo-2d数值模型模拟真实的泥石流过程。尽管这些泥石流危险性模型变得越来越完善，但是大多数模型，特别是前两类模型，主要以静态的危险性评价结果为主。这两类模型通常都没有考虑泥石流危险性的动态变化。然而，由于泥石流的诱发因素(如降雨)和基本因素(如地貌、植被覆盖、土地利用)是随时间和空间而不断变化的，因此泥石流的危险性并不是一个静态的概念，而是随时间不断发生着变化的。数值模拟方法虽然可以在泥石流危险性评价中考虑时间因素，但模型所需的参数一般难以准确获取，而且模拟过程较为复杂，不便于推广应用。flow-r是瑞士洛桑大学hortonetal(2011，2013)提出的一种基于gis的区域重力灾害危险性评价数值模型。不同于一般的因子叠加方法，该模型集成了山地灾害不同演化阶段的多个物理力学模型来对山地灾害的危险范围进行定量分析。它可以在区域尺度范围内以较少的数据量快速评估重力灾害的危险性。该模型自提出后，引起广泛关注，已被应用于瑞士、意大利、美国、韩国等国家滑坡、崩塌、泥石流等地质灾害的危险性评价中。例如，rahmanetal(2017)在研究中应用flow-r模型模拟滑坡源区的运动以评估滑坡危险范围；gladeetal(2012)建立了多灾种灾害风险评估框架，应用flow-r模型对泥石流、崩塌、雪崩等灾害的危险性进行评估；blais-stevens和behnia(2016)运用flow-r模型对公路走廊沿线泥石流的敏感性进行了分析识别。与数值模拟方法相比，flow-r模型不仅可利用地形地貌、地层岩性等较易获取的评价指标进行泥石流源区分析，并且运用能量守恒定律和传播概率原理对泥石流源区的运动范围进行分析，而且在泥石流资料不完整的区域，参数易获取、可操作性强，在应用方面具有很大的优势。但是，flow-r模型没有考虑灾害的时间效应，仍然是静态评估模型。它只能考虑受环境背景条件影响的泥石流源区分布，并不能考虑由降雨诱发滑坡为泥石流提供的动态源区分布。降雨诱发滑坡为泥石流提供的动态源区可以通过降雨诱发滑坡的确定性模型来计算。这些模型考虑了降雨入渗过程中地下水变化对边坡稳定性的影响，是目前降雨诱发滑坡稳定性分析的主要研究方向。目前，研究人员提出了各种确定性模型来分析降雨引起的浅层滑坡，例如shalstab模型(dietrichetal.1998)，sinmap模型(packetal.1998)，iverson模型(iverson，2000)，李宁模型(liningetal.2012)和trigrs模型(baumetal.2010)。在这些模型中，trigrs模型近年来得到了广泛的应用和研究。该模型不仅具有良好的理论框架，还可以考虑岩土体参数的空间变异性、降雨条件的动态变化以及土壤的初始条件。滑坡和泥石流是两种既有明显区别又有紧密联系的不同的地质灾害类型。因为滑坡可以为泥石流提供松散的物质来源。因此，研究者们也常常把滑坡和泥石流结合起来进行研究。例如，gomesetal(2013)耦合了shalstab和flo-2d模型，stancanellietal(2017)耦合了trigrs和flo-2d模型，hsuandliu(2019)耦合了trigrs和debris-2d模型。尽管这些模型的耦合研究仍存在一定的局限性和不足，但是这些耦合模型为准确评估泥石流的危险性提供了新的途径。前已述及，泥石流是一个开放的动态系统。在实践中，泥石流灾害的诱发因素、物源和危险性等都会随时间的变化而变化。因此，仅通过运用静态评估方法或片面考虑泥石流的固体物源，很难为泥石流灾害的防灾减灾提供精准的决策依据。特别是“5·12”汶川地震后，地震灾区斜坡体的地貌、地质结构和植被覆盖以及沟谷内的松散固体物质都发生了很大变化。同一条泥石流沟在地震前后的活动性有着显著差异。例如，汶川地震灾区八一沟地震发生前植被发育良好，极少发生地质灾害。但是，在地震之后，从2008年到2011年的4年中，该沟谷就发生了5次大规模的泥石流。总冲出量就达150万m3，淤埋公路长约200m，几乎摧毁了沟内全部泥石流治理工程及安置板房。八一沟位于“5·12”汶川地震极重灾区，四川省都江堰市龙池镇云华村，沟口坐标为北纬31°3′11″n，东经103°32′25″e。八一沟在地势上整体呈现北、东、西部高，中部及南部低，横断面呈“v”字型，为易汇水及易堆积松散物地形，平均河床纵坡降和大部分坡体坡度均属泥石流易发范围。八一沟流域自沟顶到沟底依次分布震旦系安山岩及安山玄武岩、澄江-晋宁期花岗岩及三叠系砂岩、泥岩，第四系以碎块石土为主分布于堆积区和流通区。“5·12”汶川地震的主发震断裂映秀-北川发震断裂带的北支和南支同时横穿八一沟。流域处于亚热带湿润季风气候区，四季分明、夏季暴雨频繁。根据《汶川8.0级地震烈度分布图》，区域地震烈度为ⅹ度。受地震的强烈影响，区内产生大量崩滑体堆积在沟谷及两岸斜坡地带。八一沟具备了泥石流发育的有利地形、丰富的物源补给及充足的降雨条件，在一定的降雨条件下，很有可能再次形成泥石流，对周边的村庄、水库、公路等造成严重的威胁。技术实现要素：为解决现有技术中存在的问题，本发明提供了一种泥石流动态危险性评价方法，解决了以往方法没有全面考虑泥石流物源分布和多侧重于泥石流危险性的静态评价研究的问题。本发明采用的技术方案是，一种泥石流实时动态危险性评价方法，包括以下步骤：s1：将实时降雨数据、岩土体物理力学性质参数和水文参数输入到trigrs模型，输出在不同降雨历时下表征斜坡稳定性状态的稳定性系数fs；s2：将trigrs模型输出的斜坡稳定性状态分为fs＜1、fs≥1这2类；s3：将fs＜1对应为泥石流源区分类标准的“是源区”，将fs≥1对应为泥石流源区分类标准的“非源区”；s4：根据泥石流的环境背景和成因条件，确定汇流累积量、坡度、平面曲率、岩性等背景因子作为泥石流危险性的评价因子，并结合flow-r模型中各因子的源区识别阈值和泥石流的实际特征，将各评价因子图层，按照分类判据划分为“是源区”、“非源区”或者“不确定”三种类型；s5：将不同降雨历时下斜坡的稳定性系数图层fs和影响泥石流物源的各背景因子图层输入到flow-r模型中进行叠加分析，并根据“至少一次被分类为是源区，而未被分类为非源区”的栅格作为泥石流源区的识别规则，识别出不同降雨历时条件下泥石流的坡面潜在物源区的分布，实现trigrs模型和flow-r模型的耦合；s6：基于野外实地调查和遥感解译，再识别出泥石流沟道中已有的松散物源区的分布；s7：将不同降雨历时条件下泥石流的坡面潜在物源区和沟道中已有的松散物源区综合，即为不同降雨条件下泥石流总物源区的分布；s8：在此基础上，基于flow-r模型中的能量守恒方法和流向概率算法计算不同降雨历时条件下泥石流的能量和危险性概率分布；s9：根据不同降雨条件下泥石流的能量和危险性概率分布，进行泥石流的实时动态危险性分区。优选地，s5包括以下子步骤：s01：斜坡的动态稳定性分析，斜坡的动态稳定性分析利用trigrs模型进行分析：s02：泥石流的源区识别，泥石流的物源区识别包括坡面潜在物源识别和沟道已有松散物源识别，前者用flow-r模型进行识别，后者基于野外实地调查和遥感解译识别；s03：泥石流的动态危险性分析，泥石流的动态危险性的计算用flow-模型；s04：泥石流的动态危险性分区。优选地，s01包括以下子步骤：s011：确定水文参数；s012：确定物理力学性质参数；s013：确定土层厚度；s014：确定实时降雨条件；s015：斜坡动态稳定性分析。优选地，s02包括以下子步骤：s021：识别坡面潜在动态物源；s022：识别沟道已有松散物源。优选地，s03包括以下子步骤：s031：计算泥石流的动态危险性概率分布；s032：计算泥石流的动态能量分布。优选地，s04具体为将不同降雨时长的泥石流危险性概率值按分级标准分为为极高危险、高度危险、中度危险和低度危险4个等级，其等级标准为：低危险区为0.0003～0.004，中度危险区为0.004～0.032，高度危险区为0.032～0.11，极高危险区为0.11～1。本发明泥石流动态危险性评价方法的有益效果如下：1.耦合模型不仅考虑了沟道内已有的松散物源以及泥石流成灾环境背景条件对泥石流坡面潜在物源的影响，而且还考虑了降雨诱发滑坡为泥石流提供物源的动态变化。因此，泥石流源区识别考虑得更加全面、更加符合实际。2.耦合模型是在实时降雨条件下对泥石流进行的危险性评价，降雨条件是动态变化的，泥石流危险性评价结果也是动态变化的。这不同于前人多侧重于泥石流危险性的静态评价研究。附图说明图1为本发明泥石流实时动态危险性评价方法的总流程框图。图2为本发明泥石流实时动态危险性评价方法的φ为定值、c取0-30时3个评估指标变化曲线图。图3为本发明泥石流实时动态危险性评价方法的总正确率和查全率同时都大于50％的岩土体参数组合图。图4为本发明泥石流实时动态危险性评价方法的2010年8月13日龙池镇主要降雨过程图。图5为本发明泥石流实时动态危险性评价方法的trigrs模型不同降雨历时斜坡稳定性状态图。图6为本发明泥石流实时动态危险性评价方法的斜坡稳定性状态、dem、汇流累计量、坡度、平面曲率以及岩性分布图。图7为本发明泥石流实时动态危险性评价方法的泥石流总源区动态分布图。图8为本发明泥石流实时动态危险性评价方法的泥石流危险性概率动态分布图。图9为本发明泥石流实时动态危险性评价方法的泥石流能量动态分布图。图10为本发明泥石流实时动态危险性评价方法的不同危险等级面积的动态变化图。图11为本发明泥石流实时动态危险性评价方法的泥石流危险性分区图。图12为本发明泥石流实时动态危险性评价方法与flow-r模型的源区、危险性概率和能量分布结果对比图。图13为本发明泥石流实时动态危险性评价方法与flow-r模型的评价效果对比图。具体实施方式下面结合附图，以汶川地震灾区的八一沟泥石流为例，对本发明的具体实施方式做进一步详细的描述，以便于本
技术领域：
的技术人员理解本发明。但应该清楚，本发明不限于具体实施方式的范围。对本
技术领域：
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。1、斜坡的动态稳定性分析这部分需要先确定岩土体的水文参数、物理力学性质参数和土层厚度等参数，各个参数的确定过程如下：(1)水文参数根据研究区土颗粒级配组成情况，用经验统计公式(1)计算八一沟碎石土的饱和渗透系数ks，取为1.30e-06m/s。ks＝-18.4 0.1836d1 0.2189d2 0.1401d3 0.1834d4(1)式中，ks为碎石土的渗透系数(10-4cm/s)，d1为大于20mm粒径的粒组含量(％)，d2为2-20mm粒径的粒组含量(％)，d3为0.1-2mm粒径的粒组含量(％)，d4为小于0.1mm粒径的粒组含量(％)，参考郭朝旭等人(2015)研究中的颗粒级配图，d1取14.83、d2取62.67、d3取20.15、d4取2.35。另据文献(parketal.2016)，trigrs模型中的水力扩散度d0和稳定渗透率iz分别按照渗透系数ks的200倍和0.01倍确定。(2)物理力学性质参数在山地灾害的危险性评价中，混淆矩阵常被用于评估模型的优劣和参数选取的好坏。混淆矩阵的评估指标及评估统计量如表1和表2。表1混淆矩阵的评估指标实况有实况无预测有真阳(tp)假阳(fp)预测无假阴(fn)真阴(tn)表2混淆矩阵统计量运用参数反演法确定岩土体的物理力学性质参数，就是根据混淆矩阵的3个评估统计量，从大量岩土参数组合中筛选出结果最优的一组作为岩土体的参数。首先，确定c、φ值的反演范围，c取0～30kpa、φ取10°～31°；然后，以1为间隔取不同的c、φ值进行两两组合。当φ为定值时、c为唯一变量在0～30kpa范围内取值分析。发现随着c值的变大，总正确率逐渐变大，而滑坡正确率和面积比逐渐变小(图2)。此时，满足总正确率和查全率同时都大于50％的岩土参数组合将大大减小，并且在矩阵中呈带状分布(图3)。最后综合面积比的取值情况，选择c＝21kpa、φ＝28°作为岩土体参数。(3)土层厚度八一沟流域土层厚度按表3设置。表3坡度与土厚关系表坡度/°土壤平均厚度/m<305.030-404.040-603.0>602.0(3)实时降雨条件八一沟所处的龙溪河流域降水主要集中于5-9月，占全年总降水量的80％，平均月降水量最大月份为8月。2010年8月13日在暴雨的激发下八一沟发生了重大泥石流，选此次实际的降雨过程作为八一沟的降雨条件(图4)。由于前期有一次降雨过程，故取饱和情况对流域内滑坡的稳定性进行分析。(4)斜坡的动态稳定性分析将所确定的岩土体参数、实时降雨强度、降雨时间等输入trigrs模型，即公式(2)，即可得出不同降雨时刻斜坡的稳定性状态(图5)。式(2)中，φ为土的内摩擦角；c为土的凝聚力；γs为是土容重；γw是地下水容重；fs为稳定性系数；α为坡度；ψ(z,t)为压力水头，是深度和时间的函数。其中，是与重力作用相关的项，是与内聚力相关的项，是与水压力相关的项。2、泥石流的源区识别泥石流总源区包括坡面潜在物源和沟道已有松散物源。(1)识别坡面潜在动态物源首先，将trigrs模型输出的斜坡稳定性状态分为fs＜1、fs≥1这2类；将fs＜1对应为泥石流源区分类标准的“是源区”，将fs≥1对应为泥石流源区分类标准的“非源区”(表5)；然后，根据泥石流的环境背景和成因条件，确定汇流累积量、坡度、平面曲率、岩性等背景因子作为泥石流危险性的评价因子(图6)，并结合flow-r模型中各因子的源区识别阈值和泥石流的实际特征，将各评价因子图层按照分类判据划分为“是源区”、“非源区”或者“不确定”三种类型(表4)；最后，将不同降雨历时下斜坡的稳定性系数图层fs和影响泥石流物源的各背景因子图层输入到flow-r模型中进行叠加分析，并根据“至少一次被分类为是源区，而未被分类为非源区”的栅格作为泥石流源区的识别规则，识别出不同降雨历时条件下泥石流的坡面潜在物源区的分布，并且实现trigrs模型和flow-r模型的有机耦合。表4耦合模型各评价因子源区识别阈值及分类判据(2)识别沟道已有松散物源泥石流沟道内已有的松散物源基于野外实地调查和遥感解译综合确定。综合坡面动态潜在物源和沟道内已有松散物源，即可得到不同降雨历时下泥石流的总动态源区(图7)。3、泥石流的动态危险性分析泥石流的动态危险性分析包括计算泥石流的动态危险性概率分布和动态能量分布，主要基于flow-r模型进行分析。选用修正的holmgren模型计算八一沟泥石流流向概率，参考melo等、horton等人研究取holmgren方法修正系数dh为2m、坡度指数x为4；flow-r模型中惯性权重分配方法采用cosines方法。选用流向概率计算方法及运动参数见表5：表5流向概率计算方法及运动参数表在flow-r模型中设置流向概率计算方法并输入运动参数，即可计算得到泥石流的危险概率分布(图8)和能量分布(图9)。4、泥石流的动态危险性分区将不同降雨时长的危险概率值按分级标准分为为极高危险、高度危险、中度危险和低度危险4个等级。具体的分级标准为：0.0003～0.004为低危险区，0.004～0.032为中度危险区，0.032～0.11为高度危险区，0.11～1极高危险区。不同降雨历时下不同危险等级面积变化情况见图10。降雨7h后的泥石流危险性分区结果见图11。5、评价效果为了验证耦合模型的评价效果，将该耦合模型的评价结果与flow-r模型的评价结果进行对比分析(图12)。两种模型的危险范围分别与基于遥感影像提取出的八一沟历时泥石流范围进行叠加分析(图13)，分别计算混淆矩阵正确率、查全率、模拟与实际泥石流面积比三个评估指标，结果见表6。表6混淆矩阵评估结果从表6中可以看出，耦合模型的总正确率、查全率和模拟与实际泥石流面积比均明显优于flow-r模型。本发明，基于trigrs和flow-r模型的有机耦合，提出了一个泥石流实时动态危险性评价的新方法。该方法的创新性在于，它不仅考虑了降雨条件和物源条件对泥石流危险性的动态影响，而且综合考虑了成灾背景条件、降雨诱发滑坡以及沟道内已有松散物质对泥石流物源的综合贡献作用。它能更科学地定量评估动态降雨条件下泥石流的危险性，对提高泥石流评价的准确性以及更精准地指导防灾减灾具有很好的理论和实际意义。当前第1页1 2 3
技术特征：1.一种泥石流实时动态危险性评价方法，其特征在于，包括以下步骤：
s1：将实时降雨数据、岩土体物理力学性质参数和水文参数输入到trigrs模型，输出在不同降雨历时下表征斜坡稳定性状态的稳定性系数fs；
s2：将trigrs模型输出的斜坡稳定性状态分为fs＜1、fs≥1这2类；
s3：将fs＜1对应为泥石流源区分类标准的“是源区”，将fs≥1对应为泥石流源区分类标准的“非源区”；
s4：根据泥石流的环境背景和成因条件，确定汇流累积量、坡度、平面曲率、岩性等背景因子作为泥石流危险性的评价因子，并结合flow-r模型中各因子的源区识别阈值和泥石流的实际特征，将各评价因子图层按照分类判据划分为“是源区”、“非源区”或者“不确定”三种类型；
s5：将不同降雨历时下斜坡的稳定性系数图层fs和影响泥石流物源的各背景因子图层输入到flow-r模型中进行叠加分析，并根据“至少一次被分类为是源区，而未被分类为非源区”的栅格作为泥石流源区的识别规则，识别出不同降雨历时条件下泥石流的坡面潜在物源区的分布，实现trigrs模型和flow-r模型的耦合；
s6：基于野外实地调查和遥感解译，再识别出泥石流沟道中已有的松散物源区的分布；
s7：将不同降雨历时条件下泥石流的坡面潜在物源区和沟道中已有的松散物源区叠加，即可得到不同降雨条件下泥石流总物源区的分布；
s8：在此基础上，基于flow-r模型中的能量守恒方法和流向概率算法计算不同降雨历时条件下泥石流的能量和危险性概率分布；
s9：根据不同降雨条件下泥石流的能量和危险性概率分布，进行泥石流的实时动态危险性分区。
2.根据权利要求1所述的泥石流实时动态危险性评价方法，其特征在于，所述s5包括以下4个子步骤：
s01：斜坡的动态稳定性分析，斜坡的动态稳定性分析利用trigrs模型进行分析；
s02：泥石流的源区识别，泥石流的物源区识别包括坡面潜在物源识别和沟道已有松散物源识别，前者用flow-r模型进行识别，后者基于野外实地调查和遥感解译识别；
s03：泥石流的动态危险性分析，泥石流的动态危险性的计算用flow-模型；
s04：泥石流的动态危险性分区。
3.根据权利要求2所述的泥石流实时动态危险性评价方法，其特征在于，所述s01包括以下子步骤：
s011：确定水文参数；
s012：确定物理力学性质参数；
s013：确定土层厚度；
s014：确定实时降雨条件；
s015：斜坡动态稳定性分析。
4.根据权利要求2所述的泥石流实时动态危险性评价方法，其特征在于，所述s02包括以下子步骤：
s021：识别坡面潜在动态物源；
s022：识别沟道已有松散物源。
5.根据权利要求2所述的泥石流实时动态危险性评价方法，其特征在于，所述s03包括以下子步骤：
s031：计算泥石流的动态危险性概率分布；
s032：计算泥石流的动态能量分布。
6.根据权利要求2所述的泥石流实时动态危险性评价方法，其特征在于，所述s04具体为将不同降雨时长的泥石流危险性概率值按分级标准分为极高危险、高度危险、中度危险和低度危险4个等级；其等级标准为：低危险区为0.0003～0.004，中度危险区为0.004～0.032,高度危险区为0.032～0.11，极高危险区为0.11～1。
技术总结本发明公开了一种泥石流实时动态危险性评价方法，包括以下步骤：S1：将降雨数据、岩土体参数输入TRIGRS模型，计算不同降雨历时下斜坡的稳定性系数Fs，并将Fs＜1判别为“是源区”，将Fs≥1判别为“非源区”；S2：将Fs和其它评价因子图层，如汇流累积量、坡度、平面曲率、岩性等，同时输入到Flow‑R模型中，识别出不同降雨历时下泥石流坡面潜在物源区的分布，实现TRIGRS和Flow‑R模型的耦合；同时基于野外实地调查和遥感解译，再识别出泥石流沟道中已有松散物源区的分布；S3:基于Flow‑R模型，计算不同降雨历时下泥石流的能量和危险性概率分布，进行泥石流实时动态危险性分区；该方法解决了以往方法没有全面考虑泥石流物源分布和多侧重于泥石流危险性静态评价研究问题。
技术研发人员：李秀珍;聂银瓶;徐瑞池;姚杰
受保护的技术使用者：中国科学院、水利部成都山地灾害与环境研究所
技术研发日：2020.01.14
技术公布日：2020.06.09
