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

本发明专利技术公开了一种基于水热组合的冰川泥石流区域预警方法，包括以下几个步骤：S1、在冰川泥石流预警建模研究区域内采集冰川泥石流相关数据；S2、构建预警阈值曲线模型；S4、根据S3构建的预警阈值曲线，结合需预警地的气象条件是否落在预警阈值曲线上的位置关系，进行预警：若气象条件落在曲线及其以上区域，则按照曲线及其以上的泥石流事件出现统计概率进行相应预警。本发明专利技术解决了现有技术一直未能解决的冰川泥石流地区定量预警的问题，提供了一种能随着泥石流数据的增多持续有效性提升预警结果的冰川泥石流区域预警曲线和预警方法，为今后对冰川泥石流区域预警提供了量化指导和理论基础，具有普适意义。

【技术实现步骤摘要】
一种基于水热组合的冰川泥石流区域预警方法
本专利技术涉及冰川泥石流领域，具体涉及一种基于水热组合的冰川泥石流区域预警方法。
技术介绍
随着经济建设和社会发展需要，人们逐步在高海拔冰川山地地区进行工程建设，在这工程建设过程中，不可避免的会遇到冰川泥石流灾害，因此有必要对冰川泥石流发生建立预警机制。相比雨季的暴雨泥石流，冰川泥石流本身属于更小概率的事件，其起动机理研究仍处于起步阶段，再加上目前冰川泥石流数据存在样本量小、样本噪声大的问题，现有冰川泥石流预警则更多是停留在局地经验取值的层面或群测群防(即监测员巡查)，因此，无论是实际应用还是工程设计规范推荐的泥石流相关计算方法都是以简单经验性的公式计算，缺乏区域普适性应用价值。复杂的方法(如机器学习)在这个方向上，由于多数研究仅为针对单个冰川或少数次泥石流事件的分析，无法建立区域普适性的冰川泥石流有效预警方法。现阶段急需一种应对冰川泥石流区域预警的普适性有效预警方法。
技术实现思路
为解决上述问题，本专利技术目的在于提供一种基于水热组合的冰川泥石流区域预警方法，该方法解决了现有技术一直未能解决的冰川泥石流地区定量预警的问题，提供了一种能随着泥石流数据的增多持续有效性提升预警结果准确性的冰川泥石流区域预警曲线和预警方法，为今后对冰川泥石流区域预警提供了量化指导和理论基础，具有普适意义。本专利技术通过下述技术方案实现：一种基于水热组合的冰川泥石流区域预警方法，其特征在于，包括以下几个步骤：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-本文档来自技高网...

【技术保护点】
1.一种基于水热组合的冰川泥石流区域预警方法，其特征在于，包括以下几个步骤：S1、在冰川泥石流预警建模研究区域内采集冰川泥石流相关数据；S2、构建预警阈值曲线模型：SS1、基于S1中采集的冰川泥石流相关数据，选取有效累积雨量R

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


2.根据权利要求1所述的基于水热组合的冰川泥石流区域预警方法，其特征在于，所述冰川泥石流预警建模研究区域的选择原则包括：(1)、冰川泥石流所处的研究区内多数流域上游物源区正处...

【专利技术属性】
技术研发人员：刘双，胡凯衡，王璋，张晓鹏，
申请(专利权)人：中国科学院，水利部成都山地灾害与环境研究所，
类型：发明
国别省市：四川;51