★姜凱華，李東輝，孫燕南（大連交通大學研究生院，遼寧 大連 116000）

基于分形理論的牽引變壓器故障診斷

★姜凱華，李東輝，孫燕南（大連交通大學研究生院，遼寧 大連 116000）

近年來分形學迅速發(fā)展，并且在故障診斷方面得到了應用。本文運用分形理論分析了牽引變壓器發(fā)生故障時的電流采樣信號的波形特性。即根據(jù)故障牽引變壓器輸出電流波形，通過計算故障前后電流采樣信號的分形盒維數(shù)，提出了一種基于分形理論的牽引變壓器故障診斷的方法。經(jīng)過MATLAB/Simulink仿真表明，這種方法能準確識別各種故障類型，具有可行性與有效性。

分形；牽引變壓器；故障診斷；分形維數(shù)；MATLAB/Simulink Abstract: Fractal theory has been developed rapidly in recent years, and it also has been applied in fault diagnosis. The curve characteristics of the current sampling signals in the case of traction transformer fault was analyzed with fractal theory. According to the output current waveforms of fault traction transformers, the current sampling signals before and after the fault were calculated. A fault diagnosis method of traction transformer based on fractal theory is put forward. The MATLAB/ Simulink shows that the method can accurately identify various types of fault with feasibility and effectiveness.

1 引言

牽引變壓器是現(xiàn)代電氣化路供電系統(tǒng)的重要部分之一。牽引變壓器安全穩(wěn)定的工作對現(xiàn)在鐵路系統(tǒng)的運行起到至關重要的作用[1]。如果牽引變壓器發(fā)生故障，就會給鐵路系統(tǒng)帶來極大的影響。所以為了保證鐵路運輸?shù)陌踩\行，必須及時準確地判斷出牽引變壓器發(fā)生故障的故障類型，并及時做出相應的維修措施。當牽引變壓器發(fā)生故障時，會出現(xiàn)電流和電壓波動、振動、氣體、電脈沖、電輻射、磁輻射和超聲波等異常的現(xiàn)象。當前國內外對牽引變壓器故障診斷的方法有兩大類，分別是基于電氣量和非電氣量。

牽引變壓器故障診斷用電氣量檢測的方法包括局部放電法、電壓電流分析法和鐵芯多點接地檢測法等。用非電量進行故障檢測的方法包括：油量分析法、油中氣體成分含量分析法和振動在線監(jiān)測法等。但每種檢測方法都存在一些不足之處。例如，電流電壓分析法比較容易受到網(wǎng)壓波動的影響。油中氣體成分含量分析法由于色譜分析非常復雜，這樣就會對故障診斷造成較大的誤差。局部放電分析法近些年來受到廣泛關注，用脈沖電流法監(jiān)測局部放電故障時，脈沖分辨率高但信號易受干擾；基于超聲波的檢測法不易受干擾，能準確定位放電點，但靈敏度低，不能定量地判斷放電的大小[2]。

近些年來，很多專家學者為了改進上述方法的不足，對牽引變壓器故障診斷研究了大量的方法。例如神經(jīng)網(wǎng)絡、灰色理論、遺傳算法、模糊理論，專家系統(tǒng)等[3～5]，并在實際應用中取得了一定的進展。本文根據(jù)近年在非線性學科中各種理論的發(fā)展和應用，提出了一種基于分形理論[6]分析牽引變壓器故障診斷的可行性研究。

2 分形盒維數(shù)的計算

分形學是上個世紀70年代中期由哈佛大學曼德爾布羅特（Mandelbrot）教授提出的，是對沒有特征長度但是還具有一定意義下的自相似圖形和結構的總稱[7]。分形理論（FT-fractal theory）主要是對自然界中不規(guī)則事物的規(guī)律性的描述的學科。分形學可以讓人們使用新的理論去處理非線性領域中的各種問題，揭示了復雜事物背后的規(guī)律、局部和整體之間本質的內在關聯(lián)性，例如用分形學可以對不規(guī)則的幾何形體進行研究[8]。分形學的特征量是分形維數(shù)，并且分形維數(shù)在故障診斷方面得到了發(fā)展與應用。

2.1 分形盒維數(shù)

