劉小磊 李冰心 趙峰

摘 要：柔性高壓直流輸電系統(tǒng)大多采用地下電纜作為輸電線路，由于地下電纜行波波頭衰減嚴(yán)重，導(dǎo)致在長距離故障測距時很難準(zhǔn)確測出故障點(diǎn)反射波到達(dá)測量點(diǎn)的時刻。為了實(shí)現(xiàn)精確的故障測距，將奇異值（SVD）分解算法應(yīng)用到固有頻率法與行波法相結(jié)合的組合算法中。首先利用故障行波的固有頻率計(jì)算出故障點(diǎn)的大概位置，算出故障點(diǎn)反射波分別到達(dá)兩端母線的時間范圍，然后對故障行波信號進(jìn)行奇異值分解，得到SVD第二分量奇異值檢測結(jié)果，依據(jù)固有頻率法算出的時間范圍確定反射波頭到達(dá)兩端母線的準(zhǔn)確時刻。最后將初始行波和故障點(diǎn)反射波分別到達(dá)兩端母線的時間代入消除波速影響雙端測距公式，計(jì)算出故障點(diǎn)的準(zhǔn)確位置。PSCAD和MATLAB仿真證明該方法測距精度較高。

關(guān)鍵詞：柔性高壓直流；行波法；固有頻率；奇異值分解；雙端故障測距；波速

DOI：10.11907/rjdk.172745

中圖分類號：TP319

文獻(xiàn)標(biāo)識碼：A 文章編號：1672-7800（2018）006-0173-05

Abstract：In most cases underground cables are utilized as the transmission lines in flexible high voltage direct current transmission （VSC-HVDC） system， but the serious attenuation of the underground cable wave head leads to great difficulties in accurately measuring the arrival time of fault point reflected waves. Therefore a new algorithm of two terminal fault location combining natural frequency and traveling wave method based on the time-frequency characteristics of fault traveling wave is proposed. Firstly， the natural frequency of the fault traveling wave is used to calculate the approximate location of the fault point， so that the time range of the reflected wave between the bus at both ends can be calculated. Then the singular values of the fault traveling wave are decomposed， and the singular points of the fault signals are detected by the pulse of second SVD components， according to the corresponding time range； the accurate time of the wave head reaching the bus at both ends is determined. VSC-HVDC transmission lines are on on both sides of the combined capacitor， so the fault traveling wave in the bus has total reflection. At the same time， we put forward a algorithm eliminating the influence of wave velocity of two terminal fault location by the initial traveling wave and time equations that reflect when wave reaches the fault point respectively at both ends of bus time equations are listed， and then calculate the distance and determine the exact location of the fault point. PSCAD and MATLAB simulations show that the proposed algorithm has higher ranging accuracy and can solve the problem of the recognition of the reflected wave head.

Key Words：flexible high voltage direct current； traveling wave method； natural frequency； SVD； two terminal fault location； wave velocity

0 引言

柔性高壓直流輸電VSC-HVDC（voltage source converter-high voltage direct current transmission）是一種以電壓源換流器和脈寬調(diào)制等技術(shù)為基礎(chǔ)的新一代直流輸電技術(shù)，具有較好的可控性和靈活性，在新能源發(fā)電并網(wǎng)、孤島供電等領(lǐng)域有很好的應(yīng)用前景[1]。直流輸電線路大多距離較長，故障發(fā)生率較高，且故障發(fā)生后故障點(diǎn)的查找比較困難。因此，研究精確可靠的故障測距技術(shù)，對于保證電力系統(tǒng)安全運(yùn)行，提高系統(tǒng)可靠性和經(jīng)濟(jì)性有著極為重要的作用。直流故障測距裝置主要采用行波法進(jìn)行測距，行波測距法響應(yīng)速度快、測距精度高，理論上不受過渡電阻、故障類型和相角等因素影響[2-6]。行波測距分為單端測距法和雙端測距法[7]。單端測距法是指根據(jù)故障初始行波從故障點(diǎn)到測量點(diǎn)往返一次所用的時間和波速計(jì)算故障點(diǎn)位置；雙端測距法是指利用故障初始行波分別到達(dá)兩端測量點(diǎn)的時間差和波速計(jì)算故障點(diǎn)位置。

行波法測距需要對故障點(diǎn)反射波頭進(jìn)行準(zhǔn)確識別，這就降低了測距的可靠性[8-9]。而基于固有頻率的測距算法則不受波頭識別困難的影響，只需在任意故障暫態(tài)數(shù)據(jù)提取固有頻率就能進(jìn)行測距。該算法測距的精度明顯低于行波法，但可靠性較高、穩(wěn)定性強(qiáng)，可以對故障點(diǎn)位置進(jìn)行粗略估算。本文將奇異值（SVD）分解算法用到固有頻率法與行波法相結(jié)合的組合算法中，利用兩種測距方法優(yōu)勢互補(bǔ)[10-11]，將奇異值分解（SVD）算法引入使得測距計(jì)算更為精確。該算法首先利用固有頻率法測出故障點(diǎn)的粗略值，計(jì)算出故障點(diǎn)反射波到達(dá)兩端測量點(diǎn)的時間范圍，從而實(shí)現(xiàn)對反射波頭精確時刻的主動識別，很大程度上提高了測距的可靠性，而奇異值（SVD）分解具有準(zhǔn)確識別故障信號的優(yōu)點(diǎn)，能得到更為精確的時間參數(shù)。經(jīng)過PSCAD仿真和MATLAB數(shù)據(jù)處理，證明該測距算法能實(shí)現(xiàn)高壓直流輸電線路可靠、精確的故障測距。

1 固有頻率法測距

輸電線路發(fā)生故障時，形成故障暫態(tài)行波，在頻域上表現(xiàn)為一系列特定頻率諧波的形式，稱為固有頻率，固有頻率值與故障距離呈函數(shù)關(guān)系。文獻(xiàn)[12]研究了故障行波頻率、邊界條件和故障距離之間的關(guān)系：

輸電線路發(fā)生金屬性接地短路時，故障行波在短路點(diǎn)的反射系數(shù)為-1，也就是說α-F為π；直流輸電線路母線處裝有平波電抗器，一般情況下故障暫態(tài)行波頻率較高，平波電抗器阻抗較大，可以近似認(rèn)為是開路，反射系數(shù)為1， α-F取π。在故障行波的頻譜中，第一個峰值對應(yīng)的頻率最大，為主頻，利用主頻對故障點(diǎn)進(jìn)行測距最為簡單易行，這時對應(yīng)的測距公式為：

固有頻率法測距是根據(jù)故障行波的固有頻率進(jìn)行測距的，不用識別故障點(diǎn)反射波頭。用向后預(yù)測的prony算法提取故障電流行波的固有頻率。故障行波波頭平緩，在時域范圍內(nèi)對行波波頭的標(biāo)定比較困難。但是經(jīng)過時頻轉(zhuǎn)變后，可以有效提取固有頻率值，實(shí)現(xiàn)故障測距。該方法抗高阻能力強(qiáng)，受噪聲干擾小，可靠性高，缺點(diǎn)是測距精度不高。

2 雙端故障測距算法

行波在線路中傳播時，在母線、故障點(diǎn)和線路中的不連續(xù)點(diǎn)處都會發(fā)生反射和折射，而對于VSC-HVDC而言，系統(tǒng)兩端并聯(lián)有大電容，行波會在母線處發(fā)生全反射。傳統(tǒng)雙端測距法只是根據(jù)行波到達(dá)兩端母線的時間差進(jìn)行故障定位，完全沒有考慮波速不確定性對測距精度帶來的影響。行波在不同的電壓等級和線路狀態(tài)以及線路參數(shù)下傳播速度是不同的，一般情況下，波速的取值范圍是0.936c（11kV）～0.987c（500kV），這給測距帶來了較大誤差。在傳統(tǒng)行波法測距的基礎(chǔ)上，利用故障初始行波，將故障點(diǎn)的反射波分別到達(dá)測量點(diǎn)的時間參數(shù)和波速、故障距離的關(guān)系列出方程組，求解得到測距公式，這種算法能夠消除波速影響[13]。

VSC-HVDC輸電線路發(fā)生接地短路的故障最多，本文對此進(jìn)行研究。圖1為輸電線路模型，線路在F點(diǎn)處發(fā)生接地短路，M端為整流側(cè)，N端為逆變側(cè)。由于VSC-HVDC輸電線路兩端并聯(lián)大電容，所以故障行波在母線M、N處發(fā)生全反射。

故障行波傳輸網(wǎng)絡(luò)如圖2所示，t-1為故障初始行波到達(dá)M端測量點(diǎn)的時間，t-2為故障行波經(jīng)M端母線和故障點(diǎn)反射后再次到達(dá)M端測量點(diǎn)的時間，t-3為故障初始行波到達(dá)N端測量點(diǎn)的時間，t-4為故障行波經(jīng)N端母線和故障點(diǎn)反射后再次到達(dá)N端測量點(diǎn)的時間。

設(shè)MN之間的長度為L，故障點(diǎn)F到M端、N端母線的距離分別為d-1和d-2，波速為v?？闪谐鍪剑?）：

測距公式（4）是關(guān)于時間和距離的式子，消除了波速不確定給測距帶來較大誤差的影響，而t-1、t-2、t-3、t-4的測量不需要時間同步，消除了時間不同步帶來的誤差，因此不需要GPS雙側(cè)時鐘同步，節(jié)約了測距裝置成本。若兩端故障點(diǎn)反射波沒有檢測到，就只能用傳統(tǒng)雙端測距進(jìn)行故障定位。

3 SVD分解與奇異性檢測

4 雙端故障測距組合算法

綜上所述，固有頻率法是在頻域上對故障點(diǎn)進(jìn)行測距，行波法是在時域上對故障點(diǎn)進(jìn)行測距，行波法測距精度高，固有頻率法可靠性好，兩種測距方法可以優(yōu)勢互補(bǔ)。本文在行波法和固有頻率法結(jié)合的基礎(chǔ)上[10-11]，用奇異值分解（SVD）算法對故障行波信號進(jìn)行處理，得到的測距結(jié)果更為精確。

4.1 固有頻率法初測

設(shè)利用固有頻率法初步測量得到的故障距離為x，大量仿真實(shí)驗(yàn)表明，在任何形式的故障下，利用固有頻率法測距得到的結(jié)果誤差最大不超過10%，則實(shí)際故障距離范圍為（0.9x，1.1x），代入傳統(tǒng)測距公式d=v（t-2-t-1）2，可以算出故障點(diǎn)反射波到達(dá)M端母線的時間t-2的范圍為t-1+1.8xv，t-1+2.2xv。同理，可得到t-4的范圍為t-3+1.8（L-x）v，t-3+2.2（L-x）v。

4.2 波頭識別

故障行波會在線路上阻抗不連續(xù)點(diǎn)處發(fā)生突變，對故障電壓行波信號進(jìn)行奇異值分解，利用SVD第二分量脈沖極大值或極小值確定故障初始行波和故障點(diǎn)反射波分別到達(dá)兩端測量點(diǎn)的位置。

4.3 雙端測距

首先根據(jù)SVD第二分量脈沖的極值確定故障初始行波分別到達(dá)兩端測量點(diǎn)的時刻t-1和t-3，然后用固有頻率法算出故障點(diǎn)反射波到達(dá)兩端母線的時間范圍，識別出故障點(diǎn)反射波經(jīng)SVD第二分量分解后對應(yīng)的準(zhǔn)確時刻t-2和t-4，將t-1、t-2、t-3、t-4帶入公式（5）計(jì)算出故障點(diǎn)的確切位置。

5 仿真驗(yàn)證

為了驗(yàn)證本文故障測距新方法的定位精度，在PSCAD中搭建如圖3所示的VSC-HVDC模型，額定電壓為±110kV，容量為75MW，線路長度為200km，采樣頻率為1MHz，兩端分別并聯(lián) 1000μF的大電容，采用PSCAD對輸電系統(tǒng)進(jìn)行仿真，并用MATLAB進(jìn)行數(shù)據(jù)處理。

設(shè)置在距離M端20km處發(fā)生接地短路故障，過度電阻為100Ω時，經(jīng)PSCAD仿真后得到的M和N端的故障電壓波形如圖4所示。

首先用向后預(yù)測的Prony算法提取故障行波的固有頻率，其中主頻f-0=7 326Hz，在電壓等級為110kV下v-0取0.936c，代入公式（3）得x=19.165km，從而得出t-2的范圍為（t-1+122.853，t-1+150.153），同理可求出t-4的范圍為（t-3+1159.199，t-3+1416.798）。對故障電壓行波進(jìn)行奇異值分解，得到SVD第二分量奇異值檢測結(jié)果如圖5、圖6所示。根據(jù)SVD檢測結(jié)果提取突變點(diǎn)極值，可以看出，利用SVD分解到的第二分量信號能夠準(zhǔn)確測出故障初始行波到達(dá)母線的時刻t-1和t-3，t-1=70μs，t-2=204μs。算出t-3和t-4范圍，從而確定其對應(yīng)的突變點(diǎn)精確時刻： t-3=605μs，t-4=1809μs，代入公式（4）得d-1=20.029 9km，誤差為0.029 9km，可知測距精度較高。

在不同故障距離和過度電阻下的測距結(jié)果如表1所示。

用傳統(tǒng)測距方法對故障點(diǎn)進(jìn)行計(jì)算，與新方法對比，測距結(jié)果如表2所示。

從表1、表2測距結(jié)果可以看出，新測距方法的測距結(jié)果比傳統(tǒng)測距方法的測距誤差要小很多，而且消除了波速影響，具有波頭識別優(yōu)勢，測距精度更高。

6 結(jié)語

本文基于VSC-HVDC線路故障測距的可靠性和精確性，將奇異值（SVD）分解算法用到固有頻率法與行波法相結(jié)合的組合算法中。行波法和固有頻率法各有優(yōu)點(diǎn)，將兩種測距方法優(yōu)勢互補(bǔ)，引入奇異值分解（SVD）算法使得測距計(jì)算更為精確。

對于故障信號的處理，利用故障電壓行波信號分量構(gòu)造Hankel矩陣，對其進(jìn)行奇異值分解得到SVD第二分量信號，可以準(zhǔn)確測量出信號的奇異點(diǎn)，并用SVD脈沖提取出奇異點(diǎn)極值，得到精確的時間信息。根據(jù)固有頻率法測出的故障點(diǎn)位置的模糊值計(jì)算出故障點(diǎn)反射波分別到達(dá)兩端測量點(diǎn)的時間范圍，識別出故障點(diǎn)反射波到達(dá)測量點(diǎn)的準(zhǔn)確時刻，將時間參數(shù)代入測距公式得到測距結(jié)果。經(jīng)仿真驗(yàn)證，該方法有效消除了波速的影響，提高了測距精度。

參考文獻(xiàn)：

[1] 湯廣福.基于電壓源換流器的高壓直流輸電技術(shù)[M].北京：中國電力出版社，2010.

[2] 葛耀中.新型繼電保護(hù)和故障測距的原理與技術(shù)[M].西安：西安交通大學(xué)出版社，2007.

[3] 李澤文，花歡歡，鄧豐，等.基于廣域行波信息的行波保護(hù)[J].中國電機(jī)工程學(xué)報(bào)，2014，34（34）：6238-6245.

[4] 束洪春，鄔乾晉，張廣斌，等.基于神經(jīng)網(wǎng)絡(luò)的單端行波故障測距方法[J].中國電機(jī)工程學(xué)報(bào)，2011，31（4）：85-92.

[5] 蔣洪濤，黃純，楊帥雄，等.基于小波暫態(tài)能量的超高壓輸電線路暫態(tài)保護(hù)新方法[J].電網(wǎng)技術(shù)，2012，36（8）：248-253.

[6] 韓志錕，饒曙勇，姜玉山，等.基于能量比法的輸電線路行波故障測距[J].電網(wǎng)技術(shù)，2011，35（3）：216-236.

[7] 宋國兵，蔡新雷，高淑萍，等.高壓直流輸電線路故障定位研究綜述[J].電力系統(tǒng)保護(hù)與控制，2012，40（5）：133-137，147.

[8] 徐敏，蔡澤祥，劉永浩，等.基于寬頻信息的高壓直流輸電線路行波故障測距方法[J].電工技術(shù)學(xué)報(bào)，2013，28（1）： 259-265.

[9] 朱永利，范新橋，尹金良.基于三點(diǎn)電流測量的輸電線路行波故障定位新方法[J].電工技術(shù)學(xué)報(bào)，2012，27（3）：260-268.

[10] 李博雅，楊耀，楊立紅.高壓直流輸電線路單端故障測距組合算法[J].電力系統(tǒng)保護(hù)與控制，2014，42（3）：116-121.

[11] 齊建平.內(nèi)蒙古500kV輸電線路行波故障測距算法研究[D].北京：華北電力大學(xué)，2013.

[12] 武曉，何正友，彭少博，等.基于行波固有頻率的特高壓直流輸電線路縱聯(lián)保護(hù)方法[J].電力系統(tǒng)保護(hù)與控制，2013，41（11）： 67-73.

[13] 向強(qiáng)銘，王茜，陳靖秋，等.基于奇異值分解理論的雙端行波故障測距的研究[J].電力系統(tǒng)保護(hù)與控制，2014，42（12）：14-18.

[14] NANAYAKKARA O. Line Fault Location in Emerging HVDC Transmission Systems [D]. Manitoba： University of Manitoba ， 2014.

[15] 趙學(xué)智，葉邦彥.基于二分遞推SVD的信號奇異性位置精確檢測[J].電子學(xué)報(bào)，2012，40（1）：53-59.

[16] 徐巖，張銳，霍福廣，等.基于MAS信息融合的電網(wǎng)故障診斷[J].電力科學(xué)與工程，2012，28（9）：24-28.

（責(zé)任編輯：杜能鋼）