楊艷紅 ，裴 瑋 ，鄧 衛(wèi) ，肖 浩 ，齊智平

（1.中國(guó)科學(xué)院電工研究所，北京 100190；2.中國(guó)科學(xué)院大學(xué)，北京 100049）

0 引言

隨著電力系統(tǒng)負(fù)荷的不斷增長(zhǎng)，單純依靠增加裝機(jī)容量來(lái)滿足用電的需求會(huì)造成發(fā)電資源巨大的浪費(fèi)，而微電網(wǎng)作為智能電網(wǎng)的重要組成部分，可以集成可再生能源發(fā)電、燃?xì)廨啓C(jī)（GT）發(fā)電、儲(chǔ)能裝置以及部分電力負(fù)荷，從中低壓配電網(wǎng)接入電力系統(tǒng)，運(yùn)行方式靈活多樣，可以緩解電力系統(tǒng)中裝機(jī)以及輸電網(wǎng)架的建設(shè)壓力，其在未來(lái)電網(wǎng)中將發(fā)揮重要作用［1-2］。如何合理地利用有限的資源，提高微電網(wǎng)系統(tǒng)的運(yùn)行效率，力爭(zhēng)達(dá)到系統(tǒng)的優(yōu)化運(yùn)行，對(duì)于微電網(wǎng)的發(fā)展以及大規(guī)模應(yīng)用具有重要意義。

由于微電網(wǎng)盡量保證能源的本地吸收利用沒(méi)有輸電環(huán)節(jié)，其中的電源和負(fù)荷直接緊密耦合在一起，負(fù)荷的波動(dòng)直接反饋到電源，同時(shí)光伏（PV）發(fā)電、風(fēng)力發(fā)電等受到環(huán)境條件影響類型的電源輸出功率隨機(jī)波動(dòng)，這樣就造成了部分機(jī)組的頻繁啟動(dòng)和調(diào)整出力，需要合理地安排機(jī)組的啟停和出力配合。文獻(xiàn)［3-6］研究了單母線微電網(wǎng)的經(jīng)濟(jì)運(yùn)行問(wèn)題，分析了可再生能源波動(dòng)、電價(jià)政策對(duì)于微電網(wǎng)運(yùn)行的影響。文獻(xiàn)［7-8］研究了微電網(wǎng)運(yùn)行中發(fā)電成本最小、系統(tǒng)網(wǎng)損最小和環(huán)境影響費(fèi)用最低等多個(gè)目標(biāo)之間的協(xié)調(diào)優(yōu)化，采用多目標(biāo)優(yōu)化方法對(duì)微電網(wǎng)優(yōu)化運(yùn)行問(wèn)題進(jìn)行了研究。文獻(xiàn)［9-11］研究了儲(chǔ)能在微電網(wǎng)潮流優(yōu)化中的作用，表明儲(chǔ)能可以很好地平抑可再生能源和公共連接點(diǎn)的能量波動(dòng)，但是在多母線的大系統(tǒng)中配置大量?jī)?chǔ)能目前還不夠現(xiàn)實(shí)，系統(tǒng)較大的能量波動(dòng)還需要機(jī)組之間的配合來(lái)解決。

上述研究都沒(méi)有涉及微電網(wǎng)中多機(jī)組的協(xié)調(diào)配合運(yùn)行問(wèn)題，目前針對(duì)小機(jī)組的機(jī)組組合UC（Unit Commitment）優(yōu)化研究也相對(duì)較少。文獻(xiàn)［12］使用場(chǎng)景分析的方法研究了含有可再生能源和熱電聯(lián)供機(jī)組的機(jī)組組合優(yōu)化問(wèn)題，無(wú)法保證系統(tǒng)全局最優(yōu)。文獻(xiàn)［13-14］分別使用動(dòng)態(tài)規(guī)劃和改進(jìn)遺傳算法對(duì)微電網(wǎng)機(jī)組組合優(yōu)化進(jìn)行了研究，但都針對(duì)單母線微電網(wǎng)系統(tǒng)，沒(méi)有考慮系統(tǒng)網(wǎng)絡(luò)拓?fù)浼s束。文獻(xiàn)［15］對(duì)3種混合整數(shù)線性規(guī)劃方法在機(jī)組組合問(wèn)題上的應(yīng)用進(jìn)行了對(duì)比分析。文獻(xiàn)［16-17］采用Benders分解方法對(duì)傳統(tǒng)電力系統(tǒng)中的機(jī)組組合問(wèn)題進(jìn)行了求解，但由于微電網(wǎng)中電源分散接入，可再生能源發(fā)電輸出功率不穩(wěn)定，使得系統(tǒng)潮流和機(jī)組的啟停更為復(fù)雜。

針對(duì)目前對(duì)于多母線多機(jī)組聯(lián)網(wǎng)微電網(wǎng)優(yōu)化運(yùn)行研究中的不足，本文綜合考慮了微電網(wǎng)運(yùn)行優(yōu)化中的機(jī)組組合問(wèn)題和優(yōu)化潮流OPF（Optimal Power Flow）問(wèn)題，建立了考慮機(jī)組組合的微電網(wǎng)優(yōu)化運(yùn)行模型。針對(duì)模型中十分龐大的變量和約束條件數(shù)目，并且變量可以分為0/1決策變量和機(jī)組出力連續(xù)變量，采用Benders分解的方法對(duì)2類變量進(jìn)行分解，形成無(wú)網(wǎng)絡(luò)約束機(jī)組組合主問(wèn)題和網(wǎng)絡(luò)約束優(yōu)化潮流子問(wèn)題，從而可以快速地對(duì)模型進(jìn)行求解。在改造后的IEEE 13節(jié)點(diǎn)配電系統(tǒng)上進(jìn)行了算例驗(yàn)證與仿真分析。

1 微電網(wǎng)優(yōu)化運(yùn)行問(wèn)題

微電網(wǎng)中負(fù)荷隨著時(shí)間、季節(jié)的變化而變化，光伏、風(fēng)電等不穩(wěn)定電源輸出功率也隨著時(shí)間、季節(jié)的變化而變化，負(fù)荷以及不穩(wěn)定電源的這種變化使得若只改變機(jī)組的出力而不改變機(jī)組的組合，不僅會(huì)造成資源的浪費(fèi)，而且可能難以滿足負(fù)荷變化的要求，甚至造成供電安全性、可靠性的問(wèn)題。所以，在進(jìn)行微電網(wǎng)優(yōu)化運(yùn)行時(shí)，需要將經(jīng)典優(yōu)化潮流進(jìn)行擴(kuò)展，同時(shí)引入機(jī)組組合的考慮，建立考慮機(jī)組組合的微電網(wǎng)優(yōu)化運(yùn)行模型。

經(jīng)典優(yōu)化潮流問(wèn)題主要是對(duì)電力系統(tǒng)運(yùn)行過(guò)程中的某個(gè)時(shí)間片段的靜止計(jì)算分析，需要詳盡的網(wǎng)絡(luò)模型，一般沒(méi)有考慮到發(fā)電機(jī)組的啟停問(wèn)題；而經(jīng)典機(jī)組組合問(wèn)題主要解決一個(gè)調(diào)度周期內(nèi)各時(shí)間段機(jī)組的啟停調(diào)度和機(jī)組出力，需要詳盡考慮機(jī)組自身的特性，使得調(diào)度周期內(nèi)發(fā)電成本最小，通常不考慮網(wǎng)絡(luò)約束；文中以微電網(wǎng)綜合運(yùn)行成本最小為優(yōu)化目標(biāo)，建立兼顧微電網(wǎng)運(yùn)行機(jī)組組合和優(yōu)化潮流的數(shù)學(xué)優(yōu)化模型如下。

目標(biāo)函數(shù)：

約束條件如下。

a.功率平衡約束。保證系統(tǒng)中每個(gè)節(jié)點(diǎn)的輸入功率總和等于輸出功率總和。由于微電網(wǎng)中電源與負(fù)荷之間的線路較短，線路損耗比較小，同時(shí)發(fā)電機(jī)組主要滿足有功需求，暫時(shí)不考慮其無(wú)功調(diào)節(jié)能力，因此采用近似直流潮流的方法分析系統(tǒng)有功分配，建立功率平衡約束如下：對(duì)于非平衡節(jié)點(diǎn)的約束如式（2）所示，等號(hào)左邊第1項(xiàng)為發(fā)電機(jī)組注入節(jié)點(diǎn)功率之和，第2項(xiàng)為與之相連線路注入功率之和，等號(hào)右邊是負(fù)荷需求。

b.可靠性備用約束。為了保證系統(tǒng)突發(fā)應(yīng)急的供電可靠性，系統(tǒng)運(yùn)行時(shí)還需要留有一定備用容量，一般為負(fù)荷容量的10%～20%，約束如下：

