苑保利 ，林其釗，蔡振宇

(1.中國科學(xué)技術(shù)大學(xué) 能源動力，安徽 合肥 230000; 2.合肥美的電冰箱有限公司，安徽 合肥 230000)

在我國機(jī)械加工、精煉生產(chǎn)過程中，工業(yè)爐作為主要的能源供 應(yīng)設(shè)備,其重要性是顯而易見的。但由于技術(shù)、設(shè)備等原因，國內(nèi)多數(shù)煉化公司對工業(yè)爐的燃燒監(jiān)測一直處于人工監(jiān)控且不安全的階段。與國外一些發(fā)達(dá)地區(qū)相比，在燃燒監(jiān)控系統(tǒng)的研究與推廣應(yīng)用方面，我國差距較大。

近年來，在安全、效率和環(huán)境保護(hù)方面，國家對主要煉化公司的要求一直在提高，并且對人力和財力的投資逐年加大。同時，一些國內(nèi)燃燒監(jiān)測設(shè)備制造商也在逐年增多，技術(shù)與設(shè)備越來越精細(xì)化與復(fù)雜化，過去只能靠進(jìn)口的火焰監(jiān)測和控制設(shè)備，現(xiàn)能實現(xiàn)國產(chǎn)化，且更精良，可以極大提高煉化企業(yè)監(jiān)控水平。這一切，都為我國石化行業(yè)全面提升工業(yè)爐燃燒監(jiān)測及控制水平，保障煉化企業(yè)的安全生產(chǎn)，實現(xiàn)節(jié)能降耗提供了有利條件[1]。于興智[2]提出利用遠(yuǎn)程監(jiān)控及計算機(jī)圖像識別技術(shù)實現(xiàn)對排放黑煙進(jìn)行自動監(jiān)測，但這一方法無法確定煙氣含量與黑度兩者之間的對應(yīng)關(guān)系，且需依靠人來判斷顆粒物排放濃度，會引入人為誤差。劉玲玲[3]提出基于WinCC的工業(yè)爐群智能監(jiān)控系統(tǒng)，集中了智能化、信息化、網(wǎng)絡(luò)化等先進(jìn)技術(shù)，但該方法涉及的系統(tǒng)較復(fù)雜，無法廣泛、短時間推廣。馮素梅[4]提出利用主元分析(Principal Component Analysis，簡稱PCA)構(gòu)建爐管故障監(jiān)測模型，PCA方法是一種成熟有效、數(shù)據(jù)處理方法，構(gòu)建的監(jiān)測模型構(gòu)建較簡單，實用性較強(qiáng)，但PCA提取特征值波動性大，診斷結(jié)果不準(zhǔn)確[5]。

針對以上方法的不足之處，本文提出一種基于多向核熵成分分析(Multiway Kernel Entropy Component Analysis，簡稱MKECA)的故障監(jiān)測方法，應(yīng)用于工業(yè)爐管的離心鑄造過程。首先，將離心鑄造過程三維數(shù)據(jù)進(jìn)行數(shù)據(jù)預(yù)處理，利用沿批次-變量展開方法，將三維數(shù)據(jù)轉(zhuǎn)換成二維數(shù)據(jù)，再將正常批次的二維數(shù)據(jù)帶入MKECA監(jiān)測模型內(nèi)計算正常批次下的角結(jié)構(gòu)相似度統(tǒng)計量，統(tǒng)計量控制限利用核密度估計值算法獲得，本文選取99%的置信區(qū)間，最后將實時數(shù)據(jù)進(jìn)行數(shù)據(jù)預(yù)處理后，代入模型內(nèi)進(jìn)行爐管故障監(jiān)測。本文通過對乙烯裂解爐管實際生產(chǎn)數(shù)據(jù)進(jìn)行分析實驗和各種對比研究，驗證了該方法的有效性。

1 離心鑄造過程

1.1 離心鑄造過程簡介

離心鑄造過程的定義是將熔融金屬液體注入到高速旋轉(zhuǎn)的模具型筒內(nèi)，且在離心力的作用下填充并固化液態(tài)金屬的一種鑄造方法[6]。離心鑄造工業(yè)主要包括6個部件，由圖1可以看出整個工作原理與各個部件位置。

圖1 離心鑄造過程

1.2 離心鑄造過程的變量特征

