于 淼, 胡敬軒, 張壽志, 魏靜靜, 孫建群, 吳屹瀟

(1.北京建筑大學(xué) 機(jī)電與車輛工程學(xué)院,北京 100044;2.北京市建筑安全檢測工程技術(shù)研究中心,北京100044; 3.清華大學(xué) 電力系統(tǒng)及發(fā)電設(shè)備安全控制和仿真國家重點(diǎn)實驗室,北京 100084)

隨著電力系統(tǒng)核心設(shè)備電力電子化發(fā)展以及新能源發(fā)電快速變革,當(dāng)前電力系統(tǒng)呈現(xiàn)出向高比例可再生能源以及電力電子設(shè)備“雙高”發(fā)展態(tài)勢.同時輸電網(wǎng)中高壓大容量變流設(shè)備持續(xù)推廣以及配用電側(cè)電力電子技術(shù)廣泛應(yīng)用,將逐步成為新一代電力系統(tǒng)重要技術(shù)特征.在新型電力系統(tǒng)中,大型同步發(fā)電機(jī)組被小型電源替代,海量分布式電源、電動汽車、分布式儲能等以電力電子為接口的新型設(shè)備從中低壓側(cè)接入系統(tǒng),使得穩(wěn)定性分析中需要考慮動態(tài)元件數(shù)目可能達(dá)到十萬到百萬量級[1],并且各類設(shè)備動態(tài)特性差異相對較大,呈現(xiàn)出高度異構(gòu)性.對新型電力系統(tǒng)小干擾穩(wěn)定性分析需要解決由于維數(shù)爆炸引發(fā)的計算負(fù)擔(dān)、通信負(fù)擔(dān)等嚴(yán)峻挑戰(zhàn)[2].此外,電力系統(tǒng)變化時間尺度在穩(wěn)定性分析模型中可以視為恒定時間,現(xiàn)有穩(wěn)定性分析理論與工具多是針對給定平衡點(diǎn)進(jìn)行,當(dāng)平衡點(diǎn)變化時,穩(wěn)定性分析需要重新進(jìn)行.新能源供能所引起的波動以及來自大量電力電子設(shè)備頻繁投切和系統(tǒng)拓?fù)浣Y(jié)構(gòu)高速變化,使得新型電力系統(tǒng)運(yùn)行位置變化更加快速[3].因此,深入研究新型電力系統(tǒng)小干擾穩(wěn)定性機(jī)理迫在眉睫.

從我國未來新型電力系統(tǒng)“雙碳”目標(biāo)方面出發(fā),肖先勇等[4]闡述未來新型電力系統(tǒng)對“雙碳”貢獻(xiàn).在此基礎(chǔ)上,論述構(gòu)建新型電力系統(tǒng)關(guān)鍵技術(shù),總結(jié)建設(shè)新型電力系統(tǒng)將要面對的主要挑戰(zhàn).文獻(xiàn)[5]中針對我國“30·60雙碳”目標(biāo)和新型電力系統(tǒng)建設(shè)背景,闡述我國燃?xì)獍l(fā)電發(fā)展前景,提出我國新型電力系統(tǒng)安全低碳發(fā)展若干建議.文獻(xiàn)[6]中闡述在“雙碳”目標(biāo)下研究清潔電力發(fā)展路徑具有重要意義,分析我國水、風(fēng)、光等清潔能源儲量及其特征,預(yù)測中長期電力需求,提出我國未來清潔電力發(fā)展的對策與建議.

