，

(1.三峽大學(xué)電氣與新能源學(xué)院，湖北 宜昌 443002；2.宜昌供電局檢修公司，湖北 宜昌 443000)

行波小波分析法對(duì)單相接地故障選線方法的研究

賈智彬1，顧佩穎2

(1.三峽大學(xué)電氣與新能源學(xué)院，湖北 宜昌 443002；2.宜昌供電局檢修公司，湖北 宜昌 443000)

論述了配電網(wǎng)中單相接地短路故障行波的特性，運(yùn)用小波變換分析中性點(diǎn)非直接接地配電網(wǎng)中單相接地故障行波信號(hào)。導(dǎo)出了一種利用暫態(tài)行波0模量的小波系數(shù)的極大值極性識(shí)別故障線路的故障選線新方法。并證明了在故障初始行波過(guò)程中，故障點(diǎn)間歇電弧及過(guò)渡電阻對(duì)本方法的影響十分微弱。通過(guò)大量仿真試驗(yàn)，證明本方法的正確性、有效性和通用性。

行波;小波分析;配電網(wǎng);故障選線;極大值極性

1 引言

我們知道，配電網(wǎng)系統(tǒng)的中性點(diǎn)通常是不直接接地的(如不接地或者經(jīng)消弧線圈接地)。當(dāng)發(fā)生單相接地故障時(shí)，配電網(wǎng)系統(tǒng)往往允許繼續(xù)運(yùn)行一段時(shí)間，但一般不超過(guò)2小時(shí)。此時(shí)，故障點(diǎn)的間歇電弧能夠在系統(tǒng)中產(chǎn)生可能高達(dá)3.5倍相電壓的過(guò)電壓，這對(duì)于運(yùn)行人員和系統(tǒng)設(shè)備都十分危險(xiǎn)，并可能留下一些事故隱患。與連接在同一母線上的其他線路相比，故障線路中的穩(wěn)態(tài)殘余電流的大小和相位并沒(méi)有顯著的特征，所以，一直以來(lái)，判別故障線路十分困難。

一些文獻(xiàn)已經(jīng)提出了幾種中性點(diǎn)不直接接地系統(tǒng)中單相接地故障選線的方法，并多半是基于零序電流。文獻(xiàn)[1]認(rèn)為電流相位有別于其他線路的就是故障線路。文獻(xiàn)[2]運(yùn)用小波分析實(shí)現(xiàn)這種方法。文獻(xiàn)[3]把電流幅值最大的判定為故障線路。文獻(xiàn)[4]同時(shí)判別電流的大小和相位以提高選線準(zhǔn)確性。但是，如果系統(tǒng)中各線路的RLC參數(shù)有些差異或者中性點(diǎn)是經(jīng)消弧線圈接地的，這些方法就可能失效。而且，故障線路中的零序電流常常很小，很難精確獲取，所以，這些選線原理總是難以令人滿意。開(kāi)發(fā)能夠應(yīng)用在各種中性點(diǎn)非直接接地系統(tǒng)進(jìn)行準(zhǔn)確故障選線的新原理新方法，特別是適用于經(jīng)消弧線圈接地的系統(tǒng)，這仍然十分迫切。

2 配電網(wǎng)系統(tǒng)中的暫態(tài)行波

在對(duì)稱(chēng)三相輸電線路中的x點(diǎn)，某時(shí)刻t的電壓、電流行波應(yīng)滿足下列波動(dòng)方程：

(1)

其中，L、C分別為單位長(zhǎng)度線路的電感和電容矩陣[6，7]。

為求解此波動(dòng)方程，通常進(jìn)行相模變換，把相互耦合的三相電壓電流變換為相互獨(dú)立的模量，即α模量、β模量和0模量。

(2)

可以采用Karen Bauer變換：

圖1 初始波過(guò)程中的電流行程

(3)

這樣，式(1)就可以變換為：

(4)

其中，[Lc]=diag[L1C1L1C1L1C1L1C1L1C1L1C1]，L1、L0、C1，C0分別為線路單位長(zhǎng)度的正序和零序電感和電容。

求解式(4)，得到α、β、0模量電壓和電流的前向行波和反向行波波動(dòng)方程：

(5)

其中，ufwd、ifwd、ubwd和ibwd分別為前向行波和反向行波電壓、電流。

當(dāng)某條線路上發(fā)生單相接地故障時(shí)，故障點(diǎn)產(chǎn)生的暫態(tài)行波由故障線路以電磁波速向變電站母線傳播，并折射到母線上的其他線路和變壓器等元件中去，在非故障線路等元件中形成前行波。在初始行波過(guò)程期間，母線電壓為ubus，各非故障線路中的電流前行波為：

(6)

n為線路編號(hào)。在初始行波過(guò)程中，各線路中電流行波的情形如圖1所示。

圖2 在初始時(shí)間t=0+變壓器繞組的近似等效

3 考慮變壓器中的影響

變壓器及其中性點(diǎn)接地方式對(duì)母線行波電壓可能會(huì)有些影響，其效果在現(xiàn)有文獻(xiàn)中都沒(méi)有分析。在初始行波過(guò)程中，由于變壓器繞組電感很大，繞組中電流不能突變，無(wú)論繞組末端是否接地，都可以認(rèn)為電感回路開(kāi)路。所以，首先可以把單相變壓器繞組近似等效為其入口電容Ct都可以認(rèn)為，如圖2所示。繞組上x(chóng)點(diǎn)對(duì)地初始電位分布可近似為：

u=ubuse-αx

(7)

在初始時(shí)刻t=0+繞組端部(x=0)的電位為ubus。由彼得森法則，在初始行波過(guò)程中，從母線到變壓器的連接線中的電流行波可以近似為：

(8)

其中，zlink為連接線波阻抗，Ct為變壓器繞組近似等效為其入口電容。

在t=0+時(shí)刻，式(8)中的電流可以進(jìn)一步等效為：

(9)

這與一條非故障線路中的電流前向行波十分相似。

在單相接地故障初始行波過(guò)程中，中性點(diǎn)接地的Y形接線變壓器中電壓、電流行波的特性與單相變壓器規(guī)律相同，中性點(diǎn)不接地或經(jīng)消弧線圈接地的Y形接線以及無(wú)中性點(diǎn)的?形接線變壓器，包括電壓互感器，其行波特性都與上述繞組末端接地時(shí)的單相變壓器規(guī)律相似[6]。

因此，可以得出一個(gè)十分慎重的結(jié)論：無(wú)論變壓器的中性點(diǎn)是否接地，在單相接地故障時(shí)的暫態(tài)行波到達(dá)變壓器瞬間的初始行波過(guò)程中，從母線到變壓器的連接線中的電流行波可以近似為一條非故障線路中的電流前向行波。此電流行波與所有非故障線路中的電流前向行波之和為故障線路的電流行波。所以如果以u(píng)bus為參考方向，故障線路中的電流行波與所有非故障線路中的電流行波的方向正好相反。這一點(diǎn)對(duì)于我們構(gòu)成單相接地故障選線的有效的通用算法十分有益，也至關(guān)重要。

4 單相接地故障選線算法

當(dāng)配電網(wǎng)系統(tǒng)中發(fā)生單相接地故障時(shí)，在故障點(diǎn)將產(chǎn)生暫態(tài)電壓電流行波。根據(jù)上文的分析，如果以母線電壓ubus為參考方向，故障線路中的電流行波與同一母線上所有非故障線路、變壓器及電壓互感器中的電流行波的方向正好相反。本文提取行波的0模量構(gòu)成單相接地故障選線算法。

暫態(tài)行波信號(hào)是一種非穩(wěn)定的奇異性信號(hào)，特別當(dāng)故障點(diǎn)存在間歇電弧時(shí)。所以，常規(guī)信號(hào)分析方法很難有效提取故障特征信息甚至完全不適用，如快速傅里葉變換FFT、FIR濾波器、Kalman算法等。小波分析是理想的奇異信號(hào)分析工具。當(dāng)故障發(fā)生時(shí)刻，母線電壓和線路電流波形都出現(xiàn)不連續(xù)點(diǎn)，所以其小波系數(shù)會(huì)在此時(shí)刻同時(shí)出現(xiàn)極大值點(diǎn)。

