王小龍，張飛龍，王 里，許鑫瑋，張儀祥，王學(xué)斌，譚厚章

(1.國(guó)能神東煤炭集團(tuán)有限責(zé)任公司，陜西 神木 719315；2.西安交通大學(xué) 能源與動(dòng)力工程學(xué)院，陜西 西安 710049)

0 引 言

據(jù)統(tǒng)計(jì)，我國(guó)目前使用的工業(yè)鍋爐有60多萬(wàn)臺(tái)，其中有接近50萬(wàn)臺(tái)為鏈條爐等燃煤工業(yè)鍋爐[1-2]。傳統(tǒng)的燃煤工業(yè)鍋爐通常采用爐排燃燒，由于有效反應(yīng)面積小，爐排鍋爐存在燃燒效率低、污染物排放高等嚴(yán)重問(wèn)題。煤粉工業(yè)鍋爐強(qiáng)化了煤與空氣的接觸與混合，能大幅提高煤燃燒效率，有效控制污染物排放[3-5]。

然而由于燃煤性質(zhì)和運(yùn)行負(fù)荷的變化，煤粉鍋爐易出現(xiàn)燃燒不穩(wěn)定和NOx排放量高的問(wèn)題。目前，普遍使用的低氮燃燒技術(shù)是利用燃燒器的特殊結(jié)構(gòu)設(shè)計(jì)改變煤粉與空氣的階段性配比，將濃淡燃燒、空氣分級(jí)等低氮燃燒技術(shù)應(yīng)用于煤粉燃燒。高濃度煤粉的低氮旋流燃燒器能創(chuàng)建一個(gè)高溫再循環(huán)區(qū)，促進(jìn)煤粉再著火，穩(wěn)定燃燒，還能降低著火區(qū)的溫度和含氧量，形成還原性氣氛，降低NOx的生成[6-9]。預(yù)燃室設(shè)計(jì)是為了保證燃用貧煤和煙煤時(shí)爐內(nèi)燃燒穩(wěn)定高效[10]。

20世紀(jì)80年代，清華大學(xué)與上海吳涇電廠共同研制了使用根部二次風(fēng)的預(yù)燃室旋流燃燒器，緩解了預(yù)燃室前端與噴火口的積碳結(jié)焦問(wèn)題[11]。王進(jìn)卿等[12]使用商業(yè)軟件對(duì)一種新開(kāi)發(fā)的采用鈍體的預(yù)燃室燃燒器進(jìn)行了數(shù)值計(jì)算，研究一、二次風(fēng)速，二次風(fēng)旋流強(qiáng)度等對(duì)回流強(qiáng)度的影響。該燃燒器在煤粉工業(yè)鍋爐應(yīng)用中的NOx排放濃度為331 mg/Nm3[13]。陳志超等[14]通過(guò)將預(yù)燃室燃燒器的二次風(fēng)分為內(nèi)外兩級(jí)和煙氣再循環(huán)來(lái)弱化前期的一、二次風(fēng)混合，降低初期NOx生成，最終降低NOx排放。

高濃度煤粉燃燒器可以直接調(diào)控煤粉濃度控制燃燒，能明顯縮短著火時(shí)間和穩(wěn)定燃燒，阻礙NOx生成，是一種經(jīng)濟(jì)環(huán)保的燃燒技術(shù)[15]。但煤粉濃度提高意味著氣固兩相混合變差以及燃盡較差問(wèn)題加劇[16]。因此一次風(fēng)速的設(shè)計(jì)對(duì)高濃度煤粉鍋爐的安全經(jīng)濟(jì)運(yùn)行非常重要[17]。一次風(fēng)速過(guò)大或過(guò)小都會(huì)造成爐內(nèi)著火延遲、NOx排放增加、燃燒效率降低和燃燒不穩(wěn)定等問(wèn)題，同時(shí)，一次風(fēng)速過(guò)大會(huì)引起爐內(nèi)結(jié)焦，過(guò)小則會(huì)引起一次風(fēng)管堵塞和噴口燒壞[18-21]。目前對(duì)于高濃度煤粉的預(yù)燃室燃燒器研究仍較少，本文對(duì)新設(shè)計(jì)的一種帶預(yù)燃室的高濃度煤粉低氮旋流燃燒器進(jìn)行數(shù)值模擬，并用一臺(tái)25 t/h 全尺寸煤粉工業(yè)鍋爐試驗(yàn)臺(tái)對(duì)數(shù)值模擬的基本工況進(jìn)行試驗(yàn)驗(yàn)證，通過(guò)分析煤粉鍋爐的速度、溫度以及煙氣組分分布，研究該燃燒器的燃燒性能。

