基于VMD 和改進(jìn)聚類算法的配電網(wǎng)故障選線方法

王遠(yuǎn)川，李澤文，夏翊翔，毛紫玲，郭欣玉

（長沙理工大學(xué)電氣與信息工程學(xué)院，長沙 410114）

目前我國中低壓配電網(wǎng)普遍采用小電流接地方式，大數(shù)據(jù)統(tǒng)計(jì)顯示，一年中單相接地故障發(fā)生的次數(shù)達(dá)到配電網(wǎng)所有發(fā)生故障的80%左右[1]，故障可能會引起系統(tǒng)局部電壓過高，進(jìn)而對電氣設(shè)備造成破壞[2]，嚴(yán)重時(shí)甚至?xí)＜叭藛T的生命安全，因此需及時(shí)找出故障線路并予以切除。

基于主動注入法[3-5]和被動檢測法的故障選線方法是目前常用的兩種方法。通常，主動注入法易受過渡電阻的影響，且成本相對較高，適用性不強(qiáng)；被動檢測法可分為基于穩(wěn)態(tài)信息[6-7]和基于暫態(tài)信息的故障選線法[8-10]。實(shí)際運(yùn)行中，諸多因素都會對傳統(tǒng)的基于穩(wěn)態(tài)信息的故障分析法造成干擾[11]，使得該方法很不可靠[12]。相比之下，基于暫態(tài)信息的故障分析法，憑借自身諸多優(yōu)點(diǎn)[13-14]，被廣泛應(yīng)用于電力系統(tǒng)故障選線、行波定位等電網(wǎng)保護(hù)領(lǐng)域。許多故障特性，例如信號的相位、幅值、極性變化等，都存在于高頻暫態(tài)信號中[15-17]，如果能被有效地分解、提取出來，就能夠?yàn)檫x線提供思路。例如，文獻(xiàn)[18]采用小波分解將故障零序電流分解后，通過比較各線路在特征頻帶中的幅值和極性來鑒別故障線路。但在特定分解尺度下，小波變換不能同時(shí)在時(shí)間和頻率上都具有很高的精度，以及小波變換是非自適應(yīng)性的，基函數(shù)選定了之后就不能在分解過程中改變，其選擇可能會對分解結(jié)果產(chǎn)生影響。針對這些問題，經(jīng)驗(yàn)?zāi)B(tài)分解EMD（empirical mode decomposition）算法被運(yùn)用在故障選線領(lǐng)域，其能依據(jù)數(shù)據(jù)自身的時(shí)間尺度特征來將信號分解為若干模態(tài)分量，如文獻(xiàn)[19]使用EMD算法對各線路導(dǎo)納進(jìn)行分解，提取出特征分量，通過比較故障線路和非故障線路的特征分量符號進(jìn)行選線。但是EMD 分解存在端點(diǎn)效應(yīng)和頻率混疊等問題。此后，又相繼誕生了集合經(jīng)驗(yàn)?zāi)B(tài)分解EEMD（ensemble empirical mode decomposition）和互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解CEEMD（complementary ensemble empirical mode decomposition）等改進(jìn)型算法，但始終存在分解效果不理想的問題。

有一類利用形態(tài)學(xué)選線的方法，其對分解后的零序電流采用相關(guān)性聚類算法來鑒別故障線路，如文獻(xiàn)[20]使用凝聚層次聚類HAC（hierarchical agglomerative clustering）算法來計(jì)算不同線路零序電流相似性測度并聚類，但HAC 算法時(shí)間復(fù)雜度高，且具有貪心算法的缺點(diǎn)。又如，文獻(xiàn)[21]使用的基于劃分算法的K-means 聚類算法具有步驟簡單、時(shí)間復(fù)雜度、空間復(fù)雜度低的優(yōu)點(diǎn)，但該算法易受初始聚類中心選取的影響[22]，且運(yùn)算效率不及在其基礎(chǔ)之上改進(jìn)的一系列算法。

針對上述情況，為了高效、準(zhǔn)確地選出故障線路，本文提出了一種基于變分模態(tài)分解VMD 和Kmeans++聚類算法的故障選線方法。首先對提取到的零序電流信號進(jìn)行VMD 分解，得到有限個(gè)自適應(yīng)頻帶特征的本征模態(tài)函數(shù)IMF（intrinsic mode function）；而后通過對其中的低頻分量作相關(guān)性計(jì)算，對高頻分量測算其初始極性，從而確定代表各出線零序電流信號的散點(diǎn)在二維平面內(nèi)的橫、縱坐標(biāo)；最后通過K-means++聚類算法檢測點(diǎn)集的分布情況，將故障線路篩選出來。

1 故障零序電流的分析計(jì)算

1.1 小電流接地系統(tǒng)解析

配電網(wǎng)發(fā)生單相接地故障示意如圖1 所示。圖1 中，rp、Lp分別為消弧線圈的等效電阻和等效電感值；Rf為接地過渡電阻；C0j為第j條線路Lj的對地電容大??；i0j為線路Lj上流過的零序電流；iL為消弧線圈上流過的感性電流；U0f為零序電壓源。

圖1 零序電流分布Fig.1 Zero-sequence current distribution

由圖1 可知，發(fā)生單相接地故障后，系統(tǒng)中各條健全線路首端流過其本身產(chǎn)生的容性電流，這些電流都將經(jīng)過接地點(diǎn)流入故障線路中。因此，零序電流在故障和健全線路中的流向剛好相反，在非故障線路中為出線端向線路末端，而在故障線路中則為線路末端向出線端，這也導(dǎo)致故障線路中的容性電流與健全線路中的容性電流初始極性理論上也是相反的。值得注意的是，消弧線圈中此時(shí)將產(chǎn)生一個(gè)與接地點(diǎn)處容性電流相位相反的感性電流，該電流經(jīng)電感線圈注入大地再經(jīng)接地故障點(diǎn)流入故障線路中。

1.2 故障零序電流的分析

建立單相接地故障時(shí)的暫態(tài)分析等值電路模型，該電路如圖2所示。

圖2 系統(tǒng)暫態(tài)分析等值電路Fig.2 Equivalent circuit of transient analysis of system

圖2 中，R0、L0分別為配電網(wǎng)的等值電阻、電感；C∑為系統(tǒng)零序電容之和；if為總故障電流；ic為暫態(tài)電容電流。

根據(jù)基爾霍夫電壓定理列寫圖2 中等效電路回路的微分方程為

式中：φ為故障時(shí)刻電壓初始相角；Um為系統(tǒng)相對地的電壓幅值；ω為工頻角頻率；t為時(shí)間。

對上述微分方程求解，可得故障點(diǎn)處流過的容性電流和感性電流的計(jì)算公式為

式中：、分別為容性電流的高頻振蕩與工頻交流分量；為感性電流的工頻交流分量；iLdc為感性電流的衰減直流分量；Icm、ILm分別為容性電流、感性電流的幅值大??；ωf為暫態(tài)自由振蕩分量的角頻率；β為自由振蕩分量的衰減系數(shù)；τL為回路的時(shí)間常數(shù)。

容性電流和感性電流相加可得故障點(diǎn)處流過的總電流if，即

由于衰減直流分量只存在于感性電流中，因此其只會隨著感性電流從接地故障點(diǎn)流入故障線路。根據(jù)接地故障點(diǎn)電流表達(dá)式可以看出，故障零序電流的特征不僅與配電線路本身一系列固有參數(shù)相關(guān)，其各組成部分的幅值也會隨著故障初相角取值大小的改變而變化。根據(jù)容性電流高頻振蕩分量的表達(dá)式可知，其波形振幅很大程度上取決于sinφ的大??；而對于感性電流中的衰減直流分量，其波形振幅大小取決于cosφ的大小。因此當(dāng)故障電流直流分量數(shù)值較大時(shí)，高頻振蕩分量一般較小；反之，當(dāng)高頻振蕩分量數(shù)值較大時(shí)，衰減直流分量一般較小。

例如，當(dāng)φ=0°時(shí)，故障電流表達(dá)式可化簡為

從式（6）可以看出，接地故障點(diǎn)處零序電流含有全部的3 個(gè)分量，且衰減直流分量幅值較大，而高頻分量幅值較小。

當(dāng)φ=90°時(shí)，故障電流表達(dá)式可化簡為

從式（7）可以看出，此時(shí)故障零序電流主要由高頻分量與工頻分量所組成，且高頻分量有著較大的幅值，而幾乎不含有直流分量。

2 變分模態(tài)分解基本原理

VMD算法具有自適應(yīng)和完全非遞歸性[23]，其將信號的分解問題轉(zhuǎn)化為構(gòu)建變分問題并求解，尋找變分問題的最優(yōu)解，保證分解后得到的序列為具有中心頻率和有限帶寬的模態(tài)分量，實(shí)現(xiàn)各頻率下本征模態(tài)函數(shù)IMF（intrinsic mode function）的有效分離，得到預(yù)設(shè)分級層數(shù)下的K個(gè)模態(tài)分量。其核心步驟由變分約束問題的構(gòu)造和求解構(gòu)成。

首先，構(gòu)造變分約束問題，其數(shù)學(xué)表達(dá)式為

式中：f為待分解的原始信號；?t為梯度運(yùn)算；δ(t)為單位脈沖；j為虛數(shù)單位；{uk} ={u1,u2,…,uK} 為分解后得到的K個(gè)本征模態(tài)函數(shù)；{ωk} ={ω1,ω2,…,ωK} 為各模態(tài)分量的中心頻率。通過引入懲罰因子α和拉格朗日算子λ，構(gòu)造增廣拉格朗日表達(dá)式，即

通過使用交替乘法算子迭代、計(jì)算式（9），得到其鞍點(diǎn)即為變分約束問題的最優(yōu)解。各個(gè)模態(tài)分量和其中心頻率的迭代公式為

如果滿足，則將目前得到的各IMF作為輸出量輸出；若不滿足，則回到式（10）和（11）繼續(xù)進(jìn)行迭代更新。

3 選線相關(guān)算法原理及判據(jù)構(gòu)造

結(jié)合第1.2節(jié)對發(fā)生單相接地故障時(shí)故障零序電流組成及其相關(guān)特征的解析可知，含有衰減直流分量的感性電流只流入故障線路中，這將使故障發(fā)生初期故障線路零序電流低頻分量的幅值遠(yuǎn)大于健全線路，從而產(chǎn)生波形上的差異；又由于零序電流在故障和健全線路上的流向相反，故障線路與健全線路零序電流信號高頻分量的初始極性也應(yīng)當(dāng)是相反的，因此無論是低頻分量還是高頻暫態(tài)分量，故障線路相比于非故障線路均存在著明顯的差異，如果能將這種差異投射到二維平面或者三維空間中，將有利于篩選出故障線路。

第1.2 節(jié)提到，故障電流衰減直流分量數(shù)值較大時(shí)，高頻振蕩分量往往較??；反之當(dāng)高頻振蕩分量較大時(shí)，衰減直流分量往往較小。換言之，一方差異較為明顯時(shí)，另一方差異有可能難以辨識。因此，本文選擇綜合考慮二者的差異，將其投射到二維平面中，通過對電流信號低頻分量相關(guān)系數(shù)與高頻分量初始極性差異構(gòu)造的二維平面點(diǎn)集進(jìn)行聚類，從而準(zhǔn)確選出故障線路。

3.1 VMD 分解相關(guān)參數(shù)取值的確定

對采樣信號做VMD 分解之前，需要提前設(shè)定好分解層數(shù)K、懲罰因子α以及預(yù)設(shè)精度ε[24-25]。

故障零序電流信號經(jīng)過VMD 分解后將得到K個(gè)子模態(tài)分量，每個(gè)子模態(tài)都有各自的中心頻率。VMD算法與EMD等其他分解方法的一個(gè)很大的不同點(diǎn)在于，其需要人為預(yù)先設(shè)定K的取值。

零序電流信號中的高頻振蕩分量的振蕩頻率浮動區(qū)間大致為1 000～3 000 Hz，與工頻分量和直流分量相差巨大；而工頻分量與直流分量之間頻率相差不大?？紤]到實(shí)際采集的信號往往含有噪聲，因此本文所有的仿真信號都將注入25 dB的噪聲作為干擾。

本文最終采取的分解方式為：先使用FIR 濾波器濾去頻率過高的噪聲信號（在仿真實(shí)驗(yàn)時(shí)取5 500 Hz 作為濾波界線），而后截取故障發(fā)生時(shí)刻后2 個(gè)信號周期（40 ms）的零序電流信號先做2 層分解，得到的模態(tài)分量分別記為IMF1和IMF2。高頻振蕩分量與過濾后剩下的噪聲信號將全部分解到IMF2中，而衰減直流、工頻交流分量和低次諧波則全部存在于IMF1中。將得到的IMF1分量再次做2層分解，得到的模態(tài)分量分別記為IMF1-1與IMF1-2，這樣衰減直流分量和工頻分量被全部分解到IMF1-1中，低次諧波分量則被分解到了IMF1-2中。由此各模態(tài)分量根據(jù)頻率由低到高（IMF1-1、IMF1-2、IMF2）的主要成分如表1所示。

表1 各模態(tài)分量的主要成分Tab.1 Main components of each modal component

懲罰因子α的取值會對模態(tài)分量的帶寬產(chǎn)生影響，研究表明α值的大小與模態(tài)分量的帶寬成反比例關(guān)系。α的值越小，信號帶寬越大，各模態(tài)分量的函數(shù)會有很大一部分頻率重疊部分，即模態(tài)混疊；α的值越大，信號帶寬越窄，模態(tài)分量向中心頻率收斂速度越快，故障波形中的重要信息存在丟失的可能性。根據(jù)以往的VMD分解經(jīng)驗(yàn)，α的取值范圍通常在2 000～4 000。本文經(jīng)過多次實(shí)驗(yàn)，結(jié)合電網(wǎng)實(shí)際運(yùn)行中的工況，選擇α的取值為4 000，預(yù)設(shè)精度ε的值取為1×10-6。圖3 是某條件下的故障線路信號波形，其分解結(jié)果如圖4所示。

圖3 某條件下的故障線路信號波形Fig.3 Fault line signal waveform under certain conditions

圖4 信號分解結(jié)果Fig.4 Signal decomposition results

3.2 二維平面散點(diǎn)圖的構(gòu)造

選擇采用數(shù)學(xué)形態(tài)學(xué)中皮爾遜相關(guān)系數(shù)來反映不同出線上零序電流低頻分量IMF1-1的相關(guān)性，其公式為

式中：ai、bi分別為a、b兩個(gè)信號的第i個(gè)采樣點(diǎn)的值；n為采樣點(diǎn)個(gè)數(shù)；ρ表示兩者的相關(guān)性，-1≤ρ≤1。ρ越趨近于-1，兩波形之間相關(guān)性越低；ρ越趨近1，兩波形之間的相關(guān)性越高。假設(shè)配電網(wǎng)共有m回出線，對各出線首個(gè)信號周期（20 ms）內(nèi)的IMF1-1分量進(jìn)行相關(guān)性計(jì)算，得到它們之間的相關(guān)系數(shù){ρij} ={ρ11,ρ12,…,ρmm} ，ρij表示出線Li與Lj低頻信號分量IMF1-1之間的相關(guān)系數(shù)。將每一個(gè)出線看作二維平面的一個(gè)散點(diǎn)，各散點(diǎn)橫坐標(biāo)為該點(diǎn)所代表的出線IMF1-1分量與其他出線（包括其自身）的IMF1-1分量相關(guān)系數(shù)之和。例如，代表線路Li的點(diǎn)的橫坐標(biāo)xi表達(dá)式為

各散點(diǎn)縱坐標(biāo)構(gòu)造方法如下。首先對各線路IMF2分量做一階差分計(jì)算，即

式中：diff為一階差分計(jì)算符號；IMFi-2表示出線Li的高頻分量；n為采樣點(diǎn)個(gè)數(shù)，當(dāng)n=1時(shí)，式（15）為信號的一階差分初值。

然后使用符號函數(shù)計(jì)算其一階差分初值的極性，確定各線路高頻分量的極性，以表示其縱坐標(biāo)。例如，代表線路Li的點(diǎn)的縱坐標(biāo)yi表達(dá)式為

式中，sgn為符號函數(shù)，當(dāng)信號初始極性為正時(shí)，sgn的值為1；當(dāng)信號極性為負(fù)時(shí)，其值為-1。

