# Based on the difference of Newton’s method integrated energy system distributed collaborative optimization

^{1}Faculty of Electrical and Control Engineering, Liaoning Technical University, Huludao, China^{2}College of Mechanical and Electrical Engineering, Xinjiang Agricultural University, Urumqi, China^{3}The Electrical Engineering College, Guizhou University, Guiyang, China

With the integration of renewable energy into the grid, the traditional power system stability faced by huge challenges, and the development of integrated energy system, it is of essence to improve the coupling of multiple integrated energy systems of different types, management in the integrated energy system and reduce the pressure of communication and computing, in this paper, we construct a distributed Newton algorithm based on Newton’s method to accelerate the solving speed, which decreases the times of iterations to reduce the pressure of communication and calculation, saving the cost of operation. Besides, privacy protection is particularly important for a distributed control system, under the premise that calculation speed is guaranteed, meanwhile, privacy protection of all agents in an integrated energy system is also critical. This study uses annular directed distributed algorithm to enhance the privacy of integrated distributed energy systems in the intelligent body, so as to fully ensure the privacy safety of all agents in the system. Moreover, the forementioned difference Newton algorithm in this study avoid the behavior of Zeno, greatly accelerating the speed of iteration and finding the best energy market price,. At the same time, the privacy safety of all agentsin the distributed energy system are ensured. Finally, a distributed integrated energy system based on the algorithm proposed by this study has went through theoretical proof and simulation experiment, whose result shows the validity of the algorithm.

## 1 Introduction

Recent years have seen the need to bring a major shift of energy source from coal to and electricity in an attempt to ensure power supply Therefore, renewable energy is connected to the power system, which can ensure sustainable and reliable power supply. Despite unprecedented challenge and change, the traditional power system gradually transformed into new power system generating clean energy. However, as the new energy power systemincreasing in scale, it still faces numerous challenges, such as new balance system and complex security mechanism. The introduction of the integrated energy system can meet the regional energy demand and the multi-source strategy development is of great significance, that said, the key still lies in safeguarding regional energy safety and stable operation of energy system.

For the traditional electric power system, renewable energy is introduced into too little, therefore, most of the scholars at home and abroad research contains only a single power network (YANG and WANG, 2021), and for the centralized algorithm. These algorithms included multi-objective optimization scheduling (Zhang et al., 2020), mixed nonlinear programming (Marty et al., 2017), (Hemamalini and Simon, 2009) Newton method, and traditional iterative methods (Lin and Viviani, 1984) (Hemamalini and Simon, 2009) to solve the non-convex participant energy management problem in centralized energy systems. The literature suggests that centralized algorithms can more accurately obtain optimal values (Hemamalini and Simon, 2009), and are fast and easy to design. However, using centralized algorithms requires significant computational resources and communication costs, and damage to the central agent can be difficult to recover, such as (Yin et al., 2018). To address the aforementioned challenges, a distributed algorithm has been proposed by relevant scholars, which effectively overcomes the issue of handling large volumes of data and prevents problems such as information processing. The distributed algorithm encompasses the following aspects: Chow initially proposed a numerical method for consistency that was utilized by foreign scholars to solve the problem of distributed energy scheduling (Zhang and Chow, 2012). For a single distribution power grid, the primary consideration is the impact of electricity prices on power consumption (Xie et al., 2022a). Moreover, the robustness and control of the micro grid must be taken into account (Xie et al., 2022b), while optimization of parallel distribution in weak power grids is achieved through the application of the method of group economical ICA and the NSGA-II (Nie et al., 2023), (Zhong et al., 2022), which aims to identify the optimal operating point for the micro grid. The optimization of the micro grid mainly includes approaches such as neural network-based methods (Zhang et al., 2023), alternating direction multiplier methods (Gao et al., 2022) (Zhu et al., 2022), and dynamic programming (Yang and Yang, 2022). Scholars, such as Yang Ping, have proposed the GPS model in light of the relationship between information flow and energy flow, and have also developed optimal control strategies to explore the micro grid (Yang and Yang, 2022). The aforementioned research effectively addresses the challenges associated with communication and computing stress in traditional micro grid control centers. In addition, due to the use of distributed methods, the absence of a control center ensures that local damage has minimal impact on the entire system. However, the research only focuses on electricity and does not consider other forms of energy. With the constant improvement of the proportion of new energy power system, balance system, the security mechanism of the new type of power system problems such as challenged, in this case, the collaborative optimization of a variety of energy for the stability and security of the power system is obviously much better than a single grid. In this context, some scholars have proposed the concept of an integrated energy system that differs from single energy networks. Considering multiple energy networks together can enhance the effect of energy optimization, but it also increases the complexity of the transformation mechanisms across various forms of energy, such as electricity, gas, and heat, in both time and space scales (Lv, 2022). Consequently, the comprehensive optimization of an integrated energy system is more challenging than optimizing a single grid (Schfer et al., 2018). To address these issues, scholars have developed distributed computing methods for finding optimal values in integrated energy systems. Many studies by domestic and foreign scholars have proposed distributed non-iterative algorithms for multi-agent coordination optimization problems (Tan et al., 2019; Tan et al., 2021). For instance, one study (Munshi and Mohamed, 2019) proposed an unsupervised algorithm that uses electric meter data to determine electrical load parameters, while another study (Zhang et al., 2017) introduced the mixed alternating direction multiplier method to solve the coupling relationship between various forms of energy in the integrated energy system. Scholars such as Zhang have applied the integrated energy unit and alternating direction multiplier method to address the multivariate coupling between electric and heat energy systems and solve coordination problems between the energy networks (Zhang et al., 2017). On the basis of previous research, scholars both at home and abroad have mainly focused on the optimization of distributed optimization methods in integrated energy systems with regards to convergence speed. This has been done by addressing issues such as information attacks (Zhao et al., 2016; Duan and Chow, 2019), non-convexity of energy systems (Chen and Zhao, 2019; Huang et al., 2019; Li et al., 2019), and so on.

