A Hybrid Energy Storage System Strategy for Smoothing Photovoltaic Power Fluctuation Based on Improved HHO-VMD

. To solve the problems of large ﬂ uctuation of photovoltaic output power a ﬀ ecting the safe operation of the power grid, a hybrid energy storage capacity con ﬁ guration strategy based on the improved Harris hawks optimization algorithm optimizing variational mode decomposition (IHHO-VMD) is proposed. In this strategy, the improved Harris hawk optimization algorithm is used to adaptively select k and α in VMD parameters and decompose the photovoltaic output power and distinguish between correlated and uncorrelated modes. Similarly, the moving average method (MA) is used to extract the continuous component signal in the uncorrelated mode, and it is reconstructed with the related mode as the grid-connected power that meets the national standard. The hybrid energy storage system (HESS) is used to stabilize the ﬂ uctuation component signal. The minimum annual con ﬁ guration cost of the energy storage system is established as the objective function. The simulation results show that the improved algorithm reduces the cost of the hybrid energy storage system by 6.15% compared with the original algorithm, suppresses the power ﬂ uctuation, and improves the economy and stability of the system.


Introduction
In order to get rid of excessive dependence on fossil energy and accelerate energy transformation, China has proposed the goal of carbon peak and carbon neutrality.Increasing the proportion of photovoltaic grid-connected power generation in the power grid is an important way to achieve this goal [1].Due to the randomness, intermittence, and periodicity of photovoltaic power generation, the power of the input grid fluctuates greatly, which challenges the stability of the grid.
The photovoltaic system with an energy storage device can effectively solve the problem of photovoltaic (PV) output power fluctuation and improve the stability of the photovoltaic grid [2].However, too large capacity configuration will increase the system cost, and too small capacity configuration will lead to frequent switching of charging and discharging states of the system and shorten the service life of the energy storage device [3].Therefore, many scholars have done a lot of research on suppressing photovoltaic fluctuations and HESS capacity configuration.Reference [4] proposed a VMD-based hierarchical configuration method for HESS power.Because the selection of VMD decomposition parameters k and α has not been studied, it is subjective and empirical.Reference [5] used regenerative boiler and battery as HESS and solved the problem of frequent electrode replacement of regenerative electric boiler in a novel way.However, due to high cost and low efficiency, it will produce economic waste.References [6][7][8] used adaptive wavelet transform to distribute the output power once and obtained the power and energy storage power which satisfied the gridconnected condition.However, the empirical mode decomposition (EMD) proposed in this paper is prone to modal aliasing and poor decomposition effect and cannot accurately separate signals with similar frequencies.In reference [9], the output power of the photovoltaic power generation system is decomposed by ensemble empirical mode decomposition (EEMD).However, after the decomposition of the EEMD algorithm, the noise component of each mode will be increased, so that the center frequency will be higher, which will lead to poor decomposition and reduction effect.References [10,11] observed whether the center frequencies of the VMD modes overlapped and then selected the number of modes based on the advantages.But there are no clear measurement criteria; they are subjective.Reference [12] proposed a grid voltage stability analysis framework considering the uncertainty of photovoltaic power generation and load demand and adopted the Monte Carlo simulation.Reference [13] used the hybrid bat search algorithm and artificial neural network method to optimize the current PI controller parameters of the hybrid energy storage system with the minimum objective function as the fitness function.But the bus voltage fluctuation optimization effect is not ideal.Reference [14] mentioned a parameter adaptive VMD method, using the grasshopper optimization algorithm parameters to optimize, by combining kurtosis and correlation coefficient as the optimization function.Reference [15] introduced a hybrid energy storage system composed of gravitational potential energy storage, a battery (BAT), and a supercapacitor (SC) to improve the economy and flexibility of a renewable energy power system.In [16], a capacity configuration method of HESS considering the time scale of decommissioned batteries was summarized.The results of capacity under different discharge times of hybrid energy storage, different time scales of decommissioned power batteries, and different control strategies were considered.In [17], a power limitation algorithm and a SOC correction strategy are proposed to correct the SOC of HESS.Limiting photovoltaic power when power fluctuations exceed the HESS power capacity can improve the net income of photovoltaic and HESS power stations and smooth the performance of PV output power fluctuations.However, due to the limited power capacity of HESS, there will be a decrease in photovoltaic power output and a decrease in the power generation of power plants.
Due to the deficiencies in the above research, this paper proposes a HESS capacity optimization configuration strategy for IHHO-VMD.Firstly, for the problem that the signal decomposition effect is poor due to the subjectivity of VMD parameter selection, IHHO is used to adaptively select k and α in VMD.Then, the Euclidean distance method [18] is used to distinguish the correlation mode and the noncorrelation mode of the photovoltaic output power, and the MA is used to extract the continuous component signal in the noncorrelation mode, which is reconstructed with the correlation mode as the grid-connected power that meets the national standard, and the fluctuation component signal in the noncorrelation mode is suppressed by HESS.Finally, the minimum annual configuration cost of the energy storage system is established as the objective function, and the IHHO-VMD is used to solve the model to obtain the energy storage capacity configuration that meets the requirements.
The main contributions of this work are as follows: (1) Compared with the original HHO algorithm, the IHHO algorithm can improve the accuracy and iteration speed of the solution process and reduce the economic cost of HESS.The results show that the power of BAT is reduced by 8.54%, the power of SC is reduced by 2.47%, and the economic cost of HESS is reduced by 6.15% (2) By comparison, the results of this method show that the phenomenon of light abandonment can be reduced, and it is not much different from the original photovoltaic output, which improves the stability of the grid connection

Photovoltaic Grid-Connected System Model
2.1.Optical Storage System Structure.The optical storage system is composed of PV, HESS, and a power grid, as shown in Figure 1.HESS consists of BAT and SC.P PV is the photovoltaic output power; P ref is the power input to the grid; P HESS is the charging and discharging power of HESS; P BAT is the charge and discharge power of BAT; P SC is the charge and discharge power of SC; C 0 , C 1 , C 2 , and C 3 are DC/DC converters.The photovoltaic array battery and the hybrid energy storage battery are connected to the DC bus through the DC/DC converter.When P HESS > 0, HESS is charged; when P HESS < 0, HESS discharges.

Photovoltaic
Volatility.The grid-connected power of a photovoltaic power station should meet the fluctuation requirements [19], and the maximum fluctuation limit of active power is shown in Table 1.
In order to quantify the effect of HESS smoothing active power, the PV output power fluctuation index is shown in the following: where V t represents the photovoltaic power fluctuation rate within the sampling time interval t, P t max represents the maximum output power within the sampling interval t, and P t min represents the minimum output power within the sampling time interval t.P capacity represents the rated installed capacity. 2 International Journal of Photoenergy mode is distributed to the power-type energy storage element, and the low frequency is distributed to the energy-type energy storage element.The photovoltaic grid-connected power fluctuation is small, but the grid-connected power is quite different from the original photovoltaic output power, and the phenomenon of light abandonment is more serious.In this paper, P PV has decomposed into IMF1 and IMF2-IMFK intrinsic mode function components by IHHO-VMD; then, the correlated mode and uncorrelated mode are distinguished by Euclidean distance, and MA is used to extract the continuous component signal in the uncorrelated mode for grid connection, and the remaining fluctuations are compensated by SC and BAT; the hybrid energy storage power allocation strategy is shown in (Figure 2).

Power Allocation Strategy Based on IHHO-VMD
3.1.Principle of VMD Algorithm.VMD was proposed by Konstantin Dragominetskiy in 2014 [20] to solve the problem that empirical mode decomposition is prone to modal aliasing by iteratively finding the optimal solution of the variational model for complex nonlinear signals.The structure of the variational problem of the VMD algorithm is set as follows.The original signal F can be decomposed into K intrinsic IMF instructions, u k ðtÞ with central frequency and limited bandwidth, and u k ðtÞ has a certain sparsity to ensure that the sum of the estimated bandwidths of each u k ðtÞ is the smallest.The constraint condition is that the sum of F, and all u k ðtÞ is equal [21].The variational constraint is established by the following equation: where u k is the decomposed modal component, ω k is the center frequency of each modal component, δðtÞ is the unit impulse function, k is the number of decomposition modes, * is a convolution symbol, and f ðtÞ is the original signal.To solve the above constrained variational problem, a quadratic penalty term α and a Lagrangian multiplier λ are used to transform the constrained problem into a nonconstrained solution problem.The augmented Lagrangian multiplier is shown in the following: The alternating multiplier algorithm ( 3) is used to iteratively update the modal components u k , the center frequency ω k , and the Lagrange multiplier λ, and finally, the optimal solution of the augmented Lagrange multiplier is obtained.The update methods of u k , ω k , and λ are shown in Equations ( 4)-( 6), respectively,  3 International Journal of Photoenergy where τ is the noise tolerance of λ.In the case that the decomposed signal does not contain strong noise, generally only the quadratic penalty factor α is used instead of the Lagrange multiplier λ, and τ should be set to 0 to disable λ.
The VMD algorithm will iterate until the convergence error ε satisfies the following equation: Due to the principle of the VMD algorithm, it is necessary to set parameters artificially when decomposing the signal: the number of modes K and the binomial penalty term α .In the traditional VMD algorithm, the values of K and α are not accurate, which will lead to underdecomposition or overdecomposition of the original signal, resulting in the unreasonable allocation of energy storage power in the system.Therefore, this paper uses intelligent algorithms to solve the parameter combination ðk, αÞ in the VMD algorithm, which can adaptively adjust the VMD parameters according to the change characteristics and complexity of the signal, so as to solve the problem of relying on subjective judgment.

The Harris Hawk Optimization Algorithm.
The HHO algorithm is a naturally inspired algorithm proposed by Heidari et al. [22] in 2019, which is derived from observing the pursuit and escape behavior between the Harris eagle and its prey (rabbit).The algorithm has been widely concerned since it was proposed.Mathematical formulas are used to simulate the strategy of capturing prey under different mechanisms.In HHO, the Harris hawk chooses the solution, and the prey approaches the optimal solution with iteration.The HHO algorithm includes four stages [23].
The Harris hawk is first in a waiting state, carefully checking and monitoring the search space (lb, ub) to find prey and then maintaining an appropriate balance between exploration and exploitation.Finally, according to the detection of the previous stage, the raid is performed to attack the expected prey, while the prey tries to escape from danger.According to the escape and chase, the HHO algorithm proposes four strategies to simulate the attack behavior, referring to the formula in [24].
3.2.1.Soft Siege.Under this search strategy, the prey has enough energy to escape, and the position is updated by the following equation: where Xðt + 1Þ denotes the position of the Harris hawk at the next moment, XðtÞ denotes the current position of the eagle, ΔXðtÞ represents the difference between the current position of the prey and the eagle, X rabbit ðtÞ represents the position vector of prey, O is the prey escape energy, J repre-sents the random jump intensity of the prey escape process and r is a random number between 0 and 1.

Hard
Enclosure.It indicates that the prey energy is not abundant, and there is no way to jump out of the enclosure, so it is surrounded by hardness.The expression is shown in the following: 3.2.3.Gradual Rapid Subduction Soft Encirclement.This indicates that the prey is energetic and escapes successfully from the enclosure.At this point, you need to combine Levy flights to update the location of the prey: where Levy is the flight expression, D and S are the number and random vector needed to solve Levy flight, Z represents the eagle based on the Levy flight mode, and Y represents the eagle's next action assessment indicator.
3.2.4.Gradual Fast Subduction Hard Encirclement.Prey energy is not abundant, and there is no way to jump out of the enclosure.Therefore, a progressive, fast dive hard encirclement is required, and the position is updated as follows: where X m ðtÞ denotes the value corresponding to each dimension of all individuals.

Multistrategy Improved the Harris Hawk Optimization Algorithm (IHHO).
Although the existing HHO improvement strategy improves the exploration and optimization performance of the algorithm to a certain extent, it generally only improves or changes the energy reduction method for a certain update strategy in the HHO algorithm and does not effectively avoid the blindness of global and local optimization.In view of the above problems, reference [25] proposed a multistrategy improved HHO that combines the Cauchy mutation, random shrinkage exponential function, and adaptive weight.
The Cauchy function decreases gently from the peak to both sides.After updating the position by Cauchy mutation, the Harris hawk is constrained by the local extreme point and can jump out of the local extreme.The mathematical model of the Cauchy mutation updates the optimal solution to obtain the current global optimal solution X best by the following formula:   In order to improve the local mining ability of the algorithm, the neighborhood of the prey location is updated to find a better solution; this paper introduces an adaptive weight method.By fusing adaptive weight factor ω, the algorithm has better local mining ability.The adaptive weight formula and the prey position update are shown in the following formulas: where T is the maximum number of iterations and t is the current number of iterations.
3.4.Selection of Adaptive Function.The parameters of VMD are self-selective, and the values of noise tolerance τ and convergence error ε have no obvious influence on the reconstruction accuracy of the decomposed signal, which is generally set as the default values τ = 0:1 and ε = 10 −6 .In order to use IHHO to optimize the VMD critical decomposition parameters, it is necessary to construct a fitness function.
The fitness function contains two elements: the number of modes k and the quadratic penalty factor α, which have a significant impact on the VMD decomposition results.The improved Harris eagle optimization algorithm can be used to find the optimal values of parameters k and α.
In this paper, envelope entropy (EE) is used as the IHHO fitness function, and the envelope entropy reflects the sparse characteristics of the signal.Strong coefficient characteristics, small envelope entropy, and the coefficient characteristic are weak, the envelope entropy is big, and envelope information entropy is shown in the following equation: where IMF i,j represents the envelope amplitude of the jth sampling point of the ith modal component after VMD decomposition, N is the length of the modal component signal p i , and j is the normalized modal component envelope: where IMF EE ðiÞ represents the MEE of the ith modal component.Figure 3 is the IHHO-VMD flow chart using MEE as the fitness function.

Verify the Effectiveness of IHHO-VMD
In order to verify the effectiveness of the proposed IHHOoptimized VMD, the original photovoltaic output power is decomposed by IHHO-VMD, VMD, and EMD, respectively.The photovoltaic power plant with installed capacity of 150 kW in Gansu is selected, and the photovoltaic power data on June 28 is taken as the research object.The photovoltaic output curve of this day is shown in Figure 4.
The analysis of Figure 4 shows that the maximum power of photovoltaic output is 145.3 kW, and the maximum fluctuation rate of photovoltaic within 1 min is 64.8%.According to China's grid-connected photovoltaic power station standard, a 1 min photovoltaic fluctuation rate limit of 20% and a 10 min photovoltaic fluctuation rate of 100% are required to meet the requirements.However, the 1 min PV fluctuation rate is much higher than the upper limit standard of 20%, and the power supply reliability is low, which affects the safe operation of the power grid.It is necessary to stabilize and decompose the P PV to meet the gridconnected requirements.

Intelligent Algorithm Comparative Analysis Experiment.
In order to verify the advantages of the IHHO algorithm, the arithmetic optimization algorithm (AOA), HHO, and IHHO algorithms are used to optimize the VMD parameters.The number of iterations is shown in Figure 5.
The data depicted in Figure 5 reveals that the IHHO algorithm exhibits a minimum fitness function value of 3.236 in the 5th iteration.In contrast, the HHO algorithm demonstrates a minimum fitness function value of 3.240 in the 13th iteration, and the AOA algorithm shows a minimum fitness function value of 3.242 in the 14th iteration.These results indicate that the IHHO algorithm outperforms the other algorithms in both speed and accuracy for solving the problem.In order to further verify the superiority of the parameter optimization VMD algorithm proposed in this paper, parameter optimization VMD, VMD, and EMD [26] are used to decompose P PV .The optimal parameter combination ½k, α = ½10, 6440 can be obtained by iterative calculation of the IHHO algorithm in Table 2.And the Hilbert transform is performed on the decomposed IMF. Figure 6 shows the time domain and frequency domain of the modal components obtained by the decomposition of P PV by EMD.Nine modal components are obtained.It can be seen that the different frequency components are mixed seriously, and there are more harmonics.IMF1 and IMF2 almost cover the whole frequency band, and the high-and low-frequency components are not easy to distinguish.
Figure 7 shows the decomposition of P PV using VMD under the default parameters K = 4, τ = 0:1, ε = 10 −6 , α = 2000.As shown in the figure, VMD decomposition is better than EMD decomposition, and the modal aliasing phenomenon is not obvious, which is better than being unable to selectK andαwhile modal aliasing still exist.Figure 8 shows that using IHHO to optimize VMD to adaptively search for the best decomposition parameters, the best parameter combination ½k, α = ½10, 6440 is obtained, the mode mixing phenomenon is less obvious, and the amplitude is stronger.Therefore, compared with EMD and VMD, the parameteroptimized VMD can better complete the reasonable allocation of HESS power by distinguishing high-and lowfrequency components.[27], the correlated modes and uncorrelated modes are selected.The smaller the d is, the more relevant the submode component is to the original signal.The correlation mode can be determined by evaluating the slope between the Euclidean distance between the two adjacent submodes and the original signal.If the slope between submode j and submode j + 1 is the largest,  7 International Journal of Photoenergy j is the decomposition point of the correlated and uncorrelated modes.The Euclidean distance formula between twodimensional vectors, ðx 11 , x 12, ⋯,x 1n Þ and ðx 21 , As can be seen from Figure 8, as the number of submodal decomposition increases, the frequency of submodal components also increases, but the IHHO-VMD decomposition effectively alleviates the submodal aliasing  International Journal of Photoenergy phenomenon and has better separation and robustness.According to Formula (20), the Euclidean distance between each submodal component and the original signal is obtained.It can be seen from Figure 9 that the slope between IMF1 and IMF2 is the largest, so IMF1 is used as the correlation mode, while the others are used as the noncorrelation mode.

