CN106643940B - 基于回波能量计算超声波流量计传播时间的方法 - Google Patents
基于回波能量计算超声波流量计传播时间的方法 Download PDFInfo
- Publication number
- CN106643940B CN106643940B CN201710090304.5A CN201710090304A CN106643940B CN 106643940 B CN106643940 B CN 106643940B CN 201710090304 A CN201710090304 A CN 201710090304A CN 106643940 B CN106643940 B CN 106643940B
- Authority
- CN
- China
- Prior art keywords
- echo
- energy
- ultrasonic
- time
- period
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/66—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种基于回波能量计算超声波流量计传播时间的方法,本发明通过超声波换能器接收到的回波信号的采样值的最大值先寻找到wh1,计算wh1以及其附近几个周期的回波的能量,然后根据回波的最大能量寻找到最高峰wh,最后通过线性插值寻找过零点,这样避免了避免了因寻找回波最高峰错误,而导致计算出的传播时间出现一个周期或几个周期的偏差,提高了传播时间的测量准确度。
Description
技术领域
本发明涉及超声流量测量技术领域,特别涉及一种基于回波能量计算超声波流量计传播时间的方法。
背景技术
时差法超声流量计工作原理如图1所示,它是根据介质流速与超声波在介质中顺、逆传播时产生的时间差存在一定的线性关系原理进行测量的,只要准确测量顺逆流时间差,再根据流速与其线性关系求出流速,进而可以求出瞬时流量以及累积流量。
图1中,S1,S2分别为2个超声波换能器,V为液体流速,D为管道直径,L为超声波的声程,θ为超声波进入液体的入射角。t1为换能器S1发射、S2接收时,超声波在管道中传播时间,即顺流时间;t2为换能器S2发射、S1接收时,超声波在管道中传播时间,即逆流时间。
超声流量计顺流传播时间t1和逆流传播时间t2分别用下式计算,即:
式中,C为超声波信号在水中的声速,设ΔT为顺逆流时间差,则:
由于超声波在流体中传播的速度C远远大于被测流体的实际流速V,即C2>>V2,所以(3)式可以简化为:
这里,超声波流量计测得的V是沿声道上的流体的线平均流速,而进行流量计算时需要的是截面平均流速,所以我们要将线平均速度V乘以流体修正系数k最后再乘以管道的截面积S便可求出瞬时流量Q。其公式如下:
Q=S·kV=(πD2/4)kV。
目前,一般采用对超声波换能器接收到的回波信号进行采样,获取整个回波波形的采样值,如图2所示,回波信号中,Y轴为信号幅值,X轴为信号坐标。根据回波波形采样值中的最大值xmax确定回波波形中的最高峰wh。由于现有算法通过超声波换能器接收到的回波信号的采样值的最大值xmax寻找到最高峰wh,而往往超声波换能器接收到的回波的最高峰wh相邻几个周期内的回波的幅值最大值都很接近,如图4所示。图4中回波wh1和w4的幅值最大值很接近,在超声波流量管道中有流量的情况下,很容易判断错误最高峰,从而导致根据最高峰计算出的过零点可能存在一个周期或几个周期的误差,使计算出的传播时间存在一个周期或几个周期的误差。
发明内容
本发明的目的是针对现有技术存在的不足,提供一种基于回波能量计算超声波流量计传播时间的方法,其提高超声波流量计传播时间测量的准确性,减少测量误差。
本发明的目的是采用下述方案实现的:一种基于回波能量计算超声波流量计传播时间的方法,超声波换能器发射端发射超声波信号,超声波换能器接收端接收超声回波信号,超声波传播时间的计算方法包括如下步骤:
1)找到超声波换能器接收端接收到的回波信号采样值内的回波幅值最大点;
2)分别求取该回波幅值最大点所处周期内的回波的能量以及与该回波幅值最大点所处周期相近的几个周期内的回波的能量,确定回波能量最大值Emax;
3)将步骤2)确定的回波能量最大值Emax所对应的周期内的回波作为最高峰wh,计算得出该最高峰wh的过零点的时间t0,根据公式T=t0+t'即可求出超声波传播时间T,式中,t'为一个与超声波流量计相关的固定值。
各个周期内的回波的能量为该周期内的所有采样点的幅值相加的总和。
在最高峰wh中,如果相邻两采样点的幅值极性相反时则进行线性插值运算,计算出最高峰wh的信号过零点时间t0。
在零流量的情况下计算t':当超声波流量计管道中的流体静止时通过测量声道的长度以及根据此时的声速计算出此时准确的传播时间T0,同时计算出零流速下的过零点对应的时间t'0,则t'=t'0-T0。
步骤2)分别求取该回波幅值最大点所处周期内的回波的能量以及与该回波幅值最大点所处周期相近的前后各两个或三个相邻周期内的回波的能量,确定回波能量最大值Emax。
本发明具有的优点是:本发明通过超声波换能器接收到的回波信号的采样值的最大值先寻找到wh1,计算wh1以及其附近几个周期的回波的能量,然后根据回波的最大能量寻找到最高峰wh,最后通过线性插值寻找过零点,这样避免了避免了因寻找回波最高峰错误,而导致计算出的传播时间出现一个周期或几个周期的偏差,提高了传播时间的测量准确度。
附图说明
图1为时差法超流量计的原理示意图;
图2为回波信号的示意图;
图3为本发明的方法流程图;
图4为回波信号与其回波能量的示意图;
图5为线性插值计算过零点t0的原理示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,给出了详细的实施方式和具体的操作过程,下面将结合附图对本发明作进一步的详细描述:
参见3至图5,一种基于回波能量计算超声波流量计传播时间的方法,超声波换能器发射端发射超声波信号,超声波换能器接收端接收超声回波信号,超声波传播时间的计算方法包括如下步骤:
1)找到超声波换能器接收端接收到的回波信号采样值内的回波幅值最大点xmax;
2)分别求取该回波幅值最大点xmax所处周期内的回波wh1的能量以及与该回波幅值最大点所处周期相近的几个周期内的回波的能量,确定上述周期内的回波中的回波能量最大值Emax。本发明一般分别求取包括回波wh1在内的多个相邻周期内的回波的能量。优选地,该步骤分别求取该回波幅值最大点所处周期内的回波的能量以及与该回波幅值最大点所处周期相近的前后各两个或三个相邻周期内的回波的能量,确定回波能量最大值Emax。如本实施例分别求取包括回波wh1在内的五个相邻周期内的回波w1、w2、wh1、w4、w5的能量E1、E2、Eh1、E4、E5。各个周期内的回波的能量为该周期内的所有采样点的幅值相加的总和。如回波能量Ei由下式计算:Ei=xi+xi+1+…+xi+n,式中,i为回波wi的坐标起始点,n为回波一个周期采样点个数,xi为回波幅值。回波w1、w2位于回波wh1前,回波w4、w5位于回波wh1后。本实施例中回波w4的能量E4最大。
3)将步骤2)确定的回波能量最大值Emax所对应的周期内的回波作为最高峰wh。本实施例的回波w4为最高峰wh。然后计算得出该最高峰wh的第一个过零点的时间t0,根据公式T=t0+t'即可求出超声波传播时间T,式中,t'为一个固定的值。在最高峰wh中,如果第n点的幅值xn与第n+1点的幅值xn+1极性相反即相邻两采样点的幅值极性相反时则进行线性插值运算,计算出最高峰wh的信号过零点时间t0。将计算出的回波信号坐标x过零点作为信号过零点时间t0,如图5所示,线性插值计算的过零点时间t0:
本实施例t'为一个与超声波流量计相关的固定值。在零流量的情况下计算t':当超声波流量计管道中的流体静止时通过测量声道的长度以及根据此时的声速计算出此时准确的传播时间T0,同时计算出零流速下的过零点对应的时间t'0,则t'=t'0-T0。
本发明通过上述方法求取的顺流时间t1和逆流时间t2具有较高的精确度,能够满足超声波流量计的测量要求。
以上所述仅为本发明的优选实施例,并不用于限制本发明,显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (5)
1.一种基于回波能量计算超声波流量计传播时间的方法,其特征在于,超声波换能器发射端发射超声波信号,超声波换能器接收端接收超声回波信号,超声波传播时间的计算方法包括如下步骤:
1)找到超声波换能器接收端接收到的回波信号采样值内的回波幅值最大点;
2)分别求取该回波幅值最大点所处周期内的回波的能量以及与该回波幅值最大点所处周期相近的几个周期内的回波的能量,确定回波能量最大值Emax;
3)将步骤2)确定的回波能量最大值Emax所对应的周期内的回波作为最高峰wh,计算得出该最高峰wh的过零点的时间t0,根据公式T=t0+t'即可求出超声波传播时间T,式中,t'为一个与超声波流量计相关的固定值。
2.根据权利要求1所述的方法,其特征在于:各个周期内的回波的能量为该周期内的所有采样点的幅值相加的总和。
3.根据权利要求1所述的方法,其特征在于:在最高峰wh中,如果相邻两采样点的幅值极性相反时则进行线性插值运算,计算出最高峰wh的信号过零点时间t0。
4.根据权利要求1或3所述的方法,其特征在于:在零流量的情况下计算t':当超声波流量计管道中的流体静止时通过测量声道的长度以及根据此时的声速计算出此时准确的传播时间T0,同时计算出零流速下的过零点对应的时间t'0,则t'=t'0-T0。
5.根据权利要求1所述的方法,其特征在于:步骤2)分别求取该回波幅值最大点所处周期内的回波的能量以及与该回波幅值最大点所处周期相近的前后各两个或三个相邻周期内的回波的能量,确定回波能量最大值Emax。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710090304.5A CN106643940B (zh) | 2017-02-20 | 2017-02-20 | 基于回波能量计算超声波流量计传播时间的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710090304.5A CN106643940B (zh) | 2017-02-20 | 2017-02-20 | 基于回波能量计算超声波流量计传播时间的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106643940A CN106643940A (zh) | 2017-05-10 |
CN106643940B true CN106643940B (zh) | 2019-05-17 |
Family
ID=58845177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710090304.5A Active CN106643940B (zh) | 2017-02-20 | 2017-02-20 | 基于回波能量计算超声波流量计传播时间的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106643940B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107655533B (zh) * | 2017-10-31 | 2019-09-10 | 合肥工业大学 | 一种基于回波能量积分的气体超声波流量计信号处理方法和系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4080574A (en) * | 1974-01-31 | 1978-03-21 | United Kingdom Atomic Energy Authority | Apparatus for providing time reference signals |
CN102297712A (zh) * | 2011-07-12 | 2011-12-28 | 北京理工大学 | 一种超声回波传播时间测量方法 |
CN103344288A (zh) * | 2013-06-25 | 2013-10-09 | 国家电网公司 | 一种基于零点分析的时差式超声波流量计测量方法 |
CN104330120A (zh) * | 2014-10-28 | 2015-02-04 | 姜跃炜 | 用于低能耗超声波流量表的流量检测方法及系统 |
CN104697593A (zh) * | 2015-03-24 | 2015-06-10 | 合肥工业大学 | 一种基于fpga和dsp的气体超声流量计 |
CN105628115A (zh) * | 2015-12-30 | 2016-06-01 | 浙江大学 | 一种应用于时差式超声波流量计的脉冲噪声滤除方法 |
-
2017
- 2017-02-20 CN CN201710090304.5A patent/CN106643940B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4080574A (en) * | 1974-01-31 | 1978-03-21 | United Kingdom Atomic Energy Authority | Apparatus for providing time reference signals |
CN102297712A (zh) * | 2011-07-12 | 2011-12-28 | 北京理工大学 | 一种超声回波传播时间测量方法 |
CN103344288A (zh) * | 2013-06-25 | 2013-10-09 | 国家电网公司 | 一种基于零点分析的时差式超声波流量计测量方法 |
CN104330120A (zh) * | 2014-10-28 | 2015-02-04 | 姜跃炜 | 用于低能耗超声波流量表的流量检测方法及系统 |
CN104697593A (zh) * | 2015-03-24 | 2015-06-10 | 合肥工业大学 | 一种基于fpga和dsp的气体超声流量计 |
CN105628115A (zh) * | 2015-12-30 | 2016-06-01 | 浙江大学 | 一种应用于时差式超声波流量计的脉冲噪声滤除方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106643940A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106643939B (zh) | 用于超声波流量计计算超声波传播时间的方法 | |
CN106768109B (zh) | 基于回波上升段峰值拟合和基于回波能量点定位的气体超声波流量计信号处理方法 | |
CN103344288B (zh) | 一种基于零点分析的时差式超声波流量计测量方法 | |
CN104501889B (zh) | 基于互相关时差法超声波流量的检测方法 | |
Zhou et al. | Multipath ultrasonic gas flow-meter based on multiple reference waves | |
CN109029598B (zh) | 一种互相关法气体超声波流量计中基于主次包络结合处的参考波形选取方法 | |
JP2010512510A (ja) | 超音波式流量計測方法及びシステム | |
CN107860430B (zh) | 一种基于时差法的超声波气体流量计时间差测量方法 | |
JP4976287B2 (ja) | パルス波形検出による超音波信号の受信点検出 | |
JPWO2005083372A1 (ja) | パルスドップラ方式と伝搬時間差方式の両方式対応型超音波流量計、同流量計において測定方式を自動選択する方法およびプログラム、同流量計用の電子装置 | |
CN106706056B (zh) | 一种大口径超声波水表流量测量的补偿方法 | |
CN106643937A (zh) | 一种基于超声波流量计的流量测量方法及装置 | |
CN108051036A (zh) | 非满管的超声波流量计及超声波流量测量系统 | |
CN106643940B (zh) | 基于回波能量计算超声波流量计传播时间的方法 | |
Han et al. | Studies on the transducers of clamp-on transit-time ultrasonic flow meter | |
US10955273B2 (en) | Extended range ADC flow meter | |
JP2004317459A (ja) | 超音波式ガス濃度流量測定方法および装置 | |
CN112147366A (zh) | 用于使用超声波来测量流体速度的方法 | |
KR20170035555A (ko) | 초음파 유량계의 신호처리시스템 | |
US20230243682A1 (en) | Ultrasonic flow measurement | |
JP2013250254A (ja) | 超音波式スパイロメータの多重反射防止整流管 | |
WO2019047634A1 (zh) | 一种直线距离上平均流体流速测量系统 | |
Li et al. | Research on transit-time ultrasonic flowmeter with signal characteristic analysis | |
Waluś | Mathematical modelling of an ultrasonic flowmeter primary device | |
US7617738B2 (en) | Method and apparatus for measuring flow rate of fluid |
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 |