吳浩然,李 澄
(國(guó)網(wǎng)無錫供電公司,江蘇無錫214061)
采用不同步數(shù)據(jù)的故障測(cè)距新原理研究
吳浩然,李 澄
(國(guó)網(wǎng)無錫供電公司,江蘇無錫214061)
為解決雙端測(cè)距中無法忽略的不同步誤差問題,提出了一種采用不同步數(shù)據(jù)的故障測(cè)距新原理。該原理對(duì)原有的雙端同步測(cè)距方程進(jìn)行不同步化處理,通過對(duì)測(cè)距方程兩邊同時(shí)取共軛復(fù)數(shù)后得到包含故障距離和不同步角2個(gè)未知量的方程組,該方程組可以轉(zhuǎn)化為一元二次方程,求解即可實(shí)現(xiàn)不同步故障測(cè)距。仿真結(jié)果表明,提出的不同步測(cè)距的方法在原理上能消去不同步角,測(cè)距精度高,基本不受過渡電阻影響,可適用于各種短路故障。
共軛復(fù)數(shù);一元二次方程;不同步;故障測(cè)距
惡劣天氣下,連接廠站的高壓線路頻繁發(fā)生短路故障,容易造成用戶失電。一般情況下,線路配置了保護(hù)裝置,可以確保系統(tǒng)穩(wěn)定和減少短路危害。從調(diào)度監(jiān)控系統(tǒng)端得到的線路保護(hù)裝置動(dòng)作情況可以獲悉具體線路故障形式以及實(shí)時(shí)電氣量數(shù)據(jù)。為了加快對(duì)失電用戶恢復(fù)送電,需要及時(shí)準(zhǔn)確地找到故障位置。以往的人工巡線工作隨著電力信息化發(fā)展逐步轉(zhuǎn)化為參照裝置提供的故障測(cè)距數(shù)據(jù)進(jìn)行有目的的巡線。因此精確的測(cè)距結(jié)果有利于迅速恢復(fù)供電系統(tǒng)。實(shí)際生產(chǎn)運(yùn)行中,連接廠、站的線路兩端數(shù)據(jù)存在不同步的問題。在故障分析的過程中,往往需要計(jì)算得出測(cè)距信息,而采用不同步的數(shù)據(jù)計(jì)算較同步的數(shù)據(jù)會(huì)產(chǎn)生相當(dāng)大的誤差,不利于得到準(zhǔn)確的故障分析結(jié)果。至此,需要研究雙端數(shù)據(jù)不同步的故障測(cè)距原理,從而改良線路保護(hù)裝置提供的測(cè)距數(shù)據(jù)。故障測(cè)距的理論和實(shí)際工程研究始終是國(guó)內(nèi)外的熱點(diǎn)。故障測(cè)距就是確定故障點(diǎn)相對(duì)于已知點(diǎn)的距離,即獲得故障位置相對(duì)于參考廠站位置的距離。諸多學(xué)者主要研究行波法和阻抗法。其中,阻抗法的故障分析原理可靠性高,精度優(yōu)良,簡(jiǎn)單易實(shí)現(xiàn),工程實(shí)施量小。近幾年,國(guó)內(nèi)外關(guān)于故障測(cè)距的研究主要集中在:(1)無需同步通信的T型輸電線路測(cè)距[1];(2)由電纜和架空線組成的混合線路測(cè)距[2,3];(3)含F(xiàn)ACTS補(bǔ)償裝置的測(cè)距方法研究[4];(4)基于行波法的全網(wǎng)廣域測(cè)距[5];(5)更具經(jīng)濟(jì)性的同桿雙回線測(cè)距[6]。
現(xiàn)今,連接于廠站的由各種電壓等級(jí)的線路兩端母線無法獲得同步的電氣量數(shù)據(jù),將會(huì)給具有明顯優(yōu)勢(shì)的雙端故障測(cè)距原理帶來不可忽略的誤差??紤]了雙端數(shù)據(jù)不同步的測(cè)距研究有故障分量法[7]、通過對(duì)復(fù)雜測(cè)距方程求導(dǎo)后得到簡(jiǎn)單線性方程[8,9]求解實(shí)現(xiàn)測(cè)距及形成非線性優(yōu)化問題來求解非線性最優(yōu)[10]方程的方法。
阻抗法故障測(cè)距主要通過構(gòu)造測(cè)距方程來實(shí)現(xiàn),其數(shù)學(xué)模型解析為簡(jiǎn)單的一元方程求解,而不同步故障測(cè)距則是對(duì)該簡(jiǎn)單方程增加了一個(gè)未知量 (不同步角),通過一些方法構(gòu)造出方程組,將測(cè)距問題轉(zhuǎn)化為線性一元方程求解的問題。本文主要采用共軛復(fù)數(shù)法,通過對(duì)測(cè)距方程等號(hào)兩邊取共軛復(fù)數(shù)的方法來構(gòu)造新的測(cè)距方程來實(shí)現(xiàn)。以下分別在集中參數(shù)阻抗模型和分布參數(shù)阻抗模型上進(jìn)行不同步測(cè)距原理的推導(dǎo)。
1.1 集中參數(shù)模型
假設(shè)線路完全換位,采用單線表示的集中參數(shù)線路故障的等值電路如圖1所示。
圖1 單線表示的集中參數(shù)線路故障等值圖
線路全長(zhǎng)為L(zhǎng),單位長(zhǎng)度阻抗Z,起始端S、終端R,故障位置F距起始端S的距離為D。在故障點(diǎn)F處的電壓為:
式(1,2)中:US為始端S電壓;IS為始端S電流;UR為終端R電壓;IR為終端R電流;UFS為由始端S推算到F點(diǎn)的電壓;UFR為由終端R推算到F點(diǎn)的電壓。
記線路的R端相對(duì)于S端偏差的不同步角δ。因此在故障點(diǎn)F處滿足:
整理式(4)得到故障距離D:
從式(5)中可以看出,US,IS,UR,IR,Z均為已知量,y是惟一的未知量。采用求解未知量y,然后將其代入故障測(cè)距解析式實(shí)現(xiàn)不同步測(cè)距。因此求解未知量即可。
對(duì)式(5)的等號(hào)左右兩邊取共軛復(fù)數(shù),得:
式(6)中:U*S,I*S,U*R,I*R,Z*分別為US,IS,UR,IR,Z的共軛復(fù)數(shù)。
注意到L*=L,D*=D,且y*=1/y,將其代入式(6)中得:
由式(6,7)得:
寫成如式(9)的一元二次方程的形式:
式(9)中:a=U*SIRZ+URI*SZ*+IRLZI*SZ*;b=U*SISZ+ USI*SZ*+URI*RZ*-U*RIRZ+LI*RZ*IRZ-IRLZI*RZ*;c=-U*RISZ-USI*RZ*-LI*RZ*I*SZ。
解式(9)得:
然后將式(10)代入式(5),求得(y1,D1)和(y2,D2)兩組解,此時(shí)需要判別偽根。
1.2 分布參數(shù)模型
輸電線路分布參數(shù)模型考慮了集中參數(shù)模型中忽略的分布電容影響,因此模型更加精確。假設(shè)長(zhǎng)距離線路完全換位,故障點(diǎn)F距離始端S距離為D,線路總長(zhǎng)度為L(zhǎng),線路發(fā)生故障如圖2所示。
圖2 單線表示的分布參數(shù)線路故障示意圖
由均勻傳輸線電壓方程得:
式(11,12)中:US,IS分別為S端電壓﹑電流;UR,IR分別為R端電壓﹑電流;UFS和UFR分別為從S端和R端推算至故障點(diǎn)的電壓;γ為傳播常數(shù);ZC為波阻抗。已知:
對(duì)式(15)兩邊取共軛復(fù)數(shù),得:
兩端電壓、電流的共軛復(fù)數(shù)形式物理量釋義同前述,且Z*C為特性阻抗ZC的共軛復(fù)數(shù);γ*C為傳播常數(shù)γC的共軛復(fù)數(shù);D*為故障距離D的共軛復(fù)數(shù)。
傳播常數(shù)γ(γ=α+jβ)在遠(yuǎn)距離、高電壓情況下,α可作為零處理,此時(shí)γ=jβ,且注意到D*=D、y*= 1/y,因此將等式(15,16)左右兩邊相乘,有e2D(γ+γ*)=1恒成立,因此:
式(17)中:該一元二次方程的系數(shù)為:a'=e-γL(UR+ZCIR)(U*S-Z*CI*S)-eγL(UR-ZCIR)(U*S+Z*CI*S);b'=eγ*L(U*R-Z*CI*R)eγL(UR-ZCIR)+(US+ZCIS)(U*S+Z*CI*S)-(US-ZCIS)(U*S-Z*CI*S)-e-γL(UR+ZCIR)e-γ*L(U*R+Z*CI*R);c'=(US-ZCIS)e-γ*L(U*R+Z*CI*R)-(US+ZCIS)eγ*L(U*R-Z*CI*R)。
解式(17)得:
最后,將由式(18)求得的代入式(15)可得到故障距離D。
記實(shí)際故障距離為Dt,則2種線路模型求得的故障距離誤差Δ:
1.3 新原理分析
實(shí)際應(yīng)用中,通常長(zhǎng)距離采用分布參數(shù)模型,而短距離輸電線路采樣集中參數(shù)模型測(cè)距即可符合現(xiàn)場(chǎng)要求。當(dāng)線路發(fā)生故障時(shí),通過求解由不同步化處理后的測(cè)距方程和經(jīng)共軛復(fù)數(shù)法構(gòu)造出的新方程組成的方程組,獲得故障位置。
按照上述測(cè)距原理,分別采用2種線路模型進(jìn)行計(jì)算,在2種線路模型上構(gòu)建一元二次方程求解均實(shí)現(xiàn)了不同步測(cè)距。構(gòu)建測(cè)距方程的過程均是對(duì)電氣量取共軛復(fù)數(shù)然后構(gòu)造了方程,且都是采用先求不同步時(shí)間,然后將其作為已知量代入故障距離函數(shù)求解??蔀樽杩狗ú煌焦收蠝y(cè)距研究提供理論基礎(chǔ)。上述原理中存在的一元二次方程求解過程,不可避免會(huì)存在偽根的問題,采用文獻(xiàn)[11]中的偽根判別法可以得到準(zhǔn)確的測(cè)距結(jié)果。
在PSCAD軟件中搭建一條長(zhǎng)為300 km的500 kV分布參數(shù)模型輸電系統(tǒng),如圖3所示。
圖3 雙端供電系統(tǒng)
集中參數(shù)模型即在長(zhǎng)線模型中忽略分布電容即可,具體系統(tǒng)參數(shù)如下:
(1)輸電線路。正序阻抗Z1=0.017 82+j0.313 9 Ω/km,正序?qū)Φ仉娙軨1=0.011 636 μF/km;零序阻抗Z0=0.295 22+j1.039 89 Ω/km,零序?qū)Φ仉娙軨0=0.007 685 6 μF/km。
(2)S和R兩側(cè)系統(tǒng)。S端正序阻抗ZS1=0.750 5+ j37.133 6 Ω,零序阻抗ZS0=0.43+j25.038 5 Ω;R端正序阻抗ZR1=18.5+j35.217 3 Ω,零序阻抗ZR0=14.2+ j24.598 7 Ω;兩端電壓 US=525∠δ kV,UR=525∠0 kV。
在PSCAD仿真軟件中,設(shè)置系統(tǒng)采樣頻率為1.8 kHz,故障時(shí)刻后60 ms以內(nèi)的數(shù)據(jù)通過差分傅式濾波,采用文中提出的測(cè)距原理進(jìn)行Matlab編程計(jì)算,即可實(shí)現(xiàn)測(cè)距。
線路發(fā)生A-G(單相接地)故障,當(dāng)采用2種不同的線路模型時(shí),不同步角δ變化時(shí)的測(cè)距結(jié)果如表1所示。當(dāng)分別為35°,100°,135°時(shí),分部參數(shù)模型測(cè)距誤差范圍為[0.19%,0.32%],集中參數(shù)模型測(cè)距誤差范圍為[0.22%,0.44%]。明顯可以看出,當(dāng)線路長(zhǎng)度大于100 km后,集中參數(shù)模型的測(cè)距誤差開始增大,變化時(shí)的測(cè)距結(jié)果良好。
表2為過渡電阻Rg對(duì)2種模型的兩相接地故障測(cè)距結(jié)果。分布參數(shù)模型的測(cè)距結(jié)果隨著過渡電阻的增大(Rg分別取80 Ω,100 Ω,200 Ω)幾乎不變,其測(cè)距結(jié)果的誤差范圍為[0.05%,0.26%],集中參數(shù)模型的測(cè)距誤差范圍為[0.08%,0.52%],綜合表中數(shù)據(jù)可以得出2種模型的不同步數(shù)據(jù)的測(cè)距原理可以認(rèn)為均不受過渡電阻影響。
表1 不同線路模型中δ對(duì)A-G故障的測(cè)距影響
表2 不同線路模型中Rg對(duì)AB-G故障的測(cè)距影響
表3示出了不同故障類型對(duì)2種模型的測(cè)距影響。在距離線路一側(cè)100 km位置故障時(shí)的測(cè)距結(jié)果中,4種短路故障情況下,集中參數(shù)模型測(cè)距最大誤差為0.52%,而分部參數(shù)模型的測(cè)距最大誤差為0.27%,2種方法都不受故障類型的影響。
表3 不同模型及不同故障類型的測(cè)距結(jié)果
本文提出了雙端數(shù)據(jù)不同步的故障測(cè)距原理,主要解決了不同步故障測(cè)距誤差的問題。該原理分別基于線路的集中參數(shù)和分布參數(shù)模型可推導(dǎo)出一元二次方程形式的測(cè)距方程,通過對(duì)測(cè)距方程取共軛復(fù)數(shù)構(gòu)建出的簡(jiǎn)單線性方程求解容易,能避免求解復(fù)雜的非線性方程的繁瑣與誤差擴(kuò)大化,僅需對(duì)有可能存在的偽根進(jìn)行判別,在兩種線路模型上均適用。不存在非線性方程的最優(yōu)求解不準(zhǔn)確問題,適用于4種基本的短路故障類型,工程應(yīng)用價(jià)值優(yōu)良。
參考文獻(xiàn):
[1]KHODADADI M,SHAHRTASH S M.A New Noncommunication-Based Protection Scheme for Three-Terminal Transmission Lines Employing Mathematical Morphology-Based Filters.IEEE Transactions On Power Delivery,2013,28(1)∶347-356.
[2]LIU C W,LIN T C,YU C S,et al.A Fault Location Technique for Two-Terminal Multisection Compound Transmission Lines Using Synchronized Phasor Measurements[J].IEEE Transactions On Smart Grid,2012,03(4)∶113-121.
[3]王 濤,李 澄.一種混合輸電線路故障測(cè)距新方法[J].江蘇電機(jī)工程,2015,34(5):17-20,24.
[4]AHSAEE M G,SADEH J.A Novel Fault-Location Algorithm for Long Transmission Lines Compensated by Series FACTS Devices [J].IEEE Transactions On Power Delivery,2011,26(4)∶2299-2308.
[5]KORKALI M,ABUR A.Optimal Deployment of Wide-Area Synchronized Measurements for Fault-Location Observability[J]. IEEE Transactions On Power Systems,2013,28(1)∶482-489.
[6]林富洪,王增平.采用同向正序基頻分量的雙回線故障測(cè)距原理[J].中國(guó)電機(jī)工程學(xué)報(bào),2011,31(4):93-98.
[7]施世鴻,何奔騰.基于分布參數(shù)模型的雙端非同步故障測(cè)距算法[J].電網(wǎng)技術(shù),2008,32(9):84-88.
[8]YU C S,CHANG L R,CHO J R.New Fault Impedance Computations for Unsynchronized Two-Terminal Fault-Location Computations[J].IEEE Transactions On Power Delivery,2011,26(4)∶2879-2881.
[9]DOUGLAS A G V,DIOGO B O,ADRIANO C L.A Closed-Form Solution for Transmission-Line Fault Location Without the Need of Terminal Synchronization or Line Parameters[J].IEEE Transactions On Power Delivery,2013,28(2)∶1238-1239.
[10]蘇進(jìn)喜,羅承沐,解子鳳,等.利用雙端電氣量的高壓長(zhǎng)線故障測(cè)距算法[J].清華大學(xué)學(xué)報(bào),2000,40(7):27-30.
[11]李 澄,王 鍇.雙端測(cè)量不同步的非對(duì)稱故障測(cè)距新方法[J].電測(cè)與儀表,2014,51(13):55-59.
Research on Novel Fault Location Principle with Asynchronous Data
WU Haoran,LI Cheng
(State Grid Wuxi Power Supply Company,Wuxi 214061,China)
A novel fault location principle with asynchronous data is proposed to solve the asynchronous problem of fault location at both ends which can't be ignored.The primary dual-terminal synchronous fault location equation is asynchronized in the principle by taking conjugate complex number at both sides of the equation simultaneously,which transforms the equation into a set of equations with two unknown variables of fault location and asynchronous angle.The set of equations can be transformed into a quadratic equation with one variable which can be solved to reach the goal of asynchronous fault location.Simulations show that the proposed asynchronous fault location method can eliminate the asynchronous angle and has high precision which is hardly affected by transition resistance,so it is adequate for all short circuit faults.
conjugate complex number;quadratic equation with one variable;asynchronous;fault location
TM77
A
1009-0665(2016)06-0073-04
吳浩然(1964),男,江蘇鎮(zhèn)江人,高級(jí)工程師,從事電力系統(tǒng)管理專業(yè)工作;
李 澄(1987),男,江蘇江陰人,工程師,從事電網(wǎng)調(diào)度運(yùn)行、繼電保護(hù)與故障測(cè)距方面的工作。?
2016-08-08;
2016-09-13