中電華創(chuàng)電力技術(shù)研究有限公司 陳小強(qiáng) 俞衛(wèi)新 東南大學(xué)能源與環(huán)境學(xué)院 單鑫晨 趙伶玲
根據(jù)國家的雙碳目標(biāo),火電機(jī)組將逐漸由發(fā)電主力轉(zhuǎn)變?yōu)檎{(diào)峰調(diào)頻的主力。而在負(fù)荷調(diào)整過程中,如何能滿足國家燃煤電廠超低排放環(huán)保要求將NOx 排放限值控制在50mg/m3以下,與SCR 系統(tǒng)內(nèi)部運(yùn)行狀況息息相關(guān),因此有必要對(duì)不同負(fù)荷時(shí)SCR 內(nèi)部的具體運(yùn)行情況進(jìn)行研究。
目前對(duì)SCR 的研究有數(shù)值模擬及實(shí)驗(yàn)兩種主要手段。數(shù)值模擬在SCR 流場(chǎng)診斷及流場(chǎng)優(yōu)化方面使用較多,但在負(fù)荷調(diào)整與SCR 反應(yīng)器運(yùn)行情況的關(guān)系方面應(yīng)用還較少。凌忠錢等人[1]對(duì)于不同負(fù)荷下SCR 運(yùn)行狀況對(duì)某300MW 鍋爐的SCR 系統(tǒng)進(jìn)行研究,分析了導(dǎo)流板布置對(duì)SCR 反應(yīng)器內(nèi)流場(chǎng)分布均勻性的影響,研究了不同布置方案下SCR 反應(yīng)器內(nèi)氨氮物質(zhì)的量比(簡稱氨氮比)分布的規(guī)律;朱天宇等人[2]對(duì)某600MW 超臨界燃煤機(jī)組SCR 系統(tǒng)進(jìn)行反應(yīng)器內(nèi)煙氣流動(dòng)均勻性和飛灰沉積的數(shù)值模擬,開展了煙道內(nèi)導(dǎo)流板布置形式和導(dǎo)流板結(jié)構(gòu)對(duì)SCR 反應(yīng)器內(nèi)流場(chǎng)以及飛灰沉積的影響研究。
鄒紅果[3]等人采用數(shù)值模擬對(duì)某660MW 機(jī)組SCR 系統(tǒng)流場(chǎng)進(jìn)行全面診斷,并通過數(shù)值模擬提出導(dǎo)流板優(yōu)化方案,優(yōu)化了噴氨性能;牛彩偉等人[4]利用數(shù)值模擬研究了某600MW 燃煤電站鍋爐SCR系統(tǒng)不同工況下內(nèi)部流場(chǎng)和濃度場(chǎng)的變化情況,研究表明低負(fù)荷下流場(chǎng)不均勻度升高導(dǎo)致SCR 脫硝效率降低;郭江源等人[5]采用控制變量的方式對(duì)SCR系統(tǒng)進(jìn)行試驗(yàn),分別對(duì)5臺(tái)火電廠機(jī)組在不同負(fù)荷條件下煙氣的溫度、氧含量、污染物含量、脫硝效率以及SO2/SO3轉(zhuǎn)化率等進(jìn)行測(cè)量。
本文的研究對(duì)象為某電廠600MW 超臨界機(jī)組SCR 脫硝系統(tǒng),包括SCR 反應(yīng)器入口煙道和SCR反應(yīng)器。
根據(jù)煙道及SCR 反應(yīng)器設(shè)計(jì)圖紙及現(xiàn)場(chǎng)測(cè)繪,本文按照1:1比例建立了煙道及SCR 反應(yīng)器含導(dǎo)流板的計(jì)算模型(圖1)。該爐SCR 入口煙道及反應(yīng)器內(nèi)共有6處導(dǎo)流板,入口1、2處為直導(dǎo)流板或斜擋板,之后的四處彎道都設(shè)有弧形導(dǎo)流板(3~6號(hào)導(dǎo)流板),數(shù)目分別為3、5、9、4塊,其中第三處同時(shí)還設(shè)有3塊直擋板。此外在圖中第3處與第4處導(dǎo)流板之間、即垂直上升煙道內(nèi)布置了噴氨格柵和靜態(tài)混合器,靜態(tài)混合器由垂直分布的細(xì)管布置而成。試驗(yàn)所使用的測(cè)孔安裝在噴氨前截面和第三層催化劑出口截面,分別有14個(gè),編號(hào)分別為A1~A14和B1~B14。但是由于現(xiàn)場(chǎng)實(shí)際情況,噴氨前截面有個(gè)別測(cè)孔被電廠安裝的分析儀表所擋,可測(cè)量測(cè)孔為11個(gè)。
圖1 SCR 反應(yīng)內(nèi)部幾何結(jié)構(gòu)及測(cè)點(diǎn)位置圖
基本方程:SCR 脫硝系統(tǒng)中煙氣進(jìn)行的是三維湍流流動(dòng),本文計(jì)算選用標(biāo)準(zhǔn)k-ε 湍流模型,其中連續(xù)性方程、動(dòng)量方程和k-ε 方程可統(tǒng)一表達(dá)如式(1),式中:ρ 為氣相密度;Γφ為擴(kuò)散系數(shù);Sφ為源項(xiàng);u、v、w 分別為x、y、z 方向上的速度分量;k、ε 為湍流動(dòng)能和耗散系數(shù)。計(jì)算中φ 取1、u、v、w、k 和ε 時(shí),方程分別代表連續(xù)性方程、各方向的動(dòng)量方程及湍流動(dòng)能和耗散率方程。
模型選擇及計(jì)算細(xì)節(jié):模型入口采用速度進(jìn)口條件,煙氣為不可壓縮流體,流速為入口面平均流速13.08m/s。模型出口為壓力出口條件。對(duì)于反應(yīng)器的催化劑層,采用多孔介質(zhì)模型來模擬。本文采用二階迎風(fēng)差分格式求解湍流動(dòng)能、湍流動(dòng)能耗散率、動(dòng)量方程、對(duì)流擴(kuò)散方程。壓力和速度耦合采用SIMPLE 算法,時(shí)間離散方式選用一階隱式時(shí)間迭代,同時(shí)近壁區(qū)域求解采用標(biāo)準(zhǔn)壁面函數(shù)法。計(jì)算中可通過適當(dāng)調(diào)節(jié)各參數(shù)的松弛因子來滿足收斂的需要,求解器的收斂標(biāo)準(zhǔn)是標(biāo)準(zhǔn)殘差之和小于給定值δ。
實(shí)際工程中,影響SCR 脫硝效果的兩個(gè)性能評(píng)價(jià)指標(biāo)分別為脫硝效率和氨逃逸率。對(duì)于脫硝系統(tǒng)數(shù)值模擬中,煙氣流場(chǎng)和濃度場(chǎng)的均勻程度決定著脫硝效率和是否有過量氨逃逸。為定量研究SCR 脫硝反應(yīng)器煙氣流場(chǎng)和濃度場(chǎng)的均勻程度,本文引入流場(chǎng)不均勻系數(shù)Cv和濃度場(chǎng)不均勻系數(shù)Cρ。
濃度場(chǎng)不均勻系數(shù):NOx 分布不均可導(dǎo)致SCR反應(yīng)器出口局部氨逃逸量過高,造成脫硝副產(chǎn)物硫酸氫銨的生成,從而產(chǎn)生空預(yù)器硫酸氫氨堵灰,嚴(yán)重危害機(jī)組運(yùn)行的安全穩(wěn)定性。另一方面,NOx 分布不均會(huì)導(dǎo)致催化劑各個(gè)區(qū)域損耗程度不均。本文基于NOx 的分布來評(píng)價(jià)脫硝濃度場(chǎng)情況,采用不均勻度(Cρ)分析測(cè)點(diǎn)平面NOx 分布差異程度,NOx分布不均勻度的定義為式中:ρi為測(cè)點(diǎn)煙氣速度,單位mg/Nm3;ρ 為測(cè)點(diǎn)平面平均氣流速度,單位mg/Nm3。
反應(yīng)器整體流場(chǎng):圖2所示為省煤器出口至SCR反應(yīng)器入口煙道及SCR 脫硝反應(yīng)器內(nèi)流場(chǎng)的模擬結(jié)果,可看出從煙氣入口到豎直煙道存在流場(chǎng)不均勻現(xiàn)象,經(jīng)過3號(hào)導(dǎo)流板后流場(chǎng)得以改善,在3號(hào)導(dǎo)流板下方灰斗處存在低速回流區(qū),到5號(hào)導(dǎo)流板下方即SCR 反應(yīng)器時(shí)流場(chǎng)分布已基本均勻。本文進(jìn)一步將450MW、300MW 工況下的模擬結(jié)果與600MW模擬結(jié)果進(jìn)行對(duì)比(圖2),可明顯看出當(dāng)機(jī)組功率從600MW 下 降 到450MW 或300MW 時(shí),SCR 反應(yīng)器入口煙速分別從13.08m/s 下降到10.37m/s 及7.58m/s。與高負(fù)荷工況相比,低負(fù)荷工況下整體速度下降,但流場(chǎng)分布規(guī)律并未發(fā)生較大的改變。
圖2 SCR 反應(yīng)器整體流場(chǎng)圖
噴氨前截面流場(chǎng):噴氨前截面流場(chǎng)計(jì)算結(jié)果如圖3。以600MW 負(fù)荷為例,在3號(hào)導(dǎo)流板的作用下,煙氣由水平煙道進(jìn)入豎直煙道,在噴氨前截面呈現(xiàn)中間流速高、兩側(cè)流速低的分布。其中兩側(cè)平均流速均為13.7m/s,中間平均流速為15.35m/s,主要原因是上升煙道導(dǎo)流板(3號(hào)導(dǎo)流板)的直板段偏向外側(cè)約5°。3號(hào)導(dǎo)流板的導(dǎo)流效果較為明顯,避免了煙氣在較短路程內(nèi)轉(zhuǎn)彎90°后由于離心力的作用貼向上升煙道內(nèi)側(cè)流動(dòng),從而引起外側(cè)產(chǎn)生回流區(qū),流場(chǎng)的不均勻度為8.73%。當(dāng)負(fù)荷降低時(shí),低負(fù)荷工況相比高負(fù)荷工況整體速度下降,流速分布趨勢(shì)并無較大差異。而450MW 和300MW 負(fù)荷下的流場(chǎng)不均勻度分別為8.77%和8.84%。
圖3 SCR 反應(yīng)器噴氨前截面流場(chǎng)
催化劑前截面流場(chǎng):催化劑前截面流場(chǎng)計(jì)算結(jié)果示于圖4??煽闯?00MW 負(fù)荷下的截面速度不均勻度為0.53%,其中該截面左側(cè)區(qū)域速度略低于其他區(qū)域,整體流場(chǎng)較為均勻。450MW 工況下低速區(qū)面積有所增大,截面不均勻度為0.55%。300MW工況低速區(qū)區(qū)域面積進(jìn)一步增大,截面不均勻度為0.57%。總體看來低負(fù)荷條件下低速區(qū)增大,截面不均勻度略有增加。根據(jù)SCR 脫硝反應(yīng)段模擬結(jié)果可得出,隨著負(fù)荷的降低,流速雖有所降低但流場(chǎng)的分布規(guī)律和不均勻度基本相似,同時(shí)反應(yīng)段安裝的幾組導(dǎo)流板也使得整體流場(chǎng)也比較均勻。
圖4 SCR 反應(yīng)器催化劑前截面流場(chǎng)
SCR 系統(tǒng)流場(chǎng)模擬結(jié)果與測(cè)量結(jié)果對(duì)比:本文將現(xiàn)場(chǎng)試驗(yàn)測(cè)得煙速與數(shù)值模擬結(jié)果進(jìn)行分析,由圖5可看出,噴氨前截面煙氣平均流速測(cè)量值為17.08m/s、不均勻度為28.89%,而模擬值為14.46m/s,相對(duì)誤差為15.3%。催化劑出口測(cè)點(diǎn)截面整體平均流速實(shí)測(cè)值為6.33m/s、不均勻度為21.79%,模擬值為7.75m/s,相差為1.42m/s。模擬值與現(xiàn)場(chǎng)試驗(yàn)結(jié)果趨勢(shì)一致,相對(duì)吻合較好。
圖5 噴氨前截面與催化劑出口截面煙氣流速試驗(yàn)與數(shù)值模擬結(jié)果對(duì)比
NOx 濃度測(cè)量結(jié)果:噴氨前和催化劑出口截面NOx 濃度的測(cè)量結(jié)果分別示于圖6和圖7。由圖6可看出,噴氨前截面NOx 含量分布較為均勻,各測(cè)點(diǎn)NOx 濃度差距不大。進(jìn)一步進(jìn)行量化分析,噴氨前截面、催化劑出口截面濃度不均勻度(%)在600MW、450MW、300MW 下分別為10.10/ 34.91、9.05/58.29、11.47/64.14;由圖7可知,SCR 出口截面濃度場(chǎng)的不均勻性大幅度增加。在不同的負(fù)荷下催化劑出口測(cè)點(diǎn)截面沿測(cè)孔方向,NOx 分布都呈現(xiàn)先增大后降低再增大的趨勢(shì)。同時(shí)催化劑出口測(cè)點(diǎn)截面的NOx 濃度不均勻度會(huì)隨負(fù)荷的變小有明顯的增加,450MW 負(fù)荷下和300MW 負(fù)荷下的SCR 出口截面的NOx 濃度不均勻度相比600MW 分別增加了69.0%和83.7%。
圖6 噴氨前截NOx 濃度場(chǎng)測(cè)量結(jié)果
圖7 催化劑出口截面NOx 濃度場(chǎng)測(cè)量結(jié)果
本文采集了試驗(yàn)時(shí)間段的機(jī)組DCS 運(yùn)行數(shù)據(jù)進(jìn)行分析。經(jīng)數(shù)據(jù)分析,SCR 反應(yīng)器進(jìn)、出口氧量與負(fù)荷的關(guān)系如圖8所示,隨著負(fù)荷的增高A 側(cè)入口平均氧量分別降低38.2%和52.5%,出口平均氧量分別降低19.9%和35.6%。B 側(cè)入口平均氧量分別降低33.4%和46.8%,出口平均氧量分別降低41.2%和54.0%。由此可見,該爐在運(yùn)行時(shí)氧含量隨著負(fù)荷增加而降低,當(dāng)負(fù)荷較高時(shí)爐膛溫度較高,爐膛內(nèi)著火條件、煤粉與空氣混合條件較好,燃燒穩(wěn)定,最佳過量空氣系數(shù)較低;而當(dāng)鍋爐低負(fù)荷運(yùn)行時(shí),最佳過量空氣系數(shù)較高。因此高負(fù)荷時(shí)的最佳氧含量要低于低負(fù)荷時(shí)的最佳氧含量。
圖8 SCR 反應(yīng)器進(jìn)出口氧量與負(fù)荷的關(guān)系
SCR 反應(yīng)器入口煙氣溫度與負(fù)荷的關(guān)系如圖9。可看出,隨著負(fù)荷的增加,A 側(cè)入口平均溫度分別增加了11.0%和21.8%,B 側(cè)入口平均溫度分別增加了10.2%和19.5%。入口煙氣溫度隨著負(fù)荷的升高而升高。當(dāng)鍋爐負(fù)荷發(fā)生變化時(shí)燃料的消耗量也發(fā)生變化。負(fù)荷升高時(shí)燃料消耗量增加,煙氣溫度升高。
圖9 SCR 反應(yīng)器入口煙氣溫度與負(fù)荷的關(guān)系
本文以某電廠600MW 機(jī)組SCR 脫硝系統(tǒng)為研究對(duì)象,應(yīng)用數(shù)值模擬、現(xiàn)場(chǎng)試驗(yàn)及DCS 數(shù)據(jù)分析相結(jié)合的方法,研究了負(fù)荷調(diào)整對(duì)SCR 運(yùn)行狀況的影響,得出以下結(jié)論:機(jī)組負(fù)荷發(fā)生調(diào)整時(shí),由于多組導(dǎo)流板的作用,除流速會(huì)隨著負(fù)荷減小而下降外,整體流場(chǎng)的分布規(guī)律和不均勻度基本相似,測(cè)點(diǎn)平面的速度場(chǎng)也比較均勻;負(fù)荷變化時(shí),噴氨前截面的NOx 濃度不均勻度變化不大,但催化劑出口截面的NOx 濃度不均勻度會(huì)隨著負(fù)荷的增大而增加,可能會(huì)導(dǎo)致SCR 反應(yīng)器出口局部的氨逃逸增加;SCR 反應(yīng)器入口和出口氧含量會(huì)隨著負(fù)荷增加而降低,而入口煙氣溫度則隨著負(fù)荷的增加而增加。