楊關(guān)春,彭 麗
(1.山東理工大學(xué)電氣與電子工程學(xué)院,山東淄博255091;2.山東省電力集團(tuán)公司臨沂供電公司,山東臨沂276000)
一種不受波速影響的單端行波故障測(cè)距新算法
楊關(guān)春1,彭 麗2
(1.山東理工大學(xué)電氣與電子工程學(xué)院,山東淄博255091;2.山東省電力集團(tuán)公司臨沂供電公司,山東臨沂276000)
采用提升小波變換和奇異性檢測(cè)原理對(duì)電流行波到達(dá)母線測(cè)量端的時(shí)刻進(jìn)行檢測(cè).首先利用初始行波、故障點(diǎn)反射波及對(duì)端母線反射波計(jì)算出行波波速,然后在常規(guī)單端行波故障測(cè)距算法的基礎(chǔ)上,通過(guò)消去行波波速參數(shù),提出一種不受行波波速影響的單端行波故障測(cè)距新算法. PSCAD和MATLAB仿真驗(yàn)證結(jié)果證明該方法具有很高的測(cè)距精度和實(shí)用性.
單端測(cè)距;行波;波速;提升小波;奇異性檢測(cè)
隨著我國(guó)經(jīng)濟(jì)的發(fā)展,電力系統(tǒng)規(guī)模逐漸加大,網(wǎng)絡(luò)結(jié)構(gòu)逐漸復(fù)雜,用戶(hù)對(duì)供電穩(wěn)定的要求也越來(lái)越高.因此,如何確保配電網(wǎng)的安全可靠運(yùn)行,快速準(zhǔn)確地查出故障線路及故障點(diǎn)位置,具有非常重要的意義.目前,輸電線路測(cè)距方法主要有行波法和阻抗法.與阻抗法相比,行波法不受系統(tǒng)運(yùn)行方式、接地故障電阻、線路參數(shù)和故障類(lèi)型等因素的影響[1-2].文獻(xiàn)[3]提出了利用阻抗和行波組合測(cè)距的算法,利用阻抗大小來(lái)確定故障區(qū)段,文獻(xiàn)[4]提出了一種利用地模分量和線模分量之間的延遲差來(lái)區(qū)分故障點(diǎn)反射波和對(duì)端母線反射波.
本文采提出了一種不受行波波速影響的單端行波測(cè)距新算法,本算法不需要計(jì)算行波波速,只需在單端檢測(cè)三個(gè)波頭的到達(dá)時(shí)間,就可以計(jì)算出故障點(diǎn)到測(cè)量端的距離,實(shí)現(xiàn)故障定位.
1.1 提升小波算法
提升小波算法的基本思想是,將現(xiàn)有的小波濾波器分解成基本的構(gòu)造模塊,通過(guò)一個(gè)預(yù)測(cè)算子來(lái)預(yù)測(cè)高頻分量,并初步預(yù)測(cè)低頻分量,然后通過(guò)更新算子,對(duì)初步預(yù)測(cè)的低頻分量進(jìn)行更新和修正,從而獲得信號(hào)的低頻分量.其提升方案一般分為三步[5].
其中:
(2)預(yù)測(cè):預(yù)測(cè)是利用偶數(shù)序列和奇數(shù)序列之間的相關(guān)性,用其中一個(gè)序列來(lái)預(yù)測(cè)另一個(gè)序列,通常是用偶數(shù)序列來(lái)預(yù)測(cè)奇數(shù)序列.預(yù)測(cè)過(guò)程如下:
(3)更新:經(jīng)過(guò)分裂步驟產(chǎn)生子集的某些整體特征(如均值)可能與原始數(shù)據(jù)并不一致,為了使原信號(hào)集的某些全局特征在其子集中繼續(xù)保持,必須進(jìn)行更新.將更新過(guò)程用算子U來(lái)代替,其過(guò)程如下:
1.2 奇異性檢測(cè)原理
小波變換具有很好的空間局部性,因此,本文利用小波變換的奇異性檢測(cè)原理來(lái)確定信號(hào)突變的位置.
由上述奇異性檢測(cè)理論可知,小波變換結(jié)果反應(yīng)信號(hào)對(duì)應(yīng)位置的變化率,小波變換的模極大值對(duì)應(yīng)該點(diǎn)具有最大的變化率[6].因此,小波變換模極大值點(diǎn)與信號(hào)的突變點(diǎn)是一一對(duì)應(yīng)的,所以可以用小波變換的模極大值來(lái)確定行波波頭達(dá)到的時(shí)刻.
2.1 相模變換
由于三相輸電線路間具有耦合作用,為了消除三相間耦合的影響,需要首先對(duì)行波分量進(jìn)行相模變換[7],將三相耦合的相分量解耦為相互獨(dú)立的0、α、β模分量,本文采用凱倫貝爾變換將三相電流解耦,則有
由于線模分量不受傳播系數(shù)、特征阻抗和大地電導(dǎo)率等因素影響,而且在各個(gè)頻率下其傳播速度基本一樣,所以本文采用線模分量來(lái)測(cè)距.
2.2 故障點(diǎn)反射波與對(duì)端母線反射波的確定
通過(guò)檢測(cè)故障點(diǎn)電流反射波與對(duì)端母線電流反射波模極大值的極性,來(lái)區(qū)分故障點(diǎn)反射波與對(duì)端母線反射波.行波在傳播的過(guò)程中,經(jīng)常會(huì)遇到結(jié)點(diǎn)(不同波阻抗的連接點(diǎn)稱(chēng)為結(jié)點(diǎn)).當(dāng)行波到達(dá)結(jié)點(diǎn)時(shí)就會(huì)在結(jié)點(diǎn)處發(fā)生折射和反射,若線路波阻抗一定,則在結(jié)點(diǎn)處的折、反射系數(shù)為:
由此可知,對(duì)于單一回路來(lái)說(shuō),母線處的反射系數(shù)ρ<0,折射系數(shù)γ>0,此時(shí),故障點(diǎn)的反射波與初始行波極性相同,對(duì)端母線上的反射波與初始行波極性相反,因此,可以利用模極大值的極性來(lái)確定故障點(diǎn)反射波與對(duì)端母線反射波[8-9].
2.3測(cè)距新算法
本文采用單端行波故障測(cè)距方法,只需在M端采集故障電流行波信號(hào),行波傳輸過(guò)程如圖1所示.
圖1 行波傳輸網(wǎng)格圖
當(dāng)線路發(fā)生故障時(shí),故障點(diǎn)就會(huì)向M、N兩端發(fā)出故障行波,設(shè)故障初始電流行波到達(dá)M端的時(shí)刻為t1,對(duì)應(yīng)的故障點(diǎn)反射波到達(dá)M端的時(shí)刻為t2,對(duì)端母線電流反射波到達(dá)M測(cè)量端的時(shí)刻為t3,L是線路的全長(zhǎng),首先利用到達(dá)M端這三個(gè)時(shí)刻計(jì)算出行波波速:
可以從上式得到
然后利用模極大值的極性來(lái)區(qū)分故障點(diǎn)反射波和對(duì)端母線反射波,若在M端檢測(cè)到是故障點(diǎn)反射波,則可以得到故障距離.若在M端檢測(cè)到是對(duì)端母線反射波,則可以得到故障距離d=.本文所提單端行波測(cè)距新算法流程圖如圖2所示.
為了驗(yàn)證本算法的正確性,采用PSCAD和MATLAB進(jìn)行線路仿真,仿真模型如圖3所示,采用簡(jiǎn)單的雙電源系統(tǒng),其中電壓為220k V,輸電線路長(zhǎng)度為L(zhǎng)=200km,采樣頻率f=1MHz.
設(shè)在輸電線路上距M端50km處發(fā)生單相接地短路故障,過(guò)渡電阻選為50Ω,通過(guò)對(duì)M端測(cè)得的電流行波經(jīng)相模變換后,再用小波變換求得突變量對(duì)應(yīng)的模極大值,從而可以準(zhǔn)確地確定初始電流行波、故障點(diǎn)反射波和對(duì)端母線反射波到達(dá)M端的時(shí)刻,如圖4所示,其對(duì)應(yīng)的時(shí)刻分別為t1=170μs, t2=505μs,t3=1173μs,由此可以得到波速
圖2 單端行波測(cè)距新算法流程圖
圖3 仿真模型圖
υ=298.954 km/ms,以及計(jì)算故障距離d1= 50.075km,實(shí)際測(cè)距誤差為75m,遠(yuǎn)小于理論測(cè)距誤差150m,仿真結(jié)果表明本文所提的單端行波故障測(cè)距算法是比較精確的.
圖4 電流線模分量和故障行波到的時(shí)刻示意圖(d=5O km)
對(duì)輸電線路上不同的故障距離進(jìn)行仿真,然后根據(jù)本文所提的算法進(jìn)行計(jì)算,計(jì)算結(jié)果見(jiàn)表1.可以看出,當(dāng)不同的距離處發(fā)生故障時(shí),經(jīng)本文的單端測(cè)距新算法計(jì)算所得到的故障距離,其最大誤差為105m,一般在75m左右,PSCAD和MATLAB仿真表明本文所提測(cè)距方法是可行的.表2為常規(guī)單端行波測(cè)距的仿真結(jié)果,取波速度υ=297.54km/ms,若在M測(cè)量端檢測(cè)到是故障點(diǎn)反射波,則利用公式來(lái)計(jì)算故障點(diǎn)到測(cè)量端的距離,若在M測(cè)量端檢測(cè)到是對(duì)端母線反射波,則利用公式d計(jì)算故障點(diǎn)測(cè)量端的距離,利用上述公式計(jì)算所得的結(jié)果,其誤差大約為200m左右,遠(yuǎn)遠(yuǎn)大于本文單端測(cè)距新算法計(jì)算結(jié)果誤差.
表1 常規(guī)單端行波測(cè)距的仿真計(jì)算表
表2 單端行波測(cè)距新算法的仿真計(jì)算表
從表1、表2可以看出常規(guī)單端行波測(cè)距需行波波速度來(lái)計(jì)算故障距離,容易受行波波速的影響,精確度不高,而本文的單端行波測(cè)距新算法將波速這個(gè)參數(shù)消去了,從而避免了行波波速的影響,使測(cè)距的精度有了明顯的提高.
本文提出了一種不受行波波速影響的新型單端行波故障測(cè)距算法,采用小波模極大值確定出行波波頭到達(dá)母線測(cè)量端的準(zhǔn)確時(shí)刻,并根據(jù)本文所提的單端行波測(cè)距新算法計(jì)算故障點(diǎn)到測(cè)量端得距離.該算法消除了行波波速的影響,減小了測(cè)距誤差,提高了測(cè)距的精確度.
[1]覃劍,彭麗萍,王和春.基于小波變換技術(shù)的輸電線路單端行波故障測(cè)距[J].電力系統(tǒng)自動(dòng)化,2005,29(19):62-65.
[2]徐丙垠,李京,陳平.現(xiàn)代行波測(cè)距技術(shù)及其應(yīng)用[J].電力系統(tǒng)自動(dòng)化,2001,25(23):62-65.
[3]鄭秀玉,丁堅(jiān)勇,黃娜.輸電線路單端故障定位的阻抗行波組合算法[J].電力系統(tǒng)保護(hù)與控制,2010,38(6):18-21.
[4]Abur A,Mangnago F H.Use of delays time between model components in wavelet based fault location[J].Electrical Power and Energy Systems,2000,22(6):397-403.
[5]周陳斌.基于提升小波的電力系統(tǒng)單端行波測(cè)距仿真[J].江蘇電機(jī)工程,2007,26(5):28-30.
[6]馬丹丹,王曉茹.基于小波模極大值的單端行波故障測(cè)距[J].電力系統(tǒng)保護(hù)與控制,2009,37(03):55-59.
[7]齊爽,翟玉成,等.單端行波測(cè)距新算法的研究[J].東北電力大學(xué)學(xué)報(bào),2011,31(1):48-51.
[8]艾斌,呂艷萍.基于小波模極大值極性的行波信號(hào)識(shí)別[J].電網(wǎng)技術(shù),2003,27(5):55-57.
[9]蔣濤,陸于平.不受波速影響的輸電線路的單端行波故障測(cè)距研究[J].電力自動(dòng)化設(shè)備,2004,24(12):29-32.
(編輯:劉寶江)
A new single termianal fault locating algorithm unaffeated by the traveling wave velocity
YANG Guan-chun1,PENG Li2
(1.School of Electrical and Electronic Engineering,Shandong University of Technology,Zibo 255091,China;2.Linyi Power Supply Company,Shandong Provincial Electric Power Company,Linyi 276000,China)
The traditional single terminal traveling wave fault location algorithm needs to consider the wave velocity,and is prone to large range error.Traditional single terminal fault location methods consider the wave velocity and tend to be inaccurate.In this paper,lifting wavelet transform and singularity detection are used to detect the moment that current wave reaches the measuring terminals.Firstly,initial current wave,fault reflected wave as well as peer bus reflected wave are utilized to calculate wave velocity.Then,based on the traditional single terminal fault location algorithm,a new single terminal fault location algorithm unaffected by traveling wave velocity is put forward by eliminating wave velocity parameters.Finally,the high accuracy and practicality are proved by the results of PSCAD and MATLAB simulation.
single terminal fault location;travelling wave;wave velocity;lifting wavelet;singularity detection
1672―6197(2013)01―0059―04
TM773
A
2012- 11- 06
楊關(guān)春,男,yangguanchun1888@163.com