陳仕龍，李朋松，莊啟康，畢貴紅，高敬業(yè)，趙四洪

（昆明理工大學(xué) 電力工程學(xué)院，云南 昆明 650500）

0 引言

近年來，電氣化高速鐵路在交通運輸中所占比例越來越大[1]。

傳統(tǒng)牽引供電系統(tǒng)存在的電流分相、電流負(fù)序等問題，對軌道交通運輸?shù)母咚僦剌d發(fā)展需求造成了很大的制約[2,3]。

電力電子技術(shù)的日漸成熟，給牽引供電系統(tǒng)快速穩(wěn)定發(fā)展創(chuàng)造了條件。比如，依托電力電子器件構(gòu)建的系統(tǒng)，可取消電流分相環(huán)節(jié)，提高系統(tǒng)運行的穩(wěn)定性[4]。

貫通式同相牽引供電系統(tǒng)因在提高電能質(zhì)量、系統(tǒng)穩(wěn)定性等方面具有明顯優(yōu)勢，所以目前已成為研究及應(yīng)用熱點[5]。

受沿線環(huán)境以及機車負(fù)荷沖擊等因素影響，牽引網(wǎng)線路故障可能性較高，亟需探討可靠的牽引網(wǎng)線路保護方案[6]。

現(xiàn)階段，牽引網(wǎng)的保護主要為距離保護。文獻(xiàn)[7]利用供電臂兩相之間的電壓以及諧波電流提出了饋線距離保護和電流速斷保護。文獻(xiàn)[8]通過串入高頻阻波器識別故障行波波頭，利用多測點測距方程進(jìn)行故障測距。該方法測距精度高，但未考慮串入阻波器對牽引供電系統(tǒng)的影響。文獻(xiàn)[9]分析了牽引網(wǎng)故障行波在不同頻段、不同路徑的傳播特性，通過Morlet小波變換對故障信號進(jìn)行處理，進(jìn)而得到故障測距方案。

對于同相牽引供電系統(tǒng)而言，現(xiàn)有研究主要集中在系統(tǒng)特性方面——針對同相牽引供電系統(tǒng)的控制策略、系統(tǒng)穩(wěn)定性以及提高電能質(zhì)量等問題的研究較為深入，而對牽引網(wǎng)保護方面的研究還相對較少[10-12]。文獻(xiàn)[13]利用高壓直流輸電線路正方向區(qū)內(nèi)外故障暫態(tài)功率的頻帶差異性，構(gòu)造了一種單端量保護；該方法為牽引網(wǎng)保護提供了研究思路。

同相牽引供電系統(tǒng)獨特的供電方式，導(dǎo)致其系統(tǒng)結(jié)構(gòu)特性與傳統(tǒng)牽引供電系統(tǒng)存在明顯差異，故需進(jìn)一步探討與之相適應(yīng)的保護方案。

本文以貫通式同相牽引供電系統(tǒng)為研究對象，提出了一種利用小波高、低頻暫態(tài)能量比值來區(qū)分牽引網(wǎng)線路區(qū)內(nèi)、外故障的單端方向保護方案。

1 系統(tǒng)結(jié)構(gòu)及邊界頻率特性分析

1.1 貫通式同相牽引供電系統(tǒng)結(jié)構(gòu)

貫通式同相牽引直接供電系統(tǒng)基本結(jié)構(gòu)如圖1所示（電流正方向與機車運動方向一致）。

圖1 牽引供電系統(tǒng)拓?fù)浣Y(jié)構(gòu)圖Fig. 1 Topological structure diagram of traction power supply system

圖1中，牽引變電所是供電系統(tǒng)主要組成部分，其使用的是PWM三相整流–單相PWM逆變對稱變換與雙閉環(huán)PI控制技術(shù)。牽引變電所的作用是：在消除公共電網(wǎng)與牽引供電系統(tǒng)之間相互干擾的同時，實現(xiàn)三相交流–直流–單相交流變換，使?fàn)恳W(wǎng)側(cè)電壓具有相同幅值和相位，實現(xiàn)無電分相環(huán)節(jié)情況下全線貫通供電。

變電所出口處并聯(lián)的 LC濾波器，能有效濾除電流諧波、提高電能質(zhì)量。

1.2 牽引網(wǎng)邊界頻率特性

牽引變電所出口處的 LC濾波器在起到濾除諧波作用的同時，其并聯(lián)電容和一段牽引網(wǎng)線路還構(gòu)成牽引網(wǎng)物理邊界，如圖2所示。

