唐振東，楊洪耕

（四川大學(xué) 電氣信息學(xué)院，四川 成都 610065）

0 引言

風(fēng)力發(fā)電是目前成本最接近常規(guī)電力、發(fā)展前景最大的可再生能源發(fā)電。近年來，隨著風(fēng)電場開發(fā)規(guī)模的擴(kuò)大以及并網(wǎng)技術(shù)的發(fā)展，風(fēng)電并網(wǎng)所引發(fā)的并網(wǎng)點(diǎn)電壓偏差與波動(dòng)、電流諧波注入、三相不平衡等電能質(zhì)量問題愈發(fā)受到關(guān)注［1-3］。

為抑制風(fēng)電場逆變器工作時(shí)產(chǎn)生的諧波，需要在并網(wǎng)逆變器與電網(wǎng)之間加入濾波器，LCL濾波器因具有較強(qiáng)的高頻抑制能力而被廣泛使用［4-5］。但LCL濾波器的高階特性使得多逆變器并網(wǎng)時(shí)，并網(wǎng)系統(tǒng)存在多個(gè)諧振點(diǎn)，嚴(yán)重劣化并網(wǎng)電能質(zhì)量［6-7］。多逆變器與電網(wǎng)相互作用而產(chǎn)生的諧振問題目前已成為國內(nèi)外研究的熱點(diǎn)。文獻(xiàn)［8-9］建立大型新能源站的諾頓等效模型，基于所推導(dǎo)的輸出電流和并網(wǎng)點(diǎn)電壓表達(dá)式，研究了電站中由無功補(bǔ)償裝置、輸電線路阻抗等組成的電網(wǎng)阻抗對(duì)逆變器并聯(lián)系統(tǒng)諧振現(xiàn)象及電能質(zhì)量的影響；文獻(xiàn)［10］為建立更逼近實(shí)際系統(tǒng)的并網(wǎng)逆變器輸出阻抗模型，將死區(qū)效應(yīng)和開關(guān)器件非理想特性考慮到模型中以增加準(zhǔn)確性，并對(duì)多逆變器并網(wǎng)的分布式阻抗網(wǎng)絡(luò)進(jìn)行建模和諧波交互分析。

目前，分析諧波諧振的方法主要有頻率掃描法和模態(tài)分析法，前者雖可準(zhǔn)確判定諧振頻率卻難以給出更多的諧振信息［11-12］。2005年，由Xu提出的模態(tài)分析法將系統(tǒng)節(jié)點(diǎn)導(dǎo)納矩陣進(jìn)行特征值分解識(shí)別出模式諧振頻率以及節(jié)點(diǎn)參與因子［13-15］，從而給出更多的諧振抑制、避免等信息。模態(tài)分析法已逐漸成為目前分析諧波諧振現(xiàn)象的有效手段，并被廣泛應(yīng)用于微網(wǎng)諧振以及高速鐵路牽引供電系統(tǒng)諧波諧振的分析中，卻尚未被應(yīng)用于有關(guān)風(fēng)電場并網(wǎng)系統(tǒng)的諧振現(xiàn)象研究中。

基于上述分析，本文建立風(fēng)電場逆變器并網(wǎng)等效輸出阻抗模型和電網(wǎng)等效模型，采用模態(tài)分析法給出并網(wǎng)系統(tǒng)各節(jié)點(diǎn)諧振參與因子值。最后，依據(jù)實(shí)例風(fēng)電場數(shù)據(jù)，在PSCAD平臺(tái)中搭建仿真模型，結(jié)果證明模態(tài)分析法分析并網(wǎng)諧振現(xiàn)象的可行性。

1 風(fēng)電場并網(wǎng)拓?fù)浣Y(jié)構(gòu)

1.1 并網(wǎng)逆變器等效模型

鑒于LCL濾波器體積較小且對(duì)高頻諧波有顯著的抑制作用［16］，含LCL濾波器的逆變器被廣泛應(yīng)用于風(fēng)電場中。單臺(tái)三相LCL逆變器并網(wǎng)控制結(jié)構(gòu)圖如圖1所示。圖中，Lg、Linv分別為網(wǎng)側(cè)和逆變器側(cè)電感；Cf、Cdc分別為濾波電容和直流側(cè)電容；iga、iinva分別為網(wǎng)側(cè)和逆變器側(cè)a相電流；ica為濾波電容a相電流；idc為直流側(cè)電流。

圖1 并網(wǎng)逆變器控制系統(tǒng)結(jié)構(gòu)Fig.1 Structure of control system for grid-connected inverter

LCL濾波器的固有諧振尖峰會(huì)導(dǎo)致并網(wǎng)逆變器處于無阻尼不穩(wěn)定狀態(tài)，為保證系統(tǒng)能夠穩(wěn)定可靠地工作，通常采用電容電流內(nèi)環(huán)反饋與并網(wǎng)電流外環(huán)反饋相結(jié)合的雙電流閉環(huán)控制策略［17-18］。依據(jù)文獻(xiàn)［19］，利用Clarke和Park矩陣將三相靜止坐標(biāo)系下的并網(wǎng)逆變器數(shù)學(xué)模型變換到dq旋轉(zhuǎn)坐標(biāo)系下，同時(shí)依據(jù)前饋解耦策略可將d軸和q軸控制回路解耦成為相互獨(dú)立且完全對(duì)稱的兩部分。逆變器并網(wǎng)系統(tǒng)控制框圖如圖2所示。其中，iref為并網(wǎng)電流參考值；Gig（s）為并網(wǎng)電流外環(huán)控制傳遞函數(shù)，采用PI控制；Gic（s）為電容電流內(nèi)環(huán)控制傳遞函數(shù)，采用比例控制；Ginv為逆變器橋路等效增益，且有Ginv=Uin／Um，Uin為輸入直流電壓，Um為PWM波幅值。

圖2 逆變器并網(wǎng)控制框圖Fig.2 Block diagram of grid-connected inverter control

根據(jù)上述雙閉環(huán)控制框圖，由Mason定理推導(dǎo)并網(wǎng)逆變器等效輸出阻抗：

由式（1）知，并網(wǎng) LCL逆變器輸出外特性可等效表示為如圖3所示的諾頓等效電路。圖中，Go為受控電流源控制系數(shù)。

圖3 LCL逆變器并網(wǎng)諾頓模型Fig.3 Norton model of gridconnected LCL inverter

1.2 電網(wǎng)模型

可控串聯(lián)補(bǔ)償器可以平滑快速地改變輸電線路的阻抗，具有提高輸送功率、改善潮流分布和改善系統(tǒng)動(dòng)態(tài)性能等優(yōu)點(diǎn)，因此110 kV及以上電壓等級(jí)的長距離輸電網(wǎng)絡(luò)及其電網(wǎng)可等效為理想電壓源［20-21］。建立電網(wǎng)等效模型如圖4所示。圖中，Rg為輸電線路等效電阻；Lg為輸電線路和變壓器等效電感之和；LTn為升壓變壓器Tn等效電感；CQn、LQn和RQn分別為第n個(gè)風(fēng)電機(jī)組集群并網(wǎng)時(shí)所聯(lián)結(jié)的無功補(bǔ)償裝置總等效電容、等效電感和等效電阻；uinvn為逆變器集群并網(wǎng)點(diǎn)電壓；ug為電網(wǎng)電壓。

圖4 電網(wǎng)等效模型Fig.4 Equivalent grid model

2 模態(tài)分析法及其指標(biāo)量

2.1 模態(tài)分析法

下面分析模態(tài)分析法的核心思想。若系統(tǒng)產(chǎn)生頻率為f的并聯(lián)諧振現(xiàn)象，其節(jié)點(diǎn)電壓、電流方程為［13］：

其中，Yf為在頻率f處的系統(tǒng)節(jié)點(diǎn)導(dǎo)納矩陣；Vf和If分別為節(jié)點(diǎn)電壓和節(jié)點(diǎn)注入電流向量。

節(jié)點(diǎn)導(dǎo)納矩陣Yf可分解如下：

其中，Λ 為對(duì)角特征值矩陣，Λ=diag（λ1，λ2，…，λn）；L= ［L1L2… Ln］、T=［T1T2… Tn］分別為左、右特征向量矩陣，且有L=T-1。

將式（3）代入式（2）得：

其中，Λ-1命名為“模態(tài)阻抗”。

定義Uf=TVf為模態(tài)電壓向量，Jf=TIf為模態(tài)電流向量，則有 Uf=Λ-1Jf，即：

由式（5）知，當(dāng)系統(tǒng)發(fā)生并聯(lián)諧振時(shí)，某些節(jié)點(diǎn)的注入電流將產(chǎn)生很高的節(jié)點(diǎn)電壓，即若λf1=0或非常小，則很小的模態(tài)1注入電流Jf1將導(dǎo)致很大的模態(tài)電壓Uf1，而其他模態(tài)電壓不受影響。即在模態(tài)域中，根據(jù)特征值倒數(shù)易于識(shí)別出諧振的位置。

2.2 模態(tài)分析法指標(biāo)量

諧波諧振實(shí)際上只在特定的模式下發(fā)生，它與某個(gè)節(jié)點(diǎn)的注入電流無關(guān)也并非由其引起?？蓪⒆钚〉奶卣髦捣Q為諧振“關(guān)鍵模式”，對(duì)應(yīng)的左、右特征向量則稱為“關(guān)鍵特征向量”。

模態(tài)電流Jf1可表示為電流在第一個(gè)特征向量上的線性映射，即：

由式（6）可知，若T1n有最大值，則節(jié)點(diǎn)電流Ifn將對(duì)模態(tài)1注入電流具有最大的貢獻(xiàn)度，即節(jié)點(diǎn)n最容易受到外界激勵(lì)而產(chǎn)生模態(tài)1諧振。

同理可得，模態(tài)電壓與節(jié)點(diǎn)電壓間存在關(guān)系Vf=LUf，則有：

由式（7）可知，若系統(tǒng)發(fā)生模態(tài)1諧振，則Uf1遠(yuǎn)大于其他模態(tài)電壓，此時(shí)，式（7）可近似僅由［L11L21… Ln1］TUf1來表示。如果Ln1的值最大，則節(jié)點(diǎn)n將會(huì)出現(xiàn)最大電壓值，即節(jié)點(diǎn)n處最容易觀測到模態(tài)1諧振。

綜上可得結(jié)論：右特征向量矩陣T反映關(guān)鍵模式下節(jié)點(diǎn)的諧振可激勵(lì)性，左特征向量矩陣L反映關(guān)鍵模式下節(jié)點(diǎn)的諧振可觀測性，在具有最大激勵(lì)性的節(jié)點(diǎn)注入相應(yīng)信號(hào)可有效抑制諧波諧振現(xiàn)象。

實(shí)際上，可將可激勵(lì)性和可觀測性結(jié)合成為一個(gè)新指標(biāo)量。當(dāng)系統(tǒng)發(fā)生模態(tài)1諧振時(shí)，模態(tài)阻抗 1／λf1的值遠(yuǎn)大于其他模態(tài)阻抗值，則式（4）可簡化為：

式（8）中矩陣對(duì)角線元素顯示了節(jié)點(diǎn)在關(guān)鍵模式下的可激勵(lì)性和可觀測性的結(jié)合，不妨稱之為關(guān)鍵模式下節(jié)點(diǎn)的“參與因子”，參與因子反映了節(jié)點(diǎn)對(duì)所發(fā)生諧振的參與度以及貢獻(xiàn)度，對(duì)某種模式具有最大可觀測性的節(jié)點(diǎn)同時(shí)也具有最大可激勵(lì)性。具有最大參與因子的節(jié)點(diǎn)則可被認(rèn)為是系統(tǒng)諧波諧振的中心。

3 風(fēng)電場并網(wǎng)諧振現(xiàn)象模態(tài)分析

3.1 模態(tài)分析法應(yīng)用算例

為驗(yàn)證模態(tài)分析法分析諧振現(xiàn)象的準(zhǔn)確性，以圖5所示三節(jié)點(diǎn)測試系統(tǒng)為例簡要說明。各元件參數(shù)（標(biāo)幺值）標(biāo)示于圖中，三節(jié)點(diǎn)均可被激勵(lì)或觀察諧波諧振。圖6為三節(jié)點(diǎn)測試系統(tǒng)的模態(tài)分析結(jié)果。由圖6可知，測試系統(tǒng)共有3個(gè)準(zhǔn)諧振點(diǎn)，分別為 145 Hz、335 Hz和 775 Hz處。

