金濤　許立彬　高偉　郭謀發(fā)　陳永往

摘要：針對當(dāng)前電力系統(tǒng)配網(wǎng)領(lǐng)域的過電壓識別問題，提出一種基于時(shí)頻分布特征的配電網(wǎng)暫時(shí)過電壓分類識別方法。在暫時(shí)過電壓中，單相金屬性接地過電壓的三相電壓能量分布均勻度最低，間歇性弧光接地零序電壓的直流分量最高，鐵磁諧振零序電壓信號在幅值和頻率集中頻帶上存在很大差異。計(jì)算零序電壓能量貢獻(xiàn)率、零序電壓平均值和三相電壓奇異譜熵，提取過電壓信號的時(shí)域能量分布特征；采用局部特征尺度分解（LCD）和Hilbert變換結(jié)合帶通濾波算法計(jì)算零序電壓重心頻帶，提取過電壓信號的頻域能量分布特征；結(jié)合閾值判別法實(shí)現(xiàn)暫時(shí)過電壓類型識別。該方法不需要分類器，算法簡單，計(jì)算時(shí)間少。仿真和實(shí)驗(yàn)表明該識別方法在不同故障條件下均有更高的識別率。

關(guān)鍵詞：中性點(diǎn)不接地系統(tǒng)；暫時(shí)過電壓；能量貢獻(xiàn)率；平均值；奇異譜熵；重心頻帶

DOI：10.15938/j.emc.2018.11.000

中圖分類號：TM 73

文獻(xiàn)標(biāo)志碼：A

文章編號：1007-449X（2018）11-0000-00

0引言

運(yùn)行經(jīng)驗(yàn)表明，過電壓是影響配電網(wǎng)安全運(yùn)行的重要因素之一[1-2]。暫時(shí)過電壓持續(xù)時(shí)間較長，易引起設(shè)備絕緣損壞，從而引發(fā)各類短路故障，危及配電網(wǎng)供電可靠性。因此，及時(shí)檢測配電網(wǎng)出現(xiàn)的暫時(shí)過電壓，準(zhǔn)確區(qū)分過電壓類型，對于配電網(wǎng)災(zāi)害預(yù)防和故障分析具有重要意義。

國內(nèi)外普遍采用的過電壓識別方法是先提取特征量，再結(jié)合支持向量機(jī)、極端學(xué)習(xí)機(jī)、徑向基核函數(shù)（radial basis function， RBF）神經(jīng)網(wǎng)絡(luò)等分類算法，得到識別結(jié)果。文獻(xiàn)[3]采用Morlet小波方法進(jìn)行故障特征提取，但對于不同種類的過電壓，小波方法存在最優(yōu)小波基選擇困難的問題。文獻(xiàn)[4]采用S變換進(jìn)行信號分解，通過核主成分分析法提取鐵磁諧振過電壓的特征量結(jié)合支持向量機(jī)進(jìn)行鐵磁諧振的識別。文獻(xiàn)[5]采用經(jīng)驗(yàn)?zāi)B(tài)分解提取電弧故障特征，并構(gòu)建極端學(xué)習(xí)機(jī)，用于配網(wǎng)不同負(fù)載電弧故障的識別，但經(jīng)驗(yàn)?zāi)B(tài)分解存在模態(tài)混疊和端點(diǎn)效應(yīng)的問題，對于包含不同時(shí)間尺度的過電壓信號，分解后的分量不唯一，該方法的適應(yīng)性有待提高。文獻(xiàn)[6]和文獻(xiàn)[7]采用RBF神經(jīng)網(wǎng)絡(luò)作為分類方法，神經(jīng)網(wǎng)絡(luò)對樣本數(shù)量要求較高，且RBF神經(jīng)網(wǎng)絡(luò)存在局部極小值問題。文獻(xiàn)[8]建立過電壓分層識別系統(tǒng)，采用多級支持向量機(jī)對特征量逐級選擇，應(yīng)用于暫時(shí)過電壓類型識別，但多級支持向量機(jī)需要大量訓(xùn)練樣本，且訓(xùn)練周期長。文獻(xiàn)[9]結(jié)合傅里葉變換和小波變換提取三類比值作為特征參量對鐵磁諧振過電壓信號進(jìn)行分類識別。文獻(xiàn)[10]利用希爾伯特黃變換（HilbertHuang transform， HHT）作為過電壓信號的時(shí)頻分析工具，基于瞬時(shí)幅值譜、Hilbert邊際譜和Hilbert時(shí)頻譜提取特征量，對于某些過電壓信號，分解時(shí)同樣會出現(xiàn)端點(diǎn)效應(yīng)、模態(tài)混疊的問題。文獻(xiàn)[11]和文獻(xiàn)[12]采用數(shù)學(xué)形態(tài)學(xué)濾波算子提取零序過電壓信號波形的突變特征，作為識別弧光接地過電壓的判據(jù)，所提方法較為新穎，但不同的信號采樣頻率和干擾信號會影響結(jié)構(gòu)元素的選取，進(jìn)而影響信號處理結(jié)果[13]。文獻(xiàn)[14]基于S變換模矩陣奇異值結(jié)合分類器識別操作及雷電過電壓，在高頻瞬態(tài)過電壓信號的識別上具有一定的參考價(jià)值。局部特征尺度分解（local characteristicscale decomposition，LCD）是一種新的非線性、非平穩(wěn)信號的自適應(yīng)時(shí)頻分析方法。該方法繼承了經(jīng)驗(yàn)?zāi)B(tài)分解法（empirical mode decomposition，EMD）的優(yōu)點(diǎn)，并改善了EMD存在的端點(diǎn)效應(yīng)、模態(tài)混疊和分解速度等問題，已廣泛應(yīng)用于機(jī)械故障診斷領(lǐng)域[15]。

