Research on Economic Operation Strategy of CHP Microgrid Considering Renewable Energy Sources and Integrated Energy Demand Response

: Renewable energy resources (RESs) play an important role in the upgrading and transformation of the global energy structure. However, the question of how to improve the utilization e ﬃ ciency of RESs and reduce greenhouse gas emissions is still a challenge. Combined heating and power (CHP) is one e ﬀ ective solution and has experienced rapid development. Nevertheless, with the large scale of RESs penetrating into the power system, CHP microgrid economic operation faces great challenges. This paper proposes a CHP microgrid system that contains renewable energy with considering economy, the environment, and system ﬂexibility, and the ultimate goal is to minimize system operation cost and carbon dioxide emissions (CO 2 ) cost. Due to the volatility of renewable energy output, the fuzzy C-means (FCM) and clustering comprehensive quality (CCQ) models were ﬁrst introduced to generate clustering scenarios of the renewable energy output and evaluate the clustering results. In addition, for the sake of improving the ﬂexibility and reliability of the CHP microgrid, this paper considers the battery and integrated energy demand response (IEDR). Moreover, the strategy choices of microgrid operators under the condition of grid-connected and islanded based on environment and interest aspects are also developed, which have rarely been involved in previous studies. Finally, this stochastic optimization problem is transformed into a mixed integer linear programming (MILP), which simpliﬁes the calculation process, and the results show that the operation mode under di ﬀ erent conditions will have a great impact on microgrid economic and environmental beneﬁts.


Introduction
With the increasing depletion of fossil energy and the deterioration of the environment becoming more and more serious, the development and utilization of RESs is the main direction of future energy development. Distributed energy based on RESs such as wind power, solar energy, and biomass widely exist as clean energy, and development technology is becoming more and more mature, which has received wide attention around the world. Many countries have set future renewable energy development goals to promote sustainable energy development [1].The share of RESs in Finland will increase to 38% by 2020 [1]. The Chinese government has formulated four energy revolutions, aimed at achieving the goal of nonfossil energy accounting for 15% and 20% of primary energy consumption in 2020 and 2030, respectively [2]. Meanwhile, in traditional power systems, a large part (40-60%) of fossil energy or renewable energy sources may be wasted because traditional power systems cannot effectively use the input energy [3]. Further, the fluctuation of renewable energy output, market electricity prices, and load demand make it hard to elevate comprehensive energy utilization efficiency and reduce greenhouse gas emissions. In this regard, one of the effective solutions to solve these problems is the combined heating and power (CHP) microgrid. As an independent and controllable energy supply system, CHP microgrids can satisfy local energy demands and guarantee the safe and stable operation of a power system, reduce extra costs, and improve the environment [4,5]. However, with the huge penetration of renewable energy, the operation mode of the CHP microgrid has changed dramatically, which needs a new mode of operation to match the current situation.
The economic dispatch of a CHP microgrid with renewable energy has been the focus of significant efforts in the literature. Reference [6] presents a multiobjective economic emission dispatch to minimize total operation cost and CO 2 emission cost for the microgrid that contains RESs. Reference [7] developed an optimization dispatching model for CHP plants to minimize system operation costs and satisfy the end-consumers' demands. For the CHP-direct heating(CHP-DH) system, the authors established energy balance, system control, and operation constraints in [8], where the optimized objective function is to minimize system operation cost and net acquisition for heat and power. A stochastic optimization model considefing the economic and environmental impacts and reliability is proposed in [9], whose ultimate goal is to minimize fuel cost. In order to solve the multiobjective economic dispatch of the microgrid, a particle swarm optimization (PSO) method and enhanced firefly algorithm (FA) have been adopted, respectively, in [10,11]. Reference [12] developed an optimization model from both the environmental and economic viewpoints and proposes optinal operation and management strategy for the mincrogrid. Reference [13] introduces an economic model and optimization method for the CHP system, including solar, wind, and fuel cell. The optimization goal is to minimize system operation costs and maintenance costs and compare the effects of genetic algorithms and particle swarm optimization algorithms in the model. In these references, although the fuel cost and the operation cost of the generator units are taken into account, the impact of environmental on the economic operation of CHP microgrids is rarely considered.
Further, several studies focus on the impact of auxiliary functions, such as demand response and batteries on CHP microgrid operation. Reference [14] satisfies the electrical and thermal load demands of end-users with minimal cost by implementing demand response (DR) programs. Reference [15] considers the energy storage system as a non-ideal battery and analyzes the impact of battery on the economic operation of CHP microgrid. Reference [16] considers DR as a virtual generator, replacing user-side resources with supply-side energy, which can effectively improve energy efficiency and reduce energy costs. References [17,18] believe that under fluctuation of electricity prices, real-time electricity prices can improve the adaptability of DR and thus propose a price-based DR management model for real-time electric and thermal loads. However, most studies about DR focus on the power load but neglect the heat load. Therefore, it is necessary to propose a novel DR model.
In addition, the impact of volatility of renewable energy on grid shocks has not been well resolved. Some of the literature has studied the uncertainty processing method. Reference [19] describes the uncertainty of wind power output through normal distribution but cannot accurately characterize wind power changes. A weighting method is applied to systematically combine multiple predictive models in [20], which can reduce the prediction bias. References [21][22][23][24] deal with wind power uncertainty through scenario generation and reduction methods. The solution to the uncertainty of renewable energy output are relatively simple and quite limited. It is urgent to find new ways to reduce the impact of renewable energy output uncertainty on the power system.
Through the analysis of the previous literature, we summarize the structure and influencing factors of the CHP microgrid system proposed in several typical studies and compare them with this paper, as shown in Table 1.
To deal with the problem mentioned above, we propose an optimal scheduling method based on stochastic optimization against the uncertainty of renewable energy output and consider the economy, environment, as well as system flexibility aspects in this paper. In addition, the operation mode of the CHP microgrid system and the behavior of system operators in different situations are discussed, and operational strategies are designed, respectively based on the interests of system operators and the environment. In short, the main contributions of this paper are included as follows.
(1) Considering system operation cost, carbon dioxide emissions cost, as well as system flexibility in the objective function. (2) Based on the fuzzy C-means (FCM) clustering method, the number of historical scenarios of wind power and solar energy output can be reasonably reduced to a fixed value. Then, a novel clustering evaluation method named CCQ is proposed to determine the number of best scenario categories, which can reduce the uncertainty of renewable energy output compared with previous literature methods. (3) A new demand side management method called integrated energy demand response is employed.
Compared with a single electrical or heat load demand response method, the flexibility of the CHP microgrid system is improved, and the system operation cost is reduced significantly. (4) The operation mode of the CHP microgrid and the behavior of the system operator in different situations are discussed, and different operational strategies are put forward in the light of the interests of system operators and environmental benefits.
The rest of proposed paper is organized as follows: Section 2 mainly describes the structure of the CHP microgrid and its mathematical model. Section 3 presents the objective functions and constraints of the CHP microgrid, and the solution to the proposed model. Case studies and simulation results as well as the corresponding operational strategy are explained in Section 4. Section 5 summarizes the main conclusions in this paper.

