傅里葉算法暫態(tài)特性以及微機保護的動作時間

劉華龍，阮宜菁

（大連船舶重工集團有限公司，遼寧大連 116000）

0 引言

微機保護算法的關(guān)鍵技術(shù)途徑可表述為：計算體現(xiàn)被保護對象運行方式及特點的物理量參數(shù)，包括電壓、電流等電氣量的有效值和矢量，或算出其基波分量、序分量或某次諧波分量的模值大小及相位大小，并與給定的閾值進行比較[1]。

在微機保護中各個繼電器都是由其相應(yīng)的算法實現(xiàn)的。目前的主流算法都是通過采集輸入的電氣量，經(jīng)算法公式計算得到邏輯原理所需的數(shù)據(jù)，并與定值進行比較。最常用的是傅里葉算法，實現(xiàn)對過流繼電器、濾序器和移相器等的數(shù)據(jù)計算。上述保護算法充分利用微處理器運算速度快、處理速度高的特點，完成傳統(tǒng)保護所不能實現(xiàn)的功能[2]。

算法研究的主要問題有2點：準確度和速度（時間窗）。準確度指的是采用微型計算機采樣所得到的離散數(shù)據(jù)點與原始信號在各個采樣時刻實際信號在數(shù)值上的逼近程度即采集精度。而速度即時間窗表示所用來執(zhí)行保護算法所需要的數(shù)據(jù)時間長度，所用時間越長，則保護動作越慢。在一定采樣頻率下進行離散采樣時，故障發(fā)生后用于執(zhí)行保護算法所需的采樣點個數(shù)被稱為此算法所需的時間窗。

不難發(fā)現(xiàn)，速度和準確度在同一個保護算法中難以兼顧，是矛盾的關(guān)系。因此，研究保護算法的核心就是在準確度和保護速度上進行權(quán)衡。此外，微機保護中通常會引入數(shù)字濾波器，某些保護僅通過算法可實現(xiàn)，而有些保護則需配合數(shù)字濾波器共同工作，數(shù)字濾波器所產(chǎn)生的時間延時也會影響保護的速度。因此評價算法時還需綜合考慮其對數(shù)字濾波器的要求[2]。

1 傅里葉算法

1.1 全周傅里葉算法

式中：N為工頻一周波內(nèi)的采樣點數(shù)。N確定后，對于特定成分的頻率，如對于基波k＝1，式中

1.2 半周傅里葉算法

全波傅里葉算法的有效數(shù)據(jù)窗為1個周波即20 ms，在某些要求保護動作速度較高的故障情況下難以滿足快速性的要求。因此，為了提高保護動作速度，可以采用半波傅里葉算法[3]。半周傅氏算法數(shù)據(jù)窗為12點，響應(yīng)速度很快，但無法濾除直流分量，需與差分濾波器一起使用。

其基本點積計算為

與全波傅氏算法相比，半波傅氏算法不能消除偶次諧波分量（包括直流分量）。因此，在工程應(yīng)用中，為兼顧繼電保護裝置的準確性和快速性，多采取如下策略：故障起始時刻進行半波傅氏算法運算提供保護邏輯所需的計算數(shù)據(jù)，一周波即20 ms后采用全波傅氏算法。

1.3 關(guān)于各類傅里葉算法的評述

綜合上述分析，在響應(yīng)時間方面，半波傅里葉算法由于所使用的數(shù)據(jù)窗口更短，判別故障所需時間就更短，因此保護動作時間更短。但其不能濾除直流分量和偶次諧波分量，對于包含偶次諧波和非周期分量的輸入信號源，采用半波算法所計算的數(shù)據(jù)精度低于全波算法。當追求保護快速性時，應(yīng)當選用半周傅氏算法（如差動保護、短延時等）；當在臨界位置附近時，為提高靈敏度，需要精確計算（如零序過流保護、零序差動保護等），則需選用全周傅里葉算法并使動作速度稍慢一些。

2 過流保護動作時間與返回時間

詳細分析由傅里葉算法影響的過流繼電器動作時間與返回時間。采用傅里葉算法的過流保護暫態(tài)特性如圖1所示。

圖1 短路電流及傅里葉算法的輸出波形

圖1為短路電流，圖2為傅里葉算法輸出。為定量分析傅里葉算法的暫態(tài)特性，圖1中傅里葉算法的輸出可用圖2予以抽象。

圖2 傅里葉算法的過流保護暫態(tài)特性圖

上述過程中，O點為故障發(fā)生時刻，A 點為計算值超過整定值時刻，B點為算法時間窗終止時刻，C點為故障切除時刻，D點為計算值返回時刻，E 點為時間窗終止時刻。A 點、B點、C點、D點對應(yīng)全周傅里葉算法，A′點、B′點、D′點對應(yīng)半周傅里葉算法或半周積分法。

圖2為近似圖形，做出如下假設(shè)：

1）去除非周期分量的震蕩過程，認為時間窗結(jié)束后直接進入穩(wěn)態(tài)。

2）認為過渡過程為線性遞增過程。不同的合閘角度導致不同的遞增過程，但總體而言，時間窗結(jié)束后都能進入穩(wěn)態(tài)。

在上述前提下進行推導，其中tB是時間窗固有常數(shù)：對于全周算法，tB＝20 ms；對于半周算法tB＝10 ms。

tA為自故障發(fā)生后到計算值超過整定值的時間，不妨稱之為上升過渡時間，其大小與故障電流有關(guān)，在穩(wěn)態(tài)故障電流為整定值n倍的情況下）

自計算值超過整定值開始計時至計算值進入穩(wěn)態(tài)的時間，不妨稱之為穩(wěn)態(tài)過渡時間：

容易知道自故障切除后起到計算值小于返回值的時間也是穩(wěn)態(tài)過渡時間，即

在上述前提下，計算保護的可返回時間tF，即自故障后開始計時，在該時間之前由別處斷路器完全切除故障電流，保護裝置不會跳閘。在此概念下，故障起始點為O 點，故障切除點位C 點，故可返回時間tF＝tC。

故障完全切除后，保護的數(shù)據(jù)窗內(nèi)還剩余一部分故障時的數(shù)據(jù)，因此計算值無法立刻變?yōu)?，存在過渡過程；在此期間，保護內(nèi)部的計時器會繼續(xù)計時。因此，欲使保護在此種情況下不動作，必須使得保護計時器起始時刻與保護計時器終止時刻之差值比整定值來得小，即

有時為使動作速度加快，在內(nèi)部計數(shù)時會減去提前量tq，整定值tzd則變?yōu)閠zd－tq。由于計數(shù)時間變短，所以即使動作時間加快了，但可返回時間也相應(yīng)縮短了

3 結(jié)論

通過分析傅里葉算法的暫態(tài)特性以及其與微機保護的動作時間之間關(guān)系可知：

1）在保護的動作延時中加入內(nèi)部提前量不能使總體指標變好，因為雖然加快了動作時間，但可返回時間卻縮短了。

2）使指標變好的唯一方法是縮短時間窗tB，由于抑制非基頻信號的能力是數(shù)據(jù)窗長度的函數(shù)，數(shù)據(jù)窗越長，抑制能力越強。因此存在精度與速度上的矛盾，必須取舍。

3）故障電流越大，動作速度越快?？己藙幼鲿r間時應(yīng)以小電流n＝1.2按標準進行考核，以免保護不能快速動作。故障電流越大，可返回時間越短，故考核可返回時間應(yīng)以大電流n＝10.0來進行，以免保護不能及時返回。在上述2種極端情況下保護仍能滿足時間指標時，可保證過流繼電器滿足使用需求。