羅 閩，趙伶玲，李偲宇

(東南大學(xué) 能源熱轉(zhuǎn)換及其過程測控教育部重點(diǎn)實(shí)驗(yàn)室，南京 210096)

?

空氣預(yù)熱器硫酸氫銨積灰的數(shù)值研究

羅 閩，趙伶玲，李偲宇

(東南大學(xué) 能源熱轉(zhuǎn)換及其過程測控教育部重點(diǎn)實(shí)驗(yàn)室，南京 210096)

建立了適用于液態(tài)黏性NH4HSO4沉積的積灰模型，對某SCR脫硝改造后的300 MW機(jī)組回轉(zhuǎn)式空氣預(yù)熱器蓄熱板的積灰過程進(jìn)行了數(shù)值模擬，分析了煙灰顆粒法向碰撞速度和角度、積灰強(qiáng)度、積灰厚度及積灰概率等積灰特性的變化規(guī)律.結(jié)果表明：在NH4HSO4沉積區(qū)，上、下板平均積灰強(qiáng)度分別為普通積灰區(qū)的7.3倍和4.5倍，且不同顆粒粒徑的積灰概率保持在90%左右，而在普通積灰區(qū)均小于35%；上板平均積灰強(qiáng)度約為下板平均積灰強(qiáng)度的2倍；積灰厚度與積灰強(qiáng)度的分布規(guī)律一致.

回轉(zhuǎn)式空氣預(yù)熱器； SCR脫硝； 積灰； NH4HSO4

選擇性催化還原法(SCR)以其高效的特點(diǎn)成為當(dāng)前尾部煙氣脫硝技術(shù)的主流，并在我國得到了廣泛應(yīng)用[1].SCR系統(tǒng)中，還原劑(如氨水、尿素等)在催化劑作用下，選擇性地將煙氣中的NOx還原為N2，降低煙氣中NOx的含量[2].SCR系統(tǒng)通常布置在省煤器與空氣預(yù)熱器之間，從SCR系統(tǒng)中逃逸的氨氣會與SO3、H2O反應(yīng)生成NH4HSO4[3].NH4HSO4在420～493 K時呈黏性液態(tài)[4]，而回轉(zhuǎn)式空氣預(yù)熱器的中低溫段正好在此溫度區(qū)間運(yùn)行，液態(tài)黏性的NH4HSO4會迅速粘在傳熱元件表面，吸附大量飛灰，進(jìn)而造成空氣預(yù)熱器堵塞，一般的吹灰器無法有效清除此積灰[5].

與普通無NH4HSO4沉積的低溫積灰相比，NH4HSO4沉積區(qū)的積灰由于黏性NH4HSO4的存在，兩者在煙灰顆粒與蓄熱板發(fā)生碰撞時的能量、速度衰減等方面均明顯不同，后者的積灰特性與松散性積灰存在顯著差異，需對二者分開進(jìn)行討論.現(xiàn)有文獻(xiàn)中，為便于建模和計算，對鍋爐低溫受熱面上的積灰大多采用松散性積灰的處理方法[6]，這種處理方法對SCR脫硝改造后回轉(zhuǎn)式空氣預(yù)熱器的積灰問題顯然不再適用.

筆者對SCR脫硝改造后回轉(zhuǎn)式空氣預(yù)熱器高溫段(無黏性NH4HSO4沉積)、中低溫段(存在黏性NH4HSO4沉積)積灰進(jìn)行建模，并使用該模型對某SCR脫硝改造后300 MW機(jī)組回轉(zhuǎn)式空氣預(yù)熱器蓄熱板的積灰過程進(jìn)行了三維數(shù)值模擬，分析了NH4HSO4對回轉(zhuǎn)式空氣預(yù)熱器積灰特性的影響，對分析和解決SCR脫硝改造后空氣預(yù)熱器的吹灰問題具有重要的理論指導(dǎo)意義.

1 計算模型及研究對象

1.1 數(shù)值計算模型

采用標(biāo)準(zhǔn)k-ε雙方程[7]模擬蓄熱板內(nèi)部湍流流動，煙灰顆粒粒徑分布采用Rosin-Rammler分布[8]，平均粒徑為25 μm，其中最大顆粒粒徑為100 μm，最小顆粒粒徑為1 μm.顆粒在煙氣中運(yùn)動時受到曳力作用會顯著改變其運(yùn)動軌跡，顆粒所受曳力采用文獻(xiàn)[9]的計算方法.

顆粒碰撞對積灰的計算十分重要，當(dāng)顆粒與蓄熱板壁面發(fā)生碰撞時，其黏彈性作用采用Maxwell模型[10]進(jìn)行描述：

(1)

式中：ξ為應(yīng)力，Pa；下標(biāo)D、S分別表示黏滯體和彈簧；θ為應(yīng)變；ζ為材料黏性系數(shù)；E為楊氏模量，Pa；t為時間，s.

本文中普通積灰區(qū)只考慮其彈性作用，忽略黏滯體作用；而在NH4HSO4沉積區(qū)既考慮彈性作用，又考慮黏滯體作用.

碰撞發(fā)生時，顆粒主要受到壁面對其法向力和切向力的作用.法向彈性系數(shù)Ψn由Hertz接觸理論[11]確定，切向彈性系數(shù)Ψt根據(jù)Mindlin-Deresiewicz接觸理論[12]確定，整個碰撞過程中，質(zhì)量為m的振子的阻尼系數(shù)χ與恢復(fù)系數(shù)e間的關(guān)系采用文獻(xiàn)[13]的計算方法.此外，碰撞過程中，在NH4HSO4沉積處，黏滯力的存在會使煙氣能量產(chǎn)生較大的衰減.在普通積灰區(qū)，恢復(fù)系數(shù)為0.9；而在NH4HSO4沉積區(qū)，由于液態(tài)NH4HSO4的強(qiáng)黏性，恢復(fù)系數(shù)會發(fā)生變化，為便于計算，將其取為定值0.4.

采用基于速度的積灰準(zhǔn)則作為顆粒與蓄熱元件碰撞時顆粒是否產(chǎn)生沉積的判定依據(jù)，臨界積灰速度由以下公式[14]確定：

(2)

(3)

式中：vcr為臨界積灰速度，m/s；μ為泊松比；dp為顆粒直徑，m；下標(biāo)p和m分別表示顆粒和受熱面.

顆粒與蓄熱板碰撞時，當(dāng)顆粒與蓄熱板的法向碰撞速度大于其在該處的臨界積灰速度時，顆粒因碰撞而產(chǎn)生的彈性儲能足以克服其附著能，產(chǎn)生反彈；反之，顆粒沉積在蓄熱板上.NH4HSO4沉積區(qū)受熱面的楊氏模量為4×105Pa，普通積灰區(qū)受熱面的楊氏模量則為7.91×109Pa，對于同一顆粒，其在NH4HSO4沉積區(qū)的臨界積灰速度比其在普通積灰區(qū)高約一個數(shù)量級.

空氣預(yù)熱器入口煙溫為612.3 K，根據(jù)NH4HSO4的物性，當(dāng)回轉(zhuǎn)式空氣預(yù)熱器受熱面溫度在420～493 K時，作為NH4HSO4沉積區(qū)，此段積灰采用NH4HSO4積灰計算方法.顆粒碰撞和積灰的各模型均通過編寫UDF程序，采用商業(yè)軟件Fluent進(jìn)行計算.

1.2 研究對象及網(wǎng)格劃分

