石油与天然气化工  2018, Vol. 47 Issue (6): 80-85, 91
天然气中C6+组成分析不同定量方式对计算烃露点的影响
曾文平 1,2, 蒲长胜 1,2, 王晓琴 1,2, 张镨 1,2, 王伟杰 1,2, 何飞 1     
1. 中国石油西南油气田公司天然气研究院;
2. 中国石油天然气质量控制和能量计量重点实验室
摘要:由天然气组成数据计算烃露点是获取烃露点的一种重要的检测方法,天然气组成分析已有相关方法标准,但对于C6+组成分析有不同的定量和数据处理方式。概述了计算烃露点的主要影响因素, 重点分析了C6+延伸组成分析不同定量方式对计算烃露点结果的影响。采用气相色谱法分析了两个实际天然气样品,根据不同的定量方法和数据处理方式获得的C6+延伸组成数据预测了天然气的烃露点。结果表明, 对于C6+延伸组成的分析,C6~C16按碳数归类,其中苯、甲苯、环己烷和甲基环己烷等组分单独定量时,由此获得完整天然气组成数据计算的烃露点和直接测定烃露点平均偏差在3 ℃以内,初步验证了计算法和直接测定法获得的烃露点结果的一致性。
关键词C6+延伸组成    架桥定量    碳数归类    烃露点    特殊组分    
Effect of natural gas C6+ composition obtained by different quantitative methods on calculation of hydrocarbon dew point
Zeng Wenping1,2 , Pu Changsheng1,2 , Wang Xiaoqin1,2 , Zhang Pu1,2 , Wang Weijie1,2 , He Fei1     
1. Research Institute of Natural Gas Technology, PetroChina Southwest Oil & Gasfield Company, Chengdu, Sichuan, China;
2. Key Laboratory of Natural Gas Quality Control and Energy Measurement, PetroChina, Chengdu, Sichuan, China
Abstract: The calculation of hydrocarbon dew point from natural gas composition data is an important detection method for obtaining hydrocarbon dew point. Natural gas composition analysis has related methods and standards, but there are different quantitative and data processing methods for C6+ composition analysis. In this paper, the main factors influencing the calculation result of hydrocarbon dew point are summarized, among which the difference between different quantitative methods of C6+ extension analysis for predicting hydrocarbon dew point is emphatically discussed. Two real natural gas samples were analyzed by gas chromatography. The experimental data were obtained by different quantitative and data processing methods for C6+ extension components, and the hydrocarbon dew point of natural gas is predicted by using composition data of natural gas. The results show that C6+ components were classified according to the carbon number from C6 to C16, in which benzene, toluene, cyclohexane and methyl cyclohexane were respectively determined, and the complete analytical data of natural gas composition were obtained. The average deviation between the hydrocarbon dew point calculated through natural gas composition data and directly measured hydrocarbon dew point is less than 3 ℃.The consistency of hydrocarbon dew point obtained by calculation method and direct determination method is preliminarily verified.

在天然气的压力、温度发生变化时,当天然气的温度低于一定压力下的烃露点后,天然气中较重烃类就会析出形成液相,这些液相组分的存在将对管道输送、天然气利用产生不同程度的不利影响[1-3]。因此,在各国天然气的气质标准及天然气管道设计规范中, 除对气体中的H2S、CO2、水露点、发热量等有具体要求外,同时也对烃露点指标进行了规定。

测定天然气的烃露点主要有仪器直接测定和通过组成数据由专用软件计算两种方法。两种方法各有优缺点,直接测定方法虽然更方便、快速,但测定时受天然气中水的干扰较大,当烃露点低于水露点时,消除水对测定烃露点的干扰是准确获得烃露点的关键因素。通过软件计算天然气的烃露点,可获得天然气的相包络线,由此能够得到不同压力下天然气的烃露点,但通过组成数据计算烃露点时,影响计算结果准确度的因素较多。本文概述了计算烃露点的主要影响因素, 重点针对分析C6+延伸组分时不同定量方法和数据处理方式获得组成结果对计算烃露点的影响开展研究,并通过直接测定的烃露点验证不同定量方法获得结果的合理性和准确性,以期获取一种相对可靠的C6+延伸组分分析定量和数据处理方式,为采用天然气组成数据准确预测烃露点或烃含量提供技术参考。

