油库质量流量计现场检定测量不确定度评定

黄智渊

江西省宜春市计量所　（江西　宜春　336000）

科里奥利质量流量计是一种直接式质量流量仪表，它具有准确度高、稳定性好、量程比大、介质适应性强等特点，且不受被测介质物理参数和管内流体流动状态的影响，能够实现对质量流量直接测量[1]，广泛应用于油库中成品油的计量。依据JJG 1038—2008《科里奥利质量流量计检定规程》，使用静态容积法油流量标准装置，利用分项组合法对成品油贸易结算用科里奥利质量流量计进行现场检定，并就被检流量计示值误差的测量不确定度进行评定。

1　概述

1）测量依据：JJG 1038—2008《科里奥利质量流量计检定规程》。

2）环境条件：环境温度一般为5~45℃，相对湿度一般为35%~95%，大气压力一般为86~106 kPa。

3）测量标准：二等标准金属量器（容积为2 000 L），数字防爆温度计（分度值为0.1℃），工作用玻璃液体温度计（分度值为0.1℃），二等标准石油密度计，秒表（分度值为 0.01 s）。

4）测量对象：科里奥利质量流量计（0.2级）。

5）测量介质：成品油。

6）测量方法：在规定的环境条件下，按照检定规程规定的测量程序，以油库工况下固定的发油流量点发油，测量3次，计算得到被检流量计的相对示值误差。

2　数学模型及不确定度传播率

2.1　数学模型

为使测量结果不确定度评定过程更具逻辑性，采用分步的方式构建数学模型[2-3]。

1）数学模型 1，计量温度下的油品体积（Vt）：

2）数学模型 2,油品在 20℃下的标准体积（V20）：

3）数学模型3，考虑空气浮力修正后的油品计量密度（ρ）：

4）数学模型4，油品在空气中的标准计重质量（Qs）：

5）数学模型5，流量计相对误差（E）：

式中：Vs为标准金属量器读出容积，L；βs为标准金属量器的体膨胀系数，℃-1；ts为标准金属量器处液体温度 （计量温度），℃；VCF为体积修正系数；ρ20为油品在20℃下的标准密度，g/cm3；Q为被检流量计累积质量流量，kg。

2.2　不确定度传播率

不确定度传播率计算如下：

由于 Q≈Qs，得

3　各输入量的标准不确定度

3.1　计量温度下油品体积的相对标准不确定度

依据数学模型1计算下列参数[4-8]。

3.1.1标准器引入的相对标准不确定度分量ur（Vs）

依据标准器的检定证书和JJG 2024—1989《容量计量器具检定系统》规定，二等标准金属量器的扩展不确定度为 2.5×10-4，k=3。 因此：

3.1.2标准器体膨胀系数引入的相对标准不确定度

分量 ur（βs）

标准金属量器的体膨胀系数βs=50×10-6℃-1，查资料，其不确定度为5×10-6℃-1，符合均匀分布。因此：ur（βs）=

依据检定经验，考虑气候条件影响，检定时油库发油温度控制在（20±10）℃范围内。取极限温度值：ts＝30 ℃。

相对灵敏系数为：

3.1.3油温引入的相对标准不确定度分量ur（ts）

依据数字温度计的校准证书，温度测量值的扩展不确定度为0.2℃，k=2。取极温度限值：ts＝30℃。因此：

相对灵敏系数为：

3.1.4相对标准不确定度ur（Vt）

由数学模型1得到相对标准不确定度ur（Vt）：

3.2　密度测量和密度相关参数不确定度分析

依据GB/T 1884—2000标准中的方法，使用二等标准石油密度计和工作用玻璃液体温度计，分别测量油品的视密度（ρ′）和试验温度（t′），并查 GB/T 1885—1998标准中表59B得到油品20℃下的标准密度（ρ20）。使用数字温度计测量标准金属量器内油温（ts），查GB/T 1885-1998标准中表60B得到体积修正系数VCF。为减小因查表所引起的误差，查表所得的标准密度（ρ20）和体积修正系数VCF均使用双线性内插法求得。

