Next Article in Journal
Approximation of Permanent Magnet Motor Flux Distribution by Partially Informed Neural Networks
Next Article in Special Issue
Integrated Carbon-Capture-Based Low-Carbon Economic Dispatch of Power Systems Based on EEMD-LSTM-SVR Wind Power Forecasting
Previous Article in Journal
Carbon Accounting Measurement with Digital Non-Financial Corporate Reporting and a Comparison to European Automotive Companies Statements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Distributionally Robust Chance-Constrained Unit Commitment with N-1 Security and Renewable Generation

Engineering Research Center of Education Ministry for Renewable Energy Power Generation and Grid Control, Xinjiang University, Urumqi 830047, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(18), 5618; https://doi.org/10.3390/en14185618
Submission received: 26 July 2021 / Revised: 27 August 2021 / Accepted: 27 August 2021 / Published: 7 September 2021

Abstract

:
With the increasing penetration of renewable energy generation, one of the major challenges is the problem of how to express the stochastic process of wind power and photovoltaic output as the exact probability density and distribution, in order to improve the security and accuracy of unit commitment results, a distributed robust security-constrained optimization model based on moment uncertainty is proposed, in which the uncertainty of wind and photovoltaic power is captured by two uncertain sets of first- and second-order moments, respectively. The two sets contain the probability distribution of the forecast error of the wind and photovoltaic power, and in the model, the energy storage is considered. In order to solve the model effectively, firstly, based on the traditional chance-constrained second-order cone transformation, according to the first- and second-order moments polyhedron expression of the distribution set, a cutting plane method is proposed to solve the distributed robust chance constraints. Secondly, the modified IEEE-RTS 24 bus system is selected to establish a simulation example, an improved generalized Benders decomposition algorithm is developed to solve the model to optimality. The results show that the unit commitment results with different emphasis on economy and security can be obtained by setting different conservative coefficients and confidence levels and, then, provide a reasonable decision-making basis for dispatching operation.

1. Introduction

In recent years, the rapid growth of renewable energy represented by wind and photovoltaic power in the energy system has greatly reduced the consumption of fossil fuels and greenhouse gas emissions [1]. However, its uncertainty and intermittence have a great impact on the security and stable operation of the power system. Therefore, how to solve the problem of unit commitment with large scale wind and photovoltaic power uncertainty has become a research hotspot.
In many unit commitment (UC) models considering the uncertainty of renewable generation, the traditional UC adopts a simple and feasible standby capacity method to deal with the uncertainty of renewable generation. This method is simple and fast, but it is blind and arbitrary. As a common method to deal with data uncertainty in optimization problems, stochastic programming has been widely explored in the past decade [2,3,4,5,6]. It mainly includes chance constraint and conditional risk method. This kind of method not only has large calculation scale, but also needs to know the deterministic distribution of random variables. Another tool for optimization against uncertainty is robust optimization. This method seeks the optimal objective values of all possible worst-case scenarios in the uncertainty set. In references [7,8,9,10], a robust optimization model for different objective functions is proposed based on the modeling of uncertain load and renewable energy through an adjustable robust optimization framework. However, in view of the fact that the worst case involved does not always happen, the optimization results are too conservative or aggressive. As a compromise method, distributed robust optimization can be used to solve the limitations of stochastic programming and robust optimization. In reference [11], a moment-based distributed robust model for UC is established, and the linear decision rule approximation method is applied to solve the model. However, only the uncertainty of wind power is considered, the impact of energy storage is not considered. In reference [12], the occurrence of accidents is regarded as a random event, and the moment uncertainty robust optimization method of unit commitmentconstrained by emergency accidents is studied without considering the uncertainty of renewable energy. The distribution set of the first- and second-order moment of wind power forecast error is given in reference [13], and the model is transformed into a bilinear matrix inequality problem to solve. Reference [14] obtained the first- and second-order moment information of wind power forecast error based on historical experience data and, then, transformed the optimal power flow model into convex second-order cone programming to solve the problem. In the above methods [11,12,13,14], DC power flow is used to replace AC power flow approximately, without considering the influence of reactive power and node voltage on the optimization results. The operation strategy obtained on this basis may lead to the risk of under voltage or over voltage on some key buses of the system. Reference [15] proposes a stochastic unit commitment (SUC) model to avoid the load loss risk caused by wind power uncertainty. In reference [16], the uncertainty of wind power is represented by the scenario tree model. Such a scenario tree is constructed by considering several stages, and each time period has several hours. The uncertainty of wind power generation in each cycle is represented by some discrete branches. This technology is suitable for modeling the recourse behavior of generator and energy storage device when the system uncertainty is realized successively. Reference [17] proposes a transmission capacity margin assessment (TCMA) model. Using dual theory and linearization technology, the model is transformed into mixed integer linear programming (MILP). In reference [18], the control characteristics of wind power plants is summarized, and the ability of wind power plants to provide fast regulation and frequency response services was determined.
Considering that the prediction accuracy of wind and photovoltaic power cannot accurately describe its probability, stochastic or robust optimization, such as widely used approximate normal distribution [19], beta distribution [20], mixed Laplace distribution [21], Cauchy distribution [22], α-steady-state distribution [23], etc., are essentially approximate methods, which is difficult to maintain the effectiveness of decision making. Distributed robust optimization combines robust optimization and stochastic optimization, two excellent methods with risk decision making, to implement the scheduling strategy under the probability distribution scenario of the worst uncertainty in the fuzzy set. For example, distributed robust unit commitment decision [24,25], reserve scheduling [26], and optimal power flow [27] can improve the effectiveness of decision making. The construction form of fuzzy sets determines different conservatism and solution efficiency. At present, the commonly used distribution robust optimization method based on probability density function is based on the probability uncertainty of Monte Carlo sampling points [28], which retains the characteristics that the probability density function of long-term fitting is similar to that of short-term fitting. However, this method does not make full use of the moment of random variables, and the variable scale is huge, it is difficult to obtain its dual problem. In a long statistical period, the prediction error of wind and photovoltaic power follows a certain distribution [29,30]. However, in a short period of time, its distribution is similar to the determined distribution, but it is not exactly the same, that is, the moment of the probability distribution is also uncertain. The distributed robust optimization under moment uncertainty (DRO-MU) method recently studied in the field of mathematics [31] can consider the similar but different characteristics of long-term fitting of random variables [32].
Given the context above, a robust optimization model of UC distribution with AC power flow constraints based on uncertainty of wind and photovoltaic power forecast error moment is proposed. Considering the network security constraints, the coupling relationship among active power, reactive power, voltage, and phase angle is considered, and the collaborative optimization is carried out. By setting conservative coefficient and confidence level, the optimal UC results with different emphasis on economy and security can be obtained, which can provide reasonable decision-making basis for dispatching operation. the main contributions of this paper are as follows: we develop a moment-based distributionally robust chance-constrained unit commitment with N-1 security and renewable generation model, while considering AC optimal power flow (OPF) and energy storage constraints; by considering an ambiguity set including the first- and second-order moment information obtained from empirical data, the wind and solar generation uncertainties are properly modeled.
The remainder of this paper is structured as follows. Section 2 first provides preliminaries on expansion time horizon and uncertainty modeling and, then, presents the proposed distributionally robust chance-constrained unit commitment (DRCCUC) model. Section 3 explains the solution methodology. Section 4 provides a comprehensive case study based on the IEEE RTS-24 node test system. Section 5 concludes the paper.

2. Distributed Robust Unit Commitment Model

2.1. Wind and Photovoltaic Fluctuations

In this paper, the active power output of each wind farm and photovoltaic power station is regarded as the sum of forecast value and forecast error value. The forecast value is regarded as the determined variable, and the forecast error value is treated as a random variable. It is assumed that the forecast error of all wind and photovoltaic power follow the normal distribution with zero mean value and are independent of each other [14]. For each wind farm, the hourly power is expressed as P b , t W + ω b , t W , where P b , t W is the forecast active power of wind farm at bus b at time t, and ω b , t W is the forecast error of active power of wind farm at bus b at time t. For each photovoltaic power station, the hourly output is expressed as P b , t PV + ω b , t PV , where P b , t PV is the forecast active power of photovoltaic power station at bus b at time t, and ω b , t PV is the forecast error of active power of photovoltaic power station at bus b at time t. As the security operation of power system requires that the power generation and power consumption should always be balanced, any deviation of wind and photovoltaic power from the forecast value must be balanced through the adjustment of controllable power generation. Since the power injection of each bus is fluctuating, a control scheme is needed to ensure that the power generation is always equal to the demand within a given time window. In this paper, it is assumed that all the generators participating in the control respond in a proportional manner to the fluctuation of the generalized load (assuming the forecast actual load minus the random wind and photovoltaic power), but the proportional coefficients may be different. Therefore, we call it initial value joint result frequency control and quadratic frequency control affine control. In this paper, the adjustment of wind and photovoltaic power fluctuation is modeled as an affine control strategy, and the activation of generation reserve is controlled by automatic generation control (AGC). The power balance is established in tens of seconds by AGC as follows:
P ^ i , t = P i , t α i , t Ω t t T ,   i G
Ω   t = Ω   t W   + Ω   t PV ;   Ω   t W   = w = 1 W ω w , t W ;   Ω   t PV = p v = 1 P V ω p v , t PV
where G is the number of generators; W is the number of wind farms; P V is the number of photovoltaic power stations; T is the scheduling cycle and the time interval is 1 h; P i , t is active power of generator; α i , t is participation coefficient of AGC controlled generator; Ω   t is the sum of forecast errors of output power of all wind farms and photovoltaic power stations; Ω   t W   is all wind farm power forecast error; Ω   t PV is all power forecast error of photovoltaic power station. When α i , t is constant or equal in the scheduling period, the calculation processability of the unit commitment model will be significantly improved. However, it is taken as an optimization variable in this paper for economically advantageous. In practical application, because of the randomness and inaccuracy of wind speed, usually only part of the probability distribution information can be obtained. Therefore, this paper uses two uncertainty intervals to describe the uncertainty of wind and photovoltaic power forecast error, μ ¯ b w , μ ¯ b w and μ ¯ b pv , μ ¯ b pv represent the mean change interval of wind and photovoltaic power forecast error at the bus b, and ( σ b w ) 2 ( σ ¯ b w ) 2 , ( σ b w ) 2 + ( σ ¯ b w ) 2 and ( σ b pv ) 2 ( σ ¯ b pv ) 2 , ( σ b pv ) 2 + ( σ ¯ b pv ) 2 represent the variance change interval of wind and photovoltaic power forecast error at the bus and assume that the random variables of wind and photovoltaic power forecast error are independent of each other. In fact, there are several probability distribution functions whose mean value and variance are consistent with these information, then the fluctuation range of mean value and variance of wind and photovoltaic power is represented by polyhedron set as follows:
U μ = μ W + P V : μ b w + μ b pv μ ¯ b w + μ ¯ b pv b W + P V μ b w + μ b pv μ ¯ b w + μ ¯ b pv Γ μ W + P V
U σ 2 = σ 2 W + P V ( σ b w ) 2 + ( σ b pv ) 2 ( σ ¯ b w ) 2 + ( σ ¯ b pv ) 2 b W + P V ( σ b w ) 2 + ( σ b pv ) 2 ( σ ¯ b w ) 2 + ( σ ¯ b pv ) 2 Γ σ W + P V
where μ b w and μ b pv are the mean variation of wind and photovoltaic power forecast error at bus b, respectively; μ ¯ b w and μ ¯ b pv are the upper limit of mean value change in wind and photovoltaic power forecast error at bus b, respectively; Γ μ is the mean interval conservative coefficient (the range limiting parameter of uncertainty set); σ ¯ b w and σ ¯ b pv are the upper limit of forecast error variance of wind and photovoltaic power at bus b, respectively; Γ σ is the conservative coefficient of variance interval.

2.2. DRCCUC Model

The objective function of DRCCUC minimizes the operation cost of generator, including no-load cost, start-up cost, shutdown cost, fuel cost, spinning reserve cost, and emergency reserve cost.
min i = 1 G t = 1 T { c i 0 x i , t + k K c i k P i , t k + C i , t U + C i , t D + [ c i r ( R i , t + + R i , t ) + c i * R i , t * ] + w = 1 W t = 1 T C w , t W + p v = 1 P V t = 1 T C p v , t PV }
where c i 0 is no load cost of generator; x i , t is the binary variable of generator start and stop state; c i k is the section k (k = 1, 2, 3) rate for the unit fuel cost function; P i , t k is the active power of section k for the fuel cost function of generator; C i , t U is the cost of generator start-up; C i , t D is the cost of generator shutdown; c i r is the spinning reserve cost for generators; R i , t + is the positive spinning reserve of generator; R i , t is the negative spinning reserve of generator; c i * is the contingency spinning reserve cost for generator accident; R i , t * is the contingency spinning reserve of generator. C w , t W is the penalty cost of wind power spillage; C p v , t PV is the penalty cost of photovoltaic power spillage. The main constraints are formulated as follows:
y i , t z i , t = x i , t x i , t 1 t T , i G
y i , t + z i , t 1   t T , i G
C i , t U C i U x i , t x i , t 1   and   C i , t U 0 t T , i G
C i , t D C i D x i , t 1 x i , t   and   C i , t D 0 t T , i G
P i , t R i , t P i min x i , t t T , i G
P i , t + R i , t + + R i , t * P i max x i , t t T , i G
i = 1 G α i , t = 1 ,   0 α i , t x i , t i G , t T  
i = G r α i , t = δ r , t = θ r , t = 0 t T
b B b r Y b , b δ b , t = i G b α i , t b B , t T
P i , t g c = P i , t + P ^ i , t g c t T ,   i , g c G
P ^ i , t g c R i , t * ,   i = 1 G P ^ i , t g c = 0 ,   P ^ g c , t g c = P g c , t t T ;   i , g c G
Pr W , P V ~ U μ , σ 2 R i , t Ω t α i t 1 ε i G μ U μ , σ 2 U σ 2
Pr W , P V ~ U μ , σ 2 R i , t + Ω t α i t 1 ε i G μ U μ , σ 2 U σ 2
where y i , t is the binary variable of generator starting state; z i , t is a binary variable of generator shutdown state; C i U is the cost of generator start-up; C i D is the Cost of generator shutdown; Constraint (6) according to the start-up and stop state between t hour and t − 1 determines whether the generator is in the start-up or stop state at t hour; Constraint (7) ensures that start-up and shut-down of any generator will not occur at the same time. P i min and P i max are the minimum and maximum active power of generator, respectively; θ r , t is the reference bus voltage phase angle; δ r , t is an auxiliary variable; P ^ i , t g c is the increased active power of generator in case of generator fault. In order to ensure the power balance, Constraint (16) forces the generator to stop unexpectedly, and the power output change is equal to its pre accident power output; Pr ( ) represents the probability of occurrence of random events; ε i G indicates the probability of insufficient reserve capacity of generator; Constraints (10) and (11) constrain the positive and negative rotating reserve and emergency reserve capacity of the generator; Constraints (17) and (18) ensure that the generator emergency reserve capacity is greater than the expected reserve activation. The combination of Constraints (11) and (16) can ensure that the total accident reserve is sufficient to cover any generator shutdown. Constraints (17) and (18) are robust chance constraints on the distribution of generator reserve capacity. Constraint (12) limits the number of generators participating in the automatic generation control based on the on–off state.
x i , t x i , 0 t t i on + t i off ,   i G
t = t ¯ t y i , t x i , t t t i on ,   t ¯ = t t i min on + 1 ,   i G
t = t _ t z i , t 1 x i , t t t i off ,   t _ = t t i min off + 1 ,   i G
P i down P i , t 1 P i , t t T ,   i G
P i up P i , t P i , t 1 t T ,   i G
Constraints (19)–(23) are the minimum start-up and shutdown time and climbing constraints of the generator, where t i on is the accumulated running time of the generator from t = 0; t i off is the accumulated shutdown time of generator from t = 0; t i min on is the minimum allowable operation time of generator; t i min off is the minimum allowable shutdown time of generator; P i up is the generator active output rise rate limit; P i down is the generator active power output reduction rate limit. Equation (19) constrains the switching state of the generator under initial conditions; Equations (20) and (21) are the constraints of minimum starting time and minimum stopping time of generator; Equations (22) and (23) are generator ramp constraints.
F m , n , t = β m , n ( θ m , t θ n , t ) ,   F m , n , t F m , n MAX { m , n } L
b B Y b , b θ b , t = i G b P i , t + P b , t W + ω b , t w + P b , t PV + ω b , t pv P b , t D + P b , t DIS P b , t CH b B , t T
Pr W , P V ~ U μ , σ 2 ( F m , n , t + β m , n Ω t ( δ n , t δ m , t ) + β m , n ω ^ t W ( π m , t π n , t ) + β m , n ω ^ t PV ( π m , t π n , t ) F m , n MAX ) 1 ε l L { m , n } L , μ U μ , σ 2 U σ 2
Pr W , P V ~ U μ , σ 2 ( F m , n , t + β m , n Ω t ( δ n , t δ m , t ) + β m , n ω ^ t W ( π m , t π n , t ) + β m , n ω ^ t PV ( π m , t π n , t ) F m , n MAX ) 1 ε l L { m , n } L , μ U μ , σ 2 U σ 2  
Equations (24)–(27) are DC power flow constraints under N state, where F m , n , t is the active power flow of transmission line; β m , n is the susceptance of the line; F m , n MAX is the maximum active power flow of transmission line; P b , t DIS is the active power of energy storage discharge; P b , t CH is the charging active power for energy storage; π m , t is the mth row of the inverse of the admittance matrix; ε l L is the probability of power overrun in normal operation.
F m , n , t c = β m , n ( θ m , t c θ n , t c ) { m , n } L
b B Y b , b c θ b , t c = i G b ( P i , t + δ i , t c ) + P b , t W + ω b , t w + P b , t PV + ω b , t pv P b , t D + P b , t DIS P b , t CH b B , t T
Pr W , P V ~ U μ , σ 2 ( F m , n , t c F m , n MAX ) 1 ε l LC   { m , n } L , μ U μ , σ 2 U σ 2  
Pr W , P V ~ U μ , σ 2 ( F m , n , t c F m , n MAX ) 1 ε l LC   { m , n } L , μ U μ , σ 2 U σ 2  
Equations (28)–(31) are DC power flow constraints under N-1 state of generator and transmission line, where F m , n , t c is the active power flow of transmission line in case of accident; θ m , t c is the bus phase angle in case of accident; ε l LC is the probability of power overrun of transmission line in case of accident. Equations (26), (27), (30), and (31), respectively, represent the robust chance constraints of line power flow distribution under normal operation and accident conditions.
P b , t G + P b , t W + P b , t PV + P b , t DIS P b , t CH P b , t D = b b V b , t V b , t G b b cos θ b b , t + B b b sin θ b b , t b , b B
Q b , t G + Q b , t W + Q b , t PV Q b , t D = b b V b , t V b , t G b b sin θ b b , t B b b cos θ b b , t b , b B
V b MIN V b , t V b MAX b B ,   t T
Q i MIN Q i , t Q i MAX i G ,   t T
F m , n , t P 2 + F m , n , t Q 2 F m , n , t S 2 { m , n } L , t T
where Q b , t G is the reactive power of thermal power unit; Q b , t W is the reactive power of wind power; Q b , t V is the photovoltaic reactive power; Q b , t D is the reactive power of the load; V b , t is the node voltage amplitude; G b b is the real part of the bus admittance matrix; B b b is the imaginary part of node admittance matrix; V b max is the upper limit of node voltage amplitude; V b min is the lower limit of node voltage amplitude; Q i , t is the reactive power output of thermal power units; Q i max is the upper limit of reactive power of thermal power unit; Q i min is the lower limits of reactive power of thermal power units; F m , n , t P and F m , n , t Q are the line active and reactive power; F m , n , t S is the line apparent power.
S e s , t CH + S e s , t DIS 1 e s E S , t T
E e s , t = E e s , t 1 + P e s , t CH η + P e s , t DIS / η t T , e s E S
0 η P e s , t CH P e s MAX S e s , t CH t T , e s E S
0 P e s , t DIS / η P e s MAX 1 S e s , t CH t T ,   e s E S
0.1 E e s MAX E e s , t 0.9 E e s MAX t T ,   e s E S
where S e s , t CH is the binary variable of energy storage charging state; S e s , t DIS is the binary variable of energy storage discharge state; E e s , t is the charge capacity of energy storage; P e s , t CH and P e s , t DIS are energy storage charging and discharging power, respectively; η is the conversion efficiency of charging and discharging for energy storage; P e s MAX and E e s MAX are the maximum rated power and capacity of energy storage, respectively; Equation (37) is used to restrict the charging and discharging state of energy storage, Equation (38) is used to calculate the charging quantity of energy storage, Equations (39) and (40) restrict whether the energy storage can be charged and discharged at the same time, and Equation (41) restricts the upper and lower limits of charging quantity of energy storage.

3. Solution Methodology

Due to the nonconvex of distributed robust chance constraints and the strong nonlinearity of AC power flow equations, the large-scale mixed integer nonlinear programming (MINLP) model established in the previous section is difficult to solve directly by using existing commercial solvers. Among many UC solving algorithms, the general Benders decomposition (GBD) algorithm is widely used to solve MINLP problems. In order to improve the efficiency of traditional GBD, according to the basic idea of variable separation in reference [33], this paper proposes an improved GBD algorithm. Compared with traditional GBD, it is mainly reflected in the following two aspects: (1) an improved cut plane method for solving distributed robust chance constraints is proposed, which reduces the number of unnecessary iterations; (2) a network security constraint reduction method is proposed. Through a large number of simulation calculations, it is known that even if there are power flow violations of security constraints in N and N-1 states, it is not necessary to add them to the model, but only a small number of security constraints need to be added to the model. Once these key security constraints are satisfied, the network security constraints can be reduced; then, a large number of other network security constraints will be automatically satisfied.

3.1. Distributed Robust Chance Constraint

The processing and solving of robust chance constraints of nonconvex distribution is the key to the effective solution of the model. The traditional fixed distribution chance constraints can finally be transformed into convex second-order cone (SOC) constraints, which can be solved by commercial solvers. The specific conversion process can be referred to [34]. In the traditional chance-constrained reconstruction process, firstly, the wind power forecast error assumption is required to obey the normal distribution with known mean and variance; secondly, the existing commercial solvers are time consuming and inefficient for solving models with a large number of SOC constraints. For the distributed robust chance constraint in the form of Equation (42), it is required to keep a high probability in all possible distribution sets U. That is to say, for each distribution in a polyhedral uncertainty sets, there is a single convex constraint in the form of traditional chance constraint, so (42) represents a set of potential infinite convex constraint sets, and there is no known deterministic reconstruction expression, which leads to the existing commercial solvers cannot solve directly. Therefore, based on reference [35], this paper proposes an improved cutting-plane approach to deal with large-scale chance constraint cases. The specific solving steps are as follows: (1) to solve the linear relaxation Formula (43) of constraint Formula (42), whether the constraint is satisfied must be confirmed in each iteration; (2) by solving Equation (43), we can judge whether Equation (42) is satisfied, that is, there is a linear objective for each fixed and all internal maximization problems. Therefore, we can judge whether a given constraint (42) is satisfied by calculating a linear optimization. If the result of Equation (44) shows that Equation (43) is true, the algorithm will terminate. Otherwise, we can obtain the mean and variance sets that make Equation (43) not true. The formative Equation (44) linearizes the constraint and returns to the first iteration.
Pr z N u , τ y T z d 1 ε u , τ U
max μ U u y T u + max τ U τ y T τ y Φ 1 1 ε d
y T u * + Φ 1 1 ε y * T Σ * y * + Φ 1 1 ε y * T τ * y * y * T τ * y y * d

3.2. Model Decomposition Strategy

Model decomposition strategy: firstly, the model is relaxed to obtain the relaxed model without (26), (27), (30), and (31) stochastic power flow security constraints. Then, according to the main idea that discrete variables and complex continuous variables are treated separately by GBD algorithm, the relaxed model is decomposed into a main problem and two sub-problems, in which sub-problem one is used to verify AC power flow. The second sub problem is used to check the generator N-1 fault. After solving the relaxed main problem and obtaining the unit commitment decision, the main problem solution is usually not feasible because the main problem does not consider the constraints of power flow, reactive power, and network security. At this time, it needs to be brought into the sub-problem for AC power flow and generator N-1 verification, and the verification results are modified in the form of benders cut. After the iteration of the main sub-problem is completed, whether the stochastic power flow of the transmission line fully meets the security constraints is checked. If there is any unsatisfied situation, the unsatisfied stochastic power flow constraints of the transmission line are added to the relaxed main problem, and the above solving process is repeated until all the constraints are met and the optimization results are output.
The Models (1)–(41) are decomposed into the following main unit commitment problems and sub-problems. Formula (45) is the main problem. For the sake of simplicity, the constraint conditions are all replaced by formula numbers.
o b j . min i = 1 G t = 1 T { c i 0 x i , t + k K c i k P i , t k + C i , t U + C i , t D + [ c i r ( R i , t + + R i , t ) ] + w = 1 W t = 1 T C w , t W + p v = 1 P V t = 1 T C p v , t PV } + ψ t s . t . 1 4 , 6 10 , 12 14 , 17 25 , 37 41
It can be seen from Equation (45) that the decomposed main problem contains all binary variables, and the constraints are linear or can be converted to SOC, so it is a typical MISOC programming problem. The last two terms of its objective function are obtained from benders cut set. Because there are two sub-problems in the proposed decomposition strategy, two alternative variables are set.
o b j . min b = 1 B Δ P b , t + + Δ P b , t + Δ Q b , t + + Δ Q b , t s . t . P b , t G + P b , t W + P b , t PV + P b , t DIS P b , t CH P b , t D + Δ P b , t + Δ P b , t = P b Q b , t G + Q b , t W + Q b , t V Q b , t D + Δ Q b , t + Δ Q b , t = Q b P i , t = P i , t ( k ) x i , t = x i , t ( k ) 32 36 Δ P b , t + 0 , Δ P b , t 0 , Δ Q b , t + , Δ Q b , t 0
The first subproblem is shown in Equation (46), where P b and Q b are the simplified form of power flow; P i , t ( k ) and x i , t ( k ) are the solution of the main problem of the kth iteration, which is treated as a constant in this paper. Δ P b , t + , Δ P b , t , Δ Q b , t + , and Δ Q b , t are the active and reactive relaxation variables introduced for bus b to ensure that the sub-problem always has a solution. It can be seen that the power flow constraint sub-problem has no coupling constraints such as climbing, so it can be directly decomposed into T smaller nonlinear programming sub-problems. Each sub-problem can be solved independently to improve the computational efficiency of the original problem. Take the sub-problem corresponding to time period t as an example to analyze: when the objective function value of the sub-problem is less than a certain threshold, for example 10 4 , the solution of the main problem is considered to meet the power flow and other operation constraints under time period t, the sub-problem is said to be feasible; otherwise, it is called infeasible. For the infeasible sub-problem, it is necessary to feed back the constraint information to the main problem in order to modify the solution of the main problem in the corresponding period. The key to the correction is benders cut.
b B Δ P b , t + + Δ P b , t + Δ Q b , t + + Δ Q b , t ( k ) + i G λ   i , t 1 ( k ) P i , t P i , t ( k ) + i G λ   i , t 2 ( k ) x i , t x i , t ( k ) 0
where λ   i , t 1 and λ   i , t 2 are Lagrange multipliers with k marked equality constraints in Equation (46), and Equation (47) is benders cut of infeasible sub-problem for time period t. The expression contains the information of constraint overrun and Lagrange multiplier, which can be regarded as the linear approximation of the sub-problem at the solution of the main problem. If it is brought into the constraint set of Equation (45), the solution of the main problem can be modified [36].
The second sub-problem is composed of the unit commitment decision and the generator output level solution obtained from the optimization of the main problem, which is responsible for calculating the new generation level for any generator outage. Each time period has a sub-problem. The objective of each sub-problem is to minimize the total reserve cost of all generator accidents.
o b j . min i = 1 G m i 2 R i , t * s . t . ( 11 ) , ( 15 ) , ( 16 ) , ( 28 ) , ( 29 )
For the convenience of explanation and explanation, only the variables of the main problem are reserved on the right side of the sub-problem constraint. After the values of P i , t , x i , t , and R i , t + are obtained from the iterative solution of the main problem. According to the result of the sub-problem solution, the feasible cut or the optimal cut are generated, respectively. These cutting planes are added to the main problem, and then the main problem is solved again. The objective value of the main problem provides a lower bound, and the upper bound is calculated by using the solution of the sub-problem. When the difference between the upper and lower limits is within the set threshold, the inner loop of the algorithm is terminated.

3.3. Algorithm

The flow chart of the algorithm is shown in Figure 1, the main steps are as follows:
  • The algorithm parameters are initialized, the number of iterations is set to zero, the upper bound (UB) and the lower bound (LB) are set to positive and negative infinity, respectively;
  • The main problem (45) is solved with the sum of substitution variables being zero and benders optimal cut set and feasible cut set being empty, and the results of generator start-up, stop, active power and rotating reserve power, energy storage charging and discharging state and charging and discharging power, and wind power photovoltaic active power output are obtained;
  • The first sub-problem is solved with the data stream obtained from the main problem. If the optimization result is feasible, benders optimal cut is generated according to Equation (47); if the optimization result is not feasible, benders feasibility cut is generated according to Equation (48);
  • The second sub-problem is solved by the data stream obtained from the main problem. If the optimization result is not feasible, benders feasibility cut is generated according to the formula; if the optimization result is feasible, benders optimal cut is generated according to the formula;
  • If all the first and second sub-problems in the scheduling cycle are feasible, then go to the next step, otherwise go to the seventh step;
  • Calculate the value of the dual gap and judge whether it meets the threshold of the end of the algorithm iteration. If it meets, go to step 8. If not, go to the next step;
  • Add the optimal cut set and feasible cut set to the main problem to get the new data stream and, then, go to the third step.
  • Calculate whether the network security constraints satisfy the constraints (26), (27), (30), and (31) in N and N-1 states. If all the constraints are satisfied, the optimization results are output and the algorithm ends. Otherwise, the unsatisfied constraints are added to the main problem, and the second step is to solve it again.

4. Numerical Case Studies

In order to verify the correctness and feasibility of the above model and algorithm, IEEE RTS-24 system [37] is selected for test. The total installed capacity of this system is 3405 MW. All thermal generator parameters contained in the test system are shown in Table 1. In addition, six wind farms with a capacity of 260 MW are connected to buses 2, 3, 5, 8, 17, and 21 of the system, three photovoltaic power stations with capacity of 400 MW are connected to buses 6, 16, and 23, and six sodium flow battery energy storage power stations with rated capacity of 100 MW are connected to buses 1, 2, 3, 5, 6, and 10. See Figure 2 for the schematic diagram of the grid placement of the power units. Nuclear power units participate in peak load regulation operation of power grid according to G mode “15-1-7-1” and maximum 50% [38]. The peak time electricity price of the system is 125 USD/MWh; the average electricity price is 78 USD/MWh; the valley time tariff is 31 USD/MWh; the energy storage conversion efficiency is 80%, and the complete charge/discharge time is 8 h; the penalty cost of wind/photovoltaic is 40 USD/MWh; let mean interval conservative coefficient be equal to variance interval conservative coefficient; the confidence level of each robust chance constraint is equal. Julia0.6 is selected as the programming language, CPLEX12.7 is used as the main problem solver, and Ipopt3.1 is used to solve the nonlinear AC sub-problems. All the calculations are completed on a computer equipped with Intel Core i5-2310 2.9GHz.

4.1. The Results of This Model

Using this method to simulate the test system, the confidence level is 0.992, and the conservative coefficient is 0.6. The results of the next day 24-hour conventional unit start-up and shutdown are shown in Figure 3. It is easy to see that two nuclear power units, 22 and 23, and three large capacity thermal power units, 20, 21, and 30, start-up at all times to provide extra base load. Other thermal power units with smaller capacity start and stop according to the changes in load and wind and photovoltaic power output.
The hourly online generator capacity (OGC), wind, and photovoltaic power output, energy storage charge and discharge power, and spinning reserve power are shown in Figure 4. It can be seen that the energy storage charge is concentrated in the valley period, and the energy storage discharge is concentrated in the peak period, indicating the feasibility of the decision results. Furthermore, by changing the parameters of the model, sensitivity analysis was used to detect the influence of the parameters on the UC results.

4.2. Influence of Conservative Coefficient on Objection

The fixed confidence level is 0.99. When the conservative coefficient changes from 0.001 to 1.0, the change details of the total generation cost of the test system are shown in Table 2. In order to compare the influence of AC and DC power flow on UC results, the table also gives the UC optimization results based on DC.
It can be seen from the Table 2 that with the gradual increase in the conservative coefficient, the spinning reserve cost, the penalty cost, and the total power generation cost of the system gradually increase, resulting in the gradual increase in the objection value. The change trend of them with the conservative coefficient is shown in Figure 5. It can be seen that when the conservative coefficient increases from 0.001 to 1, the spinning reserve power increases significantly synchronously, with a cumulative increase of 50%; the active power of the generator increases slowly with a cumulative increase of 6%; the active power output of wind and photovoltaic decreased slowly with the increase in conservative coefficient, with a cumulative decrease of 9%. The reduction in wind and photovoltaic power output indicates that the proportion of wind and photovoltaic energy consumption is reduced, and the abandoned wind and photovoltaic power is increased. The reason is that the increase in the conservative coefficient means that the fluctuation range of the forecast error becomes larger, which indicates that the accuracy of the wind and photovoltaic power output forecast decreases. At this time, in order to cope with the uncertainty of the continuous increase in wind and photovoltaic power and meet the safety requirements of power system operation as far as possible, more spinning reserve capacity is bound to be reserved, which leads to the increase in the generation cost of conventional units. In addition, the accuracy of wind and photovoltaic output forecast is reduced, which will increase the power of abandoned wind and solar power and, then, increase the cost of abandoned wind and solar power. These two reasons lead to the increase in the objection value.

4.3. Influence of Confidence Level on Objection

The change in the objection value when the confidence level changes from 0.98 to 0.998 is shown in Table 3. As a comparison, the optimization results considering only DC power flow are given.
It can be seen from Table 3 that with the gradual increase in the confidence level, the spinning reserve cost, abandonment cost, and fuel cost of the system gradually increase, resulting in the gradual increase in the objection value. The change trend of them with the confidence level is shown in Figure 6. It can be seen that when the confidence level gradually increases from 0.98 to 0.998, the synchronous increase in spinning reserve is obvious, and the cumulative increase is 33%; the active power of the generator increases slowly with a cumulative increase of 7%; the output power of wind and photovoltaic decreased slowly with the increase in conservative coefficient, with a cumulative decrease of 10%. The increase in confidence level means the increase in risk level of generator output and line power flow out of limit. When the chance constraint meets higher constraint probability, its corresponding feasible region is narrowed, and the optimization in the smaller feasible region must be at the expense of certain economic objection value.

4.4. Influence of AC Power Flow on Optimization Results

Compared with the optimization results without considering AC power flow in Table 2 and Table 3, when the parameters are consistent, the UC objection value, fuel cost, and abandoned wind and solar energy cost increase slightly with considering AC power flow, and the target value in Table 2 increases by 0.27–0.88%; the increase in fuel cost is between 0.16% and 0.81%, and the increase in abandoned wind and solar energy cost is between 2.3% and 5.4%. In Table 3, the objection value increases between 0% and 0.9%; the increase in fuel cost is between 0% and 0.59%, and the increase in abandoned wind and solar energy cost is between 0% and 11%.
The UC model based on DC does not consider the influence of reactive power and node voltage on the optimization results, which may lead to reactive power imbalance, and then, cause the risk of under voltage or over voltage on some buses. From the optimization results in Table 2 and Table 3, it can be seen that only when the confidence level is greater than or equal to 0.996, the optimization results meet the reactive power balance. All bus node voltages are within the allowable fluctuation range. In order to investigate the influence of AC on start-up, shutdown, and node voltage in detail, the conservative coefficient of 0.6 and confidence level of 0.992 are selected to optimize the calculation and track the iterative process. The results of the first solution of the main problem are substituted into the AC check sub-problem, and all bus node voltages are shown in Figure 7.
It is easy to see that there are three times of node voltage overrun in the sixth period, in which the lower limit of bus 7 is 1.74%, the lower limit of bus 8 is 0.22%, and the higher limit of bus 13 is 0.21%. In the 12th period, there were three instances of node voltage overrun, in which the lower limit of bus 1 was 0.37%, the lower limit of bus 2 was 0.56%, and the lower limit of bus 5 was 0.95%. In the 13th period, there were three instances of node voltage overrun, in which the lower limit of bus 1 was 0.53%, the lower limit of bus 2 was 0.67%, and the lower limit of bus 5 was 1.2%. In the 14th period, there were three instances of node voltage overrun, in which the lower limit of bus 1 was 0.78%, the lower limit of bus 2 was 0.99%, and the lower limit of bus 5 was 1.5%. It is easy to see that most of the voltage overruns are that the node voltage is less than the minimum limit, and the maximum overrun of 1.74% occurs at bus 7 of period 6. Usually, the node voltage is too low because of the lack of reactive power supply, which indicates that it is necessary to increase the number of generators started in these periods when the AC power flow is not feasible to increase the reactive power supply. According to the generated benders cut, it is returned to the main problem for the second solution. The comparison of unit startup and shutdown obtained from the two main problems is shown in Figure 8.
Through the comparison, it can be found that the main changes are in six periods, 11 units are added, nine units are added in 12 and 13 periods, and two units are added in 14 periods. Because these units are added with reactive power, the voltage of out-of-limit node returns to the limit range. Finally, after 14 iterations, the voltage of all nodes in the whole scheduling cycle is within the limit. It can be seen that UC considering AC power flow constraints can comprehensively consider the impact of reactive power and node voltage on system power generation capacity and line power flow when making a unit start-up and shutdown plan and ensure node voltage safety by coordinating and adjusting the unit start-up and shutdown status and output in advance.

4.5. Comparison with Stochastic Unit Commitment Model

In order to investigate the feasibility and calculation efficiency of the proposed model, the chance-constrained unit commitment [39] is selected and compared with the proposed model. In order to make the comparison more real and effective, the AC verification sub-problem is added to the CCUC model. The calculation objection value of the model is shown in Figure 9. It can be seen from the figure that the conservative coefficient is infinitely close to 0 regardless of the confidence level. It means that the fluctuation range of mean and variance is infinitesimal and close to 0. At this time, the mean and variance are close to a fixed value. The model in this paper is similar to the traditional chance-constrained model. From the calculation results of the two models, it can be seen that the results calculated by this method are very close to the results obtained by traditional chance-constrained SOC transformation, and the error is between 0.027% and 0.58%. It is shown that the proposed cut algorithm has high accuracy.
The iteration and calculation time of the main and subproblem of the model in this paper are shown in Table 4.The results show that under different confidence levels, the number of iterations of the main sub-problem is slightly larger than that of the CCUC, and the solving time of the sub-problem is consistent with that of CCUC, but the total computing time is significantly lower than that of CCUC, After analysis, it is considered that this is mainly because the cutting plane algorithm proposed for distributed robust chance constraints greatly reduces the complexity of the main problem model and the solution time of the main problem. Furthermore, it shows that this algorithm has high computational efficiency and can meet the requirements of the actual system.

5. Conclusions

In order to solve the problem of insufficient historical data and unknown real probability distribution of wind and photovoltaic power output, this paper improves the traditional chance-constrained unit commitment model. Firstly, the general moment uncertainty distribution set of wind and photovoltaic power forecast error is selected, its polyhedral expression is given, and an improved cut plane algorithm is proposed; secondly, the AC power flow model, which truly reflects the physical law of power grid operation, is adopted to ensure the node voltage security and improve the accuracy of UC results. Through comparative analysis of simulation examples, the following conclusions are obtained.
(1) In this model, the uncertainty of renewable energy output of UC is transformed into a conservative coefficient, and the risk of generator output and line power flow out of limit is transformed into the confidence level. After selecting reasonable conservative coefficient and confidence level, the UC optimization results of this model can achieve the decision goal of minimizing operation cost under certain risk.
(2) Under the premise that the objection value of UC based on AC is less than 1.0%, the safety margin of bus voltage is significantly improved, and the risk of bus voltage overrun is reduced, so as to improve the safety of system operation.
(3) The simulation results show that the larger the range of uncertainty set, the greater the total scheduling cost of the system. Under different confidence levels, the optimal scheduling results focusing on economy or security are different. Compared with the chance-constrained method with specific distribution, the results show that the moment uncertain distribution robust optimization method has good computational accuracy and higher computational efficiency.

Author Contributions

Conceptualization, Q.S. and W.W.; methodology, H.W.; software, Q.S.; validation, Q.S., W.W. and H.W.; formal analysis, H.W.; investigation, Q.S.; resources, H.W.; data curation, Q.S.; writing—original draft preparation, Q.S.; writing—review and editing, Q.S.; visualization, H.W.; supervision, W.W.; project administration, W.W.; funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Chinese National Natural Science Foundation (No. 51667020), Innovation team project of the Ministry of Education of China (No. IRT-16R63), Open project of Key Laboratory of Xinjiang Autonomous Region (No. 2018D0400).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xia, Q.; Zhong, H.; Kang, C. Review and Prospects of the Security Constrained Unit Commitment Theory and Applications. Proc. CSEE 2013, 33, 94–103. [Google Scholar]
  2. Zhao, C.; Guan, Y. Data-Driven Stochastic Unit Commitment for Integrating Wind Generation. IEEE Trans. Power Syst. 2015, 31, 2587–2596. [Google Scholar] [CrossRef]
  3. Uckun, C.; Botterud, A.; Birge, J. An Improved Stochastic Unit Commitment Formulation to Accommodate Wind Uncertainty. IEEE Trans. Power Syst. 2015, 31, 2507–2517. [Google Scholar] [CrossRef]
  4. Wang, C.; Wei, H.; Wu, S. Stochastic-Security-Constrained Unit Commitment Considering Uncertainty of Wind Power. Power Syst. Technol. 2017, 41, 1419–1427. [Google Scholar]
  5. Wang, B.; Wang, S.; Zhou, X.; Watada, J. Multi-objective unit commitment with wind penetration and emission concerns under stochastic and fuzzy uncertainties. Energy 2016, 111, 18–31. [Google Scholar] [CrossRef]
  6. Wu, Z.; Zeng, P.; Zhang, X.-P.; Zhou, Q. A Solution to the Chance-Constrained Two-Stage Stochastic Program for Unit Commitment With Wind Energy Integration. IEEE Trans. Power Syst. 2016, 31, 4185–4196. [Google Scholar] [CrossRef]
  7. Bertsimas, D.; Litvinov, E.; Sun, X.A.; Zhao, J.; Zheng, T. Adaptive Robust Optimization for the Security Constrained Unit Commitment Problem. IEEE Trans. Power Syst. 2013, 28, 52–63. [Google Scholar] [CrossRef]
  8. Wu, W.; Wang, K.; Li, G. Affinely Adjustable Robust Unit Commitment Considering the Spatiotemporal Correlation of Wind Power. Proc. CSEE 2017, 37, 4089–4097. [Google Scholar]
  9. Xiong, P.; Jirutitijaroen, P. Two-stage adjustable robust optimization for unit commitment under uncertainty. IET Gener. Transm. Distrib. 2014, 8, 573–582. [Google Scholar] [CrossRef]
  10. Fan, L.; Wang, K.; Li, G.; Wu, W.; Ge, W. Robust Unit Commitment Considering Temporal Correlation of Wind Power. Autom. Electr. Power Syst. 2018, 42, 91–97. [Google Scholar]
  11. Xiong, P.; Jirutitijaroen, P.; Singh, C. A Distributionally Robust Optimization Model for Unit Commitment Considering Uncertain Wind Power Generation. IEEE Trans. Power Syst. 2017, 32, 39–49. [Google Scholar] [CrossRef]
  12. Zhao, C.; Jiang, R. Distributionally Robust Contingency-Constrained Unit Commitment. IEEE Trans. Power Syst. 2017, 33, 94–102. [Google Scholar] [CrossRef]
  13. Wang, Z.; Bian, Q.; Xin, H.; Gan, D. A distributionally robust co-ordinated reserve scheduling model considering CVaR-based wind power reserve requirements. IEEE Trans. Sustain. Energy 2016, 7, 625–636. [Google Scholar] [CrossRef]
  14. Xie, W.; Ahmed, S. Distributionally robust chance constrained optimal power flow with renewables: A conic reformulation. IEEE Trans. Power Syst. 2018, 33, 1860–1867. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Wang, J.; Ding, T.; Wang, X. Conditional value at risk-based stochastic unit commitment considering the uncertainty of wind power generation. IET Gener. Transm. Distrib. 2018, 12, 482–489. [Google Scholar] [CrossRef]
  16. Xiong, P.; Singh, C. Optimal Planning of Storage in Power Systems Integrated with Wind Power Generation. IEEE Trans. Sustain. Energy 2015, 7, 232–240. [Google Scholar] [CrossRef]
  17. Yang, Y.; Lei, X.; Zhai, Q.; Wu, J.; Guan, X. Transmission capacity margin assessment in power systems with uncertain wind integration. In Proceedings of the 2017 13th IEEE Conference on Automation Science and Engineering (CASE), Xi’an, China, 20–23 August 2017. [Google Scholar]
  18. Tarnowski, G.C.; Kjær, P.C.; Dalsgaard, S.; Nyborg, A. Regulation and frequency response service capability of modern wind power plants. In Proceedings of the IEEE Power & Energy Society General Meeting, Minneapolis, MN, USA, 25–29 July 2010. [Google Scholar]
  19. Lin, Y.; Yang, M.; Wan, C.; Wang, J.; Song, Y. A Multi-Model Combination Approach for Probabilistic Wind Power Forecasting. IEEE Trans. Sustain. Energy 2019, 10, 226–237. [Google Scholar] [CrossRef]
  20. Luig, A.; Bofinger, S.; Beyer, H.G. Analysis of confidence intervals for the prediction of regional wind power out-put. In Proceedings of the European Wind Energy Association Conference, Copenhagen, Denmark, 1 January 2001; pp. 725–728. [Google Scholar]
  21. Tewari, S.; Geyer, C.J.; Mohan, N. A Statistical Model for Wind Power Forecast Error and its Application to the Estimation of Penalties in Liberalized Markets. IEEE Trans. Power Syst. 2011, 26, 2031–2039. [Google Scholar] [CrossRef]
  22. Hodge, B.; Milligan, M. Wind power forecasting error distributions over multiple time scales. In Proceedings of the Power and Energy Society General Meeting, Detroit, MI, USA, 24–29 July 2011. [Google Scholar]
  23. Bruninx, K.; Delarue, E. A Statistical Description of the Error on Wind Power Forecasts for Probabilistic Reserve Sizing. IEEE Trans. Sustain. Energy 2014, 5, 995–1002. [Google Scholar] [CrossRef]
  24. Liu, M.; Zeng, C.; Miao, H. Unit commitment model based on distributionally robust chance constraints. Electr. Meas. & Instrum. 2021, 58, 32–36. [Google Scholar] [CrossRef]
  25. Shi, Y.; Wang, L.; Chen, W.; Guo, C. Distributed Robust Unit Commitment with Energy Storage Based on Forecasting Error Clusting of Wind Power. Autom. Electr. Power Syst. 2019, 43, 3–12. [Google Scholar] [CrossRef]
  26. Wei, W.; Liu, F.; Mei, S. Distributionally Robust Co-Optimization of Energy and Reserve Dispatch. IEEE Trans. Sustain. Energy 2015, 7, 289–300. [Google Scholar] [CrossRef]
  27. Zhang, Y.; Shen, S.; Mathieu, J.L. Distributionally robust chance-constrained optimal power flow with uncertain renewables and uncertain reserves provided by loads. IEEE Trans. Power Systems 2017, 32, 1378–1388. [Google Scholar] [CrossRef]
  28. Liu, S.; Jian, J.; Wang, Y.; Liang, J. A robust optimization approach to wind farm diversification. Int. J. Electr. Power Energy Syst. 2013, 53, 409–415. [Google Scholar] [CrossRef]
  29. Wang, S.; Li, B.; Yu, J.; Xu, T. Analysis on Time-Varying Characteristics of Probability Error in Forecast of Wind Speed and Wind Power. Power Syst. Technol. 2013, 37, 967–973. [Google Scholar]
  30. Xie, L.; Gu, Y.; Zhu, X.; Genton, M.G. Short-Term Spatio-Temporal Wind Power Forecast in Robust Look-ahead Power System Dispatch. IEEE Trans. Smart Grid 2014, 5, 511–520. [Google Scholar] [CrossRef]
  31. Delage, E.; Ye, Y. Distributionally robust optimization under moment uncertainty with application to data-driven problems. Oper. Res. 2010, 58, 595–612. [Google Scholar] [CrossRef] [Green Version]
  32. Cheng, J.; Lisser, A.; Letournel, M. Distributionally robust stochastic shortest path problem. Electron. Notes Discret. Math. 2013, 41, 511–518. [Google Scholar] [CrossRef]
  33. Benders, J.F. Partitioning procedures for solving mixed-variables programming problems. Numer. Math. 1962, 4, 238–252. [Google Scholar] [CrossRef]
  34. Sha, Q.; Wang, W. Optimal Allocation of Energy Storage Considering N-1 Security Constraints Based on Chance Constraints. Adv. Eng. Sci. 2019, 51, 147–156. [Google Scholar]
  35. Lubin, M.; Dvorkin, Y.; Backhaus, S. A Robust Approach to Chance Constrained Optimal Power Flow with Renewable Generation. IEEE Trans. Power Syst. 2015, 31, 3840–3849. [Google Scholar] [CrossRef] [Green Version]
  36. Wang, C.; Wei, H.; Wu, S. A Decomposition coordination Algorithm Applied to Hydro-thermal Unit Commitment Problems With Power Flow Constraints. Proc. CSEE 2017, 37, 3148–3161. [Google Scholar]
  37. Grigg, C.; Wong, P.; Albrecht, P.; Allan, R.; Bhavaraju, M.; Billinton, R.; Chen, Q.; Fong, C.; Haddad, S.; Kuruganty, S.; et al. The IEEE Reliability Test System-1996. A report prepared by the Reliability Test System Task Force of the Application of Probability Methods Subcommittee. IEEE Trans. Power Syst. 1999, 14, 1010–1020. [Google Scholar] [CrossRef]
  38. Zhao, J.; Liu, D.; Lei, Q.; Du, Z.; Wang, J.; Zhou, L.; Huang, Y. Analysis of nuclear power plant participating in peak load regulation of power grid and combined operation with pumped storage power plant. Proc. CSEE 2011, 31, 1–6. [Google Scholar]
  39. Wang, B.; Yang, X.; Short, T.; Yang, S. Chance constrained unit commitment considering comprehensive modelling of demand response resources. IET Renew. Power Gener. 2017, 11, 490–500. [Google Scholar] [CrossRef]
Figure 1. Flow chart of the proposed algorithm.
Figure 1. Flow chart of the proposed algorithm.
Energies 14 05618 g001
Figure 2. Schematic diagram of the grid.
Figure 2. Schematic diagram of the grid.
Energies 14 05618 g002
Figure 3. On and off status of thermal power units.
Figure 3. On and off status of thermal power units.
Energies 14 05618 g003
Figure 4. Status of source and load and storage.
Figure 4. Status of source and load and storage.
Energies 14 05618 g004
Figure 5. Variation of generator and wind power with different conservative coefficient.
Figure 5. Variation of generator and wind power with different conservative coefficient.
Energies 14 05618 g005
Figure 6. Variation of generator and wind power with different confidence level.
Figure 6. Variation of generator and wind power with different confidence level.
Energies 14 05618 g006
Figure 7. First iteration bus voltage.
Figure 7. First iteration bus voltage.
Energies 14 05618 g007
Figure 8. On and off status of generator.
Figure 8. On and off status of generator.
Energies 14 05618 g008
Figure 9. Comparison of model calculation results.
Figure 9. Comparison of model calculation results.
Energies 14 05618 g009
Table 1. Summary of all thermal power units contained in the test system.
Table 1. Summary of all thermal power units contained in the test system.
UnitNodeTechnology P i max
(MW)
P i min
(MW)
P i up
(MW/h)
P i down
(MW/h)
UT
(h)
DT
(h)
Noload Costs (USD) c i 1
(USD/MW)
c i 2
(USD/MW)
c i 3
(USD/MW)
1–21OCGT2089010021454.628.9729.2429.70
3–41CCGT764012012032263.418.4219.2320.11
5–62OCGT2089010021454.628.9729.2429.70
7–82CCGT764012012032263.418.4219.2320.11
9–117CCGT1001042042042306.617.5918.2818.96
12–1413CCGT19710431031043482.917.2117.7118.23
15–1915OCGT125.4607021365.529.4630.1330.86
2015IGCC15554.2470802416415.523.8124.5225.25
2116IGCC15554.2470802416415.523.8124.5225.25
2218Nuclear40010028028016824188.36.967.237.50
2321Nuclear40010028028016824188.36.967.237.50
24–2922CCGT502612012021626.128.3129.2530.49
30–3123IGCC15554.2470802416415.523.8124.5225.25
3223Coal35014014014085303.826.2126.7127.20
Table 2. Object value of different conservative coefficients.
Table 2. Object value of different conservative coefficients.
Γ σ , Γ μ Power FlowReserve Cost/USDPenalty Cost/USDFuel Cost/USDObjection/USD
0.001AC91,96813,595572,691627,738
0.2119,37744,340593,714710,645
0.4130,96162,423623,802767,520
0.6134,23666,844629,503780,846
0.8137,46675,979638,049802,159
1138,75084,398640,529814,397
0DC91,96813,223568,113625,064
0.2119,37743,340592,763708,703
0.4130,96159,610620,243760,830
0.6134,23663,431626,464774,388
0.8137,46673,222633,233798,928
1138,75080,420639,155809,511
Table 3. Object value of different confidence levels.
Table 3. Object value of different confidence levels.
ε i G , ε l L Power FlowReserve COST/USDPenalty Cost/USDFuel Cost/USDObjection/USD
0.02AC120,97454,096599,869723,844
0.01134,23666,844629,503780,846
0.008138,25378,032634,329 802,046
0.006143,27589,658642,196826,545
0.004150,083105,532657,015865,024
0.002161,082135,132677,079929,561
0.02DC120,97448,721598,148720,156
0.01134,23663,431626,464774,388
0.008138,25374,450630,581794,805
0.006143,27587,215640,130822,573
0.004150,083105,532657,015865,024
0.002161,082135,132677,079929,561
Table 4. Algorithm calculation efficiency comparison.
Table 4. Algorithm calculation efficiency comparison.
ε i G , ε l L AlgorithmIterationsMaster Problem/sSubproblem/sTotal/s
0.02DRSCUC712020980
0.01513520775
0.008415020680
0.0066165201110
0.0045180201000
0.0025210201150
0.02CCSCUC5300201600
0.013330201050
0.008238020800
0.0064440201840
0.0043520201620
0.0024600202480
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sha, Q.; Wang, W.; Wang, H. A Distributionally Robust Chance-Constrained Unit Commitment with N-1 Security and Renewable Generation. Energies 2021, 14, 5618. https://doi.org/10.3390/en14185618

AMA Style

Sha Q, Wang W, Wang H. A Distributionally Robust Chance-Constrained Unit Commitment with N-1 Security and Renewable Generation. Energies. 2021; 14(18):5618. https://doi.org/10.3390/en14185618

Chicago/Turabian Style

Sha, Qiangyi, Weiqing Wang, and Haiyun Wang. 2021. "A Distributionally Robust Chance-Constrained Unit Commitment with N-1 Security and Renewable Generation" Energies 14, no. 18: 5618. https://doi.org/10.3390/en14185618

APA Style

Sha, Q., Wang, W., & Wang, H. (2021). A Distributionally Robust Chance-Constrained Unit Commitment with N-1 Security and Renewable Generation. Energies, 14(18), 5618. https://doi.org/10.3390/en14185618

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop