#### Abstract

A mass of thermal energy from the combined cooling, heating, and power (CCHP) system will be wasted by the buildings with low ratio of heat to power demand. To further improve the using efficiency of the excess thermal energy, an active method, which adds the organic Rankine cycle (ORC) equipment in the CCHP system, is configured. A complementary configuration of this CCHP-ORC is also investigated. This paper presents a matrix modeling approach to establish a mathematical model of the CCHP-ORC system. The CCHP-ORC system can be viewed as a multiple input and multiple output (MIMO) model. Energy conversion of the system components is described as the efficiency matrices. Energy flow between the system components is described as dispatch matrices. Energy conversion and flow from the system input to output are modeled by a conversion matrix. The objective function and the constraints of the system are determined. The optimal operation strategies are obtained by solving the optimization problem of minimizing the evaluation criteria function. Comprehensive case studies are conducted based on a hypothetical CCHP-ORC system. The study results reveal that the effectiveness and economic efficiency of the proposed approaches outperform those obtained from conventional CCHP systems.

#### 1. Introduction

According to the statistical data from the Energy Information Administration of the U.S. Department of Energy, in many countries, building energy consumption accounts for about 36% of the total energy consumption, and nearly 30% of greenhouse gas (GHG) emissions [1]. Particularly in China, the building energy consumption has an annual growth rate over 10% with the rapid development of economy [2, 3]. The increasing building energy consumption has led to serious building-related problems such as restricting economic development, exacerbating energy crisis, and environmental issues such as greenhouse effect and air pollution. To deal with these problems, combined cooling, heating and power (CCHP) systems which can improve energy utilization efficiency and reduce GHG emissions are developed.

A CCHP system is an integrated energy supply system, which can also be referred to as the trigeneration system [4]. It provides building users electricity, cooling energy, and heating energy by installing the power generation unit (PGU) [5, 6]. The CCHP system makes use of the rejected heat from the PGU to generate thermal energy for cooling and heating demands [7, 8]. As a result, the overall fuel energy utilization efficiency of CCHP systems can reach 70%-90%, which is significantly higher than that of separation production (SP) systems [9, 10]. As shown in the literature, CCHP systems have been broadly applied in office buildings, hotels, hospitals, commercial malls, and other building types.

Performance of the CCHP system depends on the system structure and the operation strategy. In [11], the authors propose a new CCHP system structure with both absorption and electric chillers to provide more flexibility in supplying energy. Four typical CCHP systems with existing technologies for various size ranges are illustrated in [12]. With the development of the renewable energy, solar-integrated CCHP systems [13–17], SOFC-trigeneration systems [18, 19], and biomass-trigeneration system [20] are studied.

Once the system structure is determined and implemented, the operation strategy is the crucial factor for the CCHP system design. There are two popular operation strategies [21–24]: following the electric load (FEL) and following the thermal load (FTL). However, as illustrated in the literature [11] both of the FEL and FTL strategies inherently waste a certain amount of energy. Therefore, designing an optimal operation strategy to improve the energy utilization efficiency by reducing the energy waste is urged.

Operation strategy design relies on the selection of evaluation criteria, which should be aligned with the stakeholders’ objectives. Energy savings and environmental impacts are two popular aspects when constructing evaluation criteria. Fumo* et al.* proposed the evaluation criteria considering the primary energy consumption saving (PECS) and carbon dioxide emission reduction (CDER) [25, 26]. Li* et al.* established a mix-integer nonlinear programming (MINLP) model to evaluate economic efficiency and sensitivity of CCHP systems [27]. Wang* et al.* studied the energetic and environmental benefits of the CCHP system in comparison to SP systems based on a particle swarm optimization algorithm [28]. Jing* et al.* estimated the PEC and pollution emissions of the BCHP system in FEL and FTL operation strategies based on life cycle assessment (LCA) [29, 30].

For an existing CCHP system, when the electricity to thermal energy requirement ratio of building users is exactly equal to the output ratio of the PGU, the maximum efficiency of the CCHP system is achieved [31]. However, this ideal situation rarely happens due to the random variation of energy requirements. When the supply is insufficient, adding extra energy (fuel or electricity) is a direct treatment. However, when the supply exceeds the demand, the situation becomes complicated. In the case surplus electricity is generated, it can be sold back to power grid or converted to thermal energy. Similarly, the excess thermal energy can be stored in thermal storage units or converted to electricity by some special equipment, such as organic ranking cycle (ORC). In this paper, an active excess thermal energy treatment method which uses an ORC to convert the excess thermal to electricity is investigated.

