何 意，郭 蘇，彭懷午，段楊龍，宋國濤，阿依努爾·庫爾班

(1.河海大學(xué)水利水電學(xué)院，南京 210024;2.河海大學(xué)能源與電氣學(xué)院，南京 211100;3.中國電建集團西北勘測設(shè)計研究院有限公司,西安 710065)

0 前 言

由于化石能源短缺及環(huán)境污染問題，可再生能源發(fā)電技術(shù)發(fā)展受到了廣泛關(guān)注，其中風(fēng)力發(fā)電及光伏發(fā)電是目前較為成熟的可再生能源發(fā)電形式[1]。由于可再生能源的間歇性和不穩(wěn)定性，儲能系統(tǒng)通常被用于調(diào)節(jié)風(fēng)電-光伏發(fā)電系統(tǒng)的功率波動，從而提高可再生能源的利用效率及滲透率水平[2]。常用的儲能技術(shù)包括蓄電池、抽水蓄能、壓縮儲能、氫儲能等，其中蓄電池和抽水蓄能是較為成熟、應(yīng)用廣泛的新型儲能方式[3]。

多能互補系統(tǒng)可利用不同資源之間的互補特性提高系統(tǒng)的供電可靠性，其中風(fēng)電-光伏-儲能聯(lián)合發(fā)電系統(tǒng)的規(guī)劃及調(diào)度問題受到了國內(nèi)外有關(guān)研究人員的廣泛關(guān)注。Javed等[4]評估了離網(wǎng)型風(fēng)電-光伏-蓄電池聯(lián)合系統(tǒng)的技術(shù)經(jīng)濟性能，并通過遺傳算法求解了系統(tǒng)容量優(yōu)化問題。Guo等[5]研究了并網(wǎng)型風(fēng)電-光伏-熔鹽儲熱聯(lián)合系統(tǒng)的多目標(biāo)容量優(yōu)化，并通過案例分析得出了熔鹽儲熱的技術(shù)經(jīng)濟性能優(yōu)于蓄電池儲能。Ma等[6]研究了抽水蓄能電站在獨立型微網(wǎng)中的技術(shù)經(jīng)濟性，結(jié)果表明抽水蓄能電站可有效調(diào)節(jié)風(fēng)電、光伏的出力波動并滿足負荷需求。Baghaee等[7]研究了獨立型風(fēng)電-光伏-儲氫-燃料電池聯(lián)合系統(tǒng)的容量優(yōu)化配置，并通過多目標(biāo)粒子群優(yōu)化算法求解。上述研究均聚焦于單儲能系統(tǒng)的容量優(yōu)化問題，單一儲能難以同時滿足風(fēng)電光伏功率調(diào)節(jié)過程中的功率密度及容量密度的雙重需求[8]，然而多種儲能技術(shù)的聯(lián)合運行可有效利用不同儲能的互補特性，提高儲能系統(tǒng)的調(diào)節(jié)能力及綜合性能。

國內(nèi)外學(xué)者已對聯(lián)合儲能系統(tǒng)進行了初步的研究。常見的聯(lián)合儲能形式包括蓄電池-超級電容、蓄電池-氫儲能、蓄電池-飛輪儲能以及蓄電池-抽水蓄能等。Javed等[8]研究了離網(wǎng)型抽水蓄能-電池聯(lián)合儲能的協(xié)同運行策略，但未研究聯(lián)合儲能的容量優(yōu)化配置。Zhang等[9]研究了蓄電池-氫儲能聯(lián)合系統(tǒng)的容量優(yōu)化問題，未考慮聯(lián)合儲能的協(xié)同優(yōu)化策略，然而容量配置和運行策略之間相互耦合，規(guī)劃設(shè)計階段需要對兩者進行協(xié)同優(yōu)化。本文主要研究離網(wǎng)型風(fēng)電-光伏-抽水蓄能-蓄電池聯(lián)合發(fā)電系統(tǒng)的容量及運行兩階段協(xié)同優(yōu)化。其中，一階段容量優(yōu)化模型以最小化聯(lián)合系統(tǒng)的平準(zhǔn)化度電成本為目標(biāo)，利用量子粒子群算法求解系統(tǒng)容量優(yōu)化配置。二階段運行優(yōu)化模型以最小化失負荷率為目標(biāo)，考慮不同儲能的運行約束，利用混合整數(shù)線性規(guī)劃求解聯(lián)合儲能協(xié)同運行策略。此外，本文研究多種智能優(yōu)化算法的求解性能，并通過案例分析驗證聯(lián)合儲能的有效性。

1 系統(tǒng)結(jié)構(gòu)