1 燃燒器幾何模型

基于本預(yù)燃式低氮旋流燃燒器的實(shí)際尺寸和結(jié)構(gòu)，使用Inventor軟件1∶1建立燃燒器的三維模型，燃燒器結(jié)構(gòu)如圖1所示。基本原理為：將高濃度煤粉與一次風(fēng)混合通入，同時(shí)二次風(fēng)分為2股，一股為內(nèi)二次風(fēng)，以直流方式進(jìn)入預(yù)燃室，另一股為外二次風(fēng)，經(jīng)過(guò)可調(diào)旋流葉片進(jìn)入爐內(nèi)；一次風(fēng)出口裝有一圈齒狀的穩(wěn)焰環(huán)，用于穩(wěn)定一次風(fēng)的著火。

圖1 低氮燃燒器示意Fig.1 Schematic diagram of low-nitrogen burner

爐膛燃燒區(qū)實(shí)際為一個(gè)L型矩形區(qū)域，但當(dāng)爐膛足夠大時(shí)，具體形狀對(duì)低氮燃燒器出口附近的流場(chǎng)影響不大，為簡(jiǎn)化計(jì)算，取一個(gè)直徑5 m、高10 m的圓柱體區(qū)域，并在該區(qū)域出口處設(shè)置一個(gè)圓臺(tái)的收束結(jié)構(gòu)，以縮小計(jì)算區(qū)域出口尺寸，避免過(guò)大的回流影響爐膛內(nèi)的流場(chǎng)分布。本文采用ICEM對(duì)幾何模型進(jìn)行網(wǎng)格劃分，對(duì)整個(gè)燃燒器區(qū)域使用六面體結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分，并對(duì)燃燒器出口、預(yù)燃室內(nèi)部及爐膛主要燃燒區(qū)域進(jìn)行加密(圖2)。

圖2 網(wǎng)格數(shù)劃分示意Fig.2 Schematic diagram of grid number division

2 數(shù)值模擬設(shè)置及正確性驗(yàn)證

2.1 數(shù)值模擬設(shè)置

本文使用ANSYS Fluent軟件進(jìn)行數(shù)值模擬。① 湍流模型選擇Realizablek-ε模型，該模型適用于圓形旋轉(zhuǎn)射流的模擬，更適用于反應(yīng)帶旋的湍流流動(dòng)，與本文的旋流燃燒器實(shí)際流場(chǎng)吻合；② 對(duì)煤粉顆粒的運(yùn)動(dòng)過(guò)程，本文選用拉格朗日隨機(jī)軌道模型進(jìn)行模擬；③ 煤粉燃燒過(guò)程中使用適用性廣的DO模型；④ 煤粉燃燒模型包括揮發(fā)分析出模型、揮發(fā)分氣相燃燒模型以及焦炭燃燒模型3部分：對(duì)揮發(fā)分析出過(guò)程，本文采用兩步競(jìng)爭(zhēng)反應(yīng)模型模擬揮發(fā)分析出過(guò)程，認(rèn)為揮發(fā)分的揮發(fā)速率受揮發(fā)分含量的影響；揮發(fā)分氣相燃燒過(guò)程使用渦耗散模型進(jìn)行模擬；焦炭燃燒過(guò)程選用擴(kuò)散-動(dòng)力反應(yīng)模型進(jìn)行模擬[22-23]。

預(yù)燃室采用耐火材料，預(yù)燃室壁面設(shè)置為絕熱壁面邊界，其他主要配風(fēng)參數(shù)見(jiàn)表1，煤粉參數(shù)見(jiàn)表2。此次驗(yàn)證試驗(yàn)的負(fù)荷為最大負(fù)荷的20%。邊界條件設(shè)置完成后，先對(duì)爐膛的冷態(tài)流場(chǎng)進(jìn)行迭代計(jì)算，得到充分發(fā)展的流場(chǎng)后，再耦合顆粒相，加入輻射和燃燒模型進(jìn)行熱態(tài)計(jì)算。

表1 試驗(yàn)工況主要參數(shù)

表2 煤工業(yè)分析和元素分析

2.2 網(wǎng)格無(wú)關(guān)性驗(yàn)證

