許世森，王肖肖，劉剛，李小宇，任永強，譚厚章

（1 中國華能集團清潔能源技術研究院有限公司，煤基清潔能源國家重點實驗室，北京 102209；2 西安交通大學熱流科學與工程教育部重點實驗室，陜西 西安 710049）

煤炭仍是我國目前經(jīng)濟發(fā)展的主要能源支柱，同時也是主要的碳排放來源。我國承諾在2030 年實現(xiàn)“碳達峰”，2060 年實現(xiàn)“碳中和”，為實現(xiàn)這一目標需要對煤炭更加清潔高效利用，同時盡量減少煤炭利用后的碳排放。煤氣化是實現(xiàn)煤的清潔高效利用重要技術之一，也是發(fā)展煤基化學合成（氨、甲醇、乙酸、烯烴等）、液體燃料合成（二甲醚、汽油、柴油）、聯(lián)合循環(huán)發(fā)電（integrated gasification combined cycle，IGCC）等工業(yè)過程的基礎。整體煤氣化聯(lián)合循環(huán)發(fā)電融合了化工和電力行業(yè)優(yōu)點，實現(xiàn)了對煤炭利用過程中低污染物（NO和SO）排放以及低成本CO捕集和封存的協(xié)同，將有望實現(xiàn)煤炭利用過程中CO近零排放。氣流床氣化技術因具有煤種適應性廣、處理量大、氣化效率高等優(yōu)點而得到了廣泛應用。

掌握氣流床氣化爐中煤粉顆粒燃燒和氣化的控制機理，對于設計和優(yōu)化氣化爐運行參數(shù)至關重要。氣化爐中煤粉顆粒轉(zhuǎn)化率主要取決于煤粉顆粒固有氣化速率和停留時間，因為焦炭氣化速率遠小于脫揮發(fā)分速率，脫揮發(fā)分完成后剩余的焦炭顆粒主要通過與CO和HO 的氣化反應消耗。事實上，溫度和顆粒粒徑是影響焦炭顆粒氣化速率的重要因素。氣化爐中溫度范圍分布較廣，意味著不同區(qū)域顆粒氣化過程的控制因素不相同。Kim 等研究表明隨著顆粒尺寸增大，焦炭氣化反應性降低，且高溫時影響更加明顯。主要是因為隨著溫度升高內(nèi)孔擴散主導了焦炭的固有氣化反應速率，而減小煤顆粒尺寸可以極大地降低煤粉顆粒在高溫氣化時氣體擴散阻力。已有研究表明煤粉經(jīng)過超細化后，具有優(yōu)越的理化特性和燃燒特性，其在大型燃煤機組再燃技術還原NO已經(jīng)取得了巨大成功。煤粉粒徑小到一定程度后稱其為超細粉，目前對超細煤粉粒徑還沒有嚴格的定義，針對目前氣流床氣化的工程實踐，本文將粒徑小于40μm煤粉顆粒定義為超細煤粉。此外，煤粉粒徑不但影響溫度場、組分場，而且影響氣化爐壁面固態(tài)和液態(tài)渣層分布。因為氣流床氣化爐一般均為液態(tài)排渣，氣化爐正常運行時壁面形成一層動態(tài)的固態(tài)渣層和液態(tài)渣層。合理控制氣化爐溫度和液、固渣層厚度的動態(tài)平衡，對氣化爐穩(wěn)定、可靠運行至關重要。已有學者報道煤粉粒徑對氣流床氣化性能和碳轉(zhuǎn)化率的影響，但其主要針對粒徑大于50μm 的顆粒，且沒有涉及對氣化壁面渣層流動特性研究。

