張俊峰

【摘 要】以某100t/h鍋爐為原型搭建冷態(tài)實(shí)驗(yàn)臺(tái)，研究再燃燃燒器噴口流場特性及其對(duì)爐內(nèi)上升氣流覆蓋效果的影響，并采用標(biāo)準(zhǔn)k-ε湍流模型對(duì)其進(jìn)行數(shù)值模擬。結(jié)果表明：再燃區(qū)采用八個(gè)噴口比四個(gè)噴口更能使?fàn)t膛內(nèi)氣流充滿度強(qiáng)，再燃?xì)饬鲗?duì)爐膛內(nèi)上升氣流覆蓋效果更好；在各噴口氣流速度不變的工況下，爐膛四角入射角為14.17°、側(cè)墻中心處入射角為84.57°再燃?xì)饬鲗?duì)上二次風(fēng)的覆蓋效果最好；在保持該入射角度不變的情況下，四角再燃?xì)饬魉俣葹?0m/s，側(cè)墻中心處再燃?xì)饬魉俣葹?0m/s時(shí)（工況3），再燃?xì)饬鲗?duì)上二次風(fēng)的覆蓋效果最好。

【關(guān)鍵詞】煤粉爐；再燃噴口；冷態(tài)流場；數(shù)值模擬

再燃是一種很有前途的脫氮方法[1]，因其投運(yùn)成本低、脫氮效率高等優(yōu)點(diǎn)而受到廣泛重視。天然氣由于本身不含氮、灰和硫，且能比其它燃料產(chǎn)生更多的烴根，被認(rèn)為是最理想、使用最為廣泛的再燃燃料。近年來許多國家都對(duì)天然氣再燃降低NOx機(jī)理的研究發(fā)現(xiàn)，再燃燃料中含氮組分HCN對(duì)再燃過程有著重要的影響。Patry和Engel[2]在研究氮氧化物和甲烷的反應(yīng)后認(rèn)為其生成產(chǎn)物為HCN、水和氫氣，發(fā)現(xiàn)約有70%的NOx在相對(duì)短的時(shí)間內(nèi)實(shí)現(xiàn)了轉(zhuǎn)化。該項(xiàng)技術(shù)的關(guān)鍵之一是再燃?xì)怏w燃料與爐內(nèi)煙氣的混合特性[3]，所以，研究再燃噴口射流流場特性對(duì)探討爐內(nèi)氣體混合特點(diǎn)、及其對(duì)降低NOx的影響有重要意義。

本文考慮到僅由爐膛四角射入的天然氣與爐內(nèi)煙氣混合程度較差，降低再燃還原NOx的效果，因此，設(shè)計(jì)時(shí)在四側(cè)爐墻中心處各增加一個(gè)再燃噴口以提高再燃?xì)怏w在爐內(nèi)的充滿度，實(shí)現(xiàn)再燃區(qū)八點(diǎn)噴射，達(dá)到有效降低NOx的目的。文中采用相似?；碚撚盟俣缺碚鳚舛?，將兩相流處理為單相流體，通過數(shù)值模擬與冷態(tài)實(shí)驗(yàn)研究的方法，研究再燃區(qū)噴口數(shù)目、氣流入射角度，再燃?xì)饬魉俣鹊葘?duì)爐內(nèi)流場影響，為優(yōu)化天然氣再燃技術(shù)及工程應(yīng)用提供參考。

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

1.1 實(shí)驗(yàn)裝置

以某100t/h鍋爐為原型，根據(jù)1：5的比例搭建冷態(tài)實(shí)驗(yàn)臺(tái)，實(shí)驗(yàn)系統(tǒng)由三部分組成：鍋爐爐膛，燃燒器組件和送風(fēng)管路。整個(gè)爐膛用有機(jī)玻璃制造，燃燒器組件由下至上依次為下一次風(fēng)噴口、下二次風(fēng)噴口、上一次噴口、上二次風(fēng)噴口、再燃?xì)怏w噴口和燃盡風(fēng)噴口，再燃區(qū)八個(gè)噴口在爐膛同一截面上分布示意圖如圖1所示。本實(shí)驗(yàn)的工質(zhì)為空氣，測(cè)量儀器為熱球風(fēng)速儀，實(shí)驗(yàn)中測(cè)試區(qū)及爐膛一角測(cè)點(diǎn)坐標(biāo)布置如圖2所示：x軸為以噴口中心為原點(diǎn)的氣流軸線；y軸為沿燃燒器噴口的水平軸線；O點(diǎn)為燃燒器噴口中心。

