伍 豪,韋紅旗,周 帥,趙傳輝
(東南大學 能源與環(huán)境學院,南京 210096)
石灰石-石膏濕法脫硫工藝具有技術(shù)成熟、效率高、運行可靠等特點,廣泛應用于國內(nèi)燃煤電廠。在濕法脫硫工藝中,脫硫塔入口煙氣與脫硫漿液初期接觸區(qū)域長期處于干濕交界狀態(tài),脫硫塔和入口煙道的不合理設計易造成煙氣在脫硫塔入口產(chǎn)生回流及脫硫效率下降[1-2]。同時,煙氣的回流也會引起入口煙道內(nèi)壁及內(nèi)部桁架結(jié)垢,并使得系統(tǒng)阻力增大,長期結(jié)垢將嚴重削弱煙道的載重能力[3]。
經(jīng)濕法脫硫工藝處理后的煙氣溫度將降低至50 ℃左右,低溫煙氣會使煙囪出口出現(xiàn)白色煙羽,并伴有石膏雨,而且低溫濕煙氣會對尾部煙道、煙囪造成低溫腐蝕[4]。國內(nèi)外燃煤電廠普遍采用低低溫煙氣處理技術(shù)以避免該情況的出現(xiàn),而中間熱媒介煙氣換熱器(MGGH)是對傳統(tǒng)低低溫煙氣換熱器的改進。MGGH在運行過程中將空氣預熱器或靜電除塵器出口煙氣的熱量,以水為媒介傳遞給脫硫塔出口的低溫煙氣[5-6]。由于MGGH入口前存在變截面煙道與彎頭,易導致在入口處產(chǎn)生渦流,引起入口流場分布不均勻,換熱器進口流速不均勻會造成換熱效率的下降。BLECICH P[7]在對管翅式換熱器在空氣進口流量分配不均勻狀況下的研究中發(fā)現(xiàn),空氣側(cè)的流量分布不均勻?qū)е聯(lián)Q熱器的換熱效率下降,并引起換熱器的阻力增加。另外,在換熱器入口處產(chǎn)生的渦流容易造成管束振動,產(chǎn)生噪聲污染。因此,降低MGGH入口處流場的不均勻性,對解決換熱管束的振動及換熱效率下降等問題,保證MGGH安全、高效運行,并延長使用壽命具有十分重要的工程應用前景及意義。
實際工程中常采用加裝導流裝置的方式提高煙道內(nèi)流場的均勻性,以保證換熱器高效運行,防止換熱管束及煙道出現(xiàn)大幅度振動[8]。在轉(zhuǎn)彎煙道處增添導流板能明顯抑制流場的分離,有利于減少和阻止灰塵沉積,同時避免煙氣經(jīng)過轉(zhuǎn)彎煙道時,在離心力的作用下造成后續(xù)煙道內(nèi)流場的均勻性變差[9-12]。陽君等[13]在對列管式煙氣-煙氣換熱器(GGH)煙道系統(tǒng)的研究中發(fā)現(xiàn),同時采用導流板加布風板的組合可以抑制彎頭和變截面煙道處的流場分離,提升換熱器入口流場的均勻性。在對不同形態(tài)的導流板對流場影響的研究中發(fā)現(xiàn),帶弧度并在其后帶延長直板的導流板可以更好地引導煙氣流動,減少煙氣回流[14-16]。王為術(shù)等[17]運用數(shù)值模擬方法對選擇性催化還原(SCR)脫硝煙道系統(tǒng)內(nèi)的導流板進行優(yōu)化,提高了SCR脫硝煙道系統(tǒng)的流場均勻性,并使其適用于各負荷工況。在對構(gòu)建幾何模型進行數(shù)值模擬時,網(wǎng)格數(shù)目、質(zhì)量、計算邊界條件以及計算模型的選取都會影響數(shù)值模擬結(jié)果的準確性,因此需要將數(shù)值模擬結(jié)果與現(xiàn)場試驗數(shù)據(jù)進行對比。
筆者在對某燃煤電廠MGGH煙道系統(tǒng)的研究中,使用ANSYS FLUENT軟件以及多孔介質(zhì)模型對MGGH煙道系統(tǒng)進行數(shù)值模擬,采用加裝導流板的方式對MGGH煙道系統(tǒng)內(nèi)的流場進行優(yōu)化,并對優(yōu)化前后的模擬結(jié)果進行分析。
筆者研究的MGGH煙道系統(tǒng)整體模型按1∶1構(gòu)建,模型由原煙道、凈煙道、脫硫塔、MGGH四部分組成,構(gòu)建幾何模型時忽略其內(nèi)部煙道內(nèi)部支撐結(jié)構(gòu),同時忽略導流板的厚度。為提高數(shù)值模擬的效率,將MGGH煙道系統(tǒng)整體在脫硫塔的除霧器處分為兩個計算域(分別為原煙道計算域與凈煙道計算域),具體計算模型見圖1。由于MGGH內(nèi)部存在數(shù)目龐大的換熱管和翅片,以及復雜的支撐結(jié)構(gòu)和吹灰設備,若按原模型進行網(wǎng)格劃分,所需網(wǎng)格數(shù)目龐大,需要消耗大量的計算資源,因此筆者采用PATANKAR S V等[18]提出的多孔介質(zhì)模型對其進行簡化。
圖1 原煙道與凈煙道計算域
為確保數(shù)值模擬結(jié)果的準確性,原煙道與凈煙道計算域按結(jié)構(gòu)化網(wǎng)格劃分,并保證近壁面第一層網(wǎng)格的y+值在30以下,計算模型采用realizablek-ε模型,模型進出口邊界條件分別采用速度進口及壓力出口,并以現(xiàn)場試驗測得的MGGH煙氣冷卻器壓降(485 Pa)與煙氣再熱器壓降(406 Pa)為衡量標準,進行網(wǎng)格無關(guān)性驗證,模擬結(jié)果見圖2。由圖2可知:當原煙道、凈煙道計算域網(wǎng)格數(shù)分別大于251萬和237萬,MGGH煙氣冷卻器與煙氣再熱器的壓降趨于穩(wěn)定,而且壓降模擬值與現(xiàn)場試驗數(shù)據(jù)相符。因此,筆者選取的原煙道、凈煙道計算域網(wǎng)格數(shù)分別為287萬和299萬。
圖2 網(wǎng)格無關(guān)性分析
數(shù)值模擬過程中假設煙氣為不可壓縮牛頓流體,引風機出口處速度分布均勻,且煙氣在脫硫塔內(nèi)經(jīng)噴淋漿液沖刷后均勻地進入除霧器。MGGH采用多孔介質(zhì)模型簡化。多孔介質(zhì)模型是將多孔介質(zhì)區(qū)域內(nèi)的流體和固體簡化為同一個控制體,并用孔隙率來表示換熱管束的影響。在數(shù)值模擬的過程中,將固體對流體造成的影響當作是附加在流體上的動量源項增加至控制微分方程。在FLUENT軟件中,多孔介質(zhì)的動量方程中附加的動量源項由黏性損失項和內(nèi)部損失項組成。多孔介質(zhì)的壓降模擬公式[19]為:
(1)
式中:Si為i(x,y,z)向動量源項,Pa/m;vi為i(x,y,z)向速度,m/s;α為質(zhì)滲透面積,m2;μ為動力黏度,Pa·s;C2為壓力跳躍系數(shù),m-1;ρ為密度,kg/m3。
流場的均勻性在眾多工程研究中都是十分重要的衡量指標,為對各情況下流場的均勻狀況進行量化,通常采用現(xiàn)場試驗或數(shù)值模擬的方式獲得設置的各數(shù)據(jù)采集點的速度,然后將數(shù)據(jù)帶入速度偏差系數(shù)Cv[20]計算公式,Cv越小流場越均勻。
(2)
以速度偏差系數(shù)作為流場不均勻性評價的方法,需要進行采集點的設置、獲取采集點上的數(shù)據(jù)和處理數(shù)據(jù)等一系列操作,這將導致工作量的增大,為此筆者引入速度不均勻度[21]對流場均勻性進行評價。采用速度不均勻度可以直接在模擬的結(jié)果中對流場均勻性進行展示、分析,取消數(shù)據(jù)采集點的設置,為流場的優(yōu)化提供優(yōu)化依據(jù)。速度不均勻度越小,流場越均勻。
(3)
筆者將采用速度偏差系數(shù)和速度不均勻度兩種方式對流場的均勻性進行評價。
2.1.1 初始型原煙道分析
為便于對優(yōu)化前與優(yōu)化后的原煙道、凈煙道的流場模擬結(jié)果進行分析,定義在優(yōu)化前為初始型煙道,在初始型煙道基礎(chǔ)上進行優(yōu)化后為改進型煙道。
圖3是初始型原煙道的數(shù)值模擬結(jié)果,測得MGGH煙氣冷卻器壓降為485 Pa,與現(xiàn)場試驗數(shù)據(jù)相符。
由圖3可知:初始型原煙道存在3個流動紊亂區(qū)域,區(qū)域1、2有渦流產(chǎn)生。一方面,渦流的存在會造成鍋爐尾部煙道系統(tǒng)的阻力增大,風機運行能耗增加;另一方面,渦流位于距離引風機出口后較近的彎頭處會極大地影響引風機的安全運行。在區(qū)域3處煙氣多集中于一側(cè),致使脫硫塔入口處煙氣分布極不均勻,易造成在脫硫塔入口處出現(xiàn)漿液卷吸、結(jié)垢等現(xiàn)象。為使MGGH煙氣冷卻器及脫硫塔入口處于均勻的流動狀態(tài),需要對這些區(qū)域紊亂流場進行優(yōu)化。
2.1.2 原煙道改進型方案
針對在區(qū)域1、2處產(chǎn)生的渦流,以及在脫硫塔入口處煙氣嚴重偏向一側(cè)的情況,在引風機出口彎頭1、彎頭2、煙道匯合前后處、煙氣冷卻器出口后變截面傾斜煙道及脫硫塔入口前傾斜擴張煙道處增加導流板。安裝位置見圖4。
圖4 改進型原煙道導流板位置
2.2.1 初始型凈煙道分析
圖5為初始型凈煙道數(shù)值模擬結(jié)果,測得MGGH煙氣再熱器壓降為406 Pa,與現(xiàn)場試驗相符。
圖5 初始型凈煙道速度流場分布
由圖5可知:凈煙道內(nèi)主要的不穩(wěn)定流動分布在區(qū)域1、2、3。由于經(jīng)脫硫塔排出的煙氣聚集于出口水平煙道頂部,同時在煙氣再熱器進口前存在變截面擴張煙道,造成在變截面擴張煙道處煙氣充滿度不夠,區(qū)域1處流場紊亂,煙氣再熱器入口處的流速分布明顯不均。由于在出口水平煙道前存在2個連續(xù)的彎頭,造成在區(qū)域2、3處產(chǎn)生大渦流,并伴隨著流場嚴重偏心現(xiàn)象。因此,為使煙氣冷卻器入口獲得均勻的流動狀態(tài),消除凈煙道內(nèi)存在的大渦流,需要對這2個區(qū)域的不穩(wěn)定流場進行優(yōu)化。
2.2.2 凈煙道改進型方案
針對在煙氣再熱器入口流場速度不均勻,以及在區(qū)域2處產(chǎn)生的渦流情況,在脫硫塔出口水平煙道后變截面擴張煙道、出口彎頭1處增加導流板。安裝位置見圖6。
圖6 改進型凈煙道導流板位置
為更直觀清晰地展現(xiàn)初始型與改進型原煙道、凈煙道的流場,筆者選取了相關(guān)的截面進行流場分析;同時,對初始型與改進型的MGGH煙氣冷卻器、MGGH煙氣再熱器及脫硫塔入口的速度均勻度進行對比。
3.1.1 原煙道優(yōu)化后流場分布
圖7為改進型原煙道流場整體速度流線圖。
圖7 改進型原煙道流場分布對比
改進型原煙道在彎頭1、2處沒有渦流產(chǎn)生,雖仍存在少量紊亂的流線,但該處流場整體分布較為均勻;而且改進型煙氣冷卻器出口后的煙氣不再明顯地集中于一側(cè),脫硫塔入口處的流場得到極大的改善。
同時,為清楚地展現(xiàn)原煙道局部流場分布,選取Y=16.5 m處截面作為分析截面,具體截面的位置見圖8。
圖8 原煙道選取分析截面
圖9為原煙道選取分析截面速度流場分布圖。由圖9可知:相比于初始型原煙道的流場情況,改進型原煙道在MGGH煙氣冷卻器后的變截面煙道處的煙氣充滿度更高,且消除了在脫硫塔進口處煙氣流動嚴重偏心狀態(tài)。
圖9 原煙道Y=16.5 m截面速度流場分布
3.1.2 速度均勻度對比
圖10為MGGH煙氣冷卻器入口處速度不均勻度云圖。
圖10 煙氣冷卻器入口不均勻度云圖
由圖10可知:初始型原煙道在煙氣冷卻器入口處速度不均勻度大于±0.2的區(qū)域主要集中在入口左下角處。相較于初始型原煙道,改進型原煙道在煙氣冷卻器入口處速度不均勻度大于±0.2的區(qū)域減少,該處的流場狀況良好。MGGH煙氣冷卻器入口速度偏差系數(shù)由11.6%降低至9.3%,避免換熱管束振動,有利于MGGH煙氣冷卻器安全、高效運行。滿足設計要求(速度偏差系數(shù)小于10%)。
圖11為脫硫塔入口處速度不均勻度云圖。由圖11可知:初始型原煙道在脫硫塔入口處速度不均勻度大于±0.2的區(qū)域較大,且多集中于脫硫塔入口一側(cè);改進型原煙道在脫硫塔入口處的流場得到很大的改善,改進型原煙道在脫硫塔入口處速度不均勻度大于±0.2區(qū)域大幅度減少,優(yōu)化效果明顯。速度偏差系數(shù)由21.3%降低至6.3%,避免在入口處出現(xiàn)煙氣卷吸現(xiàn)象,保證了脫硫塔安全、高效運行,滿足設計要求。
圖11 脫硫塔入口不均勻度云圖
3.2.1 凈煙道優(yōu)化后流場分布
圖12為初始型和改進型凈煙道流場整體速度流線圖。
由圖12可知:改進型凈煙道在煙氣再熱器入口前變截面擴張煙道內(nèi)的煙氣充滿度更高,煙氣再熱器入口處的流場得到改善。在連續(xù)彎頭1處沒有大渦流產(chǎn)生,相比初始型凈煙道,在該處的流場整體分布更為均勻,增加導流板后優(yōu)化效果明顯。
圖12 凈煙道流場分布對比
3.2.2 速度均勻度對比
圖13為MGGH煙氣再熱器入口速度不均勻度云圖。
圖13 煙氣再熱器入口不均勻度云圖
由圖13可知:初始型凈煙道在煙氣再熱器入口處速度不均勻度大于±0.2的區(qū)域多;改進型凈煙道在該處的流場分布保持良好,改進型凈煙道速度均勻度大于±0.2的區(qū)域較少,煙氣再熱器入口速度分布更加均勻,速度偏差系數(shù)由27.7%降低至9.6%,滿足設計要求。
MGGH煙道系統(tǒng)內(nèi)原設計的流場狀況差,MGGH煙氣冷卻器、煙氣再熱器及脫硫塔入口處流場狀況無法滿足其工作要求。原煙道在煙氣冷卻器入口前的彎頭及傾斜擴張煙道,使得流場分布較不均勻,而且煙氣冷卻器出口后的變截面傾斜煙道使得煙氣聚集于一側(cè),造成脫硫塔入口速度分布偏心嚴重。凈煙道在煙氣再熱器入口前的變截面擴張煙道,造成其進口處流態(tài)嚴重偏心,而且在出口水平煙道前的2個連續(xù)彎頭處出現(xiàn)大渦流。故需要在原煙道及凈煙道內(nèi)加裝導流板以改善其內(nèi)部流動。
(1) 在不改變煙道結(jié)構(gòu)的條件下,通過增加導流板,可有效提升流場均勻性,并抑制渦流的產(chǎn)生、擴散。速度不均勻度能在模擬結(jié)果中清楚地展示流場的均勻狀況,為改善流場均勻性提供優(yōu)化方向。
(2) 改進型原煙道(加裝導流板后)有效提高了MGGH煙氣冷卻器入口的速度均勻性,入口大部分區(qū)域內(nèi)的速度不均勻度降低至±0.2以內(nèi),速度偏差系數(shù)降低至9.3%;同時,改善脫硫塔入口煙氣偏向一側(cè)的流場狀況,脫硫塔入口的速度偏差系數(shù)降低至6.7%。
(3) 凈煙道在變截面擴張煙道處流場偏心嚴重,且在出口水平煙道處渦流明顯,導流板的加入抑制了渦流的延伸,MGGH煙氣冷卻器入口大部分區(qū)域內(nèi)的速度不均勻度降低至±0.2以內(nèi),速度偏差系數(shù)降低至9.6%。