句榮濱，宋旭日，鐘麗波，趙曉娜，鄭 璐

(1.國網(wǎng)遼寧省電力有限公司，遼寧 沈陽 110006；2.中國電力科學(xué)研究院，北京 100192； 3.國網(wǎng)沈陽供電公司，遼寧 沈陽 110003；4.國網(wǎng)遼寧省電力有限公司技能培訓(xùn)中心，遼寧 錦州 121003； 5.國網(wǎng)葫蘆島供電公司,遼寧 葫蘆島 125003)

微電網(wǎng)是具有規(guī)模小、分散式特點(diǎn)的獨(dú)立系統(tǒng)，利用豐富而先進(jìn)的現(xiàn)代電力技術(shù)，將分布式電源(燃?xì)廨啓C(jī)、風(fēng)機(jī)、光伏發(fā)電、燃料電池等)儲(chǔ)能及保護(hù)系統(tǒng)等裝置整合在一起，直接與用戶側(cè)相接。微電網(wǎng)作為我國社會(huì)經(jīng)濟(jì)發(fā)展的基礎(chǔ)設(shè)施和重要的公用事業(yè)，充分發(fā)揮其資源優(yōu)化配置，滿足了電力短缺地區(qū)用戶的特定需求，有效提高了用戶側(cè)的電力特性，例如在提高本地供電可靠性的同時(shí)降低了饋線損耗；在保持本地電壓穩(wěn)定性和提供不間斷電源的同時(shí)提高了余熱能量利用效率等[1-4]。

目前微電網(wǎng)控制特性、能量管理及控制是微電網(wǎng)集中性研究范疇，而對(duì)于微電網(wǎng)動(dòng)態(tài)運(yùn)行特性缺少更加專業(yè)且深入地發(fā)掘及研究。文獻(xiàn)[5-6]首先在不同微電源的微電網(wǎng)動(dòng)態(tài)仿真結(jié)果的基礎(chǔ)上，分析研究不同微電源配置方案對(duì)微電網(wǎng)運(yùn)行模式轉(zhuǎn)換動(dòng)態(tài)運(yùn)行特性的影響，但沒有對(duì)不同微電源的特性進(jìn)行對(duì)比研究分析；文獻(xiàn)[7]建立了風(fēng)/光/柴互補(bǔ)的微電網(wǎng)，對(duì)故障情況下微電網(wǎng)故障恢復(fù)能力進(jìn)行了仿真驗(yàn)證，但沒有對(duì)抑制微電網(wǎng)暫態(tài)特性的措施進(jìn)行研究；文獻(xiàn)[8]考慮了風(fēng)/光/燃?xì)廨啓C(jī)等微電源的動(dòng)態(tài)特性，對(duì)微電網(wǎng)運(yùn)行模式切換進(jìn)行了仿真，但忽略了模式切換過程中的暫態(tài)特性；文獻(xiàn)[9-10]闡述了基于IEC60909標(biāo)準(zhǔn)的分布式電源短路計(jì)算，用交/直流電壓源代替分布式電源，但忽略了各微電源的外在特性。

風(fēng)能和太陽能具有天然互補(bǔ)性，因此風(fēng)能與太陽能混合微電網(wǎng)可以向用戶提供更加穩(wěn)定的電能，從而提高微電網(wǎng)的供電可靠性。本文在電磁暫態(tài)仿真軟件EMTP/ATP搭建太陽能電池模型、永磁同步風(fēng)力發(fā)電機(jī)模型、負(fù)荷模型，并利用軟件自帶的變壓器、線路等一起構(gòu)成簡單的微電網(wǎng)模型。根據(jù)IEC60909標(biāo)準(zhǔn)對(duì)微電網(wǎng)暫態(tài)運(yùn)行特性(包括三相接地短路、單相接地短路等)進(jìn)行了仿真研究，結(jié)果表明本文建立的微電源及微電網(wǎng)模型能夠很好的模擬實(shí)際運(yùn)行情況，為實(shí)際工程計(jì)算提供了方法和手段。

