一种浅层崩塌灾害预警模型建立方法、装置、设备及介质制造方法及图纸

本发明专利技术公开一种浅层崩塌灾害预警模型建立方法、装置、设备及介质，方法包括获取待预警区域内的历史灾点数据；基于数字高程模型划分斜坡单元，提取多个孕灾因子；基于梯度增强决策树构建崩塌易发性预测模型并评价其准确性，得到易发性分区评价结果；根据易发性分区评价结果，统计分析不同易发区内历史灾点当日雨量及前期雨量，通过相关系数确定前期降雨天数，建立前期降雨

【技术实现步骤摘要】
一种浅层崩塌灾害预警模型建立方法、装置、设备及介质


[0001]本专利技术涉及岩土灾害预警
，尤其涉及一种浅层崩塌灾害预警模型建立方法、装置、设备及介质。

技术介绍

[0002]浅层土质崩塌多发生于汛期，具有“点多面广”及“体量小预警难”等特征，给山区村民的生命财产安全带来严重威胁，也给防灾减灾工作造成麻烦。在每年汛期到来之际，大规模转移村民的现象时有发生，不仅浪费大量的人力物力，也给村民和防灾人员造成严重困扰。其归根的原因在于浅层土质崩塌的预警阈值不准，要么不够灵敏，要么过于敏感。
[0003]目前，关于浅层土质崩塌，如滑坡、崩塌或泥石流的预警方法，主要采用：统计分析、数值模拟及物理模型。其中，数值模拟或物理模型方法主要用于灾害启动机理研究及单点灾害预警阈值的确定，更多的是关注理论层面，但在推广到工程实践及实际应用上存在着不足。而基于统计分析所建立的模型虽然缺乏机理支撑，但经过多年的发展，得到了世界各国的广泛应用。统计模型的建立通常只需要充足的雨量资料和灾点信息作为支撑，并未考虑坡度、植被和岩性等其他影响因素对斜坡失稳的影响。
[0004]综上，现有技术中的预警模型建立方法未考虑斜坡失稳的影响，导致其预警阈值的准确度和可靠性不高。

技术实现思路

[0005]本专利技术提供了一种浅层崩塌灾害预警模型建立方法、装置、设备及介质，以解决现有技术中的预警模型建立方法未考虑斜坡失稳的影响，导致其预警阈值的准确度和可靠性不高的技术问题。
[0006]第一方面，为了解决上述技术问题，本专利技术提供了一种浅层崩塌灾害预警模型建立方法，包括：
[0007]获取待预警区域内的历史灾点数据，所述历史灾点数据包括发灾位置、发灾时间和发灾时临近雨量站数据；
[0008]基于数字高程模型划分斜坡单元，并基于ArcGIS提取多个孕灾因子；
[0009]以所述历史灾点数据和所述孕灾因子作为正样本集，选取相同数量的非崩塌点和孕灾因子作为负样本集，组成总样本集并按比例划分为训练集、验证集和测试集；
[0010]基于梯度增强决策树构建崩塌易发性预测模型并评价其准确性，得到易发性分区评价结果；
[0011]根据所述易发性分区评价结果，统计分析不同易发区内历史灾点当日雨量及前期雨量，通过相关系数确定前期降雨天数，建立前期降雨
‑
当日降雨阈值模型；
[0012]根据坡体是否被治理、植被条件及现阶段变形迹象，对所述斜坡单元划分易发性等级，并修正所述斜坡单元的易发性等级及降雨阈值；
[0013]根据所述斜坡单元的易发性等级及降雨阈值，构建多等级时空预警模型。
[0014]优选地，所述基于数字高程模型划分斜坡单元，包括：
[0015]生成无洼地，计算流向和流量；
[0016]生成河网并将河流连接，生成分水岭并将分水岭转面，生成正向集水流域面，反转原DEM；
[0017]重复所述生成河网并将河流连接，生成分水岭并将分水岭转面，生成正向集水流域面的步骤，生成反向集水流域面；
[0018]合并所述正向集水流域面和所述反向集水流域面，生成初始斜坡单元；
[0019]对不合理边界进行修正，得到斜坡单元。
[0020]优选地，所述基于ArcGIS提取多个孕灾因子，包括：
[0021]利用ArcGIS中的欧几里德距离工具计算道路、断层和水系到各斜坡因子的距离；
[0022]取所有栅格单元距离的平均值作为斜坡单元的最终值；
[0023]其中，在高程、坡度、坡向、高差、平面曲率、剖面曲率、TWI和NDVI提取时，按一个斜坡单元内所有栅格单元的平均值作为各个斜坡单元的最终赋值；在利用斜坡单元进行地层岩性提取时，若一个单元内涉及多个岩性，则按照地层岩性所占面积最大的岩性作为斜坡因子的最终赋值。
[0024]优选地，多个所述孕灾因子包括：高程、坡度、坡向、高差、平面曲率、剖面曲率、地形湿度指数TWI、归一化植被指数NDVI、岩性、到道路距离、到断层距离及到水系距离。
[0025]优选地，所述以所述历史灾点数据和所述孕灾因子作为正样本集，选取相同数量的非崩塌点和孕灾因子作为负样本集，组成总样本集并按比例划分为训练集、验证集和测试集，包括：
[0026]以所述历史灾点数据和所述孕灾因子作为正样本集；
[0027]在距历史灾点半径缓冲区外，随机选取相同数量的非崩塌点，构建负样本集；
[0028]根据所述正样本集和所述负样本集组成总样本集；
[0029]按照8：1：1的比例将总样本集划分为训练集、验证集和测试集，分别用于模型的训练、调参和测试。
[0030]优选地，所述基于梯度增强决策树构建崩塌易发性预测模型并评价其准确性，得到易发性分区评价结果，包括：
[0031]初始化弱分类器：
[0032][0033]f0(x)为初始弱分类器，L(y
i
,γ)为一个可微的损失函数；
[0034]计算负梯度：
[0035][0036]在回归树中，实则就是y
‑
c，c指的是每个叶子节点样本标签的平均值；
[0037]将公式(2)中的残差作为样本的新值，将(x
m
,r
im
),m＝1,2,3,
…
,M作为下一棵树的训练数据；
[0038]根据所述样本的新值，训练得到新的弱学习器G
k
(X)为：
[0039][0040]返回至所述计算负梯度的步骤，直到满足最低误差要求，最终形成的强分类器为：
[0041][0042]式中，f
m
(x)是崩塌易发性预测模型，f0(x)为初始模型，其值为损失函数的最小值；m为最大迭代次数，J是第j棵回归树的叶节点数，γ是损失函数在叶节点区域负梯度的最优拟合值，对应的第j棵回归树的叶节点区域为R
jm
(j＝1，2，3，
…
，J)，/为指示函数，R
jm
是第t个决策树对应的叶节点区域；最大迭代次数(m)设置为150次，考虑数据时学习率为0.5，叶节点数(J)设置为20。
[0043]优选地，所述根据所述易发性分区评价结果，统计分析不同易发区内历史灾点当日雨量及前期雨量，通过相关系数确定前期降雨天数，建立前期降雨
‑
当日降雨阈值模型，包括：
[0044]分别统计不同易发性分区内已发灾点数量、灾害发生时的当日雨量及前九天雨量；
[0045]通过相关性分析，确定前期降雨天数，建立不同易发性分区的前期降雨
‑
当日降雨阈值模型；
[0046]其中，所述相关性系数的公式为：
[0047][0048]式中，r(X，Y)为两变量之间的相关程度；cov(X，Y)为X与Y的协方差；var(X)为x的方差；var(Y)为y的方差。
[0049]第二方面，本专利技术提供了一种浅层崩塌灾本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种浅层崩塌灾害预警模型建立方法，其特征在于，包括：获取待预警区域内的历史灾点数据，所述历史灾点数据包括发灾位置、发灾时间和发灾时临近雨量站数据；基于数字高程模型划分斜坡单元，并基于ArcGIS提取多个孕灾因子；以所述历史灾点数据和所述孕灾因子作为正样本集，选取相同数量的非崩塌点和孕灾因子作为负样本集，组成总样本集并按比例划分为训练集、验证集和测试集；基于梯度增强决策树构建崩塌易发性预测模型并评价其准确性，得到易发性分区评价结果；根据所述易发性分区评价结果，统计分析不同易发区内历史灾点当日雨量及前期雨量，通过相关系数确定前期降雨天数，建立前期降雨
‑
当日降雨阈值模型；根据坡体是否被治理、植被条件及现阶段变形迹象，对所述斜坡单元划分易发性等级，并修正所述斜坡单元的易发性等级及降雨阈值；根据所述斜坡单元的易发性等级及降雨阈值，构建多等级时空预警模型。2.根据权利要求1所述的浅层崩塌灾害预警模型建立方法，其特征在于，所述基于数字高程模型划分斜坡单元，包括：生成无洼地，计算流向和流量；生成河网并将河流连接，生成分水岭并将分水岭转面，生成正向集水流域面，反转原DEM；重复所述生成河网并将河流连接，生成分水岭并将分水岭转面，生成正向集水流域面的步骤，生成反向集水流域面；合并所述正向集水流域面和所述反向集水流域面，生成初始斜坡单元；对不合理边界进行修正，得到斜坡单元。3.根据权利要求1所述的浅层崩塌灾害预警模型建立方法，其特征在于，所述基于ArcGIS提取多个孕灾因子，包括：利用ArcGIS中的欧几里德距离工具计算道路、断层和水系到各斜坡因子的距离；取所有栅格单元距离的平均值作为斜坡单元的最终值；其中，在高程、坡度、坡向、高差、平面曲率、剖面曲率、TWI和NDVI提取时，按一个斜坡单元内所有栅格单元的平均值作为各个斜坡单元的最终赋值；在利用斜坡单元进行地层岩性提取时，若一个单元内涉及多个岩性，则按照地层岩性所占面积最大的岩性作为斜坡因子的最终赋值。4.根据权利要求3所述的浅层崩塌灾害预警模型建立方法，其特征在于，多个所述孕灾因子包括：高程、坡度、坡向、高差、平面曲率、剖面曲率、地形湿度指数TWI、归一化植被指数NDVI、岩性、到道路距离、到断层距离及到水系距离。5.根据权利要求1所述的浅层崩塌灾害预警模型建立方法，其特征在于，所述以所述历史灾点数据和所述孕灾因子作为正样本集，选取相同数量的非崩塌点和孕灾因子作为负样本集，组成总样本集并按比例划分为训练集、验证集和测试集，包括：以所述历史灾点数据和所述孕灾因子作为正样本集；在距历史灾点半径缓冲区外，随机选取相同数量的非崩塌点，构建负样本集；根据所述正样本集和所述负样本集组成总样本集；
按照8：1：1的比例将总样本集划分为训练集、验证集和测试集，分别用于模型的训练、调参和测试。6.根据权利要求1所述的浅层崩塌灾害预警模型建立方法，其特征在于，所述基于梯度增强决策树构建崩塌易发性预测模型并评价其准确性，得到易发性分区评价结果，包括：初始化弱分类器：f0(x)为初始弱分类器，L(y
i
,γ)为一个可微的损失函数；计算负梯度：在回归树中，实则就是y
‑
c，c指的是每个叶子节点样本标签的平均值；将公式...

【专利技术属性】
技术研发人员：梁柱，刘伟，张庆华，张守文，姜凯华，程伟，黄厚赞，
申请(专利权)人：广州市城市规划勘测设计研究院，
类型：发明
国别省市：