常永亮

(京能秦皇島熱電有限公司，河北秦皇島066300)

0 引言

火電機組熱工一次風量測點在磨煤機系統(tǒng)運行過程中起到重要的調(diào)節(jié)作用，影響磨煤機一次風量的2個調(diào)門分別是熱風調(diào)節(jié)門和冷風調(diào)節(jié)門，熱風調(diào)節(jié)門主要調(diào)節(jié)磨煤機入口一次風量，冷風調(diào)節(jié)門主要調(diào)節(jié)磨煤機出口風溫。因此，一次風量測點的準確與否直接影響鍋爐的燃燒以及磨煤機的安全運行。據(jù)統(tǒng)計，河北某熱電機組熱工專業(yè)1個月內(nèi)發(fā)生6次磨煤機一次風量測點堵塞缺陷，根據(jù)此缺陷，熱工專業(yè)已將磨煤機一次風量測點反吹掃工作作為定期工作執(zhí)行?；诖?，磨煤機一次風量在線故障預警研究具有十分重要的研究意義。

目前，故障預警系統(tǒng)主要有基于機理模型的故障預警系統(tǒng)[1]和基于數(shù)據(jù)驅(qū)動模型的故障預警系統(tǒng)。發(fā)電機組由于工況變化復雜，被控對象受外界影響大，系統(tǒng)之間相互耦合導致機理建模的模型不能十分準確地匹配現(xiàn)場實際情況。隨著智能控制的不斷成熟，基于現(xiàn)場數(shù)據(jù)建模的故障預警系統(tǒng)逐漸得到應用，并在一定程度上優(yōu)于機理建模方法。

現(xiàn)有的火電機組分散控制系統(tǒng)DCS（distributed control system）功能十分完善，擁有詳細的故障預警功能，但是火電機組中各個測點之間都是相互耦合、相互關(guān)聯(lián)的，而DCS本身自帶的故障報警功能都是基于單點上下限的報警模式，缺乏多變量的數(shù)據(jù)融合報警功能?；谏舷孪薜膱缶Ｊ酵荒芗皶r發(fā)現(xiàn)故障，提前處理?；诖?，本文采用多變量相互融合的故障預警模型進行磨煤機一次風量故障預警研究，采用多變量狀態(tài)估計MSET（multivariate state estimation techniques）的方法對磨煤機一次風量進行建模，進行故障預警。

1 多變量狀態(tài)估計

多變量狀態(tài)估計算法最初是應用于檢測核電機組傳感器、設(shè)備以及機組運行參數(shù)劣化的非參數(shù)、非機理的建模方法。近年來，隨著風電機組裝機容量不斷擴大，參考文獻[2—5]的作者將該算法逐漸應用在風電機組齒輪箱的故障預警研究應用中。MSET算法采用觀測向量來組成該算法的估計矩陣，該估計矩陣盡量包含了機組全部的運行工況。在火電機組、風電機組、核電機組中，任何一輔機或者主機的監(jiān)測參數(shù)都是相互關(guān)聯(lián)、相互耦合的，基于此，多變量狀態(tài)估計算法將這些主要的監(jiān)測變量組成向量的形式來作為估計矩陣中的其中某一工況點?；谝陨显?，估計矩陣幾乎包含了某一設(shè)備、某一系統(tǒng)、某一機組的所有工況點。MSET采用相似度函數(shù)來檢測目前采集到的觀測向量是不是在估計矩陣中的某一工況點，當相似度很低時說明某一設(shè)備、某一系統(tǒng)、某個機組正處在故障邊緣，此時需要提醒運行或者熱控人員開始排查設(shè)備、系統(tǒng)、機組的潛在故障點。相比于DCS上下限的故障報警模式，該方法不但能夠?qū)崟r監(jiān)測設(shè)備、系統(tǒng)、機組的運行狀態(tài)，還能夠?qū)⒏鱾€監(jiān)測變量融合到一起提前發(fā)現(xiàn)設(shè)備的劣化趨勢。MSET的詳細數(shù)學推導公如下，其中估計矩陣記為D，估計矩陣中的任何一列代表狀態(tài)向量。某一時刻新觀測到的狀態(tài)向量代表觀測向量，記為Xobs，則

