王 丹 楊煜普 屈衛(wèi)東

(上海交通大學(xué)電子信息與電氣工程學(xué)院自動(dòng)化系，上海 200240)

近年來，隨著現(xiàn)代化工及冶金等工業(yè)過程日益大規(guī)?；蛷?fù)雜化，工業(yè)過程的安全問題越來越受到人們的關(guān)注。基于多元統(tǒng)計(jì)分析的故障檢測與診斷方法也成為近年來故障檢測與診斷領(lǐng)域的研究熱點(diǎn)，并在工業(yè)過程中成功應(yīng)用[1～3]。偏最小二乘(PLS)技術(shù)能夠根據(jù)正常工況的生產(chǎn)數(shù)據(jù)，準(zhǔn)確捕捉質(zhì)量變量與過程變量之間的關(guān)系，對生產(chǎn)工況進(jìn)行有效監(jiān)測，且PLS統(tǒng)計(jì)檢測技術(shù)不依賴于過程機(jī)理模型，訓(xùn)練時(shí)不需要故障樣本，能夠彌補(bǔ)其他統(tǒng)計(jì)方法(例如PCA)無法考慮過程變量對質(zhì)量變量影響的不足，因此近年來在化工過程的質(zhì)量控制及在線檢測等方面得到了廣泛研究和應(yīng)用[4～6]。但是PLS方法無法反映過程的動(dòng)態(tài)時(shí)序特性，這在一定程度上影響了它的故障檢測準(zhǔn)確率?？深A(yù)測元分析[7](Forecastable Component Analysis，F(xiàn)oreCA)作為一種新的統(tǒng)計(jì)信號處理方法克服了這個(gè)不足。可預(yù)測元分析是一種全新的用于多變量時(shí)序相關(guān)信號的特征提取方法，它能從已有的數(shù)據(jù)中捕捉到系統(tǒng)的動(dòng)態(tài)特性，并以此來預(yù)測系統(tǒng)運(yùn)行變化的趨勢，因此所提取的特征能從本質(zhì)上描述工業(yè)過程。

筆者將可預(yù)測元分析方法與偏最小二乘法回歸方法相結(jié)合并用于故障檢測，通過將樣本映射到可預(yù)測子空間，使用最小二乘回歸，進(jìn)一步提高了模型的預(yù)測性能，同時(shí)構(gòu)造CUSUM和SPE統(tǒng)計(jì)量對系統(tǒng)進(jìn)行監(jiān)控，這樣能夠較好地檢測均值偏差在兩倍標(biāo)準(zhǔn)差以下的故障。該方法克服了傳統(tǒng)偏最小二乘法無法反映過程時(shí)序特性的不足，能夠預(yù)測系統(tǒng)運(yùn)行變化的趨勢，反映出系統(tǒng)的動(dòng)態(tài)特性，因此能夠提升故障檢測的準(zhǔn)確率。

1 基本算法①

1.1 可預(yù)測元分析

可預(yù)測元分析的基本思想是假設(shè)矩陣X∈Rn×m，其中n為樣本個(gè)數(shù)，m為變量個(gè)數(shù)，通過線性變換WT∈Rk×n，可得：

S=WTX

(1)

其中W為由可預(yù)測元列向量組成的可預(yù)測元矩陣，S為得分矩陣，F(xiàn)oreCA需要解決的問題即由觀測矩陣X估計(jì)S和W。

γy(k)=E(yt-μy)(yt-k-μy)T,k∈R

(2)

其中k為時(shí)延。

定義單變量平穩(wěn)過程的譜密度為對其自協(xié)方差函數(shù)的傅里葉變換，得：

(3)

(4)

由文獻(xiàn)[7]可知，一個(gè)平穩(wěn)過程的熵越大越難被預(yù)測，且白噪聲無法被預(yù)測，可得：

Hs,a(yt)≤Hs,a(白噪聲)

(5)

因此可定義平穩(wěn)過程的可預(yù)測度為：

(6)

對于多變量二階平穩(wěn)過程Xt，考慮線性變換yt=wTXt，其中w(w∈Rn)是式(1)中W的列向量，即可預(yù)測元，此時(shí)yt可以看成是一個(gè)單變量的二階平穩(wěn)過程。文獻(xiàn)[7]給出了ForeCA的最優(yōu)化問題：

(7)

s.t.wTΣXw=1

對式(7)進(jìn)行求解，首先使用加權(quán)交疊平均譜估計(jì)法對隨機(jī)過程進(jìn)行譜密度估計(jì)[8]，再使用EM-Like算法求取可預(yù)測元[7]。文獻(xiàn)[7]給出了此算法的詳細(xì)步驟，通過此算法可以得出一組按照可預(yù)測度由高到低順序排列的可預(yù)測元，即可得到式(1)中的可預(yù)測元矩陣wT。

