楊 晉 芳

(1.煤科院節(jié)能技術(shù)有限公司，北京 100013； 2.煤炭資源高效開采與潔凈利用國家重點實驗室，北京 100013；3.國家能源煤炭高效利用與節(jié)能減排技術(shù)裝備重點實驗室，北京 100013)

0 引 言

煤粉工業(yè)鍋爐是近十年來發(fā)展起來的1種新型工業(yè)鍋爐，以其高效、節(jié)能、環(huán)保等優(yōu)點代替鏈條爐并迅速進入市場。煙氣氧含量是燃煤鍋爐的重要指標(biāo)之一[1]，其對鍋爐的效率和排放有著重要的影響[2-3]，很多學(xué)者致力于鍋爐氧含量優(yōu)化的相關(guān)研究[4-6]。

目前煤粉工業(yè)鍋爐煙氣氧含量的測量主要通過直接測量方式，借助氧化鋯傳感器或熱磁式傳感器測量煙氣氧含量[7]，以上儀器可直接顯示氧含量的百分比，其反應(yīng)速度快、測量精度較高，然而也存在以下3個問題:① 價格成本高；② 氧化鋯探頭在長期使用中會被粉塵堵塞而出現(xiàn)老化現(xiàn)象，嚴(yán)重影響儀器的穩(wěn)定性和準(zhǔn)確性；③ 氧量儀還需定期進行效驗。因此，此次研究考慮采用間接測量技術(shù)對煙氣氧含量進行檢測。

間接測量采用軟測量的方法[8],即利用易于測量的一些輔助變量以估計主要變量,已在電站鍋爐上推廣應(yīng)用并取得較好的效果。采用軟測量方法間接獲得煙氣氧含量，可與儀器測量值比較并對其進行修正，同時可實現(xiàn)對氧含量的實時監(jiān)控，從而提高鍋爐的燃燒效率[9]。

軟測量的建模方法主要分為以下4種[10]，分別基于工藝機理分析方法、模式識別方法、人工神經(jīng)網(wǎng)絡(luò)方法和支持向量機方法。基于工藝機理分析方法存在滯后，很難符合實際工業(yè)現(xiàn)場對監(jiān)測數(shù)據(jù)實時性的需求；基于模式識別方法通常與人工神經(jīng)網(wǎng)絡(luò)和支持向量機等方法相結(jié)合，并應(yīng)用至實際生產(chǎn)中；基于人工神經(jīng)網(wǎng)絡(luò)方法易形成局部極小而得不到全局最優(yōu)，很多學(xué)者針對此項缺點對其進行改進?；诨貧w支持向量機方法應(yīng)用于回歸估計問題中，得到很好的效果，此方法預(yù)測精度很高且泛化能力強。目前已有較多基于此方法進行軟測量模型的建立，并已應(yīng)用至工程實際。但該方法對于大規(guī)模訓(xùn)練樣本難以實施，建模精度不高。至目前為止，最常用的軟測量方法主要包括基于人工神經(jīng)網(wǎng)絡(luò)和支持向量機的方法。

針對基于人工神經(jīng)網(wǎng)絡(luò)和支持向量機方法的煙氣氧含量間接預(yù)測技術(shù)，已有大量研究并成功應(yīng)用于電站鍋爐。胡世廣[11]利用神經(jīng)網(wǎng)絡(luò)進行火電廠煙氣含氧量的測量并進行建模研究，為火電廠的氧量測量技術(shù)提供新手段。馬良玉[12]等基于L-M算法改進神經(jīng)網(wǎng)絡(luò)模型以更準(zhǔn)確地預(yù)測煙氣含氧量,為實現(xiàn)鍋爐的燃燒優(yōu)化調(diào)整、節(jié)煤降耗奠定基礎(chǔ)。Zhenghao Tang[13]等使用DBN深度學(xué)習(xí)模型對電廠鍋爐氧含量進行預(yù)測，結(jié)果顯示可準(zhǔn)確地預(yù)測結(jié)果。Ma Liangyu[14]等基于ANN對1000MW電廠煙氣氧含量進行預(yù)測，結(jié)果顯示此模型可精確預(yù)測煙氣氧含量，對電廠改善鍋爐燃燒、節(jié)約能源和減少煤耗具有指導(dǎo)作用。楊秀等[15]采用1種基于粒子群改進的支持向量機算法對氧含量軟測量進行建模研究，仿真結(jié)果表明預(yù)測精度較高、泛化能力良好。任錦等[16]提出改進最小二乘支持向量機模型PSO-LSSVM并對氧含量進行預(yù)測，最終仿真結(jié)果顯示:改進的軟測量模型預(yù)測精度更高且泛化性更好、模型可靠性更高。Changliang Liu等[17]使用LS-SVM對煙氣氧含量進行預(yù)測，通過預(yù)測值和實際值比較可知，其提出的模型可較準(zhǔn)確地預(yù)測電廠煙氣氧含量。

以上研究均針對電站鍋爐，而目前對煤粉工業(yè)鍋爐氧含量的預(yù)測仍很稀少，由此筆者提出煤粉工業(yè)鍋爐氧含量預(yù)測模型，以期為工業(yè)煤粉鍋爐燃燒系統(tǒng)優(yōu)化提供指導(dǎo)。

1 現(xiàn)場數(shù)據(jù)采集

現(xiàn)場數(shù)據(jù)采集基于煤粉工業(yè)鍋爐系統(tǒng)的以下工藝流程:煤粉由供料器均勻送入鍋爐燃燒器并在爐膛內(nèi)燃燒，煙氣經(jīng)爐膛內(nèi)的SNCR(Selective Non-Catalytic Reduction)選擇性非催化還原脫硝和尾部的SCR(Selective Catalytic Reduction)選擇性催化還原脫硝后，進入NGD(No Gap Desulfurization)半干法脫硫除塵一體化裝置進行脫硫除塵，潔凈煙氣排至煙囪進入大氣，鍋爐飛灰經(jīng)倉泵輸送至密閉灰塔排出，經(jīng)封閉罐車運輸出廠。

煤粉鍋爐DCS(Distributed Control System)分布式控制系統(tǒng)中的歷史數(shù)據(jù)是整個鍋爐工藝過程的所有設(shè)備測點數(shù)據(jù)，而氧含量的預(yù)測只與部分工藝相關(guān)，即進料、燃燒和煙氣排放3個過程，因此采集數(shù)據(jù)時只需篩選出關(guān)于3個過程的相關(guān)變量即可。

以某礦20 t/h煤粉工業(yè)鍋爐為研究對象，以20 min為采樣時間，采集了采暖高峰期2個月的數(shù)據(jù)。根據(jù)系統(tǒng)機理和工藝流程對變量進行初選，選擇了給料頻率、引風(fēng)機頻率、爐膛負(fù)壓、前爐膛溫度、省煤器出口煙溫等20個變量。以上變量分為控制變量和狀態(tài)變量兩大類。控制變量為技術(shù)操作可控制的變量，采集的控制變量為6個，分別為給料頻率、補水閥開度、引風(fēng)機頻率、二次風(fēng)閥開度、進風(fēng)閥1開度、進風(fēng)閥2開度。狀態(tài)變量是反應(yīng)鍋爐狀態(tài)的一些變量，不能直接控制，初篩包括狀態(tài)變量14個，分別為蒸汽流量、前爐膛溫度、后爐膛溫度、一次風(fēng)壓、省煤器入口煙溫、省煤器入口水溫、省煤器出口煙溫、省煤器出口水溫、蒸汽壓力、二次風(fēng)流量、爐膛負(fù)壓、蒸汽溫度、爐尾負(fù)壓、汽包水位。

2 數(shù)據(jù)預(yù)處理和變量篩選

2.1 數(shù)據(jù)預(yù)處理

從現(xiàn)場運行的歷史數(shù)據(jù)庫中所采集的數(shù)據(jù)一般情況下都含有異常數(shù)據(jù)及干擾數(shù)據(jù)，其會影響模型的訓(xùn)練，因此需對歷史數(shù)據(jù)進行預(yù)處理[18]。此次研究刪除異常值和缺失值后，最終得到2 598條數(shù)據(jù)。

2.2 變量篩選

輸入的變量直接影響模型的預(yù)測準(zhǔn)確性，太多的輸入變量將會降低預(yù)測精度并增加模型的計算時間，因此有必要將冗余數(shù)據(jù)去除。通過歷史數(shù)據(jù)集的相關(guān)性選擇，使用F檢驗判斷每個輸入變量和煙氣氧含量的關(guān)系。F檢驗、方差齊性檢驗是用來捕捉每個特征與標(biāo)簽之間線性關(guān)系的過濾方法。使用F檢驗將篩選出與氧含量有著顯著關(guān)系的變量。通過以下公式(1)和公式(2)計算F值，F(xiàn)值越大則輸入變量與氧含量之間的相關(guān)性就越大。

(1)

(2)

式中，Xi代表所有樣本在i號特征上取值的n維列向量，y是氧含量;ri是樣本的相關(guān)系數(shù)。

經(jīng)過F檢驗，從上述20個變量中選取與氧含量關(guān)系顯著的變量:蒸汽流量、給料頻率、后爐膛溫度、一次風(fēng)壓、二次風(fēng)閥開度、進風(fēng)閥1開度、進風(fēng)閥2開度、省煤器入口煙溫、引風(fēng)機頻率、蒸汽壓力、爐膛負(fù)壓、二次風(fēng)流量、爐尾負(fù)壓、蒸汽溫度14個變量，其中控制變量5個，狀態(tài)變量9個。

2.3 數(shù)據(jù)標(biāo)準(zhǔn)化

從DCS中采集的歷史數(shù)據(jù)，其各變量有著不同的單位，一些變量大小在數(shù)值上相差很大，有的甚至達(dá)到幾個數(shù)量級。為了不影響模型的精度，需要將采集到的數(shù)據(jù)進行Z-score 標(biāo)準(zhǔn)化，以消除因為數(shù)據(jù)數(shù)量級相差大造成的變量對模型作用差異大。

Z-score標(biāo)準(zhǔn)化公式詳見如下公式(3):

(3)

經(jīng)過Z-score標(biāo)準(zhǔn)化的所有數(shù)據(jù)都聚集在0附近，方差為1，經(jīng)過處理后的數(shù)據(jù)不會造成因數(shù)據(jù)的數(shù)量級差異而影響模型精度等問題。

3 建立模型

3.1 支持向量回歸(SVR)

支持向量機是用于分類的算法，支持向量也可用于回歸，稱之為支持向量回歸(SVR)。支持向量機可分為支持向量分類和支持向量回歸[19]，支持向量回歸因其優(yōu)越的學(xué)習(xí)性能，使其在系統(tǒng)辨識、預(yù)測估計等領(lǐng)域被廣泛研究應(yīng)用。文中預(yù)測氧含量值，采用支持向量回歸算法。

給定數(shù)據(jù)集:D={(x1,y1),(x2,y2),……,(xn,yn)},yn∈R，且令:

f(x)=wTx+b

(4)

w和b是待確定的模型參數(shù)，T為轉(zhuǎn)置符號，希望學(xué)得1個回歸模型,使得f(x)與y盡可能接近。假設(shè)能容忍f(x)與y之間最多有ε的偏差，如圖1所示。

圖1 支持向量回歸模型原理圖Fig.1 Schematic diagram of support vector regression model

即僅當(dāng)f(x)與y之間的差別絕對值大于ε時才計算損失。于是，SVR問題可形式化為:

(5)

式中,C為正則化常數(shù);lε是圖1中所示的ε-不敏感損失函數(shù):

(6)

引入松弛變量和拉格朗日乘子，最終可解得SVR為:

(7)

3.2 建立模型

上文經(jīng)過F檢驗已篩選出14個變量，其中控制變量5個、狀態(tài)變量9個。模型中的變量見表1?？紤]到輸入變量的性質(zhì)，將模型訓(xùn)練分為控制變量、狀態(tài)變量以及兩者合并3種方式，并分別對3個訓(xùn)練的模型依次進行測試。

表1 模型中的變量詳情Table 1 Variables details in model

3.3 模型訓(xùn)練與預(yù)測

將經(jīng)過處理的2 598條數(shù)據(jù)的70%作為訓(xùn)練集，建立控制變量、狀態(tài)變量以及兩者合并的支持向量回歸模型并進行訓(xùn)練，30%作為測試集測試模型的性能。3種情況下煙氣氧含量軟測量模型的預(yù)測結(jié)果分別如圖2、圖3和圖4所示，藍(lán)色實線代表直接測量所得的真實值，桔色虛線代表預(yù)測值。

