Optimal Allocation of Energy Storage System Considering Multi-Correlated Wind Farms

With the increasing penetration of wind power, not only the uncertainties but also the correlation among the wind farms should be considered in the power system analysis. In this paper, Clayton-Copula method is developed to model the multiple correlated wind distribution and a new point estimation method (PEM) is proposed to discretize the multi-correlated wind distribution. Furthermore, combining the proposed modeling and discretizing method with Hybrid Multi-Objective Particle Swarm Optimization (HMOPSO), a comprehensive algorithm is explored to minimize the power system cost and the emissions by searching the best placements and sizes of energy storage system (ESS) considering wind power uncertainties in multi-correlated wind farms. In addition, the variations of load are also taken into account. The IEEE 57-bus system is adopted to perform case studies using the proposed approach. The results clearly demonstrate the effectiveness of the proposed algorithm in determining the optimal storage allocations considering multi-correlated wind farms.


Introduction
Due to the ever-increasing concern over the global climate change and fossil energy decrease, interest in wind power has drawn more attention.However, a high penetration of wind power raises a problem of system instability, caused by the nature of wind uncertainty.The goal of economic emission dispatch (EED) with the energy storage system is to minimize the power system cost and emissions meanwhile ensuring the stability and power quality with facilitating penetration of distributed wind resources.
In recent years, many researchers have focused on solving the EED problem associated with wind power [1][2][3][4][5][6][7].Kumar and Suna in [1] developed a chemical reaction optimization algorithm based on the chemical molecular reaction to optimize total cost in wind-fossil-fuel-based power systems.With the help of this method, both environmental emissions and transmission losses are reduced.Zhan et al. [2] presented a multi-objective evolutionary algorithm to solve combined economic and emission dispatch problem incorporating the wind power, which decreased both the cost and the emission.In order to reduce the emissions, Denny and O'Malley [3] developed a forecasting method for wind power, which is utilized in the dispatch decisions.In [5], Tan et al. utilized two-point estimation method combined with particle swarm optimization to solve the stochastic economic load dispatch problem in microgrid which integrates the renewable energy generation.A traditional low-carbon emission dispatch analysis combined with the impact of wind power was proposed in [7].
However, the integration of a significant amount of wind power into power systems causes crucial challenges with determining optimal emission reduction, which result from its uncertainties.Especially, many wind sites are geographically close; therefore, interdependency of wind speeds at different sites should be considered in EED [8][9][10][11][12].As mentioned in [13], various approaches such as Monte Carlo simulation and analytical solutions have been employed to model the correlation of the multiple dependent wind sources.Monte Carlo simulation technology is introduced in [14] to calculate the indices of load expectation and the loss of energy expectation considering the wind speed correlation between the wind farms.In [15], Vallée et al. used the method of Monte Carlo simulation to consider the impact of uncertain system parameters on optimal power flow, but any Monte Carlo method requires intensive simulations to reach convergence.In order to reduce the computational burden, the analytical probabilistic methods were proposed as follows.The authors in [16] made use of the first-order second-moment method (FOSMM) to account for the uncertainties and correlations of the system load.The cumulant method (CM) in [17] was developed to study probability load flow.Due to the simplification used in FOSMM and CM, these two methods may result in relatively large errors.PEM can overcome the above demerits of analytical methods, which utilizes deterministic routines in solving the probabilistic problems while requiring low computational simulations [18].In [19], an efficient point estimate method was proposed to take uncertainty of nodal data and line parameters into account in the load-flow computations and to estimate the corresponding variations in the solution.
In this paper, the wind speed correlation between the two wind farms and their joint distribution is modeled by the Clayton-Copula method combined with a new PEM to discretize the joint distribution.The main advantages of the proposed method are twofold: (1) high computational efficiency when modeling and discretizing the joint distribution and (2) easy scale-up when dealing with multiple (more than two) wind farms.Furthermore, combining the proposed modeling and discretizing method with Hybrid Multi-Objective Particle Swarm Optimization (HMOPSO) [20][21][22], this paper develops an optimum energy storage allocation method considering inter-dependent wind farms to minimize wind-penetrated power system operation cost and emission.In addition, the load variation is considered by peak and off-peak load levels.

Problem Formulation
The goal of this study is to analyze the economic allocation of ESSs by considering the joint distribution of multi-correlated wind farms.The EED problem with ESSs incorporating multi-correlated wind farms is formulated by the reduction of CO x emissions along with system cost minimization.The objective function encompasses the joint probability, power system operation cost and greenhouse gas emissions.

Modeling Multiple Dependent Wind Probability
Clayton-Copula method is used in this paper to model the joint probability of multi-correlated wind sources.Clayton-Copula method is a subset of the Copula method family which transfers multivariate marginal distributions of random variables into one-dimensional marginal [23][24][25].Among various types of copulas, the Clayton-Copula is mostly used to study correlated risks because of their ability to capture dependence and the easiness of the construction and the implementation [26].Unlike Normal-Copula and t-Copula, Clayton-Copula can be employed to joint marginal distributions together, which follows the different types of the distribution [27].
Let F XY (x, y) denote the joint cumulative distribution function with marginal F X (x) and F Y (y).Then there exists Clayton-Copula such that for all x and y: where a is the correlation coefficient.Then the bivariate probability density function f XY (x, y) can be obtained by Equation (2): With the help of the Clayton-Copula function, the marginal distribution of random variables for wind power and their dependency structure are constructed.

Discretizing Joint Distribution
Instead of using Monte Carlo method, a novel PEM is proposed in this paper to discretize the continuous multi-correlated wind power distribution.The basic concept of this discretizing method can be separated into two parts: (1) Calculate the first few moments of wind joint distribution when one of the wind farms is at its limits either zero or full wind power; (2) Find the remaining correlated finite points representing both wind farms generating power not at their limits.The procedure of the discretization are described as follows.

