顧曉利，馬 旭，李麗鮮，程康華，李忠正，肖 睿

（1.南京林業(yè)大學(xué) 化學(xué)工程學(xué)院，江蘇 南京210037；2.華南理工大學(xué) 制漿造紙工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，廣東 廣州510640；3.東南大學(xué) 能源與環(huán)境學(xué)院，江蘇 南京210096）

我國(guó)生物質(zhì)能儲(chǔ)量巨大，僅農(nóng)林廢棄物每年就約15億t，其中約5～7億t可資源化利用，折合3～4億t標(biāo)準(zhǔn)煤［1］。若能合理、充分地將這些生物質(zhì)資源以各種形式加以利用，可以大大減輕環(huán)境污染，實(shí)現(xiàn)CO2的零排放，對(duì)我國(guó)能源與環(huán)境可持續(xù)發(fā)展具有深遠(yuǎn)意義。由可再生的生物質(zhì)資源熱解氣化制取生物質(zhì)燃?xì)馐且环N較理想的替代化石燃?xì)獾姆桨?，具有污染排放少、產(chǎn)品能源品位高等優(yōu)勢(shì)，是生物質(zhì)清潔高效利用的有效途徑［2-3］。但生物質(zhì)燃?xì)獯嬖跓嶂档汀⒔褂秃看蟮热秉c(diǎn)，導(dǎo)致生物質(zhì)燃?xì)獾钠焚|(zhì)和化石燃?xì)獯嬖谳^大差距。

近年來，熱重-紅外聯(lián)用（TG-FT IR）技術(shù)逐漸成為生物質(zhì)熱解研究領(lǐng)域進(jìn)行動(dòng)態(tài)特性分析的新工具［4-6］。由于分析過程中不僅可獲得物質(zhì)熱分解的失重與溫度關(guān)系，還可實(shí)時(shí)檢測(cè)物質(zhì)熱分解氣相產(chǎn)物及部分液相產(chǎn)物的組成［7-8］，因而越來越受到研究者的重視。筆者以常見的廢松木屑為農(nóng)林廢棄物的代表底物，借助TG-FT IR手段探討其熱解反應(yīng)過程及產(chǎn)物分布規(guī)律，分別引入微分法和積分法推導(dǎo)出熱解動(dòng)力學(xué)方程，并通過固定床放大熱解實(shí)驗(yàn)進(jìn)行驗(yàn)證，為進(jìn)一步研究生物質(zhì)資源替代傳統(tǒng)石化資源技術(shù)建立數(shù)據(jù)基礎(chǔ)。

1 實(shí)驗(yàn)部分

1.1 試劑與儀器

實(shí)驗(yàn)所用的農(nóng)林廢棄物底物為松木屑，源自南京林業(yè)大學(xué)木材加工廠，其相關(guān)元素分析及工業(yè)分析結(jié)果見表1。高純N2和O2由南京通廣特種氣體股份有限公司提供。

表1 松木屑組成與特性Table 1 Composition and properties of pine sawdust

采用法國(guó)Setaram公司TG92型常壓熱重分析儀，以α-Al2O3為參比物，由程序設(shè)定升溫速率、終溫和終溫保留時(shí)間。實(shí)驗(yàn)過程中采用水冷方式使熱天平的工作溫度處在較低的狀態(tài)。采用德國(guó)Bruker公司Vector型傅里葉變換紅外光譜儀測(cè)定樣品的FT-IR譜。

采用Emerson公司組裝的多組分煤氣分析儀測(cè)定熱解氣體產(chǎn)物組成。該分析儀包括德國(guó)羅斯蒙特分析公司Rosemount NGA 2000和Rosemount Hydros 100 2種檢測(cè)器。NGA 2000主要由檢測(cè)CO2和CO的NGA 2000NDIR非散色紅外分析儀、檢測(cè)O2的NGA 2000PMD氧量分析儀、檢測(cè)CH4的NGA 2000HFID氫離子火焰檢測(cè)器以及檢測(cè)H2的Hydros 100熱導(dǎo)檢測(cè)器組成。

其他設(shè)備有北京七星華創(chuàng)股份有限公司D07-7B型質(zhì)量流量計(jì)（MFC）；上海同田生物科技股份有限公司TBP1022型高壓恒流泵，0～10mL／min，精確度1%。

1.2 熱重-紅外聯(lián)用分析法

實(shí)驗(yàn)前將松木屑底物置于80℃下烘干2h除去水分。稱取10mg底物置于熱重反應(yīng)裝置中，以10～80℃／min由環(huán)境溫度升至1000℃并保持15min，載氣為99.999%的高純N2，流量65mL／min。熱重分析儀與傅里葉變換紅外光譜儀之間通過保溫管連接，以保證反應(yīng)過程中產(chǎn)生的熱解蒸氣全部進(jìn)入傅里葉變換紅外光譜儀，在線檢測(cè)蒸氣組成。

1.3 固定床木屑熱解反應(yīng)與氣體產(chǎn)物分析

圖1為木屑固定床熱解氣化反應(yīng)實(shí)驗(yàn)裝置示意圖。固定床反應(yīng)器的材料為2520不銹鋼，長(zhǎng)450mm、內(nèi)徑32mm、外徑38mm。采用電熱爐加熱反應(yīng)器，由熱電偶控制反應(yīng)器內(nèi)溫度，由質(zhì)量流量計(jì)（MFC）控制氣體流量，工作溫度5～45℃，工作壓差范圍0.1～0.5MPa，最大工作壓力3MPa。由高壓恒流泵出來的水經(jīng)蒸氣發(fā)生器加熱至200℃進(jìn)入反應(yīng)器。

圖1 木屑固定床熱解氣化反應(yīng)實(shí)驗(yàn)裝置示意圖Fig.1 The schematic of fixed-bed device for sawdust pyrolytic gasification

向固定床反應(yīng)器內(nèi)通入流量為400mL／min的N2，同時(shí)加熱反應(yīng)器。當(dāng)爐內(nèi)溫度達(dá)到設(shè)定溫度時(shí)，打開恒流泵和蒸氣發(fā)生器，調(diào)節(jié)恒流泵的流量為0.6mL／min（熱態(tài)時(shí)約為720mL／min），與 N2充分混合后通入反應(yīng)器，混合氣體中水蒸氣體積分?jǐn)?shù)約占60%。準(zhǔn)確稱取2g松木屑底物，當(dāng)反應(yīng)器和蒸氣發(fā)生器的溫度均達(dá)到穩(wěn)定后，接入集氣袋，同時(shí)開始加料。加料完成后繼續(xù)通入N2和水蒸氣吹掃5min，打開O2閥門，控制O2流量為400mL／min并持續(xù)3min。此時(shí)爐內(nèi)殘?zhí)考按呋瘎┥戏e炭可以基本完全燃燒，并檢測(cè)到多余的O2。實(shí)驗(yàn)結(jié)束后關(guān)閉加熱爐，分析氣體組成。

2 結(jié)果與討論

2.1 木屑的TG-DTG曲線

圖2為不同升溫速率下木屑的TG-DTG曲線。一般來說，生物質(zhì)的熱解過程基本歸結(jié)于纖維素、半纖維素和木質(zhì)素3種主要高聚物熱解的綜合過程［9］。以20℃／min升溫速率為例，從圖2可以看出，木屑熱解主要分為水分析出、木質(zhì)素與纖維素的熱解和殘余焦炭的緩慢分解3個(gè)階段。其中，第2階段是木屑熱解的主反應(yīng)階段（226～471℃），373℃時(shí)達(dá)到最大失重峰，木屑熱解過程中70%左右的失重在這一階段完成。471℃以后基本穩(wěn)定，1000℃時(shí)熱裂解剩余質(zhì)量分?jǐn)?shù)為19.53%。木屑熱解主要是纖維素和木質(zhì)素的熱裂解，前者生成醛、酸等大分子有機(jī)成分和CO2、H2O等小分子產(chǎn)物。譚洪［9］推測(cè)，在木屑熱解過程中，纖維素會(huì)先發(fā)生解聚生成活性纖維素，隨著溫度的升高，纖維素分子結(jié)構(gòu)中各種鍵發(fā)生斷裂和重整的同時(shí)，開始釋放出大量的揮發(fā)份，此時(shí)纖維素大分子進(jìn)入到熱裂解的主要階段，繼而生成焦油、焦炭和小分子的輕質(zhì)氣體。而木質(zhì)素的熱解主要包括木質(zhì)素聚合體側(cè)鏈的（OH鍵、苯丙烷上的苯基C（C鍵等斷裂生成大量含氧化合物（水、CO、CO2）和甲氧基團(tuán)側(cè)鏈斷裂生成CH4等氣相產(chǎn)物［10-11］。隨著溫度的升高，低溫?zé)峤馍傻慕褂蜕疃葻崃呀?，生成CO等小分子化合物。

圖2 不同升溫速率下木屑的TG-DTG曲線Fig.2 TG-DTG curves of sawdust pyrolysis at different heating rates

2.2 木屑熱解產(chǎn)物的FT-IR分析

生物質(zhì)熱解產(chǎn)物復(fù)雜，因此其紅外光譜也復(fù)雜。紅外光譜中各種官能團(tuán)的歸屬主要有［9］：670～690cm-1處對(duì)應(yīng)CO2的C＝O鍵面內(nèi)變形振動(dòng)峰，2250～2400cm-1處對(duì)應(yīng)CO2的C＝O鍵收縮振動(dòng)峰，2015～2235cm-1處對(duì)應(yīng)CO的C≡O(shè)鍵收縮振動(dòng)峰，1650～1820cm-1處對(duì)應(yīng)醛、酮、酸類C＝O鍵伸縮振動(dòng)峰，2770～2845cm-1處對(duì)應(yīng)醛類C—H鍵伸縮振動(dòng)峰，1320～1440cm-1處對(duì)應(yīng)醇或酚類O—H鍵變形振動(dòng)峰，1045～1100cm-1處對(duì)應(yīng)芳香醚或酯類C—O—C伸縮振動(dòng)峰，3470～3760cm-1處對(duì)應(yīng)O—H鍵伸縮振動(dòng)峰。由于熱重分析儀和紅外分析儀之間采用200℃保溫管連接，生物質(zhì)熱解產(chǎn)生的酮、醛、酚等活性組分可能在保溫管內(nèi)聚合結(jié)焦或者發(fā)生二次裂解，因此對(duì)大分子的焦油進(jìn)行定量分析比較困難。

圖3為升溫速率20℃／min時(shí)木屑在不同溫度下熱解產(chǎn)物的FT-IR譜。由圖3可以看出，熱解液相產(chǎn)物主要有酚、醇、酯、酸、醛、酮、芳香烴和烷烴，其波數(shù)主要集中在1000～1800cm-1，然而各組分復(fù)雜，且大部分含量較低，屬于木屑熱解的副產(chǎn)物，其中醛／酮類物質(zhì)的含量相對(duì)較高。醛類物質(zhì)是導(dǎo)致生物質(zhì)深度熱解產(chǎn)物焦油生成的主要原因，且隨著熱解溫度的升高，醛酮類物質(zhì)深度熱解造成其含量不斷降低，焦油的生成量隨之增加［13］，因此控制醛／酮類物質(zhì)的含量對(duì)減少木屑熱解副產(chǎn)物焦油的生成很重要。

圖3 升溫速率20℃／min條件下木屑在不同溫度下熱解產(chǎn)物的紅外吸收光譜圖Fig.3 FT-IR spectra of the product from sawdust pyrolysis at different temperatures with 20℃／min heating rate

圖4 木屑?xì)庀酂峤猱a(chǎn)物的紅外吸光度隨熱解溫度的變化Fig.4 Infrared absorption value of gas products from sawdust pyrolysis vs pyrolysis temperature

