CN111256786A - 一种判断气体超声流量计双阈值触发错波的方法 - Google Patents

一种判断气体超声流量计双阈值触发错波的方法 Download PDF

Info

Publication number
CN111256786A
CN111256786A CN202010079049.6A CN202010079049A CN111256786A CN 111256786 A CN111256786 A CN 111256786A CN 202010079049 A CN202010079049 A CN 202010079049A CN 111256786 A CN111256786 A CN 111256786A
Authority
CN
China
Prior art keywords
ultrasonic receiving
receiving signal
peak value
wave
threshold
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
CN202010079049.6A
Other languages
English (en)
Other versions
CN111256786B (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN202010079049.6A priority Critical patent/CN111256786B/zh
Publication of CN111256786A publication Critical patent/CN111256786A/zh
Application granted granted Critical
Publication of CN111256786B publication Critical patent/CN111256786B/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
    • G01F25/00Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
    • G01F25/10Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of flowmeters
    • 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

Landscapes

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

Abstract

本发明涉及一种判断气体超声流量计双阈值触发错波的方法,包括下列步骤:选取峰值;计算基准值;确定判断区间的比值;进行错波判定:判定为后触发错波,前触发错波还是判定为不错波;求取修正后的传播时间:根据错波判定结果修正超声波接收信号最终的传播时间。

Description

一种判断气体超声流量计双阈值触发错波的方法
技术领域
本发明属于流量测量技术领域,涉及气体超声流量计。
背景技术
在气体超声流量计的研究过程中,获取超声波在介质中的传播时间以及时间测量方案的精度是流量计精度的保证。现阶段超声波传播时间的测量方案主要有模拟方案与数字方案两种,模拟方案的电路设计相对简单,但准确度在大流速时难以得到保证。数字方案的测量结果准确度高,但是电路结构相对模拟方案更加复杂,方案的设计难度大。
双阈值触发法是模拟方案中一种应用于接收信号准确测时的方法。其原理是设定两个阈值,第一个阈值作为触发电平用以判断信号是否到达,第二个阈值即过零点作为信号准确到达的时刻。通过电路中的高精度计时芯片测量从单片机开始电平输出到接收信号触发第二个阈值的区间即超声信号的传播时间。为保证各流速下均可在相同位置触发,通常要与自动增益控制电路(AGC)相配合使用。当管道内气体的流速较大时,超声波接收信号衰减严重,使得第一个阈值的触发位置出现向前或向后错误触发的现象,造成实际检测到的超声波传播时间产生较大的测量误差,将该情况称为错波现象。
发明内容
本发明提供一种解决模拟方案中双阈值法存在的错波问题的方法,在此将其命名为峰值采集差值法,通过该方法实现了准确测时的目的,同时可有效解决自动控制电路的滞后性问题。技术方案如下:
一种判断气体超声流量计双阈值触发错波的方法,包括下列步骤:
步骤一:选取峰值;利用ADC采集超声波接收信号阈值触发后的五个峰值k1,k2,k3,k4,k5,并求前三个峰值之和,阈值的选择要保证所采集的前三个峰值未达到超声波接收信号最大峰值,但在第四个或第五个峰值达到超声波接收信号最大峰值,求和结果设为k:
k=k1+k2+k3
步骤二:计算基准值t:设定三个长度为10的数组A1,A2,A3,依次动态采集10次超声波接收信号,并将每一次超声波接收信号阈值触发之后的三个峰值k1,k2,k3对应存入数组A1,A2,A3,分别将数组A1,A2,A3中的10个元素相加并取平均值,得到每一个峰的平均峰值设为t1,t2,t3;将其求和,结果为基准值t:
t=t1+t2+t3
步骤三:确定判断区间的比值:根据超声波接收信号的规律以及峰值保持电路的设计,设定超声波接收信号的最大峰值为kp,并求出当次超声波接收信号最大峰值kp与静态下初始超声波接收信号最大峰值k0的比值,结果设为β:
Figure BDA0002379614370000021
步骤四:进行错波判定:
当k≥βt时,若k-βt>k4-k1,则判定为后触发错波,否则判定为不错波;
当k<βt时,若βt-k>k4-k1,则判定为前触发错波,否则判定为不错波;
步骤五:求取修正后的传播时间:根据错波判定结果修正超声波接收信号最终的传播时间;设双阈值法求得的传播时间为T,超声波接收信号频率为f,则修正后的传播时间T1为:
后错波时
Figure BDA0002379614370000022
前错波时
Figure BDA0002379614370000023
附图说明
图1双阈值触发法原理图
图2本发明实施例的双阈值法采集峰值示意图
图3本发明实施例的后错波情况下两数值分布图
具体实施方式
以下将以前述DN50双声道气体超声流量计为例,结合技术方案中的步骤方法,给出各步的参数设置和实施方法。测量系统使用MSP430FR6047单片机作为主控制器,利用其ADC12采集模块采集超声波接收信号相应峰值,并用峰值采集差值法进行错波判定。具体步骤如下:
步骤一:选取峰值。
(一)使用双阈值法判断并求出超声波的传播时间T,如图1所示。第一个触发阈值作为触发电平用以判断信号是否到达,第二个特征点检测阈值即过零点作为信号准确到达的时刻。硬件电路中超声波接收信号的零点基准值为1.6V,故第一个触发阈值设为2.2V,第二个特征点检测阈值设为1.6V。通过电路中的高精度计时芯片TDC-GP21测量从单片机开始电平输出到接受信号特征点检测阈值之间的时间即为超声波的传播时间T。
(二)利用单片机ADC12采集模块采集双阈值法判断后的超声波接收信号各个峰值,并将其命名为k1,k2,k3,k4,k5,根据峰值保持电路的原理可知,当超声波接收信号峰值点达到最大值时,最大峰值点以后的峰值点将保持与最大峰值点的大小一致,如图2所示。在软件程序中将前三个峰值相加,其值设为k:
k=k1+k2+k3
步骤二:计算基准值。
求出超声波接收信号的基准值t,其是判定当次超声波接收信号是否错波的标准。设定三个长度为10的数组A1,A2,A3,将每一次超声波接收信号触发之后的三个峰值k1,k2,k3对应存入数组A1,A2,A3。数组中存入的元素个数超过10时,去掉数组的第一个元素,其余元素向前循环一位,下一次信号的峰值作为数组的第十位元素。确保每次超声波接收信号到达时,数组中都有十个峰值元素。将数组A1,A2,A3中的10个元素相加并取平均值,得到每一个峰的平均峰值设为t1,t2,t3。将其求和,结果为t:
t=t1+t2+t3
步骤三:确定判断区间。
(一)求出当此超声波接收信号的最大峰值。根据超声波接收信号的规律以及峰值保持电路的设计,k5为当此超声波接收信号的最大峰值。
(二)找出静态下超声波接收信号的最大峰值。该峰值是由自动增益控制电路手动设定的已知量,设其值为k0
(三)求出当次超声波接收信号最大峰值k5与静态下初始超声波接收信号最大峰值k0的比值,结果设为β:
Figure BDA0002379614370000031
步骤四:进行错波判定。
当k≥βt时,若k-βt>k4-k1,则判定为后触发错波,否则判定为不错波。
当k<βt时,若βt-k>k4-k1,则判定为前触发错波,否则判定为不错波。
如图3所示为后错波情况下采集峰值之和k乘系数β之后的基准值t分布关系图。由图可知,k和βt可以准确地对错波现象进行判别。
步骤五:求取修正后的传播时间。
根据错波判定结果修正双阈值法求得的传播时间T。设超声波接收信号频率为f,则修正后的传播时间T1为:
后错波时
Figure BDA0002379614370000032
前错波时
Figure BDA0002379614370000041
为了验证本方法在气体超声波流量计中时间测量的准确性,设计了与双阈值触发法的对比实验,借助Labview控制NI6110数据采集板卡多次采集1m/s、3m/s、5m/s、10m/s、15m/s、20m/s、25m/s流速点下自动增益控制电路之后的超声波接收信号,将采集得到的超声波接收信号分别用双阈值触发法与峰值采集差值法计算超声波传播时间,每个流速点下处理1000个超声波接收信号,通过对比两种方法在不同流速下测得正确的期望传播时间的个数N来比较其优劣。处理结果如表1所示。
表1两种方法在不同流速下测得正确传播时间个数
Figure BDA0002379614370000042
通过表1可以看出,一旦管道内气体流速大于等于10m/s,由于信号幅值波动较剧烈且信号信噪比变小,双阈值触发法错波且随流速变大越来越严重。峰值采集差值法测得的正确期望传播时间远多于双阈值触发法,只会产生极少量的错误识别,而少量错误数据可以通过滤波进行滤除。
由以上对比分析可知,本发明提出的峰值采集差值法在不同流速下都能解决错波问题,实现准确测时,该方法具有灵活高效、抗干扰能力强、精度高等特点。

Claims (1)

1.一种判断气体超声流量计双阈值触发错波的方法,包括下列步骤:
步骤一:选取峰值;利用ADC采集超声波接收信号阈值触发后的五个峰值k1,k2,k3,k4,k5,并求前三个峰值之和,阈值的选择要保证所采集的前三个峰值未达到超声波接收信号最大峰值,但在第四个或第五个峰值达到超声波接收信号最大峰值。求和的结果设为k:
k=k1+k2+k3
步骤二:计算基准值t:设定三个长度为10的数组A1,A2,A3,依次动态采集10次超声波接收信号,并将每一次超声波接收信号阈值触发之后的三个峰值k1,k2,k3对应存入数组A1,A2,A3,分别将数组A1,A2,A3中的10个元素相加并取平均值,得到每一个峰的平均峰值设为t1,t2,t3;将其求和,结果为基准值t:
t=t1+t2+t3
步骤三:确定判断区间的比值:根据超声波接收信号的规律以及峰值保持电路的设计,设定超声波接收信号的最大峰值为kp,并求出当次超声波接收信号最大峰值kp与静态下初始超声波接收信号最大峰值k0的比值,结果设为β:
Figure FDA0002379614360000011
步骤四:进行错波判定:
当k≥βt时,若k-βt>k4-k1,则判定为后触发错波,否则判定为不错波;
当k<βt时,若βt-k>k4-k1,则判定为前触发错波,否则判定为不错波;
步骤五:求取修正后的传播时间:根据错波判定结果修正超声波接收信号最终的传播时间;设双阈值法求得的传播时间为T,超声波接收信号频率为f,则修正后的传播时间T1为:
后错波时
Figure FDA0002379614360000012
前错波时
Figure FDA0002379614360000013
CN202010079049.6A 2020-02-03 2020-02-03 一种判断气体超声流量计双阈值触发错波的方法 Active CN111256786B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010079049.6A CN111256786B (zh) 2020-02-03 2020-02-03 一种判断气体超声流量计双阈值触发错波的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010079049.6A CN111256786B (zh) 2020-02-03 2020-02-03 一种判断气体超声流量计双阈值触发错波的方法

Publications (2)

Publication Number Publication Date
CN111256786A true CN111256786A (zh) 2020-06-09
CN111256786B CN111256786B (zh) 2021-07-06

Family

ID=70949309

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010079049.6A Active CN111256786B (zh) 2020-02-03 2020-02-03 一种判断气体超声流量计双阈值触发错波的方法

Country Status (1)

Country Link
CN (1) CN111256786B (zh)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5639971A (en) * 1996-10-04 1997-06-17 Dieterich Technology Holding Corp. Method and apparatus for detecting a signal
JP2003014515A (ja) * 2001-06-28 2003-01-15 Fuji Electric Co Ltd 超音波流量計
CN1846178A (zh) * 2003-07-03 2006-10-11 丹尼尔工业公司 用于渡越时间超声计量仪的峰值切换检测器
CN101886939A (zh) * 2010-06-10 2010-11-17 宁波大学 一种时差法超声流量计静态漂移抑制模型及抑制方法
CN102636252A (zh) * 2012-04-10 2012-08-15 吉林大学 一种超声波到达精确时刻检测的方法及装置
CN104897248A (zh) * 2015-06-12 2015-09-09 天津大学 一种噪声背景下准确估计超声流量计传播时间方法
US20150323374A1 (en) * 2014-05-07 2015-11-12 Azbil Corporation Ultrasonic flow meter and ultrasound absorbing body fault evaluating method
CN105300508A (zh) * 2015-10-30 2016-02-03 北京控制工程研究所 一种超声波信号传播时间的双阈值检测电路及检测方法
CN105628115A (zh) * 2015-12-30 2016-06-01 浙江大学 一种应用于时差式超声波流量计的脉冲噪声滤除方法
CN105698886A (zh) * 2016-05-04 2016-06-22 湖北锐意自控系统有限公司 一种基于超声波检测技术的气体流量检测方法
JP2017116458A (ja) * 2015-12-25 2017-06-29 アズビル株式会社 超音波流量計
CN109579950A (zh) * 2018-12-07 2019-04-05 天津商业大学 一种气体超声波流量计的防错波检测装置
CN109612541A (zh) * 2018-12-07 2019-04-12 天津商业大学 一种气体超声波流量计的防错波信号处理电路
CN109931996A (zh) * 2019-02-26 2019-06-25 天津大学 一种准确测量气体超声流量计信号传播时间的方法
CN112304376A (zh) * 2020-10-27 2021-02-02 浙江大学 基于数据融合的超声波流量计流量测量方法

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5639971A (en) * 1996-10-04 1997-06-17 Dieterich Technology Holding Corp. Method and apparatus for detecting a signal
JP2003014515A (ja) * 2001-06-28 2003-01-15 Fuji Electric Co Ltd 超音波流量計
CN1846178A (zh) * 2003-07-03 2006-10-11 丹尼尔工业公司 用于渡越时间超声计量仪的峰值切换检测器
CN101886939A (zh) * 2010-06-10 2010-11-17 宁波大学 一种时差法超声流量计静态漂移抑制模型及抑制方法
CN102636252A (zh) * 2012-04-10 2012-08-15 吉林大学 一种超声波到达精确时刻检测的方法及装置
US20150323374A1 (en) * 2014-05-07 2015-11-12 Azbil Corporation Ultrasonic flow meter and ultrasound absorbing body fault evaluating method
CN104897248A (zh) * 2015-06-12 2015-09-09 天津大学 一种噪声背景下准确估计超声流量计传播时间方法
CN105300508A (zh) * 2015-10-30 2016-02-03 北京控制工程研究所 一种超声波信号传播时间的双阈值检测电路及检测方法
JP2017116458A (ja) * 2015-12-25 2017-06-29 アズビル株式会社 超音波流量計
CN105628115A (zh) * 2015-12-30 2016-06-01 浙江大学 一种应用于时差式超声波流量计的脉冲噪声滤除方法
CN105698886A (zh) * 2016-05-04 2016-06-22 湖北锐意自控系统有限公司 一种基于超声波检测技术的气体流量检测方法
CN109579950A (zh) * 2018-12-07 2019-04-05 天津商业大学 一种气体超声波流量计的防错波检测装置
CN109612541A (zh) * 2018-12-07 2019-04-12 天津商业大学 一种气体超声波流量计的防错波信号处理电路
CN109931996A (zh) * 2019-02-26 2019-06-25 天津大学 一种准确测量气体超声流量计信号传播时间的方法
CN112304376A (zh) * 2020-10-27 2021-02-02 浙江大学 基于数据融合的超声波流量计流量测量方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘泽忠 等: ""超声波流量计自适应双门限触发法研究"", 《化工自动化及仪表》 *
郑丹丹 等: ""速度分布对气体超声流量计声传播规律的影响"", 《天津大学学报(自然科学与工程技术版)》 *
闵闿: ""基于双阈值电平比较法的时差法超声波流量计信号处理研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
陈洁 等: ""基于双阈值比较法超声波流量计信号处理"", 《电子测量与仪器学报》 *

Also Published As

Publication number Publication date
CN111256786B (zh) 2021-07-06

Similar Documents

Publication Publication Date Title
CN108548578B (zh) 一种基于自适应阈值的超声波回波信号特征峰识别方法
JPS62135783A (ja) 増分自動ゲイン・コントロ−ル付レ−ダ受信機及び受信方法
CN109579950B (zh) 一种气体超声波流量计的防错波检测装置
CN101179283A (zh) 脉冲噪声的检测及抑制方法
TW200919986A (en) Jamming detector and jamming detecting method
EP2912406A1 (en) Ultrasonic measurement apparatus and method
CN107271973B (zh) 韦布尔杂波环境下基于偏斜度和均值比的恒虚警检测方法
CN111256786B (zh) 一种判断气体超声流量计双阈值触发错波的方法
US7504884B2 (en) Automatic gain control apparatus of cable modem, and lock-time acquisition method using the same
CN112526602A (zh) 一种基于长短时窗和ar模型方差激增效应的p波到时拾取方法
CN110118979A (zh) 基于广义互熵的改进差分进化算法估计多径参数的方法
CN110299926B (zh) 一种面向低信噪比环境的水声信号检测方法
CN112332807A (zh) 一种弱包络信号检测方法和系统
US20170059718A1 (en) Signal processing device and radiation measurement device
RU2684440C1 (ru) Способ отождествления объектов, обнаруженных несколькими системами
CN109117698B (zh) 一种基于最小均方误差准则的噪声背景估计方法
EP3831090A1 (en) Sound sensing with time-varying thresholds
CN109931996B (zh) 一种准确测量气体超声流量计信号传播时间的方法
CN115755183A (zh) 一种基于尾波干涉的震源位置追踪方法
CN115508788A (zh) 基于稳健噪底估计的信道化检测方法
CN108627857B (zh) 多径检测方法、装置及gnss接收机
CN111879368B (zh) 一种流量测量中相移检测方法、装置及设备
CN114674384A (zh) 超声波流量计的防错波检测方法、装置、设备及流量计
CN110032758B (zh) 计算电信号的能量的方法、装置和计算机存储介质
CN107346030A (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