一种基于形变有序性的滑坡灾害多监测点融合预警方法与流程
1.本发明涉及地质灾害监测技术领域，具体涉及一种基于形变有序性的滑坡灾害多监测点融合预警方法。
背景技术：
2.滑坡的监测预警是滑坡灾害由“事后救助”向“事前预防”转变的关键环节，如何提升滑坡预警预报的可靠性是滑坡监测预警工作的重点和难点。当前主流的监测预警主要以滑坡位移作为主要参数，在预警过程中则主要采用单个监测点的形变时间曲线进行预警。由于单一曲线往往仅能代表监测点局部附近很小范围内的变形情况，而无法反应整个坡体的变形特性，从而往往出现误报、漏报的现象。因此，如何综合利用同一滑坡体上不同监测点的形变特征，做到多监测点的融合利用，从滑坡整体角度进行预警预报，从而避免现有预警预报方法的疏漏和不足，是新形势下滑坡监测预警需要解决的紧迫问题。
技术实现要素：
3.本发明意在提供一种滑坡灾害多监测点融合预警方法，以解决单个监测点的形变曲线进行预警的疏漏和不足。
4.为达到上述目的，本发明采用如下技术方案：一种基于形变有序性的滑坡灾害多监测点融合预警方法，包括以下步骤：
5.s1：获取滑坡多个监测点的监测传感器预设时刻的形变速率；
6.s2：根据形变速率计算获得形变速率信息熵hv(t)；
7.s3：根据地形地貌特征获取滑坡变形主方向，即目标形变方向；
8.s4：获取各监测点形变方向，即监测形变方向；
9.s5：计算监测形变方向和目标形变方向的形变方向关联熵hd(t)；
10.s6：根据形变速率信息熵hv(t)的状态和形变方向关联熵hd(t)的状态，获取滑坡预警等级。
11.优选的，作为一种改进，所述根据形变速率计算获取形变速率信息熵的具体步骤包括：根据形变速率的大小进行分组，并获取形变速率组间距l；获取形变速率组内每个形变速率在当前形变速率组内出现的概率pi(t)；计算多个监测传感器的预设时刻形变速率的形变速率信息熵。
12.优选的，作为一种改进，所述形变速率组间距l为：
[0013][0014]
其中，x
max
，x
min
分别为监测点形变速率的最大和最小值，k为形变速率分组数。
[0015]
优选的，作为一种改进，所述形变速率组内每个形变速率在当前形变速率组内出现的概率pi(t)为：
[0016]
[0017]
其中，ni(t)为预设时刻形变速率落在第i个分组区间内的监测传感器数量，n为监测传感器总数。
[0018]
优选的，作为一种改进，所述形变速率信息熵hv(t)具体为：
[0019][0020]
其中，i表示第i个形变速率的分组。
[0021]
优选的，作为一种改进，所述根据形变速率的大小进行分组具体包括：根据监测传感器数量确定分组数，根据分组数和形变速率的大小对形变速率进行分组。
[0022]
优选的，作为一种改进，所述形变方向关联熵hd(t)具体为：
[0023][0024]
其中，i表示第i个监测点，pi为关联系数概率，n为监测传感器总数。
[0025]
优选的，作为一种改进，所述计算监测形变方向和目标形变方向的形变方向关联熵hd(t)具体包括：根据方位角表达监测形变方向序列和目标形变方向序列；对监测形变方向序列和目标形变方向序列进行归一化处理；对归一化处理后的监测形变方向序列和目标形变方向序列进行灰关联度计算；就根据灰关联度计算关联系数概率分布；根据关联系数概率分布计算形变方向关联熵hd(t)。
[0026]
优选的，作为一种改进，所述预警等级包括：蓝预警、黄预警、橙预警和红预警。
[0027]
优选的，作为一种改进，所述形变速率信息熵hv(t)状态包括sh
v-1、sh
v-2和sh
v-3，其中，所述sh
v-1表示形变速率信息熵hv(t)随时间的变化平稳发展或围绕中轴线上下波动，无明显的增加或降低趋势，所述sh
v-2表示形变速率信息熵hv(t)随时间的变化具有明显的下降趋势，所述sh
v-3表示形变速率信息熵hv(t)随时间的变化持续保持下降趋势，且逐步逼近于零；
[0028]
所述形变方向关联熵hd(t)状态sh
d-1、sh
d-2、sh
d-3表示，其中，所述sh
d-1表示形变方向关联熵hd(t)随时间的变化平稳发展或围绕中轴线上下波动，无明显的增加或降低趋势，所述sh
d-2表示形变方向关联熵hd(t)随时间的变化具有明显的上升趋势，所述sh
d-3表示形变方向关联熵hd(t)随时间的变化持续保持上升趋势，且逐步逼近于-ln(1/n)，其中n为为同一滑坡体上的监测传感器数量；
[0029]
所述形变速率信息熵hv(t)状态为sh
v-1和形变方向关联熵hd(t)状态为sh
d-1时，所述预警等级为蓝预警；所述形变速率信息熵hv(t)状态为sh
v-1和形变方向关联熵hd(t)状态为sh
d-2时或所述形变速率信息熵hv(t)状态为sh
v-2和形变方向关联熵hd(t)状态为sh
d-1时，所述预警等级为黄预警；所述形变速率信息熵hv(t)状态为sh
v-1和形变方向关联熵hd(t)状态为sh
d-3时、所述形变速率信息熵hv(t)状态为sh
v-2和形变方向关联熵hd(t)状态为sh
d-2时或所述形变速率信息熵hv(t)状态为sh
v-3和形变方向关联熵hd(t)状态为sh
d-1时，所述预警等级为橙预警；所述形变速率信息熵hv(t)状态为sh
v-2和形变方向关联熵hd(t)状态为sh
d-3时、所述形变速率信息熵hv(t)状态为sh
v-3和形变方向关联熵hd(t)状态为sh
d-2时或所述形变速率信息熵hv(t)状态为sh
v-3和形变方向关联熵hd(t)状态为sh
d-3时，所述预警等级为红预警。
[0030]
本方案利用同一滑坡体上不同监测点的形变速率和形变方向特征，做到多监测点的融合利用，从滑坡的整体角度进行预警预报，从而避免现有预警预报方法的疏漏和不足；滑坡在进入加速变形阶段时对地表起控制作用的是滑坡的整体运动，局部微地貌控制作用明显减弱，此时地表的变形趋于有序，根据有序性演化过程可以对滑坡变形破坏阶段进行定量描述，从而可以根据有序性进行滑坡预警，本方案利用监测传感器的获取滑坡的形变速率，通过形变速率有序性演化过程来对滑坡变形阶段进行定量描述，同时选取确认主滑方向作为目标方向，根据监测点监测得到的形变方向与主滑方向的关联程度来表示形变方向有序性，再结合形变速率有序性和形变方向有序性，建立预警依据，与现有的预警模型形成互补，提升灾害预警预报的可靠性和准确性；利用滑坡形变的两个重要参数：形变速率和形变方向两个参数，参数较少，能够描述滑坡变形破坏阶段且减少预警误差。
附图说明
[0031]
图1为本发明实施例的形变速率求取示意图。
[0032]
图2为本发明实施例的形变方向计算示意图。
[0033]
图3为本发明实施例的示例滑坡的图3为形变速率信息熵随时序演化过程曲线。
[0034]
图4为本发明实施例的合位移方向计算示意图。
[0035]
图5为本发明实施例的形变速率信息熵hv(t)随时间的变化图。
[0036]
图6为本发明实施例的形变方向关联熵hd(t)随时间的变化图。
具体实施方式
[0037]
下面通过具体实施方式进一步详细说明：
[0038]
实施例：
[0039]
一种基于形变有序性的滑坡灾害多监测点融合预警方法，具体包括以下步骤：
[0040]
s1：获取滑坡多个监测点的监测传感器预设时刻的形变速率；
[0041]
完整的滑坡监测数据分为初始变形阶段、等速变形阶段、加速变形阶段。在滑坡变形破坏过程中，形变速率是表征滑坡稳定状态的一个重要参数，当某一滑坡表面所安装的位移传感设备达到一定的规模，符合统计条件，则可以从统计学角度对滑坡总体形变速率有序性进行分析。具体过程如下：
[0042]
提取某一时刻不同监测传感器的形变速率，具体过程如图1所示；
[0043]
首先在滑坡变形时序曲线上，按照固定的时间段(一般以天或小时为单位进行计算)，在时间轴上将时间位移曲线划分为若干小段，计算t时刻的形变速率时，提取(t
‑△
t,t)时间段内的形变增量
△
s，则t时刻形变速率为：
[0044][0045]
s2：利用形变速率，对形变速率大小进行分组，组数采用美国学者斯特吉斯提出的经验公式，即：
[0046]
k＝1+3.322
×
log(n)
ꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0047]
式中：k为形变速率分组数，n为同一滑坡体上的监测传感器数量。上述计算结果存在小数的，按照四舍五入规则处理。
[0048]
s21：根据形变速率组数，求组间距，即：
[0049][0050]
式中：l为形变速率组间距，x
max
，x
min
分别为监测点形变速率的最大和最小值；
[0051]
s22：根据速率分组结果，统计每一区间内，形变速率在该区间内的传感器数量，并计算该区间内的概率，即：
[0052][0053]
式中：pi(t)为t时刻第i个分组区间内形变速率出现的概率，ni(t)为t时刻形变速率落在第i个分组区间内的传感器数量，n为传感器总数。
[0054]
s23：引入信息熵理论，定量计算某一时刻形变速率有序性：
[0055][0056]
s24：对0～t整个时间段内任意一个时刻t分别计算有序性，则可获得形变速率有序性随时间演化过程曲线。
[0057]
以下是获得形变速率有序性随时间演化过程曲线其中一个示例，例如某一滑坡上设有20个形变监测点，以周为时间单位计算，统计从计算当周起前10周的形变速率情况，按照式(1)对每一天的形变速率进行计算，得到最大最小形变速率统计如下：
[0058][0059]
式中，t0表示当周，t1表示计算时间前1周，以此类推。
[0060]
(1)则根据上述信息可知n＝20，相应的按照式(2)可得到分组数k＝5.322，取k＝5。
[0061]
分别建立t0～t9共计10周内形变速率直方图，以t0为例进行说明。
[0062]
(2)首先根据式(3)及表1，计算t0时刻直方图组间距，l0：
[0063][0064]
(3)于是可以建立t0时刻的分组规则：
[0065]
[1.2～1.6)、[1.6～2.0)、[2.0～2.4)、2.4～2.8)、[2.8～3.2]
[0066]
(4)根据上述分组规则，按照式(4)统计行不速率出现在各个区间内的频率，p1(t0)、p2(t0)、p3(t0)、p4(t0)、p5(t0)，通过计算得到：p1(t0)＝0.15、p2(t0)＝0.14、p3(t0)＝
0.36、p4(t0)＝0.21、p5(t0)＝0.14。
[0067]
(5)根据式(5)计算t0时刻形变速率信息熵：
[0068][0069]
(6)依次对t0、t1、t2、t3、t4、t5、t6、t7、t8、t9时刻的形变速率信息熵进行计算得到形变速率有序性随时间演化曲线，具体如图3所示，图3为形变速率信息熵随时序演化过程曲线。
[0070]
s3：根据地形地貌特征获取滑坡变形主方向，即目标形变方向。
[0071]
对于典型地质灾害而言，在调查过程中通常会确定一个主滑方向，此主滑方向可以作为我们的目标方向。
[0072]
s4：获取各监测点形变方向，即监测形变方向；
[0073]
根据各监测点监测得到的形变方向与主滑方向的关联程度来表示形变方向有序性，显然，关联程度越大，形变方向系统有序性更强，根据灰理论，这一关联性可用灰关联系数来表示。
[0074]
s51：在进行有序性分析过程中，仅考虑平面方向，按照s1中形变速率求解时所确定的时间间隔，通过矢量计算，得到形变平面合方向，并以空间方位角表示。
[0075]
传感器上直接获取的形变增量包括x方向的形变增量δsx和y方向的形变增量δsy结合图4(合位移方向计算示意图)可以通过计算得到合位移方向角α：
[0076][0077]
当δs
x
》0且δsy≥0，则为第一象限，此时合位移方向角α＝α；
[0078]
当δs
x
《0且δsy≥0，则为第二象限，此时和位移方向角α＝180
°‑
α；
[0079]
当δs
x
《0且δsy《0，则为第三象限，此时和位移方向角α＝180
°
+α；
[0080]
当δs
x
》0且δsy《0，则为第四象限，此时和位移方向角α＝360
°‑
α；
[0081]
当δs
x
＝0且δsy》0，则α＝90
°
；
[0082]
当δs
x
＝0且δsy《0，则α＝270
°
；
[0083]
s52：形变方向序列表达：设x＝(x1,x2,
…
,xn)为某一时刻t灾害体上各监测点构成的形变方向序列，y＝(y1,y2,
…
,yn)为目标形变方向序列，此处的y1＝y2＝,
…
,yn＝y*，y*为滑坡的主滑方向。形变方向以方位角表示，由于方位角范围在0～360
°
之间，各监测点形变方向与主滑方向的差异不能完全按照普通数列来处理，具体如如图2所示。
[0084]
当方位角x未超过y*+180
°
时，x与y*的角度差可以直接将二者相减，而当x超过y*+180
°
时，x与y*的差值实际上应该为x-360
°
后再与y*作差。于是可将x形变序列划分为两个部分，即：
[0085][0086]
s53：形变方向归一化处理：对形变方向观测值序列x＝(x1,x2,
…
,xn)和目标方向
序列y＝(y1,y2,
…
,yn)作归一化处理，结果如下：
[0087][0088]
式中：xi为t时刻第i个监测点的形变方向值，x
min
，x
max
分别为t时刻监测点方向的最大最小值。显然，对于目标方向序列y，由于各监测点的目标方向均为滑坡主滑方向，因此其无量纲化后的结果为全为1,即：y
′i＝1。
[0089]
s54：灰关联度计算：t时刻序列x＇和y＇的灰关联系数表示为：
[0090][0091]
式中：δ(min)＝min(|x
′
i-y
′i|)，δ(max)＝max(|x
′
i-y
′i|)，0＜ρ＜1为分辨系数，通过调节该值可以改变关联系数的显著性，这里取ρ＝0.5。
[0092]
s55：计算关联系数概率分布：
[0093][0094]
s56：计算监测形变方向和目标形变方向的形变方向关联熵hd(t)，形变方向关联熵hd(t)具体为：
[0095][0096]
式中，pi(t)表示在t时刻的关联系数概率。
[0097]
s6：根据形变速率信息熵hv(t)的状态和形变方向关联熵hd(t)的状态，获取滑坡预警等级。
[0098]
公式(5)为形变速率有序性的最终表达式，从该表达式及其分析过程可知，速率的信息熵随着有序性的增强而减小，当所有的监测点形变速率趋于统一值时，形变速率信息熵趋于零，此时可以认为滑坡即将发生。因此，可以将形变速率趋于零作为形变速率有序性预警判据。
[0099]
对于形变方向有序性，采用式(9)计算其灰关联熵，显然，对于式(9)中所确定的灰关联熵，hd越大，形变方向的有序性越强，当所有形变方向都趋于一个方向时，ξi为常数，此时pi＝1/n，灰关联熵达到最大，其最大值为-ln(1/n)。因此可以通过实际获得的形变方向灰关联熵与最大灰关联熵-ln(1/n)的对比，来确定预警判据。
[0100]
滑坡在临近发生的时候，其形变的速率和方向均趋于一致，因此单一利形变速率和形变方向进行预警，往往容易造成误报和漏报，因此有必要将二者结合起来，通过二者的相互补充和应证，来提升滑坡预警预报的准确性。
[0101]
首先定义形变速率信息熵和形变方向关联熵的三种状态，其中形变速率信息熵状态sh
v-1、sh
v-2、sh
v-3表示，形变方向关联熵状态以sh
d-1、sh
d-2、sh
d-3表示，各状态具体定义如下：
[0102]
形变速率信息熵hv(t)随时间的变化具体如附图5(形变速率有序性等级划分示意图)所示。
[0103]
sh
v-1表示形变速率信息熵hv(t)随时间的变化平稳发展或围绕中轴线上下波动，无明显的增加或降低趋势；
[0104]
sh
v-2表示形变速率信息熵hv(t)随时间的变化具有明显的下降趋势(连续三个以上时间间隔内的变化速率基本一致或有所增加)；
[0105]
sh
v-3表示形变速率信息熵hv(t)随时间的变化持续保持下降趋势，且逐步逼近于零。
[0106]
形变方向关联熵hd(t)随时间的变化具体如附图6(形变方向有序性等级划分示意图)所示。
[0107]
sh
d-1表示形变方向关联熵hd(t)随时间的变化平稳发展或围绕中轴线上下波动，无明显的增加或降低趋势；
[0108]
sh
d-2表示形变方向关联熵hd(t)随时间的变化具有明显的上升趋势(连续三个以上时间间隔内的变化速率基本一致或有所增加)；
[0109]
sh
d-3表示形变方向关联熵hd(t)随时间的变化持续保持上升趋势，且逐步逼近于-ln(1/n)。
[0110]
基于上述两种熵值的状态定义，构建预警判据矩阵如下：
[0111][0112]
本方案利用同一滑坡体上不同监测点的形变速率和形变方向特征，做到多监测点的融合利用，从滑坡的整体角度进行预警预报，从而避免现有预警预报方法的疏漏和不足；滑坡在进入加速变形阶段时对地表起控制作用的是滑坡的整体运动，局部微地貌控制作用明显减弱，此时地表的变形趋于有序，根据有序性演化过程可以对滑坡变形破坏阶段进行定量描述，从而可以根据有序性进行滑坡预警，本方案利用监测传感器的获取滑坡的形变速率，通过形变速率有序性演化过程来对滑坡变形阶段进行定量描述，同时选取确认主滑方向作为目标方向，根据监测点监测得到的形变方向与主滑方向的关联程度来表示形变方向有序性，再结合形变速率有序性和形变方向有序性，建立预警依据，与现有的预警模型形成互补，提升灾害预警预报的可靠性和准确性；利用滑坡形变的两个重要参数：形变速率和形变方向两个参数，参数较少，能够描述滑坡变形破坏阶段且减少预警误差。
[0113]
以上所述的仅是本发明的实施例，方案中公知的具体技术方案和/或特性等常识在此未作过多描述。应当指出，对于本领域的技术人员来说，在不脱离本发明技术方案的前提下，还可以作出若干变形和改进，在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部
的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。本技术要求的保护范围应当以其权利要求的内容为准，说明书中的具体实施方式等记载可以用于解释权利要求的内容。
技术特征：
1.一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于，包括以下步骤：s1：获取滑坡多个监测点的监测传感器预设时刻的形变速率；s2：根据形变速率计算获得形变速率信息熵hv(t)；s3：根据地形地貌特征获取滑坡变形主方向，即目标形变方向；s4：获取各监测点形变方向，即监测形变方向；s5：计算监测形变方向和目标形变方向的形变方向关联熵h
d
(t)；s6：根据形变速率信息熵hv(t)的状态和形变方向关联熵h
d
(t)的状态，获取滑坡预警等级。2.根据权利要求1所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述根据形变速率计算获取形变速率信息熵的具体步骤包括：根据形变速率的大小进行分组，并获取形变速率组间距l；获取形变速率组内每个形变速率在当前形变速率组内出现的概率pi(t)；计算多个监测传感器的预设时刻形变速率的形变速率信息熵。3.根据权利要求2所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述形变速率组间距l为：其中，x
max
，x
min
分别为监测点形变速率的最大和最小值，k为形变速率分组数。4.根据权利要求2所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述形变速率组内每个形变速率在当前形变速率组内出现的概率pi(t)为：其中，n
i
(t)为预设时刻形变速率落在第i个分组区间内的监测传感器数量，n为监测传感器总数。5.根据权利要求4所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述形变速率信息熵hv(t)具体为：其中，i表示第i个形变速率的分组。6.根据权利要求2所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述根据形变速率的大小进行分组具体包括：根据监测传感器数量确定分组数，根据分组数和形变速率的大小对形变速率进行分组。7.根据权利要求1所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述形变方向关联熵h
d
(t)具体为：其中，i表示第i个监测点，p
i
为关联系数概率，n为监测传感器总数。8.根据权利要求1所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述计算监测形变方向和目标形变方向的形变方向关联熵h
d
(t)具体包括：根据
方位角表达监测形变方向序列和目标形变方向序列；对监测形变方向序列和目标形变方向序列进行归一化处理；对归一化处理后的监测形变方向序列和目标形变方向序列进行灰关联度计算；就根据灰关联度计算关联系数概率分布；根据关联系数概率分布计算形变方向关联熵h
d
(t)。9.根据权利要求1所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述预警等级包括：蓝预警、黄预警、橙预警和红预警。10.根据权利要求9所述的一种基于形变有序性的滑坡灾害多监测点融合预警方法，其特征在于：所述形变速率信息熵hv(t)状态包括sh
v-1、sh
v-2和sh
v-3，其中，所述sh
v-1表示形变速率信息熵hv(t)随时间的变化平稳发展或围绕中轴线上下波动，无明显的增加或降低趋势，所述sh
v-2表示形变速率信息熵hv(t)随时间的变化具有明显的下降趋势，所述sh
v-3表示形变速率信息熵hv(t)随时间的变化持续保持下降趋势，且逐步逼近于零；所述形变方向关联熵h
d
(t)状态sh
d-1、sh
d-2、sh
d-3表示，其中，所述sh
d-1表示形变方向关联熵h
d
(t)随时间的变化平稳发展或围绕中轴线上下波动，无明显的增加或降低趋势，所述sh
d-2表示形变方向关联熵h
d
(t)随时间的变化具有明显的上升趋势，所述sh
d-3表示形变方向关联熵h
d
(t)随时间的变化持续保持上升趋势，且逐步逼近于-ln(1/n)，其中n为为同一滑坡体上的监测传感器数量；所述形变速率信息熵hv(t)状态为sh
v-1和形变方向关联熵h
d
(t)状态为sh
d-1时，所述预警等级为蓝预警；所述形变速率信息熵hv(t)状态为sh
v-1和形变方向关联熵h
d
(t)状态为sh
d-2时或所述形变速率信息熵hv(t)状态为sh
v-2和形变方向关联熵h
d
(t)状态为sh
d-1时，所述预警等级为黄预警；所述形变速率信息熵hv(t)状态为sh
v-1和形变方向关联熵h
d
(t)状态为sh
d-3时、所述形变速率信息熵hv(t)状态为sh
v-2和形变方向关联熵h
d
(t)状态为sh
d-2时或所述形变速率信息熵hv(t)状态为sh
v-3和形变方向关联熵h
d
(t)状态为sh
d-1时，所述预警等级为橙预警；所述形变速率信息熵hv(t)状态为sh
v-2和形变方向关联熵h
d
(t)状态为sh
d-3时、所述形变速率信息熵hv(t)状态为sh
v-3和形变方向关联熵h
d
(t)状态为sh
d-2时或所述形变速率信息熵hv(t)状态为sh
v-3和形变方向关联熵h
d
(t)状态为sh
d-3时，所述预警等级为红预警。
技术总结
本发明涉及地质灾害监测技术领域，公开了一种基于形变有序性的滑坡灾害多监测点融合预警方法，包括以下步骤：获取滑坡多个监测点的监测传感器预设时刻的形变速率；根据形变速率计算获得形变速率信息熵HV(t)；计算监测形变方向和目标形变方向的形变方向关联熵H