顯而易見，木屑熱解的主產(chǎn)物為氣相產(chǎn)物，包括CO2、CH4、CO等。圖4為在20、40、80℃／min升溫速率條件下木屑熱解氣相產(chǎn)物的紅外吸光度隨熱解溫度的變化。一般來說，在一定濃度范圍內(nèi)，每種氣體的紅外吸光度值與氣體的濃度遵守朗伯-比爾定律［12］，即近似滿足正比關(guān)系。由圖4看到，隨著升溫速率的增加，同一熱解溫度下，氣相產(chǎn)物中CO、CO2和CH4的紅外吸光度也隨之增加；各升溫速率下，該3種產(chǎn)物的紅外吸光度隨熱解溫度變化的規(guī)律基本相同。木屑熱解過程中，CO紅外吸光度隨熱解溫度的升高存在2個(gè)峰值，第1個(gè)峰值出現(xiàn)在300～350℃，由纖維素和木質(zhì)素的熱解產(chǎn)生，第2個(gè)峰出現(xiàn)在700～800℃，可能是由于低溫時(shí)縮合結(jié)焦的大分子在高溫時(shí)再次發(fā)生裂解引起。作為永久性的小分子生成物，CO主要通過揮發(fā)份中不穩(wěn)定羰基的斷裂生成，木屑中除了本身含有大量羰基，熱裂解過程中普遍發(fā)生的脫水反應(yīng)也會(huì)形成羰基。可見CO可由一次反應(yīng)中產(chǎn)生，但更多的是由一次揮發(fā)分通過二次反應(yīng)形成。由于熱解反應(yīng)初期CO的形成顯得比較薄弱，因而通常將CO劃分為纖維素?zé)崃呀獾亩萎a(chǎn)物。CO2的紅外吸光度峰值出現(xiàn)在300～350℃，但從700℃以后沒有出現(xiàn)第2個(gè)峰值，僅有緩慢的增加。這是由于CO2主要在低溫段生成，高溫段的二次反應(yīng)對(duì)其貢獻(xiàn)不明顯，說明CO2的生成機(jī)理與CO并不相同［7］。通常認(rèn)為，CO2主要形成于纖維素低溫焦炭化期間，通過脫水纖維素的反應(yīng)生成［14］。此外，部分揮發(fā)分的二次反應(yīng)也會(huì)釋放少量CO2，但后者在纖維素的熱裂解過程中并不是主導(dǎo)反應(yīng)。而700℃以后CO2的緩慢增加主要由木質(zhì)素的熱解造成［15］，雖然木質(zhì)素苯環(huán)中的C—C鍵很牢固，但連接各個(gè)單體苯環(huán)上的側(cè)鏈和氧鍵橋相對(duì)脆弱，受熱易發(fā)生斷裂，形成活潑的含苯環(huán)自由基，易與其他分子或自由基發(fā)生縮合反應(yīng)生成結(jié)構(gòu)更穩(wěn)定的易于結(jié)焦的大分子，同時(shí)釋放了大量的CO、CO2及水分子［16］。CH4的紅外吸光度峰值只有1個(gè)，位于350～400℃，明顯高于CO、CO2出現(xiàn)峰值的溫度，比木屑熱解的TG-DTG曲線出現(xiàn)最大失重峰的溫度要高，可見生成CH4所需的能量比 CO 和 CO2高［17］。

2.3 木屑熱解動(dòng)力學(xué)模型的計(jì)算

生物質(zhì)熱解屬于復(fù)雜的反應(yīng)過程，求解其熱解動(dòng)力學(xué)參數(shù)的方法很多，其中以Kissinger的微分法［18］和Ozawa的積分法［19］最具代表性。一般來說，忽略反應(yīng)溫度對(duì)活化能的影響，并假設(shè)其符合簡(jiǎn)單動(dòng)力學(xué)方程，生物質(zhì)熱解動(dòng)力學(xué)反應(yīng)方程可由式（1）表示。式（1）中的α，k，f（α）分別由式（2）～（4）計(jì)算。

式（1）～（4）中，α為反應(yīng)度或失重率；k為反應(yīng)速率常數(shù)，min-1；m0、mt和mf分別為生物質(zhì)的原始質(zhì)量、在熱解時(shí)某一時(shí)刻t的質(zhì)量和熱解結(jié)束后的質(zhì)量，mg；t為反應(yīng)時(shí)間，min；T為反應(yīng)溫度，℃；A為頻率因子，min-1；E為活化能，kJ／mol；R為理想氣體常數(shù)，8.314×103kJ／（mol·K）；f（α）為與反應(yīng)機(jī)理相關(guān)的函數(shù)；n為反應(yīng)級(jí)數(shù)，生物質(zhì)熱解反應(yīng)通常為一級(jí)反應(yīng)，n＝1。

