Mathematical Problems in Engineering

Volume 2017 (2017), Article ID 9153297, 11 pages

https://doi.org/10.1155/2017/9153297

## Two-Stage Robust Security-Constrained Unit Commitment with Optimizable Interval of Uncertain Wind Power Output

^{1}College of Energy and Electrical Engineering, Hohai University, Nanjing, Jiangsu 211100, China^{2}State Grid Jiangsu Electric Power Company, Nanjing, Jiangsu 210024, China^{3}Electric Power Research Institute of State Grid Jiangsu Electric Power Company, Nanjing, Jiangsu 211103, China^{4}School of Electrical Engineering, Southeast University, Nanjing, Jiangsu 210096, China

Correspondence should be addressed to Liudong Zhang

Received 16 February 2017; Accepted 11 June 2017; Published 20 July 2017

Academic Editor: Fazal M. Mahomed

Copyright © 2017 Dayan Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Because wind power spillage is barely considered, the existing robust unit commitment cannot accurately analyze the impacts of wind power accommodation on on/off schedules and spinning reserve requirements of conventional generators and cannot consider the network security limits. In this regard, a novel double-level robust security-constrained unit commitment formulation with optimizable interval of uncertain wind power output is firstly proposed in this paper to obtain allowable interval solutions for wind power generation and provide the optimal schedules for conventional generators to cope with the uncertainty in wind power generation. The proposed double-level model is difficult to be solved because of the invalid dual transform in solution process caused by the coupling relation between the discrete and continuous variables. Therefore, a two-stage iterative solution method based on Benders Decomposition is also presented. The proposed double-level model is transformed into a single-level and two-stage robust interval unit commitment model by eliminating the coupling relation, and then this two-stage model can be solved by Benders Decomposition iteratively. Simulation studies on a modified IEEE 26-generator reliability test system connected to a wind farm are conducted to verify the effectiveness and advantages of the proposed model and solution method.

#### 1. Introduction

At present, the variability, limited predictability, and antipeaking regulation inherent in wind power have created significant challenge to power systems operation with high wind power penetration [1, 2]. Unit commitment (UC), one of the most crucial processes in power systems schedule and operation, has been evolved from deterministic formulation into uncertainty formulation [3–6] to cope with uncertainties in wind power generation and load. In the stochastic UC model, the uncertainty of wind power is represented by numerous scenarios of possible wind power output which are often required to ensure the quality of the UC solution. Hence, the limitation with stochastic UC is that the UC problem size and computational requirement increased with the number of scenarios. Another drawback of stochastic UC is that the scenario generation method usually creates all scenarios based on certain probabilistic distribution assumption without specifying what scenarios can represent the ramp events. However, the probabilistic distribution of the uncertainty parameters is difficult to be acquired in its real world applications.

Robust optimization [7, 8] is an alternative uncertainty method to account for the uncertain parameters in optimization problems. Because the bounded uncertainty set considered in robust optimization method with a high solving efficiency is easy to be obtained in reality and the robust optimization strategies can consider the worst case to guarantee the security of power systems operation under all possible scenarios within a given uncertainty set, the robust optimization method applied to UC with uncertain wind power output has become the current research hotspots [9–14]. A contingency-constrained UC with an - security criterion based on robust optimization is proposed in [9], where a single-level equivalent robust counterpart of the original problem is formulated. Bertsimas et al. [10] and Jiang et al. [11], respectively, propose a two-stage adaptive robust UC model with adjustable uncertainty budget to reduce the conservativeness of conventional robust optimization. Hu et al. [12] proposes an effective method to acquire robust solutions to the security-constrained UC (SCUC) problem, which takes wind and load uncertainties into account via interval numbers. Based on [11], the uncertainty in demand response is involved in [13] and it assumes that the price-elastic demand curve is also varying within a given range to develop a multistage robust UC model. In [14], a new concept recourse cost defining the upper bound of redispatch cost when uncertainties are revealed is proposed to reduce the conservativeness of robust optimization. All these works [9–14] indicate that robust optimization can ensure system security under the worst-case scenario, and, therefore, it is an effective approach for solving optimization problems with uncertainties.

However, in the existing works [9–14], the robust optimization approaches applied to solve UC problem mainly focus on the adjustment of conservativeness of uncertain parameters including wind power output and load as well as demand response. Wind power spillage which is also called “wind power curtailment” is rarely considered in these robust UC models which have a premise that wind power generation can be completely accommodated no matter how much wind power output and its uncertainty are. As is known, in order to more accurately analyze the impacts of wind power accommodation on on/off schedules and spinning reserve requirements of conventional generators and consider the transmission capacity limits, it is necessary to treat the boundaries of wind power prediction interval as the optimizable decision variables and not merely uncertain parameters in robust UC models to determine the optimal wind power accommodation amount [15].

Therefore, based on a given wind power prediction interval, this paper first proposes a novel double-level robust SCUC problem formulation with optimizable interval of uncertain wind power output. In the proposed model, the boundaries of wind power prediction interval are treated as the optimizable decision variables. And the lower optimization model in the proposed double-level model contains the minimum or maximization formulation which represents the worst-case scenarios for the spinning reserve constraints and the transmission flow constraints. By shrinking the wind power prediction interval, this proposed model can curtail wind power to balance the dispatch cost of conventional generators and dispatch infeasibility penalty cost for wind farms, which is beneficial for mitigating the uncertainty of wind power, especially when the spinning reserve or the capacity of transmission lines is not enough. Furthermore, based on the optimal allowable wind power output interval sent from the operator of dispatch center, the wind farm could make the wind turbines be in the maximum power point tracking operation mode within allowable interval according to the actual wind speed condition to maximize the use of wind energy.

In addition, due to the invalid dual transform in solution process caused by the coupling relation between the discrete and continuous variables, such as the on/off schedules and power output of conventional generators, the proposed model is difficult to be solved and hence this paper proposes a two-stage iterative solution method based on Benders Decomposition (BD). Specifically, the discrete and continuous variables are regarded as the first and the second-stage decision variables, respectively, to eliminate the coupling relation between the discrete and continuous variables, and then the proposed double-level robust interval SCUC problem formulation can be transformed into a single-level and two-stage model which can be solved iteratively by BD algorithm.

The remainder of this paper is laid out as follows. First of all, Section 2 describes the double-level robust SCUC model with optimizable interval of uncertain wind power output. Then, the procedure of two-stage iterative solution method based on BD is given in Section 3. In Section 4, the case studies and simulation results analysis are presented. At last, in Section 5, main conclusions are summarized.

#### 2. Double-Level Robust SCUC Model with Optimizable Interval of Uncertain Wind Power Output

##### 2.1. Objective Function

The double-level robust interval SCUC model determines the allowable wind power generation interval from system security and economy points of view. The objective function involves two parts: the dispatch cost of conventional generators and the penalty cost of possible wind power spillage over all schedule periods. The dispatch cost contains the fuel cost, start-up and shut-down cost, and spinning reserve cost of generators. The penalty cost of wind power spillage can be expressed to be proportional to the magnitude of the difference between the wind power prediction interval and the allowable wind power output interval to maximize the wind power utilization [16]. Thus, the objective function of proposed UC model is defined as where is the number of time periods in the schedule horizon; the fuel cost of conventional generators is usually expressed as a piecewise linearization of quadratic convex function, and are the number of conventional generators and the segment number of piecewise linearized function, respectively;* NW* is the number of wind farms; and are cost slope of segment and minimum fuel cost, respectively, of generator , and meets ; is the power output of generator in segment during period ; denoting the on/off schedule of generator during period is a binary variable; and are up- and down-spinning reserve amounts of generator , respectively; and are the up- and down-spinning reserve cost coefficients of generator , respectively; and are the start-up and shut-down cost coefficients of generator , respectively; , , , and are the upper/lower limits of wind power prediction interval and the upper/lower limits of allowable wind power output interval for wind farm , respectively; is the penalty coefficient of the upper/lower limit deviation of wind power output interval for wind farm .

##### 2.2. Constraints

###### 2.2.1. Power Balance Constraint

where is the number of load buses; is the predicted load of bus during period ; and are the power output during period and the minimum power output of generator , respectively; is the most economic power output for wind farm during period .

###### 2.2.2. The Upper and Lower Limits for the Power Output of Conventional Generators

where , and are the power points, respectively, of linearized power output interval of generator .

###### 2.2.3. The Minimum On/Off Schedule Time Constraints of Conventional Generators

is introduced to show the on-time of generator needed to be maintained at the initial schedule period, and can be calculated aswhere is the minimum up-time of generator ; is the on-time of generator at the initial schedule period; denotes the on/off schedule of generator at the initial schedule period. If the generator needs shut-down, it must meet the minimum up-time as long as it is in on schedule, and hence this constraint can be listed as follows:where a binary variable is introduced to judge whether the generator is in the start-up process during period .

Similarly, is introduced to show the off-time of generator needed to be maintained at the initial schedule period, and can be calculated aswhere is the minimum down-time of generator ; is the off-time of generator at the initial schedule period. If the generator needs start-up, it must meet the minimum down-time as long as it is in off schedule, and hence this constraint can be listed as follows:where a binary variable is introduced to judge whether the generator is in the shut-down process during period . , , and satisfy

In addition, the start-up cost and shut-down cost of generator in the objective function (1) can be linearized to and , respectively.

###### 2.2.4. Spinning Reserve Constraints of Conventional Generators

where and are the ramp-up and ramp-down rates of generator , respectively; is the time resolution of per schedule period.

The spinning reserve constraint violation of conventional generator may cause wind power spillage. From the dynamic response ability of systems, wind power output mutation can lead to the decrease of adjustable capacity of generators. When the adjustable capacity of systems is smaller, the security level of systems is lower, and this scenario will be worse. Therefore, the worst-case scenarios for the spinning reserve constraints should be satisfied to guarantee the system security, which can be formulated as (10)–(13): where is the set of wind farms; and are up- and down-spinning reserve, respectively, supplied by systems under the worst-case scenario during period ; and are the power output of wind farm under the worst-case scenario for the up- and down-spinning reserve constraints, respectively, during period ; and are the minimum up- and down-spinning reserve requirements of systems, respectively, during period .

###### 2.2.5. Ramp-Rate Limits for Conventional Generators

From the dynamic response ability of systems, the ramp-rate limits for conventional generators under the worst-case scenario can be formulated as follows: where and are up and down power generation adjusting amounts for generator under the worst-case scenarios (10) and (12), respectively, during period . The generators may be in off schedule in UC during periods and ; hence it assumes in (18) and (19) that the power output of generator will reach the minimum value once the generator starts up and the power output of generator is required to be the minimum value before it shuts down.

###### 2.2.6. Transmission Flow Constraints

The transmission flow constraints violation may also lead to wind power spillage. From the network security, the wind power output boundary value will result in the maximum load rate of transmission line achieved. When the load rate of transmission line is higher, the security level of systems is lower, and this scenario will be worse. Therefore, the worst-case scenarios for the transmission flow constraints should be satisfied to guarantee the network security, which can be formulated as follows:

where is the number of buses; is the number of transmission lines; referring to [17] is power transfer distribution factor of bus to line ; , , and represent the generator , wind farm , and load connected to bus , respectively; and are the maximum positive and negative power flows of line during period ; and are the decision variables of power output of wind farm under the worst-case scenario for the positive and negative transmission flow constraint; is the capacity limit of line .

###### 2.2.7. Allowable Wind Power Output Interval Constraints

The upper limits of the maximum allowable wind power output interval must be lower than the upper limits of wind power prediction interval for every wind farm. Meantime, the lower limits of maximum allowable wind power output interval must be smaller than or equal to the lower limits of wind power prediction interval for practical generation schedule.

#### 3. Two-Stage Iterative Solution Method Based on BD Algorithm

Because the lower optimization models (10), (12), (20), and (21) contain the minimum or maximization formulation, the double-level robust interval SCUC model described in Section 2 cannot be solved efficiently. By the equivalent transformation of the minimum and maximization expressions based on linear duality theory [18], the double-level model needs to be transformed into a single-level model to be solved. However, due to the coupling relation between the continuous variables and the discrete variables of generator schedule in the upper and lower models, the upper model requirements for the lower model are not merely the maximum and minimum power generation capacity constraints. Therefore, the original duality transformation of the two-level model to a single-level model applied to economic dispatch will be invalid in the UC, and the corresponding double-level interval SCUC model is difficult to be solved.

