徐笑， 李志鏗， 楊海森， 廖威， 郭俊宏

(1.深圳供電局有限公司，廣東 深圳 518100；2.中國能源建設(shè)集團廣東省電力設(shè)計研究院有限公司，廣東 廣州 510530；3.廣東電網(wǎng)有限責(zé)任公司電力調(diào)度控制中心，廣東 廣州 510600)

0 引 言

近年來分布式電源(Distributed Generation,DG)在配電網(wǎng)側(cè)得到廣泛應(yīng)用，不僅滿足用戶側(cè)需求又保證電網(wǎng)清潔，但也對傳統(tǒng)配電網(wǎng)的運行和規(guī)劃帶來嚴峻的考驗[1]。主動配電網(wǎng)(Active Distribution Network, ADN)能夠?qū)崿F(xiàn)對配電網(wǎng)內(nèi)的分布式能源主動控制和主動管理，改善潮流分布[2]。合理的ADN規(guī)劃不僅可以節(jié)省投資費用，提高供電質(zhì)量，還能消納大規(guī)模間歇式可再生能源[3]。

經(jīng)濟性和可靠性是考量主動配電網(wǎng)規(guī)劃的兩大重要指標[4]，但通常二者又相對矛盾。為了提高配網(wǎng)規(guī)劃的可靠性，一方面不可避免要增強電網(wǎng)的電氣聯(lián)系，增加建設(shè)成本；另一方面，減少投資建設(shè)費用又可能降低供電質(zhì)量和可靠性。所以如何權(quán)衡與協(xié)調(diào)二者的關(guān)系是主動配電網(wǎng)規(guī)劃研究的重點。

主動配電網(wǎng)可靠性規(guī)劃主要考慮網(wǎng)架敷設(shè)、變電站規(guī)劃、DG選址定容和無功設(shè)備配置等[5]，屬于一個復(fù)雜的混合整數(shù)非線性規(guī)劃(MINLP)問題。目前常見的求解算法包括：啟發(fā)式規(guī)劃算法、數(shù)學(xué)規(guī)劃算法[6-7]和人工智能優(yōu)化算法[8-10]。

廣義的Benders分解算法[11-12]可以將復(fù)雜的耦合約束和變量分開處理，形成相對容易解決的主問題和子問題，能夠較好地處理MINLP這類問題。因此，本文基于廣義Benders分解算法，將配網(wǎng)網(wǎng)架結(jié)構(gòu)和DG選址配置作為整數(shù)決策變量，建立以綜合費用最小為目標的規(guī)劃主問題；再對規(guī)劃后的主網(wǎng)架進行N-1安全性校驗，建立以年缺供電費最小為目標的可靠性子問題。通過迭代主問題和子問題以獲得最優(yōu)網(wǎng)架和DG規(guī)劃方案。

1 AND可靠性規(guī)劃數(shù)學(xué)模型

將AND可靠性規(guī)劃分成主問題和子問題:主問題為ADN規(guī)劃成本優(yōu)化問題，以年均綜合費用最小為目標，包括線路網(wǎng)架和DG的年等值建設(shè)費用，設(shè)備維護費用以及對尚未通電負荷節(jié)點的停電補償費用；子問題為在主問題所得的配網(wǎng)規(guī)劃下進行N-1安全性校驗，以年均缺供電費最小為目標的可靠性問題。

1.1 AND可靠性規(guī)劃主問題

主問題的目標函數(shù)表達式如下：

minT=Iv+Mv+Sv+Lv

(1)

(2)

(3)

(4)

式中:T為總投資費用;Iv為網(wǎng)架和DG的投資費用;Mv為設(shè)備維護費用;Sv為停電補償費用;Lv為ADN規(guī)劃的可靠性費用，將在1.2節(jié)中給出;r為折現(xiàn)率，取為8%；Nt表示配電網(wǎng)規(guī)劃年限，取為3年;t為規(guī)劃年的編號；Nk表示待建線路的數(shù)目;k為線路編號；Ng表示待并網(wǎng)DG的數(shù)目;g為DG編號;c1為線路建設(shè)費用；c2為DG建設(shè)費用；c3為設(shè)備的維護費用；c4為停電補償單價費用;i為節(jié)點編號;αij,t、γi,t、δi,t為01變量,αij,t=1表示第t個規(guī)劃年間線路ij連通；γi,t=1表示第t個規(guī)劃年間在節(jié)點i處接入DG；δi,t=1表示節(jié)點i負荷在第t個規(guī)劃年投運；αij,t-1、γi,t-1為t-1年相應(yīng)情況。

約束條件包括：

(1)配電網(wǎng)的輻射狀約束

設(shè)置變電站所在節(jié)點為根節(jié)點，網(wǎng)架規(guī)劃從根節(jié)點逐漸延伸而成。當(dāng)βij=1，βji=0時，表示節(jié)點i是節(jié)點j的父節(jié)點；當(dāng)二者均為0時，則代表ij兩節(jié)點沒有連通。

(5)

式中:Nr為變電站節(jié)點數(shù)目;Ωr為根節(jié)點集合。

(2)配電網(wǎng)設(shè)備建設(shè)約束

線路、DG等設(shè)備的投建具有不可逆轉(zhuǎn)性。

(6)

(3)潮流約束

基于文獻[13]的DistFlow潮流計算公式，同時為了提高算法的收斂性，將其簡化如下：

(7)

式中:Vi、Vj分別為節(jié)點i的電壓幅值；Pij，Qij為線路ij流通的有功、無功功率；Rij和Xij分別為線路ij的電阻和電抗。

結(jié)合變量αij,t、γi,t、δi,t，可得到最終潮流模型：

(8)

(4)電壓幅值與線路容量約束

(9)

式中:Vi,max，Vi,min為節(jié)點i電壓幅值的上下限;Iij,max，Iij,min為線路ij電流幅值上下限。

1.2 AND可靠性規(guī)劃子問題

在得到主動配電網(wǎng)規(guī)劃方案后進行N-1安全性校驗以驗證規(guī)劃方案的可靠性。子問題為可靠性費用，以年缺供電費最小為目標：

(10)

約束條件為：

(12)

2 求解算法

在廣義Benders分解算法的框架下，將目標函數(shù)改寫為以下緊湊形式：

minT=C(αij,t,γi,t,δi,t)+ω(n)

(13)

s.t.αij,t(n),γi,t(n),δi,t(n)∈f(αij,t,γi,t,δi,t)

(14)

(15)

ω(n)≥ω-

(16)

式中:C(αij,t，γi,t，δi,t)包括網(wǎng)架和DG的建設(shè)費用，設(shè)備維護費用以及停電補償費用;ω為一個連續(xù)變量，用以表示子問題的優(yōu)化結(jié)果;f(αij,t，γi,t，δi,t)為由主問題的約束條件確定的變量αij,t，γi,t，δi,t取值范圍;n為迭代次數(shù);D(σ(n)i,s,t)為第n次迭代后的子問題函數(shù)值;μ(n)ij,s,t、η(n)i,t、ψ(n)i,t為第n次迭代子問題得到的對偶變量;α(n)ij,s,t、γ(n)i,t、δ(n)i,t分別為第n次迭代得到變量αij,s,t、γi,t、δi,t的值;ω-為ω初始化的下限值。式(16)為Benders割集，即該算法的核心。

子問題則是在主問題所求解的網(wǎng)架規(guī)劃與DG選址方案的基礎(chǔ)上，求解最小的年缺供電費：

minD(σi,s.t)

(17)

s.t.σi,s,t∈g(σi,s,t)

(18)

(19)

(20)

(21)

式中:g(σi,s,t)為子問題的約束條件，由1.2節(jié)中約束條件確定的變量σi,s,t取值范圍。

具體步驟如下：

(1)初始化上下界，UB(0)=+∞，LB(0)=0；收斂精度ε取作10-6。

(2)求解主問題，得到一組網(wǎng)架規(guī)劃和DG選址配置的可行解以及迭代過程中的下界LB(n)，LB(n)=C(α(n)ij,,t、γ(n)i,t、δ(n)i,t)+ω(n),將這組可行解傳遞給子問題作為ADN的決策規(guī)劃。

(3)求解子問題，得到在這組可行解下的最低可靠性成本D[σ(n)i,s,t]。更新對偶變量μ(n)ij,s,t、η(n)i,t、ψ(n)i,t信息，產(chǎn)生新的Benders割集添加到主問題的約束條件中。計算本次的上界UB(n)=D[σ(n)i,s,t]+C[α(n)ij,t、γ(n)i,t、δ(n)i,t]。

