ORIGINAL RESEARCH article

Front. Energy Res., 22 September 2025

Sec. Sustainable Energy Systems

Volume 13 - 2025 | https://doi.org/10.3389/fenrg.2025.1633719

Two stage coordination planning method of wind power and storage considering uncertainty of distributed source-load

  • 1. Electric Power Research institute of Yunnan Electric Power Grid Co.Ltd., Kunming, China

  • 2. Hubei Engineering and Technology Research Center for AC/DC Intelligent Distribution Network, School of Electrical Engineering and Automation, Wuhan University, Wuhan, China

Abstract

Introduction:

With the widespread integration of distributed power sources, the power grid is facing challenges such as increased losses, rising costs, voltage fluctuations, and overload, resulting in greater operational complexity. Traditional scheduling methods are no longer adequate, making reasonable planning of distributed power generation and energy storage configurations particularly crucial.

Methods:

This article proposes a two-stage wind-storage coordination planning method that considers source-load uncertainty. The approach is based on an improved antlion algorithm and incorporates distributed energy storage charging and discharging strategies. The first stage focuses on wind power site selection and capacity determination, using voltage offset, network loss, and comprehensive system cost as evaluation indicators. A multi-objective function model is established to balance grid stability and economic efficiency. The second stage introduces distributed energy storage devices to reduce power fluctuations while minimizing the sum of operation, maintenance, and storage investment costs, thereby optimizing the energy storage charging and discharging strategy. The improved antlion algorithm, enhanced with adaptive Lévy flight and golden sine theory, is used to solve the two-stage planning model.

Results:

The proposed method effectively improved system-level voltage distribution, reduced network losses, and lowered overall system costs. Specifically, it achieved a 27.95% increase in total capacity, a reduction of 32.14 kW in active power loss, and a total cost decrease of 221,200 yuan. The improved antlion algorithm demonstrated strong search capability, fast convergence speed, and high computational accuracy.

Discussion:

The results indicate that the proposed method is better aligned with practical requirements compared to traditional approaches. The improvements in system performance and cost efficiency highlight the effectiveness of the two-stage planning framework and the enhanced optimization algorithm. The method offers a viable solution for the integrated planning of wind power and energy storage systems under uncertainty.

1 Introduction

The integration of distributed power sources injects new voltage power into the distribution network, and the network topology and power flow distribution will also change accordingly. Unreasonable integration may result in problems such as reverse transmission of branch power flow, voltage exceeding limits, and increased line losses, affecting system operation (). Meanwhile, in distributed power generation, wind and photovoltaic power generation, as the main distributed energy sources, have the advantages of being renewable and environmentally friendly. However, their output power is unstable due to changes in wind speed and light intensity, which may lead to insufficient power supply or resource waste. Therefore, optimizing the configuration of distributed power sources and utilizing energy storage technology to mitigate their adverse effects on the power grid is crucial.

In terms of distributed power generation planning models, Chu and Qiao considered the output efficiency and load rate of distributed power generation units. They formulated a planning model with the objective of minimizing the comprehensive operational cost of the distribution network. Huang et al. calculated power flow and network losses using Monte Carlo sampling and applied a genetic algorithm to optimize costs, network losses, and surplus electricity from distributed power sources. Cao et al. addressed the uncertainties associated with wind, solar, and load variations by employing Latin hypercube sampling combined with an improved synchronous substitution method to generate representative scenarios. The model was solved using an improved particle swarm optimization algorithm, aiming to minimize the annual comprehensive cost. Su et al. proposed a coordinated optimization strategy for wind power, solar power, load demand, and energy storage systems, focusing on determining the optimal power and capacity configuration of energy storage devices. Their objective function included distribution network investment costs, maintenance costs, power purchase costs, and reliability costs, which were optimized using the particle swarm optimization algorithm. Other researchers have also used variables such as network loss as objective functions for analysis and optimization. However, most of the aforementioned studies focus on single-objective optimization, which may overlook the complex interactions in system operations and deviate from practical engineering applications. To address this limitation, scholars both domestically and internationally have conducted further research into multi-objective optimization models. Mohammad et al. constructed a multi-objective function based on indicators such as network loss, voltage deviation, and short-circuit current, and solved it using optimization algorithms. Truong et al. introduced a quasi-adversarial chaotic symbiotic biological search algorithm to address multi-objective optimization problems. Banihashemi et al. developed a multi-objective optimization model aimed at minimizing voltage deviation, line loss, and operational costs, which was solved using an improved genetic algorithm. Li et al. applied the theory of chance-constrained programming and employed the non-dominated sorting genetic algorithm (NSGA) to optimize objectives including minimizing the operational risk of distribution networks and minimizing annual comprehensive costs.

Overall, many literature currently use simple deterministic models to model distributed power generation planning problems, without considering source load uncertainty or the impact of distributed energy storage (Zhenqi, 2021b; Paiva et al., 2017; ; Xu et al., 2017; Sivaram et al., 2019; ; Junyang et al., 2018; ). Therefore, this article will establish a more comprehensive distributed power generation planning model, taking into account the uncertainty of distributed power generation output and the integration of energy storage, to ensure the safety, reliability, and economy of the power system. Therefore, this article proposes a distributed wind storage coordination planning method that takes into account the uncertainty of source load. Firstly, a multi-objective model is established with the constraints of power flow, voltage, and power, aiming to minimize system network losses, voltage deviations, and overall system costs. Taking into account the uncertainty of source and load, a first stage distributed wind power coordination optimization strategy is proposed; Then, taking into account constraints such as power supply, energy storage capacity, and State of charging/discharge, combined with decision variables obtained in the previous stage, a model is established based on system economic indicators, and a second stage distributed energy storage planning method is proposed; Finally, the improved antlion algorithm with adaptive Levy flight and golden sine theory as improvement factors was used to solve the proposed two-stage wind storage coordination planning method. Through simulation verification, it was proved that the proposed method can effectively improve the system voltage distribution level, reduce network losses, and further reduce the overall system cost, bringing good stability and economy to the system operation.

2 Two stage wind storage coordination planning methods

This article considers distributed energy storage charging and discharging strategies and proposes a two-stage wind storage coordination planning method that takes into account source load uncertainty (; Weiguo et al., 2016; Zhenqi, 2021a; ). The specific model diagram is shown in Figure 1.

FIGURE 1

2.1 Distributed wind power planning model

2.1.1 Objective function

The connection of power supply to the distribution network can effectively improve the system voltage level and reduce network losses, but an unreasonable connection scheme can have a significant impact on the operation of the distribution network and disrupt the safe and reliable operation of the system. Therefore, in this section, a multi-objective system for distribution network operation is established based on three indicators of power grid stability and economy, namely, system network loss, node voltage deviation, and annual comprehensive cost, to coordinate and plan the integration of distributed power sources into the distribution network (Yurong et al., 2020). The specific model is as follows:

In the

Formula 1

,

indicates the active power loss of the system;

indicates the system voltage offset value;

indicates the comprehensive cost of the system;

indicates the maximum active power loss of the system;

indicates the maximum comprehensive cost of the system;

indicate the inertia weight coefficients of each indicator, determined by the Analytic Hierarchy Process, they are

,

。

  • 1. System active power loss

In the

Formula 2

,

indicates Line loss;

indicates the conductivity value between node i and node j; L indicates the number of system lines;

,

indicate the voltage values of nodes i and j;

indicates the voltage phase angle difference between nodes i and j (

Ping et al., 2018

).

  • 2. System voltage offset value

In the

Formula 3

,

indicates the voltage offset value;

indicates the actual voltage value of the node;

indicates the rated voltage of the system; N indicates the number of system nodes.

  • 3. Annual comprehensive cost of the system

In the Formula 4, N indicates the number of typical scenarios; indicates probability of occurrence of typical wind load scenarios; indicates operating costs in typical scenarios; indicates investment costs in typical scenarios; indicates Network loss cost in typical scenarios; indicates Environmental benefits and costs in typical scenarios; indicates under typical scenarios, the cost of purchasing electricity online from a large power company; indicates power generation subsidies in typical scenarios. Among them, the cost expressions are as follows:

  • a. Running cost:

In the

Formula 5

,

indicates operating costs of wind power;

indicates the unit cost coefficient for wind power operation is taken as 0.045 yuan per kW hour;

indicates Selected

i

th wind power generation power;

indicates Annual duration of wind power generation;

indicates A collection of nodes that can be connected to wind power (

Junqiang et al., 2016

).

  • b. Investment cost:

In the

Formula 6

,

indicates Wind power investment cost;

indicates investment coefficient for wind power, etc.,;

indicates unit capacity investment cost of wind power at 10,000 yuan per kW hour; The expression for

is as follows

Formula 7

, and in the formula, d indicates annual interest rate which is 0.8,

indicates The full life cycle of wind power is taken as 10 years:

  • c. Line loss cost:

In the

Formula 8

,

indicates unit length line loss cost, taken as 0.6 yuan per kW hour;

indicates maximum network loss on line j;

indicates the maximum annual load loss time is taken as 4200 h.

  • d. Environmental benefit cost:

In the

Formula 9

n indicates the number of environmental pollutants,

indicates Unit cost of pollutant m;

indicates under the traditional power generation mode, the emission indicators for pollutant m are as follows: sulfur dioxide emission indicator is 4.5 g per kW hour, with a cost of 7.3 yuan per kilogram; nitrogen dioxide emission indicator is 1.64 g per kW hour, with a cost of 10 yuan per kilogram; carbon dioxide emission indicator is 90 g per kW hour, with a cost of 0.8 yuan per kilogram;

indicates total load of the line;

indicates Network loss; all emission indicators are based on typical data of the Chinese power system (

;

Ministry of Ecology and Environment, 2022

;

Jinnan et al., 2019

).

  • e. Electricity purchase cost:

In the

Formula 10

,

indicates the unit cost of purchasing electricity is 0.5 yuan per kW hour;

indicates maximum system load;

indicates the maximum annual load utilization time is taken as 4,200 h.

  • f. Power generation subsidy:

In the Formula 11, indicates The unit government subsidy fee is 0.4 yuan per kW hour.

2.1.2 Constraint condition

The main constraints considered by the distributed wind power planning model are as follows (

Kaiyuan, 2023

):

  • 1. Trend constraints

In the

Formula 12

,

,

indicate separately active and reactive power of node

;

,

indicate separately actual voltage of node j;

,

indicate separately conductance and susceptance of nodes and branches between them;

indicates voltage phase angle difference between nodes; N is the number of system nodes.

  • 2. Node voltage constraint

In the

Formula 13

,

indicates the voltage amplitude of node

;

,

indicates separately the upper and lower limits of the voltage at node

.

  • 3. Branch power constraint

In the

Formula 14

,

indicates the active power of the line between nodes

and

;

indicates separately the upper and lower limits of active power of the line connecting nodes

and

.

  • 4. Wind turbine capacity constraint

In the Formula 15, indicates maximum capacity of wind turbines; indicates maximum wind power capacity; indicates minimum capacity of wind power.

2.2 Distributed energy storage planning model

2.2.1 Objective function

In order to further optimize the system operation, this section introduces energy storage devices with peak shaving, valley filling, and flat wave suppression effects (

Mohammad, 2014

;

Truong et al., 2020

;

;

Ke et al., 2017

). The reasonable introduction of it greatly improves the stability and performance of the system. However, the cost of energy storage devices is high, and a large amount of investment can also increase the economic operating costs of the system, resulting in resource losses. Therefore, this section focuses on the balance between energy storage devices and power supply and demand, considering energy storage charging and discharging strategies and the entire life cycle of the devices. With the goal of minimizing energy storage investment and operation costs, the optimal energy storage device charging and discharging strategy is obtained. The objective function is as follows (

Formula 16

) (

):

  • 1. Energy storage operation and maintenance costs

In the

Formula 17

,

indicates operating costs of energy storage devices;

indicates the unit cost coefficient for energy storage operation is set at 0.05 yuan per kW hour;

,

indicate separately Energy storage charging power, energy storage discharging power;

,

Indicates the charging and discharging efficiency of the energy storage device, taken as 0.9;

,

indicates the charging and discharging state, with a value of 0 or 1, where 1 represents charging and 0 represents discharging.

  • 2. Energy storage investment cost

In the Formula 18, indicates investment cost of energy storage devices; indicates investment coefficient for energy storage devices, etc. , indicate separately unit power and capacity cost of energy storage, the energy storage power and capacity are 5,000 yuan per kilowatt or kW hour; indicates actual power of energy storage; indicates energy storage capacity, The expression for is as follows (Formula 19), where d indicates an annual interest rate of 0.8, indicates the full life cycle of energy storage and other devices is taken as 10 years:

2.2.2 Constraint condition

For the above objective function, the constraints of this model include energy storage device charging and discharging power, capacity constraints, state constraints, as well as system power balance constraints (

Ibrahim et al., 2008

).

  • 1. Energy storage charging and discharging power constraint

In the

Formula 20

indicate separately energy storage charging and discharging power;

indicate separately maximum charging and discharging power of energy storage;

indicate separately minimum charging and discharging power for energy storage.

  • 2. Energy storage capacity constraint

In the

Formula 21

,

indicates capacity of energy storage device;

indicates maximum capacity of energy storage device (

Zechun et al., 2017

).

  • 3. The constraints on the storage charging and discharging states are shown in Formula 22.

  • 4. System power balance constraint

In the Formula 23, indicates the power of the load.

3 Solution of two-stage wind storage coordination planning method based on improved antlion algorithm

3.1 Principle of antlion algorithm

