李方舟，李文英，馮杰

（太原理工大學煤科學與技術教育部和山西省重點實驗室，山西省煤科學與技術省部共建國家重點實驗室培育基地，山西 太原 030024）

?

固體熱載體法褐煤熱解過程中的傳質傳熱特性

李方舟，李文英，馮杰

（太原理工大學煤科學與技術教育部和山西省重點實驗室，山西省煤科學與技術省部共建國家重點實驗室培育基地，山西 太原 030024）

2015-05-12收到初稿，2015-11-06收到修改稿。

聯(lián)系人：李文英。第一作者：李方舟（1990—），男，博士研究生。

Foundation item: supported by the National Natural Science Foundation of China (51276120，U1361202)，the Higher Specialized Research Fund for the Doctoral Program (20121402110016) and the National High Technology Research and Development Program of China(2011AA05A202).

摘要：為揭示在固定床反應器中固體熱載體法快速熱解褐煤工藝過程中的熱、質傳遞機理，建立了固體熱載體法褐煤熱解過程中的傳質傳熱模型。模型包括球型顆粒的一維非穩(wěn)態(tài)導熱方程和基于分布活化能模型的動力學模塊，分別采用有限容積法與Matlab軟件中遺傳算法工具箱對二者進行數(shù)值計算。通過呼倫貝爾褐煤熱重實驗數(shù)據(jù)與溫度測定實驗數(shù)據(jù)分別驗證了預測的動力學參數(shù)及顆粒傳熱模型結果。研究發(fā)現(xiàn)，熱、質變化在固體熱載體法褐煤熱解工藝中呈現(xiàn)復雜的耦合特性。此外，考察了在不同初始溫度、熱載體進料比與煤顆粒半徑條件下，褐煤在熱解過程中顆粒內部溫度場在徑向上隨時間的變化規(guī)律，并分析了產(chǎn)物釋放速率與溫度場的關聯(lián)性，結果表明熱歷程改變是工藝條件對熱解產(chǎn)物分布造成影響的根本原因。

關鍵詞：固體熱載體；褐煤熱解；傳熱；傳質；數(shù)值模擬

引　言

固體熱載體法褐煤熱解技術是指在熱解條件下，將褐煤與具有載熱功能的固體載體混合加熱，得到焦油、煤氣和半焦的技術，屬于褐煤熱解提質加工工藝。采用固體熱載體作為給熱介質，增大了傳熱面積，避免了外熱式熱解時焦油在反應器邊壁分解的問題。此外，褐煤在快速熱解條件下，所含的含氧官能團、烷基側鏈斷裂會形成更多的煤氣和焦油，焦油中的脂肪烴、芳烴和酚類物質可加工獲得化學品和燃料油，熱解半焦熱值得到提高，可作為燃燒或氣化的合格原料[1]，工藝實現(xiàn)了褐煤的分級煉制，提升了氫、氧的元素利用率，是綜合利用褐煤資源的有效方式之一。

為了掌握固體熱載體煤熱解過程中工藝條件對產(chǎn)物組成及分布的作用規(guī)律，許多學者從實驗分析與模型計算兩個角度進行了研究。其中通過建立數(shù)學模型來研究固體熱載體法煤熱解工藝引起了越來越多的關注，例如，Liang等[2]針對移動床固體熱載體熱解工藝建立了一維穩(wěn)態(tài)數(shù)值模型，包含氣相-雙固相傳熱理論與多重反應模型，模型能預測氣、煤與熱載體的軸向溫度以及在不同操作條件下產(chǎn)物的變遷，但忽略顆粒內溫度梯度；郭治等[3]用差分法解固體熱載體法熱解傳熱模型，但邊界條件僅考察熱輻射，并用一級動力學模型描述煤熱解揮發(fā)分釋放過程；王洪亮等[4]對固體間的傳熱途徑做了詳盡的分析，建立了熱載體球煤粉熱解過程的傳熱模型，計算采用Comsol軟件，但簡化了能量方程中源項和動力學部分。

盡管通過數(shù)值計算的方法來研究固體熱載體法煤熱解過程方面已有了一定的成果，但由于系統(tǒng)內熱量傳遞不僅方式多樣，且與顆粒內傳質過程相互影響，并進一步影響產(chǎn)品分布，使機理趨于復雜化，以上對該工藝的數(shù)值模擬均對源項部分或熱解動力學進行了較大的簡化，弱化模型的準確性。

近年來，通過數(shù)值模型對大顆粒煤單獨熱解過程中顆粒內的傳熱規(guī)律的研究逐漸完善。其中，Zhao等[5]采用傳熱方程和復雜煤熱解模型（FGDVC），研究了在流化床或固定床中燃燒的大顆粒煤的熱解過程，模型計算將煤熱解的熱效應設為零；劉訓良等[6-7]在Adesanya等[8]與胡國新等[9]基礎上，結合分布活化能模型（DAEM）及能量守恒方程對煤顆粒單獨熱解的傳熱傳質過程建立數(shù)學模型，將反應熱與揮發(fā)分對流散熱納入源項，提升了模型的準確度。因而有必要結合大顆粒煤熱解過程規(guī)律，建立更全面的數(shù)值模型來進一步探究固體熱載體法煤熱解過程的熱、質傳遞機理。

在褐煤固體熱載體固定床熱解工藝中，原料煤的Biot數(shù)（Bi）大于臨界值0.1[10]，研究固體熱載體法熱解的傳熱過程不宜選用集總熱容法，即大顆粒煤內的溫度梯度產(chǎn)生傳熱阻力，會影響毫米級褐煤顆粒的熱解反應與傳質速率，最終改變產(chǎn)物組成與分布。因此，探究固體熱載體與大顆粒煤混合時的質、熱傳遞規(guī)律是揭示該工藝中熱解產(chǎn)品所經(jīng)歷溫度變化的基礎，本文通過分析進料比、溫度、粒徑和停留時間對傳熱過程的擾動，選擇匹配的工藝條件參數(shù)，為反應器設計、運行提供理論依據(jù)。

1　實驗樣品及方法

1.1樣品

實驗所用原料選取呼倫貝爾褐煤，元素分析和工業(yè)分析見表1。將原煤破碎、篩分至0.15～0.18 mm及8 mm左右兩種粒徑的顆粒樣品，均置于110℃的烘箱內干燥8 h，之后，用銼刀打磨大粒徑的煤粒，制成類球體形狀，將處理后的樣品均分別放到密封容器內待用。并選取粒徑在1.8～2 mm的石英砂樣品（SiO2含量達99.5%以上）備用。

表1　呼倫貝爾褐煤的工業(yè)分析與元素分析Table 1　Physical properties analysis of Hulunbuir lignite

1.2實驗裝置

本文建立的固定床實驗裝置如圖1所示，該系統(tǒng)主要由固定床反應器、溫控電爐、測溫元件和尾氣處理裝置組成。固定床反應器材質為石英玻璃，為減小揮發(fā)分釋放過程中的二次反應，反應器的高徑比設計為1:1（?60 mm×60 mm），距反應器下端15 mm處安裝有石英燒結作為物料支撐板，高純氮（99.999%）作為載氣，從反應器底端進氣口吹進，進料口與出氣口位于頂端。熱電偶6、7分別測量反應器外壁和器內物料表面溫度，選用標號WRNK-162的K型熱電偶，規(guī)格?0.5×600×L2400，允差值為±2.5℃，靈敏度能滿足測量要求。

圖1　固定床實驗裝置流程Fig.1　Flowsheet of fixed bed reactor1—N2cylinder; 2—mass flowmeter; 3—heating furnace; 4—temperature controller; 5—quartz reactor; 6,7—thermocouple; 8—off gas collector

1.3實驗方法

1.3.1熱重實驗在HCT-1型熱重分析儀（北京恒久科學儀器廠）上，以30 K·min-1升溫速率將粒徑為0.15～0.18 mm的10.3 mg樣品從室溫程序升溫到1233 K進行熱重實驗。載氣為高純氮（99.999%），流量為100 ml·min-1。

1.3.2測溫實驗熱電偶分別用銅絲固定在單顆煤粒表層和中心，采集所在部位溫度信息。設定加熱爐升溫程序及載氣流量參數(shù)，待系統(tǒng)穩(wěn)定運行后，將石英砂按一定質量比例加入固定床反應器，并將反應器置于加熱爐中，在石英砂達到設定溫度后，將反應器迅速移出加熱爐，同時將連有熱電偶的煤顆粒插入反應器，使煤顆粒進入石英砂熱載體床層，開始記錄時間和溫度，在顯示溫度開始從峰值緩慢下降時，認為進入散熱階段，停止計時。

2　模型建立與數(shù)值分析方法

在固定床反應器中用固體熱載體使褐煤熱解的過程包含復雜的熱量傳遞機理。在固定床反應器中將預熱后的固體熱載體（石英砂）與煤快速混合，熱載體向周圍環(huán)境發(fā)散熱量，煤顆粒表面先受到氣-固間對流與固-固間輻射的傳熱作用，部分直接接觸熱載體的煤會發(fā)生導熱現(xiàn)象。熱流從煤顆粒表層經(jīng)固相向中心傳遞，整個煤顆粒處于熱的非穩(wěn)態(tài)，表現(xiàn)為顆粒徑向上的溫度梯度隨時間而變化。在顆粒逐漸升溫過程中伴隨著煤的分解，熱解的反應熱會擾動顆粒溫度梯度，同時，顆粒內反應速率、孔隙率及物性參數(shù)如熱導率、比熱容等也隨溫度變化而變動。非黏結性煤釋放的揮發(fā)分與蒸發(fā)的水分以氣體形式從孔道內向顆粒表面擴散，經(jīng)二次反應后進入外部載氣，并帶走部分熱量，如圖2所示。

圖2　褐煤與固體熱載體的傳熱途徑Fig.2　Heat transfer methods between lignite (L) and solid heat carriers (S)Q1—heat transfer interparticles L; Q2—heat transfer between L and S;Q3—heat transfer interparticles S; Q4—heat conduction in L; Q5—heat conduction in S; Q6—heat radiation between L and S; Q7—heat radiation in pores of L; Q8—heat convection between L and volatiles; Q9—heat convection between S and volatiles

針對褐煤固體熱載體法熱解工藝提出下列假設來簡化模型：

（1）假設物料與熱載體在極短時間內均勻混合，且忽略固定床內的軸向返混現(xiàn)象，即物料顆粒停留時間相等，受熱過程相同；

（2）褐煤與石英砂均認為是實心球形顆粒，在快速熱解過程中顆粒形狀、體積不變；

（3）假設顆粒內熱量傳遞是一維非穩(wěn)態(tài)過程，傳熱僅在顆粒徑向上發(fā)生；

（4）模型未考慮二次反應的影響。

2.1固體顆粒傳熱模型

褐煤和石英砂的內部由于存在孔隙及孔內的氣氛，熱量傳遞是由熱傳導、對流傳熱和輻射傳熱3種基本方式共同完成。在假設顆粒均為實心顆粒后，僅考察最主要的熱傳導方式，參照固體在球坐標系中的能量守恒方程式分別建立褐煤與石英砂熱載體的導熱方程。

2.1.1固體熱載體顆粒的導熱方程固體熱載體反應前是等溫熱源，顆粒表層熱量會以接觸煤粒傳導、相鄰煤粒輻射和與氣氛對流的方式向外發(fā)散，并從熱載體內部得到補充。僅作為供熱介質無需考慮廣義源項，其一維非穩(wěn)態(tài)導熱方程[4]

初始條件：邊界條件：

2.1.2煤顆粒的導熱方程褐煤受熱分解時，熱量以顆粒內部氣、固相為傳遞介質，主要的熱流向包含外界受熱、熱解各反應吸放熱（含蒸發(fā)熱）和揮發(fā)分逸出的對流熱散失。顆粒內部溫度梯度的能量平衡方程[6，9]

源項S包含了不能歸入非穩(wěn)態(tài)相、對流相及擴散相中的一切其他項，式(6)右邊3項分別為熱解反應熱、蒸發(fā)熱和熱解揮發(fā)分逸散攜帶的顯熱

由于認為煤顆粒內部存在熱、質傳遞阻力，揮發(fā)分釋放速率不再均勻，揮發(fā)分在r處質流率為

初始條件：

當褐煤快速熱解時，揮發(fā)分劇烈釋放，近似認為環(huán)繞顆粒周圍氣氛與煤粒表面溫度相同，忽略對流傳熱作用。

邊界條件：

2.1.3質量守恒方程由于已經(jīng)假設煤顆粒粒徑保持不變，可近似得到其質量守恒方程

2.2傳熱模型參數(shù)

2.2.1褐煤物性參數(shù)的經(jīng)驗公式由于褐煤在受熱過程中會發(fā)生復雜的熱裂解、熱縮聚反應，其物性參數(shù)在過程中不再是恒定不變的，而是隨著熱歷程與反應程度改變，采用經(jīng)驗公式來表示

比熱容[4]

2.2.2傳熱系數(shù)的經(jīng)驗公式

（1）接觸傳熱系數(shù)的經(jīng)驗公式

Sun和Chen[11]、Natarajan等[12]以及Watson 等[13]都對顆粒間接觸傳熱問題進行過研究。Watson 等[13-14]曾推導出適用于完全彈性碰撞的光滑理想球體的經(jīng)驗公式，本文傳熱模型中不考慮煤與固體熱載體間的接觸應力或摩擦力，可采用這種接觸傳熱系數(shù)的經(jīng)驗公式

復合熱導率

楊氏模量

（2）對流傳熱系數(shù)的經(jīng)驗公式

氣固相傳熱系數(shù)已有多種經(jīng)驗公式，比如Wakao等[15]和Ranz[16]以及Rowe等[17]氣固傳熱系數(shù)模型，可以看出模型中Nu=f(Re，Pr)，Re反映了顆粒周圍流場慣性力的影響，Pr則體現(xiàn)了氣體物理性質的影響，考慮到固體熱載體工藝一般包含不同粒徑顆粒混合的特點，床層的平均空隙率也應納入考察，這里采用Gunn[18]提出的經(jīng)驗公式，其偏差較小[5,19]

Nusselt數(shù)

Prandtl數(shù)

Reynolds數(shù)

（3）輻射傳熱系數(shù)的經(jīng)驗公式

②文化負載詞意義需引申：游覽故宮時，導游講解到：“每一座宮都是一個院落，墻高院深，門戶嚴密。”筆者經(jīng)過釋意分析，認為“院深”意在表現(xiàn)宮殿的封閉和安靜，因此并未直接譯出其字面含義deep，而是將其解釋為closed and quiet.

在工程實際中一般將固體材料視為灰體，引入修正參數(shù)系統(tǒng)黑度εδ，則輻射傳熱系數(shù)可表示為

其中，系統(tǒng)黑度[4]為

2.3源項中動力學參數(shù)模型

源項S公式[式（6）和式（7）]中熱解反應速率和對流散熱量都由揮發(fā)分的釋放速率反映。為得到呼倫貝爾褐煤熱解揮發(fā)分釋放的活化能等動力學參數(shù)，在模型中引入對煤種和升溫速率有寬廣適應性的分布活化能模型（DAEM）[20]，該模型通過式（24）來計算煤熱解到某一時刻揮發(fā)分的釋放量。

為了方便計算及以熱重實驗數(shù)據(jù)驗證，將DAEM方程記做非等溫條件下熱解的DAEM形式

其中，β為褐煤熱解的升溫速率。

本文所研究的固體顆粒傳熱模型包含了球形顆粒的一維非穩(wěn)態(tài)導熱方程和基于分布活化能模型的動力學模塊，其數(shù)值求解方法包括以下幾個步驟：①采用內節(jié)點法將空間區(qū)域離散化；②采用有限容積法（FVM）來離散導熱方程，離散過程為全隱式格式[21]；③采用三對角矩陣法（TDMA）對離散后的方程進行迭代求解；④DAEM公式（25）的內層積分參考文獻[22]的方法簡化

使其適用于離散方程，外層積分用Romberg算法[23]求得高精確度的解；⑤DAEM中的動力學參數(shù)采用Matlab軟件中的遺傳算法工具箱(GAOT)求得。

3　實驗結果與討論

3.1動力學參數(shù)求解與驗證

為了計算過程的簡便，假設熱解反應的DAEM中動力學參數(shù)均為恒定值。對DAEM公式（25）中的動力學參數(shù)E0、k0與σ采用Matlab軟件中的遺傳算法工具箱（GAOT）來求得最優(yōu)解[24]。其中在給定的動力學參數(shù)范圍內（150＜E0＜300；10＜σ＜70）[25]，目標函數(shù)選取各溫度下數(shù)值解與實驗數(shù)據(jù)的偏差，偏差最小時得E0、k0與σ最優(yōu)解。

目標函數(shù)：

其中，wcal,i、wexp,i分別為不同溫度點處失重的數(shù)值解和實驗值，實驗值源自熱重數(shù)據(jù)，實驗條件見1.3.1，n為選取的計算點數(shù)，程序編制流程見圖3。

圖3　呼倫貝爾褐煤的DAEM動力學參數(shù)在GAOT中程序框圖Fig.3　Computing programme of GAOT determining kinetic parameters of Hulunbuir lignite in DAEM

結合熱重實驗數(shù)據(jù)，將上述程序編制流程在Matlab軟件中運行，計算結果為E0=233.14 kJ·mol-1，k0=5.14×1012s-1，σ=29 kJ·mol-1。采用1.3.1實驗條件下的熱重實驗數(shù)據(jù)對計算結果驗證，對比結果見圖4。

圖4　DAEM預測呼倫貝爾褐煤熱解相對失重與實驗數(shù)據(jù)的對比Fig.4　Comparison of DAEM prediction and experimental data from lignite thermogravimetric analysis

從預測結果與實驗結果的對比發(fā)現(xiàn)，用DAEM來預測褐煤熱解失重時會與實驗值有一定的誤差。在DAEM的假設中認為分布活化能E近似為一個正態(tài)分布的函數(shù)，并將頻率因子k0看作定值，在GAOT的編程中接受這一假設。但有研究表明頻率因子并非常數(shù)，而是會隨著活化能的增大而增大，二者之間存在“補償效應”[20,26]。當采用恒定的頻率因子進行計算時會產(chǎn)生誤差，且誤差會隨著溫度升高而增大。但由于本文考察褐煤中低溫快速熱解的熱歷程，因此忽略上述原因造成的誤差影響。

3.2顆粒傳熱模型的驗證

為驗證傳熱模型的正確性，在固定床實驗室裝置上分別測定了熱解過程中煤顆粒表面、中心及石英砂熱載體表面的溫度，測溫實驗方法見1.3.2，實驗條件設定與表2中工況3的計算參數(shù)相同，與傳熱模型預測結果對比見圖5。

圖5　顆粒傳熱模型預測結果與測溫實驗數(shù)據(jù)對比Fig.5　Comparison of particle heat transfer model prediction and temperature measuring experiments

由圖5可見，在相同參數(shù)條件下，顆粒傳熱模型的計算結果與實驗結果吻合較好，表明建立的模型適用，二者稍有偏差的原因可能是煤顆粒與固體熱載體在固定床中接觸方式不同，煤顆粒并非均勻受熱的理想狀態(tài)，相互之間的熱量傳遞受到影響。

3.3煤顆粒內部熱量傳遞的計算結果

為了探究不同初始溫度、熱載體進料比與煤顆粒粒徑等工藝條件下，褐煤在熱解過程中顆粒內部溫度場和產(chǎn)物釋放速率隨時間在徑向上變化規(guī)律，表2列舉了4種工況下的條件參數(shù)。

表2　工況參數(shù)Table 2　Parameters of operation conditions

圖6為對應的4個工況下褐煤顆粒內部溫度場隨時間在徑向上變化的情況。均以工況1[圖6(a)]作為基準與其他工況條件進行比較。

工況2是將固體熱載體的進料溫度提高了200 K，從圖6(b)中可見，褐煤顆粒熱解最終的平衡溫度提升了近220 K，且煤粒達到表面與中心溫差低于3 K的時間縮短了約72 s，外層到達650 K僅用了約19 s，可歸結為更高初始溫度的熱載體與煤接觸，傳熱動力正比于溫差值，使得顆粒表層升溫速率加快，同時也提升顆粒內層熱量傳遞的劇烈程度，煤顆粒內層直到中心在相對較短時間達到均一溫度，熱解反應快而徹底。煤顆粒表層和中心升溫速率相差較大，在到達平衡溫度前的某一時刻煤顆粒徑向上傳熱滯后，產(chǎn)生間隔溫差很大的溫度梯度，在20 s時顆粒外層與中心最大溫差約250 K，之后差距會逐漸減小，煤顆粒自身的傳熱阻力是熱解快慢的速控步驟。

工況3是將褐煤粒徑縮小一半，從圖6(c)中可見在某一時刻，煤顆粒徑向上升溫速率差距較小，即顆粒內部各層溫度梯度不顯著，煤粒表面與中心的溫差達到小于3 K的時間縮短了約85 s。因為小體積煤顆粒在與固體熱載體堆積過程中，會有更多的表面接觸到熱源，固體間熱傳導與輻射作用增強；且熱量傳遞到顆粒中心的距離縮短，越靠近顆粒內層其升溫速率的提升效果越顯著。此外，減小顆粒體積只能促進煤近表面部分的升溫均勻程度，但對煤粒內部升溫速率以及最終平衡溫度影響小。

圖6　不同工況下褐煤顆粒內部溫度場變化情況Fig.6　Changing temperature field in inner particles of lignite under different conditions

工況4是將固體熱載體與褐煤的進料質量比增大了一倍，在圖6(d)中反映出煤顆粒熱解最終平衡溫度提升60 K以上，且煤粒表面與中心的溫差達到小于3 K的時間縮短了約45 s，煤顆粒表層初始升溫速率增大，可以反映出煤顆粒在剛接觸熱載體時受熱均勻，近表面層的熱流密度大，這是因為增加固體熱載體量加大了有效傳熱面積；顆粒表層與中心的溫度差隨時間呈先增大后減小的趨勢，但徑向上溫差減小的幅度不大，顆粒內仍會產(chǎn)生明顯的溫度梯度。由于煤顆粒表面與顆粒內的傳熱方式不同，所以顆粒內外呈現(xiàn)出不同的升溫曲線，由Fourier定律可知，煤顆粒表面主要受到熱傳導與熱輻射作用，傳熱系數(shù)較大，有效傳熱面積對熱通量影響大，而煤的熱導率較小，顆粒內熱通量主要受限于傳熱阻力，固體熱載體量增加對煤表層溫度場擾動劇烈，但對顆粒內影響較弱；此外，煤熱解平衡溫度上升是由于系統(tǒng)內初始能量總量增加，更多的熱量可供給煤的升溫與熱解過程。

3.4煤顆粒內部揮發(fā)分逸出的計算結果

與溫度場類似，顆粒較大的褐煤內部各處的熱解反應速率由于受熱程度不同步，在顆粒內距表層不同位點處也會出現(xiàn)揮發(fā)分瞬時釋放量的差異。假設不考慮氣體、焦油在逸出過程中的二次反應影響，那么揮發(fā)分的釋放速率可認為是溫度與升溫速率的函數(shù)，利用溫度場的規(guī)律就能進一步研究煤顆粒局部揮發(fā)分逸出特性。對工況1，褐煤徑向上各處揮發(fā)分隨時間的釋放速率如圖7所示。

由圖7可見，在熱解過程中，褐煤顆粒局部揮發(fā)分釋放速率隨時間變化呈拋物線形，沿煤顆粒半徑方向的揮發(fā)分生成存在“滯后效應”，煤顆粒表面達到最大釋放速率的時間比中心提前了約95 s，且越接近顆粒內部揮發(fā)分最大釋放速率值越高，時間上越滯后，這與文獻[6]結果類似。與圖6(a)對比，褐煤顆粒內局部揮發(fā)分最大釋放速率均出現(xiàn)在對應徑向上的升溫結束階段，即最終的平衡溫度區(qū)，可見最大熱解反應速率的產(chǎn)生與局部溫度和升溫速率有內在關聯(lián)。根據(jù)DAEM中活化能分布假設，升溫速率會影響煤熱解的基元反應途徑，造成各反應的活化能呈現(xiàn)某種連續(xù)分布的函數(shù)形式，使熱解反應在不同的溫度段具有不同的活化能；DAEM的建立基于Arrhenius方程，反應速率受活化能與溫度影響。褐煤在工況1下熱解30 s后，顆粒中心升溫速率超過表面，且升溫階段比表面長，造成顆粒中心揮發(fā)分最大瞬時釋放速率大于表面。

圖7　工況1中褐煤顆粒局部揮發(fā)分瞬時釋放速率Fig.7　Local transient volatile yield rate of a lignite particle under condition 1

在計算過程中，揮發(fā)分釋放速率是導熱方程源項中的一個動力學參數(shù)，反映了褐煤熱解反應和褐煤與氣相對流散熱程度，但同時它的計算又依賴于導熱方程提供的顆粒內局部的溫度數(shù)據(jù)，體現(xiàn)了質、熱傳遞之間存在復雜的耦合特性。

從計算結果分析，褐煤不斷受到固體熱載體傳遞的熱量，會發(fā)生由表及里推進的熱解過程，生成的揮發(fā)分以氣體的形式從褐煤顆粒的孔道中逸出，釋放速率受熱解反應和傳質的劇烈程度影響。固體熱載體向褐煤傳熱越多，熱解反應速率越快，生成大量揮發(fā)分，造成顆粒內外的氣相濃度差，它與冷熱氣流的對流作用一起構成了揮發(fā)分逸出的驅動力，熱量傳遞主導著這一過程。煤熱解揮發(fā)分的釋放速率規(guī)律與煤顆粒熱解過程中的熱歷程之間存在關聯(lián)性。

4　結　論

用數(shù)值模擬結合實驗數(shù)據(jù)分析方法，對固體熱載體法熱解褐煤這一特定工藝過程中的熱、質傳遞機理進行了研究，詳盡考察了煤與固體熱載體之間的熱量傳遞方式以及褐煤熱解過程中顆粒內部熱、質傳遞間的交互作用；模擬獲得了4個不同工況下，煤顆粒徑向上溫度場的變化規(guī)律。此外，還探究了揮發(fā)分釋放速率在徑向上隨時間的變化規(guī)律，結論如下。

（1）固體熱載體法褐煤熱解工藝條件的變化均會影響煤顆粒內部的溫度場，其中初始溫度升高、熱載體進料比加大及煤顆粒粒徑減小能在不同程度上縮短顆粒整體達到均一、平衡溫度的時間，提高熱解效率。

（2）煤顆粒內局部的揮發(fā)分釋放速率隨時間表現(xiàn)出類拋物線的變化規(guī)律，不同徑向上的位點到達速率峰值的時間不同，越接近顆粒中心滯后效應越明顯，且在顆粒中心位置揮發(fā)分釋放速率峰值最高。

（3）在同一工況下，褐煤顆粒內部的溫度場與揮發(fā)分釋放速率之間存在有關聯(lián)性的變化規(guī)律。熱歷程主導著煤顆粒內部熱解反應進度與方式，是工藝條件對熱解產(chǎn)物造成影響的根本原因。

對固體熱載體法熱解褐煤過程，揮發(fā)分在大顆粒煤的多孔基質內及固體顆粒間存在對流與擴散的傳質現(xiàn)象，并伴隨較劇烈的聚合、裂解等二次反應。因此，后續(xù)的研究會進一步考察褐煤熱解揮發(fā)分在顆粒內的傳質機理及二次反應作用，探究熱、質傳遞歷程中的交互影響，并得出不同工藝條件下的褐煤熱解最終產(chǎn)物組成及分布。

符號說明

A——顆粒表面積，m2

c——顆粒比熱容，J·kg-1·K-1

cp,g——流體比定壓熱容，J·kg-1·K-1

cp,v——揮發(fā)分比定壓熱容，J·kg-1·K-1

d——褐煤顆粒粒徑，m

E——顆粒楊氏模量，Pa

E0——平均分布活化能，kJ·mol-1

e——床層平均空隙率

Fn——法向應力

ΔHevapor——水分的蒸發(fā)熱，kJ·kg-1

ΔHreaction——褐煤熱解的反應熱，kJ·kg-1

R——氣體常數(shù)，8.314 J·(mol·K)-1

Rc,Rp,r——分別為煤顆粒半徑、固體熱載體顆粒半徑和

顆粒徑向上任意位置到中心的距離，m

