GB T 17747.2-1999 天然气压缩因子的计算 第2部分;用摩尔组成进行计算.pdf

《GB T 17747.2-1999 天然气压缩因子的计算 第2部分;用摩尔组成进行计算.pdf》由会员分享，可在线阅读，更多相关《GB T 17747.2-1999 天然气压缩因子的计算 第2部分;用摩尔组成进行计算.pdf（20页珍藏版）》请在麦多课文档分享上搜索。

1、GB/T 17747.2 1999 前兰主仨习本标准等效采用ISO12213-2，1997(天然气压缩因于的讨算用摩尔组成进行计算。本标准在技术内容上和编写格式与ISO12213-2，1997完全-致。本标准取消了国际标准中的附录F、附录G，本标准中高位发热量和相对密度采用的参比条件同我国石油气体所采用的标准参比条件有所差别，为方便使用，在4.4. 1和4.4.2增加了注2和注3。天然气压缩因于的计算标准包括以下3个部分t第1部分导论和指南h第2部分用摩尔组成进行计算);第3部分用物性值进行计算。本标准是第2部分(GB/T17747. 21999)，用摩尔组成进行汁算。本标准的附录A、附录B、

2、附录C、附录D是标准的附录。本标准的附录E是提示的附录。本标准由原中国石油天然气总公司提出。本标准由石油工业天然气专业标准化技术委员会归口并负责解释。本标准起草单位:中国石油天然气集团公司四川石油管理局天然气研究院本标准主要起草人:罗勤、陈腐良、曾文平、许文晓、富朝英、陈荣松。GB/T 17747.2一1999ISO前言ISO(国际标准化组织)是各国家标准化机构(lSO成员)组成的世界性的联合会。制定国际标准的工作通常由ISO技术委员会完成。对技术委员会提出的项目感兴趣的每个成员都有权参加。与ISO保持联系的各政府或非政府的国际性组织也可以参加此项工作。所有电工技术方面的标准化工作，ISO与I

3、EC(国际电工委员会)保持密切的合作。由技术委员会通过的国际标准草案交各成员进行表决投票，要求至少有75%的成员同意，才能作为国际标准正式发布。国际标准ISO12213-2是由天然气技术委员会ISO/TC193下的天然气分析分委员会制定的。ISO 12213天然气压缩因子的计算标准包括以下3个部分z第1部分g导论和指南;一第2部分:用摩尔组成进行计算;第3部分g用物性值进行计算。附录A、附录B、附录C、附录D是标准的附录。附录E、附录F、附录G是提示的附录。 353 1 范围中华人民共和国国家标准天然气压缩因子的计算第2部分:用摩尔组成进行计算Nalural gas-Calculalion o

4、f compression faclor Parl 2, Calculalion using molar-composition analysis G/T 17747.2一1999eqv ISO 12213-2 ,1997 本标准规定了天然气、含人工掺合物的天然气和其他类似混合物仅以气体状态存在时的压缩因子计算方法。该计算方法是用已知的气体的详细的摩尔分数组成和相关压力、温度计算气体压缩因子。该计算方法又称为AGA8-92DC计算方法，主要应用于在输气和配气正常进行的压力户和温度T范围内的管输气，计算不确定度约为土0.1%。也可在更宽的压力和温度范围内，用于更宽组成范围的气体，但汁算结果的不确