c.發(fā)電機(jī)組輸出功率約束。發(fā)電機(jī)組必須滿足最小負(fù)荷率要求，也不能超過(guò)額定功率，約束關(guān)系如下：

d.傳輸功率約束。2個(gè)節(jié)點(diǎn)之間傳輸?shù)墓β适艿綄?dǎo)線物理?xiàng)l件約束，不能超過(guò)最大輸電能力，約束關(guān)系如下：

e.微電網(wǎng)與外部電網(wǎng)接口功率約束。根據(jù)我國(guó)實(shí)際情況，目前聯(lián)網(wǎng)運(yùn)行的微電網(wǎng)不允許向電網(wǎng)賣電，其從電網(wǎng)吸收功率受到最大傳輸功率限制。

f.機(jī)組運(yùn)行狀態(tài)約束：

g.0/1變量約束：

其中，F(xiàn)j為機(jī)組j運(yùn)行時(shí)的固定費(fèi)用；Sj為機(jī)組j啟停費(fèi)用；Ej為機(jī)組j運(yùn)行燃耗費(fèi)用；Ee為系統(tǒng)購(gòu)電費(fèi)用；vjt為機(jī)組j時(shí)段t運(yùn)行狀態(tài)變量；yjt為機(jī)組 j時(shí)段t啟停變量；Pjt為機(jī)組j時(shí)段t有功功率；Pet為系統(tǒng)時(shí)段t與大電網(wǎng)接口功率；Bnp為線路n－p電納；θnt為節(jié)點(diǎn)n時(shí)段t電壓相角；θpt為節(jié)點(diǎn)p時(shí)段t電壓相角；Knp為線路n－p損耗系數(shù)；Dnt為節(jié)點(diǎn)n時(shí)段t負(fù)荷需求；Rt為時(shí)段t系統(tǒng)備用容量；Cnp為線路n－p傳輸容量限制；分別為機(jī)組j的最小和最大出力限制；Pemax為系統(tǒng)與大電網(wǎng)接口功率限制；N為機(jī)組數(shù)；T為時(shí)段數(shù)；Mn為節(jié)點(diǎn)n的機(jī)組個(gè)數(shù)；p為與節(jié)點(diǎn)n相連接的節(jié)點(diǎn)。

式（1）—（8）所建立的微電網(wǎng)優(yōu)化運(yùn)行模型本質(zhì)上是混合整數(shù)規(guī)劃模型，模型中0/1決策變量個(gè)數(shù)較少，而與運(yùn)行時(shí)間相關(guān)的連續(xù)運(yùn)行變量數(shù)目很多，固定其中的0/1變量以后，余下的優(yōu)化問(wèn)題可以快速求解，因此考慮采用分解的方法求解這類復(fù)雜模型。

2 Benders分解方法

Benders分解方法是J.F.Benders在1962年首先提出的，目的是用于解決混合整數(shù)規(guī)劃問(wèn)題，即連續(xù)變量與整數(shù)變量同時(shí)出現(xiàn)的極值問(wèn)題［18］。Benders分解廣泛應(yīng)用于帶有決策變量和連續(xù)運(yùn)行變量的混合規(guī)劃問(wèn)題，利用該問(wèn)題中決策變量和連續(xù)運(yùn)行變量的不同特點(diǎn)及相互關(guān)系，將其分解為決策主問(wèn)題和運(yùn)行子問(wèn)題，并行迭代求解，按普通的優(yōu)化問(wèn)題求解方法（如原-對(duì)偶內(nèi)點(diǎn)法）求解運(yùn)行子問(wèn)題后，將其解信息通過(guò)Benders割返回到主問(wèn)題，再對(duì)主問(wèn)題求解并判斷是否達(dá)到最優(yōu)化。應(yīng)用問(wèn)題不同，Benders的分解形式也不同，其核心就是把一系列的子問(wèn)題解，通過(guò)Bender割的形式返回到主問(wèn)題，從而使子問(wèn)題和主問(wèn)題分開(kāi)并行求解，此種方法已經(jīng)在電力系統(tǒng)中有部分應(yīng)用［16-17］。

Benders分解的基本原理如下，對(duì)于一個(gè)如下的混合規(guī)劃問(wèn)題P：

其中，變量y是復(fù)雜變量，整個(gè)約束矩陣中都含有y，當(dāng)y被固定后，變量x的子集是完全獨(dú)立的；A為系數(shù)矩陣；c、b 為系數(shù)向量；f（y）、F（y）為變量 y 的函數(shù)；S為y的可行集。

