CN109029602B - 基于超声波的流量测量方法及流量计 - Google Patents

基于超声波的流量测量方法及流量计 Download PDF

Info

Publication number
CN109029602B
CN109029602B CN201810987006.0A CN201810987006A CN109029602B CN 109029602 B CN109029602 B CN 109029602B CN 201810987006 A CN201810987006 A CN 201810987006A CN 109029602 B CN109029602 B CN 109029602B
Authority
CN
China
Prior art keywords
ultrasonic
frequency
signal
max
flow
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
Application number
CN201810987006.0A
Other languages
English (en)
Other versions
CN109029602A (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.)
Taihua Wisdom Industry Group Co Ltd
Original Assignee
Taihua Wisdom Industry Group 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 Taihua Wisdom Industry Group Co Ltd filed Critical Taihua Wisdom Industry Group Co Ltd
Priority to CN201810987006.0A priority Critical patent/CN109029602B/zh
Publication of CN109029602A publication Critical patent/CN109029602A/zh
Application granted granted Critical
Publication of CN109029602B publication Critical patent/CN109029602B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/66Measuring 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
    • G01F1/663Measuring 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 by measuring Doppler frequency shift

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明公开了一种基于超声波的流量测量方法及流量计,所述测量方法包括:步骤S1:设置超声波出射信号的频率f0;步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算FFT采样频率和混频基准信号fr;步骤S3:对超声波回波信号进行模数采样处理,得到采样信号;步骤S4:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱;步骤S5:根据所述功率谱计算得到流速v;步骤S6:检测v是否大于Q*vmax,其中Q小于1;步骤S7:如果v>Q*vmax,则根据流速v计算流量;步骤S8:如果v≤Q*vmax,则设置vmax=Q*vmax,返回步骤S2。通过发明提高了低流速下的流量测量精度。

Description

