Optimal Planning of Integrated Energy Systems Based on Coupled CCHP

With the widespread attention on clean energy use and energy efficiency, the integrated energy system (IES) has received considerable research and development. This paper proposed an electricity-gas IES optimization planning model based on a coupled combined cooling heating and power system (CCHP). The planning and operation of power lines and gas pipelines are considered. Regarding CCHP as the coupled hub of an electricity-gas system, the proposed model minimizes total cost in IES, with multistage planning and multi-scene analyzing. Renewable energy generation is also considered, including wind power generation and photovoltaic power generation. The numerical results reveal the replacing and adding schemes of power lines and gas pipelines, the optimal location and capacity of CCHP. In comparison with conventional separation production (SP), the optimization model which regards CCHP as the coupled hub attains better economy. At the same time, the influence of electricity price and natural gas price on the quantities of purchasing electricity and purchasing gas in the CCHP system is analyzed. According to the simulation result, a benchmark gas price is proposed, which shows whether the CCHP system chooses power generation. The model results and discussion demonstrate the validity of the model.


Introduction
Owing to low capital costs, high thermal efficiency, and relatively low fuel cost, natural gas-fired generation is more attractive than traditional fossil units.According to the Energy Information Administration [1], natural gas consumed by the electric power sector has shown a sharp growth from 22.3% in 2000 to 34.1% in 2017.At the same time, the interactions between electrical and gas systems have become tighter [2,3].Based on this background, it is necessary to regard the electrical and gas systems as an integrated energy system (IES) [4,5].
Specifically, IES refers to the energy integration system formed by the optimization of energy production, transmission, distribution, conversion, consumption, and storage in the process of planning, design, construction, and operation.IES can bring a great number of benefits, which mainly include complementary strengths of diverse energy systems for system operation and design, facilitating the integration of local renewable and sustainable energy resources, carbon emission reduction by increasing the whole system reliability and energy efficiency, improved system reliability and resilience [6].
Research of IES is presently increasing, and many scholars have explored in the aspect of energy flow, reliability and flexibility evaluation, system planning and operation schedule.Reference [7][8][9][10][11] analyzed energy flow in a coupled electricity and natural gas system.Reference [7] proposed a mixed-integer linear programming (MILP) method to calculate optimal power flow (OPF) and developed a linear energy hub model based on variables to prevent the introduction of dispatch factor variables.Reference [8] proposed a novel decentralized multi-carrier OPF of large-scale IES in a carbon-trading market to give full play to environmental and economic advantages of the system.Reference [9] concentrated on an integrated formulation for the steady-state analysis of electricity-gas coupled systems and proposed a general method to execute a single power and a single gas flow analysis in a uniform framework based on the Newton-Raphson formulation.Reference [10] presented a method for integrated optimization of coupled power flows of various energy infrastructures such as gas, electricity, and district heating systems.Reference [11] proposed a multi-agent genetic algorithm to decompose the multi-energy OPF problem into the separate traditional OPF problem in such a way that most advantages of integrated analysis of IES would be reserved.Reference [12][13][14][15][16][17][18][19] studied the interdependency of IES for system reliability evaluation and flexibility assessment.Reference [12] introduced a new reliability assessment approach to IES, which proposes a hierarchical decoupling optimization framework for both the OPF problems and the energy hub optimal dispatch and accommodates an impact-increment-based state enumeration method to accelerate the reliability assessment process.Reference [13] presented an original methodology to assess the flexibility that the gas system can provide to the power system and introduced a novel metric to evaluate the integrated electrical and gas flexibility, which is used to impose gas-related inter-temporal inter-network constraints on the power system OPF.Reference [14] developed a novel high-resolution spatial and temporal integrated electricity-gas-heat model to evaluate the impact of low-carbon heating options on gas and electricity transmission networks considering heat demand requirements.Reference [15] explored the economic and technical feasibility of deploying the electric boilers and heat storage tanks under practical operational regulations and typical power grids.Reference [16] proposed an integrated gas and electricity transportation system programming algorithm to enhance the power system resilience in extreme events.The proposed variable uncertainty set describes the influence of extreme conditions, and the power system resilience represented by a set of constraints solves the robust model for the integrated planning.Reference [17] presented a tri-level robust optimization-based network model to minimize the total weighted electricity and gas load shedding of the worst case in the integrated electricity and gas distribution systems with respect to hardening uncertain damages and budget limits by disasters.Reference [18] presented a new reliability evaluation approach to integrate system reconfiguration into the reliability evaluation process.Therefore, state evaluation with system reconfiguration can be conducted autonomously.Reference [19] developed a novel quasi-dynamic simulation model to analyze the interdependence between electricity and natural gas transmission networks and the impact on security of supply.
Regarding system planning and operation schedule, some works are found.Traditionally, power system and natural gas system are planned separately; however, with the increasing importance of IES, there is a growing interest in the integrated planning of IES with the coupled natural gas-fired generators [20][21][22][23][24][25][26][27].Reference [20] presented a new planning expansion model in the electricity and natural gas distribution system with many natural gas-fired generators.The model obtains higher economy as compared to traditional separate planning models.Reference [21] presented a comprehensive long-term planning model of natural gas-fired generators, capacitor banks and natural gas distribution pipelines.The outputs of the planning model are the best location and size of the capacitor banks and the natural gas-fired generators.Reference [22] proposed a two-stage stochastic optimization model that provides a general balance between computational tractability and accuracy and pertains uncertainty of natural gas and electricity demands and analyzed the trade-off of building natural gas-related and other facilities.Reference [23] proposed a stochastic integrated day-ahead scheduling model to dispatch hourly load and generation resources and considered a comprehensive framework to address the dispatchability of a set of fuel-constrained natural gas-fired units in the natural gas transportation network.Reference [24] proposed a unified planning and operation optimization methodology to assess flexibility embedded in both investment and operation stages subject to long-term uncertainties.Reference [25] presented a long-term robust optimization planning model for interdependent systems and considered transmission lines, pipelines, generators, gas compressor stations and gas suppliers as investment candidates.Reference [26] presented a long-term, multistage, and multiarea model for the supply/interconnections expansion planning of IES with considering the natural gas value chain and electrical systems value chain.Reference [27] developed a comprehensive and systematic planning model to understand, optimize and develop energy grids to reach higher social welfare with consideration of the expansion of natural gas-fired generators, electricity transmission lines and gas pipelines.
As other coupling technologies, combined heat and power system (CHP) and combined cooling heating and power system (CCHP) are gaining increasing importance and attention in the process of system planning and operation schedule [28][29][30][31][32][33].The CCHP integrates the prime mover, the alternator, the refrigeration equipment, and the heat recovery system to implement the cascade use of energy [34].In a CCHP system, by-product heat is as much as 60-80% of total primary energy and can be recycled for different uses [35].The integrated efficiency of IES with CCHP can reach more than 80%.Therefore, IES with the coupled CCHP is the most potential and the most promising operation mode in IES today.Reference [28] developed a novel approach for optimal economic dispatch scheduling to minimize environmental emissions and maximize economic profit based on integration of CHP systems with traditional thermoelectric units.Reference [29] proposed a coordinated operation strategy considering the gas hydraulic calculation in gas network AC and power flow in the power network and transferred the power fluctuation of renewable energy sources to cooling or heating system and gas distribution network with the coupled external characteristics modeling of CCHP system.Reference [30] proposed an optimal planning of CHPs at arbitrary nodes in IES and employed energy hub approach to calculate the energy consumption of each bus.Reference [31,32] proposed the cooling and electricity coordinated micro-grid real-time dispatching and day-ahead scheduling models considering the performance of ice-storage air-conditioners and the partial load performance of CCHPs.Reference [33] proposed a MILP model to simulate an electricity-gas IES with CCHP and wind power to find an optimal dispatch strategy.And the system security constraints are integrated into the optimal dispatch model.
The previous works about system planning of IES mainly focus on the integrated planning with the coupled natural gas-fired generators.There are only a few papers that discuss the system planning with CCHP.Motived by the issue above, the main contributions of the paper are as follows: (1) This paper proposes a MILP optimization planning model of IES based on the coupled CCHP considering the planning and operation of power lines and gas pipelines.The renewable sources, including wind and solar power, are also considered.
(2) This paper discusses the topology constraints of radical distribution power network and introduces piecewise linear approximation to deal with the nonlinear nature of gas flow.The model uses a three-stage optimization framework and each stage includes four typical days to analyze the economic dispatch strategy and an extreme day to analyze the stability of system operation.
The remainder of this paper is organized as follows.Section 2 introduces the optimization planning model framework and the mathematical model of IES.Section 3 describes the numerical results.The discussion is presented in Section 4, and conclusions in Section 5.