應(yīng)用線性規(guī)劃的對(duì)偶原理，求解P1問(wèn)題的對(duì)偶問(wèn)題P2：

其中，u為問(wèn)題P1的對(duì)偶變量。

通過(guò)求解子問(wèn)題P2后，原問(wèn)題P便可以寫(xiě)成如下問(wèn)題P31：

再通過(guò)變換將P31寫(xiě)成如下極點(diǎn)形式的問(wèn)題P32：

其中，ui（i=1，2，…，m）為對(duì)偶問(wèn)題 P2 的 m 個(gè)極點(diǎn)。

引入任意一個(gè)新的變量z，上述主問(wèn)題可以寫(xiě)成如下問(wèn)題P3：

給定下界：LB=min z。

綜上所述，應(yīng)用Benders分解方法，將原始優(yōu)化問(wèn)題P分解為子問(wèn)題P2和主問(wèn)題P3，通過(guò)反復(fù)迭代過(guò)程，修改系統(tǒng)的上、下限UB和LB，當(dāng)UB=LB時(shí)，得到最優(yōu)解。

3 微電網(wǎng)優(yōu)化運(yùn)行問(wèn)題的分解

應(yīng)用Benders分解方法將復(fù)雜的微電網(wǎng)優(yōu)化運(yùn)行問(wèn)題分解為主問(wèn)題和子問(wèn)題，Benders主問(wèn)題確定機(jī)組的啟停調(diào)度，子問(wèn)題求解網(wǎng)絡(luò)的優(yōu)化潮流，子問(wèn)題向主問(wèn)題傳遞某個(gè)確定啟停狀態(tài)在次優(yōu)附近的余量信息，主問(wèn)題根據(jù)此信息更新啟停調(diào)度，反復(fù)迭代直到找到最優(yōu)解。分解后主問(wèn)題為0/1混合整數(shù)線性規(guī)劃，子問(wèn)題為非線性規(guī)劃，子問(wèn)題決定運(yùn)行費(fèi)用，同時(shí)考慮了傳輸限制和網(wǎng)絡(luò)損耗。應(yīng)用Benders分解方法求解微電網(wǎng)優(yōu)化運(yùn)行問(wèn)題的基本原理如圖1所示。

圖1 Benders方法求解微電網(wǎng)優(yōu)化運(yùn)行問(wèn)題的基本原理Fig．1 Basic principle of microgrid operation optimization based on Benders decomposition

子問(wèn)題和主問(wèn)題的數(shù)學(xué)形式如下。

a.子問(wèn)題在第w次迭代時(shí)的模型如下。

目標(biāo)函數(shù)：

約束條件：

b.主問(wèn)題在第w次迭代時(shí)的模型如下。

目標(biāo)函數(shù)：

約束條件：

式（22）為Benders割的形式，這些割提供了子問(wèn)題中總體運(yùn)行費(fèi)用的一個(gè)較低估計(jì)，每次迭代增加一個(gè)割。式（23）和（24）2個(gè)約束稱為可行割，強(qiáng)制主問(wèn)題產(chǎn)生的解滿足負(fù)載和備用容量的需求，從而保證子問(wèn)題的可行性。并且，上述模型中只有α（w）一個(gè)實(shí)數(shù)變量。

通過(guò)子問(wèn)題與主問(wèn)題間的循環(huán)迭代求解流程如圖2所示。

綜上所述，使用Benders分解的方法可以成功將微電網(wǎng)優(yōu)化運(yùn)行問(wèn)題拆解成經(jīng)典的機(jī)組組合問(wèn)題和優(yōu)化潮流問(wèn)題，從而再分別進(jìn)行求解，降低系統(tǒng)計(jì)算復(fù)雜度，能夠快速收斂到系統(tǒng)最優(yōu)值。

4 算例分析

