一种山洪灾害风险动态预警方法及装置与流程

1.本技术涉及山洪灾害预报领域，尤其涉及一种山洪灾害风险动态预警方法及装置。


背景技术：

2.山洪灾害是自然灾害中造成人员伤亡的主要灾种，山洪灾害突发性强，破坏力大，且多发生于山丘区，涨洪历时短、流域汇流块、洪水陡涨陡落，因此预报预警难度很大。
3.目前山洪灾害预警主要包括实测雨量(水位、流量)预警、洪水预报预警、山洪灾害气象预报预警。采用实测雨量预警的难点在于如何科学、合理地确定预警指标，以及如何降低土壤湿度对预警准确性的影响，而采用实测水位和流量更加合理，但是目前能够监测水位和流量的遥测站点较少，不能有效覆盖山丘区；但由于山丘区环境的复杂性和山洪形成的多要素性，雨量和水位等预警指标的分析过程中对土壤含水量考虑不足，导致基于实测数据山洪预警存在精确度不高，预见期过短的问题。洪水预报预警是指通过水文模型以实测和预报的降雨数据为输入条件，对降雨径流过程进行模拟，将降雨数据转化为水位和流量数据进行预警。洪水模拟可以有效延长预警的预见期，但是由于模型参数确定困难，因此预报预警的准确性很难得到保证。山洪灾害气象预报预警是指以气象预报降雨为基础进行预警，预警方式与实测降雨预警类似，其优点是预见期延长，缺点是预警准确度不高。
4.现有的山洪灾害预警方法均存在预见期短或预警准确性低的问题，导致现有的山洪灾害预警方法不能满足山洪防御工作需要。


技术实现要素：

5.本技术提供了一种山洪灾害风险动态预警方法及装置，用于延长山洪灾害预警的预见期及提高山洪灾害预警的准确性。
6.本技术第一方面提供了一种山洪灾害风险动态预警方法，包括：
7.将预警区域均匀划分为若干网格，并获取每个网格的风险因子及风险等级,所述风险因子为影响所述山洪灾害风险的因素；
8.根据所述风险因子确定所述每个网格的实时土壤含水量；
9.根据所述实时土壤含水量确定所述每个网格的网格土壤下渗能力；
10.根据所述风险因子确定致灾因子及所述致灾因子对应的致灾权重，所述致灾因子包括：危险性因子、承灾体因子及易损性因子；
11.根据所述致灾因子及所述致灾权重确定所述每个网格的风险指数；
12.根据所述风险指数确定所述每个网格的不同预警时段及不同所述风险等级的网格基准指标；
13.根据所述网格土壤下渗能力及所述网格基准指标确定所述每个网格的风险预警指标；
14.获取所述每个网格的数值预报降雨插值，并根据所述数值预报降雨插值统计所述
每个网格在所述不同预警时段的最大降雨量；
15.比较所述最大降雨量及所述风险预警指标，确定所述每个网格在所述不同预警时段的风险预警等级。
16.可选地，在所述比较所述最大降雨量及所述风险预警指标，确定所述每个网格在所述不同预警时段内的风险预警等级之后，所述山洪灾害风险动态预警方法还包括：
17.将所述风险预警等级与风险实际等级进行对比，得到对比结果；
18.根据所述对比结果评估风险预警效果。
19.可选地，在所述根据所述风险因子确定所述每个网格的实时土壤含水量之前，所述山洪灾害动态预警方法还包括：
20.根据所述风险因子中的土壤类型及平均坡度通过第一公式确定所述每个网格的最大土壤含水量；
21.所述第一公式为：
22.w
mi
＝wm(1-sin ai)
0.7
；
23.其中，w
mi
为第i个网格的最大土壤含水量，wm为根据所述土壤类型确定的最大土壤含水量初值，ai为第i个网格的平均坡度；
24.所述根据所述风险因子确定所述每个网格的实时土壤含水量包括：
25.根据所述每个网格的最大土壤含水量确定所述每个网格的实时土壤含水量。
26.可选地，所述根据所述每个网格的最大土壤含水量确定所述每个网格的实时土壤含水量包括：
27.根据所述每个网格的最大土壤含水量通过第二公式确定所述每个网格的实时土壤含水量；
28.所述第二公式为：
29.p
a，t
＝k1(p
t-1
+p
a，t-1
)；
[0030][0031]
其中，p
a，t
为第i个网格的实时土壤含水量、p
t-1
为t-1时刻的降雨量、p
a，t-1
为t-1时刻的第i个网格的土壤含水量，当p
a，t
≥w
mi
时，取p
a，t
＝w
mi
，w
mi
为第i个网格的最大土壤含水量、k1为日消退系数、em为目标气象站每个月的平均日最大蒸散发量。
[0032]
可选地，所述根据所述实时土壤含水量确定所述每个网格的网格土壤下渗能力包括：
[0033]
根据实时土壤含水量通过第三公式确定所述每个网格的网格土壤下渗能力；
[0034]
所述第三公式为：
[0035]
f(p
a，t
)＝fc+f0e-μm
；
[0036]
m＝f
t
+k2p
a，t
；
[0037]
其中，f(p
a，t
)为每个网格的网格土壤下渗能力，fc为稳定下渗率、f0为初始下渗率、μ为下渗指数、m为表层土湿、k2为土壤表层厚度与包气带厚度的比值，取值范围为[0，1]、p
a，t
为第i个网格的土壤实时含水量、f
t
为累计下渗量。
[0038]
可选地，所述根据所述致灾因子及所述致灾权重确定所述每个网格的风险指数包括：
[0039]
根据所述致灾因子及所述致灾权重通过第四公式确定所述每个网格的风险指数；
[0040]
所述第四公式为：
[0041][0042]
其中，risk为风险指数、h为危险性因子、e为承灾体因子、v为易损性因子、hi为危险性因子的归一化指标、ej为承灾体因子的归一化指标、vk为易损性因子的归一化指标、m为危险性因子的指标数量、n为承灾体因子的指标数量、l为易损性因子的指标数量、wi为危险性因子的致灾权重、wj为承灾体因子的致灾权重、wk为易损性因子的致灾权重。
[0043]
可选地，所述根据所述风险指数确定所述每个网格的不同预警时段及不同风险等级的网格基准指标包括：
[0044]
对所述风险指数risk进行降序排列，获取所述每个网格的秩i
risk
；
[0045]
根据所述每个网格的秩i
risk
确定所述每个网格不同预警时段及不同风险等级的网格基准指标包括：
[0046]
根据所述每个网格的秩i
risk
通过第五公式确定所述每个网格不同预警时段及不同风险等级的网格基准指标；
[0047]
所述第五公式为：
[0048][0049]
其中，g为风险等级、h为预警时段、n为网格点总数、h
min
(g，h)为风险等级及预警时段的基准指标下限、h
max
(g，h)为风险等级及预警时段的基准指标上限。
[0050]
可选地，所述根据所述网格土壤下渗能力及所述网格基准指标确定所述每个网格的风险预警指标包括：
[0051]
根据所述网格土壤下渗能力及所述网格基准指标通过第六公式确定所述每个网格的风险预警指标；
[0052]
所述第六公式为：
[0053]r′
(g，h)＝f(p
a，t
)+r(risk，g，h)；
[0054]
其中，r
′
(g，h)为风险预警指标、f(p
a，t
)为土壤下渗能力、r(risk，g，h)为风险预警基准指标、g为风险等级、h为预警时段、risk为网格风险指数、p
a，t
为第i个网格的实时土壤含水量。
[0055]
可选地，所述获取所述每个网格的数值预报降雨插值包括：
[0056]
根据气象站点的数值预报降雨插值通过第七公式获取所述每个网格的数值预报降雨插值；
[0057]
所述第七公式为：
[0058][0059]
其中，z为每个网格的数值预报降雨插值、zi(i＝1、2、
…
n)为气象站点实测的若干个数值预报降雨插值、di为插值点与第i个气象站点之间的距离、p为距离的幂。
[0060]
本技术第二方面提供了一种山洪灾害风险动态预警装置，包括：
[0061]
划分单元，用于将预警区域均匀划分为若干网格，并获取每个网格的风险因子及
风险等级,所述风险因子为影响所述山洪灾害风险的因素；
[0062]
第一确定单元，用于根据所述风险因子确定所述每个网格的实时土壤含水量；
[0063]
第二确定单元，用于根据所述实时土壤含水量确定所述每个网格的网格土壤下渗能力；
[0064]
第三确定单元，用于根据所述风险因子确定致灾因子及所述致灾因子对应的致灾权重，所述致灾因子包括：危险性因子、承灾体因子及易损性因子；
[0065]
第四确定单元，用于根据所述致灾因子及所述致灾权重确定所述每个网格的风险指数；
[0066]
第五确定单元，用于根据所述风险指数确定所述每个网格的不同预警时段及不同所述风险等级的网格基准指标；
[0067]
第六确定单元，用于根据所述网格土壤下渗能力及所述网格基准指标确定所述每个网格的风险预警指标；
[0068]
获取单元，用于获取所述每个网格的数值预报降雨插值，并根据所述数值预报降雨插值统计所述每个网格在所述不同预警时段内的最大降雨量；
[0069]
比较单元，用于比较所述最大降雨量及所述风险预警指标，确定所述每个网格在所述不同预警时段内的风险预警等级。
[0070]
可选地，所述山洪灾害风险动态预警装置还包括：
[0071]
对比单元，用于将所述风险预警等级与风险实际等级进行对比，得到对比结果；
[0072]
评估单元，用于根据所述对比结果评估风险预警效果。
[0073]
可选地，所述洪灾害风险动态预警装置还包括：
[0074]
第七确定单元，用于根据所述风险因子中的土壤类型及平均坡度通过第一公式确定所述每个网格的最大土壤含水量；
[0075]
所述第一公式为：
[0076]wmi
＝wm(1-sin ai)
0.7
；
[0077]
其中，w
mi
为第i个网格的最大土壤含水量，wm为根据所述土壤类型确定的最大土壤含水量初值，ai为第i个网格的平均坡度；
[0078]
所述第一确定单元具体用于：
[0079]
根据所述每个网格的最大土壤含水量确定所述每个网格的实时土壤含水量。
[0080]
可选地，所述第一确定单元具体用于：
[0081]
根据所述每个网格的最大土壤含水量通过第二公式确定所述每个网格的实时土壤含水量；
[0082]
所述第二公式为：
[0083]
p
a，t
＝k1(p
t-1
+p
a，t-1
)；
[0084][0085]
其中，p
a，t
为第i个网格的实时土壤含水量、p
t-1
为t-1时刻的降雨量、p
a，t-1
为t-1时刻的第i个网格的土壤含水量，当p
a，t
≥w
mi
时，取p
a，t
＝w
mi
，w
mi
为第i个网格的最大土壤含水量、k1为日消退系数、em为目标气象站每个月的平均日最大蒸散发量。
[0086]
可选地，所述第二确定单元具体用于：
[0087]
根据实时土壤含水量通过第三公式确定所述每个网格的网格土壤下渗能力；
[0088]
所述第三公式为：
[0089]
f(p
a，t
)＝fc+f0e-μm
；
[0090]
m＝f
t
+k2p
a，t
；
[0091]
其中，f(p
a，t
)为每个网格的网格土壤下渗能力，fc为稳定下渗率、f0为初始下渗率、μ为下渗指数、m为表层土湿、k2为土壤表层厚度与包气带厚度的比值，取值范围为[0，1]、p
a，t
为第i个网格的土壤实时含水量、f
t
为累计下渗量。
[0092]
可选地，所述第四确定单元具体用于：
[0093]
根据所述致灾因子及所述致灾权重通过第四公式确定所述每个网格的风险指数；
[0094]
所述第四公式为：
[0095][0096]
其中，risk为风险指数、h为危险性因子、e为承灾体因子、v为易损性因子、hi为危险性因子的归一化指标、ej为承灾体因子的归一化指标、vk为易损性因子的归一化指标、m为危险性因子的指标数量、n为承灾体因子的指标数量、l为易损性因子的指标数量、wi为危险性因子的致灾权重、wj为承灾体因子的致灾权重、wk为易损性因子的致灾权重。
[0097]
可选地，所述第五确定单元具体用于：
[0098]
对所述风险指数risk进行降序排列，获取所述每个网格的秩i
risk
；
[0099]
根据所述每个网格的秩i
risk
确定所述每个网格不同预警时段及不同风险等级的网格基准指标包括：
[0100]
根据所述每个网格的秩i
risk
通过第五公式确定所述每个网格不同预警时段及不同风险等级的网格基准指标；
[0101]
所述第五公式为：
[0102][0103]
其中，g为风险等级、h为预警时段、n为网格点总数、h
min
(g，h)为风险等级及预警时段的基准指标下限、h
max
(g，h)为风险等级及预警时段的基准指标上限。
[0104]
可选地，所述第六确定单元具体用于：
[0105]
根据所述网格土壤下渗能力及所述网格基准指标通过第六公式确定所述每个网格的风险预警指标；
[0106]
所述第六公式为：
[0107]r′
(g，h)＝f(p
a，t
)+r(risk，g，h)；
[0108]
其中，r
′
(g，h)为风险预警指标、f(p
a，t
)为土壤下渗能力、r(risk，g，h)为风险预警基准指标、g为风险等级、h为预警时段、risk为网格风险指数、p
a，t
为第i个网格的实时土壤含水量。
[0109]
可选地，所述获取单元具体用于：
[0110]
根据气象站点的数值预报降雨插值通过第七公式获取所述每个网格的数值预报降雨插值；
[0111]
所述第七公式为：
[0112][0113]
其中，z为每个网格的数值预报降雨插值、zi(i＝1、2、
…
n)为气象站点实测的若干个数值预报降雨插值、di为插值点与第i个气象站点之间的距离、p为距离的幂。
[0114]
本技术第三方面提供了一种山洪灾害风险动态预警装置，所述装置包括：
[0115]
处理器、存储器、输入输出单元以及总线；
[0116]
所述处理器与所述存储器、所述输入输出单元以及所述总线相连；
[0117]
所述存储器保存有程序，所述处理器调用所述程序以执行第一方面以及第一方面中任一项可选的一种山洪灾害风险动态预警方法。
[0118]
本技术第四方面提供了一种计算机可读存储介质，所述计算机可读存储介质上保存有程序，所述程序在计算机上执行时执行第一方面以及第一方面中任一项可选的一种山洪灾害风险动态预警方法。
[0119]
从以上技术方案可以看出，本技术具有以下优点：本技术通过每个网格的数值预报降雨插值进行山洪预警，有效延长了山洪灾害预警的预见期，从而可以最大程度地为群众安全转移争取时间，减少人员伤亡和财产损失，并且，在制定风险预警指标时，充分考虑了每个网格的实时土壤含水量以及致灾因子对风险预警指标的影响，有效提高了山洪灾害预警的准确性。
附图说明
[0120]
为了更清楚地说明本技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0121]
图1为本技术提供的一种山洪灾害风险动态预警方法一个实施例流程示意图；
[0122]
图2为本技术提供的一种山洪灾害风险动态预警方法另一个实施例流程示意图；
[0123]
图3为本技术提供的一种山洪灾害风险动态预警装置一个实施例结构示意图；
[0124]
图4为本技术提供的一种山洪灾害风险动态预警装置另一个实施例结构示意图；
[0125]
图5为本技术提供的一种山洪灾害风险动态预警装置一个实施例结构示意图；
[0126]
图6为本技术提供的基准指标上下推荐表；
[0127]
图7为本技术提供的判断矩阵比较量化表。
具体实施方式
[0128]
本技术提供了一种山洪灾害风险动态预警方法，用于延长山洪灾害预警的预见期及提高山洪灾害预警的准确性。
[0129]
需要说明的是，本技术提供的一种山洪灾害风险动态预警方法，可以应用于终端，还可以应用于服务器上，例如终端可以是智能手机或电脑、平板电脑、智能电视、智能手表、便携计算机终端也可以是台式计算机等固定终端。为方便阐述，本技术中以终端为执行主体进行举例说明。
[0130]
请参阅图1，图1为本技术提供的一种山洪灾害风险动态预警方法的一个实施例，
该山洪灾害风险动态预警方法包括：
[0131]
101、终端将预警区域均匀划分为若干网格，并获取每个网格的风险因子及风险等级，风险因子为影响山洪灾害风险的因素；
[0132]
在本实施例中，终端根据基础数据的密度和分辨率将预警区域均匀划分为若干网格，并以每个网格为风险分析对象，提取每个网格的风险因子及风险等级，该风险为子为影响山洪灾害风险的因素，该风险因子包括：土壤类型、平均坡度、防洪能力、不同频率影响人口及不同时长设计暴雨值；该风险等级包括：蓝色风险等级、黄色风险等级、橙色风险等级及红色风险等级四个风险等级，该基础数据为预警区域的图片或照片。例如：终端根据基础数据的密度和分辨率将预警区域均匀划分为0.05
°×
0.05
°
的若干个网格。
[0133]
102、终端根据风险因子确定每个网格的实时土壤含水量；
[0134]
在本实施例中，终端根据风险因子中的土壤类型及平均坡度确定每个网格的实时土壤含水量，具体的确定过程将在下一个实施例中进行详细描述，此处不再赘述。
[0135]
103、终端根据实时土壤含水量确定每个网格的网格土壤下渗能力；
[0136]
在本实施例中，终端获取每个网格的稳定下渗率、初始下渗率、累计下渗量、下渗指数、表层土湿、土壤表层厚度与包气带厚度的比值，并根据实时土壤含水量及每个网格的稳定下渗率、初始下渗率、累计下渗量、下渗指数、表层土湿、土壤表层厚度与包气带厚度的比值确定每个网格的网格土壤下渗能力。具体的确定过程将在下一个实施例中进行详细描述，此处不再赘述。
[0137]
104、终端根据风险因子确定致灾因子及致灾因子对应的致灾权重，致灾因子包括：危险性因子、承灾体因子及易损性因子；
[0138]
在本实施例中，终端根据风险因子筛选致灾因子并确定每个致灾因子对应的致灾权重，致灾因子包括危险性因子、承灾体因子及易损性因子。例如：危险性因子包括：1h百年一遇设计暴雨值，其对应的致灾权重为0.2970、6h百年一遇设计暴雨值，其对应的致灾权重为0.5396、或平均坡度，其对应的致灾权重为0.1276。承灾体因子包括：5年一遇洪水影响人口，其对应的致灾权重为0.5539、20年一遇洪水影响人口，其对应的致灾权重为0.3185、或100年一遇洪水影响人口，其对应的致灾权重为0.1276。易损性因子包括：防洪能力，其对应的致灾权重为1。
[0139]
105、终端根据致灾因子及致灾权重确定每个网格的风险指数；
[0140]
在本实施例中，终端根据每个网格的致灾因子及致灾权重确定每个网格的风险指数，具体的确定过程将在下一个实施例中进行详细描述，此处不再赘述。
[0141]
106、终端根据风险指数确定每个网格的不同预警时段及不同风险等级的网格基准指标；
[0142]
在本实施例中，终端将每个网格的风险指数进行降序排列，并获取每个网格的秩，终端根据每个网格的秩确定每个网格的不同预警时段及不同风险等级的网格基准指标。不同预警时段可以包括1h、3h、6h、24h等预警时段，不同风险等级可以包括蓝色风险等级、黄色风险等级、橙色风险等级及红色风险等级等风险等级。
[0143]
107、终端根据网格土壤下渗能力及网格基准指标确定每个网格的风险预警指标；
[0144]
在本实施例中，终端根据网格土壤下渗能力及网格基准指标确定每个网格的不同预警时段及不同风险等级的风险预警指标。不同预警时段可以包括1h、3h、6h、24h等预警时
段，不同风险等级可以包括蓝色风险等级、黄色风险等级、橙色风险等级及红色风险等级等风险等级。
[0145]
108、终端获取每个网格的数值预报降雨插值，并根据数值预报降雨插值统计每个网格在不同预警时段的最大降雨量；
[0146]
在本实施例中，终端根据气象站点的数值预报降雨插值通过反距离权重模型计算每个网格的数值预报降雨插值，并根据每个网格的数值预报降雨插值统计每个网格在不同预警时段的最大降雨量。
[0147]
109、终端比较最大降雨量及风险预警指标，确定每个网格在不同预警时段的风险预警等级。
[0148]
在本实施例中，终端比较每个网格的最大降雨量及风险预警指标,确定每个网格在不同预警时段的风险预警等级。
[0149]
在本实施例中，终端通过每个网格的数值预报降雨插值进行山洪预警，有效延长了预警的预见期，从而可以最大程度地为群众安全转移争取时间，减少人员伤亡和财产损失，并且，在制定风险预警指标时，充分考虑了每个网格的实时土壤含水量以及致灾因子对风险预警指标的影响，有效提高了预警的准确性。
[0150]
为使本技术提供的一种山洪灾害风险动态预警方法更加的明显易懂，下面对本技术提供的一种山洪灾害风险动态预警方法进行详细说明：
[0151]
请参阅图2，图2为本技术提供的一种山洪灾害风险动态预警方法另一个实施例，该山洪灾害风险动态预警方法包括：
[0152]
201、终端将预警区域均匀划分为若干网格，并获取每个网格的风险因子及风险等级，该风险因子为影响山洪灾害风险的因素；
[0153]
本实施例中的步骤201与前述图1实施例中的步骤101类似，具体此处不再赘述。
[0154]
202、终端根据风险因子中的土壤类型及平均坡度确定每个网格的最大土壤含水量；
[0155]
在本实施例中，终端根据风险因子中的土壤类型及平均坡度通过第一公式确定每个网格的最大土壤含水量，第一公式为：w
mi
＝wm(1-sin ai)
0.7
，其中，w
mi
为第i个网格的最大土壤含水量，wm为根据所述土壤类型确定的最大土壤含水量初值，ai为第i个网格的平均坡度。
[0156]
203、终端根据每个网格的最大土壤含水量确定每个网格的实时土壤含水量；
[0157]
在本实施例中，终端根据每个网格的最大土壤含水量通过第二公式确定每个网格的实时土壤含水量，第二公式为：
[0158]
p
a，t
＝k1(p
t-1
+p
a，t-1
)，
[0159][0160]
其中，p
a，t
为第i个网格的实时土壤含水量、p
t-1
为t-1时刻的降雨量、p
a，t-1
为t-1时刻的第i个网格的土壤含水量，当p
a，t
≥w
mi
时，取p
a，t
＝w
mi
，w
mi
为第i个网格的最大土壤含水量；k1为日消退系数、em为目标气象站每个月的平均日最大蒸散发量。
[0161]
204、终端根据实时土壤含水量确定每个网格的网格土壤下渗能力；
[0162]
在本实施例中，终端根据实时土壤含水量确定每个网格的网格土壤下渗能力包
括：终端根据实时土壤含水量通过第三公式确定每个网格的网格土壤下渗能力，第三公式为：
[0163]
m＝f
t
+k2p
a，t
；
[0164]
f(p
a，t
)＝fc+f0e-μm
；
[0165]
其中，f(p
a，t
)为每个网格的网格土壤下渗能力，fc为稳定下渗率、f0为初始下渗率、μ为下渗指数、m为表层土湿、k2为土壤表层厚度与包气带厚度的比值，取值范围为[0，1]、p
a，t
为第i个网格的土壤实时含水量、f
t
为累计下渗量。
[0166]
205、终端根据风险因子确定致灾因子及致灾因子对应的致灾权重，致灾因子包括：危险性因子、承灾体因子及易损性因子；
[0167]
在本实施例中，终端根据风险因子确定致灾因子及致灾因子对应的致灾权重包括：终端利用成分分析法对风险因子进行降维处理，筛选出致灾因子，并确定致灾因子对应的致灾权重。例如：风险因子包括：土壤类型、平均坡度、防洪能力、不同频率洪水影响人口及不同时长设计暴雨值；终端通过成分分析法筛选出危险性因子为：1h百年一遇设计暴雨值(权重为0.2970)、6h百年一遇设计暴雨值(权重为0.5396)以及平均坡度(权重为0.1634)；终端通过成分分析法筛选出承灾体因子为：5年一遇洪水影响人口(权重为0.5539)、20年一遇洪水影响人口(权重为0.3185)以及100年一遇洪水影响人口(权重为0.1276)；终端通过成分分析法筛选出易损性因子为：防洪能力(权重为1)。
[0168]
在本实施例中，本技术通过一致矩阵法确定各层次各因素之间的权重，即不把所有因素放在一起比较，而是两两相互比较，采用相对尺度，尽可能减少性质不同的诸因素相互比较的困难，以提高准确度。如对某一准则，对其下的各方案进行两两对比，并按其重要性程度评定等级。aij为要素i与要素j重要性比较结果。按两两比较结果构成的矩阵称作判断矩阵。判断矩阵具有如下性质：
[0169][0170]
判断矩阵元素aij通过判断矩阵比较量化表(图7)进行标度。例如：构造危险性因子的判断矩阵如下所示，求解λmax＝3.0092,ci＝0.0046,cr《0.1，通过一致性检验，为可接受矩阵。其中1h百年一遇设计暴雨值、6h百年一遇设计暴雨值和平均坡度的权重系数分别为0.2970、0.5396和0.1634。构造承灾体的判断矩阵如下所示，求解λmax＝3.0713,ci＝0.0357,cr《0.1，通过一致性检验，为可接受矩阵。其中5年一遇、20年一遇、100年一遇洪水位下的人口数量的权重系数分别为0.5539、0.3185和0.1276。而由于易损性因子只有防洪能力这一项指标，因此防洪能力的权重为1。
[0171]
在本实施例中，采用主成分分析(pca)法通过选择不相关的度量，表示最自然的风险特征类别簇集。pca是一种统计步骤，使用正交变换来转换一组可能相关变量的观测值，到被称为主成分(pcs)的线性不相关变量的一组值。这种转换是这样定义的，第一主成分有最大的方差(因此，尽可能地解释数据中的可变性)，每个随后的成分依次是最高的，在方差可能下的约束，它是正交于前面的成分。在这项研究中，对于风险特征的每个类别，共同解释了至少每个风险特征类别方差的三分之二(66.6％)的主成分被保留。如果属于给定类别的变量不都有相同的单位，在分析发生之前，它们被缩放成单位方差。pca可以显著减少风险特征分析的要素数量，降低风险分析复杂度。
[0172]
206、终端根据致灾因子及致灾权重确定每个网格的风险指数；
[0173]
在本实施例中，终端根据致灾因子及致灾权重通过第四公式确定每个网格的风险指数，第四公式为：
[0174][0175]
其中，risk为风险指数、h为危险性因子、e为承灾体因子、v为易损性因子、hi为危险性因子的归一化指标、ej为承灾体因子的归一化指标、vk为易损性因子的归一化指标、m为危险性因子的指标数量、n为承灾体因子的指标数量、l为易损性因子的指标数量、wi为危险性因子的致灾权重、wj为承灾体因子的致灾权重、wk为易损性因子的致灾权重。
[0176]
207、终端根据风险指数确定每个网格的不同预警时段及不同风险等级的网格基准指标；
[0177]
在本实施例中，终端根据风险指数确定每个网格的不同预警时段及不同风险等级的网格基准指标包括：对每个网格的风险指数risk进行降序排列，获取每个网格的秩i
risk
；根据每个网格的秩i
risk
通过第五公式确定所述每个网格不同预警时段及不同风险等级的网格基准指标，第五公式为：
[0178][0179]
其中，g为风险等级、h为预警时段、n为网格点总数、h
min
(g，h)为风险等级及预警时段的基准指标下限、h
max
(g，h)为风险等级及预警时段的基准指标上限。
[0180]
在本实施例中，通过查询基准指标上下限推荐表(图6)获取hmin(g，h)及hmax(g，h)，例如g为3h，h为黄色预警等级，则此时通过图6查询到hmin(1h，黄色)＝多年降雨均值*0.8、hmax(1h，黄色)＝多年降雨均值；将hmin(1h，黄色)及hmax(1h，黄色)代入得到网格的r(risk，1h，黄色)的值。
[0181]
208、终端根据网格土壤下渗能力及网格基准指标确定每个网格的风险预警指标；
[0182]
在本实施例中，终端根据网格土壤下渗能力及网格基准指标通过第六公式确定每个网格的风险预警指标，第六公式为：
[0183]r′
(g，h)＝f(p
a，t
)+r(risk，g，h)；
[0184]
其中，r
′
(g，h)为风险预警指标、f(p
a，t
)为土壤下渗能力、r(risk，g，h)为风险预警基准指标、g为风险等级、h为预警时段、risk为网格风险指数、p
a，t
为第i个网格的实时土壤含水量。
[0185]
209、终端获取每个网格的数值预报降雨插值，并根据数值预报降雨插值统计每个网格在不同预警时段额最大降雨量；
[0186]
在本实施例中，终端根据气象站点的数值预报降雨插值通过第七公式获取每个网格的数值预报降雨插值，第七公式为：
[0187][0188]
其中，z为每个网格的数值预报降雨插值、zi(i＝1、2、
…
n)为气象站点实测的若干个数值预报降雨插值、di为插值点与第i个气象站点之间的距离、p为距离的幂。
[0189]
在本实施例中，距离平方反比根据插值点与样本点之间的距离计算样本点的权重系数，离插值点越近的样本点其权重越大，权重贡献与距离负相关。p是距离的幂，它显著影响插值计算的结果，在本实施例中，p取2。
[0190]
210、终端比较最大降雨量及风险预警指标，确定每个网格在不同预警时段的风险预警等级；
[0191]
本实施例中的步骤210与前述图1实施例中的步骤109类似，具体此处不再赘述。
[0192]
211、终端将风险预警等级与风险实际等级进行对比，得到对比结果；
[0193]
在本实施例中，例如终端确定某个网格在3h时的风险预警等级为红色，在6h时的风险预警等级为红色，在24h时的风险预警等级为红色，当实际灾害发生时，终端获取该网格在3h时的风险实际等级为红色、在6h时的风险实际等级为红色，在24h时的风险实际等级为红色，则终端得到的对比结果为风险预警等级与风险实际等级完全相同。
[0194]
212、终端根据对比结果评估风险预警效果。
[0195]
在本实施例中，当对比结果大于预设阈值时，终端确定风险预警效果良好。例如：预设阈值为95％，对比结果中风险预警等级与风险实际等级的相同的概率为97％，则终端确定风险预警效果良好。根据对比结果评估风险预警效果，可以迅速发现风险预警效果是否良好，当风险预警效果不好时，可以及时调整预警指标，从而提高风险预警效果。
[0196]
上述对本技术提供的一种山洪灾害风险动态预警方法进行了说明，下面对本技术提供的一种山洪灾害风险动态预警装置进行说明：
[0197]
请参阅图3，图3为本技术提供的一种山洪灾害风险动态预警装置一个实施例，该山洪灾害风险动态预警装置包括：
[0198]
划分单元301，用于将预警区域均匀划分为若干网格，并获取每个网格的风险因子及风险等级,所述风险因子为影响所述山洪灾害风险的因素；
[0199]
第一确定单元302，用于根据所述风险因子确定所述每个网格的实时土壤含水量；
[0200]
第二确定单元303，用于根据所述实时土壤含水量确定所述每个网格的网格土壤下渗能力；
[0201]
第三确定单元304，用于根据所述风险因子确定致灾因子及所述致灾因子对应的致灾权重，所述致灾因子包括：危险性因子、承灾体因子及易损性因子；
[0202]
第四确定单元305，用于根据所述致灾因子及所述致灾权重确定所述每个网格的风险指数；
[0203]
第五确定单元306，用于根据所述风险指数确定所述每个网格的不同预警时段及不同所述风险等级的网格基准指标；
[0204]
第六确定单元307，用于根据所述网格土壤下渗能力及所述网格基准指标确定所述每个网格的风险预警指标；
[0205]
获取单元308，用于获取所述每个网格的数值预报降雨插值，并根据所述数值预报降雨插值统计所述每个网格在所述不同预警时段内的最大降雨量；
[0206]
比较单元309，用于比较所述最大降雨量及所述风险预警指标，确定所述每个网格在所述不同预警时段内的风险预警等级。
[0207]
本实施例系统中，各单元所执行的功能与前述图1所示方法实施例中的步骤对应，具体此处不再赘述。
[0208]
下面对本技术提供的一种山洪灾害风险动态预警装置进行详细说明，请参阅图4，图4为本技术提供的一种山洪灾害风险动态预警装置另一个实施例，该山洪灾害风险动态预警装置包括：
[0209]
划分单元401，用于将预警区域均匀划分为若干网格，并获取每个网格的风险因子及风险等级,所述风险因子为影响所述山洪灾害风险的因素；
[0210]
第一确定单元402，用于根据所述风险因子确定所述每个网格的实时土壤含水量；
[0211]
第二确定单元403，用于根据所述实时土壤含水量确定所述每个网格的网格土壤下渗能力；
[0212]
第三确定单元404，用于根据所述风险因子确定致灾因子及所述致灾因子对应的致灾权重，所述致灾因子包括：危险性因子、承灾体因子及易损性因子；
[0213]
第四确定单元405，用于根据所述致灾因子及所述致灾权重确定所述每个网格的风险指数；
[0214]
第五确定单元406，用于根据所述风险指数确定所述每个网格的不同预警时段及不同所述风险等级的网格基准指标；
[0215]
第六确定单元407，用于根据所述网格土壤下渗能力及所述网格基准指标确定所述每个网格的风险预警指标；
[0216]
获取单元408，用于获取所述每个网格的数值预报降雨插值，并根据所述数值预报降雨插值统计所述每个网格在所述不同预警时段内的最大降雨量；
[0217]
比较单元409，用于比较所述最大降雨量及所述风险预警指标，确定所述每个网格在所述不同预警时段内的风险预警等级。
[0218]
可选地，所述山洪灾害风险动态预警装置还包括：
[0219]
对比单元410，用于将所述风险预警等级与风险实际等级进行对比，得到对比结果；
[0220]
评估单元411，用于根据所述对比结果评估风险预警效果。
[0221]
可选地，所述洪灾害风险动态预警装置还包括：
[0222]
第七确定单元412，用于根据所述风险因子中的土壤类型及平均坡度通过第一公式确定所述每个网格的最大土壤含水量；
[0223]
所述第一公式为：
[0224]wmi
＝wm(1-sin ai)
0.7
；
[0225]
其中，w
mi
为第i个网格的最大土壤含水量，wm为根据所述土壤类型确定的最大土壤含水量初值，ai为第i个网格的平均坡度；
[0226]
所述第一确定单元402具体用于：
[0227]
根据所述每个网格的最大土壤含水量确定所述每个网格的实时土壤含水量。
[0228]
可选地，所述第一确定单元402具体用于：
[0229]
根据所述每个网格的最大土壤含水量通过第二公式确定所述每个网格的实时土壤含水量；
[0230]
所述第二公式为：
[0231]
p
a，t
＝k1(p
t-1
+p
a，t-1
)；
[0232][0233]
其中，p
a，t
为第i个网格的实时土壤含水量、p
t-1
为t-1时刻的降雨量、p
a，t-1
为t-1时刻的第i个网格的土壤含水量，当p
a，t
≥w
mi
时，取p
a，t
＝w
mi
，w
mi
为第i个网格的最大土壤含水量、k1为日消退系数、em为目标气象站每个月的平均日最大蒸散发量。
[0234]
可选地，所述第二确定单元403具体用于：
[0235]
根据实时土壤含水量通过第三公式确定所述每个网格的网格土壤下渗能力；
[0236]
所述第三公式为：
[0237]
f(p
a，t
)＝fc+f0e-μm
；
[0238]
m＝f
t
+k2p
a，t
；
[0239]
其中，f(p
a，t
)为每个网格的网格土壤下渗能力，fc为稳定下渗率、f0为初始下渗率、μ为下渗指数、m为表层土湿、k2为土壤表层厚度与包气带厚度的比值，取值范围为[0，1]、p
a，t
为第i个网格的土壤实时含水量、f
t
为累计下渗量。
[0240]
可选地，所述第四确定单元405具体用于：
[0241]
根据所述致灾因子及所述致灾权重通过第四公式确定所述每个网格的风险指数；
[0242]
所述第四公式为：
[0243][0244]
其中，risk为风险指数、h为危险性因子、e为承灾体因子、v为易损性因子、hi为危险性因子的归一化指标、ej为承灾体因子的归一化指标、vk为易损性因子的归一化指标、m为危险性因子的指标数量、n为承灾体因子的指标数量、l为易损性因子的指标数量、wi为危险性因子的致灾权重、wj为承灾体因子的致灾权重、wk为易损性因子的致灾权重。
[0245]
可选地，所述第五确定单元406具体用于：
[0246]
对所述风险指数risk进行降序排列，获取所述每个网格的秩i
risk
；
[0247]
根据所述每个网格的秩i
risk
确定所述每个网格不同预警时段及不同风险等级的网格基准指标包括：
[0248]
根据所述每个网格的秩i
risk
通过第五公式确定所述每个网格不同预警时段及不同风险等级的网格基准指标；
[0249]
所述第五公式为：
[0250][0251]
其中，g为风险等级、h为预警时段、n为网格点总数、h
min
(g，h)为风险等级及预警时段的基准指标下限、h
max
(g，h)为风险等级及预警时段的基准指标上限。
[0252]
可选地，所述第六确定单元407具体用于：
[0253]
根据所述网格土壤下渗能力及所述网格基准指标通过第六公式确定所述每个网格的风险预警指标；
[0254]
所述第六公式为：
[0255]r′
(g，h)＝f(p
a，t
)+r(risk，g，h)；
[0256]
其中，r
′
(g，h)为风险预警指标、f(p
a，t
)为土壤下渗能力、r(risk，g，h)为风险预警基准指标、g为风险等级、h为预警时段、risk为网格风险指数、p
a，t
为第i个网格的实时土壤
含水量。
[0257]
可选地，所述获取单元408具体用于：
[0258]
根据气象站点的数值预报降雨插值通过第七公式获取所述每个网格的数值预报降雨插值；
[0259]
所述第七公式为：
[0260][0261]
其中，z为每个网格的数值预报降雨插值、zi(i＝1、2、
…
n)为气象站点实测的若干个数值预报降雨插值、di为插值点与第i个气象站点之间的距离、p为距离的幂。
[0262]
本实施例系统中，各单元所执行的功能与前述图2所示方法实施例中的步骤对应，具体此处不再赘述。
[0263]
本技术还提供了一种山洪灾害风险动态预警装置，请参阅图5，图5为本技术提供的一种山洪灾害风险动态预警装置一个实施例，该山洪灾害风险动态预警装置包括：
[0264]
处理器501、存储器502、输入输出单元503、总线504；
[0265]
处理器501与存储器502、输入输出单元503以及总线504相连；
[0266]
存储器502保存有程序，处理器501调用程序以执行如上任一一种山洪灾害风险动态预警方法。
[0267]
本技术还涉及一种计算机可读存储介质，计算机可读存储介质上保存有程序，当程序在计算机上运行时，使得计算机执行如上任一一种山洪灾害风险动态预警方法。
[0268]
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0269]
在本技术所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
[0270]
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0271]
另外，在本技术各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0272]
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本技术各个实施例方法的全部或部分
步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。