陶佳蘭，喻 敏，陳貴詞，王 斌

基于SWT的電力系統(tǒng)基波檢測(cè)

陶佳蘭1,2，喻 敏1,2，陳貴詞1,2，王 斌3

(1.冶金工業(yè)過(guò)程系統(tǒng)科學(xué)湖北省重點(diǎn)實(shí)驗(yàn)室(武漢科技大學(xué))，湖北 武漢 430065；2.武漢科技大學(xué)理學(xué)院，湖北 武漢 430065；3.武漢科技大學(xué)信息科學(xué)與工程學(xué)院，湖北 武漢 430081)

在噪聲混入含有基波的信號(hào)時(shí)，傳統(tǒng)的時(shí)頻分析方法在基波提取過(guò)程中易出現(xiàn)模態(tài)混疊。為了準(zhǔn)確檢測(cè)出基波分量，利用時(shí)頻分析精度較高的同步擠壓小波變換(Synchrosqueezing Wavelet Transform, SWT)實(shí)現(xiàn)基波檢測(cè)。首先，采用SWT將含有基波的信號(hào)分解為一組內(nèi)蘊(yùn)模態(tài)類(lèi)函數(shù)(Intrinsic Mode Type functions, IMTs)，第一個(gè)分量IMT1即代表基波。然后，該分量經(jīng)Hilbert變換實(shí)現(xiàn)基波頻率和幅值的測(cè)量。在諧波幅值瞬變、噪聲混入、基波頻率波動(dòng)、間諧波頻率靠近基波和諧波的情境下進(jìn)行算法驗(yàn)證。實(shí)驗(yàn)結(jié)果表明，SWT能夠準(zhǔn)確提取基波，頻率精度最高可達(dá)10-8量級(jí)，具有較強(qiáng)的抗噪性，且SWT的基波提取能力強(qiáng)于諧波和間諧波。

基波檢測(cè)；模態(tài)混疊；同步擠壓小波變換；希爾伯特變換；抗噪性

0 引言

在電力系統(tǒng)故障診斷、儀器儀表輸出信號(hào)的處理、電能質(zhì)量分析等領(lǐng)域[1]，需要從含有噪聲的信號(hào)中提取基波信息，且各種電力設(shè)備在電力系統(tǒng)非線性負(fù)荷和電力系統(tǒng)電壓作用下，易發(fā)生擾動(dòng)[2]，并伴隨大量影響電能質(zhì)量的諧波和間諧波，而部分諧波和間諧波檢測(cè)的相關(guān)研究，也需依賴(lài)于所提取到的高精度基波頻率[3]。所以，準(zhǔn)確地從電力系統(tǒng)信號(hào)中檢測(cè)出基波參數(shù)對(duì)電力系統(tǒng)的穩(wěn)定運(yùn)行具有重要的工程實(shí)用價(jià)值，同時(shí)也是部分諧波和間諧波檢測(cè)研究工作的基礎(chǔ)[4]。

目前，時(shí)頻分析方法在基波檢測(cè)中較為常見(jiàn)，但不同的時(shí)頻分析方法因電力系統(tǒng)環(huán)境的不同，基波檢測(cè)效果各有優(yōu)劣。快速傅里葉變換(Fast Fourier Transform, FFT)[5]適用于平穩(wěn)周期電力信號(hào)的基波提取，且要求滿足采樣同步，否則容易出現(xiàn)頻譜泄露和柵欄效應(yīng)[6]。而短時(shí)傅里葉變換(Short Time Fourier Transform, STFT)分析窗的長(zhǎng)度固定，對(duì)于含有基波的暫態(tài)擾動(dòng)，不易在時(shí)域和頻域同時(shí)達(dá)到較高的精度檢測(cè)要求，且計(jì)算量稍大[7]。小波變換(Wavelet Transform, WT)[8]的基波檢測(cè)效果則取決于最優(yōu)小波基本函數(shù)的選取[9]。相比傳統(tǒng)窗函數(shù)固定的STFT與基函數(shù)固定的WT，希爾伯特-黃變換(Hilbert-Huang Transform, HHT)[10]包含經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)和Hilbert變換兩部分，因EMD中的基函數(shù)具有自適應(yīng)性，所以更適用于非平穩(wěn)信號(hào)，經(jīng)EMD可得到代表基波和諧波等不同頻率分量的固有模態(tài)函數(shù)(Intrinsic Mode Functions, IMFs)，但EMD存在的模態(tài)混疊問(wèn)題[11]會(huì)影響電力系統(tǒng)信號(hào)分解的精度，繼而產(chǎn)生EMD的系列改進(jìn)方法：集合經(jīng)驗(yàn)?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition, EEMD)[12]和自適應(yīng)噪聲的完全集合經(jīng)驗(yàn)?zāi)B(tài)分解(Complete Ensemble Empirical Mode Decomposition with Adaptive Noise, CEEMDAN)[13]等，但EEMD可能使原信號(hào)被噪聲污染，計(jì)算量大，CEEMDAN計(jì)算量雖有減小，卻和EEMD一樣可能存在虛假分量[14]。而變分模態(tài)分解(Variational Mode Decomposition, VMD)[15]則是將信號(hào)分解問(wèn)題轉(zhuǎn)化為變分問(wèn)題，通過(guò)迭代搜索變分模型中的最優(yōu)解，確定出各模態(tài)分量的頻率中心和最優(yōu)帶寬，但因預(yù)先設(shè)定模態(tài)數(shù)目的差異，存在模態(tài)過(guò)分解和丟失的情況[16]。同步擠壓小波變換[17]是由Daubechies等人提出的一種以小波變換為基礎(chǔ)的時(shí)頻重排算法，該方法在時(shí)間-尺度上對(duì)信號(hào)能量進(jìn)行重新分配，經(jīng)同步擠壓銳化頻率曲線，可重構(gòu)出各頻率曲線所對(duì)應(yīng)分量的波形，具有抗噪性，且能有效改善模態(tài)混疊問(wèn)題[18]。

