CN110927723B - 毫米波雷达泥石流智能监测预警系统与方法 - Google Patents

Abstract

本公开涉及一种毫米波雷达泥石流智能监测预警系统与方法，属于地质监测领域，能够精确、及时、成本低廉地对泥石流进行监测预警。一种毫米波雷达泥石流智能监测预警方法，包括：向被监测的泥石流区域发射第一微波信号；接收被所述泥石流区域反射回来的回波信号；对反射回来的回波信号进行下变频处理，得到回波中频信号；对所述回波中频信号进行采样；基于采样后的回波中频信号对被监测的泥石流区域进行成像，得到关于被监测的泥石流区域的图像；利用所述图像来监测被监测的泥石流区域的变化；基于监测到的变化进行泥石流预警。

Description

毫米波雷达泥石流智能监测预警系统与方法

技术领域

本公开涉及地质监测领域，具体地，涉及一种毫米波雷达泥石流智能监测预警系统与方法。

背景技术

目前，一种泥石流监测技术是采用接触式预警仪，然而，由于泥石流破坏性较大，接触式预警仪很容易在泥石流发生时损坏，从而影响监测预警精度和预警效果。另一种泥石流监测技术是采用非接触式预警仪。例如采用天气雷达监测泥石流隐患区域上空的降雨云团情况，但是该方法的造价高，而且由于其预警范围较大而不适用于小流域降雨情况。再例如采用超声波泥位计对平直、规则、不易冲毁的断面进行超声波监测并根据泥位大小实现对泥石流的监测与预警，但是该方法得到的信息是延时的，导致监测和预警不及时。

发明内容

本公开的目的是提供一种毫米波雷达泥石流智能监测预警系统与方法，能够精确、及时、成本低廉地对泥石流进行监测预警。

根据本公开的第一实施例，提供一种毫米波雷达泥石流智能监测预警方法，包括：向被监测的泥石流区域发射第一微波信号；接收被所述泥石流区域反射回来的回波信号；对反射回来的回波信号进行下变频处理，得到回波中频信号；对所述回波中频信号进行采样；基于采样后的回波中频信号对被监测的泥石流区域进行成像，得到关于被监测的泥石流区域的图像；利用所述图像来监测被监测的泥石流区域的变化；基于监测到的变化进行泥石流预警。

可选地，所述基于监测到的变化进行泥石流预警，包括：在监测到的变化指示所述图像中变化的像素数量超过预设阈值时，进行泥石流滑坡预警。

可选地，所述方法还包括：估计泥石流滑坡速度。

可选地，所述估计泥石流滑坡速度包括：在下变频处理得到的回波中频信号经历残余视频相位补偿处理之后，对残余视频相位补偿处理后的回波中频信号进行一维距离向压缩处理；在得到了Q个一维距离向压缩处理数据之后，将Q个一维距离向压缩处理数据组成Q行N列矩阵，其中N为距离向采集点数，Q为总的采集周期数；对所得到的矩阵信号组进行列向FFT变换；在变换得到的频域中寻找尖峰位置；计算所述尖峰位置处的滑坡速度值作为所估计的泥石流滑坡速度。

可选地，所估计的泥石流滑坡速度为：

其中，f₀为中心频率，Br为信号带宽，V为所估计的泥石流滑坡速度，C为电磁波速度，f_v为目标的多普勒频率值。

根据本公开的第二实施例，提供一种毫米波雷达泥石流智能监测预警系统，包括：阵列天线，用于向被监测的泥石流区域发射第一微波信号，并接收被所述泥石流区域反射回来的回波信号；MIMO模块，用于对反射回来的回波信号进行下变频处理，得到回波中频信号；采样模块，用于对所述回波中频信号进行采样；成像模块，用于基于采样后的回波中频信号对被监测的泥石流区域进行成像，得到关于被监测的泥石流区域的图像；变化监测模块，用于利用所述图像来监测被监测的泥石流区域的变化；预警模块，用于基于监测到的变化进行泥石流预警。

可选地，所述预警模块用于：在监测到的变化指示所述图像中变化的像素数量超过预设阈值时，进行泥石流滑坡预警。

