空氣分級燃燒下NOx生成特性的研究

孫保民，王頂輝，段二朋，郭永紅，曹為華，張守恒

（1.華北電力大學(xué) 電站設(shè)備狀態(tài)監(jiān)測與控制教育部重點實驗室，北京102206；2.綏中發(fā)電有限責(zé)任公司，葫蘆島125222）

氮氧化物（NOx）是主要的大氣污染物之一，易造成酸雨和光化學(xué)煙霧，對環(huán)境和人體產(chǎn)生很大的危害.在我國，NOx主要來源于化石燃料的燃燒，尤其是火電廠生產(chǎn)過程.根據(jù)最新的《火電廠大氣污染排放指標(biāo)》，燃煤電廠需采取有效的措施來實現(xiàn)NOx的減排.空氣分級燃燒技術(shù)在一定程度上實現(xiàn)了NOx的減排，與SCR 煙氣脫硝技術(shù)相比，空氣分級燃燒技術(shù)的成本較低、操作簡單，在我國電站鍋爐中應(yīng)用廣泛.對空氣分級燃燒下NOx的生成特性進行研究有助于進一步了解其降低NOx排放量的機理，對完善空氣分級燃燒技術(shù)具有重要的指導(dǎo)意義.

國內(nèi)外許多學(xué)者利用數(shù)值模擬方法對不同燃盡風(fēng)率下的空氣分級燃燒進行了研究［1-2］，但這些研究主要側(cè)重于權(quán)衡燃盡風(fēng)率變化產(chǎn)生的正面和負面影響，進而選擇較合適的燃盡風(fēng)率，NOx生成特性分析也是基于NOx質(zhì)量濃度分布的，然而爐內(nèi)某處NOx質(zhì)量濃度明顯受到上游累積作用的影響，因此這種分析方法并不夠直觀.近年來，部分學(xué)者對NOx反應(yīng)速率進行了分析［3-5］，相對而言，這種方法更為科學(xué)，但這些研究只是通過NOx生成速率來反映NOx的生成區(qū)域，并沒有深入分析NOx生成特性.筆者通過爐內(nèi)NOx的反應(yīng)速率考察了各區(qū)域NOx生成與還原的快慢，比較不同燃盡風(fēng)率下NOx的反應(yīng)速率，結(jié)合其他參數(shù)的變化趨勢，分析了爐內(nèi)各區(qū)域NOx的生成特性和主要控制因素，為利用空氣分級燃燒技術(shù)更好地控制NOx排放量提供參考.

1 研究對象與計算方法

1.1 研究對象

以俄羅斯塔干羅格鍋爐廠生產(chǎn)的800 MW 超臨界一次中間再熱直流鍋爐為研究對象.該鍋爐為矩形單爐膛T 型布置，鋼架全懸吊結(jié)構(gòu)，采用平衡通風(fēng)和固態(tài)排渣.爐膛橫截面尺寸為：長30.986m，寬15.472m，頂棚管標(biāo)高為81m.燃燒方式為旋流燃燒器前后墻對沖燃燒，在爐膛標(biāo)高22.9～35.8m之間布置4層旋流燃燒器，各層間距為4.3m，單側(cè)每層布置6只，全爐共48只，配備8套直吹式制粉系統(tǒng).為解決NOx排放質(zhì)量濃度過大的問題，利用空氣分級燃燒技術(shù)對該鍋爐進行改造，即在爐膛標(biāo)高43.5m 和44.5m 處增設(shè)2層燃盡風(fēng)噴口，單側(cè)每層布置6個燃盡風(fēng)噴口，共24個.下層燃盡風(fēng)噴口與最上層一次風(fēng)噴口間距為7.7m.

所用煤種為該鍋爐實際燃用煤種，燃煤量和空氣量均與實際運行數(shù)據(jù)相吻合.當(dāng)燃盡風(fēng)率發(fā)生變化時，主燃區(qū)二次風(fēng)風(fēng)量發(fā)生相應(yīng)變化，對燃盡風(fēng)率分別為15%、20%和25%的工況進行計算.

1.2 計算區(qū)域及網(wǎng)格劃分

