基于定頻變點采樣技術(shù)的電網(wǎng)諧波檢測方法

梅 永 王柏林

(1.南京信息工程大學 電子與信息工程學院，南京 210044；2.河海大學 能源與電氣學院，南京 210098)

?

基于定頻變點采樣技術(shù)的電網(wǎng)諧波檢測方法

梅 永1王柏林2

(1.南京信息工程大學 電子與信息工程學院，南京 210044；2.河海大學 能源與電氣學院，南京 210098)

本文提出了一種新的定頻變點電網(wǎng)諧波檢測方法。首先對電網(wǎng)信號進行定頻整周采樣，即以恒定采樣頻率fS對電網(wǎng)電壓(電流)波形采樣，由過零檢測器控制每個數(shù)據(jù)窗口無限逼近整數(shù)(S)個周波，然后對采樣信號進行變點DFT計算。事先把不同點數(shù)的DFT系數(shù)預存在大容量flash存儲器中。本文采用定頻變點法實現(xiàn)同步采樣，同步偏差一定小于1/L0(L0=SfS/55)，只要選取足夠大的采樣頻率fS和采樣周期數(shù)S，采樣同步偏差就足夠小，不存在時延，從而減少了諧波測量同步誤差。

定頻整周采樣；過零檢測；變點DFT； 諧波測量

0 引言

諧波研究的意義，首先是因為諧波的危害十分嚴重。諧波使電能的生產(chǎn)、傳輸和利用的效率降低，使電氣設備過熱、產(chǎn)生振動和噪聲，并使其絕緣老化，使用壽命縮短，甚至發(fā)生故障或燒毀[1]。其次，諧波研究的意義還可以上升到從治理環(huán)境污染、維護綠色環(huán)境的角度來認識。對電力系統(tǒng)這個環(huán)境來說，無諧波就是“綠色”的主要標志之一。因此電力系統(tǒng)諧波電壓測量[2-4]、諧波電流測量[5-6]、隨機環(huán)境下的諧波分析[7-9]、諧波狀態(tài)估計[10]、電能質(zhì)量特征參量的提取[11]具有重要的意義。為此許多學者對電能質(zhì)量問題尤其是諧波測量進行廣泛研究，也作出了重要的貢獻[1-16]。提出了許多諧波測量算法，從早期的模擬濾波器法到傅里葉變換[2-3]、小波變換[14]、希爾伯特黃變換[12-13]、其他算法[14-16]等。

雖然各種新的諧波測量算法層出不窮，由于各種原因大部分算法在工程中的應用還不夠成熟和完善。在電力工程中真正廣為應用的仍然是FFT算法，IEC61000-4-7諧波測量標準推薦的測量算法也仍然是FFT算法。這一事實說明：1)基于FFT的諧波測量算法仍然具有很強的應用價值和研究意義；2)許多新的諧波測量算法還要在工程應用上狠下功夫；3)FFT算法與新算法的融合可能是拉近研究與應用的一條途徑。

1 非同步采樣下的采樣誤差分析

電力系統(tǒng)諧波測量的常用方法是FFT(快速傅立葉變換)法，其基本原理是：首先對電網(wǎng)波形采樣和ADC(模數(shù)轉(zhuǎn)換)得到采樣序列，然后對采樣序列進行FFT得到各次諧波、間諧波的幅值和相位信息。FFT是DFT(離散傅立葉變換)的快速算法，在數(shù)據(jù)點數(shù)很多時，F(xiàn)FT比標準DFT快得多，但FFT要求數(shù)據(jù)點數(shù)必須是2的整次冪。理論上，F(xiàn)FT方法要求同步采樣，即等間隔整周期采樣—— 同一數(shù)據(jù)窗口數(shù)據(jù)的采樣周期相等且每個數(shù)據(jù)窗口的時寬剛好等于整數(shù)個電網(wǎng)周波。由于電網(wǎng)頻率是不斷變化的，同步采樣很難實現(xiàn)，所以，數(shù)據(jù)總存在采樣同步偏差，采樣同步偏差直接造成各次諧波、間諧波的幅值和相位的測量誤差，這稱之為諧波測量同步誤差。

經(jīng)過理論分析和仿真驗證，同步誤差越大造成諧波幅值誤差越大；同步偏差一定的情況下，幅值較小的諧波產(chǎn)生的測量相對誤差越大[17]。

目前減少采樣采樣同步偏差的方法主要有：1)硬件同步：主要指加PLL(鎖相環(huán))的同步方法；2)軟件同步：利用DSP的定時器實時改變采樣頻率。硬件同步和軟件同步的采樣頻率都是變化的，但每一數(shù)據(jù)窗口內(nèi)的采樣點數(shù)相同且都是2的整次冪。PLL同步是IEC61000-4-30標準推薦的方法，IEC61000-4-30標準要求PLL的采樣同步偏差不大于萬分之三。電網(wǎng)頻率恒定不變時PLL做到萬分之三就不容易了，更何況電網(wǎng)頻率時刻都在波動，PLL時刻都處在閉環(huán)調(diào)節(jié)的動態(tài)過程之中，這種閉環(huán)調(diào)節(jié)還存在時延，所以實際應用中PLL很難保證采樣同步偏差不大于萬分之三，傳統(tǒng)的硬件同步和軟件同步在實際應用中都不容易獲得高精度。IEC61000-4-30(2003)標準還要求FFT加矩形窗，Matlab仿真表明：同步偏差為萬分之三時，矩形窗FFT對40次以上的諧波、間諧波的計算結(jié)果就是超標的。軟件同步可以看作是一個軟件PLL，除了上述弊端，軟件同步比硬件PLL的時延大，因為后者是鎖相、前者是鎖周。鎖周是用前一個周波的寬度確定下一個周波的采樣周期，時延至少是一個周波。軟件同步還有個問題就是定時器的分數(shù)定時誤差。

本文提出一種固定采樣頻率、采樣點數(shù)可變的采樣技術(shù)，能很好的減少同步偏差造成的諧波測量誤差。

2 定頻變點采樣方法

所謂定頻變點采樣技術(shù)是以恒定采樣頻率fS對電網(wǎng)電壓(電流)波形采樣，由過零檢測器控制每個數(shù)據(jù)窗口無限逼近S個整周波；用固定的采樣頻率fS對電網(wǎng)波形采樣，允許采樣頻率很高，而且越高越有利。

例如一個數(shù)據(jù)窗口的時間寬度規(guī)定為S個電網(wǎng)整周波，S個電網(wǎng)周波的開始(電網(wǎng)波形過零)采當前數(shù)據(jù)窗口的第1點，S個電網(wǎng)周波的結(jié)束(電網(wǎng)波形過零)停止本數(shù)據(jù)窗口的采樣，采樣同時計算采樣點數(shù)Lk。按IEC61000-4-30標準的要求采樣周波數(shù)S=10。

設Lk的可能取值是：

L0,L0+1,…,L0+n,…,L0+N-1，其中L0、n、N為正整數(shù)，N個Lk分別對應第1,2,…,n,…,N組DFT系數(shù)，根據(jù)三角函數(shù)的循環(huán)性，第n組DFT系數(shù)只需存儲L0+n個三角函數(shù)值，不需要存2M(L0+n)個三角函數(shù)值。

考慮電網(wǎng)頻率f波動的最大范圍為45～55Hz， 如果采樣頻率fS和數(shù)據(jù)窗口的周波數(shù)S都固定不變，那么f=55Hz時的采樣點數(shù)L0最小：

f=45Hz時的采樣點數(shù)L0+N-1最大：

其中，round{ }表示取整運算。

3 變點數(shù)DFT計算諧波參數(shù)

設當前數(shù)據(jù)窗口的采樣點數(shù)為Lk，且Lk個采樣數(shù)據(jù)(采樣序列)為：

y=[y(1)y(2) …y(Lk)]

當前數(shù)據(jù)窗口的采樣結(jié)束后立即計算Lk點實序列三角函數(shù)型DFT：

X=FyT

式中，yT表示向量y的轉(zhuǎn)置；F為實序列三角函數(shù)型DFT的系數(shù)矩陣：

其中，M是要求的頻點數(shù)。按IEC61000-4-30標準的要求：需要測量諧波、間諧波共500個頻點即M=500。一次DFT需完成500×Lk次乘加運算。

以前由于計算機的數(shù)據(jù)處理能力不高，DFT曾經(jīng)是難以實施的，如今高速DSP(數(shù)據(jù)處理器)的出現(xiàn)已經(jīng)使得直接用DFT作為諧波分析的工具成為現(xiàn)實。DFT的最大優(yōu)勢在于不要求一個數(shù)據(jù)窗口的采樣點數(shù)是2的整次冪，我們則進一步發(fā)展為：不要求每個數(shù)據(jù)窗口的采樣點數(shù)相同即變點DFT。

為了提高DFT的計算速度，我們采取了兩個措施：

1)不是進行Lk×Lk次乘加運算，而是進行500×Lk次乘加運算，因為高于50次的諧波、間諧波不需要測量。

2)不同點數(shù)的DFT系數(shù)(系數(shù)矩陣F的元素)用查表法得到。把不同點數(shù)的DFT系數(shù)預存在大容量flash存儲器中，設Lk的可能取值是L0,L0+1,…,L0+n,…,L0+N-1，它們分別對應第1,2,…,n,…,N組DFT系數(shù)，根據(jù)三角函數(shù)的循環(huán)性，第n組DFT系數(shù)只需存儲L0+n個三角函數(shù)值，不需要存2M(L0+n)個三角函數(shù)值。

采取了這兩個措施的DFT比FFT慢不了多少，用現(xiàn)代DSP足以從容地完成。

4 仿真研究

給出仿真模型如下：

取DSP定時器工作頻率為150MHz，取采樣率fS=60kHz、采樣周期S=10，計算出采樣點數(shù)L0=10909，采樣同步偏差小于萬分之一。表1是不同電網(wǎng)頻率下，本算法所得到的諧波幅值的相對誤差(Matlab仿真結(jié)果)：

表1中，E1、E2、E3、E25、E31、E48、E49、E50分別代表第1、2、3、25、31、48、49、50次諧波幅值的相對誤差。可見，本方法很容易使諧波精度達標。

5 結(jié)論

本文不用傳統(tǒng)的硬件同步和軟件同步，而是采用定頻變點法實現(xiàn)同步。即固定采樣頻率，不需要每個數(shù)據(jù)窗口的采樣點數(shù)相同， 解決了各個數(shù)據(jù)窗口的數(shù)據(jù)點數(shù)不同且不是2的整次冪的問題，這種方案即使在電網(wǎng)頻率不斷波動的情況下，也能實現(xiàn)同步偏差小于萬分之一、甚至更小，不存在時延。減少了采樣同步偏差就直接減少了諧波測量誤差。對于不同點數(shù)的DFT的系數(shù)利用可編程片上系統(tǒng)計算出來存在數(shù)據(jù)存儲器里，供DSP實時計算諧波測量值查表，這樣大大節(jié)省了運算時間，提高了諧波分析的實時性。本文提出的減少采樣同步偏差的方法有如下特點：

表1 諧波幅值的相對誤差matlab仿真結(jié)果表

1)結(jié)構(gòu)簡單、易實現(xiàn)，定頻整周采樣器可以用過零檢測器加ADC構(gòu)成，變點DFT計算器可以用高速DSP(或PSOC)加flash存儲器構(gòu)成，不需要常見的PLL同步環(huán)節(jié)；

2)很容易做到采樣同步偏差小于萬分之一、甚至更小，而傳統(tǒng)的硬件同步和軟件同步方法很難做到采樣同步偏差小于萬分之一；

3)電網(wǎng)頻率波動與電網(wǎng)頻率相對穩(wěn)定時的采樣同步偏差基本一致，而傳統(tǒng)的硬件同步和軟件同步在電網(wǎng)頻率波動時的采樣同步偏差比電網(wǎng)頻率相對穩(wěn)定時大得多；

4)同步不存在時延，而硬件同步和軟件同步同步都存在時延；

5)極小的采樣同步偏差保證了極小的諧波測量誤差。

[1] 林海雪.公用電網(wǎng)諧波國標中的幾個問題[J].電網(wǎng)技術(shù)，2003，27(1)

[2] 劉敏，王克英. 基于快速傅里葉變換與誤差最小原理的諧波分析方法[J].電網(wǎng)技術(shù)，2006，30(19)

[3] 薛蕙, 楊仁剛. 基于FFT 的高精度諧波檢測算法[J].中國電機工程學報, 2002, 22(12)

[4] 李明，王曉茹.采用時頻原子方法的電壓閃變檢測[J].電網(wǎng)技術(shù)，2010，34(7)

[5] 戴朝波，林海雪. 兩種諧波電流檢測方法的比較研究[J].中國電機工程學報，2002，22(1)

[6] 劉國海，呂漢聞，劉穎，等. 基于改進RLS 算法的諧波電流檢測方法[J]. 電力自動化設備，2010，30(10)

[7] L. Sainz, M. Caro, J. Pedra, Study of electric system harmonic response,IEEE Trans. Power Deliv. 19 (2) (2004) 868-874.

[8] L. Sainz, J. Pedra, M. Caro, Steinmetz circuit influence on the electric system harmonic response, IEEE Trans. Power Deliv. 20 (2) (2005)

[9] 劉應梅，高玉潔.基于Prony 法的暫態(tài)擾動信號分析[J].電網(wǎng)技術(shù)，2006，30(4)

[10] 梁志瑞, 葉慧強, 趙飛. 電力系統(tǒng)諧波狀態(tài)估計研究綜述[J].電力系統(tǒng)保護與控制, 2010, 38(15)

[11] 賈清泉，于連富，董海艷，等.應用原子分解的電能質(zhì)量擾動信號特征提取方法[J].電力系統(tǒng)自動化，2009，33(24)

[12] 張宇輝，賀健偉，李天云，等.基于數(shù)學形態(tài)學和HHT 的諧波和間諧波檢測方法[J].電網(wǎng)技術(shù)，2008，32(17)

[13] 蘇玉香，劉志剛，李科亮，等.Hilbert-Huang 變換在電氣化鐵路諧波檢測中的應用[J].電網(wǎng)技術(shù)，2008，32(18)

[14] 汪斌，王年，蔣云志，等.改進FastICA 算法在諧波檢測中的應用[J].電力自動化設備，2011，31(3)

[15] 郭成，李群湛，賀建閩，等.電網(wǎng)諧波與間諧波檢測的分段Prony算法[J].電網(wǎng)技術(shù)，2010，34(3)

[16] 趙聞蕾，鄒積巖，孔莉，等.基于小波變換的電氣化鐵道牽引負荷諧波分析[J].鐵道學報，2011，33(6)

[17] 梅永，王柏林.電力系統(tǒng)諧波分析的同步校正法.電力自動化設備， 31(1)

10.3969/j.issn.1000-0771.2015.09.03