韓 晉，黃 健

（江蘇金智科技股份有限公司，江蘇 南京 211100）

0 引言

同期并列的理想條件是并列斷路器兩側(cè)電源電壓的3個(gè)狀態(tài)量全部相等，即：待并側(cè)頻率與系統(tǒng)側(cè)頻率相等；待并側(cè)電壓與系統(tǒng)側(cè)電壓幅值相等；相角差為零。壓差和頻差的存在將導(dǎo)致并網(wǎng)瞬間，并列點(diǎn)兩側(cè)出現(xiàn)一定的無(wú)功功率和有功功率的交換。電網(wǎng)和發(fā)電設(shè)備一般都具有承受一定功率交換的能力。相對(duì)而言，相角差的存在會(huì)給斷路器兩側(cè)帶來(lái)更多的傷害，嚴(yán)重時(shí)會(huì)誘發(fā)次同步諧振。在實(shí)際并列操作中，只要并列合閘時(shí)沖擊電流較小，不危及電氣設(shè)備，壓差和頻差允許有一定的偏差，相角差要盡可能接近零。因此，應(yīng)確保同期捕捉算法在相角差為零時(shí)完成并網(wǎng)。

機(jī)組同期在壓差和頻差不滿足要求時(shí)，要對(duì)電壓、頻率進(jìn)行調(diào)節(jié)，一般采用PID調(diào)節(jié)算法或模糊控制算法。一旦壓差和頻差條件滿足，就要采用同期捕捉算法預(yù)測(cè)最佳合閘時(shí)機(jī)。本文主要研究同期捕捉算法，及對(duì)算法有直接影響的相角差的量測(cè)系統(tǒng)。

1 常用算法

并列點(diǎn)兩側(cè)的電壓矢量差Ud=Ug-Us的波形如圖1所示，其為周期的脈沖電壓，點(diǎn)B是相角差為0（即Ud為0）的時(shí)刻?？紤]到合閘回路動(dòng)作有一定延時(shí)，為了能在點(diǎn)B合閘，要提前發(fā)出合閘命令，如在點(diǎn)A發(fā)出合閘命令（tDC為合閘回路動(dòng)作時(shí)間）。同期捕捉算法的任務(wù)是選擇點(diǎn)A。Ud的每個(gè)周期有且只有1個(gè)最佳的合閘點(diǎn)，實(shí)際運(yùn)行中要求能夠捕捉到第1次出現(xiàn)的最佳合閘點(diǎn)，以實(shí)現(xiàn)準(zhǔn)確快速的并網(wǎng)。

常用的同期捕捉算法有2種，算法1是利用相角差與滑差及滑差加速度的關(guān)系的預(yù)測(cè)方法，算法2是基于相角差序列｛δi｝、運(yùn)用現(xiàn)代預(yù)測(cè)理論的預(yù)測(cè)方法。

圖1 滑差波形示意圖Fig.1 Waveform of slip frequency

1.1 算法1

考慮相角差可能存在加速度的情況，相角差隨滑差及滑差加速度變化的模型為：

其中，δ為相角差，ω為滑差角頻率，t為時(shí)間。

軟件根據(jù)合閘回路動(dòng)作時(shí)間tDC及式（1）不斷計(jì)算導(dǎo)前角 δYJ，并檢測(cè)實(shí)時(shí)相角差 δ，捕捉到 δ-δYJ＜ ε時(shí)，發(fā)出合閘信號(hào)。只有在ω及Δω/Δt始終保持為常數(shù)的前提下，才能保證該模型的預(yù)測(cè)精度。

1.2 算法2

從相角差量測(cè)系統(tǒng)獲得的相角差序列｛δi｝（i為時(shí)間遞增方向）雖然是一維序列，但其中包含了｛δi｝變化的各次函數(shù)的信息，如滑差 ωs、滑差加速度 ω˙s、加速度變化率¨s等，以及可能出現(xiàn)的量測(cè)噪聲。算法2應(yīng)用現(xiàn)代預(yù)測(cè)理論對(duì)上述信息進(jìn)行挖掘來(lái)預(yù)測(cè)合閘點(diǎn)，本文主要介紹2種預(yù)測(cè)算法：一種是對(duì)該序列采用預(yù)測(cè)的最簡(jiǎn)模型來(lái)預(yù)測(cè)同期點(diǎn)，同時(shí)對(duì)采樣數(shù)據(jù)進(jìn)行殘差遞推的分段辨識(shí)以提高預(yù)測(cè)精度；另一種是對(duì)序列采用多元線性回歸，進(jìn)行殘差分析，對(duì)預(yù)測(cè)值進(jìn)行修正。2種預(yù)測(cè)算法都是根據(jù)已有的相角差序列不斷預(yù)測(cè)tDC之后的相角差，捕捉到＜ε時(shí)，發(fā)出合閘信號(hào)。

算法1、2都是運(yùn)行于離散的數(shù)字采樣系統(tǒng)中，為了避免丟失合閘時(shí)機(jī)，ε不能過(guò)小，然而也不能過(guò)大，否則會(huì)降低同期質(zhì)量。算法2與算法1相比能避免受到太多量測(cè)噪聲和系統(tǒng)波動(dòng)的影響，早期的微機(jī)準(zhǔn)同期裝置因微機(jī)水平比較落后，采樣誤差比較大，多采用此模型。隨著目前微機(jī)技術(shù)水平的提高，采樣精度已不再是問(wèn)題。這2種算法的應(yīng)用都很成熟。

筆者所開(kāi)發(fā)的同期裝置采用了雙CPU的結(jié)構(gòu)，分別運(yùn)行以上2種算法。2個(gè)CPU的判讀結(jié)果構(gòu)成表決系統(tǒng)，只有在2種算法同時(shí)判斷為可以合閘時(shí)，才能發(fā)出合閘命令。

2 問(wèn)題的提出

由于相角差的變化是連續(xù)的，而CPU對(duì)相角差的測(cè)量和計(jì)算、算法的計(jì)算都是離散的。因此，理論上無(wú)法每次都能捕捉到相角差絕對(duì)為零的時(shí)刻。在系統(tǒng)頻率為工頻的情況下，合閘時(shí)刻的最大理論誤差為±10 ms，該誤差隨系統(tǒng)頻率的不同而不同。對(duì)于高速斷路器而言，精度偏低。為了獲得更好的同期時(shí)刻，應(yīng)插入一個(gè)更短的時(shí)間間隔，判別在本次發(fā)出合閘脈沖是否比下一次合適。

以下舉例分析離散間隔對(duì)合閘精度的影響。根據(jù)相角差與頻差的近似關(guān)系：δ=ωst=2πΔf t（ωs為兩側(cè)滑差；Δf為兩側(cè)頻率差），以Δf=0.25 Hz為例進(jìn)行計(jì)算。

由計(jì)算結(jié)果可知，以10 ms為間隔，算法本身給同期效果帶來(lái)的誤差范圍為［-0.9°，+0.9°］，而以1 ms為間隔，算法本身給同期效果帶來(lái)的誤差范圍將減小到［-0.09°，+0.09°］，接近于 0。

由以上分析可知：在相同的條件下，減小數(shù)據(jù)離散度可顯著提高合閘精度。而只有進(jìn)一步提高原始數(shù)據(jù)密度，并對(duì)算法進(jìn)行相應(yīng)的改進(jìn)，才能減小合閘離散度。以上分析是在測(cè)量無(wú)誤差、未計(jì)及滑差加速度、導(dǎo)前時(shí)間設(shè)置無(wú)誤差的前提下進(jìn)行的，以下將針對(duì)此問(wèn)題討論同期捕捉算法的改進(jìn)和實(shí)現(xiàn)。

3 算法改進(jìn)

合閘時(shí)刻是以相角差序列、滑差序列和導(dǎo)前時(shí)間為基礎(chǔ)通過(guò)算法預(yù)測(cè)的；同期捕捉算法只有在數(shù)據(jù)更新后計(jì)算才有意義；而相角差量測(cè)序列｛δi｝的更新是以10 ms為周期的。2種算法都存在相角差時(shí)間離散過(guò)大而引起誤差過(guò)大的問(wèn)題。為了將合閘動(dòng)作時(shí)間的誤差減小至±1 ms，必須每隔1 ms執(zhí)行一次捕捉算法，并且相角差要在1 ms內(nèi)得到更新。

3.1 算法1的改進(jìn)及實(shí)現(xiàn)

算法1的執(zhí)行分為根據(jù)模型求解導(dǎo)前角和根據(jù)相角差的實(shí)測(cè)值與導(dǎo)前角比較判斷合閘時(shí)機(jī)兩部分。

在已知合閘回路動(dòng)作時(shí)間的前提下，根據(jù)式（1）不斷計(jì)算合閘導(dǎo)前角δYJ（即在此時(shí)電壓條件下，為使斷路器在相角差接近零時(shí)并網(wǎng)而提前的相角差角度）。式（1）所需滑差及滑差的加速度計(jì)算如下：

式（2）和式（3）所需要數(shù)據(jù)從序列｛δi｝中取得。

由于｛δi｝是10 ms更新一次，原有算法在相角差更新時(shí)，計(jì)算一次導(dǎo)前角，并判斷一次合閘條件。本文采用插值原理對(duì)較短時(shí)間內(nèi)的相角差進(jìn)行估計(jì)，即隨后的10 ms內(nèi)每隔1 ms進(jìn)行一次估計(jì)，估計(jì)值為每 1 ms更新一次，因此可每隔 1 ms判斷是否滿足若條件滿足，則發(fā)合閘命令，否則繼續(xù)執(zhí)行算法。一旦有新的相角差產(chǎn)生，即序列｛δi｝更新后，將采用新的數(shù)據(jù)計(jì)算和 δYJ。

3.2 算法2的改進(jìn)

所有的預(yù)測(cè)算法都是以序列｛δi｝為基礎(chǔ)進(jìn)行預(yù)測(cè)的。預(yù)測(cè)算法有很多種，以下介紹2種：預(yù)測(cè)算法1［7］對(duì)序列進(jìn)行線性擬合；本文提出了預(yù)測(cè)算法2，對(duì)數(shù)據(jù)進(jìn)行拋物線擬合。2種算法都挖掘了預(yù)測(cè)算式中的有用信息，對(duì)可能出現(xiàn)的誤差進(jìn)行了修正。

3.2.1 預(yù)測(cè)算法 1［7］

對(duì)于相角差序列｛δi｝（i為時(shí)間遞增方向），其一步預(yù)測(cè)模式為：

實(shí)際中采用預(yù)測(cè)的最簡(jiǎn)模型，根據(jù)遞推殘差辨識(shí)預(yù)測(cè)原理，對(duì)預(yù)測(cè)值與實(shí)際值 δi的插值再進(jìn)行遞推高階參數(shù)，直至

3.2.2 預(yù)測(cè)算法2

對(duì)相角差序列進(jìn)行拋物線擬合，擬合基函數(shù)y=at2+bt+c，假設(shè) t=0 時(shí)，y=δi-2；t=T 時(shí)，y=δi-1；t=2T時(shí)，y=δi；對(duì) t=3T 時(shí)的時(shí)的分別進(jìn)行一步預(yù)測(cè)和二步預(yù)測(cè)。通過(guò)求解方程組

加入殘差修正后的一步預(yù)測(cè)模型和二步預(yù)測(cè)模型分別為：

以上2種預(yù)測(cè)算法對(duì)原始序列要求量小，運(yùn)算簡(jiǎn)單。實(shí)際滑差的變化為非線性的，因此預(yù)測(cè)算法2比預(yù)測(cè)算法1更能真實(shí)反映實(shí)際情況，預(yù)測(cè)精度更高。

假設(shè)在勻速滑差的情況下，δk為實(shí)際值，為估計(jì)值，εk為該點(diǎn)的估計(jì)誤差，εk的方差為σ2。第n步預(yù)測(cè)為：

不難求出其預(yù)測(cè)方差為2n2σ2，可見(jiàn)預(yù)測(cè)方差是隨預(yù)測(cè)步數(shù)的平方增加的。

因此，當(dāng)tDC一定時(shí)，算法的精度就取決于步長(zhǎng)，步長(zhǎng)越短，步數(shù)就越多，預(yù)測(cè)的準(zhǔn)確度就急劇下降。例如tDC為100 ms時(shí)，對(duì)其進(jìn)行100次步長(zhǎng)為1 ms的預(yù)測(cè)，其誤差遠(yuǎn)大于對(duì)其進(jìn)行一次步長(zhǎng)為100 ms的預(yù)測(cè)。為了提高預(yù)測(cè)的準(zhǔn)確度，采用加大預(yù)測(cè)步長(zhǎng)的方法來(lái)減少預(yù)測(cè)步數(shù)所造成的誤差。因此本文所提出的預(yù)測(cè)算法都是針對(duì)固定間隔是tDC的序列。

4 改進(jìn)算法的實(shí)現(xiàn)

以下所討論的相角差的量測(cè)系統(tǒng)適用于算法1和算法2。算法1的實(shí)現(xiàn)比較簡(jiǎn)單，不再贅述，以下主要討論算法2的實(shí)現(xiàn)。

4.1 相角差的獲取

同期捕捉算法中最主要的參量就是相角差。相角差量測(cè)可以通過(guò)測(cè)量電壓瞬時(shí)值再進(jìn)行數(shù)學(xué)處理而獲得，也可以通過(guò)2個(gè)同相電壓的正向過(guò)零點(diǎn)的時(shí)差來(lái)計(jì)算相角差而獲得。

第1種方法采用富士算法計(jì)算兩側(cè)電壓的瞬時(shí)相量而獲得相角差。由于并列雙方頻率不相等，難以同時(shí)對(duì)其進(jìn)行整周期采樣，計(jì)算結(jié)果難免會(huì)存在不同步誤差。為獲取準(zhǔn)確的合成相量，必須對(duì)采樣數(shù)據(jù)進(jìn)行同步化處理。文獻(xiàn)［5］提出利用雙線性插值平均法對(duì)采樣數(shù)據(jù)進(jìn)行同步化，然而這只適用于頻差較小的情況，對(duì)于頻差較大的情況，難以確定2路采樣數(shù)據(jù)的時(shí)間對(duì)應(yīng)關(guān)系。因此該方法不具有可行性。

第2種方法將正弦波轉(zhuǎn)變?yōu)榉讲ǎ鶕?jù)方波的過(guò)零比較計(jì)算來(lái)獲得兩側(cè)的頻率和相角差。該方法可獲得充足和高質(zhì)量的原始數(shù)據(jù)，但數(shù)據(jù)離散度較大。由于發(fā)電機(jī)和系統(tǒng)的慣性時(shí)間常數(shù)很大，可認(rèn)為在幾十毫秒內(nèi)發(fā)電機(jī)和系統(tǒng)的頻率不變，即在此段時(shí)間內(nèi)同期點(diǎn)兩側(cè)頻差是勻速運(yùn)動(dòng)。因此本文提出采用線性插值和線性外推值的方法來(lái)增加數(shù)據(jù)的密度。

線性插值原理見(jiàn)圖2，只用于為算法2生成序列｛δ′j｝。線性外推的原理見(jiàn)圖 3，在 2種算法中都有應(yīng)用，用于1 ms中斷中更新相角差。

圖2 線性插值原理圖Fig.2 Schematic diagram of linear interpolation

圖3 線性外推原理圖Fig.3 Schematic diagram of linear extrapolation

算法 2 需要｛δi｝、｛δ′j｝、｛δ″k｝3 個(gè)序列。固定間隔為10 ms的采樣序列｛δi｝由采樣計(jì)算獲得，經(jīng)過(guò)插值計(jì)算后轉(zhuǎn)變?yōu)楣潭ㄩg隔為 1 ms 的序列｛δ′j｝，從中提取固定間隔為 tDC的序列｛δ″k｝，｛δ″k｝即為算法所需序列。在1 ms中斷中，｛δi｝不會(huì)更新，序列｛δ′j｝和｛δ″k｝將會(huì)更新。算法根據(jù)序列｛δ″k｝不斷預(yù)測(cè) δ″k+1或 δ″k+2。

序列｛δi｝的時(shí)間間隔隨頻率的不同而不同，實(shí)際運(yùn)行中電壓頻率可能不是50 Hz，該序列間隔將不是10 ms，因此要根據(jù)實(shí)際的頻率進(jìn)行相應(yīng)的轉(zhuǎn)變，保證｛δ″k｝為固定間隔是tDC的序列，否則會(huì)帶來(lái)預(yù)測(cè)誤差。

本文提出的預(yù)測(cè)模型要求被測(cè)序列為遞增序列，而在一個(gè)周期內(nèi)相角差變化先遞增后遞減，因此要做適當(dāng)處理，將相角差序列轉(zhuǎn)變?yōu)閱握{(diào)遞增序列。

4.2 算法實(shí)現(xiàn)的若干問(wèn)題

算法2的實(shí)現(xiàn)主要包括相角差序列的獲取、預(yù)測(cè)隊(duì)列的創(chuàng)建、算法的計(jì)算、殘差分析等。以下主要討論判別預(yù)測(cè)點(diǎn)的單調(diào)性和針對(duì)tDC過(guò)大而需要進(jìn)行二步預(yù)測(cè)的問(wèn)題。

算法2的實(shí)現(xiàn)采用兩段預(yù)測(cè)的方式。第1段是對(duì)tDC-1 ms之后的相角差進(jìn)行預(yù)測(cè)，產(chǎn)生預(yù)測(cè)值1；第2段是對(duì)tDC之后的相角差進(jìn)行預(yù)測(cè)，產(chǎn)生預(yù)測(cè)值2。因此需要 2 個(gè)不同時(shí)間間隔的預(yù)測(cè)序列｛δ″k｝，這2個(gè)序列均由序列｛δ′j｝轉(zhuǎn)變而來(lái)。2段預(yù)測(cè)所得的2個(gè)結(jié)果用來(lái)判斷tDC之后的相角差的單調(diào)性，以確保不會(huì)逆功率合閘。只有2段預(yù)測(cè)值都遞增，且小于一定值，才認(rèn)為合閘條件滿足。其條件表達(dá)式如下：

其中，ε1和ε2為允許誤差。如圖4所示，僅滿足條件式（6）將在區(qū)域A內(nèi)合閘，條件式（7）是為避免在區(qū)域B內(nèi)合閘，只有2個(gè)條件都滿足才能在區(qū)域C合閘。

圖4 合閘區(qū)域示意圖Fig.4 Schematic diagram of closing area

大多數(shù)情況下，一個(gè)周期內(nèi)的數(shù)據(jù)是足夠使用的，這與合閘回路的動(dòng)作時(shí)間與兩側(cè)頻差有關(guān)。目前廣泛使用的SF6斷路器和真空斷路器的動(dòng)作時(shí)間一般不超過(guò)100 ms，加上中間繼電器的動(dòng)作時(shí)間等，一般不會(huì)超過(guò)200 ms。即使以最短的周期，允許頻差0.5 Hz為例，一個(gè)周期將是2 s，也有足夠的數(shù)據(jù)支持算法預(yù)測(cè)同期點(diǎn)。因此針對(duì)合閘回路動(dòng)作時(shí)間不太大的情況，取預(yù)測(cè)步長(zhǎng)為tDC，根據(jù)式（4）進(jìn)行一步預(yù)測(cè)?？紤]到極個(gè)別的對(duì)于合閘回路動(dòng)作時(shí)間比較大的情況，為了能在一個(gè)周期內(nèi)得到足夠多的數(shù)據(jù)，取預(yù)測(cè)步長(zhǎng)為tDC/2，根據(jù)式（5）進(jìn)行二步預(yù)測(cè)。

4.3 算法的實(shí)驗(yàn)驗(yàn)證

圖5 同期合閘錄波圖Fig.5 Recorded wave of synchronous closing

實(shí)驗(yàn)證明減小數(shù)據(jù)離散度后的算法也能明顯改善合閘質(zhì)量，鑒于之前的分析，不再給出錄波波形。

實(shí)驗(yàn)證明改進(jìn)的算法確實(shí)提高了合閘的精度。改進(jìn)后的算法計(jì)算量有所增加，然而現(xiàn)在CPU處理能力強(qiáng)大，計(jì)算量已不再是問(wèn)題。

5 結(jié)語(yǔ)

同期捕捉算法有2種類型：一種是對(duì)相角差隨滑差變化規(guī)律的模型的求解；一種是對(duì)相角差序列進(jìn)行分析，對(duì)合閘點(diǎn)進(jìn)行預(yù)測(cè)。預(yù)測(cè)算法有多種，本文介紹了分段遞推殘差和基于殘差分析的多元線性回歸的預(yù)測(cè)方法。針對(duì)2類算法中同時(shí)存在離散度過(guò)大而造成的誤差，本文對(duì)算法進(jìn)行了改進(jìn)，并對(duì)對(duì)算法有直接影響的相角差量測(cè)系統(tǒng)進(jìn)行了處理，即通過(guò)方波過(guò)零比較獲取最原始的數(shù)據(jù)，通過(guò)插值計(jì)算估計(jì)值來(lái)增加數(shù)據(jù)密度。實(shí)驗(yàn)證明改進(jìn)后的算法使同期合閘質(zhì)量有了明顯提高。所介紹的算法不僅適用于發(fā)電機(jī)同期和線路同期，也適用于對(duì)重合閘等情況。