洪文鵬，裴彩鋒，劉廣林

(東北電力大學(xué)能源與動(dòng)力工程學(xué)院，吉林吉林132012)

近年來(lái)，我國(guó)火力發(fā)電廠建設(shè)得到迅猛發(fā)展，隨著我國(guó)煤炭等能源資源緊缺，人們節(jié)能降耗意識(shí)的日益增強(qiáng).國(guó)家政策對(duì)環(huán)境保護(hù)的要求越來(lái)越高。目前，我國(guó)發(fā)電能源構(gòu)成中還是以煤為主(占80%)，雖然今后將大力開發(fā)西部水電資源，但預(yù)測(cè)在本世紀(jì)30-40年內(nèi)以煤為主的結(jié)構(gòu)還不會(huì)改變。到2020年全國(guó)總裝機(jī)容量將達(dá)到9.5億千瓦左右，發(fā)電量將達(dá)到42000億千瓦時(shí)左右，其中火電裝機(jī)比重仍然占70%?；痣姍C(jī)組在為人們提供電和熱的同時(shí)，產(chǎn)生了大量的廢棄物，嚴(yán)重污染環(huán)境。氨法脫硫技術(shù)是以氨作為吸收劑脫除燃煤煙氣中的二氧化硫，具有脫硫成本隨煤的含硫量增加而下降、裝置阻力小、工藝節(jié)能優(yōu)化、脫硫裝置可靠、裝置配備設(shè)備少、既脫硫又脫硝等優(yōu)勢(shì)，在電廠脫硫行業(yè)中逐漸開始應(yīng)用［1-2］。

噴淋塔是氨法脫硫設(shè)備的主體，塔內(nèi)流場(chǎng)的優(yōu)劣直接影響脫硫效率與運(yùn)行成本［3］。目前，我國(guó)尚無(wú)噴淋塔設(shè)計(jì)與制造標(biāo)準(zhǔn)，國(guó)內(nèi)的研究主要是針對(duì)截面為圓形的噴淋塔，截面為矩形的噴淋塔少見研究與應(yīng)用。由于矩形截面的噴淋在外部結(jié)構(gòu)建造及內(nèi)部設(shè)備吊運(yùn)安裝上比較簡(jiǎn)便，矩形截面的脫硫塔在國(guó)外的大機(jī)組上已有應(yīng)用。本文利用流體計(jì)算軟件(CFD)Fluent對(duì)幾種長(zhǎng)寬比不同的矩形截面噴淋塔內(nèi)部流場(chǎng)進(jìn)行數(shù)值模擬，以其為矩形截面噴淋塔的優(yōu)化設(shè)計(jì)作出研究。

1 模擬對(duì)象及基本參數(shù)

1.1 模擬對(duì)象

由當(dāng)前研究成果得知，煙氣進(jìn)入噴淋塔內(nèi)主要受進(jìn)氣口對(duì)面及兩側(cè)的墻壁限制，噴淋塔截面的形狀對(duì)煙氣流場(chǎng)有直接影響［4］;因此，本文以某電廠圓柱形噴淋塔為參照對(duì)象，建立截面積、高度相同，截面長(zhǎng)寬比不同的矩形截面噴淋塔。噴淋塔采用氨水立式噴淋，煙氣入口垂直于塔體，無(wú)導(dǎo)流板等優(yōu)化措施，煙氣從下部進(jìn)氣口進(jìn)入，經(jīng)過(guò)噴淋后從上部流出的方式。針對(duì)相同負(fù)荷下、噴淋層雙開的方式，考查截面長(zhǎng)寬比不同的噴淋塔的煙氣流場(chǎng)。噴淋塔漿液池液面以上部分見圖1［5］。

圖1 三種噴淋塔模型

1.2 噴淋塔的基本參數(shù)

表1 噴淋塔的基本參數(shù)

2 模型簡(jiǎn)化與數(shù)值模擬方法

2.1 模型簡(jiǎn)化與假設(shè)

(1)將煙氣視為不可壓縮牛頓流體;

(2)假設(shè)漿液滴為球形，不考慮漿液滴的碰撞、破碎及聚并;

(3)不考慮煙氣與液滴之間的傳質(zhì)、傳熱和化學(xué)反應(yīng);

(4)不考慮煙氣中的硫化物與氨水液滴之間的化學(xué)反應(yīng)。

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

3.1 基本控制方程

采用歐拉法將煙氣處理為連續(xù)相對(duì)其進(jìn)行描述?；谏鲜黾僭O(shè)與簡(jiǎn)化，不可壓煙氣流場(chǎng)的連續(xù)控制方程可表示為［6］:

3.2 動(dòng)量守恒方程

式(2)-式(4)中p為靜壓，ρ為煙氣密度，V為煙氣的速度，F(xiàn)x、Fy、Fz為漿液滴對(duì)煙氣流場(chǎng)的反作用力。采用時(shí)均方法可將式(2)-式(4)處理為雷諾平均的Navier-Stokes方程。

3.3 液滴顆粒運(yùn)動(dòng)方程

采用拉格朗日法來(lái)描述液滴的運(yùn)動(dòng)，即選用Discrete Phase Models(DPM)模型來(lái)跟蹤液滴運(yùn)動(dòng)。吸收塔內(nèi)的漿液液滴受到諸如重力、曳力、浮力等力的作用，本文主要考慮重力和曳力對(duì)液滴的作用，液滴顆粒的運(yùn)動(dòng)方程可表示為:

式中:FD(ug-up)為顆粒受煙氣的單位質(zhì)量曳力;ρp為液滴密度;g為重力加速度;ug、up分別為煙氣和液滴速度;FD為曳力系數(shù)。

3.4 液滴粒徑分布

采用Rosin-Rammler模型來(lái)描述液滴的粒徑分布:

3.5 煙氣與液滴間的作用

在DPM模型中煙氣(連續(xù)相)與液滴(離散相)之間的作用可表示為:

式中:Fo為其它作用力;mp為顆粒質(zhì)量。

4 數(shù)值模擬方法

計(jì)算體為漿液面以上部分，運(yùn)用Gambit建立三維模型，采用分塊劃分生成結(jié)構(gòu)化網(wǎng)格;運(yùn)用fluent6.3軟件，湍流模型采用標(biāo)準(zhǔn)k-e方程，紊流邊界條件選用水力直徑與紊流強(qiáng)度，壓力采用SIMPLE算法，近壁面處采用壁面函數(shù)修正。

5 模擬結(jié)果與討論

為比較幾種模型的流場(chǎng)，主要考察有無(wú)噴淋時(shí)三維模型內(nèi)煙氣跡線、煙氣速度、塔內(nèi)壓強(qiáng)的指標(biāo)［7］。

圖2 無(wú)噴淋時(shí)煙氣跡線

圖3 雙層噴淋時(shí)煙氣跡線

由圖2中可以看出，在模型Ⅱ中旋流最大，模型Ⅲ中旋流最小，這是由于模型Ⅱ中截面邊長(zhǎng)長(zhǎng)(Y)大于邊寬(X)，截面寬度與進(jìn)氣口寬度相近，煙氣進(jìn)入塔內(nèi)無(wú)法向兩側(cè)擴(kuò)散，直接沖擊入口對(duì)面墻壁，煙氣產(chǎn)生向上及向下的速度，因此在上、下方向產(chǎn)生的旋流較大，在模型Ⅲ中，截面邊長(zhǎng)(Y)小于邊寬(X)，煙氣可向兩側(cè)擴(kuò)散，沖壁現(xiàn)象較弱，由此產(chǎn)生的旋流較小。

雙層噴淋時(shí)，由于液滴的重力及粘性作用，煙氣流場(chǎng)受到壓制，液滴的噴淋對(duì)煙氣的分布起到重新分配的作用，由圖3可以看出三種模型中旋流區(qū)域均變小。為考察塔內(nèi)煙氣流速，在兩層噴淋層上、下及Z軸17.5米及20.5米處建立截面觀察煙氣速度分布。

圖4 雙層噴淋時(shí)煙氣速度分布

由圖4可以看出，由于液滴的噴淋作用，煙氣的動(dòng)能得到消耗，速度隨著高度的增加而減少，在速度分布方面，模型Ⅰ、Ⅲ中速度梯度較少，速度分布較均勻，延長(zhǎng)了煙氣與液滴接觸時(shí)間，提高了反應(yīng)效率;但模型Ⅲ中煙氣入口側(cè)墻壁附近煙氣速度較低，壓強(qiáng)較大。由圖5中可以看出，模型Ⅲ中壓強(qiáng)較高，阻力較大，因此會(huì)增大循環(huán)水泵與增壓風(fēng)機(jī)電耗，增加運(yùn)行成本，模型Ⅰ、Ⅱ中壓強(qiáng)相差較小。

圖5 雙層噴淋時(shí)塔內(nèi)壓強(qiáng)分布

6 結(jié) 論

(1)在同等負(fù)荷下，截面積相同的三種模型中，截面邊長(zhǎng)小于邊寬(Ⅲ)或邊長(zhǎng)等于邊寬(Ⅰ)的模型在流場(chǎng)的均勻性及塔內(nèi)煙氣速度的分布方面要優(yōu)于截面長(zhǎng)大于寬(Ⅱ)的模型;

(2)模型Ⅲ的流場(chǎng)較模型Ⅰ均勻，但進(jìn)氣口墻壁側(cè)速度較低，壓強(qiáng)較高，致使阻力增大，增加了運(yùn)行電耗，但具體的經(jīng)濟(jì)性，還需考慮考慮建材造價(jià)、運(yùn)行成本等指標(biāo);

(3)對(duì)比分析模擬結(jié)果，截面邊長(zhǎng)小于邊寬或邊長(zhǎng)等于邊寬的模型經(jīng)濟(jì)性較好，模擬結(jié)果可作為矩形截面噴淋塔的工程設(shè)計(jì)及運(yùn)行優(yōu)化的依據(jù)。

［1］鄧新云，肖懷秋，禹練英.我國(guó)燃煤煙氣脫硫技術(shù)研究進(jìn)展［J］.廣州化工，2008，36(1):24-27.

［2］Michalski J A.Aerodynamic Characteristics of flue gas Desulfurization Spray Tower Spoly Dispersity Consideration［J］.Ind.Eng.Chem.Res，2000，39(9):3314-3324.

［3］郭偉，李彩亭，曾光明，等.傘罩型濕式脫硫除塵塔入口結(jié)構(gòu)優(yōu)化模擬［J］.環(huán)境工程學(xué)報(bào)，2009，3(6):1060-1064.

［4］汪洋.濕法脫硫噴淋塔數(shù)值模擬［D］.北京:華北電力大學(xué)，2006.

［5］賈勇，鐘秦，杜冬冬，等.氨法煙氣脫硫模型的研究［J］.動(dòng)力工程，2009，29(10):960-965.

［6］林永明，高翔，施平平.大型濕法煙氣脫硫噴淋塔內(nèi)阻力特性數(shù)值模擬［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2008，28(5):28-33.

［7］肖國(guó)俊.濕法脫硫噴淋塔煙氣入口角度優(yōu)化數(shù)值模擬［J］.電站輔機(jī)，2007(3):18-21.