Optimal Model
The MILP optimization planning model of IES based on the coupled CCHP is proposed.Compared with the planning model based on the coupled natural gas-fired generators, the planning model in this paper has higher energy efficiency.
Based on the traditional expansion planning model of multistage distribution network and gas network [36] and typical CCHP system planning model [34], the IES physical programming model including distribution network, gas network and CCHP system is presented as shown in Figure 1.The CCHP physical programming model is presented as shown in Figure 2. The IES optimization programming model including model input and model output is shown in Figure 3.In the IES physical programming model, some nodes set up the CCHP system as the electric and gas coupled hub, and the other nodes are only the electric load node or the natural gas load node.The CCHP physical programming model describes the connection and energy flow relation of gas turbine, auxiliary gas boiler, heat recovery system, absorption chiller and heating coil in the CCHP system, and how to cooperate with electric air conditioner and renewable energy power generation to output cold and heat load.The optimization objective of this paper is to minimize the total investment and operation cost of IES.Transmission lines or pipes that transmit electricity, gas, cold and heat from nodes to customers are not considered.This part of the line or pipeline is common to different planning models.It is located within the user and depends on the design of the user's architecture.It usually does not change with the external planning.In the IES physical programming model, some nodes set up the CCHP system as the electric and gas coupled hub, and the other nodes are only the electric load node or the natural gas load node.The CCHP physical programming model describes the connection and energy flow relation of gas turbine, auxiliary gas boiler, heat recovery system, absorption chiller and heating coil in the CCHP system, and how to cooperate with electric air conditioner and renewable energy power generation to output cold and heat load.The optimization objective of this paper is to minimize the total investment and operation cost of IES.Transmission lines or pipes that transmit electricity, gas, cold and heat from nodes to customers are not considered.This part of the line or pipeline is common to different planning models.It is located within the user and depends on the design of the user's architecture.It usually does not change with the external planning.The IES optimization programming model contains 3 parts: (1) Model input: (a) The IES planning and construction options include: The distribution line expansion and construction options, the gas pipeline expansion and construction options and the CCHP construction options.(b) The load conditions of IES operation include: The electricity, gas, cold and heat load curves of nodes, and the electricity and gas price curves.(c) IES objective function weights include: The unit distance cost of distribution line, gas pipeline, the unit power cost of CCHP, and maintenance cost, the per kilowatt hour cost of power, the per cubic meter cost of natural gas, the unit power cost of electricity and gas load shedding with corresponding conversion coefficient.
(2) The objective function is added in accordance with each weight value.The economy of planning and construction, the economy of electricity and natural gas trading, and the economy of electricity load shedding and gas load shedding are obtained.The decision variables are divided into two categories: (a) Planning decision variables include replacing existing distribution lines ( P,R, , ), adding new distribution lines ( P,A, , ), replacing existing natural gas pipelines ( G,R, , ), adding new natural gas pipelines ( G,A, , ) and constructing CCHP systems ( CCHP, , ).
(b) Operation decision variables include the operating state ( t y ) of the power distribution lines, gas pipelines and CCHP systems in the IES, the power level ( SUB , , h s t Ε ) of electricity load in the IES system and the flow level ( SUB , , h s t Ψ ) of the natural gas load.
(3) Model output: (a) IES planning and construction plan includes distribution lines expansion and construction scheme, gas pipelines expansion and construction scheme, CCHP construction and selection scheme.The IES optimization programming model contains 3 parts: (1) Model input: (a) The IES planning and construction options include: The distribution line expansion and construction options, the gas pipeline expansion and construction options and the CCHP construction options.(b) The load conditions of IES operation include: The electricity, gas, cold and heat load curves of nodes, and the electricity and gas price curves.(c) IES objective function weights include: The unit distance cost of distribution line, gas pipeline, the unit power cost of CCHP, and maintenance cost, the per kilowatt hour cost of power, the per cubic meter cost of natural gas, the unit power cost of electricity and gas load shedding with corresponding conversion coefficient.
(2) The objective function is added in accordance with each weight value.The economy of planning and construction, the economy of electricity and natural gas trading, and the economy of electricity load shedding and gas load shedding are obtained.The decision variables are divided into two categories: (a) Planning decision variables include replacing existing distribution lines (x P,R,J j,t ), adding new distribution lines (x P,A,K k,t ), replacing existing natural gas pipelines (x G,R,J j,t ), adding new natural gas pipelines (x G,A,K k,t ) and constructing CCHP systems (x CCHP,m n,t ).
(b) Operation decision variables include the operating state (y t ) of the power distribution lines, gas pipelines and CCHP systems in the IES, the power level (E SUB h,s,t ) of electricity load in the IES system and the flow level (Ψ SUB h,s,t ) of the natural gas load.
(3) Model output: (a) IES planning and construction plan includes distribution lines expansion and construction scheme, gas pipelines expansion and construction scheme, CCHP construction and selection scheme.(b) IES operation plan includes CCHP operation scheme and IES network topology scheme under different types of daily load scenarios.

Objective
The main objective of the optimization model proposed in this paper is to get the minimum cost of the construction and operation of IES in multistage planning.Therefore, the objective function includes many aspects:

Objective
The main objective of the optimization model proposed in this paper is to get the minimum cost of the construction and operation of IES in multistage planning.Therefore, the objective function includes many aspects: The objective functions obtained from the above four costs are as follows: where, C LP represents the cost of per unit power load shedding.C LG represents the cost of per unit gas load shedding.r P e,h,T and r P s,h,T represent the h hour power load shedding under extreme day scenario and typical day scenario.r G e,h,T and r G s,h,T represent the h hour gas load shedding under extreme day scenario and typical day scenario.In the modeling, a large numerical value of the load shedding coefficient is set.Optimization program to minimize the objective function will avoid the appearance of load shedding during the operation.The appearance of the cost indicates that the construction plan does not meet the load requirements and the planning and construction plan should be readjusted.

