莊慧敏, 肖 建
(1. 成都信息工程大學控制工程學院,四川 成都610225;2. 西南交通大學電氣工程學院,四川 成都610031)
分布式發(fā)電(distributed generation,DG)技術,尤其是可再生能源(renewable energy sources,RES)的發(fā)展可以優(yōu)化能源結構,提升綠色能源利用率,推動節(jié)能減排和實現(xiàn)經濟可持續(xù)發(fā)展[1]. 近年來,配電網(wǎng)中DG 的滲透率迅速增長,傳統(tǒng)配電網(wǎng)正逐步演變?yōu)榫哂斜姸嗫烧{可控資源(可調DG、儲能裝置、柔性負荷等)的主動配電網(wǎng).主動配電網(wǎng)的優(yōu)化調度策略是主動配電網(wǎng)對可控資源實施主動管理的核心技術,也是實現(xiàn)網(wǎng)絡安全經濟運行的重要手段[2].
目前,對含風力發(fā)電的電力系統(tǒng)調度及微網(wǎng)優(yōu)化調度的研究[3-4]比較多,而對主動配電網(wǎng)的優(yōu)化調度研究較少.
現(xiàn)有的主動配電網(wǎng)優(yōu)化調度模型按目標函數(shù)可以分為配電網(wǎng)運行成本最?。?-6]、電力公司效益最大[7-8]、社會綜合效益最大[9].文獻[6]以整個調度周期內配電網(wǎng)運行成本最低為目標函數(shù),以可控分布式能源及聯(lián)絡開關作為控制手段,實現(xiàn)主動配電網(wǎng)的優(yōu)化調度. 在文獻[6]的目標函數(shù)中,將聯(lián)絡開關的狀態(tài)作為控制變量,達到優(yōu)化網(wǎng)絡結構的目的,但沒有給出發(fā)電功率關于該狀態(tài)量的具體表達式,在實施方面存在一定的困難. 文獻[7]在主動配電網(wǎng)的調度模型中考慮了輸電和高壓配電網(wǎng)網(wǎng)損的影響,以電力公司效益最大化為目標,得到了區(qū)域電網(wǎng)中含分布式電源的主動配電網(wǎng)優(yōu)化調度模型,并得出效益最大時分布式電源對主動型配電網(wǎng)的最佳滲透率.
在調度對象方面,大多以DG 為主體,部分文獻包括了儲能裝置,少數(shù)文獻[10]以可控負荷為主體.目前考慮需求響應的優(yōu)化調度策略成為研究熱點[11-12].文獻[12]針對微網(wǎng)的分布式發(fā)電調度,提出了考慮實時電價、可中斷負荷等需求側響應策略的分布式發(fā)電協(xié)調優(yōu)化調度模型.
通過對主動配電網(wǎng)優(yōu)化調度的研究現(xiàn)狀可以看出,目前的調度策略存在以下問題:
(1)多數(shù)研究只考慮了日前調度策略,沒有考慮RES 發(fā)電預測及負荷預測的偏差對調度策略的影響;
(2)主要集中于DG 的有功功率優(yōu)化,而沒有考慮DG 作為無功資源對電網(wǎng)的輔助服務作用;
(3)沒有考慮DG 在配電網(wǎng)中引起的過壓問題,過壓可能會觸發(fā)保護裝置將其從電網(wǎng)切除,浪費電網(wǎng)資源和可再生能源.
基于上述問題,本文提出主動配電網(wǎng)兩階段優(yōu)化調度策略:日前優(yōu)化調度和日內滾動修正兩階段調度.日前經濟調度模型以一個完整調度周期內的配電網(wǎng)運行成本最低為目標函數(shù),以可調度DG、可控負荷作為控制手段,同時確保儲能系統(tǒng)在整個調度周期的能量守恒及容量約束;而日內滾動修正以日前調度計劃為基礎,以DG 的有功消減量和無功功率作為控制變量,每小時根據(jù)最新的短期預測信息更新日前調度值,從而優(yōu)化網(wǎng)絡電壓,降低網(wǎng)絡損耗,提高DG 并網(wǎng)能力. 最后,通過擴展的IEEE34 節(jié)點配電網(wǎng)算例對本文建立的模型進行了驗算.
主動配電網(wǎng)的日前調度模型是對次日可調分布式能源(發(fā)電機、儲能系統(tǒng))有功功率及響應負荷減載量、減載時段的計劃.因此,主動配電網(wǎng)優(yōu)化調度的控制向量可以表示為
控制向量的前Ng個分量是可控分布式發(fā)電單元的輸出功率,中間Ns個分量是儲能單元的充放電功率,最后的Nl個分量是響應負荷的減載量.本文選取可減載負荷模型作為需求響應模型.主動配電網(wǎng)可在負荷高峰時段與用戶簽訂負荷減載合同,并給予補償.
將目標函數(shù)F 定義為整個調度周期內主動配電網(wǎng)的運行成本最小.運行成本包括分布式發(fā)電機的發(fā)電成本、從關口購電的成本以及補償需求響應的成本.
式中:t 為時段編號;
T 為調度周期的總時段數(shù),對于每個階段可認為各DG 單元功率、儲能單元功率以及負荷大小不變;
Δt 為單位階段的時長;
Ng為可調DG 單元總數(shù);
g 為可調DG 編號;
Nl為響應負荷總數(shù);
l 為可控負荷編號;
Cg(t)和PDGg(t)分別為第g 個可調DG 單元在時段t 的單位發(fā)電成本和有功功率;
Cl(t)和PRSl(t)分別為第l 個響應負荷在時段t 的減載補償價和減載量;
PEX(t)和CEX(t)分別為配變關口的交換功率及其電價.
該目標函數(shù)從本質上能反映出優(yōu)化調度分布式能源、儲能充放電策略及需求響應策略合理化帶來的收益.一般來說,分布式能源發(fā)電成本優(yōu)于電網(wǎng),在負荷一定的情況下,分布式發(fā)電利用率越高,電網(wǎng)輸送電能就越少,目標函數(shù)的值就越小.此外,儲能系統(tǒng)或需求響應的優(yōu)化調度可以實現(xiàn)削峰填谷帶來的峰谷電價差收益,并且響應用戶還會因此獲取一定的補償,也可以使目標函數(shù)值變小.
主動配電網(wǎng)優(yōu)化調度的約束條件除了常規(guī)的潮流平衡方程、節(jié)點電壓上下限約束、支路功率最大限制約束之外,還包括如下約束:
(1)功率平衡
式中:N 為所有DG 單元個數(shù);
PL(t)為時段t 的總負荷功率.
(2)儲能單元能量守恒
式中:
Es0和Es1為第0 和1 個儲能單元在時段t 的剩余能量.
(3)發(fā)電單元及儲能單元的功率和能量限制
式中:
Pmax,DGj為第j 個分布式發(fā)電單元的功率上限值(對于風力發(fā)電或光伏發(fā)電為預測值);
Pmin,ESs(小于0)、Pmax,ESs、Emin,s、Emax,s分別為第s 個儲能單元充放電功率的下限值和上限值、剩余能量的最小允許值和最大允許值.
(4)需求響應負荷的減載量限制
式中:
Pmin,RSl、Pmax,RSl分別為第l 個需求響應負荷的減載量下限值和上限值.
(5)配變關口功率約束
為了減小主動配電網(wǎng)的功率波動對輸電網(wǎng)的影響,需考慮配電網(wǎng)根節(jié)點的關口交換功率約束,即
本文建立的優(yōu)化調度模型與傳統(tǒng)模型的區(qū)別主要在于:將儲能單元及需求響應用于抑制可再生能源功率的波動,提高系統(tǒng)運行的經濟性. 所建立的模型是復雜的混合整數(shù)規(guī)劃問題.
由于RES 發(fā)電預測及負荷預測都受其自身不確定性的影響,因此在配電網(wǎng)的實際運行中,可能會出現(xiàn)節(jié)點電壓越限的情況. 此外,當配電網(wǎng)出現(xiàn)意外事故(例如某支路斷線)時還可能導致某支路過流等異常現(xiàn)象.因此,為了保證配電網(wǎng)安全、可靠和經濟地運行,還需在日前優(yōu)化調度的基礎上進行日內優(yōu)化調度.
日內優(yōu)化調度模型是根據(jù)最新的短期預測信息和負荷需求等,對日前調度計劃進行滾動調整.在滾動修正階段,由于更接近實際情況,并且合理地消減可控DG 的發(fā)電量、鼓勵DG 向配電網(wǎng)提供無功支持,可改善電壓問題、提高大規(guī)模DG 的消納能力[13].因此,引入可控DG 有功功率消減量及向配電網(wǎng)提供的無功功率,既考慮了系統(tǒng)運行的經濟性,又從數(shù)學角度把消減量和無功輸出作為松弛變量,保證模型存在可行解.
日內優(yōu)化調度的啟動周期一般為30 ~60 min,本文取調度周期為1 h,即每小時執(zhí)行一次.調度過程為:收集配電網(wǎng)系統(tǒng)信息(網(wǎng)絡狀態(tài)信息、可再生能源短期平均功率預測值和負荷需求預測值),然后計算三相潮流,若發(fā)現(xiàn)有技術條件約束不能滿足(例如節(jié)點電壓越限等),則啟動日內優(yōu)化調度程序,修正參與網(wǎng)絡電壓控制的分布式發(fā)電單元的有功功率和無功功率輸出.
在DG 穩(wěn)態(tài)運行時采用PQ 類型. DG 通過電力電子裝置或常規(guī)旋轉電機接口并網(wǎng),且其并網(wǎng)功率已能實現(xiàn)有功無功分別獨立調節(jié)[14].
考慮到有載調壓變壓器(on-load tap changers,OLTC)分接頭不應頻繁調整,本文假設在日內優(yōu)化的實時控制時段內,OLTC 分接頭固定. 當配電網(wǎng)運行過程中出現(xiàn)過壓或支路功率越限等異常情況時,可通過消減可控DG 的有功功率和控制其無功輸出進行緩解或消除異常情況.
日內優(yōu)化調度具有實時性,因此模型的控制變量不宜過多,目標主要偏向于系統(tǒng)的經濟性. 在日內優(yōu)化調度中,配電網(wǎng)將對支持優(yōu)化調度的發(fā)電單元給予一定的補償.為此,目標函數(shù)J 考慮主動配電網(wǎng)的調度成本(包括DG 消減有功功率的成本CGC和提供無功支持的成本CQ)及能量損耗成本CLOS,即
(1)消減有功功率成本
式中:
NGC為可消減有功功率的DG 個數(shù);
CGCj為第j 個DG 消減單位電量(1 kW·h)的成本;
Pset,DGj為可調DG 的日前優(yōu)化調度值,對于不可調DG(例如風力或光伏發(fā)電單元)為發(fā)電功率預測值;
PDGj為第j 個DG 的控制目標值;
PGCj為第j 個DG 的有功功率消減值.
(2)提供無功支持成本
式中:
NQ為可向配電網(wǎng)提供無功支持的DG 個數(shù);
Qk和CQk分別為第k 個DG 輸出的無功功率和提供單位無功量(1 kVar·h)的成本.
(3)能量損耗成本
式中:CLOS為單位能量損耗成本;
PLOSi為第i 條支路損失的有功功率;
Nbr為配電網(wǎng)支路總數(shù);
Gi(m,n)為連接節(jié)點m 和n 的支路i 的電導;V 和θ 分別為節(jié)點的電壓幅值和相位.
消減有功功率是一種比較昂貴的控制手段,通常只在電網(wǎng)安全受到威脅時采用. 因此,其控制成本取較大值.
約束條件除了日前優(yōu)化調度模型的約束外,還包括如下約束:
式中:
Pmax,GCj為第j 個DG 的最大允許消減量;
Qmax,k、Qmin,k分別為第k 個DG 的最大、最小無功率.
本文提出的兩階段調度優(yōu)化模型是一個混合整數(shù)規(guī)劃模型,可采用優(yōu)化規(guī)劃軟件CPLEX 對其求解.考慮到本文的日內調度中需要計算配電網(wǎng)潮流,用Matlab 編寫了線性規(guī)劃模型的混合整數(shù)規(guī)劃優(yōu)化模型求解程序,而配電網(wǎng)模型建立及日內調度的潮流計算則由軟件PSCAD/EMTDC 執(zhí)行,兩個軟件的接口程序采用C+ +語言編寫.
以擴展的IEEE34 節(jié)點配電網(wǎng)為算例進行分析.測試網(wǎng)絡結構如圖1 所示.
圖1 IEEE34 節(jié)點測試系統(tǒng)Fig.1 IEEE34-node test system
該系統(tǒng)包含33 條支路,輻射運行,電壓等級為24.9 kV,最大負荷總功率為7 500 kW,可中斷負荷為1 000 kW.配電網(wǎng)節(jié)點802、812、818 各連接一臺燃氣輪機(gas turbine,GT),節(jié)點856 與一臺柴油(diesel,DE)發(fā)電機相連,相關的發(fā)電功率極限和發(fā)電成本單價見表1. 兩個不可調的光伏(photovoltaic,PV)發(fā)電單元分別和配電網(wǎng)節(jié)點844、826 相連,峰值容量均為500 kW,節(jié)點842、824 各連接一個儲能裝置(energy storage source,ESS)(ESS842、ESS824),充放電功率上限均為120 kW,電量上限均為600 kW·h,充放電效率均為0.95.其他網(wǎng)絡和負荷參數(shù)見文獻[15].
表1 可調DERs 的模型參數(shù)Tab.1 Model parameters of the dispatchable DERs
該算例以一天24 h 為一個完整調度周期,調度間隔為1 h,配變關口的市場電價采用階梯電價(如圖 2 所示),可中斷負荷的補償價為0.9 元/(kW·h).
圖2 市場電價Fig.2 Market electrical price
為了分析簡便,兩個PV 單元采用相同的光照曲線,全天的有功功率預測見圖3. 圖3 中的功率值為相對于最大功率的標幺值.
用程序計算出測試系統(tǒng)各調度單元全天的最優(yōu)有功功率如圖4 所示.
由圖4 看出,DG 的總功率特性與負荷特性基本保持一致,當負荷較大時DG 的總功率較大,負荷較小時總功率較小.而從配電網(wǎng)關口的購電量則與負荷成反比趨勢,即在荷峰時段從關口購電量最小,在荷谷時段從關口購電量較大. 此外,GT、DE的功率特性和可中斷負荷的響應特性主要是由根節(jié)點關口的實時電價決定的,在峰時期間電價成本較高,GT、DE 基本都處于滿發(fā)狀態(tài);在荷谷時段電價成本較低,GT、DE 功率為其最小值;在最高峰期間,電價成本高于可中斷負荷的補償價,可中斷負荷供電.
圖3 光伏發(fā)電的有功功率預測Fig.3 Forecasted active power outputs of photovoltaic generation
圖4 日前調度結果Fig.4 Day-ahead scheduling solution
通過圖4 還發(fā)現(xiàn),儲能單元的充放電策略通過優(yōu)化計算也得到了合理的分配:在負荷最低谷期間(1:00 ~5:00)處于充電狀態(tài),在負荷最高峰期間進行(9:00 ~14:00)放電,在其他時間(負荷平穩(wěn)期間)充放電功率為0.
因此,通過對主動配電網(wǎng)各可控單元的優(yōu)化調度,不僅能夠減少配電網(wǎng)的運行成本,還可對大電網(wǎng)起到積極的削峰填谷作用.
取系數(shù)
α=1.5 元/(kW·h),
β=γ=0.7 元/(kW·h).
GT 和DE 的初始功率因數(shù)素均為0.85. 為驗證日內優(yōu)化調度策略的有效性,以負荷的低谷(2:00 ~3:00)與高峰(10:00 ~11:00)兩個典型時段為例.
經優(yōu)化計算得到這兩個時段的GT 和DE 的有功消減量及無功輸出,見表2 和表3. 荷谷時段的網(wǎng)絡損耗由日內優(yōu)化前負荷功率的3. 32%降為3.11%,而荷峰時段由3.24%降為3.02%.此時全天運行成本為40 926.19 元,較未采用優(yōu)化調度時的運行成本42 309. 72 元,運行成本降低了3.27%.
表2 荷谷時段日內優(yōu)化結果Tab.2 Intra-day solution in load trough period
表3 荷峰時段日內優(yōu)化結果Tab.3 Intra-day solution in load peak period
由表2 可看出,荷谷時段的GT1、GT2 的無功輸出降為零,GT3 和DE 的無功輸出為負,DE 的有功輸出消減了120 kW.這是因為在荷谷時段,母線節(jié)點電壓較高,而有功消減成本比提供無功支持的成本高得多,因此,主要借助無功調節(jié)來調壓:減少無功功率輸出(例如GT1、GT2)或發(fā)出容性無功功率(例如GT3 和DE). 由表3 可看出,荷峰時段各GT 的功率因素調整為0.8,增加了無功輸出,從而提升母線節(jié)點電壓.
日內調度調整前后的DG 母線和主要負荷母線A 相的相電壓幅值曲線如圖5、圖6 所示.
圖5 荷谷時段各母線A 相電壓曲線Fig.5 Phase-A voltage profiles at various network buses in load trough period
圖6 荷峰時段各母線A 相電壓曲線Fig.6 Phase-A voltage profiles at various network buses in load peak period
由圖5 ~6 可看出,通過日內優(yōu)化調整日前調度值后,無論是荷谷時期還是荷峰時段,母線節(jié)點電壓都由原來的越限狀態(tài)恢復到電壓規(guī)定范圍內(±5%以內),這樣可避免引起電壓越限的DG 被過壓(或低壓)保護裝置切除.
本文針對主動配電網(wǎng)的技術特點及其現(xiàn)有優(yōu)化調度方法存在的問題,提出了兩階段一體化調度方法:日前經濟優(yōu)化調度和日內滾動修正. 日前優(yōu)化調度模型以分布式能源、儲能單元及可控負荷為調度對象,以整個調度周期內主動配電網(wǎng)的運行成本最小為目標,使配電網(wǎng)的長期運行更經濟.
為了避免高滲透率的DG 引起過壓等問題而被切除,在日內實時優(yōu)化調度模型中,將可調DG的無功功率也作為優(yōu)化變量,參與配電網(wǎng)的電壓調整,實現(xiàn)可調DG 的有功功率輸出和無功功率輸出的全局最優(yōu)分配,提高配電網(wǎng)消納新能源的能力.最后通過算例結果分析,驗證了所提主動配電網(wǎng)優(yōu)化調度模型的正確性和有效性.
本文所研究的主動配電網(wǎng)優(yōu)化調度模型是建立在市場電價和可控負荷補償價格為已知的基礎上的,有待進一步將市場電價和需求響應策略融入到優(yōu)化調度模型中,使主動配電網(wǎng)獲得更好的優(yōu)化經濟運行效果.
致謝:本文工作得到成都信息工程大學引進人才項目(KYTZ201317)的資助.
[1] 王成山,李鵬. 分布式發(fā)電、微網(wǎng)與智能配電網(wǎng)的發(fā)展與挑戰(zhàn)[J]. 電力系統(tǒng)自動化,2011,34(2):10-14.WANG Chengshan, LI Peng. Development and challenges of distributed generation,micro-grid and smart distribution network[J]. Automation of Electric Power Systems,2011,34(2):10-14.
[2] 范明天,張祖平,蘇傲雪,等. 主動配電系統(tǒng)可行技術的研究[J]. 中國電機工程學報,2013,33(22):12-18.FAN Mingtian,ZHANG Zuping,SU Aoxue,et al.Enablingtechnologiesforactivedistribution systems[J]. Proceedings of the CSEE,2013,33(22):12-18.
[3] 陳道君,龔慶武,張茂林,等. 考慮能源環(huán)境效益的含風電場多目標優(yōu)化調度[J]. 中國電機工程學報,2011,31(13):10-17.CHEN Daojun,GONG Qingwu,ZHANG Maolin,et al.Multi-objective optimal dispatch in wind power integrated system incorporating energy-environmental efficiency[J]. Proceedings of CSEE,2011,31(13):10-17.
[4] 洪博文,郭力,王成山. 微電網(wǎng)多目標動態(tài)優(yōu)化調度模型與方法[J]. 電力自動化設備,2013,33(3):100-107.HONG Bowen,GUO Li,WANG Chengshan,et al.Model and method of dynamic multi-objective optimal dispatch for microgrid[J]. Electric Power Automation Equipment,2013,33(3):100-107.
[5] WIERZBOWSKI M, OLEK B. Integration of the embedded generation into distribution systems at the competitive markets[C]∥Innovative Smart Grid Technologies Conference (ISGT 2014). Washington D.C.:IEEE,2014:1-5.
[6] 尤毅,劉東,鐘清,等. 主動配電網(wǎng)優(yōu)化調度策略研究[J]. 電力系統(tǒng)自動化,2014,38(9):177-183.YOU Yi,LIU Dong,WU Qing,et al. Research on optimal schedule strategy for active distribution network[J]. Automation of Electric Power Systems,2014,38(9):177-183.
[7] 竇震海,牛煥娜,高燕,等. 主動型配電網(wǎng)日前調度策略研究[J]. 農業(yè)工程學報,2014,30(11):126-133.DOU Zhenhai,NIU Huanna,GAO Yan,et al. Study on day-ahead dispatch strategy of active distribution network[J]. Transactions of the Chinese Society of Agricultural Engineering,2014,30(11):126-133.
[8] GILL S,KOCKAR I,AULT G W. Dynamic optimal power flow for active distribution networks[J]. IEEE Transactions on Power Systems,2014,29(1):121-131.
[9] AYMAN B,SALAMA M M A. Management scheme for increasing the connectivity of small-scale renewable DG[J]. IEEE Transactions on Sustainable Energy,2014,27(1):1-8.
[10] BARTONA J,HUANGB S,DAVID I. The evolution of electricity demand and the role for demand side participation in ADN[J]. Energy Policy,2013,52(5):85-102.
[11] 劉小聰,王蓓蓓,李揚. 智能電網(wǎng)下計及用戶側互動的發(fā)電日前調度計劃模型[J]. 中國電機工程學報,2013,33(1):30-38.LIU Xiaocong,WANG Beibei,LI Yang. Day-ahead generation scheduling model considering demand side interaction under smart grid paradigm[J]. Proceedings of the CSEE,2013,33(1):30-38.
[12] 劉自發(fā),劉剛,劉幸. 基于量子差分進化算法的分布式電源協(xié)調優(yōu)化調度[J]. 電網(wǎng)技術,2013,37(7):1922-1928.LIU Zifa,LIU Gang,LIU Xing. Coordinated optimal dispatching of distributed generation based on quantum differential evolution algorithm[J]. Power System Technology,2013,37(7):1922-1928.
[13] 王彩霞,魯宗相. 風電功率預測信息在日前機組組合中的應用[J].電力系統(tǒng)自動化,2011,35(7):13-18.WANG Caixia,LU Zongxiang. Unit commitment based on wind power forecast[J]. Automation of Electric Power Systems,2011,35(7):13-18.
[14] 劉楊華,吳政球,涂有慶,等. 分布式發(fā)電及其并網(wǎng)技術綜述[J]. 電網(wǎng)技術,2008,32(15):71-76.LIU Yanghua,WU Zhengqiu,TU Youqing,et al. A Survey on distributed generation and its networking technology[J]. Power System Technology, 2008,32(15):71-76.
[15] 耿立卓. 智能配電網(wǎng)狀態(tài)估計與測量配置[D]. 天津:天津大學,2011.