可选地，所述预警模块还用于：估计泥石流滑坡速度。

可选地，所述预警模块还用于：在下变频处理得到的回波中频信号经历残余视频相位补偿处理之后，对残余视频相位补偿处理后的回波中频信号进行一维距离向压缩处理；在得到了Q个一维距离向压缩处理数据之后，将Q个一维距离向压缩处理数据组成Q行N列矩阵，其中N为距离向采集点数，Q为总的采集周期数；对所得到的矩阵信号组进行列向FFT变换；在变换得到的频域中寻找尖峰位置；计算所述尖峰位置处的滑坡速度值作为所估计的泥石流滑坡速度。

可选地，所估计的泥石流滑坡速度为：

其中，f₀为中心频率，Br为信号带宽，V为所估计的泥石流滑坡速度，C为电磁波速度，f_v为目标的多普勒频率值。

通过采用上述技术方案，能够全天候全天时地监测泥石流区域，提供泥石流区域的实时成像与变化监测，能够对泥石流隐患区域实现大面覆盖监测，而且还能够最大限度地降低系统复杂度，降低中频信号带宽，降低对数字采样的要求，提高了泥石流监测的精度。

本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。

附图说明

附图是用来提供对本公开的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本公开，但并不构成对本公开的限制。在附图中：

图1示出根据本公开一种实施例的毫米波雷达泥石流智能监测预警方法的流程图。

图2示出根据本公开一种实施例的毫米波雷达泥石流智能监测预警系统的示意框图。

图3示出了根据本公开实施例的毫米波雷达泥石流智能监测预警系统的成像几何。

图4示出根据本公开一种实施例的毫米波雷达泥石流智能监测预警方法的流程图。

具体实施方式

以下结合附图对本公开的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本公开，并不用于限制本公开。

图1示出根据本公开一种实施例的毫米波雷达泥石流智能监测预警方法的流程图。如图1所示，该方法包括：

在步骤S11中，向被监测的泥石流区域发射第一微波信号；

在步骤S12中，接收被泥石流区域反射回来的回波信号；

在步骤S13中，对反射回来的回波信号进行下变频处理，得到回波中频信号；

在步骤S14中，对回波中频信号进行采样；

在步骤S15中，基于采样后的回波中频信号对被监测的泥石流区域进行成像，得到关于被监测的泥石流区域的图像；

在步骤S16中，利用图像监测被监测的泥石流区域的变化。

通过采用上述技术方案，能够全天候全天时地监测泥石流区域，提供泥石流区域的实时成像与变化监测，能够对泥石流隐患区域实现大面覆盖监测，而且还能够最大限度地降低系统复杂度，降低中频信号带宽，降低对数字采样的要求，提高了泥石流监测的精度。

图2示出根据本公开一种实施例的毫米波雷达泥石流智能监测预警系统的示意框图。如图2所示，该系统包括：阵列天线11，用于向被监测的泥石流区域发射第一微波信号，并接收被泥石流区域反射回来的回波信号；MIMO模块13，用于对反射回来的回波信号进行下变频处理，得到回波中频信号；采样模块14，用于对回波中频信号进行采样；成像模块15，用于基于采样后的回波中频信号对被监测的泥石流区域进行成像，得到关于被监测的泥石流区域的图像；变化监测模块16，用于利用图像监测被监测的泥石流区域的变化。

第一微波信号可以由毫米波雷达泥石流智能监测预警系统外部的设备产生，也可以由MIMO模块13产生。

进一步参考图2，根据本公开实施例的毫米波雷达泥石流智能监测预警系统还可以包括微波开关网络12，用于通过不同开关网络切换逐一传输MIMO模块13产生的第一微波信号至阵列天线11，以及通过不同开关网络切换将阵列天线11接收到的反射回来的回波信号传输给MIMO模块13。毫米波雷达泥石流智能监测预警系统还可以包括系统控制器18，用于对微波开关网络12、MIMO模块13、采样模块14、成像模块15、变化监测模块16进行控制。毫米波雷达泥石流智能监测预警系统还可以包括预警模块17，用于对泥石流进行预警。

通过采用上述技术方案，能够全天候全天时地监测泥石流区域，提供泥石流区域的实时成像与变化监测，能够对泥石流隐患区域实现大面覆盖监测，而且还能够最大限度地降低系统复杂度，降低中频信号带宽，降低对数字采样的要求，提高了泥石流监测的精度。

下面结合图2所示的毫米波雷达泥石流智能监测预警系统结构图，详细描述毫米波雷达泥石流智能监测预警系统的各个组成模块以及如何进行泥石流监测与预警。

图2中的阵列天线11可以由多幅在水平、垂直方向上排列的独立天线阵元构成。阵列天线11可以包括发射阵列天线和接收阵列天线。

微波开关网络12由发射阵列开关网络、接收阵列开关网络、驱动器和微波开关集中控制器构成。发射阵列开关网络通过不同开关网络切换逐一传输MIMO模块13产生的第一微波信号至阵列天线11中的发射阵列天线，而后将阵列天线11中的接收阵列天线接收到的由被监测泥石流区域反射回来的回波信号、通过接收阵列开关网络传输给MIMO模块13进行信号接收。驱动器用于对发射阵列开关网络、接收阵列开关网络进行驱动，微波开关集中控制器用于对驱动器进行控制。

MIMO模块13可以在系统控制器18的控制下，通过MIMO收发控制器控制频率源产生微波信号，产生的微波信号经过功率放大器放大后被传输至耦合器，然后，耦合器会输出两路微波信号s₁(t)和s₂(t)。第一路微波信号s₁(t)经过功率放大器放大后被传输至微波开关网络12中的发射阵列开关网络，并通过阵列天线11中的发射阵列天线辐射出去。第二路微波信号s₂(t)用于在后续的下变频处理中使用。另外，如果阵列天线11的接收阵列天线包括多路接收阵列天线，则第二路微波信号s₂(t)经过放大器和功分器后会输出与多路接收阵列天线相对应的多路信号。例如，假设接收阵列天线包括2路接收阵列天线，这2路接收阵列天线分别接收回波信号s_re(t-τ₁)和s_re(t-τ₂)，则第二路微波信号s₂(t)经过放大器和功分器后会输出2路信号s₂₁(t)和s₂₂(t)，然后，回波信号s_re(t-τ₁)和s_re(t-τ₂)经阵列天线11的接收阵列天线和微波开关网络12的接收阵列开关网络后，其中一路回波信号s_re(t-τ₁)经过放大器后通过混频器与s₂₁(t)进行下变频处理，而后通过滤波器和放大器形成第一路中频信号“IF1”S_IF1(t)，同时，另一路回波信号s_re(t-τ₂)经过放大器后通过混频器与s₂₂(t)进行下变频处理，而后通过滤波器和放大器形成第二路中频信号“IF2”S_IF2(t)。其中，t为距离向时间变量，且t∈[-T_r/2,T_r/2]；τ为被监测目标到阵列天线11的电磁波传播时间。为了保证系统总体工作效率和实现高分辨率成像，频率源产生的微波信号为调频连续波(FMCW)信号，信号带宽为B_r，持续时间为T_r。通过MIMO模块13的如此处理，能够最大限度地降低系统复杂度，降低中频信号带宽，从而降低数字采样的要求。

在一种实施方式中，MIMO模块13发射的信号S_tr(t)为：

其中，f_c为系统工作频率，t为距离向时间变量，且t∈[-T_r/2,T_r/2]，T_r为信号持续时间，K_r为信号调频率，信号带宽为B_r＝K_rT_r，j为虚数符号。

MIMO模块13接收的反射回来的回波信号经下混频器混频后得到的中频回波信号的表达式为S_IF(t)：

其中，y_n为阵列天线11中第n个阵元的方位向坐标，阵列天线11的最大阵元数N为

d为阵元间距，R(y_n)为阵列天线11阵元y_n到被监测目标P_n的距离，

为阵列天线11阵元y_n的目标散射系数，Σ为场景中所有被监测目标的求和运算，j为虚数符号，pi为圆周率，C为光速，L为虚拟阵列长度，k为表示场景中第k个目标。

然后MIMO模块13会对回波中频信号进行残余视频相位补偿，得到：

其中，

为目标散射系数，n为第n个阵列虚拟单元。

采样模块14通过模数转换器(Analog to Digital converter,ADC对前面公式(3)中得到的中频信号进行量化，量化位数为例如8～14bit，采样率f_s为信号带宽B_r的2～3倍，通常取2.4倍。系统控制器18对采样模块14的采样操作进行控制。经过采样模块14的采样操作之后，会得到相应的数字信号，例如如果经过MIMO模块13的下变频处理之后得到2路中频信号，则经过采样模块14的采样操作之后也会得到对应的2路数字信号。然后所得到的数字信号被传输到成像模块15中进行成像处理。

成像模块15的成像处理操作主要包括以下步骤。

步骤S151：对公式(3)中的回波中频信号进行距离向逆傅里叶变换，变换后的表达式为：

其中，C为光速。

步骤152：对距离向信号进行分割。也即，对步骤S151中的S_IFFT(t,y_n)，根据判据公式

进行分割，从而将距离向逆傅里叶变换后的回波中频信号分割成近距离向回波中频信号和远距离向回波中频信号，也即：

其中，

K为距离向的采样点数。

步骤S153：基于近距离向回波中频信号执行近距离方位向叠加压缩成像，其中，图像的方位向表示为以阵列天线11的中心为原点、以垂直于阵列天线11的方向作为法向对用户需要的角度范围θ_rang进行等间隔角度划分，共划分为M份；以变量θ_m表示，其坐标刻度如下式所示：

图像的距离向变量ρ表示为ρ∈[0～R_min]，其中

刻度为

其中B为信号带宽(其与B_r含义一样，可以换成B_r)；m为距离单元索引。则，近距离方位向叠加压缩成像表达式为：

其中，H为图像距离向采样点数，h为采样点的序号，ρ_hm为图像上m行h列到阵列天线11中的阵元的距离。

步骤S154：基于远距离向回波中频信号执行远距离方位向叠加压缩成像，其中图像方位向的表示同步骤153中的近距离方位向叠加压缩成像一致。图像距离向的刻度与步骤153中的一致但取值范围为ρ∈[R_min，R_max]，方位向通过构建加权函数来实现方位向压缩成像，则加权函数表达式为：

则方位向压缩成像的表达式为：

设场景中目标P(R_P,θ_P)，根据幂级数叠加公式对上式进行合成以消除求和符号即可得：

其中，C为光速，r为实际的距离变量，R_P表示P目标距虚拟阵列天线的距离、θ_P为P目标与虚拟阵列天线法线的夹角。

步骤155：对近距离图像与远距离图像进行组合拼接，即将I_近(ρ,θ_m)图像与I_远(ρ,θ_m)图像进行拼接，即可得到一副完整关于泥石流区域的图像(例如二维图像、三维图像等)I(ρ，θ)，即为：

图3示出了根据本公开实施例的毫米波雷达泥石流智能监测预警系统的成像几何。成像模块15通过采用近距离和远距离分别成像，弥补了阵列天线11成像算法上对近距离成像散焦的缺点。

变化监测模块16利用图像监测被监测的泥石流区域的变化，包括：选取不同时相的两幅图像；基于不同时相的两幅图像，利用交叉熵进行泥石流区域的变化检测处理。

在一种实施方式中，变化监测模块16所执行的操作可以包括：

步骤S161：执行滤波处理，例如，执行Lee滤波处理流程以降低图像相干斑噪声。

步骤S162：执行变化检测提取，也即选取两幅不同时相的图像，利用交叉熵来进行变化监测处理。步骤S162可以包括以下若干步骤：

步骤S1621：构造逼近真实统计分布模型的多项式，选取两幅不同时相的图像，其窗口大小为M*N，分别选取两时相图像中某同一个像素点为中心原点，并以M*N固定窗口大小选取相应的样本数据，构造如下式所示的多项式来近似样本图像统计分布模型的概率密度函数，则表达式为：

f(x)＝G_amma(x)+E(x)G_amma(x) (11)

其中，x为概率密度的随机变量；G_amma(x)为广义GAMMA分布，E(x)如下所示：

其中，n为样本数量，H₃(x)、H₄(x)、H₆(x)为Hermite多项式，分别如下所示：

变量x为数据样本的随机变量，其如下式所示：

其中，a_i表示变量x的i阶原点矩，即

其中x₁,x₂,…x_n为样本数据。

步骤S1622：计算交叉熵，也即，通过上步骤中获取的两幅图像窗口样本的概率密度函数(分别记为

和

)带入到交叉熵的计算公式中，即为：

其中，

步骤S1623：确定变化阈值T，其具体的过程如下：

步骤S16231：建立统计直方图，根据差异图像构建统计直方图，选择直方图中峰值处的变化值T0；

步骤S16232：确定阈值T，以T0开始一直到Tmax，(注意：T0<T1<T2....<Tmax)，一直计算当前位置处Tn的峰值与Tn-1处的峰值比，选择第一个该峰值比小于1时的Tn为最佳阈值T；

步骤S1624：确定像素变化状态，根据上一步骤确定交叉熵的阈值T，用步骤S1622中求得交叉熵KLD与设定的阈值T进行比较，若大于，则将图像中某位置像素标记为变化，若小于，则将图像中某位置像素标记为未变化；

步骤S1625：循环迭代，移动下一位置像素，重新执行步骤S1621～S1624，直到将图像中的像素全部遍历完毕后跳出循环，即显示图像变化部分。

图4示出根据本公开一种实施例的毫米波雷达泥石流智能监测预警方法的流程图，其在图1的基础上，还包括步骤S17，基于监测到的变化进行泥石流预警。其中，可以采用面积预警的方式进行预警。也即选定初始图像，然后利用变化监测模块16记录并统计后续时刻获取的图像与初始图像之间的变化像素，设面积预警的预设阈值为TT，即当变化的像素数量超过TT时，则进行泥石流滑坡预警，当小于TT时，则继续监测。

另外，在发生面积预警的情况下，还需要进一步估计泥石流滑坡速度，以便增强预警效果。

在估计泥石流滑坡速度的情况下，首先，将毫米波雷达泥石流智能监测预警系统的工作模式切换为单通道模式，也即阵列开关网络中只开启某固定发射通道和接收通道，并且实现每秒f_M次的收发频率，并将信号进行存储。该收发频率直接影响速度分辨力，可根据场景的需求灵活设置。下面将具体阐述滑坡速度估计的处理过程。

首先进行一维距离向逆傅里叶变换，其中，发射的微波信号与前文公式(1)相同，该微波信号经过混频器后输出的中频信号的表达式与公式(2)类似，即为：

其中，

为目标散射系数；Σ为场景中所有被监测目标的求和运算；R_k为第k个被监测目标到阵列天线11相位中心的距离；K_r为调频率；C为光速；t为距离向时间。

然后对公式(16)执行残余视频相位补偿处理，将接收到的信号分别进行一维距离向压缩处理；其具体过程与上文所述一致，也即压缩后的表达式为：

然后，获取多普勒速度。也即，设定总的采集周期数为Q次，则每达到Q次后就对采集到的数据进行处理，将Q次采集到的信号按照如下所示的方式排列，即组成一个Q行N列矩阵，其中Q为采集次数，N为距离向采集点数；

其中，

为第n次采集的信号，n取值为[1,2,…,Q]。

然后对上述信号组进行列向FFT变换，变换到频域，通过寻找尖峰位置即可求出相应距离位置处的速度值，其速度与频率的换算关系表达式为：

其中，f₀为中心频率，Br为信号带宽，V为所估计的泥石流滑坡速度，C为电磁波速度，f_v为目标的多普勒频率值。

另外，根据本公开实施例的毫米波雷达泥石流智能监测预警系统可选择工作频率范围为77GHz～79GHz，在此工作频率下，可以充分发挥系统的全天候、全天时以及不受气候和环境等因素的影响的优势，能够实现对泥石流隐患区域的点云成像。

以上结合附图详细描述了本公开的优选实施方式，但是，本公开并不限于上述实施方式中的具体细节，在本公开的技术构思范围内，可以对本公开的技术方案进行多种简单变型，这些简单变型均属于本公开的保护范围。

另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本公开对各种可能的组合方式不再另行说明。

此外，本公开的各种不同的实施方式之间也可以进行任意组合，只要其不违背本公开的思想，其同样应当视为本公开所公开的内容。

Claims (6)

1.一种毫米波雷达泥石流智能监测预警方法，其特征在于，包括：

产生两路相同的微波信号，其中一路微波信号为第一微波信号，另一路微波信号为第二微波信号；

向被监测的泥石流区域发射所述第一微波信号；

接收被所述泥石流区域反射回来的回波信号；

将反射回来的回波信号通过混频器与所述第二微波信号进行下变频处理，并对下变频处理后的反射回来的回波信号进行滤波和放大，得到回波中频信号；

对所述回波中频信号进行采样；

基于采样后的回波中频信号对被监测的泥石流区域进行成像，得到关于被监测的泥石流区域的图像；

利用所述图像来监测被监测的泥石流区域的变化；

基于监测到的变化进行泥石流预警；

其中，所述方法还包括：估计泥石流滑坡速度；

其中，所述估计泥石流滑坡速度包括：

在下变频处理得到的回波中频信号经历残余视频相位补偿处理之后，对残余视频相位补偿处理后的回波中频信号进行一维距离向压缩处理；

在得到了Q个一维距离向压缩处理数据之后，将Q个一维距离向压缩处理数据组成Q行N列矩阵，其中N为距离向采集点数，Q为总的采集周期数；

对所得到的矩阵信号组进行列向FFT变换；

在变换得到的频域中寻找尖峰位置；

计算所述尖峰位置处的滑坡速度值作为所估计的泥石流滑坡速度。

2.根据权利要求1所述的方法，其特征在于，所述基于监测到的变化进行泥石流预警，包括：

在监测到的变化指示所述图像中变化的像素数量超过预设阈值时，进行泥石流滑坡预警。

3.根据权利要求1所述的方法，其特征在于，所估计的泥石流滑坡速度为：

其中，f₀为中心频率，Br为信号带宽，V为所估计的泥石流滑坡速度，C为电磁波速度，f_v为目标的多普勒频率值。

4.一种毫米波雷达泥石流智能监测预警系统，其特征在于，包括：

阵列天线，用于向被监测的泥石流区域发射第一微波信号，并接收被所述泥石流区域反射回来的回波信号；

MIMO模块，用于将反射回来的回波信号通过混频器与第二微波信号进行下变频处理，并对下变频处理后的反射回来的回波信号进行滤波和放大，得到回波中频信号，其中所述第二微波信号为与所述第一微波信号相同的另一微波信号；

采样模块，用于对所述回波中频信号进行采样；

成像模块，用于基于采样后的回波中频信号对被监测的泥石流区域进行成像，得到关于被监测的泥石流区域的图像；

变化监测模块，用于利用所述图像来监测被监测的泥石流区域的变化；

预警模块，用于基于监测到的变化进行泥石流预警；

其中，所述预警模块还用于：估计泥石流滑坡速度；

其中，所述预警模块还用于：

在下变频处理得到的回波中频信号经历残余视频相位补偿处理之后，对残余视频相位补偿处理后的回波中频信号进行一维距离向压缩处理；

在得到了Q个一维距离向压缩处理数据之后，将Q个一维距离向压缩处理数据组成Q行N列矩阵，其中N为距离向采集点数，Q为总的采集周期数；

对所得到的矩阵信号组进行列向FFT变换；

在变换得到的频域中寻找尖峰位置；

计算所述尖峰位置处的滑坡速度值作为所估计的泥石流滑坡速度。

5.根据权利要求4所述的系统，其特征在于，所述预警模块用于：

在监测到的变化指示所述图像中变化的像素数量超过预设阈值时，进行泥石流滑坡预警。

6.根据权利要求4所述的系统，其特征在于，所估计的泥石流滑坡速度为：

其中，f₀为中心频率，Br为信号带宽，V为所估计的泥石流滑坡速度，C为电磁波速度，f_v为目标的多普勒频率值。

Images