1 预测烃露点的影响因素

计算法获得烃露点通过现场获取实际天然气样品,然后采用气相色谱法分析天然气的组成,再由天然气组成数据通过软件计算天然气的烃露点。由于天然气是一种组分十分复杂的混合物,影响计算法获得烃露点结果的因素较多,其主要影响因素包括样品的代表性、组成分析结果的准确性以及计算模型的选择等。

采用GB/T 13609-2017《天然气取样导则》规定的方法进行取样[4],现场取样示意图见图 1。取样钢瓶在实验室经过特殊处理,在取样过程中,由于压力的突降,取样管线和接头等部位可能造成节流,节流处会急剧降温,若不采取一定的措施则会造成重烃的析出、损失,导致取样不准确,使取到的样品没有代表性。因此,取样时应对取样探头、取样管线和取样阀等部位采取加热保温措施。

图 1     用于天然气延伸分析样品的取样示意图 Figure 1     Sampling diagram for natural gas extension analysis samples

由于天然气中组分较多,要准确分析每个组分的含量十分困难。天然气组成全分析可分为两部分完成,对于N2、CO2、C1~C5烃类等主要组分(以下简称常量组分)含量的分析采用GB/T 13610-2014《天然气的组成分析气相色谱法》[5]、GB/T 27894《天然气在一定不确定度下用气相色谱法测定组成》系列标准[6-11],而对于C6及更重烃类组分(以下简称C6+延伸组分)含量的分析采用GB/T 17281-2016《天然气中丁烷至十六烷烃类的测定气相色谱法》规定方法[12]。对于常量组分的分析,由于采用了气体标准物质外标定量,一般可以获得准确可靠的组成数据,但对于C6+延伸组分含量的分析,由于C6+延伸组分较多,要准确分析每个组分的含量较困难。天然气组成分析结果的准确可靠是影响计算烃露点结果的主要因素,特别是C6+延伸组分含量的分析结果,对计算的烃露点有着重要的影响。以下将重点分析研究C6+延伸组分不同定量方法和数据处理方式获得结果对计算烃露点的影响,并通过直接测定的烃露点验证不同定量方法获得结果的合理性和可靠性。

计算烃露点的软件主要是基于SRK和PR等状态方程,曾有研究表明,分别采用CHEMCAD 5.2和HYSYS 3.0两个软件计算不同压力下气体标准物质的烃露点时,得到的烃露点结果相差均在1 ℃以内[2],说明不同软件计算的烃露点结果是基本一致的。

2 天然气中C6+组成分析
2.1 定量方法
2.1.1 常量组分

具有代表性的气样和已知组成的气体标准物质在相同的操作条件下,用气相色谱法进行分离。样品中重烃组分可以在某个时间通过改变流过色谱柱载气的方向,获得一组色谱峰,这组重烃组分是C6和更重组分的加和峰。由气体标准物质的组成值,通过对比样品气和气体标准物质的响应值,计算获得样品的相应组成。

2.1.1.1 戊烷和更轻组分

测量每个组分的峰面积,比较气样和标准气中相应组分的响应值,按式(1)计算气样中i组分的含量:

$ {y_i}{\rm{ = }}{y_{si}}\left( {{A_i}/{A_{si}}} \right) $ (1)

式中:yi为气样中i组分的摩尔分数, %; ysi为标准气中i组分的摩尔分数,%;Ai为气样中i组分的峰面积, μV·s;Asi为标准气中i组分的峰面积,μV·s。

2.1.1.2 己烷和更重组分加和峰

测量反吹的己烷及更重组分的峰面积,并在同一色谱图上测量正戊烷、异戊烷的峰面积,按式(2)计算气样中C6+的含量:

$ y\left( {{\rm{C}}_6^ + } \right) = \frac{{y\left( {{{\rm{C}}_5}} \right)A\left( {{\rm{C}}_6^ + } \right)M\left( {{{\rm{C}}_5}} \right)}}{{A\left( {{{\rm{C}}_3}} \right)M\left( {{\rm{C}}_6^ + } \right)}} $ (2)

式中:y(C6+)为气样中C6及更重组分加和峰的摩尔分数,%;y(C5)为气样中异戊烷与正戊烷摩尔分数之和,%;A(C6+)为气样中C6和更重组分加和峰的峰面积, μV·s;A(C5)为气样中异戊烷和正戊烷的峰面积之和, μV·s;M(C5)为戊烷的相对分子质量,取值为72;M(C6+)为C6和更重组分加和峰的相对分子质量,取值为92。

