一种滑坡变形数据实时过滤方法与流程

1.本发明涉及针对在突发型滑坡临滑发生时对临滑变形滑坡的数据误差过滤方法。


背景技术：

2.滑坡变形数据的误差过滤平滑是滑坡时间预报的重要前期工作，确保预报结果不因误差产生严重偏差，近年来常用的移动平均法（孟蒙,陈智强,黄达,等.基于h-p滤波法、arima和var模型的库区滑坡位移综合预测[j].岩土力学,2016,37(s2):552-560.）、最小二乘法（郭献涛,黄腾,臧妻斌,等. 基于最小二乘三维表面匹配算法的滑坡变形测量[j].岩土力学,2015,36(5):1421-1427.）、kalman滤波（彭鹏,单治钢,董育烦,等.基于kalman滤波融合算法的库区滑坡动态变形监测综合信息提取[j].岩石力学与工程学报,2014,33(8):1520-1525.）、神经网络法过滤（翟会君,翟亚锋,朱涛,等.基于回归-elm神经网络模型的滑坡变形及失稳预测模型[j].河北工业科技,2017,34(6):440-447）等多种变形数据处理方法都能在事后较好地平滑变形数据，但对于突发型滑坡临滑实时预报，由于变形监测设备获取的实时数据具有一定的误差不确定性，在误差过滤中需要实时反映波动究竟是由于滑坡进入临滑加速还是精度误差引起，而单一的过滤方法无法较好地平衡误差剔除和及时识别加速变形，在实时预报中过滤程度较低会造成误报，过滤程度过高则会导致漏报，这一问题对临滑预报的偏差影响极大，缺乏科学有效的临滑时间预报方法，对于人员紧急疏散和应急抢险非常不利。


技术实现要素：

[0003]
本发明的目的在于克服现有过滤方法的不足，适应现实需要，提供一种滑坡变形数据实时过滤方法，本发明方法避免了采用最小二乘法导致不能实时反映波动究竟是由于滑坡进入临滑加速还是精度误差引起的问题，为滑坡应急处置提供可靠的时间参考。
[0004]
为了实现本发明的目的，本发明所采用的技术方案为：一种滑坡变形数据实时过滤方法，包括以下步骤：a、通过室内试验和野外现场监测获得的大量变形曲线进行分析，根据原始数据得到的变形速率数据波动趋势特点； b、确定误差基本特征和波动趋势；c、基于实时变形数据的不确定性，采用常规移动平均法、最小二乘法、卡尔曼滤波、样条平滑等常用数据过滤方法相结合，研究不同过滤参数变化下变形数据的误差处理效果和临滑预报偏差特征；d、利用自主研发的智能监测设备在加速变形阶段可自动加密监测频率的优势，建立数据识别过滤区，进行变形数据缓冲过滤。
[0005]
优选的，步骤a具体为：步骤a具体为：依次取24个时间段的原始变形速率时间散点图为一个固定时间段长度，每次顺延一个时间点，即1-24,2-25，3-26逐步往后延后一个时间点，避免其他变量的影响，将第一个固定时间段上的散点用最小二乘法拟合成一条直线
并求出解析式:y=ax b。
[0006]
优选的，步骤c具体为：a、令时间点2即x
2，
令x=x2代入解析式y=ax b，求解y
2.，
得到点（x2,y2）；同时，用最小二乘法求出第二个时间段（2-25）上的散点拟合直线解析式y0=a0x0 b0以及该解析式的中点（x
中
，y
中
），x
中
即为第13个时间点（x
中=
x
13
）；b、将点（x2，y2）与点（x
中
，y
中
）连线，并求得该直线解析式y1=a1x1 b1c、在解析式y1=a1x1 b1的基础上求出该直线上的第二个点的坐标（x3，y3），同时，用最小二乘法求出第三个时间段（3-26）上的散点拟合直线解析式y2=a2x2 b2以及该解析式的中点（x
中1
，y
中1
）.d、将点（x3，y3）与点（x
中1
，y
中1
）连线，并求得该直线解析式y3=a3x3 b3；e、重复cd操作。
[0007]
本发明的有益效果在于：本发明方法考虑变形误差的特点和突发型滑坡的极短加速变形特征，可通过整合现有误差过滤方法，基于丰富的临滑变形监测数据上，构建更合理的突发型滑坡实时过滤方法，为后续滑坡临滑预报提供支持。
附图说明
[0008]
图1原始速率变形图；图2过滤后变形速率图；图3第一个固定时间段所拟合直线解析式示意图；图4第二个固定时间段所拟合直线解析式示意图；图5连点成线示意图；图6第三个固定时间段所拟合直线解析式示意图；图7连点成线示意图。
具体实施方式
[0009]
下面结合附图和实施例对本发明进一步说明：实施例1：一种滑坡变形数据实时过滤方法，参见图3-图7，包括以下步骤：a、依次取24个时间段的原始变形速率时间散点图为一个固定时间段长度，每次顺延一个时间点，即1-24,2-25，3-26逐步往后延后一个时间点，即1-24,2-25，3-26逐步往后延后一个时间点，避免其他变量的影响，将第一个固定时间段上的散点用最小二乘法拟合成一条直线并求出解析式:y=ax b；b、令时间点2即x
2，
令x=x2代入解析式y=ax b，求解y
2.，
得到点（x2,y2）；同时，用最小二乘法求出第二个时间段（2-25）上的散点拟合直线解析式y0=a0x b0以及该解析式的中点（x
中
，y
中
），x
中
即为第13个时间点（x
中
=x
13
）；c、将点（x2，y2）与点（x
中
，y
中
）连线，并求得该直线解析式y1=a1x b1，在解析式y1=a1x1 b1的基础上求出该直线上的第二个点的坐标（x3，y3），同时，用最小二乘法求出第三个时间段（3-26）上的散点拟合直线解析式y2=a2x b2以及该解析式的中点（x
中1
，y
中1
）；d、将点（x3，y3）与点（x
中1
，y
中1
）连线，并求得该直线解析式y3=a3x3 b3；e、重复cd操作。
[0010]
本发明方法的具体应用如下，可参见图1、图2：1、党川4#滑坡位于甘肃省永靖县盐锅峡镇黑方台西南侧黄河边，滑源区长300m，
宽20m，总滑坡体积约34万m3。自2017年8月底开始，滑坡开始产生变形且变形速率逐渐加快。2017年10月1日凌晨5时许，党川4#滑坡附近连续产生滑动，在滑源区形成了3个凹槽，并在滑坡下方形成了超过300m长的堆积体。其中，监测获取的突发型滑坡变形数据并非光滑曲线，会影响临滑预报的准确性将会出现多个误差点，从而导致多次误报的发生，将会在2017年6月3日、2017年7月5日、2017年8月7日、2017年8月20日等（如图1）多个点产生误报情况，而实际滑坡发生时间2017年10月1日凌晨5点，具有较大的误差预报情况。若采用本发明过滤方法构建可靠的实时变形数据过滤方法，对获取的监测数据进行实时过滤，很好地平衡误差剔除和及时识别加速变形（如图2）。
[0011]
综上，本发明针对现有数据过滤方法不能体现实时性、预报时间滞后两方面问题，通过本发明过滤方法结合常规移动平均法、卡尔曼滤波、样条平滑等常用数据过滤方法对监测数据进行实时过滤，并平衡和剔除误差数据，构建可靠的实时变形数据过滤方法，保证数据的实时更新和规避数据误差。本发明方法操作简单，避免了无法实时过滤数据和采用计算区域不同导致同一预报方法得到的结果有差异的问题，为滑坡应急处置提供可靠的时间参考。
[0012]
本发明的实施例公布的是较佳的实施例，但并不局限于此，本领域的普通技术人员，极易根据上述实施例，领会本发明的精神，并做出不同的引申和变化，但只要不脱离本发明的精神，都在本发明的保护范围内。