基于最小二乘曲線擬合提取基頻分量

牛智銳，宋運忠

（河南理工大學(xué) 電氣工程與自動化學(xué)院，焦作 454003)

0 引言

輸電線路監(jiān)測到的故障中80%是單相接地故障，而瞬時故障占其中的50%～90%[1]。瞬時故障是指在線路上僅持續(xù)幾個周期的暫時性故障。雖然瞬時故障能夠在一定延時后重合成功，產(chǎn)生的故障電弧熄滅后不會對電力系統(tǒng)和電氣設(shè)備造成重大損壞，但是這些故障通常是同一位置發(fā)生永久性故障的前兆[2，3]。瞬時故障造成電流巨大變化，電路中的電感會阻礙電流的變化，產(chǎn)生隨時間呈指數(shù)衰減的直流偏移分量。直流偏移分量對基頻相量的提取造成影響，導(dǎo)致故障定位產(chǎn)生較大的誤差。

如何快速消除直流偏移分量并準(zhǔn)確地提取基頻分量，一直是國內(nèi)外學(xué)者研究的熱點問題。采用全波和半波傅里葉算法提取故障信號的基頻分量，存在計算時間長、響應(yīng)速度慢、無法完全濾除直流偏移分量的缺點[4～6]。使用同步相量估計算法[7，8]修正離散傅里葉變換（Discrete Fourier Transform，DFT）可以較好的消除直流偏移分量提取基頻分量，缺點是計算量較大。文獻(xiàn)[9]將Hilbert變換與DFT相結(jié)合計算基頻相量，該算法的計算過程非常復(fù)雜。文獻(xiàn)[10]先對輸入的電流信號進(jìn)行積分，然后在每個采樣瞬間減去衰減的直流偏移分量獲得基頻相量。該算法將故障持續(xù)一個周期后的采樣數(shù)據(jù)用于計算，其難點在于必須能夠準(zhǔn)確地檢測出故障的起始時間。文獻(xiàn)[11]需要2.25個周期的波形數(shù)據(jù)，將經(jīng)驗?zāi)B(tài)分解與Hilbert變換相結(jié)合計算基頻相量。算法易受波形周期的限制，不適合應(yīng)用于繼電保護(hù)。使用最小二乘算法[12]消除直流偏移分量提取基頻相量，算法的性能會隨著時間常數(shù)的減小而降低，難以滿足電力系統(tǒng)對于穩(wěn)定性的要求。

針對上述方法存在數(shù)學(xué)模型較為復(fù)雜、運算量較大的缺點，本文構(gòu)建新的數(shù)學(xué)模型提取基頻分量。利用均方根域-小波變換方法獲取完整的故障時間段數(shù)據(jù)，使用非線性最小二乘曲線擬合提取基頻分量。通過多個算例的仿真驗證和分析，算法能夠快速準(zhǔn)確地提取基頻分量。將其應(yīng)用于繼電保護(hù)裝置開展故障定位研究工作，可以很好的與行波法[13～16]相互補(bǔ)充，更好地保障輸電系統(tǒng)的安全、可靠運行。

1 構(gòu)建故障信號數(shù)學(xué)模型

當(dāng)輸電系統(tǒng)發(fā)生單相接地故障時，電流產(chǎn)生衰減的直流偏移分量。將故障電流信號看作由一個含有諧波的交流分量iAC和一個隨時間常數(shù)呈指數(shù)衰減的直流偏移分量iDC組成[17]：

式中，直流偏移分量初始時刻的大小取決于故障時刻的入射角θn，直流偏移分量的衰減時間常數(shù)取決于系統(tǒng)在故障點的比值X/R。

實際上，諧波成分比較復(fù)雜，系統(tǒng)發(fā)生單相接地故障，其幅值與基頻分量相比很小，諧波只會造成數(shù)據(jù)的離散并不能改變基頻信號的趨勢?；谧钚《饲€擬合方法構(gòu)建擬合函數(shù)，將故障電流信號表示為一個對稱的正弦基頻交流分量和一個隨時間常數(shù)呈指數(shù)衰減的直流偏移分量：

2 故障仿真實驗與基頻分量提取算法

本節(jié)在MATLAB/Simulink仿真環(huán)境下搭建如圖1所示的雙端電源系統(tǒng)輸電線路單相接地故障仿真模型，仿真模型參數(shù)為：兩端電源電壓EG=220∠30°kV、EH=220∠0°kV，系統(tǒng)頻率50Hz，輸電線路總長度300km，設(shè)置單相接地故障點F距G端150km，兩端電源內(nèi)阻ZG=1.06+j44.92Ω、ZH=1.05+j43.18Ω，線路正序阻抗Z1=0.02083+j0.282Ω/km，線路零序阻抗Z0=0.1148+j0.7186Ω/km。仿真時間0.2s，采樣頻率為5kHz，即每個周波采樣100點。

圖1 輸電線路故障模型

過渡電阻為零時，在監(jiān)測位置G記錄的三相電流和電壓波形如圖2所示。

圖2 記錄電壓和電流波形

2.1 檢測故障信號

系統(tǒng)發(fā)生單相接地故障，微機(jī)繼電保護(hù)裝置通過阻抗繼電器的視在阻抗或過流繼電器的均方根電流與預(yù)先設(shè)定好的阻抗域或繼電器吸合電流作比較來檢測故障。上述兩種方法將窗口時域波形轉(zhuǎn)換為阻抗復(fù)數(shù)平面[18]或均方根域，能夠成功地檢測故障發(fā)生和故障清除。

基于過流繼電器的工作原理，設(shè)定兩倍的均方根電流作為過流繼電器預(yù)先設(shè)定的均方根域。如圖3所示，故障電流信號首次超過均方根域時檢測到故障發(fā)生，對應(yīng)的故障起始時刻為0.0612s；故障電流信號末次高于均方根域后故障清除，對應(yīng)的故障清除時刻為0.1404s。

圖3 均方根域檢測故障電流

2.2 獲取完整的故障時間段數(shù)據(jù)

獲取完整的故障時間段數(shù)據(jù)應(yīng)用于提取基頻分量有兩個優(yōu)點：1）當(dāng)故障持續(xù)時間較短、特別是小于一個周期時，因為故障數(shù)據(jù)點較少，所以使用完整的故障時間段數(shù)據(jù)提取基頻分量準(zhǔn)確度更高。2）能夠最大限度地減少非故障期間數(shù)據(jù)點。

本文以均方根域檢測到的故障時刻為依據(jù)，將常見的時頻域分析工具小波變換應(yīng)用于暫態(tài)電能質(zhì)量分析中，能夠準(zhǔn)確地檢測與定位故障信號的起止時刻與持續(xù)時間。當(dāng)小波變換應(yīng)用于電能質(zhì)量檢測時，小波基的選取是一項重要的工作，文獻(xiàn)[19]和文獻(xiàn)[20]提出db小波是檢測電能質(zhì)量擾動的一個很好的選擇。如圖4所示，當(dāng)故障電流信號首次超過均方根域，檢測到故障發(fā)生。選擇小波系數(shù)能量集中的db4為小波基進(jìn)行小波變換，平方小波變換系數(shù)最大值對應(yīng)的0.0608s為校正后的故障起始時刻。當(dāng)故障電流信號末次高于均方根域后，故障清除，采用同樣的方法得到校正后的故障清除時刻為0.1438s。

圖4 使用db4小波變換校正故障時間

故障持續(xù)時間為0.0830s，提取完整的故障發(fā)生時段電流信號如圖5所示。

圖5 獲取的故障電流信號

2.3 求解模型參數(shù)

根據(jù)計算機(jī)數(shù)值計算的需要，引出向量a=(β1，β2，β3，β4)表示模型參數(shù)，式(2)改寫為：

