基于風(fēng)驅(qū)動算法的鍋爐NOx排放模型優(yōu)化

牛培峰，　趙　振，　馬云鵬，　陳　科，　王丘亞，　趙慶沖

(燕山大學(xué) 工業(yè)計算機控制工程河北省重點實驗室， 河北秦皇島 066004)

?

基于風(fēng)驅(qū)動算法的鍋爐NOx排放模型優(yōu)化

牛培峰，趙振，馬云鵬，陳科，王丘亞，趙慶沖

(燕山大學(xué) 工業(yè)計算機控制工程河北省重點實驗室， 河北秦皇島 066004)

為了準(zhǔn)確地預(yù)測電站鍋爐的NOx排放量，以某300 MW亞臨界循環(huán)流化床鍋爐為研究對象，利用自適應(yīng)風(fēng)驅(qū)動優(yōu)化(AWDO)算法和極端學(xué)習(xí)機(ELM)進行綜合建模，并根據(jù)不同工況下現(xiàn)場收集的樣本數(shù)據(jù)檢驗該模型的預(yù)測能力；將該模型的預(yù)測值與基本極端學(xué)習(xí)機、差分進化算法、粒子群算法和基本風(fēng)驅(qū)動算法優(yōu)化的極端學(xué)習(xí)機模型預(yù)測值進行比較.結(jié)果表明：AWDO算法可以更好地找到優(yōu)化參數(shù)，該算法優(yōu)化的極端學(xué)習(xí)機模型具有良好的預(yù)測精度和泛化能力，可以準(zhǔn)確、有效地預(yù)測電站鍋爐的NOx排放量.

NOx排放； 循環(huán)流化床鍋爐； 極端學(xué)習(xí)機； 風(fēng)驅(qū)動優(yōu)化算法

燃料燃燒產(chǎn)生的氮氧化物(主要是NO和NO2，總稱NOx)是大氣污染的主要有害物質(zhì)之一.煤炭的燃燒是NOx的主要來源之一[1].在中國，50%以上的NOx來源于電站鍋爐[2].NOx排放的控制過程分為2個階段，首先要建立NOx排放的穩(wěn)態(tài)模型，然后再優(yōu)化鍋爐運行參數(shù)[3].因此，建立有效的NOx排放特性預(yù)測模型對電站鍋爐的燃燒優(yōu)化至關(guān)重要.

然而，影響鍋爐NOx排放量的因素有很多，且鍋爐的燃燒過程具有強耦合、非線性等特征，很難用機理模型去描述.近年來，神經(jīng)網(wǎng)絡(luò)廣泛應(yīng)用于解決工程領(lǐng)域的建模與控制問題，但傳統(tǒng)神經(jīng)網(wǎng)絡(luò)訓(xùn)練算法需經(jīng)多次迭代調(diào)整才能確定網(wǎng)絡(luò)權(quán)值[4]，存在計算時間長、易陷入“過擬合”的缺點[5].對此，筆者提出了基于改進的自適應(yīng)風(fēng)驅(qū)動優(yōu)化(Adaptive Wind Driven Optimization，AWDO)算法和極端學(xué)習(xí)機(Extreme Learning Machine，ELM)[6]的綜合建模.WDO算法是美國Bayraktar等[7]在2010年提出的一種新的尋優(yōu)算法，其核心是應(yīng)用牛頓第二定律并結(jié)合理想氣體狀態(tài)方程研究空氣質(zhì)點在大氣中的受力、運動情況[7-9].WDO算法在實際工程領(lǐng)域如電磁學(xué)[7,9-10]、衛(wèi)星圖像處理[11]、云計算[12]等的成功應(yīng)用，證明其是一種有效的優(yōu)化算法.ELM通過求解線性方程組，解析求得具有最小范數(shù)的最小二乘解作為網(wǎng)絡(luò)權(quán)值，整個訓(xùn)練過程可一次完成無需迭代，這一優(yōu)勢使得ELM在模式識別與回歸估計等問題中獲得成功的應(yīng)用[13-15].

