洪文鵬

(東北電力大學(xué)能源與動力工程學(xué)院)

目前，大部分火力發(fā)電廠脫硫系統(tǒng)采用石灰石-石膏濕法脫硫工藝，取消了氣氣熱交換器(Gas Gas Heater，GGH)裝置，直接將凈煙氣從煙囪排出，形成濕煙囪排放方案。無GGH裝置的脫硫系統(tǒng)投產(chǎn)后，雖然有效避免了GGH的堵塞問題，但由于無煙氣再熱措施，排煙溫度較低，煙氣自煙囪口排出后不能有效地抬升、擴散到大氣中，煙氣中攜帶的粉塵和液滴聚集在煙囪附近，落到地面形成石膏雨或出現(xiàn)酸雨現(xiàn)象。近年來，大型火電機組參與電網(wǎng)調(diào)峰，而現(xiàn)有的固定結(jié)構(gòu)除霧器難以適應(yīng)負荷大幅波動帶來的煙氣流量的劇烈突變，容易導(dǎo)致除霧效率下降。因此，分析除霧器在不同煙氣流量下的性能，研究提高除霧效率的方法具有重要的理論與現(xiàn)實意義。

筆者以濕法煙氣脫硫系統(tǒng)中帶鉤波紋板除霧器為研究對象，采用數(shù)值模擬技術(shù)求解鉤片直段長度H、鉤片圓弧段半徑r及鉤片圓弧段轉(zhuǎn)折角度β等因素對除霧器性能的影響[1]。提出了一種鉤片偏轉(zhuǎn)角度可調(diào)的高效帶鉤波紋板除霧器。采用數(shù)學(xué)建模及數(shù)據(jù)擬合等方法，建立機組負荷及燃料成分等參數(shù)與除霧器鉤片偏轉(zhuǎn)角度的關(guān)系[2]，并在此基礎(chǔ)上提出了除霧器的優(yōu)化運行方法。

1 模擬對象與簡化假設(shè)

1.1模擬對象

濕法煙氣脫硫過程如圖1所示，噴淋塔上部設(shè)置噴淋層，噴淋后的石膏漿液進入下部的反應(yīng)池，吸收液由漿液循環(huán)泵進行循環(huán)，噴嘴將吸收液分散成1～3mm的液滴并沿塔向下噴灑，煙氣由塔的下部引入，以一定速度向上流動，以此經(jīng)過冷卻區(qū)和吸收區(qū)，氣液兩相逆流接觸傳熱傳質(zhì)，吸收后的尾氣經(jīng)由除霧器后從塔頂排出[3，4]。筆者以除霧區(qū)的帶鉤波紋板除霧器作為研究對象。

圖1 石灰石-石膏濕法噴淋脫硫過程示意圖

模擬帶鉤波紋板除霧器的結(jié)構(gòu)如圖2所示，除霧器葉片轉(zhuǎn)折角β=90°，鉤片偏轉(zhuǎn)角度α=45°，第一、二級除霧器板間距分別為35、25mm。參考工程應(yīng)用情況[1]，筆者對鉤片直段長度H、鉤片圓弧段半徑r和鉤片圓弧段轉(zhuǎn)折角度β對除霧器特性的影響進行數(shù)值模擬研究。

圖2 除霧器結(jié)構(gòu)與計算區(qū)域網(wǎng)格劃分

1.2模型假設(shè)與簡化

根據(jù)實際WFGD工程的運行環(huán)境和滿足工程要求的條件，對除霧器內(nèi)煙氣與液滴兩相流動狀況做假設(shè)和簡化[5，6]：除霧器簡化為二維流動；不考慮漿液滴的碰撞聚并與液滴的破碎；液滴視為均勻硬質(zhì)球形；不考慮二次夾帶作用的影響。

2 數(shù)學(xué)模型與計算方法

2.1計算區(qū)域與網(wǎng)格劃分

計算區(qū)域采用結(jié)構(gòu)化網(wǎng)格，由于鉤片附近流場變化比較劇烈，因而在鉤片所在的區(qū)域網(wǎng)格劃分應(yīng)細密一些。經(jīng)過網(wǎng)格精度無關(guān)化驗證后，劃分網(wǎng)格總數(shù)約67 000，除霧器的網(wǎng)格劃分如圖2所示。

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

2.2.1煙氣流場控制方程

采用歐拉法將煙氣處理為連續(xù)相并對其進行描述?；谏鲜黾僭O(shè)與簡化，不可壓縮煙氣流場的連續(xù)性控制方程可表示為：

div(V)=0

(1)

動量守恒控制方程為：

(2)

(3)

(4)

式中Fx、Fy、Fz——漿液滴對煙氣流場的反作用力，N；

p——壓力，Pa；

V——煙氣速度，m/s；

u、v、w——煙氣在x、y、z方向上的速度，m/s；

ρ——煙氣密度，kg/m3。

采用時均方法將式(2)～(4)處理為雷諾平均Navier-Stokes方程。根據(jù)除霧器通道內(nèi)煙氣流動湍流實際情況，采用低雷諾數(shù)k-ε模型來模擬除霧器內(nèi)的湍流流動。

2.2.2液滴顆?？刂品匠?/p>

由于石膏液滴的體積分數(shù)遠小于10%，因此，采用拉格朗日法來描述液滴的運動較精確[7]，即選用離散相顆粒軌道模型來跟蹤液滴運動。液滴受到諸如重力、熱泳力及Saffman力等的作用，筆者主要考慮重力和曳力對液滴的作用[8]。根據(jù)液滴在煙氣流場中的受力分析，液滴運動方程可表示為：

(5)

式中FD(u-up)——液滴的單位質(zhì)量曳力；

Fx——熱泳力、布朗力及Saffman升力等作用力；

gx(ρp-ρ)/ρp——顆粒在運動中受到的重力和浮力。

2.3模型的邊界條件

計算平臺采用Fluent 14.0軟件，用Gambit軟件生成網(wǎng)格。連續(xù)相方面，介質(zhì)為密度1.1 kg/m3、動力粘度19.548 6μPa·s的空氣，假設(shè)氣相在入口截面的速度均勻分布，出口表壓0Pa[9，10]；離散相方面，介質(zhì)為密度1 200kg/m3、動力粘度54.9mPa·s的脫硫循環(huán)漿液，給定液滴的初始速度與煙氣入口速度一致，液滴在入口截面上均勻分布[11]。

3 數(shù)值模擬結(jié)果與分析

3.1數(shù)值計算結(jié)果準(zhǔn)確性驗證

參考楊柳的實驗數(shù)據(jù)[12]，將模擬結(jié)果與實驗結(jié)果進行對比，以驗證文中數(shù)值模型的準(zhǔn)確性。以實驗中采用的板型尺寸、煙氣流速V及霧滴粒徑d等參數(shù)作為模擬參數(shù)，如圖3所示，實驗結(jié)果與模擬結(jié)果符合較好，表明文中的模擬計算可以用于除霧器性能的模擬研究。

圖3 計算結(jié)果與實驗數(shù)據(jù)的比較

3.2結(jié)構(gòu)參數(shù)對除霧器效率的影響

隨著鉤片直段長度H的增加，除霧器效率呈現(xiàn)升高趨勢(圖4a)。原因是鉤片直段長度的增加會使氣流在通道內(nèi)發(fā)生劇烈偏轉(zhuǎn)，致使液滴在慣性力的作用下脫離連續(xù)相的流線直接撞擊鉤片表面的幾率增加，導(dǎo)致除霧效率上升。一方面鉤片可以直接攔截一定量的液滴，另一方面可以使攜帶液滴的煙氣加速，從而增加離心慣性力。鉤片直段長度直接關(guān)系到鉤片間通道攔截氣流的程度，使得除霧器通道內(nèi)湍流程度加劇，流場內(nèi)擾動的增強會提高除霧效率。相同的，圓弧段轉(zhuǎn)折角度β和圓弧段半徑r的增大也會增加鉤片對葉片通道間的攔截面積，從而使除霧效率呈現(xiàn)上升趨勢(圖4b、c)。

3.3結(jié)構(gòu)參數(shù)對壓力損失的影響

煙氣在折板式除霧器中流動時，會在流道彎曲的區(qū)域形成回流區(qū)。鉤片直段長度H的增加，使流道變得更加狹窄，氣流在通流面積減小時速度增大，離心力隨之呈平方增大，以致產(chǎn)生更大的漩渦。每一個回流區(qū)都會阻塞通道，使得主流加速，引起靜壓減小。隨著流道變寬、流體的逐漸減速和回流區(qū)下游的再附著，壓力會逐漸恢復(fù)，但是由于能量的耗散，壓力無法完全恢復(fù)，造成壓力損失。鉤片直段長度H對除霧器內(nèi)擾動影響最為明顯，隨著H的增加，壓力損失迅速上升(圖5a)。同樣，圓弧段轉(zhuǎn)折角度β和圓弧段半徑r的增加，也會引起壓力損失的上升(圖5b、c)。三者的結(jié)構(gòu)參數(shù)關(guān)系到除霧器通道內(nèi)的局部壓力損失，對除霧器總壓降貢獻較大。

圖4 不同結(jié)構(gòu)參數(shù)對除霧效率的影響

圖5 不同結(jié)構(gòu)參數(shù)對壓力損失的影響

3.4高效帶鉤波紋板除霧器葉片

在兩級除霧系統(tǒng)組合中，總壓力損失不超過150Pa，為達到理想的除霧效率，選擇了兩種波紋板除霧器葉片的組合方式。第一級除霧器：H=22mm、β=35°、r=33mm、D=35mm。第二級除霧器：H=19mm、β=30°、r=34mm、D=25mm。

在額定煙氣流速下，液滴粒徑d不同時兩種除霧器的除霧效率η見表1，第一級除霧器的壓力損失為73.14Pa，第二級除霧器的壓力損失為68.96Pa。

表1 液滴粒徑不同時兩種除霧器的除霧效率

4 典型300MW機組除霧器鉤片偏轉(zhuǎn)角度的優(yōu)化

4.1300MW機組系統(tǒng)基本參數(shù)

筆者針對某300MW亞臨界機組配備的WFGD系統(tǒng)開展除霧器的優(yōu)化研究。該WFGD系統(tǒng)可在鍋爐負荷40%～100%之間安全連續(xù)運行。液滴的粒徑按照某電廠實測數(shù)據(jù)輸入，通過Ros-in-Rammler方法來模擬，液滴的最小粒徑為10μm，最大粒徑40μm，平均粒徑21μm，分布指數(shù)3.77。在保證鍋爐正常燃燒用氧的前提下，適當(dāng)減小風(fēng)量，用C語言自編程序求解機組不同負荷與煤種下的煙氣量，當(dāng)燃用設(shè)計煤種(收到基碳、氫、氧、氮、硫、灰分、水分和低位發(fā)熱量分別為56.9%、4.4%、9.1%、1.2%、0.6%、14.8%、13.0%和22 415kJ/kg)時，額定工況下WFGD入口煙氣量為310m3/s，實際煙氣量與機組負荷之間的關(guān)系如圖6所示。

圖6 機組燃用設(shè)計煤種時煙氣量隨負荷的變化

4.2鉤片可調(diào)與固定結(jié)構(gòu)除霧器性能的對比

鉤片角度確定的相關(guān)因素較多，為了獲得最佳方案，對不同結(jié)構(gòu)參數(shù)除霧器進行了大量數(shù)值實驗。鉤片可調(diào)的除霧器模擬結(jié)果如圖7所示，與鉤片固定結(jié)構(gòu)(鉤片偏轉(zhuǎn)角度α=45°)除霧器比較可見，流場得到了明顯的加速和擾動。

圖7 固定鉤片與可調(diào)鉤片除霧器的比較

鉤片與水平位置夾角(鉤片偏轉(zhuǎn)角度α)逐漸減小的過程中，鉤片呈現(xiàn)張開的趨勢，增加了攔截面積，從而使一部分液滴直接撞擊到鉤片上而被捕集。鉤片與葉片之間流道變得更加狹窄，使攜帶液滴的煙氣加速，提高了液滴捕集效率。與此同時，鉤片背弧處產(chǎn)生的回流區(qū)也在逐漸增大，引起較大的能量損失。并且從工程實施來看，鉤片角度調(diào)節(jié)在制造和安裝階段也相對困難。

4.3除霧器的優(yōu)化運行

除霧器優(yōu)化運行的思路是：保證煙氣壓力能克服除霧器阻力的前提下，每級除霧器壓力損失不大于75Pa，合理調(diào)整鉤片偏轉(zhuǎn)角度，從而提高除霧效率。表2為案例機組在部分機組負荷下第一級除霧器的最優(yōu)運行方案(燃用設(shè)計煤種時)。

