Coordinated Optimal Operation Method of the Regional Energy Internet

The development of the energy internet has become one of the key ways to solve the energy crisis. This paper studies the system architecture, energy flow characteristics and coordinated optimization method of the regional energy internet. Considering the heat-to-electric ratio of a combined cooling, heating and power unit, energy storage life and real-time electricity price, a double-layer optimal scheduling model is proposed, which includes economic and environmental benefit in the upper layer and energy efficiency in the lower layer. A particle swarm optimizer–individual variation ant colony optimization algorithm is used to solve the computational efficiency and accuracy. Through the calculation and simulation of the simulated system, the energy savings, level of environmental protection and economic optimal dispatching scheme are realized.


Introduction
Energy shortages, environmental pollution and climate change are important factors restricting the sustainable development of the world, and the corresponding energy and environmental problems have become a major strategic issue of great concern at home and abroad.The development of the energy internet (EI) has become one of the key ways to solve these problems [1][2][3][4].Since Jeremy Rifkin put forward the concept, the EI has been widely discussed around the world.The concept highlights that future energy consumption should take full account of electricity, heat, gas and other forms of energy coupling [5].In addition to the coordination and intensive use of multiple types of energy through the information system, another important feature of the EI is that the process of energy production, flow, and consumption becomes considerable and controllable [6,7].
At present, there are some preliminary studies in the field of EI, including the research framework [8], the key technology [9], the planning and operation [10][11][12], prognostic [13][14][15], and network security [16].References 8 and 9 summarize the technical structure and the key technology of the EI, which provides some considerations for related field research; prognostic areas include power prediction, predictive diagnosis, and predictive control.In [13], predictive control is applied to the optimized control of the microgrid to achieve the optimal online scheduling.In [14,15], predictive diagnosis is modeled in uninterruptible power supplies and a transformer fault.However, there is no relevant research on the EI, which is perhaps a hot topic for further research.The network structure of the EI is very complex, and network security is another important research direction.Reference [16] has studied the network architecture, security issues, and hardware implementation of a smart grid, dealing with its security aspects and showing solutions for the realization of a wireless network.
The optimal operation strategy is the key factor to achieve the benefits of the EI.However, there is a big difference between EI scheduling and traditional power system scheduling.First of all, there are strong random factors in the power side of renewable energy output and the user side of the load demand, which increase the uncertainties of the scheduling process.Secondly, the system can provide electricity and heat for the users simultaneously, and there is a coupling relationship between the electricity and heat-how to combine the demand of heat and electric load, and the operation constraint of the energy storage (ES) will be a complex decision problem.Finally, there exists a two-way interaction between the regional EI and the external power system.Under the electricity market environment, the EI decision makers need to adjust the power generation plan according to the incentive of different electric power markets to ensure the maximum operation efficiency of the system.
To date, some scholars have carried out relevant research.In [10], the energy flow and information flow characteristics of the EI are analyzed, and the operation mechanism of a dynamic multi-agent system is established with the aim of minimizing the cost of energy supply.However, there is a lack of discussion on the combined cooling, heating and power (CCHP) unit in the model.In [17], the regional integrated energy system optimization scheduling method is established with the aim of economic and environmental protection, and the load-scene classification is optimized, but the function of the ES is not taken into account.In [18], a three-level cooperative global optimization method is proposed, but it is difficult to ensure energy efficiency and the complete decoupling of heat/electricity by using a CCHP unit with a constant heat-to-electric ratio (HTER).In [19,20], energy efficiency is used as a goal to establish an optimization model.In [19], a demand response (DR) program is considered as one of the most cost-effective energy alternatives, and in [20] the development of a CCHP strategy is established to achieve energy efficiency.
The current research is more one-sided, which is mainly reflected in the following ways: (1) The existing analog system is incomplete; a complete regional EI should include multiple types of load (especially heat and cold loads), multi-power (wind turbine (WT), photovoltaic (PV), CCHP, etc.) and ES; (2) CCHP units use a constant HTER strategy, leading to the power and the load not matching, and affecting the system's comprehensive energy efficiency.The objective function is mostly the combination of economic and environmental benefits, but also ignoring the energy efficiency; (3) when dealing with the ES model, it usually considers its operating constraints, but it cannot guarantee the number and depth of its discharge, thus affecting its life expectancy; (4) due to the complexity and multiple objectives of the model, the conventional optimization algorithm is not ideal in convergence, calculation speed and initial sensitivity.
Based on the characteristics of the EI, this paper studies the system architecture, energy flow characteristics and coordinated optimization method of the regional EI.Compared with previous studies, the main contribution of this paper is as follows: Considering a CCHP-adjustable HTER, the ES life equivalent model and real-time electricity prices, a double-layer optimal scheduling model taking into account the economic benefits, environmental benefits and energy efficiency is proposed, and a particle swarm optimizer-individual variation ant colony optimization (PSO-IVACO) algorithm is used to solve the computational efficiency and accuracy.Through the calculation and simulation of the actual system, the energy savings, level of environmental protection and economic optimal dispatching schemes are realized.

Structure
According to the scale of construction, EI can be divided into region-level, city-level, national-level and global EI.Different types of EI have similarities and differences, of which the regional EI has the most basic system architecture, with the following characteristics [18,21]: Sustainability 2017, 9, 848 3 of 14 (1) The regional EI is an organic whole combined by a power grid, a thermal network, a gas network, a regional distributed energy station and an operation center; (2) There is a variety of energy devices in the region, such as WT, PV, CCHP units, and ES equipment (including electric storage, heat storage and cold storage), to achieve a variety of energy complements; (3) There is a wide variety of loads in addition to electric load, including heat and cold load; (4) Information flow and energy flow are more closely linked and complex.A regional EI system structure is shown in Figure 1.
Sustainability 2017, 9, 848 3 of 14 (1) The regional EI is an organic whole combined by a power grid, a thermal network, a gas network, a regional distributed energy station and an operation center; (2) There is a variety of energy devices in the region, such as WT, PV, CCHP units, and ES equipment (including electric storage, heat storage and cold storage), to achieve a variety of energy complements; (3) There is a wide variety of loads in addition to electric load, including heat and cold load; (4) Information flow and energy flow are more closely linked and complex.A regional EI system structure is shown in Figure 1.In the regional EI, the electricity supply is realized by various types of distributed power supplies, storage equipment, CCHPs and electricity purchases from the outside system.The heat and cold supply is mainly supported by the CCHP, cold/heat storage equipment, or electric heating/cooling equipment.It can be seen that the ES and the CCHP are very important, as they are the key equipment to realize the mutual conversion and storage between different energy sources.Therefore, the key to realizing the optimal operation of EI lies in the development of strategies involving these two devices.The energy input and output of the EI are shown in Figure 2.

Communications network
Sustainability 2017, 9, 848 3 of 14 (1) The regional EI is an organic whole combined by a power grid, a thermal network, a gas network, a regional distributed energy station and an operation center; (2) There is a variety of energy devices in the region, such as WT, PV, CCHP units, and ES equipment (including electric storage, heat storage and cold storage), to achieve a variety of energy complements; (3) There is a wide variety of loads in addition to electric load, including heat and cold load; (4) Information flow and energy flow are more closely linked and complex.A regional EI system structure is shown in Figure 1.In the regional EI, the electricity supply is realized by various types of distributed power supplies, storage equipment, CCHPs and electricity purchases from the outside system.The heat and cold supply is mainly supported by the CCHP, cold/heat storage equipment, or electric heating/cooling equipment.It can be seen that the ES and the CCHP are very important, as they are the key equipment to realize the mutual conversion and storage between different energy sources.Therefore, the key to realizing the optimal operation of EI lies in the development of strategies involving these two devices.In the regional EI, the electricity supply is realized by various types of distributed power supplies, storage equipment, CCHPs and electricity purchases from the outside system.The heat and cold supply is mainly supported by the CCHP, cold/heat storage equipment, or electric heating/cooling equipment.It can be seen that the ES and the CCHP are very important, as they are the key equipment to realize the mutual conversion and storage between different energy sources.Therefore, the key to realizing the optimal operation of EI lies in the development of strategies involving these two devices.

Coordinated Optimal Operation Model
The heat/cold equipment in the EI includes a gas boiler, the CCHP unit, absorption chillers, an electric heating/cooling device, and cold/heat storage.The electric power supply equipment includes WT, PV, the CCHP unit, ES, and purchasing extranet.Firstly, the output model and cost model of different power supply devices should be constructed.Then, the coordinated and optimized operation model is constructed.

Power Output Model and Cost Model
(1) Gas boiler model: The amount of fuel consumed by the gas boiler is related to its thermal power: where Q boi (t) is the fuel consumed in ∆t period, H boi (t) is the thermal power, η boi (t) is the thermal efficiency of the boiler.Gas cost is: where C f uel (t) is the gas cost, K f uel is the unit price of natural gas.
(2) CCHP model: The CCHP unit consists of two parts: the gas generator and the waste heat boiler.The amount of fuel consumed by the electric power generation is: where Q G (t) is the fuel consumed in ∆t period, P G (t) and P G0 (t) are the inputted natural gas power of the generator and the waste heat boiler, respectively.η G (t) and η G0 (t) are the power conversion efficiency of the generator and the waste heat boiler, respectively.Gas cost is: At the same time, the heat recovery system in the waste heat boiler concentrates the heat energy generated by the gas generator, and provides the thermal power as follows: where H rs (t) is the thermal power and η rs (t) is the heat recovery efficiency.Define the HTER of the CCHP unit [22]: As can be seen from Equation ( 6), under normal circumstances, when the HTER is not adjustable, the conversion efficiency cannot be changed, and P G0 (t) = 0, η G0 (t) = 0.The waste heat boiler can control P G0 (t) to change the HTER.In a certain electric power output range, the CCHP's electric efficiency changes within a very small range, so the HTER will greatly affect the thermal efficiency of the system (heat-to-electric ratio adjustable principle).
However, a larger amount of fuel is not necessarily better.We need to increase the amount of fuel consumption when P G0 (t) is increased, thus increasing the cost.When the amount of fuel increases to a certain extent, the economy will be reduced.As such, in different scheduling periods, the greatest need is to find the optimal HTER.
(3) The absorption chiller converts the heat generated by the waste heat boiler and gas boiler into cold energy: where H ac (t) is the cold power of the absorption chiller, H c rs (t) and H c boi (t) are the cold power provided by the waste heat boiler and gas boiler, respectively.η ac is the conversion efficiency.
(4) Electric cooling/heating device model: The cold/heat power provided by the electric cooling/heating machine is related to the power consumption and energy efficiency ratio: where H ec (t) and H eh (t) are the cold and heat power provided by the machine, respectively.P ec (t) and P eh (t) are the corresponding power required; η ec and η eh are the corresponding efficiency.
(5) PV and WT model: WT and PV are a nonlinear DC source, and the output power changes with natural conditions, which shows strong nonlinearity.In order to reach the maximum power point, the maximum power point tracking (MPPT) is required for a WT/PV generation system.The operation cost of a WT/PV can be described as follows: where C WT (t) and C PV (t) are the operation cost of WT and PV, respectively.K WT and K PV are the operation management coefficient, P WT (t) and P PV (t) are the corresponding power at time t.
(6) ES model: The remaining capacity of the ES at time t is related to the remaining capacity at time t-1 and charge/discharge capacity from t-1 to t: where SOC(t) is the remaining capacity at time t; P es (t) is the charge/discharge power in the time period t-1, when the ES is discharging, P es (t) > 0; when the ES is charging, P es (t) < 0, η es is the efficiency.For the ES, its life depends on many factors, such as operating temperature, maximum charging current, and depth of discharge.The operating temperature and the charging current of the ES are usually related to their heat dispersion and control system.This paper focuses on the relationship between charging times and the charge/discharge depth.When the charge/discharge depth is R, the maximum cycle times of charging-discharging N ES can be expressed as [23]: Where α 1 ~α5 the characteristic parameters of the ES, which can be obtained from the service life test data provided by the manufacturer.When finishing a charge/discharge cycle, the ES life-loss accounts for 1/N ES percentage of the total loss of the lifespan of the ES, and the equivalent economic loss cost C 0 is: Sustainability 2017, 9, 848 6 of 14 where C cos t is the initial investment cost.Therefore, the equivalent total economic loss for the ES C ES is: (7) Power exchange model: Under an electricity market environment, the EI will exchange power with an external grid; the model is as follows: where C buy is the purchase cost, C sold is the selling income, and P grid (t) is the exchange power.When purchased, it is positive, and vice versa.e 1 (t) and e 2 (t) are the real-time electricity price of purchasing and selling.

