李 偉，王金星，黃江城，王春波

（華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003）

300 MW煤粉／高爐煤氣混燃鍋爐NOx優(yōu)化數(shù)值模擬

李 偉，王金星，黃江城，王春波

（華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003）

利用Fluent軟件對(duì)某鋼廠300 MW煤粉／高爐煤氣混燃鍋爐進(jìn)行數(shù)值模擬，探討了燃燼風(fēng)配比、高爐煤氣配比、過量空氣系數(shù)及運(yùn)行負(fù)荷等因素對(duì)爐膛內(nèi)CO濃度分布及NOx濃度分布的影響規(guī)律。研究表明：燃燼風(fēng)配比增加，燃燒器附近的CO濃度峰值增加且爐膛內(nèi)部的NOx濃度整體降低；高爐煤氣配比增加，爐膛內(nèi)部的CO濃度的峰值下降，同時(shí)NOx濃度呈整體下降趨勢(shì)；隨著過量空氣系數(shù)增加，燃燒器上部的CO濃度均有所下降，而NOx濃度逐漸上升；運(yùn)行工況對(duì)爐膛內(nèi)部CO濃度的分布影響不大，隨著負(fù)荷的降低，峰值略有升高，然而NOx濃度分布在燃燒器上部受負(fù)荷影響較顯著，隨著負(fù)荷下降，NOx濃度逐漸下降。

煤粉；高爐煤氣；混燃鍋爐；NOx；CO；數(shù)值模擬

0 引言

隨著鋼鐵工業(yè)的發(fā)展，生產(chǎn)中的副產(chǎn)煤氣（高爐煤氣、焦?fàn)t煤氣及轉(zhuǎn)爐煤氣等）大量增加，這些煤氣約占鋼鐵企業(yè)總能耗量的30%～40%［1］。其中焦?fàn)t和轉(zhuǎn)爐煤氣熱值較高可在生產(chǎn)或生活中利用，因此高爐煤氣的回收和利用是鋼鐵企業(yè)節(jié)能降耗的重要環(huán)節(jié)。其中高爐煤氣在鋼鐵企業(yè)中產(chǎn)量巨大，每煉1 t生鐵約可以得到高爐煤氣3 500 m3。高爐煤氣中可燃成分主要是CO，但是所占比例很小只有20%左右，不可燃成分體積分?jǐn)?shù)高達(dá)80%左右。所以高爐煤氣熱值非常低，利用率很低而且燃燒不穩(wěn)定［2］。因此，如何有效地利用高爐煤氣是目前面臨的節(jié)能問題。

高爐煤氣和煤粉在爐內(nèi)摻燒是一種有效利用途徑［3］。然而由于高爐煤氣的低熱值，混燃鍋爐的燃燒存在燃燒穩(wěn)定性差、煤粉燃盡困難等問題［4］。目前，國(guó)內(nèi)外學(xué)者針對(duì)煤粉爐內(nèi)燃燒進(jìn)行了模擬［5～10］。對(duì)煤粉／高爐煤氣混燃鍋爐模擬的研究較少，尤其針對(duì)煤粉／高爐煤氣混燃爐膛內(nèi)NOx的分布特性更為罕見［11］。本文利用Fluent軟件對(duì)煤粉／高爐煤氣混燃鍋爐進(jìn)行模擬，以便對(duì)爐膛內(nèi)NOx分布進(jìn)行分析探討，進(jìn)而期望對(duì)今后研究煤粉／高爐煤氣混燃鍋爐NOx優(yōu)化排放提供參考。

1 研究對(duì)象

選取某鋼廠300 MW機(jī)組煤粉／高爐煤氣混燃鍋爐為模擬對(duì)象。該鍋爐為亞臨界參數(shù)、一次中間再熱、自然循環(huán)汽包爐，采用四角切圓燃燒方式，設(shè)計(jì)燃料為煙煤，設(shè)計(jì)燃燼風(fēng)配比為10%，摻燒高爐煤氣比例為20%。燃用煤的實(shí)驗(yàn)分析見表1，高爐煤氣成分分析見表2。

鍋爐為單爐膛，燃燒器為四角布置的擺動(dòng)式燃燒器，切向燃燒，上下擺動(dòng)的最大角度為±30°。燃燒器出口射流中心線和水冷壁中心線的夾角分別為38°和44°，在爐膛中心形成逆時(shí)針旋向的兩個(gè)直徑不同的假想切圓。為了削弱爐膛出煙氣的旋轉(zhuǎn)強(qiáng)度，減小四角燃燒引起爐膛出

2 計(jì)算方法

表1 煤的實(shí)驗(yàn)分析Tab.1 Experimental analysis of coal

表2 高爐煤氣成分分析Tab.2 Analysis of BFG

出口煙溫偏差，燃燼風(fēng)室被設(shè)計(jì)成反切，使其噴嘴出口中心線同主噴嘴中心線成12°的夾角。鍋爐的寬、深為14 048 mm×12 468 mm，高為58 900 mm。鍋爐共布置7層燃燒器噴口，其中上5層為煤粉噴口，下2層為高爐煤氣噴口。鍋爐爐膛構(gòu)架和燃燒器布置見圖1。

模擬軟件采用Fluent，爐膛網(wǎng)格采用分段劃分，從冷灰斗至爐膛出口劃分為7部分。爐膛截面網(wǎng)格采用Paving方法生成非結(jié)構(gòu)四邊形網(wǎng)格，燃燒器體網(wǎng)格用Cooper方法沿著爐膛高度方向鋪展生成六面體網(wǎng)格。Paving方法生成的輻射狀網(wǎng)格線與四角射流的氣流軌跡基本平行，減小了網(wǎng)格線與流線的夾角，降低數(shù)值偽擴(kuò)散［12］。由于燃燒器區(qū)域的溫度梯度較高，因此為保證數(shù)值計(jì)算的精確性，對(duì)該區(qū)域劃分的網(wǎng)格較密。整個(gè)爐膛網(wǎng)格總數(shù)約為120萬。

數(shù)值計(jì)算中湍流流動(dòng)采用雙方程模型，輻射傳熱采用P-1模型［13］，離散相顆粒軌跡采用隨機(jī)跟蹤模型，焦炭燃燒采用動(dòng)力—擴(kuò)散限制模型，揮發(fā)份熱解采用兩步競(jìng)相反應(yīng)模型，氣相湍流燃燒采用混合分?jǐn)?shù)／概率密度函數(shù)方法。

熱力型NOx根據(jù)廣義的Zeldovich機(jī)理計(jì)算，燃料型NOx根據(jù)DeSoete機(jī)理分為揮發(fā)份NOx和焦炭NOx兩部分來計(jì)算，其中揮發(fā)分N占總?cè)剂螻的比例為45%～60%［14］。在計(jì)算中揮發(fā)份N轉(zhuǎn)化為中間產(chǎn)物HCN，焦炭N直接轉(zhuǎn)化為NOx。

3 模擬結(jié)果及分析

為能整體描述爐膛內(nèi)部CO和NOx的分布情況，本文CO濃度及NOx濃度所選用的為某一爐膛高度的截面平均濃度。

3.1 燃燼風(fēng)配比的影響

燃燼風(fēng)是減小NOx排放的重要參數(shù)。選取高爐煤氣摻燒比例為20%，過量空氣系數(shù)為1.2，B MCR的運(yùn)行工況下，分別模擬了燃燼風(fēng)配比（占總風(fēng)量的比例）為0%、5%、10%及15%4種工況。CO及NOx隨爐膛高度的分布規(guī)律，模擬結(jié)果見圖2。

如圖2（a）所示，4種工況下的CO濃度的變化趨勢(shì)相近。例如，4種情況最高峰都出現(xiàn)在燃燒器區(qū)域，且隨著燃燼風(fēng)比例的增加CO濃度峰值有所升高。由于燃燒器區(qū)域溫度很高，不添加燃燼風(fēng)時(shí)，爐膛氧氣量能夠充足地保證燃料充分燃燒，因此造成CO濃度較小。而隨著燃燼風(fēng)的添加，減小了燃燒器區(qū)域的氧氣量，甚至造成燃燒器區(qū)域的缺氧狀態(tài)，導(dǎo)致了燃燒緩慢，有助于CO生成，不利于已存在CO迅速氧化。此外，燃燒器上部，隨著爐膛高度升高，CO的截面濃度呈現(xiàn)出減小趨勢(shì)，最后到達(dá)爐膛出口時(shí)已經(jīng)很小了。

從圖2（b）可以看出，隨著燃燼風(fēng)比例的增加，NOx的排放量呈下降趨勢(shì)。由于是在B MCR工況保持總風(fēng)量不變的條件下進(jìn)行的，可能是因?yàn)橐徊糠侄物L(fēng)分出來改為燃燼風(fēng)加在燃燒器上面，減小爐膛下層氧氣量，形成富燃料的還原氣氛，降低了燃燒區(qū)域的溫度，減小了熱力型NOx的生成。此外，鍋爐上層氧氣含量也將有所增加，進(jìn)而降低了上層溫度，從而對(duì)NOx的形成也產(chǎn)生了一定的抑制作用。根據(jù)多次模擬經(jīng)驗(yàn)得出，分出來的燃燼風(fēng)部分不宜太大。燃燼風(fēng)過大，可能會(huì)將上層的氮進(jìn)一步氧化成NOx，反而增大NOx的排放，因此對(duì)于不同的燃燒爐需要選擇一個(gè)合適的燃燼風(fēng)比例，本爐子在燃燼風(fēng)為15%的情況下運(yùn)行NOx的排放最優(yōu)。

3.2 高爐煤氣配比的影響

煤價(jià)的上漲和鋼廠副產(chǎn)品高爐煤氣的難以二次利用使得多數(shù)電廠開始嘗試煤粉摻混高爐煤氣進(jìn)行混燃。在B MCR工況，過量空氣系數(shù)為1.2，燃燼風(fēng)比例為15%的條件下，分別模擬了高爐煤氣摻燒比為0%，10%，20%，30%4種工況。CO及NOx隨爐膛高度的分布規(guī)律，模擬結(jié)果見圖3。

從圖3（a）可以看出，純?nèi)济汗r的CO濃度最大，約為6%。隨著高爐煤氣摻燒量的增加，CO濃度峰值逐漸下降。比如，摻燒30%高爐煤氣工況降至4%左右。由于高爐煤氣的成分含有20%的CO，這可能是燃燒初期摻燒30%高爐煤氣工況比其他工況的CO濃度高的主要原因。隨著高爐煤氣摻燒比增加，CO濃度的峰值逐漸降低且向右偏移?？赡苁且?yàn)楦郀t煤氣與空氣易充分混合且相比于煤粉更易燃燒和燃盡。因此，隨著高爐煤氣摻混量的增加，爐膛內(nèi)燃料更易完全燃燒；而煤粉與空氣混合的充分性較差，致使燃燒不充分，從而CO的峰值較高。此外，高爐煤氣熱值低燃燒緩慢著火推遲導(dǎo)致峰值向右偏移，摻燒比例越大延遲越多。燃料在爐膛出口處均能基本燃盡［15］。

如圖3（b）所示，隨著高爐煤氣摻燒比的增加，爐膛內(nèi)整體NOx濃度分布明顯降低。例如，當(dāng)摻燒高爐煤氣量達(dá)到30%時(shí)，NOx排放量達(dá)到225×10－6，相對(duì)純煤粉的降低了一半。結(jié)合表2，高爐煤氣中的氮元素主要以N2的形式存在，而N2轉(zhuǎn)化為NOx的過程又相對(duì)困難，其可能是造成摻混高爐煤氣NOx濃度分布下降的主要原因。此外，由于高爐煤氣熱值較低，從而燃燒純煤時(shí)的爐膛溫度相比于摻混高爐煤氣后要高很多，致使熱力型NOx成倍增加。

3.3 過量空氣系數(shù)的影響