本文采用LCDHilbert變換結(jié)合帶通濾波算法提取零序電壓信號的頻域分布特征，將其與三相電壓、零序電壓的時(shí)域分布特征有機(jī)結(jié)合，采用閾值判別法對配電網(wǎng)暫時(shí)過電壓進(jìn)行類型識別。該方法不需要使用分類器，可有效克服樣本獲取困難、耗時(shí)長等問題，且在迭代次數(shù)上優(yōu)于EMD，用于暫時(shí)過電壓的實(shí)時(shí)類型識別具有良好的效果。

1配電網(wǎng)暫時(shí)過電壓特性分析

配電網(wǎng)暫時(shí)過電壓包括諧振過電壓和工頻電壓升高?；诖?，著重研究中性點(diǎn)不接地系統(tǒng)中出現(xiàn)的對絕緣危害較大的不對稱接地故障引起的工頻電壓升高和電磁式電壓互感器飽和引起的鐵磁諧振過電壓。下面分析過渡電阻阻值大小對不對稱接地故障過電壓倍數(shù)的影響。圖1為單相接地故障的復(fù)合序網(wǎng)圖，Rf為過渡電阻，Uf|0|為f點(diǎn)正常時(shí)的電壓；Z∑（1）、Z∑（2）、Z∑（0）分別對應(yīng)正序網(wǎng)、負(fù)序網(wǎng)、零序網(wǎng)f點(diǎn)的等值阻抗。

設(shè)故障相為A相，EA=EA∠0°為故障前A相的正常電壓，則EA=Uf|0|，計(jì)算故障后A相電壓UA、B相電壓UB為：

UA=3RfZ∑（1）+Z∑（2）+Z∑（0）+3RfEA，（1）

UB=a2Uf（1）+aUf（2）+Uf（0）=

a2（EA-Z∑（1）If（1））+a（-Z∑（2）If（2））+

（-Z∑（0）If（0））。（2）

式中：a為旋轉(zhuǎn)因子，a=1∠120°，a2=1∠240°，a2+a+1=0。X∑（1）、X∑（2）和X∑（0）分別為正序、負(fù)序和零序網(wǎng)的等效電抗。因?yàn)樵趯?shí)際系統(tǒng)中，X∑（i）R∑（i）（i=1，2，0），因此忽略系統(tǒng)的三序電阻，則X∑（i）≈Z∑（i）（i=1，2，0），且X∑（1）≈X∑（2），If（1）=If（2）=If（0）=Uf|0|Z∑（1）+Z∑（2）+Z∑（0）+3Rf，則式（2）可變?yōu)?/p>

UB=a2EA-j（X∑（0）-X∑（1））If（1）=

a2EA-j（X∑（0）-X∑（1））EAj（2X∑（1）+X∑（0））+3Rf。（3）

令k0=X∑（0）X∑（1），則式（3）可變?yōu)?/p>

UB=EB-j（k0-1）EAj（k0+2）+3RfX∑（1）=

EB-1-j3+3RfX∑（1）j（k0+2）+3RfX∑（1）EA。（4）

當(dāng)Rf=0時(shí)，UA=0，由于X∑（0）X∑（1），因此k0=

SymboleB@ ，UB=EB-EA=3EB∠-30°；當(dāng)Rf阻值增大至與零序電抗X∑（0）數(shù)量級相同時(shí)，kf=RfX∑（1）=∞，則式（4）可變?yōu)?/p>

UB=EB-1-3jk0kf+3EA。（5）

Rf越大，kf越大，|UB|越接近|EB|，因此在高阻接地故障時(shí)非故障相電壓UB的過電壓倍數(shù)很低，故僅考慮單相金屬性接地的引起的工頻電壓升高。

圖2和圖3是典型的過電壓三相電壓和零序電壓時(shí)域波形，過電壓發(fā)生時(shí)刻均為0.01 s，采樣頻率為20 kHz，以故障前半個(gè)工頻周波為起始時(shí)刻，給出5個(gè)工頻周期的過電壓波形示例。

單相金屬性接地過電壓發(fā)生時(shí)，如圖2（a）所示，故障相電壓幅值跌落至0，非故障相電壓幅值升高至原線電壓幅值，因此相電壓之間能量分布相差很大，盡管故障點(diǎn)、故障時(shí)刻和過渡電阻會發(fā)生變化，這個(gè)特征仍保持不變，因此可作為識別單相金屬性接地過電壓的依據(jù)?；」饨拥剡^電壓的三相波形如圖2（b）所示，當(dāng)C相電壓達(dá)到最大值時(shí)（0.01 s）發(fā)生第一次絕緣擊穿，故障點(diǎn)的電弧在工頻電流通過第一個(gè)過零點(diǎn)時(shí)熄滅，此時(shí)，非故障相上的自由電荷將沿三相對地電容重新分布，于是在各相上便產(chǎn)生了同等的位移電壓，熄弧時(shí)段三相電壓恢復(fù)為正弦電壓。給出其零序電壓波形如圖2（d）所示，由于位移電壓的存在，零序電壓每隔半個(gè)周波出現(xiàn)一次直流分量，可作為識別弧光接地過電壓的時(shí)域特征。

鐵磁諧振過電壓以單相故障消失作為“激發(fā)”條件，0.01 s時(shí)刻單相接地故障消失，原故障相電壓振幅逐漸增大，三相電壓的能量分布較為均勻，但不同類型的諧振過電壓振蕩頻率不同（如圖3（a）、圖3（b）、圖3（e）所示）：分頻諧振可能包含1/2次、1/3次、1/5次等分次諧波含量，振蕩頻率為工頻的分?jǐn)?shù)倍；高頻諧振的振蕩頻率為工頻的整數(shù)次倍；中性點(diǎn)位移電壓屬于基頻零序電壓時(shí)，產(chǎn)生基頻諧振，顯然其基波含量占比最大。鐵磁諧振的振蕩頻率不同必然會反映在零序電壓上，為提高識別速率，可將零序電壓能量集中的頻帶作為區(qū)分鐵磁諧振過電壓的特征量。

2配電網(wǎng)暫時(shí)過電壓識別方法分析

根據(jù)以上分析，提出一種基于LCDHilbert變換和奇異譜熵的過電壓識別方法，對典型的暫時(shí)過電壓進(jìn)行分層識別。如圖4所示，給出涉及到的幾類過電壓的分層結(jié)構(gòu)圖。間歇性弧光過電壓雖屬操作過電壓，但因電弧間歇性熄滅和重燃，持續(xù)時(shí)間較長，故將間歇性弧光接地過電壓歸為暫時(shí)過電壓一起分析。

首先，基于零序電壓在時(shí)域上的能量分布差異，采用能量貢獻(xiàn)率作為特征量，排除由斷路器操作引發(fā)的過電壓，從而獲得本文的識別對象，如圖4虛線框內(nèi)所示；第二，基于三相電壓能量分布差異，采用奇異譜熵作為特征量，識別單相金屬性接地；第三，基于零序電壓直流分量的高低，采用平均值的絕對值作為特征量，識別間歇性弧光接地，同時(shí)，基于鐵磁諧振零序電壓能量集中頻帶的不同，采用重心頻帶作為特征量，識別高頻諧振、基頻諧振和分頻諧振。下面介紹地方法的原理。

2.1時(shí)域特征量提取

在分析過程中，操作過電壓是配電網(wǎng)中較為常見的故障，在分析過程中需要先行排除。操作過電壓持續(xù)時(shí)間較短，一般小于10 ms，暫時(shí)過電壓持續(xù)時(shí)間較長，從幾個(gè)工頻周期到1.5～2 h不等。在所截取的時(shí)段內(nèi)，操作過電壓的能量分布多集中于故障后0.5個(gè)周波，而暫時(shí)過電壓能量分布較均勻，因此，暫時(shí)過電壓和操作過電壓的零序電壓信號在故障后的0.5個(gè)周波的能量與所截取時(shí)段的總能量的比值相差較大。定義能量貢獻(xiàn)率E為

E=∑N1k=1v0（k）2/∑N2k=1v0（k）2。（6）