Although about the optimal management of traditional integrated energy system research have been studied from many aspects, can meet the balance between supply and demand, energy scheduling and achievements, but the traditional studies in distributed energy system, between the iterative speed and privacy protection problems still exist, Literature (Zhang et al., 2017; Li et al., 2020) has shown that when the iteration speed is slow, it is possible to ignore the kink behavior, which can lead to a loss in wireless loop. The traditional energy management approach is thus faced with the challenge of being too slow for a long period of time, resulting in delayed energy supply during emergencies and inevitable loss. Pressure is too large and integrated energy systems of communication, iterative speed too slow, resulting in the high cost problem. Based on the above analysis, this paper combined with the ring to algorithm design difference Newton’s method, Newton method and greatly accelerate the iteration speed, reduce the number of iterations, delay to solve energy, reduce communication and computing pressure, privacy protection and integrated energy systems between adjacent agent. The main contributions of this paper are as follows.

1) The difference Newton’s method is proposed based on the Newton’s method. In comparison with the traditional approach, this method is characterized by a faster rate of iteration and a reduced number of iterations. Consequently, in the context of an integrated energy system, it can significantly reduce communication and computational overhead and mitigate energy latency.

2) The ring-based distributed algorithm is employed in the integrated energy system to ensure privacy protection among its participants. Conventional algorithms rely on participants to protect their privacy, but they require the exchange of adjacent information such as energy input/output power and energy price. In this regard, the use of the ring-based distributed algorithm effectively addresses this issue by providing complete privacy protection for all participants.

3) To avoid the divergence of Newton’s method during the downhill, a kino line may occur where the algorithm iterates indefinitely within a limited time, leading to an infinite loop. To prevent this behavior, the proposed algorithm in this paper is designed to avoid kino behavior.

The paper is structured as follows: In the first section, the integrated energy system established in this study is introduced, along with the function and constraint conditions for each participant’s cost in the system. The second section describes the Difference Newton’s method employed in this paper (DNEA) and how it addresses the energy coupling issue. The third section provides a proof of DNEA’s speed, convergence, optimality, and avoidance of kino behavior. The fourth section presents a simulation validation of the proposed DNEA algorithm in the established integrated energy system. Finally, the fifth section offers concluding remarks.

## 2 Models of IES

The composition of body, specifically its internal structure, comprises various energy devices that serve distinct functions. 1)The power generation device includes distributed coal-fired generators, gas generators, solar generators, wind turbines, and energy storage systems. 2) The heating device comprises distributed coal-fired heat production devices, distributed combustion gas heat production equipment, photovoltaic production engines, and distributed storage devices. 3) The thermal electric power plant machine. 4) The distributed gas suppliers constitute the remaining energy equipment. Each of these energy bodies accommodates three energy types, namely, electricity, gas, and heat, which need to be considered while accounting for electricity price fluctuations, load-side random scheduling, and demand response for electricity, gas, and heat. Because in traditional power system load is not adjustable, so you need to increase the “adjustable load”, such as: conversion of electricity and heat, electricity and hydrogen between mutual conversion, energy storage device, etc.,. This paper is suitable for small energy network connected to the electricity grid if there are toning, can from the grid to compensate.

The figure illustrates the encryption of data to safeguard the privacy of each pluripotent micro power grid in the integrated energy system. During each iteration process, the control center establishes encrypted data to ensure secure communication.

### 2.1 Models of renewable energy devices

This article considers five energy bodies and takes into account the changing demands and environmental factors affecting the energy efficiency of each equipment within these bodies. Based on this, the following constraints have been established.

1) Distributed energy physical quantity balance constraints

Type: *i*th a total electricity energy body; *i*th coal-fired generator set, energy storage device, gas supplier collection, cogeneration unit set, pv electricity production device set, load, and fan electricity production device under given energy body scenarios.

For each CLP energy body, the production amount of equipment for gas and heat must remain equal to the difference between the total load and the exchange of energy on the load side during operation. In addition, the exchange of distributed energy should satisfy the following constraints:

Regarding the heat transfer in type 2) variables, the concrete form of electricity exchange is as follows: the same redundancy avoidance applies, and it is not presented here. 3), *i*th energy body;

2) Considering the randomness, volatility, and renewable energy output fluctuation value limit conditions, the following expressions apply:

Type:

As a renewable energy source, solar electricity is widely utilized currently; however, its volatility and randomness necessitate the inclusion of photovoltaic (pv) motor and solar heating into the following functions:

Type:

Fan electricity production cost function:

Type:

Note: this article does not consider the relationships between kW and heat, electricity, and gas; therefore, the comprehensive energy kW should be applied to all energy units involved.

3) Considering coal heating device, coal thermal power plant and cogeneration plant output fluctuation value limit condition are:

Type:

4) Coal-fired thermal power equipment, coal heating equipment, and combined heat and power output ramping restrictions are as follows:

Type: *T* moment coal-fired thermal power equipment compared to T-1 h, the maximum change in *T* moment coal-fired heat production equipment compared to *T-1* h, and the maximum change in T moment cogeneration equipment compared to *T-1* h for electricity and heat production, respectively. *T-1* coal-fired thermal power equipment, *T-1* coal-fired heat production equipment, and *T-1* cogeneration plant for electricity and heat, respectively.

Coal is commonly used in energy bodies to produce electricity and heat, and its consumption characteristics are determined by the energy size produced at time T. Although its stability is reliable, it causes a certain level of environmental pollution and is subject to ramping constraints. The cost functions for coal-fired generators and heat production engines are as follows:

Type:

As an energy body, the CLP cogeneration unit exhibits a thermal coupling relationship with its main equipment. The cost function is determined by the simultaneous production of electric and thermal energy, and its output falls within a specific range. The cogeneration unit cost function is as follows:

Type:

5) Distributed gas supplier supply constraints are given by the following expression:

Type:

6) Resistance and energy storage restrictions for energy storage devices and heat values are as follows:

Type:

7) Electricity, heat, and energy storage device discharge and power change constraints are as follows:

Type:

Energy storage devices store energy when electricity prices are low and release energy when prices are high, playing a regulatory role. Consequently, they are essential equipment within energy bodies. Due to varying factors across different types of energy storage devices, a unified storage device cost function is established as follows:

Type:

8) Gas-to-electricity transfer and thermal conversion rates are as follows:

Type: *a* and *b* represent the gas-to-electricity and heat conversion rates, respectively.

9) Lateral load consumption of electricity and heat cost functions are as follows:

Type:

### 2.2 The interests of the function

The benefit of the energy body function revenue function and cost function are two parts, and the mathematical expression is as follows:

Type: *i*th a energy body in T time overall interests, *i*th total earnings, *i*th the total cost. The specific expressions for the profit function and cost function in Eq. 34 are as follows:

Type:

### 2.3 The objective function

This study focuses on the optimization of an integrated energy system that aims to coordinate the use of different types of energy, such as electricity, gas, and heat, with the goal of reducing production costs and improving energy efficiency. Specifically, we investigate the collaborative optimization of three energy supply and demand types, namely, electricity, heat, and gas. The energy conversion model is employed to allow for price adjustments among the three types of energy (Yu-Shuai et al., 2020). We strive to achieve a globally optimal solution that takes full advantage of the complementary characteristics of various energy sources to obtain the most economic price, while ensuring the balance between energy supply and demand. The optimization problem is formulated as (34), using a modified version of Newton’s method, known as the difference Newton’s method, to speed up the algorithm iteration and reduce communication and computational pressure. The objective function is expressed as type (40), which represents the overall benefits of the entire energy system, while ensuring that the net energy value is zero, thus guaranteeing the balance between energy supply and demand. The function expression is presented below.

## 3 The algorithm design

### 3.1 The traditional Newton iteration method

The traditional numerical method referred to as Newton’s method is commonly utilized to solve nonlinear equations. Within integrated energy systems, collaborative optimization can utilize the traditional Newton’s method to solve complex optimization models. The basic idea of this method is to iteratively solve equations of zero and determine the direction of the next iteration through the first-order approximation of the equation. In practical applications, the traditional Newton method can achieve optimal cooperation among parties based on the integrated energy system states and parameters of continuous optimization. When addressing comprehensive power, heat, and gas optimization problems in energy systems, the traditional Newton method can be applied to solve conflicts and coordinate energy sources to realize optimal system synergy (Tan and Li, 2022).

### 3.2 DNEA algorithm design

In integrated energy systems, each energy entity is equipped with its own processor. Nonetheless, the absence of a centralized system causes a distributed structure which lacks significant communication and computing power. Consequently, the computing capability of each processor in the energy entity is limited. When faced with the high-speed and high utilization of renewable energy, the conventional algorithm struggles to compute the parameters in a distributed energy system. To tackle this challenge, this study introduces the difference Newton’s method as a solution. This method is designed to resolve issues with slow calculation speed and communication and computing pressure. By applying the difference Newton’s method to every energy entity during the calculation process, the computational speed is significantly enhanced, while maintaining a balance between supply and demand and maximizing profit. The following outlines the design of the difference Newton’s method in the calculation process of the integrated energy system.

Due to the convex nature of the cost function designed in this paper, the derivative of the energy body’s benefit function for all participants is the iterative price

Type: Given the goal of finding the optimal balance between supply and demand conditions for the maximum overall energy interests, the price of electricity and heat is assumed to be the same. Hence, the price of electricity and heat are represented by

Type:

The difference between the value of

Assuming that the second energy iteration price is

Type: according to the relationship between price and the imbalance between supply and demand, if all functions are linear and the balance between supply and demand is guaranteed, the energy price is *i*th iteration imbalance between supply and demand, and

### 3.3 Ring signature algorithm design

This paper ensures privacy among all participants in the integrated energy system by initializing the encrypted energy data as

Type: In the first i-1 iterations, energy body information *i*th energy body, and *i*th energy body. The information transfer direction is shown in Chart 1. This process strengthens privacy protection among energy bodies, ensuring each energy body's privacy and fully protecting the privacy of all agents in the system. This is an improvement compared to traditional privacy protection, which does not adequately protect privacy between integrated energy system energy bodies.

## 4 DNEA algorithm theory to prove

### 4.1 DNEA iteration speed

Setting parameters *R*^{+} denotes all positive values. The traditional iterative method using Newton’s method establishes an adjustment volume of each time as

Type: The total number N for iteration,

### 4.2 DNEA convergence is proved

Starting with the first iteration,

### 4.3 DNEA optimality

In the energy body, a balance between supply and demand for electricity and heat must be maintained. Electricity and heat conversion are employed to compensate for energy deficiencies. Assuming equal prices for electricity and heat, if the energy body lacks a heat quantity of

### 4.4 DNEA kinetic behavior analysis

This section provides mathematical proof that the difference Newton iterative algorithm designed in this paper effectively avoids kinetic behavior. Based on the convergence proof, the difference Newton designed in this paper gradually tends toward a supply and demand balance of 0. Setting parameter

## 5 DNEA algorithm simulation analysis

To verify the algorithm presented in this manuscript, a testing of a distributed algorithm based on differential points and Newton’s method was performed on an energy system described in Appendix 5, wherein the fundamental parameters of the simulation device are exhibited. Figure 1 portrays the integrated energy system under consideration. A condition is stipulated in this paper, wherein the energy supply and demand must remain within 10 kW of each other to achieve a balance between energy production and consumption. The algorithm was designed to accomplish this goal, and it achieved the desired balance within five iterations. The simulation process is as follows.

### 5.1 DNEA simulation analysis iteration speed

A comparison between the iterative processes of the distributed algorithm based on differential points and Newton’s method and the traditional finite difference algorithm was conducted on five energy systems with unbalanced supply and demand. The simulation results, depicted in Figure 2, demonstrate that the differential points Newton algorithm achieved a supply and demand balance of less than 10 kW within the fifth iteration without causing any disturbances. In contrast, the traditional algorithm failed to reach the optimal balance even after 50 iterations. Dichotomy optimization for the iterative error is small, the eight time but also known from the analysis of the simulation, the final results, not avoid kino, cannot ensure the final iteration for optimal results. Simulation results depicted in Figure 2, Figure 3 indicate that all participants in the energy system, designed in this manuscript, achieved rapid convergence to the optimal forecast price point while ensuring a balance between supply and demand. These results verify the feasibility and effectiveness of the algorithm presented in this paper. The differential points Newton algorithm was employed to calculate the energy system under conditions of supply and demand balance, which yielded faster iteration speeds compared to traditional computing methods by several orders of magnitude. Consequently, energy losses, delays, and communication costs were considerably reduced, leading to more efficient energy systems.

### 5.2 DNEA simulation convergence analysis

In this paper, an integrated energy system has been designed. The simulation diagram presented in Figure 4, Figure 5, Figure 6 depicts the supply side, and reveals that the total electricity generated is 69828694.93 kW. The final values of heat and gas are 35815789.39 kW and 368579.09 kW, respectively. The relevant parameters, such as the change of price and the fast convergence to phase, contribute to achieving a stable state. As a result, the aforementioned values remain practically unchanged.

Furthermore, it is worth noting that Newton’s method, when used in this study, yields highly accurate and reliable calculation results. Moreover, it significantly reduces the communication and computation workload associated with the distributed energy system, facilitating a fast scheduling process and ensuring a steady state energy system. As shown in Figure 7, the power supply amounts to 72907191.36 kW, and the heating load reaches 32737316.53 kW, thereby achieving a balance between supply and demand.

Combining the simulation results presented in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, it can be concluded without surprise that by utilizing the difference Newton’s method, this study is capable of rapidly achieving a stable state, while simultaneously ensuring that the results remain stable post convergence.

### 5.3 Ring to privacy protection analysis

Traditional algorithm (Tan and Li, 2022) privacy often meet:

Type: *Mpl*_{T} refers to privacy, and *Ac*_{T}, represents average connectivity, the term *Na*_{T} denotes an agent. Combined with the simulation Figure 9 shows in conventional privacy protection algorithms, an increase in the number of agents and average connectivity leads to a gradual enhancement in privacy levels. However, in this study, the researchers have designed a distinct Newton algorithm with a ring structure. In this approach, the growth in the number of agents and average connectivity does not result in elevated privacy levels, but rather maintains comprehensive privacy protection for the agents involved.

### 5.4 Usage scenarios and limit analysis of the algorithm

The utilization of the DNEA algorithmic approach is not applicable within the context of an integrated energy system, as the cost function follows a convex pattern. This system is characterized by interconnections, and it should not be employed in scenarios with high-quality energy standards, such as parks, intelligent buildings, hospitals, among others.

## 6 Conclusion

