Integrated Optimization Design of Combined Cooling , Heating , and Power System Coupled with Solar and Biomass Energy

The combined cooling, heating, and power (CCHP) systems coupled with solar energy and biomass energy can meet the needs of island or rural decentralized and small-scale integrated energy use, which have become increasingly popular in recent years. This study presents a renewable energy sources integrated combined cooling, heating, and power (RES-CCHP) system, driven by a biogas fueled internal combustion engine (ICE) and photovoltaic (PV) panels, which is different from the traditional natural gas CCHP system. Owing to the solar energy volatility and the constraint of biomass gas production, the traditional optimization design method is no longer applicable. To improve the energetic, economic and environmental performances of the system, an integrated design method with renewable energy capacity, power equipment capacity and key operating parameters as optimization variables is proposed. In addition, a case study of a small farm in Jinan, China, is conducted to verify the feasibility of the proposed RES–CCHP system structure and the corresponding optimal operation strategy. The results illustrate that the implementation of the optimal design is energy-efficient, economical and environmentally-friendly. The values of primary energy saving ratio, annual total cost saving rate and carbon emission reduction ratio are 20.94%, 11.73% and 40.79%, respectively. Finally, the influence of the volatility of renewable energy sources on the optimization method is analyzed, which shows that the RES–CCHP system and the method proposed are robust.


Introduction
Exploration of renewable energy sources (RES) is an important way to address the energy and environmental crisis, which has become one of the main tasks and strategic subjects for the development of the national economy and for establishing a prosperous society in the 21st century [1].China spans a large geographical area with rural areas associated with prolonged sunshine durations and abundant renewable energy sources, such as solar and biomass resources, which have unique advantages.It is estimated that the total amount of biomass resources is approximately equivalent to 460 billion kg of standard coal in China each year.However, at present, most of the biomass is directly used for combustion, and the conversion efficiency is only 10%-20%, thus resulting in serious energy waste, environmental pollution and ecological damage.Simultaneously, the rural residents in China have a prominent feature of small dispersion, and the centralized supply model is unable to alleviate the difficulties of heating supply and local pollution in the vast rural and remote areas.Therefore, the traditional electricity-based energy supply is no longer advisable.The combined cooling, heating and power (CCHP) system, based on the principle of energy cascade utilization, can simultaneously meet the electricity, cooling, and heating loads, which improves the energy efficiency and reduces air pollutant emissions dramatically [2][3][4][5][6].Simultaneously the CCHP system provides an effective way for solar and biomass energy utilization [7][8][9][10][11], which especially has received widespread attention in Europe [12][13][14].The choices of primary energy sources and system capacity have a significant impact on the overall performance of the CCHP system, which has been extensively studied by researchers [15][16][17][18][19].
Most of the current researches choose a number of algorithms for the traditional CCHP system driven by natural gas to determine the optimal system capacity.Zeng et al. [20] selected the genetic algorithm (GA) to optimize the capacity of a combined cooling, heating and power system with ground source heat pump (CCHP-GSHP), and the integrated performances of the coupling system were improved.Liu et al. [21] presented a matrix modeling approach to optimize the CCHP system with hybrid chillers, and the optimal power generation unit (PGU) capacity was determined based on the use of the enumeration algorithm, then an illustrative case study was conducted to show that the proposed method performed much better compared to traditional strategies.Wei et al. [22] adopted the non-dominated sorting genetic algorithm-II (NSGA-II) to find a set of Pareto optimal configurations of the CCHP system.Wang et al. [23] applied a particle swarm optimization algorithm (PSOA) to obtain the optimal design of the building cooling heating and power (BCHP) system.
Several documents have focused on the CCHP system with renewable energy sources.Yousefi et al. [24] proposed a CCHP system with photovoltaic/thermal (PV/T) panels, which was optimized based on the economy, energy, and emission objective functions.The GA optimization and an analytical hierarchical process were used to find the best size of the system components.Fani et al. [25] studied a solar-assisted combined cooling, heating, and power (SCCHP) system in office buildings in Iran.Noussan et al. [26] found out the optimal configuration of a biomass-fired CHP system with thermal energy storage from the economic and energetic points of view.Gazda [27] proposed a distributed energy system driven by solar energy and biogas, then analyzed the energy and environmental performance of the system.Ameri et al. [28] developed a mixed integer linear programming (MILP) model for determining the optimal capacity of CCHP systems with photovoltaic systems so that the initial and operating costs of the energy supply system were minimized.Ebrahimi et al. [29] designed a hybrid-CCHP system that consisted of a basic CCHP system and a solar collector for a residential building in five different climates.The maximum rectangle method was applied to identify the optimal size of the prime mover.Ju et al. [30] compared the performance of a CCHP and renewable energy-based hybrid energy system driven by distributed energy resources (DERs CCHP) and a CCHP system driven by natural gas (NG CCHP), which indicated that the former had better operational performance.Li et al. [31] applied the preference-inspired coevolutionary algorithm (PICEA-g) to the configuration of the CCHP system with renewable energies.
Owing to the volatility and intermittency of renewable energy sources and the time-varying load, a source-load mismatch problem emerges, whereby the addition of RES makes the CCHP system structure and energy flow more complicated.Most of the optimization design methods only optimize the capacity of the PGU while primary energy sources have not been optimized, which is applicative for the natural gas-driven CCHP system to achieve the optimal targets since natural gas is sufficient and stable.However, biomass energy is completely different from natural gas since the utilization of biomass energy is affected by daily production, which is a constraint for optimal design.Therefore, the current optimization design of CCHP system is rendered not applicable.
This study presents a RES-CCHP system driven by a biogas-fueled internal combustion engine (ICE) and photovoltaic (PV) panels.Separate optimization of PV and PGU, would not elicit the best results.Thus, comprehensive consideration of renewable energy and load characteristics, four-decision variables (the capacity of the PGU, the number of PV panels, the on-off coefficient of PGU, and the ratio of electric cooling to cooling load) were selected to be optimized in order to achieve collaborative optimization.A multi-objective integrated optimization model that maximizes the annual total cost Energies 2019, 12, 687 3 of 21 saving rate (ATCSR), primary energy saving ratio (PESR), and carbon emission reduction ratio (CERR), and considers the off-design performance of the PGU, is proposed.The GA is adopted to determine the capacity and key operation parameters of the RES-CCHP system.In addition, a case study of a small farm in Jinan, China is conducted to ascertain the feasibility of the proposed RES-CCHP system structure and the corresponding optimization design method.
The remainder of this study is organized as follows: Section 2 describes the structure and energy flows of the RES-CCHP system.Section 3 presents the optimization model.Results of the simulation and optimization are given and analyzed in Section 4. The conclusions are summarized in the last section.
Compared with existing approaches, the main contributions of this paper are summarized as follows: (1) Considering the fluctuation of solar energy and the limitation of biomass energy, the PGU will frequently operate in a variable operating condition.This paper analyzed off-design performance of the PGU, and the daily biogas production was added into the optimization model as a constraint condition, so that the obtained equipment capacity can be more guaranteed to operate in an efficient operation area and the renewable energy has higher utilization.(2) The capacity of the PGU, the number of PV panels, the on-off coefficient of PGU, and the ratio of electric cooling to cooling load were collectively optimized to improve the energetic, economic and environmental performances of the system.(3) The influence of the volatility of renewable energy sources on the design results and system performance is analyzed to prove the proposed RES-CCHP system and the design method are robust.

