王霄翔 程 立 熊慕文

（南京南瑞繼保電氣有限公司，南京 211102）

?

一種非同步采樣下的微機測控裝置直流濾波方法

王霄翔 程 立 熊慕文

（南京南瑞繼保電氣有限公司，南京 211102）

直流分量是影響測量精度的重要因素。使用均方根算法計算有效值時，為了在非同步采樣下濾除直流分量，提出一種在非同步采樣下的濾除直流分量的方法。該方法利用梯形積分原理，在系統(tǒng)頻率和采樣頻率偏離的情況下，將一個周期內(nèi)的采樣點分為整數(shù)點和分?jǐn)?shù)點，分別計算整數(shù)點和分?jǐn)?shù)點的面積來計算直流分量。為求取分?jǐn)?shù)點時對應(yīng)的采樣值，采用了線性插值算法。所提出的算法無需硬件頻率跟蹤電路，數(shù)據(jù)窗長度僅為一個周波，適合多種微機保護測控裝置。實際裝置運行表明，在系統(tǒng)頻率在45～55Hz之間變化時，此方法可以有效濾除直流分量，提高裝置測量精度。

均方根；非同步采樣；直流分量；線性插值

大多數(shù)微機保護算法的計算就是對采樣信號進行頻譜分析而獲得其參數(shù)估算的過程，對算法性能的評價也取決于是否能在較短的數(shù)據(jù)窗中獲得精確的估計值。在各種算法中，傅里葉變換是被廣泛使用進行信號頻譜分析的重要工具。目前在微機繼電保護測控算法中，比較經(jīng)常使用的是離散傅里葉變換[1]。然而實際采樣信號中，往往還包含著一些直流分量，必須要將直流分量可靠的濾除才能得到精確的測量結(jié)果。對于完整周期信號，離散傅里葉變換可以較為容易的濾除直流分量，人們進行了廣泛的研究，如文獻[2]基于三角函數(shù)的和差公式及線性方程組的求解，提出了一種新的遞推離散傅里葉變換算法。在微機測量裝置的實現(xiàn)過程中，如果只需要計算有效值不需要計算諧波分量，離散傅里葉變換算法的計量量較大，實現(xiàn)較為復(fù)雜；使用均方根（RMS）算法較為簡單易行，但均方根算法受直流分量影響較大，需要使用可靠的方法濾除直流分量。測量裝置采樣通道較多，每一路通道的頻率不一定一致，而裝置采樣頻率只能保證和某一通道一致，這樣采樣頻率和通道頻率就不完全同步，這時實際采樣就是為非同步采樣[3]。同步采樣時可以很容易的濾除直流分量，但非同步采樣時，直流分量需要單獨處理。

本文提出一種在非同步采樣情況下濾除直流的算法，主要用于微機測控裝置濾除穩(wěn)態(tài)下的恒定直流分量。該算法基于的梯形積分原理，將系統(tǒng)頻率偏移時，一個周期內(nèi)的采樣點數(shù)分為整數(shù)點和分?jǐn)?shù)點兩部分，并利用分?jǐn)?shù)點的計算結(jié)果對整數(shù)點的計算結(jié)果進行修正。該算法僅需軟件測頻，無需修改硬件電路；算法簡單直觀，僅需一個工頻周期的采樣數(shù)據(jù)，計算量小、延遲短。實際裝置運行表明，在系統(tǒng)頻率發(fā)生變化時，依靠本算法依然能很好的濾除直流分量。

1　算法分析

1.1全波傅氏算法

全波傅氏算法可以濾除恒定的直流分量和基波整數(shù)次諧波分量。

假定被采樣的信號具有如下形式：

式中，a0為直流分量；分別為k次諧波的幅值和初相角。

根據(jù)傅氏級數(shù)原理，可以得到各次諧波分量的實部和虛部的時域表達式為

式中，T為基頻分量的周期；ω為基頻分量的角頻率（ω=2π/T）。

將分解出來的實部和虛部時域表達式離散化處理后可得式中，k為采樣點數(shù)序號；N為每周波采樣點數(shù)。

通過傅氏級數(shù)算法，可以得到直流分量、基波分量以及高次諧波分量。但是如果我們只需要基波分量時，傅氏級數(shù)算法雖然準(zhǔn)確，但計算量大，實現(xiàn)起來比較復(fù)雜，會增加軟件資源消耗。在這種情況下，一般不選用傅氏級數(shù)算法，選擇較為簡單的均方根（RMS）算法，可以很容易的得到有效值。

1.2均方根（RMS）算法

假定被采樣的信號具有式（1）的形式若每周期采樣點數(shù)為 N，可以使用均方根算法得到該信號的有效值，即

式中，x(i)為每個采樣點的值。

從均方根算法的原理可以看出，在混入直流信號時，直流信號也會混入有效值之中，影響計算精度，故必須使用有效的方法濾除直流分量。

同步采樣時，濾除直流分量的方法較多，可以將N個采樣點的算數(shù)平均得到直流分量，但非同步采樣時，所得到的采樣點并不是一個完整周期，該方法已不再適用，需要使用其他方法。

1.3算法原理

設(shè) x(i)是一個周期為 T的正弦連續(xù)信號，將連續(xù)信號離散化，每個周期的采樣點為 N，可以很容易的得到曲線一個周期內(nèi)的面積，即

對于一個不含直流分量的周期信號，正半周波和負(fù)半周波的面積相等，但符號相反，整個曲線所包圍的面積S=0。

如果該周期信號包含直流信號，此時所求得的面積應(yīng)為矩形ABCD的面積，如圖1所示。該信號的直流分量應(yīng)為

求得直流分量之后，可以對每個采樣點進行直流濾波補償，再進行其他計算，這樣就消除了直流分量的影響。

當(dāng)采樣為非同步采樣時，每個周期的采樣點并不是N或者N的整數(shù)倍，這時直接采用式（7）計算面積時就會存在偏差，存在多加了一部分或者少了加了一部分面積的情況，不能得到準(zhǔn)確的直流分量。

圖1　周期信號

當(dāng)系統(tǒng)采樣頻率 f 波動導(dǎo)致基波周期發(fā)生ΔT變化、而采樣頻率 fs固定不變時，一個周期內(nèi)的采樣點數(shù)為

對N取整得到整數(shù)采樣點為

此時分?jǐn)?shù)采樣點為

例如，設(shè)采樣頻率 fs為 4kHz。當(dāng)系統(tǒng)頻率由50Hz偏移到49.6Hz時，可得Ns=80，N=0.65。

在非整周期采樣下，將采樣點N分為整數(shù)采樣點Ns和分?jǐn)?shù)采樣點N，利用梯形積分算法分別計算其面積，然后用分?jǐn)?shù)點面積補償整數(shù)點面積，從而得到了完整的曲線面積，可以求得直流分量。

對于整數(shù)采樣點 Ns，依然可以通過式（7）求得整數(shù)部分面積。

分?jǐn)?shù)采樣點N的面積計算如圖2所示。圖中曲線1位工頻波形，曲線2位頻率偏移時的波形。線段 ANsNsN和曲線 AN所包圍的面積即為分?jǐn)?shù)點面積。為計算此面積，需要首先獲取N處采樣點的值。我們已知Ns和Ns+1處采樣點的值，通過線性插值可以求得N處采樣點的值x(N)。

進而可求得分?jǐn)?shù)采樣點的面積，即

圖2　面積補償算法

利用整數(shù)采樣點部分的面積Ns和分?jǐn)?shù)部分的面積就可以得到此時的直流分量：

至此已經(jīng)得到了直流分量，可以逐個采樣點進行濾除直流的計算。

2　應(yīng)用實現(xiàn)

在變電站中，交流電壓電流經(jīng)過二次PT/CT再經(jīng)接入測控裝置。測控裝置內(nèi)部經(jīng)過A/D轉(zhuǎn)換成數(shù)字信號，然后根據(jù)一定的采樣周期進行采樣，再對采樣點進行濾波后進行保護測量計算。

測控裝置一般有多個采樣通道，比如單間隔線路測控裝置會對母線側(cè)電壓和同期側(cè)電壓同時采樣，需要對每個通道進行濾除直流后再進行計算，我們以母線側(cè)的電壓為基準(zhǔn)進行頻率跟蹤，并以該頻率為基準(zhǔn)進行采樣。在線路開關(guān)合上時，母線側(cè)和線路側(cè)電壓處于同一個系統(tǒng)內(nèi)，他們頻率相同，故測控裝置使用同一個采樣周期對所有通道進行采樣均是同步采樣；當(dāng)線路開關(guān)斷開，需要進行同期合閘時，線路側(cè)和母線側(cè)的頻率就會存在偏差，此時同期側(cè)電壓采樣就是非同步采樣。雖然線路側(cè)和母線側(cè)頻率相差不大，但造成的計算誤差會超過測控裝置要求的0.2%精度范圍，需要單獨處理。流程圖如圖3所示。

測控裝置改進前，仍然使用將N個采樣點算數(shù)平均的方法計算直流分量進行直流濾波。在頻率偏離采樣頻率時會出現(xiàn)較大誤差，且采樣值抖動，不能滿足0.2%精度要求。測量結(jié)果見表1。

采用本文所述方法改進后，在頻率偏離采樣頻率的情況下，裝置計算有效值穩(wěn)定，不再出現(xiàn)抖動或跳動的情況，能夠滿足0.2%的精度要求。測量結(jié)果見表2。

圖3　流程圖

表1　改進前測量結(jié)果

表2　改進后測量結(jié)果

3　結(jié)論

本文提出一種在非同步采樣條件下，利用梯形積分修正方法的濾除直流算法。當(dāng)系統(tǒng)頻率和采樣頻率不一致時，利用線性插值計算分?jǐn)?shù)采樣點，利用分?jǐn)?shù)采樣點的面積和整數(shù)采樣點的面積來進行直流分量的計算。

嵌入式微機保護測控裝置，其計算能力有限，主要用來監(jiān)視正常系統(tǒng)穩(wěn)態(tài)下的運行狀態(tài)。對于系統(tǒng)故障時可能出現(xiàn)的暫態(tài)變化過程不通過測控裝置進行監(jiān)視，故本文所述方法不考慮濾除系統(tǒng)故障時可能出現(xiàn)的衰減直流分量。雖然有很多方法可以較好的濾除恒定的直流分量，但是算法復(fù)雜，運算量較大，不適合實現(xiàn)。此外，我們也可以通過插值的方法保證在不同頻率下每個整周期的采樣點一致，但這樣會消耗較多的資源。本算法既不需要硬件濾波、跟頻，也不需要軟件插值等復(fù)雜運算，在最低資源消耗的情況下，取得了較高的計算精度，非常適合運用于嵌入式微機測控裝置中。目前，本算法以及該成功應(yīng)用于某測控裝置中，并在多個變電站現(xiàn)場長期穩(wěn)定運行。

[1] 陳德樹. 計算機繼電保護原理與技術(shù)[M]. 北京: 水利電力出版社, 1992.

[2] 李晨, 張航, 張愛民, 等. 一種能濾除衰減直流分量的心遞推離散傅氏算法[J]. 繼電器, 2005, 33(17): 17-20.

[3] 邱海峰, 周浩. 非同步采樣下電網(wǎng)諧波分析方法的探討[J]. 繼電器, 2008(1).

[4] 熊崗, 陳陳. 一種能濾除衰減直流分量的交流采樣新算法[J]. 電力系統(tǒng)自動化, 1997, 21(2): 24-26.

[5] 張立華, 徐文立, 常成, 等. 一種適用于微機保護的新的遞推DFT算法[J]. 店里系統(tǒng)自動化, 2000, 24(5): 28-31.

[6] 蘇文輝, 李鋼. 一種能濾去衰減直流分量的改進全波傅氏算法[J]. 電力系統(tǒng)自動化, 2002, 26(23): 42-44.

[7] 黃純, 何怡剛, 江亞群, 等. 交流采樣同步方法的分析與改進[J]. 中國電機工程學(xué)報, 2002, 22(9): 38-42. [8] 蔡菲娜, 左伍衡. 改進的雙速率同步采樣法及其傅里葉變換[J]. 浙江大學(xué)學(xué)報: 工學(xué)版, 2005, 39(3): 414-417.

The Method for Filtering the DC Component of Microcomputer Bay Control Unit with Non-synchronous Sampling

Wang Xiaoxiang Cheng Li Xiong Muwen
(NR Electric Co., Ltd, Nanjing 211102)

DC component is a vital factor which affects the accuracy of measurement. In order to filter out the DC component in the non-synchronous sampling, a method of filtering the DC component is proposed while using the RMS algorithm to calculate the effective value. This method uses the trapezoidal integration principle to divide sampling points in one cycle into integer points and points, and calculate the area of the integer points and the points to calculate the DC component, in the case of system frequency and sampling frequency deviation. Using the linear interpolation algorithm to obtain the sampling value of the fractional point. The proposed method does not need hardware frequency tracking circuit. Data window length is only one fundamental cycle. This method is suitable for a variety of computer protection and monitoring device. Actual device operation shows that this method can effectively remove the DC component and improve the measurement accuracy when the system frequency changes between 45～55Hz.

RMS; non-synchronous sampling; DC component; linear interpolation

王霄翔（1985-），男，工程師，碩士，主要從事變電站自動化系統(tǒng)，智能變電站自動化系統(tǒng)方面的研究。