李顥然，薛屹洵*，戴鐵潮，常馨月，潘昭光，孫宏斌,

(1.太原理工大學(xué) 電氣與動(dòng)力工程學(xué)院，山西 太原 030002；2.國網(wǎng)浙江省電力有限公司，浙江 杭州 310007；3.清華大學(xué)電機(jī)工程與應(yīng)用電子技術(shù)系，北京 100084)

近年來，隨著新能源的大規(guī)模開發(fā)利用，各地都出現(xiàn)了不同程度的棄風(fēng)棄光現(xiàn)象，給社會(huì)帶來一定的經(jīng)濟(jì)損失[1–2]。電解水制氫作為消納新能源電力的一種有效途徑，已被廣泛應(yīng)用于中國的一些城市與工業(yè)園區(qū)[3]，因此出現(xiàn)了越來越多的電–氫耦合系統(tǒng)。

目前，國內(nèi)外學(xué)者已經(jīng)對(duì)電–氫耦合系統(tǒng)展開了一系列研究。一些研究聚焦于提升電–氫耦合系統(tǒng)的靈活性，以適應(yīng)高比例可再生能源接入的電網(wǎng)。Niaz等[4]提出將堿性電解槽與電池儲(chǔ)能系統(tǒng)集成，以克服可再生能源出力的波動(dòng)性；Ban等[5]提出了一個(gè)同時(shí)具有電解槽和燃?xì)獍l(fā)電設(shè)施的能源中心，并驗(yàn)證了該運(yùn)行模式的可行性及在風(fēng)電消納上的有效性。在電–氫耦合系統(tǒng)的運(yùn)行調(diào)度方面，范宏等[6]建立了一種氫儲(chǔ)能系統(tǒng)的動(dòng)態(tài)特性模型，并將其運(yùn)用到智慧園區(qū)多樓宇系統(tǒng)的經(jīng)濟(jì)調(diào)度中，以此降低多樓宇系統(tǒng)的運(yùn)行成本；顧玖等[7]提出構(gòu)建一種電–氫一體化站，通過加入一體化站的運(yùn)行模式，提高了電力系統(tǒng)削峰填谷能力及新能源消納能力；Amin等[8]提出了一種基于熱電聯(lián)產(chǎn)與儲(chǔ)氫融合的風(fēng)能電網(wǎng)日前調(diào)度方案，并通過仿真驗(yàn)證了所提出方案的有效性；Zhong等[9]提出了一種含有電制氫裝置的微電網(wǎng)實(shí)時(shí)優(yōu)化模型，并采用預(yù)測(cè)控制和短視策略解決優(yōu)化問題；Jun等[10]提出了一個(gè)考慮氫氣銷售的電–氫耦合系統(tǒng)模型，綜合考慮儲(chǔ)氫、售氫及區(qū)域電力系統(tǒng)、供熱系統(tǒng)進(jìn)行優(yōu)化調(diào)度以降低系統(tǒng)的運(yùn)行成本。在電–氫協(xié)同規(guī)劃方面，許傳博等[11]構(gòu)建了一種儲(chǔ)氫優(yōu)化配置模型，通過配置發(fā)電側(cè)新能源站場(chǎng)的儲(chǔ)氫容量，解決風(fēng)光站場(chǎng)的能源消納問題；李奇等[12]提出一種考慮電–氫耦合的混合儲(chǔ)能微電網(wǎng)容量?jī)?yōu)化配置方法，建立微電網(wǎng)優(yōu)化配置模型，并提出該微電網(wǎng)運(yùn)行控制策略；Wang等[13]提出一種輸電與輸氫網(wǎng)絡(luò)聯(lián)合規(guī)劃方法，量化輸氫系統(tǒng)的靈活性，實(shí)現(xiàn)電氫資源的互補(bǔ)。

然而，上述研究?jī)H在樓宇、城市等中小規(guī)模用氫背景下考慮，并不適用于工業(yè)園區(qū)中大規(guī)模用氫的生產(chǎn)環(huán)境，難以保證在工業(yè)背景下的成本問題；同時(shí)，上述研究?jī)H考慮了氫氣作為電能儲(chǔ)存媒介的作用，而忽略了氫氣在化工行業(yè)的重要用途。事實(shí)上，氫氣作為一種重要的化工原料，在合成氨[14]、甲醇等化工領(lǐng)域有著十分重要的作用，以氫氣為主要原料的氫化工園區(qū)也在中國北方廣泛存在，同時(shí)化工行業(yè)的負(fù)荷特性決定其具有參與需求側(cè)響應(yīng)的潛力[15]，可以為電力系統(tǒng)提供可調(diào)的靈活性資源。然而，得益于成熟的技術(shù)條件及較低的制氫成本[16]，當(dāng)前大部分氫氣仍來自于化石燃料制氫，如天然氣制氫和煤氣化制氫[17]，這種傳統(tǒng)的化石燃料制氫方式往往伴隨著嚴(yán)重的碳排放問題，不利于中國“雙碳”目標(biāo)的實(shí)現(xiàn)，亟需開展以化工園區(qū)為背景的電–氫耦合系統(tǒng)研究。

基于以上考慮，本文對(duì)化工園區(qū)場(chǎng)景下的電–氫耦合系統(tǒng)進(jìn)行了全面研究。建立了涵蓋園區(qū)源–網(wǎng)–荷–儲(chǔ)各角度的電–氫耦合系統(tǒng)各側(cè)設(shè)備的經(jīng)濟(jì)調(diào)度模型，在氫源側(cè)，將煤制氫與電解水制氫同時(shí)作為園區(qū)氫源，使電解水制氫在園區(qū)中既可以消納富余電力，又可以替代部分煤制氫，實(shí)現(xiàn)煤與可再生能源的有機(jī)協(xié)同；在氫網(wǎng)側(cè)，提出采用管道實(shí)現(xiàn)園區(qū)內(nèi)的短距離氫氣運(yùn)輸；在氫負(fù)荷側(cè)，充分挖掘氫燃料發(fā)電[18–19]與氫制氨的靈活性，使其作為可調(diào)度的資源。系統(tǒng)調(diào)度以實(shí)現(xiàn)園區(qū)收益最大化為目標(biāo)，并加入棄風(fēng)棄光懲罰以實(shí)現(xiàn)最大限度消納新能源，綜合考慮園區(qū)電力系統(tǒng)與氫系統(tǒng)的運(yùn)行約束條件，對(duì)系統(tǒng)優(yōu)化調(diào)度進(jìn)行仿真分析。

1 園區(qū)級(jí)電–氫耦合系統(tǒng)架構(gòu)

本文以實(shí)際綠氫制氨園區(qū)為基礎(chǔ)[20–21]，構(gòu)建了如圖1所示的化工園區(qū)電–氫耦合系統(tǒng)總體架構(gòu)。園區(qū)內(nèi)風(fēng)、光等可再生能源充足，同時(shí)園區(qū)電力系統(tǒng)通過變電站與外部電網(wǎng)進(jìn)行電力交易，氫源包括煤制氫與電解水制氫，氫氣通過儲(chǔ)氫罐進(jìn)行存儲(chǔ)，園區(qū)內(nèi)的氫氣使用管道進(jìn)行短距離運(yùn)輸，氫負(fù)荷包括加氫站、氫制氨，以及使用氫氣作為燃料進(jìn)行發(fā)電的氫燃料發(fā)電。

圖1 園區(qū)級(jí)電–氫耦合系統(tǒng)Fig. 1 Park-level electricity–hydrogen coupling system

整個(gè)園區(qū)通過電解水制氫與氫燃料發(fā)電使電力系統(tǒng)與氫系統(tǒng)緊密耦合在一起。在園區(qū)正常運(yùn)行時(shí)，通過外電網(wǎng)與分布式風(fēng)電、光伏對(duì)園區(qū)供電，滿足電負(fù)荷所需電能。氫氣通過煤制氫與電解水制氫制備，當(dāng)前已被廣泛使用的電解水制氫設(shè)備主要為堿性電解槽（alkaline water electrolyser，AWE）與質(zhì)子交換膜電解槽（proton exchange membrane electrolyser，PEME），而堿性電解槽在工業(yè)大規(guī)模量產(chǎn)上有更廣泛的應(yīng)用[22–23]，由氫源制備的氫氣先由儲(chǔ)氫罐進(jìn)行存儲(chǔ)，而后通過管道輸送至氫的各個(gè)負(fù)荷，當(dāng)前較為成熟的管道輸氫技術(shù)為天然氣混氫運(yùn)輸，出于安全考慮，一般將氫氣的體積比控制在5%以內(nèi)[24]。在負(fù)荷側(cè)，氫燃料發(fā)電與氫制氨通過跟隨電源與電負(fù)荷的實(shí)時(shí)變化情況，同時(shí)響應(yīng)分時(shí)電價(jià)，以達(dá)到園區(qū)整體運(yùn)行成本最低及收益的最大化的效果。

2 化工園區(qū)氫能系統(tǒng)模型

對(duì)化工園區(qū)氫系統(tǒng)關(guān)鍵設(shè)備進(jìn)行建模，包括煤制氫、電解槽、儲(chǔ)氫罐、氫制氨、氫燃料發(fā)電機(jī)，以及輸氫管網(wǎng)。

2.1 氫源模型

1）煤制氫

煤制氫是指煤在高溫、高壓下與水蒸氣反應(yīng)生成氫氣與二氧化碳的過程，其化學(xué)反應(yīng)原理如下：

2）電解水制氫

本文依據(jù)電極動(dòng)力學(xué)建立了堿性電解槽的電流–電壓經(jīng)驗(yàn)關(guān)系模型：

故電解槽電壓與電流即可簡(jiǎn)化為線性形式，根據(jù)電化學(xué)理論與法拉第定律，電解槽的實(shí)際產(chǎn)氫量可以表示為：

由于電解槽所用電為直流電，從電網(wǎng)接入的電能需通過整流裝置將交流電轉(zhuǎn)化為直流電，電的整流過程可表示為：

2.2 儲(chǔ)氫罐模型

在儲(chǔ)氫罐模型中，重點(diǎn)關(guān)注罐內(nèi)總氫氣質(zhì)量的

2.3 氫負(fù)荷模型

1）氫制氨

氫制氨是合成氨的重要來源之一，利用氫氣與氮?dú)庠诟邷馗邏合碌姆磻?yīng)生產(chǎn)氨氣，其化學(xué)反應(yīng)方程式如下：

在實(shí)際氫制氨工藝中原料往往要進(jìn)入多反應(yīng)容器，同時(shí)配有加壓、冷卻、分離等輔助系統(tǒng)以達(dá)到較高的轉(zhuǎn)化率，故氫制氨的生產(chǎn)過程往往也消耗大量電能。Li等[25]基于布勞恩軸向三床絕緣反應(yīng)器（Braun axial three-beds insulation reactor，BATIR），建立了氫制氨模型，證明反應(yīng)器中各變量，包括溫度、輸入電能、輸入氫氣、輸入氮?dú)?、輸出熱能、輸出氨氣之間存在一種線性關(guān)系。基于該模型，給出表征氫氣輸入、氮?dú)廨斎?、用電量及氨氣產(chǎn)量之間線性關(guān)系的方程：

2.4 輸氫管網(wǎng)模型

式（19）即為管道穩(wěn)態(tài)方程的線性形式，在實(shí)際管道運(yùn)行中，往往需考慮氣體流向，由此可將式（19）改寫為：

3 電–氫協(xié)同優(yōu)化調(diào)度模型

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

為實(shí)現(xiàn)新能源就地消納與園區(qū)運(yùn)行效益最大化的目標(biāo)，建立了如下的目標(biāo)函數(shù)：

3.2 約束條件

1）氫源

氫源約束除了應(yīng)滿足式（2）～（11）的模型，還應(yīng)滿足如下的運(yùn)行約束：

式（24）表示基于煤制氫容量規(guī)定了其用煤量的上限與下限，式（25）表示基于電解槽容量規(guī)定了電解槽輸入電能的上下限。

2）儲(chǔ)氫

儲(chǔ)氫約束除了應(yīng)滿足式（12）的模型，還應(yīng)滿足如下的運(yùn)行約束：

式（26）為儲(chǔ)氫罐的容量限制；式（27）表示充氫與放氫的速率限制；式（28）表示假設(shè)每天初始時(shí)刻有一個(gè)氫氣最低需求量，則在一天的初始時(shí)刻的儲(chǔ)氫量應(yīng)大于這一最低需求量。

3）氫負(fù)荷

氫負(fù)荷約束除了應(yīng)滿足式（14）～（15）的模型，還應(yīng)滿足如下的運(yùn)行約束：

4）輸氫管網(wǎng)

輸氫管網(wǎng)約束除了應(yīng)滿足式（16）～（22）的模型，還應(yīng)滿足如下運(yùn)行約束：

式（32）表示氫源出口一般通過壓縮機(jī)接入管道，因此可假設(shè)管道在氫源節(jié)點(diǎn)處的壓力為定值；式（33）表示管道承受壓力范圍。

5）電量平衡約束

園區(qū)電力系統(tǒng)應(yīng)滿足如下電量平衡約束：

4 算例分析

4.1 算例設(shè)置

本文所構(gòu)建的調(diào)度模型是含非線性約束的優(yōu)化問題，算例研究在16 GB內(nèi)存、i7–1165G7 CPU的計(jì)算環(huán)境下完成。仿真程序使用MATLAB R2021a開發(fā)并使用ipopt求解器進(jìn)行求解。

本文基于典型化工園區(qū)電氫耦合系統(tǒng)設(shè)計(jì)算例，算例電氫耦合系統(tǒng)拓?fù)浣Y(jié)構(gòu)如圖2所示，其中，P2H為電解水制氫、C2H為煤制氫、HFG為氫燃料發(fā)電、H2A為氫制氨。

圖2 算例電–氫耦合系統(tǒng)拓?fù)浣Y(jié)構(gòu)Fig. 2 Electricity–hydrogen coupling system topology in case studies

由圖2可知，園區(qū)系統(tǒng)由一個(gè)11節(jié)點(diǎn)的電力系統(tǒng)和一個(gè)7節(jié)點(diǎn)的氫能系統(tǒng)耦合而成。在氫能系統(tǒng)中節(jié)點(diǎn)2設(shè)有電解水制氫，節(jié)點(diǎn)3設(shè)有煤制氫，節(jié)點(diǎn)1、4設(shè)有儲(chǔ)氫，節(jié)點(diǎn)5設(shè)有氫制氨，節(jié)點(diǎn)7設(shè)有氫燃料發(fā)電機(jī)，節(jié)點(diǎn)6為加氫站用氫負(fù)荷，其中，電解水制氫、氫制氨與氫燃料發(fā)電機(jī)分別耦合在電力系統(tǒng)的節(jié)點(diǎn)3、6與11，風(fēng)電與光伏分別在電力系統(tǒng)的節(jié)點(diǎn)7和8，在節(jié)點(diǎn)1處通過變電站與外電網(wǎng)進(jìn)行電力交易。

分布式風(fēng)電、光伏出力數(shù)據(jù)、園區(qū)總電源、用戶用電負(fù)荷數(shù)據(jù)如圖3所示，園區(qū)中加氫站氫負(fù)荷數(shù)據(jù)如圖4所示。

圖3 園區(qū)電量數(shù)據(jù)Fig. 3 Park electricity data

圖4 加氫站用氫數(shù)據(jù)Fig. 4 Hydrogen data for hydrogen refueling station

電解槽參數(shù)采用典型堿性電解槽參數(shù)[26]，參數(shù)值如表1所示，氫制氨運(yùn)行參數(shù)如表2所示。

