Integrated Energy System Configuration Optimization for Multi-Zone Heat-Supply Network Interaction

: The integrated energy system effectively improves the comprehensive utilization of energy through cascade utilization and coordinated scheduling of various types of energy. Based on the independent integrated energy system, the thermal network interaction between different load characteristic regions is introduced, requiring a minimum thermal grid construction cost, CCHP investment operation cost and carbon emission tax as the comprehensive optimization targets, and making overall optimization to the configuration and operation of the multi-region integrated energy systems. This paper focuses on the planning of equipment capacity of multi-region integrated energy system based on a CCHP system and heat network. Combined with the above comprehensive target and heat network model, a mixed integer linear programming model for a multi-region CCHP system capacity collaborative optimization configuration is established. The integrated energy system, just a numerical model solved with the LINGO software, is presented. Taking a typical urban area in Shanghai as an example, the simulation results show that the integrated energy system with multi-zone heat-suply network interaction compared to the single area CCHP model improved the clean energy utilization of the system, rationally allocates equipment capacity, promotes the local consumption of distributed energy, and provides better overall system benefits.


Introduction
A regional integrated energy system (RIES) is a comprehensive regional energy supply network formed by the coupling of single energy systems such as electricity, gas and thermal (cold). It is a clean, economic, efficient and environmentally friendly energy supply system at the present stage [1,2].The multi-energy complementarity, synergistic optimization and energy cascade utilization of the integrated energy system have improved the comprehensive energy efficiency and reduced the emissions of air pollutants. At the same time, the integrated energy link among multiple regions balances the difference of energy use between regions, cuts peak loads and fills valleys, and improves energy supply reliability [3,4]. However, the complex structure of comprehensive energy system, the coupling of various energy sources, and the matching of installed capacity directly affect the economy of the system operation and the adjustment strategy of unit operations, which is of great theoretical significance and application value for the allocation of multiple energy resources and the optimization of operation strategies in the comprehensive energy system [5][6][7][8][9].
Reference [10] comprehensively coordinated multi-energy forms on the supply side and demand side, and carried out coordinated planning by adopting the grid method to achieve multi-energy complementarity and energy cascade utilization. Reference [11] fully considered the uncertainty of renewable energy output and terminal loads in the integrated energy system, and studied the scheduling optimization of integrated energy systems based on interval linear programming. In [12,13], under different operating modes and different load structures, the optimal configuration and comprehensive operation efficiency of the CCHP system combined with energy storage devices are studied. References [14,15] consider the coordinated planning and operation of the multi-area CCHP system of the heat-supply network model, which improves the gas turbine utilization rate, reduces the gas boiler configuration capacity and the thermal energy transmission loss, and significantly reduces the operating cost. In [16], the energy storage device concept is introduced into the distributed coordination system to optimize the configuration of the mixed integer linear programming model to realize the synchronization optimization of the system structure and operation, and the simultaneous optimization of each device and the energy storage device. References [17][18][19][20] combine solar energy with CCHP systems to optimize the number of units and gas turbine capacity of the optimal joint supply system under different operational control strategies, and achieves the maximum comprehensive benefit of the multi-target joint supply system. Reference [21] proposed a multi-objective optimization model for urban integrated electrical power, thermal and gas grids, which is used for control optimization of modified PRS integrated with thermal users, and thermally integrates the entire system with the district heating network, the effectiveness of the model in terms of economic and environmental performance is quantified by software. With the advancement of the carbon emission trading mechanism, references [22][23][24] comprehensively consider the impact of carbon emissions on the CCHP system, and establish a low-carbon scheduling multi-objective optimization model for carbon trading costs, fuel costs and environmental costs. Reference [25] integrates renewable energy (RE) into an autonomous CCHP system to simultaneously achieve zero environmental emissions and higher power generation and energy efficiency advantages, using an evolutionary particle swarm optimization algorithm to optimize the different configuration size of the autonomous RE-CCHP system. Reference [26] uses the Analytic Hierarchy Process (AHP) to optimize the configuration of hybrid CCHP systems considering three objective functions: annual operating cost ratio (AOCR), primary energy saving ratio (PESR) and carbon emission reduction rate (CERR). Coupling of cold, heat and electric power loads between multi-region CCHP systems, coordinated planning and optimized operation of multi-zone systems, higher utilization of equipment than single CCHP systems, reduced configuration capacity, and significantly reduced operating costs, enabling multiple systems "Multiple horizontal complementarity, vertical source network charge and storage coordination". However, at present, there are few studies on multi-energy networks such as inter-area multi-energy flow access to grids, heat networks, gas networks, etc., in terms of mixed power flow and optimal scheduling among different energy networks [27][28][29][30].
The above references mostly optimize the operation scheduling of integrated energy systems under the condition of known equipment capacity, and research on the collaborative planning of equipment capacity between various energy sources and between multiple regions is lacking. In this paper, a heat network is ingeniously introduced to regions of different load characteristics, In the references, the heating network construction of a single area is generally carried out. Based on the thermal energy-flow constraint, a simplified heat network model is established, which implements thermal energy coupling between different regions. The CCHP system takes into account the combination of devices with different characteristics, establishes an objective function which is considered operating cost, unit investment cost and carbon emission tax. The thermal energy interaction and pipeline flow changes of CCHP systems in various regions are analyzed with examples. The collaborative configuration optimizes regions with different load characteristics. The capacity allocation and network structure of various energy facilities are optimized separately to minimize the total energy supply cost of the system. In the aspect of configuration planning optimization, the redundancy of similar equipment between regions is reduced, and the centralized energy supply characteristics of high-efficiency equipment are reflected. The cost of inter-regional collaborative optimization operation is lower than that of the original independent system operation, the carbon emission tax is reduced, and the clean energy consumption is increased. To a certain extent, the overcapacity or shortage of individual regions is eliminated, and the complementary characteristics of supply and demand between regions are visually reflected.

