郭權(quán)利，苑舜

（1.沈陽工業(yè)大學(xué)電氣工程學(xué)院，遼寧 沈陽110023 ;2.電網(wǎng)節(jié)能與控制遼寧省重點(diǎn)實(shí)驗(yàn)室，沈陽工程學(xué)院電氣工程系，遼寧 沈陽 110136;3.國家電監(jiān)會(huì)東北監(jiān)管局，遼寧 沈陽110006）

0 引言

目前，微電網(wǎng)中絕大部分采用電力電子器件作為微電源與電網(wǎng)的接口,微電源種類不同，應(yīng)用接口方式也不同[1-4]。因此，在Matlab 仿真軟件中，需要建立不同接口仿真模型。即使同一個(gè)接口在Matlab 不同的運(yùn)行模式下，其接口仿真模型也不同[5-6]。Matlab/Simulink 所提供的接口模型較少，接口模型己成為阻礙使用Matlab/Simulink 進(jìn)行微電網(wǎng)仿真研究的主要原因。如何處理接口已經(jīng)是仿真的核心問題。接口控制部分大多采用矢量控制技術(shù)[7-8]，所以本文對(duì)矢量控制技術(shù)進(jìn)行了研究，并設(shè)計(jì)出基于矢量控制技術(shù)的適合于不同運(yùn)行模式下的微電網(wǎng)接口模型，以方便用戶在Matlab/Simulink 下進(jìn)行電力系統(tǒng)及微電網(wǎng)的建模與仿真。

1 基于矢量控制的微電網(wǎng)接口模型的建立

矢量控制技術(shù)是隨著電力電子技術(shù)的發(fā)展而出現(xiàn)的新技術(shù)。最初是為了滿足現(xiàn)代交流調(diào)速系統(tǒng)較高的動(dòng)態(tài)要求，利用旋轉(zhuǎn)坐標(biāo)變換，采用磁場(chǎng)定向控制技術(shù)，將電動(dòng)機(jī)的定子電流分解成勵(lì)磁電流和轉(zhuǎn)矩電流，簡(jiǎn)化了交流電動(dòng)機(jī)的動(dòng)態(tài)數(shù)學(xué)模型，使交流電機(jī)變得像直流電機(jī)一樣易于控制[9]。其技術(shù)的核心是采用合理的坐標(biāo)變換。本文借鑒矢量控制思想，利用不同的坐標(biāo)變換，采用電壓定向控制技術(shù)，設(shè)計(jì)了三種應(yīng)用于不同運(yùn)行模式下的微電源與電力系統(tǒng)接口模型。按照電壓定向所在的坐標(biāo)系的不同，本文將其分為基于dq0 、αβ、旋轉(zhuǎn)abc 坐標(biāo)系下的微電網(wǎng)接口模型。所有的接口模型均采用等效原則進(jìn)行坐標(biāo)變換，即：相同的接口在不同坐標(biāo)系下的模型中，其電壓、電流、磁鏈等合成矢量不變、傳輸?shù)墓β什蛔儭?/p>

1.1 微電網(wǎng)接口

圖1 為典型的微電源接口——電壓型換流器[10]。圖中：uca、ucb、ucc和usa、usb、usc分別為換流器三相電壓和微電網(wǎng)的三相電壓；ia、ib、ic為由換流器側(cè)向網(wǎng)側(cè)逆變電流；C為換流側(cè)電容；idl為微電源提供的直流電流；idc為換流器側(cè)向微電網(wǎng)轉(zhuǎn)換的直流電流。

圖1 換流器結(jié)構(gòu)圖Fig.1 Converter structure

其狀態(tài)方程為

本文以圖1 為例建立其分別基于dq0、αβ、旋轉(zhuǎn)abc 坐標(biāo)系下的微電網(wǎng)接口模型。

1.2 dq0 坐標(biāo)系下的接口模型

首先建立dq0 坐標(biāo)系。坐標(biāo)系的建立分三步進(jìn)行，第一步：將靜止的abc 坐標(biāo)系統(tǒng)轉(zhuǎn)換為靜止的αβ0 坐標(biāo)系統(tǒng)，第二步：將αβ0 坐標(biāo)系統(tǒng)轉(zhuǎn)換為旋轉(zhuǎn)的dq0 坐標(biāo)系統(tǒng)。第三步：重新定向dq0 坐標(biāo)系統(tǒng)。

abc 系統(tǒng)向αβ0 的系統(tǒng)轉(zhuǎn)換所選擇的坐標(biāo)系如圖2(a)所示。坐標(biāo)變換關(guān)系為

αβ0 系統(tǒng)向旋轉(zhuǎn)的dq0 系統(tǒng)的轉(zhuǎn)換所選擇的坐標(biāo)系如圖2(b)所示，d、q 軸和合成矢量F都以電網(wǎng)角速度sω旋轉(zhuǎn)，分量ud、uq的長短不變。α、β 軸是靜止的，α 軸與d 軸的夾角?隨時(shí)間而變化，因此F在α、β 軸上的分量的長短也隨時(shí)間變化，相當(dāng)于交流電壓的瞬時(shí)值。顯然，uα、uβ、u0和ud、uq、u0間存在下列關(guān)系：

圖2 坐標(biāo)變換Fig.2 Coordinate transformation

為了實(shí)現(xiàn)電網(wǎng)與微電源之間傳輸?shù)挠泄β屎蜔o功功率的解耦控制，微電源與電力系統(tǒng)接口通常采用基于電網(wǎng)電壓定向的矢量控制方案，使電網(wǎng)電壓的合成矢量方向與新的dq0坐標(biāo)系d軸方向相同。在圖2(b)中，將合成矢量F方向定義為新的d 軸，逆時(shí)針超前d 軸90°的方向?yàn)閝 軸。至此完成了dq0坐標(biāo)系的建立。

根據(jù)式（3）、式（4）將式（1）中的iabc、ucabc、usabc轉(zhuǎn)化到圖2（b）坐標(biāo)系下，其接口狀態(tài)方程可變?yōu)?/p>

顯然，在αβ0 系統(tǒng)向旋轉(zhuǎn)的dq0 系統(tǒng)的轉(zhuǎn)換過程中，需要較多地測(cè)量或計(jì)算轉(zhuǎn)換矩陣式（4）中的?值，也即需要使用鎖相技術(shù)，給仿真及控制系統(tǒng)的設(shè)計(jì)帶來不便，需要進(jìn)一步簡(jiǎn)化模型，可采用本文的第二種模型。

1.3 αβ坐標(biāo)系下的微電網(wǎng)接口模型

首先建立αβ坐標(biāo)系，αβ 坐標(biāo)系的建立分兩步進(jìn)行，第一步：將靜止的abc 坐標(biāo)系統(tǒng)轉(zhuǎn)換為靜止αβ的坐標(biāo)系統(tǒng)，并求出合成矢量；第二步：定向dq 坐標(biāo)。

abc 系統(tǒng)向αβ的系統(tǒng)轉(zhuǎn)換同樣選取如圖2(a)所示坐標(biāo)系。按照坐標(biāo)等效原則，轉(zhuǎn)換矩陣為

合成矢量按下式計(jì)算：

按照微電源與電力系統(tǒng)接口所采用基于電網(wǎng)電壓定向的矢量控制方案，選取合成矢量方向F為d軸，逆時(shí)針超前d 軸90°的方向?yàn)閝 軸。如圖2（c）所示。則αβu在dq 坐標(biāo)系下的坐標(biāo)udq為

根據(jù)式（6）、式（8）將式（1）中的iabc、ucabc、usabc轉(zhuǎn)化到圖2（c）坐標(biāo)系下，其方程可變?yōu)?/p>

微電源向微電網(wǎng)提供的三相電流為

式（9）、式（10）構(gòu)成了接口電路的狀態(tài)方程。

在dq 坐標(biāo)系下，微電源經(jīng)過換流器向微電網(wǎng)提供的功率可表示為

由式（11）可見，在電網(wǎng)電壓保持恒定時(shí)，通過設(shè)定id就可以調(diào)節(jié)微電源經(jīng)換流器轉(zhuǎn)換到微電網(wǎng)側(cè)的有功功率，通過設(shè)定iq就可以調(diào)節(jié)微電源經(jīng)換流器轉(zhuǎn)換到微電網(wǎng)側(cè)的無功功率。從而實(shí)現(xiàn)微電源與電網(wǎng)之間傳輸?shù)挠泄β屎蜔o功功率的解耦控制。經(jīng)式（6）～式(10)可以建立靈活多變的不同控制方式的接口模型。圖3(a)為V/F 控制的接口模型示意圖。具體原理可參閱文獻(xiàn)[10-12]。

1.4 abc 坐標(biāo)系下的微電網(wǎng)接口模型

abc 坐標(biāo)系的建立分兩步進(jìn)行，第一步：將靜止的ABC 坐標(biāo)系統(tǒng)轉(zhuǎn)換為旋轉(zhuǎn)的abc 的坐標(biāo)系統(tǒng)，并求出合成矢量；第二步：定向dq 坐標(biāo)。

圖3 接口模型Fig.3 Interface mode

靜止的ABC系統(tǒng)向旋轉(zhuǎn)的abc系統(tǒng)的轉(zhuǎn)換所選擇的坐標(biāo)系如圖4（a）所示。旋轉(zhuǎn)坐標(biāo)系的a 軸參考起始位置與靜止坐標(biāo)系A(chǔ) 軸重合。

圖4 坐標(biāo)變換Fig.4 Coordinate transformation

a、b、c 軸和合成矢量F都以電網(wǎng)角速度sω旋轉(zhuǎn)，分量ua、ub、uc的長短不變。 A、B、C 軸是靜止的，靜止的坐標(biāo)軸ABC 與旋轉(zhuǎn)坐標(biāo)的軸abc的夾角?隨時(shí)間而變化，顯然，ua、ub、uc和uA、uB、uC之間存在下列關(guān)系

合成矢量按下式計(jì)算：

按照微電源與電力系統(tǒng)接口所采用基于電網(wǎng)電壓定向的矢量控制方案，選取合成矢量方向F為d軸，逆時(shí)針超前d 軸90°的方向?yàn)閝 軸。如圖4(b)所示。則合成矢量在dq 坐標(biāo)系下的坐標(biāo)udq為

由于Matlab 中的測(cè)量模塊已經(jīng)完成了第一步（式12）的坐標(biāo)轉(zhuǎn)換，所以只需根據(jù)式（13）、式（14）將式（1）中的iabc、ucabc、usabc轉(zhuǎn)化到圖4（b）坐標(biāo)系下即可，其方程可變?yōu)?/p>

則微電源向微電網(wǎng)提供的三相電流為

由于坐標(biāo)變換采用了電網(wǎng)電壓定向的矢量控制方案，同樣可以實(shí)現(xiàn)微電源與電網(wǎng)之間傳輸?shù)挠泄β屎蜔o功功率的解耦控制。向微電網(wǎng)提供的功率的表達(dá)式與式（11）相同。經(jīng)式（11）～式(16)可以建立接口模型。圖3(b)為P/Q 控制的接口模型。

2 仿真分析

圖5 為微電網(wǎng)仿真系統(tǒng)示意圖。微電網(wǎng)由配電網(wǎng)、微型燃?xì)廨啓C(jī)[11-14]、光伏發(fā)電[15-18]、變壓器、輸電線路及一般負(fù)荷組成，具體參數(shù)如圖中所示。仿真過程采用標(biāo)幺值，基準(zhǔn)容量S=100 kVA 基準(zhǔn)電壓V=400 V，基準(zhǔn)頻率f=50 Hz。

Matlab/Simulink 中的powersystem 模塊集分為三種運(yùn)行模式：連續(xù)運(yùn)行模式、離散運(yùn)行模式、矢量運(yùn)行模式。本例分別在離散運(yùn)行模式及矢量運(yùn)行模式下進(jìn)行仿真。

圖5 微電網(wǎng)仿真系統(tǒng)圖Fig.5 Microgrid simulation system diagram

2.1 離散運(yùn)行模式下的仿真

在此模式下，微電網(wǎng)采用基于dq0 或αβ坐標(biāo)系下的接口模型進(jìn)行微電網(wǎng)與大電網(wǎng)并網(wǎng)運(yùn)行及并網(wǎng)、孤島轉(zhuǎn)換的仿真。并網(wǎng)運(yùn)行時(shí)，微型燃?xì)廨啓C(jī)采用P/Q 控制方式，在微電網(wǎng)孤島運(yùn)行時(shí)采用V/F控制方式。

微電網(wǎng)在0.0～3.0 s 時(shí)并網(wǎng)運(yùn)行，3.0～5.0 s，孤島運(yùn)行。0.5 s，燃?xì)廨啓C(jī)有功功率由初始給定0.1升至0.5，無功功率給定由0 變?yōu)?.2；1.5 s 光照強(qiáng)度由最初的500 W/m2變?yōu)?00 W/m2，光伏發(fā)電的無功功率給定與初始給定（0）相同；在t=3 s 時(shí)斷路器斷開，微電網(wǎng)變?yōu)楣聧u運(yùn)行。有功功率、無功功率仿真結(jié)果如圖6（a）；母線電壓的幅值及微電網(wǎng)的頻率在微電網(wǎng)運(yùn)行方式轉(zhuǎn)換時(shí)的變化情況如圖6(b)所示。

圖6 離散運(yùn)行模式下的仿真圖形Fig.6 Simulation of the discrete mode

2.2 矢量運(yùn)行模式下的仿真

在此模式下，微電網(wǎng)采用基于旋轉(zhuǎn)abc 坐標(biāo)系下的接口模型進(jìn)行微電網(wǎng)與大電網(wǎng)并網(wǎng)運(yùn)行及配電線路故障的仿真。并網(wǎng)運(yùn)行仿真參數(shù)及過程與離散運(yùn)行模式下的并網(wǎng)仿真相同，其有功功率、無功功率仿真結(jié)果如圖7(a)。在t=3 s 線路發(fā)生三相短路故障，持續(xù)時(shí)間為0.3 s，仿真結(jié)果如圖7(b)所示。

由仿真結(jié)果可以看出，并網(wǎng)運(yùn)行時(shí)在兩種運(yùn)行模式下的仿真結(jié)果相同；在微電網(wǎng)并網(wǎng)、孤網(wǎng)轉(zhuǎn)換過程中，仿真結(jié)果顯示了微電網(wǎng)維持頻率及電壓穩(wěn)定的過程；配電線路發(fā)生故障時(shí)，仿真結(jié)果也反映出了微電網(wǎng)電壓、電流變化。

本仿真中雖然只包含太陽能、燃?xì)廨啓C(jī)兩種微電源，但其接口模型同樣可以應(yīng)用于風(fēng)機(jī)、燃料電池、蓄電池組、儲(chǔ)能電容等微電源與微電網(wǎng)的連接。因此采用本文設(shè)計(jì)的基于矢量控制的接口模型連接組成的微電網(wǎng)進(jìn)行仿真，能夠正確反映出微電網(wǎng)在各種擾動(dòng)下的電壓、頻率、有功、無功等變化，基于矢量控制的接口模型能夠滿足微電網(wǎng)仿真要求。

圖7 矢量運(yùn)行模式下的仿真圖形Fig.7 Simulation of the phasor mode

3 結(jié)論

1）本文設(shè)計(jì)的矢量控制的接口模型，解決了微電源與微電網(wǎng)在不同模塊集、不同運(yùn)行模式下組網(wǎng)仿真的接口問題。

2）基于矢量控制的接口模型，不論在何種運(yùn)行模式、何種坐標(biāo)系下，均采用電網(wǎng)電壓定向的矢量控制方案，使電壓向量在dq 坐標(biāo)系下q 軸分量為零，實(shí)現(xiàn)微電源與電網(wǎng)之間傳輸?shù)挠泄β屎蜔o功功率的解耦控制。從而使控制變量的物理含義十分清晰，簡(jiǎn)化了接口模型。

3）從建立矢量控制的微電源接口模型的過程中可以看出，模型略去了接口中開關(guān)的暫態(tài)過程，只反映了接口兩端的轉(zhuǎn)換關(guān)系，從而極大地加快了仿真速度。

4）矢量控制技術(shù)應(yīng)用較為廣泛，基于矢量控制的接口模型的建立不僅豐富了Matlab 資源庫，而且方便了用戶在Matlab/Simulink下進(jìn)行其他方面的建模與仿真。

[1]高明遠(yuǎn).雙向 DC-DC 變換器基于切換系統(tǒng)的建模與儲(chǔ)能控制[J].電力系統(tǒng)保護(hù)與控制,2012,40(3):129- 134.

GAO Ming-yuan.Modeling and energy storage control for bi-directional DC-DC converter based on switching system[J].Power System Protection and Control,2012,40(3):129-134.

[2]王陽,魯宗相,閡勇.微電網(wǎng)中微電源接口仿真模型的研究與比較[J].電力系統(tǒng)自動(dòng)化,2010,34(1):84-88.

WANG Yang,LU Zong-xiang,MIN Yong.Research and comparison of DG interface simulation models in a microgrid[J].Automation of Electric Power Systems,2010,34(1):84-88.

[3]Suryanarayanan S,Mitra J.Enabling technologies for the customer-driven microgrid[C]// Proceedings of 2009 IEEE PES General Meeting,Calgary,Canada,July 28,2009.

