脫硫吸收塔深度除塵氣液兩相流數值模擬研究

蔣 楠，厲雄峰，葛春亮，劉文櫸，李晨朗，王先董

（1.浙江天地環保科技股份有限公司，杭州 311100；2.中國機械設備工程股份有限公司，北京 100055）

0 引言

近年來我國經濟社會快速發展，同時也產生了持續增長的能源消耗，引發了嚴重的大氣環境問題。顆粒物是我國城市大氣空氣中最主要的污染物之一，其中細顆粒物的污染尤為嚴重[1-4]。細顆粒物PM2.5 是指空氣動力學直徑小于2.5 μm的顆粒物，很容易懸浮在大氣中。由于其粒徑小，排放到大氣后隨氣流的傳輸距離更遠，在空氣中的停留時間也更長，所以其產生的污染影響范圍更廣、影響時間更長。

根據生成機理，大氣中顆粒物的來源主要可以分為一次顆粒物排放和二次顆粒物排放，其中，工業煙（粉）塵排放量為1 232.6 萬t，占全國排放總量的80.1%，工業是一次顆粒物最主要的排放源。而在工業行業中，電力、熱力生產和供應業是煙塵排放量的前3 位。因此，必須從源頭開始采取一系列手段，控制工業行業的顆粒物排放，以改善大氣環境質量。

石灰石-石膏濕法脫硫工藝是目前電力行業運用最為廣泛的一種脫硫技術[5]。研究表明[6-7]：對于粉塵而言，WFGD（濕法煙氣脫硫）系統的除塵效率能達到60%～80%，但該系統存在設備投資高、布置場地受限等不足[8]。在后超低排放時代，提高脫硫吸收塔協同除塵能力，挖掘吸收塔深度除塵潛力還需進一步研究[9-10]。脫硫吸收塔內部流場較為復雜，涉及到氣液固三相介質流動，為了解脫硫吸收塔對顆粒物的脫除機制，利用數值模擬技術進行分析是十分必要的。

1 模型建立

1.1 物模建立

噴淋脫硫吸收塔深度除塵過程包括氣液兩相耦合、氣液固三相流動以及液滴捕集細微顆粒物等復雜過程，本文基于Fluent 軟件先對脫硫塔內的氣液兩相流動進行數值模擬計算，分析塔內氣液兩相流動的分布情況，以便后續探究脫硫吸收塔協同除塵的捕集機制。

以某中試裝置中的噴淋脫硫塔為計算對象，開展氣液兩相流動模擬分析。中試脫硫塔的結構如圖1 所示。

圖1 脫硫塔結構示意

脫硫塔塔高29 m，除漿液池部分之外的塔徑為3 m，以漿液池底部平面的圓心作為坐標軸的原點；入口管道方向與水平面呈20°，入口中軸線的高度為6.46 m；均流板的設置高度為9.26 m，開孔率為0.40。塔內共設置5 個噴淋層，噴淋層截面布置如圖1 所示，8 個噴嘴均勻布置在截面的同一個圓周上，噴淋層編號為1—5 號，高度分別為11.26 m，13.26 m，15.26 m，17.26 m，19.26 m。

1.2 數模建立

1.2.1 基本假設

根據氣液兩相的流動特征和液相的體積分數，塔內氣液兩相模擬以氣相為連續相，脫硫液滴為離散相，即采用Eulerian-Lagrangian 多相流模型[11-14]。煙氣的流動視為定常不可壓縮粘性流動，以下為對氣相控制方程的描述。

質量守恒方程（連續性方程）：

動量守恒方程：

能量守恒方程：

式中：ρg為氣相密度；ug為氣相速度矢量；μg為煙氣的動力粘性系數；p 為氣相壓力；g 為重力加速度矢量；Fd為噴淋液滴對單位體積流體的作用力；e 為單位質量煙氣的內能；λ 為氣相導熱系數；Tg為煙氣的溫度；為煙氣與外界的能量交換。

噴淋塔內，漿液液滴主要受到重力、曳力、浮力、Basset 力、Saffman 力等作用力，相對于曳力和重力，其他力都很小，可以忽略[15]。液滴運動方程采用牛頓第二定律表示，氣流中單個液滴的平衡方程[16]為：

式中：md為漿液滴質量；ud為液滴速度矢量；FD為單個液滴受到的曳力；ρd為液滴密度；dd為液滴直徑；Red為液滴雷諾數；CD為曳力系數。

對式（4）進行時間積分就能得到離散液滴顆粒的速度，進而得到液滴顆粒的運動軌跡。

1.2.2 網格劃分

對吸收塔進行網格劃分，采用了非結構化網格和結構化網格相結合的方式。噴淋區由于考慮了噴淋管道的影響，結構比較復雜，故采用非結構化網格繪制，同時針對重要觀察區進行局部加密處理，其余部分均生成結構化網格。研究前期，針對該模型進行了網格無關性驗證，在本模型的網格數量下，網格的疏密程度對其模擬結果影響不大。網格劃分如圖2 所示。

圖2 網格劃分示意

1.2.3 邊界條件

表1 給出了典型工況下數值模擬計算的邊界條件和相關參數設置。

表1 氣液兩相流場邊界條件與參數設置

1.2.4 數值計算方法

本文的數值模擬計算基于Ansys Fluent 軟件進行求解。本文中，液滴的體積分數視為10%以下，可采用DPM（離散相模型）離散軌道模型在Lagrange 坐標系下計算液滴運動軌跡[17]。采用Realizable k-ε 模型在Euler 坐標系中進行計算，采用穩態計算，使用分離式求解器Pressure-Based Solver，差分方程的求解采用SIMPLEC 算法，動量、能量、湍動能等的離散方法均采用二階迎風格式[18-20]。