5、定度会增加(见附录E)。有关该计算方法应用范围和应用领域更详细的说明见GB/T17747.1。2 号|用标准下列标准所包含的条文，通过在本标准中引用而构成为本标准的条文。本标准出版时，所示版本均为有效。所有标准都会被修订，使用本标准的各方应探讨使用下列标准最新版本的可能性。GB/T 3102.3-1993 力学的量和单位GB/T 3102.4一1993热学的量和单位GB/T 11062-1998 天然气发热量、密度、相对密度和沃泊指数的计算方法(neqISO 6976 ,1 995) GB/T 17747.1-1999 天然气压缩因于的计算第1部分:导论和指南(eqvISO 12213-1 ,

6、1 997) 3 定义相关于本标准的所有定义见GB/T17747. L文中出现的符号所代表的含义及单位见附录A，4 计算方法4. 1 原理AGA8.92DC计算方法所使用的方程是基于这样的概念.管输天然气的容量性质可由组成来表征和计算。组成、压力和温度用作计算方法的输入数据。该计算方法需要对气体进行详细的摩尔组成分析。分析包括摩尔分数超过O.000 05的所有组分。对典型的管输气，分析组分包括碳数最高到C，或C，的所有短类，以及NCO，和He。对其他气体，分析需要考虑如H，O蒸气、H，S和C，H等组分。对人造气体.H，和CO也可能是重要的分析组分。4.2 AGA8.92DC方程AGA8-92D

7、C计算方法使用AGA8详细特征方程(下面表示为AGA8-92DC方程，见GB/T 17747.1);该方程是扩展的维利方程，可写作国l!t质量技术监督局1999-05-17批准1999-12-01实施，) I GB/T 17747.2-1999 Z二1+ Bpm -p，2:c:十三C:(b, 呐卢卡，)!1.exp( - c ,1;.) .=13 ,=13 式中，Z压缩因于pB 第二维利系数;Pm摩尔密度(单位体积的摩尔数hp, 对比密度sb.、k，一一常数，见表B1; (n* 温度和组成的函数的系数。对比密度p，同摩尔密度Pm相关，两者的关系由下式给出.p, = K3Pm 、2) 式中zk

8、混合物体积参数。摩尔密度表示为:Pm = p/(ZRT) . ( 3 ) 式中:户绝对压力sR 摩尔气体常数;T 热力学温度。压缩因子Z的计算方法如下:首先利用附录B给出的相关式计算出B和C:(n=13-58)。然后通过适当的数值ii-:算方法，求解联立方程(1)和(3)得到pm和Z.计算程序流程见图B1.4. 3 输入变量AGA8-92DC计算方法要求输入的变量包括绝对压力、热力学温度和摩尔组成。摩尔组成是以摩尔分数表示下列组分，NCOAr，CH.、C2Hofl, C3Ha n-C4H10-. i-C4HlO, n-CsH12, i CSH12、C6H14、C7H16、CSH18、C9H20

9、、C1oH22、H，、CO、H，S、He、0，和H，O。注1，如果C7H16、C8H18C9H20、ClOHU摩尔分数未知，允许用ct表示且的摩尔分数目应进行敏感度分析，以检验此近似法是否会使计算结果变差。摩尔分数大于o.000 05的所有组分都必须在计算中考虑。痕量组分(如C，H，等)应按表l中指定的赋值组分处理。所有组分的摩尔分数之和为1士O.000 1。如果已知体积分数组成，则应将其换算成摩尔分数组成，具体换算方法见GB/T11062。表1微量和痕量组分一览表微量和痕量组分0 , A, H,S CZH4 , CzH2 丙烯、丙二烯丁烯、丁二烯新戊饶、戊烯、苯、环戊烧C同分异构体、环己饶、

10、乙苯、二甲苯c，同分异构体、环庚烧、甲苯c.同分异构体C同分异构体c，.同分异构体和更高碳数娃类4. 4 应用范围4.4.1 管输气AGA8-92DC计算方法对管输气的应用范围如下:指定赋值组分0 , A, H,S co, C3Ha n-C4H!o n-C6H12 n-CsHa n-C1H16 n-CSH18 n-C9Hzo n-C1oH22 355 GB/T 17747.2-1999 绝对压力z热力学温度:高位发热量相对密度:o MPal 2 士0.001.TC02 土0.001XCH4 :tO.001 XC2H6 士0.001XC:3He 土。.0005 I巳4Hl0土0.0003 zcf

11、 :t0. 000 1 IH，和Xco:t o. 001 4.5.4 结果的表述压缩因子和摩尔密度计算结果应保留至小数点后四位或五位，同时给出压力和温度以及所使用的计算方法(GB/T17747. 2.AGA8-92DC计算方法。验证计算机程序时，压缩因子计算结果应给出更多的位数。:);)8 符号a., B B;i) b. C, c; Ei E EIJ E:; F Fi F, f. G G, G, Gi) G; g. H K K, K, K, k. M GB/T 17747.2-1999 附录A(标准的附录)符号和单位含义悔数(表B1)第二维利系数混合物交互作用系数C方程(B1)和方程(B2)J

