【技术实现步骤摘要】
本专利技术属于泥石流灾害评价，尤其涉及基于暴发频率和冲出方量的泥石流易发性评价方法、系统、设备及存储介质。
技术介绍
1、泥石流是山区最常见，最具破坏性的地质灾害之一，泥石流不仅严重威胁人们的生命财产安全也制约了山区的经济发展。泥石流易发性评价能反映泥石流可能发生的位置，在泥石流防治决策中有重要地位，也是泥石流领域全球性的研究热点。而选择合适的分析单元是泥石流易发性评价中重要的环节之一。关于泥石流分析单元，前人大多使用网格单元，该分析单元能充分考虑连续型和类型的影响因子。而这样的研究是因为混淆了滑坡和泥石流的概念，网格单元在滑坡易发性评价中能代表足够的信息，但在泥石流易发性评价研究中，流域单元能比网格单元获取更丰富的灾情特征和空间信息。目前已有学者使用流域单元，对不同的易发性评价模型的性能进行了评价和比较，但他们针对泥石流的取值仅限于是与非（1或0），虽然有一定进步，但仍难以体现泥石流的发育程度和灾情信息，凸显泥石流的个体差异。
技术实现思路
1、本专利技术的目的在于：基于暴发频率和冲出方量的泥石流易发性评价方法、系统、设备及存储介质，以解决现有技术难以体现泥石流的发育程度和灾情信息，无法凸显泥石流的个体差异的问题。
2、本申请实施例是这样实现的，提供基于暴发频率和冲出方量的泥石流易发性评价方法，包括：提供影响因子x值；
3、提供泥石流发育程度y值，泥石流发育程度y值的计算公式如下：
4、；
5、式中，为近50年内泥石流的暴发频率，单
6、将影响因子x值和泥石流发育程度y值输入地理探测器-频率比模型获得所有泥石流的易发性指数；
7、将获得的所有泥石流的易发性指数进行易发性等级划分，获得泥石流易发性评价结果。
8、在一些实施例中，影响因子x值包括流域平均坡度，沟道纵比降，流域面积，ndvi，降雨指数，断裂分布线密度，岩性，土地利用类型和相对高差中的一种或多种。
9、在一些实施例中，采用泥石流流域中软岩分布面积在整个流域面积中的占比来代表岩性；和/或采用耕地面积在整个流域面积中的占比代表土地利用类型。
10、在一些实施例中，将影响因子x值和泥石流发育程度y值输入到地理探测器中，提取影响因子x值的权重，且针对泥石流发育程度y值，筛选出主控因子f。
11、在一些实施例中，每个泥石流流域作为分析单元，采用地理探测器模型探测各个影响因子x值对泥石流发育的影响程度，并以q值代表影响因子x值的权重。
12、在一些实施例中，将主控因子f带入频率比模型，得到泥石流的易发性指数。
13、在一些实施例中，将主控因子f划分成n种类型或n个等级，分别是fi（i=1，2，3…，n），则fi的频率比值fri定义为：
14、；
15、式中，表示已知泥石流d中的频率，表示研究区所有拟评价的泥石流中的频率，表示已知泥石流d中的数量，表示泥石流d的总数量，表示研究区所有拟评价的泥石流中的总数量， a表示研究区所有拟评价的泥石流的总数量。
16、针对某一研究区，将该研究区某一泥石流d的不同类型主控因子的频率比相加，就得到该泥石流d的易发性指数:
17、；
18、式中，表示易发性指数，表示主控因子的数量，表示不同类型主控因子的频率比值；
19、同上，依次获得该研究区所有泥石流的易发性指数。
20、相应的，本申请实施例还提供基于暴发频率和冲出方量的泥石流易发性评价系统，包括：
21、影响因子x值模块，用于提供影响因子x值；
22、泥石流发育程度y值模块，用于提供泥石流发育程度y值，泥石流发育程度y值的计算公式如下：
23、；
24、式中，为近50年内泥石流的暴发频率，单位为次/年；为泥石流最近一次暴发冲出方量，单位为m3；
25、泥石流的易发性指数获得模块，用于将影响因子x值和泥石流发育程度y值输入地理探测器-频率比模型，获得所有泥石流的易发性指数；
26、泥石流易发性评价结果获得模块，用于将获得的所有泥石流的易发性指数进行易发性等级划分，获得泥石流易发性评价结果。
27、相应的，本申请实施例还提供计算机设备，包括储存器和处理器，所述储存器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器执行如上述方法的步骤。
28、相应的，本申请实施例还提供计算机可读存储介质，存储有计算机程序，所述计算机程序被处理器执行时，使得所述处理器执行如上述方法的步骤。
29、本申请中：
30、ndvi为归一化植被指数；
31、综上所述，由于采用了上述技术方案，本专利技术的有益效果是：
32、本申请提出利用一次性冲出方量和50年内泥石流的暴发频次的乘积来表征已暴发泥石流沟谷的发育程度，即y值，该y值物理意义为多年年均冲出的固体物质，更能体现泥石流的发育程度；本申请将影响因子x值和泥石流发育程度y值输入地理探测器-频率比模型获得所有泥石流的易发性指数，将获得的所有泥石流的易发性指数进行易发性等级划分，获得泥石流易发性评价结果，本申请获取的高易发性以上流域的空间分布结果更接近实际现场调查的结果，而传统方法为y值的模型获取的易发性结果则高易发流域数量仅有59条，远小于现场调查的结果，这也说明本申请新专利技术的方法较传统方法具有优越性。
本文档来自技高网...【技术保护点】
1.基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，包括：
2.根据权利要求1所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，影响因子X值包括流域平均坡度，沟道纵比降，流域面积，NDVI，降雨指数，断裂分布线密度，岩性，土地利用类型和相对高差中的一种或多种。
3.根据权利要求2所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，采用泥石流流域中软岩分布面积在整个流域面积中的占比来代表岩性；和/或采用耕地面积在整个流域面积中的占比代表土地利用类型。
4.根据权利要求1所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，将影响因子X值和泥石流发育程度Y值输入到地理探测器中，提取影响因子X值的权重，且针对泥石流发育程度Y值，筛选出主控因子F。
5.根据权利要求4所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，每个泥石流流域作为分析单元，采用地理探测器模型探测各个影响因子X值对泥石流发育的影响程度，并以q值代表影响因子X值的权重。
6.根据权利要求4所述的基于暴发
7.根据权利要求6所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，将主控因子F划分成n种类型或n个等级，分别是Fi（i=1，2，3…，n），则Fi的频率比值FRi定义为：
8.基于暴发频率和冲出方量的泥石流易发性评价系统，其特征在于，包括：
9.计算机设备，其特征在于，包括储存器和处理器，所述储存器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器执行如权利要求1-7中任一项所述方法的步骤。
10.计算机可读存储介质，其特征在于，存储有计算机程序，所述计算机程序被处理器执行时，使得所述处理器执行如权利要求1-7任一项所述方法的步骤。
...【技术特征摘要】
1.基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，包括：
2.根据权利要求1所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，影响因子x值包括流域平均坡度，沟道纵比降，流域面积，ndvi，降雨指数，断裂分布线密度，岩性，土地利用类型和相对高差中的一种或多种。
3.根据权利要求2所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，采用泥石流流域中软岩分布面积在整个流域面积中的占比来代表岩性；和/或采用耕地面积在整个流域面积中的占比代表土地利用类型。
4.根据权利要求1所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，将影响因子x值和泥石流发育程度y值输入到地理探测器中，提取影响因子x值的权重，且针对泥石流发育程度y值，筛选出主控因子f。
5.根据权利要求4所述的基于暴发频率和冲出方量的泥石流易发性评价方法，其特征在于，每个泥石流流域作为分析单元，采用...
【专利技术属性】
技术研发人员：张佳佳，李奕霖，田尤，高波，黄海，李元灵，陈龙，谢绮丽，
申请(专利权)人：中国地质科学院探矿工艺研究所，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。