天然气基本压缩因子计算方法

合集下载

一种新型天然气压缩因子数值计算方法

0. 99
1. 58
2. 59
4. 30
6. 00
9. 33
由表 3 可以看出，在对比温度为 1. 05 时，本文
计算方法与其他计算方法一样误差较大，最大误差
达到了 60% 。从压缩因子三维图可以看出，对比
温度在 1. 05 ～ 1. 10 之间时，Z 值曲面表现出了较
强的扭曲性，这也是造成各方法预测精度均较低的
近几年国内相继发现了一批高温高压天然气田［8］。原有压缩因子计算方法适用压力范围低的弊端逐渐暴露。石油大学郭绪强教授针对这一问题进行了相关试验，取得了丰富的高压天然气实验数据［9］。
将郭绪强教授发表的高压天然气实验数据与传统天然气压缩因子图版［10］叠加，发现天然气压缩因子在高压阶段具有较强的延展性，表现出了较好的规律。利用三维绘图软件将数据进行处理，可以发现天然气压缩因子曲面较为复杂。因此，本文利用传统压缩因子图版与郭绪强教授发表的高压天然气实验数据进行拟合，尝试找到高精度的能够
引言
目前使用较多的天然气压缩因子计算方法，包括 Dranchk － Abu － Kassem 方法（ DAK）［1 － 2］，Hankinson － Thomas － Phillips 方法（ HTP ）［3］，Dranchuk － Purvis － Robinson 方法（ DPR）［4］，以及由石油大学李相方教授根据天然气压缩因子图版拟合的李相方方法（ LXF）。这些计算公式均是根据 Standing 和 KatZ 1942 年提出的压缩因子图版［5］采用不同拟合方法拟合得到的［6］。在不同的对比压力及对比温度下，误差均较大。根据李相方教授的统计，各方法的最大误差均超过了 55%［7］。

AGA8—92DC计算方法天然气压缩因子计算(最漂亮的)

AGA8—92DC计算方法天然气压缩因子计算摘要：按照GB／T 17747．2—1999《天然气压缩因子的计算第2部分：用摩尔组成进行计算》，采用AGA8—92DC计算方法，用VB编程计算了天然气压缩因子。

用二分法求解状态方程，精度满足工程需要。

关键词：压缩因子；AGA8—92DC计算方法；二分法1概述工作状态下的压缩因子是天然气最重要的物性参数之一，涉及到天然气的勘探、开发、输送、计量和利用等各个方面。

实测天然气压缩因子所需的仪器设备价格高，不易推广，因此计算方法发展很快，主要为经验公式和状态方程计算方法。

1992年6月26日，国际标准化组织(ISO)天然气技术委员会(TC193)及分析技术分委员会(TC193/SC1)在挪威斯泰万格(Stavanger)召开了第四次全体会议，会上推荐了两个精度较高的计算工作状态下天然气压缩因子的方程，目PAGA8-92DC方程、SGERG-88方程[1]。

随后，国际标准化组织于1994年形成了国际标准草案[2]。

AGA8-92DC方程来自美国煤气协会(AGA)。

美国煤气协会在天然气压缩因子和超压缩因子表的基础上，开展了大量研究，于1992年发表了以状态方程为基础计算压缩因子的AGA No.8报告及AGA8-92DC方程[2]。

1994年，四川石油管理局天然气研究所遵照中国石油天然气总公司技术监督局的指示，对国际标准化组织1992年挪威斯泰万格会议推荐的AGA8-92DC 方程、SGERG-88方程进行验证研究，于1996年底基本完成[2]。

1999年，四川石油管理局天然气研究院(前身为四川石油管理局天然气研究所)起草的《天然气压缩因子的计算》GB/T 17747.1～3—1999被批准、发布。

《天然气压缩因子的计算》GB/T 17747.1～3—1999包括3个部分：《天然气压缩因子的计算第1部分：导论和指南》GB/T 17747.1—1999，《天然气压缩因子的计算第2部分：用摩尔组成进行计算》GB/T 17747.2—1999，《天然气压缩因子的计算第3部分：用物性值进行计算》GB/T 17747.3—1999。

天然气压缩因子的分析及其计算

天然气压缩因子的分析及其计算谢莉莉;刘劲松【摘要】根据天然气压缩因子的2种计算方法:用摩尔组成进行计算和用物怀值进行计算编制计算机程序,并运用此程序研究天然气压缩因子与温度、压力之间的关系.【期刊名称】《上海计量测试》【年(卷),期】2011(000)005【总页数】4页(P27-30)【关键词】天然气;压缩因子;计算方法【作者】谢莉莉;刘劲松【作者单位】上海公正燃气计量站;上海公正燃气计量站【正文语种】中文0 引言天然气是重要的能源之一，随着天然气贸易量的增加，其流量计量越来越被人们重视。

在天然气流量计量中，天然气压缩因子是决定其准确与否的关键因素之一。

天然气压缩因子是实际气体状态采用理想气态方程时引入的偏差修正系数。

实际上，符合理想气态方程的理想气体是不存在的，实验表明，只有在低压高温下实际气体才可以近似被看作理想气体。

由于实际气体与理想气体的差异，使得对气体流量测量的准确性和可靠性难以评价，特别是低温、高压管道气体流量的测量，在这种情况下，管道中的被测介质就不能用理想气体状态方程进行描述。

在高压、低温下，任何气体理想状态方程都会出现明显的偏差，而且压力越高，温度越低，这种偏差就越大，因而需要引入一个压缩校正因子Z来修正气体的状态方程，如式（1）所示。

因此，天然气压缩因子Z在天然气这一重要能源计量中起着举足轻重的作用。

虽然GB/T 17747-1999《天然气压缩因子的计算》对天然气压缩因子进行了详细的描述，但是国内大部分是使用超压缩因子来计算天然气流量，对于压缩因子大多是文献上查得的或是通过图表获得。

若是用图表方式，则整个计算过程不仅费时费力，而且计算误差大，结果不准确。

而国外的进口流量计，像压缩因子等技术核心不公开，因此有必要编制一套计算程序来计算天然气压缩因子，确保天然气流量计量的准确性。

本文将介绍程序的编制简要以及运用该程序研究压缩因子与温度、压力之间的关系，并对两种方法进行比较。

1 计算程序编制天然气压缩因子的计算方法有2种：用天然气的摩尔组成进行计算和用天然气的物性值进行计算。

AGA8—92DC计算方法天然气压缩因子计算（最漂亮的）

AGA8—92DC计算⽅法天然⽓压缩因⼦计算（最漂亮的）AGA8—92DC计算⽅法天然⽓压缩因⼦计算摘要：按照GB／T 17747．2—1999《天然⽓压缩因⼦的计算第2部分：⽤摩尔组成进⾏计算》，采⽤AGA8—92DC计算⽅法，⽤VB编程计算了天然⽓压缩因⼦。

⽤⼆分法求解状态⽅程，精度满⾜⼯程需要。

关键词：压缩因⼦；AGA8—92DC计算⽅法；⼆分法1概述⼯作状态下的压缩因⼦是天然⽓最重要的物性参数之⼀，涉及到天然⽓的勘探、开发、输送、计量和利⽤等各个⽅⾯。

实测天然⽓压缩因⼦所需的仪器设备价格⾼，不易推⼴，因此计算⽅法发展很快，主要为经验公式和状态⽅程计算⽅法。

1992年6⽉26⽇，国际标准化组织(ISO)天然⽓技术委员会(TC193)及分析技术分委员会(TC193/SC1)在挪威斯泰万格(Stavanger)召开了第四次全体会议，会上推荐了两个精度较⾼的计算⼯作状态下天然⽓压缩因⼦的⽅程，⽬PAGA8-92DC⽅程、SGERG-88⽅程[1]。

随后，国际标准化组织于1994年形成了国际标准草案[2]。

AGA8-92DC⽅程来⾃美国煤⽓协会(AGA)。

美国煤⽓协会在天然⽓压缩因⼦和超压缩因⼦表的基础上，开展了⼤量研究，于1992年发表了以状态⽅程为基础计算压缩因⼦的AGA No.8报告及AGA8-92DC⽅程[2]。

1994年，四川⽯油管理局天然⽓研究所遵照中国⽯油天然⽓总公司技术监督局的指⽰，对国际标准化组织1992年挪威斯泰万格会议推荐的AGA8-92DC ⽅程、SGERG-88⽅程进⾏验证研究，于1996年底基本完成[2]。

1999年，四川⽯油管理局天然⽓研究院(前⾝为四川⽯油管理局天然⽓研究所)起草的《天然⽓压缩因⼦的计算》GB/T 17747.1～3—1999被批准、发布。

