CN116125459A - 一种用于侧扫雷达的有效测速单元判定方法 - Google Patents

一种用于侧扫雷达的有效测速单元判定方法 Download PDF

Info

Publication number
CN116125459A
CN116125459A CN202310053518.0A CN202310053518A CN116125459A CN 116125459 A CN116125459 A CN 116125459A CN 202310053518 A CN202310053518 A CN 202310053518A CN 116125459 A CN116125459 A CN 116125459A
Authority
CN
China
Prior art keywords
radar
unit
flow velocity
flow
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.)
Granted
Application number
CN202310053518.0A
Other languages
English (en)
Other versions
CN116125459B (zh
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.)
Jiangsu Naiwch Cooperation
Nanjing Water Conservancy and Hydrology Automatization Institute Ministry of Water Resources
Original Assignee
Nanjing Water Conservancy and Hydrology Automatization Institute Ministry of Water Resources
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 Nanjing Water Conservancy and Hydrology Automatization Institute Ministry of Water Resources filed Critical Nanjing Water Conservancy and Hydrology Automatization Institute Ministry of Water Resources
Priority to CN202310053518.0A priority Critical patent/CN116125459B/zh
Publication of CN116125459A publication Critical patent/CN116125459A/zh
Application granted granted Critical
Publication of CN116125459B publication Critical patent/CN116125459B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
    • 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)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种用于侧扫雷达的有效测速单元判定方法，包括：获取流速单元数据；判断动态流速单元有效性，得到各流速单元的位置指数；判定水面信号质量，获取信号强度对于测流准确度的影响程度；实时判定各单元流速有效性，获取对应的流速偏离指数；结合得到的多个因子判断侧扫雷达单元流速有效性。本发明通过对测速单元进行多方面判定，量化其有效性，从而能够去除侧扫雷达干扰单元数据，提高了用于最终计算的流速样本准确性，使得侧扫雷达系统流量计算精度大大提高。

Description

一种用于侧扫雷达的有效测速单元判定方法
技术领域
本发明属于测量技术领域，涉及水文测量技术，具体涉及一种用于侧扫雷达的有效测速单元判定方法。
背景技术
侧扫雷达是进行大江大河、跨界河流等河流表面流速测量的水文雷达，采用非接触式雷达技术，实现对河流表面流场、网格点流速进行连续监测，同时提供数据上传服务；系统通过水位、过流面积、表面流速等数据交互，完成流量数据网络合成，实现全天候、连续自动河流流量监测。在水文测流行业中，近年来侧扫雷达测流技术应用越来越广泛。相对于必须安装在水流正上方的传统点雷达，侧扫雷达可采集水流表面多个单元的流速，但实际使用中我们却发现其监测结果误差较大，输出结果精确度较低。究其原因，我们认为并不是所有侧扫雷达采集的单元的流速样本均可参与最终的流量计算，但究竟其中哪些是对精确计算有正面效果，哪些是应被剔除的，现有技术中并没有精确的判定方法。
发明内容
针对现有技术中存在的问题，本发明公开了一种用于侧扫雷达的有效测速单元判定方法，有机结合了动态测速单元、水面信号强度、实时流速等多个因子，形成一套有效的测速单元有效性的判定方法，剔除了干扰单元，保留对精确计算有正面效果的流速样本。
为达到上述目的，本发明的技术方案如下：
一种用于侧扫雷达的有效测速单元判定方法，包括如下步骤：
步骤1，获取流速单元数据
侧扫雷达设备通过测量计算，生成若干沿测流断面方向的动态数量的扇形水面流速单元；
步骤2，通过各流速单元的位置指数，判定动态流速单元有效性；
步骤3，通过水面信号质量，判定动态流速单元有效性；
步骤4，获取各流速单元的流速偏离指数，判定实时流速有效性；
步骤5，结合步骤2-4中得到的多个因子判断侧扫雷达单元流速有效性。
进一步的，所述步骤2包括如下子步骤：
（1）收集河道监测站的安装信息；
（2）根据测站安装信息，构建测站三维空间模型如下：
                       （1）
式中，为流速单元起点距，为雷达安装起点距，为雷达安装高程，为实时水位值，为雷达发射天线至表面流速单元位置的直线距离；
（3）为三维空间模型接入实时水位，计算第i个流速单元的位置指数：
                 （2）
