徐健,周著斌
(西安工程大學(xué) 電信學(xué)院,陜西 西安710048)
混合電能質(zhì)量檢測方法的研究
徐健,周著斌
(西安工程大學(xué) 電信學(xué)院,陜西 西安710048)
結(jié)合電能穩(wěn)態(tài)和暫態(tài)的分析方法,提出了對兩者的混合擾動的檢測。短時傅里葉變換(STFT)提出了加短時窗的思想,小波變換改進(jìn)了STFT窗口寬度固定的缺陷,S變換改進(jìn)了小波變換窗口面積固定,不能提取任意頻次成分的缺點。并對比分析了這三種時頻分析法對電壓暫降與諧波混合的信號的檢測結(jié)果,隨著學(xué)術(shù)思想的進(jìn)步,檢測結(jié)果更加精確。
混合擾動;短時傅里葉變換;小波變換;S變換
電能質(zhì)量可以簡單的分為暫態(tài)和穩(wěn)態(tài)擾動,現(xiàn)今對這兩種擾動的研究也很多,尤其是對穩(wěn)態(tài)的分析[1]。穩(wěn)態(tài)電能質(zhì)量,如電壓波動與閃變、諧波、三相不平衡[2],其主 要檢測的是波動幅值、不平衡因子、出現(xiàn)頻率等特性;暫態(tài)電能質(zhì)量,如脈沖震蕩、電壓暫升、電壓暫降、電壓中斷[3],其主要檢測的是峰值、頻譜、起止時刻、持續(xù)時間等特性。然而在現(xiàn)實生活中,出現(xiàn)的電能質(zhì)量往往不是單一形式的,而是兩種甚至多種電能質(zhì)量的混合[4],并且還混有噪聲等影響因素,所以對混合電能質(zhì)量的研究很有必要。電壓暫降和諧波的混合擾動是典型的混合臺電能質(zhì)量。
電能質(zhì)量檢測的方法很多,如傅立葉變換,短時傅立葉變換,dq變換,小波變換,S變換等[5]。
短時傅立葉變換,小波變換,S變換雖是3種不同名稱的檢測方法,其基本思想?yún)s是相同的,3種方法是依次改進(jìn)的關(guān)系[6]。3種方法的理論思路相同,隨著時代的進(jìn)步,從短時傅立葉變換到小波變換再到S變換其學(xué)術(shù)思想也在進(jìn)一步改進(jìn),在電能質(zhì)量中的檢測范圍更廣,檢測值更加精確。用這3種方法分別對諧波與電壓暫降組成的混合態(tài)電能質(zhì)量堿性檢測,可以清晰地看出3種方法的聯(lián)系與區(qū)別。
1.1短時傅里葉變換
短時傅里葉變換(STFT)是傅里葉變換的一種改進(jìn),其在傅里葉變換的基礎(chǔ)上提出了加短時窗的概念,從而能夠確定時變信號局部信號的頻率[7]。
離散傅里葉變換數(shù)學(xué)表達(dá)式為:
短時傅里葉變換表達(dá)式為:
由表達(dá)式可知短時傅里葉變換的改進(jìn)之處主要體現(xiàn)在其的加窗性,且窗函數(shù)可以選擇漢寧窗、漢明窗、高斯窗。移動窗函數(shù),使原始信號h(t)與g(t)在不同的有限時間寬度內(nèi)是平穩(wěn)信號,從而得到原始信號的光譜。
1.2小波變換
小波變換是在短時傅里葉變換基礎(chǔ)上的改進(jìn),短時傅里葉變換的缺點是一旦窗函數(shù)確定,其窗口寬度將確定,這對非平穩(wěn)信號的研究有很大的局限性,小波變換在此進(jìn)行了提升[8]。
連續(xù)小波變換表達(dá)式:
其中,a:伸縮因子,b:平移因子
由表達(dá)式可知,小波變換窗口可以平移,即位置不固定,具有可變性,且其形狀也可以變,較短時傅里葉變換時變性更好。
1.3S變換
S變換是小波變換學(xué)術(shù)思想的進(jìn)一步提高,它克服了小波變換窗口大小固定,受噪聲影響較大等缺陷。
S變換的表達(dá)式為:
由表達(dá)式可知,S變換的窗寬由頻率決定,窗寬可變,且可以單獨(dú)提取任意頻次的信號,較小波變換實時性更好[9]。
2.1混合態(tài)電能質(zhì)量的模型建立
電壓暫降數(shù)學(xué)模型為[10]:
其中,α的范圍為0.1—0.9;t1<t2;t≥1時,u(t)=1,其余值為0
諧波數(shù)學(xué)模型為[11]:
其中0.05≤αi≤0.15
故可建立混合態(tài)模型為:
取α=0.5,t1=0.3,t2=0.7,w=100∏,α3=α5=α7=0.1
為了更接近實際效果,對該模型加入30 db的噪聲,然后分別進(jìn)行仿真處理。
2.2三種方法的matlab仿真
圖1 短時傅里葉變換三維圖
圖2 小波變換分解圖
圖3 S變換三維圖
圖4 短時傅里葉及小波變換幅頻曲線
圖1、圖2、圖3是3種方法對混合信號的一種概況描述,從這3幅仿真圖也可以看出3種方法的對波形進(jìn)行檢測的數(shù)學(xué)思維[11]。圖1中能基本確定所檢測信號的類型,但并不清晰,因為短時傅里葉是對于原始信號先乘以某確定的某種窗函數(shù),該窗函數(shù)僅在一段時間不為零其仿但一旦確定窗函數(shù)窗寬等因素將無法進(jìn)行進(jìn)一步調(diào)整;圖2是采用的小波變換多分辨率思想,通過設(shè)置采樣頻率,提取出我們需要檢測的頻率下的波形,較短時傅里葉變換更靈活些,對細(xì)節(jié)的提取較精確;圖3中S變換所做的分析圖更加清晰,由于其變換之處的優(yōu)勢在于高斯窗口的高度和寬度隨頻率而變化,這樣就克服了短時傅葉變換窗口高度和寬度固定的缺陷,也克服了小波變換窗口大小固定的缺陷。
3種變換對被檢測信號的幅頻特性分析如圖4、5所示。
圖5 S變換對幅值特性的分析
3個仿真圖都能分析出被檢測信號的幅值信息,由圖4左邊可知短時傅里葉變換得到的幅值曲線明顯不如小波變換得到的幅值精確,小波變換幅值提取過程是基于多分辨率原理[12],逐步提取的。圖4右邊可以看出a1中可以測出頻率f= 50,150,250,350的幅值,a2可以測出f=50,150,250時的幅值,a4時只能測出基頻幅值,短時傅里葉與小波變換除了在f= 50,150,250,350等處,其他頻率所對應(yīng)的幅值基本為零,即短時傅里葉變換與小波變換沒有反映出各頻率之間的變化情況。而圖5中S變換則反映出了幅值隨著頻率的增大而減小這一特性[13]。
3種變換對被檢測信號的突變特性分析如圖6~8。
圖6 短時傅里葉變換光譜圖(窗寬M=128)
圖7 小波熵對信號的檢測
圖8 S變換的幅值包絡(luò)圖
圖6中,亮色代表原始信號中所含的頻率成分,亮色程度的面積代表各個頻率成分所占的比例,可以看出短時傅里葉變換基本顯示出了混合信號的特性,但不能得到其精確值;圖7中采用小波變換與信息熵相結(jié)合的思想,由圖可知在0.3、0.7附近,熵系數(shù)明顯增大,即反映出此時被檢測信號成分發(fā)生較大變化,且由于噪聲和諧波的存在,熵系數(shù)并不穩(wěn)定于某一穩(wěn)定值。較短時傅里葉變換,小波變換的熵系數(shù)能夠反映出被檢測信號的整個非穩(wěn)態(tài)的大致變化;圖8是S變換下的檢測情況,被檢測信號的頻率成分,幅值變化一目了然,說明了S變換在非穩(wěn)態(tài)檢測中的優(yōu)勢[14]。
為了更加精確地說明3種變化的區(qū)別,分別提取信號被檢測以后的突變起止時刻,幅值,得表如表1。
表1 3種方法檢測特征量的提取值
表1數(shù)據(jù)的數(shù)據(jù)更加清楚地說明了三者對于非平穩(wěn)信號的檢測精度,但由于噪聲的存在,3種方法均不能達(dá)到最優(yōu),在實際檢測中,需要消噪處理[15]。
文章以短時傅里葉為基礎(chǔ),一步步提出了小波變換、S變換,三種算法的基本學(xué)術(shù)思想相同,但理論知識的應(yīng)用在提高。
電壓暫降與諧波的混合態(tài)擾動,更加符合實際電能質(zhì)量狀態(tài)。通過3種方法對混合態(tài)信號的檢測仿真圖對比,可得結(jié)果如下:短時傅里葉變換能基本確定混合信號的突變的起止時刻,幅頻變化信息,但誤差較大;小波變換對非平穩(wěn)信號的特征的提取明顯優(yōu)于STFT,但并不能提取任意頻次的信號,且受噪聲影響較大;S變換不管是對幅頻值的檢測,還是時間的檢測,精度都明顯高于前兩者。
[1]YUAN Shuai,Huang Runchang,TONG Weiming,et al.Voltage flicker measurement based on DSP platform and its error correction[J].The power system protection and control,2008,4(21):47-51.
[2]程浩忠,電能質(zhì)量介紹[M].北京:中國電力出版社,2008.
[3]程浩志,呂干云,周荔丹.電能質(zhì)量監(jiān)測與分析[M].北京科學(xué)出版社,2012.
[4]Youssef A M,Abdel-Galil T K,El-Saadany E F,et al.Disturbance classification utilizing dynamic time warping classifier [J].IEEE Transon Power Delivery,2004,19(1):272-278.
[5]桑松,柴玉華,孫影.基于小波變換和快速傅里葉變換的諧波檢測[J].電測與儀表,2012,49(7):29-32.
[6]Dash P K,Panigrahi B K,Sahoo D K,et al.Power quality disturbance data compression,detection,and classification using integrated spline wavelet and S-transform[J].IEEE Transactions on Power Delivery,2003,18(2):595-600.
[7]趙鳳展,楊仁剛.基于短時傅里葉變換的電壓暫降擾動檢測[J].中國電機(jī)工程學(xué)報,2007,27(10):28-34.
[8]陳祥訓(xùn).利用小波變換測量與識別輕微幅度變化與相位變化的方法[J].中國電機(jī)工程學(xué)報,2003,23(12):12-18.
[9]唐求,滕召勝,高云鵬,等.基于S變換的平方檢測法測量電壓閃變,中國電機(jī)工程學(xué)報,2012,32(7):60-67.
[10]周林,吳紅春,夢婧,等.電壓暫降分析方法研究[J].高電壓技術(shù),2008,34(5):1001-1008.
[11]張伏生,耿中行,葛耀中.電力系統(tǒng)諧波分析的高精度FFT算法[J].中國電機(jī)工程學(xué)報,1999,19(3):63-66.
[12]任 震,黃群古,黃雯瑩.多層多分辨分析算法及在電力系統(tǒng)中的應(yīng)用[J].電網(wǎng)技術(shù),2012(6):5-7.
[13]劉守亮,肖先勇,楊洪耕.基于S變換模時頻矩陣相似度的短時電能質(zhì)量擾動分類[J].電網(wǎng)技術(shù).2006(5):22-25.
[14]Stockwell R G,Mansinha L,Lowe R P.Localization of the Complex Spectrum:The S Transform[M].IEEE Transactions on Signal Processing,1996,44(4):998-1001.
[15]Pinnegar C R,Mansiha L.Time-local spectral analysis for non-stationary time series:the S-transform for noise signal[J].Fluctuation and Noise Letters,2003,3(3):357-364.
Study on the method of hybrid power quality detection
XU Jian,ZHOU Zhu-bin
(College of Electronic Information of Xi'an Polytechnic University,Xi'an 710048,China)
Combined with the analysis method of electric energy and transient state,the detection of mixed disturbance is proposed.The short-time Fourier transform(STFT)is proposed to add short time window;Wavelet transform improves the STFT's window width and fixed the defect;S transform is an improvement of wavelet transform which improves the window area of the wavelet transform and can extract any frequency components.And the detection results of the three time-frequency methods are compared and analyzed.Along with the progress of the academic thought,the test result is more accurate.
mixed disturbance;short-time Fourier transform;wavelet transform;S transform
TN801
A
1674-6236(2016)22-0053-04
2016-03-12稿件編號:201603143
陜西省科技廳工業(yè)科技攻關(guān)項目(2015GY075)
徐 ?。?963—),男,陜西西安人,副教授。研究方向:物聯(lián)網(wǎng),信號處理,電能質(zhì)量檢測。