張銳+張恒煜+周永陽(yáng)+丁鵬

摘 要：短時(shí)閃變嚴(yán)重度Pst是衡量電能質(zhì)量的一個(gè)重要指標(biāo)，為提高Pst的計(jì)算精度，提出了基于原子稀疏分解和實(shí)數(shù)編碼量子進(jìn)化算法（AtomRQEA）相結(jié)合的短時(shí)閃變嚴(yán)重度計(jì)算方法。其核心是根據(jù)電壓波動(dòng)與閃變信號(hào)的特點(diǎn)構(gòu)造了原子庫(kù)，由實(shí)數(shù)編碼量子進(jìn)化算法優(yōu)化原子特征參數(shù)，自適應(yīng)選擇最佳匹配原子重構(gòu)電壓波動(dòng)信號(hào)，進(jìn)而計(jì)算短時(shí)閃變嚴(yán)重度Pst。仿真結(jié)果表明：基于AtomRQEA方法，能夠提高短時(shí)閃變嚴(yán)重度Pst計(jì)算精度，驗(yàn)證了提出方法的有效性和適用性。

關(guān)鍵詞：電能質(zhì)量；短時(shí)閃變嚴(yán)重度；原子分解

DOI：10.15938/j.jhust.2017.06.009

中圖分類(lèi)號(hào)： TM866

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

文章編號(hào)： 1007-2683（2017）06-0046-06

Abstract：Short term flicker severity Pst is an important index of evaluating power quality. To improve precision of calculation AtomRQEA is proposed. The core is to construct atom library in terms of characteristic of voltage fluctuation and flicker， using RQEA（Realcoded Quantum Evolutionary Algorithm） to optimize characteristic parameters of atom， choosing selfadaption optimal matching atom to reconstruct voltage signal， and calculate flicker severity Pst. Simulation results show the method proposed in this paper can improve the precise of calculating short term flicker severity Pst， and prove it is available and applicable.

Keywords：power quality； short term flicker severity； atomic decomposition

0 引 言

目前我國(guó)堅(jiān)強(qiáng)智能電網(wǎng)的建設(shè)在不斷的推進(jìn)，光伏發(fā)電、風(fēng)力電源等再生能源的不斷接入對(duì)現(xiàn)有配電網(wǎng)結(jié)構(gòu)產(chǎn)生巨大的沖擊，電力系統(tǒng)正面臨著新的電能質(zhì)量問(wèn)題，對(duì)電能質(zhì)量進(jìn)行實(shí)時(shí)監(jiān)測(cè)、評(píng)估，為及時(shí)獲得現(xiàn)場(chǎng)數(shù)據(jù)，采取合理的應(yīng)對(duì)措施從而全面提高電網(wǎng)的安全性、經(jīng)濟(jì)性、適應(yīng)性具有重要意義。

IEC標(biāo)準(zhǔn)中短時(shí)閃變嚴(yán)重度Pst是衡量電壓波動(dòng)與閃變程度的主要綜合指標(biāo)[1]，通常采用快速傅里葉變換（FFT）提取電壓波動(dòng)信號(hào)的特征量，再計(jì)算Pst。然而，該方法存在頻譜泄露和柵欄效應(yīng)，導(dǎo)致計(jì)算精度較低，且實(shí)際應(yīng)用效果差。文[2-4]提出采用修正FFT算法提取電壓波動(dòng)信號(hào)的特征量，用插值FFT克服FFT缺點(diǎn)。但是，該方法不能從根本上消除頻譜泄露和刪欄效應(yīng)的影響，很難改善計(jì)算精度，且插值FFT易于受到噪聲的影響。

近年來(lái)，文[5-7]提出了原子分解方法，該方法是一種時(shí)頻分析法，通過(guò)自適應(yīng)選取過(guò)完備原子庫(kù)中的原子分解信號(hào)，實(shí)現(xiàn)對(duì)信號(hào)的表示。被廣泛地應(yīng)用在信號(hào)處理領(lǐng)域。原子分解過(guò)程實(shí)質(zhì)上是使用一系列的相關(guān)信號(hào)估計(jì)原始信號(hào)，且原始信號(hào)可以以更少的計(jì)算量重構(gòu)原始信號(hào)。信號(hào)分解采用匹配跟蹤（MP）算法[8-9]，也就是說(shuō)，每一次迭代，從原子庫(kù)選擇原子匹配信號(hào)的剩余分量，與信號(hào)最匹配原子的線(xiàn)性組合來(lái)描述剩余分量。但面對(duì)離散空間的搜索問(wèn)題，MP幾乎無(wú)法找到全局最優(yōu)解。文中提出原子稀疏分解與實(shí)數(shù)編碼量子進(jìn)化算法（atomic decomposition and realcoded quantum evolutionary algorithm，AtomRQEA）相結(jié)合的短時(shí)閃變嚴(yán)重度計(jì)算方法，主要是利用實(shí)數(shù)編碼量子進(jìn)化算法自適應(yīng)選擇與信號(hào)最佳匹配的原子，再利用最佳匹配原子的線(xiàn)性組合重構(gòu)電壓波動(dòng)與閃變信號(hào)的包絡(luò)線(xiàn)，進(jìn)而計(jì)算短時(shí)閃變嚴(yán)重度Pst。通過(guò)仿真試驗(yàn)表明：提出AtomRQEA方法能夠提高短時(shí)閃變嚴(yán)重度Pst的計(jì)算精度。

