低溫甲醇脫碳工藝熱力學(xué)模型和過(guò)程模擬研究

第39卷第7期化學(xué)工程Vol 39 No. 72011年7月CHEMICAL ENGINEERING( CHINA)低溫甲醇脫碳工藝熱力學(xué)模型和過(guò)程模擬研究肖珍平12,房鼎業(yè),應衛勇,劉森,陳健(1.華東理工大學(xué)大型工業(yè)反應器T程教育部工程研究中心化學(xué)工程聯(lián)合國家重點(diǎn)實(shí)驗室,上海200237;2.中國石化集團寧波工程有限公司,浙江寧波315103;3.清華大學(xué)化學(xué)工程聯(lián)合國家重點(diǎn)實(shí)驗室,北京10004)摘要:為了進(jìn)行含高壓氣體和強極性物質(zhì)的體系的低溫甲醇脫碳工藝的模擬研究建立相應的熱力學(xué)模型。采用Soave-Redlich- Kwong(SRK)立方型狀態(tài)方程,結合 Huron-Vidal混合規則和非隨機雙流體Non- Random-Two- Liqui(NRTL)活度系數模型建立的熱力學(xué)模型,從氣體溶解度和氣液平衡數據擬合獲得了45對活度系數模型參數。用于低溫甲醇洗脫碳工藝的過(guò)程模擬,低壓和高壓脫碳系統的模擬結果和實(shí)際T業(yè)數據符合很好。結果表明:采用熱力學(xué)模型適合于同時(shí)含有H2,N2,CH4等非極性氣體和CH2OH,H2O等強極性溶劑的低溫甲醇脫碳工藝的過(guò)程模擬關(guān)鍵詞:CO2捕集;氣體吸收;低溫甲醇;過(guò)程模擬中圖分類(lèi)號:TQ021.4文獻標識碼:A文章編號:10059954(2011)07-007404Thermodynamic model and process simulation for CO2removal with low temperature methanolXIAO Zhen-ping", FANG Ding-ye, YING Wei-yong, LIU Sen', ChEN Jian(1. Engineering Research Center of Large Scale Reactor Engineering and Technology of the Ministry of EducationState Key Laboratory of Chemical Engineering, East China University of Science and Technology, Shanghai200237, China; 2. Sinopec Ningbo Engineering Company Limited, Ningbo 315103, Zhejiang Province, China3. State Key Laboratory of Chemical Engineering, Tsinghua University, Beijing 100084, China)Abstract: For the study on simulation of CO2 capture processes with low temperature methanol, a thermodynamicmodel was established for a system with gases and strong polar components under high pressure. The SoaveRedlich-Kwong(SRK cubic equation of state, Huron- Vidal mixing rule and Non-Random-Two-Liquid (NRTLmodel for activity coefficients were used, and 45 pairs of parameters activity coefficient model were regressed withas solubility and vapor-liquid equilibrium data. The simulation of the CO2 removal pros was camed outresults show that the thermodynamic method is valid for the process simulation with nonpolar gases such as H2, N2CH4 and strong polar components such as CH3 OH and H,0.Key words: CO2 capture; gas absorption; low temperature methanol; process simulation氣候變暖已成為全世界關(guān)注的焦點(diǎn)問(wèn)題,為等)和物理吸收溶劑(低溫甲醇、聚乙二醇二甲醚、了減少CO2的排放除了節能降耗、使用新能源保碳酸丙烯酯和N甲基吡咯烷酮等)。采用低溫甲醇護環(huán)境以外,CO2的捕集和封存技術(shù)(CCS)顯得尤吸收分離脫除CO2和H2S的過(guò)程,通常稱(chēng)為低溫甲為重要。但目前CO2的捕集與封存技術(shù)的實(shí)施費醇洗工藝。目前常用于合成氨和煤化工等工業(yè)中的用和能耗仍然很高降低其能耗和成本是目前國際高壓合成氣的脫硫和脫碳。它具有吸收能力強、溶研究熱點(diǎn)(2)。劑便宜不降解的優(yōu)點(diǎn)。但是流程比較復雜,特別是吸收法有較悠久的歷史,大量用于煉油、合成由于低溫和分別脫硫脫碳的要求,物流數量較多,換氨制氫、天然氣凈化等工業(yè)過(guò)程,相對比較成熟,有熱系統比較復雜。大量經(jīng)驗可以參考。目前使用的主要溶劑有化學(xué)吸流程比較復雜多變,就需要有完整精確的流程收溶劑(有機胺水溶液和無(wú)機堿碳酸鉀、氨水溶液模擬程序。但由于體系復雜溫度和壓力范圍變化收稿日期:20110127作者簡(jiǎn)介:肖珍平(1965-)男博士研究生,主要從事化工工藝系統研究和設計工作, E-mail: xiaochu,snoc@ sinopec,com肖珍平等低溫甲醇脫碳工藝熱力學(xué)模型和過(guò)程模擬研究75大流程模擬中熱力學(xué)模型的選擇就稱(chēng)為關(guān)鍵。式中:G為分子相互作用能,G=xp[-ayT本文從體系的特點(diǎn)出發(fā),選擇了狀態(tài)方程加吉布斯(RT)];下標jk分別代表各組分;7,為溫度的函混合規則的熱力學(xué)方法獲取熱力學(xué)方程中的參數,數;c,d為溫度參數。對于每一對組分i和j,活度系從而進(jìn)行該體系的流程模擬。數的計算需要3個(gè)參數,即相互作用能參數r;,7以及非隨機因子a(a4=a)。a與溫度及溶液組成1氣體吸收溶解度計算的熱力學(xué)模型無(wú)關(guān),本文取0.25。在用低溫甲醇溶液對CO2進(jìn)行吸收分離的過(guò)程中體系中的氣體組分包括CO2,CO,H2,CH4,Ar,2模型參數的計算N2,H2S,COS,再加CHOH和H2O,同時(shí)整個(gè)體系中根據已有的汽液平衡資料,回歸了所有組分之的各個(gè)塔的壓力范圍大,從常壓到80MPa左右。對間的相互作用參數。如果某個(gè)二元體系的數據,溫于這樣一個(gè)組分沸點(diǎn)差距較大壓力變化較大,又同度范圍較寬的比較重要的2個(gè)組分之間,可以采用時(shí)含有非極性和極性締合組分的體系,采用狀態(tài)方3個(gè)溫度系數的參數,例如CO2CH2OH,是本項目程計算是最合適的。本文選擇 Soave-Redlich-kwog中最為重要的參數。H2OCH2OH參數主要是在溶(SRK)方程進(jìn)行計算劑回收塔中發(fā)揮作用,采用的是2個(gè)溫度系數的參-b v(V+(1)數。有一些二元體系的數據溫度范圍很窄主要是惰性氣體之間的,對模擬計算影響不大,本文只采用式中P為壓力,R為氣體參數,T為溫度,V為摩爾1個(gè)系數的參數。體積。純物質(zhì)參數a和b直接采用PRO/Ⅱ中給定本文回歸了45對NRTL模型的相互作用參數,的參數。有些二元體系沒(méi)有數據,采用替代組分的數據,其中在采用狀態(tài)方程法計算相平衡時(shí),如組分極性主要是COs,用CO2的代替。和體積大小差別很大,采用常用線(xiàn)性混合規則就不在不同溫度下,二元體系CO2CH3OH的計算適合。本文涉及的甲醇和水都是強極性組分,和結果和文獻[810]中數據的比較見(jiàn)圖1。CO2、H2等輕組分,所以必須采用更為精確的超額自140由能CE型混合規則。在PO/Ⅱ中只有 Huron423.15v373.15Vidal混合規則°是考慮了組分間的強烈的非理想105°323.15x298.15性的混合規則,是最合適的。同時(shí)它結合了狀態(tài)方237.15▲213.15程和活度系數的優(yōu)點(diǎn),可同時(shí)計算汽液二相的性質(zhì):6.RT ART(2)∑x(3)式中:x為摩爾分數;下標i為組分;a4,b;分別為純組分i的SRK方程參數。在 Huron-Vidal混合規則0.4中,需要有一個(gè)活度系數模型計算超額自由能G。r(CO,)圖1CO2-CHOH二元體系的汽液平衡計算在PRO/Ⅱ中,可以使用 Non-Random-Two-LiquidFig. 1 Vapor-liquid equilibrium calculation for binary(NRTL)方程。NRTL方程是由 Renon和system CO,-CH,OHPrausnitz1從局部濃度概念出發(fā)獲得的一個(gè)非隨機雙流體模型活度系數模型,特別適用于含有強極性二元體系H2SCH1OH的結果見(jiàn)圖2。從組分的體系的活度系數y,表達式如下:圖1和圖2可以看出酸性氣體CO2和H2S和主要G∑x4G溶劑甲醇的計算結果與文獻數據完全吻合,為低溫甲醇洗體系的過(guò)程模擬打下了堅實(shí)基礎。2《4啊叫2氣體和水的汽液平衡計算,用二元體系CO2(4)H2O作為例子,計算結果見(jiàn)圖3,該體系的計算結果表明1,氣體在水中的溶解度也可以得到很好的(5)模擬和計算?；瘜W(xué)工程2011年第39卷第7期二元體系CH2OH-H2O也是體系中的一個(gè)關(guān)V273.15鍵,對于溶劑回收塔具有重要意義。圖5是該二元°258.15248.15體系的計算結果,在整個(gè)溫度范圍243-423K內▲223.15計算結果和文獻[15]數據非常吻合,為溶劑甲醇的純化精度及其能耗估算提高了計算精度。10.0001.000x0x(Hrs)-日a圖2H2SCH3OH二元體系的汽液平衡計算0010Fig- 2 Vapor-liquid equilibrium calculation for binary system H,S-CH,OH040.60.81.0圖5CH,OH2O二元體系的汽液平衡計算Fig 5 Vapor-liquid equilibrium calculation for binary system CH, OH-H,0T/K對于氣體和液體的焓值和熵值的計算,本文直▲373.15接采用PRO/Ⅱ中的計算程序,即采用SRK方程和·348.15v32HV混合規則及其從汽液平衡數據擬合得到的參數308.15x298.15進(jìn)行熱量的計算。后面的計算結果表明,總體上和288.15實(shí)際情況符合。液體密度數據采用API數據庫,氣0.010.03體密度也采用SRKV方程進(jìn)行計算圖3CO2-H2O二元體系的汽液平衡計算3模擬計算結果Fig-3 Vapor-liquid equilibrium calculation for binary system CO-H,o3.1低壓系統的模擬結果惰性氣體和甲醇的汽液平衡計算,用二元體系低壓系統的數據來(lái)自于國內某廠(chǎng),表1是其脫N2CH2OH作為例子計算結果見(jiàn)圖4,由于文獻[14]硫塔的計算結果結果和現場(chǎng)數據符合很好,H2S摩數據比較分散所以直接釆用壓力的計算值和文獻值爾分數計算值為14×10°。脫碳塔的凈化氣的計進(jìn)行比較。該體系的計算結果表明惰性氣體在甲醇算結果見(jiàn)表2,結果和現場(chǎng)數據也非常符合,CO2摩中的溶解度也可以得到很好的模擬和計算。爾分數計算值為11×10-°,現場(chǎng)數據為10×10。囊1合成氣經(jīng)脫硫后的模擬計算結果(2.9MPa)2. 9 MPa after h s removal脫硫后的合成氣現場(chǎng)值本文計算值流量/(kmol·h1)9859.79868.7溫度/℃-23.2壓力/MPa2.92.9H,0.570156970.004264.238E-030.38610.3865H, s1.438E-06組分摩爾COS01.78lE-07100分數9.0E-058.188E-050.038470.0383A圖4N2CH2OH二元體系的汽液平衡計算7.6E-047.555E-04CH、OH2.0E-044.033E-04Fig4 Vapar-liquid equilibrium calculation for binary system N2-CH,OHH0肖珍平等低溫甲醇脫碳工藝熱力學(xué)模型和過(guò)程模擬研究77囊2合成氣經(jīng)脫硫磯后的模擬計算結果(28MPa)低壓28MPa和高壓77MPa條件下的低溫甲醇脫Table 2 Simulation results for synthesis gas at硫脫碳的流程進(jìn)行了模擬,模擬結果和現場(chǎng)結果也2.8 MPa after H,S and CO, removal符合很好。模擬結果表明,本項目選擇的熱力學(xué)模脫碳后的合成氣現場(chǎng)值本文計算值型及回歸得到的參數是正確的,可以用于各種壓力流量/( kmol h)5988.35989.8條件下的低溫甲醇洗脫碳系統的模擬和設計。溫度/℃45.8-45.6壓力/MPa282.8參考文獻H20.92910.9295[1] IPCC. Climate change 2007-the physical science basis6.8lE-036.770Esummary for policymakers of the working group I repo1.0E-051.148E[R]. Cambridge: Cambridge University Press, 2007H,s00[2]費維揚,艾寧陳健溫室氣體CO2的捕集和分離組分摩爾0分離技術(shù)面臨的挑戰與機遇[冂].化工進(jìn)展,2005,24分數1.3E-041.382E-04(1):18[3]陳健,于燕梅唐宏青含二甲醚二元體系相平衡計算N20.062670623及對分離流程的影響[冂.天然氣化工,2005,30(1):1.22E-031.221E-0371-78.CH, OH 8.0E-05 7. 116E-05 [4] SOAVE G. Equilibrium constants from a modified Redlich-H,Oquation of state [J]Eng Sci, 1972, 27(6):1197-203[5]CHEN Jian, LI Zongcheng. Study on the GE-type mixing3.2高壓系統的模擬結果rule for cubic equations of state[ J]. Tsinghua Science高壓系統數據來(lái)自國內某廠(chǎng)的現場(chǎng)數據,吸收and Technology, 1996, 1(4): 410-415塔壓力是7.7MPa。合成氣經(jīng)過(guò)脫硫和脫碳后的模[6] HURON M J, VIDAL J. New mixing rules in simpl擬計算結果見(jiàn)表3,模擬結果和現場(chǎng)數據也非常符equations of state forVle ofideal mixtures [J]. Fluid Phase Equilibria, 1979, 3(4)合,H2S和CO2的摩爾分數計算值均小于1×10-6。255-271[7] RENON H, PRAUSNITZ J M. Local compositions in豪3脫硫碳后的合成氣的模擬計算結果(7.7MPa)thermodynamic excess functions for liquid mixtures [J]Table 3 Simulation results for purified synthesis gas atAIChE J,1968,14(1):135-1447.7 MPa after H,s and CO2 removal[8]KATAYAMA Takashi, OHGAKI Kazunari, MAEKAWAGoro, et al. Isothermal vapor-liquid equilibria of acetone-脫硫碳后的合成氣現場(chǎng)值計算值carbon dioxide and methanol-carbon dioxide systems at流量/(kmol·h-1)3935.8high pressure[ J]. J Chem Eng Jap, 1975, 8(2):89-92.溫度/℃-56.39[9]HONG Jane H, KOBAYASHI Riki. Vapor-liquid equilib-壓力/MPa7.77.7rium studies for the carbon dioxide- methanol systemH20.95270.9527[J]. Fluid Phase Equilibria, 1988, 41(3): 269-276.CO0.03730.0371[10] FISCHER K, CHEN Jian, PETRI M, et al. Solubility ofHS and CO, in N-octyl-2-pyrrolidone and of H,S in metha-3.940E-07nol and benzene[J]. AIChE J, 2002, 48(4): 887-893H6.709E-23[11] RUPERT D, PATRICK JR, TEJA A S. High pressure組分摩爾cOSphase equilibria in the carbon dioxide n-hexadecane and分數2.38E-032.392E-03carbon dioxide-water systems[J]. Can J Chem Eng2.24E-032.266E-031988,66(2):319-3235.34E-035.503E-03[12] TOCHIGI Katsumi, KOJIMA Kazuo.Prediction oflar gas solubilities in water, alcohols and aqueousCH3OH2.266E-052507E-05solutions by the modified ASOG method[ J]. Fluid PhaseH,O6.967E-09Equilibia,1982,8(3):221-232[13]HENNI A, MATHER A E. Solubility of CO2, N,0, CHand C,H. in polar solvents [J]. J Can Petrol Tech4結論1999,38(13):4648根據低溫甲醇洗脫碳系統的特點(diǎn),選擇了SRK[14] KATAY AMA Takashi, NITTA Tomoshige. Solubilities of方程,結合采用NRTL模型的HⅤ混合規則,建立了hydrogen and nitrogen in alcohols and幾 hexane[]同時(shí)用于汽液二相的熱力學(xué)模型。根據汽液平衡數Chem Eng Data,1976,21(2):194196.據擬合回歸了45對二元相互作用參數。所有參數[15 HIRATA Makoto, OHE Shuzo, NAGAHAMA KunioComputer aided data book of va在低溫甲醇洗脫碳系統工藝的范圍內適用。分別對[M]. Tokyo: Kodansha, 1975.