所研究的對象為某SCR脫硝改造后300 MW機(jī)組回轉(zhuǎn)式空氣預(yù)熱器，其蓄熱板為CU型，上板為波紋板蓄熱板，下板為波形定位板蓄熱板，如圖1所示.該空氣預(yù)熱器的轉(zhuǎn)速為2.32 r/min，對一個旋轉(zhuǎn)周期(即25.86 s)內(nèi)空氣預(yù)熱器蓄熱板的積灰過程進(jìn)行計算，流體流動方向?yàn)閦方向.取CU型蓄熱板下部定位板的一個周期(寬度為38 mm，長度為空氣預(yù)熱器轉(zhuǎn)子高度1.6 m)的蓄熱元件作為計算域，網(wǎng)格尺寸為1.2 mm×1.2 mm×1.2 mm，總網(wǎng)格數(shù)約為38萬.

圖1 蓄熱元件幾何模型Fig.1 Geometric model of the regenerative unit

2 模型驗(yàn)證

通過對模擬結(jié)果與已有文獻(xiàn)實(shí)驗(yàn)數(shù)據(jù)[15-16]的比較，對所建NH4HSO4沉積的積灰模型進(jìn)行驗(yàn)證.首先按照文獻(xiàn)[15]的實(shí)驗(yàn)條件建立模型，邊界條件也按其條件進(jìn)行選取，其上、下板均為光滑平板，入口煙溫為600 K，流速為8.5 m/s.計算所得煙溫分布見圖2.從圖2可以看出，本文計算煙溫值與文獻(xiàn)煙溫值的誤差很小，精度較高.實(shí)驗(yàn)結(jié)果表明，煙溫超過493 K時，NH4HSO4呈氣態(tài)，NH3轉(zhuǎn)換生成NH4HSO4的轉(zhuǎn)換率超過90%，不發(fā)生沉積；在420～493 K煙溫區(qū)間內(nèi)，NH4HSO4呈黏性液態(tài)，且其轉(zhuǎn)換率達(dá)到60%～80%，NH4HSO4發(fā)生沉積，故在此煙溫區(qū)間內(nèi)，按NH4HSO4沉積區(qū)的處理方法進(jìn)行積灰計算符合實(shí)際情況.

圖2 煙溫分布與NH4HSO4轉(zhuǎn)換率的關(guān)系Fig.2 Gas temperature distribution vs. conversion rate of NH4HSO4

此外，本文計算得到的臨界積灰速度與文獻(xiàn)[16]中實(shí)驗(yàn)結(jié)果的比較見表1和表2.從表1和表2可以看出，采用本文模型計算得出的臨界積灰速度與實(shí)驗(yàn)值的誤差約為6%～15%，在允許的誤差范圍之內(nèi).因此，采用本文模型計算得到的臨界積灰速度作為判定顆粒是否沉積的準(zhǔn)則可以準(zhǔn)確反映顆粒沉積的基本特征.

表1 粒徑為54 μm顆粒的臨界積灰速度Tab.1 Critical deposition velocity for particles of size 54 μm

表2 粒徑為42 μm顆粒的臨界積灰速度Tab.2 Critical deposition velocity for particles of size 42 μm

3 計算結(jié)果與討論

3.1 溫度場

一個旋轉(zhuǎn)周期后，煙溫沿z方向的分布如圖3所示.從圖3可以看出，從煙氣入口到出口，煙溫從612.3 K逐漸下降至約440 K；在z>1.25 m時，煙溫降至493 K以下，NH4HSO4開始產(chǎn)生冷凝，黏附在蓄熱板表面，為NH4HSO4沉積區(qū)；而0

圖3 沿流動方向煙溫的變化曲線Fig.3 Gas temperature profile along flow direction

3.2 法向碰撞速度和角度

沿z方向，上、下板出現(xiàn)顆粒碰撞時，顆粒(既包括反彈顆粒，又包括沉積的顆粒)的法向碰撞速度和角度的計算結(jié)果見圖4.從圖4可以看出，沿著z方向，上板的顆粒法向碰撞速度由煙氣入口處的0.68 m/s降至出口的0.04 m/s.法向碰撞速度降低說明顆粒的動能在衰減，這主要是由顆粒與蓄熱板之間的碰撞造成的.法向碰撞速度降得越快以及越低，說明顆粒與蓄熱板碰撞造成顆粒的能量衰減越大.在上板，顆粒法向碰撞角度由煙氣入口處的3.94°減至出口處的0.24°，法向碰撞角度減小則說明隨著流動的進(jìn)行和顆粒的動能衰減，顆粒逐漸適應(yīng)煙氣流動，與主流的跟隨性越來越好.

