Abstract
In this paper, aiming to achieve the target of carbon emission orientation, a multi-objective optimization model of the multi-energy flow coupling system is proposed, in which all the environmental protection, system economy, and energy efficiency are comprehensively considered as the addressed objectives. To solve the developed model, by combining the analytic hierarchy process (AHP) and the improved entropy weight method, a so-called AHP-improved entropy weight method is proposed and utilized for weighting the considered objectives, and the model is transformed into a single objective optimization problem, namely, the collaborative optimization model. Then, to expedite the process, a simplified primal dual interior point method is proposed to solve the model. Finally, the results of a case study indicate that the proposed multi-objective collaborative optimization can obtain the optimal solution of the system. In addition, the convergence and global optimization ability of the simplified primal dual interior point method show better characteristics when solving the proposed model.
1 Introduction
Energy is the basis and important guarantee for human survival. There are many problems in traditional energy systems, such as independent energy supply, low cascade utilization level, energy waste, and environmental pollution (; ; ). The multi-energy flow coupling system (MEFCS) is an energy system form that integrates public cold, heat, electricity, and gas. Its purpose is to integrate multiple energy sources, such as electric energy, natural gas, and thermal energy in a certain area, so as to realize collaborative optimal operation, collaborative management, and complementary mutual assistance among various forms of energy subsystems (; ). In addition, under the background of “double carbon”, the transformation of clean and low-carbon energy is an inevitable trend of global energy development.
Traditional energy systems are planned and operated independently, and only a single situation needs to be considered in their optimal scheduling. However, for a multi-energy flow coupling system, the correlation among energy subsystems should be considered in planning and operation (). In ), considering the multi-timescale characteristics, an electrical and thermal energy sharing model of interconnected microgrids with combined heating and power (CHP) and photovoltaic systems was built, in which CHP could operate in a hybrid mode by selecting the operating point flexibly. In ), a multi-objective bi-level optimization model considering the total cost and carbon dioxide emission was built, while the energy efficiency of multi-energy flow coupling system was ignored. In ) and ), under the condition of meeting the basic needs of power, gas, and heat loads, the coordinated planning of multi-energy flow coupling system was considered, in order to reduce the construction cost of transmission lines, gas pipelines, and power plants as much as possible. In ), the authors presented an optimization framework for the optimal scheduling of a multi-energy microgrid, where a number of aggregated end-users were considered. In ), a novel mixed integer linear programming optimization algorithm has been developed to compute the optimal management of a micro-energy grid, where the total cost, the NOx, and the CO2 emissions of the system were taken into consideration. To meet the safety constraints, in ), an optimal coordination control strategy (OCCS) for a hybrid energy storage system was developed considering the state-space equation to describe the OCCS, the constraints of the OCCS, and the objective function to express the optimal coordination control performance. In ), the authors considered the day ahead optimal scheduling problem of electricity gas interconnected systems, where the two-way energy flow was taken as a non-convex nonlinear mixed integer linear programming problem, and a second-order cone programming (SOCP) method has been proposed. In ) and ), the uncertainty caused by renewable energy and multi-energy load was considered, and the robust optimization and stochastic optimization methods were adopted to deal with it respectively, so as to ensure that the system can still maintain stable operation under the worst conditions. In ), a grid-connected two mass DFIG and a grid-supportive single mass squirrel cage induction generator-based flywheel energy storage system model have been considered for controller design and proof-of-concept exploration. In ), the flexible resources (FRs) on both the energy supply and load sides were introduced into the optimal dispatch of the integrated electricity-heat energy system (IEHES) and further modeled to alleviate the renewable fluctuations, and the solution for FRs participating in IEHES dispatch was given, with goals of maximizing the renewable penetration ratio and lowering operation costs. It can be seen that most of the existing results consider optimization of the economic objectives of the multi-energy flow coupling system, where the index is relatively single, and less consideration is paid on the carbon emission level in the operation of the system. At the same time, the operation strategy is the lack of comprehensive comparison and verification.
In solving the MEFCS collaborative optimization model, when considering multiple optimization objectives including carbon emission, investment and operation cost, and energy utilization, the traditional single objective optimization algorithm may be difficult to ensure that the solution result is the optimal solution of the original problem. In ), the load characteristics and various constraints of the integrated community energy system were considered, and the operating model with the goal of minimizing operating costs was optimized. In ), the energy consumption cost and environmental cost of the multi-energy flow coupling system were considered comprehensively, the optimal scheduling model of multi-energy flow coupling system was proposed, and the optimal scheduling model was transformed into a mixed integer linear programming problem. In ), the method of the probability scenario had been used to model the uncertainties of the distributed renewable energies (DREs) and loads, which could better characterize the impact of uncertainty on the planning and design of the MEFCS. In ), a two-stage robust generation scheduling model was proposed for the dynamic safety constraints of the natural gas pipeline network and the uncertainty of wind power, and a new solution method was developed to avoid the nonlinearity of gas flow constraints. In ), the multi-objective optimization model was transformed into a single objective optimization model through the multi-objective programming hierarchical solution method, and the primal dual interior point method was used to solve the model. Based on the fast particle swarm optimization algorithm, in ), a dual-decomposition-based distributed algorithm was designed to address the problem that the data and information of the EHs during the operation were confidential and should be kept by each owner, where the optimal consensus problem was used for the dual problem to update the multipliers, in ), the proposed MEFCS planning model, formulated as a two-stage MILP problem, was solved by the Benders decomposition (BD) method to determine the optimal capacity of each component in MEFCS planning.
To be pointed out that, the research on the optimization of multi-energy flow coupling system at home and abroad mainly focuses on the simplification of the optimization model. However, on the one hand, it will lead to the reduction of solution accuracy, at the same time, because the models are more and more complex, which are difficult to be simplified. Therefore, the heuristic algorithm has become an important way to deal with optimization problems. However, the traditional heuristic algorithm has the problems of poor convergence and easy to fall into local optimization, and how to find a simplified and better algorithm is another motivation of this paper. Based on the above discussions, in this paper, the environmental protection goal is taken as the leading factor, the economic and energy efficiency goals are comprehensively considered, the multi-objective collaborative optimization model is developed for the multi-energy flow coupling system, which can be transformed into a single objective optimization model through the linear combination of the analytic hierarchy process and the improved entropy weight method, and then the model can be solved by using the simplified primal dual interior point method. The results avoid falling into local optimization and accelerate convergence. Case studies verify the effectiveness of the proposed algorithm.
2 Modeling of Multi-Energy Flow Coupling System
A typical multi-energy flow coupling system structure is shown in Figure 1, which is internally connected through the power grid, thermal pipe network, and cooling transmission network. The equipment involved distribution power source includes a wind turbine (WT) and photovoltaic (PV). Cogeneration includes CHP, a gas turbine (GT), a waste heat boiler (WHB), a ground source heat pump (HP), an electric refrigerator (ER), an absorption refrigerator (AR), and other energy conversion equipment, as well as electric energy storage (EES), heat energy storage (HES), and other energy storage equipment.
FIGURE 1
2.1 Modeling of Distributed Generations
2.1.1 Wind Turbine
where indicates the wind turbine generation power (kW) in time period t, is the wind energy utilization efficiency of the wind turbine, represents the blade radius (m), represents the air density (), and is the air velocity () in time period t.
2.1.2 Photovoltaic
where refers to the output power (kW) of photovoltaic equipment during the period t, represents the test power (kW) under standard conditions t, refers to the light intensity () in the period t, is the test light intensity () under standard conditions, K is the power temperature coefficient, which is taken as −0.0047; , , and represent the solar panel temperature, reference temperature, and external ambient temperature (), respectively; normally the reference temperature is taken as 25 ; and expresses the solar radiation intensity () in time period t.
2.2 Modeling of Energy Conversion Unit
2.2.1 Cogeneration Unit
The cogeneration unit generates electric energy and heat energy at the same time by consuming natural gas. Its operation mode can be expressed aswhere , , and are the electric power, thermal power, and gas power consumed by the internal cogeneration unit in scheduling period t, respectively, and and are the power generation efficiency and heating efficiency of cogeneration units, respectively.
2.2.2 Gas Turbine and Waste Heat Boiler
The gas turbine generates electric energy by consuming natural gas, and part of the discharged flue gas can be transformed into available calorific value through a waste heat boiler. Their working characteristics can be expressed aswhere and indicate the gas turbine generation power and flue gas waste heat power during the period t, respectively; represents the low calorific value of natural gas, which is set as 9.78 in this paper; expresses the natural gas consumption during the period t; t is the scheduling period; and represent the power generation efficiency and loss rate of gas turbine, respectively; is the recovery efficiency of the waste heat boiler; and is the heat recovery power of the waste heat boiler in time period.
2.2.3 Ground Source Heat Pump
The heat pump is a high-efficiency and energy-saving equipment in the multi-energy flow coupling system. It can convert low-grade heat energy into high-grade heat energy by consuming electric energy. Its operation mode is given bywhere and represent the heat energy generated and electric energy consumed of the ground source heat pump during the period t, respectively; and is the conversion efficiency of the heat pump.
2.2.4 Electric Chiller and Absorption Chiller
The electric chiller generates cold power by consuming electric power during operation, and the absorption chiller generates cold power by absorbing thermal power. Its mathematical model is as follows:where and represent the cool power generated in time period t of the electric chiller and absorption chiller, respectively; and represent the conversion efficiency of the electric chiller and absorption chiller, respectively; indicates the electric energy consumed of the electric chiller during the period t; represents the heat energy consumed of the absorption chiller during the period t; is the rated conversion efficiency of the absorption chiller; , , and are the refrigeration coefficient of the absorption chiller, respectively; and is the load rate when the absorption chiller is working.
2.3 Modeling of Energy Storage Equipment
where represents the energy storage of energy storage equipment i in time period t, and are the charging power and discharging power of energy storage equipment i in time period t, and represent the charging efficiency and discharging efficiency of energy storage equipment i, and is the consumption rate of energy storage equipment i.
3 Modeling of Multi-Objective Collaborative Optimization
In the multi-objective collaborative optimization of MEFCS considering carbon emissions, the optimization objectives considered in this paper include the environmental protection objective, economic objective, and energy efficiency objective.
3.1 Each Optimization Objective Function
3.1.1 Environmental Protection Objective
Aiming at minimizing the emission of MEFCS in 1 day, the optimization model can be established as follows:where and represent the emission coefficient corresponding to the combustion of natural gas and the consumption of electric energy; in this paper, they are taken as 184 g/kWh and 877 g/kWh, respectively.
3.1.2 Economic Objective
In order to minimize the operation cost of MEFCS in 1 day, the optimization model can be formulated aswhere and are the cost coefficients corresponding to the electric energy and natural gas consumed by the system, respectively; is the maintenance cost of equipment i; is the rated capacity of equipment i; and N represents the total amount of equipment.
3.1.3 Energy Efficiency Objective
Primary energy utilization is defined as the ratio of MEFCS load to MEFCS primary energy input in a day. Aiming at the maximum utilization of primary energy, the optimization model can be formulated aswhere , , and represent the total load of the system in a day, respectively, and represents the network loss rate of transmission line, which is usually chosen as 5%.
3.2 Constraint Condition
3.2.1 Energy Balance Constraints
1) Power balance constraint
2) Heat energy balance constraint
3) Cool energy balance constraint
3.2.2 Upper and Lower Limits of Equipment Output
where represents the output power of equipment i in period t.
3.2.3 Energy Storage Constraints
where and represent the upper limit of charging and discharging power of energy storage equipment i, respectively; represents 0–1 variable; is the energy storage of equipment i in period t; and and represent the upper and lower limits of the charging and discharging state of the energy storage equipment i, respectively.
3.3 Collaborative Optimization Objective
The developed optimization model is a multi-objective optimization problem. First, the optimal solution of each objective is obtained through single objective optimization, and then the optimization results of each objective are standardized, so the multi-objective optimization is transformed into single objective optimization with the help of the linear weighting method. Finally, the single objective optimization algorithm can be solved.
3.3.1 Normalization and Standardization
As the environmental protection goal and economic goal belong to very small goals, that is, the smaller the final result, the better, while the energy efficiency goal belongs to maximum goals, the larger the final result, the better. Therefore, before establishing the collaborative optimization objectives, each single objective should be normalized and standardized, which can be expressed as follows:where and represent the membership function of very small target and maximum target, respectively, is the ith objective function, and and are the minimum and maximum of the ith objective function, respectively.
3.3.2 Index Weighting
Generally, the methods of weighting indicators can be divided into subjective method, objective method, and the combination of subjective and objective methods. The subjective weighting method is simple to operate and does not need the support of original data, but the subjectivity of weighting results is often too large. The objective weighting method can show the relationship between indicators well, but it has high requirements for the original data. Therefore, in this paper, a new combination method based on the analytic hierarchy process (AHP) and the improved entropy weight method is adopted.
The analytic hierarchy process first judges the relative importance of each index through decision-making experts and scores each index with an integer between 1 and 9, and then the judgment matrix is obtained,where n denotes the number of indicators. A is a positive reciprocal matrix, which satisfies .
To be noted that, due to the environmental protection goal is taken as the leading factor in this paper, when forming the judgment matrix, the score of the environmental protection index is relatively high so that the final weight is relatively maximum.
Then check the consistency of the judgment matrix,where represents the consistency proportion. If , the consistency verification passes, otherwise the judgment matrix needs to be modified. and represent the consistency index and average random consistency index, respectively. is the maximum eigenvalue of judgment matrix A.
When the judgment matrix A passes the consistency check, the eigenvector corresponding to its maximum eigenvalue is obtained and normalized, that is, the weight vector is obtained by the analytic hierarchy process,
The entropy weight method reflects the amount of information contained in each index through the entropy value of each index. Generally speaking, the smaller the entropy value, the greater the amount of index information and the greater the weight should be set. Since the standard entropy weight method is mainly applied to multiple schemes, the entropy weight method can be improved bywhere represents the characteristic specific gravity of the ith target, is the entropy of the ith target, and is the membership function of the ith objective.
To be noted that, this improvement is mainly to adapt to the optimization model. The objective functions have been processed and converted into the form of membership function. Therefore, in order to adapt to this form, the index value is replaced by membership . Because this is not an evaluation problem, there are no multiple schemes to be evaluated. Therefore, the possible maximum and minimum values of each membership degree are substituted into the formula to reduce the individual deviation. In this way, there are three evaluation schemes in terms of quantity, that is, . Therefore, the above formula is obtained.
According to the calculation results of entropy value of each index, the weight of each index can be obtained by
Therefore, the weight vector is obtained by the improved entropy weight method,
In order to obtain the combined weight of AHP and improved entropy weight method, the coupling vector is taken as follows:where and represent the coupling weight of index I for weight coefficients and , respectively, which can be expressed as
Therefore, the weight after coupling is
Normalize it, one has
Then, the combined weight of the AHP improved entropy weight method can be finally expressed as
3.3.3 Collaborative Optimization Objective
After obtaining the index weight, combined with the standardized objective function in the above sections, we can obtain the comprehensive satisfaction goal, that is, the collaborative optimization objective is given as
4 Optimization Algorithm
Based on the above analysis, it can be found that the multi-objective collaborative optimization of the multi-energy flow coupling system considered in this paper is a complex nonlinear programming problem. In order to make the solution speed and convergence meet the requirements of practical problems, the simplified primal dual interior point algorithm is used in this paper. For the sake of brevity, first, the optimization model described above is transformed into the following general form:where is the state variable, including the output power, external power purchase, gas purchase, etc., of each equipment, is the equality constraint, including the power balance constraint of the system, the energy balance constraint at the beginning and end of the scheduling cycle of energy storage equipment, etc.; is the inequality constraint, including the upper and lower limits of equipment output, energy storage charge and discharge constraints, etc.; and and represent the upper and lower bounds of the inequality, respectively.
When dealing with this optimization model with the traditional interior point algorithm, relaxation variables and are introduced first, where r represents the number of inequality constraints; thus, the original inequality constraints are transformed into the equality constraints. The resulting optimization model is formulated as
At the same time, the size of relaxation variables and should be restricted to ensure that the objective function is always far away from the solution boundary so that it can be solved in the feasible domain as follows:where represents the introduced disturbance factor.
At this point, the inequality constraints contained in the optimization model described in this paper have all been converted into the equality constraints, and the Lagrange function for this optimization problem can be expressed aswhere , , and all represent the Lagrange operators, also known as the dual variables. By deriving this Lagrange function, the optimal solution to this optimization problem can be obtained.
In this paper, by simplifying the original dual interior point algorithm, the simplified original dual interior point method can be utilized to solve the optimization model, and the simplified process is to rewrite the inequality constraint towhere and represent the generalized inequality constraints and generalized upper bounds, respectively.
It can be found from the traditional interior point algorithm that in the process of dealing with inequality constraints, the upper and lower bounds of inequality constraints need to be relaxed, and then converted into equality constraints, respectively. At the same time, Lagrange operators are also introduced for equality constraints converted from upper-bound inequality constraints and lower-bound inequality constraints, respectively, which introduce more variables in the Lagrange function. This simplification algorithm greatly reduces the relaxation variables and corresponding Lagrange operators introduced in the optimization model, improves the convergence speed of the algorithm while guaranteeing the calculation accuracy, and reduces the amount of programming to a certain extent. The rest of the algorithm is handled similarly to the traditional algorithm, which are not discussed here. For ease of understanding, the multi-objective collaborative optimization calculation flow of the multi-energy flow coupling system based on the simplified primal dual interior point algorithm is shown in Figure 2, where is the dual gap, represents the convergence accuracy, which is taken as 10–6 in this paper, and represents the maximum number of iterations, normally is set as 300.
FIGURE 2
Based on the above analysis, the collaborative optimization of MEFCS dominated by the carbon emission targets in this paper can be summarized as the following steps:
Step 1. Enter parameter information of the multi-energy flow coupling system, such as system load, rated capacity of each unit, equipment parameters, and time-of-use electricity price.
Step 2. Establish the steady-state operation model of each equipment, as shown in Eqs 1–7.
Step 3. Establish the objective functions and constraints of the MEFCS, as shown in Eqs 8–15.
Step 4. Weight each objective function using the analytic hierarchy process-improved entropy weight method, as shown in Eqs 18–28.
Step 5. Convert each objective function into a collaborative optimization objective through the membership function and the obtained weight information, as shown in Eqs 16, 17, 29.
Step 6. Solve the MEFCS collaborative optimization model by the primal dual interior point algorithm until the optimal solution is obtained or the algorithm does not converge.
The above steps can be clearly represented by the flowchart shown in Figure 3.
FIGURE 3
5 Case Study
5.1 Case Description
In this paper, the typical multi-energy flow coupling system shown in Figure 1 is selected as an example. The capacity of each equipment is as follows: one photovoltaic generator unit with a rated output of 700 kW and one wind turbine generator unit with a rated output of 500 kW, one cogeneration unit with a rated output of 3 MW, one gas turbine with a rated output of 2 MW, one waste heat boiler with a rated output of 1 MW, four heat pumps with a rated output of 500 kW, four electric refrigerators and four absorption refrigerators with a rated output of 200 kW, and four batteries and heat storage equipment with a rated capacity of 500 kwh. Other economic and technical parameters of the equipment can be found in ) and ). The time of use electricity price information of the multi-energy flow coupling system is shown in Figure 4, and the price of natural gas is 2.71 yuan/m3 (). The load data of the system are shown in Figure 5.
FIGURE 4
FIGURE 5
5.2 Results Analysis
The output curve of each equipment is shown in Figure 6. It can be seen that the supply of electric energy and heat energy of the system is mainly guaranteed by a gas turbine and heat pump, but only the output of each equipment is not enough to meet the load demand during the peak load period of the system. At this time, the system needs to purchase electricity from the external power grid to jointly supply energy to the load. In addition, it can be seen from the figure that the discharge time of the power storage equipment is 03:00–21:00, and the heat release time of the heat storage equipment is 06:00–19:00. In other periods, that energy storage equipment is in the charged states.
FIGURE 6
While using the multi-objective collaborative optimization model proposed in this paper to solve the multi-energy flow coupling system, three separate objectives are solved respectively. After finding the individual optimization of each objective, its state variables are substituted into the other two objectives to obtain the respective results of the three objectives in this case. The comparison between the operation results of each scheme and the operation results of multi-objective collaborative optimization is shown in Table 1. It can be found from Table 1 that the CO2 emission under multi-objective collaborative optimization increases by 3.8% compared with that under single objective F1 optimization, the system operation cost increases by 6.1% compared with that under single objective F2 optimization, and the primary energy utilization rate decreases by 7.2% compared with that under single objective F3 optimization. Although each objective under multi-objective collaborative optimization is not the optimal solution, the contradiction and conflict between each single objective are balanced in the optimization process. On the premise of taking the minimum carbon emission of the system as the leading objective, the comprehensive satisfaction of the system is significantly higher than the solution results of each single objective, and a relatively satisfactory optimal scheduling scheme is given.
TABLE 1
| Operation form | Multi-objective collaborative optimization | Single objective optimization | ||
|---|---|---|---|---|
| F1 optimal | F2 optimal | F3 optimal | ||
| F | 0.9685 | 0.8943 | 0.8304 | 0.8612 |
| F1 (kg) | 10,403.5 | 10,022.2 | 15,113.3 | 10,543.7 |
| F2 (yuan) | 14,568.8 | 16,205.2 | 13,727.1 | 14,037.3 |
| F3 (%) | 79.04 | 81.67 | 67.94 | 86.33 |
Comparison between multi-objective collaborative optimization and single-objective optimization.
In order to further compare the differences between the multi-objective collaborative optimization proposed in this paper and the traditional single objective optimization, the optimization results of each single objective and the multi-objective collaborative optimization results are analyzed period by period, as shown in Figure 7. In the figure, , , and represent the environmental protection objective, economic objective, and energy efficiency objective under collaborative optimization, while , , and represent the environmental protection objective, economic objective, and energy efficiency objective under single objective optimization. Figure 7A shows the carbon emissions in each period of the two optimization methods. It can be seen from the figure that the two carbon emission curves cross each other. Except that the carbon emissions during collaborative optimization in 17:00–21:00 are significantly higher than those in single objective optimization, they are very close in other times. It shows that when taking the minimum carbon emission as the leading objective, the effect of collaborative optimization is not different from the single objective optimization with the minimum carbon emission, and the working state of each equipment is also relatively stable. Figure 7B shows the cost curves of the two optimization methods. During 3:00–12:00, the cost of multi-objective collaborative optimization is about 100 yuan/h higher than that of single objective optimization, and the two curves almost coincide after 12:00. It can be seen from Table 1 that the comprehensive satisfaction of multi-objective collaborative optimization is obviously higher than that of single objective optimization. On this basis, it ensures that the operation cost of the system is not too high, and it is almost the same as that of single objective optimization in most periods, indicating that the result of multi-objective collaborative optimization is ideal. Figure 7C shows the comparison of energy efficiency of the two methods in each period. Although the operation energy efficiency of multi-objective optimization in each period is not as good as that of single objective energy efficiency optimization, the overall operation result is relatively stable, indicating that each equipment can achieve stable energy supply and continuous output during the operation of the system, and the working state is not easy to fluctuate violently.
FIGURE 7
In order to highlight the effectiveness of the simplified primal dual interior point method proposed in this paper, the particle swarm optimization (PSO) algorithm is selected to compare with the algorithm proposed in this paper. The solution process curves of the simplified primal dual interior point method and particle swarm optimization algorithm for the system comprehensive satisfaction objective are shown in Figures 8A,B, respectively. They tend to converge at the 25th and 65th iterations, respectively. It can be seen that the convergence of the simplified primal dual interior point method is better than that of the particle swarm optimization algorithm. In addition, from the solution results of the two algorithms, it can be seen that the simplified primal dual interior point method finally converges near 0.9685, while the particle swarm optimization algorithm finally converges only near 0.8352, which still has a certain deviation from the global optimal solution. Therefore, the global optimization ability of the simplified primal dual interior point method proposed in this paper is also stronger than that of the particle swarm optimization algorithm.
FIGURE 8
6 Conclusion
In this paper, the multi-objective collaborative optimization model of MEFCS has been developed considering each of the environmental protection, system economy, and energy efficiency as the objectives of this study, in which the carbon emission orientation goal could be achieved. Moreover, the simplified primal dual interior point method has been used to solve the constructed model. According to the obtained results, the following findings have been concluded: 1) The contradiction and conflict between the three objectives (CO2 emission, system operation cost, and primary energy utilization rate) were relatively balanced under the proposed collaborative optimization operation, which have clearly demonstrated that the satisfaction of collaborative optimization operation considering multiple objectives could be higher than that considering a single objective of the system. 2) Each equipment of the system could achieve stable energy supply as well as continuous output throughout the whole operation process, whereas the working state was difficult to fluctuate sorely. 3) At the same time, the operator could adjust the weight of the three objectives through his own will, so as to get the best operation results that meet his requirements. 4) In addition, the simulation results have illustrated that the simplified primal dual interior point method being adopted in this paper has better convergence and global optimization ability in multi-objective collaborative optimization. However, more investigations are needed regarding the sensitivity analysis on the critical parameters of the multi-energy flow coupling system, which will be considered in our future research work.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
XZ: Methodology, software, writing—original draft, data curation. YY: Conceptualization of this study, supervision, review and editing. HW: Review and editing.
Funding
This work was supported by the National Natural Science Foundation of China (No. 51477041).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Glossary
- Acronyms
- AHP
Analytic hierarchy process
- AR
Absorption refrigerator
- BD
Benders decomposition
- CHP
Combined heating and power
- DRE
Distributed renewable energy
- EES
Electric energy storage
- ER
Electric refrigerator
- FRs
Flexible resources
- GT
Gas turbine
- HP
Heat pump
- HES
Heat energy storage
- IEHES
Integrated electricity-heat energy system
- MEFCS
Multi-energy flow coupling system
- OCCS
Optimal coordination control strategy
- PV
Photovoltaic
- SOCP
Second-order cone programming
- WHB
Waste heat boiler
- WT
Wind turbine
- Parameters
Wind energy utilization efficiency
Blade radius
Air density
Test power under standard conditions
Test light intensity under standard conditions
- K
Power temperature coefficient
Power generation efficiency
Heating efficiency of cogeneration units
Low calorific value of natural gas
Power generation efficiency
Loss rate of gas turbine
Recovery efficiency of waste heat boiler
Conversion efficiency of the heat pump
Conversion efficiency of electric chiller and absorption chiller
Conversion efficiency of absorption chiller
Refrigeration coefficients of absorption chiller
Load rate
CO2 emission coefficient corresponding to the combustion of natural gas and the consumption of electric energy
Cost coefficients corresponding to the electric energy and natural gas consumed by the system
- N
Total amount of equipment
Total load of the system in a day
Network loss rate of transmission line
Consistency proportion
Consistency index
Average random consistency index
- Variables
Wind turbine generation power in period t
Air velocity in period t
Output power of photovoltaic equipment in period t
Light intensity in period t
Solar panel temperature in period t
Reference temperature in period t
External ambient temperature in period t
Solar radiation intensity in period t
Electric power consumed by internal cogeneration unit in period t
Thermal power consumed by the internal cogeneration unit in period t
Gas power consumed by the internal cogeneration unit in period t
Gas turbine generation power in period t
Flue gas waste heat power in period t
Natural gas consumption during in period t
Heat recovery power of the waste heat boiler in period t
Heat energy generated of the ground source heat pump in period t
Electric energy consumed of the ground source heat pump in period t
Cool power generated of EC in period t
Cool power generated of AC in period t
Electric energy consumed of the electric chiller in period t
Heat energy consumed of the absorption chiller in period t
Energy storage of energy storage equipment i in period t
Energy storage of equipment i in period t
Charging power and discharging power of energy storage equipment i in period t
Output power of equipment i in period t
Upper limit of the charging and discharging power of energy storage equipment i
Energy storage of energy storage equipment i in period t
Energy storage of equipment i in period t
Consumption rate of energy storage equipment i
Maintenance cost of equipment i
Rated capacity of equipment i
Upper and lower limits of the charging and discharging state of the energy storage equipment i
Membership function of a very small target and maximum target
ith objective function
Minimum and maximum of the ith objective function
Charging efficiency and discharging efficiency of energy storage equipment i
Maximum eigenvalue
Characteristic specific gravity of the ith target
Entropy of the ith target
Membership function of the ith objective
State variable
Equality constraint
Inequality constraint
Upper and lower bounds of inequality
Generalized inequality constraints
Generalized upper bounds
References
1
BaratiF.SeifiH.SepasianM. S.NateghiA.Shafie-khahM.CatalaoJ. P. S. (2015). Multi-period Integrated Framework of Generation, Transmission, and Natural Gas Grid Expansion Planning for Large-Scale Systems. IEEE Trans. Power Syst.30 (5), 2527–2537. 10.1109/TPWRS.2014.2365705
2
CleggS.MancarellaP. (2016). Integrated Electrical and Gas Network Flexibility Assessment in Low-Carbon Multi-Energy Systems. IEEE Trans. Sustain. Energ.7 (2), 718–731. 10.1109/TSTE.2015.2497329
3
FanD.DouX.XuY.WuC.XueG.ShaoY. (2021). A Dynamic Multi-Stage Planning Method for Integrated Energy Systems Considering Development Stages. Front. Energ. Res.9, 723702. 10.3389/fenrg.2021.723702
4
GhoshS.KamalasadanS. (2017). An Energy Function-Based Optimal Control Strategy for Output Stabilization of Integrated DFIG-Flywheel Energy Storage System. IEEE Trans. Smart Grid.8 (4), 1922–1931. 10.1109/TSG.2015.2510866
5
HuQ.LiangY.DingH.QuanX.WangQ.BaiL. (2022). Topological Partition Based Multi-Energy Flow Calculation Method for Complex Integrated Energy Systems. Energy.244, 123152. 10.1016/j.energy.2022.123152
6
HuangW.ZhangN.YangJ.WangY.KangC. (2019). Optimal Configuration Planning of Multi-Energy Systems Considering Distributed Renewable Energy. IEEE Trans. Smart Grid.10 (2), 1452–1464. 10.1109/TSG.2017.2767860
7
KlyapovskiyS.YouS.CaiH.BindnerH. W. (2019). Integrated Planning of a Large-Scale Heat Pump in View of Heat and Power Networks. IEEE Trans. Ind. Applicat.55 (1), 5–15. 10.1109/TIA.2018.2864114
8
KoltsaklisN. E.KnápekJ. (2021). Optimal Scheduling of a Multi-Energy Microgrid. Chem. Eng. Trans.88, 901–906. 10.3303/CET2188150
9
LiC.YangH.ShahidehpourM.XuZ.ZhouB.CaoY.et al (2020). Optimal Planning of Islanded Integrated Energy System with Solar-Biogas Energy Supply. IEEE Trans. Sustain. Energ.11 (4), 2437–2448. 10.1109/TSTE.2019.2958562
10
LiuN.WangJ.WangL. (2019). Hybrid Energy Sharing for Multiple Microgrids in an Integrated Heat-Electricity Energy System. IEEE Trans. Sustain. Energ.10 (3), 1139–1151. 10.1109/TSTE.2018.2861986
11
LuoZ.GuW.WuZ.WangZ.TangY. (2018). A Robust Optimization Method for Energy Management of CCHP Microgrid. J. Mod. Power Syst. Clean. Energ.6 (1), 132–144. 10.1007/s40565-017-0290-3
12
MaT.WuJ.HaoL.LeeW.-J.YanH.LiD. (2018). The Optimal Structure Planning and Energy Management Strategies of Smart Multi Energy Systems. Energy.160, 122–141. 10.1016/j.energy.2018.06.198
13
NicolosiF. F.AlberizziJ. C.CaligiuriC.RenziM. (2021). Unit Commitment Optimization of a Micro-grid with a MILP Algorithm: Role of the Emissions, Bio-Fuels and Power Generation Technology. Energ. Rep.7, 8639–8651. 10.1016/j.egyr.2021.04.020
14
QuM.DingT.JiaW.ZhuS.YangY.BlaabjergF. (2021). Distributed Optimal Control of Energy Hubs for Micro-integrated Energy Systems. IEEE Trans. Syst. Man. Cybern, Syst.51 (4), 2145–2158. 10.1109/TSMC.2020.3012113
15
ShenF.ZhaoL.DuW.ZhongW.QianF. (2020). Large-scale Industrial Energy Systems Optimization under Uncertainty: A Data-Driven Robust Optimization Approach. Appl. Energ.259, 114199. 10.1016/j.apenergy.2019.114199
16
SirventM.KanelakisN.GeisslerB.BiskasP. (2017). Linearized Model for Optimization of Coupled Electricity and Natural Gas Systems. J. Mod. Power Syst. Clean. Energ.5 (3), 364–374. 10.1007/s40565-017-0275-2
17
SunY.ZhangB.GeL.SidrovD.WangJ.XuZ. (2020). Day-ahead Optimization Schedule for Gas-Electric Integrated Energy System Based on Second-Order Cone Programming. Csee Jpes.6 (1), 142–151. 10.17775/CSEEJPES.2019.00860
18
WangD.ZhiY.ZhiY.YuB.ChenZ.AnQ.et al (2018). Optimal Coordination Control Strategy of Hybrid Energy Storage Systems for Tie-Line Smoothing Services in Integrated Community Energy Systems. Csee Jpes.4 (4), 408–416. 10.17775/CSEEJPES.2017.01050
19
WangY.ZhangN.ZhuoZ.KangC.KirschenD. (2018). Mixed-integer Linear Programming-Based Optimal Configuration Planning for Energy Hub: Starting from Scratch. Appl. Energ.210, 1141–1150. 10.1016/j.apenergy.2017.08.114
20
WangL.HouC.YeB.WangX.YinC.CongH. (2021). Optimal Operation Analysis of Integrated Community Energy System Considering the Uncertainty of Demand Response. IEEE Trans. Power Syst.36 (4), 3681–3691. 10.1109/TPWRS.2021.3051720
21
WangW.HuangS.ZhangG.LiuJ.ChenZ. (2021). Optimal Operation of an Integrated Electricity-Heat Energy System Considering Flexible Resources Dispatch for Renewable Integration. J. Mod. Power Syst. Cle.9 (4), 699–710. 10.35833/MPCE.2020.000917
22
WangY.LiR.DongH.MaY.YangJ.ZhangF.et al (2019). Capacity Planning and Optimization of Business Park-Level Integrated Energy System Based on Investment Constraints. Energy.189, 116345. 10.1016/j.energy.2019.116345
23
WuJ.LiB.ChenJ.DingY.LouQ.XingX.et al (2021). Multi-objective Optimal Scheduling of Offshore Micro Integrated Energy System Considering Natural Gas Emission. Int. J. Electr. Power Energ. Syst.125, 106535. 10.1016/j.ijepes.2020.106535
24
XiaoH.PeiW.PeiW.DongZ.KongL. (2018). Bi-level Planning for Integrated Energy Systems Incorporating Demand Response and Energy Storage under Uncertain Environments Using Novel Metamodel. Csee Jpes.4 (2), 155–167. 10.17775/CSEEJPES.2017.01260
25
YangJ.ZhangN.KangC.XiaQ. (2018). Effect of Natural Gas Flow Dynamics in Robust Generation Scheduling under Wind Uncertainty. IEEE Trans. Power Syst.33 (2), 2087–2097. 10.1109/TPWRS.2017.2733222
26
ZhangS.GuW.YaoS.LuS.ZhouS.WuZ. (2021). Partitional Decoupling Method for Fast Calculation of Energy Flow in a Large-Scale Heat and Electricity Integrated Energy System. IEEE Trans. Sustain. Energ.12 (1), 501–513. 10.1109/TSTE.2020.3008189
27
ZhaoB.ConejoA. J.SioshansiR. (2018). Coordinated Expansion Planning of Natural Gas and Electric Power Systems. IEEE Trans. Power Syst.33 (3), 3064–3075. 10.1109/TPWRS.2017.2759198
28
ZhouZ.LiuP.LiZ.NiW. (2013). An Engineering Approach to the Optimal Design of Distributed Energy Systems in China. Appl. Therm. Eng.53 (2), 387–396. 10.1016/j.applthermaleng.2012.01.067
Summary
Keywords
multi-energy flow coupling system, multi-objective collaborative optimization, combined weighting method, simplified primal-dual interior point algorithm, carbon emission
Citation
Zong X, Yuan Y and Wu H (2022) Multi-Objective Optimization of Multi-Energy Flow Coupling System With Carbon Emission Target Oriented. Front. Energy Res. 10:877700. doi: 10.3389/fenrg.2022.877700
Received
17 February 2022
Accepted
28 March 2022
Published
10 May 2022
Volume
10 - 2022
Edited by
Qingxin Shi, North China Electric Power University, China
Reviewed by
Linquan Bai, University of North Carolina at Charlotte, United States
Nikolaos Koltsaklis, Czech Technical University in Prague, Czechia
Zhenkun Li, Shanghai University of Electric Power, China
Updates
Copyright
© 2022 Zong, Yuan and Wu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xuanjun Zong, zongxuanjun@163.com
This article was submitted to Smart Grids, a section of the journal Frontiers in Energy Research
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.