算例選取經(jīng)過(guò)改造后的IEEE 13節(jié)點(diǎn)［19］配電系統(tǒng)進(jìn)行優(yōu)化分析，改造后的系統(tǒng)添加了1臺(tái)索拉透平公司的Saturn 20型小型燃?xì)廨啓C(jī)，額定功率輸出1210 kW；添加3臺(tái)科勒公司的燃?xì)鈨?nèi)燃機(jī)，額定功率輸出分別為 800 kW（GG1）、480 kW（GG2）、480 kW（GG3）；在節(jié)點(diǎn)634加入400 kW裝機(jī)容量的風(fēng)力發(fā)電機(jī)組（WT），在節(jié)點(diǎn)675接入100 kW裝機(jī)容量的光伏發(fā)電。改造后的系統(tǒng)如圖3所示，機(jī)組運(yùn)行數(shù)據(jù)如表1所示，風(fēng)力發(fā)電機(jī)組和光伏發(fā)電模擬功率輸出如圖4所示。標(biāo)準(zhǔn)IEEE 13節(jié)點(diǎn)配電系統(tǒng)為不平衡系統(tǒng)，本算例將其改造為平衡系統(tǒng)，線路參數(shù)取原來(lái)參數(shù)的正序分量，系統(tǒng)容量基準(zhǔn)1000 kV·A，電壓基準(zhǔn)4.16 kV，系統(tǒng)有功和線路電抗分布見(jiàn)表2，表中有功功率與電抗為標(biāo)幺值。

圖2 Benders分解方法的求解流程圖Fig.2 Flowchart of microgrid operation optimization based on Benders decomposition

圖3 改造后的IEEE 13節(jié)點(diǎn)系統(tǒng)Fig.3 Modified IEEE 13-bus system

表1 機(jī)組運(yùn)行數(shù)據(jù)Table 1 Operational data of generator

圖4 風(fēng)力發(fā)電、光伏發(fā)電模擬輸出Fig.4 Simulative outputs of wind power and PV power

表2 系統(tǒng)有功和線路電抗分布Table 2 Distributions of system active power and transmission line reactance

選取典型辦公類負(fù)荷對(duì)上述系統(tǒng)運(yùn)行時(shí)的機(jī)組組合和有功分配進(jìn)行優(yōu)化計(jì)算，系統(tǒng)總體日負(fù)荷曲線和不同節(jié)點(diǎn)日負(fù)荷空間分布情況如圖5、6所示，系統(tǒng)白天峰值負(fù)荷大約是夜晚低谷負(fù)荷的3～4倍，節(jié)點(diǎn)671和節(jié)點(diǎn)675的負(fù)荷水平明顯高于其他節(jié)點(diǎn)，負(fù)荷總體呈現(xiàn)多峰多谷的空間分布。為了便于計(jì)算，對(duì)負(fù)荷進(jìn)行了分段線性化處理。

圖5 系統(tǒng)總體日負(fù)荷曲線Fig.5 Overall daily load curve of system

圖6 不同節(jié)點(diǎn)的日負(fù)荷空間分布Fig.6 Spatial distribution of power load among different nodes

為了計(jì)算分析系統(tǒng)運(yùn)行時(shí)的經(jīng)濟(jì)特性，天然氣價(jià)格取為2.05元/m3；微電網(wǎng)與外部電網(wǎng)接口電價(jià)取為峰-平-谷 3 段電價(jià)，分別為 1.2、0.8、0.4 元/（kW·h）；微電網(wǎng)與外部電網(wǎng)接口交換功率最大為其峰值負(fù)荷的20%。由于燃?xì)廨啓C(jī)發(fā)電效率低，為了提高其運(yùn)行經(jīng)濟(jì)性，文中考慮了其以熱電聯(lián)供方式運(yùn)行，燃?xì)廨啓C(jī)以發(fā)電調(diào)度為核心，不參與熱能的平衡。因此，將其熱效率折算到了發(fā)電成本考核當(dāng)中，具體而言，燃?xì)廨啓C(jī)余熱回收效率取其發(fā)電效率的2倍，回收的熱能價(jià)格取平均電價(jià)的1/3，因此折算后的燃?xì)廨啓C(jī)發(fā)電效率為其原始發(fā)電效率的1.67倍。

應(yīng)用文中所述的Benders分解方法將機(jī)組狀態(tài)變量、啟停變量與連續(xù)運(yùn)行變量進(jìn)行分解，經(jīng)過(guò)2.4 s 10次迭代計(jì)算，目標(biāo)函數(shù)收斂到最優(yōu)值，迭代過(guò)程如圖7所示。系統(tǒng)運(yùn)行優(yōu)化結(jié)果為：日運(yùn)行費(fèi)用總計(jì)3.448萬(wàn)元，其中機(jī)組啟停費(fèi)用0.042萬(wàn)元，約占總費(fèi)用的1.22%；從電網(wǎng)買電費(fèi)用0.59萬(wàn)元，約占總費(fèi)用的17.11%。如果系統(tǒng)按照常規(guī)方式運(yùn)行，機(jī)組GT與GG1工作在額定功率或者關(guān)機(jī)狀態(tài)，不足的功率通過(guò)啟停GG2、GG3以及通過(guò)電網(wǎng)獲得，此時(shí)的日運(yùn)行費(fèi)用總計(jì)需要3.873萬(wàn)元。優(yōu)化后的與常規(guī)未優(yōu)化的系統(tǒng)運(yùn)行經(jīng)濟(jì)性對(duì)比如圖8所示，通過(guò)優(yōu)化節(jié)省了11%的運(yùn)行費(fèi)用。

圖7 迭代收斂過(guò)程Fig．7 Iterative convergence curves

圖8 系統(tǒng)運(yùn)行經(jīng)濟(jì)性對(duì)比Fig．8 Comparison of system operating economy

如果采用常規(guī)的先分別進(jìn)行機(jī)組組合優(yōu)化和潮流優(yōu)化，再進(jìn)行滾動(dòng)迭代，其中機(jī)組組合優(yōu)化和潮流優(yōu)化使用拉格朗日松弛方法，需要計(jì)算6.2 s 14次迭代，最終系統(tǒng)運(yùn)行費(fèi)用的優(yōu)化結(jié)果為3.522萬(wàn)元，其與本文所述方法的迭代收斂過(guò)程對(duì)比如圖7所示?？梢?jiàn)本文所提出的Benders分解方法收斂速度更快，優(yōu)化結(jié)果也較好。

優(yōu)化后的機(jī)組啟停與出力情況如圖9所示，其中大機(jī)組GT、GG1保持高負(fù)荷率運(yùn)行，GT只在白天高負(fù)荷時(shí)段運(yùn)行，GG1承擔(dān)系統(tǒng)基荷；通過(guò)調(diào)節(jié)小機(jī)組的啟停和從電網(wǎng)吸收的功率來(lái)平衡負(fù)荷的變化。由圖中可以看出電價(jià)變化對(duì)于系統(tǒng)運(yùn)行的影響，夜間低谷電價(jià)0.4元/（kW·h）要優(yōu)于使用小機(jī)組GG2、GG3發(fā)電，優(yōu)先使用外部電網(wǎng)的功率輸入；在白天高峰電價(jià)，盡量避免從外部電網(wǎng)輸入功率，但是由于午后風(fēng)力發(fā)電輸出功率驟降，所有機(jī)組都達(dá)到最大功率輸出，仍然不夠的部分由外部電網(wǎng)輸入。

圖9 機(jī)組啟停及出力優(yōu)化結(jié)果Fig．9 Results of unit startup /shutdown and power output optimization

為了檢驗(yàn)系統(tǒng)運(yùn)行狀態(tài)是否滿足要求，以系統(tǒng)峰值負(fù)荷（t=11）為代表進(jìn)行時(shí)間斷面分析，檢查系統(tǒng)各支路上的傳輸功率是否越限，此時(shí)系統(tǒng)的有功潮流分布如圖10所示。

圖10 系統(tǒng)有功潮流分布Fig．10 Distribution of system active power flow

由圖 10 可見(jiàn)，支路 632-671、671-692、692-675線路傳輸?shù)墓β瘦^大，這是由于節(jié)點(diǎn)671和節(jié)點(diǎn)675的負(fù)荷水平較高，但沒(méi)有超過(guò)線路傳輸功率的限制，系統(tǒng)峰值負(fù)荷運(yùn)行時(shí)滿足要求，其他時(shí)間斷面也不會(huì)出現(xiàn)過(guò)負(fù)荷情況。

5 結(jié)論

本文建立了考慮機(jī)組組合時(shí)微電網(wǎng)的優(yōu)化運(yùn)行模型，綜合優(yōu)化了系統(tǒng)中網(wǎng)絡(luò)有功潮流分配和發(fā)電機(jī)組的啟停配合。針對(duì)模型中含有大量的0/1決策變量和機(jī)組出力連續(xù)變量，研究了采用Benders分解的方法對(duì)模型進(jìn)行求解。通過(guò)在改造后的IEEE 13節(jié)點(diǎn)系統(tǒng)上進(jìn)行算例的計(jì)算與仿真分析，結(jié)果表明采用本文所述的優(yōu)化方法可以顯著提高微電網(wǎng)運(yùn)行時(shí)的經(jīng)濟(jì)特性。

文中采用近似直流潮流的計(jì)算方法來(lái)分析系統(tǒng)有功分配和建立功率平衡約束，為了進(jìn)一步精確優(yōu)化微電網(wǎng)運(yùn)行，下一步研究工作可以將文中的近似直流潮流部分改為采用交流潮流。