考慮時變電碳因子的園區(qū)綜合能源系統(tǒng)低碳經(jīng)濟調(diào)度

翁格平，任嬌蓉，姚 艷，王 亮，黃 芳

（國網(wǎng)浙江省電力有限公司寧波供電公司，浙江 寧波 315100）

0 引言

在全球環(huán)境變暖及能源危機日益嚴重的背景下，我國于2020年9月提出“2030碳達峰、2060碳中和”的雙碳目標［1-2］。目前我國能源結(jié)構(gòu)仍以化石能源為主導，風電光伏等新能源占比較小，能源行業(yè)碳排放占我國總體碳排放量比重較大。因此充分挖掘能源行業(yè)的碳減排能力，對于我國雙碳目標的實現(xiàn)具有重要意義。

具有多能耦合互補、能源轉(zhuǎn)換形式多樣等優(yōu)勢的綜合能源系統(tǒng)是能源行業(yè)低碳經(jīng)濟運行的重要解決方案之一［3-5］?，F(xiàn)有綜合能源系統(tǒng)低碳經(jīng)濟運行的研究，大多是從提高新能源消納水平［6-7］、高碳機組低碳化改造［8］、碳交易市場與能量市場的交互［9-10］等角度出發(fā)的。

從提高新能源消納水平角度出發(fā)，降低火電機組與高碳機組出力，從而可達到低碳運行的目標。文獻［11］以碳排放最少、風電消納水平最高、系統(tǒng)運行成本最低為目標，采用蓄熱式電采暖消納富余風電，有效降低了系統(tǒng)碳排放總量。文獻［12］考慮了機組層面碳排放約束對綜合能源系統(tǒng)運行優(yōu)化，研究結(jié)果表明考慮碳排放約束與需求響應能夠在一定程度上改善風電棄風現(xiàn)象，并有效降低系統(tǒng)運行過程中產(chǎn)生的二氧化碳。文獻［13］在某區(qū)域綜合能源系統(tǒng)中引入地源熱泵與需求響應，對綜合能源系統(tǒng)運行成本及環(huán)境成本等目標進行優(yōu)化，在提高風電消納水平的同時，降低了系統(tǒng)二氧化碳、二氧化硫等污染物排放量。

上述研究中，機組層面的碳排放特性被建模為定值［14］，實際上機組碳排放量會因機組運行工況、能源占比的不同而變化。所采用的度電碳排放因子應是時變曲線，反應能源占比的變化以及運行工況的變化［15］。此外，上述研究只考慮了常規(guī)發(fā)電機組如火電、熱電聯(lián)產(chǎn)機組在發(fā)電過程中的碳排放，沒有考慮到聯(lián)絡(luò)線功能及儲能裝置在充放電過程中的等效碳排放。

從高碳機組低碳化改造的角度出發(fā)，對高碳機組以及火電機組裝設(shè)碳捕集裝置，減少直接排放到空氣中的二氧化碳量。文獻［16］構(gòu)建靈活捕獲運行模式調(diào)節(jié)的碳捕集電廠與需求響應協(xié)同優(yōu)化模型，利用碳捕集設(shè)備捕碳能力的同時，采用需求響應改善了用能曲線。文獻［17］采用碳捕集技術(shù)中的富養(yǎng)燃燒技術(shù)對燃煤機組進行低碳化改造，仿真結(jié)果表明，富氧燃燒技術(shù)減碳能力優(yōu)于傳統(tǒng)燃燒后捕集的碳捕集電廠。文獻［18］提出一種計及儲熱裝置與碳捕集電廠的電力系統(tǒng)優(yōu)化調(diào)度模型，仿真結(jié)果表明碳捕集裝置在大幅降低碳排放的同時犧牲了一定的經(jīng)濟性。

上述研究需要在機組層面進行設(shè)備改造，但碳捕集設(shè)備的投資改造及維修費用較大，其在運行過程中的能耗也相對較高，目前經(jīng)濟性較差。