通過(guò)比較母線電壓和各線路電流初始行波0模量小波系數(shù)極大值極性，本文構(gòu)成了一種新穎的通用的基于小波分析的故障選線算法。其流程圖如圖3所示。

圖3 故障選線程序流程圖

5 應(yīng)用

對(duì)圖1所示的配電網(wǎng)系統(tǒng)，本文建立了一個(gè)基于MATLAB的仿真模型，其中變壓器、線路都采用了分布參數(shù)模型，系統(tǒng)會(huì)在定義的時(shí)刻發(fā)生單相接地接地故障。圖4所示為線路1在仿真時(shí)間0.008s發(fā)生單相接地接地故障時(shí)的母線零序電壓和線路零序電流波形。由于系統(tǒng)消弧線圈為過(guò)補(bǔ)償方式，所以，可以看到，故障線路與非故障線路的穩(wěn)態(tài)零序電流波形十分相似。這一特征會(huì)使所有基于穩(wěn)態(tài)故障電壓電流的選線方法失效。

在仿真計(jì)算時(shí)，采用了定步長(zhǎng)離散解法，計(jì)算時(shí)間步長(zhǎng)為1×10-7s，線路長(zhǎng)度、RLC參數(shù)、故障線路、故障位置、故障時(shí)間、接地電阻等都可以根據(jù)要求進(jìn)行定義。從系統(tǒng)穩(wěn)態(tài)開(kāi)始運(yùn)行仿真模型，測(cè)量母線電壓和各線路電流，識(shí)別故障類(lèi)型和發(fā)生時(shí)間，并截取故障前后0.2ms的數(shù)據(jù)。

MATLAB提供了幾十種小波函數(shù)，可供分析這些故障數(shù)據(jù)。本文采用了Haar小波。圖5所示為在尺度50下的Haar連續(xù)小波變換(CWT)系數(shù)波形，第一個(gè)波形為母線電壓行波0模量CWT系數(shù)，其他為各線路及變壓器連接線中電流行波0模量CWT系數(shù)。在故障時(shí)刻，所有這些CWT系數(shù)都有一個(gè)極大值，而且線路1的電流行波0模量CWT系數(shù)極大值極性與其他的正好相反。于是，可以判定線路1為故障線路。

在采樣周期為1×10-7s時(shí)的Haar連續(xù)小波變換199.22kHz。圖6所示為在尺度1～100(相當(dāng)于頻域中的9.96MHz～110.7kHz)下的母線電壓、故障線路電流、非故障線路電流和變壓器連接線電流的行波0模量Haar CWT系數(shù)的3D圖形。其中的CWT系數(shù)極大值極性在很大的范圍里都符合上述的規(guī)律，都能夠完成故障線路的判定。

本文進(jìn)一步計(jì)算了這些行波0模量的Haar小波包系數(shù)，在小波樹(shù)的很多細(xì)節(jié)節(jié)點(diǎn)(高頻部分)上，其系數(shù)極大值極性都十分顯著地符合上述規(guī)律。所以，本文提出了一個(gè)實(shí)用的故障選線方法，我們比較多個(gè)尺度下的連續(xù)小波變換和多個(gè)節(jié)點(diǎn)的小波包的系數(shù)極大值極性，并對(duì)比較結(jié)果進(jìn)行表決得出最終的選線結(jié)果。這樣，就可以有效地避免選線錯(cuò)誤，保證選線可靠性。

在不同的條件下運(yùn)行此配電網(wǎng)系統(tǒng)仿真模型，并應(yīng)用本文的選線方法進(jìn)行故障線路的判斷。本文主要進(jìn)行了下列條件下的單相接地故障仿真試驗(yàn)：

圖4 母線零序電壓及線路零序電流

圖5 母線電壓和線路電流0模量在尺度50下的邊續(xù)小波變換系數(shù)

圖6 母線電壓和線路電流0模量在尺度1～100下的邊續(xù)小波變換系數(shù)的3D圖形

