循環(huán)流化床是一種商業(yè)化運(yùn)行較為成熟的煤清潔燃燒枝術(shù)，具有燃料適應(yīng)性強(qiáng)、負(fù)荷調(diào)節(jié)范圍大、運(yùn)營費(fèi)用低、燃燒效率高、可實(shí)現(xiàn)爐內(nèi)脫硫脫硝等優(yōu)點(diǎn)，已被廣泛應(yīng)用到工業(yè)中，在構(gòu)建新型電力系統(tǒng)中有望發(fā)揮更大作用。

目前，許多學(xué)者采用多種手段對循環(huán)流化床鍋爐爐內(nèi)流動(dòng)特性進(jìn)行研究，如冷態(tài)試驗(yàn)、數(shù)值計(jì)算等。沈志恒等采用歐拉－歐拉雙流體模型模擬1025t/h循環(huán)流化床爐內(nèi)氣固流動(dòng)行為，耦合了化學(xué)反應(yīng)模型，給出了爐內(nèi)氣固流動(dòng)及其反應(yīng)特性。王坤等基于CPFD方法模擬研究了循環(huán)流化床床內(nèi)燃燒煤與生物質(zhì)顆粒時(shí)顆粒的流動(dòng)特性。

煤在爐內(nèi)燃燒是一個(gè)復(fù)雜的傳熱傳質(zhì)過程，尤其是循環(huán)流化床鍋爐，爐內(nèi)固體燃料受底部上行流化風(fēng)的作用，相比層燃鍋爐和煤粉鍋爐其運(yùn)動(dòng)軌跡隨機(jī)性較大，這就對循環(huán)流化床鍋爐爐內(nèi)顆粒運(yùn)動(dòng)和燃燒過程的研究造成了極大的挑戰(zhàn)。循環(huán)流化床鍋爐爐內(nèi)顆粒的運(yùn)動(dòng)和燃燒過程不但影響著鍋爐的運(yùn)行效率，同時(shí)還關(guān)系到鍋爐的安全運(yùn)行。為此有必要將理論與實(shí)際相結(jié)合去分析流化床鍋爐爐內(nèi)顆粒的流動(dòng)狀態(tài)。

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

送入循環(huán)流化床鍋爐的燃料需提前經(jīng)中速或高速磨煤機(jī)進(jìn)行預(yù)處理后達(dá)到設(shè)計(jì)所要求的粒徑后送入爐膛，燃料顆粒在爐內(nèi)床層處流化燃燒，當(dāng)其自身重力小于流化風(fēng)給予的浮力時(shí)，燃料顆粒將向爐膛上方邊燃燒邊運(yùn)動(dòng)。上行的流化風(fēng)作用于顆粒形成的浮力隨著高度的增加而減小，當(dāng)燃燒受到阻斷使得其重力大于浮力時(shí)顆粒將落回床層受熱并重新被點(diǎn)燃，而隨著燃燒進(jìn)程的推進(jìn)，顆粒中的可燃物與氧氣反應(yīng)使得其質(zhì)量變輕，顆粒本身所受浮力大于重力時(shí)將隨高溫?zé)煔鈴臓t膛出口逃逸進(jìn)入旋風(fēng)分離器，質(zhì)量相對較輕的顆粒將隨煙氣從分離器出口逃逸，而顆粒質(zhì)量相對較重的將在分離器底部的返料口回到爐膛，重新參與爐膛內(nèi)的燃燒過程，直到可燃物燃盡，質(zhì)量變輕。這與循環(huán)流化床鍋爐設(shè)計(jì)的循環(huán)倍率有關(guān)。顆粒燃燒過程是一個(gè)縮核反應(yīng)，其外徑將隨著燃燒進(jìn)程的加劇而縮小，同時(shí)送入流化床內(nèi)的燃料粒徑是經(jīng)篩網(wǎng)進(jìn)行初選，其粒徑必定大小不一，對于煤顆粒燃燒過程需考慮送入爐膛的煤以及燃燒過程中煤顆粒粒徑的變化。為此，對于爐膛內(nèi)氣固兩相流動(dòng)采用歐拉－歐拉雙流體模型，在此基礎(chǔ)上建立顆粒動(dòng)理學(xué)－多組分反應(yīng)模型。

1.1 基本控制方程

本文是對于氣固兩相的動(dòng)量方程、質(zhì)量方程和能量方程進(jìn)行求解，在求解過程中將氣固兩相以混合物的方式進(jìn)行處理，在氣固兩相輸運(yùn)過程中以質(zhì)量分?jǐn)?shù)的方式考慮氣固兩相間的摩擦，以能量耗散的方式計(jì)算其曳力，得到氣固兩相整體的流動(dòng)狀態(tài)。具體計(jì)算公式如下:

1.1.1 質(zhì)量守恒方程

式中:下標(biāo)g表示氣相;s表示固相;i表示顆粒相，i=1，2;v為速度;S為質(zhì)量源項(xiàng);ε為空隙率;ρ為密度。

1.1.2 氣相動(dòng)量守恒方程

式中:βgsi為氣相–顆粒相i間曳力系數(shù);g為當(dāng)?shù)刂亓铀俣?pg為壓力;τg為氣相應(yīng)力張量。

1.1.3 氣相動(dòng)量守恒方程

式中，τsi為不同顆粒相的應(yīng)力張量。

1.1.4 能量方程

氣固兩相能量守恒方程分別為:

式中，λ為氣相或固相的導(dǎo)熱系數(shù);H為焓;hgsi和hsgi為相間換熱系數(shù);SH為非均相化學(xué)反應(yīng)熱源;T為溫度

2 計(jì)算模型及主要參數(shù)

本文以SHX116－1．6/150/90－AⅡ循環(huán)流化床熱水鍋爐為研究對象建立三維幾何模型。由于爐內(nèi)燃燒過程的復(fù)雜性，基于現(xiàn)有計(jì)算資源做適當(dāng)簡化，為多方面地反映爐膛內(nèi)氣固兩相流的具體形態(tài)，在計(jì)算過程中將旋風(fēng)分離器予以剔除，只保留爐膛主體部分，但為貼近循環(huán)流化床的正常運(yùn)行，按照設(shè)計(jì)說明書中給出的循環(huán)倍率計(jì)算旋風(fēng)分離器的返料量并采用UDF的方式設(shè)置回料。模擬對象116MW循環(huán)流化床鍋爐結(jié)構(gòu)簡圖與爐膛模型如圖1所示。

鍋爐設(shè)計(jì)容量為116MW，設(shè)計(jì)煤種AⅡ，出口介質(zhì)溫度130℃，進(jìn)口介質(zhì)溫度70℃，設(shè)計(jì)額定壓力1．6MPa，設(shè)計(jì)介質(zhì)循環(huán)量1652000kg/h，燃料消耗量23377．3kg/h，爐膛高25．7m(Z向)，寬4．43m(X向)，深8．83m(Y向)。一次風(fēng)入口位于爐膛底部布風(fēng)板處，布風(fēng)板面積為2．2×8．64m2，考慮到風(fēng)帽結(jié)構(gòu)較為復(fù)雜，將布風(fēng)板上的進(jìn)風(fēng)口設(shè)置為小孔，孔徑為100mm，爐膛腿部布置了二次風(fēng)和給料口以及回料口，為得到爐膛底部氣固兩相流流動(dòng)過程更多的細(xì)節(jié)，將該區(qū)域網(wǎng)格進(jìn)行了加密處理，采用Gambit軟件對該爐型進(jìn)行了網(wǎng)格劃分，總共劃分網(wǎng)格數(shù)為801554。

循環(huán)流化床在運(yùn)行過程中很難保證其燃燒的顆粒粒徑一致，結(jié)合設(shè)計(jì)煤種粒度的要求，采用雙組分的形式進(jìn)行模擬，其中大顆粒粒徑設(shè)置為1mm，細(xì)顆粒粒徑為0．1mm。相對于設(shè)計(jì)數(shù)據(jù)中粒徑的*大值和*小值。具體模擬的參數(shù)設(shè)置如表1所示。