1.2 偏最小二乘法

給定輸入矩陣X∈Rn×N包含n個(gè)樣本，每個(gè)樣本N個(gè)過程變量，輸出矩陣Y∈Rn×M包含n個(gè)樣本，每個(gè)樣本M個(gè)質(zhì)量變量。PLS通過隱變量對兩個(gè)數(shù)據(jù)塊的關(guān)系進(jìn)行建模，它將n×N零均值矩陣X和n×M零均值矩陣Y分解為：

(8)

(9)

式中Ek、Fk——擬合誤差矩陣；

P——X的負(fù)載矩陣；

Q——Y的負(fù)載矩陣；

T——得分矩陣，T=[t1,…,tk]；

在PLS模型中，負(fù)載向量和得分向量通過最大化解釋各自的信息，同時(shí)也使X與Y的相關(guān)程度最大來求得。最常見的計(jì)算PLS模型的算法是Nipals算法，Y的預(yù)測回歸方程為：

(10)

其中，BPLS是PLS回歸系數(shù)矩陣，權(quán)重矩陣M是由Nipals算法定義的，T=XM。

在復(fù)雜的多變量系統(tǒng)中，PLS算法將自變量X和因變量Y看成是具有線性關(guān)系的數(shù)據(jù)矩陣。沒有逐個(gè)對變量判斷其留取與舍棄，而是利用信息分解的思路將顯變量系統(tǒng)中的信息重新組合，綜合篩選，提取出既能最大程度解釋自變量信息，又能最大程度反映自變量與因變量間線性關(guān)系的互相正交的綜合變量(隱變量)。PLS用獨(dú)立的隱變量進(jìn)行建模、預(yù)測，使得該方法可以廣泛應(yīng)用于數(shù)據(jù)不完整、變量間存在多重相關(guān)性的場合。

2 基于ForePLS的故障檢測模型

2.1 CUSUM統(tǒng)計(jì)量

工業(yè)過程中存在很多慢漂移的故障，為了檢測這種微小的變化，筆者用CUSUM統(tǒng)計(jì)量對其進(jìn)行檢測?；贑USUM統(tǒng)計(jì)量的表格累加法為了檢測樣本均值向上和向下漂移，定義了兩個(gè)統(tǒng)計(jì)量，即：

SH(i)=max[0,xi-(μ0+K)+SH(i-1)],SH(0)=0

(11)

SL(i)=max[0,(μ0-K)-xi+SL(i-1)],SL(0)=0

(12)

(13)

其中μ0是樣本實(shí)際的均值，xj為第j個(gè)樣本值，筆者用訓(xùn)練樣本均值代替。K為參考值，一般取0.5Δ，Δ為期望檢測出的偏差，取值在[0.5σ,2.0σ]內(nèi)。其控制限為5倍的標(biāo)準(zhǔn)差[9]。

2.2 SPE統(tǒng)計(jì)量

首先選取一段正常工況下的觀測數(shù)據(jù)X(X∈Rn×N)，其中n為變量個(gè)數(shù)，N為采樣點(diǎn)數(shù)，對其運(yùn)用ForeCA算法，得可預(yù)測元矩陣：

WT=[w1,w2,…,wn]T∈RN×N

(14)

(15)

過程殘差可表示為：

(16)

(17)

SPE統(tǒng)計(jì)量的控制限用核密度估計(jì)確定，具體參見文獻(xiàn)[11]。

2.3 基于ForePLS的故障檢測步驟

基于ForePLS的故障檢測分為兩個(gè)階段——離線訓(xùn)練階段和在線檢測階段。

離線訓(xùn)練階段。首先采集正常工況下的訓(xùn)練數(shù)據(jù)X，對其進(jìn)行預(yù)處理后，使用ForeCA算法提取出可預(yù)測主元矩陣W，然后在可預(yù)測子空間進(jìn)行PLS回歸，再計(jì)算訓(xùn)練數(shù)據(jù)在可預(yù)測子空間的CUSUM統(tǒng)計(jì)量和SPE統(tǒng)計(jì)量，最后計(jì)算兩個(gè)統(tǒng)計(jì)量的控制限——H和SPEα。

