CN116577772A - 流速测量方法、流速流量仪、计算机设备及存储介质 - Google Patents

流速测量方法、流速流量仪、计算机设备及存储介质 Download PDF

Info

Publication number
CN116577772A
CN116577772A CN202310862467.6A CN202310862467A CN116577772A CN 116577772 A CN116577772 A CN 116577772A CN 202310862467 A CN202310862467 A CN 202310862467A CN 116577772 A CN116577772 A CN 116577772A
Authority
CN
China
Prior art keywords
flow rate
flow
path
phase difference
density
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
CN202310862467.6A
Other languages
English (en)
Other versions
CN116577772B (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.)
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 CN202310862467.6A priority Critical patent/CN116577772B/zh
Publication of CN116577772A publication Critical patent/CN116577772A/zh
Application granted granted Critical
Publication of CN116577772B publication Critical patent/CN116577772B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C13/00Surveying specially adapted to open water, e.g. sea, lake, river or canal
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/10Pre-processing; Data cleansing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Hydrology & Water Resources (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明涉及水文测量技术领域,具体涉及流速测量方法、流速流量仪、计算机设备及存储介质。包括如下步骤:S1:对一帧的ADC采样信号的I路和Q路分别进行DFT处理,分别得到I路和Q路的相位;S2:将I路的相位和Q路的相位作差,得到相位差;S3:根据目标值,对相位差进行初筛,再对初筛后的相位差进行干扰去除;S4:根据干扰去除后的相位差,计算得出流速;S5:将每一帧的流速组合,形成多帧的流速,对多帧的流速进行自动聚类;S6:流速按照密度准则分类,密度最大的一类求平均后作为此次测量的流速结果,其密度为置信度,将小密度流速滤除。从而能够提高对水流流速测量的准确性。

Description

流速测量方法、流速流量仪、计算机设备及存储介质
技术领域
本发明涉及水文测量技术领域,具体涉及流速测量方法、流速流量仪、计算机设备及存储介质。
背景技术
水文指的是自然界中水的变化、运动等的各种现象。建国60年来,水文事业随着国民经济和社会的发展而不断发展,取得了巨大成就。1949年,水位站203处、雨量站2处。我国已建成各类水文站点3万多处,形成包括水位、流量、雨量、水质、地下水、蒸发、泥沙等项目齐全、布局比较合理的水文站网。
其中,对水流流速的测量成为了研究水文的重要参照之一,而流速测量算法是实现流速测量的理论基石。而当今所采用的传统流速测量算法则使用了传统的运行框架。公开号为CN115526074A的中国发明专利,公开了基于河流表面流速的截面流速分布反演方法,该方法通过分析截面内区流速分布和外区流速分布的差异和联系获取河流截面流速分布,解决了微波雷达无法测量水面以下流速的问题。公开号为CN116222676A的中国发明专利,公开了一种精准定位的毫米波水流量监测方法及系统,通过毫米波射雷达及数据处理单元对河面某一网格内的流速进行无接触监测,由于摄像头指向与天线组的波束指向一致,因此摄像头获取采集点所在河面中的位置,并根据河道内垂线流速的分布曲线来分析计算河道整体的流速。
但在实际的工作环境中,在上述的测量中会存在一些误差因素,误差因素会与水流一起处理后求得的相对延迟量会产生较大偏差。这会出现对目标频带的识别有一定误差,进而对于幅度无法做到精确检测。因此,常规的测量方法得出的流速会存在一定的误差,降低流速测量的准确性。
上述问题是本领域亟需解决的问题。
发明内容
本申请人针对上述现有生产技术中的缺点,提供一种流速测量方法,从而能够提高对水流流速测量的准确性。
为了解决上述技术问题,本发明提供的方案是:一种流速测量方法,包括如下步骤:
S1:对一帧的ADC采样信号的I路和Q路分别进行DFT处理,分别得到I路和Q路的相位;
S2:将I路的相位和Q路的相位作差,得到相位差;
S3:根据90°目标值,对相位差进行初筛,再对初筛后的相位差进行干扰去除;
S4:根据干扰去除后的相位差,计算得出流速;
S5:将每一帧的流速组合,形成多帧的流速,对多帧的流速进行K-mean自动聚类;
S6:流速按照密度准则分类,密度最大的一类求平均后作为此次测量的流速结果,其密度为置信度,将小密度流速滤除。
作为本发明的进一步改进,所述S3中初筛的条件为,相位差在45°-135°之间的子载波。
作为本发明的进一步改进,所述S3中的所述干扰去除条件为:相位差在60°-120°之间,初筛目标频带内最大幅度大于目标值,初筛目标加权相位差后在75°-105°之间。
一种流速流量仪,流速流量仪包括所述流速测量方法并可以根据流速测量方法得出的流速数据解算出流量数据。
一种计算机设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述流速测量方法。
一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现所述流速测量方法。
本发明的有益效果:
本发明逻辑合理,通过设定目标值,目标值为90°,之后根据90°两侧的范围对筛选目标依次做初筛与干扰去除处理,以保障频带内所有被选出符合条件的数据的准确性,进而能够筛除放大器带来的误差因素,避免放大器带来的误差,使得本算法与传统的幅度检测相比更加准确。因而保障了流速计算的精确度。
附图说明
图1是本发明的流程示意图;
图2是子载波指数和幅度与相位差之间的关系图;
图3是捕获指数与流速之间的关系图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好地理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
参照图1所示,为本发明的一实施例,通过以下步骤实现对流速的测量,S1:对一帧的ADC采样信号的I路和Q路分别进行DFT处理,分别得到I路和Q路的相位;S2:将I路的相位和Q路的相位作差,得到相位差;S3:根据目标值,对相位差进行初筛,再对初筛后的相位差进行干扰去除;S4:根据干扰去除后的相位差,计算得出流速;S5:将每一帧的流速组合,形成多帧的流速,对多帧的流速进行自动聚类;S6:流速按照密度准则分类,密度最大的一类求平均后作为此次测量的流速结果,其密度为置信度,将小密度流速滤除。
在常规的水流流速的测量中,使用超声波测量水流流速时,会得到很多干扰的相位差,为了能够提高测量的准确性,需要将干扰相位差去除。在现实的测量中,ADC采样信号中会存在水流信号和实际影响信号,实际影响信号是客观存在无法消除的,因此水流信号和实际影响信号的混合信号默认为水流信号。之后水流信号的I路和Q路分别进行DFT处理,DFT处理后的I路和Q路还需要经过放大器对信号进行放大,在放大的过程中,放大器会对信号产生一定的误差因素,如何去除放大器产生的误差因素成为提高流速测量准确性的关键。因此,将这些误差因素排除,得到准确的相位差,才能够计算出更加准确的流速。由于I路和Q路是对信号的正交处理,因此I路和Q路的相位差的理论值永远保持在90°,而放大器产生的误差因素的相位差并不是90°,因此将目标值设定在90°。由于90°是理论值,在现实中很难达到90°,因此需要选取90°两侧的数值,选取靠近90°的相位差就能够确定该相位差为水流的相位差。
首先,对一帧的ADC采样信号的I路和Q路分别进行DFT处理,DFT处理为离散傅里叶变换,通过DFT处理能够分别得到I路和Q路的角度信息,对于任一复数z=a+bi,先对复数做tanz=b/a运算,再对所求结果做tan^(-1)运算以求出复数信号的相位,分别选取其中的相位信息。之后将I路的相位和Q路的相位作差,得到相位差。根据目标值90°,对相位差进行初筛,初筛的条件为,相位差在45°-135°之间的子载波,再对初筛后的相位差进行干扰去除。干扰去除包括如下条件,相位差在60°-120°之间,初筛目标频带内最大幅度大于目标值,初筛目标加权相位差后在75°-105°之间。通过上述的筛选和干扰去除,得到的相位差就为水流的相位差,之后通过相位差法求得该水流的流速。
上述步骤的计算方式是通过每抓取一帧进而算出所抓取帧的速度。由于每次采集均采用多帧采集,因此需要测量多个速度。由于环境的变化,会出现噪声干扰等影响。这会导致这些速度的被测结果的与理想值会存在差异,此时所测数据的方差很大,导致这些数据非常离散。因此在运算前需要找出离散数据中的奇异值并做筛除处理。
根据上述的步骤,求得每一帧的流速,之后将每一帧的流速组合,形成多帧的流速,接着对多帧的流速进行K-mean自动聚类,流速按照密度准则分类,密度最大的一类求平均后作为此次测量的流速结果,其密度为置信度,将小密度流速滤除。
如图2所示,是上述对一帧的I路和Q路处理后的子载波索引,左侧的纵坐标为重要性,右侧的纵坐标为相位差(度)。其中椭圆圈出的部分为干扰去除后的相位差。
如图3所示,是上述对多帧流速的K-mean自动聚类的示意图,其横坐标为帧数,纵坐标为速度(m/s),其中实线为聚类前的数据,虚线为聚类后的数据。该图为30帧流速的自动聚类效果图。流速按密度准则分类,每类流速的平均值即为其流速。因此需要找出密度最大的一类。舍弃密度较小的速度类。保留密度最大一类的帧数范围是7~30帧,对7~30帧内所有的流速值求平均后即作为本次测量的流速结果。最后求解置信度,利用公式C=T/N,其中C为置信度,T代表实际检测结果,N代表总数。舍弃流速条件小于等于0m/s的数据,根据对图表分析舍弃0-6帧的三组流速数据,保留在7-30帧的24个有效数据。所以根据图表计算出的置信度结果为80%。其中,小密度流速自动被滤除。置信度的定义为,在统计学中,一个概率样本的置信区间是对这个样本的某个总体参数的区间估计。置信区间展现的是这个参数的真实值有一定概率落在测量结果的周围的程度。置信区间给出的是被测量参数测量值的可信程度范围,即前面所要求的“一定概率”。这个概率被称为置信水平。
S5中的自动聚合采用K-mean自动聚类,K-mean算法是一种典型的基于划分的聚类算法,也是一种无监督学习算法。K-Means算法的思想很简单,对给定的样本集,用欧氏距离。作为衡量数据对象间相似度的指标,相似度与数据对象间的距离成反比,相似度越大,距离越小。
在一个实施例中,提供了一种流速流量仪,通过上述的流速测量方法实现对水流流速的测量并可以根据流速测量方法得出的流速数据解算出流量数据。
在一个实施例中,提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现以下步骤:
对一帧的ADC采样信号的I路和Q路分别进行DFT处理,分别得到I路和Q路的相位;
将I路的相位和Q路的相位作差,得到相位差;
根据目标值,对相位差进行初筛,再对初筛后的相位差进行干扰去除;
根据干扰去除后的相位差,计算得出流速;
将每一帧的流速组合,形成多帧的流速,对多帧的流速进行自动聚类;
流速按照密度准则分类,密度最大的一类求平均后作为此次测量的流速结果,其密度为置信度,将小密度流速滤除。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
对一帧的ADC采样信号的I路和Q路分别进行DFT处理,分别得到I路和Q路的相位;
将I路的相位和Q路的相位作差,得到相位差;
根据目标值,对相位差进行初筛,再对初筛后的相位差进行干扰去除;
根据干扰去除后的相位差,计算得出流速;
将每一帧的流速组合,形成多帧的流速,对多帧的流速进行自动聚类;
流速按照密度准则分类,密度最大的一类求平均后作为此次测量的流速结果,其密度为置信度,将小密度流速滤除。
以上对本发明的较佳实施例进行了描述;需要理解的是,本发明并不局限于上述特定实施方式,其中未尽详细描述的设备和结构应该理解为用本领域中的普通方式予以实施;任何熟悉本领域的技术人员,在不脱离本发明技术方案作出许多可能的变动和修饰,或修改为等同变化的等效实施例,这并不影响本发明的实质内容;因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰,均仍属于本发明技术方案保护的范围内。

Claims (6)

1.一种流速测量方法,其特征在于,包括如下步骤:
S1:对一帧的ADC采样信号的I路和Q路分别进行DFT处理,分别得到I路和Q路的相位;
S2:将I路的相位和Q路的相位作差,得到相位差;
S3:根据90°目标值,对相位差进行初筛,再对初筛后的相位差进行干扰去除;
S4:根据干扰去除后的相位差,计算得出流速;
S5:将每一帧的流速组合,形成多帧的流速,对多帧的流速采用K-mean自动聚类;
S6:流速按照密度准则分类,密度最大的一类求平均后作为此次测量的流速结果,其密度为置信度,将小密度流速滤除。
2.如权利要求1所述的一种流速测量方法,其特征在于,所述S3中初筛的条件为,相位差在45°-135°之间的子载波。
3.如权利要求1所述的一种流速测量方法,其特征在于,所述S3中的所述干扰去除条件为:相位差在60°-120°之间,初筛目标频带内最大幅度大于目标值,初筛目标加权相位差后在75°-105°之间。
4.一种流速流量仪,其特征在于,流速流量仪包括如权利要求1至3任一项所述流速测量方法并可以根据所述流速测量方法得出的流速数据解算出流量数据。
5.一种计算机设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至3任一项所述流速测量方法。
6.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至3任一项所述流速测量方法。
CN202310862467.6A 2023-07-14 2023-07-14 流速测量方法、流速流量仪、计算机设备及存储介质 Active CN116577772B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310862467.6A CN116577772B (zh) 2023-07-14 2023-07-14 流速测量方法、流速流量仪、计算机设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310862467.6A CN116577772B (zh) 2023-07-14 2023-07-14 流速测量方法、流速流量仪、计算机设备及存储介质

Publications (2)

Publication Number Publication Date
CN116577772A true CN116577772A (zh) 2023-08-11
CN116577772B CN116577772B (zh) 2023-09-12

Family

ID=87540067

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310862467.6A Active CN116577772B (zh) 2023-07-14 2023-07-14 流速测量方法、流速流量仪、计算机设备及存储介质

Country Status (1)

Country Link
CN (1) CN116577772B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6119070A (en) * 1995-09-25 2000-09-12 Schlumberger Industries, S.A. Method for acoustically measuring a fluid flow rate
CN104007286A (zh) * 2014-06-10 2014-08-27 上海航征测控系统有限公司 一种用于河道表面流速的自动监测系统
CN104865581A (zh) * 2015-03-19 2015-08-26 中国科学院空间科学与应用研究中心 一种基于gnss-r技术的河流流速测量方法与系统
CN112819817A (zh) * 2021-02-26 2021-05-18 澜途集思生态科技集团有限公司 基于图计算的河流流速估计方法
CN113012195A (zh) * 2021-03-04 2021-06-22 西安电子科技大学 基于光流计算的河流表面流速估计的方法和可读存储介质
CN114002666A (zh) * 2021-10-26 2022-02-01 南京航空航天大学 任意天线构型下星载ati-sar洋流流速提取方法及设备
CN114252111A (zh) * 2021-12-14 2022-03-29 北京中科蓝电科技有限公司 基于雷达信号的河道流量测量方法
CN116222676A (zh) * 2023-05-08 2023-06-06 成都赐华科技有限公司 精准定位的毫米波水流量监测方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6119070A (en) * 1995-09-25 2000-09-12 Schlumberger Industries, S.A. Method for acoustically measuring a fluid flow rate
CN104007286A (zh) * 2014-06-10 2014-08-27 上海航征测控系统有限公司 一种用于河道表面流速的自动监测系统
CN104865581A (zh) * 2015-03-19 2015-08-26 中国科学院空间科学与应用研究中心 一种基于gnss-r技术的河流流速测量方法与系统
CN112819817A (zh) * 2021-02-26 2021-05-18 澜途集思生态科技集团有限公司 基于图计算的河流流速估计方法
CN113012195A (zh) * 2021-03-04 2021-06-22 西安电子科技大学 基于光流计算的河流表面流速估计的方法和可读存储介质
CN114002666A (zh) * 2021-10-26 2022-02-01 南京航空航天大学 任意天线构型下星载ati-sar洋流流速提取方法及设备
CN114252111A (zh) * 2021-12-14 2022-03-29 北京中科蓝电科技有限公司 基于雷达信号的河道流量测量方法
CN116222676A (zh) * 2023-05-08 2023-06-06 成都赐华科技有限公司 精准定位的毫米波水流量监测方法及系统

Also Published As

Publication number Publication date
CN116577772B (zh) 2023-09-12

Similar Documents

Publication Publication Date Title
CN107861117B (zh) 一种适用于连续波周界监视雷达的多目标参数测量方法
CN102692627B (zh) 用于处理多普勒测量信号的仪器和方法
CN108469608B (zh) 一种运动平台雷达多普勒质心精确估计方法
CN110596787A (zh) 一种基于x波段全固态双偏振雨量雷达的降水估测方法
CN109490862B (zh) 一种基于相位差分统计谱的载频估计方法
CN110210067B (zh) 一种基于测量轨迹的阈值直线确定方法、装置
CN110837079B (zh) 一种基于雷达的目标检测方法及装置
RU2718144C1 (ru) Способ классификации, определения координат и параметров движения шумящего в море объекта в инфразвуковом диапазоне частот
RU2739000C1 (ru) Гидроакустический комплекс для обнаружения движущегося подводного источника звука, измерения пеленга на источник звука и горизонта источника звука в мелком море
CN107678025B (zh) 海浪波高计算方法和装置、存储介质及处理器
CN105204017B (zh) 一种基于正则优化的高分辨雷达角跟踪方法
CN109901147A (zh) 超高频雷达强杂区船舶目标检测的改进ca-cfar方法及系统
CN115494497A (zh) 一种相位编码脉冲雷达信号测向方法
CN116224274A (zh) 基于动态门限因子的低慢小目标检测方法、装置及设备
CN113189554B (zh) 雷达实测回波数据的处理方法、电子设备及存储介质
CN116577772B (zh) 流速测量方法、流速流量仪、计算机设备及存储介质
CN111352083A (zh) 一种高频地波雷达多接收通道增益自动校准方法及装置
Nobach A model-free noise removal for the interpolation method of correlation and spectral estimation from laser Doppler data
Williams et al. Estimation of hurricane winds from SeaWinds at ultrahigh resolution
CN112763796B (zh) 一种高精度测量lfm载频的方法
CN113887450B (zh) 基于调制谱特征的水中目标在线筛选方法
CN113552563B (zh) 垂测信息和高频地波雷达杂波信息的对应性分析方法
CN116359851A (zh) 一种基于融合网络的雷达有源干扰检测识别方法及装置
CN111983579B (zh) 一种参差脉冲重复时间消除雷达速度模糊的方法
CN114070679A (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