CN114184256A - 一种多目标背景下的水位测量方法 - Google Patents

一种多目标背景下的水位测量方法 Download PDF

Info

Publication number
CN114184256A
CN114184256A CN202111493031.1A CN202111493031A CN114184256A CN 114184256 A CN114184256 A CN 114184256A CN 202111493031 A CN202111493031 A CN 202111493031A CN 114184256 A CN114184256 A CN 114184256A
Authority
CN
China
Prior art keywords
target
water level
targets
radar
spectral line
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111493031.1A
Other languages
English (en)
Inventor
朱世平
陈德莉
袁开见
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuxi Hangzheng Technology Co ltd
Original Assignee
Wuxi Hangzheng Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuxi Hangzheng Technology Co ltd filed Critical Wuxi Hangzheng Technology Co ltd
Priority to CN202111493031.1A priority Critical patent/CN114184256A/zh
Publication of CN114184256A publication Critical patent/CN114184256A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • G01F23/22—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
    • G01F23/28—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
    • G01F23/284—Electromagnetic waves

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Thermal Sciences (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种多目标背景下的水位测量方法，包括步骤：雷达水位计发射天线以固定周期发射一组调频连续波信号；三个接收天线接收采样后的第一锯齿波回波信号和第二锯齿波回波信号；将信号处理得到离散傅里叶频谱；在频谱上用CFAR算法检测出局部幅值最大的m个目标的谱线位置，求解m个目标的四元组属性值；通过若干周期的测量，得到若干周期的检测得到目标的四元组属性值；将目标分类，并根据预设规则选择类内的目标类为待定目标类，将距离均值最小的目标类判定为水面目标。本发明可准确识别出水体目标，排除干扰目标的影响，提高水位测量的准确性，降低对雷达水位计的安装要求。

Description

一种多目标背景下的水位测量方法
技术领域
本发明属于水位测量技术领域，尤其涉及一种多目标背景下的水位测量方法。
背景技术
目前市面上的调频连续波体制雷达水位计和脉冲体制雷达水位计都是单个发射天线，单个接收天线或者收发共用一个天线的形式，只能获取目标的距离和回波强度信息。当雷达水位计波束照射范围内存在如建筑物边壁、凸台等人造物，树木、水草等自然物甚至行船等移动目标时，雷达回波中存在多目标干扰。现有的解决措施有：在自然河道和渠道场景下，通过加长雷达水位计的横臂安装支架，将雷达水位计安装位置尽量远离河岸，伸入水体区域，减小雷达波照射边壁、树木和水草的几率，这种安装方式无疑会增加水位监测设施的土建成本，且依然无法排除行船等动目标的干扰，容易出现水位跳变；在城市地下排水管网窨井水位监测场景下，通过外壳结构设计(如透镜天线结构)以减小雷达水位计的波束照射范围，或者将水位计安装在窨井正中心位置，以减小井边壁及金属踏步的反射，这种安装方式对现场安装的要求高，且会阻碍后期工人下井进行管道清淤等工作。因此，现有的雷达水位计测量方法难以实现多目标背景下的水体目标的正确识别，以及高精度、稳定准确的水位测量要求。
另外，现有的提高雷达水位测量精度的方法主要有：在回波的频率域，利用频谱差值法对目标频率进行细化处理，以提高目标空高和水位测量的精度。在多目标存在的情况下，该方法的计算效率和实时处理能力有待提升。
发明内容
为简化水位测量的计算方法，提高计算效率，增强多目标背景下水体目标的检测识别能力，提高水位测量准确度和稳定性，本发明提出了一种多目标背景下的水位测量方法。
本发明公开的一种多目标背景下的水位测量方法，使用雷达水位计测量水位，所述雷达水位计包括一个单阵元发射天线TX和三个单阵元接收天线RX1、RX2和RX3，所述水位测量方法包括以下步骤：
雷达水位计通过发射天线TX以固定周期Tframe发射一组调频连续波信号，每组调频连续波由第一锯齿波S1和第二锯齿波S2组成，两段锯齿波间隔Tshape时间发射；
三个接收天线接收采样后的第一锯齿波回波信号S1 RX1(n)，S1 RX2(n)，S1 RX3(n)和第二锯齿波回波信号S2 RX1(n)，S2 RX2(n)，S2 RX3(n)，n∈[1,N]；将S1 RX1(n)，S1 RX2(n)，S1 RX3(n)和S2 RX3(n)加N点汉明窗后，做N点离散傅里叶变换后，得到离散傅里叶频谱
Figure BDA0003399116580000021
和
Figure BDA0003399116580000022
其中Mi(k)和
Figure BDA0003399116580000023
为第i个频谱在谱线位置k的频谱幅值和相角，k∈[1,N]；
在S1 RX3(n)的频谱上用CFAR算法在k∈[1,N/2]范围内检测出局部幅值最大的m个目标的谱线位置k1、k2、...、km，求解所述m个目标的四元组属性值；
通过若干周期Tframe的测量，得到若干周期Tframe的所述m个目标的四元组属性值；
将所述若干目标分类，并根据预设规则选择类内的目标类为待定目标类，将距离均值最小的目标类判定为水面目标。
进一步的，所述三个接收天线两两构成正交关系，且天线阵元中心间距Di相等，间距小于等于雷达发射信号波长的一半，定义三个接收天线的中心点为原点O，RX2和RX3中心点连线方向为Y轴，RX1和RX3中心点连线方向为X轴，建立参考直角坐标系。
进一步的，所述目标的四元组属性为(r,θ,ρ,v)，若目标W在第一平面上的投影点为W'，在第二平面上的投影点为W"，OZ为雷达中心点到水平面的垂线，则θ为OW'与OZ的夹角，ρ为OW"与OZ的夹角，r为雷达中心点O到目标W的距离，v为目标W在OW径向上的相对速度，所述第一平面和第二平面互相垂直，并分别和水平面垂直。
进一步的，单段锯齿波形状的调频连续波发射时起始频率为f0，结束频率为f1，带宽为BW，发射时间为Tp，且Tp<Tshape，其余时间段雷达不工作。
进一步的，所述求解m个目标的四元组属性值包括以下步骤：
对于第i个目标，1<＝i<＝m，
步骤一：在S1 RX3(n)的频谱谱线ki的相邻两条谱线ki-1和ki+1中找到幅值最大的谱线，假设谱线幅值M3(ki)<M3(ki+1)，则以谱线ki和ki+1的频谱幅值M3(ki)和M3(ki+1)进行步骤二计算；
步骤二：计算目标i的距离ri为：
Figure BDA0003399116580000031
其中c为电磁波在真空中的传播速度，a为系数，Bw为带宽；
步骤三：计算目标i的角度参量(θi,ρi)：
Figure BDA0003399116580000032
Figure BDA0003399116580000033
其中λ0＝c/f0，为FMCW起始扫频频率f0对应的波长；
Figure BDA0003399116580000041
为接收天线RX1和RX3的回波信号频谱在谱线ki处的相位差，计算方法如下：
Figure BDA0003399116580000042
Figure BDA0003399116580000043
为接收天线RX2和RX3的回波信号频谱在谱线ki处的相位差，计算方法如下：
Figure BDA0003399116580000044
步骤四：计算目标i的径向相对速度vi：
根据所述系数a，先求取目标i在频谱
Figure BDA0003399116580000045
和
Figure BDA0003399116580000046
中对应频率k＝(ki-1+a)处的相角
Figure BDA0003399116580000047
和
Figure BDA0003399116580000048
Figure BDA0003399116580000049
Figure BDA00033991165800000410
再求目标i的径向相对速度vi：
Figure BDA00033991165800000411
其中
Figure BDA00033991165800000412
vi为正数时，表示目标i的运动方向靠近雷达水位计；vi为负数时，表示目标i的运动方向远离雷达水位计；Tshape是决定目标的径向速度的测量范围的参数；
重复上述步骤一至步骤四过程，得到当前frame雷达数据检测到的其余目标的四元组属性值。
进一步的，所述系数a的计算公式为：
Figure BDA0003399116580000051
其中p0,p1,p2,p3为预设的系数，M3(ki)和M3(ki+1)为谱线ki和ki+1的幅值；
Tshape通过下式计算：
Figure BDA0003399116580000052
其中vmax为检测目标的最大径向速度。
进一步的，按照以下方法完成水面目标识别：
步骤一：所有目标以距离这一属性值为条件，通过K均值聚类算法，分为若干组不同距离的目标类；
步骤二：对类内目标，计算距离均值r，空间立体角θ和ρ的均方差θstd和ρstd，径向速度的均值v和均方差vstd；
步骤三：以空间立体角均方差为判定条件，将同时满足θstd≥Th1和ρstd≥Th1的目标类作为待定目标类，其中Th1为设定的门限值；
步骤四：计算待定目标类径向速度均方差与均值之比，在满足
Figure BDA0003399116580000053
条件下，将距离均值
Figure BDA0003399116580000054
最小的目标类判定为水面目标，其中Th2为设定的门限值。
进一步的，识别出水面目标后，将识别为水面目标类的距离均值r换算为雷达水位计距离水面的空高值h:
Figure BDA0003399116580000055
其中
Figure BDA0003399116580000056
为水面目标类的θ均值，
Figure BDA0003399116580000057
为水面目标类的ρ均值。
进一步的，用设定的安装高度减去空高值得到当前水深值。
本发明的有益效果如下：
在雷达水位计已有的测量目标距离和回波强度的功能基础上，利用一发多收调频连续波体制的优势，快速解算出目标相对雷达水位计的距离、空间立体角和径向运动速度，在雷达水位计安装现场存在不可排除的多目标干扰的情况下，利用目标的距离、回波强度、空间立体角和径向运动速度等属性，准确识别出水体目标，排除干扰目标的影响，提高水位测量的准确性，降低对雷达水位计的安装要求。
附图说明
图1本发明的水位测量方法流程图；
图2本发明在典型应用场景下的雷达水位计安装示意图；
图3本发明的天线示意图；
图4目标三维示意图；
图5本发明的调频连续波发射信号形式。
具体实施方式
下面结合附图对本发明作进一步的说明，但不以任何方式对本发明加以限制，基于本发明教导所作的任何变换或替换，均属于本发明的保护范围。
本发明使用雷达水位计进行水位测量，雷达水位计包括一个单阵元发射天线TX和三个单阵元接收天线RX1、RX2和RX3，如图1所示。本实施例中以4个目标的识别任务为例进行识别，但对于2个或3个目标或超过4个目标的场景，本发明也同样适用，本实施例对此不作限制。三个接收天线两两构成正交关系，且天线阵元中心间距Di相等，间距小于等于雷达发射信号波长的一半，定义三个接收天线的中心点为原点，RX2和RX3中心点连线方向为Y轴，RX1和RX3中心点连线方向为X轴，建立参考直角坐标系。
水位测量方法包括以下步骤：
雷达水位计通过发射天线TX以固定周期Tframe发射一组调频连续波信号，每组调频连续波由第一锯齿波S1和第二锯齿波S2组成，两段锯齿波间隔Tshape时间发射；单段锯齿波形状的调频连续波发射时起始频率为f0，结束频率为f1，带宽为BW，发射时间为Tp，且Tp<Tshape，其余时间段雷达不工作。
三个接收天线接收采样后的第一锯齿波回波信号S1 RX1(n)，S1 RX2(n)，S1 RX3(n)和第二锯齿波回波信号S2 RX1(n)，S2 RX2(n)，S2 RX3(n)，n∈[1,N]；将S1 RX1(n)，S1 RX2(n)，S1 RX3(n)和S2 RX3(n)加N点汉明窗后，做N点离散傅里叶变换后，得到离散傅里叶频谱
Figure BDA0003399116580000071
和
Figure BDA0003399116580000072
其中Mi(k)和
Figure BDA0003399116580000073
为第i个频谱在谱线位置k的频谱幅值和相角，k∈[1,N]；
在S1 RX3(n)的频谱上用CFAR算法在k∈[1,N/2]范围内检测出局部幅值最大的m个目标的谱线位置k1、k2、...、km，求解m个目标的四元组属性值；
目标的四元组属性为(r,θ,ρ,v)，若目标W在第一平面上的投影点为W'，在第二平面上的投影点为W"，OZ为雷达中心点到水平面的垂线，则θ为OW'与OZ的夹角，ρ为OW"与OZ的夹角，r为雷达中心点O到目标W的距离，v为目标W在OW径向上的相对速度，所述第一平面和第二平面互相垂直，并分别和水平面垂直。
通过若干周期Tframe的测量，得到若干周期Tframe的m个目标的四元组属性值；
求解m个目标的四元组属性值包括以下步骤：
对于第i个目标，1<＝i<＝m，
步骤一：在S1 RX3(n)的频谱谱线ki的相邻两条谱线ki-1和ki+1中找到幅值最大的谱线，假设谱线幅值M3(ki)<M3(ki+1)，则以谱线ki和ki+1的频谱幅值M3(ki)和M3(ki+1)进行步骤二计算；
步骤二：计算目标i的距离ri为：
Figure BDA0003399116580000081
其中c为电磁波在真空中的传播速度，a为系数，Bw为带宽；
系数a的计算公式为：
Figure BDA0003399116580000082
其中p0,p1,p2,p3为预设的系数，M3(ki)和M3(ki+1)为谱线ki和ki+1的幅值；
Tshape通过下式计算：
Figure BDA0003399116580000083
其中vmax为检测目标的最大径向速度。
步骤三：计算目标i的角度参量(θi,ρi)：
Figure BDA0003399116580000084
Figure BDA0003399116580000085
其中λ0＝c/f0，为FMCW起始扫频频率f0对应的波长；
Figure BDA0003399116580000086
为接收天线RX1和RX3的回波信号频谱在谱线ki处的相位差，计算方法如下：
Figure BDA0003399116580000087
Figure BDA0003399116580000088
和
Figure BDA0003399116580000089
为第1个和第3个频谱在谱线位置ki的相角；
Figure BDA00033991165800000810
为接收天线RX2和RX3的回波信号频谱在谱线ki处的相位差，计算方法如下：
Figure BDA0003399116580000091
Figure BDA0003399116580000092
和
Figure BDA0003399116580000093
为第2个和第3个频谱在谱线位置ki的相角。
步骤四：计算目标i的径向相对速度vi：
根据系数a，先求取目标i在频谱
Figure BDA0003399116580000094
和
Figure BDA0003399116580000095
中对应频率k＝(ki-1+a)处的相角
Figure BDA0003399116580000096
和
Figure BDA0003399116580000097
Figure BDA0003399116580000098
Figure BDA0003399116580000099
Figure BDA00033991165800000910
和
Figure BDA00033991165800000911
为第3个和第4个频谱在谱线位置ki的相角，
Figure BDA00033991165800000912
和
Figure BDA00033991165800000913
为第3个和第4个频谱在谱线位置ki+1的相角；
再求目标i的径向相对速度vi：
Figure BDA00033991165800000914
其中
Figure BDA00033991165800000915
vi为正数时，表示目标i的运动方向靠近雷达水位计；vi为负数时，表示目标i的运动方向远离雷达水位计；Tshape是决定目标的径向速度的测量范围的参数；
重复上述步骤一至步骤四过程，得到当前frame雷达数据检测到的其余目标的四元组属性值。
将若干目标分类，并根据预设规则选择类内的目标类为待定目标类，将距离均值最小的目标类判定为水面目标。
按照以下方法完成水面目标识别：
步骤一：所有目标以距离这一属性值为条件，通过K均值聚类算法，分为若干组不同距离的目标类；
步骤二：对类内目标，计算距离均值
Figure BDA0003399116580000101
空间立体角θ和ρ的均方差θstd和ρstd，径向速度的均值
Figure BDA0003399116580000102
和均方差vstd；
步骤三：以空间立体角均方差为判定条件，将同时满足θstd≥Th1和ρstd≥Th1的目标类作为待定目标类，其中Th1为设定的门限值；
步骤四：计算待定目标类径向速度均方差与均值之比，在满足
Figure BDA0003399116580000103
条件下，将距离均值
Figure BDA0003399116580000104
最小的目标类判定为水面目标，其中Th2为设定的门限值。
识别出水面目标后，将识别为水面目标类的距离均值
Figure BDA0003399116580000105
换算为雷达水位计距离水面的空高值h:
Figure BDA0003399116580000106
其中
Figure BDA0003399116580000107
为水面目标类的θ均值，
Figure BDA0003399116580000108
为水面目标类的ρ均值。
进一步的，用设定的安装高度减去空高值得到当前水深值。
实施例
雷达水位计安装在城市排水污水管网窨井口，固定在井壁，雷达波束向下照射，如图2所示，在雷达的天线波束角照射范围内，有水面目标W，还有井壁上金属踏步A1和A2，井底的水泥留台B等干扰目标。本实施例以4个目标为例，对于多于1个的目标检测也同样适用，本实施例对此不作限制。
雷达水位计采用一发三收的天线形式(一个单阵元发射天线TX，三个单阵元接收天线RX1、RX2和RX3)，三个接收天线两两构成正交关系，且天线阵元中心间距Di相等，间距小于等于雷达发射信号波长的一半，如图2所示。定义三个接收天线的中心点为原点，RX2和RX3中心点连线方向为Y轴，RX1和RX3中心点连线方向为X轴，建立参考直角坐标系。
以图3的直角坐标系为基准，波束向下方向为Z轴方面建立三维坐标系，图2示意的场景中四个目标(W、A1、A2和B)的三维坐标关系如图4所示。以水面目标W为例，定义该目标在XOZ平面上的投影点为W'，OW'与OZ的夹角为θW；该目标在YOZ平面上的投影点为W"，OW"与OZ的夹角为ρW；定义雷达中心点O到目标W的距离为rW；另外，定义目标W在OW径向上的相对速度为vW。因此，将(rW,θW,ρW,vW)四个参数定义为水面目标W的属性，同理，目标A1的属性为(rA1,θA1,ρA1,vA1)，目标A2的属性为(rA2,θA2,ρA2,vA2)，目标B的属性为(rB,θB,ρB,vB)。
雷达水位计通过发射天线TX以固定周期Tframe发射一组调频连续波信号，每组调频连续波由两段相同锯齿波S1和S2组成，两段锯齿波间隔Tshape时间发射。单段锯齿波形状的调频连续波发射时起始频率为f0，结束频率为f1，带宽为BW，发射时间为Tp(需满足Tp<Tshape)，其余时间段雷达不工作，发射信号的形状如图5所示。
雷达水位计发射S1信号时，三个接收天线同时接收到四个目标的反射信号，经过混频下变频、中频放大和AD采样后，得到三组N点数字回波信号S1 RX1(n)，S1 RX2(n)，S1 RX3(n)，n∈[1,N]，其中N为采样点数。同理，雷达水位计发射S2信号时，三个接收天线得到的采样后的数字回波信号表示为S2 RX1(n)，S2 RX2(n)，S2 RX3(n)。
对S1 RX1(n)，S1 RX2(n)，S1 RX3(n)和S2 RX3(n)分别加N点汉明窗(hamming window)后，做N点离散傅里叶变换，得到对应的离散傅里叶频谱可表示为
Figure BDA0003399116580000111
和
Figure BDA0003399116580000112
其中Mi(k)和
Figure BDA0003399116580000113
为第i个频谱在谱线位置k的频谱幅值和相角，k∈[1,N]。
在S1 RX3(n)的频谱
Figure BDA0003399116580000121
上用CFAR算法在k∈[1,N/2]范围内检测出局部幅值最大的四个目标(A1、A2、B和W)的谱线位置k1、k2、k3和k4，下面以W目标为例，求解该目标的属性值(rW,θW,ρW,vW)。
步骤一：在谱线k4的相邻两条谱线k4-1和k4+1中找到幅值最大的谱线，假如M3(k4-1)<M3(k4+1)，则以谱线k4和k4+1的频谱值进行下一步计算。
步骤二：计算目标W的距离rW为：
Figure BDA0003399116580000122
其中c为电磁波在真空中的传播速度，系数a的计算公式为：
Figure BDA0003399116580000123
通过系数a的计算，可将现有的Chirp-Z变换、频谱差值法简化为一个与幅值之比相关的一元三次函数关系，该函数的四个待定系数p0,p1,p2,p3通过仿真数值的曲线拟合和实测验证确定。与现有方法相比，该方法可在保证测距精度的前提下，大大降低目标频率细化的计算量，保证了水位测量的实时性。
步骤三：计算目标W的角度参量(θW,ρW)：
Figure BDA0003399116580000124
Figure BDA0003399116580000125
其中λ0＝c/f0，为FMCW起始扫频频率f0对应的波长；
Figure BDA0003399116580000126
为接收天线RX1和RX3的回波信号频谱在谱线k4处的相位差，计算方法如下：
Figure BDA0003399116580000127
同理，
Figure BDA0003399116580000131
为接收天线RX2和RX3的回波信号频谱在谱线k4处的相位差，计算方法同上。
步骤四：计算目标W的径向相对速度vW：
根据步骤二求取的系数a，先求取目标W在频谱
Figure BDA0003399116580000132
和
Figure BDA0003399116580000133
中对应频率kW＝(k4-1+a)处的相角
Figure BDA0003399116580000134
和
Figure BDA0003399116580000135
Figure BDA0003399116580000136
Figure BDA0003399116580000137
再求目标W的径向相对速度vW：
Figure BDA0003399116580000138
其中
Figure BDA0003399116580000139
vW为正数时，表示目标W的运动方向靠近雷达水位计；vW为负数时，表示目标W的运动方向远离雷达水位计；需要注意的是，Tshape决定了目标的径向速度的测量范围：
Figure BDA00033991165800001310
可通过选择合适的Tshape值，在覆盖场景运动目标速度范围的条件下，保证测速的精度，以便后续的目标识别处理过程。
重复上述步骤一至步骤四过程，得到当前frame雷达数据检测到的四个目标(A1、A2、B和W)的属性值。
通过若干周期Tframe的测量，得到若干目标的测量结果(四项属性值)。按照以下方法完成水面目标识别：
步骤一：所有目标以距离这一属性值为条件，通过K均值聚类算法，分为若干组不同距离的目标类；
步骤二：对类内目标，计算距离均值
Figure BDA0003399116580000141
空间立体角θ和ρ的均方差θstd和ρstd，径向速度的均值
Figure BDA0003399116580000142
和均方差vstd；
步骤三：以空间立体角均方差为判定条件，将同时满足θstd≥Th1和ρstd≥Th1(Th1为设定的门限值)的目标类作为待定目标类；
步骤四：计算待定目标类径向速度均方差与均值之比，在满足
Figure BDA0003399116580000143
(Th2为设定的门限值)条件下，将距离均值
Figure BDA0003399116580000144
最小的目标类判定为水面目标。
结合目标的距离、回波强度、空间立体角和径向运动速度属性，将目标类分类，再从待定目标类中识别出水体目标，可排除干扰目标的影响，与现有技术相比，大大提高了水位测量的准确性。
将识别为水面目标类的距离均值
Figure BDA0003399116580000145
换算为雷达水位计距离水面的空高值h:
Figure BDA0003399116580000146
用设定的安装高度减去空高值得到当前水深值。
本发明的有益效果如下：
在雷达水位计已有的测量目标距离和回波强度的功能基础上，利用一发多收调频连续波体制的优势，快速解算出目标相对雷达水位计的距离、空间立体角和径向运动速度，在雷达水位计安装现场存在不可排除的多目标干扰的情况下，利用目标的距离、回波强度、空间立体角和径向运动速度等属性，准确识别出水体目标，排除干扰目标的影响，提高水位测量的准确性，降低对雷达水位计的安装要求。
本文所使用的词语“优选的”意指用作实例、示例或例证。本文描述为“优选的”任意方面或设计不必被解释为比其他方面或设计更有利。相反，词语“优选的”的使用旨在以具体方式提出概念。如本申请中所使用的术语“或”旨在意指包含的“或”而非排除的“或”。即，除非另外指定或从上下文中清楚，“X使用A或B”意指自然包括排列的任意一个。即，如果X使用A；X使用B；或X使用A和B二者，则“X使用A或B”在前述任一示例中得到满足。
而且，尽管已经相对于一个或实现方式示出并描述了本公开，但是本领域技术人员基于对本说明书和附图的阅读和理解将会想到等价变型和修改。本公开包括所有这样的修改和变型，并且仅由所附权利要求的范围限制。特别地关于由上述组件(例如元件等)执行的各种功能，用于描述这样的组件的术语旨在对应于执行所述组件的指定功能(例如其在功能上是等价的)的任意组件(除非另外指示)，即使在结构上与执行本文所示的本公开的示范性实现方式中的功能的公开结构不等同。此外，尽管本公开的特定特征已经相对于若干实现方式中的仅一个被公开，但是这种特征可以与如可以对给定或特定应用而言是期望和有利的其他实现方式的一个或其他特征组合。而且，就术语“包括”、“具有”、“含有”或其变形被用在具体实施方式或权利要求中而言，这样的术语旨在以与术语“包含”相似的方式包括。
本发明实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以多个或多个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。上述提到的存储介质可以是只读存储器，磁盘或光盘等。上述的各装置或系统，可以执行相应方法实施例中的存储方法。
综上所述，上述实施例为本发明的一种实施方式，但本发明的实施方式并不受所述实施例的限制，其他的任何背离本发明的精神实质与原理下所做的改变、修饰、代替、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (9)

1.一种多目标背景下的水位测量方法，使用雷达水位计测量水位，其特征在于，所述雷达水位计包括一个单阵元发射天线TX和三个单阵元接收天线RX1、RX2和RX3，所述水位测量方法包括以下步骤：
雷达水位计通过发射天线TX以固定周期Tframe发射一组调频连续波信号，每组调频连续波由第一锯齿波S1和第二锯齿波S2组成，两段锯齿波间隔Tshape时间发射；
三个接收天线接收采样后的第一锯齿波回波信号S1 RX1(n)，S1 RX2(n)，S1 RX3(n)和第二锯齿波回波信号S2 RX1(n)，S2 RX2(n)，S2 RX3(n)，n∈[1,N]；将S1 RX1(n)，S1 RX2(n)，S1 RX3(n)和S2 RX3(n)加N点汉明窗，做N点离散傅里叶变换后，得到离散傅里叶频谱
Figure FDA0003399116570000011
和
Figure FDA0003399116570000012
其中Mi(k)和
Figure FDA0003399116570000013
为第i个频谱在谱线位置k的频谱幅值和相角，k∈[1,N]；
在S1 RX3(n)的频谱
Figure FDA0003399116570000014
上用CFAR算法在k∈[1,N/2]范围内检测出局部幅值最大的m个目标的谱线位置k1、k2、...、km，求解所述m个目标的四元组属性值；
通过若干周期Tframe的测量，得到若干周期Tframe的所述m个目标的四元组属性值；
将所述若干目标分类，并根据预设规则选择类内的目标类为待定目标类，将距离均值最小的目标类判定为水面目标。
2.根据权利要求1所述的多目标背景下的水位测量方法，其特征在于，所述三个接收天线两两构成正交关系，且天线阵元中心间距Di相等，间距小于等于雷达发射信号波长的一半，定义三个接收天线的中心点为原点O，RX2和RX3中心点连线方向为Y轴，RX1和RX3中心点连线方向为X轴，建立参考直角坐标系。
3.根据权利要求2所述的多目标背景下的水位测量方法，其特征在于，所述目标的四元组属性为(r,θ,ρ,v)，若目标W在第一平面上的投影点为W'，在第二平面上的投影点为W"，OZ为雷达中心点到水平面的垂线，则θ为OW'与OZ的夹角，ρ为OW"与OZ的夹角，r为雷达中心点O到目标W的距离，v为目标W在OW径向上的相对速度，所述第一平面和第二平面互相垂直，并分别和水平面垂直。
4.根据权利要求2所述的多目标背景下的水位测量方法，其特征在于，单段锯齿波形状的调频连续波发射时起始频率为f0，结束频率为f1，带宽为BW，发射时间为Tp，且Tp<Tshape，其余时间段雷达不工作。
5.根据权利要求2所述的多目标背景下的水位测量方法，其特征在于，所述求解m个目标的四元组属性值包括以下步骤：
对于第i个目标，1<＝i<＝m，
步骤一：在S1 RX3(n)的频谱谱线ki的相邻两条谱线ki-1和ki+1中找到幅值最大的谱线，假设谱线幅值M3(ki)<M3(ki+1)，则以谱线ki和ki+1的频谱幅值M3(ki)和M3(ki+1)进行步骤二计算；
步骤二：计算目标i的距离ri为：
Figure FDA0003399116570000021
其中c为电磁波在真空中的传播速度，a为系数，Bw为带宽；
步骤三：计算目标i的角度参量(θi,ρi)：
Figure FDA0003399116570000022
Figure FDA0003399116570000023
其中λ0＝c/f0，为FMCW起始扫频频率f0对应的波长；
Figure FDA0003399116570000024
为接收天线RX1和RX3的回波信号频谱在谱线ki处的相位差，计算方法如下：
Figure FDA0003399116570000031
Figure FDA0003399116570000032
和
Figure FDA0003399116570000033
为第1个和第3个频谱在谱线位置ki的相角；
Figure FDA0003399116570000034
为接收天线RX2和RX3的回波信号频谱在谱线ki处的相位差，计算方法如下：
Figure FDA0003399116570000035
Figure FDA0003399116570000036
和
Figure FDA0003399116570000037
为第2个和第3个频谱在谱线位置ki的相角；
步骤四：计算目标i的径向相对速度vi：
根据所述系数a，先求取目标i在频谱
Figure FDA0003399116570000038
和
Figure FDA0003399116570000039
中对应频率k＝(ki-1+a)处的相角
Figure FDA00033991165700000310
和
Figure FDA00033991165700000311
Figure FDA00033991165700000312
Figure FDA00033991165700000313
Figure FDA00033991165700000314
和
Figure FDA00033991165700000315
为第3个和第4个频谱在谱线位置ki的相角，
Figure FDA00033991165700000316
和
Figure FDA00033991165700000317
为第3个和第4个频谱在谱线位置ki+1的相角；
再求目标i的径向相对速度vi：
Figure FDA00033991165700000318
其中
Figure FDA00033991165700000319
vi为正数时，表示目标i的运动方向靠近雷达水位计；vi为负数时，表示目标i的运动方向远离雷达水位计；Tshape是决定目标的径向速度的测量范围的参数；
重复上述步骤一至步骤四过程，得到当前frame雷达数据检测到的其余目标的四元组属性值。
6.根据权利要求5所述的多目标背景下的水位测量方法，其特征在于，所述系数a的计算公式为：
Figure FDA0003399116570000041
其中p0,p1,p2,p3为预设的系数，M3(ki)和M3(ki+1)为谱线ki和ki+1的幅值；
Tshape通过下式计算：
Figure FDA0003399116570000042
其中vmax为检测目标的最大径向速度。
7.根据权利要求1所述的多目标背景下的水位测量方法，其特征在于，按照以下方法完成水面目标识别：
步骤一：所有目标以距离这一属性值为条件，通过K均值聚类算法，分为若干组不同距离的目标类；
步骤二：对类内目标，计算距离均值
Figure FDA0003399116570000045
空间立体角θ和ρ的均方差θstd和ρstd，径向速度的均值
Figure FDA0003399116570000046
和均方差vstd；
步骤三：以空间立体角均方差为判定条件，将同时满足θstd≥Th1和ρstd≥Th1的目标类作为待定目标类，其中Th1为设定的门限值；
步骤四：计算待定目标类径向速度均方差与均值之比，在满足
Figure FDA0003399116570000043
条件下，将距离均值
Figure FDA0003399116570000044
最小的目标类判定为水面目标，其中Th2为设定的门限值。
8.根据权利要求1所述的多目标背景下的水位测量方法，其特征在于，识别出水面目标后，将识别为水面目标类的距离均值
Figure FDA0003399116570000051
换算为雷达水位计距离水面的空高值h:
Figure FDA0003399116570000052
其中
Figure FDA0003399116570000053
为水面目标类的θ均值，
Figure FDA0003399116570000054
为水面目标类的ρ均值。
9.根据权利要求8所述的多目标背景下的水位测量方法，其特征在于，
用设定的安装高度减去所述空高值得到当前水深值。
CN202111493031.1A 2021-12-08 2021-12-08 一种多目标背景下的水位测量方法 Pending CN114184256A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111493031.1A CN114184256A (zh) 2021-12-08 2021-12-08 一种多目标背景下的水位测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111493031.1A CN114184256A (zh) 2021-12-08 2021-12-08 一种多目标背景下的水位测量方法

Publications (1)

Publication Number Publication Date
CN114184256A true CN114184256A (zh) 2022-03-15

Family

ID=80542785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111493031.1A Pending CN114184256A (zh) 2021-12-08 2021-12-08 一种多目标背景下的水位测量方法

Country Status (1)

Country Link
CN (1) CN114184256A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117268496A (zh) * 2023-11-21 2023-12-22 无锡航征科技有限公司 基于雷达波的水位测量方法、装置及存储介质
CN117367541A (zh) * 2023-09-28 2024-01-09 深圳妙月科技有限公司 一种水位监测方法、装置、电子设备及存储介质

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117367541A (zh) * 2023-09-28 2024-01-09 深圳妙月科技有限公司 一种水位监测方法、装置、电子设备及存储介质
CN117268496A (zh) * 2023-11-21 2023-12-22 无锡航征科技有限公司 基于雷达波的水位测量方法、装置及存储介质
CN117268496B (zh) * 2023-11-21 2024-01-30 无锡航征科技有限公司 基于雷达波的水位测量方法、装置及存储介质

Similar Documents

Publication Publication Date Title
JP5460066B2 (ja) Ｖｈｆ／ｕｈｆレーダー局を使用して河川流量パラメータを監視するシステムおよび方法
CN109102702A (zh) 基于视频车辆检测和雷达信号融合的车辆测速方法
CN105445701B (zh) Ddma‑mimo雷达目标的单脉冲角度估计方法
CN114184256A (zh) 一种多目标背景下的水位测量方法
CN111006743B (zh) 一种基于平面雷达水位计的水位测量排除干扰的方法
CN110297233B (zh) Lfmcw阵列雷达信号并行流水化处理方法
CN107607936B (zh) 一种高频天地波雷达海洋表面流反演方法
CN114942428B (zh) 多目标跟踪的混合体制雷达数据关联算法
CN111025273B (zh) 一种畸变拖曳阵线谱特征增强方法及系统
CN111896913B (zh) 高频雷达单极子/交叉环天线通道增益校准方法及装置
CN110879391B (zh) 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法
CN109884337B (zh) 一种利用高频地波雷达探测海面风向的方法
CN115657007A (zh) 一种基于毫米波雷达的河流表面流速测量方法
Yang et al. Multibeam sonar bottom detection using multiple subarrays
CN110967677B (zh) 时频域分级多目标分辨与测量方法
CN113189554A (zh) 雷达实测回波数据的处理方法、电子设备及存储介质
CN111505590A (zh) 一种高频地波雷达通道校准方法及系统
CN116719023A (zh) Gpr能量聚焦和分辨率联合约束下的介质介电常数估计方法
CN110907937B (zh) 一种基于“t”型阵的掩埋物合成孔径三维成像方法
TW200819775A (en) Tracking target method of radar with frequency modulation continuous wave
CN104360335A (zh) 基于ais目标指示的电离层参数反演方法
JP3393183B2 (ja) ドップラ速度検出装置及び検出方法
CN115236664A (zh) 一种航海雷达图像反演有效波高的方法
CN113075633A (zh) 一种基于距离-俯仰联合字典的目标定位方法
CN113534130A (zh) 基于视线角度的多站雷达多目标数据关联方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination