Feasible Modified Subgradient Method for Solving the Thermal Unit Commitment Problem as a New Approach

The Lagrangian relaxationLRbased methods are commonly used to solve the thermal unit commitment UC problem which is an important subject in power system engineering. The main drawback of this group of methods is the difference between the dual and the primal solutions which gives some significant problems on the quality of the feasible solutions. In this paper, a new approach, feasible modified subgradient F-MSG method which does not require finding an unconstrained global minimum of the Lagrangian function and knowing an optimal value of the problem under consideration in order to update dual variables at the each iteration, is firstly used for solving the thermal UC problem. The major advantage of the proposed approach is that it guarantees the zero duality gap and convergence independently from the size of the problem. In order to discuss the advantages of this method, the four-unit Tuncbilek thermal plant, which is located in Kutahya region in Turkey, is chosen as a small test system. The numerical results show that F-MSG gives better solutions as compared to the standard LR method.


Introduction
Solving thermal UC problem is very important for power system optimization.The thermal UC problem, which is defined as a nonlinear, mixed-integer, and nonconvex, is considered to be an NP Nonlinear-Programming hard problem.The high dimension of the possible solution space is the real difficulty in solving the problem.There have been various methods for solving the thermal UC problem in the literature.They are based on mathematical programming and metaheuristic-based approaches.These methods are priority list 1-3 , dynamic programming 4 , simulated annealing 5-7 , tabu search 8, 9 , evolutionary programming 10, 11 , constraint logic programming 12 , expert systems 13 , fuzzy systems 14 , genetic algorithms 15-19 , LR 20-24 , interior point algorithm 25 , memetic algorithm 26 , neural network 27, 28 , and hybrid models 29, 30 .LR which is most widely used approach among the aforementioned methodologies solves the thermal UC problem by "relaxing" or temporarily ignoring the coupling constraints and solving the problem as if they did not exist.Based on the duality theory, the LR method subsequently tries to find the values of the Lagrange multipliers that maximize the dual objective function.The duality gap, a major problem in the nonlinear programming, has been long recognized as an inherent disadvantage of these methods.If the LR-based methods are used for solving the UC problem, dual solution may be far away from the optimal solution.The duality gap for the problem of thermal UC is an important measure of the quality of the solution, which is defined by the smaller gap, the better solution 31 .In 32 , different mathematical-based methods LR, penalty function, and augmented Lagrangian penalty function are compared to each other according to feasible cost, dual cost, duality gap, number of iterations, and duration time.According to 32 , it is seen that there are differences between primal value and dual value.In the nonlinear programming, the duality gap has been investigated and the theoretical tools for zero duality gap condition have been improved extensively in 33-38 .In this paper, one of the methods based on dual optimization technique, F-MSG method, which has the best performance in eliminating the duality gap in the literature, is firstly used for solving the thermal UC problem as a novel approach.A dual problem with respect to the sharp augmented Lagrangian is constructed for this problem.F-MSG method solves without any duality gap for large class of nonconvex constrained problems.This proposed approach is compared to LR method to show that F-MSG gives an optimal solution obtaining the zero duality gap value.The LR method chosen in testing of F-MSG methods' performance lies from the fact that LR method is the most widely used method in solving UC problems in the literature.In order to compare these methods, the UC problem for fourunit Tuncbilek thermal plant, which is in Kutahya region in Turkey, is solved using MATLAB and General Algebraic Modeling System GAMS .Note that there is a general acceptance that whenever the system size gets smaller, the duality gap value gets bigger.However, it is shown that this is not the case with the F-MSG method for a small size system 38 .
The remaining sections are outlined as follows.Section 2 addresses the thermal UC problem formulations.This section includes an objective function and constraints of the problem.In Section 3, usage of F-MSG method for solving the thermal UC problem is explained in detail.In Section 4, numerical results for F-MSG and LR methods are provided.Finally the conclusions are presented, the main contributions of this paper are highlighted, and possible research directions are discussed.

Review of Thermal Unit Commitment Problem
The objective of thermal UC problem is to minimize the production cost over the scheduled time horizon e.g., 8 hours, 24 hours under the generator operational and spinning reserve constraints while all generating unit constraints are satisfied.The objective function of the UC problem can be formulated as where

Mathematical Problems in Engineering 3
The minimization of the objective function is provided to the following constraints.
a Power Balance Constraint.For satisfying the load balance in each stage, the forecasted load demand should be equal to the total power of the generated power for feasible combination: b Spinning Reserve Constraints.In order to meet any abnormal operating conditions, the spinning reserve must be considered: c Generation Limits Constraint.Each unit must satisfy the generation range: d Ramp-Up and Ramp-Down Constraints.For each unit, output is limited by ramp up/down rate at each hour as follows:

2.5
e Minimum Up-Down Time Constraint.Once the unit is started up, it should not be shut down before a minimum up-time period is met as follows, and once the unit is shut down, there is a minimum down time before it can be started up again:

2.6
The problem is structured binary variables making it clearly nonconvex.These binary variables cause a great deal of trouble and difficulty in solving the UC.Load balance and spinning reserve constraints are coupling constraints for the thermal UC problem.

F-MSG Method for Thermal Unit Commitment
The F-MSG method guarantees a zero duality gap for a wide class of nonconvex optimization problems.It does not require finding an unconstrained global minimum of the Lagrangian function and knowing an optimal value of the problem under consideration in order to update dual variables at the each iteration, and it guarantees convergence.

