杜瑞嶺， 吳 鏗， 巢昌耀， 徐大安， 張 黎， 張 兵

(鋼鐵冶金新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(北京科技大學(xué))， 100083 北京)

熱重分析法研究煤粉燃燒過程動(dòng)力學(xué)的Arrhenius方程修正式

杜瑞嶺， 吳 鏗， 巢昌耀， 徐大安， 張 黎， 張 兵

(鋼鐵冶金新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(北京科技大學(xué))， 100083 北京)

煤粉的燃燒過程常被作為能量的提供者應(yīng)用于不同工業(yè)，為研究升溫速率對(duì)煤粉燃燒過程動(dòng)力學(xué)的影響. 本文采用熱重分析法對(duì)20 ℃/min、25 ℃/min、30 ℃/min和35 ℃/min升溫速率下的煤粉燃燒過程進(jìn)行分析. 結(jié)果表明：隨升溫速率的提高，煤粉燃燒過程有明顯的熱滯后現(xiàn)象. 根據(jù)煤粉燃燒過程的特點(diǎn)，以反應(yīng)速率曲線波谷對(duì)應(yīng)的溫度點(diǎn)，將著火點(diǎn)到燃盡點(diǎn)的燃燒過程劃分為兩個(gè)階段，并分別采用界面化學(xué)反應(yīng)模型和內(nèi)擴(kuò)散模型來描述相應(yīng)階段的動(dòng)力學(xué)過程. 由所獲得動(dòng)力學(xué)參數(shù)可知，不同升溫速率下活化能Eai和指前因子Ai的動(dòng)力學(xué)補(bǔ)償效應(yīng)可表示為lnAi=aEai+b，升溫速率β對(duì)活化能Eai的影響可表示為Eai=△Ealnβi+Ea0. 隨將常用的Arrhenius方程lnki=-Ea0/RT+lnA0修正為lnki=-Ea0/RT+△Ea(a-1/RT)lnβi+lnA0來描述升溫速率β對(duì)煤粉燃燒過程反應(yīng)速率常數(shù)k的影響. 而后，采用10 ℃/min、15 ℃/min、40 ℃/min和45 ℃/min升溫速率下的煤粉燃燒試驗(yàn)對(duì)Arrhenius方程修正式的外推性進(jìn)行驗(yàn)證，效果良好. 因此，Arrhenius方程修正式不僅能很好地描述升溫速率β對(duì)煤粉燃燒過程化學(xué)反應(yīng)速率常數(shù)k的影響，而且還具有一定的外推性.

煤粉燃燒；升溫速率；Arrhenius方程修正式；動(dòng)力學(xué)參數(shù)；非等溫體系

煤粉燃燒過程常被作為能量的提供者應(yīng)用于不同工業(yè)[1]. 研究發(fā)現(xiàn)升溫制度對(duì)煤粉燃燒過程有重要的影響. 然而，前人關(guān)于升溫速率對(duì)煤粉燃燒過程的影響主要是定性分析，定量研究較少[2-3]. 參考文獻(xiàn)[4]提出一種預(yù)測(cè)不同升溫速率下煤粉熱解過程動(dòng)力學(xué)常數(shù)的方法. 其表達(dá)式為k1/k2=(m1/m2)n，其中，k為熱解反應(yīng)速率常數(shù)，m為升溫速率，指數(shù)n為升溫速率對(duì)煤熱解反應(yīng)速率的影響程度. 然而，關(guān)系式只能簡(jiǎn)單計(jì)算k的平均值. 忽略溫度對(duì)速率常數(shù)的影響，沒有給出速率常數(shù)k與溫度和升溫速率的關(guān)系式. 因此，表達(dá)式不能直接用來描述煤粉燃燒動(dòng)力學(xué)過程. 所以本文重點(diǎn)研究升溫速率對(duì)煤粉燃燒過程動(dòng)力學(xué)的影響.

在等溫條件下，Arrhenius方程常用來描述反應(yīng)速率常數(shù)與溫度之間的關(guān)系. 隨著熱分析動(dòng)力學(xué)的發(fā)展，Arrhenius方程被應(yīng)用到非等溫領(lǐng)域. 然而，無論采用哪種熱力學(xué)分析方法(如：模型擬合法[5]，等轉(zhuǎn)化率法[6]，分布活化能法(DAEM)[7])，在相同溫度下活化能Ea和速率常數(shù)k都會(huì)隨著升溫速率β的變化而改變. 因此，本文提出一個(gè)Arrhenius方程修正式來描述升溫速率β對(duì)反應(yīng)速率常數(shù)k的影響.

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

實(shí)驗(yàn)過程中所用煤粉試樣取自某鋼鐵廠1#高爐的噴吹煤粉，其工業(yè)分析和元素分析見表1，粒度分布見表2.

表1 煤粉工業(yè)分析與元素分析結(jié)果

表2煤粉粒度分布

Tab.2 Particle size distribution of pulverized-coal used in the experiment %

實(shí)驗(yàn)儀器采用北京恒久科學(xué)儀器廠生產(chǎn)HCT-1綜合熱分析儀，按照設(shè)定程序自動(dòng)升溫，期間實(shí)驗(yàn)數(shù)據(jù)由計(jì)算機(jī)自動(dòng)采集. 每次實(shí)驗(yàn)所用試樣質(zhì)量為10±0.2 mg. 將煤粉試樣放入到Al2O3坩堝內(nèi)，由室溫25 ℃以10 ℃/min升溫至105 ℃恒溫10 min，以除去樣品中的吸附水；以不同的升溫速率(20 ℃/min、25 ℃/min、30 ℃/min和35 ℃/min)加熱到950 ℃. 在整個(gè)實(shí)驗(yàn)過程中采用氣流量150 mL/min的空氣作為載氣.

3.4小結(jié)中為驗(yàn)證Arrhenius方程修正式在其它升溫速率下的外推性，又對(duì)升溫速率10 ℃/min、15 ℃/min、40 ℃/min和45 ℃/min煤粉燃燒過程進(jìn)行實(shí)驗(yàn).

2 結(jié)果與討論

2.1升溫速率對(duì)煤粉燃燒過程的影響

本文采用熱重分析法對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析處理，圖1、2分別給出不同升溫速率下煤粉燃燒過程的轉(zhuǎn)化率(α)和反應(yīng)速率(dα/dT)曲線圖. 其中，反應(yīng)速率曲線是轉(zhuǎn)化率曲線的一次微分曲線.

根據(jù)實(shí)際應(yīng)用情況，本文主要研究從著火點(diǎn)到燃盡點(diǎn)(Ti～Tb)的燃燒過程動(dòng)力學(xué). 由圖1、2可知，曲線形狀相似，隨升溫速率提高整個(gè)曲線向右偏移. 整個(gè)燃燒過程表現(xiàn)出明顯的滯后性[8]. 從圖2可看出，每條反應(yīng)速率曲線都有兩個(gè)峰值和一個(gè)波谷[9]，從燃燒反應(yīng)宏觀動(dòng)力學(xué)方面入手，現(xiàn)象被認(rèn)為是由于反應(yīng)機(jī)理和控速環(huán)節(jié)的改變而導(dǎo)致. 因此，以燃燒反應(yīng)速率曲線中波谷值所對(duì)應(yīng)溫度點(diǎn)將燃燒過程中著火點(diǎn)到燃盡點(diǎn)的溫度區(qū)間(Ti～Tb)分為兩個(gè)階段，結(jié)果見表3. 由表3可知，在非等溫條件下，不同升溫速率下溫度區(qū)間劃分也不相同.

圖1 不同升溫速率下的煤粉燃燒過程轉(zhuǎn)化率曲線

Fig.1 Fractional conversion curves of pulverized-coal combustion at different heating rates

圖2不同升溫速率下的煤粉燃燒過程反應(yīng)速率曲線(dα/dT)

Fig.2 Reaction rate curves of pulverized-coal combustion at different heating rates (dα/dT)

表3 不同升溫速率下的反應(yīng)溫度區(qū)間、峰值溫度、特征溫度點(diǎn)

注：Ti為著火點(diǎn)溫度，T1為第一個(gè)峰值所對(duì)應(yīng)溫度，Tv波谷所對(duì)應(yīng)溫度，T2為第二個(gè)峰值所對(duì)應(yīng)溫度，Tb為燃盡溫度.

第一階段，主要是揮發(fā)分的析出和燃燒過程[10]. 例如：液體和氣體產(chǎn)物大量生成、脫揮發(fā)分作用，以及揮發(fā)分的燃燒過程. 此外，由表3可看出，粉煤著火點(diǎn)溫度Ti隨升溫速率增加而增加，而不是一個(gè)常數(shù). 因此，著火點(diǎn)溫度不僅與煤的性能有關(guān)，而且還受升溫速率的影響.

第二階段，主要是焦炭的氧化過程. 分為3個(gè)步驟：空氣中的氧擴(kuò)散到碳表面，碳和氧的表面化學(xué)反應(yīng)，氣體反應(yīng)產(chǎn)物的擴(kuò)散過程[11]. 隨著反應(yīng)進(jìn)行，溫度逐漸升高，化學(xué)反應(yīng)速率逐漸加快并超過擴(kuò)散速率. 因此，控速環(huán)節(jié)的改變是圖2中出現(xiàn)雙峰的主要原因.

2.2Arrhenius方程修正式

在無限小的時(shí)間范圍內(nèi)，非等溫過程被看成等溫過程. 其反應(yīng)速率為

dα/dt=k(T)f(α) .

(1)

式中：f(α)為反應(yīng)模型；α為轉(zhuǎn)化率；t為時(shí)間，單位是s；k(T)為反應(yīng)速率常數(shù)，單位是s-1；T為溫度，單位是K.

樣品的轉(zhuǎn)化率α為

α=(m0-mt)/(m0-mfinal) .

(2)

式中：m0為樣品初始質(zhì)量，單位是g；mt為樣品在t時(shí)刻的質(zhì)量，單位是g；mfinal為樣品最終反應(yīng)的剩余質(zhì)量，單位是g.

將Arrhenius方程k(T)=Aexp(-Ea/RT)和升溫速率的關(guān)系式dT=βdt代入到式(1)中. 其近似積分式為[12]

(3)

式中：G(α)為f(α)-1的積分式；Ea為活化能，單位是J·mol-1；A為指前因子，單位是s-1；R為理想氣體常數(shù)，值為8.314 J·(K·mol)-1；β為升溫速率，單位是K·s-1.

在同一升溫速率下，升溫速率β為一個(gè)定值. 對(duì)式(3)取對(duì)數(shù)得

(4)

將動(dòng)力學(xué)機(jī)理函數(shù)G(α)代入到式(4). 根據(jù)ln[G(α)/T2]與1/T的線性關(guān)系，可求出相應(yīng)的活化能和指前因子.

一般情況下，不同升溫速率下的活化能Ea和指前因子A之間存在著動(dòng)力學(xué)補(bǔ)償效應(yīng)[13]，為

lnAi=aEai+b.

(5)

式中下標(biāo)i代表不同的升溫速率. 因此，指前因子為Ai=exp(aEai+b)，Arrhenius方程改寫為

(6)

[14]從分布活化能模型中推導(dǎo)出活化能Ea與升溫速率β的關(guān)系式，dβ/dT=β(2/T+Ea/RT2). 對(duì)其進(jìn)行積分得Ea=RTln(T2)-RTln(β). 本文采用不同的動(dòng)力學(xué)模型G(α)來描述煤燃燒過程的相應(yīng)動(dòng)力學(xué)區(qū)間，研究發(fā)現(xiàn)升溫速率對(duì)活化能的影響可表示為Eai=△Ealnβi+Ea0. 因此，式(6)被修正為

(7)

式中：A0=exp(aEa0+b)，根據(jù)lnAi和Eai與Eai和lnβi擬合曲線的斜率和截距可求出參數(shù)a、b、Ea0和△Ea. 與Arrhenius方程ki(T)=Aiexp(-Eai/RT)相比，指前因子被修正為lnAi=a△Ealnβi+lnA0.

對(duì)式(7)簡(jiǎn)化后取對(duì)數(shù)，得Arrhenius方程修正式

lnki=-Ea0/RT+△Ea(a-1/RT)lnβi+lnA0.

(8)

為簡(jiǎn)化k與β的關(guān)系式，由式(7)得

(9)

將式(9)與參考文獻(xiàn)[4]中經(jīng)驗(yàn)關(guān)系式k1/k2=(m1/m2)n相比可知，指數(shù)n被修正為△Ea(a-1/RT). 即在不同溫度下反應(yīng)速率常數(shù)k的平均值被修正成一個(gè)瞬時(shí)值.

2.3相應(yīng)動(dòng)力學(xué)參數(shù)求解

粉煤燃燒過程屬于氣固多相反應(yīng). 兩個(gè)反應(yīng)區(qū)間的反應(yīng)機(jī)理不同，如果只采用單一模型來擬合整個(gè)燃燒過程顯然是不合理的[15]. 因此本文采用分段的方法來擬合煤粉燃燒的動(dòng)力學(xué)過程. 在第一階段，由于溫度相對(duì)較低，反應(yīng)產(chǎn)物較少，氣體反應(yīng)物和生成物易穿過固體產(chǎn)物層. 燃燒過程主要控速環(huán)節(jié)是化學(xué)反應(yīng)過程. 因此，第一階段采用界面化學(xué)反應(yīng)模型來進(jìn)行擬合[16]. 隨溫度升高，反應(yīng)進(jìn)行，界面化學(xué)反應(yīng)速率逐漸升高. 與此同時(shí)，固體產(chǎn)物逐漸積累，氣體反應(yīng)物和產(chǎn)物的擴(kuò)散過程受到阻礙. 控速環(huán)節(jié)由界面化學(xué)反應(yīng)變?yōu)閿U(kuò)散過程. 因此，第二階段采用內(nèi)擴(kuò)散模型來進(jìn)行擬合. 對(duì)應(yīng)的機(jī)理函數(shù)見表4.

表4兩個(gè)反應(yīng)階段的動(dòng)力學(xué)機(jī)理函數(shù)

Tab.4 Kinetic mechanism functions for two reaction regions

反應(yīng)階段f(α)G(α)第一階段f1(α)=3(1-α)2/3G1(α)=1-(1-α)1/3第二階段f2(α)=3/2[(1-α)-1/3-1]-1G2(α)=1-2/3α-(1-α)2/3

根據(jù)表3中對(duì)應(yīng)的溫度區(qū)間和表4中相應(yīng)的動(dòng)力學(xué)機(jī)理函數(shù)，將不同反應(yīng)區(qū)間的實(shí)驗(yàn)點(diǎn)依據(jù)式(4)進(jìn)行擬合，結(jié)果見圖3、4.

圖3 不同升溫速率下第一階段的擬合曲線

Fig.3 Model fitting curves for first region at different heating rates

圖4 不同升溫速率下第二階段的擬合曲線

Fig.4 Model fitting curves for second region at different heating rates

從圖3、4中可看出，實(shí)驗(yàn)點(diǎn)的線性關(guān)系很好. 由表5可知，所有的相關(guān)系數(shù)均超過0.99. 這表明，界面化學(xué)反應(yīng)模型和內(nèi)擴(kuò)散模型能很好地描述煤粉燃燒的兩個(gè)反應(yīng)階段. 根據(jù)圖3、4中擬合直線的斜率K和截距C，求得相應(yīng)的活化能Ea=-RK和指前因子A=-βKexp(C)，結(jié)果見表5.

由表5可知，在相同控速環(huán)節(jié)內(nèi)，活化能隨升溫速率增加而減小. 煤燃燒過程中存在著許多放熱反應(yīng). 在非等溫條件下，升溫速率越快，單位溫度內(nèi)積累的熱量越多. 因此，快速升溫能促進(jìn)燃燒反應(yīng)，致使活化能隨升溫速率提高而降低[17]. 升溫速率對(duì)活化能Ea和指前因子A有重要的影響. 使得反應(yīng)速率常數(shù)k隨升溫速率β的變化而改變.

根據(jù)表5中的數(shù)據(jù)，將不同階段的Eai和lnAi與Eai與lnβi分別進(jìn)行擬合，結(jié)果見圖5、6.

表5 不同升溫速率下的活化能和指前因子

圖5、6中實(shí)驗(yàn)點(diǎn)的線性關(guān)系良好. 由表6可知，擬合曲線的相關(guān)系數(shù)都超過0.97. 表明不同升溫速率下活化能與指前因子的動(dòng)力學(xué)補(bǔ)償效應(yīng)可表示為lnAi=aEai+b，升溫速率對(duì)活化能的影響可表示為Eai=△Ealnβi+Ea0. 由擬合曲線的斜率和截距求得相應(yīng)參數(shù)a、b、△Ea和Ea0，結(jié)果見表6. 將表6中相應(yīng)動(dòng)力學(xué)參數(shù)代入式(8)，得到不同溫度區(qū)間反應(yīng)速率常數(shù)ki的表達(dá)式，見表7.

將表7中Arrhenius方程修正式與等溫下的Arrhenius方程(ki(T)=Aiexp(-Eai/RT))對(duì)比可知，在非等溫條件下，Arrhenius方程修正式中反應(yīng)速率常數(shù)ki是溫度T與升溫速率β的二元函數(shù). 因此，在非等溫條件下將升溫速率β引入到Arrhenius方程中是非常有必要的.

2.4不同升溫速率下的外推性

采用10 ℃/min、15 ℃/min、40 ℃/min和45 ℃/min升溫速率下的煤粉燃燒過程數(shù)據(jù)，來驗(yàn)證Arrhenius方程修正式lnki=-Ea0/RT+△Ea(a-1/RT)lnβi+ln(A0)的外推性. 不同升溫速率下的特征溫度點(diǎn)見表8.

圖5 活化能和指前因子之間的動(dòng)力學(xué)補(bǔ)償效應(yīng)

Fig.5 Kinetic compensation effect between activation energy and pre-exponential factor

圖6 升溫速率與活化能之間的關(guān)系

Fig.6 Relationship between heating rate and activation energy

表6 不同反應(yīng)階段相應(yīng)的動(dòng)力學(xué)參數(shù)和相關(guān)系數(shù)

表7 不同溫度區(qū)間反應(yīng)速率常數(shù)ki的表達(dá)式

表8不同升溫速率下粉煤燃燒過程的特征溫度

Tab.8 Characteristic temperatures of pulverized-coal combustion at different heating rates

升溫速率/(℃/min)Ti/℃T1/℃Tv/℃T2/℃Tb/℃1044650153158562015451510543603642204585225556176632546153657063669330466543580650723354695555886657394047356259667274845476570605680757

從表8中發(fā)現(xiàn)，升溫速率對(duì)著火點(diǎn)溫度Ti、波谷所對(duì)應(yīng)溫度Tv和燃盡點(diǎn)溫度Tb的影響可表示為T=T0+△Tlnβi. 根據(jù)擬合結(jié)果求得相應(yīng)的參數(shù)，見表9.

表9 特征溫度的擬合結(jié)果

由表9可知，擬合曲線的相關(guān)系數(shù)均超過0.97. 表明關(guān)系式T=T0+△Tlnβi，能夠用來描述升溫速率對(duì)分段點(diǎn)溫度的影響，進(jìn)而可知不同升溫速率下溫度區(qū)間的劃分. 結(jié)合表7能直接預(yù)測(cè)不同升溫速率不同溫度下反應(yīng)速率常數(shù)k的值.

用表4中的機(jī)理函數(shù)根據(jù)式(4)對(duì)升溫速率10 ℃/min、15 ℃/min、40 ℃/min和45 ℃/min下煤粉燃燒過程的實(shí)驗(yàn)點(diǎn)進(jìn)行擬合，進(jìn)而求得相應(yīng)的活化能Eai和指前因子Ai. 隨后，分別從兩個(gè)溫度區(qū)間各選取一個(gè)溫度點(diǎn)500 ℃和610 ℃來對(duì)比反應(yīng)速率常數(shù)ki的實(shí)驗(yàn)值和預(yù)測(cè)值. 然后，應(yīng)用不同升溫速率下的活化能和指前因子，根據(jù)Arrhenius方程(ki(T)=Aiexp(-Eai/RT))計(jì)算出反應(yīng)速率常數(shù)k在溫度500 ℃和610 ℃的實(shí)驗(yàn)值見表10.

表10 不同升溫速率下的反應(yīng)速率常數(shù)值ki/s-1

根據(jù)表7中Arrhenius方程表達(dá)式，直接計(jì)算出相應(yīng)溫度(500 ℃和610 ℃)下反應(yīng)速率常數(shù)k與升溫速率β的預(yù)測(cè)曲線. 圖7為反應(yīng)速率常數(shù)實(shí)驗(yàn)值和預(yù)測(cè)曲線的對(duì)比圖. 其中，圖例點(diǎn)代表實(shí)驗(yàn)值，實(shí)線代表Arrhenius方程修正式的預(yù)測(cè)曲線.

圖7 反應(yīng)速率常數(shù)實(shí)驗(yàn)值和預(yù)測(cè)值的對(duì)比

Fig.7 Comparison results between experimental data and predicted curves of reaction rate constants

從圖7可看出，反應(yīng)速率常數(shù)k值隨升溫速率β的增大而減小，且反應(yīng)速率常數(shù)實(shí)驗(yàn)值與預(yù)測(cè)曲線具有良好的一致性. 表明Arrhenius方程修正式對(duì)其他升溫速率下具有良好的外推性. 因此，只要有3組不同升溫速率下煤粉燃燒過程的實(shí)驗(yàn)數(shù)據(jù)就能求得不同反應(yīng)區(qū)間的Arrhenius方程修正式，進(jìn)而求出不同溫度不同升溫速率下的反應(yīng)速率常數(shù).

3 結(jié) 論

在非等溫?zé)岱治鰟?dòng)力學(xué)過程中，本文給出一個(gè)Arrhenius方程的修正式來描述升溫速率β對(duì)煤粉燃燒過程反應(yīng)速率常數(shù)k的影響.

1)隨著升溫速率增加，煤粉燃燒過程表現(xiàn)出明顯的熱滯后現(xiàn)象. 根據(jù)煤粉燃燒過程的特點(diǎn)，以反應(yīng)速率曲線達(dá)到波谷時(shí)對(duì)應(yīng)的溫度點(diǎn)，將著火點(diǎn)到燃盡點(diǎn)的燃燒過程從劃分為兩個(gè)階段. 隨后，根據(jù)兩個(gè)階段的特點(diǎn)，分別采用界面化學(xué)反應(yīng)模型和內(nèi)擴(kuò)散模型來描述相應(yīng)的動(dòng)力學(xué)過程. 結(jié)果表明，實(shí)驗(yàn)數(shù)據(jù)與模型擬合的線性關(guān)系良好，所有相關(guān)系數(shù)均超過0.99.

2)為描述升溫速率對(duì)煤粉燃燒過程動(dòng)力學(xué)的影響，本文給出一個(gè)Arrhenius方程的修正式lnki=-Ea0/RT+△Ea(a-1/RT)lnβi+lnA0. 其中，動(dòng)力學(xué)參數(shù)a、△Ea、Ea0和A0可由3組不同升溫速率下的實(shí)驗(yàn)值求得. 研究表明修正的Arrhenius方程不僅能很好地描述升溫速率對(duì)反應(yīng)速率常數(shù)k的影響，而且還具有一定的外推性.

參考文獻(xiàn)

[1] WU Keng, DING Rucai, YANG Sen, et al. Research on unconsumed fine coke and pulverized coal of BF dust under different PCI rates in BF at Capital Steel Co. [J]. ISIJ International, 2010, 50(3): 390-395. DOI:10.2355/isijinternational.50.390.

[3] LEMAIRE R, BRUHIER C, MENAGE D, et al. Study of the high heating rate devolatilization of a pulverized bituminous coal under oxygen-containing atmospheres [J]. Journal of Analytical & Applied Pyrolysis, 2015, 114:22-31. DOI: 10.1016/j.jaap.2015.04.008.

[4] YU Juan, ZHANG Mingchuan. A simple method for predicting the rate constant of pulverized-coal pyrolysis at higher heating rate [J]. Energy & Fuels, 2003, 17: 1085-1090. DOI: 10.1021/ef020281r.

[5] SUN Xiaoyan, YIN Shumei, WANG Hongyuan, et al. Effect of the addition of cornstalk to coal powder/coal tar combustion [J]. Journal of Thermal Analysis and Calorimetry, 2012, 109(2): 817-823. DOI: 10.1007/s10973-011-1769-1.

[6] 許桂英, WATKINSON A P, ELLIS N, 等. 亞麻纖維熱解動(dòng)力學(xué)的“model free”法和Coats-Redfern模型擬合法研究[J]. 化工學(xué)報(bào), 2010, 09: 2480-2487.

XU Guiying, WATKINSON A P, ELLIS N, et al. Kinetic study on thermal decomposition of flax fibers with model-free and Coats-Redfern model fitting kinetic approaches [J]. CIESC Journal, 2010 (9): 2480-2487.

[7] MIURA K C. A new and simple method to estimatef(E) andk0(E) in the Distributed Activation Energy Model from three sets of experimental data [J]. Energy & Fuels, 1995, 9: 302-307.

[8] KOCABAS Z ?, OKYAY A F, YüRüM ? Y. Combustion characteristics of Turkish hazelnut shell biomass, lignite coal, and their respective blends via thermogravimetric analysis [J]. Journal of Thermal Analysis and Calorimetry, 2015, 119: 1723-1729. DOI: 10.1007/s10973-014-4348-4.

[9] K?K M V. An investigation into the combustion curves of lignites [J]. Journal of Thermal Analysis and Calorimetry, 2001, 64: 1319-1323. DOI: 10.1023/A:1011586105543.

[10]KOK M V. Simultaneous thermogravimetry-calorimetry study on the combustion of coal samples: Effect of heating rate [J]. Energy Conversion and Management, 2012, 53(1): 40-44. DOI: 10.1016/j.enconman.2011.08.005.

[11]XU Zhijie, SUN Xin, KHALEEL M A. A generalized kinetic model for heterogeneous gas-solid reactions [J]. Journal of Chemical Physics, 2012, 137(7): 8653-8660. DOI: 10.1063/1.4740242.

[12]VLAEV L T, MARKOVSKA I G, LYUBCHEV L A. Non-isothermal kinetics of pyrolysis of rice husk [J]. Thermochimica Acta, 2003, 406: 1-7. DOI: 10.1016/S0040-6031(03)00222-3.

[13]WANG Guangwei, ZHANG Jianliang, SHAO Jiugang, et al. Characterisation and model fitting kinetic analysis of coal/biomass co-combustion [J]. Thermochimica Acta, 2014, 591: 68-74. DOI: 10.1016/j.tca.2014.07.019.

[14]劉旭光, 李保慶. 分布活化能模型的理論分析及其在半焦氣化和模擬蒸餾體系中的應(yīng)用[J]. 燃料化學(xué)學(xué)報(bào), 2001, 29(1): 54-59.

LIU Xuguang, LI Baoqing. Theoretical elucidation of distributed activation energy model and its applications in char gasification and simulation distillation systems [J]. Journal of Fuel Chemistry and Technology, 2001, 29(1): 54-59.

[15]巢昌耀, 吳鏗, 杜瑞嶺, 等. 煤粉與半焦的混合燃燒特性及動(dòng)力學(xué)分析[J]. 工程科學(xué)學(xué)報(bào), 2016, 38(11): 1532-1538. DOI:10.13374/j.issn2095-9389.2016.11.004.

CHAO Changyao, WU Keng, DU Ruiling, et al. Combustion characteristics and kinetic analysis of pulverized coal and semicoke[J]. Chinese Journal of Engineering, 2016, 38(11): 1532-1538. DOI:10.13374/j.issn2095-9389.2016.11.004.

[16]MAITRA S, BANDYOPADHYAY N, DAS S, et al. Non-Isothermal decomposition kinetics of alkaline earth metal carbonates [J]. Journal of the American Ceramic Society, 2007, 90(4): 1299-1303. DOI: 10.1111/j.1551-2916.2007.01607.x.

[17]LIU Xiang, CHEN Meiqian, WEI Yuanhang. Kinetics based on two-stage scheme for co-combustion of herbaceous biomass and bituminous coal[J]. Fuel, 2015, 143: 577-585. DOI: 10.1016/j.fuel.2014.11.085.

AmodifiedArrheniusequationforstudyingthecombustionkineticsofpulverized-coalbythermogravimetricanalysis

DU Ruiling, WU Keng, CHAO Changyao, XU Daan, ZHANG Li, ZHANG Bing

(State Key Laboratory of Advanced Metallurgy (University of Science and Technology Beijing) 100083 Beijing, China)

The combustion process of pulverized-coal is widely used in different industries as a supplier of energy. In order to investigate the influence of heating rate on the combustion kinetics of pulverized-coal, the study of combustion pulverized-coal at different heating rates (20 ℃/min, 25 ℃/min, 30 ℃/min, and 35 ℃/min) was conducted by thermogravimetric analysis. The results show that the entire combustion process indicates a thermal hysteresis phenomenon with heating rate increasing. According to the characteristics of the pulverized-coal combustion, the combustion process from ignition temperature to burnout temperature is divided into two regions according to the temperature where the reaction rate curve reaches the valley value. Then, interfacial chemical reaction model and internal diffusion model are successfully applied to describe the two regions of combustion process. The obtained corresponding kinetic parameters show that the kinetic compensation effect between activation energyEaiand pre-exponential factorAifor different heating rates is expressed as lnAi=aEai+band the influence of heating rateβon activation energyEaican be described asEai=△Ealnβi+Ea0. Later, the Arrhenius equation (lnki=-Ea0/RT+lnA0) is modified as lnki=-Ea0/RT+△Ea(a-1/RT)lnβi+lnA0to describe the influence of heating rateβon the rate constantkof pulverized-coal combustion. Finally, the extrapolation reliability of the modified Arrhenius equation is validated by the experimental data at different heating rates(10 ℃/min, 15 ℃/min, 40 ℃/min, and 45 ℃/min). Therefore, the modified Arrhenius equation is not only the most expedient way to depict the combustion kinetic at different heating rates, but also provides extrapolation reliability over a broad range.

pulverized-coal combustion; heating rates; modified Arrhenius equation; kinetic parameters; non-isothermal system

10.11918/j.issn.0367-6234.201612074

TQ530.2

A

0367-6234(2017)11-0081-06

2016-12-14

國(guó)家自然科學(xué)基金(51274026)

杜瑞嶺(1989—)，男，博士研究生

吳 鏗，wukeng@metall.ustb.edu.cn

(編輯苗秀芝)