EE2351 Uw

EE2351
POWER SYSTEM ANALYSIS
A.Ahamed Riazudeen
EEE DEPARTMENT
EE2351 A.Ahamed Riazudeen 1

AP/EEE
UNIT I
INTRODUCTION

AP/EEE
Power system network

AP/EEE
SINGLE LINE DIAGRAM
It is a diagrammatic representation of a power system in

which the components are represented by their symbols.

AP/EEE
COMPONENTS OF A POWER SYSTEM
1.Alternator
2.Power transformer
3.Transmission lines
4.Substation transformer
5.Distribution transformer
6.Loads

AP/EEE
MODELLING OF GENERATOR AND
SYNCHRONOUS MOTOR
1Φ equivalent circuit of generator 1Φ equivalent circuit of synchronous motor

AP/EEE
MODELLING OF TRANSFORMER
E2 N I
K   2  1
E1 N1 I2
R2
R 01  R1  R 2 '  R1  =Equivalent resistance referred to 1o
K2
X
X 01  X 1  X 2 '  X 1  22 =Equivalent reactance referred to 1o
K
AP/EEE
MODELLING OF TRANSMISSION LINE
Π type T type

AP/EEE
MODELLING OF INDUCTION MOTOR
1
'
R r (  1 ) =Resistance representing load
s
R  R S  R r ' =Equivalent resistance referred to stator
X  X S  X r
'
=Equivalent reactance referred to stator

AP/EEE
per unit=actual value/base value
Let KVAb=Base KVA
kVb=Base voltage
Zb=Base impedance in Ω
 kVb   kVb 
2 2
Zb  
MVAb KVAb
1000

AP/EEE
Changing the base of per unit quantities
Let z = actual impedance(Ω)

Z b = base impedance (Ω)
Z Z Z * MVAb
Z p.u   
Zb  kVb  2
 kVb 
2
MVAb
Let kVb ,old & MVBb ,old represent old base values
kVb ,new & MVBb ,new represent new base values

AP/EEE
Z * MVAb , old
Z p .u , old   (1)
 kV 
2
b , old
Z p .u , old * MVAb , old

Z   (2)
 kV 
2
b , old
Z * MVAb , new
Z p .u , new   (3)
 kV 
2
b , new
 
2
kV MVAb , new
Z p .u , new  Z p .u , old
b , old
* *
 kV 
2
b , new
MVAb , old

AP/EEE
ADVANTAGES OF PER UNIT CALCULATIONS
 The p.u impedance referred to either side of a 1Φ

transformer is same
 The manufacturers provide the impedance value in p.u
 The p.u impedance referred to either side of a 3Φ
transformer is same regardless of the 3Φ connections Y-
Y,Δ-Y
 p.u value always less than unity.

AP/EEE
IMPEDANCE DIAGRAM
• This diagram obtained by replacing each component by

their 1Φ equivalent circuit.
Following approximations are made to draw impedance
diagram
1. The impedance b/w neutral and ground omitted.
2. Shunt branches of the transformer equivalent circuit
neglected.

AP/EEE
REACTANCE DIAGRAM
 It is the equivalent circuit of the power system in which

the various components are represented by their
respective equivalent circuit.
 Reactance diagram can be obtained after omitting all

resistances & capacitances of the transmission line from
impedance diagram

AP/EEE
REACTANCE DIAGRAM FOR THE GIVEN POWER SYSTEM
NETWORK

AP/EEE
PROCEDURE TO FORM REACTANCE
DIAGRAM FROM SINGLE DIAGRAM
1.Select a base power kVAb or MVAb
2.Select a base voltage kVb
3. The voltage conversion is achieved by means of transformer kVb on LT
section= kVb on HT section x LT voltage rating/HT voltage rating
4. When specified reactance of a component is in ohms
p.u reactance=actual reactance/base reactance
specified reactance of a component is in p.u
 