工業(yè)爐生產(chǎn)過程主體由原料熔融與離心鑄造兩大部分構(gòu)成。原料熔融階段，該過程是將原材料熔融，主要是為了獲得合格的鋼水，用于澆鑄。鋼水中含有多種化學(xué)元素，比如C、Si、Mn、P、S、Cr、Ni、W、Nb、Mo、Cu、Al、Fe。

離心鑄造工業(yè)過程是工業(yè)爐第二大組成部分，該過程存在批次、時間、變量三個維度數(shù)據(jù)，與其他間歇工業(yè)工程類似，是一個復(fù)雜的工業(yè)過程，每個工業(yè)爐管的鑄造過程可視為一個批次，每個工業(yè)爐管隨時間遷移，特征變量也會隨之變化，批次、時間、變量三者之間都有復(fù)雜的聯(lián)系。本文參考文獻(xiàn)[4]，選取工業(yè)爐離心鑄造工業(yè)過程中的變量如表1所示。

表1 離心鑄造過程監(jiān)測變量

1.3 離心鑄造過程三維數(shù)據(jù)處理

針對上節(jié)所述，離心鑄造過程表現(xiàn)為三維數(shù)據(jù)，其中，I為批量數(shù)，J為變量數(shù)，K為采樣點數(shù)(時間點)，比普通連續(xù)工業(yè)過程的數(shù)據(jù)多出一維批次數(shù)據(jù)，數(shù)據(jù)結(jié)構(gòu)更加復(fù)雜。三維數(shù)據(jù)向二維擴(kuò)展，常用的兩種方法：沿批量和變量擴(kuò)展，其中沿批次方向展開原理圖如圖2所示。

圖2 沿批次數(shù)據(jù)展開

沿批次方向展開的主要優(yōu)點是，在對數(shù)據(jù)集列向矩陣進(jìn)行標(biāo)準(zhǔn)化后，能突出顯示每個批次的變化信息，提取該過程變量的平均軌跡，并且消除或減弱非線性的影響。主要缺點是，在監(jiān)視故障時，各批次數(shù)據(jù)需要完整，每個批次數(shù)據(jù)長度一致，當(dāng)進(jìn)行故障監(jiān)測時需要估計從當(dāng)前采樣點到批次結(jié)束的數(shù)據(jù)值，該過程預(yù)估值與實際值會存在一定誤差，預(yù)估值較少,誤差可以忽略不計，但預(yù)估數(shù)值較多時，整個監(jiān)測結(jié)果誤差較大。沿變量方向展開原理如圖3所示，該方法的優(yōu)點是無需考慮批次數(shù)據(jù)長度一致性，對監(jiān)測采樣數(shù)據(jù)也無需填補(bǔ)成一個完整的批次，但缺點無法突出批次之間的關(guān)系與差異，且該方法無法消除或減弱時間變量上的非線性，故對故障敏感程度低，從而實現(xiàn)的故障監(jiān)測效果較差。

圖3 沿變量數(shù)據(jù)展開

針對上述兩種三維數(shù)據(jù)預(yù)處理的不足，本文引用文獻(xiàn)[7]所提的方法，將三維數(shù)據(jù)沿批次-變量展開，展開原理圖如圖4所示。將三維數(shù)據(jù)，先沿批次方向展開，再按變量方向展開，沿批次展開可以突出批次信息，沿變量展開無需考慮完整批次數(shù)據(jù)，可以綜合前兩者優(yōu)點，彌補(bǔ)兩者缺陷，能有效的提取數(shù)據(jù)信息，為構(gòu)建高效的故障監(jiān)測模型打好良好的數(shù)據(jù)基礎(chǔ)。

圖4 沿批次-變量數(shù)據(jù)展開

2 基于MKECA離心鑄造過程故障監(jiān)測

2.1 KECA理論

KECA與KPCA選取特征數(shù)據(jù)不同在于，KECA通過Renyi熵大小來選取特征信息。假設(shè)某一概率問題中有n個事件，第i個事件Xi產(chǎn)生的概率為Pi(i=1,2…,n)，則時間轉(zhuǎn)換成Renyi熵指標(biāo)的公式為

(1)

式(1)對數(shù)函數(shù)屬于單調(diào)函數(shù)，在取值范圍內(nèi)單增或單減。因此可以單獨分析積分部分，對積分部分進(jìn)行單獨定義，公式如下：

(2)

