胡 南，徐 夢，張 縵，楊海瑞，張守玉，趙 冰，鞏太義，王家林

(1.長春工程學(xué)院 吉林省水利電力工程物理級仿真與安全科技創(chuàng)新中心，吉林 長春 130012；2.清華大學(xué) 熱科學(xué)與動力工程教育部重點實驗室，北京 100084；3.上海理工大學(xué) 能源與動力工程學(xué)院，上海 200093；4.華電國際電力股份有限公司天津開發(fā)區(qū)分公司，天津 300270)

0 引 言

循環(huán)流化床(circulating fluidized bed，CFB)鍋爐具有燃料適用范圍廣、低成本燃燒中脫硫、低氮氧化物排放等優(yōu)點，是大規(guī)模清潔利用高硫及高灰煤、矸石、中煤、煤泥和生物質(zhì)等燃料的最佳選擇[1]。CFB燃燒技術(shù)始終向大型化、高參數(shù)方向發(fā)展。近年來，我國在超臨界循環(huán)流化床鍋爐發(fā)電技術(shù)方面取得了突破，引領(lǐng)世界在該領(lǐng)域的發(fā)展[2-5]。目前我國正在運行的超臨界CFB鍋爐已達(dá)40臺。同時，在超超臨界CFB鍋爐的研發(fā)、設(shè)計方面已開展大量工作[6-10]。2019年初，2項660 MW超超臨界CFB燃煤發(fā)電項目被正式列為國家電力示范項目，標(biāo)志著我國CFB發(fā)電開始邁向超超臨界階段。

在CFB鍋爐大型化發(fā)展歷程中，隨著鍋爐參數(shù)提高、蒸發(fā)量增大，鍋爐結(jié)構(gòu)也在不斷發(fā)展變化。首先，為了保證足夠的反應(yīng)容積和更多受熱面的布置，爐膛截面積和高度均相應(yīng)增加。600 MW超臨界CFB鍋爐爐膛高度已經(jīng)超過52 m，寬度超過25 m。清華大學(xué)[11-14]通過對超高爐膛縱向氣固流動特性的試驗和理論研究，實現(xiàn)了爐膛內(nèi)床壓的優(yōu)化，有效降低了水冷壁的磨損，提高了機組的經(jīng)濟性和可用性。其次，由于二次風(fēng)穿透能力限制，隨著爐膛截面積的增加，二次風(fēng)無法射入到爐膛中心區(qū)域，導(dǎo)致爐膛中心形成貧氧區(qū)，嚴(yán)重影響燃燒效率。為了解決該問題，可采用大寬深比的單布風(fēng)板爐膛，減少二次風(fēng)穿透所需要的射程，對于300 MW以上等級的CFB鍋爐，布風(fēng)板的寬深比通常在5以上；也可在爐膛底部從中間分開，形成雙布風(fēng)板的“褲衩腿”式爐膛結(jié)構(gòu)，二次風(fēng)從支腿內(nèi)側(cè)和外側(cè)同時注入爐膛[15]。對于更高參數(shù)的CFB鍋爐，也可考慮雙爐膛并聯(lián)和環(huán)形布風(fēng)的結(jié)構(gòu)形式[16-19]。第三，旋風(fēng)分離器的數(shù)量和布置方式不斷更新。為了保證分離器效率，分離器設(shè)計必須滿足特定準(zhǔn)則，不能簡單按比例放大，因此大型CFB鍋爐均采用多分離器并聯(lián)布置方案。對于300～350 MW的CFB鍋爐通常為3個分離器在一側(cè)并聯(lián)布置，白馬電廠600 MW超臨界CFB鍋爐為6個分離器；Foster Wheeler公司緊湊型布置的CFB鍋爐，采用較小直徑的分離器，分離器個數(shù)多于其他爐型；波蘭460 MW超臨界CFB鍋爐，在爐膛兩側(cè)共布置8個旋風(fēng)分離器。

由于鍋爐結(jié)構(gòu)上的變化，具有超大截面和多分離器并聯(lián)結(jié)構(gòu)的大型CFB鍋爐，在運行過程中發(fā)現(xiàn)的橫向參數(shù)非均勻性問題，逐漸成為關(guān)注的焦點[20-25]，主要包括熱力學(xué)橫向非均勻性、氣固組分的橫向非均勻性以及氣固流動的橫向非均勻性。熱力學(xué)橫向非均勻性主要是指溫度和傳熱系數(shù)在水平方向的偏差。張攀等[24]對國內(nèi)2臺300 MW CFB鍋爐進行了測試，在不同運行工況下爐膛內(nèi)均存在明顯的橫向溫度偏差和熱流密度偏差，同時也發(fā)現(xiàn)，在鍋爐熱負(fù)荷下降時，溫度和熱流密度橫向非均勻性更加嚴(yán)重。氣固組分橫向非均勻性主要是由于給煤進入爐膛密相區(qū)后，顆粒橫向擴散能力有限，使煤燃燒和放熱橫向分布不均。李金晶等[25]在300 MW CFB鍋爐上發(fā)現(xiàn)，沿橫向爐膛中間床溫明顯高于兩側(cè)，通過調(diào)節(jié)給煤點分配，可改變溫度分布。胡南等[26]通過總結(jié)已有橫向擴散系數(shù)研究成果，給出了大型CFB鍋爐密相區(qū)固體橫向擴散系數(shù)的取值范圍。

氣固流動不均勻性主要是由于爐膛結(jié)構(gòu)、氣固流動特性等原因?qū)е?，由于氣固流動對氣固組分分布和局部傳熱影響顯著，因此本文重點針對CFB鍋爐大型化過程中產(chǎn)生的新的氣固流動非均勻性問題，主要包括超大布風(fēng)截面的布風(fēng)非均勻性問題、多分離器并聯(lián)的氣固分配非均勻性問題以及特殊工況下產(chǎn)生的橫向波動和翻床問題。其中前二者在穩(wěn)定運行條件下，不隨時間變化，屬于靜態(tài)非均勻性問題；床壓波動和翻床問題屬于動態(tài)氣固流動橫向非均勻性問題。

1 靜態(tài)氣固流動非均勻性

1.1 多分離器并聯(lián)結(jié)構(gòu)

分離器是CFB系統(tǒng)重要的組成部分。分離器阻力特性是分離器設(shè)計的重要參數(shù)，對于多回路并聯(lián)的大型CFB鍋爐，分離器阻力特性決定各并聯(lián)回路的氣固流量分配。氣相流量一定的前提下，攜帶固相顆粒的質(zhì)量濃度將影響分離器壓降。早期試驗已發(fā)現(xiàn)，分離器壓降隨顆粒濃度增大呈先減后增的非單調(diào)特性[27-29]。當(dāng)顆粒濃度小于拐點處顆粒濃度時，固相濃度升高增加了顆粒與壁面撞擊的幾率，導(dǎo)致兩相流的動量損失增加，分離器內(nèi)湍流強度減小，導(dǎo)致分離器壓降減小。當(dāng)顆粒濃度進一步增加，大于拐點處顆粒濃度時，在分離器內(nèi)顆粒/顆粒團與氣體的相互作用增加，導(dǎo)致分離器壓降隨顆粒濃度的增加而增加[27]。

分離器的壓降ΔPc一般表述為

(1)

其中，fg、fp分別為氣相和顆粒相修正系數(shù)；Vi為分離器入口氣體速度，m/s；ρg為氣相密度，kg/m3。氣體修正系數(shù)fg主要取決于分離器入口和氣體出口的橫截面積之比。Baskakov等[27]提出了一個簡便的關(guān)聯(lián)式表征顆粒濃度對分離器壓降的影響，即

(2)

式中，Cp為分離器入口處顆粒濃度；a1、a2、a3為與分離器結(jié)構(gòu)有關(guān)的系數(shù)。

基于分離器壓降隨入口顆粒濃度增加先降后升的關(guān)系，在穩(wěn)定運行時，即使各分離器的壓降相等且氣體流量相同的情況下，入口顆粒濃度仍有可能不同，如圖1所示[30]。在理論研究方面，Grace等[31-32]提出并聯(lián)分離器系統(tǒng)內(nèi)氣固流動的多解方程，認(rèn)為系統(tǒng)穩(wěn)定運行的2個約束條件是各并聯(lián)分離器的壓降相等及氣體和固體的質(zhì)量守恒。Mo等[33]建立了并聯(lián)多分離器內(nèi)氣固流動的控制方程組。對于雙分離器并聯(lián)系統(tǒng)，氣固兩相流的控制方程組存在一個自由度，因而具有多值性，從理論上論證了分離器固體循環(huán)流率偏差的原因。

圖1 分離器壓降隨循環(huán)流率變化[30]Fig.1 Pressure drop of cyclone changing with solid circulating rate[30]

許多科研人員在多分離器的?；囼灧矫骈_展了相關(guān)工作，對于單側(cè)3分離器、4分離器，雙側(cè)布置的6分離器以及環(huán)形布置的6分離器結(jié)構(gòu)進行試驗研究，見表1。研究的核心問題是在不同風(fēng)速、不同床料量的運行條件下，各個回路的循環(huán)流率、分離器壓降以及提升管內(nèi)的壓力分布。中科院、浙江大學(xué)和重慶大學(xué)分別針對不同形式的多回路CFB鍋爐進行了數(shù)值模擬[6,39-40]。

表1 大型CFB鍋爐冷態(tài)?；囼?/p>

研究結(jié)果表明[35]，忽略并聯(lián)分離器及入口段間的微小加工差異，理論上對于純氣相的流動，各分離器的氣體流量相近，且分離器阻力越大，流量一致性越高，爐膛出口結(jié)構(gòu)偏差所引起的氣相流量偏差很小。在兩相流條件下，一側(cè)并聯(lián)3個分離器中，中間分離器循環(huán)流率和壓降均較小，且隨著風(fēng)速增加，循環(huán)流率偏差會進一步加大[6,35,38]。循環(huán)流率的偏差原因是內(nèi)墻邊壁處的顆粒濃度高，而兩側(cè)分離器附近的爐墻面積更大。王法軍等[37]對單側(cè)4分離器進行研究得到了相反的規(guī)律，即中間2個分離器的循環(huán)流率較高，循環(huán)流率偏差同樣隨風(fēng)速增加而增加，其主要原因是兩側(cè)分離器與中間2個分離器的入口角度不同。

爐內(nèi)屏式過熱器和尾部煙道對氣固流量分配的影響不可忽略。許霖杰等[7]、王超[39]等通過帶有懸吊屏的CFB鍋爐爐膛內(nèi)數(shù)值模擬發(fā)現(xiàn)，懸吊屏附近顆粒濃度較高，合理布置懸吊屏，可影響內(nèi)部流場，提高中間分離器的循環(huán)流率。郭強等[41]通過數(shù)值模擬研究了分離器出口煙道形式對分離器內(nèi)流量分配的反作用，研究發(fā)現(xiàn)，分離器出口煙道結(jié)構(gòu)對并聯(lián)三分離器內(nèi)的煙氣流動特性有較大影響。

1.2 超大截面的布風(fēng)均勻性

CFB鍋爐的布風(fēng)系統(tǒng)是CFB鍋爐的重要組成部分之一，包括風(fēng)室、布風(fēng)板以及風(fēng)帽。布風(fēng)系統(tǒng)的主要作用是使流化氣體沿截面均勻流出，同時防止顆粒落入風(fēng)室。CFB鍋爐風(fēng)室進風(fēng)方式主要包括兩側(cè)進風(fēng)和后墻進風(fēng)。相關(guān)測試和模擬研究均發(fā)現(xiàn)[42-47]，目前大型CFB鍋爐的布風(fēng)均勻性均有待提高，而進風(fēng)方式對于布風(fēng)均勻性的影響十分明顯。對于兩側(cè)進風(fēng)形式，風(fēng)量呈兩側(cè)風(fēng)帽流量高、中間風(fēng)帽流量略低的分布[44]；對于梯形風(fēng)室后進風(fēng)形式，表現(xiàn)為遠(yuǎn)離進風(fēng)口一側(cè)風(fēng)帽流量偏大[45]。布風(fēng)非均勻程度主要受布風(fēng)系統(tǒng)壓降與爐膛壓降之比(壓降比)影響。隨布風(fēng)面積的增加，布風(fēng)非均勻性進一步加劇，甚至出現(xiàn)局部床層壓力波動變大、流化死區(qū)、漏渣等現(xiàn)象，嚴(yán)重影響鍋爐的穩(wěn)定運行和燃燒效率。在布風(fēng)不均勻條件下，氣固流動橫向分布也存在偏差，同時反作用于供風(fēng)系統(tǒng)，床面氣固系統(tǒng)與布風(fēng)系統(tǒng)的耦合作用研究還較匱乏。