2 脫硫塔內氣液兩相流場分析

2.1 脫硫塔內液滴分布

圖3 給出塔內縱截面液滴體積分數分布，有助于更好地理解塔內液滴對顆粒物的捕集情況。

如圖3 所示，液滴分布在噴淋層及其下方的空間，入口附近以及出口區均不存在噴淋液滴。塔內共設置5 個噴淋層，可以看出液滴濃度的變化趨勢在高度上會呈現一定的周期性，該周期性會影響氣液兩相速度場、溫度場的分布。

圖3 液滴體積分數縱截面分布云圖

從圖4 可以看出每個分噴淋區內橫截面上液滴體積濃度分布會隨著高度不同而不斷變化。圖4（a）所示的平面緊挨噴嘴下方，可以看到8 個明顯的液滴密集圈，分別對應8 個噴嘴，平面中心較大一部分區域液滴體積分數較小或幾乎為0，說明在每個噴淋層所對應的分噴淋區內，該層的液滴基本被壁面或管道捕集。圖4（b）所示的平面，之前的濃度聚集點已擴張成較大的圓，而平面中點由于各噴嘴噴出的液滴跡線在這里有交叉，形成了局部液滴高濃度區。圖4（c）所示的平面，液滴的體積分數分布呈現錯落規則的花瓣狀，但是平面內的體積分數較小，最高只有0.23%。圖4（d）所示的平面，只有壁面附近有部分液滴。根據液滴的運動方向，該平面內液滴大部分都將要撞擊到壁面而被塔壁捕集，難以繼續對顆粒物數量的減少做出貢獻。

圖4 液滴體積分數橫截面分布云圖

在高度7～20 m 選取30 個平面，計算液滴的面平均體積分數，得到隨高度的變化曲線如圖5所示。

圖5 液滴橫截面平均體積分數隨高度的變化

液滴體積分數的變化曲線一共出現5 個峰，峰值所在的塔內高度與噴淋層的位置相對應。液滴體積分數高，才能有更多的液滴對顆粒物進行捕集，故脫硫塔除塵效率的提升需要更高的液滴體積分數。

2.2 氣液兩相速度場

煙氣從塔底部的入口進入噴淋脫硫塔后，氣相主流從下至上流動，流動過程中經過均流板。均流板可以讓煙氣分布趨于均勻，保證液滴與氣相很好地混合，有利于脫硫塔中脫硫漿液對SO2的吸收，同時也有利于液滴對顆粒物的捕集。

從圖6 可以看到，液滴在豎直方向上的運動距離較短，上層噴淋出的液滴幾乎無法運動到下層噴淋以下；所有的液滴跡線未到達煙氣入口處就已終止，液滴未能運動到塔底部的入口附近并繼續起到捕集顆粒或吸收其他污染物的作用。

圖6 液相速度分布

從圖7、圖8 來看，在每個分噴淋區的中間高度附近，近壁處區域的相對速度數值會高于其余區域的值，而兩相相對速度最大值出現在入口附近的對側壁面，此處對應的慣性捕集效率應該最高，有利于顆粒物的脫除。但是入口對側的壁面附近，液滴均集中在塔壁附近且體積分數很低，沒有足夠的液滴參與顆粒物捕集。

圖7 氣液兩相相對速度縱截面分布云圖

圖8 氣液兩相相對速度橫截面分布云圖

2.3 氣液兩相溫度場

煙氣以400 K 的初始溫度從脫硫塔入口進入，在塔內的溫度分布如圖9 所示。

圖9 氣相溫度縱截面分布云圖

可以看出，縱截面等溫線的分布呈現下凹形狀，尤其是在1—3 號噴淋層下方形成不同的分區，有明顯的分界主要是由于相應位置的液滴體積分數較高。其中，凹下去的區域為前文所述的噴淋層下方的液滴跡線交叉點，而分界線中位置較高的部分處于噴嘴所在的位置，也是液滴的聚集區，因此煙氣在這些區域與液滴發生熱交換失去較多熱量，使得煙氣溫度明顯下降。

氣液兩相溫度差在塔內的分布如圖10 和圖11 所示，塔內沒有液滴的區域兩相相對溫度設為0。

圖10 氣液兩相相對溫度縱截面分布云圖

圖11 氣液兩相相對溫度橫截面分布云圖

液滴在塔內的溫度變化小，塔內兩相溫度差的變化主要受到煙氣溫度變化的影響。由圖11 可以看到，在1 號噴淋層以下氣液兩相溫差明顯較高，主要因為該區域的煙氣溫度高，兩相溫差最高處位于入口對側壁面處，達到76.2 K。固體顆粒在有溫度梯度的流場中，由于其兩側的高溫區和低溫區給予的壓力不同，會使其向低溫區遷移，有利于液滴對顆粒的熱泳捕集。而隨著塔內高度增加，氣液兩相溫差不斷減小，捕集效率也不斷減弱。

3 結語

氣液兩相相對速度越大，則液滴對顆粒物的慣性捕集效率就越高，有利于微米顆粒物在塔內的捕集。入口高度對側壁面附近，雖然氣液相對速度較大，但是液滴數目較少，體積分數較低；而壁面附近有更多相對速度較大的區域，則液滴的慣性效率較高，這些區域中脫硫塔對微米級顆粒的脫除作用更強，顆粒數濃度的減小應該更加明顯。

氣液兩相溫差越大，則液滴對顆粒物的熱泳捕集效率就越高，有利于亞微米顆粒物在塔內的捕集。但是塔內溫差最大的區域，液滴體積分數很低，故溫差較大區域內的除塵效果還要將各影響因素結合起來進行分析。