天然气发热量间接测量不确定度评估方法初探

天然气发热量间接测量不确定度评估方法初探

蔡 黎1,2唐 蒙1,2迟永杰1,2

1.中国石油西南油气田公司天然气研究院

2.中国石油集团公司天然气质量控制与能量计量重点实验室

摘要天然气发热量作为天然气的主要使用性能指标，按现代分析计量学要求，发热量值的分析测试需要给出其不确定度。相比使用仪器直接测量天然气发热量值，使用天然气组成计算天然气发热量的间接测量法，也需要进行不确定度评估。使用不确定度传递原则，初步探讨了使用天然气组成及其不确定度进行发热量不确定度评估的方法。并以称量法制备的两元气体标准物质为例，使用气体标准物质的组成不确定度进行发热量不确定度评估，标准物质发热量不确定度为0.1%级别。

关键词发热量间接测量不确定度评估

1天然气发热量测试方法

随着天然气工业的发展，天然气的分析及其物化性质的确定已建立了一套较为完整的体系。天然气自身性质特殊，包括其组分多样、组分含量多变、组分性质相似但不等同等。发热量作为天然气的主要使用性能，其影响因素众多，须使用准确的方法测量，以保证包括发热量在内的天然气物化性质测试结果与真实值一致。在天然气发热量的测量方面，主要的手段有直接测量和间接测量两种[1-3]。

发热量直接测量技术已达到较高水平，但由于发热量直接测量仪器对环境的要求较高，维护工作量较大。因此，目前在管道现场多使用的是发热量间接测量技术，该技术以在线气相色谱仪测量天然气组成，通过计算的方式间接测量天然气发热量。两种测量方法相互验证，均有存在的必要性。ISO在2006年发布了技术报告ISO TR 24094-2006《天然气分析-气体标准物质的验证方法》[4]，用不同等级的发热量直接测量装置分别测试标准物质和天然气样品，对比验证两种方法获取的发热量值，其比对方法见图1、图2。通过此技术报告，进一步确认了发热量间接测量和直接测量两种方法测试结果的一致性，不确定度水平为同一级别。由于使用色谱分析组成数据进行计算的发热量间接测量方法较为简便易行，因而天然气发热量使用间接测量已成为一种趋势[5]。

现代计量对所有测量量值均有不确定度要求。在此要求下，发热量的间接测量结果应给出其不确定度值。尤其是在即将要实施的天然气能量计量中，天然气发热量间接测量结果需要给出其不确定度，才能最终对交接计量的天然气总能量进行不确定度评估。我国的发热量间接测量均使用GB/T 11062-1998《天然气发热量、密度、相对密度和沃泊指数的计算方法》进行[6]，此标准给出发热量间接测试的方法，以及发热量的偏差评估方法，但并未严格地按照JJF 1059.1-2012《测量不确定度评定与表示》中规定的方法进行不确定度评估[7]。所以，建立由天然气组成不确定度向发热量不确定度传递的方法迫在眉睫。本文主要探讨进行天然气发热量间接测试时天然气发热量不确定度的评估方法。

2天然气发热量间接测量不确定度评估方法

2.1天然气发热量间接测量方法

天然气发热量间接测量使用气相色谱获取天然气组成，随后由天然气组成数据进行计算。首先使用公式(1)计算气体标准物质的理想气体发热量。

(1)

对于理想气体，其摩尔分数与体积分数相同，在计算理想气体发热量时，可以使用摩尔分数进行计算。在完成理想气体发热量计算后，使用公式(2)进行天然气真实气体发热量计算。

(2)

式中，Zmix为在计量参比条件下的天然气压缩因子。

而要计算气体标准物质的压缩因子Zmix，需要使用公式(3)。

(3)

在以上计算过程中，进行天然气真实气体发热量不确定度评估时，需使用计量校准规范JJF1059.1-2012中规定的不确定度传播律，首先进行理想气体发热量不确定度的评估，再结合真实气体压缩因子的不确定度进行真实气体发热量的不确定度评估。

2.2理想气体不确定度评估

气体标准物质理想气体发热量按照公式(1)进行计算，而要评估此过程的发热量不确定度，需要各组分的纯物质发热量不确定度和摩尔分数不确定度。摩尔分数的不确定度在进行天然气组成分析的过程获取。