圖5 測試系統(tǒng)Fig.5 Test system

圖6 測試系統(tǒng)模態(tài)分析結(jié)果Fig.6 Results of modal analysis for test system

圖7為采用傳統(tǒng)頻率掃描法對(duì)各節(jié)點(diǎn)的分析結(jié)果，分別為從節(jié)點(diǎn)1、2、3看進(jìn)去的等效輸出阻抗的頻率掃描特性。圖7中第一圖顯示系統(tǒng)準(zhǔn)諧振點(diǎn)為775Hz；圖7中第二、三圖顯示系統(tǒng)在145 Hz和335 Hz均可能發(fā)生諧波諧振。對(duì)比圖6與圖7可知，兩者所得諧振點(diǎn)一致，證明了模態(tài)分析法在分析諧振現(xiàn)象時(shí)的可行性和準(zhǔn)確性。同時(shí)，模態(tài)分析法中，對(duì)每個(gè)特定的諧振頻率，只有一個(gè)諧振模式與之對(duì)應(yīng)，由此可進(jìn)一步計(jì)算各節(jié)點(diǎn)的諧振參與因子，從而得到更多的諧振相關(guān)信息。

圖7 測試系統(tǒng)頻率掃描結(jié)果Fig.7 Results of frequency scan for test system

3.2 基于模態(tài)分析法的風(fēng)電場并網(wǎng)諧振分析

建立多逆變器并網(wǎng)模型如圖8所示。圖中，Zg（s）為電網(wǎng)等效阻抗；Zfk（s）、Zok（s）和 Gkirefk（s）分別為線路阻抗、逆變器并網(wǎng)等效輸出阻抗和諾頓等效電路受控電流源，k=1，2，…，n-1。

圖8 逆變器并網(wǎng)模態(tài)分析模型Fig.8 Modal analysis model of grid-connected inverter

采用模態(tài)分析法研究并網(wǎng)諧振現(xiàn)象的具體步驟如下。

a.由圖8所示模態(tài)分析模型，可求得系統(tǒng)節(jié)點(diǎn)導(dǎo)納矩陣Y（s）。因節(jié)點(diǎn)導(dǎo)納矩陣本身具有對(duì)稱性，為簡化表達(dá)式，僅列出其上三角元素，即：

b.對(duì)式（9）所示的節(jié)點(diǎn)導(dǎo)納矩陣Y（s）進(jìn)行特征值分解，即 Y（s）=L（s）Λ（s）T（s），從而得到左、右特征向量矩陣 L（s）、T（s）及對(duì)角特征值矩陣 Λ（s）。

c.分別用 L（s）、T（s）表示模態(tài)電壓和模態(tài)電流向量，即有 U（s）=T（s）V（s）、J（s）=T（s）I（s）。 同時(shí)，計(jì)算系統(tǒng)頻率為f時(shí)各模式下的模態(tài)阻抗值。

d.重復(fù)步驟b、c直至遍歷所有可能引起系統(tǒng)諧振的諧波頻率，從而得到各諧波頻率下的系統(tǒng)模態(tài)阻抗曲線圖，進(jìn)而分析風(fēng)電場并網(wǎng)諧振現(xiàn)象。

e.模態(tài)阻抗值最大的點(diǎn)即為關(guān)鍵諧振峰值點(diǎn)，根據(jù)式（8）將其所對(duì)應(yīng)的左、右特征向量依次取對(duì)應(yīng)對(duì)角線元素，求積可得各節(jié)點(diǎn)關(guān)鍵模式下的參與因子，并根據(jù)參與因子值的大小判斷系統(tǒng)諧振中心。

4 實(shí)例分析

4.1 實(shí)例介紹

為檢驗(yàn)?zāi)B(tài)分析法分析風(fēng)電場并網(wǎng)諧振現(xiàn)象的可行性，以某風(fēng)電場為實(shí)例進(jìn)行仿真驗(yàn)證，其主接線圖如圖9所示。風(fēng)電場裝機(jī)容量為99 MW，風(fēng)電機(jī)組型號(hào)為華銳SL1500，機(jī)組類型為水平軸異步風(fēng)力發(fā)電機(jī)，單機(jī)容量為1.5 MW。風(fēng)電場分為兩期建設(shè)，每期均安裝有33臺(tái)風(fēng)電機(jī)組。為便于后文通過節(jié)點(diǎn)導(dǎo)納矩陣實(shí)現(xiàn)并網(wǎng)諧振的模態(tài)分析，對(duì)實(shí)例系統(tǒng)進(jìn)行節(jié)點(diǎn)序號(hào)的標(biāo)注，將可能存在同類型諧波電流激勵(lì)源的節(jié)點(diǎn)進(jìn)行并聯(lián)處理，共選取9個(gè)節(jié)點(diǎn)，其中，節(jié)點(diǎn)1—3與節(jié)點(diǎn)5—7為實(shí)際風(fēng)電機(jī)組集電線聯(lián)接處，每條集電線節(jié)點(diǎn)均聯(lián)接有11臺(tái)風(fēng)機(jī)。風(fēng)電機(jī)組出口電壓為690 V，經(jīng)箱式變壓器升壓至35 kV后，再通過主變變換至110 kV由單回架空線并網(wǎng)。并網(wǎng)系統(tǒng)其他參數(shù)如表1所示。

4.2 模態(tài)分析結(jié)果

根據(jù)表1所示實(shí)例系統(tǒng)并網(wǎng)參數(shù)，基于PSCAD平臺(tái)并采用平臺(tái)中的變壓器、輸電線路等模塊搭建仿真模型，風(fēng)電場并網(wǎng)諧振的模態(tài)分析結(jié)果如圖10所示（h為諧振次數(shù)，模態(tài)阻抗為標(biāo)幺值；后同）。分析運(yùn)算中各并網(wǎng)參數(shù)均采用標(biāo)幺值，并與MATLAB編程相結(jié)合，節(jié)點(diǎn)選取則以圖9所標(biāo)注的9個(gè)系統(tǒng)節(jié)點(diǎn)為準(zhǔn)。為便于后文分析，假設(shè)此并網(wǎng)風(fēng)電場為滿出力運(yùn)行狀態(tài)。

圖9 實(shí)例系統(tǒng)示意圖Fig.9 Schematic diagram of a power system

表1 系統(tǒng)參數(shù)Table 1 System parameters

圖10 模態(tài)分析結(jié)果Fig.10 Results of modal analysis

由圖10可知，并網(wǎng)系統(tǒng)可能在3種模式下發(fā)生諧波諧振現(xiàn)象，準(zhǔn)諧振點(diǎn)及其對(duì)應(yīng)的諧振模態(tài)阻抗值分別為模式 3（27.2，378.9）、模式 4（29.1，95.6）和模式5（15.3，425.7）。表2詳細(xì)給出了在3種關(guān)鍵諧振模式下各節(jié)點(diǎn)的諧振參與因子值。

表2 關(guān)鍵模式參與因子表Table 2 Participation factors for key modes

由表2可知，因搭建仿真模型時(shí)，35 kV母線I、Ⅱ段呈完全對(duì)稱狀態(tài)，故3種諧振頻率下節(jié)點(diǎn)1、2、3、4 與節(jié)點(diǎn) 5、6、7、8 的諧振參與因子值對(duì)應(yīng)相等。且此8個(gè)節(jié)點(diǎn)的參與因子值相差不多，即發(fā)生諧振時(shí)，8個(gè)節(jié)點(diǎn)將會(huì)受到同等程度的影響。節(jié)點(diǎn)9在諧振模式5下具有最大的參與因子值，則可將其作為此模式下的諧振中心。而在諧振模式3和模式4下，節(jié)點(diǎn)9均具有最小的參與因子值，可認(rèn)為此時(shí)節(jié)點(diǎn)9相較于其他節(jié)點(diǎn)最不易觀測到諧振現(xiàn)象，同時(shí)在節(jié)點(diǎn)9也最不易被激勵(lì)起諧波諧振。

4.3 動(dòng)態(tài)條件下的風(fēng)電場并網(wǎng)諧振分析

集電線作為風(fēng)電機(jī)組與電力系統(tǒng)的聯(lián)絡(luò)線，其線路距離的長短對(duì)并網(wǎng)諧振現(xiàn)象存在著不可忽略的影響。繼續(xù)以上述3種關(guān)鍵諧振模式為例，分析當(dāng)集電線路長度由5 km漸增至20 km時(shí)風(fēng)電場并網(wǎng)諧振變化情況，如圖11所示（l為集電線長度）。

圖11 集電線長度改變時(shí)的模態(tài)分析結(jié)果Fig.11 Results of modal analysis vs.tieline length

由圖11可知，模式諧振與節(jié)點(diǎn)諧振沒有直接的對(duì)應(yīng)關(guān)系，某一模式諧振可能包含有較多的參與節(jié)點(diǎn)。圖11（a）顯示，當(dāng)集電線長度在 8.2～14.7 km時(shí)，諧振曲線1出現(xiàn)，所對(duì)應(yīng)的諧振次數(shù)h由29.4變?yōu)?6.7，諧振幅值也隨之驟減；當(dāng)集電線長度在11.5～20 km時(shí)，諧振曲線2出現(xiàn)，所對(duì)應(yīng)的諧振次數(shù)h由15.8減至11.3。圖11（b）顯示，隨著集電線長度由5 km增至11 km，曲線3所對(duì)應(yīng)諧振次數(shù)h恒為29.1；當(dāng)集電線長度增至15 km時(shí)，諧振曲線4出現(xiàn)，所對(duì)應(yīng)諧振次數(shù)h為11，諧振幅值呈現(xiàn)波動(dòng)性變化。圖11（c）顯示，諧振曲線5和曲線6分別在集電線長度為8.9 km和11.2 km處諧振幅值減至零，即在此處諧波諧振現(xiàn)象將會(huì)消失。

一般而言，風(fēng)機(jī)的運(yùn)行與否和所處地帶風(fēng)速的大小有很大關(guān)系，當(dāng)風(fēng)速在3～25 m／s時(shí)，風(fēng)機(jī)進(jìn)入正常運(yùn)行狀態(tài)。以集電線1為例，考慮并非所有風(fēng)機(jī)都正常并網(wǎng)運(yùn)行，即集電線所接風(fēng)機(jī)臺(tái)數(shù)小于11臺(tái)時(shí)的情況，其并網(wǎng)諧振頻率以及各節(jié)點(diǎn)諧振參與因子如表3所示。

由表3可知，隨著集電線1所聯(lián)結(jié)風(fēng)機(jī)臺(tái)數(shù)的增加，諧振模式5和模式3所對(duì)應(yīng)諧振次數(shù)逐漸減小。諧振模式5下，節(jié)點(diǎn)9的參與因子值始終保持最大，即恒為諧振中心，在此處最易被激勵(lì)起諧振也最易觀測到諧振現(xiàn)象。而在關(guān)鍵模式3下，節(jié)點(diǎn)1—8的參與因子值相差無幾且都遠(yuǎn)大于節(jié)點(diǎn)9，即發(fā)生諧波諧振時(shí)，節(jié)點(diǎn)1—8均受影響嚴(yán)重。參與因子值同時(shí)也反映出諧波諧振的傳播距離。

通過比較可發(fā)現(xiàn)，在集電線路長度和并網(wǎng)風(fēng)機(jī)臺(tái)數(shù)改變的情況下，風(fēng)電場并網(wǎng)諧振次數(shù)和諧振幅值都在動(dòng)態(tài)變化著。而節(jié)點(diǎn)參與因子、諧振中心等信息的提取可為風(fēng)電場治理、諧振抑制提供一定的理論指導(dǎo)。

表3 并網(wǎng)風(fēng)機(jī)臺(tái)數(shù)改變時(shí)的參與因子表Table 3 Participation factors for different grid-connected generator quantities

5 結(jié)論

針對(duì)大規(guī)模風(fēng)電場并網(wǎng)時(shí)易發(fā)生的諧波諧振問題，本文建立并網(wǎng)逆變器等效輸出阻抗模型和電網(wǎng)等效模型，并基于此采用模態(tài)分析法進(jìn)行諧振研究，得出以下結(jié)論：

a.相比于傳統(tǒng)的頻率掃描法，模態(tài)分析法不僅能準(zhǔn)確判定并網(wǎng)系統(tǒng)諧振次數(shù)，還可給出各節(jié)點(diǎn)諧振參與因子、諧振中心等諧振信息；

b.集電線長度以及并網(wǎng)風(fēng)機(jī)臺(tái)數(shù)的改變都將很大程度地改變風(fēng)電場并網(wǎng)諧振的諧振次數(shù)、諧振幅值以及諧振中心的位置；

c.實(shí)例仿真驗(yàn)證了本文所用模態(tài)分析法的可行性，所得結(jié)論可為工程實(shí)際中抑制諧波諧振、評(píng)估風(fēng)電場電能質(zhì)量等方面提供理論指導(dǎo)。

參考文獻(xiàn)：

［1］MOHAMED Y A R．Mitigation of converter-grid resonance，gridinduced distortion，and parametric instabilities in converter-based distributed generation［J］．IEEE Transactions on Power Electronics，2011，26（3）：983-996．

［2］于汀，蒲天驕，劉廣一，等．含大規(guī)模風(fēng)電的電網(wǎng)AVC研究與應(yīng)用［J］．電力自動(dòng)化設(shè)備，2015，35（10）：81-86．YU Ting，PU Tianjiao，LIU Guangyi，et al．AVC research and application for grid with large-scale wind power［J］．Electric Power Automation Equipment，2015，35（10）：81-86．

［3］潘偉，李勇，曹一家，等．用于大規(guī)模集中式風(fēng)電并網(wǎng)的VSCHVDC 頻率控制方法［J］．電力自動(dòng)化設(shè)備，2015，35（5）：94-99．PAN Wei，LI Yong，CAO Yijia，et al．Frequency control of gridconnection system based on VCS-HVDC for large-scale centralized wind farm［J］．Electric Power Automation Equipment，2015，35（5）：94-99．

［4］許斌，張志強(qiáng)，李程昊，等．應(yīng)用于雙饋風(fēng)力發(fā)電機(jī)的LCL型濾波器［J］．電力自動(dòng)化設(shè)備，2015，35（5）：44-50．XU Bin，ZHANG Zhiqiang，LI Chenghao，et al．LCL filters applied in doubly fed induction generator［J］．Electric Power Automation Equipment，2015，35（5）：44-50．

［5］JALILI K，BERNET S．Design of LCL filters of active-front-end two level voltage-source converters［J］．IEEE Transactions on Industrial Electronics，2009，56（5）：1674-1689．

［6］WANG F，DUARTE J，HENDRIX M，et al．Modeling and analysis of grid harmonic distortion impact of aggregated DG inverters［J］．IEEE Transactions on Power Electronics，2011，26（3）：786-797．

［7］HE J．Investigation and active damping of multiple resonances in a parallel-inverter-based microgrid［J］．IEEE Transactions on Power Electronics，2013，28（1）：234-246．

［8］謝寧，羅安，馬伏軍，等．大型光伏電站與電網(wǎng)諧波交互影響［J］．中國電機(jī)工程學(xué)報(bào)，2013，33（34）：9-16．XIE Ning，LUO An，MA Fujun，et al．Harmonic interaction between large-scale photovoltaic power stations and grid［J］．Proceedings of the CSEE，2013，33（34）：9-16．

［9］周林，張密．大型光伏電站諧振現(xiàn)象分析［J］．電力自動(dòng)化設(shè)備，2014，34（6）：8-14．ZHOU Lin，ZHANG Mi．Analysis of resonance phenomenon in large-scale photovoltaic power plant［J］．Electric Power Automation Equipment，2014，34（6）：8-14．

［10］許德志，汪飛，毛華龍，等．多并網(wǎng)逆變器與電網(wǎng)的諧波交互建模與分析［J］．中國電機(jī)工程學(xué)報(bào)，2013，33（12）：64-71．XU Dezhi，WANG Fei，MAO Hualong，etal．Modelingand analysis of harmonic interaction between multiple grid-connected inverters and the utility grid［J］．Proceedings of the CSEE，2013，33（12）：64-71．

