天然气压缩因子计算及影响因素分析
- 格式:doc
- 大小:931.50 KB
- 文档页数:9
天然气流量计量影响因素及对策摘要:在我国当前的能源结构中,清洁能源的占比越来越提高。
其中,最显著的是天然气能源,而且天然气已经是大多数城市使用最多的能源,并且应用在日常生活的各个方面。
近年来天然气的用户开发已经达到了饱和状态,而当前天然气公司的重点应该以提高自身运营规模为主。
而对当前天然气流量进行计量,可以为天然气公司的日常经营提供一些数据上的支持,从而来进一步的为制定政策奠定数据上的基础。
因此在本文中对当前天然气流量计量中的影响因素进行了分析,并提出了相应的改进对策,从而来提高了计量结果的准确度。
关键词:天然气,流量计量,影响因素,对策1.影响天然气体流量计计量准确性的主要因素1.1测量器件自身的影响,影响计量数据的准确对当前天然气流量计量仪表计量中发现,有50%的流量仪表在计量上的误差,还有2%的仪表时不计量的情况,还有2% 的是流量计在极小范围内不进行计量等以及3%是流量计量仪表出现故障。
再者在对流量计量仪表的维修过程中可以看出,一些轴承磨损情况非常严重、叶轮形状出现变形以及一些元器件出现老化的情况,进一步严重着影响着涡轮流量计的计量准确度。
此外,还有管道所运输的天然气中含有一些细微的杂质,进一步加快了测量仪表的老化。
更重要的是涡轮流量计在所使用的环境中,温度,湿度相对较高,同时也加快了测量仪表的老化情况。
1.2 流量计未及时修正1.2.1 流量计仪表系数自身的变化在当前大多数流量计仪表系数指的是所测量的天然气自身的体积与测量零部件所产生信号之间的对应关系,其中在气体涡轮流量计中,流量计仪表系数与叶轮和轴承的精确度有着非常大的关系。
但是在实际加工的过程中,或者存在一些误差,导致每一个流量计仪表的系数存在一些细微的偏差。
此外,经过长时间的使用,流量计仪表中的叶轮会发生变形和老化,轴承出现磨损,进一步也影响着仪表系数的准确度。
在使用一定时间后,流量计的仪表系数就会出现偏差,如果没有及时进行检查并修订仪表系数,就会导致流量计在测量时测量的误差加大。
AGA8—92DC计算方法天然气压缩因子计算摘要:按照GB/T 17747.2—1999《天然气压缩因子的计算第2部分:用摩尔组成进行计算》,采用AGA8—92DC计算方法,用VB编程计算了天然气压缩因子。
用二分法求解状态方程,精度满足工程需要。
关键词:压缩因子;AGA8—92DC计算方法;二分法1概述工作状态下的压缩因子是天然气最重要的物性参数之一,涉及到天然气的勘探、开发、输送、计量和利用等各个方面。
实测天然气压缩因子所需的仪器设备价格高,不易推广,因此计算方法发展很快,主要为经验公式和状态方程计算方法。
1992年6月26日,国际标准化组织(ISO)天然气技术委员会(TC193)及分析技术分委员会(TC193/SC1)在挪威斯泰万格(Stavanger)召开了第四次全体会议,会上推荐了两个精度较高的计算工作状态下天然气压缩因子的方程,目PAGA8-92DC方程、SGERG-88方程[1]。
随后,国际标准化组织于1994年形成了国际标准草案[2]。
AGA8-92DC方程来自美国煤气协会(AGA)。
美国煤气协会在天然气压缩因子和超压缩因子表的基础上,开展了大量研究,于1992年发表了以状态方程为基础计算压缩因子的AGA No.8报告及AGA8-92DC方程[2]。
1994年,四川石油管理局天然气研究所遵照中国石油天然气总公司技术监督局的指示,对国际标准化组织1992年挪威斯泰万格会议推荐的AGA8-92DC 方程、SGERG-88方程进行验证研究,于1996年底基本完成[2]。
1999年,四川石油管理局天然气研究院(前身为四川石油管理局天然气研究所)起草的《天然气压缩因子的计算》GB/T 17747.1~3—1999被批准、发布。
《天然气压缩因子的计算》GB/T 17747.1~3—1999包括3个部分:《天然气压缩因子的计算第1部分:导论和指南》GB/T 17747.1—1999,《天然气压缩因子的计算第2部分:用摩尔组成进行计算》GB/T 17747.2—1999,《天然气压缩因子的计算第3部分:用物性值进行计算》GB/T 17747.3—1999。
管输天然气贸易计量中压缩因子的计算肖迪;巩大利【摘要】管输天然气的贸易结算经常采用体积计量和能量计量两种方式,压缩因子作为计算参数直接影响到计量准确度.国家标准GB/T 17747提供了天然气压缩因子的两种计算方法:摩尔组成法和物性值法.目前国内管输天然气压力普遍在6 MPa 以上、12 MPa以下,在这种工况条件下,物性值法计算压缩因子与摩尔组成法计算结果偏差比较大,尤其是非烃含量高(高含N2或CO2)的气体,采用物性值法更需慎重.在管输天然气贸易计量中,应采用适用范围更广,计算精度更高的摩尔组成法;物性值法是在现场增设在线物性参数测量仪器而采用的简单方法,此方法适用于无法得到气体组成且对计量准确度要求不高的情况.【期刊名称】《油气田地面工程》【年(卷),期】2011(030)009【总页数】1页(P24)【关键词】天然气;压缩因子;摩尔组成法;物性值法【作者】肖迪;巩大利【作者单位】国家石油天然气大流量计量站;国家石油天然气大流量计量站【正文语种】中文近年来,我国天然气工业迅速发展,建设了一批管道工程项目,引进了多条跨国管道。
管输天然气的贸易结算经常采用体积计量和能量计量两种方式,压缩因子作为计算参数直接影响到计量准确度。
国家标准《天然气压缩因子的计算(GB/T17747-1999)》规定了天然气压缩因子的两种计算方法,通过对两种方法比较,可明确各自的适用范围,确保国家和企业的合法权益。
国家标准《天然气压缩因子的计算GB/T 17747)》提供了天然气的压缩因子的两种计算方法:摩尔组成法和物性值法。
摩尔组成法也叫详细特征法(源自AGA8-92DC),采用已知天然气的详细摩尔组成和相关压力、温度计算压缩因子;物性值法,又称为总体特征法(源自SGERG-88),通过获取天然气的高位发热量、相对密度、CO2含量和N2含量中任意3个变量作为输入变量的压缩因子计算方法。
利用物性值计算压缩因子时,GB/T 17747不推荐采用N2含量作为输入变量之一,只给出了前3个变量作为输入变量时的压缩因子计算方法。
1.天然气相关物性参数计算密度计算: TZR PM m =ρ ρ——气体密度,Kg/m 3;P ——压力,Pa ;M ——气体千摩尔质量,Kg/Kmol ;Z ——气体压缩因子;T ——气体温度,K ;R m ——通用气体常数,8314.4J/Kmol·K 。
2.压缩因子计算:已知天然气相对密度∆时。
96.28M =∆ M ——天然气的摩尔质量。
∆+=62.17065.94pc T510)05.493.48(⨯∆-=pc P ;pc pr P P P = pcpr T T T =; P ——工况下天然气的压力,Pa ;T ——工况下天然气的温度,k ;P Pc —临界压力;T Tc ——临界温度。
对于长距离干线输气管道,压缩因子常用以下两式计算:668.34273.01--=prpr T P Z 320107.078.068.110241.01prpr pr pr T T T P Z ++--=对于干燥天然气也可用经验公式估算: 15.1117.0100100P Z +=标况流量和工况流量转换。
为了控制Welas 的5L/min 既 0.3立方米每小时的工况流量。
Q 2------流量计需要调节的流量值P 2------0.1MpaT 2------293.15K (20℃ )Z 2------标况压缩因子Q 1------0.3m 3/hP 1------ 工况压力(绝对压力MPa )T 1------开尔文KZ 1-------工况压缩因子转换公式为12221211p T Z Q Q p T Z。
天然气压缩因子及标准体积计算天然气压缩因子及标准体积计算导语:天然气是一种重要的能源资源,广泛应用于工业、民用和交通领域。
在储存和运输过程中,天然气往往会受到压缩或膨胀的影响。
为了更准确地计量天然气的数量,我们需要了解天然气压缩因子和标准体积的概念,并掌握相关的计算方法。
一、天然气压缩因子的概念天然气的体积与压力、温度以及成分有关,而天然气的压缩因子则是描述天然气体积变化的重要参数。
压缩因子是指实际天然气体积与理论天然气体积之间的比值。
天然气在不同压力和温度下的压缩因子是不同的,通常用Z表示。
当Z=1时,说明天然气符合理想状态,即PV=ZnRT,其中P是天然气的压力,V是天然气的体积,n是物质的摩尔数,R是气体常数,T是天然气的温度。
当Z小于1时,说明天然气存在压缩,体积变小;当Z大于1时,说明天然气存在膨胀,体积变大。
二、天然气压缩因子的计算方法天然气压缩因子的值受到很多因素的影响,如压力、温度、天然气的成分和田间条件等。
常见的计算方法有实验法和经验法。
实验法是通过实验测定压缩因子的值。
实验室通常使用高精度的实验装置,将天然气在不同压力和温度下进行测量,并计算出相应的压缩因子。
这种方法准确度高,但时间成本较高,不适合大规模应用。
经验法是通过统计数据建立的数学模型来计算压缩因子。
常用的经验法有很多,如Dranchuk-Abou-Kassem (DAK)模型、Peng-Robinson (PR)模型等。
这些模型基于一定的假设和实验数据,通过计算方程以及相似物性参数,预测天然气的压缩因子。
经验法计算速度较快,适用于大规模计算,但存在一定的误差。
三、标准体积的概念和计算标准体积是指天然气在标准条件下的体积,常用单位是立方米(m³)。
标准条件一般指标准大气压(101.325千帕)和摄氏度为15℃(或20℃)的状态。
天然气的实际体积与标准体积之间存在一定的关系,可以通过压缩因子进行计算。
标准体积与实际体积之间的关系可以用以下公式表示:V_std = V_actual * Z * (P_std / P_actual) * (T_actual / T_std)其中,V_std是标准体积,V_actual是实际体积,Z是压缩因子,P_std和P_actual分别是标准压力和实际压力,T_actual和T_std分别是实际温度和标准温度。
AGA8—92DC计算方法天然气压缩因子计算摘要:按照GB/T 17747.2—1999《天然气压缩因子的计算第2部分:用摩尔组成进行计算》,采用AGA8—92DC计算方法,用VB编程计算了天然气压缩因子。
用二分法求解状态方程,精度满足工程需要。
关键词:压缩因子;AGA8—92DC计算方法;二分法1概述工作状态下的压缩因子是天然气最重要的物性参数之一,涉及到天然气的勘探、开发、输送、计量和利用等各个方面。
实测天然气压缩因子所需的仪器设备价格高,不易推广,因此计算方法发展很快,主要为经验公式和状态方程计算方法。
1992年6月26日,国际标准化组织(ISO)天然气技术委员会(TC193)及分析技术分委员会(TC193/SC1)在挪威斯泰万格(Stavanger)召开了第四次全体会议,会上推荐了两个精度较高的计算工作状态下天然气压缩因子的方程,目PAGA8-92DC方程、SGERG-88方程[1]。
随后,国际标准化组织于1994年形成了国际标准草案[2]。
AGA8-92DC方程来自美国煤气协会(AGA)。
美国煤气协会在天然气压缩因子和超压缩因子表的基础上,开展了大量研究,于1992年发表了以状态方程为基础计算压缩因子的AGA No.8报告及AGA8-92DC方程[2]。
1994年,四川石油管理局天然气研究所遵照中国石油天然气总公司技术监督局的指示,对国际标准化组织1992年挪威斯泰万格会议推荐的AGA8-92DC 方程、SGERG-88方程进行验证研究,于1996年底基本完成[2]。
1999年,四川石油管理局天然气研究院(前身为四川石油管理局天然气研究所)起草的《天然气压缩因子的计算》GB/T 17747.1~3—1999被批准、发布。
《天然气压缩因子的计算》GB/T 17747.1~3—1999包括3个部分:《天然气压缩因子的计算第1部分:导论和指南》GB/T 17747.1—1999,《天然气压缩因子的计算第2部分:用摩尔组成进行计算》GB/T 17747.2—1999,《天然气压缩因子的计算第3部分:用物性值进行计算》GB/T 17747.3—1999。
编号:天然气流量计量压缩因子影响探讨张秋萍中国石化仪征化纤股份有限公司设备管理部联系电话:*************电子邮箱:************************二○一四年五月天然气流量计量压缩因子影响探讨张秋萍(仪征化纤股份有限公司设备管理部江苏仪征211900)摘要:通过分析天然气流量计量原理,利用文献资料提供的实验数据,结合理论研究结论,得出天然气压缩因子的影响因素有:工况压力、温度及气体组份,定性并部分定量给出了这些因素对压缩因子及计量值的影响,探讨了天然气流量计量的最佳解决方案。
关键词:压缩因子工况压力工况温度气体组份计量误差To investigate the influence of compressibility factor of naturalgas flow meteringZhang quipping(Equipment management department of Yizheng Chemical Fibre Co.,Ltd.,Yizheng Jiangsu211900)Abstract:Through the analysis of natural gas flow measurement principle, the experimental data provided by the literature, combined with the theoretical research conclusions, obtained the influencing factors of natural gas compression factor: working pressure, temperature and gas composition, qualitative and quantitative given the impact of these factors on the compression factor and the value of measurement, discusses the optimal solution of natural gas flow measurement.Keywords: compression factor working pressure working temperature gas composition measurement error.一、引言目前,国内天然气作为一种重要的清洁能源越来越多地被应用,管输天然气的贸易结算大都采用换算到标准状态(1个标准大气压、20℃)下的标况体积流量。
天然气压缩因子的计算天然气的压缩系数计算方法可采用GB/T 17747-1999《天然气压缩因子的计算》,或AGA NX-19方程。
当为非贸易计量场合和贸易计量中符合GB/T 18603-2001《天然气计量系统技术要求》表A1准确度为C 级要求的计量装置可考虑使用AGA NX-19方程,其它应采用SGERG-88或AGA 8-92DC 方程。
本文描述AGA NX-19和SGERG-88两种计算方法。
1.用物性值进行计算天然气压缩因子的公式本计算公式参照国家标准GB/T 17747中SGERG-88公式,该计算公式使用高位发热量、相对密度和CO 2含量作为输入变量。
在GB/T 17747中,用物性值计算天然气压缩因子公式如下:21mm C B Z ρρ++= (1 /(ZRT p m =ρ (2式中有关符号表示见本文后述的符号说明。
天然气压缩因子Z 的值由方程(1、(2联解求得(1式中:天然气第二维利系数B 由方程(B1求得B x x B x B 2111212+=12++++++233222255254424332324422B x B x B x B x x +++ (B1(B1式中:CH x x =1 (B2 22N x x = (B3 23CO x x = (B4 24H x x = (B5 CO x x =5 (B6 2(1(0([2(1(b 0(11120H0011H H H H H b T b b T b T b B +++++=CH H T ]222222]2(1(0([CH H H H H T b T b b +++ (B7B 14,B 15,B 22,B 23,B 24,B 33, B 44和B 55是温度函数的二次多项式,即:2-12B B B +×+= (B91/23311130.865(B -B B = (B10(B7和(B8式中维利系数温度展开式系数b(0,b(1和b(2的数值 ij b(0 b(1 b(2 CH H0 -4.25468×10-1 2.86500×10-3 -4.62073×10-1CH H1 8.77118×10-4 -5.56281×10-6 8.81510×10-9 CH H2 -8.24747×10-7 4.31436×10-9 -6.08319×10- 12 N 2 22 -1.44600×10-1 7.40910×10-4-9.11950×10-7CO 2 33 -8.68340×10-1 4.03760×10-3 -5.16570×10-6H 244-1.10596×10-3 8.13385×10-5-9.87220×10-8 CO 55 -1.30820×10-1 6.02540×10-4 -6.44300×10-7CH+ N 2 12 y =0.72+1.875×10-5(320-T2 CH+ CO 2 13 y =-0.865CH+ H 2 14 -5.21280×10-2 2.71570×10-4-2.50000×10-7CH+ CO 15-6.87290×10-2-2.39381×10-6 5.18195×10-7N 2+ CO 2 23 -3.39693×10-1 1.61176×10-3 -2.04429×10-6N 2+ H 2 24 1.20000×10-2 0.00000 0.00000天然气第二维利系数C 由方程(C1求得1233211222444343333323323222332222232133231333C x C x C x x C x x C x C x x ++++++ (C1 (C1式中:211120001112(1(0([2(1(0(T c T c c T c T c c C H H H H H H +++++=CH H ]22222]2(1(0([CH H H H H T c T c c +++ (C2C 222,C 333,C 444,C 113,C 223和C 233是温度函数的二次多项式,即:22(1(0(T c T c c C ijk ijk ijk ijk ++= (C3维利系数温度展开式中系数c(0,c(1和c(2的数值ijk c(0 c(1 c(2CH H0 -3.02488×10-1 1.95861×10-3 -3.16302×10-6CH H1 6.46422×10-4-4.22876×10-6 6.88157×10-9CH H2 -3.32805×10-7 2.23160×10-9 -3.67713×10-12N 2 222 7.84980×10-3-3.98950×10-5 6.11870×10-8CO 2 333 2.05130×10-3 3.48880×10-5 -8.37030×10-8H 2 444 1.04711×10-3-3.64887×10-8 4.67095×10-9CH+ CH+ N 2 112 y =0.92+0.0013(T-270 CH+ CH+ CO 2 113 y = 0.92 CH+ CH+ H 2 114 y = 1.20 CH+ CH+CO1157.36748×10-3-2.76578×10-5 3.43051×10-8CH+ N 2+ N 2 122 y =0.92+0.0013(T-270 CH+ N 2+ CO 2 123 y =1.10 CH+ CO 2+ CO 2 133 y =0.92 N 2+ N 2+ CO 2 223 5.52066×10-3 -1.68609×10-5 1.57169×10-8N 2+ CO 2+ CO 2 2333.58783×10-3 8.06674×10-6 -3.25798×10-8其他非同类交互作用维利系由方程(C4求得: ijk ijk y C =3/1(kkk jjj iii C C C (C4(C4式中ijk y 由(C5~(C8给出:92.0133113==y y (C620.1114=y (C710.1123=y (C8式中符号:H ——摩尔发热量,单位:MJ ·kmol -1x ——组分的摩尔分数CH ——等价烃类 CO ——一氧化碳 CO 2——二氧化碳 H 2——氢气 N 2——氮气m ρ——摩尔密度,单位: kmol -1·m 3p ——绝对压力,单位:MPaR ——摩尔气体常数,其值为0.008314510 m 3·kmol -1K -1 T ——热力学温度,单位:K2.用AGA NX-19公式计算天然气压缩因子的方法天然气超压缩系数Fz 是因天然气特性偏离理想气体定律而导出的修正系数,其定义为1Z ZnFz =………………………………………………………(3 式中: Zn —天然气在标准参比条件下的压缩因子;Z1 —天然气在操作条件下的压缩因子。
天然气压缩因子计算及影响因素分析王春生;徐玉建;田明磊;董国庆;徐畅;陈钊【摘要】Measurement shortage will often arise between the head and the end of nature gas pipeline which is a vital important influencing factor of transmission cost. Regarding to the phenomenon of measurement shortage, we focused on the compressibility factor and tried to solve the problem by optimizing the calculation method of the compressibilityfactor so that the phenomenon can be well control. On the basis of BWRS equation, first equation coefficients were obtain by Excel, then the gas density was calculated with these coefficients, finally all these results were put into the gas state equation to obtain the compressibility factor. By solving the gas compressibility factor, its main influencing factors were determined, which could help to correct the throughput of natural gas to keep measurement shortage to the minimum.%天然气长输管道首端与末端之间往往会出现输差,输差是影响输气成本的一个最关键的因素。
2.1 天然气临界参数计算2.1.1 天然气平均分子量天然气是混合气体,分子量不是一成不变的,其平均分子量按Key 规则计算: g i i M y M =∑ (2.1) 式中 M g —天然气的平均分子量kg/mol ;M i 、y i —天然气中i 组分的分子量和摩尔分数。
2.1.2 天然气的相对密度首先假定空气和天然气都取同一标准状态,天然气的相对密度可用下式表示:28.9729g g g g g air air M M M r M ρρ===≈(2.2) 式中 r g —天然气的相对密度;g ρair ρ—同一标准状态下,天然气、空气的密度kg/m 3;g M air M —天然气、空气的平均分子量kg/mol 。
2.1.3 拟临界压力P PC 和拟临界温度T PC① 组分分析方法pc i ci p y p =∑pc i ci T yT =∑(2.3) g i i M y M =∑式中 ci p —— 天然气组分i 的临界压力(绝),MPa ;ci T —— 天然气组分i 的临界温度,(273+t)°K 。
② 相关经验公式方法在缺乏天然气组分分析数据的情况下,可引用Standing 在1941年发表的相关经验公式对于干气2pc 2pc 4.6660.1030.2593.31817g g g gp T γγγγ=+-=+- (2.4)对于湿气2pc 2pc 4.8680.35639.7103.9183.339.7g g g gp T γγγγ=+-=+- (2.5)也可以用下面经验关系式进行计算 对于干气pc pc pc pc 4.88150.386192.2222176.66670.74.77800.248292.2222176.66670.7g g g g gg p T p T γγγγγγ=-=+≥=-=+< (2.6)对于湿气pc pc pc pc 5.10210.6895132.2222176.66670.74.77800.2482106.1111152.22220.7g g g g gg p T p T γγγγγγ=-=+≥=-=+< (2.7)注意:上式是对于纯天然气适用,而对于含非烃CO 2 、H 2S 等可以用Wichert 和Aziz 修正。
天然气压缩因子规范
天然气压缩因子规范是指在一定温度和压力条件下,天然气体积收缩率的计算规范。
根据国家标准《液化天然气技术要求》(GB/T 19856-2005)的规定,天然气的压缩因子应满足以下要求:
1. 在温度为15℃,压力为101.325 kPa (1兆帕)时,天然气的压缩因子应大于等于1.0。
2. 在温度为15℃,压力为50.0 kPa (0.5兆帕)时,天然气的压缩因子应大于等于0.8。
3. 在温度为15℃,压力为25.0 kPa (0.25兆帕)时,天然气的压缩因子应大于等于0.6。
4. 在温度为15℃,压力为10.0 kPa (0.1兆帕)时,天然气的压缩因子应大于等于0.4。
5. 在温度为15℃,压力为5.0 kPa (0.05兆帕)时,天然气的压缩因子应大于等于0.2。
浅谈压缩因子及其在燃气计量中的影响Gas Compression Factor(气体压缩因子)摘要:本文通过对压缩因子的介绍,浅要分析了其在燃气计算中的影响。
关键词:临界温度;临界压力;势能;动能;压缩因子;燃气计量1 前言随着燃气业的不断发展,工商业燃气用户也越来越多,特别是将来天然气在广东珠三角的推广将使珠三角燃气行业得到更迅猛的发展。
由于燃气公司与供应商,燃气公司与用户之间,都是存在经济的利益关系,因而在燃气计量上必然存在着或多或少的矛盾。
特别是大型工业用户与燃气公司的计量问题的更容易发生矛盾。
为了减少此类矛盾的激化,燃气公司在自身利益不受损害的情况下,都将会尽量减少计量上的误差,使两者的矛盾得到缓和。
因而燃气公司在计量上就应了解许多注意事项,其中压缩因子就是众多影响计量准确度的潜在因素之一。
2 压缩因子的产生及在计量上的影响在实际的气体计量的过程中,气体状态方程:Z=PV/RT,压缩因子z在计算中引入了临界温度Tc和临界压力Pc两种参数,其中压缩因子Z=f(Pc,Tc)是随温度及压力而变化的(关系图l,2),其中:临界温度指的是气体加压液化所允许的最高温度,一般分子间的引力越大对应的临界温度越高;如甲烷临界温度Tc为191.05K,丙烷临界温度Tc为368.85K;临界压力指的是气体在临界温度时发生液化所需要的最小压力;如甲烷临界压力Pc为4.6407MPa,丙烷临界压力Pc为4.3975MPa;在实际计算中,还要引入对比压力Pr和对比温度Tr,所谓的对比压力就是实际工作压力和临界压力Pc的比值,对比温度同理亦是实际工作温度和临界温度Tc的比值。
压缩因子的数值在不同温度压力下的也不是完全沿一个趋势变化的,如天然气是先随压力增大而变小,到达一定程度后又逐渐随压力增大而变大。
那为什么会出现这种情况呢?从微观上讲,这主要是由于分子间的作用力造成的。
一定温度下的气体在压力较小时,分子间的距离较大,分子间一般表现为引力,造成实际气体比理想气体更易于被压缩,但随着压力增加气体分子逐渐靠近,分子间的作用力表现为排斥力造成气体难于被压缩,进而形成压缩因子的数值上的变化,并且该数值随气体种类不同而不同。
天然气基本压缩因子计算方法编译:阙洪培(西南石油大学审校:刘廷元这篇文章提出一个简便展开算法:任一压力-温度的基本压缩因子的输气监测计算。
这个算法中的二次维里系数来源于参考文献1。
计算的压缩因子接近AGA 8状态方程值[2]。
1 测量在天然气工业实用计量中,压力、温度变化作为基本(或标准条件,不仅地区间有差别,而且在天然气销售合同也有不同。
在美国,通常标准参考条件是60°F和14.73 psia。
欧洲常用的基本条件是0 ℃和101.325 kPa,而标准条件是15 ℃和101.325 kPa。
阿根廷也用15 ℃和101.325 kPa,而墨西哥则用的是20 ℃和1kg/ sq cm(绝对。
计算真实气体的热值、密度、基本密度、基本体积、以及沃贝指数时要求已知基本条件的压缩因子。
表1是理想气体值。
表1中的理想气体值不能用于密闭输气,必须计算相应基本条件的压缩因子。
计算其它基本条件的压缩因子可用AGA 8 程序,但代数计算较复杂,计算机编程共有三组软件,比较耗时。
本文提出了一个展开算法,计算密闭输气基本条件(基本条件可是任何压力温度的压缩因子。
2 压缩因子接近外界条件时,即压力小于16 psia,截断维里状态方程(方程组中的方程1较好地描述了天然气的体积性质。
方程1中,各符号的物理意义是:Z = 基本条件下压缩因子B = 二次维里系数R = 气体常数P = 基本条件的绝对压力T = 温度条件的绝对压力天然气基本压缩因子接近1,如0.99,B必然为负(图1方程2是混合物的二次维里系数,式中B ij = B ji为组分i和j的二次交互维里系数,B ii为纯组分i 的二次维里系数。
二次维里系数是温度的函数。
也可用方程3求B,便于手工计算。
比较适合密闭输气计算,方程3中B i的平方根为总因子,参见参考文献1,3,4。
问题的提出:表中常见60°F总因子值,而未见有其它基本温度条件的总因子值。
天然气压缩因子计算及影响因素分析作者:王春生徐玉建田明磊董国庆徐畅陈钊来源:《当代化工》2015年第06期摘要:天然气长输管道首端与末端之间往往会出现输差,输差是影响输气成本的一个最关键的因素。
针对出现的输差问题,以天然气组分为基础,以压缩因子作为突破口,通过着重理解天然气压缩因子的解法与改进来得到控制输差。
以BWRS方程为重点,通过Excel求得方程系数,然后从中解出气体密度,再代入气体状态方程中求得压缩因子。
通过对天然气压缩因子的求解,得到影响压缩因子的主要因素,从而修正到天然气输量,以便减少输差。
关键词:输差;压缩因子;BWRS方程;影响因素中图分类号:TQ 018 文献标识码: A 文章编号: 1671-0460(2015)06-1408-04Calculation of Natural Gas Compressibility Factor and Its Influence FactorsWANG Chun-sheng1,XV Yu-jian1,TIAN Ming-lei1,DONG Guo-qing1,XV Chang1,CHEN Zhao2(1. Northeast Petroleum University, Heilongjiang Daqing 163318,China;2. 7th China Petroleum Engineering&Construction Corp, Shandong Qingdao 266300,China)Abstract: Measurement shortage will often arise between the head and the end of nature gas pipeline which is a vital important influencing factor of transmission cost. Regarding to the phenomenon of measurement shortage, we focused on the compressibility factor and tried to solve the problem by optimizing the calculation method of the compressibility factor so that the phenomenon can be well control. On the basis of BWRS equation, first equation coefficients were obtain by Excel, then the gas density was calculated with these coefficients, finally all these results were put into the gas state equation to obtain the compressibility factor. By solving the gas compressibility factor, its main influencing factors were determined, which could help to correct the throughput of natural gas to keep measurement shortage to the minimum.Key words: Measurement shortage; Compressibility factor; BWRS equation; Influence factor天然气与其他能源材料,例如煤炭和石油相比-天然气的热值较高,利用率较高,并且对环境的污染很小。
当下,国家大力提倡使用清洁能源,天然气成为能源的主流已成为不争的事实。
然而在天然气输送过程中往往会出现输差现象,而这一现象除了降低管道输送过程中的经济效益外,还会带来一些安全隐患。
因此,如何提高对管道运输过程中的输差的控制成为当下一个重要问题。
本文以庆哈管线为研究实例,针对庆哈管线存在输差现状,开展研究。
结合现场数据,逆推相应管段管存量的计算误差,以及气体组成对压缩因子的影响。
最后综合研究成果及数据结论,了解其他因素的影响。
制定具有实际效果的的可行性解决方案,从而指导现场生产实践,最终达到显著改善庆哈输气管线输差的目标。
1 压缩因子压缩因子对天然气计量起着至关重要的作用[1]。
在日常生产实验中,我们对天然气的各项参数进行计算时,往往将其看成是理想气体来适应状态方程的需要。
但是在实际中,理想状态的气体是不存在的,因此我们便需要一个纠正系数,来对我们计算的结果进行修正,使计算结果更具有准确性更接近于实际情况,而这个修正系数,便是压缩因子Z,如下面的公式所示:(1)压缩因子在天然气的各项计算当中起着至关重要的作用,特别是对天然气的计量方面。
因此压缩因子的准确与否对天然气各项参数的计算具有决定性的作用。
常规计算压缩因子的方法有通过气体压缩因子图得出的方法、通过实验求取的方法、SK 图版法、AGA8-92DC计算法以及用BWRS方程求解压缩因子[2,3]。
2 用BWRS方程求解压缩因子2.1 BWRS方程的基本理论BWRS状态方程是一个多参数状态,其基本形式为[4]:(2)式中: P —压力,kPa;T—温度,K;ρ—密度,kmol/m3 ;Rm—通用气体常数,为8.314 3 kJ/(kmol·K)。
由式(2)可得,压缩因子为:(3)由气体状态方程p=ZρRT可得:Z=p/(ρRT)(4)上述式中,A0,B0,C0,D0,E0,a,b,c,d,ɑ,γ为状态方程得11个参数。
可以根据临界温度Tci、临界密度ρci以及偏心因子wi之间的关系式来求得纯组分的11个参数,其关系式如下所示:(5)式中Ai,Bi为通用常数(i=1,2,3,...,11)见表1。
表1 通用常数Ai和Bi的值Table 1 Universal constant value of Ai and Bi当所求目标为混合物时,BWRS采用下式进行计算:(6)式中:yi—气相或液相混合物中第i组分的摩尔分数;Kij—为第i,j组分间交互作用系数(Kij=Kji)。
Kij为和理论混合物的偏差,偏差随着Kij值的怎大而增大,对于相同组分,Kij=0。
2.2 压缩因子的求解首先,通过计算得出气体常数以及天然气分子量[5],之后运用式(5)求出流体中各个纯组分所对应的11个参数值:A0i,B0i,C0i,D0i,E0i,ai,bi,ci,di,αi,γi接着由式(6)混合物的混合规则,求得流体混合物对应的11个参数:A0,B0,C0,D0,E0,a,b,c,d,α,γ。
将已知压力p,温度T代入(3)中计算出气体的密度。
可以将式(3)改为以下形式以便于求解,求解时运用弦截法。
(7)弦截法迭代公式为:(8)即为:(9)弦截法求解需设两个初值,可设ρ1=0,ρ2=p/(RT)。
3 基本概况以庆哈输气管线的数据为例,输入气体组分,纯物质的临界温度,临界压力,偏心因子,分子量,气体通用常数Aj,Bj基础数据,输入二元交互系数Kij,计算BWRS方程中所需的各纯组分所对应的11个参数值:A0i,B0i,C0i,D0i,E0i,ai,bi,ci,di,αi,γi 如表2,其中xi为摩尔比。
计算天然气混合物对应的11个参数,即待解方程的所有参数:A0,B0,C0,D0,E0,a,b,c,d,α,γ 。
根据Excel与公式可知需先计算j变量的和,再计算i变量的和。
根据BWRS方程计算气体密度。
我们假定T=303.15 K,压力p=4 000 kPa,则BWRS方程的系数为A=2 520.389 1,B=-131.91,C=10.852,D=3.685,E=0.005 51,求得气体压缩密度为:ρ=1.714969 kmol/m3Z=p/(ρRT)=4 000/(1.714969×8.314×303.15)=0.9254现场中采用的Z=1,修正了(1-0.9254)/1×100%=7.46%根据公式(10)就可以减少管道储气量计算误差。
4 计算分析4.1 温度对压缩因子的影响从可以发现影响压缩因子的主要因素是温度,随着温度的升高,压缩因子逐渐接近于1 如图1。
图1 压缩因子随温度的变化Fig.1 Compression factor variation with temperature4.2 压力对压缩因子的影响压力与压缩因子成正比例关系如图2。
图2 压缩因子随压力变化曲线Fig.2 Compression factor variation with pressure curve4.3 天然气组成对压缩因子的影响由于天然气的主要成分是甲烷,所以主要探究甲烷含量与压缩因子的关系,如图3,甲烷含量越大,压缩因子越接近于1。
图3 压缩因子随甲烷含量变化曲线Fig.3 Compression factor variation with methane content4.4 管存量及输配率随压缩因子的变化分析假设管道直径Φ711×12.5(15),管道长度为70 km,理想管存为=3.14×[(711-25)/1000/2]2×3.14×7000/4=25 859 m3。
标准压力P0为101.325 kPa,标准温度T0为293.15 K,末端压力P按4 000 kPa,末端温度T按303.15 K。
则由公式 V储=V×P×T0/(P0×T×Z)知V×P×T0/(P0×T)为一常数,即987 169.28。
所以可以得到压缩因子对管存量的影响:采用修正后的压缩因子,管存量的到了相当大的改善,如图4。
图4 压缩因子对管存量的影响Fig.4 Influence of compression factor on the tube stock图5 压缩因子对输配率的影响Fig.5 Influence of compression factor on the distribution rate未修正的压缩因子为1,即V储=987 169.28 m3。
修正后的压缩因子为0.925 4,得到:V储=1 066 748.73 m3由于损耗量=总进气量-管存量-总外输量,而损耗率=总损耗量/总进气量×100%。
可见,管存量增大,损耗量便减小了,进而损耗率减小,又输配率=100%-损耗率,所以输配率得到提高如图5阴影部分(表2)。