張龍，陳宸，韓寧，王亞慧

(1. 北京林業(yè)大學(xué) 工學(xué)院，北京 100083； 2. 美國德州大學(xué)Dallas分校 電子工程系，Richardson 75080； 3. 北京建筑大學(xué) 電氣與信息工程學(xué)院，北京 100044)

隨著城市化進(jìn)程的加速，高層和超高層建筑日益增加，人們對于建筑物安全和舒適度的要求也越來越高。在整個建筑物中，建筑電氣是關(guān)鍵技術(shù)之一，它包括了照明系統(tǒng)、供配電系統(tǒng)、動力設(shè)備系統(tǒng)、辦公及管理自動化等主要內(nèi)容。不同子系統(tǒng)間的相互關(guān)聯(lián)越來越緊密，同時也使建筑電氣故障發(fā)生幾率增加。故障智能診斷技術(shù)已在電力系統(tǒng)、機(jī)械等領(lǐng)域發(fā)展得比較成熟[1-2]，但在建筑電氣方面還處于空白階段，基本依靠人工檢測查找故障原因。目前關(guān)于故障診斷的算法多采用人工神經(jīng)網(wǎng)絡(luò)算法、支持向量機(jī)算法(SVM)、小波變換法等。壓縮感知[3](compressive sensing，CS)方法是近年來興起的一種新的理論算法，已在人臉識別分類和圖像處理上有了廣泛應(yīng)用[4-6]。該理論基于線模型，其核心假設(shè)是信號的稀疏性，只要信號是稀疏或可壓縮的，低維觀測信號就能很好地恢復(fù)到高維原始信號。將其用于分類，也就是用訓(xùn)練樣本本身作為基元素去表示測試樣本，用與測試樣本相同類的訓(xùn)練樣本的線性組合來表示輸入的待識別的樣本，達(dá)到分類目的[7]。

1 建筑電氣系統(tǒng)故障診斷原理

建筑電氣系統(tǒng)常見故障類型、表現(xiàn)特征及危害如表1所示，具體故障診斷的系統(tǒng)結(jié)構(gòu)如圖1所示。

表1 建筑電氣系統(tǒng)常見故障

建筑電氣系統(tǒng)故障診斷實(shí)質(zhì)上是基于征兆集/故障集的映射模式，即故障發(fā)生時的征兆提取和故障狀態(tài)判斷。由于建筑電氣系統(tǒng)故障種類繁多且發(fā)生機(jī)率隨機(jī)，所以本文以建筑電氣故障模擬實(shí)驗(yàn)平臺為研究對象，對建筑物中常見的電氣故障，包括絕緣故障、接地系統(tǒng)故障、配電系統(tǒng)接地故障等進(jìn)行故障診斷研究，按照不同故障工作狀態(tài)的診斷目的和對象，選擇便于診斷的狀態(tài)信號(電流、電壓及電阻值)，將傳感器加裝于系統(tǒng)關(guān)鍵回路，通過數(shù)據(jù)采集器收集故障時異常信號，提取故障特征，輸入處理后的數(shù)據(jù)，經(jīng)過故障診斷算法判別輸出故障類別及相應(yīng)位置，提示報警信息，最后根據(jù)問題所在提出控制措施和維修策略。

圖1 建筑電氣實(shí)驗(yàn)平臺故障識別系統(tǒng)框圖Fig.1 Fault identification system block diagram of the building electrical experiment platform

2 壓縮感知理論基礎(chǔ)

壓縮感知理論是由Donoho與Candes等提出的理論框架，具有廣闊的應(yīng)用前景，受到各個領(lǐng)域研究人員的廣泛關(guān)注。在壓縮感知中，信號的稀疏表示和重構(gòu)是整個理論的核心。

在此基礎(chǔ)上考慮信號重構(gòu)問題，如果構(gòu)造一個觀測矩陣φ∈Rm×n(m遠(yuǎn)小于n)，原始信號為x∈Rn×1，x在φ上的線性測量值為y∈Rm×1，即

y=φx=φΨS

(1)

式(1)展示了原始信號x在觀測矩陣φ變換下的線性投影，現(xiàn)在考慮如何從信號y中重構(gòu)出原始信號x來，理論證明，可以通過對測量值y的最優(yōu)l0范數(shù)問題的求解來實(shí)現(xiàn)重構(gòu)[8]，獲得x的精確或近似逼近解：

(2)

但由于信號y的維數(shù)遠(yuǎn)遠(yuǎn)低于原始信號x的維數(shù)(m?n)，所以式(1)的解有無窮多個，是一個NP-hard問題，考慮到S為稀疏向量，通過合理選擇觀測矩陣φ和稀疏矩陣ψ，通?？蓪⑹?2)轉(zhuǎn)換為求解l1范數(shù)下的最優(yōu)問題：

(3)

求解l1優(yōu)化問題，可以利用內(nèi)點(diǎn)法、梯度投影法、二階圓錐規(guī)劃、匹配追蹤法等方法求解[7]。

3 壓縮感知故障診斷方法

故障診斷的過程實(shí)質(zhì)上就是一個分類的過程，利用故障時的異常信號，提取故障特征，通過算法判別分類各種故障類型。

3.1 訓(xùn)練樣本的組成

假設(shè)需要對k類故障進(jìn)行分類，每個故障樣本維數(shù)為p維，組成一個p×1維的列向量v，第i類故障的訓(xùn)練樣本數(shù)為ni(i=1,2,…,k)，組成訓(xùn)練樣本矩陣如式(4)所示：

(4)

式中：vi,j為第i類故障的第j個訓(xùn)練樣本，Ai為第i類故障的訓(xùn)練樣本矩陣。

由于樣本所屬類別i未知，將所有k個類別的n個訓(xùn)練樣本拼接在一起，組成完備訓(xùn)練樣本矩陣A：

3.2 測試樣本的稀疏分解

若待分類故障的測試樣本y∈Rp屬于第i類，則y可以通過第i個故障訓(xùn)練樣本集合線性表示：

y=ai,1vi,1+ai,2vi,2+…+ai,nivi,ni

式中：ai,j為權(quán)重系數(shù)。

當(dāng)給定一個故障測試數(shù)據(jù)y時，該樣本所屬類別是未知的，需要求出它是樣本集中哪種故障。因此通過完備矩陣A來線性表示出待分類的故障y：

y=a1,1v1,1+…+a1,n1v1,n1+…+ai,1vi,1+…+

ai,nivi,ni+…+ak,1vk,1+…+ak,nkvk,nk

(5)

y=AX

(6)

利用式(5)可以解出矩陣X：

因此X的系數(shù)理論上只有ai不為0，其他k-1個系數(shù)都為0，可見X是一個稀疏向量，可看作是測試樣本y的稀疏分解。

3.3 測試樣本分類

若要求對給定的未知測試樣本進(jìn)行歸類，那么只要根據(jù)式(6)對于每一個y解出稀疏向量X，結(jié)合已知的完備矩陣A，X上只有與該被測樣本有關(guān)的ni個系數(shù)為非0值，即可知道待分類故障的類別。實(shí)際問題中通過第3節(jié)中介紹的求解l1最小化問題，獲得x的精確或近似逼近解，但實(shí)際求解結(jié)果并非如理論所述，x的非零元素將會散布于很多類間，為了通過X的值完成分類工作，需采用以下分類函數(shù)[9]：

(7)

4 稀疏表示分類算法故障診斷實(shí)驗(yàn)

4.1 建筑電氣故障模擬實(shí)驗(yàn)平臺

建筑電氣故障模擬實(shí)驗(yàn)平臺是本文實(shí)驗(yàn)室研究階段的重要試驗(yàn)對象，其原產(chǎn)于德國，集合了住宅建筑物內(nèi)部低壓配電系統(tǒng)中常見的低壓電氣裝置，如：斷路器、熔斷器、RCD(剩余電流保護(hù)器)、單向插座、三相插座等，如圖2。

圖2 建筑電氣系統(tǒng)測試平臺MA2067Fig.2 Experimental platform of building electrical system-MA2067

該實(shí)驗(yàn)平臺的內(nèi)部結(jié)構(gòu)如圖3所示，電源供電為220 V、50 Hz交流電，由變壓器轉(zhuǎn)變?yōu)?5V直流輸出，為弱電保護(hù)板供電。弱電保護(hù)板對強(qiáng)電系統(tǒng)中的單相和三相系統(tǒng)進(jìn)行保護(hù)。強(qiáng)電系統(tǒng)是該實(shí)驗(yàn)臺主體，系統(tǒng)通過故障設(shè)置面板上的22個開關(guān)的斷開閉合對強(qiáng)電系統(tǒng)中四大類阻值故障、22個故障位置進(jìn)行模擬故障設(shè)置，斷開即為通路，閉合即為相應(yīng)部位故障發(fā)生。