圖2 牽引網(wǎng)邊界等效模型Fig. 2 Equivalent model of traction network boundary

圖2中ZL為牽引網(wǎng)一段線路的等效阻抗，ZC為牽引變電所供電段區(qū)外所并聯(lián)電容的等效阻抗。

根據(jù)線路等效參數(shù)[14]和文獻(xiàn)邊界定義[15]，可計算牽引網(wǎng)邊界傳遞函數(shù)G(jω)。

定義邊界傳遞函數(shù)為U2/U1，即：

由此可得直接供電系統(tǒng)牽引網(wǎng)邊界幅頻特性G(jω)，如圖3所示。

圖3 牽引網(wǎng)邊界幅頻特性Fig. 3 Amplitude frequency characteristics of traction network boundary

由圖3可知：對于頻率在0～550 Hz頻段的信號，其G(jω)≥1；當(dāng)頻率f≥800 Hz時，其幅值遠(yuǎn)小于1，且隨頻率的增加持續(xù)呈下降趨勢?？梢姡瑺恳冸娝╇姸螀^(qū)外所并聯(lián)的電容與牽引網(wǎng)的一部分線路共同構(gòu)成的線路邊界，是信號高頻部分產(chǎn)生較大幅值衰減的主要影響因素。

2 牽引網(wǎng)故障暫態(tài)過程分析

2.1 反方向故障特征分析

正向區(qū)內(nèi)線路末端發(fā)生非金屬性接地故障產(chǎn)生的高頻分量，將會受到線路波阻抗及過渡電阻2方面的影響；而反向區(qū)外靠近線路保護安裝處，故障高頻分量只受到線路邊界衰減單方面影響，所以：有可能出現(xiàn)正向區(qū)內(nèi)線路末端故障保護裝置檢測到的高頻分量小于反向區(qū)外近保護安裝處故障產(chǎn)生的高頻分量；這將造成保護誤動作。因此，單靠邊界元件，即依靠高頻分量大小區(qū)分區(qū)內(nèi)、外故障將遇到一定困難，需要在保護安裝處增加方向元件來有效區(qū)分正方向與反方向故障，從而消除保護死區(qū)。

考慮到實際工程應(yīng)用問題，本文采用暫態(tài)功率能量極性來區(qū)分正反方向故障[16]。

規(guī)定電流正方向為母線流向線路，從而可以得到以下結(jié)論：當(dāng)發(fā)生反方向故障時，保護裝置測得的暫態(tài)電流Δi與暫態(tài)電壓Δu突變方向相同；當(dāng)發(fā)生正方向故障時，保護裝置測得的暫態(tài)電流Δi與暫態(tài)電壓Δu突變方向相反。

定義暫態(tài)功率ΔP=ΔiΔu。經(jīng)計算，得到線路不同位置故障暫態(tài)功率ΔP，如圖4所示。

圖4 正反向故障暫態(tài)功率變化曲線Fig. 4 Transient power change curve of forward and reverse faults

由圖4可以看出，暫態(tài)電流、電壓突變方向的異同造成了暫態(tài)功率幅值變化的不同，故可利用暫態(tài)功率的極性判別故障發(fā)生的方向。

2.2 正方向區(qū)內(nèi)外故障特征分析

2.2.1 高頻暫態(tài)量提取

多尺度小波分解能夠很好地表征故障信號高、低頻帶的信息：通過正交小波基對故障信號進(jìn)行逐層的快速分解，每一層分解的結(jié)果都包含了一組高頻段信息與低頻段信息，即信號的細(xì)節(jié)部分與近似部分。該方法在頻域分析方面具有優(yōu)良的性能，廣泛應(yīng)用于非平穩(wěn)信號的處理。分解過程可表示為：

設(shè)J為信號分解的最大尺度，第j尺度下的逼近信號表示為Aj，細(xì)節(jié)信號表示為Dj。選取合適的信號采樣頻率，利用高通濾波器h對逼近信號Aj進(jìn)行相關(guān)濾波處理后，采用偶次降采樣即可得到j(luò)+1尺度下的細(xì)節(jié)系數(shù)Dj+1。保留j+1尺度下的細(xì)節(jié)系數(shù)Dj+1后再對逼近信號Aj+1進(jìn)行分解，逐層分解后可得到原始信號在各個頻帶上的細(xì)節(jié)系數(shù)[17,18]。分解過程如圖5所示。

圖5 多尺度小波分解示意圖Fig. 5 Schematic diagram of multiscale wavelet decomposition

定義第j個尺度下信號的小波能量為：

式中：Ej代表j尺度下的小波能量；Dj(n)為小波分解系數(shù)；N為采樣點個數(shù)。

本文采用小波多尺度對故障電流信號的線模分量進(jìn)行分解，以得到信號的多頻帶能量信息?？紤]需要利用信號的高頻分量進(jìn)行分析，將采樣頻率設(shè)置為50 kHz，選取db5小波為小波基，多尺度分解6層。每個尺度對應(yīng)頻帶如表1所示。

表1 不同分解尺度所對應(yīng)的頻帶Tab. 1 Frequency bands corresponding to different decomposition scales kHz

據(jù)表1可知，D1尺度下的小波能量為故障電流信號的高頻能量，D6尺度下的小波能量為故障電流信號的低頻能量。

2.2.2 正方向區(qū)內(nèi)外故障特征

如圖1（a）所示，F(xiàn)1、F2、F3分別表示牽引網(wǎng)故障發(fā)生的不同位置。根據(jù)疊加原理，牽引網(wǎng)線路區(qū)內(nèi)F2點發(fā)生接地故障以及區(qū)外F3點發(fā)生接地故障后的附加網(wǎng)絡(luò)，如圖6、圖7所示。

圖6 牽引網(wǎng)區(qū)內(nèi)故障附加網(wǎng)絡(luò)Fig. 6 Fault additional network in traction network area

圖7 牽引網(wǎng)區(qū)外故障附加網(wǎng)絡(luò)Fig. 7 Fault additional network outside the traction network area

圖6—7中：ZL1、ZL2、ZC分別為牽引變出口處一部分線路的等效阻抗、牽引網(wǎng)線路剩余部分的等效阻抗、牽引變供電段區(qū)外所并聯(lián)電容的等效阻抗。

區(qū)內(nèi)故障時，故障電流行波IL只經(jīng)過線路流向檢測點，故障信號的高頻能量E1和低頻能量E6未經(jīng)過邊界的衰減，其幅值應(yīng)較大。

區(qū)外故障時，故障電流行波IL經(jīng)過線路、線路邊界流向檢測點，故障信號的高頻能量E1和低頻能量E6經(jīng)過線路和邊界的衰減，其幅值應(yīng)小于區(qū)內(nèi)故障信號高、低頻能量幅值；但對于故障信號的低頻能量而言，其衰減程度遠(yuǎn)不及高頻能量。因此，利用故障信號的高、低頻暫態(tài)能量比值，與保護裝置預(yù)先設(shè)定好的整定值進(jìn)行數(shù)值大小的對比，就可以判斷故障發(fā)生的位置。

當(dāng)牽引網(wǎng)線路發(fā)生接地故障時，對故障電流信號進(jìn)行多尺度小波分解，可得故障信號高頻、低頻暫態(tài)能量譜。

以典型故障——金屬性接地故障與高阻性接地故障為例，能量譜如圖8、圖9所示。

圖8 區(qū)內(nèi)外中點經(jīng)0.1 Ω過渡電阻接地故障時電流暫態(tài)能量譜Fig. 8 Current transient energy spectrum in case of neutral point grounding fault through 0.1 Ω transition resistance inside and outside the area

圖9 區(qū)內(nèi)外中點經(jīng)300 Ω過渡電阻接地故障時電流暫態(tài)能量譜Fig. 9 Current transient energy spectrum in case of grounding fault with 300 Ω transition resistance at the midpoint inside and outside the area

由圖8、圖9可知，無論對于金屬性接地故障還是高阻性接地故障，牽引網(wǎng)邊界都對故障信號的高頻部分具有強烈衰減作用。此特性使得區(qū)內(nèi)、外故障時故障信號的高頻能量存在明顯差異，但對于低頻能量而言，其差異相對較小。

通過上述分析可得出結(jié)論：利用小波多尺度分解提取故障電流的高、低頻暫態(tài)能量，并用其能量比值區(qū)分區(qū)內(nèi)、外故障具有可行性。

3 保護方案

當(dāng)同相牽引供電系統(tǒng)正常運行時，其電流波形為標(biāo)準(zhǔn)正弦波。發(fā)生故障后，系統(tǒng)失衡，電流幅值急劇上升，故可采用電流的變化率作為故障啟動判據(jù)。

3.1 方向判據(jù)