System Design
The RES-CCHP system consists of PGU, PV panels, a heat recovery unit, an absorption chiller (AC), an electric chiller (EC), a heat exchanger, and an auxiliary boiler, which is shown in Figure 1.Correspondingly, the system design is based on the following points: 1.
The electricity subsystem consists of PV, power grid (PG) and PGU by biogas.The PV will generate low power when the illumination is low, and the deficiency is thus supplemented by biomass power generation.Biomass reserves are abundant and biomass power generation is more controllable.This can be used with noncontrollable photovoltaic systems as an energy mix to improve system stability.The electricity generated can meet the electrical demand and drive the EC to meet the cooling demand.The RES-CCHP system is connected to grid so that the system can sell the surplus electricity to the PG or purchase the cheap power.

2.
The selection of the PGU has a decisive influence on the performance of the entire system.Turbine engines (TEs) and internal combustion engines (ICEs) are two types of power equipment normally used in current CCHP systems.In view of their capacities at a level of tens of kilowatts, ICEs are superior to TEs in terms of efficiency and cost, so a small biogas-fueled ICE is selected as a power plant for the RES-CCHP system.

3.
Because of the small amount of recovered waste heat in the case of a small ICE, the single-effect, hot water absorption chiller is adopted.In addition, a water recycling subsystem is employed to recover waste heat from jacket water (JW) and exhaust gas (EXH) so that hot water is produced for heating and driving the absorption chiller, and the excess heat is provided for biogas fermentation to maintain a constant fermentation temperature and ensure efficient biogas production.When the system waste heat is insufficient, shortfall heat can also be supplemented by the biogas-fueled boiler.

4.
As the coefficient of performance (COP) of single-effect hot water absorption chiller is low, an EC with an increased performance coefficient is used to achieve hybrid cooling.

Energy Flow Analysis
The first step of the system design is to analyze its energy flow.In Figure 1, Epv and Epgu respectively represent the electricity generated by the PV and PGU, while Egb and Egs respectively represent the electricity bought from the grid and sold electricity.Eec is the input power of the electric chiller, and Le, Lh, and Lc, represent the user demands for electricity, heating, and cooling, respectively.

Biogas ICE JW Exchanger
Absorption Chiller EXH Exchanger Electric Chiller Biogas Boiler

Island or rural area Loads
Animal waste， straws and forest wastes The electrical energy balance is expressed as: where δ(t) represents the state variable of interaction with the power grid, i.e., whereby the number 1 indicates that the RES-CCHP system purchases electricity from PG in period t, whereas 0 indicates that the system sells electricity to PG in period t.
The fuel consumption of the PGU is expressed as: where ηe,pgu is the electrical efficiency of the PGU, which are both affected by the part-load ratio (PLR).Meanwhile, the primary energy consumption of the PG is: where ηe,pg and ηd,pg are the efficiencies of the electricity generation and distribution of the PG, respectively.Correspondingly, the heating balance of the RES-CCHP system is: Considering that the system is designed according to the island and rural energy supply conditions, the dependence and impact on the PG should be reduced as much as possible, while the RES-CCHP system operates in a grid-connected mode.Therefore, the system follows the electric demand.

