CN115308435B - 一种双向测流的侧扫多普勒雷达流速仪 - Google Patents
一种双向测流的侧扫多普勒雷达流速仪 Download PDFInfo
- Publication number
- CN115308435B CN115308435B CN202211078678.2A CN202211078678A CN115308435B CN 115308435 B CN115308435 B CN 115308435B CN 202211078678 A CN202211078678 A CN 202211078678A CN 115308435 B CN115308435 B CN 115308435B
- Authority
- CN
- China
- Prior art keywords
- flow
- water
- water surface
- zero
- flow rate
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 344
- 238000005259 measurement Methods 0.000 claims abstract description 51
- 238000001914 filtration Methods 0.000 claims description 20
- 238000000691 measurement method Methods 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 238000003786 synthesis reaction Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000002156 mixing Methods 0.000 description 3
- 230000002457 bidirectional effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 239000003643 water by type Substances 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/24—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
- G01P5/241—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/581—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse modulated waves and based upon the Doppler effect resulting from movement of targets
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measuring Volume Flow (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请公开了一种双向测流的侧扫多普勒雷达流速仪，涉及水文技术领域，该侧扫多普勒雷达流速仪悬挂安装在水域岸边的支架上，向水域的表面不同位置发射两个电磁波波束并接收回波得到两个流速分量，并将两个流速分量投影到水面得到两个水面投影流速后基于卡尔曼跟踪算法进行滤波处理，对滤波处理后的两个水面投影流速进行矢量合成得到水域的流速测量结果。该侧扫多普勒雷达流速仪通过获取水域的水面上两个不同方向的水面投影流速进行合成，即便水域中的水流方向发生变化也能准确快速的得到流速测量结果，本申请结合卡尔曼跟踪算法进行滤波处理，可以降低天气等外界因素带来的误差，提高测量精度。
Description
技术领域
本申请涉及水文技术领域，尤其是一种双向测流的侧扫多普勒雷达流速仪。
背景技术
非接触式的多普勒雷达流速仪通过计算接收回波信号与发射电磁波信号的频率差值，基于多普勒原理可以得到液面流速，从而实现流速测量。多普勒雷达流速仪的操作简单，可在远端进行控制及通信，同时不受泥沙及污水腐蚀影响，使用寿命长，维护成本低，近年来已被广泛应用于水文水利等领域。
目前在使用多普勒雷达流速仪时，一般要求多普勒雷达流速仪要安装在水渠、河道、水库、湖泊等水域的中心线位置处，且要顺着水流方向发射脉冲并采集回波信号，才能实现流速测量的效果。但是实际应用环境比较复杂，受风向等其它因素影响，多普勒雷达流速仪安装好以后的工作过程中，水域中的水流方向并不固定，当水流方向发生变化时，多普勒雷达流速仪发射的脉冲不再顺着水流方向，就会导致测得的流速不准确。
发明内容
本申请人针对上述问题及技术需求，提出了一种双向测流的侧扫多普勒雷达流速仪，本申请的技术方案如下：
一种双向测流的侧扫多普勒雷达流速仪，该侧扫多普勒雷达流速仪悬挂安装在水域岸边的支架上，侧扫多普勒雷达流速仪执行的流速测量方法包括：
按照俯角向水域的表面不同位置发射两个电磁波波束，并分别接收两个电磁波波束经水面反射后返回的回波信号；
基于多普勒原理根据两个回波信号分别计算得到的两个流速分量，并将两个流速分量按照俯角投影到水面得到两个水面投影流速；
基于卡尔曼跟踪算法对两个水面投影流速进行滤波处理；
对滤波处理后的两个水面投影流速和/>进行矢量合成得到水域的流速测量结果。
其进一步的技术方案为，侧扫多普勒雷达流速仪在基于卡尔曼跟踪算法对两个流速分量进行滤波处理时：
对于每个流速分量，基于流速分量的当前时刻的测量值、前一时刻的预测值和前一时刻的测量值更新卡尔曼跟踪算法的***噪声的协方差和测量噪声的协方差，并利用更新后的卡尔曼跟踪算法对流速分量的当前时刻的测量值进行滤波处理。
其进一步的技术方案为，更新***噪声的协方差和测量噪声的协方差的方法包括按照下列公式更新：
；
其中，是更新后的测量噪声的协方差，/>是更新后的***噪声的协方差，是更新前的测量噪声的协方差，/>是更新前的***噪声的协方差；/>是补偿量权值；/>，/>是流速分量的当前时刻的测量值，/>是流速分量的前一时刻的测量值，/>是流速分量的前一时刻的预测值。
其进一步的技术方案为，对滤波处理后的两个水面投影流速和/>进行矢量合成得到水域的流速测量结果的方法包括：
根据滤波处理后的两个水面投影流速和/>的流动速率以及方向，通过矢量合成得到水域中水流的流动速率/>和流动方向，得到流速测量结果；其中，水域中水流的流动速率为/>，/>是根据滤波处理后的两个水面投影流速/>和/>的方向确定的两个水面投影流速之间的夹角。
其进一步的技术方案为，两个水面投影流速之间的夹角为：
当两个水面投影流速均沿着水流正方向或者均沿着水流负方向时，确定；
当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且两个水面投影流速与零度角方向的夹角之和超过180°时，确定；
当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且两个水面投影流速与零度角方向的夹角之和未超过180°时，确定；
其中，是水面投影流速/>与零度角方向的夹角，/>是水面投影流速/>与零度角方向的夹角，水流正方向从侧扫多普勒雷达流速仪所在一侧的岸边指向对向的另一侧岸边，水流负方向从对向的另一侧岸边指向侧扫多普勒雷达流速仪所在一侧的岸边；零度角方向平行于侧扫多普勒雷达流速仪所在位置的岸边。
其进一步的技术方案为，得到水域中水流的流动方向的方法包括：
当两个水面投影流速的方向均沿着水流正方向时，确定水域中水流的流动方向沿着水流正方向且与零度角方向的夹角为；当两个水面投影流速的方向均沿着水流负方向时，确定水域中水流的流动方向沿着水流负方向且与零度角方向的夹角为；其中，/>，/>表示两个水面投影流速中与零度角方向的夹角较小的一个水面投影流速，/>表示/>与零度角方向的夹角；/>表示两个水面投影流速中与零度角方向的夹角较大的一个水面投影流速，/>表示/>与零度角方向的夹角；
当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且两个水面投影流速与零度角方向的夹角之和超过180°时：当时确定水域中水流沿着水流正方向且与零度角方向的夹角为/>；当/>时确定水域中水流沿着水流负方向且与零度角方向的夹角为/>；
当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且两个水面投影流速与零度角方向的夹角之和未超过180°时：当时确定水域中水流沿着水流正方向且与零度角方向的夹角为/>；当/>时确定水域中水流沿着水流负方向且与零度角方向的夹角为/>；
其中，，/>是两个水面投影流速中沿着水流正方向的水面投影流速，/>是/>与零度角方向的夹角；/>是两个水面投影流速中沿着水流负方向的水面投影流速，/>是/>与零度角方向的夹角；
其中，水流正方向从侧扫多普勒雷达流速仪所在一侧的岸边指向对向的另一侧岸边，水流负方向从对向的另一侧岸边指向侧扫多普勒雷达流速仪所在一侧的岸边；零度角方向平行于侧扫多普勒雷达流速仪所在位置的岸边。
其进一步的技术方案为，侧扫多普勒雷达流速仪执行的流速测量方法还包括：
对矢量合成得到的水域的流速测量结果进行矢量分解，得到各个水面方向的流速分量，构建得到水域的二维水面流速剖面图，二维水面流速剖面图包含水域的水面任意方向上的流速分量。
其进一步的技术方案为，侧扫多普勒雷达流速仪包括主控芯片及其相连的两个射频收发芯片，侧扫多普勒雷达流速仪悬挂安装在水域岸边的支架上时，两个射频收发芯片的射频收发通道按照俯角分别朝向水域的表面不同位置；主控芯片分别通过两个射频收发芯片发射两个电磁波波束并分别接收对应的回波信号。
其进一步的技术方案为，俯角与支架的高度和水域的宽度相关，侧扫多普勒雷达流速仪发射的电磁波波束覆盖水域的有效区域。
其进一步的技术方案为，侧扫多普勒雷达流速仪执行流速测量方法得到水域的流速测量结果后进入休眠状态，直至下一个工作周期再执行流速测量方法。
本申请的有益技术效果是：
本申请公开了一种双向测流的侧扫多普勒雷达流速仪，该侧扫多普勒雷达流速仪通过获取水域的水面上两个不同方向的水面投影流速，合成得到流速测量结果，即便水域中的水流方向发生变化也能准确快速的得到流速测量结果，而且本申请结合自适应更新噪声的协方差的卡尔曼跟踪算法进行滤波处理，可以降低天气等外界因素带来的误差，提高测量精度。
相较于传统多普勒雷达流速仪只能获取固定方向上的流速的限制，本申请在得到流速测量结果后，可以计算得到各个水平水面方向上的流速分量，形成二维水面流速剖面图，实现全方向的流速测量。
该侧扫多普勒雷达流速仪只需悬挂安装在水域的岸边，而无需悬挂安装在水流中线上，因此降低了使用安装难度。
附图说明
图1是一个实施例中的双向测流的侧扫多普勒雷达流速仪的工作场景示意图。
图2是一个实施例中的侧扫多普勒雷达流速仪执行的流速测量方法的方法流程图。
图3中的（a）是一个实例中两个水面投影流速均沿着水流正方向的角度关系示意图，图3中的（b）是一个实例中两个水面投影流速均沿着水流负方向的角度关系示意图。
图4中的（a）和（b）是水面投影流速沿着水流负方向、水面投影流速/>沿着水流正方向，且/>的两个实例中的角度关系示意图。
图5中的（a）和（b）是水面投影流速沿着水流正方向、水面投影流速/>沿着水流负方向，且/>的两个实例中的角度关系示意图。
图6是一个实例中根据流速测量结果构建得到的水域的二维水面流速剖面图的示意图。
具体实施方式
下面结合附图对本申请的具体实施方式做进一步说明。
本申请公开了一种双向测流的侧扫多普勒雷达流速仪，请参考图1，该侧扫多普勒雷达流速仪1悬挂安装在水域3岸边的支架2上。侧扫多普勒雷达流速仪1内部包括主控芯片及其相连的两个射频收发芯片，当该侧扫多普勒雷达流速仪1悬挂安装在水域3岸边的支架2上时，两个射频收发芯片的射频收发通道按照俯角分别朝向水域3的表面不同位置，主控芯片分别通过两个射频收发芯片发射两个电磁波波束并分别接收对应的回波信号。如图1所示，侧扫多普勒雷达流速仪1通过两个射频收发芯片发射电磁波波束4和电磁波波束5。
安装侧扫多普勒雷达流速仪1时，保证两个射频收发芯片的射频收发通道按照俯角朝下照射水域3的水面，与水平面之间的俯角/>与支架2的高度H和水域3的宽度W相关。使得侧扫多普勒雷达流速仪1发射的电磁波波束覆盖水域3的有效区域，有效区域可以预先设定，通过调整俯角/>可以使得波束照射到对岸的位置处，覆盖面积较广。常规方法将多普勒雷达流速仪安装在水域3中心线位置时，若水域3较宽，则需要增加悬挂雷达的杆臂的长度，而本申请不会有这个问题，只需将侧扫多普勒雷达流速仪安装在岸边，可以适应不同宽度的水域3。
该侧扫多普勒雷达流速仪1执行的流速测量方法包括如下步骤，该方法由该侧扫多普勒雷达流速仪1内部的主控芯片完成，请参考图2所示的流程图：
步骤210，按照俯角向水域3的表面不同位置发射两个电磁波波束4和5，并分别接收两个电磁波波束经水面反射后返回的回波信号。
步骤220，基于多普勒原理根据两个回波信号分别计算得到的两个流速分量。对于接收到的每个回波信号，将回波信号与本振信号进行混频，混频后经过低通滤波放大电路滤除频率为的高频信号，得到频率为/>的频偏信号写为/>，其中，/>为多普勒频移角频率，/>为本振信号角频率，/>和/>为参数。对频偏信号进行采样，将频偏信号转换成离散数字形式的离散信号，去除离散信号的直流分量，经过傅里叶变换计算得到多普勒频率，然后根据已知的频率-速度转换公式即可确定多普勒频率对应的流速分量。
步骤230，将两个流速分量按照俯角投影到水面得到两个水面投影流速，得到的水面投影流速包括流动速率和方向。
两个水面投影流速的流动速率为，/>，其中，/>和/>是计算得到的两个电磁波波束对应的流速分量。
水面投影流速的方向有两种情况：一种是从射频收发芯片在水面的垂直投影位置指向电磁波波束在水面的照射位置，比如图1中水面投影流速的方向可以从A指向B，水面投影流速/>的方向可以从C指向D。另一种是从电磁波波束在水面的照射位置指向射频收发芯片在水面的垂直投影位置，比如图1中水面投影流速/>的方向可以从B指向A，水面投影流速/>的方向可以从D指向C。
本申请定义从侧扫多普勒雷达流速仪1所在侧的岸边指向另一侧岸边的方向为零度角方向，从另一侧岸边指向侧扫多普勒雷达流速仪1所在侧的岸边的方向为水流负方向，并将侧扫多普勒雷达流速仪1安装位置处的河岸区域看作是一条直线，预先定义一个平行于岸边的零度角方向，比如图1和图3所示的零度角方向。在本申请中，由于对接收到的回波信号进行了正交混频处理，对每个波束的两路IQ信号都进行采样，因此在得到两个水面投影流速可以得到其矢量角，该矢量角指示水面投影流速的方向，实际水面投影流速的矢量角指示了水面投影流速与零度角方向之间的夹角，并指示了水面投影流速沿着水流零度角方向还是水流负方向。
比如对于图3中的（a），水面投影流速的矢量角指示其与零度角方向的夹角为且沿着水流零度角方向，比如实际应用时，水面投影流速/>的矢量角为+30°度，表示其与零度角方向的夹角为30°且沿着水流零度角方向。而对于图3中的（b），水面投影流速/>的矢量角指示其与零度角方向的夹角为/>且沿着水流负方向，比如实际应用时，水面投影流速/>的矢量角为-150°，表示其与零度角方向的夹角为150°且沿着水流负方向。
步骤240，基于卡尔曼跟踪算法对两个水面投影流速进行滤波处理。
由于侧扫多普勒雷达流速仪测量的是水域3表面的流速，因此容易受风、雨等外界天气影响，导致测量结果出现异常跳动，每一个水面投影流速的异常都会导致最终得到的流速测量结果的大小及方向不准确。因此为了减小外界因素造成的测量误差，本申请通过卡尔曼跟踪算法对得到的两个水面投影流速进行滤波处理。
常规的卡尔曼跟踪算法中的***噪声的协方差和测量噪声的协方差通常为固定的预设值，这并不适合本申请这种存在较多异常数据的场景，因此在一个实施例中，侧扫多普勒雷达流速仪在基于卡尔曼跟踪算法对两个流速分量进行滤波处理时：对于每个流速分量，基于流速分量的当前时刻的测量值、前一时刻的预测值/>和前一时刻的测量值/>更新卡尔曼跟踪算法的***噪声的协方差和测量噪声的协方差。包括按照下列公式更新：
；
其中，是更新后的测量噪声的协方差，/>是更新后的***噪声的协方差，是更新前的测量噪声的协方差，/>是更新前的***噪声的协方差。/>是补偿量权值，可以根据实际使用过程中对异常值的容忍度来进行调整。/>。
然后利用更新后的卡尔曼跟踪算法对流速分量的当前时刻的测量值进行滤波处理。如此卡尔曼跟踪算法的***噪声的协方差和测量噪声的协方差不断动态更新，既保证平稳过程数据预测值对测量值的依赖占比，保证数据真实性，又可以有效的抑制异常的测量值。
步骤250，对滤波处理后的两个水面投影流速和/>进行矢量合成得到水域3的流速测量结果。根据两个水面投影流速/>和/>的流动速率以及方向，通过矢量合成得到水域中水流的流动速率/>和流动方向，也即得到的流速测量结果包括水流的流动速率/>和流动方向。
根据水面投影流速和水面投影流速/>的矢量角可以得到水面投影流速/>和水面投影流速/>之间的夹角/>，假设水面投影流速/>与零度角方向的夹角为/>，水面投影流速/>与零度角方向的夹角为/>，/>，/>，计算得到的/>：
（1）当两个水面投影流速均沿着水流正方向或者均沿着水流负方向时，确定，如图3中（a）表示两个水面投影流速均沿着水流正方向，图3中的（b）表示两个水面投影流速均沿着水流负方向。
（2）当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且时，确定/>。如图4中的（a）和（b）均以水面投影流速/>沿着水流负方向、水面投影流速/>沿着水流正方向为例。
（3）当两个水面投影流速的方向中一个为水流正方向、另一个为水流负方向，且时，确定/>。如图5中的（a）和（b）均以水面投影流速/>沿着水流正方向、水面投影流速/>沿着水流负方向为例。
得到两个水面投影流速的夹角后，即可确定水域中水流的流动速率为。
确定水域中水流的流动方向的方法为：
（1）当两个水面投影流速的方向均沿着水流正方向时，确定水域中水流的流动方向也沿着水流正方向且与零度角方向的夹角为，/>。
当两个水面投影流速的方向均沿着水流负方向时，确定水域中水流的流动方向也沿着水流负方向且与零度角方向的夹角为，/>。
其中，表示两个水面投影流速中与零度角方向的夹角较小的一个水面投影流速，/>表示/>与零度角方向的夹角。/>表示两个水面投影流速中与零度角方向的夹角较大的一个水面投影流速，/>表示/>与零度角方向的夹角。也即当/>为/>时/>为，当/>为/>时/>为/>。当/>为/>时/>为/>，当/>为/>时/>为/>。
如图3中的（a）以两个水面投影流速的方向均沿着水流正方向，且为/>、/>为/>，以及/>为/>、/>为/>为例。图3中的（b）以两个水面投影流速的方向均沿着水流负方向，且/>为/>、/>为/>，以及/>为/>、/>为/>为例。。
（2）当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且时：
当时确定水域中水流沿着水流正方向且与零度角方向的夹角为/>。
当时确定水域中水流沿着水流负方向且与零度角方向的夹角为。
其中，，/>是两个水面投影流速中沿着水流正方向的水面投影流速，/>是/>与零度角方向的夹角。/>是两个水面投影流速中沿着水流负方向的水面投影流速，/>是/>与零度角方向的夹角。也即当/>为/>时/>为/>，当/>为/>时/>为/>。当/>为/>时/>为/>，当/>为/>时/>为/>。
如图4中的（a）以为/>、/>为/>，且/>为/>、/>为/>，且/>的情况为例。如图4中的（b）以/>为/>、/>为/>，且/>为/>、/>为/>，且/>的情况为例。
（3）当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且时：
当时确定水域中水流沿着水流正方向且与零度角方向的夹角为/>。
当时确定水域中水流沿着水流负方向且与零度角方向的夹角为/>。
其中，，其余定义如上述情况（2）。
如图5中的（a）以为/>、/>为/>，且/>为/>、/>为/>，且/>的情况为例。如图5中的（b）以/>为/>、/>为/>，且/>为/>、/>为/>，且/>的情况为例。
由此可以得到水域的流速测量结果，该侧扫多普勒雷达流速仪在实际应用时一般按照周期性执行流速测量，则该侧扫多普勒雷达流速仪执行流速测量方法得到水域的流速测量结果后进入休眠状态以降低功耗，直至执行下一个工作周期中唤醒并按上述流程再次执行流速测量方法。
在得到水域的流速测量结果后，还可以对矢量合成得到的水域的流速测量结果进行矢量分解，由此得到各个水面方向的流速分量，构建得到水域的二维水面流速剖面图，二维水面流速剖面图包含水域的水面任意方向上的流速分量如图6示出了其中三个不同的水面方向上由流速测量结果矢量分解得到的流速分量/>以及其构成的圆形的二维水面流速剖面图。
以上所述的仅是本申请的优选实施方式，本申请不限于以上实施例。可以理解，本领域技术人员在不脱离本申请的精神和构思的前提下直接导出或联想到的其他改进和变化，均应认为包含在本申请的保护范围之内。
Claims (10)
1.一种用双向测流的侧扫多普勒雷达流速仪执行的流速测量方法，其特征在于，所述侧扫多普勒雷达流速仪悬挂安装在水域岸边的支架上，所述流速测量方法包括：
按照俯角向水域的表面不同位置发射两个电磁波波束，并分别接收两个电磁波波束经水面反射后返回的回波信号；
基于多普勒原理根据两个回波信号分别计算得到的两个流速分量，并将两个流速分量按照俯角投影到水面得到两个水面投影流速；
基于卡尔曼跟踪算法对两个水面投影流速进行滤波处理；
对滤波处理后的两个水面投影流速和/>进行矢量合成得到水域的流速测量结果。
2.根据权利要求1所述的流速测量方法，其特征在于，所述侧扫多普勒雷达流速仪在基于卡尔曼跟踪算法对两个流速分量进行滤波处理时：
对于每个流速分量，基于流速分量的当前时刻的测量值、前一时刻的预测值和前一时刻的测量值更新卡尔曼跟踪算法的***噪声的协方差和测量噪声的协方差，并利用更新后的卡尔曼跟踪算法对流速分量的当前时刻的测量值进行滤波处理。
3.根据权利要求2所述的流速测量方法，其特征在于，更新***噪声的协方差和测量噪声的协方差的方法包括按照下列公式更新：
；
其中，是更新后的测量噪声的协方差，/>是更新后的***噪声的协方差，/>是更新前的测量噪声的协方差，/>是更新前的***噪声的协方差；/>是补偿量权值；，/>是流速分量的当前时刻的测量值，/>是流速分量的前一时刻的测量值，/>是流速分量的前一时刻的预测值。
4.根据权利要求1所述的流速测量方法，其特征在于，对滤波处理后的两个水面投影流速和/>进行矢量合成得到水域的流速测量结果的方法包括：
根据滤波处理后的两个水面投影流速和/>的流动速率以及方向，通过矢量合成得到水域中水流的流动速率/>和流动方向，得到流速测量结果；其中，水域中水流的流动速率为/>，/>是根据滤波处理后的两个水面投影流速/>和/>的方向确定的两个水面投影流速之间的夹角。
5.根据权利要求4所述的流速测量方法，其特征在于，两个水面投影流速之间的夹角为：
当两个水面投影流速均沿着水流正方向或者均沿着水流负方向时，确定；
当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且两个水面投影流速与零度角方向的夹角之和超过180°时，确定；
当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且两个水面投影流速与零度角方向的夹角之和未超过180°时，确定；
其中，是水面投影流速/>与零度角方向的夹角，/>是水面投影流速/>与零度角方向的夹角，水流正方向从所述侧扫多普勒雷达流速仪所在一侧的岸边指向对向的另一侧岸边，水流负方向从对向的另一侧岸边指向所述侧扫多普勒雷达流速仪所在一侧的岸边；零度角方向平行于所述侧扫多普勒雷达流速仪所在位置的岸边。
6.根据权利要求4所述的流速测量方法，其特征在于，得到水域中水流的流动方向的方法包括：
当两个水面投影流速的方向均沿着水流正方向时，确定水域中水流的流动方向沿着水流正方向且与零度角方向的夹角为；当两个水面投影流速的方向均沿着水流负方向时，确定水域中水流的流动方向沿着水流负方向且与零度角方向的夹角为/>；其中，/>，/>表示两个水面投影流速中与零度角方向的夹角较小的一个水面投影流速，/>表示/>与零度角方向的夹角；/>表示两个水面投影流速中与零度角方向的夹角较大的一个水面投影流速，/>表示/>与零度角方向的夹角；
当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且两个水面投影流速与零度角方向的夹角之和超过180°时：当时确定水域中水流沿着水流正方向且与零度角方向的夹角为/>；当/>时确定水域中水流沿着水流负方向且与零度角方向的夹角为/>；
当两个水面投影流速中一个沿着水流正方向、另一个沿着水流负方向，且两个水面投影流速与零度角方向的夹角之和未超过180°时：当时确定水域中水流沿着水流正方向且与零度角方向的夹角为/>；当/>时确定水域中水流沿着水流负方向且与零度角方向的夹角为/>；
其中，，/>是两个水面投影流速中沿着水流正方向的水面投影流速，/>是/>与零度角方向的夹角；/>是两个水面投影流速中沿着水流负方向的水面投影流速，/>是/>与零度角方向的夹角；
其中，水流正方向从所述侧扫多普勒雷达流速仪所在一侧的岸边指向对向的另一侧岸边，水流负方向从对向的另一侧岸边指向所述侧扫多普勒雷达流速仪所在一侧的岸边；零度角方向平行于所述侧扫多普勒雷达流速仪所在位置的岸边。
7.根据权利要求1所述的流速测量方法，其特征在于，所述流速测量方法还包括：
对矢量合成得到的水域的流速测量结果进行矢量分解，得到各个水面方向的流速分量，构建得到水域的二维水面流速剖面图，所述二维水面流速剖面图包含水域的水面任意方向上的流速分量。
8.根据权利要求1所述的流速测量方法，其特征在于，所述侧扫多普勒雷达流速仪包括主控芯片及其相连的两个射频收发芯片，所述侧扫多普勒雷达流速仪悬挂安装在水域岸边的支架上时，两个射频收发芯片的射频收发通道按照俯角分别朝向水域的表面不同位置；所述主控芯片分别通过两个射频收发芯片发射两个电磁波波束并分别接收对应的回波信号。
9.根据权利要求1所述的流速测量方法，其特征在于，俯角与支架的高度和水域的宽度相关，所述侧扫多普勒雷达流速仪发射的电磁波波束覆盖水域的有效区域。
10.根据权利要求1所述的流速测量方法，其特征在于，所述侧扫多普勒雷达流速仪执行的流速测量方法得到水域的流速测量结果后进入休眠状态，直至下一个工作周期再执行流速测量方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211078678.2A CN115308435B (zh) | 2022-09-05 | 2022-09-05 | 一种双向测流的侧扫多普勒雷达流速仪 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211078678.2A CN115308435B (zh) | 2022-09-05 | 2022-09-05 | 一种双向测流的侧扫多普勒雷达流速仪 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115308435A CN115308435A (zh) | 2022-11-08 |
| CN115308435B true CN115308435B (zh) | 2024-04-16 |
Family
ID=83867377
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211078678.2A Active CN115308435B (zh) | 2022-09-05 | 2022-09-05 | 一种双向测流的侧扫多普勒雷达流速仪 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115308435B (zh) |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104076353A (zh) * | 2014-06-30 | 2014-10-01 | 西安空间无线电技术研究所 | 一种面目标回波波束中心速度测量方法 |
| JP2015025699A (ja) * | 2013-07-25 | 2015-02-05 | 株式会社東京建設コンサルタント | ドップラー水象レーダー |
| CN104502910A (zh) * | 2015-01-13 | 2015-04-08 | 上海航征测控***有限公司 | 一种用于测量河道表面流速的手持式雷达 |
| CN109001723A (zh) * | 2018-06-29 | 2018-12-14 | 薛志勇 | 雷达流速仪 |
| CN109143223A (zh) * | 2018-08-14 | 2019-01-04 | 中国电子科技集团公司第三十八研究所 | 一种双基地雷达的空间目标跟踪滤波装置及方法 |
| CN109596856A (zh) * | 2018-12-25 | 2019-04-09 | 北京奥特美克科技股份有限公司 | 应用在雷达波流速仪的流速测量方法及雷达波流速仪 |
| CN110221299A (zh) * | 2019-07-04 | 2019-09-10 | 南京信息工程大学 | 一种星载双通道双波束InSAR测流*** |
| CN110568215A (zh) * | 2019-09-30 | 2019-12-13 | 南京微通电子技术有限公司 | 一种雷达流水计及其测量方法 |
| CN110672875A (zh) * | 2019-10-21 | 2020-01-10 | 中国科学院电子学研究所 | 基于Chirp-Z变换的表面水流速度检测方法 |
| CN110673143A (zh) * | 2019-09-30 | 2020-01-10 | 西安电子科技大学 | 一种子孔径大斜视sar俯冲成像的两步处理方法 |
| CN112630772A (zh) * | 2021-01-15 | 2021-04-09 | 陕西长岭电子科技有限责任公司 | 侧向两波束车载多普勒测速雷达设备 |
| CN113108767A (zh) * | 2021-04-07 | 2021-07-13 | 王陶然 | 一种无人机载雷达水文信息实时监测方法 |
| CN113933827A (zh) * | 2021-10-18 | 2022-01-14 | 无锡市雷华科技有限公司 | 一种非接触式船载测速*** |
| CN114778885A (zh) * | 2022-04-28 | 2022-07-22 | 珠海微度芯创科技有限责任公司 | 一种雷达流速仪及控制方法、装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10613212B2 (en) * | 2017-08-14 | 2020-04-07 | Oculii Corp. | Systems and methods for doppler-enhanced radar tracking |
| FR3090893B1 (fr) * | 2018-12-20 | 2021-10-08 | Commissariat Energie Atomique | Procede et systeme de mesure de la vitesse d’un porteur par rapport au sol |
-
2022
- 2022-09-05 CN CN202211078678.2A patent/CN115308435B/zh active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015025699A (ja) * | 2013-07-25 | 2015-02-05 | 株式会社東京建設コンサルタント | ドップラー水象レーダー |
| CN104076353A (zh) * | 2014-06-30 | 2014-10-01 | 西安空间无线电技术研究所 | 一种面目标回波波束中心速度测量方法 |
| CN104502910A (zh) * | 2015-01-13 | 2015-04-08 | 上海航征测控***有限公司 | 一种用于测量河道表面流速的手持式雷达 |
| CN109001723A (zh) * | 2018-06-29 | 2018-12-14 | 薛志勇 | 雷达流速仪 |
| CN109143223A (zh) * | 2018-08-14 | 2019-01-04 | 中国电子科技集团公司第三十八研究所 | 一种双基地雷达的空间目标跟踪滤波装置及方法 |
| CN109596856A (zh) * | 2018-12-25 | 2019-04-09 | 北京奥特美克科技股份有限公司 | 应用在雷达波流速仪的流速测量方法及雷达波流速仪 |
| CN110221299A (zh) * | 2019-07-04 | 2019-09-10 | 南京信息工程大学 | 一种星载双通道双波束InSAR测流*** |
| CN110568215A (zh) * | 2019-09-30 | 2019-12-13 | 南京微通电子技术有限公司 | 一种雷达流水计及其测量方法 |
| CN110673143A (zh) * | 2019-09-30 | 2020-01-10 | 西安电子科技大学 | 一种子孔径大斜视sar俯冲成像的两步处理方法 |
| CN110672875A (zh) * | 2019-10-21 | 2020-01-10 | 中国科学院电子学研究所 | 基于Chirp-Z变换的表面水流速度检测方法 |
| CN112630772A (zh) * | 2021-01-15 | 2021-04-09 | 陕西长岭电子科技有限责任公司 | 侧向两波束车载多普勒测速雷达设备 |
| CN113108767A (zh) * | 2021-04-07 | 2021-07-13 | 王陶然 | 一种无人机载雷达水文信息实时监测方法 |
| CN113933827A (zh) * | 2021-10-18 | 2022-01-14 | 无锡市雷华科技有限公司 | 一种非接触式船载测速*** |
| CN114778885A (zh) * | 2022-04-28 | 2022-07-22 | 珠海微度芯创科技有限责任公司 | 一种雷达流速仪及控制方法、装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115308435A (zh) | 2022-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN201262647Y (zh) | 利用vhf/uhf雷达站监测河流流速参数的*** | |
| CN106405556B (zh) | 车辆目标信息探测识别***及其信号处理方法 | |
| CN101872020B (zh) | 基于频谱变换的水下运动目标定位导航方法和装置 | |
| EP3460417B1 (en) | Single-beam radar level and velocity sensing | |
| KR101631967B1 (ko) | 장주기 파고측정용 gnss 기반 부이에서의 수신기 위상말림 보상을 위한 ppp 및 마그네틱 컴파스 통합시스템 | |
| CN113302459B (zh) | 非侵入式敞开通道流量计 | |
| CN113108767A (zh) | 一种无人机载雷达水文信息实时监测方法 | |
| CN111766599B (zh) | 一种实时测量海浪高度的激光雷达装置及方法 | |
| CN104007286A (zh) | 一种用于河道表面流速的自动监测*** | |
| KR102555054B1 (ko) | 멀티포인트 레이더 유속계로부터 획득된 정보와 다양한 기상 관측 정보를 빅데이터화 및 데이터 마이닝하여 하천 유량을 측정하는 브레이니 rtu 및 상기 멀티포인트 레이더 유속계와 상기 브레이니 rtu를 포함하는 하천 유량 측정 시스템 | |
| CN115657007A (zh) | 一种基于毫米波雷达的河流表面流速测量方法 | |
| CN109946659A (zh) | 一种车载毫米波雷达线性调频连续波运动频率扩展校正方法 | |
| CN115308435B (zh) | 一种双向测流的侧扫多普勒雷达流速仪 | |
| CN116087908A (zh) | 一种基于协同作业的雷达高精度物位计测量方法 | |
| CN114184256A (zh) | 一种多目标背景下的水位测量方法 | |
| CN206019788U (zh) | 一种微波雷达波浪水位测量装置 | |
| JP2015111097A (ja) | 水際から離隔した位置からの水面挙動の計測方法 | |
| WO2014192528A1 (ja) | 表層潮流推定装置、レーダ装置、表層潮流推定方法、及び表層潮流推定プログラム | |
| CN112698319A (zh) | 一种雷达测量目标角度的实验方法 | |
| CN113805182B (zh) | 一种基于矢量运算的物体运动监测方法 | |
| CN110673128B (zh) | 一种基于间断上下调频波的x波段岸基雷达测流方法 | |
| CN112630772A (zh) | 侧向两波束车载多普勒测速雷达设备 | |
| CN110542902B (zh) | 一种适用于车载毫米波雷达的简易测高方法 | |
| CN109959907B (zh) | 一种对高度数据进行处理的方法 | |
| CN109959928B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |