天然气压缩因子的计算

合集下载

AGA8—92DC计算方法天然气压缩因子计算(最漂亮的)

AGA8—92DC计算方法天然气压缩因子计算(最漂亮的)

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。

75立方压缩天然气体积计算公式(一)

75立方压缩天然气体积计算公式(一)

75立方压缩天然气体积计算公式(一)
计算天然气体积的公式
在石油和天然气行业,计算天然气体积是一个重要的任务。

其中一个常用的计算公式是75立方压缩天然气体积计算公式。

这个公式可以帮助我们计算压缩天然气的体积。

公式
75立方压缩天然气体积计算公式如下:
V = (75 x P) / Z
其中, V 是天然气的体积(单位:立方米) P 是天然气的压力(单位:兆帕) Z 是压缩因子(无单位)
举例说明
让我们通过一个例子来说明这个计算公式。

假设我们有一个压力为 10 兆帕(MPa)的天然气样本,而其压缩因子为。

我们想要计算这个天然气样本的体积。

根据公式:
V = (75 x P) / Z
代入数值:
V = (75 x 10) / V = 立方米
因此,该天然气样本的体积为立方米。

这个例子展示了如何使用75立方压缩天然气体积计算公式来计算天然气的体积。

只需要知道压力和压缩因子,就可以轻松地计算出天然气的体积。

总结起来,75立方压缩天然气体积计算公式是一个在石油和天然气行业常用的公式,它可以帮助我们通过压力和压缩因子来计算天然气的体积。

压缩因子计算

压缩因子计算

天然气压缩因子的计算气田上大多数在高压下生产,为控制其流动需要安装节流阀。

当气流经过节流阀时,气体产生膨胀,其温度降低。

如果气体温度变得足够低,将形成水合物(一种固体结晶状的冰雪物质)。

这就会导致管道和设备的堵塞。

【1】从而,在天然气的集输过程当中,不管对天然气或天然气管道进行怎样的处理,都离不开气体的三个状态参数:压力P 、体积V、温度T。

而根据真实气体状态方程PV ZnRT =可知,在确定某个状态参数的时候需要先计算一个压缩因子Z。

如果能够更精确的确定压缩因子,从而确定气体的状态参数,对于研究天然气的收集、预处理和输送等问题具有重要意义。

下面简要介绍下压缩因子及其计算方法。

真实气体是实实在在的气体,它是为了区别于理想气体而引人的。

真实气体占有一定空间,分子之间存在作用力,因此真实气体性质与理想气体性质就有偏离。

压缩因子就是反映这种真实气体对理想气体的偏离程度大小。

在温度比临界温度高的多、压力很小时,偏离不太显著;反之偏离就很显著。

下面将介绍一种计算压缩因子的方法(Dranchuk-Purvis-Robinson 法)。

压缩因子的关系式如下:5635214373831()()()(1)exp()pr pr pr pr pr A A A A A Z A A T T T T A A A T =++++++++-52pr pr pr 222prpr pr ρρρρρρ (1)式中A 1到A 8都是常数,具体数据可到参考文献上查阅,ρpr 为无因次拟对比密度,它和压缩因子满足关系式: 0.27prpr pr p ZT ρ= (2)其中p pr 和T pr 分别为拟对比压力和拟对比温度。

由于式(2)为非线性方程,欲计算Z ,可采用牛顿迭代法(Newton-Raphson )。

在已知p pr 和T pr 的情况下,需经过迭代过程求解ρpr ,其公式如下:()(1)()'()()()i pr i i pr pr i pr f f ρρρρ+=- (3)迭代求得拟对比密度ρpr ,即可易求得压缩因子。

对aga8-92dc方程计算天然气压缩因子的探讨

对aga8-92dc方程计算天然气压缩因子的探讨

对aga8-92dc方程计算天然气压缩因子的探讨
对于aga8-92dc方程计算天然气压缩因子的探讨,涉及两个主要的概念,一个是aga8-92dc方程,另一个是压缩因子。

aga8-92dc方程是一种多参数热力学模型,它提供了精确地描述天然气性质的方法和参数。

它可以进一步详细地描述天然气中各种物理参数,如温度、压强、密度等,以及天然气的组分组成,并使用七个物理参数对天然气混合进行深入分析。