Integrated Energy System Structure
In the integrated energy system, the three types of loads, including electricity, cold and thermal, are mainly provided by electric power grids, gas turbines (GT), and gas boilers (GB). The photovoltaic generator set and solar collector are added to the system plan to form a comprehensive solar energy utilization and supply system (PVCU CCHP). Multi-energy can be supplied to the cold, thermal and electrical load through different energy transfer mediums "Electrical, Thermal, and Cold Bus" [13]. The multi-energy flow structure of the integrated energy system is shown in Figure 1.

Heat-Supply Network Structure
The heat-supply network is used as the link to connect the thermal load between the CCHP systems. The heat-supply network is laid along the road and arranged in a ring according to the geographical location. When any pipe in annular network is damaged, a gate valve can separate it from other pipelines for maintenance, so it has higher safety reliability. The thermal network in this paper is a simple model, which is established by referring to the network node interactive power model. Figure 2 shows the thermal flow distribution in the heat-supply network pipe. The thermal loss existing in the thermal flow through the pipe is defined as ∆ , . , and , represent the thermal energy at both ends of the pipe at time .

Heat-Supply Network Model
This model is based on the basic principle of heat transfer and pipe networks, with the interacting power ( , , ) between each area and heat-supply network, the flow of heat , , the node flow , as the main factors. The optimization variables are based on the energy conservation law (flow balance) constraint, and the following simplified linearized hot water network model is established. In the short-distance region, assuming that the temperature field of the hot network is a steady-state field, the interaction node has no heat loss, and the heat energy loss is only related to the flow distance of the heat flow in the pipeline, and the feed water temperature and the return water temperature are constant. The time horizon and operating cost model applied are based on 24 h a day use.
(2) Thermal energy loss equation: , , − ∆ , = 0 where ∆ , is the heat loss of the pipeline, , , and , , indicate the mutual thermal energy at two ends of the pipeline at time .
is the thermal energy loss rate of pipeline for unit length, taking 0.1. | , , | > | , , | represents thermal flow from node to node , and | , , | < | , , | represents thermal flow from node to node . Equation (4) means that the sum of interactive power and pipeline thermal loss is 0.

Operation Cost of Heat-Supply Network Investment
(1) Investment cost of the pipeline: where , is the cost of laying pipeline and , is model cost of pipeline. is the length of pipe from to in heat-supply network.
(2) Electric charge for pump operation: where m is the number of CCHP co-supply systems. , , is the unit price of power purchase at time for the region . h is the ratio of electricity consumption to transferred thermal quantity, which means the amount of electricity consumed in the transmission of unit of thermal, taking 0.0059. ∆ is the time interval, take ∆ = 1 h. is the number of running time periods.

CCHP Mathematical Model
The lowest total cost is just the solution of the objective function. The proposal doesn't refer to a multi-objective method. In the process of programming and solving, available energy, various types of load and the cost of unit power are input variables. The total cost of output, various types of interactive power, the rated capacity of each unit and the output power of the unit at each moment are output variables.

Objective Function
where K is the unit capacity price of equipment. and are and alternative model sets, respectively.
(2) Operating cost: , , where , is the total operating cost of the CCHP system. , , is the electricity purchase price of the electric grid. is the ratio coefficient of operation and maintenance cost, taking 0.03. is the annual maintenance cost of the system. , , is the interactive electric power between the system and the electric grid.
. , is the output power of the gas turbine, and . , is the heat energy output of the gas boiler. h , and h , are the conversion efficiency of gas turbine and gas boiler, respectively.
is the price of natural gas, taking 2.37¥/m 3 . is the calorific value of natural gas, taking 9.78 (kw · h)/m , , is the electrical storage device running loss cost, taking 0.02¥/(kw · h). , , is the charge and discharge power of the electrical storage device. (3) emission tax model: where and are electricity emission factor and use gas emission factor, respectively, taking 0.8 kg/(kw • h) and 0.19 kg/(kw • h).
is the emission tax, taking 0.2 RMB/kg. This paper will combine the above optimization objectives, investment cost, operating cost and environmental cost, as the total objective function, as shown in the following expression: where is the integrated cost, and the lowest objective function value is the highest benefit. The model does not involve a separate optimization discussion of each sub-objective functions in the comprehensive target, and does not discuss the weighting of the sub-functions in the total objective function.

Constraints
(1) Power balance constraint. The electric balance, thermal balance and cold balance constraints at each time slot of each load area are: , , = , , , , = . .
where , , , , , , , , are the user's demand for electricity, thermal and cold loads in area at time .
, , is the electric power consumed by electric refrigeration; , , is the output power of photovoltaic generator set. . . is the heat energy power consumed by absorption refrigeration; , , is the output power of solar collector.
. and . are respectively the refrigeration efficiency of electric refrigerator and absorption refrigerator for corresponding models in area .
(2) Interactive power constraint: where , and , are respectively the lower limit and upper limit for the electric power of the tie line.
, and , are the lower and upper bound value of heat-supply network circulating power, respectively.
(3) Upper and lower bound for unit output: where , , and , , are the upper and lower bound of type r GT output, respectively. , , and , , are upper and lower bound of type r GB output, respectively.
where , is the storage state of the electrical storage device.
, , is the charging and discharging electric power at the corresponding time. In order to guarantee the service life and efficiency of the electrical storage device, the lower limit of the storage capacity of the electrical storage device is 0.15 of the total capacity, and the upper limit is 0.9 of the total capacity. Charging and discharging electric power shall not exceed 30% of the total capacity. At the same time, it is guaranteed that the storage state of the initial and ending time is the same.
(5) Photovoltaic output constraint: where , , is the available solar energy for each load area and h is the efficiency of photovoltaic power generation. h is the collector conversion efficiency. The upper limit of , , and , , are the rated power of the equipment. In conclusion, all models are integer linear models, so Lingo software is used for calculation optimization [16]. The parameters involved in this paper are shown in the Appendix.

Results
This paper takes the integrated energy planning of a typical urban area in Shanghai (Caoxi area) as an example, and divides it into three load areas, including commercial area, office area and residential area according to local industrial and load characteristics, as shown in Figure 2. And all kinds of time-of-use electricity prices in Shanghai are shown in Appendix Table A2. Caoxi region is determined by grid division of the Shanghai urban electric power network. The region is an independent electric power supply ring network unit. On this basis, thermal and cold energy supply are divided into blocks and centralized scheduling, photovoltaic power generating set and photovoltaic collector are arranged on the sunny surface of each region. Among them, the residential area has relatively high heat-to-electric ratio, which is between 1.3 and 3. The heat-to-electric ratio in the office area tends to 1 and the heat-to-electric ratio in the commercial area is around 0.5. The equipment parameters are shown in Appendix Table A1, Various load values are shown in Appendix Table A4.

