泥石流沟道中山洪泥石流灾害监测方法与流程

本发明涉及一种监测方法，特别是涉及一种在泥石流沟道中利用环境噪声信号对山洪泥石流灾害的运动加以监测的方法。属于信号装置技术领域、山区环境灾害防治领域。

背景技术：

山洪泥石流灾害形成主要是由环境条件变化的激发所致。其中最重要的三个环境条件是特定地形地貌条件、丰富的松散土石碎屑固体物质来源、强降雨。此三者也是山区灾害预报研究中用于判断山洪泥石流是否发生的主要标准。从环境噪声的研究角度看，激发山洪泥石流的外部自然活动(如降水、强风、土石运动等等)都能够从环境噪声值的变化中得到体现。并且，环境噪声值的变化特征，还能够反映多山洪泥石流的诸多有研究意义的运动特征。同时，环境噪声数据不受环境光线条件、能见度等条件影响，能够在夜间监测中发挥明显优势。

目前，山洪泥石流灾害的监测主要有两类方式：第一类方式是布置超声波物位计、雷达测速仪对通过断面的泥石流进行观测，只能检测到泥位、流速(如泥石流规模较小，不走传感器下面过，无法测量)，并辅以人工观测记录。这类方式获得的监测数据较粗糙，准确性不高。而且在灾害危险较重的时段与区域，也不无实施人工观测。二是在超声波物位计、雷达测速仪的基础上增加图像监控制备，利用计算机图形分析技术提高对山洪泥石流发生运动的监测科学性。这类方式更为科学准确，获得的监测数据质量好。但为数较多的山洪泥石流发生在夜晚，光线、大降雨、起雾等自然因素会严重影响图像监控制备的成像质量，降低后期分析的有效性。并且这类监测系统往往造价高，不便于在山区大规模建设。

技术实现要素：

本发明的目的就是针对现有技术的不足，提供一种利用泥石流沟道中的环境噪声信号对沟道中的山洪泥石流灾害的发生、运动实施监测的方法。

为实现上述目的，本发明的技术方案如下：

一种泥石流沟道中山洪泥石流监测方法，其特征在于：利用环境噪声分析实现对泥石流沟道中山洪泥石流灾害的监测，依如下步骤实施：

步骤s1、选定泥石流沟道作为监测区域，设置泥石流沟道监测断面，在监测断面上的泥石流沟道岸旁安装环境噪声传感器，环境噪声传感器与上位机通信连接；在上位机中设置初始背景噪声值lbackground、时间阈值tnsl、泥石流判定阈值lnsl、山洪判定阈值lsh，开始采集分析数据；

步骤s2、山洪泥石流灾害监测预报

上位机获取环境噪声传感器采集的环境噪声信号l样本，进行分析：

步骤s21、确定校准后lbackground

自数据采集时间起，将环境噪声信号样本l与lbackground对比，采用l校准lbackground：若lbackground≤l，上位机记录lbackground，若lbackground＜l，上位机将l记录为lbackground；上位机根据至少48h无风无雨环境中的l校准lbackground，得到校准后lbackground用于后续运算；

步骤s22、滤波

滤波开始，标记当前时刻为t，

自t时刻起，取t至t－5s间环境噪声信号l的最小值记为第1个lmin，，再沿时间轴取t+1s至t－4s间环境噪声信号的最小值记为第2个lmin，以此类推得到滤波lmin；

自t时刻起，计算t至t－30s间平均lmin，记为第1个laverage，再沿时间轴取t+1s至t－29s间平均lmin，记为第2个laverage，得到滤波laverage；

自t时刻起，计算t至t－300s间平均lmin，记为第1个lthreshold，再沿时间轴取t+1s至t－299s间平均lmin，记为第2个lthreshold，得到滤波lthreshold；

步骤s23、监测泥石流沟道内自然活动变化趋势

自t时刻起沿t+1时刻方向实时分析环境噪声信号l，当laverage>lthreshold时，记录两波交叉时刻为tstart，判断自tstart时刻起沟道内自然活动趋势增大，执行步骤s24；若laverage≤lthreshold时，判断沟道内自然活动变化不大或趋势减小，继续监测直到满足条件确定tstart；

步骤s24、监测自然活动峰值到达监测断面

自tstart时刻起，将laverage波峰值记为lmax、波峰时刻记为tmax，且判定有强自然活动在时刻tmax到达监测断面，发出监测信号；