式中：β1=A1表示對稱的基頻交流分量幅值，單位A；β2=θ1表示對稱的基頻交流分量初相位，單位rad；β2=B表示指數(shù)衰減的直流偏移分量幅值，單位A；β4-1=-ω/(X/R)表示指數(shù)衰減的直流偏移衰減常數(shù)，單位s-1。

通過對故障信號電氣量特點的分析，將故障時刻ti對應(yīng)的單相接地故障電流vi進(jìn)行曲線擬合求解模型參數(shù)。式(4)定義有約束的最小二乘目標(biāo)函數(shù)：使用信賴域反射法[21～23]求解式(4)，得到基頻分量β1cos(2πfti+β2)和直流偏移分量β3e-β4ti。

2.4 提取電氣量

本文針對以下四種典型的故障情況按照上述步驟進(jìn)行分析：理想情況下的單相接地短路故障、含有諧波分量的單相接地故障、持續(xù)時間小于一個周期的單相接地短路故障和存在過渡電阻的單相接地故障。

2.4.1 單相接地短路故障

理想情況下的單相接地短路故障曲線擬合結(jié)果如圖6所示，基頻分量提取結(jié)果如表1所示。

圖6 分解故障電流

2.4.2 含有諧波分量的單相接地短路故障

電網(wǎng)中由于正弦電壓加壓于非線性負(fù)載，導(dǎo)致電流發(fā)生畸變產(chǎn)生諧波，一般以3，6n±1次諧波較為常見。其中，3次諧波是主要成分之一且幅值最大，對繼電保護(hù)和設(shè)備的沖擊最強(qiáng)，對電網(wǎng)的影響也最為嚴(yán)重。

在故障信號中分別添加不同類型的諧波分量，其幅值為基波幅值的10%，基頻分量參數(shù)提取結(jié)果如表1所示。

表1 含有不同諧波分量的計算結(jié)果

將3，5，7，11，13，17，19次諧波同時添加至故障電流信號，改變諧波幅值，基頻分量參數(shù)提取結(jié)果如表2所示。

表2 不同幅值時的計算結(jié)果

當(dāng)諧波幅值為基波幅值的30%，曲線擬合結(jié)果如圖7所示。

歐拉iQ高續(xù)航版也在車展現(xiàn)場上市，其綜合續(xù)航里程達(dá)401公里，不同配置車型補(bǔ)貼后售價分別為9.28萬元、9.98萬元和10.88萬元。與此同時，歐拉還對外發(fā)布了汽車新零售平臺“歐拉ORA”超級APP?！皻W拉ORA”超級APP可實現(xiàn)線上下單功能，資訊、話題和社區(qū)等線性連接功能，與歐拉粉絲建立起立體連接。

圖7 分解含有諧波分量的故障電流

由表1可知，當(dāng)故障信號只含有5次諧波時，基頻分量幅值的誤差最大；隨著諧波次數(shù)增大，誤差逐漸減小。當(dāng)故障信號只含有3次諧波時，基頻分量初相位的誤差最大；隨著諧波次數(shù)的增大，誤差減小。由表2可知，當(dāng)故障信號中含有多個諧波分量時，隨著諧波幅值的逐漸增大，提取的基頻分量幅值和初相位的誤差也逐漸增大。

然而，即使使用含有多個諧波分量且諧波幅值為基波幅值30%的故障信號提取基頻分量參數(shù)，基頻分量幅值和初相位的誤差也僅為-0.2521%、0.2088%。擬合出的基頻分量參數(shù)與無諧波時相比誤差很小，證明本文建立的最小二乘函數(shù)模型能夠較好的抑制諧波分量，基頻分量提取精度較高。

2.4.3 持續(xù)時間小于一個周期的單相接地短路故障

微機(jī)繼電保護(hù)裝置使用傅里葉濾波器提取基頻相量需要一個周期的波形數(shù)據(jù)。如圖8(a)所示，當(dāng)故障持續(xù)時間小于一個周期，故障期間和故障結(jié)束后的數(shù)據(jù)點都將應(yīng)用于傅里葉濾波器，導(dǎo)致提取的基頻分量精度降低、故障定位誤差增大。曲線擬合結(jié)果如圖8(b)所示，電氣量提取結(jié)果如表3所示。

圖8 故障持續(xù)時間不足一個周期

2.4.4 存在過渡電阻的單相接地故障

當(dāng)單相接地故障存在過渡電阻時，隨著過渡電阻阻值增大，系統(tǒng)中出現(xiàn)直流偏移分量的可能性逐漸減小。由于阻值的增加降低了故障期間的X/R值，所以產(chǎn)生的直流偏移分量衰減速度快。

設(shè)定過渡電阻為50Ω時系統(tǒng)發(fā)生單相接地故障，故障仿真信號如圖9(a)所示，曲線擬合結(jié)果如圖9(b)所示，電氣量提取結(jié)果如表3所示。

圖9 存在過渡電阻的單相接地故障

表3 不同情況的計算結(jié)果

3 算法驗證分析

微機(jī)繼電保護(hù)裝置基于阻抗法的Takagi方法計算從監(jiān)測點到故障位置的距離：

式中，U表示母線電壓降，I表示故障電流，Iper表示故障前的電流，I0表示零序電流，Z1為線路正序阻抗，Z0為線路零序阻抗，*表示共軛復(fù)數(shù)。

同時，故障定位結(jié)果誤差定義如下：

本節(jié)通過故障定位結(jié)果，驗證在上述四種典型的故障情況下提取基頻分量的準(zhǔn)確性。

3.1 含有諧波分量的單相接地短路故障

由表4可知，含有諧波分量的單相接地短路故障的故障定位誤差很小，說明提取的基頻分量魯棒性強(qiáng)、準(zhǔn)確度高，能夠很好的消除單相接地故障產(chǎn)生的直流偏移分量。

3.2 故障持續(xù)時間小于一個周期的單相接地短路

經(jīng)計算，與單相接地短路故障相比誤差增大。主要原因有兩點：1）采樣點數(shù)據(jù)較少會導(dǎo)致基頻相量計算結(jié)果誤差增大。2）故障數(shù)據(jù)中參雜了非故障期間的數(shù)據(jù)點。

3.3 存在過渡電阻的單相接地故障

由表4可知，單相接地故障存在過渡電阻時，與單相接地短路故障相比誤差增大。此時，誤差源自Takagi方法假設(shè)系統(tǒng)是均勻的，故障點兩端的故障電流分量相位角一致[24]。然而，在仿真環(huán)境下系統(tǒng)是非均勻的，違反假設(shè)導(dǎo)致誤差增大。但是算法可以避免提取基頻相量時造成額外誤差的情況，基頻分量提取結(jié)果準(zhǔn)確度較高。

表4 故障距離和誤差

4 結(jié)語

本文構(gòu)建的故障信號數(shù)學(xué)模型能夠快速消除單相接地故障產(chǎn)生的直流偏移分量、準(zhǔn)確提取故障信號基頻分量，故障定位結(jié)果精度較高。數(shù)學(xué)模型具有以下四個優(yōu)點：

1）運用函數(shù)逼近原理，把故障信號看作基頻交流分量和直流偏移分量相疊加，使用數(shù)學(xué)表達(dá)式直接消除直流偏移分量，提取基頻分量更加簡單。

2）數(shù)學(xué)模型通常在10～200次迭代后提取基頻分量，算法擁有快速高效的計算效率。

3）諧波分量對現(xiàn)有數(shù)學(xué)模型擬合出的基頻相量影響很小，對諧波分量有較好的抑制作用。

4）可以計算任意時長瞬時故障信號包含的基頻相量，算法不受波形周期的限制，準(zhǔn)確度較高。