3.2.1油品视密度的标准不确定度u（ρ′）

1）二等标准石油密度计的修正值引入的标准不确定度分量u（ρ1）。在油品密度测量时，密度计读数需要修正，依据检定证书（校准证书），其修正值的扩展不确定度U=0.000 20 g/cm3，k=2。因此：

2）密度计估读引入的标准不确定度分量u（ρ2）。二等标准石油密度计的分度值为0.000 5g/cm3，GB/T 1884—2000要求的读数间隔为分度值的1/5，按均匀分布计算。因此：

3）密度计倾斜读数引入的标准不确定度分量u（ρ3）。 按照 JJG 86—2011《标准玻璃浮计检定规程》的规定，浮计干管与液面间的垂直偏差不得大于0.1个分度值，服从均匀分布。因此：

4）密度计测量过程中油液温度变化引入的标准不确定度分量 u（ρ4）。 GB/T 1884—2000 规定，测量过程中油液温度最大变化为0.5℃。油液温度变化引起的密度计示值变化为：Δρ=ρ′·Δt·βM， 服从均匀分布。其中ρ′取常用石油密度计的最大刻度0.850 g/cm3，Δt取最大允许温度变化 0.5 ℃，βM为密度计体膨胀系数 25×10-6℃-1。

5）密度计修正引入的标准不确定度分量u（ρ5）。依据GB/T 1884—2000要求，对观察到的密度计读数进行修正后，记录到0.000 1 g/cm3。修约间隔为0.000 05 g/cm3，服从均匀分布。因此：

6）合成标准不确定度u（ρ′）。计算结果如下：

7）扩展不确定度U（ρ′）。油品视密度的扩展不确定度为：U（ρ′）≈0.000 3g/cm3，k=2。

3.2.2试验温度的标准不确定度u（t′）

1）工作用玻璃液体温度计引入的标准不确定度分量u（t1）。在油品试验温度测量时，温度计读数需要修正，依据检定证书（校准证书），其修正值的扩展不确定度为U=0.06℃，k=2。因此：

2）温度计读数引入的标准不确定度分量u（t2）。依据GB/T 1884—2000要求，温度计读数时记录温度接近到0.1℃。假设近似取值间隔为0.05℃，服从均匀分布。因此：

3）温度计修正引入的标准不确定度分量u（t3）。依据GB/T 1884—2000要求，对观察到的温度计读数进行修正后，记录到接近0.1℃。修约间隔为0.05℃，服从均匀分布。因此：

4）合成标准不确定度 u（t′）。 计算结果如下：

5）扩展不确定度U（t′）。试验温度的扩展不确定度为：U（t′）≈0.1 ℃，k=2。

3.2.3标准密度的标准不确定度u（ρ20）

假设油品的视密度取值与GB 19147—2016和GB 17930—2016标准中规定相同，则车用柴油和车用汽油标准密度分别为0.810~0.850 g/cm3和0.720~0.775 g/cm3。假设试验温度取值和油库发油时的温度范围相同，为（20±10） ℃。

1）试验温度测量引入的标准不确定度 ut（ρ20）。根据上述计算，试验温度的扩展不确定度为0.1℃。查GB/T 1885—1998表59B,得到由此引入的标准密度最大变化为0.000 1 g/cm3，服从均匀分布。因此：

2）视密度测量引入的标准不确定度 uρ′（ρ20）。 根据上述计算，视密度的扩展不确定度0.000 3 g/cm3。查GB/T 1885—1998表59B,得到由此引入的标准密度最大变化为0.000 3 g/cm3，服从均匀分布。因此：

3）标准密度计重复性引入的标准不确定度分量uρ20（ρ20）。 参照标准 GB/T 1884—2000，两次测量结果的极差最大为0.000 2 g/cm3。最终测量结果取两次测量结果的平均值。因此：

4）合成标准不确定度u（ρ20）。计算结果如下：

5）扩展不确定度U（ρ20）。标准密度的扩展不确定度为：U（ρ20）≈0.000 5 g/cm3， k=2。