The Ant Lion Optimization (ALO) algorithm, proposed by Mirjalili in 2015, is a metaheuristic optimization approach inspired by the hunting behavior of ant lions in nature (Jasim et al., 2023). This algorithm simulates several core processes: random walking, trap building, luring ants, capturing prey, and the implementation of an elite mechanism (). The key innovation of ALO lies in its adaptive boundary contraction strategy, wherein the radius of the ant lion trap decreases progressively with each iteration. This feature enables a smooth transition from global exploration to local exploitation and helps prevent premature convergence. The algorithm offers advantages such as fewer required parameters and a strong capacity for balance (). When applied to solve the two-stage wind-storage coordinated planning method presented in this paper, the specific solution procedure is as follows:

3.1.1 Ant random walk model

In the Formula 24, indicates ant random walk normalization route; indicates ants randomly walk along unnormalized routes; , indicate separately the -th ant randomly walks the minimum and maximum values in the -dimensional variable; , indicate the -th ant randomly walks to the minimum and maximum values after the -dimensional variable iteration.

3.1.2 Antlion trap model

Simulate the process using the roulette wheel selection mechanism and define it as:

In the Formula 25, indicates the position of the jth antlion in the -th iteration; , indicate separately the minimum and maximum values of all ants after the -th iteration.

3.1.3 Ant trapped in trap model

The trap range is defined as follows (Formula 26):

In the Formula 27, I indicates the size of the trap range; indicates the calibration coefficient is related to the number of iterations; indicates maximum number of iterations.

3.1.4 Elite antlion model

The antlion with the highest fitness during each iteration is called the elite antlion, In the iteration, the position of the ant is:

In the Formula 28, indicates the position of the ith ant in the -th iteration; , indicate separately the step size of ants randomly walking around the antlion and elite antlion in the -th iteration; (5) Antlions prey on ants and reconstruct new trap models

When the fitness of ants is higher than that of antlions, ants are captured. At this point, the antlion sets the capture location as the reconstruction trap location. The formula for this model is:

In the Formula 29, indicates the position of the -th ant in the -th iteration; indicates the position of the -th antlion in the -th iteration.

3.2 Improvement of antlion algorithm

Actual tests demonstrate that although the traditional ant lion algorithm employs a diversified search strategy, it suffers from limited local search capability and is prone to becoming trapped in local optima, thereby constraining improvements in solution quality. Furthermore, conventional algorithms exhibit inadequate convergence accuracy in high-precision application scenarios. To address these limitations, this study proposes an adaptive Levy flight mechanism. Leveraging its long-step-length jumping characteristic, this mechanism effectively enables individuals to escape local optimal regions, thereby enhancing the algorithm’s global exploration capability and mitigating premature convergence. Additionally, the golden sine theory is incorporated. Utilizing its refined search and rapid convergence properties, this approach facilitates more precise and efficient local exploitation within promising solution regions, significantly improving the algorithm’s convergence accuracy. By integrating an adaptive strategy that dynamically adjusts the balance between Levy flight-based exploration and golden sine-based exploitation, the proposed method intelligently coordinates global search and local development processes, ensuring high-precision convergence while effectively avoiding local optima (Jianfang et al., 2025).

3.2.1 Adaptive levy flight (Tanyildizi and Demir, 2017)

The random walk pattern of Levi’s flight follows the Levi distribution shown in Equation 30:

In the Formula 30, a represents the peak height of the Levi distribution. When a is a non integer positive real number, the position update is performed using the method shown in Formula 31 below:

In the Formula 31, , indicate the position of the -th ant in the -th and +1st iterations; indicates a random search path that follows the Levy distribution. conforms to the following constraints:

In the Formula 32, indicates levy parameter, usually setting at 1.5; μ and ν indicate random numbers that follow a standard normal distribution; The expression for is as follows (Formula 33):

In this section, an adaptive Levy flight mechanism is introduced to mutate ants in their position updates, making their updated positions more diverse. The introduction of this improvement factor makes the search scope more comprehensive, the overall optimization efficiency of the population higher, and effectively avoids the solution results from falling into local optima.

3.2.2 Golden sine algorithm ()

The position update process of the Golden Sine Algorithm is shown in the following Formula 34:

During the position update process, individual positions are randomly generated first, and is used to represent the position of the -th particle in the Tth iteration in -dimensional space, among them , Use to represent the historical best position of the -th particle in the Tth iteration. and are random number, , . and are coefficients obtained through the golden ratios, and golden ratios are , , .

The steps of the improved antlion algorithm based on adaptive Levy flight and golden sine algorithm are shown in Figure 2:

FIGURE 2

4 Example analysis

4.1 Validation of the superiority of the improved antlion algorithm