2.1.2 C6+延伸组分
2.1.2.1 用丁烷标准气定量(方法一)

如果使用外标法定量,应测量标准气和样品气两者中各组分的峰面积。样品气中碳数为ni组分的摩尔分数按式(3)计算。

$ {y_{i, n}} = {y_{{{\rm{C}}_4}}}\left( {\frac{{4{A_{i, n}}}}{{n{A_{\rm{C}}}_{_4}}}} \right) $ (3)

式中:yi, n为样品气中碳数为ni组分的摩尔分数, %;yC4为标准气中丁烷的摩尔分数,%;Ai, n为样品气中i组分的峰面积, pA·s;AC4为标准气体混合物中丁烷的峰面积, pA·s;n为重烃组分的碳数。

2.1.2.2 用样品气中戊烷架桥定量(方法二)

如果GB/T 17281-2016中规定的方法是和GB/T 13610-2014中所规定的方法一起使用,应把GB/T 13610-2014得到的未归一化结果与GB/T 17281-2016得到的结果相结合,体现在一个测试报告中。戊烷作为架桥组分时,样品气中碳数为n(6≤n≤16)的i组分摩尔分数按式(4)计算。

$ {y_{i, n}} = {y_{{{\rm{C}}_5}}}\left( {\frac{{5{A_{i, n}}}}{{n{A_{\rm{C}}}_{_5}}}} \right) $ (4)

式中:yi, n为样品气中碳数n(6≤n≤16)的i组分摩尔分数, %;yC5为由GB/T 13610-2014测定的样品气中正戊烷的摩尔分数,%;Ai, n为样品气中i组分的峰面积,pA·s;AC5为样品气中正戊烷的峰面积,pA·s。

2.2 数据处理方式
2.2.1 C6+所有组分均按碳数归类(处理方式1)

由于天然气组分的复杂性,计算烃露点时,C6及更重组分需按碳数归类进行定量分析,即正戊烷(n-C5)和正己烷(n-C6)之间的组分归入C6,其中不包括nC5,但包括n-C6,其他重烃组分依次类推,按式(5)计算。

$ {y_{\rm{C}}}_{_n} = \sum\limits_{i = 1}^k {{y_i}{{_{\rm{C}}}_{_n}}} $ (5)

式中:yCn为按碳数归类的重烃组分的摩尔分数,%,即在色谱图上保留时间为碳数(n-1)的正构烷烃和碳数为n的正构烷烃之间所有组分含量之和, 其中不包括(n-1)的正构烷烃, 但包括n的正构烷烃的含量;yiCn为碳数为ni组分的摩尔分数,%。

2.2.2 C6+中特殊组分单独定量(处理方式2)

当采用组成数据,通过专用软件计算天然气的烃露点时,GB/T 30492-2014《天然气烃露点计算的气相色谱分析要求》要求对于C6+组分的定量分析,当按碳数归类有交叉的情况时,应将苯、甲苯、环己烷、甲基环己烷等组分进行单独定量分析[15]

对于苯、甲苯、环己烷和甲基环己烷等组分的定量分析,由于在色谱图中,苯和环己烷出峰的保留时间在n-C6和n-C7之间,甲苯和甲基环己烷在n-C7和n-C8之间出峰,对这几个组分的定量可采用式(6)计算。

$ {y_{i, n}} = {y_{{{\rm{C}}_5}}}\left( {\frac{{5{A_{i, n}}}}{{n{A_{\rm{C}}}_{_5}}}} \right) $ (6)

式中:yC5为由GB/T 13610-2014测定的样品气中正戊烷的摩尔分数, %;Ai, n为样品气中i组分(苯、甲苯、环己烷和甲基环己烷)的峰面积, μV·s;AC5为样品气中正戊烷的峰面积, μV·s;n为重烃组分的碳数,对于苯和环己烷,n=7,对于甲苯和甲基环己烷,n=8。

3 结果与讨论
3.1 结果
3.1.1 不同定量分析方法获得C6+组成

