張華元
(云南電網(wǎng)通海供電有限公司,云南 玉溪 652700)
頻率響應(yīng)分析法是一種檢測電力變壓器繞組變形故障快速、有效、無損的方式,具有較小的輸出響應(yīng)偏差及較高的故障診斷穩(wěn)定性而被廣泛應(yīng)用。根據(jù)激勵電源的性質(zhì),頻率響應(yīng)分析法分為正弦頻率響應(yīng)法(Sweep Frequency Response Analysis,SFRA)和脈沖頻率響應(yīng)法(Impulse Frequency Response Analysis,IFRA)[1]。其中,正弦頻率響應(yīng)法應(yīng)用較多,已出現(xiàn)中國電力行業(yè)標(biāo)準(zhǔn)和IEC 標(biāo)準(zhǔn)。正弦頻率響應(yīng)法直接在頻域?qū)崿F(xiàn)繞組激勵和響應(yīng)信號的測量,并繪制頻率響應(yīng)曲線,因而不易受外界干擾和系統(tǒng)穩(wěn)定性的影響,但該方法目前離線應(yīng)用較多,檢測速度一般,并且繞組狀態(tài)評判標(biāo)準(zhǔn)多采用基于相關(guān)系數(shù)評判的中國電力行業(yè)標(biāo)準(zhǔn),具有很大的局限性[2]。
目前對于脈沖頻率響應(yīng)法的應(yīng)用中,激勵信號和響應(yīng)信號的測量在時域?qū)崿F(xiàn),頻率響應(yīng)曲線的獲取,大多采用快速傅里葉變換算法[3]。然而,快速傅里葉變換 (Fast Fourier Transform,F(xiàn)FT)適用于處理平穩(wěn)信號,對暫態(tài)信號的時頻轉(zhuǎn)換容易造成頻譜泄露和柵欄現(xiàn)象等缺陷[4]。因此,有必要開展適合處理暫態(tài)信號新的時頻分析方法以獲得脈沖頻率響應(yīng)曲線。
本文提出以短時傅里葉變換算法(Short Time Fourier Transform,STFT)處理暫態(tài)信號,以期獲取正確的繞組脈沖頻率響應(yīng)圖譜,通過理論推導(dǎo)、仿真驗證和試驗驗證,該方法能夠提高脈沖頻率響應(yīng)法在線應(yīng)用與繞組狀態(tài)評估的準(zhǔn)確性。
短時傅里葉變換中,短窗寬能夠提供較好的時域解析度,長窗寬能夠提供較好的頻域解析度,根據(jù)Heisenberg 不確定性準(zhǔn)則,短時傅里葉變換的時域解析度和頻域解析度不能同時達到最優(yōu),這就限制了它的應(yīng)用[5]。但由于本文的重點是正確獲取繞組脈沖頻率響應(yīng)曲線,獲取曲線的前提是獲得信號精確頻率分布,只要能夠在保證足夠頻率分辨率的同時,考慮一定的時域分辨率,即能不失真地獲得繞組脈沖頻率響應(yīng)曲線,又能消除采用快速傅里葉變換帶來的缺陷[6]。
信號x (t)的連續(xù)短時傅里葉變換定義如式(1),其中ω 其中ω 是窗函數(shù):
信號x (n)的離散短時傅里葉變換定義如式(),實際上,離散形式是連續(xù)的一種采樣形式。
假定脈沖頻率響應(yīng)分析法中的激勵電壓信號Vin (t)和響應(yīng)電壓/電流信號Rout (t)的短時傅里葉變換分別由Vin (里葉變換和Rout (葉變換分別表示,根據(jù)邊際譜的定義,則脈沖頻率響應(yīng)曲線可由下式(3)獲得。
離散形式如式(4)所示:
為了驗證本文提出方法的可行性,建立了一個變壓器單繞組的等效電路模型,該繞組由67 餅構(gòu)成,R、L 和C 構(gòu)成的單餅等值電路模型如圖1 所示。
圖1 模擬型變壓器
本文將Matlab 軟件和Pspice 軟件相結(jié)合,利用Pspice 進行電路瞬態(tài)仿真,利用Matlab 進行數(shù)據(jù)后處理。
需要注意的是,Pspice 軟件的瞬態(tài)分析采用了變步長的方法以減小累積誤差,這在Matlab 數(shù)據(jù)后處理時無法獲取等間隔的采樣數(shù)據(jù),因此,在進行時頻轉(zhuǎn)換前,本文采用了三次樣條插值的方式,從Pspice 仿真結(jié)果獲取等間隔采樣數(shù)據(jù)。為了檢驗本文提出方法的可行性,仿真時對繞組等效電路模型進行了掃頻分析,獲取繞組的正弦頻率響應(yīng)曲線;數(shù)據(jù)處理同時采用了傳統(tǒng)的FFT變換,進行短時傅里葉變換時采用Hamming 窗,最后將三種結(jié)果進行對比,具體仿真流程如圖2所示。
圖2 仿真分析流程圖
進行脈沖注入繞組的仿真分析時,脈沖參數(shù)選擇接近實際試驗使用的10 ns 前沿、10 ns 后沿、500 ns 脈寬、500 V 幅值的方波脈沖,脈沖頻率響應(yīng)曲線由注入繞組首端的激勵電壓信號和響應(yīng)電流信號構(gòu)成。
采用本設(shè)計提出的基于短時傅里葉變換的脈沖頻率響應(yīng)曲線獲取方法處理上述數(shù)據(jù),結(jié)果如圖3-圖5 中的(a)所示,圖中還顯示了FFT 算法處理的結(jié)果。為了更直觀的展現(xiàn)包括正弦頻率響應(yīng)分析法01MHz 頻段的低頻部分頻率響應(yīng)曲線的走勢和諧振點信息,采用繞組等效電路導(dǎo)納參數(shù)作為頻率響應(yīng)曲線的縱軸,橫軸頻率取為05MHz,分析結(jié)果如圖3-圖5 中(b)所示,其中藍色、綠色和紅色曲線分別表示A、B 和C 三相繞組脈沖頻率響應(yīng)曲線。
圖3 Y 形連接中性點接地下三相試驗曲線對比圖
圖4 Y 形連接中性點不接地下三相試驗曲線對比圖
從圖3-圖5 中可以看出,經(jīng)過STFT 變換處理后的脈沖頻率響應(yīng)曲線主要諧振點的信息未丟失,特別是反映繞組各種變形故障的頻段,在保留基本諧振點信息的同時幾乎不受噪聲影響,有利于根據(jù)頻率響應(yīng)曲線判斷繞組的變形狀況。
參考中國電力行業(yè)標(biāo)準(zhǔn)DL/2004,采用相關(guān)系數(shù)數(shù)學(xué)指標(biāo)評估兩個頻率響應(yīng)曲線的相似程度,兩個轉(zhuǎn)移函數(shù)序列X (k)和Y (k),k=0,1 系數(shù)數(shù)學(xué)的標(biāo)準(zhǔn)偏差定義如下式(5)和(6)。
圖5 三角形連接下三相試驗曲線對比圖
序列的協(xié)方差定義如下式(7):
序列的歸一化協(xié)方差系數(shù)如下式(8):
相關(guān)系數(shù)定義為式():
對FFT 和STFT 變換的結(jié)果分別進行處理,依據(jù)DL/2004 頻段劃分標(biāo)準(zhǔn)在01 范圍內(nèi)計算兩相曲線之間的相關(guān)系數(shù),如下表1 所示。從表中可以看出,STFT 處理后的兩相繞組間的相關(guān)系數(shù)均大于4,明顯大于FFT 處理后的相關(guān)系數(shù)。根據(jù)試驗設(shè)置可知,大修后的變壓器三相繞組健康狀態(tài)良好,結(jié)構(gòu)一致,頻率響應(yīng)曲線應(yīng)盡可能相似。因此,本文提出的基于STFT 的變換方法優(yōu)于FFT 算法,更適合處理暫態(tài)信號,從而獲取繞組正確的脈沖頻率響應(yīng)曲線。
表1 不同數(shù)據(jù)處理方式下相間曲線數(shù)據(jù)對比(Y 形連接躍居第一 點接地)
另外,需要說明的是,經(jīng)過STFT 處理后的三相頻率響應(yīng)曲線在02 頻段差異較小,隨著頻率的增加,曲線差異明顯變大,且較FFT 處理的結(jié)果更清晰。將頻率響應(yīng)法的檢測頻段提高至10,可以檢測繞組的微小變形,結(jié)合本文測試結(jié)果,由于頻響曲線的高頻段對繞組的微小變形較為靈敏,變壓器的三相健康繞組結(jié)構(gòu)一致,相間繞組的微小差異即在頻率響應(yīng)曲線的高頻段體現(xiàn)出來。
本文提出采用短時傅里葉算法構(gòu)建變壓器繞組的脈沖頻率響應(yīng)曲線,采用Pspice 軟件建立了單繞組等效電路模型并進行暫態(tài)時域仿真,通過仿真分析可知,經(jīng)過短時傅里葉算法處理后,頻率響應(yīng)曲線清晰度高,噪聲較小,頻響曲線更接近實際情況,通過測試結(jié)果可以看出,高頻段曲線變化對繞組的微小變形較為靈敏變壓器相間繞組的微小差異即在頻率響應(yīng)曲線的高頻段體現(xiàn)出來,因此本文提出的方法可以快速準(zhǔn)確的診斷變壓器故障。
[1]郭創(chuàng)新,朱承治,張琳,等.應(yīng)用多分類多核學(xué)習(xí)支持向量機的變壓器故障診斷方法[J].中國電機工程學(xué)報.2010 (13).
[2]楊廷方,李景祿,曾祥君,等.基于多方法組合診斷模型的大型變壓器故障診斷[J].電力系統(tǒng)自動化,2009(20).
[3]顏秋容,劉欣,尹建國.基于小波理論的電力變壓器振動信號特征研究[J].高電壓技術(shù),2007 (01).
[4]汲勝昌,王世山,李清泉,等.用振動信號分析法監(jiān)測變壓器繞組狀況[J].高電壓技術(shù),2002 (04).
[5]汲勝昌,李彥明,傅晨釗.負載電流法在基于振動信號分析法監(jiān)測變壓器鐵心狀況中的應(yīng)用[J].中國電機工程學(xué)報,2003 (06).
[6]王璋奇,王孟.電力變壓器繞組軸向振動穩(wěn)定性分析[J].中國電機工程學(xué)報,2002 (07).