對該鍋爐進行全尺寸三維模擬，將鍋爐冷灰斗底部至水平煙道入口之間的區(qū)域作為計算區(qū)域，利用Gambit軟件對爐膛整體和燃燒器附近區(qū)域進行網(wǎng)格劃分（見圖1），均采用六面體結(jié)構(gòu)化網(wǎng)格，生成大約70萬個網(wǎng)格.根據(jù)爐內(nèi)燃燒特點合理調(diào)整網(wǎng)格的疏密，即對流場變化劇烈的燃燒器和燃盡風(fēng)噴口附近區(qū)域進行局部網(wǎng)格加密.同時根據(jù)旋流燃燒器出口流場的特點，將每個燃燒器附近的網(wǎng)格劃分為與其同軸的局部輻射狀網(wǎng)格，這樣可使流線與網(wǎng)格線盡可能吻合，有效地避免了偽擴散［6］.

圖1 爐膛網(wǎng)格劃分Fig.1 Grid division of the furnace

1.3 數(shù)學(xué)模型

1.3.1 流動及燃燒模型

爐內(nèi)湍流流動采用可實現(xiàn)雙方程模型（realizablek-εmodel）計算［7］，近壁區(qū)域采用標(biāo)準(zhǔn)壁面函數(shù)（standard wall function）［8］，輻射傳熱采用P1輻射模型，采用隨機軌道模型來跟蹤煤粉顆粒的運動，煤粉燃燒過程中的化學(xué)反應(yīng)與各組分的輸運采用非預(yù)混燃燒模型，氣相湍流燃燒采用混合分數(shù)-概率密度函數(shù)法（mixture-reaction／PDF），揮發(fā)分的釋放采用單步反應(yīng)模型（single-rate model），焦炭燃燒采用動力／擴散控制燃燒模型（kinetics／diffusion-limited char combustion model）.

1.3.2 NOx生成模型

與其他組分相比，爐內(nèi)NOx的質(zhì)量濃度很小，可認為其對爐內(nèi)燃燒和其他組分分布的影響十分微弱，因此NOx的計算采用后處理方法［9］.煤粉燃燒過程中生成的NOx分為快速型NOx、熱力型NOx和燃料型NOx.由于快速型NOx生成量極少，計算中只考慮熱力型NOx和燃料型NOx.湍流脈動對NOx生成的影響則通過溫度和氧濃度的β型概率密度函數(shù)來考慮.

熱力型NOx的生成遵循擴展的Zeldovich鏈鎖反應(yīng)機制，O 與OH 自由基物質(zhì)的量濃度由部分平衡法計算得到［10］.熱力型NOx的生成速率為

燃料型NOx的N 來自于揮發(fā)分或焦炭.揮發(fā)分N 的轉(zhuǎn)化采用De＇Soete模型［11-12］，即揮發(fā)分N 先轉(zhuǎn)化為中間產(chǎn)物HCN，HCN 既能與O2反應(yīng)生成NOx，又能將NOx還原為N2.焦炭N 則直接轉(zhuǎn)化為NOx，這3個反應(yīng)的反應(yīng)速率分別為

式（1）～式（4）中：cN2、cO2、cHCN和cNO分別為N2、O2、HCN 和NO 的物質(zhì)的量濃度；MNO為NO 的相對分子質(zhì)量；MN為N 的相對原子質(zhì)量；R為通用氣體常數(shù)；T為溫度，K；b隨O2濃度的變化而變化；Rc為焦炭的燃燒速率；wNC為焦炭中N 的質(zhì)量分數(shù)；η為焦炭中N 向NO 轉(zhuǎn)化的頻率因子.

1.4 數(shù)學(xué)模型與計算方法的驗證

數(shù)學(xué)模型與計算方法選用是否恰當(dāng)直接影響到計算結(jié)果的準(zhǔn)確性.筆者測量了該鍋爐近期滿負荷運行（未改造）工況下爐膛出口的相關(guān)參數(shù)，對相應(yīng)的工況進行數(shù)值計算，并將試驗結(jié)果與數(shù)值計算結(jié)果進行比較以驗證所用模型與方法的正確性.表1給出了試驗結(jié)果與模擬結(jié)果的對比，其中CO 體積分數(shù)、O2體積分數(shù)和NOx質(zhì)量濃度、試驗值與計算值的相對誤差分別為8.1%、1.2%和4.7%.由表1可以看出，除CO 體積分數(shù)外，其他2個參數(shù)的相對誤差均不超過5%，這是因為CO 體積分數(shù)很小，數(shù)值容易波動，其相對誤差是可以接受的.因此本文所選的計算模型比較恰當(dāng)，利用相同的方法對改造工況進行計算，得到的相關(guān)數(shù)據(jù)有很大的參考價值.

