基于波形時域分布特征的變壓器勵磁涌流識別

王 雪 王增平

（華北電力大學電力工程系 保定 071003）

1 引言

差動保護仍是當前廣為應用的變壓器主保護，隨著大容量變壓器的投入運行，變壓器差動保護的快速性和可靠性面臨更加嚴峻的考驗，如何區(qū)分勵磁涌流和故障電流仍是變壓器差動保護的重要研究課題之一。工程中主要應用二次諧波制動判據、間斷角判據、波形對稱判據來防止勵磁涌流造成差動保護誤動[1-2]，但是由于大容量變壓器鐵心材料的改進以及系統中無功補償裝置的應用使得涌流和故障暫態(tài)電流中二次諧波含量的差別減小，二次諧波制動判據不正確動作的可能性提高了，波形對稱判據和間斷角判據都是依據勵磁涌流波形發(fā)生畸變并有間斷的特征，但鐵心極度飽和或產生對稱性涌流時，上述原理的保護都有可能誤動。近年來眾多學者提出很多新方法，總體可分為兩類:一類是繼續(xù)探究電流波形特征并利用不同的技術手段加以提取[3-6]；另一類是引入電壓量，構成諸如基于差有功[7]、磁通特性[8]或等效瞬時勵磁電感[9-10]的勵磁涌流識別方法，或者構成不依賴于差動保護的新型主保護原理，如基于等效瞬時漏感[11]或回路方程[12]的新型保護。這些新型主保護在理論上完全不受勵磁涌流的影響，具有很好的研究前景，但是由于引入電壓量必須考慮電壓互感器的暫態(tài)響應特性，且電壓互感器二次斷線時必須退出運行，此時只能依賴電流差動保護。因此研究僅依賴電流特征的勵磁涌流識別方法仍具有十分重要的意義。

當變壓器鐵心飽和時，勵磁涌流波形將出現嚴重的畸變，而故障電流基本符合正弦波的特征，二者波形的時域分布特征存在明顯的差異。本文根據勵磁涌流生成的機理，分析了勵磁涌流和故障電流波形的時域分布特征，并給出了相應的勵磁涌流識別判據。動模實驗和數字仿真結果均表明該方法具有原理簡單、便于實現、靈敏度高等特點，性能上優(yōu)于二次諧波制動原理和間斷角原理。

2 基本原理

2.1 勵磁涌流波形的時域分布特征

圖1給出了用直線近似表示的鐵心磁化曲線，圖中Φs為鐵心的飽和磁通。當鐵心磁通φ ＜Φs時勵磁電流為零，等值勵磁阻抗為無窮大，φ ＞Φs時鐵心飽和，產生很大的勵磁電流，等值勵磁阻抗大大減小(為常數)。

圖1 變壓器鐵心磁化曲線Fig.1 The magnetization curve of transformer core

以單相變壓器為例加以分析，設空載合閘時電源電壓為

式中，α為合閘初相角。

設變壓器鐵心剩磁為Φr，額定工作磁通為Φm，在簡化的鐵心磁化特性下可得勵磁涌流為[13]

式中，θ =ωt+α；θ1=arccos（cosα+kr-ks）；kr=Φr/Φm；ks=Φs/Φm。

令α=0、kr=0.8、ks=1.1，將各參數代入式（2）可得勵磁涌流波形如圖2a中實線所示，很明顯勵磁涌流波形不再滿足正弦波特征，且偏向時間軸一側。

直方圖技術是數字圖像處理中的一種常用方法，它能提供非常有用的圖像信息統計資料。對于某個灰度圖像，假設其灰度級為[0，L]，圖像總像素個數為 n，定義灰度直方圖為離散函數 h(r)=nr，其中 r表示某個灰度級(0≤r≤L)，nr表示圖像中灰度級為r的像素個數[14]。借鑒直方圖處理技術來描述勵磁涌流波形分布特征，首先以一個工頻周期內差動電流瞬時值的最大絕對值為基準對差電流進行歸一化處理，歸一化后的電流最大取值范圍為[-1，1]，將該取值區(qū)間等分為2N個子區(qū)間，各子區(qū)間標號分別記為-N, ··, -1，1, ··, N，每個區(qū)間的長度為1/N。分別統計一個工頻周期內電流瞬時值位于每個子區(qū)間的采樣點數目，假設位于第n個子區(qū)間的采樣點數目為mn，定義第n個子區(qū)間的分布系數為

式中，M為一個工頻周期的采樣點總數；n∈[-N，N]。

取N=10，縱向劃分涌流波形，如圖2a中網格線所示，計算每個區(qū)間的分布系數。將區(qū)間標號作為橫坐標，分布系數作為縱坐標，可得勵磁涌流波形分布系數直方圖，如圖2b所示?？梢娪捎诓ㄐ未嬖陂g斷部分，且偏向正半軸，因此標號為負的子區(qū)間的分布系數為零，第1區(qū)間分布系數最大，而第10區(qū)間分布系數相對要小。實際變壓器鐵心進入和退出飽和階段，勵磁涌流波形呈凹弧特性和明顯的尖頂波特性，如圖2a中虛線所示。顯然此時第1、10子區(qū)間的分布系數差異更大。

圖2 不對稱勵磁涌流及波形分布系數直方圖Fig.2 Asymmetric inrush current waveform and its distribution coefficient histogram

對于 Yd聯結三相變壓器，由于不能得到三角形側的繞組電流，實際上差動電流為兩相電流之差，此時往往有一相差動電流波形呈現對稱性。當兩相的暫態(tài)磁通相同時，可能出現正負峰值相等的對稱性勵磁涌流，如圖3a所示，此時的波形分布系數直方圖如圖3b所示。與不對稱涌流的情況相比，此時直方圖呈對稱分布，但仍具有第1區(qū)間分布系數最大，而第±10區(qū)間分布系數相對要小的特點。

圖3 對稱性勵磁涌流及波形分布系數直方圖Fig.3 Symmetric inrush current waveform and its distribution coefficient histogram

2.2 故障電流波形的時域分布特征

不計變壓器繞組電阻時，基于等效瞬時電感的變壓器模型滿足如下方程[10]:

式中，u為繞組端電壓；id為差電流；Lk為等效瞬時電感。

由式（4）可知，故障電流的波形變化特征由繞組端電壓和等效瞬時電感決定。端電壓僅在故障發(fā)生時刻發(fā)生突變，之后即保持不變。由于變壓器內部故障時鐵心不會飽和，其勵磁電感和各繞組漏感均可假定為常數，因此等效瞬時電感為常數，故障差電流也基本保持工頻正弦特性。圖4a給出了一個工頻周期正弦特性的故障電流波形，此時的波形分布系數直方圖如圖 4b所示?？梢娭狈綀D呈對稱分布，與勵磁涌流情況不同的是第±1區(qū)間分布系數最小，而第±10區(qū)間分布系數最大。

圖4 故障電流及波形分布系數直方圖Fig.4 Fault current waveform and its distribution coefficient histogram

實際上在各子區(qū)間內電流基本上是單調變化的，變化率越大，則處于該子區(qū)間內的采樣點數就越少，因此根據式（3）的定義可知各子區(qū)間的分布系數反比于該區(qū)間內電流的變化率。對于勵磁涌流，無論非對稱涌流還是對稱性涌流，第±1區(qū)間內電流波形處于間斷部分，電流變化率小，采樣點多，因而分布系數大，而在第±N區(qū)間內波形為尖頂波，電流變化率大，采樣點少，因而分布系數?。粚τ诠收想娏?，第±1區(qū)間內電流處于過零點附近，電流變化率大，采樣點少，分布系數小，在第±N區(qū)間內波形處于峰值附近，電流變化率相對要小，采樣點較多，分布系數相對比較大。這和前面的計算結果是完全吻合的。

2.3 保護方案的設計

為了描述勵磁涌流和故障電流的各子區(qū)間分布系數的差異，本文定義電流波形分布特征系數D如下:

式中，Kd(±1)和 Kd(±N)分別為第±1和±N 子區(qū)間的波形分布系數。

顯然，故障電流的分布特征系數較小，而勵磁涌流的分布特征系數要大得多。

綜合上述分析，本文提出如下的保護方案:

（1）差電流歸一化。以一個工頻周期內差動電流瞬時值的最大絕對值為基準對其進行歸一化處理。

（2）分別計算第±1、±N子區(qū)間的電流波形分布系數 Kd(±1)和 Kd(±N)。

（3）計算波形分布特征系數。當D≥Dset時，判斷為勵磁涌流閉鎖差動保護；當 D＜Dset時，判斷為內部故障開放差動保護跳閘。Dset為判據的門檻值。

要計算波形分布特征系數，首先需要確定劃分的子區(qū)間數目。對于固定的電流波形，如果子區(qū)間數目不同，則每個子區(qū)間內采樣點的數目也會不同，波形分布特征系數也將隨之變化。對于圖 2、圖 4給出的簡化條件下的勵磁涌流和故障電流波形，波形分布特征系數可由下面公式分別計算

從式（6）可知勵磁涌流的分布特征系數還與合閘時的狀態(tài)有關，勵磁涌流越嚴重波形分布特征系數越小。眾所周知當電壓瞬時值過零時刻合閘同時剩磁越大、飽和磁通越小，產生的勵磁涌流越嚴重，本文考慮比較極端的情況，取α=0，kr=1.0，ks=1.1。確定變壓器合閘狀態(tài)參數后，波形分布特征系數的大小就僅與子區(qū)間的數目有關。式(6)、式(7)描述的二者之間關系如圖5所示。

圖5 分布特征系數與子區(qū)間數目關系示意圖Fig.5 Relationship between the distribution characteristic coefficient and the numbers of subintervals

從圖5可以看到隨著子區(qū)間數目的增加，勵磁涌流波形分布特征系數隨之變大，故障電流波形分布特征系數隨之變小，主要是因為對于勵磁涌流，隨著子區(qū)間數目的增加，第±N子區(qū)間分布系數顯然會減小，而由于波形間斷部分采樣點始終處于第±1子區(qū)間內，因而第±1子區(qū)間分布系數減小較少，所以勵磁涌流的波形分布特征系數會變大。而對于故障電流，第±N子區(qū)間分布系數比第±1子區(qū)間分布系數要大得多，因此隨著子區(qū)間數目的增加第±N子區(qū)間分布系數相對減小較少，而第±1子區(qū)間分布系數相對減小較大，因而故障電流的分布特征系數會變小。因此從增大勵磁涌流和故障電流差異角度來看，可以劃分較多的子區(qū)間。另一方面，由于實際上勵磁涌流在變壓器鐵心不飽和階段并非完全等于零，如果子區(qū)間過多，則第±1子區(qū)間分布系數可能減小更多，那么其波形分布特征系數反而會減小。本文綜合考慮這兩方面的因素，取N=10。

3 方法驗證

為了分析提出的新型變壓器勵磁涌流識別判據的性能，本文分別利用大量動模實驗和數字仿真數據[15]對其進行了驗證。

3.1 動模實驗及數字仿真系統簡介

動模實驗系統如圖6所示，其中實驗變壓器由三個單相變壓器組成，采用 Yd11聯結。單相變壓器參數如下:額定容量為10kVA，低壓側額定電壓為 380V，額定電流為 25.3A，高壓側額定電壓為1kV，額定電流為10A，空載電流為1.45%×25.3A，空載損耗為 1%×10kVA，短路損耗為 0.35%×10kVA，短路電壓為（9.0～15.0）%×1kVA，高壓側線圈288匝，由于系統電壓的原因實驗中僅接入78匝，低壓側線圈112匝，全部接入。

圖6 變壓器實驗系統圖Fig.6 Transformer experimental system

動模實驗的內容包括:正常變壓器空載合閘、內部故障變壓器空載合閘、運行中變壓器發(fā)生內部故障等工況。

由于動模實驗中不能控制變壓器鐵心剩磁以及合閘初相角，難以獲得最嚴重的勵磁涌流的情況，同時故障位置也不能任意設置，因此動模實驗數據具有一定的局限性。作為動模實驗的重要補充，本文同時進行了數字仿真，仿真中可任意設置與運行工況有關的各參數，可使結果更全面。數字仿真中變壓器參數與實驗變壓器參數相同。

3.2 分析結果

圖7給出了變壓器正?？胀稌r的勵磁涌流。圖7a、7c分別為偏向負軸和對稱性勵磁涌流，各自的波形分布系數直方圖分別為圖7b、7d，顯然和前面的理論分析結果是相符合的。按式（5）計算得到波形分布特征系數分別為7.1和7.6。

圖7 動模實驗中勵磁涌流及波形分布特征直方圖Fig.7 Inrush currents obtained by dynamic experiment and its waveform distribution coefficient histogram

圖8給出了變壓器發(fā)生B相匝間故障時差電流的波形，其中圖8a、8c分別為帶故障變壓器空載合閘和變壓器運行中發(fā)生故障的情況。從圖8a可以看到空載合閘于故障變壓器時差電流還帶有明顯的涌流特征，如前5ms內的波形呈尖頂波特點，在之后的一段時間內鐵心退出飽和故障電流的特征才呈現出來，圖8b為此時的波形分布系數直方圖，基本符合故障的特點，與前面涌流時的結果有明顯的區(qū)別。圖8c中的故障電流基本符合正弦特征，其波形分布直方圖也與前面理論分析結果十分符合。按式（5）計算得到波形分布特征系數分別為1.0和0.23。

圖8 動模實驗中匝間故障電流及波形分布特征直方圖Fig.8 Turn-to-turn fault currents obtained by dynamic experimental and its waveform distribution coefficient histogram

表1給出了不同運行狀態(tài)下實驗變壓器差電流波形分布特征系數的計算結果。每個運行狀態(tài)下的結果均為多次實驗結果的取值范圍。從表1可得到，正常變壓器空投時波形分布特征系數最小值為2.5，帶故障變壓器空投時該系數最大值為1.38，而變壓器正常運行中發(fā)生故障時該系數最大值僅有0.91，勵磁涌流和故障電流的波形分布特征系數具有較大差異。本文取 Dset=2，即可可靠區(qū)分勵磁涌流和故障電統。

表1 實驗數據的波形分布特征系數計算結果Tab.1 Experimental results of the waveform distribution characteristics coefficient of the differential current waveforms

表2給出了仿真數據的波形分布特征系數、二次諧波含量以及間斷角的計算結果。對于正?？蛰d合閘的情況統計各相的結果，對于變壓器發(fā)生故障的情況只統計故障相計算結果中的最小值，其中對于間斷角判據不統計空投時間斷角較小但同時波寬小于140°可判斷為涌流的情況。

表2 仿真數據的波形分布特征系數（D）、二次諧波含量（I2/I1）和間斷角（θJ）的計算結果Tab.2 Simulation results of waveform distribution characteristics coefficient, second harmonic content and the dead angle

從表2可以看到，仿真數據的波形分布特征系數范圍與實驗數據的結果相符，定值取2時可完全區(qū)分勵磁涌流和故障電流。對于二次諧波制動判據如果定值取15%，對于部分情況包括正常空載合閘和變壓器帶輕微故障空載合閘(表2中用“*”標出)，保護會誤動或延時動作。對于間斷角判據，如果定值取 65°，則表 2中標“**”的情況同樣存在誤動或延時動作的可能。

從上面對比情況可知，波形分布特征判據在性能上優(yōu)于現場普遍采用的二次諧波判據和間斷角判據。主要在于二次諧波判據只能反應波形一個周期的總體特征，間斷角判據僅反應鐵心不飽和時波形出現間斷的局部特點，而波形分布特征判據綜合考慮鐵心不飽和時波形出現間斷以及鐵心飽和時波形呈尖頂波的特征，可以更全面地反應涌流和故障電流之間的差異，因此性能得以提高。

4 結論

本文利用直方圖技術分析了變壓器勵磁涌流和內部故障狀態(tài)下差動電流波形時域分布特征，在此基礎上提出一種基于波形分布特征的新型勵磁涌流識別方法，并給出了具體的實現方案。該方法可以綜合考慮變壓器差電流在時域上的局部分布特征，更好地反應了勵磁涌流和故障電流之間的差異。實驗數據和仿真數據分析結果表明本方法具有更高的可靠性和靈敏度，性能優(yōu)于二次諧波判據和間斷角判據。

[1]孫志杰, 陳云侖.波形對稱原理的變壓器差動保護[J].電力系統自動化, 1996, 20(4): 42-46.Sun Zhijie, Chen Yunlun.Transformer differential protection based on the characteristic analysis of the first and second half cycle of the magnetizing inrush current [J].Automation of Electric Power Systems,1996, 20(4): 42-46.

[2]林湘寧, 劉世明, 楊春明, 等.幾種波形對稱法變壓器差動保護原理的比較研究[J].電工技術學報,2001, 16(4): 44-49.Lin Xiangning, Liu Shiming, Yang Chunming, et al.Study on comparisons among some waveform symmetry principle based transformer differential protection[J].Transactions of China Electrotechnical Society, 2001, 16(4): 44-49.

[3]陳德樹, 尹項根, 張哲, 等.虛擬三次諧波制動式變壓器差動保護[J].中國電機工程學報, 2001,21(8): 19-23.Chen Deshu, Yin Xianggen, Zhang Zhe, et al.Virtual third harmonic restrained transformer differential protection principle and practice[J].Proceedings of the CSEE, 2001, 21(8): 19-23.

[4]鄭濤, 劉萬順, 肖仕武, 等.一種基于數學形態(tài)學提取電流波形特征的變壓器保護新原理[J].中國電機工程學報, 2004, 24(7): 18-24.Zheng Tao, Liu Wanshun, Xiao Shiwu, et al.A new algorithm based on the mathematical morphology for power transformer protection[J].Proceedings of the CSEE, 2004, 24(7): 18-24.

[5]焦邵華, 劉萬順, 劉建飛, 等.用小波理論區(qū)分變壓器的勵磁涌流和短路電流的新原理[J].中國電機工程學報, 1999, 19(7): 1-5.Jiao Shaohua, Liu Wanshun, Liu Jianfei, et al.A new principle of discrimination between inrush current and fault current of transformer based on wavelet [J].Proceedings of the CSEE, 1999, 19(7): 1-5.

[6]馬靜, 王增平, 徐巖.用相關函數原理識別變壓器勵磁涌流和短路電流的新方法[J].電網技術, 2005,29(6): 78-81.Ma Jing, Wang Zengping, Xu Yan.A novel method to identify inrush current and short circuit current based on correlation function[J].Power System Technology,2005, 29(6): 78-81.

[7]Yabe K.Power differential method for discrimination between fault and magnetizing inrush current in transformer[J].IEEE Transactions on Power Delivery,1997, 12(3): 1109-1118.

[8]孫志杰, 曾獻華, 湯漢松.磁通制動原理在變壓器差動保護中的應用[J].電力自動化設備, 2005,25(1): 79-81.Sun Zhijie, Zeng Xianhua, Tang Hansong.Application of flux restraint theory in transformer differential protection[J].Electric Power Automation Equipment, 2005, 25(1): 79-81.

[9]徐巖, 王增平, 楊奇遜, 等.基于電壓電流微分波形特性的變壓器保護新原理的研究[J].中國電機工程學報, 2004, 24(2): 61-65.Xu Yan, Wang Zengping, Yang Qixun, et al.Research on novel transformer protection based on the characteristics of voltage and differential current[J].Proceedings of the CSEE, 2004, 24(2): 61-65.

[10]葛寶明, 于學海, 王祥珩, 等.基于等效瞬時電感判別變壓器勵磁涌流的新方法[J].電力系統自動化,2004, 28(7): 44-48.Ge Baoming, Yu Xuehai, Wang Xiangheng, et al.A novel equivalent instantaneous inductance based algorithm used to distinguish inrush currents for transformers[J].Automation of Electric Power Systems, 2004, 28(7): 44-48.

[11]馬靜, 王增平, 王雪.基于等效瞬時漏感的變壓器保護新原理[J].電力系統自動化, 2006, 30(23):64-68.Ma Jing, Wang Zengping, Wang Xue.Novel principle of power transformer protection based on equivalent instantaneous leakage inductance[J].Automation of Electric Power Systems, 2006, 30(23): 64-68.

[12]王增平, 徐巖, 王雪, 等.基于變壓器模型的新型變壓器保護原理的研究[J].中國電機工程學報,2003, 23(12): 54-58.Wang Zengping, Xu Yan, Wang Xue, et al.Study on the novel transformer protection principle based on the transformer model[J].Proceedings of the CSEE,2003, 23(12): 54-58.

[13]張保會, 尹項根.電力系統繼電保護[M].北京: 中國電力出版社, 2005.

[14]岡薩雷斯, 等.數字圖像處理[M].2版.阮秋琦, 等譯.北京: 電子工業(yè)出版社, 2003.

[15]王雪, 王增平.基于有限元法的變壓器電感參數計算方法的研究[J].電力系統保護與控制, 2009,37(24): 11-14.Wang Xue, Wang Zengping.An investigation on transformer inductance calculation based on finite element method[J].Power System Protection and Control, 2009, 37(24): 11-14.