故障特征較為明顯時(shí)，故障線路零序電流低頻IMF1-1分量波形與健全線路存在明顯差異，因此它們之間的相關(guān)系數(shù)也將較小，而健全線路之間的波形差異則相對較小，相關(guān)系數(shù)較大并趨近于1。根據(jù)式（14），做相加運(yùn)算后，故障線路的散點(diǎn)橫坐標(biāo)值會較小，而健全線路的散點(diǎn)橫坐標(biāo)值則較大，且數(shù)值比較接近；又由于高頻IMF2分量初始極性存在差異，故障線路與健全線路散點(diǎn)的縱坐標(biāo)值將為相反數(shù)。即使某一頻段下故障特征不明顯導(dǎo)致單一坐標(biāo)無法區(qū)分故障線路，在二維平面中代表故障線路的散點(diǎn)依然是一個(gè)離群點(diǎn)。

3.3 K-means++聚類算法

基于劃分聚類的K-means++算法的原理簡單來說就是對一堆散點(diǎn)進(jìn)行聚類，首先確定將這些散點(diǎn)聚成幾類，然后挑選幾個(gè)點(diǎn)作為初始中心點(diǎn)，并根據(jù)一定規(guī)則迭代重置聚類中心點(diǎn)，直到達(dá)到類中的點(diǎn)都足夠近、類間的點(diǎn)都足夠遠(yuǎn)的目標(biāo)效果。傳統(tǒng)的K-means 聚類算法采取隨機(jī)選取的方式確定初始聚類中心，但最終的聚類結(jié)果可能會受到不同初始聚類中心選取的影響，從而對選線結(jié)果產(chǎn)生影響。K-means++的改進(jìn)之處在于，其在選取初始聚類中心時(shí)可以確保初始聚類中心之間的距離盡可能地遠(yuǎn)，從而更高效、準(zhǔn)確地進(jìn)行聚類分析。

在本文中故障線路和非故障線路各自的散點(diǎn)分別形成一個(gè)簇（記為故障線路簇和健康線路簇），因此需要將所有散點(diǎn)劃分為2類。由于故障線路所屬的散點(diǎn)理論上屬于離群點(diǎn)，K-means++聚類算法在首輪分類時(shí)即可將代表故障線路的散點(diǎn)單獨(dú)劃分出來，從而大大提高了運(yùn)算速度和選線準(zhǔn)確率。

該算法通過對各出線零序電流信號IMF1-1分量相關(guān)系數(shù)和IMF2分量初始極性差異構(gòu)造的點(diǎn)集P={P1,P2,…,Pm} 進(jìn)行聚類時(shí)的具體步驟為如下。

步驟1從輸入的數(shù)據(jù)點(diǎn)集P中隨機(jī)選一個(gè)點(diǎn)作為故障線路簇Q1和健康線路簇Q2的首個(gè)初始聚類中心c1。

步驟2分別計(jì)算點(diǎn)集中其他的點(diǎn)到該聚類中心的歐氏距離，根據(jù)距離盡可能遠(yuǎn)地原則選擇一個(gè)點(diǎn)作為另一個(gè)初始聚類中心c2。用xci、yci表示聚類中心ci的橫、縱坐標(biāo)，第i個(gè)點(diǎn)到其的歐氏距離dis表示為

步驟3根據(jù)點(diǎn)集中所有點(diǎn)與這兩個(gè)聚類中心的距離最短原則將它們劃分到相應(yīng)的線路簇。

步驟4重新計(jì)算兩個(gè)簇各自的聚類中心，即

式中：z表示單個(gè)樣本向量；|Qi|表示簇的大小。

步驟5將點(diǎn)集P中全部點(diǎn)根據(jù)與聚類中心歐氏距離最小原則進(jìn)行劃分。

步驟6計(jì)算標(biāo)準(zhǔn)測度函數(shù)，當(dāng)滿足收斂條件時(shí)，函數(shù)停止迭代，算法終止；如不滿足則重復(fù)步驟3～步驟5。

3.4 選線流程

選線流程如圖5所示，具體步驟如下。

圖5 故障選線流程Fig.5 Flaw chart of fault line selection

步驟1利用檢測裝置持續(xù)監(jiān)測系統(tǒng)母線零序電壓U0是否超過額定電壓UN的115%，符合條件即啟動故障選線裝置，提取各出線在故障發(fā)生時(shí)刻后2個(gè)信號周期（40 ms）內(nèi)的零序電流信號波形。

步驟2將提取到的各出線上的零序電流信號用VMD 算法按照前文所述方法先過濾后再進(jìn)行3層分解，得到預(yù)設(shè)個(gè)數(shù)的本征模態(tài)函數(shù)。

步驟3根據(jù)式（13）和（16），通過對故障發(fā)生時(shí)刻后首個(gè)信號周期（20 ms）內(nèi)的低頻分量做相關(guān)性計(jì)算，對高頻分量測算初始極性，進(jìn)而確定各點(diǎn)橫、縱坐標(biāo)，作出二維平面散點(diǎn)分布圖。

步驟4使用K-means++算法進(jìn)行聚類分析，篩選出故障線路。

4 仿真分析

4.1 仿真模型

本節(jié)將在仿真的基礎(chǔ)上分析驗(yàn)證本文所提算法的實(shí)用性。通過仿真軟件Pscad 建立配電網(wǎng)模型，獲得各出線單相接地條件下的零序電流信號波形；將零序電流信號導(dǎo)入Matlab工作區(qū)中進(jìn)行VMD算法分解，獲得故障信號不同頻段下的模態(tài)分量；使用本文所提出的方法進(jìn)行模擬選線，得到最終結(jié)果數(shù)據(jù)。

在Pscad中搭建一個(gè)6條出線6個(gè)負(fù)荷的10 kV配電網(wǎng)模型，其結(jié)構(gòu)如圖6所示。系統(tǒng)中，線路L1、L3、L4、L6為架空線路，除L1長度為15 km之外，其余長度均為12 km；線路L2為電纜線路，長度8 km；線路L5為纜線混合線路，電纜長4 km，架空線長12 km。各出線相關(guān)參數(shù)如表2所示，各負(fù)荷功率均設(shè)為2 MW，消弧線圈電感值采用過補(bǔ)償10%進(jìn)行設(shè)置。

表2 線路參數(shù)Tab.2 Parameters of lines

圖6 仿真電路模型Fig.6 Simulation circuit model

4.2 不同初始條件下的仿真

仿真實(shí)驗(yàn)選定的3 個(gè)變量分別是故障距離、故障初相角和接地過渡電阻。首先選取一個(gè)典型案例進(jìn)行詳細(xì)分析，隨后在選定兩個(gè)變量而不斷對剩余的單一變量做出改變的情況下進(jìn)行仿真實(shí)驗(yàn)，以驗(yàn)證本文所提出選線方法的適用性。

4.2.1 典型案例分析

假設(shè)距離出線L2始端4 km處發(fā)生A相單相接地短路故障，故障時(shí)刻為0.04 s，故障初相角設(shè)置為90°，接地過渡電阻大小設(shè)置為2 kΩ，通過故障錄波裝置提取出各出線上故障發(fā)生時(shí)刻后2 個(gè)信號周期（40 ms）的零序電流波形如圖7所示。

圖7 零序電流波形Fig.7 Zero-sequence current waveforms

將提取到的電流進(jìn)行濾波處理后，按照前文所述方式進(jìn)行3層VMD分解，提取故障條件下各出線零序電流分解結(jié)果中的IMF1-1分量和IMF2分量分別如圖8和圖9所示。

圖8 各出線IMF1-1波形Fig.8 IMF1-1 waveform of each line

圖9 各線路IMF2波形Fig.9 IMF2 waveform of each line

觀察圖8和圖9中各出線上零序電流分解后得到的低頻分量IMF1-1和高頻分量IMF2可知，在故障初相角為90°的情況下，故障線路與非故障線路低頻分量之間的波形差異不大，說明故障線路L2上流過的零序電流低頻分量中幾乎不含衰減直流分量；除故障線路L2的高頻分量初始極性為正外，其余各線路中零序電流高頻分量的初始極性全部為負(fù)。

算例的分解結(jié)果與本文第1 節(jié)中關(guān)于零序電流的成分及特性分析相吻合。

計(jì)算各出線IMF1-1分量之間的相關(guān)系數(shù)，所得結(jié)果如表3所示。

表3 各線路IMF1-1之間的相關(guān)系數(shù)Tab.3 Correlation coefficient between IMF1-1 of each line

計(jì)算各出線IMF2分量的初始極性，所得結(jié)果如表4所示。

表4 各線路IMF2的初始極性Tab.4 Initial polarity of IMF2 of each line

計(jì)算代表各出線的散點(diǎn)的橫、縱坐標(biāo)，并做出所有散點(diǎn)在二維平面分布，結(jié)果如圖10所示。

圖10 二維平面的散點(diǎn)分布Fig.10 Scatter distribution on two-dimensional plane

從圖中可以看到各散點(diǎn)的橫坐標(biāo)值均較為接近，代表僅通過低頻分量之間的波形相關(guān)性并不足以使故障線路的散點(diǎn)成為離群點(diǎn)，必須通過縱坐標(biāo)代表的高頻分量初始極性之間的差異構(gòu)造離群點(diǎn)。由于初相角φ=90°，高頻分量初始極性差異較為明顯，保證了故障線路與健全線路散點(diǎn)的縱坐標(biāo)值為相反數(shù)，因此故障線路L2的點(diǎn)在平面中仍然成為了一個(gè)離群點(diǎn)。

使用K-means++聚類算法對圖10 中的點(diǎn)集進(jìn)行聚類分析，得到輸出聚類標(biāo)簽為[2 1 2 2 2 2]，表示出線L2與其他線路不同，單獨(dú)劃分為了一類，因此可以得出結(jié)論為出線L2發(fā)生了單相接地故障。

4.2.2 適用性分析

選取故障初相角、故障距離以及過渡電阻為3個(gè)變量，在選定兩個(gè)變量而不斷對剩余的單一變量做出改變的情況下進(jìn)行仿真實(shí)驗(yàn)。受限于篇幅不再詳細(xì)分析，仿真結(jié)果如表5所示。

表5 仿真結(jié)果Tab.5 Results of simulation

從表5 中可以看出，不論故障初始條件如何改變，包括過渡電阻超過2 kΩ 等對選線不利的工況下，本文提出的算法依然能夠準(zhǔn)確地選擇出故障線路。除表5中列舉的情況，還進(jìn)行了大量的仿真實(shí)驗(yàn)，包括改變噪聲強(qiáng)度等，本文提出的算法均能正確地選出故障線路，表明該選線方法不受條件改變的影響，具有較強(qiáng)的適用性。

4.2.3 方法優(yōu)越性比較

為驗(yàn)證本文所提出的方法的優(yōu)越性，現(xiàn)假設(shè)在距離出線L5始端4 km 處發(fā)生單相接地故障，故障初相角設(shè)為60°，過渡電阻依次設(shè)置為100 Ω、500 Ω、1 000 Ω。分別用本文的方法和進(jìn)行EMD分解再使用K-means++聚類（方法2）以及進(jìn)行VMD分解再使用K-means聚類（方法3）進(jìn)行比較，結(jié)果如表6所示。

表6 不同方案的選線結(jié)果Tab.6 Line selection results under different schemes

通過對比3 種方法的結(jié)果可以看出，隨著過渡電阻的增大，高頻分量幅值的減小，本文的算法均正確地選擇出了故障線路，而方法2 有錯(cuò)誤記錄。方法3 雖然也能正確選線，但是相較本文算法，Kmeans 算法在執(zhí)行聚類任務(wù)時(shí)運(yùn)行更加耗時(shí)（約為40 ms），約為K-means++算法（約為10 ms）的4 倍，因此本文設(shè)計(jì)的選線方法優(yōu)于基于EMD分解或Kmeans聚類的相關(guān)性分析方法。

5 結(jié) 論

通過對諧振接地系統(tǒng)單相接地短路故障時(shí)故障零序電流組成特性的分析，本文提出了一種基于VMD 分解和改進(jìn)K-means++相關(guān)性聚類的故障選線方法，實(shí)現(xiàn)了精確故障選線的目的。主要結(jié)論如下。

（1）本文將VMD分解和K-means++聚類算法相結(jié)合，利用在構(gòu)造的散點(diǎn)圖中，故障線路信號所屬的散點(diǎn)屬于離群點(diǎn)的特點(diǎn)，使用由K-means 聚類算法改進(jìn)而來的K-means++聚類算法進(jìn)行分析，相比于K-means算法，提高了選線速率。

（2）本文將故障零序電流信號不同模態(tài)分量之間的差異投射到二維平面中，從兩個(gè)維度進(jìn)行聚類分析，極大提高了選線準(zhǔn)確率和保護(hù)裕度。

（3）加入噪聲進(jìn)行仿真實(shí)驗(yàn)，結(jié)果表明本文方法對噪聲具有一定的魯棒性。