A Spinning Reserve Allocation Method for Power Generation Dispatch Accommodating Large-scale Wind Power Integration

Spinning reserve allocation is a critical problem for active power dispatch with large-scale wind power penetration. A risk-based reserve allocation method that accounts multiple control sub-area coordination is given in this paper. And a multi-objective optimization model is constructed to schedule the spinning generation reserve for online active power dispatch. A fuzzy optimization method is used to transform the multi-objective optimization problem into a single-objective optimization one. The relationship between loss of load expectation and spinning reserve is also derived. And a particle swarm optimization (PSO) method is employed to provide a numerical solution to the problem. Numerical tests on IEEE RTS system are also given to validate the proposed method. Nomenclature Sets wind G Set of wind farms. G Set of conventional units. Parameters jt D Predicted load demand of area j during time period t. t D Predicted load demand of the whole system during time period t. 5358 n Number of control sub-areas. f jt p Predicted wind power output of wind farm j during time period t. max jt p Predicted maximum available output of wind farm j during time period t. j λ Penalty coefficient for wind power curtailment. Production cost coefficients of thermal unit i. N Number of thermal units. 0 t /T Initial and last time period of the optimization.

Upper /lower limits of the production cost or operational risk index for conventional units.
i k The preference factor (k 2 is the cost preference factor and k 3 is the risk preference factor).x Predicted error of wind power and load demand.
The forced outage rate of unit i during time period t.

Introduction
Wind power generation is one of the most important renewable energy forms and it has seen an ever-increasing development in the past decades all over the World, but especially in China, which has now the largest wind turbine installation capacity in the World.The installed wind turbine capacity in China had reached 62,660 MW by the end of 2012 and is still developing fast.However, despite its benefits of low operating cost and low pollutants emission, large-scale wind power generation introduces great technical challenges for power systems because of its intrinsic intermittent and fluctuant output characteristics.In particular, due to its limited predictability of wind power, the operational reliability cannot be guaranteed with the conventional deterministic spinning reserve method, and extra spinning generation reserves are needed to accommodate wind power integration [3,4].
Power generation reserves can be classified into spinning reserves and non-spinning reserves according to the North American Electric Reliability Corporation (NERC) definitions [5].In China, the power generation reserves are divided into load following reserves, contingency reserves and maintenance reserves, where the load following and contingency reserves constitute the spinning reserves, and the maintenance reserves constitute the non-spinning reserves, which are used for different time-scale power generation dispatch and control [6].In the past several years, much research has been done to evaluate the spinning reserve requirements of systems with high wind power penetration.A multi-time-scale reserve power dispatch method [6] and reserve allocation method [7] are proposed to accommodate the integration of wind power generation.Lee [8] presented an evolutionary iteration particle swarm optimization algorithm for spinning reserve optimization of a wind-thermal power system, in which the outage cost as well as the total operation cost of thermal units is considered to evaluate the spinning reserve capacity.In reference [9], a cost-benefit analysis, which takes the outage cost into account, is presented for spinning reserve optimization to achieve a better economical balance.Reference [10] presents a load shedding probability index which takes the wind power uncertainty into account, and the system reliability is used as the measure for the influence of wind power penetration.In [11], a risk reserve constraint-based dynamic economic dispatch model is given, which schedules wind power under the restraint of the given risk threshold.In [12], the reserve cost is considered in the objective which is comprised of the penalty cost of underestimation and the reserve cost of overestimation of the available wind power.Morales et al. [13] proposes a stochastic programming model for spinning reserve optimization in the power system with high wind power penetration, in which both load shedding cost and wind spillage cost are accounted.
Although much work has been done on spinning reserve optimization, most of these studies are focused on evaluating the spinning reserve capacity requirements, and only a little work pays attention to the spinning reserve allocation problem for a reliable operation of multi-area power system [14,15].In fact, with large-scale wind power penetration, the power flow fluctuates intensively and drastically, making the transmission interface easily heavy or over loaded, which brings more pressure on the transmission interfaces [16].A limited transmission interface capacity may mean more localized reserve requirements.In such cases, a not well-optimized spinning reserve allocation result for a multi-subarea system may make the inter-zonal reserve supplies constrained by the transmission interface limits, resulting to the consequence of local reserve inadequacy and transmission interface overload, which could seriously damage the operational security of the power system [17].Thus, the spinning reserve allocation problem becomes especially serious for the grids with large-scale wind power integration.Reference [14] presents a model which includes area demand constraints, area reserve constraints and tie line capacity constraints between the modeled areas, for solving the multi-area economic dispatch problem with tie line constraints.Ref. [15] presents a LP-based joint dispatch approach, which dispatches energy and reserve market concurrently, for solving the multi-product, multi-zone/area physical market dispatch problem.However, the wind power characteristics are not taken into account in the above two references and thus they are only suitable for conventional power grids with no wind power penetration.For the power grid with large-scale wind power integration, the uncertainty characteristic of wind power output has to be considered and a probability-based method is required for this.
In this paper, the spinning reserve allocation problem of a multi-subarea system is investigated for active power generation dispatch.A risk-based multi-objective spinning reserve optimization model, which considers the coordination of multiple control subareas, is proposed.The relationship between loss of load expectation (LOLE) and spinning reserve is also derived.The fuzzy optimization and particle swarm optimization methods are used for this multi-objective optimization model.This paper is organized as follows: the spinning reserve allocation optimization model for active power dispatch considering operational risk, penalty cost of wind power curtailment and power production cost is proposed in Section 2, and the fuzzy optimization is applied for this model.Then the relationship between loss of load expectation (LOLE) and spinning reserve is derived in Section 3. Section 4 describes the calculation procedure of the reserve capacity for a multiple sub-area system.Next, a PSO-based method is applied for the optimization of the proposed model in Section 5. Numerical results on the IEEE RTS system are given in Section 6.Finally, some conclusions are drawn in Section 7.

Proposed Model
A well scheduled multi-subarea spinning reserve allocation result should effectively handle the problem of local reserve inadequacy and transmission interface overload.Hence, a risk sharing principle is applied in the proposed spinning reserve allocation method, which means that the operational risk is shared evenly by all the subareas, to avoid the possibility of a large blackout or wind power curtailment.And a reserve balance index t I with the variance form, which is widely applied for risk management in finance and electricity market [18][19][20], is defined as follows to evaluate the degree of operational risk differences among multi-subareas: The variance index was first used to represent the risk by Markowitz, who in 1952 defined risk as the volatility of expected return.When the variance becomes smaller, the data sequence becomes more stable, which means a lower risk for the investment process.What is sought in the process is the benefit maximization under certain risk levels or the risk minimization under some expected benefit levels.As the subareas may have quite different load demand levels, the same amount of load losses may be disastrous for a small subarea, while it is only a very small portion for a large subarea.Hence, a normalized operation on LOLE is necessary for risk evaluation among subareas.When t I gets smaller, the running risk gets more balanced and the possibility of a large local blackout or wind power curtailment gets lower, and vice versa.
Meanwhile, the economic and environmental considerations are another two primary factors in determining the unit output during the active power dispatch process.Therefore, the spinning reserve allocation problem is essentially a multi-objective optimization problem which can be formulated as follows: (1) Objective The objective of the spinning reserve allocation model can be expressed as follows: where f jt p can be obtained from Numerical Weather Prediction (NWP) model or statistical model like the persistence method depending on the actual conditions of the wind farms.And j λ > 0 must be satisfied to achieve the objective of a minimum possibility of wind power curtailment, which can be deduced from the equal incremental principle.

① LOLE threshold constraint
From the operational reliability aspect, a certain amount of spinning reserve must be reserved in advance.This is achieved by Equation (3): where the expression of t E will be given in Section 3, and thresh E is usually specified by the operators [21].From the mathematical optimization aspect, constraint (3) becomes looser when thresh E becomes larger, and the objective of the model becomes smaller, and vice versa.In fact, this reflects the tradeoff between security and economics of the power grid.And when thresh E is small enough, there can be no solution for the model.Hence, a pre-calculation of minimum thresh E can be done by setting t E as the objective instead of Equation (2) if needed.