1 原子分解

與傳統(tǒng)意義上的信號(hào)分解不同，原子分解屬于非正交分解[10-12]，在信號(hào)空間中構(gòu)造了足夠密集的基即：原子，這些原子構(gòu)成的庫(kù)是過(guò)完備的，足以自適應(yīng)的選擇最匹配的原子來(lái)表示任意信號(hào)，從而獲得信號(hào)的稀疏表示，所選擇的原子體現(xiàn)了信號(hào)的特征。

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

本文仿真實(shí)驗(yàn)數(shù)據(jù)采用的是IEC提供的標(biāo)準(zhǔn)數(shù)據(jù)，分別是在單位瞬時(shí)閃變值時(shí)的正弦波波動(dòng)電壓波動(dòng)值和矩形波波動(dòng)電壓波動(dòng)值，在一定的電壓波動(dòng)條件下，即電壓波形、波動(dòng)的幅度、頻率不變，計(jì)算得到的Pst為0.7144。采樣頻率為1600Hz，采樣長(zhǎng)度為1s。為了比較文中方法的有效性，分別采用離散傅里葉分析方法和RQEAAtom方法對(duì)波動(dòng)電壓進(jìn)行頻譜分析。

圖1、圖2、圖3分別給出了利用RQEAAtom對(duì)頻率為8.8Hz正弦波、多頻閃變和8Hz矩形波波動(dòng)電壓的分解過(guò)程。電壓閃變現(xiàn)象對(duì)頻率為8.8Hz，正弦電壓波動(dòng)為0.25%，8Hz矩形波電壓波動(dòng)0.201%時(shí)較敏感。圖4、圖5對(duì)兩種算法計(jì)算Pst的相對(duì)誤差進(jìn)行了對(duì)比。

由圖1可以發(fā)現(xiàn)，利用RQEAAtom對(duì)8.8Hz正弦波包絡(luò)線(xiàn)信號(hào)進(jìn)行分解，經(jīng)過(guò)1次原子分解，得到包絡(luò)線(xiàn)信號(hào)的幅值和頻率估計(jì)值分別1.249999mV 和8.799998Hz，接近理論值。圖2原子第一次分解可得到幅度4.41mV，頻率1.99997Hz閃變信號(hào)，原子第2次分解時(shí)得到第二個(gè)頻率分量，測(cè)量值幅度A=1.639997 頻率f=6.00006Hz，很接近理論值。圖2（e）經(jīng)第1次和第2次分解后的殘余信號(hào)，利用文中方法能夠以很高的精度重構(gòu)正弦波動(dòng)閃變信號(hào)。endprint

圖3為矩形波閃變電壓包絡(luò)信號(hào)，矩形波信號(hào)包含了奇次諧波分量，由于引起閃變現(xiàn)象的電壓波動(dòng)頻率范圍是0.5-35Hz，這里提取了基波和3次諧波分量。由原子第一次分解獲得的電壓波動(dòng)幅度1.279553，理論值為1.279605，相對(duì)誤差為0.004063%；頻率為8Hz，相對(duì)誤差為0%。第2次原子分解得到了3次諧波分量幅值為0.426421，理論值為0.426535，相對(duì)誤差0.026726%；頻率為23.99998Hz，理論值為24Hz，相對(duì)誤差8.3×10-5%。通過(guò)確定閃變信號(hào)的頻域信息，計(jì)算Pst，與DFT算法相比較，計(jì)算Pst的結(jié)果如圖4、圖5所示，給出了計(jì)算Pst相對(duì)誤差曲線(xiàn)。

從圖4、圖5可以看出文中采用的AtomRQEA方法提高了頻域的分析能力，從而使計(jì)算Pst的相對(duì)誤差較小，且在對(duì)不同的頻率下計(jì)算的Pst相對(duì)誤差結(jié)果相差不大，而DFT方法在計(jì)算Pst時(shí)在不同的頻率下，相對(duì)誤差出現(xiàn)較大的振蕩?？梢?jiàn)，文中提出方法具有很好的性能。今后工作主要集中在實(shí)數(shù)編碼量子進(jìn)化算法和原子分解原理研究，更好的結(jié)合二者的優(yōu)點(diǎn)，擴(kuò)大基于原子分解和實(shí)數(shù)編碼量子進(jìn)化算法在電力系統(tǒng)中應(yīng)用的廣度和深度。

5 結(jié) 語(yǔ)