CHP Microgrid Structure Description and Mathematical Modelling
In this section, The CHP microgrid system is introduced, which contains CHP units, photovoltaic, wind turbine, and the integrated energy demand response, as well as mathematical model of units.

CHP System Composition and Structure
The structure of a CHP microgrid system in this paper is illustrated in Figure 1, which contains two energy flows (electrical energy flow and thermal energy flow). In order to describe the impact of battery, integrated energy demand response, and grid connection on the operation of the CHP microgrid, we considered a CHP microgrid that includes a micro gas turbine (GT), a wind turbine (WT) unit, a photovoltaic cell (PV), a battery (BT), a gas boiler (GB), and a waste heat boiler (HB), along with electricity and heating loads. model of units.

CHP System Composition and Structure
The structure of a CHP microgrid system in this paper is illustrated in Figure 1, which contains two energy flows (electrical energy flow and thermal energy flow). In order to describe the impact of battery, integrated energy demand response, and grid connection on the operation of the CHP microgrid, we considered a CHP microgrid that includes a micro gas turbine (GT), a wind turbine (WT) unit, a photovoltaic cell (PV), a battery (BT), a gas boiler (GB), and a waste heat boiler (HB), along with electricity and heating loads. In this CHP microgrid system, PV and WT provide electricity to the microgrid with renewable solar and wind energy, respectively. GT and GB use natural gas as fuel to convert chemical energy into electricity and thermal energy into the microgrid. For the sake of achieving the cascade utilization of energy and improving energy efficiency, a HB is used to recover the waste heat generated by the GT. The BT performs real-time charging and discharging according to the operating state of the CHP microgrid system, which can ensure the balance of supply and demand of the electricity in the system and promote the maximum consumption of solar energy and wind energy. In addition, the CHP microgrid system can also exchange power with the main network. When the system cannot meet the user's electrical load demand, the system operator can purchase power from the main network to meet the user load, otherwise, the system operation can sell excess electricity to the main network for additional profit. In this CHP microgrid system, PV and WT provide electricity to the microgrid with renewable solar and wind energy, respectively. GT and GB use natural gas as fuel to convert chemical energy into electricity and thermal energy into the microgrid. For the sake of achieving the cascade utilization of energy and improving energy efficiency, a HB is used to recover the waste heat generated by the GT. The BT performs real-time charging and discharging according to the operating state of the CHP microgrid system, which can ensure the balance of supply and demand of the electricity in the system and promote the maximum consumption of solar energy and wind energy. In addition, the CHP microgrid system can also exchange power with the main network. When the system cannot meet the user's electrical load demand, the system operator can purchase power from the main network to meet the user load, otherwise, the system operation can sell excess electricity to the main network for additional profit.

Micro Gas Turbine
The mathematical relationship between the energy output P t GT and the natural gas consumption F t GT can be modeled by the linear equation [25]: where a and b are gas-power conversion parameters. U t GT represents the operating state of the micro gas turbine in each dispatch cycle that is equal to 1 if the GT is on in period t and 0 otherwise. P t et , H t GT denotethe electrical energy and thermal energy generated by the gas turbine in interval t, respectively. Generally, there is a nonlinear relationship between input variables and output variables, which brings great difficulties to the solution of the problem. In order to optimize the solution, the thermal and electrical efficiency of the GT can be expressed as piecewise linearization, as shown in Figure 2, which can be modeled as mixed integer linear programming (MILP) [25].
where H t GT is the heat produced by the GT. N n GT is the coefficient of the thermal and electrical efficiency curve. M n GT represents the block limit of the thermal and electrical efficiency curve. C n GT represents the slope of block n of the thermal and electrical efficiency curve. Z tn GT is the binary variable encoding the thermal and electrical efficiency curve of the GT. l GT is the index set of the piecewise linear thermal and electrical efficiency curve. The mathematical relationship between the energy output t GT P and the natural gas consumption t GT F can be modeled by the linear equation [25]: where a and b are gas-power conversion parameters. t GT U represents the operating state of the micro gas turbine in each dispatch cycle that is equal to 1 if the GT is on in period t and 0 otherwise. Generally, there is a nonlinear relationship between input variables and output variables, which brings great difficulties to the solution of the problem. In order to optimize the solution, the thermal and electrical efficiency of the GT can be expressed as piecewise linearization, as shown in Figure 2, which can be modeled as mixed integer linear programming (MILP) [25].

Gas Boiler
where F t GB and H t GB denote natural gas consumption and the power output of the gas boiler in interval t, respectively. η GB denotes the efficiency of the gas boiler. LHV is the low calorific value of natural gas.

Waste Heat Boiler
where P t whb represents the heat recovered by the HB in period t. η whb represents the recovery efficiency of the HB, and P whb,min and P whb,max are the minimum and maximum power of the HB, respectively.

Exchange Power Between the Main Grid and the CHP Microgrid
where Pm 1 t and Pm 2 t represent electricity bought and sold to the main network in interval t, respectively.
U t buy and U t sell are the state of buying/selling, respectively, from/to main network in interval t.

Battery
where PC t and PD t denote charge/discharge power of battery in interval t, respectively. U t cha and U t diss represent the charge/discharge state of the battery in interval t, where if battery is charging (U t cha = 1) or discharging (U t diss = 1). δ es , ξ chr es , and ξ diss es refer to the energy loss rate and charge/discharge efficiency of the battery, respectively. SOC min and SOC max refer to the lower, upper limits of state-of-charge. S es is the rated capacity of the battery. S t es is the amount of battery stored in interval t.

Integrated Energy Demand Response Program
This paper introduces a new concept named integrated energy demand response (IEDR) to help the CHP microgrid system to reduce operating cost and increase system flexibility by changing the user's energy consumption pattern. In other words, the system operator can shift the load from peak time intervals to off-peak ones, as shown in Figure 3, to minimize the operating cost of the CHP microsystem as well as maximize the profit [14]. user's energy consumption pattern. In other words, the system operator can shift the load from peak time intervals to off-peak ones, as shown in Figure 3, to minimize the operating cost of the CHP microsystem as well as maximize the profit [26].

Electric Energy
where ELoad DR TL1 t represents the electrical load after the implementation of DR. ELoad t represents the actual electrical load without implementing DR. DR TL1 t indicates increased/decreased electrical load in DR. DRP max1 indicates the maximum limitation of the consumer's participation in DR.

Thermal Energy
where HLoad DR TL2 t is the thermal load after the implementation of DR. HLoad t is the actual thermal load without implementing DR. DR TL2 t represents increased/decreased thermal load in DR. DRP max2 represents the maximum limitation of the consumer's participation in DR.

CHP Microgrid Optimization Scheduling Strategy
Stochastic programming is a method for dealing with uncertain decision problems and has been widely used in the field of power systems. This paper takes the output power of CHP units and the renewable energy as the decision variables, via converting stochastic programming to mixed linear programming (MILP) and solving it with MATLAB R2016b (MathWorks company, Natick, MA, USA, 2016) and CPLEX Optimization Studio v12.8(IBM company, Amund, New York, NY, USA, 2017).

Objective Function
The operation cost of the CHP microgrid system, which contains the power exchange cost with the main network, the maintenance cost of each unit, the CO 2 emission cost, and natural gas procurement cost, should be minimized. The mathematical model of the objective function is shown below (14): where C exc is the exchange power cost. C gass is the natural gas consumption cost. C mgt is the gas turbine maintenance cost. C mpv , C mwp , C mgb , C mhb , and C mes represent the maintenance cost of each generating unit, respectively. C co 2 is the CO 2 emissions cost. C DR is the IEDR compensation cost.

The Cost of Exchange Power between the Main Network and the CHP Microgrid
where C buy and C sell are electricity procurement/sales costs from or to the main network. p t buy and p t sell are electricity procurement price and electricity sales price, respectively.
where p gas represents natural gas price. F t GT and F t GB represent the amount of natural gas consumed by gas turbines and gas boilers in interval t.

The Maintenance Cost of Each Unit
where R mgt , R mgb , R mhb , R mes , R mpv , and R mwt indicate the unit power maintenance cost of each unit in the CHP microgrid in interval t, respectively.

Carbon Dioxide Emission Cost
where P CO 2 represents carbon tax. C gt , C gb , and C grid indicate carbon dioxide emissions per unit of power from the GT, GB, and main network.

Integrated Energy Demand Response Compensation Cost
where P E and P H represent the compensation price of the unit electric load/thermal load in the demand response program.

CHP Microgrid System Constraints
This section mainly describes the CHP microgrid system constraints, including the balance of electrical energy and thermal energy, as well as the technical constraints of the units.

•
Electrical energy balance constraint • Thermal energy balance constraint • The technical constraints of the units.
The technical constraints of the units are described in Section 2.2 in detail. Please see Equations (1)-(11) for specific constraints.

RESs Uncertainty Sets
Since WT and PV power generation are subject to meteorological conditions, their power output has certain intermittency and volatility. In order to characterize the WT and PV, this paper adopted the improved fuzzy C-means (FCM) to cluster the historical power scenarios of wind power and solar energy out into several typical scenarios and then used the clustering comprehensive quality (CCQ) method to evaluate whether these scenarios are representative and determine the optimal number of scenarios [26].

FCM Clustering Method
In this paper, the objective function was established with all the same scenarios as the minimum total distance from the cluster center: where n represents the total number of scenarios; c is the number of clusters, u k j is the relative membership of the sample j to the category k; and Euclidean distance between sample j and category k; w i is the weight of the i-th indicator, and Usually, there are differences in the dimension and magnitude of each indicator in the system, so each indicator needs to be normalized, then d k j is updated to: where a max ij and a min ij are the maximum and minimum values of the i-th indicator eigenvalues in the scenario sets, respectively.
Finally, via constructing the Lagrange function, the cluster iteration can be obtained, and the loop iteration can realize the scenario clustering.

CCQ Evaluation Method
The appropriate number of clustering scenario categories can accurately reflect the characteristics of wind power and photovoltaic output and has a significant impact on the solution effect of the model. Thus, we adopted the CCQ method to determine the best scenario category.
The clustering effect can be expressed by the variance within the cluster. The larger the variance, the worse the clustering effect, and the smaller the variance, the better the clustering effect.
where n is the total number of scenes of set X; x is the sample mean value corresponding to the eigenvalue of an index; and x j is the value of this indicator in the j-th scenario.
The distance between members within each cluster should theoretically be as close as possible, so the smaller the cluster density, the better.
where C is the sum scenarios after clustering; var(c j ) represents the variance in category C j . The clustering proximity and the distance between clusters are inversely proportional. The smaller the clustering proximity, the more effectively each cluster can be separated, indicating that the classification effect is better.
where S C j is the clustering center in category C j ; δ is a gaussian constant and 2δ 2 = 1 for simplified calculation.
In summary, this paper linearly combines clustering density and proximity, that is, the CCQ method, which can be expressed as: where ξ ∈ [0, 1] is the weight of equilibrium clustering density and clustering proximity, which is 0.5 in this paper. The average comprehensive quality evaluation of each index C om determines the optimal number of clustering, as shown in the following formula: The solution procedure of the WT and PV scenario clustering and CHP microgrid optimization scheduling strategy is demonstrated in Figure 4.

Case Study
In this section, the optimization scheduling strategy of CHP microgrid based on MILP is presented. This paper selected a residential area in Beijing China as the experimental system. The structure of the experimental system is illustrated in Figure 1.

Data analysis
In this section, considering the time characteristics of WT and PV output and user load, as well as the fluctuation characteristics of electricity price, the optimal scheduling strategy of CHP microgrid system is presented. This paper selects the historical output power data of WT and PV in the winter in a residential area in Beijing, China, and 10 typical power output curves were obtained using the improved FCM and CCQ models, as shown in Figures 5a and 5b for PV and WT, respectively. The typical daily heat load curve and electric load curve are shown in Figure 6. The hourly price of electricity is observed in Figure 7.

Case Study
In this section, the optimization scheduling strategy of CHP microgrid based on MILP is presented. This paper selected a residential area in Beijing China as the experimental system. The structure of the experimental system is illustrated in Figure 1.

Data Analysis
In this section, considering the time characteristics of WT and PV output and user load, as well as the fluctuation characteristics of electricity price, the optimal scheduling strategy of CHP microgrid system is presented. This paper selects the historical output power data of WT and PV in the winter in a residential area in Beijing, China, and 10 typical power output curves were obtained using the improved FCM and CCQ models, as shown in Figure 5a,b for PV and WT, respectively. The typical daily heat load curve and electric load curve are shown in Figure 6. The hourly price of electricity is observed in Figure 7.  According to Figure 5, the maximum power output period of the PV is 11:00-14:00. In this time period, the user's electricity load and thermal load gradually reduce and reach the minimum load at 14:00. The maximum output power of the WT ranges from 20:00 at night to 06:00 the next day. In this time period, the user's electrical load gradually reduces, and the thermal load begins to reduce after reaching the maximum load at 20:00. Obviously, PV is the perfect complement to WT. Further analysis indicates that the abandonment of wind and solar energy is more prone to take place under the low valley period of users' electricity load. In grid-connected mode, when the user's electricity load demand reaches the peak time, the renewable energy can be fully utilized, and the shortage of power is supplemented by the battery discharge or purchasing from the main grid; when the electricity load is in a low period, the extra power can be stored in the battery or sold to the main grid, which is beneficial to reduce the rejection rate of wind and solar energy. Figure 6 shows a typical user's 24-hour electricity and thermal load in winter. There are two peak periods of energy consumption and two valley periods. Figure 7 presents 24 h purchase and sale electricity prices of the main grid, which is divided into three sections: peak, flat, and valley.  According to Figure 5, the maximum power output period of the PV is 11:00-14:00. In this time period, the user's electricity load and thermal load gradually reduce and reach the minimum load at 14:00. The maximum output power of the WT ranges from 20:00 at night to 06:00 the next day. In this time period, the user's electrical load gradually reduces, and the thermal load begins to reduce after reaching the maximum load at 20:00. Obviously, PV is the perfect complement to WT. Further analysis indicates that the abandonment of wind and solar energy is more prone to take place under the low valley period of users' electricity load. In grid-connected mode, when the user's electricity load demand reaches the peak time, the renewable energy can be fully utilized, and the shortage of power is supplemented by the battery discharge or purchasing from the main grid; when the electricity load is in a low period, the extra power can be stored in the battery or sold to the main grid, which is beneficial to reduce the rejection rate of wind and solar energy. Figure 6 shows a typical user's 24-hour electricity and thermal load in winter. There are two peak periods of energy consumption and two valley periods. Figure 7 presents 24 hours purchase and sale electricity prices of the main grid, which is divided into three sections: peak, flat, and valley.
Finally, the main parameters involved in the CHP microgrid system are shown in Table 2. The relevant parameters of carbon dioxide emissions are demonstrated in Table 3 [28].  According to Figure 5, the maximum power output period of the PV is 11:00-14:00. In this time period, the user's electricity load and thermal load gradually reduce and reach the minimum load at 14:00. The maximum output power of the WT ranges from 20:00 at night to 06:00 the next day. In this time period, the user's electrical load gradually reduces, and the thermal load begins to reduce after reaching the maximum load at 20:00. Obviously, PV is the perfect complement to WT. Further analysis indicates that the abandonment of wind and solar energy is more prone to take place under the low valley period of users' electricity load. In grid-connected mode, when the user's electricity load demand reaches the peak time, the renewable energy can be fully utilized, and the shortage of power is supplemented by the battery discharge or purchasing from the main grid; when the electricity load is in a low period, the extra power can be stored in the battery or sold to the main grid, which is beneficial to reduce the rejection rate of wind and solar energy. Figure 6 shows a typical user's 24-hour electricity and thermal load in winter. There are two peak periods of energy consumption and two valley periods. Figure 7 presents 24 hours purchase and sale electricity prices of the main grid, which is divided into three sections: peak, flat, and valley.
Finally, the main parameters involved in the CHP microgrid system are shown in Table 2. The relevant parameters of carbon dioxide emissions are demonstrated in Table 3 [28]. Finally, the main parameters involved in the CHP microgrid system are shown in Table 2. The relevant parameters of carbon dioxide emissions are demonstrated in Table 3 [27].

Analysis and Discussion of Results
According to the actual operating characteristics of the CHP microgrid system, this paper analyzed the system operating cost, demand response compensation cost, and CO 2 emission cost in five different operation modes based on various conditions affecting the economic operation of the CHP microgrid and then proposed corresponding strategies for system operators. The specific operation mode settings are shown in Table 4.
By comparing the total cost of the five operation modes, integrated demand response compensation costs, and CO 2 emission costs, the optimal economic dispatch strategy is obtained.

Operation Mode 1
In operation mode 1, the electricity load is provided by each microgenerator unit. As shown in Figure 8, only WT is working during the period from 01:00 to 07:00; at the same time, the price of electricity (0.35¥) is lower than the marginal cost (0.605¥) of gas turbine power generation. In order to satisfy the users' load demand, the system operator needs to purchase electricity from the main grid and store a portion of it into the battery to cope with the future increase in electricity load. At 08:00-11:00, the PV unit and the gas turbine start to work, and as the power supply is greater than the demand, electricity price gradually increases; at this time, extra power is sold to the main grid for additional profits. At 12:00-18:00, the internal power supply of the CHP microgrid is balanced. At 19:00-22:00, it is at the peak hours of power consumption; only the power output of the WT and the GT in the system are insufficient to meet the user's electricity load, so in order to avoid purchasing a large amount of electricity from the main grid in the high price period and reduce system operating costs, the system operator needs IEDR to transfer users' load. During the period of 23:00-24:00, the electricity load drops sharply, and the surplus electricity in the system is sold to the main grid for profit. This method effectively improves the flexibility and power supply reliability of the CHP microgrid. Figure 9 clearly illustrates the purchase and sale of electricity as well as charge and discharge behaviors of the CHP microgrid at various times. When the electricity price is low, the electricity is purchased; otherwise, the electricity is sold. In addition, the total amount of electricity sold is more than purchased, which indicates that in the grid-connected mode, the CHP microgrid can promote the consumption of renewable energy, and system operators can also obtain more revenue. Figure 10 shows that the users' thermal load in the CHP microgrid system is always in balance. Changes in users' electricity demand will not affect thermal load, because this mode of operation achieves a partially decoupling heat production from electricity generation. Compared with the traditional "heat-set" and "electricity-set" modes of operation, this mode of CHP microgrid is more flexible and more efficient in energy utilization. Finally, the operation cost of the system is 6188.1433¥, and the cost of CO 2 emission is 879.9429¥. Figure 10 shows that the users' thermal load in the CHP microgrid system is always in balance. Changes in users' electricity demand will not affect thermal load, because this mode of operation achieves a partially decoupling heat production from electricity generation. Compared with the traditional "heat-set" and "electricity-set" modes of operation, this mode of CHP microgrid is more flexible and more efficient in energy utilization. Finally, the operation cost of the system is 6188.1433¥, and the cost of CO2 emission is 879.9429¥.

Operation Mode 2
In operation mode 2, there is no battery device in the system. In order to ensure the balance of electric power in the system, the system operator either sells extra power to the main grid or shuts down the PV and WT power output. In the power shortage period, the users' power demand cannot be satisfied after implementing the demand response; the system operator needs to purchase power from the main network, thereby causing an increase in operating cost. Further, the flexibility of the CHP microgrid is greatly reduced because of the absence of a battery device. When the grid is interrupted, the system fails to respond in time and will cause serious economic losses. The operation cost of the system calculated by simulation is 6191.6159¥, and the cost of CO 2 emission is 875.2604¥.

Operation Mode 2
In operation mode 2, there is no battery device in the system. In order to ensure the balance of electric power in the system, the system operator either sells extra power to the main grid or shuts down the PV and WT power output. In the power shortage period, the users' power demand cannot be satisfied after implementing the demand response; the system operator needs to purchase power from the main network, thereby causing an increase in operating cost. Further, the flexibility of the CHP microgrid is greatly reduced because of the absence of a battery device. When the grid is interrupted, the system fails to respond in time and will cause serious economic losses. The operation cost of the system calculated by simulation is 6191.6159¥, and the cost of CO2 emission is 875.2604¥.

Operation Mode 3
IEDR is an effective solution for power demand side management that can intelligently control various energy loads within the system, increase system flexibility, and maximize energy resilience and savings on operating cost. Implementing demand response to adjust the demand side load, the user's electricity habits can be improved, and the maximum load of the microgrid grid as well as the installed capacity of the system can be reduced. In this mode of operation, the user cannot interact with the CHP microgrid system. In the peak load period, the system operator can only purchase power from the main grid to meet the user load and increase the system purchase cost. Due to the lack of timely response of users to the power grid, the operation security of the CHP microgrid system is not guaranteed, and the reliability of power supply is doubtful. Through modeling and calculation, the operation cost of the system is 6269.3646¥, and the cost of CO 2 emission is 849.3497¥.

Operation Mode 4
The battery and IEDR are considered as an auxiliary function that can improve system flexibility and energy efficiency. Without battery and IEDR, the cost of system operation is higher than that in operation modes 2 and 3. The operation cost of the system calculated by simulation is 6272.7725¥, and the cost of CO 2 emission cost is 843.9117¥.

Operation Mode 5
Under the condition of isolated main grid operation, the internal electricity load of the CHP microgrid system is completely provided by PV, WT, and GT. When the power supply is insufficient, it will cause a large-scale power outage and bring inconvenience to people's daily lives. When the power is surplus, little power is generated by the PV and WT units and massive renewable energy is wasted, which is not conducive to sustainable development. The operation cost of the system calculated by simulation is 6372.2021¥, and the cost of CO 2 emission is 697.7986¥.

Comparative Analysis of Results
The economic operation of the CHP microgrid system is very complicated and is limited by external environmental interference and its own factors. The system operation mode will vary with changes in the combination of various factors. The results of the above five different CHP microgrid operation modes are shown in Table 5.  Table 5 presents the operation cost, the CO 2 emission cost, integrated demand response compensation cost, and the proportion of electricity purchased from the main network in the above five operation modes. The comparison between operation modes 1 and 2 shows that adding a battery device to the CHP microgrid can reduce the total system cost, but the CO 2 emission cost will increase, because the system operator will purchase electricity and store it in the battery at a low electricity price, then sell it at a peak electricity price for additional profit. As a medium for storing electricity, the battery makes it possible for the system operator to purchase more power from the main grid at a lower price, which leads to an increase in the cost of CO 2 emissions. The behavior of earning a differential profit is completely in line with the idea of time-of-use price. In other words, this method is called the arbitrage behavior of the system operator. A comparison between operation modes 1 and 3 shows that the implementation of the integrated demand response program can also motivate operators' arbitrage behavior, which decreases the total system cost and has a greater impact on the operator than that of the battery device. A comparison between operation modes 1 and 5 shows that in the isolated model, the electricity exchange cannot be performed between the CHP microgrid and the main grid, and the internal electricity supply of the CHP microgrid is completely provided by the WT, PV, and GT. When the system load demand exceeds the system supply capacity, in order to maintain the safety of the power system, the system operator must perform load shedding, which may result in damage to the user's interests. When the power supply exceeds the demand, the system operator must reduce the output of the generator unit, which may lead to the abandonment of wind and solar energy. Figure 11a,b shows the relationship between the proportion of electricity purchase and system operation cost and CO 2 emission cost, respectively. It could be clearly judged from the simulation results that there is a close relationship between factors. As the proportion of electricity purchases goes up, the CO 2 emission also gradually increases, while the operating cost of the system gradually decreases and tends to stabilize. The more electricity you buy from the main grid, the more the CO 2 emission. The relationship between these factors is very complicated, involving electricity prices, carbon taxes, etc. Table 6 shows the changes in various costs in the CHP microgrid system under traditional demand response. Comparing operation mode 1 in Tables 3 and 4, it can be concluded that the implementation of the IEDR can effectively reduce the system operation cost and CO 2 emission cost, and the proportion of power purchase is also reduced by about 3%, but the cost of demand response compensation will increase. From operation mode 5, it can be concluded that in the isolated operation mode, the implementation of the IEDR does not change the operation mode of the system, and the system operation cost and CO 2 emission cost will not decrease. By contrast, the demand response compensation cost is greatly reduced. Therefore, the system operator should select the appropriate demand response mode according to the operating state of the microgrid to minimize the system cost.  Table 5 presents the operation cost, the CO2 emission cost, integrated demand response compensation cost, and the proportion of electricity purchased from the main network in the above five operation modes. The comparison between operation modes 1 and 2 shows that adding a battery device to the CHP microgrid can reduce the total system cost, but the CO2 emission cost will increase, because the system operator will purchase electricity and store it in the battery at a low electricity price, then sell it at a peak electricity price for additional profit. As a medium for storing electricity, the battery makes it possible for the system operator to purchase more power from the main grid at a lower price, which leads to an increase in the cost of CO2 emissions. The behavior of earning a differential profit is completely in line with the idea of time-of-use price. In other words, this method is called the arbitrage behavior of the system operator. A comparison between operation modes 1 and 3 shows that the implementation of the integrated demand response program can also motivate operators' arbitrage behavior, which decreases the total system cost and has a greater impact on the operator than that of the battery device. A comparison between operation modes 1 and 5 shows that  CHP microgrid operation is affected by many factors, some of which are variable, such as battery capacity, and some factors which are immutable, such as demand response plan, carbon tax, which are formulated according to policies, and once determined, they will not change for a while. Therefore, this paper studied the impact of battery capacity changes on the economic benefits of the CHP microgrid and obtained the optimal battery capacity under the optimal economic operating condition, which provides a theoretical reference for social capital investment in wind power and the photovoltaic industry.
The energy in the energy storage system can flow in both directions, which can suppress the fluctuation of the load, cut the peak, and fill the valley in the power system, so the power system can maintain system balance at all times. Figure 12 shows that as the battery capacity changes, the operation cost of the CHP microgrid changes accordingly. It can be seen from Figure 12 that as the battery capacity increases, the system cost gradually decreases. When the battery capacity reaches 39 KW, the system cost is minimized, and then the system cost increases rapidly as the battery capacity increases. In the process of battery capacity change, when the battery capacity is less than 39 KW, as the battery capacity increases, the power exchange between the system operator and the main network is promoted, and the profit obtained by the operator is higher than the maintenance cost and installation cost of the battery; when the battery capacity is greater than 39 KW, although the power exchange between the system operator and the main network is promoted, the profit obtained by the operator is less than the maintenance cost and installation cost of the battery. Therefore, it is necessary to comprehensively measure the impact of battery capacity on the economic operation of the microgrid and ameliorate the economic benefits of the CHP microgrid in practice.

Operational Strategy
In this paper, we assumed that the CHP microgrid system is consistent with the interest of users and analyzed the operation mode of the CHP microgrid system from the perspective of the environment and the microgrid operator benefits, respectively. For a microgrid system operator, the ultimate goal is to minimize the system operating cost and maximize the profits, so operation mode 1 is the best strategy choice. When there is a battery and an IEDR is implemented under grid connection, the operator can minimize the system operating cost, because the operator can sell the electricity in peak time stored in the battery and shift the peak power consumption to the low valley period to avoid buying electricity from the main grid when electricity prices are high. Therefore, the It can be seen from Figure 12 that as the battery capacity increases, the system cost gradually decreases. When the battery capacity reaches 39 KW, the system cost is minimized, and then the system cost increases rapidly as the battery capacity increases. In the process of battery capacity change, when the battery capacity is less than 39 KW, as the battery capacity increases, the power exchange between the system operator and the main network is promoted, and the profit obtained by the operator is higher than the maintenance cost and installation cost of the battery; when the battery capacity is greater than 39 KW, although the power exchange between the system operator and the main network is promoted, the profit obtained by the operator is less than the maintenance cost and installation cost of the battery. Therefore, it is necessary to comprehensively measure the impact of battery capacity on the economic operation of the microgrid and ameliorate the economic benefits of the CHP microgrid in practice.

Operational Strategy
In this paper, we assumed that the CHP microgrid system is consistent with the interest of users and analyzed the operation mode of the CHP microgrid system from the perspective of the environment and the microgrid operator benefits, respectively. For a microgrid system operator, the ultimate goal is to minimize the system operating cost and maximize the profits, so operation mode 1 is the best strategy choice. When there is a battery and an IEDR is implemented under grid connection, the operator can minimize the system operating cost, because the operator can sell the electricity in peak time stored in the battery and shift the peak power consumption to the low valley period to avoid buying electricity from the main grid when electricity prices are high. Therefore, the total operating cost of operation mode 1 is lower than that of the other three operating modes, but since the CO 2 emissions per kWh are much greater than the amount generated by the natural gas fuel from Table 1, as a result, the CO 2 emission cost of operation mode 1 is relatively high. From the perspective of environmental benefits, the CO 2 emissions cost of operation mode 5 is the lowest, because under the condition of the isolated model, there is no way for the CHP microgrid system and the main grid to trade electricity with each other, the electricity load all satisfied by the internal generator units of the system. Furthermore, WT and PV will not produce any CO 2 in the CHP microgrid system. GT uses natural gas as a fuel to generate energy, which also produces CO 2 ; despite this, it produces much less CO 2 than that of the main network. Therefore, operation mode 5 should be selected in consideration of environmental benefits. In summary, the economic operation strategy of the CHP microgrid is a complex project, an excellent operation strategy balancing interests and environment, which not only minimizes the operating cost but also considers various factors within the system (e.g., unit power limit, maintenance plan) and external environment (e.g., weather conditions, government policies, operator behavior).

Conclusion
This paper mainly studied the operation mode of the CHP microgrid system with RESs and the behavior of system operators in different scenarios to improve the efficiency of RES use and proposes the operating strategy for system operators under the power market. Considering the battery and IEDR, a MILP model was utilized to optimize the CHP microgrid with RESs. Firstly, considering the uncertainty and volatility of WT and PV output, the FCM and CCQ clustering evaluation method were used to cluster the WT and PV historical power output data into the optimal scenario number, and then the clustering results were comprehensively evaluated to obtain the best typical scenario. Secondly, the MILP method was used to simulate the optimized scheduling of the CHP microgrid in different scenarios, and the operating status of each unit, total system cost, CO 2 emission cost, and integrated demand response compensation cost in the system were obtained. Finally, the behavior of the system operator in different scenarios was analyzed. The results indicate that the battery and IEDR can increase system flexibility and reduce system operating cost. In the grid-connected mode, a shrewd system operator can profit handsomely by making an optimal operation strategy and trading electricity at the time-of-use price. In this regard, more social capital will be introduced to invest in the development of RESs in the future, which is conducive to environmental improvement and sustainable social development. Compared with previous works, the main research in this paper was as follows: (1) an economic operation strategy of the CHP microgrid with RESs considering CO 2 emission, battery, and IEDR was conducted; (2) a novel scene clustering method (FCM-CCQ) was adopted, which can not only cluster but also evaluate whether the scene clustering number is the optimal result; and (3) considering the uncertainty of renewable energy output.
Additionally, we proposed five different operating modes for the CHP microgrid and obtained the results by simulation. After analyzing the results, the strategy choices of the system operators under different circumstances were presented.
Finally, there are some shortcomings in this paper, which are that we did not consider the start-stop cost of generator units, the climbing slope constraint, and the model being idealized. However, it is hoped that the work can provide a useful reference for CHP microgrid economic optimization scheduling containing RESs. In future research, the influence of various factors will be comprehensively considered, making the model construction closer to the actual situation.
This paper can provide useful reference for CHP microgrid economic optimization scheduling containing RESs. As a future work, the influence of various factors will be given more comprehensive consideration, making the model construction more dynamic and accurate. Further, the start-stop cost of generator units and the climbing slope constraint can be employed in this model.