第27卷第10期煤氣與熱力Vol 27 No 102007年10月Gas heatOct液化天然氣相平衡計(jì)算李軍,吳洪松(1.天津城市建設(shè)學(xué)院能源與機(jī)械工程系,天津300384;2.中國(guó)市政工程華北設(shè)計(jì)研究院設(shè)計(jì)四所,天津300074)摘要:闡述了SRK狀態(tài)方程和逸度方程,探討了應(yīng)用SRK方程進(jìn)行相平衡計(jì)算的過(guò)程。結(jié)合計(jì)算實(shí)例,分析了計(jì)算相對(duì)誤差。關(guān)鍵詞:液化天然氣;液化;相平衡中圖分類(lèi)號(hào):TU996文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):1000-4416(2007)10-0018-04Phase Equilibria Calculation of LNGLI Jun, WU Hong-song(1. Dept. of Energy Source and Machine Engineering, Tianjin Institute of Urban ConstructionTianjin 300384, China; 2. Fourth Design Department, North China Municipal EngineeringDesign Research Institute, Tianjin 300074Abstract: The SRK state equation and fugacity equation are described, and the process of phaseequilibria calculation using SRK equations is discussed. The relative calculation error is analyzed with acalculation exampleKey words: LNG; liquefaction; phase equilibrium當(dāng)前,大多數(shù)天然氣液化裝置采用混合制冷劑T——天然氣溫度,K液化流程。混合制冷劑由C1~C5碳?xì)浠衔锛暗猇——天然氣摩爾體積,m/mol氣組成,通過(guò)逐級(jí)冷凝、蒸發(fā)、節(jié)流膨脹,逐步冷卻和a、b天然氣的系數(shù)液化天然氣。混合制冷劑和天然氣的相態(tài)和組成通a、b與天然氣的組成有關(guān)。a的計(jì)算式為過(guò)相平衡來(lái)確定,相平衡計(jì)算是閃蒸計(jì)算的基礎(chǔ),而(2)閃蒸計(jì)算又是各設(shè)備模塊計(jì)算的基礎(chǔ)。本文對(duì)應(yīng)用a=2(a)032(1-k)0.42747R2TSRK方程的相平衡計(jì)算進(jìn)行探討。1+m(1-r01SRK狀態(tài)方程和逸度方程-31m2=0.48+1.574a4-0.1760①SRK狀態(tài)方程對(duì)由n種組分組成的天然氣,其SRK狀態(tài)方程(5)為式中a—組分i的系數(shù)組分i的摩爾分?jǐn)?shù)Pv-6 v(V+6)e一后中國(guó)煤化工數(shù)式中天然氣壓力,P度,KR——摩爾氣體常數(shù),J/(mol·K),取CNMHG8.3145J/(mol·K)m:組分i的常數(shù)www.watergasheat.com李軍,等:液化天然氣相平衡計(jì)算第27卷第10期T—組分i的對(duì)比溫度式(1)~(13)適用于氣相和液相。對(duì)式(13)——組分i的偏心因子在已知y、x的情況下,計(jì)算組分i的氣相逸度系數(shù)二元交互作用系數(shù)k,與所選的狀態(tài)方程有關(guān),d時(shí),取組分i的氣相摩爾分?jǐn)?shù)y;計(jì)算組分i的般查表。部分氣液相平衡體系中狀態(tài)方程的二元液相逸度系數(shù)時(shí),z取組分i的液相摩爾分?jǐn)?shù)交互作用系數(shù)見(jiàn)文獻(xiàn)[2],當(dāng)表中查不到時(shí),可按下x。式計(jì)算:2相平衡計(jì)算過(guò)程Vv2.1相平衡狀態(tài)下的關(guān)系(6)當(dāng)達(dá)到相平衡時(shí),存在以下關(guān)系:式中V、V,組分i的臨界摩爾體積,m3/xd:x中(14)(15)b的計(jì)算式為(16)(7)L=9L(17)0.08664RTe=1-L式中b—組分i的常數(shù)式中K—組分i的相平衡常數(shù)天然氣壓縮因子Z的計(jì)算式為:L—天然氣的液化分率(9)q—液相摩爾流量,molsq—原料氣的摩爾流量,mo式中天然氣壓縮因子e天然氣的氣化分率p天然氣密度,kg/m2.2已知條件和待求量R氣體常數(shù),J/(kg·K)已知條件為天然氣壓力溫度和組成,天然氣組當(dāng)用于氣相時(shí),p取天然氣氣相密度;當(dāng)用于液成以氣態(tài)各組分的摩爾分?jǐn)?shù)表示。待求量為各組分相時(shí),p取天然氣液相密度。的氣相摩爾分?jǐn)?shù)y、液相摩爾分?jǐn)?shù)x1相平衡常數(shù)SRK狀態(tài)方程用壓縮因子Z表示為:K1,i=1,…,n;天然氣氣化分率e和液化分率L。Z-z2+(A-B-B)z-AB=0(10)2.3相平衡計(jì)算過(guò)程(11)2.3.1第1種方法2(RT)①由氣態(tài)各組分的摩爾分?jǐn)?shù)2、天然氣壓力pB=9(12)和溫度T,計(jì)算中、。式中A、B—天然氣的系數(shù)由z、P、T按式(2)~(8)計(jì)算出a、b式(10)適用于氣相和液相。式(10)是一元三b.按式(11)、(12)計(jì)算出A、B。次方程,可采用卡爾丹公式求解解方程(10),對(duì)單相(氣相或液相)天然氣②SRK逸度方程式(10)只有1個(gè)實(shí)根,這個(gè)實(shí)根等于該相的壓縮因SRK逸度方程為:子。對(duì)兩相天然氣,式(10)有3個(gè)實(shí)根,最大的實(shí)b根為氣相的壓縮因子z,最小的實(shí)根為液相的壓縮d=1(Z-1)-ln(2-B)因子Z,中間的實(shí)根無(wú)意義。d將氣相壓縮因子z、x:等代入式(13),計(jì)算2點(diǎn)()2(1-)-1得中國(guó)煤化工代入式(13),計(jì)算(13)得φCNMHG式中d—組分i的逸度系數(shù)按下式設(shè)定各組分相平衡常數(shù)K的初值第27卷煤氣與熱力ww.watergasheat comexp. 42(19)a由y、x1、P、T,按式(2)~(8)分別計(jì)算出氣相和液相的a、b。③求解氣化分率e。b.按式(11)、(12)分別計(jì)算出氣相和液相的在相平衡狀態(tài)下,存在關(guān)系A(chǔ)、Bc.分別解方程(10),得到氣相的壓縮因子Zx=1+e(K1-1)(20)和液相的壓縮因子z。K(21)d將氣相壓縮因子zY、y1等代入式(13),計(jì)算1+e(K:-1)得d。將液相壓縮因子z、x等代入式(13),計(jì)算由式(15)、(16)可得得d②按下式計(jì)算各組分相平衡常數(shù)K將式(21)、(20)代入式(22)得,(25)(K:-1)1+e(K1-1)(23)③由計(jì)算得到的K和液化分率L,計(jì)算新的對(duì)式(23)進(jìn)行迭代求解,可得e,進(jìn)而由式(18)液相摩爾分?jǐn)?shù)x1、新的氣相摩爾分?jǐn)?shù)y求得L。④按式(20)、(21)計(jì)算x1、y,并進(jìn)行歸一化K2+L(1-K)(26)使得:KK;+L(1-K1)(27)④判斷下式是否成立: