一種提高電網(wǎng)諧波測量精度的方法研究

薛軍

摘要：本文介紹了一種提高電能質(zhì)量中諧波精度的方法。步驟包括：（1）A/D 采樣獲得等間隔的原始數(shù)據(jù)序列α。（2）根據(jù)當(dāng)前頻率計算新的采樣間隔1/Fn。（3）在數(shù)據(jù)序列α的基礎(chǔ)上通過二階插值重構(gòu)新的數(shù)據(jù)序列β。（4）對新的采樣序列β進(jìn)行FFT變換得到各次諧波的有效值，同時可以獲得各次諧波畸變率、波形失真度等相關(guān)數(shù)據(jù)。該方法無需調(diào)整采樣頻率即可得到同步信號，節(jié)約了硬件和軟件的資源開銷，簡化了系統(tǒng)的設(shè)計。

關(guān)鍵詞：諧波測量;電能質(zhì)量;二階插值重構(gòu)

中圖分類號：TM935 文獻(xiàn)標(biāo)識碼：A 文章編號：1007-9416（2018）11-0054-02

0 引言

近三十年來，各種電力電子設(shè)備的快速發(fā)展導(dǎo)致電網(wǎng)的諧波污染日益嚴(yán)重。諧波的增加影響了電能生產(chǎn)、傳輸和利用效率，同時使得電氣設(shè)備過熱、產(chǎn)生振動和噪聲，縮減壽命，嚴(yán)重時可引起電力系統(tǒng)局部并聯(lián)諧振或串聯(lián)諧振，使諧波含量放大，造成電容器等設(shè)備損壞。由于諧波污染造成的故障和事故不斷發(fā)生，諧波的危害逐漸引起關(guān)注，諧波測量成為電能質(zhì)量中一項(xiàng)非常重要的課題[1-2]。

傳統(tǒng)的諧波測量為避免系統(tǒng)的頻率變化引起的“柵欄效應(yīng)”帶來的誤差，通常的做法是實(shí)時根據(jù)當(dāng)前頻率去重新設(shè)置采樣率，從而準(zhǔn)確計算出諧波分量，同時也帶來了軟件和硬件的額外開銷，使得系統(tǒng)更加復(fù)雜[3-4]。

1 技術(shù)方案

本文提出了一種通過重構(gòu)采樣數(shù)據(jù)，無需實(shí)時的設(shè)置采樣率計算諧波的方法。該諧波計算方法可以降低電力系統(tǒng)頻率的變化導(dǎo)致諧波測量的誤差，提高算法精度。解決了傳統(tǒng)方案中軟硬設(shè)計過于復(fù)雜的問題，具有更強(qiáng)的通用性和可移植能力。

本文的技術(shù)方案步驟包括：采樣獲得等間隔的原始數(shù)據(jù)序列，再根據(jù)當(dāng)前頻率計算新的采樣間隔;在原始數(shù)據(jù)序列的基礎(chǔ)上通過二階插值重構(gòu)新的數(shù)據(jù)序列;對新的采樣序列進(jìn)行FFT變換得到各次諧波的有效值。本文的技術(shù)方法具體包括以下步驟：

（1）錄制采樣序列α，該序列的長度需要根據(jù)實(shí)際周期決定;（2）計算重構(gòu)間隔1/Fs，該間隔需要當(dāng)前的實(shí)時頻率計算得到;（3）重構(gòu)序列β，該步驟需要查找插值位置附近采樣點(diǎn)進(jìn)行二階插值;（4）對新的采樣序列β進(jìn)行FFT變換，得到當(dāng)前頻率F的整數(shù)次方頻率所對應(yīng)的幅值。

2 方案實(shí)施和驗(yàn)證結(jié)果

步驟1：錄制采樣序列α

以額定采樣頻率Fn進(jìn)行采樣，F(xiàn)n的計算公式如（1）所示：

Fn=50*N? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （1）

其中，50（hz）為系統(tǒng)的額定頻率，N為采樣倍率，根據(jù)采樣定理，采樣頻率要大于信號頻率兩倍，采樣之后的數(shù)字信號才可以完整地保留原始信號的信息，所以N應(yīng)該大于需測諧波次數(shù)*2。例如，需測量19次諧波，N應(yīng)該大于38，采樣頻率應(yīng)大于1900。

采樣序列α的錄制長度為1個周期T，周期T需要根據(jù)實(shí)時測頻得到的當(dāng)前頻率F通過公式（2）獲得：T=1/F? ? ? ? ? ? ? ? ? （2）

則采樣序列α錄制的點(diǎn)數(shù)為T/（1/Fn），這里需要注意的是，為了滿足之后的插值運(yùn)算，錄制點(diǎn)數(shù)需要增加兩個點(diǎn)的裕度，即為T/（1/Fn）+2，如圖1。

步驟2：計算重構(gòu)間隔

這個步驟需要一個額外的測頻手段（硬件測頻或軟件測頻）輔助來獲取信號的當(dāng)前頻率F，假設(shè)使用Fs的采樣頻率來重采樣輸入信號，那么新的采樣間隔即為1/Fs。得到重構(gòu)間隔是為下一步獲得重構(gòu)序列β做準(zhǔn)備。Fs的計算公式為：Fs=F*N? ? ? ? ? ? ? ? ?（3）

步驟3：重構(gòu)序列β

這個步驟是最為關(guān)鍵的一步，在已知離散信號序列α和重構(gòu)間隔的基礎(chǔ)上，我們通過插值運(yùn)算獲得重構(gòu)序列β，這個序列等價于系統(tǒng)使用Fs采樣率對原始信號進(jìn)行采樣獲得的離散數(shù)據(jù)。在這里使用了一個二階插值算法，它的精度優(yōu)于線性插值，也具有適度的運(yùn)算量。

二階插值的公式為：

β（P）=α（P0）+α（P1）+ α（P2）? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （4）

其中，P為需要插入的點(diǎn)，即為1/Fs*n（n=1，2…N）。P0、P1、P2為曲線上的已知三個點(diǎn)，公式（3）即為P0、P1、P2所擬合的曲線方程，帶入橫坐標(biāo)P得到這條曲線上的縱坐標(biāo)β（P）。對于P0、P1、P2的選擇采用了內(nèi)插的方式，即P點(diǎn)在P0、P1、P2三個點(diǎn)中間的中間，內(nèi)插方式的精度由于外插方式。

在程序設(shè)計時，這里在序列α中遍歷出小于且最接近于P的點(diǎn)定為P0，P0之后的兩個點(diǎn)分別定為P1、P2。由于Fs=F*N，P2-P1=2/50 *N=0.04N。在系統(tǒng)頻率45hz～55hz的前提下，1/Fs在（0.0181～ 0.0227）N之間，則1/Fs<（P2-P1）恒成立，所以使用這個索引方式可以確保P處于P0～P2以內(nèi)。圖2為內(nèi)插示意圖。

依次計算Pn（n=1/Fs，2/Fs…N/Fs），其中對于（N-1）/Fn和N/Fn，將會使用到序列α中多預(yù)留的第T/（1/Fn）+1和T/（1/Fn）+2兩個點(diǎn)，這是為了滿足二階內(nèi)插的要求。最終將得到新的序列β，它的長度為1/Fs*N即為T，也就是新的序列剛好等分了一個周期的信號。

圖3列舉了一種信號頻率F>50hz的情況下，重構(gòu)序列β的情況。其中黑色為原始序列，紅色為重構(gòu)序列。

步驟4：對新的采樣序列β進(jìn)行FFT變換

新的采樣頻率為Fs，采樣點(diǎn)數(shù)為N，我們通過FFT可以得到頻率n*Fs/N（n=1，2…2/N）即n*F（n=1，2…2/N）下對應(yīng)的幅值，完成了各次諧波的計算。表1顯示了使用Matlab進(jìn)行重采樣前后的數(shù)據(jù)對比。

3 結(jié)語

本文提出的一種提高電能質(zhì)量中諧波測量精度的方法。可以通過重構(gòu)采樣數(shù)據(jù)，達(dá)到無需實(shí)時的設(shè)置采樣率進(jìn)行諧波計算。本技術(shù)方法降低了電力系統(tǒng)頻率的變化導(dǎo)致諧波測量的誤差，提高了算法精度。本文所述方法更加依賴于軟件算法，具有更強(qiáng)的通用性和可移植能力。

參考文獻(xiàn)

[1]陳向群.電力用戶用電信息采集系統(tǒng)[M].中國電力出版社，2012，（12）.

[2]王志軒.中國電力需求側(cè)管理變革[J].新能源經(jīng)貿(mào)觀察，2018，（09）：27-34.

[3]段宇.電網(wǎng)諧波治理和無功補(bǔ)償技術(shù)的研究[D].南昌航空大學(xué)，2012.

[4]張志婕.電力用戶用電信息管理系統(tǒng)的設(shè)計與實(shí)現(xiàn)[D].電子科技大學(xué)，2013.

A Method for Improving the Precision of Harmonic Measurement in Power Grid

XUE Jun

（Jiangsu Linyang Energy Co.， Ltd.， Nanjing Jiangsu? 210019）

Abstract：This paper introduces a method to improve harmonic accuracy in power quality. The steps include： （1） A/D sampling to obtain equally spaced raw data sequences α. （2） Calculate the new sampling interval 1/Fn based on the current frequency. （3） Reconstruct the new data sequence β by second-order interpolation on the basis of the data sequence α. （4） Perform FFT transformation on the new sampling sequence β to obtain the effective value of each harmonic， and obtain relevant data such as harmonic distortion rate and waveform distortion. The method can obtain the synchronization signal without adjusting the sampling frequency， which saves the hardware and software resource overhead and simplifies the system design.

Key words：harmonic measurement; power quality; second-order interpolation reconstruction