2 動態(tài)氣固流動非均勻性

2.1 鍋爐勻性動態(tài)非均勻現(xiàn)象

對于多分離器并聯(lián)的大型CFB鍋爐，分離器間氣固流動偏差是由于爐膛結(jié)構(gòu)、分離器特性產(chǎn)生，而布風(fēng)均勻性是由于風(fēng)室結(jié)構(gòu)和風(fēng)帽阻力等原因產(chǎn)生，在運行過程中偏差客觀存在并保持相對穩(wěn)定，屬于靜態(tài)橫向偏差。動態(tài)氣固流動的非均勻性，主要是指在不做干預(yù)的情況下，氣固流動的非均勻性不斷變化。對于褲衩腿型鍋爐，曾發(fā)生物料不斷在一側(cè)布風(fēng)板堆積的“翻床”現(xiàn)象，嚴(yán)重阻礙鍋爐的正常運行，隨著研究人員對該問題的深入認(rèn)識以及在運行過程中增加了優(yōu)化控制，“翻床”問題逐漸得到解決[25,48]，近期已鮮有相關(guān)研究報道。工程運行經(jīng)驗表明，大型CFB鍋爐在低負(fù)荷條件下會產(chǎn)生兩側(cè)壓力大幅度、長周期交替波動的動態(tài)橫向偏差現(xiàn)象(圖2)，并引起爐膛溫度偏差(圖3)，對鍋爐的安全運行產(chǎn)生更加嚴(yán)重的危害[49-50]。

圖2 CFB鍋爐床壓橫向波動[41]Fig.2 Lateral fluctuation of bed pressure in CFB boiler[41]

圖3 分離器出口煙氣溫度[41]Fig.3 Flue gas temperature of cyclone outlet[41]

李金晶等[21]認(rèn)為，形成床壓橫向波動的條件包括大尺度的爐膛截面、多條主循環(huán)回路并聯(lián)的結(jié)構(gòu)設(shè)計、母管制的流化風(fēng)供給系統(tǒng)和爐膛內(nèi)特定的流態(tài)。運行數(shù)據(jù)發(fā)現(xiàn)，只有當(dāng)一次風(fēng)量下降到一定值后，才會發(fā)生橫向波動現(xiàn)象。由于一次風(fēng)量的降低，布風(fēng)系統(tǒng)阻力下降，因此布風(fēng)穩(wěn)定性被破壞。另一方面，返料量等于分離器捕捉下來的相對大顆粒流量，但返料閥內(nèi)存有物料，因此返料閥對返料流量具有遲滯作用，而返料風(fēng)為母管制供風(fēng)系統(tǒng)，不同返料閥也會產(chǎn)生相互影響，這也是產(chǎn)生周期擾動的原因之一[20]。

2.2 床壓橫向波動理論

在理論方面，姜華偉等[49]建立了CFB鍋爐床壓橫向波動的非線性淺床長波駐波模型以及密相區(qū)氣固流態(tài)化數(shù)學(xué)模型，床面橢圓余弦波的波速C為

(3)

(4)

式中，Ug為流化風(fēng)速，m/s；Ug0為滿負(fù)荷運行條件下的密相床平均流化風(fēng)速，m/s；Kgs為氣固流動修正系數(shù)，無量綱；He為密相區(qū)膨脹床高，m；H為靜床高，m；X為波動幅度，m；k為橢圓積分或橢圓函數(shù)的模數(shù)K(k)和E(k)分別為第1類和第2類完全橢圓積分。

對于寬度為L的矩形截面流化床，床壓橫向波動周期T為

(5)