2 計算結(jié)果及分析

2.1 NOx 反應(yīng)速率分布

2.1.1 燃燒器附近NOx的反應(yīng)速率

為了解燃燒初期NOx的變化，對燃燒器出口附近的NOx反應(yīng)速率（見圖2）進行了分析.由圖2可知，NOx主要在燃燒初期產(chǎn)生.煤粉從燃燒器噴出后被高溫?zé)煔饧訜?，煤粉中的部分N 隨揮發(fā)分析出而釋放，轉(zhuǎn)變?yōu)镠CN 等中間產(chǎn)物，在燃燒初期O2相對充足時被氧化為NOx.此外，隨著燃燒的進行，焦炭中的N 會轉(zhuǎn)化為NOx.

表1 試驗結(jié)果與模擬結(jié)果的對比Tab.1 Comparison between actual measurements and calculated results

圖2 燃燒器附近NOx 的反應(yīng)速率Fig.2 Distribution of NOxformation rate around burner

由圖2還可以看到，NOx生成區(qū)域首先隨著二次風(fēng)氣流的旋轉(zhuǎn)而徑向擴張，燃燒過程中逐漸消耗O2，NOx生成速率減小，NOx生成區(qū)域沿氣流方向逐漸縮小.而此區(qū)域中心處NOx的生成速率提前減小，即NOx生成速率形成外圍高、中心低的分布情況.這是因為：一方面外圍二次風(fēng)較多，供氧比較充足，促進燃燒的同時也促進了NOx的生成；另一方面旋轉(zhuǎn)的二次風(fēng)使得大量高溫?zé)煔庠谥行幕亓?，此處O2體積分數(shù)較低，在高溫低氧條件下NOx易被還原.

2.1.2 NOx反應(yīng)速率沿爐膛高度的變化

圖3給出了當(dāng)燃盡風(fēng)率為20%時燃料型NOx反應(yīng)速率、熱力型NOx反應(yīng)速率、O2體積分數(shù)和溫度沿爐膛高度的變化情況.燃料型NOx與熱力型NOx反應(yīng)速率的疊加為NOx的總反應(yīng)速率.由圖3可知，燃料型NOx的反應(yīng)速率明顯大于熱力型NOx，說明燃煤鍋爐中生成的NOx主要來源于燃料型NOx.主燃區(qū)的NOx既有燃料型NOx也有熱力型NOx，而燃盡區(qū)的NOx幾乎全部來自燃料型NOx.由圖3還可以看到，主燃區(qū)NOx生成速率顯著大于燃盡區(qū)，因此控制主燃區(qū)NOx的生成速率是空氣分級燃燒降低NOx排放量的關(guān)鍵.

圖3 各參數(shù)沿爐膛高度的變化Fig.3 Variation of key parameters along furnace height

（1）燃料型NOx

由圖3可知，由于主燃區(qū)的空氣從不同層燃燒器間隔噴入，劇烈燃燒消耗O2，O2體積分數(shù)出現(xiàn)巨大波動，燃料型NOx的反應(yīng)速率也隨之變化.O2體積分數(shù)越大，燃燒過程中的含氮產(chǎn)物越容易被氧化為NOx，燃料型NOx反應(yīng)速率越大.當(dāng)O2體積分數(shù)很小時，燃料型NOx反應(yīng)速率變?yōu)樨撝?，即NOx被還原.這是由于燃燒過程中產(chǎn)生的HCN 等中間產(chǎn)物在缺氧條件下容易還原NOx.研究表明［13］，缺氧條件下焦炭的燃燒也會促進NOx的異相還原.由圖3可以看到，NOx還原速率最快的位置出現(xiàn)在上兩層燃燒器之間，因為此處不僅O2體積分數(shù)很小，還原性氛圍很強，而且劇烈燃燒釋放出大量還原性產(chǎn)物，有利于NOx的還原.

（2）熱力型NOx

與燃料型NOx相比，主燃區(qū)熱力型NOx的反應(yīng)速率也隨著O2體積分數(shù)的波動而變化，兩者區(qū)別在于燃料型NOx反應(yīng)速率與O2體積分數(shù)的整體變化趨勢一致，均隨燃燒器層高度的增加而減??；而熱力型NOx反應(yīng)速率隨燃燒器層高度的增加而增大，說明熱力型NOx對O2的敏感性較弱.此外，主燃區(qū)和燃盡區(qū)都有燃料型NOx，而熱力型NOx的生成幾乎只發(fā)生在溫度很高的主燃區(qū)，說明較高的溫度是生成熱力型NOx的必要因素.根據(jù)式（1），熱力型NOx的生成速率與溫度呈指數(shù)關(guān)系，即當(dāng)溫度較高時，溫度對熱力型NOx的影響比較明顯，而在溫度相對較低的還原區(qū)幾乎不生成熱力型NOx.