② Active power balance constraint
The power balance between active power generation and load demand must be kept all the time during the normal operation status of the power system, which can be expressed by Equation ( 4), in which the power losses are neglected: where D t can be obtained from linear extrapolation method for short term time periods and its prediction technology is relatively mature now.

③ Unit capacity constraint
The following inequalities should be satisfied: where, max jt p should be greater than f jt p with the prediction error considered.
where 15 T is adopted in this paper according to the operating rules of China's power grid [22].Even if the operating rules are different for other power grids [23], the parameter can be easily adjusted.
⑤ Spinning reserve definition 10 min( , ) 10 min( , ) It can be seen from Equation (2) to Equation ( 12) that the proposed model is an operational-reliability-constrained [Equation (3)] risk-sharing optimization model, which is also a multi-objective optimization model.For this type of optimization problem, the traditional multi-objective weighting method is often used, but the values of the weight coefficients, which have a significant impact on the results, are difficult to determine, hence, the fuzzy optimization method, which is well capable of tackling multi-objective optimization problems with disparate and contradicting objectives [2], is used for this model to transform the multi-objective optimization problem into a single-objective optimization problem which can be easily solved.

Fuzzy Optimization-Based Multi-Area Reserve Allocation Model
The fuzzy optimization method, which provides an excellent mathematical treatment for multi-objective optimization problems, has been widely used in the literature for large multi-objective optimization challenges [2,24,25].With the introduction of membership function, the fuzzy optimization method can transform the optimization model into a single-objective optimization problem with the original constraints and an objective to maximize the minimum of such membership function.A satisfactory number of multiple contradictory objectives can be simultaneously considered in the optimization problem.The "best" compromised solution among the Pareto-optimal solutions of the multi-objective optimization problem can be obtained with this method [26].
For the proposed reserve allocation model, the membership function for the penalty cost of wind power curtailment can be formulated as Equation (13) and is plotted in Figure 1.As both the production cost 2 f and the operating risk 3 f are required to be as small as possible, a monotonically decreasing quadratic polynomial membership function is adopted as follows to reflect the operator's attitude towards security and economics of the system, and its shape is shown in Figure 2: where, 1 α , 1 β , 1 γ and 2 α , 2 β , 2 γ can be obtained from equation ,min ( ) The results are as follows: ) where [ 1,1] i k ∈ − , and its value plays a major role on the shape of the membership function as shown in Figure 2. 0 i k > represents a pessimistic attitude towards the relevant objective and the optimization tends to decrease the corresponding objective value; while 0 i k < represents an optimistic attitude towards the corresponding objective.The values of coefficient i k for different objectives define the compromise degree between economics and security.For example, if the operating cost is what to be more pursued, then a larger positive value should be taken for 1 k .Otherwise, if the operational risk is emphasized more, then a larger positive value should be set for 2 k .Especially, 0 i k = stands for a neutral attitude to the relevant objective, and with this i k value, the membership functions ( 14) and ( 15) can be simplified as: It can be seen that the membership function approaches its maximum value 1 when the corresponding objective reaches its minimum value, and the membership function gets its minimum value 0 when the optimization objective reaches its maximum value.Therefore, the optimization result becomes better when the membership function gets larger, and vice versa.With the above definition of the membership function, the original multi-objective optimization problem can be transformed into a single objective optimization model as shown in Equation ( 19): 2 ( )

