一种基于水热组合的冰川泥石流区域预警方法

本发明涉及冰川泥石流领域，具体涉及一种基于水热组合的冰川泥石流区域预警方法。

背景技术：

随着经济建设和社会发展需要，人们逐步在高海拔冰川山地地区进行工程建设，在这工程建设过程中，不可避免的会遇到冰川泥石流灾害，因此有必要对冰川泥石流发生建立预警机制。

相比雨季的暴雨泥石流，冰川泥石流本身属于更小概率的事件，其起动机理研究仍处于起步阶段，再加上目前冰川泥石流数据存在样本量小、样本噪声大的问题，现有冰川泥石流预警则更多是停留在局地经验取值的层面或群测群防(即监测员巡查)，因此，无论是实际应用还是工程设计规范推荐的泥石流相关计算方法都是以简单经验性的公式计算，缺乏区域普适性应用价值。复杂的方法(如机器学习)在这个方向上，由于多数研究仅为针对单个冰川或少数次泥石流事件的分析，无法建立区域普适性的冰川泥石流有效预警方法。现阶段急需一种应对冰川泥石流区域预警的普适性有效预警方法。

技术实现要素：

为解决上述问题，本发明目的在于提供一种基于水热组合的冰川泥石流区域预警方法，该方法解决了现有技术一直未能解决的冰川泥石流地区定量预警的问题，提供了一种能随着泥石流数据的增多持续有效性提升预警结果准确性的冰川泥石流区域预警曲线和预警方法，为今后对冰川泥石流区域预警提供了量化指导和理论基础，具有普适意义。

本发明通过下述技术方案实现：

一种基于水热组合的冰川泥石流区域预警方法，其特征在于，包括以下几个步骤：s1、在冰川泥石流预警建模研究区域内采集冰川泥石流相关数据；s2、构建预警阈值曲线模型：ss1、基于s1中采集的冰川泥石流相关数据，选取有效累积雨量re和有效正积温ate作为泥石流预警的关键气象因子，然后对冰川泥石流预警建模研究区域内的每个地理位置每次泥石流发生前一个月内的关键气象因子进行计算，并对计算结果进行泥石流事件和非泥石流事件的标签分类；ss2、构造泥石流事件关键气象因子的概念函数：ate＝αre^β,----------公式(1)，其中α和β为待定参数；ss3、获得参数α和β，形成冰川泥石流区域预警阈值曲线：利用随机函数，生成均匀分布的随机初始种子，作为α和β待定参数的初值，设定拟提出的泥石流事件出现的概率值，再将标签分类后的泥石流事件和非泥石流事件的有效正积温和有效累积降雨数据代入其中，通过至少50000次以上的循环迭代，获得满足条件的最优α和β的参数，最终形成冰川泥石流区域预警阈值曲线；ss4、根据ss3构建的预警阈值曲线，结合需预警区域监测或预报的气象条件是否落在预警阈值曲线上的位置关系，进行预警：若气象条件落在曲线上方区域内，则按照曲线以上的泥石流事件出现统计概率进行相应预警。

其中，标签分类的分类标准是：有历史泥石流记录对应的气象条件被标记为泥石流事件样本，每次泥石流暴发前随机抽取的多次非泥石流事件对应的气象条件被标记为非泥石流事件样本。

所述冰川泥石流预警建模研究区域的选择原则包括：(1)、冰川泥石流所处的研究区内多数流域上游物源区现在处于冰川环境，(2)、气象观测数据或高分辨率气象再分析资料覆盖的时间段内区域发生泥石流次数不少于20次。

s2中有效累积降雨量re表示为：

其中re是有效降水量；右侧第一项为泥石流发生前n天的有效前期降水量，ri为第i天的累积降水量；r0为泥石流发生当天的日降水量；k是前期降雨作用的衰减因子，表明前期降雨对泥石流形成的影响大小，k和n可根据经验设置或实验测定；具体的，k根据邻近区域的已有研究取值，或当地开展水力学和土力学模型实验取值。

s2中冰雪区的有效正积温ate表示为：

右侧第一项为表示冰川泥石流预警建模研究区域内的每个地理位置每次泥石流发生前n天大于0摄氏度的有效积温部分，ti为第i天的正积温，而泥石流发生当天的正积温，通过天气预报获取。在本发明中，有效正积温均为预警地区的冰雪区或者冰川区域的有效累积温度。

进一步的，ss3中随机函数为python的random.uniform函数，参数α和β在既定预警概率的约束条件下，通过蒙特卡洛采样与筛选方法进行确定。本文的随机函数不局限于上述函数，实施例中以python的random.uniform函数为例进行具体的计算说明。

本发明将有效正积温和有效累积降雨作为关键气象因子，基于关键气象因子构建预警函数，通过对预警函数的创新构建和以可靠性概率为约束条件的结合，构建了能够将泥石流和非泥石流二者的气象特征分离开的分级阈值预警曲线，一方面基于对高易发泥石流气象条件的提取，解决了目前冰川泥石流样本量小、样本噪声大的问题，另一方面降低了预警误判的可能性。在本发明中，预警函数的创新构建和以可靠性概率为约束条件的结合主要是基于“舍弃对极小概率事件的预警，将重心放在相对高概率的事件”的大胆构思，具体可理解为：舍弃较低概率的有效累积雨量和有效正积温组合条件下极易与非泥石流气象条件混淆的小概率事件，而注重较高概率的有效累积雨量和有效正积温组合条件下那些不易与非泥石流气象条混淆件的大概率事件。