1.2 實(shí)驗(yàn)參數(shù)

根據(jù)相似理論，選擇鍋爐實(shí)際運(yùn)行參數(shù)進(jìn)行冷態(tài)?；?jì)算，可得到冷態(tài)實(shí)驗(yàn)再燃各噴口出口風(fēng)速，冷態(tài)實(shí)驗(yàn)參數(shù)與實(shí)際運(yùn)行參數(shù)對(duì)照表如表1所示。

根據(jù)研究內(nèi)容，本實(shí)驗(yàn)安排了6種工況進(jìn)行測(cè)試，各工況冷態(tài)實(shí)驗(yàn)參數(shù)如表2所示。

2 實(shí)驗(yàn)與數(shù)值計(jì)算理論基礎(chǔ)

2.1 實(shí)驗(yàn)原理

根據(jù)相似模化理論，冷態(tài)實(shí)驗(yàn)必須遵循以下原則[4]：（1）模型與實(shí)物保持幾何相似；（2）模型與實(shí)物在對(duì)應(yīng)的工況下，保證氣流運(yùn)動(dòng)狀態(tài)進(jìn)入自?；瘏^(qū)；（3）模型與實(shí)物的各股氣流動(dòng)量比相等。

由于模型與實(shí)物的一二次風(fēng)動(dòng)量比相等[[5-6]，即：

根據(jù)實(shí)物不同噴口的運(yùn)行參數(shù)，可得到冷態(tài)模型的各噴口參數(shù)。

2.2 數(shù)值計(jì)算數(shù)學(xué)模型

鍋爐內(nèi)的氣體流動(dòng)是一種三維湍流流動(dòng)，大量的研究表明，采用標(biāo)準(zhǔn)k-ω雙方程模型具有較好的適應(yīng)性。在三維直角坐標(biāo)系中，根據(jù)Navier-Stokes方程，可得到氣相湍流控制方程組，其通用形式為[7]：

入口邊界條件：一次風(fēng)V1=11.4m/s，二次風(fēng)V2=22.6m/s，再燃風(fēng)V3=51m/s，燃盡風(fēng)V4=29m/s；出口邊界條件：爐膛出口負(fù)壓取運(yùn)行值-50Pa；壁面邊界條件：無速度滑移和無質(zhì)量滲透條件。其中，邊界上的湍動(dòng)能k=■I■■，湍動(dòng)能耗散率ε=0.093/4■，湍流強(qiáng)度I=0.16Re-1/8，湍流長度尺寸l=0.007L。

本次冷態(tài)鍋爐數(shù)值模擬采用SIMPLE算法對(duì)各個(gè)變量的方程組進(jìn)行迭代求解[8]，SIMPLE算法是一種求解壓力耦合方程組的半隱式方法，主要用于不可壓縮流場的數(shù)值模擬計(jì)算。SIMPLE算法主要利用“猜測(cè)-修正”的方法，在交錯(cuò)網(wǎng)格的基礎(chǔ)上來計(jì)算壓力場，從而達(dá)到求解動(dòng)量方程的目的。

3 結(jié)果分析

實(shí)驗(yàn)測(cè)量了6種工況下，爐膛四角處再燃噴口和上二次風(fēng)噴口截面上5組無量綱距離X/b（X/b=1，2.5，5， 7.5，10，噴口的寬度b為2cm）測(cè)點(diǎn)上的風(fēng)速，在每組X/b測(cè)點(diǎn)上，在噴口中心左右側(cè)每隔2cm進(jìn)行速度測(cè)量，將測(cè)量的結(jié)果與數(shù)值模擬的結(jié)果進(jìn)行對(duì)比，本文重點(diǎn)分析爐膛中心處（X/b=7.5，10）流場分布情況。

3.1 再燃噴口數(shù)目對(duì)再燃?xì)饬鞲采w效果的影響

大量研究表明沿爐膛水平截面，爐內(nèi)氣流濃度剖面分布同速度剖面分布將存在相似性[9]，因此本次實(shí)驗(yàn)中可以用速度的分布來表征濃度的分布[10]。

工況1與工況2中再燃?xì)饬鞲咚賴娙霠t膛，在噴口前端（X/b=1，2.5，5），由于風(fēng)速大，氣流剛性強(qiáng)，再燃?xì)饬髋c上二次風(fēng)的速度中心軸線大致重合，再燃?xì)饬魉俣仍?2

