基于傳遞函數(shù)分析的高速鐵路自耦變壓器繞組軸向移位故障診斷研究

高仕斌， 嚴(yán)靜荷， 周利軍， 江俊飛， 李 威

(西南交通大學(xué) 電氣工程學(xué)院， 四川 成都 611756)

AT(Autotransformer)供電方式是目前高速鐵路的主要供電方式，其擁有供電距離長、對沿線通信線路干擾小等特點，自耦變壓器是AT供電系統(tǒng)最主要的設(shè)備之一[1]。由于高速鐵路一般采用高架橋設(shè)計，接觸網(wǎng)在高架橋上受雷擊、異物接觸等影響，會引起短路跳閘事故，而短路電流沖擊是導(dǎo)致變壓器繞組故障的主要原因。其中，軸向移位是常見的繞組故障之一，進而引起繞組間漏磁的增加及磁通分布不對稱，降低了現(xiàn)場變壓器運行的可靠性，從而影響牽引供電系統(tǒng)的安全運行[2-3]。

頻率響應(yīng)分析(Frequency Response Analysis，F(xiàn)RA)法所需的測量儀器簡便，測試靈敏度高，是目前現(xiàn)場測試變壓器繞組是否存在故障的主要方法[4-5]。頻率響應(yīng)法是在繞組一端施加頻率連續(xù)變化的正弦信號，在另一端測量信號并計算得到頻率響應(yīng)曲線[6]，目前利用FRA法進行繞組變形診斷時，都是通過故障前后或相間繞組的幅頻曲線對比得出結(jié)論[7]，未考慮系統(tǒng)的相頻特性。由FRA法的測試原理可知，其本質(zhì)是利用頻響曲線表征系統(tǒng)傳遞函數(shù)的變化，近年來的研究表明，傳遞函數(shù)的相位特性也可以反映繞組狀態(tài)[8]，為此國內(nèi)外學(xué)者對通過頻率響應(yīng)曲線幅值、相位特征獲取更為豐富的診斷信息進行了大量的研究。昆士蘭大學(xué)的學(xué)者通過FRA曲線中的幅值和相位值將FRA曲線轉(zhuǎn)換成奈奎斯特圖，并通過仿真和試驗計算了不同情況下奈奎斯特曲線的變化規(guī)律，研究故障與曲線特征的關(guān)系[9]；科廷大學(xué)的學(xué)者基于FRA曲線幅值和相位值繪制了極坐標(biāo)圖，并研究利用圖像分布特征進行故障診斷[10-11]。上述工作豐富了頻率響應(yīng)曲線的特征信息，促進了繞組變形診斷評估技術(shù)的發(fā)展，但奈奎斯特圖、極坐標(biāo)圖只是通過頻率響應(yīng)曲線轉(zhuǎn)換坐標(biāo)而得到，本質(zhì)上仍是對頻率響應(yīng)曲線進行直接分析，且分析指標(biāo)仍為偏移量特征，較為單一。由于缺乏一些更深層次的分析方法，基于FRA法的繞組變形診斷方法研究進展緩慢。近年來，通過變壓器繞組頻響數(shù)據(jù)辨識傳遞函數(shù)的方法受到國內(nèi)外學(xué)者的關(guān)注，傳遞函數(shù)能夠直接地反映繞組等效參數(shù)電路與物理結(jié)構(gòu)之間的關(guān)聯(lián)性，其系數(shù)和零極點能映射繞組的機械狀態(tài)。分析傳遞函數(shù)的變化可以判斷繞組變形，為頻率響應(yīng)分析法增添更深層次、豐富的評價信息。

圖1 自耦變壓器輻向分裂繞組

