大坝变形异常值识别方法及系统技术方案

本发明专利技术属于大坝安全监测技术领域，具体公开了一种大坝变形异常值识别方法及系统，解决了现有大坝变形异常值识别方法存在的识别结果不准确的技术问题

【技术实现步骤摘要】
大坝变形异常值识别方法及系统


[0001]本专利技术公开了一种大坝变形异常值识别方法及系统，具体为一种基于时空聚类和深度学习的大坝变形异常值识别方法及系统，属于大坝安全监测

。

技术介绍

[0002]安全监测是定量分析大坝安全性态的重要方法
。
大坝众多监测科目中，变形能够真实
、
直观
、
准确地反映大坝安全性态变化，是大坝最重要的监测内容
。
然而，采集的变形监测数据受仪器故障
、
人为因素或其他不确定因素的影响，不可避免地会出现各种异常值
。
[0003]异常数据可能是荷载急剧变化造成的正常结构反映，也可能是仪器故障
、
人为因素等造成的粗差错误数据
。
因此，对大坝变形异常值进行精准识别，对大坝监测数据分析及安全评价研究具有重要意义
。
[0004]现有技术通常是利用统计检验法和监控模型法识别大坝变形异常值
。
但统计检验法仅考虑变形数据本身，未结合大坝变形输入输出关系，容易错误将正常数据误诊断为异常值
。
监控模型法主要考虑单测点监测数据的异常值，忽略了与相似测点的时空关联特性，导致识别结果出错
。

技术实现思路