SWT方法已經(jīng)應(yīng)用在信號(hào)消噪、機(jī)械故障診斷等諸多領(lǐng)域。文獻(xiàn)[19]將SWT應(yīng)用于解決風(fēng)電匯集地區(qū)存在的非同步振蕩問(wèn)題，對(duì)比FFT和HHT，SWT更為準(zhǔn)確地檢測(cè)出了頻率相近的間諧波成分。文獻(xiàn)[20]將SWT用于旋轉(zhuǎn)機(jī)械振動(dòng)信號(hào)特征提取，有效避免了模態(tài)混疊，能夠精確描述諧波信號(hào)的頻率構(gòu)成。文獻(xiàn)[21]首先運(yùn)用EMD方法分解得到含有模態(tài)混疊的IMF3分量，再運(yùn)用SWT對(duì)IMF3分量繼續(xù)進(jìn)行分解，得到了與直接使用SWT分解原始信號(hào)相一致的分解結(jié)果，成功提取出諧波分量。

針對(duì)模態(tài)混疊使基波提取難度加大的問(wèn)題，本文在電力系統(tǒng)電壓存在諧波幅值瞬變、噪聲混入、基波頻率波動(dòng)、間諧波頻率靠近基波和諧波的情形下，結(jié)合SWT和Hilbert變換實(shí)現(xiàn)基波頻率和幅值的測(cè)量。實(shí)驗(yàn)對(duì)比EMD、CEEMDAN、VMD與SWT的基波檢測(cè)效果，分析驗(yàn)證SWT的模態(tài)分解能力、抗噪性能以及在基波檢測(cè)中的準(zhǔn)確性。

1 ?基于同步擠壓小波變換的基波提取原理

1.1 同步擠壓小波變換原理

同步擠壓小波變換通過(guò)對(duì)所劃分的頻率區(qū)間同步擠壓來(lái)銳化頻率曲線，依據(jù)信號(hào)在時(shí)間-尺度中的能量大小逐步分解出代表基波和諧波等不同頻率分量的內(nèi)蘊(yùn)模態(tài)類(lèi)函數(shù)(Intrinsic Mode Type functions, IMTs)。經(jīng)同步擠壓的各IMT之間不存在交叉項(xiàng)，在模態(tài)混疊時(shí)也能準(zhǔn)確提取基波等不同頻率分量。

SWT的基波提取主要包括以下5個(gè)步驟。

1) 對(duì)混有噪聲的多分量電力信號(hào)進(jìn)行連續(xù)小波變換(Continuous Wavelet Transform, CWT)。

對(duì)式(2)進(jìn)行傅里葉變換，得到連續(xù)小波變換系數(shù)在頻域上的計(jì)算公式為

bump小波基函數(shù)在頻域上可表示為

3) 劃分信號(hào)頻率區(qū)間。

5) 基波分量重構(gòu)。

1.2 Hilbert變換原理

經(jīng)1.1節(jié)同步擠壓小波變換的5個(gè)步驟，即可完成基波分量的提取，再結(jié)合1.2節(jié)的Hilbert變換，則能計(jì)算出基波分量的特征頻率和特征幅值。

