Multi-Objective Optimal Design of Renewable Energy Integrated CCHP System Using PICEA-g

: The integration of renewable energies into combined cooling, heating, and power (CCHP) systems has become increasingly popular in recent years. However, the optimization of renewable energies integrated CCHP (RECCHP) systems (i.e., optimal component conﬁgurations) is far from being well addressed, especially in isolated mode. This study aims to ﬁll this research gap. A multi-objective optimization model characterizing the system reliability, system cost, and environmental sustainability is constructed. In this model, the objectives include minimization of annual total cost (ATC), carbon dioxide emission (CDE), and loss of energy supply probability (LESP). The decision variables representing the conﬁguration of the RECCHP system include the number of photovoltaic (PV) panels and wind turbines (WTs), the tilt angle of PV panels, the height of WTs, the maximum fuel consumption, and the capacity of battery and heat storage tanks (HSTs). The multi-objective model is solved by a multi-objective evolutionary algorithm, namely, the preference-inspired coevolutionary algorithm (PICEA-g), resulting in a set of Pareto optimal (trade-off) solutions. Then, a decision-making process is demonstrated, selecting a preferred solution amongst those trade-off solutions by further considering the decision-maker preferences. Furthermore, on the optimization of the RECCHP system, operational strategies (i.e., following electric load, FEL, and following thermal load, FTL) are considered, respectively. Experimental results show that the FEL and FTL strategies lead to different optimal conﬁgurations. In general, the FTL is recommended in summer and winter, while the FEL is more suitable for spring and autumn. Compared with traditional energy systems, RECCHP has better economic and environmental advantages.


Introduction
The combined cooling, heating, and power (CCHP) system-also known as the tri-generation system-has received increasing attention in both academia and industry in recent years because of its high efficiency in the use of energy sources [1,2]. The CCHP system has shown to be capable of converting about 80% of the available fuel into usable energy, which is much more efficient than separation systems [3,4]. Given the depletion of fossil energy and degradation of environment, renewable energies are currently strongly recommended to be integrated into CCHP systems [5]. The use of renewable energies facilitates the application of CCHP systems in places where the power grid is unavailable (e.g., islands, plateaus, deserts).
The early studies of CCHP system optimization mainly focus on the system cost minimization [6][7][8][9][10]. Recently, researchers have started to consider system reliability and environmental sustainability, in addition to the system cost [4,[11][12][13][14][15]. Since the optimization model is often non-convex and/or non-linear and involves both integer and real decision variables, a number of algorithms have been proposed to solve the model. For example, Lahdelma et al. [16] modelled the hourly CCHP operation is presented; (II) a multi-objective model considering the criteria of system cost, system reliability, and environmental sustainability is built to optimize configurations of the RECCHP system; (III) PICEA-g has been applied for the first time to the design of a multi-objective CCHP system with renewable energies [46]. Lastly, the effect of improved FEL and FTL strategies is discussed for the RECCHP system.
The rest of this paper is organized as follows. Section 2 introduces the structure, mathematical model, evaluation criteria, and operational strategies of the considered RECCHP system. Section 3 elaborates the problem solver (i.e., the multi-objective evolutionary algorithm, PICEA-g). The case study and results are discussed in Section 4. Section 5 concludes this study and identifies some future research directions.

Traditional Energy System
Traditional energy systems in isolated mode provide the electrical demand by the power generated by the generator directly and meet the heat and cold demand using electric boilers and electric chillers, respectively (see Figure 1). In other words, all electric, heat, and cold demands are from the electricity generated by a generator, which means the waste heat generated during generator operation will be wasted.

Structure of RECCHP System
A typical CCHP system mainly contains four parts: energy sources, power generation units (PGUs), energy conversion devices, and user load demands, as shown in Figure 2.
Energy sources include fossil energies(e.g., natural gas, coal); and renewable energies (e.g., solar energy, wind energy). A PGU contains prime movers and heat recovery systems. The prime movers can be boilers, turbines, internal combustion engines, Stirling machines, fuel cells, and so on. Usually, the prime mover works with the heat recovery system, producing electricity and thermal energy. Conversion devices contain heat conversion devices (e.g., heat exchangers, electric boilers, and heating coils), cold conversion devices (e.g., electric chillers, absorption chillers, and adsorption chillers), and storage devices (e.g., batteries and heat storage tanks, HSTs) [25,47]. The electrical, cold, and heat energy are the main forms of user demand. The thermal demand is equal to the sum of cold and heat demand from user.
According to the composition of the CCHP system, Figure 3 illustrates the structure of the RECCHP system in isolated mode. The RECCHP system contains four parts: renewable energies devices, conversion devices, storage devices, and the user load. The renewable energies are converted into usable energy through photovoltaic (PV) panels and wind turbines (WTs). The function of the conversion device is to change the form of energy to meet the user demand. Additionally, given the stochastic feature of renewable energies, energy storage devices are required to improve the energy efficiency and system reliability. The introduction to the working procedure of the RECCHP is in Section 2.4.

Component Model
In order to formulate the performance of the RECCHP system, the mathematical models including the PV panel, the wind turbine, the battery, HST, and PGU are constructed as follows.

PV Panel
PV panels convert solar irradiation into electricity. The output electricity of a PV panel can be calculated as Equation (6) [48]. From the horizontal component of the solar radiation (S), the incident radiation on the inclined surface (S t ) can be calculated as follows [49]: As shown in Figure 4, S p represents solar radiation perpendicular to the inclined PV panel, which is used to calculate the maximum output power of PV panels. The relationship equations are shown as follows [50,51]: where FF and N pv respectively represent the fill factor and the total number of PV panels. I sc (t) and V oc (t) respectively denote the short-

Wind Turbine
Wind power generation is closely related to the wind velocity, and can be calculated as Equation (7) [52]. When the wind speed exceeds V f (m/s), wind turbines stop working to protect themselves. with: where V c and V f are the cut-in and the cut-off wind velocity, respectively; V r is the wind velocity measured at the reference height H r . η wt is the coefficient of wind power generation. The air density is marked as ρ. A wt represents the cross-section of the rotor and P wtr is the rated power of the wind turbine. H tow is the height of the wind tower. The power law coefficient is denoted as γ.

Renewable Energy Generation
As can be seen from Figure 3, the renewable energy generation (E REE ) comes from the solar energy (E pv ) and the wind power generation (E wt ). Therefore, it can be calculated by Equation (9) and Equation (10): where E REE el is renewable energy generation meeting the electrical demand. Q eb and Q ec are the heat and cold energy generated by electric boiler and electric chiller. η eb and η ec represent the efficiency of the electric boiler and electric chiller, respectively. E REE bat is excess renewable energy generation stored in the battery.

Battery
The activity procedure of the battery is described as Equation (11) [50]: S.T. Soc min Soc(t) Soc max ,

Heat Storage Tank
The function of the HST is described as Equations (13) and (14): Constraints in the process are explained as follows: where Q HST in and Q HST out denote the input and output thermal energy of HST interacted with outside, retained within [0, Q HST in,max ] and [0, Q HST out,max ], respectively; η HST in and η HST out are their efficiency. C HST represents the content of HST.

Power Generation Unit
The electricity (E pgu ) and the thermal energy (Q pgu ) generated by the PGU are calculated as follows: where F is the fuel supplied to the PGU, and it produces electricity (E pgu ) and waste heat as a byproduct. The waste heat can be reused as thermal energy (Q pgu ) by a heat recovery system, whose efficiency is η r . As shown in Figure 3, Q rc and Q rh are the recovery heat supplied to the absorption chiller and the heating coil; Q ac and Q hc respectively represent the cooling and heating from the absorption chiller and the heating coil. The efficiency of the absorption chiller and heating coil are marked as COP ac and η hc , respectively.
The electrical and thermal balance of the RECCHP system in isolated mode are defined as Equations (19) and (20): where E el , Q hl , and Q cl represent the electrical, heat, and cold demand, respectively. U BAT ch , U BAT dis , V HST in , and V HST out are binary variables denoting the operation status of the battery and HST.

Operation Strategy
Typically, the energy supply procedure in the RECCHP system is as follows. First, the power generated by renewable energies is used to meet the load demand (to maximize the usage of renewable energies). If the load demand is not satisfied, storage devices start working to supply energy. If the load demand still cannot be satisfied, the PGU starts to work. When the energy of the storage device is used up and the PGU works at rated power, the load demand is still unmet, then some load demands will be cut off and a loss of energy supply is recorded. With respect to the FEL and FTL operation strategies, the detailed RECCHP operational process is described in the following sections.

FEL Operation Strategy
In the FEL operation strategy, the electrical demand should be satisfied with a higher priority than the thermal demand. In this case, no extra electricity will be generated from the RECCHP system, since the amount of electricity from RECCHP follows the electrical demand. Besides, the thermal energy generated from the RECCHP system meets the thermal demand for the heating and cooling process. Note that the extra thermal energy might be generated by the PGU, which would be stored in the HST. Moreover, through the electric boiler and electric chiller, the electricity from the battery can be easily converted to the thermal energy used to meet the thermal demand. The FEL operational strategy is shown in Figure 5.

Case 1.
If DC 1 is TRUE, the electrical load demand is fully met by the renewable energy generation. The discriminant condition (DC) is expressed as Equation (22): Then, the RECCHP system runs as follows: 1. If DC 2 is TRUE, excess renewable energy is stored in the battery until it is fully charged. The discriminant condition is expressed as Equation (23): with: where E l is the amount of renewable energy generation meeting load demands. 2. If DC 2 is FALSE and DC 4 is TRUE, HST starts working to supply the thermal energy. The discriminant condition is expressed as Equation (25). The remaining cold and heat demand are denoted as Q R cl and Q R hl .
3. If DC 4 is FALSE and DC 6 is TRUE, PGU starts working to meet the gap of the thermal demand. Meanwhile, the unused electricity is stored into batteries. The discriminant condition is shown as Equation (26): 4. If DC 6 is FALSE, a loss of the thermal energy is recorded, that is, the count of the thermal energy loss (LossCount q ) plus one, as defined in Equation (27). In this state, HST is exhausted and PGU works at the rated power. The unused electricity is used to charge batteries.
If DC 1 is FALSE and DC 3 is TRUE, the system runs as follows. The discriminant condition is expressed as Equation (28): 1. If DC 4 is TRUE, HST starts to supply the thermal energy. The discriminant condition is expressed as Equation (29).
2. If DC 4 is FALSE and DC 6 is TRUE, PGU starts to work. Meanwhile, the unused electricity is stored in batteries. The discriminant condition is expressed as Equation (30): 3. If DC 6 is FALSE, a loss of the thermal energy is recorded and LossCount q increases one, as defined in Equation (27). In this state, HST is used up and PGU works at the rated power. Meanwhile, the unused electrical energy is stored into batteries.
Case 3. If DC 3 is FALSE and DC 5 is TRUE, the next steps are the same as 1, 2, and 3 in Case2. The discriminant condition is expressed as Equation (31): Case 4. If DC 5 is FALSE, a loss of the electrical energy is recorded; that is, the count of electrical energy loss (LossCount P ) plus one, defined as Equation (32). In this state, batteries are exhausted and PGU works at the rated power. In addition, the remaining thermal energy is stored into HST.

FTL Operation Strategy
In the FTL operation strategy, the thermal demand is satisfied with a higher priority than the electrical demand. In this case, no extra thermal energy will be generated, since the amount of thermal energy from RECCHP follows the thermal demand. Additionally, excess electricity would be stored in batteries if the electrical demand is less than the electricity from RECCHP system. The FTL operational strategy is shown in Figure 6. The operational details of the FTL strategy are analogous to the FEL strategy, and thus are not repeated here. However, the related discriminant condition and flowchart are described in Algorithm 1. The optimization problems of the above RECCHP model are seen as multi-objective problems (MOPs), which is implemented in MATLAB (2015(b), MathWorks, Natick, MA, USA) and three objective functions (i.e., annual total cost (ATC), carbon dioxide emission (CDE), and loss of energy supply probability (LESP)) are optimized simultaneously using a multi-objective evolutionary algorithm. MOPs often appear in real-world design scenarios, requiring simultaneous optimization of two or more goals. As these objectives are generally conflicting with each other, the optimal solution of MOPs is not a single one, but rather a set of Pareto optimal solutions [53][54][55]. Due to the population-based approach, the multi-objective evolutionary algorithms are well-suited to solve the MOPs, since they can naturally lead to an approximately trade-off surface (or Pareto front) in operation.

Algorithm 1: Flowchart of FTL
The multi-objective evolutionary algorithms aim to find solutions satisfying all the conflicting objectives ( f ) simultaneously. Without loss of generality, a general multi-objective evolutionary minimization problem is stated as follows: where Ω is the decision space. A solution x is a vector of n decision variables: x = x 1 , x 2 , · · · , x n . m is the number of objectives; g i (x) and h j (x) are the inequality and equality constraints, and p and q are the number of them, respectively. F: Ω → R m include m real-valued objective functions and R m represents the objective space. In this study, the state-of-the-art MOEA, PICEA-g, is applied as the problem solver. It is generally agreed that a set of preferences could orientate the search toward a subset of Pareto front interested by the decision maker; that is, different preferences sets may result in different Pareto fronts. The PICEA-g is a co-evolutionary algorithm in which the candidate solutions are co-evolved with a set of goal vectors. The formulations for calculating the fitness F s of the candidate solution s and the fitness F g of the preference g is given as follows [46]: where n g represents the number of solutions satisfying the preference g. N is the population size of candidate solutions. G denotes the initial goal vectors set and G c is the goal vectors set after genetic variation. when s does not satisfy any g, the fitness F s is defined as 0.
The candidate solutions try to satisfy (Pareto dominate) as many goal vectors as possible, being guided towards the Pareto optimal front. The goal vectors try to avoid being satisfied by candidate solutions. Thus, the main idea of PICEA-g is that if multiple sets of hypothetical preferences are specified during the co-evolution process simultaneously, a meaningful subset of the Pareto front for a decision maker can be achieved. PICEA-g has been demonstrated as effective both on benchmarks and some real-world applications [50,56,57]. For more information about PICEA-g, please refer to the literature [46,58,59]. The source code of PICEA-g is available online at: http://ruiwangnudt.gotoip3. com/optimization.html.

Optimization Objectives
Three objectives from the aspects of economy, environment, and reliability are considered in the optimization of the RECCHP system. The three objectives are annual total cost, carbon dioxide emission, and loss of energy supply probability.

Annual Total Cost
Annual total cost (ATC) contains three parts; i.e., the initial investment cost (C inv ), the operating and maintenance cost (C om ), and the fuel cost (C F ). They are calculated as follows: • The annual operating and maintenance cost can be written as: • The annual fuel cost is described as: • The initial investment cost can be written as: with: CRF denotes the capital recovery factor, written as follows: with: where C REE cap , C PGU cap , and C S cap respectively represent the capital of the renewable energy device, PGU, and the storage device. Cap bat and Cap hst denote the capital of the battery and HST; Cap wt , Cap tow , Cap ac , Cap hc , and Cap pgu are the capital of the wind turbine, the wind tower, the absorption chiller [60], the heating coil, and PGU. C PGU max represents the maximum capacity of the PGU. The operation and maintenance cost of the PV panel, the wind turbine, the wind tower, PGU, the absorption chiller, and the heating coil are respectively marked as OM pv , OM wt , OM tow , OM pgu , OM ac , and OM hc . H tow is the height of the wind tower. T pv , T wt and T pgu represent the running time of the PV panel, the wind turbine, and the PGU. i real is the annual real interest rate. The annual inflation rate and nominal interest rate, denoted as in f and i n , are 1.5% and 3.75%, respectively. L, L bat , and L hst are the lifetime of the system, the battery, and HST. F price is the fuel price.

Carbon Dioxide Emission
The carbon dioxide emission (CDE) is used to measure the environmental sustainability of the system, which can be calculated using the emission conversion factor η CO 2 .
where F(t) denotes the hourly fuel consumption.

Loss of Energy Supply Probability
Loss of energy supply probability (LESP) is defined as the probability of failing to satisfy the electrical demand and the thermal demand. It is calculated by Equation (50): where E total and Q total are the available electrical and thermal energy generated from the RECCHP system. LossCount p and LossCount q represent the count of electrical and thermal energy loss over a year, respectively.

System Simulation and Optimization Model
The decision variables representing the configuration of the RECCHP system include the number of photovoltaic (PV) panels and wind turbines (WTs), the tilt angle of PV panels, the height of WTs, the maximum fuel consumption, and the capacity of battery and heat storage tank (HST). Decision variables related to renewable energies include the number of PV panels (N pv ) and wind turbines (N wt ), the inclination angle of the tilted PV panel (β), and the height of the wind tower (H tow ), which are obtained from the PV and wind tower model. The maximum capacity of PGU (C PGU max ) impacts the selection of other devices and the capital cost [61], and thus is also considered as a decision variable obtained from the PGU model. In addition, the maximum capacity of the battery (C BAT max ) and HST (C HST max ) are two other decision variables obtained from battery and HST model. In summary, decision variables are written as follows: The simulation flowchart in the optimization process is described in Figure 7. E REE (D, H) represents the total renewable energy generation produced by PV panels and wind turbines at hour H of day D. Through this simulation, individuals can obtain the corresponding values of ATC, CDE, and LESP. In order to minimize the three objectives (i.e., ATC, LESP, CDE) simultaneously, a three-objective optimization model is constructed as follows: The parameter range is determined with the daily experience to improve solving speed.
where x is a vector consisting of decision variables; Ω is the decision space; and the objective set is defined as {F(x)|x ∈ Ω}.

Input Data
This section presents a case study to verify the proposed model and algorithm. A place (latitude 41.65 • ) in Spain is considered. Meteorological data are averaged over ten years that include the solar radiation intensity, the wind velocity at 10 m altitude, and the environmental temperature, as shown in Figure 8.
The hourly electrical, heat, and cold demand in representative days of spring, autumn, summer, and winter are taken from the study [62]. The distributions of electrical, heat, and cold load demands during a day are shown in Figure 9.

Representative day in winter
Electricity Heating Cooling Figure 9. Hourly electrical, cold, and heat demands in representative days.
From Figure 9, the electrical demand is basically stable, while cold and heat demands vary greatly throughout a year. Specifically, the heat demand is large in winter while it is small in summer. The cold demand is the opposite.

Parameters
All parameters related to the RECCHP system are summarized in Tables 1-3 [28,50]. Additionally, the service life of the RECCHP system is assumed as 15 years in this study, and those of the battery and HST are 5 years. The technical characteristics of RECCHP system components as the inputs to the simulation and optimization process are given in Tables 1-3. The specifications, lifetime, and efficiency of all components are presented in Table 1. Table 2 shows the characteristics of wind turbine modules. Additionally, the wind turbine height is limited between 5 m and 30 m, and the rated power of wind turbine (P wtr ) is 10 kW. The investment cost of a wind turbine tower depends on its length. Besides, Table 3 shows the characteristics of PV modules that encompass the short-circuit current, the open-circuit voltage, the maximum operating voltage, and current. Both short-circuit current temperature coefficient and open-circuit voltage temperature coefficient equal 0.95. The operation and maintenance cost (C om ) of battery, heat storage tank (HST), PGU, and heat recovery system are related to the operation time. The fuel is liquefied natural gas (LNG) whose price is 6.5 CNY/kg. In addition, the parameters settings of PICEA-g are shown in Table 4.

Results and Discussion
The planning time span in this article is one year (8760 h). For each system configuration, we simulated the operation of the various components of the RECCHP system every hour of the year until 8760 h. Then, the indicators of system configuration were obtained (i.e., ATC, CDE, and LESP), determining the system configuration as good or not. Through constant iteration, the system configuration was optimized toward the optimal direction of the indicator. By solving the model with PICEA-g, a set of Pareto optimal solutions was obtained as shown in Figures 10 and 11.
In both Figures 10 and 11, sub-figures (b), (c), and (d) are the left view, main view, and top view of sub-figure (a). As can be seen from sub-figure (b), the results from the optimization problem highlight a clear trade-off between indicator ATC and LESP. Reducing the value of LESP will often increase the value of ATC; that is, a high-energy reliable system is more expensive than a low-demand one. From sub-figure (d), indicators LESP and CDE are also contradictory. Reducing the value of LESP will often increase CDE. In other words, a higher energy reliable system will emit more CO 2 . If the desirable level for LESP and CDE are defined, we can find a lower cost for a better solution with the PICEA-g approach.
Among the obtained Pareto optimal solutions, the decision-maker can further introduce their preference so as to select the preferred one for implementation. For example, three solutions from the Pareto set are selected according to the desired LESP less than 0.1%, and the satisfactory configurations of the RECCHP system are shown in Table 5.
• A reliable RECCHP system should have a low LESP value, which is preferably 0. Three solutions with LESP less than 0.1% for the FEL and FTL strategies are shown in Table 5. • Furthermore, if the decision-maker is more concerned about the system cost, the solution with the smallest ATC should be chosen. Likewise, solutions with the minimum CDE are more likely to be selected if the fuel usage and CO 2 emission are considered as more important for the decision-maker. Table 5 presents the final optimal configurations of the RECCHP under the FEL and FTL strategies. It is found that the optimal system configurations for the two operation strategies are greatly different-for example, solution 1 in both FEL and FTL (marked as bold inTable 5). For FEL, 10 solar panels, 59 degree tilt angle for PV panels, 19 wind turbines, 29 m height for wind tower, 650 kW·h for the capacity of PGU, 145 kW·h for the capacity of battery, 2 kW·h for the capacity of heat storage tank, causing 1.55 (10 5 ) CNY for annual total cost, 0% for LESP, and 4.49 (10 5 ) kg, are the best components of RECCHP system in terms of decreasing cost. For FTL, 20 solar panels, 59 degree tilt angle for PV panels, 16 wind turbines, 25 m height for wind tower, 489.7 kW·h for the capacity of PGU, 0.2 kW·h for the capacity of battery, 439.5 kW·h for the capacity of the heat storage tank, causing 1.35 (10 5 ) CNY for annual total cost, 0% for LESP, and 3.83 (10 5 ) kg, are the best components of the RECCHP system in terms of decreasing cost. To further illustrate the difference, configurations with the smallest ATC are used to simulate a one-year energy supply process of the RECCHP system. The results are shown in Figures 12 and 13.       It can be clearly seen from Figures 12 and 13 that in the simulation process, all the constraints of the system are satisfied. It can be observed that when the renewable energy generation is unable to meet the higher load demand, the storage system will first discharge until the maximum depth of discharge is reached, and then the PGU starts working to supply the deficit energy.

•
In both FEL and FTL strategies, the fuel consumption is large in summer and small in other seasons. This might be caused by the reduction of renewable energy generation and the increase of the cold demand in summer, which can be observed from Figures 12a,b and 13a,b.

•
Comparing the fuel consumption of PGU under FEL and FTL strategies, it is obvious that the fuel consumption in FEL is large in summer and winter while small in spring and autumn. In addition, the difference for FTL in different seasons is as evident as that for the FEL.

•
Regarding the FEL strategy, we can observe from Figure 12c,d that the battery has a large capacity while the capacity of HST is small and sometimes even reaches 0.
• Regarding the FTL strategy, we can observe from Figure 13c,d that the capacity of the battery is small while HST has great capacity. Obviously, in spring and autumn, HST is only slightly used. The reason may be that the thermal demand in these seasons is pretty small, which can be directly met by renewable energy generation. Hence, some thermal energy generated from PGU may be wasted.  The thermal demand in summer and winter is high. The FTL aims to meet the thermal demand as much as possible. Excess electricity generated from the PGU is stored into batteries which can then be used as needed, improving the system performance even further. The FEL aims to meet the electrical demand as much as possible. This strategy is shown as suitable for spring and autumn. In these seasons, the thermal demand is relatively low, which then can be easily met by the recovery thermal energy from the PGU. Using the FEL strategy, in the course of its operation for one year, renewable energy generation is 1.67× 10 5 kW· h, which is equivalent to saving 5.5 × 10 4 · kg of LNG. LNG consumption of the RECCHP system is 2.33 × 10 5 kg, while traditional energy systems need to consume 4.63 × 10 5 kg of LNG. It can be seen that using the RECCHP system reduces the LNG consumption and CO 2 emissions by 50%. Besides, the annual total cost (ATC) of the RECCHP system is 1.55 × 10 5 CNY, while it is 3.03 × 10 5 CNY for traditional energy systems; thus, the RECCHP system has a better economic performance. This is mainly due to the savings in fuel costs. With the increase in fuel prices, the economic advantages of the RECCHP system will be more pronounced. When working in the FTL strategy, the total renewable energy generation is 1.51 × 10 5 kW· h, which is equivalent to saving 5.0 × 10 4 kg of LNG. The LNG consumption of the RECCHP system is 2.05 × 10 5 kg, which means that using the RECCHP system reduces the LNG consumption and CO 2 emissions by 55%. Similarly, the RECCHP system has good economic performance under the FTL strategy because of the reduction of fuel usage.

Conclusions
The combined cooling, heating, and power (CCHP) system has been widely used due to its high efficiency in terms of energy usage. This research studied optimal configurations of renewable energies integrated CCHP systems in isolated mode. A multi-objective optimization model aimed at the minimization of system cost, loss of energy supply, and carbon emission is constructed. It is solved by the multi-objective evolutionary algorithm, PICEA-g, obtaining a set of Pareto optimal solutions (optimal system configurations). These solutions represent different trade-offs among the three objectives. Lastly, the most satisfactory solution is finalized by introducing the decision-maker preference. In addition, it is found that the optimal system configurations vary greatly under FEL and FTL strategies. The FEL strategy is more suitable for spring and autumn, while the FTL strategy is more suitable for winter and summer. It can be seen that RECCHP system has better economic and environmental advantages. More specifically, compared to traditional energy systems, using the RECCHP system reduces the LNG consumption and CO 2 emissions by 50% and 55% in FEL and FTL, respectively. At the same time, the ATC of the RECCHP system is lower.
In terms of future studies, first we would like to explore day-ahead optimal scheduling of CCHP systems. Second, an adaptive selection of operation strategies would be quite useful, which deserves further study. Lastly, more effective algorithms (i.e., dealing with mixed categorical and real decision variables) will be investigated.