從碳交易市場與能量市場交互的角度出發(fā)，可利用固定碳交易價格、階梯型碳交易機制來約束高碳機組出力，實現(xiàn)系統(tǒng)低碳運行的目標。文獻［19］構(gòu)建了碳交易機制下的多主體均衡競價策略模型，在充分調(diào)動供能主體供能可靠性與積極性的同時，控制各供能主體的碳排放水平。文獻［20］分析電力、天然氣與碳排放市場的交互作用，通過市場手段調(diào)控不同能源交易方式，從而達到碳減排效果。文獻［21］為應對城市碳排放量與廢物處理的需求，提出了基于階梯碳交易機制的含廢物處理的城市綜合能源系統(tǒng)低碳經(jīng)濟模型。

上述研究從市場角度出發(fā)對源側(cè)機組層面的碳排放特性進行優(yōu)化，并沒有考慮到荷側(cè)可調(diào)度資源的低碳特性。荷側(cè)可調(diào)度資源在電碳因子曲線的引導下同樣可以響應調(diào)度需求，能夠?qū)崿F(xiàn)綜合能源系統(tǒng)低碳與經(jīng)濟運行的協(xié)調(diào)優(yōu)化［22］。

綜上，本文提出一種考慮時變電碳因子的園區(qū)綜合能源系統(tǒng)低碳經(jīng)濟調(diào)度模型。首先，在時變電碳因子曲線的基礎(chǔ)上，對源側(cè)機組、儲能裝置、聯(lián)絡(luò)線碳排特性進行建模；荷側(cè)利用可調(diào)度資源的削減轉(zhuǎn)移用能，構(gòu)建低碳環(huán)境下需求響應模型。最后，構(gòu)建考慮時變電碳因子的綜合能源系統(tǒng)低碳經(jīng)濟模型，在減少碳排放的同時，提高園區(qū)運行經(jīng)濟性。

1 低碳環(huán)境下綜合能源系統(tǒng)架構(gòu)

本節(jié)以浙江地區(qū)某園區(qū)型綜合能源系統(tǒng)為例進行分析，系統(tǒng)架構(gòu)如圖1 所示，其中CHP 為熱電聯(lián)產(chǎn)機組。能源供給機組包括分布式風電、分布式光伏、聯(lián)絡(luò)線輸送功率以及本地火電機組；多能耦合機組包括電轉(zhuǎn)氣、電制冷機、熱電聯(lián)供等機組；負荷側(cè)以園區(qū)內(nèi)部企業(yè)負荷為主。碳能流線條的粗細代表了能源在生產(chǎn)轉(zhuǎn)換過程中產(chǎn)生的二氧化碳的數(shù)量多少。

1.1 電碳因子曲線

目前電力系統(tǒng)全環(huán)節(jié)碳排放計量與分析的現(xiàn)狀為主要關(guān)注與發(fā)電側(cè)碳計量，且是簡單的非時變的固定碳排放系數(shù)。但發(fā)電側(cè)各類發(fā)電機組的實際運行工況并非時時處于額定運行狀態(tài)，因此每一時段產(chǎn)生的二氧化碳量應是時變曲線，傳統(tǒng)的平均電碳因子無法滿足降碳新要求，且無法體現(xiàn)時空特性。同時，上級電網(wǎng)聯(lián)絡(luò)線功率中的碳計量應考慮在低碳模型中。園區(qū)時變電碳因子則是在聯(lián)絡(luò)線中火電占比及相應碳排放數(shù)據(jù)的基礎(chǔ)上，進而耦合園區(qū)內(nèi)發(fā)電機組出力得到的。

園區(qū)時變電碳因子數(shù)值的高低表征火電出力占比，當火電占比較高時，產(chǎn)生較多的二氧化碳，對應電碳因子數(shù)值較高。因此，本文構(gòu)建了考慮聯(lián)絡(luò)線碳排放的園區(qū)內(nèi)時變電碳因子曲線模型如圖2所示。

圖2 時變電碳因子曲線

1.2 計及碳排放特性的熱電聯(lián)供建模

本文所采用的熱電聯(lián)供由燃氣輪機及余熱鍋爐組成，計及其碳排放特性的熱電聯(lián)供模型可表示為：