將式（3）和（4）代入到式（1）可得式（5）。

在底物熱解過程的TG-FTIR測(cè)定中設(shè)定為線性升溫速率，將升溫速率β＝dT／dt代入式（5）得式（6）。

式（6）中，熱解溫度T的表達(dá)式為T＝T0＋β·t，T0為熱解初始溫度，℃。

分別采用Kissinger的微分法和Ozawa的積分法對(duì)式（6）進(jìn)行數(shù)學(xué)處理，可得式（7）和式（8）。

式（7）中，Tp為失重DTG最大峰值對(duì)應(yīng)的溫度，℃；n＝1時(shí)，g（α）＝-ln（1-α）。進(jìn)而可以擬合出熱解反應(yīng)過程的活化能E和頻率因子A。

2.3.1 Kissinger法

木屑熱解的動(dòng)力學(xué)參數(shù)可由Kissinger算法推出的式（7）計(jì)算得來。由不同升溫速率條件下的Tp值的變化作出ln（β／Tp2）-1／Tp曲線，如圖5所示。由圖5得出，該擬合線的斜率為-15.65331×103，進(jìn)而推出活化能E＝130.14kJ／mol；截距為16.62308，推得頻率因子lnA＝26.28min-1。

圖5 基于Kissinger算法的木屑熱解動(dòng)力學(xué)擬合曲線Fig.5 Kinetic linearization curve of sawdust based on Kissinger method

2.3.2 Ozawa法

由Ozawa算法推出的式（8）可計(jì)算出木屑熱解的動(dòng)力學(xué)參數(shù)。由不同升溫速率條件下的各參數(shù)值的變化作出lnβ-1／Tp曲線，如圖6所示。由圖6得出，該擬合線的斜率為-16.85612×103，進(jìn)而推出活化能E＝133.21kJ／mol；截距為31.42234，推得頻率因子lnA＝26.47min-1。與Kissinger算法得出的動(dòng)力學(xué)參數(shù)基本一致。

圖6 基于Ozawa算法的木屑熱解動(dòng)力學(xué)擬合曲線Fig.6 Kinetic linearization curve of sawdust based on Ozawa method

2.4 固定床木屑熱解反應(yīng)的影響因素

采用固定床作反應(yīng)器，按1.3節(jié)設(shè)定的工況條件進(jìn)行木屑固定床熱解實(shí)驗(yàn)，其中熱解溫度選擇在680～870℃范圍變化。此部分可視作木屑熱解的放大過程。

2.4.1 熱解溫度對(duì)木屑熱解產(chǎn)物分布的影響

圖7 熱解溫度（T）對(duì)木屑熱解反應(yīng)產(chǎn)物分布的影響Fig.7 Influence of pyrolytic temperature（T）on product distribution of sawdust pyrolysis

圖7給出了不同熱解溫度下木屑熱解氣化過程中氣相、固相及焦油的質(zhì)量分?jǐn)?shù)。由圖7可以看出，680℃時(shí)氣相產(chǎn)物的質(zhì)量分?jǐn)?shù)為34.0%，隨著熱解溫度的升高，氣相產(chǎn)物所占比例逐漸增加，870℃時(shí)氣相產(chǎn)物質(zhì)量分?jǐn)?shù)達(dá)到67.5%。這是由于氣相產(chǎn)物主要由焦油和焦炭?jī)刹糠洲D(zhuǎn)化而成，低溫時(shí)隨著溫度的升高，焦油分解生成氣相產(chǎn)物，而當(dāng)溫度達(dá)到810℃以后焦油分解變得很緩慢，但焦炭的氣化反應(yīng)加快，生成大量H2、CO等氣體。

2.4.2 熱解溫度對(duì)木屑熱解氣相產(chǎn)物分布的影響