In recent years, more attention has been put on the ORC system. A complementary structure of CCHP-ORC system is configured in [31, 32]. A biomass fuelled trigeneration system integrated with ORC is studied in [33]. Three trigeneration systems, combined SOFC with ORC, combined biomass combustor with ORC, and combined solar collectors with ORC, are analyzed in [34]. Liu* et al.* researched five kinds of building in Shanghai and evaluated the CCHP-ORC system [35]. The aforementioned literature shows that the ORC technology has advantages in using the thermal energy and its combination with trigeneration system has good application prospects. So we are going to adopt the CCHP-ORC structure which uses the gas turbine as PGU, an ORC as waste heat recycle equipment, and electric and absorption chillers as hybrid chillers. The advantages of our CCHP-ORC system are reflected by the higher energy utilization, and the adjustability of the electricity to thermal energy output ratio.

When the structure of the CCHP-ORC system is configured, the operation strategy should be designed. Firstly, the mathematical model of the CCHP-ORC system should be established. From the perspective of control theory, the CCHP-ORC system can be treated as the connection between the energy input and the building users’ demand. We consider using transfer function in control theory to describe this input-output relation. Because the energy input includes fuel and electricity from the power grid, and the building users’ demand includes cooling, heating, and power, the CCHP-ORC system is a multiple input and multiple output (MIMO) system. In [36], the concepts of dispatch factors and coupling matrix are used to describe the system, and the power flow and operation strategy are optimized using Karush-Kuhn-Tucker (KKT) conditions. In [37], the authors model the system by introducing the concepts of bifurcations, junctions, and the backtracking. However these modeling methods are too complicated.

The main contribution of this paper is that a more intuitive and comprehensive approach, a matrix transfer function modeling approach (referred to as matrix modeling approach), is proposed to model the CCHP-ORC system. Based on the matrix modeled system, sequential quadratic programming (SPQ) algorithm is adopted to optimize the linear objective function, which is subjected to the nonlinear equalities and inequalities. The result of the optimization problem is the optimal flow and operation strategy for the CCHP-ORC.

This paper is organized as follows. The configuration and matrix modeling of the CCHP-ORC system are described in Section 2. In Section 3, the performance criteria for the CCHP-ORC system are presented and the optimization problem is established. Finally, the effectiveness of the system configuration and operation strategy are verified through a case study of a hypothetical office building in Section 4. Section 5 concludes this paper.

#### 2. System Configuration and Matrix Modeling

In this section, an optimal CCHP-ORC structure and an intuitive matrix modeling of the CCHP-ORC system will be introduced. In order to make full use of energy and improve the flexibility of energy supply, the configuration of the CCHP-ORC system is presented in Figure 1.

In addition to traditional equipment, a complementary ORC and an electric chiller are added as important components which can not only utilize the excess thermal energy but also dynamically adjust the electricity to thermal energy output ratio of the entire system. In Figure 1, when , the excess thermal energy exists. The ORC is used to convert the excess thermal to electricity. The cooling load is provided by absorption and electric chillers. The electric chiller has a high coefficient of performance (COP) which leads to a high cooling efficiency. The output of the ORC and the electric chiller can be adjusted dynamically to keep the CCHP-ORC system running at the matched status according to the users’ energy requirements.

##### 2.1. Efficiency Matrices of System Components

The energy conversion of the components in the CCHP-ORC system can be described using efficiency matrices, namely, coupling matrices. In this paper, is an equipment set and equipment , , , , and denote the fuel, electricity, cooling energy, and heating energy, respectively. The input and the output of the component are defined as and , respectively. So the input-output relation of the component can be described aswhere is the efficiency matrix of the component .

Consequently, the input-output relation of the ORC can be represented in the following matrix form:

Following the same procedure, efficiency matrices for the PGU, auxiliary boiler, heat recovery system, heating unit, electric chiller, and absorption chiller can be obtained as

Each efficiency matrix represents the performance of the energy conversion of each component in the CCHP-ORC system.

##### 2.2. Dispatch Matrices of System

The dispatch matrices indicate the energy flow between the components in the system. Usually, the dispatch factors exist at the bifurcations of the system. For example in Figure 1, the thermal energy produced by the heat recovery system and the auxiliary boiler will be separated into three parts: one for the absorption chiller, one for the heating unit, and the last one for the ORC. Let , , and denote the dispatch factors for the absorption chiller, the heating unit, and the ORC, respectively. Then we havesubject to

At the ORC side, we havewhere and denote the dispatch factors for the PGU and the auxiliary boiler. To facilitate the later optimization-based operation strategy design, (8) can be rewritten in a matrix form aswhere is the dispatch matrix for the ORC and is the system input.

By following the similar procedure, the dispatch matrices for the PGU, the auxiliary boiler, electric chiller, heat recovery system, heating unit, and absorption chiller can be defined as

Similarly, in addition to (7a), (7b), (7c), and (7d), we have

##### 2.3. Conversion Matrix of the CCHP-ORC System

The CCHP-ORC system modeled in a matrix form can be viewed as an input-out model. The conversion matrix of the CCHP-ORC system describes the relationship between the input and the output and represents the energy flow procedure of the system. From Figure 1, the system input and output are defined asrespectively.

Then the conversion matrix of the CCHP-ORC system, , can be defined as

In order to obtain the conversion matrix , the expression of the output elements, , , and , should be written as functions of the system input. Without the loss of generality, the parasitic electricity can be assumed to be a constant and further assumed to be zero. The output element can be represented as

For the cooling part, we havewhere the electric cooling to cool load ratio, which is inherently included in the modeling procedure and can be obtained through the optimization, is defined as .

The heating demand is solely provided by the heating unit. Thus, we have

The conversion matrix of the CCHP-ORC system can be obtained from (15)–(17) aswhere

The main objective of this paper is to determine the dispatch factors and the system input to optimize the performance criteria objective function.

#### 3. Optimal Operation Strategy Design

Having the whole system modeled in a matrix form, we can formulate the operation strategy design as an optimization problem. The objective function and associated constraints will be illustrated in the rest of this section.

##### 3.1. Performance Criteria for the CCHP-ORC System

As mentioned in Section 1. Evaluation of a CCHP system may contain many factors, including primary energy consumption, energy conversion efficiency, total cost of the system, environmental impact, life cycle assessment, and combination of multiple factors. In this paper, three popular criteria, i.e., primary energy consumption, the total cost of the system, and the environmental impact, are adopted to evaluate the performance of the CCHP-ORC system.

###### 3.1.1. Primary Energy Consumption (PEC)

PEC is a direct and common indicator to reflect the operating efficiency of a CCHP-ORC system. PEC is defined aswhere is the parasitic electricity, is the generation efficiency of the grid power plant, and is the transmission efficiency of grid.

###### 3.1.2. Hourly Total Cost (HTC)

In order to obtain the optimal energy flow and operation strategy for every hour, HTC is adopted as a key element in the operation strategy optimization. HTC of the CCHP-ORC system is defined aswhere and are the unit prices of electricity and the fuel, respectively; and are the carbon conversion factors of the the electricity and fuel, respectively; is the carbon tax rate; is the installed capacity of each equipment and is the initial capital cost of each equipment; is the number of equipment included in the CCHP-ORC system; and is the life of each equipment, where we assume that all equipment have the same life time. 8760 (h) represents the hours number of the whole year.

###### 3.1.3. Carbon Dioxide Emission (CDE)

Greenhouse gas emission control is an emerging topic in the CCHP system design. In a CCHP-ORC system equipped with fossil-fuel combustion units, CDE serves as an indicator of environment impacts and can be calculated aswhere and are the carbon conversion factors of the the electricity and fuel, respectively.

###### 3.1.4. Performance Evaluation Criteria (EC) Function

To comprehensively evaluate the system performance, we define a performance evaluation criteria (EC) function as a weighted summation of the three above mentioned criteria. The most commonly used way to evaluate the performance of a CCHP-ORC system is to compare those factors with those of the SP system. Thus, the hourly EC function can be redefined aswhere , , and are the weights of PEC, HTC, and CDE, respectively. Three weights satisfy, , and are PEC, HTC, and CDE of the SP system and can be, respectively, calculated as (20), (21), and (22).

Consequently, the annual function can be written aswhere is the hourly function value on day at hour .

##### 3.2. Optimization Problem Formulation

Formulating a complete optimization problem includes determining decision variables, designing objective function, and setting constraints. The whole procedure will be elaborated on this part.

###### 3.2.1. Decision Variables

As mention in Section 2, dispatch factors and energy inputs are the variables that need to be adjusted to coordinate the energy flow of the system, meet the building demand, and minimize the objective function. Thus, the dispatch factors and system energy inputs can be selected as decision variables.

It can be readily observed from (7a), (7b), (7c), (7d) and (11a), (11b), (11c), (11d) that the dispatch factors , , , , , , and are dependent variables. Thus, the dispatch factors can be reduced to be , , , and , which are independent of each other.

The system energy inputs include and . Thus, the optimizer consisting of all the above-mentioned decision variables can be defined as

###### 3.2.2. Objective Function

By using the defined optimizer , the function in (23) can be written as a linear function of the optimizer aswhere

will serve as the objective function of the optimization-based operation strategy design. , , and can be set according to different evaluation requirements. For example, would be set relatively higher if more emphases are put on the environmental aspect.

###### 3.2.3. Equality Constraint

In the sense of balance, (14) will serve as the equality constraint of the optimization problem. Further, since some elements in are represented by the decision variables, it is necessary to write as a function of . Thus we havewhereand

Similarly, the system input should also be represented by aswhere

Then the nonlinear equality constraint in (14) becomes

###### 3.2.4. Inequality Constraints

The inequality constraints include the bound of the decision variables, various capacities of system components, and thresholds of the components output. The decision variables can be separated into two parts: one consists of dispatch factors, and the other consists of system input. To facilitate the later derivation, we define dispatch factor vector as

According to (7a), (7b), (7c), (7d) and (11a), (11b), (11c), (11d), we have

In this paper, we assume that no electricity would be sold back; thus the constraints for fuel and electricity input can be represented by

Inequalities in (36a), (36b), and (38) are linear inequality constraints seeing as the bounds of the decision variables.

Each component in the CCHP-ORC system has a specific capacity, which is the output upper bound. Denote the capacity of ORC; then the upper bound of the ORC output can be represented as

The output upper bound of the PGU, heat recovery system, auxiliary boiler, electric chiller, absorption chiller, and heating unit can be denoted as , , , , , and and defined in a similar way as in (39), respectively.

In addition, the lower bound of the ORC, PGU, heat recovery system, auxiliary boiler, electric chiller, absorption chiller, and heating unit can be denoted as , , , , , , and . The lower bound of a component implies that the output will be cut down when the expected output of the component is lower than this bound. In this paper, the lower bound can be set as 5 percent of the maximum capacity of each component.

Based on the above discussions, the bounding conditions for component output can be represented by

Then from (9) and (32), (40a), (40b) can be rewritten aswhere is the dispatch matrix for the component in the CCHP-ORC system and can also be represented by the function of by following the similar procedure of deriving (29) aswhereand , , and .

Therefore, (36a), (36b), and (38) are the linear inequality constraints, and (41a) and (41b) are the nonlinear inequality constraint of the optimization problem. The optimal solution can be obtained by solving

The objective function integrates three aspects, i.e., PEC, HTC, and CDE, into it, and it is comprehensive. Three weights , , and , can be selected according to different requirements. For example, if more priority is need to be given to environmental factors, the carbon dioxide emission (CDE) should be put on more emphasis, then is raised, and and are accordingly decreased.

#### 4. Case Study

##### 4.1. Description of a Hypothetical Building

In this section, a CCHP-ORC system is installed for a hypothetical office building in Wuhan, Hubei, China, because the CCHP-ORC system is more suitable to the commercial building and the office building with low ratio of heat to power demand [35]. We chose EnergyPlus [38, 39] to analyze energy consumption of this building. This building is assumed to operate all year round and has five floors with a total construction area of . The first floor consists of exhibition halls and office rooms. The second to the fifth floors are office rooms. Some important construction parameters of the hypothetical office building are listed in Table 1.

The hourly cooling, heating, and power requirements in one year of the hypothetical office building are shown in Figure 2.

To ease the performance comparison, the hourly cooling, heating, and power requirements in representative days of spring, summer, autumn, and winter are shown in Figures 3–6, respectively.

From Figures 2–6, it can be noticed that the cooling requirements in summer are the most, the heating requirements in winter are more than those of other seasons, and the electricity requirements are substantially the same for each season.

##### 4.2. Simulation Parameters

Table 2 shows a set of technical parameters of the installed CCHP-ORC system for the hypothetical building. In this paper, the fuel is the natural gas which is widely used in China. The gas turbine is chosen as PGU. For performance comparison, some parameters of a conventional CCHP system, coefficient of performance, and efficiency of the components in the CCHP-ORC system, conversion factors, and weight factors are included in Table 2. The weighting factors of the evaluation criteria can be chosen to satisfy different requirements. Some parameters listed in Table 2 fluctuate with the market, such as price of electricity, natural gas, and carbon tax. Most of the system parameters can be gotten from equipment manual and some references.

##### 4.3. Test Result

SQP algorithms have been proved highly effective for solving general constrained problems with smooth objective and constraint functions [40]. In this paper, we chose SQP algorithms to solve the optimization problem established in Section 3.2. In order to accelerate the convergence, some feasible initial points which are calculated from FEL and FTL strategies are manually designated. So the local optimal solution can be avoided and the convergence time is significantly reduced.

The constructed optimization problem is solved in MATLAB by using fmincon.m. The computer used to solve this problem is configured with 2.6 GHz Intel Core i5-3230M processor and 4 GB 1067 MHz memory. Time elapsed for each step of optimization is about . This time is much shorter than one hour; in other words, the optimal energy flow and operation strategy for the next hour can be quickly obtained at the beginning of this hour. Thus, we can say that the optimization method is appropriate for practical use.

The performance comparison of FEL, FTL, and the proposed optimal operation strategy is shown in Figures 7–10 which includes spring, summer, autumn, and winter, respectively. From Figures 8 and 10, all values of the optimal operation strategy are obviously less than those of the FEL and FTL strategies. From Figures 7 and 9, most values are also obviously less than those of FEL and FTL except 8:00-10:00. Because it is the beginning working time during 8:00-10:00, electricity requirements are very large, and there are almost no heating and cooling requirements. Taking into account the full day values, the performance of the optimal operation strategy is much better than that of FEL and FTL strategies. In spring, summer, and autumn, especially in the morning and in the evening, the optimal energy flow has excellent performance, while in the middle of a day, the advantage is not that evident. In winter, the performance in the middle of a day is much better than that in the morning and in the evening.

The reason for the proposed optimal operation strategy performing better than FEL and FTL strategies is that both FEL and FTL strategies are inevitably wasting energy. When operating at FEL mode, electricity for the building is firstly satisfied, and then if the thermal energy provided by PGU cannot meet the thermal demand, the auxiliary boiler will be activated; if the thermal supply exceeds the thermal demand, the excess thermal energy will be wasted. When operating at FTL mode, the thermal requirements are firstly satisfied, and then if the electricity provided by PGU cannot meet the electricity demand, additional electricity should be purchased from the power grid; if the electricity provided by PGU is sufficient, a certain amount of electricity will be wasted.

The CCHP-ORC system includes ORC component as a power generator to make full use of excess thermal energy, and its integration with CCHP system has good application prospects. According to the energy requirements in representative days shown in Figures 3–6, the hourly outputs of the electric chiller and the ORC in representative days can be calculated as Figures 11–14.

From Figures 11–14, we can find that the electric chiller has high utilization rate in summer because of the high temperature in Wuhan, and the ORC plays a significant role in other three seasons.

Substituting the primary energy consumptions, the daily total cost, and the carbon dioxide emission into (20)–(22), the daily performance criteria (i.e., PEC, DTC, CDE) for the CCHP-ORC and the CCHP system can be calculated. Capacities of all equipment in the CCHP-ORC and CCHP system are listed in Table 3.

In Table 3, the capacity of PGU in the CCHP system is 100 kW, which is searched in with the whole year (8760h) data using enumeration algorithm. In the CCHP-ORC system, ORC is combined as power generator. In order to satisfy the power requirements, the sum capacity of PGU and ORC should be equal 100kW which is the optimal capacity in the CCHP system. In addition, efficiency of PGU is higher than that of the ORC. For the office building with low ratio of heat to power demand, power requirements are more than other requirements, so the capacity of PGU should be far larger than that of ORC. Through enumeration, the optimal capacity of PGU and ORC is 80kW and 20kW, respectively.