風(fēng)電-光伏-抽水蓄能-蓄電池聯(lián)合發(fā)電系統(tǒng)主要包括風(fēng)電場、光伏場、直流母線、蓄電池、水庫、水泵及水輪機，其簡化結(jié)構(gòu)如圖1所示。當(dāng)風(fēng)電光伏的聯(lián)合出力大于負荷需求時，多余電能可直接以電能形式儲存在蓄電池中，也可通過水泵抽水工況以水勢能形式儲存在水庫中；當(dāng)風(fēng)電光伏的聯(lián)合出力小于負荷需求時，不足電能可通過蓄電池放電補充，也可利用水庫放水工況經(jīng)水輪機發(fā)電補充。由于水泵及水輪機的技術(shù)運行特性，20%額定工況下水泵及水輪機無法正常運行，該階段通過蓄電池單獨調(diào)節(jié)，其他工況下可通過抽水蓄能-蓄電池協(xié)同調(diào)節(jié)。

圖1 風(fēng)電-光伏-抽水蓄能-蓄電池聯(lián)合發(fā)電系統(tǒng)結(jié)構(gòu)圖

2 組件數(shù)學(xué)模型

2.1 風(fēng)電模型

風(fēng)電的輸出功率取決于風(fēng)機輪轂高度處的實際風(fēng)速以及風(fēng)功率曲線，其出力計算方法如下所示[5]。

(1)

式中:PW、PWR分別為風(fēng)機理論出力和風(fēng)機額定功率，kW；v、vci、vci、vr分別為輪轂處實際風(fēng)速、風(fēng)機切入風(fēng)速、切出風(fēng)速、額定風(fēng)速，m/s。

2.2 光伏模型

光伏的輸出功率取決于光伏傾斜面上的輻照度以及運行溫度，其出力計算方法如下所示[5]。

(2)

(3)

式中:PPV、PPVR分別為光伏理論出力、光伏額定功率，kW；IT、Istc、Iref分別為實際輻照強度、標(biāo)準(zhǔn)測試條件輻照強度、參考輻照強度，W/m2；TPV、Tstc、Tamb、NOCT、Tref分別為光伏實際運行溫度、標(biāo)準(zhǔn)測試條件溫度、環(huán)境溫度、光伏額定運行溫度、參考溫度，℃；εT為溫度修正系數(shù)。

2.3 蓄電池模型

蓄電池具有靈活的調(diào)節(jié)能力，能夠有效消納風(fēng)電光伏的多余電能并用于補充不足負荷。蓄電池經(jīng)儲放電工況后的容量如下所示[4]。

Ebat(t)=Ebat(t-1)+Pbat.c(t)·ηbat.c△t

(4)

(5)

式中:Ebat(t-1)和Ebat(t)分別為儲放電前、后蓄電池容量，kWh；Pbat.c(t)和Pbat.d(t)分別為蓄電池充電功率及放電功率，kW；ηbat.c和ηbat.d分別為蓄電池充電效率及放電效率。

2.4 抽水蓄能模型

抽水蓄能電站具有技術(shù)成熟，響應(yīng)速度快等特點，抽水蓄能電站經(jīng)儲放電工況后的容量如下所示[6]。

(6)

(7)

式中:VUR(t-1)和VUR(t)分別為儲放能前后的水庫容量，m3；PP(t)和Pt(t)分別為抽水蓄能電站儲能功率及發(fā)電功率，kW；ηP和ηt分別為水泵效率及水輪機效率；ρ為水密度，取1 000 kg/m3；g為重力加速度，取9.81 m/s2；h為額定水頭，取40 m。

3 優(yōu)化模型

本文的兩階段優(yōu)化模型分別為容量配置優(yōu)化模型及聯(lián)合儲能協(xié)同運行優(yōu)化模型。容量優(yōu)化配置模型以系統(tǒng)經(jīng)濟性最優(yōu)為目標(biāo)，通過量子粒子群優(yōu)化算法進行隨機尋優(yōu)求解。聯(lián)合儲能協(xié)同運行優(yōu)化模型以系統(tǒng)可靠性最優(yōu)為目標(biāo)，通過混合整數(shù)線性規(guī)劃進行精確求解。其中，一階段模型的容量配置作為二階段模型的邊界條件；二階段的目標(biāo)函數(shù)作為返回值嵌入到一階段目標(biāo)函數(shù)中。詳細優(yōu)化目標(biāo)函數(shù)及約束條件如下所示。

3.1 運行優(yōu)化模型

下層運行優(yōu)化模型以最小化失負荷率為目標(biāo)，考慮不同儲能技術(shù)的運行約束，優(yōu)化聯(lián)合儲能的協(xié)同運行策略，即每時刻各儲能系統(tǒng)的充放電功率。

