CN106706056A - 一种大口径超声波水表流量测量的补偿方法 - Google Patents

一种大口径超声波水表流量测量的补偿方法 Download PDF

Info

Publication number
CN106706056A
CN106706056A CN201710131154.8A CN201710131154A CN106706056A CN 106706056 A CN106706056 A CN 106706056A CN 201710131154 A CN201710131154 A CN 201710131154A CN 106706056 A CN106706056 A CN 106706056A
Authority
CN
China
Prior art keywords
ultrasonic
time
propagation time
temperature
accurate
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
CN201710131154.8A
Other languages
English (en)
Other versions
CN106706056B (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.)
JINAN RUIQUAN ELECTRONICS CO Ltd
Original Assignee
JINAN RUIQUAN ELECTRONICS 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 JINAN RUIQUAN ELECTRONICS CO Ltd filed Critical JINAN RUIQUAN ELECTRONICS CO Ltd
Priority to CN201710131154.8A priority Critical patent/CN106706056B/zh
Publication of CN106706056A publication Critical patent/CN106706056A/zh
Application granted granted Critical
Publication of CN106706056B publication Critical patent/CN106706056B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F15/00Details of, or accessories for, apparatus of groups G01F1/00 - G01F13/00 insofar as such details or appliances are not adapted to particular types of such apparatus
    • G01F15/02Compensating or correcting for variations in pressure, density or temperature

Landscapes

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

Abstract

本发明涉及自来水设施领域,特别公开了一种大口径超声波水表流量测量的补偿方法。它包括如下步骤:1)通过时间测量模块的第一波检测功能检测超声波在基表管段内的传播时间,通过换算得到当前流体的温度;2)通过大量实验数据拟合温度和超声波声波传输的时间的关系公式;3)利用关系公式来补偿TDC_GP22寄存器2中的屏蔽时间,来优化采样时机,获得准确的传播时间;4)通过错波滤波算法,对传播时间进行滤波,最终得到稳定准确的超声波传播时间;5)通过插值算法,积分得到准确的累计流量。通过此方法可以在任意一个温度下,只需要一次合格的校正,就可以实现整个温度要求范围内的准确测量。

Description