采用气相色谱法分析了两个实际天然气样品,常量组分采用GB/T 13610-2014,C6+延伸组分采用GB/T 17281-2016。其中,对于C6+各组分的定量,分别采用气体标准物质中n-C4架桥定量分析(方法一)和样品气中n-C5架桥定量分析(方法二),表 1为C6+延伸组分采用两种定量分析方法获得的结果。

表 1    两种定量分析方法获得1号和2号样品组成结果比较(C6+组分按碳数归类) Table 1    Comparison of the results of composition of sample No.1 and No.2 obtained by two quantitative analysis methods (C6+ components classified by carbon number)

结果表明, 分别采用气体标准物质中n-C4架桥定量和样品气中n-C5架桥定量分析两个样品C6+延伸组分时,其绝对偏差最大值为C6和C7的分析结果,均为0.000 5%,说明两种定量分析方法获得的结果是一致的。

此外,对于C6+延伸组分除只按碳数进行归类外(处理方式1),一般还将苯、甲苯、环己烷、甲基环己烷等组分进行单独定量分析(处理方式2),两个实际天然气样品按照数据处理方式2获得的结果见表 2

表 2    按照两种数据处理方式获得的1号和2号样品分析结果(C6+特殊组分单独定量) Table 2    Comparison of the results of composition analysis of sample No.1 and No.2 obtained by two data processing methods(C6+ special components single quantification)

结果表明, 将苯、甲苯、环己烷、甲基环己烷等组分进行单独定量分析时,分别采用气体标准物质中n-C4架桥定量分析和样品气中n-C5架桥定量分析时,其中甲基环己烷分析结果的偏差最大,为0.001 2%。

3.1.2 由不同定量和数据处理方式获得C6+组成计算烃露点

采用专用软件计算了两个实际天然气样品不同压力下的烃露点,见表 3表 4

表 3    不同压力下1号样品计算的烃露点 Table 3    Results of hydrocarbon dew point for sample No.1 at different pressure

表 4    不同压力下2号样品计算的烃露点 Table 4    Results of hydrocarbon dew point for sample No.2 at different pressure

结果表明:C6+组成按碳数归类时,采用标准气中n-C4架桥定量和样品气中n-C5架桥定量时,计算的烃露点是一致的,两者相差在0.2 ℃以内;而C6+中特殊组分单独定量分析时,采用标准气中n-C4架桥定量和样品气中n-C5架桥定量,计算的烃露点也是一致的,两者相差也在0.2 ℃以内;但采用特殊组分单独定量时,计算的烃露点比C6+组成只按碳数归类时计算的烃露点低,最大偏低值达到6 ℃左右,说明采用两种数据处理方式获得的C6+组成结果对计算的烃露点有较大影响。

3.2 讨论
3.2.1 不同定量方法获得C6+组成结果对计算烃露点的影响

采用标准气中n-C4架桥定量和样品气中n-C5架桥定量两种方法获得的天然气中C6+延伸组成分析结果是一致的,而采用不同定量方法获得的组成结果计算获得的烃露点也是一致的,两者差值在0.2 ℃以内。由于采用样品气中n-C5架桥定量时,不需要使用气体标准物质。因此,采用样品气中n-C5架桥定量分析更为方便适用。

3.2.2 烃露点计算值与直接测定值的一致性比较

现场实验时,采用增压的方式将天然气样品的压力升高到10 MPa左右,然后采用GB/T 27895-2011规定的冷却镜面目测法直接测定不同压力下的烃露点,表 5为直接测定的两个天然气样品烃露点。

表 5    直接测定两个天然气样品的烃露点 Table 5    Direct determination results of hydrocarbon dew point of two natural gas samples

在现场测定增压后样品气烃露点的同时,采用1 L钢瓶取增压后的样品,在实验室分析其组成,然后通过组成数据,采用SRK状态方程预测天然气的相包络线,并与直接测定烃露点结果绘制相包络线比较,见图 2图 3

图 2     1号样品直接测定烃露点相态曲线和预测相态曲线的比较 Figure 2     Comparison of phase behavior curves between direct determination and predictive hydrocarbon dew point for sample No.1

图 3     2号样品直接测定烃露点相态曲线和预测相态曲线的比较 Figure 3     Comparison of phase behavior curves between direct determination and predictive hydrocarbon dew point for sample No.2