失負荷率(Loss of Power Supply Probability，LPSP)表示系統(tǒng)未滿足負荷需求的概率，其計算方法如下所示。

(8)

式中：Pw(t)、Ppv(t)、PES(t)、Pload(t)分別為t時刻風(fēng)電出力，光伏出力，聯(lián)合儲能出力及負荷需求，kW；T為仿真時長，本文取24 h。

蓄電池及抽水蓄能的運行約束如下所示。其中，公式(9)～(10)表示蓄電池及水庫每時刻容量約束；公式(11)～(13)表示蓄電池，水泵及水輪機每時刻功率約束；公式(14)～(15)表示蓄電池及水庫初始末時刻容量約束。

(1-DOD)·Ebat.max≤Ebat(t)≤Ebat.max

(9)

VUR.min≤VUR(t)≤VUR.max

(10)

0≤Pbat.c(t),Pbat.d(t)≤Pbat.max

(11)

Pp.min≤Pp(t)≤Pp.max

(12)

Pt.min≤Pt(t)≤Pt.max

(13)

(1-ε)·Ebat(0)≤Ebat(T)-Ebat(0)≤(1+ε)·Ebat(0)

(14)

(1-ε)·VUR(0)≤VUR(T)-VUR(0)≤(1+ε)·VUR(0)

(15)

式中:Ebat.max為蓄電池最大容量，kWh；Pbat.max為蓄電池額定功率，kW；DOD為放電深度；VUR.min和VUR.max分別為水庫最小容量及最大容量，m3；Pp.max及Pt.max分別表示水泵及水輪機的額定功率，kW；Pp.min及Pt.min分別表示水泵及水輪機的運行下限，本文取20%；ε為儲能容量允許偏差，本文取10%。

3.2 容量優(yōu)化模型

上層容量優(yōu)化模型以最小化平準(zhǔn)化度電成本為目標(biāo)，考慮二階段優(yōu)化結(jié)果LPSP作為罰函數(shù)，優(yōu)化各組件的容量配置。

平準(zhǔn)化度電成本(Levelized Cost of Energy，LCOE)表示系統(tǒng)全生命周期總成本與總發(fā)電量的比值，其計算方法如下所示。

(16)

式中:IC和AC分別為系統(tǒng)初始投資成本及年運行維護成本，元；Ew、Epv、EES分別為風(fēng)機、光伏、儲能的首年發(fā)電量，kWh；dpv為光伏衰減率，本文風(fēng)機和儲能衰減率忽略不計；NS為系統(tǒng)設(shè)計壽命，本文取20年；i為折現(xiàn)率，本文取6.4%。

本文的優(yōu)化模型如下所示,其中公式(17)是下層運行優(yōu)化目標(biāo)函數(shù),公式(18) 是上層容量優(yōu)化目標(biāo)函數(shù)。

f1=minLPSP[Pbat.c(t),Pbat.d(t),Pp(t),Pt(t)]

(17)

f2=minLCOE[Cw,Cpv,Pp.max,Pt.max,VUR.max,Pbat.max,Ebat.max]

+αf1

(18)

式中:Cw、Cpv分別為風(fēng)電及光伏裝機容量，kW；上層決策變量為系統(tǒng)各組件的容量及額定功率；下層優(yōu)化的決策變量為每時刻蓄電池及抽水蓄能的儲放電功率。

3.3 優(yōu)化算法

下層聯(lián)合儲能協(xié)同運行優(yōu)化模型采用混合整數(shù)線性規(guī)劃進行精確求解，上層容量優(yōu)化模型利用量子粒子群算法(Quantum Particle Swarm Optimization，QPSO)進行隨機尋優(yōu)求解。PSO具有計算簡單易于實現(xiàn)的優(yōu)點，但其通過個體最優(yōu)以及全局最優(yōu)學(xué)習(xí)的方式更新粒子位置缺少隨機性，易陷入局部最優(yōu)。QPSO引入量子行為理論，提高算法的隨機搜索能力[10]。PSO及QPSO的位置更新如公式(19) ～ (20)所示。

xi+1=xi+ω·vi-1+rand·c1·(pbest-xi)

+rand·c2·(gbest-xi)

(19)

xi+1=[φ·pbest+(1-φ)·gbest]

±β|Mbest-xi|In(1/μ)

(20)

式中：pbest、gbest、Mbest為個體最優(yōu)，全局最優(yōu)及個體最優(yōu)均值；ω、c1、c2是PSO的控制參數(shù)慣性權(quán)重和學(xué)習(xí)因子;β是QPSO的控制參數(shù)創(chuàng)新因子。

4 實驗結(jié)果及討論

4.1 數(shù)據(jù)準(zhǔn)備

本文以某地區(qū)(25°04′N,67°56′E)風(fēng)光資源及負荷數(shù)據(jù)為案例，研究風(fēng)電-光伏-抽水蓄能-蓄電池聯(lián)合發(fā)電系統(tǒng)的容量及運行兩階段協(xié)同優(yōu)化。系統(tǒng)組件經(jīng)濟性參數(shù)如表1所示，代表日風(fēng)電光伏出力及負荷曲線如圖2所示。

表1 系統(tǒng)組件經(jīng)濟性參數(shù)表

圖2 代表日風(fēng)電光伏出力及負荷曲線圖

4.2 算法性能分析

為了研究QPSO算法的性能，QPSO、PSO及遺傳算法(Genetic Algorithm，GA)多種單目標(biāo)進化算法分別求解該容量優(yōu)化問題。不同算法隨迭代過程變化的收斂曲線如圖3所示。根據(jù)圖3可知，QPSO的收斂速度最快，且收斂結(jié)果最優(yōu)；PSO的收斂速度以及收斂結(jié)果均次之;GA的收斂速度最慢且收斂結(jié)果最差。因此，QPSO在求解該容量優(yōu)化問題時具有最佳的算法性能。

圖3 算法收斂曲線圖

4.3 優(yōu)化結(jié)果

該案例分析的容量優(yōu)化最優(yōu)解如表2所示，聯(lián)合儲能協(xié)同優(yōu)化運行結(jié)果如圖4～5所示。根據(jù)表2可知，當(dāng)風(fēng)電、光伏、蓄電池及水庫的容量分別為207 MW、156 MW、354 MWh、1 500 m3時，聯(lián)合儲能系統(tǒng)可完成凈負荷的完全調(diào)節(jié)，此時平準(zhǔn)化度電成本及負荷缺電率分別為0.2018 元/kWh和0%；根據(jù)圖4～5可知，23:00-4:00，夜間風(fēng)電出力較大且負荷需求較小導(dǎo)致風(fēng)電光伏總出力大于負荷需求，因此抽水蓄能-蓄電池聯(lián)合儲能充電消納多余電能，此時抽水蓄能及蓄電池剩余容量均上升；5:00-7:00，風(fēng)電出力減少且負荷需求增大導(dǎo)致風(fēng)電光伏總出力小于負荷需求，因此聯(lián)合儲能放電補充不足負荷，此時抽水蓄能及蓄電池剩余容量均下降；8:00-17:00，光伏出力增大導(dǎo)致風(fēng)電光伏總出力大于負荷需求，聯(lián)合儲能充電；18:00-22:00，夜間無光伏出力且為負荷峰值階段，聯(lián)合儲能放電。此外，抽水蓄能作為主要調(diào)節(jié)單元，蓄電池用于補充抽水蓄能運行，且運行周期末儲能容量滿足初始末容量約束，可保證多時段的連續(xù)調(diào)節(jié)。以上結(jié)果表明：通過聯(lián)合儲能的協(xié)同運行可有效調(diào)節(jié)可再生能源波動出力與負荷需求之間的差值，實現(xiàn)無棄風(fēng)棄光損失且無不足負荷，充分提高可再生能源的供電可靠性及消納率水平。

圖4 抽水蓄能-蓄電池聯(lián)合儲能充放電功率圖

圖5 抽水蓄能-蓄電池聯(lián)合儲能剩余容量圖

表2 容量優(yōu)化結(jié)果表

5 結(jié) 論

本文研究了風(fēng)電-光伏-抽水蓄能-蓄電池聯(lián)合發(fā)電系統(tǒng)的容量及運行兩階段協(xié)同優(yōu)化，其以系統(tǒng)經(jīng)濟性和可靠性為優(yōu)化目標(biāo)，利用量子粒子群算法和混合整數(shù)線性規(guī)劃方法求解容量配置及聯(lián)合儲能運行策略。通過案例分析得出以下結(jié)論：

(1) 量子粒子群算法的收斂速度及收斂結(jié)果均優(yōu)于常規(guī)粒子群算法及遺傳算法，說明量子粒子群算法在求解該容量優(yōu)化問題時具有最佳的算法性能。

(2) 當(dāng)風(fēng)電、光伏、蓄電池及水庫的容量分別為207 MW、156 MW、354 MWh和1 500 m3時，聯(lián)合儲能系統(tǒng)可完成凈負荷的完全調(diào)節(jié)。

(3) 抽水蓄能-蓄電池聯(lián)合儲能可通過協(xié)同優(yōu)化運行有效消納多余電能并補充不足負荷，實現(xiàn)零棄風(fēng)棄光損失及零不足負荷。