胡 光 莊可好 高暉勝 辛煥海 汪海蛟 張雅君

低慣量交流系統(tǒng)并網(wǎng)變流器次/超同步振蕩分析

胡 光1莊可好1高暉勝1辛煥海1汪海蛟2張雅君3

（1. 浙江大學(xué)電氣工程學(xué)院 杭州 310027 2. 中國(guó)電力科學(xué)研究院有限公司 北京 100192 3. 國(guó)網(wǎng)上海市電力公司電力科學(xué)研究院 上海 200437）

低慣量電力系統(tǒng)中跟網(wǎng)型變流器可能會(huì)與同步機(jī)轉(zhuǎn)子動(dòng)態(tài)產(chǎn)生交互，給系統(tǒng)穩(wěn)定性分析引入新的問(wèn)題，僅靠跟網(wǎng)型設(shè)備動(dòng)態(tài)和短路比衡量交流系統(tǒng)強(qiáng)度不再準(zhǔn)確。為此，該文綜合分析系統(tǒng)慣量和短路比對(duì)跟網(wǎng)型變流器小干擾同步穩(wěn)定的影響，提出適用于分析低慣量系統(tǒng)中并網(wǎng)變流器穩(wěn)定性的方法，所提方法計(jì)算量低、物理意義明確。首先，推導(dǎo)變流器和同步機(jī)混聯(lián)系統(tǒng)的閉環(huán)模型，并提取出能表征低慣量動(dòng)態(tài)的相互作用系數(shù)矩陣，定性分析兩類(lèi)設(shè)備間的交互關(guān)系。其次，基于模態(tài)攝動(dòng)理論，將原多機(jī)系統(tǒng)解耦為穩(wěn)定性等效的子系統(tǒng)。然后，選取適用于分析小干擾同步穩(wěn)定問(wèn)題的相位主導(dǎo)回路，根據(jù)該回路提出綜合考慮慣量和短路比的穩(wěn)定性分析方法。最后，基于仿真算例驗(yàn)證了所提低慣量系統(tǒng)中并網(wǎng)變流器小干擾同步穩(wěn)定性分析方法的有效性。

低慣量系統(tǒng) 小干擾同步穩(wěn)定 廣義短路比 最小子系統(tǒng) 系統(tǒng)強(qiáng)度

0 引言

“雙碳”背景下，隨著新能源的大規(guī)模并網(wǎng)，電力系統(tǒng)逐漸演變?yōu)榈蛻T量、低短路比的弱交流系統(tǒng)[1]（以下簡(jiǎn)稱(chēng)為“弱交流系統(tǒng)”）。相較于傳統(tǒng)電力系統(tǒng)，弱交流系統(tǒng)中源網(wǎng)耦合加劇，動(dòng)力學(xué)行為更加復(fù)雜[2-4]，并網(wǎng)變流器振蕩失穩(wěn)問(wèn)題凸顯。

在研究系統(tǒng)振蕩失穩(wěn)問(wèn)題時(shí)，系統(tǒng)電壓支撐強(qiáng)度（以下簡(jiǎn)稱(chēng)為“系統(tǒng)強(qiáng)度”）因可用于評(píng)估系統(tǒng)的穩(wěn)定性和響應(yīng)性能而備受關(guān)注[5]。系統(tǒng)強(qiáng)度由電網(wǎng)強(qiáng)度和設(shè)備動(dòng)態(tài)共同決定，其中，電網(wǎng)強(qiáng)度描述了不考慮接入設(shè)備動(dòng)態(tài)的開(kāi)環(huán)交流電網(wǎng)性能，可通過(guò)短路比（Short Circuit Ratio, SCR）指標(biāo)進(jìn)行量化評(píng)估[6]。對(duì)于設(shè)備動(dòng)態(tài)方面對(duì)系統(tǒng)強(qiáng)度的影響，鑒于目前我國(guó)交流系統(tǒng)中新能源機(jī)組的鎖相同步控制技術(shù)較為成熟，在未來(lái)一段時(shí)間，新能源設(shè)備主要以跟網(wǎng)型控制并網(wǎng)（不特殊說(shuō)明，本文所研究的變流器特指跟網(wǎng)型變流器），組網(wǎng)功能仍主要由同步發(fā)電機(jī)組承擔(dān)[7]，而同步機(jī)的慣量在很大程度上影響著其組網(wǎng)的能力[8]。由此可見(jiàn)，除交流系統(tǒng)短路比外，系統(tǒng)慣量也是衡量系統(tǒng)強(qiáng)度、評(píng)估并網(wǎng)變流器穩(wěn)定性（后文穩(wěn)定性均指小干擾同步穩(wěn)定性）的重要指標(biāo)。

短路比常被用于表征電力電子設(shè)備饋入交流系統(tǒng)時(shí)交流電網(wǎng)的相對(duì)強(qiáng)度，可用于分析交流網(wǎng)絡(luò)連接特征對(duì)系統(tǒng)穩(wěn)定性的影響[9-12]。變流器單饋入系統(tǒng)中，隨著并網(wǎng)線路阻抗增大，短路比減小，鎖相環(huán)動(dòng)態(tài)和線路動(dòng)態(tài)耦合加劇[13]，導(dǎo)致由鎖相環(huán)主導(dǎo)的系統(tǒng)小干擾同步失穩(wěn)問(wèn)題頻發(fā)。文獻(xiàn)[9]提出廣義短路比（generalized Short Circuit Ratio, gSCR）的概念，建立了系統(tǒng)穩(wěn)定性和電網(wǎng)強(qiáng)度之間的顯性關(guān)系，實(shí)現(xiàn)了多變流器并網(wǎng)系統(tǒng)穩(wěn)定裕度的量化。文獻(xiàn)[10]基于模態(tài)攝動(dòng)理論，將廣義短路比理論拓展至異構(gòu)多機(jī)系統(tǒng)。然而，上述研究重點(diǎn)關(guān)注系統(tǒng)低短路比特性對(duì)系統(tǒng)穩(wěn)定性的影響，分析時(shí)通常將同步機(jī)等效為理想電壓源節(jié)點(diǎn)，沒(méi)有考慮弱交流系統(tǒng)的低慣量特性。

次/超同步頻段下，弱交流系統(tǒng)的低慣量特性使得同步機(jī)無(wú)法等效為理想電壓源，可能導(dǎo)致同步機(jī)與變流器兩類(lèi)設(shè)備的同步環(huán)節(jié)產(chǎn)生交互耦合[14-19]，由此可見(jiàn)，交流系統(tǒng)對(duì)并網(wǎng)變流器穩(wěn)定性的支撐強(qiáng)度還受到系統(tǒng)慣量的影響。針對(duì)兩機(jī)互聯(lián)系統(tǒng)，文獻(xiàn)[16-17]通過(guò)時(shí)域仿真和特征值分析，研究了同步機(jī)慣量常數(shù)和網(wǎng)絡(luò)阻抗對(duì)系統(tǒng)穩(wěn)定性的影響；文獻(xiàn)[18]分析了同步機(jī)的物理特性和控制參數(shù)對(duì)電力電子設(shè)備穩(wěn)定性的影響，揭示了兩類(lèi)設(shè)備動(dòng)態(tài)的交互路徑。針對(duì)多機(jī)系統(tǒng)，文獻(xiàn)[19]建立了系統(tǒng)的詳細(xì)狀態(tài)空間模型，通過(guò)仿真分析發(fā)現(xiàn)，與全電力電子設(shè)備電力系統(tǒng)相比，含有小容量同步機(jī)的弱交流系統(tǒng)更容易失穩(wěn)。文獻(xiàn)[20]給出了一種便于拓展的多設(shè)備統(tǒng)一坐標(biāo)系建模方法，保留了各設(shè)備同步環(huán)節(jié)的動(dòng)態(tài)，為探索弱交流系統(tǒng)穩(wěn)定機(jī)理提供了一個(gè)分析方法。文獻(xiàn)[21]針對(duì)同步穩(wěn)定問(wèn)題，兼顧了各類(lèi)設(shè)備模型的魯棒性和穩(wěn)定分析的便捷性，給出了電流源視角下各類(lèi)設(shè)備的統(tǒng)一建模方法。文獻(xiàn)[22]針對(duì)系統(tǒng)小干擾穩(wěn)定性提升問(wèn)題，啟發(fā)式地分析了慣量對(duì)系統(tǒng)穩(wěn)定性的影響。上述文獻(xiàn)基于仿真法或模型解析法探索了低慣量特性對(duì)系統(tǒng)穩(wěn)定性的影響，為分析弱交流系統(tǒng)穩(wěn)定性提供了一定的模型和理論基礎(chǔ)。但基于仿真分析所得結(jié)論可能僅適用于特定工況，不利于揭示失穩(wěn)機(jī)理，而且全系統(tǒng)電磁暫態(tài)模型仿真計(jì)算量過(guò)高，尚缺乏有效的全電磁仿真手段[23]。基于模型的解析方法又存在系統(tǒng)全階模型過(guò)于復(fù)雜，容易出現(xiàn)“維數(shù)災(zāi)”問(wèn)題[24]。此外，上述研究尚未分析弱交流系統(tǒng)多機(jī)間耦合機(jī)理，也缺少綜合考慮弱交流系統(tǒng)低短路比和低慣量特性，以及對(duì)系統(tǒng)穩(wěn)定性的定量化分析。