中國華能集團于2012 年在天津建立了中國首座HNCERI（Huaneng Clean Energy Research Institute）兩段干粉加壓氣化爐用于整體煤氣化聯(lián)合循環(huán)發(fā)電系統(tǒng)（IGCC）示范電站，穩(wěn)定持續(xù)運行至今。氣化爐一段反應室產(chǎn)生的高溫合成氣（1400～1600℃）被二段噴入的煤和水蒸氣進行化學降溫，提升了合成氣產(chǎn)量和冷煤氣效率。比國外先進干粉氣化技術冷煤氣效率可提高2%～3%，比氧耗降低10%～15%。然而，該氣化爐實際運行過程中也存在一些問題。比如氣化爐下端黑水過濾收集的濾餅中殘?zhí)剂扛撸瑲饣癄t在不同負荷下長周期運行時飛灰含碳量可達27.82%～44.9%，這主要由于氣化爐二段碳轉(zhuǎn)化率較低。掌握煤粉顆粒在氣化爐中轉(zhuǎn)化控制機理的空間分布對于進一步提高兩段干粉加壓氣化爐碳轉(zhuǎn)化率至關重要，同時也可以為其他同類型干粉氣流床氣化爐提供參考。由于在承受高溫高壓爐壁上開孔取樣和測溫極易引起合成氣的泄露和爆炸，使用實驗的方法很難對氣化爐內(nèi)部流場以及氣化過程控制因素進行研究，而計算流體動力學（CFD）數(shù)值計算具有成本低、靈活性高的優(yōu)點，已被廣泛用于大型工業(yè)設備內(nèi)部流場結(jié)構(gòu)研究。

本文使用CFD 數(shù)值模擬方法，采用考慮氣體在焦炭顆粒表面擴散的隨機孔模型和Seggiani提出的熔渣模型分別計算氣化爐中焦炭顆粒固有氣化反應速率和壁面液、固渣層分布特性，通過改變煤粉粒徑（20μm、40μm、70μm、100μm、200μm），研究煤粉粒徑對兩段干粉加壓氣化爐顆粒壁面沉積特性、液固渣層分布特性以及碳轉(zhuǎn)化率的影響，探討氣化爐中不同區(qū)域煤粉顆粒轉(zhuǎn)化速率的控制機理，旨在為進一步提高兩段干粉加壓氣化爐碳轉(zhuǎn)化率，實現(xiàn)氣化爐長周期穩(wěn)定運行提供理論依據(jù)。同時，也可為其他干粉氣流床氣化爐進一步提升氣化性能提供參考。

1 數(shù)值計算模型和研究對象

1.1 物理模型

兩段干粉加壓氣化爐由兩個直徑相等的圓柱組成，通過一個縮口連接，其主要結(jié)構(gòu)參數(shù)如圖1所示。氣化爐下半部分圓柱是一段燃燒室，有4個對置安裝的燒嘴；上半部分圓柱是二段還原室，有兩個對置安裝的燒嘴。為了研究氣化爐壁面溶渣沉積特性，對氣化爐一段燃燒室壁面進行了分區(qū)，如圖1所示，燃燒室內(nèi)噴嘴下方直段均分為5層，噴嘴上方直段分13層，上下錐形縮口各均分為3層。

圖1 氣化爐主要尺寸及分層示意圖（單位：mm）

1.2 數(shù)學模型

在實際氣化爐運行過程中，氣化爐中氣固兩相會發(fā)生一些幾乎同時進行的復雜物理化學過程，主要包括：①氣體與煤粉顆粒的湍流混合；②氣相各組分之間的擴散、混合及復雜而劇烈的化學反應；③煤粉顆粒脫揮發(fā)分及揮發(fā)分裂解燃燒過程；④煤顆粒與氣體之間的異相反應過程；⑤對流與輻射換熱過程。本文采用一系列子模型描述以上過程，具體模型如表1 所示。共考慮的13 個氣相反應，如表2所示。下文介紹計算焦炭固有氣化速率的隨機孔模型和熔渣沉積模型。

表1 數(shù)值模擬所采用的子模型

表2 氣相反應反應動力學參數(shù)[14]

1.2.1 焦炭氣化模型

準確計算焦炭氣化反應速率對于準確的模擬氣化爐內(nèi)各個場分布和氣化爐碳轉(zhuǎn)化率至關重要。一般認為焦炭異相反應速率分為三個控制區(qū)域：區(qū)域Ⅰ（化學反應速率控制，低溫）、區(qū)域Ⅱ（孔隙擴散和化學反應速率共同控制，中高溫）和區(qū)域Ⅲ（氣體擴散控制，高溫）。本文使用隨機孔模型，并考慮氣體在孔擴散的影響，使用合適的有效因子η來修正不同區(qū)域的總體反應速率，每種氣化劑的氣化速率方程如式(1)所示，每種氣化劑與焦炭氣化反應動力學參數(shù)采用Kajitani實驗研究數(shù)據(jù)，如表3所示。

表3 焦炭氣化反應動力學參數(shù)[19]

有效因子η的確定如式(2)。

其中，?為Thiele 模量，對于每種氣化劑?表示如式(3)。

對于每種氣化劑的有效擴散率如式(4)。

假設煤的孔隙度為0.5，彎曲度=1/。

1.2.2 渣層流動模型

使用Seggiani提出的熔渣模型計算氣化爐壁面渣層分布和熱損失。圖2(a)為某一層壁面渣單元熱量質(zhì)量傳遞截面示意圖，對每一個渣層控制單元建立質(zhì)量、能量、動量守恒方程。質(zhì)量守恒、能量守恒、動量守恒分別如式(5)～式(7)。

邊界條件為：=0，d/d=0；=，=0。

熔渣流動模型與氣化爐模型的集成以及模型的輸入輸出變量如圖2(b)所示，CFD數(shù)值計算為熔渣模型提供、和，熔渣模型使用MATLAB軟件進行求解計算，將計算得到熔渣表面溫度返回給CFD繼續(xù)求解，CFD計算得到新的、和，傳遞給熔渣模型求解出新的，依此交互，直至熔渣模型兩次求解得到的壁面溫度相差在5%之內(nèi)。熔渣模型同時得到氣化爐壁面固、液渣層厚度分布特性。煤灰黏性隨溫度變化的實驗測量數(shù)據(jù)如圖3 所示，計算中所使用的煤灰渣其他參數(shù)如表4所示。

表4 煤渣物性參數(shù)和計算參數(shù)

圖2 熔渣熱量傳遞和模型變量傳輸示意圖

圖3 煤灰黏度隨溫度變化特性

由于氣化爐模型對稱性，采用四分之一模型進行數(shù)值模擬，網(wǎng)格和邊界設置如圖4所示，網(wǎng)格均采用六面體結(jié)構(gòu)化網(wǎng)格，網(wǎng)格數(shù)量為66 萬。入口邊界條件如表5 所示；表6 和表7 為基礎工況計算所用煤粉入爐基工業(yè)分析、元素分析和粒徑分布，其中煤粉顆粒平均粒徑為35μm。使用Simple算法對壓力和速度進行耦合求解；壓力使用PRESTO 格式；動量、、方程、組分輸運方程、能量方程均使用二階迎風格式；收斂標準除能量和DO方程設置為10，其余均為10。

表5 入口邊界條件

表6 煤工業(yè)分析（質(zhì)量分數(shù)）

圖4 網(wǎng)格及邊界示意圖

一般情況下氣化爐中煤粉顆粒粒徑分布不是均勻分布的，但為了更加清楚地理解氣化過程中煤粉粒徑變化對氣化爐中碳轉(zhuǎn)化率和壁面渣層沉積特性的影響，模擬中保持其他邊界條件不變，將入口煤粉粒徑設置為統(tǒng)一直徑，具體設置工況如表7所示。

表7 基礎工況粒徑分布

2 模擬結(jié)果

2.1 基礎工況

表9為本文所使用數(shù)學模型預測的氣化爐出口主要氣體組分摩爾分數(shù)和碳轉(zhuǎn)化率與工業(yè)運行數(shù)據(jù)和樊強等模擬結(jié)果的對比?？梢钥闯霰疚乃褂脭?shù)學模型對氣化爐出口主要氣體組分摩爾分數(shù)和碳轉(zhuǎn)化率的預測結(jié)果與工業(yè)運行值一致，微量氣體甲烷的預測值與工業(yè)運行值存在一定偏差。樊強等使用簡化的概率密度（PDF）模型考慮湍流和化學反應相互作用，模擬結(jié)果代表了平衡狀態(tài)下氣化爐出口各個組分的含量，可以看出氣化爐出口組分濃度基本達到平衡狀態(tài)。相比樊強等使用的Wen等提出的焦炭反應模型，本文所使用的考慮氣體在焦炭顆粒表面擴散的隨機孔模型計算焦炭氣化速率，并結(jié)合Seggiani提出的灰渣模型預測氣化爐壁面熱損失，更準確地預測出了氣化爐一段和二段的碳轉(zhuǎn)化率。圖5(a)顯示了基礎工況對應的氣化爐中心截面溫度場，可以看出氣化爐一段溫度很高，射流區(qū)和上下流股區(qū)溫度在1900K之上，這可以保證煤粉顆粒在較短停留時間內(nèi)達到高碳轉(zhuǎn)化率，二段燒嘴噴入水蒸氣和煤粉對合成氣進行化學降溫導致二段合成氣溫度持續(xù)下降。從圖5(b)和(c)可以看出氣化爐一段四股高速氧氣射流撞擊后形成兩個大的上下循環(huán)區(qū)，這可以增加顆粒的停留時間，氣化爐二段氣體流動主要以平推流為主。

圖5 基礎工況溫度和速度矢量分布

表8 工況設置

表9 氣化爐出口組分模擬值與工業(yè)值對比

圖6為氣化爐一段壁面氣體溫度、渣層表面溫度和顆粒在壁面沉積率分布。從圖6(a)可以看出，氣化爐上縮口處壁面氣體溫度約為1800K。隨著與縮口距離增大，壁面氣體溫度逐漸降低，至燒嘴平面上部最低，約為1750K，隨后迅速升高，至下端出口處又迅速降低。這是因為燒嘴平面上部上循環(huán)區(qū)的存在，導致上流股區(qū)的高溫氣體沿著壁面回流被壁面逐漸冷卻，而在燒嘴平面以下也存在類似的下循環(huán)區(qū)，且燒嘴平面以下區(qū)域相當于一個高溫煙氣循環(huán)的死區(qū)，因此溫度較燒嘴上部高。壁面渣層表面溫度分布趨勢和氣體溫度分布趨勢相同，但溫度值約低100K，這與Seggni對Shell 氣化爐壁面液渣表面溫度模擬結(jié)果相似。圖6(b)為灰顆粒在壁面沉積率沿著壁面分布，顆粒沉積率在氣化爐底部縮口斜面最大，這是因為大量顆粒向下運動因慣性撞擊到底部斜面導致。顆粒沉積率在燒嘴平面附近最小，燒嘴平面以上區(qū)域分布比較均勻，不同于Shell氣化爐壁面顆粒沉積率的雙峰分布。

圖6 壁面氣體溫度、熱流度和顆粒沉積率分布

圖7為基礎工況壁面液固渣層厚度和壁面熱流密度分布，可以看出液態(tài)渣層厚度沿著氣化爐壁面從1.5mm 逐漸增加到5mm，但在燒嘴平面下部略有減小，之后在底部斜面迅速增加。不同于四燒嘴切圓布置的Shell 氣化爐和單噴嘴旋流的西門子（GSP）氣化爐壁面液態(tài)渣層厚度分布，其特征均是沿著氣化爐壁面從上至下持續(xù)增加，但三種粉煤氣流床氣化爐液渣厚度分布均在1～5mm 之間。造成兩段干粉加壓氣化爐燒嘴平面下部液渣厚度減小的原因在于此區(qū)域壁面附近氣體溫度急劇升高，導致液態(tài)渣層黏度減小，流動速度增大，同時顆粒沉積量在此區(qū)域較小。而在氣化爐底部因為斜面的存在導致渣層流速減小，同時顆粒沉積量在此區(qū)域最大，兩者共同導致液渣厚度在氣化爐底部急劇增大。從圖7(b)可以看出，氣化爐壁面固態(tài)渣層厚度在縮口處約為12mm，沿著壁面逐漸增加，同樣是到燒嘴平面上部達到最大值30mm，之后急劇減小，到底部斜面又急劇增加。這種分布規(guī)律與壁面氣體溫度和顆粒沉積率分布息息相關。GSP氣化爐壁面固態(tài)渣層厚度分布是先減小后增加，厚度在5～40mm 之間。三種氣流床氣化爐壁面液態(tài)渣層分布規(guī)律的不同，反映了三種氣流床氣化爐內(nèi)部氣流組織方式和煤粉顆粒運動規(guī)律存在根本差異。