Selection Analysis of Uncorrelated Modes.
In this paper, for the processing of uncorrelated modes, reference [28] moving average method is used to extract the continuous component signals in uncorrelated modes.On the basis of the simple average method, the moving average is calculated by increasing and decreasing the old and new data in sequence.The calculation formula of the moving average method is as follows: where P st represents the continuous component signal of the uncorrelated mode at time t, P umt represents the uncorrelated modal signal at time t, P qt represents the wave component signal of the uncorrelated mode at time t, N represents the length of moving average period, and M represents the number of points measured.

Capacity Optimization Model of Hybrid Energy Storage System
5.1.Objective Function.A function is established to minimize the annual comprehensive cost of the hybrid energy storage system, and Formulas ( 22) and ( 23) are calculated in the following: where C is the annual comprehensive cost of the hybrid  9 International Journal of Photoenergy output power and grid input power should be consistent with the hybrid energy storage system output: where P HESS is the charge and discharge power of hybrid energy storage, P W is the output power of a photovoltaic power generation system, and P O represents gridconnected power.

Energy Storage Charging and Discharging Power
Constraint.The real-time output power of the battery and supercapacitor should be within the rated power range: where P B and P C are the rated power of the battery and supercapacitor; P B ðtÞ and P C ðtÞ are the charging and discharging power of the battery and supercapacitor at time t; P B ðtÞ ≥ 0 and P C ðtÞ ≥ 0 mean the charging of the battery and supercapacitor, respectively; and P B ðtÞ ≤ 0 and P C ðtÞ ≤ 0 mean battery and supercapacitor discharge, respectively.

SOC Constraints of Energy Storage.
The SOC of the battery and supercapacitor should be within the prescribed range: where SOC B ðtÞ and SOC C ðtÞ are the states of charge of the battery and the supercapacitor at time t, SOC B−min and SO C B−max are the upper and lower limits of battery state of charge, SOC C−min and SOC C−max are the upper and lower limits of the supercapacitor state of charge, η is the charge and discharge efficiency of battery and supercapacitor, PðtÞ is the power at time t, Δt is the sampling interval time, and E is the energy storage capacity.

Example Analysis
In this paper, the Euclidean distance method is used to select the uncorrelated mode, and the uncorrelated mode is processed by Formula (21).By reconstructing the uncorrelated mode continuous component signal and the correlated mode signal, the photovoltaic grid-connected power is obtained.
Figure 10(a) shows the method used in reference [29] to directly use the relevant modes as photovoltaic gridconnected power.Therefore, the maximum fluctuation rate of 1 min obtained by IHHO-VMD is 0.92%, which meets the national photovoltaic grid-connected standard.Although the grid-connected curve is relatively smooth, it is quite different from the original photovoltaic output, and the phenomenon of light abandonment is more serious.The method in this paper can be seen in Figure 10(b).The maximum power fluctuation rate of 1 min is 18.2%, which meets the requirements of grid-connected power stipulated by the state, and the grid-connected curve is more consistent with the original photovoltaic output, and the photovoltaic output is more fully utilized.The spectrum of the uncorrelated modal fluctuation component signal is shown in Figure 10.
By observing the spectrum diagram in Figure 11, it can be seen that the low-frequency fluctuation is mainly concentrated in IMF2 ~IMF4, the high frequency is mainly concentrated in IMF5 ~IMF10, and the modal aliasing part in IMF4~IMF5 is the least.Therefore, IMF4 is used as the demarcation point.The low frequency is compensated by the battery, and the high frequency is compensated by the supercapacitor.
It can be seen in Figure 12 that BAT compensates for low-frequency power and SC compensates for highfrequency power, and the distribution results also conform to the response characteristics of BAT and SC.At the same time, the switching times of BAT between charge and discharge are significantly less than that of SC, which plays a significant role in prolonging the battery life.10 International Journal of Photoenergy The parameters calculated by the hybrid energy storage system [30,31] are shown in Table 3.According to the annual average configuration cost and constraints of the energy storage system, IHHO is used to solve the annual average cost of a hybrid energy storage capacity configuration system.In order to verify the effectiveness of the proposed method, this paper compares it with the Harris hawk algorithm to optimize VMD (HHO-VMD-MA) and the arithmetic optimization algorithm to optimize VMD (AOA-VMD-MA).The different method configuration results are shown in Table 4.
As shown in Figure 13 and Table 4, the IHHO-VMD-MA method adopted in this paper can reduce the capacity configuration of HESS and the configuration cost of the energy storage power station.Compared with the HHO-VMD-MA method, the power configuration of BAT decreased by 8.54%, the power configuration of SC decreased by 2.47%, and the annual configuration decreased by 6.15%.Compared with AOA-VMD-MA, the power configuration of BAT decreased by 18.34%, the power configuration of SC decreased by 6.90%, and the annual configuration decreased by 19.15%.The results show that the method adopted in this paper can effectively reduce the economic cost.

Conclusions
Aiming at the fluctuation of photovoltaic output power, this paper proposes IHHO to optimize VMD parameters to allocate hybrid energy storage power and draws the following conclusions: (1) By dividing the photovoltaic output power into correlated mode and uncorrelated mode and using MA to connect the continuous component and correlated mode in the uncorrelated mode to the grid, the power-type energy storage material SC bears the high-frequency part of the power, and the energy-type energy storage material BAT bears the low-frequency part of the power, which can reduce the number of battery charging and discharging times and effectively prolong the service life of the battery (2) The combination of VMD and MA smooths the power of unrelated modes in P P and decomposes on the basis of meeting the upper limit of gridconnected volatility stipulated by national standards.
Compared with the traditional variational mode decomposition, it is more effective to improve the stability of photovoltaic grid-connected power (3) Compared with the traditional HHO, the improved IHHO can effectively jump out of the local optimal solution and improve the accuracy and convergence speed of the algorithm.The annual configuration cost is reduced from 76126 CNY to 71447 CNY, a decrease of 6.15% The proposed MA can effectively suppress the fluctuation of photovoltaic grid-connected power.However, there may be a phenomenon of limiting photovoltaic output power and reducing the power generation of photovoltaic power stations.
It is recommended that further research be undertaken in the following aspects: (1) In the part of photovoltaic power decomposition, better methods can be explored and adopted to reduce the phenomenon of light abandonment.For example, the application of the Kalman filter in suppressing photovoltaic fluctuations (2) In this paper, BAT and SC are used as HESS, and there are some limitations in energy storage capacity.
In the future, electric hydrogen can be used as part of HESS to reduce environmental pollution and reduce economic costs

Figure 2 :
Figure 2: Hybrid energy storage power allocation strategy.

3. 3 . 1 .
The Cauchy Variation.Aiming at the problem that the HHO optimization algorithm is easy to fall into local optimal value, the global search ability of the HHO algorithm is improved.Combined with the Cauchy operator, the Cauchy distribution function mutation effect is used to optimize 4International Journal of Photoenergy the global object.The standard Cauchy distribution function formula is as follows:

Figure 4 :
Figure 4: The original photovoltaic output curve of the day.

Figure 5 :
Figure 5: Comparison of fitness function values.

4. 3 .
Selection of Correlated and Uncorrelated Modes 4.3.1.Selection Analysis of Related Modes.According to the Euclidean distance theorem in reference

Figure 13 :
Figure 13: Convergence curves of different methods.

Table 1 :
Active power fluctuation requirements of the photovoltaic power station.

Table 2 :
The optimum parameters of each algorithm.
energy storage system; C B , C C , C Bom , and C Com are equal annual investment costs and annual operation and maintenance costs of battery and supercapacitor; P B , P C , E B , and E C are rated power and rated capacity of the battery and supercapacitor; K BP , K CP , K BE , and K CE are power cost coefficient and capacity cost coefficient of battery and supercapacitor; K Bom and K Com are the operation and maintenance cost coefficients of battery and supercapacitor; γ is the discount rate; and T B and T C are the service life of the battery and supercapacitor.

Table 4 :
The results obtained by configuring different methods.

Table 3 :
Related parameters of the energy storage system.