(4)判斷UB(n)-LB(n)是否小于收斂精度ε，如果是,則將當(dāng)前的可行解作為最優(yōu)解輸出；反之則返回步驟(3)，令n=n+1。

(5)算法迭代結(jié)束。

3 算例分析

原始網(wǎng)架拓撲圖如圖1所示，共有54個節(jié)點，61條可選規(guī)劃線路。規(guī)劃年限為3 a，電壓等級為10 kV，容量基準值為1 MVA，λs均為0.01。規(guī)劃期內(nèi)并網(wǎng)DG的數(shù)量為6臺，每臺DG的功率為PG=1.5 MW，QG=0.5 Mvar。DG投入成本c2=1萬元/臺，設(shè)備維護成本c3=0.5萬元/臺，停電補償費用單價c4=4萬元/MW，缺供電費用單價c5為1萬元/MW。

對上述算例進行求解得到最優(yōu)的網(wǎng)架規(guī)劃和DG選址配置方案如圖2所示，其中圖2(a)、圖2(b)、圖2(c)分別為第1～第3年的規(guī)劃結(jié)果圖。

圖2 所提模型規(guī)劃結(jié)果圖

對比算例2，利用遺傳算法求解，得到的網(wǎng)架規(guī)劃和DG選址配置方案如圖3所示。

為了驗證所提的網(wǎng)架、DG協(xié)調(diào)規(guī)劃與先確定DG后規(guī)劃網(wǎng)架的縱向規(guī)劃兩種方案的優(yōu)劣，算例3選取規(guī)劃網(wǎng)架中負荷較大的6個節(jié)點作為已選定的DG安裝位置，再對網(wǎng)架進行規(guī)劃。3個算例的各項費用計算結(jié)果列于表1，求解時間列于表2。

圖3 遺傳算法規(guī)劃結(jié)果圖

對比可知：算例1和算例2的DG最終安裝位置相同，但是每臺DG的具體安裝年限卻不盡相同，這也會導(dǎo)致網(wǎng)架線路的規(guī)劃方案不同；含有DG的線路往往較長，連接更多的負荷節(jié)點，這是因為DG不僅能就近供應(yīng)負荷，而且當(dāng)線路發(fā)生故障時能形成孤島，恢復(fù)部分負載，降低可靠性費用；算例1的可靠性費用為3.81萬元，總費用為401.86萬元，算例2的可靠性費用為4.12萬元，總費用為408.73萬元。由此可見，算例1的ADN規(guī)劃不僅經(jīng)濟性優(yōu)于算例2，且具有更高的可靠性。算例1較算例2盡管在計算時間上不占優(yōu)勢，但能夠得到更理想的優(yōu)化結(jié)果。

表1 算例計算結(jié)果對比

表2 算例求解時間對比

算例3綜合費用較算例1多41.04萬元，這表明先確定DG位置再進行網(wǎng)架優(yōu)化規(guī)劃并不能獲得最佳的ADN規(guī)劃方案，只有將DG選址規(guī)劃和網(wǎng)架規(guī)劃協(xié)同考慮，才能更好地發(fā)揮DG在減緩ADN規(guī)劃建設(shè)費用方面的作用，同時能夠充分利用DG就地供給電能的優(yōu)勢，減少因尚未建設(shè)而造成的停電補償費用。算例3的求解時間較算例1少了0.503 s，這是因為算例3不需要對DG選址進行最優(yōu)決策，減少了ADN規(guī)劃的變量數(shù)量,自然減少了求解時間。

4 結(jié)束語

利用廣義Benders分解算法將主動配電網(wǎng)可靠性規(guī)劃分成ADN規(guī)劃成本決策優(yōu)化主問題和可靠性效益優(yōu)化子問題，通過算例分析驗證了所提方法的有效性和可行性：相較于遺傳算法，所提的模型求解方法能夠更好地收斂到主動配電網(wǎng)的最優(yōu)規(guī)劃方案，在保證求解準確性的基礎(chǔ)上仍有較快的求解速度；DG的接入位置影響配網(wǎng)的潮流分布，在主動配電網(wǎng)規(guī)劃中協(xié)同考慮DG選址規(guī)劃和網(wǎng)架規(guī)劃，才能充分體現(xiàn)DG并網(wǎng)對ADN規(guī)劃的影響，從而獲得整體最優(yōu)的主動配電網(wǎng)規(guī)劃方案。