何 銘，郭雪巖，楊 帆，呂宏紅

（上海理工大學(xué) 能源與動(dòng)力工程學(xué)院，上海 200093）

火電廠燃煤過程中產(chǎn)生的二氧化硫是導(dǎo)致酸雨污染的重要因素，因此國(guó)內(nèi)燃煤電廠普遍加裝脫硫裝置。其中石灰石?石膏濕法煙氣脫硫技術(shù)是目前應(yīng)用最為廣泛的一種煙氣脫硫工藝。大部分濕法脫硫工藝采用氧化噴槍裝置將氧化空氣鼓入脫硫吸收塔攪拌池內(nèi)，進(jìn)而實(shí)現(xiàn)將亞硫酸鈣氧化為硫酸鈣，并生成石膏晶體的過程。該過程中，攪拌槳承擔(dān)著分散氧化空氣，促進(jìn)池內(nèi)氧化反應(yīng)的作用，空氣在攪拌池內(nèi)的分散程度會(huì)影響系統(tǒng)的脫硫效率。

關(guān)于氣液攪拌池的研究多集中在進(jìn)氣方式優(yōu)化、氣泡尺寸及分布的測(cè)量等方面。Chaumat等利用雙光學(xué)探針對(duì)氣液兩相流場(chǎng)中氣泡尺寸及氣泡平均速度進(jìn)行了測(cè)量，為氣液攪拌池的研究提供了一定的參考；宋月蘭等對(duì)氣液兩相攪拌槽內(nèi)的流體運(yùn)動(dòng)進(jìn)行了實(shí)驗(yàn)及數(shù)值模擬研究，結(jié)果表明，雙歐拉方法能較好地模擬槽內(nèi)氣液兩相流的流動(dòng)狀況；王鐵峰等利用光纖探頭技術(shù)對(duì)三相循環(huán)流化床中的氣泡大小及其分布進(jìn)行了研究，發(fā)現(xiàn)漿液速度對(duì)氣泡平均直徑的影響較小；鐘允攀等通過改變攪拌槳轉(zhuǎn)速來研究攪拌射流對(duì)氧化空氣分布的影響，發(fā)現(xiàn)氧化空氣分布與攪拌射流的長(zhǎng)度及射流出口速度正相關(guān)；張朝平等研究了不同氧化噴槍結(jié)構(gòu)對(duì)氣液攪拌池氣含率分布的影響；周屈蘭等研究了空氣流量對(duì)脫硫塔漿液池顆粒分布的影響，發(fā)現(xiàn)空氣流量對(duì)漿液池中的顆粒分布影響較為顯著。

現(xiàn)有文獻(xiàn)中，有關(guān)氣液攪拌池內(nèi)氧化空氣進(jìn)口位置對(duì)池內(nèi)氣體分布影響的研究較為少見，因此本文對(duì)5種不同氧化噴槍布局中氣液攪拌池兩相流場(chǎng)進(jìn)行數(shù)值模擬，以研究氧化噴槍位置對(duì)氣液攪拌池空氣分布的影響，從而為脫硫吸收塔氣液攪拌池的設(shè)計(jì)提供一定的參考。

1 模型與網(wǎng)格無(wú)關(guān)性驗(yàn)證

1.1 物理模型及數(shù)值方法

攪拌池直徑D為11 m，液位高度H為7 m。4臺(tái)側(cè)入式攪拌槳的安裝高度h = 3 m，攪拌槳水平安裝角θ = 5°，豎直安裝角α = 10°，攪拌軸長(zhǎng)度L = 1 m。攪拌池結(jié)構(gòu)示意圖如圖1所示。攪拌槳采用直徑為1.2 m的KCX1200型三葉寬旋槳式攪拌槳。

圖1 攪拌池結(jié)構(gòu)示意圖Fig.1 Structure of the stirred tank