表1 堿性電解槽參數(shù)Tab. 1 Parameters of alkaline electrolyzer

表2 氫制氨參數(shù)Tab. 2 Parameters of hydrogen to ammonia

本文采用的電價(jià)為分時(shí)電價(jià)，具體如表3所示。

表3 分時(shí)電價(jià)數(shù)據(jù)Tab. 3 Time-of-use electricity price

4.2 算例結(jié)果分析

為驗(yàn)證所提方法的有效性，設(shè)置3個(gè)典型場(chǎng)景如下：

場(chǎng)景1：電解槽與氫制氨的運(yùn)行功率均為可調(diào)，即電解槽跟隨新能源出力，同時(shí)考慮氫制氨的需求側(cè)響應(yīng)；

場(chǎng)景2：電解槽無法跟蹤新能源出力，即電解槽處于定功率運(yùn)行狀態(tài)，氫制氨運(yùn)行功率可調(diào)；

場(chǎng)景3：不考慮氫制氨的需求側(cè)響應(yīng)，即氫制氨處于定功率運(yùn)行狀態(tài)，電解槽運(yùn)行功率可調(diào)。

化工園區(qū)電–氫耦合系統(tǒng)在3種場(chǎng)景下的運(yùn)行結(jié)果分別如圖5、6及7所示，包括電解水制氫的運(yùn)行曲線、氫燃料發(fā)電機(jī)的運(yùn)行曲線與氫制氨的運(yùn)行曲線。

圖5 場(chǎng)景1運(yùn)行結(jié)果Fig. 5 Operation results of case 1

從圖5可以看出：新能源出力的峰值位于上午的07：00—09：00時(shí)與下午的15：00—17：00時(shí)，在這些時(shí)段，園區(qū)供大于需，若將大量的富余電力送到電網(wǎng)，將帶來巨大的懲罰成本。因此，在場(chǎng)景1中，富余電力通過電解槽轉(zhuǎn)化為氫氣，實(shí)現(xiàn)新能源的消納，尤其是在新能源總出力最大的為09：00，電解槽的產(chǎn)氫量也達(dá)到了最大值，同時(shí)氫制氨也在最大功率運(yùn)行，盡可能地消納風(fēng)電與光伏。在14：00、19：00時(shí)，園區(qū)內(nèi)新能源發(fā)電難以滿足負(fù)荷需求。考慮到此時(shí)較高的購電電價(jià)，從外電網(wǎng)購電將帶來巨大的購電成本，不利于園區(qū)的收益最大化，因此，園區(qū)在這兩個(gè)時(shí)刻，利用氫燃料發(fā)電機(jī)進(jìn)行發(fā)電滿足負(fù)荷需要。同時(shí)，為緩解這些時(shí)刻電負(fù)荷及氫負(fù)荷壓力，氫制氨也處于最小功率運(yùn)行，從而大幅減少園區(qū)的外購電成本。

從圖6可以看出：在場(chǎng)景2中，電解槽處于定功率運(yùn)行狀態(tài)時(shí)。當(dāng)新能源出力最大時(shí)，電解槽運(yùn)行功率過低而無法消納過剩的新能源電力；當(dāng)園區(qū)負(fù)荷最大的時(shí)刻，電解槽也無法降低自身功率，反而增大了園區(qū)的購電負(fù)擔(dān)。在這種運(yùn)行方式下，園區(qū)與外界的電力交易更加頻繁且交易電量明顯增加，不利于園區(qū)整體的發(fā)展。

圖6 場(chǎng)景2運(yùn)行結(jié)果Fig. 6 Operation results of case 2

從圖7可以看出：在場(chǎng)景3中，氫制氨處于定功率運(yùn)行狀態(tài)，無法為園區(qū)提供可調(diào)的靈活型資源。例如：在可再生電源出力出現(xiàn)低谷的04：00、15：00及19：00時(shí)，適當(dāng)降低氫制氨的生產(chǎn)功率，可以減小園區(qū)總負(fù)荷，降低購電成本；在可再生電源出力的高峰09：00時(shí)，適當(dāng)增大氫制氨的產(chǎn)能，有利于促進(jìn)新能源的消納，最大限度地減少棄風(fēng)棄光。

圖7 場(chǎng)景3運(yùn)行結(jié)果Fig. 7 Operation results of case 3

3種場(chǎng)景下園區(qū)與電網(wǎng)的交易電量對(duì)比如圖8所示。表4展示了園區(qū)級(jí)電–氫耦合系統(tǒng)在不同運(yùn)行方式下的運(yùn)行成本及總收益。

圖8 3種場(chǎng)景的交易電量對(duì)比Fig. 8 Comparison of transacted electricity quantity for three cases

表4 3種場(chǎng)景的運(yùn)行成本對(duì)比Tab. 4 Comparison of operation cost for three cases

結(jié)合圖8和表4可以看出，在場(chǎng)景2中，由于電力交易的頻繁，導(dǎo)致園區(qū)購電成本大幅增加，總收益也大幅下降。場(chǎng)景3與場(chǎng)景1的對(duì)比表明，在考慮氫制氨的需求側(cè)響應(yīng)時(shí)，制氨收益基本不變，同時(shí)由于加入了可調(diào)的靈活性資源，可以使購電成本與棄風(fēng)棄光懲罰成本進(jìn)一步降低，進(jìn)而使園區(qū)總收益提升約21.5%。

5 結(jié) 論

由于新能源電力的裝機(jī)容量日益提高，導(dǎo)致大量富余電力從配電網(wǎng)注入輸電網(wǎng)絡(luò)或出現(xiàn)嚴(yán)重的棄風(fēng)棄光現(xiàn)象。本文以化工園區(qū)為例，通過電力系統(tǒng)與氫系統(tǒng)的耦合，來解決過剩電力難以消納的問題，在氫源側(cè)，通過電解水制氫來消納富余電力；在氫負(fù)荷側(cè)，通過氫燃料發(fā)電機(jī)在負(fù)荷高峰時(shí)刻進(jìn)行發(fā)電，減少在電價(jià)峰時(shí)的外購電量。本文從系統(tǒng)的源–網(wǎng)–荷–儲(chǔ)角度，建立了各側(cè)設(shè)備的物理模型，在此基礎(chǔ)上，建立了以化工園區(qū)為例的園區(qū)級(jí)電–氫耦合系統(tǒng)經(jīng)濟(jì)調(diào)度模型。從案例研究中可以看出，電解槽、氫燃料發(fā)電機(jī)、氫制氨通過跟蹤新能源出力與園區(qū)負(fù)荷的變化，實(shí)現(xiàn)對(duì)新能源電力的削峰填谷；同時(shí)可以得出，電–氫耦合系統(tǒng)協(xié)同運(yùn)行的基礎(chǔ)是電解槽跟隨源荷功率變化，可大幅減少園區(qū)購電與棄風(fēng)棄光，在此基礎(chǔ)上，考慮氫負(fù)荷側(cè)的需求側(cè)響應(yīng)，充分挖掘園區(qū)的靈活性資源，可以進(jìn)一步降低系統(tǒng)運(yùn)行成本，提升園區(qū)總體運(yùn)行效益。

由于考慮負(fù)荷側(cè)的需求側(cè)響應(yīng)可為園區(qū)帶來直觀的經(jīng)濟(jì)效益，因此，下一步將對(duì)氫負(fù)荷側(cè)進(jìn)行細(xì)致研究，建立更精細(xì)的模型來進(jìn)行優(yōu)化調(diào)度。此外，氫燃料電池的熱電聯(lián)供技術(shù)的應(yīng)用也逐漸廣泛，將出現(xiàn)越來越多的電–氫–熱耦合系統(tǒng)，今后可在本研究的基礎(chǔ)上，融合供熱系統(tǒng)，進(jìn)一步提升能效，降低成本。

（編輯 趙 婧）