式中，为测流断面近雷达边界处起点距，为测流断面远雷达边界处起点距，由水位和测流断面的大断面资料可求得；为直线的水面投影，为流速单元近雷达边界与测流断面之间的垂直距离；为第i个流速单元处水深，为雷达测速所需的最小水深。
进一步的，的计算方法如下：
。
进一步的，所述步骤2中的安装信息至少包括：雷达安装起点距、雷达安装高程、大断面信息。
进一步的，所述步骤3包括如下子步骤：
（1）根据侧扫雷达的测速单元布设信息，在侧扫雷达工作的断面采用历史资料或比测仪器多次测量流速信息；
（2）采集侧扫雷达相应单元格的信号质量信息；
（3）将同一时刻同一起点距单元格流速与侧扫雷达单元流速对比，分析出能准确反映流场状态的。
进一步的，所述步骤3中，通过侧扫雷达设备采集处理海量信号得出，其具体值由雷达设备在一次测量完成后发出。
进一步的，所述步骤4包括如下子步骤：
（1）根据侧扫雷达的测速单元布设信息，在侧扫雷达工作的断面采用采用历史资料或流速比测仪器多次测量流速，率定出第i个流速单元处的水位-水面流速关系；
（2）将同一时刻同一起点距流速比测仪器得出的流速与侧扫雷达单元流速对比，分析侧扫雷达自然工况中的测速有效性；
（3）根据率定出的第i个流速单元处的水位-水面流速关系，计算当前水位第i个流速单元的流速值偏离指数：
               （4）
式中，为当前水位第i个流速单元的流速值，为侧扫雷达设备所能测到的最大流速，为侧扫雷达设备所能测到的最小流速，为由断面测流资料率定出的第i个流速单元处的水位-水面流速关系。
进一步的，所述步骤4中，若测量断面无率定成果，将直接取值。
进一步的，所述步骤5包括如下过程：
通过下列公式计算每个单元的流速值有效性：
                                   （5）
式中，为第i个流速单元的有效性指数，取值范围0-1；为第i个流速单元的位置指数，取值0或1；为第i个流速单元的信号可信度指数，取值范围0-1；为第i个流速单元的流速值偏离指数。
本发明的有益效果为：
本发明通过对测速单元进行多方面判定，量化其有效性，从而能够去除侧扫雷达干扰单元数据，提高了用于最终计算的流速样本准确性，使得侧扫雷达系统流量计算精度大大提高。
附图说明
图1为测站空间示意图。
图2为测站空间俯视图。
具体实施方式
以下将结合具体实施例对本发明提供的技术方案进行详细说明，应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
本发明提供了一种用于侧扫雷达的有效测速单元判定方法，包括如下步骤：
步骤1，获取流速单元数据
侧扫雷达设备通过测量计算，生成一系列沿测流横断面方向的动态数量的扇形水面流速单元；
步骤2，判定动态流速单元有效性，包括以下子步骤：
（1）收集河道监测站的安装信息，包括雷达安装起点距、雷达安装高程、大断面信息等。
（2）根据测站安装信息，构建测站三维空间模型（如图1、图2所示），为后续计算提供支撑：
                         （1）
假设起点距的零点位于雷达安装岸边，起点距自雷达安装岸边向过水断面另一侧岸边逐渐增大，式中，
为流速单元起点距，m；
为雷达安装起点距，m；
为雷达安装高程，m；
为实时水位值，m；
为雷达发射天线至表面流速单元位置的直线距离，由雷达波在空气中的传播速度及时间算得，m；
（3）为三维空间模型接入实时水位（由现场设立的水位计采集得来），计算第i个流速单元的位置指数作为判定测速单元是否有效的一个因子，计算公式如下：
                 （2）
为测流断面近雷达边界处起点距，由水位和测流断面的大断面资料可求得，m；
为测流断面远雷达边界处起点距，由水位和测流断面的大断面资料可求得，m；
为第i个流速单元处水深，由水位和测流断面的大断面资料可求得，m；
为雷达测速所需的最小水深，为侧扫雷达一项技术参数，m；
为直线的水面投影；；
为流速单元近雷达边界与测流断面之间的垂直距离。
步骤3，判定水面信号质量，包括如下子步骤：
（1）根据侧扫雷达的测速单元布设信息，在侧扫雷达工作的断面采用走航式ADCP（也可以采用其他流速比测仪器，如流速仪等）多次测量流速信息。此外，也可根据实际需要采用历史资料中的流速数据。
（2）采集侧扫雷达相应单元格的信号质量信息。
（3）将同一时刻同一起点距ADCP第一层单元格流速与侧扫雷达单元流速对比，分析出能准确反映流场状态的，作为判定侧扫雷达测速单元是否有效的一个因子。
通过侧扫雷达设备采集处理海量信号得出，其具体值由雷达设备在一次测量完成后发出，介于0和1之间，雷达设备提供每个测速单元的信号强度数据。通过信号强度与现场的流场资料的相关分析，输出每个单元格的可信指数。
步骤4，判定实时流速有效性：
（1）根据侧扫雷达的测速单元布设信息，在侧扫雷达工作的断面采用走航式ADCP（或流速仪等其他流速比测仪器）多次测量流速资料，率定出第i个流速单元处的水位-水面流速关系。此外，也可根据实际需要采用历史资料中的流速数据。
（2）根据率定出的第i个流速单元处的水位-水面流速关系，计算当前水位第i个流速单元的流速值偏离指数，作为判定测速单元是否有效的一个因子：
                    （4）
式中，
为当前水位第i个流速单元的流速值，m/s；
为侧扫雷达设备所能测到的最大流速，为侧扫雷达一项技术参数，m/s；
为侧扫雷达设备所能测到的最小流速，为侧扫雷达一项技术参数，m/s；
为由断面测流资料率定出的第i个流速单元处的水位-水面流速关系，若测量断面无率定成果，可暂时将直接取值即可。
步骤5，根据多个因子判断侧扫雷达单元流速有效性
其中每个单元的流速值有效性可通过下列公式来评判：
                                   （5）
式中，
为第i个流速单元的有效性指数，取值范围0-1，1表示有效性最好；
为第i个流速单元的位置指数，取值0或1，0表示不可用，1表示可用；
为第i个流速单元的信号可信度指数，取值范围0-1，1表示信号可信度最好；
为第i个流速单元的流速值偏离指数，非负数；
由式（5）算得的流速单元流速值有效性指数一般应在0.65以上可认为目标流速单元流速值是有效的，也可根据测流断面具体情况分析有效性指数的门槛值。
应用上述用于侧扫雷达的有效测速单元判定方法后得到各单元的流速有效性值，剔除流速值有效性指数低于给定数值的流速单元，基于有效的目标流速单元进行流速计算，能够得到更为精确的计算结果。结合前述步骤及本段内容，能够得到改进的流速计算方法。
需要说明的是，以上内容仅仅说明了本发明的技术思想，不能以此限定本发明的保护范围，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰均落入本发明权利要求书的保护范围之内。

Claims (9)