在新型電力系統(tǒng)穩(wěn)定性分析研究方面,文獻(xiàn)[7]中從分析“雙高”電力系統(tǒng)特點(diǎn)出發(fā),對經(jīng)典穩(wěn)定性進(jìn)行概述以及簡要?dú)w納,探討經(jīng)典穩(wěn)定性分類適用性,最后嘗試提出一種新的穩(wěn)定性分類方法,以求在經(jīng)典分類邏輯上兼容“雙高”電力系統(tǒng)穩(wěn)定性分析.文獻(xiàn)[8]中對電力系統(tǒng)由于小干擾產(chǎn)生失穩(wěn)現(xiàn)象進(jìn)行總結(jié)與歸納,再結(jié)合“雙高”系統(tǒng)特性與需求,探討現(xiàn)有研究理論方法不足.文獻(xiàn)[9]中基于電力系統(tǒng)提出寬頻振蕩廣域監(jiān)測與預(yù)警系統(tǒng)分析方法,并通過模擬振蕩場景驗證該系統(tǒng)有效性.文獻(xiàn)[10-12]中主要對非線性系統(tǒng)以及新型并網(wǎng)電力系統(tǒng)穩(wěn)定平衡點(diǎn)、穩(wěn)定域及建模進(jìn)行分析和研究,對直接法研究結(jié)果具有一定參考價值.文獻(xiàn)[13-16]中以新型電力系統(tǒng)為研究對象,基于李雅普諾夫直接法理論對系統(tǒng)不同擾動形式穩(wěn)定性進(jìn)行分析,并對所提方法進(jìn)行仿真驗證.文獻(xiàn)[17]中基于系統(tǒng)穩(wěn)定性分析方法,提出適用于新型電力系統(tǒng)穩(wěn)定判據(jù)以及穩(wěn)定性分類思路.文獻(xiàn)[18]中基于等系統(tǒng)模型分析策略,依據(jù)所提方法進(jìn)行阻尼穩(wěn)定控制器參數(shù)設(shè)計,提出優(yōu)化函數(shù)進(jìn)行參數(shù)計算新思路,在特定工況下進(jìn)行阻尼穩(wěn)定控制器參數(shù)配置.文獻(xiàn)[19]中綜述目前大型風(fēng)力發(fā)電(簡稱風(fēng)電)機(jī)組實現(xiàn)慣量控制,從穩(wěn)定性分析、控制系統(tǒng)優(yōu)化、性能分析與評估等方面,介紹慣量控制在雙饋風(fēng)電機(jī)組以及全功率風(fēng)電機(jī)組的研究與工程現(xiàn)狀.文獻(xiàn)[20]中研究隨機(jī)擾動下電力系統(tǒng)強(qiáng)迫振蕩特點(diǎn),基于同步機(jī)機(jī)械功率上隨機(jī)擾動功率譜密度模型,為多機(jī)系統(tǒng)建立模態(tài)疊加形式頻率響應(yīng)模型,由功率譜密度分析方法得到同步機(jī)功角的功率譜密度表達(dá)式,進(jìn)而分析隨機(jī)擾動下系統(tǒng)強(qiáng)迫振蕩的4個影響因素.

綜上文獻(xiàn)研究發(fā)現(xiàn),當(dāng)前新型電力系統(tǒng)小干擾穩(wěn)定性研究普遍使用基于增廣狀態(tài)矩陣描述線性化模型.此種模型實現(xiàn)較為復(fù)雜,且編程繁瑣,為后續(xù)計算分析帶來極大不便.現(xiàn)有新型電力系統(tǒng)穩(wěn)定性分析方法對應(yīng)對運(yùn)行工作點(diǎn)快速變化尚存在亟需解決問題,對小擾動穩(wěn)定分析提出更高要求.傳統(tǒng)小干擾穩(wěn)定性分析特征值法會產(chǎn)生由于維數(shù)過高而引起的計算負(fù)擔(dān),對電力系統(tǒng)運(yùn)行工作點(diǎn)快速變化情況無法做到快速穩(wěn)定性分析,并且單純求解特征值無法為之后控制提供理論指導(dǎo)依據(jù).為解決上述問題,提出一種基于李雅普諾夫函數(shù)直接分析法的新型電力系統(tǒng)快速穩(wěn)定性研究方法,利用同步相量測量裝置(PMU)廣域監(jiān)測系統(tǒng)數(shù)據(jù),在理論上考慮運(yùn)行工況不確定穩(wěn)定性,為實現(xiàn)小干擾穩(wěn)定性增強(qiáng)控制策略提供基礎(chǔ),可實現(xiàn)對新型電力系統(tǒng)小干擾穩(wěn)定性快速分析,并可在實際應(yīng)用中改善系統(tǒng)穩(wěn)定性.

1 基本理論

1.1 梯度下降法

梯度下降法[21]是經(jīng)典優(yōu)化方法之一,通常用于無約束問題.梯度下降法的依據(jù)是對某個具體函數(shù)而言,沿梯度正方向函數(shù)值增加速度最快,與之相反,函數(shù)值減少最快.

(1)

式中:

(2)

在達(dá)到收斂條件時,迭代結(jié)束.由梯度下降法迭代公式可知,下一個點(diǎn)選擇與當(dāng)前點(diǎn)位置和梯度相關(guān);反之,如果要計算函數(shù)的最大值,沿著梯度反方向前進(jìn)即可.從整體看,無論是計算函數(shù)的最大值或最小值,都需要構(gòu)建一個迭代關(guān)系λ′:

(3)

即對所有的i,都滿足迭代關(guān)系.從而可以得到梯度下降法函數(shù)表達(dá)式:

λ(x)=x-ηf(x)

(4)

1.2 新型電力系統(tǒng)穩(wěn)定性直接法分析模型

對于實際大型電力系統(tǒng)性能分析而言,要建立整個新型電力系統(tǒng)高階數(shù)學(xué)模型難度較大,即便建立出模型,在實際應(yīng)用過程中也會產(chǎn)生諸如維數(shù)爆炸、模型過于復(fù)雜、計算過于繁瑣等問題.為規(guī)避電力系統(tǒng)模型建立過于復(fù)雜,文獻(xiàn)[22]中提出將大型電力系統(tǒng)中每個發(fā)電機(jī)模型假定為一個子系統(tǒng),各子系統(tǒng)通過傳輸網(wǎng)絡(luò)互相連接,在考慮子系統(tǒng)模型及其相互作用條件的情況下,通常采用子系統(tǒng)模型代替整個系統(tǒng)模型進(jìn)行穩(wěn)定性分析.對于整個系統(tǒng)統(tǒng)一的慣性中心,在某一區(qū)域發(fā)生振蕩時,總存在表現(xiàn)為兩組發(fā)電機(jī)之間相對搖擺的減速與加速機(jī)組,上述機(jī)組可以簡化為一個由兩臺等值發(fā)電機(jī)構(gòu)成互聯(lián)的電力系統(tǒng)等效研究[23],等值系統(tǒng)機(jī)間振蕩對應(yīng)原系統(tǒng)的區(qū)間振蕩.當(dāng)發(fā)生低頻振蕩時,兩臺發(fā)電機(jī)轉(zhuǎn)速增量和功角增量具有同頻反向正弦振蕩的特點(diǎn).

首先對于PMU所持續(xù)采集到的原始大數(shù)據(jù),采用文獻(xiàn)[24]中的方法進(jìn)行篩選與降維預(yù)處理,建立數(shù)據(jù)矩陣,利用數(shù)據(jù)降維方法將高維矩陣降階到低維矩陣,代入電力系統(tǒng)狀態(tài)方程中,從而可以利用PMU所采集到的數(shù)據(jù)代入數(shù)據(jù)矩陣中對新型電力系統(tǒng)穩(wěn)定性進(jìn)行分析.基于同步發(fā)電機(jī)模型進(jìn)行推導(dǎo),利用接入雙饋異步風(fēng)力發(fā)電機(jī)組(DFIG)的電力系統(tǒng)狀態(tài)矩陣代替新型電力系統(tǒng)模型進(jìn)行分析,并依據(jù)其已被證明有效性的降階模型進(jìn)行進(jìn)一步穩(wěn)定性分析.

發(fā)電機(jī)勵磁繞組方程:

(5)

定子電壓方程:

ud=Xqiq-raid

(6)

(7)

式中:ud為d軸電壓;uq為q軸電壓;Xq為q軸同步電抗;iq為q軸電流;ra為電樞電阻.

轉(zhuǎn)子運(yùn)動方程:

(9)

式中:TJ為發(fā)電機(jī)組的慣性時間常數(shù);D為轉(zhuǎn)軸定阻尼常數(shù);Tm為原動機(jī)輸出機(jī)械力矩.

進(jìn)一步忽略f繞組的暫態(tài)過程,上述模型可以降階為發(fā)電機(jī)經(jīng)典二階實用模型即負(fù)荷采用恒阻抗模型:

(10)

當(dāng)發(fā)電機(jī)組采用二階經(jīng)典模型時,前文所述等效兩臺等值發(fā)電機(jī)系統(tǒng)轉(zhuǎn)子運(yùn)動方程組為

(11)

式中:Δδ1和Δδ2分別為兩臺發(fā)電機(jī)的功角增量;Δω1和Δω2分別為兩臺發(fā)電機(jī)的轉(zhuǎn)速增量;M1和M2、D1和D2分別為兩臺發(fā)電機(jī)慣性時間常數(shù)和阻尼力矩系數(shù);ΔP1e和ΔP2e為電磁功率增量.

考慮新型電力系統(tǒng)中電力系統(tǒng)核心設(shè)備的變化以及新能源發(fā)電大規(guī)模的應(yīng)用,風(fēng)電已經(jīng)成為現(xiàn)今主流的分布式發(fā)電形式之一,本文應(yīng)用接入DFIG電力系統(tǒng)矩陣進(jìn)行分析,給出新型電力系統(tǒng)狀態(tài)矩陣:

x=[Δδ1Δδ2Δω1Δω2]T

(12)

(13)

式中:x為狀態(tài)偏差列向量;A為狀態(tài)空間矩陣.

(14)

電磁功率方程可表示為

(15)

風(fēng)電場的動態(tài)有功表示為

(16)

式中:k1=(E1Vw0/X13)cos(δ10-δw0)∝1/X13;k2=(-E2Vw0/X23)cos(δw0-δ20)∝1/X23;k3=g1g2∝ΔPw;k4=g1g3∝ΔPw.其中:E1和E2分別為區(qū)域1和區(qū)域2的電動勢;Vw0為風(fēng)電場接入點(diǎn)電壓初始值,δw0為風(fēng)電場接入點(diǎn)相位增量Δδw的初始值;X13和X23分別為兩臺發(fā)電機(jī)到風(fēng)電場接入點(diǎn)的電抗;g1,g2,g3為風(fēng)電場動態(tài)頻率特性模值,各變量下標(biāo)中的0表示初始值.

2 基于梯度動態(tài)偏差李雅普諾夫直接法穩(wěn)定性分析

2.1 李雅普諾夫直接法

李雅普諾夫直接法借助一個李雅普諾夫能量函數(shù)V(x(t))及根據(jù)微分方程所計算得到的V沿著軌跡導(dǎo)數(shù)符號性質(zhì)直接推斷穩(wěn)定性問題[25].考慮以下一般非線性自治動力系統(tǒng):

(17)

系統(tǒng)存在以下狀態(tài)方程式:

(18)

對于線性定常系統(tǒng)而言,找到一個使得狀態(tài)在原點(diǎn)平衡漸近穩(wěn)定充分必要條件為:對于任意給定的一個對稱正定矩陣Q,一定存在唯一一個正定對稱矩陣P,使得原線性定常系統(tǒng)狀態(tài)方程成立.對于其平衡狀態(tài)達(dá)到大范圍漸近穩(wěn)定充分必要條件以代數(shù)方程表示為

ATP+PA=-Q

(19)

并且有

V(x)=xTPx

(20)

V(x)稱作該系統(tǒng)的李雅普諾夫能量函數(shù),被用作系統(tǒng)穩(wěn)定性分析.由于正定實對稱矩陣Q任意給定,為計算簡便,下文計算取Q為單位矩陣.將矩陣A代入系統(tǒng)狀態(tài)式(18),建立李雅普諾夫方程式(19),求解方程式可得對稱矩陣P,通過對稱矩陣P正定性判定可以對系統(tǒng)穩(wěn)定性進(jìn)行判斷.

2.2 新型電力系統(tǒng)小干擾狀態(tài)變量偏差值

利用以下對角矩陣代入求解計算系統(tǒng)某個狀態(tài)變量偏差:

QJ=diag[1,…,1,0,…,0]

(21)

若求解某個狀態(tài)變量的偏差量,取相對應(yīng)的對角矩陣元素為1,其余為0即可.當(dāng)式(21)中描述系統(tǒng)已被判定為穩(wěn)定后,借助李雅普諾夫方程所得矩陣求取狀態(tài)變量積分值[16]如下:

(22)

QJ=diag(q1,q2,…,qn)

(23)

在計算中,用單位矩陣I代入QJ矩陣,可得下式:

(24)

J=xTPJx

(25)

式中:PJ為對稱矩陣,將QJ和A代入式(19)可求得.J值能夠反映系統(tǒng)狀態(tài)偏差量的動態(tài)穩(wěn)定性,J值越小,表明狀態(tài)偏差量衰減速度快且振蕩幅度小,動態(tài)過程短,系統(tǒng)小干擾穩(wěn)定性好[16].

2.3 基于梯度下降法的狀態(tài)響應(yīng)曲線時間加權(quán)動態(tài)偏差值

當(dāng)分別分析每個發(fā)電機(jī)組時,為計算簡便,且起到與上述積分狀態(tài)偏差量近似效果,提出一種基于狀態(tài)量曲線梯度下降偏差的時間加權(quán)值,用于描述單個狀態(tài)量的低頻振蕩過程,對狀態(tài)量的動態(tài)過程進(jìn)行分析評價,且可以對穩(wěn)定過程進(jìn)行分析.首先以振蕩曲線為梯度下降的目標(biāo)函數(shù):

x(t)=(t1,t2,…,tn)

(26)

對該函數(shù)進(jìn)行如下梯度下降迭代過程,當(dāng)i>0時有:

(27)

…

(28)

一旦達(dá)到收斂條件,在梯度為零點(diǎn)時求得該點(diǎn)極值,梯度下降迭代結(jié)束.由梯度下降法迭代公式可知,下一個點(diǎn)選擇和當(dāng)前點(diǎn)位置與其梯度相關(guān).反之,如果計算函數(shù)最大值,沿著梯度反方向前進(jìn)即可.如此反復(fù),即可計算出整條曲線每次振蕩波峰以及波谷極值點(diǎn)數(shù)值.從整體來看,無論是計算函數(shù)最大值或最小值,都需要構(gòu)建如下迭代關(guān)系h:

h(t)=t-ηf(t)

(29)

在經(jīng)過本文直接法判定系統(tǒng)穩(wěn)定性步驟后,如果系統(tǒng)已經(jīng)可被判定為穩(wěn)定,則記錄整個振蕩過程時間,對該狀態(tài)變量曲線使用梯度下降法,迭代計算出該曲線每次振蕩波峰以及波谷極值點(diǎn)數(shù)值,依據(jù)每次振蕩過程時間在整個振蕩過程占比加權(quán)計算出整個振蕩過程極值偏差量,并以該值來具體評估單個機(jī)組單一狀態(tài)量振蕩劇烈程度.設(shè)ti為波峰到波谷之間的時間間隔,αi為波峰/波谷,具體加權(quán)偏差值計算式如下:

(30)

2.4 基于PMU數(shù)據(jù)的梯度動態(tài)偏差李雅普諾夫直接法實現(xiàn)

提出基于梯度動態(tài)偏差李雅普諾夫直接法對新型電力系統(tǒng)穩(wěn)定性分析,具體實現(xiàn)流程圖如圖1所示.算法步驟如下:

圖1 基于梯度動態(tài)偏差李雅普諾夫直接法流程圖

步驟1將PMU數(shù)據(jù)進(jìn)行匯總形成電網(wǎng)數(shù)據(jù),對PMU所實時采集到的數(shù)據(jù)進(jìn)行降維預(yù)處理,生成小干擾穩(wěn)定性分析模型需要的數(shù)據(jù)矩陣.

步驟2采用二階同步發(fā)電機(jī)模型代替整個電力系統(tǒng)模型進(jìn)行小干擾穩(wěn)定性分析,并給出前文已推導(dǎo)接入風(fēng)電機(jī)組新型電力系統(tǒng)小干擾穩(wěn)定性分析所需要的狀態(tài)空間矩陣A.

步驟3將經(jīng)過處理的PMU數(shù)據(jù)帶入狀態(tài)空間矩陣A,利用MATLAB中命令求解李亞普諾夫方程,并利用其求解式(19).在求解方程式時,由于正定實對稱矩陣Q任意給定,為計算簡便,故可以取Q為單位矩陣I,可以快速解出正定對稱矩陣P.通過矩陣P的正定性來判定電力系統(tǒng)穩(wěn)定性,借助李雅普諾夫方程分析系統(tǒng)穩(wěn)定性比較簡捷快速.基于已推導(dǎo)出接入DFIG電力系統(tǒng)狀態(tài)空間矩陣A,采用直接法已求取P矩陣對新型電力系統(tǒng)進(jìn)行小干擾穩(wěn)定性分析.

步驟4借助李雅普諾夫方程求出的矩陣P獲取所求狀態(tài)偏差量平方積分值,實際上是系統(tǒng)各狀態(tài)偏差量平方在正無窮時間軸上的積分之和;在系統(tǒng)被判定穩(wěn)定之后,根據(jù)矩陣對角元素取值不同求得各狀態(tài)偏差量平方積分值加權(quán)和,可以對已被判定穩(wěn)定系統(tǒng)進(jìn)行動態(tài)過程描述.如果僅關(guān)心部分狀態(tài)變量動態(tài)過程,可以利用求取偏差變量平方積分值方法分析該系統(tǒng)達(dá)到穩(wěn)定的動態(tài)過程,且某個單獨(dú)動態(tài)偏差量可以通過對角矩陣代入求解.

步驟5為了單個狀態(tài)量偏差計算方便,且起到與積分狀態(tài)偏差量近似作用,在經(jīng)過直接法判定系統(tǒng)小干擾穩(wěn)定性步驟后,如果系統(tǒng)已被判定為穩(wěn)定,此時記錄下整個振蕩過程時間,對該狀態(tài)變量曲線使用梯度下降法,迭代計算出該曲線每次振蕩波峰以及波谷極值點(diǎn)數(shù)值.依據(jù)式(30),每次振蕩過程時間在整個振蕩過程的占比加權(quán)計算出整個振蕩過程極值偏差量.

3 算例仿真

3.1 小干擾穩(wěn)定性直接法算例分析

將IEEE 10機(jī)39節(jié)點(diǎn)新英格蘭系統(tǒng)[26]接入風(fēng)電場,如圖2所示.圖2中每個節(jié)點(diǎn)處均安裝 PMU,系統(tǒng)電壓等級為 345 kV,頻率為 60 Hz,采樣頻率為 3 kHz,系統(tǒng)中各母線負(fù)荷設(shè)置為恒功率模型.首先,在節(jié)點(diǎn)28和節(jié)點(diǎn)29處加入PMU數(shù)據(jù),并分出一條支路,由多臺DFIG風(fēng)力機(jī)等效為一個風(fēng)電場.然后,借助文獻(xiàn)[22]中實測PMU數(shù)據(jù),選擇其中 0～179.9 s時間段所采集到的電流電壓幅值、相角,采用文獻(xiàn)[24]中方法對數(shù)據(jù)進(jìn)行預(yù)處理,將高維矩陣降維生成穩(wěn)定性分析需要的參數(shù)矩陣,將PMU數(shù)據(jù)代入狀態(tài)空間A矩陣中,運(yùn)用前文所提李雅普諾夫直接法對算例系統(tǒng)進(jìn)行小擾動穩(wěn)定性分析計算.使用MATLAB設(shè)計穩(wěn)定性判定以及計算指標(biāo)函數(shù)值的算法程序.

圖2 風(fēng)電場接入新英格蘭10機(jī)39節(jié)點(diǎn)系統(tǒng)結(jié)構(gòu)

利用2.4節(jié)所述算法步驟分別對WSCC 3機(jī)系統(tǒng)、中國電科院6機(jī)系統(tǒng)以及新英格蘭10機(jī)系統(tǒng)應(yīng)用特征值法和直接法進(jìn)行穩(wěn)定性判定,分析結(jié)果如表1所示.可知,兩種方法呈現(xiàn)判定結(jié)果一致,表明采用李雅普諾夫直接法可以準(zhǔn)確可靠地判定電力系統(tǒng)穩(wěn)定性,且不會因為某些參數(shù)省略導(dǎo)致判定錯誤結(jié)果.

表1 特征值法與直接法穩(wěn)定性判定結(jié)果

3.2 不同風(fēng)電場出力下新型電力系統(tǒng)低頻特性影響

參照文獻(xiàn)[27],設(shè)置區(qū)域1向區(qū)域2送電,風(fēng)電機(jī)組在送電側(cè),輸送距離為25 km,在增大風(fēng)電場輸出功率同時調(diào)整同步發(fā)電機(jī)的輸出功率,保持聯(lián)絡(luò)線傳輸功率為400 MW不變,其中風(fēng)電場輸出功率占區(qū)域1總輸出功率比例用風(fēng)力場出力表示,新型電力系統(tǒng)區(qū)域間振蕩模式變化情況如表2所示.可知,隨著風(fēng)電場出力的增加,區(qū)域間振蕩模式1的振蕩頻率和阻尼比在總體上均表現(xiàn)為增大趨勢;區(qū)域間振蕩模式2的振蕩頻率逐漸增大,但是阻尼比卻逐漸減小.

