黃南天，張衛(wèi)輝，徐殿國，蔡國偉，劉 闖，張書鑫

（1.東北電力大學(xué)電氣工程學(xué)院，132012吉林吉林；2.哈爾濱工業(yè)大學(xué)電氣工程及自動(dòng)化學(xué)院，150001哈爾濱）

采用多分辨率廣義S變換的電能質(zhì)量擾動(dòng)識(shí)別

黃南天1，張衛(wèi)輝1，徐殿國2，蔡國偉1，劉 闖1，張書鑫1

（1.東北電力大學(xué)電氣工程學(xué)院，132012吉林吉林；2.哈爾濱工業(yè)大學(xué)電氣工程及自動(dòng)化學(xué)院，150001哈爾濱）

為提高電能質(zhì)量復(fù)合擾動(dòng)識(shí)別能力，提出一種采用多分辨率廣義S變換（multiresolution generalized S-transform，GST）的擾動(dòng)識(shí)別方法.首先，將信號(hào)頻譜分為低頻、中頻、高頻3個(gè)頻域，分別設(shè)定窗寬調(diào)整因子，使其在各個(gè)頻域具有不同的時(shí)-頻分辨率，滿足不同擾動(dòng)信號(hào)識(shí)別要求.并針對(duì)高頻振蕩識(shí)別問題，設(shè)計(jì)基于基頻傅里葉譜特征的自適應(yīng)窗寬調(diào)整方法.在此基礎(chǔ)上，提取6種特征用于構(gòu)建決策樹.最后，提出最小分類損失原則，確定決策樹節(jié)點(diǎn)分類閾值，設(shè)計(jì)擾動(dòng)分類器.仿真與實(shí)測(cè)信號(hào)實(shí)驗(yàn)證明，新方法能夠準(zhǔn)確識(shí)別含5種復(fù)合擾動(dòng)在內(nèi)的13種擾動(dòng).相較于S變換、廣義S變換和Hyperbolic S變換，新方法具有更好的特征表現(xiàn)能力，分類效果好，抗噪聲干擾能力強(qiáng).

電能質(zhì)量；電能質(zhì)量暫態(tài)擾動(dòng)；S變換；多分辨率；決策樹

電能質(zhì)量暫態(tài)擾動(dòng)分類是電能質(zhì)量分析與控制的重要基礎(chǔ)，對(duì)暫態(tài)治理、電力電子設(shè)備狀態(tài)監(jiān)控、擾動(dòng)源定位等工作具有重要意義［1］.隨著智能電網(wǎng)將“安全、經(jīng)濟(jì)、電能質(zhì)量”作為建設(shè)核心［2］，電能質(zhì)量的監(jiān)控與分析對(duì)擾動(dòng)識(shí)別的準(zhǔn)確性、高效性提出了更高的要求.常用的擾動(dòng)識(shí)別方法一般包括信號(hào)處理與模式識(shí)別兩個(gè)步驟［1］.電力系統(tǒng)發(fā)生的擾動(dòng)中含大量復(fù)合擾動(dòng)，這對(duì)信號(hào)處理提出了更高的要求［3］.同時(shí)，電能質(zhì)量監(jiān)測(cè)系統(tǒng)需要處理海量監(jiān)測(cè)數(shù)據(jù)，因此，擾動(dòng)識(shí)別方法需具有良好的分類效率.

現(xiàn)有信號(hào)處理方法多使用時(shí)-頻分析方法，包括希爾伯特-黃變換（Hilbert-Huang transform，HHT）［4］、小波變換（wavelet transform）［5］和 S變換（S-transform，ST）［6-7］等.S變換具有較好的抗噪性，采用的高斯窗函數(shù)具有可變的時(shí)-頻分辨率，能夠滿足不同頻域信號(hào)的分析要求；其結(jié)果為二維的時(shí)-頻矩陣，便于提取信號(hào)的時(shí)域和頻域特征.廣義S變換（generalized S-transform，GST）還可以通過改進(jìn)窗函數(shù)，提高其時(shí)-頻特征表現(xiàn)能力［7］.

常用的模式識(shí)別方法包括神經(jīng)網(wǎng)絡(luò)（neural networks，NNs）［8］、支持向量機(jī) （support vector machine，SVM）［6］與決策樹（decision tree，DT）［9］等.相比其他方法，決策樹分類效率高、實(shí)現(xiàn)簡單，適用于對(duì)分類效率要求高的應(yīng)用環(huán)境.