式中：N1為0.5個(gè)周波的采樣點(diǎn)數(shù)；N2為截取的總時(shí)段的采樣點(diǎn)數(shù)；v0（k）為零序電壓信號序列。

“熵”最早出現(xiàn)在熱力學(xué)研究中，是用于衡量一個(gè)系統(tǒng)中分子狀態(tài)的均勻狀況，分子狀態(tài)越均勻，熵值越高，當(dāng)系統(tǒng)中分子分布最均勻時(shí)，熵值達(dá)到最大[16]?；谙嚯妷褐g的能量分布差異，由于接地瞬間單相金屬性接地三相電壓可能有高頻振蕩，因此以故障發(fā)生后2個(gè)工頻周波為起始點(diǎn)，截取5個(gè)工頻周波的三相電壓信號，計(jì)算其奇異譜熵作為識別單相金屬性接地過電壓的特征量。

設(shè)截取的信號采樣點(diǎn)數(shù)為n，三相電壓信號矩陣為U=[ua，ub，uc]T，對矩陣U的每一行作奇異值分解，可得奇異譜Λ=diag（μ1，μ2，μ3），三相電壓信號奇異譜熵定義為

S=-∑3i=1qilogqi。（7）

式中：qi=μi/∑3i=1μi（i=1，2，3）表征了μi對整個(gè)奇異譜的貢獻(xiàn)率[17]；S可有效表征三相電壓信號奇異值分布的不均勻度，進(jìn)而反映三相電壓之間幅值分布的差異。

間歇性弧光接地的零序電壓包含較高的直流分量，采用零序電壓信號的平均值絕對值表征這個(gè)特性。零序電壓信號的平均值絕對值M定義為

M=|∑nk=1v0（k）|n。（8）

式中：v0（k）是各種過電壓零序電壓信號；n為故障時(shí)刻開始5個(gè)周波的零序電壓信號的采樣點(diǎn)數(shù)。

2.2LCDHILBERT變換

LCDHilbert變換由兩部分組成：首先，通過LCD將原信號分解成有限個(gè)內(nèi)稟尺度分量（intrinsic scale component，ISC）；其次，對各個(gè)內(nèi)稟尺度分量進(jìn)行Hilbert變換，進(jìn)而得到瞬時(shí)頻率函數(shù)。

LCD分解計(jì)算得到的ISC應(yīng)滿足以下條件：1）信號數(shù)據(jù)的任意兩個(gè)相鄰極值點(diǎn)符號互異；2）信號數(shù)據(jù)的任意兩個(gè)相鄰極大（或小）值點(diǎn)確定的直線在二者之間的極?。ɑ虼螅┲迭c(diǎn)所對應(yīng)的橫坐標(biāo)處的函數(shù)值與該極?。ɑ虼螅┲迭c(diǎn)的比值保持不變[14]。

根據(jù)上述兩個(gè)條件，將過電壓信號進(jìn)行局部特征尺度分解，分解成若干個(gè)ISC和一個(gè)殘余項(xiàng)之和為

x（t）=∑ni=1ISCi（t）+r（t）。（9）

式中：x（t）為原始信號；n為ISC的個(gè)數(shù)；ISCi（t）為第i個(gè)ISC；r（t）為剩余分量?；l諧振的零序電壓信號如圖1所示，其經(jīng)LCD分解后的ISC如圖5所示。

2.3頻域特征量提取

利用LCDHilbert變換結(jié)合帶通濾波算法可實(shí)現(xiàn)過電壓信號在頻帶上的分離并濾去高頻噪聲干擾，增強(qiáng)信號可比性，提高識別成功率。帶通濾波算法的基本原理是：對故障波形進(jìn)行LCDHilbert變換，得到ISC分量和瞬時(shí)頻率函數(shù)之后，根據(jù)瞬時(shí)頻率范圍選擇合適的頻帶寬度Δf，以此進(jìn)行頻帶劃分，則第i個(gè)頻帶的頻率范圍為[（i-1）Δf ，iΔf]，將所有ISC分量不在該頻率范圍內(nèi)的瞬時(shí)能量點(diǎn)都置為零，其余在該頻率范圍內(nèi)的瞬時(shí)能量點(diǎn)保持不變，即可得到所有ISC分量在該頻帶的成分。最后將經(jīng)第i個(gè)頻帶濾波后的所有ISC分量疊加，可得到原始信號在該頻帶內(nèi)的分量，同理可得到其他頻帶內(nèi)的分量。

不同頻帶重構(gòu)波形的能量大小可反映過電壓信號的頻域分布情況。以各頻帶重構(gòu)信號的平方和表征該頻帶信號能量，第i頻帶的能量為

ei=∑mk=1c2i（k）。（15）