Equipment Capacity Planning
The CCHP system without thermal network interaction is set as mode 1, and the CCHP system with thermal network interaction is set as mode 2. In the planning, considering the actual load and operation characteristics, the equipment with high utilization rate or high economic efficiency and meeting the load requirements is selected according to the optimization results, the multi-energy flow structure of the integrated energy system is redesigned. The optimized configuration and one-time investment cost contrast are shown in Table 1.
It can be seen from Table 1, after considering the heat-supply network, the configuration capacity of gas turbine of users (commercial) with low heat-to-electric ratio increases significantly, while users (residents) with high heat-to-electric ratio do not allocate gas turbine and gas boiler. When considering the heat-supply network, the users with relatively low heat-to-electric ratio can increase the output of the gas turbine to meet their own demand and reduce the electric power purchase, while the surplus thermal energy can be used to subsidize other users with relatively high heat-to-electric ratio through the heat-supply network. Therefore, the gas turbine capacity configuration in the commercial area under mode 2 increases. At the same time, for the residents with high heat-to-electric ratio in mode 1, the insufficient thermal energy can only be reignited through the gas boiler. While combined with the heat-supply network, the thermal energy can be supplied by offices and business districts, so as to reduce the output of gas boilers and reduce the configuration of gas boilers. On the other hand, the users with low heat-to-electric ratio are given priority to deploy gas turbines, and the users with high heat-to-electric ratio are given priority to deploy gas boiler. With the increase of heating network, the overall thermal power ratio of the original high thermal power ratio area will be reduced, and gas turbines will be selected for economy. Photovoltaic output and load characteristics are highly coupled, so the number of photovoltaic generating units has been greatly increased to promote the consumption of clean energy. After connecting to the heating network, the non-gas boiler configuration in residential areas increases the capacity of the photovoltaic collector to supplement the increase of daytime heat load; the gas turbine output in office and commercial areas is increased and the thermal energy supply is sufficient. Under the economic target, the two areas omit the photovoltaic collector configuration. When planning the pipeline layout, the heat-supply network fluid flow , is obtained using Equation (1). Then, the cross sectional area of the pipeline (to limit the fluid flow velocity in the pipeline) can be calculated from , < , , . The temperature of feed water is taken as 100 ℃, the backwater is 70 ℃, and the peak velocity is 0.6 m/s. The planning results are shown in Table  2.
By comparing the investment costs of the two modes, it can be seen that mode 2 saves capacity allocation of gas turbine and gas boiler, and mainly increases the cost of setting photovoltaic generating and laying heat-supply network pipelines. The total investment cost of mode 2 is 6859 million yuan, which is lower than that of mode 1 about 6969.7 million yuan.