结果表明:C6+组成采用特殊组分单独定量与只按碳数归类数据处理方式获得的组成结果,对烃露点预测值有较大影响。C6+中特殊组分单独定量分析时,计算的烃露点比只按碳数归类时计算的烃露点低,最大偏低值达到6 ℃左右;而采用特殊组分单独定量时,计算预测值与实际测试值更为接近,两者平均绝对偏差在3.0 ℃以内,基本能够满足生产上对预测烃露点准确度的要求。

4 结语

(1) 天然气中C6+延伸组成分析主要有两种定量方法,即采用标准气中n-C4架桥定量和样品气中n-C5架桥定量,两种方法获得的组成结果对计算的烃露点是一致的,说明选用两种定量方法均可。但在实际使用中,由于不需要使用气体标准物质,因此采用样品气中n-C5架桥定量分析更为方便适用。

(2) 采用不同数据处理方式获得组成数据对计算的烃露点有较大影响,实际应用中,C6+组成采用苯、甲苯、环己烷和甲基环己烷等组分单独定量,而其余组分按照碳数归类分析时,获得的完整组成数据计算的烃露点和直接测定结果更为接近,两者平均偏差在3 ℃以内。因此,在需要使用组成数据计算烃露点时,宜采用苯、甲苯、环己烷和甲基环己烷等组分单独定量,而其余组分按照碳数归类的分析和数据处理方式,由此获得的烃露点更为准确可靠。

参考文献
[1]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气词汇: GB/T 20604-2006[S].北京: 中国标准出版社, 2007.
[2]
曾文平, 熊钢. 计算法获得天然气烃露点影响因素探讨[J]. 石油与天然气化工, 2011, 40(5): 510-513. DOI:10.3969/j.issn.1007-3426.2011.05.020
[3]
周理, 张镨, 蔡黎, 等. 天然气烃露点预测研究进展[J]. 石油与天然气化工, 2017, 46(4): 87-92. DOI:10.3969/j.issn.1007-3426.2017.04.017
[4]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气取样导则: GB/T 13609-2017[S].北京: 中国标准出版社, 2017.
[5]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气的组成分析气相色谱法: GB/T 13610-2014[S].北京: 中国标准出版社, 2015.
[6]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气在一定的不确定度下用气相色谱法测定组成第1部分: 分析导则: GB/T 27894.1-2011[S].北京: 中国标准出版社, 2012.
[7]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气在一定不确定度下用气相色谱法测定组成第2部分: 测量系统的特性和数理统计: GB/T 27894.2-2011[S].北京: 中国标准出版社, 2012.
[8]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气在一定不确定度下用气相色谱法测定组成第3部分: 用两根填充柱测定氢、氦、氧、氮、二氧化碳和直至C8的烃类: GB/T 27894.3-2011[S].北京: 中国标准出版社, 2012.
[9]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气在一定不确定度下用气相色谱法测定组成第4部分: 实验室和在线测量系统中用两根色谱柱测定氮、二氧化碳和C1至C5及C6+的烃类: GB/T 27894.4-2012[S].北京: 中国标准出版社, 2013.
[10]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气在一定不确定度下用气相色谱法测定组成第5部分: 实验室和在线工艺系统中用三根色谱柱测定氮、二氧化碳和C1至C5及C6+的烃类: GB/T 27894.5-2012[S].北京: 中国标准出版社, 2013.
[11]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气在一定不确定度下用气相色谱法测定组成第6部分: 用三根毛细管色谱柱测定氢、氦、氧、氮、二氧化碳和C1至C8的烃类: GB/T 27894.6-2012[S].北京: 中国标准出版社, 2013.
[12]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气中丁烷至十六烷烃类的测定气相色谱法: GB/T 17281-2016[S].北京: 中国标准出版社, 2016.
[13]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.气体分析词汇: GB/T 14850-2008[S].北京: 中国标准出版社, 2008.
[14]
张娅娜, 陈赓良, 曾文平, 等. 用气相色谱法测定天然气组成国内外标准对比研究[J]. 石油与天然气化工, 2007, 36(1): 74-78. DOI:10.3969/j.issn.1007-3426.2007.01.019
[15]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气烃露点计算的气相色谱分析要求: GB/T 30492-2014[S].北京: 中国标准出版社, 2014.
[16]
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.天然气烃露点的测定冷却镜面目测法: GB/T 27895-2011[S].北京: 中国标准出版社, 2012.