2.2 燃盡風(fēng)率對NOx 生成特性的影響

2.2.1 燃盡風(fēng)率對NOx反應(yīng)速率的影響

圖4～圖6分別給出了不同燃盡風(fēng)率下NOx反應(yīng)速率、截面平均O2體積分數(shù)（即O2體積分數(shù)）和焦炭燃燒速率沿爐膛高度的變化.由圖4～圖6可以看到，各圖中3條曲線的變化趨勢相似.燃盡風(fēng)率增大，主燃區(qū)NOx反應(yīng)速率、O2體積分數(shù)和焦炭燃燒速率均有所減小，而在燃盡區(qū)均增大.不同的燃盡風(fēng)率導(dǎo)致燃燒過程在主燃區(qū)和燃盡區(qū)的分配情況不同，NOx反應(yīng)速率也隨之變化.

不同燃盡風(fēng)率下主燃區(qū)焦炭燃燒速率均較大，即主燃區(qū)的燃燒比較劇烈，釋放的含氮產(chǎn)物較多，因此限制NOx生成量的主要因素不是氮的釋放過程，而是其轉(zhuǎn)化過程.在氮的轉(zhuǎn)化過程中，O2體積分數(shù)為主要影響因素.因為在O2體積分數(shù)較大時，燃燒釋放的HCN 等中間產(chǎn)物易被氧化為NOx，而在O2體積分數(shù)很小時則易還原NOx.由燃盡風(fēng)率變化引起的主燃區(qū)O2體積分數(shù)的差異對NOx反應(yīng)速率產(chǎn)生了明顯的影響.對比圖4和圖5中主燃區(qū)的部分曲線，NOx反應(yīng)速率與O2體積分數(shù)的變化趨勢十分一致，說明O2體積分數(shù)是控制主燃區(qū)NOx反應(yīng)速率的關(guān)鍵因素.

圖4 不同燃盡風(fēng)率下NOx 反應(yīng)速率沿爐膛高度的變化Fig.4 Distribution of NOxformation rate along furnace height at different over-fire air rates

圖5 不同燃盡風(fēng)率下O2 體積分數(shù)沿爐膛高度的變化Fig.5 Distribution of O2volume fraction along furnace height at different over-fire air rates

圖6 不同燃盡風(fēng)率下焦炭燃燒速率沿爐膛高度的變化Fig.6 Distribution of coke burn-out rate along furnace heightat different over-fire air rates

為了保證主燃區(qū)未燃盡煤粉在燃盡區(qū)完全燃燒，燃盡區(qū)O2體積分數(shù)均較大，煤粉在O2充足且混合較充分的條件下燃燒，此時限制NOx生成量的主要因素不是O2體積分數(shù).根據(jù)前文分析，燃盡區(qū)生成的NOx幾乎全部來自燃料型NOx，限制NOx生成的主要因素則是燃料燃燒釋放的氮.對比圖4和圖6中燃盡區(qū)的曲線，燃盡區(qū)NOx反應(yīng)速率曲線與焦炭燃燒速率曲線的變化十分吻合.通過改變風(fēng)量分配，空氣分級燃燒技術(shù)將一部分煤粉的燃盡過程由主燃區(qū)轉(zhuǎn)移到燃盡區(qū).燃盡風(fēng)率越大，燃盡區(qū)燃燒的煤粉越多，燃燒速率越大，隨之釋放的氮也就越多，因此NOx的反應(yīng)速率越大.

2.2.2 燃盡風(fēng)率對NOx生成區(qū)域和還原區(qū)域的影響

要研究燃盡風(fēng)率對NOx生成區(qū)域和還原區(qū)域的影響，需要了解爐內(nèi)具體位置NOx的反應(yīng)速率.圖7為某列燃燒器所在縱截面上NOx反應(yīng)速率分布云圖.為了方便分析，只畫出了3 條典型的NOx反應(yīng)速率等值線，其中反應(yīng)速率為0 的等值線為NOx生成區(qū)域和還原區(qū)域的分界線.由于NOx的生成速率顯著大于還原速率，因此將生成速率大于1×10-3mol／（m3·s）的區(qū)域作為NOx的主要生成區(qū)域，而將還原速率大于1×10-4mol／（m3·s）的區(qū)域作為NOx的主要還原區(qū)域.