本文首先針對(duì)擾動(dòng)信號(hào)特征頻域分布特點(diǎn)開展分析，設(shè)計(jì)具有多重窗口寬度的多分辨率廣義S變換（multiresolution generalized S-transform，MGST）；然后，從原始信號(hào)及MGST時(shí)-頻模矩陣中提取6種特征，構(gòu)建決策樹，并建立最小分類損失原則，確定含交叉樣本的節(jié)點(diǎn)分類閾值；最后，通過對(duì)仿真與實(shí)測(cè)數(shù)據(jù)進(jìn)行分類，對(duì)比不同方法的分類準(zhǔn)確率，驗(yàn)證新方法的有效性與魯棒性.

1 多分辨率廣義S變換基本原理

1.1 廣義S變換

Stockwell等［10］于1996年提出了 S變換方法.設(shè)輸入信號(hào)為h（t），經(jīng)過S變換后為

S變換結(jié)果為一個(gè)二維復(fù)矩陣，稱 S矩陣（S-matrix）.對(duì)矩陣各元素求模后得到S模矩陣，列向量反映某時(shí)刻信號(hào)幅頻特性，行向量描述信號(hào)在特定頻率下的時(shí)域分布.

非平穩(wěn)信號(hào)中，不同頻率成分在發(fā)生畸變時(shí)的時(shí)-頻分布特點(diǎn)不同.其中，信號(hào)的高頻部分變化劇烈，而低頻部分變化相對(duì)平穩(wěn)［11］.國內(nèi)學(xué)者為獲得更好的時(shí)-頻分辨率，引入窗寬調(diào)整因子 λ，使通過調(diào)整λ的值，使窗寬隨頻率成反比變化的速度發(fā)生改變，改進(jìn)后的S變換形式為

由式（3）可以得到廣義S變換的離散表達(dá)式為

其中：f→ n／NT，τ→ jT；G（m，n）由高斯窗函數(shù)w（t，f）經(jīng)快速傅里葉變換獲得.

1.2 不同擾動(dòng)信號(hào)時(shí)-頻特征表現(xiàn)能力要求分析

從擾動(dòng)信號(hào)能量的時(shí)-頻分布特性分析，電壓暫降、暫升、中斷、閃變4類擾動(dòng)信號(hào)能量集中于基頻50 Hz附近；諧波信號(hào)能量分布于基頻和諧波頻率（一般考慮13次以下奇次諧波）附近；暫態(tài)振蕩信號(hào)擾動(dòng)成分能量分布頻域范圍高于諧波頻率；電壓尖峰、電壓切痕的擾動(dòng)成分能量分布在基頻和諧波頻率附近，但諧波頻域部分信號(hào)能量較小，容易受噪聲干擾.從擾動(dòng)分析的需要看，低頻擾動(dòng)分析內(nèi)容包括信號(hào)幅值變化及起止點(diǎn)定位等，需要GST具有更高的時(shí)間分辨率；諧波等中頻擾動(dòng)的分析目的為確定信號(hào)是否含有諧波成分，需要更高的頻域分辨率；高頻特征用于識(shí)別振蕩與含振蕩的復(fù)合擾動(dòng)，要避免其他類型擾動(dòng)信號(hào)的高頻能量和噪聲的影響.因此，需要分別考慮不同頻域的窗寬調(diào)整因子設(shè)定.

1.3 多分辨率廣義S變換

現(xiàn)有廣義S變換在分析復(fù)合擾動(dòng)時(shí)，如果信號(hào)為含有諧波的復(fù)合擾動(dòng)信號(hào)（如諧波含暫降等），GST采用較小λ值進(jìn)行分析，對(duì)復(fù)合擾動(dòng)中的暫降、暫升等分析能力明顯下降，影響分類準(zhǔn)確率.

