張志禹，朱宗耀，滿蔚仕

　　(西安理工大學 自動化與信息工程學院，陜西 西安 710048)

       摘要：S變換是一種由小波變換及短時傅里葉變換改進后的時頻分析工具，由于其良好的自適應分辨率近年來已被用于電力系統故障研究。但由于S變換運算過程存在窗函數之外的冗余運算，計算量較大，影響運算速度。為此提出了一種基于快速S變換的小電流接地方式下的故障選線新算法。利用快速S變換后各線路零序電流在特征頻帶內的模值進行分析選線。最后通過大量仿真算例證明了此方法能夠有效地進行故障選線，正確率高，并且節省了運算時間。

　　關鍵詞：小電流接地系統；時頻分析；快速S變換

　　中圖分類號：TM77文獻標識碼：ADOI： 10.19358/j.issn.1674-7720.2016.24.021

　　引用格式：張志禹，朱宗耀，滿蔚仕. 基于快速S變換的小電流單相接地故障選線［J］.微型機與應用，2016,35（24）：73-76.

0引言

　　我國配電網主要采用小電流接地方式，包括中性點不接地、中性點經消弧線圈接地和中性點經高阻抗接地方式。采用小電流接地方式的優勢是供電可靠，當發生單相接地故障后可以繼續運行1~2小時［1］。但若不能及時準確地檢測出故障線路進行檢修可能造成故障進一步擴大，對人身及設備安全造成危害，因此對小電流接地故障選線進行研究是非常必要的。

　　目前故障選線方法主要有基于穩態量的方法、基于暫態量的方法、穩態與暫態法相結合的選線方法和注入信號法等。故障發生后暫態量幅值大，易檢測，并且包含豐富的故障信息。由于近年來信號分析技術廣泛地應用于故障檢測與分析，因此基于暫態量的小電流接地系統故障選線技術得到了快速的發展。故障信號的能量譜在不同的故障發生時刻具有不同的特點，文獻［2］通過小波變換得到不同頻段的能量極大值，從而判定故障線路。文獻［3］對各線路的零序電流進行小波包分解，計算各線路故障后1/4周期的時頻特征量，通過相關系數來判定故障線路，仿真證明此法準確、可靠。文獻［4］利用S變換進行故障選線，通過對零序電流進行S變換提取各頻率點的模值和相角信息作為選線判據，但由于S變換過程中存在大量冗余計算，增大了運算量使運算時間變長，影響選線速度。

　　快速S變換在繼承了S變換良好的時頻分析能力的基礎上又解決了S變換運算效率低的問題，速度快，實時性好。對于饋線較多的系統采用此法能大幅度減少運算量。因此本文提出了一種利用快速S變換進行故障選線的新方法。利用MATLAB構建配電線路故障仿真模型，對此選線方法進行了大量仿真驗證，證明了此方法能準確、可靠地對小電流接地系統進行故障選線。

1S變換與快速S變換

　　1.1S變換

　　S變換是由小波變換與短時傅里葉變換（STFT）結合演變而來的，由STOCKWELL R G等人［5］首次提出。其數學表達式如下：

　　式（1）中：τ和v是S域的時間和頻率，ω(τ－t,ν)=是單位區域內所有頻率ν 的廣義窗函數。

　　由式（1）可看出S變換相對于STFT的優點在于高斯窗口的高度和寬度隨頻率而變化，這就克服了STFT窗口高度和寬度固定的缺陷［6］。從而對低頻信號具有較高的頻率分辨率，對高頻信號具有較高的時間分辨率。

　　1.2快速S變換

　　由于在S變換的過程中，采樣頻率是根據傅里葉變換來確定的，沒有考慮窗函數的影響，從而使變換處理過程中存在大量的冗余數據［7］。因此快速S變換將窗函數之外的冗余運算去除了，從而提高了運算速度。

　　S變換的時間復雜度為O(N3)，而快速S變換的時間復雜度為O(NlogN)，相比之下快速S變換的計算速度更快［6］，且采樣點數N越大，快速S變換的速度優勢體現越明顯。

　　1.3S變換與快速S變換效果對比

　　圖1為一個1 024點的測試信號，圖2為S變換（ST）及快速S變換（KST）的重構圖。可以看出S變換重構信號與快速S變換重構信號基本一致。

　　圖3、圖4分別為對測試信號進行S變換及快速S變換后的模系數時頻圖，橫軸對應采樣點，縱軸對應頻率。可以看出，相比較S變換，快速S變換由于去除了冗余信息使得時頻圖不再連續、平滑。但這些數據也完全足夠重構出原始信號。

2仿真模型

　　圖5所示為一個4饋線的10 kV工頻50 Hz的配電系統。線路L1、L2及L3的長度分別為40 km、30 km和50 km。消弧線圈按8%的過補償方式整定。采樣頻率取6 400 Hz。

　　中性點經消弧線圈接地系統后，由于消弧線圈電感電流的補償作用，當發生單相接地故障時，故障線路的暫態零序電流幅值最大這一結論不再成立。因此必須確定特征頻帶，特征頻帶的上限頻率即為所有線路暫態零序電流相位一致的首段頻帶［8］。可計算得本系統特征頻帶的上限頻率為1 200 Hz，即24倍工頻。而中性點經消弧線圈接地系統與中性點不接地系統的特征頻帶的下限頻率不一樣，中性點不接地時為0 Hz，而中性點經消弧線圈接地時為150 Hz，即3倍工頻。

　　因此取特征頻帶為3~24倍工頻。在此頻帶內選線算法不受中性點接地方式的影響。

3基于快速S變換的故障選線方法

　　3.1數據處理

　　中性點經消弧線圈接地系統，線路L1發生A相接地故障，故障點距離母線5 km，接地電阻20 Ω，故障合閘角為90°，對應第288個采樣點，檢測線路L1的零序電流如圖6所示。

　　當系統正常運行時，由于負荷電流包含零序分量、負荷存在三相不平衡以及受到電流互感器勵磁電流分量的影響，電流互感器輸出的零序電流并不為零（由于很小需放大才可以觀察得到），因此為消除干擾，需對暫態零序電流進行處理［9］：

　　i0i(t)=io0i［tf:(tf+T)］－io0i［(tf－T):tf］（2）

　　式（2）中tf為故障時刻，T為周期，i代表故障線路，io0i［tf:(tf+T)］表示故障后一周期的零序電流，io0i［(tf－T):tf］代表故障前一周期的零序電流。i0i(t)則為消除干擾后用于后續分析的零序電流。

　　采樣頻率設為6 400 Hz，因此一個周期T包含128個點，圖6中L1線路在288點處發生故障，因此對故障點后一個周波的數據進行快速S變換處理，并取模，得到如圖7所示的模時頻圖，橫軸為采樣點，縱軸為頻率，1~64行平均等分3 200 Hz。

　　頻率間隔△f=fs/N=6 400/128=50 Hz，由第2節的分析知特征頻帶取3~24行。由于暫態能量在故障后1/4的T內快速衰減，因此只取1~32列，又由于邊緣效應，因此去掉前6列，故最后只取7~32列。所以最終取3~24行，7~32列的模值矩陣A(22×26）。

3.2選線判據

　　定義mi（i=1，2，3代表線路）為對應線路的模值矩陣A(22×26）中所有元素的和。分別計算三條線路對應m1、m2和m3，比較大小，m值最大的即判定為故障線路。

4仿真驗證

　　在圖5的仿真模型中，分別對中性點不接地和經消弧線圈接地兩種方式進行了仿真。對所有線路距離出線端不同的位置設置A相接地故障，分別設置不同的故障合閘角（本文分別取0°、45°、90°），取不同的接地電阻（本文分別取0.01 Ω、20 Ω、200 Ω、2 000 Ω）進行仿真分析。經過大量的仿真驗證證明，在線路不同故障位置、不同故障初始角及不同接地電阻的情況下均能準確判斷出故障線路。以下列舉幾種典型情況進行驗證說明。

　　算例1：中性點不接地系統，故障線路為L1，故障位置距離出線端25 km，故障合閘角45°，接地電阻2 000 Ω。

　　分別計算出L1、L2和L3對應的m1=1579.6，m2=951.4681，m3=896.4882。因為m1最大，故判定故障線路為L1，選線正確。

　　算例2：中性點不接地系統，故障線路為L2，故障位置距離出線端10 km，故障合閘角0°，接地電阻20 Ω。

　　分別計算出L1、L2和L3對應的m1=2 109.4, m2=5 265.3，m3=3 286.2。m2最大，故判定故障線路為L2，選線正確。

　　算例3：中性點不接地系統，故障線路為L3，故障位置距離出線端40 km，故障合閘角90°，接地電阻0.01 Ω。

　　分別計算出L1、L2和L3對應的m1=12 943，m2=12 582，m3=23 325。m3最大，故判定故障線路為L3，選線正確。

　　算例4：中性點經消弧線圈接地系統，故障線路L1，故障位置距離出線端5 km，故障合閘角90°，接地電阻20 Ω。

　　分別計算出L1、L2和L3對應的m1=37 466，m2=17 144，m3=25 776。m1最大，故判定故障線路為L1，選線正確。

　　算例5：中性點經消弧線圈接地系統，故障線路為L2，故障位置距離出線端20 km，故障合閘角45°，接地電阻200 Ω。

　　分別計算出L1、L2和L3對應的m1=4 487.5，m2=11 645，m3=6 205.7。m2最大，故判定故障線路為L2，選線正確。

　　算例6：中性點經消弧線圈接地系統，故障線路為L3，故障位置距離出線端30 km，故障合閘角0°，接地電阻200 Ω。

　　分別計算出L1、L2和L3對應的m1=1 960.2，m2=1 529.7，m3=2 555.7。m3最大，故判定故障線路為L3，選線正確。

5結論

　　本文將快速S變換（FST）應用于小電流接地故障選線，提出了一種新的故障選線方法，克服了傳統依靠S變換進行故障選線時運算量大、速度慢的問題。尤其當饋線較多，處理的數據量大時，本方法的優勢體現更明顯。經過大量仿真，證明本方法可以準確、可靠地進行故障選線。

　　參考文獻

　　［1］ 于永源,楊綺雯. 電力系統分析［M］. 北京:中國電力出版社,2007.

　　［2］ 吳樂鵬,黃純,林達斌,等. 基于暫態小波能量的小電流接地故障選線新方法［J］.電力自動化設備,2013,33（5）:7075.

　　［3］ 束洪春,朱夢夢,黃文珍,等.基于暫態零序電流時頻特征量的配網故障選線方法［J］.電力自動化設備,2013,33（9）:16.

　　［4］ 張鈞,何正友,賈勇.基于S變換的故障選線新方法［J］.中國電機工程學報,2011,31（10）:109115.

　　［5］ STOCKWELL R G，MANSINHA L，LOWE R P． Localization of the complex spectrum：the Stransform［J］．IEEE Transactions on Signal Processing，1996，44(4)：9981001．

　　［6］ 張志禹,滿蔚仕,郗壘,等. 快速S變換在電能質量分析中的應用［J］. 電網技術，2013,37（5）：12851290.

　　［7］ 吳彥華. 改進的離散S變換快速算法與連續小波變換算法性能分析［J］. 信號處理，2012，28（7）：973979.

　　［8］ 張保會,趙慧梅,張文豪,徐靖東. 基于特征頻帶內暫態零序電流特點的配電網單相接地故障選線方法［J］. 電力系統保護與控制，2008，36（13）：510.

　　［9］ 蘇杰鋒. 基于S變換的小電流接地故障選線方法的研究［D］.廣州：華南理工大學,2012.