GBT17747.2-2011天然气压缩因子的计算用摩尔组成进行计算.pdf 24页

'lcs75．060E24囝园中华人民共和国国家标准GB／T17747．2—2011代替GB／T17747．2—1999天然气压缩因子的计算第2部分：用摩尔组成进行计算Naturalgas--Calculationofcompressionfhcto卜Part2：Calculationusingmolar-compositionanalysis2011-12-30发布(ISO12213-2：2006，MOD)2012．06一01实施宰瞀鳃紫瓣訾襻赞星发布中国国家标准化管理委员会况19 前言⋯⋯⋯⋯⋯⋯··1范围··⋯⋯⋯⋯·-2规范性引用文件·3术语和定义⋯⋯·4计算方法⋯⋯⋯-附录A(规范性附录)附录B(规范性附录)附录C(规范性附录)附录D(规范性附录)附录E(资料性附录)目次GB／T17747．2—2011⋯·⋯⋯··-⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯·--⋯⋯⋯··-⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯Ⅲ-⋯··⋯⋯--⋯·--·⋯⋯··⋯···⋯··-·⋯···⋯⋯⋯⋯--·⋯-···⋯⋯⋯··⋯··--⋯·-··⋯·1⋯⋯⋯··-·⋯··⋯⋯⋯⋯···⋯⋯·⋯一·⋯⋯·-·-⋯·--⋯····⋯··-⋯·⋯⋯··⋯·-·⋯···-1⋯··-·-⋯-⋯⋯⋯⋯···-⋯⋯⋯··⋯⋯⋯⋯··-··⋯···⋯·⋯⋯--⋯·⋯⋯⋯·⋯··⋯···1⋯-·⋯···⋯·⋯·⋯···-·⋯⋯⋯····⋯·--⋯·····⋯···⋯⋯·⋯--··⋯--⋯-·⋯···⋯⋯-··1符号和单位⋯⋯⋯⋯⋯·⋯⋯·⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯·7AGA8—92DC计算方法描述⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯·⋯⋯⋯⋯一9计算示例⋯⋯⋯⋯⋯⋯···⋯⋯⋯⋯⋯⋯⋯⋯·⋯⋯···⋯⋯⋯··⋯⋯⋯⋯⋯⋯16压力和温度的换算因子·⋯⋯⋯⋯·⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯·⋯⋯⋯⋯17更宽范围的应用效果⋯⋯·⋯⋯⋯⋯·⋯⋯··⋯⋯⋯⋯⋯⋯·⋯⋯⋯⋯·⋯⋯⋯18 刖菁GB／T17747．2—2011GB／T17747《天然气压缩因子的计算》包括以下3个部分：——第1部分：导论和指南；——第2部分：用摩尔组成进行计算；——第3部分：用物性值进行计算。本部分是第2部分。本部分按照GB／T1．1—2009给出的规则起草。本部分代替GB／T17747．2--1999(天然气压缩因子的计算第2部分：用摩尔组成进行计算》。本部分与GB／T17747．2—1999相比主要变化如下：一一按ISO12213．2：2006修改了表1“微量和痕量组分一览表”中的内容；一一修改了图中的符号和图注；——删除正文中不确定度数值前的“士”号。本部分使用重新起草法修改采用ISO12213—2：2006((天然气压缩因子的计算第2部分：用摩尔组成进行计算》。本部分与IsO12213-2：2006的主要差异是：——第2章规范性引用文件中，将一些适用于国际标准的表述修改为适用于我国标准的表述，ISO标准替换为我国对应内容的国家标准，其余章节对应内容也作相应修改；——在4．4．1和4．4．2增加了将高位发热量和相对密度换算为我国天然气标准参比条件下相应值的注；——将B．2．1．1中注的内容移至4．3最后一段；——删除正文中不确定度数值前的“士”号；一一删除ISO标准前言，重新起草本部分前言；——删除第5章的内容；一一删除附录F和参考文献；本部分由全国天然气标准化技术委员会(sAc／Tc244)归口。本部分起草单位：中国石油西南油气田分公司天然气研究院、中国石油西南油气田分公司安全环保与技术监督研究院。本部分主要起草人：罗勤、许文晓、周方勤、黄黎明、常宏岗、陈赓良、李万俊、曾文平、富朝英、陈荣松、丘逢春。Ⅲ 1范围天然气压缩因子的计算第2部分：用摩尔组成进行计算GB／T17747．2—2011GB／T17747的本部分规定了天然气、含人工掺合物的天然气和其他类似混合物仅以气体状态存在时的压缩因子计算方法。该计算方法是用已知气体的详细的摩尔分数组成和相关压力、温度计算气体的压缩因子。该计算方法又称为AGA8—92DC计算方法，主要应用于在输气和配气正常进行的压力户和温度T范围内的管输气，计算不确定度约为0．1％。也可在更宽的压力和温度范围内，用于更宽组成范围的气体，但计算结果的不确定度会增加(见附录E)。有关该计算方法应用范围和应用领域更详细的说明见GB／T17747．1。2规范性引用文件下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本(包括所有的修改单)适用于本文件。GB3102．31993力学的量和单位GB3102．4—1993热学的量和单位GB／T11062--1998天然气发热量、密度、相对密度和沃泊指数的计算方法(ISO6976：1995，NEQ)GB／T17747．1—2011天然气压缩因子的计算第1部分：导论和指南(IsO12213-1：2006，MOD)3术语和定义GB／T17747．1给出的术语和定义适用于本文件。文中出现的符号所代表的含义及单位见附录A。4计算方法4．1原理AGA8—92DC计算方法所使用的方程是基于这样的概念：管输天然气的容量性质可由组成来表征和计算。组成、压力和温度用作计算方法的输人数据。该计算方法需要对气体进行详细的摩尔组成分析。分析包括摩尔分数超过0．00005的所有组分。对典型的管输气，分析组分包括碳数最高到ct或cs的所有烃类，以及N。、CO：和He。对其他气体，分析需要考虑如HzO蒸气、HzS和CzH。等组分。对人造气体，Hz和CO也可能是重要的分析组分。4．2AGA8-92DC方程AGA8-92DC计算方法使用AGA8详细特征方程(下面表示为AGA8—92DC方程，见GB／T17747．1)；1 GB／T17747．2—2011该方程是扩展的维利方程，可写作方程(1)：1858Z一1+BP。一Pr∑c：+∑c：(6。一c。6妒：·)硭·exp(一‘中：一)⋯⋯⋯⋯(1)_一15^一13式中：Z——压缩因子；B——第二维利系数；P。——摩尔密度(单位体积的摩尔数)；pr——对比密度；b。，c。，k。——常数(见表B．1)；c：——温度和组成函数的系数。对比密度m同摩尔密度m相关，两者的关系由方程(2)给出p。=K3P。式中：K——混合物体积参数。摩尔密度表示为方程(3)：口。一p／(ZRT)⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯(3)式中：p——绝对压力；R——摩尔气体常数；T——热力学温度。压缩因子z的计算方法如下：首先利用附录B给出的相关式计算出B和c：(n一13～58)。然后通过适当的数值计算方法，求解联立方程(1)和(3)得到舳和z。计算程序流程见图B．1。4．3输入变量AGA8—92DC计算方法要求输入的变量包括绝对压力、热力学温度和摩尔组成。摩尔组成包括下列组分的摩尔分数：N：、C02、Ar、CH。、C2H。、C。H。、n—c。H。。、i—c。H。。、n-C；H。：、i—C5H12、C6HH、C7H16、C。H18、C9Hzo、C10H22、Hz、CO、H2S、He、Oz和H20。注：如果GH。叫CH·㈨CHz州C。Hn摩尔分数未知，允许用C。+表示总的摩尔分数．应进行灵敏度分析，以检验此近似法是否会使计算结果变差．摩尔分数大于0．00005的所有组分都必须在计算中考虑。痕量组分(如c：H。等)应按表1中指定的赋值组分处理。所有组分的摩尔分数之和为l士o．0001。如果已知体积分数组成，则应将其换算成摩尔分数组成，具体换算方法见GB／T11062--1998。当压力和温度不以MPa和K表示时，必须将其分别换算成以MPa和K表示的值(有关换算因子觅GB3102．31993和GB3102．4—1993以及附录D)。表1微量和痕量组分一览表微量和痕量组分指定赋值组分02Ar、Ne、Kr、XeAmH2S一氧化二氨CO= 表1(续)GB／T17747．2—2011微量和痕量组分指定赋值组分氨CH．乙烯、乙炔、甲醇、氢氰酸czH6丙烯、丙二烯、甲硫醇C3H8丁烯、丁二烯、硫氧碳、二氧化硫nC⋯H新戊烷、戊烯、苯、环戊烷、二硫化碳n-C5H1zC。同分异构体、环己烷、甲苯、甲基环戊烷n—C6H【‘C，同分异构体、乙基环戊烷、甲基环己烷、环庚烷、乙苯、二甲苯n—C7Hl5C。同分异构体、乙基环己烷nCsHlIc。同分异构体H—CeH∞C。。同分异构体和更高碳数烃类”⋯CH224．4应用范围4．4．1管输气AGA8·92DC计算方法对管输气的应用范围如下：绝对压力：0MPa≤p≤12MPa热力学温度：263K≤丁≤338K高位发热量：30MJ·m13≤Hs≤45MJ·m_3相对密度：0．55≤d≤o．80洼：将本条中的高位发热量和相对密度换算为我国天然气标准参比条件下的高位发热量和相对密度，则高位发热量范围为27．95MJ·m-3～41．93MT·m，相对密度范围为0．550～0．800．天然气中各组分的摩尔分数应在以下范围以内：CH。O．7≤z(CH。)≤1．OON2O≤x(N2)≤O．20C020≤x(CO2)≤0．20CzH6O≤x(C2H6)≤0．10C3H8O≤x(C3H8)≤0．035C．Hl。0≤z(QHlo)≤O．015CsH12O≤x(C5H12)≤0．005C6HHO≤x(C6Hl‘)≤O．oolC7H。。O≤z(C7H16)≤0．0005C8+o≤z(cs+)≤O．0005H20≤_z(H2)≤O．103 GB／T17747．2—2011CO04z(CO)40．03He04z(He)40．005H：O04z(H20)40．00015所有摩尔分数大于0．00005的组分都不可忽略。微量和痕量组分见表l，并按指定的赋值组分处理。AGA8—92DC计算方法仅适用于单相气态(高于露点)混合物在操作压力和操作温度下压缩因子计算。4．4．2更宽的应用范围超出4．4．1所给出范围的应用范围如下：绝对压力：0MPa≤p≤65MPa热力学温度：225K≤T≤350K相对密度：0．55≤d40．90高位发热量：20MJ·m1≤Hs≤48MMJ·m_3注：将本条中的高位发热量和相对密度换算为我国天然气标准参比条件下的高位发热量和相对密度，则更宽的高位发热量范围为18．64MJ-m1～44．73MJ·m，相对密度范围为0．550--0．900。天然气中主要组分摩尔分数允许范围如下：CH‘0．50≤z(CH；)≤1_00N204z(N2)40．50C0204z(CO：)40．30C2H604z(C2H6)40．20C3H804x(C3H8)40．05H204z(H：)40．10管输气中微量和痕量组分含量范围见4．4．1。在超出以上范围应用时，AGA8—92DC方法的计算性能见附录E。4．5不确定度4．5．1管输气压缩因子计算的不确定度AGA8—92DC计算方法在4．4．1给出的管输气应用范围(温度为263K～350K，压力最大为12MPa)内，计算结果的不确定度为0．1％(见图1)。当温度高于290K，压力在最大为30MPa的范围内时，计算结果的不确定度也为0．1％。温度低于263K时，仅当压力在最高至10MPa的范围内，计算结果的不确定度才能保持在0．1％内。不确定度水平是通过将天然气压缩因子计算值与实验值数据库相比较而得到的(天然气压缩因子计算示例见附录c)。另外还同由称量法配制的模拟天然气混合物的压缩因子实验数据作了详细比较。用于试验本计算方法的两个数据库中实验测定值的不确定度在0．1％以内。 芒=爸GB／T17747．2—2011TfKAGA8-DC92方程P——压力}了’——温度}1——△Z≤0．1％；2——O．1％～O．2％；3——O．2％～O．5％．注：给出的不确定度范围仅适合于满足下面条件的天然气和类似气体：x(N2)≤O．20、x(C02)≤O．20、I(c2H6)≤0．10、x(H：)≤O．10、30MJ·m_3≤Hs≤45Ⅶ·m～、0．55≤丞≤0，80图1压缩因子计算的不确定度范围4．5．2更宽应用范围压缩因子计算的不确定度超出4．4．1给出气质范围的气体压缩因子计算的预期不确定度见附录E。4．5．3输入变量不确定度的影响表2列出的是相关输入变量的典型不确定度值，这些值可在最优操作条件下获得。根据误差传播分析，输入变量的不确定度会对压力为6MPa，温度在263K～338K范围内的压缩因子计算结果产生约0．1％的附加不确定度。当压力大于6MPa时，附加不确定度会更大，且大致与压力成正比例增加。表2相关输入变量的典型不确定度值输人变量绝对不确定度绝对压力o．02MPa温度0．15K情性组分摩尔分数0．001x(N2)0．001 GB／T17747．2—2011表2(续)输入变量绝对不确定度x(C02)0．001z(CH．)0．001z(C2H5)0．001x(C3H5)0．0005z(C‘-110H)0．0003x(C5_)0．0001z(Hz)和x(CO)o．0014．5．4结果的表述压缩因子和摩尔密度计算结果应保留至小数点后四位或五位，同时给出压力和温度以及所使用的计算方法(GB／T17747．2，AGAS-92DC计算方法)。验证计算机程序时，压缩因子计算结果应给出更多的位数。 符号口^BB二，b。附录A(规范性附录)符号和单位含义常数(表B．1)第二维利系数混合物交互作用系数[方程(B．1)和方程(B．2)]常数(表B．1)与温度和组成相关的系数组分i的特征能量参数(表B．2)组分J的特征能量参数第二维利系数的二元能量参数第二维利系数的二元能量交互作用参数混合物高温参数组分i的高温参数(表B．2)组分j的高温能量参数常数(表B．1)混合物定位参数组分i的定位参数(表B．2)组分j的定位参数二元定位参数二元定位交互作用参数(B．3)常数(表B．1)高位发热量体积参数组分i的体积参数(表B．2)组分J的体积参数二元体积交互作用参数(表B．3)常数(表B．1)摩尔质量组分i的摩尔质量气体混合物的组分数整数(1～58)绝对压力GB／T17747．2—2011单位m3·kmo[KMJ·m。3(m3／km01)1”7iTl3／km01)1“‘m3／km01)1／3kg·kmol一1MPa“CE马岛乓FRB^GGq岛瞄&魄KK巧％kM舰Nnp GB／T17747．2—20”四级参数组分i的四级参数组分J的四级参数常数(表B．1)气体常数(一o．008314510)组分i的偶极参数(表B．2)组分j的偶极参数常数(表B．1)绝对温度混合物能量参数混合物二元能量交互作用参数(表B．3)常数(表B．1)组分i的组合参数(表B．2)组分j的组合参数常数(表B．1)气体混合物中组分i的摩尔分数气体混合物中组分J的摩尔分数压缩因子质量密度气体的对比密度摩尔密度MJ·(km01·K)Kkg·m一3kraol·m3Qqt；RsS“Tu砜‰眠Ⅳ．‰五qzPp阳 B．1概述附录B(规范性附录)AGA8-92DC计算方法描述GB／T17747．2—2011使用4．2给出的AGA8—92DC方程计算气体混合物的压缩因子z。本附录将详细描述用AGA8—92DC方程进行压缩因子计算的有关方法和计算机执行程序，并给出必要的常数值。验证计算机程序用的压缩因子数据见附录c。如果计算机程序能够得到与附录C中数据相等的计算结果，则可使用。B．2AGA8—92Dc计算方法的计算机执行程序B．2．1计算程序概述B．2．1．1输入热力学温度T，绝对压力P和混合物中各组分的摩尔分数z，；B．2．1．2计算状态方程系数B和cj(n一13--58)，两者均取决于丁和z．；B．2．1．3利用改写的状态方程，迭代求解摩尔密度P。，以得到压力p；B．2．1．4当由B．2．1．3计算出的压力与B．2．1．1的输入压力，在规定的收敛范围内(如1×10“)相一致时，即得到压缩因子计算值。计算程序流程见图B．1。图B．1AGA9-92DC计算方法的程序流程图 GB／T17747．2—2011B．2．2计算步骤B．2．2．1输入热力学温度T，绝对压力P和混合物中各组分的摩尔分数z。；B．2．2．2根据B．2．2．1中输入的热力学温度丁和天然气组分摩尔分数z。，计算与组成和温度有关的状态方程系数B和c：(n一13～58)。第二维利系数B由方程(B．1)计算：混合交互作用系数B二由方程(B．2)计算：B二一(Gj+1一g_)h(QQ，+1一吼)k(F肛F72+1一^)‘(s；s，+1一s。)’一(w,wi+1一姒)、⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯(B．2)二元参数Ei和Gi由方程(B．3)和(B．4)计算：E#一E；(E。E，)“2⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯(B．3)G。一G：(Gt+G，)／2⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯(B．4)注B．1：除了表B．3中给出的目和G：外，所有其他二元交互作用参数耳和G：的值都是1．0．系数c：(n一13～58)由方程(B．5)计算：c：一n。(G+l—g。)‘(Q2+1一q。)％(F+1一，_)7-u～T一～⋯⋯⋯⋯(B．5)用共形求解混合方程(B．6)～(B．9)计算混合物参数u、G和Q；二重加和时，i从1～N—l变化，而相对每一个i值，j从i+1～N变化。Q一∑X。Q1)(E。E．)572⋯⋯⋯⋯⋯⋯(B．6)1)(G+G，)⋯⋯⋯⋯⋯⋯⋯(B．7)⋯⋯⋯⋯⋯⋯⋯⋯(B．8)⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯(B．9)注B．2：除了表R3给出的％、写、q和U,i外，所有其他的二元交互作用参数Kd、乓、q和U的值都为1．0。表＆2中H：的F(H)一1．0，而其他组分的Fi值都为0；水的w(H20)一1．0，而其他组分的Wi值都为0．B．2．2．3在压缩因子z的计算中，气体的组成z。、热力学温度T和绝对压力p都是已知的，问题在于要用表示压力P的状态方程来计算摩尔密度P。。将定义压缩因子z的方程(1)代人方程(3)(见4．2)，获得压力的状态方程(B．10)：1858P一阻RT[1+马D。一Pr∑c：+∑c：(6。一“五。P：一)砧-exp(一“P：t)]⋯⋯(B．10)^=13^一13方程(B．10)用标准状态方程密度检索法求解。由于已获得压力P的表达式[方程(B．10)3，则求解摩尔密度风使计算出的压力与输入的压力两者的差值在预先设定的范围(如1×10“)以内。对比密度B通过混合物体积参数K与摩尔密度P。相关联[见4．2中方程(2)]。混合物体积参数K由方程(B．11)计算：NN-1NK5一(∑z，K；72)2+2∑∑z．z，(K；一1)(K。K，)s／z⋯⋯⋯⋯⋯(B．11)I—It—l，一斗1注：求和时，下标i指的是气体混合物中第i个组分；下标J指的是气体混合物中第J个组分}N指的是混合物中的】0)K(～vE．～暮_B工z。∑H。∑㈣一下口。∑H—Bz。∑一¨∑Ⅲ+、)Ez。∑Ⅲ，【一U*一9@z工。∑=}M∑㈦+Gz。∑㈦一一GFz。∑Ⅲ一F GB／T17747．2—2011组分数。单重求和中，i是l～N间的整数值。例如，对含12个组分的混合物，N—12，单重求和中将有12项；二重求和中i从1～Ⅳ一1变化，而相对每一个i值，j从i+l～N变化。例如，对含12个组分的混合物，如果％的值都不为1．0，则二重求和将有66项。由于许多K。的值都为1．0，因此对许多天然气混合物，二重求和中非零项的数目很少。除了表B．3给出的K。外，所有其他的K。的值都为1．o．B．2．3求出摩尔密度m后，利用压力、温度、摩尔密度和摩尔气体常数计算压缩因子，见方程(B．12)。Z—p／(口。RT)⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯(B．12)密度p可由方程(B．13)计算：p=』Ⅵp。⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯(B．13)式中，M根据方程(B．14)计算：密度值应保留至小数点后第三位。表B．1状态方程参数6Ik，g-“f．1o．153832600100．O0O21．341953000100．50O03—2．998583000lO1．0O0O4一o．048312280lO03．50O0O050．37579650010O一0．510O6—1．589575000104．51O0O7—0．05358847010o．501O080．8865946301O7．50O1O9一O．710237040109．5O0O1O10—1．47l72200010O6．0O0l111．3218503501012．00112—0．78665925010O12．50O0O1132．291290×10—91l3—6．OO0l0O14o．1576724001l22．o0O015—0．436386400l123．0O0O0O16—0．044081590l122．00l0O17—0．003433888142．0O0O180．0320590501l411．00O0O190．02487355020—0．50O0O200．07332279020o．50O02l一0．00l6005732120．0O0O0O220．6424706002124．O0O23—0．4162601002l26．oO0O0O24J—O．0668995702l421．OOMz。∑㈨一M GB／T17747．2—201表B．1(续)巩^。g一幽f．250．27917950021423．01026—0．69660510021422．001O027—0．002860589214—1．00lO028—0．0080988363O0—0．5O1O0293．150547000317．010O0300．00722447931—1．001O031—0．7057529003126．0O0320．5349792003l24．01O0O33—0．0793149103131．01O34—1．4184650003139．O1O0O35—5．99905×10_17314—13．0O10O360．10584020031421．00O0O370．0343172903l48．001038—0．00702284740O～0．50390．024955870400．00O0400．0429681804122．00O410．7465453004127．00O0O42一O．2919613004l29．0O10O437．29461600041422．0O044—9．93675700041423．0O045—0．00539980850O1．0046—0．2432567005129．0O0O0470．0498701605l23．0010O480．0037337975148．0O0O0O491．87495100051423．0010O0500．00216814460O1．5O051—0．6587164006l25．01O0520．0002055187O一0．501O0530．0097761957124．00O054—0．020487080817．010O550．0155732208123．0O560．0068624158l20．01O57—0．00l2267529121．0O0O0580．0028509089120．00l0O12 表B．2特征参数GB／T17747．2—201摩尔质量能量参数体积参数识别化合物ME．K．定位参数四极参数高位参数偶极参数组合参数号G。QlF．S：W．kg·kmol一1K(m3·kmol1)1n1CH-16．0430151．3183000．46192550．02N228．013599．7377800．44791530．0278150．03CO：44．0100241．9606000．45574890．1890650．6900000．04C2H630．0700244．1667000．52792090．0793000．05C3H544．0970298．1183000．58374900．1412390．06H2018．0153514．0156000．38258680．3325001．0677500．01．5822001．07HzS34．0820296．3550000．46182630．0885000．6332760．00．3900000．08H22．015926．957940O．35149160．034369O．01．00．09CO28．0100105．5348000．45338940．0389530．010Ot31．998812Z．7667000．41869540．0210000．00．O0．011i—C‘Hlo58．1230324．0689000．64069370．2566920．012n—CIHlo58．1230337．6389000．63414230．2818350．013tC5H1272．1500365．5999000．67385770．3322670．014n—C㈣H72．1500370．6823000．67983070．3669110．015卅C6Hl‘86．1770402．6362930．71751180．2897310．016nC7Hl‘100．2040427．7226300．75251890．3375420．017卅C8H¨114．2310450．3250220．78495500．3833810．018卅C9H2口128．2580470．840891O．81527310．4273540．019卅C10H22142．2850489．5583730．84378260．4696590．020He4．00262．6101110．35898880．021Ar39．9480119．629900O．421655l0．0表B．3二元交互作用参数识别号化合物对彤UdK。《J12CH．+N2o．9716400．8861061．0036303C020．9606440．9638270．9959330．8076534C2H‘5C3Hl0．9946350．9908771．0076196H200．7082187H2S0．9314840．7368331．00008013 GB／T17747．2—2011表B．3(续)识别号化合物对耳UdKd《J8Hz1．1705201．1563901．0232601．9573109C00．99012610Oz11i—C‘Hlol_01953012水C‘HloO．9898440．99229l0．99759613{一C5Hlz1．00235014nCsHl20．9992681．0036701．00252915加C6H141_1072741．3025760．98296216n—C7H160．8808801．1919040．98356517肛c8H¨0．8809731．2057690．98270718n—c9H200．8810671．219634o．98184919卅CIoH220．8811611．2334980．98099123N=+CO：1．0227400．835058o．9823610．9827464C2H‘0．970120O．8164311．0079605C3Ha0．9459390．9155026H20o．7469547H。S0．9022710．9934760．9425968H21．086320o．4088381．0322709CO1．00571010021．021000110C4H100．94691412nC㈣HO．973384O．99355613i—C5H120．95934014n-C㈣H0．94552034CO：+C2H60．9250530．9698701．0085100．3702965C3H8o．9602376H200．8494081．6730907H2S0．9550521．0452901．0077908Hz1．2817909C01．5000000．900000100211二C‘Hl。0．90684912n-C㈣H0．89736214 表B．3(续)6B／"F17747．2—2011识别号化合物对E：U"K“G：J13i—C5H120．72625514n—C⋯Ho．85976415n—C⋯H0．8551341．0666380．91018316n—C7H“0．8312291．0776340．89536217卅GHla0．8083101．0881780．88115218n-C9H200．7863231．0982910．86752019”ⅦCH2z0．76517l1．1080210．85440645C2H6+C3HI1．0225601．065173O．9868936H200．6931687H2S0．9468710．9719260．9999698H21．1644601．6166601．0203409CO100211fC‘Hlo1．25000012n—aHlo1．0130601．250ooo13i—C5H121．25000014n—C㈣H1．0053201．25000058CaH。+H21．03478712n一巴Hl。1．004900715H2S+n—C6H1‘1．0086921．0289730．96813016卅GHl61．0101261．0337540．96287017月一C㈣H1．011501I．0383380．95782818n—C9H201．0128211．0427350．95244l19n—C㈣H21．0140891．0469660．9483388gH2+C01．100000lO021li—C‘Hlo1．30000012n-C㈣H1．30000015 GB／T17747．2—2011附录C(规范性附录)计算示例用AGA8号报告中描述的经过验证的计算机程序对下述示例进行了压缩因子计算，该经验证的计算机程序包含附录B所描述的子程序。表C．1以摩尔分数表示的气体组成分析数据气体组成1’气样24气样38气样48气样54气样6。气样x(C02)o．006o．005o．015o．016o．076o．011z(N2)0．003o．0310．0100．1000．0570．117z(I-12)O．OOO．095O．OOz(C0)0．00o．0100．00z(CH．)0．9650．9070．8590．7350．8120．826z(C2H‘)0．018o．04500．085o．0330．043o．035I(CaH^)0．00450．00840．0230．00740．0090．0075z(i-C4H10)0．00100．00100．00350．00120．00I50．0012x(n—C4Hlo)0．00100．00150．00350．00120．00150．0012z(i-C5H12)0．00050．00030．00050．00040．000．0004z(”一C5lll2H)0．00030．00040．00050．00040．000．0004z(C6Hl‘)0．00070．00040．000．00020．000．0002z(C7Hl6)0．000．00010．000．0001z(CoHll)0．000．00010．00表C．2压缩因子计算结果条件Pt14气样28气样3o气样4”气样5’气样64气样bar℃60一3．150．840530．833480．793800．885500．826090．85380606．850．86l990．855960．822060．901440．849690．873706016．850．880060．874840．845440．9150l0．869440．890526036．850．908670．904660．881830．936740．90052O．917235056．850．93011O．926960．908680．95318O．923680．93730120—3．150．721330．710440．641450．810240．695400．750741206．850．760250．750660．689710．837820．737800．7858612016．850．793170．784750．731230．861370．77369O．8156912036．850．845150．838630．796970．899130．830220．8631112056．850．883830．878700．845530．927660．872110．89862注：1bar=0．1～fLP& 附录D(规范性附录)压力和温度的换算因子GB／T17747．2—2011如果输入压力和温度的单位不是MPa和K，为了使用附录B描述的计算机执行程序进行计算，应作单位换算。部分换算因子如下：压力：p(MPa)一Ep(bar)]×10-1p(MPa)一D(atra)]X0．101325p(MPa)一p(psia)J／145．038p(MPa)一D(psig)+14．6959]／145．038温度：T(K)一f(℃)+273．15T(K)一[t(下)--32]／1．8+273．15丁(K)=Et(。R)]／1．817 GB／T17747．2—2011附录E(资料性附录)更宽范围的应用效果在温度263K～338K、压力最高至30MPa的范围内，利用实验测定值数据库，对AGA8—92DC计算方法进行了全面检验，这些实验数据取自给定组成范围内的管输气(见4．4．1)，在此组成范围内压缩因子计算的不确定度已在4．5中给出。对更宽范围(对组成而言)天然气(见4．4．2)压缩因子计算的不确定度的粗略估计，见图E．1～图E．4。图中以压力为纵坐标，横坐标分别为Nz、COz、CzHe和CsHe含量。图E．1～图E．4显示了在压力最高为30MPa的范围内，AGA8—92DC计算方法的应用效果。不确定度范围取决于压力、温度和组成，也强烈地受到相界位置的影响。图E．1～图E．4中给出的不确定度极限是在最坏情况下获得的结果，也就是说它们不是最佳选择。实验数据不足以决定不确定度范围的边界位置时，用虚线将所估计的不确定度区划为两个区域。气体全组成对相界位置会有强烈的影响，使用者应当进行相界计算。AGA8—92DC方程(丁=263K～338K)P——压力；XN：——Nz的摩尔分数#1——△Z≤O．1％；2——△Z0．1％～O．2％．图E．1计算高含N2天然气压缩因子时不确定度的估计范围 曼瓷GB／T17747．2—2011：。＼p——压力；z(COz)——cO=的摩尔分数；1——△Z≤0．1％；2——△Zo．1％～0．2％；3——△zo．2％～o．5％；4——△Zo．5％～3．o％。01j(c02)0’2AGA8-92DC方程(T=263K～338K)图E．2计算高含CO：天然气压缩因子时不确定度的估计范围AGA8—92DC方程(T=263K～338K)p——压力；z(c2H。)——GHs的摩尔分数；l——△Z≤O．1％；2——△Zo．1％no．2％；3——△Z0．2％～O．5％。图E．3计算高含c：116天然气压缩因子时不确定度的估计范围19 GB／T17747．2—2011皇≥一2一j＼＼＼00．05o．1叫C3Hs)AGA8-92DC方程(Tffi263K一338K)P——压力；x(CaHa)——csHs的摩尔分数；1——△Z≤O．1％；2——O．1％t0．2％。圉E．4计算高含G珥天然气压缩因子肘不确定度的估计范围压力最高至10MPa，温度在263K～338K时，压缩因子计算的综合结果总结如下：在给定的压力温度范围之内，只有组分摩尔分数在表E．1列出值范围内的气体，其压缩因子计算的不确定度才分别在0．1％，0．2％和0．5％以内。表E．1AGA8—92DC计算方法的计算综合结果不确定度范围内的摩尔分数组分o．I％0．2％0．5％Nz≤0．50CO=≤0．23≤0．26≤0．28c2H。≤0．13≤0．20c3H8≤0．06≤O．10'