1.一种用于侧扫雷达的有效测速单元判定方法，其特征在于，包括如下步骤：
步骤1，获取流速单元数据
侧扫雷达设备通过测量计算，生成若干沿测流断面方向的动态数量的扇形水面流速单元；
步骤2，通过各流速单元的位置指数，判定动态流速单元有效性；
步骤3，通过水面信号质量，判定动态流速单元有效性；
步骤4，获取各流速单元的流速偏离指数，判定实时流速有效性；
步骤5，结合步骤2-4中得到的多个因子判断侧扫雷达单元流速有效性。
2.根据权利要求1所述的用于侧扫雷达的有效测速单元判定方法，其特征在于，所述步骤2包括如下子步骤：
（1）收集河道监测站的安装信息；
（2）根据测站安装信息，构建测站三维空间模型如下：
（1）
式中，为流速单元起点距，为雷达安装起点距，为雷达安装高程，为实时水位值，为雷达发射天线至表面流速单元位置的直线距离；
（3）为三维空间模型接入实时水位，计算第个流速单元的位置指数：
（2）
式中，为测流断面近雷达边界处起点距，为测流断面远雷达边界处起点距，由水位和测流断面的大断面资料可求得；为直线的水面投影，为流速单元近雷达边界与测流断面之间的垂直距离；为第个流速单元处水深，为雷达测速所需的最小水深。
3.根据权利要求2所述的用于侧扫雷达的有效测速单元判定方法，其特征在于，的计算方法如下：
。
4.根据权利要求2所述的用于侧扫雷达的有效测速单元判定方法，其特征在于，所述步骤2中的安装信息至少包括：雷达安装起点距、雷达安装高程、大断面信息。
5.根据权利要求1所述的用于侧扫雷达的有效测速单元判定方法，其特征在于，所述步骤3包括如下子步骤：
（1）根据侧扫雷达的测速单元布设信息，在侧扫雷达工作的断面采用历史资料或比测仪器多次测量流速信息；
（2）采集侧扫雷达相应单元格的信号质量信息；
（3）将同一时刻同一起点距单元格流速与侧扫雷达单元流速对比，分析出能准确反映流场状态的。
6.根据权利要求5所述的用于侧扫雷达的有效测速单元判定方法，其特征在于，所述步骤3中，通过侧扫雷达设备采集处理海量信号得出，其具体值由雷达设备在一次测量完成后发出。
7.根据权利要求1所述的用于侧扫雷达的有效测速单元判定方法，其特征在于，所述步骤4包括如下子步骤：
（1）根据侧扫雷达的测速单元布设信息，在侧扫雷达工作的断面采用历史资料或流速比测仪器多次测量流速，率定出第个流速单元处的水位-水面流速关系；
（2）根据率定出的第个流速单元处的水位-水面流速关系，计算当前水位第个流速单元的流速值偏离指数：
 （4）
式中，为当前水位第个流速单元的流速值，为侧扫雷达设备所能测到的最大流速，为侧扫雷达设备所能测到的最小流速，为由断面测流资料率定出的第个流速单元处的水位-水面流速关系。
8.根据权利要求6所述的用于侧扫雷达的有效测速单元判定方法，其特征在于，所述步骤4中，若测量断面无率定成果，将直接取值。
9.根据权利要求1所述的用于侧扫雷达的有效测速单元判定方法，其特征在于，所述步骤5包括如下过程：
通过下列公式计算每个单元的流速值有效性：
 （5）
式中，为第个流速单元的有效性指数，取值范围0-1；为第个流速单元的位置指数，取值0或1；为第个流速单元的信号可信度指数，取值范围0-1；为第个流速单元的流速值偏离指数。
CN202310053518.0A 2023-02-03 2023-02-03 一种用于侧扫雷达的有效测速单元判定方法 Active CN116125459B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310053518.0A CN116125459B (zh) 2023-02-03 2023-02-03 一种用于侧扫雷达的有效测速单元判定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310053518.0A CN116125459B (zh) 2023-02-03 2023-02-03 一种用于侧扫雷达的有效测速单元判定方法

Publications (2)

Publication Number Publication Date
CN116125459A true CN116125459A (zh) 2023-05-16
CN116125459B CN116125459B (zh) 2023-07-11

Family

ID=86309685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310053518.0A Active CN116125459B (zh) 2023-02-03 2023-02-03 一种用于侧扫雷达的有效测速单元判定方法

Country Status (1)