[0005]本申请的目的在于，提供一种大坝变形异常值识别方法及系统，以解决现有技术中大坝变形异常值识别方法存在的识别结果不准确的技术问题
。
[0006]本专利技术的第一方面提供了一种大坝变形异常值识别方法，包括：
[0007]对大坝上多个测点进行聚类分区，并获取每个分区内每个测点的预测值和实测值之间的残差；
[0008]根据每个分区内所有测点的残差，确定每个分区的综合残差；
[0009]结合所述预测值和所述综合残差，确定每个分区内测点的实测值中的异常值
。
[0010]优选地，对大坝上多个测点进行聚类分区，具体包括：
[0011]获取大坝上每一个测点与其他测点之间的综合距离；
[0012]根据所述综合距离，利用最小离差平方和系统聚类法对多个测点进行聚类分区
。
[0013]优选地，获取大坝上每一个测点与其他测点之间的综合距离，具体包括：
[0014]获取大坝上每一个测点与其他测点之间的绝对距离
、
增量距离和增速距离；
[0015]根据所述绝对距离
、
增量距离和增速距离，确定每一个测点与其他测点之间的综合距离
。
[0016]优选地，所述绝对距离根据一个测点的变形值与另一个测点的变形值之间的差值确定；
[0017]所述增量距离根据一个测点相较于前一时刻的变形增量与另一个测点相较于前一时刻的变形增量之间的差值确定；
[0018]所述增速距离根据一个测点的变形增速与另一个测点的变形增速之间的差值确
定
。
[0019]优选地，获取每个分区内每个测点的预测值和实测值之间的残差，具体包括：
[0020]优化每个分区的预测模型的超参数；
[0021]根据每个分区的环境影响因子
、
每个测点的历史测值和所述超参数，训练对应分区的所述预测模型；
[0022]获取预测模型对对应每个分区内每个测点的预测值，并确定每个测点的预测值和实测值之间的残差
。
[0023]优选地，优化每个分区的预测模型的超参数，具体包括：
[0024]利用蜜獾算法优化每个分区的预测模型的超参数
。
[0025]优选地，所述预测模型为双向长短时记忆神经网络模型
。
[0026]优选地，根据每个分区内所有测点的残差，确定每个分区的综合残差，具体包括：
[0027]获取每个分区内每一个测点的残差与其余测点的残差之间的相关系数；
[0028]根据所述分区内的所有所述相关系数，确定每个分区的综合残差
。
[0029]优选地，结合所述预测值和所述综合残差，确定每个分区内测点的实测值中的异常值，具体包括：
[0030]利用所述综合残差确定每个分区的基准值；
[0031]将每个分区内的每个测点的预测值与实测值之间的差值与所述基准值进行对比，确定每个分区内测点的实测值中的异常值
。
[0032]本专利技术的第二方面提供了一种大坝变形异常值识别系统，包括：
[0033]残差模块，所述残差模块用于对大坝上多个测点进行聚类分区，并获取每个分区内每个测点的预测值和实测值之间的残差；
[0034]综合残差模块，所述综合残差模块用于根据每个分区内所有测点的残差，确定每个分区的综合残差；
[0035]异常值确定模块，所述异常值确定模块用于结合所述预测值和所述综合残差，确定每个分区内测点的实测值中的异常值
。
[0036]本专利技术的大坝变形异常值识别方法及系统，相较于现有技术，具有如下有益效果：
[0037]本专利技术方法利用最小离差平方和系统聚类法完成测点变形聚类分区，采用蜜獾算法对预测模型进行超参数优化，基于最优超参数进行建模预测
。
然后，在相似测点残差相关性的基础上构建异常值判别指标，当研究测点的监测数据序列残差大于指标标准，即可识别为异常值
。
[0038]本专利技术方法考虑大坝变形的时空关联特性并提出高精度的大坝变形预测模型，可以更加合理地识别出大坝的异常测值，提升了对变形监测数据的异常值识别的精度，为大坝变形提供了一种高精度预测及异常值识别方法
。
附图说明
[0039]图1为本专利技术大坝变形异常值识别方法的流程图；
[0040]图2为本专利技术实施例中气温与上游水位时间序列过程线；
[0041]图3为本专利技术实施例中坝体测点
Ward
聚类分区结果图；
[0042]图4为本专利技术实施例中坝体各分区测点径向水平位移过程线，其中
(a)
为图3中的
分区
I
内各测点径向水平位移过程线；
(b)
为图3中的分区Ⅱ内各测点径向水平位移过程线；
(c)
为图3中的分区Ⅲ内各测点径向水平位移过程线；
(d)
为图3中的分区Ⅳ内各测点径向水平位移过程线；
(e)
为图3中的分区
Ⅴ
内各测点径向水平位移过程线；
[0043]图5中的
(a)
为本专利技术实施例中各模型在
PL3
‑1测点预测结果与实际变形序列对比图，
(b)
为本专利技术实施例中各模型在
PL3
‑1测点的各评价指标分析对比图；
[0044]图6中的
(a)
为本专利技术实施例中各模型在
PL4
‑1测点预测结果与实际变形序列对比图，
(b本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种大坝变形异常值识别方法，其特征在于，包括：对大坝上多个测点进行聚类分区，并获取每个分区内每个测点的预测值和实测值之间的残差；根据每个分区内所有测点的残差，确定每个分区的综合残差；结合所述预测值和所述综合残差，确定每个分区内测点的实测值中的异常值
。2.
根据权利要求1所述的大坝变形异常值识别方法，其特征在于，对大坝上多个测点进行聚类分区，具体包括：获取大坝上每一个测点与其他测点之间的综合距离；根据所述综合距离，利用最小离差平方和系统聚类法对多个测点进行聚类分区
。3.
根据权利要求2所述的大坝变形异常值识别方法，其特征在于，获取大坝上每一个测点与其他测点之间的综合距离，具体包括：获取大坝上每一个测点与其他测点之间的绝对距离
、
增量距离和增速距离；根据所述绝对距离
、
增量距离和增速距离，确定每一个测点与其他测点之间的综合距离
。4.
根据权利要求3所述的大坝变形异常值识别方法，其特征在于，所述绝对距离根据一个测点的变形值与另一个测点的变形值之间的差值确定；所述增量距离根据一个测点相较于前一时刻的变形增量与另一个测点相较于前一时刻的变形增量之间的差值确定；所述增速距离根据一个测点的变形增速与另一个测点的变形增速之间的差值确定
。5.
根据权利要求1所述的大坝变形异常值识别方法，其特征在于，获取每个分区内每个测点的预测值和实测值之间的残差，具体包括：优化每个分区的预测模型的超参数；根据每个分区的环境影响因子
、
每个测点的历史测值...

【专利技术属性】
技术研发人员：费新峰，刘海敏，张治，张小宁，岳子，
申请(专利权)人：中国电建集团西北勘测设计研究院有限公司，
类型：发明
国别省市：