圖3 建筑電氣故障模擬實(shí)驗(yàn)平臺內(nèi)部結(jié)構(gòu)Fig.3 Physical model of electrical test platform

4.2 故障特征量的選擇與故障分類

根據(jù)實(shí)驗(yàn)平臺能夠模擬的實(shí)際住宅建筑物中的常見故障，故障類型可分為線路阻抗故障(E1)、連續(xù)性故障(E2)、接地電阻異常(E3)、絕緣電阻過小(E4)共4種，再加上正常狀態(tài)(E5)，所以本文的診斷狀態(tài)共有5類。通過采集實(shí)驗(yàn)平臺10個不同測試位置的故障信息值(電阻值)作為算法輸入的特征分量，位置信息如表2所示。

表2 建筑電氣實(shí)驗(yàn)平臺故障特征與對應(yīng)位置

4.3 實(shí)驗(yàn)方案設(shè)計(jì)

為驗(yàn)證本文提出的故障診斷分類方法的有效性，設(shè)計(jì)了以下3個實(shí)驗(yàn)，分別采用支持向量機(jī)、基于稀疏表達(dá)分類算法的l1分類器和l2分類器。在壓縮感知理論中，對于信號的重建，要求稀疏矩陣構(gòu)成正交基底。利用稀疏表達(dá)(sparse representation)來做分類(見式(6))，由于完備矩陣A通常是奇異矩陣(不可逆)，所以求解系數(shù)矩陣X的時候，需要采用正規(guī)化手段。本文采用了l1和l22種正規(guī)化方法分別來求解系數(shù)矩陣X，這2種方法的不同點(diǎn)在于它的目標(biāo)函數(shù)[10]。

式中：L采用Tikhonov正則化矩陣，λ則為正則化參數(shù)。

本文通過實(shí)驗(yàn)平臺實(shí)測收集了樣本數(shù)據(jù)共55組，5種狀態(tài)模式(線路阻抗故障、連續(xù)性故障、接地電阻異常、絕緣電阻過小、正常)，每個樣本含10個不同位置故障信息特征分量。3次實(shí)驗(yàn)時，每一類故障均隨機(jī)選取一個樣本作為測試樣本，剩下的50個樣本為訓(xùn)練樣本，重復(fù)50組實(shí)驗(yàn)，然后取平均值作為最終分類的準(zhǔn)確度，并計(jì)算診斷運(yùn)行時間。本文實(shí)驗(yàn)都是運(yùn)行在2.13 GHz的雙核處理器上。

4.4 診斷結(jié)果及分析

根據(jù)實(shí)驗(yàn)方案及步驟，此時式(6)中的矩陣A的維度為10×50，測試樣本y的維度為10×1。實(shí)驗(yàn)結(jié)果如圖4、圖5。圖4是由l2分類器求解的x的稀疏系數(shù)(本圖采用屬于線路阻抗故障(E1)的測試樣本)，圖5即是由式(7)計(jì)算出的最終5個殘差項(xiàng)。

圖4 稀疏表示系數(shù)xFig.4 The sparse representation coefficients x

圖5 殘差項(xiàng)ri(y)Fig.5 The residuals ri(y)

從圖4中可以明顯看出，當(dāng)輸入第1類故障數(shù)據(jù)用于測試時(不在訓(xùn)練樣本中)，所得到的稀疏表示向量的確是非常稀疏的，最大系數(shù)與第1類故障訓(xùn)練數(shù)據(jù)相對應(yīng)，則從直觀上判斷該測試樣本應(yīng)屬于第1類故障。圖5則是通過數(shù)學(xué)計(jì)算方法來精確判斷待測故障屬于哪一類別。類別1的殘差值最小，因此該測試樣本屬于第1類。

表3展示了l1分類器、l2分類器以及SVM在相同實(shí)驗(yàn)條件和方法下的故障診斷結(jié)果。

表3 實(shí)驗(yàn)結(jié)果對比

通過表3中數(shù)據(jù)可以看出，這3種分類算法對測試樣本的識別率(準(zhǔn)確度)都是令人滿意的。但是l1分類器在運(yùn)算時間上遠(yuǎn)遠(yuǎn)大于后2種方法，這在電氣系統(tǒng)運(yùn)行過程中將造成來不及使保護(hù)器動作的安全隱患，而同樣是基于稀疏表示的l2分類器不論是從分類準(zhǔn)確率和診斷時間上都優(yōu)于l1分類器，可以認(rèn)為本文提出的基于壓縮感知理論的稀疏表示分類算法的診斷效果與支持向量機(jī)SVM算法不相上下，也可以作為建筑電氣系統(tǒng)故障診斷的一種新方法加以推廣。

