葉洪波， 姜陽， 陳雪梅， 崔勇， 俞越， 陸超

(1. 國網(wǎng)上海市電力公司，上海 200122；2. 清華大學(xué)電機(jī)工程與應(yīng)用電子技術(shù)系，北京 100084；3. 上海電力設(shè)計(jì)院有限公司，上海 200025)

0 引言

電力系統(tǒng)的慣量響應(yīng)對(duì)頻率調(diào)節(jié)有重要影響[1—4]。實(shí)時(shí)監(jiān)測(cè)同步發(fā)電機(jī)慣量是電力系統(tǒng)頻率安全穩(wěn)定在線分析與控制的重要基礎(chǔ)。一般采用甩負(fù)荷實(shí)驗(yàn)離線測(cè)量發(fā)電機(jī)慣量，但轉(zhuǎn)子渦動(dòng)和測(cè)量噪聲會(huì)降低結(jié)果的可信度[5—6]。而在線辨識(shí)能更好地反映發(fā)電機(jī)在實(shí)際工況下的動(dòng)態(tài)特性[7]，因此采用在線辨識(shí)方法估計(jì)發(fā)電機(jī)慣量。同步相量測(cè)量裝置(phasor measurement unit,PMU)可在電力系統(tǒng)正常運(yùn)行中實(shí)時(shí)測(cè)量豐富的類噪聲數(shù)據(jù)，反映電力系統(tǒng)動(dòng)態(tài)特性[8]，為在線監(jiān)測(cè)發(fā)電機(jī)慣量提供數(shù)據(jù)條件[9]。使用系統(tǒng)辨識(shí)，以搖擺方程為基礎(chǔ)，實(shí)時(shí)估計(jì)PMU功率和頻率擾動(dòng)信號(hào)的動(dòng)態(tài)模型，并從中獲取發(fā)電機(jī)慣量的方法已得到較多研究[10—12]。

這類辨識(shí)方法對(duì)于仿真系統(tǒng)中產(chǎn)生的數(shù)據(jù)有較好的效果[10，12]，但在實(shí)測(cè)數(shù)據(jù)中的適用性和可靠性一直難以保證[13—14]。其原因一是實(shí)際電網(wǎng)高度復(fù)雜，而現(xiàn)有模型完備性有所欠缺；二是數(shù)據(jù)質(zhì)量存在問題[15]，使用壞數(shù)據(jù)或精度較低的數(shù)據(jù)進(jìn)行辨識(shí)會(huì)得到錯(cuò)誤的結(jié)果；三是類噪聲數(shù)據(jù)信噪比較低[16]，導(dǎo)致部分?jǐn)?shù)據(jù)信息不足而無法辨識(shí)。

處理實(shí)測(cè)數(shù)據(jù)是提升辨識(shí)效果的重要手段。數(shù)據(jù)的修復(fù)和篩選是處理實(shí)測(cè)數(shù)據(jù)的2種途徑?，F(xiàn)有基于類噪聲的參數(shù)辨識(shí)工作多采用異常數(shù)據(jù)檢測(cè)、恢復(fù)、降采樣和去噪等方法[17—21]，然而，在數(shù)據(jù)質(zhì)量較低時(shí)，不能保證修復(fù)后數(shù)據(jù)的可靠性。數(shù)據(jù)的篩選多采用頻域分析方法，其目的在于篩選出可以反映系統(tǒng)動(dòng)態(tài)特性的數(shù)據(jù)段。文獻(xiàn)[22]指出，輸入信號(hào)頻譜覆蓋振蕩頻段時(shí)，才能較準(zhǔn)確地辨識(shí)負(fù)荷參數(shù)。用類噪聲信號(hào)辨識(shí)系統(tǒng)振蕩模式時(shí)，文獻(xiàn)[9，23]基于類噪聲信號(hào)的功率譜在關(guān)心頻段的峰度指標(biāo)進(jìn)行數(shù)據(jù)篩選，但關(guān)心頻段的選取主觀性較強(qiáng)，且可解釋性較差。慣量辨識(shí)領(lǐng)域尚無有關(guān)數(shù)據(jù)篩選的文獻(xiàn)。

文中針對(duì)功頻一階模型，基于類噪聲信號(hào)的譜分析和系統(tǒng)辨識(shí)理論，建立慣量估計(jì)的先驗(yàn)方差統(tǒng)計(jì)量作為可辨識(shí)指標(biāo)，反映不同數(shù)據(jù)窗用于辨識(shí)的效果，實(shí)現(xiàn)對(duì)數(shù)據(jù)的評(píng)估和篩選。文中和現(xiàn)有數(shù)據(jù)篩選文獻(xiàn)的思路均為類噪聲信號(hào)的頻域分析，但文中方法的理論性、可解釋性與適應(yīng)性較強(qiáng)，仿真和實(shí)測(cè)數(shù)據(jù)的算例驗(yàn)證了文中方法的有效性。

1 辨識(shí)模型和算法與評(píng)價(jià)指標(biāo)

針對(duì)發(fā)電機(jī)機(jī)端量測(cè)的類噪聲數(shù)據(jù)，采用基于系統(tǒng)辨識(shí)的方法實(shí)現(xiàn)慣量監(jiān)測(cè)。辨識(shí)的模型是發(fā)電機(jī)穩(wěn)定工作點(diǎn)附近建立的小信號(hào)模型，其輸入為離散時(shí)間的標(biāo)幺有功擾動(dòng)，記為u(n)；輸出為標(biāo)幺頻率擾動(dòng)，記為y(n)。文中該模型結(jié)構(gòu)固定為一階的輸出誤差(output error,OE)模型，如式(1)所示，各項(xiàng)的具體形式如式(2)所示。

y(n)=G(z,θ)u(n)+H(z,θ)e(n)

(1)

(2)

式中：z-1為移位算子；A，B為2個(gè)模型參數(shù)，將其合并記為θ；e(n)為白噪聲。對(duì)PMU量測(cè)的類噪聲數(shù)據(jù)截取長度為N的數(shù)據(jù)窗。對(duì)于式(1)、式(2)所述OE模型，采用預(yù)報(bào)誤差算法進(jìn)行優(yōu)化求解。

ypre(n)=-Aypre(n-1)+Bu(n-1)

(3)

(4)

式中：ypre(n)為對(duì)輸出的預(yù)報(bào)；θN為長度為N的數(shù)據(jù)窗辨識(shí)出的模型參數(shù)。忽略噪聲，將系統(tǒng)傳遞函數(shù)G(z,θ)雙線性變換連續(xù)化，得到連續(xù)系統(tǒng)傳遞函數(shù)，最終得到慣量M：

(5)

式中：h(0+)為連續(xù)系統(tǒng)沖激響應(yīng)在0+時(shí)刻的值；TS為PMU中2幀數(shù)據(jù)的時(shí)間間隔；M(θ)為數(shù)據(jù)窗的慣量估計(jì)。盡管一階OE模型可能不足以準(zhǔn)確建模發(fā)電機(jī)的調(diào)速器特性，但該模型結(jié)構(gòu)的辨識(shí)結(jié)果數(shù)值穩(wěn)定性較高。而采用更高階的模型較易出現(xiàn)數(shù)值問題，得到大量錯(cuò)誤的辨識(shí)結(jié)果，導(dǎo)致整個(gè)數(shù)據(jù)集的辨識(shí)成功率降低。

定義數(shù)據(jù)窗的慣量辨識(shí)誤差為ε，擬合度為η。2個(gè)量都描述一組數(shù)據(jù)窗的辨識(shí)效果，擬合度越高，模型預(yù)測(cè)輸出和實(shí)際輸出越接近。ε和η的計(jì)算如下：

ε=(|M-M0|/M0)×100%

(6)

(7)

式中：M0為慣量真值；M為該數(shù)據(jù)窗慣量辨識(shí)值。ypre(k)根據(jù)辨識(shí)出的θN使用式(3)進(jìn)行遞推得到。設(shè)某數(shù)據(jù)集中數(shù)據(jù)窗總數(shù)為NT。定義有效誤差門限為εTH，并認(rèn)為滿足ε<εTH的數(shù)據(jù)窗為有效辨識(shí)。統(tǒng)計(jì)數(shù)據(jù)集有效辨識(shí)的數(shù)據(jù)窗個(gè)數(shù)NS，并定義辨識(shí)有效率ξ。ξ的計(jì)算如下：

ξ=(NS/NT)×100%

(8)

對(duì)數(shù)據(jù)集中各個(gè)數(shù)據(jù)窗的辨識(shí)誤差和擬合度取平均，得到平均誤差εav和平均擬合度ηav。平均誤差、平均擬合度和辨識(shí)有效率可描述辨識(shí)方法對(duì)數(shù)據(jù)集總體的辨識(shí)效果。

2 先驗(yàn)方差與數(shù)據(jù)篩選

系統(tǒng)辨識(shí)理論分析需要以一個(gè)理論上的“真實(shí)系統(tǒng)”為基礎(chǔ)，該系統(tǒng)應(yīng)反映物理上功率和頻率的關(guān)系。調(diào)速器動(dòng)作時(shí)，真實(shí)系統(tǒng)結(jié)構(gòu)如圖1所示[12]。其中，Gp(s)為調(diào)速器的動(dòng)態(tài)特性傳遞函數(shù)；Δpm(s)為發(fā)電機(jī)機(jī)械功率擾動(dòng)；M為發(fā)電機(jī)慣量；D為發(fā)電機(jī)阻尼系數(shù)。

圖1 調(diào)速器-轉(zhuǎn)子系統(tǒng)結(jié)構(gòu)

該系統(tǒng)為閉環(huán)系統(tǒng)，考慮了調(diào)速器的動(dòng)態(tài)特性，即頻率對(duì)機(jī)械功率的影響，而搖擺方程反映了功率對(duì)頻率的影響。電機(jī)正常運(yùn)行時(shí)，頻率波動(dòng)較小，可認(rèn)為調(diào)速器不動(dòng)作，搖擺方程中機(jī)械功率不變，發(fā)電機(jī)的有功功率擾動(dòng)和頻率擾動(dòng)滿足一階微分方程，因此假設(shè)離散化的真實(shí)系統(tǒng)結(jié)構(gòu)為一階OE模型是合理的。采用OE模型求解時(shí)，若辨識(shí)的模型結(jié)構(gòu)和真實(shí)系統(tǒng)結(jié)構(gòu)相同，則數(shù)據(jù)窗長度N趨于無窮時(shí)，按照式(3)和式(4)估計(jì)的模型參數(shù)θN依概率1收斂于真實(shí)系統(tǒng)的參數(shù)θ0；同時(shí)參數(shù)辨識(shí)θN服從漸近正態(tài)分布，其漸近期望為θ0，漸近方差記為cov(θN)[19]。cov(θN)的計(jì)算如下：

(9)

其中：

(10)

式中：Φu(ω)為u(n)的功率譜；Φue(ω)，Φeu(ω)為u(n)和e(n)的互功率譜；λ0為白噪聲的方差[19]。

因此可認(rèn)為所述的一階模型辨識(shí)結(jié)果是依概率1收斂至真實(shí)系統(tǒng)參數(shù)的。而不同數(shù)據(jù)窗得到辨識(shí)結(jié)果誤差的波動(dòng)較大這一現(xiàn)象理解為相應(yīng)數(shù)據(jù)得到的慣量估計(jì)量方差不同。產(chǎn)生較大方差的輸入數(shù)據(jù)，其辨識(shí)結(jié)果偏差較大的可能性更高。

每個(gè)獨(dú)立的數(shù)據(jù)窗的量測(cè)數(shù)據(jù)均可看作具有相應(yīng)功率譜的隨機(jī)過程的一個(gè)實(shí)現(xiàn)，進(jìn)而估計(jì)該隨機(jī)過程各個(gè)不同實(shí)現(xiàn)的數(shù)據(jù)用于辨識(shí)后得到的慣量估計(jì)的方差。該方差體現(xiàn)了辨識(shí)結(jié)果的分散性，分散性越大，隨機(jī)過程的每個(gè)特定實(shí)現(xiàn)進(jìn)行參數(shù)辨識(shí)得到誤差較大的值的可能性越高。慣量估計(jì)的先驗(yàn)方差由估計(jì)參考系統(tǒng)的參數(shù)、估計(jì)模型參數(shù)的方差和估計(jì)慣量的方差3個(gè)步驟實(shí)現(xiàn)。

2.1 估計(jì)參考系統(tǒng)的參數(shù)

由式(9)可知，求取模型參數(shù)的漸近方差時(shí)，需要已知真實(shí)系統(tǒng)的參數(shù)，因此可以通過預(yù)辨識(shí)估計(jì)一個(gè)參考系統(tǒng)參數(shù)θ1來代替真實(shí)系統(tǒng)，以估計(jì)方差。具體過程為：對(duì)一段數(shù)據(jù)集中的各個(gè)數(shù)據(jù)窗，進(jìn)行慣量辨識(shí)并統(tǒng)計(jì)擬合度。一般認(rèn)為，擬合度越高，辨識(shí)結(jié)果越準(zhǔn)確。因此設(shè)置一個(gè)擬合度的門限值ηTH，對(duì)數(shù)據(jù)集中擬合度高于ηTH的數(shù)據(jù)窗的模型參數(shù)各取平均值。為實(shí)現(xiàn)ηTH選取的自適應(yīng)性，ηTH可以選取為各數(shù)據(jù)窗辨識(shí)擬合度的某分位數(shù)。

θ1=mean({θNi|ηi>ηTH})

(11)

式中：θNi，ηi分別為數(shù)據(jù)集中第i個(gè)數(shù)據(jù)窗所得模型參數(shù)和擬合度；mean(·)為對(duì)各向量取平均值。

2.2 模型參數(shù)的方差

假設(shè)真實(shí)系統(tǒng)結(jié)構(gòu)為一階OE模型，則u(n)與e(n)是獨(dú)立的隨機(jī)變量，Φue(ω)≡Φeu(ω)≡0。數(shù)據(jù)窗長度N相同時(shí)，可不考慮N對(duì)參數(shù)估計(jì)方差的影響。由于OE模型的特殊噪聲結(jié)構(gòu)，T′(z,θ)第二列為0，式(9)化簡(jiǎn)得到R矩陣，可反映數(shù)據(jù)窗辨識(shí)得到的模型參數(shù)θN的方差。

(12)

其中：

(13)

各組數(shù)據(jù)中輸入功率譜可用周期圖估計(jì)：

(14)

使用2.1節(jié)估計(jì)出的θ1，可估計(jì)出模型參數(shù)的方差。實(shí)際計(jì)算時(shí)，式(12)由數(shù)值積分得到。

2.3 慣量的方差

由于模型參數(shù)θ和慣量估計(jì)M具有式(5)所述的函數(shù)關(guān)系，因此可以根據(jù)辨識(shí)得到的θN的方差，進(jìn)一步估計(jì)慣量估計(jì)M的方差。假設(shè)辨識(shí)出的θN≈θ1，則根據(jù)高斯近似表達(dá)式可得到慣量先驗(yàn)方差的估計(jì)量σ[19]。

(15)

σ和慣量的方差成正比關(guān)系，在各數(shù)據(jù)窗長度N相同時(shí)，可直接用σ衡量一個(gè)數(shù)據(jù)窗慣量估計(jì)方差的大小。

2.4 篩選方法

由于不同機(jī)組先驗(yàn)方差變化較大，因此采用自適應(yīng)方法。對(duì)一段較長時(shí)間(例如15 min)的PMU數(shù)據(jù)集，先均勻取等長數(shù)據(jù)窗，計(jì)算各數(shù)據(jù)窗先驗(yàn)方差并排序，取出先驗(yàn)方差較小的那部分?jǐn)?shù)據(jù)窗(經(jīng)驗(yàn)上取前10%或25%)，得到一個(gè)新數(shù)據(jù)集。一般而言，新數(shù)據(jù)集的總體辨識(shí)效果優(yōu)于原數(shù)據(jù)集。

2.5 討論

事實(shí)上，2.1節(jié)中估計(jì)參考系統(tǒng)參數(shù)的步驟與2.2～2.3節(jié)中估計(jì)先驗(yàn)方差的步驟相獨(dú)立，用于估計(jì)參考系統(tǒng)的數(shù)據(jù)集與篩選的數(shù)據(jù)集可以不一致。參考系統(tǒng)更新的頻度低于求取慣量估計(jì)方差的頻度。在參考系統(tǒng)固定且已經(jīng)先驗(yàn)求出的條件下，式(12)和式(15)中帶有θ1的項(xiàng)均已固定，不同組數(shù)據(jù)的先驗(yàn)方差σ只與其模型輸入數(shù)據(jù)u的功率譜Φu(ω)相關(guān)。每個(gè)數(shù)據(jù)窗求先驗(yàn)方差時(shí)的主要運(yùn)算為該數(shù)據(jù)窗功率譜的估計(jì)和數(shù)值積分的計(jì)算，算法實(shí)現(xiàn)時(shí)那些固定的值可事先算好并儲(chǔ)存。Φu(ω)表征擾動(dòng)的能量，一般而言，擾動(dòng)的能量越大，辨識(shí)的方差越小，辨識(shí)結(jié)果越準(zhǔn)確，這與式(12)相一致，即Φu(ω)的值越大，θN的方差越小。定義加權(quán)矩陣W(ω)：

W(ω)=K(ejω,θ1)(K(e-jω,θ1))T

(16)

式(12)可理解為采用W(ω)給Φu(ω)在不同頻率加權(quán)并積分。經(jīng)簡(jiǎn)易的推導(dǎo)可知，在頻率較高時(shí)，W(ω)各項(xiàng)的值迅速衰減，因此先驗(yàn)方差對(duì)各個(gè)頻段的敏感性不一致，對(duì)較低頻段尤為敏感。而類噪聲信號(hào)的能量也多集中于低頻，因此只需要關(guān)心輸入功率譜在較低頻的幅值。文中方法是頻域篩選方法的延續(xù)與規(guī)范化。與現(xiàn)有文獻(xiàn)中的方法相比，通過估計(jì)參考系統(tǒng)與加權(quán)矩陣可自適應(yīng)地確定關(guān)心頻段的位置與不同頻率的權(quán)值。文中方法的特點(diǎn)為：基于先驗(yàn)方差，在辨識(shí)前即可對(duì)數(shù)據(jù)進(jìn)行篩選，從而不辨識(shí)那些被去除的數(shù)據(jù)窗，提升辨識(shí)的效率與準(zhǔn)確度。

3 數(shù)據(jù)評(píng)估篩選算法的驗(yàn)證

將上述數(shù)據(jù)篩選方法用于處理仿真和實(shí)測(cè)數(shù)據(jù)，以檢驗(yàn)該方法的有效性。

3.1 仿真數(shù)據(jù)算例

使用基于Matlab的電力系統(tǒng)時(shí)域仿真軟件PSAT得到電力系統(tǒng)類噪聲工況的數(shù)據(jù)。仿真系統(tǒng)為標(biāo)準(zhǔn)的IEEE 10機(jī)39節(jié)點(diǎn)系統(tǒng)，系統(tǒng)結(jié)構(gòu)如圖2所示，其中G1—G10為發(fā)電機(jī)[24—25]。對(duì)26號(hào)節(jié)點(diǎn)的有功功率施加標(biāo)準(zhǔn)差為0.05倍均值的典型類噪聲擾動(dòng)，該類噪聲信號(hào)為電網(wǎng)中實(shí)測(cè)有功功率信號(hào)。

圖2 IEEE 39節(jié)點(diǎn)系統(tǒng)結(jié)構(gòu)

測(cè)量有功擾動(dòng)下各個(gè)發(fā)電機(jī)端的有功功率和頻率，共進(jìn)行6 min仿真。固定每個(gè)數(shù)據(jù)窗長1 min，實(shí)時(shí)傳送速率為25次/s，即N為1 500。辨識(shí)前，對(duì)各數(shù)據(jù)窗進(jìn)行標(biāo)幺化、去均值等預(yù)處理得到u(n)與y(n)。為近似真實(shí)系統(tǒng)中PMU頻率量測(cè)精度不高的情況，對(duì)頻率數(shù)據(jù)進(jìn)行了舍入，將頻率的精度設(shè)為0.001 Hz(標(biāo)幺值2×10-5)，與實(shí)測(cè)數(shù)據(jù)情況相符。先取仿真數(shù)據(jù)中前80 s的數(shù)據(jù)集用于參考系統(tǒng)參數(shù)估計(jì)，之后使用該參考系統(tǒng)對(duì)后280 s的數(shù)據(jù)集進(jìn)行篩選，對(duì)各數(shù)據(jù)集中每隔2 s取一個(gè)數(shù)據(jù)窗進(jìn)行辨識(shí)。估算參考系統(tǒng)參數(shù)，其中ηTH取各組數(shù)據(jù)擬合度的95%分位數(shù)，事實(shí)上相當(dāng)于取擬合度最高的一組數(shù)據(jù)窗模型參數(shù)為θ1。各機(jī)組的慣量設(shè)定值為M0，θ1，表1為使用θ1經(jīng)式(5)得到的慣量M1及其誤差ε。

表1 各機(jī)組參考系統(tǒng)參數(shù)及其慣量誤差

