Study on the Performance of a Solar Heating System with Seasonal and Cascade Thermal-Energy Storage

: Seasonal solar thermal-energy storage systems used for space heating applications is a promising technology to reduce greenhouse gas emissions. A novel solar heating system with seasonal and cascade thermal-energy storage based on zeolite water is proposed in this study. The system’s efﬁciency is improved through cascade storage and the release of solar energy. The energy storage density is improved through the deep coupling of daily energy storage and cross-seasonal energy storage. A mathematical model of the system-performance analysis is established. The system performances in the non-heating and heating seasons and throughout the year are analyzed by considering the Chifeng City of China as an application case. The results indicate that the average collection efﬁciency of the proposed system is 2.88% higher in the non-heating season and 7.4% higher in the heating season than that of the reference system. Furthermore, the utilization efﬁciency of the proposed system is 37.16%, which is 3.26% higher than that of the reference system. Further, the proposed system has a supply heat of 2135 GJ in the heating season, which is 9.66% higher than the reference system. This study provides a solution for the highly efﬁcient solar energy utilization for large-scale space-heating applications.


Introduction
Renewable energy can contribute to reducing primary energy consumption and environmental impacts. However, the use of renewable energy is limited by its intermittency and volatility. Thermal energy storage (TES) technology can be used to address the mismatch between renewable-energy supply and fluctuating demand, and greatly increase the efficiency of renewable energy [1].
In the past decade, TES systems have been extensively applied and studied [2][3][4][5][6][7][8]. TES technology can be classified into sensible heat storage, latent heat storage, and thermochemical heat storage according to different heat-storage mechanisms. It can be further classified into short-term and long-term heat storage based on the heat-storage period. The seasonal TES (STES) system can solve the time-and space-discrepancy problems of solar energy utilization, and thus, it has been used in space heating [9][10][11][12][13][14][15][16].
STES systems possess many patterns, such as seasonal hot water tank TES, seasonal underground TES, and seasonal thermochemical heat storage. Solar heating systems with STES have generally been adopted and investigated because they can solve the time and location mismatch between solar energy and heat load. Kulkarni [17] investigated the thermal responses of seasonal heat storage and solar thermal systems. The results indicated an enormous waste of exergy in the system based on a natural gas boiler, compared with the systems based on heat pumps with seasonal heat storage and solar thermal collectors with an efficiency of 0.47. Zhao [18] simulated the feasibility of seasonal underground thermal-energy storage in Tianjin, China. The results indicated a thermal recovery ratio of less than 67% and different temperature distributions under the three modes. Dada [19] simulated a novel seasonal heat-storage system dedicated to storing heat in the ground during the hot period and then recovering it during the cold period. This approach led to a remarkable increase in temperature levels as well as heat-energy recovery efficiency, which increased from 41% to 71%. Liu [20] investigated the thermal equilibrium of a solar seasonal storage system coupled with a ground-source heat pump. The results showed that the solar-energy utilization efficiency was 50.2%. The feasibility of seasonal thermochemical heat storage has been investigated [21][22][23]. It was found that STES overcomes the drawback of the intermittency of solar energy and contributes to storing heat in summer for use in winter. Therefore, the seasonal solar TES (SSTES) system is a promising technology that minimizes greenhouse gas emissions by harnessing solar energy for space-heating applications.
Thermochemical heat-storage technology, which has the advantages of high energy density and long duration, is recognized as a promising TES technology for short-and long-term storage applications. Sorption heat storage, a classical thermochemical heatstorage technology, is divided into different types according to various sorption materials such as zeolite, silica gel, and hygroscopic salts. Zeolite 13X has many advantages, such as sufficient heat of adsorption as a hydrophilic material, high thermal-storage density, tolerance to water, no-swelling property, utilization of low-temperature thermal energy, and environmental protection. Therefore, zeolite 13X has been adopted and investigated for many adsorption thermal-storage systems [24][25][26][27].
However, previous studies on STES have always focused on the recovery quantity of heat storage without the recovery quality of heat storage. The temperature recovery must be less than 100% in sensible and latent heat storage. However, temperature recovery can be more than 100% in thermochemical heat-storage technology, for example, adsorption heat pumps [25][26][27].
This study aims to build a novel solar heating system with seasonal and cascade thermal-energy storage based on zeolite water and verify the feasibility of this system by selecting the case of a district heating network in Chifeng, China. Furthermore, sensitivity analyses of the proposed system are analyzed in detail. This study provides a solution for more efficient solar energy utilization for large-scale space-heating applications. Figure 1 illustrates a solar heating system with seasonal and cascade TES based on zeolite water-comprised of solar thermal collectors, a zeolite adsorption reactor, a condenser (C), water tanks (W1 and W2), pumps (P1-P4), heat exchangers (EXC1-EXC3) and valves (V0-V10-3). The zeolite adsorption reactor consisted of many reaction units, and each reaction unit could run independently. The zeolite adsorption reactor was filled with zeolite and was set by embedded heat exchangers.

System Description
In the non-heating season, the solar heat-collection, zeolite adsorption, and mediumtemperature TES subsystems work, while the low-temperature TES, energy release, and heating subsystems do not work. Based on the desorption dynamic of zeolite, the desorption reaction occurred quickly when the zeolite moisture content was higher and did not occur easily when the zeolite moisture content was lower. To improve the efficiency of solar-energy utilization, the system adopted solar-energy cascade heat storage and zeolite cascade desorption. Medium-temperature heat was collected and stored in tank W2 when the solar radiation intensity was lower, and the stored heat was used to desorb zeolite with a higher moisture content. Here, tank W2 also played a role in regulating the desorption time. In this operating mode, valves V2-1, V2-2, V4-1, V4-2, and V0 were opened. High-temperature heat was collected and used to desorb zeolite with a lower moisture content when the solar radiation intensity was higher. Moreover, in this operating mode, valves V3-1, V3-2, V4-1, V4-2, and V0 were opened. Thus, the amount of solar energy collected via all operating modes was eventually stored in the zeolite adsorption reactor. In the reactor, an endothermic desorption reaction occurred, and the wet zeolite became dry to complete the regeneration of zeolite; at the same time, water vapor was generated during the endothermic reaction, which was removed from the reactor and condensed in condenser C. In the non-heating season, the solar heat-collection, zeolite adsorption, and mediumtemperature TES subsystems work, while the low-temperature TES, energy release, and heating subsystems do not work. Based on the desorption dynamic of zeolite, the desorption reaction occurred quickly when the zeolite moisture content was higher and did not occur easily when the zeolite moisture content was lower. To improve the efficiency of solar-energy utilization, the system adopted solar-energy cascade heat storage and zeolite cascade desorption. Medium-temperature heat was collected and stored in tank W2 when the solar radiation intensity was lower, and the stored heat was used to desorb zeolite with a higher moisture content. Here, tank W2 also played a role in regulating the desorption time. In this operating mode, valves V2-1, V2-2, V4-1, V4-2, and V0 were opened. High-temperature heat was collected and used to desorb zeolite with a lower moisture content when the solar radiation intensity was higher. Moreover, in this operating mode, valves V3-1, V3-2, V4-1, V4-2, and V0 were opened. Thus, the amount of solar energy collected via all operating modes was eventually stored in the zeolite adsorption reactor. In the reactor, an endothermic desorption reaction occurred, and the wet zeolite became dry to complete the regeneration of zeolite; at the same time, water vapor was generated during the endothermic reaction, which was removed from the reactor and condensed in condenser C.
In the heating season, the system adopted cascade heat storage and release to improve energy efficiency. During the charging process, medium-temperature heat was stored in tank W2 when the solar radiation intensity was higher, and valves V2-1 and V2-2 were opened in the operating mode. Low-temperature heat was stored in tank W1 when the solar radiation intensity was low, and valves V1-1 and V1-2 were opened in operating mode. During the discharging process, the heating backwater was heated sequentially by stored low-temperature and medium-temperature heat. It was then used for heating when the heating load was low or when solar energy was abundant. Moreover, V5-1, V5-2, V6-1, V6-2, V7-1, V7-2, V8-1, V8-2, and V10-1 were opened in this operating mode. The back- In the heating season, the system adopted cascade heat storage and release to improve energy efficiency. During the charging process, medium-temperature heat was stored in tank W2 when the solar radiation intensity was higher, and valves V2-1 and V2-2 were opened in the operating mode. Low-temperature heat was stored in tank W1 when the solar radiation intensity was low, and valves V1-1 and V1-2 were opened in operating mode. During the discharging process, the heating backwater was heated sequentially by stored low-temperature and medium-temperature heat. It was then used for heating when the heating load was low or when solar energy was abundant. Moreover, V5-1, V5-2, V6-1, V6-2, V7-1, V7-2, V8-1, V8-2, and V10-1 were opened in this operating mode. The backwater was heated by the stored low-temperature, medium-temperature, and zeolite adsorption heat alternatively. Subsequently, it was used for heating when the heating load was high, or the solar energy was insufficient. An exothermic adsorption reaction occurred in the reactor, and the dry zeolite became wet. Furthermore, V5-1, V5-2, V6-1, V6-2, V7-1, V7-2, V8-1, V8-2, V9-1, and V9-2 were opened in this operating mode.

Mathematical Model
To simplify the proposed model, the following assumptions were considered: (1) The radiation loss of any component and pipeline of the system were ignored; (2) The average irradiation intensity of solar energy and the average atmospheric temperature were specified for any hour; (3) The space-heating load on any day of each natural month was assumed to be constant and stable; (4) The temperature and mass ratio of adsorbed water on the zeolite were uniform at any point in the reaction unit; Energies 2022, 15, 7733 4 of 21 (5) The zeolite particles were packed uniformly in the reactor, and the void fraction was uniform; (6) The thermodynamic properties of zeolites were stable.

Collection of Solar Energy
Equation (1) [17,20] describes the efficiency of solar collection. Meteorological parameters and hourly daily irradiation intensity are shown in Figures 2-4 for the Chifeng City of China. The irradiation intensity increased and then decreased on any day. Furthermore, Figures 5 and 6 show that solar collection efficiency varies with different irradiation intensities and heat-collection temperatures. The efficiency increased with an increase in the irradiation intensity when the heat-collection temperature was constant. The efficiency decreased with increased heat-collection temperature when the irradiation intensity was constant. The heat-collection temperature greatly affected the efficiency under a lower irradiation intensity. To improve the collection efficiency of solar energy, the system collected high-and low-temperature heat at high and low irradiation intensities, respectively. The energy of the solar collection can be calculated using Equation (2). The average collection efficiency of solar energy can be calculated using Equation (3).

Adsorption and Desorption
The rates of desorption and adsorption were calculated using the desorption and adsorption dynamic models expressed in Equations (4) and (5) [28,29]. The equilibrium mass ratio of adsorbed water on the zeolite was obtained during the charging process using Equation (6). During the discharging process, the zeolite was in contact with water, and the equilibrium mass ratio of adsorbed water on the zeolite was set as the maximum mass ratio of 0.261 [28].
(4) Figure 6. Variation of solar-energy collection efficiency in the heating season.

Adsorption and Desorption
The rates of desorption and adsorption were calculated using the desorption and adsorption dynamic models expressed in Equations (4) and (5) [28,29]. The equilibrium mass ratio of adsorbed water on the zeolite was obtained during the charging process using Equation (6). During the discharging process, the zeolite was in contact with water, and the equilibrium mass ratio of adsorbed water on the zeolite was set as the maximum mass ratio of 0.261 [28].
As shown in Figure 7, the transient physisorption capacity of the zeolite in the desorption process obtained using the mathematical model and experiment [29] was in good agreement. Thus, it can be seen that the proposed model is adaptable. As shown in Figure 8, the transient physisorption capacity of the zeolite during the adsorption process obtained using the mathematical model and experiment [30] were in good agreement. Thus, it can be validated that the proposed model is adaptable. The energy balance for the charging process of the system can be expressed using Equation (7), where the term on the left-hand side represents the heat transferred from the heat source in the embedded heat exchanger. Furthermore, the first and second terms on the right-hand side of the equation represent the sensible heat change of the wet zeolite and desorption heat stored by the zeolite, respectively. The heat of charging was obtained As shown in Figure 8, the transient physisorption capacity of the zeolite during the adsorption process obtained using the mathematical model and experiment [30] were in good agreement. Thus, it can be validated that the proposed model is adaptable. As shown in Figure 8, the transient physisorption capacity of the zeolite during the adsorption process obtained using the mathematical model and experiment [30] were in good agreement. Thus, it can be validated that the proposed model is adaptable. The energy balance for the charging process of the system can be expressed using Equation (7), where the term on the left-hand side represents the heat transferred from the heat source in the embedded heat exchanger. Furthermore, the first and second terms on the right-hand side of the equation represent the sensible heat change of the wet zeolite and desorption heat stored by the zeolite, respectively. The heat of charging was obtained  The energy balance for the charging process of the system can be expressed using Equation (7), where the term on the left-hand side represents the heat transferred from the heat source in the embedded heat exchanger. Furthermore, the first and second terms on the right-hand side of the equation represent the sensible heat change of the wet zeolite and desorption heat stored by the zeolite, respectively. The heat of charging was obtained using Equation (8).
The energy balance for the discharge process of the system can be expressed using Equation (9), where the term on the left-hand side represents the generated adsorption heat, and the first and second terms on the right-hand side represent the sensible heat changes of the wet zeolite and supplied water in the embedded heat exchanger, respectively. The heat of discharging is calculated using Equation (10).