Energy Flow Analysis
The first step of the system design is to analyze its energy flow.In Figure 1, E pv and E pgu respectively represent the electricity generated by the PV and PGU, while E gb and E gs respectively represent the electricity bought from the grid and sold electricity.E ec is the input power of the electric chiller, and L e , L h , and L c , represent the user demands for electricity, heating, and cooling, respectively.
The electrical energy balance is expressed as: where δ(t) represents the state variable of interaction with the power grid, i.e., whereby the number 1 indicates that the RES-CCHP system purchases electricity from PG in period t, whereas 0 indicates that the system sells electricity to PG in period t.
The fuel consumption of the PGU is expressed as: where η e,pgu is the electrical efficiency of the PGU, which are both affected by the part-load ratio (PLR).Meanwhile, the primary energy consumption of the PG is: where η e,pg and η d,pg are the efficiencies of the electricity generation and distribution of the PG, respectively.Correspondingly, the heating balance of the RES-CCHP system is: Energies 2019, 12, 687 5 of 21 where Q h,ac is the input power of absorption chiller, Q ex is heat loss of the system, Q fe is the heat required for biogas fermentation, Q re represents the recovered waste heat from the PGU, and Q b represents the heating power from the auxiliary boiler.The recovered waste heat from the PGU in period t can be expressed as: where Q jw and Q exh are the jacket water heat and the exhaust heat, respectively, and η jw,he and η exh,he are the heat exchange efficiencies of the jacket water and exhaust, respectively.The input power of the auxiliary gas boiler at time t can be expressed as: where η b is the thermal efficiency of the biogas boiler.Thus, the biogas consumption of the RES-CCHP system is: The cooling balance of the CCHP system is: where Q c,ac and Q c,ec are the cooling output of the absorption and electric chillers, respectively.The cooling capacity of the absorption chiller can be expressed as follows: where COP ac is the coefficient of performance of the absorption chiller.
The electricity consumption of the electric chiller is given as: where COP ec is the COP of the electric refrigerator.
The electric cooling ratio α can be expressed as:

Off-Design Performance of the PGU
PGU is the core power equipment of the system, whose off-design characteristics are obvious owing to the fluctuation of load and renewable energy sources, and determine the entire system operation.The off-design performance is simulated by the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) data for a naturally aspirated ICE, as shown in Table 1 [32].
In the table, η th,ice and η e,eg denote the thermal efficiency of the ICE and the electrical efficiency of the electrical generator (EG), respectively.Therefore, the electrical efficiency of the PGU can be expressed as: Energies 2019, 12, 687 6 of 21 According to the above data, the electrical efficiency of the PGU at different PLR can be obtained as follows: where r(t) is the part-load ratio of the PGU.Equivalently, f jw , f exh , and f ex , correspond to the proportion coefficients of the total waste heat of the unit from the jacket water, exhaust gas, and heat loss, respectively, which satisfy the following relationship: The heat coefficient of the jacket water can also be obtained based on data fitting, as shown in Equation ( 15): The heat coefficient of exhaust gas is shown in Equation ( 16): ). ( 16)

Objective Function
The multi-objective (MO) function, including the environment, economic, and energy objectives, is used to evaluate the overall performance of the RES-CCHP system, compared to the separated production (SP) system.In the SP system, the PG is employed to supply the electricity, and the electric chiller and gas boiler are used to produce the cooling and heating energy, respectively.

Energy Objective
It is vital for RES-CCHP to improve the ratio of the energy usage which is the basic criterion of the newly developed system.PESR is the ratio of the energy-saving of the RES-CCHP system compared to the SP system, which can be illustrated as: where F SP is the yearly energy consumption of the SP system, and F CCHP is the energy consumption of the RES-CCHP system in a year.
The energy consumption of the RES-CCHP system is composed of the amount of biogas consumption of the PGU and the boiler and the consumption of the grid.Thus, F CCHP can be calculated in accordance to: where F CCHP,grid is the prime energy consumption of the grid, F CCHP,gas is the biogas consumption of the PGU and the boiler, and T and v are the time interval and standard coal conversion factor (8.13 kWh/kg), respectively.In the SP system, the primary energy is mainly consumed by the grid and auxiliary boilers, so FSP can be calculated in accordance to: where G SP,b is the heating power from the auxiliary boiler in the SP system.