The calculation speed for the traditional iterative algorithm was evaluated multiple times, and each participant achieved the optimal convergence value. This confirmed the feasibility and stability of the proposed algorithm. However, it should be noted that this algorithm is not suitable for participants in the energy system with convex cost functions and constraints. Future scholars may explore non-convex cost functions and constraints to further develop the proposed algorithm. One limitation of this paper is that it did not address the non-convex cost function and constraints for the proposed algorithm.

## 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

XC and XK has corrected the format and simulation prove. Other study is carried out by XL. All authors contributed to the article and approved the submitted version.

## 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.

## Supplementary material

The supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenrg.2023.1215786/full#supplementary-material

## References

Chen, G., and Zhao, Z. Y. (2019). Delay effects on consensus-based distributed economic dispatch algorithm in microgrid. *IEEE Trans. Power Syst.* 33 (1), 602–612. doi:10.1109/tpwrs.2017.2702179

Duan, J., and Chow, M. Y. (2019). A novel data integrity attack on consensus based distributed energy management algorithm using local information. *IEEE Trans. Industrial Inf.* 15 (3), 1544–1553. doi:10.1109/tii.2018.2851248

Gao, S., He, J., Yang, S., and Xiao, B. (2022). En ergy management of interconnected microgrid bas ed on alternating direction method of multipliers method [J]. *Grid Clean Energy* 38 (06), 1 13–120.

Hemamalini, S., and Simon, S. P. (2009). Maclaurin series-based Lagrangian method for economic dispatch with valve-point effect. *IET Generation, Transm. Distribution* 3 (9), 859–871. doi:10.1049/iet-gtd.2008.0499

Huang, B. N., Liu, L. N., Zhang, H. W., Li, Y. S., and Sun, Q. Y. (2019). Distributed optimal economic dispatch for microgrids considering communication delays. *IEEE Trans. Syst. Man, Cybern. Syst.* 49 (8), 1634–1642. doi:10.1109/tsmc.2019.2900722

Li, F. Y., Qin, J. H., and Kang, Y. (2019). Multi-agent system based distributed pattern search algorithm for non-convex economic load dispatch in smart grid. *IEEE Trans. Power Syst.* 34 (3), 2093–2102. doi:10.1109/tpwrs.2018.2889989

Li, Y., Gao, D. W., Gao, W., Zhang, H., and Zhou, J. (2020). Double-mode energy management for multi-energy system via distributed dynamic event-triggered Newton-raphson algorithm. *IEEE Trans. Smart Grid* 11 (6), 5339–5356. doi:10.1109/tsg.2020.3005179

Lin, C. E., and Viviani, G. L. (1984). Hierarchical economic dispatch for piece wise quadratic cost functions. *IEEE Trans. Power Apparatus Syst.* 6, 1170–1175. doi:10.1109/tpas.1984.318445

Lv, H. (2022). *Research on coordination method of multi-type demand response considering electrica l thermal coupling [D]*. Jinan, China: Shandong University.

Marty, F., Serra, S., Sochard, S., and Reneaume, J. M. (2017). Economic optimization of a combined heat and power plant: Heat vs electricity. *Energy Procedia* 116, 138–151. doi:10.1016/j.egypro.2017.05.062

Munshi, A. A., and Mohamed, Y. A. R. I. (2019). Unsupervised nonintrusive extraction of electrical vehicle charging load patterns. *IEEE Trans. Industrial Inf.* 15 (1), 266–279. doi:10.1109/tii.2018.2806936

Nie, Y., Peng, C., Hu, Y., He, Y., Ma, G., and Huang, C. (2023). Parallel distribute d optimal economic dispatch of high penetration microgrid based on edge computing [J/OL]. *So uthern Power Grid Technol.* 1-11. [2023-01-20]].

Schfer, B., Beck, C., Aihara, K., Witthaut, D., and Timme, M. (2018). Non-Gaussian power grid frequency fluctuations characterized by Lévy-stable laws and superstatistics. *Nat. Energy* 3, 119–126. doi:10.1038/s41560-017-0058-z

Tan, J., and Li, H. (2022). Annular directed distrib uted algorithm for energy internet. *Int. Trans. Electr. Energy Syst.* 2022. doi:10.1155/2022/7717605