(1)系統(tǒng)處于不同的運(yùn)行方式：變壓器中性點(diǎn)經(jīng)電阻接地、不接地或經(jīng)消弧線圈接地(過(guò)補(bǔ)償方式、欠補(bǔ)償方式、完全補(bǔ)償方式)；系統(tǒng)負(fù)荷輕載或者重載。

(2)輸電線路條件：線路數(shù)量不同(2、3或者更多條線路)；不同的導(dǎo)線RLC參數(shù)(架空線或者電纜線)及長(zhǎng)度。

(3)故障點(diǎn)條件：不同的故障時(shí)間、故障線路、故障點(diǎn)位置(從首端到末端)，并帶有不同的過(guò)渡電阻(從0.0001～500Ω)。

6 結(jié)論

本文以小波分析為數(shù)學(xué)工具，分析中性點(diǎn)非直接接地的配電網(wǎng)中單相接地故障時(shí)的行波信號(hào)。小波分析對(duì)于暫態(tài)奇異性信號(hào)的適用性是毋庸置疑的。本文導(dǎo)出了一種新穎的基于行波和小波的單相接地故障選線方法，證明了在行波初始過(guò)程中，變壓器(包括電壓互感器)及其中性點(diǎn)接地方式的影響是很小的。這一故障選線方法運(yùn)用故障前后0.2ms的數(shù)據(jù)，此時(shí)間十分短暫，所以，故障點(diǎn)間歇電弧及其可變的電阻的影響可以不必考慮。

此故障選線方法適用于各種接地方式的配電網(wǎng)，大量仿真試驗(yàn)證明了其準(zhǔn)確性和可靠性。此故障選線方法完全能夠設(shè)計(jì)為實(shí)時(shí)執(zhí)行的程序，只是因?yàn)椴蓸宇l率和行波的頻率都很高，就需要很大的系統(tǒng)內(nèi)存開(kāi)銷(xiāo)和較高的電磁兼容性能。

[1] Chaari,P.Bastard,and M.Meunier,“Prony′s method:an efficient tool for the analysis of earth fault currents in Petersen-coil-protected networks.” IEEE Trans.Power Delivery,1995,10(3):1234-1241.

[2] Chaari,M.Meunier,and F.Brouaye,“Wavelet:A New Tool for the Resonant Grounded Power Distribution Systems Relaying.” IEEE Trans.on Power Delivery,1996,11(3):1301-1308.

[3] Liang Jun,et al,“A method of fault line detection in distribution systems based on wavelets”,Proc.2002 International Conference on Power System Technology,2002,4(13-17):2635-2639.

[4] Jia Qingquan,et al,“Multi-criteria relaying strategy for single phase to ground fault in MV power systems”,in Proc.2002 International Conference on Power System Technology,2002,2(13-17):683-687.

[5] 崔錦泰,程正興.小波分析導(dǎo)論[M].西安交通大學(xué)出版社,1995.1.

[6] 周澤存,高電壓技術(shù)[M].水利電力出版社,1988.6.

[7] 解廣潤(rùn),電力系統(tǒng)過(guò)電壓[M].水利電力出版社,1985.6.

2013-07-02

賈智彬(1970-)，男，工程師，主要從事于絕緣與過(guò)電壓研究；

顧佩穎(1974-)，女，工程師，主要從事于電力系統(tǒng)過(guò)電壓研究。

ResearchonthePhaseEarthFaultFeederDetectionbyTravellingWaveandWavelets

JIAZhi-bin1,GUPei-ying2

(1.Three Gorges University,Yichang 443002,China;2.Power Supply Bureau of Yichang，Yichang 443000,China)

In the paper the characteristic of single phase fault traveling waves in distribution grid is discussed,wavelets have been presented to analyze single phase fault traveling wave signals in an ineffectively grounded distribution grid.A novel earth failure detection technique that distinguishes the fault feeder by the wavelets coefficient maximum polarity of the traveling wave zero mode is introduced.In the primary process the influences of fault intermittent arc,transformer and its neutral point are proved be feeble.The exactness,availability and universality of the technique are proved by a great deal of simulating test.

traveling wave；wavelets；earth fault detection；maximum polarity

1004-289X(2013)06-0019-05

TM72

B