一种基于数据压缩变换的矿区边坡预测方法技术

本发明专利技术公开了一种基于数据压缩变换的矿区边坡预测方法，利用以往的滑坡累积位移观测数据，利用变分模态分解的方法进行滑坡数据本身的分解，期望通过分解的方法赋予每组分解数据独特的物理意义，以方便长短时记忆神经网络的预测。再通过分析每组数据的波动程度，针对波动较大的数据进行沃尔什—阿达玛变换，将改进后的数据进行长短时记忆神经网络的预测，再经过沃尔什—阿达玛逆变换将预测后的数据还原，最终将各组的预测值进行累积相加得出滑坡累积位移。该方法可在有限的观测值下，得出较高的预测精度。高的预测精度。高的预测精度。

【技术实现步骤摘要】
一种基于数据压缩变换的矿区边坡预测方法


[0001]本专利技术涉及露天矿边坡监测领域，具体为一种基于数据压缩变换的矿区边坡预测方法。

技术介绍

[0002]依据国家统计局的数据，我国近五年的滑坡灾害，约占地质灾害的50％以上，因此滑坡是我国主要的地质灾害，每年造成的大量人员伤亡和难以统计的财产损失。而边坡累积位移预测是滑坡灾害的预测的关键部分，因此如果能够露天矿边坡监测运用到生产实践中，将极大的提高滑坡灾害的预测能力。
[0003]现阶段的边坡累积位移预测模型分为四种，分别是确定性模型、统计模型、非线性模型和组合模型。近几年来，随着研究的深入，单一模型已经不能满足实际生产的需要，越来越多的学者，通过深度学习与其它学科相结合，来获取滑坡累积位移预测值。同时现有的滑坡位移预测模型往往需要很多的观测数据和考虑很多的滑坡影响因子，例如：像地质的岩层，断层，节理，高程，坡度，降水量，归一化植被指数，滑坡的水平累积位移等。上述的这些数据有时往往在生产实践中很难获取，从而导致预测模型构建的困难。
[0004]因此本文选取了在实际生产中最易获得的数据，利用滑坡的水平累积位移作为基础来构建模型，提出一种WHT
‑
VMD
‑
LSTM组合模型来降低上述各种模型对数据种类和数量的要求，选取了在实际生产中最易获得的数据，利用滑坡的水平累积位移作为基础来构建WHT
‑
VMD
‑
LSTM组合模型。

技术实现思路