The Boundary Points
(1) Apply PEM such as described in [28,29] to the two wind farms separately.Assume that K and L discrete probability points are obtained for wind farms 1 and 2, separately.
(2) Calculate the joint probabilities ( (5) Calculate the probabilities ( K ∑ i=1 P i,max ) when the wind power in wind farm 2 achieves maximum combining with K different discretized wind power in wind farm 1.K boundary points are attained.
Note that the total number of discretized wind power in the boundary is 2K + 2L − 4 because 4 points are duplicating points.

The Interior Points
The remaining correlated discrete points (interior points) are calculated by integrating Clayton-Copula method with a PEM, which is given in below: where M 1 and M 2 are the rated wind power of wind farm 1 and wind farm 2, respectively; µ XY , σ XY , λ ij represents the mean, the standard deviation, and the central moment of (X, Y); P ij is the probability of interior points; P inner denotes the total probability of interior points, which is equal to one minus the total probability of boundary points; z = (X,Y)−µ XY σ XY denotes the standardized value of (X, Y).Consequently, altogether (K − 2) × (L − 2) interior points are finally determined.It was noticed that these interiors points are calculated not based on the discrete points obtained in the boundary point calculation and the total sum of all probabilities is equal to one, which is as shown in Equation (7):

Objective Function
The goal of this paper is optimally allocating the ESSs and generator outputs in order to reduce the total expected system operation cost and emissions, while considering the uncertainties and correlations of multiple wind farms and the variations of loads.The multi-objective functions are given by (8): where N is the total number of points of the discretized wind power; Prob peak−load and Prob o f f peak−load denote the probability of peak-load condition and off peak-load condition; Prob i is the probability of i-th discrete point for the joint distribution of two wind farms; Cost i or k is the total operation cost at the i-th or k-th discrete point ($/h) with the peak load or off peak load; Emission i or k is the total emission at the i-th or k-th discrete point (kg/h) with the peak load or with the off peak load: where NG is the number of generators; C(P Gj ) is the fuel cost of generator j ($/h); E CO x (P Gj ) is the CO x emission of generator j (kg/h); C s is the cost of ESS ($/h); a j , b j , c j are the fuel cost coefficients of generator j; θ j , δ j , γ j represent the CO x emission coefficients of generator j; c ops is the operation and maintenance cost of ESS, which is taken to be $29/MWh [30] in this paper; P storage denotes the power capacity of installed ESS (MW).
The first objective function in ( 8) is to calculate the total expected operation cost by optimally allocating ESS and determining the outputs of all the different types of generators factoring in the wind distribution as well as considering peak and off-peak loads.Due to the environmental concern, the emission equation in the second objective function of ( 8) is proposed to decrease the gaseous pollutants.

Problem Constraints
There are two types of constraints considered in this research: equality and inequality constraints.

Equality Constraints
These constraints (11) are related to the nonlinear power flow equations:

Inequality Constraints
The inequality constraints contain the bus voltages, the reactive power of generation and the tap of the transformer.The inequality constraints are shown as follows: where V i is the RMS value of the bus i voltage; Q Gi represents the reactive power of generator i, and T i denotes the tap of transformer i.

Solution Method
In order to solve the EED problem considering multiple dependent wind farms, the copula method combined with the proposed PEM is implemented to calculate the discretized joint probability of two inter-dependent wind farms, one in Madison and the other in Milwaukee.The algorithm of HMOPSO is applied to minimize the operation cost and emission and the algorithm consists of MOPSO with NSGA-II [31] for solving the multi-objective problems, and the probability power flow for obtaining the system status considering the probabilistic wind distribution.Additionally, the MOPSO with NSGA-II is detailed in the part of Appendix B.

Wind Distribution
A 10-year daily wind speed data for the cities of Madison and Milwaukee in USA is used to fit the Weibull distribution, which is shown in Equation ( 13): where k is called the shape parameter, and λ is the scale parameter.
To obtain the wind power distribution, a linear approximation equation is utilized as follows: where W is the injected power.X is actual wind speed.M is the maximum power of wind turbine.
α and β are the linear coefficients.V ci , V co and V no , respectively, denote the cut-in wind speed, cut-out wind speed and normal wind speed.The rated power for the cities of Madison and Milwaukee are 45 MW and 330 MW, respectively derived from 30% penetration.The cut-in, cut-out, and rated speed are 3.5 m/s, 40 m/s, and 13.5 m/s, respectively.The parameters are summarized in Table 1.Since Madison and Milwaukee are relatively close in geographical and metrological terms, the wind distributions of the two sites cannot be assumed to be independent so it is necessary to consider the wind speed dependencies.In this paper, a bivariate distribution between these two wind farms is constructed by Clayton-Copula method described in Section 2. By applying the Weibull distribution (Equation ( 13)) into Clayton-Copula method (Equation ( 2)), the bivariate probability density function f XY (x, y) can be rewritten as shown in Equation ( 15): where P 11−33 is the sum of the probabilities of the interior points.

Discretizing Wind Power Joint Distribution
Following the procedure given in Section 2, the joint distribution can be discretized as follows.

•
Applying the 5-PEM [28,29] to the wind farms separately, five discrete probability points are firstly calculated separately; • Calculate the probabilities ( (2) The interior points

•
The interior discrete probability points of the joint distribution are decided by Equations ( 3)-( 6) and (13).Three points are selected for each wind farm distribution in the calculation and a total of 9 discrete probability points are obtained.
The total discrete probability points including both boundary and interior points are 25.The sum of the 25 discrete probabilities is equal to one.Table 2 shows the results of the 25 discrete probabilities and their corresponding wind power contribution from each.More specifically, the first 16 points are boundary points, and the others are interior points.The Clayton-Copula method can be used to build multivariate distribution of dimensions more than two.The proposed PEM can also discretize the joint distribution of higher dimension easily.Supposing that there are three dependent wind farms and their joint probability can be calculated by Equation ( 16): Then the joint probability density function can be obtained: The joint probability of correlated wind power can be redefined as: Similarly, the new continuous distribution can be discretized using the new PEM.

Hybrid MOPSO for Economic Emission Dispatch
By incorporating a fast and elitist multi-objective algorithm and probabilistic load flow calculation, a HMOPSO algorithm is developed to search for the best combination of the placements and sizes of energy storage devices in power systems.A mixed-integer MOPSO [32][33][34] is applied to optimally locate and rate the storage.As a part of the probabilistic load flow, Newton-Raphson power flow calculates the power equations for the constraints.
The procedure of HMOPSO to solve the EED problem associated with multi-correlated wind farms can be summarized as follows: (1) Randomly generate population P with N particles for initializing all generators' voltage, output power and position and size of ESS.The random selections of swarm of particles considering constraints and corresponding velocity for each particle are initialized.(2) Model joint probability of multi-wind farms by Clayton-Copula method.
(3) Discretize the joint wind power distribution into 25-point distribution by the new proposed estimation method, which has been discussed in Section 3, 25 scenarios in two different load conditions for each study case are created.(4) Through probabilistic power flow, evaluate the particles by fitness function and recall their best positions associated with the best fitness value.( 5) Check and preserve the p best and g best , if the algorithm has not yet converged, update the p best and g best .(6) Duplicate population P to population Q to form a combined population R and update the position and velocity of each particle.( 7) Sort the members in population R through NSGA-II with elitism algorithm for selecting N best solutions to renew population P. ( 8) Repeat Steps 4-7 until all the scenarios are considered.

System Configuration
The proposed HMOPSO algorithm has been applied to IEEE 57-bus system and compared with several other methods in order to test its quality and robustness.The system consists of 7 generations and 42 loads, where bus 1 is the slack bus, buses 3, 6, 9, 12 are defined as power-voltage (PV) nodes, and other buses are real and reactive power (PQ) nodes [35].According to the sensitivity analysis, wind generations are added to bus 2 and bus 8 with the rated size of 45 MW and 330 MW, respectively.
A typical daily load profile for a community in the city of Madison is scaled.The peak-load condition, which is 1250.80MW, and the off peak-load condition, which is 750.48MW (60% of maximum load) are selected as variable load conditions to conduct HMOPSO algorithm in case studies.

Economic and Emission Analysis
The impacts of the integration of wind power, two different load conditions, and ESSs in five cases are studied and compared to demonstrate the effectiveness of the proposed method.
Case 1: A regular probabilistic load flow analysis for the system considering the entire joint wind power distribution, but without ESS installation; Case 2: An optimal load flow analysis to determine the best ESS allocation under the worst-case scenario assuming zero wind power for both wind farms; The resultant costs and emissions of Cases 1-3 are listed in Table 3.It can be seen from Table 3 that the total operation cost and emissions are reduced with the increase of wind power in Case 1.The operation cost in peak load condition varies from $108,816/h to $56,937/h according to the changes of total wind power from 0 to 375 MW and the emissions decrease from 86,162 kg/h to 29,556 kg/h.Similarly, the operation cost in off peak-load condition has the same pattern, which ranges from $23,660/h to $35,280/h with the emission from 24,289 kg/h to 6027.7 kg/h.However, since no ESS was considered in this case, the system confronts with both low voltage and high power loss problems.Specifically, when the total wind power is 0, the voltage at the wind generator buses and load area experiences a low voltage problem associated with a high power loss, which is 128.85MW.Under this situation, even if the outputs of the diesel generators reach the maximum, the low voltages still appear at buses 20, 30, 31, 32, 33 and 34, which drop to 0.96 p.u. below.

Case 2
In this case, the optimal load flow analysis is implemented to determine the locations and the corresponding sizes of ESSs under the worst operating condition, considering the zero wind power situation.In this case, buses 14, 45, 47 and 52 are found to be the best places to install ESSs with sizes of 29.98 MW, 26.30MW, 12.66 MW and 15.14 MW, separately.The system is then operated with the given allocation of ESSs, considering twenty-five different wind power situations in peak-load and off peak-load conditions; even though the total operation cost and emissions are significantly reduced comparing with Case 1, no matter the total size of the installed ESSs or the cost and emissions are higher than that in Case 3.
In addition, the total expected operation cost in Case 2 is equal to $44,309.38/h and the expected emission is calculated to 23,454.15kg/h by Equation ( 8) and the fuel cost and emission coefficients are given in Table 4. Case 3 In this case, the optimal ESS allocation considers the entire joint wind distribution.The proposed HMOPSO is implemented to minimize the operation cost and emissions through searching the optimal ESS allocation and generator outputs.In HMOPSO, the algorithm starts with randomly generating a swarm of particles and each particle initializes a random ESS size within constraints at each bus.The size of ESSs is evaluated by the cost and emissions function and will be updated by HMOPSO.In the end, the size of ESSs at some buses becomes zero, which means that these buses do not need to install any ESS.The remaining ESSs converge to their optimal allocations.As a result, the best placements for ESSs in Case 3 are buses 10, 32, 46 and 52 with the size of 12.89 MW, 13.12 MW, 16.50 MW and 29.87 MW, respectively.Table 5 describes optimization results for expected values of operation cost, greenhouse gas emissions and power loss in four cases.Similar to Case 2, the total operation cost and emissions are reduced, and voltage profiles are improved.The results from Table 5 show that the proposed method in Case 3, which considered the entire wind distribution and employed the HMOPSO, delivered the best results in all categories.Notice the total expected operation cost is calculated to be $38,522.24/h in Case 3, which is much less than that in Case 2, indicating that the allocation of ESSs selected by the proposed HMOPSO algorithm is much better than only using the worst-case scenario.Supposing that the system operated in one year (8760 h), the system in Case 3 will save $506,953, 46.4/year compared to the conventional optimization result, which is determined by the worst case.Furthermore, even though the total size of ESS in Case 3 is less than that in Case 2, the CO x emissions decreased more, which is 17.59% lower than that in Case 2.

Case 4
In order to demonstrate the effectiveness of HMOPSO, the method of MOGA is compared to solve EED problem in the same operation situation as Case 3. As shown in Table 5, the total operation cost and emissions are $40,754.70/hand 20,939.28kg, respectively, which are both higher than that in Case 3. Furthermore, the total size of ESSs found by MOGA is larger than HMOPSO, which is 80.78 MW.However, the system has to confront with a higher voltage deviation than that in Case 3 and the voltage deviation is calculated by , which is shown in Figure 1.
In addition, the computing speed of the HMOPSO is faster.Specifically, the time to operate the algorithm in Case 4 once is 224.52 s which is 12.35 s slower than that in Case 3.
Energies 2017, 10, 625 11 of 16 The results from Table 5 show that the proposed method in Case 3, which considered the entire wind distribution and employed the HMOPSO, delivered the best results in all categories.Notice that the total expected operation cost is calculated to be $38,522.24/h in Case 3, which is much less than that in Case 2, indicating that the allocation of ESSs selected by the proposed HMOPSO algorithm is much better than only using the worst-case scenario.Supposing that the system operated in one year (8760 h), the system in Case 3 will save $506,953, 46.4/year compared to the conventional optimization result, which is determined by the worst case.Furthermore, even though the total size of ESS in Case 3 is less than that in Case 2, the COx emissions decreased more, which is 17.59% lower than that in Case 2.

Case 4
In order to demonstrate the effectiveness of HMOPSO, the method of MOGA is compared to solve EED problem in the same operation situation as Case 3. As shown in Table 5, the total operation cost and emissions are $40,754.70/hand 20,939.28kg, respectively, which are both higher than that in Case 3. Furthermore, the total size of ESSs found by MOGA is larger than HMOPSO, which is 80.78 MW.However, the system has to confront with a higher voltage deviation than that in Case 3 and the voltage deviation is calculated by , which is shown in Figure 1.In addition, the computing speed of the HMOPSO is faster.Specifically, the time to operate the algorithm in Case 4 once is 224.52 s which is 12.35 s slower than that in Case 3. Case 5 In this case, the two wind distributions are assumed to be independent.The probability comparison between with and without joint wind distribution is shown in Figure 2. By utilizing the HMOPSO in the same way, the best allocation for ESSs without considering joint wind distribution

Case 5
In this case, the two wind distributions are assumed to be independent.The probability comparison between with and without joint wind distribution is shown in Figure 2. By utilizing the HMOPSO in the same way, the best allocation for ESSs without considering joint wind distribution are found to be 10, 32, 46 and 52 with the size of 26.37 MW, 17.01 MW, 9.35 MW and 25.86 MW, respectively.The total expected operation cost is calculated to be $39,455.60/hby Equation ( 8), which is $933.36/hhigher than that considering joint wind distribution in Case 3 and the total expected emission without joint wind distribution is 21,976.54kg/h, which is 12.05% more than that with joint wind distribution.Furthermore, all values obtained in this case are larger than those in Case 3 for the sake of ignoring the interdependency of multi-correlated wind farms.8), which is $933.36/hhigher than that considering joint wind distribution in Case 3 and the total expected emission without joint wind distribution is 21,976.54kg/h, which is 12.05% more than that with joint wind distribution.Furthermore, all values obtained in this case are larger than those in Case 3 for the sake of ignoring the interdependency of multi-correlated wind farms.
(a) (b) Case 6 With the rapid development of the renewable technology, three correlated wind farms are considered in this case.These wind farms are installed in buses 2, 3 and 8 with the rated power of 45 MW, 125 MW and 330 MW, respectively.By applying the HMOPSO in the same way, the best allocation for ESSs in Case 6 are found to be 10, 32, 46, 49 and 52 with the size of 24.88 MW, 23.64 MW, 26.10 MW, 7.76 MW and 28.49MW, respectively.The total expected operation cost is optimized to be $62,008.81/hwith 42,967 kg/h CO2 emissions.In addition, according to the hardware condition of Intel Core i5-4460, 3.20 GHz CPU and a 4.00 GB RAM, it takes 1053.86 s to operate the optimization algorithm once.

Conclusions
As the penetration of wind power continues to increase in the power grids, it becomes important to consider the uncertainties of multi-correlated wind power generations when optimizing the placements and sizes of ESSs.In this paper, the Clayton-Copula method incorporating with a new PEM is presented to calculate the discretized joint probability for two dependent wind farms.A HMOPSO algorithm is applied to solve economic emission problems associated with optimal energy storage allocation, where the variation of loads is also taken into account.Furthermore, unlike many other optimization methods which only consider the worst-case scenario or the expected wind power scenario, the entire multi-correlated wind power distribution is taken into account in the optimization.The approach of how to scale up to more interdependent wind farms is also discussed.An IEEE 57-bus system is adopted to perform case studies.The results illustrate the importance of considering the wind speed correlations between the wind farms and also show that the proposed approach is able to find optimal placements and sizes of ESSs for reduction of total cost and emissions.In the future study, both the annualized investment costs of the renewable energy and ESSs and the load profile along the year (8760 h) as well as the single cycle of the ESSs will be considered.

Case 6
With the rapid development of the renewable technology, three correlated wind farms are considered in this case.These wind farms are installed in buses 2, 3 and 8 with the rated power of 45 MW, 125 MW and 330 MW, respectively.By applying the HMOPSO in the same way, the best allocation for ESSs in Case 6 are found to be 10, 32, 46, 49 and 52 with the size of 24.88 MW, 23.64 MW, 26.10 MW, 7.76 MW and 28.49MW, respectively.The total expected operation cost is optimized to be $62,008.81/hwith 42,967 kg/h CO 2 emissions.In addition, according to the hardware condition of Intel Core i5-4460, 3.20 GHz CPU and a 4.00 GB RAM, it takes 1053.86 s to operate the optimization algorithm once.

Conclusions
As the penetration of wind power continues to increase in the power grids, it becomes important to consider the uncertainties of multi-correlated wind power generations when optimizing the placements and sizes of ESSs.In this paper, the Clayton-Copula method incorporating with a new PEM is presented to calculate the discretized joint probability for two dependent wind farms.A HMOPSO algorithm is applied to solve economic emission problems associated with optimal energy storage allocation, where the variation of loads is also taken into account.Furthermore, unlike many other optimization methods which only consider the worst-case scenario or the expected wind power scenario, the entire multi-correlated wind power distribution is taken into account in the optimization.The approach of how to scale up to more interdependent wind farms is also discussed.An IEEE 57-bus system is adopted to perform case studies.The results illustrate the importance of considering the wind speed correlations between the wind farms and also show that the proposed approach is able to find optimal placements and sizes of ESSs for reduction of total cost and emissions.In the future study, both the annualized investment costs of the renewable energy and ESSs and the load profile along the year (8760 h) as well as the single cycle of the ESSs will be considered.the size of P t+1 is larger than N.The new population P t+1 is used for the next iteration and, in every iteration, the individuals compete with each other by binary tournament selection operator.

( 3 )( 4 )
j ) when the output power of wind farm 1 is 0 MW combining with L different discretized wind powers in wind farm 2. L points in boundary are obtained.Calculate the probabilities ( 0 ) when the output power in wind farm 2 is 0 MW combining with K different discretized wind powers in wind farm 1.K points in boundary can be obtained.Calculate the probabilities ( L ∑ j=1 P max,j ) when the wind power in wind farm 1 achieves maximum combining with L different discretized wind power in wind farm 2. L boundary points are attained.
1 ,j ) when the output power of wind farm 1 is 0 MW or reaches the maximum combining with five different discretized wind powers of wind farm 2. Consequently, 10 discrete points are obtained;•Apply the same algorithm into wind farm 2 and the probabilities that are 5 determined.10 discrete points are obtained herein.