压缩因子则是指在维持温度和体积恒定的情况下,压缩比的大小。

使用aga8-92dc方程,可以根据天然气各种物理参数和组分组成,推算出不同的压缩因子,从而更好地描述天然气的变化情况。

天然气压缩因子计算

天然气压缩因子计算

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。

天然气与CNG(压缩天然气)之间的换算比例

天然气与CNG(压缩天然气)之间的换算比例

天然气车的计算方法——CNG车气瓶加气量A(m3)= L×N×P/1000其中:L(升)=气瓶标定容积,N(个)=气瓶个数,P (大气压)=气瓶内气体压强(一般为20MPa,200个大气压);气瓶加气量:G(Kg)= A×ρ,其中:ρ-天然气密度(0.716 Kg /m3)。

例1)车辆CNG气瓶为8×120 + 4×80L。

则加气量=(8×120 + 4×80L)×200/1000 = 256(立方米)= 256×0.716 = 183.3(kg)——LNG车气瓶加气量A(m3)= L×N×600/1000例2)车辆LNG气瓶为450L。

则加气量= 450 ×600/1000 = 270(立方米)= 270 × 0.716 = 193.3简易算法:CNG加气量(立方米)=气瓶总容积×0.2;LNG加气量(立方米)=气瓶总容积×0.6——CNG车在驾驶室的仪表板上有剩余气量显示系统——显示气瓶内的剩余气体压力为多少MPa,根据气瓶加气量计算方法:加气量(m3)=气瓶总容积L×气瓶个数N ×气瓶内气体压强P/ 1000,从而算出气瓶内还剩多少气量,提示司机及时加气。

例1)车辆CNG气瓶为8×120 + 4×80L,驾驶室仪表板上显示气瓶内剩余气体压力为5MPa。

则剩余气量=(8×120 + 4×80L)×50/1000 = 64(立方米)。

一般来说,当CNG车驾驶室仪表板上显示气瓶内的剩余气体压力为5MPa时,司机就应及时给车辆加气;当剩余气体压力为2MPa时,司机就必须给车辆加气。

定义:是指以天然气为燃料的一种气体燃料汽车(有压缩天然气CNG和液化天然气LNG)。

分类:按燃料使用状况的不同,可分为:(1)单燃料天然气汽车:发动机只使用CNG或LNG作为燃料。

Matlab编程天然气压缩因子计算模型

Matlab编程天然气压缩因子计算模型

1程序目的利用AGA8-92DC模型计算天然气的压缩因子,该程序主要应用于在输气和配气正常进行的压力P和温度T范围内的管输气的压缩因子计算2数学模型:AGA8-92DC模型2.1模型介绍此模型是已知气体详细的摩尔分数组成和相关压力、温度来计算气体压缩因子。

输入变量包括绝对压力、热力学温度和摩尔组成。

摩尔组成是以摩尔分数表示下列组分:CO2、N2、H2、CO、CH4、C2H6、C3H8、i-C4H10、n-C4H10、i-C5H12、n-C5H12、n-C6H14、n-C7H16、n-C8H18。

2.2 模型适用条件绝对压力:0MPa<P<12MPa热力学温度:263K≤T≤338K高位发热量:30MJ·m-3≤HS≤45 MJ·m-3 相对密度:0.55≤d≤0.80天然气中各组分的摩尔分数应在以下范围内:CH4:0.7≤xCH4≤1.0N2:0≤xN2≤0.20CO2:0≤xCO2≤0.20C2H6:0≤xC2H6≤0.10C3H8:0≤xC3H8≤0.035C4H10:0≤xC4H10≤0.015C5H12:0≤xC5H12≤0.005C6H14:0≤xC6H14≤0.001C7H16:0≤xC7H16≤0.0005C8H18和更高碳数烃类:C8H18:0≤xC8H18≤0.0005H2:0≤xH2≤0.10CO :0≤x CO ≤0.03如果已知体积分数组成,则应将其换算成摩尔分数组成。

所有摩尔分数大于0.00005的组分都不可忽略。