在下板，顆粒法向碰撞速度由煙氣入口處的0.3 m/s降至出口處的0.1 m/s，沿z方向也存在由于顆粒碰撞蓄熱板造成的動能衰減，但動能衰減幅度較上板?。欢ㄏ蚺鲎步嵌扔蔁煔馊肟谔幍?°減至出口處的1°，其減小幅度僅為上板的1/4，說明在下板附近顆粒在入口處隨煙氣運(yùn)動的跟隨性較好，這主要是由于上、下板波紋面的結(jié)構(gòu)不同造成的.

圖4 顆粒法向碰撞速度和角度分布Fig.4 Profiles of normal impact velocity and angle

3.3 積灰強(qiáng)度

采用單位面積上沉積的煙灰顆粒質(zhì)量來衡量蓄熱板表面的積灰嚴(yán)重程度，即

(4)

式中：md為蓄熱板某位置處的積灰質(zhì)量；A為相應(yīng)位置的面積；O為積灰強(qiáng)度，積灰強(qiáng)度的增長速率dO/dt也可看成是積灰速率.

進(jìn)一步分析數(shù)據(jù)，得到一個旋轉(zhuǎn)周期后上、下板的積灰強(qiáng)度(見圖5).從圖5(a)可以看出，上板積灰強(qiáng)度明顯大于下板積灰強(qiáng)度，上板的平均積灰強(qiáng)度為198 μg/mm2，下板的平均積灰強(qiáng)度為98 μg/mm2，上板平均積灰強(qiáng)度約為下板平均積灰強(qiáng)度的2倍.與0

圖5(b)給出了不同時刻上、下板積灰強(qiáng)度沿z方向的分布.從圖5(b)可以看出，積灰強(qiáng)度隨時間增加而增大，在0

圖5 積灰強(qiáng)度分布Fig.5 Distribution of deposition intensity

3.4 積灰厚度

當(dāng)運(yùn)行一個旋轉(zhuǎn)周期時，計算所得上、下板的積灰厚度分布見圖6.從圖6可以看出，上板平均積灰厚度為112 μm，下板平均積灰厚度為55 μm.在0

圖6 積灰厚度分布Fig.6 Distribution of deposition thickness

3.5 積灰概率

引入積灰概率來衡量顆粒與蓄熱板碰撞時，顆粒沉積于蓄熱板表面的可能性：

ηd=nd/nimp

(5)

式中：nd為蓄熱板某位置處沉積顆粒的數(shù)量；ηd為積灰概率；nimp為與對應(yīng)位置發(fā)生碰撞的顆粒的總數(shù)，nimp值越大，說明顆粒越易沉積，反之，說明顆粒越難沉積.

對不同粒徑范圍內(nèi)上、下板普通積灰區(qū)(0

此外，從圖7還可以看出，在普通積灰區(qū)，顆粒的積灰概率隨著顆粒粒徑的增大而減小，粒徑越大越不容易沉積，當(dāng)粒徑為100 μm時，積灰概率僅為2%，而在粒徑為25 μm時，積灰概率為35%.這是因?yàn)轭w粒臨界積灰速度隨著顆粒粒徑的增大而降低，相對而言，小顆粒與蓄熱板碰撞時，其法向碰撞速度更易滿足積灰條件而黏附在蓄熱板表面.但總體而言，在普通積灰區(qū)，各粒徑顆粒的臨界積灰速度均處于較低水平，不易沉積.而在NH4HSO4沉積區(qū)，粒徑為100 μm和25 μm時的顆粒積灰概率分別為92%和90%，積灰概率幾乎與顆粒粒徑無關(guān)，均保持在90%左右的較高水平.這說明在此區(qū)域無論是大顆粒還是小顆粒，碰撞時均極易產(chǎn)生灰分沉積.這是由于黏性NH4HSO4的存在會導(dǎo)致臨界積灰速度比普通積灰區(qū)高約一個數(shù)量級.因此在NH4HSO4沉積區(qū)，顆粒與蓄熱板碰撞時的法向碰撞速度低于其臨界積灰速度，使得產(chǎn)生沉積的可能性較高(即積灰概率較大)，較大的積灰概率也會導(dǎo)致此區(qū)域的積灰強(qiáng)度和積灰厚度明顯大于普通積灰區(qū).

