流體模擬試驗技術(shù)在鍋爐SCR反應(yīng)器設(shè)備優(yōu)化中的應(yīng)用

張峰

（華電新疆發(fā)電有限公司昌吉熱電廠，新疆 昌吉 831100）

1 SCR反應(yīng)器內(nèi)流體模擬系統(tǒng)

昌吉熱電廠1、2號機組原脫硝裝置采用選擇性催化還原法(SCR)脫硝技術(shù)，單臺機組脫硝催化劑分為“2+1”層設(shè)計，初裝兩層，脫硝裝置入口NOx濃度按450mg/Nm3(標(biāo)干，6%O2)，出口NOx濃度不高于100mg/m3（標(biāo)態(tài)，干基，6%O2）設(shè)計。CFD模型的建立主要是從省煤器下部到空預(yù)器入口法蘭，如下圖1所示為初步設(shè)計階段的SCR系統(tǒng)形狀示意圖。如圖2所示氨氣及空氣混合其他主要通過噴氨格柵噴入SCR脫硝系統(tǒng)煙道并與其內(nèi)部的煙氣混合，最后進入SCR反應(yīng)器內(nèi)，在催化劑的作用下NH3與NOx發(fā)生化學(xué)反應(yīng)生產(chǎn)水和氮氣，煙氣脫硝后經(jīng)SCR出口煙道排出，然后進入空氣預(yù)熱器等下游各級輔機設(shè)備。

圖1 CFD模型圖

通過利用加裝導(dǎo)流板的形式可以將煙氣進入催化劑層前盡可能分布均勻，為確保氣流的均勻性需要在每一個煙氣流向改變處及變徑段設(shè)置相應(yīng)的導(dǎo)流板。如圖1所示，導(dǎo)流板統(tǒng)一安裝在調(diào)整管道的頂部及反應(yīng)器入口位置。另外為了確保脫硝的效率，通常需要在噴氨格柵前加裝導(dǎo)流板來使氨和煙氣能夠充分混合，如圖2，同時導(dǎo)流板的安置還要充分考慮到系統(tǒng)阻力的影響情況。

圖2 噴氨格柵、上升煙道頂部導(dǎo)流板

2 工況模擬分析

（1）數(shù)學(xué)模型。在數(shù)值模擬實驗中，科學(xué)合理地選取數(shù)學(xué)模型對后期的計算結(jié)果至關(guān)重要。本次模擬對象主要針對流場的改進設(shè)計，所以在確保高精度的條件下，為了減少計算量，主要應(yīng)用氣相湍流模型。依據(jù)SCR裝置內(nèi)煙氣流動時湍流的情況，主要采取標(biāo)準(zhǔn)k-ε湍流模型來模擬系統(tǒng)內(nèi)煙氣的湍流運動。在正交直角坐標(biāo)系下，k-ε湍流模型基本控制方程可以表示：

其中，φ為因變量，Γφ為因變量φ的擴散系數(shù)，Sφ為因變量守恒方程中所對應(yīng)的源項。

（2）多孔介質(zhì)模型。對于SCR反應(yīng)器內(nèi)的催化劑層壓降，通過將催化劑層看作多孔介質(zhì)進行模擬。其壓降損失模擬公式如下：

式中：Si——i方向上動量源項，Pa/m；μ——流動動力粘度，Pa·s；α——介質(zhì)滲透性；vi——i向速度分量，m/s；ρ——密度，kg/m3；C2——內(nèi)部阻力因子，1/m。

（3）物質(zhì)運輸模型。由于在進行催化還原過程中的流動介質(zhì)存在混合情況，因此選用混合物物質(zhì)輸運模型進行模擬。通過求解描述每種組成物質(zhì)的對流、擴散和反應(yīng)源的守恒方程來模擬混合和輸運，可以模擬多種同時發(fā)生的化學(xué)反應(yīng)。當(dāng)選擇解化學(xué)物質(zhì)的守恒方程時，通過第i種物質(zhì)的對流擴散方程預(yù)估每種物質(zhì)的質(zhì)量分?jǐn)?shù)Yi守恒方程采用以下的通用形式：

3 模擬結(jié)果分析

（1）標(biāo)準(zhǔn)偏差計算分析。隨著對氮氧化合物排放的控制及人們對SCR系統(tǒng)性能要求的不斷提高，目前我們主要通過偏差系數(shù)的大小來衡量SCR系統(tǒng)的設(shè)計性能。偏差系數(shù)的計算公式如下式：

式中：Cv——標(biāo)準(zhǔn)偏差系數(shù)；σ——標(biāo)準(zhǔn)偏差；x——平均值。

在進行數(shù)值模擬計算時，一般常用鍋爐最大連續(xù)運行工況時，SCR反應(yīng)器內(nèi)第一層催化劑上面的速度及NH3濃度分布等偏差(Cv值)是否達到標(biāo)準(zhǔn)來考核系統(tǒng)性能的優(yōu)劣。

（2）數(shù)值模擬結(jié)果分析。根據(jù)SCR脫硝裝置的幾何尺寸及運行參數(shù)等，對初步設(shè)計的方案進行模擬及分析，經(jīng)過反復(fù)的模擬和調(diào)整，得到優(yōu)化后的導(dǎo)流板布置方式和技術(shù)指標(biāo)。經(jīng)過校驗和對比，發(fā)現(xiàn)原有煙道導(dǎo)流板流場分布基本合理，建議改進方案如下：將導(dǎo)流板3間距減少一半，進行加密；調(diào)整導(dǎo)流板4角度；調(diào)整導(dǎo)流板5兩側(cè)各五塊導(dǎo)流板長度；加密導(dǎo)流板6，間隔減少一倍，并減少噴氨側(cè)前五塊導(dǎo)流板長度。其余原有導(dǎo)流板及整流格柵可以利舊，保持不變。研究結(jié)果表明，該建議設(shè)計方案能取得很好地?zé)煔饬鲃臃植己途S持系統(tǒng)較低的壓降。實驗在100%BMCR工況下測定，圖3為優(yōu)化方案100%BMCR工況下第一層催化劑入口截面速度分布，將該截面的各點速度大小導(dǎo)出代入公式（4），計算速度分布的標(biāo)準(zhǔn)偏差系數(shù)為9.08%，滿足技術(shù)指標(biāo)關(guān)于“速度最大偏差為±10%”的要求。

圖4為優(yōu)化方案100%BMCR工況下SCR反應(yīng)器內(nèi)第一層催化劑入處的氨濃度分布圖，將該截面的各點氨濃度分布大小導(dǎo)出代入公式（4），求得標(biāo)準(zhǔn)偏差系數(shù)為4.81%，符合技術(shù)指標(biāo)要求。

圖3 100%BMCR工況下第一催化劑床層速度云圖

圖4 100%BMCR工況下催化劑第一床層氨濃度云圖

圖5 為建議優(yōu)化方案100%BMCR第一層催化劑上游速度偏角余弦值分布云圖，將該截面的各點余弦值分布大小導(dǎo)出代入公式（4），求得標(biāo)準(zhǔn)偏差系數(shù)為0.20%，速度偏角平均值為3.51o（即cos（θ）平均值為0.99813，θ=3.51o），滿足技術(shù)指標(biāo)關(guān)于“速度偏角小于±5°”的要求。

圖5 100%BMCR工況下第一層催化劑入口速度偏角

4 結(jié)語

通過對華電新疆發(fā)電有限公司昌吉熱電廠1、2號機組煙氣脫硝超低排放改造工程的CFD研究分析發(fā)現(xiàn)，原始方案模擬結(jié)果顯示流場分布中速度場大小和方向基本合理，能夠達到第一層催化劑入口處關(guān)于煙氣速度和煙氣入射角度不滿足相關(guān)要求。改造方案模擬結(jié)果顯示第一層催化劑入口處的速度分布更加均勻，速度標(biāo)準(zhǔn)偏差系數(shù)降至9.08%；氨濃度分布更均勻，氨濃度標(biāo)準(zhǔn)偏差系數(shù)下降至4.81%，達到了技術(shù)指標(biāo)關(guān)于“氨濃度偏差小于±5%”的要求。通過優(yōu)化設(shè)計后系統(tǒng)計算數(shù)據(jù)與物理模擬測試的結(jié)果相吻合，實現(xiàn)了優(yōu)化設(shè)計的目標(biāo)要求。該項目導(dǎo)流板方案的成功設(shè)計對于其他機組脫硝SCR反應(yīng)器的設(shè)計具有很好的應(yīng)用價值。