表2 新型電力系統(tǒng)區(qū)域間振蕩模式隨風(fēng)電場出力變化

3.3 狀態(tài)偏差平方積分值的低頻振蕩算例分析

在電力系統(tǒng)發(fā)生小干擾時,機(jī)電模式最小阻尼比常作為低頻振蕩評價指標(biāo),文獻(xiàn)[16]中按照動態(tài)最小阻尼比配置準(zhǔn)則對電力系統(tǒng)進(jìn)行阻尼穩(wěn)定控制器(PSS)配置,并給出特征值法所得最小阻尼比.本文采用的系統(tǒng)動態(tài)偏差平方積分值法利用系統(tǒng)狀態(tài)矩陣A建立李雅普諾夫方程組,調(diào)用 MATLAB 中李雅普諾夫方程函數(shù)解出矩陣P,判定矩陣P正定性,從而利用直接法判斷系統(tǒng)穩(wěn)定性.根據(jù)前文所述動態(tài)偏差值方法,借助李亞普諾夫方程求得指標(biāo)函數(shù)值即各狀態(tài)偏差量加權(quán)平方積分值,動態(tài)偏差方積分值能夠反映系統(tǒng)狀態(tài)偏差量動態(tài)性能.

表3列出文獻(xiàn)[16]中最小阻尼比以及動態(tài)偏差量平方積分值計算結(jié)果,最小阻尼比是廣為應(yīng)用的低頻振蕩評價指標(biāo),而動態(tài)平方積分值也可以對系統(tǒng)性能進(jìn)行分析.如表3所示,隨著PSS在系統(tǒng)中配置數(shù)量越來越多,最小阻尼比也越來越大,表明系統(tǒng)穩(wěn)定性越來越好.

表3 算例系統(tǒng)穩(wěn)定性分析過程

3.4 依據(jù)狀態(tài)變量曲線時間加權(quán)梯度動態(tài)偏差值PSS配置

利用所得指標(biāo)函數(shù)值對算例系統(tǒng)小擾動穩(wěn)定性進(jìn)行分析,以驗證方法的可行性及正確性.由文獻(xiàn)[27]可知,在新英格蘭10機(jī)39節(jié)點(diǎn)系統(tǒng)中6#機(jī)組為等值機(jī),故不需要設(shè)計PSS.動態(tài)阻尼比配置方法可用于配置PSS位置確定,從而達(dá)到優(yōu)化及電力系統(tǒng)小干擾穩(wěn)定性目的.在MATLAB軟件Simulink中搭建含風(fēng)電新英格蘭10機(jī)39節(jié)點(diǎn)系統(tǒng),在該系統(tǒng)中6#機(jī)組構(gòu)建小干擾信號模塊,模擬系統(tǒng)受到小干擾低頻振蕩,從而達(dá)到對系統(tǒng)進(jìn)行小干擾穩(wěn)定性分析目的.利用所求取狀態(tài)變量曲線時間加權(quán)梯度動態(tài)偏差值描述系統(tǒng)受到小干擾后發(fā)電機(jī)組低頻振蕩動態(tài)過程,通過對系統(tǒng)中不同機(jī)組低頻振蕩動態(tài)偏差值比較,對系統(tǒng)中振蕩幅度大且振蕩時間長機(jī)組進(jìn)行PSS配置,發(fā)電機(jī)未配置PSS功角振蕩情況如圖3所示,各發(fā)電機(jī)組狀態(tài)曲線時間加權(quán)梯度偏差和阻尼比如表4所示.

表4 各發(fā)電機(jī)組狀態(tài)曲線時間加權(quán)梯度偏差和阻尼比

3.5 依據(jù)狀態(tài)變量曲線時間加權(quán)梯度動態(tài)偏差值與動態(tài)最小阻尼比配置方法PSS配置

