竹材熱解失重動(dòng)力學(xué)模型研究

第25卷第2期高?；瘜W(xué)工程學(xué)報No, 2 WoL 252011年4月Journal of Chemical Engineering of Chinese UniversitiesApr.2011文章編號:1003-901502011)022023105竹材熱解失重動(dòng)力學(xué)模型研究谷雙,陳紀忠(浙江大學(xué)化學(xué)工程與生物工程學(xué)系化學(xué)工程聯(lián)合國家重點(diǎn)實(shí)驗室,浙江杭州310027)摘要:采用熱重分析儀在溫度范圍313~873K和N2氣流保護下,以不同升溫速率(10、20、30、50Kmin2)對兩種竹材(#1和#2)進(jìn)行了熱解失重動(dòng)力學(xué)研究,并分別采用非模型動(dòng)力學(xué)法( model-free kinctics,MFK)和模型擬合法model- fitting method分析了熱重分析結果,結果表明,#樣品在不同升溫速率下的DTG曲線(xiàn)表現為單峰,非模型動(dòng)力學(xué)法計算獲得的表觀(guān)活化能E值隨樣品失重變化很小,說(shuō)明#1樣品的熱解過(guò)程比較簡(jiǎn)單,當采用動(dòng)力學(xué)模型進(jìn)行擬合時(shí),一級簡(jiǎn)單反應動(dòng)力學(xué)模型即可表征其熱解失重過(guò)程;而#2樣品的熱解行為相對比較復雜,其DTG曲線(xiàn)表現為雙峰,由非模型動(dòng)力學(xué)計算得到的E值隨a增大而單調上升,采用動(dòng)力學(xué)模型擬合可得到較復雜的獨立兩步反應模型第一步為一步簡(jiǎn)單反應,生成揮發(fā)物;第二步為兩步平行反應,分別生成揮發(fā)物和炭關(guān)鍵詞:竹材:熱重分析:熱解;動(dòng)力學(xué)中圖分類(lèi)號:Q5222;TQ0132文獻標識碼:AStudy on Pyrolysis Kinetics of BambooGU Shuang, CHEN Ji-zhong(State Key Laboratory of Chemical Engineering, Department of Chemical and Biological Engineering,Zhejiang University, Hangzhou 310027, China)Abstract: The pyrolysis of two different kinds of bamboo were carried out by thermogravimetry (TG)atdifferent heating rates(10, 20, 30 and 50 K min" )within the temperate range of 313-873 K and under Natmosphere; And the pyrolysis kinetics of bamboo were analyzed by both model-free-kinetics(MFK) andmodel-fitting-method. The results show that pyrolysis of sample I is a quite simple process with single DtGpeak, and the apparent activation energy e calculated by MFK varies little with the variation of conversion aThis means that the pyrolysis of sample 1 is a simple reaction and can be explained as one step first-orderreaction by model-fitting-method. The pyrolysis of sample #2 is a more complex process with two distinct DTGpeaks; Value of E calculated by MFK increases monotonously with the increase of a. According to themodel-fitting-method, the pyrolysis process of sample #2 may be described as two independent steps. The firststep is quite simple and can be described by one- step simple reaction model, while the second step is more complexand can be explained as existing two parallel first-order reactions forming volatile matter and char separatelKey words: bamboo; thermogravimetry analysis; pyrolysis; kinetics1引言竹子是重要的森林資源之一。在我國,竹材分布廣,資源蘊藏量大,是一種理想的可再生資源。對竹材熱解過(guò)程進(jìn)行系統、深入的研究,并深入認識竹材熱解的動(dòng)力學(xué)規律,將有利于開(kāi)發(fā)新型竹類(lèi)產(chǎn)品,對產(chǎn)品加工工藝設計、優(yōu)化也是十分必要的,具有極大的經(jīng)濟價(jià)值和社會(huì )價(jià)值。但是,目前國內關(guān)于竹產(chǎn)品研究和開(kāi)發(fā)很少,竹類(lèi)產(chǎn)品的生產(chǎn)也多停留在原始的初加工階段?，F有文獻對竹材熱解的研究多為竹醋液組分分析及竹炭性質(zhì)分析及應用等,對竹材熱解過(guò)程及其動(dòng)力學(xué)的中國煤化工收稿日期:20100907:修訂日期:20101025CNMHG作者簡(jiǎn)介:谷雙(1982-),女,遼寧確蘆島人,浙江大學(xué)博土生,通訊聯(lián)系人:陳紀忠,Em高?；瘜W(xué)工程學(xué)報2011年4月研究鮮見(jiàn)報道。杜瑛等對毛竹進(jìn)行主要化學(xué)成分分析,并給出了其慢速熱解的液體產(chǎn)物的詳細組成,探討了熱解溫度對產(chǎn)物組成的影響③;蔣劍春等對孝順竹進(jìn)行常規及快速熱解,研究了竹材在不同熱解條件下液相、固相及氣相產(chǎn)品的組成及得率差異;邵千釣等對兩種竹質(zhì)材料(毛竹和孝順竹的熱解失重行為進(jìn)行了研究,并采用一級反應過(guò)程描述其失重動(dòng)力學(xué),得到主要熱解反應的平均活化能;鄧天昇通過(guò)假設熱解反應級數,建立竹材熱解動(dòng)力學(xué)模型,分析了竹材熱解反應的反應級數,認為熱解反應級數與升溫速率有著(zhù)一定的關(guān)系,一般為15或2級;曾凱斌對竹材的非等溫動(dòng)力學(xué)過(guò)程進(jìn)行數據擬合,認為半纖維素在較低溫度下熱解,僅生成揮發(fā)性組分,表現為第一步熱解,可以用簡(jiǎn)單一級反應描述;纖維素和木質(zhì)素在較高溫度下熱解,表現為第二步熱解,是競爭的平行一級反應,產(chǎn)物分別是揮發(fā)性組分和竹炭。第一熱解過(guò)程與熱解升溫速率有關(guān),而第二熱解過(guò)程與升溫速率無(wú)關(guān)。目前,對于采用熱分析法研究固相反應動(dòng)力學(xué)已有多種數據處理方法。從操作方式上可以分為單掃描速率法和多重掃描速率法?，F有文獻對于生物質(zhì)熱解動(dòng)力學(xué)的研究多采用單一掃描速率法。研究表明,使用單一掃描速率法得到的動(dòng)力學(xué)結果并不可靠,往往不能反映固態(tài)反應的本質(zhì):應該采用多重掃描速率法來(lái)測定熱分析數據,并通過(guò)用等轉化率法確定活化能隨轉化率的變化情況。本文將分別采用非模型動(dòng)力學(xué)法和模型擬合法分析竹材熱解的非等溫動(dòng)力學(xué)特征,比較該兩種方法的差異和適用范圍,為熱重分析的動(dòng)力學(xué)分析與研究提供新思路。2實(shí)驗部分21實(shí)驗原料實(shí)驗原料毛竹分別采自浙江杭州老和山(#)及遂昌(#2),竹齡3~5年取竹子的下部材料為實(shí)驗樣品,預處理為粉末狀,在353K下干燥8h22實(shí)驗儀器與實(shí)驗過(guò)程熱解研究采用熱重分析儀(瑞士梅特利一托利多公司, TGA/SDTA851°,溫度靈敏度:001K,質(zhì)量靈敏度:0.01mg)。在高純氮氛下進(jìn)行熱解實(shí)驗,升溫速率分別為10、20、30及50Kmin-,熱解溫度為313-873K。固定天平載氣流量:40 mL min-,樣品保護氣流量:30 mL min223熱解動(dòng)力學(xué)數據處理方法當升溫速率恒定時(shí)dr動(dòng)力學(xué)方程的微分式表示為d7se"(-a)”積分式G(a)=l(A/B)exp(E/RT)dT=So(A/B)exp(E/RT)dT當采用“模型擬合法”計算時(shí),需事先假設熱解反應的動(dòng)力學(xué)模型,嘗試將各種動(dòng)力學(xué)模型函數的微分式∫(a)或積分式G(a)代入,選擇使方程獲得最佳線(xiàn)性者為所求動(dòng)力學(xué)模型。同時(shí)通過(guò)非線(xiàn)性擬合計算得到動(dòng)力學(xué)參數活化能E及速率指前因子A非模型動(dòng)力學(xué)MFK法采用的是 Vyazovka的等轉化率原理,即在一定的轉化率下反應速率只是反應溫度的函數,不必確定動(dòng)力學(xué)方程。因此,MFK法是在相同的轉化率a下,認為G(a)為定值,利用不同加熱速率B下所測得的TG曲線(xiàn)數據,計算獲得熱解反應的表觀(guān)活化能E隨轉化率a的變化曲線(xiàn)可見(jiàn),模型擬合法需要豐富的動(dòng)力學(xué)模型知識,且所九命的證化能E和速率指前因子A越多,因此模型擬合法是很復雜的動(dòng)力學(xué)解析工中國煤化工的動(dòng)力學(xué)方程,只是通過(guò)不同加熱速率下的TG數據計算獲得熱解反應的表CNMH化關(guān)系,解析工作簡(jiǎn)單易行。如果MFK法計算的表觀(guān)活化E隨轉化率a變化,說(shuō)明熱解反應是復雜的反應過(guò)程,反之第25卷第2期谷雙等熱解失重動(dòng)力學(xué)樸型研究表觀(guān)活化能E不隨轉化率a變化,熱解反應為簡(jiǎn)單反應過(guò)程。本文利用梅特利一托利多公司的熱解動(dòng)力學(xué)分析軟件 STARSW920,嘗試將非模型動(dòng)力學(xué)法和模型擬合法結合,解析竹材熱解反應的動(dòng)力學(xué)模型3竹材熱解動(dòng)力學(xué)1.13.1竹材熱解失重曲線(xiàn)不同升溫速率下1#和2#兩種竹材的熱重分- 50 K min -#1-+10 K- min-析實(shí)驗的TG和DTG曲線(xiàn)如圖1和圖2所示。由圖可知,在實(shí)驗條件下,隨著(zhù)溫度的升高,竹材熱解30Kmin4-#2→50Kmin-#2主要經(jīng)歷了三個(gè)階段第一階段是從加熱開(kāi)始到480K,主要是竹材中的水分蒸發(fā)即竹材的干燥階段,可認為竹材的化學(xué)成分沒(méi)有變化;第二階段從480K到約630K,為竹材熱解的主要失重過(guò)程;第三階段為熱解殘留圖1不同升溫速率下兩種竹材熱解的TG曲線(xiàn)Fig. I tG curves of both samples at different heating rate物的緩慢炭化,失重較慢。在計算竹材熱解動(dòng)力學(xué)時(shí),主要針對第二階段失重,并忽略第一步干燥階001210 Kmin"-#段對后期熱解的影響。.2竹材熱解動(dòng)力學(xué)模型3.2.1非模型動(dòng)力學(xué)法利用梅特利一托利多公司的熱解動(dòng)力學(xué)分析軟000450 K min件 STARSW920,以非模型動(dòng)力學(xué)法(MFK)對兩種竹材樣品在不同升溫速率下的失重行為進(jìn)行計算分析,得到的活化能EMK-a曲線(xiàn)如圖3所示。從圖3500600700800900可以看出,對#1竹材樣品,當a<0.55時(shí),即竹材處于熱解失重的主要階段,EwBK值約為140圖2不同升溫速率下兩種竹材熱解的DTG曲線(xiàn)k]. mol,且隨a變化很小,說(shuō)明熱解反應為一簡(jiǎn)單g 2 DTG curves of both samples at different heating rates反應過(guò)程,活化能E基本為一常數。當失重進(jìn)一步增加,竹材進(jìn)入緩慢炭化階段,Ewκ值隨a急劇增sample #1加,說(shuō)明熱解反應為復雜反應。而對#2竹材樣品來(lái)講,在整個(gè)熱解過(guò)程中計算得到的表觀(guān)活化能EwBK值均隨α增大而單調增加,且EM值均大于#竹材樣品的E值,說(shuō)明#2竹材的熱解反應為復雜反應過(guò)程。322模型擬合法0102030405060708對#竹材,熱解得到的DTG曲線(xiàn)為單峰,且由非模型動(dòng)力學(xué)法計算得到的E值隨a變化很小圖3兩種竹材的非模型函數法結果Fig 3 MFK results of both sample因此,采用一級簡(jiǎn)單反應動(dòng)力學(xué)模型表征竹材的熱解失重過(guò)程,得到的動(dòng)力學(xué)參數如表1所示表1不同升溫速率下#竹材熱解的動(dòng)力學(xué)參數可以看出,不同升溫速率下計算得到的動(dòng)力學(xué)Table 1 Kinetic parameters of sample #l atdifferent heating rates參數十分接近,其活化能值在120kmo左右,與非模型動(dòng)力學(xué)法得到的活化能EMx值略有差異主要是由計算過(guò)程中采用不同方法(分別為積分法中國煤化工m和微分法造成的,但這一差異在允許范圍之內。因THiCNMHG97X101.11×10此對#1樣品,熱解過(guò)程比較簡(jiǎn)單,升溫速率對熱11年4月解過(guò)程的影響并不明顯。表2不同升溫速率下#2樣品熱解的動(dòng)力學(xué)參數Table 2 Kinetic parameters of sample #2 at different beating rates對#竹材,DTG曲線(xiàn)為雙峰,非模型動(dòng)力學(xué)法計算得到的/Kmin194731.10x10l417.89表觀(guān)活化能E值隨a變化十分明15815157×101319907408×10266311.28×10°160.58249x10134931.16×101280.30206×10顯,因此參照文獻[刁],提出如35858x101318755287×101下動(dòng)力學(xué)模型第一步; Solid,h→ Volatile表3不同升溫速率對#樣品各步反應失重率的影響Table3 Effect of heating rates on weight loss of sample #2ate First DTG First step Second DTG Two steps Total weig第二步: Solor〃 olatile2n-peak/k weight loss/% peak/K weight loss/% loss/%54998602.156071554687101模型認為,分解過(guò)程發(fā)生的是前5606521.13610.157131后兩步獨立的失重反應。第一步606.82為一步簡(jiǎn)單反應,生成揮發(fā)物;第二步為兩步平行反應,分別生成揮發(fā)物和炭。第一步反應得到的揮發(fā)物得率即第一步反應的反應物占整體反應物的質(zhì)量分率作為模型參數之一,與三個(gè)反應各自的活化能指前因子同時(shí)由計算擬合得到。計算得到的不同升溫速率下反應的活化能、指前因子如表2所示由表2可以看出,升溫速率由10Kmn升至30Kmin4時(shí),各步反應的E及A均隨升溫速率的增大而增大;對第二步熱解過(guò)程,生成揮發(fā)性組分的E、A要大于生成焦炭的E、A,因此,升溫速率增大有利于生成揮發(fā)性組分,這與實(shí)際相符;但當升溫速率進(jìn)一步升至50Kmin時(shí),得到的E及A明顯下降。已有文獻指出,當升溫速率較快時(shí),熱解反應受傳熱影響,計算得到的活化能將明顯小于真實(shí)值由此可見(jiàn),升溫速率逐漸增大過(guò)程中,反應逐漸由動(dòng)力學(xué)控制轉變?yōu)閭鳠峥刂?。同時(shí)擬合得到的不同升溫速率下第一步獨立反應所得到的揮發(fā)物得率結果如表3所示?？梢钥闯?各步反應的DTG峰溫及得到的揮發(fā)物得率即反應失重率隨升溫速率大致呈上升趨勢,加快升溫速率有利于揮發(fā)類(lèi)產(chǎn)物的生成將兩種竹材的動(dòng)力學(xué)計算結果對比可以看出,升溫速率對其失重行為的影響表現出不同的規律。有文獻指出,相近產(chǎn)地的毛竹其纖維素、半纖維素、木質(zhì)素等含量變化很小,產(chǎn)地主要影響竹材的微量元素含量,由此造成兩種樣品對加熱條件的敏感程度不同。3.23結果討論綜合TG曲線(xiàn)及動(dòng)力學(xué)分析結果可以看出,對#1樣品,DTG曲線(xiàn)表現為單峰,升溫速率對其主要失重階段的熱解行為的影響并不明顯,對不同升溫速率下的熱解曲線(xiàn)進(jìn)行非模型動(dòng)力學(xué)計算得到的E值隨a變化很小,相應的,從模型動(dòng)力學(xué)的擬合結果可以看到,采用一級簡(jiǎn)單反應動(dòng)力學(xué)模型即可表征其熱解失重過(guò)程。而#2樣品的熱解行為相對比較復雜,其DTG曲線(xiàn)表現為雙峰,由非模型動(dòng)力學(xué)計算得到的E值隨α增大而單調上升,并且隨a增大,Ea曲線(xiàn)斜率亦逐漸增大。在采用模型動(dòng)力學(xué)進(jìn)行擬合時(shí),得到較復雜的獨立兩步反應模型,第一步為一步簡(jiǎn)單反應,生成揮發(fā)物;第二步為兩步平行反應,分別生成揮發(fā)物和炭。從模型動(dòng)力學(xué)的結果還可以看出升溫速率對#2樣品熱解的影響:隨升溫速率逐漸增大,反應逐漸由動(dòng)力學(xué)控制轉變?yōu)閭鳠峥刂?。由此可?dòng)力進(jìn)行分析得到的BQ曲線(xiàn)可明顯示反過(guò)10二濫m程是否可采用簡(jiǎn)單動(dòng)力學(xué)模型進(jìn)行擬合:而采用模型動(dòng)力學(xué)進(jìn)行計算,可進(jìn)一步得到關(guān)于反應的動(dòng)力學(xué)信0.5息,包括具體的動(dòng)力學(xué)參數活化能及指前因子的值并可在一定程度上揭示可能的反應機理33模型可靠性驗證中國煤化工曲線(xiàn)的比較由于模式函數法在計算過(guò)程中采用非線(xiàn)性擬合方CNMHPIG curve with法,因此將所求得的竹材熱解過(guò)程中各步的活化能和第25卷第2期雙等:竹材熱解失重動(dòng)力學(xué)模型研究指前因子回代,對竹材熱解過(guò)程進(jìn)行擬合,并與實(shí)驗TG值比較。由圖4可知,兩者吻合較好,表明所提出的竹材熱解機理模型很好地預測了竹材的熱解失重與溫度的關(guān)系。4結論本文分別采用非模型動(dòng)力學(xué)法( model-free kinetics,MFK)和模型擬合法(mode- fitting method)對兩種不同產(chǎn)地竹材在不同升溫速率(10、20、30、50Kmin2)下的熱解失重動(dòng)力學(xué)進(jìn)行研究,結果表明,不同產(chǎn)地毛竹表現出不同的熱解行為。通過(guò)非模型動(dòng)力學(xué)法計算分析,發(fā)現兩種竹材的熱解表觀(guān)活化能隨轉化率的變化關(guān)系也不同,#1竹材樣品的熱解表觀(guān)活化能隨轉化率變化很小,為簡(jiǎn)單反應過(guò)程,可用一級簡(jiǎn)單反應動(dòng)力學(xué)模型進(jìn)行描述;而#2竹材樣品的熱解表觀(guān)活化能隨轉化率增大而單調上升,為復雜反應,以模型擬合法得到較復雜的獨立兩步反應模型,第一步為一步簡(jiǎn)單反應,生成揮發(fā)物;第二步為兩步平行反應,分別生成揮發(fā)物和炭。符號說(shuō)明A一指前因子,s摩爾氣體常數,83145 J-mol 'K活化能,k一轉化率K反應速率常數升溫速率,Kmin-1溫度,K模式函數1/f(a)的積分函數參考文獻:[l] Zhejiang Forest Resources Monitoring Centre(浙江省森林資源監測中心), Forestry Survey and Design Institute of ZhejiangProvince(浙江省林業(yè)勘察設計院). The status and role of forestry in the economical and social development of Zhejiangprovince(林業(yè)在浙江國民經(jīng)濟與社會(huì )發(fā)展中的地位和作用)R]1996[2] ZHOU Shan(周珊),HUZe-you(胡澤友) YU Jing-quan(喻景權) Biodegradation of phenol wastewater by pseudomonas即.mobilized on bamboo-carbon(竹炭固定化假單胞菌處理含酚廢水的研究) J Chem Eng of Chinese Univ(高?；瘜W(xué)工程學(xué)報)2008,22(5)889-8943] DU Ying(杜瑛,QIwi-yan(齊衛艷) Miao Xia(菌霞),eal. Chemical analysis of pubescens and its pyrolysis(毛竹的主要化學(xué)成分分析與熱解 J Chem Ind and eng( Chinak化工學(xué)報,2004,15(12)209921024] JIANG Jian-chun(蔣劍春), DENG Xian-lun(鄧先倫) ZHANG Yan-ping(張燕萍,era. Study on the characteristic feature ofbamboo pyrolysis(材熱解特性研究 Chemistry and Industry of Forest Products(林產(chǎn)化學(xué)與工業(yè))2005,25(B10):15-18S] SHAO Qian-jun(邵千鈞, PENG Jin-xng彭錦星)ⅪUQun-ang(徐群芳,eral. Study on pyrolysis characteristics and kinetics ofbamboo under liner heating temperature(竹質(zhì)材料熱解失重行為及其動(dòng)力學(xué)研究)[Act" Energy Solaris Sinica(太陽(yáng)能學(xué)報)006,27()671466] DENG Tian-sheng(鄧天昇). Investigation of Bamboo Pyrolysis Process竹材熱解動(dòng)力學(xué)的研究)[D, Hangzhou(杭州) ZhejiangUniversity(浙江大學(xué))20041 ZENG Kai-bin(曾凱斌, JIANG Bin-bo(蔣斌波), CHEN Ji-zhong(陳紀忠) Kinetic of bamboo pyrolysis(材熱解過(guò)程的動(dòng)力學(xué)研究. J Chem Ind and Eng( China)化工學(xué)報,2006,57(2)318-323.[8] LU Zhen-ong(陸振榮). New progress in thermal analysis kinetics studies((熱分析動(dòng)力學(xué)的新進(jìn)展)[門(mén) Chinese Journal oforganic Chemistry(無(wú)機化學(xué)學(xué)報,1998,14(2):19-1269] REN Ning(任寧, ZHANG Jian-jun(張建軍) Progress in datum treatment methods of themal analysis kinetics(熱分析動(dòng)力學(xué)數據處理方法的研究進(jìn)展) Progress in Chemistry(化學(xué)進(jìn)展)2006,18(4)410416[10] Gabor Varhegyi, Michael Jerry Antal, Emma Jakab, er al. Kinetic modeling of biomass pyrolysis []. Journal of Analytical andApplied Pyrolysis, 1997, 42(1): 73-87中國煤化工CNMHG