The performance comparisons of the CCHP-ORC system and the CCHP system are shown in Table 4. It is clear that all the performance values of the CCHP-ORC system are smaller than those of the CCHP system and the performances of the CCHP-ORC system are better than those of the CCHP system. The integration of ORC and PGU can improve the energy efficiency and reduce the capacity requirements of the auxiliary boiler and PGU. With the reduction of the ORC’s price, the superiority of the CCHP-ORC system will be more and more outstanding.

If we assume that the price of ORC is about 100000 Yuan, from Table 4, we only consider the PEC and DTC, we can find the whole year’s PEC saving is 58190.13 Yuan, the total consumption reduction is 16123.88 Yuan, and the sum of the whole year’s PEC saving and the total consumption reduction is 74314 Yuan. So the CCHP-ORC system runs 1 year and 4 months, and the cost of ORC is recovered.

#### 5. Conclusion

A complementary configuration of the CCHP-ORC system is investigated to further improve the utilization efficiency of the excess thermal energy from the CCHP system. A comprehensive and systematic matrix modeling approach is proposed to describe the CCHP-ORC system. The matrix model includes efficiency matrices of system components, system conversion matrix, and dispatch factors. The optimal operation strategy is obtained by solving the optimization problem. Both economic and environmental factors are considered in the objective function of the optimal operation strategy design. The weights of economic and environmental factors can be changed according to different requirements. At last, an illustrative case study is conducted to show that the effectiveness and economic efficiency of the proposed approaches outperform those obtained from conventional CCHP systems.

The main challenges of applying this method to real system are that some uncertainties in the CCHP-ORC system exist; for example, a likely source of uncertainty in the heat exchanger components of the CCHP-ORC system is fouling. The efficiencies of heat exchanger components are changed. In real system, we should consider these uncertainties by adding an error of less than 5% of the nominal value of the parameters.

The main problem of the CCHP-ORC system is high investment cost of the ORC. However the unit price of the ORC will continually decline because of the technology maturing in recent years. Therefore the CCHP-ORC system will be more attractive and valuable. In addition, all calculation in this paper is based on an accurate load profile. While, in practical applications, accurate load profile does not exist, therefore, in future work, accurate load forecasting should be considered. In a further step, excess electricity sold back feature can be considered. In recent years, with the development of the renewable energy including solar energy, wind energy, and tide energy, more and more renewable energy should be incorporated into the CCHP system to make the trigeneration system more efficient.

#### Nomenclature

CCHP: | Combined cooling, heating, and power |

CDE: | Carbon dioxide emission |

COP: | Coefficient of performance |

DTC: | Daily total cost |

EC: | Evaluation criteria |

FEL: | Following the electric load |

FTL: | Following the thermal load |

HTC: | Hourly total cost |

ORC: | Organic Rankine cycles |

PEC: | Primary energy consumption |

PGU: | Power generation unit |

SP: | Separation production. |

*Symbols*

C: | Cost |

E: | Electricity |

F: | Fuel |

L: | Facility life |

N: | Installed power |

Q: | Thermal energy |

: | Efficiency |

: | Carbon conversion factor. |

*Subscripts*

ac: | Absorption chiller |

annual: | Annual value |

: | Auxiliary boiler |

: | Cooling |

ca: | Carbon |

: | Electricity |

ec: | Electric chiller |

: | Fuel |

grid: | Local power grid |

: | Heating |

hour: | Hourly value |

hrc: | Recovered heat cooling |

hrh: | Recovered heat heating |

: | Total consumption |

orc: | Organic Rankine cycles |

: | Parasitic |

pgu: | Power generation unit |

pgum: | Capacity of PGU |

: | Recovered. |

*Superscripts*

: | SP system. |

#### Data Availability

All research data reported in this article are available. All of these data can be obtained free of charge. If you want to get these data, you can directly contact Yufang Chang; her email is [email protected].

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work has been supported in part by the Green Industry Leading Program of Hubei University of Technology (CPYF2017003), the National Nature Science Fund of China (61601176) and the National Nature Science Fund of China (41601394).