《天然⽓压缩因⼦的计算》GB/T 17747.1～3—1999包括3个部分：《天然⽓压缩因⼦的计算第1部分：导论和指南》GB/T 17747.1—1999，《天然⽓压缩因⼦的计算第2部分：⽤摩尔组成进⾏计算》GB/T 17747.2—1999，《天然⽓压缩因⼦的计算第3部分：⽤物性值进⾏计算》GB/T 17747.3—1999。

压缩因子计算

天然气压缩因子的计算气田上大多数在高压下生产，为控制其流动需要安装节流阀。

当气流经过节流阀时，气体产生膨胀，其温度降低。

如果气体温度变得足够低，将形成水合物（一种固体结晶状的冰雪物质）。

这就会导致管道和设备的堵塞。

【1】从而，在天然气的集输过程当中，不管对天然气或天然气管道进行怎样的处理，都离不开气体的三个状态参数：压力P 、体积Ｖ、温度Ｔ。

而根据真实气体状态方程PV ZnRT =可知，在确定某个状态参数的时候需要先计算一个压缩因子Ｚ。

如果能够更精确的确定压缩因子，从而确定气体的状态参数，对于研究天然气的收集、预处理和输送等问题具有重要意义。

下面简要介绍下压缩因子及其计算方法。

真实气体是实实在在的气体，它是为了区别于理想气体而引人的。

真实气体占有一定空间，分子之间存在作用力，因此真实气体性质与理想气体性质就有偏离。

压缩因子就是反映这种真实气体对理想气体的偏离程度大小。

在温度比临界温度高的多、压力很小时，偏离不太显著；反之偏离就很显著。

下面将介绍一种计算压缩因子的方法（Dranchuk-Purvis-Robinson 法）。

压缩因子的关系式如下：5635214373831()()()(1)exp()pr pr pr pr pr A A A A A Z A A T T T T A A A T =++++++++-52pr pr pr 222prpr pr ρρρρρρ （1）式中A 1到A 8都是常数，具体数据可到参考文献上查阅，ρpr 为无因次拟对比密度，它和压缩因子满足关系式： 0.27prpr pr p ZT ρ= （2）其中p pr 和T pr 分别为拟对比压力和拟对比温度。

由于式（2）为非线性方程，欲计算Z ，可采用牛顿迭代法（Newton-Raphson ）。

在已知p pr 和T pr 的情况下，需经过迭代过程求解ρpr ，其公式如下：()(1)()'()()()i pr i i pr pr i pr f f ρρρρ+=- （3）迭代求得拟对比密度ρpr ，即可易求得压缩因子。

天然气压缩因子计算

１．天然气相关物性参数计算密度计算： TZR PM m =ρ ρ——气体密度，Kg/m 3；P ——压力，Pa ；M ——气体千摩尔质量，Kg/Kmol ；Z ——气体压缩因子；T ——气体温度，K ；R m ——通用气体常数，8314.4J/Kmol·K 。

２．压缩因子计算：已知天然气相对密度∆时。

96.28M =∆ M ——天然气的摩尔质量。