由圖3（a）可知，當(dāng)氣流接近爐膛中心處（X/b=7.5）時(shí)，工況1由于再燃?xì)饬魉俾食霈F(xiàn)衰減，再燃?xì)饬髟跔t內(nèi)旋轉(zhuǎn)氣流沖擊下向鍋爐爐墻一側(cè)偏斜[11]，其速度最大值出現(xiàn)在Y/b=1附近；相比于工況2的氣流偏斜率要小。由圖3（b）可知，工況1再燃?xì)饬鞯淖畲笾禐?3.1m/s，不能明顯高于上二次風(fēng)速，表示該處再燃?xì)饬鞲采w效果較差。而在工況2中，由于鍋爐四側(cè)墻中心處增設(shè)了再燃噴口，通過周圍氣流的卷吸作用，從而保證了四角側(cè)再燃風(fēng)有足夠大的動(dòng)量，能夠射入爐膛旋轉(zhuǎn)氣流中心，對(duì)上二次風(fēng)的覆蓋率更高。

3.2 側(cè)墻中心處再燃?xì)饬魅肷浣嵌葘?duì)覆蓋效果的影響

對(duì)于燃煤鍋爐，切圓直徑是爐內(nèi)氣流流動(dòng)狀態(tài)的重要特征[8]，理想的爐內(nèi)氣流流動(dòng)狀態(tài)是在爐膛中心形成合適的圓形旋轉(zhuǎn)火焰[12-13]。

在保證再燃風(fēng)速及上二次風(fēng)速不變的情況下，改變側(cè)墻中心處再燃?xì)饬魅肷浣嵌?，?0°逐漸調(diào)整至74.57°，由圖4（a）可以看出，在遠(yuǎn)離噴口的位置（x/b=7.5）處，隨著入射角度的降低，工況2，3，4的再燃?xì)饬髦行奶幾畲笏俣赛c(diǎn)分別出現(xiàn)在y/b=0.5，1，1.5左右，該現(xiàn)象表明在不同的入射角度下，由于四側(cè)墻中心處再燃?xì)饬髋c四角側(cè)的再燃?xì)饬骰旌蠌?qiáng)度不同，從而引起四角側(cè)再燃?xì)饬髌甭实牟町?，最終形成不同直徑的切圓。

工況4再燃?xì)饬髦行奶幾畲笏俣赛c(diǎn)出現(xiàn)在y/b=1.5，再燃?xì)饬髌甭首畲?，在爐內(nèi)形成較大的切圓，切圓直徑過大，使再燃?xì)饬魅菀踪N墻，在實(shí)際運(yùn)行中易造成水冷壁的結(jié)渣[14]。

根據(jù)圖4（b）可知，工況2，3中再燃?xì)饬魉俣仍?2

3.3 側(cè)墻中心處再燃?xì)饬魉俣葘?duì)覆蓋效果的影響

通過研究側(cè)墻中心處再燃噴口不同的入射角度對(duì)爐內(nèi)再燃?xì)饬鞲采w效果的影響發(fā)現(xiàn)，當(dāng)四角側(cè)再燃?xì)饬魅肷浣嵌葹?1.17°，四墻側(cè)再燃?xì)饬魅肷浣嵌葹?4.57°時(shí)，再燃?xì)饬髂芎芎玫母采w上二次風(fēng)射流。在該入射角度不變的前提下，改變其入射氣流速度，其結(jié)果如下。

由圖5可知，工況5下，四角再燃噴口（X/b=1，2.5，5））氣流高速噴入爐膛，剛性較強(qiáng)，計(jì)算值與模擬值吻合度較好，但是由于側(cè)墻中心處再燃?xì)饬饕?4.57°射入爐膛，對(duì)四角側(cè)再燃?xì)饬鳟a(chǎn)生了一個(gè)橫向推動(dòng)力，兩股氣流相互卷吸，抵消了部分動(dòng)量，故四角再燃?xì)饬餮刂鴩娍谳S線衰減速度加快，在爐膛中間區(qū)域（x/b=7.5，10），再燃?xì)饬魉俣葻o法完全明顯的高于上二次風(fēng)氣流流速，在實(shí)際運(yùn)行時(shí)，爐膛的高溫火焰處NOX生成含量相對(duì)較高，如果再燃?xì)饬髟谠搮^(qū)域無法完全包圍住爐膛內(nèi)旋轉(zhuǎn)氣流，容易造成脫硝效率的下降。

工況3下，在噴口出口處（x/b=1，2.5，5）四角再燃?xì)饬髋c上二次風(fēng)氣流的最高速度中軸線大致保持吻合，并且在-2