实施预警时，α和β确定进行拟合曲线后，当预警地的气象条件落到拟合的曲线以上区域内，则按照曲线以上区域的泥石流事件出现统计概率进行相应预警。倘若拟合了不同概率下的曲线模型，则可实现分级预警。如同时拟合50％和80％概率的两条曲线模型，则根据预警区域内的实际气象条件落在这两条曲线划分的几个区域位置关系进行分级预警。

本发明与现有技术相比，具有如下的优点和有益效果：

1、本发明能够根据泥石流样本的可获得性，自定义寻找不同概率的冰川泥石流气象预警曲线，为当地气象、地灾部门的分级预警提供参考；

2、本发明能够突破以往冰川泥石流预警时简单经验式公式预测的局限性，具有普遍应用价值。

附图说明

此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：

图1为实施例中不同气象组合条件下的泥石流概率发生分析图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。

实施例1

冰川泥石流预警流程具体实施如下：

1、选择冰川泥石流频发的藏东南帕隆藏布中下游区域开展实施例演示。研究区受印度洋夏季风影响，沿海拔高度有较大的降水梯度。该流域的年降雨量随着海拔高度的增加而减少，波密气象站的降雨量约为900mm。由于强烈的夏季风和孟加拉湾潮湿的气候，帕隆藏布成为高原上气候最湿润的地区。由于高山环境的降温作用，形成了中国最大的海洋冰川分布区。大量冰川泥石流严重威胁着川藏国道、川藏铁路及少数民族村镇。

2、以往由于当地许多泥石流未造成较严重社会经济损失，因此没有得到及时记录，泥石流发生的准确信息很难获得。在这个实例中，收集了24个灾害性冰川泥石流时间，包括地点和精确到日的日期。收集了高分辨率网格化气温和降水数据集，记录跨度为1979-2018年，时空分辨率为3小时和0.1°。此外，结合中国第二次冰川资源清查数据集和遥感解译，提取了包括面积、边界和位置在内的冰川信息，最终按照上述的泥石流预警建模研究区域选择原则，选择如表1中展示的地理位置作为泥石流预警建模研究区域。

表1泥石流资料

3、在上述泥石流预警建模研究区域内，收集每个地理位置处每次泥石流发生前一个月的冰川泥石流降水量和温度数据，选取有效累积雨量re和有效正积温ate作为泥石流预警的关键气象因子，用计算有效累积雨量re，其中re是有效降水量；右侧第一项为泥石流发生前n天的有效前期降水量，ri为第i天的累积降水量；r0为泥石流发生当天的日降水量；k是前期降雨作用的衰减因子，表明前期降雨对泥石流形成的影响大小，k和n可根据经验设置或实验测定；用计算有效正积温ate，右侧第一项为表示冰川泥石流预警建模研究区域内的每个地理位置每次泥石流发生前n天大于0摄氏度的有效积温部分，ti为第i天的正积温，而泥石流发生当天的正积温，通过天气预报获取，结果如表2。

表2泥石流事件7日有效雨量和3个月有效累积温度

4、对上述计算得到的关键气象因子：有效累积雨量re和有效正积温ate进行泥石流事件和非泥石流事件的标签分类。标签分类的分类标准是：有历史泥石流记录对应的气象条件被标记为泥石流事件样本，每次泥石流暴发前随机抽取的多次非泥石流事件对应的气象条件被标记为非泥石流事件样本。

5、构造泥石流事件关键气象因子的概念函数的有效正积温与有效累积雨量的概念函数：

ate＝αre^β,如图1所示，设定拟提出的预警对应的所有事件中泥石流出现的概率值为80％，用random.uniform函数，生成均匀分布的随机初始种子，作为α和β为待定参数的初值，再将所有记录的泥石流事件和非泥石流事件的有效正积温和有效累积降雨数据代入其中，通过50000次的循环迭代，获得满足条件的最优α＝2423.40，β＝-0.479。

6、预警：如图1所示，已经拟合出泥石流爆发率为80％的曲线，预警时，若当前获得的关键气象因子落到曲线以上区域(包括曲线上)，则按照曲线以上的泥石流事件出现统计概率进行相应预警，即进行泥石流发生概率在80％以上的预警。

7、分级预警：拟合了不同泥石流出现概率下的曲线模型，则可实现分级预警。如拟合泥石流出现概率分别为80％、50％的两条曲线模型，如图1所示，预警时，根据计算得出的关键气象因子落在图上的位置关系，可进行分级预警。

值得说明的是，后续随着日后泥石流事件样本量的增加，泥石流事件与非泥石流事件的可区分性会越来越强，分级预警程度将更为精细化，预报的准确性将得到持续提升。

本发明解决了现有技术一直未能解决的冰川泥石流地区定量预警的问题，提供了一种能随着泥石流数据的增多持续有效性提升预警结果的冰川泥石流区域预警曲线和预警方法，为今后对冰川泥石流区域预警提供了量化指导和理论基础。

本发明中，未详细描述的均是现有技术。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。