付光晶, 于躍, 張峰, 張士文
(上海交通大學(xué) 電子信息與電氣工程學(xué)院,上海 200240)
GIS作為氣體絕緣全封閉組合式電器設(shè)備,與普通的高壓變電站設(shè)備相比,具有體積小、空間利用率高、故障率小、運行安全性高的特點,因此目前在高壓、超高壓和特高壓領(lǐng)域被廣泛應(yīng)用。但是由于其密封式的結(jié)構(gòu)設(shè)計,也帶來了故障定位困難,檢修流程繁雜,事后斷電時間久,斷電范圍大等缺點。對于GIS頻發(fā)的分合閘故障,傳統(tǒng)的對設(shè)備進行定期定點故障檢測和修理的方法已經(jīng)不能滿足設(shè)備自身的工作要求,因此設(shè)計出一種能夠快速準確分析GIS分合閘狀態(tài)的方法便具有了重要的意義。
目前對于GIS分合閘故障檢測的方法主要是在設(shè)備分合閘動作時對激勵線圈中的驅(qū)動電流進行采樣,然后再利用相應(yīng)的算法對采樣的電流波形進行故障分析。例如華中科技大學(xué)的易慧于2006年提出了最小二乘法和小波理論可以對分合閘電流進行故障分析[1]。浙江大學(xué)的曹飛于2008年采用分散包絡(luò)統(tǒng)計法對分合閘線圈電流波形的輪廓進行綜合分析和研究[2]。江蘇財經(jīng)學(xué)院的宋錦剛于2011年提出了將分合閘信號進行小波分解后對固定頻段的信號重構(gòu),然后利用相似性原則對故障信號進行分類的方法[3]。上海交通大學(xué)的黃志于2015年利用核主元分析法對分合閘線圈電流波形特征值進行降維,再利用支持向量機對故障波形進行分類[4]。本文借鑒和參考了國內(nèi)外學(xué)者對同類課題的研究成果,針對以往算法中模型復(fù)雜,預(yù)測精度差等缺點,依據(jù)GIS分合閘經(jīng)典波形的特征點搭建了BP神經(jīng)網(wǎng)絡(luò),并利用遺傳算法和大量的樣本對神經(jīng)網(wǎng)絡(luò)進行優(yōu)化和訓(xùn)練,實現(xiàn)了GIS分合閘故障的有效分析和精準分類。
圖1顯示了GIS正常分合閘時驅(qū)動線圈上的經(jīng)典電流波形[5]285,整個分合閘過程可以被分為4個階段,其中T1~T2階段為線圈初始通電階段,此時電流按照指數(shù)形式上升,但由于鐵芯受到的機械阻力大于電磁力,鐵芯無法運動。隨著通電線圈內(nèi)的電流的不斷增加,到達T2時刻時,鐵芯的電磁力將和其受到的阻力相等,此時鐵芯處于運動的臨界狀態(tài),隨后在T2~T3階段鐵芯在電磁力的作用下開始做加速運動并產(chǎn)生反向感應(yīng)電動勢,使得線圈電流減小,直至鐵芯在T3時刻觸碰到分合閘觸發(fā)器后停止運動,此后線圈電流會在T3~T4階段再次以指數(shù)形式上升,直至T4時刻鐵芯能夠推動脫扣裝置脫扣,脫扣后,分合閘彈簧將釋放彈性勢能,帶動斷路器本體迅速動作,在T4~T5階段輔助開關(guān)將開始動作,產(chǎn)生電弧,此時電壓迅速升高,電流將迅速減小到最初值。
圖1 正常分合閘情況下線圈電流波形
從圖中可以看到,分合閘過程中線圈電流波形存在三個極值點,分別出現(xiàn)在T2,T3,T4三個時刻,其中T2時刻是鐵芯動作的臨界點,T3時刻鐵芯撞擊到分合閘觸發(fā)器,T4時刻脫扣裝置脫扣,這三組極值點的數(shù)值大小對于分析分合閘故障具有非常重要的意義。
圖2顯示了出現(xiàn)脫扣卡滯或者脫扣過軟情況時線圈電流的變化趨勢,大量試驗證明,機構(gòu)老化和潤滑過度經(jīng)常是導(dǎo)致分合閘出現(xiàn)脫扣問題的主要原因[6]。當(dāng)GIS設(shè)備出現(xiàn)老化時,脫扣裝置卡滯會延長T3~T4區(qū)間的時間,此時線圈電流達到第二個峰值的時間會加長,也會間接導(dǎo)致電流的峰值增大。相反,當(dāng)GIS設(shè)備出現(xiàn)潤滑過度時,脫扣裝置過軟會縮短T3~T4區(qū)間的時間,此時線圈電流達到第二個峰值的時間會縮短,其電流的峰值也會小于正常分合閘時的情況。
圖2 脫扣卡滯或過軟情況下線圈電流波形
圖3顯示了當(dāng)供電電壓高于或低于正常值時線圈電流的變化趨勢,大量試驗表明,當(dāng)供電電壓過高時,由于設(shè)備的等效電阻沒有變化,此時激勵電流的增長速度會快于正常情況,鐵芯產(chǎn)生的電磁力會大于正常值,從而導(dǎo)致鐵芯的臨界運動時間和脫扣裝置的脫扣時間提前到達,使得分合閘線圈電流的波形被整體抬升。相反,當(dāng)供電電壓不足的時候,激勵電流的增長速度會慢于正常情況,鐵芯產(chǎn)生的電磁力會小于正常值,從而導(dǎo)致鐵芯的臨界運動時間和脫扣裝置的脫扣時間延緩到達,使得分合閘線圈電流的波形被整體拉低。
圖3 供電電壓過高或過低情況下線圈電流波形
BP神經(jīng)網(wǎng)絡(luò)是目前應(yīng)用最為普遍的一種神經(jīng)網(wǎng)絡(luò),其網(wǎng)絡(luò)的拓撲結(jié)構(gòu)一般包含輸入層,隱藏層和輸出層,其中隱藏層又以單隱層應(yīng)用最多。BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練由輸入信號的正向計算和輸出信號的反向回饋組成,當(dāng)輸入信號經(jīng)過網(wǎng)絡(luò)計算后,其網(wǎng)絡(luò)輸出與期望輸出之間有誤差時,我們可以將誤差定義為:
(1)
式中:dk代表了輸出層第k個神經(jīng)元的期望輸出;ok代表了輸出層第k個神經(jīng)元的實際輸出。將上式展開到隱藏層,有:
(2)
式中:函數(shù)f代表了從隱藏層到輸出層的轉(zhuǎn)移函數(shù);wik代表了隱藏層第i個神經(jīng)元到輸出層第k個神經(jīng)元的網(wǎng)絡(luò)權(quán)值;yj代表了隱藏層第j個神經(jīng)元的網(wǎng)絡(luò)輸出。
BP神經(jīng)網(wǎng)絡(luò)對權(quán)值調(diào)整的原則是使輸出誤差不斷減小,因此權(quán)值調(diào)整量被定義為與輸出誤差的負梯度成正比關(guān)系,即:
(3)
(4)
式中:w和v分別代表了隱藏層到輸出層的權(quán)值和輸入層到隱藏層的權(quán)值;負號代表了權(quán)值按照負梯度下降;常數(shù)η代表了神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)速率,其一般定義為(0,1)中的一常數(shù)。對上述兩式繼續(xù)推導(dǎo),可得:
(5)
(6)
式中:對輸出層均有j=0,1,…,m;k=1,2,…,l;對隱藏層均有i=0,1,…,n;j=1,2,…m;f′(netk)代表了輸出層轉(zhuǎn)移函數(shù)對輸出層第k個神經(jīng)元的偏導(dǎo)數(shù);f′(netj)代表了隱藏層轉(zhuǎn)移函數(shù)對隱藏第j個神經(jīng)元的偏導(dǎo)數(shù)。
利用上述公式即可對網(wǎng)絡(luò)權(quán)值和神經(jīng)元閾值進行一次調(diào)整,反復(fù)利用樣本數(shù)據(jù)對神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練即可使神經(jīng)網(wǎng)絡(luò)誤差不斷收斂以致達到理想狀態(tài)。
在本設(shè)計中,根據(jù)GIS分合閘電流波形的特點,建立三層且隱藏層神經(jīng)元個數(shù)為10的神經(jīng)網(wǎng)絡(luò)。圖4展示了本神經(jīng)網(wǎng)絡(luò)系統(tǒng)的整體結(jié)構(gòu), 其中神經(jīng)網(wǎng)絡(luò)輸入層對應(yīng)分(合)閘線圈電流波形中三組極值點處的6個特征值,輸出層對應(yīng)分合閘動作后設(shè)備的5種狀態(tài)類型。設(shè)備的狀態(tài)類型和其對應(yīng)的標簽值如表1所示。
圖4 BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)
故障類型電壓過高電壓過低正常情況脫扣過軟脫扣卡滯標簽值-4-2024
本設(shè)計中輸入層到隱藏層和隱藏層到輸出層的轉(zhuǎn)移函數(shù)都設(shè)置為Sigmoid函數(shù),學(xué)習(xí)效率η設(shè)置為0.2,權(quán)值引用動量表達式ΔW(t)=ηδX+αΔW(t-1),其中α=0.2,利用神經(jīng)網(wǎng)絡(luò)對GIS分合閘故障判斷的算法流程圖如圖5所示,其一共分為以下幾個步驟:
第一步,搭建三層BP神經(jīng)網(wǎng)絡(luò),并為神經(jīng)網(wǎng)絡(luò)權(quán)值和神經(jīng)元閾值設(shè)置初始值,設(shè)置轉(zhuǎn)移函數(shù)和學(xué)習(xí)效率。
第二步,將GIS所有的分合閘線圈電流樣本輸入到神經(jīng)網(wǎng)絡(luò)中,經(jīng)過神經(jīng)網(wǎng)絡(luò)拓撲結(jié)構(gòu)計算得到輸出值,并利用式(1)計算出樣本標簽值和實際輸出之間的誤差。
第三步,利用式(2)~式(6)和第二步得到的誤差值對網(wǎng)絡(luò)的權(quán)值和神經(jīng)元閾值進行負梯度調(diào)整,并反復(fù)迭代,直至達到訓(xùn)練次數(shù)。判斷誤差值是否達到精度要求,如果沒有達到則重新進行訓(xùn)練,如果達到了則表明訓(xùn)練過程結(jié)束。
第四步,將測試樣本代入訓(xùn)練后的神經(jīng)網(wǎng)絡(luò)中,利用輸出值與設(shè)定標簽值之間誤差判斷神經(jīng)網(wǎng)絡(luò)的性能。
圖5 BP神經(jīng)網(wǎng)絡(luò)整體流程圖
對于上節(jié)利用神經(jīng)網(wǎng)絡(luò)進行故障分類的算法,雖然其結(jié)構(gòu)簡單且易于操作,但是卻具有一些天生的缺陷,例如誤差經(jīng)常僅能收斂到局部極小點而無法達到全局最優(yōu)點,這種問題經(jīng)常是由于神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和神經(jīng)元閾值設(shè)置不當(dāng)而導(dǎo)致的。針對這些問題,可以用遺傳算法來進行優(yōu)化。遺傳算法的原理是將問題的待求量轉(zhuǎn)化成二進制染色體,然后在問題求解的過程中,反復(fù)根據(jù)“適者生存”的原則對染色體進行選擇,交叉和變異,以此來交換染色體的基因,使其能夠達到最優(yōu)解。
圖6展示了如何利用遺傳算法對基于神經(jīng)網(wǎng)絡(luò)的GIS故障分類方法進行優(yōu)化,其主要算法流程分為以下幾步。
第一步:建立神經(jīng)網(wǎng)絡(luò),設(shè)置神經(jīng)網(wǎng)絡(luò)參數(shù),將GIS分合閘線圈電流樣本歸一化。
第二步:為神經(jīng)網(wǎng)絡(luò)的連接權(quán)值和神經(jīng)元閾值進行二進制編碼,并將其連接成串組成一個個體,隨機生成40組這樣的個體組成一個種群,針對種群中的每個個體對應(yīng)的神經(jīng)網(wǎng)絡(luò),分別將歸一化后的樣本輸入并計算網(wǎng)絡(luò)輸出值。
第三步:利用網(wǎng)絡(luò)輸出值計算適應(yīng)度函數(shù)值,根據(jù)適應(yīng)度函數(shù)值對樣本重新選擇,對于適應(yīng)度不高的個體進行淘汰,隨后對種群中的個體進行交叉和變異以豐富種群特性。
第四步:將進化過的離散種群中的每個個體的數(shù)值重新賦值給神經(jīng)網(wǎng)絡(luò)權(quán)值和神經(jīng)元閾值,針對新種群中的由每個個體組成的神經(jīng)網(wǎng)絡(luò),再次輸入樣本計算子代的目標函數(shù)值。
第五步:根據(jù)子代的目標函數(shù)值,將子代個體選擇性的插入到父代樣本中。
第六步:遺傳代數(shù)計數(shù)器加1,判斷遺傳代數(shù)是否達到設(shè)定數(shù)值,如果沒有達到,返回步驟三繼續(xù)進行進化,如果達到了,輸出每次進化后的誤差值。
圖6 遺傳算法優(yōu)化設(shè)計流程圖
對于GIS分合閘過程中線圈電流的一些經(jīng)典特征值,經(jīng)過對文獻[5]289的查閱將其列在表2中,從表2的數(shù)據(jù)可以發(fā)現(xiàn)其特征值的變化趨勢與第1節(jié)所描述的電流波形特性基本一致。其中(T2,A2)對應(yīng)鐵芯動作臨界點的時間和線圈電流,(T3,A3)對應(yīng)鐵芯撞擊到分合閘觸發(fā)器的時間和線圈電流,(T4,A4)對應(yīng)脫扣裝置脫扣的時間和線圈電流。
表2 GIS分合閘線圈電流經(jīng)典特征值
利用表2中的經(jīng)典特征數(shù)據(jù),針對每種分合閘情況擬合40組特征點,并將前20組用于神經(jīng)網(wǎng)絡(luò)的訓(xùn)練,后20組用于訓(xùn)練后神經(jīng)網(wǎng)絡(luò)的性能測試,表3為抽調(diào)的100組訓(xùn)練樣本中的5組,并分別將不同的分合閘電流特征值賦予相應(yīng)標簽值。
表3 神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本抽調(diào)
將100組數(shù)據(jù)點輸入到神經(jīng)網(wǎng)絡(luò)中經(jīng)過100次訓(xùn)練后,可以得到相應(yīng)的網(wǎng)絡(luò)參數(shù),輸入層和隱藏層之間的權(quán)值系數(shù)V為:
其中權(quán)值矩陣V的第i行第j列元素對應(yīng)的輸入層的第i個神經(jīng)元和隱藏層的第j個神經(jīng)元之間的權(quán)值大小。而隱藏層和輸出層之間的權(quán)值系數(shù)W為:
W=[0.407 -0.138 -0.059 -0.369 0.014 -0.501 0.079 -1.180 0.005 0.208]T
權(quán)值矩陣W的第k列元素對應(yīng)的隱藏層的第k個神經(jīng)元和輸出層神經(jīng)元之間的權(quán)值大小。此外,經(jīng)過神經(jīng)網(wǎng)絡(luò)訓(xùn)練的還有隱藏層和輸出層的閾值,經(jīng)過訓(xùn)練后其值為:
T1=[-2.352 2.112 -1.784 -0.719 1.031 -1.140 -0.289 -2.303 -2.057 -1.530]TT2=-0.235
至此,經(jīng)過100次訓(xùn)練的神經(jīng)網(wǎng)絡(luò)已經(jīng)形成,圖7展示了在100次神經(jīng)網(wǎng)絡(luò)訓(xùn)練的過程中,輸出誤差的變化情況。
圖7 BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練過程誤差變化
從圖7中可以看到,神經(jīng)網(wǎng)絡(luò)誤差值在前10次訓(xùn)練過程中收斂迅速,但是在第10次到第50次之間速度變緩,直到訓(xùn)練到53次時,神經(jīng)網(wǎng)絡(luò)的網(wǎng)絡(luò)權(quán)重和神經(jīng)元閾值進入到極小值區(qū)域,網(wǎng)絡(luò)趨于穩(wěn)定。網(wǎng)絡(luò)經(jīng)過100次訓(xùn)練后其輸出與期望值的誤差為1.824。圖8顯示了經(jīng)過50次進化后,每次進化過程中種群中適應(yīng)性最優(yōu)個體經(jīng)過100次訓(xùn)練后對應(yīng)的誤差值。
圖8 遺傳算法進化過程中誤差的變化
在擬定遺傳算法參數(shù)時,將子代和父代之間的代溝設(shè)置為0.95,這樣可以保證誤差隨著遺傳代數(shù)的變化持遞減性的變化。圖8顯示出當(dāng)種群進化到22代時,其誤差值已經(jīng)達到0.119,在隨后的的進化過程中誤差值沒有再縮小,對比單獨利用神經(jīng)網(wǎng)絡(luò)進行故障分類時,當(dāng)訓(xùn)練到53次便進入極小值區(qū)域且誤差值為1.824,其性能有了極大的提高。
針對訓(xùn)練結(jié)束后的神經(jīng)網(wǎng)絡(luò),本試驗利用剩余的100組樣本對網(wǎng)絡(luò)性能進行測試,表4抽調(diào)了100組測試樣本中的5組。
表4 神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本抽調(diào)
從表4中可以看出,單單使用神經(jīng)網(wǎng)絡(luò)進行故障分類,當(dāng)判據(jù)苛刻時,成功率只有80%,而應(yīng)用遺傳算法優(yōu)化后的算法,其成功率為100%。
此外,圖9展示了100個測試樣本在遺傳算法優(yōu)化前后輸出誤差的對比,圖中清晰地表明了經(jīng)過了遺傳算法優(yōu)化的神經(jīng)網(wǎng)絡(luò),在50次進化后,測試樣本的輸出誤差明顯減小。
圖9 引入遺傳算法前后誤差對比
本文簡介了GIS分合閘時的基本特性和典型故障類型,并提出了基于BP神經(jīng)網(wǎng)絡(luò)的和遺傳算法的故障分類方法,試驗表明,利用遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)具有結(jié)構(gòu)簡單,分類精度高的優(yōu)點,能夠很好地對GIS分合閘故障進行分類。