選取數(shù)據(jù)窗內(nèi)電流突變量Δi及電壓突變量Δu，經(jīng)計算得暫態(tài)功率ΔP，再對暫態(tài)功率ΔP進(jìn)行積分運算，便可得到暫態(tài)功率能量EP，即：

式中：N為采樣點個數(shù)。

當(dāng)發(fā)生反方向故障時，暫態(tài)功率能量EP>0；發(fā)生正方向故障時，暫態(tài)功率能量EP<0。

3.2 正向區(qū)內(nèi)外故障判據(jù)

對保護裝置采集的故障電流信號進(jìn)行6尺度小波分解，求取第1尺度高頻能量E1和第6尺度低頻能量E6。定義E1與E6的比值KP為：

式中：K為比例系數(shù)。

設(shè)置K目的是保證暫態(tài)能量能夠較好地被保護裝置測量到。綜合考慮采樣頻率及時間窗長度，取K=1 000。

利用此高低頻能量比值作為區(qū)內(nèi)外保護動作判據(jù)：當(dāng)KP滿足式（6）時，故障發(fā)生在保護區(qū)內(nèi)，保護動作；反之則為區(qū)外發(fā)生故障，保護不動作。

式中：Kset為保護動作閾值。

經(jīng)過仿真驗證并考慮保留一定裕度，取Kset為1。

保護流程如圖10所示。

圖10 保護方案流程圖Fig. 10 Flow chart of protection scheme

4 仿真研究

本文利用PSCAD/EMTDC仿真軟件搭建了貫通式同相牽引直接供電系統(tǒng)仿真模型。仿真模型參數(shù)：公用電網(wǎng)側(cè)電壓為 220 kV，牽引網(wǎng)側(cè)電壓為27.5 kV，系統(tǒng)頻率為50 Hz，電源容量為6 GVA；牽引變之間的線路長度為35 km；采樣時間選取5 ms。

4.1 不同故障初始角仿真實驗

選取不同故障初始角，分別在牽引網(wǎng)線路的不同位置設(shè)置相同過渡電阻（50 Ω）接地故障。

表2為保護方向判別及高低頻能量比值計算結(jié)果。由表2數(shù)據(jù)可見，保護對于不同故障初始角接地故障均能準(zhǔn)確動作。

表2 不同故障初始角保護判斷結(jié)果Tab. 2 Judgment results of different fault initial angle protection

4.2 不同故障距離仿真實驗

選定相同故障初始角（45°），在線路不同位置，以5 km為間隔距離，設(shè)置相同過渡電阻（0.1 Ω）接地故障進(jìn)行仿真。

仿真結(jié)果如表3所示。

表3 不同故障距離保護判斷結(jié)果Tab. 3 Judgement results of different fault distance protection

從表3可以看出，線路能夠得到可靠保護。

4.3 不同過渡電阻仿真實驗

仿真故障設(shè)置：相同故障初始角（60°），正方向不同位置，過渡電阻分別為 0.1 Ω、10 Ω、50 Ω、100 Ω、300 Ω；故障距離同樣以首端、中點、末端為例。

仿真結(jié)果表4所示。

表4 不同過渡電阻保護判斷結(jié)果Tab. 4 Judgment results of different transition resistance protection

由表4所示仿真結(jié)果可知，對于不同的故障位置、不同故障初始角、不同故障距離、不同過渡電阻，本文所提保護方案均能準(zhǔn)確可靠判別。

5 結(jié)論

本文利用牽引網(wǎng)邊界對故障電流信號高頻部分的衰減作用以及正反方向故障時暫態(tài)功率能量極性的差異，提出了一種基于小波暫態(tài)能量的牽引網(wǎng)單端量方向保護方案。

方案可快速準(zhǔn)確地區(qū)分牽引網(wǎng)線路區(qū)內(nèi)、外故障，在不同過渡電阻、不同故障初始角情況下均能可靠判別。

（1）牽引變電所供電段區(qū)外所并聯(lián)的電容與牽引網(wǎng)的一部分線路共同構(gòu)成的線路邊界是故障電流信號高頻部分產(chǎn)生較大幅值衰減的主要影響因素。

（2）線路反方向故障時，暫態(tài)功率能量的極性與正方向故障時的極性相反。

（3）線路正方向故障時，利用小波多尺度分解可提取故障電流信號的高低頻暫態(tài)能量，通過比較區(qū)內(nèi)外故障時的高頻暫態(tài)能量與低頻暫態(tài)能量比值可構(gòu)造正方向區(qū)內(nèi)外故障判據(jù)。