3.2.4油品计重密度的相对标准不确定度ur（ρ）

由数学模型 3，可知：u（ρ）=u（ρ20），因此测量点的最大相对标准不确定度为：

3.2.5体积修正系数的相对标准不确定度ur(VCF)

依据油品的标准密度（ρ20）和计量温度（ts），查GB/T 1885—1998表60B得到体积修正系数VCF。

1）试验温度测量引入的标准不确定度ut(VCF)。数字温度计的扩展不确定度0.2℃。查GB/T 1885—1998表60B得到由此引入的体积修正系数的最大变化为0.000 32，服从均匀分布。因此：

2）标准密度测量引入的标准不确定度uρ(VCF)。根据上述计算，标准密度的扩展不确定度0.000 5 g/cm3。查GB/T 1885—1998表60B得到由此引入的体积修正系数的最大变化为0.000 03，服从均匀分布。因此：

3）合成标准不确定度u(VCF)，计算结果如下：

4）体积修正系数的相对标准不确定度ur(VCF)。体积修正系数的最大变化是在标准密度0.720 g/m3，计量温度10℃附近取得。因此：

3.3　油品在空气中的标准计重质量引入的相对标准不确定度 ur（Qs）

由数学模型 2、3、4，可得：

Qs=V20·ρ=Vt·VCF·ρ。

相对标准不确定度 ur（Qs） 为：

3.4　被检流量计累积质量流量测量结果的重复性引入的相对标准不确定度ur（Q）

选取一台常规的科里奥利质量流量计（0.2级），以油库工况下固定的发油流量点发油，测量3次，测量结果如表1所示。

4　标准不确定度汇总

根据上述测量不确定度评定过程，将各输入量的标准不确定度进行汇总，如表2所示。

表1　科里奥利质量流量计（0.2级）工况发油流量点的测量结果

表2　标准不确定度汇总表

5　合成标准不确定度

各项相对标准不确定度分量彼此相互独立，依据不确定度传播率，求得合成标准不确定度为：

6　扩展不确定度

取k=2，则流量计累积流量相对误差的扩展不确定度为：U=uc×k=0.090%。

7　测量结果不确定度报告

被检流量计工况下固定发油流量点的累积流量相对误差为0.005%，其扩展不确定度为：U=0.090%,k=2。

为简化评定过程，在计算灵敏系数时，未根据实测温度计算，而是直接取极限温度值；在密度测量和密度相关参数不确定度分析过程中，未具体区分油品，而是直接取极限值进行评定。在实际过程中，可根据具体情况评定上述分量，获得更准确的测量结果不确定度。

8　结束语

由测量不确定度评定过程可知，现场检定的测量结果不确定度小于被检流量计最大允许误差的1/2。该静态容积法油流量标准装置符合国家计量标准JJG 2063—2007《液体流量计器具检定系统表检定规程》和JJG 1038—2008《科里奥利质量流量计检定规程》的相关要求，可以用于现场检定准确度等级0.2级及以下的科里奥利质量流量计。

参考文献：

[1]苏彦勋，梁国伟，盛健．流量计量与测试[M]．2版.北京：中国计量出版社，2007．

[2]中国石油天然气股份有限公司．成品油计量员手册[M]．北京：中国石油天然气股份有限公司，2009．

[3]叶培德．测量不确定度理解评定与应用[M]．北京：中国质检出版社，2013．

[4]纪建英．流量容量计量标准技术报告[M]．北京：中国计量出版社，2009．

[5]黄锋，谭山．力学计量器具建标指南[M]．北京：中国质检出版社，2013．

[6]杨有涛，陈梅．流量标准装置建标指南[M]．北京：中国质检出版社，2013．

[7]中国石油化工集团公司石油化工科学研究院技术.原油和液体石油产品密度实验室测定法（密度计法）：GB/T 1884—2000[S]．北京：中国标准出版社，2004.

[8]石油静态和轻烃计量标准化技术归口单位．石油计量表产品部分（GB/T 1885—1998 专用）[M]．北京：中国标准出版社，2005．