T,T0,TR——分別為顆粒任意位置溫度、顆粒初始溫度和

顆粒表面溫度，K

Tg——流體溫度，K

t——時間，s

V,V*——分別為煤熱解揮發(fā)分析出量、平衡析出量，ml

VM——工業(yè)分析中水分(ad)，ml

v——顆粒泊松比

|vg?vs|——氣固相相對速率，m·s-1

wC,wH——分別為干燥基下煤中C、H元素的質量分數(shù)，%

β ——升溫速率，K·s-1

ε ——顆粒黑度

λ,λg——分別為顆粒熱導率、流體熱導率，W·(m·K)-1

μg——流體動力黏度，kg·(m·s)-1

ρ,ρ0,ρg,ρw——分別為顆粒密度、煤顆粒初始密度、流體密

度和水密度，kg·m-3

σ0——黑體輻射（斯蒂芬‐玻爾茲曼）常數(shù)，5.67×

10-8W·m-2·K-4

σ ——活化能分布的標準偏差，kJ·mol-1

φ ——顆粒間角系數(shù)

下角標

p ——固體熱載體

t,t+Δt ——分別為某時間節(jié)點、與時間節(jié)點t相差Δt的時間節(jié)點

References

[1] 李文英，鄧靖，喻長連. 褐煤固體熱載體熱解提質工藝進展[J]. 煤化工，2012，40(1): 1-5. LI W Y，DENG J，YU C L. Development of lignite pyrolysis with solid heat carrier [J]. Coal Chem. Ind.，2012，40(1): 1-5.

[2] LIANG P，WANG Z F，BI J C. Simulation of coal pyrolysis by solid heat carrier in a moving-bed pyrolyzer [J]. Fuel，2008，87(4/5): 435-442.

[3] 郭治，杜銘華，杜萬斗. 固體熱載體褐煤熱解過程的數(shù)學模型與模擬計算[J]. 神華科技，2010，8(2): 69-72. GUO Z，DU M H，DU W D. The mathematical model and simulation of solid heat carrier pyrolysis of lignite [J]. Shenhua Sci. Technol.，2010，8(2): 69-72.

[4] 王洪亮，蒙濤，張華，等. 球型固體熱載體煤粉熱解過程傳熱計算及分析[J]. 潔凈煤技術，2014，20(3): 90-94. WANG H L，MENG T，ZHANG H，et al. Heat transfer calculation and analysis of coal particle pyrolysis with solid ball heat carrier [J]. Clean Coal Technol.，2014，20(3): 90-94.

[5] ZHAO Y X，SERIO M A，SOLOMON P R. A general model for devolatilization of large coal particles [J]. Symp. (Int.) Combust.，1996，26(2): 3145-3151.

[6] 劉訓良，曹歡，王淦，等. 煤顆粒熱解的傳熱傳質分析[J]. 計算物理，2014，31(1): 59-66. LIU X L，CAO H，WANG G，et al. Numerical analysis of heat and mass transfer during pyrolysis of coal particle [J]. Chin. J. Comput. Phys.，2014，31(1): 59-66.

[7] LIU X L，WANG G，PAN G，et al. Numerical analysis of heat transfer and volatile evolution of coal particle [J]. Fuel，2013，106: 667-673.

[8] ADESANYA B A，PHAM H N. Mathematical modelling of devolatilization of large coal particles in a convective environment [J]. Fuel，1995，74(6): 896-902.

[9] 胡國新，田學偉，許偉，等. 大顆粒煤在移動床中的熱解模型[J].上海交通大學學報，2001，35(5): 733-736. HU G X，TIAN X W，XU W，et al. Devolatilization model of large coal particles in moving bed [J]. J. Shanghai Jiaotong Univer.，2001，35(5): 733-736.

[10] CHERN J S，HAYHURST A N. A simple theoretical analysis of the pyrolysis of an isothermal particle of coal [J]. Combust. Flame，2010，157(5): 925-933.

[11] SUN J，CHEN M M. A theoretical analysis of heat transfer due to particle impact [J]. Int. J. Heat Mass Transfer，1988，31(5): 969-975.

[12] NATARAJAN V V R，HUNT M L. Heat transfer in vertical granular flows [J]. Exp. Heat Transfer，1997，10(2): 89-107.

[13] WATSON L V，MCCARTHY J J. Heat conduction in granular materials [J]. AIChE J.，2001，47(5): 1052-1059.

[14] BHARADWAJ R，KETTERHAGEN W R，HANCOCK B C. Discrete element simulation study of a Freeman powder rheometer [J]. Chem. Eng. Sci.，2010，65(21): 5747-5756.

[15] WAKAO N，KAGUEI S，F(xiàn)UNAZKRI T. Effect of fluid dispersion coefficients on particle-to-fluid heat transfer coefficients in packed beds: correlation of Nusselt numbers [J]. Chem. Eng. Sci.，1979，34(3): 325-336.

