孫 杰，尹太元，王 躍，劉 普，王曉雷，樂(lè) 波

（1.西安交通大學(xué)電氣工程學(xué)院，西安 710049；2.鄭州輕工業(yè)大學(xué)電氣與信息工程學(xué)院，鄭州 450000；3.中原工學(xué)院電子信息學(xué)院，鄭州 450007；4.國(guó)網(wǎng)經(jīng)濟(jì)技術(shù)研究院有限公司，北京 102209）

模塊化多電平換流器MMC（modular multi-level converter）由于其模塊化、開關(guān)頻率低、輸出電壓波形畸變小等優(yōu)點(diǎn)，逐漸成為柔性直流輸電系統(tǒng)中廣泛采用的換流器[1-3]。近年來(lái)，分布式能源大規(guī)模并網(wǎng)，電網(wǎng)的慣性和調(diào)節(jié)能力隨著分布式能源的高滲透率而降低[4]，虛擬同步發(fā)電機(jī)VSG（virtual syn-chronous generator）控制技術(shù)由于具有類似同步發(fā)電機(jī)的調(diào)頻和調(diào)壓能力，并能夠?yàn)橄到y(tǒng)提供一定的慣性和阻尼支撐等優(yōu)點(diǎn)，為解決新能源大規(guī)模并網(wǎng)所帶來(lái)的穩(wěn)定性問(wèn)題提供了新的解決方案，逐漸得到國(guó)內(nèi)外學(xué)者的高度重視和廣泛研究[5-6]。

基于端口阻抗的穩(wěn)定性分析法是目前被學(xué)術(shù)界和工業(yè)界廣泛采用的穩(wěn)定性分析方法，多個(gè)柔直工程均采用了此方法，并且其準(zhǔn)確性和有效性也得到了充分的驗(yàn)證[7-9]。目前采用VSG 控制策略的并網(wǎng)逆變器的阻抗建模多集中在兩電平換流器，文獻(xiàn)[10-11]研究了電壓型VSG 和電流型VSG 的小信號(hào)阻抗建模及其與弱電網(wǎng)的交互穩(wěn)定性。而MMC 由于內(nèi)部動(dòng)態(tài)特性復(fù)雜，在子模塊電容電壓和橋臂電流中含有很多諧波，且不同諧波之間存在復(fù)雜的耦合關(guān)系[12-15]，相較于VSC 的阻抗建模要復(fù)雜得多，且其阻抗模型跟VSC 有較大差別。文獻(xiàn)[16]首次基于諧波線性化的建模方法建立了包含環(huán)流抑制和鎖相環(huán)的單電流環(huán)MMC 整流器的阻抗模型。文獻(xiàn)[17]基于諧波狀態(tài)空間法建立了MMC 交流側(cè)阻抗，分析了級(jí)聯(lián)系統(tǒng)的穩(wěn)定性。文獻(xiàn)[18]對(duì)具有完整控制環(huán)路的MMC-HVDC 系統(tǒng)進(jìn)行了序阻抗建模并進(jìn)行了穩(wěn)定性分析。

但是上述文獻(xiàn)均未涉及VSG 控制策略下的MMC 阻抗建模，其端口阻抗特性尚不明確，也鮮有基于阻抗分析法對(duì)采用VSG 控制策略的MMC 進(jìn)行并網(wǎng)穩(wěn)定性研究，而由于VSG 控制策略的MMC在接入弱電網(wǎng)時(shí)對(duì)電網(wǎng)具有慣性和阻尼支撐的優(yōu)點(diǎn)，未來(lái)在柔直領(lǐng)域中有廣闊的應(yīng)用前景，因而十分有必要建立其阻抗模型，并研究其與弱電網(wǎng)交互時(shí)的穩(wěn)定性問(wèn)題，這對(duì)基于MMC 拓?fù)湓诓煌刂颇Ｊ较碌淖杩鼓Ｐ秃驮谌嶂惫こ讨械姆€(wěn)定性研究是一項(xiàng)重要補(bǔ)充，能夠?yàn)槠洳⒕W(wǎng)安全運(yùn)行提供理論依據(jù)，并對(duì)其參數(shù)設(shè)計(jì)提供理論指導(dǎo)。

本文首先研究MMC 的VSG 控制策略，然后基于諧波線性化方法詳細(xì)推導(dǎo)VSG 控制策略下的MMC 阻抗模型，并進(jìn)行驗(yàn)證。然后根據(jù)得到的MMC 阻抗伯德圖，分析其阻抗特性，研究VSG 控制參數(shù)對(duì)MMC 阻抗特性的影響以及VSG 控制策略下MMC 的阻抗靈敏度。再結(jié)合電網(wǎng)阻抗特性，并根據(jù)奈奎斯特穩(wěn)定判據(jù)分析在弱電網(wǎng)條件下VSG控制策略下MMC 的并網(wǎng)穩(wěn)定性。最后通過(guò)實(shí)驗(yàn)驗(yàn)證上述穩(wěn)定性分析的正確性，本文研究結(jié)果可為采用VSG 控制策略的MMC 并網(wǎng)系統(tǒng)的控制和穩(wěn)定性分析以及穩(wěn)定性優(yōu)化提供參考。

1 MMC 拓?fù)浼捌銿SG 控制策略

MMC 的結(jié)構(gòu)框圖如圖1 所示。其由6 個(gè)三相橋臂組成，每相上下2 個(gè)橋臂完全相同，每個(gè)橋臂由N 個(gè)半橋子模塊和1 個(gè)橋臂電感L 組成，L 的等效寄生電感為rL。圖1 中，vj（j=a，b，c）為MMC 交流端口電壓，ij（j=a，b，c）為MMC 交流端口電流，iju（j=a，b，c）為流過(guò)三相上橋臂的電流，ijl（j=a，b，c）為流過(guò)三相下橋臂的電流，Udc為MMC 直流側(cè)電壓。通過(guò)控制橋臂子模塊投入的數(shù)量改變交流端口電壓，橋臂電感L 可以抑制相間環(huán)流并減少內(nèi)部電流的諧波含量。

圖1 MMC 的結(jié)構(gòu)框圖Fig.1 Structural diagram of MMC

圖2 為包含環(huán)流抑制環(huán)節(jié)的VSG 控制框圖，包括功率計(jì)算部分、功率控制部分和環(huán)流抑制部分。功率計(jì)算部分及功率控制部分主要是為了模擬同步發(fā)電機(jī)的特性，其又分為有功-頻率控制和無(wú)功-電壓控制。有功-頻率控制主要是為了模擬同步發(fā)電機(jī)的一次調(diào)頻特性和轉(zhuǎn)子慣性；無(wú)功-電壓控制是為了模擬電壓調(diào)節(jié)特性。

圖2 VSG 控制框圖Fig.2 Block diagram of VSG control

VSG 輸出有功功率Pe和輸出無(wú)功功率Q 由瞬時(shí)功率理論計(jì)算，其計(jì)算公式為

式中：vα、vβ為由三相輸出電壓經(jīng)Clark 變換得到的輸出電壓；iα、iβ為由三相輸出電流經(jīng)Clark 變換得到的輸出電流。

VSG 有功控制器的數(shù)學(xué)方程為

式中：J 為虛擬轉(zhuǎn)動(dòng)慣量；ω0為電網(wǎng)的額定角頻率；ω 為VSG 有功控制環(huán)路的輸出角頻率；Pm為虛擬機(jī)械功率；Pe為電磁功率；D 為阻尼系數(shù)；Pref為有功功率指令；kp為有功-頻率下垂系數(shù)；θ 為VSG 輸出角度。

VSG 無(wú)功控制器的數(shù)學(xué)方程為

式中：E 為虛擬內(nèi)電勢(shì)的峰值；Eref為內(nèi)電勢(shì)峰值指令；Qref為無(wú)功功率指令；Q 為輸出無(wú)功功率；kq為無(wú)功-電壓下垂系數(shù)。

VSG 的有功-頻率控制器和無(wú)功-電壓控制器的輸出共同決定其產(chǎn)生的調(diào)制波，數(shù)學(xué)方程為

式中，vma、vmb、vmc為VSG 生成的調(diào)制波。

環(huán)流抑制控制器生成的調(diào)制波為mcir，將其與vmabc疊加，即可得到MMC 的調(diào)制波信號(hào)m。

2 VSG 控制策略的MMC 阻抗建模