圖7 積灰概率分布Fig.7 Distribution of deposition probability

4 結(jié) 論

(1)NH4HSO4沉積區(qū)的積灰強(qiáng)度、積灰厚度和積灰概率等均明顯大于普通積灰區(qū)，積灰更加嚴(yán)重.

(2)空氣預(yù)熱器上板的法向碰撞速度和角度均高于下板的法向碰撞速度和角度.

(3)液態(tài)NH4HSO4的黏性作用對空氣預(yù)熱器積灰特性的影響十分顯著，在進(jìn)行積灰計算時，須予以考慮.

[1] 董建勛, 李永華, 馮兆興, 等. 選擇性催化還原煙氣脫硝反應(yīng)器的變工況運(yùn)行分析[J]. 動力工程, 2008, 28(1): 142-146.

DONG Jianxun, LI Yonghua, FENG Zhaoxing,etal. Operation analysis of flue gas denitrification reactors with selective catalytic reduction under variable working conditions[J]. Journal of Power Engineering, 2008, 28(1): 142-146.

[2] 李德波, 廖永進(jìn), 徐齊勝, 等. 我國SCR脫硝催化劑服役過程中運(yùn)行規(guī)律的研究[J]. 動力工程學(xué)報, 2014, 34(10): 808-813.

LI Debo, LIAO Yongjin, XU Qisheng,etal. Experimental study on running condition of the catalyst in a domestic SCR denitrification system[J]. Journal of Chinese Society of Power Engineering, 2014, 34(10): 808-813.

[3] 周洪煜, 張振華, 張軍, 等. 超臨界鍋爐煙氣脫硝噴氨量混結(jié)構(gòu)-徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)最優(yōu)控制[J]. 中國電機(jī)工程學(xué)報, 2011, 31(5): 108-113.

ZHOU Hongyu, ZHANG Zhenhua, ZHANG Jun,etal. Mixed structure-radial basis function neural network optimal control on spraying ammonia flow for supercritical boiler flue gas denitrification[J]. Proceedings of the CSEE, 2011, 31(5): 108-113.

[4] 曾壁群, 姚友工, 楊博, 等. SCR煙氣脫硝過程中硫酸氫銨的形成原因及清洗方法[J]. 華電技術(shù), 2015, 37(1): 40-41, 44.

ZENG Biqun, YAO Yougong, YANG Bo,etal. Formation reasons of ammonium bisulfate in SCR flue gas deNOxprocess and cleaning method[J]. Huadian Technology, 2015, 37(1): 40-41, 44.

[5] 解濤, 金鵬. 300 MW機(jī)組空氣預(yù)熱器冷端腐蝕堵灰原因及對策[J]. 華北電力技術(shù), 2010，40(3): 47-50.

XIE Tao, JIN Peng. Cause of corrosion and blocking ash at cold junction of air preheater of 300 MW unit[J]. North China Electric Power, 2010, 40(3): 47-50.

[6] 陳衡, 王云剛, 趙欽新, 等. 燃煤鍋爐低溫受熱面積灰特性實(shí)驗(yàn)研究[J]. 中國電機(jī)工程學(xué)報, 2015, 35(增刊1): 118-124.

