佟 曦, 陳 衡, 茍凱杰, 徐 鋼, 劉文毅, 張國強

(華北電力大學(xué) 能源動力與機械工程學(xué)院, 北京 102206)

隨著我國“雙碳”目標(biāo)的提出,各行各業(yè)都在積極開展碳減排工作。電力行業(yè)的碳排放量約占我國能源活動二氧化碳量的40%,減排潛力巨大,其低碳轉(zhuǎn)型發(fā)展對于我國實現(xiàn)“碳中和”戰(zhàn)略目標(biāo)至關(guān)重要。為實現(xiàn)該戰(zhàn)略目標(biāo),我國電力行業(yè)一方面進行火電機組的靈活性改造,提高機組調(diào)峰能力[1];另一方面大力發(fā)展風(fēng)能、光伏等新能源發(fā)電技術(shù),并且大規(guī)模接入電網(wǎng),以減少化石能源等消耗[2-5]。由于風(fēng)能、光伏等新能源具有出力不確定性,以及發(fā)電的間歇性和波動性,其大規(guī)模接入電網(wǎng)后,給電網(wǎng)的電壓水平、運行控制等帶來巨大困擾[6-10]。儲能設(shè)備可以在時間和空間的尺度上快速進行能量調(diào)度,解決調(diào)峰難的問題,從而維持系統(tǒng)的穩(wěn)定性,因此近年來備受關(guān)注和使用[11-12]。

研究碳交易機制對新能源電力系統(tǒng)運行優(yōu)化的影響至關(guān)重要,陳錦鵬等[13]構(gòu)建了綜合能源系統(tǒng)(IES)模型,并且引入階梯式碳交易機制,構(gòu)建以購能成本、碳排放成本、棄風(fēng)成本最小的低碳經(jīng)濟運行目標(biāo)進行優(yōu)化求解,驗證了所提策略的有效性。陳彥奇等[14]提出了計及階梯型碳交易機制的綜合能源系統(tǒng)優(yōu)化模型,以系統(tǒng)運行維護成本和碳交易成本之和最小為目標(biāo),分析典型日下不同基準(zhǔn)碳價對系統(tǒng)運行方式和經(jīng)濟性的影響。王俐英等[15]將需求響應(yīng)和碳交易機制引入綜合能源系統(tǒng)運行調(diào)度中,并基于信息間隙決策理論(IGDT)構(gòu)建了考慮階梯式碳交易機制以及需求響應(yīng)的綜合能源系統(tǒng)雙層隨機優(yōu)化調(diào)度模型,最后通過求解證明了所提系統(tǒng)的低碳環(huán)保。Cao等[16]提出了面向運營商和多個虛擬電廠(MVPP)的Stackelberg博弈策略,考慮階梯式碳價格機制和雙補償需求響應(yīng)機制,采用克里金元模型擬合虛擬電廠(VPP)的能量調(diào)度模型,最后分析了多個算例,驗證了所提策略能夠有效促進MVPP的能量互補和碳減排。

平準(zhǔn)化度電成本(LCOE),也叫電力平均化成本,是用于衡量發(fā)電項目在其整個生命周期內(nèi)每單位電力產(chǎn)生的平均成本,在研究系統(tǒng)經(jīng)濟性中起到了至關(guān)重要的作用。劉漂[17]從全生命周期角度研究了風(fēng)電成本的構(gòu)成,并針對現(xiàn)有LCOE模型提出了改進的LCOE分析計算模型,通過物理建模驗證了所提出的LCOE模型的有效性。Mostafa等[18]建立了電池儲能成本模型,并考慮了長期、中期和短期應(yīng)用以及技術(shù)特征,使用儲能平準(zhǔn)化度電成本(LCCOS)和LCOE進行技術(shù)經(jīng)濟評估,此外,還對影響存儲成本的控制因素進行了敏感性分析。Li等[19]建立了實現(xiàn)碳捕集與封存(CCS)一體化氣化聯(lián)合循環(huán)(IGCC)建模、預(yù)測和優(yōu)化的計算方法,選擇碳捕集率(CCR)和LCOE作為優(yōu)化目標(biāo),評估了多目標(biāo)優(yōu)化的最優(yōu)解集和碳交易市場的效果。

當(dāng)前大多數(shù)研究只考慮單一儲能系統(tǒng)接入電網(wǎng)的運行優(yōu)化及配置優(yōu)化,考慮儲能系統(tǒng)LCOE的運營策略研究很少。因此,筆者構(gòu)建了含火力發(fā)電、新能源(風(fēng)能、光伏)發(fā)電,以及抽水蓄能和電化學(xué)儲能的多能互補系統(tǒng)數(shù)學(xué)模型;同時,考慮LCOE及碳排放成本,引入階梯式碳交易機制,以總成本之和最小為優(yōu)化目標(biāo),利用MATLAB的Cplex求解器進行優(yōu)化求解,并且設(shè)置不同場景及不同碳交易機制,以比較不同調(diào)度方式的優(yōu)劣,最終提出成本最低、碳排放最少、收益最大化的電力系統(tǒng)下電化學(xué)儲能及抽水蓄能的協(xié)同優(yōu)化調(diào)度策略,為新型電力系統(tǒng)的儲能規(guī)劃提供參考,助力電網(wǎng)高效清潔化及經(jīng)濟最優(yōu)化的發(fā)展。

1 階梯式碳交易機制模型

階梯式碳交易是在碳市場中設(shè)立多個不同層級的碳排放配額,每個層級對應(yīng)不同的行業(yè)或行業(yè)細(xì)分。不同行業(yè)的企業(yè)擁有不同數(shù)量的碳排放配額,企業(yè)在交易市場上可以自由買賣碳排放配額。與普通碳交易模式相比,階梯式碳交易模式的優(yōu)勢在于更精準(zhǔn)地針對不同行業(yè)的減排需求,可以更有效地推動高排放行業(yè)的減排工作[20]。階梯式碳交易機制模型主要包含無償碳排放配額模型、實際碳排放模型和階梯式碳交易模型[21]。

1.1 碳排放配額模型

所構(gòu)建的系統(tǒng)中,碳排放源有燃煤機組及系統(tǒng)外購電力。所構(gòu)建系統(tǒng)的碳排放配額為

ECEQ=ECFP+ESOP

(1)

式中:ECEQ為系統(tǒng)碳排放配額,t;ECFP為燃煤機組碳排放配額,t;ESOP為系統(tǒng)外購電力碳排放配額,t。

燃煤機組碳排放配額為

(2)

式中:Ee(t)為t時刻機組供電CO2配額,t;Eh(t)為t時刻機組供熱CO2配額,t;T為周期;Δt為時間長度,h。

其中,t時刻機組供電CO2配額為

Ee(t)=Qe(t)×Be×F1×Fr×Ff

(3)

式中:Qe(t)為t時刻機組供電量,MW;Be為機組所屬類別的供電基準(zhǔn),t/MW;F1為機組冷卻方式修正系數(shù);Fr為機組供熱量修正系數(shù),燃煤機組供熱量修正系數(shù)為1-0.22×供熱比;Ff為機組負(fù)荷修正系數(shù)。

t時刻機組供熱CO2配額為

Eh(t)=Qh(t)×Bh

(4)

式中:Qh(t)為t時刻機組供熱量,GJ;Bh為機組所屬類別的供熱基準(zhǔn),t/GJ。

系統(tǒng)外購電力碳排放配額為

(5)

式中:ηe為單位電量碳排放配額系數(shù);Pbuy(t)為t時刻系統(tǒng)外購電量,MW。

1.2 實際碳排放模型

燃煤機組的碳排放量為

(6)

式中:EACE為燃煤機組的碳排放量,t;βe為燃煤機組單位電量碳排放系數(shù),t/MW。

系統(tǒng)外購電力碳排放量為

(7)

式中:EPCE為系統(tǒng)外購電力碳排放量,t;βs為系統(tǒng)外購電力單位電量碳排放系數(shù),t/MW。

不考慮新能源(風(fēng)能、光伏)發(fā)電及各類型儲能放電時的碳排放量,認(rèn)為其碳排放量為0,則當(dāng)其出力時,相當(dāng)于給系統(tǒng)減少碳排放量。減少的碳排放量為

(8)

因此,系統(tǒng)的實際碳排放量EACT為

EACT=EACE+EPCE-ERCE

(9)

1.3 階梯式碳交易模型

實際參與碳交易市場的碳交易額ECET為

ECET=EACT-ECEQ

(10)

當(dāng)ECET大于0時,表示主體當(dāng)前碳排放量高于碳排放配額,需要到碳市場購買碳排放配額;當(dāng)ECET小于0時,表示主體當(dāng)前碳排放量低于碳排放配額,可以將剩余碳排放配額拿到碳市場進行出售[22-23]?？紤]經(jīng)濟效益最大化以及提高電力企業(yè)減排的積極性,構(gòu)建階梯式碳交易計算模型,即

(11)

式中:Ccar為系統(tǒng)的碳交易成本;c為市場的碳交易基準(zhǔn)價格;d為碳排放區(qū)間長度;α為碳交易價格增長幅度。

傳統(tǒng)碳交易機制下,其碳交易模型為

(12)

2 協(xié)同調(diào)度優(yōu)化模型

2.1 LCOE模型

LCOE和全生命周期成本(LCC)都是在能源經(jīng)濟學(xué)和工程經(jīng)濟學(xué)中常用的概念,LCOE是衡量發(fā)電項目在其整個生命周期內(nèi)每單位電力產(chǎn)生的平均成本[24],而LCC是衡量整個項目或設(shè)備在其壽命周期內(nèi)的總成本。LCOE通常用于比較不同類型發(fā)電項目的經(jīng)濟性,而LCC用于做出長期投資決策,確保項目或設(shè)備的經(jīng)濟可行性。較低的LCOE意味著單位電力成本更低,使得該發(fā)電項目在市場上更具競爭力。在發(fā)電產(chǎn)業(yè)中,LCOE的下降通常與技術(shù)進步、規(guī)模效應(yīng)和政策支持等因素有關(guān)。因此,LCOE被廣泛用于評估各種發(fā)電技術(shù)的經(jīng)濟性,例如光伏、風(fēng)能、煤炭、天然氣等,以確定最具成本效益的能源解決方案。

LCOE的計算考慮了項目的總成本,包括建設(shè)、運營、維護、折舊成本等,以及項目預(yù)期的發(fā)電量[25-26]。這些成本在項目的壽命周期內(nèi)被折算,然后除以預(yù)期的發(fā)電量,從而得到每單位電力的平均成本。

LCOE可以表示為

(13)

式中:CLCOE為LCOE成本,元;I為初始投資,包括項目前期規(guī)劃建設(shè)以及設(shè)備購置安裝等費用,元;On為第n年的運行成本,包括在使用階段產(chǎn)生的費用,如能源消耗、人工勞動力、原材料等,元;Mn為第n年的維護成本,包括在使用過程中確保設(shè)備正常運轉(zhuǎn)和良好狀態(tài)的費用以及設(shè)備發(fā)生故障或損壞時所需的修理和維修費用,元;Gn為第n年的設(shè)備發(fā)電量,MW;N為年數(shù)。

考慮折現(xiàn)率時,LCOE的表達(dá)式為

(14)

式中:rLCOE為考慮折現(xiàn)率時的LCOE,元;r為折現(xiàn)率。

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

本文調(diào)度優(yōu)化以系統(tǒng)購能成本CB、CLCOE、碳排放成本CC以及棄風(fēng)棄光成本CWP之和C最小為目標(biāo)函數(shù),其數(shù)學(xué)表達(dá)式為

minC=CB+CLCOE+CC+CWP

(15)

系統(tǒng)購能成本CB主要包括外購電力成本、風(fēng)光儲的上網(wǎng)成本,以及火電上網(wǎng)成本,其表達(dá)式為

(16)

式中:Pg(t)為t時刻燃煤機組出力,MW;Cbuy(t)、Cpv、Cwt、Cs,dis,i(t)、Cs,cha,i(t)和Cg分別為t時刻系統(tǒng)購電價格、光伏上網(wǎng)電價、風(fēng)電上網(wǎng)電價、儲能上網(wǎng)電價、儲能充電電價以及燃煤機組上網(wǎng)電價,元/(kW·h)。

LCOE包括火電、風(fēng)電、光伏以及儲能的單位標(biāo)準(zhǔn)度電成本之和,其表達(dá)式為

(17)

碳排放成本主要包括系統(tǒng)外購電力碳排放成本和燃煤機組碳排放成本。階梯式碳交易下的碳排放計算公式參考式(11);傳統(tǒng)碳交易下的碳排放的計算公式參考式(12)。

棄風(fēng)棄光成本主要為風(fēng)能及光伏所產(chǎn)生的電能不能被完全利用的部分,其成本公式為

(18)

式中:a為風(fēng)光懲罰系數(shù),元/(kW·h);Ppv,max(t)為t時刻光伏最大出力,MW;Pwt,max(t)為t時刻風(fēng)能最大出力,MW。

2.3 約束條件

燃煤電廠發(fā)電功率約束為

0.3Pg,max≤Pg(t)≤0.8Pg,max

(19)

式中:Pg,max為火電機組額定出力,MW。

光伏發(fā)電出力約束為

0≤Ppv(t)≤Ppv,max(t)

(20)

風(fēng)力發(fā)電出力約束為

0≤Pwt(t)≤Pwt,max(t)

(21)

抽水蓄能系統(tǒng)約束為

應(yīng)該說，“首都餐飲業(yè)品質(zhì)提升工作”是北京市餐飲業(yè)以習(xí)近平總書記關(guān)于食品安全系列講話精神為指導(dǎo)，全面實施國家市場監(jiān)管總局提出的《餐飲服務(wù)食品安全操作規(guī)范》的自律表現(xiàn)，是繼北京餐飲業(yè)實現(xiàn)“明廚亮灶”、“陽光餐飲”后的一次全面提升，是為迎接2022年冬奧會打造餐飲業(yè)的首都標(biāo)準(zhǔn)、北京品牌，使北京餐飲業(yè)成為彰顯首都文化魅力、良好生態(tài)環(huán)境、和諧文明社會、安定富裕生活的載體，成為體現(xiàn)首都城市內(nèi)在品質(zhì)的亮麗名片，在全國餐飲行業(yè)中做出了表率。下一步，中國烹飪協(xié)會將在北京市市場監(jiān)督管理局的指導(dǎo)下，繼續(xù)深入落實“首都餐飲業(yè)品質(zhì)提升工作”，力爭取得更大的成效，惠及更多的企業(yè)和消費者。

(22)

Pp,min

(23)

式中:Et為上水庫初始能量,MW·h;Emin為上水庫最小能量,MW·h;ηh為發(fā)電轉(zhuǎn)化效率;Ph(t)為t時刻抽水蓄能系統(tǒng)發(fā)電功率,MW;Ph,min為抽水蓄能系統(tǒng)發(fā)電功率最小值,MW;Ph,max為抽水蓄能發(fā)電功率最大值,MW;Pp(t)為t時刻水泵抽水功率,MW;Pp,min為水泵抽水功率最小值,MW;Pp,max為水泵抽水功率最大值,MW。

Eu,min

(24)

Ed,min

(25)

式中:Eu為上水庫儲存能量,MW·h;Eu,min為上水庫儲存能量最小值,MW·h;Eu,max為上水庫儲存能量最大值,MW·h;Ed為下水庫儲存能量,MW·h;Ed,min為下水庫儲存能量最小值,MW·h;Ed,max為下水庫儲存能量最大值,MW·h。

電化學(xué)儲能系統(tǒng)約束為

Pcha,min≤Pcha(t)≤Pcha,max

(26)

Pdis,min≤Pdis(t)≤Pdis,max

(27)

式中:Pcha,max、Pcha,min分別為儲能電池充電功率上、下限,MW;Pdis,max、Pdis,min分別為儲能電池放電功率上、下限,MW;Pcha(t)、Pdis(t)分別為t時刻儲能電池充放電功率,MW。

考慮到儲能系統(tǒng)的壽命問題,在儲能系統(tǒng)的充放電過程中,儲能系統(tǒng)的電池荷電狀態(tài)(SOC)不能超過設(shè)定的上下限,即

RS,min≤RSOC≤RS,max

(28)

式中:RS,min為最小荷電狀態(tài),取0.2;RS,max為最大荷電狀態(tài),取0.9;RSOC為儲能系統(tǒng)的電池SOC。

系統(tǒng)外購電力約束為

0≤Pbuy(t)≤Pbuy,max

(29)

式中:Pbuy,max為系統(tǒng)最大外購電量,MW。

功率平衡約束為

(30)

式中:Pload(t)為t時刻系統(tǒng)負(fù)荷,MW。

3 算例分析

采用MATLAB Yalmip工具箱的Cplex優(yōu)化求解器進行優(yōu)化調(diào)度求解。選取24 h為1個周期,電力負(fù)荷及風(fēng)光出力數(shù)據(jù)參考我國西北某區(qū)域電網(wǎng)夏季典型日的實際數(shù)據(jù)。其中,火電裝機3 000 MW、風(fēng)電裝機1 500 MW、光伏裝機1 000 MW、抽水蓄能裝機1 000 MW,電化學(xué)儲能的功率/容量為200 MW/500 MW·h。根據(jù)2021年、2022年度全國碳排放權(quán)交易配額總量設(shè)定與分配實施方案,各參數(shù)取值為:Be=0.815 9,Fr=0.934,F1=1,Ff=1.087,Bh=0.110 4,ηe=0.728,βe=0.8,βs=1.08,c=0.25,α=0.25,d=5 000,a=0.3。該區(qū)域的電力負(fù)荷及風(fēng)光出力數(shù)據(jù)如圖1所示。

圖1 負(fù)荷出力曲線

系統(tǒng)外購電力價格按照當(dāng)?shù)叵募痉謺r電價并且執(zhí)行尖峰電價機制,分時電價曲線如圖2所示。火電上網(wǎng)電價為0.25元/(kW·h),風(fēng)光上網(wǎng)電價為0.262元/(kW·h)。電化學(xué)儲能充放電價按照分時電價執(zhí)行,根據(jù)最新的儲能政策,補貼設(shè)置為:棄風(fēng)棄光段(00:00:00—05:00:00、13:00:00—18:00:00)充電補貼0.55元/(kW·h),放電補貼0.25元/(kW·h);平段充放電均為0.2元/(kW·h)。抽水蓄能上網(wǎng)電價與燃煤機組上網(wǎng)電價相等,抽水電價為燃煤機組上網(wǎng)電價的75%?？紤]折現(xiàn)率及機組壽命,計算得出各類型能源的LCOE為:燃煤機組的單位度電成本為0.15元/(kW·h)、光伏機組的單位度電成本為0.18元/(kW·h)、風(fēng)電機組的單位度電成本為0.16元/(kW·h)、抽水蓄能系統(tǒng)的單位度電成本為0.25元/(kW·h)、電化學(xué)儲能系統(tǒng)的單位度電成本為0.65元/(kW·h)。

圖2 夏季電網(wǎng)分時電價

為比較不同碳交易機制不同儲能參與下的系統(tǒng)優(yōu)化調(diào)度情況以及碳排放情況,提出6種場景,具體場景情況見表1。

表1 場景設(shè)置

按照所構(gòu)建的6種場景及所給定的參數(shù),以系統(tǒng)總成本最小為目標(biāo),進行調(diào)度優(yōu)化計算,得出的優(yōu)化調(diào)度結(jié)果如圖3～圖8所示。

圖3 場景一

圖4 場景二

圖5 場景三

圖6 場景四

圖7 場景五

圖8 場景六

不同場景下的碳排放情況見表2。比較不同碳交易機制下的碳排放情況,場景四相較于場景一,碳排放量減少1 502.65 t;場景五相較于場景二,碳排放量減少1 288.47 t;場景六相較于場景三,碳排放量減少1 634.68 t。可以看出,階梯式碳交易機制下,其碳排放量小于傳統(tǒng)碳交易機制下的碳排放量,證明了階梯式碳交易的優(yōu)勢。比較不同儲能配置方法調(diào)度優(yōu)化后的碳排放情況,可以看出抽水蓄能與電化學(xué)儲能協(xié)同調(diào)度下的碳排放量低于抽水蓄能和電化學(xué)儲能單獨調(diào)度下的碳排放量,證明了電化學(xué)儲能與抽水蓄能協(xié)同調(diào)度下參與碳市場的可行性。

表2 不同場景下的碳排放情況

為驗證不同碳交易機制下不同儲能配置方式的經(jīng)濟效益優(yōu)勢,計算其電力市場以及碳市場的收益情況,結(jié)果匯總于表3。其中,電力市場收益的計算公式為:負(fù)荷×分時電價-購能成本-LCOE。由于本系統(tǒng)新能源比例較高,因此最終的碳排放配額高于實際碳排放量,二者之差即為可以進行售賣的碳配額量,其收益計算公式參考式(11)和式(12)。

由表3可以得出,傳統(tǒng)碳市場交易機制下的售碳收益低于階梯式碳市場交易機制下的售碳收益,這將更有利于鼓勵高碳排企業(yè)進行低碳轉(zhuǎn)型,從而售賣多余碳排放配額獲得更高收益;抽水蓄能與電化學(xué)儲能協(xié)同調(diào)度下的電力市場收益和碳市場收益均高于抽水蓄能和電化學(xué)儲能單獨調(diào)度下的電力市場收益和碳市場收益,進一步證明了階梯式碳交易機制下電化學(xué)儲能與抽水蓄能協(xié)同調(diào)度下參與電力市場及碳市場的可行性。

4 結(jié)論

(1) 各時段各類型能源的配置情況為:抽水蓄能主要通過谷時蓄能、峰時放能來實現(xiàn)對電網(wǎng)的調(diào)峰,并且獲得收益;電化學(xué)儲能主要通過峰谷分時電價實現(xiàn)高低電價間套利,以及通過容量成本回收機制獲得補償收益。二者協(xié)同發(fā)展的經(jīng)營模式,有利于維護電網(wǎng)的穩(wěn)定,并且提高系統(tǒng)整體收益。

(2) 傳統(tǒng)碳排放機制下,系統(tǒng)的整體碳減排積極性難以調(diào)動,引入階梯式碳交易機制后,系統(tǒng)的碳減排積極性有明顯的提高。傳統(tǒng)碳市場交易機制下的售碳收益低于階梯式碳市場交易機制下的售碳收益,這將更有利于鼓勵高碳排企業(yè)進行低碳轉(zhuǎn)型。

(3) 抽水蓄能與電化學(xué)儲能協(xié)同調(diào)度優(yōu)化下的碳排放量低于抽水蓄能和電化學(xué)儲能單獨調(diào)度下的碳排放量,證明了抽水蓄能與電化學(xué)儲能協(xié)同調(diào)度下參與碳市場的可行性;抽水蓄能與電化學(xué)儲能協(xié)同調(diào)度下的電力市場收益高于抽水蓄能和電化學(xué)儲能單獨調(diào)度下的電力市場收益,證明了電化學(xué)儲能與抽水蓄能協(xié)同調(diào)度下參與電力市場的可行性。所提出的電化學(xué)儲能與抽水蓄能協(xié)同參與電力市場與碳市場的配置方案,驗證了電化學(xué)儲能與抽水蓄能同時接入電網(wǎng)后進行系統(tǒng)調(diào)度優(yōu)化的可行性。