過量空氣系數(shù)是制約爐膛內(nèi)燃燼程度的主要因素。在高爐煤氣摻燒比例為20%、燃燼風(fēng)比例為15%、運(yùn)行工況為B MCR的情況下，分別對(duì)過量空氣系數(shù)1，1.1，1.15及1.2共4種工況進(jìn)行了模擬，其結(jié)果見圖4。

如圖4（a）4條CO濃度分布曲線所示，隨著過量空氣系數(shù)的增加，爐膛整體的CO濃度分布呈下降趨勢(shì)。這是因?yàn)?，過量空氣的增加使?fàn)t膛內(nèi)的氧量更為充足，從而有助于燃料的充分燃燒，最終表現(xiàn)為爐膛內(nèi)部整體CO濃度的下降。此外，過量空氣系數(shù)是1時(shí)爐膛出口CO濃度偏大，說明未能充分燃燒，而其余工況均較小。進(jìn)而也說明了選擇適當(dāng)過量空氣系數(shù)的必要性。從圖4（b）可以看出，隨著過量空氣系數(shù)的增加，爐膛出口NOx濃度排放量逐漸增加。例如，NOx濃度整體分布在過量空氣系數(shù)為1.2時(shí)最大，過量空氣系數(shù)為1時(shí)最小。由于過量空氣系數(shù)的增加有利于燃料燃燒更充分，進(jìn)而推測(cè)，在選取的過量空氣范圍內(nèi)，過量空氣系數(shù)越大，燃料燃燒越充分，使得爐膛內(nèi)部的溫度更高，從而促進(jìn)了熱力型NOx的生成。此外，過量空氣系數(shù)的增加，有利于NHi等氣體氧化成NOx，也可能是造成NOx增加的重要因素。由于燃料中的N生成NOx的轉(zhuǎn)化率隨空氣過量系數(shù)的增加呈上升趨勢(shì)，因此從降低污染的角度來講，在考慮最佳運(yùn)行條件的同時(shí)，也應(yīng)考慮過量空氣系數(shù)的選擇。

3.4 運(yùn)行工況的影響

電力需求的波動(dòng)性決定了運(yùn)行負(fù)荷的可調(diào)節(jié)性。因此對(duì)高爐煤氣比例為20%、燃燼風(fēng)為15%、過量空氣系數(shù)為1.2時(shí)，選取不同運(yùn)行工況（100%B MCR 鍋 爐 最 大 蒸 發(fā) 量、85%B MCR、70%B MCR）進(jìn)行模擬，其結(jié)果見圖5。

從圖5（a）可以看出，負(fù)荷變化對(duì)CO濃度影響不明顯，CO濃度峰值隨運(yùn)行負(fù)荷減小略有上升且均為燃燼風(fēng)出口附近。推測(cè)其原因可能為風(fēng)量隨著運(yùn)行負(fù)荷增大而增加，從而對(duì)CO有一定的稀釋作用，表現(xiàn)為CO濃度下降。如圖5（b）所示。隨著負(fù)荷的降低，NOx濃度呈下降趨勢(shì)?？赡苁且?yàn)?，隨著負(fù)荷降低，燃料量有所減少，燃料NOx排放減少，同時(shí)爐內(nèi)溫度有所降低，熱力NOx降低，所以總的氮氧化物減少，表現(xiàn)為NOx濃度分布下降。文獻(xiàn) ［16］也得到類似結(jié)論。

5 結(jié)論

通過對(duì)某鋼廠300 MW四角切圓煤粉／高爐煤氣混燃鍋爐進(jìn)行數(shù)值模擬，探討了燃燼風(fēng)配比、高爐煤氣配比、過量空氣系數(shù)及運(yùn)行負(fù)荷等因素對(duì)爐膛CO分布及NOx分布的影響規(guī)律。結(jié)果表明：燃燼風(fēng)比例增加，燃燒器附近的CO濃度峰值增加且爐膛內(nèi)部的NOx濃度整體降低；高爐煤氣摻混比增加，爐膛內(nèi)部的CO濃度的峰值下降，同時(shí)NOx濃度呈整體下降趨勢(shì)；隨著過量空氣系數(shù)增加，燃燒器上部的CO濃度均有所下降，而NOx濃度逐漸上升；運(yùn)行工況對(duì)爐膛內(nèi)部CO濃度的分布影響不大，隨著負(fù)荷的降低，峰值略有升高，然而NOx濃度分布在燃燒器上部受負(fù)荷影響較顯著，隨著負(fù)荷的下降，NOx濃度逐漸下降。

［1］汪洪濤，程艷玲，楊天助.220 t／h全燒高爐煤氣高溫高壓電站鍋爐的研究 ［J］.節(jié)能技術(shù)，2005，23（5）：422-425.

Wang Hongtao，Chen Yanling，Yang Tianzhu.Research of 220t／h high temperature and pressure plant boiler with combusting pure BFG ［J］.Energy Conservation Technology，2005，23 （5）：422-425.

［2］趙洪波，付麗梅，莊緒平.220 t／h煤粉鍋爐摻燒高爐煤氣改造 ［J］.電站系統(tǒng)工程，2003，19 （6）：43.

Zhao Hongbo，F(xiàn)u Limei，Zhuang Xuping.Retrofit of 220t／h coal-fired boiler blended with blast furnace Gas［J］.Power SystemEngineering，2003，19 （6）：43.

［3］劉定平.低熱值高爐煤氣與煤粉混燒技術(shù)的探討 ［J］.熱科學(xué)與技術(shù)，2003，2 （1）：74-79.

Liu Dingping.Study of mixing combustion techniques for coal and law heat value blast-furnace gas ［J］.Journal of Thermal Science and Technology，2003，2 （1）：74-79.

［4］湛志鋼，陳剛.煤粉鍋爐摻燒高爐煤氣對(duì)煤粉燃盡影響的研究 ［J］.動(dòng)力工程，2004，24 （2）：179-182.

Zhan Zhigang，Chen Gang.Influence of pulverized coal combustion boiler adding blast furnace gas on coal complete combustion ［J］.Power Engineering，2004，24（2）：179-182.

［5］Choeng R C，Chang N K.Numerical investigation on the flow，combustion and NOxemission characteristics in a 500 MWe tangentially fired pulverized-coal boiler［J］.Fuel，2009，88 （9）：1720-1731.

［6］李永華，陳鴻偉，劉吉臻，等.800 MW鍋爐混煤燃燒數(shù)值模擬 ［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2002，22（6）：101-104.

Li Yonghua，Chen Hongwei，Liu Jizhen，et al.Number simulation of blending coals combustion of 800 MWboiler［J］.Proceedings of the CSEE，2002，22 （6）：101-104.

［7］Yin Chungen，Caillat S，Harion J，et al.Investigation of the flow，combustion，heat-transfer and emissions from a 609 MWutility tangentially fired pulverized-coal boiler［J］.Fuel，2002，81 （8）：997-1006.

［8］Park H Y，F(xiàn)aulkner M，Turrell MD，et al.Coupled fluid dynamics and whole plant simulation of coal combustion in a tangentially-fired boiler ［J］.Fuel，2010，89 （8）：2001-2010.

［9］Belosevic S，Sijercic M，Oka S，et al.Three-dimensional modeling of utility boiler pulverized coal tangentially fired furnace ［J］.International Journal of Heat and Mass Transfer，2006，49 （19-20）：3371-3378.

［10］王政允，孫保民，郭永紅，等.330 MW前墻燃燒煤粉鍋爐爐內(nèi)溫度場(chǎng)的數(shù)值模擬及優(yōu)化 ［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2009，29 （20）：18-24.

Wang Zhengyun，Sun Baomin，Guo Yonghong，et al.Numerical simulation and optimization on temperature field of 330 MWfront wall fired furnace ［J］.Proceedings of the CSEE，2009，29 （20）：18-24.

［11］吳彥坤，高正陽，崔偉春，等.燃燼風(fēng)對(duì) NOx爐內(nèi)分布影響的數(shù)值模擬研究 ［J］.電力科學(xué)與工程，2008，24 （3）：63-66.

Wu Yankun，Gao Zhengyang，Cui Weichun，et al.Study on the effect of over fired air on NOxdistribution in furnace by numerical simulation ［J］.Electric Power Science and Engineering，2008，24 （3）：63-66.

［12］潘維，池作和，斯東坡，等.200 MW四角切圓燃燒鍋爐改造工況數(shù)值模擬 ［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2005，25 （8）：110-115.

Pan Wei，Chi Zuohe，Si Dongpo，et al.Numerical Simulation of combustion process in a 200 MWtangentially fired furnace to study furnace reconstruction［J］.Proceedings of the CSEE，2005，25 （8）：110-115.

［13］Sivathanu Y R，F(xiàn)aeth G M.Generalized state relationships for scalar properties in non-premixed hydrocar-bon／ air flames ［J］.Combustion and Flame，1990，82 （2）：211-230.

［14］高正陽，宋瑋，方立軍，等.1000 MW超超臨界機(jī)組雙切圓鍋爐NOx排放特性的數(shù)值模擬 ［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2009，29 （32）：12-18.

Gao Zhengyang，Song Wei，F(xiàn)ang Lijun，et al.Numerical simulation of NOxemission of a dual circle tangential firing boiler of a 1 000 MWultra supercritical unit［J］.Proceedings of the CSEE，2009，29 （32）：12-18.

［15］楊軼，陳剛.煤粉和高爐煤氣混燒鍋爐燃燒問題的分析及改 造 ［J］.電站系統(tǒng)工程，2003，19 （2）：36-38.

Yang Yi，Chen Gang.Technical reform of combustion system in pulverized coal and BFG combustion boiler ［J］Power System Engineering，2003，19（2）：36-38.

［16］王學(xué)棟，辛洪昌，欒濤，等.330 MW機(jī)組鍋爐燃燒調(diào)整對(duì)NOx排放濃度影響的試驗(yàn)研究 ［J］.電站系統(tǒng)工程，2007，23 （3）：7-10.

Wang Xuedong，Xin Hongchang，Luan Tao，et al.Research and test on Influence of boiler combustion adjusting on NOxemission of 330 MWUnit ［J］.Power System Engineering，2007，23 （3）：7-10.

Numerical Simulation of 300 MWBlast Furnace Gas and Pulverized Coal Blended Combustion Boiler on the Optimizing of NOx

Li Wei，Wang Jinxing，Huang Jiangcheng，Wang Chunbo
（School of Energy Power and Mechanical Engineering，North China Electric Power University，Baoding 071003，China）

300 MWblast furnace gas and pulverized coal blended combustion Boiler was numerical simulated by the fluent software.The effects of OFA ratio，blast furnace gas ratio，excess air coefficient and operating duty on CO concentration distribution and NOxconcentration distribution in the furnace cavity were studied.Some results and conclusions were illustrated through some curves.The more OFA ratio，the higher peak value of CO concentration in the burner area and the less NOxconcentration in the whole furnace cavity；The increase of blast furnace gas ratio made the peak value of CO concentration descend and meanwhile made the NOxconcentration decline integrally；The increase of excess air coefficient always can lead to the decreasing of CO concentration over the burner，and give rise to the gradual rising of NOxconcentration.The effect of operating duty on CO concentration distribution of the furnace cavity is tiny.With the operating duty decreases，the peak value rises weakly.Besides，the effect of operating duty on NOxconcentration distribution over the burner is obvious，the NOxconcentration decreases along with the reducing of operating duty.

coal；BFG；blended combustion boiler；NOx；CO；numerical simulation.

TK 229

A

2012－05－23。

教育部新世紀(jì)優(yōu)秀人才支持計(jì)劃 （NCET－08－0770）；中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng) （重點(diǎn)）資金資助項(xiàng)目（11ZG12）。

李偉 （1987-），男，碩士研究生，研究方向?yàn)闈崈裘喝紵?，E-mail：snzx198748＠126.com。