2003BWRS*():BWRS,,2003,22(10)16~21在输气管道的水力热力计算中需要计算许多热物性参数,详细地介绍了应用BWRS方程求解各参数的方法过程以及物性计算程序编制方法,给出了相应的计算实例通过分析实例计算和比较得出结论,应用BWRS方程计算天然气物性,其计算精度高,适用范围广天然气物理性质BWRS方程应用,,ZhscpcvkTkv,kDi,,一BWRS方程BWRS,1!:p=RT+B0RT-A0-C0T2+D0T3-E0T42+bRT-a-dT3+a+dT6+c3T21+2exp-2(1)p,kPa;T,K;,kmol/m;R,R=8.3143kJ/(kmol∀K)1(1)A,B,C,D,E,a,b,c,d,,11i,Tcici!i1!:ci∀B0i=A1+B1!iciA0iRTci=A2+B2!iciC0iRT3ci=A3+B3!i2cii=A4+B4!i2cibi=A5+B5!i2ciaiRTci=A6+B6!i3cii=A7+B7!i2ciciRT3ci=A8+B8!iciD0iRT4ci=A9+B9!i2cidiRT2ci=A10+B10!iciE0iRT5ci=A11+B11!iexp(-3.8!i)(2)(2)A1~A11B1~B111,StarlingKE(PVT,)Tcici!i2,*113001,;:(0413)6686544∀16∀1AjBj(j)AjBj10.44369000.11544921.2843800-0.92073130.35630601.70871040.5449790-0.27089650.52862900.34962160.48401100.75413070.0705233-0.04444880.50408701.32245090.03074520.179433100.07328280.463492110.0064500-0.0221432Tci(K)pci(MPa)CH4190.694.604C2H6305.384.880C3H8369.894.250iC4H10408.133.648nC4H10425.183.797iC5H12460.373.374nC5H12469.493.369C6H14507.283.012C7H16540.282.736C8H18568.582.487N2126.153.394CO2304.097.376H2O647.3022.048ci(kmol/m3)!i∀i10.05000.013016.0426.75660.101830.0684.99940.157044.0943.80120.183058.1203.92130.197058.1203.24690.226072.1463.21490.252072.1462.71670.302086.1722.34670.3530100.1982.05680.4120114.22411.09900.035028.01610.63800.210044.01017.85700.344018.0152,1!B0=ixiB0iA0=ijxixjA0i1/2A0j1/2(1-kij)C0=ijxixjC0i1/2C0j1/2(1-kij)3=ixii1/22b=ibixi1/33a=ixiai1/33=ixii1/33c=ixici1/33D0=ijxixjD0i1/2D0j1/2(1-kij)4d=ixidi1/33E0=ijxixjE0i1/2E0j1/2(1-kij)5(3),xixjij,kijij,,kij,,kij=0Starling18kijkij,BWRS二各热力学参数的求解(温度T和压强P已知)1,,,,AiBi,c0pih0is0i,2BWRS11(2)B0iA0iC0iibiaiiciD0idiE0i11,(3)()B0A0C0bacD0dE0113∀17∀2210吴玉国等:BWRSBWRS:F()=RT+B0RT-A0-C0T2+D0T3-E0T42+bRT-a-dT3+a+dT6+c3T21+2exp-2-p=0(4)Tpxi,F()=0,,2!:k+1=k-1F(k)-kF(k-1)F(k)-F(k-1)k,,,1=0,2=pRTk+1-k###=10-4,3~6(:kmol/m3),p=ZRT,Z4hs,h=(h-h0)+h0(5)(5)h0,(h-h0),,s=(s-s0)+s0(6)(6)s0,(s-s0),(1)等温焓(h-h0)差的计算,T0p0Tp:dh=cpdT+u-TvTpdp(7):dhT=v-TvTpTdp(8)BWRSp,(v/T)p,,d(pv)=pdv+cdp,:∃p2p1vdpT=p2v2-p1v1-∃v2v1pdvT(9),p0%0p:(h-h0)T=∃pp0%0v-TvTpdpT(10),∃p2p1vTpdpT=-∃v2v1pTvdvT,v=1p0v0=RT,:h-h0=p-RT+∃0p-TpTd2(11)BWRS(11),h-h0=B0RT-2A0-4C0T2+5D0T3-6E0T4+122bRT-3a-4dT2+156a+7dT5+cT23-3+22-24exp-2(12)h0i(13)h0i=h00i+∃TT0c0pidT(13)c0pi,kJ/(kmol∀K);h00iiT0p0,kJ/kmol,,,,h0i=Ai+BiT+CiT2+DiT3+EiT4+FiT5(14)(14)AiBiCiDiEiFic0pih0is0i(2),T=0,Pa=0h0i=0,,,API-129&∀18∀20030(),:B=xi∀iBi∀(15)xii;∀ii;Bii;∀()ACDEFc0ps0,(14),kJ/kg,(12),kJ/kmol(2)等温熵差s-s0的计算2,3!:ds=cvTdT+pTvdvds=cpTdT+pTvdv(16),ds=pTvdv,v=1/,:ds=-pTd2(17)(17),:s-s0=-∃101.325RTpTd2=-RlnRT101.325+∃0R-pTd2(18)BWRS,:s-s0=-RlnRT101.325-B0R+2C0T3-3D0T4+4E0T5∋-0-12bR+dT22-20+5dT25-50-2cT31+122exp-2-1+122exp-2(19)(19)0101.325/RT,kJ/(kmol∀K)s0(p0=101325Pa)(16),s0i(20):s0i=∃TT0cpdTT+s00i(20)(20)s00iiT0=0,p0=101325Pa,s0i:s0i=BilnT+2CiT+32DiT2+43EiT3+54FiT4+Gi(21)(21)s0ikJ/(kg∀K)5,cpcv:c0p=yiC0pi(22):c0p=B+2CT+3DT2+4ET3+5FT5(23),,:c0p-c0v=R(24):cv=c0v+∃0-T22pT2d(25)BWRS,:cv=c0v+6C0T3-12D0T4+20E0T5+dT22-2ad5T25+3cT32+2exp-∃2-2(26),(24)(26)cv():1kJ/(kmol∀K)=1∀kJ/(kg∀K)(27):cp-cv=T2∀pT2pT(28)pT=RB0R+2C0T3-3D0T4+4E0T52+bR+dT23-dT26∀19∀2210吴玉国等:BWRS-2c6T31+2exp-2(29)pT=RT+2B0RT-A0-C0T2+D0T3-E0T4+3bRT-a-dT2+6a+dT5+3c2T21+2-2324exp-2(30)pTpTcpcvkJ/(kmol∀K),k=cp/cv,kTkv,k=cp/cv,-,,-2,3!,Di=lim%p%0%t%ph=tph(31)Di:Di=1cpT2∀pTpT-1(32)(32)cpkJ/(kmol∀K),Di&/kPa,,,Tp,,,Ts,1,三计算实例10&101325PaZ(2.0685)Tc=126.15K,c=11.099kmol/m3,!=0.03511,BWRS,=20.6895kmol/m3,Z=p/RT,Z=2.15645,T=273.15K,p=101325Pa()=2.15645-2.06852.0685∋100%=4.25%1Ts20.755()CH4CO297.65&17.48MPam=0.144g/cm3,[4]33Tc(K)VC(cm3/mol)MPc(MPa)190.699.016.0434.54304.294.044.0107.28,CH4CO2m=0.1432g/cm3,VCc=0.1432-0.1440.144∋100%=-0.556%∀20∀2003hsZcvcpkkvkTDi4(1)p=10∋105Pa,t=10&(2)p=10∋105Pa,t=35&(3)p=10∋105Pa,t=50&(4)p=50∋105Pa,t=35&(5)p=100∋105Pa,t=50&BWRS,54mol%mol%CH486.43nC60.05C2H61.83nC70.09C3H80.49nC80.02iC40.12N20.62nC40.13CO210.10iC50.06H2O0.01nC50.055(kmol/m3)Zh(kJ/(kmol∀K))s(kJ/(kmol∀K))cp(kJ/(kmol∀K))cv(kJ/(kmol∀K))kkvkTDi(&/MPa)(1)0.43610.974111922.3252.3538.455729.08911.3221.28781.29785.34(2)0.39800.980612893.8255.6439.280830.13961.3031.27831.28524.42(3)0.37840.983613487.5257.5239.878430.83891.2931.27231.27773.97(4)2.14870.908312168.4240.5544.493031.00911.4351.31481.31054.14(5)4.27080.871512020.3234.9949.206532.23581.5261.41031.30632.82四结论BWRS1970StarlingKEBWR11,BWR,Tr=0.3,3PVTCO2H2SN2,0.5%~2%5!(2)StarlingKEHanMSBWRS,11cTc,!1C1~C8cTc!1,BWRS(2),4.25%,,RK-5.1%PR-10.37%(1)BWRS,,BWRSStarlingKEHanMSBWRS-175~237.78&101.3528~34473.72kPa14,1.16%-156.67~360&344.7372~17236.86kPa23,5.117kJ/kg94.8%5.2%,72-156.67~148.89&1723.686~17389.49kPa,0.0218kJ/(kg∀&),,-2,(3),BWRS,,-0.556%,RK-1.597%(2),BWRS,,BWRS1,StarlingKE,HanMS:ThermodatarefinedforLPG(Part14:Mixture),HydrocarbonProcessing,197251(5),129.2,:,(),19893,:,(),19874,:,(),19935,:,(),1996(:20021113):刘春阳∀21∀2210吴玉国等:BWRS,1965,1987,1966,1988(),1991(),(),1970,2001(),,()1978,2001,,1977,2000,,1974,2002(),,(),1963,1986,,,1960,1993,,1968,1992,,1972,1995,,,1964,1983,1987,,,,1966,1989,1979,2001,,1970,1992,∀60∀2003emphasesofslugflowinthefuture.Theworkaimsatful