本文分別設(shè)置了前置型、后置型、下置型、左置型和右置型5種氧氣噴槍布局。氧化噴槍布局三維示意圖及二維俯視圖如圖2所示，其中：三維示意圖中虛線框標(biāo)記的面為氧化空氣出口位置；二維俯視圖中灰色圓為氧化空氣出口位置。對(duì)氧化噴槍作簡(jiǎn)化處理，將其設(shè)置為長(zhǎng)度為0.5 m的圓柱形。前置型氧化空氣出口面圓心距底面2.5 m，距壁面1.5 m，方向?yàn)榇怪毕蛳?；后置型氧化空氣出口面圓心距底面2.5 m，距壁面0.63 m，方向?yàn)檠財(cái)嚢栎S方向；下置型氧化空氣出口面圓心距底面2 m，距壁面1 m，方向?yàn)榇怪毕蛏稀Ｒ陨?種布局中噴槍出口面圓心均能映射在攪拌軸安裝點(diǎn)與攪拌池中心點(diǎn)的連線上。而左置型和右置型布局中噴槍的位置則分別為由下置型布局中噴槍以攪拌軸為旋轉(zhuǎn)軸順時(shí)針和逆時(shí)針旋轉(zhuǎn)90°而成。

圖2 氧化噴槍布局三維示意圖及二維俯視圖Fig.2 Three-dimensional diagram and two-dimensional top view of air-lance layout

參照陳佳等的研究，采用周期性邊界條件處理交界面的1/4計(jì)算域，并對(duì)氣液攪拌池流場(chǎng)進(jìn)行數(shù)值模擬。使用多重參考系（MRF）方法模擬攪拌槳的旋轉(zhuǎn)，湍流模型采用standard k?ε方法，多相流采用雙歐拉方法，氣液間曳力采用Tomiyama模型，噴槍進(jìn)口采用速度進(jìn)口邊界。頂部氣液界面采用degassing脫氣邊界。根據(jù)Barigou等的實(shí)驗(yàn)研究，氣液攪拌釜大部分區(qū)域內(nèi)氣泡的尺寸在3～5 mm之間，氣泡直徑設(shè)為4 mm。脫硫吸收塔攪拌池運(yùn)行溫度為45 ℃。工況參數(shù)如表1所示。

表1 工況參數(shù)
Tab.1 Working conditions

參數(shù)名稱 數(shù)值空氣進(jìn)口速度/（m·s?1） 10空氣密度/（kg·m?3） 1.225水密度/（kg·m?3） 988.04水動(dòng)力黏度/（Pa·s） 0.005 988

本文中未考慮氧化反應(yīng)過程對(duì)模型的影響，原因有以下兩點(diǎn)：①本文主要研究目標(biāo)為通過優(yōu)化氧化噴槍布局使攪拌池內(nèi)空氣分布更均勻，以提高攪拌池內(nèi)脫硫效率，因此主要研究攪拌池流場(chǎng)而非氧化反應(yīng)過程；②考慮化學(xué)反應(yīng)模型時(shí)需加入亞硫酸鈣及硫酸鈣固相顆粒，并將空氣分為氧化氣體和非氧化氣體，多相流模型內(nèi)相數(shù)將增加到5相，從而增加計(jì)算量，因此本文中計(jì)算氣液攪拌池流場(chǎng)時(shí)未考慮化學(xué)反應(yīng)模型。當(dāng)未考慮所加入的固相顆粒對(duì)攪拌池內(nèi)化學(xué)反應(yīng)的影響時(shí)：①將無(wú)法獲取脫硫塔攪拌池中的顆粒沉降情況；②氧氣與亞硫酸鈣發(fā)生化學(xué)反應(yīng)會(huì)使氣相體積分?jǐn)?shù)降低，且這個(gè)過程將無(wú)法被觀察到。因此，可在下一步研究中使用考慮了化學(xué)反應(yīng)的模型，以獲取更準(zhǔn)確的模擬結(jié)果。

1.2 網(wǎng)格無(wú)關(guān)性驗(yàn)證

采用四面體非結(jié)構(gòu)網(wǎng)格對(duì)計(jì)算區(qū)域進(jìn)行離散，對(duì)前置型布局流場(chǎng)進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證，分別計(jì)算網(wǎng)格數(shù)量為250萬(wàn)、360萬(wàn)、600萬(wàn)和860萬(wàn)時(shí)的氣相分布。圖3給出了網(wǎng)格數(shù)量對(duì)氣相分布的影響。從圖中可知，網(wǎng)格數(shù)量為600萬(wàn)時(shí)的氣相分布與860萬(wàn)時(shí)的較為相似。同時(shí)，本文統(tǒng)計(jì)了各模型中氣相體積分?jǐn)?shù)大于5%區(qū)域體積，網(wǎng)格無(wú)關(guān)性驗(yàn)證結(jié)果如表2所示。綜合對(duì)比后確定，合適的網(wǎng)格數(shù)量為600萬(wàn)。

表2 網(wǎng)格無(wú)關(guān)性驗(yàn)證
Tab.2 Grid independence verification

網(wǎng)格數(shù)量/萬(wàn) 氣相體積分?jǐn)?shù)大于5%區(qū)域體積/m3 250 13.87 變化率/%55.30 360 7.68 19.27 600 6.42 3.41 860 6.20 0（基準(zhǔn)值）

圖3 網(wǎng)格數(shù)量對(duì)氣相分布的影響Fig.3 Influence of grid number on the distribution of gas phase

2 結(jié)果分析

2.1 氧化噴槍布局對(duì)攪拌池流場(chǎng)的影響

5種布局中氧化噴槍均位于攪拌槳附近，在此區(qū)域內(nèi)氣液兩相間相互作用較強(qiáng)，局部流場(chǎng)會(huì)對(duì)氧化空氣在整個(gè)攪拌池中的分布產(chǎn)生較大影響。5種布局中攪拌槳附近流場(chǎng)如圖4所示，其中攪拌槳安裝高度為3 m。由速度分布可看出，在攪拌槳推進(jìn)流體運(yùn)動(dòng)的過程中形成了一股由攪拌槳向攪拌池中心運(yùn)動(dòng)的高速攪拌射流。該射流使氧化空氣隨漿液向攪拌池中心運(yùn)動(dòng)，對(duì)氧化空氣的擴(kuò)散有增益作用，因此，應(yīng)盡可能使氧化空氣被攪拌射流所影響。前置型和下置型布局中，氧化空氣均向攪拌池中心運(yùn)動(dòng)，說明這兩種布局均使氧化空氣能夠有效地被攪拌射流所推進(jìn)。后置型布局中，氣體從攪拌槳后方進(jìn)入攪拌區(qū)域。由于攪拌槳的高速旋轉(zhuǎn)對(duì)空氣射流存在阻擋作用，較少氣體能夠直接通過攪拌槳，因此大部分氣體從攪拌槳側(cè)方流出后向上運(yùn)動(dòng)。對(duì)比其他4種布局中的攪拌射流方向可以看出，后置型布局中匯集在此處的氧化空氣對(duì)攪拌射流產(chǎn)生了不利影響。左置型和右置型布局中，由于氧化空氣出口在攪拌槳側(cè)方，受攪拌射流影響較小，因此，氧化空氣射流方向?yàn)樨Q直向上，鮮有空氣隨攪拌射流向水平方向運(yùn)動(dòng)。綜合以上分析可知，合理的氧化噴槍布局應(yīng)滿足兩點(diǎn)要求：①氧化空氣出口位置應(yīng)位于攪拌射流路徑上，以使氧化空氣能夠有效地被攪拌射流所影響；②鼓出的氧化空氣不能匯集在攪拌槳附近，以免對(duì)攪拌射流產(chǎn)生不利影響。綜合比較后可知，前置型和下置型布局符合這兩點(diǎn)要求。