Tan, Z., Zhong, H., Wang, J., Xia, Q., and Kang, C. (2019). Enforcing intra-regional constraints in tie-line scheduling: A projection based framework. *IEEE Trans. Power Syst.* 34 (6), 4751–4761. doi:10.1109/tpwrs.2019.2913876

Tan, Z., Zhong, H., Xia, Q., and Kang, C. (2021). Non-iterative multiarea coordinated dispatch via condensed system representation. *IEEE Trans. Power Syst.* 36, 1594–1604. doi:10.1109/tpwrs.2020.3019828

Xie, Z., Lu, C., Wang, X., Wen-tao, Y., and Qi, C. (2022a). Microgrid based on distribute d resource cooperative contro research and appl ication of bidding strategy [J]. *Electr. Age* 10, 37–42.

Xie, Z., Lu, C., wang, X., Wen-tao, Y., and Qi, C. (2022b). Microgrid based on distribute d resource cooperative contro research and appl ication of bidding strategy [J], the electrical age, 2022(10):37-42 control in dcmicrogrid based on linear matrix I nequalities [J], header-icon. *J. Sol. Energy* 43 (05), 45–52.

Yang, M., and Wang, J. (2021). Optimal scheduling of islanded microgrid considering uncertain output of renewable energy[J]. *Chin. J. Electri cal Eng.* 41 (03), 973–985.

Yang, Y., and Yang, P. (2022). The hierarchical modeling approach for centralized control microgrid cyber physical system [J]. *Chin. J. Electr. Eng.* 42 (19), 7088–7102.

Yin, S., Qian, A., Zeng, S. Q., Wu, Q., Ran, H., and Jiang, D. (2018). Challenges and prospects of multi-energy distributed optimization for energy internet. *Power Syst. Technol.* 42 (5), 1–10.

Yu-Shuai, L., Tian-yi, L., Gao, W., and Gao, W. (2020). Distributed collaborative optimization operati on approach for integrated energy system based on asynchronous and dynamic event-triggering C ommunication strategy [J]. *J. Automation* 46 (09), 1831–1843.

Zhang, H., Li, Y., Gao, D. W., and Zhou, J. (2017). Distributed optimal energy management for energy internet. *IEEE Trans. Industrial Inf.* 13 (6), 3081–3097. doi:10.1109/tii.2017.2714199

Zhang, N., Sun, Q., and Yang, L. (2020). A two-stage multi-objective optimal scheduling in the integrated energy system with we energy modeling. *Energy* 215, 119121. doi:10.1016/j.energy.2020.119121

Zhang, Y., Lin, Y., Huang, G., Yang, X., Weng, G., and Zhou, Z. (2023). Rev iew on applications of deep reinforcement learning in regulation of microgrid systems [J/OL]. *Power Grid Technol.* 1-15. [2023-01-20].

Zhang, Z., and Chow, M. Y. (2012). Convergence analysis of the incremental cost consensus algorithm under different communication net work topologies in a smart grid. *IEEE Trans. Power Syst.* 27 (4), 1761–1768. doi:10.1109/tpwrs.2012.2188912

Zhao, C. C., He, J. P., Cheng, P., and Chen, J. M. (2016). Analysis of consensus based distributed economic dispatch under stealthy attacks. *IEEE Trans. Industrial Electron.* 64 (6), 5107–5117. doi:10.1109/tie.2016.2638400

Zhong, Y., Zhou, Y., Gao, Z., Ton, J., and Yu, D. (2022). Optimal scheduling of micro grid based on ICA and NSGA-II compound algori thm [J]. *Appl. Electr.* 41 (11), 78–8 4.

Keywords: difference Newton method, integrated energy system, energy scheduling, optimization, distributed algorithm

Citation: Liu X, Chen X and Ke X (2023) Based on the difference of Newton’s method integrated energy system distributed collaborative optimization. *Front. Energy Res.* 11:1215786. doi: 10.3389/fenrg.2023.1215786

Received: 02 May 2023; Accepted: 31 May 2023;

Published: 13 June 2023.

Edited by:

Yushuai Li, University of Oslo, NorwayReviewed by:

Xiaofeng Li, Anhui University, ChinaBingyu Wang, North China Electric Power University, China

Copyright © 2023 Liu, Chen and Ke. 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: Xinyu Ke, 848618850@qq.com