茍凱杰，呂鳴陽，高悅，陳衡，張國強(qiáng)，雷兢

(華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，北京 102206)

隨著世界能源緊缺、環(huán)境污染等問題的日益突出，積極響應(yīng)“雙碳”戰(zhàn)略目標(biāo)[1]，建設(shè)一個(gè)節(jié)約能源、減少污染、綠色低碳、可持續(xù)發(fā)展的智能化社會(huì)，大力發(fā)展可再生能源，規(guī)劃建設(shè)風(fēng)力發(fā)電、光伏發(fā)電項(xiàng)目，形成綠色清潔的發(fā)電方式成為研究熱點(diǎn)[2]。虛擬電廠(virtual power plant，VPP)概念的提出為達(dá)成這些目標(biāo)給出了嶄新的想法[3-5]。

目前，國內(nèi)外學(xué)者已經(jīng)開展了許多關(guān)于VPP運(yùn)行結(jié)構(gòu)和優(yōu)化方面的研究。文獻(xiàn)[6]對(duì)VPP在包含風(fēng)電、光伏等的情況下，分別考慮可中斷負(fù)荷和備用的協(xié)同調(diào)度模型采用經(jīng)典場景集成的算法進(jìn)行了研究和分析。文獻(xiàn)[7]在考慮需求響應(yīng)和用戶用電滿意度的情況下，搭建了優(yōu)化調(diào)度模型。文獻(xiàn)[8]在允許棄風(fēng)棄光的情況下，通過改進(jìn)量子粒子群算法，構(gòu)建了一個(gè)優(yōu)化VPP收益的模型。文獻(xiàn)[9]分析了VPP內(nèi)部不同組成成分的特性，得到了包含需求響應(yīng)的VPP內(nèi)部優(yōu)化調(diào)度模型。文獻(xiàn)[10]構(gòu)建了包含風(fēng)電、光伏、水電、儲(chǔ)能和生物質(zhì)等分布式能源的VPP調(diào)度模型，通過對(duì)該模型的研究，將多種分布式能源進(jìn)行組合聯(lián)用可以消除分布式能源本身不確定性的影響。文獻(xiàn)[11]定義了商業(yè)型和技術(shù)型VPP，以用戶側(cè)收益最大化為目標(biāo)，綜合利用激勵(lì)型和價(jià)格型需求響應(yīng)的調(diào)度作用，對(duì)比不同場景下VPP的經(jīng)濟(jì)效益。文獻(xiàn)[12]將大規(guī)模儲(chǔ)能技術(shù)運(yùn)用于VPP中，得到了可以有效提升VPP經(jīng)濟(jì)性的方法。文獻(xiàn)[13]利用算例仿真，考慮風(fēng)電輸出的不確定性因素，成功降低了由于風(fēng)電預(yù)測誤差而導(dǎo)致的系統(tǒng)運(yùn)行成本增加的影響。這些文獻(xiàn)均對(duì)VPP的調(diào)度優(yōu)化進(jìn)行了研究分析，但是大部分仍然側(cè)重于研究風(fēng)電的影響，綜合考慮風(fēng)-光-儲(chǔ)、需求響應(yīng)和電價(jià)變化的多重影響比較少。本文將這些影響因素分別組合進(jìn)行計(jì)算分析，經(jīng)過多種情況模型的比較相對(duì)優(yōu)良的VPP調(diào)度方案。

本文旨在對(duì)VPP日前經(jīng)濟(jì)調(diào)度進(jìn)行優(yōu)化，將不同情景運(yùn)行成本作為目標(biāo)函數(shù)，同時(shí)考慮功率平衡、設(shè)備運(yùn)行等約束條件。進(jìn)一步研究VPP在可以和主網(wǎng)進(jìn)行功率交換背景下的調(diào)度方法，采用基于反向?qū)W習(xí)的混沌映射自適應(yīng)粒子群算法，研究不同情景下分布式能源的出力。最終，經(jīng)過對(duì)實(shí)際案例的分析，得出不同情況下VPP的經(jīng)濟(jì)優(yōu)化和最小峰谷差運(yùn)行方法。

1 計(jì)及風(fēng)-光-儲(chǔ)和需求響應(yīng)的日前經(jīng)濟(jì)調(diào)度模型

1.1 風(fēng)力發(fā)電模型

風(fēng)力發(fā)電機(jī)的出力狀況與風(fēng)速變化緊密相關(guān)，風(fēng)力機(jī)發(fā)電的出力[14]

(1)

式中：PE為風(fēng)電機(jī)組的額定功率；νin為切入速度；νout為切出速度；νe為額定速度；ν為風(fēng)速，服從雙參數(shù)威布爾分布。

1.2 光伏發(fā)電模型

光伏發(fā)電系統(tǒng)采用光伏電池板進(jìn)行電能轉(zhuǎn)換[15]，其輸出功率受到日照強(qiáng)度和氣候條件的影響，可以表示為

Pvi=riηA.

(2)

式中：ri為i時(shí)段的光照強(qiáng)度，服從Beta分布；A為太陽能電池陣列板的總面積；η為光電轉(zhuǎn)化效率。

1.3 需求側(cè)響應(yīng)

a)可再生能源消納率?？稍偕茉聪{率Mre是在可再生能源發(fā)電量最高的時(shí)候?qū)⒁恍┴?fù)載轉(zhuǎn)移，以便更好地消耗可再生能源，可表示為[16]

(3)

式中：Pgr為可再生能源系統(tǒng)輸出的有功功率；Pg為所有電源提供的有功功率。

b)用戶舒適度。采取需求側(cè)響應(yīng)措施會(huì)使整體的負(fù)荷結(jié)構(gòu)發(fā)生變化，從而影響用戶的舒適度，其可以表示為[17]：

(4)

QGi=|Pwimax-Pwi|+|Pvimax-Pvi|+|Pzimax-Pzi|，

(5)

QLi=Pwimax+Pvimax+Pzimax.

(6)

式(4)—(6)中：Ccom為用戶的用戶舒適度；QGi為i時(shí)段采用需求側(cè)響應(yīng)后用電量與采用需求側(cè)響應(yīng)前用電量之差的絕對(duì)值的和；QLi為i時(shí)段采用需求側(cè)響應(yīng)前的總用電量；Pzi為i時(shí)段與主網(wǎng)交換有功功率；Pzimax為i時(shí)段采用需求側(cè)響應(yīng)前與主網(wǎng)交換有功功率；Pwi為i時(shí)段風(fēng)力機(jī)有功功率；Pwimax為i時(shí)段采用需求側(cè)響應(yīng)前風(fēng)力機(jī)有功功率；Pvi為i時(shí)段光伏有功功率；Pvimax為i時(shí)段采用需求側(cè)響應(yīng)前光伏有功功率。

c)電價(jià)的影響?？蛻粲秒娏颗c電價(jià)之間的關(guān)系用電價(jià)伸縮系數(shù)ξij表示。當(dāng)前時(shí)段電價(jià)和其他時(shí)段電價(jià)都將影響當(dāng)前時(shí)段用電量[18]，因此在分析價(jià)格伸縮性時(shí)，需要綜合考慮上述兩種影響因素。ξij可以定義為客戶用電量波動(dòng)率與價(jià)格波動(dòng)率的比值[19]：

(7)

式中：Qi、Cj分別為i(j)時(shí)段初始的用電量、電價(jià)；ΔQi、ΔCj分別為i(j)時(shí)段考慮電價(jià)影響后用電量、電價(jià)的改變量。i=j時(shí)為i時(shí)段的自伸縮系數(shù)ξii，i≠j時(shí)為i時(shí)段對(duì)j時(shí)段的互伸縮系數(shù)ξij。

由自伸縮系數(shù)和互伸縮系數(shù)可以組成需求側(cè)響應(yīng)電價(jià)伸縮系數(shù)矩陣ζ：

(8)

式中：下標(biāo)f、p、g分別對(duì)應(yīng)伸縮系數(shù)歸屬的峰值、平段和峰谷時(shí)段。采用分時(shí)電價(jià)后，用電變化量[20]

(9)

綜上，采用需求側(cè)響應(yīng)之后i時(shí)段的電量

Q′i=Qi+ΔQ.

(10)

2 目標(biāo)函數(shù)和約束條件

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

以系統(tǒng)最終的運(yùn)行成本最小作為研究目標(biāo)，針對(duì)不同方案的機(jī)組運(yùn)行情況進(jìn)行分析，總的供電費(fèi)用可以表示為：

(11)