[0005]本专利技术目的是为了提供一种在数据量较少时，滑坡累积位移预测的方法和提高预测精度。
[0006]为实现上述目的，本专利技术提供了一种基于数据压缩变换的矿区边坡预测方法，包括以下步骤：
[0007]建立WHT
‑
VMD
‑
LSTM新模型，对滑坡进行定期监测获取原始位移观测数据；
[0008]根据VMD分解后的平均残差，确定出分解层数；
[0009]根据VMD分解的原理确定出其带宽设置、对噪声的容忍程度、第一个模态是不是直流分量、中心频率和收敛精度5个参数；
[0010]根据确定出VMD分解的参数，然后进行原始观测数据的分解；
[0011]对分解后的数据画出趋势图，分析其波动的频率；
[0012]通过观察数据趋势图单调性改变的次数，如果大于改变次数大于10则波动较大，该组需要进行沃尔什—阿达玛变换压缩的分组；
[0013]根据分组后的数据分别利用长短时记忆神经网络进行预测；
[0014]将其中经过沃尔什—阿达玛变换压缩的预测值进行反变换，最后将各组的预测值相加，得出滑坡累积位移预测值。
[0015]本专利技术的核心在于利用信号分解与重构和图像压缩的思想，对原始累积滑坡位移观测数据进行处理，使其能够更好的利用长短时记忆神经网络进行预测滑坡累积位移值，解决了长短时记忆神经网络在预测时需要大量训练值的问题，有利于当观测数据不足时，进行滑坡累加位移的预测。
附图说明
[0016]图1是本专利技术的预测滑坡累积位移流程图；
[0017]图2是本专利技术的长短时记忆神经网络的结构图；
[0018]图3是经过VMD分解后的各分量的第一组趋势图；
[0019]图4是经过VMD分解后的各分量的第二组趋势图；
[0020]图5是经过VMD分解后的各分量的第三组趋势图；
[0021]图6是经过VMD分解后的各分量的第四组趋势图；
[0022]图7是经过VMD分解后的各分量的第五组趋势图；
[0023]图8是经过VMD分解后的各分量的第六组趋势图；
[0024]图9是经过沃尔什—阿达玛变换的第三组与原始个分量的趋势比较图；
[0025]图10是经过沃尔什—阿达玛变换的第四组与原始个分量的趋势比较图；
[0026]图11是经过沃尔什—阿达玛变换的第五组与原始个分量的趋势比较图；
[0027]图12是经过沃尔什—阿达玛变换的第六组与原始个分量的趋势比较图。
具体实施方式
[0028]下面详细通过实例介绍一下本专利技术的实现过程，通过实例介绍本专利技术是为了有助于理解本专利技术，而不是局限于该实例。
[0029]本专利技术的实验原始数据来源于实地的观测数据。
[0030]本专利技术的实验流程请参考图一，具体有以下步骤：
[0031]建立WHT
‑
VMD
‑
LSTM新模型，对滑坡进行定期监测获取原始位移观测数据；
[0032]通过多次VMD分解结果累加值与原始数据进行平均残差与最大残差的比较，从表2的结果来看当分解层数为6时其平均残差值最小，且其最大残差也在1mm内，因此本实验取VMD的分解层数为6。
[0033]最终将VMD分解的参数设为：alpha带宽设置为104，tau对噪声的容忍程度为0，k需要分解的层数为6，DC第一个模态是不是直流分量设为0，init设为0，即每个分量的中心频率为0，tol收敛的精度为设为1e
‑
10。
[0034]输入VMD的参数带入到VMD处理代码中即可得出分解后的六个分量如图3所示，其VMD的基本原理公式如下所示：
[0035][0036][0037][0038]其中：k为模态分解个数、{u
k
}为第k个模态分量、{w
k
}为模态分量的中心频率、f(t)为狄拉克函数、*为卷积运算、r为拉格朗日乘法算子、a为二次惩罚因子以保证信号重构的精度、h(t)用来保持约束条件的严格性。
[0039]表1原始累计位移数据
[0040][0041]表2 VMD分解结果
[0042]VMD分解层数三层四层五层六层平均残差值0.0170680.0167890.0000790.000077最大残差值2.1794811.2700001.2605320.641517。
[0043]观察图3
‑
图8中的六组分量趋势图，其经过VMD分解后使其每个分组具有其独特的意义，第一二组是趋势分量，第三四组是周期分量，第五六组是随机分量。除了第一二组分量外，其余组的分量波动频率都高，趋势复杂，不利于后期长短时记忆神经网络的学习，因此需要对第三
‑
六组的分量使用图像压缩的一种方法沃尔什—阿达玛变换压缩变量，使其波动频率减少，有利于后期长短时记忆神经网络的学习。
[0044]利用沃尔什—阿达玛变换压缩后的各组分量的图像如图9
‑
图12所示，从图中明显的可以看出经过压缩后的图像其数据变小了，趋势变缓了，比原来的数据更有利于长短时记忆神经网络的预测，其主要利用的公式为：
[0045][0046]其中WAL(n,i)是矩阵变换核。
[0047]然后将各组分量分别进行预测，将第一，二组的分量数据的前90％为训练集，后10％为测试集，由于沃尔什—阿达玛变换的特性，对变换数据本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于数据压缩变换的矿区边坡预测方法，其特征在于包括以下步骤：建立WHT
‑
VMD
‑
LSTM新模型，对滑坡进行定期监测获取原始位移观测数据；根据VMD分解后的平均残差，确定出分解层数；根据VMD分解的原理确定出其带宽设置、对噪声的容忍程度、第一个模态是不是直流分量、中心频率和收敛精度5个参数；根据确定出V...

【专利技术属性】
技术研发人员：肖玉强，李东印，赵鑫，
申请(专利权)人：中铁十九局集团矿业投资有限公司，
类型：发明
国别省市：