基于能源集線器矩陣化建模的冷熱電多能流優(yōu)化調(diào)度

陸懷谷，莊重，陳亮，張偉，尹俊杰，王建華

(1. 國(guó)網(wǎng)江蘇省電力有限公司 常州供電分公司，江蘇 常州 213164；2. 東南大學(xué) 電氣工程學(xué)院，江蘇 南京 210018)

0 引言

在過(guò)去較長(zhǎng)一段時(shí)間內(nèi)，由于各類(lèi)能源系統(tǒng)彼此之間存在較少的耦合環(huán)節(jié)，因此實(shí)行不同能源系統(tǒng)相對(duì)獨(dú)立管理的模式，如電網(wǎng)、熱網(wǎng)、天然氣網(wǎng)等屬于不同公司獨(dú)立規(guī)劃、獨(dú)立管理和獨(dú)立維護(hù)[1]。面臨新形勢(shì)，過(guò)去運(yùn)行的管理模式會(huì)導(dǎo)致能源使用效率總體不高，而且面臨極端天氣等自然災(zāi)害(如臺(tái)風(fēng)、冰雹、地震)時(shí)，各類(lèi)供能系統(tǒng)(如供電、供氣、供熱、供冷等)之間缺乏必要的協(xié)調(diào)配合[2]，從而導(dǎo)致單方面的能源故障存在引發(fā)復(fù)雜連鎖反應(yīng)的風(fēng)險(xiǎn)。

為了將具有不同物理本質(zhì)的多種能源系統(tǒng)統(tǒng)一規(guī)劃調(diào)度，當(dāng)前學(xué)者已進(jìn)行較為深入的研究和探討。文獻(xiàn)類(lèi)比電路模型，并借鑒例如時(shí)域頻域變換、戴維南等值定理、集中參數(shù)等效等傳統(tǒng)電路分析方法[3]，將多能流的傳輸網(wǎng)絡(luò)進(jìn)行統(tǒng)一規(guī)劃，取得了較好的數(shù)學(xué)形式上的統(tǒng)一，亦有專(zhuān)家學(xué)者指出可能由于追求統(tǒng)一模型而引入了一些不必要的人為假設(shè)，值得進(jìn)一步深入探討研究以實(shí)現(xiàn)較為統(tǒng)一的模型形式和較為靈活的調(diào)度效果。部分研究著眼于提出“統(tǒng)一能路”的概念及其方法論，并以天然氣網(wǎng)絡(luò)、熱力網(wǎng)絡(luò)為研究對(duì)象推導(dǎo)了統(tǒng)一能路理論中的氣路、水路與熱路模型[4]，對(duì)于未來(lái)大規(guī)模的多能流耦合規(guī)劃提供了一定的指導(dǎo)意見(jiàn)。

進(jìn)一步地，國(guó)內(nèi)外在綜合能源系統(tǒng)矩陣化建模研究方面也進(jìn)行了較為豐富的討論研究。蘇黎世聯(lián)邦理工學(xué)院FAVRE-PERROD P教授團(tuán)隊(duì)提出能源集線器(energy hub, EH)這一概念，是用以表征不同能源載體的輸入、輸出、轉(zhuǎn)換、存儲(chǔ)的多端口裝置[5]。意大利學(xué)者率先提出在冷熱電三聯(lián)供(combined cooling, heating and power，CCHP)系統(tǒng)中根據(jù)輸出能流情況，基于“回溯”算法探求能流輸入情況和EH內(nèi)部的分布情況[6]，較為細(xì)致地刻畫(huà)了三聯(lián)供系統(tǒng)內(nèi)部的能流分布及其轉(zhuǎn)化關(guān)系。能源集線器標(biāo)準(zhǔn)化矩陣建模在定義了節(jié)點(diǎn)-端口能量轉(zhuǎn)化和端口-支路關(guān)聯(lián)關(guān)系等基礎(chǔ)上，能夠?qū)崿F(xiàn)計(jì)算機(jī)自動(dòng)化建模，并實(shí)現(xiàn)綜合能源系統(tǒng)線性化、標(biāo)準(zhǔn)化建模，適用于多種結(jié)構(gòu)的綜合能源系統(tǒng)[7-8]。然而也存在能流之間部分耦合關(guān)系被忽略簡(jiǎn)化，導(dǎo)致矩陣中出現(xiàn)大量“0”元素和為了追求矩陣描述形式而把簡(jiǎn)單問(wèn)題復(fù)雜化的問(wèn)題，尤其是忽略了能源集線器內(nèi)部元件運(yùn)行工況的實(shí)時(shí)變化。在規(guī)劃調(diào)度方面，文獻(xiàn)[9-10]在已知各類(lèi)能源年度典型日負(fù)荷曲線、輸入能源最大容量與單位價(jià)格、可選設(shè)備及其投資價(jià)格等數(shù)據(jù)情況下，對(duì)區(qū)域綜合能源系統(tǒng)的設(shè)備選擇、連接方式以及典型日的設(shè)備運(yùn)行方式進(jìn)行優(yōu)化。

在綜合能源系統(tǒng)中，多能流耦合關(guān)系較為緊密、耦合環(huán)節(jié)較為復(fù)雜、耦合形式較為典型的是冷熱電三聯(lián)供系統(tǒng)[11]。該系統(tǒng)的設(shè)計(jì)運(yùn)行模式建立在能量梯級(jí)利用的概念之上，首先通過(guò)燃?xì)廨啓C(jī)燃燒天然氣產(chǎn)生高品位電能，再將低品位熱能通過(guò)余熱鍋爐、制冷機(jī)等設(shè)備轉(zhuǎn)化用于供熱和制冷。根據(jù)數(shù)據(jù)統(tǒng)計(jì)，傳統(tǒng)的集中發(fā)電、遠(yuǎn)程送電模式的發(fā)電效率較低，且在輸配電過(guò)程中存在較大的電能損耗。類(lèi)似地，集中輸送熱能每公里有一定的溫度降落，因此一般規(guī)劃過(guò)程中傳輸距離設(shè)計(jì)僅為數(shù)公里。然而，經(jīng)過(guò)能量梯級(jí)利用的冷熱電三聯(lián)供系統(tǒng)使能源利用效率從傳統(tǒng)發(fā)電模式的40%左右提高到80%～90%，且輸電損耗相對(duì)較少[12]。其優(yōu)勢(shì)在于可以在系統(tǒng)內(nèi)引進(jìn)轉(zhuǎn)化效率更高的機(jī)組，實(shí)現(xiàn)電、氣、熱、冷的最優(yōu)分配轉(zhuǎn)化模式，顯著提高多種能源利用效率；可以強(qiáng)化居民用戶之間、建筑樓宇之間、生產(chǎn)企業(yè)之間的有效連接，較大程度上實(shí)現(xiàn)能源共享；可以有效融入太陽(yáng)能發(fā)電、光伏集熱、風(fēng)能發(fā)電、地?zé)崂玫瓤稍偕茉矗瑥亩行p少二氧化碳排放量，推動(dòng)能源結(jié)構(gòu)的升級(jí)轉(zhuǎn)型[13]。

針對(duì)上述問(wèn)題，本文構(gòu)建了能源集線器的通用能量交換分析模型，通過(guò)矩陣一體化描述能源集線器的能量輸入端、能量輸出端以及集線器內(nèi)部的拓?fù)潢P(guān)系，集中明確地表征多能流之間的供能方式、用能形式和耦合關(guān)系。隨后對(duì)冷熱電三聯(lián)供系統(tǒng)進(jìn)行詳細(xì)闡述，建立燃?xì)廨啓C(jī)、電制冷機(jī)、輔助鍋爐等內(nèi)部元件矩陣化模型，完整描述多能流的輸入、輸出耦合關(guān)系并進(jìn)行一定程度的解耦，最后通過(guò)算例分析驗(yàn)證了所述矩陣化模型在多能流耦合優(yōu)化調(diào)度方面的作用。

1 通用能量交換分析模型

當(dāng)前綜合能源系統(tǒng)方面的研究多專(zhuān)注于考慮更為復(fù)雜的情況和創(chuàng)新優(yōu)化調(diào)度理論，而瑞士學(xué)者較早提出能源集線器的概念：它被認(rèn)為是可以轉(zhuǎn)換、調(diào)節(jié)和存儲(chǔ)多個(gè)能量載體的單元，代表了不同能源基礎(chǔ)設(shè)施和負(fù)載之間的接口。圖1表述了能源集線器概念。

圖1 通用能量交換分析模型

在研究過(guò)程中，穩(wěn)態(tài)流量模型相對(duì)暫態(tài)流量模型來(lái)說(shuō)更為簡(jiǎn)潔，且針對(duì)較小范圍的能源集線器來(lái)說(shuō)也更為適用。對(duì)流經(jīng)能量轉(zhuǎn)換設(shè)備的流量進(jìn)行簡(jiǎn)單分析，將其能量效率定義為穩(wěn)態(tài)輸出與輸入之比。通過(guò)多個(gè)輸入和輸出，可以定義一個(gè)元素與功率流一一對(duì)應(yīng)的轉(zhuǎn)換矩陣。耦合矩陣描述了能流從集線器的輸入到輸出的轉(zhuǎn)換，元素可從集線器的轉(zhuǎn)換器結(jié)構(gòu)和轉(zhuǎn)換器的效率特性中得出。從通用網(wǎng)絡(luò)流量到更詳細(xì)的穩(wěn)態(tài)功率流量模型，電、氣、熱等各種流量均適用，其適用的程度取決于所進(jìn)行的研究類(lèi)型。

能源樞紐左端的P向量表示多能源系統(tǒng)的原始能源輸入，下標(biāo)m∈M表示輸入能源種類(lèi)；右端的L向量表示經(jīng)過(guò)轉(zhuǎn)換后的能源輸出，下標(biāo)n表示輸出能源種類(lèi)。因此數(shù)學(xué)層面上的能源集線器是一個(gè)映射函數(shù)：

L=f(P)

(1)

式中函數(shù)f(·)可以考慮到各種形式能源的傳輸、轉(zhuǎn)換、存儲(chǔ)等環(huán)節(jié)。因此可以利用一個(gè)耦合矩陣來(lái)描述輸入-輸出之間的關(guān)系：

(2)

式中cij為耦合因子，表示第j種形式能源輸出與第i種形式能源輸入的比值。各種形式能源的輸入到輸出可以分為兩個(gè)步驟：能源分配和能源傳輸或轉(zhuǎn)化。能源分配是指各種能源以一定的比例分配到不同能源傳輸或轉(zhuǎn)化設(shè)備。能源傳輸或轉(zhuǎn)化是指能源輸入該設(shè)備后通過(guò)機(jī)械、化學(xué)等途徑進(jìn)行轉(zhuǎn)化，具有一定的轉(zhuǎn)化效率。所以可以將式(2)中的耦合矩陣作進(jìn)一步的分解：

(3)

可以簡(jiǎn)記為

L=ηNP=CP

(4)

式中：vij為分配因子，表示第i種形式能源輸入分配到用于轉(zhuǎn)化第j種形式能源的比例，對(duì)應(yīng)的N為分配矩陣；ηij為效率因子，表示第i種形式能源轉(zhuǎn)化為第j種形式能源的效率，對(duì)應(yīng)的η為轉(zhuǎn)化矩陣。

2 CCHP系統(tǒng)矩陣化建模

2.1 CCHP系統(tǒng)及其內(nèi)部元件

圖2展示了基于CCHP的能源集線器，該系統(tǒng)用于詳細(xì)說(shuō)明如何建立能源集線器模型，它由熱電聯(lián)供 (combined heat and power，CHP) 單元、輔助鍋爐(auxiliary boiler，AB)、電制冷機(jī)(electric refrigerator，ER)和吸收式制冷機(jī)(absorption refrigerator，AR)組成，其中熱電聯(lián)供單元包括燃?xì)廨啓C(jī)(gas turbine，GT)和余熱鍋爐(waste-heat boiler，WB)兩部分。

圖2 CCHP系統(tǒng)及其內(nèi)部元件示意圖

能源輸入端為電能和天然氣兩種，分別記為Pe、Pg；能源輸出端連接電負(fù)荷Le、熱負(fù)荷Lh、冷負(fù)荷Lc三種負(fù)荷形式。

輸入電能Pe從電力傳輸系統(tǒng)獲得，進(jìn)入能源集線器，一部分直接用于電負(fù)荷，其電能損耗忽略不計(jì)，另一部分進(jìn)入電制冷機(jī)ER轉(zhuǎn)化為冷能，滿足部分冷負(fù)荷需求。

(5)

輸入電能Pe經(jīng)分配矩陣Ne作用后，分為Pe1和Pe2兩部分，其中ve1、ve2分別表示分配比例，且ve1+ve2=1。

天然氣Pg通過(guò)氣網(wǎng)輸送至冷熱電三聯(lián)供系統(tǒng)中，Pg1、Pg2部分用于熱電聯(lián)供單元CHP(包括燃?xì)廨啓C(jī)GT和余熱鍋爐WB)，另一部分Pg3用于輔助鍋爐AB。

(6)

輸入天然氣Pg經(jīng)分配矩陣Ng作用后，分為Pg1、Pg2和Pg3，其中vg1、vg2、vg3表示分配比例，且vg1+vg2+vg3=1。需要注意的是，在實(shí)際生產(chǎn)過(guò)程中，天然氣Pg1、Pg2首先全部注入燃?xì)廨啓C(jī)GT燃燒轉(zhuǎn)化，隨后利用余熱鍋爐WB進(jìn)行熱能回收，并非圖2中所示天然氣Pg1、Pg2分別進(jìn)入燃?xì)廨啓C(jī)GT、余熱鍋爐AB。如果按照實(shí)際生產(chǎn)次序進(jìn)行建模過(guò)于復(fù)雜，為簡(jiǎn)化計(jì)算采用上述表達(dá)方式；且從能量轉(zhuǎn)化形式這一本質(zhì)上來(lái)說(shuō)，將Pg1視為注入燃?xì)廨啓C(jī)GT產(chǎn)生電能，Pg2視為注入余熱鍋爐AB產(chǎn)生熱能也有一定道理。

對(duì)于燃?xì)廨啓C(jī)GT來(lái)說(shuō)，燃燒生產(chǎn)的電能一部分直接用于電負(fù)荷，另一部分經(jīng)過(guò)電制冷機(jī)ER生成冷能供應(yīng)給冷負(fù)荷。由圖2可知，燃?xì)廨啓C(jī)為輸入一端口、輸出二端口形式：

(7)

對(duì)于余熱鍋爐WB來(lái)說(shuō)，收集的熱能部分用于熱負(fù)荷，其余部分用于吸收式制冷機(jī)轉(zhuǎn)化為冷能。

(8)

對(duì)于輔助鍋爐AB來(lái)說(shuō)，

(9)

對(duì)于電制冷機(jī)ER來(lái)說(shuō)，

(10)

對(duì)于吸收式制冷機(jī)AR來(lái)說(shuō)，

(11)

2.2 輸入、輸出矩陣化模型

由圖2所示的基于CCHP系統(tǒng)能源集線器，根據(jù)式(2)-式(4)可得輸入矩陣如下：

(12)

經(jīng)分配矩陣作用后的分配矩陣P′如下：

(13)

從輸出端向能源集線器內(nèi)部看，輸出矩陣即負(fù)載矩陣L：

(14)

聯(lián)立式(12)-式(14)，并根據(jù)式(4)的化簡(jiǎn)形式L=CP即可求得輸入、輸出的耦合矩陣C。

2.3 優(yōu)化目標(biāo)與約束條件

優(yōu)化目標(biāo)：對(duì)能源服務(wù)商而言，能源集線器向用能負(fù)荷輸出的電能、熱能、冷能，分別賦予“售能單價(jià)”；輸入的電能、天然氣從電力系統(tǒng)和天然氣系統(tǒng)購(gòu)入，分別賦予“購(gòu)能單價(jià)”。在每個(gè)時(shí)隙內(nèi)，滿足用戶用能需求的前提下，實(shí)現(xiàn)盈利最大化。

約束條件：1)每一時(shí)隙分配矩陣中的元素都在各自合理區(qū)間內(nèi)動(dòng)態(tài)變化，且滿足不同種類(lèi)輸入能源集線器的能源分配比例之和為1；2)燃?xì)廨啓C(jī)GT、電制冷機(jī)ER等單元設(shè)備的輸入能源量應(yīng)處于(min, max)合理的范圍之內(nèi)；3)轉(zhuǎn)化單元的轉(zhuǎn)化效率與設(shè)備性能、環(huán)境因素等有緊密關(guān)系，同樣應(yīng)處于合理區(qū)間之內(nèi)；4)還應(yīng)滿足能源集線器內(nèi)部能流轉(zhuǎn)化的等式約束，如式(14)所示。

公式化表達(dá)如下：

s.t. ∑vm=1, ?m∈M,

Pmin≤Pi≤Pmax, ?i,

ηmin≤ηi≤ηmax, ?i,

&式(14)

(15)

功率乘時(shí)間等于能量，即P×t=E。本文所考慮的時(shí)間間隙均為1 h，因此公式(15)表示利潤(rùn)π等于盈利R減去成本C，省略時(shí)間常數(shù)。其中，α為服務(wù)商向用戶出售能源的單價(jià)；β為服務(wù)商從系統(tǒng)購(gòu)入能源的單價(jià)。

已知服務(wù)商從系統(tǒng)購(gòu)入能源單價(jià)β根據(jù)提前簽訂的合同，向用戶出售能源單價(jià)α在基準(zhǔn)價(jià)格的基礎(chǔ)上產(chǎn)生一定波動(dòng)。本算例中的變量包括：①可控變量，分配比例v和購(gòu)入能源量P，兩者為優(yōu)化變量；②不可控變量，轉(zhuǎn)化效率η，僅可通過(guò)分配比例v控制CCHP系統(tǒng)內(nèi)各單元的啟停狀態(tài)。

3 算例分析

3.1 原始數(shù)據(jù)

為了驗(yàn)證本文方法的有效性，使用MATLAB -YALMIP平臺(tái)開(kāi)發(fā)了上述優(yōu)化程序，該算例分析中所用的硬件環(huán)境是Intel(R) Core(i) i7-6500 M CPU@2.50 GHz，8 GB RAM，Win10 64位(操作系統(tǒng))，MATLAB R2019 b(開(kāi)發(fā)環(huán)境)，YALMIP版本是R20200116。

fmincon函數(shù)可以被用來(lái)求解帶約束條件的非線性多變量函數(shù)(constrained nonlinear multivariable function)的最小值,即可以用來(lái)求解帶約束的非線性規(guī)劃問(wèn)題。其非線性規(guī)劃模型表達(dá)式中b、beq是向量，A、Aeq是矩陣，c(x)和ceq(x)是返回向量的函數(shù)，而f(x)是返回標(biāo)量的函數(shù)，f(x)、c(x)和ceq(x)可以是非線性函數(shù)。

(16)

能源服務(wù)商購(gòu)入和售出能源單價(jià)如表1所示。在實(shí)際情況中，價(jià)格隨供需等多種因素而波動(dòng)。因此，在算例分析中單價(jià)隨機(jī)波動(dòng)±5%，每小時(shí)更新一次。在未來(lái)研究中，可將負(fù)荷響應(yīng)、價(jià)格杠桿納入考量范圍之內(nèi)。

表1 能源服務(wù)商購(gòu)入和售出能源單價(jià)

根據(jù)表2可以得知，本算例所考慮的CCHP系統(tǒng)中轉(zhuǎn)化單元輸入功率和轉(zhuǎn)化效率范圍。根據(jù)統(tǒng)計(jì)計(jì)算結(jié)果給定t=0時(shí)刻的初始數(shù)據(jù)后，為后續(xù)每一時(shí)刻的轉(zhuǎn)化效率增加random隨機(jī)函數(shù)，用于模擬各個(gè)元件運(yùn)行轉(zhuǎn)化效率的不確定性。

表2 轉(zhuǎn)化單元輸入功率和轉(zhuǎn)化效率范圍

夏季、冬季典型日負(fù)荷需求24 h內(nèi)隨時(shí)間變化如圖3所示。根據(jù)“需求決定供給”的原則，綜合能源服務(wù)商從系統(tǒng)中購(gòu)入能源量由該時(shí)隙的能源需求量及其種類(lèi)所決定。由圖3可見(jiàn)，電負(fù)荷需求基本保持穩(wěn)定，符合峰谷波動(dòng)規(guī)律，且夏季用電需求較大；熱負(fù)荷需求在夏季7時(shí)、20時(shí)附近出現(xiàn)峰值，在冬季中午前后出現(xiàn)谷值但基本平穩(wěn)；冷負(fù)荷在夏季15時(shí)前后出現(xiàn)峰值，較大需求范圍保持6 h左右，日出前后需求最少，而冬季對(duì)冷負(fù)荷需求整體維持在較低水平。

圖3 夏季、冬季典型日負(fù)荷需求曲線圖

3.2 仿真結(jié)果

夏季、冬季典型日利潤(rùn)曲線隨時(shí)間波動(dòng)情況如圖4所示。對(duì)比圖3的夏季典型日負(fù)荷需求可知，利潤(rùn)曲線走勢(shì)與冷負(fù)荷需求走勢(shì)基本一致，均為10時(shí)-20時(shí)出現(xiàn)較大值，且18時(shí)附近利潤(rùn)最大，此時(shí)電、熱、冷負(fù)荷需求量之和也接近峰值。類(lèi)似地，將圖3的冬季典型日負(fù)荷需求曲線與圖4冬季利潤(rùn)曲線對(duì)比可知，利潤(rùn)基本保持平穩(wěn)，與冬季典型日負(fù)荷需求走勢(shì)類(lèi)似。從圖4夏季、冬季利潤(rùn)曲線對(duì)比可知，夏季能源收益在白天時(shí)段高于冬季能源收益，而夜間基本保持一致，甚至低于冬季能源收益。

圖4 夏季、冬季典型日利潤(rùn)曲線圖

為了直觀地體現(xiàn)本文所提調(diào)度方案在能源服務(wù)商收益方面的優(yōu)勢(shì)，在算例分析中采用了3種不同的方案進(jìn)行對(duì)比。方案1)本文所提出的CCHP系統(tǒng)能源集線器調(diào)度方法，根據(jù)情況實(shí)時(shí)調(diào)整購(gòu)入電、氣能源量以及各轉(zhuǎn)化單元的分配量；方案2)電能由電力系統(tǒng)直接滿足負(fù)荷需求，熱能和冷能由天然氣系統(tǒng)燃燒轉(zhuǎn)化以滿足熱負(fù)荷、冷負(fù)荷需求；方案3)天然氣直接通過(guò)燃?xì)廨啓C(jī)等設(shè)備滿足電、熱、冷負(fù)荷需求，不再?gòu)碾娏ο到y(tǒng)購(gòu)入。

以夏季典型日為例，分別采用以上3種方案進(jìn)行優(yōu)化調(diào)度，能源服務(wù)商利潤(rùn)如圖5所示。整體上來(lái)說(shuō)，3種方案的利潤(rùn)走勢(shì)基本相近。將13時(shí)-18時(shí)的曲線細(xì)節(jié)放大可知，方案1明顯優(yōu)于其他兩種方案，每個(gè)時(shí)段不同程度地獲得更多收益，約為6%左右。進(jìn)一步來(lái)說(shuō)，除了考慮能源服務(wù)商的收益之外，能源消耗總量也應(yīng)納入考量范圍之內(nèi)，這與溫室氣體的排放、不可再生能源儲(chǔ)量等密切相關(guān)。通過(guò)計(jì)算可知，方案1所消耗的能源總量也較少，合理采用不同種類(lèi)的能源以滿足用能需求，一定程度上實(shí)現(xiàn)了綠色節(jié)能生產(chǎn)供應(yīng)。進(jìn)一步地，將冬季、夏季典型日各分配單元在不同時(shí)隙的分配比例用累積柱形圖表示，從電力系統(tǒng)、天然氣系統(tǒng)購(gòu)入的電能、天然氣用梯形圖表示，如圖6、圖7所示。

圖5 夏季不同供能方案利潤(rùn)對(duì)比圖

根據(jù)圖6(a)柱形圖可知，從電力系統(tǒng)獲得的電能經(jīng)分配后主要用于電負(fù)荷需求，少部分用于電制冷機(jī)ER轉(zhuǎn)化為冷能；從天然氣系統(tǒng)獲得的天然氣絕大部分注入熱電聯(lián)供系統(tǒng)(約2/3用于燃?xì)廨啓C(jī)GT，約1/3用于余熱鍋爐WB)，極少部分用于輔助鍋爐AB，當(dāng)熱電聯(lián)供系統(tǒng)滿載運(yùn)行仍無(wú)法滿足用能需求時(shí)，才啟用輔助鍋爐AB進(jìn)行補(bǔ)充，符合系統(tǒng)設(shè)計(jì)初衷；燃?xì)廨啓C(jī)GT所發(fā)出的電能少部分注入電制冷機(jī)ER中進(jìn)行電-冷轉(zhuǎn)化，大部分直接用于電負(fù)荷供應(yīng)，這是由于天然氣成本較低，雖然經(jīng)燃?xì)廨啓C(jī)GT、電制冷機(jī)ER等設(shè)備轉(zhuǎn)化過(guò)程中有一定的能量損失，相對(duì)于從電力系統(tǒng)直接購(gòu)能來(lái)說(shuō)更為經(jīng)濟(jì)，一定程度上避免了遠(yuǎn)距離輸電的能量損耗；余熱鍋爐WB絕大部分時(shí)間直接供應(yīng)熱負(fù)荷，在冷負(fù)荷需求高峰時(shí)段會(huì)將注入吸收式制冷機(jī)AR的分配比例增大，以滿足冷負(fù)荷需求；由于輔助鍋爐AB的轉(zhuǎn)化效率較高，主要用于余熱鍋爐WB供應(yīng)冷負(fù)荷時(shí)，余熱鍋爐會(huì)投入生產(chǎn)。

類(lèi)似地，圖7(a)柱形圖一定程度上驗(yàn)證了上述分析。縱向?qū)Ρ认募尽⒍镜姆峙浔壤齺?lái)看，購(gòu)入電能、余熱鍋爐WB和輔助鍋爐AB在時(shí)間變化上峰谷基本相反。這是因?yàn)橄募纠湄?fù)荷需求較大，且時(shí)段性較強(qiáng)，在中午時(shí)刻附近各設(shè)備會(huì)充分分配能源注入電制冷機(jī)ER和吸收式制冷機(jī)AR中，以滿足陡然增加的冷負(fù)荷需求。

對(duì)比圖6(b)、圖7(b)的梯形圖可知，從系統(tǒng)購(gòu)入的電力和天然氣對(duì)系統(tǒng)的峰谷調(diào)節(jié)也有一定的積極作用。尤其是冬季典型日的天然氣購(gòu)入量基本保持平穩(wěn)，這對(duì)整個(gè)天然氣系統(tǒng)的平穩(wěn)運(yùn)行具有重要意義；夜間購(gòu)入較多天然氣以滿足供熱需求，也能保證整體系統(tǒng)的穩(wěn)定出氣量。夏季典型日的電能購(gòu)入量對(duì)電力系統(tǒng)削峰填谷作用更為明顯，白天時(shí)段正常是用電高峰期，通過(guò)CCHP系統(tǒng)的調(diào)節(jié)使中午前后的6 h高峰時(shí)間購(gòu)入電能減少，并用天然氣代替。由于本算例僅考慮單個(gè)用戶的用能需求，因此CCHP系統(tǒng)購(gòu)入能源波動(dòng)起伏較大，如果將此方案推廣至園區(qū)級(jí)甚至跨區(qū)級(jí)范圍，其削峰填谷、節(jié)能高效的特點(diǎn)應(yīng)當(dāng)更為明顯。

圖6 冬季典型日分配參數(shù)累積柱形圖、輸入能源梯形圖

圖7 夏季典型日分配參數(shù)累積柱形圖、輸入能源梯形圖

4 結(jié)語(yǔ)

本文提出的冷熱電三聯(lián)供系統(tǒng)能源集線器矩陣化建模方法，可以得到如下結(jié)論。

1)通過(guò)矩陣一體化描述能源集線器的能量輸入端、能量輸出端以及集線器內(nèi)部的拓?fù)潢P(guān)系，集中明確地表征多能流之間的供能方式、用能形式和耦合關(guān)系。

2)實(shí)現(xiàn)對(duì)能源梯級(jí)利用，有效提高綜合利用率。冷熱電三聯(lián)供系統(tǒng)的綜合能源利用率比傳統(tǒng)供能方式有顯著提升，而且由于靠近用能負(fù)荷可以減少能量傳輸損失。

3)對(duì)電力系統(tǒng)和天然氣系統(tǒng)有削峰谷的作用，合理優(yōu)化多能流之間的轉(zhuǎn)化關(guān)系，增強(qiáng)不同供能系統(tǒng)的互濟(jì)協(xié)調(diào)能力，提高用能安全性和可靠性。

4)對(duì)比其他供能方案，該方法既實(shí)現(xiàn)了節(jié)能環(huán)保、減少碳排放的綠色發(fā)展目標(biāo)，保證了用戶的多樣化用能需求，又使能源服務(wù)商利益最大化。

該類(lèi)問(wèn)題研究的難點(diǎn)在于耦合矩陣的分離和解耦，研究更為快速準(zhǔn)確的解耦方式對(duì)于能源集線器優(yōu)化速度的提升具有顯著效果, 有待在未來(lái)工作中深入探討。