1 IEC60909標(biāo)準(zhǔn)概述

IEC60909標(biāo)準(zhǔn)是2001年由國際電工委員會(huì)IEC(International Electrotechnical Commission)制定的新的關(guān)于短路計(jì)算的國際標(biāo)準(zhǔn)，標(biāo)準(zhǔn)明確指出三相交流系統(tǒng)短路電流計(jì)算使用的元件數(shù)學(xué)模型、計(jì)算方法、計(jì)算假定條件等，可以為暫態(tài)穩(wěn)定計(jì)算提供支持。

IEC60909標(biāo)準(zhǔn)規(guī)定電力系統(tǒng)內(nèi)部某點(diǎn)發(fā)生短路故障時(shí)，在該點(diǎn)引入1個(gè)虛擬電壓源作為網(wǎng)絡(luò)的唯一電壓源。其他電源電勢為零，并用內(nèi)阻抗代替，包括同步發(fā)電機(jī)、同步電動(dòng)機(jī)、異步電動(dòng)機(jī)和其他饋電網(wǎng)絡(luò)的電勢等。短路電流計(jì)算如圖1所示。

對(duì)稱三相短路初始電流計(jì)算為

(1)

IEC 60909對(duì)于各元件通用模型處理方法基本是以元件銘牌參數(shù)作為變量，以相應(yīng)系數(shù)進(jìn)行計(jì)算修正。通過不同系數(shù)值的變量變換，得出模型系統(tǒng)最大及最小短路電流值。各元件數(shù)學(xué)模型主要分析如下。

a.同步發(fā)電機(jī)(適用于永磁同步風(fēng)力發(fā)電機(jī))

系統(tǒng)發(fā)生短路時(shí)，同步發(fā)電機(jī)存在1個(gè)暫態(tài)過程，尤其靠近發(fā)電機(jī)機(jī)端短路瞬間發(fā)電機(jī)會(huì)產(chǎn)生較大的次暫態(tài)電勢，該暫態(tài)過程不可以忽視。假定短路瞬間E″的縱分量保持不變，如圖2所示，可以算出發(fā)電機(jī)短路阻抗修正公式，見式(2)、式(3)。

(2)

(3)

b.異步電動(dòng)機(jī)(適用于太陽能電池)

IEC60909規(guī)定通過DC/AC并網(wǎng)的微電源，忽略其原動(dòng)機(jī)類型，可將其等效為異步電動(dòng)機(jī)，等效電路為接地阻抗，短路阻抗值為

(4)

式中：UrM、IrM和SrM分別為異步電動(dòng)機(jī)的額定電壓、額定電流和額定視在功率；ILR/IrM為異步電動(dòng)機(jī)的堵轉(zhuǎn)電流與額定電流比值，一般取3。

c.變壓器

IEC60909標(biāo)準(zhǔn)規(guī)定變壓器短路阻抗可以根據(jù)式(5)進(jìn)行修正。

(5)

式中：xT為短路電抗，對(duì)于三繞組變壓器，把xT改為xAB，xBC，xAC。

2 微電源模型

2.1 永磁同步風(fēng)力發(fā)電機(jī)模型

a.風(fēng)輪機(jī)模型

風(fēng)輪機(jī)簡化數(shù)學(xué)模型為

(6)

式中：Pg為風(fēng)輪機(jī)獲得風(fēng)能的功率，W;ρ為空氣密度(標(biāo)準(zhǔn)狀況為1.29 kg/m3)；A為掃風(fēng)面積，m2；R為風(fēng)輪機(jī)葉片半徑，m；Cp為風(fēng)能利用系數(shù)；Vw為風(fēng)速，m/s。

風(fēng)輪機(jī)捕獲風(fēng)能的大小是由風(fēng)能利用系數(shù)CP決定，CP是槳距角β和葉尖速比λ的函數(shù)，葉尖速比λ計(jì)算為

(7)

式中：ω為風(fēng)輪角速度,rad/s；R為風(fēng)輪半徑,m；Vw為風(fēng)速,m/s。

b.永磁同步風(fēng)力發(fā)電機(jī)模型

永磁同步風(fēng)力發(fā)電機(jī)仿真模型如圖3所示[11]。

永磁同步風(fēng)力發(fā)電機(jī)在dq同步旋轉(zhuǎn)坐標(biāo)系下電壓方程為

(8)

式中：id和iq分別為發(fā)電機(jī)的d軸和q軸電流；Ld和Lq分別為發(fā)電機(jī)的d軸和q軸的電感；Ra為定子電阻；ωe為電角頻率，ωe=npωg；np為發(fā)電機(jī)轉(zhuǎn)子極對(duì)數(shù)；λ0為永磁鐵磁鏈，λ0=0.175 Wb；ud和uq分別為發(fā)電機(jī)出口電壓的d軸和q軸分量。

永磁同步風(fēng)力發(fā)電機(jī)電磁轉(zhuǎn)矩為

Te=1.5np[(Ld-Lq)idiq+iqλ0]

(9)

永磁同步風(fēng)力發(fā)電機(jī)傳動(dòng)模型為

(10)

式中：ω為風(fēng)輪角速度,rad/s；R為風(fēng)輪半徑,m；Vw為風(fēng)速,m/s。

2.2 太陽能發(fā)電單元模型

2.2.1 太陽能電池?cái)?shù)學(xué)模型

本文采用的太陽能電池?cái)?shù)學(xué)模型為工程用模型，模型計(jì)算技術(shù)參數(shù)為工程技術(shù)參數(shù)，包括短路電流Isc、開路電壓Uoc、最大功率點(diǎn)電流Im、最大功率點(diǎn)電壓Um、最大功率點(diǎn)功率Pm，從而可以在一定精度的計(jì)算模型基礎(chǔ)上，研究分析太陽能電池特性。太陽能電池模型描述如下[12-13]。

a.參考條件下，輸出電壓為Uo、對(duì)應(yīng)電流為I時(shí)，太陽能電池的I-U特性方程為

(11)

在最大功率點(diǎn)處，式(11)仍然成立。

式中:參數(shù)C1、C2分別為

(12)

(13)

式中:Np為光伏電池并聯(lián)數(shù);Ns為光伏電池串聯(lián)數(shù)。

b.考慮不同光照強(qiáng)度、不同溫度的太陽能電池輸出特性變化

首先通過已知的參考日照強(qiáng)度和參考電池溫度的Isc、Uoc、Im、Um計(jì)算出新的日照強(qiáng)度和電池溫度的I＇sc、U＇oc、I＇m、U＇m,然后代入實(shí)用表達(dá)式得到新的太陽能電池的I-U特性曲線。

(14)

式中：Sref為參考太陽能輻射強(qiáng)度(1000 W/m2)；Tref為參考電池溫度(25 ℃)；ΔT為實(shí)際電池溫度與參考電池溫度的差值；ΔS為實(shí)際光強(qiáng)與參考光強(qiáng)差值；太陽光強(qiáng)為S和電池溫度為T時(shí)，I＇sc為太陽能電池短路電流；U＇oc為開路電壓；I＇m為最大功率點(diǎn)電流；U＇m為最大功率點(diǎn)電壓；a、b、c的典型值為a=0.0025/℃、b=0.5/℃、

c=0.00 288/℃。

2.2.2 太陽能電池控制策略

太陽能電池輸出電流為直流，輸出電壓受溫度、光強(qiáng)等因素影響，所以輸出頻率并不穩(wěn)定，因此需要通過整流逆變裝置進(jìn)行電路調(diào)節(jié)后對(duì)負(fù)載供電，從而提高穩(wěn)定的供電特性。太陽能電池并網(wǎng)控制如圖4所示。

太陽能電池采用下垂控制，高壓系統(tǒng)的功率因數(shù)受控制頻率影響，因此可以通過改變逆變器輸出電壓頻率控制有功功率的動(dòng)態(tài)表現(xiàn)，同理控制逆變器輸出電壓幅值可以控制其無功功率的流動(dòng)。

3 風(fēng)/光互補(bǔ)微電網(wǎng)電磁暫態(tài)模型

微電網(wǎng)仿真系統(tǒng)電壓等級(jí)為10 kV；變壓器容量為400 kVA；變比為10 kV/400 V；輸電線路為等值模型，其中Lc=0.264 mH/km,Rc=0.28 Ω/km，為方便計(jì)算負(fù)荷模型由RL串、并聯(lián)組成；系統(tǒng)功率因數(shù)為0.8。具體參數(shù)如圖5所示，為方便計(jì)算，供電線路均采用三相三線制，所有的微電源和負(fù)荷為三相電源和負(fù)荷，仿真主要設(shè)計(jì)三相對(duì)稱接地短路和單相接地短路，利用電力系統(tǒng)分析軟件EMTP/ATP進(jìn)行仿真計(jì)算。

4 仿真分析

當(dāng)公共耦合母線發(fā)生故障，與其對(duì)應(yīng)的PCC(公共耦合點(diǎn))并網(wǎng)斷路器檢測到系統(tǒng)電壓降落，并有較大的故障電流持續(xù)從微電網(wǎng)涌入大電網(wǎng)，此時(shí)斷路器能夠及時(shí)動(dòng)作。微電網(wǎng)內(nèi)部各微電源在合理的控制策略下，維持微電網(wǎng)電壓和頻率穩(wěn)定，保證對(duì)微電網(wǎng)內(nèi)部負(fù)荷持續(xù)供電。

4.1 發(fā)生三相接地短路故障

假設(shè)圖5中公共耦合母線t=0.8 s時(shí)發(fā)生三相接地短路故障，故障持續(xù)0.4 s，t=1.2 s時(shí)故障解除。微電網(wǎng)暫態(tài)特性分析如下。

微電網(wǎng)各微電源出口端電壓值如圖6所示，0.8 s前，系統(tǒng)正常運(yùn)行，端電壓為220 V；0.8 s時(shí)系統(tǒng)發(fā)生三相接地短路故障，端電壓下降為正常值的20%左右；1.2 s后故障解除，端電壓恢復(fù)到短路前正常水平。由于風(fēng)機(jī)、太陽能電池等效短路阻抗不同，同時(shí)由于線路阻抗的存在，故障時(shí)其端電壓并不相同。

如圖7所示，0.8 s前系統(tǒng)電流處于正常水平，0.8 s系統(tǒng)發(fā)生三相接地短路故障。短路初期，圖7(a)中太陽能電池出口短路沖擊電流高達(dá)68 A，5個(gè)周波之后，由于直流分量迅速衰減到零，短路電流穩(wěn)定值為38.3 A左右，為正常值的6～7倍，短路沖擊系數(shù)為1.78；1.2 s后故障切除，系統(tǒng)恢復(fù)正常。圖7(b)中風(fēng)機(jī)出口短路沖擊電流高達(dá)120 A，3個(gè)周波之后，由于直流分量迅速衰減到零，短路電流穩(wěn)定值為67.54 A左右，為正常值的4～5倍，短路沖擊系數(shù)為1.75；1.2 s后故障切除，系統(tǒng)恢復(fù)正常。由于風(fēng)機(jī)、太陽能電池內(nèi)阻抗值不同，短路電流直流分量衰減速度存在一定的差距，阻抗L/R比值越低，直流分量衰減越快，短路電流達(dá)到暫態(tài)穩(wěn)定的速度也就越快。

4.2 發(fā)生單相接地短路故障

假設(shè)圖5中公共耦合母線t=0.8 s時(shí)發(fā)生A相接地短路故障，系統(tǒng)檢測到故障發(fā)生，斷路器斷開微電網(wǎng)與大電網(wǎng)之間連接，故障持續(xù)0.4 s，t=1.2 s后故障解除，微電網(wǎng)重新并網(wǎng)運(yùn)行。以風(fēng)機(jī)為例進(jìn)行系統(tǒng)故障暫態(tài)特性分析。

如圖8所示，0.8 s之前系統(tǒng)正常運(yùn)行，風(fēng)機(jī)各相端電壓均為220 V；0.8 s時(shí)系統(tǒng)發(fā)生A相接地短路故障，風(fēng)機(jī)A相端電壓下降為正常值的10%左右，而B、C兩相則上升為線電壓；1.2 s之后故障切除，A、B、C各相端電壓恢復(fù)到短路前正常水平。

如圖9所示，0.8 s之前，系統(tǒng)電流處于正常水平；0.8 s系統(tǒng)發(fā)生A相接地短路故障，與B、C相相比，A相短路電流幅值增加較大，為正常值的4倍左右；1.2 s之后系統(tǒng)恢復(fù)正常。

5 算例驗(yàn)證

為驗(yàn)證在EMTP/ATP平臺(tái)上搭建的模型及基于IEC60909標(biāo)準(zhǔn)的短路計(jì)算有效性，本部分采用手算的方法來驗(yàn)證，包括相關(guān)計(jì)算模型參數(shù)和最大化的三相接地短路故障結(jié)果，系統(tǒng)接線如圖5所示，建模過程中對(duì)系統(tǒng)進(jìn)行如下處理。

a.將太陽能發(fā)電單元、風(fēng)力發(fā)電單元以等值電壓源加內(nèi)阻抗串聯(lián)代替，其他參數(shù)與圖5相同。

b.忽略線路正序和零序的對(duì)地電容，采用等值線路參數(shù)。

c.負(fù)荷采用定負(fù)荷模型，以恒定阻抗代替。

表1 風(fēng)力發(fā)電單元三相短路計(jì)算結(jié)果對(duì)比

表2 太陽能發(fā)電單元三相短路計(jì)算結(jié)果對(duì)比

由于在數(shù)學(xué)模型和計(jì)算方法上的差異，EMTP/ATP和手算模型的計(jì)算結(jié)果不完全一致，但差別不大。從表1和表2的計(jì)算結(jié)果來看，短路電流交流分量初始值誤差較大，為13%左右，沖擊電流誤差值較小，基本可以驗(yàn)證模型的有效性。

引起計(jì)算結(jié)果誤差的原因有以下3個(gè)方面。

a.模型元素自身誤差，例如風(fēng)力發(fā)電單元、太陽能發(fā)電單元。

b.應(yīng)用典型參數(shù)帶來的誤差，例如輸電線路采用的典型參數(shù)與實(shí)際參數(shù)有一定偏差。

c.運(yùn)行方式選取帶來的誤差，IEC60909的結(jié)果實(shí)際是對(duì)最大短路電流的估算，并沒有提供系統(tǒng)運(yùn)行方式的詳細(xì)數(shù)據(jù)。

6 結(jié)束語

本文在電磁暫態(tài)仿真軟件EMTP/ATP環(huán)境中搭建風(fēng)/光互補(bǔ)發(fā)電的微電網(wǎng)模型，并基于IEC60909標(biāo)準(zhǔn)對(duì)微電網(wǎng)暫態(tài)短路特性進(jìn)行仿真，詳細(xì)分析了微電網(wǎng)公共耦合母線處三相接地短路和單相接地短路故障。研究表明在故障情況下，微電源出口電壓會(huì)降為正常值的20%左右，短路電流水平為正常值的6～10倍；故障解除后微電網(wǎng)能夠快速恢復(fù)穩(wěn)定運(yùn)行，滿足用戶的電能質(zhì)量要求；最后通過與人工計(jì)算的短路電流結(jié)果進(jìn)行對(duì)比，結(jié)果證明了模型的有效性。