本文選擇了21.8萬(wàn)、35.9萬(wàn)、118.2萬(wàn)、181.7萬(wàn)4種不同網(wǎng)格數(shù)量的網(wǎng)格進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證，得到不同網(wǎng)格數(shù)量下燃燒器出口軸線上的速度分布如圖3所示。計(jì)算結(jié)果表明，當(dāng)網(wǎng)格數(shù)目超過(guò)118.2萬(wàn)后，燃燒器出口軸線上的速度分布基本不變，綜合考慮計(jì)算效率與計(jì)算精度，本文選擇數(shù)目為118.2萬(wàn)的網(wǎng)格進(jìn)行模擬計(jì)算。

圖3 不同網(wǎng)格數(shù)模擬結(jié)果比較Fig.3 Simulation results of different grid numbers

2.3 模型正確性驗(yàn)證

進(jìn)行不同工況下燃燒器運(yùn)行模擬的數(shù)值計(jì)算前，首先對(duì)25 t/h煤粉工業(yè)鍋爐研究平臺(tái)上的典型煤粉燃燒工況進(jìn)行了模擬計(jì)算，將模擬計(jì)算結(jié)果與試驗(yàn)測(cè)量的溫度場(chǎng)和NOx濃度等關(guān)鍵數(shù)據(jù)進(jìn)行對(duì)比，以驗(yàn)證本研究所選取的數(shù)值模型的有效性，其燃燒器實(shí)際爐膛結(jié)構(gòu)與溫度測(cè)點(diǎn)如圖4所示。驗(yàn)證試驗(yàn)工況下的內(nèi)二次風(fēng)葉片角度為45°，該試驗(yàn)系統(tǒng)實(shí)際運(yùn)行時(shí)維持爐膛負(fù)壓穩(wěn)定在-100 Pa左右。

圖4 實(shí)際爐膛結(jié)構(gòu)與溫度測(cè)點(diǎn)Fig.4 Actual furnace structure and temperaturemeasurement points

測(cè)點(diǎn)溫度的模擬結(jié)果與實(shí)際測(cè)量結(jié)果對(duì)比如圖5所示。

由圖5可知，試驗(yàn)測(cè)量溫度與數(shù)值計(jì)算結(jié)果的最大偏差為52 ℃，誤差范圍±4.5%，小于工業(yè)上常規(guī)要求的10%。因此，可認(rèn)為本研究所選取的模型能夠較準(zhǔn)確預(yù)測(cè)爐膛內(nèi)的溫度分布。

圖5 測(cè)點(diǎn)溫度試驗(yàn)結(jié)果與模擬結(jié)果對(duì)比Fig.5 Comparison between experimental results andsimulation results of measuring point temperature

該工況下?tīng)t膛出口參數(shù)的模擬值與試驗(yàn)值對(duì)比見(jiàn)表3。模擬爐膛出口截面NOx濃度誤差為3.8%；模擬焦炭轉(zhuǎn)化率95.1%，折算飛灰含碳量16.2%，與實(shí)際值的誤差為8.7%，因此可認(rèn)為熱態(tài)模擬結(jié)果與測(cè)量結(jié)果符合較好，選擇設(shè)置的計(jì)算模型與邊界條件合理，模擬結(jié)果能夠較好地與實(shí)際燃燒情況吻合。

表3 飛灰含碳量及NOx模擬結(jié)果與試驗(yàn)結(jié)果對(duì)比

3 計(jì)算工況

本文模擬的3種工況中的煤粉參數(shù)與第2部分驗(yàn)證試驗(yàn)中所用煤種相同，具體見(jiàn)表2，所有模擬工況運(yùn)行負(fù)荷為最大負(fù)荷的70%。在保證空氣過(guò)量系數(shù)不變情況下，改變基礎(chǔ)工況的一次風(fēng)速，形成新的工況1和工況2，3種工況主要參數(shù)見(jiàn)表4。

表4 模擬工況主要參數(shù)

4 結(jié)果與分析

4.1 速度分析

中心截面軸向速度(vz)分布云圖如圖6所示，可知該燃燒器的流場(chǎng)特點(diǎn)為：該燃燒器內(nèi)二次風(fēng)射流具有較大的旋轉(zhuǎn)慣性和較高的速度，能夠形成較大的回流區(qū)?；亓鲄^(qū)可以分成2部分，一部分是預(yù)燃室內(nèi)的環(huán)形回流區(qū)，另一部分是爐膛中心位置的大回流區(qū)。前者主要作用是保證煤粉穩(wěn)定燃燒，后者主要作用是卷吸高溫?zé)煔饧訜崦悍?，同時(shí)提供低氧的還原性氣氛，將部分燃燒生成的NO還原成N2，抑制NOx生成。

圖6 中心截面軸向速度分布Fig.6 Axial velocity distribution in the central section