5 結(jié)束語

本文以建筑電氣系統(tǒng)故障模擬實(shí)驗(yàn)平臺數(shù)據(jù)為研究對象，提出了基于壓縮感知理論的稀疏表示分類算法作為故障診斷方法引入建筑電氣系統(tǒng)。通過實(shí)驗(yàn)對比，分析了l1分類器、l2分類器以及SVM在故障診斷分類中的性能。結(jié)果表明，本文提出的算法分類準(zhǔn)確率達(dá)到了96.4%，診斷運(yùn)行時間0.260 1 s，得到了較好的效果。稀疏表示分類算法不僅是對現(xiàn)有故障診斷方法的有益補(bǔ)充，也可將其進(jìn)一步應(yīng)用于智能型建筑電氣故障診斷系統(tǒng)中，對于及早發(fā)現(xiàn)和排除故障發(fā)揮重要作用。

參考文獻(xiàn):

[1]邵曉非, 寧媛, 劉耀文, 等. 電力系統(tǒng)故障診斷方法綜述與展望[J]. 工業(yè)控制計(jì)算機(jī), 2012, 25(12): 4-7.

SHAO Xiaofei, NING Yuan, LIU Yaowen, et al. Review and prospects of fault diagnosis in power system[J]. Industrial Control Computer, 2012, 25(12): 4-7.

[2]黃建, 胡曉光, 鞏玉楠, 等. 高壓斷路器機(jī)械故障診斷專家系統(tǒng)設(shè)計(jì)[J]. 電機(jī)與控制學(xué)報, 2011, 15(10): 43-49.

HUANG Jian, HU Xiaoguang, GONG Yunan, et al. Machinery fault diagnosis expert system for high voltage circuit breaker[J]. Electric Machines and Control, 2011, 15(10): 43-49.

[3]DONOHO D. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289-1306.

[4]CHEN C, TRAMEL E W, FOWLER J E. Compressed-sensing recovery of images and video using multihypothesis predictions[C]//Proceedings of the 45th Asilomar Conference on Signals, Systems, and Computers. Pacific Grove, CA, 2011: 1193-1198.

[5]CHEN C, FOWLER J E. Single-image super-resolution using multihypothesis prediction[C]//Proceedings of the 46th Asilomar Conference on Signals, Systems, and Computers. Pacific Grove, CA, 2012： 608-612.

[6]CHEN C, LI W, TRAMEL E W, et al. Reconstruction of hyperspectral imagery from random projections using multihypothesis prediction[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014， 52(1)： 365-374.

[7]韓安太, 郭小華, 廖忠, 等. 基于壓縮感知理論的農(nóng)業(yè)害蟲分類方法[J]. 農(nóng)業(yè)工程學(xué)報, 2011, 27(6): 203-207.

HAN Antai, GUO Xiaohua, LIAO Zhong, et al. Classification of agricultural pests based on compressed sensing theory[J]. Transactions of the CSAE, 2011, 27(6): 203-207.

[8]蔡騁, 張明, 朱俊平.基于壓縮感知理論的雜草種子分類識別[J].中國科學(xué)， 2010, 40: 160-172.

CAI Cheng, ZHANG Ming, ZHU Junping. Weed seeds classification based on compressive sensing theory[J]. Science China， 2010, 40: 160-172.

[9]WRIGHT J, YANG A, GANESH A, et al, Robust face recognition via sparse representation[J]. IEEE Trans Pattern Anal, 2009, 31: 210-227.

[10]CHEN Chen, LIU Kui , KEHTARNAVAZ N. Real-time human action recognition based on depth motion maps[J]. Journal of Real Time Image Processing, 2013， doi：10.1007/s11554-013-0370-1.

[11]王亞慧,張龍,韓寧.建筑電氣系統(tǒng)故障診斷方法研究[J].計(jì)算機(jī)仿真, 2014, 31(2): 436-440.

WANG Yahui, ZHANG Long, HAN Ning. The research on fault diagnosis of building electrical system[J]. Computer Simulation, 2014, 31(2): 436-440.