Based on the objective function and corresponding constraints mentioned in Section 2, this section adopts the antlion algorithm, improved antlion algorithm, and particle swarm algorithm for optimization and solution. The number of iterations for each algorithm mentioned above is set to 100, with a population size of 40. This example uses the MATLAB 2021a simulation platform, with a computer model of Thinkpad X13, a processor model of Intel (R) Core (TM) i5-10210U CPU @ 1.60 GHz, and a memory capacity of 8 GB. The resulting running results are shown in the following figure (Zhou et al., 2024).

From the comparative analysis of Figure 3, it can be seen that compared to the other two algorithms, the improved Antlion algorithm has a longer computation time due to the addition of adaptive Levy flight and golden sine algorithm modules. The initial iteration curves of the three algorithms all approximate a linear descent, indicating that all three algorithms have fast optimization speeds in the initial stage. However, after multiple iterations, the convergence speed of particle swarm optimization algorithm and antlion algorithm is significantly slower than that of the improved antlion algorithm. The improved antlion algorithm obtained the optimal solution in the 18th iteration, the antlion algorithm obtained the optimal solution in the 44th iteration, and the particle swarm optimization algorithm obtained the optimal solution in the 56th iteration. Therefore, it can be seen that the introduction of the improvement factor in the improved antlion algorithm accelerates its convergence speed, enhances its optimization ability, and correspondingly increases its convergence accuracy. Therefore, it can be verified that the improved antlion algorithm used in this article has the advantages of strong search ability, fast convergence speed, and high accuracy in dealing with distributed power planning problems.

FIGURE 3

4.2 Example parameter settings

This section refers to the IEEE 33 node standard distribution network model, and the power flow calculation adopts the forward backward method. The system consists of 33 nodes, 32 branches, a reference voltage level of 12.66KV, a three-phase reference power value of 10MVA, and a balance node that is not connected to distributed power sources. The total load size of the system is 3,715 + j2300KVA, with nodes 1–8 and 17–26 selected for industrial load, and the total load size is 2,150 + j1045KVA; The residential load is taken from nodes 9–16 and 27–33, with a total load size of 1,565 + j1255KVA. Typical wind load scenarios are generated using the couple function and k-means clustering method, as shown in Figure 4.

FIGURE 4

The probabilities of each scenario are shown in Figures 5, 6, and the typical wind load scenario is taken as 10 for solving wind power planning.

FIGURE 5

FIGURE 6

4.3 Phase 1 distributed wind power planning and solution

Based on the above distribution network model and the output probability of the source load scenario, this section adopts the improved antlion algorithm mentioned in Section 2.2, selects 40 individual antlions, and performs optimization operations with a maximum iteration of 100 times. To verify the effectiveness of the proposed method, this section introduces three additional planning methods for comparative analysis. The results of the four solutions are shown in Table 1 and Figure 7.

TABLE 1

CategoryImproved antlion algorithmAntlion algorithmParticle swarm optimization algorithmNot considering uncertainty
Wind power access location and capacity12(722 KW)
31(427 KW)
8(567 KW)
24(406 KW)
14(450 KW)
26(496 KW)
4(410 KW)
19(488 KW)
Total capacity1149 KW973 KW946 KW898 KW

Power planning methods under different methods.

FIGURE 7

Compared with the other three schemes, the planning method solved by the improved antlion algorithm has the highest power access capacity and penetration rate, laying a good foundation for the economic and reliable operation of the system. To further verify the effectiveness of the method proposed in this article, this section will analyze it one by one from three aspects: voltage distribution, active power loss and voltage deviation, and system comprehensive cost.

4.3.1 System voltage distribution

Figure 8 shows the system voltage distribution diagrams corresponding to the four schemes. As shown in the figure, the improved antlion algorithm results in the most significant voltage increase effect, further verifying the superiority of the planning method based on the improved antlion algorithm.

FIGURE 8

4.3.2 System active power loss and voltage deviation

The active losses and voltage deviations of each scheme are shown in Table 2. From the table, it can be seen that the improved antlion algorithm has the smallest active power loss and voltage deviation compared to the other three schemes, which can effectively improve the economy and stability of power system operation compared to other schemes。

TABLE 2

CategoryImproved antlion algorithmAntlion algorithmParticle swarm
Optimization algorithm
Not considering uncertainty
Active loss62.17 KW86.09 KW89.24 KW94.31 KW
Voltage deviation0.1464p.u0.1572p.u0.1597p.u0.1684p.u

Active loss and voltage deviation.

4.3.3 System comprehensive cost

The system cost diagrams for different schemes are shown in Figure 9:

FIGURE 9