式中：FCHP，CO2（t）為熱電聯(lián)產(chǎn)機組碳排放量；PGT，e（t）和PGT，h（t）分別為t時段燃氣輪機電、熱出力；PCHP，h（t）為t時段熱電聯(lián)產(chǎn)機組熱出力；ηGT、ηloss、ηrec、ηh分別為燃氣輪機發(fā)電效率、熱損失參數(shù)、余熱鍋爐熱回收效率、產(chǎn)熱效率；為t時段熱電聯(lián)產(chǎn)機組等效單位碳排放數(shù)據(jù)。

1.3 電池儲能裝置碳排放模型

電池儲能裝置的能量流動如圖3所示。本文認為電池儲能裝置在充放電過程中會產(chǎn)生一定的二氧化碳，且由分布式風電及光伏給儲能充電時，由于新能源的零碳、清潔能源特性并不會產(chǎn)生二氧化碳的排放，因此儲能裝置的碳排放特性僅需要計及聯(lián)絡(luò)線充能及熱電聯(lián)產(chǎn)機組充能帶來的等效碳排放量。電池儲能裝置碳排放模型如式（5）—（7）所示［23］。

圖3 電池儲能充放電碳流示意圖

式中：Sess（t）和Sess（t-1）為t與t-1時段電池儲能設(shè)備的能量狀態(tài)；ηess為電池儲能裝置的自損率；Pess，cha（t-1）、Pess，dis（t-1）、ηess，cha、ηess，dis為t-1 時段電池儲能裝置的充、放電功率及充、放電效率；為儲能裝置的等效碳排放量；Pwvess（t）為t時段由新能源供給電池儲能的充電功率；A和B為儲能充、放電標志，充電時A=1，放電時B=1。

2 綜合需求響應低碳特性

2.1 激勵型需求響應低碳特性

激勵型需求響應通過對負荷側(cè)可調(diào)度資源進行調(diào)度的同時，對調(diào)度量進行補償，從而引導用戶改變用能曲線，響應系統(tǒng)調(diào)度需求。低碳環(huán)境下的激勵型需求響應具有較大的低碳潛力。本文擬對園區(qū)內(nèi)部高耗能產(chǎn)業(yè)進行激勵型需求響應。

在高電碳因子時段發(fā)電機組進行發(fā)電時產(chǎn)生的二氧化碳遠高于低電碳因子時段的二氧化碳排放量，因此可通過價格補償手段激勵用戶削減用能，在挖掘荷側(cè)可調(diào)度資源低碳屬性的同時，實現(xiàn)系統(tǒng)經(jīng)濟運行。激勵型需求響應低碳特性建模如下：

式中：Pcut（t）為t時段可削減負荷削減量；Pcut，max（t）和Pcut，min（t）分別為t時段可削減負荷削減量的上、下限；xt為t時段可削減負荷是否削減的0-1變量；FIDR，cut為可削減負荷調(diào)用費用；δcut為單位可削減負荷激勵補償價格。

2.2 價格型需求響應低碳特性

激勵型需求響應是在價格補償引導下響應低碳運行指標，而價格型需求響應低碳特性則是在分時電價以及時變電碳因子的雙重激勵引導下進行的響應系統(tǒng)低碳運行項目，雙重激勵引導作用如圖4所示。本文擬對園區(qū)內(nèi)部高耗能產(chǎn)業(yè)以及汽車制造業(yè)采取價格型需求響應項目。

圖4 雙重激勵下價格型需求響應特性

負荷側(cè)部分可調(diào)度資源具有轉(zhuǎn)移用能潛力，既可以在高電碳因子時段減少用能、在低電碳因子時段增加用能，也可以在價格信號引導下在電價較高時段減少用能、在電價較低時段增加用能。

式中：Fmov，co2為可轉(zhuǎn)移負荷需求響應二氧化碳減少量；FIDR，mov為可轉(zhuǎn)移負荷調(diào)用費用；δmov為單位可轉(zhuǎn)移負荷調(diào)用價格；Pmov（t）為t時段可轉(zhuǎn)移負荷轉(zhuǎn)移量；ΔPmov，max（t）和ΔPmov，min（t）為t時段可轉(zhuǎn)移負荷轉(zhuǎn)移量上、下限；yt為t時段可轉(zhuǎn)移負荷是否轉(zhuǎn)移的0-1 變量；可轉(zhuǎn)移負荷需要保證在調(diào)度周期內(nèi)用戶的用能總量保持不變。

