牛培峰, 李進柏, 劉 楠, 李國強, 王榮彥
(1.燕山大學(xué) 工業(yè)計算機控制工程河北省重點實驗室,河北秦皇島 066004; 2.燕山大學(xué) 材料科學(xué)與工程學(xué)院,河北秦皇島 066004)
電廠鍋爐是大氣污染物NOx的主要排放源之一[1]。NOx易形成霧霾和酸雨,不僅污染環(huán)境而且會造成一定的經(jīng)濟損失。因此,如何降低NOx排放是電廠面臨的重要課題。目前,電廠控制NOx排放的方法主要有尾部煙氣脫硝和燃燒優(yōu)化技術(shù)[2]。煙氣脫硝技術(shù)雖然在一定程度上能降低NOx排放質(zhì)量濃度,但是需要安裝額外設(shè)備和后期維護,操作成本較高。燃燒優(yōu)化技術(shù)運行成本低,通過優(yōu)化鍋爐燃燒過程,使NOx排放質(zhì)量濃度達到期望目標(biāo)。
鍋爐燃燒優(yōu)化技術(shù)的基礎(chǔ)是建立一個精確的燃燒特性模型,并以此模型為前提,利用智能算法優(yōu)化鍋爐運行時的可調(diào)參數(shù),給出操作量最佳值,實現(xiàn)鍋爐清潔高效運行。很多學(xué)者致力于研究電站鍋爐燃燒優(yōu)化技術(shù),并取得了不錯的成果。李慶偉等[3]利用集成支持向量機建立600 MW四角切圓煤粉爐NOx排放特性模型,并采用改進的粒子群算法優(yōu)化配風(fēng)方式以降低NOx排放質(zhì)量濃度。張文廣等[4]采用自適應(yīng)模糊推理辨識方法建立循環(huán)流化床鍋爐NOx排放特性模型,然后利用果蠅優(yōu)化算法對鍋爐運行工況進行優(yōu)化,降低了NOx排放質(zhì)量濃度。余廷芳等[5]利用BP神經(jīng)網(wǎng)絡(luò)模型建立燃煤鍋爐NOx排放和熱效率特性模型,然后利用遺傳算法優(yōu)化燃燒工況,降低了NOx排放質(zhì)量濃度并提高了熱效率。
極限學(xué)習(xí)機(ELM)[6]是一種新型單層前饋神經(jīng)網(wǎng)絡(luò),只要確定輸入權(quán)值和隱層閾值,就可以通過最小二乘法確定網(wǎng)絡(luò)參數(shù),避免了因反復(fù)迭代陷入局部最小的問題,具有很強的擬合能力和泛化能力,能用于燃煤鍋爐NOx排放質(zhì)量濃度預(yù)測建模。但是ELM隨機生成的輸入權(quán)值和隱層閾值會影響模型性能,可以將權(quán)值和閾值的選取作為一個優(yōu)化問題,借助智能優(yōu)化算法進行尋優(yōu)。花授粉算法(FPA)[7]是一種新型啟發(fā)式算法,因結(jié)構(gòu)簡單、可調(diào)參數(shù)少、魯棒性好被應(yīng)用于無線電定位[8]、特征選擇[9]和光伏模塊參數(shù)估計[10]等工程領(lǐng)域。但算法中固定轉(zhuǎn)換概率和局部隨機搜索也會導(dǎo)致算法全局搜索與局部搜索能力失衡,搜索精度低。
筆者首先針對FPA收斂精度、搜索過程中全局搜索與局部搜索能力失衡等不足,在局部搜索中引入一種基于適應(yīng)度值的步長,在全局搜索和局部搜索轉(zhuǎn)化過程中構(gòu)造動態(tài)自適應(yīng)轉(zhuǎn)換概率,提出改進的花授粉算法(AFPA)。然后以某電廠330 MW燃煤鍋爐運行數(shù)據(jù)為基礎(chǔ),采用ELM建立 NOx排放特性模型,由于隨機權(quán)值和閾值的影響,利用AFPA對ELM的權(quán)值和閾值進行優(yōu)化,建立更加精確的NOx排放特性模型。在該模型的基礎(chǔ)上,利用AFPA對鍋爐運行時可調(diào)參數(shù)進行調(diào)整,以降低燃煤鍋爐NOx排放質(zhì)量濃度。
ELM的核心就是將單隱層神經(jīng)網(wǎng)絡(luò)的訓(xùn)練學(xué)習(xí)過程轉(zhuǎn)化為求解線性最小二乘問題,然后通過MP廣義逆計算出輸出權(quán)值。簡單表述如下:
每一個觀測樣本的輸出為:
(1)
上述方程可簡寫為:
Hβ=T
(2)
其中,
(3)
(4)
(5)
式中:H為隱層輸出矩陣;w1,…,wm為隨機輸入權(quán)值;b1,…,bm為隱層閾值。
輸出權(quán)值β可以通過下式計算得到:
(6)
式中:H+為H的MP廣義逆。
ELM的學(xué)習(xí)過程如下:(1)確定隱層節(jié)點數(shù),并選擇隱層激勵函數(shù);(2)隨機初始化輸入權(quán)值W和隱層閾值B;(3)按照式(3)計算出隱層輸出矩陣H;(4)按照式(6)計算出輸出權(quán)值β。
FPA受啟發(fā)于自然界中植物花朵傳粉過程,其搜索過程分為全局搜索和局部搜索,二者轉(zhuǎn)換概率由常數(shù)P控制。
全局搜索的位置更新變化公式為:
(7)
由于昆蟲能以各種步長移動一大段距離,Levy飛行可以有效模擬這種特性。L(λ)>0,并且
(8)
式中:Γ(λ)為標(biāo)準伽馬函數(shù);s和s0均為步長。
局部搜索的位置更新變化公式為:
(9)
FPA雖然具有良好的尋優(yōu)能力,但也存在尋優(yōu)精度不足和易陷入局部最優(yōu)等缺點。針對FPA的不足,進行了如下2處改進。
(1)FPA中P為固定常數(shù),易導(dǎo)致在迭代過程中算法的開采能力與探索能力失衡。為了平衡FPA全局搜索與局部搜索的能力,提出自適應(yīng)轉(zhuǎn)換概率:
(10)
式中:N為最大迭代次數(shù),F(xiàn)PA中固定概率取0.8時,算法表現(xiàn)良好,故常數(shù)仍取0.8。P值隨著算法迭代次數(shù)的增加而增大,使得算法前期傾向于全局搜索,而后期傾向于局部搜索,有效地平衡了算法的開采能力與探索能力。
(2)FPA中局部搜索的實質(zhì)是在解空間的隨機游走,易陷入局部最優(yōu)。為改善局部搜索能力,提出一種基于適應(yīng)度值的步長對局部搜索方案進行改進的方法。該步長表達式如下:
(11)
式中:ft和ft-1分別為當(dāng)前適應(yīng)度值和前一次迭代的適應(yīng)度值,且當(dāng)ft-1為0時,W=1。
于是,局部搜索的位置更新變化公式可改寫為:
(12)
W反映了前一次迭代對下一次迭代的修正。當(dāng)前一次迭代適應(yīng)度值改變較大時,提高前一次解在當(dāng)前迭代中所占比重;當(dāng)前一次適應(yīng)度值改變較小時,降低前一次解修正比重。W隨著適應(yīng)度值的變化不斷調(diào)整,避免算法在局部迭代中陷入局部最小。
鍋爐燃燒系統(tǒng)存在多個變量,各個變量之間相互耦合,存在復(fù)雜的非線性關(guān)系,很難用傳統(tǒng)的方法建立精確的鍋爐燃燒特性模型。ELM是一種新型的建模方法,無需了解鍋爐燃燒系統(tǒng)的運行狀況,只要根據(jù)鍋爐運行的歷史數(shù)據(jù)就可建立鍋爐燃燒特性的黑箱模型。ELM隨機給定的輸入權(quán)值和隱層閾值會影響模型的預(yù)測精度,權(quán)值和閾值的選擇可以看成是一個簡單的優(yōu)化問題,利用AFPA對輸入權(quán)值和隱層閾值進行尋優(yōu),以確保模型的性能。
在燃煤鍋爐燃燒過程中,影響NOx排放的因素很多,需要根據(jù)NOx排放質(zhì)量濃度與各個變量的關(guān)聯(lián)度來確定模型輸入?yún)?shù)。以某電廠330 MW八角雙切圓燃煤鍋爐為研究對象,該鍋爐配備4臺磨煤機、4臺排粉機和32臺給粉機,燃燒器采用三級點火方式。以集散系統(tǒng)(DCS)采集的20組數(shù)據(jù)為基礎(chǔ),選用鍋爐負荷(1個)、給煤機轉(zhuǎn)速(4個)、二次風(fēng)速(5個)、燃盡風(fēng)擋板開度(3個)、爐膛溫度(1個)、空氣預(yù)熱器出口煙氣平均含氧體積分數(shù)(以下簡稱含氧量,1個)、煤質(zhì)特性(7個)和一次風(fēng)速(4個)共26個對NOx排放影響關(guān)聯(lián)度較大的變量作為輸入,以NOx排放質(zhì)量濃度作為輸出。20組數(shù)據(jù)分為訓(xùn)練集和測試集2個部分,訓(xùn)練集有17個樣本(即樣本1~樣本17)和測試集有3個樣本(即樣本18~樣本20)。
在利用AFPA對ELM參數(shù)進行尋優(yōu)過程中,AFPA參數(shù)設(shè)置如下:種群數(shù)設(shè)置為20,迭代次數(shù)為200,由于輸入權(quán)值和隱層閾值的最小值和最大值均為-1和1,所以解的上下限設(shè)置為[-1,1]。ELM模型神經(jīng)元個數(shù)為30,隱層激勵函數(shù)為“sigmoid”,詳細表達式為:
(13)
由于數(shù)據(jù)樣本較小,為了最大程度地利用樣本信息和防止發(fā)生過擬合現(xiàn)象,優(yōu)化ELM參數(shù)時,對訓(xùn)練集的17個樣本進行“留一法”處理,即每次選擇16個樣本作為訓(xùn)練數(shù)據(jù),剩下一個樣本作為預(yù)測數(shù)據(jù),循環(huán)17次遍歷所有訓(xùn)練集樣本。以訓(xùn)練數(shù)據(jù)和預(yù)測數(shù)據(jù)的均方根誤差之和作為優(yōu)化目標(biāo)函數(shù),在200次迭代中得到最小優(yōu)化目標(biāo)函數(shù)的解,作為最優(yōu)權(quán)值和閾值參數(shù)帶入ELM。優(yōu)化目標(biāo)函數(shù)定義如下:
f=f1+f2
(14)
(15)
(16)
在確定最優(yōu)參數(shù)后,先將訓(xùn)練集17個樣本代入ELM,建立鍋爐NOx排放特性模型;然后將測試集3個樣本輸入所建模型,得到NOx排放質(zhì)量濃度的預(yù)測值,并與實際NOx排放質(zhì)量濃度進行對比,測試模型的預(yù)測精度。
圖1給出了訓(xùn)練集17個樣本NOx排放質(zhì)量濃度實際值與預(yù)測值的對比。由圖1可知,預(yù)測值能夠很好地跟蹤實際值,模型擬合效果較好。圖2給出了測試集3個樣本NOx排放質(zhì)量濃度實際值與預(yù)測值的對比。由圖2可知,對于未知樣本,所建模型能夠?qū)Ox排放質(zhì)量濃度進行預(yù)測,并且預(yù)測值與實際值相差不大,模型具有很好的預(yù)測能力和泛化能力。
為了進一步驗證模型性能,采用BP神經(jīng)網(wǎng)絡(luò)和基本ELM根據(jù)相同的數(shù)據(jù)建立鍋爐NOx排放特性模型,并進行對比。ELM模型參數(shù)與第3.2節(jié)相同,輸入權(quán)值和隱層閾值是隨機給定的。BP神經(jīng)網(wǎng)絡(luò)模型采用26-26-35-1的結(jié)構(gòu),各層傳遞函數(shù)均采用tansig函數(shù),表1給出了3個模型對測試集樣本NOx排放質(zhì)量濃度的預(yù)測結(jié)果。
圖1 訓(xùn)練集預(yù)測值與實際值的對比
Fig.1 Comparison of training set between predicted results and actual measurements
圖2 測試集預(yù)測值與實際值的對比
Fig.2 Comparison of test set between predicted results and actual measurements
表1 測試集樣本預(yù)測結(jié)果對比Tab.1 Comparison of test set predicted results
由表1可以看出,AFPA-ELM模型、ELM模型和BP神經(jīng)網(wǎng)絡(luò)模型都能對燃煤鍋爐NOx排放質(zhì)量濃度進行預(yù)測,就預(yù)測結(jié)果的精度而言,AFPA-ELM模型和ELM模型的預(yù)測精度優(yōu)于BP神經(jīng)網(wǎng)絡(luò)模型,AFPA-ELM模型的預(yù)測值相比于ELM模型更逼近實際值。AFPA-ELM模型對樣本18~樣本20的預(yù)測相對誤差分別為-2.80%、0.64%和2.02%,相對誤差比另外2種模型小。結(jié)合圖1、圖2和表1可以看出,AFPA-ELM模型對訓(xùn)練集樣本和測試集樣本的預(yù)測效果均良好,為實現(xiàn)鍋爐燃燒優(yōu)化提供了基礎(chǔ)。
以所建NOx排放特性模型為基礎(chǔ),選用DCS中可以直接調(diào)節(jié)的變量:給煤機轉(zhuǎn)速(4個)、二次風(fēng)速(5個)、燃盡風(fēng)擋板開度(3個)、含氧量(1個)和一次風(fēng)速(4個)共17個參數(shù)作為可調(diào)參數(shù),利用AFPA對這17個參數(shù)進行尋優(yōu),剩下的9個變量作為不可調(diào)節(jié)參數(shù)保持不變。將這17個參數(shù)定義為一個待優(yōu)化變量x=[x1,x2,…,x17],其中x1~x4為4個給煤機轉(zhuǎn)速,x5~x9為5個二次風(fēng)速,x10~x12為3個燃盡風(fēng)擋板開度,x13為含氧量,x14~x17為4個一次風(fēng)速。優(yōu)化目標(biāo)為降低NOx排放質(zhì)量濃度,優(yōu)化目標(biāo)函數(shù)描述為:
minfNOx(x)=f(x)
(17)
為了保證優(yōu)化結(jié)果的可靠性和可行性,根據(jù)現(xiàn)場采集數(shù)據(jù)的最大值和最小值及文獻[11]確定17個優(yōu)化參數(shù)的約束范圍,各參數(shù)約束范圍見表2。
表2 參數(shù)約束范圍Tab.2 Parameter constrains
隨機抽出樣本3和樣本13,采用AFPA對其可調(diào)參數(shù)進行優(yōu)化,以降低NOx排放質(zhì)量濃度。為了驗證AFPA的有效性,利用FPA對相同樣本進行優(yōu)化,2種算法的參數(shù)設(shè)置完全相同,經(jīng)過多次實驗,最大迭代次數(shù)設(shè)置為100,以100次尋到的最小值作為NOx排放質(zhì)量濃度的最優(yōu)值,根據(jù)對應(yīng)的解確定17個可調(diào)參數(shù)的值。2種算法對樣本3的迭代曲線見圖3和圖4,其優(yōu)化結(jié)果見表3。
圖3 AFPA對樣本3的優(yōu)化過程
Fig.3 Optimization process for sample 3 by AFPA
圖4 FPA對樣本3的優(yōu)化過程Fig.4 Optimization process for sample 3 by FPA
表3 2種算法的優(yōu)化結(jié)果Tab.3 Optimization results of two algorithms
由圖3和圖4可以看出,隨著迭代次數(shù)的增加,NOx排放質(zhì)量濃度逐漸下降,最后達到穩(wěn)定。但AFPA所需迭代次數(shù)比FPA少,說明AFPA比FPA尋優(yōu)速度快。由表3可知,相比于原始數(shù)據(jù),采用FPA和AFPA優(yōu)化后的NOx排放質(zhì)量濃度都有一定幅度降低,其中采用AFPA優(yōu)化后NOx排放質(zhì)量濃度的降幅更加明顯,這就說明相比于FPA,AFPA具有更好的尋優(yōu)性能。對比優(yōu)化結(jié)果可以發(fā)現(xiàn),降低給煤機轉(zhuǎn)速可以降低NOx排放質(zhì)量濃度,但是若給煤量波動較大會影響鍋爐負荷,在實際操作過程中,可以對給煤量進行微調(diào)。降低一次風(fēng)速和二次風(fēng)速使得NOx排放質(zhì)量濃度降低,一次風(fēng)量和二次風(fēng)量減少使得煤粉充分燃燒的可能性降低,可以降低NOx排放質(zhì)量濃度。燃盡風(fēng)擋板開度增大,增加燃盡區(qū)風(fēng)率使得燃燒器內(nèi)過量空氣系數(shù)降低,爐內(nèi)平均火焰溫度下降,減少了熱力型NOx的生成,這與文獻[12]現(xiàn)場實驗結(jié)果一致。降低含氧量可以在很大程度上抑制煤中揮發(fā)分N向NO的轉(zhuǎn)換過程,減少燃料型NOx的生成,這符合燃料型NOx產(chǎn)生機理。
為降低燃煤鍋爐NOx排放質(zhì)量濃度,首先針對花授粉算法的不足,提出改進的花授粉算法,并利用該算法優(yōu)化ELM參數(shù),建立330 MW燃煤鍋爐NOx排放特性模型,相比于ELM模型和BP神經(jīng)網(wǎng)絡(luò)模型,所建模型具有更高的預(yù)測精度,可以對鍋爐NOx排放質(zhì)量濃度進行預(yù)測,為預(yù)測電站鍋爐NOx排放質(zhì)量濃度提供了一種新的方法。在該模型的基礎(chǔ)上,分別利用花授粉算法和改進的花授粉算法對鍋爐可調(diào)參數(shù)進行優(yōu)化,2種算法優(yōu)化之后NOx排放質(zhì)量濃度都有所降低,但是改進花授粉算法優(yōu)化后的效果更加明顯,證明該算法是一種高效的算法。