其中，n為觀測向量變量個數(shù)；i為采樣時刻。估計矩陣總的狀態(tài)向量個數(shù)記作m，所以多變量狀態(tài)估計矩陣維數(shù)為n×m，估計矩陣記為

多變量狀態(tài)估計模型將觀測向量Xobs與估計矩陣D中的狀態(tài)向量進行對比，然后計算出估計向量Xest，估計向量Xest為估計矩陣D與權(quán)值向量W的積，記為

其中，權(quán)值向量W通過最小化殘差方差ε來進行計算求出，ε＝Xest－Xobs，在保證殘差方差最小的情況下求出權(quán)值向量W，則

基于以上計算原理，可以得知觀測向量與估計向量之間的殘差就可以代表當前磨煤機一次風量系統(tǒng)是否存在故障，此時數(shù)值是否準確。對于兩個狀態(tài)向量之間的殘差大小，本文采用式（7）進行計算。

當數(shù)值S越大，則當前設(shè)備出現(xiàn)故障的概率越大。為了保證故障預警精度更好，更貼近現(xiàn)場實際情況，針對模型預警殘差進行滑動窗口平滑處理，保證預警更為有效，從一定程度上避免了該模型出現(xiàn)誤報現(xiàn)象。假設(shè)多變量狀態(tài)預警模型預測殘差點為n個，則計算n個殘差點的平均值代表此時故障預警殘差值，然后按照滑動窗口大小為n的窗口依次移動一個點，同理下一時刻的殘差點則為移動后的n個殘差的平均值，其原理如圖1所示。

圖1 滑動窗口原理圖

本文采用MATLAB軟件進行仿真實驗。估計矩陣的主要數(shù)據(jù)來源由機組正常運行過程中機組各個運行工況點的數(shù)據(jù)組成，該數(shù)據(jù)從DCS歷史數(shù)據(jù)中導出使用。

2 磨煤機一次風量測量裝置

本文研究對象為河北某電廠正壓直吹式磨煤機一次風量，該機組為超臨界35萬kW熱電聯(lián)產(chǎn)機組，單臺機組配置5臺正壓直吹式磨煤機。一次風量測量裝置采用差壓式原理進行測量，每臺磨煤機一次風量測點配置3個獨立的取樣裝置，經(jīng)過流量計算書得出一次風流量。其中，3個一次風流量取樣后作為磨煤機熱一次風調(diào)節(jié)門的被調(diào)節(jié)對象，其準確度直接影響風煤比、機組燃燒、磨煤機安全運行等。在該機組中經(jīng)常發(fā)生一次風量測量裝置取樣管積灰現(xiàn)象，致使風量波動較大，經(jīng)常導致磨煤機熱一次風調(diào)門由自動狀態(tài)切換至手動狀態(tài)。

3 預警模型建立