The improved antlion algorithm planning method has the highest power supply capacity and penetration rate, resulting in higher power investment and operating costs compared to other schemes. However, due to the reasonable planning of this scheme, the system line losses, environmental pollution, and power purchase costs are significantly reduced, effectively offsetting the high investment and operation costs. Table 3 shows the comparison of comprehensive costs under different schemes.

TABLE 3

CategoryImproved antlion algorithmAntlion algorithmParticle swarm optimization algorithmNot considering uncertainty
Operating cost (10,000 yuan/year)88.7369.9767.1063.56
Investment cost (10,000 yuan/year)145.86128.96126.72121.04
Network loss (RMB 10000/year)18.7129.2930.1627.43
Environmental pollution (10,000 yuan/year)56.4867.3162.1369.18
Government subsidy (10,000 yuan/year)68.1957.3853.0951.34
Electricity purchase cost (10,000 yuan/year)257.85273.71284.73291.65
Comprehensive cost (10,000 yuan/year)499.44511.69517.75521.52

Comprehensive costs under different schemes.

The cost of line loss for improving the antlion algorithm planning method is 187,100 yuan, which is the smallest compared to the other three schemes. This indicates that the proposed improved antlion algorithm can effectively assist distributed power generation in reducing line network loss and improving the system’s economic efficiency. In terms of purchasing electricity costs, due to the integration of distributed power sources, the self supply level of the system has increased, and the demand for purchasing electricity from the power grid has correspondingly decreased, thereby reducing the cost of purchasing electricity. Compared with the other three schemes, the introduction of the improved antlion algorithm further reduces the system’s electricity purchase cost, verifying the superiority of this algorithm.

4.4 Phase 2 distributed energy storage planning and solution

Based on the planning model described in Section 1, combined with the improved Antlion algorithm, the following energy storage charging and discharging and configuration planning strategies can be obtained. Table 4 and Figure 10 show the layout results of the energy storage device.

TABLE 4

CategoryImproved antlion algorithmAntlion algorithmParticle swarm optimization algorithmNot considering uncertainty
Energy storage access location and capacity33 (670 kW)4 (732 kW)10 (547 kW)24 (494 kW)

Distributed energy storage configuration schemes under different methods.

FIGURE 10

According to the particle swarm algorithm and planning methods that do not consider source load uncertainty, the energy storage installation capacity is relatively small. In the solution based on the antlion algorithm, although the energy storage access capacity is large, due to the proximity of the 4 nodes to the front end of the line and the higher economic cost of large energy storage capacity, the antlion algorithm solution is not suitable as the optimal processing method compared to the improved antlion algorithm. The introduction of improved antlion algorithm makes the energy storage installation capacity moderate, with good economic benefits, and the installation location close to distributed power sources, which is in line with the original intention of energy storage installation, that is, to suppress real-time power fluctuations caused by the integration of distributed power sources.

Figure 11 shows the charging and discharging status of energy storage and the charging and discharging power situation. It can be seen from the figure that energy storage is not always working, but indirectly charging and discharging based on the real-time operating status of the system. In order to verify the effectiveness of energy storage access, this section selects typical days as shown in Figure 12 for validity verification.

FIGURE 11

FIGURE 12

As shown in Figure 13, with the addition of distributed energy storage devices, the load curve fluctuation of the system on typical days shows a significant easing trend compared to before the addition of energy storage devices, and the peak valley difference of load output is significantly reduced, suppressing power fluctuations.

FIGURE 13

Based on the energy storage charging and discharging strategy in Figure 11, it can be seen that during peak load operation, the energy storage device serves as a power source to assist distributed wind power in providing electricity to the system. However, during low load periods, the energy storage device acts as a load and plays a charging role. The integration of distributed energy storage devices is coordinated with distributed wind power generation. When wind power output is low, energy storage serves as a power source to supply power to the load, achieving reliable operation of the system.

5 Summary

This article introduces a distributed energy storage charging and discharging strategy, and proposes a two-stage wind storage coordination planning method, which is solved by an improved antlion algorithm. The main conclusions are as follows: (1) The improved antlion algorithm with adaptive Levy flight and golden sine algorithm can effectively improve the voltage distribution level of the system, reduce network losses, and further reduce the overall cost of the system, bringing good stability and economy to the system operation. (2) The integration of distributed energy storage effectively suppresses power fluctuations caused by the uncertainty of distributed power generation output, and cooperates with distributed wind power to provide power supply. Although the cost of the method proposed in this paper is not optimal, it combines economy and reliability, providing decision-makers with more diverse and reasonable planning methods to choose from.