式中：C為全天總供電費(fèi)用；Czi為購電價(jià)格；Cwi為風(fēng)力機(jī)單位運(yùn)營成本；Cvi為光伏單位運(yùn)營成本；yi為狀態(tài)變量，當(dāng)系統(tǒng)從主電網(wǎng)購買電力時(shí)，其值為1，當(dāng)系統(tǒng)向主電網(wǎng)售出電力時(shí)，其值為0；Cgi為從主網(wǎng)購電價(jià)格；Cmi為向主網(wǎng)售電價(jià)格；Cci為將蓄電池從充電狀態(tài)切換到放電狀態(tài)所需要的成本；Pci為蓄電池的放電功率；CDR為激勵(lì)型需求響應(yīng)成本(補(bǔ)償費(fèi)用)。

(12)

式中：λ為衡量可轉(zhuǎn)移負(fù)荷所獲得經(jīng)濟(jì)補(bǔ)償?shù)南禂?shù)；Pin，i為轉(zhuǎn)入的負(fù)荷量；Pout，i為轉(zhuǎn)出的負(fù)荷量。當(dāng)系統(tǒng)總出力大于負(fù)荷時(shí)，Pin，i≥0，Pout，i=0；當(dāng)系統(tǒng)總出力小于負(fù)荷時(shí)，Pin，i=0，Pout，i≥0。

2.2 約束條件

2.2.1 有功功率平衡約束

有功功率平衡約束為：

Pwi+Pvi+Pzi=Qi.

(13)

考慮需求響應(yīng)后，此約束可表示為：

Pwi+Pvi+Pzi=Q′i.

(14)

轉(zhuǎn)入轉(zhuǎn)出負(fù)荷、與主網(wǎng)交易功率約束分別見式(15)、(16)：

(15)

0≤|Pzi|≤Pzimax.

(16)

式中：Pin，max為最大轉(zhuǎn)入負(fù)荷；Pout，max為最大轉(zhuǎn)出負(fù)荷；Pzimax為與主網(wǎng)交換最大功率限制。

2.2.2 機(jī)組運(yùn)行約束

風(fēng)電機(jī)組、光伏機(jī)組、蓄電池運(yùn)行約束為：

0≤Pwi≤Pwimax，

(17)

0≤Pvi≤Pvimax，

(18)

Smin≤Si≤Smax.

(19)

式(17)—(19)中：Si為蓄電池荷電狀態(tài)；Smin為荷電狀態(tài)下限；Smax為荷電狀態(tài)上限。

3 算例分析

3.1 算法流程

采用反向?qū)W習(xí)的混沌映射自適應(yīng)粒子群算法對(duì)風(fēng)-光-儲(chǔ)能和需求響應(yīng)不同組合搭配的5種調(diào)度方案進(jìn)行探討，算法的具體流程[21]如圖1所示。

圖1 算法流程Fig.1 Algorithm flowchart

圖2 系統(tǒng)結(jié)構(gòu)Fig.2 System structure

圖3 新能源預(yù)測出力和凈負(fù)荷Fig.3 Predicted output and net load of new energy

圖4 風(fēng)-光參與運(yùn)行結(jié)果Fig.4 Wind and photovoltaic participating in the operation

3.2 案例分析

在本案例中，采集河北邢臺(tái)某地區(qū)在夏季典型日的風(fēng)-光和負(fù)荷情況作為數(shù)據(jù)來源。風(fēng)力機(jī)的裝機(jī)容量為 360 kW，運(yùn)維成本為 0.52元/kWh，切入風(fēng)速為5 m/s，切出風(fēng)速為25 m/s，額定風(fēng)速為13 m/s。光伏發(fā)電設(shè)備的裝機(jī)容量為260 kW，運(yùn)維成本為0.75元/kWh。蓄電池的額定容量為 700 kWh，電池的荷電狀態(tài)運(yùn)行范圍為 0.4～0.9，初始荷電狀態(tài)為 0.4，充電狀態(tài)和放電狀態(tài)相互轉(zhuǎn)變的成本為 0.1元/kWh。蓄電池在1 h內(nèi)最多可以充放電功率為儲(chǔ)能容量的20%。VPP與主網(wǎng)之間的允許交換功率不超過 200 kW。系統(tǒng)所涉及的設(shè)備如圖 2所示，售電和購電價(jià)格在一天內(nèi)的波動(dòng)見表1。

表1 售電和購電價(jià)格Tab.1 Electricity sale and purchase prices

3.3 基于各種條件不確定的影響分析

3.3.1 主網(wǎng)供電的影響

