CN109470320A - 一种超声波水表阶跃滤波方法 - Google Patents

一种超声波水表阶跃滤波方法 Download PDF

Info

Publication number
CN109470320A
CN109470320A CN201811296784.1A CN201811296784A CN109470320A CN 109470320 A CN109470320 A CN 109470320A CN 201811296784 A CN201811296784 A CN 201811296784A CN 109470320 A CN109470320 A CN 109470320A
Authority
CN
China
Prior art keywords
time difference
current time
difference value
water meter
ultrasonic
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
CN201811296784.1A
Other languages
English (en)
Other versions
CN109470320B (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.)
NANJING LINYANG ELECTRIC POWER TECHNOLOGY CO LTD
Jiangsu Linyang Solarfun Co Ltd
Original Assignee
NANJING LINYANG ELECTRIC POWER TECHNOLOGY CO LTD
Jiangsu Linyang Solarfun 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 NANJING LINYANG ELECTRIC POWER TECHNOLOGY CO LTD, Jiangsu Linyang Solarfun Co Ltd filed Critical NANJING LINYANG ELECTRIC POWER TECHNOLOGY CO LTD
Priority to CN201811296784.1A priority Critical patent/CN109470320B/zh
Publication of CN109470320A publication Critical patent/CN109470320A/zh
Application granted granted Critical
Publication of CN109470320B publication Critical patent/CN109470320B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/667—Arrangements of transducers for ultrasonic flowmeters; Circuits for operating ultrasonic flowmeters
    • G01F1/668—Compensating or correcting for variations in velocity of sound

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Measurement Of Levels Of Liquids Or Fluent Solid Materials (AREA)

Abstract

一种超声波水表阶跃滤波方法，该方法对超声波水表的顺、逆流时差值进行首波跃变滤波，首先，获取超声波在水中逆流传播时间t逆和超声波在水中顺流传播时间t顺，计算当前时刻的顺、逆流时差值Δti；其次，设置最大时差阈值ΔTmax，以及相邻时差的差值阈值ΔTdmax；对当前时差值和相邻时差值进行判断，循环查找首波跃变阶数n；对当前时刻的时差值进行滤波。本发明弥补了普通滤波算法的不足，能够兼顾温度反算与正反向同时计量的情况；有效的减少因杂质、气泡、污垢、过载流量导致的超声波首波跃变、幅值下降，以致计量精度受到影响。

Description