For this reason, in this section, firstly, the discrete variables are regarded as the first-stage decision variables, and the continuous variables are regarded as the second-stage decision variables to eliminate the decouple relationship between the discrete and continuous variables. Then according to the linear duality theory, the double-level robust interval SCUC model described in Section 2 can be transformed into an equivalent single-level and two-stage robust interval SCUC model. Finally, the two-stage SCUC model can be solved iteratively by BD algorithm.

##### 3.1. Construction of a Single-Level and Two-Stage Robust Interval SCUC Model

In order to facilitate the description of the establishment process of the two-stage robust interval SCUC model and the solution method based on BD algorithm, firstly, the double-level model in a compact form is described as follows:where (23) is the objective function of the upper model, where is the column vector of decision variables composed of the discrete variables , , and , is the column vector of decision variables composed of the continuous variables , , , , , , , , , , , and , the coefficient column vectors and are acquired based on (1); the constraint conditions and of the upper model represent (5) and (7)-(8), represents (16), (17), and (22), represents (3), (9), (18), and (19), represents (2); the extremal problem of the lower model represents (10), (12), (14), and (15), represents (10)–(13), (20), and (21), , and are the coefficient matrixes, and , and are the coefficient column vectors.

The specific construction process of the single-level and two-stage robust interval SCUC model can be described as follows.

*Step 1. *The discrete variable vector is regarded as the first-stage decision variable, and the continuous variable vector is regarded as the second-stage decision variable. Then the double-level models (23) and (24) can be equivalent to the double-level and two-stage models (25)-(26)where

*Step 2. *In the second-stage model is a given vector, and can be regarded as a double-level robust interval economic dispatch model which cannot be solved directly due to the coupling relationship between the upper and lower model. However, according to a linear duality theory, from the goal of minimizing in the upper model, the constrained extremal problem in the lower model can be equivalent to the duality problem of the lower model. Hence, the double-level optimization model can be transformed to a single-level model to be solved.

In the constraint conditions of the lower models (10), (12), (20), and (21), the left-hand and right-hand constraints of , , , and are all decision variables, and the direct duality transformation may cause a nonlinear term of the product of two decision variables. Hence, these constraints need firstly simplification treatment. Take , for example, introduce a new continuous variable (), and let . Accordingly, , , and are substituted by , , and , respectively.

The lower optimization models (10), (12), (20), and (21) are, respectively, substituted by the corresponding duality problem to ensure that the dual objective function of the lower optimization model is the upper or lower bound of the original optimization model. Before the duality transformation, the decision variables of the lower optimization model are , , and . Correspondingly, after the duality transformation the duality variables are , , , and , respectively, and the duality model can be formulated asThe general form of this transformation based on linear duality theory is given in Appendix.

*Step 3. * and in (14) and (15) are substituted by the left side term in (27) and (28), respectively. Then the second-stage and double-level optimization model can be transformed to a single-level model , where satisfieswhere is the column vector of decision variables composed of the continuous variables , , , , , , , , , , , and ; denotes the second constraints of (16), (17), (22), and (27)–(30); denotes the first constraints of (3), (9), (18), (19), and (27)-(30); denotes (2), (14), and (15); , , , , and are coefficient matrixes; , , and are coefficient column vectors.

Finally the single-level and two-stage robust interval SCUC model equivalent to the original double-level models (23) and (24) can be obtained:

##### 3.2. Procedure of the Two-Stage Iterative Solution Method Based on BD Algorithm

The duality subproblem of the second-stage optimal subproblem in the two-stage robust interval SCUC model (32) iswhere the variables , , and are the duality variables of the variable . The feasible solution set (34) of the duality subproblem is independent of the discrete variable , and, based on the duality theory, the two-stage robust interval SCUC model can be expressed as

The pole set and the polar direction set of the feasible solution set (34) for the Benders subproblem are and , respectively. Introduce the relaxed variable , and construct the Benders main problem of the two-stage robust interval SCUC model (35) which can be described as follows: where (38) and (39) are the Benders optimality and feasibility cutting planes, respectively, connecting the master problem and subproblem.

Referring to [19, 20], the flow chart of the solution method based on BD is shown in Figure 1. Initialization firstly: set the iteration count , the given initial upper bound and lower bound , the convergence tolerance , = , and = ; construct a feasible solution . The detailed procedure of the solution method based on BD can be listed as follows.