中國儲能網(wǎng)訊：
摘 要 為了改善新能源發(fā)電波動對電網(wǎng)的影響，提出了一種光熱-跨臨界壓縮二氧化碳儲能（transcritical compressed carbon dioxide energy storage，TC-CCES）循環(huán)集成熱力系統(tǒng)，采用模塊化機理建模方法，基于能質(zhì)平衡關(guān)系分別建立TC-CCES系統(tǒng)與光熱系統(tǒng)的動態(tài)數(shù)學(xué)模型，獲取TC-CCES系統(tǒng)在儲釋能階段關(guān)鍵參數(shù)的動態(tài)響應(yīng)曲線。研究結(jié)果表明，系統(tǒng)的儲能密度達到28.43 kW/m3，儲能效率與循環(huán)效率分別為58.01%和60.85%，動態(tài)數(shù)學(xué)模型的最大誤差均小于5%。此外，太陽直射輻射變化促使系統(tǒng)熱源溫度變化，而系統(tǒng)負荷對熱源溫度變化非常敏感，熱源溫度升高2.29%，換熱器負荷升高3.36%，而且在某地區(qū)四季典型日冬季比秋季機組負荷低了23.9%。提出的動態(tài)數(shù)學(xué)模型可用于分析太陽能發(fā)電的動態(tài)特性，可為控制系統(tǒng)的設(shè)計提供理論參考。
關(guān)鍵詞 光熱；T-CO2儲能循環(huán)；模型；動態(tài)特性；建模仿真
Research of dynamic characteristics of photothermal coupled transcritical compressed carbon dioxide energy storage cycle
WANG Di CHEN WeiqianSUN Lingfang ZHOU Yunlong
（School of Automation Engineering, Northeast Electric Power University, Jilin 132012, Jilin, China）
Abstract: In order to improve the impact of fluctuations in new energy generation on the power grid, this paper proposes an integrated thermal system called the photothermal coupled transcritical compressed carbon dioxide energy storage (TC-CCES) cycle. The dynamic mathematical models of TC-CCES system and photothermal system were established based on the energy and mass balance relationship, and the dynamic response curves of key parameters of TC-CCES system in energy storage and release stage were obtained. The research results show that the energy storage density of the system reaches 28.43 kW/m3, the energy storage efficiency and cycle efficiency are 58.01% and 60.85% respectively, and the maximum error of the dynamic mathematical model is less than 5%. In addition, changes in direct solar radiation cause the system heat source temperature to change, and the system load is very sensitive to changes in heat source temperature. The heat source temperature increases by 2.29%, and the heat exchanger load increases by 3.36%. In a typical day of four seasons in a certain area, the unit load in winter is 23.9% lower than that in autumn. The dynamic mathematical model presented in this paper can be used to analyze the dynamic characteristics of solar power generation, and lays a theoretical foundation for the design of the control system.
Keywords: photothermal；T-CO2 energy storage cycle；model；dynamic characteristic；modeling and simulation
引 言
面對能源安全與環(huán)保問題，積極發(fā)展可再生能源已成為我國實現(xiàn)能源可持續(xù)發(fā)展的必然選擇[1]。然而，可再生能源發(fā)電的間歇性和不穩(wěn)定性給其有效并網(wǎng)帶來了嚴(yán)峻的挑戰(zhàn)[2]。為了解決該問題，儲能技術(shù)應(yīng)運而生[3]，即在低谷期儲能裝置進行能量儲存，而在用電高峰期儲能則轉(zhuǎn)化為電能輸出，有效地平衡電力需求。
目前，可成熟應(yīng)用的大型儲能方式主要包括抽水蓄能（PHS）和壓縮空氣儲能（CAES）[4]。但PHS因地理條件的限制，難以廣泛應(yīng)用；而傳統(tǒng)CAES依賴燃氣輪機技術(shù)、化石燃料的補充燃燒以及大規(guī)模適合壓縮空氣儲存的天然洞穴，同時也存在污染環(huán)境的排放問題[5]。與空氣相比，CO2具有較低的黏度、較大的擴散系數(shù)和密度，而且臨界點更低，故壓縮CO2儲能系統(tǒng)結(jié)構(gòu)更加緊湊、儲能密度更高[6]。眾多學(xué)者對壓縮CO2儲能系統(tǒng)進行了廣泛研究。Ahmadi等[7]研究了一種新型壓縮二氧化碳儲能系統(tǒng)，并根據(jù)數(shù)學(xué)模型分析了不同熱力學(xué)參數(shù)對系統(tǒng)性能的影響。Kim等[8]從系統(tǒng)設(shè)計角度提出了跨臨界壓縮二氧化碳儲能（TC-CCES）的概念系統(tǒng)，該系統(tǒng)通過等溫壓縮和等溫膨脹過程降低壓縮機單元功耗，提高膨脹機單元輸出功率，從而顯著提升系統(tǒng)的循環(huán)效率。Hao等[9]提出了一種基于熱泵的TC-CCES系統(tǒng)，該系統(tǒng)利用熱泵提升膨脹機初始溫度，無須補充燃燒，使系統(tǒng)具有較好的綜合熱力性能。Zhang等[10]對基于布雷頓循環(huán)的TC-CCES系統(tǒng)進行熱力學(xué)性能分析，結(jié)果表明，與先進絕熱CAES系統(tǒng)相比，該TC-CCES系統(tǒng)具有更高的系統(tǒng)熱力性能。Liu等[11]提出了雙地下儲氣室的CCES系統(tǒng)，并對比分析了S-CO2和T-CO2儲能系統(tǒng)的熱力學(xué)特性，得出CCES系統(tǒng)具有更高能量存儲密度。Zhang等[12]提出了一種采用填充床蓄熱器的TC-CCES系統(tǒng)，研究結(jié)果表明該系統(tǒng)的往返效率比液態(tài)CO2系統(tǒng)高4.05%，具有良好的應(yīng)用前景。吳毅等[13]提出了一種新型TC-CCES系統(tǒng)，將CO2進行液態(tài)存儲，實現(xiàn)冷能和熱能的循環(huán)使用，研究表明該系統(tǒng)具有較好的熱力性能和較高的儲能密度。Liu等[14]提出了帶有冷回?zé)崞鞯囊簯B(tài)CO2儲能系統(tǒng)，通過系統(tǒng)優(yōu)化和改進提高了性能。
Manikumar等[15]提出了一個由太陽能集熱系統(tǒng)驅(qū)動的跨臨界CO2動力循環(huán)，并研究了聚光面積、聚光功率和CO2流量對集熱器熱性能的影響。Yamaguchi等[16]研究了不同季節(jié)太陽能輸入對CO2電力系統(tǒng)性能的影響，研究表明太陽能集熱器是影響系統(tǒng)性能的決定性因素。Alzahrani等[17]設(shè)計了一種兩級壓縮的太陽能跨臨界CO2發(fā)電循環(huán)，但其儲能效率僅為34%。Zhang等[18]探討了太陽能發(fā)電塔內(nèi)熔鹽流速控制系統(tǒng)，并通過實驗獲得了系統(tǒng)的動態(tài)特性曲線。Singh等[19]通過模擬冬夏兩季太陽輻射輸入對CO2循環(huán)的影響研究了太陽能直接加熱CO2循環(huán)的熱力學(xué)特性。一方面太陽能的間歇性和短時性會對發(fā)電系統(tǒng)產(chǎn)生負面影響，另一方面動態(tài)特性是控制系統(tǒng)設(shè)計的基礎(chǔ)，因此研究系統(tǒng)的動態(tài)特性具有重要意義。
綜上所述，盡管國內(nèi)外研究學(xué)者已經(jīng)提出了一些光熱耦合壓縮CO2發(fā)電循環(huán)方案，但尚未見有深入研究TC-CCES系統(tǒng)與光熱耦合后的動態(tài)特性及系統(tǒng)性能。本研究提出一種光熱-跨臨界二氧化碳儲能循環(huán)系統(tǒng)的復(fù)合方案，采用建模仿真的方法構(gòu)建復(fù)合系統(tǒng)的動態(tài)數(shù)學(xué)模型，并通過不同工況模擬獲得TC-CCES系統(tǒng)儲能和釋能階段關(guān)鍵參數(shù)的動態(tài)特性。以期深入理解TC-CCES系統(tǒng)的動態(tài)特性及控制系統(tǒng)的設(shè)計提供理論參考。
1 研究對象介紹
光熱耦合TC-CCES循環(huán)集成熱力系統(tǒng)結(jié)構(gòu)圖如圖1(a)所示，圖1(b)為光熱系統(tǒng)接收器的俯瞰結(jié)構(gòu)圖。TC-CCES系統(tǒng)由壓縮機、透平、換熱器、儲冷/儲熱罐等部件構(gòu)成。復(fù)合系統(tǒng)的T-s圖如圖2所示。
圖1 系統(tǒng)結(jié)構(gòu)圖
Fig.1 System structure drawing
圖2 光熱耦合TC-CCES循環(huán)系統(tǒng)的T-s圖
Fig.2 T-s diagram of photothermal coupled TC-CCES cycle system
儲能過程：從低壓罐LPT中流出的液態(tài)CO2流經(jīng)節(jié)流閥TV（1→2），在蒸發(fā)器EV中蒸發(fā)后成為氣態(tài)（2→3），然后CO2經(jīng)壓縮機壓縮（3→4，5→6），在壓縮過程中同時產(chǎn)生大量的熱，壓縮熱通過間冷器被來自冷罐的水吸收，吸收壓縮熱后的熱水儲存在儲熱罐中（4→5，6→7），經(jīng)冷凝器高溫超臨界CO2被冷卻為液態(tài)CO2（7→8），這種高壓液態(tài)CO2儲存在高壓罐HPT中。
釋能過程：從高壓罐HPT中流出的高壓液態(tài)CO2經(jīng)低溫泵升壓（9→10），儲存在儲熱罐中的高溫儲熱介質(zhì)進入加熱器加熱高壓CO2（10→11，13→14），預(yù)熱的CO2進入換熱器與吸收太陽能的蓄熱介質(zhì)交換熱量（11→12），然后高溫高壓CO2進入透平膨脹做功。透平排氣經(jīng)冷凝器冷卻至液態(tài)（15→16），儲存在低壓罐LTP中。
2 數(shù)學(xué)模型
本研究建模平臺采用Matlab/Simulink，CO2物性參數(shù)通過調(diào)用工質(zhì)物性查詢軟件REFPROP獲得。
2.1 TC-CCES系統(tǒng)模型
2.1.1 壓縮機模型
壓縮機是T-CO2儲能循環(huán)系統(tǒng)的關(guān)鍵部件之一，考慮此系統(tǒng)所需壓力比較小，為使結(jié)構(gòu)緊湊，選用單級離心式壓縮機。由于CO2的密度、熱導(dǎo)率及比熱容等參數(shù)在臨界點附近會發(fā)生突變，導(dǎo)致對壓縮機性能要求更高，在建模中使用壓縮機性能曲線能更好地反映壓縮機的性能特性[20]。
圖3 壓縮機性能
Fig.3 Compressor performance
結(jié)合上述壓縮機性能圖，壓縮機出口工質(zhì)焓的求取如式（4）所示：
2.1.2 轉(zhuǎn)軸模型
轉(zhuǎn)軸數(shù)學(xué)模型主要用以計算轉(zhuǎn)速的變化，根據(jù)角動量平衡方程，轉(zhuǎn)軸動態(tài)模型如下[21]：
式中，ω為角速度，rad/s；t為時間，s；J為轉(zhuǎn)動慣量，kg?m2，是整個轉(zhuǎn)子系統(tǒng)的轉(zhuǎn)動慣量；Pec為電機輸出功率，W；Pc為壓縮機耗功，W。
2.1.3 透平模型
徑向透平理想效率與速比的關(guān)系可由式(9)~式(11)求得，如圖4所示。
圖4 徑向透平理想效率與速比的關(guān)系
Fig.4 Efficiency of an ideal radial turbine as a function of velocity ratio
結(jié)合透平進口參數(shù)，透平的出口工質(zhì)焓的求取如式（12）所示：
2.1.4 換熱設(shè)備模型
T-CO2儲能循環(huán)系統(tǒng)的換熱設(shè)備包括T-CO2-熔鹽換熱器、熔鹽-T-CO2換熱器以及預(yù)冷器。換熱器、預(yù)冷器作為T-CO2儲能循環(huán)系統(tǒng)的主要換熱部件，對整個系統(tǒng)的效率有著很大的影響。換熱設(shè)備的動態(tài)數(shù)學(xué)模型主要計算冷端和熱端的換熱量以及冷端和熱端的溫度變化[22]。
根據(jù)能量守恒方程，冷、熱流體比焓隨時間的變化可表示為[22]：
2.1.5 泵模型
泵主要是用來升高工質(zhì)的壓力，使工質(zhì)在管道內(nèi)流動的過程中可以達到一定的流速。由于泵的流動響應(yīng)速度快且復(fù)雜，為了便于計算，本研究采用基于經(jīng)驗數(shù)據(jù)的簡化靜態(tài)模型[24]。泵的性能如圖5所示。
圖5 泵的性能
Fig.5 Pump performance
泵的流量如式（20）所示[24]：
2.1.6 儲熱/儲冷罐模型
儲熱/儲冷罐模型的研究主要是分析罐內(nèi)換熱流體的溫度變化，其質(zhì)量守恒與能量守恒方程如下[25]：
2.1.7 節(jié)流閥模型
T-CO2儲能循環(huán)系統(tǒng)中節(jié)流閥主要起穩(wěn)壓作用。
因工質(zhì)流動快，其閥門內(nèi)的減壓過程可視作絕熱等焓過程，如式(26)所示[26]：
式中，hin,value為節(jié)流閥的入口焓值，kJ/kg；hout,value為節(jié)流閥的出口焓值，kJ/kg。
2.2 光熱模型
2.2.1 定日鏡場
定日鏡場主要由若干定日鏡組成，通過自動跟蹤裝置將太陽能匯集到太陽能接收器中。
定日鏡接收到的太陽總輻射能與定日鏡面積和太陽直接輻射量（direct normal irradiance，DNI）有關(guān)，如式(27)所示[27]：
式中，PDNI為定日鏡場接收到的太陽總輻射能，kW；DNI為太陽直接輻射量，kW/m2；Amirror為定日鏡場的面積，m2。
接收器上的太陽輻射功率主要受定日鏡場效率影響，如式（28）所示[27]：
式中，圖片為接收器上的太陽輻射功率，kW；ηhel為定日鏡場效率，根據(jù)文獻[28]計算。
2.2.2 接收器模型
塔式太陽能發(fā)電站的集熱子系統(tǒng)主要由接收器和接收塔兩部分組成。接收器內(nèi)傳熱工質(zhì)為60%硝酸鈉與40%硝酸鉀的混合熔鹽。接收器由很多吸熱管組成，吸熱管材質(zhì)為316H不銹鋼，具體尺寸參數(shù)見表1[29]。
表1 接收器系統(tǒng)參數(shù)
Table 1 Receiver system parameter
對接收器的動態(tài)數(shù)學(xué)模型進行簡化，并做出如下假設(shè)[28]：
（1） 假設(shè)接收器每塊面板上每一點接收到的太陽能流密度相同[29]；
（2） 接收器管壁傳熱過程中局部達到熱力學(xué)平衡。
在這些假設(shè)的基礎(chǔ)上，可以認為位于同一面板的每條管道都是等效的，并且每條管道的出入口熔鹽溫度是相等的。同時，由于每條熔鹽回路分別經(jīng)過12個面板，每條回路的模型是將12個面板連接起來，接收器的最終模型是兩個回路的疊加。因此，面板上的管道建模就相當(dāng)于接收器模型的建模。
對接收器管道的管壁建立熱平衡方程，管壁能量平衡方程為[29]：
2.2.3 光熱換熱裝置模型
2.3 系統(tǒng)儲能效率
儲能效率ηES用來衡量一個儲能周期內(nèi)存儲的電能對整個儲能周期外部產(chǎn)能的凈貢獻能力，將ηES定義為系統(tǒng)能量釋放過程的凈電能生產(chǎn)與儲存過程的凈電能消耗之比，如式（37）所示[31]：
2.4 系統(tǒng)循環(huán)效率
循環(huán)效率ηCY用來衡量系統(tǒng)循環(huán)運行的效率，將ηCY定義為發(fā)電系統(tǒng)所做功與儲能系統(tǒng)所消耗功之比，如式（38）所示：
2.5 系統(tǒng)儲能密度
系統(tǒng)儲能密度定義為：
3 模型驗證
3.1 TC-CCES模型驗證
3.1.1 靜態(tài)驗證
本研究基于文獻[32]所搭建的循環(huán)系統(tǒng)對模型進行靜態(tài)驗證，在相同的輸入條件下對比分析文獻與仿真結(jié)果，從而驗證模型的可靠性。對比結(jié)果見表2。
表2 TC-CCES模型靜態(tài)驗證
Table 2 Steady-state verification of TC-CCES model
注：表中所列參數(shù)均為出口參數(shù)。
由表2可知，系統(tǒng)各參數(shù)的相對誤差均小于2%，其中溫度的最大誤差為1.31%，壓力的最大誤差為0.09%，計算誤差較小。因此所建模型具有良好的精度，滿足系統(tǒng)動態(tài)特性的分析需求。
3.1.2 動態(tài)驗證
采用文獻[33]對換熱設(shè)備進行動態(tài)驗證，對換熱設(shè)備冷端入口溫度施加正弦信號，其動態(tài)響應(yīng)過程與文獻[33]基本一致，其各項參數(shù)相對誤差均在5%以內(nèi)。其具體對比結(jié)果如圖6所示。
圖6 換熱設(shè)備出口溫度動態(tài)驗證
Fig.6 Dynamic verification of outlet temperature of heat exchange equipment
3.2 光熱系統(tǒng)模型驗證
3.2.1 靜態(tài)驗證
基于文獻[34]所搭建的循環(huán)系統(tǒng)對所建立模型進行靜態(tài)驗證，把位于34°52′18″N、116°50′3″W的Solar Two電站的運行參數(shù)作為參考對象，在相同的輸入條件下對比分析文獻與仿真結(jié)果，從而驗證模型的可靠性。對比結(jié)果見表3。
表3 光熱模型靜態(tài)驗證
Table 3 Steady-state verification of photothermal model
由表3可見，溫度的最大誤差在1%以內(nèi)，計算誤差較小。因此所建模型具有良好的精度，滿足系統(tǒng)動態(tài)特性的分析需求。
3.2.2 動態(tài)驗證
以文獻[35]中的相關(guān)數(shù)據(jù)對光熱系統(tǒng)接收器展開動態(tài)驗證，針對文獻所構(gòu)建系統(tǒng)，參照其系統(tǒng)運行中各部件的運行數(shù)據(jù)，保證接收器與文獻的入口條件相同，通過出口條件的對比驗證模型的準(zhǔn)確性。
首先，在接收器穩(wěn)定運行情況下，將太陽輻射能階躍增大5%后，第一塊板出口溫度從323.5℃上升到324.7℃，溫度上升1.2℃，經(jīng)過約800 s后溫度趨于穩(wěn)定；第十二塊板升溫幅度較大，出口溫度從566.44℃上升到579.39℃，溫度上升12.95℃，經(jīng)過約2500 s后溫度趨于穩(wěn)定。將兩塊板的出口參數(shù)與文獻對比，其相對誤差均小于5%。接收板熔鹽出口溫度動態(tài)響應(yīng)變化曲線如圖7所示。
圖7 太陽輻射能階躍擾動對太陽能接收器出口的影響
Fig.7 Variation in solar radiation energy step perturbation on solar receiver outlet
其次，在接收器穩(wěn)定運行情況下，改變?nèi)埯}入口流速，使其階躍增大5%后，第一塊板出口溫度從323.5℃降低到322.4℃，溫度降低1.1℃，經(jīng)過約900 s后溫度趨于穩(wěn)定；第十二塊板出口溫度從565.76℃降低到555.43℃，溫度降低10.33℃，經(jīng)過約2600 s后溫度趨于穩(wěn)定。將兩塊板的出口參數(shù)與文獻對比，其相對誤差均小于5%。接收板熔鹽出口溫度動態(tài)響應(yīng)變化曲線如圖8所示。
圖8 流量階躍擾動對太陽能接收器出口的影響
Fig.8 Effect of flow step perturbation on solar receiver outlet
4 動態(tài)特性分析
基于已驗證模型，根據(jù)系統(tǒng)設(shè)計結(jié)構(gòu)搭建整個系統(tǒng)的仿真模型。表4和表5分別為光熱耦合TC-CCES系統(tǒng)的設(shè)計參數(shù)和系統(tǒng)性能指標(biāo)計算結(jié)果。由表5可得，系統(tǒng)儲能效率、循環(huán)效率、儲能密度分別為58.01%、60.85%、28.43 kW/m3，具有較高的儲能性能。
表4 系統(tǒng)設(shè)計參數(shù)
Table 4 System design parameter
表5 系統(tǒng)性能指標(biāo)計算結(jié)果
Table 5 Calculation result of system performance indicators
4.1 單日DNI對系統(tǒng)運行參數(shù)的影響
首先選取單日DNI變化，如圖9所示，其中最大輻射值點在12時，而且DNI在17時后接近0。
圖9 單日太陽直接輻射量
Fig.9 Direct normal irradiance in single day
選取某日從早上7:00到晚上17:00的太陽直接輻射變化，如圖9中局部放大圖所示。在開始運行前必須避免熔鹽在管道內(nèi)凝固的現(xiàn)象發(fā)生，故而需要先預(yù)熱接收器，用以保證管壁的溫度高于573.15 K，此后開始正常運行。本研究設(shè)定整體模型的仿真時間為7:30~16:30，共9 h。而后對整體模型進行仿真，得到主要輸出變量的變化情況，如圖10所示。接收器與換熱器出口溫度變化如圖10（a）所示，DNI變化促使系統(tǒng)熱源溫度變化，熔鹽出口溫度直接影響換熱器冷端出口溫度即透平入口溫度。換熱器負荷變化如圖10（b）所示，對比分析換熱與集熱時間，在所選取的時間段里換熱器負荷剛呈減緩趨勢且仍在變化，故換熱時間比集熱時間更長。DNI波動過程中，高于平均值的能量被儲存在儲熱罐中，從而維持換熱器中熔鹽與水的換熱量的穩(wěn)定，接收器停止工作后罐內(nèi)儲存的熱量可以用來繼續(xù)加熱水，從而使換熱器的工作時間得到延長。
圖10 接收器-換熱器出口參數(shù)變化曲線
Fig.10 Receiver-heat exchanger outlet parameter change curve
4.2 不同DNI對系統(tǒng)運行參數(shù)的影響
本研究設(shè)定的不同工況下DNI特性曲線如圖11所示。其中，C工況為實際收集的DNI的變化情況，設(shè)其為標(biāo)準(zhǔn)的DNI；B工況代表1.1倍的標(biāo)準(zhǔn)DNI；A工況代表1.2倍的標(biāo)準(zhǔn)DNI。
圖11 不同工況下DNI特性曲線
Fig.11 DNI characteristic curves under different working conditions
選取某日從早上7:00到晚上17:00的3種工況下太陽輻射強度變化，如圖12所示。
圖12 不同工況下DNI變化
Fig.12 DNI changes under different working conditions
對比3種工況，系統(tǒng)所接收到的太陽熱能不同，其中A工況下系統(tǒng)接收的最多，故而A工況擁有最長的加熱時間。隨著DNI比例的逐漸降低，接收器出口熔鹽溫度緩慢下降[圖13（a）]。從9：00到12:00，A工況的DNI從1097.51 W/m2上升到1229.89 W/m2，升高了12.06%；熔鹽出口溫度從887.95 K上升到908.34 K，升高了2.29%；換熱器負荷從16.36 MW上升到16.91 MW，升高了3.36%。換熱器出口參數(shù)動態(tài)響應(yīng)主要由熱端入口溫度變化引起，熱端入口溫度降低導(dǎo)致冷熱兩側(cè)出口溫度均降低[圖13（b）]。故換熱量減少，換熱負荷也隨之降低[圖13（c）]。
圖13 不同工況下接收器-換熱器出口參數(shù)變化曲線
Fig.13 Change curve of outlet parameters of receiver and heat exchanger under different working conditions
4.3 四季不同DNI對系統(tǒng)運行參數(shù)的影響
本研究選取某地一年四季典型日太陽輻射，如圖14所示。
圖14 DNI四季變化曲線
Fig.14 DNI variation curve in four seasons
為研究機組關(guān)鍵參數(shù)的動態(tài)特性曲線，在配備儲冷/儲熱罐的情況下對系統(tǒng)進行模擬，如圖15所示。由圖可知，換熱器的負荷變化秋季最高，冬季最低，機組運行更加穩(wěn)定。這些結(jié)果主要是由于冬季太陽輻射較低造成的。以夏季為例，換熱器負荷在9：00前下降，是由于這段時間太陽輻射較低，機組主要依靠熱罐中儲存的能量。太陽輻射繼續(xù)增加，換熱器負荷逐漸上升，在約18：00時達到峰值。隨后，由于太陽輻射逐漸減少，太陽能集熱器熔鹽出口溫度下降，換熱器負荷開始緩慢下降。
圖15 換熱器負荷變化曲線
Fig.15 Heat exchanger load change curve
由于太陽輻射冬季較低、秋季較高，造成換熱器冷端出口溫度變化秋季最高、冬季最低，故而透平入口溫度變化趨勢與換熱器冷端出口溫度相同，如圖16（a）所示。在四季典型日不同DNI下，透平入口溫度在9：00前呈降低趨勢，透平入口溫度降低導(dǎo)致等熵出口焓降低，而且入口焓降低頻率大于出口焓，流量不變，故導(dǎo)致透平做功功率也呈降低趨勢，如圖16（b）所示。太陽輻射繼續(xù)增加，換熱器冷端出口溫度逐漸上升，透平入口溫度也逐漸上升，在18：00時達到峰值。隨后，由于太陽輻射逐漸減少，太陽能接收器熔鹽出口溫度下降，換熱器熱端入口溫度開始緩慢下降，熱端與冷端出口溫度均開始緩慢下降，繼而透平入口溫度也開始緩慢下降，流量不變，等熵出口焓降低，則透平做功功率也開始下降。綜上所述，秋季機組負荷維持在2.34 MW左右，冬季機組負荷維持在1.78 MW左右，故而在四季典型日該地區(qū)冬季負荷比秋季低23.9%。
圖16 不同DNI透平出口參數(shù)變化曲線
Fig.16 Variation curves of turbine outlet parameters of different DNI
結(jié)合以上研究，TC-CCES系統(tǒng)在不同工況下儲能密度、循環(huán)效率對比仿真結(jié)果見表6。
表6 不同工況下（四季）系統(tǒng)性能指標(biāo)
Table 6 System performance index under different working conditions （four seasons）
圖17為四季不同DNI工況下TC-CCES系統(tǒng)循環(huán)效率變化曲線。儲能循環(huán)系統(tǒng)運行效率是衡量太陽能發(fā)電性能的重要指標(biāo)，綜合上述分析可知，本研究選取的四季典型日冬季太陽輻射最低，秋季太陽輻射最高，故而系統(tǒng)循環(huán)效率春夏兩季低于秋季、高于冬季，系統(tǒng)循環(huán)效率約為45.53%，如圖17所示。
圖17 系統(tǒng)循環(huán)效率變化曲線
Fig.17 Variation curves of system cycle efficiency
圖18為壓縮機出口參數(shù)響應(yīng)曲線。以春季為例，9：00前太陽輻射較低，換熱器換熱量減少，使得透平入口溫度在9：00前下降，導(dǎo)致泵入口壓力升高，繼而導(dǎo)致壓縮機流量升高，因為壓縮機出口焓降低頻率與入口焓基本一致，故壓縮機耗功有些許升高；壓縮機入口溫度升高，導(dǎo)致入口熵降低，入口密度上升，等熵出口焓降低，故出口壓力上升；又因為出口焓降低頻率大于出口壓力上升頻率，故出口溫度呈現(xiàn)降低趨勢。太陽輻射繼續(xù)增加，壓縮機出口溫度逐漸上升，在約19：00時達到峰值。隨后，由于太陽輻射逐漸減少，太陽能接收器熔鹽出口溫度下降，換熱器熱端入口溫度開始緩慢下降，換熱器換熱量緩慢下降，繼而透平入口溫度也開始緩慢下降，泵入口壓力升高，導(dǎo)致壓縮機流量升高，入口熵降低，入口密度上升，等熵出口焓降低，故出口壓力開始上升，則壓縮機出口溫度開始下降，耗功升高。
圖18 壓縮機出口參數(shù)響應(yīng)曲線
Fig.18 Response curve of compressor outlet parameters
圖19為透平出口參數(shù)響應(yīng)曲線。以春季為例，9：00前太陽輻射較低，造成透平入口溫度在9：00前下降，等熵出口焓降低，而且入口焓降低頻率大于出口焓，流量不變，故導(dǎo)致透平做功功率也呈降低趨勢；透平入口溫度降低，流量不變，入口熵降低，入口密度上升，等熵出口焓降低，故導(dǎo)致出口壓力上升；又因為出口焓降低頻率大于出口壓力上升頻率，故出口溫度呈現(xiàn)降低趨勢。太陽輻射繼續(xù)增加，透平出口溫度逐漸上升，在約19：00時達到峰值。隨后，由于太陽輻射逐漸減少，太陽能接收器熔鹽出口溫度下降，換熱器熱端入口溫度開始緩慢下降，熱端與冷端出口溫度均開始緩慢下降，繼而透平入口溫度也開始緩慢下降，流量不變，入口熵降低，入口密度上升，等熵出口焓降低，故出口壓力開始上升，則透平出口溫度下降，同時做功下降。
圖19 透平出口參數(shù)響應(yīng)曲線
Fig.19 Response curve of turbine outlet parameters
5 結(jié)論
（1）本研究構(gòu)建了光熱-跨臨界二氧化碳儲能系統(tǒng)模型，研究結(jié)果表明，系統(tǒng)的儲能效率、循環(huán)效率分別為58.01%、60.85%，系統(tǒng)儲能密度為28.43 kW/m3，具有較好的系統(tǒng)熱力性能和研發(fā)前景。
（2）對機理模型進行了變工況研究，對比3種工況，其中A工況加熱溫度最高，最高可使透平入口溫度達到535.83 K。而且在某地四季典型日，該地區(qū)冬季負荷比秋季低23.9%。
（3）跨臨界CO2儲能系統(tǒng)處于熱力學(xué)研究的基礎(chǔ)階段，本研究結(jié)果可為TC-CCES系統(tǒng)動態(tài)特性的認識提供理論參考，并為光熱-TC-CCES循環(huán)的控制系統(tǒng)設(shè)計奠定理論基礎(chǔ)。
