Unit Commitment Solution Methods Explained
Unit Commitment Solution Methods Explained
International Scholarly and Scientific Research & Innovation 1(11) 2007 1638 [Link]/1307-6892/11207
World Academy of Science, Engineering and Technology
International Journal of Energy and Power Engineering
Vol:1, No:11, 2007
single unit dynamic programming and inclusion of ramp rate such as simulated annealing [17], tabu search [18], genetic
constraints are proposed by Guan et al. [5]. Kuloor et al. [6] algorithms [19]-[21], and greedy randomized adaptive search
show environmental constraint incorporation technique in unit procedure [22].
commitment problem formulation.
Due to the imperfections of the dynamic programming III. MIXED INTEGER PROGRAMMING
algorithm, the application of a unit commitment expert system Integer programming optimizes integer function of integer
[7], [8] has been proposed to supplement this method. The variables. A modification of standard integer programming
constraints, which are difficult or impractical to be that allows non-integer function is known as mixed-integer
implemented in unit commitment algorithm, can be handled programming (MIP). MIP treates the objective and constraint
by this expert system. In this paper, various proposed methods functions as continuous and the variables as integers.
for unit commitment have been described. Finally a Branch and bound is one of the techniques used for the
comprehensive algorithm [9], [10] for solving industry-grade solution of the integer problem. It is a technique to solve a
unit commitment problem is described. discrete variable problem by solving a sequence of simpler
problems derived from the original problem. In solving using
II. HEURISTIC METHODS the branch and bound method, one needs to define 1) the
problems in the branch and bound three, 2) the method of
Heuristic methods are non-rigorous computer aided
solving each problem on the tree, and 3) the method of
empirical methods, which make the unit commitment
searching the tree. The branch and bound tree is fully
International Science Index, Energy and Power Engineering Vol:1, No:11, 2007 [Link]/Publication/11207
decisions according to a pre-calculated priority list and described if one defines the problem corresponding to the top
incorporate all the operating constraints heuristically. node of the tree and the method of obtaining the children of
Baldwin et al. [11] have used a heuristic approach for unit any node [1].
commitment. All units are shut down and started up in strict Dillon et al. [23] have formulated the unit commitment
priority order. The priority list is prepared on the basis of the problem as a linear MIP problem. Then they have used
average full load cost of each unit, where the average full load standard integer programming algorithm for solving the
cost is the net heat rate at full load times the fuel cost. commitment schedule.
Kerr et al. [12] have proposed heuristic approach that One of the proposed MIP methods [24] transforms the
begins with an initial feasible schedule and then follows the linear optimization problems that arise during the search
sequence of steps for adjusting the starting and stopping times procedure in the branch and bound algorithm into capacitated
accordingly to reduce the cost of operation. transhipment problems. These are then solved efficiently by a
The heuristic method proposed by Happ et al. [13] uses a network-based solution procedure.
sub-optimizer to get a feasible and near optimal commitment. Bond and Fox [25] presented an algorithm based on a
Then that method employs an optimizer to further optimize combination of mixed integer-linear and dynamic
the schedule. The optimizer uses a sequence of steps programming. Mixed integer-linear programming is used to
repeatedly until no further reduction in cost is observed. determine feasible combinations of units at each scheduling
point, while a novel dynamic programming approach
Heuristic approach to the short-term unit commitment
identifies promising scheduling routes in the time domain.
problem has also been implemented using an expert system
MIP models proposed so far have employed linear cost
[14], [15].
functions although more accurate cost models are available.
Heuristic methods have the following advantages [16]:
To date, branch and bound techniques have only been
• are flexible and allow for the consideration of employed on moderate sized systems using linear models [1].
practical operating constraints On the contrary, many available economic dispatch algorithms
use quadratic cost curves. Moreover, the present trend is
• feasible solutions if there are any are usually
toward improved modeling of unit input/output characteristics
obtained
with more detailed non-linear models. The MIP based
• computational requirements in terms of memory approach, using only linear models, was found to take long
and running time are modest computation time. The MIP formulation of unit commitment
would become a very large problem demanding extremely
The main shortcoming of heuristic methods is that they
long computation time if applied to a typical generation mix
cannot guarantee the optimal solutions or even furnish an
with more detailed non-linear models.
estimate of the magnitude of their sub-optimality. This aspect
becomes rather significant in large-scale power systems, as a
IV. BENDERS DECOMPOSITION
small percentage, e.g. 0.5%, in the costs of unit commitment
schedules represents a substantial financial annual saving. In Benders Decomposition method [26]-[29], the unit
Therefore, it is advantageous to employ more rigorous commitment problem is decomposed into a master problem
methods compared with heuristics methods to generate more involving only the discrete commitment variables and a
economical solutions as the size of a system grows, despite the subproblem involving the continuous generation variables.
requirement of comparatively large computational efforts. The subproblem corresponds to the economic dispatch
problem with a given commitment. The marginal costs, for
More recently, metaheuristic approaches have been used
International Scholarly and Scientific Research & Innovation 1(11) 2007 1639 [Link]/1307-6892/11207
World Academy of Science, Engineering and Technology
International Journal of Energy and Power Engineering
Vol:1, No:11, 2007
each hour, from the subproblem are used to constrain the units, then the state space size is 31 which is a reasonable
allowed commitments in the master problem. The master number.
problem supplies commitments to the economic dispatch The approaches [31], [36] used selection techniques for
problem. The master problem and the economic dispatch choosing the most promising states from all possible states
subproblem are solved iteratively until the solution converges. and implemented approximate economic dispatch subroutines
Any economic dispatch routine can be applied for solving the to reduce computer running time requirement. Variable
subproblem. window truncated dynamic programming that adjusts the
The major difficulty in Benders decomposition approach is window size according to the incremental load demands in
the determination of the solution of the master problem, which adjacent hours and controls the program execution time have
is still regarded as a large scale integer optimization problem. been proposed. Kumar and Palanisamy [37] proposed a two-
Turgeon [26] has solved the master problem by a variational step process that uses a direct computation Hopfield neural
approach and a branch and bound algorithm whereas network to generate economic dispatch. Then using dynamic
Baptistella and Geromel [27] have solved it by relaxation in programming the generator schedule is produced.
the master level of Benders decomposition approach. In order The approaches [7], [8] have an integrated expert system
to improve the efficiency, some of the constraints which are into the truncated dynamic programming based unit
difficult to handle, such as nonlinear minimum up and down commitment program to check and modify commitment
time constraints, are replaced by simpler constraints in the results.
actual formulation of the scheduling algorithm. For instance,
International Science Index, Energy and Power Engineering Vol:1, No:11, 2007 [Link]/Publication/11207
Habibollahzadeh and Bubenko [28] did not use the minimum VI. LAGRANGIAN RELAXATION
up and down time constraints in their mathematical model but, In Lagrangian relaxation approaches, unit commitment
instead, included a constraint that allowed only one problem is formulated in terms of 1) a cost function, that is the
commitment per day for each unit. Ma and Shahidehpour [29] sum of terms each involving a single unit, 2) a set of
dealt with transmission-constraint by introducing a proper constraints involving a single unit, and 3) a set of coupling
constraint called Benders cut. constraints (the generation and reserve requirements), one for
each hour in the study period, involving all units. Cohen and
V. DYNAMIC PROGRAMMING Sherkat [1] have reported that an approximate solution to this
problem can be obtained by adjoining the coupling constraints
Dynamic programming acts as an important optimization onto the cost using Lagrangian multipliers. The cost function
technique with broad application in many areas [30]. Dynamic (primal objective function) of the unit commitment problem is
programming decomposes a problem into a series of smaller relaxed to the power balance and the generating constraints
problems, solves the small problems, and develops an optimal via two sets of Lagrangian multipliers to form a Lagrangian
solution to the original problem step-by-step. The optimal dual function. The dual problem is then decoupled into small
solution is developed from the subproblem recursively. subproblems which are solved separately with the remaining
In its fundamental form, the dynamic programming constraints. Meanwhile, the dual function is maximized with
algorithm for unit commitment problem examines every respect to the Lagrangian multipliers, usually by a series of
possible state in every interval. Some of these states are iterations.
rejected instantly because they are found infeasible. But even, Unfortunately, duality theory has shown that for nonconvex
for an average size utility, a large number of feasible states problems there will be a duality gap between the cost obtained
will exist and the requirement of execution time will stretch by solving the relaxed problem and the optimal cost of the
the capability of even the largest computers. Hence many original problem. Since the commitment decision variables are
proposed techniques used some sort of simplification and discrete, the unit commitment problem is nonconvex. Due to
approximation to the fundamental dynamic programming the nonconvexity of the unit commitment problem, the dual
algorithm [31]. solution seldom satisfies the power balance constraints and the
The approach, first used by Lowery [32], and later refined reserve constraints. Hence, in addition to solving the dual
by Ayoub and Patton [33], selected unit generation output as a problem, a suboptimal feasible solution is usually searched
state variable and on-line capacity as the stages. Ayoub and near the dual optimal point. Even though the dual optimal
Patton included probabilistic techniques for reserve solutions may violate the feasibilities of the original problem,
determination in the developed code. they may usually provide sharp lower bounds.
A typical approach [34], [35] determines some nominal Muckstadt and Koenig [38] employed Lagrangian
commitment, which is determined to be good for each hour. relaxation to replace the common linear programming
Choices that have to be considered are minimum number of relaxation approach, which dropped the integrality
units, in the priority ordering, needed to meet the reserve requirements of the variables, in the fathoming process of
constraints and the result of the above priority list branch and bound algorithms. This causes a significant
optimization. A set of units in the priority list about the improvement of computational efficiency compared with
nominal commitment are then chosen for optimization - the previous branch and bound algorithms.
units below that set are assumed to be committed while the Based on the sharp bound provided by the Lagrangian dual
units above that set are assumed to be off. If this set contains 5 optimum, it is expected that a suboptimal feasible solution
near the dual optimal point can be accepted as a proper
International Scholarly and Scientific Research & Innovation 1(11) 2007 1640 [Link]/1307-6892/11207
World Academy of Science, Engineering and Technology
International Journal of Energy and Power Engineering
Vol:1, No:11, 2007
solution for the primal problem. Merlin and Sandrin [39] generation levels of two consecutive hours is difficult to deal
presented a more direct and fairly efficient technique based on with. A straightforward application of dynamic programming
this idea. In this algorithm, a modified subgradient method technique may lead to suboptimal results as illustrated by
was used which incorporated the search of the suboptimal Guan et al. [5].
feasible solution along the direction for maximizing the dual Guan et al. [5] have presented an approach that does not
function. discretize generation levels and handles ramp rate constraints
Some methods have been suggested which can be systematically. The thermal subproblem without ramp rate
implemented in systems with fuel constraints. An additional constraint is solved by first constructing a state transition
set of multipliers has been introduced to associate the fuel diagram where the optimal generation levels of all up states
constraints with the primal function. Cohen and Wan [40] are computed without discretizing generation levels. Dynamic
proposed a successive approximation approach for altering the programming technique is then applied with only a few well-
three sets of multipliers. Three iteration loops were structured states. This eliminates the difficult trade-off
constructed to update these multipliers independently, between computational requirements and the accuracy as
according to their corresponding constraints. The iteration of needed by most approaches that discretize generation levels.
these three sets of multipliers was done by enclosing the three Ramp rate constraints are relaxed by introducing an additional
iteration loops with an external loop. Once the dual optimal set of multipliers for a unit with the constraints. The
had been found, the commitment of the fuel-constrained units subproblem is then solved as if there were no ramp rate
would be fixed. Then, any violations of the constraints were constraint. An intermediate level is introduced to update this
International Science Index, Energy and Power Engineering Vol:1, No:11, 2007 [Link]/Publication/11207
International Scholarly and Scientific Research & Innovation 1(11) 2007 1641 [Link]/1307-6892/11207
World Academy of Science, Engineering and Technology
International Journal of Energy and Power Engineering
Vol:1, No:11, 2007
the least rigorous of the mathematical programming as cycling of gas turbine and steam turbine units, group
techniques but offers the best performance in computing constraints, etc.
requirements. Under the current computing technology, it
seems that Lagrangian relaxation is the viable mathematical VIII. CONCLUSION
programming technique to solve large-scale unit commitment.
The savings potential of unit commitment in absolute terms
This is reflected by the fact that they are widely reported in
increases with the increase of power system size. The effort to
the literature.
develop a unit commitment approach that is able to handle
To get an industry-grade efficient algorithm based on
Lagrangian relaxation, the techniques like non-discretization large systems consisting of both thermal and hydro generating
of generation levels, handling of ramp rate and environmental units therefore offers a large profitable return. In order to be
constraints [5], [6] can be incorporated in the algorithm [2], feasible, the method to be developed must be flexible,
[4]. Transmission loss can be incorporated using a general efficient and reliable. Many proposed methods have been
transmission loss formula [42] whose expression has a similar described along with their strengths and weaknesses. It seems
quadratic form to the B matrix loss formulation. The that the comprehensive algorithm that combines the strengths
generation cost and water discharge rate can be represented as of different methods and overcomes each other’s weaknesses
continuously increasing quadratic function of the generation. would be a suitable approach for solving industry-grade unit
Finally, an efficient unit commitment expert system [8] can be commitment problem.
developed as a supplement to the Lagrangian relaxation
International Science Index, Energy and Power Engineering Vol:1, No:11, 2007 [Link]/Publication/11207
method. REFERENCES
To fulfill the aim cited above, hydrothermal scheduling [1] A. I. Cohen and V. R. Sherkat, “Optimization based methods for
based Lagrangian relaxation (HTSLR) approach to solve the operations scheduling”, Proceedings of the IEEE, vol. 75, no. 12, pp.
unit commitment problem for a large practical system 1574-1591, 1987.
comprising both thermal and hydro generating units is [2] S. K. Tong and S. M. Shahidehpour, “An innovative approach to
generation scheduling in large-scale hydro-thermal power systems with
proposed [9]. Commitment states of thermal units are obtained fuel constrained units”, IEEE Trans. on Power Systems, vol. 5, no. 2, pp.
by solving thermal subproblems only. To achieve the output 665-673, 1990.
levels of hydro units, the hydrothermal scheduling is [3] N. P. Padhy, “Unit commitment - a bibliographical survey”, IEEE Trans.
performed with a thermal unit commitment schedule obtained on Power Systems, vol. 19, no. 2, pp. 1196-1205, 2004.
[4] A. Aoki, T. Satoh, M. Itoh, T. Ichimori, and K. Masegi, “Unit
by solving only thermal subproblems. Extensive constraints commitment in a large scale power system including fuel constrained
are taken into account such as status restriction of individual thermal and pumped storage hydro”, IEEE Trans. on Power Systems,
generating units i.e. must run, must out, base load, cycling and vol. 2, no. 4, pp. 1077-1084, 1987.
peaking, power balance, spinning reserve, minimum up/down [5] X. Guan, P. B. Luh, H. Yan, and J. A. Amalfi, “An optimization-based
method for unit commitment” International Journal of Electrical Power
time, capacity limits, ramp rate, limited generation for the first and Energy Systems, vol. 14, no. 1, pp. 9-17, 1992.
and last hour, sulfur oxide emission and hydro constraints. [6] S. Kuloor, G. S. Hope, and O. P. Malik, “Environmentally constrained
Non-linear functions are employed for thermal generation cost unit commitment”, IEE Proceedings- Generation, Transmission &
and water discharge rate. A general transmission loss formula Distribution, vol. 139, no. 2, pp. 122-128, 1992.
[7] S. Mokhtari, J. Singh, and B. Wollenberg, “A unit commitment expert
[42] whose expression has a similar quadratic form to the B system”, IEEE Trans. on Power Systems, vol. 3, no. 1, pp. 272-277,
matrix loss formulation has been utilized for incorporating 1988.
transmission loss. [8] M. S. Salam, A. R. Hamdan, and K. M. Nor, “Integrating an expert
In the HTSLR approach, the commitment schedule may be system into a thermal unit commitment algorithm”, IEE Proceedings-
Generation, Transmission & Distribution, vol. 138, no. 6, pp. 553-559,
so sensitive to the variations of the Lagrange multipliers that a 1991.
slight modification of the multipliers may change the status of [9] M. S. Salam, K. M. Nor, and A. R. Hamdan, “Hydrothermal scheduling
several units. This sensitivity problem is more serious for based Lagrangian relaxation approach to hydrothermal coordination”,
systems having several groups of identical units. Even though IEEE Trans. on Power Systems, vol. 13, no. 1, pp. 226-235, 1998.
[10] M. S. Salam, K. M. Nor, and A. R. Hamdan, “Comprehensive algorithm
fuel costs of identical units can be slightly modified to make
for hydrothermal co-ordination”, IEE Proceedings- Generation,
small differences among cost characteristics, this sensitivity Transmission & Distribution, vol. 144, no. 5, pp. 482-488, Sept. 1997.
problem still exists. In other words, unnecessary commitment [11] C. J. Baldwin, K. M. Dale, and R. F. Dittrich, “A study of the economic
of some units may be possible in the solution given by this shutdown of generating units in daily dispatch”, AIEE Trans. Part III,
Power Apparatus and Systems, vol. 78, pp. 1272-1284, 1959.
method. In order to overcome this difficulty, a refinement
[12] R. H. Kerr, J. L. Scheidt, A. J. Fontana, and J. K. Wiley, “Unit
process similar to that proposed by Tong and Shahidehpour commitment”, IEEE Trans. on Power Apparatus and Systems, vol. 85,
[2] has been developed for HTSLR approach. This refinement no. 5, pp. 417-421, 1966.
process inspects some candidate units whose shutdown may [13] H. H. Happ, R. C. Johnson, and W. J. Wright, “Large scale hydro-
result in additional reduction of the operating cost. thermal unit commitment: method and results”, IEEE Trans. on Power
Apparatus and Systems, vol. 90, no. 3, pp. 1373-1384, 1971.
A unit commitment expert system [8] has also been [14] S. K. Tong, S. M. Shahidehpour, and Z. Ouyang, “A heuristic short-term
developed [10]. It was employed as a preprocessor as well as a unit commitment”, IEEE Trans. on Power Systems, vol. 6, no. 3, pp.
postprocessor to the HTSLR based unit commitment program 1210-1216, 1991.
[15] S. Li, S. M. Shahidehpour, and C. Wang, “Promoting the application of
to check and alter commitment results by adjusting the input expert systems in short-term unit commitment”, IEEE Trans. on Power
data if necessary. It handles constraints which are difficult or Systems, vol. 8, no. 1, pp. 286-292, 1993.
impractical to be implemented in commitment algorithm such [16] E. Khodaverdian, A. Brameller, and R. M. Dunnett, “Semi-rigorous
thermal unit commitment for large scale electrical power systems”, IEE
International Scholarly and Scientific Research & Innovation 1(11) 2007 1642 [Link]/1307-6892/11207
World Academy of Science, Engineering and Technology
International Journal of Energy and Power Engineering
Vol:1, No:11, 2007
Proceedings- Generation, Transmission & Distribution, vol. 133, no. 4, [40] A. I. Cohen and S. H. Wan, “A method for solving the fuel constrained
pp. 157-164, 1986. unit commitment problem” IEEE Trans. on Power Systems, vol. 2, no. 3,
[17] F. Zhuang and F. D. Galiana, “Unit Commitment by Simulated pp. 608-614, 1987.
Annealing,” IEEE Trans. on Power Systems, vol. 5, no. 1, pp. 311-317, [41] H. Yan, P. B. Luh, X. Guan, and P. M. Rogan, “Scheduling of
1990. hydrothermal power systems”, IEEE Trans. on Power Systems, vol. 8,
[18] B. Xiaomin, S. M. Shahidehpour, and Y. Erkeng, “Constrained unit no. 3, pp. 1358-1365, 1993.
commitment by using tabu search algorithm,” in Proc. Int Conf. on [42] O.I. Elgerd, Electric energy systems theory: an introduction. McGraw-
Electrical Engineering, vol. 2, 1996, pp.1088-1092. Hill Book Company, New York, 1971, pp. 294-296.
[19] S. A. Kazarlis, A. G. Bakirtzis, and V. Petridis, “A genetic algorithm
solution to the unit commitment problem,” IEEE Trans. on Power
Systems, vol. 11, no. 1, pp. 83-90, 1996.
[20] S. O. Orero and M. R. Irving, “A genetic algorithm modeling framework
and solution technique for short term optimal hydrothermal scheduling,”
IEEE Trans. on Power Systems, vol. 13, no. 2, pp. 501-516, 1998.
[21] S. Liyong, Z. Yan, and J. Chuanwen, “A matrix real-coded genetic
algorithm to the unit commitment problem”, Electric Power Systems
Research, vol. 76, no. 9-10, pp. 716-728, 2006.
[22] A. Viana, J. P. de Sousa, and M. Matos, “A new metaheuristic approach
to the unit commitment problem”, 14th Power Systems Computation
Conference, Sevilla, Spain, 24-28 June 2002, Session 05, Paper 5.
[23] T. S. Dillon, K. W. Edwin, H. D. Kochs, and R. J. Taud, “Integer
programming approach to the problem of optimal unit commitment with
probabilistic reserve determination”, IEEE Trans. on Power Apparatus
International Science Index, Energy and Power Engineering Vol:1, No:11, 2007 [Link]/Publication/11207
International Scholarly and Scientific Research & Innovation 1(11) 2007 1643 [Link]/1307-6892/11207