(國(guó)網(wǎng)天津?yàn)I海供電公司，天津 300450)

0 引 言

當(dāng)輸電線路發(fā)生故障后，在故障電氣量中會(huì)包含一種頻率很高的暫態(tài)量信號(hào)，該信號(hào)被稱為行波信號(hào)，其包含了與故障情況相關(guān)的豐富信息。正確地從行波信號(hào)中識(shí)別和提取與故障相關(guān)的信息可以為快速準(zhǔn)確地識(shí)別故障相、選擇性地切除故障相提供依據(jù)。在繼電保護(hù)現(xiàn)有工程實(shí)際和研究中，為了保證選相的快速性，故障選相方法采用行波信號(hào)的特征直接構(gòu)成判別方法。目前對(duì)于采用行波信號(hào)進(jìn)行故障相別判別的方法主要包括以下幾種：1)采用行波方向判別方法作為故障相別判別的方法，當(dāng)根據(jù)方向判別方法判斷某相為正方向時(shí)，該相即判斷為故障相。這種選相判別方法在RALDA型行波保護(hù)中[1]得到了應(yīng)用。2)基于各相電流行波峰值進(jìn)行故障選相的方法，通過(guò)比較檢測(cè)期間三相行波電流故障分量峰值與電流最大的一相的故障分量最大值進(jìn)行故障相判別。當(dāng)任意一相電流的最大值大于電流最大一相的故障分量最大值的1/2時(shí)，判斷該相為故障相。這種選相方法是對(duì)選相方法1)的改進(jìn)，在工程實(shí)踐中已經(jīng)有應(yīng)用[2]。3)基于電流行波比較原理的選相方法，通過(guò)比較三相行波電流之差(IA-IB，IB-IC，IC-IA)的極性進(jìn)行判別。當(dāng)兩兩之間的極性相反時(shí)，相鄰元件輸出正值并判定為故障相[2]。

上述方法通過(guò)利用故障時(shí)行波信號(hào)中的暫態(tài)量信息，為故障時(shí)選相提供了思路，但在選相正確性方面仍存在一些問(wèn)題。其中，方法1)和方法2)利用行波方向判別方法進(jìn)行選相，無(wú)法避免因不同相間的電磁耦合引起的選相不正確問(wèn)題。同時(shí)，與方法3)一樣，由于選相元件采用故障時(shí)行波瞬時(shí)值的特征進(jìn)行故障相選擇，而瞬時(shí)值受波形畸變、瞬變干擾的影響較大，容易造成分析結(jié)果的不可靠。此外，由于故障發(fā)生地點(diǎn)具有隨機(jī)不可控的特性，這將造成行波波頭在瞬時(shí)采樣過(guò)程中采樣值的數(shù)目不易確定。

小波變換因其在信號(hào)處理過(guò)程中具有使信號(hào)時(shí)頻特性局部化的優(yōu)點(diǎn)，近年來(lái)在數(shù)學(xué)和數(shù)字信號(hào)處理領(lǐng)域得到了充分的應(yīng)用和發(fā)展[3-5]。在電力行業(yè)應(yīng)用方面，小波變換被廣泛應(yīng)用于行波信號(hào)的分析。這里的基本出發(fā)點(diǎn)是充分利用小波變換在暫態(tài)信號(hào)處理上的優(yōu)點(diǎn)，避免選相受電磁耦合、波形畸變、瞬變等干擾的影響?？紤]到行波信號(hào)在到達(dá)測(cè)量點(diǎn)時(shí)的突變性質(zhì)以及小波變換具有將行波信號(hào)的特征進(jìn)行時(shí)頻局部化的作用[6-7]，因此結(jié)合行波信號(hào)和B樣條二進(jìn)小波變換，對(duì)電流行波信號(hào)進(jìn)行時(shí)頻分析，通過(guò)模故障分量構(gòu)成更完善、更快速、更準(zhǔn)確的故障選相方法。首先對(duì)B樣條二進(jìn)小波的優(yōu)點(diǎn)進(jìn)行分析，然后運(yùn)用B樣條二進(jìn)小波對(duì)解耦變換后的模電流分量進(jìn)行分析，計(jì)算得到模電流分量初始行波的模極大值，并根據(jù)模電流分量極大值在不同故障類型下的不同特征進(jìn)行對(duì)比分析從而構(gòu)建出故障選相原理。最后，通過(guò)仿真驗(yàn)證選相原理的正確性。

1 B樣條小波及其優(yōu)點(diǎn)

在小波分析中，由于B樣條小波具有優(yōu)良的光滑性和緊支撐性，因此，在現(xiàn)代行波故障分析中一般使用B樣條函數(shù)的導(dǎo)函數(shù)作為基小波函數(shù)。B樣條函數(shù)是一種半正定的、分段光滑的函數(shù)，m次B樣條函數(shù)的表達(dá)式為

βm(x)=βm-1(x)·β0(x)，m≥1

(1)

其中：

(2)

由表達(dá)式可知B樣條函數(shù)是分段多項(xiàng)式，這使得其處理非常簡(jiǎn)化。此外B樣條函數(shù)具有對(duì)稱性，可以由它構(gòu)造出對(duì)稱或反對(duì)稱的樣條小波，且可以達(dá)到非常高的近似精度。

以B樣條函數(shù)構(gòu)造的二進(jìn)小波變換具有平移不變性，使信號(hào)奇異點(diǎn)的檢測(cè)變得簡(jiǎn)單。當(dāng)行波到達(dá)信號(hào)測(cè)量點(diǎn)時(shí)，測(cè)量點(diǎn)將測(cè)得一發(fā)生“突變”的信號(hào)，對(duì)該信號(hào)進(jìn)行小波變換，即可獲得對(duì)應(yīng)于信號(hào)突變的模極大值。因此，借助B樣條二進(jìn)小波變換即可將行波信號(hào)的分析轉(zhuǎn)化為其小波變換下的模極大值的分析，從而使基于行波信號(hào)的故障選相簡(jiǎn)單易行。

根據(jù)已有文獻(xiàn)的分析結(jié)論有：三次B樣條函數(shù)對(duì)提取含噪聲信號(hào)的奇異性是漸進(jìn)最優(yōu)的。因此，采用三次樣條函數(shù)作為二進(jìn)小波變換的基函數(shù)，其對(duì)應(yīng)表達(dá)式為

(3)

對(duì)故障時(shí)的行波信號(hào)進(jìn)行三次B樣條二進(jìn)小波變換即可獲得故障時(shí)行波信號(hào)特征。

2 故障電流的解耦變換與小波變換

當(dāng)輸電線路上發(fā)生故障時(shí)，可運(yùn)用小波變換對(duì)量測(cè)點(diǎn)的信號(hào)進(jìn)行處理得到該信號(hào)的模極大值，通過(guò)模極大值即可對(duì)該信號(hào)進(jìn)行表征。運(yùn)用模極大值的表征作用即可通過(guò)工頻電流信號(hào)的分析實(shí)現(xiàn)基于模故障分量的故障選相，從而將工頻量信號(hào)應(yīng)用于行波選相中，實(shí)現(xiàn)僅通過(guò)故障時(shí)電流的初始行波判斷故障相別的目的。

由于三相輸電線路中各相電流不獨(dú)立，因此無(wú)法完全消除非故障相對(duì)故障相電流的影響，需要通過(guò)解耦變換得到獨(dú)立的模電流分量進(jìn)行故障選相。以下假定在故障發(fā)生后母線流向故障點(diǎn)的三相故障電流分別為IA、IB、IC，計(jì)算在不同類型故障條件下的母線流向故障點(diǎn)的行波電流。運(yùn)用相模變換對(duì)三相故障電流進(jìn)行解耦變換與分析，建立非獨(dú)立的三相電流與解耦的模電流之間的數(shù)量關(guān)系。相模變換可采用的變換方法主要包括Karenbauer變換、Clarke變換等幾種常用的解耦變換方法。以Karenbauer變換為例建立故障后三相電流的模電流分量，Karenbauer變換因子的表達(dá)式為

(4)

利用上述Karenbauer變換因子對(duì)故障電流進(jìn)行變換可得三相故障電流與各模電流的關(guān)系如下：

(5)

式中：零模電流分量I0通過(guò)大地和導(dǎo)線構(gòu)成電流回路，該回路的波阻抗相對(duì)較大，波傳播的速度相對(duì)較小，傳播過(guò)程中的衰耗相對(duì)較大；線模電流分量Iα、Iβ通過(guò)導(dǎo)線構(gòu)成電流回路，該回路的波阻抗相對(duì)較小，波傳播的速度與光速基本接近，且比較穩(wěn)定，傳播過(guò)程中的衰耗相對(duì)較小。為便于故障選相，引入“γ模”：

(6)

γ模不是獨(dú)立的，它由α模和β模的線性組合而得到。需要注意的是，γ模的定義是建立在三相平衡系統(tǒng)上的。根據(jù)式(5)、式(6)和故障邊界條件可以分別得到對(duì)應(yīng)于接地型短路、非接地型短路和三相短路時(shí)各電流模量的情況，如表1所示。

在獲得解耦的各模量電流后即可進(jìn)行B樣條二進(jìn)小波變換，對(duì)故障后的模量電流進(jìn)行特性分析。

3 故障特征分析及選相原理

根據(jù)上述分析可以得到發(fā)生不同類型故障及不同相故障時(shí)的故障電流模分量的特征。

1)故障情況下的零模分量

在無(wú)干擾的情況下，如果線路發(fā)生接地故障，此時(shí)的零模分量I0不等于0，而如果發(fā)生的故障為相間或者三相短路時(shí)零模分量I0等于0。因而可以根據(jù)零模電流的大小判斷故障為接地故障還是相間或三相短路。

2)各種情況下的線模分量

兩相相間短路時(shí)，線模電流滿足關(guān)系式max(|3Iα|，|3Iβ|，|3Iγ|)=2min(|3Iα|，|3Iβ|，|3Iγ|)，再根據(jù)最大模量電流及對(duì)應(yīng)線模電流可以確定故障相。

單相短路時(shí)，3個(gè)線模電流中有1個(gè)為0，其他2個(gè)線模電流相等。根據(jù)模量電流為0的線模電流可以確定故障相。

兩相接地故障時(shí)，零模電流與線模電流存在3種情況：當(dāng)滿足I0=Iβ+Iγ時(shí)，兩相接地故障的故障相為AB相；當(dāng)滿足I0=Iα-Iγ時(shí)，兩相接地故障的故障相為CA相；當(dāng)滿足I0=-Iβ-Iα?xí)r，兩相接地故障的故障相為BC相。

綜上，基于B樣條二進(jìn)小波變換模極大值的故障相判別原理及其過(guò)程可總結(jié)為：通過(guò)三相電流計(jì)算得到故障情況下各個(gè)模量電流初始行波的B樣條二進(jìn)小波模極大值；通過(guò)電流模量中零模電流的小波模極大值與0的大小關(guān)系判斷是否為接地故障，不為0時(shí)為接地故障，否則為非接地故障；通過(guò)電流模量中線模分量的關(guān)系可以判斷故障時(shí)的相別。

表1 不同故障類型和故障相下的模量電流

故障類型、故障相與零模、線模電流的對(duì)應(yīng)關(guān)系如式(7)、式(8)、式(9)所示。

(7)

(8)

(9)

不滿足上述3種條件者，判為三相故障。故障類型的判別與選相過(guò)程如圖1所示。

圖1 故障類型判別與選相流程

理論上可以根據(jù)上述等式條件準(zhǔn)確判別出故障類型及故障相，然而實(shí)際信號(hào)的提取和處理過(guò)程中存在各種誤差和近似，因而需將上述判別故障類型和故障相的等式條件變?yōu)椴坏仁綏l件，即當(dāng)?shù)仁絻蛇叺牟钪翟谝欢ǖ恼`差范圍內(nèi)即認(rèn)為等式成立，從而可在實(shí)際工程中實(shí)現(xiàn)上述判別方法。

4 仿真驗(yàn)證

為了驗(yàn)證了上述故障選相方法在不同故障類型和故障相別情況下的準(zhǔn)確性和優(yōu)越性，針對(duì)圖2所示的仿真系統(tǒng)，采用PSCAD軟件對(duì)其進(jìn)行分析，獲得系統(tǒng)故障時(shí)測(cè)量端處的暫態(tài)電流行波信號(hào)。并將PSCAD中獲得的數(shù)據(jù)導(dǎo)入Matlab中進(jìn)行數(shù)據(jù)處理，根據(jù)處理后的結(jié)果進(jìn)行故障選相。數(shù)據(jù)處理采用基于三次B樣條的小波函數(shù)進(jìn)行小波變換，獲得所需各個(gè)模量的模極大值。

圖2 仿真系統(tǒng)

仿真中每隔0.000 002 s取一個(gè)點(diǎn)，即采樣頻率為0.5 MHz。以A相接地短路為例，故障發(fā)生在距離M側(cè)母線90 km處，過(guò)渡電阻為0.01 Ω。保護(hù)安裝處獲得的A相接地短路電流及其模電流、小波變換結(jié)果如圖3所示(故障發(fā)生在第1096個(gè)采樣點(diǎn)處)。根據(jù)零模電流不為0，可知故障類型為接地故障，又由于I0≈Iα≈Iβ，Iγ≈0，對(duì)比表1和式(7)可知該故障類型為單相接地短路，故障相為A相。

同樣地，利用B樣條小波變換對(duì)其他單相接地短路、相間短路、相間接地短路和三相短路進(jìn)行分析，可以得到對(duì)應(yīng)不同故障類型下的零模電流和線模電流的小波模極大值如表2所示。對(duì)比表1和式(7)、式(8)、式(9)的故障模量特征即可由模電流的小波模極大值判斷故障類型。

5 結(jié) 語(yǔ)

基于B樣條二進(jìn)小波變換，提出了一種根據(jù)暫態(tài)行波信號(hào)特征進(jìn)行故障相別判斷的方法，該方法利用模電流極大值數(shù)據(jù)進(jìn)行故障相別判別，充分顯示了選相時(shí)故障特征量的顯著性，使選相的物理概念更加清晰；故障選相判據(jù)執(zhí)行流程簡(jiǎn)單，選相邏輯清晰、速度快、可靠性高。該選相方法對(duì)故障點(diǎn)的過(guò)渡電阻、故障前的負(fù)荷電流、系統(tǒng)阻抗等因素具有很強(qiáng)的魯棒性。此外，由于選相過(guò)程中合理利用了故障電流中的零模分量，使該方法在接地故障電流波形嚴(yán)重畸變的情況下也能夠正確選相，克服了現(xiàn)有故障選相判別方法因相間電磁耦合作用造成選相錯(cuò)誤的問(wèn)題。

圖3 A相故障電流及各模電流小波變換結(jié)果

故障類型故障相模電流小波模極大值/kA3I03Iα3Iβ3Iγ單相接地A-0.081 8-0.076 8-0.075 50B0.335 7-0.312 600.312 8C-0.254 11.406 0×10-30.236 20.237 3兩相接地AB3.442 0×10-13-0.451 3-0.224 90.226 4BC4.530 0×10-13-0.319 40.318 60.638 0AC2.761 4×10-130.092 20.186 90.094 6兩相接地AB0.223 0-0.451 3-0.121 40.329 9BC-0.462 70.219 40.247 3-0.101 6AC-0.294 6-0.045 70.186 90.231 2三相接地ABC5.587 0×10-13-0.451 30.186 90.638 0