(遼寧工程技術(shù)大學 電氣與控制工程學院，遼寧 葫蘆島 125105)

電氣設(shè)備中的線路在運行中如果存在絕緣碳化，外界引起的空氣電離、短路等情況，均有可能導致故障電弧，而故障電弧是引發(fā)電氣火災的主要原因。據(jù)統(tǒng)計，2016年我國共發(fā)生31.2萬起火災，其中由電氣原因引發(fā)的火災9.5萬起，占火災總數(shù)的30.40%[1]。由此可見，準確地識別故障電弧對用電安全、預防電氣火災具有非常重要的意義。

現(xiàn)階段故障電弧識別方法有：基于小波變換能量與神經(jīng)網(wǎng)絡(luò)結(jié)合的方法，對采集的信號進行小波分解，將各層細節(jié)信號能量的平均值和標準差輸入BP神經(jīng)網(wǎng)絡(luò)后構(gòu)成小波神經(jīng)網(wǎng)絡(luò)，實現(xiàn)不同負載的測試樣本識別[2]；基于支持向量機(SVM)的故障電弧識別，將采集到的電流數(shù)據(jù)訓練基于支持向量機的故障電弧識別模型，對線性負載和非線性負載回路中的故障電弧的特征識別[3]；基于遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的故障電弧識別[4]；基于電弧電磁信號的故障電弧識別方法，通過故障電弧時域信號的模最大值提取，得到故障電弧的電磁輻射信號特征值進行識別[5]；采用離散小波變換來獲得線電流波形的時頻域特征，以反映串聯(lián)型故障電弧，然后用一些實測數(shù)據(jù)對徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)進行訓練識別[6]；基于電設(shè)備放電與紫外線輻射的關(guān)系，提出了用紫外脈沖法檢測開關(guān)齒輪電弧故障的方法，這種新技術(shù)的開關(guān)柜保護系統(tǒng)通過分析電弧產(chǎn)生的紫外線，檢測出電弧故障[7]等。

本文的方法在故障電弧識別方面是第一次應(yīng)用，將實驗數(shù)據(jù)提取特征后，訓練經(jīng)過遺傳優(yōu)化EM的高斯混合模型，得到故障特征模型和正常特征模型，然后，將要識別的數(shù)據(jù)輸入到已經(jīng)得到的模型中，利用概率最大的分類方法，判斷數(shù)據(jù)的故障與否。實驗結(jié)果表明該方法在故障電弧識別方面有很大的應(yīng)用價值。

1 實驗

1.1 實驗設(shè)置

為了得到故障電弧電流數(shù)據(jù)，按照UL1699搭建串聯(lián)型故障電弧實驗裝置，實驗裝置圖如圖1所示，其原理圖如圖2所示。兩個電極材質(zhì)分別是銅和碳，通過步進電機來轉(zhuǎn)動控制兩個電極，實現(xiàn)靜止電極和固定電極的距離從相接觸到漸漸分離，從而產(chǎn)生串聯(lián)型故障電弧，在此過程中通過實時數(shù)據(jù)采集卡對電路中的實時電流進行采集。

圖2 實驗原理圖

1.2 實驗方案

實驗在實驗室環(huán)境下進行，負載采用RL負載柜，電壓為36 V，在不同電流數(shù)值下，通過調(diào)整負載柜的功率因數(shù)，獲得不同狀態(tài)下的故障電流和正常電流數(shù)據(jù)。實驗方案條件如表1所示。

表1 實驗條件

1.3 實驗結(jié)果分析

典型電流波形如圖3所示。負載回路中發(fā)生故障電弧時，電流幅值發(fā)生了明顯的波動，并出現(xiàn)了短暫的“零休”現(xiàn)象；在不同實驗條件和故障電弧的穩(wěn)定性下，故障電流圖像會有所不同，但均出現(xiàn)了波峰幅值減小、波動等故障電弧的典型特征。

由于采用的是訓練模型的方法，以及故障電弧的隨機性和隱蔽性，所以要采集大量實驗數(shù)據(jù)。實驗中共采集總數(shù)據(jù)3231組，其中正常電流848組，故障電流2383組。

圖3 典型電流波形圖

2 電流特征提取

由于發(fā)生故障電弧時電流信號夾雜了大量噪聲和諧波信號，難以準確地提取特征，所以必須對采集到的原始數(shù)據(jù)進行數(shù)據(jù)預處理。對于故障電流信號的特征提取的主要方法是小波變換和傅里葉變換。傅里葉變換忽略了信號的時域特性，而小波變換兼顧了信號的時域和頻域變化，故使用小波變換做數(shù)據(jù)預處理[8]。

利用MATLAB中的ddencmp函數(shù)指令對電流信號進行閾值降噪，再使用wdencmp函數(shù)指令對電流信號所選用的db4小波包進行3層分解，圖4為降噪前后的電流波形圖。經(jīng)過降噪處理后的故障電流信號中明顯噪聲得到了有效的消除，并且其原有特征(比如“平肩部”)得到了完好的保留。

圖4 降噪前后電流波形圖

然后將db5分解為4層信號，從分解的多層信號中提取出既包含電流信號的主要趨勢又包含電流信號中干擾噪聲信號的分量來進行特征值提取，圖5為小波4層變換細節(jié)信號。取高頻信號A4和低頻信號D4、D3、D2、D1來進行特征提取[9]。

由于電流信號的變化是不規(guī)則的，所以引入了近似熵去量化其不規(guī)則程度。圖6為近似熵計算流程圖。其中：

① 向量之間的關(guān)聯(lián)程度Ci(m,r)為

(1)

式中，h為Heaviside函數(shù)；r為相似容限。

② 矢量序列{X(n)}的平均自相關(guān)程度Φ(m,r)為

(2)

③ 近似熵Sr為

Sr=Φ(m,r)-Φ(m+1,r)

(3)

通過上述操作，將提取到的各層電流細節(jié)信號的近似熵構(gòu)造成故障電弧的特征向量，并輸入到遺傳優(yōu)化EM的高斯混合模型。

3 故障電弧識別

3.1 高斯混合模型

高斯混合模型(Gaussian Mixture Model，GMM)是一種很常用的聚類分析方法，它認為樣本的概率密度是由幾個高斯模型加權(quán)之和組成的。每個高斯模型就代表了一個類，然后把樣本數(shù)據(jù)在各個高斯模型上分別進行投影，繼而得到樣本數(shù)據(jù)在各個類別上的概率情況。得到概率后可用各種不同的模型對觀察到的事

圖5 4層小波變換示意圖

圖6 近似熵計算流程圖

物和現(xiàn)象做判決選取概率最大的類所為判決結(jié)果[10]。GMM執(zhí)行的運算過程就是訓練出幾個概率分布。GMM的定義為

(4)

式中，K為模型的個數(shù)；πk為第k個高斯的權(quán)重；P(x|k)為第k個高斯概率密度，其均值為μk,方差為σk。對此概率密度的估計就是要求出πk、μk、σk。當求出P(x)的表達式后，求和式的各項結(jié)果就分別代表樣本x屬于各個類的概率。

對于GMM，在k已知的情況下，EM(Expectation-Maximization,期望最大)算法是一種有效的估計參數(shù)方法[11]，算法流程如下：

(1) 初始化。

協(xié)方差矩陣Cj0設(shè)為單位矩陣，每個模型比例的先驗概率αj0=1/M;均值μj0設(shè)為隨機數(shù)。

(2) 估計步驟(E-step)。

令αj的后驗概率為

(5)

(3) 最大化步驟(M-step) 更新權(quán)值、均值和方差矩陣。

(6)

(7)

(8)

(4) 收斂條件。

不斷地迭代步驟(2)、步驟(3)，重復更新上面3個值，直到P(X|Ф)-P(X|Ф)′<ε,其中P(X|Φ)′為更新參數(shù)后計算的值，即前后兩次迭代得到的結(jié)果變化小于一定程度則終止迭代，通常ε=10-5。

然而，EM算法的最大缺點就是在很多情況下只能收斂到局部最優(yōu)解，所以將遺傳算法引入到了EM算法中。

3.2 遺傳優(yōu)化EM算法

遺傳算法(Genetic Algorithm，GA)是一種參考生物界自然選擇和自然遺傳機制的隨機化全局搜索方法。GA主要包含遺傳變異、適者生存、進化概念3方面內(nèi)容，其主要執(zhí)行選擇、交叉和變異3個步驟[12]。理論研究表明，在進化過程中使用保留最佳個體策略的GA能收斂于最優(yōu)解的概率是趨向于1。所以將GA中的復制變異應(yīng)用在優(yōu)化EM參數(shù)估計上。

算法的過程是將遺傳算法應(yīng)用在GMM上，采用和諧函數(shù)作為適應(yīng)度的度量，與此同時選用復制(保留沒有的高斯)和變異(將高斯執(zhí)行過程)兩種遺傳算子。

其基本思想為:訓練所有樣本，從一個高斯分布開始，然后進行高斯，直到滿足達到設(shè)定好的k值終止。

主要有兩個步驟:一是通過一種全局搜索方法，找到現(xiàn)有GMM中的某個高斯分布進行，找到新分量合適的初始參數(shù);二是用EM估計，重新定義新分量的參數(shù)并使其達到最大似然值[13-15]。參數(shù)優(yōu)化過程為：

(1) 初始化，賦值K=1；對高斯模型隨機進行一次分裂，得到兩個新的初步的高斯分布；

(2) 在兩個高斯分布中隨機選取一個中心點，比較該類內(nèi)每個點與新舊兩個中心點的歐氏距離，與距離近的形成一類，這就形成了兩個更加準確的高斯分布。樣本x歸為第j個高斯的條件是，當且僅當q(j|x)=maxiq(i|x),其中：

后驗概率為

(9)

樣本xi和xj之間的歐氏距離為

(10)

(3) 使用EM算法對新產(chǎn)生的高斯進行參數(shù)優(yōu)化；判斷是否達到了迭代終止條件P(X|Ф)-P(X|Ф)′<10-5，如沒達到返回執(zhí)行(2)，依次循環(huán)下去。此為遺傳過程；

(4) 選擇具有最佳適應(yīng)度的后代進行復制并加以保留，這樣，每次結(jié)束后K值加1。當然，新增高斯之后的模型適應(yīng)度要大于之前的結(jié)果，否則取消[16]，輸出參數(shù)。

3.3 整體算法流程

整體算法流程圖如圖7所示。

圖7 整體算法流程圖

4 實驗結(jié)果

4.1 識別過程

4.1.1 訓練過程

訓練模型的目的是在識別程序輸入運行時的電流數(shù)據(jù)(故障電流數(shù)據(jù)和正常電流數(shù)據(jù))，然后按照特征的不同形成分類，即故障模型和正常模型。訓練指的是利用收集的電流樣本數(shù)據(jù)讓計算機學習如何對電流值進行分類。

4.1.2 測試過程

測試是將要識別的數(shù)據(jù)輸入已經(jīng)訓練過的模型，然后模型按照訓練時的分類特點對所要識別的數(shù)據(jù)進行概率最大分類，按照測試的順序輸出結(jié)果：正常，輸出為norm current，故障，輸出為fault current。

4.2 識別結(jié)果及分析

識別結(jié)果如表2所示。

表2 識別結(jié)果

由表2容易看出，遺傳優(yōu)化的GMM算法對故障電弧的識別效果顯著，識別率在91%～99%之間，實現(xiàn)了算法高效的識別故障電弧的目的。采集的數(shù)據(jù)特征和訓練時采用的數(shù)據(jù)數(shù)量不同，都會影響到識別率。

① 數(shù)據(jù)特征：在穩(wěn)定和不穩(wěn)定時采集到的數(shù)據(jù)其特征相差較大，穩(wěn)定時的特征具備故障電弧的大多數(shù)特征，而不穩(wěn)定時的特征變異較大，所以會影響識別率。

② 訓練數(shù)據(jù)的數(shù)目：訓練數(shù)目越多，尤其是故障電流數(shù)據(jù)，會訓練更多電流特征，建立更完善的訓練模型，分類會更加準確，識別率會更高。

5 結(jié)束語

根據(jù)故障電弧的電流信號的特點使用小波變換和近似熵計算來提取特征向量，再把提取到的特征向量輸入遺傳優(yōu)化EM后的GMM，對故障電弧進行識別。由結(jié)果可見，此方法對于特征明顯的故障電弧信號識別率可高達99.5%，而對特征特別不明顯的故障電弧信號也能達到91.1%，實現(xiàn)了故障電弧的有效識別。隨著樣本數(shù)據(jù)的增多，訓練數(shù)目增加，識別率還會進一步提高。所以，此方法能夠?qū)收想娀∵M行有效識別。