BWRS方程在天然气物性计算中的应用 BWRS方程在天然气物性计算中的应用

BWRS方程在天然气物性计算中的应用

  • 期刊名字:油气储运
  • 文件大小:
  • 论文作者:吴玉国,陈保东
  • 作者单位:辽宁石油化工大学油气储运工程系
  • 更新时间:2020-03-24
  • 下载次数:
论文简介

1603年设计计算BWRS方程在天然气物性计算中的应用吴玉国陈保东(辽宁石油化工大学油气储运工程系)吴玉国陈保东:BwRS方程在天然气物性计算中的应用,油气储运,2003,22(10)16~21。摘要在输气管道的水力、热力计算中需要计算许多热物性参数,详细地介绍了应用BwRS方程求解各参数的方法、过程以及物性计算程序编制方法,给出了相应的计算实例。通过分析、实例计算和比敦得出结论,应用BWRS方程计算天然气物性,其计算精度高,适用范围广。主题词天然气物理性质BWRS方程应用在输气管道的工艺计算中,需要使用的热物性由它的临界参数Ta、a及偏心因子a;从下列公式参数有,密度ρ压缩因子Z、焓h、熵s、定压比热c求得:和定容比热cυ、温度绝热指数kr和容积指数k。比B01=A1+B1a1热比k以及焦耳一汤姆逊效应系数D4等。在计算Az+ B2a天然气物性时,必须考虑高压、低温的影响,因此必须采用适用于实际气体的状态方程。A,+B,o一、BWRS方程odri=A+ B2a,b: As +bsaBWRS方程是一个多参数状态方程,其形式As+ bse&at a,+ b, wiP=7+(RF=异+一pacL- As+ bso+(r--号)+(+号p段r=A2+B3a1+(1+2)exp(-m2)(1)=A10+B10a1式中p——系统的压力,kPa;T—系统的温度,K;rTa=A1+B1aexp(-3.8c;)气相或液相的密度,kmol/R—气体常数,R=8.3143式(2)中参数A1~Au和B1~B1的值见表1,这是 Starling K E通过正构烷烃采用多种热力性质kJ/(kmol·K)1、方程式中各参数的求法分析(PVT焓和蒸气压)关联得到的。部分天然气组分的Ta、p和a;值列于表2。由于不同来源的数式(1)中的A,B,C,D,Ea,,a,a,y为状态据通常有一定的出人,在使用本方程时不宜使用其方程式的11个参数。对于纯组分i各参数,均可它来源的数据。113001,辽宁省抚顺市;电话:(0413)6686544第22卷第10期吴玉国等:BWRS方程在天然气物性计算中的应用17表1通用挚数A和B的值2、状态方程的混合法则参数的下标参数的值在计算混合物时,求解参数可按以下的混合规则进行计算。0.44369000.1154491.284380A6=∑∑xx1A0.35630600.54497900.270896C=∑∑xxCC12(1一k)30.52862900.3496210,48401100.754130y=[2x20,07052330.0444480,5040870(3)90,03074520.1794330.0732828a=[∑xa,3]30.00645000.022143[2x ci3]衰2天然气纯物质的物理参数物质名称分子式临界温度T。临界压力Pd=[∑xd,434.604E。=∑2xxEm12E1(1一k)5甲烷CH.乙烷在以上各式中,x、x分别为气相或液相混合物中i和j组分得到的摩尔分数,k为i组分间的异丁烷i-C HIo408.133.648交互作用系数,它表示和理想混合物所发生的偏差,n-C Huk越大,说明偏离越远。对于同一组分,显然k=异戊烷i-CH,460.370。 Starling给出了18种常见组分间的k数据。对正戊烷rCH正己烷3.012于每一个状态方程的k值是互不相同的,在使用正庚烷CHi540.282.736BwRS方程时不应选用其它来源的数据正辛烷CH12.487氮气3.394二、各热力学参数的求解氧化碳CO7.376水蒸气HO(温度T和压强P已知)物质名称临界密度P偏心因子分子量(knol/m)1、计算所需的基础数据甲烷计算所需的数据包括所计算流体的组分,各组0.0130分的摩尔分数和分子量,临界点压强、温度和密度丙烷4.99940.157044,094偏心因子,通用常数A1和B,流体各组分中c、h异丁烷3.80120.1830和s°的常数值,流体中各组分间二元作用系数正丁烷3.92130.197058.1202、BwRS方程的11个参数计算异戊烷3.24690,226072.146首先利用式(2)计算流体中各个纯组分对应的正戊烷0,252072.14B、A0n己烷2.71670.3020、Da、d,和E11个参数86.172正庚烷2.34670.3530然后再利用式(3)的混合规则求出所求流体(混合正辛烷0.4120物)对应的B、A、C0、y、b、a、a、c、D、d和E011个氮气11.09900.0350参数氧化碳10.63800.2100求解流体密度水蒸气18.015BwRS方程的形式可改写成下列函数形式:18·2003年FD-RT+(R-A一导+-裂)为+(r-a-÷y+a(+)p.-()又根据偏导数的循环关系,得到(+》)(=)p=0(4)[(),小]=[)+],并且在给定T、px的情况下求解方程,F(PD)=0,1,pw=RT,所以可得到得出ρ值。用正割法求解密度迭代公式为2:p+,=-(3=a(e2h-)=2-R7x+p-7(B)12(1)式中的下角标k表示迭代序号,再用正割法求将BWRS状态方程代入式(11),即可导出气相解时应先设两个密度初值,定出区间。对于气体初或液相等温焓差的计算公式。值可按理想气体考虚设=0一,选代计算(-)=(BRr2A到|p+1-p|≤c为止6 Eo+1(2T-3a当收敛指标为c=10-‘时,一般P的迭代次数只需要3~6次(注意:此时求出的密度单位为)y2+3a(6a+r)kmol/m3)。同时,由公式p=pZRT,可求得压缩因子Z+[3-(3+2-yp)计算流体的焙h和实际气体的焓不同于理想气体。在给定的压力exp(- yo)(12)和温度下,实际气体的焓一般由在该温度下的理想而纯物质理想气体的焓h可根据其定压比热气体的焓再加上一个修正项求得按式(13)计算。式(5)中的h°为系统温度下理想气体的焓,而?=h%i+c,dT(13)(h-h°)即为修正项,称为等温焓差。同理,可求得式中c—理想气体i的定压比热,实际气体的熵。kJ/kmol·K);5=(s-50)+s°h—理想气体i在参考状态T。、户下的焓式(6)中的。为系统温度下理想气体的熵,而值,kJ/kmol(s-s0)即为修正项,称为等温熵差参考状态及其焓值的选择原则上是任意的,但(1)等温焓(h-h°)差的计算常用绝对温度和绝对压力都等于零的状态作为基由热力学基本理论可知,当单位质量气体从一准。一些手册中列有大量烃类和有关组分的理想气个参考状态T和p到另一个状态T、P时:体焓的数据,为方便使用,将这些数据回归成多(7)项式。h=A,+BT+CT+DT对于等温过程+Er+FiTS(14)式(14)中的A、B1、C、D、E和F是求解cP因BWRS方程以p为显函数,不方便求导数(ah8和s的常数值(一些物质的具体常数值列于表U/aT),为此,可作如下变换,因为d(p)=pdv+2)。对于非烃类气体取T=0,P=0状态下的hcdp,所以0为基准,对于非烃类气体选用该基准时,液体的-)2=(P一)-(0py)()焓常为负值。为避免出现这种情况,AP1手册中烃类组分的焙均以-129℃温度下饱和液体的焙作为由此,恒温下由压力p→0积分至p的焓0基准。第22卷第10期吴玉国等:BWRS方程在天然气物性计算中的应用19对于流体(混合物),其计算常数为:+s8(20)B=(∑xpB,)式中x组分的摩尔分数;式(20)中的s表示纯物质i在参考状态T。=组分的分子量;0,p°=101325Pa时的绝对熵值。同样,s也可用B,—i组分的常数多项式表达为:p—流体(混合物)的平均分子量。s:= BInT+2 Cr+oDT同理可解出其它常数A、C、D、E和F的值。上述方法同样适用于求解理想气体混合物的c。和s+Er +FiT+g(21)的计算。在计算焓值的过程中应特别注意单位的换算,因为在式(14)中,焓的单位是kJ/kg,而在式式(21)中s的单位也为kJ/(kg·K)(12)中,焓差的单位是kJ/kmol。5、气体的比热、绝热指数和节流效应系数(2)等温熵差(s-s9)的计算对于气体而言,通常要计算定压比热c。和定容根据热力学关系式2,3):比热cds- dT+ at).du实际气体混合物在低压下的定压比热为:(22)(16)T+(1)dv也可用下述拟合方程计算低压下的定压比热:cp =B+2 CT+3 DT+4 ET3+5 FT对于等温差,有d=(器)dc,又因为v(23)1/p,所以得到在低压时,可将气体混合物看作理想气体所以:(17)cp-co=R对式(17)进行积分,得到:高压下的定容比热为(-1)-=()(25)Rln/四T代入BWRS方程,得(18)c=2+(-122+2952)代入BWRS方程,得(-)=-R135)[(y2+2)exp(-p2)-2](26)BR+2C-3B+4)这样,由式(24)和式(26)就可以求得c(计算时要注意单位的换算)(-)-号(+共)(一)1 kJ/(kmol K)=lk/(kg.K)(27)号共(-)-3[(+)高压下定压比热与定容比热之差为:cx(-2)-(1+2)p(y)](28)(19)式(19)中的p为101.325/RT,等温熵差的单位是k/(kmol·K)。s°为同一温度及单位压力PR BR+-+华)p=101325Pa)下理想气体混合物的熵。根据式(16),纯物质的s可由式(20)求得:+(+÷)油气储运(1+y2)exp(-y2)(29)子Z的数值(试验值是2.0685)查得氮的临界常数为T=126.15K,RRT+2 BoRT-Ao11.099kmol/m3,a=0.035。按照上述方法计算出状态方程的11个参数,然后代人BWRS方程,用正学)+3(RT一a割法求出密度。经计算得到P=20.6895kmol/m3,+6a(a+号+(+再由公式2=p/RT,求得Z=2.15645,T=273.15K,p=101325Pa(计算时注意单位换算)。3)x(-x)(30)计算的相对误差=215645-2.0685×100%=2.0685以下用到的()(⑦)相同,比热45%的单位均为kJ/( kmol. K)。在低压下作热力计算时,通常所用的气体绝热输入已知的数据T,s指数k=cp/c给定任意的两压强P1,P在高压下求解绝热过程中的状态参数需要使用并按程序计算对应的s,不同状态下的绝热指数,如温度绝热指数kr或容积由点(s1,P1)、(521P2)按线性绝热指数k,而此时的k=c/c,只能称为定压、定关系求出s对应的压强P容比热比,而不能称为绝热指数。按序计算TP对应的s节流效应又称焦耳一汤姆逊效应,温度下降的数值与压力下降数值的比值称为节流效应系数,又P,=P, P2=P,s:=s2,8;=S称焦耳一汤姆逊效应系数(2,即.=(a2)P即为Ts下的压强aP-0(△p(31)输出结果由热力学关系式可以导出节流效应系数D,的结束」计算式:图1已知温度T和熵s求解其它参数的计算程序图算例2。求甲烷含量为0.755(分子分数)的D4=(32)CHCO2混合气体在97.65℃和17.48MPa时的密度。式(32)中c,的单位为kJ/( kmol. K),D,的单已知试验值为m=0.144g/cm3,查文献[4]得位为℃/kPa。到的有关数据见表3。按照上述计算过程绘制计算机程序框图,并编3天然气各组分的已知蕃碱参数制物性计算程序。在工程实际问题中,上述情况只是其中一种,即已知条件可能是另外的两个热物理组分(K) (cm/mol)(MPa)参数。在已知温度T和压强p中的一个及另外的甲烷190.6任意一个热物性参数的情况下,仍然可以应用上述二氧化碳304,294.044.0107.28的物性计算程序。例如在一些工程实际问题中已知温度T和熵s来求解其它参数,这种问题的求解将上述数据代入物性计算程序,求得CH4CO2过程可采用图1的方式,其余情况解法相同。混合气体的密度为m=0.1432g/cm3。计算时要注意单位换算并将vc数据转换成相应的p值。三计算实例计算结果的=0.1432=0.14×100%相对误差0.144算例1。计算0℃及101325Pa时氮的压缩因第22卷第10期吴玉国等:BWRS方程在天然气物性计算中的应用各工况下h、5、Z、c、、k、k、kr及D:热物性见表5。参数的计算如下。天然气组成见表4衰4天然气各组分的摩尔分各工况数据如下。然气组分mo%‖天然气组分ml%(1)p=10×105Pa,t=10℃86.43(2)p=10×10°Pa,t=35℃。CHe(3)p=10×105Pa,t=50℃CsHsccc090.12(4)p=50×10°Pa,t=35℃10.10(5)p=100×105Pa,t=50℃。H2O应用BWRS方程物性计算程序计算,计算结果5各工况下天然气物性参败计算结果况(km/m3)Z(k/(kml·K))(kJ/(kmol.K)(k/(kmol·K))(kJ/(kmol·K)kk(℃/MPa)43610.97411922.338.455729.08911.3221,28781.29785.34(2)0.39800,980612893.8255.6439.280813961.3031.27831.28524.42(3)0,37840.983613487.5257.5239.87810.83891.2931.27231.27773.97(4)2.14870,908312168.4240.5544.493031.00911.4351.31481.31054.145)4.27080.871512020.3156.67~360℃、压力为344.7372~17236.86四结论kPa下的23种混合物的焙差时同试验值的平均绝对偏差为5.117kJ/kg.当预测摩尔分数分别为BWRS方程是在1970年由 Starling K E等人94.8%和5.2%的甲烷丙烷混合物的熵时,72个检在BwR方程基础上提出的一个具有11常数的测点覆盖了温度-156.67~148.89℃,压强状态方程,其目的是拓宽BWR方程的应用范围。1723.686~17389.49kPa的范围,所得到的熵差经改进后,对比温度可低至T=0.3,且在比临界密同文献值的平均绝对偏差为0.0218k/(kg·℃)度高3倍的条件下也能计算气体的PVT关系。在由于状态方程在热力学的一致性,其它容积性质,如计算轻烃气体、CO2、H2S和N2的容积性质时,其误热容焦耳一汤姆逊系数等都具有较高的计算精度。差在0.5%~2%之间,文中使用的式(2)实际上算例2是对于混合物的情况按照式(3)的混合法是 Starling KE和 Han Ms在纯物质的BWRS方则,利用BwRS方程计算,结果表明计算误差很小程的基础上,把方程中11个常数分别和临界参数仅为-0.556%,远远低于RK方程的-1.597%(同A、T构作成无量纲的组合数,并将各组数与偏心为计算例2)。总之,BWRS方程的优点是计算精度因子a进行关联得到的。表1中的数值是由C1~较高,适用范围较广。因此,BWRS方程值得在工C正构烷烃的p、T和a回归得到的。在算例1程计算中推广应用。中,计算BWRS方程的参数时就是应用式(2)计算的。计算结果表明误差为4.25%精度仍然较高高于RK方程的-5.1%及PR方程的-10.37%1, Starling K E, Han M S, Thermo data refined for LPG(Part 1〔同为计算例1)。经过改进的BwRS方程是一个普2,姚光镇输气管道设计与管理石袖大学出版杜(东营),1989.遍化的状态方程,不仅使用方便,而且其使用范围进3,苏长蒋等离等工程热力学,高尊教育出版社(北京),1987步扩大。BWRS方程同样能很好地预测混合物4,薰景山等化工热力学滑华大学出版社(北京),193的性质。 Starling K E和 Han M S用BWRS方程5,童景山;流体的热物理性质,中国石化出版社(北京),198预测温度为-175~237.78℃、压力为101.352834417372kPa下的14种混合物的密度,结果同试(收稿日期:2002-11-13)编辑:刘春阳验值的平均绝对偏差为1.16%。当预测温度为

论文截图
版权:如无特殊注明,文章转载自网络,侵权请联系cnmhg168#163.com删除!文件均为网友上传,仅供研究和学习使用,务必24小时内删除。