PSOC Unit 4
PSOC Unit 4
com
UNIT –IV
www.Vidyarthiplus.com
www.Vidyarthiplus.com
PART-A
PART-B
www.Vidyarthiplus.com
www.Vidyarthiplus.com
UNIT – IV
1.Mention the constraints in unit commitment problem.
1.Spinning reserve:
The change in generation required to meet power demand is called as participation factor.
Spinning reserve is the total amount of generation available from all units synchronized
on the system minus the present load and losses being supplied.
www.Vidyarthiplus.com
www.Vidyarthiplus.com
7.Write the coordination equation for (i) With loss (ii) Without loss.
Refer PART-B
8.What is meant by incremental fuel cost curve?
Refer PART - B
9.What are the factors affecting the cost of generation?
(i) Equality constraint
(ii) Inequality constraints
(a) Generator Constraints
(b) Voltage Constraints
(c) Running space capacity constraints
(d) Transformer tap settings.
(e) Transmission line constraints
www.Vidyarthiplus.com
www.Vidyarthiplus.com
ECONOMIC DISPATCH
Economic Operation of Power Systems
One of the earliest applications of on-line centralized control was to provide a central
facility, to operate economically, several generating plants supplying the loads of the system.
Modern integrated systems have different types of generating plants, such as coal fired
thermal plants, hydel plants, nuclear plants, oil and natural gas units etc. The capital
investment, operation and maintenance costs are different for different types of plants.
During operation of the plant, a generator may be in one of the following states:
Performance Curves
Input-Output Curve
This is the fundamental curve for a thermal plant and is a plot of the input in British
thermal units (Btu) per hour versus the power output of the plant in MW as shown in Fig 2.1
www.Vidyarthiplus.com
www.Vidyarthiplus.com
The incremental cost is the product of incremental fuel rate and fuel cost (Rs / Btu or $ /Btu).
The curve in shown in Fig. 4. The unit of the incremental fuel cost is Rs / MWh or $ /MWh.
In general, the fuel cost Fi for a plant, is approximated as a quadratic function of the
generated output PGi.
2
Fi = ai + bi PGi + ci PG Rs / h --------------------------------- (4)
The incremental fuel cost is a measure of how costly it will be produce an increment of
power. The incremental production cost, is made up of incremental fuel cost plus the
incremental cost of labour, water, maintenance etc. which can be taken to be some percentage
of the incremental fuel cost, instead of resorting to a rigorous mathematical model. The cost
www.Vidyarthiplus.com
www.Vidyarthiplus.com
curve can be approximated by a linear curve. While there is negligible operating cost for a
hydel plant, there is a limitation on the power output possible. In any plant, all units normally
operate between PGmin, the minimum loading limit, below which it is technically infeasible
to operate a unit and PGmax, which is the maximum output limit.
Solution Methods:
1. Lagrange Multiplier method
2. Lamda iteration method
3. Gradient method
4. Dynamic programming
5. Evolutionary Computation techniques
--------------------------------------------------------------------------------------
1.. The Economic Dispatch Problem without Loss.
------------------------------------------ (6)
This is a constrained optimization problem that may be attacked formally using advanced
www.Vidyarthiplus.com
www.Vidyarthiplus.com
calculus methods that involve the Lagrange function. In order to establish the necessary
conditions for an extreme value of the objective function, add the constraint function to the
objective function after the constraint function has been multiplied by an undetermined
multiplier. This is known as the Lagrange function and is shown in Eq(7)
------------------------------------------------ (7)
The necessary conditions for an extreme value of the objective function result when we
take the first derivative of the Lagrange function with respect to each of the independent
variables and set the derivatives equal to zero. In this case,there are N+1 variables, the N
values of power output, Pi, plus the undetermined Lagrange multiplier, λ. The derivative of
the Lagrange function with respect to the undetermined multiplier merely gives back the
constraint equation. On the other hand, the N equations that result when we take the
partial derivative of the Lagrange function with respect to the power output values one at a
time give the set of equations shown as Eq. 8.
When we recognize the inequality constraints, then the necessary conditions may be
expanded slightly as shown in the set of equations making up Eq. 9
--------------------------------------------- (9)
EXAMPLE 1
Suppose that we wish to determine the economic operating point for these three units
when delivering a total of 850 MW. Before this problem can be solved,
the fuel cost of each unit must be specified. Let the following fuel costs are in effect.
Unit 1: Coal-fired steam unit: Max output = 600 MW Min output = 150 MW
Input-output curve:
www.Vidyarthiplus.com
www.Vidyarthiplus.com
Input-output curve:
Unit 3: Oil-fired steam unit: Max output = 200 MW, Min output = 50 MW
Input-output curve:
Then
www.Vidyarthiplus.com
www.Vidyarthiplus.com
EXAMPLE 2
Suppose the price of coal decreased to 0.9 P/MBtu. The fuel cost function for unit 1 becomes
If one goes about the solution exactly as done here, the results are
This solution meets the constraint requiring total generation to equal 850 MW, but units 1 and
3 are not within limit. To solve for the most economic dispatch while meeting unit limits,
suppose unit 1 is set to its maximum output and unit 3 to its minimum output. The dispatch
becomes
Next, calculate the incremental cost for units 1 and 3 to see if they meet the conditions.
Note that the incremental cost for unit 1 is less than λ, so unit 1 should be at its maximum.
However, the incremental cost for unit 3 is not greater than ,λ so unit 3 should not be forced
to its minimum. Thus, to find the optimal dispatch, allow the incremental cost at units 2 and 3
to equal λ as follows.
www.Vidyarthiplus.com
www.Vidyarthiplus.com
------------------ (10)
The same procedure is followed in the formal sense to establish the necessary
conditions for a minimum-cost operating solution, The Lagrange function is shown in Eq.11.
In taking the derivative of the Lagrange function with respect to each of the individual power
outputs, Pi, it must be recognized that the loss in the transmission network, P loss is a function
of the network impedances and the currents flowing in the network. For our purposes, the
currents will be considered only as a function of the independent variables P i and the load
Pload taking the derivative of the Lagrange function with respect to any one of the N values of
Pi results in Eq. 11. collectively as the coordination equations.
----------------------------------- (11)
It is much more difficult to solve this set of equations than the previous set with no
losses since this second set involves the computation of the network loss in order to establish
the validity of the solution in satisfying the constraint equation. There have been two general
approaches to the solution of this problem. The first is the development of a mathematical
expression for the losses in the network solely as a function of the power output of each of
the units. This is the loss-formula method discussed at some length in Kirchmayer’s
Economic Operation of Power Systems. The other basic approach to the solution of this
problem is to incorporate the power flow equations as essential constraints in the formal
establishment of the optimization problem. This general approach is known as the optimal
power flow.
www.Vidyarthiplus.com
www.Vidyarthiplus.com
www.Vidyarthiplus.com
www.Vidyarthiplus.com
Lambda-iteration method
We could, for example, store tables of data within the computer and interpolate
between the stored power points to find exact power output for a specified value of
incremental cost rate. Another approach would be to develop an analytical function for the
power output as a function of the incremental cost rate, store this function (or its coefficients)
in the computer, and use this to establish the output of each of the individual units.
www.Vidyarthiplus.com
www.Vidyarthiplus.com
Example:3
Given the generator cost functions found in Example 2.1, solve for the economic
0
dispatch of generation with a total load of 800 MW.Using α = 100 and starting from P1 = 300
0 0
MW, P2 = 200 MW, and P3 =300 MW, we set the initial value of λ. equal to the average of
the incremental costs of the generators at their starting generation values. This value is
9.4484.The progress of the gradient search is shown in Table 3.2. The table shows that the
iterations have led to no solution at all. Attempts to use this formulation
will result in difficulty as the gradient cannot guarantee that the adjustment to the generators
will result in a schedule that meets the correct total load of 800 MW.A simple variation of
this technique is to realize that one of the generators is always a dependent variable and
remove it from the problem. In this case, we pick P3 and use the following:
Note that this function stands by itself as a function of two variables with no
load-generation balance constraint (and no λ). The cost can be minimized by
a gradient method and in this case the gradient is:
\
Note that this gradient goes to the zero vector when the incremental cost at generator 3 is
equal to that at generators 1 and 2. The gradient steps are performed in the same manner as
previously, where:
www.Vidyarthiplus.com
www.Vidyarthiplus.com
------------------------------------ (12)
Each time a gradient step is made, the generation at generator 3 is set to 800 minus the sum
of the generation at generators 1 and 2. This method is often called the “reduced gradient”
because of the smaller number of variables.
----------------------------------------- (13)
This is true for each of the N units on the system, so that
---------------------------------------(14)
The total change in generation (=change in total system demand) is, of course, the sum of the
individual unit changes. Let Pd be the total demand on the generators (where Pload+Ploss&),
then
www.Vidyarthiplus.com
www.Vidyarthiplus.com
---------------------------------------- (15)
The earlier equation, 15, can be used to find the participation factor for each unit as follows
--------------------------------------- (16)
The computer implementation of such a scheme of economic dispatch is straightforward. It
might be done by provision of tables of the values of FY as a function of the load levels and
devising a simple scheme to take the existing load plus the projected increase to look up these
data and compute the factors. somewhat less elegant scheme to provide participation factors
would involve a repeat economic dispatch calculation at. The base-point economic generation
values are then subtracted from the new economic generation values and the difference
divided to provide the participation factors. This scheme works well in computer
implementations where the execution time for the economic dispatch is short and will always
give consistent answers when units reach limits, pass through break points on piecewise
linear incremental cost functions, or have nonconvex cost curves.
EXAMPLE 2.4
Starting from the optimal economic solution found in Example 2A; use the participation
factor method to calculate the dispatch for a total load of 900 MW.
www.Vidyarthiplus.com
www.Vidyarthiplus.com
UNIT COMMITMENT
The life style of a modern man follows regular habits and hence the present society
also follows regularly repeated cycles or pattern in daily life. Therefore, the consumption of
electrical energy also follows a predictable daily, weekly and seasonal pattern. There are
periods of high power consumption as well as low power consumption. It is therefore
possible to commit the generating units from the available capacity into service to meet the
demand. The previous discussions all deal with the computational aspects for allocating load
to a plant in the most economical manner. For a given combination of plants the
determination of optimal combination of plants for operation at any one time is also desired
for carrying out the aforesaid task. The plant commitment and unit ordering schedules extend
the period of optimization from a few minutes to several hours. From daily schedules weekly
patterns can be developed. Likewise, monthly, seasonal and annual schedules can be prepared
taking into consideration the repetitive nature of the load demand and seasonal variations.
Unit commitment schedules are thus required for economically committing the units in plants
to service with the time at which individual units should be taken out from or returned to
service.
www.Vidyarthiplus.com
www.Vidyarthiplus.com
from all units synchronized (i.e., spinning) on the system, minus the present load and losses
being supplied. Spinning reserve must be carried so that the loss of one or more units does not
cause too far a drop in system frequency. Quite simply, if one unit is lost, there must be
ample reserve on the other units to make up for the loss in a specified time period.Spinning
reserve must be allocated to obey certain rules, usually set by regional reliability councils (in
the United States) that specify how the reserve is to be allocated to various units. Typical
rules specify that reserve must be a given percentage of forecasted peak demand, or that
reserve must be capable of making up the loss of the most heavily loaded unit in a given
period of time.Others calculate reserve requirements as a function of the probability of not
having sufficient generation to meet the load.Not only must the reserve be sufficient to make
up for a generation-unit failure, but the reserves must be allocated among fast-responding
units and slow-responding units. This allows the automatic generation control system to
restore frequency and interchange quickly in the event of a generating-unit outage. Beyond
spinning reserve, the unit commitment problem may involve various classes of “scheduled
reserves” or “off-line” reserves. These include quick-start
diesel or gas-turbine units as well as most hydro-units and pumped-storage hydro-units that
can be brought on-line, synchronized, and brought up to full capacity quickly. As such, these
units can be “counted” in the overall reserve assessment, as long as their time to come up to
full capacity is taken into account. Reserves, finally, must be spread around the power system
to avoid transmission system limitations (often called “bottling” of reserves) and to allow
various parts of the system to run as “islands,” should they become electrically disconnected.
www.Vidyarthiplus.com
www.Vidyarthiplus.com
operating temperature. There are two approaches to treating a thermal unit during its down
period. The first allows the unit’s boiler to cool down and then heat back up to operating
temperature in time for a scheduled turn on. The second (called banking) requires that
sufficient energy be input to the boiler to just maintain operating temperature. The costs for
the two can be compared so that, if possible, the best approach (cooling or banking) can be
chosen.
Cc = cold-start cost (MBtu)
F = fuel cost
Cf= fixed cost (includes crew expense, maintenance expenses) (in R)
α = thermal time constant for the unit
t = time (h) the unit was cooled
Start-up cost when banking = Ct x t x F+Cf
Where
Up to a certain number of hours, the cost of banking will be less than the cost of cooling, as is
illustrated in Figure 5.3.Finally, the capacity limits of thermal units may change frequently,
due to maintenance or unscheduled outages of various equipment in the plant; this must also
be taken
Other Constraints
Hydro-Constraints
Unit commitment cannot be completely separated from the scheduling of hydro-units. In this
text, we will assume that the hydrothermal scheduling (or “coordination”) problem can be
separated from the unit commitment problem. We, of course, cannot assert flatly that our
treatment in this fashion will always
result in an optimal solution.
Hydro-Constraints
www.Vidyarthiplus.com
www.Vidyarthiplus.com
Must Run
Some units are given a must-run status during certain times of the year for reason of voltage
support on the transmission network or for such purposes as supply of steam for uses outside
the steam plant itself.
Fuel Constraints
We will treat the “fuel scheduling” problem system in which some units have limited fuel, or
else have constraints that require them to burn a specified amount of fuel in a given time,
presents a most challenging unit commitment problem.
-----------------------------------------(19)
For the total period of M intervals, the maximum number of possible combinations is (2N -
l)M, which can become a horrid number to think about.
For example, take a 24-h period (e.g., 24 one-hour intervals) and consider systems with 5, 10,
24
20, and 40 units. The value of (2N - 1) becomes the following.
www.Vidyarthiplus.com
www.Vidyarthiplus.com
These very large numbers are the upper bounds for the number of enumerations required.
Fortunately, the constraints on the units and the load-capacity relationships of typical utility
systems are such that we do not approach these large numbers. Nevertheless, the real
practical barrier in the optimized unit commitment problem is the high dimensionality of the
possible solution space. The most talked-about techniques for the solution of the unit
commitment problem are:
1. Priority-list schemes,
2. Dynamic programming (DP),
3. Lagrange relation (LR).
www.Vidyarthiplus.com
www.Vidyarthiplus.com
5. Dynamic-Programming Solution
Dynamic programming has many advantages over the enumeration scheme, the chief
advantage being a reduction in the dimensionality of the problem. Suppose we have found
units in a system and any combination of them could serve the (single) load. There would be
a maximum of 24 - 1 = 15 combinations to test. However, if a strict priority order is imposed,
there are only four combinations to try:
Priority 1 unit
Priority 1 unit + Priority 2 unit
Priority 1 unit + Priority 2 unit + Priority 3 unit
Priority 1 unit + Priority 2 unit + Priority 3 unit + Priority 4 unit
The imposition of a priority list arranged in order of the full-load averagecost rate would
result in a theoretically correct dispatch and commitment only if:
www.Vidyarthiplus.com
www.Vidyarthiplus.com
Forward DP Approach
One could set up a dynamic-programming algorithm to run backward in time starting
from the final hour to be studied, back to the initial hour. Conversely, one could set up the
algorithm to run forward in time from the initial hour to the final hour. The forward approach
has distinct advantages in solving generator unit commitment. For example, if the start-up
cost of a unit is a function of the time it has been off-line (i.e., its temperature), then a
forward dynamic-program approach is more suitable since the previous history of the
unit can be computed at each stage. There are other practical reasons for going forward. The
initial conditions are easily specified and the computations can go forward in time as long as
required. A forward dynamic-programming algorithm is shown by the flowchart
The recursive algorithm to compute the minimum cost in hour K with combinati
Fcost(K,I)= min[Pcost(K,I)+Scost(K-1,L:K,I)+Fcost(K-1,L)] ---------------------------------- (20)
Where
Fcost(K, I ) = least total cost to arrive at state ( K , I )
Pcost(KI, ) = production cost for state ( K ,I )
Scost(K - 1, L: K , I)= transition cost from state (K - 1, L) to state ( K , I )
State (K, 1) is the Zth combination in hour K. For the forward dynamic programming
approach, we define a strategy as the transition, or path, from one state at a given hour to a
state at the next hour.
Note that two new variables, X and N, have been introduced in Figure 2.11
X = number of states to search each period
N = number of strategies, or paths, to save at each step
These variables allow control of the computational effort (see below Figure).For complete
n
enumeration, the maximum number of the value of X or N is 2 – 1
www.Vidyarthiplus.com
www.Vidyarthiplus.com
www.Vidyarthiplus.com
www.Vidyarthiplus.com
www.Vidyarthiplus.com
www.Vidyarthiplus.com
Example:2.5
The costs of two units at the busses connected through a transmission line are (with P1and P2
in MW): IC1=15+0.125 P1; IC2=20+0.05 P2, If 125 MW is transmitted from unit-1 to the
load at bus-2, at which the unit-2 is present, a line loss of 15.625 MW is incurred. Find the
required generation for each of the units and the power received by the load when the system
lambda is Rs.24.0 per MWHr. Use Penalty Factor method.
Solution:
With unit-2 not contributing to the line loss, it is due to the unit-1 alone, and hence,
2 2 -3 -1
DPL/dP2 = ITL2 =0; where, PL=B11P1 ; i.e., B11= PL/ P1 = 15.625/1252 = 10 MW
-3 2 -3
Thus, PL=10 P1 so that dPL/dP1= ITL1 = 2(10 ) P1 MW
Hence we have,
www.Vidyarthiplus.com
www.Vidyarthiplus.com
www.Vidyarthiplus.com