圖7 不同燃盡風(fēng)率下縱截面上NOx 反應(yīng)速率的變化Fig.7 Distribution of NOxformation rate on vertical section at different over-fire air rates

由圖7可知，不同燃盡風(fēng)率下NOx的主要生成區(qū)域和還原區(qū)域在該截面上的分布位置基本一致.NOx的主要生成區(qū)域位于燃燒器出口附近.在燃盡區(qū)隨著“后期補燃”的進行，燃盡風(fēng)噴口高度附近的爐膛中部出現(xiàn)較小的NOx生成區(qū)域，由于溫度較低，此區(qū)域內(nèi)NOx生成速率相對較小.而NOx的主要還原區(qū)域則出現(xiàn)在主燃區(qū)的爐膛中部和主燃區(qū)與燃盡區(qū)之間，在這些區(qū)域內(nèi)主燃區(qū)送入的O2被消耗殆盡，而燃盡風(fēng)還未送入，較強的還原性氣氛促進了NOx的還原.

由圖7還可以看到，不同燃盡風(fēng)率主要影響了NOx生成區(qū)域與還原區(qū)域的大小.在爐膛下部區(qū)域，隨著燃盡風(fēng)率的增大，NOx生成區(qū)域減小，還原區(qū)域增大.爐膛上部區(qū)域燃盡風(fēng)率對NOx生成區(qū)域的影響十分明顯，當(dāng)燃盡風(fēng)率由15%增至25%時，燃盡區(qū)NOx的生成區(qū)域明顯增大，燃盡風(fēng)噴口高度附近的爐膛中部NOx生成速率明顯增大，這是由于煤粉在燃盡區(qū)燃燒份額增大的緣故.

2.2.3 燃盡風(fēng)率對NOx質(zhì)量濃度分布的影響

圖8給出了不同燃盡風(fēng)率下水平截面NOx平均質(zhì)量濃度（即NOx質(zhì)量濃度）沿爐膛高度的變化.由于NOx主要在主燃區(qū)生成，因此圖8中不同燃盡風(fēng)率下NOx質(zhì)量濃度在主燃區(qū)均急劇增大.燃盡風(fēng)率增大，主燃區(qū)NOx生成速率和生成區(qū)域減小，NOx質(zhì)量濃度明顯減小.主燃區(qū)與燃盡區(qū)之間的NOx被還原，NOx質(zhì)量濃度均明顯減小.燃盡區(qū)會有少量NOx生成，NOx質(zhì)量濃度均略有增大，燃盡風(fēng)率越大，燃盡區(qū)NOx生成速率越大，NOx質(zhì)量濃度增大的幅度也越大，因此3條曲線間的差距減小.燃盡風(fēng)率為15%、20%和25%時，爐膛出口的NOx質(zhì)量 濃 度 分 別 為449.8 mg／m3、422.1 mg／m3和395.7mg／m3，可見燃盡風(fēng)率增大明顯將低了NOx的排放量.

圖8 不同燃盡風(fēng)率下NOx質(zhì)量濃度沿爐膛高度的變化Fig.8 Distribution of NOx mass concentration along furnace height at different over-fire air rates

3 結(jié) 論

（1）NOx主要在燃燒初期生成，但當(dāng)燃料與O2混合不充分時會發(fā)生NOx的還原反應(yīng).

（2）主燃區(qū)燃料型NOx和熱力型NOx的反應(yīng)速率均隨著O2體積分數(shù)的變化而波動，熱力型NOx對O2體積分數(shù)的敏感性弱于燃料型NOx，而對溫度的敏感性強于燃料型NOx.

（3）燃盡風(fēng)率增大，主燃區(qū)NOx生成速率減小，其主要控制因素為O2體積分數(shù)；而燃盡區(qū)NOx生成速率增大，其主要控制因素為焦炭燃燒速率.

（4）控制主燃區(qū)NOx的生成速率是空氣分級燃燒降低NOx排放的關(guān)鍵.燃盡風(fēng)率增大，主燃區(qū)NOx生成區(qū)域減小，還原區(qū)域增大，NOx排放質(zhì)量濃度減小.

