Module 5
Module 5
Module 5
Unit 1 Unit 2
400 MW
1200 MW
SECURE POST CONTINGENCY STATE
Ires
PTDFij = AP; = ai
where I=line index, i= bus where power is
injected. SIMULATE AN OUTAGE
j= bus where power is taken out OF LINE RUSING THE
SYSTEM MODEL
Afj= change in megawatt power flow on line 1
when a change in generation, APi, occurs at bus i
ANY LINE FLOWS YES DISPLAY
AP;= change in generation at bus i EXCEED ;IMIT ALARM MESSAGE
The PTDF factor then represents the sensitivity of ANY BUS VOLTAGES) YES OISPLAY
the flow on line to a shift of power from ito j. ALARM MESSAGE
oUTSIOE LIMITS /
Suppose one wanted to study the outage of a large
generating unit and it was assumed that all the NO
Yij Sk pmax
k#i
4
Basavarajappa S R, E& E, BIET, Davangere
where Pax = maximum MW rating for
Yij proportionality factor for pickup on generator k
generating unit jwhen unit 11a1lS
Ihen, to test for the flow on line 1, under the
participate in making up the assumntion that all the generators in the interconevuou
loss, use the following:
, =f+ PTDFj AP -
PTDF1j4 Yj AP:)
Note that this assumes that no unit will
detailed generation pickup algorithm that took actuallyofhit its maximum, If this is apt to be the case, a more
account generation limits would be required.
2. Line Outage Distribution
LODF factors are Factors: The
used in a similar manner, START
only they apply to the testing for overloads P AT EACH
GEN, BUSs
when transmission circuits are lost. By ON ALL
READ EXISTING
SYSTEM CONDITIONS
definition, the line outage LINES
has the following meaning:distribution factor
Af
LODF,k = d,k
where LODFk= line outage distribution
factor when monitoring line I after an outage AP, =
on line k CHECK ALL LINES
FOR OVERLOAD
Af; = change in MW flow on line DI^PLAY AFTER GENERATOR
NO
ALARM OUTAGES
f =original flowon line k before it MESSAGE
was outaged (opened) YES
used to give rapidanalysis of the system, but Fig. 3 AC power flow security analysis
they cannot give information about MVAR
flows and voltages. Slower, full AC power
flow methods give full accuracy but take too
long.
Because of the way the power system Select te bad casc trom he ul
case tst and store ln a short st
is designed and operated, very few of the
outages will actually cause trouble. That is,
most of the time spent running AC power Shot es o mot
flows willgo for solutions of the power flow Ikaly bsd cascs
model that discover that there are no
problems. Only a few of the power flow
solutions will conclude that an overload or
voltage violation exists. The solution to find a ouages
way to select contingencies in sucha way that Pick outago Ifrom he short ist and
remova that conponent rom he power
only those that are likely to result in an ow modet
the operators being warned. Fig. 4 AC power flow security analysis with
contingency case selection
6
Basavarajappa S R, E& E, BIET, Davangere
Contingency Selection and Ranking: An over-load nerlomance index is used to find how much aparticular
ouiage might atfeet the power system. The definition for the overload performance index (PI) is as follows:
211
Pl = (Pnow.l
pmax as nT, PAow) < Pmax ’ PIJ and Paow) > P’PIT
Allbranches
lf n (a positive integer) is a large number, the P will be a
and it will be large if one or more lines are overloaded, The small thennumber if all flows are within limit,
index. problem is how to use this performance
Various techniques have been tried to obtain the
calculations can be made exactly if n= 1:that is,a tableof PIvalue of PI when a branch is taken out. These
values, one for each line in the network, can be
calculated quite quickly. The selection procedure then involves ordering the PI
least. The lines corresponding to the top of the list are then table from largest value to
when n =1, the PI does not change suddenly from near zero to picked and placed them on the short list. However,
near infinity
Instead. it rises asa quadratic function. A line that is just below its limit as the branch exceeds its limit.
that is just over its limit. contributes to Pl almost equal to one
The result is a PI that may be large
when many lines are loaded just below their
limit. Thus the Pl's ability to distinguish or Begin power flow solution
detect bad cases is limited when n= 1. Trying
to develop an algorithm that can quickly
calculate PI when n =2 or larger has proven
extremely difficult. Buid B' and B' matrices
One way to perform an outage case
Full outage
selection is to perform what has been called case list
the 1P1Q method. Here, a decoupled power Model outage casse
flow is used. As shown in Fig.5. the solution
procedure is interrupted after one iteration Solve the P-theta equation for
(one P - 0 calculation and one Q - V the A9'S
calculation; thus, the name 1P1Q). With this
procedure, the PI can use as large an n value
as desired, say n = 5. There appears to be eew = A,
sufficient information in the solution at the
end of the first iteration of the decoupled Solve the O-V equation for the
power flow to give a reasonable PI.
Another advantage to this procedure is the
fact that the voltages can also be included in
the PI. Thus,a different PIcan be used, such
as:
2n Calculate flows and voltages for this
PI = /Pnow,!
pmax
case then calculate the Pl
Allbranches
2m
A|E| Pick next outage case
PIL0st
(one entry for
AJJ buses
(4|E|max each outage
i case)
Where AIElis the difference between Fig.5 1P1Q Contingency Selection Procedure
the voltage magnitude as solved at the end of
the 1P1Q procedure and the base-case voltage magnitude. AEma is a value set by utility engineers indicating
how much they wish to limit a bus voltage from changing on one outage case.
To complete the security analysis, the PI list is sorted so that the largest PI appears at the top. The
security analysis can then start by executing full power flows with the case which is at the top of the list, then
solve the case which is second, andso on down the list. This continues until either a fixed number of cases is
solved, or until a predetermined number of cases are solved which do not have any alarms.
Basavarajappa S R, E& E, BIET, Davangere
State Estimation of Power Systems: Introduction:
State estimation is the process of assigning a value to an unknown system state variable based on
measurements from that system according to some criteria. Usually, the process involves imperfect
measurements that are redundant, and the process of estimating the system states is based on a statistical
criterion that estimates the true value of the state variables to minimize or maximize the selected criterion. A
commonlyused and familiar criterion is that of minimizing the sum of the squares of the differences between
the estimated and true" (i.e., measured) values of a function. State estimators maybe both staticand dynamic.
Both types of estimators have been developed for power systems.
In a power system, the state variables are the voltage magnitudes and relative phase angles at the
System nodes. Measurements are required in order to estimate the system performance in real time for both
system security control and constraints on economic dispatch. The inputs to an estimator are imperfect power
system measurements of voltage magnitudes and power, VAR, or ampere-flow quantities. The estimator
isdesigned to produce the "best estimate" of the system voltage and phase angles, recognizing that there are
errors in the measured quantities and that there may be redundant measurements. The output data are then
used in system control centers in the implementation of the security-constrained dispatch and control of the
system.
When normally distributed, unbiased meter error distributions are assumed, each of these approaches
results in identical estimators. Thus, utilize the maximum likelihood approach because the method introduces
the measurement error weighting matrix R]in a straightforward manner.
The maximum likelihood procedure asks the following question:What is the probability (or likelihood)
that I will get the measurements I have obtained? This probability depends on the random error in the
measuring device (transducer) as well as the unknown parameters to be estimated. Therefore, a reasonable
procedure would be one that simply chose the estimate as the value that maximizes this probability. The
maximum likelihood estimator assumes that know the probability density function (PDF) of the random errors
in the measurement.
The least-squares estimator does not require that know the PDF for the sample or measurement errors.
Assume that the PDF of sample or measurement error is anormal (Gaussian) distribution, it will end up with
the same estimation formula. Hence, proceed to develop estimation formula using the maximum likelihood
criterion assuming normal distributions for measurement errors. The result will be a least-squares or more
precisely a weighted least-squares estimation formula.
First, we introduce the concept of random measurement error. The measurements are assumed to be in
error: that is, the value obtained from the measurement device is close to the true value of the parameter being
measured but differs by an unknown eror. Mathematically, this can be modelled as follows.
Let zmeas be the value of a measurement as received from a measurement device.
Let z'rue be the true value of the quantity being measured. Finally, let n bethe random measurement
error. Then represent measured value as
zmeas = ztrue + n ’ (5.1)
Basavarajappa SR, E&E, BIET, Davangere 8
ne random number,n, serves to model the uncertainty in the measurements. If the measurement error
1S unbiased, the PDF of n is usually chosen as a normal distribution with zero mean. Note that OLnel
measurement PDFs will also work in the maxiMum likelihood method as well. The PDF of nis
1
PDF(n) = oV2 exp ’ (5.2)
where o is called the standard deviation and g² is called the
variance of the random number. PDF(n)
behavior of n. A plot of PDFM) isshown in describes
Fig. 6.
the PDF (n)
Note that o, the standard
to model the seriousness of the deviation, provides a way
random measurement error.
If o is large, the measurement is relatively
inaccurate (i.e.,
apoor-quality measurement device), whereas a small value
of o denotes a small error spread (i.e., a
measurement device). The normal distributionhigher-quality
is commonly 20 30
used for modelling measurement errors since it is the
Fig.6 The normal distribution
distribution that willresult when many factors contribute to
the overall error.