對(duì)比3個(gè)工況預(yù)燃室內(nèi)軸向速度分布，發(fā)現(xiàn)3個(gè)工況預(yù)燃室內(nèi)的環(huán)形回流區(qū)面積沒(méi)有明顯變化，但隨著一次風(fēng)速的增加，預(yù)燃室附近的環(huán)形回流區(qū)逆向速度增大，即回流流量增大，這有利于煤粉的穩(wěn)定著火。發(fā)生這種現(xiàn)象的原因是隨著一次風(fēng)射流速度的增大，環(huán)形回流區(qū)內(nèi)的負(fù)壓增大，環(huán)形回流區(qū)對(duì)內(nèi)二次風(fēng)的卷吸能力增強(qiáng)。預(yù)燃室內(nèi)流場(chǎng)示意如圖7所示，該回流區(qū)卷吸部分二次風(fēng)及爐膛中心的高溫?zé)煔?，并攜帶部分一次風(fēng)及煤粉回到預(yù)燃室內(nèi)劇烈燃燒，燃燒產(chǎn)生的煙氣一部分與內(nèi)二次風(fēng)混合加熱內(nèi)二次風(fēng)，剩下的與一次風(fēng)混合，加熱一次風(fēng)射流內(nèi)的煤粉，保證煤粉的穩(wěn)定著火與燃燒。

圖7 基礎(chǔ)工況預(yù)燃室內(nèi)流場(chǎng)示意Fig.7 Schematic diagram of the flow field in thepre-combustion chamber under basic conditions

對(duì)比3個(gè)工況爐膛中心大回流區(qū)面積發(fā)現(xiàn)：隨著一次風(fēng)速的增加，爐膛中心回流區(qū)面積逐漸減小，工況二的一次風(fēng)剛性過(guò)強(qiáng)，煤粉射流越深入爐膛內(nèi)部，使?fàn)t膛中心的大回流區(qū)遭到破壞。燃燒器出口軸線上軸向速度分布如圖8所示，工況1和基礎(chǔ)工況軸向速度衰減快，有明顯的回流區(qū)，而工況2的中心區(qū)域沒(méi)有回流。

圖8 軸線上軸向速度分布對(duì)比Fig.8 Comparison of axial velocity distribution on the axis

4.2 溫度分析

燃燒器出口中心線的溫度分布如圖9(a)所示，爐膛沿程截面平均溫度分布如圖9(b)所示。可知隨著一次風(fēng)速增大，預(yù)燃室內(nèi)的高溫區(qū)逐漸向預(yù)燃室深處移動(dòng)，分析該現(xiàn)象的成因?yàn)椋弘S著一次風(fēng)速的增大，卷吸的內(nèi)二次風(fēng)量增大，使得著火位置向一次風(fēng)靠近；一次風(fēng)增大時(shí)爐膛的高溫區(qū)域向下游延遲，爐膛出口溫度從799 ℃增大到869 ℃。這是因?yàn)轱L(fēng)速大時(shí)，一次風(fēng)射流能穿透爐膛內(nèi)的大回流區(qū)，使得煤粉與O2的混合被削弱，焦炭燃燒速率降低，主體火焰燃燒推后。

圖9 溫度分布對(duì)比Fig.9 Comparison of temperature distribution

通過(guò)預(yù)燃室出口溫度判斷煤粉在預(yù)燃室的著火情況(表5)。雖然工況2卷吸的二次風(fēng)量很大，預(yù)燃室內(nèi)煤粉燃燒更加完全，但工況2的預(yù)燃室出口溫度卻最小，這是因?yàn)橐淮物L(fēng)速過(guò)大，剛性增強(qiáng)，卷吸回到預(yù)燃室內(nèi)的煤粉量降低，預(yù)燃室內(nèi)釋放的熱量降低，使預(yù)燃室出口溫度降低。3種情況下預(yù)燃室出口平面(圖2)的平均溫度均達(dá)到1 000 ℃以上，該溫度下煤粉可以著火。

表5 各工況計(jì)算結(jié)果

4.3 各組分濃度分析

爐膛內(nèi)的組分分布能夠直觀反映膛內(nèi)的燃燒和NO生成情況，通過(guò)分析組分濃度能夠?qū)t膛內(nèi)的燃燒情況做出合理預(yù)測(cè)。