甲烷：GB/T11062-1998中给出甲烷发热量为37.044MJ/m3，但并未给出甲烷发热量的分布范围。德国联邦物理技术研究院(PTB)使用重复性为0.01%的量热计进行纯甲烷发热量测试，甲烷发热量扩展相对不确定度为Urel=0.052%(k=2)，则甲烷发热量的标准绝对不确定度为0.010MJ/m3。本文研究的气体标准物质在进行发热量不确定度评估时，甲烷的发热量也使用此不确定度值。

氮气：氮气作为惰性气体，在天然气燃烧的过程中，不发生任何化学变化，故而氮气的发热量为0，并且氮气为0的发热量不存在分布范围，也就是说，氮气发热量对不确定度的贡献可忽略。

乙烷：GB/T11062-1998给出乙烷发热量为64.91MJ/m3，但未给出乙烷发热量的分布范围。ISO/CD6976-2013《天然气发热量、密度、相对密度和沃泊指数的计算方法》中给出乙烷发热量相对不确定度urel=0.032%[8]。则乙烷发热量的绝对不确定度为0.021MJ/m3。本文研究的气体标准物质在进行发热量不确定度评估时，乙烷的发热量使用此不确定度值。

除以上组分外，常用天然气组分中各类烷烃、氢气的发热量和发热量的不确定度在ISO/CD6976-2013中均有介绍，而二氧化碳和氦气等组分无法燃烧，发热量为零，发热量的不确定度贡献可忽略。在进行理想气体发热量不确定度评估的过程中，使用JJF1059.1-2012中规定的不确定度传播律进行。由于各量值之间不相关，合成发热量不确定度时，不确定度传播律公式可简化为：

(4)

2.3真实气体压缩因子不确定度

表1 双组分模拟气体压缩因子计算结果Table1 Compressionfactoroftwo-componentsimulationgas组成摩尔分数/%参比条件下压缩因子氮气/甲烷氮气甲烷6.70393.2970.998256.46493.5360.9982416.37683.6240.9984617.36282.6380.99848乙烷/甲烷乙烷甲烷6.04593.9550.997856.02393.9770.9978511.44488.5560.9976211.62488.3760.99761 注:参比条件为20℃,101.325kPa。

2.4真实气体发热量不确定度

由于气体标准物质真实气体发热量按照公式(2)来执行，在计算的过程中，所使用的计算公式只存在乘积关系，故不确定度评估公式可以简化为公式(5)，真实气体发热量的相对不确定度等于理想气体发热量相对不确定度和参比条件下真实气体压缩因子相对不确定度的方和根：

(5)

式中，urelH真实为真实气体发热量相对不确定度；urelH理想为理想气体发热量相对不确定度；urelZ为真实气体压缩因子相对不确定度。

3天然气发热量间接测量不确定度评估实例

本文以称量法制备的双组分气体标准物质为例，进行真实气体发热量不确定度评估，评估结果见表2。

表2 真实气体发热量评估结果Table2 Resultsofuncertaintyestimationforrealgascalorificvalue组分摩尔分数/%urelx/%H0/(MJ·m-3)uH0/(MJ·m-3)urelH0/%ZmixuZmix/%H/(MJ·m-3)urelH/%N26.7030.095CH493.2970.01334.5610.0110.0320.998250.0534.6220.06 注:urelx为摩尔分数相对不确定度;H0为理想气体发热量;uH0为理想气体发热量不确定度;urelH0为理想气体发热量相对不确定度;uZmix为真实气体压缩因子相对不确定度;H为真实气体发热量;urelH为真实气体发热量相对不确定度。

评估过程如下：

氮气摩尔分数的绝对不确定度为：uN2=0.095%×6.703%≈0.006 4%

甲烷摩尔分数的绝对不确定度为：uCH4=0.013%×93.297%≈0.012%

≈0.011 MJ/m3

≈0.06%

按照以上的方法，可以进行天然气发热量间接测量标准不确定度的评估，在进行扩展不确定度评估时，仅需要乘以扩展因子(包含因子)即可。

