1. Introduction
There is no doubt that China is the country with the highest energy consumption in the world. Unlike other countries with slow energy growth, China’s energy growth accelerated in 2019 [
1]. Within this growth, the proportion of coal consumption showed a downward trend, while the total consumption of renewable energy continued to increase [
2].
However, the utilization of renewable energy was restricted by environmental and technological factors. Therefore, a complementary system of non-renewable energy and renewable energy came into being [
3]. Additionally, solar energy is a clean and inexhaustible renewable energy, while natural gas is a non-renewable energy with the characteristics of high calorific value, convenient transportation, stable and continuous [
4,
5]. In the early and mid 21st century, the combination of solar energy and natural gas has been listed as the development goal of solar thermal utilization [
6].
The concentrating solar photovoltaic/thermal (CPV/T) system uses water as a circulating cooling medium. The cooling water reduces the working temperature of photovoltaic (PV) cells and improves photovoltaic conversion efficiency [
7]. The CPV/T system significantly improves the comprehensive utilization rate of solar energy and realizes high efficient power generation [
8]. Due to the concentrator, the efficiency of the CPV/T system is higher than that of the conventional PV/T system. In particular, the thermal and electrical efficiency of the LCPV/T system has been shown to reach 54.48% and 14.49%, respectively [
9]. Haiping Chen et al. designed a distributed compound parabolic concentrator (WD-CPC). The result showed that the light uniformity value was 0.153, which decreased the light non-uniformity about 0.13 or less compared with the flat-plate compound parabolic concentrator [
10]. T. T. Chow focused on the development trend of the PV/T hybrid solar technology. In this framework, a review about the type and potential of PV/T technology was reported [
11]. Moreover, Heng Zhang et al. proposed a LCPV/T system which included a low-concentrating compound parabolic concentrator. A baffle heat exchanger channel was applied to decrease the temperature of the PV cells. The thermal efficiency reached 55.11%, while the electrical efficiency reached 12.50% [
12]. Guiqiang Li et al. established a model of PV/T system with static miniature solar concentrator. The simulated result and the experimental result were consistent. Compared with the flat plate PV/T system, the heat loss coefficient of the proposed system was lower [
13]. Zexin Wang et al. designed a hybrid CPV/T system based on an unsteady-state thermal model. The result of the model was almost the same as the experimental output of the system. The thermal efficiency of the unsteady-state thermal model was 55.3%, and the measurement was 55.8% [
14]. Bennett K. Widyolar et al. proposed a hybrid CPV/T system applied gallium arsenide solar cells. The simulated and experimental result showed that the exergy efficiency of the PV/T system could reach 37% when the thermal absorber ran at 500 °C [
15].
Ramos A et al. [
16] proposed a hybrid PV/T system for which the overall efficiency was above 70%, the electrical efficiency was 15–20%, and the thermal efficiency could exceed 50%. After reasonable design of combining with heat pump or chiller, the hybrid PV/T system could satisfy more than 60% of the heating demand and 50% of the cooling demand of urban households. In addition, Ramos A simulated the specific requirements of Seville, Rome, Madrid and Bucharest in a TRNSYS environment. The results showed that the PV/T system satisfied 60% of the heating demand and almost all of the cooling demand. The cost of this system achieved 30–40% lower than the equivalent pure photovoltaic system. Moreover, Braun R et al. [
17] applied PV/T collectors as a heat source and heat sink for a reversible heat pump. Compared with traditional solar cooling system which used a reversible air-water heat pump as a cold or heat source, the overall efficiency and economic efficiency of this system were improved due to the reduction of electricity consumption for heat rejection. The auxiliary energy for this system was electricity from the grid. Bosheng Su et al. [
18] proposed a novel combined cooling, heating and power system that could realize independent control of temperature and humidity. Compared with the traditional vapor compression air conditioning system, this novel system saved 50.41% of electrical energy when the same amount of cooling energy was generated.
In fact, only the CPV/T system cannot satisfy the demand of the user throughout a year. The main reason for this is the random fluctuation and intermittency of solar energy. To make up for the shortage, the CPV/T system needs to be combined with other stable energy supply systems [
19]. As the most important fossil fuel, natural gas has the advantages of high calorific value and convenient transportation. Natural gas has been applied in distributed systems complementary to solar energy, such as the single-effect absorption air-conditioning system that used solar energy and natural gas as heat sources to provide heating and cooling to villas [
20], the solar-natural gas household heating system that formed by the solar hot water system and natural gas wall-hung boiler system combined in series [
21], and power generation system of solar energy combined with natural gas [
22]. Francesco Calise et al. [
23] analyzed the performance of a trigeneration system which consisted of CPV/T, biomass heater, single-effect absorption chiller and multi-effect distillation system. The exergy, economic and environmental performance of such system was studied in a TRNSYS environment. The results showed that the exergy destruction and exergy output increased with the increase of the PV/T area. Ighball Baniasad Askari et al. [
24] studied and compared the economic performance of two kinds of solar systems that provided electricity, space cooling and domestic hot water. Recently, Christos Tzivanidis et al. [
25] compared the performance of three different parabolic trough solar collectors (PTC) poly-generation systems. The energy efficiency of these three systems were 78.17%, 43.30% and 37.45%, respectively. The exergy efficiencies were 15.94%, 13.08% and 12.25%, respectively. Moreover, the payback times were 5.62 years, 7.82 years and 8.49 years, respectively. In particular, Zhang Heng et al. [
26] proposed and analyzed a novel LCPV/T system which combined with absorption chiller. The inlet temperature, cooling water flow rate and refrigerant water flow rate influenced the cooling capacity of the chiller. The average COP of the chiller was 0.52. Furthermore, a LCPV/T triple-generation system was proposed by Liu Yang et al. [
27]. This system could produce hot water at 45–90 °C. Meanwhile, the electrical efficiency of the LCPV/T system was approximately 10%. The COP of the single-effect lithium bromide absorption chiller was above 0.5. To summarize, there are few studies on the dynamic characteristics of the trigeneration system that based on a CPV/T system and natural gas.
As discussed above, the trigeneration system has attracted the attention of a lot of scientists. However, research on the multi-source complementary system combining CPV/T technology with natural gas is extremely scarce. Therefore, this paper proposed a NG-LCPV/T-TG system in a TRNSYS environment. The electrical and thermal energy are produced by the LCPV/T system. This system uses different heat pumps as a primary auxiliary device and natural gas heater as a secondary auxiliary device. A novel form of energy cascade utilization is proposed in this paper. In addition, this system achieves a balance between supply and demand according to the requirements of the building. The temperature of the experimental room was controlled at 20 °C in winter and 25 °C in summer. Subsequently, the COP value of each equipment was improved. This research not only proposes a new energy utilization method, but also proves the feasibility of this novel system by one-year data. In particular, this research has a crucial effect on the promoting the comprehensive utilization of traditional energy and renewable energy.
4. Results and Discussion
Since this research is to investigate the full-condition performance of the NG-LCPV/T-TG system, the data of continuous operation in one year is simulated. The experimental data in reference [
26] is used for the verification of the model used in this article. Take the COP value of the lithium bromide absorption chiller, as a parameter for comparison. The average deviation value of the modeling results with the experimental results is 10%, and the maximum value is 18%. These errors are mainly caused due to the start and stop of the equipment in the experiment, which causes the COP value to be lower than the simulated value. However, an acceptable range is observed among the results. Therefore, the models presented in this paper are feasible and reliable to simulate the performance of the trigeneration system.
The location simulated by the proposed system is Beijing, China. Beijing is a city with four distinct seasons throughout the year. Outdoor meteorological parameters mainly focus on solar radiation intensity, wind speed and outdoor ambient temperature. The annual variation of these three parameters are manifested in
Figure 4. This meteorological parameter uses the CN-Beijing-545110.tm2 from Meteonorm in TRNSYS software. Both load of the experimental room and LCPV/T performance are affected by outdoor meteorological parameters.
Moreover, the average beam radiation intensity of Beijing in summer is lower than that in winter. This is because Beijing is frequently rainy in summer. The daul-axis tracking system is applied such that the beam radiation is always normal to the surface of PV cells. Because the low-concentrating technology is adopted, only beam radiation intensity has an impact on the performance of the proposed system. In
Figure 4, the average value of beam radiation intensity per month is calculated according to the sum of hours that has solar radiation each day. The maximum and minimum amount of the average beam radiation intensity are 391.8 W/m
2 in February and 185.7 W/m
2 in July. In addition, the hours for calculating average value of ambient temperature and wind velocity are the sum of hours in a whole day per month. The detailed number of hours used for calculating the average value per month are listed in
Appendix A-
Table A3. The ambient temperature follows a seasonal law such that is high in summer and low in winter. In winter, it is maintained at around 0 °C. The height measurement for wind velocity is 10 m. The wind speed is slightly higher in March and April.
In order to maintain comfortable, the temperature of the experimental room is controlled around 20 °C during winter and around 25 °C during summer. The allowable fluctuation range is 1 °C. A reference room that is not provided space heating and cooling is set for comparison. As demonstrated in
Figure 5, the temperature of the reference room is associated with the outdoor ambient temperature. In contrast, the temperature of the experimental room is comfortable throughout the year due to the effect of the NG-LCPV/T-TG system.
Figure 6 depicts variations of the electrical and thermal energy produced by the LCPV/T system during a day. In this paper, the 240th hour to 264th hour (24 h) is randomly selected as a typical day. The maximum amount of the electrical and thermal energy are 2.4 kW and 11.1 kW, respectively. The cumulative value of the electrical energy is 16.4 kWh, and the cumulative thermal energy is 73.5 kWh. In summary, the LCPV/T system has excellent capacity for producing electrical and thermal energy simultaneously.
As shown in
Figure 7 and
Figure 8, the electrical and thermal energy produced by the LCPV/T system during one year is presented. At the 7883th hour, when the radiation intensity is 776 W/m
2, the maximum amount of electrical energy is 2.53 kW. At the 7885th hour, when the radiation intensity is 805 W/m
2, the maximum amount of thermal energy is 13.85 kW. The electrical energy and thermal energy are affected by the radiation intensity.
The outlet water temperature on a day and a year of the LCPV/T system are obviously obtained in
Figure 9 and
Figure 10, respectively. The supreme outlet temperature in a day is detected at nearly 45 °C. In addition, the outlet temperature gradually rises from 08:00 and increases rapidly at 10:00, then decreases until 16:00. The outlet temperature in one year is around 45 °C, which only reaches 65 °C for a few days in summer. The lowest temperature in winter is relatively low, while that in summer is about 20 °C.
In particular, the temperature of the storage tank is not below zero in the non-operation stage of LCPV/T system (shown in
Figure 11) due to the heat storage property. Because of the thermal inertia, the rising time of water temperature in the storage tank is later than that the rising time of the outlet temperature of the LCPV/T system. Moreover, the temperature of storage tank 1 after 18:00 exceeds that before 08:00. As shown in
Figure 12, the temperature of tank 1 is maintained between 10–45 °C throughout the year. Consequently, storage tank 1 can be treated as a stable cryogenic heat source of the NG-LCPV/T-TG system.
In the heating mode, it is clearly visible that the LHP is on when the temperature of storage tank 1 is between 15 °C and 40 °C. Additionally, when the temperature of storage tank 1 is below 15 °C, the auxiliary natural gas heater is applied. As shown in
Figure 13, the time nodes of heating mode are the 1761th hour and 7328th hour, respectively. Additionally, the COP of the LHP is around 5, which indicates a significant improvement compared with the rated value of 3.8. Since the terminal of heat dissipation adopts capillary network system, the outlet water temperature at the load side of the LHP is maintained at about 35 °C to satisfy the temperature requirements of the experimental room. Obviously, such a combination takes sufficient advantage of the low temperature heat source, and makes a remarkable contribution to saving energy.
Figure 14 clearly shows the proportion of the thermal energy produced by LHP and NGAH in heating mode, respectively. When the thermal energy exceeds the demand of space heating, it will preheat domestic hot water. Obviously, the produced thermal energy is higher than the heat load of the experimental room. Furthermore, the produced thermal energy of NGAH is higher than that of the LHP in January and December. The reason for this is that the outdoor environment temperature is low, and the radiation intensity is weak in January and December. Under the above environment conditions, the thermal energy produced by LCPV/T system is slighter, which leads to the low temperature of storage tank 1. In February and November, the produced thermal energy of the LHP is significantly higher than that of the NGAH, because the outdoor environment temperature and the radiation intensity increase significantly in February and November.
In addition, the monthly data relevant to the accumulated thermal energy under the heating mode is listed in
Table 9. In January, the heat load is the highest, with the accumulated value up to 906 kWh. The thermal energy produced by the LHP and the NGAH are 311 kWh and 664 kWh, respectively. Besides, the thermal energy to preheat domestic hot water is 573 kWh. In particular, the maximum amount of preheating domestic hot water is 928 kWh in February. Furthermore, LCPV/T system provides 131 kWh thermal energy for heating directly. The thermal energy produced by NGAH is only 101 kWh and 17 kWh in February and November.
In cooling mode, the HHP provides 85–90 °C hot water as the heat source of the chiller.
Figure 15 refers to the performance when outlet temperature of HHP maintained at about 90 °C in a stable manner. In addition, the COP of HHP can reach 2.2. Because the outlet temperature on the load side of HHP is high, the COP is lower than that of LHP.
Figure 16 depicts that the single-effect absorption chiller can run continuously with the COP up to 0.7 in a year. The COP for this model is defined as the energy transferred from the chilled water divided by the sum of the electrical energy required by the chiller and the energy provided to the chiller by the inlet hot water. Since the energy consumed by an absorption chiller is mainly the heat energy, energy efficiency ratio (EER) is often used as the economic evaluation index. The EER is the ratio of the cooling capacity obtained by the absorption chiller to the heat consumed. The maximum amount of the EER is 0.69 during a year. The chiller is closed when the experimental room falls below 26 °C. Moreover, in the cooling mode, the capillary network system is adopted as well. The outlet temperature of the chiller is generally maintained between 10–20 °C, which can fulfil the requirement of space cooling.
Figure 17 demonstrates the outlet temperature curve of the cooling tower. The cooling tower is started only when the chiller is running. Obviously, the outlet temperature of the cooling tower is higher when the demand of the cooling load is large. Cooling tower outlet water temperature is up to 30 °C, which can fulfill the demand of the chiller.
Table 10 lists the cooling capacity and thermal energy produced in the cooling mode. In March, the cooling load is relatively low, so the energy demand is small. At the same time, the radiation intensity increases compared with that in winter. Therefore, the thermal energy to preheat domestic hot water is the highest in March, up to 1295 kWh. As the ambient temperature rises gradually, the cooling load increases, then the cooling capacity needed is increased. It is also worth noting that the maximum amount of cooling load is 1655 kWh in July. The cooling capacity provided by the chiller is 2389 kWh. The auxiliary natural gas heater is the main heat source of the chiller which provided 3297 kWh thermal energy. In addition, the minimum preheated thermal energy of domestic hot water is 104 kWh in July.
Figure 18 clearly draws the cooling load and the cooling capacity provided by the chiller under the cooling mode. The cooling capacity of the chiller is greater than the cooling load of the room. Due to the fluctuation of the indoor temperature and thermal loss of the water tanks, the cooling capacity provided by the chiller is greater than the cooling load. It is also apparent from
Figure 18 that most of the thermal energy required during summer is supplied by the NGAC. This is because the temperature of storage tank 1 is set at 40 °C when the system determines the mode of heat source supply. In Beijing, although the ambient temperature is high, the radiation intensity is not high enough in summer. Thus, the thermal energy produced by the LCPV/T system cannot satisfy the temperature requirement of the chiller. At this time, the heat source of the chiller is mostly supported by the natural gas heater, which has more economic and energy saving benefits.
Conspicuously,
Figure 19 demonstrates the power consumption of all the electrical equipment of the proposed system and produced electrical energy of LCPV/T system in one year. In the heating mode, electrical energy produced by the LCPV/T system is sufficient to meet the power consumption of the LHP and circulating water pumps (pump 1 and pump 2) in the proposed system. Nevertheless, the electrical energy produced by the LCPV/T system can only satisfy the power consumption of the chiller and the HHP in the cooling mode. The power consumption of the circulating water pumps (pump 3, pump 4, and pump 5) are supported by external grid power. In addition, since the selection of the three pumps is based on the data of previous experiment, the energy consumption of the pumps is too large. Future research will further optimize the energy consumption of the pumps.
In particular, the specific data of the power consumption and electrical energy produced by the LCPV/T system are listed in
Table 11. According to the data, the maximum amount of the electrical energy produced by LCPV/T is 403.6 kWh in April. In July, the maximum power consumption of the circulating pumps, chiller and HHP are 3370.5 kWh, 245.5 kWh, and 30.3 kWh, respectively. The maximum amount of the power consumption of HHP is 64.9 kWh in May. In February, March and November, the power consumption and the electrical energy produced by the LCPV/T system are balanced. Nevertheless, some power needs to be purchased from the grid in other months. The proposed system can simultaneously provide the space cooling and heating, part of the electrical energy required in the building, throughout the year.