∆+=62.17065.94pc T510)05.493.48(⨯∆-=pc P ;pc pr P P P = pcpr T T T =； P ——工况下天然气的压力，Pa ；T ——工况下天然气的温度，k ；P Pc —临界压力；T Tc ——临界温度。

对于长距离干线输气管道，压缩因子常用以下两式计算：668.34273.01--=prpr T P Z 320107.078.068.110241.01prpr pr pr T T T P Z ++--=对于干燥天然气也可用经验公式估算： 15.1117.0100100P Z +=标况流量和工况流量转换。

为了控制Welas 的5L/min 既 0.3立方米每小时的工况流量。

Q 2------流量计需要调节的流量值P 2------0.1MpaT 2------293.15K (20℃ )Z 2------标况压缩因子Q 1------0.3m 3/hP 1------ 工况压力（绝对压力MPa ）T 1------开尔文KZ 1-------工况压缩因子转换公式为12221211p T Z Q Q p T Z。

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

相对密度：
０ａ镇６ＭＰＭＰ＜ｐ５ａ２５＜３０２Ｋ＜Ｔ５Ｋ２ＭＪｍ－０・＇＜４ＭＪ＇＜Ｈ，８＂ｍ－
注２将本条中的高位发热量和相对密度换算为我国石油气体标准参比条件下的高位发热量和相对密度，：则高位发热１范围为７９－４．＂，２．１９ＭＪ＇相对密度范围为０５０－．，１５３ｍ－．－０８０５０
天然气中各组分的摩尔分数应在以下范围以内：
— 第３部分：用物性值进行计算。
附录Ａ、附录Ｂ附录Ｃ附录Ｄ是标准的附录。附录Ｅ附录Ｆ附录Ｇ是提示的附录。、、、、
中华人民共和国国家标准
天然气压缩因子的计算第２部分：用摩尔组成进行计算
Ｎｔａｇｓａｃａｏｏｏｒｓｉｆｃｏ－ｒｌ－ＣｌｌｉｎｃｍｐｅｓｎｔｒａｕａｕｔｆｏａＰｒ２Ｃｌｌｉｎபைடு நூலகம்ｉｇｌｒｏｐｓｉｎａｓａｔａｕｔｕｎｍａ－ｍｏｉｏａｌｉ：ｃａｏｓｏｃｔｎｙｓ
似法是否会使计算结果变差．
摩尔分数大于０００．００５的所有组分都必须在计算中考虑。痕量组分（Ｃ等）如Ｍ；应按表１中指定的赋值组分处理。所有组分的摩尔分数之和为１．１士００００，
如果已知体积分数组成，则应将其换算成摩尔分数组成，具体换算方法见Ｇ／１０２ＢＴ６０１
ｃ／１７７２９９ＳＴ４．－１９７
绝对压力：
热力学温度：
０ａ－１ＭＰＭＰ＜ｐ２ａ＜
２３６Ｋ镇Ｔ蕊３８３Ｋ

天然气压缩因子的计算

天然气压缩因子的计算天然气的压缩系数计算方法可采用GB/T 17747-1999《天然气压缩因子的计算》,或AGA NX-19方程。

当为非贸易计量场合和贸易计量中符合GB/T 18603-2001《天然气计量系统技术要求》表A1准确度为C 级要求的计量装置可考虑使用AGA NX-19方程,其它应采用SGERG-88或AGA 8-92DC 方程。

本文描述AGA NX-19和SGERG-88两种计算方法。

1.用物性值进行计算天然气压缩因子的公式本计算公式参照国家标准GB/T 17747中SGERG-88公式,该计算公式使用高位发热量、相对密度和CO 2含量作为输入变量。

在GB/T 17747中,用物性值计算天然气压缩因子公式如下:21mm C B Z ρρ++= (1 /(ZRT p m =ρ (2式中有关符号表示见本文后述的符号说明。

天然气压缩因子Z 的值由方程(1、(2联解求得(1式中:天然气第二维利系数B 由方程(B1求得B x x B x B 2111212+=12++++++233222255254424332324422B x B x B x B x x +++ (B1(B1式中:CH x x =1 (B2 22N x x = (B3 23CO x x = (B4 24H x x = (B5 CO x x =5 (B6 2(1(0([2(1(b 0(11120H0011H H H H H b T b b T b T b B +++++=CH H T ]222222]2(1(0([CH H H H H T b T b b +++ (B7B 14,B 15,B 22,B 23,B 24,B 33, B 44和B 55是温度函数的二次多项式,即:2-12B B B +×+= (B91/23311130.865(B -B B = (B10(B7和(B8式中维利系数温度展开式系数b(0,b(1和b(2的数值 ij b(0 b(1 b(2 CH H0 -4.25468×10-1 2.86500×10-3 -4.62073×10-1CH H1 8.77118×10-4 -5.56281×10-6 8.81510×10-9 CH H2 -8.24747×10-7 4.31436×10-9 -6.08319×10- 12 N 2 22 -1.44600×10-1 7.40910×10-4-9.11950×10-7CO 2 33 -8.68340×10-1 4.03760×10-3 -5.16570×10-6H 244-1.10596×10-3 8.13385×10-5-9.87220×10-8 CO 55 -1.30820×10-1 6.02540×10-4 -6.44300×10-7CH+ N 2 12 y =0.72+1.875×10-5(320-T2 CH+ CO 2 13 y =-0.865CH+ H 2 14 -5.21280×10-2 2.71570×10-4-2.50000×10-7CH+ CO 15-6.87290×10-2-2.39381×10-6 5.18195×10-7N 2+ CO 2 23 -3.39693×10-1 1.61176×10-3 -2.04429×10-6N 2+ H 2 24 1.20000×10-2 0.00000 0.00000天然气第二维利系数C 由方程(C1求得1233211222444343333323323222332222232133231333C x C x C x x C x x C x C x x ++++++ (C1 (C1式中:211120001112(1(0([2(1(0(T c T c c T c T c c C H H H H H H +++++=CH H ]22222]2(1(0([CH H H H H T c T c c +++ (C2C 222,C 333,C 444,C 113,C 223和C 233是温度函数的二次多项式,即:22(1(0(T c T c c C ijk ijk ijk ijk ++= (C3维利系数温度展开式中系数c(0,c(1和c(2的数值ijk c(0 c(1 c(2CH H0 -3.02488×10-1 1.95861×10-3 -3.16302×10-6CH H1 6.46422×10-4-4.22876×10-6 6.88157×10-9CH H2 -3.32805×10-7 2.23160×10-9 -3.67713×10-12N 2 222 7.84980×10-3-3.98950×10-5 6.11870×10-8CO 2 333 2.05130×10-3 3.48880×10-5 -8.37030×10-8H 2 444 1.04711×10-3-3.64887×10-8 4.67095×10-9CH+ CH+ N 2 112 y =0.92+0.0013(T-270 CH+ CH+ CO 2 113 y = 0.92 CH+ CH+ H 2 114 y = 1.20 CH+ CH+CO1157.36748×10-3-2.76578×10-5 3.43051×10-8CH+ N 2+ N 2 122 y =0.92+0.0013(T-270 CH+ N 2+ CO 2 123 y =1.10 CH+ CO 2+ CO 2 133 y =0.92 N 2+ N 2+ CO 2 223 5.52066×10-3 -1.68609×10-5 1.57169×10-8N 2+ CO 2+ CO 2 2333.58783×10-3 8.06674×10-6 -3.25798×10-8其他非同类交互作用维利系由方程(C4求得: ijk ijk y C =3/1(kkk jjj iii C C C (C4(C4式中ijk y 由(C5~(C8给出:92.0133113==y y (C620.1114=y (C710.1123=y (C8式中符号:H ——摩尔发热量,单位:MJ ·kmol -1x ——组分的摩尔分数CH ——等价烃类 CO ——一氧化碳 CO 2——二氧化碳 H 2——氢气 N 2——氮气m ρ——摩尔密度,单位: kmol -1·m 3p ——绝对压力,单位:MPaR ——摩尔气体常数,其值为0.008314510 m 3·kmol -1K -1 T ——热力学温度,单位:K2.用AGA NX-19公式计算天然气压缩因子的方法天然气超压缩系数Fz 是因天然气特性偏离理想气体定律而导出的修正系数,其定义为1Z ZnFz =………………………………………………………(3 式中: Zn —天然气在标准参比条件下的压缩因子;Z1 —天然气在操作条件下的压缩因子。

天然气压缩因子计算及影响因素分析

天然气压缩因子计算及影响因素分析王春生;徐玉建;田明磊;董国庆;徐畅;陈钊【摘要】Measurement shortage will often arise between the head and the end of nature gas pipeline which is a vital important influencing factor of transmission cost. Regarding to the phenomenon of measurement shortage, we focused on the compressibility factor and tried to solve the problem by optimizing the calculation method of the compressibilityfactor so that the phenomenon can be well control. On the basis of BWRS equation, first equation coefficients were obtain by Excel, then the gas density was calculated with these coefficients, finally all these results were put into the gas state equation to obtain the compressibility factor. By solving the gas compressibility factor, its main influencing factors were determined, which could help to correct the throughput of natural gas to keep measurement shortage to the minimum.%天然气长输管道首端与末端之间往往会出现输差，输差是影响输气成本的一个最关键的因素。

天然气流量计量中确定压缩因子的方法

Ｔ。
Ｐ
。
量的惰性气体。由于其热值高、价格低廉、安全性
好、燃烧物对环境污染少，被公认是优质的节能环保燃料 … 。如今天然气不仅作为居民生活燃料，
同时被广泛应用于工业生产、交通运输及航
气体，其拟临界温度、压力的定义为：
ｌ２
化
工
自动
化及
仪
表
第４Ｏ卷
天然气流量计量中确定压缩因子的方法
李言伟
（福建省安然燃气投资有限公司．福州３５００００）
摘要以天然气流量计量为背景，对天然气压缩因子的修正进行了分析，并基于数理统计方法对压缩
。 —
ｚｃ ‘ Ｔ
Ｒ—— 通用气体常数；Ｐ —Leabharlann 对比体积、温度和压力；
ｚ一
．
婷
临界压缩因子。
根据对比态的原理，真实气体在相同的对比
压力和对比温度条件下，应当具有相同的对比体积，所以压缩因子与对比状态参数应遵循以下关
２．２曲面拟合
为：设定实验数据（，Ｙ，ｚ）（ｉ＝１，２，３， …，ｎ），寻找函数，（，Ｙ，三）使得在点（，Ｙ）（＝１，２，３， …，ｎ）处的函数值与观测数据偏差的平方和最小，即求得满足误差项平方和ＳＳＥ最小的函数，（，Ｙ，

用于计量的天然气压缩因子计算方法比较[1]

N G03 0. 25 0. 60 96. 50 1. 75 0. 40 0. 10 0. 10 0. 10 0. 10 0. 10
N G04 0. 56 0. 52 94. 53 0. 96 1. 55 0. 30 0. 79 0. 22 0. 19 0. 24 0. 14
N G05 N G062)
第 20 卷第 5 期天然气工业集输工程
用于计量的天然气压缩因子计算方法比较
张福元3
(西南油气田公司天然气计量检测中心)
张福元. 用于计量的天然气压缩因子计算方法比较. 天然气工业 ,2000 ;20 (5) :73～76 摘要天然气压缩因子或超压因子计算结果的准确性直接影响天然气流量计量的准确性。当前国内天然气计量界广泛使用 A GANX219〔1〕,A GA8 号报告 , ISO 1221321997 三种天然气压缩因子计算方法标准。文章研究了这三种天然气压缩因子计算方法标准 ,并编写了 N GZCWIN 天然气压缩因子计算软件 ,通过对不同气样和不同温度、压力条件的计算 ,比较了三种计算方法的差别 ,并对这些计算方法的应用范围和不确定度提出了看法。主题词天然气计量压缩系数计算方法分析
表 2 A GA8 号报告和 ISO 12213 的适用范围
项目
A GA8 号报告
ISO 12213
管输范围扩展范围管输范围扩展范围
压力 (MPa) 0～12 0～280 0～12 0～651)
温度 ( ℃) - 8～65
相对密度
0. 554 ～0. 87
高位发热量 18. 7 (MJ/ m3) ～45. 1
天然气压缩因子计算方法简介
1. A GA8 号报告在 A GA8 号报告 1994 年版中〔2〕,提供了以组成

天然气压缩因子计算及影响因素分析

天然气压缩因子计算及影响因素分析作者：王春生徐玉建田明磊董国庆徐畅陈钊来源：《当代化工》2015年第06期摘要：天然气长输管道首端与末端之间往往会出现输差，输差是影响输气成本的一个最关键的因素。

针对出现的输差问题，以天然气组分为基础，以压缩因子作为突破口，通过着重理解天然气压缩因子的解法与改进来得到控制输差。

以BWRS方程为重点，通过Excel求得方程系数，然后从中解出气体密度，再代入气体状态方程中求得压缩因子。

通过对天然气压缩因子的求解，得到影响压缩因子的主要因素，从而修正到天然气输量，以便减少输差。

关键词：输差;压缩因子;BWRS方程;影响因素中图分类号：TQ 018 文献标识码： A 文章编号： 1671-0460（2015）06-1408-04Calculation of Natural Gas Compressibility Factor and Its Influence FactorsWANG Chun-sheng1，XV Yu-jian1，TIAN Ming-lei1，DONG Guo-qing1，XV Chang1，CHEN Zhao2（1. Northeast Petroleum University， Heilongjiang Daqing 163318，China;2. 7th China Petroleum Engineering&Construction Corp， Shandong Qingdao 266300，China）Abstract： Measurement shortage will often arise between the head and the end of nature gas pipeline which is a vital important influencing factor of transmission cost. Regarding to the phenomenon of measurement shortage， we focused on the compressibility factor and tried to solve the problem by optimizing the calculation method of the compressibility factor so that the phenomenon can be well control. On the basis of BWRS equation， first equation coefficients were obtain by Excel， then the gas density was calculated with these coefficients， finally all these results were put into the gas state equation to obtain the compressibility factor. By solving the gas compressibility factor， its main influencing factors were determined， which could help to correct the throughput of natural gas to keep measurement shortage to the minimum.Key words： Measurement shortage; Compressibility factor; BWRS equation; Influence factor天然气与其他能源材料，例如煤炭和石油相比-天然气的热值较高，利用率较高，并且对环境的污染很小。

GB T 17747.1-1999 天然气压缩因子的计算第1部分：导论和指南

积的比值：
ｚ＝ｖ彩真实）ｖ／彩理想）・・・・・・・・……（・・・・・・・・・・・・・・・・１）Ｚ（，ｙ＝ｐ｝真实）（Ｔ）・・・・・・・・・……（３）ｐＴ，）Ｖ，（／Ｒ・・・・・・・・・・・・・・・・・・
为国际标准正式发布。国际标准Ｓ１１－是由天然气技术委员会ＩＯＴ１３ＩＯ２３１２Ｓ／Ｃ下的“ ９天然气分析’ ，分委员会制定的。Ｉ１２３天然气压缩因子的计算》《Ｓ２１Ｏ标准包括以下３个部分： — 第１部分：导论和指南； — 第２部分：用摩尔组成进行计算； — 第３部分：用物性值进行计算。附录Ａ是标准的附录。附录Ｂ和附录Ｃ是提示的附录。
ｚ无量纲， — 压缩因子，值通常接近于１。
３２密度ｄｎｉ，－ｅｓｙｐｔ
见Ｇ／１６－１９１０２９８中２３ＢＴ．．
３３摩尔．组成ｍｌｃｐｓｉｏｒｏｉｎａｏｍｔｏ
用摩尔分数或摩尔百分数表示的均匀混合物中每种组分的比例。给定体积的混合物中１组分的摩尔分数２是１ ‘ 组分的摩尔数与混合物中所有组分的总摩尔数（即所有组分摩尔数之和）之比。摩尔任何化合物所含物质的量等于以克为单位的相对摩尔质量。１相对摩尔质量的推荐值见Ｇ／１０２ＢＴ６，１对于理想气体，摩尔分数或摩尔百分数与体积分数或体积百分数值完全相等。对真实气体，两者一
产物的温度降至与规定的反应物温度ｔ，相同的温度，并且除燃烧生成的水在温度ｔ下全部冷凝为液态，
外，其余所有燃烧产物均为气态。高位发热量包含天然气中所有可燃组分。
燃烧参比条件：温度ｔ为２８１Ｋ２０）压力Ｐ为１１３５，，９．（Ｃ，５５，０．ｋａ２Ｐ体积计量参比条件：温度ｔ为２３１Ｋ０）压力Ｐ为１１３５，：７．５０，（Ｃ｝０．ｋａ２Ｐ

压缩因子计算方法

压缩因子计算方法
汇报人：齐少鹏
压缩因子简介
• 由于理想气体作了两个近似：忽略气体分子本身的体积和分子间的相互作用力，所以实际气体都会偏离理想气体。
• 压缩因子Z被引用来修正理想气体状态方程:PV=nRT 。
• 压缩因子的定义式为：Z=PV/nRT ，压缩因子的量纲为一。
• 很显然，Z的大小反映出真实气体对理想气体的偏差程度，即Z等Vm（真实）除以Vm（理想）。由于Z反映出真实气体压缩的难易程度，所以将它称为压缩因子。
Tp3r
(1.18 2.82 ) Tpr
r
• 特殊定义的对比密度：
r

0.06125
p pr ZTpr
exp1.2(1

1 Tpr
)2

• 方法的适用范围：Tpr>1
AGA公式
•
美国加利福尼亚天然气协会（CNGA）公式前苏联气体研究所公式
低压下压缩因子的确定
P<35MPa
• Z=1+(0.31506-1.0467/Tpr-0.5783/T3pr)ρ +pr
(0.5353-0.6123/Tpr)ρ
2pr+0.6815ρ
/T 2pr
3pr
•
ρ
=0.27
p /ZT pr
pr
• • 迭代
RK公式
RK(Redlich—Kwong)方程是 1949 年提出的二参数状态方程，它在范德瓦尔斯基础上引入了温度对引力的修正，多用于计算干气的压缩因子，不适于计算非极性分子的压缩因子
实验方法求取天然气压缩因子
实验测定天然气压缩因子方法是将一定质量的天然气样品装入高压物性实验装置的PVT筒中,在恒温条件下测定天然气的压力与实际体积V的关系。

matlab基于bwrs方程的天然气压缩因子计算

一、概述天然气压缩因子是表征天然气在不同压力和温度条件下物性的重要参数之一。

而计算天然气压缩因子则是石油工程中的一个重要问题。

matlab基于bwrs方程的天然气压缩因子计算是使用matlab编程语言结合bwrs方程进行天然气压缩因子计算的一种方法。

本文将介绍该方法的原理和实现过程。

二、天然气压缩因子计算原理天然气的压缩因子是一个无量纲数，用于描述天然气在不同压力和温度下的实际体积与理想气体体积之间的关系。

天然气压缩因子随着压力和温度的变化而变化，因此在石油工程中，需要对天然气的压缩因子进行精确的计算。

bwrs方程是一种常用的天然气压缩因子计算方法，其表达式如下：Z = 1 + (B1 - B2/T + B3/T^3)P + (B4/T^2 - B5/T^3)P^2 +B6/T^2P^3其中，Z为天然气压缩因子，P为压力，T为温度，B1、B2、B3、B4、B5、B6为实验参数。

通过bwrs方程，可以利用给定的压力和温度条件，计算出相应的天然气压缩因子，为后续天然气输运和储存等工程设计提供重要的物性参数。

三、matlab编程实现1. 输入压力和温度在matlab中，首先需要输入天然气的压力和温度，作为bwrs方程的输入参数。

这可以通过matlab的命令窗口或者在编写的程序中实现。

2. 编写bwrs方程计算函数接下来需要编写一个函数，用于实现bwrs方程的计算。

可以通过编写一个独立的.m文件，命名为bwrs_equation.m，其中包含bwrs方程的计算逻辑。

3. 计算压缩因子在bwrs_equation.m中，利用输入的压力和温度参数，结合bwrs 方程的表达式，计算出对应的天然气压缩因子。

这包括了对bwrs方程中的系数B1至B6进行定义和赋值，以及根据输入的压力和温度计算出压缩因子Z。

4. 输出计算结果将计算得到的天然气压缩因子Z输出到matlab的命令窗口或者保存到文件中，以便后续的使用和分析。

天然气基本压缩因子计算方法

天然气基本压缩因子计算方法编译:阙洪培(西南石油大学审校:刘廷元这篇文章提出一个简便展开算法:任一压力-温度的基本压缩因子的输气监测计算。

这个算法中的二次维里系数来源于参考文献1。

计算的压缩因子接近AGA 8状态方程值[2]。

1 测量在天然气工业实用计量中,压力、温度变化作为基本(或标准条件,不仅地区间有差别,而且在天然气销售合同也有不同。

在美国,通常标准参考条件是60°F和14.73 psia。

欧洲常用的基本条件是0 ℃和101.325 kPa,而标准条件是15 ℃和101.325 kPa。

阿根廷也用15 ℃和101.325 kPa,而墨西哥则用的是20 ℃和1kg/ sq cm(绝对。

计算真实气体的热值、密度、基本密度、基本体积、以及沃贝指数时要求已知基本条件的压缩因子。

表1是理想气体值。

表1中的理想气体值不能用于密闭输气,必须计算相应基本条件的压缩因子。

参考文献提供的一些数据表和获取基本条件压缩因子方法,基本条件只能是60°F,14.73或14.696 psia。

计算其它基本条件的压缩因子可用AGA 8 程序,但代数计算较复杂,计算机编程共有三组软件,比较耗时。

本文提出了一个展开算法,计算密闭输气基本条件(基本条件可是任何压力温度的压缩因子。

2 压缩因子接近外界条件时,即压力小于16 psia,截断维里状态方程(方程组中的方程1较好地描述了天然气的体积性质。

方程1中,各符号的物理意义是:Z = 基本条件下压缩因子B = 二次维里系数R = 气体常数P = 基本条件的绝对压力T = 温度条件的绝对压力天然气基本压缩因子接近1,如0.99,B必然为负(图1方程2是混合物的二次维里系数,式中B ij = B ji为组分i和j的二次交互维里系数,B ii为纯组分i 的二次维里系数。

二次维里系数是温度的函数。

也可用方程3求B,便于手工计算。

比较适合密闭输气计算,方程3中B i的平方根为总因子,参见参考文献1,3,4。

天然气基本压缩因子计算方法

一种新型天然气压缩因子数值计算方法

AGA8—92DC计算方法天然气压缩因子计算(最漂亮的)

天然气压缩因子的分析及其计算

AGA8—92DC计算方法天然气压缩因子计算（最漂亮的）

压缩因子计算

天然气压缩因子计算

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

天然气压缩因子的计算

天然气压缩因子计算及影响因素分析

天然气流量计量中确定压缩因子的方法

用于计量的天然气压缩因子计算方法比较[1]

天然气压缩因子计算及影响因素分析

GB T 17747.1-1999 天然气压缩因子的计算 第1部分：导论和指南

压缩因子计算方法

matlab基于bwrs方程的天然气压缩因子计算

天然气基本压缩因子计算方法

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

GB T 17747.1-1999 天然气压缩因子的计算第1部分：导论和指南