CN109029602A - 基于超声波的流量测量方法及流量计 - Google Patents
基于超声波的流量测量方法及流量计 Download PDFInfo
- Publication number
- CN109029602A CN109029602A CN201810987006.0A CN201810987006A CN109029602A CN 109029602 A CN109029602 A CN 109029602A CN 201810987006 A CN201810987006 A CN 201810987006A CN 109029602 A CN109029602 A CN 109029602A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- flow
- max
- flow velocity
- 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
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
- G01F1/663—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 by measuring Doppler frequency shift
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,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,
一种基于超声波的流量计,所述流量计包括:
中心处理模块以及分别与所述中心处理模块连接的超声接收换能器和模数转换芯片;
所述超声接收换能器用于接收超声波回波信号;
所述模数转换芯片用于对所述超声波回波信号进行模数采样;
所述中心处理模块用于执行以下步骤计算得到流速v:
步骤S1:设置超声波出射信号的频率f0;
步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算混频基准信号fr以及FFT采样频率,其中,采用以下公式计算所述FFT采样频率fs:
fs=4fd,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=4fd,c为超声波在水中传播速度,θ为声波入射角;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,
需要说明的是:则在超声波多普勒流量计的测量范围0~vmax内,包含流速信息的多普勒频移范围0~fd,即为ADC需要采样的对象;混频后的中频信号范围为fd~2fd,即为FFT采样及算法分析的对象。
根据采样定理,FFT采样频率至少为被采样信号频率的两倍,则令FFT的采样频率fs=4fd,采样点数为N,则FFT算法进行频谱分析的频谱分辨率为假设根据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,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的平均值。
具体地,采用以下公式计算所述平均值
具体地,x为测量次数,测量次数可以达到预定测量次数,预定测量次数可以根据流量测量精度通过标定确定,比如,预定测量次数为:10~100之间的值。
步骤210:根据所述平均值计算流量。
具体地,采用以下公式计算所述流量:
所述流量=流水截面积*其中流水截面积为已知量。
步骤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,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,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°,则根据公式可以得到回波信号的最大频率为:
其中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算法频谱分析的频谱分辨率为根据公式可以得到多普勒超声波流量计的流速分辨率vh:
(5)假设经过FFT算法的频谱分析,得到功率谱上功率最大点的频率为f′1=15k,则多普勒回波信号频率为f1=fr+f′1=992×103+15×103=1.007×106,可以得到流速v1:
由于则重复以上步骤,采样计算x次,得到平均流速值根据公式流量=流水截面积*即可以得到测量的流量值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算法频谱分析的频谱分辨率为且流量计的流速分辨率
可以看到,在保持模数转换芯片采样点数N不变,不增加FFT算法计算规模的情况下,相对于第一阶采样,第二阶采样的流速分辨率约增加了一倍,有效提高了低流速时流量计的分辨率和采样精度。
(8)根据FFT算法频谱分析的结果以及步骤(5)的方法,得到本次测量的流速v2,保持第二阶的采样参数不变,连续采样并求得平均值流速依据公式(9)得到流量值Q。
进一步,如果两阶测量不能满足设备的精度需求,则可以在v2<1.5m/s时,进入第三阶测量,进一步提高在更小流速下的测量精度和分辨率。
通过上述实施例可知,本发明提供的基于超声波的流量测量方法及流量计,至少实现了如下的有益效果:
1)在不增加硬件开销的前提下,提高了流量计在低流速情况下的测量精度。
2)通过对流量平均值的计算进一步提高了流量计算精度以及频谱分辨率。
虽然已经通过例子对本发明的一些特定实施例进行了详细说明,但是本领域的技术人员应该理解,以上例子仅是为了进行说明,而不是为了限制本发明的范围。本领域的技术人员应该理解,可在不脱离本发明的范围和精神的情况下,对以上实施例进行修改。本发明的范围由所附权利要求来限定。
Claims (10)
1.一种基于超声波的流量测量方法,其特征在于,包括:
步骤S1:设置超声波出射信号的频率f0;
步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算FFT采样频率和混频基准信号fr,
其中,采用以下公式计算所述FFT采样频率fs:
fs=4fd,fd=fmax-f0,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。
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.根据权利要求1所述的基于超声波的流量测量方法,其特征在于,在步骤S5中,根据所述功率谱计算计算所述流速v包括:
选取功率谱上幅值最大的点的频率作为计算多普勒频移信号f′1,则多普勒回波信号的频率f1为f1=fr+f′1;
采用以下公式计算所述流速v,
7.一种基于超声波的流量计,其特征在于,所述流量计包括:
中心处理模块以及分别与所述中心处理模块连接的超声接收换能器和模数转换芯片;
所述超声接收换能器用于接收超声波回波信号;
所述模数转换芯片用于对所述超声波回波信号进行模数采样;
所述中心处理模块用于执行以下步骤计算得到流速v:
步骤S1:设置超声波出射信号的频率f0;
步骤S2:根据所述出射信号的频率f0和超声波流量计当前可测量的最大流速vmax计算混频基准信号fr以及FFT采样频率,其中,采用以下公式计算所述FFT采样频率fs:
fs=4fd,fd=fmax-f0,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。
8.根据权利要求7所述的基于超声波的流量计,其特征在于,所述中心处理模块执行的步骤S2还包括:设置测量次数x为初始值,其中,所述初始值为0;
执行的步骤S7还包括:
步骤S71:记录v至流速组,设置x=x+1;
步骤S72:判断x是否达到预定测量次数,其中,当x达到预定测量次数时,执行步骤S73,当未达到所述预定测量次数时,返回步骤S3;
步骤S73:计算所述流速组中所有v的平均值,并根据所述平均值计算所述流量。
9.根据权利要求7所述的基于超声波的流量计,其特征在于,还包括:与所述中心处理模块连接的超声发射换能器;
所述中心处理模块还用于在执行步骤S2之前,向所述超声发射换能器发送PWM信号,所述PWM信号的频率为所述出射信号的频率;
所述超声发射换能器用于根据所述PWM信号生成并发射所述超声波出射信号。
10.根据权利要求7所述的基于超声波的流量计,其特征在于,还包括:信号处理电路,所述信号处理电路连接在所述超声接收换能器与所述模数转换芯片之间,用于在所述模数转换芯片采样之前对所述超声波回波信号进行滤波和放大。
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 true CN109029602A (zh) | 2018-12-18 |
CN109029602B 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) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111650397A (zh) * | 2020-04-02 | 2020-09-11 | 浙江清环智慧科技有限公司 | 超声波回波信号的处理方法、装置、流量计及存储介质 |
CN112639412A (zh) * | 2019-03-14 | 2021-04-09 | 欧姆龙株式会社 | 流量测量装置 |
CN113108853A (zh) * | 2021-04-08 | 2021-07-13 | 西南石油大学 | 一种提高低流速流体流量测量精度的方法 |
CN113687345A (zh) * | 2021-07-23 | 2021-11-23 | 山东省科学院自动化研究所 | 一种多普勒雷达水流测速的方法及装置 |
CN115824331A (zh) * | 2023-02-22 | 2023-03-21 | 青岛鼎信通讯科技有限公司 | 一种适用于超声水表的低功耗测量方法 |
Citations (4)
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 | 上海航征测控系统有限公司 | 一种非接触式的排水管道流体流量的测量系统及方法 |
-
2018
- 2018-08-28 CN CN201810987006.0A patent/CN109029602B/zh active Active
Patent Citations (4)
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 | 上海航征测控系统有限公司 | 一种非接触式的排水管道流体流量的测量系统及方法 |
Non-Patent Citations (2)
Title |
---|
曹新亮: "超声波流速测量的多普勒频移提取方法研究", 《应用声学》 * |
邱立存: "采用频谱分析技术的多普勒超声波流量计的研究", 《微计算机信息》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112639412A (zh) * | 2019-03-14 | 2021-04-09 | 欧姆龙株式会社 | 流量测量装置 |
CN111650397A (zh) * | 2020-04-02 | 2020-09-11 | 浙江清环智慧科技有限公司 | 超声波回波信号的处理方法、装置、流量计及存储介质 |
CN113108853A (zh) * | 2021-04-08 | 2021-07-13 | 西南石油大学 | 一种提高低流速流体流量测量精度的方法 |
CN113108853B (zh) * | 2021-04-08 | 2022-05-06 | 西南石油大学 | 一种提高低流速流体流量测量精度的方法 |
CN113687345A (zh) * | 2021-07-23 | 2021-11-23 | 山东省科学院自动化研究所 | 一种多普勒雷达水流测速的方法及装置 |
CN113687345B (zh) * | 2021-07-23 | 2023-09-08 | 山东省科学院自动化研究所 | 一种多普勒雷达水流测速的方法及装置 |
CN115824331A (zh) * | 2023-02-22 | 2023-03-21 | 青岛鼎信通讯科技有限公司 | 一种适用于超声水表的低功耗测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109029602B (zh) | 2020-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109029602A (zh) | 基于超声波的流量测量方法及流量计 | |
Garbini et al. | Measurement of fluid turbulence based on pulsed ultrasound techniques. Part 2. Experimental investigation | |
CN106643939B (zh) | 用于超声波流量计计算超声波传播时间的方法 | |
JPWO2005083372A1 (ja) | パルスドップラ方式と伝搬時間差方式の両方式対応型超音波流量計、同流量計において測定方式を自動選択する方法およびプログラム、同流量計用の電子装置 | |
JPH11505020A (ja) | 空気データ装置 | |
CN106706056B (zh) | 一种大口径超声波水表流量测量的补偿方法 | |
CN109798944A (zh) | 基于微机械压电超声波换能器的流量计及渡越时间测量装置 | |
CN108917865B (zh) | 一种多普勒流速方向识别装置和控制方法 | |
WO2018233198A1 (zh) | 一种基于环状交织阵列的悬浮物动态监测方法与装置 | |
CN109579950A (zh) | 一种气体超声波流量计的防错波检测装置 | |
WO1988008516A1 (en) | Ultrasonic fluid flowmeter | |
Tian et al. | Energy peak fitting of echo based signal processing method for ultrasonic gas flow meter | |
CN104764522A (zh) | 一种超声功率测量方法及装置 | |
CN202066962U (zh) | 多功能超声波测量仪 | |
CN209166556U (zh) | 基于微机械压电超声波换能器的流量计渡越时间测量装置 | |
US5113867A (en) | Fluid flow characterizing | |
WO2005083371A1 (ja) | ドップラー式超音波流量計 | |
CN101793908A (zh) | 超声波烟气流速计量仪 | |
Haugwitz et al. | Multipath flow metering of high-velocity gas using ultrasonic phased-arrays | |
CN109001487A (zh) | 一种基于相敏检波时延测量的声波测风系统及方法 | |
Gerasimov et al. | Applications of digital signal processing methods in ultrasonic flowmeters | |
CN109632025A (zh) | 针对声学多普勒流量计的流速的数据处理方法 | |
JP2010256075A (ja) | 流量計及び流量計測方法 | |
JP4904099B2 (ja) | パルス状信号の伝搬時間測定装置及び超音波式流量測定装置 | |
Padilla et al. | Acoustic backscattering observations from non-spherical gas bubbles with ka between 0.03 and 4.4 |
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 |