分形學研究的對象的維數(shù)是非整數(shù)的。為了更好的理解盒維數(shù)的概念我們先從一個簡單問題著手：例如要計算出一個面積為S的平面區(qū)域Z的維數(shù)的問題[9]。假設要用小正方形去覆蓋這個全部平面區(qū)域Z，所用的正方形的邊長為ε，這樣所用到邊長為ε的正方形個數(shù)為N（ε），這樣我們會得到：

兩邊取對數(shù)，得：

由此可以得到Z的維數(shù)是：

分形盒維數(shù)就是用這種覆蓋法的思想進一步求出的。把平面區(qū)域擴展集合體，把正方形擴展到超立方體。所以分形盒維數(shù)就是利用長度為ε的超立方體去覆蓋所要求的集合體X來定義的。

例如X是Rn的非空有界子集，現(xiàn)在我們用N（ε）個邊長為ε的超立方體把X覆蓋，這樣就會有：

式（4）即為盒維數(shù)的定義表達式。

2.2 分形盒維數(shù)的計算方法

設離散信號 ，Y是n維歐式空間Rn上的閉集。用盡可能細的ε網(wǎng)格劃分Rn，N(ε)是集合Y的網(wǎng)格計數(shù)。因為式（4）中的極限無法按定義求出，所以在計算時需要采用近似的方法[10]。以ε網(wǎng)格作為基準，逐步放大到kε網(wǎng)格，其中k∈Z+。這樣，令NkεN(kε)為離散空間上的集合Y的網(wǎng)格計數(shù)，則由式（2～5）和（2～6）可以計算得到。

式中：，為采樣點數(shù)。

網(wǎng)絡計數(shù) 為：

其中，。

在 圖中確定線性較好的一段為無邊度區(qū)，設無標度區(qū)起點和終點分別為 ，則：

最后，可用最小二乘法確定該直線的斜率為：

盒維數(shù) 為：

3 牽引變壓器故障類型及仿真

3.1 牽引變壓器故障

變壓器的故障可分為油箱外和油箱內兩種故障。郵箱外發(fā)生故障時，一般是引出線和套管發(fā)生相與相之間的短路或者接地短路。而郵箱內部發(fā)生故障一般是繞組的相間短路、匝間短路、接地短路或者鐵芯發(fā)生了燒損壞等[11]。

本文主要以某處牽引變電所內的牽引變壓器為對象，進行牽引變壓器的相間短路與匝間短路的仿真與研究。此變電所有兩路電源進線，一路使用時，另一路則處于熱備用的狀態(tài)，本所采用AT供電方式[12]。根據(jù)在某處牽引變電所中的調研結果，其變壓器采用的雙繞組變壓器，共兩組四臺牽引變壓器。

3.2 牽引變壓器仿真

牽引變壓器故障類型多樣，本文以相間短路和閘間短路為例，搭建了牽引變電所matlab仿真。下面對牽引變電所的仿真主要模擬了牽引電源和牽引變壓器的工作情況。圖1是對牽引變電所以及接觸網(wǎng)的仿真模型，主要包括電源模塊、牽引變壓器模塊、故障模塊、測量模塊、機車模塊、線路模塊與自耦變壓器模塊[13]。本文主要研究的是對牽引變壓器的故障類型診斷，所以下面對牽引變壓器模塊和故障模塊作簡單介紹。

圖1 牽引變電所牽引網(wǎng)模型

在本模型中電源采用Simpowersystems中三相交流電源電壓參數(shù)為220KV，頻率為50Hz；圖2是對牽引變壓器的仿真，其中牽引變壓器的容量為40MVA。其參數(shù)為一次側的電壓V1為220kV，電阻R1為2.63175Ω，電感L1為1.92577H。二次側的電壓V2為27.5kV，電阻R2為0.04112Ω，電感L2為0.03009H。激勵阻抗的電阻Rm為22224.5Ω，電感Lm為545.65H。

圖2 牽引變壓器模塊

圖3 故障模塊

圖3 是故障模塊，利用MATLAB中的斷路器（Breaker）元件以及時間（Timer）元件再串聯(lián)RLC元件行程短路故障模塊[14]。其中，利用時間元件來控制發(fā)生故障的時間，而串聯(lián)的RLC元件可以仿真發(fā)生故障時有過渡電阻的情形，可以選擇不同的電阻值來仿真不同過渡電阻下的短路故障波形[15]。設置過渡電阻為零時是理想短路模塊。本文主要討論常見的牽引變壓器相間短路故障。

