閆 睿 ，鄧雨生 ，梁樹雄 ，王世川 ，鄒斌斌 ，樓 波

（1.華南理工大學(xué) 電力學(xué)院，廣東 廣州 510641；2.中國石化集團(tuán)茂名石油化工有限公司，廣東 茂名525000）

循環(huán)流化床鍋爐的燃燒效率以及NOx的排放和一次風(fēng)量、二次風(fēng)量、燃料量、煤料種類等[1]有關(guān)。 電廠面臨著提高發(fā)電效率以及降低廢氣排放的兩大任務(wù)。建立鍋爐燃燒預(yù)測模型，對(duì)不同工況下鍋爐燃燒效率以及NOx排放進(jìn)行預(yù)測以及優(yōu)化調(diào)參， 用來指導(dǎo)鍋爐的經(jīng)濟(jì)運(yùn)行成為當(dāng)下研究的重點(diǎn)和熱點(diǎn)之一。近年來，計(jì)算機(jī)技術(shù)不斷發(fā)展，同時(shí)，人工智能算法也在工業(yè)中很多領(lǐng)域得到了發(fā)展， 這使得燃燒優(yōu)化運(yùn)行變?yōu)榱丝尚小?/p>

很多學(xué)者對(duì)傳統(tǒng)燃煤鍋爐以及循環(huán)流化床鍋爐燃燒系統(tǒng)進(jìn)行了建模， 并對(duì)不同鍋爐參數(shù)進(jìn)行了預(yù)測。葉海文[2]、馬寶萍[3]、康之霞[4]等人對(duì)循環(huán)流化床鍋爐以人工智能算法為基礎(chǔ)進(jìn)行建模， 成功預(yù)測了鍋爐床溫， 證明了人工智能算法技術(shù)可以在鍋爐燃燒系統(tǒng)中得到應(yīng)用；吳恒運(yùn)[5]、彭道剛[6]、曹慶才[7]等人以神經(jīng)網(wǎng)絡(luò)和其他人工智能算法為基礎(chǔ)對(duì)傳統(tǒng)煤粉鍋爐燃燒系統(tǒng)不同參數(shù)進(jìn)行了建模并做出了優(yōu)化；李鵬輝[8]采用了基于BP 神經(jīng)網(wǎng)絡(luò)的方法對(duì)燃煤鍋爐的NOx排放進(jìn)行了預(yù)測；王凱[9]、白建云[10]等人使用BP 神經(jīng)網(wǎng)絡(luò)方法對(duì)CFB 鍋爐NOx排放量進(jìn)行了預(yù)測建模；李國強(qiáng)[11]、李少華[12]、王杰[13]、王耿杰[14]等人利用不同的人工智能算法對(duì)鍋爐燃燒系統(tǒng)進(jìn)行了優(yōu)化控制。

盡管前人利用神經(jīng)網(wǎng)絡(luò)和其他算法對(duì)燃煤鍋爐效率和NOx排放進(jìn)行了研究， 但是對(duì)循環(huán)流化床鍋爐效率和污染物的排放同時(shí)進(jìn)行預(yù)測并優(yōu)化運(yùn)行參數(shù)研究不多， 且本文中在優(yōu)化時(shí)的目標(biāo)函數(shù)中用飛灰含碳量值代替了傳統(tǒng)的鍋爐效率，簡化了模型。本文中選取的12 個(gè)輸入?yún)?shù)均為某廠循環(huán)流化床鍋爐實(shí)際運(yùn)行參數(shù)，預(yù)測輸出兩個(gè)，分別是飛灰含碳量和NOx排放量， 建立了GA-BP 鍋爐燃燒預(yù)測模型，在預(yù)測模型的基礎(chǔ)上， 使用遺傳算法對(duì)鍋爐參數(shù)進(jìn)行了優(yōu)化，目的在于使鍋爐高效低污染運(yùn)行。

1 研究對(duì)象

研究對(duì)象為美國福斯特·惠勒公司制造的緊湊型CFB 鍋爐，額定蒸發(fā)量410 t/h，型號(hào)為FW-410-12.5-530-CCFB，采用自然循環(huán)，單汽包、單爐膛、平衡通風(fēng)、半露天的布置。 單爐膛整體為“褲腳”型，從下到上爐膛結(jié)構(gòu)可依次分為密相區(qū)、 過渡區(qū)和稀相區(qū)。密相區(qū)下部水冷壁上，開有許多循環(huán)流化床鍋爐所需的特殊門孔，其中包括前墻4 個(gè)給煤口，4 個(gè)石灰石噴入口， 前后墻共計(jì)13 個(gè)一次風(fēng)口和16 個(gè)二次風(fēng)口，另一部風(fēng)從底部布風(fēng)板風(fēng)帽進(jìn)入；過渡區(qū)至稀相區(qū)爐膛內(nèi)部有七道防磨梁； 密相區(qū)鍋爐尾部出口位置靠近爐膛兩側(cè)墻，為矩形，兩個(gè)旋風(fēng)分離器把爐膛出口與后墻回料口連接起來。 本循環(huán)流化床簡單物理結(jié)構(gòu)模型如圖1 所示。

本研究樣本容量、 樣本工況的遍歷性直接決定模型的真實(shí)性和準(zhǔn)確性， 故在收集樣本數(shù)據(jù)時(shí)選擇該電廠中額定負(fù)荷條件下450 組工況數(shù)據(jù)， 其中400 組工況為訓(xùn)練集，50 組工況為測試集，以確保數(shù)據(jù)量豐富且能夠反映實(shí)際工況。

圖1 循環(huán)流化床鍋爐物理結(jié)構(gòu)示意圖

2 CFB 鍋爐燃燒系統(tǒng)預(yù)測模型

2.1 飛灰含碳量以及NOx 排放的預(yù)測模型

影響鍋爐效率的因素很多， 但是較為重要的參數(shù)是飛灰含碳量， 可以用這個(gè)參數(shù)作為衡量鍋爐燃燒效率的一個(gè)重要指標(biāo)，NOx是礦物燃料在空氣中高溫燃燒過程所產(chǎn)生的， 鍋爐燃燒過程中NOx的排放和煙氣含氧量、風(fēng)量、給煤量等鍋爐運(yùn)行的參數(shù)有關(guān)。 在鍋爐效率提高的同時(shí)，考慮降低NOx的排放，是循環(huán)流化床鍋爐燃燒的控制目標(biāo)。

BP 神經(jīng)網(wǎng)絡(luò)是一種多層的前饋的神經(jīng)網(wǎng)絡(luò)，特點(diǎn)在于誤差能夠逆向傳播進(jìn)行訓(xùn)練， 具有很強(qiáng)的多維函數(shù)映射能力， 循環(huán)流化床鍋爐則具有多變量以及多耦合的特點(diǎn)，這一特點(diǎn)恰好能得到解決。本文燃燒預(yù)測模型中，輸入?yún)?shù)m 共計(jì)12 個(gè)，分別是主蒸汽壓力、給水溫度、主蒸汽溫度、一次風(fēng)溫度、一次風(fēng)風(fēng)壓、二次風(fēng)風(fēng)壓、二次風(fēng)溫度、一次風(fēng)總量、二次風(fēng)總量、省煤器出口煙溫、煤料量以及固定碳，輸出參數(shù)n=2 為飛灰含碳量和NOx排放， 燃燒系統(tǒng)神經(jīng)網(wǎng)絡(luò)模型如圖2 所示。

2.2 神經(jīng)網(wǎng)絡(luò)模型的遺傳算法優(yōu)化

雖然BP 神經(jīng)網(wǎng)絡(luò)在化工、能源等諸多領(lǐng)域應(yīng)用廣泛， 但同時(shí)存在著學(xué)習(xí)速度慢以及容易陷入局部最優(yōu)的缺點(diǎn)， 神經(jīng)網(wǎng)絡(luò)中的隱含層的確定可以由經(jīng)驗(yàn)公式作為參考（在本文中 m=12，n=2，a 可以取1~9 的常數(shù)）， 但是并不能作為絕對(duì)的依據(jù)， 模型預(yù)測準(zhǔn)確和網(wǎng)格訓(xùn)練成功依賴于整個(gè)網(wǎng)絡(luò)結(jié)構(gòu)的確定、初始的閾值和權(quán)值。遺傳算法可以對(duì)上述問題解決， 遺傳算法是模擬自然界選擇的一種優(yōu)化算法， 神經(jīng)網(wǎng)絡(luò)的權(quán)值初始值和閾值初始值可以由遺傳算法進(jìn)行優(yōu)化，模型預(yù)測也會(huì)較為精確。遺傳算法對(duì)神經(jīng)網(wǎng)絡(luò)權(quán)值和閾值的優(yōu)化如圖3 所示。

圖2 鍋爐燃燒系統(tǒng)神經(jīng)網(wǎng)絡(luò)模型示意圖

圖3 遺傳算法對(duì)神經(jīng)網(wǎng)絡(luò)的權(quán)值和閾值的優(yōu)化過程示意圖

2.3 GA-BP 神經(jīng)網(wǎng)絡(luò)模型仿真

在Matlab 軟件中進(jìn)行模擬仿真之前要對(duì)輸入?yún)?shù)以及輸出參數(shù)進(jìn)行歸一化處理， 使得所有參數(shù)的值都在[0，1]之間，原因在于這些參數(shù)在數(shù)量級(jí)上差異很大，歸一化處理后讀取數(shù)據(jù)進(jìn)行仿真操作，遺傳算法優(yōu)化時(shí)選擇輪盤賭法以及最優(yōu)個(gè)體保留法，輸入節(jié)點(diǎn)至隱含節(jié)點(diǎn)選擇tansig 函數(shù)， 隱含節(jié)點(diǎn)至輸出節(jié)點(diǎn)選擇purelin 函數(shù)，此外訓(xùn)練誤差為0.001，學(xué)習(xí)速率為0.1。通過上述經(jīng)驗(yàn)公式不斷測試不同隱含節(jié)點(diǎn)下的誤差情況， 發(fā)現(xiàn)隱含層為22 時(shí)誤差最小，于是確定該BP 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)為12-22-2；遺傳算法優(yōu)化神經(jīng)網(wǎng)絡(luò)中，種群大小規(guī)模為100，最大遺傳代數(shù)為100，交叉概率為0.95，變異概率為0.1，適應(yīng)度函數(shù)是誤差平方和的倒數(shù)，如圖4 所示。可以看出適應(yīng)度值在不斷增加， 且在代數(shù)70 后趨于穩(wěn)定，即誤差平方和不斷減小，不斷趨于平穩(wěn)，同時(shí)說明遺傳算法的代數(shù)等設(shè)置合理。

圖4 適應(yīng)度函數(shù)值和代數(shù)之間的關(guān)系圖

預(yù)測仿真結(jié)果如圖5、圖6 所示，仿真結(jié)果參數(shù)飛灰含碳量和NOx排放均進(jìn)行了相對(duì)誤差對(duì)比和絕對(duì)誤差對(duì)比， 飛灰含碳量預(yù)測輸出和期望輸出最大絕對(duì)誤差為0.31%，相對(duì)誤差均在5%以內(nèi)，其中最大值為4.8%， 平均相對(duì)誤差3.3%；NOx排放量預(yù)測輸出和期望輸出絕對(duì)誤差最大值為3.2 mg/m3，相對(duì)誤差均在5%以內(nèi)，其中最大值為4.7%，平均相對(duì)誤差3.2%。 由此可見該GA-BP 網(wǎng)絡(luò)模型擬合度很好，能夠較準(zhǔn)確地預(yù)測飛灰含碳量以及NOx排放量，為下一步鍋爐燃燒參數(shù)的優(yōu)化打下了基礎(chǔ)。

3 鍋爐燃燒參數(shù)優(yōu)化的實(shí)現(xiàn)

3.1 設(shè)計(jì)優(yōu)化函數(shù)

上述預(yù)測模型的建立為本節(jié)優(yōu)化打下了基礎(chǔ)，循環(huán)流化床鍋爐多目標(biāo)優(yōu)化的目的是提高鍋爐燃燒效率和降低NOx排放， 考慮鍋爐熱效率的同時(shí)也要注重NOx的排放，本文采用公式（1）作為目標(biāo)函數(shù)：

式中：D 為飛灰含碳量的歸一化值；CNOx為NOx的歸一化值，所取范圍均為[0，1]，二者都由上述GA-BP 神經(jīng)網(wǎng)絡(luò)燃燒預(yù)測模型得到；a、b 表示權(quán)重， 取決于燃燒優(yōu)化過程中對(duì)燃燒效率和低污染排放的關(guān)注度。

圖5 GA-BP 網(wǎng)絡(luò)飛灰含碳量仿真結(jié)果及其誤差圖

圖6 GA-BP 網(wǎng)絡(luò)NOx 排放量仿真結(jié)果及其誤差圖

由于影響循環(huán)流化床鍋爐效率的主要因素之一是飛灰含碳量， 選擇飛灰含碳量作為尋優(yōu)的評(píng)價(jià)參量代替了傳統(tǒng)的鍋爐效率，這一舉措可以簡化模型，優(yōu)化函數(shù)中的飛灰含碳量參數(shù)可以直接使用上述GA-BP 神經(jīng)網(wǎng)絡(luò)模型的輸出值，循環(huán)流化床鍋爐燃燒系統(tǒng)優(yōu)化如圖7 所示。

3.2 遺傳算法優(yōu)化

圖7 循環(huán)流化床鍋爐燃燒優(yōu)化示意圖

對(duì)于上述設(shè)計(jì)的優(yōu)化函數(shù)， 本文采用遺傳算法進(jìn)行參數(shù)尋優(yōu)。 同時(shí)兼顧爐效和低污染，故式（1）中取權(quán)重a=b=0.5，設(shè)置整個(gè)種群規(guī)模數(shù)量100，最大進(jìn)化代數(shù)為50，變異概率為0.1，交叉概率為0.4。 選擇GA-BP 神經(jīng)網(wǎng)絡(luò)中的9 個(gè)輸入?yún)?shù)作為燃燒優(yōu)化參量（由于試驗(yàn)時(shí)以及鍋爐日常運(yùn)行時(shí)均為額定負(fù)荷運(yùn)行，故主蒸汽壓力、主蒸汽溫度、給水溫度不在優(yōu)化運(yùn)行參數(shù)內(nèi)），優(yōu)化過程中的數(shù)據(jù)均為歸一化以后的數(shù)據(jù)，優(yōu)化結(jié)束后轉(zhuǎn)換為實(shí)際值，9 個(gè)輸入?yún)⒘康膬?yōu)化區(qū)域?yàn)樵摌颖局暗?0 個(gè)樣本和未來5 個(gè)樣本所構(gòu)成的集合。 圖8 為遺傳算法優(yōu)化循環(huán)流化床鍋爐運(yùn)行參數(shù)示意圖。

選取50 組測試集樣本數(shù)據(jù)進(jìn)行優(yōu)化，結(jié)果對(duì)比如圖9 所示。可以看出，雖然有部分的優(yōu)化結(jié)果并不完美（飛灰含碳量或NOx排放的優(yōu)化值大于實(shí)際輸出值），但是從整體上分析，循環(huán)流化床飛灰含碳量和NOx排放都得到了優(yōu)化， 可以通過調(diào)整系數(shù)a、b 來實(shí)現(xiàn)飛灰含碳量和NOx排放的不同權(quán)重。

4 工業(yè)試驗(yàn)

圖8 遺傳算法優(yōu)化循環(huán)流化床鍋爐運(yùn)行參數(shù)示意圖

為了能夠驗(yàn)證優(yōu)化的結(jié)果， 選取一組高飛灰含碳量（8.79%）和高 NOx排放量（74.1 mg/m3）的優(yōu)化結(jié)果在本循環(huán)流化床鍋爐進(jìn)行了工業(yè)試驗(yàn)， 在額定負(fù)荷的條件下，運(yùn)行參數(shù)為優(yōu)化后的輸入?yún)⒘浚玫斤w灰含碳量以及NOx排放量的數(shù)據(jù)。 通過表1 可以看出：一次風(fēng)風(fēng)溫、二次風(fēng)風(fēng)溫、總一次風(fēng)、總二次風(fēng)均有下降，其中總一次風(fēng)下降比較明顯；一次風(fēng)、二次風(fēng)風(fēng)壓有小幅上升?？傄淮物L(fēng)風(fēng)量下降后，會(huì)導(dǎo)致爐內(nèi)密相區(qū)溫度升高，間接造成燃燒效率提高；總二次風(fēng)量下降，導(dǎo)致爐內(nèi)氧量減少，NOx生成量也減少。

圖9 飛灰含碳量以及NOx 排放原期望輸出與優(yōu)化輸出對(duì)比圖

表1 優(yōu)化前后結(jié)果以及和試驗(yàn)值的對(duì)比

5 結(jié)論

（1）建立了以飛灰含碳量和NOx排放量為組合的CFB 鍋爐燃燒預(yù)測模型，以飛灰含碳量代替循環(huán)流化床鍋爐效率作為參數(shù)簡化了模型， 以飛灰含碳量和NOx排放量的線性組合為優(yōu)化目標(biāo)。

（2）對(duì)該 410 t/h 循環(huán)流化床（CFB）鍋爐的燃燒系統(tǒng)進(jìn)行了建模， 優(yōu)化后一組試驗(yàn)運(yùn)行結(jié)果飛灰含碳量從 8.79%降低至 7.98%，NOx排放量降低了15.4%，實(shí)現(xiàn)了高效低污染燃燒的目標(biāo)。

（3）該優(yōu)化可以在保證鍋爐發(fā)熱量的基礎(chǔ)上，以飛灰含碳量和NOx排放量的線性組合為優(yōu)化目標(biāo)，使用遺傳算法對(duì)所建立的循環(huán)流化床鍋爐燃燒模型輸入量中的可調(diào)變量在可變化的范圍內(nèi)進(jìn)行尋優(yōu)，得到參數(shù)最終的組合， 使得最后優(yōu)化的目標(biāo)函數(shù)值最小，即可以滿足鍋爐效率達(dá)到最大的同時(shí)，NOx的排放量最小， 可以在一定程度上指導(dǎo)該鍋爐的優(yōu)化運(yùn)行，具有實(shí)際意義。