建立磨煤機一次風量故障預警模型，首先選取建模變量。在選取建模變量前通過機理分析需要的相關(guān)建模變量。為了能夠全面準確地表征磨煤機一次風量，建模變量的選取需要考慮從一次風機入口到磨煤機出口的整個工藝流程，選取典型的代表變量進行建模分析[7]?，F(xiàn)選取主要相關(guān)變量為：11號磨煤機一次風機電流IA、11號磨煤機一次風機變頻器頻率HA、12號磨煤機一次風機電流IB、12號磨煤機一次風機變頻器頻率HB、11號磨煤機電流IM、11號給煤機瞬時煤量F、11號磨煤機冷一次風調(diào)門開度FB1、11號磨煤機熱一次風調(diào)門開度FB2、11號磨煤機進出口一次風差壓DP1、進口一次風溫度T1，11號磨煤機一次風量CF。其中1號機組11、12號磨煤機一次風機為變頻調(diào)節(jié)，事故工況下切換至工頻方式，配合一次風機入口電動調(diào)節(jié)擋板調(diào)節(jié)一次風母管壓力。磨煤機一次風量故障預警模型能否準確地進行判斷，關(guān)鍵在于估計矩陣中數(shù)據(jù)的選取。本文仿真所用數(shù)據(jù)均來自于DCS歷史數(shù)據(jù)，所以不需要再次進行濾波等數(shù)據(jù)處理。由于不同變量之間的量綱均不同，為了保證模型預測更為準確，本文對所有數(shù)據(jù)進行歸一化處理，刪除量綱對估計模型計算的影響，歸一化區(qū)間為[0 1]，本文采用MATLAB中自帶函數(shù)Mapminmax進行數(shù)據(jù)歸一化處理[8]。為了保證估計矩陣中數(shù)據(jù)能夠覆蓋機組所有運行工況且數(shù)據(jù)量盡可能小的情況以保證計算迅速和節(jié)省電腦內(nèi)存使用率，參考文獻[9]中基于相似度的方法來動態(tài)更新估計矩陣數(shù)據(jù)，從而保證模型計算既準確又運算速度快?；谙嗨贫确椒ㄟM行數(shù)據(jù)篩選的基本思想是在假設(shè)機組所有運行工況正常的情況下，利用相似度函數(shù)計算2個狀態(tài)向量之間的相似程度。針對相似度很高的2組狀態(tài)向量只保留其中1個狀態(tài)向量，這樣既保證了工況點比較全面又保證了狀態(tài)向量之間沒有冗余情況。相似度公式為

其中，Rij代表2個狀態(tài)向量的相似度大小，xi代表第i個狀態(tài)向量，xj代表第j個狀態(tài)向量，δ代表相似度函數(shù)歸一化參數(shù)，||||則代表矩陣的2個范數(shù)。

基于相似度數(shù)據(jù)選取的原則為：Rij≤ε，保留兩組狀態(tài)向量，其中ε為兩組狀態(tài)向量的相似度；Rij＞ε，刪除其中一組狀態(tài)向量。參數(shù)δ的計算公式為

其中，Di為第i個參數(shù)的數(shù)值，m為變量個數(shù)。

4 故障預警模型仿真

本次仿真數(shù)據(jù)從DCS歷史站中導出某天24 h變量的歷史數(shù)據(jù)，該數(shù)據(jù)段包含了機組整個負荷變化范圍，機組負荷變化范圍在175～350 MW之間，幾乎涵蓋了該機組全部正常運行范圍，以這段數(shù)據(jù)進行仿真驗證。其中估計矩陣數(shù)據(jù)使用以上24 h數(shù)據(jù)建立模型，隨機截取其中50個向量作為觀測向量進行仿真驗證。為了檢驗該模型能夠準確識別觀測向量中的異類向量，本文采用在以上50個正常觀測向量中的11號磨煤機一次風量上逐漸添加偏移量來檢測該模型的有效性。首先檢驗該模型檢測正常的50個觀測向量的仿真，然后通過觀察估計向量與觀測向量的殘差檢驗該模型的有效性，其殘差曲線如圖2所示，擬合對比曲線如圖3所示。

圖2 檢驗正常觀測向量殘差曲線

圖3 檢驗正常觀測向量擬合對比曲線

在正常觀測向量中的一次風流量變量上逐步添加步距為0.05的誤差，再次對MSET故障預警模型檢測故障的有效性進行仿真實驗，然后根據(jù)仿真實驗和現(xiàn)場實際情況確定報警閥限值。

5 結(jié)論

通過以上仿真可以得知基于多變量狀態(tài)估計的磨煤機一次風量系統(tǒng)故障預警模型能夠及時判斷該系統(tǒng)當前的健康狀況，能夠根據(jù)趨勢及時判斷和發(fā)現(xiàn)該系統(tǒng)當前存在的安全隱患，相比于DCS越線的報警模式能夠提前進行故障預判，并且能夠有效地融合其他相關(guān)變量進行綜合預警，基于此可知，該模型在磨煤機一次風量故障預警上是有效的。