V(p)是對概率值的平方求積分，需計算出具體值才能確定Renyi熵值，本文引入Parzen概率密度算子[8]，根據(jù)高斯函數(shù)卷積理論，再結(jié)合函數(shù)的單調(diào)性，可以將公式化簡得到

(3)

(4)

式中：I表示矩陣內(nèi)所有元素均為1的N維向量；K為樣本核矩陣。

通過式(5)可以實現(xiàn)Renyi熵的核矩陣表達(dá)，該矩陣分解公式表示為

K=EDλET

(5)

式中：Dλ為特征值λ1，…,λn組成的對角矩陣；E則為由特征向量e1,…,en組成的矩陣。

因此，式(3)進(jìn)一步可表示為

(6)

式中：1為n×1的向量，式(6)可以表示為每個特征值，和特征值對應(yīng)的特征向量對Renyi熵產(chǎn)生不同的貢獻(xiàn)，得到的轉(zhuǎn)換數(shù)據(jù)：

(7)

式中：Φeca為KECA映射；D是特征值為λ1,λ2,…,λn的對角矩陣；E=[e1,e2,…en]。

2.2 角結(jié)構(gòu)相似度統(tǒng)計量

由于采樣數(shù)據(jù)經(jīng)過KECA轉(zhuǎn)換后，數(shù)據(jù)信息內(nèi)有較好的角結(jié)構(gòu)信息，因此本文引入一種新型故障監(jiān)測統(tǒng)計量-角結(jié)構(gòu)相似度統(tǒng)計量CV，用于故障監(jiān)控，定義如下[9]：

(8)

角結(jié)構(gòu)相似度統(tǒng)計量CV是計算矩陣與矩陣之間的相似度，對比矩陣則是其他矩陣的平均值。相似度越高，CV值就越高，但為在圖紙和結(jié)果表現(xiàn)更直觀，將該數(shù)據(jù)值用1減去相似度值得到統(tǒng)計量值。因此在正常批次數(shù)據(jù)下，數(shù)據(jù)與原歷史數(shù)據(jù)相似度高，就會低于控制限，反之就會高于控制限，表示為故障，可以實現(xiàn)高效故障監(jiān)測。為避免選取控制限方法而引入的誤差，本文采用核密度估計計算CV統(tǒng)計量的控制限(CK)，該方法在計算控制限時，無須假設(shè)數(shù)據(jù)遵從高斯分布，數(shù)據(jù)處理更加偏于實際生產(chǎn)。

3 實驗分析與討論

3.1 實驗數(shù)據(jù)來源

為驗證本文提出方法的有效性，共收集到32根實際工業(yè)爐生產(chǎn)的爐管數(shù)據(jù)，其中24根合格,8根不合格。本文實驗分析主要從兩部分分析，一是主元個數(shù)選取，二是整體監(jiān)測效果對比。為對比本文提出方法的有效性，兩部分都引入MKPCA方法作對比實驗。

3.2 主元個數(shù)選擇MKECA與MKPCA對比分析

在構(gòu)建MKECA故障監(jiān)測模型時，數(shù)據(jù)主元個數(shù)的選取與監(jiān)測效果有很大關(guān)系。數(shù)據(jù)主元選取過多，數(shù)據(jù)處理時間長，虛假、雜質(zhì)信息過多，影響監(jiān)測效果；主元個數(shù)過少，故障特征信息表現(xiàn)少，削弱故障表現(xiàn)能力，也會影響監(jiān)測效果，因此故障監(jiān)測效果與MKECA處理后的主元個數(shù)選取至關(guān)重要。本文將正常三維數(shù)據(jù)預(yù)處理成二維數(shù)據(jù)后，代入MKECA監(jiān)測模型內(nèi)，獲得該模型不同主元個數(shù)下的控制限和原始數(shù)據(jù)矩陣，再將預(yù)處理后的驗證數(shù)據(jù)代入監(jiān)測模型內(nèi)，計算在不同主元個數(shù)模型下，各個監(jiān)測模型的診斷率和誤報率，本文規(guī)定誤報率為合格爐管被誤診出故障的概率，診斷率為不合格爐管被診斷出故障的概率。該過程驗證的合格數(shù)據(jù)和故障數(shù)據(jù)都是同一組，為對比研究，本文增加MKPCA進(jìn)行對比驗證。表2和圖5列出了驗證數(shù)據(jù)在不同主元個數(shù)時，MKECA與MKPCA兩種監(jiān)測模型誤報率診斷率的結(jié)果。