［1］方慶艷，汪華劍，陳剛，等.燃盡風(fēng)對超超臨界鍋爐燃燒特性影響的數(shù)值模擬［J］.華中科技大學(xué)學(xué)報：自然科學(xué)版，2010，38（11）：92-95. FANG Qingyan，WANG Huajian，CHEN Gang，et al.Numerical simulation on effects of over fire air on the combustion characteristics in an ultra-supercritical boiler［J］.Journal of Huazhong University of Science and Technology：Natural Science，2010，38（11）：92-95.

［2］劉敦禹，秦明，孫巧群，等.600 MW 超超臨界墻式切圓鍋爐燃燒過程數(shù)值模擬［J］.熱能動力工程，2011，26（1）：89-93. LIU Dunyu，QIN Ming，SUN Qiaoqun，etal.Numerical simulation of the combustion process of a 600 MW ultra-supercritical wall type tangentially-fired boiler［J］.Journal of Engineering for Thermal Energy and Power，2011，26（1）：89-93.

［3］劉敦禹，秦明，劉輝，等.深度空氣分級條件下爐內(nèi)氮氧化物生成的數(shù)值模擬［J］.燃燒科學(xué)與技術(shù)，2011，17（5）：469-475. LIU Dunyu，QIN Ming，LIU Hui，etal.Numerical simulation of NOxproduction under deep air staged combustion in furnace［J］.Journal of Combustion Science and Technology，2011，17（5）：469-475.

［4］CHOENG R C，CHANG N K.Numerical investigation on the flow，combustion and NOxemission characteristics in a 500 MW tangentially-fired pulverizedcoal boiler［J］.Fuel，2009，88（9）：1720-1731.

［5］白濤，孫保民，郭永紅，等.旋流燃燒鍋爐的NOx排放數(shù)值模擬研究［J］.熱力發(fā)電，2012，41（7）：13-21. BAI Tao，SUN Baomin，GUO Yonghong，etal.Numerical simulation of NOxemission in a boiler with swirl burners［J］.Thermal Power Generation，2012，41（7）：13-21.

［6］蘇勝，蔡興飛，呂宏彪，等.采用雙混合分數(shù)／概率密度函數(shù)方法模擬混煤燃燒［J］.中國電機工程學(xué)報，2012，32（2）：45-52. SU Sheng，CAI Xingfei，Lü Hongbiao，etal.Simulation of coal blended combustion with the two-mixture-fraction／PDF method［J］.Proceedings of the CSEE，2012，32（2）：45-52.

［7］郭永紅，白濤，孫保民，等.降低鍋爐爐渣含碳量問題的研究［J］.動力工程學(xué)報，2012，32（2）：106-111. GUO Yonghong，BAI Tao，SUN Baomin，etal.Study on reducing carbon content in slag of a utility boiler［J］.Journal of Chinese Society of Power Engineering，2012，32（2）：106-111.

［8］王福軍.計算流體動力學(xué)分析［M］.北京：清華大學(xué)出版社，2005：205-206.

［9］RYOICHI K，HISAO M，AKIRA S.Numerical analysis of pulverized coal combustion characteristics using advanced low-NOxburner［J］.Fuel，2004，83（6）：693-703.

［10］劉建全，孫保民，陳曉萌，等.旋流燃燒器給粉方式對鍋爐NO 排放特性影響的數(shù)值模擬［J］.機械工程學(xué)報，2011，47（10）：119-125. LIU Jianquan，SUN Baomin，CHEN Xiaomeng，et al.NO distribution numerical simulation on different way of fuel feeding for swirl combustion boiler［J］.Journal of Mechanical Engineering，2011，47（10）：119-125.

［11］COELHO P J，CARVALHO M G.Mathematical modeling of NO formation in a power station boiler［J］.Combustion Science and Technology，1995，108（4／5／6）：363-382.

［12］THOMAS L B，F(xiàn)RANCISCO C，SEBASTIEN C，etal.Coal combustion modeling of large power plant for NOxabatement［J］.Fuel，2007，86（9）：2213-2220.

［13］周俊虎，劉廣義，劉海峰，等.神華煤燃燒再燃中NOx生成與還原試驗研究［J］.浙江大學(xué)學(xué)報：工學(xué)版，2007，41（3）：499-503. ZHOU Junhu，LIU Guangyi，LIU Haifeng，etal.Experimental study on NOxformation and reduction by Shenhua coal during combustion and reburning［J］.Journal of Zhejiang University：Engineering Science，2007，41（3）：499-503.