吳 彤，陳洪濤，趙建明，李國(guó)慶

(1.東北電力大學(xué)電氣工程學(xué)院，吉林吉林132012;2.國(guó)網(wǎng)吉林省電力有限公司松原供電公司，吉林松原138000)

隨著分布式風(fēng)電電源和儲(chǔ)能裝置在電網(wǎng)中滲透率的增加，使得配電網(wǎng)從傳統(tǒng)的輻射狀網(wǎng)絡(luò)轉(zhuǎn)變?yōu)槎嚯娫垂╇娔Ｊ?，進(jìn)而改變了配電系統(tǒng)的結(jié)構(gòu)和運(yùn)行模式，這無(wú)疑會(huì)對(duì)配電網(wǎng)的電壓分布、潮流走向等方面產(chǎn)生影響［1－3］。由于傳統(tǒng)潮流計(jì)算算法無(wú)法對(duì)分布式電源進(jìn)行處理，國(guó)內(nèi)外研究學(xué)者開(kāi)始著手研究更高效、實(shí)用的方法來(lái)滿足復(fù)雜配電網(wǎng)快速計(jì)算和分析的需求。

文獻(xiàn)［4］提出了一種基于回路分析法的簡(jiǎn)化潮流計(jì)算方法，將回路電壓方程轉(zhuǎn)化為線性代數(shù)方程進(jìn)行求解;文獻(xiàn)［5］對(duì)網(wǎng)絡(luò)中PV節(jié)點(diǎn)進(jìn)行優(yōu)化，進(jìn)而提高潮流計(jì)算效率;文獻(xiàn)［6］采用改進(jìn)的前推回代法進(jìn)行含分布式電源的配電網(wǎng)潮流計(jì)算。這些計(jì)算方法都是對(duì)傳統(tǒng)的潮流算法進(jìn)行改進(jìn)，使其能夠更好的滿足電網(wǎng)發(fā)展的需求，但在計(jì)算過(guò)程中都忽略了儲(chǔ)能裝置對(duì)潮流計(jì)算結(jié)果的影響。

為此，本文提出了一種簡(jiǎn)化實(shí)用的潮流計(jì)算方法，在綜合考慮風(fēng)電隨機(jī)性和儲(chǔ)能裝置雙向潮流特性的基礎(chǔ)上，對(duì)傳統(tǒng)潮流計(jì)算方法加以改進(jìn)，使潮流計(jì)算的迭代次數(shù)更少，收斂效果更好。

1 風(fēng)力發(fā)電機(jī)隨機(jī)出力模型

1.1 風(fēng)力發(fā)電機(jī)模型

目前風(fēng)速模型普遍采用Weibull分布，其概率密度函數(shù)為

式中:v為風(fēng)速;α為形狀參數(shù);β為尺度參數(shù)。

風(fēng)力發(fā)電機(jī)的出力與風(fēng)速的關(guān)系基本服從圖1所示的出力曲線。圖中:vci為切入風(fēng)速;vr為額定風(fēng)速;vco為切出風(fēng)速;Pr為風(fēng)機(jī)額定輸出功率。

圖1 風(fēng)力發(fā)電機(jī)出力與風(fēng)速關(guān)系曲線

如采用異步發(fā)電機(jī)，則風(fēng)力發(fā)電機(jī)的簡(jiǎn)化等值電路如圖2所示。由圖中關(guān)系可推導(dǎo)出風(fēng)力發(fā)電機(jī)的機(jī)端電壓表達(dá)式和無(wú)功功率表達(dá)式。

式中:R2為轉(zhuǎn)子電阻;X為定、轉(zhuǎn)子電抗之和;Xm為勵(lì)磁電抗;Xc為機(jī)端補(bǔ)償電容電抗;s為轉(zhuǎn)差率。

在風(fēng)速一定的條件下，風(fēng)力發(fā)電機(jī)的有功出力可視為固定值，則風(fēng)力發(fā)電機(jī)的無(wú)功輸出僅與轉(zhuǎn)差率和機(jī)端電壓有關(guān)。根據(jù)公式(2)和公式(3)消去轉(zhuǎn)差率s即可得出無(wú)功與電壓的關(guān)系Q=f(U)為:

式中，各變量含義與公式(2)－(3)中所給出的變量含義相同。

圖2 異步發(fā)電機(jī)簡(jiǎn)化等值電路

1.2 風(fēng)機(jī)出力概率分布情況

將風(fēng)速概率密度函數(shù)帶入風(fēng)機(jī)出力表達(dá)式即可得到風(fēng)機(jī)出力Pw的概率表達(dá)式［7］。

(1)Pw=0時(shí)，0 ＜ v≤vci∪v＞vco。

則有零出力情況下風(fēng)機(jī)出力概率為:

(2)0＜Pw＜Pr時(shí)vci＜v＜vr，第k段概率為

概率密度函數(shù)可由式(7)得到

欠出力離散段風(fēng)機(jī)出力Pk采用該段有功出力的中間值，即期望值。

(3)Pw=Pr時(shí)，vr＜ v≤ vco。

在式(4)－(8)中，v為風(fēng)速;α為形狀參數(shù);β為尺度參數(shù);vci為切入風(fēng)速;vr為額定風(fēng)速;vco為切出風(fēng)速;Pr為風(fēng)機(jī)額定輸出功率。

2 潮流計(jì)算中風(fēng)機(jī)和儲(chǔ)能節(jié)點(diǎn)的處理

2.1 潮流計(jì)算中風(fēng)機(jī)節(jié)點(diǎn)的處理

由于異步風(fēng)力發(fā)電機(jī)組并網(wǎng)點(diǎn)無(wú)功需同時(shí)滿足公式(4)和潮流方程，如果只將其處理成一個(gè)簡(jiǎn)單的PQ節(jié)點(diǎn)進(jìn)行潮流計(jì)算會(huì)與實(shí)際情況相差較大。為此，本文采用如下處理方法將風(fēng)機(jī)出力看成一個(gè)輸出無(wú)功隨電壓變化的PQ(V)節(jié)點(diǎn)，方法如下:

對(duì)于雙饋式和直驅(qū)式風(fēng)機(jī)，由于其機(jī)組電壓相對(duì)穩(wěn)定且對(duì)系統(tǒng)無(wú)功需求較小，故在潮流計(jì)算中可將其視為PV節(jié)點(diǎn)。

2.2 潮流計(jì)算中儲(chǔ)能節(jié)點(diǎn)的處理

對(duì)于儲(chǔ)能裝置在電網(wǎng)中既可作為負(fù)載儲(chǔ)存電能，又可以作為分布式電源向電網(wǎng)供電。當(dāng)工作在整流狀態(tài)時(shí)，儲(chǔ)能裝置處于充電狀態(tài)，能量從電網(wǎng)側(cè)流向直流側(cè);當(dāng)工作在逆變狀態(tài)時(shí)，儲(chǔ)能裝置處于放電狀態(tài)，將直流側(cè)的能量回饋給電網(wǎng)。當(dāng)儲(chǔ)能裝置采用電流型逆變器時(shí)，做PI節(jié)點(diǎn)處理;采用電壓控制逆變器時(shí)，做PV節(jié)點(diǎn)處理，本文所采用的儲(chǔ)能裝置選用電壓控制逆變器。［8－10］

另外，由于儲(chǔ)能裝置容量受限，不可能一直處于充電或者放電狀態(tài)，故需要相應(yīng)的控制策略來(lái)實(shí)現(xiàn)其自動(dòng)充放電功能。

2.3 儲(chǔ)能裝置控制策略

為確保電網(wǎng)能夠穩(wěn)定、可靠運(yùn)行，儲(chǔ)能裝置必須依照一定的控制策略進(jìn)行充放電。本文儲(chǔ)能裝置控制策略的數(shù)學(xué)模型如下:

式中:P(t－1)為風(fēng)力發(fā)電機(jī)在t－1時(shí)刻平抑后的輸出的功率;Pmax(t)、P(t)分別為t時(shí)刻輸出的最大輸出功率和平抑后的值;Pess(t)為儲(chǔ)能裝置的充放電功率;S、λ分別為裝機(jī)容量和允許的輸出功率的最大變化率。若Pess(t)＞0，則儲(chǔ)能裝置放電;若Pess(t)＜0，則儲(chǔ)能裝置充電。

也就是說(shuō)，當(dāng)風(fēng)電輸出的有功功率波動(dòng)值大于其允許的最大波動(dòng)值時(shí)，系統(tǒng)中會(huì)有能量剩余，此時(shí)將多余的能量送至儲(chǔ)能裝置進(jìn)行充電，當(dāng)儲(chǔ)能裝置充滿后，停止充電，剩余的有功功率按照能源不浪費(fèi)原則輸出;當(dāng)風(fēng)電輸出的有功功率波動(dòng)值小于其允許的最小波動(dòng)值時(shí)，系統(tǒng)中會(huì)產(chǎn)生功率缺額，此時(shí)儲(chǔ)能裝置進(jìn)行放電，當(dāng)儲(chǔ)能裝置電量放盡后，下一時(shí)刻開(kāi)始不再放電。

3 潮流算法的實(shí)現(xiàn)

(1)讀取原始數(shù)據(jù)，設(shè)定電壓初始值。

(2)根據(jù)風(fēng)力發(fā)電出力的概率分布和儲(chǔ)能裝置的控制策略確定系統(tǒng)最具代表性時(shí)刻，以此時(shí)刻下風(fēng)機(jī)出力和儲(chǔ)能裝置的帶電狀態(tài)為初始狀態(tài)。

(3)判定系統(tǒng)中節(jié)點(diǎn)類型，將風(fēng)力發(fā)電機(jī)節(jié)點(diǎn)轉(zhuǎn)化為PQ節(jié)點(diǎn)，將儲(chǔ)能節(jié)點(diǎn)轉(zhuǎn)化為PV節(jié)點(diǎn)，判斷儲(chǔ)能節(jié)點(diǎn)是否存在無(wú)功越限的情況。若出現(xiàn)無(wú)功越限，將儲(chǔ)能節(jié)點(diǎn)轉(zhuǎn)化為PQ節(jié)點(diǎn)參與潮流計(jì)算;若無(wú)功沒(méi)有越限，將儲(chǔ)能節(jié)點(diǎn)以PV節(jié)點(diǎn)參與潮流計(jì)算。

(4)由線路末端電流和功率前推計(jì)算各節(jié)點(diǎn)電流、功率。

(5)以平衡節(jié)點(diǎn)電壓回代計(jì)算線路末端節(jié)點(diǎn)電壓。

(6)計(jì)算系統(tǒng)中PV節(jié)點(diǎn)相鄰兩次迭代電壓的差值，若滿足精度，則結(jié)束潮流計(jì)算。

4 仿真驗(yàn)證

4.1 系統(tǒng)基本資料

本文在Matlab環(huán)境下進(jìn)行程序編譯來(lái)實(shí)現(xiàn)文中所提出的潮流計(jì)算方法，并采用PSAT(Power System Analysis Toolbox)搭建仿真模型，通過(guò)潮流計(jì)算結(jié)果對(duì)比來(lái)驗(yàn)證本文所提算法的有效性。

以IEEE33節(jié)點(diǎn)系統(tǒng)為例，系統(tǒng)電壓基準(zhǔn)值為12.66 kV，視在功率基準(zhǔn)值為10 MVA。系統(tǒng)中包含4臺(tái)額定功率為300 kW的風(fēng)機(jī)和一額定容量為480 kW·h的儲(chǔ)能蓄電池。其中1節(jié)點(diǎn)為平衡節(jié)點(diǎn)，風(fēng)機(jī)接入節(jié)點(diǎn)分別為7節(jié)點(diǎn)、14節(jié)點(diǎn)、25節(jié)點(diǎn)和32節(jié)點(diǎn)，儲(chǔ)能接入節(jié)點(diǎn)為11節(jié)點(diǎn)。仿真模型如圖3所示:

圖3 IEEE33節(jié)點(diǎn)測(cè)試系統(tǒng)仿真模型

其中，風(fēng)機(jī)的額定電壓為0.69 kV，經(jīng)0.69/12.66變壓器與配電網(wǎng)相連，儲(chǔ)能裝置初始化處理為恒定PV節(jié)點(diǎn)，但同時(shí)受控制策略限制其進(jìn)行充放電。

4.2 仿真計(jì)算結(jié)果

采用Matlab編程進(jìn)行計(jì)算，將計(jì)算所得結(jié)果與PSAT潮流計(jì)算結(jié)果進(jìn)行對(duì)比，結(jié)果如表1所示。

表1 Matlab編程與PSAT仿真結(jié)果對(duì)比

由表1中數(shù)據(jù)對(duì)比可以看出，經(jīng)Matlab編程求得的各節(jié)點(diǎn)電壓與PSAT仿真模型求得的節(jié)點(diǎn)電壓基本相同，且誤差在允許范圍內(nèi)。說(shuō)明本文提出的潮流計(jì)算方法可以有效的對(duì)含分布式風(fēng)電和儲(chǔ)能裝置的配電網(wǎng)進(jìn)行潮流計(jì)算。

5 結(jié) 論

本文根據(jù)風(fēng)機(jī)出力特性和儲(chǔ)能裝置的雙向潮流特性，提出了一種簡(jiǎn)化的潮流計(jì)算方法，用于解決傳統(tǒng)潮流計(jì)算方法不能處理分布式風(fēng)電和儲(chǔ)能裝置的問(wèn)題。采用Matlab編程來(lái)實(shí)現(xiàn)本文所提出的潮流計(jì)算方法，并以PSAT工具箱計(jì)算的結(jié)果進(jìn)行對(duì)比，結(jié)果表明，本文所提出的潮流計(jì)算方法可以有效對(duì)含分布式風(fēng)電和儲(chǔ)能裝置的配電網(wǎng)進(jìn)行潮流計(jì)算。

［1］孔瑾.智能配電網(wǎng)的電網(wǎng)分區(qū)及孤島辨識(shí)研究［D］.山東:山東大學(xué)，2013.

［2］楊華，吳政球，涂有慶，等.分布式發(fā)電及其并網(wǎng)技術(shù)綜述［J］.電網(wǎng)技術(shù)，2008，32(15):71－76.

［3］莫穎濤.影響我國(guó)分布式發(fā)電發(fā)展的關(guān)鍵因素［J］.電網(wǎng)技術(shù)，200832(S1):176－178.

［4］牛煥娜，井天軍，李漢成，等.基于回路分析的含分布式電源配電網(wǎng)簡(jiǎn)化潮流計(jì)算［J］.電網(wǎng)技術(shù)，2013，37(4):1032－1038.

［5］楊超，沈聰，李睿.含分布式電源配電網(wǎng)的前推回代潮流算法中PV節(jié)點(diǎn)處理方法［J］.電網(wǎng)技術(shù)，2012，36(9):238－243.

［6］王建勛，呂群芳，劉會(huì)金，等.含分布式電源的配電網(wǎng)潮流快速直接算法［J］.電力自動(dòng)化設(shè)備，2011，31(2):17－21.

［7］馬瑞，蔣斌.含風(fēng)電場(chǎng)并考慮風(fēng)速隨機(jī)性的電力系統(tǒng)電壓控制區(qū)域劃分［J］.長(zhǎng)沙理工大學(xué)學(xué)報(bào)，2008，5(3):71－76.

［8］Anju G Pillai，Thomas P C，Sreerenjini K，et al.Transient Stability Analysis of Wind Integrated Power Systems with Storage Using Central Area Controller［J］.2013 Annual International Conference on Emerging Research Areas and 2013 International conference on Microelectronics.2013:1－5.

［9］S Pérez－Londno，L Rodriguez－Garcia，Y Ulianov López.Effects of Doubly Fed Wind Generators on Voltage Stability of Power Systems［J］.2012 Sixth IEEE/PES Transmission and Distribution:Latin America Conference and Exposition(T＆D－LA).2012:1－6.

［10］陳厚合，呂志遠(yuǎn)，朱景明.基于混合連續(xù)蟻群算法的最優(yōu)潮流研究［J］.東北電力大學(xué)學(xué)報(bào)，2009，29(6):54－59.