This study is subject to certain limitations. First, the cost-benefit trade-off of the proposed method in specific application scenarios needs to be more accurately quantified. Second, the algorithm’s solution efficiency and its application potential in large-scale systems require further in-depth validation. Future research will focus on the coordinated planning of wind power and energy storage systems, incorporating multi-time-scale and multi-type flexible resources. Additionally, the adaptability of the planning framework to uncertainties and extreme operating conditions will be enhanced.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

SS: Funding acquisition, Resources, Writing – original draft, Writing – review and editing, Project administration, Validation, Conceptualization, Methodology, Supervision. JZ: Data curation, Formal Analysis, Conceptualization, Writing – review and editing, Writing – original draft. QX: Formal Analysis, Investigation, Writing – original draft. LS: Visualization, Writing – review and editing, Validation, Supervision. CW: Methodology, Software, Writing – review and editing. SX: Data curation, Formal Analysis, Writing – review and editing. MD: Methodology, Writing – review and editing, Software.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Science and Technology Project of China Yunnan Power Grid Corporation (YPGC) (Project Number: YNKJXM20222105).

Conflict of interest

Authors SS and QX were employed by Electric Power Research institute of Yunnan Electric Power Grid Co. Ltd.

The remaining 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.

The authors declare that this study received funding from China Yunnan Power Grid Corporation. The funder had the following involvement in the study: High proportion clean energy county multi energy collaborative distribution system operation control and resilience enhancement technology.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