Economy Objective
Economy is a prerequisite for the popularization and application of a large area of the system.The full life cycle cost of the system is adopted in the economy objective, which includes the investment cost, the maintenance cost, and the operation cost of the system equipment.According to the life cycle cost analysis (LCCA), it is converted into the annual cost of the system.The annual total cost saving rate (ATCSR) is used as the economy objective and expressed as [20,33]: where C CCHP and C SP are the annual costs of the RES-CCHP system and the annual cost of the SP system, respectively.The C CCHP and C SP can be described as: where C CCHP,EQ , C CCHP,OM and C CCHP,EC are the annual investment cost, maintenance cost and operation expenses of the RES-CCHP system, respectively, while C SP,EQ , C SP,OM and C SP,EC are the annual investment cost, maintenance cost and operation expenses of the SP system, respectively.The values of C CCHP,EC , C CCHP,EQ and C CCHP,OM can be calculated as: where P gb , P gs , are the purchase and sale prices of the grid, respectively, C gas and C CCHP,IN , are the annual fuel and total investment costs, respectively, R is the capital recovery factor, P gas is the price of biomass gas, σ is the system maintenance cost coefficient.The values of C SP,EC , C SP,EQ and C SP,OM can be calculated as: where E SP,gb denotes the purchase power from the grid at time t, and C SP,IN is the total investment cost of the SP system.The capital recovery factor, R, can be defined as: where r is the rate of interest and k is the service life period of each component.In this study, the values of r and k are 0.08 and 15, respectively.

Environment Objective
Environmental issues are hot issues in the context of the social development, which is an important factor in the development of the RES-CCHP system.The carbon emission reduction ratio (CERR) is expressed as [33,34]: where CE CCHP and CE SP are the carbon dioxide emissions of the RES-CCHP and SP systems, respectively.The carbon dioxide emission of the RES-CCHP can be calculated as: where µ grid and µ gas are the emission factors of the grid and biomass gas, respectively.The carbon dioxide emission of the SP system can be calculated as: where µ ng is the emission conversion factor of natural gas.Thus, the optimization objective function of the RES-CCHP system can be defined as: where 0 ≤ ω 1 , ω 2 , ω 3 ≤ 1, and ω 1 + ω 2 + ω = 1.In this case, ω 1 , ω 2 , ω 3 , are the corresponding weight coefficients.Their values respectively represent the importance of the energy, economic, and environmental influences, respectively.To attain the best performance in the case of the RES-CCHP system, the objective is to maximize MO.The weight coefficients represent the weight of each sub-objective in the optimization solution.The value of each weight coefficient has an important influence on the result of the optimal configuration.In this study, the judgment matrix method was used to classify the three sub-objectives according to their importance in the optimization configuration.In order to maximize the advantage of the system in terms of its energy saving potential, the energy indicator is selected as the first-level indicator in the optimization configuration.The economic indicator represents the economic performance of the joint supply system.This indicator has a greater impact on application promotion and serves as the second-level indicator for optimal allocation.The environmental protection index Energies 2019, 12, 687 9 of 21 represents the reduction in the operation of the CCHP system.Compared to the SP system, the CCHP system has a clear advantage in terms of pollutant emissions, and so, it is used as a third-level indicator.The weight coefficient corresponding to each index was determined by constructing a judgment matrix [35,36].

Optimization Variable
The RES-CCHP system contains a variety of equipment types, and the selection of the optimization variables should thus follow the principles of independence, importance, and decisiveness.
The equipment capacity configuration and optimization operating parameters to be optimized in this article include the following sections.

PV Capacity
As an important source of electrical energy in this system, the capacity of the PV component has a great influence on the system performance.Considering that a photovoltaic battery group is composed of a certain number of photovoltaic battery modules in series and parallel, the relationship between its capacity and the capacity of a single module can be expressed as: where N pv is the capacity of PV, M is the number of photovoltaic cell modules, and N s,pv is the capacity of a single photovoltaic cell module.Therefore, the number of photovoltaic cell modules represents the capacity of the PV.For general buildings, photovoltaic cells are usually placed in open space with good lighting conditions, such as the roof of a building.However, since the photovoltaic cell modules need to occupy a large area, the available area must be considered as a constraint in the configuration process for the capacity of the photovoltaic cell.The relationship between the PV capacity and its usable area can be expressed as: where S pv is the area of PV, S s,pv is the area of a single PV, and S max is the available area.

Capacity and Minimum Load Coefficient of PGU
PGU is the most important power equipment in this integrated power supply system, and many parameters of other equipment types, such as the capacities of the auxiliary boiler, absorption chiller, and electric chiller, are determined by the capacity of the PGU.Therefore, the capacity of PGU should be an important variable in the optimal configuration of the system.Meanwhile, in order to increase the operating efficiency of the PGU, the generator set should not be operated at low-loads for too long.Thus, N pgu and θ are selected as the optimization variables.

Electric Cooling Ratio
In general, the COP of an electric chiller is greater than that of an absorption chiller.Therefore, setting up the electric chiller in the system not only provides supplementary cooling when the absorption chiller is insufficient for cooling, but can also improve the operating efficiency of the RES-CCHP system by rationally setting the electric cooling ratio.Thus, α is selected as the optimization variable.
In combination with the above, N pv , N pgu , θ, and α were chosen in this study.These four variables were solved by the genetic algorithm to improve the comprehensive performance of the RES-CCHP system.However, when the values of θ and α are not constant, the operational strategy would become more complex, which is difficult to achieve in the real application.Therefore, the electric cooling ratio and the minimum load factor of PGU are fixed to simplify the operation strategy in this study.