Country Link
CN (1) CN116125459B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103792533A (zh) * 2014-01-21 2014-05-14 北京艾力泰尔信息技术有限公司 基于固定点的河道断面多点测流方法
CN109270295A (zh) * 2018-08-20 2019-01-25 南京世海声学科技有限公司 一种基于自相关估计及有效数据筛选的水声多普勒流速测量方法
CN111474383A (zh) * 2020-04-23 2020-07-31 水利部南京水利水文自动化研究所 一种基于大数据的河流在线流量计算方法及系统
CN111487439A (zh) * 2019-01-25 2020-08-04 交通运输部天津水运工程科学研究所 一种基于无人船的声学多普勒流速剖面仪校准装置及方法
CN114414835A (zh) * 2022-01-26 2022-04-29 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） 一种固定式adcp测流抗干扰数据处理方法
EP4078095A1 (en) * 2019-12-16 2022-10-26 Flow-Tronic S.A. Non-invasive method and device to measure the flow rate of a river, open channel or fluid flowing in an underground pipe or channel

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103792533A (zh) * 2014-01-21 2014-05-14 北京艾力泰尔信息技术有限公司 基于固定点的河道断面多点测流方法
CN109270295A (zh) * 2018-08-20 2019-01-25 南京世海声学科技有限公司 一种基于自相关估计及有效数据筛选的水声多普勒流速测量方法
CN111487439A (zh) * 2019-01-25 2020-08-04 交通运输部天津水运工程科学研究所 一种基于无人船的声学多普勒流速剖面仪校准装置及方法
EP4078095A1 (en) * 2019-12-16 2022-10-26 Flow-Tronic S.A. Non-invasive method and device to measure the flow rate of a river, open channel or fluid flowing in an underground pipe or channel
CN111474383A (zh) * 2020-04-23 2020-07-31 水利部南京水利水文自动化研究所 一种基于大数据的河流在线流量计算方法及系统
CN114414835A (zh) * 2022-01-26 2022-04-29 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） 一种固定式adcp测流抗干扰数据处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
赵正军: "研究侧扫雷达测流系统功能与应用——以允景洪水文站为例", 《水利科学与寒区工程》, pages 139 - 142 *
韩宝豆: "基于雷达侧扫技术在高安水文站的应用", 《中国水运》, pages 102 - 104 *

Also Published As

Publication number Publication date
CN116125459B (zh) 2023-07-11

Similar Documents

Publication Publication Date Title
CN108254032B (zh) 河流超声波时差法流量计算方法
CN110146895B (zh) 基于倒置式多波束回声仪的声速剖面反演方法
Banner The influence of wave breaking on the surface pressure distribution in wind—wave interactions
US7530266B2 (en) Acoustic profiler for wind, temperature, and turbulence
CN113124941B (zh) 一种河道流量非接触式测量及精确计算方法
WO2020215458A1 (zh) 一种基于声学多普勒测流仪的流场智能计算方法及系统
CN106845018B (zh) 风电场对气象雷达降雨量影响的分析与定量化评估方法
CN110906992A (zh) 基于水平adcp施测垂线流速分布的河流流量测量方法
CN110672875B (zh) 基于Chirp-Z变换的表面水流速度检测方法
CN112711899A (zh) 一种蒸发波导高度的融合预测方法
WO2023050564A1 (zh) 风电机组的布局位置检测方法、模型训练方法及装置
CN110243423A (zh) 河道流量计算方法及系统
CN113155107A (zh) 一种不规则河道断面流量测量装置及方法
CN108663727B (zh) 利用蒸发率在世界海域范围内估算蒸发波导高度的方法
CN113626990A (zh) 基于风功率预测测风塔的风电机组功率曲线验证方法
CN106014878B (zh) 风力发电机组偏航系统动作误差的测试方法及系统
CN115876288A (zh) 一种基于大数据的电子仪表故障分析方法及系统
CN116125459B (zh) 一种用于侧扫雷达的有效测速单元判定方法
CN103852813A (zh) 雨滴三维尺度检测装置及利用该装置计算雨滴体积的方法
CN113256990A (zh) 基于聚类算法的雷达采集道路车辆信息的方法及系统
CN116222676B (zh) 精准定位的毫米波水流量监测方法及系统
CN102298072A (zh) 一种微差压式高精度测风装置及其测量方法
CN203745670U (zh) 雨滴三维尺度检测装置
CN113759387B (zh) 一种基于三维激光雷达的海岸防浪建筑物越浪量测量方法
CN111487616B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240229

Address after: 210 000 No. 95 Tiexinqiao Street, Yuhuatai District, Nanjing, Jiangsu Province

Patentee after: NANJING AUTOMATION INSTITUTE OF WATER CONSERVANCY AND HYDROLOGY, MINISTRY OF WATER RESOURCES

Country or region after: China

Patentee after: JIANGSU NAIWCH COOPERATION

Address before: 210 000 No. 95 Tiexinqiao Street, Yuhuatai District, Nanjing, Jiangsu Province

Patentee before: NANJING AUTOMATION INSTITUTE OF WATER CONSERVANCY AND HYDROLOGY, MINISTRY OF WATER RESOURCES

Country or region before: China