李恩騰，徐英杰，謝小東，范偉

（浙江工業(yè)大學(xué)機(jī)械工程學(xué)院，浙江杭州310023）

隨著人們生活品質(zhì)的提高，家用熱水的需求量也在不斷增加，然而在國(guó)內(nèi)熱水超過90%都是直接燃燒化石燃料加熱，造成了嚴(yán)重的環(huán)境污染與破壞[1-2]。因此，高效清潔的跨臨界CO2熱泵熱水器具有較大開發(fā)價(jià)值。CO2具有的天然環(huán)保特性以及跨臨界CO2熱泵系統(tǒng)對(duì)于加熱高溫?zé)崴憩F(xiàn)出來的巨大優(yōu)勢(shì)，使得跨臨界CO2熱泵熱水器在住宅中的應(yīng)用具有廣闊的前景[3]。然而跨臨界CO2熱泵熱水器也存在一定的問題；首先，跨臨界CO2熱泵熱水器的投資成本太高，遠(yuǎn)遠(yuǎn)高于電、燃?xì)?、太?yáng)能熱水器的成本，很難被普通用戶所接受[4]；其次，跨臨界CO2熱泵熱水器的儲(chǔ)熱水箱會(huì)占用較大的室內(nèi)空間，給用戶的使用帶來很大不便。因此，CO2熱泵熱水器需要綜合考慮經(jīng)濟(jì)性和實(shí)用性，對(duì)設(shè)備的性能和設(shè)計(jì)進(jìn)行優(yōu)化改進(jìn)。

Lorentzen 與Pettersen[5-6]率先對(duì)CO2跨臨界循環(huán)進(jìn)行了研究，研究表明CO2熱泵系統(tǒng)在性能上可以達(dá)到較高的性能系數(shù)（COP），能滿足高溫?zé)崴a(chǎn)的需要，并且比電熱水器和燃?xì)鉄崴鳒p少約75%的能源消耗；此外對(duì)于未來環(huán)境問題而言，也避免許多不確定性因素的產(chǎn)生。Saikawa 和Koyama[7]從理論上研究了各種制冷劑單級(jí)壓縮熱泵循環(huán)的COP 上限，發(fā)現(xiàn)制冷劑CO2具有最高的COP。Qi 等[8]進(jìn)行了一項(xiàng)評(píng)估跨臨界CO2空氣源熱泵熱水器最佳排氣壓力的實(shí)驗(yàn)研究，以實(shí)現(xiàn)系統(tǒng)的最佳性能。Nawaz等[9]研究了CO2熱泵熱水器水循環(huán)速率對(duì)水箱中水溫分層的影響，表明這是可以提高CO2熱泵系統(tǒng)性能的關(guān)鍵因素。Fernandez 等[10]發(fā)現(xiàn)跨臨界CO2熱泵熱水器的熱水在水箱中長(zhǎng)期儲(chǔ)存，由于熱傳導(dǎo)的影響會(huì)使中溫水區(qū)域擴(kuò)大導(dǎo)致系統(tǒng)性能下降。為了解決該問題，Ohkura等[11]采用了從跨臨界CO2熱泵水箱側(cè)面提取溫水的方法來克服系統(tǒng)性能下降的缺點(diǎn)，同時(shí)使用數(shù)值模擬分析了這種跨臨界CO2熱泵改進(jìn)水箱的加熱系統(tǒng)性能，并通過實(shí)驗(yàn)驗(yàn)證其有效性。以上研究的重點(diǎn)主要集中在優(yōu)化跨臨界CO2熱泵系統(tǒng)的性能系數(shù)，然而單從性能這一參數(shù)上進(jìn)行優(yōu)化，可能會(huì)造成?損和經(jīng)濟(jì)成本的增加，所以對(duì)跨臨界CO2熱泵進(jìn)行多目標(biāo)優(yōu)化具有重要的現(xiàn)實(shí)意義。

Rezayan 和Behbahaninia[12]對(duì)CO2/NH3級(jí)聯(lián)熱泵循環(huán)進(jìn)行了?分析和經(jīng)濟(jì)性優(yōu)化，在制冷量、環(huán)境溫度和制冷空間溫度約束下，以最大限度地降低該級(jí)聯(lián)系統(tǒng)的總年度成本。結(jié)果表明，與基礎(chǔ)的設(shè)計(jì)系統(tǒng)相比，該優(yōu)化系統(tǒng)的總成本降低了9.34%。雖然該系統(tǒng)對(duì)?進(jìn)行了分析，但僅對(duì)單一的經(jīng)濟(jì)成本進(jìn)行優(yōu)化，無法滿足用戶的多元化需求。Fazelpour和Morosuk[13]從跨臨界CO2熱泵的性能和成本角度進(jìn)行了優(yōu)化分析，采用經(jīng)濟(jì)器作為單級(jí)跨臨界CO2熱泵的輔助裝置降低了14% 的最終成本。Aminyavari 等[14]采用遺傳算法對(duì)CO2/NH3級(jí)聯(lián)制冷系統(tǒng)的?效率、經(jīng)濟(jì)性和環(huán)境影響進(jìn)行了系統(tǒng)優(yōu)化。優(yōu)化結(jié)果表明，對(duì)于制冷量為50kW的工廠案例，Pareto 前沿選擇的最佳設(shè)計(jì)?效率為45.89%，年度總成本為277070USD/a；但是該優(yōu)化過程的熱力學(xué)模型復(fù)雜，優(yōu)化的計(jì)算量大。因此，Khanmohammadi等[15]采用了人工神經(jīng)網(wǎng)絡(luò)的方法來簡(jiǎn)化熱力學(xué)模型，以制冷量和年度成本為目標(biāo)對(duì)跨臨界CO2制冷循環(huán)進(jìn)行熱經(jīng)濟(jì)優(yōu)化，最終增加了4.6%的制冷量，降低了60.6%的年度成本。

如上所述，跨臨界CO2熱泵無論是對(duì)性能目標(biāo)進(jìn)行改進(jìn)，還是對(duì)性能系數(shù)、?效率、經(jīng)濟(jì)成本等多目標(biāo)進(jìn)行優(yōu)化都未涉及實(shí)用性的分析，然而，對(duì)于實(shí)用性的優(yōu)化有重要的現(xiàn)實(shí)意義。因此，本文基于BP 神經(jīng)網(wǎng)絡(luò)與NSGA-Ⅱ算法以跨臨界CO2熱泵的總經(jīng)濟(jì)成本與水箱體積為目標(biāo)函數(shù)建立了跨臨界CO2熱泵的優(yōu)化設(shè)計(jì)模型，并通過對(duì)常用CO2熱泵住宅案例進(jìn)行優(yōu)化為用戶提供最優(yōu)設(shè)計(jì)方案。

1 系統(tǒng)介紹及模型建立

1.1 跨臨界CO2熱泵模型建立

跨臨界CO2熱泵系統(tǒng)主要由壓縮機(jī)、氣冷器、蒸發(fā)器、膨脹閥、閥門、水泵、水箱等部件組成，結(jié)構(gòu)如圖1 所示。其循環(huán)工作原理是，CO2在蒸發(fā)器中吸收空氣熱量被蒸發(fā)為低溫低壓的氣體，然后在壓縮機(jī)中被加壓成高溫高壓的超臨界狀態(tài)，通過氣冷器將熱量釋放出來用于熱水加熱，加熱的熱水儲(chǔ)存在保溫水箱中，供給于用戶使用，在氣冷器中被冷卻了的CO2經(jīng)由膨脹閥節(jié)流降壓再進(jìn)入蒸發(fā)器中吸收空氣熱量蒸發(fā)，從而構(gòu)成了跨臨界CO2熱泵循環(huán)。CO2跨臨界壓縮循環(huán)的溫熵圖如圖2 所示。其熱力循環(huán)數(shù)學(xué)模型見式(1)～式(6)。

圖1 跨臨界CO2熱泵循環(huán)原理

圖2 跨臨界CO2循環(huán)溫熵圖

其中，CO2壓縮機(jī)的等熵效率ηS是關(guān)于壓比的多項(xiàng)式函數(shù)[式(3)][16]。

其中，p1、p2為吸排氣壓力。制熱量、COP、制冷量分別由式(4)～(6)表示。

氣冷器換熱溫差由式(7)計(jì)算。

對(duì)于該CO2熱泵熱水加熱系統(tǒng)，基于電網(wǎng)的峰谷電價(jià)制，熱泵機(jī)組可在夜間（低電價(jià)時(shí)段）運(yùn)行生產(chǎn)熱水儲(chǔ)存在水箱中用于用戶的熱水使用，因此體積越大的水箱能儲(chǔ)存更多的熱水，意味著運(yùn)營(yíng)成本就越低。然而更大的水箱也會(huì)導(dǎo)致占用更大的室內(nèi)空間、更大的水箱投資成本以及水箱散熱面積增大所帶來的漏熱損失的增加，從而導(dǎo)致用戶安裝的不便以及成本上升；并且，熱泵自身的設(shè)計(jì)參數(shù)，如氣冷器換熱器溫差，也會(huì)對(duì)熱泵系統(tǒng)的投資成本、運(yùn)營(yíng)成本產(chǎn)生影響。因此，更小的經(jīng)濟(jì)成本與更小的水箱占用空間是存在沖突關(guān)系的兩個(gè)目標(biāo)。所以，對(duì)于CO2熱泵熱水加熱系統(tǒng)而言，其設(shè)計(jì)工作還存在很大的優(yōu)化空間。

1.2 神經(jīng)網(wǎng)絡(luò)模型的建立

由于跨臨界CO2熱泵熱物理模型的復(fù)雜性，本文將采用神經(jīng)網(wǎng)絡(luò)的方法建立CO2熱泵的熱力學(xué)性能預(yù)測(cè)模型。人工神經(jīng)網(wǎng)絡(luò)技術(shù)可以被認(rèn)為是一種回歸技術(shù)，它可以表示自變量和因變量之間更高的非線性關(guān)系[17]。本研究采用的神經(jīng)網(wǎng)絡(luò)技術(shù)是前饋反向傳播（BP）神經(jīng)網(wǎng)絡(luò)。BP 神經(jīng)網(wǎng)絡(luò)模型通常由輸入層、隱藏層和輸出層組成，如圖3所示。隱藏層可以存在多個(gè)，然而，已經(jīng)證明在有足夠數(shù)量的神經(jīng)元情況下，一個(gè)隱藏層足以逼近任何連續(xù)的函數(shù)[18]，并且輸入或隱藏層可具有偏置神經(jīng)元。在不同層的神經(jīng)元之間具有連通性，每一層中的神經(jīng)元與其輸入層和輸出層的所有神經(jīng)元節(jié)點(diǎn)相連。本研究的輸入層為氣冷器換熱溫差（Δtm），輸出層為氣冷器壓力（pgc）、性能系數(shù)（COP）、壓縮機(jī)功率（Wcom）以及制冷量（Qe）。隱藏層的節(jié)點(diǎn)值由以下傳遞函數(shù)[式(8)]計(jì)算。

圖3 BP神經(jīng)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)

其中，Ni為隱藏層第i個(gè)神經(jīng)元節(jié)點(diǎn)；X為輸入層神經(jīng)元節(jié)點(diǎn)，即氣冷器換熱溫差（Δtm）；ω1i和b11i為輸入層與隱藏層之間的權(quán)重與閾值。同樣輸出層的節(jié)點(diǎn)值由以下傳遞函數(shù)[式(9)]計(jì)算。

其中，Yj為輸出層第j個(gè)神經(jīng)元節(jié)點(diǎn)；ωij和b2ij為隱藏層與輸出層之間的權(quán)重與閾值。每個(gè)輸出神經(jīng)元的誤差是預(yù)期輸出和實(shí)際輸出之間的差值平方，定義如式(10)所示。

其中，D為實(shí)際的目標(biāo)值，輸出層與隱藏層的誤差信號(hào)分別如式(11)、式(12)所示。

輸出層與隱藏層之間的權(quán)重以及輸入層與隱藏層之間的權(quán)重如式(13)、式(14)求得。

其中，α為學(xué)習(xí)率，權(quán)重的更新如式(15)、式(16)。

其中，k為誤差傳播的迭代數(shù)；β為權(quán)重更新的動(dòng)量參數(shù)。同理可得偏置的更新如式(17)、式(18)。

首先，根據(jù)熱力循環(huán)數(shù)學(xué)模型，在不同的輸入?yún)?shù)下對(duì)跨臨界CO2熱泵進(jìn)行性能模擬分析以獲得大量的驅(qū)動(dòng)數(shù)據(jù)，然后經(jīng)由BP 神經(jīng)網(wǎng)絡(luò)構(gòu)建跨臨界CO2熱泵的熱力學(xué)預(yù)測(cè)模型，訓(xùn)練流程如圖4 所示，初始參數(shù)見表1。通過網(wǎng)絡(luò)的迭代訓(xùn)練得到預(yù)測(cè)模型，具體情況如圖5所示。

采用均方根相對(duì)誤差（RMSE）作為神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型效果的評(píng)判標(biāo)準(zhǔn)[式(19)]。

圖4 BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練流程

表1 神經(jīng)網(wǎng)絡(luò)建模初始參數(shù)

對(duì)氣冷器壓力、COP、壓縮機(jī)耗功以及制冷量各參數(shù)測(cè)試集的預(yù)測(cè)結(jié)果進(jìn)行誤差分析，其結(jié)果見表2。

通過誤差分析，氣冷器壓力的均方根相對(duì)誤差最低，僅為0.65%；COP、壓縮機(jī)耗功與制冷量雖然都存在個(gè)別樣本點(diǎn)的誤差偏大，但均方根相對(duì)誤差也能達(dá)到較高的精度。因此，該神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型能滿足預(yù)測(cè)精度要求，且沒有產(chǎn)生過擬合現(xiàn)象，因此其泛化能力較高。

表2 輸出參數(shù)均方根相對(duì)誤差比較

圖5 輸出參數(shù)測(cè)試數(shù)據(jù)預(yù)測(cè)值與目標(biāo)值的比較

1.3 經(jīng)濟(jì)評(píng)估模型的建立

1.3.1 投資維護(hù)成本

熱泵系統(tǒng)的總投資成本等于各部件的投資成本之和，對(duì)于水泵的投資成本設(shè)定為定值Zhump=240USD，其余每個(gè)部件成本（Zk）由以下成本函數(shù)[式(21)～式(26)]計(jì)算[19]。

其中，D、H、δ分別為水箱的直徑、高度以及保溫層的厚度；pgc與Wcom來源于神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)結(jié)果，換熱器的成本為換熱面積的函數(shù)，換熱面積可通過式(27)求出。

其中，U是整體的換熱系數(shù)，一般估計(jì)在1.5～2.5kW/(m2·K)的范圍內(nèi)[20]；Q和Δtm分別為換熱量與換熱溫差。

在工程經(jīng)濟(jì)學(xué)中，投資成本的評(píng)估通常要進(jìn)行年化，時(shí)間間隔為一年，并且使用資本回收系數(shù)（CRF）獲得相應(yīng)的成本，其可以由式(28)[21]確定。

其中，i和n分別是年利率與系統(tǒng)的運(yùn)行壽命，見表3?？紤]維護(hù)成本，即維護(hù)系數(shù)φ=1.06，總投資維護(hù)成本為式(29)。

表3 經(jīng)濟(jì)參數(shù)

1.3.2 運(yùn)營(yíng)成本

跨臨界CO2熱泵系統(tǒng)的運(yùn)營(yíng)成本主要來源于壓縮機(jī)與水泵的電能消耗，而以上的功耗與水箱的尺寸以及熱損失有關(guān)，水箱結(jié)構(gòu)為圓柱形，其水箱的容積由式(30)求得。

在一天中水箱熱量的損失量由式(31)求得。

其水箱采用的保溫材料為聚氨酯，熱導(dǎo)率為0.022W/(m·K)，對(duì)于一天所需的熱水熱量由式(32)～式(34)求得。

其中，Vdemand為一天熱水所需的體積；ρwater和cm分別為水的密度與比熱容；Δt為水加熱的溫差。在本文中，水泵的功率被設(shè)定為Wpump=657W，因此白天與夜晚的電力成本由式(35)、式(36)求得。

1.3.3 環(huán)境成本

由于熱泵系統(tǒng)的電力消耗所造成的CO2排放，將會(huì)帶來排放罰金，由式(38)求得。

其中，C(CO2)為CO2排放罰金的單價(jià)，通常取在0.022～0.130USD/kg 的范圍[22]，在本文中設(shè)定為0.05USD/kg；m(CO2)為CO2的排放量，由式(39)求得。

其中，μ(CO2)是電網(wǎng)發(fā)電產(chǎn)生CO2的轉(zhuǎn)換系數(shù)，取為0.968kg/(kW·h)。

該經(jīng)濟(jì)模型的建立基于以下假設(shè)：

①系統(tǒng)運(yùn)行穩(wěn)定，且熱水需求量保持不變；

②所有的連接管路沒有熱量損失；

③整個(gè)模型的計(jì)算周期為一天（24h），且環(huán)境溫度保持不變；

④系統(tǒng)優(yōu)先在夜晚（低電價(jià)時(shí)期）運(yùn)行，并裝滿整個(gè)水箱，水箱儲(chǔ)存熱水不夠再在白天運(yùn)行。

2 優(yōu)化方法與案例分析

NSGA-Ⅱ是Deb 等[23]提出的一種基于Goldberg思想的帕累托最優(yōu)概念多目標(biāo)優(yōu)化遺傳算法。與傳統(tǒng)遺傳算法相比，該算法引入了一種基于Pareto和非支配解概念的新選擇機(jī)制。選擇機(jī)制通過個(gè)體之間的支配與非支配關(guān)系對(duì)種群進(jìn)行排序，具有較高非支配等級(jí)的個(gè)體被選擇的概率越高。對(duì)于同一非支配等級(jí)的一系列個(gè)體，擁擠距離越大，被選擇的概率越大。最后，根據(jù)非支配序與擁擠度算子得到了多目標(biāo)優(yōu)化的Pareto解集。一般來說，多目標(biāo)優(yōu)化數(shù)學(xué)模型可以用式(40)～式(43)描述。

其中，x∈Rn為決策變量；Ω為可行域；hi(x)和gj(x)均為約束函數(shù)；f(x)為p個(gè)目標(biāo)組成的目標(biāo)函數(shù)向量。

本文中，考慮的目標(biāo)函數(shù)為最小化的總經(jīng)濟(jì)成本目標(biāo)[式(20)]以及最小化的水箱容積目標(biāo)[式(30)]。并且為了選擇對(duì)目標(biāo)函數(shù)具有實(shí)質(zhì)影響的主要參數(shù)，進(jìn)行了參數(shù)研究，將氣冷器換熱溫差（Δtm）、水箱的直徑（D）、水箱的高度（H）以及水箱保溫層的厚度（δ）作為優(yōu)化的設(shè)計(jì)參數(shù)，其各自的變化范圍見表4。以及本案例的研究是對(duì)于6.45kW的家用CO2熱泵熱水器最優(yōu)化設(shè)計(jì)，其案例數(shù)據(jù)見表5。

表4 設(shè)計(jì)參數(shù)變化范圍

表5 案例數(shù)據(jù)

3 結(jié)果與討論

為了解決優(yōu)化問題，隨機(jī)產(chǎn)生一個(gè)父代種群，計(jì)算所有個(gè)體的目標(biāo)函數(shù)值，再進(jìn)行非支配排序和擁擠度計(jì)算，將非支配序與擁擠度距離組成擁擠度算子進(jìn)行選擇，得到非支配前沿個(gè)體；然后對(duì)種群進(jìn)行交叉、變異操作，產(chǎn)生子代種群；最后將父代與子代種群結(jié)合，選出最優(yōu)秀的種群替換父代種群。具體的優(yōu)化過程如圖6所示。

通過NSGA-Ⅱ算法程序的優(yōu)化求解，得到最終的優(yōu)化結(jié)果，其程序的調(diào)試參數(shù)見表6。研究的跨臨界CO2熱泵系統(tǒng)總經(jīng)濟(jì)成本與水箱容積兩個(gè)目標(biāo)存在明顯的沖突關(guān)系如圖7 表示，從Pareto 前沿圖中可知，水箱容積從0.02m3變化到0.52m3，導(dǎo)致了總年度成本從914USD/a 變化到984USD/a。而圖中A選擇點(diǎn)為僅考慮水箱容積目標(biāo)最小化時(shí)的方案點(diǎn)，其水箱容積只有0.02m3，然而相應(yīng)的總年度成本高達(dá)984USD/a；同樣B點(diǎn)為僅考慮總年度成本目標(biāo)最小化時(shí)的方案點(diǎn)，其總年度成本為914USD/a，而導(dǎo)致的水箱容積為0.52m3。之所以在Pareto 非劣解集中出現(xiàn)了高于案例需求量0.5m3的解方案，是因?yàn)樗浯嬖跓崃康膿p失。

圖6 NSGA-Ⅱ優(yōu)化算法流程

表6 優(yōu)化程序的調(diào)試參數(shù)

圖7 跨臨界CO2熱泵系統(tǒng)多目標(biāo)優(yōu)化的Pareto前沿

圖8 跨臨界CO2熱泵系統(tǒng)多目標(biāo)優(yōu)化的量綱為1化Pareto前沿

對(duì)于多目標(biāo)優(yōu)化問題，從Pareto前沿獲得的所有非支配解均可以被選為最優(yōu)的方案解，然而對(duì)于實(shí)際的工程問題只能選擇一個(gè)最優(yōu)的解方案。近年來，在相應(yīng)領(lǐng)域的研究中，各種決策方法被用于選擇Pareto 前沿中的最優(yōu)解，其中TOPSIS 決策法具有魯棒性強(qiáng)的特點(diǎn)得到廣泛應(yīng)用[24]。TOPSIS 決策法的中心理論是找到距離理想方案點(diǎn)最近且距離非理想方案點(diǎn)最遠(yuǎn)的點(diǎn)[24]，理想點(diǎn)為點(diǎn)A與點(diǎn)B在各自最優(yōu)方向上的交點(diǎn)，如圖8中所示，同樣非理想點(diǎn)為點(diǎn)A與點(diǎn)B在各自最差方向上的交點(diǎn)。

對(duì)于大多數(shù)的多目標(biāo)優(yōu)化問題，其目標(biāo)的物理量綱不同（本文中總年度成本是USD/a，而水箱容積是m3），因此在對(duì)優(yōu)化結(jié)果進(jìn)行TOPSIS方法優(yōu)選之前，還應(yīng)對(duì)目標(biāo)函數(shù)的值進(jìn)行量綱為1化。采用歐式方法[25]對(duì)本文的Pareto 非劣解集進(jìn)行量綱為1化，其結(jié)果如圖8 所示。然后通過計(jì)算得到Pareto前沿上所有解到理想點(diǎn)的距離（d+），以及到非理想點(diǎn)的距離（d-），通過其接近度[d-/(d-+d+)]標(biāo)準(zhǔn)，得到接近度最大的方案解，即綜合考慮了兩個(gè)優(yōu)化目標(biāo)得到的最優(yōu)解，如圖7中所示。其該點(diǎn)的水箱容積與總年度成本分別為0.235m3與958.1USD/a，相應(yīng)的最優(yōu)設(shè)計(jì)參數(shù)見表7，以及相應(yīng)的性能經(jīng)濟(jì)優(yōu)化參數(shù)見表8。

為了分析設(shè)計(jì)參數(shù)對(duì)于優(yōu)化目標(biāo)的影響，得到如圖9 所示的Pareto 種群設(shè)計(jì)參數(shù)分布。從圖中可知，保溫層厚度大致平均分布在65～75mm較小范圍內(nèi)，這是因?yàn)楸貙拥暮穸葘?duì)于水箱的設(shè)計(jì)有最佳的厚度區(qū)域，保溫層過薄保溫效果差，保溫層過厚對(duì)于保溫效果提升不大，反而造成了投資成本快速上升。而水箱的高度分散的分布在0.50～0.83m的范圍內(nèi)，水箱的直徑也均勻分布在0.30～0.93m的范圍內(nèi)，說明了以上兩個(gè)設(shè)計(jì)參數(shù)對(duì)于優(yōu)化目標(biāo)的影響較大，具有直接的作用。而氣冷器的換熱溫差集中在14.5℃的左右，說明了氣冷器換熱溫差對(duì)于優(yōu)化目標(biāo)而言影響較小。

表7 TOPSIS 優(yōu)選點(diǎn)的系統(tǒng)設(shè)計(jì)參數(shù)

表8 TOPSIS優(yōu)選點(diǎn)對(duì)應(yīng)的性能經(jīng)濟(jì)參數(shù)

4 結(jié)論

本文對(duì)跨臨界CO2熱泵系統(tǒng)進(jìn)行了實(shí)用性與經(jīng)濟(jì)性為目標(biāo)的優(yōu)化設(shè)計(jì)分析，通過BP 神經(jīng)網(wǎng)絡(luò)構(gòu)建跨臨界CO2熱泵的預(yù)測(cè)模型，再使用NSGA-Ⅱ優(yōu)化算法對(duì)以上目標(biāo)的跨臨界CO2熱泵系統(tǒng)進(jìn)行優(yōu)化，可獲得以下結(jié)論。

（1）獲得了精度高、無過擬合現(xiàn)象的跨臨界CO2熱泵神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型，能快速、準(zhǔn)確地進(jìn)行數(shù)據(jù)的預(yù)測(cè)。

（2）對(duì)本文中案例進(jìn)行優(yōu)化分析得到Pareto前沿，通過TOPSIS 決策法選出最優(yōu)的設(shè)計(jì)方案為水箱體積0.235m3、總年度成本958.1USD/a。結(jié)果表明，所選的最優(yōu)方案具有占用空間小且總經(jīng)濟(jì)成本小等優(yōu)點(diǎn)。

圖9 Pareto種群設(shè)計(jì)參數(shù)分布

（3）水箱保溫層厚度對(duì)于優(yōu)化目標(biāo)的影響主要集中在一個(gè)較優(yōu)區(qū)域，水箱的直徑與高度對(duì)于優(yōu)化目標(biāo)的影響較大，而氣冷器換熱溫差對(duì)于優(yōu)化目標(biāo)的影響較小。之后的研究可對(duì)設(shè)計(jì)參數(shù)的選取再進(jìn)行優(yōu)化改進(jìn)以獲得影響更有效的參數(shù)。

符號(hào)說明

COP—— 性能系數(shù)

D—— 水箱直徑，m

H—— 水箱高度，m

h—— 焓，J/kg

mr—— CO2流量，kg/s

pgc—— 氣冷器壓力，kPa

Qe—— 制冷量，W

t—— 溫度，℃

Δtm—— 氣冷器換熱溫差，℃

V—— 水箱容積，m3

Z—— 成本，USD

δ—— 水箱保溫層厚度，mm

下角標(biāo)

com—— 壓縮機(jī)

day—— 白天

demand—— 需求

e—— 蒸發(fā)器

env—— 環(huán)境

exv—— 電子膨脹閥

gc—— 氣冷器

in—— 進(jìn)口

ins—— 水箱保溫層

inv—— 投資

loss—— 損失

main—— 維護(hù)

night—— 夜晚

op—— 運(yùn)營(yíng)

out—— 出口

pump—— 水泵

tac—— 總年度成本

water—— 水

wt—— 水箱