圖4 5種布局中攪拌槳安裝高度處水平面流場(chǎng)Fig.4 Horizontal flow field under five types of stirrer layout

前置型和下置型布局中攪拌池垂直方向速度分布如圖5所示。前置型布局中，氧化空氣進(jìn)入攪拌池后，受浮力的影響其向上運(yùn)動(dòng)的速度較高，因此在攪拌池中的停留時(shí)間較短，不利于其在攪拌池中的氧化反應(yīng)。而下置型布局中，由于氧化空氣射流經(jīng)過攪拌槳區(qū)域，攪拌槳的轉(zhuǎn)動(dòng)對(duì)空氣射流的運(yùn)動(dòng)存在一定的阻礙作用，使其垂直方向速度更低，延長(zhǎng)了其在攪拌池中的停留時(shí)間，使攪拌池內(nèi)氧化反應(yīng)更加充分。為了完整展示攪拌池中氣相的速度分布，本文提取氣相體積分?jǐn)?shù)為10%的等值面，并在等值面上顯示垂直方向速度分布，結(jié)果如圖6所示。由圖中可見，下置型布局中流場(chǎng)氣相垂直方向速度更低，氧化空氣的分散程度遠(yuǎn)優(yōu)于前置型布局。因此，降低氧化空氣在氣液攪拌池中垂直方向的速度可以延長(zhǎng)其在攪拌池中的停留時(shí)間，亦可增加氧化空氣在攪拌池中的分布。對(duì)比垂直方向速度分布可知，下置型布局優(yōu)于前置型布局。

圖5 前置型和下置型布局中攪拌池垂直方向速度分布Fig.5 Vertical velocity distribution of the prefixal layout and underneath-type layout

圖6 前置型和下置型布局中氣相體積分?jǐn)?shù)為10%的等值面垂直方向速度分布Fig.6 Vertical velocity distribution of the prefixal layout and underneath-type layout in the isosurfacewith a gas volume fraction of 10%

2.2 氧化噴槍布局對(duì)空氣分布的影響

2.1 節(jié)分析了5種布局中攪拌槳附近流場(chǎng)，并對(duì)比了前置型和下置型布局中氧化空氣速度分布，得出下置型布局優(yōu)于其他4種布局的結(jié)論。本節(jié)將對(duì)比各布局中氧化空氣在攪拌池內(nèi)的分布，從空氣分布的角度評(píng)價(jià)5種布局的優(yōu)劣。

5種布局中對(duì)應(yīng)的氧化空氣分布如圖7所示。將氧化噴槍置于攪拌槳前方的前置型布局中，氧化空氣分布受攪拌槳射流影響偏向攪拌池中心，在攪拌池中部及上部區(qū)域，沿弧形路徑運(yùn)動(dòng)至攪拌池氣液界面；后置型布局中，攪拌槳的轉(zhuǎn)動(dòng)對(duì)氣流的水平方向運(yùn)動(dòng)存在一定的阻礙作用，氣流通過攪拌區(qū)域后主要沿垂直方向運(yùn)動(dòng)，因此氧化空氣主要分布在攪拌槳附近及其上方區(qū)域；左置型和右置型布局中，氧化空氣射流進(jìn)入攪拌區(qū)域后，大部分氣體沿垂直方向運(yùn)動(dòng)，少量氣體沿葉梢路徑旋轉(zhuǎn)至槳葉頂部后繼而沿垂直方向運(yùn)動(dòng)，其余氣體隨攪拌射流運(yùn)動(dòng)，最終沿S型路線垂直方向上運(yùn)動(dòng)至攪拌池氣液邊界；而在下置型布局中，氣體與攪拌射流混合，被打散成多股氣流后隨攪拌射流運(yùn)動(dòng)，且因攪拌射流強(qiáng)烈擾動(dòng)，其水平方向速度比前置型布局中更高，空氣射流更偏向攪拌池中心，因此，氧化空氣分布范圍更廣，提高了脫硫吸收塔攪拌池的氧化反應(yīng)效率。

