第12卷第6期燃燒科學與技術Vol 12 No 6006年12月Journal of Combustion Science and TechnologyCaO與CO2循環(huán)反應動力學特性李振山,蔡寧生,趙旭東,黃煜煜(清華大學熱能工程系熱科學與動力工程教育部重點實驗室,北京10004)摘要:為分析CaO與CO2循環(huán)反應過程動力學參數(shù)的變化規(guī)律,利用縮核模型對CaO與CO2的循環(huán)反應特性進行了研究結(jié)果表明:CaO與CO2的化學反應速率常數(shù)在590-743℃為一常數(shù);化學反應速率常數(shù)和產(chǎn)物層擴散系數(shù)隨循環(huán)次數(shù)的增加而降低,溫度對二者的降低規(guī)律有一定的影響;由產(chǎn)物層擴散階段所引起的活性降低主要發(fā)生在前幾次循環(huán)內(nèi),化學反應動力學控制階段的活性降低是引起CaO整體轉(zhuǎn)化率降低的主要因素關鍵詞:動力學參數(shù);產(chǎn)物層擴散;循環(huán)反應中圖分類號:TQ028.17文獻標志碼:A文章編號:1006-8740(2006)06-0481-05Kinetic Characteristic of the Multiple CarbonationCycle Reactions of Cao with CO2LI Zhen-shan, CAI Ning-sheng, ZHAO Xu-dong, HUANG Yu-yuKey Laboratory for Thermal Science and Power Engineering, Ministry of Education of ChinaDepartment of Thermal Engineering, Tsinghua University, Beijing 100084, China)Abstract: In order to analyze the changing trend of kinetic parameters of carbonation/ calcination cycle of CaO, the cyclicvestigated using the shrink core model. It is found that the chemical reaction rate constantof the carbonation reaction is nearly a constant at 590--743 C; chemical reaction rate constant and product layer diffusioncoefficient decrease with the increase of carbonation/ calcination cycle numbers and are both influenced considerably by tem-perature; the reactivity loss in product layer diffusion stage only occurs at the first several cycles and the reactivity loss inchemical reaction controlling stage is the main factor leading to the entire decrease of Cao cyclic carbonation conversionKeywords: kinetic parameters; product layer diffusion; carbonation/ calcination cycleCaO與CO2的碳酸鹽化反應可以應用在很多方LANL提出的零排放的煤制氫/發(fā)電技術和日本新能面,如利用CaO與CO2的碳酸鹽反應及其逆反應(即源綜合開發(fā)機構(gòu)提出的稱為 Hy Pr-RING的煤利用系CaCO3的分解反應)來吸收CO2,增加碳與蒸汽氣化反統(tǒng)4,另外,CaO與CO2的碳酸鹽化反應及其碳酸鹽應過程,而利用CaO循環(huán)吸收CO2以減少溫室氣體排的煅燒反應,也被應用于能源存儲系統(tǒng)和化學熱泵系放也越來越受到人們的重視同時,在天然氣蒸汽重整統(tǒng)制氫方面,利用CaO循環(huán)吸收CO2,以使重整反應、水在應用CaO來吸收CO2過程中,其反應主要為混氣置換反應和CO2的分離過程在一個反應器內(nèi)進行,合氣中的CO2和顆粒中CaO之間的氣固反應,而且也受到了很多研究者的關注13.此外,在各種各樣的CaO在過程中一般都是循環(huán)使用的,即CaO先與CO2零排放煤制氫系統(tǒng)方案中,以CO2接受體氣化法(即反應生成CaCO3,反應一段時間后,對生成的CaCO利用CaO來循環(huán)吸收CO2)為基礎的無氧氣化煤制氫加中國煤化工利用生成Ca0來吸收零排放利用系統(tǒng)得到了許多研究機構(gòu)的重視如美國Co2CNMHG,Ca0與CO2的反應收稿日期:200603-02基金項目:國家高技術研究發(fā)展計劃(863)資助項目(2003AA501330)雋者蘑山(93).男,博土研究生,① mails. tsinghua. edu.en482燃嬈科學與技術第12卷第6期動力學特性數(shù)據(jù)對于上述幾種應用過程而言是十分重1-(1-X)5=k·(ca2-c)t要的,動力學參數(shù)隨循環(huán)反應/煅燒次數(shù)的變化關系對3.566×10exp(-19680/7于相關過程的分析和后續(xù)的反應器設計都十分關鍵RT本文在前期CaO與CO2循環(huán)反應特性實驗研究6的當整個反應過程由產(chǎn)物層內(nèi)擴散速度控制時,反基礎上,主要考察循環(huán)次數(shù)對反應過程動力學參數(shù)的應方程為影響,為后續(xù)的研究和應用提供參考z-[z+(1-z)(1-x)實驗裝置和測量系統(tǒng)Dofr'(Cco, -.)(4)式中k為化學反應速率常數(shù),m3·s1·moll;co為1.1實驗系統(tǒng)介紹CO2反應氣體濃度,mol·m-3;c。為CO2平衡濃度,使用熱重分析儀 TA Instrument1200進行CaO與mol·m-3;t為反應時間,s;z為產(chǎn)物CaCO3與反應物O2反應實驗,TGA的升溫速率為50℃/min,當溫度達到設定溫度并穩(wěn)定后,開始通入含有CO2的氣體CaO的摩爾體積比,2.8;Dm為產(chǎn)物層擴散系數(shù),ms·m進行反應,恒溫反應時間為30min,然后停止通入含CO2的混合氣,通入N2,TCA以50℃/min的升溫速圖1為TGA實驗給出的樣品反應20min后的轉(zhuǎn)率升溫,待溫度穩(wěn)定在設定的煅燒溫度水平時,恒溫化率曲線.從圖中可以看出,CaO的碳酸鹽化過程首10min,然后TGA以50℃/mn的降溫速率降溫,待溫是反應速率非?？斓幕瘜W反應控制階段,然后反應轉(zhuǎn)化率的增加速率明顯降低,最后變化趨勢趨于平緩度穩(wěn)定在設定的反應溫度水平時,開始通人含有CN物層擴散階段8.因此,如圖1所示,對化學的氣體進行反應,至此,完成一個循環(huán)樣品采用分析純CaO,粒徑為20~60Hm,CaO與CO2反應時間約30反應控制階段使用式(3)進行描述,而使用式(4)對產(chǎn)min,CaCO3煅燒溫度870℃,時間10min.CO2體積分物層擴散控制過程進行計算數(shù)分別為14%25%和35%1.2實驗數(shù)據(jù)的整理假定TGA原始曲線起始增重點處樣品質(zhì)量為m0,原始樣品的純度為P,假定雜質(zhì)不發(fā)生反應,且反應過程中雜質(zhì)的質(zhì)量不發(fā)生變化,反應進行到t時刻,樣品的質(zhì)量為m2,則固體反應物轉(zhuǎn)化率用式(3)(1)反應時間/min對應于本文CaO吸收CO2實驗,P=0.98,CaCO圖1caO轉(zhuǎn)化率隨時間變化曲線相對分子量為Mco2=100,Ca0相對分子量為Mc0=56,所以有2實驗結(jié)果與分析X=1.299(2)21CaO與CO2反應的動力學特性假設CaO試樣是由具有一定粒度的微粒組成,并CaO在不同溫度下與CO2反應的轉(zhuǎn)化率曲線可且CaO微粒具有規(guī)則的幾何形狀和各向同性的反應參見文獻[6].值得指出的是,CaO的轉(zhuǎn)化率并不是隨活性在反應過程中,CO2氣體在CaO微粒間的空隙向溫度的增加而一直增加的溫度增加時,cn增加,使得固態(tài)試樣內(nèi)部擴散的同時,也與其周圍的微粒反應由(cc0-c)的差值減小,當溫度增加到一定值時于與固態(tài)試樣相比,CaO微粒本身比較致密,因此收縮k(eo中國煤化工0的轉(zhuǎn)化率增加變得核模型可以用來直接描述每一個微粒的動力學行緩慢為.CNMHG,小于或等于cn此時a便不冉習UU2灰生反應;與此相反,此時當整個反應過程由化學反應速度控制時氣體反CaCO的分解反應開始進行從圖2中可以看出,在應物濃度在顆粒內(nèi)均勻分布,反應方程為CO2體積分數(shù)為14%條件下,當溫度從690℃開始增加時,CaO轉(zhuǎn)化率增加速度開始下降,當溫度達到206年12月李振山等:CaO與CO2循環(huán)反應動力學特性743℃時,CaO轉(zhuǎn)化率增加得非常緩慢,而繼續(xù)升高溫度至大約755℃時,CaO已基本不與CO2發(fā)生反應720℃—690℃640℃590℃由此可以看出,CaO+CO2=CaCO3是放熱反應,溫度過高顯然不利于反應的進行圖3不同溫度下CaO最終轉(zhuǎn)化率隨循環(huán)次數(shù)的變化反應溫度/℃2CaO轉(zhuǎn)化率隨溫度變化的曲線利用式(3)對各溫度下化學反應動力學控制階段的轉(zhuǎn)化率曲線(其中的近似直線段)進行計算,求出的化學反應速率常數(shù)見表1.從表1中可見,Ca0與CO2反應過程的反應速率常數(shù)可以認為在590~743℃范圖4不同循環(huán)次數(shù)下的CaO比表面積圍內(nèi)是不變的,這一點與文獻[8]所得的結(jié)果是一致的60min,鍛燒溫度950℃,鍛燒時間20min,CaO粒徑20~60μm,CO2體積分數(shù)為35%時,不同循環(huán)次數(shù)下表1不同溫度下CaO與CO2化學反應速率常數(shù)的Ca0比表面積結(jié)果,從中可以看出,隨循環(huán)反應次化學反應常數(shù)溫度/中(CO2)顆粒粒徑數(shù)的增加,CaO比表面積快速減少,因此會導致CaO轉(zhuǎn)化率隨循環(huán)次數(shù)的增加而快速降低2.2492.6392.164%445420~60雖然隨循環(huán)次數(shù)的增加,CaO最終轉(zhuǎn)化率都會降20-60低,但在不同的溫度下,CaO最終轉(zhuǎn)化率的降低趨勢有2.2820-60所不同.低溫下caO最終轉(zhuǎn)化率隨循環(huán)次數(shù)的增加而2.225快速降低,最后趨于一穩(wěn)定值后不再變化.較高溫度下2.17820-60CaO最終轉(zhuǎn)化率的降低比較緩慢,見圖3.2.18523CaO與CO2循環(huán)反應動力學特性因為在每次循環(huán)的煅燒過程中,CaO顆粒要發(fā)生22CaO循環(huán)吸收CO2特性燒結(jié),致使CaO顆粒微觀孔隙結(jié)構(gòu)發(fā)生變化和比表面圖3是不同溫度下Ca0循環(huán)吸收Co2的能力變積減少,而二者的這些變化會使化學反應速率常數(shù)和化情況從圖中可以明顯看出,隨著CaO碳酸鹽化/煅產(chǎn)物層擴散系數(shù)降低,圖5和圖6給出了化學反應速燒循環(huán)次數(shù)的增加,CaO吸收CO2的能力顯著降低率常數(shù)和產(chǎn)物層擴散系數(shù)隨循環(huán)次數(shù)的變化關系反應溫度不同,變化趨勢有所不同.從圖3也可以明顯從圖5中可見,反應速率常數(shù)隨循環(huán)次數(shù)的增加看出,CaO循環(huán)轉(zhuǎn)化率隨溫度的升高而增加,主要是因而降低,當溫度在640~720℃范圍內(nèi)時,反應速率常為不同溫度下Ca0與CO2反應的臨界產(chǎn)物層厚不同數(shù)降低趨勢基本相同;而當溫度為590℃時,反應速率而致常中國煤化工4次循環(huán)內(nèi)然后隨CNMHG物層擴散系數(shù)的降低顆粒此時會發(fā)生燒結(jié)現(xiàn)象6燒結(jié)影響CaO吸收緩慢下降,劑的表面積而比表面積對CaO吸收CO2有著非常重圖6所最后保持在一定值且基本不再變化,如要的影響0.圖4為反應溫度700℃,反應時間484燃燒科學與技術第12卷第6期0620.1循環(huán)次數(shù)循環(huán)次數(shù)圖5不同溫度下化學反應速率常數(shù)隨循環(huán)次數(shù)的變化圖8各反應階段轉(zhuǎn)化率隨循環(huán)次數(shù)的變化(590℃)制階段的CaO吸收CO2起作用在產(chǎn)物層控制階段720℃CO2通過CaO顆粒表面的CaCO3產(chǎn)物層向內(nèi)部擴散CaO顆粒比表面積變化對CO2的產(chǎn)物層擴散影響主590℃要發(fā)生在前幾次循環(huán)過程中,對其后循環(huán)過程的擴散日.系數(shù)影響很小,而且即使在前幾次循環(huán)過程中這種影響也會隨溫度的降低而減弱,所以循環(huán)次數(shù)對產(chǎn)物層擴散控制階段的CaO轉(zhuǎn)化率的影響較小圖6不同溫度下產(chǎn)物層擴散系數(shù)隨循環(huán)次數(shù)的變化3結(jié)論(1)在CO2含量一定時,CaO與CO2反應轉(zhuǎn)化率為了進一步分析反應不同階段對CaO反應活性隨溫度的增加先增加而后降低,存在最佳反應溫度降低的影響給出了反應不同階段CaO轉(zhuǎn)化率隨循環(huán)次數(shù)的變化情況,見圖7和圖8.圖中X1代表化學反區(qū)間(2)CaO與CO2的化學反應速率常數(shù)在590應控制階段氧化鈣的轉(zhuǎn)化率;x2為整個反應階段的最743℃為一定值不隨反應溫度的變化而改變終轉(zhuǎn)化率;△X1為從化學反應控制階段結(jié)束后到產(chǎn)物(3)燒結(jié)使得CaO顆粒比表面積隨循環(huán)次數(shù)的層控制階段開始時的轉(zhuǎn)化率△x2為從產(chǎn)物層控制階增加而減小,導致CaO與CO2的化學反應速率常數(shù)和段開始時到最終的轉(zhuǎn)化率.從圖7和圖8可以看出產(chǎn)物層擴散系數(shù)隨循環(huán)反應次數(shù)的增加而降低,致使當溫度在590℃和720℃時,Ca0的活性降低主要由Ca0吸收CO2的能力隨著CaO碳酸鹽化/煅燒循環(huán)次X1的下降所引起數(shù)的增加而顯著降低從上述的分析可知,CaO活性降低主要發(fā)生在化(4)CaO活性隨循環(huán)次數(shù)的增加而降低主要發(fā)生學反應控制階段;在產(chǎn)物層擴散控制階段循環(huán)次數(shù)對在化學反應控制階段;在產(chǎn)物層擴散控制階段循環(huán)次CaO轉(zhuǎn)化率的影響較小這主要是因為在每次循環(huán)的數(shù)對Ca0轉(zhuǎn)化率的影響較小煅燒過程中,CaO顆粒要發(fā)生燒結(jié),致使CaO顆粒比表面積減少(見圖4),而這些變化主要對化學反應控參考文獻[1 Han C, Harison DP. Simultaneous shift and carbon dioxideseparation for the direct production of hydrogen[J].ChemEng Sci,1994,49(24):5875-588△[2] Balasubramanian B, Lopez A, Kaytakoglu S, et al. Hydro.000methane in a single-step process[J]. Chem Eng中國煤化工55CNMH Ghio Yoshizawa. Study onregenerative fuel reformer for a zero-emission vehicle system[J]. J Chem Eng Japan,2003,36(7):860-866[4] Lin Shiying, Yoshizo Suzuki, Hiroyuki Hatano, et al. De-圖7各反應階段轉(zhuǎn)化率隨循環(huán)次數(shù)的變化(720℃)veloping an innovative method, HyPr-RING, to produce2006年12月李振山等:CaO與CO2循環(huán)反應動力學特性hydro-gen from hydrocarbons[ J ]. Energy Conversion andery,1976,59(9/10):379383Managemen,2002,43(9/12):1283-1290[10] Borgwardt R H. Calcium oxide sintering in atmospheres con[5 Barker R. The reversibility of reaction CaCO, Ca0+CO[J]. J App! Chem Biotechnol, 1973, 23: 733--7421989,28(4):493-500[6]李振山蔡寧生,黃煜煜,等.CaO循環(huán)吸收CO2的實驗[11] Abanades J C. The maximum capture efficiency of CO2 using研究[J].燃燒科學與技術,2005,11(4):379-383.a carbonation/ calcination cycle of CaO/CaCO, [J].ChemLi Zhenshan, Cai Ningsheng, Huang Yuyu, et al. Exper-EngJ,2002,90(3):303-306mental research of Ca0 cyclic absorbing CO2[J]. Journal of [12] Li Z S, Cai N S, Huang YY, et al. Synthesis, experimentalCombustion Science and Technology, 2005, 11(4):379-studies, and Analysis of a new calcium-based carbon dioxide383(in Chinese)sorbent[J]. Energy and Fuels, 2005, 19(4): 1337--1452[7]陳小華,鄭瑛,鄭楚光,等.CaO再碳酸化的研究13]LiZs, CaiN s, Huang Y Y. Effect of preparation tempera[J].華中科技大學學報:自然科學版,2003,31(4):54ture on cyclic CO, capture and multiple carbonation-calcina[J]. Ind EnChen Xiaohua, Zheng Ying, Zheng Chuguang, et al. Re-Chem Res,2006,45(6):1911-1917search on the carbonation reaction of CaO[J]. J Huazhong [14] Abanades J C, Alvarez D. Conversion limits in the reactionand Technotogy: Nature Science Ediof CO2 with lime[ J]. Energy and Fuels, 2003, 17(2): 308tion,2003,31(4):54-55( in Chinese)[8] Bhatia S K, Perlmutter D D. Effect of the product layer on [15] Iyer M V, Gupta H, Sakadjian B B, et al. Multicyclic studythe kinetics of the CO,-Lime reaction[ J]. AICHE, 1983the simultaneous carbonation and sulfation of high-reactiv-29(1):7986ity Cao [J]. Ind Eng Chem Res, 2004, 43 (14): 3939--[9] German R M, Munir Z A. Surface area reduction during iso-3947thermal sintering[ J]. Journal of the American Ceramic Soci中國煤化工CNMHG