在頻響曲線擬合方法上，上海交大的學(xué)者最先通過普通多項式進行最小二乘辨識，但隨著曲線擬合階次和頻率點數(shù)的增加，該方法可能出現(xiàn)病態(tài)方程組，降低了辨識精度[12]；武漢大學(xué)的學(xué)者將普通多項式替換成正交多項式，從而防止了病態(tài)方程組的出現(xiàn)，增加了辨識的精度[13]。但應(yīng)用在高階系統(tǒng)上時，上述兩類方法的辨識效果仍舊不夠理想，而且目前大多數(shù)基于傳遞函數(shù)的診斷均是針對雙繞組變壓器，高低壓繞組間是相互獨立的。而高鐵自耦變壓器繞組采用了輻向分裂式結(jié)構(gòu)(見圖1)，其中O、A、B分別為繞組接線端口，O端口與軌道相連，A端口與接觸網(wǎng)連接，B端口與饋線連接，各分裂繞組(串聯(lián)繞組1、2)在變壓器的內(nèi)部連接，無法在外側(cè)直接測得各個分裂繞組的FRA曲線，這使得內(nèi)部傳遞函數(shù)更為復(fù)雜，增大了傳遞函數(shù)辨識的難度。而矢量匹配法擁有較高的計算效率和計算精度，在電力系統(tǒng)頻率特性擬合中應(yīng)用較為廣泛[14]，該方法將普通多項式代替為有理基函數(shù)進行擬合，可以顯著地提高擬合的精度，為本文自耦變壓器繞組傳遞函數(shù)的高精度辨識提供了較好的方案，也為基于傳遞函數(shù)進行繞組狀態(tài)的診斷提供了較為可靠的基礎(chǔ)。

本文以一臺實際自耦變壓器為研究對象，基于驗證的分布式集總參數(shù)電路模型，計算得到了正常和軸向移位故障下的繞組頻率響應(yīng)幅頻和相頻曲線，并利用快速松弛矢量匹配法直接擬合計算得到高精度的傳遞函數(shù)多項式，通過計算得到繞組傳遞函數(shù)系數(shù)序列的Kendall Rank相關(guān)系數(shù)，提出了自耦變壓器分裂式繞組軸向移位故障診斷方法，最后通過實際案例驗證了所提方法的有效性，為現(xiàn)場自耦變壓器的繞組狀態(tài)分析提供了參考。

1 變壓器模型

本文以一臺高速鐵路2×27.5 kV的單相自耦變壓器為研究對象，由于牽引供電系統(tǒng)通常規(guī)定自耦變壓器27.5 kV側(cè)的阻抗為0.45 Ω，該值較常規(guī)電力變壓器標(biāo)準(zhǔn)阻抗低，因此自耦變壓器一般采用輻向繞組分裂結(jié)構(gòu)以滿足阻抗要求[15]。

圖2 自耦變壓器結(jié)構(gòu)及頻率響應(yīng)測試

自耦變壓器中間繞組為公共繞組，串聯(lián)繞組由2個分裂繞組構(gòu)成，分別安置在公共繞組兩側(cè)，每個繞組均由68個普通連續(xù)式線餅組成。在現(xiàn)場開展頻響測試時，本文采取了圖2中的接線方式，該方法從自耦變壓器繞組A端口注入掃頻正弦信號Uin，在B端口測量電壓信號Uout，Rm為測量電阻，通過計算式( 1 )得到所有繞組的傳遞函數(shù)H(f)[15]，該測試方式能直接測得3個繞組的頻率響應(yīng)曲線，反映整個系統(tǒng)的傳遞函數(shù)特性。

( 1 )

變壓器在進行FRA測試時，繞組可以等效成圖3中的集總參數(shù)電路模型[16-18]，模型以繞組雙餅為電路單元，分別考慮了雙餅的縱向等值電容(Cs)、繞組對地電容(Cg1、Cg2)、各繞組間的耦合電容(Ck1、Ck2)、電容的并聯(lián)電導(dǎo)(Gs、Gg1、Gg2、Gk1、Gk2)、電阻(R)、自感(L)、互感(M)及各個錯位線餅間的電容(Cm)。

圖3 變壓器集總參數(shù)電路模型

為有效計算電路模型中的參數(shù)，基于有限元分析軟件構(gòu)建了自耦變壓器三維有限元仿真模型,見圖4[19]。其中電容計算結(jié)果見表1，Ck1、Ck2分別為公共繞組與串聯(lián)繞組1、串聯(lián)繞組2之間同一高度線餅之間的電容；Cg為繞組對地電容，因此通過有限元模型計算得到矩陣。

