1、第一作者简介:权亚强，男，高级工程师，本科，主要从事燃气计量技术标准制定研究以及智能燃气表开发设计工作。收稿日期:2022 03 01;修回日期:2022 08 20燃气表检定气路系统泄漏率的不确定度评定权亚强，侯雪丹(成都秦川物联网科技股份有限公司，四川 成都 610100)摘要:给出临界流文丘里喷嘴式燃气表检定装置检定气路系统组成，提出在燃气表示值相对误差试验中，检定气路系统泄漏率不确定度最大允许值的确定方法。基于检定气路系统泄漏率测量过程建立测量模型，通过对泄漏率的不确定度进行 B 类评定，分析泄漏率测量过程中引入的不确定度分量，对其进行评定合成。给出泄漏率的不确定度评定算例及检定装置的
2、密封性检测建议。关键词:检定气路系统;示值相对误差;泄漏率;不确定度中图分类号:TU996 7文献标志码:B文章编号:1000 4416(2023)02 0B30 041概述概述临界流文丘里喷嘴式燃气表检定装置一般由检定气路系统、信息采集系统、操作控制系统和显示系统组成1，其中检定气路系统见图 1。临界流文丘里喷嘴式燃气表检定装置有正压法和负压法之分。图 1 为负压法。滞止容器（装有一组喷嘴）大气进口串联被检表组的进口阀组串联被检表组串联被检表组的出口阀组各喷嘴支路阀组喷嘴支路汇入的背压容器真空泵大气出口标准表图 1检定气路系统检定气路系统由大气进口、串联被检表组的进口阀组、串联被检表组、串联
3、被检表组的出口阀组、滞止容器(装有一组并联的且可实现特定流量组合的临界流文丘里喷嘴)、各喷嘴支路阀组、喷嘴支路汇入的背压容器、真空泵和大气出口组成。2示值相对误差及检定气路系统泄漏分析示值相对误差及检定气路系统泄漏分析示值相对误差是指燃气表显示的气体体积(称为示值)和实际通过燃气表的气体体积(称实际值)之差与实际值的百分比。示值相对误差是燃气表最重要的计量特性之一，本质上反映了燃气表准确度。示值相对误差试验是指在相同时间间隔内，气体连续流过串联被检表组与滞止容器中的临界流文丘里喷嘴组成的串联管路系统。根据检定流量点选择开启临界流文丘里喷嘴，根据临界流文丘里喷嘴处的气体滞止温度、滞止压力等参数计
4、算出实际值，根据实际值和示值计算得到示值相对误差。若检定气路系统无任何泄漏，根据连续性方程，流过管路任意截面的质量流量处处相等，所以采用临界流文丘里喷嘴的喉部质量流量可以作为标准质量流量对被检燃气表进行检定和校准。若检定气路系统有泄漏时，临界流文丘里喷嘴的喉部质量流量与被检燃气表示值不相等。03B第 43 卷第 2 期2023 年 2 月煤 气 与 热 力GAS HEATVol 43 No 2Feb 2023DOI:10.13608/ki.1000-4416.2023.02.008由示值相对误差测量过程可知，其测量不确定度分量来源主要有:测量重复性引起的不确定度分量 urep、检定装置本身的不
5、确定度分量 ueqp、检定气路系统泄漏率的不确定度分量 uleak。GB/T 69682019膜式燃气表 规定在流量高区(分界流量 qt到最大流量 qmax的区间内，包含 qt)，示值相对误差绝对值最大允许值 IMPEV为 1 5%。在流量低区(最小流量 qmin到分界流量 qt的区间内)，IMPEV为 3%。由JJF 10942002测量仪器特性评定 第 5.3.1.4 款可知，示值相对误差的扩展不确定度 U(包含因子 k为 2)应小于或等于 IMPEV的 1/3。也就是说将 urep、ueqp、uleak合成得到的合成标准不确定度乘以包含因子 k 后，得到的值应不超过相应流量点 IMPEV
6、的1/32，即:ku2rep+u2eqp+u2leak13IMPEV(1)式中k 包含因子urep 测量重复性引起的不确定度分量ueqp 检定装置本身的不确定度分量uleak 检定气路系统泄漏率的不确定度分量IMPEV 示值相对误差绝对值最大允许值若检定流量点为流量高区，IMPEV为 1.5%。假设此时测量重复性引起的不确定度分量 urep为0.12%，检 定 装 置 本 身 的 不 确 定 度 分 量 ueqp为0.20%，则检定气路系统泄漏率的不确定度分量uleak应不大于 0.09%。3泄漏率的不确定度评定泄漏率的不确定度评定本文的泄漏率的不确定度评定包含泄漏率的 B类相对标准不确定度
7、ur(q)、泄漏率的相对扩展不确定度 Ur、泄漏率的扩展不确定度 U。3 1泄漏率测量模型在对临界流文丘里喷嘴式燃气表检定装置进行密封性检测时，其泄漏率测量模型如下1:q=3 600pVfbpjct(2)式中q 泄漏率，m3/hp 泄漏引起封闭系统压力变化量，PaVfb 封闭系统总净容积，m3pjc 封闭系统检测压力，Pat 检测时间，s由式(2)可以看出，泄漏率与封闭系统总净容积、检测压力、检测时间有关。相同结构的被测封闭系统压力变化量越大，泄漏率也越大。3 2泄漏率的 B 类相对标准不确定度ur(q)计算泄漏率的 B 类相对标准不 确 定 度 计 算 如下3 18:u2r()q=cr()p
8、 ur()p2+crV()fburV()fb2+crp()jcurp()jc2+cr()t ur()t2(3)式中ur(q)泄漏率的 B 类相对标准不确定度cr(p)、ur(p)泄漏引起封闭系统压力变化量的相对灵敏系数、B 类相对标准不确定度cr(Vfb)、ur(Vfb)封闭系统总净容积的相对灵敏系数、B 类相对标准不确定度cr(pjc)、ur(pjc)封闭系统检测压力的相对灵敏系数、B 类相对标准不确定度cr(t)、ur(t)检测时间的相对灵敏系数、B 类相对标准不确定度根据式(2)，定义相对灵敏系数 cr(xi)计算式为3 24:crx()i=xiqqxi(4)式中cr(xi)输入量 xi
9、的相对灵敏系数xi 输入量qxi 泄漏率 q 对输入量 xi的偏导数由式(4)可得 cr(p)、cr(Vfb)、cr(pjc)、cr(t)的值分别为 1、1、1、1。将上述值代入式(3)，可得泄漏率的 B 类相对标准不确定度 ur(q)计算式为3 18:ur(q)=u2r()p+u2rV()fb+u2rp()jc+u2r()t(5)泄漏率不确定度的 B 类评定本文以不确定度的 B 类评定方法为例对泄漏率相对标准不确定度进行计算分析。根据 JJF1059 12012测量不确定度评定与表示，B 类评定是用非统计分析方法进行的相对标准不确定度评定。此方法需根据检定装置、测量方法、校准检定证13Bww
10、w gasheat cn权亚强，等:燃气表检定气路系统泄漏率的不确定度评定第 43 卷第 2 期书等有关的信息或经验，判断被测量可能值区间的半宽度 a、包含因子 k。B 类相对标准不确定度 ur计算式为3 15:ur=ak(6)式中ur B 类相对标准不确定度a 被测量可能值区间的半宽度区间半宽度 a 一般根据以下信息确定:以前测量的数据;对有关技术资料和测量仪器特性的了解和经验;生产厂提供的技术说明书;校准证书、检定证书或其他文件提供的数据;手册或某些资料给出的参考数据;检定规程、校准规范或测试标准中给出的数据;其他有用的信息等。对于封闭系统中的气体，温度是影响气密性检测的重要因素。由于检定
11、装置位于恒温(20)实验室中，气密性检测期间封闭系统内外温度近似相同，故温度变化对气密性检测的影响可忽略不计。a p 的 B 类相对标准不确定度p 是由压力变送器测量的，装置上使用的压力变送器经中国测试技术研究院检定合格，证书给出其准确度等级为 0 1 级，取 a=0 1%，按均匀分布考虑，k=3。由式(6)可得 ur(p)为 0.057 7%。b Vfb的 B 类相对标准不确定度串联被检表组、管路和滞止容器组成的封闭系统可认为是刚性封闭容器，认为在气密性检测前后，封闭系统的总净容积不变4。故封闭系统总净容积 Vfb的 B 类相对标准不确定度可忽略不计。c pjc的 B 类相对标准不确定度装置
12、上使用的压力变送器经中国测试技术研究院检定合格，证书给出其准确度等级为 0 1 级，取 a=0 1%，按均匀分布考虑，k=3。由式(6)可得ur(pjc)为 0 057 7%。d t 的 B 类相对标准不确定度装置上使用的晶振经中国测试技术研究院检定合格，证书给出其相对扩展不确定度为 7 105，取 a=7 105，k=2。由式(6)可得 ur(t)为0003 5%。ur(q)结果将以上计算结果代入式(5)，可得到泄漏率的 B类相对标准不确定度 ur(q)为 0 081 7%。检定气路系 统 泄 漏 率 的 不 确 定 度 分 量 uleak应 不 大 于0.09%。3 3泄漏率的相对扩展不确
13、定度泄漏率的相对扩展不确定度 Ur计算见式(7)3 23。按均匀分布考虑，取包含因子 k=3，则Ur为 0 14%。Ur=kur(q)(7)式中Ur 泄漏率的相对扩展不确定度3 4泄漏率的扩展不确定度泄漏率的扩展不确定度计算式为3 23:U=qUr(8)式中U 泄漏率的扩展不确定度4泄漏率的不确定度评定算例泄漏率的不确定度评定算例以某地一台负压法临界流文丘里喷嘴式燃气表检定装置为例，设定密封性检测时该装置封闭系统的净容积 Vfb为 0 05 m3。对于流量低区，规定在系统的检测压力为 94 kPa 的条件下封闭 400 s，p 不超过 10 Pa;对于流量高区，规定在系统的检测压力为 90 k
14、Pa 的条件下封闭 50 s，p 不超过 15 Pa。假设该检定装置泄漏率的相对扩展不确定度 Ur为0.14%，根据式(2)、(8)，密封性检测时，计算泄漏率及其扩展不确定度，见表 1。表 1密封性检测时各流量点泄漏率及其扩展不确定度流量点流量点流量/(m3h1)检测压力pjc/kPa检测时间t/s系统压力变化量 p/Pa泄漏率 q/(m3h1)泄漏率的扩展不确定度U/(m3h1)qmin0 01694400104 8 1056 7 108qt0 6009050156 0 1048 4 107qmax6 0009050156 0 1048 4 1075检定装置的密封性检测建议检定装置的密封性检
15、测建议由上述分析和计算可知，泄漏是导致检定装置检测不准的直接原因，泄漏率可直接反映检定装置密封性的可靠程度。为了保证检测的准确性，需采用一组容积与燃气表容积相近的标准密封器对检定装置进行密封性23B第 43 卷第 2 期煤 气 与 热 力www gasheat cn检测，以下情况时检定装置应进行密封性检测。检定装置的温度变送器、压力变送器、湿度变送器等配套仪器拆卸送检后重新装回检定装置时;检定装置的环境温、湿度或其他测量条件发生了大的变化，刚刚恢复;检定装置发生了碰撞、跌落等意外事件后;检定装置停用一段时间后，重新启用时;对检定装置的检定性能有怀疑时等。6结语结语基于检定气路系统泄漏率的测量过
16、程建立测量模型，通过对泄漏率的不确定度进行 B 类评定，全面分析了各个分量来源的不确定度，并给出泄漏率的不确定度评定过程及算例。对检定装置的密封性检测提出了建议，可为燃气表检定装置泄漏率的分析研究及不确定度评定等工作提供参考。参考文献:1 权亚强，侯雪丹 临界流文丘里喷嘴式燃气表检定装置泄漏分析 J 煤气与热力，2022(10):B30 B34 2 中国计量科学研究院，中国航天机电集团二院203 所，广东省计量科学研究所，等 测量仪器特性评定:JJF10942002 S 北京:中国计量出版社，2003 3 江苏省计量科学研究院，中国计量科学研究院，北京理工大学，等 测量不确定度评定与表示:JJF 1059 12012 S 北京:中国质检出版社，2013:9 4 李宏宇，张静，彭光东，等 航天器推进系统漏率测试不确定度评定J 真空科学与技术学报，2020(7):625 629Uncertainty Evaluation of Leakage ateof Verification Gas Path Systemfor Gas MeterQUAN Yaqiang，HOU XuedanAbst