未來24 h內(nèi)風(fēng)-光總預(yù)測出力、傳統(tǒng)負(fù)荷和凈負(fù)荷變化情況如圖 3所示。只有主網(wǎng)供電時(shí)，在13 時(shí)段和19時(shí)段達(dá)到用電高峰，在16時(shí)段出現(xiàn)用電峰谷，沒有風(fēng)、光能可利用時(shí)，VPP需要從主網(wǎng)中獲取所需電能，此時(shí)電廠所需的功率與負(fù)荷相等。

3.3.2 風(fēng)-光參與調(diào)度的影響

在存在風(fēng)-光供電時(shí)，分別根據(jù)雙參數(shù)威布爾分布函數(shù)和Beta分布函數(shù)計(jì)算風(fēng)-光出力，電量分布如圖 4所示。由圖 4可以看出：將計(jì)算后的風(fēng)力機(jī)、光伏出力與主網(wǎng)出力相加后和負(fù)荷需求量相比，超出負(fù)荷部分將出售給電網(wǎng)，未滿足部分將從主網(wǎng)購買；風(fēng)-光機(jī)組參與供電后，1—8 時(shí)段風(fēng)-光出力能夠滿足負(fù)荷需求，此時(shí)功率全部來源于風(fēng)力機(jī)；在8—15時(shí)段風(fēng)-光尚有剩余，此時(shí)以0.42元/kWh出售給主網(wǎng)。在15—24時(shí)段，光伏和風(fēng)力機(jī)出力都呈現(xiàn)下降趨勢(shì)，但負(fù)荷仍保持在較高水平，此時(shí)需要從主網(wǎng)購電來滿足負(fù)荷。與只有主網(wǎng)參與時(shí)相比，加入風(fēng)力機(jī)、光伏后，從主網(wǎng)購電量明顯減少，并且可以有一定的售電收入。

3.3.3 儲(chǔ)能參與調(diào)度的影響

將計(jì)算后的風(fēng)力機(jī)、光伏出力與主網(wǎng)出力、儲(chǔ)能設(shè)備充、發(fā)電量相加后和負(fù)荷需求量相比，超出負(fù)荷部分將出售給電網(wǎng)，未滿足部分將從主網(wǎng)購買，數(shù)據(jù)結(jié)果如圖5所示。在風(fēng)-光-儲(chǔ)都參與的情況下，在1—6時(shí)段負(fù)荷仍然可以由風(fēng)力機(jī)出力來滿足，此時(shí)儲(chǔ)能設(shè)備作為負(fù)荷，這部分由主網(wǎng)滿足；負(fù)荷需求在16—24時(shí)段達(dá)到高峰，此時(shí)儲(chǔ)能作為供電方參與供電。與只有風(fēng)-光參與時(shí)相比，儲(chǔ)能設(shè)備在風(fēng)-光出力不能滿足負(fù)荷需求時(shí)，可以作為后備能源來填補(bǔ)。

圖5 風(fēng)-光-儲(chǔ)參與運(yùn)行結(jié)果Fig.5 Wind，photovoltaic and energy storage participating in the operation

3.3.4 需求響應(yīng)參與調(diào)度的影響

對(duì)于需求響應(yīng)參與的模型，先對(duì)負(fù)荷進(jìn)行K-Means聚類，將負(fù)荷分為高峰、平段、低谷3個(gè)時(shí)段，見表2。在能源供應(yīng)調(diào)度期間，通過基于價(jià)格型需求策略得到需求響應(yīng)后的負(fù)荷，可以使用激勵(lì)型需求響應(yīng)直接控制的方式來實(shí)現(xiàn)負(fù)荷轉(zhuǎn)移，并向參與者提供經(jīng)濟(jì)或其他形式的補(bǔ)償。

表2 K-Means聚類對(duì)負(fù)荷進(jìn)行分類結(jié)果Tab.2 Load classification results of K-means clustering

在考慮風(fēng)-光出力對(duì)系統(tǒng)運(yùn)行的基礎(chǔ)上，分別考慮電價(jià)型需求響應(yīng)和激勵(lì)型需求響應(yīng)后的系統(tǒng)運(yùn)行數(shù)據(jù)分別如圖6和圖7所示?？紤]了需求響應(yīng)的影響后負(fù)荷的變化情況如圖8所示。由圖 6可以看出：在風(fēng)-光和需求響應(yīng)參與的情況下，考慮電價(jià)后，在10—21時(shí)段用電量明顯降低，并保持在一個(gè)相對(duì)平穩(wěn)的用電范圍，起到了削峰填谷作用。由圖 7可以看出：在風(fēng)力機(jī)、光伏和需求響應(yīng)都參與的情況下，經(jīng)過可控負(fù)荷調(diào)整后，未滿足負(fù)荷量幾近平穩(wěn)。由圖 8可以看出：經(jīng)過需求響應(yīng)的調(diào)整變化后，在不影響用戶使用的基礎(chǔ)上將17—22時(shí)段的高負(fù)荷利用部分轉(zhuǎn)移到了24—4 時(shí)段的低負(fù)荷時(shí)段，負(fù)荷的整體變化情況更趨于平穩(wěn)，緩解了電網(wǎng)供電的壓力。