在工況6下，由于側(cè)墻中心處再燃?xì)饬黠L(fēng)速大，剛性強(qiáng)，在沿著噴口軸線上對(duì)四角再燃?xì)饬鳟a(chǎn)生了劇烈的沖擊作用，雖然在此工況下，四角側(cè)再燃?xì)饬髟?2

4 結(jié)論

（1）天然氣再燃技術(shù)中，再燃燃燒器采用分別在爐膛四角和側(cè)墻中心處布置八個(gè)噴口比僅在爐膛四角布置四個(gè)噴口的降低NOx效果好，八點(diǎn)噴射的再燃?xì)饬鲗?duì)爐內(nèi)上升氣流覆蓋效果更好，有利于脫硝效率的提高；

（2）再燃燃燒器采用八個(gè)噴口時(shí)，再燃?xì)饬魅肷浣嵌葘?duì)于爐內(nèi)上升氣流的覆蓋效果有較大的影響，其中，在各噴口再燃?xì)饬魉俣炔蛔兊那闆r下，爐膛四角氣流以41.17°射入，側(cè)墻中心處氣流以84.57°射入，再燃?xì)饬鲗?duì)上二次風(fēng)的覆蓋效果最好。

（3）再燃?xì)饬魉俣鹊拇笮?duì)降低NOx也有一定的影響，在保持四角噴口入射角為41.17°、側(cè)墻噴口入射角為84.57不變的情況下，四角側(cè)再燃?xì)饬魉俣葹?1m/s、四側(cè)墻中心處再燃?xì)饬魉俣葹?1m/s時(shí)，再燃?xì)饬鲗?duì)上二次風(fēng)的覆蓋效果最好。

【參考文獻(xiàn)】

[1]Mereb J， Wendt J O L. Reburning mechanisms in a pulverized coal combustor[C]//Twenty-third Symposium （International） on Combustion. The Combustion Institute， 1990，1273-1279.

[2]Patry M.，Engel G. Formation of HCN by the action of nitric oxide on methane atatmospheric pressure， 1. General conditions of formation[Z]. Compt. Rend.， 1950，231：1302-1304.

[3]姚向東，張忠孝，邱莉莉，等.天燃?xì)庠偃冀档蚇Ox的化學(xué)動(dòng)力學(xué)及影響因素分析[J].上海理工大學(xué)學(xué)報(bào)，2004，26（1）：62-65.

[4]李之光.相似與?；痆M].北京：國防工業(yè)出版社，1982.

[5]岑可法.鍋爐燃燒試驗(yàn)研究方法及測(cè)量技術(shù)[M].水利電力出版社，1987.

[6]閻維平，劉亞芝，黃景立.2000t/h四墻切圓鍋爐爐內(nèi)空氣動(dòng)力場試驗(yàn)研究[J].電站系統(tǒng)工程，2010（09）.

[7]陶文銓.數(shù)值傳熱學(xué).2版[M].西安：西安交通大學(xué)出版社，2001.

[8]周托，張忠孝，烏曉江，等.氣體再燃燃燒器噴口流動(dòng)特性冷態(tài)實(shí)驗(yàn)研究[J].上海理工大學(xué)學(xué)報(bào)，2007（2）.

[9]董志勇.射流力學(xué)[M].北京：科學(xué)出版社，2005.

[10]吳雙應(yīng)，李友榮，盧嘯風(fēng).采用噴氣再燃降低NOX排放時(shí)單噴嘴射流的冷態(tài)分析[J].工業(yè)加熱，2001（6）：4-7.

[11]張澤，吳少華，秦裕琨，等.爐內(nèi)噴嘴射流剛性的理論模型及試驗(yàn)研究[J].中國電機(jī)工程學(xué)報(bào)，2001，21（9）：104-109.

[12]王瑩，秦裕琨，吳少華，邢繼春.爐內(nèi)流場對(duì)水冷壁高溫腐蝕影響的數(shù)值模擬分析[J].熱能動(dòng)力工程，2000（03）.

[13]朱彤，范衛(wèi)東，信偉，陳力哲，吳少華，秦裕琨.雙爐膛爐內(nèi)空氣動(dòng)力場的數(shù)值模擬研究[J].熱能動(dòng)力工程，1997（06）.

[14]李永華，司金茹.不同擺角下四墻切圓燃燒器的數(shù)值模擬[J].中國電機(jī)工程學(xué)報(bào)，2011（11）.

[責(zé)任編輯：湯靜]