一种大口径超声波水表流量测量的补偿方法
(一)技术领域
本发明涉及自来水设施领域,特别涉及一种大口径超声波水表流量测量的补偿方法。
(二)背景技术
目前技术下,大口径超声波水表的流量补偿方法是利用PT100或者热敏电阻测到的温度,然后再对流量进行补偿,最终通过简单的滤波方法来积分流量。这种方法在理想状态下是可以实现精确的流量计算,但是由于受限于加工精度,主要的测量模块(流量传感器和超声波基表)无法做到硬件参数的完全一致,这样就导致仅通过滤波算法无法在整个温度要求范围内得到精确的流量测量。
(三)发明内容
本发明为了弥补现有技术的缺陷,提供了一种宽范围的、低加工精度要求、自动温度补偿、高精度测量的大口径超声波水表流量测量的补偿方法。
本发明是通过如下技术方案实现的:
一种大口径超声波水表流量测量的补偿方法,其特征是:包括如下步骤:
1)通过时间测量模块的第一波检测功能检测超声波在基表管段内的传播时间,通过换算得到当前流体的温度;
2)通过大量实验数据拟合温度和超声波声波传输的时间的关系公式;
3)利用关系公式来补偿TDC_GP22寄存器2中的屏蔽时间,来优化采样时机,获得准确的传播时间;
4)通过错波滤波算法,对传播时间进行滤波,最终得到稳定准确的超声波传播时间;
5)通过插值算法,积分得到准确的累计流量。
其中,所述时间测量模块为TDC_GP22。
具体地,包括以下步骤:
1)根据TDC_GP22的第一波检测功能检测超声波在基表管段内的传播时间 t和t,即为超声波从超声波换能器a到超声波换能器b的飞行时间和从超声波换能器b到超声波换能器a的飞行时间;
2)通过公式T=C0²*Δt/(2L*k)计算得到流体的温度,其中,
C0为超声波在水中的传播速度,
k为温度修正经验系数,
Δt=t-t=2L*cosθ/(C0²-u²*cosθ),u为流体的流动速度,θ为超声波换能器a和超声波换能器b之间之间的连线与流体的流动方向之间的夹角,
L为超声波换能器a和超声波换能器b之间的中心距离;
3)根据流体温度的变化来动态计算屏蔽时间,得到T和t之间的曲线;计算公式:
t=(74.606T²-8828.3T+(3E+06)*K
式中:t为屏蔽时间,T为流体温度,K为修正系数
通过精确的屏蔽时间,来优化采样时机,最终得到精确的传播时间。
4)采用错波滤波算法过滤掉错误的回波,得到稳定的超声波传播时间;
5)采用插值法对超声波水表进行校正:通过预先设定好的数据表,分别对流体温度和流量点进行插值计算,对累计流量机型补偿。
进一步地,在步骤3)中,用示波器来测量在不同流体温度下,超声波在流体中的传播时间。
本发明的有益效果是:
本发明可以在任意一个温度下,只需要一次合格的校正,就可以实现整个温度要求范围内的准确测量,方法简单,方便、快捷。
(四)附图说明
下面结合附图对本发明作进一步的说明。
图1为本发明的补偿过程示意图。
(五)具体实施方式
附图为本发明的具体实施例。如图1所示,该种大口径超声波水表流量测量的补偿方法,包括如下步骤:
如图1所示,a,b为超声波换能器;L为a和b之间的中心距离,θ为a,b之间的连线与流体流动方向之间的夹角;
1)首先根据TDC_GP22的第一波检测功能检测超声波在基表管段内的传播时间t和t,即为超声波从a到b的飞行时间和从b到a的飞行时间。
2)当流体以速度u流动时,超声波的实际传播速度C是声速C0和流体在声道方向上的速度分量的叠加,即:
C=C0+u*cosθ (2.1)
超声波在流体中顺流和逆流的传播时间分别为:
t=L/(C0+u*cosθ) (2.2)
t=L/(C0-u*cosθ) (2.3)
式中:C0为超声波在静止流体中的传播速度。
所以:t=t-t=2L*cosθ/(C0²-u²*cosθ) (2.4)
一般情况下,在多数工业系统中流体流速远小于声波传播的速度,即u²<<C0²,所以:
Δt=t-t=2L*u*cosθ/C0² (2.5)
即:u=C0²*Δt/(2L*cosθ) (2.6)
代入公式(2.1)可得到公式:
C=C0+C0²*Δt/2L (2.7)
根据在不同水温下超声波传播时间的不同得到经验公式
C=C0+k*T (2.8)
式中C0为超声波在水中的传播速度,k为温度修正经验系数,T为流体温度,
代入公式(2.8)中可以得到公式:
T=C0²*Δt/(2L*k) (2.9)
3)根据流体温度T的变化来动态计算屏蔽时间t(用来设置TDC-GP22中的Reg2寄存器),得到T和t之间的曲线公式。具体的做法是,用示波器来测量在不同温度下,超声波在流体中的传播时间,
通过大量的实验数据,得到以下公式:
t=(74.606T²-8828.3T+(3E+06)*K
式中:t为屏蔽时间,T为流体温度,K为修正系数
通过精确的屏蔽时间,来优化采样时机,最终得到精确的传播时间。
4)在长期使用过程中,由于杂质的沉积,反射面反射的信号会散射并减弱,超声波换能器接收到的信号强度也会变弱,这样在测量回波时,很容易出现错波的现象,我们采用错波滤波算法过滤掉错误的回波,得到稳定的超声波传播时间。
错波滤波算法:从理论分析可以知道,在流动的流体中,t>t,以此为判断依据对采样数据进行滤波。
通过对顺流时间和逆流时间的测量,得到6组时间数据,顺流时间(t顺1,t顺2,t顺3)和逆流时间(t逆1,t逆2,t逆3)各3组,t0为经验系数。
首先对前两组进行滤波计算:即t顺1,t顺2,t逆1,t逆2
a)首先验证t顺1和t逆1是同一组测量出来的顺流和逆流时间,即要符合 t顺1≤t逆1并且t顺1+t0>t逆1,则时差为t逆1-t顺1
b)不满足a的条件,则说明出现了错波现象,此时要用t顺2和t逆1进行匹配。即符合t顺2≤t逆1并且t顺2+t0>t逆1,则时差为t逆1-t顺2
c)不满足b的条件,此时用t顺1和t逆2进行匹配,即符合t顺1≤t逆2并且t顺1+t0>t逆2,则时差为t逆2-t顺1
d)以上条件均不符合,则时差默认设定为无意义。
然后再对后两组进行滤波计算:即t顺2、t逆2、t顺3和t逆3
e)首先验证t顺2和t逆2是同一组测量出来的顺流和逆流时间,即要符合 t顺2≤t逆2并且t顺2+t0>t逆2,则时差为t逆2-t顺2
f)不满足e的条件,则说明出现了错波现象,此时要用t顺3和t逆2进行匹配。即符合t顺3≤t逆2并且t顺3+t0>t逆2,则时差为t逆2-t顺3
g)不满足f的条件,此时用t顺2和t逆3进行匹配,即符合t顺2≤t逆3并且t顺2+t0>t逆3,则时差为t逆3-t顺2
h)以上条件均不符合,则时差默认设定为无意义。
通过这种方法就可以得到两组稳定时差,再换算成流量。
5)在对超声波水表进行校正时,不可能对每一个温度点和流量点都进行校正,此时采用插值法,通过预先设定好的数据表,分别对温度和流量点进行插值计算,对累计流量机型补偿。通过此方法可以在任意一个温度下,只需要一次合格的校正,就可以实现整个温度要求范围内的准确测量。
除说明书所述技术特征外,其余技术特征均为本领域技术人员已知技术。

Claims (4)

1.一种大口径超声波水表流量测量的补偿方法,其特征是:包括如下步骤:
1)通过时间测量模块的第一波检测功能检测超声波在基表管段内的传播时间,通过换算得到当前流体的温度;
2)通过大量实验数据拟合温度和超声波声波传输的时间的关系公式;
3)利用关系公式来补偿TDC_GP22寄存器2中的屏蔽时间,来优化采样时机,获得准确的传播时间;
4)通过错波滤波算法,对传播时间进行滤波,最终得到稳定准确的超声波传播时间;
5)通过插值算法,积分得到准确的累计流量。
2.根据权利要求1所述的大口径超声波水表流量测量的补偿方法,其特征是,所述时间测量模块为TDC_GP22。
3.根据权利要求2所述的大口径超声波水表流量测量的补偿方法,其特征是,包括以下步骤:
1)根据TDC_GP22的第一波检测功能检测超声波在基表管段内的传播时间 t和t,即为超声波从超声波换能器a到超声波换能器b的飞行时间和从超声波换能器b到超声波换能器a的飞行时间;
2)通过公式T=C0²*Δt/(2L*k)计算得到流体的温度,其中,
C0为超声波在水中的传播速度,
k为温度修正经验系数,
Δt=t-t=2L*cosθ/(C0²-u²*cosθ),u为流体的流动速度,θ为超声波换能器a和超声波换能器b之间之间的连线与流体的流动方向之间的夹角,
L为超声波换能器a和超声波换能器b之间的中心距离;
3)根据流体温度的变化来动态计算屏蔽时间,得到T和t之间的曲线;计算公式:
t=(74.606T²-8828.3T+(3E+06)*K
式中:t为屏蔽时间,T为流体温度,K为修正系数;
通过精确的屏蔽时间,来优化采样时机,最终得到精确的传播时间;
4)采用错波滤波算法过滤掉错误的回波,得到稳定的超声波传播时间;
5)采用插值法对超声波水表进行校正:通过预先设定好的数据表,分别对流体温度和流量点进行插值计算,对累计流量机型补偿。
4.根据权利要求3所述的大口径超声波水表流量测量的补偿方法,其特征是,在步骤3)中,用示波器来测量在不同温度下,超声波在流体中的传播时间。
CN201710131154.8A 2017-03-07 2017-03-07 一种大口径超声波水表流量测量的补偿方法 Active CN106706056B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710131154.8A CN106706056B (zh) 2017-03-07 2017-03-07 一种大口径超声波水表流量测量的补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710131154.8A CN106706056B (zh) 2017-03-07 2017-03-07 一种大口径超声波水表流量测量的补偿方法

Publications (2)

Publication Number Publication Date
CN106706056A true CN106706056A (zh) 2017-05-24
CN106706056B CN106706056B (zh) 2019-07-26

Family

ID=58918144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710131154.8A Active CN106706056B (zh) 2017-03-07 2017-03-07 一种大口径超声波水表流量测量的补偿方法

Country Status (1)

Country Link
CN (1) CN106706056B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107607168A (zh) * 2017-07-25 2018-01-19 辽宁航宇星物联仪表科技有限公司 一种户用超声波水表自适应方法
CN107621292A (zh) * 2017-07-25 2018-01-23 辽宁航宇星物联仪表科技有限公司 一种户用超声波水表错波补偿方法
CN109782024A (zh) * 2018-12-29 2019-05-21 杭州先锋电子技术股份有限公司 基于定态参数进行波数判断的超声波气表防计量错误方法
CN112179428A (zh) * 2020-09-15 2021-01-05 济南瑞泉电子有限公司 一种基于线-相关算法的超声波水表测量模型
CN113899418A (zh) * 2021-12-09 2022-01-07 西安安森智能仪器股份有限公司 一种超声波换能器完全不嵌入管道内部安装流量计算方法
CN114543949A (zh) * 2022-02-24 2022-05-27 安徽汉威电子有限公司 一种超声波水表自适应计量方法
CN117470405A (zh) * 2023-12-27 2024-01-30 青岛积成电子股份有限公司 基于超声波飞行时间测量流体温度的方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0791996A (ja) * 1993-09-22 1995-04-07 Mitsubishi Heavy Ind Ltd 超音波流量計
US5831175A (en) * 1996-06-12 1998-11-03 Welch Allyn, Inc. Method and apparatus for correcting temperature variations in ultrasonic flowmeters
CN1720428A (zh) * 2002-12-06 2006-01-11 恩德斯+豪斯流量技术股份有限公司 过程测量仪表
CN102147275A (zh) * 2011-03-18 2011-08-10 中国计量学院 低功耗高精度超声波流量测量方法及装置
CN102713531A (zh) * 2010-01-07 2012-10-03 松下电器产业株式会社 超声波流量计
WO2014039873A1 (en) * 2012-09-07 2014-03-13 Daniel Measurement And Control, Inc. Ultrasonic flow metering using compensated computed temperature
CN104390680A (zh) * 2014-12-22 2015-03-04 湖南威铭能源科技有限公司 时差式超声波热量表流量误差修正方法
CN105091970A (zh) * 2015-08-07 2015-11-25 北京控制工程研究所 超声波流量计动态补偿方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0791996A (ja) * 1993-09-22 1995-04-07 Mitsubishi Heavy Ind Ltd 超音波流量計
US5831175A (en) * 1996-06-12 1998-11-03 Welch Allyn, Inc. Method and apparatus for correcting temperature variations in ultrasonic flowmeters
CN1720428A (zh) * 2002-12-06 2006-01-11 恩德斯+豪斯流量技术股份有限公司 过程测量仪表
CN102713531A (zh) * 2010-01-07 2012-10-03 松下电器产业株式会社 超声波流量计
CN102147275A (zh) * 2011-03-18 2011-08-10 中国计量学院 低功耗高精度超声波流量测量方法及装置
WO2014039873A1 (en) * 2012-09-07 2014-03-13 Daniel Measurement And Control, Inc. Ultrasonic flow metering using compensated computed temperature
CN104390680A (zh) * 2014-12-22 2015-03-04 湖南威铭能源科技有限公司 时差式超声波热量表流量误差修正方法
CN105091970A (zh) * 2015-08-07 2015-11-25 北京控制工程研究所 超声波流量计动态补偿方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107607168A (zh) * 2017-07-25 2018-01-19 辽宁航宇星物联仪表科技有限公司 一种户用超声波水表自适应方法
CN107621292A (zh) * 2017-07-25 2018-01-23 辽宁航宇星物联仪表科技有限公司 一种户用超声波水表错波补偿方法
CN109782024A (zh) * 2018-12-29 2019-05-21 杭州先锋电子技术股份有限公司 基于定态参数进行波数判断的超声波气表防计量错误方法
CN112179428A (zh) * 2020-09-15 2021-01-05 济南瑞泉电子有限公司 一种基于线-相关算法的超声波水表测量模型
CN113899418A (zh) * 2021-12-09 2022-01-07 西安安森智能仪器股份有限公司 一种超声波换能器完全不嵌入管道内部安装流量计算方法
CN114543949A (zh) * 2022-02-24 2022-05-27 安徽汉威电子有限公司 一种超声波水表自适应计量方法
CN117470405A (zh) * 2023-12-27 2024-01-30 青岛积成电子股份有限公司 基于超声波飞行时间测量流体温度的方法

Also Published As

Publication number Publication date
CN106706056B (zh) 2019-07-26

Similar Documents

Publication Publication Date Title
CN106706056A (zh) 一种大口径超声波水表流量测量的补偿方法
CN106643939B (zh) 用于超声波流量计计算超声波传播时间的方法
CN103808381B (zh) 一种时差式超声波流量计的温度影响消除方法
CN108020282B (zh) 基于复系数滤波的科氏质量流量计信号处理方法
CN105841762B (zh) 超声波水表的流量计量方法和系统
CN107860430B (zh) 一种基于时差法的超声波气体流量计时间差测量方法
US11199429B2 (en) Method for measuring a speed of a fluid
CN104913830A (zh) 一种基于曲线拟合的雷达物位计信号处理方法
CN107478282A (zh) 超声流量检测信号处理方法及装置、时差法超声检测系统
CN203287060U (zh) 双声道超声波流量检测系统
US20220244082A1 (en) Flow measurement method based on dynamic optimization of three pressure sensors
CN103913203B (zh) 超声波水表流量系数处理方法
CN104748704A (zh) 薄壁结构超声共振测厚频谱分析内插校正方法
JP2014516162A (ja) 超音波伝搬時間法による流体の流量検出方法
KR101764870B1 (ko) 초음파 유량계의 신호처리시스템
Rychagov et al. Multipath flowrate measurements of symmetric and asymmetric flows
CN109324208A (zh) 一种基于超声波声速法的汽液两相流密度、质量流量及相含率一体化分析仪
CN114877975A (zh) 一种基于管路入口雷诺数的超声水表流量误差修正方法
Li et al. A novel differential time-of-flight algorithm for high-precision ultrasonic gas flow measurement
CN105372493B (zh) 基于三条dft复数谱线的信号幅值和相位测量方法
CN110081943B (zh) 一种科氏力质量流量计温度补偿的方法
CN107064546B (zh) 一种基于非圆信号频率估计算法的流体速度估计方法
CN107796492B (zh) 一种超声波燃气表在线式校准方法
CN110749357A (zh) 基于max35103和max31865测量系统的双通道液体超声波流量计测量方法
CN205300668U (zh) 双声道超声波时间检测装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP02 Change in the address of a patent holder

Address after: Room 2103 and 2104, 21 stories, Building 157 A, Shanda Road, Lixia District, Jinan City, Shandong Province

Patentee after: JINAN RICHNES ELECTRONIC Co.,Ltd.

Address before: Room 2103-2104, 21 storeys, Building 157 A, Shanda Road, Licheng District, Jinan City, Shandong Province

Patentee before: JINAN RICHNES ELECTRONIC Co.,Ltd.

CP02 Change in the address of a patent holder
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A compensation method for flow measurement of large diameter ultrasonic water meter

Effective date of registration: 20210629

Granted publication date: 20190726

Pledgee: Bank of Beijing Co.,Ltd. Jinan Branch

Pledgor: JINAN RICHNES ELECTRONIC Co.,Ltd.

Registration number: Y2021370000078

PC01 Cancellation of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20230621

Granted publication date: 20190726

Pledgee: Bank of Beijing Co.,Ltd. Jinan Branch

Pledgor: JINAN RICHNES ELECTRONIC Co.,Ltd.

Registration number: Y2021370000078

PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A Compensation Method for Flow Measurement of Large Caliber Ultrasonic Water Meters

Effective date of registration: 20230628

Granted publication date: 20190726

Pledgee: Bank of Beijing Co.,Ltd. Jinan Branch

Pledgor: JINAN RICHNES ELECTRONIC Co.,Ltd.

Registration number: Y2023980046322

PE01 Entry into force of the registration of the contract for pledge of patent right