黃鎮(zhèn)柯，金志遠(yuǎn)，朱恒毅，劉宇濃，譚 鵬

（1.廣東省能源集團(tuán)有限公司沙角C電廠，廣東 東莞 523900；2.華中科技大學(xué)能源與動力工程學(xué)院煤燃燒國家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430074）

0 引言

為控制氮氧化物（NOx）排放，選擇性催化還原（SCR）技術(shù)在燃煤電站獲得廣泛應(yīng)用［1，2］。電站SCR脫硝系統(tǒng)通過控制噴氨閥門開度來調(diào)節(jié)出口NOx的排放濃度，是一個典型的大慣性延遲系統(tǒng)。過量噴氨會造成氨逃逸量增加，與煙氣中的硫發(fā)生反應(yīng)產(chǎn)生硫酸氫銨，導(dǎo)致下游空預(yù)器堵塞，威脅電廠安全生產(chǎn)；噴氨不足會造成SCR 脫硝系統(tǒng)脫硝效率降低，影響電廠運(yùn)行環(huán)保性［3，4］。然而隨著可再生能源入網(wǎng)比例的升高［5］，燃煤機(jī)組的調(diào)峰定位越來越清晰［6-8］，頻繁的機(jī)組負(fù)荷大范圍變化使得采用傳統(tǒng)PID 控制不能滿足SCR脫硝系統(tǒng)的高效安全運(yùn)行［9，10］，因此，研究開發(fā)更精確、更快速的噴氨控制方法具有重要意義［11-13］。

建立準(zhǔn)確可靠的SCR 模型是控制的基礎(chǔ)，為實(shí)現(xiàn)SCR脫硝系統(tǒng)的精確噴氨控制必須掌握系統(tǒng)的動態(tài)特性，建立合適的動態(tài)模型。目前，SCR脫硝建模方法可分為基于催化化學(xué)反應(yīng)機(jī)理的機(jī)理建模［14］和基于運(yùn)行數(shù)據(jù)的黑箱建模［15］，其中黑箱建模方法主要是利用統(tǒng)計方法或機(jī)器學(xué)習(xí)，如神經(jīng)網(wǎng)絡(luò)［16，17］、支持向量機(jī)［18，19］等方法，采用大量運(yùn)行數(shù)據(jù)進(jìn)行模型訓(xùn)練，由于其相對機(jī)理建模具有模型精度高、建模周期短的優(yōu)勢，得到了較多的關(guān)注。

劉吉臻等［20］將多尺度核函數(shù)、滑動窗口法和偏最小二乘法結(jié)合起來，提出了自適應(yīng)多尺度核偏最小二乘法，并將其利用在SCR 脫硝系統(tǒng)穩(wěn)態(tài)建模上。秦天牧等［21］分別使用機(jī)理建模、支持向量機(jī)、反向傳播（back propagation，BP）神經(jīng)網(wǎng)絡(luò)和偏最小二乘法對SCR 脫硝系統(tǒng)進(jìn)行動態(tài)建模，并比較了4 種模型的性能差別。周洪煜等［22］引入動態(tài)結(jié)構(gòu)的RBF 神經(jīng)網(wǎng)絡(luò)在SCR 脫硝建模問題上，并提出了相應(yīng)的噴氨控制系統(tǒng)。

本文以某660 MW亞臨界燃煤機(jī)組SCR脫硝系統(tǒng)為研究對象，利用歷史運(yùn)行數(shù)據(jù)，建立基于多層感知器神經(jīng)網(wǎng)絡(luò)（MLP）的SCR 脫硝系統(tǒng)預(yù)測模型，旨在為SCR噴氨控制器設(shè)計與優(yōu)化提供依據(jù)。

1 研究對象與方法

1.1 研究對象

本文研究對象為某660 MW 亞臨界燃煤機(jī)組，由美國ABB-CE公司設(shè)計生產(chǎn)，類型為亞臨界參數(shù)、單汽包中間再熱、控制循環(huán)鍋爐，鍋爐最大連續(xù)蒸發(fā)量為2 100 t/h，膜式水冷壁，單爐膛四角雙切圓燃燒。SCR系統(tǒng)采用液氨作為還原劑，液氨氣化后經(jīng)稀釋風(fēng)機(jī)稀釋，由噴氨格柵送入SCR脫硝反應(yīng)器入口，與煙氣混合后，在催化劑作用下與煙氣中的NOx發(fā)生選擇性催化還原反應(yīng)。該脫硝系統(tǒng)采用V2O5-WO3/TiO2型的蜂窩式催化劑，催化劑分3層布置反應(yīng)器內(nèi)。

SCR 脫硝系統(tǒng)出口NOx值受很多因素影響，主要因素有負(fù)荷、煙氣流量、給煤量、SCR 脫硝系統(tǒng)入口NOx值、噴氨閥門開度、SCR 脫硝系統(tǒng)溫度、煙氣含氧量等［23，24］，其中控制變量為噴氨閥門開度，最主要的外部擾動包括SCR 脫硝系統(tǒng)入口NOx值以及煙氣流量、SCR 脫硝系統(tǒng)溫度、催化劑活性等。催化劑活性是一個長周期衰變過程，難以監(jiān)測；煙氣流量、SCR 脫硝系統(tǒng)溫度與負(fù)荷密切相關(guān)。綜合考慮以上因素，本文所建立的SCR 脫硝系統(tǒng)預(yù)測模型輸入為：當(dāng)前及過去時刻的①負(fù)荷，②SCR 脫硝系統(tǒng)入口NOx值，③噴氨閥門開度，④SCR 脫硝系統(tǒng)出口NOx值；模型輸出為：下一時刻的①SCR脫硝系統(tǒng)出口NOx值。本文共采集了1個月的運(yùn)行數(shù)據(jù)，前一周數(shù)據(jù)用于訓(xùn)練，后三周數(shù)據(jù)用于驗(yàn)證模型穩(wěn)定性。數(shù)據(jù)采樣周期為10 s，訓(xùn)練數(shù)據(jù)集總長度為69 121，各參數(shù)變化范圍如表1所示。

表1 各參數(shù)單位及變化范圍Table 1 Each parameter unit and variation range

1.2 多層感知器神經(jīng)網(wǎng)絡(luò)

多層感知器神經(jīng)網(wǎng)絡(luò)［25，26］（MLP）的概念最先由Rosenblatt 于1952 年提出。多層感知器神經(jīng)網(wǎng)絡(luò)是人工神經(jīng)網(wǎng)絡(luò)的一種，一般由一個輸入層、幾個隱藏層和一個輸出層組成。人工神經(jīng)網(wǎng)絡(luò)［27，28］由一組突觸連接構(gòu)成，每個連接都有對應(yīng)的權(quán)重，所有輸入經(jīng)過加權(quán)之后與偏置相加，通過激活函數(shù)來轉(zhuǎn)化輸出，使得其幅度范圍縮小到有限值。其數(shù)學(xué)表達(dá)式為：

式（1）中，x為輸入值，和為第1層的權(quán)值，σ1為第1層的激活函數(shù)，神經(jīng)網(wǎng)絡(luò)的輸出是最后一層的值zm。

1.3 基于多層感知機(jī)的SCR脫硝系統(tǒng)預(yù)測建模

為了避免各輸入?yún)?shù)由于數(shù)量級和變化范圍不同所導(dǎo)致的建模效果不佳，對輸入?yún)?shù)進(jìn)行歸一化處理。引入EarlyStopping［29］方法避免過擬合，該方法可根據(jù)指定的評價標(biāo)準(zhǔn)及時停止模型的訓(xùn)練，當(dāng)連續(xù)n代訓(xùn)練模型的評價低于之前訓(xùn)練的最好結(jié)果時，立即停止模型的訓(xùn)練，這樣就可以避免繼續(xù)訓(xùn)練導(dǎo)致誤差繼續(xù)增大，激活函數(shù)選用線性整流函數(shù)（ReLU），采用Adam 優(yōu)化器對MLP 神經(jīng)網(wǎng)絡(luò)模型進(jìn)行訓(xùn)練［30-34］。Adam 是一種基于低階矩自適應(yīng)估計的一階梯度隨機(jī)目標(biāo)函數(shù)優(yōu)化算法，該方法易于實(shí)現(xiàn)，計算效率高，內(nèi)存需求少，非常適合于數(shù)據(jù)或參數(shù)較大的問題，同時也適用于非平穩(wěn)目標(biāo)和非常嘈雜或稀疏梯度的問題。

使用均方誤差（MSE）評價標(biāo)準(zhǔn)對模型進(jìn)行評價，其公式為：

2 結(jié)果與討論

2.1 基于多層感知機(jī)的SCR脫硝系統(tǒng)預(yù)測建模

采用基于多層感知機(jī)對SCR脫硝系統(tǒng)預(yù)測建模需對時間步長、網(wǎng)絡(luò)層數(shù)、隱藏層節(jié)點(diǎn)數(shù)及訓(xùn)練批次數(shù)據(jù)量大小進(jìn)行優(yōu)化。

時間步長調(diào)參結(jié)果如圖1所示，計算時模型采用2個隱藏層，每個隱藏層節(jié)點(diǎn)數(shù)為5。

由圖1 可以看出，時間步長設(shè)置太小或太大都會導(dǎo)致MSE較大，當(dāng)時間步長為3時，訓(xùn)練集和測試集的MSE 較小。在保證模型精度的情況下，考慮計算所需時間和資源，最終選擇模型的時間步長為3。

圖1 時間步長調(diào)參結(jié)果Fig.1 Parameter adjustment results of time step-length

神經(jīng)網(wǎng)絡(luò)隱藏層數(shù)量優(yōu)化如圖2 所示，模型的隱藏層節(jié)點(diǎn)數(shù)為5，時間步長為3。

圖2 隱藏層個數(shù)調(diào)參結(jié)果Fig.2 Parameter adjustment results of the number of hidden layers

由圖2 可以看出，隱藏層個數(shù)除了1 層之外，其他2至5層MSE誤差均相差不大，在保證模型精度的情況下，考慮計算所需時間和資源，最終選擇隱藏層個數(shù)為2。

神經(jīng)網(wǎng)絡(luò)每個隱藏層節(jié)點(diǎn)數(shù)優(yōu)化如圖3所示，隱藏層個數(shù)為2個，時間步長為3。由圖3可以看出，每個隱藏層個數(shù)為1個和2個時，模型效果很差，驗(yàn)證集MSE誤差分別為31和29，當(dāng)隱藏層節(jié)點(diǎn)數(shù)在3個至5個時，模型的效果提升明顯，驗(yàn)證集MSE 誤差分別為0.189、0.191和0.193，因此確定最終每個隱藏層節(jié)點(diǎn)數(shù)為3個。

圖3 不同節(jié)點(diǎn)數(shù)調(diào)參結(jié)果Fig.3 Parameter adjustment results of the number of different nodes

經(jīng)過對模型參數(shù)的調(diào)整，最終確定多層感知器結(jié)構(gòu)為：時間步長為3，隱藏層數(shù)為2，隱藏層節(jié)點(diǎn)數(shù)為3，激活函數(shù)選用ReLU，其訓(xùn)練過程如圖4 所示，可以看出，在第12 代-15 代的評價指標(biāo)已不下降，并在50 代時停止了訓(xùn)練。

圖4 MLP訓(xùn)練過程Fig.4 Training process of MLP

圖5顯示了模型在測試集上的預(yù)測效果和誤差分布，可以看出，預(yù)測值與實(shí)際值基本吻合，相關(guān)系數(shù)達(dá)到了0.99，絕對誤差在0 附近波動，波動范圍不大，說明訓(xùn)練出來的模型有很好的預(yù)測能力。

圖5 SCR出口NOx值預(yù)測效果Fig.5 Prediction effect of NOx concentration at SCR outlet

圖6為測試集的預(yù)測值與實(shí)際值之間誤差的頻率分布圖，可以很明顯地看出測試集的誤差呈期望值為0 的正態(tài)分布，約有78.78%測試集樣本的誤差絕對值位于±0.5 mg/Nm3之間，有18.5%測試集樣本的誤差絕對值位于±0.5 mg/Nm3至±1 mg/Nm3之間。

圖6 MLP模型訓(xùn)練誤差頻率Fig.6 Training error frequency of MLP model

2.2 基于多層感知機(jī)的SCR脫硝系統(tǒng)模型階躍測試

為進(jìn)一步驗(yàn)證模型的效果，針對多層感知機(jī)模型還做了階躍測試，以驗(yàn)證不同特征變化下，SCR 出口NOx值的響應(yīng)情況，具體階躍效果如圖7所示。

圖7 展示了SCR 噴氨閥門開度特征階躍效果，針對不同的SCR 噴氨閥門開度階躍幅度，模型均做出合理且平滑的響應(yīng)。SCR噴氨閥門開度階躍+0.5%，SCR出口NOx值降低-0.36 mg/Nm3，SCR 噴氨閥門開度階躍+1%，SCR 出口NOx值降低0.72 mg/Nm3，SCR 噴氨閥門開度階躍+1.5%，SCR 出口NOx值降低0.36 mg/Nm3，SCR 噴氨閥門開度階躍+0.5%，SCR 出口NOx值降低1.09 mg/Nm3，SCR 噴氨閥門開度階躍-0.5%，SCR 出口NOx值升高0.37 mg/Nm3，SCR 噴氨閥門開度階躍-1%，SCR 出口NOx值升高0.74 mg/Nm3，SCR 噴氨閥門開度階躍-1.5%，SCR出口NOx值升高1.1 mg/Nm3。

圖7 MLP模型特征階躍響應(yīng)Fig.7 Features step response of MLP model

3 結(jié)論

本文以某660 MW 亞臨界燃煤鍋爐為研究對象，建立了基于多層感知器神經(jīng)網(wǎng)絡(luò)的SCR脫硝系統(tǒng)預(yù)測模型。結(jié)果表明，MLP 訓(xùn)練模型在測試集上的均方誤差為0.189，相關(guān)系數(shù)達(dá)到了0.99，同時針對SCR 噴氨閥門開度特征的階躍變化也作出了合理的響應(yīng)，所建立的模型可用于SCR噴氨控制器設(shè)計與優(yōu)化。