［11］張民，何正友，高仕斌．不同負(fù)荷模型下高速鐵路牽引供電系統(tǒng)諧波諧振敏感度分析［J］．電力系統(tǒng)保護(hù)與控制，2015，43（20）：21-28．ZHANG Min，HE Zhengyou，GAO Shibin．Harmonic resonance modal sensitivity analysis for high-speed railway traction supply system considering different harmonic load models［J］．Power System Protection and Control，2015，43（20）：21-28．

［12］胡海濤，何正友，錢澄浩，等．基于模態(tài)分析的全并聯(lián)AT網(wǎng)動(dòng)態(tài)諧波諧振研究［J］．電網(wǎng)技術(shù)，2012，36（11）：163-169．HU Haitao，HE Zhengyou，QIAN Chenghao，et al．Modal analysis based research on dynamic harmonic resonanceofallparallel autotransformer traction supply network［J］．Power System Technology，2012，36（11）：163-169．

［13］徐文遠(yuǎn)，張大海．基于模態(tài)分析的諧波諧振評(píng)估方法［J］．中國電機(jī)工程學(xué)報(bào)，2005，25（22）：89-93．XU Wenyuan，ZHANG Dahai．A modal analysis method for harmonic resonance assessment［J］．Proceedings of the CSEE，2005，25（22）：89-93．

［14］HUANG Z，CUI Y，XU W．Application of modal sensitivity for power system harmonic resonance analysis［J］．IEEE Transactions on Power Systems，2007，22（2）：222-231．

［15］CUI Y，XU W．Harmonic resonance mode analysis using real symmetrical nodal matrices［J］．IEEE Transactions on Power Delivery，2007，22（3）：1989-1990．

［16］鄭昕昕，肖嵐，王長寶，等．三相變流器LCL濾波器參數(shù)優(yōu)化新方法［J］．中國電機(jī)工程學(xué)報(bào)，2013，33（12）：55-63．ZHENG Xinxin，XIAO Lan，WANG Changbao，et al．A new parameter optimization method of LCL filters in three-phase converters［J］．Proceedings of the CSEE，2013，33（12）：55-63．

［17］許德志，汪飛，阮毅．LCL、LLCL和LLCCL濾波器無源阻尼分析［J］．中國電機(jī)工程學(xué)報(bào)，2015，35（18）：4725-4735．XU Dezhi，WANG Fei，RUAN Yi．Passive damping of LCL、LLCL and LLCCL filters［J］．Proceedings of the CSEE，2015，35（18）：4725-4735．

［18］胡偉，孫建軍，馬謙，等．多個(gè)逆變器間的交互影響分析［J］．電網(wǎng)技術(shù)，2014，38（9）：2511-2518．HU Wei，SUN Jianjun，MA Qian，et al．Analysis on interactive influences among multi grid-connected inverters［J］．Power System Technology，2014，38（9）：2511-2518．

［19］許德志，汪飛，阮毅，等．多逆變器并網(wǎng)系統(tǒng)輸出阻抗建模與諧波交互［J］．電機(jī)與控制學(xué)報(bào)，2014，18（2）：1-7．XU Dezhi，WANG Fei，RUAN Yi，etal．Outputimpedance modeling and harmonic interactions of multiple inverters gridconnected system ［J］．Electric Machines and Control，2014，18（2）：1-7．

［20］鄭旭，裴志宏，孫海順，等．可控串補(bǔ)動(dòng)模裝置研制及阻抗特性研究［J］．電力自動(dòng)化設(shè)備，2005，25（5）：50-53．ZHENG Xu，PEI Zhihong，SUN Haishun，et al．Design and implement of dynamic simulation system of TCSC and study on resistance characteristics［J］．Electric Power Automation Equipment，2005，25（5）：50-53．

［21］楊勇，趙春江．分布式發(fā)電系統(tǒng)中并網(wǎng)逆變器比例諧振控制［J］．電力自動(dòng)化設(shè)備，2011，31（11）：51-55．YANG Yong，ZHAO Chunjiang．Proportionalresonancecontroller of grid-connected inverter for distributed generation system［J］．Electric Power Automation Equipment，2011，31 （11）：51-55．