12、常数(表B1) 常数(表B1) 与温度和组成相关的系数组分z的特征能量参数(表B2)组分J的特征能量参数(表B2)第二维利系数的二元能量参数第二维利系数的二元能量交互作用参数(表B3)混合物高温参数组分z的高温参数(表B2)组分J的高温参数(表B2)常数(表B1)混合物定位参数组分z的定位参数(表B2)组分J的定位参数(表B2)二元寇位参数二元定位交互作用参数(表B2)常数(表B1)高位发热量体积参数组分z的体积参数(表B2)与雹组分j的体积参数(表B2)二元体积交互作用参数(表B3)常数(表B1)摩尔质量单位m3 kmol-1 K K K MJ m-3 (m3 kmol-1)1/3 (m k

13、mol-1 )1/3 (m3 kmol-1)1/3 kg kmol-1 359 GB/T 17747.2-1999 M, 组分z的摩尔质量kg kmol-1 N 气体混合物的组分数n 整数(158) 户绝对压力MPa Q 四极参数Q, 组分z的四极参数Q, 组分J的四极参数q 常数(表Bl)R 气体常数(=0.008314510)MJ kmol-1 K 1 5, 组分2的偶极参数(表B2)5, 组分1的偶极参数(表B2)5 常数(表B1)T 热力学温度K U 混合物能量参数K UJ 混合物二兀能量交互作用参数(表B3)U 常数(表Bl)W, 组分z的组合参数(表B2)Wj 组分1的组合参数表B

14、2)W 常数(表B1)x , 组分z在气体混合物中的摩尔分数X , 组分J在气体混合物中的摩尔分数Z 压缩因子p 质量密度kg m-3 p, 气体的对比密度p. 摩尔密度kmol m 附亵B(标准的附录)AGA8-92DC计算方法描述B1 概述计算气体混合物压缩因子Z使用4.2给出的AGA8-92DC方程。本附录将详细描述用AGA892DC方程进行压缩因子计算的有关方法和计算机执行程序，并给出必要的常数值。验证计算机程序用;.HO GB/T 17747.2-1999 的压缩因子数据见附录C。如果计算机程序能够得到与附录C中数据相等的计算结果，则可使用。B2 AGA8-92DC计算方法的计算机执

15、行程序B2.1 计算概述B2. 1. 1 输人热力学温度T.绝对压力户和混合物中各组分的摩尔分数x;注B，当压力和温度用非MPa和K的其他单位表示时，必须分别换算成以MPa和K表示的值，有关换算因于见GB/T 3102. 31993和GB/T3102. 41993以及附录D，B2. 1. 2 计算状态方程系数B和C:(n二1358).两者均取决于T和Xi;B2. 1. 3 利用改写的状态方程，迭代求解摩尔密度队，以得到压力p;B2. 1.4 当由B2.1. 3计算出的压力与B2.1. 1的输人压力，在规定的收敛范围内(如1X10-)相一致时，即得到压缩因子计算值。计算程序流程见图B10输入T,