3 低碳經(jīng)濟調(diào)度模型

3.1 目標函數(shù)

本文所提綜合能源系統(tǒng)低碳經(jīng)濟調(diào)度模型是基于時變電碳因子的以系統(tǒng)運行總成本和二氧化碳排放最少的多目標最優(yōu)模型，運行總成本包括機組運維成本、購能成本以及需求響應調(diào)用費用。

式中：FIES為園區(qū)綜合能源系統(tǒng)總目標函數(shù)；Feco和為經(jīng)濟性指標與低碳性指標；α1和α2分別為經(jīng)濟性指標與低碳指標對應的權(quán)重系數(shù)，由模糊層次分析法得到［24］。

式中：FIDR、Fbuy、Fmc分別為園區(qū)綜合能源系統(tǒng)調(diào)度周期內(nèi)的需求響應調(diào)用費用、購能成本、運維費用；Pbuyg（t）為t時段購氣量；δg和δe分別為相應的購氣、購電單價；Pi（t）和δi為第i類機組出力和運維單價。

3.2 約束條件

1）功率平衡約束

式中：Pel（t）、Pgl（t）、Phl（t）、Pcol（t）分別為t時段電、氣、熱、冷負荷；Phd（t）、Pwl（t）、Pvl（t）、PCHP，el（t）分別為t時段火電出力、風電供給負荷量、光伏供給負荷量、熱電聯(lián)產(chǎn)供給負荷量；Pe，co（t）為t時段電制冷機出力；Pbuygl（t）和Pp2g（t）分別為t時段購氣供給負荷量、電轉(zhuǎn)氣裝置出力；ηp2g為電轉(zhuǎn)氣運行效率。

2）聯(lián)絡(luò)線購能約束

式中：Psell（t）為t時段售電功率；Pbuyg，max（t）、Pbuyg，min（t）、Pbuy，max（t）、Pbuy，min（t）、Psell，max（t）、Psell，min（t）分別為t時段購氣功率上、下限，購電功率上、下限以及售電功率上、下限。

3）多能耦合機組運行約束

式中：Pw（t）、Pv（t）、PCHP，e（t）分別為t時段風電、光伏、熱電聯(lián)產(chǎn)機組電出力；Phd，max（t）、Phd，min（t）、Pw，max（t）、Pw，min（t）、Pv，max（t）、Pv，min（t）、PCHPe，max（t）、PCHPe，min（t）、Pp2g，max（t）、Pp2g，min（t）、Sess，max（t）、Sess，min（t）分別為t時段火電機組、風電、光伏、熱電聯(lián)產(chǎn)機組、電轉(zhuǎn)氣裝置出力的上限和下限，儲能荷電狀態(tài)上、下限。

4）荷側(cè)可調(diào)度資源響應約束

