Optimal Conﬁguration of Energy Storage Systems in High PV Penetrating Distribution Network

: In this paper, a method for rationally allocating energy storage capacity in a high-permeability distribution network is proposed. By constructing a bi-level programming model, the optimal capacity of energy storage connected to the distribution network is allocated by considering the operating cost, load fluctuation, and battery charging and discharging strategy. By constructing four scenarios with energy storage in the distribution network with a photovoltaic permeability of 29%, it was found that the bi-level decision-making model proposed in this paper saves 2346.66 yuan and 2055.05 yuan, respectively, in daily operation cost compared to the scenario without energy storage and the scenario with single-layer energy storage. After accessing IEEE-33 nodes for simulation verification, it was found that the bi-level decision-making model proposed in this paper has a good inhibition effect on voltage fluctuation and load fluctuation after energy storage configuration. In addition, this paper analyzes the energy storage that can be accessed by photovoltaic distribution networks with different permeability and finds that when photovoltaic permeability reaches 45% and corresponding energy storage is configured, the economic and energy storage benefits of the system are the best.


Introduction
In recent years, with global carbon dioxide emissions hitting record highs, China has proposed a "two-carbon" target to tackle environmental problems.Promoting the development of new energy and the transformation of energy structures has become an important part of global development.Due to abundant reserves and easy access, solar energy has been developing rapidly in recent years, and its proportion in the power grid has been increasing year by year [1].While improving energy utilization, this has brought a lot of trouble to the power distribution network.With the continuous increase in the penetration rate of photovoltaics integrated into the power distribution network, problems such as voltage collapse may occur, which has a serious impact on the safe and stable operation of the system [2].
Studies have shown that a large number of photovoltaics connected to the distribution network will also increase the number of system equipment, which will bring a burden to the system and easily generate harmonic interference.In addition, the retrograde power generated by the grid connection is prone to exceed the limit of the system node voltage, which not only reduces the power quality but also deteriorates the user experience.When high-penetration photovoltaics are connected to the grid, the uncertainty of output cannot be matched with the load of the distribution network in real-time, which will affect the power balance of the system.When the photovoltaic output fluctuates greatly due to the change in environment and climate, the stability of the system will be affected [3].In addition, the high-penetration photovoltaic grid connection requires a large number of power electronic equipment to join the distribution network, which leads to the existence of harmonics and affects the power quality [4], and the dispatching flexibility of the distribution network is greatly reduced, which makes it more difficult for the power system to develop a power generation plan.If such problems cannot be properly solved, it will not only seriously threaten the safe and stable operation of the system but also cause a waste of energy and limit the future development of photovoltaic power generation [5].
The authors in [6][7][8] analyze the influence of photovoltaic systems from the aspects of voltage fluctuation, voltage amplitude, and frequency.From the perspective of stability, Rasoul proposed a new framework to analyze the influence of different photovoltaic permeability on voltage stability [9].During the study, Zetty found that in a high permeability renewable energy distribution network, load fluctuation is the main factor leading to the voltage fluctuation of the system, and the realization of various fluctuations in the high-light voltage permeability distribution network is important content to achieve the increase of photovoltaic permeability in the distribution network [10].The introduction of energy storage devices improves the power quality while improving the photovoltaic stable output [11][12][13].Through reasonable regulation and control of a BESS, the absorption of new energy on the power generation side can be completed, the permeability of distributed power supply on the transmission side and distribution side can be improved, and the safe, stable, and economic operation of the system can be ensured [14,15].The authors in [16][17][18] studied the working principle and characteristic analysis of different types of energy storage devices and different types of BESSs and discussed the practicability of combining BESS energy storage and generation measurement.From the perspective of photovoltaic and load output prediction, Rahman and Zhao verified the feasibility of combining energy storage optimization configuration with the prediction by comparing scenarios with or without prediction [19,20].In order to meet the photovoltaic energy storage demand in the distribution network, Wang's multiple operation scenarios of energy storage were divided into grid scenarios to obtain the demand relationship of energy storage capacity under different operating conditions and to complete the calculation of energy storage capacity [21].
Access to energy storage equipment requires considerable capital investment in actual project construction and operation and maintenance.Therefore, the demand response for energy storage capacity is important content in optimizing energy storage configuration.In [22], Balouch proposed an optimization goal of matching demand and supply.Based on the analysis of line planning, low-cost scheduling, and demand response, the energy utilization efficiency and comprehensive operating cost of a smart grid were optimized.The authors in [23,24] introduced the improved optimization algorithm to improve the optimization ability so as to determine the optimal scheme of energy storage optimization configuration and realize a higher degree of response between demand and supply by analyzing various indicators of access to the power grid.In [23], Balouch optimized a response scheduling scheme by introducing the GWCSO algorithm.Higher robustness and computational efficiency of the algorithm make the optimization results more advantageous in power cost and peaking ratio.In [24], Mostafa improves the PSO algorithm, improves the accuracy and effectiveness of the algorithm, and optimizes the location and capacity allocation of energy storage in distributed networks.While the optimization objects are complex and diverse when connecting to the power system, the choice and update of the Pareto optimal solution will determine the quality of the final optimization result [25].
In the existing studies, it seems obvious for everyone to apply energy storage in highpermeability photovoltaic distribution networks [26][27][28][29][30][31][32].In the case of low photovoltaic permeability, access to energy storage can indeed improve photovoltaic output and power quality.However, few researchers have analyzed whether energy storage can still meet expectations in the scenario of high photovoltaic permeability, and how to rationally allocate energy storage in a distribution network with high photovoltaic permeability.In this paper, the application of energy storage in a high permeability photovoltaic scenario is analyzed, and the energy storage in a high-light volt distribution network is configured by establishing a two-layer planning model of the distribution network.The optimal size of energy storage was configured considering the fluctuation of power grid voltage and load, economic benefits and energy storage benefits, and the working condition of energy storage in the scenario of high-light voltage permeability, and the improvement of benefits in all aspects of the distribution network were studied.Finally, the feasibility of the proposed method was verified in the IEEE-33 node system.
The main contributions of this study are summarized below: • Proposed a method for optimal allocation of energy storage capacity of a distribution network based on a two-layer programming model and verified its feasibility.