16、 , Xi 计算状事方程罩黠B相C:(.-13-58) 草帽曹E不成立输出Z=p/(现RT)图B1AGA8-92DC计算方法的程序流程图B2.2 计算步骤B2. 2. 1 输入绝对压力户、热力学温度T和天然气混合物中各组分的摩尔分数X;B2. 2. 2 根据B2.2. 1中输人的热力学温度T和天然气组分摩尔分数品，计算与组成和温度有关的状态方程系数E和c:(n=1358); 第二维利系数B按方程(81)计算: 18 N N B = L;a.T-. L; L;XiXjB;jEi;(K凡)1i) l B ( -l ;=1 j=l B:i, = (Gi，十1 g.),.(Q,Q, + 1 q.队(F

17、!Fj+ 1 /.)f,(SiSj + 1 s.),(WiWj + 1四.)飞( B2 ) 二元参数Ei，和Gij由方程(B3)和(84)计算gEij = Eij (EiEj)比Gij = Gij (Gi + G) /2 . ( B3 ) . ( B4 ) 注日2，除了表B3中给出的E;和G;外，所有其他二元交互作用参数E;和G，j的值都是1.0, 系数c:(n=1358)由方程(B5)计算c; = a.(G + 1 g.).(Q + 1 q.).(F + 1 /.)U . T-(B5 ) 用以下共形求解混合方程计算混合物参数U、G和Q;二重加和时.1从lNl，变化，而相对每一个t值.从z十l

18、N变化。-lfl GB/T 17747.2-1999 N N-1 N U .L;x，E盯2十2.L;.L;x;x,(uij -l)(E盯%. . . .( B6 ) ,=1 ;=1 j=i+l N N-l N G .L;x十2.L;.L; X;Xj(G;j一l)(G.+Gj). ( B7 ) Q .L;x;Q; ( B8 ) i= 1 F三点，( B9 ) ;-1 注B3，除了表B3给出的K川，E，;.G，和U中H矶2的F(H比2)=1.O.而其他组分的F，值都为0，革的W(H)=1.0，而其他组分的W，值都为O.B2. 2. 3 在压缩因子Z的计算中，气体的组成岛、热力学温度T和绝对压力户都

19、是已知的，问题在于要用表示压力户的状态方程来计算摩尔密度Pm.将定义压缩因子Z的方程(1)代人方程(3)(见4.2)，获得压力的状态方程(810)，p PmRTl + Bpm - p, .L;C: +三C:(b.仇fI，.)p;exp(-c.仲)J .( B10 ) 方程(B10)用标准状态方程密度检索法求解。由于已获得压力户的表达式方程(B10汀，则求解摩尔密度pm使计算出的压力与输入的压力两者的差值在预先设定的范围(如1X10-)以内。对比密度冉通过混合物体积参数K与摩尔密度pm相关联C见2中方程(2)J。混合物体积参数K由方程(811)计算N N-l N K5 .L;x;K严J+ 2.L

20、; .L; x;xj(K!j一l)(K几;-1 j-i+l . ( B11 ) 注B4，求和时，下标z指的是气体混合物中第z个组分g下标J指的是气体混合物中第J个组分，N指的是混合物中的组分数。单重求和中，I是1-N阔的整数值.例如，对吉12个组分的混合物.N需=-12，单重求和中将有12项，二重求和中z从1-N-1变化，而相对每一个a值.J从z十1-N变化.例如，对吉12个组分的混合物，如果托，的值都平为1.0，则二重求和将有66项.由于许多瓦，的值都为1.0，因此对许多天然气混合物，二重求和中非零项的数目很少.除了表B3给出的Koj外，所有其他的瓦，的值都为1.O. B2. 2. 4 求出

21、摩尔密度pm后，利用压力、温度、摩尔密度和摩尔气体常数计算压缩因子。Z = p/(PmRT) ( B12 ) 密度p可由方程(B13)计算zp=M卢刚. ( B13 ) 式中，M根据方程(B14)计算zM = .L;x;M; ( B14 ) 密度值应保留至小数点后第三位。表B1状态方程参数n a. b. C. 是.U. g瞻q. f. S. 四1 O. 153 832 600 l 。0.0 。2 1. 341 953 000 1 。0.5 。3 2. 998 583 000 l 。1. 0 。4 一0.018312 280 1 。3. 5 。5 0.375 796 500 1 。一0.51 。

22、6 - 1. 589575000 1 。4. 5 1 。7 -0.053588470 l 。O. 5 。I 。丛(2nud nd n叫Jw-1-一一;一续7-( an-I T-B 1一表T-,ia nu G-n a. b. C. 是.u饨g. q. f . 飞足.8 0.886 594 630 1 。7.5 。1 。9 0.710237 040 l 。9. 5 。1 。1. 471 722000 I 。6.0 。1 10 11 1. 321 850 350 1 。12.0 。l 12 0.786 659 250 l 。12. 5 。1 13 2.291 290 X 10 1 l 3 一6.0。

23、1 。14 0.157672 400 1 1 2 2.0 。l 2 3.0 。1 5 0.436 386 400 1 16 一0.044081 590 1 1 2 2.0 。1 。17 0.003 433 888 1 1 4 2.0 。18 0.032 059 050 1 1 4 11. 0 。19 0.024 873 550 2 。o. 5 。20 。.073322 790 2 。0.5 。21 一0.001600 573 2 1 2 0.0 。22 O. 642 470 600 2 1 2 4.0 。23 -0.416260100 2 1 2 6.0 。24 0.066 899 570 2

24、 1 4 21. 0 。25 o. 279 179 500 2 1 4 23.0 1 。26 0.696 605 100 2 1 4 22.0 。l 。27 O. 002 860 589 2 l 4 一1.0。l 。28 O. 008 098 836 3 。一0.5。1 。29 3.150547000 3 1 1 7.0 1 。30 O. 007 224 479 3 1 l 一1.0。l 。31 0.705 752 900 3 l 2 6.0 。32 O. 534 979 200 3 1 2 4.0 l 。33 0.079 314 910 3 1 3 1. 0 1 。34 一1.4184650

25、003 1 3 9.0 1 。35 5.999 05X 10- 3 1 4 一13.0 。1 。36 0.105 840 200 3 1 4 21. 0 。37 0.034 317 290 3 1 4 8.0 。1 。38 一O.007 022 847 4 。一0.5。39 0.024 955 870 4 。0.0 。40 0.042 968 180 4 1 2 2.0 。41 O. 746 545 300 4 1 2 7. 0 。42 0.291 961 300 4 1 2 9.0 。1 。36:-: GB/T 17747.2-1999 表Bl(完) a b 是U g q匍f 四13 7.2

26、94616000 4 1 4 22.0 。44 一一9.936757 000 4 1 4 23.0 。45 0.005 399 808 5 。1. 0 。46 0.243 256 700 3 1 2 9.0 。17 0.049 870 160 5 l 2 3. 0 。1 。48 0.003 733 797 5 1 4 8.0 。49 1. 874 951 000 5 I 4 23.0 。1 。50 0.002168144 6 。1. 5 。51 u. 658 716 400 6 2 5.0 1 。52 0.000205 518 7 。一O.5 。1 。53 0.009 776 195 7 1

27、z 4.0 。54 。020487 080 8 1 1 7. 0 1 。55 0.015 573 220 8 2 3. 0 。一56 0.000 862 415 8 1 2 。1 。57 0.001 226 752 9 1 2 1. 0 。8 O. 002 850 908 9 2 。1 。表B2特征参数识摩尔质量能量参数别化告物M, E , 体积垂数K, 定位参数四极参数高温参数偶极参数组合垂数号kg kmol-I K (m3 kmol-1)1/ G, Q, F S, W, 1 CH, 16.0430 151. 318300 0.461 925 5 。0, 0 0.0 。1- -2 I N,

28、28.013 5 99.737 780 0.4479153 0.027 815 0.0 0.0 0.0 0, 0 3 CO, 44.0100 24 1. 960 600 0.455748 9 0.189 065 0.690 000 0.0 。0, 0 4 CZH6 30. 070 0 244.166700 0.527 920 9 0, 079 300 0.0 0.0 0.0 。5 C3H8 44.097 0 298.118 300 0.5837490 0.141 239 。0.0 。0.0 。H,O 18.0153 514.015 600 o. 382 586 8 0, 332 500 1.

29、067 750 0.0 1.582 200 1. 0 7 HS 34, 082 0 296.355 000 0.461 826 3 O. 088 500 0, 633 276 0.0 0, 390 000 O. 0 8 H, 2.0159 26.957 940 0.351 491 6 0, 034 369 。1. 0 0.0 O. 0 9 CO 28.0100 105. 534 800 0.453 389 4 0, 038 953 。0.0 O. 0 O. 0 10 0 , 31. 998 8 122.766700 0.418 695 4 O. 021 000 。0.0 0.0 。11 i-C

30、HIU 58. 123 0 324.068 900 O. 640 693 7 O. 256 692 0.0 。0.0 12 II-C .j H10 58.123 0 337.638900 。.6341423O. 281 835 。.0 。0, 0 L一13 t于C5H1272.150 0 365.599 900 0.673 857 7 O. 332 267 0.0 0, 0 。0.0 L一L一14 lI-C,H12 72.1500 370, 682 300 0.679 830 7 0.366 911 。0.0 。0.0 15 lI-C而H86.1770 402, 636 293 0.71751

31、18 0.289 731 0.0 0.0 O. 0 O. 0 GB/T 17747.2-1999 表B2(完)识摩尔质量能量参数体积参数定位参数四极参数高温参数别化合物M. E K , G, Q, F , s W , l芝kg. kmol-1 ! K (m3 kmor-1)11 100.2川1427.7川O16 n宁C;H160.7525189 o. 337 542 0.0 。0.0 。17 n-CHI8 114.231 0 450.325 022 0.784 955 0 0.383 381 0.0 。0.0 18 n-C9H211 128.2580 470.840891 0.815 273

32、1 0.427 354 0.0 。0.0 。19 n.CJuHn 142.285 0 489.558 373 0.843 782 6 0.469 659 0.0 。0.0 20 Hc 4.002 6 2.610111 0.358 988 8 。0.0 0.0 0.0 0.0 A, 39.948 0 119.629 900 O. 421 655 1 0.0 0.0 0.0 。O. 0 表B3二元交互作用参数值识别号化合物对E,; Uj K CJ CH,+ N, 0.971 640 0.886 106 1. 003 630 CO, O. 960 644 0.963 827 0.995 933 O.

33、 807 653 唱.4 CZH6 U C3Hs 0.994 635 0.990 877 1. 007619 6 H,O 0.708 218 7 H,S 0.931 484 0.736 833 1. 000 080 卜8 H, 1. 170520 1. 156 390 1. 023 260 1.957 310 日CO 0.990 126 卜10 0 , 11 i-C.H 10 1. 019 530 12 n-C4H10 O. 989 844 0.992 291 0.997 596 13 i-CsH12 1. 002 350 14 n-CSHI2 0.999.268 1. 003 670 1.

34、002 529户15 n-C 6HH 1. 107274 1. 302 576 0.982 962 16 n骨C1H160.880 880 1. 191 904 0, 983 565 17 n-CSHI8 0.880 973 1. 205 769 0.982 707 18 n-C9H20 0.881 067 1. 219 634 O. 981 849 19 n-C1oHn 0.881 161 1. 233 498 0.980 991 2 3 N,+ CO, 1. 022740 O. 835 058 O. 982 361 0.982746 -一4 C2H5 0.970 120 0.816 431

35、 1. 007 960 5 C3H8 O. 945 939 0.915 502 6 H,O 0.746954 7 H,S 0.902 271 0.993476 0.942 596 8 H, 1. 086 320 0.408 838 1. 032 270 l(; 识别号化合物对, J 9 CO 10 0 , 11 i-CH10 12 n-CHIO 13 i-CsH12 14 n-CsH12 3 4 CO,+ C2He 5 C3H8 6 H,O 7 H,S 8 H, 9 CO 10 0 , 11 i-C4H10 12 n-C4H10 13 i-CsH12 14 n-CSH12 15 n-CsHI4

36、 16 n-C1H16 17 n-CSHI8 18 n-C9Hzo 19 n-C1oHn 4 5 CzHs+ CsH 6 H,O 7 H,S 8 H, 9 CO 10 0 , 11 i-CtH10 12 n-CtH10 13 i-C5Ha 14 n-CsH12 5 自C3Ha十H, 12 n-CH10 :3G GB/T 17747.2-1999 表B3(续E,j Uj 1. 005710 1. 021 000 0.946 914 0.973 384 O. 993 556 O. 959 340 O. 945 520 O. 925 053 0.969 870 0.960 237 0.849 408

37、 O. 955 052 1. 045 290 1. 281 790 1. 500 000 0.900 000 0.906 849 O. 897 362 0.726 255 0.859 764 0.855 134 1. 066 638 O. 831 229 1. 077 634 。.808310 1. 088 178 0.786 323 1. 098 291 0.765171 1. 108 021 1. 022 560 1. 065173 0.693 168 0.946 871 0.971 926 1. 164460 1. 616 660 1. 250 000 1. 013 060 1. 250

38、 000 1. 250 000 1. 005 320 1. 250 000 1. 034 787 1. 004 900 K ij Gii 1. 008 510 0.370 296 1. 673 090 1. 007 790 0.910 183 O. 895 362 0.881 152 0.867 520 0.854 406 O. 986 893 0.999 969 1. 020 340 识别号化告物对g J 7 15 H ,S+ n-CsH 14 16 n-C?HI6 17 n-CSH 18 18 n-C9H20 19 n-C1oHU 8 9 H ,+ co 10 0 , 11 i-C4H10

39、 12 n-CtH10 GB!T 17747.2-1999 表B3(完)E,; U ij 1. 008 692 1. 028 973 1. 010 126 1. 033 754 1. 011 501 1. 038 338 1. 012 821 1. 042 735 1. 014 089 1. 046 966 1. 100 000 1. 300 000 1. 300 000 附录C(标准的附录计算示例K;j G;j 0.968 130 0.962 870 0.957 828 。.952441 o. 948 338 用经过验证的计算机程序，对下述实例(见表Cl)进行了压缩因子计算，所得出的压缩因子

40、计算结果(见表C2)可作为编制的计算机程序验证用数据。如果编制的计算机程序能够得到与之相等的计算结果，则可使用。该经验证的计算机程序包含附录B中所描述的计算机执行程序。表Cl以摩尔分数表示的气体组成分析数据气体组成l气样2气样3气样俨气样5气样6气样1COE 0.006 0.005 0.015 0.016 0.076 0.011 XNZ 0.003 0.031 0.010 0.100 0.057 0.117 .rH, 0.00 0.00 0.00 0.095 0.00 0.00 .rco 0.00 0.00 0.00 0.010 0.00 0.00 .rCH4 0.965 0.907 0.85

41、9 0.735 0.812 0.826 XCZH6 0.018 0.045 0.085 0.033 0.043 0.035 .xC3HS 0.004 5 0.008 4 0.023 0 0.007 4 0.009 0 O. 007 5 .r, -C4H10 0.0010 0.001 0 0.003 5 O. 001 2 O. 001 5 0.001 2 Xn-C,H10 O. 001 0 0.001 5 0.003 5 0.001 2 O. 001 5 O. 001 2 Xi-C5H12 0.000 5 0.000 3 0.000 5 0.000 4 0.00 0.000 4 .:r-，SHI

42、2 0.000 3 0.000 4 O. 000 5 0.000 4 0.00 O. 000 4 XCSIl14 0.000 7 0.000 4 0.00 0.000 2 0.00 0.000 2 -XC,H16 0.00 0.00 0.00 0.000 1 0.00 O. 000 1 rCSH18 0.00 0.00 0.00 0.000 1 0.00 0.00 :l ( 7 GB/T 17747. 21999 表C2压缩因子Z计算结果输配条件户7 1气样f气样俨气样俨气样5气样俨气样肌1PaK 6 270 0.840 53 0.833 48 O. 793 80 o. 885 50 O. 8

43、26 09 0.853 80 6 280 0.861 99 O. 855 96 O. 822 06 O. 901 44 0.849 69 。.873 70 6 290 O. 880 06 0.87484 O. 845 44 0.915 01 0.869 44 O. 890 52 6 310 0.908 67 O. 904 46 0.881 83 0.936 74 O. 900 52 0.91723 6 330 0.930 11 O. 926 96 O. 908 68 0.953 18 0.923 68 0.937 30 12 270 0.721 33 0.71044 O. 641 45 0.8

44、10 24 0.695 40 0.750 74 12 280 O. 760 25 O. 750 66 0.689 71 0.837 82 0.737 80 O. 785 86 12 290 0.79317 。8475 0.731 23 0.86137 。773690.81569 12 310 0.845 15 0.838 63 0.796 98 0.899 13 O. 830 22 0.863 11 12 330 O. 883 83 0.878 70 0.845 53 0.927 66 。.87211 O. 898 62 附录D(标准的附录)压力和温度的换算因子如果输入压力和温度的单位不是MP

45、a和K，为了使用附录B描述的计算机执行程序进行计算，应作单位换算。部分换算因子如下:压力:户(MPa)= 户(bar)J X 10-1 p(MPa)=户(alm)J X 0.101 325 户(MPa)=p(psia)J/145.038p(MPa)=户(psig)+14.695 9J/145. 038 温度T(K)=t(C)+273.15 T(K)=t(于)-32J/1.8十273.15 T(K)=t(OR)J/l. 8 附亵E(提示的附录)更宽范固的应用效果在温度263-338K、压力最高至30MPa的范围内，利用实验测定值数据库，对AGA8-92DC计算方法进行了全面检验，这些实验数据取自

46、给定组成范围内的管输气(见4.4.1)，在此组成范围内压缩因于计算的不确定度已在4.5中给出。对更宽范围(对组成而言)天然气i见4.4. 2)压缩因子计算的不确定度的粗略估计，见图E1-图E4。图中以压力为纵坐标，横坐标分别为N2.C02.C2H，和C，H，含量。.; (; s GB/T 17747.2-1999 30 ZE. AGA8.92DC方程T263-338K A飞nR国B 20 A 10 A 6.z豆士0.1%B土o.1问士0.2%。o. 1 0.2 o. 3 O. 4 O. 5 N，摩尔分量XNZ)图El计算高含N，天然气压缩因子时不确定度的估计范围图El图E4显示了在压力最高为3

47、0MPa的范围内.AGA8-92DC计算方法的应用效果。不确定度范围取决于压力、温度和组成，也强烈地受到相界位置的影响。图El图E4中给出的不确定度极限是在最坏情况下获得的结果，也就是说它们不是最佳选择。实验数据不足以决定不确定度范围的边界位置时，用虚线将所估计的不确定度区划为两个区域。气体全组成对相界位置会有强烈的影响，使用者应当进行相界计算。30 L E AGA8-92DC方程T=263-338K e、4国哥B A 6.z士0.1%B土Q.l%-士0.2%C土O.2Yo-土0.5%D土O.SYo-士3.0%A 10 0 0 O. 1 O. 2 O. 3 CO，_尔分量(XC02) 图E2计算高含CO，天然气压缩因子时不确定度的估计范围:.; fj 9 GB!T 17747.2-1999 0 3 占24刊出AGA8-92DC方程T263-338K B 20 ，卢C 10 A B A .Z$.三土0.1%B土0.1%-土Q.2%C士0.2%-士0.5%。O. 05 o. 1 0.15 0, 2 CzHs摩如分Ik(xC2H) 图E3计算高含C2H，天然气压缩因子时不确定度的估计范围0 3 tE4刊出AGA8.92DC方程T263-338K 20 B 10 飞飞飞飞A A 6z运士0.1%B士0.1%-士0.2%。O. 05 o