圖6 風(fēng)-光和電價(jià)型需求響應(yīng)參與運(yùn)行結(jié)果Fig.6 Wind，photovoltaic and electricity price demand response participating in the operation

圖7 風(fēng)-光和激勵(lì)型需求響應(yīng)參與運(yùn)行結(jié)果Fig.7 Wind and incentivized demand response participating in the operation

圖8 需求響應(yīng)調(diào)整前后負(fù)荷變化Fig.8 Load changes before and after demand response adjustment

3.3.5 綜合因素對(duì)調(diào)度的影響

綜合考慮風(fēng)-光-儲(chǔ)和需求響應(yīng)對(duì)系統(tǒng)運(yùn)行的整體運(yùn)行結(jié)果如圖9所示。由圖 9可以看出：在風(fēng)-光-儲(chǔ)和需求響應(yīng)都參與的情況下，與風(fēng)-光-儲(chǔ)參與運(yùn)行結(jié)果圖 5數(shù)據(jù)相比，負(fù)荷缺額量明顯降低，儲(chǔ)能也得到了更好的利用，其他方面效果也有明顯的改善。

圖9 全因素參與運(yùn)行結(jié)果圖Fig.9 All factors participating in the operation

將主網(wǎng)和風(fēng)-光參與確定為情景1；主網(wǎng)和風(fēng)-光-儲(chǔ)參與的情況為情景2；主網(wǎng)和風(fēng)-光參與的情況下考慮需求響應(yīng)的作用為情景3；主網(wǎng)、風(fēng)-光-儲(chǔ)都參與的情況下考慮需求響應(yīng)的作用為情景4。不同情景下的VPP運(yùn)行成本、用戶舒適度等數(shù)據(jù)見表3。

表3 不同情景的運(yùn)行結(jié)果Tab.3 Operating results of different scenarios

從成本方面來看，考慮多種供電方式的協(xié)同作用能夠有效減少運(yùn)行成本的投入，利用發(fā)電設(shè)備及儲(chǔ)能設(shè)備比從主網(wǎng)直接購電更加優(yōu)惠。將情景1與情景3、情景2與情景4分別進(jìn)行對(duì)比，改變的運(yùn)行條件只是加入了需求響應(yīng)的調(diào)度作用，2組數(shù)據(jù)對(duì)比之后可以看出需求響應(yīng)的調(diào)度作用能有效提升用戶用電滿意度。情景1與情景2對(duì)比，加入儲(chǔ)能設(shè)備后能有效減少運(yùn)行時(shí)未滿足負(fù)荷。情景4與情景1相比成本降低了9.87%、用戶舒適度提高了17.44%，情景4與情景2相比成本降低了4.47%，用戶舒適度提高了3.51%。綜合考慮各方面運(yùn)行結(jié)果，可以得出考慮風(fēng)-光-儲(chǔ)和需求響應(yīng)的多方面協(xié)同作用，無論是在成本還是用戶舒適度方面都能得到更優(yōu)的結(jié)果。

4 結(jié)束語

為了降低VPP的運(yùn)行成本，提高用戶舒適度，本文建立了一個(gè)考慮風(fēng)-光-儲(chǔ)和需求響應(yīng)的VPP優(yōu)化調(diào)度模型。引入儲(chǔ)能設(shè)備和激勵(lì)型需求響應(yīng)，讓用戶在考慮分時(shí)電價(jià)的影響后主動(dòng)調(diào)整用電時(shí)間，從而起到消峰填谷的作用。利用反向?qū)W習(xí)的混沌映射自適應(yīng)粒子群算法對(duì)運(yùn)行成本進(jìn)行仿真分析，經(jīng)過計(jì)算風(fēng)-光-儲(chǔ)和需求響應(yīng)都參與的方案與情景1相比不僅可以將VPP的運(yùn)行成本降低9.87%，還可以將用戶舒適度提高17.44%。