模型分析了靜床高對床壓橫向波動時爐內(nèi)床壓分布、床壓波動幅度和周期、活性區(qū)和非活性區(qū)床高、非活性區(qū)空隙率分布的影響。認(rèn)為隨著靜床高增加，床壓橫向波動幅度隨靜床高的增加呈正比增大趨勢，從而床面波速增加，床壓橫向波動往返周期隨之減小[49]。

文獻[50]建立了基于氣固系統(tǒng)固有頻率的壓力波動簡化模型，預(yù)測了平均床壓和波動周期的關(guān)系。但目前對于橫向波動僅停留在大型CFB鍋爐運行過程中實測數(shù)據(jù)和理論模型分析上，缺少系統(tǒng)的試驗數(shù)據(jù)支撐。文獻[45]中基于橫向波動固有頻率的假設(shè)給出了微分方程，在擾動條件下方程為

(6)

式中，β為阻尼系數(shù)，s-1；F0為橫向擾動強度，其物理意義為擾動使密相區(qū)所有顆粒產(chǎn)生的最大加速度，m/s2；ω0為氣固系統(tǒng)橫向波動的固有角頻率，s-1；Ω為擾動角頻率，s-1。

綜合考慮床尺寸、顆粒性質(zhì)、床壓以及流化風(fēng)速等因素，波動固有頻率的計算關(guān)聯(lián)式為

(7)

式中，ρs為顆粒密度，kg/m3；cs為顆粒體積濃度，無量綱；a為無量綱系數(shù)；L′為橫向波動發(fā)生方向幾何尺寸，m；Pb為平均床層壓降，Pa。

式(6)的近似解析式為

(8)

由式(8)可以看出，擾動F0越強，波動幅度越大；阻尼β越小、擾動頻率和固有頻率越接近，波動幅度越大。模型關(guān)鍵是選擇合理選擇阻尼系數(shù)β，理論上阻尼系數(shù)與流化風(fēng)速和顆粒性質(zhì)有關(guān)。由式(8)可以獲得波動幅度隨擾動頻率變化的趨勢，如圖4所示，擾動頻率與固有頻率相等時，波動幅度達(dá)到最大值。

圖4 波動幅度與擾動頻率的關(guān)系Fig.4 Relationship between pressure fluctuation and disturbance frequency

3 結(jié)語與展望

氣固流態(tài)化，從小截面的化工氣固反應(yīng)裝置，發(fā)展到具有超大截面的大型CFB鍋爐，眾多科研和工程技術(shù)人員通過理論、試驗和實爐測試開展了相關(guān)工作。目前中國在CFB燃燒大型化、高參數(shù)方面取得了世界領(lǐng)先。但CFB鍋爐大型化過程中，氣固流動的橫向非均勻性問題仍有部分工作需要進一步開展。

1)多分離器并聯(lián)的大型CFB鍋爐，氣固流量分配的均勻性主要取決于分離阻力特性、分離器布置、分離器入口煙道結(jié)構(gòu)、爐膛出口附近的懸吊屏布置以及分離器出口的煙道布置，在鍋爐設(shè)計過程中需綜合考慮，根據(jù)氣固流動橫向分布特性更合理安排受熱面的布置。

2)超大截面的流化布風(fēng)系統(tǒng)的進風(fēng)方式對于布風(fēng)均勻性具有顯著影響，大型CFB鍋爐由于布置和成本原因，側(cè)進風(fēng)和后墻進風(fēng)的方式均無法達(dá)到均勻性。非均勻布風(fēng)條件下氣固流動特性以及床層與布風(fēng)系統(tǒng)的耦合作用，還有待進一步研究。

3)大型CFB鍋爐低負(fù)荷條件下的橫向波動問題對鍋爐的安全運行威脅極大。大截面的布風(fēng)系統(tǒng)、一次風(fēng)降低導(dǎo)致的布風(fēng)失穩(wěn)以及多回路并聯(lián)的不均勻性是導(dǎo)致床壓橫向波動的主要原因，但目前在實驗室尺度還缺少系統(tǒng)研究。波動產(chǎn)生機理以及由此引申的橫向波動基礎(chǔ)理論仍需進一步深入分析。