濕法煙氣脫硫塔內(nèi)流場模擬及優(yōu)化

藺海艷，雷本喜

（江蘇科行環(huán)保科技有限公司，江蘇鹽城224051）

濕法煙氣脫硫塔內(nèi)流場模擬及優(yōu)化

藺海艷，雷本喜

（江蘇科行環(huán)保科技有限公司，江蘇鹽城224051）

以某熱電廠脫硫裝置為模擬對象，采用CFD流場模擬軟件對脫硫裝置進行幾何建模、網(wǎng)格劃分、數(shù)值模擬及計算等，分析脫硫裝置內(nèi)的速度分布、溫度分布、煙氣入射角及壓損等，并進行優(yōu)化，得到最佳結(jié)構(gòu)及內(nèi)部導流板布置形式。

脫硫；CFD；流場分布；托盤

引言

電廠燃煤鍋爐要達到ρ（SO2）≤35mg/Nm3的超低排放標準，石灰石-石膏濕法脫硫技術(shù)仍將是首選技術(shù)。該技術(shù)脫硫裝置內(nèi)煙氣流場既能影響脫硫效率，也會影響除霧效率。

本文運用CFD數(shù)值模擬技術(shù)對脫硫系統(tǒng)進出口煙道、脫硫塔內(nèi)部煙氣分布進行數(shù)字模擬分析，預測脫硫系統(tǒng)的煙氣分布，優(yōu)化設(shè)計導流板、聚氣環(huán)等裝置，從而使煙氣能夠充分均勻混合，保證流場的均勻性，并有效降低系統(tǒng)的阻力。同時對吸收塔內(nèi)部托盤、噴淋層、除霧器等進行近似模擬計算，探討吸收塔內(nèi)部流場分布，為脫硫系統(tǒng)的設(shè)計提供參考。

1　模型建立與計算

1.1幾何建模及網(wǎng)格劃分

按照某熱電廠脫硫系統(tǒng)裝置建立1∶1比例的三維幾何模型。模型建立之后，進行網(wǎng)格劃分將模型離散化，即通過有限的網(wǎng)格節(jié)點來描述實際的空間連續(xù)實體。考慮到模型復雜的內(nèi)部結(jié)構(gòu)，計算區(qū)域需采用結(jié)構(gòu)化和非結(jié)構(gòu)化的混合網(wǎng)格進行劃分，整個模型總網(wǎng)格數(shù)約460萬。

模型結(jié)構(gòu)復雜，做如下簡化：

1）以空氣代替煙氣，進入除霧器通道可把氣體視為不可壓縮氣體。

2）以水代替漿液。

3）忽略漿池部分，忽略煙氣流經(jīng)漿池液面時產(chǎn)生的擴容流動，將漿池液面作為固定、絕熱壁面來處理。

4）不考慮漿液池上方氣液兩相區(qū)域的化學反應(yīng)。

5）將托盤簡化為多孔介質(zhì)模型。

6）除霧器部分簡化為多孔介質(zhì)模型。

7）將各噴淋層中的主管和支管簡化為粗細均勻的圓柱管體。

8）忽略漿池部分，忽略煙氣流經(jīng)漿池液面時產(chǎn)生的擴容流動，將漿池液面作為固定、絕熱壁面來處理。

1.2數(shù)學模型及計算方法

1.2.1湍流模型

脫硫系統(tǒng)內(nèi)部的流動是一個復雜的湍流流動過程，考慮到數(shù)學模型的可靠性和工程應(yīng)用的可行性，本文選用標準k-?雙方程模型來模擬脫硫裝置內(nèi)部的煙氣流動情況。在直角坐標系下，基本控制微分方程如下：

以上控制方程求解時，采用了求解壓力耦合方程的全隱算法SIMPLE算法。

1.2.2多孔介質(zhì)模型

多孔介質(zhì)模型可以應(yīng)用于很多問題，如通過充滿的介質(zhì)、通過濾紙、穿孔圓盤、流量分配器以及管道堆的流動。對于本模型中托盤和除霧器的模擬，將由多孔介質(zhì)模型來實現(xiàn)。

當使用這一模型時，就定義了一個具有多孔介質(zhì)fluid單元區(qū)域，而且流動的壓力損失由多孔介質(zhì)的動量方程中所輸入的內(nèi)容來決定。事實上多孔介質(zhì)不過是在動量方程中具有了附加的動量損失而已。因此，多孔介質(zhì)模型具有以下限制：

1）流體通過介質(zhì)時不會加速，因為事實上出現(xiàn)的體積阻塞并沒有在模型中出現(xiàn)。這對于過渡流是有很大影響的，因為它意味著FLUENT不會正確地描述通過介質(zhì)的過渡時間。