2.3 模型描述2.3.1 已知条件绝对压力P 、热力学温度T 、组分数N ; 各组分的摩尔分数X i ,i = 1~N ; 查附表1、2、3得到的以下数据:58种物质的状态方程参数a n ,b n , c n ,k n ,u n ,g n ,q n ,f n ,s n ,w n ; 14种识别组分的特征参数M i ,E i ,K i ,G i ,Q i ,F i ,S i ,W i ;14种识别组分的二元交互作用参数E ij ∗,U ij ,K ij ,G ij ∗。

对aga8-92dc方程计算天然气压缩因子的探讨

对aga8-92dc方程计算天然气压缩因子的探讨

对aga8-92dc方程计算天然气压缩因子的探

aga8-92dc是一个理想气体压缩因子模型,用于描述天然气的状
态变化,它可以根据温度和压强来确定天然气的压缩因子。

它由一个
三重导数加上近似调和项构成,把古典的Van der Waals方程的系数
进行了改进,而且它的表达式比Van der Waals方程更小,不仅简单
易懂,而且几乎没有实验误差。

在aga8-92dc模型中,其基本原理是
经过多次压缩后,介绍天然气压缩因子和温度、压强之间的变化关系,从而对边界值中的压缩指数、比容指数和热压缩系数进行逆推,建立
相应的函数关系。

由此,利用aga8-92dc模型可以计算天然气压缩因子,应用于天然气田开采地质、采收率计算及天然气运输管网分析等
领域。

天然气压缩因子的计算

天然气压缩因子的计算

天然气压缩因子的计算天然气的压缩系数计算方法可采用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 —天然气在操作条件下的压缩因子。

天然气压缩系数计算

天然气压缩系数计算

天然气压缩系数计算天然气压缩系数是一个重要的物理参数,它描述了天然气在压缩过程中体积变化的程度。

在石油天然气的开采、储存和运输过程中,了解和计算天然气的压缩系数对于技术人员具有重要的指导意义。

首先,我们来了解一下什么是天然气的压缩系数。

天然气是一种可燃气体混合物,主要成分是甲烷(CH4)。

在常规条件下,天然气的体积与压力成反比,即压力越高,体积越小。

天然气压缩系数则是用来描述这种压力和体积之间的关系的物理量,通常用字母Z表示。

天然气的压缩系数与温度和压力有关。

随着温度的升高,压缩系数会下降;随着压力的升高,压缩系数会增加。

压缩系数的计算公式如下所示:Z = PV / RT其中,Z表示压缩系数,P表示压力,V表示体积,R表示气体常数,T表示温度。

通过测量或计算压缩系数,我们可以进一步了解天然气在不同条件下的体积变化情况。

天然气压缩系数的计算对于天然气工程领域具有重要的指导意义。

首先,它可以帮助工程师们了解天然气储量的变化情况。

通过测量地下储气库中天然气的压缩系数,可以推算出储气库中存储的天然气量。

这对于储气库的规划和管理至关重要。

其次,天然气压缩系数的计算对于天然气的运输也非常重要。

在长距离管道输送过程中,天然气会经历压缩和脉动,对管道的设计和管道压力的控制提出了要求。

根据天然气的压缩系数,工程师可以确定所需的管道直径和压力设定,以确保天然气的安全运输。

此外,天然气压缩系数的计算还可以用于天然气的流量计算。

在煤层气开采和油气田开发过程中,天然气的流量是一个重要的指标。

通过测量天然气的压缩系数,可以准确计算出天然气的流量,为天然气的评估和开发提供依据。

综上所述,天然气压缩系数的计算对于天然气工程领域具有重要的指导意义。

它不仅可以帮助我们了解天然气在不同条件下体积的变化情况,还可以在天然气储存、运输和开发等方面提供科学的依据。

通过进一步研究和应用天然气压缩系数,我们可以更好地利用和管理天然气资源,促进能源的可持续发展。

天然气压缩因子计算及影响因素分析

天然气压缩因子计算及影响因素分析

天然气压缩因子计算及影响因素分析王春生;徐玉建;田明磊;董国庆;徐畅;陈钊【摘要】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.%天然气长输管道首端与末端之间往往会出现输差,输差是影响输气成本的一个最关键的因素。