[16] RANZ W E. Friction and transfer coefficients for single particles and packed beds [J]. Chem. Eng. Prog.，1952，48(5): 247-253.

[17] ROWE P N，CLAXTON K T，LEWIS J B. Heat and mass transfer from a single sphere in an extensive flowing fluid [J]. Trans. Inst. Chem. Eng.，1965，43(1): T14-T31.

[18] GUNN D J. Transfer of heat or mass to particles in fixed and fluidised beds [J]. Heat Mass Transfer，1978，21(4): 467-476.

[19] 郭雪巖，柴輝生，晁東海. 大顆粒流化床傳熱數(shù)值模擬與氣固傳熱模型比較[J].上海理工大學學報，2012，34(1): 81-87. GUO X Y，CHAI H S，CHAO D H. Numerical simulation of large particle fluidized bed and comparison of gas-particle heat transfer models [J]. J. Univ. Shanghai Sci. Technol.，2012，34(1): 81-87.

[20] 楊景標，張彥文，蔡寧生. 煤熱解動力學的單一反應模型和分布活化能模型比較[J]. 熱能動力工程，2010，25(3): 301-305. YANG J B，ZHANG Y W，CAI N S. A comparison of a single reaction model with a distributed activation energy one based on coal pyrolysis kinetics [J]. J. Eng. Therm. Energy Power，2010，25(3): 301-305.

[21] 陶文銓. 數(shù)值傳熱學[M]. 上海: 上海交通大學出版社，2001: 79-103. TAO W Q. Numerical Heat Transfer [M]. Shanghai: Shanghai Jiaotong University Press，2001: 79-103.

[22] CAI J M，YAO F S，YI W M，et al. New temperature integral approximation for nonisothermal kinetics [J]. AIChE J.，2006，52(4): 1554-1557

[23] GAO F. Applications of matlab in mathematical analysis [J]. J. Software，2011，6(7): 1225-1229.

[24] 孫玉鳳，高虹，王通洲. 生物質熱解分布活化能模型的遺傳算法實現(xiàn)[J]. 沈陽理工大學學報，2010，29(3): 63-66. SUN Y F，GAO H，WANG T Z. Application of genetic algorithm in distributed activation energy model of biomass pyrolysis [J]. J. Shenyang Ligong Univ.，2010，29(3): 63-66.

[25] MUSTAFA G，SEMMIN F. A direct search method for determination of DAEM kinetic parameters from nonisothermal TGA data (note) [J]. Appl. Math. Comput.，2002，130: 619-628.

[26] MIURA K，MAKI T. A simple method for estimating f(E) and k0(E) in the distributed activation energy mode [J]. Energy Fuels，1998，12: 864-869.

[27] YI Q，F(xiàn)ENG J，LU B C，et al. Energy evaluation for lignite pyrolysis by solid heat carrier coupled with gasification [J]. Energy Fuels，2013，27(8): 4523-4533.

Characteristics of mass and heat transfer in lignite pyrolysis with solid heat carrier

LI Fangzhou，LI Wenying，F(xiàn)ENG Jie
(Key Laboratory of Coal Science and Technology of Ministry of Education and Shanxi Province，Training Base of State Key Laboratory of Coal Science and Technology Jointly Constructed by Shanxi Province and Ministry of Science and Technology，Taiyuan University of Technology，Taiyuan 030024，Shanxi，China)

Abstract:A comprehensive numerical model coupled two correlative one-dimensional unsteady heat conduction equations of spherical particle with distributed activation energy model has been developed for heat and mass transfer mechanism in lignite pyrolysis with solid heat carrier. The finite volume method and the genetic algorithm optimization toolbox based on Matlab software were employed to calculate the thermal and dynamic parameters，separately，and the reliability of the predictions was further validated by thermogravimetric data of the Hulunbuir lignite and temperature measuring experiment on a laboratory-scale fixed bed reactor，respectively. It was found that the variations in mass and heat transfer during lignite pyrolysis with solid heat carrier showed a complex coupling characteristic. The time-dependent rules of temperature field in radial direction have been obtained by varying operation conditions，such as coal particle radius，the initial temperature and feed amount of solid heat carrier. Besides，the relationship between the releasing rate of pyrolytic products and temperature field revealed that the change of temperature field with heating time in lignite pyrolysis was the primary cause of differentdistribution of the pyrolytic products.

Key words:solid heat carrier; lignite pyrolysis; heat transfer; mass transfer; numerical simulation

基金項目：國家自然科學基金項目（51276120，U1361202）；高等學校博士學科點專項科研基金（20121402110016）；國家高技術研究發(fā)展計劃項目（2011AA05A202）。
Received date: 2015-05-12.

Corresponding author:LI Wenying，ying@tyut.edu.cn

DOI：10.11949/j.issn.0438-1157.20150598

中圖分類號：TQ 021

文獻標志碼：A

文章編號：0438—1157（2016）04—1136—09