圖7(c)顯示了從穩(wěn)態(tài)能量守恒方程中計算得到的兩段干粉加壓氣化爐一段熱流密度分布?？梢钥闯鰺崃髅芏仍跉饣癄t縮口處為152kW/m，沿著壁面向下逐漸減小，至燒嘴平面上部最小，這是因為壁面氣體溫度沿著此方向逐漸降低[圖6(a)]，導致溫差降低，同時液態(tài)渣層和固態(tài)渣層厚度逐漸增大[圖7(a)、(b)]，導致熱阻升高。而在燒嘴平面下部直段熱流密度逐漸增加，至底部斜面熱流密度再次降低?？梢钥闯鰵饣癄t壁面熱流密度分布與爐膛溫度和壁面渣層厚度分布息息相關。實際工業(yè)氣化爐運行過程中一段產(chǎn)生5MPa高壓蒸汽量為14～17t/h，計算得到所需熱量為9～10.3MW，而數(shù)值計算結(jié)果表明，氣化爐一段總散熱量9.28MW左右，說明本文所使用的的模型可以較好的預測氣化爐一段壁面散熱。

圖7 基礎工況壁面液固渣層厚度和壁面熱流密度分布

2.2 粒徑對氣化爐壁面渣層的影響

圖8為煤粉粒徑對氣化爐壁面沉積率和液態(tài)渣層表面溫度分布的影響，從圖8(a)可以看，出煤粉粒徑主要影響氣化爐底部斜坡附近壁面顆粒沉積率分布，而對氣化爐一段縮口和直段顆粒沉積量影響較小，表現(xiàn)為隨著煤粉粒徑從20μm 增加到200μm，底部斜面顆粒沉積率最大值從0.12kg/m增大到0.39kg/m。這主要是因為大顆粒質(zhì)量較大，轉(zhuǎn)化速率較慢，在氣流撞擊后受重力影響易向氣化爐下部運動，而向下運動過程中由于慣性較大更易沉積在下端斜坡壁面。從圖8(b)可以看出，隨著煤粉粒徑增大液態(tài)渣層表面溫度呈上升趨勢，其中粒徑小于100μm 時，渣層表面溫度從上至下分布趨勢與基礎工況一致。當粒徑增大到200μm 時，渣層表面溫度比其他粒徑下高出很多，且從氣化爐縮口至底部斜面一直增大。這是因為煤粉粒徑增大，一段碳轉(zhuǎn)化率降低，而氣化爐內(nèi)焦炭的轉(zhuǎn)化主要通過吸熱氣化反應完成，導致了氣化爐總體爐溫升高。

圖8 煤粉粒徑對顆粒沉積率和液態(tài)渣層表面溫度影響

圖9顯示了煤粉粒徑對氣化爐液、固渣層厚度的影響。從圖9(a)可以看出，煤粉粒徑的改變對氣化爐液態(tài)渣層沿壁面分布規(guī)律基本沒有影響，但隨著煤粉粒徑增大，液態(tài)渣層厚度在氣化爐上縮口和直段減小，變化量均很小。這主要是因為隨著煤粉粒徑增大液渣表面溫度升高，導致渣層黏度降低，流速增大。而在底部斜坡和下渣口附近液渣厚度隨著粒徑增大而增加，這是因為此區(qū)域粒徑增大灰渣沉積率增加所致。分析圖9(b)可知，煤粉粒徑增大，氣化爐上縮口和直段壁面固態(tài)渣層厚度降低，而在底部斜坡壁面固態(tài)渣層厚度基本相同，這主要與液態(tài)渣層表面溫度和熱流密度有關。

圖9 粒徑對壁面固液渣層影響

2.3 粒徑對氣化爐氣化性能的影響