基于SWT的基波檢測(cè)流程如圖1所示。

圖1 基于SWT的基波檢測(cè)流程圖

2 ?信號(hào)仿真實(shí)驗(yàn)與分析

為了驗(yàn)證SWT方法用于基波檢測(cè)的性能，本文在諧波幅值瞬變、噪聲混入、基波頻率波動(dòng)、間諧波頻率靠近基波和諧波的情境下分別對(duì)暫態(tài)電壓信號(hào)[24-25]開(kāi)展仿真實(shí)驗(yàn)與分析。

2.1 諧波幅值瞬變的仿真信號(hào)

圖2 SWT時(shí)頻分析

將SWT的分解結(jié)果與EMD、CEEMDAN、VMD的分解結(jié)果對(duì)比，其前4項(xiàng)模態(tài)分量波形如圖3所示。4種方法檢測(cè)出在后6個(gè)周期(圖3中箭頭所截取的周期數(shù))的基波分量與原始基波的波形吻合效果，如圖4所示。

由圖3可知，EMD、CEEMDAN、VMD、SWT分解得到的IMF2、IMF3、IMF1、IMT1分量分別對(duì)應(yīng)基波；CEEMDAN的基波提取效果較EMD稍差，波形相比之下略不整齊，兩者都出現(xiàn)了模態(tài)混疊現(xiàn)象，導(dǎo)致不能從信號(hào)()中完全提取出基波；而VMD雖能較為完整地提取出基波，但相比SWT方法的分解結(jié)果，VMD存在邊界效應(yīng)，兩端峰值較小，且波形不整齊。由此可知，SWT具有更好的模態(tài)分解能力，且因SWT檢測(cè)出后3個(gè)IMT模態(tài)分量的頻率依次為89.999 999 Hz、149.718 474 Hz和210.000 000 Hz，表明其也適用于諧波和間諧波檢測(cè)。

圖3 模態(tài)分解結(jié)果

圖4 基波波形吻合效果1

由圖4可知，因3次諧波的幅值在0.6 s和0.9 s時(shí)刻發(fā)生階躍變化， CEEMDAN和EMD的分解結(jié)果受到一定程度影響，頻率和幅值不穩(wěn)定，波形不整齊，它們?cè)诜逯堤幍牟ㄐ闻c原始基波波形的吻合效果較差；而SWT和VMD在峰值處的整體吻合度較好，但對(duì)比之下，SWT在波形右端點(diǎn)處與原始基波波形的吻合效果更好，VMD在最后2個(gè)峰值處的幅值明顯小于基波幅值，與基波波形的吻合效果略差，由此可知，4種方法的基波波形吻合度排序?yàn)镾WT＞VMD＞EMD＞CEEMDAN。

為分析4種方法測(cè)量出的頻率和幅值隨時(shí)間變化的數(shù)值波動(dòng)穩(wěn)定性，列出了4種方法經(jīng)Hilbert變換得到的基波頻率、幅值及其標(biāo)準(zhǔn)差的具體數(shù)值，如表1所示，其與原始基波頻率標(biāo)準(zhǔn)值(50 Hz)和幅值標(biāo)準(zhǔn)值(120 V)之間所存在的相對(duì)誤差如表2所示。

由表1和表2可以看出，SWT結(jié)合Hilbert變換計(jì)算得到的基波頻率和幅值隨時(shí)間波動(dòng)的標(biāo)準(zhǔn)差最小，分別為0.0464和0.2755；對(duì)比4種方法檢測(cè)出的基波頻率和幅值所存在的相對(duì)誤差，精度由高到低排序?yàn)镾WT＞VMD＞EMD＞CEEMDAN；其中，SWT的基波頻率檢測(cè)精度達(dá)到了10-8量級(jí)，相比EMD、CEEMDAN和VMD，SWT的檢測(cè)精度更高，幅值檢測(cè)精度也達(dá)到了10-4量級(jí)，同時(shí)，對(duì)比數(shù)字萬(wàn)用表FLUKE15B的技術(shù)指標(biāo)對(duì)基波測(cè)量參數(shù)的最大允許誤差(±1.0%)[28]，SWT滿足該精度要求，表明該方法能夠較為準(zhǔn)確地實(shí)現(xiàn)基波檢測(cè)，且具有較強(qiáng)的頻率測(cè)量能力。

表1 4種方法的檢測(cè)結(jié)果

表2 4種方法的基波檢測(cè)相對(duì)誤差