Constraints
Considering the operational status of the RES-CCHP system, the optimization process should both satisfy the equality constraints expressed by Equations ( 1)-( 11) and the following inequality constraints: where N pgu,max and N pv,max are the upper capacity bounds of the PGU and PV, respectively, F pgu,gas and F b,gas are the amount of biomass gas consumed by the PGU and the boiler, respectively, and F gas is the daily production of the biomass gas.Equations ( 33) and ( 34) ensure that N pgu and N pv are within a reasonably feasible range.Under normal circumstances, the configuration capacity of PGU and the electric chiller in the RES-CCHP system should not exceed the maximum electricity and cooling load values.Therefore, the maximum value of the electricity and cooling load is selected as the upper limit of the capacity configuration of PGU and the electric chiller.The capacity of PV is limited by the site area for the layout.The upper limit of the capacity of PV is determined by the area of the layout site.The daily production of the biogas is determined by the biomass energy available to users.From Table 1, it can be seen that when the electrical load rate is less than 20%, the power generation efficiency of PGU is lower than 14.11%.Therefore, in order to ensure the efficient operation of the PGU, the lower limit of the minimum load rate of the PGU is set to 0.2.

Optimization Algorithm
Since the CCHP system is a complex energy system that can operate at various operating conditions, it is difficult to optimize it using simple mathematical methods.Therefore, this study adopted a genetic algorithm to solve the multi-objective optimization model.The operating process of the solution method is shown in Figure 2. First, relevant parameters need to be set in advance, including the GA parameters, cost parameters, building load, and performance parameters of the RES-CCHP and SP systems.The optimization variables were coded in binary form.Then, the fitness of the objective function was calculated.The searching process will not stop until the convergence condition is met.After that, the optimization results will be obtained, including N pv , N pgu , θ, and α.The calculation process was implemented in MATLAB.

Building Load Simulations
Farms are typical end-users of the RES-CCHP system, with stable annual electric, heating, and cooling demands.In order to verify the correctness of the proposed optimal configuration method, this study considered a suburban farm in Jinan, China, as an example for simulating the year-by-year cold, heat, and electrical loads of the farm's building, based on its energy analysis software eQUEST (DOE-2.2,LBNL, Berkeley, CA, USA), as shown in Figures 3-5.The situation of the farm and biogas are shown in Table 2 and Table 3, respectively.This study establishes a photovoltaic cell module model using the TRNSYS software simulation platform.The annual output power of the PV in the climate environment of Jinan, in China is shown in Figure 6.

Building Load Simulations
Farms are typical end-users of the RES-CCHP system, with stable annual electric, heating, and cooling demands.In order to verify the correctness of the proposed optimal configuration method, this study considered a suburban farm in Jinan, China, as an example for simulating the year-by-year cold, heat, and electrical loads of the farm's building, based on its energy analysis software eQUEST (DOE-2.2,LBNL, Berkeley, CA, USA), as shown in Figures 3-5.The situation of the farm and biogas are shown in Tables 2 and 3, respectively.This study establishes a photovoltaic cell module model using the TRNSYS software simulation platform.The annual output power of the PV in the climate environment of Jinan, in China is shown in Figure 6.

Building Load Simulations
Farms are typical end-users of the RES-CCHP system, with stable annual electric, heating, and cooling demands.In order to verify the correctness of the proposed optimal configuration method, this study considered a suburban farm in Jinan, China, as an example for simulating the year-by-year cold, heat, and electrical loads of the farm's building, based on its energy analysis software eQUEST (DOE-2.2,LBNL, Berkeley, CA, USA), as shown in Figures 3-5.The situation of the farm and biogas are shown in Table 2 and Table 3, respectively.This study establishes a photovoltaic cell module model using the TRNSYS software simulation platform.The annual output power of the PV in the climate environment of Jinan, in China is shown in Figure 6.In accordance to Figures 3-5, the time distribution characteristics of the building load can be obtained.

•
The loads of the building have obvious seasonal characteristics.According to the time distributions of electric, cooling, and heating loads, it can be divided into three typical seasonal load conditions corresponding to the winter, summer, and transition season.

•
The quarterly alternation exhibits no significant effect on the distribution of electric loads, but the electric load variation is associated with distinct time periods.The specific performance is assessed assuming that the electric load remains unchanged during the working hours, and the electric load is significantly reduced during the nonworking hours and holidays.

•
The heating load is mostly increased in the winter and transition seasons, and there is no heating load demand in the summer.
It should be pointed out that in order to avoid repeated load calculations, the electrical load of the building simulated in this study only included electricity consumed by lighting and electric machinery and does not include the electric energy load generated by driving the electric refrigerator.

Simulation Parameter Settings
Known conditions, such as the characteristic parameters of RES-CCHP and SP system, are shown in Table 4.The SP system, which serves as a reference in the optimized configuration only contains the electric chiller and the auxiliary boiler, and its capacity could be determined according to the maximum electric and thermal load throughout the year.Equipment costs of RES-CCHP and SP systems are listed in Table 5.In this study, biomass gas was used as the primary energy in the system.Compared with natural gas, various pollutant emission coefficients are different, as shown in Table 6 [31,37].The electricity price has a significant influence on the economic efficiency of the system.Different electricity prices are adopted in different regions, which would affect the operation of the system and the capacity allocation.This study fully considered the electricity price structure of Jinan in the test site and adopts three types of electricity price structures: peak, valley, and flat.The time-of-use price policy is currently applied in Jinan, China, where three price stages matched with peak-flat-valley value periods.The time-of-use price and timetable are shown in Table 7 [22,38].The weight coefficients corresponding to each index are shown in Table 8.The GA parameters and variable search range are shown in Tables 9 and 10, respectively.

Simulation Results Analysis
Analyzed results are presented according to the optimization flow, based on the simulated annual load data obtained in Jinan.These were obtained in combination with the above parameter settings and in accordance to the optimization of the MATLAB software platform used to solve the equipment capacity and key operating parameters, as shown in Table 11.
The performance values of the optimized configuration are shown in Table 12.
From the above data, we can see that the values of primary energy saving ratio, annual total cost saving rate and carbon emission reduction ratio are 20.94%,11.73% and 40.79%, respectively.These values fully reflect that the optimization design method can advance the RES-CCHP system in energy conservation, cost saving and emission reduction.Meanwhile, the RES-CCHP system responds to the country's new energy policy, and if the state subsidies are considered, the economic performance will improve.In order to further analyze the operating characteristics of the optimally configured RES-CCHP system at different time periods, the indicators of the co-supply system were analyzed on a monthly basis.As shown in Figure 7, the changing trends of the environment and energy indicators were consistent with the economic indicator, thus showing that the overall trends of the lower winter conditions and higher summer conditions.The main reason was attributed to the significant increase of the light intensity in the summer, whereby the output power of the PV had increased significantly.This in turn significantly reduced the amount of electricity purchased by the system and the power generated by the generator set, thereby reducing CO 2 emissions and the primary energy consumption.

Item
Value Primary energy saving ratio 20.94% Annual total cost saving rate 11.73% Carbon emission reduction ratio 40.79%Mo 20.66% From the above data, we can see that the values of primary energy saving ratio, annual total cost saving rate and carbon emission reduction ratio are 20.94%,11.73% and 40.79%, respectively.These values fully reflect that the optimization design method can advance the RES-CCHP system in energy conservation, cost saving and emission reduction.Meanwhile, the RES-CCHP system responds to the country's new energy policy, and if the state subsidies are considered, the economic performance will improve.
In order to further analyze the operating characteristics of the optimally configured RES-CCHP system at different time periods, the indicators of the co-supply system were analyzed on a monthly basis.As shown in Figure 7, the changing trends of the environment and energy indicators were consistent with the economic indicator, thus showing that the overall trends of the lower winter conditions and higher summer conditions.The main reason was attributed to the significant increase of the light intensity in the summer, whereby the output power of the PV had increased significantly.This in turn significantly reduced the amount of electricity purchased by the system and the power generated by the generator set, thereby reducing CO2 emissions and the primary energy consumption.The electric load data of buildings throughout the year is shown in Figure 3.In addition, the electricity consumption of the electric chiller should be considered in the summer.In Figure 8, the green line represents the electricity obtained from the PV, the blue line represents the electricity output of the PGU, and the red line shows the electricity from the grid.During the daytime, especially when the sunlight is particularly strong in the summer, most of the electric load is satisfied by the PV.At the same time, the output of PGU is less.At night, or when the light is weak, it is supplemented by PGU.If the ratio of the required electricity to the rated capacity of PGU is less than the optimized key value (θ = 0.2), the PUG stopped working and all the electricity was The electric load data of buildings throughout the year is shown in Figure 3.In addition, the electricity consumption of the electric chiller should be considered in the summer.In Figure 8, the green line represents the electricity obtained from the PV, the blue line represents the electricity output of the PGU, and the red line shows the electricity from the grid.During the daytime, especially when the sunlight is particularly strong in the summer, most of the electric load is satisfied by the PV.At the same time, the output of PGU is less.At night, or when the light is weak, it is supplemented by PGU.If the ratio of the required electricity to the rated capacity of PGU is less than the optimized key value (θ = 0.2), the PUG stopped working and all the electricity was provided by the grid.On the contrary, the situation that the required electricity was larger than the rated electrical capacity of PGU will cause the shortfall electricity, and that will occur in the summer days with relatively high cooling demand.The surplus electricity demand will be purchased from the PG. Figure 9 shows the heat distribution of the system.The heat load of the building is represented by the black line, and the blue line is the inadequate heat from the boiler which primarily appears in the winter when the heating load exceeds the residual heat of the PGU, while in the summer days, the absorption chiller cannot provide enough cooling load.Concluded from the results, the equipment capacity and operating parameters obtained by using the optimization design method can match efficiently with the energy demands.The system can make full use of solar energy and biomass energy, and can benefit from the time-of-use electricity price, resulting in the reduction of the operating cost.
provided by the grid.On the contrary, the situation that the required electricity was larger than the rated electrical capacity of PGU will cause the shortfall electricity, and that will occur in the summer days with relatively high cooling demand.The surplus electricity demand will be purchased from the PG. Figure 9 shows the heat distribution of the system.The heat load of the building is represented by the black line, and the blue line is the inadequate heat from the boiler which primarily appears in the winter when the heating load exceeds the residual heat of the PGU, while in the summer days, the absorption chiller cannot provide enough cooling load.Concluded from the results, the equipment capacity and operating parameters obtained by using the optimization design method can match efficiently with the energy demands.The system can make full use of solar energy and biomass energy, and can benefit from the time-of-use electricity price, resulting in the reduction of the operating cost.

Sensitivity Analysis
Since the system designed in this study utilized solar and biomass energy, the randomness and volatility encountered in the process of renewable energy utilization were inevitable (in which biomass can utilize constant temperature fermentation technology to achieve stable gas production).Thus, the climate, and in particular the prediction of the illumination intensity throughout the year, had certain influences on the system design.As shown in Figure 10, the illumination intensity changed from 60% to 150% of the test data (shown in Figure 6).Each interval is 10%.It can be seen that the number of optimized PV increased with increasing light intensity, from 159 to 245.The carbon dioxide emission reduction rate and primary energy saving rate show an increasing trend.Among them, the primary energy saving rate increased from 14.21% to 29.70% and the economic indicators increased from 6.17% to 16.71%.Meanwhile, the rate of reduction of carbon dioxide emissions increased from 39.06% to 44.07%.
The changes of biogas production will also influence the system performance.Therefore, it is necessary to study how the biogas production affects the system performance.As shown in Figure 11, the daily production of the biogas changed from 140 to 240 m 3  provided by the grid.On the contrary, the situation that the required electricity was larger than the rated electrical capacity of PGU will cause the shortfall electricity, and that will occur in the summer days with relatively high cooling demand.The surplus electricity demand will be purchased from the PG. Figure 9 shows the heat distribution of the system.The heat load of the building is represented by the black line, and the blue line is the inadequate heat from the boiler which primarily appears in the winter when the heating load exceeds the residual heat of the PGU, while in the summer days, the absorption chiller cannot provide enough cooling load.Concluded from the results, the equipment capacity and operating parameters obtained by using the optimization design method can match efficiently with the energy demands.The system can make full use of solar energy and biomass energy, and can benefit from the time-of-use electricity price, resulting in the reduction of the operating cost.

Sensitivity Analysis
Since the system designed in this study utilized solar and biomass energy, the randomness and volatility encountered in the process of renewable energy utilization were inevitable (in which biomass can utilize constant temperature fermentation technology to achieve stable gas production).Thus, the climate, and in particular the prediction of the illumination intensity throughout the year, had certain influences on the system design.As shown in Figure 10, the illumination intensity changed from 60% to 150% of the test data (shown in Figure 6).Each interval is 10%.It can be seen that the number of optimized PV increased with increasing light intensity, from 159 to 245.The carbon dioxide emission reduction rate and primary energy saving rate show an increasing trend.Among them, the primary energy saving rate increased from 14.21% to 29.70% and the economic indicators increased from 6.17% to 16.71%.Meanwhile, the rate of reduction of carbon dioxide emissions increased from 39.06% to 44.07%.
The changes of biogas production will also influence the system performance.Therefore, it is necessary to study how the biogas production affects the system performance.As shown in Figure 11, the daily production of the biogas changed from 140 to 240 m 3

Sensitivity Analysis
Since the system designed in this study utilized solar and biomass energy, the randomness and volatility encountered in the process of renewable energy utilization were inevitable (in which biomass can utilize constant temperature fermentation technology to achieve stable gas production).Thus, the climate, and in particular the prediction of the illumination intensity throughout the year, had certain influences on the system design.As shown in Figure 10, the illumination intensity changed from 60% to 150% of the test data (shown in Figure 6).Each interval is 10%.It can be seen that the number of optimized PV increased with increasing light intensity, from 159 to 245.The carbon dioxide emission reduction rate and primary energy saving rate show an increasing trend.Among them, the primary energy saving rate increased from 14.21% to 29.70% and the economic indicators increased from 6.17% to 16.71%.Meanwhile, the rate of reduction of carbon dioxide emissions increased from 39.06% to 44.07%.
The changes of biogas production will also influence the system performance.Therefore, it is necessary to study how the biogas production affects the system performance.As shown in Figure 11, the daily production of the biogas changed from 140 to 240 m 3 /d.Each interval is 10 m 3 /d.It can be seen that the three performances show increasing trends when the daily production was less than 200 m 3 /d.However, the ATCSR decreased when the daily production exceeded 200 m 3 /d.That's because the high daily output that cannot be utilized completely means huge investment.
Therefore, the fluctuation of the illumination and the changes of biogas production had minor effects on the optimization design method proposed in this study, which shows that the optimization method is robust.
Energies 2018, 11, x FOR PEER REVIEW 17 of 22 be seen that the three performances show increasing trends when the daily production was less than 200 m 3 /d.However, the ATCSR decreased when the daily production exceeded 200 m 3 /d.That's because the high daily output that cannot be utilized completely means huge investment.Therefore, the fluctuation of the illumination and the changes of biogas production had minor effects on the optimization design method proposed in this study, which shows that the optimization method is robust.

Conclusions
This study presented a RES-CCHP system driven by a biogas-fueled ICE and PV panels, which capitalized on the respective advantages of the RES and CCHP systems.Comprehensive consideration of renewable energy and load characteristics, four-decision variables (the capacity of  be seen that the three performances show increasing trends when the daily production was less than 200 m 3 /d.However, the ATCSR decreased when the daily production exceeded 200 m 3 /d.That's because the high daily output that cannot be utilized completely means huge investment.Therefore, the fluctuation of the illumination and the changes of biogas production had minor effects on the optimization design method proposed in this study, which shows that the optimization method is robust.

Conclusions
This study presented a RES-CCHP system driven by a biogas-fueled ICE and PV panels, which capitalized on the respective advantages of the RES and CCHP systems.Comprehensive consideration of renewable energy and load characteristics, four-decision variables (the capacity of

Conclusions
This study presented a RES-CCHP system driven by a biogas-fueled ICE and PV panels, which capitalized on the respective advantages of the RES and CCHP systems.Comprehensive consideration of renewable energy and load characteristics, four-decision variables (the capacity of power generation unit (PGU), the number of PV panels, the on-off coefficient of PGU, and the ratio of electric cooling to cooling load) were selected to be optimized in order to achieve collaborative optimization.A multi-objective optimization model was proposed that aimed to maximize the PESR, ATCSR, and CERR of the RES-CCHP system.Given the proposed optimization model, a GA was adopted to determine the capacity and operation of the RES-CCHP system.
A case study of a small farm in Jinan, China, was conducted to ascertain the feasibility of the proposed RES-CCHP system structure and the corresponding optimal operation strategy.The results illustrated that the implementation of the optimization design was beneficial to energy conservation, cost saving and emission reduction.The values of PESR, ATCSR, and CERR are 20.94%,11.73% and 40.79%, respectively, which demonstrates the advantages of the RES-CCHP system are clear compared with the SP system.The annual energy distribution results show that equipment capacity and operating parameters obtained by using the optimization design method can match efficiently with the energy demands and reduce the operating cost.Finally, the influences of the volatility of solar energy and the changes of biogas production on the optimization design method are analyzed, which show that the RES-CCHP system and the method proposed are robust.
The RES-CCHP system can effectively stabilize the randomness and volatility of the solar energy and biomass energy through the optimization design proposed in this study.It can significantly improve the energy efficiency and the solar and biomass energy utilization rate, and is suitable for supplying energy to a decentralized rural area.
To capitalize on the advantages of the RES-CCHP system in energy saving and emission reduction, and to improve the economic index, future research will combine the system design optimization and the output plan of key equipment in accordance to load changing, so as to further enhance the system performance.

Figure 1 .
Figure 1.Structure and energy flows of renewable energy sources integrated combined cooling, heating, and power (RES-CCHP) system.

Figure 2 .
Figure 2. Flow diagram of the execution of the genetic algorithm (GA).

Figure 2 .
Figure 2. Flow diagram of the execution of the genetic algorithm (GA).

Figure 2 .
Figure 2. Flow diagram of the execution of the genetic algorithm (GA).

Figure 3 .
Figure 3. Annual electricity load of building.

Figure 3 .
Figure 3. Annual electricity load of building.

Figure 4 .
Figure 4. Annual heating load of building.

Figure 5 .
Figure 5. Annual cooling load of building.

Figure 6 .
Figure 6.Annual output of single photovoltaic plane.

Figure 4 .
Figure 4. Annual heating load of building.

Figure 5 .
Figure 5. Annual cooling load of building.

Figure 6 .
Figure 6.Annual output of single photovoltaic plane.

Figure 6 .
Figure 6.Annual output of single photovoltaic plane.

Figure 7 .
Figure 7. Monthly performance of the RES-CCHP system.

Figure 7 .
Figure 7. Monthly performance of the RES-CCHP system.

Figure 8 .
Figure 8. Annual electricity distribution of the RES-CCHP system.

Figure 9 .
Figure 9. Annual heat distribution of the RES-CCHP system.

Figure 8 .
Figure 8. Annual electricity distribution of the RES-CCHP system.

Figure 8 .
Figure 8. Annual electricity distribution of the RES-CCHP system.

Figure 9 .
Figure 9. Annual heat distribution of the RES-CCHP system.

Figure 9 .
Figure 9. Annual heat distribution of the RES-CCHP system.

Table 1 .
Performance factors of naturally aspirated internal combustion engine (ICE).

Table 2 .
Description of the simulated building.

Table 3 .
Production information of the biogas.

Table 2 .
Description of the simulated building.

Table 3 .
Production information of the biogas.

Table 2 .
Description of the simulated building.

Table 3 .
Production information of the biogas.

Table 4 .
Characteristic parameters of the renewable energy sources integrated combined cooling, heating, and power (RES-CCHP) system and SP system.

Table 5 .
Unit prices of used equipment.

Table 7 .
Time-of-use electricity prices.

Table 11 .
Optimal results of the RES-CCHP system.

Table 12 .
Annual performance value of RES-CCHP system.

Table 12 .
The performance values of the optimized configuration are shown in Table12.Annual performance value of RES-CCHP system.