針對ELM模型的性能依賴隨機初始的權(quán)值和閾值的不足，筆者應(yīng)用AWDO算法對ELM模型的權(quán)值和閾值進行優(yōu)化，建立準(zhǔn)確、有效的NOx排放預(yù)測模型，并與其他幾種ELM模型預(yù)測結(jié)果進行了對比研究.

1　極端學(xué)習(xí)機

ELM模型首先隨機生成輸入權(quán)值和隱層閾值，然后該模型就變成一個線性系統(tǒng)，輸出權(quán)值則通過最小二乘計算得到，其簡單表述[16]如下.

ELM模型的每個觀測樣本輸出tk可表示為

(1)

上述N個方程可以簡寫為矩陣形式：

(2)

(3)

(4)

(5)

式中：H為隱層輸出矩陣；W=[w1,…,wm]T，為輸入權(quán)值；B=[b1,…,bm]T，為隱層閾值.

輸出權(quán)值β可通過MP廣義逆分析計算得到，由最小范數(shù)二乘解的求法可得：

(6)

式中：H+為H的MP廣義逆.

2　風(fēng)驅(qū)動算法

WDO算法是對抽象的空氣質(zhì)點在大氣中運動軌跡的模擬，簡單通用，容易實現(xiàn)，具有較好的全局收斂能力和魯棒性，可用于處理多維和多模問題，并對連續(xù)和離散優(yōu)化問題也適用[9].

2.1WDO算法的基本原理

在WDO算法中，以抽象的單位體積的空氣質(zhì)點為研究對象，并假設(shè)其處于流體靜力學(xué)平衡狀態(tài)，且滿足理想氣體狀態(tài)方程.由大氣動力學(xué)可知，大氣運動中的水平運動強于垂直運動，因此WDO算法只研究風(fēng)的水平運動.

對簡化的空氣質(zhì)點運用牛頓第二定律，并由理想氣體狀態(tài)方程可得：

(7)

(8)

式中：ρ為空氣質(zhì)點的密度；a為空氣質(zhì)點的加速度；Fi為作用在空氣質(zhì)點上的力；p為壓強；R為理想氣體常數(shù)；T為溫度.

WDO算法沒有完全模擬大氣運動，如將紊流阻力歸為摩擦力的一部分、忽略了平流層離心力等微小影響的作用力等，其主要考慮4個作用力：(1) 氣壓梯度力FPG，方向由高壓區(qū)指向低壓區(qū)，這是啟動空氣質(zhì)點運動的基本作用力；(2) 摩擦力FF，方向與FPG相反，這里用的是簡化的摩擦力表達(dá)式；(3) 重力FG，在物理學(xué)三維坐標(biāo)系中，其方向垂直指向地心，假如將地心考慮為直角坐標(biāo)系的原點，則可簡單地認(rèn)為重力是使質(zhì)點移向坐標(biāo)系原點的力，同樣地，將問題映射到N維中，重力表示指向坐標(biāo)系原點的力；(4) 科氏力FC，由地球旋轉(zhuǎn)產(chǎn)生，在WDO算法迭代中代表其他任一維度對空氣質(zhì)點當(dāng)前所在維度速度和位置的影響.各簡化表達(dá)式如下：

(9)

(10)

(11)

(12)

(13)

(14)

式中：p′為最優(yōu)壓力值；x′為最優(yōu)位置.

將式(8)、式(14)帶入式(13)，速度改變量Δu用U-u表示，U為空氣質(zhì)點下一次迭代中的速度矢量，并用綜合系數(shù)c來代替其他的系數(shù)，即c=-2ΩRT，可得：

(15)

式(15)中引入了壓力值p和p′，這可能會使得更新后的速度變得很大，從而降低WDO算法的可執(zhí)行性.因此按壓力值(或適應(yīng)度值)升序(或降序)排列后的序號i來代替壓力值，而在最優(yōu)位置x′時壓力值(或適應(yīng)度值)最小(或最大)，設(shè)為1.并令γ=1-α，且有γ∈(0,1)，為慣性權(quán)重，則速度公式更新為

(16)

則位置公式更新為

X=x+U

(17)

式中：X為空氣質(zhì)點下一次迭代的位置.

2.2AWDO算法

Shi等[17]通過實驗研究表明，較大的γ有助于跳出局部最優(yōu)，進行全局搜索；較小的γ有助于局部搜索，加速算法收斂，獲得更好的精度.為了克服WDO算法固定γ導(dǎo)致的算法易早熟、陷入局部最優(yōu)和收斂精度不高的缺點，筆者引入上述思想，對WDO算法的慣性權(quán)重γ進行改進，提出了隨機擾動的自適應(yīng)風(fēng)驅(qū)動優(yōu)化算法.

改進后的慣性權(quán)重表達(dá)式為

(18)

式(18)中的前2項使γ隨著迭代次數(shù)的增加而線性微分遞減，這使得迭代前期利于全局搜索，迭代后期利于局部尋優(yōu)；第3項是隨機擾動因子，其隨著迭代次數(shù)的增加而呈現(xiàn)非線性變化，保持了種群的多樣性，增加了跳出局部最優(yōu)的可能性.

AWDO算法流程如下：(1) 初始化.設(shè)置種群規(guī)模、最大迭代次數(shù)、相關(guān)參數(shù)(g，c，RT等)、搜索邊界和壓力函數(shù)(或適應(yīng)函數(shù))等；(2) 隨機設(shè)置空氣質(zhì)點的起始位置和速度；(3) 計算當(dāng)前空氣質(zhì)點的壓力值(或適應(yīng)度值)，并按照壓力值(或適應(yīng)度值)排序；(4) 按照式(16)更新空氣質(zhì)點的速度；(5) 按照式(17)更新空氣質(zhì)點的位置；(6) 若達(dá)到終止條件，則輸出最優(yōu)值，否則返回步驟(3).

2.3算法性能測試

為了驗證AWDO算法的性能，將其應(yīng)用到6個經(jīng)典的基準(zhǔn)優(yōu)化問題上(見表1)，其中函數(shù)f1、f2為高維的單峰基準(zhǔn)函數(shù)，函數(shù)f3、f4為高維的多峰基準(zhǔn)函數(shù)，函數(shù)f5、f6為固定維數(shù)的多峰基準(zhǔn)函數(shù).并將運行結(jié)果與差分進化(DE)算法、粒子群(PSO)算法和基本W(wǎng)DO算法進行比較.DE算法參數(shù)設(shè)置為比例系數(shù)F0=0.6，交叉率Cr=0.85.PSO算法參數(shù)設(shè)置為學(xué)習(xí)因子c1=c2=2，慣性權(quán)重γ′=1.4，γ″=0.4.基本W(wǎng)DO算法參數(shù)設(shè)置為重力加速度g=0.6，綜合系數(shù)c=0.7，其中RT=1，慣性權(quán)重γ=0.4.AWDO算法參數(shù)設(shè)置為慣性權(quán)重γ′=0.9，γ″=0.4，其他參數(shù)與基本W(wǎng)DO算法一樣.

為了比較的公平性，上述4種優(yōu)化算法的相關(guān)參數(shù)設(shè)置相同：最大迭代次數(shù)為1 000、種群規(guī)模為40、每個仿真實驗均獨立運行20次，每次均以不同的隨機初始種群開始.此外，對高維單峰和高維多峰函數(shù)的維數(shù)均依次設(shè)為10維、20維和50維3種.以表格的形式記錄各優(yōu)化算法在設(shè)定參數(shù)條件下達(dá)到最優(yōu)值的平均值Tm、均方差Tsd和運行一次的平均時間trun，如表2所示.

表1　6個經(jīng)典基準(zhǔn)測試函數(shù)

注：1) 其中n、2表示維數(shù).

從表2可以看出，在設(shè)定參數(shù)條件下，AWDO算法能夠找到f1、f3～f6的理論最優(yōu)值(或者非常接近，可以近似看成理論最優(yōu)值)，其中不包括f3在維數(shù)為10的時候.另外對f2來說，雖然AWDO算法精度沒有達(dá)到理論最優(yōu)，但與其他3個算法相比，其尋優(yōu)質(zhì)量提高了4個以上的數(shù)量級.雖然AWDO算法對f3的部分維數(shù)測試結(jié)果不如DE算法，但AWDO算法的trun卻遠(yuǎn)小于DE算法.因此，對于絕大部分測試函數(shù)來說，與其他3個算法相比，AWDO算法不僅找到最優(yōu)解的精度更高，而且收斂速度也比較快.

表2　4種算法對基準(zhǔn)測試函數(shù)的運行結(jié)果

3　AWDO-ELM模型設(shè)計及仿真

3.1數(shù)據(jù)樣本

以某熱電廠300 MW亞臨界循環(huán)流化床鍋爐為研究對象，其中300組測試數(shù)據(jù)是由集散控制系統(tǒng)(DCS)數(shù)據(jù)庫中每隔30 s采樣一次而得，其中包括了負(fù)荷在50%、70%和100%左右的數(shù)據(jù)各100組.選取了影響鍋爐NOx排放特性的參數(shù)(如鍋爐負(fù)荷、給煤量、床溫、氧量、一次風(fēng)量和二次風(fēng)量等22個參數(shù))作為輸入變量，以NOx排放質(zhì)量濃度作為輸出，應(yīng)用ELM建立模型.隨機選取270組樣本(每個負(fù)荷下各90組)作為訓(xùn)練樣本，剩下的30組樣本作為測試樣本，用于對訓(xùn)練好的模型進行預(yù)測.為使每個輸入變量對ELM的影響作用相同，在樣本數(shù)據(jù)輸入到ELM之前，對其進行了歸一化處理，處理后的輸入和輸出數(shù)據(jù)范圍均為[-1,1].

3.2仿真結(jié)果及分析

從圖1可以看出，除個別工況外，AWDO-ELM模型能夠?qū)τ?xùn)練樣本進行很好的擬合.從圖2可以看出，AWDO-ELM模型能夠?qū)y試樣本進行很好的預(yù)測.綜上所述，AWDO-ELM模型的辨識和預(yù)測能力很強，能夠精確地對循環(huán)流化床鍋爐NOx排放質(zhì)量濃度進行預(yù)測.

為了體現(xiàn)AWDO-ELM模型的預(yù)測效果，同時采用ELM模型、DE-ELM模型、PSO-ELM模型和WDO-ELM模型對循環(huán)流化床鍋爐NOx排放量進行預(yù)測.為了綜合評價各模型預(yù)測值與實際值的準(zhǔn)確度，引入以下3個性能指標(biāo)：均方差(MSE)、平均絕對誤差(MAE)和平均相對誤差率(MAPE).各模型輸入變量、輸出變量和數(shù)據(jù)樣本均與AWDO-ELM模型相同，各算法參數(shù)的設(shè)置與前文一致，預(yù)測結(jié)果如表3和表4所示.

圖1　訓(xùn)練樣本預(yù)測值與實際值的對比Fig.1　Comparison of training data between predicted results and actual measurements

圖2　測試樣本預(yù)測值與實際值的對比Fig.2　Comparison of test data between predicted results and actual measurements

從表3可以看出，對于訓(xùn)練樣本，AWDO-ELM模型的MSE為22.123 3，MAE為3.577 1，MAPE為0.027 0，3項指標(biāo)均小于其他4種模型.因此，對于訓(xùn)練樣本而言，AWDO-ELM模型的預(yù)測能力和精度均優(yōu)于其他4種模型.

表3　訓(xùn)練樣本的準(zhǔn)確度對比

從表4可以看出，對于測試樣本，AWDO-ELM模型的MSE為21.215 7，MAE為3.739 4，MAPE為0.028 6，3項指標(biāo)均遠(yuǎn)小于其他4種模型.另外，AWDO-ELM模型的最大MAE為9.390 0，最大MAPE為0.115 6，也均遠(yuǎn)小于其他4種模型.因此，對于測試樣本而言，AWDO-ELM模型的泛化能力和模型預(yù)測精度都遠(yuǎn)優(yōu)于其他4種模型.

表4　測試樣本的準(zhǔn)確度對比

結(jié)合表3和表4，對比各項指標(biāo)可知，雖然其他4種模型也能夠?qū)τ?xùn)練樣本進行很好的預(yù)測，但預(yù)測效果都沒有AWDO-ELM模型好.

