1. Introduction
In the northwest and northeast regions of China, which are rich in wind energy resources, the wind power generation has experienced rapid development in recent years, which brought a need of optimal planning of the access location and access capacity of large-scale integrated wind power plants to improve wind power consumption and to ensure the power system safely operation [
1,
2].
Power plant planning aims to solve the issue about when and where to build what kind of power plants with an overall consideration of power grid conditions, such as grid network structure, load development, power system reliability, economic efficiency, operation efficiency, and so on [
3]. The common optimal planning of conventional power plants considers the smallest system economic cost as the optimization objective with the load demand supplied and the power system operating safely. The system’s total economic cost usually comprises of two parts: the power plants’ investment costs, which are related to the installed capacity, and the operating costs, which are related to the actual generated capacity.
Regarding the issue of optimal planning of wind power plants, according to wind power’s generation characteristics, the peaking ability of other types of power generation will be influenced by the integrated large-scale wind power generation. In solving the access location and access voltage level of integrated wind power plants, the total capacity of integrated wind power, the constraints of grid network structure, and the power system’s reliability and economy need to be taken into account [
4,
5,
6]. In [
7], based on a specific power operation mode, the capable access capacity of wind power generation was solved by the interior-point method based on the affine scaling theory. In [
8], considering the constraints of the transmission lines’ limited capacity, system’s spinning reserve requirements, and regular power units’ output limits, the wind farm’s capable access capacity was calculated using the chance-constrained method. The literature [
9] proposed an optimal method based on a market constraints optimization algorithm to solve multi-wind farms’ access locations and access capacities considering quarterly load demand and wind speed distribution characteristics. The literature [
10] also analyzed the transmission network’s reliability with integration of large-scale wind power, based on which of the relationships between wind farms’ installed capacity and the transmission system’s expansion planning, as well as the balance between system reliability and investment benefit, were investigated. In [
11], based on the evaluation of annual electricity demand, average wind speed, and illumination intensity, an optimization model using a genetic algorithm was proposed to calculate the total access capacity, the optimal capacity proportion, and the installation position of a wind-photovoltaic power hybrid system.
Up to now, among the studies focused on the optimal planning of wind power plants at home and abroad, most research are carried out based on distribution networks [
12,
13,
14,
15]. Though abundant achievements have been made, the existing methods show limitations facing the optimal planning of large-scale wind power plants which are centrally integrated into the transmission network. In China, according to the regional distribution characteristic of wind resource, wind power based generation is mostly exploited in a large scale more than hundred-MW-level in the northwest and northeast areas, however, which are a thousand kilometers far away from China’s east load center. In the meantime, China’s current new energy development policy encourages wind power energy to be exploited vigorously under the condition that the installed wind power can be effectively consumed. In this situation, there may be dozens or hundreds of wind farms integrated in a central area, and the power generated by wind farms need to be sent to distant load centers through the main transmission network. Hence, the optimal planning of access location and access capacity of wind power plants group is facing the challenges of large amount of total access capacity, widely available exploited areas, and variously-capable access locations.
In this paper, considering the mutual restriction between the planning of large-scale wind power plants and the planning of the power system’s network, an optimal planning model of access location and access capacity of large-scale wind power plants is put forward. The optimization objective takes the power flow equilibrium degree, investment costs, and active network loss as multi-optimization goals, and the proposed model is solved by the improved differential evolution (IDE) algorithm. With the solution results described by the Pareto pattern, the weights of the multi-objective attributes of each solution are calculated based on the entropy weight method, according to which all the alternative solutions are ranked and the final compromised solution is selected by the method of technique for order preference by similarity to ideal (TOPSIS). Based on a practical planning need of large-scale integrated wind power in an actual power grid in China in 2020, the best wind power access solution is successfully identified under the complex actual grid circumstances using the proposed model, which verifies the model’s effectiveness.
2. Optimal Planning Model of Large-Scale Wind Power Generation
The plan of access location and access capacity of newly planned wind power plants influences the corresponding expansion of power grid networks, the operation reliability, and the economic efficiency of the power system. To avoid the probable investment losses caused by wind farm construction and network construction, and to reach the optimal overall social benefits, the planning of integrated wind farms’ access locations and access capacities is better carried out with the consideration of the planning of network expansion. However, in this situation, the issue about optimal planning of wind power plants turns to be a multi-objective and nonlinear planning problem.
2.1. Objective Function
(1) Power flow equilibrium degree
With large-scale wind power integrated into the power system, if the wind power’s access has not been well planned, the power flow disequilibrium in different transmission lines may happen due to the volatility and randomness of wind power output characteristics. When one transmission line’s power flow reaches its capacity limit while the other lines are still lightly loaded, the power system’s security and reliability declines. However, the balance degree of system power flow on different transmission lines can be improved by better general planning schemes of wind farms and system network expansion. Hence, the power flow equilibrium degree can be an effective index indicating the effectiveness of the planning scheme of wind power and network expansion.
The power flow equilibrium degree can be described by the standard deviation of different line’s power flow ratios, as written as below:
where
l is the amount of the transmission lines in the system,
rij is the power flow ratio between the practical power flow and the rated capacity of line from bus
i to bus
j, |
Pij| is the actual power flow,
Pijmax is the rated capacity of line, and
is the average of
rij.(2) Economic efficiency
On this issue about the optimal planning of large-scale centrally integrated wind power plants, the economic efficiency is mainly considered in two parts: the investment costs and the operation costs. Investment costs includes the expenses of wind farms’ construction, corresponding connection lines’ construction, and main power grid network’s expansion. Operation costs is mainly caused by system operation loss which changes for different planning schemes since the power flow distribution will be different. In the composition of system operation losses, the active network loss is mainly considered in this issue.
Taking the minimum sum of the investment costs and the active network loss as the economic optimization objective, the economic efficiency is described as:
where
nij is the amount of planning lines between bus
i and bus
j, and
cij is the price of each line,
Ploss is the active network loss of power system,
Closs is the cost coefficient of
Ploss, determined by the electricity price and annual hours of maximum load loss.
2.2. Constraints
The constraints of the proposed model mainly include the wind power access capacity constraint, the power and electricity balance constraint, the power system operation constraint, and the power grid network constraint, etc.
(1) Constraint of wind power access capacity:
where
Wins is the limit of total capacity of integrated wind power.
(2) Power flow constraint:
where
Y0 is the initial node admittance matrix, Δ
Y is the change of node admittance matrix with the planning lines.
yij is the admittance of each planning line.
(3) Constraint of conventional power plants’ output:
(4) Node voltage constraint:
(5) Constraint of transmission line power flow:
(6) Constraint of network expansion:
where,
and
are the amounts of existed and planned transmission lines between bus
i and bus
j, respectively.
(7) Constraint of wind power distribution at booster station:
where
Wi,max is the maximum limit of bus
i’s capable access capacity.
(8) Security and stability constraint
The planning of large-scale wind power access locations and capacity should satisfy the requirement in the standard
Guidelines for Power System Security and Stability [
16] written by the State Grid Corporation of China. This standard stipulates the reliability requirement for power systems under different stability levels including: (і)
N-1 condition, when the power system loses any one operating electric equipment; (ii)
N-2 condition, when the system loses any two operating electric equipment; and (iii) serious fault conditions.
3. Multi-Objective Optimization Process
3.1. Pareto-Based Optimal Solution Set
Since the proposed optimization model of large-scale wind power’s access location and access capacity is a multi-objective and nonlinear programming problem, the improved differential evolution (IDE) algorithm [
17] is adopted here to solve the Pareto optimal solutions. The basic process of the DE algorithm consists three main steps: (і) Mutation: for two random vectors in the population, their difference is set as the variable source to generate a new vector; (ii) Crossover: based on mutation, a new generation is evolved from the current population each round; (iii) Choice: the population generation evolves until the evolutionary population contains the optimal solution.
Considering the calculation complexity and the DE algorithm’s characteristics, the original DE algorithm’s process is improved as below:
(1) Improved mutation process
For two random vectors
p and
q in the original population, if every attribute in vector
p is not poor and at least one attribute is better than those in vector
q, it is called that
p and
q obey Pareto dominance and
p is the dominated vector. In the
Gth mutation process, a new temporary population will be composed by the way that vector
ui,G+1 begins to compare with its parent vector
xi,G: (і) If the two vectors satisfy the Pareto dominance principle, the dominated vector is chosen as the new temporary population; (ii) If not, the two vectors are both chosen as the new temporary population. In this competing way, the amount of vectors in the temporary population is between
NP and 2
NP, and the former
NP vectors become the next population generation by non-dominated sorting and crowded distance sorting based on the improved non-dominated sorting genetic algorithm [
18].
(2) Adaptive adjustment of control parameters
Using the IDE algorithm, the convergence and robustness of the solving process are directly affected by the control parameters
F and
CR. In the standard DE algorithm,
F and
CR are constants parameters, and it always takes a long time to find the most appropriate value. Therefore, the adaptive adjustment is adopted to calculate the value of
F and
CR as the population evolves, which can eliminate the inconvenience of manual calculation and increase the diversity of the population. The adaptive adjustment strategies of
F and
CR are described as follows [
19]:
where
r1,
r2,
r3,
r4 are random numbers in [0,1] obeying the uniform distribution.
τ1 and
τ2 are the adaptive adjustment possibility of
F and
CR, respectively. In this work, the values of
Fmin,
Fmax,
τ1, and
τ2 are set as 0.1, 0.9, 0.1, and 0.1, respectively.
3.2. Acquisition of the Final Compromised Solution
With the Pareto optimal solution set solved by the IDE algorithm, the attribute weights of optimization objectives are then calculated, and the weighted Pareto decision-making matrix is composed using the information entropy method [
20]. In the weighted Pareto decision-making matrix, all the alternate solutions are sorted using the TOPSIS [
21], based on the distances from each alternative solution to the ideal solution (IS) and the negative solution (NS), which are both calculated. After the comparison, the solution closest to the IS and at the mean time farthest away from NS, is chosen as the compromised Solution. The decision-making steps based on the TOPSIS method are as follows:
(1) Establishment of the weighted-standardized decision matrix
Firstly, establish the standardized decision matrix
R, of which the square sum of a same attribute of all solutions adds up to 1:
Based on Equation (13), the weighted-standardized decision matrix V = {vij} is then established through the way that each column of matrix R multiply by its corresponding weight wj.
(2) Calculation of IS and NS
Second, provided that a lower value of the objective attribute indicates a better solution, then the ideal solution
A+ and the negative solution
A− can be expressed as:
(3) Calculation of distance
Then the distances of every alternative solution to IS and NS are calculated based on
N-dimensional Euclid Norm, expressed as:
(4) Sorting of alternative solutions and decision-making
Finally, the approximate degree of each solution to the ideal solution is described by:
According to Equation (16), it satisfies 0 ≤ Ci ≤ 1, and if vi is the IS, Ci = 1, if vi is the NS, Ci = 0. Since the value of Ci directly indicates the approximate degree of each alternative solution to the IS, all the solutions in Pareto optimal solution set can be sorted according to the value of Ci, and the one which has largest value of Ci is chosen as the final compromised solution.
3.3. Overall Optimization Process
Matched with the proposed optimal planning model of large-scale wind power, the access location and access capacity of wind farms are calculated and optimized in the following steps:
- (1)
Initialization of the calculation model. With the original system data, parameters, and algorithms input, the original generation of wind farm planning solutions are generated. Then the system power flow for each original solution is calculated by the Newton-Raphson method, based on which of the adaptive degrees of each solution is fetched.
- (2)
Acquisition of the first test generation of wind power plants planning solutions. Based on the results of the adaptive degree calculation of each original solution, the adaptive adjustment is applied calculating the value of F and CR as the population evolves. After mutation and crossover processes, the first test generation of the wind power plants planning solutions is generated, and the adaptive degree of each solution is calculated as well.
- (3)
Calculation of the Pareto-based optimal solution set. The IDE algorithm is applied to calculate the Pareto optimal solution set, as described in
Section 3.1.
- (4)
Acquisition of the final compromised solution. The TOPSIS method is applied to rank all the alternate solutions and reveal the final compromised solution out as described in
Section 3.2. The overall optimization process is shown in
Figure 1.
4. Simulation Analysis
4.1. Simulation Models
4.1.1. Power Grid Conditions
The simulation analysis is carried out based on a practical regional power grid in China which already contains 7300 MW first-phase integrated wind power generation. A second-phase access project of 4350 MW wind power generation in the year 2020 is ready for planning. The power system circumstances in 2020 are given in
Table 1 in which the data are given according to the future forecasting of the constructions of power plants, substations, transmission lines, and the development of load demands based on the current data, and the corresponding diagram is shown in
Figure 2.
As shown in
Figure 2, the basic information and constraints of this power gird are as follows:
- (1)
This regional power system consists of three 750 kV substations: T3, T5, and T9, of which each 750 kV substation is able to connect up to sixteen 330 kV lines, and this regional system consists of seven 330 kV substations: T1, T2, T4, T6, T7, T8, and T10, of which each 330 kV substation is able to connect up to eight 330 kV lines.
- (2)
Wind farms in the first-phase and second-phase are all integrated at 330 kV substation. The economic transport capacity of 330 kV transmission line is taken as 309 MW and the limited continuous transport capacity to be 782 MW. The power output of each wind farm from its booster substation to 330 kV transmission line is assumed as 600 MW.
- (3)
Through HVDC (high voltage direct current transmission) lines, this regional power system transmits power with other regional power grids with amount of −5000 MW, 8000 MW, and 2000 MW to District 1, District 2, and District 3, respectively.
- (4)
The locations of wind farms booster substations has been determined, and the capable exploited capacity of wind power around each booster substations is also certain. The network structure of wind farms’ internal collection and the corresponding network losses are not taken into consideration in this simulation.
4.1.2. Wind Farms Access Conditions and Boundary Conditions
As shown in
Figure 2, the detailed information of installed capacity and the transmission scheme of the first-phase 7300 MW wind power plants are given in
Table 2.
Affected by the wind resources distribution limit, there are only eight alternate wind farm’s booster substations, named by W2-1–W2-8, capable for integration of second-phase 4350 MW wind power generation. The 330 kV substations capable for wind power integration are the substations T1–T6. The transmission line’s length from the capable wind-power booster stations W2-1–W2-8 to the capable 330 kV substation T1–T6 are given in
Table 3.
4.2. Optimal Planning of the Second-Phase Wind Farms
In optimization process of the access location and access capacity planning of the second-phase wind power generation, some key variables are set as: (1) the unit cost of 750 kV line and 330 kV line are 2.44 million RMB and 1.5 million RMB per kilometer, respectively; (2) the cost coefficient of active network loss is 1.5 million RMB per MW; (3) the scale amount of the original population is 200; and (4) the maximum number of calculation iteration rounds is 1000.
In the early stage of the practical optimization process, the power system’s power flow equilibrium degree changes only in a small fluctuation degree, which is far less than the variation scale of the economic efficiency. It is suggested that the economy of transmission line’s cost plays a leading role among the optimization objectives in the early stage of optimization process. Hence, the optimal access locations of wind power booster substations and the corresponding transmission lines’ expansion plan are preliminarily selected.
In the later stage of the optimization process, since the transmission lines’ costs are the same based on the acquired wind farm access location solution, the capacity distribution scheme of wind farms are optimized according to the optimization objectives of power flow equilibrium degree and the active network loss.
Eventually the Pareto optimal solution set is generated with 30 solutions regarding wind farms’ access locations and access capacity combinations, of which the target cost of transmission lines is 504 million RMB, the target distribution of power flow equilibrium degree is between 0.3578 and 0.3725, and the target distribution of economic efficiency is between 711.0~863.3 million RMB. The distribution trend of all the 30 solutions in the Pareto optimal solution set are shown in
Figure 3. The optimization objective’s attribute values of each solution are given in
Table 4.
Based on the Pareto optimal solution set, as shown in
Figure 3 and
Table 4, the attribute weights of the power plow equilibrium degree and active network loss are then, respectively, calculated as 0.0357 and 0.9643 by the information entropy method. This means that the economic efficiency plays a more decisive role in this optimization planning process of the second-phase wind farms’ access location and access capacity.
Based on the TOPSIS method, all of the alternate solutions in the Pareto optimal solution set are sorted, and Solution 3 turns to be the optimal compromise solution which has the smallest distance to the ideal solution and the longest distance to the negative ideal solution. In Solution 3, (1) the power plow equilibrium degree is 0.3725; (2) the economic efficiency is 711 million RMB. The 750 kV network in this regional power system does not need strengthening, the total length of the planned transmission lines from second-phase wind farms to the 330 kV substations is up to 336 km with a total investment of 504 million RMB; (3) the active network loss is 138 MW; and (4) Solution 3 also satisfies the power flow check that no transmission line’s power flow exceeds the rated capacity.
The detailed access location and access capacity information of the second-phase wind farms in Solution 3 are given in
Figure 4 and
Table 5.