齊 鄭， 趙傳宗， 紀(jì) 鵬， 張 海， 王 瑩

(1. 華北電力大學(xué) 電氣與電子工程學(xué)院， 北京 102206； 2. 國(guó)網(wǎng)遼寧省電力有限公司 a. 撫順供電公司， b. 人董部， 遼寧 撫順 113008； 3. 國(guó)家電投集團(tuán)東北電力有限公司 撫順撫電能源分公司， 遼寧 撫順 113008)

小電流接地系統(tǒng)發(fā)生單相接地故障時(shí)，其暫態(tài)電壓電流信號(hào)含有豐富的時(shí)域、頻域特征.暫態(tài)接地電流包括了容性分量和感性分量，由等效電路圖易得暫態(tài)容性電流流經(jīng)故障路徑和非故障路徑，感性電流只分布在故障路徑上.在配電網(wǎng)各個(gè)節(jié)點(diǎn)可以快速準(zhǔn)確地從暫態(tài)接地電流中提取感性電流分量，結(jié)合拓?fù)浣Y(jié)構(gòu)能夠?qū)崿F(xiàn)故障點(diǎn)的在線區(qū)段定位.感性電流的主要表現(xiàn)形式是衰減的直流分量，信號(hào)特征明顯不易與豐富的諧波相互混淆，因此在故障定位上具有一定的優(yōu)勢(shì).對(duì)于小電流接地系統(tǒng)，中性點(diǎn)采用了非有效接地方式導(dǎo)致接地電流較小，從而對(duì)實(shí)現(xiàn)暫態(tài)信號(hào)的有效信息提取提出了較高的要求.

從復(fù)雜暫態(tài)信號(hào)中提取目標(biāo)信號(hào)的常用方法有傅里葉變換[1-2]、小波變換[3]和普羅尼算法等.這些經(jīng)典的非自適應(yīng)算法往往采用有限正交基函數(shù)展開形式來表現(xiàn)實(shí)際的復(fù)雜信號(hào)，因此非自適應(yīng)算法具有局限性.對(duì)于接地暫態(tài)信號(hào)這種具有大量豐富時(shí)域、頻域特征的信號(hào)，自適應(yīng)原子稀疏分解算法可以實(shí)現(xiàn)精確的稀疏表達(dá)[4-5]，其運(yùn)用過完備原子庫(kù)代替?zhèn)鹘y(tǒng)算法的正交基函數(shù)，避免了表現(xiàn)形式的局限.

基于原子稀疏分解理論，本文采用匹配追蹤算法分解接地暫態(tài)電流，自適應(yīng)地提取衰減分量成分，并結(jié)合拓?fù)浣Y(jié)構(gòu)定位故障點(diǎn).

1 衰減直流分量法定位原理

諧振接地系統(tǒng)單相接地故障零序等值回路如圖1所示.圖1中U0為零序電壓，L0為零序等值電感；R0為接地電阻零序等值電阻；C0為三相對(duì)地電容；RL和L為消弧線圈的有功損耗電阻和電感.

圖1 單相接地故障零序等值回路Fig.1 Zero-sequence equivalent circuit ofsingle-phase grounding fault

根據(jù)接地零序等值電路圖，可以得到單相接地故障的接地電流及衰減直流分量表達(dá)式為

(1)

(2)

式中：iC為容性電流分量；iL為感性電流分量；A、Af、ADC分別為各項(xiàng)幅值；ωf、ω分別為各項(xiàng)頻率；τL、τC為電感和電容回路的時(shí)間常數(shù)；φ為故障發(fā)生時(shí)刻故障相電壓相角；Uphm為非故障條件下電源相電壓幅值.

由零序等值網(wǎng)絡(luò)可知，以衰減直流分量為代表的感性分量存在于故障路徑上，區(qū)別于存在故障路徑和非故障路徑的容性分量，因此可以判定故障路徑并進(jìn)一步定位故障點(diǎn).具體策略為：在配電網(wǎng)母線和各個(gè)分支節(jié)點(diǎn)處測(cè)定衰減直流分量的有無，含有衰減直流分量的節(jié)點(diǎn)位于故障路徑上，故障路徑終點(diǎn)節(jié)點(diǎn)的負(fù)荷側(cè)區(qū)段為故障區(qū)段.以圖2為例，假設(shè)D點(diǎn)發(fā)生了單相接地故障，分別在A～H節(jié)點(diǎn)測(cè)定衰減直流分量的有無，若提取策略有效可以測(cè)算出OABCD節(jié)點(diǎn)存在衰減直流分量，E、F、G、H節(jié)點(diǎn)不存在衰減直流分量，則可確定故障路徑為OABCD，進(jìn)一步可以確定故障路徑終點(diǎn)區(qū)段DE為接地故障所在區(qū)段.

圖2 故障定位原理圖Fig.2 Principle diagram of fault location

2 接地暫態(tài)信號(hào)的稀疏分解

當(dāng)系統(tǒng)發(fā)生單相接地故障時(shí)，暫態(tài)信號(hào)通常由工頻分量、諧波分量、雜音和直流分量組成，其信息量巨大，且含有豐富的頻域、時(shí)域特征且隨時(shí)間變化易出現(xiàn)非正弦畸變[6-7]，這使得實(shí)際的暫態(tài)信號(hào)分解難度較大.傳統(tǒng)的非自適應(yīng)算法通常采用單一正交的信號(hào)投影方法，即

y=ψx

(3)

式中：x為原始信號(hào)；ψ為固定正交映射；y為信號(hào)的投影.單獨(dú)運(yùn)用某一個(gè)固定的正交基ψ表達(dá)任意的復(fù)雜信號(hào)難以實(shí)現(xiàn)精確的稀疏表達(dá)[8]，導(dǎo)致故障路徑上感性和容性分量成分難以區(qū)分.在非自適應(yīng)算法中，具有代表性的有傅里葉變換、小波分析和PRONY算法，其中，傅里葉實(shí)現(xiàn)了信號(hào)從時(shí)域向頻域的變換，可以在一定精度上提取工頻分量和諧波，但是難以辨識(shí)衰減直流分量；小波變換適用于非平穩(wěn)信號(hào)的分解，具備對(duì)暫態(tài)電流的濾波去噪的能力，但是對(duì)于給定的信號(hào)特征進(jìn)行提取存在一定難度[9]；PRONY算法適用于特征信號(hào)的提取，但是其龐大的運(yùn)算量和維度限制了其應(yīng)用范圍.

2.1 信號(hào)稀疏分解理論

(4)

常見的構(gòu)建方法是用多種核函數(shù)來表示原子，本文采用了GABOR原子庫(kù)，GABOR原子由調(diào)制高斯窗函數(shù)組成[11]，即

(5)

式中：G(t)為高斯窗函數(shù)；s為伸縮因子；u為平移因子；v為頻率；w為相位.對(duì)上述分組進(jìn)行進(jìn)一步調(diào)制，形成對(duì)應(yīng)原子庫(kù).

稀疏分解的自適應(yīng)性對(duì)解決實(shí)際暫態(tài)信號(hào)成分復(fù)雜，衰減直流分量難以提取的問題有促進(jìn)作用，可以更好地刻畫暫態(tài)信號(hào)特征.算法的自適應(yīng)性可以提高暫態(tài)信號(hào)分解的精確性和分解過程的靈活性，更加準(zhǔn)確地提取衰減直流分量以實(shí)現(xiàn)故障定位.

2.2 運(yùn)用匹配追蹤算法實(shí)現(xiàn)暫態(tài)信號(hào)提取

多重子信號(hào)逼近是復(fù)雜暫態(tài)信號(hào)成分分析和處理的核心問題之一，n原子逼近問題可以描述為：在一個(gè)龐大的原子庫(kù)中，按照既定策略選擇n個(gè)子信號(hào)(原子)，實(shí)現(xiàn)復(fù)雜暫態(tài)信號(hào)在原子子空間中投影能量最大.但在實(shí)際應(yīng)用時(shí)，對(duì)n的數(shù)值有一定限制，實(shí)際暫態(tài)接地電流信號(hào)由工頻電流、3、5、7、9次諧波和直流成分組成，使得n的數(shù)值有限避免了計(jì)算復(fù)雜度組合爆炸.匹配追蹤基于“貪婪”算法的核心思想，每次搜索和剩余信號(hào)最接近的原子，通過局部最優(yōu)解逐步逼近原始信號(hào).

對(duì)于原子庫(kù)D={gk，k=0，1，…，K}和實(shí)際單相接地暫態(tài)電流信號(hào)F，將暫態(tài)電流反映到原子庫(kù)中的原子g0上，即

s.t.F=〈F，g0〉g0+R

|〈F，g0〉|≥asup|〈F，g0〉|

a∈(0，1]

(6)

通過“貪婪策略”選擇原子g0，使得暫態(tài)剩余信號(hào)R最小，|〈F，g0〉|最大，并將R視為新的F帶入下一次迭代，繼續(xù)選擇新的原子g1直到

(7)

3 配電網(wǎng)單相接地故障區(qū)段定位

對(duì)采集的不同區(qū)段接地暫態(tài)電流信號(hào)進(jìn)行稀疏分解，可提取最優(yōu)原子中符合衰減直流分量特性的原子，記為衰減直流.計(jì)算衰減直流幅值，同時(shí)選取故障發(fā)生在相位角φ時(shí)的接地暫態(tài)電流的有效值為閾值.由于衰減直流只存在于故障路徑上，因此可以區(qū)分故障路徑和非故障路徑，故障路徑的終端為接地故障點(diǎn).故障定位的具體步驟為：