荷側(cè)可調(diào)度資源在響應系統(tǒng)調(diào)度指令時，需要滿足一定的約束條件，如調(diào)度次數(shù)限制、調(diào)度時間段限制等?？烧{(diào)度資源轉(zhuǎn)移時間約束為：

式中：ymax和ymin分別為轉(zhuǎn)移時間上、下限約束；λmov為t時段可調(diào)度資源占比，由負荷側(cè)上傳可調(diào)度資源占比。

可調(diào)度資源削減時間約束為：

式中：xmax和xmin分別為可削減負荷削減次數(shù)上、下限約束。

4 算例分析

4.1 算例系統(tǒng)架構(gòu)及參數(shù)

本文以浙江地區(qū)某園區(qū)型綜合能源系統(tǒng)為研究對象，其架構(gòu)如圖1所示。園區(qū)內(nèi)部多能負荷預測曲線、分布式風光出力預測曲線及能源價格如圖5—7 所示。多能耦合機組裝機容量等參數(shù)如表1—2所示。

表1 多能耦合機組參數(shù)

圖5 園區(qū)風光出力及負荷曲線

圖6 負荷類型及用電量曲線

圖7 分時電價曲線

本文以調(diào)度周期24 h 進行仿真，調(diào)度時間步長為1 h。系統(tǒng)內(nèi)部電價采用浙江某地現(xiàn)有分時電價進行仿真。

為驗證本文所提低碳經(jīng)濟調(diào)度模型的有效性，設(shè)置如下場景進行仿真驗證，其中場景3為本文所提出的考慮時變電碳因子的源荷儲協(xié)調(diào)綜合能源系統(tǒng)低碳經(jīng)濟調(diào)度模型。

表2 電池儲能裝置參數(shù)

場景1 為基礎(chǔ)場景，不涉及儲能裝置及需求響應。

場景2為儲能場景，裝設(shè)儲能裝置且計及儲能充放電與聯(lián)絡(luò)線功率的碳排放。

場景3為本文所提考慮時變電碳因子的園區(qū)綜合能源系統(tǒng)低碳經(jīng)濟運行場景，在時變電碳因子曲線、價格信號引導下，源側(cè)機組、儲能、荷側(cè)可調(diào)度資源協(xié)調(diào)配合。

4.2 場景調(diào)度結(jié)果分析

本文所設(shè)定的各場景調(diào)度結(jié)果如表3中數(shù)據(jù)所示。場景1為基礎(chǔ)場景，由于未考慮儲能裝置的能量時移作用以及需求響應，系統(tǒng)運行經(jīng)濟性較差，聯(lián)絡(luò)線購能功率較多導致碳排放量較多。場景2中裝設(shè)了儲能裝置，計及儲能裝置能量時移特性及碳排放特性，系統(tǒng)運行總成本相較于場景1 下降1.06%，二氧化碳排放水平相較于場景1 下降2.81%，減少運行成本的同時降低碳排放，驗證了儲能裝置在碳排放特性上的低碳性。場景3為本文所提出的考慮時變電碳因子的園區(qū)綜合能源系統(tǒng)低碳經(jīng)濟調(diào)度模型，通過合理控制儲能充放電以及荷側(cè)可調(diào)度資源的調(diào)用計劃，使得系統(tǒng)運行總成本下降6.61%，而系統(tǒng)碳排放水平下降6.05%，源荷儲協(xié)調(diào)配合，提高運行經(jīng)濟性、降低系統(tǒng)運行過程中的碳排放水平。

表3 低碳經(jīng)濟調(diào)度數(shù)據(jù)

進一步分析3個場景的聯(lián)絡(luò)線交互功率與對應的碳排放水平，如圖8所示。場景3相較于場景1，高電碳因子時段、高用能電價時段的聯(lián)絡(luò)線交互功率下降19.9%，通過需求響應使得用電高峰、電碳因子高時段用電負荷下降，說明用戶在電碳因子以及分時電價雙重激勵下的引導作用降低了系統(tǒng)碳排放水平及運行成本。

圖8 場景聯(lián)絡(luò)線功率對比

場景1的各機組出力調(diào)度結(jié)果如圖9所示。由圖9可知本園區(qū)型綜合能源系統(tǒng)供能主要與外部聯(lián)絡(luò)線供能，而園區(qū)內(nèi)分布式風電與光伏得到全額消納。熱電聯(lián)產(chǎn)機組因需承擔園區(qū)內(nèi)熱負荷需求，全天處于工作狀態(tài)，且因其“以熱定電”運行約束，熱電聯(lián)產(chǎn)機組出力相對固定隨著用熱負荷變動。場景1 聯(lián)絡(luò)線功率占全天總用電負荷的53.14%，且在17—18 時段、20—22 時段聯(lián)絡(luò)線功率較高，因此可見場景1整體碳排放水平與運行成本較高。

圖9 場景1機組出力

場景2為儲能場景且計及儲能充放電過程中的碳排放。場景2 調(diào)度結(jié)果如圖10 所示。由于本文裝設(shè)的儲能裝置容量較小，儲能裝置僅在18、20—21 時段對用電負荷進行放電，此外，儲能裝置還在14—16 時段、20—21 時段對電轉(zhuǎn)氣裝置放電。將電價較低時段儲存電能釋放出來，在起到能量時移作用的同時降低高電碳因子曲線時段的外購電，降低園區(qū)二氧化碳排放水平。高電碳因子時段園區(qū)外購電功率相較于場景1下降1.75%。

圖10 場景2機組出力

場景3為本文所提考慮時變電碳因子的綜合能源系統(tǒng)低碳經(jīng)濟調(diào)度模型，調(diào)度結(jié)果如圖11所示，荷側(cè)可調(diào)度資源調(diào)用計劃如圖12所示。由圖11可知，通過需求響應對于總用電負荷曲線的改善作用，需求響應后峰谷差下降48%的同時，緩解了用電高峰供需緊張關(guān)系，減少了高電碳因子時段用能，促進園區(qū)低碳經(jīng)濟運行。高電碳因子時段出力相較于場景2 進一步下降了19.37%，而總體二氧化碳排放水平相應地下降了3.33%。

圖11 場景3機組出力

由圖12 可知，高耗能企業(yè)通過用能削減與用能轉(zhuǎn)移，削減了午高峰及晚高峰用電量，在緩解晚高峰用電矛盾的同時，降低了自身企業(yè)用電成本，同時響應系統(tǒng)調(diào)度需求，達到低碳運行的目的。晚高峰削減負荷量占晚高峰電負荷總量的2.82%。此外，汽車制造業(yè)由于具有負荷可轉(zhuǎn)移特性且調(diào)度時間較自由，因此在電價低谷時段、低電碳因子時段增加用能，而在電價平時段、電價高峰時段、高電碳因子時段減少用能。晚高峰時段汽車制造業(yè)轉(zhuǎn)移用能量占晚高峰電負荷總量的6.46%。而高新技術(shù)產(chǎn)業(yè)作為基礎(chǔ)負荷，可調(diào)度特性較差，因此不對高新技術(shù)產(chǎn)業(yè)進行調(diào)度。

圖12 荷側(cè)可調(diào)度資源調(diào)用計劃

4.3 多目標優(yōu)化調(diào)度模型分析

本文所設(shè)置需求響應場景為低碳性與經(jīng)濟性均衡下的最優(yōu)解，現(xiàn)研究不同權(quán)重下低碳性與經(jīng)濟性的變化情況。不同權(quán)重下的低碳指標與經(jīng)濟指標參數(shù)如圖13所示。由圖13可知，隨著低碳指標權(quán)重的降低，系統(tǒng)碳排放水平處于先增加后減少的趨勢，而經(jīng)濟性指標則隨著權(quán)重系數(shù)的增加處于先減少后增加的趨勢，最終在本文均衡解處取得多目標優(yōu)化調(diào)度的最優(yōu)解。

圖13 多目標調(diào)度結(jié)果

因為本文外購電水平較高，各個場景外購電均占總用電負荷的50%以上，因此碳排放水平主要取決于外購電占比，低碳指標與經(jīng)濟性指標具有相同的數(shù)據(jù)表征形式，即通過需求響應與儲能裝置的合理充放電使得高電碳因子時段聯(lián)絡(luò)線購能功率減少，低電碳因子時段外聯(lián)絡(luò)線購能功率增加，多余電能通過儲能裝置達到能量時移、低碳經(jīng)濟的作用。因此本質(zhì)上低碳與經(jīng)濟指標具有相同的變化趨勢。

5 結(jié)論

在我國“雙碳”目標背景下，針對能源行業(yè)碳排放水平較高問題，本文基于時變電碳因子構(gòu)建了園區(qū)綜合能源系統(tǒng)低碳經(jīng)濟多目標優(yōu)化調(diào)度模型，通過算例分析可得到如下結(jié)論：

1）構(gòu)建考慮聯(lián)絡(luò)線碳排放水平、儲能裝置碳排放特性的電力行業(yè)全環(huán)節(jié)碳流排放模型進行低碳經(jīng)濟調(diào)度，能夠有效優(yōu)化負荷曲線、降低碳排。

2）在調(diào)度過程中，充分挖掘荷側(cè)可調(diào)度資源的低碳特性，在高電碳因子時段削減轉(zhuǎn)移用能，使得園區(qū)運行總成本下降6.61%，而二氧化碳排放水平下降6.05%，驗證了本文源荷儲協(xié)調(diào)優(yōu)化模型在降低碳排、減少運行成本的有效性。

3）本文所提出的綜合能源系統(tǒng)低碳經(jīng)濟調(diào)度模型能夠兼顧系統(tǒng)兩個方面的優(yōu)化目標，實現(xiàn)低碳性與經(jīng)濟性的統(tǒng)一。