可見，估計(jì)的參考系統(tǒng)慣量與設(shè)定值較為接近，參考系統(tǒng)能夠近似真實(shí)系統(tǒng)。之后對(duì)各機(jī)組后280 s的數(shù)據(jù)集各數(shù)據(jù)窗進(jìn)行慣量辨識(shí)。統(tǒng)計(jì)每個(gè)數(shù)據(jù)窗的辨識(shí)誤差ε，并進(jìn)行先驗(yàn)方差σ的估計(jì)。G5的先驗(yàn)方差和辨識(shí)誤差關(guān)系如圖3所示?？梢娤闰?yàn)方差和辨識(shí)誤差有明顯正相關(guān)關(guān)系，先驗(yàn)方差σ較小的數(shù)據(jù)組，辨識(shí)誤差ε一般也較小。對(duì)G1，G5，G6，G9，篩選出先驗(yàn)方差較低的50%數(shù)據(jù)窗，得到新的數(shù)據(jù)集。仿真數(shù)據(jù)辨識(shí)誤差較小，因此取有效誤差門限εTH為5%。

圖3 G5各組慣量辨識(shí)誤差和先驗(yàn)方差關(guān)系

表2為篩選前后數(shù)據(jù)集的數(shù)據(jù)窗數(shù)、平均擬合度、平均辨識(shí)誤差和辨識(shí)有效率。以G5為例，圖4為篩選前后模型參數(shù)分布與相對(duì)應(yīng)的誤差，圖中等高線為辨識(shí)誤差。數(shù)據(jù)篩選前后，使用的辨識(shí)模型與算法相同。由表2和圖4可得，篩選后，各機(jī)組平均誤差、辨識(shí)有效率和平均擬合度都有所提升，辨識(shí)效果提高。

表2 仿真數(shù)據(jù)篩選前后慣量辨識(shí)結(jié)果對(duì)比

圖4 G5模型參數(shù)與辨識(shí)誤差

算例還展示了文中方法的先驗(yàn)性，估計(jì)參考系統(tǒng)的數(shù)據(jù)集和用于篩選檢驗(yàn)的數(shù)據(jù)集相互分離。用文中方法篩選一個(gè)數(shù)據(jù)集時(shí)，參考系統(tǒng)可以根據(jù)先于該數(shù)據(jù)集的辨識(shí)經(jīng)驗(yàn)求得。求取參考系統(tǒng)的辨識(shí)過程與篩選過程是獨(dú)立的，篩選時(shí)完全不用對(duì)該數(shù)據(jù)集中的各數(shù)據(jù)窗進(jìn)行辨識(shí)操作。

3.2 實(shí)測(cè)數(shù)據(jù)算例

用文中方法處理電力系統(tǒng)中實(shí)測(cè)數(shù)據(jù)，選取實(shí)際電網(wǎng)中的5臺(tái)發(fā)電機(jī)的數(shù)據(jù)用于處理，分別記作1號(hào)至5號(hào)。

1號(hào)發(fā)電機(jī)典型的有功功率、頻率類噪聲信號(hào)波形如圖5所示，兩圖選取時(shí)間段相同，展示了功率和頻率同步變化?？梢杂^察到頻率信號(hào)的毛刺缺陷。為解決這一缺陷，數(shù)據(jù)預(yù)處理時(shí)在標(biāo)幺化、去均值前增加零相移低通濾波操作[20]。由于信號(hào)能量主要分布在較低的頻段，濾波操作對(duì)頻率信號(hào)有效成分的影響較弱，因此能有效地濾除波形中的毛刺與高頻噪聲。文中處理實(shí)測(cè)數(shù)據(jù)時(shí)，零相移低通濾波器的截止頻率選為0.5 Hz。

圖5 1號(hào)發(fā)電機(jī)典型有功功率和頻率波形

頻率數(shù)據(jù)的精度、實(shí)時(shí)傳送速率、數(shù)據(jù)窗長分別為0.001 Hz，25次/s，1 min。以2 s為間距均勻取數(shù)據(jù)窗。以1號(hào)發(fā)電機(jī)為例，總長為1 000 s的類噪聲數(shù)據(jù)，取前200 s的數(shù)據(jù)估計(jì)參考系統(tǒng)。后800 s中各數(shù)據(jù)窗的先驗(yàn)方差和辨識(shí)誤差如圖6所示。其中ηTH為各組辨識(shí)擬合度的90%分位數(shù)。根據(jù)圖6，先驗(yàn)方差σ和慣量辨識(shí)誤差ε具有正相關(guān)性。且實(shí)測(cè)數(shù)據(jù)的辨識(shí)中較多辨識(shí)結(jié)果較差的數(shù)據(jù)組，辨識(shí)誤差最大可以達(dá)到11 859.9%。在計(jì)算辨識(shí)誤差時(shí)，采用實(shí)際調(diào)度系統(tǒng)使用的仿真模型中相應(yīng)機(jī)組慣量作為真值。

圖6 1號(hào)發(fā)電機(jī)各組慣量辨識(shí)誤差與先驗(yàn)方差關(guān)系

每個(gè)機(jī)組選取出先驗(yàn)方差較小的前10%的數(shù)據(jù)窗作為新數(shù)據(jù)集，表3為篩選前后數(shù)據(jù)集的窗數(shù)、平均誤差、辨識(shí)有效率和平均擬合度。

表3 實(shí)測(cè)數(shù)據(jù)篩選前后慣量辨識(shí)結(jié)果

由于實(shí)測(cè)數(shù)據(jù)辨識(shí)誤差較大，因此εTH選取為30%?？梢缘玫胶头抡鏀?shù)據(jù)類似的結(jié)論。可見用文中方法先對(duì)各數(shù)據(jù)組進(jìn)行篩選，再進(jìn)行辨識(shí)可以有效提高辨識(shí)的質(zhì)量，其中1號(hào)發(fā)電機(jī)和2號(hào)發(fā)電機(jī)篩選后辨識(shí)精度得到顯著提升。

加權(quán)矩陣W(ω)是一個(gè)2×2的復(fù)矩陣，記其第i行第j列的元素為wij(ω)。圖7為W(ω)矩陣各個(gè)元素的虛實(shí)部值與ω的關(guān)系。

圖7 1號(hào)發(fā)電機(jī)先驗(yàn)加權(quán)矩陣W的各元素

可見，對(duì)1號(hào)發(fā)電機(jī)，在辨識(shí)之前應(yīng)該關(guān)注的是輸入數(shù)據(jù)中非常低頻的有功擾動(dòng)，該頻段有功擾動(dòng)越大，數(shù)據(jù)窗辨識(shí)結(jié)果較好；ω大于0.01π rad/s以后，加權(quán)矩陣的值幾乎為0，考慮實(shí)時(shí)傳送速率為25次/s，該ω對(duì)應(yīng)的頻率為0.125 Hz。上述結(jié)果表明，對(duì)發(fā)電機(jī)，低于0.125 Hz的頻段是頻率擾動(dòng)u的有效頻段。為了節(jié)省計(jì)算時(shí)間，在計(jì)算式(12)時(shí)可以只在低頻段積分(如[-0.02π, 0.02π])，而不用在[-π,π]整個(gè)區(qū)間積分。

實(shí)測(cè)數(shù)據(jù)的辨識(shí)效果較差，有辨識(shí)誤差極大的組。在這種情況下，若不進(jìn)行數(shù)據(jù)篩選，直接進(jìn)行辨識(shí)，則第1章的辨識(shí)方法不能實(shí)際應(yīng)用。但若采用第2章所述篩選方法，可在辨識(shí)之前，有效地篩選出辨識(shí)結(jié)果較好的數(shù)據(jù)窗。

4 結(jié)論

文中在已有的類噪聲數(shù)據(jù)辨識(shí)慣量方法的基礎(chǔ)上，從概率的角度建立了先驗(yàn)方差統(tǒng)計(jì)量，并將其應(yīng)用于數(shù)據(jù)篩選。在進(jìn)行辨識(shí)運(yùn)算之前，即可先驗(yàn)地從大量數(shù)據(jù)中選取出辨識(shí)效果較好的一些數(shù)據(jù)窗，為辨識(shí)算法的實(shí)施提供重要幫助。文中將所提方法應(yīng)用于處理仿真和實(shí)測(cè)數(shù)據(jù)，檢驗(yàn)了方法的有效性。但由于實(shí)測(cè)數(shù)據(jù)的復(fù)雜性及概率的不確定性等原因，篩選后辨識(shí)結(jié)果雖有所提升，但仍有不盡人意的地方，這也是下一步研究中亟需解決的問題。

本文得到國網(wǎng)上海市電力公司科技項(xiàng)目(SGSH0000DKJS1900681)資助，謹(jǐn)此致謝！