由圖5可以看出，當(dāng)NR=20 dB時(shí)，因噪聲混入，3次諧波的頻率曲線波動(dòng)加劇，210 Hz的頻率分量受此影響很大，其頻率曲線出現(xiàn)波動(dòng)和間斷的情況，但基波的頻率曲線未受明顯影響，表明SWT在噪聲混入時(shí)仍能夠檢測(cè)出基波，具有抗噪性。對(duì)比SWT與EMD、CEEMDAN、VMD所分解的不同項(xiàng)模態(tài)分量波形，分析比較4種方法抗噪的能力，4種方法的模態(tài)波形如圖6所示。

由圖6可知，噪聲混入后，EMD和CEEMDAN的虛假分量變多，代表基波的IMF分量延后出現(xiàn)，它們的 IMF4對(duì)應(yīng)基波分量，但存在嚴(yán)重的模態(tài)混疊，不能完全提取出基波；而VMD分解出的IMF1分量即對(duì)應(yīng)基波，但相比SWT分解出的基波IMT1，VMD的整體波形不整齊，存在一定程度的模態(tài)混疊，且圖6(c)中后3個(gè)IMF分量的部分頻率和幅值信息也已經(jīng)被噪聲淹沒(méi)；而SWT不僅能較為完整地提取出基波，且其后3個(gè)IMT分量的頻率分別為90.000 022 Hz、149.708 364 Hz和209.999 452 Hz，表明在噪聲混入時(shí)，該方法對(duì)諧波、間諧波也具有識(shí)別能力，突出了SWT較好的抗噪能力。

圖5 SWT時(shí)頻分析(SNR=20 dB)

圖6 模態(tài)分解結(jié)果(SNR=20 dB)

2.2 基波頻率波動(dòng)的仿真信號(hào)

設(shè)置采樣時(shí)長(zhǎng)為1.0 s，采樣頻率為4096 Hz，截取基波在前0.12 s的頻率波動(dòng)情況，如圖7所示。

圖7 基頻波動(dòng)(0～0.12 s)

由圖7可以看出，圖7(a)箭頭指在基波波形的0.0591 s時(shí)刻位置，由圖7(b)中基波在理想狀態(tài)下的時(shí)間-頻率特性曲線表明，該基波頻率波動(dòng)的信號(hào)分量在前0.0591 s時(shí)段內(nèi)存在50.2 Hz的基波頻率成分，在0.0591～1.0000 s時(shí)段存在50.0 Hz的基波頻率成分，基波頻率在前0.0591時(shí)段內(nèi)存在波動(dòng)。

分別用EMD、CEEMDAN、VMD和SWT對(duì)含有基波頻率波動(dòng)的暫態(tài)信號(hào)進(jìn)行基波檢測(cè)，所分解出的模態(tài)分量如圖8所示。前6個(gè)周期(圖8中箭頭所截取的周期數(shù))提取出的分量與原始基波的波形吻合效果如圖9所示。

由圖8和圖9可知，EMD、CEEMDAN、VMD和SWT都能分解出基波，每種方法的IMF2、IMF3、IMF1、IMT1對(duì)應(yīng)基波，但I(xiàn)MT1分量波形受基波波動(dòng)導(dǎo)致的非同步采樣影響更小，SWT和VMD所分解出的3次諧波在起止時(shí)刻存在抖動(dòng)，SWT抖動(dòng)更輕微；VMD所分解出的基波存在邊界效應(yīng)，起止峰值小于原始基波峰值，CEEMDAN和EMD出現(xiàn)了不同程度的模態(tài)混疊，峰值處存在不同大小的幅值，整體波形不整齊，基波波形吻合效果較差。

圖8 模態(tài)分解結(jié)果(基頻波動(dòng))

圖9 基波波形吻合效果2

2.3 間諧波靠近基波和諧波的仿真信號(hào)

由圖10可知，加入20 dB的高斯白噪聲，EMD和CEEMDAN所分解的基波分量IMF4出現(xiàn)了不同程度的模態(tài)混疊，不能完全提取基波；VMD分解得到的基波分量IMF1也存在較明顯的模態(tài)混疊，且左右端點(diǎn)處存在邊界效應(yīng)，其后3個(gè)模態(tài)分量信息已經(jīng)被噪聲淹沒(méi)；而SWT分解的分量按能量從大到小排列，基波能量最大，所以第1個(gè)分量IMT1為基波，不像EMD系列是從高頻到低頻排列，必須對(duì)所有分解項(xiàng)進(jìn)行參數(shù)分析，直至找到基波分量，SWT所分解的IMT1分量的頻率為49.999 990 Hz。此外，在諧波和間諧波出現(xiàn)模態(tài)混疊的情形下，SWT依然能準(zhǔn)確提取基波，表明其基波提取能力強(qiáng)于諧波和間諧波。4種方法抗噪的能力由強(qiáng)到弱排序?yàn)镾WT＞VMD＞CEEMDAN＞EMD。