2）多孔介質(zhì)對于湍流的影響只是近似的。在多孔介質(zhì)中，默認的情況下FLUENT會求解湍流量的標準守恒方程，因此介質(zhì)中的湍流被這樣處理：固體介質(zhì)對湍流的生成和耗散速度沒有影響。

多孔介質(zhì)的動量方程具有附加的動量源項。該源項由兩部分組成，一部分是粘性損失相（Darcy），另一部分是內(nèi)部損失項，如下：

其中：Si是i向（x、y或z動量源項），D和C是規(guī)定的矩陣。在多孔介質(zhì)單元中，動量損失對于壓力梯度有貢獻，壓降和流體速度（或速度平方）成比例。

對于簡單的均勻多孔介質(zhì)：

其中：α是滲透因子，C2是內(nèi)部阻力因子。

根據(jù)設(shè)計要求，本次模擬托盤阻力設(shè)定△P=500 Pa，除霧器阻力設(shè)定為△P=200 Pa，進而得到計算時內(nèi)部阻力系數(shù)（y方向即為煙氣流動方向）：

2　計算結(jié)果分析

2.1速度分布

對比脫硫裝置內(nèi)導流板前后煙氣速度分布，分析可知，通過在入口煙道拐角處增設(shè)導流板，可合理調(diào)節(jié)煙氣的分布，改善煙氣速度流場，減緩高速煙氣流對煙道的沖刷。

脫硫裝置內(nèi)煙氣速度分布，在脫硫塔內(nèi)部增設(shè)托盤之后，塔盤下方的煙氣在托盤壓力作用下自動調(diào)整，使得塔盤下煙氣速度分布均勻。

除霧器內(nèi)流速不均會造成局部堵塞或液滴逃逸，嚴重影響除霧器的效果，所以頂吸收塔內(nèi)除霧器上游的速度均勻性是影響脫硫性能的重要指標。頂層除霧器上游截面煙氣流速均勻，流速處于3.2～4 m/s之間，經(jīng)計算，該截面平均煙氣流速約3.44m/s，標準偏差為0.47，相對標準偏差為13.66%。一般認為速度相對標準偏差小于15%即滿足行業(yè)內(nèi)速度均勻指標。

2.2壓力分布

增加導流板后，煙氣速度分布更加均勻，高速流區(qū)得到緩解，整個系統(tǒng)阻力減小，經(jīng)計算，增設(shè)導流板后煙道系統(tǒng)阻力減小約100 Pa。

主要的壓力損失集中在彎道、托盤和除霧器層阻力。經(jīng)計算，100%負荷時，總壓降約1 045 Pa，其中托盤層壓降約占500 Pa，除霧器層的阻力約150 Pa，進口煙道阻力約100 Pa，出口煙道阻力180 Pa，脫硫塔本體約115 Pa，4層噴淋層約1 000 Pa，可算出脫硫塔煙道系統(tǒng)阻力為280 Pa，脫硫塔阻力為1 750 Pa，滿足該項目壓力損失要求。

2.3煙氣流動方向

在導流板的合理布置條件下，各段煙道內(nèi)的煙氣均能保證良好的方向均勻性，回流漩渦少，可有效避免煙氣沖刷和局部積灰現(xiàn)象。煙氣經(jīng)過托盤后，煙氣流向得到明顯矯正。經(jīng)計算，來流煙氣入射角與豎直方向夾角均小于10°，有利于煙氣與噴淋漿液形成對流，充分接觸混合，提高反應(yīng)效率。

3　結(jié)語

1）在脫硫系統(tǒng)煙道內(nèi)優(yōu)化設(shè)置導流板可優(yōu)化流場，減小系統(tǒng)阻力。

2）脫硫塔內(nèi)增設(shè)托盤，可使塔內(nèi)流場均勻，煙氣流向得到矯正，對提高脫硫效率和除霧效率有很大的幫助。

（編輯：賈娟）

The Flow Field Simulation and Optim ization of theW et Flue Gas Desulfurization

Lin Haiyan，Lei Benxi
（Jiangsu Cohen Environmental Protection of Science and Technology Co.，Ltd.，Yancheng Jiangsu 224051）

Based on a thermal power plant desulphurization device as simulation objects，the CFD flow field simulation softwarewas carried outon the desulfurization device geometrymodeling，meshing，numerical simulation and calculation，the velocity distribution in the analysis of desulfurization device，the temperature distribution and pressure loss of the incidence Angle and the flue gas，and optimize it，get the bestguide plate structure and the internal layout.

desulfurization；CFD；flow field distribution；tray

X701.3

A

2095-0748（2016）17-0026-02

10.16525/j.cnki.14-1362/n.2016.17.10

2016-07-29

藺海艷（1986—），女，山西呂梁人，研究生，助理工程師，研究方向：大氣污染物治理。