各模型對測試樣本的NOx排放量預(yù)測相對誤差見圖3.從圖3可以看出，與其他4種模型相比，AWDO-ELM模型的相對誤差曲線比較平穩(wěn)，波動較小，最大相對誤差也比其他模型小，這說明AWDO-ELM模型對NOx排放量的預(yù)測結(jié)果更準(zhǔn)確.

綜上所述，AWDO-ELM模型的泛化能力更強，預(yù)測精度更高，非常適合循環(huán)流化床鍋爐NOx排放量的工程預(yù)測.

圖3　測試樣本的相對誤差對比Fig.3　Comparison of relative error of test data among different models

4　結(jié)　論

以某熱電廠300 MW亞臨界循環(huán)流化床鍋爐為研究對象，以現(xiàn)場采集的NOx燃燒特性數(shù)據(jù)為樣本，利用改進的自適應(yīng)風(fēng)驅(qū)動優(yōu)化算法優(yōu)化的極端學(xué)習(xí)機AWDO-ELM模型建立了循環(huán)流化床鍋爐NOx排放量預(yù)測模型，并將該模型的預(yù)測值與ELM模型、DE-ELM模型、PSO-ELM模型和WDO-ELM模型的預(yù)測值進行比較.結(jié)果表明，AWDO-ELM模型可以更準(zhǔn)確、有效地預(yù)測NOx排放量，為電站鍋爐預(yù)測NOx排放提供了一種新的方法.

[1]歐陽子區(qū), 朱建國, 呂清剛. 無煙煤粉經(jīng)循環(huán)流化床預(yù)熱后燃燒特性及NOx排放特性實驗研究[J]. 中國電機工程學(xué)報, 2014, 34(11): 1748-1754.

OUYANG Ziqu, ZHU Jianguo, Lü Qinggang. Experimental study on combustion and NOxemission of pulverized anthracite coal preheated by a circulating fluidized bed[J]. Proceedings of the CSEE, 2014, 34(11): 1748-1754.

[2]肖海平, 張千, 王磊, 等. 燃燒調(diào)整對NOx排放及鍋爐效率的影響[J]. 中國電機工程學(xué)報, 2011, 31(8): 1-6.

XIAO Haiping, ZHANG Qian, WANG Lei,etal. Effect of combustion adjustment on NOxemission and boiler efficiency[J]. Proceedings of the CSEE, 2011, 31(8): 1-6.

[3]牛培峰, 麻紅波, 李國強, 等. 基于支持向量機和果蠅優(yōu)化算法的循環(huán)流化床鍋爐NOx排放特性研究[J]. 動力工程學(xué)報, 2013, 33(4): 267-271.

NIU Peifeng, MA Hongbo, LI Guoqiang,etal. Study on NOxemission from CFB boilers based on support vector machine and fruit fly optimization algorithm[J]. Journal of Chinese Society of Power Engineering, 2013, 33(4): 267-271.

[4]張弦, 王宏力. 限定記憶極端學(xué)習(xí)機及其應(yīng)用[J]. 控制與決策, 2012, 27(8): 1206-1210.

ZHANG Xian, WANG Hongli. Fixed-memory extreme learning machine and its applications[J]. Control and Decision, 2012, 27(8): 1206-1210.

[5]魏輝, 陸方, 羅永浩, 等. 燃煤鍋爐高效、低NOx運行策略的研究[J]. 動力工程, 2008, 28(3): 361-366.

WEI Hui, LU Fang, LUO Yonghao,etal. Research on operation strategy on reducing NOxemissions and improving efficiency of coal-fired boiler[J]. Journal of Power Engineering, 2008, 28(3): 361-366.

[6]HUANG Guangbin, ZHU Qinyu, SIEW C K. Extreme learning machine: theory and applications[J]. Neurocomputing, 2006, 70(1/2/3): 489-501.

[7]BAYRAKTAR Z, KOMURCU M, WERNER D H. Wind Driven Optimization (WDO): a novel nature-inspired optimization algorithm and its application to electromagnetics[C]//2010 IEEE Antennas and Propagation Society International Symposium. Toronto, Canada: IEEE, 2010: 1-4.

