王 飛, 潘衛(wèi)國, 王文歡, 韓 濤
(1.上海電力學院, 上海 200090; 2.上海發(fā)電環(huán)保工程技術研究中心, 上?!?00090)
?
墻式切圓鍋爐燃燒過程的數(shù)值模擬
王飛1, 潘衛(wèi)國1, 王文歡1, 韓濤2
(1.上海電力學院, 上海200090; 2.上海發(fā)電環(huán)保工程技術研究中心, 上海200090)
摘要:為研究墻式切圓鍋爐爐內(nèi)的燃燒問題,采用Fluent模擬軟件,分析了爐內(nèi)的速度場、溫度場和NOspan的分布規(guī)律.結果表明,爐膛燃燒器區(qū)橫截面的切圓直徑隨著爐膛高度的增加逐漸增加;在爐膛中心縱截面上,爐內(nèi)中心溫度低于兩側溫度,而且隨著爐膛高度的增加,溫度先增后降,爐壁附近出現(xiàn)局部高溫;NO濃度沿高度方向先升后降.
關鍵詞:超臨界鍋爐; 燃燒; 數(shù)值模擬
目前,國內(nèi)使用的煤粉鍋爐80%以上是四角切圓鍋爐,旋流燃燒鍋爐的使用率不到10%[1].四角切圓鍋爐具有結構簡單、燃料的適應性強、風粉混合均勻等明顯特點,在中國乃至全世界的應用范圍非常廣泛.而在常見的切圓鍋爐中,又存在一種燃燒器布置于4面墻壁上的墻式切圓鍋爐.墻式切圓鍋爐充分利用了爐膛的容積,因而具有避免水冷壁結渣、熱負荷均勻以及爐膛出口溫度偏差小等顯著特點,被廣泛應用在實際的鍋爐運行中.
國內(nèi)外的學者對切圓鍋爐進行了大量研究,采用的是試驗、理論和數(shù)值模擬等方法.張開利等人[2]采用磨煤機系統(tǒng)調(diào)整試驗,優(yōu)化調(diào)整了機組的運行參數(shù);秦明等人[3]通過實驗室冷態(tài)?;瘜嶒?得出了墻式切圓燃燒噴口采用分組布置的形式可以有效減輕氣流貼壁的結論;潘維等人[4]采用數(shù)值模擬的方法進行了切圓鍋爐改造工況的研究;潘維和李彥鵬等人[5-6]針對切圓鍋爐爐膛網(wǎng)格的生成方法進行了研究,發(fā)現(xiàn)pave 方法畫出的輻射狀網(wǎng)格能夠有效地抑制偽擴散;呂太等人[7]對燃燒器噴嘴射入的角度進行了研究,通過Fluent數(shù)值模擬多種不同的射入角度,從而分析對比不同射入角度前后所形成的溫度場、速度場及組分場的情況,分析所選角度的合理性,為同類機組的抑制結焦提供了有效參考.文獻[8]采用一種新網(wǎng)格對609 MW的四角切圓鍋爐爐內(nèi)流動、傳熱、燃燒及NOx排放等情況進行了模擬,并進行實驗驗證,結果表明,模擬結果與實驗相吻合,數(shù)值計算可行、可靠.文獻[9] 對500 MW四角切圓鍋爐爐內(nèi)流動、燃燒、傳熱及NOx生成過程等進行了研究,著重研究了空氣分級燃燒、二次風的配置、煙道氣循環(huán),以及燃料分級燃燒(再燃燒)降低NOx排放等技術.
本文基于Fluent軟件,以某電廠墻式切圓鍋爐為研究對象,對其爐內(nèi)的燃燒過程進行數(shù)值模擬,得出速度場、溫度場和NOx的分布情況及其燃燒規(guī)律,以期為研究該類型鍋爐的安全運行提供理論依據(jù).
1研究對象
本文主要以某電廠HG-1950/25.4-HM15超臨界直流鍋爐為研究對象.該鍋爐為一次中間再熱、超臨界壓力變壓運行的直流鍋爐,橫截面為(爐膛寬度)20 402.3 mm×(深度)20 027.3 mm,采用中速磨直吹式制粉系統(tǒng).
鍋爐采用墻式切圓新型燃燒方式,主燃燒器布置在水冷壁的4面墻上,每層4只燃燒器對應1臺磨煤機,燃燒器分組布置,共有3組,每組由下往上依次是二次風、一次風、二次風、一次風、二次風.SOFA燃燒器布置在主燃燒器區(qū)上方水冷壁的4個角,以實現(xiàn)分級燃燒降低NOx排放.鍋爐的結構如圖1所示.
2網(wǎng)格劃分
本文基于CFD軟件,對鍋爐爐內(nèi)的燃燒過程進行數(shù)值模擬.為研究方便,采用非結構化網(wǎng)格和結構化網(wǎng)格相結合的方法對爐膛進行網(wǎng)格劃分,將爐膛整體劃分為冷灰斗區(qū)、燃燒器區(qū)、上爐膛區(qū)和爐膛出口4大部分.
燃燒器部分是煤粉和空氣噴入的主要部分,會產(chǎn)生劇烈的湍流和燃燒等現(xiàn)象,因此為減少模擬過程出現(xiàn)偽擴散對計算準確性的影響,對燃燒器區(qū)適當?shù)剡M行了網(wǎng)格局部加密,同時采用非結構化網(wǎng)格進行劃分.而其他幾部分結構相對均勻,可采用結構化網(wǎng)格進行劃分.所生成的爐膛網(wǎng)格模型如圖2所示.
3數(shù)學模型及研究方法
鍋爐的爐內(nèi)燃燒過程是湍流流動與化學反應相互耦合的過程,它不僅涉及湍流和多組分的流動,而且相關的化學反應包括氣相和顆粒相的燃燒,因此整個過程非常復雜[10].為簡化整個燃燒過程,本文采用標準k-?雙方程湍流模型[11]模擬氣相湍流的輸送;運用隨機軌道模型(Stochastic Tracking Model)跟蹤煤粉顆粒;采用混合分數(shù)概率密度函數(shù)法(PDF)模擬氣相湍流的燃燒;采用動力-擴散模型(Kinetic/Diffusion Limited Model)模擬焦炭的燃燒[12];運用雙競爭反應模型(The Two Competing Rates Model)進行煤粉揮發(fā)分的釋放;用于采用P-1輻射模型[13]用于輻射傳熱.計算時采用Simple算法[14].
模擬時,一次風、二次風和燃盡風入口的邊界條件設置為速度入口,煙道出口設置為壓力出口邊界.爐膛壁面采用標準壁面方程、無滑移邊界條件,即溫度邊界條件,給定壁面溫度為600 K.運行參數(shù)見表1.
煤粉顆粒采取Rosin-Rammler分布,其最小顆粒為5 μm,最大顆粒為310 μm,平均顆粒為60 μm,分布指數(shù)取為1.5.煤質(zhì)分析見表2.
4計算結果與分析
4.1速度場
圖3為各層一次風口的速度分布圖.從圖3可以看出,各個截面的速度很好地形成了切圓的形式,符合設計要求.煤粉從一次風口進入爐膛以后,在湍流的作用下立即與空氣混合.
同時,隨著爐膛高度的增加,各個截面所形成的切圓直徑逐漸增大.這是因為,隨著爐膛高度的增加,爐內(nèi)旋轉(zhuǎn)強度也逐漸增加,從而使得形成的切圓直徑增大.
4.2溫度場
圖4為不同截面的溫度分布示意.從圖4可以看出,煤粉從一次風口進入爐膛后被迅速點燃,然后在4股射流的相互干擾影響下,形成穩(wěn)定的切圓.這也與上述得出的速度分布情況相一致.煤粉在離開壁面一定的距離時,揮發(fā)分快速析出并開始燃燒,釋放出大量的熱量,從而形成局部高溫區(qū).因此,高溫區(qū)主要發(fā)生在煤粉的射流軌跡上,不在一次風口處,而是離水冷壁壁面有一定的距離,這符合鍋爐的實際運行情況.同時,隨著爐膛高度的增加,環(huán)狀高溫區(qū)不斷增大,這也與爐膛高度增加導致旋流強度逐漸增加以及所形成的切圓直徑不斷增大有關.
圖5為爐膛中心縱截面的溫度分布.
從圖5可以看出,從燃燒器噴出的煤粉立即與空氣混合,燃燒釋放大量的熱,使得爐膛兩側的溫度高于爐膛中心溫度,但高溫區(qū)不在壁面處,也不在爐膛中心處,而是在離水冷壁有一定的距離處.而且隨著爐膛高度的增加,爐內(nèi)溫度先升后降,在燃盡風區(qū)域,隨著溫度較低的燃盡風的大量進入,使得燃燒不充分的焦炭充分燃燒,進而使得該區(qū)域的溫度不斷升高.隨著爐膛高度的進一步增加,出現(xiàn)了折焰角以后,溫度開始降低.
4.3NOx分布
鍋爐燃燒中產(chǎn)生的NOx有3種,分別為燃料型NOx,快速型NOx和熱力型NOx.由于快速型NOx的量極少,因此鍋爐燃燒中所說的NOx指的就是燃料型NOx和熱力型NOx兩種.而所生成的NOx中,90%以上是NO,故以下以NO濃度代表NOx濃度.圖6為z=10中心截面NO濃度分布圖.
由圖6可以看出,高濃度的NO主要分布在燃燒器區(qū)域的兩側靠近水冷壁的壁面處.對于燃料型的NO來說,煤粉顆粒揮發(fā)出的含氮的化合物被氧化成HCN,隨后被氧化成NO,這部分氧氣充足,NO的生成量增加.過了主燃區(qū)后,由于氧氣被大量消耗,氧氣含量大幅下降,部分NO被還原,使得NO濃度也略微開始下降,而燃盡風的加入,使得焦炭開始燃燒,NO的濃度也開始有所回升.過了燃盡風后,此時NO處于還原性氣氛,NO濃度開始大幅下降,最終趨于穩(wěn)定.而對于熱力型的NO來說,它是由空氣中的N2在高溫下直接被氧氣氧化而生成的,這與圖6中的高NO區(qū)域正好相對應.
5結語
本文對墻式切圓鍋爐爐內(nèi)的燃燒過程進行了數(shù)值模擬,較為準確地描述了爐內(nèi)的速度場、溫度場和NOx的分布規(guī)律,此研究結果可以為該類型鍋爐的設計和運行提供一定的理論依據(jù).
參考文獻:
[1]吳治永,李文彥.某新型旋流燃燒鍋爐的數(shù)值模擬[J].現(xiàn)代電力,2004(2):11-16.
[2]張開利,李樹海.印度嘉佳電廠660 MW超臨界鍋爐燃燒優(yōu)化試驗研究[J].鍋爐技術,2015(2):39-43.
[3]秦明,吳少華,郭建明,等.四墻切圓濃淡燃燒器射流特性的實驗研究[J].節(jié)能技術,2015(1):66-68.
[4]潘維,池作和,斯東波,等.200 MW 四角切圓燃燒鍋爐改造工況數(shù)值模擬[J].中國電機工程學報,2005,25(8):110-115.
[5]潘維,池作和,斯東波,等.四角切圓燃燒鍋爐爐膛網(wǎng)格生成方法的研究[J].動力工程,2005,25(3):359-363.
[6]李彥鵬,張強,徐益謙,等.四角切圓爐三維流動模擬中減小偽擴散的新方法[J].能源研究與利用,2001(2):30-32.
[7]呂太,趙凱,許賀.600 MW墻式切圓鍋爐燃燒器數(shù)值模擬優(yōu)化及改造研究[J].電站系統(tǒng)工程,2015(2):1-4.
[8]YIN Chungen,CAILLAT S ébastien,HARION Jean-Luc,etal.Investigation of the flow,combustion,heat-transfer and emissions from a 609 MW utility tangentially fired pulverized-coal boiler[J].Fuel,2002,81(5):997-1006.
[9]CHOI Choeng Ryul,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):1 720-1 731.
[10]EATONETAL A M,SMOOT L D,HILL S C,etal.Components,formulations,solutions,evaluation and application of comprehensive combustion models[J].Progress in Energy and Combustion Science,1999,25(4):387-436.
[11]林海,馬曉茜,余昭勝.大型城市生活垃圾焚燒爐的數(shù)值模擬[J].動力工程學報,2010,30(2):128-132.
[12]呂太,王澤民,吳紅峰,等.600 MW機組鍋爐分層摻燒神華煤的混煤燃燒數(shù)值模擬[J].電站系統(tǒng)工程,2015(1):29-31.
[13]鄧中乙,余永生,謝紅,等.二次風偏轉(zhuǎn)對減輕超超臨界墻式鍋爐高溫腐蝕數(shù)值模擬[J].中國電力,2015(1):47-51.
[14]劉沛奇,韓冰,李壯,等.采用近流線網(wǎng)格對NOx生成過程的數(shù)值研究[J].發(fā)電與空調(diào),2014(6):1-5.
(編輯胡小萍)
Numerical Simulation of Combustion Characteristics in a Wall-tangentially Fired BoilerWANG Fei1, PAN Weiguo1, WANG Wenhuan1, HAN Tao2
(1.ShanghaiUniversityofElectricPower,Shanghai200090,China; 2.ShanghaiTechnologyResearch
CenterinPowerEnvironmentalProtectionEngineering,Shanghai200090,China)
Abstract:The combustion characteristics in the furnace of wall-tangentially fired boiler are numerically simulated by using Fluent.The velocity field,the temperature field and the NOspandistribution are analysised.The results show that on the cross section of the burner area in the furnace,the diameter of the circle increases with the increase of the furnace elevation.In the longitudinal section of the furnace center,the central temperature is lower than that in both sides.The furnace temperature firstly incresaes and then decreases with the increase of the furnace elevation.Also,the local high temperature occurs near the wall.The fraction of NO increases at first,and then decreases with the increase of the furnace elevation.
Key words:supercritical boiler; combustion; numerical simulation
中圖分類號:TK224.11
文獻標志碼:A
文章編號:1006-4729(2016)01-0005-04
基金項目:上海發(fā)電環(huán)保工程技術研究中心項目(11dz2281700); 上海市科學技術委員會科技攻關項目(13dz1202703,14dz1200200).
通訊作者簡介:王飛(1990-),男,在讀碩士,陜西延安人.主要研究方向為鍋爐安全運行.E-mail:sgjlwf@163.com.
收稿日期:2015-07-19
DOI:10.3969/j.issn.1006-4729.2016.01.002