Relationship between LOLE and Spinning Reserve
One of the key problems for the proposed optimization model is how to derive the expression for LOLE.In this section, the relationship between spinning reserve and LOLE is investigated considering wind power integration.For the power system penetrated with wind power, the contingency considered for LOLE should include forced outages of generators and predicted errors of load demand and wind power.Generally, more than triple generator failures are negligible as their probability is very small.Hence, only three scenarios are considered for LOLE calculation: (1) No generator outage and only the predicted error of load demand and wind power are accounted; (2) one generator outage, with the predicted error of load demand and wind power accounted; (3) two generators outage, with the predicted error of load demand and wind power accounted.The corresponding expressions for LOLE in these scenarios are derived as follows: (1) For the scenario of no generator outage, only the predicted error of wind power and load demand are accounted for LOLE calculation: As the predicted error variable x only takes the value larger than R t in Equation ( 20), the integral of − gives the reserve inadequacy expectation.Meantime, the product , gives the probability of no generator outage.Hence, 1t E gives the LOLE with no generator outage.It can be seen that Equation ( 20) is only correlated with the variable R t after integral.When R t becomes larger, 1t E becomes lower.And this is also applicable for Equations ( 22) and (23), which means more reserve, less risk in the power system.
It is indicated from the statistical data of wind power output that the wind power prediction error gets larger when the forecast period becomes further, and it approximately follows a β distribution for a specific area [27].However, there is generally a wide geographical expansion for large-scale wind farms with numerous wind turbines, so their output usually follows a normal distribution according to the central-limit theorem.Hence, a normal distribution with mean 0 and standard deviation 1 ( ) t σ is assumed for wind power forecast error in this paper.And load forecast error can also be formulated as a normal distribution with standard deviation 2 ( ) t σ .As wind power and load demand forecast errors are usually considered uncorrelated, the sum of the standard deviation of the two forecast errors can be expressed as Meanwhile, the instantaneous forced outage rate for a unit during time period t can be expressed as [28]: where A(0) and U(0) are the status of the unit during the initial time period.And if it is on, A(0) = 1, U(0) = 0, and vice versa.Obviously the steady-state forced outage rate is , then y will also follow normal distribution, and Equation ( 20) can be simplified as: where Q function is defined as: The LOLE expressions for the other two scenarios can be derived in a similar way and are described as follows.
(2) For the scenario of one unit outage, LOLE can be expressed as follows: , , , Similar to Equation ( 20), the integral of gives the reserve inadequacy expectation with one unit outage.Meantime, the product gives the probability of one unit outage.
Hence, 2t E gives the LOLE with one unit outage.Meanwhile, it can be seen that Equation ( 22) is not only correlated with the variable R t , but also correlated with , i t P and , i t R .
(3) For the scenario of two units outage, LOLE can be expressed as follows: ) Similar to Equation ( 20), the integral of gives the reserve inadequacy expectation with two units outage.Meantime, the product , , , gives the probability of two units outage.Hence, 3t E gives the LOLE with two units outage.
Summarizing the above three scenarios together, the overall LOLE expression can be obtained, which is required in the former part to formulate constraint (3) and objective t I .Meanwhile, it can be seen from the expression that the optimization model is a non-linear and non-convex one, from which only a sub-optimal solution can be obtained with state-of-the-art technology.
The above three scenarios have covered most of the operating states of the system.If higher-order outages need to be considered, a fast sorting method [29] can be used to rank the probability of failure.And the first k failures with relatively high probabilities can then be selected and their corresponding impact on LOLE can be derived in a similar way.
Taken the above LOLE expression of Equations ( 21)-( 23) into 3 μ and Equation (3) of the spinning reserve allocation optimization model (19), an explicit joint optimization model for the active power schedule P i,t and reserve R i,t can be formulated.For this model, the particle swarm optimization method depicted in Section 5 is used, with which the optimum reserve allocation results can then be obtained.

Calculation of t R for Multi-Sub-Area System
When there are multiple sub-areas in the system, R t for each sub-area will take a different form as the inter-zonal supplies of reserve may be constrained by the transmission interface limits.To illustrate this, a two sub-area system shown in Figure 3 In Figure 3, MN t R should satisfy the following constraints: , , If we denote t R for sub-area M and sub-area N respectively as , What should be noted is that t R is different in Equations ( 21)-( 23) as MN t R is different under different circumstances.If two different outage events with units i and j occur in sub-area M and sub-area N, respectively, LOLE for these two sub-areas can be respectively expressed as follows: , , , It can be seen that the only difference between Equations ( 27), ( 28) and ( 23) is that , .The LOLE expressions for scenarios of more than two sub-areas can be derived in a similar way.

Solution
The proposed spinning reserve allocation model ( 19) is hard to solve by conventional mathematical programming algorithms as its constraints are nonlinearly coupled and not continuously differentiable.Hence, a particle swarm optimization method, which doesn't presuppose any conditions for the model and no approximations and limitations are required for the characteristic of the objective function and search space, is used in this paper.The particle swarm optimization, first proposed by Kennedy and Eberhart [30,31], is a population-based optimization method which can generate a high-quality solution within shorter computation time and exhibits a more stable convergence characteristic than other stochastic methods [8,[32][33][34].Due to the relative simplicity, faster convergence and less strict requirements for parameter tuning, PSO has been successfully applied to various fields of power system optimization such as optimal power flow [35], economic dispatch [36][37][38] and reactive power and voltage control [39].The main procedure of this algorithm is described as follows: (1) With the position of the particle swarm as the generator output, and the velocity of the particle swarm as the adjustment, the following matrix is stochastically generated as the initial value of the particle swarm: A step-by-step initialization is adopted in the procedure to account for the ramp rate constraint of the generator output between any two consecutive time intervals.For the individual j of particle i, the following formula is used to generate its position during time period t: (2) The fitness function for every particle is then calculated.For the equality constraints like = g(x) b ，the following penalty function is adopted: [ max( ( ), 0) ]) in which the loss of load expectation is calculated according to Equations ( 21)-( 23).
(3) The best fitness position for particle i is calculated as follows: ( ) ( ( 1)) ( ( )) ( 1) ( ( 1)) ( ( )) The global optimal position is calculated for the whole particle swarm: The correction terms are calculated during this step for the position and velocity of particle, i.e.: where, [31]; ω max and ω min are commonly set to 0.9 and 0.4 respectively; (6) When the solution variation in K consecutive iterations is less than the given threshold or the maximum iteration number has been reached, the iteration terminates.Otherwise, returning to step (3) for the next iteration.Taken the obtained global optimal position ( ) Gbest k back into the equality constraints such as the power balance constraint (4) and the inequality constraints like constraints (3) and ( 5)-( 12), the feasibility of the solution can be checked.And if it is not satisfied, go back to step (1) to repeat the iteration.

Numerical Tests
The IEEE RTS system is used in this part as the test system to validate the proposed model, in which the units with the same type and on the same bus are merged as shown in Figure A1.The unit parameters are listed in Table A1 and unit #14 is set to be a wind farm with 700 MW capacity.The wind prediction error follows the normal distribution with mean 0 and standard deviation (0.044 0.006 ) The load forecast error also follows the normal distribution with mean 0 and standard deviation 2% of the predicted value.The failure rate for generator i is 0.005-0.0045i/13,and the threshold value for load shedding due to reserve inadequacy is LSI (load shedding incidents) = 4.The ramp rate for thermal units is set to 2% of the unit capacity.The load and wind forecast values are respectively shown in Tables 1 and 2. Each time interval is 15 min and the optimization horizon is 4 h.

Relationship between Reserve Capacity and LOLE
Taken the units' output fixed to their initial values, the reserve demand under different LOLE and LOLP (loss of load probability) values are shown in Figure 4.It can be seen that the reserve demand increases with time, which is mainly due to the time-increasing characteristics of the wind forecast error and the forced outage rate of units.The gradient of the reserve demand curve with LOLP index is smaller than that with LOLE index with time increases.Meantime, when LOLE decreases, and the reserve demand increases correspondingly.
The relationship between the reserve and LOLP/LOLE during the 16th time period is shown in Figure 5.
It can be seen from Figure 5 that when LOLE decreases, the gradient of the curve is correspondingly reduced, indicating a varying important degree of the reserve capacity on system security.When the security level is low, the reserve has a relatively larger effect on the system, and a remarkable improvement of the system reliability would be seen with only a small reserve increase.This effect is weakened when the security level is relatively high.However, the relationship between the reserve and LOLP is approximately linear in this regime.Therefore, it can be concluded that LOLE is a better index for system reliability than LOLP in this respect.

Fuzzy-Optimization-Based Multi-Sub-Area Reserve Allocation
The test system is divided into two sub-areas with the dark dashed line, as shown in Figure A1 in this part, in which the unit capacity of sub-area 2 takes up about 20% of the total generation capacity, while the load demand took about 47% of the total load demand.Hence, it is a power receiving subarea with the transmission interface between the two sub-areas heavy-loaded.The LOLE threshold is set to 250 MW, and the transmission interface limit is set to 600 MW. .The production cost with different operational risk indexes 3,max   f is shown in Figure 6.It can be seen that when the optimistic attitude towards system reliability is taken, i.e., k 2 takes a positive value and k 3 takes negative value, the production cost is the lowest.When a pessimistic attitude towards system reliability is taken, i.e., k 2 takes negative value and k 3 takes a positive value, the production cost is the highest.Particularly, the production cost is in the middle when a neutral attitude towards system reliability is taken (i.e., k 2 = k 3 = 0).(2) Effect of LOLERD Keeping the LOLE threshold and the transmission interface limit constant, the production cost and LOLE under different LOLERD values are given in Figure 7 in which the linear membership function is adopted.It can be seen from Figure 7 that when LOLERD decreases, the production cost increases while LOLE for the power receiving sub-area 2 decreases and the power sending sub-area 1 remains relatively constant.They gradually become consistent, indicating the validity of the proposed method on balancing the sub-area operational risk.The relationships between the reserve demand of sub-area 2, the production cost, LOLE and LOLERD are shown in Figure 9.
As can be seen, when the transmission interface capacity increases, the production cost decreases and the reserve demand of power receiving sub-area 2 increases.Meantime, the LOLE of sub-area 2 is remarkably reduced with only a small change of the power flow.The results indicate a high sensitivity of transmission interface capacity to both the economy and security, which also means that some measurements should be taken to reinforce the interconnection.It can be seen that when the proposed reserve balance factor is accounted for, the active power output of unit #5, which has a larger capacity and has reached its upper limit, is curtailed.The outputs of units #1 and #3, which have relatively smaller capacities, are increased.Thus, the reserve capacity of sub-area 2 is increased and the operational risk is correspondingly decreased.Overall, more reserve is allocated for the power receiving sub-area to cope with the wind power fluctuation or the generator outage when the transmission interface is heavy-load.
The proposed model has a heavier computational burden than a conventional economic dispatch problem, because of the nonlinear and non-convex constraints in the model.A Matlab program is developed for the proposed method, and tested on an Intel(R) Core(TM) 2 Duo CPU computer with 2 GB RAM.The average CPU time on IEEE RTS system is about 356.1 s, but its performance heavily depends on the initial values of the particle swarms.For the well-generated initial values, its computation time can be reduced to 56 s.However, in the worst cases, the calculation time is 756.07 s.As shown in Figure 12, we believe that there are still lots of works to do for improving the performance of the solution in the future.Meanwhile, as the main focus of this paper is to evaluate the effects of wind power generation on the allocation of spinning reserve for a multi-subarea power system, where the results can further be summarized to develop expert knowledge for real-time spinning reserve determination.And some optimization techniques like the successive-linearization technique can be investigated if a fast speed is required in the future work.(6) Reserve allocation results on a three-subarea test system In this part, the IEEE RTS test system is further divided into a three sub-area system with the red dashed line, as shown in Figure A1.The load demand of sub-area 3 took about 13.3% of the total load demand, and the capacity of the transmission interface was set to be 100 MW.The reserve allocation results comparison of the reserve capacity between the proposed method and conventional probabilistic method during the 14th time period with 3,max f = 0.5 is shown in Figure 13.It can be seen from Figure 13 that when the proposed reserve balance is considered, the reserve capacity of unit #10 which has the largest installed capacity in sub-area 3 is reduced.The reserve capacity of unit #11, which is the largest thermal unit in sub-area 1, is increased.Meanwhile, the reserve capacity in sub-area 2 keeps constant.Accordingly, the reserve capacity of sub-area 1, which has large-scale wind power integration, is increased and the operational risk is correspondingly decreased.

Conclusions
Due to the large-scale penetration of wind power and its high fluctuation characteristics, the power flow of the grid varies intensively and drastically, which often makes the transmission interface margin insufficient and results in problems of local reserve inadequacy.Thus, the spinning reserve allocation problem for multi-area system becomes especially prominent for operational security of wind power systems.And it is investigated in this paper.Firstly, a risk-based multi-objective spinning reserve allocation optimization model is proposed.The relationship between the loss of load expectation and spinning reserve taken large-scale wind power integration into account is derived.Then, a fuzzy optimization method is applied to transform the multi-objective optimization problem into a single-objective optimization one.Finally, a particle swarm optimization method is applied to solve it.
From the numerical tests on IEEE RTS system, it can be seen that when the proposed reserve allocation method is considered, more reserve capacity will be reserved for the sub-areas with large-scale wind power integration or high load shedding losses when the transmission interface is heavy-load.With more reserve capacity, the scenarios of local reserve inadequacy or transmission interface overload due to the variation of wind power or forced outage of conventional units can be handled well.In this procedure, the unit output is reallocated compared with the results of the conventional method, and the reserve of those generators with relatively large capacities decreases while those with relatively small capacities increases.Overall, the sub-area operational risk is balanced by the reserve allocation method, and the effectiveness of the proposed method is validated.

Figure 1 .
Figure 1.Membership function for the penalty cost of wind power curtailment.

Figure 2 .
Figure 2. Membership functions with different k values.
is used as the example.

Figure 3 .
Figure 3. Two sub-area example for t R calculation.


Adding the above two penalty functions into the objective yields the final optimization objective: represents respectively the impact of the local best and the global best position on the individual particles in the swarm.

Figure 5 .
Figure 5. Relationship between reserve and loss of load probability and expectation.
respectively considered in this part, in which the upper limit of production cost is set to 2

Figure 6 .
Figure 6.Production cost under different preference factor values.

Figure 7 .
Figure 7. Effect of LOLERD on the production cost and LOLE.

Figure 9 .( 5 )
Figure 9. Results with different transmission interface limits.(a) Relationship between reserve proportion of sub-area 2 and LOLERD; (b) Relationship between the production cost and LOLERD; (c) Relationship between sub-area LOLE and LOLERD. .

Figure 10 .
Figure 10.Active power output of the units in sub-area 2.

Figure 11 .
Figure 11.Spinning reserve capacity of the units in sub-area 2.

Figure 12 .
Figure 12.Calculation time of the PSO method.

Figure 13 .
Figure 13.Spinning reserve capacities of all the units in the system.

Figure A1 .
Figure A1.IEEE RTS system with two/three sub-area division.
Upward/downward spinning reserve contribution of unit i during time period t.
n Number of control sub-areas.fjtp Predicted wind power output of wind farm j during time period t. max jt p Predicted maximum available output of wind farm j during time period t. j λ Penalty coefficient for wind power curtailment.N Number of thermal units.0t/T Initial and last time period of the optimization.threshEThreshold of LOLE.i p / i p Maximum/minimum generation output of unit i. i pd Δ / li k Generation distribution shift factor of unit i to transmission interface l. l TL / l TL Positive/negative power flow limit of transmission interface l with bus loads excluded.λ/μ The failure rate and the repair rate of a unit.H The swarm size.iγ/ j ζ Penalty coefficient for equality/inequality constraints of the PSO optimization.fpTotal predicted wind power output results.Variables t E LOLE of the whole system during time period t.E LOLE with no/one/two generator outage(s) accounted during time period t.E LOLE of sub-area M/N during time period t. t I Loss of load expectation ratio deviation (LOLERD).ltTLPower flow of transmission interface l during time period t, which can be calculated from the output of generators and loads multiplied with their corresponding sensitivity.jtEGLOLE of area j during time period t. jt r LOLE ratio (LOLER) of area j during time period t. t r The average of LOLER of all the control sub-areas.wjtp Wind power generation schedule of wind farm j during time period t. it p Generation schedule of conventional thermal unit i during time period t. 1 f Penalty term for wind power curtailment.2fProduction cost of conventional thermal units.3fOperating risk index of the system.
Probability density function of the forecast error.
t σ Standard deviation of the forecast error.Spinning reserve of sub-area M /sub-area N during time period t. t ij x Position of unit j of particle i during time period t. t ij v Velocity of unit j of particle i during time period t.P Upper/lower limits of unit i with ramp rate constraints accounted.( ) i x k Position of particle i during iteration k. ω The inertia weight. 1 μ / 2 μ / 3 μ Membership function of wind curtailment/economics/security

Table 1 .
Predicted load demand data.

Table 2 .
Predicted wind power output data.