Distribution Network Constraints
Reference [37] proposed a distribution network expansion planning model with distributed power supplying and energy storage based on Kirchhoff's law.The problem of this paper is related to the corresponding distribution network planning constraints.Based on the above model, the conventional constraints of distribution network are modified and expanded accordingly.
(1) Power balance and power limits (KCL) In ( 12), f SUB , f CCHP pgu , r P , f rep , f ec , f eh represent the current vector form of substation injection power, CCHP generation power, node power load shedding, node electric load power, electric air conditioner refrigeration power and electric air conditioner heat power.The above variables form an equation through the nodal incidence matrix (S EL , S RL , S AL ) and the branch current vector (f E , f R , f A ).
(2) Nodal voltage calculation and limits (KVL) Z and f are the impedance and current values of the corresponding type distribution lines, respectively.V is the node voltage column vector.The subscript row j represents the j column.M is a large enough positive number.When the distribution line does not run (y = 0), the inequality is relaxed.When the distribution line runs (y = 1), the inequality becomes the equation to calculate the corresponding node voltage.The node voltage is constrained in (17).
(3) Maximum load constraint of distribution line (5) Power load shedding constraint All constraints mentioned above are set up at any planning time t, any typical day scenario s, and any time point h in formula ( 12)-( 22).( 6) Logic constraint of distribution line construction During all the planning periods, only one expansion or new construction is considered for the same line, and the following constraints are obtained.
After the expansion of the replacement line and the construction of the new line, the corresponding lines will be put into use, thus the following constraints are obtained.
After the expansion line is completed, the original lines are no longer used, and the following constraints are obtained.

) Spanning tree constraint
The distribution network in this paper is a radial network, so it is necessary to consider the spanning tree constraint to ensure that the topology of the distribution network is a radial network.The spanning tree has the following 3 kinds of characters.
α P,R j,t = y P,R,0 α P,E e,t = β P,+,E e,t α P,R j,t = β P,+,R j,t ∑ n∈P LN (β P,back,E e,t The 0-1 variable α in Equations ( 28)-( 30) is used to indicate whether the line exists in spanning tree topology.The 0-1 variables β P,+ and β P,− in Equations ( 31)- (33) indicate "parent node" relationship at both ends of the branch.β P,+ represents a positive direction along the branch.Equations ( 34) and (35) indicate the parent node of the load node exists and is unique.There is not parent node in substation node.β P,back represents the back to this node in the branch.Its value is β P,+ or β P,− .
If the network has 1, 2, 3 three nodes, we suppose the branch direction is 1-2-3.The node 3 is the root node, and the node 2 is the parent node of node 1.We suppose the direction of the branch is the positive direction.The branch of 1-2 is the branch 1, and the branch of 2-3 is the branch 2. For node 2, the direction of the back to node 2 is opposite to the positive direction of branch 1 and is the same direction as the positive direction of branch 2. That is, β P,back

Gas Pipe Network Constraint
Gas pipe network constraints mainly consider natural gas node flow balance constraints, gas pipeline flow and pressure relationship constraints, flow constraints of gas pipeline, natural gas station, gas load shedding and construction logic constraints of gas pipeline.
(1) Flow balance constraint of natural gas node In (36), Ψ SUB , Ψ G r , Ψ CCHP pgu , Ψ CCHP boiler , Ψ rep represent the vector form of natural gas station injection flow, node gas load shedding flow, gas flow consumed by CCHP gas turbine, gas flow consumed by CCHP auxiliary gas boiler, node gas load flow.The above variables form an equation through the nodal incidence matrix (S EP , S RP , S AP ) and the branch flow vector (Ψ E , Ψ R , Ψ A ).
(2) Natural gas flow-pressure relation and pressure amplitude constraint The relationship between the flow of natural gas pipeline and the pressure of the node is In (37), Ψ ij represents natural gas flow from node i to node j.C 2 ij represents pipeline constant of i-j gas pipeline.p i represents gas absolute pressure of node i. p j represents gas absolute pressure of node j.P i represents the square of gas absolute pressure of node i. P j represents the square of gas absolute pressure of node j.
Because MILP is adopted in the optimization model, there are only linear constraints in the constraint conditions.For nonlinear constraints in (37), piecewise linearization is required, as shown in Figure 4.If the network has 1, 2, 3 three nodes, we suppose the branch direction is 1-2-3.The node 3 is the root node, and the node 2 is the parent node of node 1.We suppose the direction of the branch is the positive direction.The branch of 1-2 is the branch 1, and the branch of 2-3 is the branch 2. For node 2, the direction of the back to node 2 is opposite to the positive direction of branch 1 and is the same direction as the positive direction of branch 2. That is,

Gas Pipe Network Constraint
Gas pipe network constraints mainly consider natural gas node flow balance constraints, gas pipeline flow and pressure relationship constraints, flow constraints of gas pipeline, natural gas station, gas load shedding and construction logic constraints of gas pipeline.
(1) Flow balance constraint of natural gas node

, h s t h s t h s t h s t h s t h s t h s t h s t
In (36), represent the vector form of natural gas station injection flow, node gas load shedding flow, gas flow consumed by CCHP gas turbine, gas flow consumed by CCHP auxiliary gas boiler, node gas load flow.The above variables form an equation through the nodal incidence matrix ( EP S , RP S , AP S ) and the branch flow vector ( E Ψ , R Ψ , A Ψ ).
(2) Natural gas flow-pressure relation and pressure amplitude constraint The relationship between the flow of natural gas pipeline and the pressure of the node is In (37), ij Ψ represents natural gas flow from node i to node j. 2 ij C represents pipeline constant of i-j gas pipeline.i p represents gas absolute pressure of node i. j p represents gas absolute pressure of node j. i P represents the square of gas absolute pressure of node i. j P represents the square of gas absolute pressure of node j.
Because MILP is adopted in the optimization model, there are only linear constraints in the constraint conditions.For nonlinear constraints in (37), piecewise linearization is required, as shown in Figure 4.In (37), the function on the leftmost side of the equation is a nonlinear function.Suppose f = sign(Ψ ij )Ψ 2 ij and the function curve is divided into a linear combination of N segments.The points of each subsection are Introduce a new variable z k .The piecewise linear function f linear and natural gas pipeline flow Ψ are respectively expressed as Reintroduce a new 0-1 variable y k .
The relationship between z k and y k meets the following constraints.
According to the above steps, the nonlinear constraints can be transformed into linear constraints.In (37), the relationship between natural gas flow and pressure can be transformed into the following constraints.
(3) Maximum flow constraint of gas pipeline (5) flow constraint of natural gas load shedding The extreme load scenario usually refers to the heaviest load period of each planning period and is used to reserve the safe boundary of the system operation.The constraints are the same as the typical day scenario.

CCHP Constraint
The relationship between devices and energy flow in CCHP is shown in Figure 2.
(1) Energy flow relationship of CCHP internal equipment The following formulas are set up under the condition of ∀n ∈ Π CCHP .

Renewable Energy Generation
Renewable energy generation includes wind power generation and photovoltaic power generation.We need to consider the complementarity of wind and light and timing characteristic.This paper studies the 4 timing characteristics of spring, summer, autumn and winter, and models two kinds of power generation output.
The operation of wind power generation and photovoltaic power generation depends on local weather temperature and geographical environment and has high randomness and volatility.The output of wind turbines is mainly related to wind speed.The output of photovoltaic power is mainly related to the temperature and the radiance of the sun, both of which have obvious timing characteristics.Thus, the season and the weather have a great influence on wind power generation and photovoltaic power generation.There is the highest temperature and the strongest solar radiation in summer in China, so the photovoltaic power is the most.However, the wind is the weakest in summer, so the output of the wind turbine is limited.In winter, both temperature and solar radiation are low, so the output of photovoltaic power is the smallest.However, the wind is the strongest, so the output of the wind turbines is the largest.This shows that wind power and PV have very good complementarity.The photovoltaic output is greatly affected by the weather.The radiance will be greatly reduced on cloudy and rainy days, so the photovoltaic output will also be reduced at the same time.The seasonal timing characteristics of wind turbine and PV in a typical area are shown in Figures 5 and 6.

Renewable Energy Generation
Renewable energy generation includes wind power generation and photovoltaic power generation.We need to consider the complementarity of wind and light and timing characteristic.This paper studies the 4 timing characteristics of spring, summer, autumn and winter, and models two kinds of power generation output.
The operation of wind power generation and photovoltaic power generation depends on local weather temperature and geographical environment and has high randomness and volatility.The output of wind turbines is mainly related to wind speed.The output of photovoltaic power is mainly related to the temperature and the radiance of the sun, both of which have obvious timing characteristics.Thus, the season and the weather have a great influence on wind power generation and photovoltaic power generation.There is the highest temperature and the strongest solar radiation in summer in China, so the photovoltaic power is the most.However, the wind is the weakest in summer, so the output of the wind turbine is limited.In winter, both temperature and solar radiation are low, so the output of photovoltaic power is the smallest.However, the wind is the strongest, so the output of the wind turbines is the largest.This shows that wind power and PV have very good complementarity.The photovoltaic output is greatly affected by the weather.The radiance will be greatly reduced on cloudy and rainy days, so the photovoltaic output will also be reduced at the same time.The seasonal timing characteristics of wind turbine and PV in a typical area are shown in Figures 5 and 6    According to the timing characteristic curves of wind turbine and PV, we can see that the two outputs have seasonal complementarity.In the time of the PV output of 0, the wind turbine works well.The PV output is stronger in the daytime, while the wind turbine output is weaker.So, the daily output of the two kinds of generation is complementary.Due to the randomness and volatility of photovoltaic and wind turbines, the impact of this fluctuation on the power grid can be balanced through the CCHP system.
In this paper, fixed capacity wind power generation and photovoltaic power generation are set up at relevant nodes as part of grid input.CCHP system is used to stabilize its volatility.

Model Solution
The optimization model proposed in this paper is a mixed-integer programming model.Mixedinteger programming including MILP and mixed-integer nonlinear programming (MINLP).MINLP is one of the most difficult problems in the field of planning.It belongs to the non-deterministic polynomial (NP)-hard problem.Although many scholars have carried out extensive research on MINLP problem, there is still a big gap in theory and algorithm compared with MILP problem.In addition, there are still some bottlenecks that are difficult to solve.It is difficult to determine whether the current solution is optimal or not.Therefore, this paper uses a more widely applicable MILP model, whose logic is clear, and it is easy to converge.The good scalability is beneficial to the global optimization of IES.
The optimization model is solved by MATLAB + YALMIP + GUROBI, that is, the mathematical model is established by programming tool software YALMIP in MATLAB, and the model is solved by commercial optimization solver GUROBI.The versions of software are MATLAB 2017b, YALMIP 20180612 and GUROBI 7.5.2.The CPU of desktop is Intel Core i5-4440 3.10 GHz with 8.00 GB RAM.The computational time of the easy example is about 60 min, and the relative error between the optimal feasible solution and the optimal solution is 1%.According to the timing characteristic curves of wind turbine and PV, we can see that the two outputs have seasonal complementarity.In the time of the PV output of 0, the wind turbine works well.The PV output is stronger in the daytime, while the wind turbine output is weaker.So, the daily output of the two kinds of generation is complementary.Due to the randomness and volatility of photovoltaic and wind turbines, the impact of this fluctuation on the power grid can be balanced through the CCHP system.
In this paper, fixed capacity wind power generation and photovoltaic power generation are set up at relevant nodes as part of grid input.CCHP system is used to stabilize its volatility.

Model Solution
The optimization model proposed in this paper is a mixed-integer programming model.Mixed-integer programming including MILP and mixed-integer nonlinear programming (MINLP).MINLP is one of the most difficult problems in the field of planning.It belongs to the non-deterministic polynomial (NP)-hard problem.Although many scholars have carried out extensive research on MINLP problem, there is still a big gap in theory and algorithm compared with MILP problem.In addition, there are still some bottlenecks that are difficult to solve.It is difficult to determine whether the current solution is optimal or not.Therefore, this paper uses a more widely applicable MILP model, whose logic is clear, and it is easy to converge.The good scalability is beneficial to the global optimization of IES.
The optimization model is solved by MATLAB + YALMIP + GUROBI, that is, the mathematical model is established by programming tool software YALMIP in MATLAB, and the model is solved by commercial optimization solver GUROBI.The versions of software are MATLAB 2017b, YALMIP 20180612 and GUROBI 7.5.2.The CPU of desktop is Intel Core i5-4440 3.10 GHz with 8.00 GB RAM.The computational time of the easy example is about 60 min, and the relative error between the optimal feasible solution and the optimal solution is 1%.

Case Conditions
In this paper, an IES system is simulated with a network coupled with an 18-node distribution network and a 11-node low pressure natural gas network.The IES network topology is shown in Figure 7.The distribution network has 16 load nodes and 2 substation nodes, 24 planning lines of 20kV.Node 17 and node 18 are substation nodes, and the other nodes are electricity load nodes.
The natural gas network has 9 load nodes and 2 natural gas station nodes, 15 planning pipelines.Node 10 and node 11 nodes are natural gas station nodes, and the other nodes are gas load nodes.
The double solid line shown in the figure represents a fixed line or pipe.The single solid line represents the line or pipe to be replaced.The dashed line represents a new line or pipe to be added.
As shown in Figure 8, each natural gas load node covers natural gas, cold and heat loads in the adjacent electric load node.The natural gas node is set at one point in the coverage area, and the CCHP system is set at the node.
The two cases are listed as follows: Case 1: CCHP is the couple of IES with electric air conditioner cooperated to supply electricity, gas, cold and heat load.
Case 2: SP cooperated with electric air conditioner to supply electricity, gas, cold and heat load.

Case Conditions
In this paper, an IES system is simulated with a network coupled with an 18-node distribution network and a 11-node low pressure natural gas network.The IES network topology is shown in Figure 7.The distribution network has 16 load nodes and 2 substation nodes, 24 planning lines of 20kV.Node 17 and node 18 are substation nodes, and the other nodes are electricity load nodes.
The natural gas network has 9 load nodes and 2 natural gas station nodes, 15 planning pipelines.Node 10 and node 11 nodes are natural gas station nodes, and the other nodes are gas load nodes.
The double solid line shown in the figure represents a fixed line or pipe.The single solid line represents the line or pipe to be replaced.The dashed line represents a new line or pipe to be added.
As shown in Figure 8, each natural gas load node covers natural gas, cold and heat loads in the adjacent electric load node.The natural gas node is set at one point in the coverage area, and the CCHP system is set at the node.
The two cases are listed as follows: Case 1: CCHP is the couple of IES with electric air conditioner cooperated to supply electricity, gas, cold and heat load.
Case 2: SP cooperated with electric air conditioner to supply electricity, gas, cold and heat load.

Case Conditions
In this paper, an IES system is simulated with a network coupled with an 18-node distribution network and a 11-node low pressure natural gas network.The IES network topology is shown in Figure 7.The distribution network has 16 load nodes and 2 substation nodes, 24 planning lines of 20kV.Node 17 and node 18 are substation nodes, and the other nodes are electricity load nodes.
The natural gas network has 9 load nodes and 2 natural gas station nodes, 15 planning pipelines.Node 10 and node 11 nodes are natural gas station nodes, and the other nodes are gas load nodes.
The double solid line shown in the figure represents a fixed line or pipe.The single solid line represents the line or pipe to be replaced.The dashed line represents a new line or pipe to be added.
As shown in Figure 8, each natural gas load node covers natural gas, cold and heat loads in the adjacent electric load node.The natural gas node is set at one point in the coverage area, and the CCHP system is set at the node.
The two cases are listed as follows: Case 1: CCHP is the couple of IES with electric air conditioner cooperated to supply electricity, gas, cold and heat load.
Case 2: SP cooperated with electric air conditioner to supply electricity, gas, cold and heat load.The physical programming model of SP, as shown in Figure 9, shows the connection and energy flow relation of gas boiler and heat exchanger, and how to cooperate with the electric air conditioner to output cold and heat load.
The CCHP and SP options are shown in Table 1.The total construction cost and total operating cost in the CCHP selection scheme include the cost of gas turbine, gas boiler, heat recovery system and other equipment.The total construction cost and total operating cost in the SP selection scheme include the cost of gas-fired boilers, heat exchangers and other equipment.The three different capacity schemes of CCHP system and the three different capacity schemes of SP system are related to the characteristics of gas turbine and the gas boiler and their functions in their respective systems.Different capacity gas turbines and gas-fired boilers correspond to the different capacity schemes of the other equipment of their respective systems.
The physical programming model of SP, as shown in Figure 9, shows the connection and energy flow relation of gas boiler and heat exchanger, and how to cooperate with the electric air conditioner to output cold and heat load.
The CCHP and SP options are shown in Table 1.The total construction cost and total operating cost in the CCHP selection scheme include the cost of gas turbine, gas boiler, heat recovery system and other equipment.The total construction cost and operating cost in the SP selection scheme include the cost of gas-fired boilers, heat exchangers and other equipment.The three different capacity schemes of CCHP system and the three different capacity schemes of SP system are related to the characteristics of gas turbine and the gas boiler and their functions in their respective systems.Different capacity gas turbines and gas-fired boilers correspond to the different capacity schemes of the other equipment of their respective systems.

Case Results
The cost comparison of each case is shown in Table 2.
The time price curve of electricity load is derived from the official website of PJM company of the United States.The time-sharing price curve of natural gas comes from the official website of Henry trading center of the United States.

Case Results
The cost comparison of each case is shown in Table 2.
The time price curve of electricity load is derived from the official website of PJM company of the United States.The time-sharing price curve of natural gas comes from the official website of Henry trading center of the United States.

Result Analysis
(1) Total cost analysis

Case1 < Case2
The total of CCHP system is lower than that of SP system, which mainly comes from the decline in the cost of purchasing electricity from the power grid.
(2) Cost analysis of distribution network With CCHP coupling electric-gas system, gas turbine power generation can supply local electric load, which reduces transmission power of distribution network.Distribution network can choose distribution line model with lower maximum load flow and more flexible network frame construction.The two cases of third stage line planning results are shown in Figure 10.

Result Analysis
(1) Total cost analysis

Case1 < Case2
The total cost of CCHP system is lower than that of SP system, which mainly comes from the decline in the cost of purchasing electricity from the power grid.
(2) Cost analysis of distribution network With CCHP coupling electric-gas system, gas turbine power generation can supply local electric load, which reduces transmission power of distribution network.Distribution network can choose distribution line model with lower maximum load flow and more flexible network frame construction.The two cases of third stage line planning results are shown in Figure 10.(3) Cost analysis of gas pipeline network In addition to the same natural gas load supplied by the CCHP and SP systems, gas turbines and gas-fired boilers in CCHP also consume natural gas and supply electric, cold and heat loads.In SP, only gas-fired boilers consume natural gas and only supply heat load.Therefore, the natural gas flow of the Case1 is larger, and the selected gas pipeline price is also higher.The two cases of third stage pipeline planning results are shown in Figure 11.(3) Cost analysis of gas pipeline network In addition to the same natural gas load supplied by the CCHP and SP systems, gas turbines and gas-fired boilers in CCHP also consume natural gas and supply electric, cold and heat loads.In SP, only gas-fired boilers consume natural gas and only supply heat load.Therefore, the natural gas flow of the Case1 is larger, and the selected gas pipeline price is also higher.The two cases of third stage pipeline planning results are shown in Figure 11.

(4) Cost of CCHP or SP
There are more internal devices in CCHP, and the construction operation cost under the same capacity is higher than that of SP.At the same time, it is considered that the electric air conditioner built by the user side can supply cold and heat load.Some nodes in the SP system do not have to build the gas boiler for heating.The CCHP system has the high efficiency of comprehensive energy use.compared with the electric air conditioner, the cost of cold and heat load in CCHP is low, so most of the nodes will build a combined cooling heating and power system.Thus, the construction and operation cost of the CCHP system will be higher than that of the SP system.
(5) Cost of purchasing electricity and natural gas On the one hand, the gas turbine in CCHP can use waste heat to supply cold and heat load, on the other hand, the gas turbine generation can supply electricity load.Although more natural gas is used, the cost of gas purchase will be improved.The cost of electricity purchase will be greatly reduced, and the total cost will be reduced.
The cost of purchasing electricity from the power grid and the cost of purchasing gas from the natural gas network in stage 3, typical day 2 of the two cases are shown in Figure 12.The total cost of purchasing electricity and purchasing natural gas in stage 3, typical day 2 of the two cases is shown in Figure 13.From 23 h to 8 h, the cost of purchasing electricity and the cost of purchasing gas in Case1 and Case2 are the same, which indicates that the gas turbine and auxiliary boilers in the CCHP system are not working at this time, and the gas boilers in the SP system are not working.It shows that the electricity price is low, and the cold and heat load are provided by the electric air conditioner.From 9 h to 22 h, the cost of purchasing electricity in the SP system is greater than that of CCHP system and the cost of gas in the SP system is less than that of CCHP system.It indicates that the electricity price is higher, and the absorption chiller and the heat recovery device in the CCHP system begin to supply cool and heat load.This paper does not give unnecessary details on the specific replacing and adding schemes of power lines and gas pipelines, the optimal location and capacity of CCHP.
From the above discussion, it can be concluded that the case without coupled CCHP reveals higher economic costs.Although the cost of electricity is lower, the cost of gas is higher.Therefore, the total investment cost is higher.Obviously, the planning with CCHP as coupled hub attains higher comprehensive use rate of energy and promotes energy conservation and environmental protection.
In view of too many parameters in the model, the computation speed is slow.Thus, we consider adopting the Benders decomposition algorithm to simplify the presented MILP into a two-level optimization framework including master and sub problems.(

4) Cost of CCHP or SP
There are more internal devices in CCHP, and the construction operation cost under the same capacity is higher than that of SP.At the same time, it is considered that the electric air conditioner built by the user side can supply cold and heat load.Some nodes in the SP system do not have to build the gas boiler for heating.The CCHP system has the high efficiency of comprehensive energy use.compared with the electric air conditioner, the cost of cold and heat load in CCHP is low, so most of the nodes will build a combined cooling heating and power system.Thus, the construction and operation cost of the CCHP system will be higher than that of the SP system.
(5) Cost of purchasing electricity and natural gas On the one hand, the gas turbine in CCHP can use waste heat to supply cold and heat load, on the other hand, the gas turbine generation can supply electricity load.Although more natural gas is used, the cost of gas purchase will be improved.The cost of electricity purchase will be greatly reduced, and the total cost will be reduced.
The cost of purchasing electricity from the power grid and the cost of purchasing gas from the natural gas network in stage 3, typical day 2 of the two cases are shown in Figure 12.The total cost of purchasing electricity and purchasing natural gas in stage 3, typical day 2 of the two cases is shown in Figure 13.From 23 h to 8 h, the cost of purchasing electricity and the cost of purchasing gas in Case1 and Case2 are the same, which indicates that the gas turbine and auxiliary boilers in the CCHP system are not working at this time, and the gas boilers in the SP system are not working.It shows that the electricity price is low, and the cold and heat load are provided by the electric air conditioner.From 9 h to 22 h, the cost of purchasing electricity in the SP system is greater than that of CCHP system and the cost of gas in the SP system is less than that of CCHP system.It indicates that the electricity price is higher, and the absorption chiller and the heat recovery device in the CCHP system begin to supply cool and heat load.This paper does not give unnecessary details on the specific replacing and adding schemes of power lines and gas pipelines, the optimal location and capacity of CCHP.
From the above discussion, it can be concluded that the case without coupled CCHP reveals higher economic costs.Although the cost of electricity is lower, the cost of gas is higher.Therefore, the total investment cost is higher.Obviously, the planning with CCHP as coupled hub attains higher comprehensive use rate of energy and promotes energy conservation and environmental protection.
In view of too many parameters in the model, the computation speed is slow.Thus, we consider adopting the Benders decomposition algorithm to simplify the presented MILP into a two-level optimization framework including master and sub problems.

Impact Analysis of Price Factors
It is known from the previous section that the cost of purchasing electricity and gas are the most in the total cost of CCHP.Under the condition of constant electric, gas, cold and heat load, the electricity price, and natural gas price will affect the cost of purchasing electricity and the cost of purchasing gas.At the same time, because CCHP is the coupled hub, the relationship between the electricity price and the natural gas price will directly affect the relationship between the quantity of purchasing electricity and the quantity of purchasing gas.
For the specific typical day, the time price curve of electricity price and gas price will affect the quantity of purchasing electricity and gas.For any hour of the typical day, if the electricity price is

Impact Analysis of Price Factors
It is known from the previous section that the cost of purchasing electricity and gas are the most in the total cost of CCHP.Under the condition of constant electric, gas, cold and heat load, the electricity price, and natural gas price will affect the cost of purchasing electricity and the cost of purchasing gas.At the same time, because CCHP is the coupled hub, the relationship between the electricity price and the natural gas price will directly affect the relationship between the quantity of purchasing electricity and the quantity of purchasing gas.
For the specific typical day, the time price curve of electricity price and gas price will affect the quantity of purchasing electricity and gas.For any hour of the typical day, if the electricity price is

Impact Analysis of Price Factors
It is known from the previous section that the cost of purchasing electricity and gas are the most in the total cost of CCHP.Under the condition of constant electric, gas, cold and heat load, the electricity price, and natural gas price will affect the cost of purchasing electricity and the cost of purchasing gas.At the same time, because CCHP is the coupled hub, the relationship between the electricity price and the natural gas price will directly affect the relationship between the quantity of purchasing electricity and the quantity of purchasing gas.
For the specific typical day, the time price curve of electricity price and gas price will affect the quantity of purchasing electricity and gas.For any hour of the typical day, if the electricity price is higher and the gas price is lower, the price of the gas turbine generation will be low.If the price is lower than the purchasing electricity price, the gas turbine will generate.If the electricity price is lower and the gas price is higher, the price of gas turbine generation is higher.If the price is higher than the purchasing electricity price, the gas turbine will not generate.
Different natural gas prices and fixed electricity price in the stage 3, typical day 2 are shown in Figure 14.The electricity price is the electricity price in the upper section case.The gas price 2 is the gas price in the upper section case.The gas price 1, the gas price 3, and the gas price 4 are 0.73 times, 1.2 times and 1.47 times of the gas price 2, respectively.The electricity price of 4 h is the low, and the electricity price of 17 h and 21 h is high.Compared with electricity price, the price of gas fluctuates less in the whole day.
Energies 2018, 11, x FOR PEER REVIEW 20 of 27 higher and the gas price is lower, the price of the gas turbine generation will be low.If the price is lower than the purchasing electricity price, the gas turbine will generate.If the electricity price is lower and the gas price is higher, the price of gas turbine generation is higher.If the price is higher than the purchasing electricity price, the gas turbine will not generate.Different natural gas prices and fixed electricity price in the stage 3, typical day 2 are shown in Figure 14.The electricity price is the electricity price in the upper section case.The gas price 2 is the gas price in the upper section case.The gas price 1, the gas price 3, and the gas price 4 are 0.73 times, 1.2 times and 1.47 times of the gas price 2, respectively.The electricity price of 4 h is the low, and the electricity price of 17 h and 21 h is high.Compared with electricity price, the price of gas fluctuates less in the whole day.The quantity of purchasing electricity for different gas prices in stage 3, typical day 2 is shown in Figure 15.The quantity of purchasing gas for different gas prices in stage 3, typical day 2 is shown in Figure 16.
Different gas prices and benchmark gas price in stage 3, typical day 2 are shown in Figure 17.The benchmark natural gas price is the natural gas price that matches the given electricity price, which is, based on the electricity price at every moment, the price of the critical natural gas for the gas turbine to choose power generation at this time.When the actual gas price is higher than benchmark gas price, the gas turbine does not generate, and when the actual natural gas price is lower than the price, the gas turbine starts to generate.
From the gas price 1 curve and the benchmark natural gas price curve of Figure 17, the gas price 1 is higher than the benchmark natural gas price during 3 h to 5 h, indicating that the actual gas price is higher, and the gas turbine does not generate at this time.Therefore, the quantity of purchasing electricity of gas price 1 in Figure 15 is the largest, and the quantity of purchasing gas of gas price 1 in Figure 16 is the smallest.For other moments, gas price 1 is lower than the benchmark natural gas price, indicating that the actual gas price is lower and gas turbine generates at the time.Therefore, the quantity of purchasing electricity of gas price 1 in Figure 15 reaches the minimum, and the quantity of purchasing gas of gas price 1 in Figure 16 reaches the maximum.In addition to the fixed electricity load supplied by the gas turbine, all the other electricity loads come from the upper power grid.Thus, the quantity of purchasing electricity fluctuates with the electricity load curve.The quantity of purchasing electricity for different gas prices in stage 3, typical day 2 is shown in Figure 15.The quantity of purchasing gas for different gas prices in stage 3, typical day 2 is shown in Figure 16.
Different gas prices and benchmark gas price in stage 3, typical day 2 are shown in Figure 17.The benchmark natural gas price is the natural gas price that matches the given electricity price, which is, based on the electricity price at every moment, the price of the critical natural gas for the gas turbine to choose power generation at this time.When the actual gas price is higher than benchmark gas price, the gas turbine does not generate, and when the actual natural gas price is lower than the price, the gas turbine starts to generate.
From the gas price 1 curve and the benchmark natural gas price curve of Figure 17, the gas price 1 is higher than the benchmark natural gas price during 3 h to 5 h, indicating that the actual gas price is higher, and the gas turbine does not generate at this time.Therefore, the quantity of purchasing electricity of gas price 1 in Figure 15 is the largest, and the quantity of purchasing gas of gas price 1 in Figure 16 is the smallest.For other moments, gas price 1 is lower than the benchmark natural gas price, indicating that the actual gas price is lower and gas turbine generates at the time.Therefore, the quantity of purchasing electricity of gas price 1 in Figure 15 reaches the minimum, and the quantity of purchasing gas of gas price 1 in Figure 16 reaches the maximum.In addition to the fixed electricity load supplied by the gas turbine, all the other electricity loads come from the upper power grid.Thus, the quantity of purchasing electricity fluctuates with the electricity load curve.From the gas price 2 curve and the benchmark natural gas price curve of Figure 17, the gas price 2 is lower than the benchmark natural gas price during 10 h to 22 h, indicating that the actual gas price is lower, and the gas turbine generates at this time.Therefore, the quantity of purchasing electricity of gas price 2 in Figure 15 reaches the minimum, and the quantity of purchasing gas of gas price 2 in Figure 16 reaches the maximum.For other moments, gas price 2 is higher than the benchmark natural gas price, indicating that the actual gas price is higher, and the gas turbine does    From the gas price 2 curve and the benchmark natural gas price curve of Figure 17, the gas price 2 is lower than the benchmark natural gas price during 10 h to 22 h, indicating that the actual gas price is lower, and the gas turbine generates at this time.Therefore, the quantity of purchasing electricity of gas price 2 in Figure 15 reaches the minimum, and the quantity of purchasing gas of gas price 2 in Figure 16 reaches the maximum.For other moments, gas price 2 is higher than the benchmark natural gas price, indicating that the actual gas price is higher, and the gas turbine does From the gas price 2 curve and the benchmark natural gas price curve of Figure 17, the gas price 2 is lower than the benchmark natural gas price during 10 h to 22 h, indicating that the actual gas price is lower, and the gas turbine generates at this time.Therefore, the quantity of purchasing electricity of gas price 2 in Figure 15 reaches the minimum, and the quantity of purchasing gas of gas price 2 in Figure 16 reaches the maximum.For other moments, gas price 2 is higher than the benchmark natural gas price, indicating that the actual gas price is higher, and the gas turbine does not generate at this time.Therefore, the quantity of purchasing electricity of gas price 2 in Figure 15 reaches the maximum, and the quantity of purchasing gas of gas price 2 in Figure 16 reaches the minimum, which only meets the basic gas load.
Energies 2018, 11, x FOR PEER REVIEW 22 of 27 not generate at this time.Therefore, the quantity of purchasing electricity of gas price 2 in Figure 15 reaches the maximum, and the quantity of purchasing gas of gas price 2 in Figure 16 reaches the minimum, which only meets the basic gas load.From the gas price 3 curve and the benchmark natural gas price curve of Figure 17, the gas price 3 is lower than the benchmark natural gas price during 11 h to 18 h and 21 h to 22 h, indicating that the actual gas price is lower, and the gas turbine generates at this time.Therefore, the quantity of purchasing electricity of gas price 3 in Figure 15 reaches the minimum, and the quantity of purchasing gas of gas price 3 in Figure 16 reaches the maximum.However, the quantity of purchasing electricity is not the lowest, and the quantity of purchasing gas is not the highest.Because the benchmark gas price is only a measure of whether the gas turbine is generating.When the gas turbine chooses to generate, the proportion of the actual generation capacity to the rated power generation is related to the actual natural gas price.According to the simulation analysis, when the actual gas price is not higher than the gas price 2, the actual generation capacity reaches the rated power generation.When the actual gas price is higher than the gas price 2, the actual power generation is lower than the rated power generation, and the higher the actual natural gas price is, the lower the proportion of the actual power generation is.For other moments, gas price 3 is higher than the benchmark natural gas price, indicating that the actual gas price is higher, and the gas turbine does not generate at this time.Therefore, the quantity of purchasing electricity of gas price 3 in Figure 15 reaches the maximum, and the quantity of purchasing gas of gas price 3 in Figure 16 reaches the minimum.
From the gas price 4 curve and the benchmark natural gas price curve of Figure 17, the gas price 4 is higher than the benchmark natural gas price at any time, indicating that the gas turbine does not generate.Therefore, the quantity of purchasing electricity of gas price 4 in Figure 15 reaches the maximum, and the quantity of purchasing gas of gas price 4 in Figure 16 reaches the minimum.

Conclusions
In this paper, an optimal planning model of an electricity-gas IES based on coupled CCHP is proposed.The distribution lines expansion and construction scheme, gas pipelines expansion and From the gas price 3 curve and the benchmark natural gas price curve of Figure 17, the gas price 3 is lower than the benchmark natural gas price during 11 h to 18 h and 21 h to 22 h, indicating that the actual gas price is lower, and the gas turbine generates at this time.Therefore, the quantity of purchasing electricity of gas price 3 in Figure 15 reaches the minimum, and the quantity of purchasing gas of gas price 3 in Figure 16 reaches the maximum.However, the quantity of purchasing electricity is not the lowest, and the quantity of purchasing gas is not the highest.Because the benchmark gas price is only a measure of whether the gas turbine is generating.When the gas turbine chooses to generate, the proportion of the actual generation capacity to the rated power generation is related to the actual natural gas price.According to the simulation analysis, when the actual gas price is not higher than the gas price 2, the actual generation capacity reaches the rated power generation.When the actual gas price is higher than the gas price 2, the actual power generation is lower than the rated power generation, and the higher the actual natural gas price is, the lower the proportion of the actual power generation is.For other moments, gas price 3 is higher than the benchmark natural gas price, indicating that the actual gas price is higher, and the gas turbine does not generate at this time.Therefore, the quantity of purchasing electricity of gas price 3 in Figure 15 reaches the maximum, and the quantity of purchasing gas of gas price 3 in Figure 16 reaches the minimum.
From the gas price 4 curve and the benchmark natural gas price curve of Figure 17, the gas price 4 is higher than the benchmark natural gas price at any time, indicating that the gas turbine does not generate.Therefore, the quantity of purchasing electricity of gas price 4 in Figure 15 reaches the maximum, and the quantity of purchasing gas of gas price 4 in Figure 16 reaches the minimum.

Conclusions
In this paper, an optimal planning model of an electricity-gas IES based on coupled CCHP is proposed.The distribution lines expansion and construction scheme, gas pipelines expansion and construction scheme, CCHP construction and selection scheme are obtained through case analysis, which verifies the validity of the model.
The optimization model proposed in this paper is to consider how to plan and run IES rationally from the perspective of regional energy suppliers.In addition, the planning results bring many effective suggestions to regional energy suppliers on system planning and operation schedule.With the continuous deepening of the electric power system reform and the continuous development of IES, the role of regional energy suppliers is bound to appear.How to run IES rationally and efficiently will become a research focus in the future.Node-branch incidence matrix of existing power lines S RL  Node-branch incidence matrix of power lines to be replaced S AL  Node-branch incidence matrix of power lines to be added S EP  Node-branch incidence matrix of existing gas pipelines S RL  Node-branch incidence matrix of gas pipelines to be replaced S AP  Node-branch incidence matrix of gas pipelines to be added Z Impedance of power line C 2  Gas pipeline constant of gas pipeline from node i to node j L cv lower calorific value

Energies 2018 ,
11, x FOR PEER REVIEW 6 of 27 (b) IES operation plan includes CCHP operation scheme and IES network topology scheme under different types of daily load scenarios.

( 1 )
IES construction costs include the distribution line construction cost C LI , the gas pipeline construction cost C PI and the CCHP construction cost C CCHP .(2) IES operation costs include the distribution line operation cost O LI , the gas pipeline operation cost O PI and the CCHP operation cost O CCHP .(3) Energy transaction costs include electricity transaction cost C PT and natural gas transaction cost C GT .(4) Shedding load costs include shedding electricity load cost C CLPL and shedding gas load cost C CLGL .

2 =
For node 3, the direction of the back to node 3 is opposite to the positive direction of branch 2. That is, β P,back

Energies 2018 ,
11, x FOR PEER REVIEW 10 of 27 For node 3, the direction of the back to node 3 is opposite to the positive direction of branch 2. That is,

Figure 4 .
Figure 4. Linearized schematic diagram of relationship between natural gas flow and pressure.

Figure 4 .
Figure 4. Linearized schematic diagram of relationship between natural gas flow and pressure. .

Figure 5 .
Figure 5. Output curve of wind turbine.Figure 5. Output curve of wind turbine.

Figure 5 .
Figure 5. Output curve of wind turbine.Figure 5. Output curve of wind turbine.

Figure 8 .
Figure 8.The corresponding graph of natural gas load node and electric load node.Figure 8.The corresponding graph of natural gas load node and electric load node.

Figure 10 .
Figure 10.Power lines planning topologies of different cases in stage 3.

Figure 10 .
Figure 10.Power lines planning topologies of different cases in stage 3.

Figure 11 .
Figure 11.Gas pipelines planning topologies of different cases in stage 3.

Figure 11 .
Figure 11.Gas pipelines planning topologies of different cases in stage 3.

Figure 12 .
Figure 12.Cost curves of purchasing electricity and purchasing gas for different cases in stage 3, typical day 2.

Figure 13 .
Figure 13.Total cost curves of purchasing electricity and purchasing gas for different cases in stage 3, typical day 2.

Figure 12 .
Figure 12.Cost curves of purchasing electricity and purchasing gas for different cases in stage 3, typical day 2.

Figure 12 .
Figure 12.Cost curves of purchasing electricity and purchasing gas for different cases in stage 3, typical day 2.

Figure 13 .
Figure 13.Total cost curves of purchasing electricity and purchasing gas for different cases in stage 3, typical day 2.

Figure 13 .
Figure 13.Total cost curves of purchasing electricity and purchasing gas for different cases in stage 3, typical day 2.

Figure 14 .
Figure 14.Different gas prices and electricity price in stage 3, typical day 2.

Figure 14 .
Figure 14.Different gas prices and electricity price in stage 3, typical day 2.

Figure 15 .
Figure 15.Quantity of purchasing electricity for different gas prices in stage 3, typical day 2.

Figure 16 .
Figure 16.Quantity of purchasing gas for different gas prices in stage 3, typical day 2.

Figure 15 .
Figure 15.Quantity of purchasing electricity for different gas prices in stage 3, typical day 2.

Figure 15 .
Figure 15.Quantity of purchasing electricity for different gas prices in stage 3, typical day 2.

Figure 16 .
Figure 16.Quantity of purchasing gas for different gas prices in stage 3, typical day 2.

Figure 16 .
Figure 16.Quantity of purchasing gas for different gas prices in stage 3, typical day 2.

Figure 17 .
Figure 17.Different gas prices and benchmark gas price in stage 3, typical day 2.

Figure 17 .
Figure 17.Different gas prices and benchmark gas price in stage 3, typical day 2.
Construction cost of replacing power line j, type J C K RL,k Construction cost of adding power line k, type K O EL,e Operation cost of existing power line e O J RL,j Operation cost of replacing power line j, type J O K AL,k Operation cost of adding power line k, type K C J RP,j Construction cost of replacing gas pipeline j, type J C K RP,k Construction cost of adding gas pipeline k, type K O EP,e Operation cost of existing gas pipeline e O J RP,j Operation cost of replacing gas pipeline j, type J O K AP,k Operation cost of adding gas pipeline k, type K i Discount rate of interest of stage λ Discount rate of interest of year C CCHP P,n,m Construction price of CCHP systems of type m, node n O CCHP P,n,m Operation price of CCHP systems of type m, node n ρ s,T Proportion of s in the planning year T P SUB E,h,s,T Price of purchasing electricity from power grid to power substation P SUB G,h,s,T Price of purchasing gas from natural gas network to gas substation C LP Coefficient of power load shedding cost C LG Coefficient of gas load shedding cost f SUB max Maximum current limits of power substation Ψ SUB max Maximum gas flow limits of gas substation f max Maximum current limits of power lines Ψ max Maximum gas flow limits of gas pipelines S EL

Table 1 .
Candidate options of CCHP and SP.

Table 1 .
Candidate options of CCHP and SP.

Table 2 .
Cost comparison for cases.

Table 2 .
Cost comparison for cases.
Electricity required by electric air conditioner to supply cold load E AC eh Electricity required by electric air conditioner to supply heat load Q b Heat provided by the boiler Q ch Heat required by the absorption chiller to handle the cooling requirement of users Q hc Heat required by heating coil to handle the heating requirement of users