諧波線性化是一種獲取換流器阻抗的簡(jiǎn)便方法，其應(yīng)用范圍比較廣。相比于傳統(tǒng)dq 坐標(biāo)系中的建模，諧波線性化的穩(wěn)態(tài)運(yùn)行軌跡可以包含任意數(shù)量的諧波，計(jì)算出來(lái)的阻抗有清晰的物理解釋并可以直接測(cè)量[19]。下面基于諧波線性化思想對(duì)基于VSG控制策略下的MMC 進(jìn)行阻抗建模。

2.1 MMC 電氣系統(tǒng)小信號(hào)模型

基于諧波線性化的阻抗建模首先需要在MMC交流側(cè)注入頻率為fp的擾動(dòng)電壓，以a 相為例，加入擾動(dòng)電壓之后的a 相端口電壓va（t）為

式中：V0為a 相電網(wǎng)電壓的幅值；f1為a 相電網(wǎng)電壓的頻率；Vp為在MMC 交流端口注入的擾動(dòng)電壓幅值；φvp為其初相位，對(duì)于b 相和c 相，只須根據(jù)注入擾動(dòng)信號(hào)的正負(fù)序?qū)⑾辔辉赼 相相位的基礎(chǔ)上超前或滯后2π/3 即可。

MMC 由于子模塊數(shù)目很多，等效開關(guān)頻率極高，可以認(rèn)為所有子模塊電容電壓始終保持均衡[20]，根據(jù)MMC 的對(duì)稱性，以a 相上橋臂為例，可以得到MMC 的功率級(jí)小信號(hào)模型[16]為

對(duì)于給定頻率為fp的電壓擾動(dòng)將會(huì)導(dǎo)致系統(tǒng)產(chǎn)生一系列頻率為fp±kf1（k=1，2，…，n）的小信號(hào)諧波分量，雖然理論上會(huì)產(chǎn)生無(wú)限多頻次的諧波，但這些諧波的幅值會(huì)隨著k 的增加而減小[21-22]，考慮數(shù)學(xué)運(yùn)算的復(fù)雜性和阻抗建模精度的要求，通常只須考慮3 次以內(nèi)的諧波即可,表1 給出了各頻次諧波正負(fù)序和差模共模的信息[16]。

表1 小信號(hào)諧波之間的關(guān)系Tab.1 Relationship among small-signal harmonics

2.2 MMC 控制部分的小信號(hào)模型

根據(jù)式（6）和式（7）可以看出，Iau、uau和Mau是穩(wěn)態(tài)量，可以從MMC 輸出功率計(jì)算得到，可以消掉，此時(shí)如果能用和表示，那么就可以通過(guò)計(jì)算和的比值進(jìn)而得到MMC 的阻抗。下面對(duì)MMC 交流端口注入正序小擾動(dòng)電壓時(shí)的阻抗進(jìn)行詳細(xì)推導(dǎo)，MMC 交流端口注入負(fù)序電壓擾動(dòng)時(shí)的推導(dǎo)過(guò)程與正序時(shí)類似，本文不再贅述。

在MMC 交流端口注入fp頻次的擾動(dòng)電壓，所產(chǎn)生的小信號(hào)電流響應(yīng)中只有差模量會(huì)流到交流端口，根據(jù)表1 中的正負(fù)序和差模共模關(guān)系，對(duì)于正序電壓擾動(dòng)而言，其產(chǎn)生的差模量的頻次分別為fp和fp-2f1頻次，分別為差模正序信號(hào)和差模負(fù)序信號(hào)，其三相表達(dá)式為

對(duì)式（8）和式（9）進(jìn)行Clark 變換，將小擾動(dòng)電壓和2 個(gè)頻次的小信號(hào)電流轉(zhuǎn)換到αβ 坐標(biāo)系下，Clark變換矩陣Tαβ和變換后小信號(hào)電壓、電流分別為

設(shè)V0和I0為穩(wěn)態(tài)電壓和電流的幅值，Vα和Vβ為αβ 坐標(biāo)系下MMC 的穩(wěn)態(tài)電壓，Iα和Iβ為αβ 坐標(biāo)系下MMC 的穩(wěn)態(tài)電流，表達(dá)式分別為

對(duì)于VSG 功率控制環(huán)路，無(wú)功控制環(huán)的輸出電壓可視為定值[23]，根據(jù)式（1）可以得到小信號(hào)形式的有功功率為

由式（16）可以看出，小信號(hào)有功功率中包含小信號(hào)電壓和2 個(gè)頻次小信號(hào)電流的影響，因此小信號(hào)有功功率可以分為3 部分，其中，對(duì)于fp頻次電流部分為

對(duì)于fp-2f1頻次電流部分為

對(duì)于fp頻次電壓部分為

由圖2 所示的VSG 控制框圖可得VSG 控制部分的傳遞函數(shù)M（s）為

M（s）代表了從電磁功率Pe經(jīng)有功-頻率控制環(huán)節(jié)得到頻率后，再積分成角度θ 的傳遞函數(shù)。

燃?xì)馄髽I(yè)因?yàn)槭艿阶陨硖匦缘挠绊?，在開展財(cái)務(wù)信息化建設(shè)工作時(shí)，根據(jù)需要可以將不同地點(diǎn)、區(qū)域具有共性較強(qiáng)的、重復(fù)率高的、易于標(biāo)準(zhǔn)化和流程化的財(cái)務(wù)業(yè)務(wù)或職能從分、子公司中剝離出來(lái)，之后將其融合到信息平臺(tái)中，實(shí)現(xiàn)信息傳遞和共享。只有這樣，才能將規(guī)模效應(yīng)以及協(xié)同效應(yīng)充分發(fā)揮，減少財(cái)務(wù)運(yùn)行成本，提升財(cái)務(wù)運(yùn)營(yíng)效率以及服務(wù)質(zhì)量，從而降低財(cái)務(wù)運(yùn)行風(fēng)險(xiǎn)，給燃?xì)馄髽I(yè)的健康發(fā)展?fàn)I造良好的條件。

VSG 控制環(huán)路輸出的小信號(hào)相角為

此時(shí)在VSG 輸出的三相調(diào)制波相角中引入了相角擾動(dòng)Δθ。假設(shè)沒(méi)有擾動(dòng)時(shí)的相角為θ0，則有

將式（21）、式（22）代入式（4）中并忽略大信號(hào)項(xiàng)，可以得到VSG 控制環(huán)路生成的小信號(hào)調(diào)制波為

綜上并利用三角函數(shù)積化和差化簡(jiǎn)，由fp頻次電流所影響產(chǎn)生的a 相小信號(hào)調(diào)制波為

結(jié)合式（8）和式（24）可以看出，fp頻次的小信號(hào)電流在調(diào)制波中產(chǎn)生fp和fp-2f1兩個(gè)頻次的諧波，且與原來(lái)的小信號(hào)電流相比均有π/2 的相位延遲。因此，fp頻次的小信號(hào)電流對(duì)VSG 控制環(huán)路生成調(diào)制波的系數(shù)矩陣表示為

同理，fp-2f1頻次的小信號(hào)電流對(duì)VSG 控制環(huán)路生成調(diào)制波的系數(shù)矩陣表示為

其中

因此，VSG 控制環(huán)路中小信號(hào)電流對(duì)調(diào)制波的系數(shù)矩陣表示為

fp頻次的小信號(hào)電壓對(duì)VSG 控制回路生成調(diào)制波的系數(shù)矩陣表示為

環(huán)流抑制環(huán)對(duì)調(diào)制波的系數(shù)矩陣可以表示為[10]

式中，Hc為環(huán)流抑制環(huán)的傳遞函數(shù)。

2.3 MMC 阻抗模型的仿真驗(yàn)證

各控制環(huán)路調(diào)制波的影響可以表示為

將式（34）代入式（6）和式（7），可得到小擾動(dòng)電壓與橋臂小信號(hào)電流的關(guān)系為

式中，U 為7 階單位矩陣。

由此可得小擾動(dòng)電壓與橋臂小信號(hào)電流的7階導(dǎo)納矩陣為

由于fp頻次的相電流是橋臂電流的2 倍[16]，因此MMC 在fp頻次的阻抗模型為

為了驗(yàn)證上述建立阻抗模型的正確性，在Matlab/Simulink 中搭建仿真模型，模型參數(shù)如表2所示，通過(guò)頻率掃描法獲取仿真阻抗幅頻特性和相頻特性并與阻抗模型的解析值進(jìn)行對(duì)比驗(yàn)證。

表2 MMC 的系統(tǒng)和控制參數(shù)Tab.2 System and control parameters of MMC

MMC 阻抗模型掃頻值與解析值對(duì)比結(jié)果如圖3 所示。由圖3 可見(jiàn)，解析值與掃頻值吻合，驗(yàn)證了上述所建立阻抗模型的正確性。

圖3 MMC 阻抗掃頻值與解析值對(duì)比結(jié)果Fig.3 Results of comparison between frequency scanning and analytical values of MMC impedance

3 MMC 并網(wǎng)穩(wěn)定性分析

3.1 VSG 控制參數(shù)對(duì)MMC 阻抗的影響

3.1.1 控制參數(shù)對(duì)阻抗曲線的影響

有功-頻率下垂系數(shù)kp對(duì)MMC 阻抗的影響如圖4 所示，對(duì)比3 組參數(shù)（kp=1，10，100）得到的阻抗曲線，kp增大會(huì)使50 Hz 附近的2 個(gè)阻抗尖峰變小，而其他頻段的阻抗特性幾乎沒(méi)有變化，說(shuō)明kp的改變對(duì)其他頻段的阻抗特性影響不大。

圖4 kp=1、10、100 時(shí)的MMC 阻抗Fig.4 Impedance of MMC when kp equals 1,10 and 100,respectively

阻尼系數(shù)D 對(duì)MMC 阻抗的影響如圖5 所示，從3 組不同的參數(shù)（D=200，300，400）得到的阻抗曲線可以看出，阻尼系數(shù)對(duì)MMC 阻抗的影響主要在45～55 Hz，阻尼系數(shù)越大，50 Hz 附近的2 個(gè)尖峰越小，系統(tǒng)阻尼越強(qiáng)。而其他頻段的阻抗特性幾乎沒(méi)有變化，說(shuō)明D 的改變對(duì)其他頻段的阻抗特性影響不大。

圖5 D=200、300、400 時(shí)的MMC 阻抗Fig.5 Impedance of MMC when D equals 200,300 and 400,respectively

虛擬轉(zhuǎn)動(dòng)慣量J 對(duì)MMC 阻抗的影響如圖6 所示，從3 組不同的參數(shù)（J=0.025 4,0.050 9,0.102 0）得到的阻抗曲線可以看出，虛擬轉(zhuǎn)動(dòng)慣量對(duì)MMC阻抗的影響主要在40～55 Hz。而其他頻段的阻抗特性幾乎沒(méi)有變化，說(shuō)明J 的改變對(duì)其他頻段的阻抗特性影響不大。

圖6 J=0.025 4、0.050 9、0.102 0 時(shí)的MMC 阻抗Fig.6 Impedance of MMC when J equals 0.025 4,0.050 9 and 0.102 0,respectively

3.1.2 控制參數(shù)的阻抗靈敏度分析

控制參數(shù)的阻抗靈敏度能夠反映當(dāng)控制器不同控制參數(shù)變化時(shí)對(duì)換流器端口阻抗特性的影響大小，當(dāng)一個(gè)控制參數(shù)在某個(gè)頻段得到的靈敏度值越大，說(shuō)明該控制參數(shù)在該頻段對(duì)換流器阻抗特性的影響越大，也即阻抗特性在此頻段對(duì)該參數(shù)的變化越敏感；反之，當(dāng)一個(gè)控制參數(shù)在某個(gè)頻段得到的靈敏度值越小時(shí)，說(shuō)明該控制參數(shù)在該頻段對(duì)換流器阻抗特性的影響越小，也即阻抗特性在此頻段對(duì)該參數(shù)的變化越不敏感?？刂茀?shù)的阻抗靈敏度包括阻抗幅值靈敏度和阻抗相角靈敏度，定義為

式中：p 為原本控制參數(shù)；Δp 為控制參數(shù)變化量；S為控制參數(shù)阻抗靈敏度。

有功-頻率下垂系數(shù)kp的阻抗靈敏度如圖7 所示，可見(jiàn)，在頻率為47、49、51 Hz 時(shí)幅值靈敏度曲線較高，說(shuō)明在這3 個(gè)頻率附近的阻抗幅值對(duì)kp的變化比較敏感，而1～36 Hz、60～1 000 Hz 頻段的阻抗幅值靈敏度曲線幾乎與橫軸重合，說(shuō)明在這些頻段kp的變化不會(huì)使阻抗幅值產(chǎn)生明顯變化；在頻率為42、48、52 Hz 時(shí)相角靈敏度較高，說(shuō)明在這3 個(gè)頻率附近的阻抗相角對(duì)kp的變化比較敏感；而1～40 Hz、54～1 000 Hz 頻段的阻抗相角靈敏度曲線幾乎與橫軸重合，說(shuō)明在這些頻段kp的變化不會(huì)使阻抗相角產(chǎn)生明顯變化。上述分析與kp變化時(shí)阻抗曲線發(fā)生變化的頻段一致。

圖7 kp 的阻抗靈敏度Fig.7 Impedance sensitivity of kp

阻尼系數(shù)D 的阻抗靈敏度如圖8 所示，在頻率為47、49、51 Hz 時(shí)幅值靈敏度曲線較高，說(shuō)明在這3 個(gè)頻率附近的阻抗幅值對(duì)D 的變化比較敏感；而1～36 Hz、60～1 000 Hz 頻段的阻抗幅值靈敏度曲線幾乎與橫軸重合，說(shuō)明在這些頻段D 變化來(lái)使阻抗幅值明顯變化；相角靈敏度在頻率為42、48、52 Hz時(shí)較高，說(shuō)明在這3 個(gè)頻率附近的阻抗相角對(duì)D 變化比較敏感；而1～40 Hz、54～1 000 Hz 頻段的阻抗相角靈敏度曲線幾乎與橫軸重合，說(shuō)明在這些頻段D 變化來(lái)使阻抗相角明顯變化。上述分析與D 變化時(shí)阻抗曲線發(fā)生變化的頻段一致。

圖8 D 的阻抗靈敏度Fig.8 Impedance sensitivity of D

虛擬轉(zhuǎn)動(dòng)慣量J 的阻抗靈敏度如圖9 所示，幅值靈敏度曲線在頻率為43、48、52 Hz 時(shí)較高，說(shuō)明在這3 個(gè)頻率附近的阻抗幅值對(duì)J 的變化比較敏感；而1～30 Hz、70～1 000 Hz 頻段的阻抗幅值靈敏度曲線幾乎與橫軸重合，說(shuō)明在這些頻段J 的變化不會(huì)使阻抗幅值產(chǎn)生明顯變化；相角靈敏度在頻率為41、47 Hz 時(shí)較高，說(shuō)明在這2 個(gè)頻率附近的阻抗相角對(duì)J 的變化比較敏感；而1～30 Hz、60～1 000 Hz 頻段的阻抗相角靈敏度曲線幾乎與橫軸重合，說(shuō)明在這些頻段J 的變化不會(huì)使阻抗相角產(chǎn)生明顯變化。上述分析與J 變化時(shí)阻抗曲線發(fā)生變化的頻段一致。

圖9 J 的阻抗靈敏度Fig.9 Impedance sensitivity of J

3.2 穩(wěn)定性分析

圖10 為MMC 并網(wǎng)系統(tǒng)示意，MMC 和電網(wǎng)參數(shù)見(jiàn)表2。輸電線路和變壓器阻抗等效為電網(wǎng)阻抗，包括等效電網(wǎng)電感Lg1、Lg2和等效電網(wǎng)電容Cg。基于VSG 控制策略的MMC 可以視為電壓源，此時(shí)MMC 并網(wǎng)系統(tǒng)的簡(jiǎn)化拓?fù)淙鐖D11 所示。

圖10 MMC 并網(wǎng)系統(tǒng)示意Fig.10 Schematic of grid-connected MMC system

圖11 MMC 并網(wǎng)系統(tǒng)簡(jiǎn)化拓?fù)銯ig.11 Simplified topology of grid-connected MMC system

在圖11 中，流過(guò)傳輸線路電流I 為

式中：Vc為VSG 策略下MMC 的等效電壓源；Vg為電網(wǎng)電壓；Zg為電網(wǎng)阻抗；ZMMC為MMC 交流端口小信號(hào)阻抗。在式（39）中，Vc和Vg是穩(wěn)態(tài)量，而Zg是無(wú)源阻抗，因此（Vc-Vg）/Zg也是穩(wěn)態(tài)量。此時(shí)系統(tǒng)穩(wěn)定性取決于1/[1+（ZMMC/Zg）]。根據(jù)奈奎斯特穩(wěn)定判據(jù)[10]，當(dāng)MMC 和電網(wǎng)都是獨(dú)立穩(wěn)定的條件下，通過(guò)阻抗比ZMMC/Zg即可判斷系統(tǒng)的穩(wěn)定性。在MMC 與電網(wǎng)的阻抗交點(diǎn)圖中，若MMC 與電網(wǎng)阻抗幅值交點(diǎn)對(duì)應(yīng)的相角裕度較大，則系統(tǒng)是穩(wěn)定的，若相角裕度很小，則系統(tǒng)將在該頻率附近振蕩，若相角裕度小于等于0°，則系統(tǒng)將發(fā)生振蕩并發(fā)散。