Case 3 :
HMOPSO with ESS considering the entire joint wind distribution; Case 4: Multi-objective genetic algorithm (MOGA) with ESS considering the entire joint wind distribution; Case 5: HMOPSO with ESS without considering the joint wind distribution; Case 6: HMOPSO with ESSs considering three dependent wind farms.

Figure 1 .
Figure 1.The expected voltage deviation for different cases.

Figure 1 .
Figure 1.The expected voltage deviation for different cases.

Figure 2 .
Figure 2. The probability comparison between independent and dependent wind distributions.(a) Independent probability; and (b) dependent probability.

Figure 2 .
Figure 2. The probability comparison between independent and dependent wind distributions.(a) Independent probability; and (b) dependent probability.

Table 1 .
Parameters of two farms.

Table 2 .
Probability of multi-correlated wind farms.

Table 3 .
Operation cost and emissions in different load conditions.

Table 4 .
Cost and emission coefficients of generators.

Table 5 .
Expected value in four cases.

Table 5 .
Expected value in four cases.
Energies 2017, 10, 625 12 of 16 are found to be 10, 32, 46 and 52 with the size of 26.37 MW, 17.01 MW, 9.35 MW and 25.86 MW, respectively.The total expected operation cost is calculated to be $39,455.60/hby Equation (