Operation Optimization
Through the above optimized configuration, the capacity of each regional equipment in different cases is determined. Based on the previous research, the operating output of each unit is analyzed as follows.
(1) Thermal network interaction energy and flow direction As shown in Figure 3, after adding the heat-supply network, the commercial area with low heat-to-electric ratio generates a large amount of excess thermal energy to supply thermal to the heat-supply network while ensuring enough power supply for local area. The resident users with higher heat-to-electric ratio absorb thermal energy from the heat-supply network. The office area is close to the gas turbine heat-to-electric ratio at 9:00-13:00, and the thermal energy surplus is less and there is no thermal load during the 20:00-8:00 period. The gas turbine keeps running due to the electrical load, and supplies thermal to the heat-supply network. At 14:00-17:00, the office area absorbs thermal from the heat-supply network. The interaction between CCHP system and the heat supply network accurately reflects the configuration and operation of the energy supply equipment. It is stipulated that the flow direction of commercial and office districts to residential areas should be positive, and the flow direction of commercial districts to office areas should be positive. In the heat-supply network planning, about 30% allowance of pipeline radius is left. It can be seen from Figure 4a that the flow direction is positive in each time period, that is, the pipeline flows from commercial area to residential area in 24 h. The flow rate of the pipeline is 0 during the period 13:00-18:00 with no thermal flow, and the other times are positive. The thermal energy is transported to residential areas from the office area during this period. The pipeline is only positive at time period 13:00-18:00, and the thermal energy flows from commercial area to office area, and there is no thermal interaction at other time. Combined with the energy flow interaction analysis at each time period in Figure 3, the thermal transfer capacity of commercial and residential pipelines is much higher than that of other pipelines. However, by planning the radius of thermal network pipelines, the velocity of each pipeline is maintained at a reasonable value. This laying method reduces the inert effects and thermal loss of thermal transfer, while improving the efficiency of pipeline utilization and investment economy.
It can be seen from above that the utilization rate of commercial-office heat-supply network pipe is low, and the pipe laying is far away and the cost is high. Therefore, considering that only the commercial-residential-office pipe should be laid, the pipeline velocity is shown in Figure 4b. After omitting the pipe , the radius of pipe increases from 5 to 6, the radius of pipe remains unchanged, the pipe utilization rate increases, and the laying cost decreases.  (2) Comparison of interaction with electric power grid Figure 5 shows that when the residential area is in the time period 0:00-9:00 electrical load valley section in mode 1, the gas turbine operation meets the regional electric power demand. In mode 2, the electric energy is provided by the electric power grid and photovoltaic unit. During the period of 12:00-15:00, due to the increase in PV configuration capacity and the increase in photovoltaic power generation, the purchase of electricity is slightly reduced.
While the office area with relatively low heat-to-electric ratio is supplied to the thermal energy of the region, except for the period 14:00-17:00 office area gas turbine decline in output, it has increased the purchase of electricity to the electric grid. During the period from 23:00-8:00, the residential area has less electricity and high thermal load. While the commercial area has a certain electrical load and the thermal load tends to zero. At this time, the gas turbine output increases and complements the energy characteristics of the residential area. Therefore, the utilization rate of the gas turbine is increased, and it is more economical to use the distributed power source preferentially, and the difference in energy consumption is supplemented by the electric power grid.
(3) Electric power output distribution In Figure 6a, residential area system power load is mainly supplied by the electric grid, and makes full use of photovoltaic power generation during the period 8:00-16:00, while the electricity price in the valley section 21:00-7:00, 12:00-17: 00 charges the electrical storage device. During the peak period of 8:00-11:00 and 18:00-20:00, the absorption refrigeration unit replaces the electric refrigeration unit, and the residential area absorbs more thermal from the heat-supply network. At the same time, the charged electrical storage device discharges during this period to reduce the purchase of electricity from the electric power grid and effectively reduce the peak-filling effect.
In Figure 6b, the difference between day and night electrical load of the office system is obvious. The output of the gas turbine is almost the same as that of the electric power grid. The waste thermal is recovered and supplied to the thermal load or the thermal load is supplied to the residential area through the heat-supply network. Only absorption refrigeration units and no electric refrigeration units are installed in the planning, and the output of absorption refrigeration units is maintained at a higher level during the operation. We fully utilize photovoltaic power generation from 7:00-17:00, and charge the electrical storage device at 23:00-5:00 and 13:00-18:00. During the period 8:00-12:00 and 19:00-21:00, the electrical storage device is discharged during the peak period of electricity price. Similarly, in the commercial area with high heat-to-electric ratio in Figure 6c, the system electric power load is mainly satisfied by the purchase of electricity from the gas turbine and electric power grid. The gas turbine assumes the main power supply task. The electrical storage device charge and discharge and photovoltaic power generation operation conditions in this area are similar to the above two regions. The refrigeration unit operates at 23:00-5:00 during the electricity price period of the valley section, and the cooling load is independently satisfied by the absorption refrigeration unit during other periods. The load characteristics of residential and commercial areas and the analysis of the operating conditions of the equipment after planning show that the load characteristics of the two areas have strong complementary correlation characteristics. Compared with the operation of the original independent CCHP system, the joint regional planning for joining the heat-supply network reduces the repeatability of the equipment in each area while enhancing the multi-energy complementarity. The division of energy between the areas with load characteristics is (4) Sensitivity of natural gas price Figure 7 shows that as the price of natural gas rising, the operating costs of various regions also increases. The operating costs of commercial areas with increased gas turbine capacity after adding heat-supply networks increased obviously. In the later period of rising natural gas prices, the cost-saving curve is flat and the savings are decreasing which means the increase in natural gas prices has weakened the advantages of the heat-supply network. Due to the lack of electric heat-supply equipment in the system, even if the price of natural gas rises, the gas turbine must operate to maintain the energy demand of the thermal load. Therefore, when the gas turbines of each system maintain low-line operation, the utilization efficiency of the equipment is low, and the effect of increasing the comprehensive benefit of the heat-supply network to the system is limited.