Two-Layer Optimization Objective Function
The two-layer optimization model is a system optimization model with double hierarchical structure.The lower layer optimization problem is based on the scheme given by the upper decision.The optimal value of the lower decision is fed back to the upper layer and the upper layer makes a matched global optimal benefit decision according to the optimal value of the lower layer.In this paper, the upper model describes the EI energy allocation strategy, with the best economic and environmental benefits as the goal; the lower model describes the HTER strategy, with the best energy efficiency as the goal.The objective function is shown in the following Figure 3: Therefore, the equivalent total economic loss for the ES ES C is: (7) Power exchange model: Under an electricity market environment, the EI will exchange power with an external grid; the model is as follows: When purchased, it is positive, and vice versa.e1(t) and e2(t) are the real-time electricity price of purchasing and selling.

Two-Layer Optimization Objective Function
The two-layer optimization model is a system optimization model with double hierarchical structure.The lower layer optimization problem is based on the scheme given by the upper decision.The optimal value of the lower decision is fed back to the upper layer and the upper layer makes a matched global optimal benefit decision according to the optimal value of the lower layer.In this paper, the upper model describes the EI energy allocation strategy, with the best economic and environmental benefits as the goal; the lower model describes the HTER strategy, with the best energy efficiency as the goal.The objective function is shown in the following Figure 3: The upper objective function is a combination of economic and environmental benefits: where gass K is the cost conversion coefficient of polluted gas, () env C t is the pollution discharge costs, and f3 is a combination of economic and environmental benefits.
The optimal energy efficiency described in this paper refers to the ratio of the output energy to the input energy, indicating the utilization of the energy.The lower objective function is:  The upper objective function is a combination of economic and environmental benefits: where K gass is the cost conversion coefficient of polluted gas, C env (t) is the pollution discharge costs, and f 3 is a combination of economic and environmental benefits.The optimal energy efficiency described in this paper refers to the ratio of the output energy to the input energy, indicating the utilization of the energy.The lower objective function is: where ξ is the energy efficiency, L e (t), L h (t), L c (t) are the total output electric, heat and cold power, respectively.P e (t), P h (t), P c (t) are the total input of electric, heat and cold power, respectively.Constraints include: (1) As one of the key concerns regarding the implementation of the proposed method is the grid's capability to handle the operation, the condition of the power system is very important.The operation needs to satisfy the power flow constraints and static stability constraints: (a) Power flow constraints.That is, the power output of the heat/cold/electric power and load demand should match: L e (t) (b) In addition to the power flow constraints, we should also consider the grid to be in a stable operating range, which is mainly reflected in the bus voltage: where ∆u i (t) is the voltage deviation of ith bus at time t (per unit value).
(2) Gas generator operating constraints.The power generation efficiency of the gas generator decreases with the decrease of the power, and the power consumption is reduced below the cutoff power.As such, it is necessary to limit the power generation power: where ϕ cut is the resection factor, P max G is the maximum power of the gas generator.
(3) Gas boiler output constraints: where H max boi is the maximum power of the gas boiler, H h boi (t) is the thermal power output by the gas boiler.(4) Heat recovery system constraints: where H max rs is the maximum power of the recovery system, H h rs (t) is the thermal power output by the recovery system.
(5) Absorption chillers constraints: where H max ac is the maximum power of the absorption chillers.(6) Electric cooling/heating constraints: where H max ec and H max eh are the maximum power of the electric cooling/heating devices, respectively.(7) ES constraints: The ES working cycle is one day in duration; the energy stored on that day will be released on the same day.This is because the loss of heat and cold would be huge if long-term storage is used.Notwithstanding, the main goal of the ES is to achieve peak load in one day.So there are: where t 1 is the starting scheduled time, A is the current capacity.At the same time, as with all types of ES equipment, the charge and discharge power should meet the maximum constraints.Take heat storage as an example, when discharge, where S max is the maximum power of discharge/charge, µ is the maximum charge/discharge efficiency.Cold storage and electric storage are similar.At the same time, SOC(t) should meet the constraints: where SOC min and SOC max are the minimum constraint and maximum constraint of SOC, respectively.(8) Tie line power constraint: where P grid,max is the transmission power limit.

Solving Algorithm
The current optimization algorithms include particle swarm optimization (PSO), ant colony optimization (ACO), artificial neural network (ANN) and genetic algorithm (GA).However, these algorithms have some flaws.For example, the ANN algorithm's generalization ability is not strong, sometimes creating overlearning; the GA is easily premature, or it is over calculated.Referring to the method in [24], we used the PSO-IVACO method to solve computational efficiency and accuracy.
Ant colony optimization (ACO) has the characteristics of self-organization and versatility, but it has the disadvantages of poor global search ability and slow searching speed.The IVACO algorithm introduces the individual population difference based on the ACO algorithm, which enhances the global search ability of the ACO algorithm.As a result of the multi-selection strategy, the individual ant mixes the cooperative path selection to avoid the premature convergence of the solution.
Due to the lack of pheromone in the initial search, the accumulation of pheromone is larger, and the search speed is slower in the IVACO algorithm.The PSO algorithm is a kind of stochastic global optimization technique; the required parameters are few, but there are early and local convergence problems.In this paper, the optimal solution of the PSO algorithm is taken as the initial solution of the IVACO algorithm and searches further for the optimal solution.The PSO-IVACO algorithm flow-chart is shown in Figure 4.
Sustainability 2017, 9, 848 8 of 14 where ,max grid P is the transmission power limit.

Solving Algorithm
The current optimization algorithms include particle swarm optimization (PSO), ant colony optimization (ACO), artificial neural network (ANN) and genetic algorithm (GA).However, these algorithms have some flaws.For example, the ANN algorithm's generalization ability is not strong, sometimes creating overlearning; the GA is easily premature, or it is over calculated.Referring to the method in [24], we used the PSO-IVACO method to solve computational efficiency and accuracy.
Ant colony optimization (ACO) has the characteristics of self-organization and versatility, but it has the disadvantages of poor global search ability and slow searching speed.The IVACO algorithm introduces the individual population difference based on the ACO algorithm, which enhances the global search ability of the ACO algorithm.As a result of the multi-selection strategy, the individual ant mixes the cooperative path selection to avoid the premature convergence of the solution.
Due to the lack of pheromone in the initial search, the accumulation of pheromone is larger, and the search speed is slower in the IVACO algorithm.The PSO algorithm is a kind of stochastic global optimization technique; the required parameters are few, but there are early and local convergence problems.In this paper, the optimal solution of the PSO algorithm is taken as the initial solution of the IVACO algorithm and searches further for the optimal solution.The PSO-IVACO algorithm flowchart is shown in Figure 4.
Start PSO parameter initialization.Including the search space, the initial position of the particle, the initial velocity of the particle, the number of iterations Solve the individual optimal position and population optimal position under the current iteration Global optimal value changes within control error?

Update the particle velocity and position
The   The solution of the model is as follows: (1) Input power matrix of load, PV, WT, and other constants.The results of the optimal scheduling at t-1 are taken as the initial value variables at t (if t = 0, then the initial value takes the last time in the previous day); (2) PSO is used to solve iteratively.Calculate the lower target first, obtain the position and velocity of the particles, and feed back to the upper layer.Then, calculate the upper target, update the velocity and position, and feed back to the lower layer.After each iteration, perform a power flow calculation to check whether the stability constraint is met.If it is not satisfied, return to the first step, adjust the initial value.If satisfied, carry on to the next step; (3) Repeat iterations until the objective function value changes within the control range; (4) The optimized variable results from (2) are taken as the initial value of the IVACO iteration.
Calculate by following the steps above; then, output the final result.

Basic Data
The example system in this paper is based on the regional EI shown in Figure 1.We assume that the similar type of devices contained in the EI are same model.The simulation parameters used in this paper are shown in Table 1 [21][22][23][24], and the real-time electricity prices are shown in Appendix A. Typical load scenarios can be divided into the heating, cooling or transition period, in which the cooling period is the most typical.At this time, the electric (E) heating (H)/cooling (C) loads both exist, and the load level is higher than the other period of this year.We selected typical days in the cooling period as an example to optimize the calculation.WT, PV and load are shown in Figure 5.Typical load scenarios can be divided into the heating, cooling or transition period, in which the cooling period is the most typical.At this time, the electric (E) heating (H)/cooling (C) loads both exist, and the load level is higher than the other period of this year.We selected typical days in the cooling period as an example to optimize the calculation.WT, PV and load are shown in Figure 5.As can be seen from Figure 5, the load of the EI exhibits obvious peak-valley characteristics (E load and C load).The peak appears at 14:00 and 19:00, which corresponds to real-time electricity prices.This feature will also profoundly affect the ES strategy.

Tie line power
Set the following three simulation modes: 1. Mode 1 (M1): In this situation, the CCHP unit adopts a constant HTER to supply the electric energy and heat/cold energy; 2. Mode 2 (M2): The model proposed in this paper.The CCHP unit adopts an adjustable HTER.3. Mode 3 (M3): On the basis of M2, the ES operating cost is not included in the objective function, only taking into account its operational constraints.

Simulation Results
After calculating using the PSO-IVACO method, the optimal HTER is shown in Figure 6; the cost calculation results are shown in Figure 7; the results of energy efficiency are shown in Figure 8; and the results of the operation in one day are shown in Table 2.As can be seen from Figure 5, the load of the EI exhibits obvious peak-valley characteristics (E load and C load).The peak appears at 14:00 and 19:00, which corresponds to real-time electricity prices.This feature will also profoundly affect the ES strategy.
Set the following three simulation modes: 1. Mode 1 (M1): In this situation, the CCHP unit adopts a constant HTER to supply the electric energy and heat/cold energy; 2.
Mode 2 (M2): The model proposed in this paper.The CCHP unit adopts an adjustable HTER.

3.
Mode 3 (M3): On the basis of M2, the ES operating cost is not included in the objective function, only taking into account its operational constraints.

Simulation Results
After calculating using the PSO-IVACO method, the optimal HTER is shown in Figure 6; the cost calculation results are shown in Figure 7; the results of energy efficiency are shown in Figure 8; and the results of the operation in one day are shown in Table 2.
Typical load scenarios can be divided into the heating, cooling or transition period, in which the cooling period is the most typical.At this time, the electric (E) heating (H)/cooling (C) loads both exist, and the load level is higher than the other period of this year.We selected typical days in the cooling period as an example to optimize the calculation.WT, PV and load are shown in Figure 5.As can be seen from Figure 5, the load of the EI exhibits obvious peak-valley characteristics (E load and C load).The peak appears at 14:00 and 19:00, which corresponds to real-time electricity prices.This feature will also profoundly affect the ES strategy.
Set the following three simulation modes: 1. Mode 1 (M1): In this situation, the CCHP unit adopts a constant HTER to supply the electric energy and heat/cold energy; 2. Mode 2 (M2): The model proposed in this paper.The CCHP unit adopts an adjustable HTER.3. Mode 3 (M3): On the basis of M2, the ES operating cost is not included in the objective function, only taking into account its operational constraints.

Simulation Results
After calculating using the PSO-IVACO method, the optimal HTER is shown in Figure 6; the cost calculation results are shown in Figure 7; the results of energy efficiency are shown in Figure 8; and the results of the operation in one day are shown in Table 2.     Firstly, the influence of HTER on the optimal operation is analyzed: (1) It can be seen from Figure 6 that the optimal HTER at each time is affected by the change of the heat-electric load ratio in the EI.This shows that the optimal HTER is related to the load distribution in the region.(2) From Figures 7 and 8, and Table 2, we can determine that in the case of adjustable HTER, the CCHP unit enhances the heat power output limit, which leads to an increase in the output of the CCHP unit in the source; its operation cost and energy efficiency have improved compared with the direct use of boiler heating.(3) From Figures 6~8, it can be found that there is a certain positive correlation between the optimal HTER and the energy efficiency, and there is a certain inverse correlation with the operation cost.When the HTER is relatively high, the energy efficiency is high and the cost is also lower.It can be concluded that by further analysis, the HTER and efficiency exist in two valley areas, namely 13:00-14:00 and 18:00-19:00.These two periods are the electric peak periods; the electric power output is high, affecting the CCHP unit thermoelectric balance ratio.(4) The load situation has a greater impact on the operation of the EI.When the load is large, the economy and energy efficiency of the EI will be affected.Therefore, it is effective to improve the comprehensive benefit of the EI by a demand response strategy to reduce the peak and valley difference of the load.Firstly, the influence of HTER on the optimal operation is analyzed: (1) It can be seen from Figure 6 that the optimal HTER at each time is affected by the change of the heat-electric load ratio in the EI.This shows that the optimal HTER is related to the load distribution in the region.(2) From Figures 7 and 8, and Table 2, we can determine that in the case of adjustable HTER, the CCHP unit enhances the heat power output limit, which leads to an increase in the output of the CCHP unit in the source; its operation cost and energy efficiency have improved compared with the direct use of boiler heating.(3) From Figures 6~8, it can be found that there is a certain positive correlation between the optimal HTER and the energy efficiency, and there is a certain inverse correlation with the operation cost.When the HTER is relatively high, the energy efficiency is high and the cost is also lower.It can be concluded that by further analysis, the HTER and efficiency exist in two valley areas, namely 13:00-14:00 and 18:00-19:00.These two periods are the electric peak periods; the electric power output is high, affecting the CCHP unit thermoelectric balance ratio.(4) The load situation has a greater impact on the operation of the EI.When the load is large, the economy and energy efficiency of the EI will be affected.Therefore, it is effective to improve the comprehensive benefit of the EI by a demand response strategy to reduce the peak and valley difference of the load.Firstly, the influence of HTER on the optimal operation is analyzed: (1) It can be seen from Figure 6 that the optimal HTER at each time is affected by the change of the heat-electric load ratio in the EI.This shows that the optimal HTER is related to the load distribution in the region.(2) From Figures 7 and 8, and Table 2, we can determine that in the case of adjustable HTER, the CCHP unit enhances the heat power output limit, which leads to an increase in the output of the CCHP unit in the source; its operation cost and energy efficiency have improved compared with the direct use of boiler heating.(3) From Figures 6-8, it can be found that there is a certain positive correlation between the optimal HTER and the energy efficiency, and there is a certain inverse correlation with the operation cost.When the HTER is relatively high, the energy efficiency is high and the cost is also lower.It can be concluded that by further analysis, the HTER and efficiency exist in two valley areas, namely 13:00-14:00 and 18:00-19:00.These two periods are the electric peak periods; the electric power output is high, affecting the CCHP unit thermoelectric balance ratio.(4) The load situation has a greater impact on the operation of the EI.When the load is large, the economy and energy efficiency of the EI will be affected.Therefore, it is effective to improve the comprehensive benefit of the EI by a demand response strategy to reduce the peak and valley difference of the load.
(5) From some of the calculation results (e.g., 09:00 and 14:00), energy efficiency and economic cost are not achieving the best results simultaneously.This is due to the two-layer objective function of this paper.However, from a longer time scale, the model proposed in this paper strikes a reasonable balance between economic cost and energy efficiency, so that the calculation result of the whole day is satisfactory.
Next, we analyzed the impact of ES on scheduling.It can be seen from the curves that, due to the different ES scheduling strategies of M2 and M3, the optimal HTER, operating cost and energy efficiency are different, and the effect of energy efficiency is especially obvious.In order to better explain this problem, we took electricity storage as an example, and drew output curves within a day of M2 and M3, as shown in Figure 9.
(5) From some of the calculation results (e.g., 09:00 and 14:00), energy efficiency and economic cost are not achieving the best results simultaneously.This is due to the two-layer objective function of this paper.However, from a longer time scale, the model proposed in this paper strikes a reasonable balance between economic cost and energy efficiency, so that the calculation result of the whole day is satisfactory.
Next, we analyzed the impact of ES on scheduling.It can be seen from the curves that, due to the different ES scheduling strategies of M2 and M3, the optimal HTER, operating cost and energy efficiency are different, and the effect of energy efficiency is especially obvious.In order to better explain this problem, we took electricity storage as an example, and drew output curves within a day of M2 and M3, as shown in Figure 9.It can be seen from Figure 9, compared to M3, that the ES output in M2 is more reasonable, which is mainly reflected in: (1) The output changes in M2 are more stable, and the charge/discharge interval is more stable and works more efficiently.This can extend its life expectancy, which also can be seen from Table 2, (the cost of ES is lower in M2); (2) The ES strategy and real-time electricity price are more closely aligned in M2.ES discharges in the peak price and charges in the valley price.This maximization of the benefits reduces the overall operating costs; (3) Finally, a reasonable charge/discharge of ES will ease the CCHP unit's output pressure, and optimize its operation, resulting in an improvement in energy efficiency.
In summary, the strategy proposed in this paper is more reasonable than the traditional operation method, and provides economic efficiency and high levels of environmental benefits and energy efficiency.

Conclusions
On the basis of the traditional distribution network optimization model, this paper studies the coordinated optimization operation strategy of the regional EI, and uses the improved PSO-IVACO algorithm to solve the issue.The research results show that compared to the traditional distribution network, the EI scheduling process is more complex, and the ES system is particularly important.At the same time, objective factors have become the key components influencing the economic operation in the system, which needs to determine the most reasonable choice that will ensure maximum benefit.
Future research directions include: (1) Make appropriate simplifications for the optimal operation model of CCHP, and study the thermodynamic properties of the system to ensure that the optimization results will be closer to the actual project; It can be seen from Figure 9, compared to M3, that the ES output in M2 is more reasonable, which is mainly reflected in: (1) The output changes in M2 are more stable, and the charge/discharge interval is more stable and works more efficiently.This can extend its life expectancy, which also can be seen from Table 2, (the cost of ES is lower in M2); (2) The ES strategy and real-time electricity price are more closely aligned in M2.ES discharges in the peak price and charges in the valley price.This maximization of the benefits reduces the overall operating costs; (3) Finally, a reasonable charge/discharge of ES will ease the CCHP unit's output pressure, and optimize its operation, resulting in an improvement in energy efficiency.
In summary, the strategy proposed in this paper is more reasonable than the traditional operation method, and provides economic efficiency and high levels of environmental benefits and energy efficiency.

Conclusions
On the basis of the traditional distribution network optimization model, this paper studies the coordinated optimization operation strategy of the regional EI, and uses the improved PSO-IVACO algorithm to solve the issue.The research results show that compared to the traditional distribution network, the EI scheduling process is more complex, and the ES system is particularly important.At the same time, objective factors have become the key components influencing the economic operation in the system, which needs to determine the most reasonable choice that will ensure maximum benefit.
Future research directions include: (1) Make appropriate simplifications for the optimal operation model of CCHP, and study the thermodynamic properties of the system to ensure that the optimization results will be closer to the actual project; (2) The operational state of the regional EI is complex and changeable, and there may be some deviation in the offline scheduling strategy, and further research on the online scheduling method is needed; (3) The proposed scheduling strategy can be used as a reference, which can be used to study the actual ES scheduling strategy.

Figure 1 .
Figure 1.The structure of regional energy internet.The energy input and output of the EI are shown in Figure2.

Figure 2 .
Figure 2. The flow-chart of energy flow.

Figure 1 .
Figure 1.The structure of regional energy internet.

Figure 1 .
Figure 1.The structure of regional energy internet.The energy input and output of the EI are shown in Figure2.

Figure 2 .
Figure 2. The flow-chart of energy flow.

Figure 2 .
Figure 2. The flow-chart of energy flow.

Figure 5 .
Figure 5. Power inputs of WT, PV and load.

Figure 5 .
Figure 5. Power inputs of WT, PV and load.

Figure 5 .
Figure 5. Power inputs of WT, PV and load.

Figure 9 .
Figure 9. Charge and discharge power. )) A is the current capacity.At the same time, as with all types of ES equipment, the charge and discharge power should meet the maximum constraints.Take heat storage as an example, when discharge,

Table 2 .
The results of the operation in one day.

Table 2 .
The results of the operation in one day.

Table 2 .
The results of the operation in one day.

Table A2 .
The spot sale price.