圖7 氧化噴槍位置對(duì)氧化空氣分布的影響Fig.7 Effect of air-lance position on the air distribution

氧化空氣分布主要集中在攪拌池中部及上部區(qū)域。以攪拌槳中心所在高度（0 m）為基準(zhǔn)高度，在氣液攪拌池0～5 m高度區(qū)間內(nèi)每隔0.5 m截取一水平截面，統(tǒng)計(jì)各平面內(nèi)氣相體積分?jǐn)?shù)大于10%區(qū)域（定義該區(qū)域?yàn)檠趸諝夥植紖^(qū)域）面積，結(jié)果如圖8（a）所示。由圖中可知，下置型布局中氧化空氣分布區(qū)域面積明顯大于其他布局，下置型布局中1～3 m高度處的中部區(qū)間的氧化空氣分布區(qū)域平均面積為1.5 m左右，下部區(qū)間及上部區(qū)間氧化空氣分布區(qū)域面積在2 m以上，其他4種布局中氧化空氣分布區(qū)域面積總體均在0.25～1 m。統(tǒng)計(jì)攪拌池中氣相體積分?jǐn)?shù)大于10%區(qū)域體積，結(jié)果如圖8（b）所示。下置型布局中氣相體積分?jǐn)?shù)大于10%區(qū)域體積約分別為前置型和后置型布局中的2倍，且遠(yuǎn)大于左置型和右置型布局中的。因此，下置型布局中空氣擴(kuò)散程度最佳。

圖8 5種布局中氧化空氣分布Fig.8 Oxidation air distribution of five types of air-lance layout

3 結(jié) 論

采用設(shè)置了周期性邊界條件的1/4計(jì)算域模型對(duì)氣液兩相攪拌池進(jìn)行數(shù)值模擬。根據(jù)模擬結(jié)果，從攪拌槳附近流場(chǎng)、氧化空氣垂直方向速度分布及攪拌池內(nèi)空氣分布等3個(gè)方面分析了5種噴槍布局對(duì)氧化空氣在氣液攪拌池內(nèi)分布的影響，得出以下結(jié)論：

（1）前置型和下置型布局使氧化空氣能夠有效地被攪拌射流所影響，增強(qiáng)了氧化空氣的水平方向運(yùn)動(dòng)能力，使其分布范圍更廣；而后置型布局對(duì)攪拌射流產(chǎn)生了阻礙作用，左置型和右置型布局中氧化空氣無(wú)法被攪拌射流所推進(jìn)，這3種布局中，氧化空氣難以運(yùn)動(dòng)至攪拌池中心，分散程度較差。

（2）前置型布局中氧化空氣向上運(yùn)動(dòng)的速度較大，在攪拌池中停留時(shí)間較短，而下置型布局中，攪拌槳對(duì)氧化空氣垂直方向運(yùn)動(dòng)的趨勢(shì)有阻礙作用，延長(zhǎng)了其在攪拌池中的停留時(shí)間，使池內(nèi)氧化反應(yīng)更加充分。對(duì)比氧化空氣垂直方向速度分布可知，下置型布局優(yōu)于前置型布局。

（3）下置型布局中，氧化空氣與攪拌射流混合后被打散成多股氣流，其分散程度更高，其氣相體積分?jǐn)?shù)大于 10% 區(qū)域體積約分別為前置型和后置型布局中的2倍，且遠(yuǎn)大于左置型和右置型布局中的。可見，下置型布局中空氣擴(kuò)散程度最佳。采用將氧化噴槍置于攪拌槳下方的下置型布局能有效擴(kuò)散氧化空氣，提升脫硫吸收塔攪拌池內(nèi)的氧化效率。