式中：ei為第i頻帶的信號能量；m為重構(gòu)信號的數(shù)據(jù)長度；ci（k）為第i頻帶的重構(gòu)信號。

物理學(xué)中用物理重心來反映物體受重力作用所集中到的點(diǎn)，將重心的概念應(yīng)用到頻域分析上，選取各個(gè)頻帶中能量最大的頻帶作為頻域的“重心”，可反映波形諧波含量的集中位置。定義重心頻帶為

Ng=imax。（16）

式中：Ng為重心頻帶；imax為能量最大的頻帶編號。

2.4方法的計(jì)算流程

確定了特征量E、S、M和Ng后，針對中性點(diǎn)不接地系統(tǒng)，給出配電網(wǎng)暫時(shí)過電壓類型的識別流程如圖6所示。暫時(shí)過電壓類型識別步驟如下：

1）獲取配電網(wǎng)母線處的三相電壓和零序電壓信號，對三相電壓序列求取奇異譜熵S，判斷其是否小于閾值，若是，判為單相金屬性接地過電壓；

2）若奇異譜熵S大于設(shè)定的閾值，求取零序電壓序列平均值的絕對值M，判斷其是否大于閾值，若是，判為間歇性弧光接地過電壓；

3）若M小于設(shè)定的閾值，采用LCDHilbert變換結(jié)合帶通濾波算法重構(gòu)相應(yīng)的零序電壓波形，求取重心頻帶Ng，按照劃定好的范圍進(jìn)行分類。

3仿真實(shí)驗(yàn)分析和討論

根據(jù)以上分析，采用ATPEMTP電磁暫態(tài)仿真軟件搭建中性點(diǎn)不接地配網(wǎng)模型，如圖7所示。其中：G為110 kV無窮大系統(tǒng)電源；T1為110 kV/10.5 kV主變壓器，T2為10 kV/0.4 kV配電變壓器；負(fù)荷采用感性負(fù)載R+jX=（20+j40） Ω等值替代。饋線共有全架空線路（OL）、全電纜線路（CL）和線-纜混合線路3種類型，線路參數(shù)如表1所示。K1為時(shí)控開關(guān)，用于模擬接地故障，Rf為過渡電阻。母線PT為電磁式電壓互感器，以“F”開頭的編號均為仿真實(shí)驗(yàn)設(shè)置的故障點(diǎn)。

根據(jù)工頻熄弧理論分析所得的弧光接地過電壓倍數(shù)較接近配電網(wǎng)出現(xiàn)的弧光接地的過電壓倍數(shù)，該理論認(rèn)為故障相每隔一個(gè)工頻周期電弧就會熄滅和重燃一次，并且假定電弧在工頻電流過零時(shí)熄滅，達(dá)到最大恢復(fù)電壓時(shí)重燃[12]，因此基于該理論，通過周期性投切時(shí)控開關(guān)K1進(jìn)行弧光接地過電壓的模擬。采用單相接地故障消失作為鐵磁諧振的“激發(fā)”條件，為了實(shí)現(xiàn)鐵磁諧振的參數(shù)匹配要求[18]，通過投切不同長度的饋線形成不同的匹配電容，達(dá)到改變諧振回路三相導(dǎo)納的目的，以產(chǎn)生基頻諧振、分頻諧振和高頻諧振。圖7中，投入饋線L3、饋線L6進(jìn)行分頻諧振的仿真，投入饋線L5、饋線L6進(jìn)行基頻諧振的仿真，投入饋線L6進(jìn)行高頻諧振的仿真。其中線路的參數(shù)如表1所示。

將過渡電阻設(shè)置為1、2、5、10、20 Ω，故障初相角選擇±90°（對應(yīng)正負(fù)電壓峰值時(shí)刻，若分析對象為鐵磁諧振，則是指單相接地故障消除時(shí)刻初始相角選擇±90°），每種過電壓記錄10組波形，計(jì)算其奇異譜熵。如圖9所示，每組數(shù)據(jù)的1～5個(gè)點(diǎn)對應(yīng)-90°初相角、過渡電阻分別為1～20 Ω時(shí)的奇異譜熵，第6～10個(gè)點(diǎn)對應(yīng)90°初相角、過渡電阻分別為1～20 Ω時(shí)的奇異譜熵?？梢钥闯觯簡蜗嘟饘傩越拥剡^電壓的奇異譜熵最低，均小于1.2；其余四類過電壓的奇異譜熵均大于1.5，隨著過渡電阻的增大，單相金屬性接地的奇異譜熵呈上升趨勢。計(jì)算在20 Ω過渡電阻情況下，不同故障點(diǎn)、不同故障初相角的100組三相電壓數(shù)據(jù)的奇異譜熵，判別單相金屬性接地過電壓的閾值可取1.4。