計(jì)算的初始和邊界條件設(shè)置如下:

(1)進(jìn)出口邊界條件:送風(fēng)和給煤均采用速度入口，出口處設(shè)置為壓力出口。

(2)壁面邊界條件:氣相采用無滑移壁面，顆粒相與壁面之間考慮碰撞和摩擦，選用Johnson和Jackson模型。

3  結(jié)果與討論

3.1 顆粒濃度

圖2為循環(huán)流化床爐內(nèi)細(xì)顆粒濃度分布云圖，圖3展示的是循環(huán)流化床爐內(nèi)細(xì)顆粒沿高度方向顆粒濃度分布圖。從圖中可以明顯地發(fā)現(xiàn)細(xì)顆粒較多地分布于爐膛底部，受上行的流化風(fēng)影響，顆粒隨流化風(fēng)上行，隨高度的增加顆粒濃度不斷降低。由于給料口以及返料口持續(xù)地向爐內(nèi)送入小顆粒，致使在該水平截面上顆粒濃度有個(gè)較大的突變，但隨高度的增加小顆粒的濃度基本趨于穩(wěn)定。從圖3中可以發(fā)現(xiàn)小顆粒分布的細(xì)節(jié)，即受一次風(fēng)影響，貼壁處風(fēng)壓較小，顆粒多聚集于壁面附近，而在爐膛上部，由于壁面附近顆粒的環(huán)—核流動(dòng)，受貼壁處下行風(fēng)的影響，爐膛中心顆粒濃度較高而近壁處顆粒濃度較低。從圖3可以明確地看出二次風(fēng)對顆粒濃度沿高度影響不明顯，但在顆粒濃度云圖中可以發(fā)現(xiàn)二次風(fēng)對細(xì)顆粒起到了攪拌作用，對其顆粒濃度的影響主要表征在其水平截面上顆粒的濃度分布狀況。

圖4為循環(huán)流化床爐內(nèi)粗顆粒濃度分布云圖，圖5展示的是循環(huán)流化床爐內(nèi)粗顆粒沿高度方向顆粒濃度分布圖。其分布形態(tài)與小顆粒基本一致，主要區(qū)別為大顆粒為爐膛內(nèi)預(yù)設(shè)的，無其它進(jìn)料口，為此其隨爐膛高度的增加濃度降低。與小顆粒相比大顆粒質(zhì)量較大，與小顆粒濃度分布圖相比可以發(fā)現(xiàn)更多的大顆粒聚集于爐膛底部，同時(shí)爐膛上方大顆粒濃度明顯小于小顆粒。與小顆粒相似，二次風(fēng)的射入既起到了對爐內(nèi)物料的摻混，同時(shí)其橫向風(fēng)也間接地阻斷了顆粒上行的動(dòng)力，從圖3和圖5可以明顯地發(fā)現(xiàn)該處顆粒濃度的曲線斜率梯度較大。

3.2  速度分布

圖6為爐內(nèi)氣相速度云圖和矢量圖，從圖中可以明顯地看到在一次風(fēng)、二次風(fēng)以及播煤風(fēng)的聯(lián)合作用下爐膛底部氣相運(yùn)動(dòng)較為激烈，其運(yùn)動(dòng)狀態(tài)呈現(xiàn)不規(guī)則狀，而在Z=10m以上爐膛內(nèi)氣相速度方向相對一致均為向上運(yùn)動(dòng)，受近壁風(fēng)壓作用，貼壁處下行的氣相速度較大，而爐膛中心氣相速度相對較小。在爐膛出口處由于爐膛形狀發(fā)生了突變，氣相運(yùn)動(dòng)軌跡亦發(fā)生變化，在該處爐膛中心區(qū)域氣相速度相對紊亂，氣相將在設(shè)置在爐膛兩側(cè)的出口逃逸，從而導(dǎo)致其在出口處中心部分氣相速度較大。