•
Used the K-means method to complete an analysis of the uncertain photovoltaic output into the deterministic scenario.

•
The multi-objective particle swarm optimization algorithm was improved to solve the optimal configuration, and the advantages of the improved algorithm were compared.

•
By constructing different scenarios, it was verified that energy storage can still improve the power quality in the distribution network with high-light voltage and permeability.

•
Through analysis of the optimal configuration of energy storage in the distribution network with different photovoltaic permeabilities, the optimal economic photovoltaic permeability was concluded.
The rest of this paper is organized as follows: The treatment method for the PV uncertainty and the selection of the PV working curve is introduced in Section 2. In Sections 3 and 4, the bi-level decision-making programming model is constructed and solved to realize the addressing and capacity selection of the energy storage device.At the same time, in the fourth section, the specific content of particle swarm optimization is described.In Section 5, four scenarios are constructed to discuss the benefits generated by energy storage configuration and optimization benefits brought by algorithm improvement.Finally, in Section 6, we summarize the content of the thesis.

Analysis of Photovoltaic Output Characteristic
Due to the great influence of light and the environment, photovoltaic power generation is full of uncertainties.For further analysis, we collected the annual daily output data of photovoltaic power stations (annual output of an operational photovoltaic power station in Henan Province from June 2019 to July 2020), as shown in Figure 1.When considering energy storage benefits, excessive uncertainty in output will lead to uncertainty in energy storage benefits.In order to avoid this influence, this paper will process various output curves by clustering the division method, summarizing the photovoltaic output with high uncertainty into six typical output scenarios and analyzing them, and transforming the uncertainty into a deterministic analysis.