為此，本文綜合分析了系統(tǒng)廣義短路比和慣量對(duì)跟網(wǎng)型變流器小干擾同步穩(wěn)定的影響，提出了適用于低慣量多機(jī)系統(tǒng)的并網(wǎng)變流器穩(wěn)定裕度量化方法。首先，推導(dǎo)了變流器和同步機(jī)混聯(lián)運(yùn)行系統(tǒng)的閉環(huán)模型，利用相互作用系數(shù)矩陣定性分析了兩類(lèi)設(shè)備間的耦合關(guān)系；其次，進(jìn)一步基于模態(tài)攝動(dòng)理論，將多機(jī)原系統(tǒng)解耦為穩(wěn)定性等效的子系統(tǒng)；然后，根據(jù)適用于分析系統(tǒng)小干擾同步穩(wěn)定問(wèn)題的主導(dǎo)回路，提出了綜合考慮慣量和短路比的穩(wěn)定性分析方法。最后，通過(guò)仿真算例驗(yàn)證了所提穩(wěn)定性分析方法的有效性。

1 同步機(jī)和變流器混聯(lián)系統(tǒng)建模

本節(jié)首先對(duì)要研究的問(wèn)題及使用的假設(shè)條件進(jìn)行了描述。然后以問(wèn)題為導(dǎo)向，建立了同步機(jī)和變流器混聯(lián)系統(tǒng)的閉環(huán)模型。進(jìn)一步通過(guò)數(shù)學(xué)變換，分離出系統(tǒng)中兩類(lèi)設(shè)備耦合的動(dòng)態(tài)，形成單位相互作用系數(shù)矩陣（Unit Interaction Factor Matrix, UIFM），并根據(jù)UIFM定性分析了兩類(lèi)設(shè)備間的 交互。

1.1 問(wèn)題描述及假設(shè)

考慮圖1所示的含臺(tái)變流器和臺(tái)同步機(jī)的多機(jī)混聯(lián)系統(tǒng)，其中節(jié)點(diǎn)1～連接新能源變流器，節(jié)點(diǎn)+1～+連接同步機(jī)。

圖1 多機(jī)系統(tǒng)示意圖

分析并網(wǎng)變流器鎖相環(huán)主導(dǎo)的系統(tǒng)小干擾穩(wěn)定性時(shí)，由于變流器鎖相環(huán)動(dòng)態(tài)與同步機(jī)搖擺方程動(dòng)態(tài)處于不同時(shí)間尺度，可基于奇異攝動(dòng)定理假設(shè)同步機(jī)節(jié)點(diǎn)電壓幅值和相位恒定，即認(rèn)為同步機(jī)節(jié)點(diǎn)為無(wú)窮大節(jié)點(diǎn)[25]。然而低慣量系統(tǒng)中，兩類(lèi)同步環(huán)節(jié)的帶寬可能產(chǎn)生交互耦合，即在鎖相環(huán)的時(shí)間尺度上，同步機(jī)無(wú)法提供理想的電壓支撐能力，不能再簡(jiǎn)單地將同步機(jī)節(jié)點(diǎn)視為無(wú)窮大節(jié)點(diǎn)，需要考慮其動(dòng)態(tài)對(duì)并網(wǎng)變流器穩(wěn)定性的影響。本文聚焦弱交流系統(tǒng)中并網(wǎng)變流器鎖相環(huán)主導(dǎo)的小干擾同步穩(wěn)定問(wèn)題，擬研究問(wèn)題可歸納為：

（1）弱交流系統(tǒng)中的低慣量特性，是如何誘發(fā)并網(wǎng)變流器小干擾失穩(wěn)的？

（2）如何綜合短路比和慣量?jī)蓚€(gè)維度，定量化分析弱交流系統(tǒng)的小干擾穩(wěn)定性、評(píng)估系統(tǒng)強(qiáng)度？

為了更好地闡述和論證上述問(wèn)題，考慮以下合理簡(jiǎn)化：①網(wǎng)絡(luò)中電容動(dòng)態(tài)對(duì)所研究的鎖相環(huán)主導(dǎo)的小擾動(dòng)穩(wěn)定問(wèn)題影響較小，各線路的阻感比均相同[5]；②考慮負(fù)荷為恒電流負(fù)荷，在恒功率或恒阻抗負(fù)荷時(shí)，可以通過(guò)端口等值等方法，將負(fù)荷等效至設(shè)備端口處，進(jìn)而將其動(dòng)態(tài)歸至設(shè)備側(cè)進(jìn)行解 耦[26]；③新能源變流器間的動(dòng)態(tài)特性和同步機(jī)間的動(dòng)態(tài)特性相似（兩類(lèi)設(shè)備各自同構(gòu)），且同類(lèi)設(shè)備間的相位差較小[9]。當(dāng)考慮兩類(lèi)設(shè)備各自都異構(gòu)時(shí)，可以參考文獻(xiàn)[11]的處理方法，本文不再贅述。

1.2 混聯(lián)系統(tǒng)閉環(huán)模型

按圖1中的虛線可將系統(tǒng)模型劃分為網(wǎng)絡(luò)和設(shè)備兩部分，對(duì)應(yīng)的導(dǎo)納矩陣分別為N()和G()。根據(jù)多變量頻域理論，系統(tǒng)閉環(huán)模型sys()表示[27]為

1.2.1 設(shè)備側(cè)導(dǎo)納模型

對(duì)于小干擾同步穩(wěn)定問(wèn)題，設(shè)備以電流源形式建模并兼顧了各類(lèi)設(shè)備模型的魯棒性和穩(wěn)定分析的便捷性[21]，可建立簡(jiǎn)潔統(tǒng)一的閉環(huán)模型。因此，本文設(shè)備側(cè)建模以電流源形式建模。

首先建立變流器設(shè)備線性化的導(dǎo)納模型。全局坐標(biāo)系下，第（=1～）臺(tái)變流器導(dǎo)納模型VSC()可以表示為

其中

式中，yg1()、yg2()、yg3()和yg4()為變流器在同步dq坐標(biāo)系下以自身容量標(biāo)幺后的導(dǎo)納矩陣元素，推導(dǎo)過(guò)程及具體表達(dá)式見(jiàn)文獻(xiàn)[13]；為穩(wěn)態(tài)時(shí)坐標(biāo)系間的相位差；Bi為第臺(tái)變流器經(jīng)系統(tǒng)基準(zhǔn)容量標(biāo)幺后的自身額定容量。值得一提的是，同步dq坐標(biāo)系是穩(wěn)態(tài)時(shí)與局部dq坐標(biāo)系重合的同步旋轉(zhuǎn)坐標(biāo)系，該坐標(biāo)系下的導(dǎo)納模型包含了同步環(huán)節(jié)（鎖相環(huán)）的動(dòng)態(tài)，因此僅通過(guò)穩(wěn)態(tài)的坐標(biāo)變換[21]和容量歸一化，即可得到全局坐標(biāo)系下變流器設(shè)備（=1～）導(dǎo)納矩陣動(dòng)態(tài)。

建立同步機(jī)的導(dǎo)納模型。為了重點(diǎn)分析慣量的耦合作用，本文建立諾頓等效后內(nèi)電流源端口處的模型，將自身容量標(biāo)幺下的定子感應(yīng)電動(dòng)勢(shì)的動(dòng)態(tài)SG歸至網(wǎng)絡(luò)，如附圖1所示。由于同步機(jī)諾頓等效后為內(nèi)電流源并聯(lián)導(dǎo)納的形式，將SG歸至網(wǎng)絡(luò)時(shí)增加了同步機(jī)設(shè)備節(jié)點(diǎn)的自導(dǎo)納，增加量為對(duì)角矩陣SG=diag(Bj/SG)，其中對(duì)角元素Bj/SG為比值的形式，表征了考慮容量歸算后各節(jié)點(diǎn)導(dǎo)納增加量的有名值。式（4）給出了同步全局坐標(biāo)系下第（=+1～+）臺(tái)同步機(jī)內(nèi)電源端口處的導(dǎo)納模型，具體推導(dǎo)過(guò)程及各元素的表達(dá)式如附錄式（A1）～式（A8）所示。

結(jié)合跟網(wǎng)型變流器和同步機(jī)的導(dǎo)納模型，全局坐標(biāo)系下設(shè)備側(cè)動(dòng)態(tài)模型可表示為

其中

1.2.2 交流電網(wǎng)模型

全局坐標(biāo)系下，交流電網(wǎng)動(dòng)態(tài)模型[9]（以電流流出變流器為正方向）為

其中

1.3 基于UIFM的設(shè)備耦合特性分析

系統(tǒng)閉環(huán)框圖如圖2所示，對(duì)應(yīng)的矩陣形式的系統(tǒng)閉環(huán)模型可表示為

求解系統(tǒng)閉環(huán)特征方程det(sys0())=0，即可根據(jù)系統(tǒng)是否有特征值實(shí)部大于0判斷系統(tǒng)穩(wěn)定性。為了進(jìn)一步分析同步機(jī)動(dòng)態(tài)對(duì)并網(wǎng)變流器穩(wěn)定性的影響，根據(jù)網(wǎng)絡(luò)特征保持的Schur分解[28]，可以得到關(guān)注變流器穩(wěn)定問(wèn)題的閉環(huán)模型為

為表達(dá)簡(jiǎn)潔，式（9）及本節(jié)后文各式隱去了“()”。針對(duì)本文關(guān)注的低慣量場(chǎng)景，考慮系統(tǒng)中的同步機(jī)容量較小，根據(jù)Neumann級(jí)數(shù)性質(zhì)，式（9）可進(jìn)一步表示為

其中

由此可見(jiàn)，當(dāng)不考慮同步機(jī)動(dòng)態(tài)的交互耦合時(shí)（對(duì)應(yīng)節(jié)點(diǎn)變?yōu)闊o(wú)窮大或接地節(jié)點(diǎn)），系統(tǒng)的閉環(huán)特征方程中僅包含變流器和網(wǎng)絡(luò)動(dòng)態(tài)的交互，與變流器多饋入系統(tǒng)的閉環(huán)特征方程相同[10]。此時(shí)相互作用矩陣UIF=0×n。由此可見(jiàn)，UIF項(xiàng)包含了同步機(jī)和互聯(lián)兩類(lèi)設(shè)備的網(wǎng)絡(luò)動(dòng)態(tài)，定性地描述了低慣量系統(tǒng)中同步機(jī)和變流器動(dòng)態(tài)的交互作用，為低慣量系統(tǒng)的穩(wěn)定性分析提供了模型基礎(chǔ)。此外，整理為UIF形式的閉環(huán)模型式（10），可以通過(guò)攝動(dòng)理論和模態(tài)解耦進(jìn)行降階，從而能通過(guò)低階等效模型量化分析低慣量對(duì)并網(wǎng)變流器穩(wěn)定性的影響。

2 穩(wěn)定性等效的最小子系統(tǒng)

在第1節(jié)建立的系統(tǒng)閉環(huán)模型的基礎(chǔ)上，本節(jié)首先利用矩陣攝動(dòng)理論，構(gòu)造了一個(gè)與原系統(tǒng)穩(wěn)定性分析等效的同構(gòu)系統(tǒng)；其次通過(guò)模態(tài)解耦，將構(gòu)造出的等效同構(gòu)系統(tǒng)解耦為多個(gè)子系統(tǒng)，并給出了主導(dǎo)系統(tǒng)穩(wěn)定性的子系統(tǒng)（即最弱子系統(tǒng)）的小信號(hào)電路模型。由此可提取出影響原復(fù)雜系統(tǒng)穩(wěn)定性的關(guān)鍵因素，并將其具象為子系統(tǒng)中底層元件的 交互。

2.1 等效同構(gòu)系統(tǒng)

文獻(xiàn)[10]基于廣義短路比的概念，給出了一系列將多機(jī)同構(gòu)的系統(tǒng)解耦為等效子系統(tǒng)的方法，進(jìn)而可以通過(guò)分析低階的子系統(tǒng)，表征原高階系統(tǒng)的穩(wěn)定性。根據(jù)上述思路，首先需要構(gòu)造一個(gè)與原系統(tǒng)穩(wěn)定性等效的同構(gòu)系統(tǒng)。根據(jù)矩陣攝動(dòng)理論[30]，可以構(gòu)造出與式（10）所示原系統(tǒng)穩(wěn)定性等效的同構(gòu)系統(tǒng)為

其中

表1 等效同構(gòu)系統(tǒng)中各參量含義及選取方法

2.2 穩(wěn)定性等效子系統(tǒng)

對(duì)于式（14）中的等效同構(gòu)系統(tǒng)，可將其解耦為多個(gè)獨(dú)立的子系統(tǒng)[10]，重點(diǎn)關(guān)注其中主導(dǎo)系統(tǒng)穩(wěn)定性的最弱子系統(tǒng)，該子系統(tǒng)與原系統(tǒng)穩(wěn)定性等效，其閉環(huán)模型可表示為

式中，DVSC=VSCeq()Dc為等效變流器模型的電流小信號(hào)動(dòng)態(tài)；DSG=SGeq()Dc為等效同步機(jī)模型的電流小信號(hào)動(dòng)態(tài)；Dnet=gSCR()Dc為流過(guò)等效網(wǎng)絡(luò)模型電流的小信號(hào)動(dòng)態(tài)；Dload=LeqDc為等效負(fù)荷模型的電流小信號(hào)動(dòng)態(tài)；Dc為虛擬公共連接點(diǎn)的電壓小信號(hào)模型。根據(jù)基爾霍夫定律，圖3給出了解耦子系統(tǒng)對(duì)應(yīng)的小信號(hào)電路模型。

圖3 等效子系統(tǒng)的小信號(hào)電路模型

3 基于主導(dǎo)回路的穩(wěn)定性分析

本節(jié)首先選取了適用于分析小干擾同步穩(wěn)定問(wèn)題的相位主導(dǎo)回路，得到與最小子系統(tǒng)穩(wěn)定性等效的單輸入單輸出（Single-Input Single-Output, SISO）系統(tǒng)，可以通過(guò)奈奎斯特判據(jù)分析其穩(wěn)定性；然后建立了衡量系統(tǒng)穩(wěn)定裕度的指標(biāo)，并給出了綜合考慮系統(tǒng)慣量和短路比的二維穩(wěn)定運(yùn)行域計(jì)算方法。

3.1 主導(dǎo)回路的選取

第2節(jié)建立的子系統(tǒng)模型考慮了設(shè)備的詳細(xì)動(dòng)態(tài)，其中一些環(huán)節(jié)對(duì)于本文關(guān)注的小干擾穩(wěn)定問(wèn)題影響較小。

首先，為便于揭示耦合機(jī)理并分離出慣量和電網(wǎng)強(qiáng)度對(duì)變流器小干擾穩(wěn)定性的影響，對(duì)兩類(lèi)設(shè)備模型進(jìn)行一些合理簡(jiǎn)化。本文關(guān)注的變流器小干擾同步失穩(wěn)問(wèn)題主要是鎖相環(huán)與弱電網(wǎng)（低短路比低慣量）相互作用導(dǎo)致的次/超同步振蕩。對(duì)于變流器，在該振蕩模態(tài)所在的頻段，根據(jù)奇異攝動(dòng)理論，保留決定穩(wěn)定性的鎖相環(huán)動(dòng)態(tài)進(jìn)行建模，忽略功率外環(huán)等其他環(huán)節(jié)的動(dòng)態(tài)，僅保留式（2）中yg4()的動(dòng)態(tài)[13]；對(duì)于同步機(jī)，為重點(diǎn)關(guān)注慣量的影響，采用經(jīng)典二階模型[20, 31]。結(jié)合附錄式（A7）的具體表達(dá)式，保留了兩類(lèi)設(shè)備基本同步動(dòng)態(tài)性質(zhì)的子系統(tǒng)簡(jiǎn)化模型可表示為

選取式（21）中簡(jiǎn)化系統(tǒng)的主導(dǎo)回路，對(duì)于某個(gè)環(huán)節(jié)主導(dǎo)的穩(wěn)定問(wèn)題，基于該環(huán)節(jié)對(duì)應(yīng)的主導(dǎo)回路模型得到的穩(wěn)定性分析結(jié)果的魯棒性和標(biāo)稱(chēng)性較好[32]。本文側(cè)重于分析跟網(wǎng)型變流器鎖相環(huán)動(dòng)態(tài)與弱交流系統(tǒng)耦合導(dǎo)致的小干擾同步穩(wěn)定問(wèn)題，變流器模型中的yg4()元素包含了鎖相環(huán)動(dòng)態(tài)，反映了變流器的同步動(dòng)態(tài)，因此選取yg4()同步回路作為主導(dǎo)回路進(jìn)行分析，對(duì)應(yīng)SISO系統(tǒng)閉環(huán)特征方 程為

式中，()、()分別為系統(tǒng)的前向通路和反饋回路。以變流器端口的視角觀察，反饋回路()包含了除變流器外系統(tǒng)其余環(huán)節(jié)的動(dòng)態(tài)，反映了變流器所接入系統(tǒng)的系統(tǒng)強(qiáng)度。由于()包含了系統(tǒng)的慣量和廣義短路比，因此系統(tǒng)強(qiáng)度受到等效慣量和廣義短路比的共同影響，其中廣義短路比反映了系統(tǒng)的電網(wǎng)強(qiáng)度，慣量項(xiàng)反映了系統(tǒng)的慣量水平。

至此，本節(jié)通過(guò)數(shù)學(xué)變換將式（8）所示的多輸入多輸出（Multiple-Input Multiple-Output, MIMO）系統(tǒng)，轉(zhuǎn)化為穩(wěn)定性等效的SISO 系統(tǒng)，進(jìn)而可以通過(guò)奈奎斯特穩(wěn)定判據(jù)判斷系統(tǒng)的穩(wěn)定性并量化穩(wěn)定裕度。

3.2 穩(wěn)定性分析

根據(jù)上述分析，圖4給出了多機(jī)系統(tǒng)的穩(wěn)定性分析及穩(wěn)定裕度評(píng)估流程。具體地，首先，計(jì)算系統(tǒng)的穩(wěn)定運(yùn)行域，由于系統(tǒng)的穩(wěn)定性受到gSCR和eq兩個(gè)指標(biāo)影響，在合理范圍內(nèi)遍歷兩指標(biāo)的參數(shù)值，在不同參數(shù)取值的情況下對(duì)式（22）應(yīng)用奈奎斯特判據(jù)判斷子系統(tǒng)是否穩(wěn)定，進(jìn)而形成兩個(gè)指標(biāo)共同約束的二維穩(wěn)定運(yùn)行域（Feasible Region, FR）；其次，對(duì)于某運(yùn)行工況下的多機(jī)系統(tǒng)，通過(guò)圖4所示方法得到該工況下的等效子系統(tǒng)，當(dāng)該子系統(tǒng)的gSCR和eq落在可行域中時(shí)，即(gSCR,eq)∈FR，系統(tǒng)穩(wěn)定；最后，在系統(tǒng)穩(wěn)定時(shí)，可以根據(jù)穩(wěn)定裕度函數(shù)s評(píng)估系統(tǒng)強(qiáng)度。

4 算例分析

為驗(yàn)證上述理論分析的有效性，基于Matlab/ Simulink平臺(tái)搭建如附圖2所示的兩區(qū)4機(jī)系統(tǒng)和附圖3所示的10機(jī)39節(jié)點(diǎn)系統(tǒng)，選取系統(tǒng)基準(zhǔn)容量為100 MV·A，設(shè)備關(guān)鍵控制參數(shù)、設(shè)備容量和交流網(wǎng)絡(luò)線路參數(shù)見(jiàn)附表1～附表5，同步發(fā)電機(jī)的勵(lì)磁等控制環(huán)節(jié)所用參數(shù)見(jiàn)文獻(xiàn)[31]。

圖4 多機(jī)系統(tǒng)穩(wěn)定性分析流程

4.1 頻域理論分析

本節(jié)基于兩區(qū)4機(jī)系統(tǒng)驗(yàn)證所提理論分析方法的有效性。

首先，驗(yàn)證第2節(jié)中解耦子系統(tǒng)與原系統(tǒng)小干擾穩(wěn)定的等效性。本文主要關(guān)注不同類(lèi)型設(shè)備間的交互，為簡(jiǎn)化分析，負(fù)荷采用恒電流負(fù)荷。修改系統(tǒng)網(wǎng)絡(luò)阻抗和慣量的值，觀察原兩區(qū)4機(jī)系統(tǒng)和解耦子系統(tǒng)的主導(dǎo)特征根的變化，其中原系統(tǒng)的特征根通過(guò)Matlab/Simulink中搭建模型的狀態(tài)空間方程計(jì)算得到，子系統(tǒng)的特征根通過(guò)Matlab中的m文件編程計(jì)算得到。兩區(qū)4機(jī)系統(tǒng)等效前后主導(dǎo)特征根對(duì)比如圖5所示。

圖5 兩區(qū)4機(jī)系統(tǒng)等效前后主導(dǎo)特征根對(duì)比

圖5a給出了eq=20的高慣量系統(tǒng)中，網(wǎng)絡(luò)線路電感36、46和56由各自的初始值逐漸增加至初始值的600%，對(duì)應(yīng)廣義短路比gSCR由2.876減小至1.913時(shí)，原兩區(qū)4機(jī)系統(tǒng)和解耦子系統(tǒng)的主導(dǎo)特征根的變化，可以看出，電網(wǎng)強(qiáng)度降低，惡化了跟網(wǎng)型變流器小干擾同步穩(wěn)定性。圖5b給出了gSCR= 2.1的低短路比系統(tǒng)中，各同步機(jī)慣量由5降至1.9時(shí)，系統(tǒng)主導(dǎo)特征根的變化，隨著系統(tǒng)等效慣量降低，圖5b中同步機(jī)轉(zhuǎn)子運(yùn)動(dòng)方程對(duì)應(yīng)模態(tài)的振蕩頻率逐漸增高，增至8 Hz附近（對(duì)應(yīng)特征根虛部為50 rad/s）時(shí)，兩類(lèi)設(shè)備同步環(huán)節(jié)動(dòng)態(tài)的帶寬產(chǎn)生交互耦合，導(dǎo)致系統(tǒng)發(fā)生小干擾失穩(wěn)。由此可見(jiàn)，低慣量系統(tǒng)的相位動(dòng)態(tài)特性，將會(huì)增加鎖相環(huán)型變流器的失穩(wěn)風(fēng)險(xiǎn)。在系統(tǒng)強(qiáng)度（網(wǎng)絡(luò)阻抗和慣量）變化的過(guò)程中，原始兩區(qū)4機(jī)系統(tǒng)和等效子系統(tǒng)的主導(dǎo)特征根基本保持一致，驗(yàn)證了子系統(tǒng)和原兩區(qū)4機(jī)系統(tǒng)小擾動(dòng)穩(wěn)定的等效性，因此通過(guò)模態(tài)解耦得到的子系統(tǒng)可以用于分析弱交流系統(tǒng)的穩(wěn)定性問(wèn)題。

值得一提的是，上述算例驗(yàn)證了子系統(tǒng)的穩(wěn)定性與原系統(tǒng)等效，在該方法正確的情況下，可以利用解耦子系統(tǒng)分析各類(lèi)穩(wěn)定性問(wèn)題，例如，同步機(jī)調(diào)速器與電力電子設(shè)備外環(huán)的低頻段交互[18]、電力電子設(shè)備間的交互[33]等。由于本文關(guān)注低慣量系統(tǒng)中同步機(jī)轉(zhuǎn)子運(yùn)動(dòng)動(dòng)態(tài)和跟網(wǎng)型變流器鎖相環(huán)動(dòng)態(tài)的耦合，對(duì)于其他參數(shù)的影響不再展開(kāi)分析。此外，本文所提簡(jiǎn)化模型式（21）的合理性和準(zhǔn)確性驗(yàn)證詳見(jiàn)附錄第3節(jié)。

接下來(lái)，驗(yàn)證第3節(jié)中選取的同步主導(dǎo)回路分析低慣量系統(tǒng)穩(wěn)定問(wèn)題的有效性。在m文件中編寫(xiě)式（22）、式（23）所示SISO系統(tǒng)的閉環(huán)模型，并繪制其開(kāi)環(huán)傳遞函數(shù)的奈奎斯特曲線。圖6a和圖6b分別給出了gSCR由2.876減小至1.913和慣量由5降至1.9時(shí)的奈奎斯特曲線。可以看出，主導(dǎo)回路的穩(wěn)定性判斷結(jié)果與系統(tǒng)極點(diǎn)分析所得結(jié)果（見(jiàn)圖5）一致。進(jìn)一步地，根據(jù)3.2節(jié)中式（24）計(jì)算圖6中穩(wěn)定系統(tǒng)的裕度，表2給出了系統(tǒng)強(qiáng)度量化結(jié)果，s的值越小，穩(wěn)定裕度越大，系統(tǒng)強(qiáng)度越強(qiáng)。

圖6 兩區(qū)4機(jī)系統(tǒng)主導(dǎo)回路的奈奎斯特曲線

表2 兩區(qū)4機(jī)系統(tǒng)強(qiáng)度量化結(jié)果

Tab.2 Quantitative results of the two area four machine power system strength

4.2 時(shí)域仿真驗(yàn)證

進(jìn)一步地，從時(shí)域角度驗(yàn)證所提理論的準(zhǔn)確性。首先遍歷系統(tǒng)等效慣量由1.5～9變化和系統(tǒng)廣義短路比由1.9～3.1變化的參數(shù)域，根據(jù)如圖4所示的穩(wěn)定性分析流程，計(jì)算各工況下系統(tǒng)是否穩(wěn)定，并給出了如圖7所示的系統(tǒng)二維穩(wěn)定運(yùn)行域。

圖7 兩區(qū)4機(jī)系統(tǒng)二維穩(wěn)定運(yùn)行域

在Matlab/Simulink中搭建附圖2中的兩區(qū)4機(jī)系統(tǒng)的電磁暫態(tài)仿真模型。系統(tǒng)的初始工況為圖7中的A點(diǎn)，在1.5 s和3 s分別切換至B和C點(diǎn)。系統(tǒng)運(yùn)行在平衡點(diǎn)的情況下，于0.25 s時(shí)在同步機(jī)1的轉(zhuǎn)子處施加一轉(zhuǎn)速小擾動(dòng)時(shí)，變流器1端口的有功功率時(shí)域仿真波形如圖8所示。

圖8 兩區(qū)4機(jī)系統(tǒng)時(shí)域仿真波形

可以看出，在A、B和C三點(diǎn)處對(duì)應(yīng)的工況下，時(shí)域波形分別振蕩發(fā)散、等幅振蕩和振蕩收斂，與理論計(jì)算結(jié)果一致。其中，振蕩發(fā)散的特征根為0.042±44.33i，振蕩收斂特征根為-0.091±44.33i，振蕩收斂對(duì)應(yīng)特征根實(shí)部的絕對(duì)值大于振蕩發(fā)散對(duì)應(yīng)特征根的絕對(duì)值，因此收斂速度快于發(fā)散的速度，與圖8中的仿真結(jié)果一致。此外，從曲線中得到的原兩區(qū)4機(jī)系統(tǒng)的振蕩頻率為7.12 Hz，這與圖5b中根據(jù)子系統(tǒng)計(jì)算得到的臨界穩(wěn)定特征根0.009± 44.33i的振蕩頻率7.05 Hz基本一致，說(shuō)明了原始兩區(qū)4機(jī)系統(tǒng)和等效子系統(tǒng)的最弱振蕩模式相同。

4.3 適用性驗(yàn)證

本節(jié)基于IEEE 10機(jī)39節(jié)點(diǎn)系統(tǒng)，旨在驗(yàn)證本文所提穩(wěn)定性分析方法的適用性，因此重點(diǎn)對(duì)比了最終得到的SISO簡(jiǎn)化模型和原系統(tǒng)的穩(wěn)定性等效性。具體地，在模態(tài)攝動(dòng)解耦和主導(dǎo)回路選取后，得到穩(wěn)定性等效的SISO系統(tǒng)，并對(duì)比了Nyquist穩(wěn)定性判斷結(jié)果和時(shí)域仿真結(jié)果。

首先，根據(jù)附表5中系統(tǒng)交流網(wǎng)絡(luò)線路參數(shù)，計(jì)算出系統(tǒng)的廣義短路比gSCR=2.18。根據(jù)圖4中的穩(wěn)定性分析流程，將原多機(jī)系統(tǒng)模態(tài)解耦為等效子系統(tǒng)并選取同步主導(dǎo)回路進(jìn)行分析。圖9給出了保持網(wǎng)絡(luò)參數(shù)不變（gSCR=2.18）情況下，等效慣量eq為2、2.82和4時(shí)的系統(tǒng)Nyquist曲線。由圖可知，在三種工況下，基于主導(dǎo)回路的穩(wěn)定性判斷結(jié)果分別為振蕩發(fā)散、等幅振蕩和振蕩收斂。

圖9 10機(jī)39節(jié)點(diǎn)系統(tǒng)主導(dǎo)回路的奈奎斯特曲線

在系統(tǒng)位于初始工況（gSCR=2.18,eq=2）平衡點(diǎn)的情況下，于1 s時(shí)在設(shè)備1（變流器）的端口處施加一q軸電壓小擾動(dòng)，并在8 s和11 s時(shí)將等效慣量設(shè)置為eq=2.82和eq=4，圖10給出了變流器5端口的有功功率時(shí)域仿真波形。由圖可知，時(shí)域仿真結(jié)果與理論分析結(jié)果一致，這驗(yàn)證了本文所提穩(wěn)定性分析方法在多機(jī)系統(tǒng)中的適用性。

圖10 10機(jī)39節(jié)點(diǎn)系統(tǒng)時(shí)域仿真波形

綜上所述，本節(jié)從模態(tài)理論分析和時(shí)域仿真兩方面驗(yàn)證了所提模態(tài)攝動(dòng)解耦和基于主導(dǎo)回路的穩(wěn)定性分析方法的有效性。

5 結(jié)論

本文針對(duì)弱交流系統(tǒng)中并網(wǎng)變流器次/超同步振蕩問(wèn)題，綜合分析了廣義短路比和等效慣量對(duì)并網(wǎng)變流器穩(wěn)定性的影響。主要結(jié)論如下：

1）同步機(jī)容量小的低慣量系統(tǒng)可以解耦為個(gè)單變流器和單同步機(jī)互聯(lián)的最小子系統(tǒng)，其中最弱子系統(tǒng)與原系統(tǒng)的穩(wěn)定性等效。

