變壓器繞組振蕩波建模及軸向移位故障分析

周利軍，周祥宇，吳振宇，林 桐，徐肖偉，張陳擎宇

（1. 西南交通大學(xué) 電氣工程學(xué)院，四川 成都611756；2. 云南電網(wǎng)有限責(zé)任公司電力科學(xué)研究院，云南 昆明650217）

0 引言

變壓器是電力系統(tǒng)的核心設(shè)備之一，特別是大型電力變壓器結(jié)構(gòu)復(fù)雜、造價昂貴，一旦發(fā)生故障就會導(dǎo)致大面積的停電，造成巨大的經(jīng)濟(jì)損失［1］。據(jù)統(tǒng)計，對于110 kV 等級以上的電力變壓器，因繞組變形而導(dǎo)致的變壓器故障比例為50%～60%［2］。

受運(yùn)輸、地震等不可抗拒因素的影響，變壓器繞組易發(fā)生機(jī)械變形。另外由于電力系統(tǒng)運(yùn)行環(huán)境復(fù)雜，繞組易發(fā)生由絕緣老化、短路故障引起的電氣故障［3-4］。繞組微小形變對變壓器前期正常運(yùn)行的影響有限，然而隨著繞組絕緣的老化，其抗短路能力降低，若再次經(jīng)受短路電流可能會產(chǎn)生嚴(yán)重變形，準(zhǔn)確、有效的變壓器繞組變形檢測方法對確保電力系統(tǒng)的安全可靠運(yùn)行至關(guān)重要［5］。短路阻抗法、頻率響應(yīng)分析FRA（Frequency Response Analysis）是目前應(yīng)用最為廣泛的繞組變形檢測方法。近年來，為了能進(jìn)一步提高現(xiàn)場繞組變形檢測的可靠性，國內(nèi)外學(xué)者針對不同的檢測方法進(jìn)行了大量的研究。劉勇等結(jié)合頻率響應(yīng)法和阻抗法，提出了掃頻阻抗法，通過在實(shí)驗(yàn)室模擬繞組故障驗(yàn)證了該方法的有效性［6］；姚陳果等提出一種從變壓器套管末屏注入脈沖信號在線檢測繞組變形的方法，建立脈沖頻率響應(yīng)模型研究脈沖頻率響應(yīng)的影響因素，取得了一定的效果［7-8］；MITCHELL S D 等建立了變壓器繞組模型，闡述了不同測試方法對參數(shù)變化靈敏性的機(jī)理，用于分析頻率響應(yīng)和繞組故障之間的對應(yīng)關(guān)系［9］。上述方法對于變壓器繞組狀態(tài)檢測有一定的效果，但是由于其輸入激勵信號等級小，現(xiàn)場環(huán)境復(fù)雜易受到噪聲干擾，需要經(jīng)驗(yàn)豐富的工程師對繞組故障進(jìn)行識別。吳振宇等提出了高壓振蕩波法，通過實(shí)驗(yàn)驗(yàn)證了振蕩波用于繞組故障診斷的可行性，為了進(jìn)一步分析故障下繞組參數(shù)和振蕩波的關(guān)聯(lián)性，需要對振蕩波進(jìn)行建模分析。國內(nèi)外學(xué)者針對變壓器繞組進(jìn)行了大量建模分析：HASHEMNIA N 等建立了集總參數(shù)電路模型，得到了軸向移位以及徑向變形故障下繞組電氣參數(shù)與頻率響應(yīng)的關(guān)系［11-12］；SHINTEMIROV A等提出了適用于高、低頻的繞組混合多導(dǎo)體傳輸線模型，分析了高頻下頻率響應(yīng)的變化規(guī)律，可用于識別繞組輕微變形［13-14］；ABUSIADA A針對V-I 軌跡繞組在線檢測方法，建立了變壓器分布式參數(shù)電路模型，仿真得到了繞組在不同故障下的V-I 軌跡圖，驗(yàn)證了該方法能夠有效識別繞組故障［15］。上述建模方法都忽略了繞組間的耦合電容及互感的影響，對于繞組檢測方法的建模精度有一定的限制。

