劉斯偉，李庚銀，周 明

?

雙饋風電機組對接入?yún)^(qū)域系統(tǒng)暫態(tài)功角穩(wěn)定性的影響分析

劉斯偉，李庚銀，周 明

(新能源電力系統(tǒng)國家重點實驗室(華北電力大學)，北京 102206)

系統(tǒng)慣性中心的動態(tài)行為可以有效地反映系統(tǒng)的整體動態(tài)趨勢。基于系統(tǒng)慣性中心的動態(tài)數(shù)學模型，推導了雙饋風電機組(Doubly Fed Induction Generator, DFIG)接入后的系統(tǒng)慣性中心運動方程。并通過分析DFIG接入對系統(tǒng)慣性中心運動方程中關鍵因素的影響，揭示了DFIG以不同方式接入系統(tǒng)時對系統(tǒng)慣性中心功角加速度暫態(tài)響應的影響機理。當DFIG直接接入系統(tǒng)時，系統(tǒng)慣性中心的暫態(tài)功角加速度的變化趨勢主要取決于DFIG的接入容量和暫態(tài)功率響應特性。當DFIG等容替代同步發(fā)電機組接入系統(tǒng)時，DFIG的暫態(tài)功率響應和被替代同步發(fā)電機組的功率和慣量在原系統(tǒng)中的比重大小將共同決定系統(tǒng)慣性中心的暫態(tài)功角加速度的變化情況。仿真結果驗證了理論分析的有效性和正確性。

雙饋風電機組(DFIG)；系統(tǒng)慣性中心；暫態(tài)功角穩(wěn)定性；功角加速度；影響機理

0 引言

當今由于能源危機、環(huán)保壓力以及風電制造技術的迅速發(fā)展，風力發(fā)電在世界范圍內(nèi)得到了高度重視，在我國，風力發(fā)電主以大規(guī)模集中的方式并入電網(wǎng)，遠離負荷中心，傳輸距離較遠[1-3]。雙饋風電機組(Doubly Fed Induction Generator，DFIG)因其變速運行范圍寬、變流器容量小等優(yōu)點，成為我國風電場中應用最廣泛的風力發(fā)電機型[4-5]。通過引入背靠背的雙PWM(Pulse Width Modulation)變流器組實現(xiàn)交流勵磁控制，DFIG的轉(zhuǎn)速與系統(tǒng)頻率實現(xiàn)了機電解耦，并且DFIG的動態(tài)響應特性將主要取決于變流器的控制策略。DFIG與傳統(tǒng)同步發(fā)電機組在工作原理上的區(qū)別使其具有不同的動態(tài)行為特征，隨著其裝機容量的不斷增長，如何準確分析并網(wǎng)DFIG對系統(tǒng)暫態(tài)穩(wěn)定性的影響成為風電發(fā)展的一個難題[6-11]。

很多文獻關于DFIG對接入系統(tǒng)暫態(tài)穩(wěn)定性的影響展開了研究[12-18]，文獻[15]基于CIGRE 39節(jié)點系統(tǒng)模型對DFIG替代同步發(fā)電機組的情況進行了仿真，通過比較系統(tǒng)的極限切除時間(Critical Cut Time，CCT)，指出DFIG替代同步發(fā)電機組可以提高系統(tǒng)的魯棒性。文獻[16]考慮DFIG的并網(wǎng)接入點位置、聯(lián)絡線阻抗比和接入容量比作為對比條件，通過比較系統(tǒng)內(nèi)任意同步發(fā)電機組間的最大功角差，指出DFIG對系統(tǒng)內(nèi)同步發(fā)電機組暫態(tài)穩(wěn)定性的影響既可能是正面也可能是負面的。文獻[17]則研究了DFIG與系統(tǒng)中同步電機的交互作用，認為DFIG可能對電力系統(tǒng)穩(wěn)定性不利。文獻[18]通過建立含DFIG的電力系統(tǒng)暫態(tài)能量函數(shù)(Transient Energy Function, TEF)，從系統(tǒng)整體的角度分析，經(jīng)不同仿真條件的計算比較，認為同步發(fā)電機組接入相較于DFIG接入時，系統(tǒng)擁有更高的暫態(tài)穩(wěn)定性。在關于DFIG對系統(tǒng)暫態(tài)穩(wěn)定性影響的研究中，以往文獻基于不同的研究方法和系統(tǒng)模型開展了大量仿真算例和實例分析，然而得出的結論并不完全一致，雖然大部分研究并對仿真結果進行了詳細分析，然而仍未能有效地指出DFIG對接入系統(tǒng)暫態(tài)穩(wěn)定性的本質(zhì)影響機理。

系統(tǒng)慣性中心(Center of Inertia, COI)的思想廣泛應用于電力系統(tǒng)的暫態(tài)穩(wěn)定性分析[19]，其功角和轉(zhuǎn)速的動態(tài)方程能夠有效反映區(qū)域系統(tǒng)的整體動態(tài)穩(wěn)定趨勢。本文基于系統(tǒng)慣性中心的動態(tài)數(shù)學模型，對DFIG接入后的系統(tǒng)慣性中心運動方程進行了推導，并通過分析DFIG接入對運動方程中關鍵因素的影響方式，深入研究了系統(tǒng)慣性中心功角角速度的暫態(tài)變化的影響機理。最后應用時域仿真對理論分析進行有效驗證。

1 ?系統(tǒng)慣性中心動態(tài)模型

系統(tǒng)的慣性中心可以被描述為系統(tǒng)內(nèi)所有同步發(fā)電機的實時功角和轉(zhuǎn)速以轉(zhuǎn)動慣量為權重的加權平均值，定義為[20]

其中，

(2)

系統(tǒng)慣性中心的運動方程可以表示為

其中，

(4)

2 ?DFIG對系統(tǒng)暫態(tài)功角穩(wěn)定性的影響

2.1 含DFIG的系統(tǒng)慣性中心動態(tài)模型

在傳統(tǒng)電力系統(tǒng)的暫態(tài)過程中，各同步發(fā)電機組的功角和轉(zhuǎn)速通過彼此間的電氣聯(lián)系相互牽制、相互影響，共同維持著系統(tǒng)的暫態(tài)功角穩(wěn)定。而DFIG的勵磁控制過程中實現(xiàn)了機電解耦，不具備類似于同步發(fā)電機的功角特性和慣量特性，這使得DFIG既不存在功角失穩(wěn)機制，也無法直接參與同步發(fā)電機組之間的互同步穩(wěn)定機制。

因此當DFIG接入系統(tǒng)后，系統(tǒng)慣性中心的運動方程仍由系統(tǒng)內(nèi)所有同步發(fā)電機的運動方程構成，表示為

式中，上標“＇”表示DFIG接入后系統(tǒng)各電氣量的狀態(tài)。

式(5)中，DFIG接入后系統(tǒng)慣性中心的慣性時間常數(shù)、電磁功率和機械功率將發(fā)生以下改變。

式中，下標“dis”表示被替代同步發(fā)電機組的相關變量。

式中，下標“|0|”表示電氣量的穩(wěn)態(tài)值。

2.2 DFIG對系統(tǒng)暫態(tài)功角穩(wěn)定性的影響

本節(jié)進一步分析當DFIG直接接入系統(tǒng)或替代同步發(fā)電機組接入系統(tǒng)情況下，DFIG對系統(tǒng)暫態(tài)功角穩(wěn)定性的影響機理。式(3)和式(5)分別描述了DFIG接入前后系統(tǒng)慣性中心的功角加速度，令表示兩者的差值，則有：

式(11)描述了DFIG接入后系統(tǒng)慣性中心功角加速度動態(tài)響應的變化情況，下面結合DFIG的不同接入方式展開分析。

(1)?DFIG直接接入系統(tǒng)，系統(tǒng)的總負荷不變，同步發(fā)電機組數(shù)量不變，部分機組減發(fā)。此時系統(tǒng)慣性中心的慣性時間常數(shù)在DFIG接入系統(tǒng)前后保持不變，即

首先將式(7)、式(8)和式(10)代入式(9)，得到系統(tǒng)慣性中心功角加速度的變化為

(11)

考慮到系統(tǒng)的總負荷不變，同時假設DFIG的接入對系統(tǒng)暫態(tài)電壓穩(wěn)定性影響極小可忽略不計，即系統(tǒng)負荷的動態(tài)響應影響保持基本不變，則DFIG接入前后系統(tǒng)的總功率需求將滿足：

將式(12)進一步代入式(11)，得到系統(tǒng)慣性中心功角加速的變化為

(13)

分析式(13)可知，當DFIG直接接入系統(tǒng)時，系統(tǒng)慣性中心的暫態(tài)功角加速度的增減變化將主要取決于DFIG的接入容量和其暫態(tài)功率響應特性。

(2)?DFIG以等容量替代同步發(fā)電機組的方式接入系統(tǒng)，系統(tǒng)的總負荷不變，同步發(fā)電機組數(shù)量減少。設此時被替代的同步發(fā)電機組的慣性時間常數(shù)和機械功率分別為

將式(14)代入式(6)，可得系統(tǒng)慣性中心的慣性時間常數(shù)為

(15)

由于系統(tǒng)的總負荷不變，即系統(tǒng)內(nèi)所有發(fā)電機機械功率不變，根據(jù)式(14)可得DFIG和系統(tǒng)慣性中心的機械功率為

在此基礎上，假設通過采取無功補償措施，使DFIG的接入對系統(tǒng)暫態(tài)電壓穩(wěn)定性影響極小可忽略不計，從而保證系統(tǒng)負荷的動態(tài)響應維持不變，則對系統(tǒng)總功率需求的描述仍然滿足式(12)。

將式(7)、式(8)、式(12)、式(15)和(16)代入式(9)，可得系統(tǒng)慣性中心功角加速度的變化為

為進一步梳理式(17)的結構，定義：

(18)

式中，下標“p.u.”表示相關變量以自身故障前穩(wěn)態(tài)值為基準的標幺值。

將式(18)代入式(17)，則DFIG替代同步發(fā)電機組接入系統(tǒng)后，系統(tǒng)慣性中心功角加速度的變化可以進一步表達為

分析式(19)，令：

(20)

和以標幺值形式分別表示了系統(tǒng)總功率需求和DFIG輸出功率的暫態(tài)響應與自身穩(wěn)態(tài)值的差值。由此可以得到，當DFIG替代同步發(fā)電機組接入系統(tǒng)且不影響系統(tǒng)的暫態(tài)電壓穩(wěn)定性時，系統(tǒng)慣性中心的暫態(tài)功角加速度的增減變化將共同取決于DFIG的暫態(tài)功率響應，和被替代同步發(fā)電機的功率和慣量在原系統(tǒng)中的比重大小。

3 ?算例分析

3.1 算例系統(tǒng)

為了驗證上述關于DFIG對接入系統(tǒng)暫態(tài)功角穩(wěn)定性的影響分析，基于PSCAD/EMTDC仿真平臺搭建了IEEE-3機9節(jié)點仿真系統(tǒng)用于算例分析，如圖1所示。

圖1 IEEE-3機9節(jié)點系統(tǒng)示意圖

系統(tǒng)中的3臺發(fā)電機初始均由同步發(fā)電機組表示，配有調(diào)速器模型和勵磁系統(tǒng)模型；負荷均采用恒阻抗模型；同步發(fā)電機組和負荷的詳細參數(shù)如表1和表2所示。接入系統(tǒng)的DFIG風電場模型采用加權等值的方法，用1臺詳細的DFIG模型表示。系統(tǒng)4?s時于BUS6節(jié)點發(fā)生三相對稱短路故障，持續(xù)0.15?s后故障切除。基于本文的研究目的，暫態(tài)過程中假設風速保持額定風速12?m/s不變。

表1 IEEE-3機9節(jié)點系統(tǒng)發(fā)電機數(shù)據(jù)

表2 IEEE-3機9節(jié)點系統(tǒng)負荷數(shù)據(jù)

3.2 DFIG直接接入系統(tǒng)

DFIG風電場分別以30?MW和60?MW的額定容量接入至BUS9節(jié)點，系統(tǒng)中的同步發(fā)電機組不被替代，并保持系統(tǒng)中負荷不變。圖2展示了兩組不同容量DFIG風電場直接接入系統(tǒng)后，系統(tǒng)故障情況下DFIG電磁功率和機械功率的差值以及系統(tǒng)慣性中心功角加速度變化情況，圖2(a)中為DFIG電磁功率和機械功率的差值，有。

圖2 DFIG直接接入時的系統(tǒng)暫態(tài)響應

由圖2可見，故障期間，系統(tǒng)慣性中心的功角轉(zhuǎn)速處于加速狀態(tài)，DFIG電磁功率和機械功率的差值為負；系統(tǒng)慣性中心的功角轉(zhuǎn)速處于減速狀態(tài)，DFIG電磁功率和機械功率的差值為正。相對應的，根據(jù)式(13)，可以看到隨著DFIG接入容量的增加，的絕對值增大，而系統(tǒng)慣性中心的功角轉(zhuǎn)速加速度減小，功角轉(zhuǎn)速在故障期間的加速過程和故障后的減速過程都更為緩慢。

3.3 DFIG替代同步發(fā)電機組接入系統(tǒng)

DFIG風電場分別以替代同步發(fā)電機組GEN1和GEN3的方式接入系統(tǒng)，系統(tǒng)負荷保持不變：

CASE A：DFIG風電場等容量替代同步發(fā)電機組GEN1，配置合適容量的靜態(tài)無功補償器(static VAR compensator SVC)對DFIG并網(wǎng)點提供暫態(tài)電壓支撐，保證系統(tǒng)的暫態(tài)電壓穩(wěn)定性不改變。

CASE B：DFIG風電場等容量替代同步發(fā)電機組GEN3，同樣配置合適容量的SVC，保證系統(tǒng)的暫態(tài)電壓穩(wěn)定性不改變。

圖3和圖4分別展示了DFIG風電場替代同步發(fā)電機組GEN1和GEN3情況下的系統(tǒng)暫態(tài)響應和慣性中心功角加速度的暫態(tài)變化趨勢。

圖3 DFIG替代GEN1時的系統(tǒng)暫態(tài)響應

圖4 DFIG替代GEN3時的系統(tǒng)暫態(tài)響應

根據(jù)式(17)和式(18)可知，當DFIG替代同步發(fā)電機組接入系統(tǒng)時，系統(tǒng)慣性中心功角加速度的增減與否取決于式(17)中(/)–的正負號；變化幅度的大小則同時受影響于被替代同步發(fā)電機組的不同慣性時間常數(shù)。由圖3(a)、3(b)和圖4(a)、4(b)可以看到(/)、和(/)–的暫態(tài)變化趨勢：故障期間，兩組仿真的(/)–均小于零，此時系統(tǒng)慣性中心功角處于加速階段，因此增速加速度減小；故障切除后，系統(tǒng)慣性中心功角開始減速，在前0.3?s內(nèi)，CASE A的(/)–小于零，而CASE B的(/)–大于零，因此系統(tǒng)慣性中心功角的暫態(tài)響應分別有如圖3(c)和圖4(c)的變化趨勢。

4 ?結論

雙饋風電機組既不存在功角失穩(wěn)機制，也無法直接參與同步發(fā)電機組之間的互同步穩(wěn)定機制。本文基于系統(tǒng)慣性中心的動態(tài)數(shù)學模型，推導了雙饋風電機組接入后的系統(tǒng)慣性中心運動方程，并通過分析雙饋風電機組接入對系統(tǒng)慣性中心運動方程中關鍵因素的影響，揭示了雙饋風電機組以不同接入方式系統(tǒng)時，對系統(tǒng)慣性中心的暫態(tài)功角加速度的影響機理。

當雙饋風電機組直接接入系統(tǒng)時，系統(tǒng)慣性中心功角加速度的暫態(tài)變化趨勢主要取決于雙饋風電機組的接入容量和其暫態(tài)功率響應特性；當雙饋風電機組以等容替代同步發(fā)電機組的方式接入系統(tǒng)時，雙饋風電機組的暫態(tài)功率響應情況和被替代同步發(fā)電機組的功率和慣量在原系統(tǒng)中的比重大小將共同決定系統(tǒng)慣性中心功角加速度的暫態(tài)變化。本文工作對于進一步研究雙饋風電機組接入對區(qū)域系統(tǒng)間暫態(tài)功角穩(wěn)定性具有重要幫助。

[1] JIANG Liping, CHI Yongning, QIN Haiyan, et al. Wind energy in China[J]. IEEE Power and Energy Magazine, 2011, 9(6): 36-46.

[2] 王境彪, 希望?阿不都瓦依提, 晁勤, 等. 雙饋風電場接入交流或交直流混合系統(tǒng)低頻振蕩模態(tài)對比分析[J]. 電力系統(tǒng)保護與控制, 2015, 43(13): 23-29.

WANG Jingbiao, XI Wang?Abuduwayiti, CHAO Qin, et al. Modal contrastive analysis of the low-frequency oscillation of AC and AC/DC hybrid system integrated by DFIG-based wind farms[J]. Power System Protection and Control, 2015, 43(13): 23-29.

[3] 齊尚敏,李鳳婷,何世恩,等. 具有低電壓穿越能力的集群接入風電場故障特性仿真研究[J].電力系統(tǒng)保護與控制,2015,43(14): 55-62.

QI Shangmin, LI Fengting, HE Shien, et al. Research on fault characteristics of grid-connected large scale wind farms with LVRT capability[J]. Power System Protection and Control, 2015, 43(14): 55-62.

[4] PENA R, CLARE J C, ASHER G M. Doubly fed induction generator using back-to-back PWM converters and its application to variable-speed wind-energy generation[J]. Electric Power Applications, 1996, 143(3): 231-241.

[5] 馮雙磊, 趙海翔, 任普春, 等. 基于PSCAD/EMTDC的雙饋式變速恒頻風電機組動態(tài)模型仿真[J]. 電網(wǎng)技術, 2007, 31(17): 30-35.

FENG Shuanglei, ZHAO Haixiang, REN Puchun, et al. PSCAD/EMTDC based simulation study on dynamic model of doubly-fed variable speed wind turbine[J]. Power System Technology, 2007, 31(17): 30-35.

[6] 李巖松,郭世繁,任國威. 雙饋式風力發(fā)電系統(tǒng)的機電暫態(tài)建模與運行分析[J].電力系統(tǒng)保護與控制, 2013, 41(13): 101-109.

LI Yansong, GUO Shifan, REN Guowei. Modelling and operation simulation of doubly-fed type wind power generation system electromechanical transient model[J].Power System Protection and Control, 2013, 41(13): 101-109.

[7] 秦超,劉艷麗,余貽鑫,等.含雙饋風機電力系統(tǒng)的動態(tài)安全域[J]. 電工技術學報, 2015, 30(18): 157-163.

QIN Chao, LIU Yanli, YU Yixin, et al. Dynamic security region of power systems with double fed induction generator[J]. Transactions of China Electrotechnical Society, 2015, 30(18): 157-163.

[8] 王振樹, 劉巖, 雷鳴, 等. 基于Crowbar的雙饋機組風電場等值模型及并網(wǎng)仿真分析[J]. 電工技術學報, 2015, 30(4): 44-51.

WANG Zhenshu, LIU Yan, LEI Ming, et al. Doubly-fed induction generator wind farm aggregated model based on Crowbar and integration simulation analysis[J]. Transactions of China Electrotechnical Society, 2015, 30(4): 44-51.

[9] 林俐, 楊以涵. 基于擴展等面積定則的含大規(guī)模風電場電力系統(tǒng)暫態(tài)穩(wěn)定性分析[J]. 電力系統(tǒng)保護與控制, 2012, 40(12): 105-110.

LIN Li, YANG Yihan. Analysis of transient stability of power system including large scale wind power based on the extended equal area rule[J]. Power System Protection and Control, 2012, 40(12): 105-110.

[10] 孫勇, 陳晨, 王瑞明, 等. 電壓兩次跌落對雙饋型風電機組的影響分析[J]. 電網(wǎng)與清潔能源, 2013, 29(8): 109-112.

SUN Yong, CHEN Chen, WANG Ruiming, et al. Impacts on DFIG wind turbine during the voltage double-dip[J]. Power System and Clean Energy, 2013, 29(8): 109-112.

[11] 劉鋮, 蔡國偉, 楊德友, 等. 計及風機無功約束的電壓穩(wěn)定性分析[J]. 電力系統(tǒng)及其自動化學報, 2013, 25(1): 1-6.

LIU Cheng, CAI Guowei, YANG Deyou, et al. Analysis of voltage stability considering reactive power limitation of wind turbine[J]. Proceedings of the CSU-EPSA, 2013, 25(1): 1-6.

[12] 李東東, 梁自超, 周玉琦. 含風電場的受端系統(tǒng)暫態(tài)電壓穩(wěn)定性評估[J]. 電力系統(tǒng)保護與控制, 2015, 43(13): 8-14.

LI Dongdong, LIANG Zichao, ZHOU Yuqi. Transient voltage stability assessment of receiving end grid incorporating with wind farm[J]. Power System Protection and Control, 2015, 43(13): 8-14.

[13] HAO Zhenghang, YAO Zhiqing, LI Shaohua, et a1. The contribution of double-fed wind farms to transient voltage and damping of power grids[J]. Tehnicki Vjesnik-Technical Gazette, 2015, 22(1): 43-49.

[14] 時銳, 焦然, 陳艷霞. 含大規(guī)模風電電網(wǎng)的靜態(tài)電壓仿真分析[J]. 高壓電器, 2014, 50(12): 105-109.

SHI Rui, JIAO Ran, CHEN Yanxia. Static voltage simulation analysis of grid with large-scale wind power[J]. High Voltage Apparatus, 2014, 50(12): 105-109.

[15] LIN Li, HOU Hongjing, Li Wenxia, et al. Simulation and comparison of transient stability of power system includingDFIGs wind farm based on detailed model[C]?// Sustainable Power Generation and Supply, Nanjing, 2009: 1-6.

[16] SHI Libao, DAI Shiqiang, NI Yixin, et al. Transient stability of power systems with high penetration of DFIG based wind farms[C] // Power & Society General Meeting, Calgary, 2009: 1-6.

[17] 郝正航, 余貽鑫, 曾沅. 雙饋風力發(fā)電機功角暫態(tài)行為及其控制策略[J]. 電力自動化設備, 2011, 31(2): 79-83.

HAO Zhenghang, YU Yixin, ZENG Yuan. Transient performance of DFIG power angle in wind farm and its control strategy[J]. Electric Power Automation Equipment, 2011, 31(2): 79-83.

[18] CHOWDHURY M A, HOSSEINZADEH N, SHEN W X, et al. Comparative study on fault responses of synchronous generators and wind turbine generators using transient stability index based on transient energy function[J]. International Journal of Electrical Power & Energy Systems, 2013, 51(10): 145-152.

[19] KUNDUR P. Power system stability and control[M]. New York, USA: Mc Graw Hill, Inc, 1992: 461.

[20] 倪以信, 陳壽孫, 張寶霖. 動態(tài)電力系統(tǒng)的理論和分析[M]. 北京:清華大學出版社, 2002: 186.

(編輯 葛艷娜)

Impact analysis of doubly-fed induction generator on the transient angle stability of the region with wind power integrated

LIU Siwei, LI Gengyin, ZHOU Ming

(State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources (North China Electric Power University), Beijing 102206, China)

The dynamic motion of system center of inertia (COI) can effectively indicate the general trend of a power system dynamics. Based on this, the developed COI dynamic model considering integrated doubly-fed induction generators (DFIGs) is proposed. With the impacts of DFIGs on the key factors in COI dynamic motions deduced, the influence mechanism of DFIGs on the power angle accelerating speed (PAAS) of system COI is then deeply investigated. If DFIGs are integrated without displacing synchronous generators, the change in PAAS of system COI is mainly determined by the capacity and transient power response of integrated DFIGs. And when synchronous generators are displaced, the change in PAAS of system COI is affected jointly by the DFIG transient power response and the proportions of the inertia and capacity of displaced synchronous generators in the origin system.The theoretical analysis is effectively verified and illustrated by simulation studies. This work is supported by the Major Program of the National Natural Science Foundation of China (No. 51190103) and the “111” Project (No. B08013).

doubly-fed induction generator (DFIG); system center of inertia (COI); transient angle stability; power angle accelerating speed; influence mechanism

10.7667/PSPC151725

國家自然科學基金重大項目課題(51190103)；高等學校學科創(chuàng)新引智計劃項目(B08013)

2015-09-25；

2015-11-17

劉斯偉(1988-)，男，博士，研究方向為雙饋風電機組的建模與控制、風電接入對系統(tǒng)暫態(tài)穩(wěn)定性的影響分析；E-mail: liusw@necpu.edu.cn 李庚銀(1964-)，男，博士，教授，博士生導師，主要研究方向為新能源電力系統(tǒng)分析與控制，先進輸變電技術，電能質(zhì)量等；E-mail: ligy@ncepu.edu.cn 周 明(1967-)，女，博士，教授，博士生導師，主要研究方向為新能源電力系統(tǒng)分析與控制，電力市場，電能質(zhì)量等。E-mail: zhouming@ncepu.edu.cn