中心截面CO濃度分布如圖10所示?？芍S著一次風(fēng)速的增加，CO的分布面積在預(yù)燃室內(nèi)減小，在爐膛內(nèi)增加。這是因?yàn)橐淮物L(fēng)增大能卷吸更多內(nèi)二次風(fēng)進(jìn)入預(yù)燃室，使煤粉燃燒完全，所以預(yù)燃室內(nèi)CO濃度下降，而一次風(fēng)速大會(huì)破壞爐膛內(nèi)的回流區(qū)，與二次風(fēng)的混合效果減弱，爐膛中心的煤粉不完全燃燒生成CO，由于焦炭的燃燒速率降低，CO在爐膛中分布面積增大，煤粉整體燃盡效果變差。表5統(tǒng)計(jì)了爐膛出口焦炭的燃盡程度，可見(jiàn)增大一次風(fēng)速后，焦炭轉(zhuǎn)化率下降明顯，煤粉燃盡效果不佳。

圖10 中心截面CO濃度分布Fig.10 CO concentration distribution in central section

煤粉揮發(fā)分(Vol)析出位置的氣氛對(duì)控制NO的排放非常重要。在氧化性氣氛下，揮發(fā)分氮極易被O2氧化成NO，導(dǎo)致較高的NOx排放；而在還原性氣氛下，揮發(fā)分中的NO前驅(qū)物如HCN及NH3將更傾向于被還原成N2，揮發(fā)分產(chǎn)生的大量CH自由基也能在還原性氣氛及焦炭的作用下將已經(jīng)生成的NO還原成N2。

中心截面的Vol濃度分布如圖11所示。與CO分布相似，一次風(fēng)與二次風(fēng)的混合變差使得爐膛內(nèi)Vol分布面積隨一次風(fēng)速的增大而增大。結(jié)合圖10分析，該燃燒器出口形成的流場(chǎng)能夠保證Vol大量析出的位置處于CO包圍的還原性氣氛中，可以抑制揮發(fā)分氮向NO轉(zhuǎn)移，同時(shí)將預(yù)燃室及爐膛大回流區(qū)內(nèi)卷吸煙氣中的NO還原成N2。

圖11 中心截面的Vol濃度分布Fig.11 Vol concentration distribution of the central section

NOx生成速率分布如圖12所示，出口截面NOx濃度對(duì)比見(jiàn)表5。隨著一次風(fēng)速?gòu)?7 m/s增加到20 m/s，出口截面NO濃度提高約10%，這是由于一次風(fēng)速過(guò)大卷吸了大量?jī)?nèi)二次風(fēng)，預(yù)燃室內(nèi)還原性氣氛被破壞，且大的一次風(fēng)速破壞了爐膛內(nèi)的大回流區(qū)，使其無(wú)法卷吸并還原煙氣中的NO，導(dǎo)致出口截面NOx排放偏高。一次風(fēng)速0從17 m/s降到14 m/s時(shí)出口截面NO濃度也略有升高，主要原因是一次風(fēng)速過(guò)低，一、二次風(fēng)在爐膛后期混合減弱，卷吸的煙氣量不足，還原的NO量較少。

圖12 中心截面的NO生成速率分布Fig.12 Distribution of NO production rate in the central section

綜上，14～17 m/s存在最優(yōu)的一次風(fēng)速，使得爐膛出口的NOx排放最低。

5 結(jié) 論

1)一次風(fēng)速過(guò)低時(shí)風(fēng)速衰減過(guò)快，一、二次風(fēng)后期混合減弱，不利于煤粉穩(wěn)定著火，同時(shí)使得爐膛中心卷吸的煙氣量變少，未參與到中心回流區(qū)還原其中NO的煙氣量變多，造成NO排放量升高。

2)適當(dāng)增大一次風(fēng)速可增大預(yù)燃室環(huán)形回流區(qū)的回流量，有利于煤粉著火和穩(wěn)定燃燒。

3)一次風(fēng)速過(guò)大時(shí)會(huì)破壞爐膛內(nèi)大回流區(qū)，使得二次風(fēng)與煤粉的混合被削弱，不利于爐膛內(nèi)煤粉穩(wěn)定著火；同時(shí)爐膛內(nèi)焦炭氧化速率降低，火焰被拉長(zhǎng)，煤粉整體燃盡效果差，焦炭轉(zhuǎn)化率下降明顯；同時(shí)一次風(fēng)速過(guò)大則無(wú)法卷吸并還原爐膛煙氣中的NO，將導(dǎo)致出口截面NOx排放偏高。一次風(fēng)速?gòu)?7 m/s增加到20 m/s，出口截面NO濃度提高約10%。因此本文研究的預(yù)燃式低氮燃燒器的最佳一次風(fēng)速在14～17 m/s。