為此，本文基于一種新型變壓器繞組故障診斷方法——振蕩波法，以一臺10 kV 的三繞組變壓器為研究對象，通過有限元建模計算獲得繞組等效電氣參數(shù)，搭建了繞組全電容、電感矩陣參數(shù)狀態(tài)空間方程，獲得振蕩波數(shù)學(xué)模型。最后通過試驗(yàn)驗(yàn)證了本文振蕩波模型的有效性，并重點(diǎn)剖析了繞組軸向移位故障下電氣參數(shù)變化和振蕩波之間的關(guān)聯(lián)性，為變壓器繞組振蕩波建模提供了理論依據(jù)，為繞組故障識別提供參考。

1 振蕩波原理分析

當(dāng)頻率大于1 kHz 時變壓器繞組可以等效為無源網(wǎng)絡(luò)，由電阻、電感、電容構(gòu)成集總參數(shù)電路模型［16-17］。振蕩波是在變壓器繞組的一端輸入高壓直流電源，經(jīng)高壓可控電力電子開關(guān)（簡稱高壓可控開關(guān)）開斷形成暫態(tài)激勵信號，在繞組另一端獲取的高壓暫態(tài)響應(yīng)信號［10］。振蕩波原理如圖1 所示，在穩(wěn)態(tài)時高壓可控開關(guān)處于斷開狀態(tài)，高壓直流電源通過限流電阻給繞組前端的高壓電容充電，變壓器繞組相當(dāng)于短路狀態(tài)。當(dāng)高壓可控開關(guān)閉合時，將形成2 個回路，即直流電源經(jīng)限流電阻和開關(guān)形成的直流回路，以及高壓電容與變壓器繞組之間形成的暫態(tài)振蕩回路。高壓電力電子開關(guān)閉合后，高壓電容相當(dāng)于激勵，對繞組形成的電阻-電感-電容回路放電，由于繞組電感及電容特性，高壓電容不能立即釋放電荷，存在一個反復(fù)充放電的過程，而又由于繞組的電阻特性，充放電時電壓幅值是不斷衰減的，從而形成振蕩波信號。變壓器繞組構(gòu)成的振蕩回路是變壓器的固有屬性，當(dāng)變壓器繞組狀態(tài)改變時，振蕩波的振蕩頻率與幅值也會發(fā)生相應(yīng)的改變，兩者之間具有一定的映射關(guān)系。根據(jù)標(biāo)準(zhǔn)IEC60071-2 可知，相同激勵下的暫態(tài)響應(yīng)（即振蕩波）與變壓器繞組自身的屬性相關(guān)，這也驗(yàn)證了振蕩波信號的特征可以用于變壓器繞組的故障診斷［18］。

圖1 振蕩波原理圖Fig.1 Principle diagram of oscillation wave

為了驗(yàn)證變壓器繞組振蕩波產(chǎn)生的可行性，本文針對一臺220 kV 的大型三相電力變壓器進(jìn)行了現(xiàn)場試驗(yàn)，見附錄中圖A1。變壓器高壓側(cè)三繞組采用Y 連接，從高壓側(cè)繞組的中性點(diǎn)輸入高壓暫態(tài)激勵，在高壓側(cè)三相相繞組的末屏測試響應(yīng)信號，測試結(jié)果見圖2。圖中，由上至下3 組波形分別對應(yīng)10、20、30 kV暫態(tài)激勵的測試結(jié)果。

在暫態(tài)激勵下變壓器繞組形成的振蕩回路能夠獲得振蕩波，且振蕩頻率是變壓器繞組自身的固有屬性，和輸入激勵的幅值無關(guān)。由于三相繞組的結(jié)構(gòu)相同，在不同激勵幅值下輸出的A、B、C 相三繞組振蕩波具有良好的重復(fù)性，因此利用振蕩波識別繞組故障有一定的可行性。

圖2 現(xiàn)場試驗(yàn)振蕩波形Fig.2 Oscillation wave obtained from field test

2 有限元建模及參數(shù)計算

