Integrated Energy Micro-Grid Planning Using Electricity , Heating and Cooling Demands

Many research work has demonstrated that taken the Combined Cooling Heating and Power system (CCHP) as the core equipment, the integrated energy system (IES) can bring obvious benefit to energy efficiency, CO2 emission reduction and operation economy in urban areas. Compared with isolated IES, integrated energy micro-grid (IEMG) which is formed by connecting multiple regions IES together, through distribution and thermal network, can further improve the reliability, flexibility, cleanliness and economy of regional energy supply. Based on the existing IES model, this paper describes the basic structure of IEMG and built a IEMG planning model. The planning based on the mixed integer linear programming, and economically construction planning scheme are calculated by using known electricity, heating and cooling loads information and the given multiple equipment selection schemes. At last, the model is validated by a case study. The results show that the application of IEMG can effectively improve the economy of regional energy supply.


Introduction
Increasing pressure on energy resources endowment and environmental problems resulting from the use of the energy internet (EI) are a major focus of energy researchers and practitioners [1].Combined cooling, heating, and power unit (CCHP) technology integrates production of power from electrical and thermal systems, and solves problems caused by their separate decision-making frameworks.With support of a CCHP, an integrated energy systems (IES) can provide multiple energy flows (electricity, heating, steam, cooling, and desalination) by combining different energy production equipment (natural gas, solar, wind, etc.).It has become widely accepted as one of the most efficient examples of integration of multiple energy sources [2].The biggest benefit of IES is that different kinds of energy production systems are no longer planned separately or operated independently.It thus takes the overall process of energy production-from generation and transmission to consumption-into full consideration during the stages of planning, construction, and operation [3,4].Numerous cases testify to the strength of IES (with a CCHP at its core) in urban areas where it improves primary energy use efficiency, CO 2 emission reduction, and the operational economy [5,6].Its position in the energy network has become increasingly important.
IES theories can be summarized as electricity-heat IES and electricity-gas-heat IES [7].The author of [8] expounded the importance of electricity and heating networks combination analysis.The physical Energies 2018, 11, 2810 3 of 20 strategy, the model divides the energy cycle into winter, summer, and an interim period, according to changes in load demands.An extreme load scenario was added to further guarantee the accuracy and reliability of the planning results.Finally, this paper verifies the IEMG planning model with practical examples to prove its significance in guiding the construction and operation of IEMG for regional integrated energy suppliers.

The Integrated Energy Micro-Grid
2.1.Structure of the Integrated Energy System (IES) and the Integrated Energy Micro-Grid (IEMG) IES usually consists of a CCHP, distributed power sources (adjustable and/or non-adjustable), the electrical load, heating load, cooling load, thermal network, and the electrical network.It can also be connected to the external power grid by a transformer substation.The Busbar structure for IES is shown in Figure 1 [27].
Energies 2018, 11, x FOR PEER REVIEW 3 of 19 operational strategy, the model divides the energy cycle into winter, summer, and an interim period, according to changes in load demands.An extreme load scenario was added to further guarantee the accuracy and reliability of the planning results.Finally, this paper verifies the IEMG planning model with practical examples to prove its significance in guiding the construction and operation of IEMG for regional integrated energy suppliers.

Structure of the Integrated Energy System (IES) and the Integrated Energy Micro-Grid (IEMG)
IES usually consists of a CCHP, distributed power sources (adjustable and/or non-adjustable), the electrical load, heating load, cooling load, thermal network, and the electrical network.It can also be connected to the external power grid by a transformer substation.The Busbar structure for IES is shown in Figure 1 [27].In IES, the demands of electricity, heating, and cooling loads can be satisfied simultaneously.The electricity demand is generated by the CCHP, distributed power resources (photovoltaic (PV) cells, for example), and an external grid (when the load demand exceeds the total electricity capacity), and so on.The heating demands are fulfilled by the CCHP and the gas boiler, and the cooling demands are met by both an electric chiller and an absorption chiller.The heat recovery boiler (HRB) acts as a waste heat recovery facility, and can collect the waste heat generated by both the CCHP and the gas boiler, which can significantly improve the amount of heat used in the system.The power grid, the thermal equipment linking to the link energy production equipment, and the different loads operate together to achieve energy circulation through the whole system.It is necessary to emphasize that the interaction between the CCHP and the thermal network is bidirectional, so there is a switch apparatus between them to achieve directional selectivity.The heat networks of different regions can also transfer heat through a switch apparatus between the heat exchanger and the heat load.
A schematic diagram of the IEMG is shown in Figure 2. The IEGM connects several regions' IES (here called the subarea of the IEMG) by a micro-grid, a heat network, and a natural gas pipeline network, to make a scheduling balance within the whole region possible.Thus, the IEMG regards the multiple IES as a controllable whole-they can be safely connected to the low-voltage distribution grid and operate in a flexible manner.Meanwhile, through the coordinated control of the equipment in these regions, the IEMG can provide a more economical, efficient, and reliable supply of energy for different kinds of loads.Furthermore, it can connect to the external power grid and the thermal network, and can thus purchase electricity and heat energy when the overall output in the region is insufficient, or sell surplus energy to external buyers.It can, therefore, significantly improve the In IES, the demands of electricity, heating, and cooling loads can be satisfied simultaneously.The electricity demand is generated by the CCHP, distributed power resources (photovoltaic (PV) cells, for example), and an external grid (when the load demand exceeds the total electricity capacity), and so on.The heating demands are fulfilled by the CCHP and the gas boiler, and the cooling demands are met by both an electric chiller and an absorption chiller.The heat recovery boiler (HRB) acts as a waste heat recovery facility, and can collect the waste heat generated by both the CCHP and the gas boiler, which can significantly improve the amount of heat used in the system.The power grid, the thermal equipment linking to the link energy production equipment, and the different loads operate together to achieve energy circulation through the whole system.It is necessary to emphasize that the interaction between the CCHP and the thermal network is bidirectional, so there is a switch apparatus between them to achieve directional selectivity.The heat networks of different regions can also transfer heat through a switch apparatus between the heat exchanger and the heat load.
A schematic diagram of the IEMG is shown in Figure 2. The IEGM connects several regions' IES (here called the subarea of the IEMG) by a micro-grid, a heat network, and a natural gas pipeline network, to make a scheduling balance within the whole region possible.Thus, the IEMG regards the multiple IES as a controllable whole-they can be safely connected to the low-voltage distribution grid and operate in a flexible manner.Meanwhile, through the coordinated control of the equipment in these regions, the IEMG can provide a more economical, efficient, and reliable supply of energy for different kinds of loads.Furthermore, it can connect to the external power grid and the thermal network, and can thus purchase electricity and heat energy when the overall output in the region is insufficient, or sell surplus energy to external buyers.It can, therefore, significantly improve the economic efficiency for integrated energy suppliers.To summarize, due to the multi-region and multi-energy complementation, the electricity reliability, economical efficiency, and comprehensive utilization rate of energy within the region of the IEMG is effectively improved.

The Principle Elements of the IEMG
The IEMG model in this paper includes three main parts: energy generation, the transformation network, and loads.Generation refers to the electricity, heat, and other energy supply equipment, including the gas combustion engine, CCHP, photovoltaics, substation, and boiler.The network denotes the electrical power grid and the thermal network.The loads are the electrical, heating, and cooling loads.This section models these devices mathematically.

The Generation Equipment
(1) Combined Cooling, Heating and Power System (CCHP) In this paper, the CCHP is built using a constant efficiency model, and the relationship between its thermoelectric power and fuel consumption is established through approximation of a linear function [28,29].
The fuel combustion of the CCHP can be calculated by the following: , , = , , /3.6 Accordingly, its electrical output is: The net calorific value of the waste heat is: The net calorific value of the jacket-cooling water is: The expression of fuel combustion describes the total calorific value , , converted from the combustion inflow , , (unit: m 3 /h) per unit time, where is the average calorific value (which is a constant, 32.967 MJ/m 3 for natural gas).The other three expressions describe the electrical power, and the available calorific value in the waste heat generated by the CCHP.In these equations, , , is the power output of the CCHP; , , and , , are the net calorific values of the waste heat and the jacket-cooling water, respectively, in kW.The parameters and are two known coefficients that are used to fit the performance of the CCHP.In practice, the output of the CCHP is subject to other technical constraints, including the service life of the unit, the maximum and minimum output

The Principle Elements of the IEMG
The IEMG model in this paper includes three main parts: energy generation, the transformation network, and loads.Generation refers to the electricity, heat, and other energy supply equipment, including the gas combustion engine, CCHP, photovoltaics, substation, and boiler.The network denotes the electrical power grid and the thermal network.The loads are the electrical, heating, and cooling loads.This section models these devices mathematically.

The Generation Equipment
(1) Combined Cooling, Heating and Power System (CCHP) In this paper, the CCHP is built using a constant efficiency model, and the relationship between its thermoelectric power and fuel consumption is established through approximation of a linear function [28,29].
The fuel combustion of the CCHP can be calculated by the following: Accordingly, its electrical output is: The net calorific value of the waste heat is: The net calorific value of the jacket-cooling water is: The expression of fuel combustion describes the total calorific value Q CCHP s,d,D converted from the combustion inflow V CCHP s,d,D (unit: m 3 /h) per unit time, where θ NG is the average calorific value (which is a constant, 32.967 MJ/m 3 for natural gas).The other three expressions describe the electrical power, and the available calorific value in the waste heat generated by the CCHP.In these equations, p CCHP s,d,D is the power output of the CCHP; q GAS s,d,D and q WA s,d,D are the net calorific values of the waste heat and the jacket-cooling water, respectively, in kW.The parameters α and β are two known coefficients that are used to fit the performance of the CCHP.In practice, the output of the CCHP is subject to other technical constraints, including the service life of the unit, the maximum and minimum output limitation, the ramp rate of the unit, the minimum continuous running time, and the minimum continuous downtime.These constraints and parameters are selected differently during the different optimizing purposes.The main purpose of this paper is the planning and optimization of the IEMG.Therefore, the specific constraints and parameter selection will be detailed later.
(2) Distributed Generation The distributed generation in this paper is the PV power generation system, and the mathematical model applied is as follows.The output of a PV system is affected by weather, temperature, and solar illumination.If the PV output is p PV s,d,D , it can be modeled as: where ζ t is the actual illumination intensity during the tth hour (kW); ζ t,s is the illumination intensity under standard conditions; A is the total area of the PV panels, which is: A = ∑ M m=1 A m , in which, A m is the area of a single panel; f pv is the power derating factor of the PV system, denoting the ratio of the actual output power to the rated output power, which is used to represent the power loss caused by dirt, rainwater, or snow on the PV panels, and by the aging of the panels (its value here is taken to be 0.9); η is the overall conversion efficiency of the PV panels: in which η m is the conversion efficiency of a single panel (kW); α p is the power temperature coefficient (%/ • C) (which is generally −0.47); T stc is the reference temperature of the PV generation system measured under standard conditions (25 • C here); T s can be calculated by: where T a is the ambient temperature (T a ) and v pv is the wind velocity (m/s).
After building the model of the PV system, it is necessary to describe the relationship between solar radiation and the output of the PV system.The beta-distribution probability function can be used to express the output variation of the PV generation system.According to [30]: where α and β are the shape parameters of the beta-distribution; P M is the total power of the PV array, and R M is the maximum power that the PV array can output.
(3) Gas Boiler (GB) The heat required in the IEMG is mainly produced by two devices, a CCHP and a gas boiler.During production, the thermal energy is mainly provided by the CCHP.Once it cannot meet the heating load, the gas boiler (GB) can convert the chemical energy of the fuel into thermal energy with a high conversion efficiency, to achieve the thermal balance of the system.Assuming that the gas boiler converts the energy of the natural gas into heat at a constant conversion efficiency, the thermal power of the natural gas that is consumed by the gas boiler is: and the heat supply efficiency is: where Q CB h,d,D is the thermal power of the fuel consumed by the boiler (kW); q CB h,d,D is the boiler's heat output (kW); and η GB d,D is the heating efficiency coefficient of the boiler.Similar to the CCHP, some waste heat of the gas boiler can be reused.The formula for the available waste heat efficiency is: where q CBr h,d,D is the total available waste heat of the boiler, and η GBr d,D is the waste heat energy efficiency coefficient of the boiler.
(4) Heat Recovery Boiler The recycling of waste heat is an important means for improving energy efficiency.In IES and IEMGs, a heat recovery boiler (HRB) is used to collect waste heat from the system.The model of the heat recovery boiler is: where q re h,d is the output of the heat recovery boiler (kW), q GAS s,d,D is the waste heat from CCHP (kW); q CBr h,d,D is the waste heat from gas boiler (kW); η re is the thermal efficiency of the equipment.

(5) Chiller
There are two kinds of chillers commonly used in IES and IEMGs.These are the absorption chiller (AC) and the electric chiller (EC).ACs are driven by a thermal medium, such as lithium bromide or ammonia solution; during operation, the working medium vaporization absorbs a lot of heat from the refrigerant water, so as to achieve cooling.The refrigeration principle of the EC involves, first, compressing the gas refrigerant by electricity, then discharging the refrigerant into a condenser.Under set pressure and temperature conditions, the low temperature and low pressure refrigerant cools the air or the condensed water in the condenser to achieve a cooling effect.
The models of the two chillers are as follows.For an AC: where q AC s,d is the cooling output (kW); q AC,in h,d,D is the heat input (kW); and η AC d,D is the refrigeration coefficient, which is the ratio of the heat input to the cooling output, and it is usually used to measure the performance of an AC.For an EC: where q EC s,d is the cooling output (kW); p EC,in h,d,D is the electric power input (kW); and η COP,EC d,D is the refrigeration coefficient of the EC.