Performance Evaluation Indexes
Several indexes were adopted in this study to evaluate the performances of the SSTES system based on zeolite water.
(1) Efficiency of the total system The efficiency of the total system is defined as the amount of obtained space heating divided by the amount of solar heat collection and can be expressed using Equation (11).
(2) Utilization efficiency of solar energy Solar energy utilization efficiency is defined as the amount of obtained space heating divided by the amount of solar radiation heat and can be expressed using Equation (12).
(3) Efficiency of TES based on zeolite water. The efficiency of TES based on zeolite water is defined as the heat of discharging divided by the heat of charging and can be expressed using Equation (13).
(4) Energy density: The mass-energy density of zeolite is defined as the heat of discharging divided by the mass of zeolite and can be expressed using Equation (14).
The mass-energy density of zeolite: The volume energy density of zeolite is defined as the heat of discharging divided by the volume of zeolite and can be expressed using Equation (15). The volume energy density of zeolite: The mass-energy density of the system is defined as the amount of obtained space heating divided by the mass of zeolite and can be expressed using Equation (16).
The mass-energy density of the system: The volume energy density of the system is defined as the amount of obtained space heating divided by the volume of zeolite and can be expressed using Equation (17).
The volume energy density of the system: (5) Entransy of supply water The Entransy of the supply water is defined by Equation (18). Space heating can be expressed using Equation (19).
Entransy of supply water: Space heating:

Case Study
The physical properties and operational conditions adopted in the calculations were used to simulate the characteristics of the proposed system and are listed in Table 1. The non-heating season was defined as the period from 16 April to 14 October, and the remaining months were defined as the heating season.

Performance of the Proposed System during the Non-Heating Season
A typical day in each month was selected to show the characteristics of the charging process in detail. Figure 9 shows a schematic of the temperature of the solar collection on a typical day in July. The system collected high-(135 • C) and medium-temperature (105 • C) heat under high and low irradiation intensities, respectively. The durations for solar collection at 135 • C and 105 • C were planned based on half of the heat of the collection because this value best reflects the characteristics of the proposed collection mode. a typical day in July. The system collected high-(135 °C) and medium-temperature (105 °C) heat under high and low irradiation intensities, respectively. The durations for solar collection at 135 °C and 105 °C were planned based on half of the heat of the collection because this value best reflects the characteristics of the proposed collection mode. The simulation results are listed in Table 2 under the specified conditions on a typical day of each month. The accumulation of heat energy for the entire non-heating season can be obtained by summing the daily amount. Table 3 lists the simulation results of the proposed system during the non-heating season under the specified conditions. As observed from the table, the average collection efficiency of solar energy is 0.69. The ratio of desorption energy to The simulation results are listed in Table 2 under the specified conditions on a typical day of each month. The accumulation of heat energy for the entire non-heating season can be obtained by summing the daily amount. Table 3 lists the simulation results of the proposed system during the non-heating season under the specified conditions. As observed from the table, the average collection efficiency of solar energy is 0.69. The ratio of desorption energy to collection heat is 0.90. The ratio represents the energy-conversion efficiency of the energy-storage process.

Performance of the Proposed System during the Heating Season
A typical day in each month was selected to reveal the characteristics of the discharging process in detail. Here, Figure 10 shows a schematic of the temperature of the solar collection on a typical day in January. The system collected medium-(100 • C) and lowtemperature (70 • C) heat under high and low irradiation intensities, respectively. The duration of solar collection at 100 • C and 70 • C was planned based on half of the heat of the collection.
Energies 2022, 15, x FOR PEER REVIEW 12 of collection heat is 0.90. The ratio represents the energy-conversion efficiency of the energ storage process.

Performance of the Proposed System during the Heating Season
A typical day in each month was selected to reveal the characteristics of the dischar ing process in detail. Here, Figure 10 shows a schematic of the temperature of the sol collection on a typical day in January. The system collected medium-(100 °C) and low temperature (70 °C) heat under high and low irradiation intensities, respectively. The d ration of solar collection at 100 °C and 70 °C was planned based on half of the heat of t collection. The simulation results are listed in Table 4 under the specified conditions on a typic day of each month. The simulation results are listed in Table 4 under the specified conditions on a typical day of each month. The accumulation of heat energy for the entire heating season can be obtained by summing the daily amounts. Here, Table 5 lists the simulation results of the proposed system during the heating season under the specified conditions. As observed from the table, the average collection efficiency of solar energy is 0.68. The proportion of heat supplied at 120 • C reached 0.95. This is because the adsorption heat was insufficient to heat the supply water at 85 • C.

Annual Performance of the Proposed System
The proposed system exhibits a good performance because it adopts solar cascade heat collection, cascade storage of heat energy, and cascade heating. To better illustrate the advantages of the proposed system, we simulated a conventional system and compared the performance of the two systems.
Here, Figure 11 shows a conventional solar-heating system with STES based on zeolite water. In the non-heating season, the solar heat-collection and zeolite adsorption subsystems worked, while the TES, energy release, and heating subsystems did not work. In this operating mode, valves V3-1, V3-2, and V0 were opened. Collected solar energy was eventually stored in the zeolite adsorption reactor. In the heating season, the system adopted heat storage and heat release. During the charging process, high-temperature heat was stored in tank W1, and valves V1-1 and V1-2 were opened in this operating mode. During the discharging process, the backwater was heated using the stored heat, and then, it was used for heating when the heating load was low or when the solar energy was abundant. Valves V2-1, V2-2, V4-1, and V4-2 were opened in this operating mode. The backwater was heated using the stored heat and zeolite adsorption heat simultaneously, and then it was used for heating when the heating load was high, or the solar energy was insufficient. An exothermic adsorption reaction occurred in the reactor, and the dry zeolite became wet. Valves V2-1, V2-2, V4-1, V4-2, V5-1, and V5-2 were opened in this operating mode. Energies 2022, 15, x FOR PEER REVIEW dry zeolite became wet. Valves V2-1, V2-2, V4-1, V4-2, V5-1, and V5-2 were opene operating mode. Figure 11. The conventional solar-heating system with STES based on zeolite water.
The operating conditions of the proposed and reference systems are listed in Table 6. Operating conditions of the proposed and reference systems.

Item
Reference System The performances of the two systems under the specified conditions are liste ble 7. The proposed system had several advantages. It had a supply heat of 2135 G is 188 GJ higher than that of the reference system. The efficiency of the total syst 54.18%, which is 1.28% higher than that of the reference system, and the utilizat ciency of solar energy was 37.16%, which is 3.26% higher than that of the reference  The operating conditions of the proposed and reference systems are listed in Table 6. The performances of the two systems under the specified conditions are listed in Table 7. The proposed system had several advantages. It had a supply heat of 2135 GJ, which is 188 GJ higher than that of the reference system. The efficiency of the total system was 54.18%, which is 1.28% higher than that of the reference system, and the utilization efficiency of solar energy was 37.16%, which is 3.26% higher than that of the reference system.

Sensitivity Analyses
Typical days in July and January were selected, and sensitivity analyses were performed to investigate the characteristics of the proposed system in detail.

Medium Temperature of Solar Collection
Here, Figure 12 shows the variations in the final mass ratio of adsorbed water on zeolite, the mass of zeolite, the heat of solar collection and the desorption heat of zeolite with the medium temperature of solar collection during the non-heating season. As observed from the figure, the mass of zeolite decreased with the increase in the medium temperature of solar collection in the non-heating season. This is because the desorption depth of zeolite at the medium temperature of solar collection increased. The higher the temperature of zeolite, the lower the equilibrium mass ratio of adsorbed water on zeolite was. As the desorption depth of zeolite increased, the mass of zeolite decreased. Similarly, the final mass ratio of adsorbed water on zeolite decreased.

Sensitivity Analyses
Typical days in July and January were selected, and sensitivity analyses were performed to investigate the characteristics of the proposed system in detail.

Medium Temperature of Solar Collection
Here, Figure 12 shows the variations in the final mass ratio of adsorbed water on zeolite, the mass of zeolite, the heat of solar collection and the desorption heat of zeolite with the medium temperature of solar collection during the non-heating season. As observed from the figure, the mass of zeolite decreased with the increase in the medium temperature of solar collection in the non-heating season. This is because the desorption depth of zeolite at the medium temperature of solar collection increased. The higher the temperature of zeolite, the lower the equilibrium mass ratio of adsorbed water on zeolite was. As the desorption depth of zeolite increased, the mass of zeolite decreased. Similarly, the final mass ratio of adsorbed water on zeolite decreased.  As observed from the figure, the heat of solar collection decreased with the increase in the medium temperature of solar collection in the non-heating season. This is because the higher the temperature of solar collection, the lower the efficiency of solar collection and the heat of solar collection. Similarly, the desorption heat of zeolite exhibited the same trend. This is because only the sensible loss was included in thermal loss, and the desorption heat of zeolite was equal to the difference between the heat of solar collection and the thermal loss. Figure 13 shows the variations in the final mass ratio of adsorbed water on zeolite, the mass of zeolite, the heat of solar collection and the desorption heat of zeolite with the high temperature of solar collection during the non-heating season. As observed from the figure, the mass of zeolite decreased with the increase in the high temperature of solar collection in the non-heating season. This is because the heat of solar collection at medium temperature decreased. The higher the temperature of solar collection, the lower the heat of solar collection. At half of the total heat of solar collection, the heat of solar collection at medium temperature decreased. Meanwhile, the mass ratio of the adsorbed water on zeolite, which could be desorbed, was unchanged because the equilibrium mass ratio of the adsorbed water on the zeolite remained constant at medium temperature. However, the final mass ratio of the adsorbed water on the zeolite decreased and then remained unchanged. When the high temperature of solar collection in the non-heating season was less than 135 • C, the mass ratio of the adsorbed water on the zeolite that can be desorbed decreased. This is because the equilibrium mass ratio of the adsorbed water on the zeolite decreased. When the high temperature of solar collection in the non-heating season was greater than 135 • C, the mass ratio of the adsorbed water on the zeolite that can be desorbed remained unchanged. This is because the heat of solar collection at high temperature was limited and was equal to half of the total heat of solar collection. As observed from the figure, the heat of solar collection decreased with the increase in the medium temperature of solar collection in the non-heating season. This is because the higher the temperature of solar collection, the lower the efficiency of solar collection and the heat of solar collection. Similarly, the desorption heat of zeolite exhibited the same trend. This is because only the sensible loss was included in thermal loss, and the desorption heat of zeolite was equal to the difference between the heat of solar collection and the thermal loss. Figure 13 shows the variations in the final mass ratio of adsorbed water on zeolite, the mass of zeolite, the heat of solar collection and the desorption heat of zeolite with the high temperature of solar collection during the non-heating season. As observed from the figure, the mass of zeolite decreased with the increase in the high temperature of solar collection in the non-heating season. This is because the heat of solar collection at medium temperature decreased. The higher the temperature of solar collection, the lower the heat of solar collection. At half of the total heat of solar collection, the heat of solar collection at medium temperature decreased. Meanwhile, the mass ratio of the adsorbed water on zeolite, which could be desorbed, was unchanged because the equilibrium mass ratio of the adsorbed water on the zeolite remained constant at medium temperature. However, the final mass ratio of the adsorbed water on the zeolite decreased and then remained unchanged. When the high temperature of solar collection in the non-heating season was less than 135 °C, the mass ratio of the adsorbed water on the zeolite that can be desorbed decreased. This is because the equilibrium mass ratio of the adsorbed water on the zeolite decreased. When the high temperature of solar collection in the non-heating season was greater than 135 °C, the mass ratio of the adsorbed water on the zeolite that can be desorbed remained unchanged. This is because the heat of solar collection at high temperature was limited and was equal to half of the total heat of solar collection.  As observed from the figure, the heat of solar collection decreased with the increase in the high temperature of solar collection in the non-heating season. This is because the higher the temperature of solar collection, the lower the efficiency of the solar collection and the heat of the solar collection. However, the desorption heat of zeolite first increased rapidly and then decreased gradually. This is because the desorption heat of zeolite was equal to the difference between the heat of solar collection and the thermal loss. When the high temperature of solar collection in the non-heating season was less than 135 • C, the thermal loss significantly decreased because the mass ratio of the adsorbed water on the zeolite that can be desorbed decreased. In addition to the sensible heat loss, much of the heat of solar collection could not be converted into desorption heat of zeolite because of the limitation of the equilibrium mass ratio of the adsorbed water on the zeolite. When the high temperature of solar collection in the non-heating season was greater than 135 • C, only sensible loss was included in the thermal loss.  Figure 14 shows the variations in the temperature, the mass, the heat and the entransy of supply water with the low temperature of solar collection during the heating season. As observed from the figure, the mass of supply water decreased with the increase in the low temperature of solar collection in the heating season. This is because the heat of solar collection at a low temperature decreased with the increase in the low temperature of solar collection in the heating season. However, the temperature of the supply water increased and then remained unchanged. When the low temperature of solar collection in the heating season was less than 70 • C, the temperature of the supply water increased because of the decrease in the mass of the supply water. When the low temperature of solar collection in the heating season was greater than 70 • C, the supply water temperature remained unchanged. This is because the medium temperature of solar collection was constant during the heating season. As observed from the figure, the heat of solar collection decreased with the increase in the high temperature of solar collection in the non-heating season. This is because the higher the temperature of solar collection, the lower the efficiency of the solar collection and the heat of the solar collection. However, the desorption heat of zeolite first increased rapidly and then decreased gradually. This is because the desorption heat of zeolite was equal to the difference between the heat of solar collection and the thermal loss. When the high temperature of solar collection in the non-heating season was less than 135 °C, the thermal loss significantly decreased because the mass ratio of the adsorbed water on the zeolite that can be desorbed decreased. In addition to the sensible heat loss, much of the heat of solar collection could not be converted into desorption heat of zeolite because of the limitation of the equilibrium mass ratio of the adsorbed water on the zeolite. When the high temperature of solar collection in the non-heating season was greater than 135 °C, only sensible loss was included in the thermal loss. Figure 14 shows the variations in the temperature, the mass, the heat and the entransy of supply water with the low temperature of solar collection during the heating season. As observed from the figure, the mass of supply water decreased with the increase in the low temperature of solar collection in the heating season. This is because the heat of solar collection at a low temperature decreased with the increase in the low temperature of solar collection in the heating season. However, the temperature of the supply water increased and then remained unchanged. When the low temperature of solar collection in the heating season was less than 70 °C, the temperature of the supply water increased because of the decrease in the mass of the supply water. When the low temperature of solar collection in the heating season was greater than 70 °C, the supply water temperature remained unchanged. This is because the medium temperature of solar collection was constant during the heating season.  As observed from the figure, the heat of the supply water first decreased gradually. It then reduced rapidly with the increase in the low temperature of solar collection in the heating season. When the low temperature of solar collection in the heating season was less than 70 • C, the heat of the supply water decreased gradually because the heat of solar collection decreased. When the low temperature of solar collection in the heating season was greater than 70 • C, the heat of supply water decreased rapidly. This is because much of the heat of solar collection at medium temperature could not be converted into the heat of the supply water owing to the limitation of the temperature difference in heat exchange. However, the entransy of the supply water first increased gradually and then decreased rapidly. The maximum entransy of supply water was 875.1 MJK/(m 2 *d) when the low temperature of solar collection in the heating season was 70 • C.

Medium Temperature of Solar Collection
Here, Figure 15 shows the variations in the temperature, the mass, the heat and the entransy of supply water with the medium temperature of solar collection during the heating season. As observed from the figure, the mass of the supply water decreased with the increase in the medium temperature of solar collection in the heating season. This is because the heat of solar collection at low temperatures decreased. At half of the total heat of solar collection, the heat of solar collection at low temperatures decreased. However, the supply water temperature first increased and then remained unchanged. When the medium temperature of solar collection in the heating season was less than 100 • C, the temperature of the supply water increased. This is because of the limitation of temperature difference of heat exchange. When the medium temperature of solar collection in the heating season was greater than 100 • C, the temperature of the supply water remained unchanged. This is because the heat of solar collection at medium temperature was limited and was equal to half of the total heat of solar collection in the heating season. As observed from the figure, the heat of the supply water first decreased gradually. It then reduced rapidly with the increase in the low temperature of solar collection in the heating season. When the low temperature of solar collection in the heating season was less than 70 °C, the heat of the supply water decreased gradually because the heat of solar collection decreased. When the low temperature of solar collection in the heating season was greater than 70 °C, the heat of supply water decreased rapidly. This is because much of the heat of solar collection at medium temperature could not be converted into the heat of the supply water owing to the limitation of the temperature difference in heat exchange. However, the entransy of the supply water first increased gradually and then decreased rapidly. The maximum entransy of supply water was 875.1 MJK/(m 2 *d) when the low temperature of solar collection in the heating season was 70 °C.

Medium Temperature of Solar Collection
Here, Figure 15 shows the variations in the temperature, the mass, the heat and the entransy of supply water with the medium temperature of solar collection during the heating season. As observed from the figure, the mass of the supply water decreased with the increase in the medium temperature of solar collection in the heating season. This is because the heat of solar collection at low temperatures decreased. At half of the total heat of solar collection, the heat of solar collection at low temperatures decreased. However, the supply water temperature first increased and then remained unchanged. When the medium temperature of solar collection in the heating season was less than 100 °C, the temperature of the supply water increased. This is because of the limitation of temperature difference of heat exchange. When the medium temperature of solar collection in the heating season was greater than 100 °C, the temperature of the supply water remained unchanged. This is because the heat of solar collection at medium temperature was limited and was equal to half of the total heat of solar collection in the heating season.  As observed from the figure, the heat of supply water first increased rapidly and then decreased gradually with the increase in the medium temperature of solar collection in the heating season. When the medium temperature of solar collection in the heating season was less than 100 • C, the heat of supply water increased rapidly. This is because an increasingly lesser amount of heat from solar collection at medium temperature could not be converted into the heat of the supply water owing to the limitation of the temperature difference in the heat exchange rate. When the medium temperature of solar collection in the heating season was greater than 100 • C, the heat of the supply water decreased gradually because the heat of solar collection decreased. The maximum heat of the supply water was 4.889 MJ/(m 2 *d) when the medium temperature of solar collection in the warm season was 100 • C. Meanwhile, the entransy of the supply water first increased rapidly and then decreased gradually. The maximum entransy of the supply water was 875.1 MJK/(m 2 *d) when the medium temperature of solar collection in the warm season was 100 • C.