2
kV MVAb , new
X p .u , new  X p .u ,old
b , old
* *
 kV 
2
b , new
MVAb , old

AP/EEE
p.u. calculation of 3 winding transformer
Zp=Impedance of primary winding

Zs’=Impedance of secondary winding
Zt’=Impedance of tertiary winding
Short circuit test conducted to find out the above 3
impedances

AP/EEE
1
Zp  Zps  Zpt  Zst ' 
2
1
Zs  Zps  Zst '  Zpt 
'
2
1
Zt  Zps  Zpt  Zst ' 
'
Z ps = Leakage impedance measured in 1o with 2o short circuited

and tertiary open.
Z pt = Leakage impedance measured in 1o with tertiary short
circuited and 2o open.
Z st ' = Leakage impedance measured in 2o with tertiary short

circuited and 1o open and referred to primary

AP/EEE
PRIMITIVE NETWORK
It is a set of unconnected elements which provides information
regarding the characteristics of individual elements. it can be
represented both in impedance & admittance form

AP/EEE
BUS ADMITTANCE(Y BUS) MATRIX
Y BUS can be formed by 2 methods

1.Inspection method
2.Singular transformation
 Y11 Y12   Y1n 

 
Y BUS =
 Y21 Y22   Y2n 
Y   
 n1 n 2
Y Ynn 

AP/EEE
INSPECTION METHOD
For n bus system
Diagonal element of Y BUS
n
Yii   yij
j 1
Off Diagonal element of Y BUS
Yij yij

AP/EEE
SINGULAR TRANSFORMATION METHOD
Y BUS = A  y A
T
Where [y]=primitive admittance
A=bus incidence matrix

AP/EEE
ALGORITHM FOR FORMATION OF THE
BUS IMPEDANCE MATRIX
• Modification of Zbus matrix involves any one of the following 4 cases
Case 1:adding a branch impedance zb from a new bus p to

the reference bus
Addition of new bus will increase the order the Zbus matrix by 1
 zoriginal 0 
Zbus,new   
 0 zb 
(n+1)th column and row elements are zero except the diagonal
diagonal element is zb

AP/EEE
Case 2: adding a branch impedance zb from a new bus p
to the existing bus q
Addition of new bus will increase the order the Zbus matrix by 1
The elements of (n+1)th column and row are the elements of
qth column and row and the diagonal element is Zqq+Zb
Case 3:adding a branch impedance zb from an existing bus p to

the reference bus
The elements of (n+1)th column and row are the elements of
qth column and row and the diagonal element is Zqq+Zb and
(n+1)th row and column should be eliminated using the following
formula
Zj(n1)Z(n1)k
Zjk,act  Zjk  j 1,2...n;k 1,2..n
Z(n1)(n1)

AP/EEE
Case 4:adding a branch impedance zb between existing buses h and q
elements of (n+1)th column are elements of bus h column –
bus q column and elements of (n+1)th row are elements of
bus h row – bus q row the diagonal element= Z b  Z hh  Z qq  2 Z hq
and (n+1)th row and column should be eliminated using the following
formula
Zj(n1)Z(n1)k
Zjk,act  Zjk  j 1,2...n;k 1,2..n
Z(n1)(n1)

AP/EEE
UNIT II
POWER FLOW ANALYSIS

AP/EEE
BUS CLASSIFICATION
1.Slack bus or Reference bus or Swing bus:

|V| and δ are specified. P and Q are un specified, and to be
calculated.
2.Generator bus or PV bus or Voltage controlled bus:

P and |V| are specified. Q and δ are un specified, and to be
calculated
3.Load bus or PQ bus:

P and Q are specified. |V| and δ are un specified, and to be
calculated
AP/EEE
ITERATIVE METHOD
n
I p   YpqVq
q 1
S p  Pp  jQ p  VP I p *
Pp  jQ p n
*
  YpqVq
VP q 1
The above Load flow equations are non linear and

