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 修正。
修正常数的计算公式为:()()()()()2222222pc pc 4.75460.21020.03 1.158310 3.06121084.9389188.49440.9333 1.4944g CO N H S gCO N p T γφφφγφφ--=-+-⨯+⨯=+-- (2.8)2.1.4 拟对比压力P Pr 和拟对比温度T Pr 的计算对比参数就是指某一参数与其应对应的临界参数之比:即pr pcp p p =Pr pcTT T =(2.9)2.2 天然气的偏差因子Z 计算天然气偏差因子Z 的计算是指在某一压力和温度条件下,同一质量气体的真实体积与理想体积之比值。
V Z V 实际理想=V nRTp实际=(2.10) 计算天然气偏差因子方法较多,下面主要介绍几种常用的计算方法2.2.1 Pong.Robinson 方程法()()()RT a T p V b V V b b V b =--++- (2.11) 式中0.511()()(1)nni j i j i j ij i j a T x x a a K αα===-∑∑ (2.12)1ni i j b x b ==∑ (2.13)220.45724cri crR T a p = (2.14) 0.0788cri crRT b p = (2.15) ()20.5ri 11i i m T α⎡⎤=+-⎣⎦ (2.16)20.37464 1.54220.26992i i i m ωω=+- (2.17)式中 K ij —天然气的交互作用参数; p cr —组分i 的气体临界压力; T cr —组分i 的气体临界温度; T r —组分i 的对比温度; ωi —组分i 的偏心因子。
由方程可得到关于Z 的方程()()023)1(32223=-----+--B B AB Z B B A Z B Z (2.18)22TR aPA =RT bPB =(2.19)2.2.2 Cranmer 方法3pr pr 223pr pr 1.04670.578310.315060.61230.53530.6815pr pr pr Z T T T T ρρρ⎛⎫=+-- ⎪ ⎪⎝⎭⎛⎫⎛⎫+-+⨯ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭ (2.20) ()()pr pr 0.27/pr p Z T ρ=⨯⨯ (2.21)式中 pr ρ——拟对比密度。
已知P 、T 求Z ,计算步骤如下: 第一步 计算pc p ,pc T ; 第二步 计算pr p ,pr T ;第三步 对Z 赋初值,取Z o=1,利用式(3.12)计算pr ρ 第四步 将pr ρ值代入式(3.11),计算Z 该方法适用于p<35MP 的情况。
2.2.3 DPR 法1974年,Dranchuk ,Purvis 和Robinson 等人在拟合Standing.Katz 图版的基础上,提出了计算偏差因子Z 的牛顿迭代公式。
3212pr 3pr r 45pr r 5232256r pr 7r pr 8r 8r 1//++A /)()/(/)(1)exp()Z T A T A T A A T A T A A ρρρρρρ=++++++-(A A )( (2.22)pr r pr027p ZT ρ⋅=(2.23)()k 32r r pr pr 12pr 3pr r 3645pr r 56r pr 33227r pr 8r 8r 0.27///++A /)()(/)(1)exp()0f p T T A T A T A A T A T A A ρρρρρρρρ=-++++++-=(A A )( (2.24)()28r k 3r 12pr 3pr r 2545pr r 56pr r 324267pr r 8r 8r 1//)+A /)(3)(/)(6)(/)3(3)(2)A f T A T A T A A T A T A A eρρρρρρρρ-'=+++++⎡⎤++-⎣⎦(A A )(2( (2.25)()()k r k 1krr k r f f ρρρρ+=-' (2.26)123456780.315062371.04670990.57832290.535307710.612320320.104888130.681570010.68446549A A A A A A A A ==-=-==-=-== (2.27)在已知pr p 和pr T 的情况下,由(2.2)式求解Z 时,采用迭代法。
即首先给定的Z 的一个初值Z o(例如Z o=1.0),由(2.23)式求出r ρ,作为(2.24)式迭代的初值。
比较r ρ与用(2.26)式计算所得的r ρk+1之值,如r ρ-r ρk+1<0.00001,则可将求得了r ρ值代入 (2.2)式求得Z 值。
否则,用最后求出的r ρ继续循环,直到r ρ-r ρk+1<0.00001为止。
2.2.4 DAK 法该方法发表于1975年,方程如下:34512pr 3pr 4pr 5pr r 222567pr 8pr r 97pr 8pr r 232210r pr 11r 11r 1//+/+/++A /+A /)A /+A /)(/)(1)exp()Z T A T A T A T A T T A T T A T A A ρρρρρρ=+++-++-(A A )(( (2.28)10.326521.07003 0.533940.0156950.0516560.54757 0.736180.184490.1056100.6134110.7210A A A A A A A A A A A ==-=-==-==-==== (2.29)解题方法和DPR 法步骤思路一样,但所用公式不同:3452r r pr pr 12pr 3pr 4pr 5pr r 232667pr 8pr r 97pr 8pr r 332210r pr 11r 11r ()0.27///+/+/++A /+A /)A /+A /)(/)(1)exp()F p T T A T A T A T A T T A T T A T A A ρρρρρρρρ=-+++-++-(A A )(( (2.30)211r 345r 12pr 3pr 4pr 5pr r 222567pr 8pr r 97pr 8pr r 324610pr r 11r 11r `()1//+/+/++A /+A /)3A /+A /)6(/)3(32)A F T A T A T A T A T T A T T A T A A e ρρρρρρρρ-=+++⨯⨯-⨯+⎡⎤+-⎣⎦(A A )2(( (2.31) 此法适用于1.0<Tpr<3.0,0.2<Ppr<3.02.2.5 平均值法将以上计算方法结果累加除以计算方法的个数2.3 天然气压缩因子计算天然气的压缩系数就是指在恒温条件下,随压力变化的单位体积变化量,即g 1TV C V p ⎛⎫∂=-⎪∂⎝⎭ (2.32)Cg —— 气体压缩系数,1/MPa ;TV P ∂⎛⎫ ⎪∂⎝⎭——温度为T 时气体体积随压力的变化率,m3/Mpa ; V —— 气体体积,m 3;(负号说明气体压缩系数与压力变化的方向相反。
) 由真实气体的PVT 方程,得下式:/V nRTZ p = (2.33)经过一系列的推导及换算,得到天然气压缩系数表达式,如下所示:12g 23pr 12()2()10.271(/)()2()pr pr r pr pr r pr pr r f T f T C P Z T P Z f T f T ρρρ⎧⎫+⎪⎪=-⎨⎬⎡⎤++⎪⎪⎣⎦⎩⎭ (2.34)式中:32113()pr pr pr A A f T A T T =-- 56243()pr pr prA A f T A T T =-- (2.35) 2.4 天然气体积系数计算天然气的体积系数就是指:地层条件下某一摩尔气体占有的实际体积与地面标准条件下同样摩尔量气体占有的体积之比,由下式表示:sc R g sc sc scfp ZT V B V pZ T ==(2.36) 式中 R V —— 地层条件下气体的体积,m3;sc V —— 地面标准状态下气体的体积,m3;g B ——天然气的体积系数,m3/m3(标)。