K-Means Cluster Analysis Method
The K-means clustering method is a classical clustering analysis method based on the iterative method which has the advantages of high efficiency and convenience in processing large-scale uncertain data [33].Through the K-means clustering method, a large number of output data can be refined and extracted, and fewer typical output scenarios can be obtained that can represent the output of photovoltaic power stations.
The K-means algorithm sets an initial cluster center in all scenes and iterates clustering for a large number of scenes based on the optimal distance.The iteration is not finished until the clustering presents a steady-state equilibrium.The iteration results are shown in Figure 2.After the whole process is complete, the center of each cluster scene is set as the partition scene, and the probability of each cluster scene is set as the required partition probability P r (s).Set the number of initial scenarios ξ s (s = 1, 2 . . .N) to N s .The number of target scenarios is M s , and the entire calculation procedure is as follows: 1 M s target scenes with random data are set as the cluster center, and the set of these scenes is Center = ξ Center s (s = 1, 2 . . .M s ).
2 Excluding the cluster center set, the other scene set is set as Member = {ξ Center s } (s = 1, 2 . . .M s ), and the distance from the other scene set to the cluster center scene set is calculated: The other scene sets excluding the cluster center set are divided into the nearest cluster center according to the distance calculated in 2 .We obtain the cluster set Cluster = {C i }, i = 1, 2 . . .M s , where C i is a set of similar scenarios.4 Set the same cluster C i including L s scenarios.Add the distances from each scenarios to the others: and scene ξ k in CT s = Min(CT s ) is selected as the clustering center of the next iteration.This is used to calculate the next iteration cluster center set. 5 At this point, stable cluster centers and clustering results can be obtained by repeating steps 2 -4 .The probability number of each type of scenario is the probability number of a single scenario in that type of scenario.
The process of the clustering algorithm to reduce the scene is shown in Figure 3.

Selection of Typical Output Scenarios
After data processing and division by the K-means clustering method, six output scenarios as shown in Figure 4 can be obtained.The occurrence probability and the number of curves of each output scenario are shown in Table 1.As can be seen from Figure 4 and Table 1, the photovoltaic output power in Scenario 1 is low, while the output in Scenario 2 has great fluctuation and uncertainty.The intermediate level of output in Scenarios 3 and 4 cannot represent photovoltaic output, and Scenarios 5 and 6 have a high probability and good output curve.In contrast, Scenario 5 with the maximum annual output is selected as the typical photovoltaic output curve, which can better reflect the output characteristics of photovoltaic power generation.In order to facilitate the analysis of the combined effect of photovoltaic and energy storage under different permeability, in this paper, we will only select Scenario 5, which is the most representative and has the highest probability of occurrence at the same time, as the analysis object to study the influence of energy storage access to the power grid during daily operation on voltage fluctuation, operation cost, and other benefits of the distribution network.

BESS Bi-Level Decision-Making Model Configuration
Due to the mutual influence between the optimal configuration of the energy storage system and the stable operation of the distribution network, this will bring difficulties to the dispatching of the energy storage devices and will cause the operation stability of the distribution network to decline.Therefore, it is necessary to consider a reasonable location and capacity while taking into account the operation economy of the distribution network.The bi-level decision-making model relies on its own two-level hierarchical structure to optimize the system objectives hierarchically.The upper and lower levels influence each other and seek the overall optimal solution according to the independent objective function and the corresponding constraints [34,35].

Upper-Level Model Objective Function
In the upper-level optimization, energy storage configuration location, rated power, and installed capacity are considered to reduce the total cost of the energy storage system and distribution network investment and maintenance.The installation location and capacity of the BESS are optimized.After the optimal configuration of energy storage is obtained, the information is transmitted to the lower level to adjust the charge and discharge power of energy storage.
where F min is the minimum daily total cost after the energy storage is connected; f ope is the total cost of distribution network operation investment.f sto is the input costs for energy storage construction.X = [x 1 , x 2 , x 3 ], and x 1 , x 2 , x 3 , respectively, represent the BESS input node, power, and capacity.

Lower Objective Function
In the lower-level optimization, due to the influence of the energy storage installation location and capacity selection on the energy storage life, the lower-level decision-making model fully considers the change of the energy storage charging and discharging power to realize the economical operation mode of the distribution network and achieve the smallest fluctuation range of node voltage and load.The lower optimization objective function is as follows: where f 2 is the amplitude of the voltage fluctuation of the distribution network node caused by the access to energy storage, f 3 is the amplitude of the load fluctuation of the distribution network, and y = [y 1 , y 1 , y 2 . . .y 24 ] is the average hourly charging and discharging power of the energy storage system throughout the day.
(1) The voltage fluctuation of distribution network nodes caused by energy storage access can be expressed as: where N indicates the number of system nodes, and u k (t) indicates the voltage value of node k at time t.
(2) The load fluctuation of the distribution network caused by access to energy storage can be expressed as: where P ave represents the average load in a period of time when energy storage is connected.

Constraints
(1) Constraints at the BESS access node where x 1 are nodes invested in energy storage, and N max is the maximum number of nodes expressed as energy storage input.(2) The constraints of power rating and capacity energy storage devices can be expressed as where P max ess and P min ess is the maximum and minimum value of the rated output of the energy storage, and E max ess and E min ess are the maximum and minimum values of the energy storage input capacity.
(3) Power balance constraints P grid + P gv = P load + P loss + P ess (7) where P grid is the power value received by the grid, P gv is the PV output power, P load is the output power, P loss is the network loss, and P ess is the BESS input power.
(  (10) where SOC 0 is the initial state of the energy storage system, including power and capacity.E ssn is the rated capacity of the energy storage battery.( 7) Supplementary constraints 1 Due to the limitation of the SOC range of the BESS, there will be a large number of infeasible solutions during the recovery of its all-day charging and discharging power.If its charge and discharge power is processed, this will greatly improve the convergence rate in the solution process and reduce the amount of calculation.
where P ess (t) is the energy storage charge and discharge power value that has been processed at time t.In this way, the infeasible solution is transformed into an effective feasible solution, and the charging and discharging power that is not within the SOC range of the energy storage is changed to 0.
2 Using the penalty function method to deal with the constraints that are not within the valid range: (12)   where M is the penalty coefficient.g i (x) indicates the negative inequality constraint.h i (x) indicates the positive inequality constraints.k i (x) indicates the equality con- straint at zero.r, s, t are the number of constraints.

Solution of Model
For the bi-level programming model, this paper selects the genetic algorithm (GA) for the optimization of the upper layer and improved multi-objective particle swarm optimization (IMOPSO) for the optimization of the lower level.The calculation process is as follows: 1 Input the demand parameters of the distribution network into the system. 2 Initialize the decision variables for the upper level (including BESS installation location, power rating, and capacity).Under the constraint, the population and other parameters of the GA algorithm are initialized.3 Initialize the decision variables of the lower level, including the BESS charge-discharge method.Under its constraints, the IMOPSO algorithm population and other related parameters are initialized to solve the initial fitness of each optimization objective.4 After the optimization of the lower layer is completed, the TOPSIS multi-attribute decision-making method is used to select the upper Pareto solution set obtained, and the best scheme is selected and fed back to the upper layer to solve the fitness of the upper layer target.
The upper-level GA algorithm population is updated, and the third and fourth steps are continuously executed until the upper-level optimization is completed.6 The optimal BESS configuration scheme of the upper layer, the corresponding optimal charge-discharge method of the lower layer, and the optimal Pareto solution set are obtained.
The calculation process is shown in Figure 5.In the following sections, we will give a specific description of the improvement content of the multi-objective particle swarm optimization.

Improved PSO Algorithm
The particle swarm optimization (PSO) algorithm is an iterative optimization algorithm.By converting the target into a certain number of particles, the position is updated in each iteration, and the optimal solution is searched through continuous iteration.The updating method is as follows: x where w represents inertia weight; c 1 , c 2 are acceleration coefficients of particle motion respectively.r 1 , r 2 are randomly selected in the range of 0 to 1; id represents the ddimensional component of the particle numbered i in the optimal position vector at time k; and g

(k)
d represents the d-dimensional component of the optimal position of all particles at time k.
In the process of iterative optimization, the traditional muti-objective particle swarm optimization algorithm is prone to fall into local optimal and appear 'Premature conver-gence'.Therefore, this paper cross-mutates the prescribed bit vectors of particles to prevent them from falling into the local optimum.
X i,d = X min + (X max − X min )r (15) where X min represents the minimum position variable of the particle; r represents any value between 0 and 1. Random cross variation is carried out when the probability distribution P < P m is satisfied.When the fitness of particles tends to be the same or locally optimal, the w will increase.Otherwise, the w will decrease as the particles tend to disperse [36].
where f i is the fitness value of the particle i, w max ,w min are the maximum and minimum of w, and f max , f min , f avg are the average and maximum, minimum, and average fitness of all particles at present, respectively.
In order to make the optimal solution distributed evenly in a certain range, the Pareto solution set should be optimized step by step.In this paper, a dynamic image Pareto solution set updating method is adopted, as is shown in Figure 6.In the initial distribution stage of particles, referring to all particles of the original comparison rule, an image V iD is established with the optimal particle as an optimal Pareto solution set.With the progress of iteration, the average value of particles is obtained between every two images, and the optimal solution is used to make the speed and position of particles updated at this time so that the particle population moves towards the target direction.In the whole process, new particles constantly exchange information with image particles and update data in image particles continuously until the end of iteration: where t is the number of iterations; D is the dimension of the decision variable; p iD is the best historical value of single particle i; and pbest p is the optimal particle of the current Pareto solution set in the image.Average pbest p is the superior particle of the optimal Pareto solution set in the two images.There are differences in the selection of the optimal solutions between the multiobjective particle swarm optimization and single-objective particle swarm optimization, and the results obtained are complementary dominated Pareto solutions, which cannot be obtained by direct comparison of the particle fitness function.Therefore, in this section, the first 20% Pareto solutions with lower crowding distances and higher priority orders are randomly selected to guide the iterative updating of the particle population.

Multi-Attribute Decision-Making Based on TOPSIS Method
After solving the lower multi-objective optimization problem, the optimal solution obtained by the IMOPSO algorithm is a set of Pareto solutions, and the selection of the optimal solution is essentially a multi-attribute decision problem.TOPSIS is a method for ordering by similarity to an ideal solution, it selects the optimal solution set and the worst solution set through the established initial decision matrix, then compares the distance between the two solution sets and the evaluation index with the optimal solution set and finally sorts them to evaluate the pros and cons of the scheme.The TOPSIS method has high strict requirements in selecting weights.In this paper, the information entropy method is used to determine the weight of each target value.The information entropy method determines w by the difference of the target value in the Pareto solution, improves the accuracy of the final decision, reduces the difference, and ensures the objectivity of the decision.By using the TOPSIS method, we can determine a set of optimal Pareto solutions to guide us to choose an energy storage configuration scheme.
The optimal solution of the Pareto solution set obtained is selected from X 1 ∼ X N and combined into N alternative schemes.The scheme X i is selected from N records.It is the composition of some optimal solutions in the Pareto solution set.g m (X i ) represents the value of the mth attribute of the scheme X i .Since each attribute is different, it should be unified and changed into the same type.The new attribute value is G m (X i ) , which can be expressed as: where d(x i ) is the relative distance of scheme x i ; d + (x i ) represents the distance between scheme x i and the optimal solution.d − (x i ) represents the distance from solution x i , the negative worst solution.λ m indicates the weight value of g m (X i ), which is randomly set between 0 and 1. g m+ and g m− indicate the optimal and worst values of all schemes g m .

Case Description
In this paper, the proposed scheme is tested on the IEEE-33 node distribution network [37].In addition, the structure of the system is shown in Figure 7.
In this paper, the rated voltage of the selected distribution network is 12.66 kV, and the total load is 3715 kW + j2300 kvar.The upper and lower limits of the node voltage are specified as not exceeding ±5% of the rated voltage.Node 1 is a balance node, which is connected to the upper-level distribution network for power transmission.Taking into account the actual work and construction of photovoltaic power generation connected to the distribution network, photovoltaic power generation is connected to Node 9, and its installed power generation capacity is 1.077 MW (29% penetration rate).The typical PV output curve selected by the method in Section 2 is shown in Figure 8.The typical daily load curve in this area is shown in Figure 9.In this paper, the battery is used as the energy storage system for research and introduces the time-of-use pricing strategy proposed in [38].The specific time-of-use price is shown in Table 2.The energy storage control parameters are shown in Table 3.The specific setting parameters of the energy storage configuration optimization simulation are shown in Table 4.In order to study the actual effect of energy storage configuration, we first analyzed the specific benefits of a photovoltaic distribution network connecting to energy storage configuration and demonstrated that energy storage still has good benefits in the high-light volt distribution network.Then, we compared the photovoltaic distribution network scenarios under different permeability and analyzed and compared the change of photovoltaic permeability with the corresponding change of optimal energy storage configuration scheme.The specific analysis content is introduced in the following section.

Energy Storage Optimization Scenario Division
Analyze the effectiveness of the method proposed in this paper, set different conditions, divide it into four scenarios, and compare them one by one to verify the feasibility of the method:   By comparing Scenario 1 and Scenario 2, it can be found that the voltage amplitude curve of the photovoltaic distribution network is smoother after the energy storage is connected, and the voltage fluctuation and load fluctuation are reduced to a large extent, which indicates that the BESS plays a good role in suppressing the node voltage fluctuation and load fluctuation when it is connected to the distribution network.Compared with Scenario 3, the load fluctuation range of Scenario 2 is smaller, and the load smoothing capacity is better.At the same time, the total cost of Scenario 2 is 803.27 yuan lower than that of Scenario 3 economically, which verifies the good characteristics of the model proposed in this paper.
The optimal Pareto solution set distribution of Scenario 2 and Scenario 4 is shown in Figure 12.Scenario 2 adopts the improved IMOPSO algorithm in this paper to solve the inner model, and Scenario 4 adopts the unimproved MOPSO algorithm to solve it.The Pareto solution set in Scene 2 is more evenly distributed than that in Scene 4 due to the introduction of particle cross mutation, adaptive inertia weight, and the Pareto solution set update method of the dynamic image.Moreover, Scene 2 adopts multi-attribute decisions based on the TOPSIS method, resulting in a more diverse solution set.Through the comparison of scenes, it is obvious that the optimization results of the IMOPSO algorithm are obviously better than the MOPSO algorithm, and the search accuracy is higher.In order to compare the performance of the two algorithms, the external solution set and the spacing S are used in this paper to measure the optimization performance of the two algorithms.The S index refers to whether the particles in the Pareto solution set are evenly distributed in space.The mean variance of the particle density distance is used in this paper to characterize the uniformity and global nature of the population particles, as shown in Equation (23).
where I represents the average of all particles I(x i ) in the Pareto solution set.
According to the different internal environments of the two algorithms, after 20 cycles, the node voltage and load fluctuations in the optimization target are taken as the research object, as shown in Figure 13 and Table 6.By combining and comparing the charts, it was concluded that the IMOPSO algorithm proposed in this paper reduces the number of iterations in the node voltage fluctuation and load fluctuation, and the convergence performance is obviously better than the MOPSO algorithm.In addition, the improved algorithm and Pareto solution set update strategy make the solution set distribution more uniform and the type of solution set more diverse, and the improved MOPSO has better robustness and convergence.

Energy Storage Benefit Analysis under Different Photovoltaic Permeability
In order to verify the effectiveness of the dual-layer multi-objective optimal configuration model of the energy storage system proposed in this paper in the high-light volt permeability distribution network, the upper limit of photovoltaic power generation permeability was set at 60%, and the verification started from 30% permeability.Using the optimal configuration strategy of the BESS, the curve as shown in Figure 14 was obtained.When the photovoltaic permeability increases from 30% to 60%, the capacity and power of the energy storage system have an obvious rising trend.When the photovoltaic permeability reaches 50%, the growth slows down and tends to remain unchanged.In other words, it is of little significance to increase the capacity of the energy storage system when the permeability reaches a certain level.
Figure 15 below shows the variation trend of energy storage investment and the total cost of distribution network operation under different photovoltaic permeability.It can be clearly seen that the total cost of the system decreases first and then increases when the photovoltaic permeability increases, and the total cost is the minimum when the permeability is 45%.As the cost of photovoltaic power generation decreases with the continuous increase of the permeability but is limited by the load level, the cost of the energy storage system increases with the increase of the capacity.The interaction between the two makes the total cost of the system decrease to the minimum when the photovoltaic permeability is 45%.When the permeability increases again, the system's total cost will keep rising, and the system operation economy will be seriously affected.
Based on the discussion of the above two legends, it is found that the total capacity of the BESS should be controlled in the optimal range according to the actual situation, and the photovoltaic permeability should also be controlled at a certain value so as to ensure the system operation economy while ensuring the safe and stable operation of the system.In order to improve the overall economy of the system, this paper selected 45% photovoltaic permeability to verify and analyze the two-layer programming model of the energy storage system proposed in this paper.
As shown in Figure 16 below, after optimizing the configuration of the energy storage system with 45% photovoltaic permeability, the load curve of the distribution network presents an obvious smoothing trend, and the peak-valley difference decreases.The sufficiency proves that the two-layer optimal configuration model of energy storage can still effectively improve the off-peak load, reduce the peak load of the distribution network, and increase the scheduling flexibility of the distribution network under the condition of high photovoltaic permeability.

Conclusions
In order to ensure the power quality and scheduling flexibility of the photovoltaic distribution network with increasing permeability, this paper proposes a joint optimization operation mode of optical storage.Firstly, the PV output model was analyzed, and the scenario planning method was applied.The K-means clustering algorithm was used to divide the output scenarios, and the typical output scenarios were selected for analysis.A BESS two-layer decision model was established, and the improved IMOPSO algorithm was used to solve the two-layer model.The IEEE-33 node example was adopted, and the simulation verification was carried out based on the current feed-in price, selected energy storage parameters, and other parameters.The simulation analysis results are as follows: 1 Access to energy storage can effectively smooth the load fluctuation and voltage fluctuation of system nodes in a photovoltaic distribution network.In a distribution network with high-light volt permeability, energy storage can effectively improve the off-peak load of the distribution network and reduce the peak load, thus increasing the scheduling flexibility of the distribution network. 2 The bi-level programming model proposed in this paper has a good optimization ability for the rational allocation of energy storage.The improved IMOPSO in this paper has good convergence performance and robustness and has good applicability in application optimization 4 When the optimal energy storage capacity under different photovoltaic permeability is configured, the total cost of the system is optimal when the photovoltaic permeability is 45%, and when the permeability increases again, the total cost of the system will keep rising and seriously affect the operation economy of the system.The analysis of this paper provides a theoretical basis for the optimal configuration of the energy storage system and an important reference for the safe, stable, and economic operation of a high permeability photovoltaic distribution network.

Future Work
In this paper, only batteries are considered in the selection of batteries in the energy storage system.However, with a wider application of energy storage, a single energy storage system may not be able to meet the actual demand in the future.In subsequent research, we will combine other types of energy storage for optimization analysis of hybrid energy storage.
At the same time, because of the variety of renewable energy, more and more distributed power is connected to the distribution network.This paper only analyzes access to photovoltaic power generation.In a follow-up study, we will conduct a further study on scenarios with access to various energy sources.
In addition to the voltage fluctuation and load fluctuation considered in this paper, the power system with energy storage access has more indicators to measure security.In a follow-up study, we will also analyze the improvement and influence of energy storage access on various indicators.

Figure 1 .
Figure 1.Annual daily output curve of photovoltaic power plants.

Figure 7 .
Figure 7.The topology diagram of IEEE-33 bus system.

Figure 9 .
Figure 9.Typical daily characteristic curve of load.

Scenario 1 :
No energy storage.Scenario 2: With access to energy storage, use the IMOPSO algorithm in this paper to solve the optimization objective of lower-level model in the bi-level decision-making model; introduce the charging and discharging strategy of the energy storage system to simulate and analyze it.Scenario 3: When solving its single-level model, ignore the charging and discharging management strategy of energy storage in the lower model, and only the energy storage system and distribution network are considered to have the lowest total cost.At the same time, in the time-of-use electricity price model, the energy storage system is charged and discharged at a constant power regardless of the high or low electricity price.Scenario 4: The optimal configuration result of energy storage in Scenario 2 is used as the constraint condition of this scenario, and the traditional multi-objective PSO algorithm is used to simulate and analyze the lower model in the optimal configuration model of the energy storage double-level.Node voltage curves and load curves in different scenarios are shown in Figures 10 and 11 below, and Table5shows the optimization results of different scenarios.

Figure 10 .
Figure 10.System node voltage curve in different scenarios.

Figure 13 .
Figure 13.Convergence curves of external solutions for different objectives.

Figure 14 .
Figure 14.Changes in energy storage capacity and power under different photovoltaic penetration rates.

Figure 15 .
Figure 15.Total cost curve under different PV penetration rates.

Figure 16 .
Figure 16.Distribution network load curve before and after energy storage configuration optimization under 45% photovoltaic penetration rate.

Table 2 .
Time-of-use electricity price table.

Table 3 .
Energy storage control parameter table.

Table 5 .
Optimization results in different scenarios.
Figure 11.Load curve in different scenarios.