4 分形盒維數(shù)的應用

牽引變壓器故障信號波形通過MATLAB仿真模型獲取，并且通過（To File）元件導出.mat文件用以分形盒維數(shù)的計算。為了避免正常狀態(tài)下的電壓電流數(shù)據(jù)過多，影響故障點左右時刻的分形計算，所以只選擇故障點左右時刻進行仿真。仿真的開始時間設為0s，結束時間設為0.04s，其它的參數(shù)保持為系統(tǒng)默認參數(shù)。在故障模塊中設置時間為0.02s時發(fā)生故障，故障類型為三相短路。然后調整三相線路故障模塊，進行不同情況的短路故障。并使用Fraclab工具箱進行分形的計算。

4.1 Fraclab工具箱

Fraclab是一種可以對信號或者圖像進行分形分析與處理的工具箱。其功能是可以計算一維或者二維信號的各種分形量，如維度，霍爾德指數(shù)或多重分形譜等。利用Fraclab工具箱還可以對信號進行評估、檢測、去噪、造型、分割、分類和合成等。本文主要由仿真中得出正常與不同故障的數(shù)據(jù)通過工具箱計算得出每個數(shù)據(jù)的分形盒維數(shù)。

4.2 仿真結果與故障分析

4.2.1 正常波形與分形維數(shù)

利用電壓波形或者電流波形判斷牽引變壓器故障。以電流波形為例，圖4為正常情 況下牽引變壓器電

圖4 正常電壓電流波形

壓與電流波形，波形呈現(xiàn)出規(guī)律的周期函數(shù)變化。利用MATLAB仿真結果提取出電流波形的數(shù)據(jù)并生成.mat文件以便用于Fraclab工具箱的使用。圖5為利用工具箱計算的分形盒維數(shù)為2.05。

圖5 正常電壓波形分形數(shù)

4.2.2 故障波形與分形維數(shù)

圖6～圖9分別為ABC三相短路、A相接地短路、AB兩相短路的波形圖與匝間短路故障波形圖，其中匝間短路故障使牽引變壓器由27.5kV變?yōu)?0kV。

圖6 ABC三相短路波形

圖7 A相接地短路波形

圖8 AB兩相短路波形

圖9 匝間短路波形

由圖8、圖9可看出，當0.02s時發(fā)生故障。電流波形發(fā)生明顯的變化，三相短路故障波形變化最大。然后提取出變化的電壓波形，利用Fraclab工具箱進行分形盒維數(shù)的計算得到下面圖10～圖13。

圖10 ABC三相短路波形分形數(shù)

圖11 A相接地短路波形分形數(shù)

圖12 AB兩相短路波形分形數(shù)

圖13 匝間短路波形分形數(shù)

此外還可以進行多種不同的故障設定，并對不同故障數(shù)據(jù)逐一進行盒維數(shù)計算，其部分數(shù)據(jù)的計算結果如表1所示。

表1 部分故障分形盒維數(shù)

由表1故障波形分形盒維數(shù)計算結果可知，正常狀態(tài)下的分形盒維數(shù)為2.05，當牽引變壓器發(fā)生ABC三相短路故障時分形盒維數(shù)為1.53。當牽引變壓器發(fā)生兩相短路故障時計算的分形盒維數(shù)分別是：AB兩相短路故障為1.20，BC兩相短路故障為1.21，AC兩相短路故障為1.19。當牽引變壓器發(fā)生單相接地短路故障時計算的分形盒維數(shù)分別是：A相接地短路故障時為1.13，B相接地短路故障時為1.14，C相接地短路故障時為1.13。當發(fā)生匝間短路故障時，牽引變壓器二次側由27.5kV下降到了10kV，分形盒維數(shù)為2.01。

由此可得出不同的故障分形維數(shù)有明顯的變化，可以明顯地區(qū)分出故障類別，也由此說明了分形學監(jiān)測牽引變壓器是否發(fā)生故障并且能判明故障類型，是可行的。而且分形學運算簡單，只需要分形維數(shù)就可以辨別故障類型，比目前常規(guī)的方法，簡潔明了。

5 結論