F-MSG Algorithm
Minimize f 0 x subject to x ∈ S satisfying f x 0, where S is a subset of a metric space X, and f 0 : X → R and f : X → R p are given functions.
Sharp Augmented Lagrangian is defined as follows: and defining the dual function as The dual problem Using the above definitions, the F-MSG is as follows 38 .

F-MSG Algorithm
Step 1. Choose positive numbers ε 1 ε 2 , Δ 1 , M and a number H 1 .Set n 1, t 0, and q 0. Step Step 3. Given v k , c k , solve the following constraint satisfaction problem P H n : find an element x ∈ S such that If a solution to 3.4 does not exist, then go to Step 6; otherwise, if a solution x k exists, then check whether f x k 0 or not.If f x k 0 or f x ≤ ε 1 , then go to Step 5; otherwise go to Step 4.
Step 4. Update dual variables as

3.5
where s k is a positive step size parameter defined as with α > 0, and 0 < δ < 2. In this study, 3.7 is used for updating.The appropriate values are replaced for updating in this equation.We will also require the step size s k corresponding to the dual variables v k , c k to satisfy the following property: Set k k 1, update l k in such a way that l k → ∞ as k → ∞, and repeat Step 3.
Step 5. Let x k be a solution to 3.4 with f x k 0. In this case L x k , v k , c k f 0 x k .Set q q 1 and check t.If t 0, then set Δ n 1 Δ n ; otherwise set Δ n 1 1/2 Δ n .Check Δ n 1 < ε 2 , and then stop.f 0 x k is an approximate optimal value, x k is an approximate primal solution, and v k , c k is an approximate dual solution; otherwise set H n 1 min{f 0 x k , H n − Δ n 1 }, n n 1, and go to Step 2.
Step 6. Set t t 1.If q 0, then set

and go to Step 2.
According to F-MSG method for solving the UC problem, the objective function and constraints are defined as follows.

Object Function
Slack variables are added to inequality constraints.Equality and inequality constraints can be defined as follows.

3.13
Minimum Up and Down Constraints 3.14 According to 3.9 , sharp augmented Lagrangian is defined as follows:

3.15
The epsilon value is chosen as 0.95 z j in the thermal UC problem solution.After 3.15 is constructed using the objective function and constraints, all the steps of the F-MSG method are applied to the thermal UC problem.

Numerical Results
F-MSG method is applied to the system to find the feasible results and feasible cost, dual cost, and duality gap values.This proposed approach is compared to LR method to show that F-MSG gives an optimal solution obtaining the zero duality gap value.The data used in the analysis is taken from Turkish Electric Power Company and Electricity Generation Company.The system parameters and reserve requirements of four units are given in Table 1.
In this study, the 24-hour day is subdivided into 8 discrete stages.The load demands for the stages are given in Table 2.The spinning reserve demand is set as 10% of the system demand for each stage.LR method is coded in MATLAB, and F-MSG method is coded in GAMS which is a high-level modeling system for mathematical programming problems 40 .
MINOS is used as a GAMS solver.These computer programs are run on Pentium core-1.6GHz personal computer.Optimal solution of the dual model set up by the sharp augmented Lagrange function, which is proposed by Gasimov, is equal to the primal solution of the system 37, 38 .The F-MSG method is run until the norm is equal to 7.36 × 10 −12 , epsilon is equal to 1.72 × 10 7 , then all the constraints reduce to zero, and feasible solution is obtained.Primal value is equal to the dual value yielding zero duality gap value.The value for the c parameter is 27182.94and the UC schedule, feasible cost, dual cost, and duality gap values for the F-MSG method are given in Table 3.
The UC schedule, feasible cost, dual cost, and duality gap values for the LR method are given in Table 4.
Feasible cost, dual cost, and duality gap values are found by F-MSG and LR method.These numerical results are given in Table 5.
It can be seen from Table 5 that LR has the duality gap value and F-MSG performs much better than LR, since duality gap value is reduced to zero.The quality of the solution is improved when the duality gap is decreased.In this paper zero duality gap is achieved and feasible solution is attained by applying the F-MSG method.

Stage
Unit Combination

Conclusion
The F-MSG method is firstly applied to solve the thermal UC problem which is an important and hard-solving problem in power system engineering.Feasible cost, dual cost, and duality gap values from F-MSG and LR methods are compared.The most attractive feature of the proposed approach is that the duality gap value of F-MSG method over the scheduled time horizon is zero.Also, the total production cost over the scheduled time horizon of the F-MSG method is less than the ones of the standard LR methods.The results of the F-MSG method for solving the thermal UC problem are very impressive and the quality of feasible solution is significantly improved.This approach can be applied to the thermal UC problem for any size of systems to obtain the feasible solutions.

Table 1 :
Unit characteristics for four-unit Tuncbilek thermal plant.Unit No P min MW P max MW SC $ a $/h b $/MWh c $/MW 2 h RU $ RU $ T up h T down h

Table 3 :
UC schedule for F-MSG method.

Table 4 :
UC schedule for LR method.

Table 5 :
The numerical results for F-MSG and LR methods.
Maximum available capacity of unit i at hour t MW P t min,i : Minimum available capacity of unit i at hour t MW P t r : System spinning reserve at hour t MW RU i : Ramp-up rate of unit i MW/h RD i : Ramp-down rate of unit i MW/h Continuously off-time of unit i up to hour t, h T t i,on : Continuously on-time of unit i up to hour t, h T i,up : Minimum up-time of unit i, h