表2 燃用設(shè)計煤種時部分負荷下第一級除霧器最優(yōu)運行方案

4.4運行優(yōu)化結(jié)果

綜合以上計算結(jié)果，提出除霧器優(yōu)化運行方案。圖8和表3分別給出了除霧器運行優(yōu)化方案和運行優(yōu)化效果，并與原方案進行對照?？梢钥闯觯涸?0%～100%負荷工況下，僅通過調(diào)節(jié)鉤片偏轉(zhuǎn)角度，除霧效率就可以平均提高4.10%，并且在負荷較低時，除霧效率提高明顯。

5 結(jié)論

5.1固定結(jié)構(gòu)的除霧器不具有負荷適應(yīng)性，隨著機組負荷的降低，除霧效率下降。改變鉤片直段長度、圓弧段半徑、圓弧段偏轉(zhuǎn)角度和鉤片偏折角度能夠?qū)ΤF器內(nèi)氣液兩相流場進行調(diào)節(jié)，使之達到理想狀態(tài)，提高除霧效率。

圖8 原方案與優(yōu)化方案的比較

機組負荷L%鉤片偏轉(zhuǎn)角度α/(°)效率提高值Δη/%第一級除霧器第二級除霧器第一級除霧器第二級除霧器9043441．470．328041422．161．397039401．311．116035391．341．815032372．961．374029345．433．66

5.2以300MW機組WFGD系統(tǒng)為例，當(dāng)燃用設(shè)計煤種時，在40%～100%額定工況下，采用筆者提出的優(yōu)化方法，除霧效率平均提高4.10%。但可調(diào)鉤片的制造和加裝在工程上相對困難。

[1] 姚杰，仲兆平，周山明.濕法煙氣脫硫帶鉤波紋板除霧器結(jié)構(gòu)優(yōu)化數(shù)值模擬[J].中國電機工程學(xué)報，2010，30(14)：61～67.

[2] 任浩仁，李蔚，盛德仁，等.火電機組變工況下運行指標(biāo)應(yīng)達值的分析[J].中國電機工程學(xué)報，1999，19(9)：3～8.

[3] 林永明.大型石灰石-石膏濕法噴淋脫硫技術(shù)研究及工程應(yīng)用[D].杭州：浙江大學(xué)，2006.

[4] 祝杰，葉世超，吳振元，等.石灰石-石膏脫硫噴淋塔吸收區(qū)高度計算模型[J].化工學(xué)報，2014，65(8)：2896～2900.

[5] 林永明，高翔，施平平，等.大型濕法煙氣脫硫噴淋塔內(nèi)阻力特性數(shù)值模擬[J].中國電機工程學(xué)報，2008，28(5)：29～30.

[6] 趙毅，華偉，王小明，等.濕式煙氣脫硫設(shè)備中折線型擋板除霧器壓降的數(shù)值模擬[J].中國電機工程學(xué)報，2004，24(8)：215～218.

[7] 孫志春，郭永紅，肖海平，等.鼓泡脫硫塔除霧器除霧特性數(shù)值研究及實驗驗證[J].中國電機工程學(xué)報，2010，30(8)：68～73.

[8] Wang Y，James P W.Assessment of an Eddy-interaction Model and Its Refinements Using Predictions of Droplet Deposition in a Wave-Plate Demister[J].Institution of Chemical Engineers，1999，77(8)：692～698.

[9] 魏星，李偉力，凡鳳仙，等.脫硫塔氣固兩相流場優(yōu)化的數(shù)值模擬研究[J].中國電機工程學(xué)報，2006，26(7)：15～17.

[10] 徐鋼，袁星，楊勇平，等.火電機組煙氣脫硫系統(tǒng)的節(jié)能優(yōu)化運行[J].中國電機工程學(xué)報，2012，32(32)：22～28.

[11] 趙毅，華偉，王亞君，等.濕式煙氣脫硫塔中折線型擋板除霧器分離效率的數(shù)值模擬[J].中國動力工程學(xué)報，2005，25(2)：294～295.

[12] 楊柳.燃煤煙氣脫硫脫氮一體化工藝及技術(shù)研究[D].南京：東南大學(xué)，2006.