天然气物性参数(新)

1、。2.1 天然气临界参数计算2.1.1天然气平均分子量天然气是混合气体，分子量不是一成不变的，其平均分子量按Key 规则计算 :M gyi M i（ 2.1 ）式中 M g天然气的平均分子量 kg/mol ；Mi、yi 天然气中 i 组分的分子量和摩尔分数。2.1.2天然气的相对密度首先假定空气和天然气都取同一标准状态，天然气的相对密度可用下式表示：gM gM gM g（2.2 ）rgM air28.9729air式中 r g天然气的相对密度；gair 同一标准状态下，天然气、空气的密度3；kg/mM g M air 天然气、空气的平均分子量kg/mol 。2.1.3拟临界压力 PPC和拟临界

2、温度 TPC 组分分析方法ppcyi pciTpcyTi ci（2.3 ）M gyi M i式中pci 天然气组分 i 的临界压力 ( 绝),MPa；Tci 天然气组分 i 的临界温度 ,(273+t)°K。 相关经验公式方法在缺乏天然气组分分析数据的情况下，可引用 Standing 在 1941 年发表的相关经验公式对于干气。1。ppc4.6660.1030.252ggTpc93.3181 g72g对于湿气ppc4.8680.35639.72ggTpc103.9183.339.72gg也可以用下面经验关系式进行计算对于干气ppc4.88150.3861 gTpc92.2222176

3、.6667ggppc4.77800.2482 gTpc92.2222176.6667gg（ 2.6 ）对于湿气（2.4 ）（2.5 ）0.70.7ppc5.10210.6895 gTpc132.2222 176.6667 gg0.7（2.7 ）ppc4.77800.2482 gTpc106.1111 152.2222 gg0.7注意：上式是对于纯天然气适用，而对于含非烃CO2 、 H2 S 等可以用 Wichert 和Aziz 修正。修正常数的计算公式为 :p4.75460.2102g0.031.158310 2N23.0612102pcCO2H 2 STpc84.9389188.4944 g

4、 0.9333 CO2 1.4944(2.8)N 22.1.4拟对比压力 PPr 和拟对比温度 TPr 的计算对比参数就是指某一参数与其应对应的临界参数之比：即pprpppcTPrT（2.9）Tpc2.2天然气的偏差因子Z 计算天然气偏差因子 Z 的计算是指在某一压力和温度条件下，同一质量气体的真实体积。2。与理想体积之比值。ZV实际 V实际V理想nRTp计算天然气偏差因子方法较多，下面主要介绍几种常用的计算方法2.2.1 Pong.Robinson方程法RTa(T )pb) b(V b)V b V (V式中nna(T )xi x j ( ai a j i j ) 0.5 (1 Kij )i

5、1j 1nbxbiij 1（ 2.13 ）aiR2Tcr20.45724pcrbi0.0788 RTcrpcr2i1mi 1 Tri0.5mi0.374641.5422i0.26992i2式中Kij 天然气的交互作用参数；pcr 组分 i 的气体临界压力；Tcr 组分 i 的气体临界温度；Tr 组分 i 的对比温度；i 组分i 的偏心因子。由方程可得到关于Z 的方程Z 3(1 B)Z 2A 3B22BZ AB B2B30aPA2 T2R（ 2.10 ）（ 2.11 ）（ 2.12 ）（ 2.14 ）（ 2.15 ）（ 2.16 ）（ 2.17 ）（ 2.18 ）。3。bPB（2.19 ）RT2

6、.2.2 Cranmer方法1.04670.5783Z 1 0.31506TprTpr3pr0.6123220.6815pr0.5353prTpr3Tpr（ 2.20 ）pr 0.27 ppr / Z Tpr（2.21 ）式中pr 拟对比密度。已知 P、T 求 Z，计算步骤如下：第一步计算 ppc ， Tpc ；第二步计算 ppr ， Tpr ；第三步 对Z赋初值，取 Z，利用式（3.12）计算pro=1第四步将pr 值代入式（ 3.11 ），计算 Z该方法适用于 p<35MP的情况。2.2.3 DPR法1974 年,Dranchuk ， Purvis和 Robinson 等人在拟合 S

7、tanding.Katz图版的基础上, 提出了计算偏差因子 Z 的牛顿迭代公式。Z 1 （ A1 A2 / TprA3 / Tpr 3） r +（ A4 +A 5 / Tpr ) r 2（2.22 ）( A5 A6 r5 ) / Tpr( A72 / Tpr3)(1 A8 r2 )exp( A8 r2 )r0 27ppr（2.23 ）rZTpr。4。frkr0.27 ppr / Tpr （ A1A2 / Tpr A3 / Tpr3） r2（+A 5 / Tpr )3( A5 A6 r6)Tpr+ A4r( A7r3/ Tpr3)(1A8r2 )exp(A8r2 ) 0f rk1 （A1A2 /

8、 Tpr A3 / Tpr 3）(2 r )（ 2.24 ）（A4+A5/ Tpr )(3r2 ) ( A5 A6 / Tpr )(6r5 )( A7 / Tpr 3 ) 3 r 2A8 (3 r 4 ) A822(2 r 6 ) e A8 rfk（ 2.25）k 1krrrfkrA10.31506237A21.0467099A30.5783229（ 2.26）A40.53530771A50.61232032A60.10488813A70.68157001A80.68446549（ 2.27 ）在已知 ppr 和 Tpr 的情况下，由 (2.2) 式求解 Z 时，采用迭代法。即首先给定的Z的一

9、个初值 Zo( 例如 Zo=1.0) ，由 (2.23) 式求出r，作为 (2.24)式迭代的初值。比较 r 与用 (2.26) 式计算所得的 rk+1 之值，如 r -rk+1<0.00001，则可将求得了 r 值代入 (2.2) 式求得 Z 值。否则，用最后求出的r 继续循环，直到r -rk+1<0.00001为止。2.2.4 DAK法该方法发表于 1975 年，方程如下：Z 1（A1A2 / TprA3 / Tpr3 +A4 / Tpr 4 +A5 / Tpr5） r +（A6 +A 7 / Tpr +A 8/ Tpr22（/ Tpr +A2)5（2.28 ）) rA9 A7

10、8 / TprrA10( r2 / Tpr3)(1A11r2 )exp(A11r2 )。5。A10.3265A21.0700A30.5339A40.01569A50.05165A60.5475（2.29 ）A70.7361A80.1844A90.1056A100.6134A110.7210解题方法和 DPR法步骤思路一样，但所用公式不同：F ( r )r0.27 ppr / Tpr （ A1A2 / Tpr A3 / Tpr 3 +A4 / Tpr 4 +A5 / Tpr 5） r 2 +（2) r3（A2)6A6 +A 7 / Tpr +A 8 / TprA97 / Tpr +A 8 / T

11、prrA10 ( r3 / Tpr3 )(1 A11 r2 )exp(A11r2 )（ 2.30 ）F ( r ) 1 （ A1A2 / TprA3 / Tpr 3 +A4 / Tpr 4 +A5 / Tpr 5） 2 r +（A6 +A 7 / Tpr +A 82) 32（2) 65（2.31 ）/ TprrA9 A 7 / Tpr +A 8/ Tprr( A10 / Tpr 3) 3 r 2A11(3 r 422 A11 r6 ) e A11 r此法适用于 1.0<Tpr<3.0 ，0.2<Ppr<3.02.2.5平均值法将以上计算方法结果累加除以计算方法的个数2.

12、3天然气压缩因子计算天然气的压缩系数就是指在恒温条件下，随压力变化的单位体积变化量，即Cg1V（2.32）Vp TCg气体压缩系数 ,1/MPa；V温度为 T 时气体体积随压力的变化率， m3/Mpa；P TV气体体积 ,m3；( 负号说明气体压缩系数与压力变化的方向相反。 )由真实气体的 PVT方程，得下式：。6。V nRTZ / p（2.33 ）经过一系列的推导及换算，得到天然气压缩系数表达式，如下所示：Cg10.27f1(Tpr ) 2 f 2(Tpr ) r（2.34）Ppr2Z 3Tpr Ppr1 ( r / Z) f1 (Tpr ) 2 f 2 (Tpr ) r式中：A2A3A5A6f1 (Tpr )A1 TprTpr3f 2(Tpr ) A4TprTpr3（ 2.35 ）2.4天然气体积系数计算天然气的体积系数就是指： 地层条件下某一摩尔气体占有的实际体积与地面标准条件下同样摩尔量气体占有的体积之比，由下式表示：VpscZTfBgR（2.36 ）VscpZscTsc式中VR 地层条件下气体的体积 ,m3；Vsc 地面标准状态下气体的体积 ,m3；Bg 天然气的体