本文設(shè)計(jì)一種分頻域優(yōu)化的多分辨率廣義S變換方法（multiresolution generalized S-transform，MGST）.將GST模矩陣按照擾動(dòng)信號(hào)擾動(dòng)成分所在頻域分為低頻、中頻、高頻3個(gè)頻域，針對(duì)不同頻域范圍內(nèi)擾動(dòng)類型識(shí)別需要，定義不同的λ值.低頻頻域?yàn)?～100 Hz，主要用于分析擾動(dòng)信號(hào)中包含的暫降、暫升、中斷、閃變、尖峰、切痕6種擾動(dòng)現(xiàn)象；中頻部分為101～700 Hz，主要識(shí)別擾動(dòng)信號(hào)中的諧波成分；700 Hz以上為高頻頻域，主要識(shí)別振蕩成分.通過分別設(shè)定不同頻域的窗寬調(diào)整因子，以滿足復(fù)合擾動(dòng)分析要求.

1.3.1 低頻與中頻部分窗寬調(diào)整因子設(shè)定

由于將單一信號(hào)進(jìn)行分頻域處理，不需要考慮不同頻域時(shí)-頻分辨率難以兼顧的問題，因此新方法低頻部分窗寬可以比GST更窄，取窗寬調(diào)整因子λLF=2；中頻部分比 GST更寬，取窗寬調(diào)整因子.窗寬調(diào)整因子通過統(tǒng)計(jì)實(shí)驗(yàn)確定.

1.3.2 高頻部分窗寬調(diào)整因子確定

為區(qū)分電壓暫降、中斷和含振蕩擾動(dòng)，應(yīng)保留較高的振蕩信號(hào)高頻能量，MGST需要具有更高的時(shí)間分辨率；與噪聲信號(hào)區(qū)分，則需要較高的頻率分辨率，以抑制噪聲能量對(duì)分類結(jié)果的干擾.因此，高頻部分采用單一窗寬調(diào)整因子不能滿足振蕩識(shí)別的需要.

因此，MGST對(duì)高頻部分窗寬調(diào)整因子進(jìn)行自適應(yīng)設(shè)定.通過基頻快速傅里葉譜值A(chǔ)F判斷信號(hào)是否含有基頻擾動(dòng)（暫降或中斷）.為避免噪聲干擾，通過統(tǒng)計(jì)實(shí)驗(yàn)，設(shè)定當(dāng)0.997 pu≤AF≤1.003 pu時(shí)，判斷信號(hào)基頻無擾動(dòng)，取否則，判斷信號(hào)基頻有擾動(dòng)，取

1.3.3 MGST運(yùn)算流程

MGST運(yùn)算流程如圖1所示，虛線框內(nèi)部分為較S變換增加的流程，新方法僅在S變換的基礎(chǔ)上增加了基頻譜值判定環(huán)節(jié)，未增加S變換的復(fù)雜度.圖中，F(xiàn)FT為快速傅里葉變換，IFFT為逆傅里葉變換.

圖1 MGST計(jì)算流程

比較采用MGST、ST、GST、HST 4種方法處理100組擾動(dòng)信號(hào)所用時(shí)間，結(jié)果如圖2所示.HST需要處理更復(fù)雜的窗函數(shù)，GST需要對(duì)諧波頻域是否存在諧波進(jìn)行判斷，因此，計(jì)算時(shí)間略高于ST與MGST.MGST與ST所用計(jì)算時(shí)間接近.與HST與GST相比較，新方法運(yùn)算效率更高.

圖2 采用不同方法計(jì)算100組數(shù)據(jù)時(shí)間比較

1.4 MGST性能分析

圖3比較了 MGST與 ST［7］、GST［8］、HST［14］4種方法在分析相同擾動(dòng)信號(hào)時(shí)，中頻與低頻部分的時(shí)-頻分析能力.圖3（a）為不同變換后的基頻幅值曲線.分析信號(hào)類型為諧波含暫降，信號(hào)含0.1 pu下跌幅度的電壓暫降和3、5、7次諧波，暫降發(fā)生的起始時(shí)刻為第1 168采樣點(diǎn)，結(jié)束時(shí)間為第1 200采樣點(diǎn)，持續(xù)0.5周期.由圖3（a）可知，MGST針對(duì)基頻幅值下跌程度、時(shí)間定位等的分析結(jié)果相對(duì)準(zhǔn)確，性能優(yōu)于ST、GST與HST.圖3（b）比較了MGST與其他方法的諧波分析能力.所分析諧波信號(hào)為h（t）= sin ω0t+0.7sin 3ω0t+0.08sin 5ω0t+0.4sin 7ω0t.新方法分析結(jié)果受旁瓣效應(yīng)影響更小，諧波分析更準(zhǔn)確.為避免振蕩信號(hào)特征受噪聲影響，采用信號(hào)的高頻頻域能量特征識(shí)別振蕩.當(dāng)基頻存在擾動(dòng)時(shí)，MGST在高頻部分窗寬調(diào)整因子取值與GST相同，具有相同的特征表現(xiàn)能力.當(dāng)基頻無擾動(dòng)時(shí)，MGST頻率分辨率高于其他方法.