弧光接地熄弧時(shí)段零序電壓存在不為零直流分量，計(jì)算弧光接地與鐵磁諧振的零序電壓平均值如圖10所示。可以看出，弧光接地的零序電壓平均值的絕對值最大。因此，采用零序電壓平均值的絕對值M作為識別弧光接地過電壓的特征量，計(jì)算不同過渡電阻、不同故障點(diǎn)的100組零序電壓數(shù)據(jù)的M值，確定其閾值為2 000。

以鐵磁諧振的零序電壓作為分析對象，采用帶通濾波算法對LCD分解獲得的各個(gè)ISC依據(jù)Hilbert瞬時(shí)頻率函數(shù)分頻帶多支重構(gòu)：通過式（10）～式（14）計(jì)算不同過渡電阻、不同故障初相角等因素影響下的5個(gè)工頻周波的零序電壓波形的瞬時(shí)頻率向量，共200組數(shù)據(jù)的仿真結(jié)果，發(fā)現(xiàn)信號的瞬時(shí)頻率集中在0～600 Hz范圍內(nèi)，因此將此范圍作為帶通濾波算法的總頻帶寬度；將總頻帶劃分為6個(gè)頻帶：0～10 Hz、10～20 Hz、20～30 Hz、30～40 Hz、40～80 Hz、80～600 Hz，對各頻帶依次編號，0～10 Hz分量為頻帶1，10～20 Hz為頻帶2，以此類推；設(shè)當(dāng)前濾波頻帶為i，若瞬時(shí)頻率不在頻帶頻率范圍內(nèi)，則將該點(diǎn)對應(yīng)的ISC瞬時(shí)能量置零，最后將經(jīng)第i頻帶濾波后的所有ISC疊加，可得原始信號在該頻帶內(nèi)的分量，同理可得到其他頻帶內(nèi)的分量。基頻諧振零序電壓重構(gòu)后的波形如圖11所示，幅值最高的重構(gòu)信號出現(xiàn)在基波頻帶，與基頻諧振的原理分析一致。

計(jì)算不同過渡電阻、不同故障初相角影響下鐵磁諧振零序電壓的重心頻帶，如表2所示。可以看出：分頻諧振重心頻帶為2或3，這是因?yàn)榉诸l諧振1/2次或1/3次諧波含量最高；基頻諧振重心頻帶為5，這是因?yàn)榛l諧振的基波含量最高；高頻諧振重心頻帶為6，這是因?yàn)楦哳l諧振的高次諧波含量最高。因此重心頻帶Ng可作為識別分頻諧振、基頻諧振、高頻諧振的特征量。判別高頻諧振的Ng取為6，判別基頻諧振的Ng取為5，判別分頻諧振的Ng閾值取為4。

通過改變故障點(diǎn)、故障相別、故障初相角、過渡電阻4個(gè)因素，對圖7所示的中性點(diǎn)不接地系統(tǒng)發(fā)生的暫時(shí)過電壓進(jìn)行仿真，以驗(yàn)證識別方法的有效性，仿真數(shù)據(jù)的實(shí)驗(yàn)條件如表3所示。同時(shí)給出其余仿真數(shù)據(jù)的識別統(tǒng)計(jì)結(jié)果，如表4所示。將LCD應(yīng)用于所提方法對于暫時(shí)過電壓類型識別的正確率可達(dá)到94.98%，將EMD應(yīng)用于所提方法對暫時(shí)過電壓類型識別也是有效的，準(zhǔn)確率達(dá)92%以上，但是低于LCD方法，原因在于其混疊效應(yīng)較為嚴(yán)重，而將頻帶范圍為40～80 Hz的基頻諧振誤識別為30～40 Hz的分頻諧振。采用LCD方法的識別速度總體優(yōu)于EMD方法。所提方法不僅可以實(shí)現(xiàn)鐵磁諧振的分類，還可獲取諧振頻率范圍，對于鐵磁諧振的抑制具有參考價(jià)值。

為了模擬現(xiàn)場電磁環(huán)境的噪聲干擾，抽取表4中已正確識別的若干組過電壓數(shù)據(jù)加入20 dB的高斯白噪聲，驗(yàn)證噪聲對識別方法的影響。識別結(jié)果如表5所示，識別正確率為100%，這是由于帶通濾波算法本身具有良好的濾波效果，因此噪聲對于所提方法的影響十分微弱，說明了所提方法具有良好的抗噪性能。