Energy Network Model
The mathematical expression for the energy network, which connects different devices and different regions in the IEMG, can be represented by network topology.Its connection mode is described by the incidence matrix of the topology structure.It consists of the combination of several IES, as shown in Figure 1, with the framework of Figure 2 as an example.As in Figure 1, suppose that the equipment and energy network connection (both from the power grid and the thermal network) in the IEMG are nodes.Each pipe or power line serves as a branch, taking the flow direction of the working medium as the branch direction.A basic model of the energy network expressed by the incidence matrix is by the following.Assume that V is the node set of the network.and E is the set of the power lines or pipes (referred to in the following as the set of edges).The energy network can Energies 2018, 11, 2810 7 of 20 be represented by the v by e incidence matrix A = A ij .Thus, each node of the network is a row of the incidence matrix, and each edge is a column [31].The relationship between the nodes and edges can be indicated by the sign of A ij .When A ij = 1, the node v i is linked with edge e j and the direction points away from v i .When A ij = −1, v i is also linked with e j and the direction points towards v i .When A ij = 0, v i and e j are not linked.
Through this approach, the incidence matrix can represent any connection modes of the network system, as with the energy network.However, in order to describe the energy network more accurately, the incidence matrix needs to be further expanded.Based on the form of the incidence matrix A, it is split into the start incidence matrix A 1 and the end point incidence matrix A 2 , to represent the node set of the starting or ending points of the power and pipe lines, respectively.Therefore, A 1 and A 2 are defined as follows: Hence, suppose that the basic loop set of G(V, E) is L containing p elements, and its basic loop matrix is B = [B hk ].Thus, in the matrix B, each element B hk describes the relationship of the loop L h (L h ∈ L, h = 1, 2, . . ., p) with edge k (a branch of the grid or pipeline in the thermal network).When B hk = 1, the loop L h is in the same direction as edge k; when B hk = −1, the loop L h is opposite to the edge k.If B hk = 0, the edge k is not in the loop L h .This method can be used to describe most of the network system.The energy network based on the incidence matrix is summarized by [32], and the matrix can express the energy network as: where A is the incidence matrix of the energy network, B is the basic loop matrix, H is the energy extensive flow matrix, and ∆X is the energy-intensive difference matrix.The equivalent transfer characteristics for incompressible fluids in the energy network is: where H is the flow of energy transferred in the network, H * is the equivalent extensive energy flow in the transfer process, R is the transmission resistance, L and S are the length and cross-sectional area of the transmission line, and K is the transfer coefficient of extensibility.By combining these two equations, the energy transfer characteristic equation set for the energy network can be established in order to describe the energy transfer state at each node.The advantage of the incidence matrix is that it turns the topology relationship, and the structure of the nodes and edges in the network, into variables in a matrix which is convenient for calculations.It is also helpful for making real-time variations of the connection mode during network analysis and, therefore, it simplifies and expands the analysis of the energy network system.In addition, the incidence matrix can be used to calculate the power flow at any position of the energy network, which can improve the accuracy of energy network planning and operation optimization.On the basis of the incidence Energies 2018, 11, 2810 8 of 20 matrix expression for the energy network, models of the electrical network and the thermal network, and their constraints can be determined, as follows.
(1) Electrical Network Model Based on the above, by taking into account Kirchhoff's current law (KCL), Kirchhoff's voltage law (KVL), and electricity flow constraints for network systems, the power grid model of the IEMG containing the distributed generation sources can be derived.
a.The Kirchhoff's current constraint (KCL): This equation reflects the equilibrium relationship between the inlet current and the outlet current at any node in scenario t.In the equation, S EA t is the node-branch incidence matrix of the power grid in scenario t, f EA t is the branch current, G E t is the input power of the node from the power generation equipment, r t is the lost electrical load, and d E t is the electrical load at the node.b.Kirchhoff's voltage constraint (KVL) and the voltage magnitude constraint are: Here, row j is the j-th column, T refers to the transpose, V is the column vector of the node voltage, and Z is the line impedance.The inequality defines the magnitudes of the node voltage in which V min and V max are the maximum and the minimum voltage magnitudes, respectively.
(2) Thermal Network Model Similar to the above, using the energy network incidence matrix and Kirchhoff's laws, a model describing the working principles of the thermal network can be established.Three functions are used to define the transmission flow of the working medium, the relationship of heat and flow, and the change of the transmission pressure.In addition, apart from the relationship between the heat transfer and the mass flow, it is necessary to consider the corresponding heat loss [33].The thermal network model is as follows.
a.The transmission flow constraint: where S H A t is the node-branch incidence matrix of the thermal network in scenario t, q H A t is the energy flow in the branch pipe line in scenario t, G E t is the input power of the node from heat generation equipment, and d H t is the thermal load at the node.
b.The heat-flow constraint: According to the equivalent energy transfer characteristic equation for the thermal network, the relationship between the available heat and the flow is: where P t is the energy (heat) intensity in the pipeline section, k is the specific heat capacity of the working medium, T A1 and T A2 are the feed-water temperature and return-water temperature at the node, respectively.After the relationship between the heat and the flow in the thermal network has been determined, the heat intensity at the node conforms to the following heat balance constraint: where Z H A j,t is the demand for the thermal intensity at node j during period t; q H A j,t is the energy flow at node j; P j,t is the heat column vector in the pipeline section connected to node j.

Energy Balance of the IEMG
From the configuration and structure of the IEMG, as shown in Figure 3, the electrical/ cooling/heating loads should be balanced in any district, or balanced over the whole region, and in all of the different scenarios this principle is basically the same.The energy balance equations in the IEMG are introduced below.a. Balance of the electrical load (in all scenarios): The electrical load balance shall be satisfied in any scenario s (s = c, h, t, e).p SUB s is the power from the external grid, p PV s,d is the PV generation power, p CCHP s,d is the power of the CCHP (kW), l s,d is the pure electrical load, and p EC s,d is the power of the electric chiller (kW).b.Balance of the cooling load (in scenarios of a cooling supply period): The demand for cooling is satisfied by two devices: the electric chiller and the heat adsorption chiller: where q CA h,d is the total load of cooling, q EC h,d is the input power of the electric chiller, and q AC h,d is the input power of the adsorption chiller.The demand for cooling is satisfied by two devices: the electric chiller and the heat adsorption chiller: where , is the total load of cooling, , is the input power of the electric chiller, and , is the input power of the adsorption chiller.c.The balance of heating load (in the scenario of a heating supply period): The heating load is satisfied by the CCHP and the gas boiler: where , is the total heating load, , is the heat supplied by the CCHP, and , is the heat from the gas boiler.

Planning Process and Framework
In this study, the procedure for IEMG planning is summarized as follows: (1) Extract the regional division, loads, and other planning-related data and information, and carry out the overall regional energy supply equipment configuration (for electricity, heating, and cooling).c.The balance of heating load (in the scenario of a heating supply period): The heating load is satisfied by the CCHP and the gas boiler: where q H h,d is the total heating load, q CCHP h,d is the heat supplied by the CCHP, and q GB h,d is the heat from the gas boiler.

Planning Process and Framework
In this study, the procedure for IEMG planning is summarized as follows: (1) Extract the regional division, loads, and other planning-related data and information, and carry out the overall regional energy supply equipment configuration (for electricity, heating, and cooling).( 2) Obtain the overall configuration capacity of the energy supply equipment from step 1, combined with the load characteristics of each region, and do the regional equipment type selection and capacity optimization.(3) According to the equipment selection and capacity optimizing results, deduce the electricity, heating and the cooling load balance operation simulation of each region, and output the results.(4) Deduce the load balance operation simulation on the basis of quarterly and extreme scenarios, and output the results.( 5) Test and determine whether the regional and quarterly simulation results conform to the energy flow and all other constraints.(6) If the result does not satisfy all the constraints, adjust the selection and capacity results until all constraints are met, then output the corresponding configuration.
In the steps above, Step (1) is preparatory work.Its main function is to determine the macroscopic capacity of the whole region on the basis of the known information, in order to narrow the scope of subsequent optimization.Step (2) determines the equipment selection and the installed capacity of each region on the basis of the macro-planning results.Steps (3) and (4) formulate the operation strategy and calculate the system operation cost.This is done through regional and situational operation simulation, using the planning scheme and the scenarios determined in the previous steps.Steps (5) and (6) ensure that the results meet the requirements of the constraints, and improve the accuracy of the optimization.The process is also shown in Figure 3.
Assume that an IEMG satisfies its electrical/heating/cooling loads through PV panels, natural gas (for the CCHP and gas boiler), as well as by purchasing electricity from the external grid.It is then necessary to consider the IEMG plan from the aspect of expanding its original capacity or building a new transformer substation and PV system to supply electricity.A new CCHP construction can meet the electrical/cooling/heating energy demands.Increasing the number of gas boilers compensates for heating between the CCHP heating output and the heating demand.Adding more chilling equipment can satisfy the cooling load.Accordingly, the decision process of the IEMG planning model is given in Figure 4.
Thus, the decision variables in the model can be classified into two types: construction and operation.
The constructional decision variables are mainly binary, where '0' and '1' mean to undo/do a decision, respectively.To be specific, in the type selection option The operation decision variables are continuous and include: the electrical generation output of the CCHP, g CCHP s ; the heating generation output of the CCHP, q CCHP s ; the heating power of the gas boiler, q GB s ; the power of the electricity purchased from the external gird to the substation, g SUB s ; and the power of the cooling generation equipment, q C s .In addition, there are four typical load periods mentioned during the optimization, which measure the economics of the operational strategies.These are the transitional period (s = t), the cooling supply period (s = c), the heating supply period (s = h), and the extreme period (s = e).The extreme period indicates the unusual and sudden situation in which high cooling supplementation is required in summer, and makes sure that the results of the planning are reliable under extreme conditions.
In the steps above, Step ( 1) is preparatory work.Its main function is to determine the macroscopic capacity of the whole region on the basis of the known information, in order to narrow the scope of subsequent optimization.Step (2) determines the equipment selection and the installed capacity of each region on the basis of the macro-planning results.Steps ( 3) and ( 4) formulate the operation strategy and calculate the system operation cost.This is done through regional and situational operation simulation, using the planning scheme and the scenarios determined in the previous steps.Steps (5) and ( 6) ensure that the results meet the requirements of the constraints, and improve the accuracy of the optimization.The process is also shown in Figure 3.
Assume that an IEMG satisfies its electrical/heating/cooling loads through PV panels, natural gas (for the CCHP and gas boiler), as well as by purchasing electricity from the external grid.It is then necessary to consider the IEMG plan from the aspect of expanding its original capacity or building a new transformer substation and PV system to supply electricity.A new CCHP construction can meet the electrical/cooling/heating energy demands.Increasing the number of gas boilers compensates for heating between the CCHP heating output and the heating demand.Adding more chilling equipment can satisfy the cooling load.Accordingly, the decision process of the IEMG planning model is given in Figure 4.The constructional decision variables are mainly binary, where '0' and '1' mean to undo/do a decision, respectively.To be specific, in the type selection option d, , is the decision variable for whether to invest in the new CCHP in district D or not; similarly, , is the decision variable of any other power generator in district D. , is the decision variable of the heating generator, , is the decision variable of the cooling generator, and is the decision variable of a new or expanded transformer substation.
The operation decision variables are continuous and include: the electrical generation output of the CCHP, ; the heating generation output of the CCHP, ; the heating power of the gas boiler, ; the power of the electricity purchased from the external gird to the substation, ; and the power of the cooling generation equipment, .In addition, there are four typical load periods mentioned during the optimization, which measure the economics of the operational strategies.These are the transitional period (s = t), the cooling supply period (s = c), the heating supply period (s = h), and the extreme period (s = e).The extreme period indicates the unusual and sudden situation in which high cooling supplementation is required in summer, and makes sure that the results of the planning are reliable under extreme conditions.

The Objective Function
The overall objective of the planning model is to meet the maximum energy needs of the whole region while minimizing the sum of the construction cost and the operation cost.Hence, the objective function in this model consists of three parts: the planned construction cost, the planned operation cost, and the value of lost loads, minC where C I NV is the planned construction cost, C OPE is the planned operation cost, and C VOLL is the value of lost load.The calculation of each part is as follows.
(1) The planned construction cost: The function of the planning model is to select the economic optimal among several construction schemes, then it should be noticed that there are many different supplement construction portfolios.The equation below gives the cost which is determined for construction.Now, assume that plan D is known to be executed, thus, the planned construction cost C I NV of plan D includes the construction cost of the CCHP, the construction cost of the electrical/cooling/heating generator, and the expansion cost of the transformer substation: of expanding the transformer substation (J represents plan J); and x is a binary decision variable, where the cost is taken into account when the value is 1.
(2) The planned operation cost: Here, equipment maintenance and depreciation costs are put aside.The planned operation cost C OPE includes the cost of the fuel for the CCHP and the gas boiler, and the cost of purchased electricity: where r represents the system run cycle, ∑ r r (1+i) r is the total net cost of the annual operation, i is the discount rate, and ε s is the proportional contribution of scenario s to the entire planning period.For instance, when the planning period is one year (12 months), if the heating supply period contains four months from 15 November to 15 March, the proportion is 4/12 = 0.333; if the cooling supply period contains three months from 15 June to 15 September, the proportion is 3/12 = 0.25, then the transitional period contains the other five months and the proportion is 5/12 = 0.417.Pr GAS and Pr SUB are the prices of natural gas and external electricity, respectively.F f uel s,d is the fuel consumption per unit time in district d, which consists of the fuel consumed by the CCHP and the gas boiler: g SUB s,d is the quantity of the electricity purchased from the external grid by a substation in district d.
(3) The value of the lost load: In this part, C VOLL indicates the compensation cost for unsatisfied electrical/heating/cooling loads, which are not supplied during scenario s in district d.
Here, R d,s is the capacity of the lost loads and P VOLL is the cost coefficient of the lost loads.It should be pointed out that P VOLL is set to a relatively high value, in order to avoid load loss during operation.

Constraint Conditions
In planning the IEMG, the variation of generation equipment parameters in the model should be within a certain range.Their constraint conditions are given in the following.
a.The permeability constraint on the distributed generation (DG): The proportion of the DG to the total installed capacity should be within a certain range:

b. Operational constraints for the CCHP:
The operational constraints for the CCHP include the maximum and the minimum outputs, the ramp rate, and the maximum and the minimum continuous running times.They are given by: where p CCHP min,d,D and p CCHP max,d,D are the lower and the upper limits of the CCHP's output (kW).The number of generation units in each district needs to be limited, due to geographical factors.In this model, we only limit the total number of CCHP units, and allow one per district: is the total number of CCHP generation units in district d.This constraint can be adjusted according to planning requirements.c.Operation constraints on the gas boiler: The output heating power of the gas boiler during operation should be no larger than its rated power: where q GB max,d,D is the rated power.Similarly, the construction constraint for the gas boiler is: This ensures that the boilers in district d are of the same capacity.

d. Constraints on the chillers:
The power of the adsorption chiller and the electrical chiller during operation should be no larger than their rated power: 0 ≤ q AC/EC s,d ≤ q AC/EC max,s,d where q AC/EC max,s,d is the rated power of the chiller (kW).e.The power flow constraint: The power flow of the system should be limited according to the magnitude of the current in the electrical network: where f EA max is the upper limit of the current magnitude and y EA is the conductance value at the corresponding node.
f.The balance constraint for heat loss in the thermal network: If there is too much heat loss in the pipelines of the thermal network, the temperature of the working medium in the pipelines will become lower than the temperature of the working medium in the return-water system.As a result, the thermal network will be ineffective.To ensure the efficiency of the thermal network, we therefore need to ensure that the power (temperature) of the usable heat in the pipelines is higher than a critical value and lower than the maximum power that can be transferred in the pipelines: P min t ≤ P i,t * ≤ P max t where P i,t * is the power of the usable heat in the working medium at node j, P min t is the lower critical value of the power of the usable heat, and P max t is the maximum power of the usable heat.If the working medium flows away from node i then the value of P i,t * is positive, otherwise it is negative.
g.The supply capacity constraint for the transformer substation (in the external grid): In the transformer substation, the total power supply capacity should not be greater than the product of the load and the capacity-load ratio, expressed as: The first inequality describes the relationship between the original power supply capacity p SUB 0 and the expanded capacity ∑ J ∅ SUB x SUB J p SUB J .The sum is the expanded total power supply capacity.To ensure accuracy and effectiveness in planning, the product of the total power supply capacity and capacity-load ratio γ should be valid for the extreme load scenario (s = e).The range of the capacity-load ratio γ in this model is 1.8~2.1.The second inequality means that only one transformer substation expansion plan in the set ∅ SUB will be carried out.

Calculation Method
Generally, dynamic programming (DP) algorithms could be implemented to energy planning optimization [34].Here, we used the mixed integer linear programming (MILP) method to solve the IEMG planning model.The model involves the following decision variables: the output of the PV power system, the input and outputs of the CCHP, the electricity purchased from or sold to the external power grid, the input of the conversion equipment, and the input and output of the gas boiler.The model can be solved using a mature algorithm, or directly by commercial software, such as CPLEX, GUROBI, and LINGO [35].In this study, the model was built by the software MatLab and Yalmip, and solved by the optimization software GUROBI.The process of the optimization algorithm is as in Figure 5.

Case Description
In this work, a new development area of a municipality was taken as a case study.The planning data and the predicted annual saturated electrical/cooling/heating load data were already known, as

Case Description
In this work, a new development area of a municipality was taken as a case study.The planning data and the predicted annual saturated electrical/cooling/heating load data were already known, as shown in Figure 6.

Case Description
In this work, a new development area of a municipality was taken as a case study.The planning data and the predicted annual saturated electrical/cooling/heating load data were already known, as shown in Figure 6.From Figure 5, the data of the loads were classified under four scenarios: cooling period, heating period, transitional period, and extreme cooling period.It should be noted that during the calculation the load on the air conditioner should be subtracted during the cooling period.
Alternative planning options for the energy supply in the region are listed in Table 1.The parameters of the different CCHP units are listed in Table 2.The parameters used for the 9.5 MW CCHP unit are the same as those for the 5 MW unit.

6.
Map showing the energy supply planning information of the case study area.
From Figure 5, the data of the loads were classified under four scenarios: cooling period, heating period, transitional period, and extreme cooling period.It should be noted that during the calculation the load on the air conditioner should be subtracted during the cooling period.
Alternative planning options for the energy supply in the region are listed in Table 1.The parameters of the different CCHP units are listed in Table 2.The parameters used for the 9.5 MW CCHP unit are the same as those for the 5 MW unit.Other operational parameters included the calorific value and price of gas: 32.967 MJ/m 3 and 0.5 USD/m 3 , respectively; the external electricity price in this area was a commercial price, which was 0.16 USD/kWh (the purchase price from the weighted average of the peak-valley electricity prices).The total planning period was 10 years and the annual discount rate i was 5%.To minimize the lost load, the value of P VOLL was set to 150,000 × 10 4 USD/MW.
Four cases were set to make the calculation of the planning more accurate: CASE 1 Electricity is only purchased from the external power grid, without consideration of PV and CCHP; CASE 2 A 4 MW PV system at least is built in each district and electricity can be purchased from the external power grid, without consideration of the CCHP; CASE 3 CCHP construction is considered and electricity can be purchased from the external power grid when the output of the CCHP is insufficient, without consideration of PV; CASE 4 A 4 MW PV generation system at least is built in each district, a CCHP is considered, and electricity can be purchased from the external power grid when the outputs of the CCHP and the PV are insufficient.

Results and Analysis
The IEMG planning result for each region as in Figure 7, and the overall economic results are shown in Table 3.    From Figure 7 and Table 3, the results showed that: in CASE 1, the capacity of the transformer substation is 6 × 50 MVA, the construction cost is 0.194 billion USD, and the total planning cost, including the operation cost, is 3.954 billion USD.In CASE 2, a 4 MW PV generation source is constructed in each of the seven districts, providing 28 MW in total, and the capacity of the substation is 5 × 50 MVA, the construction cost is 0.230 billion USD, and the total planning cost is 3.954 billion USD.In CASE 3, a 9.5 MW CCHP is constructed in each district (66.5 MW in total), the capacity of the substation is 3 × 50 MVA, the construction cost is 0.249 billion USD, and the total planning cost is 3.759 billion USD.In CASE 4, a 4 MW PV source and a 9.5 MW CCHP are constructed in each district, the capacity of the substation is 2 × 50 MVA, the construction cost is 0.285 billion USD, and the total planning cost is 3.762 billion USD.
In CASE 1, the demands for electricity are satisfied by purchasing electricity from the external power grid and heating is from the gas boiler.This scheme involved the least equipment.The system's structure was relatively simple, and the construction cost was, therefore, the lowest, but the energy supply form was simple and the operational cost was relatively high: 3.225 billion and 0.535 billion USD for electricity and heat, respectively, and 3.760 billion USD in total (which is the highest of the four CASES).In CASE 2, due to PV sources in each district (28 MW in total), the construction cost of the whole system was increased by 18.52% compared to CASE 1.However, the operation cost of the external electricity was lower, and the thermal generation's cost remained unchanged, leading to CASE 2 having a 5.31% reduction in operation cost and a 4.13% decrease in total planning costs.In CASE 3, a CCHP system was added to supply electricity and thermal energy.Compared to CASE 1, the construction cost was 28.22% higher, but due to the application of the CCHP the cost of purchasing electricity was 31.71%lower.The cost of the heating supply was 16.70% lower, and thus the total cost decreased by 4.86%.In CASE 4, both PV and CCHP are constructed, which meant that the scheme combines the characteristics of CASES 2 and 3. Therefore, compared to CASE 1, the construction cost was 46.73%, higher but the operation cost decreased by 9.37%.Lower operating costs offset the higher construction costs, making CASE 4 the lowest costing among the four cases.Therefore, in CASES 2, 3, and 4, the total planning-operation costs were decreased by 3.74%, 4.86%, and 6.63% compared to CASE 1, respectively.To sum up, the application of PV and the CCHP in the IEMG is efficient in reducing the total planning-operation cost as a whole.

Conclusions
This paper presents an IEMG planning model with distributed photovoltaic by MILP.First of all, the model determines the capacity construction configuration of the energy production equipment by known electricity, heating and cooling loads.Second, to further improve the feasibility of the planning results, the calculated capacity allocation is put into operational cost analysis of heating, cooling, transitional and extreme load scenarios.The model takes district energy suppliers as the main investors; the optimized capacity configuration can meet the overall energy demand of the region in different scenarios and, at the same time, give the construction and operation cost of different sub-regions.A case study is given to prove the validity of the model.The case study is in a seven sub-district development zone, and four comparison schemes are given: CASE 1 (electricity supplied by an external power grid), CASE 2 (supplied by an external grid and PV), CASE 3 (supplied by an external grid and CCHP), and CASE 4 (supplied by all above equipment).The calculations show that the ranking of the total costs is CASE 1 > CASE 2 > CASE 3 > CASE 4. Compared to CASE 1, the total planning-operation costs in the other three cases are decreased by 3.74%, 4.86%, and 6.63%, respectively, which reflects the fact that the construction of the distributed PV and CCHP generation sources are beneficial for reducing the total planning-operation costs.From the results of the model calculation, the model we have proposed can be seen as a theoretical reference for the planning of multi-district IES (an IEMG in this paper).
The study could be further improved in the following aspects: the electrical network model and the thermal network model are relatively simple, and only the constraints of the power flow in the power grid and the thermal network are discussed, without consideration of the variation of energy quantity flow rate, the variations of temperature and pressure in the thermal network, or time delay in the thermal network's heat transmission.In addition, the planning and construction of the framework of the energy network is not fully investigated.

Figure 1 .
Figure 1.The Busbar structure for the integrated energy systems (IES).

Figure 1 .
Figure 1.The Busbar structure for the integrated energy systems (IES).
transmission pressure constraint:

Energies 2018 ,
11, x FOR PEER REVIEW 9 of 19 b.Balance of the cooling load (in scenarios of a cooling supply period):

Figure 4 .
Figure 4.The main variables and decisions of the IEMG planning model.

Figure 4 .
Figure 4.The main variables and decisions of the IEMG planning model.
d,D where M PV d,D , M CCHP d,D , M GB d,D , M GH d,D , and M GC d,D are the construction costs of PV panels, the CCHP, the gas boiler, the heating equipment, and the cooling equipment, respectively, in region d.M SUB J is the cost Energies 2018, 11, 2810 12 of 20 plan D for region d, p PV,min d,D and p pv,max d,D are the lower and upper limits of the active power of the DG (kW).v PV,min d,D and v PV,max d,D are the lower and upper limits of the reactive power (kW) of DG, respectively.p PV d,D and v PV d,D are the actual active power and reactive power of DG in district d (kW).N PV,min d,D and N PV,max d,D are the lower and upper limits of the number of the DG in the system.NG PV d,D is the number of DGs in district d.

Figure 5 .
Figure 5.The optimization algorithm process of IEMG planning.

Figure 5 .
Figure 5.The optimization algorithm process of IEMG planning.

Figure 6 .
Figure 6.Map showing the energy supply planning information of the case study area.

Energies 2018 ,
11, x FOR PEER REVIEW 16 of 19the substation is 3 × 50 MVA, the construction cost is 0.249 billion USD, and the total planning cost is 3.759 billion USD.In CASE 4, a 4 MW PV source and a 9.5 MW CCHP are constructed in each district, the capacity of the substation is 2 × 50 MVA, the construction cost is 0.285 billion USD, and the total planning cost is 3.762 billion USD.

Figure 7 .
Figure 7.The IEMG planning result of each region in different cases (Unit: MW).

Figure 7 .
Figure 7.The IEMG planning result of each region in different cases (Unit: MW).
Energies 2018, 11, x FOR PEER REVIEW 4 of 19economic efficiency for integrated energy suppliers.To summarize, due to the multi-region and multi-energy complementation, the electricity reliability, economical efficiency, and comprehensive utilization rate of energy within the region of the IEMG is effectively improved.

Table 2 .
Performance parameter of the CCHP.

Table 3 .
Cost comparison of the different case study scenarios.