can be solved by following iterative methods.
1.Gauss seidal method

2.Newton Raphson method
3.Fast Decoupled method
AP/EEE
GAUSS SEIDAL METHOD
For load bus calculate |V| and δ from Vpk+1 equation
1  P  jQ p1 n
k
k 1
Vp    YpqVq   YpqVq 
k 1
p p
k *
Ypp  (VP ) q1 q p1 
For generator bus calculate Q from QPK+1 equation

  p1 n
k 

Qp  1*Im (VP ) YpqVq  YpqVq  
k 1 k * k 1
  q1 q p  
AP/EEE
• Check Qp,calk+1 with the limits of Qp
• If Qp,calk+1 lies within the limits bus p remains as PV bus
otherwise it will change to load bus
• Calculate δ for PV bus from Vpk+1 equation
• Acceleration factor α can be used for faster convergence
• Calculate change in bus-p voltage

k1 k1
V p V p V
p
k
• If |ΔVmax |< ε, find slack bus power otherwise increase

the iteration count (k)
• Slack bus power=
S SG L

AP/EEE
NEWTON RAPHSON METHOD
n
Pi  Q i  
j 1
V i V j Y ij  ij   i   j
n
Pi  
j 1
V i V j Y ij c o s ( ij   i   j )
n
Qi  j 1
V i V j Y ij s in ( ij   i   j )
P   J1 J 2   
Q    J  
J 4   V 
   3
 Pi k  Pi sc h  Pi k
 Q i k  Q i sc h  Q i k
AP/EEE
• Calculate |V| and δ from the following equation
i k 1  i k   k
Vi k 1  Vi k   Vi k
Pi k  
• If
Qi k  
• stop the iteration otherwise increase the iteration count (k)

AP/EEE
FAST DECOUPLED METHOD
 J2 & J3 of Jacobian matrix are zero
P   J1   
0
Q    0  
     V 
J4
 P 
 P  J 1     
   
 Q 
Q  J 4 V    V
  V 
P
  B ' 
 Vi
Q
  B ''  V
 Vi
1 P
     B ' 
 V
1 Q
 V    B '' 
 V

AP/EEE
ik1 ik  k
Vik1  Vik  Vik
This method requires more iterations than NR

method but less time per iteration
It is useful for in contingency analysis

AP/EEE
COMPARISION BETWEEN ITERATIVE
METHODS
Gauss – Seidal Method
1. Computer memory requirement is less.
2. Computation time per iteration is less.
3. It requires less number of arithmetic operations to
complete an iteration and ease in programming.
4. No. of iterations are more for convergence and rate of
convergence is slow (linear convergence
characteristic.
5. No. of iterations increases with the increase of no. of
buses.

AP/EEE
NEWTON – RAPHSON METHOD
 Superior convergence because of quadratic convergence.

 It has an 1:8 iteration ratio compared to GS method.
 More accurate.
 Smaller no. of iterations and used for large size systems.
 It is faster and no. of iterations is independent of the no. of buses.
 Technique is difficult and calculations involved in each iteration are
more and thus computation time per iteration is large.
 Computer memory requirement is large, as the elements of jacobian
matrix are to be computed in each iteration.
 Programming logic is more complex.

AP/EEE
FAST DECOUPLED METHOD
 It is simple and computationally efficient.

 Storage of jacobian matrix elements are60% of NR
method
 computation time per iteration is less.
 Convergence is geometric,2 to 5 iterations required for
accurate solutions
 Speed for iterations is 5 times that of NR method and 2-3
times of GS method

AP/EEE
UNIT III
FAULT ANALYSIS-BALANCED FAULT

AP/EEE
Need for fault analysis
 To determine the magnitude of fault current throughout
the power system after fault occurs.
 To select the ratings for fuses, breakers and switchgear.
 To check the MVA ratings of the existing circuit breakers
when new generators are added into a system.

AP/EEE
BALANCED THREE PHASE FAULT
 All the three phases are short circuited to each other and to earth.
 Voltages and currents of the system balanced after the symmetrical fault
occurred. It is enough to consider any one phase for analysis.
SHORT CIRCUIT CAPACITY
 It is the product of magnitudes of the prefault voltage

and the post fault current.
 It is used to determine the dimension of a bus bar and the
interrupting capacity of a circuit breaker.

AP/EEE
Short Circuit Capacity (SCC)
SCC  V 0 I F
VT
IF 
ZT
2
VT S b ,1
SCC 1   MVA / 
ZT ZT p .u
S b ,3
SCC 3
 MVA
ZT p .u
SCC 3
*10 6
If 
3 * V L ,b *10 6

AP/EEE
Procedure for calculating short circuit capacity
and fault current
 Draw a single line diagram and select common base Sb
MVA and kV
 Draw the reactance diagram and calculate the total p.u
impedance from the fault point to source (Thevenin
impedance ZT)
 Determine SCC and If

AP/EEE
ALGORITHM FOR SHORT CIRCUIT ANALYSIS
USING BUS IMPEDANCE MATRIX
• Consider a n bus network. Assume that three phase fault
is applied at bus k through a fault impedance zf
• Prefault voltages at all the buses are
V1 (0) 
V (0) 
 2 
. 
Vbus (0)   
V
 k (0)
. 
 
Vn (0) 
• Draw the Thevenin equivalent circuit i.e Zeroing all voltage sources
and add voltage source V (0) at faulted bus k and draw the reactance
k
diagram
AP/EEE
• The change in bus voltage due to fault is
  V1 
. 
 
. 
 Vbus  

 k V
. 
 
  V n 
• The bus voltages during the fault is

Vbus (F) Vbus (0) Vbus
• The current entering into all the buses is zero.the current entering
into faulted bus k is –ve of the current leaving the bus k

AP/EEE
Vbus  Zbus Ibus
 Z11 . Z1k . Z1n  0 
 
 . . . . .  . 

Vbus   Zk1 . Zkk . Zkn  Ik (F )
  
 . . . . .  . 
Z . Znn  0
 n1 . Znk 
Vk (F )  Vk (0)  Zkk Ik (F )
Vk (F )  Z f Ik (F )
Vk (0)
Ik (F ) 
Zkk  Z f
Vi (F )  Vi (0)  Zik Ik (F )
AP/EEE
UNIT IV
FAULT ANALYSIS – UNBALANCED

FAULTS

AP/EEE
INTRODUCTION
UNSYMMETRICAL FAULTS
o One or two phases are involved
o Voltages and currents become unbalanced and each phase is to be
treated individually
o The various types of faults are
Shunt type faults
1.Line to Ground fault (LG)
2. Line to Line fault (LL)
3. Line to Line to Ground fault (LLG)
Series type faults
Open conductor fault (one or two conductor open fault)

AP/EEE
FUNDAMENTALS OF SYMMETRICAL
COMPONENTS
 Symmetrical components can be used to transform
three phase unbalanced voltages and currents to
balanced voltages and currents
 Three phase unbalanced phasors can be resolved into
following three sequences
1.Positive sequence components
2. Negative sequence components
3. Zero sequence components

AP/EEE
Positive sequence components
Three phasors with equal magnitudes, equally displaced from one
another by 120o and phase sequence is same as that of original
phasors. Va1,Vb1,Vc1
Negative sequence components

Three phasors with equal magnitudes, equally displaced from one
another by 120o and phase sequence is opposite to that of original
phasors. V , V , V
a2 b2 c2
Zero sequence components

Three phasors with equal magnitudes and displaced from one another
by 0o Va0 ,Vb0 ,Vc0

AP/EEE
RELATIONSHIP BETWEEN UNBALANCED VECTORS
AND SYMMETRICAL COMPONENTS
Va  Va 0  Va1  Va 2
Vb  Vb 0  Vb1  Vb 2
Vc  Vc 0  Vc1  Vc 2
Va   1 1 1  Va 0 
V    1 a 2 
a  Va1 
 b 
Vc   1 a a 2  Va 2 
1 1 1
 
A  1 a 2 a
1 a a 2 

Similarly we can obtain for currents also

AP/EEE
SEQUENCE IMPEDANCE
Impedances offered by power system components to positive,

negative and zero sequence currents.
Positive sequence impedance
The impedance of a component when positive sequence
currents alone are flowing.
Negative sequence impedance
The impedance of a component when negative sequence
currents alone are flowing.
Zero sequence impedance
The impedance of a component when zero sequence currents
alone are flowing.
AP/EEE
SEQUENCE NETWORK
SEQUENCE NETWORK FOR GENERATOR
negative sequence network Zero sequence network

positive sequence network

AP/EEE
SEQUENCE NETWORK FOR
TRANSMISSION LINE
positive sequence network negative sequence network Zero sequence network

AP/EEE
SEQUENCE NETWORK FOR
TRANSFORMER

AP/EEE
SEQUENCE NETWORK FOR LOAD

AP/EEE
SINGLE LINE TO GROUND FAULT
Ib  0
Ic  0
Va  Zf Ia
Ia1  Ia2  Ia0  Ia /3
Ea
Ia1 
Consider a fault between phase a and Z1  Z2  Z0 3Z f
ground through an impedance zf

AP/EEE
LINE TO LINE (LL) FAULT
Ia  0
I c   Ib
Vb  Vc  Ib Z f
Ia2  Ia1
Ia0  0
Consider a fault between phase b and c Va1  Va2  Zf Ia1
through an impedance zf Ea
Ia1 
Z1  Z2  3Z f
 jEa
Ib   I c 
Z1  Z2  3Z f

AP/EEE
DOUBLE LINE TO GROUND (LLG) FAULT
Ia0  0
Ia1  Ia2  Ia0  0
Vb =Vc  Z f (Ib  Ic )  3Zf Ia0
Va0 Va1  Vb  3Zf Ia0
Ea
Ia1 
Z1  Z2 (Z0 3Z f )/(Z2  Z0 3Z f )
Consider a fault between phase b and c
through an impedance zf to ground

AP/EEE
UNBALANCED FAULT ANALYSIS USING
BUS IMPEDANCE MATRIX
SINGLE LINE TO GROUND FAULT USING Zbus
 Consider a fault between phase a and ground through
an impedance zf at bus k
For a fault at bus k the symmetrical
components of fault current
Vk (0)
Ik 0  Ik 1  Ik 2 
Zkk1  Zkk 2  Zkk 0  3Z f

AP/EEE
LINE TO LINE (LL) FAULT
Consider a fault between phase b and c through an impedance zf
Ik 0  0
Vk (0)
Ik  Ik
1 2

Z kk 1  Z kk 2  Z f

AP/EEE
DOUBLE LINE TO GROUND (LLG) FAULT
Consider a fault between phase b and c through an impedance zf to ground
Vk (0)
Ik1 
Z 2
( Z 0
 3 Z f
)
Z kk 1  kk2 kk 0
Z kk  Z kk  3Z f
Vk (0)  Z kk 1I k 1
Ik 2 
Z kk 2
Vk (0)  Z kk 1I k 1
Ik 0

Z kk 0  3Z f
Ik ( F )  Ik b  Ik c

AP/EEE
BUS VOLTAGES AND LINE CURRENTS
DURING FAULT
Vi 0 ( F )  0  Z ik 0 I k 0
Vi1 ( F )  Vi 0 (0)  Z ik 1 I k 1
Vi 2 ( F )  0  Z ik 2 I k 2
Vi 0 ( F )  V j 0 ( F )
I ij 0 
Z ij 0
Vi 1 ( F )  V j 1 ( F )
I ij1 
Z ij1
Vi 2 ( F )  V j 2 ( F )
I ij 2 
Z ij 2

AP/EEE
UNIT V
STABILITY ANALYSIS

AP/EEE
STABILITY
 The tendency of a power system to develop restoring
forces equal to or greater than the disturbing forces to
maintain the state of equilibrium.
 Ability to keep the machines in synchronism with another
machine

AP/EEE
CLASSIFICATION OF STABILITY
 Steady state stability

Ability of the power system to regain synchronism after small and
slow disturbances (like gradual power changes)
 Dynamic stability
Ability of the power system to regain synchronism after small
disturbances occurring for a long time (like changes in turbine
speed, change in load)
 Transient stability
This concern with sudden and large changes in the network
conditions i.e. . sudden changes in application or removal of loads,
line switching operating operations, line faults, or loss of excitation.

AP/EEE
 Steady state limit is the maximum power that can be
transferred without the system become unstable when
the load in increased gradually under steady state
conditions.
 Transient limit is the maximum power that can be

transferred without the system becoming unstable when
a sudden or large disturbance occurs.

AP/EEE
Swing Equation for Single Machine Infinite
Bus System
• The equation governing the motion of the rotor of a
synchronous machine
d 2m
J 2  Ta  Tm Te
dt
where
J=The total moment of inertia of the rotor(kg-m2)
 m =Singular displacement of the rotor
Tm=Mechanical torque (N-m)
Te=Net electrical torque (N-m)
Ta=Net accelerating torque (N-m)

AP/EEE
m  smt   m
dm d m
 sm 
dt dt
d 2m d 2 m
2
 2
dt dt
d 2 m
J m 2  pa  pm  pe
dt
• Where pm is the shaft power input to the machine
pe is the electrical power
pa is the accelerating power

AP/EEE
Jm  M
d 2 m
M  pa  pm  pe
dt 2
2H
M  S machine
 sm
2 H d 2 m pa pm  pe
 
 sm dt 2 S machine S machine
2 H d 2 H=machine inertia constant
 pa  pm  pe
 s dt 2
 s  2 f
H d 2
 pa   pm   pe
 f 0 dt 2
d 2  f0  f
2
  p m  p 2 max sin    0 p a p.u δ and ωs are in electrical
dt H H
d radian
 
dt
d   f0 d 2
 pa  p.u
dt H dt 2

AP/EEE
Swing Equation for Multimachine System
S m a c h in e =machine rating(base)
S system =system base
Hsystem d 2
 pa  pm  pe p.u
f dt 2
Smachine
Hsystem  Hmachine
Ssystem

AP/EEE
Rotor Angle Stability
• It is the ability of interconnected synchronous machines

of a power system to maintain in synchronism. The
stability problem involves the study of the electro
mechanical oscillations inherent in power system.
• Types of Rotor Angle Stability
1. Small Signal Stability (or) Steady State Stability
2. Transient stability

AP/EEE
Voltage Stability
 It is the ability of a power system to maintain steady

acceptable voltages at all buses in the system under
normal operating conditions and after being subjected to
a disturbance.
 The major factor for instability is the inability of the
power system to meet the demand for reactive power.

AP/EEE
• Mid Term Stability
It represents transition between short term and long
term responses.
Typical ranges of time periods.
1. Short term : 0 to 10s
2. Mid Term : 10 to few minutes
3. Long Term : a few minutes to 10’s of minutes
• Long Term Stability
Usually these problem be associated with
1. Inadequacies in equipment responses.
2. Poor co-ordination of control and protection equipment.
3. Insufficient active/reactive power reserves.

AP/EEE
Equal Area Criterion
• This is a simple graphical method to predict the transient

stability of two machine system or a single machine
against infinite bus. This criterion does not require swing
equation or solution of swing equation to determine the
stability condition.
• The stability conditions are determined by equating the
areas of segments on power angle diagram.

AP/EEE
Power-angle curve for equal area criterion
multiplying swing equation by on both sides
Multiplying both sides of the above equation by dt and then integrating between two arbitrary angles δ0 and δc

AP/EEE
Once a fault occurs, the machine starts accelerating. Once the fault is cleared, the
machine keeps on accelerating before it reaches its peak at δc ,
The area of accelerating A1
The area of deceleration is given by A2
If the two areas are equal, i.e., A1 = A2, then the power system will be stable

AP/EEE
Critical Clearing Angle (δcr) maximum allowable value of the clearing time and angle for
the system to remain stable are known as critical clearing time and angle.
δcr expression can be obtained by substituting δc = δcr in the equation A1 = A2
Substituting Pe = 0 in swing equation
Integrating the above equation

AP/EEE
Replacing δ by δcr and t by tcr in the above equation, we get
the critical clearing time as

AP/EEE
Factors Affecting Transient Stability
• Strength of the transmission network within the system

and of the tie lines to adjacent systems.
• The characteristics of generating units including inertia of
rotating parts and electrical properties such as transient
reactance and magnetic saturation characteristics of the
stator and rotor.
• Speed with which the faulted lines or equipments can be
disconnected.

AP/EEE
Numerical Integration methods
 Modified Euler’s method
 Runge-Kutta method

AP/EEE
MODIFIED EULER’S METHOD
• Using first derivative of the initial point next point is obtained

X1p  X0  t the step t1  t0 t
• dX
dt
• Using this x1p dx/dt at x1p=f(t1, x1p)
• Corrected value is
  dx   dx  
  dt     p 
   d t  X1
X 1P  X   t
X0
0
 2 
 
 
  dx   dx  

 dt     p 
   d t  X i 1
 Xi   t
c Xi
X i 1
 2 
 
 
AP/EEE
Numerical Solution of the swing equation
• Input power pm=constant

• At steady state pe=pm,
 pm 
 0  s in  1  
 p1 m ax 
E' V
p1 m ax 
X1
• At synchronous speed
0  0
E' V
p 2 max 
X2

AP/EEE
The swing equation
H d 2
 pa   pm   pe
 f 0 dt 2
d 2  f 0  f0
2
  pm  p2max sin    pa
dt H H
d
 
dt
d   f 0 d 2
 pa  2
dt H dt
Applying Modified Eulers method to above equation
t1  t0  t
 d 
 i p1   i    t
 dt  i
 d  
  ip1    i    t
 dt  i

AP/EEE
• The derivatives at the end of interval
 d 
 p    i 1
p

 dt    i 1
 d     f0 
  p   pa 
 dt   i 1  H  p
i 1
The corrected value

  d   d  
  dt     p 
   dt  i1
 i    t
 
 ic1 i
 2 
 
 
  d    d   

 dt     p 
   dt  i1
  i    t

ic1 i
 2 
 
 

AP/EEE
Runge-Kutta Method
• Obtain a load flow solution for pretransient conditions
• Calculate the generator internal voltages behind transient reactance.
• Assume the occurrence of a fault and calculate the reduced admittance matrix
• Initialize time count K=0,J=0
• Determine the eight constants
K 1k  f 1 ( k ,  k )  t
l1k  f 2 ( k ,  k )  t
K 1k lk
K 2k  f 1 ( k
 ,  k  1 )t
2 2
k k
K l
l 2k  f 2 ( k  1 ,  k  1 )  t
2 2
k
K2 l 2k
K 3  f 1 ( 
k k
,   )t
k
2 2
k k
K l
l 3k  f 2 ( k  2 ,  k  2 )  t
2 2
Kk lk
K 4k  f 1 ( k  3 ,  k  3 )  t
2 2
k k
K l
l 4k  f 2 ( k  3 ,  k  3 )  t
2 2
  k  K 1  2 K 2  2 K 3k  K 4k 
k k
  k  l1  2 l 2  2 l 3k  l 4k 
k k
6 EE2351 A.Ahamed Riazudeen 86

AP/EEE
• Compute the change in state vector
 k

 K
1
k
 2 K 2k  2 K3k  K 4k 
6
 k 
l  2l  2l3k  l4k 
1
k k
2
6
• Evaluate the new state vector
 k1  k  k
k1 k k
• Evaluate the internal voltage behind transient reactance using the relation
Epk 1  Epk cos pk 1  j Epk sin pk 1

• Check if t<tc yes K=K+1
• Check if j=0,yes modify the network data and obtain the new reduced
admittance matrix and set j=j+1
• set K=K+1
• Check if K<Kmax, yes start from finding 8 constants

AP/EEE

EE2351 Uw

Uploaded by

Copyright:

Available Formats

EE2351 Uw

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

EE2351 Uw

Uploaded by

Copyright:

Available Formats

EE2351

POWER SYSTEM ANALYSIS

EE2351 A.Ahamed Riazudeen 1

EE2351 A.Ahamed Riazudeen 2

EE2351 A.Ahamed Riazudeen 3

It is a diagrammatic representation of a power system in

EE2351 A.Ahamed Riazudeen 4

EE2351 A.Ahamed Riazudeen 5

1Φ equivalent circuit of generator 1Φ equivalent circuit of synchronous motor

EE2351 A.Ahamed Riazudeen 6

EE2351 A.Ahamed Riazudeen 8

EE2351 A.Ahamed Riazudeen 9

EE2351 A.Ahamed Riazudeen 10

Let z = actual impedance(Ω)

kVb ,new & MVBb ,new represent new base values

EE2351 A.Ahamed Riazudeen 11

Z p .u , old * MVAb , old

EE2351 A.Ahamed Riazudeen 12

 The p.u impedance referred to either side of a 1Φ

EE2351 A.Ahamed Riazudeen 13

• This diagram obtained by replacing each component by

EE2351 A.Ahamed Riazudeen 14

 It is the equivalent circuit of the power system in which

 Reactance diagram can be obtained after omitting all

EE2351 A.Ahamed Riazudeen 15

EE2351 A.Ahamed Riazudeen 16

EE2351 A.Ahamed Riazudeen 17

Zp=Impedance of primary winding

EE2351 A.Ahamed Riazudeen 18

Z ps = Leakage impedance measured in 1o with 2o short circuited

Z st ' = Leakage impedance measured in 2o with tertiary short

EE2351 A.Ahamed Riazudeen 19

EE2351 A.Ahamed Riazudeen 20

Y BUS can be formed by 2 methods

 Y11 Y12   Y1n 

EE2351 A.Ahamed Riazudeen 21

Off Diagonal element of Y BUS

EE2351 A.Ahamed Riazudeen 22

Where [y]=primitive admittance

A=bus incidence matrix

EE2351 A.Ahamed Riazudeen 23

Case 1:adding a branch impedance zb from a new bus p to

EE2351 A.Ahamed Riazudeen 24

Case 3:adding a branch impedance zb from an existing bus p to

EE2351 A.Ahamed Riazudeen 25

EE2351 A.Ahamed Riazudeen 26

POWER FLOW ANALYSIS

EE2351 A.Ahamed Riazudeen 27

1.Slack bus or Reference bus or Swing bus:

2.Generator bus or PV bus or Voltage controlled bus:

3.Load bus or PQ bus:

The above Load flow equations are non linear and

1.Gauss seidal method

For load bus calculate |V| and δ from Vpk+1 equation

For generator bus calculate Q from QPK+1 equation

• Calculate change in bus-p voltage

• If |ΔVmax |< ε, find slack bus power otherwise increase

EE2351 A.Ahamed Riazudeen 31