為進(jìn)一步驗證所提狀態(tài)變量曲線時間加權(quán)梯度動態(tài)偏差值對PSS配置可行性和有效性,針對算例仿真系統(tǒng),按照表4數(shù)據(jù)采用本文方法所得功角動態(tài)偏差平方積分值對系統(tǒng)PSS進(jìn)行配置,配置PSS機(jī)組結(jié)果分別為4#、5#、7#、8#、10#機(jī)組.同時根據(jù)文獻(xiàn)[16],利用動態(tài)最小阻尼比法對發(fā)電機(jī)組配置PSS用于應(yīng)用效果對比,配置PSS機(jī)組結(jié)果分別為1#、2#、4#、7#、9#機(jī)組.因此,在新型電力系統(tǒng)中配置PSS改善系統(tǒng)小干擾穩(wěn)定性,對同一系統(tǒng)采用不同PSS配置方案,可通過振蕩曲線直觀反映振蕩動態(tài)過程以及電力系統(tǒng)小干擾穩(wěn)定性改善性能,已配置PSS發(fā)電機(jī)功角振蕩如圖4所示.

圖4 已配置PSS發(fā)電機(jī)功角振蕩

由圖3與圖4對比可知,與配置前相比,算例系統(tǒng)機(jī)組配置PSS后穩(wěn)定性有所提升,功角波動程度明顯減小,說明依據(jù)動態(tài)阻尼比配置方法以及依據(jù)狀態(tài)曲線時間加權(quán)偏差值配置方法的有效性,但本文提出的功角偏差時間加權(quán)值可反映出所有機(jī)組在時間區(qū)間上振蕩幅度,因此可更準(zhǔn)確全面地反映系統(tǒng)動態(tài)過程.進(jìn)一步對比圖4(a)和圖4(b)可知,狀態(tài)曲線時間加權(quán)梯度偏差PSS配置系統(tǒng)低頻振蕩過程時間顯著縮短,并且振蕩幅度也得到降低,故本文所提依據(jù)狀態(tài)曲線時間加權(quán)梯度偏差PSS配置法在特定情況下具有顯著有效性和可行性.

4 結(jié)語

在“30·60”碳達(dá)峰碳中和背景下,未來電力系統(tǒng)將會發(fā)生結(jié)構(gòu)性變化,研究新型電力系統(tǒng)小干擾穩(wěn)定性機(jī)理意義重大.針對目前新型電力系統(tǒng)中存在影響穩(wěn)定性問題,提出基于PMU數(shù)據(jù)的梯度動態(tài)偏差李雅普諾夫直接分析機(jī)制,在建立新型電力系統(tǒng)穩(wěn)定性直接法分析所需二階模型基礎(chǔ)上,將李雅普諾夫直接法、新型電力系統(tǒng)小干擾狀態(tài)變量偏差值以及基于梯度下降法的狀態(tài)響應(yīng)曲線時間加權(quán)動態(tài)偏差值有機(jī)結(jié)合,完成基于PMU數(shù)據(jù)的梯度動態(tài)偏差李雅普諾夫直接法的實現(xiàn).根據(jù)含風(fēng)電新英格蘭10機(jī)39節(jié)點(diǎn)系統(tǒng),從小干擾穩(wěn)定性直接法算例分析、不同風(fēng)電場出力下新型電力系統(tǒng)低頻特性影響、狀態(tài)偏差平方積分值的低頻振蕩算例分析、依據(jù)狀態(tài)變量曲線時間加權(quán)梯度動態(tài)偏差值PSS配置以及依據(jù)狀態(tài)變量曲線時間加權(quán)梯度動態(tài)偏差值與動態(tài)最小阻尼比配置方法PSS配置5個方面驗證方法的正確性和有效性.結(jié)果表明,本文方法可縮短系統(tǒng)振蕩過程,在PSS配置數(shù)量相對有限情況下能夠提升新型電力系統(tǒng)小干擾穩(wěn)定性.此外,新能源引起發(fā)電波動性以及來自大量電力電子設(shè)備頻繁投切和系統(tǒng)拓?fù)浣Y(jié)構(gòu)改變,導(dǎo)致新型電力系統(tǒng)的運(yùn)行點(diǎn)變化更快速,本文方法對系統(tǒng)動態(tài)過程評估還沒有統(tǒng)一評價指標(biāo),后續(xù)將進(jìn)一步研究新型電力系統(tǒng)小干擾穩(wěn)定方法適應(yīng)性.