本文識別方法在表4所示的1 434組數(shù)據(jù)中的平均識別時(shí)間是0.13 s。在實(shí)驗(yàn)室進(jìn)行算法的工程驗(yàn)證，采用數(shù)據(jù)采集卡作為在線采集終端，在福州大學(xué)配電網(wǎng)物理仿真系統(tǒng)[19]上進(jìn)行接地性實(shí)驗(yàn)，處理器是電腦主機(jī)，數(shù)據(jù)采集卡內(nèi)嵌在電腦主機(jī)中，識別算法可以在線識別。

在實(shí)際電力現(xiàn)場應(yīng)用中，已在福建省數(shù)個(gè)變電站安裝了10 kV過電壓在線采集終端，其實(shí)物圖如圖12所示，其主要功能是過電壓啟動后，將過電壓波形通過GPRS通信上傳到后臺主站，由主站負(fù)責(zé)離線分析，當(dāng)前還無法應(yīng)用于實(shí)時(shí)控制。

4結(jié)論

操作過電壓的持續(xù)時(shí)間較暫時(shí)過電壓短，單相金屬性接地過電壓的三相能量分布差異較其余四種暫時(shí)過電壓大，間歇性弧光接地的零序電壓包含最高的零序分量，鐵磁諧振過電壓的能量集中頻帶均不相同。采用零序電壓能量貢獻(xiàn)率、零序電壓平均值、三相電壓奇異譜熵作為時(shí)域特征量，通過LCDHilbert變換和帶通濾波算法重構(gòu)零序電壓波形，以重心頻帶作為頻域特征量，結(jié)合閾值判別法構(gòu)成暫時(shí)過電壓的類型識別方法。該方法所提取的特征量，綜合利用了零序電壓和三相電壓信號在頻域和時(shí)域的整體分布規(guī)律，完整地刻畫了過電壓信號特征，對不同的故障點(diǎn)、故障相別、故障過渡電阻、故障初相角情況下發(fā)生的暫時(shí)過電壓均能正確識別。

參 考 文 獻(xiàn)：

[HT6SS]

[1]黃艷玲，司馬文霞，楊慶，等. 電力系統(tǒng)實(shí)測過電壓信號的特征量提取與驗(yàn)證[J].高電壓技術(shù)，2013，39（1）：60.

HUANG Yanling， SIMA Wenxia， YANG Qing，et al. Features extraction and validation from measured overvoltage signals of power system[J]. High Voltage Engineering， 2013， 39（1）： 60.

[2]CERRETTI A， GATTA M F， GERI A， et al. Ground fault temporary overvoltages in MV networks： evaluation and experimental tests[J]. IEEE Transaction Power Delivery， 2012， 27（3）： 1592.

[3]侯新國， 牛超， 楊忠林. 基于最優(yōu)Morlet小波自適應(yīng)包絡(luò)解調(diào)的弱故障特征提取方法[J]. 電機(jī)與控制學(xué)報(bào)， 2016， 20（10）： 88.

HOU Xinguo， NIU Chao， YANG Zhonglin. Method to extract weak fault feature based on optimal Morlet wavelet adaptive envelope demodulation[J]. Electric Machines and Control， 2016， 20（10）： 88.

[4]MOKRYANI G， SIANO P， PICCOLO A. Identification of ferroresonance based on Stransform and support vector machine[J]. Simulation Modelling Practice and Theory， 2010， 18（9）： 1412.

[5]張麗萍， 繆希仁， 石敦義. 基于EMD和ELM的低壓電弧故障識別方法的研究[J]. 電機(jī)與控制學(xué)報(bào)， 2016， 20（9）： 54.

ZHANG Liping， MIAO Xiren， SHI Dunyi. Research on low voltage arc fault recognition method based on EMD and ELM [J]. Electric Machines and Control， 2016， 20（9）： 54.

[6]管春，周雒維，盧偉國.基于多標(biāo)簽RBF神經(jīng)網(wǎng)絡(luò)的電能質(zhì)量復(fù)合擾動分類方法[J].電工技術(shù)學(xué)報(bào)，2011（8）：198.

GUAN Chun， ZHOU Luowei， LU Weiguo. Recognition of multiple power quality disturbances using multilabel RBF neural networks [J].Transactions of China Electrotechnical Society， 2011（8）： 198.

[7]司馬文霞，王荊，楊慶，等.HilbertHuang變換在電力系統(tǒng)過電壓識別中的應(yīng)用[J].高電壓技術(shù)，2010，36（6）：1480.

SIMA Wenxia， WANG Jing， YANG Qing， et al. Application of HilbertHuang transform to power system overvoltage recognition[J]. High Voltage Engineering， 2010， 36（6）： 1480.

[8]杜林，李欣，王麗蓉，等.電力系統(tǒng)暫時(shí)過電壓多級支持向量機(jī)分層識別[J].電力系統(tǒng)保護(hù)與控制，2012，40（4）：26.

DU Lin， LI Xin， WANG Lirong， et al. Temporary overvoltage layered pattern identification based on multistage support vector machine in powersystem[J]. Power System Protection and Control， 2012， 40（4）： 26.

[9]杜林，李欣，吳高林，等. 采用3類特征參量比值法的鐵磁諧振過電壓識別[J].高電壓技術(shù)，2011，37（9）： 2241.

DU Lin，LI Xin， WU Gaolin， et al. Ferroresonance overvoltage identification using three feature parameters of ratio method[J]. High Voltage Engineering， 2011， 37（9）： 2241.

[10]WANG Jing， YANG Qing， CHEN Lin， et al. Classification and identification of overvoltage based on HHT and SVM[J]. High Voltage Engineering， 2012， 38（8）： 2068.

[11]司馬文霞，冉銳，袁濤，等.采用數(shù)學(xué)形態(tài)學(xué)的弧光接地過電壓識別方法[J].高電壓技術(shù)，2010，36（4）：835.

SIMA Wenxia， RAN Rui， YUAN Tao， et al. Identification of arc grounding overvoltage using mathematical morphology transform[J]. High Voltage Engineering， 2010， 36（4）： 835.

[12]羅寧昭，張曉鋒，劉佳.基于GradientClosing變換的接地故障識別[J].電力系統(tǒng)保護(hù)與控制，2011，39（23）：74.

LUO Ningzhao， ZHANG Xiaofeng， LIU Jia. Ground fault identification based on GradientClosing transform[J]. Power System Protection and Control， 2011， 39（23）： 74.

[13]王健民， 張春東. 電網(wǎng)諧波檢測的數(shù)學(xué)形態(tài)學(xué)優(yōu)化算法[J]. 哈爾濱理工大學(xué)學(xué)報(bào)， 2013， 18（2）：114.

WANG Jianmin， ZHANG Chundong. The mathematical morphology optimization algorithm of harmonic detection[J]. Journal of Harbin University of Science and Technology， 2013， 18（2）：114.

[14]杜林，李欣，司馬文霞，等. S變換模矩陣和最小二乘SVM在雷電及操作過電壓識別中的應(yīng)用[J].電力自動化設(shè)備，2012，32（8）：35.

DU Lin， LI Xin，SIMA Wenxia， et al. Application of Stransform modular matrix and LSSVM in identification of lightning and switching overvoltages[J]. Electric Power Automation Equipment， 2012，32（8）：35.

[15]程軍圣，鄭近德，楊宇.一種新的非平穩(wěn)信號分析方法——局部特征尺度分解法[J].振動工程學(xué)報(bào)，2012，25（2）： 215.

CHENG Junsheng， ZHENG Jinde， YANG Yu. A nonstationary signal analysis approach： the local characteristic scale decomposition method[J]. Journal of Vibration Engineering， 2012， 25（2）： 215.

[16]陳仕龍，曹蕊蕊，畢貴紅，等.利用多分辨奇異譜熵和支持向量機(jī)的特高壓直流輸電線路區(qū)內(nèi)外故障識別方法[J].電網(wǎng)技術(shù)，2015，39（4）： 991.

CHEN Shilong， CAO Ruirui， BI Guihong， et al. Distinguish faults located inside/outside protection zone of UHVDC transmission line by multiresolution singular spectrum entropy and support vector machine[J]. Power System Technology， 2015， 39（4）： 991.

[17]楊豐源，宋輝，程序，等.基于多分辨高階奇異譜熵分析的局部放電信號特征提取[J].電網(wǎng)技術(shù)，2016， 40（10）：3266.

YANG Fengyuan， SONG Hui， CHENG Xu， et al. Feature extraction of partial discharge signal based on multiresolution and highorder singular spectrum entropy analysis[J]. Power System Technology， 2016， 40（10）： 3266.

[18]司馬文霞，郭飛，楊慶，等.變電站鐵磁諧振仿真分析及抑制措施研究[J].電力自動化設(shè)備，2007，27（6）： 23.

SIMA Wenxia， GUO Fei， YANG Qing， et al. Simulative analysis of substation ferroresonance and its suppression measures[J]. Electric Power Automation Equipment， 2007， 27（6）： 23.

[19]郭謀發(fā)，游林旭，洪翠，等.基于LCDHilbert譜奇異值和多級支持向量機(jī)的配電網(wǎng)故障識別方法[J].高電壓技術(shù)，2017，43（4）：1246.

GUO Moufa， YOU Linxu， HONG Cui， et al. Identification method of distribution network faults based on singular value of LCDHilbert spectrums and multilevel SVM[J]. High Voltage Engineering， 2017， 43（4）：1246.

（編輯：邱赫男）