[8]任作琳, 張儒劍, 田雨波. 風(fēng)驅(qū)動優(yōu)化算法[J]. 江蘇科技大學(xué)學(xué)報(自然科學(xué)版), 2015, 29(2): 153-158.

REN Zuolin, ZHANG Rujian, TIAN Yubo. Wind driven optimization algorithm[J]. Journal of Jiangsu University of Science and Technology (Natural Science Edition), 2015, 29(2): 153-158.

[9]BAYRAKTAR Z, KOMURCU M, BOSSARD J,etal. The wind driven optimization technique and its application in electromagnetics[J]. IEEE Transactions on Antennas and Propagation, 2013, 61(5): 2745-2757.

[10]BAYRAKTAR Z, KOMURCU M, JIANG Zhihao,etal. Stub-loaded inverted-F antenna synthesis via wind driven optimization[C]//Proceedings of IEEE International Symposium on Antennas and Propagation (APSURSI). Spokane, USA: IEEE, 2011: 2920-2923.

[11]BHANDARI A K, SINGH V K, KUMAR A,etal. Cuckoo search algorithm and wind driven optimization based study of satellite image segmentation for multilevel thresholding using Kapur's entropy[J]. Expert Systems with Applications, 2014, 41(7): 3538-3560.

[12]SUN Jiajia, WANG Xingwei, HUANG Min,etal. A cloud resource allocation scheme based on microeconomics and wind driven optimization[C]//Proceedings of the 8th China Grid Annual Conference. Changchun, China: IEEE, 2013: 34-39.

[13]MICHE Y, SORJAMAA A, BAS P,etal. OP-ELM: optimally pruned extreme learning machine[J]. IEEE Transactions on Neural Networks, 2010, 21(1): 158-162.

[14]LAN Yuan, SOH Y C, HUANG Guangbin. Constructive hidden nodes selection of extreme learning machine for regression[J]. Neurocomputing, 2010, 73(16/17/18): 3191-3199.

[15]MALATHI V, MARIMUTHU N S, BASKAR S. Intelligent approaches using support vector machine and extreme learning machine for transmission line protection[J]. Neurocomputing, 2010, 73(10/11/12): 2160-2167.

[16]李國強. 新型人工智能技術(shù)研究及其在鍋爐燃燒優(yōu)化中的應(yīng)用[D]. 秦皇島: 燕山大學(xué), 2013.

[17]SHI Y, EBERHART R. A modified particle swarm optimizer[C]//Proceedings of the 1998 IEEE International Conference on Evolutionary Computation Proceedings. Anchorage, USA: IEEE, 1998: 69-73.

Model Improvement for Boiler NOxEmission Based on Wind Driven Optimization Algorithm

NIUPeifeng,ZHAOZhen,MAYunpeng,CHENKe,WANGQiuya,ZHAOQingchong

(Key Lab of Industrial Computer Control Engineering of Hebei Province, Yanshan University,Qinhuangdao 066004, Hebei Province, China)

To accurately predict the NOxemission of power plant boilers, a unified model was established using adaptive wind driven optimization (AWDO) algorithm and extreme learning machine (ELM) by taking the 300 MW subcritical circulating fluidized bed boiler as an object of study. Prediction ability of the model was then verified with sample data obtained under different working conditions, and its prediction results were subsequently compared with that of basic ELM models and the ELM models optimized by differential evolution algorithm, particle swarm optimization and wind driven optimization, respectively. Results show that the AWDO algorithm has a strong capability in parameter optimization, and the ELM model optimized by AWDO algorithm has a higher prediction accuracy and generalization ability, which therefore is able to predict the NOxemission of power plants accurately and effectively.

NOxemission; circulating fluidized bed boiler; extreme learning machine; wind driven optimization algorithm

2015-10-26

國家自然科學(xué)基金資助項目(61573306，61403331)

牛培峰(1958-)，男，吉林舒蘭人，教授，博士生導(dǎo)師，研究方向為復(fù)雜工業(yè)系統(tǒng)的智能建模與智能控制和流程工業(yè)綜合自動化.電話(Tel.):0335-8072979；E-mail：npf882000@163.com.

1674-7607(2016)09-0732-07

TK224.9

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