CN110646042A - 一种用于低功耗超声流量计飞行时间差计算的互相关插值方法 - Google Patents
一种用于低功耗超声流量计飞行时间差计算的互相关插值方法 Download PDFInfo
- Publication number
- CN110646042A CN110646042A CN201910981870.4A CN201910981870A CN110646042A CN 110646042 A CN110646042 A CN 110646042A CN 201910981870 A CN201910981870 A CN 201910981870A CN 110646042 A CN110646042 A CN 110646042A
- Authority
- CN
- China
- Prior art keywords
- cross
- flight
- correlation
- low
- time difference
- 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.)
- Pending
Links
Images
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
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
- G05B19/042—Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
- G05B19/0423—Input/output
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Electromagnetism (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种用于低功耗超声流量计飞行时间差计算的互相关插值方法。在进行流量测量时,由于低功耗嵌入式***的数字模拟转换器(ADC)采集回波信号的采样频率不够高导致采集到的回波信号的时间分辨率有限,进而导致计算互相关函数后得到的飞行时间差分辨率过低,使得气体超声波流量计的测量精度不高。本发明提出了一种用于飞行时间差计算的互相关插值方法。该发明对计算得到的互相关函数进行插值处理,可以极大地提高互相关函数的时间分辨率。本发明有效地解决了基于互相关法计算飞行时间差的低功耗超声流量计采样频率不够高的问题,提高了飞行时间差的测量精度和分辨力。
Description
技术领域
本发明涉及流量测量领域,尤其涉及一种用于低功耗超声流量计飞行时间差计算的方法。
背景技术
随着电子技术和微处理器技术的快速发展,超声波测量在过去的十年中取得了显着的进步,尤其是碳氢化合物的测量。超声流量计作为一种非接触式流量计,与传统的计量设备相比具有许多优点:不同管径的兼容性,无流阻,无压力损失,测量精度高,功耗低等。
超声气体流量计通常安装在沙漠等野外,必须具有良好的抗干扰能力。它通常有一个独立的电源,如电池。因此,低功耗设计对于工业应用中的气体流量计至关重要。
飞行时间差法是超声波气体流量计中常用的直接测量方法,其原理是通过测量超声波在一定距离内的传播时间来计算上游和下游之间的飞行时间差。超声波的传播速度会叠加一个流体流速在超声波传播方向的分量,顺逆流超声波之间存在一个速度差,导致它们的转播时间存在差异。通过互相关算法可以计算出超声波在流体中传播的顺流和逆流飞行时间差,即可进一步计算得到流体的流速信息。因此,精确测量飞行时间差尤为重要。
目前,常用的飞行时间差测量方法主要有两种:阈值检测法和互相关法。由于超声波在气体中传播时衰减比较严重,因此接收到的超声波信号非常微弱,信噪比低,并且会受温度、压力等环境因素影响产生幅值波动和波形畸变,这些干扰会对阈值检测法会带来比较大的误差。在互相关法中,首先,顺流和逆流回波信号由高速ADC采样。然后计算两个回波信号的互相关函数。飞行时间差取决于与互相关函数的最大值对应的横坐标。与阈值检测法相比,互相关法仅与顺逆流回波信号的相关性有关,不易受噪声干扰,在复杂环境中具有更高的精度与稳定性。
发明内容
本发明针对现有超声流量计的不足以及实际应用的需要,提出一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,意在用一种有效的方法提高互相关法计算飞行时间差的分辨力和测量精度。
本发明一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,该方法的包括以下步骤:
S1,在管道中有气体流动时,用超声波流量计分别采集顺流超声波回波信号和逆流超声波回波信号;
S2,对采集到的顺流超声波回波信号和逆流超声波回波信号进行互相关计算;
S3,将计算得到的离散互相关函数的峰值附近进行插值处理;
进一步地,步骤S1所述的超声波回波信号采集中仅使用集成在低功耗单片机内部的数字模拟转换器。
进一步地,飞行时间差的确定方法具体为:
x(k)是逆流超声波回波信号的第k个采样点,y(k+τ)是顺流超声波回波信号的第k个采样点;Rxy(τ)是互相关函数计算结果;τ是互相关函数计算结果中幅值最大点对应的横坐标。
进一步地,步骤S3所述的离散互相关函数的峰值附近的插值处理,具体为:取在步骤S2计算得到的离散互相关函数的幅值最大值点以及该最大值点的前后各一点,在这3点之间进行插值。
进一步地,步骤S3所述的插值处理,具体为三次样条插值。
进一步地,还包括步骤S4,将插值后得到的函数曲线最大值作为互相关函数的最大值。
进一步地,还包括进而得到对应的超声波飞行时间差。
进一步地,步骤S1满足超声流量计的低功耗设计需求。
进一步地,对τ进一步计算,得到当前流速下超声波的顺逆流飞行时间差。
进一步地,取插值后的最大值点作为互相关函数的估计最大值点。
一种实现所述的用于飞行时间差计算的互相关插值方法的高精度低功耗气体超声波流量计,其特征在于气体超声波流量计采用单个低功耗单片机作为核心,由内部集成的ADC负责采集超声波回波信号,MCU负责互相关计算得到飞行时间差并转换为气体流量,MCU通过GPIO口负责流量计切换电路的控制、换能器激励信号发射、回波信号增益控制、通信以及显示屏的显示。
本发明具有的有益效果是:提出了一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,在满足低功耗设计的条件下,提高了互相关计算飞行时间差的分辨力和测量精度,进而提高了超声波流量计的测量精度。
附图说明
图1是本发明一个较佳实施方法的流程图;
图2是本发明一个较佳实施例中的回波信号;
图3是本发明一个较佳实施例中的计算得到的互相关函数;
图4是本发明一个较佳实施例中的对互相关函数得插值;
图5是本发明一个较佳实施例中的气体超声波流量计的整体结构图。
具体实施方式
以下通过具体实施例对本发明做进一步详细说明。本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本发明一个较佳实施例中,用于低功耗超声流量计飞行时间差计算的互相关插值方法流程大致如图1所示,具体步骤如下:
1)在管道中有气体流动时,用超声波流量计分别采集顺流超声波回波信号和逆流超声波回波信号,顺流超声波回波信号采样如图2所示,逆流超声波回波信号和图2类似。
2)对采集到的顺流超声波回波信号和逆流超声波回波信号进行互相关计算,计算结果如图3所示。
3)由于单片机ADC的采样频率还不满足高精度飞行时间差的测量分辨力要求,因此,需要在计算得到的离散互相关函数的峰值附近进行插值处理。取在步骤2)计算得到的离散互相关函数的幅值最大值点以及该最大值点的前后各一点,在这3点之间进行插值,取插值后的最大值点作为互相关函数的估计最大值点,如图4所示。本实施例中采用的插值方式为三次样条插值。
4)将插值后得到的函数曲线最大值作为互相关函数的预计最大值,进而得到其对应的横坐标,然后计算得到对应的超声波飞行时间差。
5)进一步地,飞行时间差的确定方法为:
x(k)是逆流超声波回波信号的第k个采样点,y(k+τ)是顺流超声波回波信号的第k个采样点;Rxy(τ)是互相关函数计算结果;τ是互相关函数计算结果中幅值最大点对应的横坐标;对τ进一步计算,得到当前流速下超声波的顺逆流飞行时间差。
图5是本发明一个较佳实施例中的气体超声波流量计的整体结构图。气体超声波流量计采用单个低功耗单片机作为核心,由内部集成的可编程脉冲发生器产生激励信号,激励信号经升压模块升压后激励发射换能器发出超声波信号,接收换能器接收到的回波信号经过滤波、放大后由单片机内部集成的ADC采集超声波回波信号,MCU负责互相关计算得到飞行时间差并转换为气体流量,MCU通过GPIO口负责流量计声道切换电路的控制、换能器激励信号发射、回波信号增益控制、通信、显示屏的显示以及其它功能。
以上详细描述了本发明的较佳具体实施例。应当理解,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化,形成本领域技术人员可以理解的其他实施方式。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (10)
1.一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,包括以下步骤:
S1,在管道中有气体流动时,用超声波流量计分别采集顺流超声波回波信号和逆流超声波回波信号;
S2,对采集到的顺流超声波回波信号和逆流超声波回波信号进行互相关计算;
S3,将计算得到的离散互相关函数的峰值附近进行插值处理。
2.根据权利要求1所述的一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,步骤S1所述的超声波回波信号采集中仅使用集成在低功耗单片机内部的数字模拟转换器。
4.根据权利要求1所述的一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,步骤S3所述的离散互相关函数的峰值附近的插值处理,具体为:取在步骤S2计算得到的离散互相关函数的幅值最大值点以及该最大值点的前后各一点,在这3点之间进行插值。
5.根据权利要求1所述的一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,步骤S3所述的插值处理,具体为三次样条插值。
6.根据权利要求1所述的一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,还包括步骤S4,将插值后得到的函数曲线最大值作为互相关函数的最大值。
7.根据权利要求6所述的一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,还包括进而得到对应的超声波飞行时间差。
8.根据权利要求2所述的一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,步骤S1满足超声流量计的低功耗设计需求。
9.根据权利要求3,所述的一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,对τ进一步计算,得到当前流速下超声波的顺逆流飞行时间差。
10.根据权利要求4所述的一种用于低功耗超声流量计飞行时间差计算的互相关插值方法,其特征在于,取插值后的最大值点作为互相关函数的估计最大值点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910981870.4A CN110646042A (zh) | 2019-10-16 | 2019-10-16 | 一种用于低功耗超声流量计飞行时间差计算的互相关插值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910981870.4A CN110646042A (zh) | 2019-10-16 | 2019-10-16 | 一种用于低功耗超声流量计飞行时间差计算的互相关插值方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110646042A true CN110646042A (zh) | 2020-01-03 |
Family
ID=69012831
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910981870.4A Pending CN110646042A (zh) | 2019-10-16 | 2019-10-16 | 一种用于低功耗超声流量计飞行时间差计算的互相关插值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110646042A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111220222A (zh) * | 2020-02-27 | 2020-06-02 | 成都千嘉科技有限公司 | 一种超声波燃气表流量的测定算法 |
CN112414481A (zh) * | 2020-10-16 | 2021-02-26 | 成都皓超科技有限公司 | 超声波燃气表超声波模块优化方法及其结构 |
CN113124948A (zh) * | 2021-05-20 | 2021-07-16 | 中国计量大学 | 一种基于fpga与互相关法的高精度时差测量方法 |
CN114152298A (zh) * | 2021-11-26 | 2022-03-08 | 浙江威星智能仪表股份有限公司 | 一种用于超声波计量飞行时间采样的方法 |
CN114397475A (zh) * | 2022-03-25 | 2022-04-26 | 青岛鼎信通讯股份有限公司 | 一种适用于超声波水表的水流流速测量方法 |
CN116878599A (zh) * | 2023-09-06 | 2023-10-13 | 青岛鼎信通讯科技有限公司 | 一种超声水表的流量计量方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0797105A3 (de) * | 1996-03-21 | 1999-09-22 | Siemens Aktiengesellschaft | Verfahren zur Laufzeitmessung eines elektrischen, elektromagnetischen oder akustischen Signals |
CN103542901A (zh) * | 2012-07-09 | 2014-01-29 | 德克萨斯仪器股份有限公司 | 一种流量计 |
DE102013106544A1 (de) * | 2013-06-24 | 2014-12-24 | Sick Ag | Ultraschallmessvorrichtung |
CN109029598A (zh) * | 2018-06-15 | 2018-12-18 | 浙江大学 | 一种互相关法气体超声波流量计中基于主次包络结合处的参考波形选取方法 |
CN109506727A (zh) * | 2018-12-24 | 2019-03-22 | 西安安森智能仪器股份有限公司 | 一种超声波流量测量方法及低功耗超声波流量计 |
US10309813B2 (en) * | 2015-05-15 | 2019-06-04 | Reliance Worldwide Corporation | Method and system for fluid flow rate measurement |
-
2019
- 2019-10-16 CN CN201910981870.4A patent/CN110646042A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0797105A3 (de) * | 1996-03-21 | 1999-09-22 | Siemens Aktiengesellschaft | Verfahren zur Laufzeitmessung eines elektrischen, elektromagnetischen oder akustischen Signals |
CN103542901A (zh) * | 2012-07-09 | 2014-01-29 | 德克萨斯仪器股份有限公司 | 一种流量计 |
DE102013106544A1 (de) * | 2013-06-24 | 2014-12-24 | Sick Ag | Ultraschallmessvorrichtung |
US10309813B2 (en) * | 2015-05-15 | 2019-06-04 | Reliance Worldwide Corporation | Method and system for fluid flow rate measurement |
CN109029598A (zh) * | 2018-06-15 | 2018-12-18 | 浙江大学 | 一种互相关法气体超声波流量计中基于主次包络结合处的参考波形选取方法 |
CN109506727A (zh) * | 2018-12-24 | 2019-03-22 | 西安安森智能仪器股份有限公司 | 一种超声波流量测量方法及低功耗超声波流量计 |
Non-Patent Citations (1)
Title |
---|
桂永芳 等: "基于互相关理论的超声波气体流量测量电路***", 《工业仪表与自动化装置》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111220222A (zh) * | 2020-02-27 | 2020-06-02 | 成都千嘉科技有限公司 | 一种超声波燃气表流量的测定算法 |
WO2021169883A1 (zh) * | 2020-02-27 | 2021-09-02 | 成都千嘉科技有限公司 | 一种超声波燃气表流量的测定算法 |
AU2021201596B2 (en) * | 2020-02-27 | 2022-04-21 | Chengdu Qianjia Technology Co. Ltd | Method for measuring flow rate of ultrasonic wave gas meter |
CN112414481A (zh) * | 2020-10-16 | 2021-02-26 | 成都皓超科技有限公司 | 超声波燃气表超声波模块优化方法及其结构 |
CN113124948A (zh) * | 2021-05-20 | 2021-07-16 | 中国计量大学 | 一种基于fpga与互相关法的高精度时差测量方法 |
CN114152298A (zh) * | 2021-11-26 | 2022-03-08 | 浙江威星智能仪表股份有限公司 | 一种用于超声波计量飞行时间采样的方法 |
CN114152298B (zh) * | 2021-11-26 | 2024-05-17 | 浙江威星智能仪表股份有限公司 | 一种用于超声波计量飞行时间采样的方法 |
CN114397475A (zh) * | 2022-03-25 | 2022-04-26 | 青岛鼎信通讯股份有限公司 | 一种适用于超声波水表的水流流速测量方法 |
CN116878599A (zh) * | 2023-09-06 | 2023-10-13 | 青岛鼎信通讯科技有限公司 | 一种超声水表的流量计量方法 |
CN116878599B (zh) * | 2023-09-06 | 2024-01-09 | 青岛鼎信通讯科技有限公司 | 一种超声水表的流量计量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110646042A (zh) | 一种用于低功耗超声流量计飞行时间差计算的互相关插值方法 | |
CN107290564B (zh) | 一种基于相位差的超声波流速测量方法 | |
CN109579950B (zh) | 一种气体超声波流量计的防错波检测装置 | |
CN109506727B (zh) | 一种超声波流量测量方法及低功耗超声波流量计 | |
CN109029598B (zh) | 一种互相关法气体超声波流量计中基于主次包络结合处的参考波形选取方法 | |
CN206930321U (zh) | 非满管超声波流量计 | |
CN107860430B (zh) | 一种基于时差法的超声波气体流量计时间差测量方法 | |
CN205037931U (zh) | 一种高精度超声波气体流量计 | |
CN104330120A (zh) | 用于低能耗超声波流量表的流量检测方法及*** | |
CN107727176A (zh) | 一种超声波气体流量计及其流量测量方法 | |
CN114397475A (zh) | 一种适用于超声波水表的水流流速测量方法 | |
CN105181046A (zh) | 一种基于双核心架构的数字化超声波流量计***及方法 | |
CN107655533A (zh) | 一种基于回波能量积分的气体超声波流量计信号处理方法和*** | |
CN110987102A (zh) | 一种高抗干扰高精度气体超声流量计及其测量方法 | |
CN102590598B (zh) | 基于多阈值比较的超声波信号过零点预测方法 | |
CN203037849U (zh) | 一种超声波测厚仪 | |
CN104897219A (zh) | 高精度低功耗超声流量计 | |
CN103869096B (zh) | 一种超声波风速仪量程拓宽方法 | |
CN102967334B (zh) | 利用对信号包络线处理测量流体流量的***及方法 | |
CN104596601A (zh) | 八声道超声波流量计传感器 | |
CN116027311A (zh) | 超声波测量回波波谷同步触发电路以及方法 | |
CN116046116A (zh) | 一种基于Lobe算法的超声波水表流量校准方法 | |
CN115808213A (zh) | 基于tdc方案的改进型气体计量方法 | |
CN212206222U (zh) | 一种高抗干扰高精度气体超声流量计 | |
CN204165597U (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200103 |
|
RJ01 | Rejection of invention patent application after publication |