文中將原子稀疏分解方法與實(shí)數(shù)編碼量子進(jìn)化算法相結(jié)合用于分析電壓波動(dòng)與閃變信號(hào)，從而計(jì)算電力系統(tǒng)短時(shí)閃變嚴(yán)重度Pst。仿真試驗(yàn)表明，采用所提方法計(jì)算Pst具有很高的計(jì)算精度，且非常穩(wěn)定。基于原子分解和實(shí)數(shù)編碼量子進(jìn)化算法不僅能夠?yàn)殡娏ο到y(tǒng)Pst計(jì)算提供了新途徑，也為相關(guān)領(lǐng)域類(lèi)似問(wèn)題的解決提供了一個(gè)工具。

參 考 文 獻(xiàn)：

[1] 趙剛， 施圍， 林海雪. 閃變值計(jì)算方法的研究[J]. 電網(wǎng)技術(shù)， 2001， 25（11）：15-18.

[2] 李志剛，汪濤，何怡剛，等.基于改進(jìn)余弦窗的加窗插值FFT諧波分析[J].電氣工程學(xué)報(bào)，2016，11（4）：15-21.

[3] 陳子珍， 夏冰冰，閻威武. 基于改進(jìn)加窗插值FFT的高精度諧波與間諧波檢測(cè)算法[J].中國(guó)電力，2015，48（9）：73-79.

[4] 郝柱，顧偉，褚建新，等.基于四譜線(xiàn)插值FFT的電網(wǎng)諧波檢測(cè)方法[J].電力系統(tǒng)保護(hù)與控制，2014， 42（19）：107-113.

[5] 張磊，黃純，江亞群，等. 應(yīng)用HS改進(jìn)原子分解的電能質(zhì)量擾動(dòng)辨識(shí)分析[J]. 電網(wǎng)技術(shù)， 2015，（01）：194-201.

[6] JIA Q. Q， YAO R.， WANG N.， MENG R. L. An Approach to Detect Harmonics/Interharmonics with Closing Frequencies Using Atomic Decomposition and Windowed Frequency Shifting Algorithm [J]. Proceedings of the CSEE， 2014， 34（27）：4605-4612.

[7] 于浩明，黃純，江亞群等.采用快速原子分解的電能質(zhì)量擾動(dòng)信號(hào)參數(shù)辨識(shí)[J].電網(wǎng)技術(shù)，2014，38（8）：2237-2243.

[8] ZHANG Jun，ZENG Pingping，MA Jianhua. Classification of Power Quality Disturbances Using PSOMP and Parametric Dictionaries[C]// 2015 International Conference on Intelligent Computing and Internet of Things， 2015：21-25.

[9] YU FENG QIN， CHEN YING. Musical Instrument Classification Based on Improved Matching Pursuit with InstrumentSpecific Atoms[C]// 2015 IIAI 4th International Congress on Advanced Applied Infformatics，2015：506-510.

[10]SUNIL BELUR NAGARAJ， NATHAN J. STEVENSON，WILLIAM P. MARNANE， et al. Neonatal Seizure Detection Using Atomic Decomposition with a Novel Dictionary [J]. IEEE Transactions on Biomedical engineering， 2014， 61 （11）：2724-2732.

[11]曲正偉，郝婉茹，王寧.原子分解快速算法在電能質(zhì)量擾動(dòng)分析中的應(yīng)用[J].電力自動(dòng)化設(shè)備，2015，35（10）：145-150.

[12]HE Qingbo， SONG Haiyue， DING Xiaoxi. Sparse Signal Reconstruction Based on Time – frequency Manifold for Rolling Element Bearing Fault Signature Enhancement[J].IEEE Transaction on Instrumentation and Measurement， 2016， 65（2）： 482-491.

[13]賈秀芳，岳娜，安海清.基于間諧波電壓的閃變傳遞計(jì)算方法[J].華北電力大學(xué)學(xué)報(bào)，2015，42（4）： 15-21.

[14]李蘭芳，馬明，盛超，等.基于旋轉(zhuǎn)因子變換插值的電壓閃變檢測(cè)算法[J].電工技術(shù)學(xué)報(bào)， 2016，31 （22）：178-185.

[15]楊爍，曹思揚(yáng)，戴朝華，等.電能質(zhì)量擾動(dòng)信號(hào)時(shí)頻原子分解的進(jìn)化匹配追蹤算法[J].電力系統(tǒng)保護(hù)與控制，2015，43（16）：80-86.

[16]LI MING. An Adaptive Quantum Evolutionary Algorithm and its Application to Path Planning [J].2015IEEE International Conference on System， Man， and Cybermetics，2015，2067-2017.

[17]CHE ADA， WU PENG，CHU FENG， ZHOU MENG CHU. Improved QuantumInspired Evolutionary Algorithm for LargeSize Lane Reservation [J].IEEE Transactions on System， Man， and Cybernetics： System， 2015， 45 （12）： 1535-1548.

[18]劉振，胡云安，彭軍，協(xié)同進(jìn)化擴(kuò)展緊致量子進(jìn)化算法[J].控制與決策，2014，29（2）：320-326.

[19]Flicker MeterFunctional and Design Specifications [S]. EN61000-4-15，2003-01.120

[20]孫樹(shù)勤. 電壓波動(dòng)與閃變[M]. 北京：中國(guó)電力出版， 1999：41-55.

（編輯：王 萍）endprint