分形學具有自然界及人類社會中大部分信號普遍具備的一種特性，其基本理論都是構建在自相似性與無標度性基礎上的。本文主要是用MATLAB/Simulink軟件三相牽引變壓器進行仿真，并對三相牽引變壓器進行各種故障的設置。進而得到三相牽引變壓器正常運行狀態(tài)下與發(fā)生各種故障下的電壓與電流波形數(shù)據(jù)。使用Fraclab工具箱計算正常狀態(tài)和故障狀態(tài)下變壓器的分形盒維數(shù)，各類變壓器故障類型通過分形盒維數(shù)可以明確區(qū)分，可以運用到變壓器故障診斷方面，仿真結果證明了該方法的有效性和準確性。

[1] 高松, 劉志剛, 徐建芳, 何士玉, 高仕斌. 基于模型診斷和專家系統(tǒng)的牽引變壓器故障診斷研究[J]. 鐵道學報, 2013, ( 7 ) ： 42 - 49.

[2] 朱佼佼, 陳特放, 付強. 車載牽引變壓器智能故障診斷技術新研究[J]. 計算機工程與應用, 2012, ( 23 ) ： 27 - 35.

[3] 楊啟平, 薛五德, 蘭之達. 變壓器故障人工智能診斷系統(tǒng)的研究[J]. 高電壓技術, 2003, ( 3 ) ： 10 - 60.

[4] 汪曉明, 何萍, 吳花, 陳振剛, 歐陽瑾, 李彥明. CP組合神經(jīng)網(wǎng)絡在基于DGA的變壓器絕緣故障診斷中的應用[J]. 高壓電器, 2008,( 6 ) ：543 - 547.

[5] 楊兵, 丁輝, 羅為民, 張小飛. 基于知識庫的變壓器故障診斷專家系統(tǒng)[J]. 中國電機工程學報, 2002, ( 10 ) ： 122 - 125.

[6] 孫雅明，王俊豐. 基于分形理論的輸電線路故障類型識別新方法[J]. 電力系統(tǒng)自動化,2005,12 ： 23 - 28.

[7] Zhang X, Wang N, Lan G, et al. Tangential Damping and its Dissipation Factor Models of Joint Interfaces Based on Fractal Theory With Simulations[J]. Journal of Tribology, 2013, 136 ( 1 ) ： 249 - 256.

[8] 王芳, 李井永. 淺談分形理論及其應用[J]. 黑龍江科技信息, 2012, ( 7 ) ： 70.

[9] Chaudhuri B B, Sarkar N. Texture segmentation using fractal dimension [J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1995, 17 (1 ) ： 72 - 77.

[10] 郝研. 分形維數(shù)特性分析及故障診斷分形方法研究[D]. 天津大學, 2012.

[11] 張保會, 尹項根. 電力系統(tǒng)繼電器保護[M]. 北京：中國電力出版社, 2005.

[12] 李宏淵. 電力機車牽引變壓器故障診斷的技術研究[D]. 西南交通大學, 2009.

[13] 楊嶸春. AT牽引網(wǎng)吸上電流比測距精度改進的研究[D]. 西南交通大學, 2009.

[14] 張蕊, 劉新正. 平衡牽引變壓器的MATLAB建模仿真[J]. 變壓器, 2006, ( 11 ) ： 20 - 25.

[15] Yang Z, Yang H, Jiang Z, et al. A new method for calculation of flocculation kinetics combining Smoluchowski model with fractal theory [J]. Colloids & Surfaces A Physicochemical & Engineering Aspects, 2013, 423 ( 423 ) ： 11 - 19.

Fault Diagnosis for Traction Transformer Based on Fractal Theory

Fractal; Traction transformer; Fault diagnosis; Fractal dimension; MATLAB/Simulink

姜凱華（1989-），男，山東濟寧人，在讀碩士，現(xiàn)就讀于大連交通大學電氣工程及自動化專業(yè)，研究方向為電力電子技術在電力系統(tǒng)中的應用。

B

1003-0492(2017)04-0092-05

TP206 .3

李東輝（1970-），男，陜西西安人，副教授，博士，現(xiàn)就職于大連交通大學，研究方向為電力系統(tǒng)分析與控制。

孫燕楠（1979-），女，山東濱州人，博士，現(xiàn)就職于大連交通大學，研究方向為牽引變流與智能控制。