表2 驗證數(shù)據(jù)選擇主元個數(shù)的結(jié)果 %

圖5 MKECA主元選取個數(shù)診斷率與誤報率結(jié)果圖

由表2數(shù)據(jù)可以看出，在MKPCA數(shù)據(jù)結(jié)果中，統(tǒng)計量T2在主元個數(shù)為6個時，效果最佳，誤報率在同模型下也最低，SPE的統(tǒng)計量在主元個數(shù)為7時綜合檢測效果最好，但MKPCA是綜合T2與SPE兩種統(tǒng)計量來算，因此綜合兩者分析，主元個數(shù)為6，效果最佳；相比MKPCA，MKECA監(jiān)測效果由圖5可以看出，雖然當(dāng)主元個數(shù)為7時，診斷率最高，但是綜合誤報率來看，主元個數(shù)為6時，MKECA的CV統(tǒng)計量監(jiān)測效果最好，因此當(dāng)取主元個數(shù)為6時，整體故障監(jiān)測效果最佳，構(gòu)建的該模型效果最好。

3.3 MKECA與MKPCA故障監(jiān)測對比研究

為研究MKECA在離心鑄造過程中的故障監(jiān)測能力，利用16組正常數(shù)據(jù)構(gòu)建MKECA離線監(jiān)測模型，再將8組合格數(shù)據(jù)集與8組故障數(shù)據(jù)集代入構(gòu)建好的故障監(jiān)測模型內(nèi)，實現(xiàn)在線故障監(jiān)測誤報率和故障診斷率。為對比研究本文提出方法的有效性，利用MKPCA進(jìn)行對比，兩種模型的主元個數(shù)都選6。

表3和圖6～圖11是取99%的置信區(qū)間取CV、T2和SPE的控制限，表2是8組合格與8組故障數(shù)據(jù)的整體監(jiān)測效果，從數(shù)據(jù)中可以看出，MKECA誤報率為6.75%，誤報率相比MKPCA的T2與SPE效果最低，說明MKECA的誤報最少，效果最好，反之診斷率須越高越好，MKECA的CV診斷率整體效果最高，說明識別故障效率最好，因此驗證了MKECA統(tǒng)計監(jiān)測效果最好，驗證了該方法的有效性。圖6～圖11展現(xiàn)的是某個合格批次和不合格批次誤報率、合格率效果圖，從圖6～圖11可以直觀看出MKECA的監(jiān)測效果更佳。

圖11 某不合格數(shù)據(jù)集MKPCA診斷率SPE效果圖

表3 故障監(jiān)測誤報率與故障診斷率對比 %

圖6 某合格數(shù)據(jù)集MKECA誤報率效果圖

圖7 某合格數(shù)據(jù)集MKPCA誤報率T2效果圖

圖8 某合格數(shù)據(jù)集MKPCA誤報率SPE效果圖

圖9 某不合格數(shù)據(jù)集MKECA診斷率CV效果圖

4 結(jié) 論

針對目前離心鑄造工業(yè)爐管生產(chǎn)過程監(jiān)控方法復(fù)雜、耗費(fèi)高等特點，本文提出基于MKECA進(jìn)行離心鑄造過程故障監(jiān)測方法。工業(yè)爐管的離心鑄造過程是一種間歇工業(yè)過程，間歇過程表現(xiàn)是三維數(shù)據(jù)集，首先將三維數(shù)據(jù)集沿批次-變量二維數(shù)據(jù)展開，再利用MKECA建立離線故障監(jiān)測模型，得到CV控制限和歷史正常矩陣數(shù)據(jù)，再代入驗證數(shù)據(jù)計算驗證數(shù)據(jù)的CV統(tǒng)計量，判斷是否超出控制限，以此方法進(jìn)行故障監(jiān)測。本文利用乙烯裂解爐管實際生產(chǎn)過程數(shù)據(jù)，代入MKPCA監(jiān)測模型內(nèi)與本文提出的MKECA監(jiān)測模型進(jìn)行對比驗證，從結(jié)果來看，本文的方法誤報率6.75%，故障診斷率是94.53%，遠(yuǎn)高于MKPCA的方法，驗證了該方法的有效性與可行性，為實際生產(chǎn)提供一條可行的故障監(jiān)測思路。

圖10 某不合格數(shù)據(jù)集MKPCA診斷率T2效果圖