圖7為小顆粒相速度云圖和矢量圖。從圖中可以看出小顆粒速度矢量圖和云圖與氣相的基本一致。其值略小于氣相速度，這是由于小顆粒是在氣相夾帶下運(yùn)動(dòng)。從圖中可以明顯地發(fā)現(xiàn)二次風(fēng)對顆粒運(yùn)行速度影響較大，在二次風(fēng)區(qū)域顆粒運(yùn)動(dòng)速度略有減低。

圖8為不同高度上顆粒速度沿爐膛寬度以及深度分布圖。由于爐膛沿寬度方向?qū)ΨQ，深度方向上除爐膛出口其他基本相同，因此從原則上來說顆粒速度應(yīng)當(dāng)為對稱型。但由于統(tǒng)計(jì)時(shí)間較短，不同高度上顆粒速度基本呈現(xiàn)對稱形態(tài)。本文分別截取了Z=6m(流化床段)，Z=10m和18m(提升管段)以及Z=25m(爐膛出口段)4個(gè)水平截面上顆粒速度分布情況圖。從圖中可以發(fā)現(xiàn)在流化床段顆粒速度波動(dòng)較大，尤其是爐膛中心區(qū)域顆粒運(yùn)動(dòng)速度較大，這是由于貼壁處氣壓較低致使大量的顆粒在此處聚集，從而導(dǎo)致更多的流化風(fēng)從爐膛中心區(qū)域穿過流化床，氣流速度和剛度較大，從而帶動(dòng)顆粒速度變大。

提升管段Z=10m處沿爐膛深度方向上明顯呈現(xiàn)中間顆粒速度高而兩邊速度低的形態(tài)，而沿爐膛寬度方向，此處爐膛呈現(xiàn)漸擴(kuò)形態(tài)，爐膛橫截面積逐漸變大從而導(dǎo)致氣相中心氣壓變高，貼壁處氣壓低，從而使得此處氣相下行的流速大，亦帶動(dòng)顆粒速度變快，而中心區(qū)域速度較低。當(dāng)進(jìn)入直管段Z=18m時(shí)，氣流已經(jīng)過了較長時(shí)間的重整，氣流上行的形態(tài)相對平穩(wěn)，中心氣流流速較大，從而帶動(dòng)顆粒速度形成中心高、兩邊低的形態(tài)。

在爐膛出口處，受出口形態(tài)的影響，以及顆粒運(yùn)行狀態(tài)已趨于穩(wěn)定，使得顆粒速度沿深度方向上除貼壁處受貼壁風(fēng)影響外，其他區(qū)域顆粒速度基本相同。在寬度方向上，由于其長度為寬度的一半，其形態(tài)亦呈現(xiàn)中間高四周低，這與提升管處基本一致。

4 結(jié)論

(1)粗細(xì)顆粒在垂直方向上均呈現(xiàn)隨著高度的提升濃度降低的形態(tài)，無論是粗顆粒還是細(xì)顆粒在爐膛高度為5m以下時(shí)濃度較高，而在流化床上部爐膛濃度較低，尤其是粗顆粒，顆粒濃度低于細(xì)顆粒1～2個(gè)數(shù)量級，這與其設(shè)計(jì)初衷基本一致，即粗顆粒起流化作用而小顆粒參與上部爐膛燃燒過程。

(2)爐內(nèi)氣相速度受爐膛結(jié)構(gòu)的影響較大，爐膛底部區(qū)域橫截面積小但受顆粒阻力的作用，需要一定的風(fēng)室風(fēng)壓，在上升過程中隨著截面積的增大，風(fēng)速降低，受顆粒流動(dòng)的影響，將在爐墻附近形成下行的貼壁風(fēng)，使得顆粒在爐膛內(nèi)部呈現(xiàn)環(huán)－核流動(dòng)。

(3)為了進(jìn)一步分析爐內(nèi)氣固兩相的流動(dòng)特性，考察了沿高度方向上不同截面氣固兩相的速度分布，發(fā)現(xiàn)隨著高度的增加，氣固兩相的速度分布更趨于一致。