天然气物性参数新.doc

2.1 天然气临界参数计算 天然气平均分子量天然气是混合气体，分子量不是一成不变的，其平均分子量按Key规则计算: 2.1式中 Mg天然气的平均分子量kg/mol；Mi、yi天然气中i组分的分子量和摩尔分数。 天然气的相对密度首先假定空气和天然气都取同一标准状态，天然气的相对密度可用下式表示： 2.2式中 rg天然气的相对密度；同一标准状态下，天然气、空气的密度kg/m3；天然气、空气的平均分子量kg/mol。 拟临界压力PPC和拟临界温度TPC 组分分析方法 2.3 式中 天然气组分i的临界压力(绝),MPa； 天然气组分i的临界温度,(273+t)°K。 相关经历公式方法在缺乏天然气组分分析数据的情况下，可引用Standing在1941年发表的相关经历公式对于干气 2.4对于湿气 2.5也可以用下面经历关系式进展计算对于干气 2.6对于湿气 2.7注意：上式是对于纯天然气适用，而对于含非烃CO2 、H2S 等可以用Wichert和Aziz修正。修正常数的计算公式为:2.1.4 拟比照压力PPr和拟比照温度TPr的计算比照参数就是指*一参数与其应对应的临界参数之比：即 2.92.2 天然气的偏差因子Z计算天然气偏差因子Z的计算是指在*一压力和温度条件下，同一质量气体的真实体积与理想体积之比值。 2.10计算天然气偏差因子方法较多，下面主要介绍几种常用的计算方法2.2.1 Pong.Robinson方程法 2.11式中 2.12 2.13 2.14 2.15 2.16 2.17式中 Kij天然气的交互作用参数；pcr组分i的气体临界压力；Tcr组分i的气体临界温度；Tr组分i的比照温度；i组分i的偏心因子。由方程可得到关于Z的方程 2.18 2.19 Cranmer方法 2.20 2.21式中 拟比照密度。P、T求Z，计算步骤如下：第一步 计算，；第二步 计算，；第三步 对Z赋初值，取Zo=1，利用式3.12计算第四步 将值代入式3.11，计算Z该方法适用于p<35MP的情况。 DPR法1974年,Dranchuk，Purvis和Robinson等人在拟合Standing.Katz图版的根底上,提出了计算偏差因子Z的牛顿迭代公式。 2.22 2.23 2.24 2.25 2.26 2.27在和的情况下，由(2.2)式求解Z时，采用迭代法。即首先给定的Z的一个初值Zo(例如Zo=1.0)，由(2.23)式求出，作为(2.24)式迭代的初值。比拟与用(2.26)式计算所得的k+1之值，如-k+1<0.00001，则可将求得了值代入 (2.2)式求得Z值。否则，用最后求出的继续循环，直到-k+1<0.00001为止。 DAK法该方法发表于1975年，方程如下： 2.28 2.29解题方法和DPR法步骤思路一样，但所用公式不同： 2.30 2.31此法适用于1.0<Tpr<3.0，0.2<Ppr<3.0 平均值法将以上计算方法结果累加除以计算方法的个数2.3 天然气压缩因子计算天然气的压缩系数就是指在恒温条件下，随压力变化的单位体积变化量，即 2.32Cg 气体压缩系数,1/MPa；温度为T时气体体积随压力的变化率，m3/Mpa；V 气体体积,m3；(负号说明气体压缩系数与压力变化的方向相反。)由真实气体的PVT方程，得下式： 2.33经过一系列的推导及换算，得到天然气压缩系数表达式，如下所示： 2.34式中： 2.352.4 天然气体积系数计算天然气的体积系数就是指：地层条件下*一摩尔气体占有的实际体积与地面标准条件下同样摩尔量气体占有的体积之比，由下式表示： 2.36式中 地层条件下气体的体积,m3； 地面标准状态下气体的体积,m3；天然气的体积系数,m3/m3(标)。在实际计算时，通常取=1.0，而当=0.101MPa，=293K时，由上式得： 2.372.5 天然气膨胀系数的计算 2.38式中 天然气膨胀系数