從圖2～圖4可看出，3種軟測量模型的煙氣氧含量預(yù)測值均能較為準(zhǔn)確地反映真實值的變化趨勢。

將圖2與圖3和圖4分析比較可知，圖2預(yù)測值與真實值差距偏大，控制變量調(diào)節(jié)后可即時反饋到從DCS中取得的數(shù)據(jù)上，例如給料頻率從20調(diào)整為10，取得的歷史數(shù)據(jù)顯示為10，而實際鍋爐工況需要一定的時間才能基本穩(wěn)定地運行在給料頻率為10的工況。因此若取到此組數(shù)據(jù)，瞬時取到的氧含量值有一定的誤差。故若只基于控制變量來預(yù)測模型，結(jié)果和直接測量之間的誤差較大。而基于狀態(tài)變量的預(yù)測結(jié)果符合最佳，出現(xiàn)以上現(xiàn)象是由于狀態(tài)變量較控制變量有一定的延遲，和直接測量得到的氧含量較為同步。對于基于控制變量和狀態(tài)變量兩者合并的SVR模型，其準(zhǔn)確度介于兩者之間。

圖2 基于控制變量的預(yù)測值與真實值比較Fig.2 Comparison of prediction values and real values based on control variables

圖3 基于狀態(tài)變量的預(yù)測值與真實值比較Fig.3 Comparison of prediction values and real values based on state variables

圖4 控制變量和狀態(tài)變量合并的預(yù)測值與真實值比較Fig.4 Comparison of prediction values and real values based on both control and state variables

3.4 模型評價

為了評價預(yù)測模型，使用均方誤差(MSE)、平均絕對誤差(MAE)和決定系數(shù)(R-Square)來衡量預(yù)測模型的性能[20]，其計算見公式(8)～公式(10):

(8)

(9)

(10)

式中m——測試數(shù)據(jù)的樣本數(shù)量，個；

yi——煙氣氧含量真實值，%；

對測試結(jié)果分析詳見表2，MSE和MAE的值越小，說明模型精度越高。而R2取值范圍為[0,1]，R-Squared越大，表示模型擬合效果越好。評價指標(biāo)進一步說明基于狀態(tài)變量的SVR模型效果最佳。

表2 模型評價指標(biāo)Table 2 Model evaluation indicators

4 結(jié) 語

煤粉工業(yè)鍋爐的燃燒工況較為復(fù)雜，而直接測量氧含量儀器成本高，還受到鍋爐運行積灰等影響，因此軟測量越來越受到人們的重視，也逐漸應(yīng)用至實際工業(yè)生產(chǎn)中。煙氣氧含量預(yù)測技術(shù)已在電站鍋爐上得以推廣應(yīng)用，對在工業(yè)鍋爐上的應(yīng)用具有借鑒意義。

為了解決工業(yè)鍋爐的氧含量軟測量問題，筆者提出了支持向量回歸模型對工業(yè)鍋爐的煙氣進行氧含量測量。整個預(yù)測模型分為以下3個部分:

(1)選取煙氣氧含量的主要影響變量并將其對應(yīng)的數(shù)據(jù)標(biāo)準(zhǔn)化，使用了F檢驗篩選出對氧含量有顯著影響的變量，最終選取14個變量作為輸入，為了防止因不同的變量數(shù)值上數(shù)量級相差較大進而造成訓(xùn)練模型時各變量作用相差巨大，對輸入數(shù)據(jù)進行Z-score標(biāo)準(zhǔn)化。

(2)建立模型。將第一步選出來的變量分為控制變量和狀態(tài)變量，并運用SVR理論構(gòu)建分別基于控制變量、狀態(tài)變量和兩者合并的氧含量預(yù)測模型，結(jié)果表明3種模型的趨勢與測量值一致，基于狀態(tài)變量SVR模型優(yōu)于其它兩者，基于狀態(tài)變量的模擬結(jié)果能較好地預(yù)測實際值的趨勢，為工業(yè)煤粉鍋爐燃燒系統(tǒng)優(yōu)化提供指導(dǎo)。

(3)模型評價。為了評價預(yù)測模型，使用均方誤差、平均絕對誤差和決定系數(shù)來衡量預(yù)測模型的性能，進一步證明了基于狀態(tài)變量的SVR模型效果最佳。