第1 節(jié)分析了振蕩波的基本原理，為進(jìn)一步研究振蕩波和繞組故障之間的關(guān)聯(lián)性，需要對其建模分析。本文以一臺10 kV／50 kV·A 的三繞組變壓器模型為研究對象，繞組結(jié)構(gòu)是分裂式繞組，由串聯(lián)繞組1（S1）、公共繞組（C）、串聯(lián)繞組2（S2）構(gòu)成，其主要參數(shù)見附錄中表A1。變壓器繞組可以等效為附錄中圖A2 所示的集總參數(shù)電路模型。本文基于圖A2，以單個線餅為單元計算變壓器繞組的電氣參數(shù)，每一單元由電阻Rh和電感Lh以及縱向等值電容Ch和電導(dǎo)Gh構(gòu)成，此外三繞組電路單元之間由耦合電容Ct、等效電導(dǎo)Gt構(gòu)成，并且考慮S1、C、S2對地電容Ck和電導(dǎo)Gk，以及每個單元之間的互感M。當(dāng)變壓器繞組發(fā)生故障時，等效電路模型中的電氣參數(shù)也會改變，從而影響振蕩波。

變壓器繞組電氣參數(shù)計算方法有公式法和有限元數(shù)值分析法。但是由于實(shí)際的變壓器內(nèi)部磁場分布不均且存在部分漏磁場［3］，采用公式法計算所得的電路模型參數(shù)精確度不高且計算量大。本文選擇Ansoft Maxwell 有限元分析軟件（下文簡稱Maxwell軟件）對變壓器繞組進(jìn)行等比例建模，從而計算電氣參數(shù)，所建模型見附錄中圖A3。

2.1 電阻參數(shù)的計算

在集總參數(shù)電路模型中，每個單元的電阻為繞組每個線餅所有匝串聯(lián)在一起的總阻值。考慮高頻下繞組導(dǎo)線的集膚效應(yīng)，每個基本單元的阻值為：

其中，f 為頻率；μ為導(dǎo)體的磁導(dǎo)率；σ為導(dǎo)體的電導(dǎo)率；a 為導(dǎo)體橫截面的徑向厚度；b 為導(dǎo)體橫截面的軸向高度；l為單個線餅繞組的總長度。

2.2 電感參數(shù)計算

利用Maxwell 軟件的靜磁場求解器計算電感矩陣。由于變壓器繞組電感僅與變壓器繞組材料、尺寸以及導(dǎo)磁媒介有關(guān)，而與所加激勵無關(guān)，本文設(shè)置變壓器繞組激勵源為1 A 電流源，利用靜磁場能量原理計算變壓器繞組電感參數(shù)矩陣。

其中，Wm為每一個線餅單元產(chǎn)生的磁場能量；L為單個線餅單元的導(dǎo)體電感；I 為繞組輸入電流值；H 為每一個線餅單元產(chǎn)生的磁通密度；B 為每一個線餅單元產(chǎn)生的磁感應(yīng)強(qiáng)度；V 表示變壓器模型整個區(qū)域空間。由式（2）可以推導(dǎo)出L為：

則2個電流回路系統(tǒng)之間儲存的磁場能量為：

故繞組線餅間的等效互感M為：

其中，I1、I2為2個線餅上的電流。

2.3 電容參數(shù)計算

利用Maxwell 軟件的靜電場求解器求線餅間的電容，在變壓器線餅上設(shè)置1 V 電壓源為激勵，繞組間和繞組自身的電容參數(shù)只與等效介電常數(shù)有關(guān)，而與變壓器繞組上的激勵大小無關(guān)。具體計算公式［19］如下：

其中，We為靜電場能量；D 為電通密度；E 為靜電場的電場強(qiáng)度。根據(jù)有限元計算繞組間總的電場能量進(jìn)而可以計算電容的大小。通過電容C 與電壓U 可將We表示為：

則在三維靜電場求解器中可以利用式（8）求解Ch、Ct、Ck。

圖3為繞組的對位餅間及錯位餅間耦合電容，包括繞組S1和繞組C的耦合電容Ca（ij）以及繞組S2和繞組C之間的耦合電容Cb（ij），其中i、j 分別表示第i、j個單餅繞組單元，當(dāng)i=j時表示對位餅間耦合電容，i≠j時表示錯位餅間耦合電容。在Maxwell 軟件的靜電場求解器中基于上述餅間電場能量計算公式獲得繞組餅間耦合電容矩陣。通過Maxwell 軟件計算獲得繞組正常情況下部分電氣參數(shù)值如表1所示。

圖3 繞組間耦合電容Fig.3 Coupling capacitance between windings

表1 繞組正常情況下的電氣參數(shù)值Table 1 Electrical parameter values of normal winding

3 振蕩波數(shù)學(xué)模型分析

由于變壓器繞組線餅數(shù)較多，為了獲得準(zhǔn)確的振蕩波波形，基于附錄中圖A2，考慮繞組S1、C、S2對位餅間以及錯位餅間的耦合電容和繞組間互感參數(shù)，在MATLAB 中建立變壓器繞組狀態(tài)空間方程模型，代入所求電路參數(shù)矩陣得到振蕩波波形?；诨鶢柣舴蚨山⒁怨?jié)點(diǎn)i 的電感電流Ii、節(jié)點(diǎn)電壓Ui為狀態(tài)變量的狀態(tài)方程，根據(jù)電路模型得到n 個節(jié)點(diǎn)的狀態(tài)方程為：

其中，T、T1為由-1、0、-1 組成的系數(shù)矩陣，兩者關(guān)系可以表示為TT=-T1；C、L 分別為基于有限元計算所得的全電容、電感參數(shù)矩陣［16］，G 為電導(dǎo)矩陣，G =2πfC tan δ，tan δ為導(dǎo)體絕緣介質(zhì)損耗角的正切；R 為電阻矩陣；U、I 分別為包含所有節(jié)點(diǎn)電壓、電流的列向量。當(dāng)激勵從某一節(jié)點(diǎn)輸入Ui時，為了計算集總參數(shù)模型輸出響應(yīng)將其單獨(dú)提取出來，然后對式（9）進(jìn)行變換后得到狀態(tài)方程如式（10）所示。

其中，T2為T 去掉第i 行后的系數(shù)矩陣；U'為對U 求導(dǎo)后的列向量；C1、G1分別為C、G去掉第i行和第i列后的參數(shù)矩陣；Ci為C 去掉第i行之后的第i列向量；Ui為只包含Ui的列向量；Gi為G去掉第i行之后的第i列向量；P為T去掉第i行后的第i列向量；T3為T1去掉第i 列向量后的系數(shù)矩陣。對式（10）進(jìn)行變換求解得到：

根據(jù)式（11）可以化簡得到：

其中，U'i為Ui的導(dǎo)數(shù)列向量。

式（12）是非齊次矩陣常微分方程，求解時域下的矩陣微分方程得到集總參數(shù)電路中每個節(jié)點(diǎn)的電壓列向量函數(shù)如式（14）所示。

其中，τ為積分變量；e?At0F(t0)為初始條件，利用振蕩波對變壓器繞組進(jìn)行離線測試時初始條件函數(shù)值為0，進(jìn)行在線測試時初始條件函數(shù)值不為0。針對三繞組變壓器以單個線餅為單元在MATLAB 中搭建狀態(tài)空間方程，變壓器的繞組S1、C、S2通過雙餅首尾連接構(gòu)成單餅，三繞組實(shí)際各有16 餅，故變壓器繞組總餅數(shù)為144，在集總參數(shù)電路中將每個線餅作為一個電路節(jié)點(diǎn)，建立144 階狀態(tài)空間方程，對其進(jìn)行解析獲得振蕩波仿真信號，通過實(shí)驗(yàn)驗(yàn)證振蕩波模型的正確性。實(shí)驗(yàn)采用10 V 等級的可控暫態(tài)激勵源，輸出信號為上升沿和下降沿1 μs，由現(xiàn)場測試獲得繞組振蕩波信號，實(shí)驗(yàn)示意圖見附錄中圖A4。

仿真與實(shí)測獲得的振蕩波時域信號見圖4，圖中的極值點(diǎn)表示不同時間節(jié)點(diǎn)下振蕩波的波峰和波谷。從圖4 中選取4 個近似重合的極值點(diǎn)1—4，將其相關(guān)信息列在表2 中。由表可見，仿真與實(shí)測結(jié)果的幅值最大相對誤差為7.6%，表明時域下仿真和實(shí)測的振蕩波吻合度較高。

圖4 振蕩波的仿真與實(shí)測結(jié)果Fig.4 Simulative and measured results of oscillation wave

表2 時域下實(shí)測和仿真振蕩波的極值點(diǎn)對比Table 2 Comparison of extreme points between measured and simulative oscillation waves in time domain

為了進(jìn)一步分析仿真和實(shí)測振蕩波的吻合程度，從頻域的角度分析仿真和實(shí)測振蕩波頻譜特性的差異性，本文通過傅里葉變換將振蕩波時域信號轉(zhuǎn)換至頻域獲得振蕩波的頻譜特性如圖5 所示。由圖可見，仿真和實(shí)測振蕩波的頻率都主要集中在600 kHz 左右；在300～500 kHz 范圍內(nèi)，仿真和實(shí)測振蕩波的頻譜幅值均較小，能量分布有一定的差異性。從圖5中選取3個相似的頻譜特性曲線極值點(diǎn)，相關(guān)信息如表3 所示。由表可見，在頻譜幅值最大的頻率點(diǎn)600 kHz 處，仿真和實(shí)測振蕩波之間的頻率相對誤差僅為1.2%，表明頻域下仿真和實(shí)測的振蕩波吻合度較高。

圖5 仿真和實(shí)測振蕩波的頻譜特性Fig.5 Spectrum characteristics of simulative and measured oscillation waves

表3 頻域下實(shí)測和仿真振蕩波的極值點(diǎn)對比Table 3 Comparison of extreme point between measured and simulative oscillation waves in frequency domain

仿真和實(shí)測振蕩波之間存在一定的偏差，這主要是由于建立的變壓器繞組模型進(jìn)行了一定程度的簡化和等效，有限元仿真獲得的繞組電氣參數(shù)和實(shí)際值有一定的偏差。實(shí)際的變壓器單餅繞組是由雙盤線圈首尾相連連續(xù)繞制構(gòu)成的，且每盤繞組由多匝線圈繞制而成，單匝線圈外側(cè)涂有絕緣，而在繞組三維有限元建模中，為了減小計算量將單餅繞組等效為一個整體，不同繞組間設(shè)置的等效介電常數(shù)和實(shí)際值有一定的差異，因此利用有限元計算的繞組餅間耦合電容和實(shí)際值會有一定的偏差；在計算縱向等值電容時，實(shí)際變壓器繞組S1、C、S2餅間的油道高度會有所不同，為了簡化計算，選取餅間油道高度的平均值作為等效油道高度，因此縱向等值電容的計算值會和實(shí)際值有一定的偏差。在計算對地電容時，繞組對地電容包括繞組S1對油箱壁的電容、繞組S2對鐵芯的電容，而頂部和底部位置的線餅對地電容同時包含線餅和上側(cè)油箱壁以及鐵軛之間的電容，此時電容效應(yīng)較為復(fù)雜，因此利用同軸圓柱電容計算公式進(jìn)行簡化計算時，所得電容也有一定的誤差。由圖5 可看出，振蕩波的振蕩頻率主要分布在600 kHz 左右，在高頻下受鐵芯渦流效應(yīng)，以及導(dǎo)體渦流效應(yīng)和肌膚效應(yīng)的影響，繞組電感參數(shù)將發(fā)生頻變［20］，因此仿真電感參數(shù)和實(shí)際值也有一定的偏差。電容和電感參數(shù)的計算誤差導(dǎo)致仿真振蕩波頻率和實(shí)測結(jié)果有一定的偏差。而由于400 kHz 下，仿真和實(shí)測振蕩波的頻譜能量分布有一定差異，時域下由不同頻率諧波疊加構(gòu)成的振蕩波會有一定差異，仿真振蕩波會由于和實(shí)測振蕩波相異的低頻諧波疊加而多出現(xiàn)1個波峰。

4 軸向移位故障下的振蕩波仿真分析

4.1 電氣參數(shù)變化規(guī)律分析

變壓器繞組處于時變的磁場中，當(dāng)繞組發(fā)生短路故障時，其瞬態(tài)電流較大，所受到的軸向電磁力會從幾磅增加到噸級別，可能造成繞組的整體軸向移位［21］。為探究軸向移位故障下繞組電氣參數(shù)的變化規(guī)律，本文在Maxwell軟件中針對繞組S1設(shè)置了1%、3%、5%軸向移位故障，計算不同故障程度下繞組電氣參數(shù)矩陣，為軸向移位故障下繞組振蕩波的變化規(guī)律提供理論依據(jù)。圖6 為變壓器繞組軸向移位故障示意圖，軸向移位百分比Has的計算公式如式（15）所示。

其中，Δh為軸向移位的高度；h為繞組的總高度。繞組S1的總高度約為460 mm，1%軸向移位故障下，繞組軸向移動高度約為4.6 mm。

圖6 軸向移位故障示意圖Fig.6 Schematic diagram of axial shift fault

不同程度的軸向移位故障下，繞組電氣參數(shù)變化規(guī)律如圖7所示。圖中，S1-1表示繞組S1的第1餅單元，其他依此類推。由圖7 可知，當(dāng)繞組發(fā)生軸向移位故障時，繞組的電容發(fā)生了明顯的改變，而電感參數(shù)變化不大，這是由于發(fā)生軸向移位故障時，繞組間的距離發(fā)生了改變從而對電容產(chǎn)生較大的影響，而電感與繞組自身材料屬性有關(guān)，距離變化對于電感變化的影響較小。

圖7 不同程度的軸向移位故障下，繞組電氣參數(shù)的變化Fig.7 Change of electrical parameters of windings under axial shift fault of different degrees

由圖7（a）可知，軸向移位故障下，繞組S1的頂部線餅對地電容呈現(xiàn)增大趨勢，底部線餅對地電容逐漸減小，這是由于軸向移位故障下，頂部繞組對油箱的間距減小所以頂部線餅對地電容增大，而底部繞組對油箱的間距增大所以底部線餅對地電容相應(yīng)減小。中部繞組對于油箱的距離基本不變，所以中部繞組對地電容基本不變。對比圖7（b）可知，隨著故障程度的加深，繞組間耦合電容逐漸下降，且其變化幅度遠(yuǎn)大于繞組對地電容，表明繞組耦合電容的變化是振蕩波變化的主導(dǎo)因素之一。繞組中部第8、9線餅耦合電容在1%、3%軸向移位故障下的變化相似，這主要是由于繞組S1、C 的油道有一定的高度差，從而形成繞組間的錯位，在發(fā)生1%、3%軸向移位故障時，繞組耦合間距對稱變化，其間距基本保持不變，繞組間耦合電容的變化也相對較小。由圖7（c）可知，當(dāng)繞組S1發(fā)生軸向移位故障時，繞組S1、C、S2的電感均隨著軸向移位程度的加深而減小，但是其變化幅度較小，其中繞組C的電感在5%軸向移位故障下的減小幅度相對較大，但是其電感值與正常時的電感值也非常接近，兩者比值為0.975。不同的繞組以及不同位置的線餅電感變化幅度有一定的差異，這主要是由于發(fā)生軸向移位故障時繞組磁場對稱性發(fā)生了改變，不同位置的磁場強(qiáng)度在軸向移位故障下的變化不一致，因此其電感減小幅度也有輕微的差異［11］。

4.2 振蕩波變化規(guī)律分析

將軸向移位故障下的電容、電感參數(shù)矩陣代入本文所建立的振蕩波數(shù)學(xué)模型，可以得到軸向移位故障下的仿真振蕩波波形如圖8 所示。由圖可見，隨著故障程度的加深，波形整體向上偏移；波峰、波谷的變化趨勢相似，均隨著故障程度的加深向右上方偏移。當(dāng)繞組發(fā)生軸向移位故障時，繞組間耦合電容隨著故障程度的加深逐漸減小，而振蕩波的振蕩頻率主要分布在600 kHz 附近，此時電容減小將導(dǎo)致繞組容抗變大，因此輸出的振蕩波幅值是逐漸增大的。

圖8 不同程度軸向移位故障下的振蕩波Fig.8 Oscillation waves under axial shift fault of different degrees

根據(jù)圖7 知，隨著軸向移位故障程度的加深，繞組電感逐漸減小，與繞組間耦合電容變化一致。電感、電容參數(shù)的變化導(dǎo)致振蕩波整體波形的振蕩頻率發(fā)生偏移，在時域中體現(xiàn)為振蕩波波形整體時間點(diǎn)改變，因此波峰、波谷在時間軸上的偏移方向也是一致的。為了量化軸向移位故障下振蕩波波峰、波谷的偏移特征，計算得到振蕩波波峰F1、F2、F4、F5和波谷G1、G2、G4、G5的幅值和時間偏移量，分別如表4和表5所示。表中，ΔT、ΔD分別為時間偏移量、幅值偏移量。

由表4 可見，隨著故障程度的加深，ΔT 和ΔD 均逐漸增大，驗(yàn)證了振蕩波波峰、波谷變化的趨勢相似，均是向右上方偏移。當(dāng)發(fā)生1%軸向移位故障時，波峰F4、F5的幅值偏移量相對較大，分別為6.82%和15.45%，波谷G1、G2的幅值偏移量相對較大，分別為14.48%和24.63%，因此1%軸向移位故障下，這4 個極值點(diǎn)偏移相對較大，而正常波形的最小偏移量約為7%，這表明針對較小的移位故障振蕩波也可以有效檢測繞組故障。振蕩波頻率主要集中在600 kHz 左右，但是除了頻譜幅值較高的頻帶，振蕩波也包含了幅值較小且頻率不一的諧波分量，由傅里葉變換原理可知時域下振蕩波曲線可以由600 kHz基波以及不同頻率的諧波構(gòu)成。因此1%軸向移位故障下，振蕩波的波峰F4、F5和波谷G1、G2的幅值變化相對較大，這可能是由于在不同的波峰、波谷時間段內(nèi)不同頻率的諧波分量不一致，電容以及電感參數(shù)變化在不同諧波分量下響應(yīng)變化也不一致，在時域中體現(xiàn)在波峰、波谷對于故障響應(yīng)靈敏性有一定的差異。

表4 軸向移位故障下振蕩波波峰、波谷的偏移量Table 4 Offset of peaks and valleys of oscillation wave under axial shift fault

5 結(jié)論

基于變壓器繞組振蕩波法，本文通過有限元建立了振蕩波數(shù)學(xué)模型，仿真得到了軸向移位下繞組電氣參數(shù)及振蕩波的變化規(guī)律，所得結(jié)論如下。

（1）變壓器繞組振蕩波仿真與實(shí)測下的波形相似度較高，仿真下振蕩波第1—4極值點(diǎn)與實(shí)測相比基本一致，其幅值誤差最大為7.6%，驗(yàn)證了模型的正確性及有效性。

（2）繞組軸向移位下電感參數(shù)變化較小，不同位置下繞組對地電容變化規(guī)律不同，中部繞組對地電容基本不變，頂部與底部繞組對地電容變化相反；繞組間耦合電容隨故障程度加深而減小，頂部耦合電容減小趨勢最大。

（3）基于振蕩波模型，仿真得到軸向移位故障下振蕩波曲線呈現(xiàn)整體向上偏移的趨勢，其波峰極值點(diǎn)偏移量相對較大，針對繞組輕微故障也具有較高的辨識性?；诒疚恼袷幉〝?shù)學(xué)模型，可以獲取不同故障下的振蕩波變化規(guī)律，為繞組故障診斷提供參考。

附錄見本刊網(wǎng)絡(luò)版（http：//www.epae.cn）。