构造湖漫溢溃坝隐患预警方法、系统、介质及设备与流程

1.本发明属于灾害预警技术领域，尤其涉及构造湖漫溢溃坝隐患预警方法、系统、介质及设备。


背景技术：

2.构造湖漫溢溃坝灾害是随全球气候变化而在青藏高原出现的新型灾种。近年来，在全球变暖的背景下，青藏高原极端气候频发、冰雪消融加剧、湖泊扩显著，卓乃湖、阿克苏库勒湖、贝勒克勒克湖等先后发生漫溢溃坝，严重威胁下游生态环境、人民生命财产和重要基础设施安全，亟需引起关注。目前对青藏高原湖泊灾害隐患的研究多集中于冰川湖，而构造湖漫溢溃坝作为一类全新灾害类型，还未引起足够重视。
3.高原冰川湖与构造湖发生溃决的成因类型具有本质差异，冰川湖多存在溢流口，常规情况下较为稳定，其溃决多是由冰滑坡或冰崩入湖形成的瞬时涌浪击穿冰碛坝体所导致的结果；构造湖漫溢溃坝是指具备天然阻水坝体和溢流通道的闭流湖，上涨水量超出临界库容导致湖水外溢并沿溢流通道形成冲沟，经由溯源侵蚀进一步破坏天然阻水坝体形成溃决洪水，造成下游基础设施毁坏或生态损毁的现象。因此，以往相关学者针对冰川湖溃决隐患提出的成灾条件、主要溃决特征、影响因素分析、易发性评价等方法并不适用于构造湖(地理学报,1989(03):343-351+385-352；冰川冻土,1996(04):61-66；山地学报,2003(s1):128-132；灾害学,2008(01):55-60；自然资源学报,2014,29(08):1377-1390；国土资源遥感,2016,28(03):110-115)；而针对目前已溃决构造湖的隐患分析还停留在溃堤成因、溃决过程及机理研究等方面(地理学报,2012,67(05):689-698；冰川冻土,2017,39(05):949-956；冰川冻土,2018,40(01):47-54；冰川冻土,2019,41(06):1467-1474；冰川冻土,2020,42(04):1344-1352.)，未见有漫溢溃坝隐患的全流程早期预警方法提出。随着遥感技术的快速发展及对高原构造湖漫溢溃坝成灾条件认识的加深，为解决上述问题提供了新的思路。


技术实现要素：