References

  • 1

    AbualigahL.ShehabM.AlshinwanM.MirjaliliS.ElazizM. A. (2020). Ant lion optimizer: a comprehensive survey of its variants and applications. Archives Comput. Methods Eng.28 (3), 1397–1416. 10.1007/s11831-020-09420-6

  • 2

    AssiriA. S.HussienA. G.AminM. (2020). Ant lion optimization: variants, hybrids, and applications. IEEE Access8, 77746–77764. 10.1109/access.2020.2990338

  • 3

    BanihashemiF.LesaniH.MakhzaniA. S. (2011). Locating and capacity determination of distributed generations using none-dominated sorting genetic algorithm international conference in swarm intelligence. Springer, 211–298.

  • 4

    ChengshanW.KaiC.YinghuaX. (2006). Location and capacity determination of distributed power sources in distribution network expansion planning. Power Syst. Autom.03, 38–43.

  • 5

    China Electricity Council (2023). “China's electricity industry annual development report,”. Beijing: China Building Materials Industry Press, 56–59.

  • 6

    ChuzhuangQ. (2017). Fuquan distributed wind power source site selection and capacity determination considering wind power and load timing. J. Power Syst. Automation29 (10), 85–90.

  • 7

    DingQ.ZengP.SunY. (2020). A planning method for the placement and sizing of distributed energy storage system considering the uncertainty of renewable energy sources. Energy Storage Sci. Technol.9 (1), 162.

  • 8

    DeyiY.ZhengyouH.TianleiZ. (2011). Distributed power generation site selection and capacity determination based on adaptive mutation particle swarm algorithm. Power Grid Technol.35 (06), 155–160. 10.13335/j.1000-3673.pst.2011.06.032

  • 9

    FeiM. (2020). Analysis of the influence of the development of distributed power on the existing power grid. Int. J. Soc. Sci. Univ.3 (2), 98–105.

  • 10

    GangulyS.SamajpatiD. (2015). Distributed generation allocation on radial distribution networks under uncertainties of load and generation using genetic algorithm. IEEE Trans. Sustain. Energy6 (3), 688–697. 10.1109/tste.2015.2406915

  • 11

    GhaemiM.ZabihinpourZ.AsgariY. (2009). Computer simulation study of the levy flight process. Physicia A Stat. Mech. Its Appl.388 (8), 1509–1514. 10.1016/j.physa.2008.12.071

  • 12

    HaifengS.MengjinH.ZhiruiL. (2016). Distributed power generation planning with energy storage devices based on temporal characteristics power automation equipment. 36 (06) 56–63.

  • 13

    IbrahimH.IlincaA.PerronJ. (2008). Energy storage systems—Characteristics and comparisons. Renew. Sustain. Energy Rev.12 (5), 1221–1250. 10.1016/j.rser.2007.01.023

  • 14

    JasimA. M.JasimB. H.BaiceanuF. C.NeaguB. C. (2023). Optimized sizing of energy management system for off-grid hybrid solar/wind/battery/biogasifier/diesel microgrid system. Mathematics11 (5), 1248. 10.3390/math11051248

  • 15

    JianfangY.ShengL.JunjieW.XiaoyiL. (2025). Ant lion optimization algorithm integrating levy flight and golden sine. Comput. Appl. Res., 2349–2353.

  • 16

    JinnanW.HongqiangJ.Fang. MethodsY. (2019). Research on environmental pollution cost assessment in China. China Environ. Sci.39 (5), 2229–2238.

  • 17

    JunqiangW.BoZ.JianhuaZ. (2016). Adjustable robust optimization for distributed wind power planning in distribution network. Power Syst. Technol.40 (7), 1989–1995.

  • 18

    JunyangM.TaoM.HangY. (2018). Distributed wind power source location and capacity determination in distribution network based on improved genetic algorithm. Electr. Autom.40 (06), 38–41.

  • 19

    KaiyuanX. U. (2023). Experience summary on general arrangement of distributed wind power project. Eng. Technol. Res.6 (5), 115–118.

  • 20

    KeL.NenglingT.ShenxiZ. (2017). Multi objective planning method for distributed power generation considering correlation. Power Syst. Autom.41 (09), 51–57.

  • 21

    Ministry of Ecology and Environment. (2022). Announcement of carbon dioxide emission factors for provincial power grids in China [EB/OL].

  • 22

    MohammadH. (2014). “Moradi, S.M. reza Tousi,Mohammad abedini,” in Multi-objective PFDE algorithm for solving the optimal siting and sizing problem of multiple DG sources. Elsevier Ltd, 56.

  • 23

    PaivaF.CostaJ.SilvaC. (2017). A serendipity-based approach to enhance particle swarm optimization using scout particles. IEEE Lat. Am. Trans.15 (6), 1101–1112. 10.1109/tla.2017.7932698

  • 24

    PingJ.TiejiangY.YongC. (2018). Study on capacity planning of distributed wind-storage integration system. Power Capacitor and React. Power Compens.39 (1), 135–141.

  • 25

    SivaramM.BatriK.Amin SalihM. (2019). Exploiting the local optima in genetic algorithm using tabu search. Indian J. Sci. Technol.1, 1–13. 10.17485/ijst/2019/v12i1/139577

  • 26

    TanyildiziE.DemirG. (2017). Golden sine algorithm: a novel math-inspired algorithm. Adv. Electr. and Comput. Eng.17 (2), 71–78. 10.4316/aece.2017.02010

  • 27

    TruongK. H.NallagowndenP.ElamvazuthiI.VoD. N. (2020). Aquasi-oppositional-chaotic symbiotic organisms search algorithm for optimal allocation of DG in radial distribution networks. Appl. Soft Comput.88, 106067. 10.1016/j.asoc.2020.106067

  • 28

    WeiguoH.JunyongL.ZhenboW. (2016). Distributed power generation site selection and capacity determination considering redundant power and timing high voltage electrical appliances. 52 (03), 93–99.

  • 29

    XuM.LiS.GuoJ. (2017). Optimization of multiple traveling salesman problem based on simulated annealing genetic algorithm. ed. MATEC Web Conf. EDP Sci.100, 02025. 10.1051/matecconf/201710002025

  • 30

    YurongW.RuolinY.SixuanX. (2020). Capacity planning of distributed wind power based on a variable-structure copula involving energy storage systems. Energies13 (7), 3450. 10.3390/en13143602

  • 31

    ZechunH.YonghuaS.ZhiweiX. (2017). Multi-objective optimization planning of large-scale distributed energy storage in distribution networks. Proc. CSEE37 (1), 82–91.

  • 32

    ZhenqiC. (2021a). Peng minfang, shen Mei'e distributed power generation site selection and capacity determination considering source load uncertainty. J. Power Syst. Automation33 (02), 59–65.

  • 33

    ZhenqiC. (2021b). Site selection and capacity determination of wind and solar energy storage based on two-stage coordinated optimization Hunan university.

  • 34

    ZhouM.MaF.JinW. (2024). Long-term cost planning of data-driven wind-storage hybrid systems. Renew. Energy223, 120073–120082. 10.1016/j.renene.2024.120073

Summary

Keywords

distributed power generation, energy storage, adaptive levy flight, golden sine theory, improved the antlion algorithm

Citation

Su S, Zhao J, Xie Q, Shang L, Wang C, Xiao S and Deng M (2025) Two stage coordination planning method of wind power and storage considering uncertainty of distributed source-load. Front. Energy Res. 13:1633719. doi: 10.3389/fenrg.2025.1633719

Received

23 May 2025

Accepted

26 August 2025

Published

22 September 2025

Volume

13 - 2025

Edited by

Haoran Zhao, Shandong University, China

Reviewed by

Rudrodip Majumdar, National Institute of Advanced Studies, India

Xinshou Tian, North China Electric Power University, China

Updates

Copyright

*Correspondence: Jie Zhao,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics