一种基于多源监测数据的滑坡稳定性智能动态预测方法

1.本发明涉及地质灾害防治方法领域，具体涉及一种基于多源监测数据的滑坡稳定性智能动态预测方法。


背景技术：

2.滑坡是斜坡上的土体或者岩体，受河流冲刷、地下水活动、雨水浸泡、地震及人工切坡等因素影响，在重力作用下，沿着一定的软弱面或者软弱带，整体地或者分散地顺坡向下滑动的自然灾害。滑坡是我国最严重的自然灾害之一，滑坡自然灾害每年造成人员伤亡超过400人，经济损失超过10亿元。因此，对滑坡进行有效防控具有重要意义。
3.滑坡防控的一个重要内容，是准确判定滑坡的稳定状态。针对滑坡的稳定状态判定是通过已有的边坡稳定性计算方法进行，边坡稳定性计算方法主要包括极限平衡法、极限分析法、强度折减法和数值模拟方法等，这些方法所采用的参数主要为岩土体的静态物理力学参数例如坡高、坡脚、弹性模量、粘聚力、内摩擦角等，无法反映滑坡体稳定性的动态发展过程。现有的对于滑坡稳定性进行预测的方法主要借助数个单一的机器学习模型，而未进一步融合不同机器学习模型的预测性能。而现有以极限平衡法等传统稳定性计算方法为基础提出的滑坡稳定性预测方法，也未考虑滑坡稳定性的动态变化，导致预测结果不够准确。


技术实现要素：

4.本发明意在提供一种基于多源监测数据的滑坡稳定性智能动态预测方法，以解决现有方法未考虑滑坡稳定性的动态变化及没有融合滑坡现场实时的多源监测数据，导致预测结果不够准确的问题。
5.本方案中的基于多源监测数据的滑坡稳定性智能动态预测方法，包括以下步骤：
6.步骤1，获取目标滑坡隐患点地质背景条件和岩土体相关的物理学参数；
7.步骤2，获取目标滑坡隐患点传感器的历史监测数据，对历史监测数据进行预处理操作，形成目标滑坡隐患点多源监测数据集合，所述预处理操作包括去除数据噪声、删除异常数据和空缺值填补；
8.步骤3，基于步骤1获取的目标滑坡隐患点地质背景条件和岩土体相关的物理学参数，通过数值模拟手段建立数值模型；
9.步骤4，获取目标滑坡隐患点的历史工况信息，基于步骤3建立的数值模型，根据工况信息设定滑坡数值模拟过程的边界条件进行数值分析计算，得到目标滑坡隐患点历史的渗流场和应力场的变化，基于计算得到的目标滑坡隐患点不同时间阶段的渗流场和应力场，结合滑坡稳定性计算方法计算获得目标滑坡隐患点不同时间阶段的稳定系数；
10.步骤5，将步骤2中目标滑坡隐患点的多源监测数据集合与步骤4中的目标滑坡隐患点在不同时间阶段的稳定系数根据时间一一对应，形成供输入机器学习模型进行训练和验证的数据集，其中，步骤2中目标滑坡隐患点的多源监测数据集合作为机器学习模型输入
的特征变量，与之对应的步骤4中目标滑坡隐患点稳定系数作为机器学习模型的标签；
11.步骤6，选取机器学习任务中广泛采用的基础机器学习模型作为集成学习框架个体学习器的备选，利用步骤5中机器学习模型训练和验证的数据集，对n个基础机器学习模型进行训练、验证和测试，同时生成供后续集成学习堆叠模型训练和测试的训练集和测试集；
12.步骤7，从步骤6中n个基础机器学习模型中选取对稳定系数预测精度最高的m个基础机器学习模型作为集成学习的基学习器，用于后续集成学习元学习器的融合；
13.步骤8，基于集成学习堆叠策略的思想，通过元学习器对步骤7中选取的m个基学习器进行融合，形成集成学习堆叠模型，并利用步骤6中生成的训练集和测试集对集成学习堆叠模型进行训练和测试，得到性能最优的集成学习堆叠模型；
14.步骤9，将目标滑坡隐患点各个传感器的实时多源监测数据输入步骤8中训练好的集成学习堆叠模型，得到多源监测数据融合下目标滑坡隐患点实时的稳定系数。
15.本方案的有益效果是：
16.融合滑坡现场实时的多源监测数据与数值模拟技术，并基于集成学习理论，实现了基于滑坡现场实时的多源监测数据对滑坡稳定性状态的动态预测，充分考虑了滑坡体的动态演化过程，其利用集成学习的堆叠策略对基础机器学习模型的预测性能进行融合增强，耦合了不同机器学习模型的非线性表征能力，有效提升了滑坡稳定预测的准确度，从而实现了滑坡稳定性实时、高效、准确的预测，为滑坡灾害防治提供支撑。
17.进一步，所述步骤3中，所述数值模拟手段为有限元、离散元、物质点法中的一种。
18.有益效果是：通过数值模拟手段中的一种，适用范围更大。
19.进一步，所述步骤4中，所述滑坡稳定性计算方法为极限平衡法、极限分析法、强度折减法、数值分析法中的一种。
20.有益效果是：通过滑坡稳定性计算方法得到滑坡的不同时间段的稳定系数，多种滑坡稳定性计算方法，适用范围更大。
21.进一步，所述步骤6中，所述广泛采用的基础机器学习模型包括经典的回归模型和基于决策树的回归模型，所述经典的回归模型包括最小二乘回归、岭回归、核岭回归、支持向量回归、最近邻回归，所述基于决策树的回归模型包括决策树、随机森林、梯度提升回归树、极限梯度提升、轻度梯度提升机。
22.有益效果是：广泛考虑了各种的常见的机器学习模型的回归性能，为集成学习框架的搭建奠定基础。
23.进一步，所述步骤6中，对n个基础机器学习模型进行训练、验证和测试，其中模型超参数的优化采用k折交叉验证和连续减半搜索方法，包括以下子步骤：
24.子步骤6.1，将步骤5中机器学习模型训练和验证的数据集按照预设比例划分为训练集和测试集；
25.子步骤6.2，将子步骤6.1中的训练集打乱并平均分为数量相同的k组，对于一个基础机器学习模型，每轮利用其中一组数据作为验证，其余k-1组作为训练，并将机器学习模型在测试集上预测一次，记录每轮的预设参数，不断调整基础机器学习模型的超参数组合，所述机器学习模型的超参数组合调整的调优方法为连续减半搜索；
26.子步骤6.3，针对一个基础机器学习模型，将子步骤6.2重复k次，得到k组验证数据
和k组测试数据，将k组验证数据纵向排列，将k组测试数据取平均值，在k个不同超参数的基础机器学习模型中选择预设参数达到指定条件的模型作为该模型的最优；
27.子步骤6.4，重复子步骤6.2和6.3直到n个基础机器学习模型全部调优完毕，将子步骤6.3中每个基础机器学习模型纵向排列的k组验证数据和取平均的测试数据分别横向排列，作为元学习器新的训练集和测试集。
28.有益效果是：通过k折交叉验证能够综合比较同一基础机器学习模型在不同超参数组合情况下的预测性能，再利用连续减半搜索方法进行较少次数的计算，并为最优的基础机器学习模型提供较多的计算资源分配。
29.进一步，所述子步骤6.1中，所述预设比例按照数据量大小进行确定，所述预设比例包括7:3和8:2。
30.有益效果是：通过预设比例的设置，能够对数据集进行合理的划分，提高数据处理的速度，同时提高预测结果的准确性。
31.进一步，所述子步骤6.2中，所述预设参数包括验证集误差、准确率或决定系数r2，所述验证集误差最小即为预设参数达到指定条件，或者所述准确率最高即为预设参数达到指定条件，或者所述决定系数r2最高即为预设参数达到指定条件。
32.有益效果是：通过多个预设参数和对应的指定条件设置，能够准确获取到机器学习模型超参数的最优组合，让预测的稳定系数更加准确。
33.进一步，所述步骤8中，所述元学习器包括线性回归和多层感知回归。
34.有益效果是：通过元学习器的设置，能够增强集成学习堆叠模型的泛化性能并减少过拟合的风险。
附图说明
35.图1为本发明基于多源监测数据的滑坡稳定性智能动态预测方法实施例的流程框图。
具体实施方式
36.下面通过具体实施方式进一步详细说明。
37.实施例一
38.基于多源监测数据的滑坡稳定性智能动态预测方法，如图1所示，包括以下步骤：
39.步骤1，获取目标滑坡隐患点地质背景条件和岩土体相关的物理学参数，基于勘察和试验手段获取物理学参数，物理学参数包括弹性模量、内摩擦角、粘聚力、饱和/残余含水量、渗透系数。
40.步骤2，获取目标滑坡隐患点现有传感器的历史监测数据，包括降雨量、库水位高程、地表位移、深部位移、裂缝宽度、孔隙水压力、地下水位、应力。对历史监测数据进行预处理操作，形成目标滑坡隐患点多源监测数据集合，预处理操作包括去除数据噪声、删除异常数据和空缺值填补。
41.步骤3，基于步骤1获取的目标滑坡隐患点地质背景条件和岩土体相关的物理学参数，通过数值模拟手段建立1:1的数值模型，数值模拟手段为有限元、离散元、物质点法中的一种，地质背景条件包括环境诱发因素和关键因素，环境诱发因素如降雨和库水位等，环境
诱发因素用ei表示，i表示第i个独立监测参数，ei可以为环境诱发因素的实时数据，如实时的降雨量、降雨强度、温度、库水位高程等，用来作为数值模拟过程中的输入和边界条件，该类数据为对目标滑坡隐患点环境条件的直接观测，关键因素是一组用于表征目标滑坡隐患点稳定状态的关键参数，用fj表示，j表示模型的第j个关键参数，如位移、孔隙水压力、地下水位等，关键因素如地表位移、深部位移、裂缝宽度、孔隙水压力、地下水位、应力等，将该类数据将作为表征目标滑坡隐患点实时稳定状态的关键参数，为后续机器学习模型的训练做准备。
42.步骤4，获取目标滑坡隐患点的历史工况信息，基于步骤3建立的数值模型，根据工况信息设定滑坡数值模拟过程的边界条件进行数值分析计算，得到目标滑坡隐患点历史的渗流场和应力场的变化，基于计算得到的目标滑坡隐患点不同时间阶段的渗流场和应力场，结合滑坡稳定性计算方法计算获得目标滑坡隐患点不同时间阶段的稳定系数，滑坡稳定性计算方法为极限平衡法、极限分析法、强度折减法、数值分析法中的一种。
43.步骤5，将步骤2中目标滑坡隐患点的多源监测数据集合与步骤4中的目标滑坡隐患点在不同时间阶段的稳定系数根据时间一一对应，形成供输入机器学习模型进行训练和验证的数据集，其中，步骤2中目标滑坡隐患点的多源监测数据集合作为机器学习模型输入的特征变量，与之对应的步骤4中目标滑坡隐患点稳定系数作为机器学习模型的标签。
44.步骤6，选取机器学习任务中广泛采用的基础机器学习模型作为集成学习框架个体学习器的备选，利用步骤5中机器学习模型训练和验证的数据集，对n个基础机器学习模型进行训练、验证和测试，同时生成供后续集成学习堆叠模型训练和测试的训练集和测试集，广泛采用的基础机器学习模型包括经典的回归模型和基于决策树的回归模型，所述经典的回归模型包括最小二乘回归、岭回归、核岭回归、支持向量回归、最近邻回归，所述基于决策树的回归模型包括决策树、随机森林、梯度提升回归树、极限梯度提升、轻度梯度提升机。
45.对n个基础机器学习模型进行训练、验证和测试，其中模型超参数的优化采用k折交叉验证和连续减半搜索方法，包括以下子步骤：
46.子步骤6.1，将步骤5中机器学习模型训练和验证的数据集按照预设比例划分为训练集和测试集，所述预设比例包括7:3和8:2。
47.子步骤6.2，将子步骤6.1中的训练集打乱并平均分为数量相同的k组，对于一个基础机器学习模型，每轮利用其中一组数据作为验证，其余k-1组作为训练，并将机器学习模型在测试集上预测一次，记录每轮的预设参数，不断调整基础机器学习模型的超参数组合，所述机器学习模型的超参数组合调整的调优方法为连续减半搜索。
48.子步骤6.3，针对一个基础机器学习模型，将子步骤6.2重复k次，得到k组验证数据和k组测试数据，将k组验证数据纵向排列，将k组测试数据取平均值，在k个不同超参数的基础机器学习模型中选择预设参数达到指定条件的模型作为该模型的最优。
49.子步骤6.4，重复子步骤6.2和6.3直到n个基础机器学习模型全部调优完毕，将子步骤6.3中每个基础机器学习模型纵向排列的k组验证数据和取平均的测试数据分别横向排列，作为元学习器新的训练集和测试集。
50.步骤7，从步骤6中n个基础机器学习模型中选取对稳定系数预测精度最高的m个基础机器学习模型作为集成学习的基学习器，用于后续集成学习元学习器的融合。元学习器
包括线性回归和多层感知回归。
51.步骤8，基于集成学习堆叠策略的思想，通过元学习器对步骤7中选取的m个基学习器进行融合，形成集成学习堆叠模型，并利用步骤6中生成的训练集和测试集对集成学习堆叠模型进行训练和测试，得到性能最优的集成学习堆叠模型。
52.骤9，将目标滑坡隐患点各个传感器的实时多源监测数据输入步骤8中训练好的集成学习堆叠模型，得到多源监测数据融合下目标滑坡隐患点实时的稳定系数。
53.本实施例对于滑坡稳定性的动态预测方法相较于现有的预测方法来说，融合滑坡现场实时的多源监测数据与数值模拟技术，并基于集成学习理论，实现了基于滑坡现场实时的多源监测数据对滑坡稳定性状态的动态预测，充分考虑了滑坡体的动态演化过程，其利用集成学习的堆叠策略对基础的机器学习模型的预测性能进行融合，耦合了不同机器学习模型的非线性表征能力，能够及时、连续地预测到滑坡产生的动态变化，有效提升了滑坡稳定预测的准确度，从而实现了滑坡稳定性实时、高效、准确的预测，能够根据滑坡动态演化的估计对滑坡实际产生后的影响范围进行预测，为滑坡灾害防治提供支撑。
54.针对现有对滑坡的预测方法，对预测结果稳定性问题的改进设计过程中，因考虑到各项参数的普通处理器的处理能力、处理速度和处理计算量，通常是对单一的机器学习模型进行相应改进，或者是稳定性计算方法对滑坡是否稳定的最终结果进行预测，而不会将多方面的监测数据并结合集成学习理论，进行动态预测滑坡稳定性状态，本实施例在此基础上，通过设置预设比例进行机器学习训练，能够在多个参数以及结合集成学习理论后降低处理所需处理器的需求，减小计算量，提高预测结果的速度和及时性。
55.以上所述的仅是本发明的实施例，方案中公知的具体结构及特性等常识在此未作过多描述。应当指出，对于本领域的技术人员来说，在不脱离本发明结构的前提下，还可以作出若干变形和改进，这些也应该视为本发明的保护范围，这些都不会影响本发明实施的效果和专利的实用性。本技术要求的保护范围应当以其权利要求的内容为准，说明书中的具体实施方式等记载可以用于解释权利要求的内容。