Analyses of Annual Performance
Here, Figure 16 shows the variations in the efficiency of the total system and the utilization efficiency of solar energy full year. As observed from Figure 16a, the efficiency of the total system increased, and the utilization efficiency of solar energy decreased with the increase in the medium temperature of solar collection in the non-heating season. As observed from Figure 16b, the utilization efficiency of solar energy increased rapidly. Subsequently, they decreased gradually, and the efficiency of the total system first increased rapidly and then increased gradually with the high temperature of solar collection during the non-heating season. As observed from Figure 16c, the utilization efficiency of solar energy and the efficiency of the total system first decreased gradually and then decreased rapidly with the increase in the low temperature of solar collection in the heating season. As observed from Figure 16d, the utilization efficiency of solar energy and the efficiency of the total system first increased rapidly. It then decreased gradually with the increase in the medium temperature of solar collection in the heating season.
As observed from the figure, the heat of supply water first increased rapidly and then decreased gradually with the increase in the medium temperature of solar collection in the heating season. When the medium temperature of solar collection in the heating season was less than 100 °C, the heat of supply water increased rapidly. This is because an increasingly lesser amount of heat from solar collection at medium temperature could not be converted into the heat of the supply water owing to the limitation of the temperature difference in the heat exchange rate. When the medium temperature of solar collection in the heating season was greater than 100 °C, the heat of the supply water decreased gradually because the heat of solar collection decreased. The maximum heat of the supply water was 4.889 MJ/(m 2 *d) when the medium temperature of solar collection in the warm season was 100 °C. Meanwhile, the entransy of the supply water first increased rapidly and then decreased gradually. The maximum entransy of the supply water was 875.1 MJK/(m 2 *d) when the medium temperature of solar collection in the warm season was 100 °C.

Analyses of Annual Performance
Here, Figure 16 shows the variations in the efficiency of the total system and the utilization efficiency of solar energy full year. As observed from Figure 16a, the efficiency of the total system increased, and the utilization efficiency of solar energy decreased with the increase in the medium temperature of solar collection in the non-heating season. As observed from Figure 16b, the utilization efficiency of solar energy increased rapidly. Subsequently, they decreased gradually, and the efficiency of the total system first increased rapidly and then increased gradually with the high temperature of solar collection during the non-heating season. As observed from Figure 16c, the utilization efficiency of solar energy and the efficiency of the total system first decreased gradually and then decreased rapidly with the increase in the low temperature of solar collection in the heating season. As observed from Figure 16d, the utilization efficiency of solar energy and the efficiency of the total system first increased rapidly. It then decreased gradually with the increase in the medium temperature of solar collection in the heating season.

Conclusions
Based on the cascade utilization of energy, a novel solar-heating system with seasonal and cascade TES based on zeolite water is proposed. Furthermore, a mode is proposed in which high-temperature heat is collected under a higher irradiation intensity and lowtemperature heat is collected under a lower irradiation intensity. The proposed system is helpful in improving the efficiency of solar collection. In the non-heating season, solar energy was collected and stored in the zeolite. In the heating season, space heating was supplied through cascade heat release. The following conclusions can be drawn from the numerical simulations performed: (1) Compared with the reference system, the solar cascade heat collection could improve the average efficiency of heat collection in the proposed system. The average heat collection efficiency is 69.08%, which is 2.88% higher than that of the reference system during the non-heating season, and 67.71%, which is 7.4% higher than that of the reference system during the heating season. Therefore, the proposed system has a high solar-energy utilization efficiency of 37.16%, which is 3.26% higher than the reference system's. (2) The supply heat generated with the proposed system during the heating season is 2135 GJ, which is 9.66% higher than the reference system. The system's mass-energy density and volume energy density are 586 kJ/(kg zeolite) and 614 MJ/(m 3 zeolite), which are 5.02% and 5.14% higher than that of the reference system, respectively. The system performance depends on the physical properties of the zeolite and the proposed system. (3) According to the sensitivity analyses of the proposed system, the utilization efficiency of solar energy for a full year is generated when the high temperature of the solar collection is 135 • C during the non-heating season under the specified conditions. Similarly, the utilization efficiency of solar energy full year is generated when the medium temperature of solar collection is 100 • C during the heating season when the other parameters remained unchanged.