步骤s25、判别山洪泥石流灾害类别

若tmax≥tstart+tnsl，则进一步判别：

若lmax-lbackground≥lnsl，判断在tmax时刻到达监测断面的是泥石流；若lnsl>lmax-lbackground≥lsh，判断在tmax时刻到达监测断面的是山洪；若lsh>lmax-lbackground，判断在tmax时刻到达监测断面的是常流水；

若tmax＜tstart+tnsl，放弃当前分析，继续步骤s23。

本发明前期研究发现，山洪泥石流在泥石流沟道中行进时形成的环境噪声数据波形在噪声值大小、衰减特征方面都具体独特性，能够使其自身与其它环境噪声相区分。例如山洪泥石流沿沟道从上游流来噪声越来越大，然后向下游流走噪声越来越小，历时几十上百秒。其噪声中最突出的时间变化特征与背景中的风雨雷电等自然气象形成的环境噪声有明显不同。同时，不同类型的山洪泥石流形成的环境噪声信号波形又具有一定程度上的相异性。例如阵性山洪泥石流形成的环境噪声波形峰值突出，易于识别，且阵与阵之间存在明显的声波中断，并且中断时间存在规律性。基于这些发现，本发明上述山洪泥石流监测方法在泥石流沟道中设立监测断面，在断面上安装环境噪声传感器，采集环境噪声信号样本l。通过滤波分析有效地提取信号中的噪声值与时间特征信息，在解读环境噪声信号样本l沿时间轴的多项衰减、振荡、激发特征后，辨识出山洪泥石流的主要运动特征及其到达监测断面的时间。

上述方法中，上位机中预置的各指标值一般可根据观测经验确定。

具体可采用：

初始背景噪声值lbackground：背景噪声值lbackground是取监测沟道中环境噪声的最低值，因而，在上位机中预置的初始值一般根据经验或前期观测的初步分析结果确定，选择一个比监测沟道中通常环境噪声值略高的值，然后再在监测过程中对其加以替换校准。

时间阈值tnsl：时间阈值tnsl是估测泥石流能够被传感器“听”到后到达监测断面所需的时间。一般沟道中泥石流在300m外产生的噪声已经很大，只要传感器无遮蔽，均可以检测到。如果以泥石流流速10m/s估算，30s会到达监测断面，所以可预置为30s。

泥石流判定阈值lnsl与山洪判定阈值lsh可根据现有观测资料选定。本发明根据在蒋家沟研究地的长期观测数据(常流水时环境噪声值在65db～75db之间，山洪发生时环境噪声在75db～90db之间，泥石流发生时环境噪声在90db以上)，提供参考值：泥石流判定阈值lnsl预置为30db、山洪判定阈值lsh预置为15db。

上述山洪泥石流监测方法在监测装置布设时的优化可以包括：一、将环境噪声信号频率设置为不低于每秒生成2个环境噪声值，即信号频率不低于2hz。二、将监测断面设置在泥石流沟道转弯处，以有效利用山洪泥石流转弯时形成的环境噪声更突出的特征，降低监测结果误判可能性。三、环境噪声传感器安装位置在距离沟道边线30m～100m范围内。该距离能够有效避免距离沟边岸太近造成传感器被冲毁，并且能够有效防止由于沟岸边噪声较大，导致山洪泥石流发出的噪声衰减和沟边背景噪声差值较小，以及距离太远造成环境噪声整体衰减特征突出，掩盖了某些细微的信号变化致使小规模山洪泥石流无法监测。

在上述山洪泥石流监测方法的基础上，可以增设监测断面，通过分析沿泥石流沟道布设的上、下二监测断面上采集到的环境噪声信号，能够进一步监测山洪泥石流灾害的两个运动参数：其一是，山洪泥石流在两监测断面间的位置；其二是山洪泥石流的流速。

与现有技术相比，本发明的有益效果是：本发明提供了一种全新的利用环境噪声信号监测泥石流沟道中山洪与泥石流运动的方法。通过将监测断面设置在保护地上游，本发明方法可以用于泥石流沟道区域的山洪泥石流灾害预警。由于环境噪声的监测不受环境光线条件、能见度等因素影响，因而本发明方法能够很好地解决夜间突发山洪泥石流时，基于图像分析的泥石流沟通监测技术应用受限的问题。并且，由于山洪泥石流在夜间发生的危害性远高于白天发生的危害性，因为本发明方法也具有更突出的山区减灾防灾效果。再加之环境噪声传感器价格相对低廉，信号传输成本小，因而通过在泥石流沟道沿岸更密集地布设监测断面，本发明可以实现对山洪泥石流的全程运动监测，确定开始结束时间、计算山洪泥石流的多种运动参数，从而可以为山洪泥石流的野外研究提供更多的山洪泥石流运动要素数据。

附图说明

图1泥石流沟道内监测断面布置示意图。

图2是原始数据、laverage、lthreshold波形图(la部分数据)。

具体实施方式

下面结合附图，对本发明的优选实施例作进一步的描述。

实施例一

如图1所示，用本发明方法对某泥石流沟道中山洪泥石流实施监测。

步骤s1、监测装置布置

图1泥石流沟道内监测断面布置示意图。图中横实线示泥石流沟道，以泥石流沟道作为监测区域，沿泥石流沟道方向设置上下二监测断面，分别记为监测断面a、监测断面b，二监测断面相距sab。在二监测断面上分别安装环境噪声传感器a与环境噪声传感器b。二环境噪声传感器可以安装在泥石流沟道同岸，也可以分别安装在泥石流沟道两岸。环境噪声传感器a、b均与上位机通信连接。

sab一般设置为100m～300m，以避免二断面距离太近差值太小不易区分，或者距离太远衰减太大误差增大。本实施方式中sab＝180m。环境噪声传感器a与泥石流沟道上游岸边垂直距离ha＝260m、环境噪声传感器b与泥石流沟道下游岸边垂直距离hb＝30m。环境噪声信号频率均为2hz。

根据前期观测经验值，在上位机中分别设置如下阈值：

传感器a：初始背景噪声值labackgroud＝100db、时间阈值tansl＝30s、泥石流判定阈值lansl＝10db、山洪判定阈值lash＝5db；

传感器b：初始背景噪声值lbbackground＝100db、时间阈值tbnsl＝30s、泥石流判定阈值lbnsl＝30db、山洪判定阈值lbsh＝15db。

设置完成，开始采集数据。

步骤s2、山洪泥石流灾害监测预报

上位机分别获取环境噪声传感器a、b采集到的环境噪声信号样本la、lb，分别进行分析：

步骤s21、校准背景噪声lbackground

自数据采集时间起，将环境噪声信号样本la与初始背景噪声labackground相对比，采用la校准labackground，具体是：若labackground≤la，则上位机记录labackground；若labackground＞la，则上位机将la记录为labackground。上位机根据至少48h无风无雨环境中的la校准labackground，得到校准后labackground，用于后续运算。校准期间若遇风雨天气则延长校准时间，直到满足上位机记录无风无雨48h的labackground值为止。类似地，初始lbbackground也经环境噪声信号样本lb校准，得到校准后lbbackground，用于后续运算。

得到校准后labackground＝44.5db。环境噪声信号样本lb分析过程略。

步骤s22、滤波

滤波开始，标记当前时刻为t，自t时刻起，取t至t－5s间环境噪声信号l的最小值记为第1个lmin，再沿时间轴取t+1s至t－4s间环境噪声信号的最小值记为第2个lmin，以此类推。每采集一个环境噪声信号l，就会得到一个lmin，形成先进先出的min堆栈。min堆栈可存储5秒10笔。

例如，由于传感器以2hz采样，所以每秒都有2个数据。若t＝5s，则回溯取第5s至第t－5s＝0s间(含第1秒、2秒、3秒、4秒、5秒)10个数据中的最小环境噪声信号l记为第1个lmin(即5s时刻的lmin)，然后沿时间轴取t+1s＝6s至t－4s＝1s间(含第2秒、3秒、4秒、5秒、6秒)10个数据中的最小环境噪声信号l记为第2个lmin(即6s时刻的lmin)，再然后沿时间轴取t+2s＝7s至t－3s＝2s间(含第3秒、4秒、5秒、6秒、7秒)10个数据中的最小环境噪声信号l记为第3个lmin(即7s时刻的lmin)。由此得到滤波lmin。

自t时刻起，计算t至t－30s间平均lmin，记为第1个laverage，再沿时间轴取t+1s至t－29s间平均lmin，记为第2个laverage。即沿时间轴形成先进先出的average堆栈。average堆栈可存储30秒60笔。

由此得到滤波laverage。

自t时刻起，计算t至t－300s间平均lmin，记为第1个lthreshold，再沿时间轴取t+1s至t－299s间平均lmin，记为第2个lthreshold。即形成先进先出的threshold堆栈。threshold堆栈可存储300秒600笔。由此得到滤波lthreshold。

上位机处理时，每采集到一个l，即进入min堆栈，同时min堆栈最先进入的l自动出栈，取min堆栈内最小值即lmin。得到lmin将lmin同时压入average堆栈与threshold堆栈。分别计算average堆栈与threshold堆栈内数据的平均值即为laverage、lthreshold。图2是原始数据、laverage、lthreshold波形图(la部分数据)

步骤s23、监测泥石流沟道内自然活动变化趋势

自t时刻起沿t+1时刻方向实时分析环境噪声信号l，当laverage>lthreshold时，记录两波交叉时刻为tstart，判断自tstart起，沟道内自然活动趋势增大，即山洪泥石流灾害可能性增加，执行步骤s24。若laverage≤lthreshold时，判断自tstart起，沟道内自然活动变化不大或趋势减小，继续监测直到满足条件确定tstart。

以环境噪声信号样本la分析为例(图2)，当第1次laverage>lthreshold时记为tastart1＝00:38:16.0；第2次laverage>lthreshold时记为tastart2＝00:44.46.0。这2个时刻沟内自然活动趋势出现增大趋势，应该注意监测沟道内环境噪声信号样本变化。

步骤s24、监测自然活动峰值到达监测断面

自tstart时刻起，将laverage波峰值记为lmax、波峰时刻记为tmax，且判定有强自然活动(即疑似山洪泥石流)在时刻tmax到达监测断面，发出监测信号。

以环境噪声信号样本la分析为例(图2)，有自tastart1起，laverage波峰值记为lamax1＝66.3db，波峰时刻记为tamax1＝00:39:00.0；自tastart2起，laverage波峰值记为lamax2＝71.8db，波峰时刻记为tamax2＝00:45:34.0。

对于样本lb，监测计算确定tbmax1＝00:39:25.5。

步骤s25、判别山洪泥石流灾害类别

若tmax≥tstart+tnsl，则进一步判别：

若lmax-lbackground≥lnsl，判断在lmax时刻到达监测断面的是泥石流；若lnsl>lmax-lbackground≥lsh，判断在tmax时刻到达监测断面的是山洪；若lsh>lmax-lbackground，判断在tmax时刻到达监测断面的是常流水。根据判断结果，发出对应的监测信号，沿时间轴继续监测。

若tmax＜tstart+tnsl，放弃当前分析，继续步骤s23。

以环境噪声信号样本la分析为例(图2)，tastart1＝00:38:16.0、tamax1＝00:39:00.0、lamax1＝66.3db、labackground＝44.5db、tansl＝30s，同时满足tmax≥tstart+tnsl和lmax-lbackground≥lnsl两项条件，表明有一阵泥石流在tamax1＝00:39:00.0通过监测断面a；tastart2＝00:44:46.0、tamax2＝00:45:34.0、lamax2＝71.8db、labackground＝44.5db、tansl＝30s，同时满足tmax≥tstart+tnsl和lmax-lbackground≥lnsl两条件，表明有另一阵泥石流在tamax2＝00:45:34.0通过监测断面a。对应发出泥石流监测信号。

实际结果：根据当日人工观测泥石流记录，在00:39和00:45各有一阵泥石流流过监测断面a。

步骤s3、计算泥石流在二监测断面间的位置

在tastart时该后，但尚未监测到泥石流到达监测断面b之前的任一时间t，沟道中泥石流处于监测断面a与监测断面b之间。计算此时泥石流在监测断面a、b间的位置。

在信号样本la、lb中，分别提取t＝00:39:10.0时的环境噪声信号la＝59.2db、lb＝66.7db；

将la、lb、sab＝180m、ha＝260m、hb＝30m代入式1，计算泥石流距离监测断面b的距离sbc：

计算有sbc≈109.56m。

步骤s4、计算沟道泥石流的流速

将步骤s2确定的泥石流到达监测断面a的时间tamax1＝00:39:00.0、泥石流到达监测断面b的时间tbmax1＝00:39:25.5，以及步骤s1中sab＝180m代入式2：

计算有泥石流的流速vab≈7.06m·s^-1。