1) 基于GABOR方法構(gòu)建原子庫(kù)；

2) 配電網(wǎng)某處發(fā)生了單相接地故障，現(xiàn)場(chǎng)PMU記錄信號(hào)上傳主站；

3) 對(duì)記錄的暫態(tài)接地電流信號(hào)設(shè)置分解參數(shù)；

4) 尋找最優(yōu)原子，計(jì)算新的信號(hào)殘差；

5) 迭代計(jì)算，實(shí)現(xiàn)對(duì)暫態(tài)信號(hào)的稀疏表達(dá)；

6) 對(duì)于若干個(gè)原子，選取符合衰減直流特征的對(duì)應(yīng)原子，計(jì)算幅值大小；

7) 根據(jù)衰減直流分量的有無確定故障路徑，路徑終端為故障點(diǎn).

4 仿真波形分析

某實(shí)際的35 kV中性點(diǎn)不接地配電網(wǎng)結(jié)構(gòu)如圖3所示，D1～D10為配電網(wǎng)不同區(qū)段，每個(gè)區(qū)段進(jìn)線端裝設(shè)有FTU裝置，A、B、…、I為饋線分段開關(guān).該配電網(wǎng)在D4區(qū)域發(fā)生單相接地故障，當(dāng)故障發(fā)生時(shí)，系統(tǒng)帶故障運(yùn)行，各個(gè)區(qū)段內(nèi)FTU裝置對(duì)暫態(tài)信號(hào)進(jìn)行采樣.

圖3 某35 kV配網(wǎng)結(jié)構(gòu)示意圖Fig.3 Schematic structure of 35 kV distribution network

以D4區(qū)段為例，F(xiàn)TU裝置記錄的接地零序電流暫態(tài)波形如圖4所示，采用匹配追蹤算法分解現(xiàn)場(chǎng)電流得到的結(jié)果如圖5所示.

圖4 D4區(qū)段零序電流Fig.4 Zero-sequence current of D4 section

對(duì)匹配的最佳原子分析可得，分解的工頻子信號(hào)、3、5、7、9次諧波子信號(hào)和衰減直流子信號(hào)符合實(shí)際.其中，諧波分量大約經(jīng)過15 ms左右衰減為零；剩余的未完全匹配的信號(hào)殘差視為噪音.通過計(jì)算噪音有效值，其低于所有分量的有效值，提取的衰減直流分量為

(8)

其衰減時(shí)間常數(shù)為1.6 ms，最大幅值為26.4 A，幅值特征明顯.由于衰減時(shí)間常數(shù)較短，可以被電流互感器準(zhǔn)確采集，同時(shí)，中性點(diǎn)非有效接地配電網(wǎng)，其暫態(tài)接地電流強(qiáng)度有限，不會(huì)出現(xiàn)因衰減直流產(chǎn)生的磁鏈積累導(dǎo)致互感器二次側(cè)出現(xiàn)局部暫態(tài)飽和的情況[12].

同理，對(duì)其余區(qū)段的暫態(tài)信號(hào)作稀疏分解，提取對(duì)應(yīng)的直流分量原子并進(jìn)行分析對(duì)比，對(duì)比結(jié)果如圖6所示.依據(jù)直流分量的幅值特性劃分故障路徑，可以得出D4區(qū)段位于故障路徑，其相鄰區(qū)段D5不屬于故障路徑，進(jìn)而定位接地點(diǎn)位于D4區(qū)段內(nèi).

圖5 現(xiàn)場(chǎng)數(shù)據(jù)分解結(jié)果Fig.5 Decomposition results of field data

圖6 衰減直流分量對(duì)比圖Fig.6 Comparison diagram of attenuatedDC component

5 結(jié) 論

針對(duì)小電流接地系統(tǒng)單相接地故障定位問題，本文提出了一種在線區(qū)段定位方法.通過原子稀疏分解理論對(duì)暫態(tài)接地電流進(jìn)行稀疏表達(dá)，提取衰減直流分量劃分故障路徑，并進(jìn)一步定位故障點(diǎn).自適應(yīng)原子稀疏分解理論可以避免傳統(tǒng)非自適應(yīng)性算法的局限，擴(kuò)大了信號(hào)描述范圍，更加適用于復(fù)雜的接地暫態(tài)信號(hào)分解.通過匹配追蹤算法，可快速準(zhǔn)確提取衰減直流成分，結(jié)合配電網(wǎng)拓?fù)浣Y(jié)構(gòu)實(shí)現(xiàn)在線區(qū)段定位.當(dāng)前，根據(jù)直流分量進(jìn)行故障定位的方法尚不適用于故障合閘角為90°的狀態(tài)，如何解決該問題需要進(jìn)一步的研究.