圖8為熱解溫度對(duì)木屑熱解氣相產(chǎn)物分布的影響。從圖8可知，680℃時(shí)氣相中CO、CO2、CH4、H2的體積分?jǐn)?shù)分別為48.4%、12.9%、10.7%和25.3%。隨著溫度的升高，CO體積分?jǐn)?shù)迅速下降，840℃時(shí)達(dá)到最低；H2體積分?jǐn)?shù)迅速上升，840℃時(shí)達(dá)到最高值；CO2體積分?jǐn)?shù)稍微上升，但基本保持在12%左右。680℃時(shí)，由于溫度較低，水蒸氣催化裂解焦油或水蒸氣與焦炭的反應(yīng)都不明顯；720℃時(shí)，水蒸氣的作用開始顯現(xiàn)，即發(fā)生式（9）的反應(yīng)，但生成的CO繼續(xù)與水蒸氣反應(yīng)生成H2（見式（10）），因此該溫度下CO體積分?jǐn)?shù)下降，而CO2和H2體積分?jǐn)?shù)上升；當(dāng)溫度繼續(xù)升高，發(fā)生式（11）的反應(yīng)，CO體積分?jǐn)?shù)開始上升。因此，生物質(zhì)熱解氣化較理想的溫度為700℃以上，這一結(jié)論與吳創(chuàng)之等［20］得到的結(jié)果相一致。

圖8 熱解溫度（T）對(duì)木屑熱解反應(yīng)氣相產(chǎn)物分布的影響Fig.8 Influence of pyrolytic temperature（T）on gas product distribution of sawdust pyrolysis

從圖8還可以看出，CH4的體積分?jǐn)?shù)從680℃的10.7%緩慢下降至840℃的6.9%，840℃以后CH4體積分?jǐn)?shù)繼續(xù)減少。這可能是因?yàn)闅饣磻?yīng)中（見式（9）～（13）），反應(yīng)式（12）和（13）是生成 CH4的重要反應(yīng)，式（12）適宜的反應(yīng)溫度一般認(rèn)為在800℃左右［15］，隨著溫度的升高，逆反應(yīng)速率上升很快，不利于CH4的生成；式（13）既是氣化過程中CH4生成的一個(gè)主要反應(yīng)，也是氣化后對(duì)生成氣進(jìn)行甲烷化反應(yīng)以增加CH4收率的重要反應(yīng)，隨著溫度的上升，其逆反應(yīng)加速，不利于CH4的生成。

3 結(jié) 論

（1）木屑熱解過程由水分析出，木質(zhì)素、纖維素的熱解和殘余焦炭的緩慢分解3個(gè)階段組成。其中第2階段是木屑熱解的主反應(yīng)階段，木屑熱解過程中約70%的失重在這一階段發(fā)生。熱解主反應(yīng)階段的失重主要由纖維素的熱裂解，生成醛、酸等大分子有機(jī)成分和CO2、H2O等小分子產(chǎn)物引起。

（2）木屑熱解產(chǎn)物以氣體產(chǎn)物為主，包括CO2、CH4、CO等。液相產(chǎn)物成分復(fù)雜，包括酚、醇、酯、酸、醛、酮等芳香烴和烷烴，且含量較低，僅醛、酮類物質(zhì)含量相對(duì)較高。隨熱解過程的加劇，上述液相產(chǎn)物深度熱解，最終生成氣態(tài)小分子化合物、液相焦油或固相焦炭。

（3）將熱解過程視為一級(jí)反應(yīng)過程，分別采用Kissinger微分法和Ozawa積分法確定了木屑熱解動(dòng)力學(xué)參數(shù)，得到熱解反應(yīng)活化能E分別為130.14和 133.21kJ／mol， 頻 率 因 子 lnA為 26.28 和26.47min-1。2種算法的結(jié)果吻合。

（4）木屑熱解較理想的熱解溫度應(yīng)控制在700℃以上，此時(shí)熱解產(chǎn)物以CO、CO2、CH4、H2氣體為主，且隨著溫度的升高，CO收率迅速下降，H2收率迅速上升，CO2和CH4收率基本不變。木屑固定床熱解實(shí)驗(yàn)結(jié)果與TG-FT IR的結(jié)論基本一致。

［1］中國(guó)可再生資源發(fā)展戰(zhàn)略研究項(xiàng)目組.中國(guó)可再生能源發(fā)展戰(zhàn)略研究叢書綜合卷［M］.北京：中國(guó)電力出版社，2008.

［2］RAO R T，SHARMA A.Pyrolysis rates of biomass materials［J］.Energy，1998，23（11）：973-978.

［3］YANG H，YAN R，CHIN T，et al.Thermogravimetric analysis fourier transform infrared analysis of palm oil waste pyrolysis［J］.Energy Fuel，2004，18（6）：1814-1821.

［4］LIOU T H.Pyrolysis kinetics of electronic packaging material in a nitrogen atmosphere［J］.Journal of Hazardous Materials，2003，103（1）：107-123.

［5］FONT R，F(xiàn)ULLANA A，CONESA J.Kinetic models for the pyrolysis and combustion of two types of sewage sludge［J］.Journal of Analytical and Applied Pyrolysis，2005，74（2）：429-438.

［6］OSVALDA S.Kinetics of pyrolysis，combustion and gasification of three biomass fuels［J］.Fuel Process Technology，2007，88（1）：87-97.

［7］LI S，LYONS H J，BANYASZ J，et al.Real-time evolved gas analysis by FT IR method：An experimental study of cellulose pyrolysis［J］.Fuel，2001，80（12）：1809-1817.

［8］BECIDAN M， SKREIBERG Φ， HUSTAD J E.Products distribution and gas release in pyrolysis of thermally thick biomass residues samples［J］.Journal of Analytical and Applied Pyrolysis，2007，78（1）：207-213.

［9］譚洪.生物質(zhì)熱解機(jī)理實(shí)驗(yàn)研究［D］.杭州：浙江大學(xué)，2005.

［10］KUOFOPANOS C，MASCHIO G，LUCCHESI A.Kinetics modeling of the pyrolysis of biomass and biomass components［J］.Canadian Journal of Chemical Engineering，1989，67（1）：75-84.

［11］NGUYEN T，ZAVARIN E，BARRALL E M.Thermal analysis of lignocellulosic materials ［J］. Polymer Reviews，1981，20（1）：1-65.

［12］華中師范大學(xué)等編.分析化學(xué)（上冊(cè)）［M］.北京：高等教育出版社，2005

［13］LU Q，TANG Z，ZHANG Y，et al.Catalytic upgrading of biomass fast pyrolysis vapors with Pd／SBA-15 catalysts［J］.Ind Eng Chem Res，2010，49（6）：2573-2580.

［14］王樹榮，劉倩，駱仲泱，等.基于熱重紅外聯(lián)用分析的纖維素?zé)崃呀鈾C(jī)理研究［J］.浙江大學(xué)學(xué)報(bào)（工學(xué)版）2006，40（7）：1154-1158（WANG Shurong，LIU Qian，LUO Zhongyang，et al.Mechanism study of cellulose pyrolysis using thermogravimetric analysis coupled with infrared spectroscopy［J］.Journal of Zhejiang University（Engineering Science），2006，40（7）：1154-1158.）

［15］何光瑩.生物質(zhì)催化氣化制取富甲烷生物質(zhì)燃?xì)庠囼?yàn)研究［D］.南京：東南大學(xué)，2011.

［16］李雪瑤，應(yīng)浩.生物質(zhì)熱解氣化機(jī)理研究進(jìn)展［J］.工程熱物理學(xué)報(bào)，2006，27（2）：45-50.（LI Xueyao，YING Hao.Research advance on mechanism of biomass pyrolysis and gasification ［J］. Advances in Fine Petrochemicals，2006，27（2）：45-50.）

［17］BROIDO A，NOLSON M A.Char yield on pyrolysis of cellulose ［J］. Combustion Flame， 1975， 24（2）：263-268.

［18］KISSINGER H E.Reaction kinetics in differential thermal analysis［J］. Anal Chem，1957，29（11）：1702-1706.

［19］OZAWA T. A new method of analyzing thermogravimetric data［J］.Bull Chem Soc Jpn，1965，38（11）：1881-1886.

［20］吳創(chuàng)之，徐冰燕，羅曾凡，等.生物質(zhì)中熱值氣化技術(shù)的分析及探討［J］.煤氣與熱力，1995，15（2）：8-14.（WU Chuangzhi，XU Bingyan，LUO Zengfan，et al.Analysis of biomass gasification for MHV fuel gas［J］.Gas ＆ Heat，1995，15（2）：8-14.）