2）弱交流系統(tǒng)中，在設(shè)備動(dòng)態(tài)固定的前提下，并網(wǎng)變流器小干擾穩(wěn)定性由線路連接阻抗和同步機(jī)慣量時(shí)間常數(shù)共同決定，可由廣義短路比和系統(tǒng)等效慣量共同約束的二維穩(wěn)定域刻畫(huà)，當(dāng)廣義短路比或等效慣量減小時(shí)，系統(tǒng)穩(wěn)定性惡化；反之，則改善。

3）隨著同步機(jī)容量減小，系統(tǒng)等效慣量水平逐漸降低，兩類(lèi)設(shè)備同步環(huán)節(jié)動(dòng)態(tài)的帶寬產(chǎn)生交互耦合，可能誘發(fā)并網(wǎng)變流器鎖相環(huán)主導(dǎo)的次/超同步振蕩失穩(wěn)。

如何綜合分析不同環(huán)節(jié)交互下的系統(tǒng)穩(wěn)定性問(wèn)題，是未來(lái)主要的工作之一。

1. 同步機(jī)建模及其簡(jiǎn)化模型

同步發(fā)電機(jī)定子感應(yīng)電動(dòng)勢(shì)動(dòng)態(tài)可以通過(guò)等效電感SG的動(dòng)態(tài)表示為

內(nèi)電動(dòng)勢(shì)的幅值動(dòng)態(tài)為

內(nèi)電動(dòng)勢(shì)的相位動(dòng)態(tài)為

式中，、分別為同步機(jī)慣量和阻尼系數(shù)；、分別為轉(zhuǎn)子角與轉(zhuǎn)子轉(zhuǎn)速；0為系統(tǒng)標(biāo)稱(chēng)頻率；e為有功功率；m為調(diào)速系統(tǒng)機(jī)械功率；gov()為調(diào)速器傳遞函數(shù)，可以由一階慣性環(huán)節(jié)近似表示；gov為其增益；gov為其等效時(shí)間常數(shù)。

根據(jù)式（A4）對(duì)同步機(jī)諾頓等效，附圖1給出了等效前后的同步機(jī)模型電路。將虛線右側(cè)的定子電感動(dòng)態(tài)歸至網(wǎng)絡(luò)側(cè)。

附圖1 同步機(jī)內(nèi)部諾頓等效示意圖

App.Fig.1 Norton equivalence diagram of a synchronous generator at internal port

線性化式（A1）～式（A4）并聯(lián)立求解，得到內(nèi)電流源處的導(dǎo)納模型為

其中

dq同步坐標(biāo)系和全局同步坐標(biāo)系存在穩(wěn)態(tài)相位差時(shí)，僅需進(jìn)行式（A7）坐標(biāo)變換即可。

其中

第3節(jié)中僅關(guān)注相位動(dòng)態(tài)的簡(jiǎn)化模型為

2. 仿真系統(tǒng)介紹

附表1 設(shè)備主要參數(shù)

App.Tab.1 Main parameters of devices

設(shè)備類(lèi)型參 數(shù)數(shù) 值 變流器功率外環(huán)PI參數(shù)0.5, 40 電流內(nèi)環(huán)PI參數(shù)0.5, 40 鎖相環(huán)PI參數(shù)4, 2 100 同步機(jī)慣量和阻尼8, 0 定子等效電感(pu)0.1

附圖2 兩區(qū)4機(jī)系統(tǒng)示意圖

App.Fig.2 Topology diagram of a two area four machine power system

附表2 兩區(qū)4機(jī)系統(tǒng)設(shè)備容量

App.Tab.2 Device capacity in two area four machine power system

設(shè) 備容量(pu) 變流器11.75 變流器21.65 同步機(jī)10.25 同步機(jī)20.35

附表3 兩區(qū)4機(jī)系統(tǒng)交流網(wǎng)絡(luò)線路參數(shù)

App.Tab.3 Network data of two area four machine power system

線路編號(hào)數(shù)值(pu) L150.013 8 L250.014 7 L360.018 3 L460.022 0 L560.024 4 線路阻感比0.1

附圖3 IEEE 10機(jī)39節(jié)點(diǎn)系統(tǒng)拓?fù)?/p>

App.Fig.3 Topology diagram of IEEE 10-machine 39-bus power system

附表4 10機(jī)39節(jié)點(diǎn)系統(tǒng)設(shè)備容量

App.Tab.4 Capacity of device in 10-machine 39-bus power system

設(shè) 備容量(pu) 變流器11.2 變流器21.4 同步機(jī)30.6 變流器41.1 變流器51.1 同步機(jī)60.7 變流器71.5 同步機(jī)80.6 變流器91.6 變流器101.1

3. 簡(jiǎn)化模型驗(yàn)證

首先，驗(yàn)證所提簡(jiǎn)化模型的合理性，式（21）中的簡(jiǎn)化模型忽略了對(duì)所關(guān)注的主導(dǎo)特征根影響較小的環(huán)節(jié)。各動(dòng)態(tài)環(huán)節(jié)對(duì)系統(tǒng)小干擾穩(wěn)定性的影響，可以通過(guò)主導(dǎo)特征根的參與因子刻畫(huà)。具體見(jiàn)附表6，附表6中各環(huán)節(jié)參與因子為同類(lèi)型設(shè)備對(duì)應(yīng)環(huán)節(jié)參與因子求和得到。由附表6可知，鎖相環(huán)動(dòng)態(tài)、轉(zhuǎn)子運(yùn)動(dòng)動(dòng)態(tài)和網(wǎng)絡(luò)動(dòng)態(tài)主導(dǎo)了低慣量系統(tǒng)的小干擾穩(wěn)定性，因此可以根據(jù)3.1節(jié)中分析對(duì)詳細(xì)模型進(jìn)行合理簡(jiǎn)化。

附表5 10機(jī)39節(jié)點(diǎn)系統(tǒng)交流網(wǎng)絡(luò)線路參數(shù)

App.Tab.5 Network data of 10-machine 39-bus power system

線路編號(hào)數(shù)值(pu)線路編號(hào)數(shù)值(pu) L1,20.0370L15,160.008 5 L1,390.022 5L16,170.008 0 L2,30.013 6L16,190.017 5 L2,250.007 7L16,210.012 1 L2,300.016 3L16,240.005 3 L3,40.019 2L17,180.007 4 L3,180.012 0L17,270.015 6 L4,50.011 5L19,330.012 8 L4,140.011 6L19,200.012 4 L5,80.010 1L20,340.016 2 L6,50.002 3L21,220.012 6 L6,70.008 3L22,230.008 6 L6,110.007 4L22,350.030 9 L7,80.004 1L23,240.031 5 L8,90.032 7L23,360.024 5 L9,390.022 5L25,260.029 1 L10,110.003 9L25,370.020 9 L10,130.003 9L26,270.013 2 L10,320.036 0L26,280.042 7 L12,110.039 1L26,290.056 3 L12,130.039 1L28,290.013 6 L13,140.009 1L29,380.032 0 L14,150.019 5L31,60.022 5 線路阻感比0.1——

附表6 主導(dǎo)模態(tài)參與因子

App.Tab.6 Participating factors of the dominant modes

動(dòng)態(tài)環(huán)節(jié)鎖相環(huán)模態(tài)參與因子搖擺方程模態(tài)參與因子 變流器鎖相環(huán)1.050.73 變流器外環(huán)0.0275.7×10-3 變流器電流內(nèi)環(huán)2.60×10-31.16×10-5 同步機(jī)轉(zhuǎn)子動(dòng)態(tài)0.611.03 同步機(jī)勵(lì)磁4.2×10-30.011 同步機(jī)調(diào)速器3.7×10-30.031 網(wǎng)絡(luò)動(dòng)態(tài)0.210.14

驗(yàn)證簡(jiǎn)化模型的準(zhǔn)確性。附圖4進(jìn)一步給出了根據(jù)所建立的簡(jiǎn)化模型和詳細(xì)模型得到的系統(tǒng)主導(dǎo)特征根對(duì)比，其中簡(jiǎn)化模型主導(dǎo)特征根是通過(guò)編寫(xiě)m文件程序計(jì)算得到，詳細(xì)模型主導(dǎo)特征根通過(guò)Simulink線性化得到?？梢钥闯?，簡(jiǎn)化模型與詳細(xì)模型的主導(dǎo)特征根基本一致，綜合上述分析，驗(yàn)證了本文所提簡(jiǎn)化模型的合理性和準(zhǔn)確性。

附圖4 簡(jiǎn)化模型與詳細(xì)模型主導(dǎo)極點(diǎn)對(duì)比

App.Fig.4 Dominant pole comparison diagram of simplified model and detailed model

[1] 袁小明, 張美清, 遲永寧, 等. 電力電子化電力系統(tǒng)動(dòng)態(tài)問(wèn)題的基本挑戰(zhàn)和技術(shù)路線[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2022, 42(5): 1904-1917.

Yuan Xiaoming, Zhang Meiqing, Chi Yongning, et al. Basic challenges of and technical roadmap to power- electronized power system dynamics issues[J]. Proceedings of the CSEE, 2022, 42(5): 1904-1917.

[2] Hatziargyriou N, Milanovi? J, Rahmann C, et al. Stability definitions and characterization of dynamic behavior in systems with high penetration of power electronic interfaced technologies[R]. New York, USA: Institute of Electrical and Electronics Engineers, 2020.

[3] Peng Qiao, Jiang Qin, Yang Yongheng, et al. On the stability of power electronics-dominated systems: challenges and potential solutions[J]. IEEE Transa- ctions on Industry Applications, 2019, 55(6): 7657- 7670.

[4] 邵冰冰, 趙崢, 肖琪, 等. 多直驅(qū)風(fēng)機(jī)經(jīng)柔直并網(wǎng)系統(tǒng)相近次同步振蕩模式參與因子的弱魯棒性分析[J]. 電工技術(shù)學(xué)報(bào), 2023, 38(3): 754-769.

Shao Bingbing, Zhao Zheng, Xiao Qi, et al. Weak robustness analysis of close subsynchronous oscillation modes’ participation factors in multiple direct-drive wind turbines with the VSC-HVDC system[J]. Transactions of China Electrotechnical Society, 2023, 38(3): 754-769.

[5] 周瑀涵, 辛煥海, 鞠平. 基于廣義短路比的多饋入系統(tǒng)強(qiáng)度量化原理與方法: 回顧、探討與展望[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2023, 43(10): 3794-3811.

Zhou Yuhan, Xin Huanhai, Ju Ping. System strength quantification principle and method of multi-infeed systems based on generalized short-circuit ratio: reviews, discussions and outlooks[J]. Proceedings of the CSEE, 2023, 43(10): 3794-3811.

[6] 于琳, 孫華東, 趙兵, 等. 新能源并網(wǎng)系統(tǒng)短路比指標(biāo)分析及臨界短路比計(jì)算方法[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2022, 42(3): 919-929.

Yu Lin, Sun Huadong, Zhao Bing, et al. Short circuit ratio index analysis and critical short circuit ratio calculation of renewable energy grid-connected system[J]. Proceedings of the CSEE, 2022, 42(3): 919-929.

[7] 郭劍波. 構(gòu)建新型電力系統(tǒng)是實(shí)現(xiàn)能源轉(zhuǎn)型、達(dá)成“雙碳”目標(biāo)的有效途徑[N]. 國(guó)家電網(wǎng)報(bào), 2021- 09-07(5).

[8] 齊磊, 趙巍, 孫孝峰, 等. 低調(diào)頻容量高比例新能源主導(dǎo)局部電網(wǎng)低頻失穩(wěn)分析[J]. 電力系統(tǒng)自動(dòng)化, 2022, 46(12): 174-183.

Qi Lei, Zhao Wei, Sun Xiaofeng, et al. Analysis on low-frequency instability of local power grid with low-frequency regulation capacity dominated by high proportion of renewable energy[J]. Automation of Electric Power Systems, 2022, 46(12): 174-183.

[9] 辛煥海, 董煒, 袁小明, 等. 電力電子多饋入電力系統(tǒng)的廣義短路比[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2016, 36(22): 6013-6027.

Xin Huanhai, Dong Wei, Yuan Xiaoming, et al. Generalized short circuit ratio for multi power electronic based devices infeed to power systems[J]. Proceedings of the CSEE, 2016, 36(22): 6013-6027.

[10] 辛煥海, 甘德強(qiáng), 鞠平. 多饋入電力系統(tǒng)廣義短路比: 多樣化新能源場(chǎng)景[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2020, 40(17): 5516-5527.

Xin Huanhai, Gan Deqiang, Ju Ping. Generalized short circuit ratio of power systems with multiple power electronic devices: analysis for various renewable power generations[J]. Proceedings of the CSEE, 2020, 40(17): 5516-5527.

[11] 卓振宇, 張寧, 謝小榮, 等. 高比例可再生能源電力系統(tǒng)關(guān)鍵技術(shù)及發(fā)展挑戰(zhàn)[J]. 電力系統(tǒng)自動(dòng)化, 2021, 45(9): 171-191.

Zhuo Zhenyu, Zhang Ning, Xie Xiaorong, et al. Key technologies and developing challenges of power system with high proportion of renewable energy[J]. Automation of Electric Power Systems, 2021, 45(9): 171-191.

[12] 馬富藝龍, 辛煥海, 劉晨曦, 等. 新能源基地柔性直流送出系統(tǒng)小擾動(dòng)電壓支撐強(qiáng)度評(píng)估[J]. 電工技術(shù)學(xué)報(bào), 2023, 38(21): 5758-5770, 5938.

Ma Fuyilong, Xin Huanhai, Liu Chenxi, et al. Small-disturbance system voltage support strength assessment method for renewables VSC-HVDC delivery system[J]. Transactions of China Electro- technical Society, 2023, 38(21): 5758-5770, 5938.

[13] 辛煥海, 李子恒, 董煒, 等. 三相變流器并網(wǎng)系統(tǒng)的廣義阻抗及穩(wěn)定判據(jù)[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2017, 37(5): 1277-1293.

Xin Huanhai, Li Ziheng, Dong Wei, et al. Generalized-impedance and stability criterion for grid-connected converters[J]. Proceedings of the CSEE, 2017, 37(5): 1277-1293.

[14] Pulgar-Painemal H, Wang Yajun, Silva-Saravia H. On inertia distribution, inter-area oscillations and location of electronically-interfaced resources[J]. IEEE Transactions on Power Systems, 2018, 33(1): 995-1003.

[15] Dong Dong, Wen Bo, Boroyevich D, et al. Analysis of phase-locked loop low-frequency stability in three- phase grid-connected power converters considering impedance interactions[J]. IEEE Transactions on Industrial Electronics, 2015, 62(1): 310-321.

[16] Ding Lizhi, Lu Xiaonan, Tan Jin. Small-signal stability analysis of low-inertia power grids with inverter-based resources and synchronous con- densers[C]//2022 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT), New Orleans, LA, USA, 2022: 1-5.

[17] Ding Lizhi, Lu Xiaonan, Tan Jin. Comparative small-signal stability analysis of grid-forming and grid-following inverters in low-inertia power systems[C]//IECON 2021-47th Annual Conference of the IEEE Industrial Electronics Society, Toronto, ON, Canada, 2021: 1-6.

[18] 徐李清, 郭春義, 楊碩. 聯(lián)接低慣量交流系統(tǒng)的MMC與發(fā)電機(jī)之間的低頻交互振蕩模式研究[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2023, 43(9): 3402-3415.

Xu Liqing, Guo Chunyi, Yang Shuo. Research on low frequency interactive oscillation mode between generator and MMC connected to low inertia AC system[J]. Proceedings of the CSEE, 2023, 43(9): 3402-3415.

[19] Markovic U, Stanojev O, Aristidou P, et al. Understanding small-signal stability of low-inertia systems[J]. IEEE Transactions on Power Systems, 2021, 36(5): 3997-4017.

[20] Gu Yunjie, Li Yitong, Zhu Yue, et al. Impedance- based whole-system modeling for a composite grid via embedding of frame dynamics[J]. IEEE Transactions on Power Systems, 2021, 36(1): 336- 345.

[21] 莊可好, 辛煥海, 高暉勝, 等. 雙軸勵(lì)磁電流源同步機(jī)視角下跟網(wǎng)/構(gòu)網(wǎng)設(shè)備建模及其互聯(lián)系統(tǒng)同步穩(wěn)定性分析[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2023, 43(20): 7759-7773.

Zhuang Kehao, Xin Huanhai, Gao Huisheng, et al. Modeling of grid-forming and grid-following devices from the perspective of dual excitation current sources synchronous generator and synchronous stability analysis of interconnected systems[J]. Proceedings of the CSEE, 2023, 43(20): 7759-7773.

[22] 花赟玥, 楊超然, 何國(guó)慶, 等. 考慮小干擾穩(wěn)定和頻率穩(wěn)定的虛擬慣量配置分析[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2021, 61(5): 437-445.

Hua Yunyue, Yang Chaoran, He Guoqing, et al. Virtual inertia configuration analysis considering small-signal stability and frequency stability[J]. Journal of Tsinghua University (Science and Technology), 2021, 61(5): 437-445.

[23] 陳國(guó)平, 李明節(jié), 董昱, 等. 構(gòu)建新型電力系統(tǒng)仿真體系研究[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2023, 43(17): 6535-6551.

Chen Guoping, Li Mingjie, Dong Yu, et al. Research on the simulation technology architecture for the new-type power system[J]. Proceedings of the CSEE, 2023, 43(17): 6535-6551.

[24] Du Wenjuan, Dong Wenkai, Wang Haifeng. A method of reduced-order modal computation for planning grid connection of a large-scale wind farm[J]. IEEE Transactions on Sustainable Energy, 2020, 11(3): 1185-1198.

[25] Yang Chaoran, Huang Linbin, Xin Huanhai, et al. Placing grid-forming converters to enhance small signal stability of PLL-integrated power systems[J]. IEEE Transactions on Power Systems, 2021, 36(4): 3563-3573.

[26] Tu Jing, Zhou Ming, Cui Hantao, et al. An equivalent aggregated model of large-scale flexible loads for load scheduling[J]. IEEE Access, 2019, 7: 143431- 143444.

[27] Skogestad S, Postelethwaite I. Multivariable feedback control[M]. New York: Wiley Publishing, 1996.

[28] 袁輝, 辛煥海, 王冠中, 等. 含SVG的新能源多饋入系統(tǒng)振蕩分析和廣義短路比計(jì)算[J]. 電力系統(tǒng)自動(dòng)化, 2021, 45(14): 38-46.

Yuan Hui, Xin Huanhai, Wang Guanzhong, et al. Analysis on oscillation of multi-infeed system with renewable energy and static var generator and calculation of its generalized short-circuit ratio[J]. Automation of Electric Power Systems, 2021, 45(14): 38-46.

[29] Gao Huisheng, Wang Guanzhong, Li Zhiyi, et al. Assessing interaction strength between synchronous generator and LCC-HVDC in AC/DC systems[C]// 2020 IEEE Power & Energy Society General Meeting (PESGM), Montreal, QC, Canada, 2020: 1-5.

[30] Stewart G W, Sun Jiguang. Matrix perturbation theory[M]. Boston: Academic Press, 1990.

[31] 高暉勝. 新能源電力系統(tǒng)模態(tài)頻率: 定義、分析與控制[D]. 杭州: 浙江大學(xué), 2022.

Gao Huisheng. Modal frequency in renewable- integrated power system: definition, analysis and control[D]. Hangzhou: Zhejiang University, 2022.

[32] 宮澤旭, 艾力西爾·亞爾買(mǎi)買(mǎi)提, 辛煥海, 等. 新能源電力系統(tǒng)并網(wǎng)設(shè)備小擾動(dòng)穩(wěn)定分析(一): 機(jī)理模型與穩(wěn)定判據(jù)適用性[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2022, 42(12): 4405-4419.

Gong Zexu, Yaermaimaiti Ailixier, Xin Huanhai, et al. Small signal stability analysis of equipment in renewable energy power system (part I): mechanism model and adaptation of stability criterion[J]. Proceedings of the CSEE, 2022, 42(12): 4405-4419.

[33] 高本鋒, 王義, 范輝, 等. 基于阻尼路徑的新能源經(jīng)LCC-HVDC送出系統(tǒng)次同步交互作用分析方法[J]. 電工技術(shù)學(xué)報(bào), 2023, 38(20): 5572-5589.

Gao Benfeng, Wang Yi, Fan Hui, et al. A sub- synchronous interaction analysis method of renewable energy generations integrated with LCC-HVDC system based on damping path[J]. Transactions of China Electrotechnical Society, 2023, 38(20): 5572- 5589.

Sub/Super Synchronous Oscillation Analysis of Grid-Connected Converter in Low Inertia AC System

111123

（1. College of Electrical Engineering Zhejiang University Hangzhou 310027 China 2. China Electric Power Research Institute Beijing 100192 China 3. Electric Power Research Institute State Grid Shanghai Municipal Electric Power Company Shanghai 200437 China）

With the widespread integration of renewable energy sources, the power system is transitioning into a weak AC system characterized by low inertia and a reduced short-circuit ratio. Compared to conventional power systems, the interaction between sources and grids becomes more intricate in weak AC systems, resulting in complex dynamic behaviors. The challenge of oscillation instability in grid-connected inverters is accentuated in low-inertia power systems. In such systems, grid-following inverters may interact with the dynamics of synchronous machine rotors, posing new challenges for stability analysis. Relying solely on the dynamics and short-circuit ratio of grid-following devices is insufficient to assess the strength of the AC system. This paper addresses the small-signal oscillation instability in the sub-synchronous frequency range. The impact of the system's generalized short-circuit ratio and inertia on the small-disturbance synchronous stability of grid-following inverters are analyzed. The proposed stability margin quantification method is designed for grid-connected inverters in low-inertia multi-machine systems with low computational complexity and clear physical significance.

This paper initially establishes a closed-loop model for the hybrid system of synchronous machines and inverters. The dynamic coupling of the two types of devices is separated by mathematically transforming the closed-loop model. A unit interaction factor matrix (UIFM) is generated, including the dynamics of synchronous machines and the interconnected network of the two types of devices. This matrix qualitatively describes the interaction between the dynamics of synchronous machines and inverters in low-inertia systems, providing a model foundation for stability analysis.

The matrix perturbation theory is employed for the closed-loop model, constructing an isomorphic system equivalent to the original system's stability. Modal decoupling is then applied to break down the constructed equivalent isomorphic system into multiple subsystems. The small-signal circuit model of the dominant subsystem is provided. Key factors affecting the stability of the original complex system are extracted, concretizing the interaction of bottom-layer components in the subsystem.

Since the minimal subsystem remains a two-dimensional multiple-input, multiple-output (MIMO) system, this paper selects a phase-dominant loop suitable for analyzing small-disturbance synchronous stability problems. The result in a single-input, single-output (SISO) system is equivalent to the stability of the minimal subsystem. The stability of this system is analyzed using the Nyquist criterion, and stability margin indicators are established. A two-dimensional stable operating region calculation method, considering both system inertia and short-circuit ratio, is provided. Finally, using the Matlab/Simulink platform, the effectiveness of the proposed analysis method is validated using a two-area, four-machine system and a 10-machine 39-bus system.

Low inertia system, small-disturbance synchronization stability, generalized short-circuit ratio, minimum subsystem, system strength

國(guó)家電網(wǎng)有限公司科技資助項(xiàng)目（52094023001W）。

2023-09-01

2023-12-05

10.19595/j.cnki.1000-6753.tces.231434

TM712

胡 光 男，1999年生，博士研究生，研究方向?yàn)樾履茉床⒕W(wǎng)穩(wěn)定分析。E-mail: 12110056@zju.edu.cn

辛煥海 男，1981年生，教授，博士生導(dǎo)師，研究方向?yàn)榻恢绷飨到y(tǒng)穩(wěn)定分析與控制、新能源并網(wǎng)穩(wěn)定分析與控制。E-mail: xinhh@zju.edu.cn（通信作者）

（編輯 陳 誠(chéng)）