Kothari
Kothari
3.1 INTRODUCTION
Electricalenergyc4nnotbe sto , but is generatedfrom natural sourcesand delivered as demrand
arises.A transmissionsystemi usedfor the delivery of bulk power over considerabledistances,
and a distributionsvstemis for local deliveries.As depictedin Figure 3.1, an intercornerct€d
power systemconsistsof mainl three parts: the generatorswhich produce the electrical energy,
the transmlssionlines which it it to far away places, and the loads which use it. Such a
configuration applies to all i nnectednetworks (regional, national, int€rnational),where the
qumberof elpmentsmey vary. he transmissionnetworks are interconnectedthrough ties so that
utilities q3n exchangepower, s reservesand render assistanceto one anotherin times of need.
Generator
Tiansmissionlines
Bus2
Figure 3.1 A simpl configuration powersystem.
of an interconnected
Since the souucesof en are so diverse (coal, oil or gas, river water, marine tide, a
radioactivematter, sun power), e choice of one or the other is made on economic,technicalor
geographicbasis. As there are few facilities to store electrical energy, the net production of a
utility (generationplus the i cws over its ties) must clearly track its total load. For an
interconnectedsystem,the fu mental problem is one of minimtzing the source expenses.The
economicdispatchingproblem is to define the production level of each plant so that the l;otal
cost of generationand transmi ion is minimum for a prescribedscheduleof loads.
. Forecastingincludes ining the peak rate of supply (power demandand volume), i.e.
energy dernandfor long-terminvestmentdecisionsand short-termoperatingdecisions.
131
_-_-J
t32 PowerSystemO Zation
O RAilNG COST
3.2 GENERATOR
The majority of generators I extant systemsare of three types-nuclear, hydro, and fossil (coal'
oil or gases).Nuclear p tend to be operated'atconstantoutput levels and hydro plants have
essentiallyno variable ting costs. Therefore, the componentsof cost that fall under th-e
category of dispatching dures are the costs of the fuel burnt in the fossil plants. The total
cost of operationincludest fuel cost,costsof laboua supplies'andmaintenance. Generally,costs
of labour, supplies and ma ntenanceare fixed percentagesof incoming fuel costs. Figure 3.2
shows a simple model of a fossil plant. The power output of the fossil plant is inpreased
sequentiallyby openinga of valves at the inlet to its steamturbine. The throttling lossesin a
valve are large when it is j st openedand small when it is fully opened.
Fuel input
As a result, the operating of the plant has the form shown in Figure 3.3. For dispatching
purposes,this cost is usuall approximatedby one or more quadratic segments.So, the fuel cost
curve is modelled as a q ic in the active power generation
Fi(Ps)
u)
&
ct)
o
O
bo
tr
cE
k
C)
a
o
'p-in
8i OutPutPower -pmu
(MW1 8;
Figure3.3 Operati
^ llll0
costs of a fossil-fired generator, f 8i
r
anct Pil'* arethe lower
and the er limits on its outPut.
Econonic Inad Dispatch of rhernml Generatilry IJnits r33
The PTin
6i
is themini um loadinglimit belowwhichit is uneconomical (or may be techni-
cally infeasible)
to operate e unit and P;:* is the maximum output limit.
The fuel cost curve y have a number of discontinuities.The discontinuitieseccur when
the output power has to extended by using additional boilers, steam condensers,or other
equipment.Discontinuities a lso appear if the cost representsthe operation of an entire power
station,so that cost has di tinuities on parallelingof generators.Within the continuity range
the incrementalfuel cost ay be expressedby a number of short line segmentsor piece-wise
linearizations.
3.3 ECONOMTC
DIS ATCH PROBLEMON A BUS BAR
Let us assumethat it is k n a pnon which generatorsare to be run to meet a particular load
demand on the station. S there is a station with NG generatorscommitted and the active
power load Pp is given, the I power generationPr, for eachgeneratorhas to be allocatedso as
to minimize the total cost. optimization problem can thereforebe statedas
NG
Minimize F(Pr) = I F;(Pgr) (3.2a)
i=l
( N G
L(Pei,1,)= F(Ps)+ hl Po - I (3.3)
\ i=l
{9 =zarPs,
+ bi (3.7)
dP,,
8i
-- L-u, ( i = I , 2 , . . . ,N G ) (3.e)
2o,
y+3 -Pp
i=l LAi
or
)"=
NG
(3.l0)
I
i=l' 6
Thus, )" can be calculatedusin Eq. (3.10)and Pr, can be calculatedusing Eq. (3.9),Now
consider the effect of the ge tor limits given by the inequality constraintof Eq. (3.2c), If a
Economic Load Dispatch of Thermal Generating units 135
particulargeneratorloading Pr, reachesthe limit P#" or Pffx, its loading is held fixed at this
valueand threbalanceload ir sftaredbetweenthe remaining generatorson an equal incremental
costbasis.
where
hf* =,Psi - Pf* ( l = l, 2, ..., R1 upper bound violations)
hY" =P#n - Pr, (j= l, 2, ..., Rz lower bound violations)
Determine the new which is original Pp minus the sum of fixed generationlevels,
*'ft"r,
Pff* - Pp-
i=l
The new demandis alloca to various committed generatorson an equal incrementalcost basis.
EXAMPLE 3.L Two units of the systemhave the following cost curves:
F(Ps)= o.o5Pg1
+ 22Pn+ l2o Rs/tr
F(Pe)= o.o6P?,
+ l6Pr, + l2o Rs/tr
whereP, is in MW. Both units operateat all times and the maximum and minimum loads on
each unit afe 100 MW and MW respectively.Determine the economic operating scheduleof
',
the plants for lsads of 80 120MW, and 180MW neglectingthe transmission line losses.
Solution Using Eq. ( .10)to calculateL
22 16
Pn* -.------- + -
--1 2 x 0 . 0 5 2 x 0 . 0 6Rs/IVIWh
(i)
f
,136 Power Sys'ternOptimizat
When PD - 80 MW
Substitutingthe value of Pp in . (i) and solving for h, we get
)" = 23.6368 Rs/]vfWh
Using Eq. (3.9) to calculateg lons
I
n = -23.6369 - 22 23.6368- 16
_- 6.36MW; tDg z - = 63.64MW
51 1lx 0.05 2 x 0.06
But Ps = 16.36IvtW < 20 MW fix Pst at the lower limit 20 Mw and the rest of the demancl
will be met by the secondg
So Pr, and Pr, =80-20=60MW
). = 25.81B18 Rs/IvIWh
From Eq. (3.9),
-22 2 5 . 8 1 8 1-8 1 6
tDr s _
r =2 5 . g l g l g _ .' , 1818MW
= ; Psz= = 8 1 . 8 1 8M2W
2*o-05 2 x 0.06
When Po i 180 M'W
Substitutethe value of Pp in Eq. ( ) and solve for ),, i.e.
), = 29.0914Rs/lvIWh
From Eq. (3.9),
- 22
29.10914 29.0914- 16
P
or
= -- -1f\
9142MW; 'D9 2 - = 109.095
MW
2l x 0.05 2 x 0.06
But Pez= 109.095MW > 100 , so fix Pr, at the upper limit 100 MW and the resr of rhe
demandwill be met by rhe first g
EXAMPLE 3.2 Inrcrementalfuel ts in rupeesper MWh for a plant consistingof two units are
40 30
P P+
),= 0 2 * o A Rs/IvIWh
1 1
0.2 0.4
-- Po +275 Rs/rvIWh
7.5
So, there are three ng conditions for L, which are obtained from minimum and
maximum operatingvalues, i.e
(r) 38 < X"s 46
(ii)M < 2vs75
(iii) 75 < r, s 80
The range of load can be obtainedas given below:
(i) For 38 < )" S 46, the lower power limit of unit I is violated because ), = 46 Rs/IvIWh
when A*in = 30 MW.
The mininnumdemand the system that can be met will be, Pp - P,'"tn+ Pfin , i.e. Po =
3 0 + 2 0 = 5 0 M W w h e2n" = 3 8Rs/lvtWh.For 1,= 46, the value of Pp can be determinedfrom Eq.
(i), i.e.
Pn = 7.5 x 46 - 275 = 70 MW
So, the range of demand 50< PD<70.
beyond50
fu Pp increases the load increments unit 1 is
are placedon unit 2 because
fixed at 30 MW (minimum ng limit). So, the incrementalcost equation of unit 2 can be
written in terms of Pp as
h=0.4x(Po-30)+30
On simplification
L = (0.4 o + 18) Rs/TvIWhfor 50 < PD < 70
i'
(ii) For 46 < )" 3 75, the ation of generatorswill be within the operatinglimits' When
- 75, the value of Pp can be ined from Eq. (i), i.e.
= 7.5 x 75 - 275= 287.5MW
Pp
Rs/IdWhfor 70 3 Po S 287.5
- 75 Rs/lvIWh
(iii) For 75 upper timit of unit 1 is violated because )"
when 1. = 80 Rs/
when Prnax= 175 MW. The max mum demandof the systemthat can be met
MWh, iS
- 175 + 125 = 300 MW
P n = P m a x * P2'"*, i.e Pp
Thble3.1. Generationschedule(Example3.2)
NG NG
Pt=7 ZPs,BijPsi MW (3.13)
i=l j=l
P* and Pr, *" the real po injections at the ith and 7th buses,respectively
Bii are the loss coefficien which are constantunder certain assumedconditions
NG is numberof generati buses.
The transmissionloss f, la of Eq. (3.13)is knownas the George'sformula.
140 Power System Optimi rcn
where
Pr- and P, ue the real injections at ith and 7th buses,respectively.
Boo,Bio,and B,; are the coefficientswhich are constantunder certain assumedconditions
NG is numberof generati buses.
The aboveconsffainedop mization problem is convertedinto an unconstrainedoptimization
problem.Lagrangemultiplier thod is used in which the function is minimized (or maximized)
with side conditionsin the form of equality constraints.Using Lagrangemultipliers, an augmented
function is defined as
l r 'L ) = F(Pe,'t
* t(r, +Pr .r,) (3.1s)
i:
where ,?,is the plier.
Necessary le opti mization problem are
t(Pr,
)
. ^(W , )= o ( i =1 , 2 , . . , N G )
)P,,
aF(P8
( i = 1, 2 , . , . ,N G ) (3.16)
dPs,
where
aFest)
is the incremen cost of the ith generator(Rs/IvIWh)
at;
dP,
= u is the irrcremental missionlosses.
dP,
Equatipn(3.16) is known the exact coordinationequation,and
L(Pe,,L) NG
-=' P
I )n +' 'P
Lt P r ,= Q
ah ?2 r ' 8 i (3.17)
)
=Za;Pr, + 6i ( f = 1, 2 , . . . ,N G ) (3.le)
=h
(3.20)
NG
ZaiPsi + bi t n,lijPe)
,>, 22 ( i = 1 , 2 , - .N
.{G
, )
,
Ji r= l )
NG
t
)
2(ai + LBii)Pe, Bi,
iPt
iPci l-bi ( r = 1 2,
, ,, NG)
j=l
j+i"2r )
'['
P$ ( l = 1, 2 , . . . ,N G ) (3.2r)
Z(ai + )"9 ii)
If the initial values of b,(i = 1,2,..., NG) and 2uare known, the above equationcan
be solved iteratively until . (3.17) is satisfied by modifying )". This technique is known
as successiveapproximation. The stepwise procedure is explained below. For simplicity it is
considered that the solution ins within limits.
142 Power System
EXAMPLE 3.3 The tuel i per hour of two plants are given as
Fz(P' ) = (0.0074r
P:r+ 10.833
Pn +240) Rs/h
Determine the economic to meet the demandof 150 MW and the correspondingcost of
generation.The transmission are given by
+ o.wzr!, - 2 x 0.w2PstP82
0.001Pr2
10.333 10.833
1 5 0 +- - L -
2 x .00889 2 x 0.00741
Compute P^ and Psz usingEq. (3.9)
11.81812- 10333 =
i P8 83.5276MW
2 x 0.00889
r 1.818t2 - 10.833=
Pe 66.47239MW
2x O.O074l
144 Power System Optimi
Step 3: Set IT = I
Step 4: Compute Pr, Prr, usingEq. (3.21)
1 r . 8 1 8U
12 - 10.333
0 x (-0.0002) x 66.472391
Pgt = = 43.44557MW
2t0. + (11.818r2
x 0.001)l
"11.81812U x (-0.0002)x $.5n61 - 10.833
tDg z - - 22.22453MW
210. 4r+(1r.fr8r2x0.002)l
Step 5: Computetransmi ion loss,
P t = ( .001X43.44s5T2+ (0.002)(22.224fi)z
_ 2.489154MW
2(O.OO02)(43.44ss7)(22.224s3)
Step 6: Compute the balancerequirementto be met
2
A , P = P p + p r - Z Pr, = 150.0+ 2.489154- 65.6701=+ 86.81905MW
i=l
Table 3 2 Generationschedule(Example3.3)
IT Pe, L LP PL
(Mw) (Rs/IvIWh) (lvIW') MW)
I 43.4456 2.2245 1r.8t8t2 + 86.81905 2.489rs
2 I 19.1606 0.5657 t6.1,5907 - 18.93145 20.79482
3 110.l3t6 7.4586 15.21250 -g .3316A4 19,,25956
4 r0r.7692 1.8196 r4.74592 r.894872 t 5.49376
5 102.7077 2.1681 14.84067 + 0.848730 15.72457
6 r03.4826 2.6953 1 4 . 8 8 31r -0.203037 15.97491
7 r03.3739 2.6494 14.87295 -0.077766 15.94553
8 rc33Ar2 2.5993 14.86907 + 0.021338 15.92185
9 103.3133 2.6048 14.870r3 + 0.007057 15.92521
l0 103.320r 2.6096 14.87049 -0.002229 15.92743
11 103.3188 2.6089 14.87037 - 0.000624 15.92705
t2 103.3181 2.6085 t4.87034 + 0.000239 15.92684
l3 103.3183 .6086 14.8703s + 0.000041 15.92689
EXAMPLE 3.4 For a enerator system, the fuel cost coefficients and the operating
generatorIimits are givenin le 3.3(a).The B-coefficients
for transmission
loss aregivenin
Table 3.3(b). Determine the ic schedulefor loads 160 MW and 210 MW.
Solution Algorithm 3 is followed to get the optimal generation schedule. The achieved
generationscheduleis given Table3.4. The numberof iterationstakenby the algorithmare
given in Thble 3.5. The val of steplength a chosen,and the achievedconvergence are given
in Thble 3.5.
The method is very sen ve to the value of a,, i.e. the step size. The number of iterations
dependupon the assumedval of u. An incorrect valueof a,, sets the solution procedurein the
oscillations.To avoid this pro lem the Newton-Raphsonmethod can be implementedto get the
solution.
Generator a; Ci
r (R$/]WW2h) (Rs/h)
I 0.00608s 10.04025 t36.9r2s 5.0 150.0
2 0.00s915 9.t60576 s9.1550 15.0 100.0
3 0.005250 8.662500 328.t250 s0.0 250.0
PD Pe F PL
Mw) O4w) (Rs/h) (Mw)
No generationlimits imposed
160.0 57.5577 7 .5238 37.9172 11.09701 2176.023 s.998648
210.0 83.4010 9 .6169 39.4862 1r.52315 2741.473 8.503935
Generationlimits imposed
160.0 53.3906 .6094 50.0000 I 1.08013 2179.r59 7.999945
210.0 79.9043 .5531 50.0000 11.51164 2743.905 10.457380
(Exarnple3.4)
3.5 Optimalschedule
PD Iteratiorts a
Mw)
No generation
li its imposed
160.0 20 0.7009506X l0-4 0.00s
210.0 20 0.7820129X 104 0.005
Generation
limits imposed
160.0 l3 0 . 5 5 3 1 3 1X l l0-4 0.00s
210.0 T4 0.2t93451X l0-4 0.005
a2L
" = L P ,
:rD r L ) " =- a L (3.22)
. O t L dP.*,
8,
#
HffiLPsi-.#L)'= (3.23)
Der,frVStr )s cal
lval ives ca be
)e(
f;r:Y:lt*1=[-r;1
lows:
(3.24)
AL _ dF, NG
^[ 2a1Pr, + I zB,Ptt
IP,, aP_
aP_ aP- o,
j=l
AL NG
I P*, (3.2sb)
AL i=l
( i = 1 , 2 , . . . , N G ;j = 1 , 2 , . . . , N G ;i * j ) (3.26b)
.--c
Economic Load Dispatch of Thernnl' Generating Units 147
Takingderivatives
of Eqr. (3. ) and (3.25b),withrespecrto i.,
d2L d2L
aTaPs aPstaA
*,{ 2Brip,
-l
#-r-4,
-d2t
-- =Q (3.26d)
a1'
Equations(3.22)and (3.23)[or . t . (3.24)l are iteratedtill no further improvement is obtained,
or
single derivativeswith respectto control variablesbecome zero. The stepwise procedure is
outlinedhere.
,[arnJ
IHll l=-tJl
L^^)
Deactivaterow and col mn of Hessianmatrix and row of Jacobianmatrix representing
the generatorwhose ge tion is fixed either at lower limit or at upper limit. This is
done so that tixed g tors cannotparticipatein allocation.
6. Gausseliminationmeth is employedin which trian_gularization and back-substitution
processesare perf, to find Pr,(i = L,2,..., R and A,l.).Here R is the numberof
_qeneratorswhich can icipate in allocation.
= Ps,* LPr,
,P,,ilt* (i = I ,2, ...,R and inew- ), + Ltr)
14E PowerSystem
EXAMPLE 3.5 Determine economicscheduleto meet the demand of 150 MW using the
Newton-Raphsonmethod. Use the data of Example3.3.
Solution Giv6n PD = 150 MW. Initial values are presented below as calculated in
Example3.3.
a2L
0
*
The Jacobian matrix ele ts are computedand are shown below:
dt
+ )u(2BrrPr,* 2B21Ps,- 1) = 1.66N47 Rs/IvIWh
+2822Pr,- 1)=2.747459Rs/lvIWh
bz + 1"(zBzLPs,
Economic Load Dispatch of Thermat Generating
Units 149
P:, + 2 B z P r , P g =
, I 3 . 5 9 3 1M
I W
.767522
.066314
IT PIl )" AP PL
(Mw) (Rs/LIWh) (Mw) (Mw)
2 103.2330 2.56359 r4.86603 0.1043829 15.9020r
3 1 0 3 .134 0 2.60577 r4.870r2 0.0578308 15.92552
4 103.3183 2.60843 14.81034 0.0031089 15.92682
s 1 0 3 .r38 3 2.60858 14.87035 0.0001812 rs.92690
= 2309.771Rs/h
Power System Opti
AzL
* a z L -'
M=_aL (3.27)
lLPr' arr.M" dPr,
.[#-')^^=-
[#.^(# '))
or
=[, E:)e+^^)-t#)
^W)*,, (3.2e)
Let
AF, -
# 2atPs,+ bi Rs/TvIWh (3.30)
ut gi
azFi ^ d2P1-=Zar+Z)"Bii
+ IW Rs/lvIWh (3.31)
q
= )" + L)', Ki =
Anew (3.32)
#
FromEq. (3.29),
( 1- K ) L * -Xi
LPr,=T (3.33)
NG
^ , P s i = - [+* P r
j=l
#-r)
.*d
Economic l-oad Dtspcttchof Thermal Generating [Jnits 151
NG
j=l
r(l- K)Mr, =P; (3.34)
NG
where P; = P p + P y LPs, (3.3s)
i=l
(1 - Ki)L""*
I (1- Ki)
j=i Yj
NG (1 - Kj)Xj
P; + I
j=L Yj
(3.36)
2. Computethe initial v uesof PsiQ-- 1,2, ...,NG) and)" by assumingP1= 0. The values
of ,1,and Pr,(i = l, be obtaineddirectly using Eqs. (3.10) and (3.9),
respectively.
3 . Set iteration counter, - l.
EXAMPLE 3.6 Determine he economicscheduleto meet the demandof 150 MW using the
approximateNewton-Raphs method.Use the data of Exarnple3.3.
Kr =2&1Pg,+2BpPr, - 0.1404663
Kz =ZBzrPgt2B22Pg,= 0.23324785
Calculate ffansmission I
_B
Py ,P!, + B22P!,+ZBpPrrPr,= 13.59311
MW
- 2arPr,+ b, - 11.81812
Rs/lvIWh
= 2a2Prr+bz = 11.81812
RsAvIWh
= 2(ar + hBLr)= 0.041416Rs/IvIWh
= Z(az+ )"Bz) = 0.06209Rs/lvIWh
14.81919Rs/lvlWh
2 | - K;
t '
iir Yj
IT Pg h AP lLnew- 11
(Mw) (Rs/ivlWh) (Mw) (Rs/It'IWh)
dFi
(3.37)
dPr,
NG
ZPs,-Py+Po (3.38)
i=i
The initial values po and A.oat calculatedby assumingthat the transmissionlossesare absent
(r.e.Py- 0) and usirfbeq. (3.10 and Eq. (3.9).
The nonlinearequationsin P* andX, canbe solvediteratively.L,et P:t(i = 1,2,..., NG) and
Lo be approximatedsolutionsto Eq. (3.37). To find the new approximation,let
Lnew _ 9+A,L
pnew - (i = l, 2, ...,NG)
'gi + LPr,
8i
( i = 1 , 2 , . . .N, G ) (3.3e)
rr- -2
t54 Power System O
and
NG
t*:'*- P/'* * Po (3.40)
,l
Taylor's expressionto first r of Ple* is
Plt* = Pot^
* H d P | . M*'
np (3'41)
,.?,3,po
where Pf is the initial tra ission losses.
loss from the initial valuecan be
The modified value f the incrementaltransmission
obtainedA S
dpT* azPottp_
-= d P o , *' Y -''sl
dPr, dPr, i"=, aPsiaPsj
To retain the classicaform and at the sametime to improve the convergence,it is possible
to include only the trh it of the summationof the aboveequation,i.e.
Similarly, the modified i tal cost from the initial value can be obtainedas
(P;.* - (3.43)
t { p *-w * ( a r ? - a'r?-oo)=
-='v7 ^ " * (Pr -) - (Ao + LA,)
aP? s' [Er*, aPr' r' ) (.- dPr, )
aPg)
+ (w
)new - ; - q p :- - ^ o (3.44)
8i
[dPt'
aPz
8i
8, a'.)
^ 1 dt Po'
Here the term un
6rtgt ! is ignored,being small.
Economic Inad Dispatch of Thermal Generating Units 15s
SubstitutingEq. (3.41)ir
NG
i=l
By separatingthe Prl"*'
Ncr
>l (3.4s)
i=l \
(
tDg+i -- l l rr
\
w h e r eP i = P p + P f - Y
a4
i=L aPr
Substitutingthe value of
I 1,,r_ -ll
l + + L o azI
dp
\- si - inew (3.47)
(
1 1 -
\
The aboveEqs. (3.47) ) can be rewritten of equalincrementalcost,
when the transmission of considered,as
(3.48a)
(3.48b)
(3.48c)
Power System O
(3.48d)
P; - P D + P tf P " :'
iJ dpr, si (3.48f)
In view of the defi problem in Eq. (3.12),the problemin Eq. (3.48) is redefinedas
zai f|, + bi = Lnew (3.49a)
and
NG
P;, = P;
i--L (3.4eb)
where
NG
= 2, ...,NG)
j=l Q.a9c)
* ar + ho Bii
a;= (3.4ed)
(1- KT,
,* bi - 1oBiiP!,
D ; =
(3.49e)
I l - Kr. 0,
p*
'8i - (1 :P)
KI r Til'*
)P (3.4ef)
\IC
NG N
p 0 -
t L \ IP!,
I >, r, 11.
B ..Po
BuFi
I gi
f = l j=
;- 1 Q.aee)
NG
D * -
t D - PD +r P f - IIr(,9
r P:,
i == ll
(3.4eh)
The problemin Eq. (3. 9) can be solvedusingEqs. (3.9) and (3.10)
as explainedin Section
3.3. The stepwiseprocedure implement the efficient method is outlined
below (Algorithm 3.5.).
Algorithm 3.5: Economic ispatch Using Efficient Method
1. Read,data,namelya br, ci (costcoefficients);Bu(B-coefficients)
(i = l, Z, ..., NG; j = l,
2, ..., NG); e (conv )nce tolerance);and ITMAX (maximum
allowableiterations),etc.
*-d
Economic Load Dispatch of Therrnal Generating units
2. Compurethe inirial ues of Ps,(i- 1, 2, .-., NG) and A, by assumingthat pt= 0. The
valuesof 2" and psie , 1,2,..., NG) can be obtaineddirectlyusing
Eqs.(3.lOiand (3.9),
respectively.
3. Set iteration counter
4. ComputeKl, ai, bi, , ilrd Pj usingEqs.(3.49c),(3.49d\,(3.49e),(3.499),and(3.49h),
respectively.
5. ComputeL** , Pi,, PrI* usingEqs.(3.10),(3.49a),and (3.49f),respectively.
6. Check I Lo* - Lo I s if 'yes' thenGOTOStep8.
7. Modify Pro = Ps:"* (i 1,2,...,NG)and1.0= i.,"*
IT=IT+lGOTOS 4 and repeat.
8. Compute the optimal I cost, and the transmissionloss.
9. Stop.
The limits of generationscan added by the procedureoutlined in Section 3.4.
EXAMPLE 3.7 Determine economic scheduleto meet the demand of 150MW usingthe
Efficient method. Use the data f Example3.3.
ho = 11.81812
Rs/IvI , Prl = 83.5276
MW P:, = 66.47239
MW
Find the incrementaltran ssion loss with initial values of generation,i.e.
8=2BrrPrl*znzzP!,= 0.233247gs
Calculate transmission
= 1j * 7!!r: =o.o11oze
Rsrrrw2h
(1 - Kl)z
az * LoBn
= 0.232479RsfirdW2h
(1 - KD'
bt - to nnP!,
=9.724716Rs/IUWh
1-r(f
158 Power System
bz - lo nzz$
= 10.02016RsA'IWh
,1 - v t t 2
0
z
- P p + P f t rc?P!,
= 136'40690
IvIw
, i=l
2 u l
P;+L3
i =t l Zai
LUi
=-=l4.8 l9l9Rs/IVIWh
2 I
I, o *
i=l LAi
Calculate new ge NS AS
*
P,, =Ln*
- bi
=90.84408MW ; = hnew- b;
Ptr, = 45.5298IvIW
zai : zai
The required new mod ed generationsare
Pi, Pi, =
P;"* = 105.7283
IvtW; P;:* = I\dW
59.32055
I - ,Kto r-K8
Checklrl'* - Lol= J. 107> 0.001
For iterationIT = 2 ::
ru Psl L P L LP lLnew- Ll
(Mw) (Rs/IaWh) ' (MW) (Mw) (Rs/lvIWh)
.-J
Load Dispatch of Thermal Generating Units
c)
o
9l 1,,
.9
o
N.
Select the reference bus and unt elementsto bus L to be part of the network. Form fr-bus
using the Z-bus algorithm and zero row and column in order to obtain (NL+NG+I) x
(NL+NG+l) matrix, where NL is number of loads, NG is the number of generatorsand NB is
the number of buses.
Es = \Is (3.s0)
where /s is a vector of load and r currentsand
The voltage vector transform into the new reference Frame I (Figure 3.5) by the conjugate
of Tr, i.e.
frunspose
I 'Eo=(f\r4T)lr=ZJr
I h-(ri (3.51)
I
Iwhere
I (3.s2)
4Tt
0
I
I
pvith
i
II
i
II
L-- _1
In
+
7
#Lg , B- l )
( k= 1 , 2 ,. . . N (3.53)
| - L; ( k = l , 2 , . . . ,N B - l ) (3.54)
NB -I
\ Loto+t, (3.ss)
k=l
Now,
(3.56)
zu
where /n is the neuffal cu
f1 is complex calcul off nominal turns ratio betweenbus k and the reference bus.
21 is calculated with the fo owing rule [Kirchamayer,1958]:
where
ft is the complex
/7 is the total load c
Economic Inad Dispatch of rhermal Generating (Inits
€ 1 -t p , = z t ; @ * - t * e , ) (3.60)
k=l
Transmission
lossis as
Pt+ = (Ihrq - Kl;r;fzrlrrlr)l = (r)r4rz (3.61)
=l'o
I
l
LO J
I = llt, 12,...,
/Nilr
irnplement on compute\ 22 be obtainedas (Figure 3.6)
a
4- (3.62)
21
l i r
.+
et- tF,
1,,
NL NL
w=>zt ztT,ilt
i=l j=L
NL
oi= Z fi zr&,j) l, 2, ...,NB + 1)
/c=l
. . J
,.ry
i
162 PowerSystemO
NL
b,= 2 Zre,k)I* (l = NL * 1,..,NB + l)
&=l
where
Iz = T{t, Ez = h. Iz
Q.e)
with
/ = [fNl* y 2, ..., fNI-*Nc*t]
. t*io I
( k = 1, 2 , . . . ,N G ) (3.6s)
tk= - , INL+NG+I - - - F
ty ty
NL
tL= \ forr
k=l
2 3 - l-t (3.66)
where
'na^
Zt(n, m) = TrW; - brT; - , ) ( n = 1 , 2 , . . . ,N G + I ; m = 1 , 2 , . . . ,N G + 1 ) ( 3 . 6 7 )
* Z 1 ( nm
INL+l
L+t- ',
,,
En=-2
tt
A
conomic Inad Dispatch of Thermal Generating Units. 163
It is assumedthat the real and tive powers at each generatorbus are linearly related as shown
in Figure3.8.
*.
€*t (k = NL * 1, ...,NB) (3.68)
rfi
= & m P k +F * (m= k-It[L) (3.6e)
€1a
Qr
Q * = Q o r +s f r
l-{
0)
'
o
P.
(.)
o
Cg
a)
&
Real power
where * ftuo]t
= la1P1* Fr,...,aNcPNc
I
3rj=R" lrwi aj ztL,i) ait 4(i, j)l *o' i = r,2, ...,NG)
(i= 1,2,..1' (3.72)
[ -l
Bro= Reldt"z't(i,NG +' 1) a,ti zt(Nc + l, i)l (3.73)
L.*-
l& Power System )pt
EXAMPLE 3.8 Use the clas method to determine the B-coefficients for a 5-bus system
shownin Figure3.9.Bus 5 is as the slack bus. The seriesimpedanceand line charging of
each line is taken liqn Table
0.08+70.20
I .<f
o
O
o
+ +
N
O C\
O
o
0.02+70.06 o
4 3
Flg 3.9 Powersystem(Exampleg.g).
Table 3.9 Power systemnetwork data (Example 3.g)
The total line charging ittance to ground at each bus is given in Thble 3 . 11 .
I j0.0s5 - jl8.1818
2 j0.020 -jsO.0
3 y0.055 -jl8.1818
4 j0.04 - j25.0
5 jo.M - j2s.0
The bus impedancealg hm is applied to obtain the Zo for the power system shown
in Figure 3.10. All line chargin are fixed to neutral and treated as bw Z (common point). The
bus 5 is a reference bus. The to build the Za-busis given in Table 3.12.
R
0.08+/0.20
oo
@
From Eqs. (3.53)and (3. , the values of ),1,and tp arraobtainedand are givenin Table3.13.
)"k tp
-_-{
Ecorwmic Load Dispatch of Thermal Generating Llnits 167
The load current is defi AS the sum of bus injection culrents.Using Eq. (3.58)'
Ir -r.426628- i0.903495
I1 = -0.524929- i0.3646r2
I2 -0.4284r - i0.241396
,: = -0.473258 - i0.297487
UsingEq. (3.57),we
f,
- j0'016095
Ir = + = 0.378144
rT
+ j0.014g75
- 1 9 7 = 0.290831
b=
lz=
- 1! 7 = 0.331025+ j0.001116
The ratios I; are used in trans tion T2.Using Eq. (3.62), 4 ls obtainedas
.,7
z4-
L
168 Power System
':
r1<--
Network
NL
0R-
r
Ltt (3.75)
i=l
NL NL
Iq=I
d .tr
- €iIt
vl(
(3.76)
i=l i=l
^SR
.+ (3.77)
Ip
Ynr
.si
tt (3.78)
€ReR
f -
IBUS - Esus (3.7e)
Load Dispatch of Thermal Generating Units
+iNL
Network
Figure3.1 Remoteequivalentindependent
of a network.
wherc
tu[-NL+NG+2
/sus = bus injection = [itr.., fp;, tL, jNL*I,...., iNr*Nc, inJr
Y22is partitionedmatrix f s i z e ( N G +I x N G + 1 )
Y21is partitioned matrix f s i z e ( N G + lx N L + 1 )
Yp is partitioned matrix f s i z e ( N L + 1 x N G + 1 )
hr is partitioned matrix size(NL+lxNL+1)
Another admittance Y3 is consideredby eliminating bw R elementsfrom the admitt-
ancematrix of Eq" (3.80). inverseof Y3is the required impedancematrix, i.e.
4= Ytr (3.81)
Re(Ir,) = 2,..,NG)
,Solution Performing I flow, the voltageand power at all busesare calculatedand given
in Thble 3.14.
P+JQ
I -0.50 - j0.2s - j0.r2s030
0.865666 -0.524929+ j0.364612
2 -0.40- 70.15 - 70.133531
0.858383 -0.42844r+ j0.241396
3 -0.45 - j0.20 - p.L2ss07
0.871962 -0.n3258 + j0.297487
4 - 0 . 1 s- j 0 . 1 0 - p.008090
0.984324 -0.151544+ j0.102838
L= -0.653587 + j0.326793
= -0.530f,/;6
+ j0.198767
= -0.579845+ j0.257709
= -0.154805+ /0.103204
IJsingEq. (3.78),
= 1 . 9 1 7 8 8-1 j 0 . 8 9 5 0 1 1
lfhe network will be as shownin Figure 3.13.
ifhe system admittance x is obtainedusing the l/-bus algorithm as
v
IBUS -
^d
Economic l-oad Dispatch of Thermal Generating Units 171
Ynr
FromEq. (3.80),the
+ 70.
lo.o+ozl5 + i0.027e74f
66 0.034555
4- + 70.036764)
fo.orotts+70.74 0.M6476
Using Eq. (3.82), B-coeffi ients are real valuesof Z matrix, i.e.
Bn = 0.M1'5240 B 0.035t046
Bzt = 0.035t046 0.0464764
-. l-v, vil--
sr=v,LTl (3.83)
where
'Z^ is the impedanceof branch connecting the ith and 7th buses
% i s the voltage of the i bus
7; is the voltage of 7th
'Ihus,
= * t l v ? l - l v ,l l v ,|/ ( o i - o i ) l
where
llv; I is the magnitudeof Itage of the ith bus
llv; I is the magnitudeof oltage of the 7th bus.
Therelore,
','.1- 14 V;
| ll I {cos(0i - 0) + 7 sin(0i - gj)}] (3.84)
Economic Inad Dispatch of Thermal Generating Units
R^+iX^
= { i Ix * =Rm
*+iX^
W
s,Su
..= = l y r i l v , t ( c o (s0 i - 0 , ) + 7 s i n ( 0 i - 0 i ) ) ] ( & , + i x ^ ) (3.85)
_Tllv,'l
The real part of Eq. (3.85) is activepowerflow from the sendingbus i to thereceivingbus7,
i.e.
P^=Pi t v i z t - t v , n v j l ( R - c o s( 0 i - 0 ) - X - s i n( 0 i - e j ) ) ] (3.86)
Assume
(i)lV;l=lV1l=1
(ii) X*D R*
(iii) (0i - 0) is small such cos (9, - 0) = 1 and sin (0i - 0i) = ei - ej
1 - -
P^- Pri= A LRrn
tLml
o,-oi (3.87)
x^
Generalizing,we have
) (3.88)
J
3.10t.2PowerLoss in Line
The active power flow from e sendingbus 7 to the receivingbus i can be expressedin the
similar way as Eq. (3.8O,i.e.
= f f |i | (,l, , Y
,,
t +l t v j l 2 ) R ^- 2 l v i l l v j l R , c o s( 0 i- e ) l (3.e1)
Equation(3.91)is rewrittenas
f z n ^ - 2 R ^ c o s ( g- 0; ) l ('.' It v i I = l V i l =11)
lz^l
2R
Rk+ (
c Q = =')z, s
1 -' ccos ,in2
$ r)
-E
4R^
l x* R^ 2',
on 3.10.t
l _ , Sectior
_4R^
E t'.' (0i - e)n is small,Section
3.10.2f
_ R .(
I\
The above equation can be tten as a function of real power generation by substituting
Eq. (!i.87) into it. Thus,
A(m, i) = =*t#
+=#(+) #)
From Eq. (3.89),it is concl that and
W=xii W=xr,r.rhus,
A(m,rr=+ (3.e3)
where,
P, is the real power on transmissionline m from sendingbus 7 to receivingbus k.
X7i and Xpi are the el of the X matrix.
X. is the reactanceof li m.
Economic Inad Dispatch of Thermal Generating Units 175
In Eq. (3.93), since all changesare compensatedby the reference bus, the total
systemgeneration
is assumed
t be unchanged.
NG NLOAD
where
litlG is number of genera
I\ILOAD is numberof I
Using GSD factors, the n line flows, after rescheduling generation, can be expressedas
NG
+) A(m,i) Ms, ( m = 1 , 2 , . . ,N L ) (3.e4)
i=l
where
P: is the basecase line fl
I,IL is the ntrmberof lines in the system.
Ilv this incrennentalsu tion, any changein line flow due to changein generationcan
easily be calculated.
NG
P* D (m, i) Ps, ( m= 1 , 2 , . .N
, L) (3.es)
i=l
where
Lt(m,f) is the GGDF
IttG is number of generati buses
IttL is number of lines.
3.10.5i Derivationof GG F
From GGDF, if a particular generatork is increasedin generationby
some ne ru will be
NG
P'rn=I O i) Ps,+ D(m, k) Msr (m = 1,2,..,NL) (3.e6)
f=l
I
L
t76 Power SystemOpt nn
The changein line flow with pect to the generationchangeis termed the GSD factor.So, the
above equationcan be rewri n a s
NC
p P^- - L ot*,i) Ps, ( 3 .1 01 )
d=l
i*R
FromiEq. (3.95),
NG
when:
Pr, is final generationfi m generatori, which is now reducedto zero
Pr* is final generationfrom the reference generator which contains the total system
generation.
Equation (3.102) is red
NG
'*= D(m. R)Psn= D(m, R) Pr, (3.103)
i=l
----da
Economic Load Dispatch of Thermal Generating Units
SubstitutingEq. (3.103)into Eq ( 3 . 1 0 1 ) ,
NG NG
D(m,R) = pm- I A(m,i)Ps,
I d=l
i=l
i*R
or
NG
P*-) a(m,i)Pr,
i=l
. d*R
D(m, (3.104)
where
P. is basecaseflow of m
Pr, is generationof the i generator.
{+ = y 2B,iP,, (3.106a)
dPr, fr
Taking the derivativesof Eq. (3
(3.106b)
NL
Pt= (3.108)
m=l
178 Power SystemOptimi.
NL NG
EP
AP
= I ^ '\ tzr(m,i)D(m,j) Pr,l (3.10e)
m=l j=l
and
NL
a2p
'r,aPr, 2 R^ D(*, i) D(m, j) (3.I 10)
m=l
NL
s,j= I *,, ,i) D(*, j) (3.11)
m=l
where
,R, is the,resistanceof li
.D(m,i) are Generalized ion Shift DistributionFactors.
EXAIIPLE 3.10 Using the DF method, determinethe B-coefficients for a 5-bus system
shownin Figure 3.14. Bus 2 is taken as the slack bts. The line data is given in Table 3.15.
2 5 4
Figu 3.14 PowerSystem (Example 3.10).
L---,,- -,. -- . 2
_,.4
Economic Inad Dispatch of Thermal Generating Units t7g
llolution From the base ase load flow, the line flows are obtainedas given in Thble 3.16.
P^ '^ + jG^
PL* + jQL* tv^l fl
I 58.00 + j39.023 86.0 + j2O.0 1.015 rl.l139
2 12.78 10. + jr3.549 46.8 + /5.8 0.980 7.9999
3 -04.77 34. + jl0.528 0.0 + J0.0 0.953 15.6829
4 8.00 0. + 710.884 36.2 + /3.0 0.981 5.8130
s -33.07 0. + j25.343 17.4+ j12.0 0.972 6.3267
lb obtain the Z-bus,the Z- us algorithmis used.The data for Z-busis describedin Thble3.17.
For generatorsl, 2 and 3, the G factors are calculatedand are tabulatedin Thble 3.18.
Power System Opti
Table3.18 GSDfactors
3 . 1 1 T R A N S M I S S I O N L SS COEFFICIENTS
An exact transmissionloss fo ula has beenderivedby Dopazo et al. fl9671,usingbus powers
and sy'stemparameters.The I systemloss is the sum of the bus powers,i.e.
N B N B
Pr+jQr= Ii=lt ' v,Ii (3.r12)
i=l
Economic Load Dispatch of Thennal Generating units
NB
+j 2 tr l /, ll 1, | (cosg; sinei - coso, sin0,)
I .l=l
) NB /Ns )
z i j t l j sl i n g r l =) r / ; l c o s r r l l z i ; t t i sl i n o i l ( 3 . 1 1 6 a )
) j=t \i=t )
NB NB
Pr+ jQt= (3.117)
i=l ,l=l
-P.
f
l f - l =
' -
(3.123)
l V ,l c o s @ ,
Figure 3.16 Powertriangletor ith bus.
" --d
onomrc Inad Dispatch of Therrnal Generating Units lg3
FromEqs.(3.120)and(3.tZZ),
r{l=eil (3.r24)
Substiturting
Eq. (3.123)into Eq. ( 3 . 11 8 ) ,
S$ D 4Pi cos(o,-oi)
^'i (3.125)
* fr wltvi I t* 4'r* oI
Equation(3.125)can be rewrittenin terms of B-coefficientsas
NB NB
PL- I I Pi Bij Pj (3.126)
i=l j=l
c o s ( 0 ;- 0 i )
B,
lvi ll V, I cosQ, cosQ;
c o s ( 0 ;- e i )
er= sinp, sin@r.
(3.r27)
Q r -I I Q,cuQi (3.128)
f=l j=I
c o s ( 0 ,- 0 i )
gii
sin @,sin p,
Let NG be the number of erating busesand for rest of buses, Pr,= 0.0.Equation(3.126)
be rewritten as
NB NB
P L -I > (ps,-po)B,i
i=l .l=l
NB NB NC INc
=I L po,Bih,*Z
i=l j=l i=l
IT
L,t=t
Pr,BiPr,
184 Power System O
NB NB Nc lrvs I NG NG
= L.tT LtY Po, B,j -I
i=l j=l l I @ u + B i Pr,
i=l Lr=t
) h*I, l
j=l
I
j=l
pr,Bipr,
J
The above equation be written in the form of Kron'slossformula
NG NC NG
,-= Boo+
I B,oPr,
+f Z ,r,Bipr, (3.tzg)
i=l r'=l ,l=l
where
NB NB
Boo=I I Po,B,iPo,
i=l ./=l
NB
Bio=-I @u+ 81) pa,
j=l
et- coo
+ X o, er +f X e* c,ier, (3.130)
j=l i=l j=l
where
NBNB
Coo =
i=l j=l
NB
cio= -I (cu+ ci) ea,
j=l
The system power I are based on the assumptionthat (i) the generatorbus-voltage
magnitudesand angles are c stantand (ii) the power
factor of each sourceis constant.However,
the: use of loss coefficien can accountfor any change
in load demand at the buses while
schredulingthe generations the system.
3.112TRANSMISSION
LOSSFORMULA:F UNCTION
OF GENERATIO AND LOADS
. (3.r2s),
(di-d;)-@i-Q)l
cos@,cos@,
- 6 ; ) s i n ( Q ;-
Qi)
cos@;cospy (3.131)
Economic l-oad Dispatch of Thermal Generating Units 185
R,',
cos (6; - 6); 6ri=ffi sin(6;-6,
can be rewritte , consideringai1 and b;;, as
( l + tan Q i Q)+bilGanQi-tanQ)f
Substitutingtan Q, = * in the
ri
PL ,,,,1",['.
)iI w + bij
On simplification,
NB NB
= -
t"i1(4Pi+ QiQ) + i1(Q;P1
)i \ 48j)l (3.r32)
i==l j=7
NG
where
sQi
Ai= 6i- Qi ( i = I 2, ...,NB)
Pr=Pr,-Po, . (t 1 , 2 , . . .N
, B)
Qr= tan-'(Q/Pt)
4 and P1 are the real power injections at ith and 7th buses,respectively
Qi ffid Q are the ive powerinjectionsat ittr and7th buses,respectively
NB is the number of ses in the network
Zi = Rq + jXii (e of impedancematrix).
The above constrain optimization problem is converted into an unconstrainedone.
Lagrange multiplier meth is used in which a function is minimized (or maximized) with
conditions in the form of ity constraints.Using Lagrangemultipliers, an augmentedfunction
is defined as
, 4)= F(Ps)
+4[X (3.r37)
"-
where 4 it the Lagrangian multiplier.
conomic l-oad Dispatch of Thermal Generating Units 187
The ndcessary conditions fo optimization Problem, Eq. (3.137), btatethat the derivatives,
the
- l, 2,...,N G) and Lo, are equal to zero, i.e.
with reslpect to control/decision v iables Psr(i
aL(Pvi' ),py a + ^o
2 ( u r ,- r')) = o ( f = I , 2 , . . . ,N G ) (3.l38a)
EP,, t , [a"o
NG
h ,)
aL(Ps,, (3.138b)
D
+Pr
dho i=l
where
AL (3.139a)
"^ *-^i" *
t = L P8i" . *t y aPs,aPsi i
"6p2 dPr,
NG
dzr (3.l3eb)
j=l
M" odrr
d, P, - ( i = l , 2 , . . . ,N G ;j = 1 , 2 , - - .N
, G) (3.141)
Bii + B1i
aPsiaPs
and available
All the,derivatives for exPressi s given in Eqs. (3.139a)and (3.139b) are known
(3 141). Using Gausseliminationmethod, Eqs. (3-139a)and
from tiqs. (3.26a-3.26d)and
The
(3.139b)are solved to find chan e in controlvariables,namely PsiQ= 1,2,..., NG) and,ln-
is
controll variables are updated. ne above procedure is repeatedtill no further improvement
achieved. The detailed stePwise procedureis outlined in Algorithm 3.8-
L
Power System
Qs'=
R
16. Check +( A,r)z s e,
Ii=l ro+,)z
if 'yes' then GOTO 20.
17. Modify P;"* = Ps,*
, ( i= 1 , 2 , . . .R, )^ d A * = 4 + L 4 .
18. CheckIII > Rr,
if 'yes' then GOTO S 20 (without convergence),
elseIII = III + l, ps.= t s : : "(*t = 1 , 2 , . . . ,R ) , =
4 Ltr*
,
GOTO Step t4 and
Economic Load Dispatch of Thermal Generating tJnits 189
GOIO Step13.
',20.Computethe optimal loss P2, atc.
total cost 4 transmission
'.21.CheckcostlFPrcv- 'yes'GOTO Step23.
Fl S E, if
122.CheckII > R2
if 'yes' then GOTOStep23 (withoutconvergence)
elseII = II + l, F Prev= 4 GOTOStep7 andrepeat.
:23.Stop.
0.02+j0.06
N
o
+
t+
O
0.01+70.03 0.08+.10.24 o
Figure3.17 Powersystemnetwork(Example3.11).
5 5
Total real demand, Pa.= 1.65 p.u. and total reactivedemand, Z Qo,= o'40P'u'
i=l l'=l
The initial I, is obtained as
4= = 253.0147
The initial PEiQ= L, 2,3) and Qr,{t = l, 2, 3) are calculatedfollowing Step2 andStep 3 of
Algo;rithm 3.8 and are tabulatedin Table 3.21 along with loads at each bus. The load flow
solutironis obtainedusing the decoupledload flow method (see Section2.11). The convergence
0.00001 is achievedin six iterations.The voltage magnitudes,voltage angles,real and reactive
pow€rs injected at each bus are obtained and are given in Table 3.22. Loss coefficients are
calculateduging Eq. (3.126). The valuesof angles 4 0, and 0 are tabulatedin Table 3.23. For
exam-ple,Qr = tan-L(QrlPrl= -0.8224847, 0r = 6, - Qr = 0.8224847;and so on for Qz, Qt, Q+
and Qs.
.ra
I Dispatch of Thermal Generaing Units 191
Pa Qa Type
(p-u.) (p.u.)
0.00 0.00 Slack
0.20 0.lo PQ
0.45 0.15 PQ
0.40 0.05 P8
0.60 0.10 PQ
6 P 0
ad) (p.u.) (P.u.)
0.4600578 -0.4955128i
08256 0.48674M 0.0664834
03014 0.0715068 -0.023574r
95604 -0.4000043 -0.M99996
28022 -0.6000013 -o.099999E
r 0 and 0 (Example3.11)
g(rad)
Otr"Ol
-0.8224847 0.8224847
0.1357489 -0.1565744
-0.3184559 0.2781545
0.1243527 -0.1739131
0.1651480 -0.237954t
low:
I"r,-IPo,-P, = 0.45626248-03
> 0.0001
i=l d=l
or
lt7;lt
l> | #l
/.,r,
*l'\ar,)
# | =-'.'-
7.7a826000>
o:oool
llfi[aP,,
)
or
-
,l-.........--...-
,/I f**l' +(M,")2= 3.76299591
> o.ooor
I i=r
No convergencecriteria is satisfied,so to go for next iteration,
i.e. 2, the valuesare updated8s
Pr,= pr, + dPr,. (i = l, 2, ...,NG) and 4= 4* &
Po = 0.441756_ O.V27|3O
= 0.414626p.tr.
Pc,= 0.6g6Vit + 0.010359= 0,697097p.u.
Pg:= 0.521505+ 0.034212= 0.555717p.u.
1, = 253.0107+ 3.7628= 256.7735Rs/p.u.h
' The aboveprocedureis repeatedtill any one convergence
criterionis satisfied.After the
fourth iteration, one of the criteria is satisfied.
.,"d
Economic Inad Dispatch of Thermal Generating Units 193
The real generationsare given in Thble 3.24 which are used as initial values for the next
iteration. The cost and transmissionloss after the first overall iteration is
Cost = 696.1357 Ryh, PL - 0.01829998p.u. and Lp = 256.9161Rs/p.u.h
The above procedure is repeatedfollowing Algorithm 3.8 till the slack bus balance is
achieved.After three iterations, the convergenceis obtained and the final solution is given in
Tables3,,25and 3.26. The slack bus mismatchis
- Pa,-4 = 3's077338-05
< o'cool
lt, |
4 = 256.9665 Cost= 695.8972Rs/tr,Py - 0.0L739662
Rs/p.u.h, p.u.
Bus V d
(p.u.) (rad)
1.060000 0.0000000
1.071436 -0.0191460
r.069649 -0.o37t636
1.066864 -0.0466896
1.052804 -o.w06622
t94 Power System Optimization
, , Using the Lagrange multiplier method, the constrained,optimization problem given ,in
Eq. (3.135) is convertedinto unconstrainedone and is given in,Eq. (3,L37).
The necessaryconditions for the optimization problem,:given by Eq. (3.137), state that the
variables( Pr,,Qr, (; = lrZ,..,NG) and 4), *"equal to
O"riuitiueswith r"rp"r, to controVdecision
: -
zero,1.9.
=Q (f=1,2,...,NG) (3.L43a)
#=#+ho[#-')
(i = l, 2, ...,NG) (3.143b)
6Qgt
NB NG
q = I P a . + P L - I P r ,= 0
AL
(3.1a3c)
i=l i=l
L .. .-.'.---_. .--l
Economic Load Dispatch of Thermal Generati4g Units 195
NB
=; hi,Q,+
I l@u+ ai) Qi + @u- bi) Pil (3.145)
'l=l
j*i
AL
(i = 1,2, ...,NG) (3.146a)
dP,,
AL ( i = 1 , 2 , . . . ,N G ) ( 3 . 1 4 6 b )
dQr,
= _ - AL Q.ra6c)
dLo
,
L, d'P,
(i = 1,2,...,NG;i = 1,2,...,NG;' * i) (3.r47
a)
Affi
From Eq. (3.143b),the following can be obtained:
d, p, d, p, (3.r47
(r = 1,2, ..,NG, j = L, 2, ..,NG) e)
dPrdQr, aQs,aPt
81
l* --;,
196 Power System Optimization
EXAMPLE 3.12 Considerthe 5-bus systemof Example 3.1I and obtain the optimum schedule.
4lt
Using the Newton-Raphsonrnethod elaboratedin Algorithm 3.8, the real and reactive power
generationsare obtained in six iterations
The Hessianmatrix elementsare obtainedas given below:
l,2,
Ht" = += ?-ar+ Z)ra;; (i = 1, 2, 3)
dPr,
,1ii=&=A'o(aii+ai) ( r = 1 , 2 , 3 j; = 1 , 2 , 3 i a n d i * j )
rr dzl
trr;1s- = lvp(bii-bi) (f = 1, 2, 3; j = l, 2, 3; k - - 4 , 5 , 6 )
ffi
i --------._* ..-l
-d
Economic Load Dispatch of Therrnal Gencrating Units 797
lrvg
r r * = f faiz=t L } ( { o , * a i ) P i + ( b i ,i ._^ b
H , , i. ll n a
) ). - , , l . 0 ( f = 1 , 2 , 3 ; k = 7 )
rr
H * i = azt A n ( b ibi -1 ) ( r ' =1 , 2 , 3 ;j = 1 , 2 , 3 i k = 4 , 5 , 6 )
ffi=
t,
H u = f fd zi t_ \ ( a i 1 + a 1 i ) ( l = 1 , 2 , 3 ; j = I , 2 , 3 i k = 4 , 5 , 6 ; I = 4 , 5 , 6 )
NB
rr azr
H *t=ffi=Zk,*ai)Q1+@u-bi)Pi)(i=l,2,3;k:4,5,6;t=7)
, p 1 -dzr
H W = f f i azt (f=1,2,3
k=; 7)
rr p - dzr
H azt (i=1,2,3
k=; 1 tl = 4 , 5 , . . . , 7 )
m= ffi
tr azt =0.0
Hn= ^
EL",
l'xs . )
#-=ti = 4lI +ai)Q;+@1;-b1ltP)
(<o, | ti= r, z,r)
\.Fr )
rr NB fNs Nr \ Nc
= >Pr, *l I \ aa(4pi - 4e) I
+ee)+ b;;(QtP1
#- " p l-
i=l \i=r;=r ) i = r "r,
3 5
I ar,- I Pd,:P, = 0.00000002
< 0.00001
i=l i=l
P8 Qs Pd Qa
(p.u.) (p.u.) (p.u.) (p.u.)
I o.426539 -0.151419 0.00 0.00 Slack
2 0.691533 -0.116962 0.20 0.10 PQ
3 0.547509 -0.062548 0.45 0.15 PQ
4 0,0 0.0 0.40 0.05 PQ
5 0.0 0.0 0.60 0.10 P8
- Par-4 6'437302E-06
< o'ooool
It, l=
The final solution is given in Thble 3.30. The convergenceis obtainedin each iteration and
is tabulatedin Table 3.31.
*-*
Ecorromic Load Dispatch of Thermal Generating Units
lossaregivenbelow:
Theoptimumvaluesof 4, cost,andtransmission
NG NB
where
ar, bi, and c; are cost coefficients
Po, is the real power load demand at the fth bus
Pr, is the real power generation(decisionvariable)
Qa, is the reactive power load demand at the ith bus
Qr, is the reactive power generation(decisionvariable)
Ps is transmissionreal power loss
Qr is transmissionreactiygpower loss
NG is the number of generationbuses
NB is the number'of busesin the network. : .
This methrod uses the fact that under normal operating conditions, the transmission loss
is quadratic in the injected bus real powers. The general form of the loss formula obtained in
Eqs. (3.132)and (3.134)is
NBNB
p2- -
I It, i1(4P,+ QiQi\+ ba(QtPi 4Q)l (3.149a)
i=l j=l
NBNB
- (3.14eb)
Qr- I It' u(4r, + QiQil+ d;i(QiPj 48)l
i=l j=l
#
200 Power SysternOptimization
where
R,,
sii= cos(d;- dr)
ffi
R,,
bii=ffisin(d;-,t,
X
=N,ffi cos(4-6,
X' i i
d
) ri=ffi
s- : i n. / c( d ; - 6 ,
Pi + iQi= (Ps,- Pa) * j(Qs, - Qa) (i = I ,2, ..., NB)
P; and Pi are ttre real power injections at the ith and 7th buses,respectively
Q; and Q1 are the reactive power injections at the ith and 7th buses,respectively
Zi = Rij + jXil @lementsof impedancematrix).
Using the Lagrange multiplier method, the constrainedoptimization problem given by
Eq. (3.148) is convertedinto an unconstrainedone.
Lr) = F(ps,)
L(ps,,er,, +A, p4, L $ p ] *,ta
+ ' p-- -i
[i t'X ea+er \,'=, A'r')' \i=r i=l
n*I (3.1s0)
)
where 4 it the Lagrangianmultiplier.
Necessaryconditions for optimization problem statedby Eq. (3.150) are
_
a L = a F *, ,Lo ' ) * , d Q r-_ n
, ( r r r _-r)*
q Aoq 0 (i = 1'2'"'' NG) (3.151a)
aP& ttr*
Nts NG
AL
= Pa+PL-I =o (3.151c)
4 I "''
-\r NB NG
wherethe incrementaltransrnission
loss expressioniare expressed
hereas
NG
dP,
=Za;iP;+ t i = 1 ,2 , . . . , N G ) (3.r52a)
dr. Zlf"u*aii)Pi+(bi,-bilQi]
j*i
NG
q
,dP, = 2aiiQ;
+ * ai)Qi+@u-bi)Pj] t; = t, z, ...,NG) (3.tszb)
Zrr,
j+i
-
EconomicLoad Dispatchof ThennalGeneratiryS
_Utt!!!
NG
=2ciiPi* - (3.152c)
l@u + c i) Pi + (d ii di)Q il
W .i=l
j*i
NG
dQr
dQr,
- 2ci1Qi+
I
j--l
-
[(",i + ci)Qi + @u d i)P;]
(i = 1,2,...N
, G) (3.152d)
j*i
EXAMPLB 3.13 Considera 5-bus systemof Example 3.ll and obtain the optimum schedule.
Solution Ireuselements and Z,u'5elementsare sameas those of Example 3.11. The initial
'values
and load flow solutions will also come out same as in Example 3.11 and are given
in Table 3.2I and,Table 3.22, respectively.The loss coefficients a;i and b;i are tabulatedin
Tables 3.27 and 3.28. (see Exarnple 3.I2). The loss coefficientscy and dii are tabulatedbelow in
Tables3.32 and 3.33
gri= (d;- 6,
ffilcos
X,,
dii= sin(d,- 6r) (i = 1,2,...,5;i = 1,2,...,5)
ffi
Table 3.32 c-coefficientsfor reactivetransmissionloss Q1 @xample 3.13)
Ct Cp C6 C6
I
L
tr
'202 Power System Optimization
loss Q1 (Example3.13)
Table 3.33 d-coefficientsfor reactivetransmission
t, =
H* azr - du)
= )r(b1'- b') + Lo(d1i ( f= 1 , 2 , 3 ; j = I , 2 , 3 ; k=4,5,6)
ffi
Hn=
azr -1.0; ( i = 1 ,2 , , 3 )
)Pi+
aPsiaLP
dzr Pi+( ( i = 1 ,2 , 3 )
.Hn=
aPsia)"q
, azt = 1 /, - b1)
,- + ).n(ds- d1i)
H *3i = Tm Lo(b6 (i- 1,2,3; j = I,2,3)
Hi*3,j*3 =
azr = 4@U + a1) + I'oGi1+c;) ( i = 1 , 2 , 3 ; j = I,2, 3)
dQrdQs
NB
a2r = ( i = I,2,3)
Hi*3,8 = ( t " u* c i ) Q i + @ u - d j ) P j )
aesM,q Ij=l
Hlt =
dzr a2t 3)
aPsiahp
H7,i+3 =
d2t ( l = 1 ,2 , 3 )
aQs.M"
e
.*rd
Economic Inail Dispatch of Thermnl Generating Units 203
onofi=ffi (i=r'2'3)
Hsi= -Plt dzt
-H
- o'ttJ
E , i * 3#=-
azL dzr (i= 1,2,3)
u, qaQs aQ\.a;,q
Hn= *= F/s=
--oo s *- = H 'u
tB=-# = Ht.=
dto dL'n dlv pdLq
aL -- ^?n;pr, l-frr-B
lfltr ) I
+ b,+^"l[](tou + a1) p1+ (bii- bil O)l- t.o
q I
) J
I
+ LofX(tr,* c1)p1+ (d1i - dij)o)
| ,, = r,2,3)
L.Ft J
:
I
dt = LLE({"u*ai)Q1
^[mn
+@u-bi,'t,).1 r . f, , \r
ft
[(rn .'\ I
) i + @ ,-i d i , pl ; ) - t o
+ L o l| ( t ' u * c 1 e = r , 2 ,3 )
I | |,t
L\r=r ) J
-x
ioi)l
A1
u "
P
z-r
i=l \;=t j=l ) ",,
i=l
1r NB / r'rg NB ) NG
ctL -- s ea +l c;1e;pi+erei+dii(eipi-pieill-ZO*
'
dhn > > u\ri'.r wiYi
# \i=r j=t ) i=l
The modifiedgenerations
which are usedas initial valuesfor the next iteration'are given below:
Pr,= 0.4265079p.u., Pr, = 0.6916468 p.u.
P.u., Prr= 0.5479347
Qr, = -0.005769
p.u., Qrr= 0.22L1785p.u., Qrr= 0.2322112
p.u.
The convergence in six interations,
obtainedduringNewton-Raphson is givenbelow:
3 5 -
L,,,
i=I i=l
= -{.00000019< 0.00001
The values of Lp, Lq, Py and Qp for the obtained schedule are given as
L, = 255.9044,L, = 0.34058850
Pr= 0.01608950p.u., Qr = 0.04761850
p.u.
Power System Optimization
In two iterations, the ovbrall convergencevalue obtained is 4.0252518-03. The final solution
achievedis given in Tables3.34 and 3.35.
b PL - 0.01589258
= 255.8771,Lq = 0.3357378, p.u., er - 0.046g0209
p.u.,
Cost = 695.4353 Rs/h
3.16 EVALUATION
OF INCREMENTAL
TRANSMISSION
LOSS
The transmission loss can be expressedin terms of B-coefficients. Hence penalty factors or
incrementallossescan be evaluatedin terrnsof B-coefficients.The transmissionlossescan also be
expressedin terms of power flow equations.The transmissionloss in terms of power injection at
various buses is
NB NG NB
P r =I e (3.154)
i=l i=l i=l
where
Pr, is the power generatedat the ith unit
Po, is the bus power demandat the ith bus
& is the bus power at the ith bus
NB is the total number of busesin the power system network
NG is the number of generatingbusessuppryingreal power.
The incrementaltransmissionloss for the ith generatingunit is
( i = 1 , 2 , . . . ,N G )
-.{
Economic Load Dispatch of Thermal Generating Units
NB
Qi= I I V,llVillcusin(d;- 6j) - Bucos
(d,- dj)l (3.1s6)
j=l
where
Y,j = Gu + jB,i are the elementsof the bus admittancematrix
I Vi I is the voltage magnitude at ttre rth bus.
4 is the voltage angle at the ith bus.
For slack or referencebus, angle is zero or fixed one in some cases.
SubstitutingEq. (3.155) into Eq. (3.154),
NB NB
PL- I I lq ll vj llcucos(6,- di) + Bijsin(6,- 6j)l
i=l j=L
On rearrangement,
NB NB NB
Pr= I I v , l 'G , , +) I t V l l v i l f c u c o(sd ,- 6 i )+ B , i s i n (-d 6; j ) ] (3.157)
i=r i=l j=l
Equations(3.155) and (3.L57) show that the distribution of 4 and P, dependon rhe bus voltage
magnitude and angle.
(3.160)
or
(3.161)
-
^ l
!
(3'157):
Expressionsfor the elemenrsof the precedingmatrix can be obtainedfrom Eq-
NB
dPt (3.162a)
E6t
.l=l
j+i
NB
ly, l Gii+ )
-6r)]
[ Z lV j I G Uc o s ( d ; ( i = M B + 1 ,M B * 2 , - - . , N B ) ( 3 . 1 6 2 b )
j=l
j*i
= cos(4 - 6j)l
(d,- dr)+ -Bu (3.r62c)
d d , $ v,llvj ll-Gusin
33
#j+i '
- l v i ll v i l [ c u s i(nd -; 6 j) - Bu cos(6, - tj)]
# (3.r62d)
, B ; i = 2 , 3 , " ' ,N B ; i + i )
( i - 2 , 3 , . . .N
NB
- z l v ; l G i iI+ t vi llcti cos(d,- 6i)- Busin(6,- 6i )]
# j:l
j*i
(i = MB + 1; MB + 2, ...,NB) (3.r62e)
r'
( i = M B + 1 , . . - , N B/;= M B + 1 , . . ' , N B ; i * i ) (3.r62t)
= y. v, vj (6;- 6j)l
(d;- 6;) * Biisin (i = 2, 3, ',NB) (3.162e)
d 6 ' A ' I ll llcucos
P
j*i
= lvillviII-G,icos
(6;- 6 j ) - B u s i n ( 6-;d j ) ]
#
(i -- 2, 3,-.-,NB; i = 2, 3,-..,NB; i * i) (3.r62h)
NB
- 2l v , l n , I, +tvj llcu sin(6,- 6i)- Bucos(6,- dr)l;
#J- j:l
j+i
(i = MB + l, MB * 2,...,NB) (3.1621)
.^.{
Economic Inad Dispatch of Thermal Generating Units
(3.163)
=[#]
t#lt#l (3.164)
The elementsfor the abovematrix are alreadyexpressedin Eqs. (3.162a),(3.I62c), and (3.162d).
The (NB - 1) simultaneousequationscan be solved by any suitable numerical technique,viz.
Gausselimination,GaussJordon,etc.
(3.165)
FromEqs.(3.165)and(3.166),
(3.167)
(3.168a)
IL
Power System Optimization
- l:]' f*l
APr (3.l68b)
LpI Lnol
where
lal
I l =
Lp)
and
A = La2,d3, ..., Gr.rg]r
F = [F*s+b Fus*Z,..., fin"]t
Eq. (3.168b)can be rewrittenas
NB NB
APr- (3.16e)
\a,Lp,+ \F,tO,
i=2 i=MB+l
From Eq. (3. 154),
L,P2=AP1
+ X oO (3.170)
i=2
To get the changein total transmission
loss,substitureEq. (3.169)into Eq. (3.170)
NB NB
Lpy-Ift*d)A4+ (3.171)
i=2
\F,tO,
i=MB+l
The incrementaltransmissionloss is definedas
the changein transmissionloss due to a changein
generationi, keeping all other generatorsconstant,
i.e.
dP,
= dP, = 1 + d i (i -2,'.., NG)
;4 (3.r72a)
{
For slackbus
dP, dP,
=o'o (3.r72b)
tro=;4
(3.r73)
. --. . ...-..c
Economic Load Dispatch of Thermal Generating units 209
( 3 . 147)
(3.r74a)
or
dPt_ _ r _ l
q-t-T Q.r74b)
Incrementaloperatingcost is defined by Eq. (3.19),
AF = Z a ; P r , +
q 'bi
( f,,,\
z a i P r ,b+i - L l t - # | f3 . r 7 s )
( dPr,
)
From Eqs. (3.174a)and (3.175),
?n,;Po,
6i
* bi= +
Li
or
'\
h ( l' ,
' s ' = [ e - o ' ) -rlo- , (r= I ,2, ...,NG) (3.176a)
or
P-si
-= -- -L (i = 1,2, "', NG) (3.176b)
?-,L,- ?4
The total transmissionloss is expressedas a function of initial power loss Pf and,the changein
total transmissionloss A,Pp,i.e.
Pr= Pf + A'P-
Since power changesin load busesare zero
P7--Pf y
L/+ rt- +) l^+ - L- P-el
^
\. L,)Lh, L, 2o, t'J
fr
L - j
210 Power SysqemOptimization
pL-pr*}^W (3
17E)
+G+.";)]
#) }l(++p,e)
The power balance equation is
NG NB
NG NB
When Pt= 0, then Z,r, = Z Po, (3.t79a)
j=l i=l
yl:r +l=xpa+pf
*},ffi hl }l(+*p:)
fr LZa,L," +(+.4)]
On simplification and cancelling the equal terms,
NG ^ l-NB NG NG - / t \
P a i - r "+' ?P : . I+ l * * r ! , )
Zm=L:
ConsideringEq. (3.L79a),
p!*;+w.",:)
)"- (3.180)
Given the numberof busesNB as 4, the number of lines NL as 4, and the number of voltage-
conffolledbusesMB as 2, find the economicgenerationschedule.
0.08+70.20
c!
o
+
.+
a
o
0.M+70.14
L- 608.2354
[*Pa+**,)/*
P r ,= ( L - b i l ( 2 x a ) (l= I,2,3)
Bus P8 Pa Qa Typt
(p.u.) (p.u.) (p.u.)
1 2.572354 1.85 Slack
2 2.192353 1.45 PV
3 2.235295 2.10 PV
4 0.0 1.60 0.80 PQ
Pt = 1.079061E-01
p.u.
Economic Load Dispatch of Thermal Generating Units 213
Ps
(p.u.)
Pd
(p.u.)
Qa
(p.u.)
P o
(p.u.) (p.u.)
I 2.572354 1.85 0.8302581 0.2982M9
2 2.192353 r.45 0.7423536 -0.3050831
3 2.235295 2.r0 0.1352946 1.0886370
4 0.0 1.60 0.80 -1.60 -0.7999994
lvl d
(p.u.) (raO
I t.02 0.0
2 1.04 0.03898002
3 1.06 -0.03029456
4 0.94239460 -0.w235717
-8.0574 0.0
14.75836 - 6.467n62
-6.7ffi96r 13.76972
- 1.56199 3.849879 I
dP, = . ap, =
o'o (slack,bus),
OP,
0.089306, = 0.079922
E E ?Pg
FromEq.(3.174a),
we get
L--,_:.
214 . Power System Qptimization
)"= = 646.711900
+ l
k h,4
From Eq. (3.176a),generationsare computedas
Pr, = 2.957119
p.u.
-P9 2 = 1.999563
p.u.
'
-P8 3 = ( n - k, ) r Lrzszro
p.u.
tu )tr=
The overall computed cost is
F - 3571.481'
Rs/h
a 6, . 6;nu*
d,lntn (i = 2,3,..., NB) (3.l8lf)
(d) Functional constraint which is a function of control variables.
- limits on reactive power
NB
Qi(v,d) = v,Z vj(Gijsin(d,- Aj)- Br cos(4 - dj)) (3.l8li)
j=l
I
hvhere
i NG is the numberof generatorbuses.
NB is the number of buses
NV is the number of voltage controlled buses
Pi is the active power injection into bus i
Qt is the reactive power injection into bus i
Pa,is the active load on bus f
Qa,is the reactive load on bus I
Pr, is the active generation on bus i
Or, is the reactive generation on bus i
Vi is the magnitude of voltage at bus i
4 is the voltage phase angle at bus i
Yij = G,i + jBii @re the elementsof admittancematrix).
The constraint minimization probleq can be transformed into an unconstrainedone by
lugmenting the load flow constraints into the objective function. The additional variables
me known as Lagrange multiplier functions or incremental cost functions in power system
lerminology.The cost function becomes:
; N B N B
L(Ps,v,A = F(Ps)+I )'r, lPi(y,d)- Pr,+ Po)+ Z ^0,[Qi(vd)- er +ed,l (3.182)
i=l i=NV+l
r
F 216 Power System Optimization
# =Xl^,,#]*,=H.,I^,,H]
(i=2,3,
NB) (3.183b)
# =X[^,,#]*,=H.,#l
L^,,
(r=NV
+,,NV
i 2,,NB)(3.183c,[
- Pt(%O - pr, po, (r= 1,2,...,N8)
! (3.183d)
#
=Qi(V,D-Qs,+Qa ( f= N V + l , N V * 2 , * . , N 8 )
T .li
(3.183e)
Any small variation in control variablesabout their initial values is
obtainedby formingthe
total
#ft;av,+Idha,pi
,=H.,
+ y, ;a,*- s'vq'
LL-=-- L (l=1,2,...,NG) (3.184a)
drrM"n dPr,
7=fi*r
X#b '.r,.X#,#46,
+,=H.,
#; ^vi+I ^hpi
"h
*,=$.,#;qLho, (3.184b)
$ dzt NB
) s dzt
aviaPsj LL,
fr j=2 aviaLej
+ $ --a'r
LLo' =-- dt (i = NV + 1, NV * 2,..,NB)
.,=fi*rav1i- dvt (3.184c)
$ a ' r &w.Z
n * u
fi dtr,a\
NB ^"
-r=fr*,
+ ! d'L ^^ dr (i = l, 2, ...,NB) (3.184d)
Walni= dL,,
Economic Load Dispatch of rhermal Generating units 217
NG
T
j=l
a z t - Zna t
(f = 1,2,..,NG) (3.185a)
q
(3.l86a)
dzr
ad,ayk
(3.186b)
azr dPi
(3.186c)
e5;
zlE OPtimization
Power SYstem
t
(3.186d)
...,NB; j = NV + 1, NV + 2, .-.,NB)
obtained by differentiating
Second-order partial derivatives required iot Eq- (3-183c) are
Eq.(3.183c) with respectto control variables,i'e'
&=v(l=NV+1,NV*2,...,NB;/=NV+1,NV*2,...,NB)(3.18
by differentiating;
Second-order partial derivatives required for Eq. (3.183d) are obtained
Eq. (3.1S3d) with respectto control variables,i'e'
( i = 1 ,2 , . . . ,N B ; i = 2 , 3 , . ' . ,N B ) (3.188a,)
azn = 0 ( i = 1 ,2 , . . . ,N B ; j = t , 2 , . . . ,N B ) (3.188c)
7Lo,r^o,
. , : l
NB
+' L\ it
AL
d6t
,t=NV+l
(3.1e0b)
NB
AL
+ t
'
'Lrl dvi
./=NV+l
(3.190c)
AL ( r = 1 , 2 , . . . ,N B ) ( 3 . 1 9 0 d )
dL,,
z#A6,+,=H-,ffi^vj=#;(i=NV+1,NV*2,...,NB)(3
3.18.1 Limits on Variables
Kuhn-Tuckerconditions
inequality constraints are assumedto be
The Kuhn-Tucker approach can be adopted and thus'
t: its lower or upper limit
inactive during the,initial solution.process.Any variable gun^9r,1red.
-trre can be deactlvatedso
values. Then, partial derivative terms with Tqpect to^the fixld variables
of rest of variables' A variable should be
that these terms do not participate in the evaluution
correctionof all other variables'
moved to its relevant limit in coordinationwith the simultaneous'
if x; = x,ltin
(3.191b)
<0
>0 if xi=xl*
Power System Optimization
Penaltyfunction method
Penalty functions ideally fulfit the requirementsfor inequality constraintsenforcement.
The
penalty function is describedbelow which is quadratic.
f r R r t
1-t \-r J.
-t) (3.r92a)
where
* " ' =I | &
yi is the target value
y; is the current value
S; is the weighting facror.
The weighting factor is automaticallycontrolled to give the appropriateamount of hardness
of enforcement.The first and secondderivativesof a; ari:
s '( y i - y i ) (3.rezb)
#=
(3.r92c)
The proper value of S, is automaticallycontrolled.For the larger values
of ^t, the target limit acts
as the rigid4rard limit and for small value of .Si,the target limit acts as
the soft limit.
The detailed algorithm for optimal load flow is outlined below.
,
Algorithm 3.10: optimal power FIow Based on Newton Method
1' Read dataai, bi, andc; (i = 1,2,..", NG), Ioad on eachbus,
line data for the power system
network.
2. Obtain Irsususing the I-bus algorithm.
3' Calculatethe initial valuesof Pr,(i = 1,2,...;NG) and2,by assuming
-(3.2a) that pt=g. Then the
problem can be statedby Eqs. and,(3.2b) and the solution can be obtaineddirectly
usingEg.(3.10)andEq. (3.9).Initializeall Io,= A, = 1,2,...,
1i NB), 2q- = 0 (j = NV + l,
N V + 2 , . . . ,N B ) , V r = | p . u . ( t= 2 , 3 , . . . , N B ) a n d d ; = 0 ( j = 2 , 3 , . . . ,
NB).
4' Calculatethe Jacobianand Hessianmatrix elements
from Eqs.(3.1g3)to (3.1g9)
AL
LPg
a"'
dt
Ad e-
AL
LHI aAe
AV
q
AL
Llq fr
dt
dLu
The Gausseliminationmethodis usedto find LP*
A6,A4, AV, andLLn.
Economic Load Dispatch of rhermal Generating (Jnits ?:21
5. Chpck convergence
fg NB NB NB r*
lI (^4,)2 +l{uo1:t+ I ( A v , ) 2\+t t t , r f l ( s
+f rm,)z FI
7. Check the limits, if any limit of a variableis violated,then imposeor removepower florv
equationor a penalty for inequality.Add or remove derivativesfor penalty or equation
changeand GOTO Step 4 to updatethe solutlon.
8. Calculatethe total cost.
9. Stop.
0:02+70.08
o.o2+jo91
L Z
222 Power SystemOptimization
)'=[*Pa,*E*)/Er')
= 291.1111
-p8"d, - * (i=1,2)
Zxa,
Ps, Qs,P, and O are tabulatedin Table 3.42 which are obtained with the initial values
tabulatedin Thble\3.41.
Thble3.42 Othervaluescalculated
with initial values(Example3.15)
P, (n.u.) 0r (n.u.) Pa (p.u.) 0a Q.u.) P (p.u.) 2 (p.u.)
I 0.7592590 0.546t460 o.20 0.10 0.1223528 0.M6r460
2 0.9407406 0.4461460 0.00 0.00 0.1223528 0.4461460
3 0.0 0.0 1.50 0.06 -0.2352939 -0.981176/.
Cost = 653.9258Rs/h
a
With the above-mentionedvalues,the following equation is solved to initiate the iterative
processto implement the Newton-Raphsonmethod.
aL_- ^ dP,*,
= + ^o';y.* 7"'dQt= -3561
'82e
fr h At;
wherethe partial derivativeof P1 with respectto d2is
#doz = VzL
' - 6 + Bzicos(62
v;(-Gzisin(d2 - dj)) = 24.960
)
F,
j*2
?P3
+A*aQ3
= Z ho,F,f aq
#
where the partial derivative of Pt with respect to 63 is
fp
=VrVE(Gn sin (d1- ds) - Bn cos (d1 - dJ) =-12.23529
a6,-
Partial derivative of P2 with respect to Q:
a n= J -
V3Lvi?Gti s i n ( d 3- 6 j ) + h i c o s ( d 3- 6 j ) ) = 2 4 . 4 7 0 5 9
E- -*r j=l
j*3
Partial'derivative
of 0i'with respect
to d3:
E P . n '
=2VtGt *I V i ( q i c o s ( d 3- 6 i ) + B t i s i n ( d 3 ' - d i ) ) = 5 i / 7 0 5 9
;%
i=\
Partial derivative of p, with respectto V3:
3
dQl
- -2vtBt
.
.I vi(Gzisin(63 6c j)-
r
4i cos(63-6)) = 22.,508240
rr /? c \\ .r/r
M
j*t
dr - Pt - Ps * Par=o'4369062
q
AL - '
P2' Ps'* Paz= o'8183877 :':
q
a L A - - .
= QE-,Qr, * Qo,- 0.9211764
A
q
Hessianmatrix elementsare obtainedfrom Eqs. (3.185-3.189).
Economic had Dispatch of Therrnal'Generating Units ?;25
4 = 2 . o a tr2o.o;
= =oo
AP; #
azt = o.o;+ = Z.oaz=
t5o.o
aPsr
aPsz dPi,
Second partial derivatives of L with respectto )"u:
azt
wl
=;^',W*u'# = 2742.609
d, P, - -VtVr(Grrcos(di -
6r)+ Bnsin (dt . 6r)) = 1.181176
a67
3
d, P,
= vrL viFG4 sos(62- di) - Bzisin(62- 6;)) = 6.240000
asi j=l
j*2
d, P, = -vtvz(Gncos (63-
6) + Bn sin (d3 - 6r)) = 3.058824
a63
a'(h (6r - 6J + 832cos (6r - 6r )) - 12-23529
as3.
Second partial with respectto 62 and 63:
A,Q,
= -890.4575
ad2a63
where
= 0.0
d'g:-
= \v2(G32sin (d'3- 6r) - Bn cos (d_j- d, = -t2.23529
at?adr ))
= -
#--r t rn,(-cysin (d3 6i)* Bticos
(d'3- dr)) = 24.4705e
j=\
second partial derivativesof L with respectto
d,2and v3:
.l- P
*3* = vzFG*
' cos(63- 6r) - Bt sin(63- drD = -2.94rr77
a62a% z\
Secondpartial derivativesof L with respectto 4 and V3:
j?'-
= -5 -di)) = 6.117647
Btisin(63
ad3ay3 Lr,(o",cos(,83 j)+
i=l
Secondparrial derivativesof L with respectto d and \:
# =oo
=p.^,,#*',W
r
i;
i:
x28 Fower System Optimization
where
#
d,e,
= o.o;
W=
o.o;
#
- 2.0G33=
n .7647r
dzt a4 =-3.05gg24;
: = a2t ,oP', dzt dP,
= 5'64705e
w=ait dw=')fr=_3.0588240; ffi= ft
Secondpartialtderivativesof Z with respectto V3 and Ln:
=22'50824
&=W
Using the Gausselimination method,the changein variablesis obtainedand the updated
values of variablesare shown in Thbles3.43 and 3.M.
The final solution after four iterationsis given in Thbles 3.45 and 3.46.
Bus P8
(p.u.)
QB
(p.u.)
Pd
(p.u.)
Qa
(p.u.)
P o
(p.u.) (p.u.)
1 0.7767037 0.1521580 o.20 0.10 0.5767038 0.0521580
2 0.9452465 -0.0327489 0.00 0.00 0.9452462 -0.0327489
3 0.0 0.0 1.50 0.06 - 1.5000000 -0.0s99991
Economic Load Dispatch of Thermal Generating Units
Cost = 660.3356Rs/h
NB
+f dzr Ad; = AL (i=1,2,...,N8) (3.193c)
j=2
),p.a6j dL o,
dzt AL
( i = N V + 1, NV * 2, ...,NB) (3.193e)
LVi =
L,rtavj dL n,
The above equationscan be upled and can be solved separately.In the matrix notation these
can be representedas
Y ,ru
I
vao (3.r94a)
pPs
Y ^ru
ThesizeoftheHessianmatfixis(NG+2xNB-l)x(NG+2xNB-l)andthesizeofJacobi
is (NG + 2 x NB - l) x l.
lI:. =F;;l
I:;,ltil,l (3.le4b)
n a l t r i ixs [ 2 x ( N B - N V -
T h e s i z e o f t h e H e s s i am l)] x[2 x (NB -NV- l)] andthesizeof
Jacobianis [2 x (NB NV - - I )] x l. The elementsof Hessianand Jacobianmatriceshas already
beendescribedin Eqs. (3.1F3)-(3.189).
The detailed algorithrfr is 6utlined below.
I n i t i a l i zael l ) v o .= ) v ( i = 1 , 2 , . . . , N B ) ,r , q -= 0 ( i = N V + l , N V * 2 , . . . , N B ) , V i = | p.u.
( i - - 2 , 3 , . . . ,N B ) a n d 4 = 0 ( i = 2 , 3 , . . . , N B ) .
4. Calculate the Jacobipnand Hessianmatrix elementsfrom Eqs. (3.183a) to (3.189) and
solve Eqs. (3.194a)qnd (3.194b).
AL
4
dr and [H2)
a6
AL
il P
Gauss elimination m$thod is used to solve these simultaneous equations separately to find
LPs, L6, L2,,,A% and LLo.
5. Checkconvergence
NB
(s and l - N B (LV,)z+ t
1g
>
ll_i=NV+l i=NV+l I
and optinnalityconditions.If convergencecondition is not satisfiedthen GOTO Step 6 else
GOTO Srep 8. i
6. Modify the variablesl
7. Check the limits, if limit of a variable is violated then impose or remove power.flow
equationor a penalt for inequality. Add or remove derivatives for penalty or equation
change and GOTO 4 to updatethe solution.
8 . Calculate the total
9. Stop.
EXAMPLE 3.16 Consi the 3-bus systemof Figure 3.19 (see Example 3.15). The series
impedance
of eachline is in Table 3.40. The systemhas two generators.The operatingcost
characteristics of two are given below.
Bus V d hP hq
(p.u.) Gad) @Vp.u.h) (RVp.u.h)
1 1.04 0.0 292.2733
2 1.04 0.00985892 291.0713
3 1.O+0361 -0.04675326 295.4336 -6.128246
-.2
Power System O,
4(Ps)= aiP?,+b,Pr,+c,
If the systemreal power loss is be minimized,the objectivefunctionis
F = P(V, 6)
In this casethe net injected real wers are fixed, the minimization of the real injected power P1
at the slack bus is equivalent minimization of total system loss. This is known as optimal
reactive power flow problem, to the load flow equations:
(a) Real power balancein the for each PV bus
,^8
Economic Load Dispatch of Thermal Generating Units 233
ion on bus i
ion on bus i
(3.1e6)
(3.r97a)
{} } ror slackbus
LorJ
l'p.l
-'
y = { I for PV bus
l.Y'J
=[l (3.1e7b)
fp.I
PQbus
E,lfor
In the above formulation, I objective function must include the slack bus power. The vector of
independent variablesy can be partitioned into two parts-a vector u of control variables which
are to be varied to achie optimum value of objective function and a vector p of fixed or
disturbance or uncontrol partrmeters.
Control parameters be voltage magnitudeson PV buses,or Pr, at buses with control-
lable power. Slack bus vol and'regulating'transformertap-settingrnay be employed as additional
control variables. Qr, may be used as control variables on buses with reactive power control
[Singh, 1993].
The optinrLization em can be restated as
min F(x, u) (3.198a)
a L=
_ E F . [ & l t1 .- =''',
E; E'LEJ 'v o (3.resd)
#=#.[#]']'=o (3.198e)
y = g(x,u,p)= o (3.1e80
Equations (3.198d),(3.198e)a (3.198f) are nonlinear algebraic equationsand can only be
solvediteratively.A simpleyet efficient iteration scheme,that can be employed, is the steepest
descentmethod (also called gradient method). The basic technique is to adjust the control
vector u, so as to move from feasible solution point (a set of values of x which satisfies
Eq. (3.198f) for given u and p; it indeed is the load flow solution) in the direction of steepest
descent (negative gradient) to a new feasible solution point with lower value of objective
function.
The computational for the gradient method with relevant details is given below:
1. Make an initial guess z, the control variables.
2. Find a feasible load fl solution from Eq. (z.ltla and b) by the Newton-Raphson
iterative rnethod. The improves the solution .r as follows:
*r+l={+L,x
where Ax is obtained by solving the set of linear equationsgiven below:
Lx= -s(t'v)
[#,"'r,]
Lx=-f+ (",
' y)l
- s(x',y)
Ldx J
The end result of Step 2 is a feasiblesolution of x and the Jacobianmatrix.
3. Solve Eq. (3.198d)for )"
(3.199a)
YL= (3.1eeb)
o-od
Economic Load Dispatch of rhermal Generating Units z3S
In the algorithm, the choice of is very critical. Too small a value of a guarantees the convergence,
but slows down the rate of ; too high a value causesoscillationsaroundthe minimum.
Severalmethodsareavailable optimumchoiceof stepsize.
Prftn
sPr,<P#"*
These inequality constraintso control variables can be easily handled. If the correction L,u; in
Eq. (3.199c)causesui to ex one of the limits, ui is set equal to the corresponding limit, i.e.
(3.200a)
After a control variablesreac any of the limits, its componentin the gradientshould continue
to be computedin later i s, as the variablemay come within limits at some later stage.
In accordancewith the hn-Tucker theorem, the necessaryconditions for minimization of L
under constraint are:
=o if Jzlnin
<ui< uf*
ft
if ui= uln (3.200b)
#=o
i f u i= u l *
#.0
Therefore, now, in Step 5 of computationalalgorithm, the gradient vector has to satisfy the
optimality condition (3.200b).
The penalty method ls for augmentationof the objectivefunction so that the new
objective function becomes
F =F(x,u)*2wi (3.201)
j
where the penalry W1is ucedfor eachviolatedinequalityconstraint.
A suitablepenalty ion is defined as
(3.202)
where y is a real positive n ber which controlsthe degreeof penalty and is called
the penalty
factor.
The necessaryconditi ns (3.198d)and (3.198e)would now be modified as given
below,
while the condition(3.198 i.e. load flow equation,remainsunchanged.
(3.203a)
(3.203b)
The vector dW/dx obtained om Eq. (3.202)would contain only one non-zeroterm
corresponding
to the dependent variable xii while d$du - 0, as the penalty functions
on dependentvariables
are independentof the I variables.
By choosing a higher alue of
[, the penalty function can be made steeper so that the
solution lies closer to the ri ;id limits; the convergence,however, will become poorer.
A good
schemeis to start with a low value of and to increase it during the optimi zation process,if the
A
solution exceedsa certain to limit.
EXAMPLE 3,17 Considert he 5-bus system of Figure 3.17 (see Example 3.11).The
series
impedanceof each line is ei in Thble 3.20. The system has three generators.Find the economic
generationschedule.The ing cost characteristicsof the three generatorsare given below:
z = 75Ps2,
+ 150Pr,+ 120.0 Rsftr
Solution Considering 1 as slack bus, buses2 and 3 as pv buses and4 and 5 busesare
consideredas Pp buses.For samplesystem,the optimization problem is statedas
3
Minimize P- \ f",r!, + b,Pr,+c;) Rsftr (3.204a)
i=l
subject to the load flow equ ons
5
Q{V,o - - - cos(6i - dr))
Z ,,(cu rin(d; 6r) ^Bu (3.204e)
i=l
J 5 5
L o , ( 4 ( v-,P
d )r ,* Pa). > l"u,(ei(y,
d)- er,* eo,)
S't | .t
L(x, u, p) = L l",Pi, + b,Pr,+
i=1
).> i=4
i=2
(3.205)
The necessary
conditions minimization are
= Z a ; P r , + b , + ) " 0( i .= 2 , 3 ) (3.206a)
#
# = (hfs,*a,)#
I W
( i=2 , 3 , 4
s ), (3.206b)
aL _ 14_ f, , , ,, d4
(bfs,+b)#, (i=4,5) (3.206c)
fr= Y
p i ( v , 6p) -r . * ' d ,- o (i--2,3,4,5) (3.206d)
#=
AL =
=_
dL
Qi(V,d)_ Q, + =Q (i=4,5) (3.206e)
q,
V,'*t= V{ + LV, (i = 4, 5)
Initial data to impl ent the Newton-Raphsonmethod is given in Table 3.51.
Bus Typt P8 Pa Qa V 6 L P Aq
(p.u (p.u.) (p.u.) (p.u.) (rad) (Rs/p.u.h) (RVp.u.h)
Bus p
" 8 Qs Pa Q a P o
(p.u.) (p-u.) (p.u.) (p.u.) (p.u.) (p.u.)
The voltage and the Itage angle at each bus are given below:
p.u.
V t = 1.0600 4 = 0'0 rad
v 2 - 1.0593p.u. 6z= -0.01739534
rad
l=-l.ffiirlT',:l
v 3 - 1.0525p.u.
p.
V + = 1.048837
v s - 1.034887
p.
Equations(3.206b) d (3.206c) can be rewritten in matrix form as
I
Economic Inad Dispatch of Thermal Generating Units
From the above eq the values of 4 ^d A.oare obtainedby implementing the Gauss
elimination method.
L r, = 258'oo3, L o, = 261.2L59, L oo= 262'8730, L o, = 267'5160
Lno = 0j1273 , Lo, - 1.167928
llvPsll - 9.604976
> 0.r
The final solution as ined after eight iterationsis given in Thble 3.54.
- Za2Prrt bz *
L p, = 5.9l5l l7E_02
dP,,
AL =
2qzPe,* bz * Lp, = -5.g525g0E_02
dP,,
l l v P sI = 8.3211368-02< 0.1
The total opeidting cost at is schedule
is = 695.5009Rs/tr
Economic Inad Dispatch of Thennal Generating Units ut
REFERENCES
Books
Arrillaga, J. and C.P.Arnold, 'omputerAnalysis of Power Systems,John Wiley & Sons, Singapore,
1990.
Elgerd, O.I., Electric Energy Theory:An Intrcduction,Znded.,ThtaMcGrawHill, 1983.
El-Hawary,
M.E.,andG.S. tensen,Optirnal Economic Operation of Power Systems,Academic
Press,NewYork, L979.
Gross,C.A., PowerSystem , Wiley, NewYork, 1979.
Kirchamayer,L.K., Economi Operation of Power Systems,Wiley EasternLtd., New Delhi, 1958.
Kusic, G.L., ComputerAided Power SystemsAnalysis, Prentice-Hall of India, New Delhi, 1986.
Mahalanabis,A.K,, D.P ari, and SJ. Ahson, Computer Aided Power System Analysis and
Control, Thta McGraw-Hi New Delhi, L99I.
Nagrath, I.J. and D.P. Modern Power SystemAnalysis, Thta McGraw-Hill, New Delhi,
1989.
Nagrath, I.J. and D.P. Power SystemErigineertng,Tata.McGraw-Hill, New Delhi, 1994.
Singh, L.P.,AdvancedPower tem Analysis and Dynamics,Znd ed.,Wiley EasternLimited, New
Delhi, L993.
Stagg,G.W.andA.H. Ei- ComputerMethods in Power Systen;tsAnalysis,McGraw-Hill, New
Delhi, 1968.
Stevenson,W.D., Elements Power SystemAnalysis, 4th ed., McGraw-Hill, New York, \982.
Wood,A.J. and ts. lV g, Power Generation, Operation and Control, John V/iley, New
York, 1984,
Papers
Alguacil, N. and A. J., C , Multiperiod optimal power flow using benders decomposition,
IEEE Transactions on r systems,PWRS-15(1),pp. L96-201,February2000.
Alvarado, EL., Penalty fac from Newton'srnethod,IEEE Trans.,YoL PAS-97,No. 6, 2031-
2037, 1979.
Aoki,.K. and T. Satoh,New gorithms for classiceconomicload dispatch, IEEE Trans. on Power
Apparatus and Systems, 103, No. 6, pp. 1423-L431,1984.
Burchett,R.C., H.H. Happ, Verath, and K.A. Wiqgau,Developmentsin optimal power flow,
IEEE Trans.,Vol. PAS-101No. 5, pp. 406414, t982.
Carpentier,J.L., Optimal flow, Electric Power & Energy Systems,Vol. 1, No. 1, pp. 13-15,
1979.
Carpentier, J.L., Optimal er flow: uses, methods and developments,Proceedingsof IFAC
conference,RI Brazil, pp. tr-25, 1985.
Chen,L., S., Matoba,H. In and T., Okabe, Surrogateconstraint method for optimal power flow,
IEEE Transactionson Pow r system,r,PWRS-13(3),pp. 1084-1089,August 1998.
Power SYstem
B.H. andS.
ChowdhurY, A reviewof recentadvancesin economicdispatch,IEEE Trans.
Vol. PW
on PowerSYstems, , No. 4, pp. L248-1259,1990.
Chun-LungChen and Nanmi Chen, Direct searchmethod for solving economicdispatch
problem considering capacity constraints, IEEE Transactions on Power Systems,
PWRS-I6(4), pp. 76L7 69, 2001.
da Costa,G.R.M., C.E.U., and A.M., de Souza, Comparative studies of optimization
methods for the optimal flow problem, Electric Power SystemsResearch,Vot 56,
pp. 249-254, 2000.
Da Costo,V. M., N. Martins, J.L.R., Pereira,Developmentsin the Newton-Raphsonpower flow
formulation based on injections, IEEE Transactionson Power systems,PWRS'14(4),
pp. 1320-1326, 1999.
Da Costa,V.M., N., Martins, nd J.L.R., Pereira, An augmentedNewton-Raphsonpower flow
formulation basedon culren injections, Int. J. Electrical Power and Energy .Sysrens,VoL 23,
pp.305-312,2001.
Dommel, H.W. and W.F. Ti Optimal Power Flow Soutions,IEEE Trans.,Vol. PAS'87, No.10,
pp. 1866-1876,1968.
Dopazo, J.F., O.A. Klitin, G. Stagg, and M. Watson, An optimization technique for real and
reactive power allocation, P edingsof IEEE, YoL 55, No. 11, pp. 1877-1885,1967.
Early, E.D., RE. Watson, and G.L. Smith, A general transmis-sionloss equation, NEE Trans,
Vol. PAS-74,pp. 510-520, 5.
Ekwue,A.O. and J.F. , Comparisonof load flow solution methods,Electric Power
SystemsResearch,Vol. 22, o. 3, pp.2L3-222,199L.
Famideh-Vojdani, A,R. and F. Galiana, Economic dispatch with generation constraints, IEEE
Traw. on Automatic Vo[ 25, pp.213-217,1980.
Fang,R.S., and A.K., David, imal dispatch under transmissioncontracts,IEEE Transactionson
Power Systems,PWRS-I pp. 732-737,May 1999.
Happ,H.H., Analysisof netw IEEE Trans,Yol. PAS'82,
with complexauto-transformers,
pp. 75-81, t963.
Happ, H.H., Analysis of ne with complex auto-tansformers, II: Relations between all open-
path and open-path-closed- impedance matrices, IEEE Trans, YoL PAS-82, pp. 958-965,
L9;63.
Happ, H.H., Analysis of with complex auto-transformers, III: Invariant reduction for
purposeof loss formula, I Trans,Vol. PAS-83,pp. 707:714,1964.
Happ, H.H., Optimal power IEEE Trans.,Vol. PAS-93,No. 3, pp. 820-830,1974.
Happ, H.H., Optimal power tch-A comprehensivesurvey, IEEE Trans.,Vol. PAS'96, No. 3,
pp. 841-854, 1977.
Happ, H.H., J.F.Hohenstein,L. Kircham ayer,and G.W. Stagg,Direct calculation of transmission
loss formula-Il, IEEE Trans. Vol. PAS-83,pp. 702:707,l9&.
Hill, E.F. and W.D. Stev , A new methodof determiningloss coefficients,IEEE Trans.
Vol. PAS-87,No. 7, pp. I 1553,July,1968.
Economic IAad Dispatch of Thermal Generating Units 243
'metric
Housoi, E;C. and G.D: Irisarri A sparsevariable optimisation method applied to the
solutionof powersystem IEEE Trans.,Vol. PAS'101, No. l, Pp. 195-202, fan, 1982'
Huneault,M. and F.D. Galiana,A survey of the optimal power flow literature, IEEE Trans. on
PowerSysrens,Vol. PW L99L
No.2, pp.762-770,
Part II
IEEE working grouP, DescriPti and bibliography of major economic security function,
and
and lll-BibliograPhy (195 1972 and lg73-lg7g), IEEE Trans. on Power Apparatus
Vol. 100, No. 1, PP. 2 L 5 - 2 3 5 , 1 9 8 1 .
Sysreru.t,
reference
Kron, G., TensorialanalYsisof integratedtransmissionsystems,Part-I: The six basic
frame, AIEE Trans.,Vol. 70, . 1239-1248,195i.
turn ratios,
Kron, G., TensorialanalYsisof integratedtransmissionsystems,Part-II: Off-nominal
AIEE Trans.Vol. 71, pp. stz, L952.
IEEE
Liang, Z. and J.D. Glover, Im ved cost functions for economic dispatch computations,
Trans on Power System.s, Vol. PWRS-6,No. 2, PP.82I-829, L99L
real-time line flow
Lin, C.8., S.T. Chen, and C.L. Huang, A two-step sensitivity approachfor
calculation, Elecffic Power stemsResearch,Vol. 2L, pp. 6349, 1991'
ordered
Meyer, W.S. and V.D. Al , Improved loss formula corirputation by optimally
eliminationtechniques,IEEE Trans,Vol. PAS'92,PP.6249, l97I-
Mohamed-Nor, K. and A.H.A. Rashid, Efficient economic dispatch algorithm for thermal unit
commitment, IEE Proceedin s C Vol. 138, No, 3, PP.2I3-2L7, l99t'
and P. Ristaanovic,
, M.S. Stott, D. Sun, A. Papalexopoulos,
Momoh J.A., R.J. Koessler,B.
Challenges to Optimal er Flow, IEEE Transactions on Power Systems, PWRS-Ia(L),
pp. 44M47, February1997
Monticelli, A. and W.H.E. Liu, taptivemovementpenalty method for the Newton optimal power
flow, IEEE Trans, on Power stems,Vol. PWRS'7, No. 1, pp.33F342, t992'
dispatch algorithm
Nanda, J., L. Hari, M.L. Koth , and J. Henry Extremely fast economic load
Vol' 139, No' 1, pp' 3946'
through modified coordinati equations, IEE Proceedings,Part C
t992.
evaluations,IEEE
}[g, W.Y., Generalizedgenerati distribution factors for power system security
Trans. on Powe4 APParatus Sysfems,Vol. 100, No. 3, pp. 1001-1005,1981'
Electric Power
Palanichamy,C. and K. Srikris na, A method for short-term generationredispatch,
SystemsResearch,Vol. L7, 129-138,1989.
Palanichamy,C. and K. Srikri Simplealgorithmfor economicpowerdispatch,ElectricPower
SystemsResearch,Vol. 21, 147-153,199r.
S.C. Tripathy,Modificaions to Newton RaPhsonload flow for ill-
Prasad,C. Durga,A.K. Jana, 'nt.
SYsrerns,Vol. 12, No. 3,
conditionedpower sYstems, Joumal of Electrical Power & Energy
pp. 1,92-L96,1990.
Ramaraj, N., R, Rajaram,and A new analytical approachto optimize a generation
Parthasarathy,
schedule,Electric Power tems Research,Vol. 11, pp. L47-I52, 1986'
to power system
Sasson,A.M. and H.M. Merri l, Some applications of optimization techniques
problems,Proceedingsof I , Vol. 62, No. 7, pP.959-972, 1974.
Power System O,
Shoults,R.R., W.M. G and s. Helmick, An efficient method for computing loss formula
coefficientsbasedupon method of least squares,IEEE Trans.,vol. pAS-9g, pp. 2144 zl5z,
1979.
Sun,D.I., B. r\shley,B. rewer,A. Hughes,and W.F.Tinney,Optimalpowerflow
by Newton
approach,IIiEE Trans., PAS-103,No. 10,pp. 286+2880,1984.
Talukdar,S.N. and F.F. W , computer-aideddispatchfor electric power systems,
Ptoceedings of
IEEE, Yol. 69, No. 10, p . 1 2 1 2 - 1 2 3119, 8 1 .
Tao Guo, Mark I Henw d, and Mariekevan Oaijen, An algorithm for combined heat power
economic dispatch, Transactions on Power Systems,PWRS-11(4), pp. t77}-t794,
November 1996.
Wang, C. and S.M. S hopur, Optimal generationscheduling with ramping costs, IEEE
Transactions on Power stems,PWRS-10(1),pp. 60-67, February1995.
Whei-Min Lin, Fu-Sheng ng and Ming-Tong Tsay, An improved Tabu search for economic
dispatchwith multiple rnrma,IEEE Transactionson power systems,pwRs-l7(l),
pp. l0g-
112, February2002.
Wong, K.P. and K. Doan, recursiveeconomic dispatch algorithm for assessingthe costs of
thermal generatorsched , IEEE Trans.on Power system.l,vol. pwRs-7, No. 2,, pp. s77-5g3,
1992.
Yung-ChungChang,Wei-T Yang, and chun chang Liu, A new method for calculatingloss
coefficients,IEEE ctions on Powersystems,PWRS-9(3),pp. 1665-1671,Augusr lgg4.
StochqsticMultiobiective
GenerqtionScheduling
6.1 INTRODUCTION
Optimal economicdispatchin electric power systemshas gainedincreasingimportanceas the cost
associatedwith generationand transmissionof electric energy keepson increasing.The p
involves the allocation of total generationrequirementsamong the available generatingunits in
the system in such a manner that the constraints imposed on different system vari les are
adequatelysatisfied and the achieved overall cost associatedwith it is a minimum.
Despite extensive research focussing on thermal power dispatch problem' muc of the
-state
effort todate has involved the developmentof deterministic models applicable to ste
conditions. Most of these attempts assumeth-esystem data to be deterministic. It m s that
all input information is known with complete certainty and the optimal plans of disl tch are
always realized exactly. In practice, there are several inaccuracies and uncertainties in input
information (Figure 6.1), which lead to deviationsfrom optimal operation.
The operating cost functions representingthe perfortnance characteristicsof then plants
are computed by calculating the overall thermodynamicperformanceof a unit consisting f boiler,
turbine, condenser,heat cycle, and associatedptant auxiliaries. Such cost functions are rate
in most cases.The inaccuraciesmay be viewed due to the following reasons[Kirchma ,
19581.
. Inaccuracies in the process of measuring the basic data' used for compu tion of
thermodynamic performance of the unit
Deviations frorn-the computed thermodynamicperformance of the unit becauseof
errors
encounteredin operation due to opdratingat other than standard'pressureand
Effect of time on equipment conditions which influences some of its operating
stics, notably its efficiency
Inaccuraciesresulting from inability to hold generationat exact desired
. Fuel cost variations
. Load forecasting erfors
. Inaccuraciesintroduced by various types of transmission'lossequations.
Unit data
\ _ ' {
I
I
I
I
I
\ - l
|- ------- - - ' : r - - L -
- - - - -; - -: -- -- - &- - - - - : - -- - - - - - 1
1--
! ,Inaccuracy: : I Uncertainty:
I tnfonnation i i tnformation
L------- _-__l t----
Unit commitmentand
economicdispatch
'-------1 procedure
'
i Inaccuracy: 1,,'
i Modelling (
r-____ .._______--J
,i Inaccuracy:
i Control
t-------
If all these factors are taken together.these will causeinaccuraciesof great magnit in the
steady-stateoperation.The effect of inaccuraciesis in an increasein the'overall cost. Vi i and
Heydt [1981] have outlined the computationaldetails of the stochasticoptimal energy ispatch
problem. The stochasticoptirnal energy dispatch algorithm employed the multivaria Gram-
Charlier series to statisticallymodel the probability density function of the control v or. The
applicability of the series has been limited by the high computationalrequirementsof c ulating
high order statisticalmoments.The method obviatessome of the difficulty through a
transformationof the variatesto be modelled in order to enhancenormality. The main ai of the
method was to produce a tool which would be useful from an operationalstandpointbu fails to
considerstochasticcost function.
Yakin [1985] has articulated an approachto the optimal generationschedulingof power
systemby treating the electricity demandat a node as a random variable with a known I
lity distribution. Particularly, a two-stagestochasticprogramming with recoursemodel been
developedfor stochasticeconomic dispatch An equivalentproblem to this two-stage m I has
been defined. The penalties for discrepanciesin the generation have been incl in the
objective function of the equivalent problem. The major difficulty during the implemen tion is
to draw out the exact values of thesepenalties.
El-Hawary and Mbamalu [1988] have investigatedthe perturbationsin the systemthermal
fuel cost and the system equality constraintsas stochasticand normally distributed w th zero
mean and a given variance.In an anotherattempt,El-Hawary and Mbamalu [1989] in uced a
method in which the system power demand was assumedrandom with zero mean unit
variance.In the third attempt,El-Hawary and Mbamalu (1991) consideredthe pertu lons ln
system power demand as random and normally distributedwith zero mean and some v
stochastic Multiobjective Generation scheduling 389
They observedthat optimality conditionsin terms of the active power generations
were b by
parametersobtalned from the variancesof active power generations.
But these do rpt
provide trade-off betweeneconomy and risk measuresdue to uncertainties
in system tion
cost and randomnessof demand.
Parti [1987] has expounded an gconomic dispatch of thermal generation while
incor ng
the randomness in system production cost and system load through
lenerator outputs, whi were
treatedas random variables.He appendedthe traditional objective firnction
of economic patch
with a penalty term accountingfor the possibledeviationsproportional
to the expectation of the
squareof unsatisfiedload becauseof randomness of generatorpower.This upp.ourhsuppr the
ffue characterof the problem by consideringonly the mon etary aspects
and fails to ex the
non-commensurabilityof the conflicting'objectives.
Besides electric energy, power plants also produce sizeable quantities
of solid wastes,
sludge,and pollutantsthat affect air and water quality.The pollutants
affectingair qualit are of
the greatestinterest.They include particulates,NO' CO' So, and other sundry
oxides of lphur
that can travel over considdrabledistances,and have long-term effects
both in spacear time.
The traditional meansfor controlling emissions,such as precipitatorsand
scrubbers,are h ware
intensive,relatively inflexible and limit ttre ratio of emissionsto energy produced
in eac plant,
but not the total emissionproducedin a region. In contrast,dispatchingr"quires
little I ware
and is flexible and effectiveat the regionallevel. Here, the obje.liu" function
used in dis
can be changedin a few moments.
In general,a large-scalesystemas typified by an electric power system,possesses
r ultiple
objectivesto be achieved,namely economicoperation,reliability, securityand minimal
im on
environment.It may be obvious that tracie-offsamong these objectivesare difficult
of
their different nature.This implies that objectivesare non-commensurable.
6.2 MULTIOBJECTIVE
STOCHASTICOPTIMALTHERMAL
PowER DISPATOH-e,-OONSTRA|NT
METHOD
Extensivestudies,associatedwith the optimal power dispatch,have been centred on
ma ing it
more efficient in algorithm and applicableto online with deterministicdata. In actual
nractice, it
is a misleading assumptionthat data is known with completecertainty. In spite
of this, it is also
,hT the dispatch is to optimize just one specific objective, or single perfo
Iu,9 .ontimal nce
index. Now the trend is to formulate multiobjective optimization pioblem with
due consi on
of uncertaintiesfor a more realistic approach.The multiobjective stochasticoptimization
lem
is describedin the subsequentsections.
6.2.1 StochasticProblemFormulation
The objective function fo be minimized is the total operatingcost for thermal generating
Its in
the system. The operatingcost curve is assumedto 6e approximatedby a quadratic fu I on
of
generatoractive power output as
NG
F r =I (o,4' + biPi+ ) (6.1)
i=l
",
i
i where
, N G is the total number of generators
I
I
I
L
System Optimization
where
NG
s r . - t
Ft= ) tft* c'ilo,1'*6,P,+ (6.4)
f=l
",1
The load demand constraint is
F, +F, (6.s)
NG NG NG NG
NG
Fz= I c'oF,' (6 1lb)
i=l
A multiobjectiveoptimization problem c
operatingcost and (b) the risk associatedwith 1
their expected values while satisfying the el
generationlimits. The multiple objectiveoptimi
Minimize fi (6.13a)
subjecrro Fz S ez (O.tgU)
NG
S
\ D - Pr+Pt
Zy,i- (6.13c)
i=l
where s2 is interpretedas the maximum tolerableobjective level. The values of r are chosenfor
which the objectiveconstraintsin Eq. (6.13) are binding at the optimal solution. As a corfstraint
is varied parametrically,a set of non-inferior solutions (with their correspondingtrade-dffs)is
generated.
The well-known method of Lagrangemultipliers is quite popular in the power ystem
planningstudies.The LagrangianL formed for the systemis given by Eq. (6.14)
NG
_)
L- r; + Ltz(Fz + F o I- ,,)
i=1
dr=;t'.
F, ^"fr.,L'+'-t1
aE ^ AF, faE I
=o (6.15a)
dt = Fzttz=O
dAn (6.15b)
#=Fr+Fo-I"--=o i=l
(6.15c)
(6.16)
R R
+ ( tr)z
tol )z+ (LJ"rz)z + (Y)'p)z+ (Yi'
f, tvr-.)z
I i=l
i=l )='"-[
then GOTO Step 14.
12. Modify control variables,
!F* = p+ Lp
13. Update iteration counter,IT = IT + 1,
Aisign new values to old variablesto continue the process'
J J
I I 0.0w2725 I 2 5l
2 2 0.0003090 1 3 679
3 3 0.0032295 2 3 65
Using this data, the resulting expected generation schedules with expected cost and risk
in Thble 6.4 for various values of e. The deterministic results are shown in Table6.3 (C
j = I, 2 , 3 ) .
Sn PD Fl Pr P2 P3
no. (Mw) ($/tr) (Mw) (MV/) (Mw)
I 140.0 361.263 46.146 54.78r .182
2 180.0 484.010 64.417 68.727 1
220.0 621.695
?'
Sr n
($rn;
F2 FL
(Mw)
Pr
(Mw) (Mw)
P2 F3
(Mw)
no. (MW2)
Fo = 140MW
I 362.2143 7t.4W 5.7648r 47.982 53.900 43.960
2 362.2M0 71.428 5.79242 47.835 53.975 M.060
3 362.1948 71.456 5.81940 47.692 54.048 M.158
4 362.1868 71.484 5.84578 47.554 54.r18 44.252
5 362.1797 71.512 5.87163 47.419 54.r87 M.345
6 362.1736 71.540 5.89695 47.287 54.253 M.436
7 362.1684 7r . 5 6 8 5.92182 47.t59 54.318 M.524
8 362.t64r 71.596 5.94625 47.034 54.381 M.611
9 362.1604 71.624 5.97028 46.9tr 54.443 44.696
10 362.1575 71.652 5.99390 46.791 54.503 M.780
PD= 180MW
11 485.8152 119.880 8.35676 67.572 67.980 52.921
12 485.7668 119.916 8.4331 67.297 68.054 53.210
l3 485.727r 1r9.952 8.52199 67.05r 68.120 s3.470
l4 485.6942 119.988 8.59468 66.821 68.178 53.709
15 485,6663 120,024 8.66279 66.619 68.232 53932
t6 485.e27 120.060 8.72711 66,425 68.282 54.141
t7 485.6222 120.096 8.78810 66.238 68.335 54.337
18 485.6051 r20.r32 8.84665 66.066 68.377 54.526
t9 485.5905 120.168 8.90289 65.902 68.4r8 54.706
20 485.5782 t20.204 8.95712 65.745 68.456 54.879
PD = 220 MW
-
Stochastic Multiobjective Generation
0.45
I
v)
0.4
o
C)
0.35
q) 0.3
\e
0.25
0.2
0.51 0.51060.51120.5il8 0.66640.667 0.68760.83020.8308
0.8314
d +
bi fi.max p.mm
($/Mwh) (Mw) (Mw)
I 0.004100 1.280 200.0 50.0
2 0.002200 0.795 210.0 210.0
3 0.000950 1.809 200.0 10.0
4 0.002145 0.657 400.0 150.0
5 0.001110 0.889 310.0 310.0
6 0.006000 0.300 200.0 100.0
7 0.010400 0.635 100.0 50.0
8 0.006350 0.572 150.0 50.0
The deterministic generationschedulesobtained are given in Table 6.7 for various demands.
Generators2 and 5 remain at their maximum loads of 210 MW and 310 MW respectivel; and are
not included.The expectedcost and risk with expectedtransmissionlossesfor vari expected
demandsare presentedin Table 6.8 which are in the non-inferior set. The correspondin expected
generationschedulesare given in Table 6.9 for various valuesof t.
F D n F r F 3 P 4 F 6 P 7 P8
(Mw) ($nr) (Mw) (Mw) (Mw) (Mw) (Mw) (Mw)
t240 1731.691 137j99 61.047 272.756 145.826 65.254 01.165
1320 1909.131 153.874 84.645 292.947 157.232 71.308 w.476
1400 2w1.048 170.938 107.708 313.463 169.036 77.534 r7.916
Sr Fl F2
no. ($n1 (MW2)
I 175r.636 2554.OO .87011
2 1751.4M 2555.4 .89693
3 175r.175 2556.88 .9238r
4 1750.951 2558.32 .95073
5 1750.731 2559.76 .97770
6 1750.51,4 256t.20 .00471
7 1750.302 2562.64 .03177
8 1750.093 2564.08 05887
9 1749.887 2565.52 .08601
t0 1749.685 2566.96 .11319
Stochastic Multiobjective Generation Sche
Sr Fr F2 FL
no. ($n; (MWz) (Mw)
1 1993.124 2682.000 64J5024
2 1990.328 2683.600 64.30595
3 1987.785 2685.200 64.27084
4 1935/48 2686.800 64.24320
5 hBst.zso 2688.400 .22176
6 t98r.271 2690.000 .20558
7 r975.903 2694.822 .17964
8 r974.342 2696.400 .1803
I
Sr Fl FL
no. ($rn; (Mw)
I 2195.333 2950.80 .79;334
2 2191.831 2952.56 72950
3 2188.678 2954.32 .67876
4 2185.806 2956.08 .63847
5 2183.t64 2957.84 70
6 2180.7t4 2959.60 .58204
7 2r78.429 2961.36 .56348
8 2174.271 296/..88 54130
9 2r72.365 2966.64 36lt6
Sr Pr P3 P4 P6 P7
no. (Mw) (Mw) (Mw) (Mw) (Mw)
194.355 r48.276 r9t.225 169.693 127.r85 1 3 .552
1 9 4 . 1I I r48.103 183.135 169.944 125.534 l 3 .393
193.873 147.927 t84.929 170.t63 124.039 l 3 .237
t93.639 r47.748 186.629 170.357 122.668 l 3 .085
193.410 r47.566 t88.252 t70.529 121.400 I .93s
193.184 147.383 189.809 170.683 t20.2r9 I .787
192.962 147.197 191:310 170.822 119.113 I .642
192.526 146.820 194.167 1,7t.O62 117.088 I .357
192.3t2 146.629 1,95.534 t7t.t67 t l:6.t55 I .2t8
OPTIMALTHERMALPOWER
STOCHASTIC
6.3 MULTIOBJECTIVE
DISPATCH-THESURROGATEWORTHTRADE-OFFMETHOD
In this section,the effect of uncertainsystemparametersis incorporated,explicitly in the
multiobjective power dispatch. Multiobjective problem is stated b$ considering (i) the ex
operatingcost, (ii) the expectedminimum NO, emission,(iii) the expectedtransmissionI , and
(iv) the expected dertiations becauseof the unsatisfied demand. The surrogate worth tr ff
techniqueis discussedto find the compromisedsolution.
Stochastic Multiob.jective Generation Scheduling
4.5
t 3.5
I
g 3
o
Er.t
' E z
GI
T,
s
1.5
0.5
1.6 t.602 1.604 1.606 1.608 1.61 t.6r2 t.6t4 1.6t6 1.618
d, -+
4.6
A
II
8
o
4.t
.E
tr
o
. 9 4
'oq)
Bq
3.7
# 1400MW
Figure 6.4 Percentage deviation in cost vs. q,.
6.3.i MultiobjectiveoptimizationProblemFormulation
The multiobjective optimization problem is viewed as a stochasticmultiobjective optimization
problem by consideringthe systempower demand,cost coefficients,NO, emission coefficients
Power System Optimization
NG
F2 = r(1+4) a,1'*a,F,+j,1 (6.18b)
i=l
where 4min and 4t"* are expectedlower and upper limits of generatoroutputs, respectively.
Since generatoroutputs Pi are norrnally distributedindependentrandom variables,so the
expecteddeviations are proportional to the expectationof the square of the unsatisfiedload
demand.These expecteddeviations are consideredas another objective to be minimized and the
objective is given as
F4 (6.2r)
Stochastic Multiobjective Generation Scheduling
var(4 )
or
NC
Fq= Z c'oF,' (6.22)
i=l
6.3.2 SolutionProcedure
To generatenon-inferior solutions to a multiobjective optimization problem, ffre e-constraint
method is utilized [Haimes, 1977].The t-constraintapproachreplacesthree objective functionsto
consffaintsas given below.
Minimize Fr 6.Zaa)
subject to F,3e, ( / = 2 , 3 ,. . . 4, ) (6.24b)
NG
FT
L P , - P D- 0 (6.24c)
f=l
where €; is the maximum tolerable objective level for the 7th objective.
Generationof non-inferiorsolutions I
(6.26a)
AL FJ ; - t ;J = 0 (j=2,3,4) (6.26b)
dhtj
N
AL \-r
PD-LP,-o (6.26c)
=--
dl,L
i=l
where
a4 = 2(r+ cf;,)aiFi
+ ui (6.26d)
dP,
9 5 =2C2nFt (6.269)
dP,
The Newton-Raphsonmethodcan be appliedto solve the abovenonlinearequations.To implement
the Newton-Raphsonrnethodthe following equationis solved iteratively.
(6.21)
R 4
13. Check either
I ro",lt * I (aLr)z+ (Ltt)z
i=l j=2
f,new= 4+44 ( i = I , 2 , . . . ,R )
L#'"-)W+L,hry ( / = 2 , 3 , 4 )
[tn* = pt+ L,yt
15. Update the iteration counter,IT = IT + l.
Assign new values to old value variables
Pi = Pin"*(r = I , 2, ..., R)
Lu=Lfr"nand p-Ix**
GOTO Step I I and repeat.
16. check the limits of generatorsand fix up as following
If
4.4minthen 4=4min
If 4 , Pr** ,h"n Pi = P,**
The systemgiven by Eq. (6.25) is solvedfor K valuesof €2,say, e1,,.. ., where elo and
.8f,,
Ff
ef; are held at some level sr9.S"t initial values of e; such that e; > ry" and er'<
= et & - 1,
only those values of 2!rz> 0, which correspondto the active constraint F!
Z, ..., K) are consideredsince they belong to the non-inferior solution. Similarly the trade-
=
off function 213is generated,where Eq. (6.25) is solved for K' different valuesof ef (k
is generated.Regressionanalysisis
1,2, ..., K'), wittr fixed level ef and e!. Simil*!, .1,1a
performedto yield the trade-off functions hnlFzl, hy,IQl and )'slFal.
3. Generation of SWT function: SWT function assignsa scalar value (on an ordinal scale)
to any given non-inferior solution. One way of specifying non-inferior solution is by
trade-off functions.Moreover, there is close.relationshipbetweenthe SWT function 1V1.1 and
the partial derivativesof the utility functions. In multiobjectiveanalysis,it is assumed
implicitly that the DM maximizeshis utility which is a monotonic decreasingfunction
of the objectivefunctions.Given a decisionvector F and the associatedconsequencesFi,
the utility is given by
( J= [ J [ 4 , F r , F, F o ] (6.2e)
By linearizing the utility function for a small changein 4, ,h" following can be obtained
[Haimes,19771.
(6.30)
The SWT function Wry is a monotonic function of Uy, with the property that WU = 0
<-> UU = 0 and is written as
1VU= hi Uti (i = 2, 3, 4) (6.31)
-'id
Stochaltic Multiobjective Gencration Scheduling ili l
5. The optimalset of decisionvectoris found by solvingthe following problem.
l--3 I
Minimize 632a)
la.Idr4 l
L i = 2 J
NC
subject to
In-Fp=o
i=l
(6.32b)
L
5F
=- 4 * $ ^ . & - ( r =l ' 2 " "N' G ) (6.34a)
tr. htrifr-F=Q
a =
=ruL FD-I
; N q =o (6.34b)
d=l
=[-;;]
l1rrI';J[il] (6.35)
=
#=#.L^"# ll ' ' 2
" '
N G) (6.36a)
dzt g
\( ie= -1 , z , : . . , N Gj *; i ; j = t , 2 , . . . , N G ) (6.36b)
W=*n@
azt azt (r=1'2""'NG) (6.36c)
ffi=ffi=-l
,'L^= o (6.36d)
dp'
,fOS- Power SYstem
Utllity functlon
generalutility function for a given DM that can predict
Here, it is assumedthat there exists a very
utility function be defined for each objective function
his behaviour and interest. Let the DM's
function F, to the other objective functions' So,
dependingon the importance of the objective
.rotaror overau utility function is defined as [osyczka and Davies, 1984].
1
rr= t k!F, (6-37)
Maximize L'-
t1-t
' L
t=1
to the
The solution vector P is then found bY maximizing the total utilitY subjected
technology constantsas defined below.
NG
tL P,-Pp=o (6.38)
i=l
Further define
(6.a0a)
such that
4
I o; = t; ki>o
(6.40b)
i=l
the aboveoptimiza-
The DM gives the weight w;on the attributebetween(0,99). The solutionof
tion probiem can be obtained as explainedin Section 5'4'4'
6.3.4' SamPleSYstemStudY
its applicability' The'
The method is applied to a six-generatorsample system to demonstrate
expected fuel cost characteristics($ltrl undertakenfor study are as:
at
Table 6.1I shows the conflicting objectives, trade-off functions, utility and SWT function
eabh non-inferior set is shown in
each non-inferior set. The decision vector P; correspondingto
Table 6.ll Expectedcost, emission,risk and transmissionloss along with utility and SWT
function, when demand is 200 MW
Sr. Fl F2 F3 F4 )"t, Wn
no. ($/h; (ke/h) (Mw) (MW2) ($/kg;
74.2663 37r.5998 542.1307 -5
I 1306.917 20.6380 r.9145
1305.927 20.&34 1.9096 73.9515 395.297r 54r.6717 4
2
1304.951 20.6559 1.9071 73.7545 429.8842 541.V+39 4
3
4 1303.951 20.6776 r.9076 73.7r25 no.ss21 540.8398 4
73.9955 500.6307 540.4&5 :7
5 1302.847 20.7r8r 1.9140
74.1388 5W.7724 540.0906 :7
6 1301.824 20.7492 1.9r78
74.0979 500.6307 539.6841 :7
7 1300.818 20.7684 1.9187
1299.788 20.7950 r.9219 74.1986 473.6248 539.2983 4
8
74.1169 415.0005 538.9263 -5
I 1298.892 20.8104 t.9215
20.8202 r.9207 73.935r 328.1359 538.5624 4
10 1298.068
73.8146 22r.8559 538.z;07 -3
ll 1297.315 20.8361 t.922r
73,8709 99.3539 s38.0527 -l
t2 1296.801 20.8&9 1.9263
(Contd.)
410 Power System Optimizntion
Sr.
no.
n
($/h1
F2
(ke/h)
F3
(Mw)
n
(MW2)
LB
($nt{w1
wn
F1 Fz F3 F4 )qt Wro
($/h1 Csft) (Mw) (MW2) ($A{w1
I t287.16r 2r.6ffi3 2.&Or 88.7065 27.O5r7 537.4r8r 90
2 1286.748 2t.6ffi 2.4000 88.8003 2r.7904 537.2714 50
3 t286.421 2r.6600 2.4000 88.9000 r8.0675 537.r605 3
4 t286.r56 2r.6ffi0 2.40m, 89.0000 15.2434 537.0745 2
5 1285.938 21.6600 2.40W 89.1000 t2.9794 $7.m71 I
6 t285.757 2r.6ffi0 2.4W 89.2000 r r.0854 536.9547 I
7 r285.607 2r.6ffi 2.4000 89.3000 9.M7r 536.9r47 0
8 r285.485 ?1.6600 2.4p/00 89.4000 7.9907 536.8859 0
9 1285.387 21.6600 2.4000 89.5000 6.6679 536.8669 0
Table 6.12. Iri this'case, the coefficientsof variations of cost, emission and B-coefficientsare
assumedlcfo.The scalar weights are 40.0, 2O.0,20.0 and 20.0 for cost, emission,and power lot;s
and risk objectives in sequence.
By regressionanalysis,the trade-off functions for 20O MW are representedas
Sr Pl Pz P3 P4 Ps P6
no. (Mw) (Mw) (Mw) (Mw) (Mw) (Mw)
I 18.202 28.223 35.16s 33.Mr 29.t02 55.267
J 18.599 28.1I7 35.050 3 5 .I 6 5 28.432 54.637
3 18.684 28.080 34.243 37.t03 27.990 53.900
4 18.458 28.073 33.313 39.318 27.781 53.057
5 17.918 28.O49 32.223 41.989 27.569 52.252
6 17.803 28.291 30.952 44.467 27.387 51.099
7 18.316 28.726 29.443 46.652 27.210 49.654
8 19.r49 29.088 27.750 48.885 26.953 48.174
9 20.870 29.488 25.949 50.623 26.512 46.558
10 23.612 29.855 24.W6 51.751 25.690 45.027
ll 27.294 29.855 21.944 52.352 24.977 43.579
12 3t.205 29.465 19.976 52.523 24.633 42.198
F2 Ft F4
FD FL
(Mw) ($/h) Cs/h) (M!U Ct"tW'l
Minimumcostdispatch
200.0 rzu.365 21.8230 2.4749 93.8997
400.0 1788.286 u.8757 8.7164 355.I 199
600.0 2386.325 37.0950 19j952 813.4999
Minimumemission disPatch
1318.079i 20.2491 t.7949 68.6702
200.0
1878.364 20.6196 6.8843 268.3033
400.0
600.0 2577.228 27.7313 15.3691 605.5116
Ma,rimumutility approach
200.0 r3w.79l 20.3M7 r.7755 6692t5
400.0 1856.690 20.9385 7.2144 268.r749
600.0 2575.991 28.2854 15.3633 600.4974
PD PL P2 P3 P4 Fs P6
made it possible
technique.The novel formulation as a multiobjective optimization problem has
to quantitatively grasp trade-off relations among conflicting objectives.
The trade-off approachis effective only up to two objectives;as the number of objectives
SWT has
increases,the selectionof the best solution becomescumbersome.An interactivemethod
power dispatch
been applied ro identify the best compromised solution for multiobjective
problem, when conflicting objectives are more than two. The major characteristicshnd advantagen
maker's
of the SWT method are that the surrogate worth functions, which relate the decision
preferencesto the non,inferior solutions through the trade-off functions, are constructed in ttre
iunctional spaceand only then are transformedinto the decision space.
Stochastic Multiobjective Generation Scheduling 413
OPTIMAL
STOCHASTIC
6.4 MULTIOBJECTIVE POWER
THERMAL
METHOD
DISPATCH_WEIGHTING
The economic dispatch problem was defined so as to determine the allocation of electricity
demand among the committed generatingunits to minimize the operating costs subject to physical
and technologicalconstraints.Most of the existing formulations of the economic dispatch are
solved as static deterministicoptimizationproblems.Actually, there are many inaccuraciesand
uncertaintiesin the input information which lead to deviationsfrom optimal operation and cause
an increasein the cost over the optimal value [Edwin and Machate, 1980]. As a result of the rise
in production costs due to uncertainfactors,the electric energy systemhas been representedas a
network characterizedby random variablesand investigatedby numerousresearchersat various
levels [Dopazo et al., 1975;Parti et al., 1983;and El-Hawaryand Mbamalu, 1991]. Although,
these approacheshave been successfulin applicationsinvolving stochasticeconomic dispatch,
but all the methods do not provide trade-off between economy and risk measuresdue to
uncertaintiesin systemproductioncost and random natureof demand.Typically, such conflicts
exist becauseno such feasiblesolution has been found which would minimize them all-
The pollution minimization problem has attracted a lot of attention due to the public
demandfor clean air. Thermal power stationsare major causesof atmosphericpollution, because
not
of high concentration of pollutants they cause. Since optimum economic dispatch is
environmentallythe best solution, many organrzations in their fight against air pollution have
come up with a new method,the so-calledminimum emission dispatch(MED). MED is used to
minimize the total stack emission(NOr) for the entire system,althoughthis may be controlled
either through post-combustioncleaningsystems(electrostaticprecipitators,stack gas scrubbers)
set
or automatically(controllingunit loading).MED rnay be obtainedby introducinga different
of generatorrepresentations into the economicdispatchproblem. The MED generatorequationfor
($fn)
each unit is a function of stack emission(NOr) versusmegawattoutput, insteadof input
versus megawatt output for the economicdispatch.
6.4.1 StochasticMultiobiectiveOptimizationProblemFormulation
units in
The objective function to be minimized is the total operatingcost for thermal generating
the system and a quadraticoperatingcost curve is assumed'
NG
F r =I ( a i P i 2+ b i P i * c ; ) (6.41)
i=1
where
a;, bi, and c; uto the cost coefficients
NG is the total number of generatingunits'
and load
A stochasticmodel of function F1, is formulated by consideringcost coefficients
be convertedinto
demand as random variables.By taking expectation,the stochasticmodel can
distributed and
its deterministicequivalent.The random variablesare assumedto be normally
function may be
statistically dependenton each other. The expectedvalue of operating cost
1987]' By (aking
obtainedthroughexpandingthe function using Taylor'sseries,aboutmean [Rao,
by
the expectationof the expandedform, the eipected operatingcost function is represented
NG
F =I (a;,
[(o,F,'+6iPi+ c, +d, var(4)+ cov(b,,4) + 2F, cov P;))
(6'42)
i=l
414 Power System OPtimization
where
Ai, bi, and Ei are expectedcost coefficients
P, is the expectedvalue of the ith generatoroutput.
In this study, varianceand covarianceare replacedby coefficientof variation (CV) and
correlationcoefficient(CC), respectively.In general,varianceand covarianceare defined as
var()Q= C'* X' (6.43)
cov(X,If- RxvCxCyXY (6.44)
where
Cy and Cy are the CV of random variablesX and I, respectively.
Rlgyis the CC of random variablesX and Y.
The value of CC is positive or negativedependingupon the sign of the covariance and itri
value lies between-1.0 and 1.0.
Using Eqs. (6.43) and (6.44),Eq. (6.42) can be rewritten in the simplified form as
NG
Fr = fft *C'n*ZRo,r,Co,Cn)o,F,'
+(1 + Ru,r,Cb,Cp)biFi+c;) (6.45)
I
i=l
For the fixed network configuration and random load demand, the equality constraint in the
classical dispatch problem is representedby the expected power balance equation stated as
NG
= Po+Pt (6.46"\
It-
i=l
where P, and FL are the expected load demand and the expected transmission loss, respectively.
Expectedtransmissionloss
The transmissionpower loss expressedthroughthe simplified well known loss formula expressiort
as a quadraticfunction of power generationsis given by [Kusic, 1986]:
NC NG
Pr= I I PiBuPj (6.47"1
i=l j=l
Power generations Pi are dependentrandom variables . Bil are also consideredas inaccuratt:
B-coefficients.The expectedtransmissionlossesusing Taylor'sseries are representedas
NG NG NG-I NG NG
FL= I > FiEuFj+)E,,vu(4)+
I I 28, cov(Pi,P)+ > 2\ cov(Pi,Bii)
i=l i=7 j=i+l i=l
NG NG
+> I 2F,cov(pi,Bu) (6.48)
i=l ,t=l
j#i
Stochastic Multiobjective Generation Scheduling 415
NG NG NG
Fr= I ft * c7,+ 2Rrr.cncr,,lE,,F,'*
f lft * Rn,p,cncr,
+ 2R4Buc4gn)FiEuFj
6.49)
i=t i=l j;i
Expected deviations
Since generatoroutputs P;s are treated as random variables,the expecteddeviations are propor-
tional tb the expectationof the squareof the unsatisfiedload demand.Theseexpecteddeviations
are consideredas the secondobjective to be minimized. The secondobjective function Fz is
representedas
NG
Fz= +Fr-I (6.s0)
i=l
which on simplificationreducesto
NG NG_I NG
Fz=Iur(4)+I \zcov(P,,Pj) (6.51a)
i=l i=l j=i+L
NG NG NG
F2 = Z r'oF,'
*> I Rnr,cncr,F,Fi (6.s1b)
i=l j=l i=l
j*i
NG
subjectto
I "_'= Fo+F, (6.s2b)
i=l
<Fi <fimax
P,min ( f = 1, 2 , . . . ,N G ) (6.52c)
NG
subjectto
>i = l 1 = Fo+F, (6.s3b)
where lr1 &r€ the levels of the weighting coefficients.This approachyields meaningfulresult to
the decision maker when solved many times for different valuesof w1r,k = I,2. The valuesof
weighting coefficientsvary from 0 to l.
To solve the scalaroptimizationproblem (6.53), the Lagrangianfunction is defined as
(_ $q _)
L ( P i , , Ar)r=F r + * 2 F 2 * 2 l F r + P L -pI ,I f6.54)
\ i = l )
where 2 is Lagrangianmultiplier.
The necessaryconditions to minimize the unconstrainedLagrangianfunction are:
The scalar optimization problem is solved using the Newton-Raphsonalgorithm. The size of the
formulated Hessian matrix is the same as that for the deterministicproblem, becausesingle
objective function is solved in both caseswith the same number of constraints.To implement the
Newton-Raphsonmethod,the following equationis solved iterativelytill no further improvement
in decision variablesis achieved.
=[-l;1
l1r^Y,:)[il] (6.s6i)
with respect
Hessianmatrix elementscan be obtainedfrom the above equations by differentiating
to control variables,one by one.
( i = 1 , 2 , . . .N
, G) (6.57ar)
Stochastic Multiobjective Generation Scheduting 417
aPtaPj anaPj
:
dzt (6.57c)
anu"
azt (6.s7d)
fr
Algorithm 6.3: Non-Inferior Solutions by the Weighting Method
l. Read data, namely cost coefficients,emissioncoefficientsand B-coefficients,demand,Err
(convergencetolerance) and ITMAX (maximum allowed iterations), M (number of
objectives),NG (numberof generators)and K (numberof inferior solutions),etc.
2. Set iteration for non-inferiorsolutions,k = t.
3. Incrementcount of non-inferiorsolutions,k = k + l.
4. If (/<> Kl GOTO Step 17.
5. Feed or generateweights,w; [i = I ,2, ..., M (hereM = 2)]
6. Computethe initial valuesof P; (i = I ,2, ..., NG) and 1, by presumingthat Fr = 0. The
valuesof ,1.and 4 O = 1,2,..., NG) canbe computeddirectlyusingEqs.(3.10)and (3.9),
respectively.
7. Assume that no generatorhas been fixed either at lower limit or at upper limit.
8. Set iteration counter, IT = l.
9. Compute Hessian and Jacobianmatrix elementsusing Eqs.(3.56a) and (3.56e) and
Eqs. (3.57a) to Eq. (3.57d),respectively.Deactivaterow and column of Hessianmatrix
and row of Jacobianmatrix representingthe generatorwhose generationis fixed either at
lower limit or at upper limit. This is done so that fixed generatorscannot participatein
allocation.
10. Gauss elimination method is employed in which triangularizationand backsubstitution
processesare performedto find L4 (i = 1,2,..., R) and A,?".Here R is the number of
generatorsthat can participate in allocation.
4 n ' *= 1 * t E ( i= 1 , 2 , . . . , R )
2,re*=),+ LJ"
6.4,3 DecisionMaking
Consideringthe imprecisenatureof the DM's judgement,it is naturalto assumethat the DM may
have fuzzy or imprecisegoals for each objectivefunction. The fuzzy sets are defined by equations
called membershipfunctions.Thesefunctionsrepresentthe degreeof membershipin certainfuzzy
setsusing valuesfrom 0 to 1 [Klir and Folger, 1993].The membershipvalue 0, indicatesincom-
patibility with the sets,while 1 denotesfull compatibility.By taking accountof the minimum and
maximum valuesof each objectivefunction togetherwith the rate of increaseof membership
satisfaction,the DM must determinemembershipfunction tt(F) in a subjectivemanner.Here it is
assumedthat p(4) is a strictly monotonicdecreasingand continuousfunction defined as
I ;4<4*"
ri
Fmax - r ir.l
. f , m i n <t lttri ,. < f l m a x
It(Ft) = .-p , . m a-x ; 4ttn (6.s8)
' t '' lm m
0 ;4>Pmax
The value of membershipfunction indicateshow much (in the scale from 0 to 1) a non-inferior
(non-dominated)solution has satisfiedthe Fi objective. The sum of the membershipfunction
values (tt(F) (l = ,1,2, ..., It[) for all the objectivescan be computedin order to measurethe
accomplishmentof each solution in satisfyingthe objectives.The accomplishmentof each non-
dominatedsolutioncan be ratedwith respectto all the K non-dominatedsolutionsby normalizing
its accomplishmentover the sum of the accomplishments of the K non-dominatedsolutionsas
follows [Thpiaand Murtagh, l99l]:
Ii=l
F(F)K
tt|= (6.se)
F(F)K
,_..J
Stochastic Multiobjectivq Generation Scheduting 4le
The function IrB in Eq. (6.59) can be treatedas a membershipfunction for non-dominated
solutions in a fuzzy set and representedas fwzzy cardinal priority ranking of the non-dominated
solutions.The solution that attainsthe maximum membershipttt in the fuzzy set so obtainedcan
be chosenas the best solution or that having the highestcardinal priority ranking.
M a x { p * t k = 1 , 2 , . . . ,K l (6.60)
Generator q. bi c; pmax
f,min
no. ($/IvIW2tr; ($/MWh) ($ltr; (Mw) (Mw)
I 0.012 6.51 135.33 105.0 30.0
2 0.005 5.627 261.19 225.0 50.0
3 0.006 5.506 264.63 250.0 70.0
The parameterstatisticsis known from past history. OtherwiseMonte Carlo simulation technique
is a useful tool in simulating the parameterstatistics.Analysis of variation in expectedcost seems
necessarybecauseof the existenceof covarianceof two random variablesin the problem formula-
tion. In the study, covarianceof bivariaterandom variablesis consideredpositive (increasing)or
negative(decreasing)pairwise. Covarianceof one pair of random variablesis consideredat a time,
whereasthe rest of the random variablesare consideredindependent(uncorrelated).All correlation
coefficientsare variedfrom:-1.0 to 1.0 in stepsof 0.5. Assumingboth the weightsw1 and w2 as 0.5,
the percentagerelative changesin Fi from deterministic FL with respectto Rqn, Ru,n,RrB- and
Rnr, (i # i) are shown in Figure 6.5 for the expecteddemandof 350 IvtW. It can be observedthat
(i) there is increasein the percentagerelative cost as Rnp, is changedfrom positive value to
negativevalue, and (ii) there is decreasein the percentagerelative cost as CC of rest of the
random variables are changed accordingly.
420 Power System Optimization
Minimum valuesof the objectives F1and F, *" obtainedby settingweights ln1 to 1.0 and wy to
0.0 and vice versa [Dhillon et al., lg93). Owing to the conflicting natureof the objectives,I will
have maximum value, when 4 ir minimurn. The respectiveminimum and maximum values of
objectivesthus obtained are given below when PD is 350 MW.
max
pimin= 3013.519($/h), F = 3063.917($fn;
8
(J
L
o
()
a)
O.
x
o
c l.)
'5
't
q.)
\J
6 R 1
-1 -0.5 0 0.5 1
Correlation coefficients
--# C a s e1 --l r-' C aseZ " ' {-" ' C ase3 -' + ' - C ase4
Over and above, an experienceddecision maker can adopt any other suitablemethod to selectthe
minimum and maximum values of the objectiveswithin which he is expectingthe compromisecl
solution, like minimizing and maximizing each objective function separatelysubject to tht:
required constraints.
of optimalpreferredsolution
Determination
Four different casesare consideredto realize the effect of covarianceof the random variablestcr
eachother (pairwise).
Case 1: In this case all random variablesare considereddependenton each other pairwise.By'
varying weights w1 (from 1.0 to 0.0) and w2 (from 0.0 to 1.0) with a regular decrementandl
incrementof 0.1 respectively,4 and Fz as computedfor non-inferior schedulesare shown inL
StochasticMultiobiective Generation Scheduling 421
Sr W1 Fr p(n) It(Fz) PB
no. ($rn1
1 1.0 0.0 3013.519 1267.589 1.00000 0.00000 0.07702
2 0.9 0.1 3013.612 1265.897 0.99815 0.06044 0.08153
3 0.8 0.2 3013.963 1263.94r 0.99119 0.13030 0.08638
4 0.7 0.3 3014.734 t26t.658 0.97589 0.21185 0.09148
5 0.6 0.4 3016.204 1258.96r 0.94672 0.30818 0.09665
6 0.5 0.5 3018.876 1255.734 0.89371 0.42345 0.10145
7 0.4 0.6 3023.724 1251.824 0J975i 0.56312 0.1047e
8 0.3 o.7 3032.785 1247.033 o.61773 0.73425 0.10413
9 0.2 0.8 3050.861 1241.166 0.25906 0.94381 0.09264
10 0.1 0.9 3056.911 1239.891 0.13902 0.98936 0.08691
l1 0.0 1.0 3063.9r7 1,239.593 0.00000 1.00000 0.vlT02
x
'3, z
€
F
C!
v)
o
; r.s
C)
(.)
o
p
xo
. E 1
E
o
.E
€ o.s
bq
0
0.6 0.7 0.8 0.9
Weightageto cost
-f-- Cost --*-- Risk
Figure 6.5 Percentage deviations in expected cost and risk with respect to weight, w1.
FD FL F2 FL PL P2 P3
(Mw) ($n; (MW2) (Mw) (Mw) (Mw) (Mw)
CaseI
350.0 3023.724 1,251.824 3.8113 82.756 t29.523 t4t.532
400.0 3412.855 1641.559 5.1616 93.337 150.823 161.001
450.0 3816.886 2084.501 6.5633 105.000 169.608 181.956
500.0 4232.638 2595.893 9.4991 105.000 187.694 2r6.805
CaseII
350.0 3022.618 125r.642 3.7855 82.r60 r29.996 141.630
400.0 3418.477 1635.635 4.4299 99.79r 145.026 159.6t3
450.0 38r5.815 2083.319 6.4339 105.000 t69.769 181.665
500.0 4230.941 2593.945 9.3079 105.000 187.889 216.419
CaseIII
350.0 2996.655 1251.488 3.7637 82.396 r29.705 r41.663
400.0 3380.846 1641.037 5.097r 92.917 151.043 161.138
450.0 3779.968 2083.282 6.4298 105.000 169:582 181.848
500.0 4189.928 2593.896 9.303r 105.000 r87.639 216.664
CaseIV
350.0 2996.306 436.72r 3.9232 82.219 t35.354 136.351
400.0 3379.355 577.308 5.4529 90.929 t57.r?L 157.403
450.0 3776.r31 739.849 7.3192 99.080 179.538 178.702
500.0 4186.835 929.535 9.7847 105.000 203.656 20t.128
-.-*{
Stochastic Multiobjective Generation ^Sc
Case 3: Cost coefficients and power generationsare consideredindependentto each othgr (Ro,h
= 0, Ru,n=0; f = 1,.2,3) along with the consideration that power generationsand B-coefficients
are also independent -
( Rntu 0; l = 1, 2,3, j = 1,2,3), whereasexpectedpower generations EIre
dependenton each other. 'Best' solutions for this case are given in Table 6.L9 for the same
demandsas in the case 2.
Case 4: All random variablesare independentof each other. In this case (R*n - 0, fa,4 = 0,
R n u ,=, 0 ; i = 1 , 2 , 3 , i = I , 2 , 3 ) a n d R n r ,= 0 ; i = 1 , 2 , 3 , i = L , 2 , 3 , i * i ) a r et a k e nT. h e
obtained 'best' solutionsare given along wiih the expectedschedulesand transmissionlossesin
Table 6.19 for various demandsunder consideration.It can be concluded from Figure 616 that a
decreasein the weightageto expectedoperatingcost leads to a decreasein the expectedrisk and
an increasein the expectedoperatingcost. Comparisonof four casesclearly-.showsan increasein
expectedcost as the correlation (dependence)of bivariate random variablescomes into qxistence
(CasesI, 2, and 3) comparedto independentrandom variables(case4). Expectedrisk is always in
conflict with the expectedoperating cost in all the cases.Expected risk is lorver in the case 4
where random variables are independentof each other. Case 2 gives very small decreasein
expectedoperating cost as well as in expectedrisk comparedto the Case 1. This is bepauseof
fixed network configuration. Case 3 gives significant decreasein the expectedcost compared
to the Cases 1 and 2. However,there is insignificant decreasein expectedrisk comparedto the
ones obtained in Cases L and 2, respectively.This is attributed to the fact that the objective
function Fr is a separablefunction and Ro,n, R4Oand Rnr, have no effect on the objective Fz
Considering,the results of all four cases,Table 6.19 shoWsan increasein the expectedoperating
cost as the correlation (dependence)of each random variablesincreases.
LOAD DISPATCH
6.5 STOCHASTICECONOMIC-E]VIISSION
The economic-emissionload dispatch(EELD) problem is a multiple non-commensurable objective
problem that minimizes both cost and emission together.In this section, a stochasticEELD
problem is formulated with the considerationof the uncertaintiesin the system production cost
and random nature of load demand.In addition,risk is consideredas anotherconflicting objective
to be minimized becauseof random load and uncertain system production cost. Validity of the
method has be6n demonstratedby analyzinga sample system consisting of six generators.
ProblemFormulation
6.5.1 StochasticEconomic-Emission
In this section, the multiobjectives with the equality and inequality constraintspertaining to the
power system optimization problem are described.The important non-commensurable objectives
taken into account here are:
. Economic operation
. Minimal impacts on environment
. Expected deviations due to unsatisfiedload
4U Power Systgm Optimization
where Co, and Cr, are the coefficientsof variation of random variablesd, and ei respectively.
R4p,is the correlationcoefficientof randomvariablesdi andP;. Rr,4 is the correlationcoefficient
of the random variablesei and Pi.
The expecteddeviationsdue to variancein power mismatchare consideredas the third
objectiveFr, to be minimizedand are given by Eqs.(6.51a)and (6.51b),
Equalityand inequalityconstraints
When the network configurationis fixed and the load demand is random, then the equality
constraintis imposedto ensurereal power balance.This is statedas
NG
Fo+FL-Zn-o (6.64)
t:l
f l m i n- <' I E- '<I f i m a x
- t ( i = 1 , 2 , ,. . . ,N G ) (6.6s)
where 4'in and 4maxare expectedlower and upper limits, respectively,of generatoroutput.
Stochastic Multiobje ve Generation Sch
Expectedtransmissionloss
The transmissionpower lossesexpressedthroughthe simplified well known loss fi la expres-
sion as a quadratic function of the power generationare given by
NG NG
Py- I > PiBuPj (6.66)
i=l j=l
where Ft, Fz, and 4 are the expectedvaluesof the objective functions to be mini zedoverthe
set of admissibledecision variable F,
6.5.2 SolutfonApproach
The weighting method is used to generatethe non-inferior solutions of the tiobjective
optimization problem. In this method the problem is convertedinto a scalar optimi on problem
and is given below as
Minimize (6.70a)
subject to (6.70b)
Z * o =l , w * >o- (6.70d)
k=l
426 Power System Optimization
L(4, L) = (6.7r)
* * F , , + ^ ( r , +;Ft l,
where ,1,is Lagrangianmultiplier.
The necessaryconditionsto minimize the unconstrainedLagrangianfunction
(6.72a)
(6.72b)
where
G + C24+ ZRo,oCo,C
4) + (l + Ru,,Cu,C
e bi (6.72c)
",4
a E
= (t+
-
- - Ciri +2R44Co,C4)7,4
+ (t + R,,nC,,Cr,)
ei (6.72d)
dP,
+d 4 = zc,on.X
'i' R4r,c4cpjFj
'! ' (6.72e)
E,
aF, NG
= 2(r* cA)8,,1+ [t* R4,r,
c4cr,]EuFi (6.72f)
i )j=l
j*i
The scalar optimization problem is solved using the Newton-Raphson algorithm. The
si of the
formulated Hessian matrix is the same as that for the deterministic problem,
becau single
objective function is solved in both caseswith the same number of constraints.To impl
the
Newtonr-Raphsonmethod, the following equationis solved iteratively till no further im1
vement
in decision variablesis achieved.
# = * . r W . ^ ( # - r ) ( , = , ,,2N, G ) (6.74a)
Stochastic Multiobjective Generation Sched, ling 427
dzt ( i = I , 2 ,. . . , N Gj ;= r , z ,. . . , N G ;* t) (p.7+a1
+ '+nI +
a F 4 = w z* @ ffii
a2r dFt -r ( i=1 , 2 , " ' , N G ) (6.74c)
ffi,= fr
* = Q (6.74d,)
ar
where
+ czflE;i (6.741)
a2F =
4 2Q+Rnr,Cncp)Eij Q*j) (6.74i)
E4E,;
Algorithm 6.3 can be used to generatethe non-inferiorsolutions.Use Eqs. (6.72) d (6.74)to
obtain the Hessianand Jacobianmatrix elements.
Ftt = 0.15247
4' + 38.53973
Pr + 756.79886
4, = 0.ros87F] + 46.r5s16F2+ 4s1.32513
Fir = 0.028$F? + 40.396.554+ r049.9s770
4o = o.o3s46F?+ 38.305fiF4+ 1243.53110
F,r = O.OzrttFl + 36.32782
Ps + 1658.5
6960
4u = o.or7zgF? + 38.2704r
F6 + t3s6.6sgzo
428 Power SystemOptimization
Fzt = o.oo4
D n, + 0.32761P1 + 13.85932
Fzz - o.oo4rs
F] + 0.32767P2 + 13.85932
Fzt = 0.C1683
Fr2 + 0.5a551P3 + 40.2669
Fz+ = 0.00683
F+2+ 0.54551Pa + 40.2669
Fzs = 0.00461
Ps2+ 0 . 5 1 1 1 P
6 s + 42.89553
Fru = 0.00461
Fo2+ 0 . 5 1 1 1 P
6 6 + 42.89553
I 10 125
2 10 150
3 35 225
4 35 210
5 130 325
6 125 315
In addition, the following values of the coefficientsof variation and correlationcoeffici nts are
assumedin the study:
C n ,= 0 . 1 , C b i= 0 . 1 , C p i = 0 . 1 ; (t = I,2,..., 6)
R p , p=i 1 . 0 ( i = 1 , 2 , . . . , 6 i j = 1 , 2 , . . . ,6 i i * j )
Owing to the existenceof variancesand covariancesof the randomvariablesin the formul ion of
the probllem,an analysisof the variationsin objectives Fr and Fz seemsnecessary.In study,
the covarianceof bivariate random variablesis consideredpositive or negative. Si e the
covarianr;eis representedby correlationcoefficients,the correlationcoefficientsare vari from
-1.0 to 1.0 in stepsof 0.5. one pair of randomvariablesis consideredat a time, while rest of
random 'variables are considered independent of each other (uncorrelated). By taking the ights
'tu1,
w2, and b'3 &s 0.4, 0.3, and 0.3, respectively, the percentage relative deviations in FL F2
l--' - ,. .
."r-t
Stochastic Multiobjective Generation Schedulifg
from their deterministicvalues with respectto Rn,n,Ru,n,Ro,n,R;4, srtd Rnr, Q*j) shown,
one b,y one, in Figures 6.7 and 6.8 respectively,for an expecteddemand bf ZOO w. It is
observ'edthat (i) there is an increasein the percentagerelative deviation in F, and
valueof Rnp,(i * /) is changedfrom a positivevalue to a negativevalue, (ii) there is decrease
1.5
U'
O
'o
q , ) 1
6q) r
a.
x
(.)
c
.9 0.5
(g
l<
C)
bo
E
o
o
o
0.)
tu
-0.5 0 0.5
Correlationcoefficients
---F Case1 ""4"' Case2
-'-c'-' C a s e3 --*-- C a s e4
Figure 6.7 Percentagedeviationin expectedcost with respectto correlationcoeffi
E
'A
16
E
()
!
(,)
()
()
a
5 1 4
Cg
-
1 a
€-) L /-
bo
Gl
!)
(u
O
c)
-1 -0.5 0 0.5
Correlation coeffi cient
4min = 5291.854
MW2 4** MW 2
= 5613.872
of optimal'best'solution
Deterrnination
(a) Case with dependent variables: In this case all the random variablesare nsidered
dependenton eachother and the weights w1, w2, arfldtil3 zta varied in the range 0.0 to in such
a way that their sum is 1.0. The percentagedeviation in objectives F1,F2, and F3 their
expectedminimum values are shownin Figures6.9, 6.10, and 6.11, respectively, with spectto
0.9
0.8
E 0.7
.5 0.6
.E o.s
a)
:@ 0.4
E
c)
0.3
o
L
d 0.2
0.1
16
l4
.e4 t z
q)
.E l0
' € 8
o 6
bo
(.) 4
o
fr
o
O.
a
L
0
0.1 0.2 0.3 0.4 0.5 0.6 0.7
Weight,w,
* 0.6 --^-.- 0.5 ----+---' 0.4
-.-(}'-. 0.3 --r(-- 0.2 ---.+-- 0 . 1
Figure 6.10 Percentagedeviationin expected
emissionwith respectto we hts.
3.5
J4
cA
Lr x
.= 2.5
o
C\'
O L
E
C)
bo
cl
c t s
c)
o
Lr
a)
the various combinationsof weights when the expecteddemand is 700 MW. T conflicting
nature of the objectives can be observedfrom these diagrams.The normalizedmembership
funr:tionp[, obtainedfrom membershipfunctions p(F), p(F), and p(U of objectve functions
432 Power System OBtiriization
fr, F;, and 4 respectively,is shownin Figure 6.12 for each non-inferior solution. h e ' b e s t '
solutions are given in Table 6.24 and corresponcingto the schedules obtained, th expected
power generationsare given in Table 6.25 for expecteddemandsof 500, 700, and lvlW.
0.019
() 0 . 0 1 8
-q
U)
tr 0.017
()
a.)
F
0.016
c.)
N
tr
r<
0.015
z
0.014
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Weight,w,
--{- 0.6 ---*-- 0.5 + 04
- . - . . } . _ . 0 . 3- - r ( - - 0.2 + 0 1
Figure 6.12 Variationin normalized
m embershipfunctioriwith respectto we
(b) Case with independent variables: In this case all the random variables are consi red inde-
p e n d e n r teoaf c h o t h e r R 4 4 = R a , n =R " , n = R n r , = 0 . 0 ( =
i 1 , 2 , . . . , 6 ;j = 1 , 2 , ., 6). The
'best' lr,,O=
r;olutionsso obtainedare given in Tables6.24 and,6.25 for expecteddemandsof 700,
and 9Ct0MW respectively.
Since mathematicalmodels are an idealizationof actual system models, the vari system
responsesare bound to deviate. The magnitude of such deviations can be best ev luated by
sensitirvityanalysisbecause:
(rDThe stability of the optimal solution rnay be critically dependenton es in the
model parameters.
(t') Some parametersmay be controllableand, therefore,it is importantto know t effects
may result from changingtheir values.
G) Other parametersmay be estimatedmore accuratelyif the solution is critically pendent
on one of them.
'fhe
multiobjectiveformulationis amenableto sensitivityanalysisby specifyinga sensltlvlty
index, that is, a function which indicatesthe relative size of the perturbationsin the so ution due
to var:iationsin the parameters.The specifiedindex may be includedas an additional 1ecilve to
be minimized in a multiobjectivesystem.Then the problem is solved to find the solution
[Osyc:zka and Davies, 1984]. A,second way of applying sensitivity analysisto a mul iobjective
problem is to solve the system;without Eeatingthe sensitivityindex as an objective d then to
evalua,tethe sensitivity after a preferred solution of the original system is found ig and
Haime,s,1983;Kaunasand Haimes,1985].
.,4.classicaleconomic dispatchproblem is formulated as a multiobjective optimi ion prob-
lem consideringtwo non-commensurable objectivesto be minimized, namely the rting'cost
and the impacts on the environment.The formulated problem adoptsan g-constraint , which
allows explicit trade-offs betweenobjective levels for each non-inferior solution t rig and
Haimes, 19831.The effects of random variationsin the model parametersof the opti I solution
to nonlinear programming lead to a sensitivity measurecalled dispersion (O). Th index is
interpretedas a first-order approximationto the standarddeviation in the optimal sol on of the
nonlinearprograms.A sensitivitytrade-offV"Q (AppendixD) is consideredthat gives explicit
represr3ntationof the trade-dffs between the sensitivity and the objective levels. sults are
obtainr:d for two sample systems having three and six generators,respectively.
6.6.1 MultiobjectiveOptimizationProblemFormulation
In the multiobjective optimization problem formulation, two important non-co
objectivesin an electrical thermal power systemare considered.These are economy
mental impacts.The multiobjectiveoptimizationproblem is defined as
N
Ivlinimize F{P'1= I
i=l
(aiPiz+ bi4 + c;) (6.75a)
Ivlinimize F}(P) = I
i=l
(d,Pi'*e,P,+fi) (6.7sb)
subjecil to
I e -(Po+P) =0
i=l
(6.75c)
where
A p O i, and c; 31tacost coefficientsof the ith generatinggnit
ffid fi are emissioncoefficients
d; ' €i'
Pp is the power demandto be met
as
P1 is the transmissionlosses,which are approximatedin terms of B-coefficients
N N
PL- PiBijPj
:, I
i=l j=l
that 2 is grange
where Lrz and lt are the Lagrangemultipliers. The subscript t2 denotes
function is r(P).
multipl:ierassociatedwith the secondconstraint,where the prime objective
The necessaryconditionsfor the optimal solution of Eq. (6.77) are
AL = F z ( P ) - t z = O (6.78b)
oLt,
aL D +h (6.78c)
f r = P o + ' Pd L
=l
-LP,=0
The optimal solution must satisfy the above Kuhn-Tuckerconditionsbesidesthe follow ng main
condition:
)42(Fz(P) - ez) - 0.0; )42 > 0.0 (6.78d)
The Lagrange multipliers related to the objectives, as constraints may be zeto or non- . The
set of rlon-zero Lagrange multipliers correspond to the non-inferior set of solutions. e set of
non-zerc Lagrangemultipliers representsthe set of trade-off ratios betweenthe principle bjective
and each of the constrainingobjectivesrespectively. The systemgiven by Eq. (6.77) is ved for
R values of E2 using the Newton-Raphson method. Only those values of L$ > 0 which
corespond to the activeconstraintsFi(P) = ei, r = L,2, ..., R are considered,sincethey ong to
the norr-inferiorsolution [Haimes, t977).
Irr the problem, the initial value of s2 is taken such that tz ) Fjtn and E2 1 F2 . Since
cost and emission are of conflicting nature,the value of objective Fy will be maximum, hen the
valueof F, obiectiveis minimum and vice versa.So, minimum and maximum valuesof 1 and F2
are obtained by performing economic dispatchand minimum emission dispatch separa ly.
Similarly, for more than two objectives,the trade-off functions LU Q + il, can be enerated
-
with respect to the fth prime objective, by varying e; while other e* (k /) are fixed at some
level. l'he ith prime objective can be replacedby the yth objective and the solution rocedure
is repeated for more information. Algorithm 6.2 can be implemented for the ge tion of
non-infbrior solution.
6.6.3 ParameterSensitivity
In practice, the values of the objective and the constraintsof the problem defined by . (6.76)
will actually depend on the decision variable P;. Usually each equation is formulated in terms of
the decision variablesand other parameters.Therefore,the values of the objective and straint
functions depend not only on the values of the decision variables but also on thc alues of
coefficirentsdi, br, ci, dr, €i,,f, and Bii (i ='1, 2, or these
coeffic:ientsare approximatedby estimation, statisticalaveraging,past experience,and there-
fore pr,oneto error.
Irr the light of the above consideration,-themultiobjeCtiveoptimizattonproblem as given by
the prc,blemdefined by Eqs.(6.76) is reformulatedby assumingthat, in addition to decision
variables Pi(i = I, 2, ..., N), the objective function is dependenton a set of'cost x;(i =
1,2,..., 3M) and the constrainedobjectivealong with the problem constrain.r,s is depen t o n a
technology parameter set yi! = I,2, ..., m).
I\4inimize lFr(P,x)l (6.79a)
subject to \(P,y)sei U=2) (6.7eb)
4(!, Y)= u (6.79c)
.d
Stoclmstic Multiobjective'Generatiort Schedulilry
where
ct,..,cNlT
lQ-..,ay, b1,..,b7r1,
v- fr,..rftu, Btt,
ld1,..dy, €1,..,€p1, 812,.-.,Bruru]'
N
q ( P , ) ) I= a . - P L (6.7ed)
i=l
lerlu)').1lhe dependenceof the optimal solutionF* to the problemof Eqs. (6.80)on the pa eter
setsx, y, z ca.nbe expressedby writing
F* = F* (*, y, z) = F* (w, z) (6.81)
wherew = f, * I y lr"
The parametersetsr and y are the approximatedvaluesof the actual unknown
Keeping this in view, x and y can be treatedas random variables.Therefore,wi dte
random variableswith meansand variancesas
definedby Eqs. (6.80) can be consideredas a function of the random variableW ; , 1 . 9 . f = F. (w)-
By ignoring higher terms, Taylor seriesexpansionabout the mean is
var(F*)= o! (6.84)
l[#)'
where thrl partial derivatives are evaluated at the point w i = wi
Power System Optimization
O- (6.8s)
ao= [ ar".)l f
_ tT ^?(aF.')
a, l*"tlm;JtrrilJj
d (6.86a)
where
V*F* = LL* Vrg I VrFl (6.86b)
(6.86c)
a4
a4 = -
dti
where P *o (6.87)
acl aC, dtj
dti
As the numberof objectivesincreases, the trade-offsof sensitivitywith respectto obj, tive levels,
ti, ire determined.Thesetrade-offsare the componentsof the vector z;, which are cal ulatedfrom
Eq. (6.86a).The calculatedtrade-offsare submittedto the decisionmaker to select preferred
solution.
V r F *= Y p F Y r P * (6.e2)
Next difflerentiating
the expressiong(P*, !) = z yields
VpgVyP*+Vr8=0 (6.e3)
Multiply Eq. (6.89a)by VuP*to obtain
YpF YrP*+ L.Yrg VrP* = 0 (6.e4)
MultiplyEq. (6.93)by t to obtain
L. YrF YrP*+ t Vrg = o (6.es)
After solving Eqs. (6.94) and (6.95), the following can be obtained
YpF YrP* - A) Yyg (6.e6)
EquatingEqs. (6.96) and Eq. (6.92), gives
V;F*=YpF%P-+VrF
Power SYstemOPtimization
problem efined b,y
Kuhn-Tucker necessary conditions for an optimal solution to the
(6.88) are
YpF+1,V6=9 (6.100)
8(P) = 0 (6.101)
(6.107)
The implicit function theorem assuresthat in some neighbourhoodof zero, there ex a unique
L* - h*(z),P* = P.(z)to thesystemEqs.in (6.100)and(6.101).Moreover, *(z)
solurfion exists a
neighbourhood the functionP-(z) and
of zero suchthat for any z in neighbourhood satisfy
the second order sufficiency conditionsfor a strict local minimum.
So V.P*and YrAlcan be expressedas
[o.".-l_l v r r (vps)rl[ol
- Lvrrr (6.108)
Lv.r.J o llrl
PurringEq.(6.106) in Eq. (6.107),the requiredexpressionis obtained.
Stochastic Multiobjective Generation Scheduling 441
F r r = 0 . 0 1 0 P+? 2 . 0 P 1
+ 10.0
F r z =0 . 0 1 2 P+1 1 . 5 P+2 2 0 . 0
Ftt= 0.00aPl+ 1.8P3 + 20.0
NtC,,emission characteristics(kg/h) are:
F z t = 0 . 0 0 0 6 5 7 0 -P0?. 0 5 4 9 7 P+1 4 . 1 1 1
Fzz= 0.00059 I6P3- 0.05880P2 + 2.593
Fzt= 0.0004906P? - 0.05014h+ 4.268
The B-coefficients of transmissionlosses(MW-l) are:
Generationschedules(MW)
Demand- Sr.
(Mw) no. Pl P2 P3
Generationschedules(MW)
Demand Sr.
(Mw) no. Pl P2
6 72.337 76.175 63 7 r
7 72.591 76.t35 63 37
8 72.836 76.094 l8
9 70.576 84.195 53 29
250.0 I 92.289 99.35r 75 75
2 96.169 94.r76 77 40
3 96.077 94.456 76 52
4 96.025 94.723 76 l 6
5 96.803 92.978 78 25
6 97.251 92.792 77 32
7 97.804 92.894 76 39
of results
Table 6.29 Comparison
Demand Dispatch Fl F2
(Mw) ($n; (kg/tt)
Dem"and Pl P2 P3 P4 Ps
Minimumcost
200.0 29.886 r5.560 8.193 55.416 28.222 6 5 . 18 6
400.0 78.397 40.818 2r.732 133.535 45.451 88.959
600.0 12i.878 66.914 36.060 2t2.450 63.M9 113,216
IMinimumemission
200.0 30.638 37.356 35.994 22.T33 36.194 39.514
400.0 57.015 66.398 7r.747 68.100 71.222 72.634
600.0 84.025 95.916 108.806 rr4.795 106.622 105.969
dispatch
Combinred
200.0 31.863 21.725 12.158 45.M5 30.2r5 60.856
400.0 63.349 5r.969 31,.384 120.118 52.952 88.463
600.0 rt2.823 79.026 48.674 187.398 74.r79 116.628
6.7 ISTOCHASTICMULTIOBJECTIVESHORT.TERM
IIIYDROTHERMAL
SCHEDULING
The sr:hedulingproblem requiresan appropriateobjective function such as minimum neration
cost or transmission lossesor minimum pollutionlevel, etc. In essence, the hydrothermalcheduling
problem can be visualized as a multiobjectiveone. A simple and robust solution me ology for
a class of multiobjectiveshort-rangefixed-headhydrothermalproblems is undertakn which
makes use of fuzzy set theory.First of all, a short-rangefixed-headhydrothermals uling
problem is formulatedin a unified multiobjectiveframework consideringstochastic st, NO,
emission, SO2 emission,and COz emissioncurves for thermal power generation nits and
uncertaintyin systemload demand.The expectedvaluesof thermalt'uelcosts,NO, emis o n s , S O 2
emissionand CO2 emission,over whole of the planningperiod are the four conflicting bjectives
to be rninimized.Further,the expectationof the squareof the unsatisfiedload becauseo possible
varian,ce of generator outputs over whole of the planning period, is incorporated another
objective to be minimized. Basically, the solution procedure for the multiobjective prob is based
on the generationof non-inferiorsolutions.The weighting techniqueis used to gen ate non-
inferior solutions,which allows explicit trade-offbetweenobjectivelevels.Exploiting fuzzy set
theory [18], a cardinalpriority rankingof the non-inferiorsolutionsis definedthat maxi izes the
satisfactionof all the objectivesand is utilized to find the best compromisingsolution from the
non-inferior solution set. The practical viability of the multiobjective hydrothermal heduling
proble:mhas been demonstrtedon three samplesystems.
6.7.1 StochasticMultiobjectiveOptimizationProblemFormulation
Consideran electric power systemnetwork having N thermal generatingplants and M hydro
plants,where M + N is the total numberof generatingplants. The basic problem is find the
activepower generationof eachplant in the systemas a function of time over a finite ti period
from 0tto T.
Stochrasticthermal model
The objective function to be minimized is the total systemoperatingcost, representedb the fuel
cost of thermal generation,over the optimizattoninterval.
possible deviation in coefficients of input-output characteristics and load demand fro their
,*p""t"d values are manipulated through the randomness of random variable, P;. A str tic
model is reducedto its deterministicequivalentby taking its expectedvalue. Presumin that
Jr = a'
+d,+ a; vartnrl)
@,1'+ biFi . rl 0 )
I [*
where
Qi, ,bp and di, are expected cost coefficients of the ith thermalunit
p
t I
jrsthe expectedpower gcnerationby the ith thermalunit.
rt = (d,,P,, +rzi)]d'
+e,,F, 6.1r4)
I [*
where d.2i,€2;,nd fzi are SO2emissioncoefficientsof the ith thermal unit.
u8 Power System OPtimization
T(N ,,
. )
rq= J > @2,42
+ ezi4+ fti) ( 6 . 1l s )
| la,
o \i=t )
wherc dti, (ti, andft;, are CO2 emissioncoefficientsof the ith thermal unit.
Development of stochasticmodels of functions Jr, ,I3, and J4 are adopted after much
cogitation that thermal generationsand load demandduring each subintervalare variables.
Any possible deviation of NO' SOz, and CO2 emission coefficients and load deman from their
expected values is managed through the random power generation P;. Presuming t random
varierblesare norrnally distributed and statistically independent, the expected value of N r emlsslon
is esrtimatedas
J2 = + aLiFi
@,,F,' + Iti+ d,, var(P,ll]a (6.116a)
I [*
wherc d1;,V1;,and fr; are expectedNO.. emissioncoefficientsof the ith thermalunit
Substituting Eq. (6.111)into Eq. (6.116a),
T(u )
i, = J- o | > [(r* cl) duFtz
+a,,F,
+frilla, (6.116b)
\i=t )
The expectedvalue of SO2emissionbecomes
it = (d,,F,'
+d,,8+Iz,+d, var(piD)" (6.117a)
I [X
where d2;, e2i, and j2,, are expectedSO2 emis5ioncoefficientsof the lth thermal uni
Equation (6.Ll7a) can be rewrittenas
i4= +voF,.,,,])
* cA)i,,F|
[rr (6.1r 8b)
I [l
Stochastic Multiobjective Generation Scheduli
Q j= K Q ( h )r ( P ) ( / = 1 , 2 , - . - ,l A (6.11e)
where
Q, T are functions of head and hydro generations,respectively
K' is a constant
Izf is total number of hydro units.
For a large capacity reservoir,it is practical to assumethat the effective head is constant
over th,eoptimizationinterval.In the caseof constanthead, fih) becomesconstant,Eq. . 1 1 9 i)s
rewritternas
Qj = K' r(P) ( j = I , 2 , . . .M
, ) (6.r20)
where K' becomesthe new constantand is formed by the multiplication of K and ). Each
hydro plant is constrained by the amount of water available for the optimization interv i.e.
J
0
Q 1d t = R i (j = L,2, ...,Itf) (6.r2r)
Q i = x 1 P , 2 + Y+
f 12 1 ( / = L , 2 , . . . 1, 4 ) (6.r22)
where xi, !i, and zi are the dischargeccefficientsof the 7th hydro plant.
Slincethe thermal generationsand load demand are random, the hydro gene ons also
become random in view of the load demand constraint as given by Eq. (6.124). A tochastic
model of function Qi is developedby deeminghydro generationand dischargecoeffici ts during
each subintervalas iandom variables.Any possibledeviation in dischargeand load de and from
the expected value is manipulatedthrough the randomnessof hydro generator Pj. ince the
randonnvariables are assumednormally distributed and independent,the expected value of
dischar:gebecomes
( / = 1 , 2 ,- . .M
, ) (6.123a)
4 i = T t F l+ y i F i + 2 , + 7 , v w ( P )
where ii, yj, and Zj are the expecteddischargecoefficientsof the 7th hydro plant.
Eq. (6.1'23a)can be rewrittenas
4 i = . Q + C f ; ) i i F+ty' r F , + 2 , ( / = 1 , 2 , . .W
., (6.r23b)
45tD Power System OPtimization
Equalityand inequalityconstraints
(i) The expected load demand equality constraint is
M+N
sr
= (6.r24)
L P' Po+Pt
i=I
"'"
,ois the expected load demand during the interval
P--min
<Fi <P.max (i = 1, 2, ..., N + IuI) (6.12s)
where
^tn interval.
P, is the expectedlower limit of the. ith generatoroutput during
Prt* is the expectedupper limit of the ith generatoroutput during interval.
Expectedtransmissionlosses
A ,coffirrlonapproach to model transmissionlosses in the system is to use the Kron's pproximated
loss formula through B-coefficients
gg
Pr= L PiBUPj (6.126)
L
i= 1 j=L
or
N+M N+M
(6.12e)
Pt=
I \F,ruF,
j=l
i=l
Stochastic Miltiobjective Generation Schedul 451
Expecteddeviations
Generatoroutputs P; are treatedas random variables,and the stochasticmodel is co erted into
its deterministicequivalentby taking its expectedvalue. So, the solutionwill provi only the
expectedvaluesof power generations. By virtue of the aboveconsideration,
therewill mismatch
in load demand. The variance of a random variable quantifies the degree of uncertaint associated
with the mean value of the random variable.The active power loss, the systemfu I cost and
emissioncurves are quadraticfunctionsof decision variable P;, and their variancesq ntify the
degree of uncertainty associated with their expected values. So, the expected mi h can be
estimiltedthrough minimization of the squarederror of the unsatisfiedpower demand, l e
E- (6.130)
where P; is the actual power generationrequiredto meet the load, which is consi random.
Using Eq. (6.124),Eq. (6.130) can be rewritten as
E_ (6.131)
M+N-I M+N
1-
") var (P,) + (6.r32)
i=| j=i+l
I _ (6.133)
"5
where;
S;i= e,
SU= Rr,r,Cr,Cr, (i #j)
452 Power SYstemOptimization
optimizationformulation
Multiorbjective
(a) expectedoperatingc t, (b) the
Multiolbjective optimtzationproblem is framed considering: the
NO, ernission,(c) the expectedSO2emission,(d) the expectedCO2 emission thermal
"*p""t.ri variables P; from
unitr, and (e) the expectedrisk associatedwith possibledeviationof random
loar demand
their respectiveexpectedvaluesover the optimizationinterval to meet the expected
in eactr interval. Each hydro plant is constrained by expected amount of water available or draw-
down 1n the interval. Mathematically, the multiobjective optimization problem is defin d a s
M
t -
I q' J , d t - RJ ; ( i = 1 , 2 , . . . ,l v l ) (6.13ac)
J
0
T
lylinimize J4 T
.L
'- [r,* c?,) * v,,F,*,f
d,,P,' (6.135d)
k=I [*
lvlinimize Js= (6.135e)
M+N
(6.13s0
subject to
Lj = l F,o= Fp1,-F11
I
k=l
t*Tir = Rj (/ = 1,2, ..., It4) (6.1359)
p_?t"sF,r<F,.flu* (6.135h)
Stochastic Multiobjective Generation Scheduli
6.7.2 SolutionProcedure
To generate
the non-interior optimizationproblem,the weighting
solutionsto the multiobjective
methodis applied,In this method,the problemis converted into a scalaroptimizatio as given
below':
Minimize (6,136a)
subjectto (6.136b)
T
( j = 1 , 2 ,. . . ,I u [ ) (6.136c)
L ' 0 4 , 0 =Rj
k=l
5 , t M
L- Z *oio*) +F'o
viair+Lr'[u'- I,,4 (6.137)
k=l k=l [,f, Zal ,/=N+l
where:
vr.is water conversionfactor of the 7th plant
/4 is incrementalcost of power deliveredin the system.
The <iptimality conditidns are describedby taking the partial derivatives ugmented
objectivefunction with respectto the decisionvariables.
i , , % * L , f % - t )^ ) =o ( i=r, 2 , . N
. . ), (6.138a)
A,,rdF,o,,"*[a";
I
k=L
ttQ i * -R j = Q ( J = 1 , 2 ,. . . M
, (6.138c)
Optimization
For+r*-ZF,o-o (6.138d)
i=l
'rffi+Looffi')
ouu.[H
(6.139a)
a2F,
aP*ka
y[&.-,) (-
= -l %o+F&
N+M \
(6.139c)
?r 1EP,* )*tr i=l )
zruFf
5
n i - . f
F i = w 1 b , *L * i t i t
j=2
( i - 1 , 2 ,. ,
To simptify calculat{ons,neglect all the terms with S,i and Tij, i * i. So, neglectig the term,
N+M
t( M+N
(2toa,+ ztowrf Ti) lrik = -l tolza,rfl+ Fi I ws
i, + 2fuor ( ik
L \. j=l
1,
LPu ( i = L , 2 ,. . . ,M
or
LF,* (6.140)
where
A4= t{Au,Plr+gt+ Li*) (6.1a0a)
rpsSi,+ trff,,1
X4= 2(t11ai+ (6.140b)
J=T
456 Power SYstem OPtirnization
Ti= 1t+Cr;*,
M+N
K|o=| zr,,,rfl
l=l
n) LF*r
rr + 2vlt oC'r,ii + zLDoTn
(2t1,w5S
( / = 1 ,2 , . . . , M i * = j + N ; ft=1,
or
-L^r,
^n 0-KI*; 1;"* -t*(ZTiF,9,o+!)v]"*
Af ^1=
or
e - K?,i h";* - E1,v]"* - L^k (6.141)
M^k
Yio
where
Ej* = tr(2TiFlo+ y,) 6.l{la)
6.141b)
vt
jk - 1rc7,ii + hooT^*)
5s*^ + voit
2(t1rw
Stochastic Multiobjective Generation Schedul
N+ M (t' t+ u \ ( N+M \
=-[";**F:o
zriErt-r)M'r Ff ( f t= 1, 2 , . . . ,D
E t; ] )
N+M N+M
- =
.I tt KlolLFp For,* F&-I &-' ( k - 1 , 2 ,. . . ,D
,l=l i=l
i ,t - K,oo,
r =
tt - I - o-
[rT" ",g
, \ x i k
'1
)
.X,,-"n,I L"fn (l - Kl"i - E jkvlew-
Yio
( N + M \
m=j+N)
co L"f* - I
M
Di1,v]* - Fr (6.r42)
j=l
where
(r - Kl,r)z (6.r42a)
C*=
Yio
Dj* = (6.r42b)
( 1- , K , ? ) A i r
Ft= (6.r42c)
x,*
System Optinization
( 6 .l 3 9 d ) ,
Substituting the values
(/=
or
T
(6.r43)
T Dir,L"fn-Hi']"*
ik=l
=oi (i=
where
T F n
Lr --
,'j S "ik"ik 6.143a)
.L Y,,.
k=l J'r
oj=vj-i,ralo+iry 'ik
6.143b)
k=l k=l
K = +.$
(6.144)
" Ly* ?yle*
ck'?rck'l
D i * Dyi x1 . * l - " r r t =
eo* j - Z + ( i = r , 2 , . .M. ,; * = i + M (6.145)
c k ' )
Qii =
E ( j = I,2, "', Il'l) (6.145a)
Stochastic Multiobjective Generation Schedulin
e' i tL= ft W
r
( j = t , 2 , . . . , ,IM
= ;1 , 2 ,. . . ,M ; j * l) (6.145b)
tv
k=l k
T
R j =o j - > T ( j = r , 2 , "w' , (6Ja5c)
k=l
Here only matrix of M x M is to be solved to calculatevjnt*, and A,t'* can be com from
Eq. (6.144).When the valuesof url'* and 21"* .areknown then n4olt - 1, 2,..., M) a LF*
(/ = 1,,2, ..., M; m = / + /f) can be computedtiom Eqs. (6.140)and (6.141),respectiely. The
detailed algorithm iiss elaborated below:
=*.i+,;"*
L+"* tvk
j=l ^
; if Pf"* t P#u*
tDnew
ik _
; if P;f"*< P,.fln
; otherwise
460 Power SYstemOPtimization
whose limits have been set either to lower or
Disallow to participatethe generation'
to zero.
upperlirnitsbysettingtherelatingcoefficients
15. i , > IT, then GOTO SteP5,
e l s er = r * 1 ,
;0
rik r[-new
ik
= Llr"*
Lo* (k = 1, 2, ..., D
vf = Vjn"* (J = 1,2, "', ll)
pB= K i=1 5
T pO!) (6.r47)
I
k=l
\ u<i!'t
i=l
The function 1tp in Eq. (6.147) can be treatedas a membershipfunction for noR inated
solutions in a fuzzy set and representedas fuzzy cardinal priority ranking of the non minated
solutions.The solutionthat attainsthe maximummembershippB, in the fuzzy set so ed can
be chosenas the best solution or the one having the highestcardinal priority ranking
Maxtp* , k = L,2, ...,Kl (6.148)
of thermalplant ($ft):
Expected,costcharacteristic
/r r(r) = 0.001991
4' + 9.606Fr@ + 373.7
Expected .NO, emission characteristic of thermal plant (kg/h):
jrr(t) = 0.006483 - 0.7902t
4' ?ft> + 28.82488
of thermalplant(kgA):
ExpectedSO2emissioncharacteristic
jrr(t) = 0.00232Pr'+ 3.84632Pr(t)
+ 182.2605
of thermalplant (kg/h):
ExpectedCOz emissioncharacteristic
jor?) = 0.084025F12
- 2.9445484
Pr(r)+ r37.7043
(Mm3ztr):
Expectedhydroplant characteristics
(MW-l):
ExpectedB-coefficients
If the coefficient of variation is zero and randomvariablesare uncorrelatedto each othe , then the
problem is considereddeterministic.The random variablesare uncorrelatedonly if the orrelation
coefficient is zero.
5
4.5 t
4
c
3.5
/ /
r ./
/ /
3
0.) .J .I-
o
bo
2.5 a /
. . ] t ' '
By st:tting all weights equal to 0.20, the effect of varianceon all objectivesis :ved
efVet
simultaneot rsly- The percentagerelative deviationsin total expectedcost, NO" emiss.on, S Soz
O
emissionarrd CO2 emissionfrom their respectivedeterministicvalues, with respectto co :ffi ents
cient
of variation are shown in Figure 6.13. It is observedfrom Figure 6.14 that the relative p )rce;ntag tage
on cost bec )ome
mes
deviationsof all the objectivesincreaseas the varianceincreases.The effect
smallerconrparedto NO' SO2and COz emissionswhen equal importanceis given to all c bje< :tivet
tves.
Further, the: effect of varianceon water conversionfactor is also considerable.
The r:ovarianceof bivariaterandom variablescan be consideredpositive or negalve The
covarianceis representedby correlationcoefficients.The correlation coefficients are va ied fron from
-1 .0 tc l.t) in steps of 0.2. The percentagerelative deviations in total expected c NO
NOr
ernission,SiO2 emissionand COzemissionfrom their respectivedeterministicvalueswil r()spec ;pect
to correlation coefficient(Rr,r, (i f i) are shownin Figure 6.15. The weights wy w2, w2 w 4r Bllt and
w5 are takr)n as 0.25,0.25,'(i.25,0.25 and 0.0 respectively.It is examinedthat (i) th re is rs an
a,
increasein the percentagerelative deviationsin total expectedcost Jr tr the value of R, Pj i i ++j t)
is changedfrom a negativevalug to a positive u+t, (ii) there is a decreasein the p :ntag
rtage
relative de'riations in total expectedNO, emission Jz as the value of Rnp,Q/ il is charse! fror from
a negativevalue to a positive value, and (iii) there is a very small effect on SO2emissio:r . /3r, ZlIl
and
Stochastic Multiobjective Generation Scheduli
15
t4
t3
t2
11
q
g
l0
Cg
9
o)
!
()
8
@
6J 7
O
(-) 6
(.) 5
4
3
2
1 -:l=-fi:#+-:+-+:1:
0.05 0.07 0.09 0.11 0.13 0.15 0.17
Coefficient of variation
COz emission J 4 objectives.From this study,it has been observedthat the existence f random
variablesgives a significanteffect to each objective either consideredindividually or in the
multiobjectiveframework.
1.5
tr
'E
o
r.zs
'5
a)
o l
bo
cd
d
$ o.zs
0)
or
7y" = 14591.82
kg Jfl'* = 15169.05
kg
tj-" = 44337.36
kg /j"o = 44665.53kg
J;'" = 245810.70
kg J fu* - 252472.20
kg
j{ * = 46805.98Mw2
"ry" = 46044.05
lvrw'
Sr w1 Jl J2 t t ( Jr ) tt(Jz)
no. ($) (ke)
1 1.0 0.0 96028.66 14839.02 1.0000 0.0000 0.07057
2 0.9 0.r 96032.16 r4774.06 0.9966 0 . 11 2 5 0.07827
3 0.8 0.2 96044.09 r4707.rr 0.9851 0.2284 0.08564
4 0.7 0.3 96067.08 14638.60 0.9628 0.3471, 0.09244
5 0.6 0.4 96104.& 14569.26 0.9265 0.467r 0.09835
6 0.5 0.5 9616r.45 14500.20 0.8716 0.5867 0.10292
7 0.4 0.6 96243.90 lM33.l3 0.79r9 0.7029 0.10549
8 0.3 0.7 96360.7r 14370,68 0.6789 0.8110 0 . 1 0 5l 5
9 0.2 0.8 96523.80 r4316.77 0.5213 0.9044 0.10061
10 0.1 0.9 96750.05 14277.39 0.3026 o.9726 0.08999
11 0.0 1.0 97063.07 14261.55 0.0000 1.0000 0.07057
3 0.6
E
a
v)
q)
2 0.4
c)
-1
conflict or vice-versa.Figures6.L7 and 6.18 show that the expectedcost is in conflict with the
expectedSO2 emissionand risk level, respectively.
Power System OPtimization
0 .9 8 5
0.97
(J
0 .9 5 5
.a
() 0 .9 4
L
C)
0.925
0.91
o.8esd 0.1 0.2 0.3 0.4 0.5 0.6 o.1 0.8 0.9
Weight,w,
o
3 0.6
(a
c)
€ 0.4
q)
ii.-. .
stochastic Multiobjective Generation scheduling 467
The non-inferiorsolution that attainsthe maximum membershippB, is
distinguishecas the
best solution among the non-inferior solutions.The weight combination,presented
a serial
number7 in Table 6.33 gives the maximum valueof i.e. 0. I0504g and thereforeprovi
lt\, the
best or preferredweight combination.
The non-inferior solutions for 126 differentsimulatedweight combinationsare ge
consideringall, the objectivessimultaneously.Non-inferior solution that acquiresthe
ma lmum
membershippj, is chosenas the best solutionand is furnishedin Thble6.34. Ttre best
lutions
are securedfor distinct values of coefficientsof variationand correlationcoefficientsin lverse
situationsand are conferredin Table 6.34.
The weight combinations and the water conversion factor v2, cotrssponding to best sch ules
are depicted in Thble 6.36. For case one, the achieved expected generation schedules of 24 hours
have been furnishedin Table 6.35. For each sub-interval,expeotedtransmissionloss FL, incre-
mental cost L, and expecteddischarge qz, are exhibitedin Table 6.35. The attained uality
constraintLPD,during each sub-interval,is also providedin Table6.35 which showsthe
of the obtainedsolutions.
Conventional economic short-term fixed-head hydrothermal power dispatch method al tes
generationscheduleto the individual generatingunits basedupon deterministiccost functi n and
load demand,ignoring inaccuraciesand uncertainties.Such generationschedulesresult n the
lowest expectedtotal cost, but this cost is also associatedwith a relatively large varianceth can
be interpretedas risk measure.Moreover,in power systemoperationplanning, there exist m Itiple
objectives to be attained, which conflict with each other and ard subject to a mutual interf It
means that any one objective can be improved only at the expenseof other objectives. n the
multiobjective framework, the analysisof hydrothermalshort-rangefixed-headis undertakr with
explicit recognition of uncertaintiesin productioncost, NO' SO2 and COz emissionsa load
demand.
6.8 STOCHASTIC
MULTIOBJECTIVE
LONG.TERIJI
HYDROTHERMAL
SCHEDULING
A modern power system may consist of several thermal, conventional hydro power lants
l
connected to various load centres through a lossy transmission network With the insigni cant
incrementalcost involved in hydro generation,the problem of minim izing the operational of
a hydrothermalsystemcan be reducedessentiallyto that of minimizing the iuel cost for al
Power SystemOptimization
whereA,Pp= (4 + F ) - Ft-F,
Table 6.36 Weight combination and water conversion factor corresponding to the sc edule
given in Table 6.34
plants under the constraintsof the water availablefor hydro generationin a planned period.
Mostly, hydrothermaloptimal schedulingis achieved,with .the asiumption,that the water nflows
to the reservoirsand the load demandsare known with completecertainty.However,thi is not
true.
The availability of limited amount of hydroelectric energy, as stored water in the system
reservoirs,
makesthe optimaloperationcomplex,because
it createsa link betweenan rating
decision in a given stage and the future consequencesof this decision. Further, it is im ible to
have perfect forecasts of the future inflow sequenceand the load variation during a given riod.
Therefore, for long-term storage regulation, it becomes necessaryto account for the rando nature
of the load and river inflow and so a stochastic representationof these must be used.
Most of the algorithmsincorporateuncertaintiesin the system load demand an water
inflows, but choosea deterministiccost function for thermalgeneratingunits. A major s rce of
uncertaintyin optimal dispatchis that associatedwith cost coefficients[Dhillon er al., 19931.
Howeverwith the increasing,concernrecentlygiven to the environmentalconsiderations illon
et al., 1993; Dhillon et al., 19941,a revised generation scheduling for the hydrothermal power
systemis requiredthat meetsthe constraintsof availablewater at hydro plants and load emand
for power while accountingfor both cost and NO, emission.
Fuzzy sets were first introducedin solving power systemlong-rangedecisionmakin prob-
lems. Fuzzy decisionmaking theoriesattemptto deal with the vaguenessor fuzzinessinh nt in
subjective or imprecise determinations of preferences,constraints, and goals. Thpia and urtagh
[1991] put up a methodologyfor solving a decisionmaking problem involving a multipli ity of
objectivesand selectioncriteria for the best compromisedsolution.
The intent of this section is to provide a technique that allows scheduling of lon -range
hydrotherrnalsystemprobabilisticallyconsideringstochasticcost and NO, emissioncur es for
thermal power generationunits and uncertaintyin load demand and reservoir water i flows.
However,there is a growing trend towardsformulatinga multiobjectiveoptimizationp lem
[El-Hawary and Ravindranath, 1991], so, the approach is developed by formuiating hydrot ermal
scheduling as multiobjective optimization problem. The expected fuel costs and NO, e i s s i o n
over whole of the planning period are consideredas two conflicting objectives.The form lation
also incorporatesany possibledeviationsin generationsover whole of the planning period as the
third objectiveto be minimized.The weightedminimax techniquelKlir and Folger, 1993] i used
to generatethe non-inferiorset by convertingthe problem into a scalaroptimization To
reduce the complexity of the problem, interval-wisedecompositionis carried out. Eac sub-
problem is separatelysolved by using the conjugategradient method to obtain the timal
discharge[Parti, 1987]. In each subinterval, thermal generationsare calculatedby a sim lified
technique,which reducesthe economicdispatchproblem into an equivalentlosslessproble . The
methodis lesstime consuming.A numericalexampleof a power systemconsistingof th hydro
and four-thermal plants is solved and the results are presented.
6.8.1 StochasticMultiobjectiveOptimizationProblemFormulation
In this section,the multiobjectiveswith equalityand inequalityconstraintsconcerningthe power
systemoptimizationproblem are described.The importantobjectivesare consideredhere, e
1. Economicoperations
2. Minimal impacts on environment
3. Expecteddeviationsdue to unsatisfiedloads.
470 Power System Optimization
Fr = (ai(Pi^)z
+ b,P,^*',,], (6.14e)
E []
where a;, b;, and c; &ra cost coefficients. P;' is the thermal power generationduring he mth
subinterval.
A stochasticmodel of function F1 during the mth subintervalis formulatedby co idering
cost coefficientsand load demand,during the nth subintervalas random variables.The xpected
value of fuel cost function rnay be obtainedthroughexpandingthe function using Tayl s series
about the mean. The obtainedexpectedfuel cost during the nth subintervalis represen by
N
Fr"= I la,1F,\' +6,P,^*V, + a, vw (P,^)+2P,^ cov(ai,Pi^)+ cov(b,,P,^)] Sltr (6.150)
i=l
where
Pi^ is the expectedvalue of thermalgeneratoroutput during the nth subinterval.
d,, 6;, and d, *. the expectedcost coefficients.
Equation(6.n50)can be rewrittenas
1V
Fr^= lei (F,^)' + Bi F,^* q ] slrt (6.151)
j=l
where
At' = [1.0 + (C(Pi\)2 + 2R(ai,Ph C(a) C(Pi\]Ai
Bi" = [ .0 + R(bi, Pi) C(bi) C(Pi) bi
C(Pi^), C(a), and C(b) are the coefficientsof variation of random variables P{, ai and br
respectively.R(a;, Pin')is the correlationcoefficientof randomvariablesai and P;^, and b i , P i )
is the correlationcoefficientof randomvariablesb; and P;''.
ExpectedNOxemission
Only thermalpowerstationsaremajorcausesof atmospheric of high oncent-
pollutionbecause
ration of pollutants caused by them. -{-heNO, emission curve for thermal power plan can be
Stochastic Multiobjective Generation Scheduling 47r
directly relatedto the cost curve through the emissionrate per MBtu, which is a consta factor
for a given type of fuel. The aim is to optimize the NO, emissionof thermal stations ith full
use of water availableduring the optimizationperiod.The amountof NO, emissionis gi n a s a
function of the generatoroutput P;^, which is quadratic.
Fz- f ti (di(p;*)z
+eip,^*
rll re (6.rs2)
;L,J I
where di, ei, ffid fi are emission coefficients.
A stochastic model is formulated by considering emission coefficients and load and as
random.Using Thylor's seriesand taking expectations,
the expectedNO, emissionfor mth
subintervalis obtainedas
N
- L
Fi - >Wt(1^f+a,F,^
+ fi + cl, var(Pi*)+21^ cov(d,,p,*)+ cov (e;,pi^;] tgltr (6.153)
i=l
where
Expected deviations
Since generatoroutputs P!' of hydro and thermal plants are treated as random variabos, the
expecteddeviationsare proportionalto the expectationof the squareof unsatisfiedload mand
during the nth subinterval.These expecteddeviationson the whole of the planning pe od are
conSideredas an objectiveto be minimized.The expecteddeviationsduring the rzth su nterval
are representedas
MW2 6.1ss)
where
PI is the expectedload demandduring the ruth subinterval
F; are the expectecitransmissionlossesduring the mth subinterval
Z i s the total number of hydro and thermalplants.
472- Power SYstemOPtimization
This on simPlificationreducesto:
T T T
u* (Pi^) (6.1s6)
F;'= i
i= I i=l j=l
j*i
where
'T"
;::;:ffT:'L
:i:ff#:::i ,ffi T:,'l:il:'
11,'::ffi
The corresPondingvariancesare:
- EKXf*t - x;*')'l
var(xl"*t) U = 1,2,..,
L)
Stochastic Multiobjective Generation Scheduling 473
= EIQfrw-F,T*)'l
var(Pfrls) (i = 1,2,...,L)
On simplification,the aboveequationcan be written as
M M
Minimize ,I Fy,\ (6. 67a)
m=l m=I
-474 Power System OPti^rzotign
M M M
ytn+l -x|-Lti' * I e | ' * I t i ' = o (j = 1 , 2 ," ' , L ) (6.l67d)
m=l nt=l m=l
( i = I , 2 , - . . ,M (6.167f)
4 m i n{ F , u ' < 4 t t *
Minimize (6.168a)
:[]
3
subjectto (6.l68b)
L*o-t
ft=l
L-ili,*ru['+Ll'
.L,l I Lt
nr=l L t=l
L
- I tii [Fji* -t,j (r + o.ssi eV|" *ij" - Q';'- s;r'>)<Q'j'- ! 1))
j=l
During each subinterval,the control variablesare the water dischargesthrough the turbi es of
hydro plants.The reservoirstorageand the hydro power generationat the end of the subi terval
are obtainedfrom Eqs.(6.160)and (6.163),respectively. Irrespective of the hydro generatr s the
thermal generationssatisfythe power transferconstraintto achievethe minimal fuel c t . T h e
dual variablesLf , ryj, and ).fi are obtainedby equatingthe partial derivativesof the I-a gian
function with respectto the dependentvariablesto zero.
* - # + r ; ( # - t ) = o ( i =r , 2 , 1 r ) r70)
#=l
= o ( i = t ' 2 '' L ) 6 . t 7t )
#--tii+fi[+-')
t\;' - h\i + Mi hi siQi' - t-r) = o 6.n2)
aor\ \'lPfirn,pillr,)l
[*(u*\- ei]il)*I
*wt'LaQi' z*(cou . 4'(*\
',!.rji,)r,...,,,
k=t J
6.n3)
where
The dual variablesLai are to be adjustedto maximize the Lagrangianfunction under constraints
of other optimality conditions.The correspondinggradientvector is
. M M M
AL = x f * t- x j - I
6L ai
i f * Z O f* I s . f ( i = r , 2 , . .L. , (6.t77)
m=l m=l m=l
The upper and lower limits on the control variables are taken care of by making se variables
equal to the respective bounded values whenever such limits are violated. For he dependent
variables, these limits can be considered by augmenting the cost function throug the Powell's
Penalty function.
L** ( i = 1 , 2 ,. . . ,N
k=l
The above set of equations is solved by a simplified method [Osyczka and avies,1984]
which is a less time-consuming algorithm. The method reduces the econ mic dispatch
problem into an equivalent losslessproblem.
7. M *d Mi t* j = 7,2, ...,L aresolvedusingEqs.(6.111)and(6.172), tively.
8. Calculatethe gradientvectorfrom Eq. (6.173).
If the optimality conditionsare achievedwithin the prescribedaccuracy then GOTO
Step 9, else adjust the expected water discharges using the conjugate g ent method,
and GOTO Step 5.
9 . Calculatethe expectedvalues,variancesand covariancesof storagefrom Eqs. (6.160),
(6.161)and (6.162),respectivelyfor all the hydro generarors.
1 0 . If (m > IvDthen GOTO Srrep11,
else set m = m + I and GOTO Step 4.
11. Check the convergence using the gradientvector given by Eq. (6.177).It nvergence ls
not achieved,then adjust 2a; using the steepestascentmethod.
t 2 . Calculate the overall cost, emission and risk objectives.Calculate the me bership
function lr(Fr)with k = l, 2, ..., K, from Eq. (6.149).
1 3 .Check if the schedulewith new weight combinationsis required.
If 'Yes' then modify the weightsand GOTO Step 2 and repeat
else stop.
where
A t = ( A ; + D i + (c(Pi))2
P i - ( B i + Ei)
From Eq. (6.178)
Lr-F,
Pi (6.180)
2a,
After summing Eq. (6.180) over i, and on simplification,
Lt= N
(6.18r)
s.(-/ 2d,,
r
'
i=l
I
k=L
* oFt"* - Ll ( i = 1, 2 , . . . ,M (6.r82)
N
new= F, + F:"* (6.183)
4
i=1
These are nonlinear equationsin P,s and /,1 and can only be solved iteratively,.Let 21old P-loto
=
(f 1,2, ..., N) be approximated solutionsto Eqs.(6.182)and (6.183).Herethe aim is to nd new
approximations
1,1"*- llta + 6)4
f,new= P--old
+6F,
Power Syste
The expressions ,l,1n"f and flnew (i = 1,2, ..., M) depend on the transmission loss a used.By
tp first order
Thylor'sexpression
Frl"* - o raPTto*.'S$
62rpotd
o L 16
L
'''f . n e w-*
1
40to)
dP, a4
rtt
fr aPiaPi
.L .lFtE
a 4 a 4 ' 3 ry(4n"*-4o"
ry=ugo.$
T=ff.hffi(4"'*-P,"'")
aE',
Putting thesevaluesin Eqs. (6.182) and (6.183),the losslepsproblem can be
N
t.{-/ 'rF,*:=F; (6.18s)
d=l
where
"/V
Pi=F;'f+F,
i=l
Pi= 1"'* Q
a' =dL - 44
To retain the classicf fora and at the same time to improw convergence,it is possi e to include
only the ith term of sumnaption,i.e.
.. [;2P-old
^^r^l d'F,o'u .. II
1'l'o
^L AP:
| +(4n"* -4"t0)l
J
Stochastic Multiobjective Generation Scheduling 479
6.8.3 SampleSystemStudy
A long range hydrothermal problem with four thermal and three hydro electric p ants is
consideredhere. Ttre expectedincremerrtalfuel cost and NO, emissionequationof plants
consideredhere are given below:
Expectedincrementalcost equations($/TvtW-h)
are:
Fu(P() = 0.00767
(,.P1")2
+ 0.80507
Fr^ + 363.7048
Fzz(P{) = 0.00167
(:,Pz^)z
+ 0.80507Fr^+ 363.7048
- r.24885Fr^ + L37.370r
F23(\*) = 0.01378(,P3'")2
- r.24885
Frq(P{) = 0.01378(:,P;')2 P{ + 137.370t
E r , = 0 . 0 0 0 9\ i, i = 0 r . 0 ( i = L , 2 , . . . , 7j;= 1 , 2 , . . . , 7 ; i * j )
Hydro plants
2.8
o
2.6
ttt-ttttI
$?
q)
oo 2.4
tr
c)
(J
L.
o
2.2
2
0.0 0.4
Wbightvu3, whenwz=0.0
Thble 6.40 Expected generation schedule when weights 1ttr w1= 0.8, w2 = 0.0, and w
Thble 6.41 Expected generation schedule when weights ire ld1 = 0.6, w2= 0.2, and
Pr P2 P\ P4 Ps P6 P7
1 52.2826 44.2199 35.1298 35.1298 9.6734 20.3345 .0045
2 50.9489 43.0630 34.2878 34.2878 20.0005 16.1455 1 .r094
3 49.4997 41.8055 33.3729 33.3729 21.5401 12.8879 I .1457
4 44.4189 37.4065 30.1650 30.1650 16.8690 10.1635 I .M46
5 43.1654 36.3222 29.3'735 29.3735 24.4937 16.2737 1 .9634
6 40.2784 33.8265 27.5.503 27,5503 25.8547 27.9998 .9763
7 37.8151 31.6991 25.9,946 25.9946 37.9125 35.7750 .8340
8 34.3518 28.7rr2 23.81J70 23.8070 37.8655 30.7354 .07&
9 24.8012 20.4902 17.7'73r 17.773r 46.6259 28.4318 .4912
10 27.956r 23.2028 19.71665 19.7665 52.8726 23.7285 .5602
11 21.5810 t7.7245 15.7"382 15.7382 66.5194 33.6934 .3091
L2 20.1380 r6.4862 14.8'263 14.8263 65.5869 25.4539 .8206
Sub int. wt = 0.8, ,Nz=0.0, w3= 0.2 wt = 0.6, W2= 0.2, w3= 0.2
Qt Qz Qt Qr Qz
I 6.5585 31 . 1 8 3 6 7.2424 7.1148 3r.4895 7.
2 14.77r0 25.2078 15.6412 14.9532 25.3325 15.739
3 16.2284 20.3914 17.0027 16.3512 20.4918 17. I
4 12.5828 16.rrzt 14.329r 12.9281 16.3529 t4. I
5 1,8.5799 26.2228 :20.1538 18.8530 26.4567 20. 98
'.2t.4299 20.0018 45.985s 2t.
6 19.7593 45.7307
'30.9220 29.4813 59.3567 31.
7 29.296r 59.1508
8 29.0027 50.7031 28.8706 50.7470 29.7 2
45.6405 "29.8922
:i6.4336 33.6552 45.8687 36.46 I
9 33.5294
10 37.5595 38.24M 35.1948 36.7977 37.4751 34.223
11 44.865r 52.2372 38.3460 47.1684 52.6787 39.87 l
t2 50.6312 39.7570 ,+6.9679 47.7347 39.4613 44. 37
0.8
ct)
-Y
t.
-t-
H 0.6 / z - \
\
,.i \\
(A
3tr 0.4
a)
.l'
G4
0.2
0
0.8 0.7 0.6 0.5
W:ight rll1, when wl=0.2
natureof objective,
Figure6.20 Conflicting cost,and expectedNO' emission
expected
solutions. The weight combination presentedat serial number 5 in the Table 6.39 giv the
maximumvalue of p#.The numberof iterationsfor the discussedalgorithmdependson the i tial
guessof dischargeand dual variable)"ai.t\n initial''guessof dischargecan be made by usin the
local variation method.
Power System Optimization
--l--
---lrf
at, , \
tr
o
- a - '
H 0.6 --
,7 x
Ea
Lr
-8 0.4
E
o
r<<
a
0.2
0.4
w3,whenwr-- 0.0
Weight,
ProblemFormulation
6.9.1 StochasticEconomic-Emission
In this section,ttre multiobjectiveswith the equality and inequality constraintspertain g to the
power systemoptimizationproblem are described.The importantnon-commensurable o jectives,
taken into accountare economicoperation,minimal impactson environment,and expect devia-
tions due to unsattisfied
load. The stcchasticmultiobjectiveformulationis adoptedby c sidering
fuel cost coefficients,NO, emission coefficients, and load demand as random varia les. The
stochasticmodels are converte$to their deterministicequivalentsby taking their expec values,
486 Power SYstem OPtimization
;t l-
tic
on
ler
Transmissionline security
bus I to bus 7 must satisfythe foll wing
The active power flow in transmissionline ,n connecting
constraint:
.186)
F#" 3 Pr^sP# ( m = I , 2 , . . . ,N L )
NG
Fr^ .187)
i=l
and
Dt-j,r= AI-j,i * Dt-j,n
.187a)
where
D pi,n= GGDF for line /-7, due to the slack generatorR
NG
Pr,,-\ e,-1,,1
i=l
i*R .187b)
NG
I1
i=l
Xt_t-Xj_, ( .187c)
Xri
and &; is the reactanceof line r-j.
where xp;, xi_iare erementsof trre bus reactancematrix,
to consider three obj ctives:
The stochastic economic emission problem is extended
generationmi hto
(i) expectedfuel cost,(ii) expectedNo, emission,and (iii) varianceof
meerrhe expectedloaddemandwithin (a) ttreexpectedgeneration limits and (b) the tra slon
as
line flow bounds. The multiobjective optimization problem is defined
Stochastic Multiobiective Generation Schedulin
(6.188a)
Minimize
(6.188b)
subject to 1
( t = 1 , 2 ," ' , N G ) (6.188c)
4 m i n< 1 < f , n a x
NG
Fz= I + I,J
(1+ R,,nc,,cflai\
ttt + C, +ZRo,oco,c4)V,1'+
(6.188f)
i=l
NG NG NG
j*i
NG NG NG
(6.188h)
FL= 0+c?)8,,1'+) I tt +Rp,Pic4cr)FiEtiFi
i=1 j=l j=l
j*i
where
d,, b;, and ci are the expected cost coefficients
6.9.2 MembershiPFunctions
by fuzzy sel using the
In an approximatereasoning,logical decision making can be defined
the membershi' functions.
operating condirtions.The fuzzy iets are defined by equationscalled
sets using the ues from
These functions representthe degree of membership in some fazzy
while I mea full com-
0 to l. The mernbershipvalue 0 indicatesincompatibility with the sets,
patibility. when neither is true, a value between 0 and I is taken.
488 Power SystemOptimization
Objectivefunctions tt(Fx)
By takingaccount,of the calculatedindividual minimum and maximum of each objectiv turr-
tion togetherwith the rate of increaseof membershipof satisfaction(Figure 6-22), the d slon
maker must determine his membership function lt(F), in a subjective manner. Here, it is a umed
ttratLt(F) is a strictlymonotonic decreasing and conttnuous function with respect to 4 and is
givenbelow:
I ; Fr 3 4min
0 : F r 2 4max
Fk
Fmrn
rk F,lu*
P,^ m ( Pt^n
'm
p-min> P- >F:
rm rm rm
[t(Pr^) = 6.1e0)
F;^rFr^rFt*
-P-to
P^
^ >
- ' l ^
I M
where Pf^ is the mean of the lower and upper limits of line flows.
Stochastic Multiobjective Generation Schedul
E min
T ? " ' T
Ec
t m t m
Generatorlimits IL(P)
The fuzzy sets more accuratelyrepresentthe operationalinequality constraintsof the power
system.Fuzzy membershipmay have a variety of shapesbut for simplicity, here th generator
limits are representedby a rectangularmembershipfunction (Figure 6.24). Mathe tically, the
membershipfunction is defined as
0
(
4 <4'"*
P,mintP-ltPrmin+AP-'
6.9.3 Performancelndex
The objective functions are reformulatedas fuzzy sets and eacn line flow constrain defines a
fuzzy region of acceptability.Keeping in mind that higher the membershipvalue the ter is the
solutiort,the line flow constraintsare viewed as objectivesthat maximize the membershi functions
for each line. The multiobjectiveoptimizationproblem of Eqs. (6.188) can be rewritt n a s
problt by
A multiobjective optimization problem can be changed to a scalar optimization
reference( get)
defining a global function which minimizes the deviationsfrom the so-called
in space of the ob.ective
objective [Osyczka and Davies, 1984]. Any reasonableor desirable point
reference objective. The scalar
is chosen by the decision maker and may be considered as the
is given bel w :
optimizationproblem obtainedfrom the above statedmultiobjectiveproblem
NL+NG+3
NL+NG+3
E = + L w , ( . )u- ) 2 + 6.194)
j=l
where h, it the target membershipfunction for the 7th function and R is a constant'
the output layer, send information to the user. The hidden neuronsform the hidden I yers in
between, and store the information obtainedthrough training.
A schematicof a three-layeredfeedforwardneural network model is shown in Fig 6.25.
There is a connectionstrength,synapsesor weight associatedwith each connection.Each neuron
can have multiple inputs while there can be only one output. The inputs to a neuron uld be
from externalstimuli or could be from the output of the other neurons.In the simplest each
neuron producesits output by computingthe inner product of its input and
signe.l(s) ociated
weightswhich is passedthrougha nonlinearfunction as shown in Figure 6.26. One monly
used nonlinearmonotonicfunction is the sigmoidal one, which can be defined as follo
1.0 (6.1es)
flx) - 1.0+ exp (-x)
Outputs
Outputlayer
frt, wLz
Hidden
layer
Input layer
Inputs I
neural network.
Figure 0.25 Schematicof feedforwardthree-layered
n"ti =,!r*iiPi
model of neuron.
Figure 6.26 Mathematical
F-
4gZ Power SYstem OPtimization
I
by modifying the interconnection strengthsamong n( rons i
information.This is accomplished
permits htgh computa ional
according ro some prescribed rules. Further, the network architecture
processing. It is also due t this
rare ro be obtained through the rnassively parallel distributed
kind of processingthat such networkshavea greatdegreeof robustness or fault tolerancetr local
famousbei g the
damages.Variouslearningalgorithmshave been recentlyintroduced,the most
backpropagationalgorithm.
Algorithm
6.9.5 BackproPagation
Backpropagation learningalgorithrnfinds the valuesof all of the weightsthat minimi the
prese'ntt, the
function using a methodof gradientdescent.That is, after each patternhas been
error gradient tc ards
error on that pattern is computed and each werght rs moved down the
value for rhat pattern. The error is actually defined by Eq. (6.194), wh.ere h'
is the
its minimum
component f the
target objective for the 7th component of the output patteru and /rl.) is the 7th
pattern. The nett rk is
actual objective produced by the network representationwith the input
specifiedas
Fj=fi(net1)- .1e6)
_
I+e
+
ne!
\rN wirPr .197)
^L
k=l
/ to be adusted.
where 141is the weight on connectionfrom unit 7 at layer (/-1) to unit k atlayer
is used is
To obtain a rule for adjustingweights,the gradientof E with respectto )t r
representedas follows:
AE = .1e8)
6iFo
dwp
6.9.6 SamPleSYstemStudY
The appticabilityof the methodis demonstratedon a samplethree-generatorpower s whose
expected cost and emission characteristicsare given in Tables 6.43 and 6.44, pectively.
ExpectedB-coefficientsfor transmissionloss are depictedin Table 6.45. Table 6.' showsthe
exiected GGDFs (Secrion3.10.5).In addition,the valuesof the CVs and CCs are tak as 0.1 and
1.0, respectively,for all ranciomvariables.
Plant Ai bi ci P,t* 4*
no. ($/Iuw2h) ($fvtwtrl ($n; (Mw) (MW)
PIarrt di €i fi
(kslTr4w2h) (kg/Iawh) (kg/h)
no.
-0.003506 -0.036788
0.02725r
-0.003506 0.030896 -0.005653
-0.036788 -0.005653 0.32295
the range of objective levels may be decided by the experienceddecision maker the
expectedsolution trend.
Neuralnetworkdesign
A three-layeredANN is formed with three neurons in input, hidden, and output layers net-
work is trained with backpropagationalgorithm to achieve the target values given in e 6.49,
whereasthe input given is ob',ainedfrom the minimum cost dispatch and is depicted in le 6.48.
The valuesof n and a arechosenas 0.000005and 0.5 respectively. The valueof R is its
5m. The strengthof connections as weightsis gtven in Table6.50.
The choiceof the number of neuronsis a difficult task.But goodresultsareachie if the
numberof neuronsis equalto the numberof patterns.The numberof iterationsdependon initial
weights and the value of n. The choice of membershipfunction is again a crucial point be the
objective functions are in conflict. Choice becomeseasy if the solution trend is known pnon.
A theofetical basis and methodology for optimal dispatch problem in a uni multi-
objective frar.neworkis established.The chapter also investigatesthe feasibility of q itative
representationof inaccuraciesand uncertaintiesof the input data and power demandfor power
dispatch problem in terms of probability and statistics.An artificial neural network lis
establishedto capturethe optimal generationdispatchfor power system operationsrvith ultiple
conflicting objectives
Stochast,ic Multiobiective Generation Scheduling 495
J wjt W;'t
Layer 1
I 72.r7046 -32.42737 - t 7 .75860
2 12.1,5827 -90.55445 .64217
3 14.84998 -r07.18450 .38823
Layer 2
I 2.62108 -40.59315 .89007
2 1.22125 42.6s099 .73229
3 0.73992 -59.48106 .13109
-
determining the optimal alternativecan be executedon personal computers. This featu also
facilitatesthe applicationsof the new technologyto me decisions.
REFERENCES
Books
Anillaga, J. and C.P.Arnold, ComputerAnalysisof PowerSystems,JohnWiley & Sons,Sin
I 990.
Christensen,G.S. and S.A. Soliman,Optimal Long-TermOperationof Electric Power S
PlenumPress,New York, 1988.
Elgerd,O.I., Electric EnergySystemsTheory:An Introduction,2ndedn.,Tata McGraw-Hill 1 9 8 3 .
El-Hawary,M.E. and G.S. Christensen,OptimalEconomicOperationof Power Systems,Ac tc
Press,New York, 1979.
Frederick Soloman, Probability and StochasticProcesses,PrenticeHall, New Jersey,1987-
Haimes,y.y., Hierarchical Analysisof Water ResourceSystem:Modeling and Optimi on of
I-arge Scale Systems,McGraw-Hill, New York, 1977'
Kirchmayer,L.K., EconomicOperationof PowerSystems,Wiley EasternLtd., New Delhi, 9 5 8
Klir, G.J. and B.yuan, Fuzzy Setsand Fuzzy Logic: Theory and Applications, Prentice- l o f
India, New Delhi, 1997.
Klir, G.J. and T.A. Folger, Fuu! Sets,(Jncertaintyand Information, Prentice-Hallof India, t993.
Kosko, B., Neural Networks and Fuzzy Systems:Dynamical Systems Approach to chine
Intelligence, Prentice-Hallof India, New Delhi, 1994.
Kusic, G.L., ComputerAided Power SystemsAnalysis,Prentice-Hallof India, New Delhi, 986.
i.s and
Mahalanabis,A.K., D.P. Kothari and S.I. Ahson, ComputerAided Power SystemAnaly
Control,TataMcGraw-Hill, New Delhi, 1991.
Modern PowerSystemAnalysis,TataMcGraw-Hill, New Del 1989.
Nagrath,I.J. and D.P.Koth a"J^,
Nagrath,I.J. and D"P.Kothari, PowerSystemEngineering,TataMcGraw-Hill, New Delhi, 994.
Osyczka,A. and B.J. Davres,Multicriterion Optimizationin Engineeringwith FORTRANP
Ellis Horwood Ltd., 1984.
papoulis,A., probability, Random Variablesand StochasticProce.tses,Tata McGraw-Hil , New
Delhi, L99l
Rao, S.S., Optimization:TheoryM Applications,Znd edn., Wiley EasternLimited, New Delhi,
1987.
Sen Gupta,J.K., StochasticProgramming,North Holland, 1972.
Analysis,McGraw-Hi l, New
Stagg,G.W. and A.H. Ei-Abiad, ComputerMethodsin Power'systems
Delhi, 1968.
Stochastic Multiobiective Generation Sc
Papers
Blaszczynski, G.M., Sensitivity study of the economic dispatch, Proceedings of PICA rence,
New Orleans,LA, LT.S.A.,L975.
Brar, Y.S., J.S. Dhillon, and D.P. Kothari, Genetic-fuzzylogic basedweightagepatternfi multi-
objective load dispatch, Asian Journal of Information Technology,Vol. 2, No. 4, pp.
-373,
Oct.-Dec.,2003.
Brar, Y.S., J.S. Dhitlon, and D.P. Kothari, Multiobjective load dispatch by fuzzy log c based
weightagepattern, Electric Power SystemResearch.Vol. 63, pp. 149-L60, 2002.
Chowdhury, N. and R. Billinton, Risk constrainedeconomic load dispatch in intet nnected
generatingsystems, IEEE Trans.on Power Systems,Vol. 5(4), pp. L239-L247, 1990.
Cohon, J.L. and D.H. Marks, A review and evaluationof multiobjectiveprogramming hniques,
WaterResourcesResearch,Vol. ll(2), pp. 208-220, 1975.
David, A.K. andZ. Rongda,An expert systernwith fuzzy setsfor optimal planning, IEEE rans. 04
Power Sysrens,Vol. 6(1), pp.59-65, 1991.
Dhillon, J.S. and D.P. Kothari, The surrogateworth tradeoff approachfor multiobj thermal
power dispatchproblem, Electric Power SystemResearch,Vol. 56, No. 02, pp. 103-1 0, Sept.,
2000.
Dhillon, J.S., S.C. parti, and D.P. Kothari, Multiobjective optimal thermal power dispatc , Int. J.
Electric Power & EnergySystens,Vol. 16, No. 6, pp. 383-390, Dec., L994'
Dhillon, J.S., S.C. parti, and D.P. Kothari, Stochasticeconomicemissionload dispatch,Electric
Power SystemResearch,Vol. 26, pp. 179-186, 1993.
-term
Dhillon, J.S., S.C.parti, and D.P. Kothari, Fuzzy decision making in multiobjective I
schedulingof hydrothermalsystem,Int. J. Electric Power & Energy Systems, Vsl. 23 No. 01,
pp. L9-29, Jan., 2001.
Dhillon, J.S., S.C. Parti, anciD.P. Kothari, Multiobjective Decision Making in Stochastic
Dispatch, Electric Machines an^dPower ,Sys/ems, Vol. 23, pp. 289-30I, 1995.
Dillon, T.S., R.W. Martin, and D. Sjelvgren,Stochasticoptimizationand modellingof large
hydrothermalsystemsfor long-termregulation,Int. I. Electrical Power & Energy Systems,
Vol. 2(l), pp. 2-20, 1980.
Dillon, T.S., S. Sestito,and S. Leung, Short terrn load forecastingusing an adaptiveneuralnetwolh
Int. Journal of Electrical Power & Energy $ysrems,Vol. 13(4), pp. I8GI92. 1991.
Dopazo,J.F.,O.A. Klitin, and A.M. Sasson,Stochasticload flows, IEEE Trans.on Power
Vol. 94(2),pp. 299409, L975.
and Systerns,
495 Power SYstem OPtimization
Edwin, K.W. and R.D. Machate, Influence of inaccurateinput data on the optimal sh -term
operation of power generationsystems,IFAC Symposiumon Automatic Control in Power
Generation,Distribution and Protection,Fretona, Souilr Africa, 1980.
El-Hawary,M.E. and G.A.N. Mbamalu,A comparisonof probabilisticperturbationand de nistic
basedoptimal power flow solutions,IEEE Trans.on PowerSystem.s, Vol. 6(3),pp.l l 105,
r99r.
El-Hawary,M.E. and G.A.N. Mbamalu,Stochasticoptimal load flow using Newton-
iterative technique,Electric Machinesand Power,Sysrerus,
Vol. 15, pp. 371-380, 1988
El-Hawary,M.E. and G.A.N. Mbamalu,Stochasticoptimal load flow using a combined quasi-
Newton and conjugategradienttechnique,Int. Journal of Electrical Power & Energy stems,,
Vol. ll(2), pp. 85-93, 1989.
El-Hawary,M.E. and K.M. Ravindranath, A generaloverviewof multiple objectiveoptima power
flow in hydrothermalelectric power systems,Electric Machinesand Power Systems, t9,
pp. 313-327, 1991.
Glimn, A.F., L.K. Kirchmayer,G.W. Stagg,and V.R. Peterson,
Accuracyconsiderations
in nomic
dispatchof powensystems,A.I.E.E.Trans.on PowerApparatusand Systems, Vol. 75, ( -ilD,
pp. ll25-ll3l, 1956.
Haimes, Y.Y. and W.A. Hall, Multiobjectivesin water resourcesystemsanalysis:The S gate
worth trade-off method, WaterResourcesResearch,Yol. l0(4), pp. 61 5424, 1974.
Hannan, E.L., Linear programmingwith multiple fuzzy goals,Fuuy setsSystem,Yol. 6, p 235-
249,1991.
Hill, E.F. and W.D. Stevenson,A new method of determiningloss coefficients,IEEE
Power Apparatusand systems,Yol. 87(7), pp. 1548-1553,1969.
Hsu, Y.Y. and C.C. Yang, Designof artificial neural networksfor short-termload forecastin . Part
II: Multilayer feedforwardnetworksfor peak load and valley load forecasting, IEE Proc eding,
Part C, Vol. 138(5),pp. 41H18, 1991.
Kaunas,J.R. and Y.Y. Haimes, Risk managementof groundwatercontaminationin a multi-
objectiveframework,WaterResourcesResearch,Yol. 21(11), pp. 172l-1730, 1985.
Kothari, D.P. and I.J. Nagrath, Optimal stochasticscheduling of hydrothermal systems using
discretemaximumprinciple,JournAIInstitutionof Engineers(lndia),Yol.6L,pp. 22-26 1980.
Leberling,H., On finding compromisesolution in multicriteriaproblemsusing the fuzz
operator,Fuzzy,SersSystem,YoI.6, pp. 105-118,1981.
Lee, K.Y., Y.T. Cha, and J.H. Park, Short-termload forecastingusing an artificial neural n ork,
IEEE Trans.on PawerSystems, Vol. 7(1), pp. lZ+132, lg9L.
Leite da Silva,A.M., R.N. Allan, S.M. Soares,and V.L. Arienti, Probabilisticload flow c
network outages,IEE Proceedings,Part C, Vol. 132(3),pp. 139-14s, 19g5.
Mazumdar,M. and C.K. Yin, Varianceof power generatingsystemproductioncosts,IEEE
on Power System.s, Vol. 4(2), pp. 662467, 1989.
Meliopoulos,A.P., X. Chao, G.J. Cokkinides,and R. Monsalvatge,Transmissionloss eval
ation
basedon probabilistic power flow, IEEE Trans. on power System^r,Vol. 6(1), pp. 3 7t ,
1991.
h
Stochastic Multiobjective Generation Schedu
Miranda, V. and J.T. Saraiva,Fuzzy nrodellingof power systemoptimal load flow, IEEE Trans. on
Power Syslerns,Vol. 7(2), pp. 843-849, 1992.
Stochasticgenerationexpansionplanning by meansof
Mo, B., J. Hegge,and I. Wange.nsteen,
stochasticdynamic programming,IEEE 'Trans.on Power Systems,Vol. 6(2), pp. , lggl.
Nanda, J., D.P. Kothari, and K.S. Lingamurthy,A new approachto economic and lnlmum
emissiondispatch,Journal Indian Institute of Science,Yol. 67, pp. 249-256, 1987.
stochasticdynamic pro
Neto, T.A.A., M.V.F. Pereira,and J. Kelman,A risk-constrained g
approachto the operationplanning of hydrothermalsystems,IEEE Trans, on Power pparatus
and Systerus,Vol. 104(2), pp. 273-279, 1985.
Ouyang, Z. and S.M. Shahidehpour,A hybrid artificial neural network-dynamicpro mmlng
approachto unit commitment,IEEE Trans.on Power System^s, Vol. 7(l), pp.236-24 1992.
Park, D,C., M.A. El-Sharkawi, and R.J. Marks II, An adaptively trained neural network, I Trans.
on Neural Networks,Vol. 2(3), pp. 334-345, 1991.
Parti, S.C., D.P. Kothari, and P.V.Gupta,Economicthermalpower dispatch,Joumal of nstitution
of Engineers (India), Vol. 63(EL-Z), pp. 126-132, 1983.
Parti, S.C., Stochasticoptimal power generationscheduling,Ph.D. (Thesis),TIET, Pati 1987.
Pereira,M.V.F., Optimal stochasticoperationsschedulingof large hydroelectricsy ms, Int.
Journal of Electrical Power & Energy Systems,Yol. 11(3), pp.16l-169, 1989.
Rarig, H.M. and Y.Y. Haimes, RisklDispersionindex method, IEEE Trans. on Systems,Man and
Cybernetics,Vol. L3(3), pp. 317-328, 1983.
Rashid,A.H.A. and K.M. Nor, An efficientmethodfor optimal schedulingof fixed-head ydro and
thermal plants, IEEE Trans.on Power Systems,Vol. 6(2), pp.632436, 1991.
Rau, N.S. and C. Necsulescu,Probabilitydistributionsof incrementalcost of prod tion and
production cost, IEEE Trans. on Power Apparatus and Sysfens, Vol. 104(12), pp. 3 93-3499,
1985.
Sasso4,A.M., Combined use of the Powell and Fletcher-Powellnonlinear programmin methods
for optimal load flows, IEEE Trans.on PowerApparatusand Systems,Vol. 88, pp. 1 30-1535,
1969.
Sasson,A.M., Non-linear programming solutions for load-flow, minimum loss and onomic
dispatchingproblems,IEEE Trans.on Power Apparatus and Systems,Yol.88(4), pp. 399409,
1969,
Sherkat, V.R.,R. Campo,K. Moslehi,andE.O.Lo, Stochastic long-termhydrothermal mization
for a multireservoirsystem, IEEE Trans. on Powq Apparatus and Systems, 104(8),
pp. 204V2050,1985.
Su, C.C.and Y.Y. Hsu, Fuzzydynamicprogramming:an applicationto unit commi nt, IEEE
Trans.on PowerSystem.s, Vol. 6(3),pp. L23t-1237,1991.
Tapia, C.G. and B.A. Murtagh, Interactive fuzzy programming with preferencecriteri in multi-
objective decision making, ComputersOperationsResearch,Vol. 18(3), pp. 307-31 , 1991.
Tsuji, A., Optimal fuel mix and load dispatchingunder environmentalconstraints,IEEE Trans.on
Power Apparatus and Sysrens,Vol. 100(5), pp. 2357-2364, 1981.
500 PowerSy opt
with severalobjectivefut
linearprogramming
8.