[4]XIONG Xiao-fu,OUYANG Jin-xin.Modeling and transient behavior analysis of an inverter-based microgrid[J].Electric Power Components and Systems,2011,40(1):112-130.

[5]張琨,高思超,畢靖.MATLAB 2010 從入門到精通[M].北京:電子工業(yè)出版社,2011.

[6]薛定宇,陳陽泉.基于MATLAB/Simulink 的系統(tǒng)仿真技術(shù)與應(yīng)用[M].北京:清華大學(xué)出版社,2011.

[7]陳海榮,徐政.基于同步坐標(biāo)變換的 VSC-HVDC 暫態(tài)模型及其控制器[J].電工技術(shù)學(xué)報(bào),2007,22(2):121-126.

CHEN Hai-rong,XU Zheng.Transient model and controller design for VSC-HVDC based on synchronous reference frame[J].Transactions of China Electrotechnical Society,2007,22(2):121-126.

[8]Chung I,Park S,Kim H,et al.Operating strategy andcontrol scheme of premium power supply interconnected with electric power systems[J].IEEE Trans on Power Delivery,2005,20(3):2281-2288.

[9]羅飛,郗曉田,文小玲,等.電力拖動(dòng)與運(yùn)動(dòng)控制系統(tǒng)[M].北京:化學(xué)工業(yè)出版社,2011.

[10]陳謙,陶國慶,胡銘.采用dq0 坐標(biāo)的VSC-HVDC 穩(wěn)態(tài)模型與控制器設(shè)計(jì)[J].電力系統(tǒng)自動(dòng)化,2004,28(16):61-66.

CHEN Qian,TAO Guo-qing,HU Ming.Steady-state model and controller design of a VSC-HVDC converter based on dq0-axis[J].Automation of Electric Power Systems,2004,28(16):61-66.

[11]張建華,黃偉.微電網(wǎng)運(yùn)行控制與保護(hù)技術(shù)[M].北京:中國電力出版社,2010.

[12]Sao C K,Lehn P W.Control and power management of converter fed microgrids[J].IEEE Trans on Power Systems,2008,23(3):1088-1098.

[13]Sao C K,Lehn P W.Intentional islanded operation of converter fed microgrids[C]// IEEE Power Eng Soc General Meeting,2006:18-22.

[14]劉君,穆世霞,李巖松,等.微電網(wǎng)中微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)整體建模與仿真[J].電力系統(tǒng)自動(dòng)化,2010,34(7):85-89.

LIU Jun,MU Shi-xia,LI Yan-song,et al.Overall modeling and simulation of microturbine generation systems in microgrids[J].Automation of Electric Power Systems,2010,34(7):85-89.

[15]才利存,常忠廷.太陽能并網(wǎng)發(fā)電系統(tǒng)的 MPPT 及孤島檢測(cè)技術(shù)[J].電力系統(tǒng)保護(hù)與控制,2012,40(4):119-124.

CAI Li-cun,CHANG Zhong-ting.MPPT and anti-islanding of photovoltaic grid-connected system[J].Power System Protection and Control,2012,40(4):119-124.

[16]邱培春,葛寶明,畢大強(qiáng).基于擾動(dòng)觀察和二次插值的光伏發(fā)電最大功率跟蹤控制[J].電力系統(tǒng)保護(hù)與控制,2011,39(4):62-67.

QIU Pei-chun,GE Bao-ming,BI Da-qiang.MPPT control for PV power generation system based on P&O algorithms and quadratic interpolation[J].Power System Protection and Control,2011,39(4):62-67.

[17]李晶,許洪華,趙海翔,等.并網(wǎng)光伏電站動(dòng)態(tài)建模及仿真分析[J].電力系統(tǒng)自動(dòng)化,2008,32(24):83-87.

LI Jing,XU Hong-hua,ZHAO Hai-xiang,et al.Dynamic modeling and simulation of the grid-connected PV power station[J].Automation of Electric Power Systems,2008,32(24):83-87.

[18]廖志凌,阮新波.任意光強(qiáng)和溫度下的硅太陽電池非線性工程簡(jiǎn)化數(shù)學(xué)模型[J].太陽能學(xué)報(bào),2009,30(4):430-435.

LIAO Zhi-ling,RUAN Xin-bo.Non-linear engineering simplification model of silicon solar cells in arbitrary solar radiation and temperature[J].Acta Energiae Solaris Sinica,2009,30(4):430-435.