圖3 低頻與中頻時(shí)-頻分析能力比較

2 擾動(dòng)信號(hào)特征提取

2.1 擾動(dòng)信號(hào)建模

擾動(dòng)信號(hào)的分類研究需要覆蓋所有擾動(dòng)參數(shù)范圍與不同信噪比的大量擾動(dòng)數(shù)據(jù).由于實(shí)測(cè)數(shù)據(jù)一般僅包含一定范圍內(nèi)的擾動(dòng)參數(shù)，不能滿足全面分析的需要，因此，一般采用仿真方法獲得擾動(dòng)信號(hào)，采用MATLAB7.0生成擾動(dòng)信號(hào)，相關(guān)參數(shù)參考文獻(xiàn)［11-13］設(shè)計(jì).其中，設(shè)定電力系統(tǒng)基頻為50 Hz，標(biāo)準(zhǔn)信號(hào)幅值取標(biāo)幺值1 pu，信號(hào)采樣率為3.2 kHz.

2.2 基于MGST的擾動(dòng)信號(hào)分析

從單一擾動(dòng)看，電壓暫降、暫升、中斷、閃變、尖峰、切痕6類擾動(dòng)的特征主要集中于低頻部分，且基頻畸變特性各有不同；諧波主要集中于中頻部分；振蕩主要集中于高頻部分.復(fù)合擾動(dòng)同時(shí)兼具不同擾動(dòng)類型特征，分別體現(xiàn)在不同頻域范圍內(nèi).因此可以從不同頻域范圍內(nèi)提取特征用于擾動(dòng)識(shí)別.

與其他方法相比，分頻域提取特征可以通過MGST的多分辨率特性獲得更具針對(duì)性的特征.

2.3 擾動(dòng)信號(hào)特征提取

綜上分析，確定6種特征構(gòu)建分類決策樹.從原始信號(hào)中提取信號(hào)1／4周期能量跌落幅度D（特征1，F(xiàn)1）與上升幅度R（特征2，F(xiàn)2），共2個(gè)特征.從MGST模矩陣中提取基頻對(duì)應(yīng)幅值標(biāo)準(zhǔn)差σFstd（特征3，F(xiàn)3）、基頻歸一化幅值因子 Af（特征 4，F(xiàn)4）、中頻各個(gè)頻率對(duì)應(yīng)平均幅值最大值A(chǔ)Mmax（特征5，F(xiàn)5）與高頻能量EHF（特征6，F(xiàn)6），共4個(gè)特征.

F1與F2的特征計(jì)算方法分別為

F3、F4、F5、F6的特征計(jì)算方法分別為

其中：n0為基頻；AMax為基頻各采樣點(diǎn)對(duì)應(yīng)幅值中的最大值；AMin為基頻各采樣點(diǎn)對(duì)應(yīng)幅值中的最小值；nH為中頻部分頻率，101 Hz≤nH≤700 Hz；Sth為高頻頻域內(nèi)幅值大于特定閾值的元素，經(jīng)統(tǒng)計(jì)實(shí)驗(yàn)得出，當(dāng)閾值設(shè)為0.02時(shí)，降噪效果較好.

3 基于決策樹的擾動(dòng)信號(hào)分類器設(shè)計(jì)

決策樹通過將復(fù)雜問題轉(zhuǎn)化為若干二分類問題，實(shí)現(xiàn)擾動(dòng)信號(hào)的識(shí)別.分類效率高，實(shí)現(xiàn)簡單［9］.但其分類效果依賴于特征的分類能力，且不同噪聲水平下，最優(yōu)分類閾值變化較大.電力系統(tǒng)中的實(shí)測(cè)信號(hào)信噪比難以確定，因此，應(yīng)該在隨機(jī)噪聲環(huán)境下，確定節(jié)點(diǎn)閾值.

3.1 采用最小分類損失原則確定分類閾值

本文識(shí)別包括電壓暫降C1、電壓暫升C2、電壓中斷C3、閃變C4、暫態(tài)振蕩C5、諧波C6、電壓切痕C7、電壓尖峰C8共8種單一擾動(dòng)，以及諧波含暫降C9、諧波含暫升C10、諧波含閃變C11、諧波含振蕩C12、暫降含振蕩C13共5種復(fù)合擾動(dòng).仿真生成覆蓋不同擾動(dòng)參數(shù)且信噪比為30～50 dB隨機(jī)仿真信號(hào)，每類500組，共6 500組開展統(tǒng)計(jì)分析.并設(shè)計(jì)最小分類損失原則，確定分類閾值.閾值的確定過程如圖4所示.

圖4 基于最小分類損失原則的閾值確定過程

圖4（a）為暫降和暫降含振蕩兩類擾動(dòng)的100組樣本特征值分布.圖4（b）為樣本F6特征值交叉部分.為了確定最優(yōu)閾值，定義誤識(shí)別率為

其中：p為誤識(shí)別待分類樣本數(shù)，M為樣本總數(shù).

設(shè)交叉范圍內(nèi)樣本特征值為分類閾值，在其中尋找Er最小值對(duì)應(yīng)的閾值，即確定具有最小分類損失的閾值.閾值確定過程：1）確定交叉樣本區(qū)域，以交叉樣本特征的最大值與最小值作為閾值上限與下限，如圖4（b）所示；2）交叉區(qū)域樣本特征值分別作為分類閾值，并計(jì)算取不同閾值時(shí)的Er；3）選擇對(duì)應(yīng)Er最小值Ermin的特征值作為分類閾值，如果樣本多于1個(gè)，則取最優(yōu)閾值范圍中間值為分類閾值.

由圖4可知，當(dāng)訓(xùn)練樣本較多時(shí)，閾值更加精確.因此，決策樹各節(jié)點(diǎn)分類閾值均通過每類擾動(dòng)樣本500組分析確定.

3.2 基于決策樹的擾動(dòng)信號(hào)分類器設(shè)計(jì)

決策樹分類器結(jié)構(gòu)如圖5所示.其中，無交叉樣本節(jié)點(diǎn)閾值由無交叉特征值范圍的中間值確定；含交叉樣本的節(jié)點(diǎn)閾值由最小分類損失原則確定.

圖5 決策樹結(jié)構(gòu)圖

4 仿真與實(shí)測(cè)數(shù)據(jù)驗(yàn)證

4.1 仿真實(shí)驗(yàn)

仿真生成信噪比分別為30、40、50 dB以及30～50 dB隨機(jī)值的仿真信號(hào)，每類各500組，驗(yàn)證新方法有效性.并構(gòu)建基于ST［6］、GST［7］、HST［14］的決策樹，開展比較試驗(yàn).表1為隨機(jī)噪聲環(huán)境下，應(yīng)用4種分類方法分別進(jìn)行分類時(shí)的準(zhǔn)確率比較表.由表1可見，新方法在識(shí)別不同類別擾動(dòng)信號(hào)時(shí)，分類準(zhǔn)確率均為最高.在識(shí)別復(fù)合擾動(dòng)C9～C13時(shí)，優(yōu)勢(shì)尤其明顯.表2分別為30、40、50 dB噪聲環(huán)境下4種分類方法的分類準(zhǔn)確率.由表2可見，新方法在不同噪聲水平下的分類總準(zhǔn)確率均高于其他方法，且保持在97.66%以上，具有良好的抗噪性和魯棒性.

4.2 實(shí)測(cè)信號(hào)實(shí)驗(yàn)

采用葡萄牙某電網(wǎng)2006年11月間實(shí)測(cè)單相電能質(zhì)量信號(hào)952組開展分析［15］.根據(jù)文獻(xiàn)［15］提出的擾動(dòng)識(shí)別方法，該組實(shí)測(cè)信號(hào)中25組擾動(dòng)信號(hào)被識(shí)別為暫降，5組被識(shí)別為中斷，910組被識(shí)別為暫態(tài)振蕩，8組被識(shí)別為諧波，4組信號(hào)擾動(dòng)類型未能確定.采用新方法識(shí)別以上實(shí)測(cè)信號(hào)，驗(yàn)證新方法有效性.實(shí)驗(yàn)過程中取采樣率為10 kHz，并對(duì)信號(hào)電壓幅值進(jìn)行歸一化預(yù)處理.

表1 分類準(zhǔn)確率比較

表2 不同噪聲環(huán)境下分類準(zhǔn)確率比較

實(shí)驗(yàn)結(jié)果顯示，新方法對(duì)中斷與暫態(tài)振蕩信號(hào)的識(shí)別結(jié)果與原系統(tǒng)相同.由于對(duì)暫降幅值變化程度定義不同，原系統(tǒng)識(shí)別的25組暫降信號(hào)中，19組信號(hào)幅值下跌幅度不足0.1 pu，不滿足IEEE 1159—1995標(biāo)準(zhǔn)定義，因此不予考慮；在下跌幅值超過0.1 pu的6組樣本中，3組識(shí)別為電壓暫降；3組含有諧波成分，識(shí)別為諧波含暫降.新方法將原方法的8組諧波信號(hào)中的5組識(shí)別為諧波，3組識(shí)別為諧波含閃變.原系統(tǒng)無法識(shí)別的4種擾動(dòng)，由新方法分析后得到的特征值及識(shí)別結(jié)果見表3.綜合以上分析可知，文獻(xiàn)［15］采用的識(shí)別系統(tǒng)未考慮復(fù)合擾動(dòng)情況，而新方法提高了復(fù)合擾動(dòng)的識(shí)別能力，成功的識(shí)別出原系統(tǒng)未發(fā)現(xiàn)的復(fù)合擾動(dòng)，具有更好的復(fù)合擾動(dòng)識(shí)別能力，能夠滿足實(shí)際電能質(zhì)量信號(hào)分析的要求.

表3 未知類型實(shí)測(cè)信號(hào)特征值與識(shí)別結(jié)果

5 結(jié) 論

1）針對(duì)不同擾動(dòng)信號(hào)的特征表現(xiàn)要求，自適應(yīng)定義不同頻域范圍內(nèi)的窗函數(shù)寬度調(diào)整因子，提高了廣義S變換的復(fù)合擾動(dòng)識(shí)別能力.

2）從不同頻域提取擾動(dòng)特征，降低了特征運(yùn)算量，避免了不同頻域信號(hào)成分對(duì)特征值的干擾.同時(shí)，從原始信號(hào)中提取的特征，能夠有效提高參數(shù)范圍相鄰的暫降與中斷兩類擾動(dòng)的識(shí)別準(zhǔn)確率.

3）設(shè)計(jì)了隨機(jī)噪聲環(huán)境下的最小分類損失原則，用于確定含交叉樣本的節(jié)點(diǎn)的分類閾值.

4）通過仿真與實(shí)測(cè)信號(hào)實(shí)驗(yàn)，驗(yàn)證了新方法的有效性與實(shí)用性.將新方法與現(xiàn)有研究成果中的不完全S變換、并行S變換等相結(jié)合，能夠提高M(jìn)GST運(yùn)算速度.此外，現(xiàn)有研究未考慮對(duì)高次諧波的識(shí)別，如何區(qū)分高次諧波和暫態(tài)振蕩是未來工作的難點(diǎn).

［1］LIEBERMAND G，ROMERO-TRONCOSO R J，OSORNIORIOS R A，et al.Techniques and methodologies for power quality analysis and disturbances classification in power systems：a review ［J］.IET Generation，Transmission＆Distribution，2011，5（4）：519-529.

［2］張伯明，孫宏斌，吳文傳，等.智能電網(wǎng)控制中心技術(shù)的未來發(fā)展［J］.電力系統(tǒng)自動(dòng)化，2009，33（17）：21-28.

［3］FARZANEHRAFAT A，WATSON N R.Power quality state estimatorfor smart distribution grids ［J］. IEEE Transactions on Power System，2013，28（3）：2183-2191.

［4］BISWALM，DASH P K.Detection and characterization of multiple power quality disturbances with a fast S-transform and decision tree based classifier［J］.Digital Signal Processing，2013，23：1071-1083.

［5］李天云，趙研，李楠，等.基于HHT的電能質(zhì)量檢測(cè)新方法［J］.中國電機(jī)工程學(xué)報(bào)，2005，25（17）：52-56.

［6］李濤，夏浪，張宇，等.基于提升復(fù)小波的暫態(tài)電能質(zhì)量擾動(dòng)的檢測(cè)與定位［J］.中國電機(jī)工程學(xué)報(bào)，2011，25（9）：66-72.

［7］黃南天，徐殿國，劉曉勝.基于S變換與SVM的電能質(zhì)量復(fù)合擾動(dòng)識(shí)別［J］.電工技術(shù)學(xué)報(bào)，2011，10（10）：23-30.

［8］徐方維，楊洪耕，葉茂清，等.基于改進(jìn)S變換的電能質(zhì)量擾動(dòng)分類［J］.中國電機(jī)工程學(xué)報(bào)，2012，32（4）：77-84.

［9］林圣，何正友，臧天磊，等.基于粗神經(jīng)網(wǎng)絡(luò)的輸電線路故障分類方法［J］.中國電機(jī)工程學(xué)報(bào)，2010，30（28）：72-79.

［10］STOCKWELL R G，MANSINHA L，LOWE R P. Localization of the complex spectrum：the S transform［J］. IEEE Trans on Signal Process，1996，44（4）：998-1001.

［11］王麗霞，何正友，趙靜.一種基于線性時(shí)頻分布和二進(jìn)制閾值特征矩陣的電能質(zhì)量分類方法［J］.電工技術(shù)學(xué)報(bào)，2011，26（4）：185-191.

［12］占勇，程浩忠，丁屹峰，等.基于S變換的電能質(zhì)量擾動(dòng)支持向量機(jī)分類識(shí)別［J］.中國電機(jī)工程學(xué)報(bào)，2005，25（4）：51-56.

［13］HOOSHMAND R， ENSHAEE A.Detection and classification ofsingle and combined power quality disturbances using fuzzy systems oriented by particle swarm optimization algorithm ［J］.Electric PowerSystems Research，2010，80：1552-1561.

［14］ROBERT P C，MANSINHA L.The S-transform with windows of arbitrary and varying shape［J］.Geophysics，2003，68（1）：381-385.

［15］RADIL T，RAMOS P M，JANEIRO F M，et al.PQ monitoring system for real-time detection and classification of disturbances in a single-phase power system［J］.IEEE Transactions on Instrumentation and Measurement，2008，57（8）：1725-1733.

（編輯 魏希柱）

Classification of power quality disturbances utilizing multiresolution generalized S-transform

HUANG Nantian1，ZHANG Weihui1，XU Dianguo2，CAI Guowei1，LIU Chuang1，ZHANG Shuxin1

（1.School of Electrical Engineering，Northeast Dianli University，132012 Jilin，Jilin，China；2.School of Electrical Engineering and Automation，Harbin Institute of Technology，150001 Harbin，China）

In order to improve the ability of complex power quality disturbances recognition，a new type of complex disturbances recognition approach based on Multiresolution Generalized S-transform（MGST）is proposed.Firstly，the spectrum of original signals is segmented into 3 frequency areas including low frequency area，medium frequency area and high frequency area.The width factor of window function in S-transform is defined respectively in different frequency areas.MGST has different time-frequency resolution in each frequency area in order to satisfy the recognition requirements of different disturbances in each frequency area.Otherwise，the width factor of window function in the high frequency area is adaptively adjusted according to the value of Fourier spectrum of the fundamental frequency.On this basis，the decision tree based on 6 features is constructed to recognize disturbance signals.Finally，the minimum classification faults rule is designed to get the optimum threshold of each node.The simulation and real signals experiments show that 13 types of disturbances including 5 types of complex disturbances are recognized accurately by the new approach.The new approach has better classification accuracy and noise immunity than other methods such as S-transform，generalized S-transform and Hyperbolic S-transform.

power quality；power quality disturbance；S-transform；multiresolution；decision tree

TM714.3

A

0367-6234（2015）09-0051-06

10.11918／j.issn.0367-6234.2015.09.010

2014-07-05.

國家自然科學(xué)基金 （51307020）；吉林省科技發(fā)展計(jì)劃（20150520114JH）；吉林市科技發(fā)展計(jì)劃（201464052）.

黃南天（1980—），男，博士，副教授；徐殿國（1960—），男，教授，博士生導(dǎo)師；蔡國偉（1968—），男，教授，博士生導(dǎo)師.

黃南天，huangnantian＠126.com.