CHEN Heng, WANG Yungang, ZHAO Qinxin,etal. Experimental investigation of ash deposition on low temperature heating surfaces in a coal-fired boiler[J]. Proceedings of the CSEE, 2015, 35(S1): 118-124.

[7] 秦赟, 關(guān)欣, 頓喆, 等. 波紋板式空氣預(yù)熱器內(nèi)流動換熱過程的數(shù)值模擬[J]. 動力工程學(xué)報, 2015, 35(3): 213-218.

QIN Yun, GUAN Xin, DUN Zhe,etal. Numerical simulation on fluid flow and heat transfer in a corrugated plate air preheater[J]. Journal of Chinese Society of Power Engineering, 2015, 35(3): 213-218.

[8] ROSIN P, RAMMLER E. Die kornzusammensetzung des mahlgutes im lichte der wahrscheinlichkeitslehre[J]. Kolloid-Zeitschrift, 1934, 67(1): 16-26.

[9] ZHU C, LIANG S C, FAN L S. Particle wake effects on the drag force of an interactive particle[J]. International Journal of Multiphase Flow, 1994, 20(1): 117-129.

[10] MAXWELL J C. On the dynamical theory of gases[J]. Philosophical Transactions of the Royal Society of London, 1867, 157: 49-88.

[11] HERTZ H. Ueber die beruehrung elastischer koerper(on contact between elastic bodies), gesammelte werke (collected works)[M]. Leipzig, Germany: Johann Ambrosius Barth, 1895.

[12] MINDLIN R D, DERESIEWICZ H. Thickness-shear and flexural vibrations of a circular disk[J]. Journal of Applied Physics, 1954, 25(10): 1329-1332.

[13] WEBER M W. Simulation of cohesive particle flows in granular and gas-solid systems[D]. Colorado, USA: University of Colorado at Boulder, 2004.

[14] BRACH R M, DUNN P F. A mathematical model of the impact and adhesion of microsphers[J]. Aerosol Science and Technology, 1992, 16(1): 51-64.

[15] MENSHA J, DUNN-RANKIN D, MUZIO L,etal. Ammonium bisulfate formation temperature in a bench-scale single-channel air preheater[J]. Fuel, 2011, 90(7): 2445-2453.

[16] ABD-ELHADY M S, RINDT C C M, WIJERS J G,etal. Particulate fouling in waste incinerators as influenced by the critical sticking velocity and layer porosity[J]. Energy, 2005, 30(8): 1469-1479.

Numerical Simulation of Ash Deposition with Adhesion of NH4HSO4in an Air Preheater

LUOMin,ZHAOLingling,LISiyu

(Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education,Southeast University, Nanjing 210096, China)

A model was built for ash deposition with adhesion of viscous liquid NH4HSO4, based on which numerical simulations were conducted to study the ash deposition process on the regenerative unit of a rotary air preheater in a 300 MW set after retrofit of its SCR denitrification system, so as to analyze the variation law of the ash deposition characteristics, such as the normal impact velocity and angle of particles, and the deposition intensity, thickness and probability, etc. Results show that the average deposition intensity in the area with NH4HSO4adhesion is 7.3 and 4.5 times of that without NH4HSO4respectively for the upper and lower wall; the deposition probability in the area with adhesion of NH4HSO4is around 90%, whereas the value in areas without NH4HSO4is lower than 35%. The average deposition intensity on the upper wall is 2 times of that on the lower wall. The distribution of deposition thickness is in the same pattern as that of deposition intensity.

rotary air preheater; SCR denitrification; ash deposition; NH4HSO4

2015-12-23

2016-02-04

羅 閩(1991-)，男，湖南婁底人，碩士研究生，研究方向?yàn)榛痣姀S低溫受熱面?zhèn)鳠崤c積灰. 趙伶玲(通信作者)，女，教授，電話(Tel．)：025-83793453；E-mail：zhao_lingling@seu.edu.cn.

1674-7607(2016)11-0883-06

TK224.9

A 學(xué)科分類號：470.30