一種交直流混合微網(wǎng)的優(yōu)化運行方法
【專利摘要】一種交直流混合微網(wǎng)的優(yōu)化運行方法，針對含有多種微源的混合微網(wǎng)，建立了計及經(jīng)濟成本和環(huán)境效益的交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型?；诮恢绷骰旌衔⒕W(wǎng)的同時含有交流區(qū)域與直流區(qū)域的結(jié)構(gòu)特性，通過協(xié)同進化改進黑洞算法將復(fù)雜的混合微網(wǎng)拆分為兩個相對對立的子區(qū)域進行協(xié)調(diào)計算，可有效解決變量較多、結(jié)構(gòu)較復(fù)雜的交直流混合微網(wǎng)的求解問題。以包含了風、光、燃料電池、蓄電池和微燃機的具體交直流混合微網(wǎng)為實施例，進行分析。本發(fā)明的一種交直流混合微網(wǎng)的優(yōu)化運行方法，從微網(wǎng)的結(jié)構(gòu)特性出發(fā)，基于協(xié)同進化框架對經(jīng)濟運行問題進行分拆和協(xié)調(diào)，同時使用改進黑洞搜索算法進行子區(qū)尋優(yōu)，降低了交直流混合微網(wǎng)的經(jīng)濟運行問題的求解難度。
【專利說明】
-種交直流混合微網(wǎng)的優(yōu)化運行方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明設(shè)及一種微網(wǎng)的優(yōu)化運行方法。特別是設(shè)及一種基于協(xié)同進化黑桐算法的 交直流混合微網(wǎng)的優(yōu)化運行方法。
【背景技術(shù)】
[0002] 當今能源與環(huán)境問題日益凸顯，發(fā)展利用可再生能源已成為共識。微網(wǎng)作為一種 包含可再生能源等分布式發(fā)電技術(shù)的綜合集成技術(shù)，憑借其對可再生能源的高度兼容性和 對分布式發(fā)電的靈活調(diào)控能力，在日益強調(diào)節(jié)能環(huán)保的今天得到了業(yè)內(nèi)的普遍重視。交直 流混合微網(wǎng)可綜合發(fā)揮交流微網(wǎng)與直流微網(wǎng)的互補優(yōu)勢，擁有更加廣泛的適用領(lǐng)域。
[0003] 交直流混合微網(wǎng)經(jīng)濟運行問題是一個高維度、非線性、多目標的復(fù)雜優(yōu)化問題。傳 統(tǒng)數(shù)學優(yōu)化方法及其改進算法計算精確，但通常對問題模型有很高要求，而且求解難度大； 智能算法，如遺傳算法、粒子群算法等，通過群體捜索W快速找到工程可行解，此類方法對 模型要求較低，亦可保證滿足要求的精度，但是維數(shù)較多時求解依然會變得困難。智能算法 在求解的問題維數(shù)較多時，求解難度加大，運一弱點可通過對問題合理分解并整體協(xié)調(diào)的 方法克服。

【發(fā)明內(nèi)容】

[0004] 本發(fā)明所要解決的技術(shù)問題是，提供一種能夠降低交直流混合微網(wǎng)的經(jīng)濟運行問 題求解難度的交直流混合微網(wǎng)的優(yōu)化運行方法。
[0005] 本發(fā)明所采用的技術(shù)方案是：一種交直流混合微網(wǎng)的優(yōu)化運行方法，包括如下步 驟：
[0006] 1)采集交直流混合微網(wǎng)的交流區(qū)域和直流區(qū)域的負荷信息數(shù)據(jù)、氣象信息數(shù)據(jù)， 綜合微網(wǎng)運行的歷史數(shù)據(jù)，對未來一天的數(shù)據(jù)進行預(yù)測，得到未來一天內(nèi)交直流混合微網(wǎng) 的交流區(qū)域負荷、直流區(qū)域負荷、風能和太陽能功率預(yù)測數(shù)據(jù)；
[0007] 2)統(tǒng)計交直流混合微網(wǎng)的微電源特性，建立微網(wǎng)內(nèi)所有可控型微電源的出力數(shù)學 模型；
[000引3)將交直流混合微網(wǎng)未來一天內(nèi)的經(jīng)濟運行分為24個時段，W交直流混合微網(wǎng)全 天運行燃料成本、運行維護成本、購電成本和環(huán)境成本最小為目標函數(shù)，考慮微網(wǎng)內(nèi)部的各 時段電能平衡約束、并網(wǎng)聯(lián)絡(luò)線容量約束、可控微源的出力限制、可控微源的最短啟停時間 約束、蓄電池的電量容量約束、蓄電池的功率容量約束、蓄電池的容量連續(xù)性約束和計算蓄 電池的周期始末電能平衡約束，建立交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型；
[0009] 4)基于步驟1)中的交流區(qū)域負荷、直流區(qū)域負荷、風能和太陽能功率預(yù)測數(shù)據(jù)，采 用協(xié)同進化改進黑桐優(yōu)化算法對步驟3)中的交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型進行求解， 得到各時段可控型微電源的出力及混合微網(wǎng)購電量優(yōu)化方案。
[0010] 步驟2)所述的出力數(shù)學模型包括：
[0011] (1)燃料電池的出力模型：
[0012]
[001引其中，各符號的定義如下：ffuel,FC為燃料電池的燃料成本，Cgas,FC為燃料電池使用 的燃料價格，LHVgas,FC為燃料的低熱值，Pfc為燃料電池的出力大小，郵為燃料電池的效率。
[0014] (2)微型燃氣輪機的出力模型：
[0015]
[0016] 其中，各符號的定義如下MT為微型燃氣輪機的燃料成本，Cgas,MT為微型燃氣 輪機使用的燃料價格，LHVgas,MT為微型燃氣輪機的燃料的低熱值，Pmt為微型燃氣輪機的的出 力大小，TlMT為微型燃氣輪機的的效率。
[0017] 步驟3)所述的的交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型包括目標函數(shù)和約束條件，其 中，目標函數(shù)為：
[001 引
[0019] 其中，各符號定義如下:F為目標函數(shù)，Cfuel為燃料成本，Com為運行維護成本，Cbuy為 購電成本，Cw為污染折算成本;t為時段，N為微源總個數(shù)，ft( ?)為第i個微源在t時段的燃料 成本函數(shù)，Pi(t)為第i個微源在t時段的實際輸出功率，ki為第i個微源的運行維護成本系 數(shù)，a%t時段時的購電電價，Pbuy(t)為混合微網(wǎng)在t時段時的購電量，M為污染物總數(shù)，^3J為 第巧巾污染物的環(huán)境評價標準，、為第i個機組的第巧巾污染物的排污系數(shù)；
[0020] 約束條件為：
[0021] (1)發(fā)電單元功率限額：
[0022]
[0023] (2)最短啟停時間約束：
[0024]
[0025] 若在改進黑桐算法的迭代趨于收斂時，部分機組沒有滿足啟停時間的約束，對沒 有滿足啟停時間約束的機組按如下方法加 W調(diào)整：如果機組連續(xù)啟動時間過短，則強制延 長運行時間；如果停機時間過短，則將違規(guī)停機時間強制改為W最低功率運行，再調(diào)整其他 機組出力W滿足系統(tǒng)的其他約束；
[0026] (3)蓄電池電量容量和功率容量約束：
[0027]
[002引
[0029]
[0030] (5)計算蓄電池周期始末電能平衡約束：
[0031] E(0)=E(24 化（t)
[0032] (6)全系統(tǒng)實時電能平衡約束：
[0033] EPi(t)+Pbuy(t)+扣 S(t)=時（t)+PLoss(t)
[0034] (7)并網(wǎng)聯(lián)絡(luò)線容量約束：
[0035]
[0036] 其中和/f"分別為輸出上限和下限，Pi(t)為第i個微源在t時段的實際輸出 功率;T。。為可控型微電源的連續(xù)開啟時間，Toff為可控型微電源的連續(xù)停機時間，可控型微 電源的連續(xù)開啟和停機時間的下限分別為寫產(chǎn)與：Cm;EES(t)為t時段的蓄電池的電量容 量，Emin、Emax分別為蓄電池電量容量的下限和上限;PES(t)為t時段的蓄電池輸出量，W輸出 為正方向，Pes,min、扣S,max分別為蓄電池輸出量的下限和上限；；E(t)為蓄電池 t時段的電量，E (t-1)為蓄電池 t-1時段的電量巧ESl為蓄電池充電效率巧ES2為蓄電池放電效率；AT為單位 時段，本發(fā)明中AT = 1; Pbuy (t)為交直流混合微網(wǎng)t時段的購電功率，Pd (t)為交直流混合微 網(wǎng)t時段的電負荷需求，PlDSS(t)為交直流混合微網(wǎng)t時段的有功網(wǎng)損;瑞T為聯(lián)絡(luò)線容量約 束下的最大購電功率。
[0037] 步驟4)所述的對交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型進行求解包括：
[0038] (1)初始化:輸入微電源、負荷參數(shù);輸入算法參數(shù)包括星體范圍、最大尋優(yōu)次數(shù)、 星體的個數(shù)、新產(chǎn)生星體距離黑桐的最小距離和星體的維數(shù);各時段交、直流區(qū)域間的交互 功率;交、直流區(qū)域間的交互功率計算迭代次數(shù)PP=I;設(shè)定交、直流區(qū)域間的交互功率改變 量e ;
[0039] (2)將交直流混合微網(wǎng)的各優(yōu)化變量根據(jù)位置劃分到交流區(qū)域和直流區(qū)域中；
[0040] (3)將交互功率計入各區(qū)域功率平衡約束中，結(jié)合步驟3)所述的約束條件，給出功 率平衡約束如下：
[0041]
[0042] 其中，各符號定義如下:PAC(t)為t時段交流區(qū)域微電源的發(fā)電功率輸出總量，包括 購電功率；巧f的為t時段交流區(qū)域蓄電池的輸出功率；巧為t時段交流區(qū)域負荷； 巧f如為t時段交流區(qū)域網(wǎng)損;PM(t)為t時段交、直流區(qū)域的交互功率;巧1(0為t時段交、直 流區(qū)域間的雙向換流器的損耗;UO為雙向換流器的損耗系數(shù)，當PM(t)〉0時，110=1，否則UO = 0;PDC(t)為t時段直流區(qū)域微電源的發(fā)電功率輸出總量;巧f的為t時段直流區(qū)域蓄電池的 輸出功率;巧f(〇為t時段直流區(qū)域負荷；記非)為t時段直流區(qū)域網(wǎng)損；
[0043] (4)設(shè)置改進黑桐算法的迭代次數(shù)Uer = I;
[0044] (5)在區(qū)域P中隨機產(chǎn)生N個星體；
[0045] (6)計算初始星體的自適應(yīng)度值，選擇最優(yōu)的自適應(yīng)度值確定初始黑桐；
[0046] (7)更新星體的位置；
[0047] (8)計算更新后星體的自適應(yīng)度值，若星體的自適應(yīng)度值優(yōu)于初始黑桐的自適應(yīng) 度值，交換所述星體與初始黑桐的位置，所述星體為新的黑桐，其他星體向著新的黑桐移動 更新，否則繼續(xù)靠近初始黑桐位置更新星體的位置，星體向著黑桐移動更新的速度和慣性 為：
[004引
[0049] 其中，各符號定義如下:Vit為星體t時刻的速度，Vit^i為星體t-1時刻的速度，定義t =1時速度為0 ; XiP為第i個星體的歷史最優(yōu)適應(yīng)值;Xit為第i個星體在t時刻的位置;祉H為黑 桐的位置;ri和n為[0山之間的隨機數(shù)；W為星體向著黑桐移動的慣性參量，Wmin、COmax為 星體向著黑桐移動的慣性參量的下限和上限，取Wmax = O.9, COmin = O.4;Niter為當前改進黑 桐算法的迭代次數(shù)，Nitermax為改進黑桐算法的最大迭代次數(shù)；
[0050] 星體向著黑桐移動的速度約束為：
[0化1 ]
[0052] 其中，各符號定義如下:Vid為星體i向著黑桐移動的速度，Vidmin與Vidmax分別為星體 向著黑桐移動的最小與最大速度;Xidmax、Xidmin分別為星體i向著黑桐移動一次的最小、最大 距離;定值參數(shù)L根據(jù)算例取定，用來限制速度的幅度；
[0053] (9)判斷星體是否符合吞隧條件，若符合，黑桐將星體吞隧，并在可行空間內(nèi)隨機 位置重新生成新星體，在可行空間內(nèi)隨機位置重新生成新星體，具體是采用Logistic映射 產(chǎn)生混濁變量，映射公式如下：
[0054] 化+I=Mk(I-Xk)
[0055] 其中，各符號定義如下:祉、xk+i為星體在第k、k+l次尋優(yōu)中的位置;y為常數(shù)，取4令 映射進入混濁狀態(tài);所述映射有不動點0、0.25、0.50、0.75、1，若星體的初始位置為所述的 不動點，則選用普通隨機數(shù)生成新星體的位置；
[0056] 若不符合吞隧條件，星體繼續(xù)向靠近黑桐位置移動更新星體的位置；
[0057] (10)判斷當前的迭代次數(shù)是否達到最大迭代次數(shù)，若達到則結(jié)束計算，輸出計算 結(jié)果，進入步驟（11);若不滿足則設(shè)定改進黑桐算法的迭代次數(shù)iter=iter+l，返回第(7) 步；
[0058] (11)交、直流子區(qū)域分別產(chǎn)生若干個最優(yōu)個體；
[0059] (12)將交、直流子區(qū)提交的個體互相組合，得到多種全網(wǎng)機組出力計劃，由于可能 存在通過跨子區(qū)的合作互補，故逐時段進行如下優(yōu)化：
[0060] min( Cfuel ( t)+Com( t)+Cbuy ( t)+Cw( t))
[0061 ]其中，各符號定義如下:Cfuel(t)為t時段的燃料成本，Com(t)為t時段的運行維護成 本，Cbuy(t)為t時段的購電成本，Cw(t)為t時段的污染折算成本，約束條件與步驟3)-致；
[0062] (13)將所有24個時段都進行第(12)步所述的優(yōu)化處理，得到全周期總成本；
[0063] (14)在第（13)步中得到的全周期總成本中找出綜合成本最低的組合，計算各時段 交、直流區(qū)域間的交互功率，當各時段交、直流區(qū)域間的交互功率的計算迭代次數(shù)PP = I時， 若各時段交、直流區(qū)域間的交互功率相對于初始的各時段交、直流區(qū)域間的交互功率改變 量均值小于給定的e，則轉(zhuǎn)第（15)步；當各時段交、直流區(qū)域間的交互功率的計算迭代次數(shù) PP>1時，若各時段交、直流區(qū)域間的交互功率相對于上一次的若各時段交、直流區(qū)域間的 交互功率改變量均值小于給定的6,則轉(zhuǎn)第(15)步；
[0064] 否則將綜合成本最低的組合作為新的邊界傳輸功率值，賦給交流、直流區(qū)域，計算 各時段交直流子區(qū)間交互功率，公式如下：
[00 化]
[0066] 其中，各符號定義如下:PM(t)為t時段交、直流區(qū)域的交互功率;PAC(t)為t時段交 流區(qū)域微電源的發(fā)電功率輸出總量，包括購電功率;巧f價為t時段交流區(qū)域蓄電池的輸出 功率;巧產(chǎn)的為*時段交流區(qū)域負荷;自的為t時段交流區(qū)域網(wǎng)損；巧:,的為t時段交、直流 區(qū)域間的雙向換流器的損耗;UO為雙向換流器的損耗系數(shù)，當PM(t)〉0時，110=1，否則UO = O; 交、直流區(qū)域間的交互功率計算迭代次數(shù)PP = PP+1，轉(zhuǎn)第(3)步繼續(xù)計算；
[0067] (15)停止計算，輸出結(jié)果。
[0068] 本發(fā)明的一種交直流混合微網(wǎng)的優(yōu)化運行方法，從微網(wǎng)的結(jié)構(gòu)特性出發(fā)，基于協(xié) 同進化框架對經(jīng)濟運行問題進行分拆和協(xié)調(diào)，同時使用改進黑桐捜索算法進行子區(qū)尋優(yōu)， 降低了交直流混合微網(wǎng)的經(jīng)濟運行問題的求解難度。本發(fā)明具有的有益效果是：
[0069] (1)根據(jù)交直流混合微網(wǎng)的結(jié)構(gòu)特性，引入?yún)f(xié)同進化框架，將復(fù)雜的交直流混合微 網(wǎng)優(yōu)化運行問題拆分成兩個相對簡單的交流區(qū)域優(yōu)化問題和直流區(qū)域優(yōu)化問題，降低了問 題求解的難度。
[0070] (2)將復(fù)雜的交直流混合微網(wǎng)優(yōu)化運行問題運用協(xié)同進化框架進行拆分后，運用 改進黑桐算法對其進行求解。由于傳統(tǒng)的黑桐算法容易陷入局部最優(yōu)而導(dǎo)致算法早熟，且 在產(chǎn)生新星體的時候如果直接按原始方法隨機生成新的星體常常不滿足約束。為改善運些 缺點，在生成新星體時加入混濁機制、加入星體被吸引的慣性和速度約束，并且對不滿足啟 停約束的機組調(diào)整其啟停時間，提升了求解效率。
【附圖說明】
[0071 ]圖1是本發(fā)明一種交直流混合微網(wǎng)的優(yōu)化運行方法的流程圖；
[0072] 圖2是本發(fā)明交直流混合微網(wǎng)結(jié)構(gòu)圖；
[0073] 圖3是本發(fā)明實施例地區(qū)典型日風力發(fā)電機組、光伏電池發(fā)電預(yù)、交流負荷、直流 負荷測及分時電價曲線圖；
[0074] 圖4是本發(fā)明中可控微源的燃料成本-機端功率曲線圖；
[0075] 圖5是本發(fā)明優(yōu)化后交直流混合微網(wǎng)交流區(qū)24個時段的優(yōu)化結(jié)果圖；
[0076] 圖6是本發(fā)明優(yōu)化后交直流混合微網(wǎng)直流區(qū)24個時段的優(yōu)化結(jié)果圖；
[0077] 圖7是本發(fā)明優(yōu)化后交直流混合微網(wǎng)整體24個時段的優(yōu)化結(jié)果圖；
[0078] 圖8是本發(fā)明優(yōu)化后交直流混合微網(wǎng)交流區(qū)與直流區(qū)24個時段的交互功率優(yōu)化結(jié) 果圖。
【具體實施方式】
[0079] 下面結(jié)合實施例和附圖對本發(fā)明的一種交直流混合微網(wǎng)的優(yōu)化運行方法做出詳 細說明。
[0080] 本發(fā)明的一種交直流混合微網(wǎng)的優(yōu)化運行方法要克服傳統(tǒng)方法及智能算法在復(fù) 雜問題求解上的弊端，將協(xié)同進化法與黑桐算法相結(jié)合，提出協(xié)同進化黑桐算法。
[0081] 本發(fā)明W交直流混合微網(wǎng)的經(jīng)濟成本和環(huán)境成本最小為目標，考慮了交直流混合 微網(wǎng)內(nèi)部的各時段電能平衡約束、并網(wǎng)聯(lián)絡(luò)線容量約束、可控微源的出力限制、可控微源的 最短啟停時間約束、蓄電池的電量容量約束、蓄電池的功率容量約束、蓄電池的容量連續(xù)性 約束、蓄電池的計算周期始末電能平衡約束等，建立交直流混合微網(wǎng)優(yōu)化運行模型，并提出 一種協(xié)同進化黑桐算法，求解各時段微電源出力和交、直流兩區(qū)域的交互功率。W包含了 風、光、燃料電池、蓄電池和微燃機的具體交直流混合微網(wǎng)為實施例，針對具體實施例，對優(yōu) 化結(jié)果進行了分析。
[0082] 如圖1所示，本發(fā)明的一種交直流混合微網(wǎng)的優(yōu)化運行方法，包括如下步驟：
[0083] 1)采集交直流混合微網(wǎng)的交流區(qū)域和直流區(qū)域的負荷信息數(shù)據(jù)、氣象信息數(shù)據(jù)， 綜合微網(wǎng)運行的歷史數(shù)據(jù)，對未來一天的數(shù)據(jù)進行預(yù)測，得到未來一天內(nèi)交直流混合微網(wǎng) 的交流區(qū)域負荷、直流區(qū)域負荷、風能和太陽能功率預(yù)測數(shù)據(jù)；
[0084] 2)統(tǒng)計交直流混合微網(wǎng)的微電源特性，建立微網(wǎng)內(nèi)所有可控型微電源的出力數(shù)學 模型；
[0085] 所述的出力數(shù)學模型包括：
[0086] (1)燃料電池的出力模型：
[0087]
[008引其中，各符號的定義如下FC為燃料電池的燃料成本，Cgas,FC為燃料電池使用 的燃料價格，LHVgas,FC為燃料的低熱值，Pfc為燃料電池的出力大小，郵為燃料電池的效率。
[0089] (2)微型燃氣輪機的出力模型：
[0090]
[0091] 其中，各符號的定義如下：ffuel,MT為微型燃氣輪機的燃料成本，Cgas,MT為微型燃氣 輪機使用的燃料價格，LHVgas,MT為微型燃氣輪機的燃料的低熱值，Pmt為微型燃氣輪機的的出 力大小，riMT為微型燃氣輪機的的效率。
[0092] 3)將交直流混合微網(wǎng)未來一天內(nèi)的經(jīng)濟運行分為24個時段，W交直流混合微網(wǎng)全 天運行燃料成本、運行維護成本、購電成本和環(huán)境成本最小為目標函數(shù)，考慮微網(wǎng)內(nèi)部的各 時段電能平衡約束、并網(wǎng)聯(lián)絡(luò)線容量約束、可控微源的出力限制、可控微源的最短啟停時間 約束、蓄電池的電量容量約束、蓄電池的功率容量約束、蓄電池的容量連續(xù)性約束和計算蓄 電池的周期始末電能平衡約束，建立交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型；
[0093] 所述的的交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型包括目標函數(shù)和約束條件，其中，目 標函數(shù)為：
[0094] C=I J-L
[00M]其中，各符號定義如下:F為目標函數(shù)，Cfuel為燃料成本，Com為運行維護成本，Cbuy為 購電成本，Cw為污染折算成本;t為時段，N為微源總個數(shù)，ft( ?)為第i個微源在t時段的燃料 成本函數(shù)，Pi(t)為第i個微源在t時段的實際輸出功率，ki為第i個微源的運行維護成本系 數(shù)，a%t時段時的購電電價，Pbuy(t)為混合微網(wǎng)在t時段時的購電量，M為污染物總數(shù)，^3J為 第巧巾污染物的環(huán)境評價標準，、為第i個機組的第巧巾污染物的排污系數(shù)；
[0096] 約束條件為：
[0097] (1)發(fā)電單元功率限額：
[009引
[0099] (2)最短啟停時間約束：
[0100]
[0101] 若在改進黑桐算法的迭代趨于收斂時，部分機組沒有滿足啟停時間的約束，對沒 有滿足啟停時間約束的機組按如下方法加 W調(diào)整：如果機組連續(xù)啟動時間過短，則強制延 長運行時間；如果停機時間過短，則將違規(guī)停機時間強制改為W最低功率運行，再調(diào)整其他 機組出力W滿足系統(tǒng)的其他約束；
[0102] (3)蓄電池電量容量和功率容量約束：
[0103]
[0104] (4)蓄電池容量連續(xù)性約束：
[0105]
[0106] (5)計算蓄電池周期始末電能平衡約束：
[0107] E(0)=E(24 化（t)
[0108] (6)全系統(tǒng)實時電能平衡約束：
[0109] EPi(t)+Pbuy(t)+扣 S(t)=時（t)+PLoss(t)
[0110] (7)并網(wǎng)聯(lián)絡(luò)線容量約束：
[0111]
[0112] 其中：巧和ifn分別為輸出上限和下限，Pi(t)為第i個微源在t時段的實際輸出 功率;T。。為可控型微電源的連續(xù)開啟時間，Toff為可控型微電源的連續(xù)停機時間，可控型微 電源的連續(xù)開啟和停機時間的下限分別為7；f與r。^;垃s(t)為t時段的蓄電池的電量容量， Emin、Emax分別為蓄電池電量容量的下限和上限;PES(t)為t時段的蓄電池輸出量，W輸出為正 方向，Pes,min、扣S,max分別為蓄電池輸出量的下限和上限；；E(t)為蓄電池 t時段的電量，E(t- 1)為蓄電池 t-1時段的電量;riEsi為蓄電池充電效率;riES2為蓄電池放電效率；AT為單位時 段，本發(fā)明中AT=l;Pbuy(t)為交直流混合微網(wǎng)t時段的購電功率，PD(t)為交直流混合微網(wǎng) t時段的電負荷需求，PlDss(t)為交直流混合微網(wǎng)t時段的有功網(wǎng)損;巧為聯(lián)絡(luò)線容量約束 下的最大購電功率。
[0113] 4)基于步驟1)中的交流區(qū)域負荷、直流區(qū)域負荷、風能和太陽能功率預(yù)測數(shù)據(jù)，采 用協(xié)同進化改進黑桐優(yōu)化算法對步驟3)中的交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型進行求解， 得到各時段可控型微電源的出力及混合微網(wǎng)購電量優(yōu)化方案。所述的對交直流混合微網(wǎng)優(yōu) 化運行數(shù)學模型進行求解包括：
[0114] (1)初始化:輸入微電源、負荷參數(shù);輸入算法參數(shù)包括星體范圍、最大尋優(yōu)次數(shù)、 星體的個數(shù)、新產(chǎn)生星體距離黑桐的最小距離和星體的維數(shù);各時段交、直流區(qū)域間的交互 功率;交、直流區(qū)域間的交互功率計算迭代次數(shù)PP=I;設(shè)定交、直流區(qū)域間的交互功率改變 量e;
[0115] (2)將交直流混合微網(wǎng)的各優(yōu)化變量根據(jù)位置劃分到交流區(qū)域和直流區(qū)域中；
[0116] (3)將交互功率計入各區(qū)域功率平衡約束中，結(jié)合步驟3)所述的約束條件，給出功 率平衡約束如下：
[0117]
[0118] 其中，各符號定義如下:PAC(t)為t時段交流區(qū)域微電源的發(fā)電功率輸出總量，包括 購電功率；巧f的為t時段交流區(qū)域蓄電池的輸出功率；巧f(。為t時段交流區(qū)域負荷； 巧己的為t時段交流區(qū)域網(wǎng)損;PM(t)為t時段交、直流區(qū)域的交互功率；巧1(0為t時段交、直 流區(qū)域間的雙向換流器的損耗;UO為雙向換流器的損耗系數(shù)，當PM(t)〉0時，110=1，否則UO = 〇;PDc(t)為t時段直流區(qū)域微電源的發(fā)電功率輸出總量；塔(Y)為t時段直流區(qū)域蓄電池的 輸出功率;巧f(〇為t時段直流區(qū)域負荷；巧苦(0為t時段直流區(qū)域網(wǎng)損；
[0119] (4)設(shè)置改進黑桐算法的迭代次數(shù)Uer = I;
[0120] (5)在區(qū)域P中隨機產(chǎn)生N個星體；
[0121] (6)計算初始星體的自適應(yīng)度值，選擇最優(yōu)的自適應(yīng)度值確定初始黑桐；
[0122] (7)更新星體的位置；
[0123] (8)計算更新后星體的自適應(yīng)度值，若星體的自適應(yīng)度值優(yōu)于初始黑桐的自適應(yīng) 度值，交換所述星體與初始黑桐的位置，所述星體為新的黑桐，其他星體向著新的黑桐移動 更新，否則繼續(xù)靠近初始黑桐位置更新星體的位置，星體向著黑桐移動更新的速度和慣性 為：
[0124]
[0125] 其中，各符號定義如下:Vit為星體t時刻的速度，Vit-I為星體t-1時刻的速度，定義t =1時速度為0 ; XiP為第i個星體的歷史最優(yōu)適應(yīng)值;Xit為第i個星體在t時刻的位置;祉H為黑 桐的位置;ri和n為[(U]之間的隨機數(shù)；《為星體向著黑桐移動的慣性參量，Wmin、COmax為 星體向著黑桐移動的慣性參量的下限和上限，取《 max = 0.9，O min = 0.4 ; Niter為當前改進黑 桐算法的迭代次數(shù)，Nitermax為改進黑桐算法的最大迭代次數(shù)；
[0126] 星體向著黑桐移動的速度約束為：
[0127]
[01%]其中，各符號定義如下:Vid為星體i向著黑桐移動的速度，Vidmin與Vidmax分別為星體 向著黑桐移動的最小與最大速度;Xidmax、Xidmin分別為星體i向著黑桐移動一次的最小、最大 距離;定值參數(shù)L根據(jù)算例取定，用來限制速度的幅度，本發(fā)明的實施例中，L = 5;
[0129] (9)判斷星體是否符合吞隧條件，若符合，黑桐將星體吞隧，并在可行空間內(nèi)隨機 位置重新生成新星體，僅使用隨機數(shù)生成新星體的遍歷性不夠理想，使用混濁映射在可行 空間內(nèi)隨機位置重新生成新星體，具體是采用Logistic映射產(chǎn)生混濁變量，映射公式如下：
[0130] 化+I=Ji 化 Q-Xk)
[0131] 其中，各符號定義如下:祉、xk+i為星體在第k、k+l次尋優(yōu)中的位置;y為常數(shù)，取4令 映射進入混濁狀態(tài);所述映射有不動點0、0.25、0.50、0.75、1，若星體的初始位置為所述的 不動點，則選用普通隨機數(shù)生成新星體的位置；
[0132] 若不符合吞隧條件，星體繼續(xù)向靠近黑桐位置移動更新星體的位置；
[0133] (10)判斷當前的迭代次數(shù)是否達到最大迭代次數(shù)，若達到則結(jié)束計算，輸出計算 結(jié)果，進入步驟（11);若不滿足則設(shè)定改進黑桐算法的迭代次數(shù)iter=iter+l，返回第(7) 步；
[0134] (11)交、直流子區(qū)域分別產(chǎn)生若干個最優(yōu)個體；
[0135] (12)將交、直流子區(qū)提交的個體互相組合，得到多種全網(wǎng)機組出力計劃，由于可能 存在通過跨子區(qū)的合作互補，故逐時段進行如下優(yōu)化：
[0136] min( Cfuel ( t)+Com( t)+Cbuy ( t)+Cw( t))
[0137] 其中，各符號定義如下:Cfuel(t)為t時段的燃料成本，Com(t)為t時段的運行維護成 本，Cbuy(t)為t時段的購電成本，Cw(t)為t時段的污染折算成本，約束條件與步驟3)-致;所 述的優(yōu)化只對時段t的可調(diào)機組進行計算，不改變機組啟停狀態(tài)和儲能的功率，因此該步的 優(yōu)化問題變量和約束都較少。
[0138] (13)將所有24個時段都進行第(12)步所述的優(yōu)化處理，得到全周期總成本；
[0139] (14)在第（13)步中得到的全周期總成本中找出綜合成本最低的組合，計算各時段 交、直流區(qū)域間的交互功率，當各時段交、直流區(qū)域間的交互功率的計算迭代次數(shù)PP = I時， 若各時段交、直流區(qū)域間的交互功率相對于初始的各時段交、直流區(qū)域間的交互功率改變 量均值小于給定的e，則轉(zhuǎn)第（15)步；當各時段交、直流區(qū)域間的交互功率的計算迭代次數(shù) PP>1時，若各時段交、直流區(qū)域間的交互功率相對于上一次的若各時段交、直流區(qū)域間的 交互功率改變量均值小于給定的6,則轉(zhuǎn)第(15)步；
[0140] 否則將綜合成本最低的組合作為新的邊界傳輸功率值，賦給交流、直流區(qū)域，計算 各時段交直流子區(qū)間交互功率，公式如下：
[0141]
[0142] 其中，各符號定義如下:PM(t)為t時段交、直流區(qū)域的交互功率;PAC(t)為t時段交 流區(qū)域微電源的發(fā)電功率輸出總量，包括購電功率；的為t時段交流區(qū)域蓄電池的輸出 功率;傳K'的為t時段交流區(qū)域負荷；巧.芳的為t時段交流區(qū)域網(wǎng)損；巧1的為t時段交、直流 區(qū)域間的雙向換流器的損耗;UO為雙向換流器的損耗系數(shù)，當PM(t)〉0時，110=1，否則UO = O; 交、直流區(qū)域間的交互功率計算迭代次數(shù)PP = PP+1，轉(zhuǎn)第(3)步繼續(xù)計算；
[0143] (15)停止計算，輸出結(jié)果。
[0144] 下面給出實施例
[0145] 考慮如圖2所示的交直流混合微網(wǎng)，微網(wǎng)工作在正常狀態(tài)下。交直流混合微網(wǎng)中的 微源有微型燃氣輪機(MT),燃料電池(FC)，蓄電池化S)，風力發(fā)電(WT)，光伏發(fā)電（PV)。其 中，風電只接入交流區(qū)域，光伏只接入直流區(qū)域。大電網(wǎng)通過PCC將電能送入交流區(qū)域，限定 購電功率上限為lOOkW，下限為0。此外，交流區(qū)域配備微燃機怡、燃料電池共始、儲能裝置 1臺。直流區(qū)域配有同型號的微燃機2臺、燃料電池3臺、儲能裝置1臺。儲能裝置最大電量 ISOkWh，最小電量20kWh，最大充放電功率為40kW，充放電效率均取92 %。微燃機和燃料電池 使用天然氣作燃料，天然氣價格取2.80元/m3。
[0146] 表1各微源的參數(shù)
[0147]
[014引表2分時電價
[0149]
[0150]
[0151]
[0152]
[0153]
[0154] (I)采集交直流混合微網(wǎng)中的氣象信息數(shù)據(jù)、負荷信息數(shù)據(jù)，綜合交直流混合微網(wǎng) 運行的歷史數(shù)據(jù)，對未來一天的數(shù)據(jù)進行預(yù)測，得到未來一天內(nèi)24個時段交直流混合微網(wǎng) 的風力發(fā)電、光伏發(fā)電、交流負荷和直流負荷預(yù)測數(shù)據(jù)。實施例地區(qū)典型日風力發(fā)電預(yù)測曲 線、光伏發(fā)電預(yù)測曲線、交流負荷預(yù)測曲線和直流負荷預(yù)測曲線如圖3所示；
[0155] (2)統(tǒng)計交直流混合微網(wǎng)中可控微源運行特性，建立所有可控微源的燃料成本-機 端功率函數(shù)，得到可控微源的燃料成本-機端功率曲線如圖4所示；
[0156] (3) W-個小時為優(yōu)化時段，將交直流混合微網(wǎng)未來一天內(nèi)的優(yōu)化運行分為24個 時段，W交直流混合微網(wǎng)全天運行的燃料成本、運行維護成本、購電成本和環(huán)境成本最小為 目標函數(shù)，考慮交直流混合微網(wǎng)內(nèi)部的各時段功率平衡、發(fā)電單元功率限額、發(fā)電單元最短 啟停時間、蓄電池電量容量和功率容量約束、蓄電池容量連續(xù)性、并網(wǎng)聯(lián)絡(luò)線容量限額為約 束條件，建立交直流混合微網(wǎng)優(yōu)化運行模型；
[0157] (4)采用協(xié)同進化改進黑桐算法對步驟(3)中的交直流混合微網(wǎng)優(yōu)化運行的數(shù)學 模型進行求解，得到各時段各區(qū)域可控型微電源的出力和兩區(qū)域間交互功率的優(yōu)化方案如 圖5-8所示。
[0158] 在本實施例中，設(shè)定改進黑桐算法的最大迭代次數(shù)100,且連續(xù)10代進化停滯則認 為收斂;各區(qū)域?qū)?yōu)個體數(shù)50;協(xié)同進化框架中T取13。
[0159] 計算結(jié)果為總成本4483.84元。從圖5-8可知，交直流混合微網(wǎng)運行過程中，由于燃 料電池發(fā)電成本較低，優(yōu)先使用燃料電池發(fā)電，在燃料電池不能滿足供能需求時，再使用微 燃機發(fā)電；蓄電池在低電價時充電，在高電價時放電;微網(wǎng)在外網(wǎng)電價低時優(yōu)先使用外網(wǎng)購 電，外網(wǎng)電價高時主要使用內(nèi)部發(fā)電單元，降低了交直流混合微網(wǎng)的運行成本。
[0160]綜上所述，通過本實施例的測試結(jié)果，說明本發(fā)明提出的一種基于協(xié)同進化改進 黑桐算法的交直流混合微網(wǎng)優(yōu)化運行方法可有效實現(xiàn)交直流混合微網(wǎng)的優(yōu)化運行，充分發(fā) 揮微網(wǎng)的經(jīng)濟效益、環(huán)保效益。
【主權(quán)項】
1. 一種交直流混合微網(wǎng)的優(yōu)化運行方法，其特征在于，包括如下步驟： 1) 采集交直流混合微網(wǎng)的交流區(qū)域和直流區(qū)域的負荷信息數(shù)據(jù)、氣象信息數(shù)據(jù)，綜合 微網(wǎng)運行的歷史數(shù)據(jù)，對未來一天的數(shù)據(jù)進行預(yù)測，得到未來一天內(nèi)交直流混合微網(wǎng)的交 流區(qū)域負荷、直流區(qū)域負荷、風能和太陽能功率預(yù)測數(shù)據(jù)； 2) 統(tǒng)計交直流混合微網(wǎng)的微電源特性，建立微網(wǎng)內(nèi)所有可控型微電源的出力數(shù)學模 型； 3) 將交直流混合微網(wǎng)未來一天內(nèi)的經(jīng)濟運行分為24個時段，W交直流混合微網(wǎng)全天運 行燃料成本、運行維護成本、購電成本和環(huán)境成本最小為目標函數(shù)，考慮微網(wǎng)內(nèi)部的各時段 電能平衡約束、并網(wǎng)聯(lián)絡(luò)線容量約束、可控微源的出力限制、可控微源的最短啟停時間約 束、蓄電池的電量容量約束、蓄電池的功率容量約束、蓄電池的容量連續(xù)性約束和計算蓄電 池的周期始末電能平衡約束，建立交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型； 4) 基于步驟1)中的交流區(qū)域負荷、直流區(qū)域負荷、風能和太陽能功率預(yù)測數(shù)據(jù)，采用協(xié) 同進化改進黑桐優(yōu)化算法對步驟3)中的交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型進行求解，得到 各時段可控型微電源的出力及混合微網(wǎng)購電量優(yōu)化方案。2. 根據(jù)權(quán)利要求1所述的一種交直流混合微網(wǎng)的優(yōu)化運行方法，其特征在于，步驟2)所 述的出力數(shù)學模型包括： (1) 燃料電池的出力模型：其中，各符號的定義如下:ffuel,FC為燃料電池的燃料成本，Cgas,FC為燃料電池使用的燃料 價格，LHVgas,FC為燃料的低熱值，PfC為燃料電池的出力大小，郵為燃料電池的效率。 (2) 微型燃氣輪機的出力模型：其中，各符號的定義如下:ffuel,MT為微型燃氣輪機的燃料成本，Cgas,MT為微型燃氣輪機使 用的燃料價格，LHVgas,MT為微型燃氣輪機的燃料的低熱值，PmT為微型燃氣輪機的的出力大 小，riMT為微型燃氣輪機的的效率。3. 根據(jù)權(quán)利要求1所述的一種交直流混合微網(wǎng)的優(yōu)化運行方法，其特征在于，步驟3)所 述的的交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型包括目標函數(shù)和約束條件，其中，目標函數(shù)為：其中，各符號定義如下:F為目標函數(shù)，Cfuel為燃料成本，C?為運行維護成本，Cbuy為購電 成本，α為污染折算成本;t為時段，N為微源總個數(shù)，ft( ·)為第i個微源在t時段的燃料成本 函數(shù)，Pi(t)為第i個微源在t時段的實際輸出功率，ki為第i個微源的運行維護成本系數(shù)，at 為t時段時的購電電價，Pbuy(t)為混合微網(wǎng)在t時段時的購電量，Μ為污染物總數(shù)，&為第巧中 污染物的環(huán)境評價標準，為第i個機組的第巧巾污染物的排污系數(shù)； 約束條件為： (1) 發(fā)電單元功率限額： Pimin《Pi(t)《Pimax (2) 最短啟停時間約束：若在改進黑桐算法的迭代趨于收斂時，部分機組沒有滿足啟停時間的約束，對沒有滿 足啟停時間約束的機組按如下方法加 W調(diào)整：如果機組連續(xù)啟動時間過短，則強制延長運 行時間；如果停機時間過短，則將違規(guī)停機時間強制改為W最低功率運行，再調(diào)整其他機組 出力W滿足系統(tǒng)的其他約束； (3) 蓄電池電量容量和功率容量約束：(5) 計算蓄電池周期始末電能平衡約束： E(0)=E(24 化(t) (6) 全系統(tǒng)實時電能平衡約束： EPi(t)+Pbuy(t)+PES(t)=&(t)+PLoss(t) (7) 并網(wǎng)聯(lián)絡(luò)線容量約束：其中：Ρι"3Χ和Ρι"ιη分別為輸出上限和下限，Pi(t)為第i個微源在t時段的實際輸出功率； Τηη為可控型微電源的連續(xù)開啟時間，Toff為可控型微電源的連續(xù)停機時間，可控型微電源的 連續(xù)開啟和停機時間的下限分別為式嚴與石f ;EES(t)為t時段的蓄電池的電量容量，Emin、 Emax分別為蓄電池電量容量的下限和上限；PES(t)為t時段的蓄電池輸出量，W輸出為正方 向，PES,min、PES,max分別為蓄電池輸出量的下限和上限；；E(t)為蓄電池 t時段的電量，E(t-l) 為蓄電池 t-1時段的電量;riESl為蓄電池充電效率;T1ES2為蓄電池放電效率；ΛΤ為單位時段， 本發(fā)明中AT=l;Pbuy(t)為交直流混合微網(wǎng)t時段的購電功率，PD(t)為交直流混合微網(wǎng)t時 段的電負荷需求，扣。ss(t)為交直流混合微網(wǎng)t時段的有功網(wǎng)損;巧Γ為聯(lián)絡(luò)線容量約束下的 最大購電功率。4.根據(jù)權(quán)利要求1所述的一種交直流混合微網(wǎng)的優(yōu)化運行方法，其特征在于，步驟4)所 述的對交直流混合微網(wǎng)優(yōu)化運行數(shù)學模型進行求解包括： (1) 初始化:輸入微電源、負荷參數(shù);輸入算法參數(shù)包括星體范圍、最大尋優(yōu)次數(shù)、星體 的個數(shù)、新產(chǎn)生星體距離黑桐的最小距離和星體的維數(shù);各時段交、直流區(qū)域間的交互功 率;交、直流區(qū)域間的交互功率計算迭代次數(shù)PP = 1;設(shè)定交、直流區(qū)域間的交互功率改變量 ε； (2) 將交直流混合微網(wǎng)的各優(yōu)化變量根據(jù)位置劃分到交流區(qū)域和直流區(qū)域中； (3) 將交互功率計入各區(qū)域功率平衡約束中，結(jié)合步驟3)所述的約束條件，給出功率平 衡約束如下：其中，各符號定義如下:PAC(t)為t時段交流區(qū)域微電源的發(fā)電功率輸出總量，包括購電 功率;/皆(0為t時段交流區(qū)域蓄電池的輸出功率;巧^'(0為t時段交流區(qū)域負荷;巧己(0為t 時段交流區(qū)域網(wǎng)損;PM(t)為t時段交、直流區(qū)域的交互功率;巧1的為t時段交、直流區(qū)域間 的雙向換流器的損耗;U0為雙向換流器的損耗系數(shù)，當Pm( t)〉0時，U0 = 1，否則U0 = 0; Pdc (t) 為t時段直流區(qū)域微電源的發(fā)電功率輸出總量；P#'(〇為t時段直流區(qū)域蓄電池的輸出功 率;巧嚴(。為t時段直流區(qū)域負荷;巧.志的為t時段直流區(qū)域網(wǎng)損； (4) 設(shè)置改進黑桐算法的迭代次數(shù)iter = 1; (5) 在區(qū)域P中隨機產(chǎn)生N個星體； (6) 計算初始星體的自適應(yīng)度值，選擇最優(yōu)的自適應(yīng)度值確定初始黑桐； (7) 更新星體的位置； (8) 計算更新后星體的自適應(yīng)度值，若星體的自適應(yīng)度值優(yōu)于初始黑桐的自適應(yīng)度值， 交換所述星體與初始黑桐的位置，所述星體為新的黑桐，其他星體向著新的黑桐移動更新， 否則繼續(xù)靠近初始黑桐位置更新星體的位置，星體向著黑桐移動更新的速度和慣性為：其中，各符號定義如下:vit為星體t時刻的速度，v^i為星體t-1時刻的速度，定義t = l 時速度為0;xip為第i個星體的歷史最優(yōu)適應(yīng)值;xi%第i個星體在t時刻的位置;祉Η為黑桐 的位置;ri和η為[0，1]之間的隨機數(shù)；ω為星體向著黑桐移動的慣性參量，COmin、《max為星 體向著黑桐移動的慣性參量的下限和上限，取Wmax = 0.9, C〇min = 0.4;Niter為當前改進黑桐 算法的迭代次數(shù)，Nitermax為改進黑桐算法的最大迭代次數(shù)； 星體向著黑桐移動的速度約束為：其中，各符號定義如下:Vid為星體i向著黑桐移動的速度，Vidmin與Vidmax分別為星體向著 黑桐移動的最小與最大速度；Xidmax、Xidmin分別為星體i向著黑桐移動一次的最小、最大距 離;定值參數(shù)L根據(jù)算例取定，用來限制速度的幅度； (9) 判斷星體是否符合吞隧條件，若符合，黑桐將星體吞隧，并在可行空間內(nèi)隨機位置 重新生成新星體，在可行空間內(nèi)隨機位置重新生成新星體，具體是采用Logistic映射產(chǎn)生 混濁變量，映射公式如下： χ??+ι=μχ??(1-χι<) 其中，各符號定義如下:Xk、Xk+l為星體在第k、k+l次尋優(yōu)中的位置;μ為常數(shù)，取4令映射 進入混濁狀態(tài);所述映射有不動點0、〇. 25、0.50、0.75、1，若星體的初始位置為所述的不動 點，則選用普通隨機數(shù)生成新星體的位置； 若不符合吞隧條件，星體繼續(xù)向靠近黑桐位置移動更新星體的位置； (10) 判斷當前的迭代次數(shù)是否達到最大迭代次數(shù)，若達到則結(jié)束計算，輸出計算結(jié)果， 進入步驟(11);若不滿足則設(shè)定改進黑桐算法的迭代次數(shù)iter = iter+l，返回第(7)步； (11) 交、直流子區(qū)域分別產(chǎn)生若干個最優(yōu)個體； (12) 將交、直流子區(qū)提交的個體互相組合，得到多種全網(wǎng)機組出力計劃，由于可能存在 通過跨子區(qū)的合作互補，故逐時段進行如下優(yōu)化： rnin ( Cfuel (t) +C〇m( t ) ~l~Cbuy (t ) +Cw( t)) 其中，各符號定義如下:Cfuel(t)為t時段的燃料成本，Ccm(t)為t時段的運行維護成本， Cbuy(t)為t時段的購電成本，Cw(t)為t時段的污染折算成本，約束條件與步驟3)-致； (13) 將所有24個時段都進行第(12)步所述的優(yōu)化處理，得到全周期總成本； (14) 在第（13)步中得到的全周期總成本中找出綜合成本最低的組合，計算各時段交、 直流區(qū)域間的交互功率，當各時段交、直流區(qū)域間的交互功率的計算迭代次數(shù)PP=1時，若 各時段交、直流區(qū)域間的交互功率相對于初始的各時段交、直流區(qū)域間的交互功率改變量 均值小于給定的ε，則轉(zhuǎn)第（15)步；當各時段交、直流區(qū)域間的交互功率的計算迭代次數(shù)PP >1時，若各時段交、直流區(qū)域間的交互功率相對于上一次的若各時段交、直流區(qū)域間的交 互功率改變量均值小于給定的ε，則轉(zhuǎn)第(15)步； 否則將綜合成本最低的組合作為新的邊界傳輸功率值，賦給交流、直流區(qū)域，計算各時 段交直流子區(qū)間交互功率，公式如下： Ρ、, (η 二 Ρ,八η-聽；（!)-咕（0 - .蠟J0 其中，各符號定義如下:PM(t)為t時段交、直流區(qū)域的交互功率;PAC(t)為t時段交流區(qū) 域微電源的發(fā)電功率輸出總量，包括購電功率；/3,.;.<(〇為t時段交流區(qū)域蓄電池的輸出功 率；f/·的為t時段交流區(qū)域負荷；巧志的為t時段交流區(qū)域網(wǎng)損;的為t時段交、直流區(qū) 域間的雙向換流器的損耗;U0為雙向換流器的損耗系數(shù)，當PM(t)〉0時，110 = 1，否則u〇 = 0; 交、直流區(qū)域間的交互功率計算迭代次數(shù)PP=PP+1，轉(zhuǎn)第(3)步繼續(xù)計算； (15)停止計算，輸出結(jié)果。
【文檔編號】H02J5/00GK106100002SQ201610604292
【公開日】2016年11月9日
【申請日】2016年7月28日 公開號201610604292.9, CN 106100002 A, CN 106100002A, CN 201610604292, CN-A-106100002, CN106100002 A, CN106100002A, CN201610604292, CN201610604292.9
【發(fā)明人】李鵬, 狄開麗, 李鑫明, 陳安偉, 周金輝, 趙波
【申請人】華北電力大學（保定）, 國網(wǎng)浙江省電力公司電力科學研究院