袁天夢(mèng)，侯鑫垚，廉杰

(國(guó)網(wǎng)冀北電力有限公司唐山供電公司，唐山 063000)

隨著能源全球化進(jìn)程的不斷推進(jìn)，高效、清潔、多樣化的現(xiàn)代能源系統(tǒng)逐漸成為中國(guó)能源發(fā)展轉(zhuǎn)型的主要途徑。當(dāng)今世界多數(shù)國(guó)家均使用大功率的單一供電系統(tǒng)為電力用戶(hù)提供能源，這種集中式供電方式為95%左右的用戶(hù)提供服務(wù)。隨著電力用戶(hù)的逐漸增加，如何保證能源供應(yīng)與安全，提高能源利用率，降低污染物排放成為電力系統(tǒng)發(fā)展過(guò)程中急需解決的問(wèn)題[1-2]。針對(duì)上述問(wèn)題，能源互聯(lián)網(wǎng)應(yīng)運(yùn)而生，此網(wǎng)絡(luò)中不僅包含多種能源單元，更能滿(mǎn)足用戶(hù)的多元化應(yīng)用需求，為集中式供電方式提供強(qiáng)有力的能源保證。分布式能源系統(tǒng)是能源互聯(lián)網(wǎng)中的重要組成部分，可為區(qū)域性配電網(wǎng)的能源提供供給[3]。分布式能源系統(tǒng)的發(fā)展不僅提高了清潔能源的利用率，更為中國(guó)電力事業(yè)與環(huán)保事業(yè)的發(fā)展提供了幫助。

國(guó)內(nèi)外一些學(xué)者對(duì)此展開(kāi)了研究，王永利等[4]提出了一種考慮需求響應(yīng)的分布式能源系統(tǒng)多能協(xié)同與儲(chǔ)能規(guī)劃方法，通過(guò)建立由風(fēng)機(jī)、光伏、熱電聯(lián)產(chǎn)和電鍋爐、儲(chǔ)能系統(tǒng)組成的綜合能源系統(tǒng)多能協(xié)同模型，通過(guò)規(guī)劃各設(shè)備容量情況，滿(mǎn)足電熱負(fù)荷需求，并采用遺傳算法對(duì)該規(guī)劃方案進(jìn)行優(yōu)化。易靈芝等[5]提出了一種含微電網(wǎng)群的主動(dòng)配電網(wǎng)雙層聯(lián)合優(yōu)化調(diào)度方法，上層模型主要是以配電網(wǎng)為研究對(duì)象，應(yīng)用遺傳算法來(lái)降低配電網(wǎng)線(xiàn)損，提高電能質(zhì)量；下層模型主要是以微電網(wǎng)為研究對(duì)象，應(yīng)用混合整數(shù)線(xiàn)性規(guī)劃來(lái)降低運(yùn)行成本。張學(xué)軍等[6]提出了一種含風(fēng)儲(chǔ)聯(lián)合單元的電力系統(tǒng)兩階段日前調(diào)度模型，針對(duì)電力市場(chǎng)環(huán)境下含風(fēng)儲(chǔ)聯(lián)合單元的電力系統(tǒng)日前調(diào)度問(wèn)題，建立了兩階段調(diào)度優(yōu)化模型。第一階段主要是以風(fēng)儲(chǔ)聯(lián)合單元的期望收益最大為目標(biāo)，對(duì)日前投標(biāo)量進(jìn)行優(yōu)化；第二階段主要是以系統(tǒng)綜合運(yùn)行成本最小為目標(biāo)，對(duì)日前機(jī)組組合進(jìn)行優(yōu)化。王子馳等[7]提出了一種基于多能互補(bǔ)分布式能源與綜合能源管理系統(tǒng)優(yōu)化調(diào)度方法，通過(guò)構(gòu)建某一區(qū)域的多能互補(bǔ)分布式能源系統(tǒng)，采用負(fù)荷預(yù)測(cè)的方法進(jìn)行系統(tǒng)的優(yōu)化，實(shí)現(xiàn)系統(tǒng)調(diào)度的最大經(jīng)濟(jì)效益。劉華等[8]提出了一種天然氣冷熱電聯(lián)供分布式能源系統(tǒng)優(yōu)化調(diào)度方法，基于改進(jìn)粒子群算法，提出天然氣冷熱電聯(lián)供系統(tǒng)中各設(shè)備的優(yōu)化調(diào)度模型，并對(duì)設(shè)備運(yùn)行進(jìn)行條件約束，確保運(yùn)行成本最低。以上方法都可以進(jìn)行分布式能源調(diào)度，并且可進(jìn)行多角度調(diào)度。但在主動(dòng)配電網(wǎng)下的分布式能源系統(tǒng)雙層兩階段調(diào)度中，優(yōu)化模型分層調(diào)度能力不高，使得能源利用率不高。

為解決以上方法在主動(dòng)配電網(wǎng)下的分布式能源系統(tǒng)雙層兩階段調(diào)度中的問(wèn)題，在考慮能源利用率和運(yùn)行成本的基礎(chǔ)上，現(xiàn)構(gòu)建一種主動(dòng)配電網(wǎng)下分布式能源系統(tǒng)雙層兩階段調(diào)度優(yōu)化模型。首先，構(gòu)建主動(dòng)配電網(wǎng)雙層能量管理下分布式能源系統(tǒng)日前與日內(nèi)兩階段調(diào)度模型，并以提高能源系統(tǒng)的利用率作為目標(biāo)函數(shù)，設(shè)定雙層兩階段調(diào)度優(yōu)化模型的約束條件，控制能源系統(tǒng)的調(diào)度過(guò)程。其次，使用帝國(guó)競(jìng)爭(zhēng)優(yōu)化算法對(duì)能源系統(tǒng)調(diào)度優(yōu)化模型進(jìn)行求解，獲取最佳調(diào)度方案。最后，構(gòu)建算例分析環(huán)節(jié)，將該模型和算法引入IEEE-33節(jié)點(diǎn)配電網(wǎng)中，進(jìn)行測(cè)試，驗(yàn)證所提模型的有效性。

1 調(diào)度優(yōu)化模型構(gòu)建與求解

1.1 主動(dòng)配電網(wǎng)下分布式能源系統(tǒng)雙層兩階段調(diào)度優(yōu)化模型構(gòu)建

主動(dòng)配電網(wǎng)下分布式能源系統(tǒng)調(diào)度優(yōu)化的目的是在確保配電網(wǎng)穩(wěn)定運(yùn)行的前提下，使可再生能源能夠得到最大化利用，系統(tǒng)的總體運(yùn)行成本最低。由于傳統(tǒng)方法下的優(yōu)化模型對(duì)于日前與日內(nèi)兩階段模型構(gòu)建結(jié)果較好，但模型分層效果不佳[9-10]。因此，將優(yōu)化模型構(gòu)建為上、下兩層，上層負(fù)責(zé)節(jié)點(diǎn)功率的調(diào)度，以網(wǎng)損最小為目標(biāo)，進(jìn)行全局優(yōu)化，下層以運(yùn)行總成本最低為目標(biāo)，局部?jī)?yōu)化[11-12]。在調(diào)度的過(guò)程中考慮日前與日內(nèi)兩個(gè)階段的能源系統(tǒng)作用與表現(xiàn)，優(yōu)化模型可表示為

(1)

式(1)中：b為下層調(diào)度變量；a為上層調(diào)度變量；f(a,b)為目標(biāo)函數(shù)。

根據(jù)式(1)可得到a、b為指定數(shù)據(jù)時(shí)，調(diào)度的最優(yōu)解為

(2)

已知上層目標(biāo)函數(shù)為滿(mǎn)足多種約束條件的能源系統(tǒng)運(yùn)行成本，在研究中，將分布式能源系統(tǒng)中的微型燃?xì)廨啓C(jī)運(yùn)行成本作為主要研究對(duì)象。通過(guò)文獻(xiàn)[13-14]分析可知，能源系統(tǒng)供能具有時(shí)段性，為保證此系統(tǒng)的正常運(yùn)轉(zhuǎn)，需在其中安裝相應(yīng)的負(fù)荷穩(wěn)定器，則此電源的運(yùn)行成本Fl可表示為

Fl=ui(Pi)2+viPi+si

(3)

式(3)中：ui、vi和si為燃料費(fèi)用系數(shù)；Pi為微型燃?xì)廨啓C(jī)組正常運(yùn)行功率。

除此成本外，風(fēng)電機(jī)組與光伏機(jī)組不計(jì)入燃料成本，其發(fā)電成本函數(shù)可使用Ji、Ki表示，則有

Ji=PjiRJ

(4)

Ki=PkiRK

(5)

式中：Pji、RJ分別為風(fēng)電機(jī)組的有效出力與維護(hù)費(fèi)用；Pki、RK分別為光伏機(jī)組的有效出力與維護(hù)費(fèi)用。將上述計(jì)算內(nèi)容融入當(dāng)前上層目標(biāo)函數(shù)中。

在下層調(diào)度模型的構(gòu)建過(guò)程中，需對(duì)系統(tǒng)中的各個(gè)電源分配情況進(jìn)行分析，以提高能源系統(tǒng)的利用率作為目標(biāo)函數(shù)。則此目標(biāo)函數(shù)[15]可表示為

minOz=αiOi+αjOj

(6)

式(6)中：Oi為調(diào)度過(guò)程中，能源系統(tǒng)的總網(wǎng)損函數(shù)；Oj為能源系統(tǒng)中各電壓節(jié)點(diǎn)偏移函數(shù)；αi與αj分別為上述兩系數(shù)在下層調(diào)度函數(shù)中的權(quán)重系數(shù)。

為保證上下層目標(biāo)函數(shù)具有科學(xué)性，將其約束條件設(shè)定如下。

(7)

式(7)中：Pimin為微型燃?xì)廨啓C(jī)組i的出力下限；Pjmax、Pkmax分別對(duì)應(yīng)能源系統(tǒng)中的風(fēng)電機(jī)組與光伏機(jī)組的出力上限。

將上述公式作為雙層兩階段調(diào)度優(yōu)化模型的約束條件，控制能源系統(tǒng)的調(diào)度過(guò)程。

1.2 主動(dòng)配電網(wǎng)下分布式能源系統(tǒng)雙層兩階段調(diào)度優(yōu)化方案求解

使用差分帝國(guó)競(jìng)爭(zhēng)算法對(duì)構(gòu)建的調(diào)度優(yōu)化模型進(jìn)行求解。在此算法的使用過(guò)程中，其最小的個(gè)體單位是國(guó)家[16]，根據(jù)相應(yīng)計(jì)算原理，可將國(guó)家視作N×1的矩陣，則有

c=[d1,d2,…,dn]

(8)

根據(jù)以往研究結(jié)果，將國(guó)家的價(jià)值與勢(shì)力通過(guò)代價(jià)函數(shù)體現(xiàn)，則有

cost=f(c)=f(d1,d2,…,dn)

(9)

在模型求解的過(guò)程中，尋找分布式能源系統(tǒng)運(yùn)行的最小成本。在計(jì)算前，需給定初始化的國(guó)家數(shù)量ncou以及原始帝國(guó)數(shù)量nemp。經(jīng)過(guò)多次計(jì)算得到初始國(guó)家后，選擇初始國(guó)家nemp中代價(jià)函數(shù)的取值最小但力量最強(qiáng)大的國(guó)家作為初始化帝國(guó)主義國(guó)家，剩余數(shù)據(jù)作為此部分帝國(guó)主義國(guó)家的殖民地。

每個(gè)國(guó)家對(duì)應(yīng)的殖民地個(gè)數(shù)的具體計(jì)算過(guò)程如下。

(10)

(11)

N(Yn)=round{dn，N}

(12)

式中：yn為第n個(gè)帝國(guó)類(lèi)型國(guó)家的代價(jià)函數(shù)取值。

使用式(10)將其轉(zhuǎn)化為標(biāo)準(zhǔn)代價(jià)取值結(jié)果，而后再將標(biāo)準(zhǔn)代價(jià)取值結(jié)果代入式(11)中，得到各個(gè)國(guó)家的勢(shì)力大小標(biāo)準(zhǔn)值dn。最后，使用式(12)得到各個(gè)國(guó)家的初始殖民地個(gè)數(shù)。在分布式能源系統(tǒng)的調(diào)度過(guò)程中，更具有優(yōu)勢(shì)的方案會(huì)通過(guò)整合較差的調(diào)度方案提高此方法的使用效果[17]，在帝國(guó)競(jìng)爭(zhēng)算法中通過(guò)競(jìng)爭(zhēng)機(jī)制模擬此過(guò)程。在研究中，帝國(guó)的總成本表達(dá)式為

U(Yn)=f(empn)+φmean[f(cempn)]

(13)

式(13)中：第n個(gè)帝國(guó)的總代價(jià)使用U(Yn)表示，帝國(guó)總體勢(shì)力與殖民地之間的關(guān)系通過(guò)φ體現(xiàn)。通過(guò)計(jì)算可知，當(dāng)φ取值越大時(shí)，兩者之間更易出現(xiàn)相互影響的情況，但其整體取值范圍控制在(0,1)。每個(gè)帝國(guó)都希望通過(guò)殖民掠奪擴(kuò)大版圖，并控制其他的國(guó)家。在此行為開(kāi)始后，較弱的帝國(guó)會(huì)成為競(jìng)爭(zhēng)計(jì)算的核心，則每個(gè)帝國(guó)被占領(lǐng)的可能性可表示為

NU(Yn)=U(Yn)-max[U(Yi)]

(14)

(15)

式中：U(Yn)與NU(Yn)分別為第n個(gè)國(guó)家的代價(jià)函數(shù)以及代價(jià)函數(shù)標(biāo)準(zhǔn)值；Pn為此帝國(guó)成為勢(shì)力最弱帝國(guó)的可能性。

在帝國(guó)互相競(jìng)爭(zhēng)侵占的過(guò)程中，勢(shì)必會(huì)出現(xiàn)帝國(guó)大量消亡的情況[18]。與此同時(shí)，殖民地的數(shù)量將會(huì)隨著計(jì)算迭代次數(shù)不斷減少。當(dāng)大部分帝國(guó)均消失后，所剩下的最強(qiáng)帝國(guó)所代表的方案，則為本次調(diào)度優(yōu)化模型的最優(yōu)解。其所侵占的殖民地位置將會(huì)收斂于帝國(guó)位置的周?chē)?/p>

2 算例測(cè)試分析

完成了新型分布式能源系統(tǒng)雙層兩階段調(diào)度優(yōu)化模型的構(gòu)建過(guò)程后，為驗(yàn)證模型具有一定的使用效果。在理論部分設(shè)計(jì)完成后，構(gòu)建算例測(cè)試分析環(huán)節(jié)，對(duì)此模型的應(yīng)用結(jié)果與性能加以分析。

2.1 算例參數(shù)設(shè)定

對(duì)比多個(gè)案例后，選擇IEEE-33節(jié)點(diǎn)配電網(wǎng)測(cè)試系統(tǒng)搭建測(cè)試環(huán)境，使用仿真實(shí)驗(yàn)的方式得到模型應(yīng)用結(jié)果。此系統(tǒng)拓?fù)浣Y(jié)構(gòu)為輻射狀三相平衡系統(tǒng)。根據(jù)算例測(cè)試要求，將系統(tǒng)電壓設(shè)定為11 V，可承載有功負(fù)荷閾值設(shè)定為3 000 kW，可承載無(wú)功負(fù)荷設(shè)定為2 000 kVar，測(cè)試系統(tǒng)基本結(jié)構(gòu)如圖1所示。

節(jié)點(diǎn)0為變電站節(jié)點(diǎn)，其余節(jié)點(diǎn)為工作節(jié)點(diǎn)圖1 IEEE-33節(jié)點(diǎn)配電網(wǎng)測(cè)試系統(tǒng)結(jié)構(gòu)圖Fig.1 IEEE-33 node distribution network test system structure diagram

為使其更加貼近真實(shí)主動(dòng)配電網(wǎng)應(yīng)用環(huán)境，將電源分別接入節(jié)點(diǎn)9與節(jié)點(diǎn)2，微型電網(wǎng)安裝在節(jié)點(diǎn)32與節(jié)點(diǎn)23的位置。選用微型燃?xì)廨啓C(jī)作為分布式電源，此電源額定功率為300 kW。引入系統(tǒng)中的微電網(wǎng)中包含燃?xì)廨啓C(jī)、光伏設(shè)備以及風(fēng)力機(jī)組。所研究的配電測(cè)試系統(tǒng)各支路的年最大負(fù)荷利用時(shí)間Tmax=150 h，單位售電價(jià)格ci=1.0元/(kW·h)，能源系統(tǒng)規(guī)劃年限設(shè)定為5年。此外，為提升算例的真實(shí)性，在本次測(cè)試中主要考慮日前與日內(nèi)兩種優(yōu)化調(diào)度模式。在算例測(cè)試中將上述設(shè)定結(jié)果作為測(cè)試的基礎(chǔ)，以此分析優(yōu)化模型的使用效果。

2.2 算例數(shù)據(jù)設(shè)定

引入MATLAB軟件中的3次樣條差值法，為算例分析提供數(shù)據(jù)基礎(chǔ)。在測(cè)試前，獲取日前24 h調(diào)度計(jì)劃，將日前24 h劃分為96個(gè)數(shù)據(jù)點(diǎn)，每個(gè)單位時(shí)間周期設(shè)定為15 min，則采樣時(shí)間間隔可表示為t=15 min。日內(nèi)調(diào)度時(shí)間周期與優(yōu)化時(shí)域均設(shè)定為10 min，則日內(nèi)采樣時(shí)間間隔可表示為t=5 min。與此同時(shí)，將此測(cè)試系統(tǒng)各節(jié)點(diǎn)參數(shù)設(shè)定如表1所示。

表1 IEEE-33節(jié)點(diǎn)配電網(wǎng)測(cè)試系統(tǒng)相關(guān)參數(shù)設(shè)定結(jié)果Table 1 Parameter setting results of IEEE-33 node distribution network test system

除上述數(shù)據(jù)外，將光伏發(fā)電成本設(shè)定為0.015 0元/(kW·h)；風(fēng)力發(fā)電成本設(shè)定為0.005 8元/(kW·h)；微電網(wǎng)發(fā)電成本為0.041 50元/(kW·h)；微型燃?xì)廨啺l(fā)電成本為0.005 14 元/(kW·h)。將此部分?jǐn)?shù)據(jù)作為算例測(cè)試中的初始成本參數(shù)，使用此參數(shù)輔助完成計(jì)算過(guò)程，并對(duì)不同調(diào)度方法應(yīng)用前后的成本變化加以分析。

2.3 所提方法調(diào)度前功率測(cè)試結(jié)果及分析

根據(jù)算例設(shè)定結(jié)果，得到分布式能源系統(tǒng)日前調(diào)度結(jié)果，如圖2所示。

圖2 所提方法應(yīng)用前功率測(cè)試結(jié)果Fig.2 Power test results before application of the proposed method

根據(jù)圖2的曲線(xiàn)變化趨勢(shì)可知，在0：00—6：00，整個(gè)分布式能源系統(tǒng)處于充電狀態(tài)，此時(shí)段內(nèi)能源系統(tǒng)中的儲(chǔ)能容量較低，各功率值均較低，系統(tǒng)中的剩余功率將通過(guò)聯(lián)通線(xiàn)路傳輸?shù)较到y(tǒng)的儲(chǔ)能單元，為后期的電能需求高峰期提供供給保障。在12：00左右，逐漸達(dá)到用電高峰期，配電網(wǎng)的功率大幅度提升，使用提出的調(diào)度優(yōu)化模型可以很好地實(shí)現(xiàn)功率的消納，并充分利用此部分能量完成其他工作項(xiàng)目。

2.4 所提方法調(diào)度后功率測(cè)試結(jié)果及分析

根據(jù)已設(shè)定的算例參數(shù)以及數(shù)據(jù)，得到主動(dòng)配電網(wǎng)分布式能源系統(tǒng)日內(nèi)調(diào)度結(jié)果，具體如圖3所示。

圖3 所提方法應(yīng)用后功率測(cè)試結(jié)果Fig.3 Power test results after application of the proposed method

根據(jù)圖3可知，所提方法應(yīng)用后功率總和得以下降。由此數(shù)據(jù)進(jìn)行分析可知，日內(nèi)調(diào)度在一定程度上縮小了全天電力需求的最大值與最小值之間的差距，使得能源系統(tǒng)電能輸出更加穩(wěn)定。與此同時(shí)，對(duì)燃?xì)廨啓C(jī)的有效出力進(jìn)行計(jì)算可知，其出力大小發(fā)生了明顯的變化，比照日前調(diào)整階段同比增長(zhǎng)了1.5%。由此可知分布式電源的出力大小均發(fā)生了較大的變化，其利用率得到了明顯的提升。

2.5 方法應(yīng)用前后的成本對(duì)比

為了突出本文方法的實(shí)用性，對(duì)比所提方法應(yīng)用前和應(yīng)用后的光伏發(fā)電成本、風(fēng)力發(fā)電成本、微電網(wǎng)發(fā)電成本、微型燃?xì)廨啺l(fā)電成本。根據(jù)實(shí)驗(yàn)初始參數(shù)設(shè)置可知，單位功率對(duì)應(yīng)的光伏發(fā)電成本設(shè)定為0.015 0元/(kW·h)；風(fēng)力發(fā)電成本設(shè)定為0.005 8元/(kW·h)；微電網(wǎng)發(fā)電成本為0.041 50元/(kW·h)；微型燃?xì)廨啺l(fā)電成本為0.005 14 元/(kW·h)。成本的具體對(duì)比結(jié)果如表2所示。

表2 所提方法應(yīng)用前后的成本對(duì)比Table 2 Cost comparison before and after application of the proposed method

根據(jù)表2可知，所提方法應(yīng)用下，光伏發(fā)電成本、風(fēng)力發(fā)電成本、微電網(wǎng)發(fā)電成本、微型燃?xì)廨啺l(fā)電成本均明顯下降，說(shuō)明所提方法具有理想的實(shí)用性，優(yōu)化主動(dòng)配電腦網(wǎng)下分布式電源的應(yīng)用效果和經(jīng)濟(jì)效益。

2.6 算例測(cè)試結(jié)果討論

在當(dāng)前分布式能源系統(tǒng)調(diào)度優(yōu)化模型的基礎(chǔ)上進(jìn)行了研究，并提出了可分階段分層調(diào)度的分布式能源系統(tǒng)調(diào)度優(yōu)化模型。根據(jù)日前與日內(nèi)兩階段調(diào)度結(jié)果以及其分析內(nèi)容，確定研究中提出模型的使用效果。算例測(cè)試結(jié)果表明：此模型具有一定的使用價(jià)值，通過(guò)對(duì)目前階段的研究可以確定，調(diào)度優(yōu)化模型可提升分布式電源的利用率。與此同時(shí)，此模型使用后，分布式能源系統(tǒng)運(yùn)行狀態(tài)更加穩(wěn)定。由此可知，此模型可降低分布式能源運(yùn)行的波動(dòng)性以及不確定性。綜上所述，提出的新型調(diào)度優(yōu)化模型具有較高的使用價(jià)值，可將其應(yīng)用到實(shí)際的電力系統(tǒng)調(diào)度管理工作中。

3 結(jié)論

分布式能源系統(tǒng)的經(jīng)濟(jì)調(diào)度問(wèn)題一直都是電力系統(tǒng)發(fā)展研究中的關(guān)鍵問(wèn)題，其中合理的調(diào)度優(yōu)化模型是降低電力系統(tǒng)運(yùn)行成本和提升系統(tǒng)運(yùn)行性能的關(guān)鍵，為此提出了雙層兩階段調(diào)度模型。經(jīng)算例測(cè)試證實(shí)，此調(diào)度優(yōu)化模型具有較好的使用效果以及良好的應(yīng)用前景。在以后的研究中可對(duì)模型的計(jì)算部分進(jìn)行更為詳細(xì)的分析，不斷完善優(yōu)化模型結(jié)果，提高模型使用效果。