由评估结果可以看到，发热量间接计算按照本文给出的方法标准，其相对不确定度为0.06%，绝对不确定度约为0.02 MJ/m3。GB 11062-1998中仅对发热量计算的准确度进行描述，其准确度不可能高于0.01 MJ/m3，但在没有实际推导过程的情况下，0.01 MJ/m3准确度是否能够保证还是未知数。并且本文在评估过程中，参比条件下压缩因子作为最终评估结果的主要贡献，其不确定度的评估相对保守。通过更为精确的压缩因子数据，最终获取的发热量不确定度可能更低。

4结论与建议

(1) ISO VAMGAS 项目已证明天然气发热量直接测量和间接测量结果一致。本文论证了天然气发热量间接测量不确定度为0.06%，与Cutler-Hammer发热量直接测定装置不确定度0.25%(k=2)处于同一级别。

(2) 我国现有标准进行发热量间接测量时，发热量不确定度评估方法有所缺失。本文描述的不确定度评估方法可以替代现有国家标准GB/T 11062-1998中规定的发热量偏差评估方法。

(3) 本文描述的不确定度评估方法适用于天然气分析用气体标准物质发热量的不确定度评估。用于评估使用称量法获取的气体标准物质时，最终发热量不确定度水平仅为0.1%级别。

参 考 文 献

[1] 李克, 潘春锋, 张宇, 等. 天然气发热量直接测量及赋值技术[J]. 石油与天然气化工, 2013, 42(3): 297-301.

[2] 王海峰， 李佳， 孙国华， 等. 基准气体热量计研究进展[J]. 石油与天然气化工, 2014, 43(2): 196-199.

[3] 陈赓良, 唐飞. 管输压力下实时记录式天然气发热量测定设备的开发[J]. 石油与天然气化工, 2013, 42(2): 181-187.

[4] ISO TR 24094: 2006.Analysis of natural gas-Validation methods for gaseous reference materials[S]. 2006.

[5] 全国天然气标准化技术委员会(SAC/TC 244). GB 17820-2012天然气[S]. 北京： 中国标准出版社, 2012-09-01.

[6] 中国石油天然气总公司. GB 11062-1998天然气发热量、密度、相对密度和沃泊指数的计算方法[S]. 北京： 中国标准出版社， 2004-04-09.

[7] 全国法制计量管理计量技术委员会. JJF 1059.1-2012测量不确定度评定与表示[S]. 北京： 中国标准出版社， 2013-06-03.

[8] ISO/CD 6976: 2013 Natural gas-Calculation of calorific values, density, relative density and Wobbe index from composition[S]. 2013.

Discussion on uncertainty estimation of indirectly measuring natural gas calorific value

Cai Li1,2, Tang Meng1,2, Chi YongJie1,2

(1.ResearchInstituteofNaturalGasTechnology,PetroChinaSouthwestOil&GasfieldCompany,Chengdu610213,China; 2.KeyLaboratoryofNaturalGasQualityandEnergyMeasurement,CNPC,Chengdu610213,China)

Abstract:Calorific value is one of the most important indexes of natural gas. Modern metrology requires all the testing results including calorific value exhibited with uncertainty. Comparing direct measurement of calorific value by heat-determiner, the indirect measurement of natural gas calorific value by calculation from the composition of natural gas requires the uncertainty estimation. In the article, the uncertainty estimation of calorific value from the uncertainty of composition is discussed. Using reference gas prepared by gravimetric method with two components as an example, the uncertainty of the calorific value is about 0.1%.

Key words:calorific value, indirect measurement, estimation of uncertainty

收稿日期：2014-10-13；编辑：钟国利

中图分类号：TE642

文献标志码：A

DOI:10.3969/j.issn.1007-3426.2015.02.022

作者简介：蔡黎(1983-)，男，四川大学物理化学专业毕业，博士。现就职于中国石油天然气股份有限公司西南油气田公司天然气研究院，主要从事天然气分析用标准物质和非常规天然气气质方面的研究。E-mail:cai_li@petrochina.com.cn

基金项目：中国石油西南油气田公司科研项目“甲烷中含多组分混合气体及发热量直接测定用气体标准物质研究”(20130308-05)。