一种超声波水表阶跃滤波方法
技术领域
本发明属超声波水表计量技术领域，特别涉及一种超声波水表阶跃滤波方法。
背景技术
超声波换能器浸泡在水中结垢后对换能器接收信号幅值影响特别大，当水中混杂气泡时，超声波在不同介质中的传播速度受影响，换能器的接收信号幅值也会有较大变化，这就会间接导致超声波接收信号发生整体错位，即首波发生跃变，严重影响计量精确度。所以对于长期使用于计量的超声波水表，既要保证计量精度，又要保证使用可靠性，保障用户利益。要保证在使用年限内的计量准确度，急需一种能够避免首波跃变保证计量精度的方法。
超声波水表基表由于加工及系统环境因素，静水下顺逆流时间并不能完全遵循t顺≤t逆的关系，采样到的顺流时间t顺与逆流时间t逆在一定范围内跳动，而时差Δt将会存在正负交替变化。因此，使用t顺≤t逆的判断依据进行首波跃变滤波将会带来问题；同时，反向流计量与正向计量同时存在的超声波水表越来越广泛应用，因此无法使用以上判据进行滤波。
使用跃变调整首波检测阈值的方法，能够有效的解决因首波跃变带来的计量误差，但存在一个较为严重的缺陷，即当增加或减小首波检测阈值后，超声波传播顺逆流传播绝对时间将变大或变小。虽然采用时差法原理计量对流速并无影响，但对于无温度传感器超声波水表而言，此方法无法有效解决温度反算问题。因此，一种针对无温度传感器超声波水表正向、反向计量的首波跃变滤波解决方法呼之欲出。
首波跃变现象非常普遍，因此选择一种合适的过滤方法至关重要。
发明内容
本发明目的是为弥补现有技术缺陷，提供一种简单易懂，节省硬件成本，保证计量精度的超声波水表阶跃滤波方法。
本发明的技术方案是：
一种超声波水表阶跃滤波方法，该方法对超声波水表的顺、逆流时差值进行首波跃变滤波。
进一步地，首波跃变滤波包括以下步骤：
S1、获取超声波在水中逆流传播时间t逆和超声波在水中顺流传播时间t顺，计算当前时刻的顺、逆流时差值Δti，填入时差FIFO缓冲区；
S2、设置最大时差阈值ΔTmax，以及相邻时差的差值阈值ΔTdmax；
S3、对当前时差值和相邻时差值按照下述公式进行判断，当其中之一符合时，发生首波跃变；
其中，i表示当前时刻，i-1表示当前时刻之前一时刻，Δti-1表示当前时刻之前一时刻的顺、逆流时差值；
S4、使用迭代循环检测方法，循环查找首波跃变阶数n；
S5、根据前述首波跃变阶数n，对当前时刻的时差值进行滤波，采用下述公式获取滤波后的当前时差值Δti′；
Δt′i＝Δti+n×To
其中：To表示一个波形宽度的时间，f表示超声波水表的换能器中心频率。
进一步地，步骤S1中，对当前时刻的顺、逆流时差值Δti进行零漂补偿处理，得到校准时差，之后填入时差FIFO缓冲区。
进一步地，步骤S4中，首波跃变阶数n的取值范围是n∈(-5，5)。
进一步地，步骤S4中，循环查找首波跃变阶数n的具体方法为：对于当前时刻i，从n＝-5开始进行查找，当满足以下全部条件时，当前n的取值为首波跃变阶数；
本发明的有益效果：
本发明弥补了普通滤波算法的不足，能够兼顾温度反算与正反向同时计量的情况。
本发明有效的减少因杂质、气泡、污垢、过载流量导致的超声波首波跃变、幅值下降，以致计量精度受到影响。
本发明从算法角度改进计量精度，实时简单，节省硬件成本，提升超声水表计量可靠性。
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
图1是静水原始时差散点图。
图2是滤波前时差散点图。
图3是滤波后时差散点图。
具体实施方式
下面将参照附图更详细地描述本发明的优选实施方式。虽然附图中显示了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
一种超声波水表阶跃滤波方法，该方法对超声波水表的顺、逆流时差值进行首波跃变滤波，包括以下步骤：
S1、获取超声波在水中逆流传播时间t逆和超声波在水中顺流传播时间t顺，计算当前时刻的顺、逆流时差值Δti，进行零漂补偿处理，得到校准时差，之后填入时差FIFO缓冲区；
S2、设置最大时差阈值ΔTmax，以及相邻时差的差值阈值ΔTdmax；
S3、对当前时差值和相邻时差值按照下述公式进行判断，当其中之一符合时，发生首波跃变；
其中，i表示当前时刻，i-1表示当前时刻之前一时刻，Δti-1表示当前时刻之前一时刻的顺、逆流时差值；
S4、使用迭代循环检测方法，循环查找首波跃变阶数n；
首波跃变阶数n的取值范围是n∈(-5，5)，循环查找首波跃变阶数n的具体方法为：对于当前时刻i，从n＝-5开始进行查找，当满足以下全部条件时，当前n的取值为首波跃变阶数；
S5、根据前述首波跃变阶数n，对当前时刻的时差值进行滤波，采用下述公式获取滤波后的当前时差值Δti′；
Δt′i＝Δti+n×To
其中：To表示一个波形宽度的时间，f表示超声波水表的换能器中心频率。
本滤波方法既能满足时差因首波跃变的影响带来的误差，同时又能满足反向计量的首波跃变滤波处理，不区分正反向计量，同理，顺流与逆流时间可通过此滤波方法滤波，达到温度反算的目的与准确度，本发明简单明了，无需任何前提条件，使用简单可靠，能够有效去除首波跃变影响，节省硬件成本，提高产品质量。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (5)

1.一种超声波水表阶跃滤波方法，其特征在于：该方法对超声波水表的顺、逆流时差值进行首波跃变滤波。
2.根据权利要求1所述的一种超声波水表阶跃滤波方法，其特征在于首波跃变滤波包括以下步骤：
S1、获取超声波在水中逆流传播时间t逆和超声波在水中顺流传播时间t顺，计算当前时刻的顺、逆流时差值Δti，填入时差FIFO缓冲区；
S2、设置最大时差阈值ΔTmax，以及相邻时差的差值阈值ΔTdmax；
S3、对当前时差值和相邻时差值按照下述公式进行判断，当其中之一符合时，发生首波跃变；
其中，i表示当前时刻，i-1表示当前时刻之前一时刻，Δti-1表示当前时刻之前一时刻的顺、逆流时差值；
S4、使用迭代循环检测方法，循环查找首波跃变阶数n；
S5、根据前述首波跃变阶数n，对当前时刻的时差值进行滤波，采用下述公式获取滤波后的当前时差值Δti′；
Δt′i＝Δti+n×To
其中：To表示一个波形宽度的时间，f表示超声波水表的换能器中心频率。
3.根据权利要求2所述的一种超声波水表阶跃滤波方法，其特征在于步骤S1中，对当前时刻的顺、逆流时差值Δti进行零漂补偿处理，得到校准时差，之后填入时差FIFO缓冲区。
4.根据权利要求2所述的一种超声波水表阶跃滤波方法，其特征在于步骤S4中，首波跃变阶数n的取值范围是n∈(-5，5)。
5.根据权利要求4所述的一种超声波水表阶跃滤波方法，其特征在于步骤S4中，循环查找首波跃变阶数n的具体方法为：对于当前时刻i，从n＝-5开始进行查找，当满足以下全部条件时，当前n的取值为首波跃变阶数；
CN201811296784.1A 2018-11-01 2018-11-01 一种超声波水表阶跃滤波方法 Active CN109470320B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811296784.1A CN109470320B (zh) 2018-11-01 2018-11-01 一种超声波水表阶跃滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811296784.1A CN109470320B (zh) 2018-11-01 2018-11-01 一种超声波水表阶跃滤波方法