用于计量的天然气压缩因子计算方法比较[1]

用于计量的天然气压缩因子计算方法比较[1]

N G03 0. 25 0. 60 96. 50 1. 75 0. 40 0. 10 0. 10 0. 10 0. 10 0. 10
N G04 0. 56 0. 52 94. 53 0. 96 1. 55 0. 30 0. 79 0. 22 0. 19 0. 24 0. 14
N G05 N G062)
第 20 卷第 5 期 天 然 气 工 业 集输工程
用于计量的天然气压缩因子计算方法比较
张 福 元3
(西南油气田公司天然气计量检测中心)
张福元. 用于计量的天然气压缩因子计算方法比较. 天然气工业 ,2000 ;20 (5) :73~76 摘 要 天然气压缩因子或超压因子计算结果的准确性直接影响天然气流量计量的准确性 。当前国内天然 气计量界广泛使用 A GANX219〔1〕,A GA8 号报告 , ISO 1221321997 三种天然气压缩因子计算方法标准 。文章研究了 这三种天然气压缩因子计算方法标准 ,并编写了 N GZCWIN 天然气压缩因子计算软件 ,通过对不同气样和不同温 度 、压力条件的计算 ,比较了三种计算方法的差别 ,并对这些计算方法的应用范围和不确定度提出了看法 。 主题词 天然气 计量 压缩系数 计算 方法 分析
表 2 A GA8 号报告和 ISO 12213 的适用范围
项 目
A GA8 号报告
ISO 12213
管输范围 扩展范围 管输范围 扩展范围
压 力 (MPa) 0~12 0~280 0~12 0~651)
温 度 ( ℃) - 8~65
相对密度
0. 554 ~0. 87
高位发热量 18. 7 (MJ/ m3) ~45. 1
天然气压缩因子计算方法简介
1. A GA8 号报告 在 A GA8 号报告 1994 年版中〔2〕,提供了以组成

Matlab编程天然气压缩因子计算模型

Matlab编程天然气压缩因子计算模型

1程序目的利用AGA8-92DC模型计算天然气的压缩因子,该程序主要应用于在输气和配气正常进行的压力P和温度T范围内的管输气的压缩因子计算2数学模型:AGA8-92DC模型2.1模型介绍此模型是已知气体详细的摩尔分数组成和相关压力、温度来计算气体压缩因子。

输入变量包括绝对压力、热力学温度和摩尔组成。

摩尔组成是以摩尔分数表示下列组分:CO2、N2、H2、CO、CH4、C2H6、C3H8、i-C4H10、n-C4H10、i-C5H12、n-C5H12、n-C6H14、n-C7H16、n-C8H18。

2.2 模型适用条件绝对压力:0MPa<P<12MPa热力学温度:263K≤T≤338K高位发热量:30MJ·m-3≤HS≤45 MJ·m-3 相对密度:0.55≤d≤0.80天然气中各组分的摩尔分数应在以下范围内:CH4:0.7≤xCH4≤1.0N2:0≤xN2≤0.20CO2:0≤xCO2≤0.20C2H6:0≤xC2H6≤0.10C3H8:0≤xC3H8≤0.035C4H10:0≤xC4H10≤0.015C5H12:0≤xC5H12≤0.005C6H14:0≤xC6H14≤0.001C7H16:0≤xC7H16≤0.0005C8H18和更高碳数烃类:C8H18:0≤xC8H18≤0.0005H2:0≤xH2≤0.10CO :0≤x CO ≤0.03如果已知体积分数组成,则应将其换算成摩尔分数组成。

所有摩尔分数大于0.00005的组分都不可忽略。

2.3 模型描述2.3.1 已知条件绝对压力P 、热力学温度T 、组分数N ; 各组分的摩尔分数X i ,i = 1~N ; 查附表1、2、3得到的以下数据:58种物质的状态方程参数a n ,b n , c n ,k n ,u n ,g n ,q n ,f n ,s n ,w n ; 14种识别组分的特征参数M i ,E i ,K i ,G i ,Q i ,F i ,S i ,W i ;14种识别组分的二元交互作用参数E ij ∗,U ij ,K ij ,G ij ∗。

天然气压缩因子计算及影响因素分析