(5) Operating expenses
It can be seen from Table 3, in the planned residential areas, gas turbines and gas boilers were abandoned, the photovoltaic capacity allocation was increased, and the carbon emission tax was greatly reduced. In the office area, the difference in equipment capacity is small under mode 1 and mode 2, and the output of the equipment is similar. The carbon emission tax fluctuations in different modes are extremely small. After adding the thermal network, the gas turbine capacity in the commercial area increased, the output of the equipment increased greatly during operation, and reduced the electricity purchase from the electric grid. Compared with the large increase in output, the carbon emission tax increased slightly which indicates the cleanliness of the gas turbine operation. The capacity allocation of photovoltaic generator sets and photovoltaic collectors have a certain impact on the comprehensive benefits of the carbon emission tax. However, the two projects have a large investment cost and the unit output has uncertainty. After the addition of the heat-supply network, the total operating costs of the three regions were reduced by 8.2% compared with the independent operating costs and the carbon emission tax was reduced by 3.32%.

Discussion
According to the operating characteristics, it can be seen that the gas equipment in the independent area maintains low-load operation, and the economic dispatching is limited. The joint planning of joining the heating network reduces the repeatability of equipment in each area, and the  operating cost (¥) division of energy between the areas with load characteristics is more clear. It embodies the operating characteristics of centralized energy supply and high utilization of equipment. The multi-region comprehensive energy system based on thermal network interaction not only saves the one-time investment cost, but also improves the overall operation efficiency while promoting photovoltaic absorption. Of course, the application of the heat network has its limitations. If the load characteristics between the regions are similar, the effect of increasing the heat network is not obvious, the heat storage device can be considered. And the distance between the regions should not be too long to reduce the heat energy transmission loss. In cold places, it can be used as a supplement to the distributed thermal energy system as a thermal power plant system, and the benefits brought by the construction of the heat network in the different load characteristics areas are particularly obvious.
This paper does not discuss the impact of surplus electric power inputted to electric grid and prospective annual load change on the comprehensive energy capacity allocation plan, and it is also simple to model the thermal network. The evaluation of the operation state of the comprehensive energy system, such as the primary energy utilization efficiency and environmental benefits, can be taken as the research direction of the next stage.

Conclusions
In this paper, the integrated energy system of typical urban residential, office and commercial mixed areas in the city is connected through the heat-supply network for joint planning, and compared with the original independent regional integrated energy system planning and operation, the following conclusions are drawn: (1) When natural gas prices fluctuate greatly, it is not recommended to configure large-capacity gas equipment, and the appropriate capacity can be configured to meet the minimum load requirements. The load characteristics of residential and commercial areas and the analysis of the operating conditions of the equipment after planning show that the two areas are highly complementary. For areas with low heat-to-electric ratio, the combination of gas turbine and absorption chiller is prioritized to meet the cold and hot load. For areas with high heat-to-electric ratio, the output of gas turbine and gas boiler should be balanced to achieve optimal working coordination. Inter-regional heat energy interacts through the heat network, enabling multi-energy flow to achieve cross-regional distribution, coordinated planning of multi-regional systems, Although the cost of heat network pipe laying has increased about 2.44 million yuan, The total investment cost of mode 2 is lower than that of mode 1 about 1.1 million yuan, reduce proportion is 1.6%, especially the investment cost of gas turbine has been reduced about 7.3 million yuan, reduce proportion is 13.3%. Significantly reduced the configuration and investment of gas turbine and gas boiler, resulting in lower equipment costs.
(2) The primary investment cost of the electrical storage device is expensive, and the lower economic benefit brought by the electrical storage device in the area to adjust the output of the gas turbine. Therefore, the electrical storage device is mainly used to interact with the power grid, and the peak value is filled with the price response. The integration of energy storage and heat network into multi-regional joint planning can achieve multi-energy scheduling in time and space between regions. In the modeling process, a carbon emission model was added. The investment cost of photovoltaic generating set has been added about 4.4 million yuan, increased by 66.7%, it means the power of PV output increased and promotes renewable energy consumption. After the addition of the heat-supply network, the total operating costs of the three regions were reduced by 8.2% compared with the independent operating costs and the carbon emission tax was reduced by 3.32%. The simulation showed that the input of the heating network reduced the overall energy consumption of the system, promoted the use of clean energy, and improved the local consumption capacity of distributed energy.

Conflicts of Interest:
The authors declare no conflicts of interest Appendix A