Publications (2)

Publication Number Publication Date
CN109470320A true CN109470320A (zh) 2019-03-15
CN109470320B CN109470320B (zh) 2020-11-06

Family

ID=65666417

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811296784.1A Active CN109470320B (zh) 2018-11-01 2018-11-01 一种超声波水表阶跃滤波方法

Country Status (1)

Country Link
CN (1) CN109470320B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113010141A (zh) * 2021-03-04 2021-06-22 湖北三赫智能科技有限公司 一种信号首波信息提取方法、系统、装置及可读存储介质
CN113865659A (zh) * 2021-09-17 2021-12-31 湖南威铭能源科技有限公司 一种高精度户用超声波水表抗干扰滤波方法及设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003194604A (ja) * 2001-12-26 2003-07-09 Babcock Hitachi Kk 音波式ガス温度・流量計とガス中音波の伝播時間の検出方法
CN101886939A (zh) * 2010-06-10 2010-11-17 宁波大学 一种时差法超声流量计静态漂移抑制模型及抑制方法
CN104034375A (zh) * 2014-06-27 2014-09-10 苏州东剑智能科技有限公司 用于采用超声波进行流量测量的方法
WO2015055209A1 (en) * 2013-10-14 2015-04-23 Flonidan A/S Flow meter algorithm
JP2018138891A (ja) * 2017-02-24 2018-09-06 株式会社東芝 超音波流量計

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003194604A (ja) * 2001-12-26 2003-07-09 Babcock Hitachi Kk 音波式ガス温度・流量計とガス中音波の伝播時間の検出方法
CN101886939A (zh) * 2010-06-10 2010-11-17 宁波大学 一种时差法超声流量计静态漂移抑制模型及抑制方法
WO2015055209A1 (en) * 2013-10-14 2015-04-23 Flonidan A/S Flow meter algorithm
CN104034375A (zh) * 2014-06-27 2014-09-10 苏州东剑智能科技有限公司 用于采用超声波进行流量测量的方法
JP2018138891A (ja) * 2017-02-24 2018-09-06 株式会社東芝 超音波流量計

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
殷光: "超声波流量测量技术研究", 《中国优秀硕士学位论文全文数据库》 *
王慧: "时差法超声波流量计数据处理方法的研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113010141A (zh) * 2021-03-04 2021-06-22 湖北三赫智能科技有限公司 一种信号首波信息提取方法、系统、装置及可读存储介质
CN113865659A (zh) * 2021-09-17 2021-12-31 湖南威铭能源科技有限公司 一种高精度户用超声波水表抗干扰滤波方法及设备
CN113865659B (zh) * 2021-09-17 2024-03-22 湖南威铭能源科技有限公司 一种高精度户用超声波水表抗干扰滤波方法及设备

Also Published As

Publication number Publication date
CN109470320B (zh) 2020-11-06

Similar Documents

Publication Publication Date Title
CN105841762B (zh) 超声波水表的流量计量方法和系统
CN106643939B (zh) 用于超声波流量计计算超声波传播时间的方法
CN108548578B (zh) 一种基于自适应阈值的超声波回波信号特征峰识别方法
CN104501889B (zh) 基于互相关时差法超声波流量的检测方法
WO2014205884A1 (zh) 一种基于零点分析的时差式超声波流量计测量方法
CN109470320A (zh) 一种超声波水表阶跃滤波方法
CN112484798B (zh) 一种基于改进卡尔曼滤波器的时差超声波流量测量方法
JP6709846B2 (ja) メーターの相互オペレーションのための超音波フローメーター自動調整
CN107860430B (zh) 一种基于时差法的超声波气体流量计时间差测量方法
US20210003436A1 (en) Time-of-flight generating circuit and chip, flow meter and method of the same
CN106706056A (zh) 一种大口径超声波水表流量测量的补偿方法
CN115824331B (zh) 一种适用于超声水表的低功耗测量方法
CN110346004A (zh) 一种双声道超声时差法的流量测量数据融合方法
US6816808B2 (en) Peak switch detector for transit time ultrasonic meters
CN112304376B (zh) 基于数据融合的超声波流量计流量测量方法
KR101764870B1 (ko) 초음파 유량계의 신호처리시스템
CN205246150U (zh) 一种自动校正零点漂移的超声水表
CN105785123A (zh) 一种基于apFFT相位差的雷达信号频率计算方法
CN111964738B (zh) 一种基于超声回波信号相似度的动态阈值设置方法
CN111473827B (zh) V形声道零飘消除方法
JPH11304559A (ja) 流量計測装置
CN106643940B (zh) 基于回波能量计算超声波流量计传播时间的方法
CN104965401A (zh) 一种多点定位相关监视系统的到达时间测量方法及系统
CN112147366A (zh) 用于使用超声波来测量流体速度的方法
CN104713631B (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