基于超声波的流量测量方法及流量计
技术领域
本发明涉及流量测量技术领域,更具体地,涉及一种基于超声波的流量测量方法及流量计。
背景技术
多普勒超声波流量计应用声学多普勒效应,根据反射声波和发射声波之间的频率差,即多普勒频移来计算被测液体的流速和流量。多普勒超声波流量计具有流速响应快,对流体压力、密度等因素不敏感,无零点偏移等优点,但是在低流速情况下,多普勒超声波流量计测量精度较低,测量误差较大。
在运用FFT算法(快速傅里叶变换)进行多普勒信号频谱分析的流量计中,往往通过提高FFT频谱分辨率的方法提高测量精度,但是在满足多普勒信号采样频率的前提下,只能通过提高采样点数来提高FFT频谱分辨率,这样不仅增加了软件处理的工作,也要消耗更多的硬件存储空间和计算单元,从而增加了软硬件开销,不利于成本控制。
发明内容
有鉴于此,本发明提供了一种基于超声波的流量测量方法及流量计,以提高低流速下的流量测量精度。
一种基于超声波的流量测量方法,包括:
步骤S1:设置超声波出射信号的频率f0
步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算FFT采样频率和混频基准信号fr
其中,采用以下公式计算所述FFT采样频率fs
fs=4fd
Figure BDA0001779925480000011
c为超声波在水中传播速度,θ为声波入射角;
采用以下公式计算所述混频基准信号fr
fr=f0-fc,fc=fd
步骤S3:对超声波回波信号进行模数采样处理,得到采样信号;
步骤S4:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱;
步骤S5:根据所述功率谱计算得到流速v;
步骤S6:检测v是否大于Q*vmax,其中Q小于1;
步骤S7:如果v>Q*vmax,则根据流速v计算流量;
步骤S8:如果v≤Q*vmax,则设置vmax=Q*vmax,返回步骤S2。
可选地,步骤S2还包括:设置测量次数x为初始值,其中,所述初始值为0;
步骤S7包括:
步骤S71:记录流速v至流速组,设置x=x+1;
步骤S72:判断x是否达到预定测量次数,其中,当x达到预定测量次数时,执行步骤S73,当未达到所述预定测量次数时,返回步骤S3;
步骤S73:计算所述流速组中所有v的平均值,并根据所述平均值计算所述流量。
可选地,在步骤S3中,以1/2倍的所述FFT采样频率fs对超声波回波信号进行模数采样处理。
可选地,在步骤S2之前,所述方法还包括:
步骤S9:向超声发射换能器发送PWM信号,所述PWM信号的频率为所述出射信号的频率。
可选地,在步骤S7中,采用以下公式计算所述流量:
所述流量=流水截面积*v,其中流水截面积为已知量。
可选地,在步骤S5中,根据所述功率谱计算计算所述流速v包括:
选取功率谱上幅值最大的点的频率作为计算多普勒频移信号f′1,则多普勒回波信号的频率f1为f1=fr+f′1
采用以下公式计算所述流速v,
Figure BDA0001779925480000021
一种基于超声波的流量计,所述流量计包括:
中心处理模块以及分别与所述中心处理模块连接的超声接收换能器和模数转换芯片;
所述超声接收换能器用于接收超声波回波信号;
所述模数转换芯片用于对所述超声波回波信号进行模数采样;
所述中心处理模块用于执行以下步骤计算得到流速v:
步骤S1:设置超声波出射信号的频率f0
步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算混频基准信号fr以及FFT采样频率,其中,采用以下公式计算所述FFT采样频率fs
fs=4fd
Figure BDA0001779925480000031
c为超声波在水中传播速度,θ为声波入射角;
采用以下公式计算所述混频基准信号fr:
fr=f0-fc,fc=fd
步骤S3:控制所述模数转换芯片对所述超声波回波信号进行模数采样处理,得到采样信号;
步骤S4:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱;
步骤S5:根据所述功率谱计算得到流速v;
步骤S6:检测v是否大于Q*vmax,其中Q小于1;
步骤S7:如果v>Q*vmax,则根据v计算流量;
步骤S8:如果v≤Q*vmax,则设置vmax=Q*vmax,返回步骤S2。
可选地,所述中心处理模块执行的步骤S2还包括:设置测量次数x为初始值,其中,所述初始值为0;
执行的步骤S7还包括:
步骤S71:记录v至流速组,设置x=x+1;
步骤S72:判断x是否达到预定测量次数,其中,当x达到预定测量次数时,执行步骤S73,当未达到所述预定测量次数时,返回步骤S3;
步骤S73:计算所述流速组中所有v的平均值,并根据所述平均值计算所述流量。
可选地,还包括:与所述中心处理模块连接的超声发射换能器;
所述中心处理模块还用于在执行步骤S2之前,向所述超声发射换能器发送PWM信号,所述PWM信号的频率为所述出射信号的频率;
所述超声发射换能器用于根据所述PWM信号生成并发射所述超声波出射信号。
可选地,还包括:信号处理电路,所述信号处理电路连接在所述超声接收换能器与所述模数转换芯片之间,用于在所述模数转换芯片采样之前对所述超声波回波信号进行滤波和放大。
与现有技术相比,本发明提供的基于超声波的流量测量方法及流量计,至少实现了如下的有益效果:
1)在不增加硬件开销的前提下,提高了流量计在低流速情况下的测量精度。
2)通过对流量平均值的计算进一步提高了流量计算精度以及频谱分辨率。
当然,实施本发明的任一产品必不特定需要同时达到以上所述的所有技术效果。
通过以下参照附图对本发明的示例性实施例的详细描述,本发明的其它特征及其优点将会变得清楚。
附图说明
被结合在说明书中并构成说明书的一部分的附图示出了本发明的实施例,并且连同其说明一起用于解释本发明的原理。
图1是本发明实施例基于超声波的流量测量方法的一种流程图;
图2是本发明实施例中换能器发射的声波入射角示意图。
图3是本发明实施例基于超声波的流量计的一种结构示意图;
图4是本发明实施例基于超声波的流量计的另一种结构示意图。
具体实施方式
现在将参照附图来详细描述本发明的各种示例性实施例。应注意到:除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。
以下对至少一个示例性实施例的描述实际上仅仅是说明性的,决不作为对本发明及其应用或使用的任何限制。
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论,但在适当情况下,所述技术、方法和设备应当被视为说明书的一部分。
在这里示出和讨论的所有例子中,任何具体值应被解释为仅仅是示例性的,而不是作为限制。因此,示例性实施例的其它例子可以具有不同的值。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步讨论。
实施例1
如图1所示是本发明实施例基于超声波的流量测量方法的一种流程图,具体包括以下步骤:
步骤100:开始。
步骤101:设置超声波出射信号的频率。
步骤102:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算FFT采样频率和混频基准信号fr,其中,采用以下公式计算所述FFT采样频率fs:fs=4fdc为超声波在水中传播速度,θ为声波入射角;fd为多普勒频移;fmax为最大流速下的回波频率。具体地,如图2所示,声波入射角θ为超声发射换能器H发射的超声波与被测液体流速v之间的夹角,其中,超声发射换能器可以由中心处理模块控制,也可以由中心处理模块之外的其他控制器控制。
采用以下公式计算所述混频基准信号fr:fr=f0-fc,fc=fd,fc为混频信号的中频基准频率。
步骤103:对超声波回波信号进行模数采样处理,得到采样信号。
步骤104:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱。
步骤105:根据所述功率谱计算得到流速v。
具体地,根据所述功率谱计算计算所述流速v包括:
选取功率谱上幅值最大的点的频率作为计算多普勒频移信号f′1,则多普勒回波信号的频率f1为f1=fr+f′1
采用以下公式计算所述流速v,
Figure BDA0001779925480000061
需要说明的是:则在超声波多普勒流量计的测量范围0~vmax内,包含流速信息的多普勒频移范围0~fd,即为ADC需要采样的对象;混频后的中频信号范围为fd~2fd,即为FFT采样及算法分析的对象。
根据采样定理,FFT采样频率至少为被采样信号频率的两倍,则令FFT的采样频率fs=4fd,采样点数为N,则FFT算法进行频谱分析的频谱分辨率为
Figure BDA0001779925480000062
假设根据FFT算法功率谱得到的幅值最大点为f′1,则多普勒回波信号为f1+fr+f′1根据公式可以得到本次测量的流速v。
进一步,当超声波发射频率f0不变时,频谱分辨率fh与流速分辨率vh呈正比关系。
步骤106:检测v是否大于Q*vmax,其中Q小于1;如果是,执行步骤107;否则,执行步骤109。
需要说明的是,Q理论上凡小于1的数均可以,均可以提高分辨率。在应用中,当Q为1/2可操作性更强,应用也最多。
步骤107:根据流速v计算流量。
具体地,采用以下公式计算所述流量:
所述流量=流水截面积*v,其中流水截面积为已知量。
步骤108:结束。
步骤109:设置vmax=Q*vmax,返回执行步骤102。
本发明实施例提供的基于超声波的流量测量方法,中心处理模块设置超声波出射信号的频率f0,根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算FFT采样频率和混频基准信号fr,对超声波回波信号进行模数采样处理,得到采样信号;将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱,根据所述功率谱计算得到流速v;如果流速v小于或等于Q*vmax,其中Q小于1,则重新对最大流速vmax进行设置,从而保证了本申请流量计算的精度。
实施例2
步骤200:开始。
步骤201:设置超声波出射信号的频率。
步骤202:设置测量次数x为初始值,其中,所述初始值为0;根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算FFT采样频率和混频基准信号fr,其中,采用以下公式计算所述FFT采样频率fs:fs=4fd
Figure BDA0001779925480000071
c为超声波在水中传播速度,θ为声波入射角;
采用以下公式计算所述混频基准信号fr:fr=f0-fc,fc=fd
步骤203:以1/2倍的所述FFT采样频率fs对超声波回波信号进行模数采样处理,得到采样信号。
步骤204:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱。
步骤205:根据所述功率谱计算得到流速v。
步骤206:检测流速v是否大于Q*vmax,其中Q小于1;如果是,执行步骤207;否则,执行步骤212。
步骤207:记录流速v至流速组,设置x=x+1。
步骤208:判断x是否达到预定测量次数;如果是,执行步骤209;否则,执行步骤203。
具体地,预定测量次数可以根据流量测量精度通过标定确定,比如,预定测量次数为:10~100之间的值。
步骤209:计算所述流速组中所有流速v的平均值。
具体地,采用以下公式计算所述平均值
Figure BDA0001779925480000072
Figure BDA0001779925480000081
具体地,x为测量次数,测量次数可以达到预定测量次数,预定测量次数可以根据流量测量精度通过标定确定,比如,预定测量次数为:10~100之间的值。
步骤210:根据所述平均值计算流量。
具体地,采用以下公式计算所述流量:
所述流量=流水截面积*
Figure BDA0001779925480000083
其中流水截面积为已知量。
步骤211:结束。
步骤212:设置vmax=Q*vmax,返回执行步骤202。
本发明实施例提供的基于超声波的流量测量方法,当流速v小于Q*vmax时,将每次计算得到的流速记录到流速组中,并在记录次数达到预设定测量次数后,计算流速组中所有流速v的平均值,并根据平均值计算流量,从而进一步保证了流量的测试精度。
实施例3
步骤300:开始。
步骤301:设置超声波出射信号的频率。
步骤302:向超声发射换能器发送PWM信号,所述PWM信号的频率为所述出射信号的频率。
步骤303:设置测量次数x为初始值,其中,所述初始值为0;根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算FFT采样频率和混频基准信号fr,其中,采用以下公式计算所述FFT采样频率fs:fs=4fd
Figure BDA0001779925480000082
c为超声波在水中传播速度,θ为声波入射角;
采用以下公式计算所述混频基准信号fr:fr=f0-fc,fc=fs
步骤304:以1/2倍所述FFT采样频率fs对超声波回波信号进行模数采样处理,得到采样信号。
步骤305:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱。
步骤306:根据所述功率谱计算得到流速v。
步骤307:检测流速v是否大于Q*vmax,Q小于1;如果是,执行步骤308;否则,执行步骤313。
步骤308:记录流速v至流速组,设置x=x+1。
步骤309:判断x是否达到预定测量次数;如果是,执行步骤310;否则,执行步骤304。
具体地,预定测量次数可以根据流量测量精度通过标定确定,比如,预定测量次数为:10~100之间的值。
步骤310:计算所述流速组中所有流速v的平均值。
步骤311:根据所述平均值计算流量。
步骤312:结束。
步骤313:设置vmax=Q*vmax,返回执行步骤303。
本发明实施例提供的基于超声波的流量测量方法,中心处理模块可以只接收超声波回波信号,可以发射完出射信号后,再接收此出射信号的回波信号,即中心处理模块向超声发射换能器发送PWM信号,所述PWM信号的频率为所述出射信号的频率。通过中心处理模块控制信号的发送与接收进一步保证了流量计算的精度。
实施例4
针对上述方法实施例,本发明还提供了一种基于超声波的流量计,如图3所示,所述流量计包括:中心处理模块以及分别与所述中心处理模块连接的超声接收换能器和模数转换芯片;所述超声接收换能器用于接收超声波回波信号;所述模数转换芯片用于对所述超声波回波信号进行模数采样;所述中心处理模块用于执行以下步骤计算得到流速v:
步骤S1:设置超声波出射信号的频率f0
步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算混频基准信号fr以及FFT采样频率,
其中,采用以下公式计算所述FFT采样频率fs:
fs=4fd
Figure BDA0001779925480000091
c为超声波在水中传播速度,θ为声波入射角;
采用以下公式计算所述混频基准信号fr
fr=f0-fc,fc=fd
步骤S3:控制所述模数转换芯片对所述超声波回波信号进行模数采样处理,得到采样信号。
步骤S4:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱。
步骤S5:根据所述功率谱计算得到流速v。
步骤S6:检测v是否大于Q*vmax,其中Q小于1。
步骤S7:如果v>Q*vmax,则根据v计算流量。
步骤S8:如果v≤Q*vmax,则设置vmax=Q*vmax,返回步骤S2。
进一步,所述中心处理模块执行的步骤S2还包括:设置测量次数x为初始值,其中,所述初始值为0;
执行的步骤S7还包括:
步骤S71:记录v至流速组,设置x=x+1。
步骤S72:判断x是否达到预定测量次数,其中,当x达到预定测量次数时,执行步骤S73,当未达到所述预定测量次数时,返回步骤S3。
步骤S73:计算所述流速组中所有v的平均值,并根据所述平均值计算所述流量。
进一步,本发明的一个实施例中,所述中心处理模块可以以1/2倍的所述FFT采样频率fs对超声波回波信号进行模数采样处理。
更进一步,本发明的另一实施例中,如图4所示,所述流量计还可以包括:
与所述中心处理模块连接的超声发射换能器;
所述中心处理模块还用于在执行步骤S2之前,向所述超声发射换能器发送PWM信号,所述PWM信号的频率为所述出射信号的频率;
所述超声发射换能器用于根据所述PWM信号生成并发射所述超声波出射信号。
进一步,本发明的一个实施例中,所述流量计还可以包括:信号处理电路,所述信号处理电路连接在所述超声接收换能器与所述模数转换芯片之间,用于在所述模数转换芯片采样之前对所述超声波回波信号进行滤波和放大。
进一步,本发明的另一个实施例中,所述流量计还可以包括:信号放大电路,所述信号放大电路连接在所述中心处理模块与所述超声发射换能器之间,用于对所述中心处理模块发送的PWM信号进行放大。
结合上述实施例,下面通过具体实施例对本发明基于超声波的流量测量方法及流量计进行详细介绍:
(1)首先通过流量计的中心控制模块产生频率为f0=1M的PWM信号,PWM信号经过信号放大电路放大,驱动超声发射换能器向水中发射超声波,由超声接收华能器对反射声波信号进行接收,并通过信号处理电路对信号进行放大滤波处理。
(2)将处理过的回波信号送入模数转换芯片进行采样,若多普勒超声波流量计测量的最大流速为vmax=6m/s,超声波的入射角θ=10°,则根据公式
Figure BDA0001779925480000111
可以得到回波信号的最大频率为:
Figure BDA0001779925480000112
其中c=1500m/s,为常温下声波在水中的传播速度。
(3)由上可得多普勒频移为fd=fmax-f0=1.008×106-1×106=8k,令混频信号的中频基准频率fc=fd=8k,则混频基准信号为fr=f0-fd=1×106-8×103=992×103,则经过混频后,得到的中频信号f′=fc+fd=2fd=16k,即为最大流速vmax时混频后得到的多普勒频移信号。同样的,可以计算得到当流速为0时,混频后的多普勒频移信号为f′0=fd=8k,即包含流速信息的多普勒信号频移范围为0~8k,为模数转换芯片采样对象,经过混频的中频频率范围为8k~16k,为FFT算法频谱分析的对象。
(4)根据采样定理,令fs=4fd=32k,ADC采样点数N=4096,则FFT算法频谱分析的频谱分辨率为
Figure BDA0001779925480000113
根据公式
Figure BDA0001779925480000114
可以得到多普勒超声波流量计的流速分辨率vh
(5)假设经过FFT算法的频谱分析,得到功率谱上功率最大点的频率为f′1=15k,则多普勒回波信号频率为f1=fr+f′1=992×103+15×103=1.007×106,可以得到流速v1
Figure BDA0001779925480000122
由于
Figure BDA0001779925480000123
则重复以上步骤,采样计算x次,得到平均流速值
Figure BDA0001779925480000124
根据公式流量=流水截面积*
Figure BDA0001779925480000125
即可以得到测量的流量值Q。
(6)如果经过FFT算法频谱分析,得到功率谱上的功率最大点f′1=9k,则根据第(5)步的方法,可以得到回波信号为f1=1.001×106,流速v1≈0.76m/s,由于v1<3m/s,为了提高测量精度,我们进入第二阶测量。
(7)根据第(2)步的方法,当vmax≤3m/s时,其最大回波频率fmax≈1.004×106,则其多普勒频移fd=1.004×106-1×106=4k,令混频信号的中频基准信号fc=fd=4k,则混频基准信号fr=1×106-4×103=996×103,混频后得到的中频信号f′=2fd=8k。
当流速为0时,混频后的多普勒频移信号为f′0=fd=4k,则包含流速信息的多普勒信号频移范围为0~4k,为模数转换芯片采样对象,经过混频的中频频率范围为4k~8k,为FFT算法频谱分析的对象。令FFT采样频率fs=4fd=16k,采样点数N=4096,则FFT算法频谱分析的频谱分辨率为且流量计的流速分辨率
Figure BDA0001779925480000127
可以看到,在保持模数转换芯片采样点数N不变,不增加FFT算法计算规模的情况下,相对于第一阶采样,第二阶采样的流速分辨率约增加了一倍,有效提高了低流速时流量计的分辨率和采样精度。
(8)根据FFT算法频谱分析的结果以及步骤(5)的方法,得到本次测量的流速v2,保持第二阶的采样参数不变,连续采样并求得平均值流速
Figure BDA0001779925480000128
依据公式(9)得到流量值Q。
进一步,如果两阶测量不能满足设备的精度需求,则可以在v2<1.5m/s
Figure BDA0001779925480000131
时,进入第三阶测量,进一步提高在更小流速下的测量精度和分辨率。
通过上述实施例可知,本发明提供的基于超声波的流量测量方法及流量计,至少实现了如下的有益效果:
1)在不增加硬件开销的前提下,提高了流量计在低流速情况下的测量精度。
2)通过对流量平均值的计算进一步提高了流量计算精度以及频谱分辨率。
虽然已经通过例子对本发明的一些特定实施例进行了详细说明,但是本领域的技术人员应该理解,以上例子仅是为了进行说明,而不是为了限制本发明的范围。本领域的技术人员应该理解,可在不脱离本发明的范围和精神的情况下,对以上实施例进行修改。本发明的范围由所附权利要求来限定。

Claims (9)

1.一种基于超声波的流量测量方法,其特征在于,包括:
步骤S1:设置超声波出射信号的频率f0
步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算FFT采样频率和混频基准信号fr
其中,采用以下公式计算所述FFT采样频率fs
fs=4fd,fd=fmax-f0
Figure FDA0002269553220000011
c为超声波在水中传播速度,θ为声波入射角,fd为多普勒频率,fmax为最大回波频率;
采用以下公式计算所述混频基准信号fr
fr=f0-fc,fc=fd,fc为中频基准频率;
步骤S3:对超声波回波信号进行模数采样处理,得到采样信号;
步骤S4:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱;
步骤S5:选取功率谱上幅值最大的点的频率作为计算多普勒频移信号f′1,则多普勒回波信号的频率f1为f1=fr+f′1
采用以下公式计算所述流速v,
步骤S6:检测v是否大于Q*vmax,其中Q小于1;
步骤S7:如果v>Q*vmax,则根据流速v计算流量;
步骤S8:如果v≤Q*vmax,则设置vmax=Q*vmax,返回步骤S2。
2.根据权利要求1所述的基于超声波的流量测量方法,其特征在于,
步骤S2还包括:设置测量次数x为初始值,其中,所述初始值为0;
步骤S7包括:
步骤S71:记录流速v至流速组,设置x=x+1;
步骤S72:判断x是否达到预定测量次数,其中,当x达到预定测量次数时,执行步骤S73,当未达到所述预定测量次数时,返回步骤S3;
步骤S73:计算所述流速组中所有v的平均值,并根据所述平均值计算所述流量。
3.根据权利要求1或2所述的基于超声波的流量测量方法,其特征在于,在步骤S3中,以1/2倍的所述FFT采样频率fs对超声波回波信号进行模数采样处理。
4.根据权利要求1或2所述的基于超声波的流量测量方法,其特征在于,在步骤S2之前,所述方法还包括:
步骤S9:向超声发射换能器发送PWM信号,所述PWM信号的频率为所述出射信号的频率。
5.根据权利要求1所述的基于超声波的流量测量方法,其特征在于,在步骤S7中,采用以下公式计算所述流量:
所述流量=流水截面积*v,其中流水截面积为已知量。
6.一种基于超声波的流量计,其特征在于,所述流量计包括:
中心处理模块以及分别与所述中心处理模块连接的超声接收换能器和模数转换芯片;
所述超声接收换能器用于接收超声波回波信号;
所述模数转换芯片用于对所述超声波回波信号进行模数采样;
所述中心处理模块用于执行以下步骤计算得到流速v:
步骤S1:设置超声波出射信号的频率f0
步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算混频基准信号fr以及FFT采样频率,其中,采用以下公式计算所述FFT采样频率fs
fs=4fd,fd=fmax-f0
Figure FDA0002269553220000021
c为超声波在水中传播速度,θ为声波入射角,fd为多普勒频率,fmax为最大回波频率;
采用以下公式计算所述混频基准信号fr
fr=f0-fc,fc=fd,fc为中频基准频率;
步骤S3:控制所述模数转换芯片对所述超声波回波信号进行模数采样处理,得到采样信号;
步骤S4:将所述采样信号与所述混频基准信号混频后进行频谱分析,得到功率谱;
步骤S5:选取功率谱上幅值最大的点的频率作为计算多普勒频移信号f′1,则多普勒回波信号的频率f1为f1=fr+f′1
采用以下公式计算所述流速v,
Figure FDA0002269553220000031
步骤S6:检测v是否大于Q*vmax,其中Q小于1;
步骤S7:如果v>Q*vmax,则根据v计算流量;
步骤S8:如果v≤Q*vmax,则设置vmax=Q*vmax,返回步骤S2。
7.根据权利要求6所述的基于超声波的流量计,其特征在于,所述中心处理模块执行的步骤S2还包括:设置测量次数x为初始值,其中,所述初始值为0;
执行的步骤S7还包括:
步骤S71:记录v至流速组,设置x=x+1;
步骤S72:判断x是否达到预定测量次数,其中,当x达到预定测量次数时,执行步骤S73,当未达到所述预定测量次数时,返回步骤S3;
步骤S73:计算所述流速组中所有v的平均值,并根据所述平均值计算所述流量。
8.根据权利要求6所述的基于超声波的流量计,其特征在于,还包括:与所述中心处理模块连接的超声发射换能器;
所述中心处理模块还用于在执行步骤S2之前,向所述超声发射换能器发送PWM信号,所述PWM信号的频率为所述出射信号的频率;
所述超声发射换能器用于根据所述PWM信号生成并发射所述超声波出射信号。
9.根据权利要求6所述的基于超声波的流量计,其特征在于,还包括:信号处理电路,所述信号处理电路连接在所述超声接收换能器与所述模数转换芯片之间,用于在所述模数转换芯片采样之前对所述超声波回波信号进行滤波和放大。
CN201810987006.0A 2018-08-28 2018-08-28 基于超声波的流量测量方法及流量计 Active CN109029602B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810987006.0A CN109029602B (zh) 2018-08-28 2018-08-28 基于超声波的流量测量方法及流量计

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810987006.0A CN109029602B (zh) 2018-08-28 2018-08-28 基于超声波的流量测量方法及流量计

Publications (2)

Publication Number Publication Date
CN109029602A CN109029602A (zh) 2018-12-18
CN109029602B true CN109029602B (zh) 2020-01-21

Family

ID=64625644

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810987006.0A Active CN109029602B (zh) 2018-08-28 2018-08-28 基于超声波的流量测量方法及流量计

Country Status (1)

Country Link
CN (1) CN109029602B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020183719A1 (ja) * 2019-03-14 2020-09-17 オムロン株式会社 流量測定装置
CN111650397A (zh) * 2020-04-02 2020-09-11 浙江清环智慧科技有限公司 超声波回波信号的处理方法、装置、流量计及存储介质
CN113108853B (zh) * 2021-04-08 2022-05-06 西南石油大学 一种提高低流速流体流量测量精度的方法
CN113687345B (zh) * 2021-07-23 2023-09-08 山东省科学院自动化研究所 一种多普勒雷达水流测速的方法及装置
CN115824331B (zh) * 2023-02-22 2023-06-16 青岛鼎信通讯科技有限公司 一种适用于超声水表的低功耗测量方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5777892A (en) * 1992-03-30 1998-07-07 Isco, Inc. Doppler shift velocity measuring system with correction factors
WO2006028990A2 (en) * 2004-09-03 2006-03-16 Nephros, Inc. Doppler flow measurement apparatus
CN104677437A (zh) * 2015-02-12 2015-06-03 延安大学 一种超声波液相流量精确化测量方法
CN105067058A (zh) * 2015-08-19 2015-11-18 上海航征测控系统有限公司 一种非接触式的排水管道流体流量的测量系统及方法

Also Published As

Publication number Publication date
CN109029602A (zh) 2018-12-18

Similar Documents

Publication Publication Date Title
CN109029602B (zh) 基于超声波的流量测量方法及流量计
CN100455999C (zh) 一种超声波测量液位的装置及方法
CN108802735B (zh) 一种用于未知声速环境的水下目标定位及测速方法和装置
US9689726B2 (en) Flow meter
US8234933B2 (en) Ultrasonic meter for determining physical quantity of a fluid based on a determined plurality of correlation coefficients
CN107576371A (zh) 一种超声波液位测量方法和超声波液位测量装置
US11480510B2 (en) Method and device for dynamically monitoring suspended matter based on annular interleaving array
JP3581676B2 (ja) 血流速度測定向け超音波診断装置
JP3016511B1 (ja) 超音波流速測定方法と装置
CN109579950B (zh) 一种气体超声波流量计的防错波检测装置
US20200200898A1 (en) Acoustic distance measuring circuit and method for low frequency modulated (lfm) chirp signals
CN106643939A (zh) 用于超声波流量计计算超声波传播时间的方法
JP2012509473A (ja) 超音波流量測定ユニットの測定変換器を校正するための方法及び装置
CN104764803A (zh) 基于超声波波长变化的材料应力检测技术
Tian et al. Energy peak fitting of echo based signal processing method for ultrasonic gas flow meter
CN105738469B (zh) 基于超声功率谱估计的液体悬浮物浓度测量方法与装置
CN202066962U (zh) 多功能超声波测量仪
CN103075981A (zh) 一种超声波测厚方法
Zhou et al. High-accuracy ultrasonic temperature measurement based on MLS-modulated continuous wave
CN108896664A (zh) 一种聚合物中声速与衰减系数一体化检测方法
CN110296913B (zh) 一种可燃粉尘扩散动态浓度的检测系统及其检测方法
CN203037849U (zh) 一种超声波测厚仪
CN111189513B (zh) 一种超声波液位测量方法、装置及超声波液位计
CN209043421U (zh) 超声波定距多点测量装置
WO2023130682A1 (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