随着我国天然气生产、消费量的快速增长[1],不同气源的天然气已经进入国内长输管网,混合输送的情况越来越多[2],不同气质、发热量的天然气以体积为基础进行贸易结算显然是不合理的。我国已经建立了完善的天然气能量计量技术、标准体系[3],计量设备基本满足实施能量计量的要求,国内天然气贸易方式转向能量计量势在必行。天然气分析系统的性能评价方法和标准在实施能量计量的情况下将得到广泛应用。
在天然气贸易过程中,常用在线或离线气相色谱仪测定天然气的组成[4-6],用于天然气质量控制、物性参数计算、贸易计量数据计算等。只有测得准确的天然气组成,才能保证贸易计量过程中数据的可靠性,从而保证天然气贸易能够顺利进行,保护贸易各方的利益。要保证测量数据的准确性,首先需确保气相色谱仪处于良好的运行状态,能够测得真实的气体组成。因此,天然气气相色谱分析系统性能评价方法、技术及标准在天然气输送、贸易过程中具有重要的应用价值[7]。
用于天然气气相色谱性能评价的现行标准为GB/T 28766-2012《天然气-在线分析系统性能评价》[8],该标准是对ISO 10723-1995《天然气-在线分析系统性能评价》的修改采用, 本文主要对GB/T 28766-2012的应用进行介绍。
性能评价用气体标准物质(GB/T 28766-2012中称为“试验气体”)是用来获取分析仪器对天然气中单独组分响应值的气体混合物,以便用被评价仪器得到的测量值与标准气中已知的组分含量进行比较,考察仪器性能。计量站现场所安装的性能评价用气体标准物质可采用已较为成熟的分压法制备[10],其不确定度水平基本可满足要求[11]。
在理想情况下,应制备尽量多的气体标准物质,以便在性能评价过程中得到仪器对所测量各组分的响应函数。对新型号的仪器应进行7个含量值试验,对气源气质相对稳定的仪器应进行至少5个含量值试验。这些含量值应尽可能覆盖被评价仪器所测气质组成的变化范围,且其中有1个点低于、1个点高于规定的范围。如果分析系统日常所测样品各组分的最低和最高含量分别为xL%和xU%,则标准物质的含量计算公式如表 1所列。这里需要注意的是,编号1的公式计算结果可能为负,此时将其置为0即可。
GB/T 13610-2014《天然气的组成分析气相色谱法》对气体标准物质组分含量的选择做了规定,CH4采用差减法定量,其他组分的含量范围一般较窄,检测器对这些组分的检测响应值与含量之间满足线性关系。因此,在现场应用时,一般可省去对气体标准物质含量和仪器检测响应值之间线性关系的考察,其余评价项目及指标要求见表 2。
将含N2、CO2、C1~C6全组分的气体标准物质作为样品气,通入在线气相色谱仪,待分析完毕,观察仪器是否能识别气体标准物质所有组分,且对样品中不存在的组分不应给出虚假的色谱峰。另外,在分析过程中电磁阀的切换信号不应对待测组分产生干扰。
连接气体标准物质至分析系统进样口,经充分吹扫进样系统后,连续分析7~10次试验气体,记录每次分析的各组分峰面积,按式(1)计算各组分分析结果的相对标准偏差(RSD),用以评价在线气相色谱仪分析结果的定量重复性。
式中:RSD为相对标准偏差,%;n为测量次数;xi为第i次分析的某组分峰面积;x为某组分峰面积平均值;i为进样分析序号。不同类型仪器峰面积单位表示方法不同,需要注意不同类型仪器的峰面积数据不可直接比较。
根据有效性试验时得到的色谱峰,通过计算分离度,考察各个相邻组分的干扰,计算如式(2)。
式中:R为分离度;tR1, tR2为相邻2个组分的保留时间,s;w1, w2为相邻2个组分色谱峰的峰宽,s。
首先,分析校正用气体标准物质3次,得到各组分的校正因子;然后,分别连续分析不同含量的气体标准物质7~10次;按以下步骤计算气体标准物质中各组分分析结果的不确定度。
(1) B类不确定度:由气体标准物质证书上给出的相对扩展不确定度计算获得。
当气体标准物质证书给出的不确定度为相对扩展不确定度(k=2)时,需换算为绝对不确定度,换算公式如式(3)。
式中:u1为气体标准物质的绝对标准不确定度,%(摩尔分数,下同);mi为组分i的摩尔分数,%;pi为组分i的相对扩展不确定度,%。
(2) A类不确定度:校正用气体标准物质和试验气体重复分析带来的不确定度分别为u2和u3。作为样品气的气体标准物质不确定度u2与作为标准气的气体标准物质不确定度u3计算公式相同,以样品气计算为例。
首先计算样品气标准偏差,计算公式如式(4)。
式中:s为标准偏差,%;n为测量次数;xi为第i次测量的样品气摩尔分数,%;x为样品气摩尔分数平均值,%;i为测量序号。
分析重复性带来的不确定度计算见式(5)。
式中:u3为分析重复性带来的相对不确定度,%;s为单次测量结果的标准偏差,%;n为测量次数。
(3) 分析结果的不确定度。
分析结果的不确定度计算公式如式(6)。
式中:u测为样品气分析结果的标准不确定度,%;u1为标准气的标准不确定度,%;u2为标准气分析重复性引入的标准不确定度,%;u3为样品气分析重复性引入的标准不确定度,%。
根据GB/T 10628-2008《气体分析校准混合气组成的测定和校验比较法》的规定,比较组分分析结果与试验气体参考值之间的一致性,计算公式如式(7)。
式中:x1为试验气体由仪器分析得到的摩尔分数,%;x2为试验气体由标准物质证书中给出的摩尔分数,%;u测为试验气体分析结果的标准不确定度,%;u试为试验气体证书上给出的标准不确定度,%。
当分析结果满足式(7)的要求时,表明在线色谱分析仪的测定结果与试验气体的参考值是一致的;反之,则认为测定结果与参考值不一致。
对某长输管道计量站所配备的在线气相色谱仪开展现场性能评价,以7瓶气体标准物质中3瓶组成有代表性的检测结果为例进行说明。
将作为载气的He通入色谱仪,仪器分析得到的谱图不出现不应有的色谱峰。将含N2、CO2、C1~C6全组分的气体标准物质作为样品气,通入在线气相色谱仪。仪器由两个流路检测所有组分,TCD1检测N2、CO2、C1~C2组分,TCD2检测C3~C6组分;仪器能识别气体标准物质所有组分,且对样品中不存在的组分没有给出虚假的色谱峰。另外,在分析过程中,电磁阀的切换信号也没有对待测组分产生干扰。因此,该分析系统的有效性满足要求。
连接气体标准物质1至分析系统进样口,经充分吹扫进样系统后,连续分析试验气体,记录每次分析的各组分峰面积,并计算各组分的相对标准偏差(RSD),以评价在线气相色谱仪分析结果的定量重复性。根据JJG 1055-2009 《在线气相色谱仪》检定规程中定量重复性RSD≤2%要求,用气体标准物质1通入仪器进行分析,所有组分定量重复性都满足RSD≤2%指标,试验数据见表 3。用气体标准物质2通入仪器进行分析,所有组分定量重复性都满足RSD≤2%指标,试验数据见表 4。用气体标准物质3通入仪器进行分析,所有组分定量重复性都满足RSD≤2%指标,试验数据见表 5。
由表 3~表 5中数据可以看到,该在线色谱仪测量3瓶气体标准物质的重复性均满足标准要求,说明仪器测量的重复性较好。
根据有效性试验得到的色谱峰,通过计算分离度,考察正丁烷、异丁烷相邻组分的干扰。从试验数据可以得出该系统的分离度达到的要求(见表 6),判断标准为JJG 1055-2009要求的分离度R≥1。
由表 6中数据可以看到,该在线色谱仪测量3瓶标准气的分离度均满足标准要求,说明仪器测量的分离度较好。
根据重复性实验中的两瓶气体标准物质,以其中一瓶作样品气,另一瓶作为标准物质进行分析结果的一致性判断依据。气体标准物质各组分的相对不确定度见表 7,其中平衡气CH4的不确定度由各个组分的不确定度按摩尔分数加权计算得到。
以气体标准物质1为样品气,气体标准物质2为校准用气体标准物质,进行试验结果的一致性判断。从试验数据可以看出(见表 8),丙烷和正己烷的一致性没有达到要求。
以气体标准物质2为样品气,气体标准物质1为标准物质,进行试验结果的一致性判。从实验数据可以看出(见表 9),丙烷和正己烷的一致性没有达到要求。
以气体标准物质3为样品气,气体标准物质1为标准物质,进行试验结果的一致性判。从实验数据可以看出(见表 10),氮气、乙烷和正己烷的一致性没有达到要求。
在分析结果与气体标准物质证书上给出的组成的一致性(考虑不确定度)方面,气体标准物质1的丙烷和C6+一致性超标,气体标准物质2的丙烷和C6+一致性超标,气体标准物质3的氮气、丙烷和C6+一致性不满足要求。3瓶气体标准物质检测结果中C6+的一致性均未达到标准要求,经了解设备运行情况后分析,其原因应为该台色谱仪色谱柱超期使用,性能下降所导致。为保证组成分析设备给出准确的分析结果,应定期开展性能评价,并按照设备管理要求更换关键配件,保证设备正常运行。
由现场色谱仪性能评价的情况可见,试验站点的在线色谱仪的有效性、分离度、重复性都能满足标准要求,但在分析结果的一致性方面,氮气、丙烷和C6+存在超标现象。在当前国内天然气贸易以体积计量为主的情况下,这种分析结果的不一致对物性参数和体积计算的影响约为0.05%。C6+的摩尔分数相差0.1%其单位发热量差别为0.16%,如果在实施能量计量的情况下,以天然气的能量作为计价和结算的依据,将对贸易双方产生较大影响。为尽量提高贸易计量数据的准确度,应进一步提高天然气计量站组成分析设备的管理水平,建议至少每2年对气相色谱仪开展一次性能评价,以确保仪器正常运行。