4.本发明所要解决的技术问题是提供构造湖漫溢溃坝隐患预警方法、系统、介质及设备。
5.本发明解决上述技术问题的技术方案如下：一种高原构造湖漫溢溃坝隐患预警方法，包括：
6.步骤1，获取构造湖遥感影像；
7.步骤2，对所述构造湖遥感影像进行分析处理，筛选出具有成灾条件的目标构造湖；
8.步骤3，确定所述目标构造湖所属气候区的类型；
9.步骤4，基于所述类型确定该目标构造湖的模拟年内预测水量增量；
10.步骤5，基于所述模拟年内预测水量增量确定所述目标构造湖的漫溢溃坝概率，根
据所述漫溢溃坝概率生成预警信息。
11.本发明的有益效果是：高原构造湖漫溢溃坝作为一种新型灾害类型，由于青藏高原湖泊众多、湖区环境地质条件多样、时空气候变化情况复杂，目前暂无有效早期识别监测预警方法。本专利提出一种高原构造湖漫溢溃坝隐患的早期预警方法，有益效果包含以下四个方面：一是提出基于构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征的漫溢溃坝条件判断方法，可准确筛选隐患分析目标对象；二是综合考虑引发构造湖扩张的流域蒸散发及补给条件，划分气候区类型，可为准确开展构造湖漫溢渍坝隐患评估提供有价值的分析依据；三是全面考虑不同时间尺度下流域降水特征的延续性和突变性，提出多时一空气候特征的年内水量增量模拟方法，可有效提升预测精度；四是提出判断构造湖年内发生漫溢渍坝风险等级的概率区间，避免主观因素干预，从而具有普遍适用性，适合在青藏高原地区湖泊分布区的大范围应用。
12.在上述技术方案的基础上，本发明还可以做如下改进。
13.进一步，所述步骤1具体为：
14.获取高原上面积在预设值以上的所有湖泊的遥感图像，并对所述所有湖泊的遥感影像进行筛选，提出冰川湖以及堰塞湖，得出构造湖遥感影像。
15.进一步，所述步骤2具体为：
16.通过漫溢渍坝条件判断公式，将μ不小于预设指数的构造湖判定为具有成灾条件的目标构造湖，所述漫溢渍坝条件判断公式具体为：
17.μ＝0.4798
·
w1+0.1857
·
w2+0.280θ
·
w3+0.0536
·
w418.其中，μ为任意一个构造湖漫溢渍坝条件指数，w1为任意一个构造湖连通性指数，w2为任意一个构造湖坝体特征指数，w3为任意一个构造湖外泄通道指数，w4为任意一个构造湖形态指数；μ∈[0，1]，形态指数；μ∈[0，1]，q
in
为任意一个构造湖入湖水量，q
out
为任意一个构造湖排泄水量，k
vn
为任意一个构造湖构造湖天然坝体新鲜岩石压缩波速度，k
vw
为任意一个构造湖风化岩石压缩波速度，r为任意一个构造湖构造湖最小外接圆半径，从湖岸处以r/4为步长依次计算河道比降，ρi为湖岸外第i处河道比降，c
l
为任意一个构造湖构造湖岸线周长，s
l
为任意一个构造湖构造湖水域面积。
[0019]
进一步，所述步骤3具体为：
[0020]
根据类型判断公式，确定所述目标构造湖所属气候区的类型，所述类型判断公式为：
[0021][0022]
其中，q为气候区类型，a
t
为给定时间尺度内的气温变化量a
p
为同时间尺度流域降水量变化量，δsg为同时间尺度冰川面积变化量，q1，q2，q3为第一类型、第二类型以及第三类型。
[0023]
进一步，所述步骤4具体为：
[0024]
通过增量公式，确定所述类型对应的模拟年内预测水量增量，所述增量公式具体为：
[0025][0026]
其中，δv
fi
为目标构造湖年内预测水量增量，θy为年时间尺度常规降水概率，θm为月时间尺度极端降水概率，θd为日时间尺度极端降水概率，δv
yr
为年常规降水条件下的湖泊水量增量，δv
mr
为月度极端降水条件下的湖泊水量增量，δv
dr
为日强降水条件下的湖泊水量增量，δv
dr
＝sb×
p
max
×
λ
×
k+s
t
×
p
max
×
λ，β为月度汛期极端降水系数，k为流域产流系数，为流域月均降水量，为流域内冰川月均变化量，et0为气象站潜在蒸散量，p
max
为历史最大单日降水量，λ为历史最大单日降水时长。
[0027]
进一步，所述步骤5具体为：
[0028]
通过概率公式，确定所述目标构造湖的漫溢渍坝概率，根据所述漫溢渍坝概率生成预警信息，所述概率公式具体为：
[0029][0030]
其中，p
l
为构造湖年内漫溢渍坝概率，为任意一个目标构造湖外围最大古湖岸线平均高程，h
l
为任意一个目标构造湖湖面高程，s
old
为任意一个目标构造湖最大古湖岸线面积。
[0031]
本发明解决上述技术问题的另一种技术方案如下：一种高原构造湖漫溢溃坝隐患预警系统，包括：
[0032]
获取模块，用于获取构造湖遥感影像；
[0033]
筛选模块，用于对所述构造湖遥感影像进行分析处理，筛选出具有成灾条件的目标构造湖；
[0034]
确定模块，用于确定所述目标构造湖所属气候区的类型；
[0035]
预测模块，用于基于所述类型确定该目标构造湖的模拟年内预测水量增量；
[0036]
预警模块，用于基于所述模拟年内预测水量增量确定所述目标构造湖的漫溢溃坝概率，根据所述漫溢溃坝概率生成预警信息。
[0037]
本发明的有益效果是：高原构造湖漫溢溃坝作为一种新型灾害类型，由于青藏高原湖泊众多、湖区环境地质条件多样、时空气候变化情况复杂，目前暂无有效早期识别监测预警方法。本专利提出一种高原构造湖漫溢溃坝隐患的早期预警方法，有益效果包含以下四个方面：一是提出基于构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征的漫溢溃坝条件判断方法，可准确筛选隐患分析目标对象；二是综合考虑引发构造湖扩张的流域蒸散发及补给条件，划分气候区类型，可为准确开展构造湖漫溢溃坝隐患评估提供有价值的分析依据；三是全面考虑不同时间尺度下流域降水特征的延续性和突变性，提出多时—空气候特征的年内水量增量模拟方法，可有效提升预测精度；四是提出判断构造湖年内发生漫溢溃坝风险等级的概率区间，避免主观因素干预，从而具有普遍适用性，适合在青藏高原地区湖泊分布区的大范围应用。
[0038]
进一步，所述获取模块具体用于：
[0039]
获取高原上面积在预设值以上的所有湖泊的遥感图像，并对所述所有湖泊的遥感
影像进行筛选，提出冰川湖以及堰塞湖，得出构造湖遥感影像。
[0040]
进一步，所述筛选模块具体用于：
[0041]
通过漫溢渍坝条件判断公式，将μ不小于预设指数的构造湖判定为具有成灾条件的目标构造湖，所述漫溢渍坝条件判断公式具体为：
[0042]
μ＝0.4798
·
w1+0.1857
·
w2+0.2809
·
w3+0.0536
·
w4[0043]
其中，μ为任意一个构造湖漫溢渍坝条件指数，w1为任意一个构造湖连通性指数，w2为任意一个构造湖坝体特征指数，w3为任意一个构造湖外泄通道指数，w4为任意一个构造湖形态指数；μ∈[0，1]，形态指数；μ∈[0，1]，q
in
为任意一个构造湖入湖水量，q
out
为任意一个构造湖排泄水量，k
vn
为任意一个构造湖构造湖天然坝体新鲜岩石压缩波速度，k
vw
为任意一个构造湖风化岩石压缩波速度，r为任意一个构造湖构造湖最小外接圆半径，从湖岸处以r/4为步长依次计算河道比降，ρi为湖岸外第i处河道比降，c
l
为任意一个构造湖构造湖岸线周长，s
l
为任意一个构造湖构造湖水域面积。
[0044]
进一步，所述确定模块具体用于：
[0045]
根据类型判断公式，确定所述目标构造湖所属气候区的类型，所述类型判断公式为：
[0046][0047]
其中，q为气候区类型，a
t
为给定时间尺度内的气温变化量a
p
为同时间尺度流域降水量变化量，δsg为同时间尺度冰川面积变化量，q1，q2，q3为第一类型、第二类型以及第三类型。
[0048]
进一步，所述预测模块具体用于：
[0049]
通过增量公式，确定所述类型对应的模拟年内预测水量增量，所述增量公式具体为：
[0050][0051]
其中，δv
fi
为目标构造湖年内预测水量增量，θy为年时间尺度常规降水概率，θm为月时间尺度极端降水概率，θd为日时间尺度极端降水概率，δv
yr
为年常规降水条件下的湖泊水量增量，δv
mr
为月度极端降水条件下的湖泊水量增量，δv
dr
为日强降水条件下的湖泊水量增量，δv
dr
＝sb×
p
max
×
λ
×
k+s
t
×
p
max
×
λ，β为月度汛期极端降水系数，k为流域产流系数，为流域月均降水量，为流域内冰川月均变化量，et0为气象站潜在蒸散量，p
max
为历史最大单日降水量，λ为历史最大单日降水时长。
[0052]
进一步，所述预警模块具体用于：
[0053]
通过概率公式，确定所述目标构造湖的漫溢渍坝概率，根据所述漫溢渍坝概率生成预警信息，所述概率公式具体为：
[0054]
[0055]
其中，p
l
为构造湖年内漫溢渍坝概率，为任意一个目标构造湖外围最大古湖岸线平均高程，h
l
为任意一个目标构造湖湖面高程，s
old
为任意一个目标构造湖最大古湖岸线面积。
[0056]
本发明解决上述技术问题的另一种技术方案如下：一种存储介质，所述存储介质中存储有指令，当计算机读取所述指令时，使所述计算机执行如上述任一项所述的一种高原构造湖漫溢渍坝隐患预警方法。
[0057]
本发明的有益效果是：高原构造湖漫溢渍坝作为一种新型灾害类型，由于青藏高原湖泊众多、湖区环境地质条件多样、时空气候变化情况复杂，目前暂无有效早期识别监测预警方法。本专利提出一种高原构造湖漫溢渍坝隐患的早期预警方法，有益效果包含以下四个方面：一是提出基于构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征的漫溢溃坝条件判断方法，可准确筛选隐患分析目标对象；二是综合考虑引发构造湖扩张的流域蒸散发及补给条件，划分气候区类型，可为准确开展构造湖漫溢溃坝隐患评估提供有价值的分析依据；三是全面考虑不同时间尺度下流域降水特征的延续性和突变性，提出多时—空气候特征的年内水量增量模拟方法，可有效提升预测精度；四是提出判断构造湖年内发生漫溢溃坝风险等级的概率区间，避免主观因素干预，从而具有普遍适用性，适合在青藏高原地区湖泊分布区的大范围应用。
[0058]
本发明解决上述技术问题的另一种技术方案如下：一种电子设备，包括上述存储介质、执行上述存储介质内的指令的处理器。
[0059]
本发明的有益效果是：高原构造湖漫溢溃坝作为一种新型灾害类型，由于青藏高原湖泊众多、湖区环境地质条件多样、时空气候变化情况复杂，目前暂无有效早期识别监测预警方法。本专利提出一种高原构造湖漫溢溃坝隐患的早期预警方法，有益效果包含以下四个方面：一是提出基于构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征的漫溢溃坝条件判断方法，可准确筛选隐患分析目标对象；二是综合考虑引发构造湖扩张的流域蒸散发及补给条件，划分气候区类型，可为准确开展构造湖漫溢溃坝隐患评估提供有价值的分析依据；三是全面考虑不同时间尺度下流域降水特征的延续性和突变性，提出多时—空气候特征的年内水量增量模拟方法，可有效提升预测精度；四是提出判断构造湖年内发生漫溢溃坝风险等级的概率区间，避免主观因素干预，从而具有普遍适用性，适合在青藏高原地区湖泊分布区的大范围应用。
附图说明
[0060]
图1为本发明一种高原构造湖漫溢溃坝隐患预警方法的实施例提供的流程示意图；
[0061]
图2为本发明一种高原构造湖漫溢溃坝隐患预警系统的实施例提供的结构框架图。
具体实施方式
[0062]
以下对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
[0063]
如图1所示，一种高原构造湖漫溢溃坝隐患预警方法，包括：
[0064]
步骤1，获取构造湖遥感影像；
[0065]
步骤2，对所述构造湖遥感影像进行分析处理，筛选出具有成灾条件的目标构造湖；
[0066]
步骤3，确定所述目标构造湖所属气候区的类型；
[0067]
步骤4，基于所述类型确定该目标构造湖的模拟年内预测水量增量；
[0068]
步骤5，基于所述模拟年内预测水量增量确定所述目标构造湖的漫溢溃坝概率，根据所述漫溢溃坝概率生成预警信息。
[0069]
在一些可能的实施方式中，高原构造湖漫溢溃坝作为一种新型灾害类型，由于青藏高原湖泊众多、湖区环境地质条件多样、时空气候变化情况复杂，目前暂无有效早期识别监测预警方法。本专利提出一种高原构造湖漫溢溃坝隐患的早期预警方法，有益效果包含以下四个方面：一是提出基于构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征的漫溢溃坝条件判断方法，可准确筛选隐患分析目标对象；二是综合考虑引发构造湖扩张的流域蒸散发及补给条件，划分气候区类型，可为准确开展构造湖漫溢溃坝隐患评估提供有价值的分析依据；三是全面考虑不同时间尺度下流域降水特征的延续性和突变性，提出多时—空气候特征的年内水量增量模拟方法，可有效提升预测精度；四是提出判断构造湖年内发生漫溢渍坝风险等级的概率区间，避免主观因素干预，从而具有普遍适用性，适合在青藏高原地区湖泊分布区的大范围应用。
[0070]
需要说明的是，步骤1，获取构造湖遥感影像。针对该步骤的具体解释参考如下示例：
[0071]
利用高分辨率遥感数据获取青藏高原面积在1km2以上的湖泊，根据成因类型筛选出构造湖，剔除冰川湖、堰塞湖等其他类型湖泊。
[0072]
步骤2，对所述构造湖遥感影像进行分析处理，筛选出具有成灾条件的目标构造湖。针对该步骤的具体解释参考如下示例：
[0073]
根据对已漫溢渍决构造湖成因分析，获取构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征，综合构建成灾决定因素的漫溢渍坝条件判断公式，筛选出具成灾条件的目标构造湖。
[0074]
μ＝0.4798
·
w1+0.1857
·
w2+0.2809
·
w3+0.0536
·
w4[0075]
其中，μ为任意一个构造湖漫溢渍坝条件指数，w1为任意一个构造湖连通性指数，w2为任意一个构造湖坝体特征指数，w3为任意一个构造湖外泄通道指数，w4为任意一个构造湖形态指数；μ∈[0，1]，无量纲；形态指数；μ∈[0，1]，无量纲；以上各指数均在0到1范围内，无量纲；指数前系数为权重值，由yaahp方法确定；当μ≥0.75时，认为满足漫溢渍坝成灾条件，q
in
为任意一个构造湖入湖水量，q
out
为任意一个构造湖排泄水量，k
vn
为任意一个构造湖构造湖天然坝体新鲜岩石压缩波速度，k
vw
为任意一个构造湖风化岩石压缩波速度，r为任意一个构造湖构造湖最小外接圆半径，从湖岸处以r/4为步长依次计算河道比降，ρi为湖岸外第i处河道比降，c
l
为任意一个构造湖构造湖岸线周长，s
l
为任意一个构造湖构造湖水域面积。
[0076]
涉及漫溢渍坝条件判断公式的构建，具体过程如下：
[0077]
结合已渍决湖泊分析，高原构造湖发生漫溢渍坝的先决条件包括四个方面，即属
于闭流型湖泊，具备天然阻水坝体和存在外泄通道，并且在形态结构上较为特殊。根据上述要素之间的相关性，建立漫溢渍坝条件指数的一般表达式：
[0078]
式中，μ＝α
·
w1+β
·
w2+γ
·
w3+δ
·
w4；
[0079]
式中，μ为构造湖漫溢渍坝条件指数，w1为连通性指数，w2为坝体特征指数，w3为外泄通道指数，w4为形态指数，上述各指数均∈[0，1]，无量纲；α，β，γ，δ系为经验公式未知数，无量纲。
[0080]
为有效表征湖泊连通性特征对漫溢渍坝的影响，与常规判定湖泊连通性的方法不同，本发明创新型的引入连通性指数概念，通过构造湖排泄水量与入湖水量特征对比分析，实现其连通特征的量化。
[0081][0082]
式中，q
in
为入湖水量，q
out
为排泄水量，分式内单位一致以消除量纲。当构造湖连通性指数趋近于0时，说明其为吞吐湖，不具备漫溢渍坝条件；当指数趋近于1时，说明其为闭流湖，存在蓄水后漫溢渍坝条件。
[0083]
为有效表征坝体物质结构特征对漫溢渍坝的影响，本发明选用压缩波速度比这一指标，来量化坝体结构类型，以区分易渍坝的松散土质坝体及较难渍决的岩质坝体。
[0084][0085]
式中，k
vn
为构造湖天然坝体新鲜岩石压缩波速度，k
vw
为风化岩石压缩波速度，其比值表示坝体物质结构特征，分式内单位一致以消除量纲。当坝体物质结构指数趋近于1时，其结构接近于强风化松散土质坝体，当趋近于0时，其结构接近于未风化稳定岩质坝体。
[0086]
为量化外泄通道特征对漫溢溃坝的影响，本发明选用河道平均比降这一指标，对不同类型的外泄通道进行量化描述。
[0087][0088]
式中，r为构造湖最小外接圆半径，从湖岸处以r/4为步长依次计算河道比降，ρi为湖岸外第i处河道比降，分式内单位一致以消除量纲。当外泄通道特征指数趋近于0时，其湖岸外侧均为高地，不存在外泄条件；当指数趋近于0.5时，其存在漫溢外流条件；当指数趋近于1时，其湖岸外侧存在较陡外泄通道，具备引发溃坝条件。
[0089]
为量化构造湖形态特征及湖岸线发育复杂程度对漫溢溃坝的影响，本发明引入构造湖形态指数这一概念，通过对已出现漫溢溃坝的构造湖岸线的几何形态分析构建计算公式，以表征其对漫溢溃坝的影响程度。
[0090][0091]
式中，c
l
为构造湖岸线周长，s
l
为构造湖水域面积，其比值表示构造湖形态的复杂程度，分式内单位一致以消除量纲。当构造湖形态指数趋近于1时，其形态相似性越强，岸线形态越规则，易出现溢流口，已发生漫溢溃坝湖泊均具备这一特征；当趋近于0时，其相似性越差，岸线形态越复杂，一般不存在溢流口，位于流域内终点的尾闾湖多具备此特征。
[0092]
为确定上述4个指数因子对判定构造湖漫溢溃坝条件的影响权重，本发明建立分层交错评价指标的层次结构，利用yaahp软件求解判断矩阵特征向量(表1)，求解α，β，γ，δ在μ的经验公式中重要性权重。
[0093]
表1指数因子的配对比较矩阵
[0094][0095][0096]
指数因子的配对比较矩阵一致性cr＝0.0304＜0.1000，通过一致性检验。
[0097]
据此，求得连通性指数权重α＝0.4798，坝体特征指数β＝0.1857，外泄通道指数γ＝0.2809，形态指数δ＝0.0536。
[0098]
将前述指数公式及权重因子带入漫溢渍坝条件指数的一般表达式，得出构造湖漫溢渍坝条件指数μ的计算公式。
[0099][0100]
式中，μ∈[0，1]，无量纲；当μ≥0.75时，认为满足漫溢渍坝条件。
[0101]
步骤3，确定所述目标构造湖所属气候区的类型。针对该步骤的具体解释参考如下示例：
[0102]
根据气象站数据及高分辨率遥感数据，获取给定时间尺度下可能出现扩张的目标构造湖流域范围内的历史气温、降水量、冰川面积变化特征，划定目标构造湖所属气候区类型；
[0103][0104]
式中，q为气候区类型，a
t
为给定时间尺度内的气温变化量a
p
为同时间尺度流域降
水量变化量，δsg为同时间尺度冰川面积变化量。
[0105]
步骤4，基于所述类型确定该目标构造湖的模拟年内预测水量增量。针对该步骤的具体解释参考如下示例：
[0106]
根据气象站数据统计目标构造湖流域内年时间尺度常规降水概率、月时间尺度极端降水概率、日时间尺度极端降水概率及不同降水条件下的水量变化量，针对不同气候变化区分类模拟年内预测水量增量。
[0107][0108]
式中，δv
fi
为目标构造湖年内预测水量增量，θy为年时间尺度常规降水概率，θm为月时间尺度极端降水概率，θd为日时间尺度极端降水概率，δv
yr
为年常规降水条件下的湖泊水量增量，δv
mr
为月度极端降水条件下的湖泊水量增量，δv
dr
为日强降水条件下的湖泊水量增量。
[0109]
不同气候区内目标构造湖水量增量模拟，具体过程如下：
[0110]
年常规气候条件下，限定目标构造湖年内水量变化量将延续往年趋势，流域范围内的气候、水文地质等条件不会出现异常变化。据此，对区内满足成灾条件的30个构造湖的水域面积、流域面积及历史年均水量增量，构建面积流域比-历史年均水量增量关系曲线，进而模拟年常规气候条件下构造湖的年内水量增量：
[0111][0112]
式中，δv
ry
为目标构造湖年常规气候条件下的预测年内水量增量，s
l
为湖泊当前面积，sb为湖泊流域面积。
[0113]
月极端降水条件下，限定目标构造湖年内月度汛期最大降水不超过历史极值，综合考虑目标构造湖流域范围内的陆面降水、湖面降水、汛期降水系数、冰川融水及蒸发等条件，进而模拟月度极端降水条件下湖泊水量增量：
[0114][0115]
式中，δv
mr
为月度极端降水条件下的湖泊水量增量，β为月度汛期极端降水系数，k为流域产流系数，为流域月均降水量，为流域内冰川月均变化量，et0为气象站潜在蒸散量。
[0116]
日极端降水条件下，限定目标构造湖年内最大单次降水不超过历史极值，短时降水无需考虑蒸发损耗及冰川融水补给，据此结合陆面降水、湖面降水等条件，构建预测日常规气候条件下的水量增量模拟方法：
[0117]
δv
dr
＝sb×
p
max
×
λ
×
k+s
t
×
p
max
×
λ
[0118]
式中，δv
dr
为日强降水条件下的湖泊水量增量，p
max
为历史最大单日降水量，λ为历史最大单日降水时长。
[0119]
步骤5，基于所述模拟年内预测水量增量确定所述目标构造湖的漫溢渍坝概率，根据所述漫溢渍坝概率生成预警信息。针对该步骤的具体解释参考如下示例：
[0120]
基于上一步得到的目标构造湖年内预测水量增量，结合剩余库容进行综合判断年内发生漫溢渍坝概率，依照概率区间划分风险等级，进而实现隐患早期识别。
[0121][0122]
其中，p
l
为构造湖年内漫溢渍坝概率，为任意一个目标构造湖外围最大古湖岸线平均高程，h
l
为任意一个目标构造湖湖面高程，s
old
为任意一个目标构造湖最大古湖岸线面积。当p
l
∈(0.80，0.90]时，为低风险；当p
l
∈(0.90，0.95]时，为中风险；当p
l
＞0.95时，为高风险。
[0123]
优选地，在上述任意实施例中，所述步骤1具体为：
[0124]
获取高原上面积在预设值以上的所有湖泊的遥感图像，并对所述所有湖泊的遥感影像进行筛选，提出冰川湖以及堰塞湖，得出构造湖遥感影像。
[0125]
优选地，在上述任意实施例中，所述步骤2具体为：
[0126]
通过漫溢渍坝条件判断公式，将μ不小于预设指数的构造湖判定为具有成灾条件的目标构造湖，所述漫溢渍坝条件判断公式具体为：
[0127]
μ＝0.4798
·
w1+0.1857
·
w2+0.2809
·
w3+0.0536
·
w4[0128]
其中，μ为任意一个构造湖漫溢渍坝条件指数，w1为任意一个构造湖连通性指数，w2为任意一个构造湖坝体特征指数，w3为任意一个构造湖外泄通道指数，w4为任意一个构造湖形态指数；μ∈[0，1]，形态指数；μ∈[0，1]，q
in
为任意一个构造湖入湖水量，q
out
为任意一个构造湖排泄水量，k
vn
为任意一个构造湖构造湖天然坝体新鲜岩石压缩波速度，k
vw
为任意一个构造湖风化岩石压缩波速度，r为任意一个构造湖构造湖最小外接圆半径，从湖岸处以r/4为步长依次计算河道比降，ρi为湖岸外第i处河道比降，c
l
为任意一个构造湖构造湖岸线周长，s
l
为任意一个构造湖构造湖水域面积。
[0129]
优选地，在上述任意实施例中，所述步骤3具体为：
[0130]
根据类型判断公式，确定所述目标构造湖所属气候区的类型，所述类型判断公式为：
[0131][0132]
其中，q为气候区类型，a
t
为给定时间尺度内的气温变化量a
p
为同时间尺度流域降水量变化量，δsg为同时间尺度冰川面积变化量，q1，q2，q3为第一类型、第二类型以及第三类型。
[0133]
需要说明的是，q1为第一种气候区类型，对应暖湿化气候区，主要分布于羌塘高原、祁连山、昆仑山和冈底斯山等地，由于气温升高、降水增加及冰川融水补给，导致区内构造湖扩张比例达70％；q2为第二种气候区类型，对应暖干化气候区，主要分布于横断山、念青唐古拉山、唐古拉山和喜马拉雅山等地，由于气温升高、降水减少及冰川融水补给，区内构造湖扩张比例不足10％；q3为第三种气候区类型，对应冷湿化气候区，主要分布喀喇昆仑山、昆仑山及冈底斯山的部分区域，由于气温降低、降水增加、冰川融水补给，区内构造湖扩张比例仅为30％。
[0134]
优选地，在上述任意实施例中，所述步骤4具体为：
[0135]
通过增量公式，确定所述类型对应的模拟年内预测水量增量，所述增量公式具体为：
[0136][0137]
其中，δv
fi
为目标构造湖年内预测水量增量，θy为年时间尺度常规降水概率，θm为月时间尺度极端降水概率，θd为日时间尺度极端降水概率，δv
yr
为年常规降水条件下的湖泊水量增量，δv
mr
为月度极端降水条件下的湖泊水量增量，δv
dr
为日强降水条件下的湖泊水量增量，δv
dr
＝sb×
p
max
×
λ
×
k+s
t
×
p
max
×
λ，β为月度汛期极端降水系数，k为流域产流系数，为流域月均降水量，为流域内冰川月均变化量，et0为气象站潜在蒸散量，p
max
为历史最大单日降水量，λ为历史最大单日降水时长。
[0138]
优选地，在上述任意实施例中，所述步骤5具体为：
[0139]
通过概率公式，确定所述目标构造湖的漫溢渍坝概率，根据所述漫溢渍坝概率生成预警信息，所述概率公式具体为：
[0140][0141]
其中，p
l
为构造湖年内漫溢渍坝概率，为任意一个目标构造湖外围最大古湖岸线平均高程，h
l
为任意一个目标构造湖湖面高程，s
old
为任意一个目标构造湖最大古湖岸线面积。
[0142]
如图2所示，一种高原构造湖漫溢渍坝隐患预警系统，包括：
[0143]
获取模块100，用于获取构造湖遥感影像；
[0144]
筛选模块200，用于对所述构造湖遥感影像进行分析处理，筛选出具有成灾条件的目标构造湖；
[0145]
确定模块300，用于确定所述目标构造湖所属气候区的类型；
[0146]
预测模块400，用于基于所述类型确定该目标构造湖的模拟年内预测水量增量；
[0147]
预警模块500，用于基于所述模拟年内预测水量增量确定所述目标构造湖的漫溢溃坝概率，根据所述漫溢溃坝概率生成预警信息。
[0148]
在一些可能的实施方式中，高原构造湖漫溢溃坝作为一种新型灾害类型，由于青藏高原湖泊众多、湖区环境地质条件多样、时空气候变化情况复杂，目前暂无有效早期识别监测预警方法。本专利提出一种高原构造湖漫溢溃坝隐患的早期预警方法，有益效果包含以下四个方面：一是提出基于构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征的漫溢溃坝条件判断方法，可准确筛选隐患分析目标对象；二是综合考虑引发构造湖扩张的流域蒸散发及补给条件，划分气候区类型，可为准确开展构造湖漫溢溃坝隐患评估提供有价值的分析依据；三是全面考虑不同时间尺度下流域降水特征的延续性和突变性，提出多时—空气候特征的年内水量增量模拟方法，可有效提升预测精度；四是提出判断构造湖年内发生漫溢溃坝风险等级的概率区间，避免主观因素干预，从而具有普遍适用性，适合在青藏高原地区湖泊分布区的大范围应用。
[0149]
优选地，在上述任意实施例中，所述获取模块100具体用于：
[0150]
获取高原上面积在预设值以上的所有湖泊的遥感图像，并对所述所有湖泊的遥感影像进行筛选，提出冰川湖以及堰塞湖，得出构造湖遥感影像。
[0151]
优选地，在上述任意实施例中，所述筛选模块200具体用于：
[0152]
通过漫溢溃坝条件判断公式，将μ不小于预设指数的构造湖判定为具有成灾条件的目标构造湖，所述漫溢渍坝条件判断公式具体为：
[0153]
μ＝0.4798
·
w1+0.1857
·
w2+0.2809
·
w3+0.0536
·
w4[0154]
其中，μ为任意一个构造湖漫溢渍坝条件指数，w1为任意一个构造湖连通性指数，w2为任意一个构造湖坝体特征指数，w3为任意一个构造湖外泄通道指数，w4为任意一个构造湖形态指数；μ∈[0，1]，形态指数；μ∈[0，1]，q
in
为任意一个构造湖入湖水量，q
out
为任意一个构造湖排泄水量，k
vn
为任意一个构造湖构造湖天然坝体新鲜岩石压缩波速度，k
vw
为任意一个构造湖风化岩石压缩波速度，r为任意一个构造湖构造湖最小外接圆半径，从湖岸处以r/4为步长依次计算河道比降，ρi为湖岸外第i处河道比降，c
l
为任意一个构造湖构造湖岸线周长，s
l
为任意一个构造湖构造湖水域面积。
[0155]
优选地，在上述任意实施例中，所述确定模块300具体用于：
[0156]
根据类型判断公式，确定所述目标构造湖所属气候区的类型，所述类型判断公式为：
[0157][0158]
其中，q为气候区类型，a
t
为给定时间尺度内的气温变化量a
p
为同时间尺度流域降水量变化量，δsg为同时间尺度冰川面积变化量，q1，q2，q3为第一类型、第二类型以及第三类型。
[0159]
优选地，在上述任意实施例中，所述预测模块400具体用于：
[0160]
通过增量公式，确定所述类型对应的模拟年内预测水量增量，所述增量公式具体为：
[0161][0162]
其中，δv
fi
为目标构造湖年内预测水量增量，θy为年时间尺度常规降水概率，θm为月时间尺度极端降水概率，θd为日时间尺度极端降水概率，δv
yr
为年常规降水条件下的湖泊水量增量，δv
mr
为月度极端降水条件下的湖泊水量增量，δv
dr
为日强降水条件下的湖泊水量增量，δv
dr
＝sb×
p
max
×
λ
×
k+s
t
×
p
max
×
λ，β为月度汛期极端降水系数，k为流域产流系数，为流域月均降水量，为流域内冰川月均变化量，et0为气象站潜在蒸散量，p
max
为历史最大单日降水量，λ为历史最大单日降水时长。
[0163]
优选地，在上述任意实施例中，所述预警模块500具体用于：
[0164]
通过概率公式，确定所述目标构造湖的漫溢渍坝概率，根据所述漫溢渍坝概率生成预警信息.所述概率公式具体为：
[0165][0166]
其中，p
l
为构造湖年内漫溢渍坝概率，为任意一个目标构造湖外围最大古湖岸线平均高程，h
l
为任意一个目标构造湖湖面高程，s
old
为任意一个目标构造湖最大古湖岸线面积。
[0167]
本发明解决上述技术问题的另一种技术方案如下：一种存储介质，所述存储介质中存储有指令，当计算机读取所述指令时，使所述计算机执行如上述任一项所述的一种高原构造湖漫溢渍坝隐患预警方法。
[0168]
在一些可能的实施方式中，高原构造湖漫溢渍坝作为一种新型灾害类型，由于青藏高原湖泊众多、湖区环境地质条件多样、时空气候变化情况复杂，目前暂无有效早期识别监测预警方法。本专利提出一种高原构造湖漫溢渍坝隐患的早期预警方法，有益效果包含以下四个方面：一是提出基于构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征的漫溢渍坝条件判断方法，可准确筛选隐患分析目标对象；二是综合考虑引发构造湖扩张的流域蒸散发及补给条件，划分气候区类型，可为准确开展构造湖漫溢渍坝隐患评估提供有价值的分析依据；三是全面考虑不同时间尺度下流域降水特征的延续性和突变性，提出多时—空气候特征的年内水量增量模拟方法，可有效提升预测精度；四是提出判断构造湖年内发生漫溢溃坝风险等级的概率区间，避免主观因素干预，从而具有普遍适用性，适合在青藏高原地区湖泊分布区的大范围应用。
[0169]
本发明解决上述技术问题的另一种技术方案如下：一种电子设备，包括上述存储介质、执行上述存储介质内的指令的处理器。
[0170]
在一些可能的实施方式中，高原构造湖漫溢溃坝作为一种新型灾害类型，由于青藏高原湖泊众多、湖区环境地质条件多样、时空气候变化情况复杂，目前暂无有效早期识别监测预警方法。本专利提出一种高原构造湖漫溢溃坝隐患的早期预警方法，有益效果包含以下四个方面：一是提出基于构造湖连通性、坝体特征、外泄通道特征、湖泊形态特征的漫溢溃坝条件判断方法，可准确筛选隐患分析目标对象；二是综合考虑引发构造湖扩张的流域蒸散发及补给条件，划分气候区类型，可为准确开展构造湖漫溢溃坝隐患评估提供有价值的分析依据；三是全面考虑不同时间尺度下流域降水特征的延续性和突变性，提出多时—空气候特征的年内水量增量模拟方法，可有效提升预测精度；四是提出判断构造湖年内发生漫溢溃坝风险等级的概率区间，避免主观因素干预，从而具有普遍适用性，适合在青藏高原地区湖泊分布区的大范围应用。
[0171]
读者应理解，在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0172]
在本技术所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的方法实施例仅仅是示意性的，例如，步骤的划分，仅仅为
一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个步骤可以结合或者可以集成到另一个步骤，或一些特征可以忽略，或不执行。
[0173]
上述方法如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0174]
以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。