在線檢測階段。首先根據(jù)實(shí)時(shí)采集的未知狀態(tài)的數(shù)據(jù)集，將此可預(yù)測模型運(yùn)用于在線數(shù)據(jù)，分別計(jì)算每個(gè)樣本數(shù)據(jù)的CUSUM和SPE統(tǒng)計(jì)量，最后比較兩個(gè)統(tǒng)計(jì)量與其對應(yīng)控制限的大小，通過比較確定系統(tǒng)是否發(fā)生故障。如果檢驗(yàn)結(jié)果在控制限以內(nèi)，則說明目前系統(tǒng)工作在可預(yù)測模型所預(yù)測的變化范圍之內(nèi)，即系統(tǒng)工作正常；反之，則說明目前系統(tǒng)的工作狀態(tài)已經(jīng)偏離可預(yù)測模型所預(yù)測的變化范圍，判斷系統(tǒng)已經(jīng)出現(xiàn)了故障。

3 TE實(shí)驗(yàn)平臺故障分析

TE實(shí)驗(yàn)平臺是Downs和Vogel根據(jù)Eastman化學(xué)公司的世界工藝流程做了少許修改于1993年提出的[12]，其中包含21個(gè)預(yù)設(shè)故障。TE過程由連續(xù)攪拌式反應(yīng)釜、分凝器、氣液分離塔、汽提塔、再沸器及離心式壓縮機(jī)等多個(gè)操作單元組成，其流程如圖1所示。

圖1 TE流程

TE過程共有A、C、D、E 4種氣體進(jìn)料，G和H兩種反應(yīng)產(chǎn)物，F(xiàn)一種副產(chǎn)品。系統(tǒng)中存在的化學(xué)反應(yīng)如下：

以上各式中，g代表氣體，liq代表液體。所有的反應(yīng)都是不可逆放熱反應(yīng)，反應(yīng)速度取決于溫度和反應(yīng)物的氣相濃度。

TE模型用于訓(xùn)練的樣本數(shù)據(jù)為500個(gè)52維向量，用于測試的樣本數(shù)據(jù)為960個(gè)52維向量，其中故障從第161個(gè)樣本點(diǎn)開始引入。筆者選擇過程中的G和H(即MEAS35和MEAS36)作為ForePLS模型的質(zhì)量變量Y；選取22個(gè)過程變量MEAS1～22和11個(gè)操作變量MV1～11作為X。采用ForePLS模型對TE過程的反應(yīng)產(chǎn)物G的含量的預(yù)測結(jié)果如圖2所示，可以看出ForePLS有很好的預(yù)測能力。

圖2 故障10發(fā)生時(shí)產(chǎn)品中組分G的含量曲線

下面以隨機(jī)變化故障中典型的故障IDV(10)為例加以分析。故障IDV(10)發(fā)生時(shí)，供料C的溫度產(chǎn)生了隨機(jī)變化。為了驗(yàn)證ForePLS的有效性，將其與PCA和PLS兩種方法進(jìn)行對比。實(shí)驗(yàn)中，F(xiàn)orePLS的隱變量個(gè)數(shù)為6，PCA的主元個(gè)數(shù)為15，PLS的隱變量個(gè)數(shù)為9，期望檢測到的均值偏離為0.5倍的標(biāo)準(zhǔn)差。圖3顯示了PCA、PLS和ForePLS 3種方法對故障IDV(10)的檢測效果。可以看出，PCA的T2統(tǒng)計(jì)量和SPE統(tǒng)計(jì)量的準(zhǔn)確率分別為45.6%和53.9%；PLS的兩個(gè)統(tǒng)計(jì)量的檢測準(zhǔn)確度都較低，分別為18.8%和27.8%；ForePLS的CUSUM和SPE統(tǒng)計(jì)量的準(zhǔn)確率為96.5%和52.9%。由此說明，筆者所提出的基于ForePLS的故障檢測方法檢測隨機(jī)變化的故障準(zhǔn)確率比PCA和PLS方法更好。

圖3 IDV(10)發(fā)生時(shí)PCA、PLS和ForePLS方法的故障檢測性能比較

4 結(jié)束語

介紹了一種基于可預(yù)測元分析和最小二乘回歸法相結(jié)合的故障檢測方法。該方法克服了傳統(tǒng)最小二乘法無法反映過程時(shí)序特性的不足，能夠有效預(yù)測系統(tǒng)運(yùn)行變化的趨勢，反映出系統(tǒng)的動(dòng)態(tài)特性。通過檢測可預(yù)測空間上的CUSUM統(tǒng)計(jì)量和SPE統(tǒng)計(jì)量，以達(dá)到檢測慢漂移等微小故障和隨機(jī)變化故障的目的。在TE模型上的仿真表明：該方法比傳統(tǒng)的PCA、PLS方法檢測精度更高，效果更好。