CN115655423A - 一种超声波燃气表的校准方法 - Google Patents
一种超声波燃气表的校准方法 Download PDFInfo
- Publication number
- CN115655423A CN115655423A CN202211653357.0A CN202211653357A CN115655423A CN 115655423 A CN115655423 A CN 115655423A CN 202211653357 A CN202211653357 A CN 202211653357A CN 115655423 A CN115655423 A CN 115655423A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic gas
- gas meter
- instantaneous flow
- estimated
- formula
- 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
Images
Landscapes
- Measuring Volume Flow (AREA)
Abstract
本发明涉及一种超声波燃气表的校准方法,包括步骤:将所有超声波燃气表置于零流量下,获取平均个体时差误差;根据个体时差误差,获取瞬时流速公式;通过最小二乘法将瞬时流速公式转换为瞬时流量公式;计算误差权重函数;并基于实际流量获得目标函数,根据目标函数得到待估计参数;将待估计参数代入瞬时流量公式,通过瞬时流量公式判断各个超声波燃气表是否测试合格;若不合格,则针对不合格的超声波燃气表单独重复实施上述步骤,以获取该超声波燃气表的待估计参数。本方案可先消除各个超声波燃气表的个体误差,而后批量化校准超声波燃气表,在保证1.5级计量精度的前提下大幅提升生产效率。
Description
技术领域
本发明涉及超声波校准技术领域,特别涉及一种超声波燃气表的校准方法。
背景技术
在燃气计量领域中,目前使用最多的远传燃气表为膜式燃气表,但其存在结构复杂、易产生机械磨损、计量精度逐渐劣化等缺点。超声波燃气表作为一种纯数字式的计量仪器,因其工作性能稳定、远程抄表安全高效、计量精度可靠、易实现数字化,而成为一种必然的发展趋势。
然而,由于超声波燃气表中的计量模组存在个体差异,如超声换能器、气体流道、电子元器件等,导致每只超声波燃气表的表内计算流量与实际流量存在误差,需通过校验后方可投入时长。
发明内容
本发明的目的在于对批量化生产的超声波燃气表的流量测量精度进行校准,提供一种超声波燃气表的校准方法。
为了实现上述发明目的,本发明实施例提供了以下技术方案:
一种超声波燃气表的校准方法,包括以下步骤:
步骤1,将所有超声波燃气表置于零流量下,获取平均个体时差误差;根据个体时差误差,获取瞬时流速公式;通过最小二乘法将瞬时流速公式转换为瞬时流量公式;
步骤2,计算误差权重函数;并基于实际流量获得目标函数,根据目标函数得到待估计参数;
步骤3,将待估计参数代入瞬时流量公式,通过瞬时流量公式判断各个超声波燃气表是否测试合格;
步骤4,若合格,则结束;若不合格,则针对不合格的超声波燃气表单独重复实施步骤1~步骤3,以获取该超声波燃气表的待估计参数。
步骤1中所述将所有超声波燃气表置于零流量下,获取平均个体时差误差的步骤,包括:
测量超声波信号的上飞信号的飞行时间tu,和下飞信号的飞行时间td,计算上下飞时间差∆t=td-tu;
步骤1中所述根据个体时差误差,获取瞬时流速公式的步骤,包括:
获取上飞信号的飞行时间估算式:
获取下飞信号的飞行时间估算式:
其中,vz为瞬时流速。
步骤1中所述通过最小二乘法将瞬时流速公式转换为瞬时流量公式的步骤,包括:
其中,v为瞬时流量;a、b为待估计参数。
步骤2中所述计算误差权重函数的步骤,包括:
计算误差权重函数:
步骤2中所述并基于实际流量获得目标函数,根据目标函数得到待估计参数的步骤,包括:
其中,minL表示使目标函数L达到最小值;
目标函数L对待估计参数a、b的偏导数分别为:
通过上式,计算出待估计参数a、b。
步骤3中所述将待估计参数代入瞬时流量公式,通过瞬时流量公式判断各个超声波燃气表是否测试合格的步骤,包括:
与现有技术相比,本发明的有益效果:
本发明首先对每只超声波燃气表进行零流量校准,再使用加权的最小二乘法确定超声波燃气表的待估计参数a、b(也即校准系数),最后对测试不合格的超声波燃气表做单独的系数校验。本方案可先消除各个超声波燃气表的个体误差,而后批量化校准超声波燃气表,在保证1.5级计量精度的前提下大幅提升生产效率。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例校准方法流程图;
图2为本发明实施例流道中换能器设置示意图;
图3为本发明实施例流道中换能器另一种设置示意图;
图4为本发明实施例误差权重函数图像示意图。
具体实施方式
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。同时,在本发明的描述中,术语“第一”、“第二”等仅用于区分描述,而不能理解为指示或暗示相对重要性,或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。另外,术语“相连”、“连接”等可以是元件之间直接相连,也可以是经由其他元件的间接相连。
本发明通过下述技术方案实现,如图1所示,一种超声波燃气表的校准方法,包括以下步骤:
步骤1,将所有超声波燃气表置于零流量下,获取平均个体时差误差;根据个体时差误差,获取瞬时流速公式;通过最小二乘法将瞬时流速公式转换为瞬时流量公式。
超声波燃气表的原理是依靠一对换能器采集声波信号,根据飞行时差法测量流体的估算瞬时流速vz0,而后可换算为估算瞬时流量v0。以反射式的流道为例,如图2所示的A和B是一对换能器,发射的超声波信号与流道的水平方向夹角为,流道高度为D。另一种直射式流道,如图3所示的A和B是一对换能器。
测量时,换能器A和换能器B分别交替发送和接受超声波信号,其中由换能器A发射,换能器B接收的超声波信号沿待测量气体的顺流方向传递,称为“上飞信号”,图2或图3中水平箭头方向为待测量气体的顺流方向;反之,由换能器B发射,换能器A接收的超声波信号沿待测量提起的逆流方向传递,称为“下飞信号”。记上飞信号的飞行时间为tu,下飞信号的飞行时间为td,而对应的飞行时间估算式为:
可推导出估算瞬时流速vz0为:
而后,将估算瞬时流速vz0换算为估算瞬时流量v0时,需考虑流道横截面积的值,由于超声波燃气表受到如模具光洁面、换能器角度之类的个体差异影响,实际的瞬时流量v0应记为:
其中,k为待估计参数,且个体之间存在差异。
为保证计量精度同时提高生产效率,本方案首先进行零流量校准,将超声波燃气表置于标准测试温度下并保证处于零流量进行校准,校准过程持续5min以上,在此期间以125ms/次的频率不断测量上下飞时间差∆t=td-tu,接着计算平均个体时差误差:,其中n为测量总次数。
其中,vz为瞬时流速。
为将瞬时流速vz精准映射到瞬时流量v,通过加权的最小二乘法拟合校准系数,为进一步缩小超声波燃气表个体差异造成的影响,对式(2)进行改变:
其中,v为瞬时流量;a、b为待估计参数。
步骤2,计算误差权重函数;并基于实际流量获得目标函数,根据目标函数得到待估计参数。
其中,minL表示使目标函数L达到最小值;
目标函数L对待估计参数a、b的偏导数分别为:
通过上式,计算出待估计参数a、b。
步骤3,将待估计参数代入瞬时流量公式,通过瞬时流量公式判断各个超声波燃气表是否测试合格。
步骤4,若合格,则结束;若不合格,则针对不合格的超声波燃气表单独重复实施步骤1~步骤3,以获取该超声波燃气表的待估计参数。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (7)
1.一种超声波燃气表的校准方法,其特征在于:包括以下步骤:
步骤1,将所有超声波燃气表置于零流量下,获取平均个体时差误差;根据个体时差误差,获取瞬时流速公式;通过最小二乘法将瞬时流速公式转换为瞬时流量公式;
步骤2,计算误差权重函数;并基于实际流量获得目标函数,根据目标函数得到待估计参数;
步骤3,将待估计参数代入瞬时流量公式,通过瞬时流量公式判断各个超声波燃气表是否测试合格;
步骤4,若合格,则结束;若不合格,则针对不合格的超声波燃气表单独重复实施步骤1~步骤3,以获取该超声波燃气表的待估计参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211653357.0A CN115655423B (zh) | 2022-12-22 | 2022-12-22 | 一种超声波燃气表的校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211653357.0A CN115655423B (zh) | 2022-12-22 | 2022-12-22 | 一种超声波燃气表的校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115655423A true CN115655423A (zh) | 2023-01-31 |
CN115655423B CN115655423B (zh) | 2023-03-21 |
Family
ID=85022397
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211653357.0A Active CN115655423B (zh) | 2022-12-22 | 2022-12-22 | 一种超声波燃气表的校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115655423B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117029975A (zh) * | 2023-10-09 | 2023-11-10 | 成都千嘉科技股份有限公司 | 一种超声波燃气表的多环境自适应方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103074873A (zh) * | 2013-02-07 | 2013-05-01 | 天津大学 | 一种非理想流场下的多声道超声流量计声道布置方法 |
US20140083202A1 (en) * | 2011-05-17 | 2014-03-27 | Endress + Hauser Flowtec Ag | Ultrasonic, Flow Measuring Device |
CN103913211A (zh) * | 2014-04-14 | 2014-07-09 | 姜跃炜 | 超声波水表的时间系数标定方法 |
US20150013471A1 (en) * | 2012-03-16 | 2015-01-15 | Hitachi Aloka Medical, Ltd. | Fluid flow rate detection device |
CN105222839A (zh) * | 2015-08-21 | 2016-01-06 | 宁波市鹏盛科技发展有限公司 | 超声波流量计仪表系数的非线性修正方法 |
CN110646041A (zh) * | 2019-09-29 | 2020-01-03 | 宁夏隆基宁光仪表股份有限公司 | 一种多声道超声流量计的数据融合计算方法 |
CN111323101A (zh) * | 2020-03-24 | 2020-06-23 | 成都千嘉科技有限公司 | 超声波表自适应的自动标定方法 |
CN112964337A (zh) * | 2021-03-17 | 2021-06-15 | 成都千嘉科技有限公司 | 一种应用于超声波燃气表的校准系统及方法 |
CN113029263A (zh) * | 2021-03-24 | 2021-06-25 | 成都千嘉科技有限公司 | 基于样本方差剔除误差的燃气损失流量测定方法及系统 |
CN113624305A (zh) * | 2021-08-27 | 2021-11-09 | 成都千嘉科技有限公司 | 一种超声波流量计校准方法及系统 |
CN114001804A (zh) * | 2021-11-01 | 2022-02-01 | 成都千嘉科技有限公司 | 一种基于时间差法超声波计量装置的校准方法及系统 |
CN114935388A (zh) * | 2022-07-26 | 2022-08-23 | 成都千嘉科技股份有限公司 | 超声波燃气表的校正系统及方法 |
-
2022
- 2022-12-22 CN CN202211653357.0A patent/CN115655423B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140083202A1 (en) * | 2011-05-17 | 2014-03-27 | Endress + Hauser Flowtec Ag | Ultrasonic, Flow Measuring Device |
US20150013471A1 (en) * | 2012-03-16 | 2015-01-15 | Hitachi Aloka Medical, Ltd. | Fluid flow rate detection device |
EP2826425A1 (en) * | 2012-03-16 | 2015-01-21 | The University of Tokyo | Device for detecting fluid flow rate |
CN103074873A (zh) * | 2013-02-07 | 2013-05-01 | 天津大学 | 一种非理想流场下的多声道超声流量计声道布置方法 |
CN103913211A (zh) * | 2014-04-14 | 2014-07-09 | 姜跃炜 | 超声波水表的时间系数标定方法 |
CN105222839A (zh) * | 2015-08-21 | 2016-01-06 | 宁波市鹏盛科技发展有限公司 | 超声波流量计仪表系数的非线性修正方法 |
CN110646041A (zh) * | 2019-09-29 | 2020-01-03 | 宁夏隆基宁光仪表股份有限公司 | 一种多声道超声流量计的数据融合计算方法 |
CN111323101A (zh) * | 2020-03-24 | 2020-06-23 | 成都千嘉科技有限公司 | 超声波表自适应的自动标定方法 |
CN112964337A (zh) * | 2021-03-17 | 2021-06-15 | 成都千嘉科技有限公司 | 一种应用于超声波燃气表的校准系统及方法 |
CN113029263A (zh) * | 2021-03-24 | 2021-06-25 | 成都千嘉科技有限公司 | 基于样本方差剔除误差的燃气损失流量测定方法及系统 |
CN113624305A (zh) * | 2021-08-27 | 2021-11-09 | 成都千嘉科技有限公司 | 一种超声波流量计校准方法及系统 |
CN114001804A (zh) * | 2021-11-01 | 2022-02-01 | 成都千嘉科技有限公司 | 一种基于时间差法超声波计量装置的校准方法及系统 |
CN114935388A (zh) * | 2022-07-26 | 2022-08-23 | 成都千嘉科技股份有限公司 | 超声波燃气表的校正系统及方法 |
Non-Patent Citations (1)
Title |
---|
郑安芳 等: ""改进的BP算法在多声道超声波流量计中的应用"" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117029975A (zh) * | 2023-10-09 | 2023-11-10 | 成都千嘉科技股份有限公司 | 一种超声波燃气表的多环境自适应方法 |
CN117029975B (zh) * | 2023-10-09 | 2023-12-15 | 成都千嘉科技股份有限公司 | 一种超声波燃气表的多环境自适应方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115655423B (zh) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105841762B (zh) | 超声波水表的流量计量方法和系统 | |
CN115655423B (zh) | 一种超声波燃气表的校准方法 | |
CN107655552B (zh) | 一种基于非测量段传播时差及零漂误差的超声水表流速修正方法 | |
CN114964429A (zh) | 一种超声波水表检定方法 | |
US5827963A (en) | System and method for determining a density of a fluid | |
CN105403265A (zh) | 一种自动校正零点漂移的超声水表及其校正方法 | |
CN114001804B (zh) | 一种基于时间差法超声波计量装置的校准方法及系统 | |
CN114777887A (zh) | 一种超声波水表的零点漂移消除方法 | |
CN111238587A (zh) | 一种超声波换能器传输时间校准方法 | |
CN109612607B (zh) | 一种温度传感器反应速度测试方法 | |
CN117405187A (zh) | 一种四声道气体超声波流量计的声道异常筛选方法 | |
CN205246150U (zh) | 一种自动校正零点漂移的超声水表 | |
CN114580214B (zh) | 基于物联网的智能型蒸汽流量计校核测试系统及方法 | |
CN110081943B (zh) | 一种科氏力质量流量计温度补偿的方法 | |
CN114812707A (zh) | 一种超声波计量表计量方法 | |
CN114485863A (zh) | 超声波水表的流量误差校正方法、系统、计算机及介质 | |
CN107796492B (zh) | 一种超声波燃气表在线式校准方法 | |
CN117309075A (zh) | 基于矫正后飞行时间的超声波燃气表检测修正系统和方法 | |
CN111561888A (zh) | 一种超声波测厚值的修正方法 | |
CN117029975B (zh) | 一种超声波燃气表的多环境自适应方法 | |
CN113569202B (zh) | 统一基准订正的河流流量全量程测量不确定度计算方法 | |
CN117407636B (zh) | 一种水下成像折射率动态修正方法 | |
CN116858346B (zh) | 基于超声波流量计的校准方法和校准装置 | |
CN116147741B (zh) | 一种基于NB-IoT的超声波水表流量校准方法及系统 | |
CN114235106A (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 |