天然气压缩因子计算及影响因素分析
c o mp r e s s i b i l i t y f a c t o r S O t h a t t h e p h e n o me n o n c a n b e we l l c o n t r o 1 .0n t h e b a s i s o f BW RS e q u a t i o n ,f i r s t e q u a t i o n c o e f ic f i e n t s we r e o b t a i n b y E x c e l , t h e n t h e g a s d e n s i t y wa s c a l c u l a t e d wi t h t h e s e c o e f i f c i e n t s , i f n a l l y a l l t h e s e r e s u l t s we r e p u t i n t o t h e g a s s t a t e e q u a t i o n t o o b t a i n t h e c o mp r e s s i b i l i t y f a c t o r . By s o l v i n g t h e g a s c o mp r e s s i b i l i t y f a c t o r , i t s ma i n i n f l u e n c i n g f a c t o r s we r e d e t e r mi n e d .wh i c h c o u l d h e l p t o c o r r e c t t h e t h r o u g h p u t o f n a t u r a l g a s t o k e e p
W A NG C h u n — s h e n g ,XVY u - j i a n ,T I A N Mi n g — l e i ,DO NG G u o — q i n g ,X VC h a n g ,C HE NZ h a o

压缩因子计算方法

压缩因子计算方法
压缩因子计算方法
汇报人:齐少鹏
压缩因子简介
• 由于理想气体作了两个近似:忽略气体分子本身的体积和分 子间的相互作用力,所以实际气体都会偏离理想气体。
• 压缩因子Z被引用来修正理想气体状态方程:PV=nRT 。
• 压缩因子的定义式为:Z=PV/nRT ,压缩因子的量纲为一。
• 很显然,Z的大小反映出真实气体对理想气体的偏差程度, 即Z等Vm(真实)除以Vm(理想)。由于Z反映出真实气 体压缩的难易程度,所以将它称为压缩因子。
Tp3r
(1.18 2.82 ) Tpr
r
• 特殊定义的对比密度:
r

0.06125
p pr ZTpr
exp1.2(1

1 Tpr
)2


• 方法的适用范围:Tpr>1
AGA公式

美国加利福尼亚天然气协会(CNGA)公式 前苏联气体研究所公式
低压下压缩因子的确定
P<35MPa
• Z=1+(0.31506-1.0467/Tpr-0.5783/T3pr)ρ +pr
(0.5353-0.6123/Tpr)ρ
2pr+0.6815ρ
/T 2pr
3pr

ρ
=0.27
p /ZT pr
pr
• • 迭代
RK公式
RK(Redlich—Kwong)方 程 是 1949 年 提 出 的 二 参数状态方程,它在范德瓦尔斯基础上引入了温度 对引力的修正,多用于计算干气的压缩因子,不适于 计算非极性分子的压缩因子
实验方法求取天然气压缩因子
实验测定天然气压缩因子方法是将一定质量的天然气 样品装入高压物性实验装置的PVT筒中,在恒温条件下测定 天然气的压力与实际体积V的关系。

天然气管存量的两种计算公式

天然气管存量的两种计算公式

天然气管存量计算公式1、第一种计算公式 Q=*V*P 均/(T 均**Z)其中V 是该管段内容积(即管段管容),Z 是压缩因子,Z=1/(1+*1000000*P 均*10^^C 2/T 均^, P 均=2/3[P 1++(P 2+)2/(P 1+P 2+2*] T 均=(T 1+T 2)/2+P 1、P 2、T 1、T 2分别为管段起、终点压力和温度;C 2是天然气相对密度(注:一定周期内会有小调整)。

总管存Q n 为各分段管存的求和。

2、第二种计算公式 (1)管段管存计算公式:式中:0V ——管段在标准状态下的管存量,单位为立方米(m 3) ;1V ——管段的设计管容量,单位为立方米(m 3) ,计算公式为:4V 21Ld ⨯⨯=π式中:π=;d ——管段的内直径,单位为米(m ); 10001pj pj V P T Z V P T Z ⨯⨯⨯=⨯⨯L ——管段的长度,单位为米(m ); pj P ——管段内气体平均压力(绝对压力),单位为兆帕(MPa );0T ——标准参比条件的温度,数值为; 0Z ——标准参比条件下的压缩因子,数值为; 0P ——标准参比条件的压力,数值为;pj T ——管段内气体平均温度,单位为开尔文(K );1Z ——工况条件下的压缩因子,根据GB/T 《天然气压缩因子的计算 第2部分:用摩尔组成进行计算》计算求得。

(2) 平均压力计算公式:式中:1P ——管段起点气体压力,单位为兆帕(MPa);2P ——管段终点气体压力,单位为兆帕(MPa)。

(3) 平均温度计算公式:123132T T T pj ⨯+⨯=式中:1T ——管段起点气体温度,单位为开尔文(K );2T ——管段终点气体温度,单位为开尔文(K )。

注:气体体积的标准参比条件是p 0=,T 0=12121223pj P P P P P P P ⎡⎤⨯=⨯+-⎢⎥+⎣⎦。

天然气基本压缩因子计算方法(最全)word资料

天然气基本压缩因子计算方法(最全)word资料

天然气基本压缩因子计算方法编译:阙洪培(西南石油大学审校:刘廷元这篇文章提出一个简便展开算法:任一压力-温度的基本压缩因子的输气监测计算。

这个算法中的二次维里系数来源于参考文献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总因子值,而未见有其它基本温度条件的总因子值。

天然气管存量的两种计算公式

天然气管存量的两种计算公式

天然气管存量计算公式1、第一种计算公式Q=293.15*V*P 均/(T 均*0.101325*Z)其中V是该管段内容积(即管段管容),Z是压缩因子,Z=1/(1+5.072*1000000*P 均*10A1.785A C 2/T 均A3.825),P 均=2/3[P 1+0.101325+ (P2+0.101325)2/(P 1+P2+2*0.101325)]T 均=(「+T2)/2+273.15P 1、P2、「、T2分别为管段起、终点压力和温度;G是天然气相对密度(注:一定周期内会有小调整)。

总管存Q为各分段管存的求和。

2、第二种计算公式(1)管段管存计算公式:V0二0沪丁刖乙式中:V。

管段在标准状态下的管存量,单位为立方米(m i);V1 —管段的设计管容量,单位为立方米(mb,计2算公式为:V厂亠^4式中:=3.1415926 ;d――管段的内直径,单位为米(m;L—管段的长度,单位为米(m;P pj ――管段内气体平均压力(绝对压力),单位为兆帕(MPa;T o ――标准参比条件的温度,数值为293.15K ;Z o――标准参比条件下的压缩因子,数值为0.9980 ;P o ――标准参比条件的压力,数值为0.101325MPa;T pj ――管段内气体平均温度,单位为开尔文(K);Z i ――工况条件下的压缩因子,根据GB/T 17747.2《天然气压缩因子的计算第2部分:用摩尔组成进行计算》计算求得(2)平均压力计算公式:式中:P i――管段起点气体压力,单位为兆帕(MPa);P2 ------------- 管段终点气体压力,单位为兆帕(MPa)(3)平均温度计算公式:T1 ――管段起点气体温度,单位为开尔文(K);T2 ――管段终点气体温度,单位为开尔文(K)。

注:气体体积的标准参比条件是P o =0.101325MPa T o =293.15KT1式中:3欢迎您的下载,资料仅供参考!致力为企业和个人提供合同协议,策划案计划书,学习资料等等打造全网一站式需求。

天然气压缩因子的计算

天然气压缩因子的计算

天然气压缩因子的计算
唐蒙
【期刊名称】《石油与天然气化工》
【年(卷),期】1987(016)001
【摘要】天然气压缩因子的计算,过去大多使用Standing-Katz图,此法比较费时,本文介绍了速度快,精度高的DPR法通过PC-1500电子计算机在计算天然气及含水量酸性气体的气体混合物压缩因子时的应用,文章认为,DPR法在对比温度Tr≥1.2时优于Standing-katz图,对于含酸性气体的气体混合物,将Wichert的修正原则在38-39℃和10-112atm的范围内与DPR法结合使用,效果较好。

【总页数】9页(P66-74)
【作者】唐蒙
【作者单位】四川石油局天然气研究所
【正文语种】中文
【中图分类】TE64
【相关文献】
1.AGA8-92DC计算方法天然气压缩因子计算 [J], 李丹华;姜东琪
2.基于AGA8-92DC模型的天然气压缩因子计算软件开发与应用 [J], 何冬辉
3.天然气压缩因子计算的仿真和实施 [J], 吕秦;杨惠德
4.GERG-2008方程计算天然气压缩因子及烃露点 [J], 张镨;周理;鲁春;何力;王晓莉;
赵红
5.天然气压缩因子计算的仿真和实施 [J], 吕秦;杨惠德
因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

114
CH+ CH+CO 115
y=1.20 7.36748×10-3
-2.76578×10-5
CH+ N2+ N2
122
y=0.92+0.0013(T-270)
CH+ N2+ CO2 CO2 133
N2+ N2+ CO2
223
N2+ CO2+ CO2 233
y=0.92 5.52066×10-3 3.58783×10-3
Z
= 1+
Bρ m
+

2 m
(1)
ρm = p /(ZRT )
(2)
式中有关符号表示见本文后述的符号说明。天然气压缩因子 Z 的值由方程(1)、(2)联解
求得
(1)式中:天然气第二维利系数 B 由方程(B1)求得
B = x12 B11 + 2x1 x2 B 12 + 2x1 x3 B13 + 2x1 x4 B14 + 2x1 x5 B15 + x22 B22 + 2x2 x3 B 23 +
2x2 x4 B24
+
x32 B33
+
x
2 4
B44
+
x52 B55
(B1)
(B1)式中:x1 = xCH
(B2)
x2 = xN2 x3 = xCO2 x4 = xH2 x5 = xCO
(B3) (B4) (B5) (B6)
B11 = bH 0 (0) + bH0 (1)T + bH 0 (2)T 2 + [bH1 (0) + bH1 (1)T + bH1 (2) T 2 ]H CH
天然气压缩因子的计算
(C4)
(C5) (C6) (C7) (C8)
2.用 AGA NX-19 公式计算天然气压缩因子的方法
天然气超压缩系数 Fz 是因天然气特性偏离理想气体定律而导出的修正系数,其定义为
Fz = Zn ………………………………………………………(3) Z1
式中: Zn — 天然气在标准参比条件下的压缩因子; Z1 — 天然气在操作条件下的压缩因子。
N2
222
7.84980×10-3
-3.98950×10-5
CO2
333
2.05130×10-3
3.48880×10-5
H2
444
1.04711×10-3
-3.64887×10-8
CH+ CH+ N2
112
y=0.92+0.0013(T-270)
CH+ CH+ CO2 113
y=0.92
CH+ CH+ H2
-10
1.0009 1.0041 1.0074 1.0107 1.0141 1.0312 1.0492 1.0679 1.0875 1.1079 1.1293 1.1517 1.1750 1.1993 1.2246 1.2507 1.2777 1.3054 1.3335 1.3619
-5
1.0007 1.0038 1.0069 1.0100 1.0131 1.0292 1.0459 1.0633 1.0814 1.1002 1.1198 1.1401 1.1612 1.1831 1.2056 1.2288 1.2526 1.2769 1.3014 1.3260
b=
9n-2mn3 54mH 3

E 2mH2
当 1.4>τ≥1.09 且 P1(绝)≤13.79MPa 时 E=1-0.00075H2.3e-20(τ-1.09)-0.0011(τ-1.09)0.5H2[2.17+1.4(τ-1.09)0.5-H]2
当 0.88≤τ<1.09 且 13.79MPa ≥P1 (绝)>8.963 MPa 时, E=1-0.00075H2.3(2-e-20( 1.09-τ))+0.455[200(1.09-τ)6-0.03249(1.09-τ)+2.0167(1.09-τ)2
SGERG-88 或 AGA 8-92DC 方程。
本文描述 AGA NX-19 和 SGERG-88 两种计算方法。
1.用物性值进行计算天然气压缩因子的公式
本计算公式参照国家标准 GB/T 17747 中 SGERG-88 公式,该计算公式使用高位发热量、
相对密度和 CO2 含量作为输入变量。 在 GB/T 17747 中,用物性值计算天然气压缩因子公式如下:
0
1.0005 1.0034 1.0063 1.0093 1.0122 1.0273 1.0429 1.0591 1.0758 1.0932 1.1111 1.1297 1.1488 1.1685 1.1887 1.2094 1.2305 1.2519 1.2735 1.2951
5
1.0004 1.0031 1.0058 1.0086 1.0113 1.0255 1.0401 1.0551 1.0707 1.0867 1.1031 1.1201 1.1375 1.1553 1.1735 1.1921 1.2109 1.2299 1.2490 1.2680
10
1.0002 1.0028 1.0053 1.0079 1.0105 1.0238 1.0374 1.0515 1.0659 1.0806 1.0958 1.1113 1.1272 1.1434 1.1599 1.1765 1.1934 1.2104 1.2274 1.2442
附表 1:用 SGERG-88 计算的压缩因子表格,表内数据按天然气相对密度为 0.581,二氧化 碳摩尔分数为 0.006,氢气摩尔分数为 0,高位发热量为 40.66MJ·m—3 计算。
温度℃ Zn/Zg 绝对压力(MPa)
0.10 0.20 0.30 0.40 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 5.50 6.00 6.50 7.00 7.50 8.00
天然气压缩因子的计算
天然气压缩因子的计算
天然气的压缩系数计算方法可采用 GB/T 17747-1999《天然气压缩因子的计算》,或 AGA
NX-19 方程。当为非贸易计量场合和贸易计量中符合 GB/T 18603-2001《天然气计量系统技
术要求》表 A1 准确度为 C 级要求的计量装置可考虑使用 AGA NX-19 方程,其它应采用
式中符号: H——摩尔发热量,单位:MJ·kmol-1
x ——组分的摩尔分数
CH——等价烃类 CO——一氧化碳 CO2——二氧化碳 H2——氢气 N2——氮气
ρm ——摩尔密度,单位: kmol-1·m3
p——绝对压力,单位:MPa R——摩尔气体常数,其值为 0.008314510 m3·kmol-1K-1 T——热力学温度,单位:K
(C1)
(C1)式中:
C111 = cH 0 (0) + cH 0 (1)T + cH 0 (2)T 2 + [cH1 (0) + cH1 (1)T + cH1 (2)T 2 ]H CH
+ [cH 2 (0)
+
cH 2 (1)T
+
cH 2 (2)T
2
]H
2 CH
(C2)
C222,C333,C444,C113,C223 和 C233 是温度函数的二次多项式,即:
(B8) (B9)
B13 = -0.865(B11B33 )1/2
(B10)
1
天然气压缩因子的计算
(B7)和(B8)式中维利系数温度展开式系数 b(0),b(1)和 b(2)的数值
ij
b(0)
b(1)
CH
H0
-4.25468×10-1 2.86500×10-3
CH
H1
8.77118×10-4
-5.56281×10-6
CH+ N2
12
y=0.72+1.875×10-5(320-T)2
CH+ CO2
13
y=-0.865
CH+ H2
14
-5.21280×10-2 2.71570×10-4
CH+ CO
15
-6.87290×10-2 -2.39381×10-6
N2+ CO2
23
-3.39693×10-1 1.61176×10-3
-18.028(1.09-τ)3+42.844(1.09-τ)4](H-1.3)(4.01952-H2)
当 0.84≤τ<1.09 且 P1 (绝)≤8.963 MPa 时, E=1-0.00075H2.3(2-e-20(1.09-τ))-1.317(1.09-τ)4H(1.69-H2)

156.47 Fp = 160.8-7.22Gr+Kp
m=0.0330378τ-2 –0.0221323τ-3 +0.0161353τ-5
– 0.133185τ-1 + 0.265827τ-2 +0.0457697τ-4
n=
m
H=
Pj+14.7 1000
τ=
tj+460 500
Pj = 145.04P1FP
tj =(1.8t1+492)Ft—460
D = (b + b2+B3 )1/3
Xc — 天然气中二氧化碳含量的摩尔分数(应小于 0.15),由气分析给出;
Xn — 天然气中氮气含量的摩尔分数(应小于 0.15),由气分析给出。
相关文档
最新文档