圖12 為MMC 阻抗與電網(wǎng)的阻抗伯德圖。二者的幅頻曲線在頻率為32、89、92 和110 Hz 有交點(diǎn)，穩(wěn)定裕度分別為11.1°、3.0°、165.8°和0.1°，其中，在110 Hz 穩(wěn)定裕度很小，系統(tǒng)在110 Hz 附近將會(huì)振蕩。

圖12 MMC 并網(wǎng)系統(tǒng)中MMC 阻抗與交流電網(wǎng)阻抗Fig.12 MMC impedance and AC grid impedance in grid-connected MMC system

3.3 實(shí)驗(yàn)驗(yàn)證

通過(guò)硬件在環(huán)實(shí)驗(yàn)驗(yàn)證上述分析。硬件在環(huán)實(shí)驗(yàn)系統(tǒng)如圖13 所示，主要由主計(jì)算機(jī)、目標(biāo)機(jī)、DSP、I/O 卡和示波器等組成。首先在主計(jì)算機(jī)上用Matalab/Simulink 的工程模式進(jìn)行編譯，將代碼下載到目標(biāo)機(jī)中進(jìn)行實(shí)時(shí)仿真（步長(zhǎng)為25 μs）。目標(biāo)機(jī)通過(guò)與其連接的多功能I/O 卡（型號(hào)Humusoft MF 624）與DSP 進(jìn)行通信。包含環(huán)流抑制控制環(huán)節(jié)的VSG 控制策略在DSP（型號(hào)TMS320F28335）中單獨(dú)運(yùn)行，目標(biāo)機(jī)中的實(shí)時(shí)模型生成的電壓和電流信號(hào)通過(guò)嵌入在目標(biāo)機(jī)器中的多功能I/O 卡適當(dāng)放大并輸出。利用示波器實(shí)時(shí)監(jiān)測(cè)和記錄目標(biāo)機(jī)提供的電壓和電流波形。PCC 點(diǎn)a 相電流實(shí)驗(yàn)波形如圖14 所示。

圖13 硬件在環(huán)實(shí)驗(yàn)平臺(tái)示意Fig.13 Schematic of hardware-in-the-loop experimental platform

圖14 PCC 點(diǎn)a 相電流實(shí)驗(yàn)波形Fig.14 Experimental waveforms of PCC current in phase a

由圖14 可以看出，加入電網(wǎng)阻抗前，系統(tǒng)運(yùn)行穩(wěn)定，當(dāng)2 s 加入電網(wǎng)阻抗后，PCC 點(diǎn)電流出現(xiàn)振蕩，說(shuō)明此時(shí)系統(tǒng)不穩(wěn)定，對(duì)PCC 點(diǎn)a 相電流進(jìn)行FFT 分析，結(jié)果如圖15 所示。從頻譜圖可以看出，PCC 點(diǎn)a 相電流在110 Hz 處諧波含量大，這與圖12 阻抗交點(diǎn)伯德圖的分析結(jié)果一致，說(shuō)明上述穩(wěn)定性分析是正確的。

圖15 PCC 點(diǎn)a 相電流諧波頻譜Fig.15 Harmonic spectrum of PCC current in phase a

4 結(jié)論

（1）基于VSG 控制策略的MMC 的阻抗特性在低頻段較為復(fù)雜，在容性和感性多次變化，高頻段主要呈感性。

（2）VSG 控制參數(shù)中的虛擬轉(zhuǎn)動(dòng)慣量和阻尼系數(shù)對(duì)50 Hz 附近的MMC 阻抗有影響，有功-下垂系數(shù)較大時(shí)對(duì)50 Hz 附近的MMC 阻抗有比較明顯的影響。

（3）當(dāng)電網(wǎng)等效阻抗在高頻段呈容性時(shí)，可能會(huì)與MMC 阻抗存在諧振點(diǎn)，導(dǎo)致MMC 并網(wǎng)系統(tǒng)在中高頻段出現(xiàn)諧振現(xiàn)象。