表1 電容計算結(jié)果

繞組電感隨頻率的增加而降低，這是由于隨激勵頻率增加鐵心渦流效應(yīng)更為明顯，同時也進一步降低鐵心導(dǎo)磁能力[20]，本文通過與測試頻率響應(yīng)曲線的對比，通過主要諧振點處匹配，得到鐵心磁導(dǎo)率隨頻率變化的趨勢，并計算得到電感矩陣隨頻率變化的系數(shù)，見圖5。

圖4 自耦變壓器三維有限元模型

圖5 電感矩陣隨頻率變化系數(shù)

由于本文自耦變壓器采用半硬自粘換位導(dǎo)線，可降低徑向變形的概率，但繞組軸向移位仍然難以避免，為此針對串聯(lián)繞組2、公共繞組、串聯(lián)繞組1(根據(jù)圖1，以各繞組對應(yīng)的端口對繞組進行標(biāo)識，分別簡稱為O、A、B繞組)上分別發(fā)生1%、3%、5%軸向移位故障情況下的FRA曲線進行研究，軸向移位的故障程度定義為

( 2 )

式中：Δh為繞組移位距離；h為繞組整體高度。

根據(jù)圖3中的集總參數(shù)電路模型，在數(shù)值計算軟件中建立狀態(tài)空間模型計算頻率響應(yīng)，該模型的可行性已在文獻[15, 21]中通過了驗證。通過有限元模型，計算得到了各繞組發(fā)生移位時的等效電路參數(shù)，并代入狀態(tài)空間模型中計算O、A、B繞組移位時的幅頻和相頻曲線圖，見圖6～圖8。由上述幾個圖可以看到不同的繞組發(fā)生移位時，幅頻和相頻曲線均發(fā)生了顯著的偏移，若能將幅頻和相頻曲線特征結(jié)合，診斷結(jié)果將更為準(zhǔn)確。

圖6 O繞組軸向移位時頻率響應(yīng)曲線

圖7 A繞組軸向移位時頻率響應(yīng)曲線

圖8 B繞組軸向移位時頻率響應(yīng)曲線

2 矢量匹配與傳遞函數(shù)計算

為了能將繞組頻率響應(yīng)幅頻和相頻曲線相結(jié)合，本文擬將兩者結(jié)合用于推導(dǎo)出系統(tǒng)傳遞函數(shù)，并通過比較傳遞函數(shù)變化特征來反映繞組的變化特征。

2.1 矢量匹配

為進一步提高傳遞函數(shù)計算的精度，本文采用快速松弛矢量匹配方法進行傳遞函數(shù)計算[22-23]，其中頻域傳遞函數(shù)有理分式近似式一般表示為

h(s)≈

( 3 )

式中：a、b為系數(shù)，均為實數(shù)；h(s)為頻域傳遞函數(shù)；s為拉普拉斯算子；n為擬合階數(shù)。

擬合的過程分為兩步：

(1) 極點定位

設(shè)置初始極點為p0,p1,…,pn，并引入一個極點替換函數(shù)σ(s)，使之滿足如下表達式

( 4 )

式中：pi為傳遞函數(shù)的極點；ri為傳遞函數(shù)的留數(shù)；ui為σ(s)的余數(shù)向量；d為實數(shù)。其中σ(s)在高頻率下趨向于1，將式( 4 )中的第二行左右兩邊乘以h(s)

( 5 )

因此在整個頻率段上，該問題轉(zhuǎn)換成求解如下超定線性方程組

Ax=b

( 6 )

式中：A為對角線上元素為初始極點的對角矩陣；b為值為1的列向量；x為待求解向量。

(2) 零點計算

σ(s)的零點計算等價于求解如下特征矩陣H的特征值

H=A-b·uT

( 7 )

式中：A為對角線元素為初始極點的對角矩陣；b為所有值為1的列向量;uT為σ(s)的余數(shù)向量的轉(zhuǎn)置。σ(s)的零點即是h(s)的極點，為了使該方法收斂，把新得到的極點作為初始極點繼續(xù)進行迭代。

2.2 傳遞函數(shù)計算

為了驗證本文所用快速松弛矢量匹配方法的精度，以現(xiàn)場某臺自耦變壓器實測FRA曲線為參考，進行了擬合計算，整個頻段的擬合階數(shù)為40階，以共軛極點為初始極點，極點間隔為對數(shù)間隔。計算結(jié)果見圖9，擬合結(jié)果精度較高，且擬合的均方根誤差小于10-6，可見該方法極大地提高了傳遞函數(shù)計算的精度。

圖9 頻率響應(yīng)曲線擬合

由于不同頻率下頻率響應(yīng)曲線電感、電容諧振點的分布變化所反映的信息是不同的，因此，在1 kHz～1 MHz全頻域上對網(wǎng)絡(luò)的傳遞函數(shù)進行分析，難以有效地反映繞組變形在不同頻率區(qū)間的信息量。為此本文采用了分頻段擬合的方法，在電力行業(yè)標(biāo)準(zhǔn)DL/T 911—2016[24]中頻段劃分的基礎(chǔ)上，將中頻段繼續(xù)細分為2個頻段，以更加突出中頻段的特征，其中劃分的4個頻段：1～100 kHz為頻段1， 100～300 kHz為頻段2，300～600 kHz為頻段3，600 kHz～1 MHz為頻段4。

3 傳遞函數(shù)系數(shù)對比研究

變壓器繞組實際上是一個非線性系統(tǒng)，其非線性特性主要受電感的頻變特性影響。但繞組的非線性程度實際上較小，因此可以將其近似為一個線性系統(tǒng)，通過擬合頻響曲線辨識得到傳遞函數(shù)是可行的。同時本文繞組故障仿真表明：不同故障情況下電感的頻變特性基本一致，而電感等非線性元件的頻變特性對FRA曲線偏移的影響，較繞組故障引起的電感、電容量變化引起的曲線偏移可以近似忽略不計。為了進一步增強各個頻段的線性化程度，采用分頻段擬合的方法。其中頻段1的非線性特征較其他3個頻段要強，擬合的效果不如其他頻段好，因此相應(yīng)地增加了該頻段擬合所需要的階數(shù)來提高精度，保證該方法的適用性。為了保證后續(xù)診斷的有效性，本文在擬合時確保了繞組正常和不同故障情況下各個頻段擬合的精度一致。

基于矢量匹配算法，針對圖6～圖8中的故障仿真曲線，分別在4個頻段內(nèi)進行了擬合，得到了相應(yīng)的傳遞函數(shù)多項式系數(shù)a0,a1，…,an-1和b0,b1,…,bn，并將其進行對數(shù)歸一化，為了區(qū)分，統(tǒng)一將系數(shù)歸一化后的值稱為參數(shù)。由圖6～圖8可知，頻率響應(yīng)曲線在頻段2和4內(nèi)的變化較為明顯，故重點分析頻段2和4內(nèi)不同繞組在不同故障程度下的各階分子、分母參數(shù)，見圖10～圖12，其中橫坐標(biāo)從左至右依次為ai和bi的不同階數(shù)，縱坐標(biāo)為對數(shù)歸一化后的參數(shù)值。

圖10 O繞組移位時頻段2內(nèi)各階參數(shù)變化

圖11 A繞組移位時頻段2內(nèi)各階參數(shù)變化

圖12 B繞組移位時頻段2內(nèi)各階參數(shù)變化

由圖10～圖12可以看到，隨著故障程度的增加，傳遞函數(shù)各階分子分母系數(shù)均發(fā)生了顯著偏移，且都向上進行了偏移。當(dāng)O繞組和A繞組發(fā)生移位時，參數(shù)a和參數(shù)b的值偏移較為接近，而B繞組發(fā)生移位時，參數(shù)b的變化較參數(shù)a要顯著，傳遞函數(shù)變化更為明顯，這是由于B繞組靠近油箱，移位時對地電容和互感變化更為明顯，但從幅頻曲線中難以直接發(fā)現(xiàn)。

頻段4內(nèi)的傳遞函數(shù)參數(shù)變化見圖13～圖15，與頻段2相反，當(dāng)不同繞組發(fā)生不同程度地偏移時，傳遞函數(shù)各階分子分母系數(shù)均向下發(fā)生了顯著偏移。O繞組和B繞組移位時參數(shù)a較參數(shù)b更為明顯，而公共繞組發(fā)生移位時，參數(shù)a和b變化較為接近，這是由于串聯(lián)繞組分別靠近鐵心和油箱，移位時對地電容變化更為顯著，傳遞函數(shù)變化更為明顯。不論是頻段2還是頻段4都對應(yīng)頻率響應(yīng)曲線上的諧振點偏移，說明諧振點幅值頻率的偏移與傳遞函數(shù)參數(shù)的變化呈現(xiàn)一定的相關(guān)性。

圖13 O繞組移位時頻段4內(nèi)各階參數(shù)變化

圖14 A繞組移位時頻段4內(nèi)各階參數(shù)變化

圖15 B繞組移位時頻段4內(nèi)各階參數(shù)變化

為能進一步分析繞組故障在不同頻段上的表現(xiàn)，表2中列出了不同繞組發(fā)生不同程度故障時各個頻段內(nèi)的參數(shù)變化規(guī)律，由于各個頻段內(nèi)的參數(shù)呈現(xiàn)出比較均勻的變化趨勢，因此參數(shù)平均變化率在一定程度上可表征參數(shù)總體的偏移規(guī)律。

表2 參數(shù)平均變化率對比表

由表2可知，頻段1內(nèi)的參數(shù)平均變化率要遠小于在其他頻段上的參數(shù)平均變化率，對應(yīng)于頻率響應(yīng)曲線在頻段1內(nèi)基本沒有發(fā)生什么偏移，而在其他頻段內(nèi)的偏移量要大，因此傳遞函數(shù)參數(shù)變化率的大小可以用來表征諧振點的偏移程度。

由于自耦變壓器特殊的繞組結(jié)構(gòu)，若串聯(lián)繞組發(fā)生故障，很難準(zhǔn)確判斷變壓器哪一個串聯(lián)繞組發(fā)生了故障，由圖10～圖15可以看出，隨著故障程度的加深，在頻段2、3、4內(nèi)，B繞組(靠近油箱)參數(shù)平均變化率變化的速率要明顯大于O繞組(靠近鐵心)參數(shù)平均變化率變化的速率，因此B繞組相應(yīng)參數(shù)的變化對軸向移位故障更加敏感。

為了定量地評價繞組故障的程度，引入Kendall Rank(肯德爾秩)相關(guān)系數(shù)來評價繞組故障與正常情況下傳遞函數(shù)參數(shù)序列之間的相關(guān)性。其中Kendall Rank相關(guān)系數(shù)τ的計算式為

( 8 )

式中：C為兩組序列中存在一致性的元素對數(shù)；D為兩組序列中存在不一致性的元素對數(shù)；N為每組序列的總元素個數(shù)。將4個頻段內(nèi)的參數(shù)合并到同一組序列Xa、Xb中

( 9 )

通過上述序列，計算出Kendall Rank相關(guān)系數(shù)τa、τb，繞組正常時τa、τb均為零，對其整理變形得到能表征繞組故障程度的參數(shù)τtotal，見表3，其中τtotal的計算式為

(10)

表3 特征參數(shù)τtotal對比表

由表3可知，當(dāng)繞組發(fā)生1%至5%移位時，τtotal值呈現(xiàn)顯著增加，并且可以看到B繞組的變化較其他2個繞組更為明顯。為能進一步分析特征參數(shù)τtotal與繞組故障程度之間的關(guān)系，根據(jù)標(biāo)準(zhǔn)DL/T 911—2016相關(guān)系數(shù)計算方法[24]，計算得到低、中、高頻的相關(guān)系數(shù)RLF、RMF、RHF，見表4。

表4 O、A、B 3個繞組不同故障程度下相關(guān)系數(shù)

由表3和表4對比可知，當(dāng)繞組發(fā)生1%的軸向移位故障時，通過相關(guān)系數(shù)法計算得到的結(jié)果為繞組正常，為此本文將Kendall Rank特征參數(shù)作為判斷繞組故障程度判據(jù)的同時，定義了一個判定區(qū)間。當(dāng)0.4≤τtotal≤0.6時，繞組有可能發(fā)生了輕度故障；當(dāng)0.6<τtotal≤0.75時，繞組有可能發(fā)生了明顯故障；當(dāng)τtotal>0.75時，繞組有可能發(fā)生了嚴(yán)重故障，這也進一步增加了該特征參數(shù)對故障程度較輕時的靈敏度。

4 實驗驗證

為能進一步驗證本文的方法，本文對一臺用于高速鐵路的32 MVA/55 kV自耦變壓器進行了分析，該變壓器短路沖擊試驗前后的頻率響應(yīng)曲線見圖16。由圖16可知，變壓器經(jīng)過短路沖擊試驗后，高頻段的曲線發(fā)生了一定的偏移。利用快速松弛矢量匹配算法分頻段對其擬合，擬合采用的階數(shù)為50階，擬合的均方根誤差為10-6，其中由于頻段4內(nèi)的幅頻曲線產(chǎn)生了明顯的變化，因此僅列出頻段4內(nèi)傳遞函數(shù)系數(shù)的變化，見圖17。

圖16 繞組短路試驗前后頻響曲線

圖17 頻段4內(nèi)各階參數(shù)變化

結(jié)果表明，繞組遭受短路沖擊后，頻段4內(nèi)傳遞函數(shù)參數(shù)變小，向下偏移，對應(yīng)于該頻段內(nèi)頻響曲線發(fā)生了較大的偏移。根據(jù)標(biāo)準(zhǔn)DL/T 911—2016中給出的相關(guān)系數(shù)計算方法，得出繞組在遭受短路沖擊后其繞組正常，而本文的Kendall Rank特征參數(shù)τtotal的計算結(jié)果為0.403，根據(jù)定義范圍，表明繞組有可能已經(jīng)發(fā)生了微小的移動，經(jīng)過吊罩檢查，發(fā)現(xiàn)O繞組端部墊塊有部分松動，并采取了緊固措施，驗證了本文所提方法的有效性。

5 結(jié)論

本文基于快速松弛矢量匹配算法和頻率響應(yīng)幅頻、相頻曲線對高速鐵路自耦變壓器的傳遞函數(shù)進行了高精度擬合計算，得到了正常情況和不同故障程度下的傳遞函數(shù)系數(shù)序列，通過計算其系數(shù)序列的Kendall Rank相關(guān)系數(shù)對繞組故障程度進行診斷。結(jié)果表明：

(1) 基于快速松弛矢量匹配算法提出的頻率響應(yīng)曲線分析方法，結(jié)合繞組的幅頻和相頻曲線變化特性，可實現(xiàn)高精度的系統(tǒng)傳遞函數(shù)計算，比單純利用幅頻曲線變化趨勢進行判斷更加全面，有助于增加變壓器繞組故障診斷的準(zhǔn)確性。

(2) 通過計算故障情況相對于正常情況的傳遞函數(shù)參數(shù)序列的Kendall Rank相關(guān)系數(shù)，定義了能夠表征繞組故障程度的特征參數(shù)τtotal，基于本文中的仿真和實測案例分析，表明該特征參數(shù)的大小能夠較為準(zhǔn)確地對繞組變形的程度進行衡量，并對微小形變表現(xiàn)出更高的靈敏度。

(3) 從本文的案例中可以看出，與靠近鐵心的分裂繞組O(串聯(lián)繞組2)相比，遠離鐵心的分裂繞組B(串聯(lián)繞組1)隨著故障程度加深，其參數(shù)平均變化率變化的速率要大于O繞組，因此可以通過判斷參數(shù)變化率變化的速率來區(qū)分是串聯(lián)繞組O還是串聯(lián)繞組B發(fā)生了故障。