溫度是影響焦炭氣化反應速率的重要因素，而氣化爐中溫度分布范圍較廣。本文定義氣化爐中某空間位置單位體積內(nèi)焦炭反應到原來濃度1/所需要時間為焦炭氣化反應的特征時間，大小為氣化爐空間某位置單位體積網(wǎng)格內(nèi)焦炭濃度的1/除以當?shù)氐臍饣磻俾?。氣化爐中氣固混合也是影響焦炭氣化速率的重要因素，氣化爐中宏觀混合主要由氣相湍流脈動控制，根據(jù)湍流能譜理論，積分尺度渦團表征了流場中完成大部分能量和動量輸運渦團的特性，同湍流混合密切相關，定義積分尺度渦團時間尺度τ=k/ε。HNCERI 氣化爐一段射流區(qū)高溫區(qū)可達2300K以上，而二段合成氣經(jīng)過化學降溫后溫度只有不到1400K。為了清晰地認知兩段干粉加壓氣化爐內(nèi)焦炭轉(zhuǎn)化的控制機理，圖10 顯示了基礎工況下氣化爐一段和二段積分渦時間尺度和氣化反應R15特征時間尺度分布。因為反應R16特征時間分布和R15 類似因此不再列出。從圖10(a)和(b)可以看出，一段射流區(qū)積分渦時間尺度數(shù)量級為10，撞擊區(qū)和靠近燒嘴平面的上下流股區(qū)積分時間尺度數(shù)量級為10，這與此區(qū)域R15特征時間尺度在同一數(shù)量級，說明此區(qū)域內(nèi)宏觀尺度的湍流脈動對焦炭顆粒異相氣化反應有強烈的影響。而氣化爐一段燃燒室內(nèi)其他區(qū)域R15特征時間尺度在10～10。說明氣化爐一段焦炭氣化區(qū)主要為射流區(qū)和靠近燒嘴平面上下流股區(qū)。R15特征時間尺度分布也說明氣化爐一段碳轉(zhuǎn)化率主要受到固有氣化速率的影響，因此可以通過提升顆粒固有氣化反應速率來提升一段碳轉(zhuǎn)化率。圖10(c)和(d)可以看出氣化爐二段射流區(qū)積分時間尺度數(shù)量級為10～10，而焦炭氣化反應R15的特征時間尺度數(shù)量級為10，說明二段宏觀尺度的湍流脈動對焦炭異相反應影響很小，此時可以依靠提升顆粒停留時間來提高二段碳轉(zhuǎn)化率。

圖10 氣化爐一段和二段氣相積分渦時間尺度和氣化反應R15特征時間分布

從圖11(a)可以看出當煤粉粒徑為20μm 和40μm時，氣化爐一段碳轉(zhuǎn)化率分別可達99.68%和99.64%，說明煤粉粒徑達到超細粉有利于氣化爐一段焦炭顆粒轉(zhuǎn)化；當煤粉粒徑增加到70μm 時，一段碳轉(zhuǎn)化率降低至99.26%；繼續(xù)增大顆粒粒徑，一段碳轉(zhuǎn)化率明顯降低，至粒徑為200μm 時，一段碳轉(zhuǎn)化率僅有95.06%。這是因為氣化爐內(nèi)焦炭主要通過反應R15和R16消耗，氣化爐一段溫度很高，其中絕大部分區(qū)域溫度大于1900K。此時孔耗散對于焦炭氣化反應速率的影響變得比化學反應速率重要。而焦炭氣化反應模型中有效因子η與顆粒粒徑成反比，因此隨著粒徑增大氣化反應速率減小。此外隨著顆粒粒徑增大，顆?？偙砻娣e也減小，同時從氣相到顆粒相的輻射傳熱也會減小，這都會降低氣化反應的速率。再者，從圖11(b)可以看出，隨著顆粒粒徑從20μm增大到200μm，一段顆粒平均停留時間從2.5s降低到了1.5s，說明一段煤顆粒停留時間也是焦炭轉(zhuǎn)化的影響因素。以上原因共同導致了隨著粒徑增加一段碳轉(zhuǎn)化率降低。

圖11 粒徑碳轉(zhuǎn)化率和停留時間的影響

從圖11(a)可以看出，氣化爐二段碳轉(zhuǎn)化率在粒徑為20μm時為69.03%，隨著粒徑增加到70μm，二段碳轉(zhuǎn)化率緩慢增加到了71.56%；當顆粒粒徑增大到100μm時，二段碳轉(zhuǎn)化率提升到了78.75%；當顆粒粒徑進一步增大到200μm 時，二段碳轉(zhuǎn)化率增加到了89%。這與隨著粒徑增大，一段碳轉(zhuǎn)化率降低，導致進入二段的氣體溫度更高，氣化劑CO和HO 濃度也增加有一定關系，但根據(jù)上文分析可知，氣化爐二段顆粒轉(zhuǎn)化率主要受到停留時間的控制，從圖11(b)也可以看出，顆粒在二段的停留時間和碳轉(zhuǎn)化率分布趨勢非常一致。與此同時，煤顆粒粒徑增加的同時會增加氣固相之間的滑移速度，這有利于增加氣固相之間的傳質(zhì)和傳熱。在氣化爐二段較長停留時間作用下，其有利于大顆粒氣化反應的進一步發(fā)生。對于新設計的兩段氣化爐可以考慮將二段噴嘴布置為切圓形式，這樣可以進一步增加煤粉顆粒在氣化爐二段中的停留時間。