圖10 模態(tài)分解結(jié)果(SNR=20 dB)

3 ?實(shí)驗(yàn)驗(yàn)證與分析

根據(jù)某牽引供電系統(tǒng)高壓側(cè)220 kV的A相瞬時(shí)電壓實(shí)測(cè)數(shù)據(jù)，采樣頻率為20 kHz，采樣時(shí)間間隔為50ms，信號(hào)長(zhǎng)度為1 s，取其前0.2 s的信號(hào)長(zhǎng)度，將SWT與EMD、CEEMDAN、VMD應(yīng)用于實(shí)際信號(hào)，分解得到前4項(xiàng)模態(tài)分量波形如圖11所示。

由圖11可知，4種方法均能分解出基波，EMD、CEEMDAN、VMD、SWT分解得到的IMF4、IMF4、IMF1、IMT1分別對(duì)應(yīng)基波，但EMD和CEEMDAN所分解出的第4個(gè)分量代表基波，存在分解分量的延遲，而SWT和VMD分解出的第1個(gè)分量即代表基波，基波檢測(cè)效率更高；SWT分解出的IMT3和IMT4的頻率分別為145.0145 Hz和250.0054 Hz，對(duì)應(yīng)3次諧波和5次諧波，而剩余由其他3種方法分解出的模態(tài)分量則近似噪聲。

圖11 模態(tài)分解結(jié)果(實(shí)測(cè)電壓)

經(jīng)Hilbert變換，得到4種方法的基波時(shí)間-頻率特性曲線如圖12所示。

由圖12可知，SWT方法提取到的基波頻率波動(dòng)較其他3種方法更小，基波頻率波動(dòng)范圍控制在49.7～50.3 Hz，波動(dòng)幅度小于1 Hz，而其他3種方法的基波頻率波動(dòng)幅度則超過(guò)了4 Hz，表明本文方法適用于實(shí)際信號(hào)的頻率檢測(cè)。

圖12 基波時(shí)間-頻率特性曲線

EMD、CEEMDAN、VMD和SWT的基波測(cè)量頻率、標(biāo)準(zhǔn)差及其與原始基波頻率標(biāo)準(zhǔn)值(50 Hz)之間所存在的相對(duì)誤差如表3所示。

表3 4種方法的檢測(cè)結(jié)果(實(shí)測(cè)電壓)

由表3可知，實(shí)際A相電壓信號(hào)經(jīng)本文方法測(cè)量得到的基波頻率與原始基波頻率的相對(duì)誤差5.2427×10-7較其他3種方法更小，精度達(dá)到10-7量級(jí)，其標(biāo)準(zhǔn)差為0.036 553，較其他3種方法的標(biāo)準(zhǔn)差也更小，測(cè)量結(jié)果更穩(wěn)定，表明SWT能更為準(zhǔn)確地從實(shí)際A相電壓信號(hào)中檢測(cè)出基波頻率，另外，該方法測(cè)量出的基波幅值為1.872 433×105V。

4 ?結(jié)論

本文在諧波幅值瞬變、噪聲混入、基波頻率波動(dòng)、間諧波頻率靠近基波和諧波的情況下，基于SWT進(jìn)行電力信號(hào)的基波提取，再通過(guò)Hilbert變換計(jì)算出基波頻率和幅值，從而完成基波檢測(cè)。

實(shí)驗(yàn)結(jié)果表明，在諧波幅值發(fā)生階躍變化時(shí)， SWT通過(guò)同步擠壓來(lái)避免頻率曲線交叉，能夠準(zhǔn)確提取出基波分量，模態(tài)分解能力強(qiáng)，可以有效避免模態(tài)混疊；對(duì)于混有噪聲的暫態(tài)信號(hào)，該方法仍能準(zhǔn)確地測(cè)量基波特征參數(shù)，表明SWT具有較強(qiáng)的抗噪性，且基波的提取能力強(qiáng)于諧波和間諧波；本文也考慮了因基波頻率波動(dòng)導(dǎo)致采樣不同步的情況，將所采用的SWT與EMD、CEEMDAN、VMD相對(duì)比，得到SWT測(cè)量出的基波頻率和幅值受采樣不同步的影響更小，與原始基波的波形吻合度更高的結(jié)論。

綜上，SWT能準(zhǔn)確地測(cè)量出基波參數(shù)，其分解出的第1個(gè)IMT分量即是基波，具有更強(qiáng)的抗噪能力，基波檢測(cè)精度高于傳統(tǒng)的時(shí)頻分析方法，如EMD、CEEMDAN和VMD，頻率精度最高可達(dá)10-8量級(jí)，且SWT基波提取的能力強(qiáng)于諧波和間諧波，是提取基波的有力工具。

[1] 秦毅, 秦樹(shù)人, 毛永芳. 小波變換中經(jīng)驗(yàn)?zāi)B(tài)分解的基波檢測(cè)及其在機(jī)械系統(tǒng)中的應(yīng)用[J]. 機(jī)械工程學(xué)報(bào), 2008, 56(3): 135-142.

QIN Yi, QIN Shuren, MAO Yongfang. Fundamental wave detection based on wavelet transform and empirical mode decomposition with application in mechanical system[J]. Journal of Mechanical Engineering, 2008, 56(3): 135-142.

[2] 蘭金晨, 朱革蘭, 李松奕. 基于廣義S變換突變識(shí)別的暫穩(wěn)態(tài)諧波檢測(cè)方法[J]. 電力系統(tǒng)保護(hù)與控制, 2020, 48(23): 40-46.

LAN Jinchen, ZHU Gelan, LI Songyi. Detection method of transient and steady state harmonics based on generalized S-transform mutation recognition[J]. Power System Protection and Control, 2020, 48(23): 40-46.

[3] 李寧, 左培麗, 王新剛, 等. 基于改進(jìn)DFT和時(shí)域準(zhǔn)同步的間諧波檢測(cè)算法[J]. 電力自動(dòng)化設(shè)備, 2017, 37(4): 170-178.

LI Ning, ZUO Peili, WANG Xingang, et al. Inter-harmonic detection based on improved DFT and time-domain quasi-synchronization[J]. Electric Power Automation Equipment, 2017, 37(4): 170-178.

[4] 孫仲民, 何正友, 臧天磊. 基于滑窗頻譜分離算法的基波參數(shù)測(cè)量[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2017, 37(19): 5604-5612, 5837.

SUN Zhongmin, HE Zhengyou, ZANG Tianlei. Estimating fundamental parameters by using sliding-window spectrum separation algorithm[J]. Proceedings of the CSEE, 2017, 37(19): 5604-5612, 5837.

[5] CHEN G L, TSAI S H, YANG K J. On performance of sparse fast Fourier transform and enhancement algorithm[J]. IEEE Transactions on Signal Processing, 2017, 65(21): 5716-5729.

[6] 陳至豪, 王立德, 王沖, 等. 基于組合余弦優(yōu)化窗四譜線插值FFT的電力諧波分析方法[J]. 電網(wǎng)技術(shù), 2020, 44(3): 1105-1113.

CHEN Zhihao, WANG Lide, WANG Chong, et al. An approach for electrical harmonic analysis based on optimized composite cosine window four-spectrum-line interpolation FFT[J]. Power System Technology, 2020, 44(3): 1105-1113.

[7] 王燕. 電能質(zhì)量擾動(dòng)檢測(cè)的研究綜述[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(13): 174-186.

WANG Yan. Review of research development in power quality disturbance detection[J]. Power System Protection and Control, 2021, 49(13): 174-186.

[8] DESAI J P, MAKWANA V H. A novel out of step relaying algorithm based on wavelet transform and a deep learning machine model[J]. Protection and Control of Modern Power Systems, 2021, 6(4): 500-511.

[9] VINAYAK B A, ANAND K A, JAGADANAND G. Wavelet-based real-time stator fault detection of inverter-fed induction motor[J]. IET Electric Power Applications, 2020, 14(1): 82-90.

[10] HUANG N E, SHEN Z, LONG S R. New view of nonlinear water waves: the Hilbert spectrum[J]. Annual Review of Fluid Mechanics, 1999, 31: 417-457.

[11] Lü Chenhuan, ZHAO Jun, WU Chao, et al. Optimization of the end effect of Hilbert-Huang transform (HHT)[J]. Chinese Journal of Mechanical Engineering, 2017, 30(3): 732-745.

[12] WU Z, HUANG N E. Ensemble empirical mode decomposition: a noise-assisted data analysis method[J]. Advances in Adaptive Data Analysis, 2009, 1(1): 1-41.

[13] 溫志平, 方江雄, 劉軍, 等. 基于CEEMDAN和SWT的地震信號(hào)隨機(jī)噪聲壓制[J]. 地球物理學(xué)進(jìn)展, 2019, 34(3): 949- 960.

WEN Zhiping, FANG Jiangxiong, LIU Jun, et al. Seismic random noise attenuation based on CEEMDAN and SWT[J]. Progress in Geophysics, 2019, 34(3): 949-960.

[14] 張慧, 孟克其勞, 張占強(qiáng), 等. 基于CEEMDAN和HHT的諧波檢測(cè)方法[J]. 電子測(cè)量技術(shù), 2020, 43(6): 83-87.

ZHANG Hui, MENG Keqilao, ZHANG Zhanqiang, et al. Harmonic detection method based on CEEMDAN and HHT[J]. Electronic Measurement Technology, 2020, 43(6): 83-87.

[15] DRAGOMIRETSKIY K, ZOSSO D. Variational mode decomposition[J]. IEEE Transactions on Signal Processing, 2014, 62(3): 531-544.

[16] 郎憲明, 朱永強(qiáng), 李金娜, 等. 基于時(shí)頻分析和改進(jìn)差分進(jìn)化的多點(diǎn)泄漏定位方法[J]. 振動(dòng)與沖擊, 2022, 41(1): 39-45, 87.

LANG Xianming, ZHU Yongqiang, LI Jinna, et al. Multi-point leakage locating method based on TFA and improved DE[J]. Journal of Vibration and Shock, 2022, 41(1): 39-45, 87.

[17] DAUBECHIES I, LU J, WU H. Synchrosqueezed wavelet transforms: an empirical mode decomposition-like tool[J]. Applied and Computational Harmonic Analysis, 2010, 30(2): 243-261.

[18] 喻敏, 王斌, 王文波, 等. 基于SST的間諧波檢測(cè)方法[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2016, 36(11): 2944-2951.

YU Min, WANG Bin, WANG Wenbo, et al. An inter- harmonic detection method based on synchrosqueezing wavelet transform[J]. Proceedings of the CSEE, 2016, 36(11): 2944- 2951.

[19] 劉博文, 張新燕, 常喜強(qiáng), 等. 基于SWT的風(fēng)電匯集地區(qū)非同步振蕩PMU暫態(tài)錄波數(shù)據(jù)分析[J]. 太陽(yáng)能學(xué)報(bào), 2020, 41(1): 217-224.

LIU Bowen, ZHANG Xinyan, CHANG Xiqiang, et al. Analysis of non-synchronous oscillation based on SWT and transient wave recorded data of PMU in wind farm integration area[J]. Acta Energiae Solaris Sinica, 2020, 41(1): 217-224.

[20] 熊炘, 占銳, 王小靜. 同步壓縮小波與希爾伯特-黃變換性能對(duì)比[J]. 振動(dòng)·測(cè)試與診斷, 2015, 35(6): 1103-1109.

XIONG Xin, ZHAN Rui, WANG Xiaojing. Performance comparison between synchronous compressed wavelet and Hilbert-Huang Transform[J]. Vibration·Testing and Diagnosis, 2015, 35(6): 1103-1109.

[21] 喻敏, 王斌, 王文波, 等. 基于同步擠壓小波變換的電力系統(tǒng)時(shí)變諧波檢測(cè)[J]. 電工技術(shù)學(xué)報(bào), 2017, 32(增刊1): 50-57.

YU Min, WANG Bin, WANG Wenbo, et al. Power system time-varying transient harmonics detection based on SWT[J]. Transactions of China Electrotechnical Society, 2017, 32(S1): 50-57.

[22] 劉景良, 任偉新, 王佐才, 等. 基于同步擠壓小波變換的結(jié)構(gòu)瞬時(shí)頻率識(shí)別[J]. 振動(dòng)與沖擊, 2013, 32(18): 37-42, 48.

LIU Jingliang, REN Weixin, WANG Zuocai, et al. Instantaneous frequency identification based on synchrosqueezing wavelet transformation[J]. Journal of Vibration and Shock, 2013, 32(18): 37-42, 48.

[23] IATSENKO D, MCCLINTOCK P V E, STEFANOVSKA A. Nonlinear mode decomposition: a noise-robust, adaptive decomposition method[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 2015, 92(3): 1-25.

[24] 李欣悅, 李鳳婷, 尹純亞, 等. 直流雙極閉鎖故障下送端系統(tǒng)暫態(tài)過(guò)電壓計(jì)算方法[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(1): 1-8.

LI Xinyue, LI Fengting, YIN Chunya, et al. Transient overvoltage calculation method of HVDC sending-end system under DC bipolar blocking[J]. Power System Protection and Control, 2021, 49(1): 1-8.

[25] 毛宇陽(yáng), 楊志. V2X雙向充電機(jī)在微電網(wǎng)中的控制策略[J]. 新能源汽車(chē)供能技術(shù), 2021, 5(1): 19-22.

MAO Yuyang, YANG Zhi. Control strategy of V2X bidirectional charger in microgrid[J]. Energy Supply Techniques of New Energy Vehicles, 2021, 5(1): 19-22.

[26] 馬也馳, 陳隆道. 應(yīng)用于頻率波動(dòng)電網(wǎng)的改進(jìn)相位差校正法[J]. 電力系統(tǒng)保護(hù)與控制, 2018, 46(9): 136-145.

MA Yechi, CHEN Longdao. An improved phase correction difference method applied to the power systems with frequency fluctuation[J]. Power System Protection and Control, 2018, 46(9): 136-145.

[27] 電能質(zhì)量電力系統(tǒng)頻率偏差: GB/T 15945—2008[S]. 北京: 中國(guó)標(biāo)準(zhǔn)出版社, 2008.

Power quality-frequency deviation for power system: GB/T 15945—2008[S]. Beijing: Standards Press of China, 2008.

[28] 倪大志. 電源質(zhì)量影響儀表測(cè)量準(zhǔn)確度[J]. 計(jì)量與測(cè)試技術(shù), 2010, 37(12): 53-55.

NI Dazhi. Power quality affects the meter measuring accuracy[J]. Metrology & Measurement Technique, 2010, 37(12): 53-55.

[29] XU D, LIU P F, CHEN Z P. Damage mode identification and singular signal detection of composite wind turbine blade using acoustic emission[J]. Composite Structures, 2021, 255: 1-16.

[30] 赫玉瑩, 王學(xué)華, 劉福鑫. LCL型并網(wǎng)逆變器分裂電流單環(huán)控制及其強(qiáng)魯棒性權(quán)值設(shè)計(jì)[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2021, 41(10): 3536-3546, 3678.

HE Yuying, WANG Xuehua, LIU Fuxin. Splitting current single-loop control and robust weight value design for LCL-type grid-connected inverter[J]. Proceedings of the CSEE, 2021, 41(10): 3536-3546, 3678.

Fundamental detection for a power system based on SWT

TAO Jialan1, 2, YU Min1, 2, CHEN Guici1, 2, WANG Bin3

(1.Hubei Province Key Laboratory of Systems Science in Metallurgical Process, Wuhan University of Science and Technology,Wuhan 430065, China; 2. College of Science, Wuhan University of Science and Technology, Wuhan 430065, China; 3. College of Information and Engineering, Wuhan University of Science and Technology, Wuhan 430081, China)

When noise is mixed into a signal containing the fundamental component, the traditional time-frequency analysis method is prone to modal aliasing in the fundamental component extraction process. To accurately detect the fundamental component, synchrosqueezing wavelet transform (SWT) with high accuracy of time-frequency analysis is applied to realize the fundamental detection. First, the signal containing the fundamental component is decomposed into a set of intrinsic mode type functions (IMTs) by SWT. The first component IMT1 represents the fundamental. Then, the fundamental frequency and amplitude are measured by Hilbert Transform. Finally, the algorithm is verified in the situation of harmonic amplitude transient, noise mixing, fundamental frequency fluctuation and interharmonic frequencies near the fundamental and harmonic frequencies. The experimental results show that SWT can accurately extract the fundamental component, and the frequency accuracy can reach 10-8orders of magnitude. The method has strong noise resistance and is better at extracting the fundamental component than harmonics and interharmonics.

fundamental detection; mode aliasing; SWT; Hilbert transform; noise resistance

10.19783/j.cnki.pspc.211660

2021-12-06；

2022-03-08

陶佳蘭(1997—)，女，碩士研究生，研究方向?yàn)殡娔苜|(zhì)量擾動(dòng)檢測(cè)；E-mail: jlantao_21@163.com

喻 敏(1975—)，女，通信作者，博士，講師，研究方向?yàn)榉中闻c小波應(yīng)用、電能質(zhì)量擾動(dòng)檢測(cè)等；E-mail: yufeng3378@163.com

陳貴詞(1978—)，男，博士，教授，研究方向?yàn)殡S機(jī)系統(tǒng)、魯棒控制等。E-mail: gcichen@yahoo.com.cn

國(guó)家自然科學(xué)基金項(xiàng)目資助(51877161，61671338)；冶金工業(yè)過(guò)程系統(tǒng)科學(xué)湖北省重點(diǎn)實(shí)驗(yàn)室基金重點(diǎn)項(xiàng)目資助(Y202007，Z201901)

This work is supported by the National Natural Science Foundation of China (No. 51877161 and No. 61671338).

(編輯 許 威)