3 結(jié)論

本文采用考慮了焦炭表面氣體組分擴散效應的隨機孔模型和Seggiani 提出的熔渣模型分別計算氣化爐內(nèi)焦炭氣化速率和壁面熔渣分布特性，研究了煤粉顆粒粒徑對HNCERI氣化爐壁面固、液渣層厚度分布和碳轉(zhuǎn)化率的影響，探討了氣化爐一段和二段煤粉顆粒轉(zhuǎn)化的控制因素。通過預測顆粒粒徑對氣化爐液固渣層厚度和碳轉(zhuǎn)化率影響，得到以下結(jié)論。

（1）本文所構(gòu)建的的數(shù)學模型對氣化爐出口碳轉(zhuǎn)化率、主要氣體組分的預測值與工業(yè)運行值一致，同時也較好地預測了氣化爐一段水冷壁散熱量。

（2）隨著煤粉粒徑增大，氣化爐底部斜坡顆粒沉積量增加。但液態(tài)渣層厚度分布基本不受粒徑變化的影響，固態(tài)渣層厚度分布受到粒徑變化影響顯著。

（3）氣化爐一段宏觀混合湍流脈動對焦炭氣化速率有很大影響，焦炭固有氣化速率和顆粒停留時間共同控制一段煤顆粒碳轉(zhuǎn)化率。模擬結(jié)果表明，隨著顆粒粒徑從20μm增加到200μm，一段碳轉(zhuǎn)化率從99.64%降低到了95.06%，即煤粉粒徑減小至超細粉有利于提高一段碳轉(zhuǎn)化率。

（4）氣化爐二段焦炭氣化速率受宏觀混合湍流脈動影響較小，二段焦炭轉(zhuǎn)化率主要受顆粒停留時間控制。模擬結(jié)果表明，隨著顆粒粒徑從20μm增加到200μm，二段碳轉(zhuǎn)化率69.03%增加到了89%，即煤粒徑適當增大有利于提高二段碳轉(zhuǎn)化率。

A—— 壁面面積，m

—— 指數(shù)前因子

c—— 渣層比熱容，J/(kg·K)

—— 有效擴散率，m/s

—— Knudsen擴散系數(shù)，m/s

—— 分子擴散系數(shù)，m/s

—— 焦炭顆粒的直徑，m

—— 活化能，J/kmol

—— 重力加速度，m/s

—— 碳的分子量，kg/kmol

—— CO的分子量，kg/kmol

—— 流出控制單元渣層質(zhì)量流量，kg/s

—— 渣沉積到壁面的質(zhì)量流量，kg/s

—— 焦炭顆粒質(zhì)量，kg

—— 初始狀態(tài)焦炭顆粒質(zhì)量，kg

—— 反應級數(shù)

—— 總壓，Pa

P—— 組分的分壓，Pa

—— 單位時間流出控制單元的渣帶走的熱量，W

—— 灰渣熔化熱，J/kg

—— 煙氣進入渣層表面的熱流密度，W/m

—— 耐火材料與金屬壁面間的熱流密度，W/m

—— 灰渣與耐火材料間的熱流密度，W/m

—— 溫度，K

—— 臨界黏度對應的溫度，K

—— 時間，s

—— 氣化劑的化學計量因子

—— 焦炭顆粒的總碳轉(zhuǎn)化率

Y—— 物質(zhì)的質(zhì)量分數(shù)

—— 渣層總厚度，m

—— 液態(tài)渣層厚度，m

—— 固態(tài)渣層厚度，m

—— Thiele模量

—— 渣黏度，Pa·s

η—— 有效因子

—— 炭顆粒的孔隙率

—— 渣密度，kg/m

—— 氣體密度，kg/m

—— 焦炭顆粒的密度，kg/m

—— 孔隙的彎曲度

—— 孔隙結(jié)構(gòu)參數(shù)