Modelling and Analysis of Electric Power Systems

Modelling and Analysis of Electric Power Systems
Power Flow Analysis Fault Analysis Power Systems Dynamics and Stability
Lecture 227-0526-00, ITET ETH Zrich u
Gran Andersson o
EEH - Power Systems Laboratory
ETH Z rich u March 2006
ii
Contents
Preface vii
Static Analysis
1
1 1 3 3 5 6 9 10 12 14 16 16 17 18 19 20 21 21 23 24 25 27
1 Introduction 1.1 Power Flow Analysis . . . . . . . . . . . . . . . . . . . . . . . 1.2 Fault Current Analysis . . . . . . . . . . . . . . . . . . . . . . 1.3 Literature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 Network Models 2.1 Lines and Cables . . . . . . . . . . . . . 2.2 Transformers . . . . . . . . . . . . . . . 2.2.1 In-Phase Transformers . . . . . . 2.2.2 Phase-Shifting Transformers . . . 2.2.3 Unied Branch Model . . . . . . 2.3 Shunt Elements . . . . . . . . . . . . . . 2.4 Loads . . . . . . . . . . . . . . . . . . . 2.5 Generators . . . . . . . . . . . . . . . . 2.5.1 Stator Current Heating Limit . . 2.5.2 Field Current Heating Limit . . 2.5.3 Stator End Region Heating Limit 3 Active and Reactive Power Flows 3.1 Transmission Lines . . . . . . . . . . . . 3.2 In-phase Transformers . . . . . . . . . . 3.3 Phase-Shifting Transformer with akm = 3.4 Unied Power Flow Equations . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . 1 .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
4 Nodal Formulation of the Network Equations
5 Basic Power Flow Problem 31 5.1 Basic Bus Types . . . . . . . . . . . . . . . . . . . . . . . . . 31 5.2 Equality and Inequality Constraints . . . . . . . . . . . . . . 32 iii
iv 5.3
Contents Problem Solvability . . . . . . . . . . . . . . . . . . . . . . . . 34 37 37 39 40 41 42 44 45 48 48 52 55 57 61 63 66 66 66 67 68 71
6 Solution of the Power Flow Problem 6.1 Solution by Gauss-Seidel Iteration . . . . . . . . . . . . . 6.2 Newton-Raphson Method . . . . . . . . . . . . . . . . . . 6.2.1 Unidimensional case . . . . . . . . . . . . . . . . . 6.2.2 Quadratic Convergence . . . . . . . . . . . . . . . 6.2.3 Multidimensional Case . . . . . . . . . . . . . . . . 6.3 Newton-Raphson applied to the Power Flow Equations . . 6.4 P QU Decoupling . . . . . . . . . . . . . . . . . . . . . 6.5 Approximative Solutions of the Power Flow Problem . . . 6.5.1 Linearization . . . . . . . . . . . . . . . . . . . . . 6.5.2 Matrix Formulation of DC Power Flow Equations 6.5.3 DC Power Flow Model . . . . . . . . . . . . . . . . 7 Fault Analysis 7.1 Transients on a transmission line . . . . . 7.2 Short circuit of a synchronous machine . . 7.3 Algorithms for short circuit studies . . . . 7.3.1 Generator model . . . . . . . . . . 7.3.2 Simplications . . . . . . . . . . . 7.3.3 Solving the linear system equations 7.3.4 The superposition technique . . . . 7.3.5 The Takahashi method . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
II
Power System Dynamics and Stability

. . . . . . .
77
79 80 80 81 82 82 84 87 89 89 90 91 94 95 95
8 Classication and Denitions of Power System Stability 8.1 Dynamics in Power Systems . . . . . . . . . . . . . . . . . . 8.1.1 Classication of Dynamics . . . . . . . . . . . . . . . 8.1.2 Modelling . . . . . . . . . . . . . . . . . . . . . . . . 8.2 Power System Stability . . . . . . . . . . . . . . . . . . . . . 8.2.1 Denition of Stability . . . . . . . . . . . . . . . . . 8.2.2 Classication of Power System Stability . . . . . . . 8.3 Literature on Power System Dynamics and Stability . . . . 9 Synchronous Machine Models 9.1 Design and Operating Principle . . . . . . 9.1.1 Rotor Types . . . . . . . . . . . . 9.1.2 Stator Field . . . . . . . . . . . . . 9.1.3 Magnetic Torque . . . . . . . . . . 9.2 Stationary Operation . . . . . . . . . . . . 9.2.1 Stationary Single Phase Equivalent
. . . . . . . . . . . . . . . . . . . . . . . . . Circuit
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
Contents 9.2.2 Phasor diagram . . . . . . . . . . . . . . . 9.2.3 Operational Limits . . . . . . . . . . . . . Dynamic Operation . . . . . . . . . . . . . . . . 9.3.1 Transient Single Phase Equivalent Circuit 9.3.2 Simplied Mechanical Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
v . 97 . 98 . 100 . 100 . 100
9.3
10 The 10.1 10.2 10.3
Swing Equation 103 Derivation of the Swing Equation . . . . . . . . . . . . . . . . 103 Analysis of the Swing Equation . . . . . . . . . . . . . . . . . 105 Swing Equation as System of First Order Dierential Equations106 109 . 109 . 111 . 113 . 120 . 122 . 123 . 126
11 Power Swings in a Simple System 11.1 The Swing Equation and its Solutions 11.1.1 Qualitative Analysis . . . . . . 11.1.2 Stable and Unstable Solutions . 11.2 Equal Area Criterion . . . . . . . . . . 11.3 Lyapunov Stability Criterion . . . . . 11.4 Small Signal Analysis . . . . . . . . . 11.5 Methods to Improve System Stability
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
12 Oscillations in Multi-Machine Systems 129 12.1 Simplied Model . . . . . . . . . . . . . . . . . . . . . . . . . 129 12.2 Detailed Model . . . . . . . . . . . . . . . . . . . . . . . . . . 130 13 Voltage Stability 13.1 Mechanisms of Voltage Instability . . . 13.1.1 Long Term Voltage Instability 13.1.2 Short Term Voltage Instability 13.2 Simple Systems for Analysis of Voltage 13.3 Analysis of Voltage Stability . . . . . . 13.3.1 Stability Indicators . . . . . . . 13.3.2 Analysis of Simple System . . . 133 . 133 . 134 . 134 . 135 . 138 . 139 . 140 149 151 152 154 155 155 155 156 158 161
. . . . . . . . . . . . . . . . . . Stability . . . . . . . . . . . . . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
14 Control of Electric Power Systems 14.1 Control of Active Power and Frequency . . . . 14.1.1 Spinning reserve . . . . . . . . . . . . . 14.1.2 Supplementary Reserves . . . . . . . . . 14.1.3 Back-Up Reserves . . . . . . . . . . . . 14.2 Control of Reactive Power and Voltage . . . . . 14.2.1 Reactive Power Control . . . . . . . . . 14.2.2 Voltage Control . . . . . . . . . . . . . . 14.3 Supervisory Control of Electric Power Systems A Phase-Shifting Transformers
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
vi B Protections in Electric Power Systems B.1 Design of Protections . . . . . . . . . . B.2 Distance Protections . . . . . . . . . . B.2.1 General Principles . . . . . . . B.2.2 Automatic Re-Closure . . . . . B.3 Out of Step Protections . . . . . . . . B.4 System Protections . . . . . . . . . . .
Contents 165 . 165 . 167 . 167 . 169 . 170 . 170
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
Preface
These lectures notes are intended to be used in the Modellierung und Analyse elektrischer Netze (Vorlesungsnummer ETH Zrich 227-0526-00) lectures u given at ETH Zrich in Information Technology and Electrical Engineering. u In the lectures three main topics are covered, i.e. Power ow analysis Fault current calculations Power systems dynamics and stability In Part I of these notes the two rst items are covered, while Part II gives an introduction to dynamics and stability in power systems. In appendices brief overviews of phase-shifting transformers and power system protections are given. The notes start with a derivation and discussion of the models of the most common power system components to be used in the power ow analysis. A derivation of the power ow equations based on physical considerations is then given. The resulting non-linear equations are for realistic power systems of very large dimension and they have to be solved numerically. The most commonly used techniques for solving these equations are reviewed. The role of power ow analysis in power system planning, operation, and analysis is discussed.1 The next topic covered in these lecture notes is fault current calculations in power systems. A systematic approach to calculate fault currents in meshed, large power systems will be derived. The needed models will be given and the assumptions made when formulating these models discussed. It will be demonstrated that algebraic models can be used to calculate the dimensioning fault currents in a power system, and the mathematical analysis has similarities with the power ow analysis, so it is natural to put these two items in Part I of the notes. In Part II the dynamic behaviour of the power system during and after disturbances (faults) will be studied. The concept of power system stability is dened, and dierent types of power system instabilities are discussed. While the phenomena in Part I could be studied by algebraic equations,
As stated above, the presentation of the power ow problem is in these lecture notes based on physical considerations. However, for implementation of these equations in computer software, an object oriented approach is more suitable. Such an approach is presented in the lecture notes Modellierung und Analyse elektrischer Netze, by Rainer Bacher, is also to be dealt with in the lectures. In those lecture notes the techniques for handling large but sparse matrices are also reviewed.
1
vii
viii
Preface
the description of the power system dynamics requires models based on dierential equations. These lecture notes provide only a basic introduction to the topics above. To facilitate for readers who want to get a deeper knowledge of and insight into these problems, bibliographies are given in the text. I want to thank numerous assistants, PhD students, and collaborators of Power Systems Laboratory at ETH Zrich, who have contributed in various u ways to these lecture notes. Zrich in March 2006 u Gran Andersson o
Part I
Static Analysis
1
Introduction
This chapter gives a motivation why an algebraic model can be used to describe the power system in steady state. It is also motivated why an algebraic approach can be used to calculate fault currents in a power system.
POWER SYSTEM is predominantly in steady state operation or in a state that could with sucient accuracy be regarded as steady state. In a power system there are always small load changes, switching actions, and other transients occurring so that in a strict mathematical sense most of the variables are varying with the time. However, these variations are most of the time so small that an algebraic, i.e. not time varying model of the power system is justied. A short circuit in a power system is clearly not a steady state condition. Such an event can start a variety of dierent dynamic phenomena in the system, and to study these dynamic models are needed. However, when it comes to calculate the fault currents in the system, steady state (static) models with appropriate parameter values can be used. A fault current consists of two components, a transient part, and a steady state part, but since the transient part can be estimated form the steady state one, fault current analysis is commonly restricted to the calculation of the steady state fault currents.
1.1
Power Flow Analysis
It is of utmost importance to be able to calculate the voltages and currents that dierent parts of the power system are exposed to. This is essential not only in order to design the dierent power system components such as generators, lines, transformers, shunt elements, etc so that these can withstand the stresses they are exposed to during steady state operation without any risk for damages. Furthermore, for an economical operation of the system the losses should be kept at a low value taking various constraints into account, and the risk that the system enters into unstable modes of operation must be supervised. In order to do this in a satisfactory way the state of the system, i.e. all (complex) voltages of all nodes in the system, must be known. With these known, all currents, and hence all active and 1
1. Introduction
reactive power ows can be calculated, and other relevant quantities can be calculated in the system Generally the power ow, or load ow, problem is formulated as a nonlinear set of equations f (x, u, p) = 0 (1.1) where f is an n-dimensional (non-linear) function x is an n-dimensional vector containing the state variables, or states, as components. These are the unknown voltage magnitudes and voltage angles of nodes in the system u is a vector with (known) control outputs, e.g. voltages at generators with voltage control p is a vector with the parameters of the network components, e.g. line reactances and resistances The power ow problem consists in formulating the equations f in eq. (1.1) and then solving these with respect to x. This will be the subject dealt with in the rst part of these lectures. A necessary condition for eq. (1.1) to have a physically meaningful solution is that f and x have the same dimension, i.e. that we have the same number of unknowns as equations. But in the general case there is no unique solution, and there are are also cases when no solution exists. If the states x are known, all other system quantities of interest can be calculated from these and the known quantities, i.e. u and p. System quantities of interest are active and reactive power ows through lines and transformers, reactive power generation from synchronous machines, active and reactive power consumption by voltage dependent loads, etc. As mentioned above, the functions f are non-linear, which makes the equations harder to solve. In the solution of the equations the linearized equation f x = y (1.2) x is quite often used and solved. This equations gives also very useful informaf tion about the system. The Jacobian matrix , whose elements are given x by f fi = . (1.3) x ij xj can be used for many useful computations, and it is an important indicator of the system conditions. This will also be elaborated on.
1.2. Fault Current Analysis
1.2
Fault Current Analysis
In the lectures Elektrische Energiesysteme it was studied how to calculate fault currents, e.g. short circuit currents, for simple systems. This analysis will now be extended to deal with realistic systems including several generators, lines, loads, and other system components. Generator (synchronous machines) are important system components when calculating fault currents and their modelling will be elaborated on and discussed.
1.3
Literature
The material presented in these lectures constitutes only an introduction to the subject. Further studies can be recommended in the following text books: 1. Power Systems Analysis, second edition, by Artur R. Bergen and Vijay Vittal. (Prentice Hall Inc.,2000, ISBN 0-13-691990-1, 619 pages) 2. Computational Methods for Large Sparse Power Systems, An object oriented approach by S.A. Soma, S.A. Khaparde, Shubba Pandit (Kluwer Academic Publishers, 2002, ISBN 0-7923-7591-2, 333 pages)
1. Introduction
2
Network Models
In this chapter models of the most common network elements suitable for power ow analysis are derived. These models will be used in the subsequent chapters when formulating the power ow problem.
LL ANALYSIS in the engineering sciences starts with the formulation of appropriate models. A model, and in power system analysis we almost invariably then mean a mathematical model, is a set of equations or relations, which appropriately describes the interactions between dierent quantities in the time frame studied and with the desired accuracy of a physical or engineered component or system. Hence, depending on the purpose of the analysis dierent models of the same physical system or components might be valid. It is recalled that the general model of a transmission line was given by the telegraph equation, which is a partial dierential equation, and by assuming stationary sinusoidal conditions the long line equations, ordinary dierential equations, were obtained. By solving these equations and restricting the interest to the conditions at the ends of the lines, the lumped-circuit line models (-models) were obtained, which is an algebraic model. This gives us three dierent models each valid for dierent purposes. In principle, the complete telegraph equations could be used when studying the steady state conditions at the network nodes. The solution would then include the initial switching transients along the lines, and the steady state solution would then be the solution after the transients have decayed. However, such a solution would contain a lot more information than wanted and, furthermore, it would require a lot of computational eort. An algebraic formulation with the lumped-circuit line model, would give the same result with a much simpler model at a lower computational cost. In the above example it is quite obvious what model that is the appropriate one, but in many engineering studies the selection of the correct model is often the most dicult part of the study. It is good engineering practice to use as simple models as possible, but of course not too simple. If too complicated models are used, the analysis and computations would be unnecessarily cumbersome. Furthermore, generally more complicated models need more parameters for their denition and to get reliable values of these require often extensive work. 5
2. Network Models
i u
Rdx
Ldx Gdx
i+di Cdx u+du
dx
Figure 2.1. Equivalent circuit of a line element of length dx Figure 2.1. Equivalent circuit of a line element of length dx
In the subsequent sections algebraic models of the most common power system components suitable for power ow calculations will be derived. If not explicitly stated, symmetrical three-phase conditions are assumed in the following.
2.1
Lines and Cables
The equivalent -model of a transmission line section was derived in the lectures Elektrische Energiesysteme, 35-505. The general distributed model is characterized by the series parameters R = series resistance/km per phase (/km) X = series reactance/km per phase (/km) and the shunt parameters B = shunt susceptance/km per phase (siemens/km) G = shunt conductance/km per phase (siemens/km) as depicted in Figure 2.1. The parameters above are specic for the line or cable conguration and are dependent on conductors and geometrical arrangements. From the circuit in Figure 2.1 the telegraph equation is derived, and from this the lumped-circuit line model for symmetrical steady state conditions, Figure 2.2. This model is frequently referred to as the -model, and it is characterized by the parameters Zkm = Rkm + jXkm = series impedance ()
sh sh Ykm = Gsh + jBkm = shunt admittance (siemens) km 1
1 In Figure 2.2 the two shunt elements are assumed to be equal, which is true for homogenous lines,i.e. a line with equal values of the line parameters along its lengths, but this might not be true in the general case. In such a case the shunt elements are replaced sh sh sh sh by Ykm and Ymk with Ykm = Ymk with obvious notation.
2.1. Lines and Cables
I km
k
zkm
I mk
m
sh ykm
sh ykm
Figure 2.2. Lumped-circuit model (-model) of a transmission line between nodes k and m.
Note. In the following most analysis will be made in the p.u. system. For impedances and admittance capital letters indicate that the quantity is expressed in ohms or siemens, and lower case letters that they are expressed in p.u. Note. In these lecture notes complex quantities are not explicitly marked as underlined. This means that instead of writing z km we will write zkm when this quantity is complex. However, it should be clear from the context if a quantity is real or complex. Furthermore, we will not always use specic type settings for vectors. Quite often vectors will be denoted by bold face type setting, but not always. It should also be clear from the context if a quantity is a vector or a scalar. When formulating the network equations the node admittance matrix will be used and the series admittance of the line model is needed
1 ykm = zkm = gkm + jbkm
(2.1) (2.2)
with gkm = and
rkm 2 rkm + x2 km
2 rkm
bkm =
xkm + x2 km
(2.3)
For actual transmission lines the series reactance xkm and the series resistance rkm are both positive, and consequently gkm is positive and bkm sh sh is negative. The shunt susceptance ykm and the shunt conductance gkm are sh is so both positive for real line sections. In many cases the value of gkm small that it could be neglected.
2. Network Models
The complex currents Ikm and Imk in Figure 2.2 can be expressed as functions of the complex voltages at the branch terminal nodes k and m:
sh Ikm = ykm (Ek Em ) + ykmEk sh Imk = ykm(Em Ek ) + ykmEm
(2.4) (2.5)
where the complex voltages are
Ek = Uk ejk Em = Um ejm This can also be written in matrix form as Ikm Imk =
sh ykm + ykm ykm sh ykm ykm + ykm
(2.6) (2.7) Ek Em
(2.8)
As seen the matrix on the right hand side of eq. (2.8) is symmetric and the diagonal elements are equal. This reects that the lines and cables are symmetrical elements. Example 2.1. The series impedance of a 138 kV transmission line section is z = r + jx = 0.0062 + j0.0360 p.u. The total shunt susceptance (double the susceptance that appears in the equivalent -model) is bsh = 0.0105 p.u. and the shunt conductance is ignored. Calculate the series conductance and series susceptance and the ratio between the series and shunt susceptances. Solution The series conductance is given by g= r2 r 0.0062 = = 4.64 p.u. 2 +x 0.00622 + 0.03602
and the series susceptance by b= r2 x 0.0360 = = 27.0 p.u. 2 +x 0.00622 + 0.03602 27.0 = 2596 0.0105 0.036 = 5.8 0.0062
The b/bsh ratio is given by b/bsh = Of interest is also the ratio x/r x/r =
2.2. Transformers
Example 2.2. The series impedance and the total shunt impedance of a 750 kV line section are z = r + jx = 0.00072 + j0.0175 p.u. bsh = 8.77 p.u. Calculate the series conductance and series susceptance and the ratio between the series and shunt susceptances and the x/r ratio. Solution The series conductance and susceptance are given by g= 0.00072 = 2.35 p.u. 0.000722 + 0.01752 0.0175 = 57 p.u. 0.000722 + 0.01752
b=
The x/r ratio and b/bsh ratio are x/r = 0.0175 = 24.3 0.00072 57 = 6.5 8.77
b/bsh =
Note. The 750 kV line has a much higher x/r ratio than the 138 kV line and, at the same time, a much smaller (in magnitude) b/bsh ratio. (Why?) Higher x/r ratios mean better decoupling between active and reactive parts of the power ow problem, while smaller (in magnitude) b/bsh ratios may indicate the need for some sort of compensation, shunt or series, or both.
2.2
Transformers
We will start with a simplied model of a transformer where we neglect the magnetizing current and the no-load losses 2 . In this case the transformer can be modelled by an ideal transformer with turns ratio tkm in series with a series impedance zkm which represents resistive (load dependent) losses and the leakage reactance, see Figure 2.3. Depending on if tkm is real or non-real (complex) the transformer is in-phase or phase-shifting.
2 Often the magnetizing current and no-load losses are modelled by a shunt impedance, with much higher impedance than the leakage impedance. The inductive part of this impedance is then determined by the value of the magetizing current and the resistive part by the no load losses.
10
2. Network Models
Uk e jqk U jqk k ke k Ikm Ikm Uk e jqk U jqk k ke k Ikm
1:tkm 1:tkm
Up e jqp Up e jqp p p
zkm zkm
Um e jqm Um e jqm m Imk Imk m
tkm:1 tkm:1
zkm zkm
Um e jqm Um e jqm m Imk m
Imk Ikm Figure 2.3. Transformer model with complex ratio tkm = akm ejkm (tkm = a1 ejkm ) km Uk e jqk U jqk k ke k Ikm Ikm 1:akm 1:akm Up e jqp Up e jqp p p zkm zkm Um e jqm Um e jqm m Imk Imk m
Figure 2.4. In-phase transformer model
2.2.1
A In-Phase Transformers k m A Ikm mk Figure 2.4 shows an in-phase transformer model indicatingIthe voltage at the k m
internal - non-physical - node p. In this model the ideal voltage magnitude Imk Ikm ratio (turns ratio) is B C Up = akm (2.9) B C Uk Since k = p , this is also the ratio between the complex voltages at nodes k and p, Ep Up ejp = = akm (2.10) Ek Uk ejk There are no power losses (neither active nor reactive) in the ideal transformer (the k - p part of the model), which yields
Ek Ikm + Ep Imk = 0
(2.11)
Then applying eqs. (2.9) and (2.10) gives Ikm |Ikm | = = akm , Imk |Imk | (2.12)
Ikm
Imk
2.2. Transformers
11
A k Ikm B C Imk m
Figure 2.5. Equivalent model for in-phase transformer
which means that the complex currents Ikm and Imk are out of phase by 180 since akm R. Figure 2.5 represents the equivalent -model for the in-phase transformer in Figure 2.4. Parameters A, B, and C of this model can be obtained by identifying the coecients of the expressions for the complex currents Ikm and Imk associated with the models of Figures 2.4 and 2.5. Figure 2.4 gives Ikm = akm ykm (Em Ep ) = (a2 ykm )Ek + (akm ykm )Em km Imk = ykm (Em Ep ) = (akm ykm )Ek + (ykm )Em or in matrix form Ikm Imk = a2 ykm akm ykm km akm ykm ykm Ek Em (2.15) (2.13) (2.14)
As seen the matric on the right hand side of eq. (2.15) is symmetric, but the diagonal elements are not equal when a2 = 1. Figure 2.5 provides mow km the following: (2.16) Ikm = (A + B)Ek + (A)Em Imk = (A)Ek + (A + C)Em or in matrix form Ikm Imk = A+B A A A+C Ek Em (2.18) (2.17)
Identifying the matrix elements from the matrices in eqs. (2.15) and (2.18) yields A = akm ykm (2.19) B = akm (akm 1)ykm C = (1 akm )ykm (2.20) (2.21)
12
2. Network Models
Example 2.3. A 138/69 kV in-phase transformer with a series resistance of zero, a 0.23 p.u. series reactance, and p.u. turns ratio of 1 : 1.030 (from the model in Figure 2.4) Calculate the equivalent -model parameters. Solution B = akm (akm 1)ykm = 1.030(1.030 1)(j0.230)1 = j0.13 p.u. C = (1 akm )ykm = (1 1.030)(j0.230)1 = j0.13 p.u. Hence, since A, B, and C denote admittances, A and B are inductive, and C is capacitive. Example 2.4. A 500/750 kV in-phase transformer with a series resistance of zero, a 0.00623 p.u. series reactance, and p.u. turns ratio of 1 : 0.950 (from the model in Figure 2.4) Calculate the parameters for the equivalent -model. (100 MVA base) Solution As in the previous example one obtains: A = j152.5 p.u. B = j7.62 p.u. C = j8.03 p.u. i.e., parameter B is capacitive and parameter C is inductive (762 Mvar and -803 Mvar, respectively, assuming nominal voltage magnitudes and a 100 MVA base). A = akm ykm = 1.030(j0.230)1 = j4.48 p.u.
2.2.2
Phase-Shifting Transformers
Phase-shifting transformers, such as the one represented in Figure 2.6 , are used to control active power ows; the control variable is the phase angle and the controlled quantity can be, among other possibilities, the active power ow in the branch where the shifter is placed. In Appendix A the physical design of phase-shifting transformer is described. A phase-shifting transformer aects both the phase and magnitude of complex voltages Ek and Ep , without changing their ratio, i.e., Ep = tkm = akm ejkm Ek Thus, p = k + km and Up = akm Uk , using eqs. (2.11) and (2.22), Ikm = t = akm ejkm km Imk (2.23) (2.22)
2.2. Transformers
13
Uk e jqk k Ikm
1:akme jjkm p qp = qk + jkm Up = akmUk
zkm
Um e jqm m Imk
Figure 2.6. Phase-shifting transformer with tkm = akm ejkm .
Uq e jqk Um e jqm zkm As with in-phase transformers, the complex currents Ikm and Imk can k be expressed in terms of complex voltages at the shifter terminal nodes: m q p Imk Ikm 1:tkm t 1 2 (2.24) I = t y (E E ) = (a y )E + (tmk:y )E Uk e jqk
km km km m p km km k km km m
Up e jqp
Imk = ykm(Em Ep ) = (tkm ykm)Ek + (ykm )Em or in matrix form Ikm Imk = a2 ykm t ykm km km tkm ykm ykm Ek Em
(2.25)
(2.26)
m p As seen thisk matrix is not symmetric if tkm is qnon-real, :1 tmk and the diagonal 1:tkm q matrix elements are not equal p a2 =km There is no way to determine if km z 1. k A, B, and C of the equivalent -model from these equations, m parameters y since the coecient tkm km of Em in eq. (2.24) diers from tkm ykm Ikm Imk in eq. (2.25), as long as there is nonzero phase shift, i.e. tkm R. A / phase-shifting transformer can thus not be represented by a -model. sh sh
U e jqk
U e jqp
U e jqq
U e jqm
ykm
ymk
Example 2.5. A 230/138 kV transformer (Figure 2.6) has a series resistance of zero, a 0.0127 p.u. series reactance, and a complex turns ratio of 1 : 1.007ej30 (Y- connection). Show that this transformer can be seen to consist of a series connection of two transformers: an ideal in-phase transformer with a turns ratio of 1 : 1.007 (constant voltage phase) and a phase-shifting transformer with a complex turns ratio of 1 : ej30 (constant voltage amplitude) and a series reactance of 0.0127 p.u. Note. If no parallel paths exist, the phase-shifting has no signicance. The introduced phase-shift can in such a case been seen as a shift of the phase angle of the reference node. Y- connected transformers are often used to provide zero-sequence de-coupling between two parts of the system, and not
14
2. Network Models
1:1 XtA = 0.2 p.u. A 1:t12 XtB = 0.4 p.u B E2 Load Iload
Figure 2.7. System for Example 2.6.
for active power ow. For active power ow control usually phase-shifting much lower than 30 is needed. Often the phase-shifting could be varied to cope with dierent loading situations in the system. Example 2.6. Consider two parallel connected transformers according to Figure 2.7. Transformer A has a turns-ratio of 1:1 (p.u./p.u.) while the turns-ratio of transformer B will be varied as described below. The transformers are feeding a load at bus 2, Iload = 1.05 45 p.u. and E2 = 10 p.u. The reactances of the transformers are given in the gure. Calculate the complex power through the transformers when the turnsratios of transformer B are (t12 = a12 ej12 ) (a) a12 = 1, 12 = 0 (b) a12 = 1.05, 12 = 0 (c) a12 = 1, 12 = 3 Comments!
2.2.3
Unied Branch Model
The expressions for the complex currents Ikm and Imk for both transformers and shifters derived above, depend on which side the tap is located; i.e., they are not symmetrical. It is however possible, to develop unied complex expressions which can be used for lines, transformers, and phase-shifters, regardless of the side on which the tap is located (or even in the case when there are taps on both sides of the device). Consider initially the model
2.2. Transformers
15
in Figure 2.8 in which shunt elements have been temporarily ignored and tkm = akm ejkm and tmk = amk ejmk . In this case Ikm = t Ipq = t (Ep Eq )ykm = t (tkm Ek tmk Em )ykm km km km and Imk = t Iqp = t (Eq Ep )ykm = t (tmk Em tkm Ek )ykm mk mk mk which together yield Ikm = (a2 Ek t tmk Em )ykm km km and Imk = (a2 Em t tkm Ek )ykm mk mk (2.30) (This expressions are symmetrical in the sense that if k and m are interchanged, as in the expression for Ikm , the result is the expression for Imk , and vice-versa.) Figure 2.9 shows the unied branch model. All the devices studied above can be derived from this general model by establishing the appropriate definitions of the parameters that appear in the unied model. Thus, for instance, if tkm = tmk = 1 is assumed, the result is an equivalent -model of a transmission line; or, if the shunt elements are ignored, and tkm = 1 and tmk = amk ejmk is assumed, then the result is a phase shifting transformer with the tap located on the bus m side. The general expressions for and Imk can be obtained from the model in Figure 2.9:
sh Ikm = (a2 Ek t tmk Em )ykm + ykm a2 Ek km km km
(2.27)
(2.28)
(2.29)
(2.31) (2.32)
and
sh Imk = (a2 Em t tkm Ek )ykm + ymk a2 Em mk mk mk
or Ikm Imk =
sh a2 (ykm + ykm ) t tmk ykm km km sh t tkm ykm a2 (ykm + ykm ) mk mk
Ek Em
(2.33)
The transformer modelled above were all two-winding transformers. In power systems there are also three-windings, and N -windings N > 3, transformers, and they can be modelled in a similar way. Instead of linear relationships between two complex currents and two complex voltages, we will have linear relationships between these quantities modelled by N N matrices. It should be notes that for N -winding transformers there are leakage reactances between each pair of windings, i.e. in total N (N 1)/2 reactances. Since the power ratings of the dierent windings are not necessarily equal, as in the two-winding case, it should be noted that the power bases for expressing the reactances in p.u. are not always equal.
p k Up = akmUkm k Up = akmUk
16
2. Network Models
Uk e jqk U jq k ke k k Ikm Ikm 1:tkm 1:tkm
zkm zkm
Uq e jqk Uq e jqk q q tmk:1 tmk:1
Figure 2.8. Transformer symmetrical model.
Uk e jqk Uk e jqk 1:tkm k 1:tkm Ikm k Ikm
Up e jqp U p jqp e
p
Uq e jqq zkm zkm
t :1 Uqq jqq mk e tmk:1 q
sh ykm sh ykm
sh ymk sh ymk
Figure 2.9. Unied branch model extended (-model).
2.3
Shunt Elements
The modelling of shunt elements in the network equations is straightforward and the main purpose here is to introduce the notation and the sign convention to be used when formulating the network equations in the coming chapters. As seen from Figure 2.10 the current from a shunt is dened as positive when injected into the bus. This means that
sh sh Ik = yk Ek
(2.34)
with Ek being the complex voltage at node k. Shunts are in all practical cases either shunt capacitors or reactors. From eq. (2.34) the injected complex power is 2 S sh = P sh + jQsh = y sh |Ek |2 = y sh Uk (2.35) k k k k k
2.4
Loads
Load modelling is an important topic in power system analysis. When formulating the load ow equations for high voltage systems, a load is most
2.5. Generators
17
Ik k
sh
Ik
Ikm
sh yk
sh yk
Ik
sh
Figure 2.10. A shunt connected to bus k.
Ik
load
k often the infeed of power to a network at a lower voltage level, e.g. a distribution network. Often the voltage in the distribution systems is kept constant by controlling the tap-positions of the distribution transformers which means that power, active and reactive, in most cases can be regarded as independent of the voltage on the high voltage side. This means that the Load load complex power Ek (Ik ) is constant, i.e. independent of the voltage magnitude Uk . Also in this case the current is dened as positive when injected into the bus, see Figure 2.11. In the general case the complex load current can be written as load load Ik = Ik (Uk ) (2.36)
load where the function Ik () describes the load characteristics 3 . More often gen the load characteristics are given for the active and reactive powers Ik
load load Pk = Pk (Uk )
(2.37) (2.38)
Qload = Qload (Uk ) k k
Generator
2.5
Generators
Generators are in load ow analysis modelled as current injections, see Figure 2.12. In steady state a generator is commonly controlled so that the active power injected into the bus and the voltage at the generator terminals are kept constant. This will be elaborated later when formulating the load ow equations. Active power from the generator is determined by the turbine control and must of course be within the capability of the turbinegenerator system. Voltage is primarily determined by the reactive power
This refers to the steady state model of the load. For transient conditions other load models apply. These are usually formulated as dierential equations and might also involve the frequency
3
sh yk
18
2. Network Models
Ik k k Ik
load
load
Load Load
Figure 2.11. Model of a load connected to bus k. gen Ik k I gen k
k Generator Generator
Figure 2.12. Model of a generator connected to bus k.
injection into the node, and since the generator must operate within its reactive capability curve it is not possible to control the voltage outside certain limits. The reactive capability of a generator depends on a number of quantities, such as active power, bus voltage and other operating conditions, and a typical example is shown in Figure 2.13. The shape of the generator capability curve is specic for each generator and depends on design characteristics, type of generator, hydro or steam turbine, stability constraints, etc. In Figure 2.13 it is also indicated what imposes the dierent limits for this particular generator. These are briey discussed below. These limits are also discussed in Chapter 9.
2.5.1
Stator Current Heating Limit
2 The losses in the armature windings are given by Rt It , with obvious notation. These losses result in a temperature rise in the armature windings, and this must be limited to a given value otherwise the generator is
2.5. Generators
19
Field current heating limit
Overexcited
Stator current heating limit
P
Underexcited
End region heating limit

Figure 2.13. Reactive capability curve of a turbo generator.
damaged or its life time is reduced. Since the complex power is given by S = P + jQ = Ut It it means that for a given terminal voltage, Ut , circles in the P Q-plane with centre at the origin correspond to constant value of the magnitude of the armature current It . The stator current limit for a given terminal voltage is thus a circle with the centre at the origin. At high loading of the generator, this is usually determining the reactive capability of the synchronous machine.
2.5.2
Field Current Heating Limit
The reactive power that can be generated at low load is determined by the eld current heating limit. It can shown that the locus for constant eld 2 current is a circle with the centre on the Q-axis at Et /Xs , where Et is the terminal voltage and Xs is the synchronous reactance. The radius is given by machine parameters and typical behaviour is shown in Figure 2.13. The eld current heating is usually limiting at overexcited operation at low load.
20
2. Network Models
2.5.3
Stator End Region Heating Limit
When the synchronous machine is underexcited the armature end leakage ux is increased. This ux enters and leaves in a direction perpendicular to the stator laminations causing eddy currents in the laminations, and hence heating. This can limit the capability, particularly for round rotor machines. A typical example of such a limitation is shown in Figure 2.13. The limits on reactive power generation and absorption given in this subsection and the two previous ones are imposed by internal synchronous machine design. As seen they are dependent on the terminal voltage of the generator. There are also other considerations that can further limit the reactive capability of the synchronous machine. These are stability limits that could further limit the operation of the machine in underexcited operation.
3
Active and Reactive Power Flows
In this chapter the expressions for the active and reactive power ows in transmission lines, transformers, phase-shifting transformers, and unied branch models are derived.
HE SYSTEM COMPONENTS dealt with in this chapter are linear in the sense that the relations between voltages and currents are linear1 . However, since one usually is interested rather in powers, active and reactive, than currents, the resulting equations will be non-linear, which introduces a complication when solving the resulting equations.
T
3.1
Transmission Lines
Consider the complex current Ikm in a transmission line Ikm = ykm (Ek Em ) + jbsh Ek km (3.1)
with quantities dened according to Figure 2.2. The complex power, Skm = Pkm + jQkm , is
2 Skm = Ek Ikm = ykm Uk ejk (Uk ejk Um ejm ) jbsh Uk km
(3.2)
sh where the conductance of ykm has been neglected. The expressions for Pkm and Qkm can be determined by identifying the corresponding coecients of the real and imaginary parts of eq. (3.2), which yields 2 Pkm = Uk gkm Uk Um gkm cos km Uk Um bkm sin km
(3.3) (3.4)
2 Qkm = Uk (bkm + bsh ) + Uk Um bkm cos km Uk Um gkm sin km km
where the notation km = k m is introduced.

1 This is at least true for the models analysed here. Dierent non-linear phenomena, e.g. magnetic saturation, can sometimes be important, but when studying steady state conditions the devices to be discussed in this chapter are normally within the region of linearity.
21
22
3. Active and Reactive Power Flows
The active and reactive power ows in opposite directions, Pmk and Pmk , can be obtained in the same way, resulting in:
2 Pmk = Um gkm Uk Um gkm cos km + Uk Um bkm sin km 2 Um (bkm
(3.5) (3.6)
Qmk =
+ bsh ) + km
Uk Um bkm cos km + Uk Um gkm sin km
From these expressions the active and reactive power losses of the lines are easily obtained as:
2 2 Pkm + Pmk = gkm (Uk + Um 2Uk Um cos km )
= gkm |Ek Em |2
(3.7)
2 2 2 2 Qkm + Qmk = bsh (Uk + Um ) bkm (Uk + Um 2Uk Um cos km ) km 2 2 = bsh (Uk + Um ) bkm |Ek Em |2 km
(3.8)
Note that |Ek Em | represents the magnitude of the voltage drop across the line, gkm |Ek Em |2 represents the active power losses, bsh |Ek Em |2 km 2 2 represents the reactive power losses; and bsh (Uk + Um ) represents the km reactive power generated by the shunt elements of the equivalent -model (assuming actual transmission line sections, i.e. with bkm < 0 and bsh > 0). km Example 3.1. A 750 kV transmission line section has a series impedance of 0.00072 + j0.0175 p.u., a total shunt impedance of 8.775 p.u., a voltage magnitude at the terminal buses of 0.984 p.u. and 0.962 p.u., and a voltage angle dierence of 22 . Calculate the active and reactive power ows. Solution The active and reactive power ows in the line are obtained by applying eqs. (3.3) and (3.4), where Uk = 0.984 p.u., Um = 0.962 p.u., and km = 22 . The series impedance and admittances are as follows: zkm = 0.00072 + j0.0175 p.u.
1 ykm = gkm + jbkm = zkm = 2.347 j57.05 p.u.
The -model shunt admittances (100 MVA base) are: bsh = 8.775/2 = 4.387 p.u. km and Pkm = 0.9842 2.3470.9840.9622.347 cos 22 +00.9840.96257.05 sin 22 p.u. Qkm = 0.9842 (57.05+4.39)0.9840.96257.05 cos 22 00.9840.9622.347 sin 22 p.u. which yield Pkm = 2044 MW Qkm = 8.5 Mvar
3.2. In-phase Transformers In similar way one obtains: Pmk = 2012 MW Qmk = 50.5 Mvar
23
It should be noted that powers are positive when injected into the line.
3.2
In-phase Transformers
The complex current Ikm in an in-phase transformer is expressed as in eq. (2.13) Ikm = akm ykm(akm Ek Em ) The complex power, Skm = Pkm + jQkm , is given by
Skm = Ek Ikm = ykm akm Uk ejk (akm Uk ejk Um ejm )
(3.9)
Separating the real and imaginary parts of this latter expression yields the active and reactive power ow equations: Pkm = (akm Uk )2 gkm akm Uk Um gkm cos km akm Uk Um bkm sin km (3.10) Qkm = (akm Uk )2 bkm + akm Uk Um bkm cos km akm Uk Um gkm sin km (3.11) These same expressions can be obtained by comparing eqs. (3.9) and 2 (3.2); in eq. (3.9) the term jbsh Uk is not present, and Uk is replaced by km akm Uk . Hence, the expressions for the active and reactive power ows on in-phase transformers are the same expressions derived for a transmission line, except the for two modications: ignore bsh , and replace Uk with km akm Uk . Example 3.2. A 500/750 kV transformer with a tap ratio of 1.050:1.0 on the 500 kV side, see Figure 2.4, has neglible series resistance and a leakage reactance of 0.00623 p.u., terminal voltage magnitudes of 1.023 p.u. and 0.968 p.u., and an angle spread of 5.3 . Calculate the active and reactive power ows in the transformer. Solution The active and reactive power ows in the transformer are given by eqs. (3.10) and (3.11), where Uk = 1.023 p.u., Um = 0.968 p.u., km = 5.3 , and akm = 1.0/1.05 = 0.9524. The series reactance and susceptance are as follows: xkm = 0.00623 p.u. bkm = x1 = 160.51 p.u. km
24
The active and reactive power ows can be expressed as Pkm = 0.9524 1.023 0.968 (160.51) sin 5.3 p.u. Qkm = (0.95241.023)2 (160.51)+0.95241.0230.968(160.51) cos 5.3 p.u. which yield Pkm = 1398 MW Qkm = 163 Mvar The reader is encouraged to calculate Pmk and Qmk . (The value of Pmk should be obvious.)
3.3
Phase-Shifting Transformer with akm = 1
The complex current Ikm in a phase shifting transformer with akm = 1 is as follows, see Figure 2.6: Ikm = ykm (Ek ejkm Em ) = ykm ejkm (Ek ejkm Em ) and the complex power, Skm = Pkm + jQkm , is thus
Skm = Ek Ikm = ykm Uk ej(k +km ) (Uk ej(k +km ) Um ejm )
(3.12)
(3.13)
Separating the real and imaginary parts of this expression, yields the active and reactive power ow equations, respectively:
2 Pkm = Uk gkm Uk Um gkm cos(km + km )
Uk Um bkm sin(km + km )
(3.14)
2 Qkm = Uk bkm + Uk Um bkm cos(km + km )
Uk Um gkm sin(km + km )
(3.15)
As with in-phase transformers, these expressions could have been obtained through inspection by comparing eqs. (3.2) and (3.13): in eq. (3.13), 2 the term jbsh Uk is not present, and km is replaced with km + km . Hence, km the expressions for the active and reactive power ows in phase-shifting transformers are the same expressions derived for the transmission line, albeit with two modiactions: ignore bsh and replace km with km + km . km Example 3.3. A - Y , 230/138 kV transformer presents a 30 phase angle shift. Series resistance is neglected and series reactance is 0.0997 p.u. Terminal voltage magnitudes are 0.882 p.u. and 0.989 p.u., and the total angle dierence is 16.6 . Calculate the active and reactive power ows in the transformer.
3.4. Unied Power Flow Equations
25
Solution The active and reactive power ows in the phase-shifting transformer are given by eqs. (3.14) and (3.15), where Uk = 0.882 p.u., Um = 0.989 p.u., km = 16.6 , and km = 30 . The series reactance and susceptance are as follows: xkm = 0.0997 p.u. bkm = x1 = 10.03 p.u. km The active and reactive power ows can be expressed as Pkm = 0.882 0.989 (10.03) (160.51) sin(16.6 + 30 ) p.u. Qkm = 0.8822 (10.03) + 0.882 0.989 (10.03) cos(16.6 + 30 ) p.u. which yield Pkm = 203 MW Qkm = 70.8 Mvar The reader is encouraged to calculate Pmk and Qmk . (The value of Pmk should be obvious.)
3.4
Unied Power Flow Equations
The expressions for active and reactive power ows on transmission lines, in-phase transformers, and phase shifting transformers, see Figure 2.9, can be expressed in the following unied forms: Pkm = (akm Uk )2 gkm (akm Uk )(amk Um )bkm sin(km + km mk ) Qkm = (akm Uk )2 (bkm + bsh ) km + (akm Uk )(amk Um )bkm cos(km + km mk ) (akm Uk )(amk Um )gkm sin(km + km mk ) (3.17) (akm Uk )(amk Um )gkm cos(km + km mk ) (3.16)
Where, for the transmission lines like the one represented in Figure 2.2, akm = amk = 1 and km = mk = 0; for in-phase transformers such as the sh sh one represented in Figure 2.4, ykm = ymk = 0, amk = 1 and km = mk = 0; sh and for a phase-shifting transformer such as the one in Figure 2.6, ykm = sh = 0, a ymk mk = 1 and mk = 0.
26
4
Nodal Formulation of the Network Equations
In this chapter the basic network equations are derived from Kirchho s Current Law (KCL) and put into forms that are suitable for the formulation of the power ow equations in the subsequent chapter
HE NET COMPLEX current injection at a network bus, see Figure 4.1, is related to the current ows in the branches incident to the bus. Applying Kirchhos Current Law (KCL) yields
sh Ik + Ik =
from generators where k is a generic node, sh and loads, Ik is the current injection from shunts, m is a node adjacent to k, k is the set of nodes adjacent to k, and N is the number of nodes in the network. The complex current Ikm in the unied branch model, Figure 2.9, is
sh Ikm = (a2 Ek t tmk Em )ykm + ykm a2 Ek km km km
Ikm ; for k = 1, . . . , N mk Ik is the net current injection
(4.1)
(4.2)
where tkm = akm ejkm and tmk = amk ejmk .
Ik k
sh
Ik
Ikm
sh yk
sh yk
Ik
sh
Figure 4.1. Generic bus with sign conventions for currents and power ows.
27
Ik k
load
Load
28
4. Nodal Formulation of the Network Equations Equations (4.1) and (4.2) yield
sh Ik = yk + mk sh a2 (ykm + ykm ) Ek km
(t tmk ykm )Em km

mk
(4.3)
for k = 1, . . . , N . This expression can be written as I = YE where I is the injection vector with elements Ik , k = 1, . . . , N E is the nodal voltage vector with elements Ek = Uk ejk Y = G + jB is the nodal admittance matrix, with the following elements (4.5) Ykm = t tmk ykm km
sh Ykk = yk + mk sh a2 (ykm + ykm ) km
(4.4)
(4.6)
We see that the nodal admittance matrix dened by eqs. (4.5) and (4.6) is modied as compared with the nodal admittance matrix without transformers. Particularly it should be noted that Y as dened above is not necessarily symmetric. For large practical networks this matrix is usually very sparse. The degree of sparsity (percentage of zero elements) normally increases with the dimensions of the network: e.g., a network with 1000 buses and 1500 branches, typically presents a degree of sparsity greater than 99 %, i.e. less than 1 % of the matrix elements have non-zero values. The kth component of I, Ik , dened in eq. (4.3)can by using eqs. (4.5) and (4.6) be written as Ik = Ykk Ek +
mk
Ykm Em =
mK
Ykm Em
(4.7)
where K is the set of buses adjacent to bus k, including bus k, and k is the set of buses adjacent to bus k, excluding bus k. Now considering that Ykm = Gkm + jBkm and Em = Um ejm , eq. (4.7) can be rewritten as Ik =
mK
(Gkm + jBkm )(Um ejm )
(4.8)
The complex power injection at bus k is

Sk = Pk + jQk = Ek Ik
(4.9)
29 and by applying eqs. (4.8) and (4.9) this gives Sk = Uk ejk

mK
(Gkm jBkm )(Um ejm )
(4.10)
The expressions for active and reactive power injections are obtained by identifying the real and imaginary parts of eq. (4.10), yielding Pk = Uk
mK
Um (Gkm cos km + Bkm sin km) Um (Gkm sin km Bkm cos km )
(4.11)
Qk = Uk
mK
(4.12)
30
4. Nodal Formulation of the Network Equations
5
Basic Power Flow Problem
In this chapter the basic power ow problem is formulated and the basic bus types are dened. Also, the conditions for solvability of the problem are discussed
HE POWER FLOW PROBLEM can be formulated as a set of nonlinear algebraic equality/inequality constraints. These constraints represent both Kirchhos laws and network operation limits. In the basic formulation of the power ow problem, four variables are associated to each bus (network node) k: Uk - voltage magnitude k - voltage angle Pk - net active power (algebraic sum of generation and load) Qk - net reactive power (algebraic sum of generation and load)
5.1
Basic Bus Types
Depending on which of the above four variables are known (given) and which ones are unknown (to be calculated), two basic types of buses can be dened: PQ bus: Pk and Qk are specied; Uk and k are calculated PU bus: Pk and Uk are specied; Qk and k are calculated PQ buses are normally used to represent load buses without voltage control, and PU buses are used to represent generation buses with voltage control in power ow calculations1 . Synchronous compensators2 , are also treated as PU buses. A third bus is also needed:
Synchronous machines are often equipped with Automatic Voltage Regulators (AVRs), which controls the excitation of the machine so that the terminal voltage, or some other voltage close to the machine, is kept at the set value. 2 Synchronous compensators, sometimes also called synchronous condensers, are synchronous machines without any active power generation or load (except for losses) used for reactive power and voltage control.
1
31
32
5. Basic Power Flow Problem U bus: Uk and k are specied; Pk and Qk are calculated
The U bus, also called reference bus or slack bus, has double functions in the basic formulation of the power ow problem: 1. It serves as the voltage angle reference 2. Since the active power losses are unknown in advance, the active power generation of the U bus is used to balance generation, load, and losses In normal power systems PQ-buses or load buses are the far most common, typically comprising more than 80% of all buses. Other possible bus types are P, U, and PQU, with obvious denitions. The use of multiple U buses may also be required for certain applications. In more general cases, the given values are not limited to the specic set of buses (P, Q, U, ), and branch related variables can also be specied. Example 5.1. Figure 5.1 shows a 5-bus network with four transmission lines and two transformers. Generators, with voltage control, are connected at buses 1, 3, and 5, and loads are connected at buses 4 and 5, and at bus 4 a shunt is also connected. Classify the buses according to the bus types PU, PQ and U. Solution Buses 1,3, and 5 are all candidates for PU or U bus types. Since only one could be U bus, we select (arbitrarily) bus 5 as U. In a practical system usually a generator, or generator station, that could produce power within a large range is selected as reference or slack bus. It should be noted that even if a load is connected to bus 5 it can only be a PU or U bus, since voltage control is available at the bus. The reference angle is set at bus 5, usually to 0. Bus 2 is a transition bus in which both P and Q are equal to zero, and this bus is consequently of type PQ. Bus 4 is a load bus to which is also connected a shunt susceptance: since shunts are modelled as part of the network, see next section, the bus is also classied as a PQ bus.
5.2
Equality and Inequality Constraints
Eqs. (4.11) and (4.12) can be rewritten as follows Pk =

mk
Pkm (Uk , Um , k , m ) Qkm (Uk , Um , k , m )

mk
(5.1)
Qk + Qsh (Uk ) = k where
(5.2)
Pkm = (qk - q m ) x km
5.2. Equality and Inequality Constraints
33
PU 1
PQ 2
PU 3
Uq
PQ
Figure 5.1. 5-bus system
k = 1, . . . , N (N is the number of buses in the network) k : set of buses adjacent to k Uk , Um : voltage magnitudes at the terminal buses of branch k m k , m : voltage angles at the terminal buses of branch k m Pkm : active power ow from bus k to bus m Qkm : reactive power ow from bus k to bus m Qsh : component of reactive power injection due to the shunt element k 2 at bus k (Qsh = bsh Um ), where bsh is the shunt impedance3 k k k A set of inequality constraints imposes operating limits on variables such as the reactive power injections at PU buses (generator buses), se section 2.5, and voltage magnitudes at PQ buses: Qmin Qk Qmax k k
min max Uk Uk Uk
3
(5.3) (5.4)
It is here assumed that all shunts are reactive without losses. If shunts with resistive components should be included, then eq. (5.1) must be modied accordingly.
34
5. Basic Power Flow Problem
When no inequality constraints are violated, nothing is aected in the power ow equations, but if a limit is violated, the bus status is changed and it is enforced as an equality constraint at the limiting value. This normally requires a change in bus type: if, for example, a Q limit of a PU bus is violated, the bus is transformed into an PQ bus (Q is specied and the U becomes a problem unknown). A similar procedure is adopted for backing-o when ever appropriate. What is crucial is that bus type changes must not aect solvability. Various other types of limits are also considered in practical implementations, including branch current ows, branch power ows, active power generation levels, transformer taps, phase shifter angles, and area interchanges.
5.3
Problem Solvability
One problem in the denition of bus type (bus classication) is to guarantee that the resulting set of power ow equations contains the same number of equations as unknowns, as are normally necessary for solvability, although not always sucient. Consider a system with N buses, where NP U are of type PU, NP Q are of type PQ, and one is of type U. To fully specify the state of the system we need to know the voltage magnitudes and voltage angles of all buses, i.e. in total 2N quantities. But the voltage angle and voltage magnitude of the slack bus are given together with the voltage magnitudes of NP U buses. Unknown are thus the voltage magnitudes of the PQ buses, and the voltage angles of the PU and the PQ buses, giving a total of NP U + 2NP Q unknown states. From the PU buses we get NP U balance equations regarding active power injections, and from the PQ buses 2NP Q equations regarding active and reactive power injections, thus in total NP U + 2NP Q equations, and hence equal to the number of unknowns, and the necessary condition for solvability has been established. Similar necessary conditions for solvability can be established when other types of buses, such as P, U, and PQU buses are used in the formulation of the power ow problem. Example 5.2. Consider again the 5-bus in Figure 5.1. Formulate the equality constraints of the system and the inequality constraints for the generator buses.
Solution In this case N = 5, NP Q = 2, NP U = 2, and of course NU = 1. The number of equations are thus: NP U + 2NP Q = 2 + 2 2 = 6, and these
5.3. Problem Solvability are: P1 = P12 + P15 P2 = P21 + P23 + P25 Q2 = Q21 + Q23 + Q25 P3 = P32 + P34 P4 = P43 + P45 Q4 + Qsh = Q43 + Q45 4
35
In the above equations P1 , P2 , P3 , P4 , Q2 , and Q4 are given. All the other quantities are functions of the bus voltage magnitudes and phase angles, of which U1 , U3 , and U5 and 5 are given. The other six, i.e. U2 , U4 , 1 , 2 , 3 , 4 , in total 6 unknowns, can be solved from the above equations, and from these all power ows and injections can be calculated. The inequality constraints of the generator buses are: Qmin Q1 Qmax 1 1 Qmin Q3 Qmax 3 3 Qmin Q5 Qmax 5 5 The reactive limits above are derived from the generator capability curves as explained in section 2.5. For the slack bus it must also be checked that the injected active and reactive powers are within the range of the generator, if not the power generation of the other generators must be changed or the voltage settings of these.
36
5. Basic Power Flow Problem
6
Solution of the Power Flow Problem
In this chapter the basic methods to solve the non-linear power ow equations are reviewed. Solution methods based on the observation that active and reactive power ows are not so strongly coupled are introduced.
N ALL REALISTIC CASES the power ow problem cannot be solved analytically, and hence iterative solutions implemented in computers must be used. In this chapter we will review two solutions methods, Gauss iteration with a variant called Gauss-Seidel iterative method, and the NewtonRaphson method.
6.1
Solution by Gauss-Seidel Iteration
Consider the power ow equations (5.1) and (5.2) which could be written in complex form as Sk = Ek
mK Ykm Em
k = 1, 2, . . . N
(6.1)
which is a the same as eq. (4.10). The set K is the set of buses adjacent (connected) to bus k, including bus k, and hence shunt admittances are included in the summation. Furthermore Ek = Uk ejk . This equation can be rewritten as
Ek =
1 Sk Ykk Ek
Ykm Em mk
k = 1, 2, . . . N
(6.2)
where k is the set of all buses connected to bus k excluding bus k. Taking the complex conjugate of eq. (6.2) yields
1 Sk Ykk Ek
Ek =
Ykm Em
mk
k = 1, 2, . . . N
(6.3)
37
38
6. Solution of the Power Flow Problem
Thus we get N 1 algebraic (complex) equations in the complex variables Ek in the form E2 = h2 (E1 , E2 , . . . , EN ) E3 = h3 (E1 , E2 , . . . , EN ) . . . EN = hN (E1 , E2 , . . . , EN ) where the functions hi are given by eq. (6.3). It is here assumed that bus number 1 is the U bus, and hence is E1 given and we have no equation for node 1. For PQ buses both the magnitude and angle of Ek are unknown, while for PU buses only the angle is unknown. For PQ buses Sk is known, while for PU buses only Pk is known. This will be discussed below in more detail. In vector form eq. (6.4) can be written as x = h(x) and based on this equation the following iterative scheme is proposed x+1 = h(x ) = 0, 1, . . . (6.6) where the superscript indicates the iteration number. Thus starting with an initial value x0 , the sequence x0 , x1 , x2 , . . . is generated. If the sequence converges, i.e. x x , then x = h(x ) (6.8) and is a solution of eq. (6.5). In practice the iteration is stopped when the changes in x becomes suciently small, i.e. when the norm of x = x+1 x is less than a pre-determined value . To start the iteration a rst guess of x is needed. Usually, if no a priori knowledge of the solution is known, one selects all unknown voltage magnitudes and phase angles equal to the ones of the reference bus, usually around 1 p.u. and phase angle = 0. This start solution is often called a at start. The dierence between Gauss and Gauss-Seidel iteration can be explained by considering eq. (6.6) with all components written out explicitly1 x+1 = h2 (x1 , x , . . . , x ) 2 N 2 x+1 = h3 (x1 , x , . . . , x ) 2 N 3 . . . x+1 = hN (x1 , x , . . . , x ) 2 N N
In this particular formulation x1 is the value of the complex voltage of the slack bus and consequently known. For completeness we have included it as a variable in the equations above, but it is actually known
1
(6.4)
(6.5)
(6.7)
(6.9)
6.2. Newton-Raphson Method
39
In carrying out the computation (normally by computer) we process the equations from top to bottom. We now observe that when we solve for x+1 3 we already know x+1 . Since x+1 is presumably a better estimate than x , 2 2 2 it seems reasonable to use the updated value. Similarly when we solve for x+1 we can use the values of x+1 and x+1 . This is the line of reasoning 4 2 3 called the Gauss-Seidel iteration: x+1 = h2 (x1 , x , . . . , x ) 2 N 2 x+1 = h3 (x1 , x+1 , . . . , x ) N 3 2 . . . x+1 = hN (x1 , x+1 , . . . , x+1 , x ) 2 N N 1 N It is clear that the convergence of the Gauss-Seidel iteration is faster than the Gauss iteration scheme. For PQ buses the complex power Sk is completely known and the calculation of the right hand side of eq. (6.3) is well dened. For PU buses however, Q is not dened but is determined so that the voltage magnitude is kept at the specied value. In this case we have to estimate the reactive power injection and an obvious choice is
Q = Ek k Ykm (Em ) mK
(6.10)
(6.11)
In the Gauss-Seidel iteration scheme one should use the latest calculated values of Em . It should be clear that also for PU buses the above iteration scheme gives an solution if it converges. A problem with the Gauss and Gauss-Seidel iteration schemes is that convergence can be very slow, and sometimes even the iteration does not converge despite that a solution exists. Furthermore, no general results are known concerning the the convergence characteristics and criteria. Therefore more ecient solution methods are needed, and one such method that is widely used in power ow computations is discussed in the subsequent sections.
6.2
Newton-Raphson Method
Before applying this method to the power ow problem we review the iteration scheme and some of its properties. A system of nonlinear algebraic equations can be written as f (x) = 0 (6.12)
where x is an n vector of unknowns and f is an n vector function of x. Given an appropriate starting value x0 , the Newton-Raphson method solves this
40
f (x) f ( x0)
f ( x1) x x2 x1 x0
Figure 6.1. Newton-Raphson method in unidimensional case
f (x)
vector equation by generating the following sequence:
f ( x0)
J(x )x = f (x ) x
1
+1
(6.13)
= x + x
f (x ) where J(x ) = f (x)/x is the Jacobian matrix with elements

Jij = fi xj
x x2 x1 x0
(6.14)
6.2.1
Unidimensional case
To get a better feeling for the method we rst study the one dimensional case, and eq. (6.12) becomes f (x) = 0 (6.15) where x is the unknown and f (x) is a scalar function. Figure 6.1 illustrates a simple case in which there is a single solution to eq. (6.15). Under these circumstances, the following algorithm can be used to nd the solution of eq. (6.15): 1. Set = 0 and choose an appropriate starting value x0 ; 2. Compute f (x );
x x2
x1
x0
41
f (x) f ( x0)
f ( x1) x x2 x1 x0
Figure 6.2. Dishonest Newton-Raphson method in unidimensional case
3. Compare f (x ) with specied tolerance ; if |f (x )| , then x = x is the solution to eq. (6.15); Otherwise, |f (x )| > , go to the next step; 4. Linearize f (x) at the current solution point [x , f (x )], as shown in Figure 6.1. That is, f (x + x ) f (x ) + f (x )x , where f (x ) is calculated at x 5. Solve f (x ) + f (x )x = 0 for x , and update the solution estimate, x+1 = x + x , where x = f (x )/f (x ); 6. Update iteration counter + 1 and go to step 2. The dishonest Newton-Raphson method is illustrated in Figure 6.2. In this case at Step 4 of the algorithm, a constant derivative is assigned and f (x ) = f (x0 ). Although the number of iterations required for convergence usually increases, it is not necessary to recalculate the derivatives for each iteration and hence the computation burden at each iteration is lower. When only limited accuracy is needed, the overall performance of the dishonest version may be better than that of the full Newton-Raphson method.
6.2.2
Quadratic Convergence
Close to the solution point x , the Newton-Raphson method normally presents a property called quadratic convergence. This can be proved for the unidimensional case discussed above if it is assumed that x is a simple (not a multiple) root and that its rst and second derivatives are continuous.
42
Hence, f (x ) = 0, and for any x in a certain neighbourhood of x , f (x) = 0. If denotes the error at the -th iteration, i.e. = x x() the Taylor expansion about x yields f (x ) = f (x() + ) x = f (x() ) + f (x() ) + 1/2f ()2 =0 where x [x() , x ]. Dividing by f (x() ), this expression can be written as f () f (x() ) x + + 1/2 () 2 = 0 (x() ) f f (x ) Since, f (x() ) f (x() ) + = () + x x() = x x(+1) = +1 f (x() ) f (x ) the following relationship between and +1 results: +1 1 f () x = 2 2 f (x ) (6.20) (6.19) (6.18) (6.17) (6.16)
In the vicinity of the root, i.e. as x x , x x , and we thus have |+1 | = 1 |f (x )| 2 2 |f (x )| (6.21)
From eq. (6.21) it is clear that the convergence is quadratic with the assumptions stated above.
6.2.3
Multidimensional Case
Reconsider now the n-dimensional case f (x) = 0 where f (x) = (f1 (x), f2 (x), . . . , fn (x))T and x = (x1 , x2 , . . . , xn )T (6.24) Thus f (x) and x are n-dimensional (column) vectors. The Newton-Raphson method applied to to solve eq. (6.22) follows basically the same steps as those applied to the unidimensional case above, (6.23) (6.22)
43
except that in Step 4, the Jacobian matrix J(x ) is used, and the linearization of f (x) at x is given by the Taylor expansion f (x + x ) f (x ) + J(x )x where the Jacobian matrix has the general form f1 x1 f2 f J= = x x .1 . . fn x1 f1 x2 f2 x2 . . . fn x2 .. . f1 xn f2 xn . . . fn xn (6.25)
(6.26)
The correction vector x is the solution to f (x ) + J(x )x = 0 (6.27)
Note that this is the linearized version of the original problem f (x +x ) = 0. The solution of eq. (6.27) involves thus the solution of a system of linear equations, which usually is done by Gauss elimination (LU factorization). The Newton-Raphson algorithm for the n-dimensional case is thus as follows: 1. Set = 0 and choose an appropriate starting value x0 ; 2. Compute f (x ); 3. Test convergence: If |fi (x )| for i = 1, 2, . . . , n, then x is the solution Otherwise go to 4; 4. Compute the Jacobian matrix J(x ); 5. Update the solution x = J1 (x )f (x ) x+1 = x + x 6. Update iteration counter + 1 and go to step 2.
(6.28)
44
6.3
Newton-Raphson applied to the Power Flow Equations
In this section we will now formulate the Newton-Raphson iteration of the power ow equations. Firstly, the state vector of unknown voltage angles and magnitudes is ordered such that x= U (6.29)
and the nonlinear function f is ordered so that the rst components correspond to active power and the last ones to reactive power: f (x) = with P2 (x) P2 . . . (s) Pm (x) Pm f (x) = Q2 (x) Q(s) 2 . . . (s) Qn (x) Qn
(s)
P(x) Q(x)
P(x) P(s) Q(x) Q(s)
(6.30)
(6.31)
In eq . (6.31) the functions Pk (x) are the active power ows out from bus k (s) given by eq. (4.11) and the Pk are the known active power injections into bus k from generators and loads, and the functions Qk (x) are the reactive (s) power ows out from bus k given by eq. (4.12) and Qk are the known reactive power injections into bus k from generators and loads. The rst m 1 equations are formulated for PU and PQ buses, and the last n 1 equations can only be formulated for PQ buses. If there are NP U PU buses and NP Q PQ buses, m 1 = NP U + NP Q and n 1 = NP Q . The load ow equations can now be written as f (x) = P(x) Q(x) =0 (6.32)
and the functions P(x) and Q(x) are called active and reactive (power) mismatches. The updates to the solutions are determined from the equation J(x ) U + P(x ) Q(x ) =0 (6.33)
6.4. P QU Decoupling The Jacobian matrix J can be written as P P U J= Q Q U which is equal to P(x) J= Q(x)
45
(6.34)
or simply
P(x) U Q(x) U
(6.35)
P P U J= (6.36) Q Q U In eq. (6.34) the matrices P/ and Q/U are always quadratic, and so is of course J.
6.4
P QU Decoupling
The ac power ow problem above involves four variables associated with each network node k: Uk , the voltage magnitude k , the voltage angle Pk , the net active power (generation - load) Qk , the net reactive power (generation - load) For transmission systems, a strong coupling is normally observed between P and , as well as between Q and U . This property will in this section be employed to simplify and speed up the computations. In the next section we will derive a linear approximation called dc power ow (or dc load ow). This linear model relates the active power P to the bus voltage angle . Let us consider a -model of a transmission line, where the series resistance and the shunt admittance both are neglected and put to zero. In this case, the active and reactive power ows are given by the following simplied expressions of eqs. (3.3) and (3.4) Pkm = Uk Um sin km xkm (6.37)
46 Qkm =
2 Uk Uk Um cos km (6.38) xkm where xkm is the series reactance of the line. The sensitivities between power ows Pkm and Qkm and the state variables U and are for this approximation given by
Uk Um cos km Pkm = k xkm
Pkm Um sin km = Uk xkm
(6.39) (6.40)
Uk Um sin km 2Uk Um cos km Qkm Qkm = = k xkm Uk xkm When km = 0, perfect decoupling conditions are observed, i.e. Pkm Uk Um = k xkm Pkm =0 Uk
(6.41)
Qkm 2Uk Um Qkm =0 = (6.42) k Uk xkm As illustrated in Figure 6.3, in the usual range of operations (relatively small voltage angles), a strong coupling between active power and voltage angle as well as between reactive power and voltage magnitudes exists, while a much weaker coupling between reactive power and voltage angle, and between voltage magnitude and active power exists. Notice, however, that for larger angles this is no longer true. In the the neighbourhood of km = 90 , there is strong coupling between P and U as well as between Q and . Example 6.1. A 750 kV transmission line has 0.0175 p.u. series reactance (the series reactance and the shunt admittance are ignored in this example). The terminal bus voltage magnitudes are 0.984 and 0.962 p.u. and the angle dierence is 10 . Calculate the sensitivities of the active and reactive power ows with respect to voltage magnitude and phase angle. Solution The four sensitivities are calculated by using eqs. (6.39) and (6.40): Uk Um cos km 0.984 0.962 cos 10 Pkm = = = 54.1 k xkm 0.0175 Pkm Um sin km 0.962 sin 10 = = = 9.5 Uk xkm 0.0175 Qkm Uk Um sin km 0.984 0.962 sin 10 = = = 9.4 k xkm 0.0175 Qkm 2Uk Um cos km 2 0.984 0.962 cos 10 = = = 58.3 Uk xkm 0.0175
6.4. P QU Decoupling
47
P,Q
VkVm / xkm Operation range Pkm qkm
-p/2
p /2
Qkm
Figure 6.3. P and Q curves for a line with a series resistance and a shunt admittance of zero and considering terminal voltages Uk = Um = 1.0 p.u.
f (x)
As seen the P and Q U couplings are much greater than the other couplings.f ( x0) If the couplings Q and P U are neglected the Newton-Raphson iteration scheme can be simplied. With this assumption the Jacobian Matrix can be written as P 1) f (x 0 JDEC = (6.43) Q x 0 U
2 1 0
x x x Thus there is no coupling between the updates of voltage magnitudes and angles and eq. (6.33) can be written as two uncoupled equations:
P + P( , U ) = 0 f ( x ) (6.44)
f ( x0)
Q U + Q( +1 , U ) = 0 U
(6.45)
f ( x1) x
48
In this formulation two systems of linear equations have to be solved instead of one system. Of course, the total number of equations to be solved is the same, but since the needed number of operations to solve a system of linear equations increases more than linearly with the dimension, it takes less operations to solve eqs. (6.44) and (6.45) than eq. (6.33) with the complete Jacobian matrix, J. It should be noted that if the iterations of eqs. (6.44) and (6.45) converge, it converges to a correct solution of the load ow equations. No approximations have been introduced in the functions P(x) or Q(x), only in the way we calculate the updates. The convergence of the decoupled scheme is somewhat slower than the full scheme, but often the faster solution time for the updates compensates for slower convergence, giving as faster overall solution time. For not too heavily loaded systems a faster overall solution time is almost always obtained. The two equations in eqs. (6.44) and (6.45) are solved sequentially, and then the updated unknowns of the rst equation, eq. (6.44), can be used to calculate the mismatches of the second system of equations, eq. (6.45), resulting in an increased speed of convergence. A number of approximations can be made to calculate the matrix elements of the the two sub-matrices of the Jacobian in eqs. (6.44) and (6.45). This will only inuence the speed of convergence of the solution. If the method converges to a solution, this is the correct solution as long as the accurate expression for P and Q are used. If approximation regarding the active and reactive power mismatches are used, the solution can be even faster, but then only a approximative solution will be obtained. This is further elaborated in next section.
6.5
Approximative Solutions of the Power Flow Problem
In the previous section the exact expressions of the power ow equations were used. However, since the power ow equations are solved frequently in the operation and planning of electric power systems there is a need that the equations can be solved fast, and for this purpose the approximations introduced in this chapter have proved to be of great value. Often the approximations described here are used together with exact methods. Approximative methods could be used to identify the most critical cases, which then are further analysed with the full models. The fast, approximative methods can also used to provide good initial guesses for a complete solution of the equations.
6.5.1
Linearization
In this subsection the linearized dc power ow equations will be derived.
6.5. Approximative Solutions of the Power Flow Problem Transmission Line
49
Consider again expressions for the active power ows (Pkm and Pmk ) in a transmission line:
2 Pkm = Uk gkm Uk Um gkm cos km Uk Um bkm sin km
(6.46) (6.47)
Pmk =
2 Um gkm
Uk Um gkm cos km + Uk Um bkm sin km
These equations can be used to determine the real power losses in a transmission line
2 2 Pkm + Pmk = gkm (Uk + Um 2Uk Um cos km)
(6.48)
If the terms corresponding to the active power losses are ignored in eqs. (6.46) and (6.47), the result is Pkm = Pmk = Uk Um bkm sin km (6.49)
The following additional approximations are often valid, particularly during light load conditions Uk Um 1 p.u. (6.50) sin km km And since bkm = 1/xkm we can simplify the expression for the active power ow Pkm to Pkm = km /xkm = k m xkm (6.53) (6.52) (6.51)
This equation is analogous to Ohms law applied to a resistor carrying a dc current: Pkm is the dc current; k and m are the dc voltages at the resistor terminals; xkm is the resistance. This is illustrated in Figure 6.4. Series Capacitor For a given voltage angle spread, the active power ow in a transmission line decreases with the line reactance (and series reactance normally increases with line length). Series compensation aims at reducing the eective electric length of the line: a series capacitor connected in series with the line. If, for example, a 40% compensation corresponds to a capacitor with a reactance
50
m rkm + jxkm
jbkm
sh
jbkm
sh
(a)
xkm
m (b)
Pkm = (qk - q m ) x km
Figure 6.4. Transmission line. (a) Equivalent -model. (b) DC power ow model
PU PQ of 40% of the original line reactance, but with opposite PU the resulting sign, 1 reactance of the compensated line becomes 60% of the original value. Thus 2 3
xcomp = xkm xsc km (6.54)
with obvious notation. In the dc power ow model the series capacitor can thus be regarded as a negative resistance inserted in series with the equivalent line resistance. In-Phase Transformer
Uq
PQ
The active power ows, Pkm and Pmk , in an in-phase transformer are given by eq. (3.10) Pkm = (akm Uk )2 gkm akm Uk Um gkm cos km akm Uk Um bkm sin km (6.55) Neglecting the terms associated with losses and introducing the same approximations used for transmission lines yields Pkm = km xkm /akm (6.56)
where further approximating akm 1, i.e. the transformer tap ratio is close to the relation between the nominal voltages of the two sides, yields the
6.5. Approximative Solutions of the Power Flow Problem
51
1:akm
jxkm
m (a)
xkm akm Pkm = (qk - q m ) (xkm akm)
m (b)
Figure 6.5. In-phase transformer. (a) Transformer comprising ideal transformer and series reactance. (b) DC power ow model
same expression as for transmission lines Pkm = This is illustrated in Figure 6.5. Phase Shifter Let us consider again the expression for the active power ow Pkm in a phaseshifting transformer of the type represented in Figure 2.6 with akm = 1 (eq. (3.14)):
2 Pkm = Uk gkm Uk Um gkm cos(km + km ) Uk Um bkm sin(km + km )
km xkm
(6.57)
As with transmission lines and in-phase transformers, if the terms associated with active power losses are ignored and Uk = Um = 1 p.u. and bkm = x1 , the result is km Pkm = sin(km + km ) xkm (6.58)
and if (km + km ) << /2, then linear approximation can be used, giving Pkm = (km + km ) xkm (6.59)
52
qk +jkm jjkm
jxkm
m (a)
1:e k
xkm
m (b)
Pkm = (qk + jkm - qm ) /xkm xkm m (c) jkm / xkm

Figure 6.6. Phase-shifting transformer. (a) Phase-shifting transformer model (b) Thvnin dc power ow model. (c) Norton dc power e e ow model
Note that Pkm has two components, the rst depending on the terminal bus voltage angles, km/xkm , and the other depending only on the phaseshifting transformer angle, km /xkm . If km is considered to be a constant, eq. (6.59) can be represented by the linearized model shown in Figure 6.6, where the constant part of the active power ow, km /xkm , appears as an extra load on the terminal bus k and an extra generation on the terminal bus m if km > 0, or vice-verse if km < 0.
6.5.2
Matrix Formulation of DC Power Flow Equations
In this section, the dc model developed above is expressed in the form I = YE. According to the dc model, the active power ow in a branch is given by Pkm = x1 km km (6.60)
where xkm is the series reactance of the branch (parallel equivalent of all the circuits existing in the branch).
6.5. Approximative Solutions of the Power Flow Problem The active power injection at bus k is thus given by Pk =
mk
53
x1 km = ( km
mk
x1 )k + km
mk
(x1 m ) km
(6.61)
for k = 1, 2, . . . , N , where N is the number of buses in the network. This can be put into matrix form as follows: P = B where P is the vector of the net injections Pk B is the nodal admittance matrix with the following elements:
Bkm = x1 km Bkk = mk
(6.62)
x1 km
is the vector of voltage angles k The matrix B in eq. (6.62) is singular, i.e. with a determinant equal to zero. This means that the system of equations in eq. (6.62) has no unique solution. Since the power losses have been ignored, the sums of the components of P is equal to zero. This means that the rows of B are linearly dependent. To make the system solvable, one of the equations in the system is removed, and the bus associated with that row is chosen as the angle reference, i.e. ref = 0. In forming the matrix B , in-phase and phase-shifting transformers are treated like transmission lines. The phase-shifting transformers also contribute to the construction of the independent vector P with the Norton equivalent injections shown in Figure 6.6(c). Example 6.2. Consider the network given in Figure 6.7 in which the reference angle is 1 = 0. Use the dc power ow method to calculate the power ows in the lines. Solution In this case, the elements of the matrix B are calculated as B22 = x1 + x1 = (1/3)1 + (1/2)1 = 5 21 23 B23 = x1 = (1/2)1 = 2 23 B32 = x1 = (1/2)1 = 2 32 B33 = x1 + x1 = (1/2)1 + (1/2)1 = 4 31 32
54
P1 = 1.5 1 P12 x12 = 1/3 P13 x13 = 1/2
P2 = -0.5 2
P23 x23 = 1/2
3 P3 = -1.0
Figure 6.7. 3-bus network. (Active power in p.u.; branch reactances in p.u.)
and thus B = and (B )1 =
5 2 2 4 1/4 1/8 1/8 5/16
The nodal voltage angles (in radians) can now easily be calculated = 1 2 = (B )1 P = 1/4 1/8 1/8 5/16 0.5 1.0 = 0.250 0.375
The power ows in the transmission lines are according to the dc power ow model P12 = x1 12 = 3 0.25 = 0.75 p.u. 12
P13 = x1 13 = 2 0.375 = 0.75 p.u. 13 P23 = x1 23 = 2 0.125 = 0.25 p.u. 23
6.5. Approximative Solutions of the Power Flow Problem
55
6.5.3
DC Power Flow Model
The linearized model P = B can be interpreted as the model for a network of resistors fed by dc current sources where P is the vector of nodal current injections, is the nodal vector of dc voltages, and B is the nodal conductance matrix, as illustrated in Figure 6.8.
56
1 P1 x12
P2
2 x23
3 P3
x15
x25
x34
(a)
P5 x45 5 4
P4
I2 I1 r12 I3
r23
r15
r25
r34
(b)
I4 r45
Figure 6.8. 6-bus network. (a) power network. (b) dc power ow model.
7
Fault Analysis
This chapter presents computation techniques for the currents that may occur in the network under symmetrical short circuit conditions. These short circuit currents determine the rating of circuit breakers that must be able to clear the fault in order to avoid permanent damage to the equipment.
O FAR we have dealt with steady state behavior of power systems under normal operating conditions. This chapter is devoted to abnormal system behavior under conditions of faults. Such conditions are caused in the system accidentally through insulation failure of equipment or ashover of lines initiated by a lightning stroke or through accidental faulty operation. In high voltage networks, short circuits are the most frequent type of faults. Short circuits may be solid or may involve an arc impedance. Figure 7.1 illustrates dierent types of short circuits. Depending on the location, the type, the duration, and the system grounding short circuits may lead to electromagnetic interference with conductors in the vicinity (disturbance of communication lines), stability problems, mechanical and thermal stress (i.e. damage of equipment, personal danger) danger for personnel The system must be protected against ow of heavy short circuit currents by disconnecting the faulty part of the system by means of circuit breakers operated by protective relaying. The safe disconnection can only be guaranteed if the current does not exceed the capability of the circuit breaker. Therefore, the short circuit currents in the network must be computed and compared with the ratings of the circuit breakers at regular intervals as part of the normal operation planning. As illustrated in Figure 7.2, the short circuit currents at network nodes are generally increasing over the years due to more generators, 57
58
a) Symmetrical three-phase short circuit a b c
7. Fault Analysis
b) Two-phase without ground contact a b c
c) Two-phase with ground contact a b c
d) Single-phase earth fault a b c
SC current
partial SC current in conductor or ground
Figure 7.1. Examples for dierent types of short circuits.
new lines in existing networks, interconnection of isolated networks to an integrated one. This is primarily a problem for the expansion planning, where the impacts of long-term changes on the short circuit currents have to be assessed. If the short circuit current exceeds the admissible limit at a network node, the circuit breakers have to be replaced by breakers with higher ratings. Alternatively, the impedance between feeder and fault location can be increased in order to reduce the short circuit current. This is commonly achieved by introducing a higher voltage level while splitting the existing lower voltage network (Figure 7.3), use of multiple busbars (Figure 7.4), fast decoupling of busbars during faults (Figure 7.5). Since changing the circuit breakers involves very high costs, the proposed means to reduce the short circuit currents are the generally preferred solution. However, this results in a more complex network structure. It
59
Figure 7.2. Development of short circuit currents over the years.
Figure 7.3. Introduction of a higher voltage level.
also leads to more possibilities to recongure the network topology during operation. When a multiple busbar is introduced, for example, a line can be switched from one busbar to another. Switching actions have a signicant inuence on the short circuit cur-
60
7. Fault Analysis
Figure 7.4. Multiple busbar operation.
ts = 70 . . . 80ms
Figure 7.5. Fast busbar decoupling.
rents, but not all possible topologies can be studied during network planning. Therefore, calculating the short circuit currents has become more and more a problem for the network operation planning. Prior to each switching action all short circuit currents of the new topology must be calculated in order to decide if the switching action may be carried out. This requires computation algorithms that are sucient fast for real time applications. The majority of system faults are not three-phase faults but faults involving one line to ground or occasionally two lines to ground. These are unsymmetrical faults requiring special computational methods like symmetrical components. Though symmetrical faults are rare, symmetrical short circuit analysis must be carried out, as this type of fault generally leads to the most severe fault current ow against which the system must be protected. Symmetrical fault analysis is, of course, simpler to carry out. A power network comprises synchronous generators, transformers, lines, and loads. Though the operating conditions at the time of fault are important, the loads can usually be neglected during short circuits, as voltages dip very low so that currents drawn by loads can be neglected in comparison with short circuit currents. The synchronous generator during short circuit has a characteristic timevarying behavior. In the event of a short circuit, the ux per pole undergoes dynamic change with associated transients in damper and eld windings. The reactance of the circuit model of the machine changes in the rst few
7.1. Transients on a transmission line
61
cycles from a low subtransient reactance to a higher transient value, nally settling at a still higher synchronous (steady state) value. Depending upon the arc interruption time of the circuit breakers, an appropriate reactance value is used for the circuit model of synchronous generators for the short circuit analysis. In a modern large interconnected power system, heavy currents owing during a short circuit must be interrupted much before the steady state conditions are established. Furthermore, from the considerations of mechanical forces that act on the circuit breaker components, the maximum current that a breaker has to carry momentarily must also be determined. For selecting a circuit breaker we must, therefore, determine the initial current that ows on occurrence of a short circuit and also the current in the transient that ows at the time of circuit interruption. We distinguish between two dierent approaches to calculate the short circuits in a power system: Calculation of transient currents Calculation of stationary currents First, we will focus on the calculation of transient currents since this will help us to understand the physical phenomena during short circuits. However, for large power systems, the computation of transient currents is not feasible. For this reason simplied techniques for short current computation will be presented that are based on stationary models.
7.1
Transients on a transmission line
Let us consider the short circuit transient on a transmission line. Certain simplifying assumptions are made at this stage: 1. The line is fed from a constant voltage source. 2. Short circuit takes place when the line is unloaded. 3. Line capacitance is negligible and the line can be represented by a lumped RL series circuit.
R
v= 2U sin(t + )
Figure 7.6. Transmission line model.
62
7. Fault Analysis
With the above assumptions the line can be represented by the circuit model of Figure 7.6. The short circuit is assumed to take place at t = 0. The parameter controls the instant on the voltage wave when short circuit occurs. It is known from circuit theory that the current after short circuit is composed of two parts, i.e. i = is + it where is represents the steady state alternating current 2U is = sin(t + ) |Z| and it represents the transient direct current 2U (R/L)t it = is (0) e = sin( ) e(R/L)t |Z| with Z= R2 + 2 L2 = tan1 L R . (7.1)
(7.2)
(7.3)
(7.4)
A plot of i = is + it is shown in Figure 7.7. In power system terminology, the sinusoidal steady state current is called the symmetrical short circuit current and the unidirectional transient component is called the DC o-set current, which causes the total short circuit current to be unsymmetrical till the transient decays.
i imm i = is + it
Figure 7.7. Waveform of a short circuit current on a transmission line.
It follows easily from Figure 7.7 that the maximum momentary short circuit current imm corresponds to the rst peak. If the decay of transient current in this short time is neglected, 2U 2U imm = sin( ) + (7.5) |Z| |Z|
7.2. Short circuit of a synchronous machine Since transmission line resistance is small, 90 . 2U 2U imm = cos + |Z| |Z|
63
(7.6)
This has the maximum possible value for = 0, i.e. short circuit occurring when the voltage wave is going through zero. Thus imm may be a high as twice the maximum of the symmetrical short circuit current: 2U imm 2 (7.7) |Z| For the selection of circuit breakers, momentary short circuit current is taken corresponding to its maximum possible value. The next question is what is the current to be interrupted ? As has been pointed out earlier, modern circuit breakers are designed to interrupt the current in the rst few cycles (ve cycles or less). With reference to Figure 7.7 it means that when the current is interrupted, the DC o-set it has not yet died out and contributes thus to the current to be interrupted. Rather than computing the value of the DC o-set at the time of interruption (this would be highly complex in a network of even moderately large size), the symmetrical short circuit current alone is calculated. This current is then increased by an empirical multiplying factor to account for the DC o-set current.
7.2
Short circuit of a synchronous machine
Under steady state short circuit conditions, the armature reaction of a synchronous generator produces a demagnetizing ux. In terms of a circuit this eect is modelled as a reactance Xa in series with the induced emf. This reactance when combined with the leakage reactance Xl of the machine is called synchronous reactance Xd . The index d denotes the direct axis. Since the armature reactance is small, it can be neglected. The steady state short circuit model of a synchronous machine is shown in Figure 7.8.
Xl Xa
Eg
Xd Synchronous reactance
Figure 7.8. Steady state short circuit model of a synchronous machine.
Consider now the sudden short circuit of a synchronous generator that has initially been operating under open circuit conditions. The machine
64
7. Fault Analysis
undergoes a transient in all the three phases nally ending up in the steady state condition discribed above. The circuit breaker must interrupt the current long before the steady condition is reached. Immediately upon short circuit, the DC o-set currents appear in all three phases, each with a different magnitude since the point on the voltage wave at which short circuit occurs is dierent for each phase. These DC o-set currents are accounted for separately on an empirical basis. Therefore, for short circuit studies, we need to concentrate our attention on the symmetrical short circuit current only. In the event of a short circuit, the symmetrical short circuit current is limited initially only by the leakage reactance of the machine. Since the air gap ux cannot change instantaneously, to counter the demagnetization of the armature short circuit current, currents appear in the eld winding as well as in the damper winding in a direction to help the main ux. These currents decay in accordance with the winding time constants. The time constant of the damper winding which has low X/R-ratio is much less than the one of the eld winding, which has high leakage inductance with low resistance. Thus, during the initial part of the short circuit, the damper and eld windings have transformer currents induced in them. In the circuit model their reactancesXf of eld winding and Xdw of damper winding appear in parallel with Xa as shown in Figure 7.9.
Xdw Xl Xf Xa Eg Subtransient reactance
Figure 7.9. Approximate circuit model during subtransient period of short circuit.
Xf Xl Xa Eg Transient reactance
Figure 7.10. Approximate circuit model during transient period of short circuit.
As the damper winding currents are rst to die out, Xdw eectively becomes open circuited and at a later stage Xf becomes open circuited.
7.2. Short circuit of a synchronous machine
65
The machine reactance thus changes from the parallel combination of Xa , Xf , and Xdw during the initial period of the short circuit to Xa and Xf in parallel (Figure 7.10) during the middle period. The machine reactance nally becomes Xa in steady state (Figure 7.8). The reactance presented by the machine in the initial period of the short circuit, i.e.
Xd = Xl +
1 1/Xa + 1/Xf + 1/Xdw
(7.8)
is called the subtransient reactance of the machine; while the reactance effective after the damper winding currents have died out, i.e.
Xd = Xl +
1 1/Xa + 1/Xf
(7.9)
is called the transient reactance. Of course, the reactance under steady conditions is the synchronous reactance. Obviously Xd < Xd < Xd . The machine thus oers a time-varying reactance which changes from Xd to Xd and nally to Xd .
i
2Isc
Subtransient period Transient period Steady state
2Isc
2Isc
0 t
Figure 7.11. Symmetrical short circuit armature current in synchronous machine.
If we examine the oscillogram of the short circuit current of a synchronous machine after the DC o-set has been removed, we will nd the current wave shape as given in Figure 7.11. The short circuit current can be divided into three periodsthe initial subtransient period when the current is large as the machine oers subtransient reactance, the middle transient period where the machine oers transient reactance, and nally the steady state period when the machine oers synchronous reactance. Hence we distinguish between
66 the steady state current Isc ,

the transient current Isc , and the initial subtransient current Isc .
7. Fault Analysis
None of the three currents includes a DC component. The steady state and transient currents are obtained by extrapolating their envelopes backwards in time. Since we are normally interested in maximum short circuit current, most short circuit computation is based only on the subtransient current. In some systems a weighted average of the subtransient and the transient reactance is used.
7.3
7.3.1
Algorithms for short circuit studies

Generator model
In the power ow analysis we are modeling generators as constant power sources (PU characteristics) which leads to nonlinear power ow equations. However, for the short transient period during a short circuit this model is no longer valid since the power and the voltage regulators operate with much large time constants. Here, the linear model outlined in Figure 7.12 is more appropriate. It is used for generators and infeeds, where E represents the subtransient emf and Xd represents the subtransient reactance of generators or the internal grid impedance of a feeder. This model may also be applied to large motor loads.
R E
Xd
Figure 7.12. Linear generator model for short circuit computation.
7.3.2
Simplications
When computing short circuits in an power system further simplications can be made. The following simplications are also proposed in the German standard VDE 0102: All line capacitances are ignored. All non-motor shunt impedances are ignored; motor loads are treated the same way as generators.
7.3. Algorithms for short circuit studies
67
The voltage magnitude and phase angle of generators and infeeds are all set to the same value
Ei = c Unom
(7.10)
where Unom is the nominal voltage of the system in which the short circuit occurs. In high voltage systems ( 35 kV) c = 1.1 represents the dierence between the eective voltage and the system voltage. If the minimum rather than the maximum initial AC short circuit current is to be calculated, c is set to c = 0.95. All tap changing transformers are in middle position. These simplications are indicated for studies regarding medium- and long-term network planning. In the planning stage, the calculations are based on estimated and hence inaccurate data. Therefore, the demands on the short circuit computation algorithm are lower than for real-time applications in the network operation, where accurate results are desired. Studies have shown that the shunt elements and loads have little inuence on the short circuit currents (0.5 % . . . 4 %) and may compensate each other. However, disregarding the actual generator pole voltages and the actual positions of tap changing transformers may sometimes lead to errors of up to 30 %.
7.3.3
Solving the linear system equations
With the linear models of the network elements, a system of linear equations can be set up for a short circuit at any node i. Figure 7.13 illustrates the network model. The voltage sources of the generators (Figure 7.12) must be transformed into equivalent current sources. The admittance of the current source is considered in the respective element of the admittance matrix.
Generator nodes I1 =
E1 Xd
Load nodes
j
Ij = 0
Network Ih =
Eh Xd
h i
Isc
In = 0
Short circuit node
Figure 7.13. Network with short circuit at node i.
68
7. Fault Analysis
In order to calculate the short circuit current Isc at node i the equation
YU=I must be solved, where Y is the nodal admittance matrix U is the nodal voltage vector U= Uj 0 if j = i if j = i
(7.11)
(7.12)
I is the vector of injection currents Ej for generator nodes j = i X dj 0 for load nodes j = i I= Isc for short circuit node without generation j = i I + Ej sc for short circuit node with generation j = i Xdj (7.13) All elements Uj of vector U are unknown variables with the exception of voltage Ui = 0 at the short circuit node i. On the contrary, all elements of vector I are known except for the current Ii at the short circuit node i. To solve the system of equations (7.11), the ith row and the ith column are removed. The remaining admittance matrix is then factorized into triangular matrices. Backward substitution yields the unknown voltages. The unknown current Ii at the short circuit node i can be calculated from the voltages using the ith row of eq. (7.11). If the short circuit currents Isci are calculated for all nodes i = 1, . . . , n the equation system has to be set up and solved for every short circuit node. Since always another row and another column are removed the matrix must be factorized each time anew. Therefore, with the proposed method the computation of all short circuit currents is considerably expensive.
7.3.4
The superposition technique
If two equal voltage sources in reverse direction are introduced into the short circuit path at node i according to Figure 7.14, neither currents nor voltages in the network are aected. Since the network consists only of linear elements the short circuit calculation is a linear problem. Therefore, the superposition principle can be applied, that means the computation of the short circuit current may be performed in two steps. These steps
7.3. Algorithms for short circuit studies
69
are illustrated in Figure 7.15 and Figure 7.16. The superposition, i.e. the addition of the results nally yields the desired values. The advantage of inserting two reverse voltage sources becomes clear if both steps are regarded separately.
Generator nodes I1
1
Load nodes
j
Network Ih
h i n
Short circuit node
Uli Ui = 0 Uli
Isc
Figure 7.14. Network with short circuit at node i.
Generator nodes I1
1
Load nodes
j
Generator nodes
1
Load nodes
j
Network Ih
h i n h
Network
n i
Short circuit node
Short circuit node
Uli I=0
Isc
Uli
Figure 7.15. Superposition 1st step.
Figure 7.16. Superposition 2nd step.
First step Principally, the voltage Uli can be chosen arbitrarily. If it is chosen in such (1) a way that Isc = 0, the rst step corresponds to the computation of node voltages and branch currents in a power ow analysis. The power ow analysis yields the voltage Uli , i.e. the voltage at the short circuit node i
70
7. Fault Analysis
before the fault occurs. If the above simplications are applied, the rst step can be omitted. For all generators, infeeds and loads the voltage is then uniformely presumed
Ei = c Unom .
(7.14)
As already mentioned, disregarding the actual generator voltages may lead to considerable errors. Therefore, this approximation is not always recommended. The rst step, i.e. calculating the pre-fault voltages, is identical for all short circuit nodes i = 1, . . . , n and needs to be done only once. Thus, the additional eort for an exact calculation of the node voltages is comparatively small. Second step In the second step the voltage Uli is inserted in reverse direction at the short circuit node i. All other sources like generators and infeeds are disregarded. This yields the equation system Y U(2) = I(2) with Elements of vector U(2) : U(2) = Elements of vector I(2) : I(2) = 0 Isc
(2)
(7.15)
Uj Uli
(2)
if j = i if j = i
(7.16)
if j = i if j = i
(7.17)
Vector U(2) contains the node voltages Uj , all of which are unknown with the exception of Ui = Uli . In vector I(2) the only non-zero current is Ii . Dividing the system of equations (7.15) by Isc yields the modied system (2) contains only unknown and (2) only given of equations (7.18) where U I values. Y U(2) = (2) I (7.18) with Elements of vector U(2) : (2)
(2) (2)
(2) Uj I
sc
if j = i (7.19) if j = i
U li Isc
7.3. Algorithms for short circuit studies Elements of vector (2) : I (2) = I 0 1 if j = i if j = i
71
(7.20)
The transformation from eq. (7.15) to eq. (7.18) leaves the admittance matrix Y unchanged. As a consequence, the triangular factorization of the admittance matrix must be done only once, even if multiple short circuit cases are investigated. The superposition technique should therefore be preferred to the direct solution of the linear system of equations (7.11). However, if we compute the solution of eq. (7.18) with the known technique, i.e. triangular factorization and forward-backward-substitution, a new forward-backward-substitution is required for every short circuit case. This disadvantage can be avoided with the Takahashi method.
7.3.5
The Takahashi method
The principle idea behind the Takahashi method is the solution of eq. (7.18) by means of inversion. This yields Z (2) = U(2) I where Z = Y 1 (7.22) is the impedance matrix of the system. The unknown initial short circuit current Isc at node i can be found in the ith element of vector U(2) Uli (2) Ui = Isc resp.
Isc =
(7.21)
(7.23)
Uli . (2) U
i
(7.24)
From eqs. (7.20) and (7.21) follows (2) Ui = zii (7.25)
where zii is an element on the main diagonal of matrix Z. Inserting eq. (7.25) in eq. (7.16) nally yields: Ui Isc = (7.26) zii For the calculation of the short circuit current only the main diagonal element zii of the impedance matrix Z is required. Knowing all main diagonal elements zii (i = 1, . . . , n) is therefore sucient to compute the short circuit currents Isci at all nodes i.
72
7. Fault Analysis
Often, also the partial short circuit currents Iji are to be determined. These are the currents that ow from adjacent nodes j to the short circuit node i. With the simplication that all generator voltages are identical the partial short currents are Iji = Uj = Uj
(2) (2)
Ui
(2) (2)
yji yji (7.27)
Uli
where yji is the admittance of the branch between j to i. If this simplication is not permitted the terms from the rst step must be added here. With eq. (7.21) follows:
Iji = zji Isci + zii Isci yji = Isci (zii zji ) yji
(7.28)
For the calculation of the partial short circuit currents from node j to the short circuit node i in addition to the main diagonal element zii all elements zji are required. Because node j is connected directly to node i, the element yji in the admittance matrix is non-zero. Although the impedance matrix Z is not sparse like the admittance matrix Y relatively few elements of Z must be known to calculate the short circuit currents and the partial short circuit currents for all nodes. Apart from the main diagonal elements zii these are the elements zji that are non-zero in the corresponding admittance matrix. In the Takahashi method only those elements of Z are determined that have non-zero counterparts in the triangularly factorized admittance matrix. Except for the few ll-ins that emerge during factorization these are the elements of Z essential for the short circuit computation. The method shall be demonstrated here: Eq. (7.22) can be written as YZ=E (7.29)
where E is the identity matrix. Y is factorized into triangular matrices and eq. (7.29) can be transformed in the following way: LDRZ=E DRZ=L | L1 (7.30)
1 1
RZ=D
Z = D1 L1 + (E R) Z
Z = D1 L1 + Z R Z
|D
|+ZRZ (7.31)
The inverse of diagonal matrix D is again a diagonal matrix. Each element on the main diagonal contains the reciprocal value of the element before the
7.3. Algorithms for short circuit studies inversion: d11 d22 .. 0 . dnn 0 1
1 d11
73
0
1 d22
D1 =
..
.
1 dnn
(7.32)
The inverse of the lower triangular matrix L is a triangular matrix with the same structure. Since all elements on the main diagonal of L are 1, the determinant is also 1. Thus, the main diagonal remains unchanged after inversion. The values of the remaining elements are of no interest due to the nesting of the system of equations (7.31). 1 = 1 0 1 1 = .. . 1 0 1 .. . 1
L1
(7.33)
The product of the diagonal matrix D1 with the triangular matrix L1 is a triangular matrix with modied main diagonal. Therefore, eq. (7.31) has the following form: z11 z12 z21 z22 . . . . . . zn1 zn2 1 z1n d11 1 z2n d22 . . = . . . . znn 0 r12 0 0 . . . . . . 0 0 0 0 0 .. . .. .
1 dnn
z11 z12 z21 z22 . . . . r(n1)n . . 0 zn1 zn2 r1n r2n . . .
(7.34)
z1n z2n . . . . . . znn
With the system of equations (7.34) all elements of the impedance matrix Z can be determined recursively. We start with the main diagonal element znn on the nth row and continue with the row (n 1). For each row, rst the secondary diagonal elements zji and then the main diagonal element zii are calculated. Here we take advantage of the fact that the admittance matrix Y and also their inversethe impedance matrix Zare symmetrical. This is due to the simplifying assumption that all tap changing transformers are in middle
74 position.
7. Fault Analysis
zij
i<j
= zji = 1 dii
not calculated if rij = 0

n
zii =
ril zlj
(7.35)
l=i+1
ril zlj
(7.36)
l=i+1
u d d d d
d d du d
3
Figure 7.17. Sample network.
Example 7.1. The whole procedure shall be demonstrated in single steps on the sample network given in Figure 7.17. With the arbitrarily chosen node numbering the admittance matrix Y has the following structure.
y11 y12 y21 y22 Y = 0 y32 0 0 y51 0
0 0 y23 0 y33 y34 y43 y44 y53 y54
y15 0 y35 y45 y55
(7.37)
Solution During the triangular factorization of Y a ll-in element c25 is generated in the upper triangular matrix. The recursive solution of the
7.3. Algorithms for short circuit studies system of equations (7.34) is then z55 = z45 z44 z35 1 d55 = z54 = r45 z55 1 = r45 z54 d44 = r34 z45 r35 z55
75
(7.38) (7.39) (7.40) (7.41) (7.42) (7.43) (7.44) (7.45) (7.46) (7.47) (7.48) (7.49)
z34 = r34 z44 r35 z54 1 z33 = r34 z43 r35 z53 d33 z25 = r23 z35 r25 z55 z23 = r23 z33 r25 z53 1 z22 = r23 z32 r25 z52 d22 z15 = r12 z25 r15 z55
z12 = r12 z22 r15 z52 1 z11 = r12 z21 r15 z51 . d11 Figure 7.18 illustrates this sequence.
u' uC e u I ' u u uC e B u' u' u e u' u T u u Calculated Element e Element not required
Figure 7.18. Solution sequence with Takahashi method.
Until now we have assumed that the admittance matrix Y is symmetric. This assumption is valid only if all tap changing transformers are in middle position. However, the Takahashi method can also be applied to unsymmetric matrices. In this more general case we get instead of eq. (7.31) two similar matrix equations. With one of these equations the secondary diagonal elements above the main diagonal are calculated, with the other
76
7. Fault Analysis
equation the elements below the main diagonal. The elements on the main diagonal can be calculated with either equation. Unsymmetric matrices require about double the eort of symmetric matrices. Since neglecting the actual tap positions can lead to errors of 30 %, considering the actual tap positions is indicated for network operation applications, that require short circuit computations with high accuracy. Concluding remarks The Takahashi method is an algorithm especially designed for the demands of short circuit computations. Having the same accuracy it needs only a fraction of the computation time of conventional methods. This advantage is important in real-time applications. However, with the Takahashi method only the initial short circuit currents can be calculated. This includes the partial short circuit currents Iji that ow from the adjacent nodes to the short circuit node i. It is not possible to determine the short current shares delivered from individual generators or infeeds. Thus the breaking current at the generators can not be determined with the Takahashi method but requires Gauss elimination.
Part II
Power System Dynamics and Stability
77
8
Classication and Denitions of Power System Stability
This chapter contains a brief overview of some of the basic concepts in power system dynamics. A classication of dierent dynamic phenomena is given and the requirements on modelling are discussed. Stability is a fundamental concept and a denition of stability in power systems is introduced. Instabilities in power systems can be classied according to their physical origins, and this classication is reviewed. Finally, a list of further reading in the subject is given.
N A ELECTRIC POWER SYSTEM a great variety of dierent dynamics occurs. These dynamic phenomena have dierent physical origin and they occur in dierent time scales. A system is in a dynamic states if the time derivative of any system quantity is non zero, and to describe a dynamic system mathematically dierential or dierence equations must be employed. As an example, the electric consumer loads in a power system vary spontaneously all the time, and consequently one could state that the power system is never in steady state in a strict mathematical sense. Such a point of view is however unpractical in most cases, and one must, as when it comes to most large and complex systems, study and analyse more or less simplied models of the system. To do relevant and adequate simplications is benecial for the analysis, and contributes also to achieving results that are easy to understand and easy to interpret. In addition, when deriving simplied models one is forced to identify the most important processes and phenomena in the studied systems, which provides insight and deeper understanding. Thus, as an example, the assumption that the loads are constant and that the electric power output from the generators are constant are relevant approximations for many studies of a power system, and the power system is under these assumptions described by a solution to the power ow equations. In the following we will refer to this state as the steady state of the power system. 1
In order to emphasise that this state is an approximation or idealisation, it is sometimes referred to as quasi steady state. More generally, a quasi steady state is a state that strictly is a dynamic state, but could with sucient accuracy be described by algebraic
1
79
80
8. Classication and Denitions of Power System Stability
8.1
Dynamics in Power Systems
There are a lot of dierent dynamical phenomena with dierent characteristics in a power system. The phenomena could be local, in which case they only involve a minor part of the system or a single component. But they can also involve interactions between dierent parts of the system that might be geographically far from each other. In many cases these system wide interactions are initiated by a local disturbance causing e.g. an earth fault with subsequent change in network topology. In this compendium interactions and phenomena that involve many power system components, e.g. generators and loads, or parts of the system are dealt with. These interactions have in common that they can cause system instabilities that can lead to black outs in large parts of the system, i.e. to interruptions of power supply for many consumers. Dynamics can also be initiated from actions from dierent controllers or from switchings of lines or other components by system operators. Such disturbances should be regarded as normal and should consequently not endanger the stability of the system.
8.1.1
Classication of Dynamics
Dynamic phenomena in power systems are usually classied as 1. Fast (electro-magnetic) transients (100 Hz - MHz) 2. Electro-mechanical swings (rotor swings in synchronous machines) (0,1 3 Hz) 3. Non-electric dynamics, e.g. mechanical phenomena and thermodynamics (up to tens of Hz) Approximate values of typical frequencies are given in brackets. One single initial event in the power system can give rise to dynamics in all the three groups above. A lightning stroke in a power line can induce so high over-voltages that the insulation fails, causing a earth fault. The earth fault can cause rotor swings in synchronous machines with high amplitudes. This can trigger protections to disconnect generators, so that an unbalance between produced and consumed power in the system arises. The frequency in the system drops and generators participating in the frequency control compensate this by increasing their power outputs. Thus the initial lightning stroke has initiated dynamics in all the three groups above. Another way to classify dynamic phenomena is given in Figure 8.1. This classication is based on the time scale of the phenomenon and the (mathematical) models used in analysis.
equations. One such example is fault current calculations, which often is done by using algebraic models even if it is a true dynamic phenomenon, see Figure 8.1.
8.1. Dynamics in Power Systems

Power Flow Harmonic Distortion C. Long Term Stability B.
(Differential and Algebraic Equations) (Algebraic equations)
81
D.
Fault Currents
Short Term Stability Stator Transients Saturation and Resonances
A.
Switching Transients Lightning Transients 10-6 10-3
(Differential Equations)
(Travelling Wave Theory)
103
Time (s)
Figure 8.1. Dynamic phenomena in a power system. In the gure approximate time scales are given and types of mathematical models used. The dierent groups are called: A. Fast (electro-magnetic) transients. B. Synchronous machine dynamics. C. Quasi steady state phenomena D. Steady state phenomena.
8.1.2
Modelling
It is of course almost impossible to develop models that can describe all dynamics in a power systems and still being of practical use. Often one has to utilise a model that captures correctly the specic dynamic phenomenon or interaction that is the aim of the particular investigation. Depending on the purpose of the study the appropriate model of a given power system component could vary signicantly. It is obvious that if the aim is to study relatively slow power oscillations between generators in the system, completely dierent models are required as compared with if one wants to analyse the inuence of lightning impulses in the windings of the synchronous machine. Even if it were theoretically possible to develop a complete model of all the dynamics in the power system, it is questionable if such a model were particularly useful. Firstly, such a model would require an enormous amount if parameter data to be uniquely specied. Secondly, the results obtained from such a model would be very hard to analyse and interpret. Critical review and understanding of obtained results is a necessary prerequisite
82
for sound engineering. When making simulations and computations, which all are done with the help of computers nowadays in system sciences, it is important to have an expectation of what are reasonable outputs. Thus trivial errors due to wrong input data les, mistakes in modelling, etc. can be eliminated to a large extent. The human factor is of utmost importance in computer based analysis and simulation. Models can in principle be erroneous in two dierent ways. Firstly, it can have the wrong structure. It can be too simple overlooking important interactions and processes or modelling them incorrectly. This is of course very serious and might give rise to detrimental consequences. But it is also very serious if wrong parameter data is used in a models of the correct structure. This latter shortcoming occurs not seldom in technical systems, which might look surprising at rst sight. Since technical systems are man made, one should in principle have access to all design parameters dening the system. But it turns out that many parameters, e.g. the gain in the controller, could easily be changed after the system has been commissioned and such changes are not always reported to system analysts. It is obvious that the consequences could be very serious. In technical systems there are of course parameters that are genuinely unknown, e.g. the ground resistivity under a power line. In a large system like a power system, thousands of parameters are needed to dene the system completely. It is a very dicult, but also very important, task to maintain and keep the data bases where all these parameter values are stored updated. This is now a special activity usually referred to as data engineering. In this compendium models needed for the problems to be analysed are developed. Due to space limitations more detailed and elaborate derivation cannot be presented, but the reader is advised to consult other sources, e.g. books in electrical machines or the books listed at the end of this chapter.
8.2
8.2.1
Power System Stability

Denition of Stability
A dynamic phenomenon in a power system is, as said above, initiated by a disturbance in the system. Such a disturbance could as an example be that a line impedance is changed due to an external cause. The behaviour of the system after this disturbance depends of course on a how large this disturbance is. A small disturbance results usually in small transients in the system that are quickly damped out, while a larger disturbance will excite larger oscillations. We all have a intuitive feeling for what is meant with stability in this context. Stability is associated with that the system oscillations decay and that the operation of the power systems can continue without any major impacts for any of the consumers. But, and this is very important, as the power system is a nonlinear system (this will be
8.2. Power System Stability
83
elaborated on later) system stability depends on the kind and magnitude of the disturbance. 2 This distinguishes nonlinear systems from linear systems that can be classied as stable or unstable independent of the disturbance, i.e. stability is a property of the linear system as such. As will be shown later, stability of a power system is strongly coupled to both the magnitude and character of the disturbance as well as to the initial operating point. Over the years many dierent denitions of power system stability have been proposed. The most recent one, which will adopted in these lecture notes, was the result of a joint IEEE/CIGRE working group. In the report in reference 7 of section 8.3, the following denition is given: Denition 8.1. Power system stability is the ability of an electric power system, for a given initial operating condition, to regain a state of operating equilibrium after being subjected to a physical disturbance, with most system variables bounded so that practically the entire system remains intact. The following two comments elaborate on some important aspects of this denition. Comment 8.1. It is not necessary that the system regains the same steady state operating equilibrium as prior to the disturbance. This would be the case when e.g. the disturbance has caused any power system component (line, generator, etc.) to trip. Voltages and power ows will not be the same after the disturbance in such a case. Most disturbances that are considered in stability analyses incur a change in system topology or structure. Comment 8.2. It is important that the nal steady state operating equilibrium after the fault is steady state acceptable. Otherwise protections or control actions could introduce new disturbances that might inuence the stability of the system. Acceptable operating conditions must be clearly dened for the power system under study. As mentioned above there are also other denitions of power system stability in the literature. They are not all identical, but could dier in some details, but most of them are in line with the denition proposed above. Power system dynamics can be modelled by systems of dierential and algebraic equations, and the mathematics of those are studied in the theory of dynamic systems. It is therefore desirable that the stability denitions introduced in a more mathematical stringent way are compatible with the more practically oriented denition above. In the report where the above denition is proposed and motivated, these aspects are further discussed.
To give an example of what is meant with magnitude in this context one could consider a three phase to earth fault on a power line that is cleared by disconnecting the line. The magnitude of this disturbance increases with the fault clearing time. If the fault clearing time is suciently small, the system will remain stable, while a longer fault clearing might cause system instability. Similarly, a high impedance fault gives a smaller disturbance than a solid fault.
2
84
8.2.2
Classication of Power System Stability
To achieve a better overview and structure of stability analyses of power systems, it is of great help to classify possible power system stability. The classication to be introduced here is based on the physical mechanism being the main driving force in the development of the associated instability. It could be either the active or the reactive power that is the important quantity. A common characteristic of the instabilities to be discussed here is that they have their origin in too large an imbalance of active or reactive power in the system, locally or globally. This imbalance can then develop in dierent ways and cause unstable behaviour depending on system characteristics. Exactly what is meant with imbalance will be elaborated in the following.
Rotor Angular or Synchronous Stability The total active electrical power fed into the power system by the generators is always equal to the active power consumed by the loads including the losses in the system. On the other there is not always a similar balance between the loads and the power fed into the generators by the prime movers, e.g. the hydro and steam turbines. If such an imbalance develops the rotating parts of the generators, and other rotating machines, will act as energy buer, and the kinetic energy stored in these will decrease or increase as a result of the imbalance. Rotor angle stability refers to the ability of synchronous machines of a power system to remain in synchronism after a disturbance. If the disturbance is local and substantial, e.g. an earth fault close to an generator, the generator can fall out of step since it has been accelerated during the fault. As quite big currents will ow in the generator windings in such a case, it must be disconnected to avoid that it is damaged. Typical time scale for such an instability to develop is a second to a couple of seconds. This kind of instability is called transient instability, and, instability is usually in form of aperiodic angular separation due to lack of synchronizing torque. This form of instability is also referred to as large-disturbance rotor angle instability. Small-disturbance (or small-signal) rotor angle stability is concerned with the ability of the power system to maintain synchronism under small disturbances. These disturbances are considered to be suciently small that linearization of the system equations is permissible for purposes of analysis. Usually small-disturbance rotor angle stability is associated with insucient damping of oscillations.
8.2. Power System Stability Frequency Stability
85
A third variety of active power imbalance, which is dierent from the ones above, is when the imbalance is not local but global. In the cases of the proceeding cases the sum of active power infeed was enough but there was an imbalance locally. But if the total power fed into the system by the prime movers is less than what is consumed by the loads, including losses, this imbalance will inuence the frequency of the whole system. As explained above the kinetic energy stored in rotating parts of the synchronous machines, and other rotating electrical machines, will compensate for the imbalance resulting in a frequency deviation. If the imbalance is not too large the generators participating in the frequency control will regulate the active power input from their prime movers, and bring back the frequency deviation to acceptable values. If the imbalance is too large, the frequency deviation will be signicant with possible serious consequences. Particularly thermal power plants are sensitive to large frequency drops of long durations, since detrimental oscillations could be excited in the turbines. As a last resort the generators are disconnected, making the situation even more serious. This type of instability is called frequency instability and the time scale could be from a few seconds up to several minutes. Since the involved mechanisms could be quite dierent, one often distinguishes between short-term and long-term frequency instability. In the latter, the control and protections characteristics of turbines, boilers, and reactors play important roles. Voltage Stability When it comes to reactive power balance the situation is not as clear and simple as concerning active power. There is always a balance between produced and consumed reactive power in every node of a network. This is in fact a direct consequence of Kirchos rst current law. When one talks about imbalance in this context we mean that the injected reactive power is such, normally too small, that the voltage in the node cannot be kept to acceptable values. (At low load the injected reactive power could be high resulting in a too high voltage, possibly higher than the equipment might be designed for. This is of course not desirable but it could usually be controlled in such a way that no instabilities develop.) When we talk about imbalance in this case we thus mean that the injected reactive power diers from the desired injected reactive power, needed to keep the desired voltage. If this imbalance gets too high, the voltages become outside the acceptable range. Reactive power is a more local quantity than active power since it cannot be transported as easily in power system where normally X >> R. This explains why voltage problems often are local, and often only occur in part of the system. When the imbalances (voltage problems) develop into insta-
86

Power System Stability
Rotor Angle Stability
Frequency Stability
Voltage Stability
Small-Disturbance Angle Stability
Transient Stability
LargeDisturbance Voltage Stability
SmallDisturbance Voltage Stability
Short Term Short Term Long Term
Short Term
Long Term
Figure 8.2. Classication of power system stability.
bilities these are called voltage instabilities or voltage collapses. In the latter case the instability develops into very low voltages in the system. In principle too high voltages can also occur at a voltage instability. Low voltages arise at high load conditions, while high voltages are associated with low load conditions. Depending on the time scale the voltage instabilities are classied as short-term, a couple of seconds, or long-term, tens of seconds to minutes. The short-term voltage instability involves dynamics of fast acting components such as induction motors, electronically controlled loads, and HVDC converters, while the long-term voltage instability involves slower acting equipment such as tap-changing transformers, thermostatically controlled loads, and generator current limiters. As for rotor angle stability one distinguishes between large-disturbance and small-disturbance voltage stability. The classication of power system instabilities is summarised as in Figure 8.2.
Connection between Instabilities and System Components As explained above the generators, i.e. the synchronous machines, are very important in angular instabilities, and it is sometimes said that these are the driving force in this instability. A more detailed analysis shows that the loads are very often the driving force when it comes to voltage instability, which consequently sometimes is called load instability.
8.3. Literature on Power System Dynamics and Stability Instabilities in Real Systems
87
The classication above is based in simplied an ideal conditions in the system. In a real system these assumptions might not be valid. In real system instabilities it is not seldom a combination of active and reactive power imbalances that trigger an instability. However, in many case it is possible to identify which is the dominating processes in the beginning of the instability. During the course of the dynamics new consequential imbalances might occur, resulting in a combined angular and voltage instability in the nal phase. There are examples of black outs in power systems that have started as slow voltage instabilities, which through low voltages have reduced the power transfer capability resulting in rotor angular instabilities causing the nal collapse of the system. On the other hand, rotor angular instabilities can cause generators to trip, which most systems are designed to cope with, but it can eect the reactive power balance in such a way that voltage instabilities can develop. The purpose of a classication is to dene a structure for a complicated problem and thereby better understanding it. Furthermore it often helps to identify important and critical quantities, processes and components in the system. Classications of this kind should not be driven too far. Most important is always that useful and adequate results are obtained from realistic models of the power system.
8.3
Literature on Power System Dynamics and Stability
These lecture notes should be seen as an introduction to power system dynamics and stability. For those that want to get deeper knowledge in the subject there are a number of books than can be recommended. They have all their strong and weaker sides, and one has often to consult several books to get a complete view of a problem. Below are only references listed that are fairly modern and focus on power system dynamics and stability. 1. Power System Stability and Control by Prabha Kundur. (McGrawHill Inc., 1994, ISBN 0-07-035958-X, 1176 pages) This is the most complete modern book on the subject and is already a classic textbook. It covers most subtopics and its approach is rather practical, but it contains a fair amount of theory also. The book contains a lot of references to other books and published papers. 2. Power System Dynamics and Stability by Jan Machowski, Janusz W. Bialek and James R. Bumby. (John Wiley & Sons Ltd, 1997, ISBN 0-471-97174-X, 461 pages) This book does not contain as many applications as the previous one.
88
8. Classication and Denitions of Power System Stability On the other hand it contains rather detailed motivations and derivations of many of the basic assumptions in power system analysis. A large part of the book is devoted to control and supervision of power systems. 3. Power Systems Dynamics. Stability and Control by K. R. Padiyar. (John Wiley & Sons Ltd, 1996, ISBN 0-471-19002-0, 629 pages) The book gives a good overview with many solved problems. The focus is on angular stability, while voltage stability is on briey dealt with. 4. Power System Dynamics and Stability by Peter W. Sauer and M. A. Pai. (Prentice Hall, 1998, ISBN 0-13-678830-0, 357 pages) The mathematical level in this book is higher than the previous ones and to fully appreciate it knowledge in the theory of nonlinear systems is needed. Only angular stability is included. 5. Power System Voltage Stability by Carson W. Taylor. (McGraw-Hill Inc., 1994, ISBN 0-07-063184-0, 273 pages) This book deals only with voltage stability. The approach is rather practical. Many examples from voltage instabilities in real systems are reviewed and analysed. 6. Voltage Stability of Electric Power Systems by Thierry Van Cutsem and Costas Vournas. (Kluwer Academic Publishers, 1998, ISBN 07923-8139-4, 378 pages) Voltage stability is the subject of this book also. However the approach is more mathematical than in the previous one, and to fully appreciate it, knowledge about nonlinear systems is required. 7. Denition and Classication of Power System Stability IEEE/CIGRE Joint Task Force on Stability Terms and Denitions, 2002.
9
Synchronous Machine Models
This chapter starts with a short discussion that motivates why synchronous machines are important when electro-mechanical oscillations are studied. A brief overview of the operating principles and design of the synchronous machine is given. This forms the basis of simplied models that are derived rst for stationary, then for dynamic operation. Even if the dynamic model is very simple and neglects a number of features, it captures the most important properties to model electro-mechanical oscillations in power systems.
YNCHRONOUS MACHINES , i.e. practically all generators together with synchronous motors and synchronous compensators, are the most important power system components in the analysis of electro-mechanical oscillations in power systems. The oscillations are manifested in that the rotors of the synchronous machines do not rotate with constant angular velocity corresponding to system frequency, but superimposed are low frequency oscillations, typically 0.1 4 Hz. It is important that these superimposed oscillations are not too large, because then the stability of the power system can be endangered. A correct description of these oscillations requires often detailed models of many dierent system components, but to get an understanding of and insight into the basic physical phenomena and processes that determine the stability it is often sucient to employ the simple model that will derived and motivated in this chapter. This simplied model of the synchronous machine together with a simple model of the power transmission system provide a description that will be used in chapter 11. As the name electro-mechanical oscillations suggests, both electrical and mechanical phenomena are involved, i.e. both currents in and voltages across dierent windings in the machines but also the mechanical motion of the rotor. Therefore, models of both electrical and mechanical parts of the synchronous machine are needed.
9.1
Design and Operating Principle
In the synchronous machine a magnetized rotor creates a rotating magnetic eld in the air gap. If the rotor eld is ideally sinusoidal and if the rotor rotates at constant speed, this will induce ideally sinusoidal voltages in the 89
90
R
9. Synchronous Machine Models

R T T F S S S T
F D
T S T R
R
Round rotor p = 2, n = 3000 min1 for f = 50 Hz
Salient pole rotor p = 4, n = 1500 min1 for f = 50 Hz
Figure 9.1. Cross-sections through dierent rotor types.
stator windings. If the machine terminals are connected, the currents owing in the stator windings create a second rotating magnetic eld which causes a torque on the rotor. In a synchronous motor, this torque drives a mechanical load; in a synchronous generator, the magnetic torque opposes the mechanical driving torque of the prime mover (e.g. a turbine). Under balanced, steady-state conditions the magnetic torque is equal to the mechanical torque, and so the rotor continues to rotate at constant speed.
9.1.1
Rotor Types
The magnetic rotor eld is generated by a eld winding F on the rotor which is fed with an adjustable direct current. In addition, the rotor has a shortcircuited damper winding D at the surface. This winding serves to dampen electrical and mechanical oscillations and to shield the eld winding from inverse rotating elds in case of asymmetries or harmonics in the stator currents. (In rotors without an explicitly realized damper winding eddy currents in the rotor iron can have a similar eect.) Depending on the application of the generator, two dierent types of rotors are used that are shown in Figure 9.1. Round Rotor Round rotors are used with high-speed turbines such as steam or gas turbines. For this reason, generators with round rotors are also called turbo generators. They can have ratings as high as 1800 MVA per unit. Due to the large centrifugal forces, the rotor consists of a long, narrow, solid steel cylinder. The eld windings are mounted in slots that are mill-cut into about 2/3 of the perimeter. Because of the discrete distribution of the windings on the
9.1. Design and Operating Principle
91
Figure 9.2. Simplied arrangement of a stator with a solid, unwound rotor.
rotor surface, the magnetic ux density in the air gap always has a stair-step form. Through proper distribution of the windings these stair-steps can be made approximately sinusoidal!
Salient Pole Rotor Salient pole rotors are used with low-speed hydro turbines with rated powers of up to 800 MVA per unit. In order to obtain the appropriate electrical power frequency in spite of the low rotor speed, salient pole rotors typically have multiple pole pairs. For run-of-river power stations the number of poles can be as high as p = 200! Such rotors have very large diameters (several meters) and short lengths. The eld windings are mounted on the individual poles. By properly designing the geometric form of the poles, the magnetic ux density in the air gap at the stator surface can be made ideally sinusoidal!
9.1.2
Stator Field
The magnetic stator eld is generated by an alternating three-phase current owing in the stator windings. The variation of this eld in space and time can best be studied using the simplied arrangement shown in Figure 9.2. In this arrangement, the stator has only one pole pair and the slots are assumed to be innitely narrow. The rotor is replaced by a solid, unwound cylinder. The stator eld is derived here for a stator with one single pole pair. The transfer to the case of a stator with p poles is left as an exercise to the reader. The magnetic eld strength in the iron is neglected and the induction in the air gap is assumed to be radial. Due to the -symmetry of the set-up
92
Amp`res law gives for the dotted path shown in Figure 9.2 e () =
A()
i da =
1 0
S()
B ds =
B() 2 0
(9.1)
where A() is the area enclosed by the path, is the width of the air gap, B() is the ux density in the air gap at the angle , and () is the total ux linkage (or number of ampere windings) enclosed by the path. One single-phase winding If only phase R is fed by an alternating current, the total enclosed ux linkage is either +w iR (t) or w iR (t), depending on whether the path encloses the upper slot (0 < < ) or the lower one ( < < 2). w is the number of windings in one slot. The left side of Figure 9.3 shows the ux linkage of a single-phase winding plotted against the angle for dierent times. Since B() is directly proportional to () as shown in eq. (9.1), the magnetic eld in the air gap is a rectangular, standing wave with a variable amplitude. Three-phase windings If the three phases R, S and T are fed by an alternating three-phase current, the integration path in Figure 9.2 always encloses three slots that carry windings with dierent momentary currents. The resulting total ux as a function of the angle is thus a superposition of three rectangular curves that are shifted by 2/3 and have dierent amplitudes. This is shown on the right side of Figure 9.3. It is evident that the fundamental component of the total ux is a traveling wave with a constant amplitude which moves with the angular velocity . This angular velocity of the stator eld is the same as the angular frequency of the alternating currents. In a real stator, the windings are of course not concentrated in innitely narrow slots. Instead, they are distributed among several slots of nite extent possibly with a dierent number of windings in dierent slots. Through proper design the magnetic eld in the air gap can be made approximately sinusoidal. It should be noted that this simplied treatment of the stator eld is only applicable for generators with round rotors where the air gap has a constant width. In the case of a salient pole rotor, the width of the air gap is variable, and therefore the magnetic stator eld depends not only on the momentary winding currents but also on the momentary rotational angle of the rotor. In some cases, the dierences between round rotor and salient pole rotor can be neglected (as will be done frequently in the remainder of this chapter). However, it should be kept in mind that this is a simplication which needs to be justied from case to case.
9.1. Design and Operating Principle
93
t=0
N
N
N 2
N
2N
2N
2 2
N
2 2
N
t = /6
N 2N
N 2N
2 2N 2N
t = /3
N
N
N 2
N
2N
2N
2 2N 2N
Figure 9.3. Magnetic eld in the air gap of a stator with a solid, unwound rotor at dierent times. Left: one single-phase stator winding; right: three-phase stator windings. The top plots show the ampere windings in the individual slots; the bottom plots show the total linkage.
94
d
Figure 9.4. Geometric quantities used to calculate the magnetic energy stored in the air gap. The width of the air gap is assumed to be much smaller than the rotor diameter: D. Rotor and stator have the length l.
9.1.3
Magnetic Torque
The magnetic elds generated by the stator and rotor currents are superposed to a total eld B = Bs + Br (9.2) with p Bs = Bs sin( ( + t)) 2 r sin( p ( + t + )) Br = B 2 (9.3)
Both elds are assumed to be ideally sinusoidal and to rotate with the same angular velocity ; the rotor eld (together with the rotor itself) leads the stator eld by the angle . p is the number of poles. Using the geometry shown in Figure 9.4, the energy stored in the magnetic eld of the air gap is given by B 2 (, t, ) 1 dV = 20 20
2
Wm =
V
B 2 (, t, ) l
D d 2
(9.4)
The torque on the rotor is obtained by calculating the derivative of the energy with respect to the displacement angle between the rotor and the stator eld Wm lD M= = 40
2
B 2 (, t, ) d
(9.5)
Using several trigonometric identities together with the fact that only Br depends on , this can be simplied to p M = Mmax sin( ) 2 (9.6)
9.2. Stationary Operation
95
1.2 1 0.8 Rotor torque 0.6 0.4 0.2 0 20 0.2 40 60 80 100 120 140 160 180
Figure 9.5. Rotor torques in function of the rotor angle . Solid line: round rotor torque; dotted line: reluctance torque; dashed line: salient pole rotor torque.
where Mmax =
lD Bs Br 40
(9.7)
is the maximum torque that the machine can exert on the rotor. If the mechanical torque from the turbine exceeds this limit, the generator will fall out of step. If no protective measures are taken to reduce the turbine power, the rotor will accelerate until the centrifugal forces eventually destroy the generator! This is particularly dangerous in the case of turbo generators, that have a high speed to begin with. As noted before, the stator eld can only be regarded as sinusoidal for a turbo generator. In the salient pole generator, the inhomogeneous air gap gives rise to a reluctance torque Mr sin(p) which is added to the torque given in eq. (9.6). This is shown in Figure 9.5 for p = 2.
9.2
9.2.1
Stationary Operation
Stationary Single Phase Equivalent Circuit
If the rotor rotates at constant speed the rotor eld will induce sinusoidal voltages in the stator windings that are shifted by 120 with respect to each other. The frequency of the voltages is equal to the rotor frequency (multiplied by the number of pole pairs). The magnitude of the voltages depends on both the rotor speed and the magnitude of the rotor eld; it can therefore be controlled via the rotor eld current. Due to the symmetry of the voltages, the generator can be represented by the single phase equivalent circuit shown in Figure 9.6. In this circuit the induced voltages are modeled
96
jX d E
I U
Figure 9.6. Stationary single phase equivalent circuit.
by the voltage source E; I is the stator current, and U is the terminal voltage of the generator. With the losses neglected, the total impedance of the generator as seen from the stator can be modeled with a single lumped inductance Xd . This is strictly only valid for the turbo generator with its homogeneous air gap. In the salient pole generator, the same eects that give rise to the reluctance torque mentioned earlier also cause the active part of the stator current (with respect to the induced voltage E), to see a smaller reactance Xq than the reactive part, for which the reactance Xd is eective. This is illustrated in Figure 9.7. Typical values for Xd and Xq are listed in Table 9.1. For our purposes in this lecture we will neglect these dierences (but keep in mind that they exist). If the terminal voltage U is assumed to lie on the real axis, the phasor equation is (9.8) E = U + jXd I or in polar form E (cos + j sin ) = U + Xd I (j cos + sin ) (9.9)
where is the angle between E and the real axis counted counter-clockwise, and the angle between I and the real axis counted clockwise. Looking only at the imaginary part gives E sin = I cos Xd (9.10)
The active power of the machine (which is directly proportional to the mechanical torque) is given by P = U I cos Using eq. (9.10) this can be written as P = UE sin Xd (9.12) (9.11)
A comparison with eq. (9.6) reveals that for a machine with a single pole pair the mechanical angle between the magnetic elds of rotor and stator is identical to the electrical angle between the phasors of the induced voltage E and the terminal voltage U .

daxis
Iq jX I I
d
97
qaxis
qq
jX I
dd
Figure 9.7. Phasor diagram of voltages and currents of the salient pole machine.
Round Rotor Xd (p.u.) Xq (p.u.) 1.0 2.3 1.0 2.3
Salient Pole Rotor 0.6 1.5 0.4 1.0
Table 9.1. Typical values of reactances for synchronous machines.
9.2.2
Phasor diagram
Figure 9.8 illustrates the relation between the three phasors for dierent modes of operation. They can be categorized as follows. Generator operation (right column): || < 90 , > 0 The machine delivers active power. The stator current phasor lies in the left half-plane; the induced voltage phasor leads the terminal voltage phasor. Motor operation (left column): || > 90 , < 0 The machine consumes active power. The stator current phasor lies in the right half-plane; the induced voltage phasor lags the terminal voltage phasor. Overexcitation (top row): |E| cos > |U | The machine delivers reactive power (i.e. it acts like a capacitor). The stator current phasor leads the terminal voltage phasor; the real part of the induced voltage phasor is larger than the terminal voltage. Underexcitation (bottom row): |E| cos < |U | The machine consumes reactive power (i.e. it acts like an inductor). The stator current phasor lags the terminal voltage phasor; the real part of the induced voltage phasor is smaller than the terminal voltage.
98
Motor Im

Generator Im
E overexcited jX I
d
U E jXdI
U Re I Re
Im
Im
underexcited
E I U Re E jX I
d
jXdI I U Re
Figure 9.8. Phasor diagrams for dierent modes of operation.
As a special case consider a synchronous machine where the rotor is running idle. In steady state the mechanical torque is zero and therefore the machine neither consumes nor delivers active power ( = 0). Nevertheless, depending on the excitation it can either consume (underexcitation) or deliver (overexcitation) reactive power, i.e. = 90 . Such a set-up is called synchronous compensator and is used for variable reactive power compensation.
9.2.3
Operational Limits
Solving eq.(9.8) for I gives I= U E + jXd jXd (9.13)
If the terminal voltage U is considered constant, the stator current is determined by two parameters. According to eq.(9.11), the real part of the stator current is given by the turbine power. The imaginary part can then be controlled by adjusting the magnitude of the induced voltage via the rotor eld current. Eq. (9.13) is illustrated in Figure 9.9 in the P-Q plane. The limits of the operating area can be constructed as follows.

Im(I) Q Stator current thermal limit
99
Rated operating point Field current thermal limit
U/jXd
E/jXd
Re(I) P Practical stability margin (20)
Safety margin against 10% voltage drops Theoretical stability limit
Figure 9.9. Operational limits of the synchronous generator.
Stator current thermal limit: |I| < Imax The stator current limit is a circle around the origin with radius Imax . Field current thermal limit: |E| = f (If ) |E| < Xd U The eld current limit is a circle around jXd .
f (If,max ) Xd
Rotor angle stability limit: The theoretical rotor angle stability limit is determined by < 90 . It U is a parallel of the real axis that passes through jXd . In practice, a is kept to this limit. safety margin of 20 Safety margin against voltage drops: If the generator delivers little or no active power and operates close to the stability limit (underexcitation), even small voltage drops at the machine terminals could cause the rotor angle to exceed 90 . To allow for 10 % voltage drops, the operating point must stay above a parallel U of the real axis passing through 0.9 jXd . Depending on the design of the generator, there may be other limiting factors (e.g. for turbo generators the stator end region heating limit mentioned in section 2.5).
100
9.3
9.3.1
Dynamic Operation
Transient Single Phase Equivalent Circuit
During network transients, the reactance of the synchronous generator is not constant. For symmetrical transients, as discussed in section 7.2 on short circuits of synchronous machines, the machine reactance itself undergoes transient changes as the machine passes through the sub-transient, transient, and steady-state stages. In the fault analysis in chapter 7, only the rst few cycles after a fault were of interest. The generator was therefore modeled with the sub-transient reactance Xd . Compared to the time scale of short-circuit overcurrents, electro-mechanical oscillations are much slower (0.1 4 Hz as mentioned in the introduction to this chapter). For the study of these phenomena, the sub-transient phase is therefore neglected and the generator is modeled with the transient reac tance Xd . Similarly, during transients, the magnitude of the induced voltage drops to the value E which can be calculated as follows E = or where the superscript denotes the steady-state quantities prior to the transient. Thus, the phasor equation valid during network transients is
E = U + jXd I E = E 0 (Xd Xd ) I 0 sin 0 0 Xd Xd Xd E0 + U 0 cos 0 Xd Xd
(9.14) (9.15)
(9.16)
9.3.2
Simplied Mechanical Model
As stated above the rotor of a synchronous machine rotates with synchronous speed in steady state. If the electrical frequency of the system is e , the mechanical angular velocity m of the rotor is given by e m = (9.17) p/2 where p is the number of poles of the machine. An important parameter in the analysis of rotor oscillations is the total moment of inertia of the synchronous machine J. This is the sum of all moments of inertia of all rotating parts of the synchronous machine, i.e. the sum of the moments of inertia of the rotor, turbines, shafts and other devices on the shaft system, e.g. generator feeding the eld winding. As for electrical quantities it is practical to express J in a suitable p.u. base and therefore the inertia constant of the synchronous machine H is dened as H=
2 0.5Jmo S
(9.18)
9.3. Dynamic Operation Type of Synchronous Machine Thermal Power Steam Turbine Gas Turbine Hydro Power Slow (< 200 min1 ) Fast ( 200 min1 ) Synchronous Compensators Synchronous Motors Inertia Constant H (s) 49 7 10 23 24 1 1.5 2
101
Table 9.2. Typical values of H for dierent types of synchronous machines.
where S is the MVA rating of the machine. In eq. (9.18) the numerator is an expression for the total kinetic energy stored in the synchronous machine in steady state and the unit for H is thus seconds. (If there is a gear box in the system, it is of course dicult to dene one mechanical angular velocity. In such a case H is calculated as the ratio between total stored kinetic energy and the MVA rating of the machine.) The inertia constant states how long time it would take to bring the machine from synchronous speed to standstill if rated power is extracted from it while no mechanical power is fed into it. The value of the inertia constant will vary within a much smaller range than the value of J for dierent machines. Table 9.2 shows typical values of H for dierent types of synchronous machines. It can be concluded that the value is higher for thermal units as compared with hydro units. Typically 30 60% of the total moment of inertia comes from the turbines and shafts for thermal units, while the corresponding value for hydro units is 5 15%. For synchronous motors the inertia constant depends to a high degree of what kind of load that is connected.
102
10
The Swing Equation
This chapter contains a derivation and a discussion of the swing equation, which is the basic mathematical relation describing how the rotor of a synchronous machine will move (swing) when there is an unbalance between mechanical power fed into the machine and the electrical power extracted from it. Even if the basic physical relations are very simple, a number of simplications must be done to make the equation useful in practice.
HE SWING EQUATION is of fundamental importance in the study of power oscillations in power systems. Electro-mechanical oscillations is an important phenomenon that must be considered in the analysis of most power systems, particularly those containing long transmission lines, as e.g. the Swedish system. In normal steady state operation all synchronous machines in the system rotate with the same electrical angular velocity, but as a consequence of disturbances one or more generators could be accelerated or decelerated and there is risk that they can fall out of step. i.e. lose synchronism. This could have a large impact on system stability and generators losing synchronism must be disconnected otherwise they could be severely damaged. (A short description of power system protections is given in Appendix B.)
10.1
Derivation of the Swing Equation
The synchronous machine models developed in the previous chapter will form the basis for the derivation of the swing equation describing the electroTe Pm Tm m Generator Pe Pm m Te Tm Motor Pe
Figure 10.1. Schematic description of powers and torques in synchronous machines.
103
104
10. The Swing Equation
mechanical oscillations in a power system. Schematically the dierent torques and powers of a synchronous machine can be depicted as in Figure 10.1. Index m denotes in the following a mechanical quantity, and index e denotes an electrical quantity. The dierential equation describing the rotor dynamics is d2 m J = Tm Te (10.1) dt2 The quantities in eq. (10.1) are: J = The total moment of inertia of the synchronous machine (kgm2 ) m = The mechanical angle of the rotor (rad) Tm = Mechanical torque from turbine or load (Nm). Positive Tm corresponds to mechanical power fed into the machine, i.e. normal generator operation in steady state. Te = Electrical torque on the rotor (Nm). Positive Te in normal generator operation. If eq. (10.1) is multiplied with the mechanical angular velocity m one gets d2 m = Pm Pe (10.2) m J dt2 where Pm = Tm m = mechanical power acting on the rotor (W) Pe = Te m = electrical power acting on the rotor (W) If the angular acceleration should be expressed in electrical angle instead, eq. (9.17) is used to give 2 d2 e m J 2 = Pm Pe p dt where the left hand side can be re-arranged: 2 2 1 2 d2 e ( J) = Pm Pe pm 2 m dt2 (10.4) (10.3)
If eq. (10.4) is divided by the rating of the machine S, and eq. (9.17) is utilised again, the result is
2 2 ( 1 m J) d2 e Pm Pe 2 = 2 e S dt S
(10.5)
Observations and experiences from real power systems show that during disturbances, the angular velocity of the rotor will not deviate signicantly from the nominal values, i.e. from m0 and e0 ,respectively. This implies that eq. (10.5) together with the denition eq. (9.18) can be written as 2H d2 e pu pu = Pm Pe e0 dt2 (10.6)
10.2. Analysis of the Swing Equation
105
where superscript pu indicates that the mechanical and electrical powers should be expressed in p.u. of the rating of the synchronous machine. It is of course possible to use another base power than the rating of the machine in eq. (10.6), but this must then be made consistently regarding the denition of H and when calculating the power values in the right hand side. If another base power than the rating is used, the physical interpretation of H made in section 9.3.2 will not be valid anymore. Furthermore the typical values of H in Table 9.2 do not apply. In the following it will be assumed, if not otherwise explicitly stated, that electrical angles and electrical angular velocities are considered, and consequently the index e in the left hand side of eq. (10.6) can be omitted. It is also assumed that powers are expressed on the same base power as H, and the superscript pu can also be omitted in eq. (10.6). We are thus going to use the following form of the swing equation: 2H d2 = Pm Pe (10.7) 0 dt2
10.2
Analysis of the Swing Equation
Before applying the swing equation to a specic system, it could be of value to briey discuss the dierent terms of eq. (10.7) and their inuence. This will give an insight into the fundamental relations governing the dynamics during rotor oscillations. The dierence between mechanical power fed into the machine and the electrical output power will cause a motion of the rotor relative to a rotation with constant angular velocity 0 . 1 Therefore it is of interest to discuss in some detail how Pe and Pm will vary during a rotor swing. The mechanical power Pm is in most case provided by a hydro, steam or gas turbine. This power is determined by the gate opening of the turbine, and the time constant for changing the mechanical power is in most cases several seconds. There are dierent reasons why the mechanical power should be changed. For each generator there is a dispatch plan made for each hour. This plan is determined by the expected load in the power system. These plans are often the result of extensive optimisations of the resources available with due consideration to load, prices, cots, etc. The mechanical input powers are changed in accordance with these dispatch plans. But it is not possible to exactly forecast the power consumption in the system, Furthermore, generators could suddenly be disconnected due to faults, resulting in a shortage of power in the system. These unbalances will cause a frequency
In the following it is assumed that the synchronous machine is a generator. The same discussion could be applied to a synchronous motor, but then the signs of the electrical and mechanical powers must be changed. But since the dynamics of the system is dominated by the generators, it is natural to use these as the basis for type discussion. However, there is no principal dierence in the dynamics between a synchronous motor and generator.
1
106
deviation in the system. This frequency deviation can be used as an input signal to selected generators, which will change their mechanical input power and thereby their electrical output so that the frequency deviation is brought down to acceptable levels, typically 0.1 Hz in interconnected systems. The frequency controller could be rather fast, but since it takes several seconds to change the mechanical power Pm it will take some time for the frequency control to act. Generally it can be said that Pm = Pm () for those generators that are equipped with frequency controllers, but if time scales up to some seconds are considered it is a fair approximation to say that Pm = constant, at least for studies and discussion of more principal nature. In section 9.2 an expression for the active power of a synchronous machine connected to an innite bus was derived, eq. (9.12). As can be seen the power depends in the angle , and the rotor angle occurs thus on both sides of the swing equation. This applies for a synchronous machine connected to an innite bus, but it can be shown to be the case also for more complex systems with more synchronous machines. In this latter case, the electrical power of a given machine will not only depend on the value of its own rotor angle, but also on the rotor angles of other machines. The result is thus a system of coupled dierential equations. A more detailed analysis of the electrical power shows that it depends not only on but also on , that is on the relative angular velocity as compared with a synchronously rotating system. This contribution, which only occurs during transients, is due to currents induced in the rotor circuits, and it tends to damp out the oscillations. A voltage controller might also give a contribution to the electrical power depending on . With fairly good accuracy the electrical power can written as Pe = Pe (, ) = Ps () + Pd () (10.8)
with Ps being the synchronising power and Pd being the damping power. This is further explained in section 11.4.
10.3
Swing Equation as System of First Order Dierential Equations
The swing equation (10.7) is an ODE of the second order. 2 ODEs of higher orders can be written as a system of rst order ODEs, which in many cases is practical. Particularly for multi-machine systems this turns out to be an attractive approach. In most cases it is not possible to solve the swing equation analytically, but one has to use numerical integration in solving it. Numerical integration requires usually that the ODE is in the form of ODEs of rst order.
2
The abbreviation ODE for Ordinary Dierential Equation is here introduced.
10.3. Swing Equation as System of First Order Dierential Equations 107 By introducing = eq. (10.7) can be written as = 0 (Pm Pe ) 2H (10.9)
This form of the swing equation will also be used in the following. The quantities and in eq. (10.9) are named state variables, or just states, and eq. (10.9) describes the system in state form or standard form. The vector [ ]T is the state vector of the system. Some clarications will be made concerning the angular velocity introduced in eq. (10.9) above. This angular velocity denotes the frequency with which the rotor oscillates relative to a system rotating with the synchronous and constant angular velocity = 0 . The absolute angle of the rotor abs relative to a reference system at time t is given by abs = 0 t + + 0 (10.10)
where 0 is the angle of the rotor at t = 0. Dierentiation of eq. (10.10) with respect to time gives abs = 0 + (10.11) from which it can be concluded that the absolute angular velocity of the rotor abs is given by abs = abs (10.12) and the angular velocity relative the synchronously rotating system is given by = (10.13) It is this latter relative angular velocity that is of interest when rotor oscillations are studied, while the absolute angular velocity is of interest when studying frequency stability.
108
11
Power Swings in a Simple System
In this chapter the models and equations earlier derived are applied to the simple system of a synchronous machine connected to an innite bus. Despite the simplicity of the system, useful and general information can be gained from this analysis. Through qualitative discussions the features of the equilibrium points of the system can be determined. Time simulations of the system show the principal behaviour of stable and unstable solutions. The equal area criterion is derived. This is a very powerful method for determining the stability of a system. Another powerful tool is provided by small signal analysis, or linear analysis, which gives valuable information about the local behaviour of the system. The chapter is concluded by a discussion of dierent ways of improving the angular stability of a power system.
COMPLETE STABILITY analysis of a power system is an extensive and complicated task. However, it turns out that many of the most important phenomena and mechanisms can be found in very simple systems, where they can be seen very clearly. In large and complicated systems it is often hard to distinguish the fundamental and decisive phenomena from the more irrelevant ones. It is therefore of importance to study simple systems to get an insight into and understanding of the basics, that can be used in the analysis of more complex systems. This chapter focuses on power oscillations in the simple system of Figure 11.1. This system can be a model of a synchronous machine or group of synchronous machines connected to a larger system through a one or more power lines. The reactance Xe in Figure 11.1 is here an equivalent reactance including transformers and parallel lines. Even if this system is simple, a number of simplication are required to get simple solutions to the system. The reader is encouraged to solve problems from the problem set to get a better understanding of the issues discussed here.
11.1
The Swing Equation and its Solutions
The solutions to the swing equations of the simple system in Figure 11.1 will be analysed and discussed in this sub-section. Despite its simplicity, a number of important conclusions concerning the angular stability in large 109
110
11. Power Swings in a Simple System
systems can be drawn from this system. Earlier a number of simplications and assumptions have been made, more or less explicitly, and for the sake of completeness and clarity they are summarised here: 1. The synchronous machine is modelled as a constant emf behind the transient reactance Xd . The angle of the emf is assumed to coincide with the rotor angle. 2. Resistances in lines, transformers, and synchronous machines are neglected. 3. Voltages and currents are assumed to be perfectly symmetrical, i.e. pure positive sequence. 4. The angular velocity is close to the nominal one. 5. Static models for lines are used. 6. The mechanical power Pm , i.e. the power from the prime mover, is constant during the transient under study. Furthermore, it is assumed that the damping power of the system can written as Pd = D (11.1) An electrical equivalent to Figure 11.1 is shown in Figure 11.2. Together with the assumptions above the swing equation can now be formulated as 2H d2 = Pm Pe 0 dt2 with Pe = and Pm = Pm0 = konstant Now we introduce Pe,max =
Eq UN Xd + Xe Eq UN sin + D Xd + Xe
(11.2)
(11.3)
(11.4) (11.5)
Xe
Infinite bus U N 0
Figure 11.1. Synchronous machine connected to innite bus.
11.1. The Swing Equation and its Solutions
111
Xd E q
Xe U N 0
Figure 11.2. Equivalent electric circuit of a synchronous machine connected to an innite bus.
and eq. (11.2) can be written as 2H d2 = Pm0 Pe,max sin D 0 dt2 (11.6)
which is the complete swing equation of the simplied system with the introduced assumptions.
11.1.1
Qualitative Analysis
In order to make a qualitative analysis of the solutions to eq. (11.6) a further simplication will be made, i.e. the damping is neglected and D is thus set to zero : 2H d2 = Pm0 Pe,max sin (11.7) 0 dt2 Of fundamental importance to a non-linear system is its equilibrium points, i.e. the points in state space where all time derivatives vanish. For the equilibrium points the right hand side of eq. (11.7) is thus zero resulting in the following equation for 0 Pm0 = Pe,max sin 0 (11.8)
Figure 11.3 shows how Pe and Pm vary with rotor angle , and the following conclusions about the equilibrium points can be drawn from this gure: 1. If Pm0 < Pe,max there are two equilibrium points, i.e. 0 and 0 for 0 . 2. If Pm0 = Pe,max there is exactly one equilibrium point 0 = /2 for 0 . 3. If Pm0 > Pe,max there is no equilibrium point. (The angle 0 is given by 0 = arcsin(Pm0 /Pe,max ).)
112

P max sin e,
P() e P m0
Figure 11.3. Diagram showing the variation of electric and mechanical power for the system in eq. (11.7).
It is clear that if Pm0 > Pe,max the system is unstable and further analysis is superuous. In this case the rotor will be accelerated until the protections trip the generator and turbine. A necessary condition for stability is that at least one equilibrium point exists. In the case of Pm0 < Pe,max the equilibrium point = 0 is stable for (suciently) small disturbances, which implies that that if the system is moved away from 0 the dynamics of the system tends to bring it back to 0 . This is veried by the following qualitative reasoning: If the rotor has been accelerated so that > 0 , the right hand side of eq. (11.7) will be negative and the system is decelerated and starts to move back to 0 . Correspondingly a deviation < 0 will give an acceleration so that is brought back to 0 . (As seen from Figure 11.3 these arguments are valid only for suciently small deviations from 0 och = 0.) It is seen that if > 0 and if 0 then the system will never return to = 0 . For a state vector [ ]T in the two dimensional state space, or phase space as it is sometimes referred to, that initially is close to [ ]T = [0 0]T , will be close to this point for all t > 0. As in a real system there is positive (hopefully) damping, so the solution will converge (asymptotically) to [0 0]T . The region that above has been described as close to [0 0]T is
113
called the region of attraction of this equilibrium point, and the curve that separates this region from points that will not converge to [0 0]T , is called the separatrix. This curve separates the stable and unstable solutions from each other. A similar discussion concerning the equilibrium point = 0 shows that this point is unstable. A stability criterion for the system in Figure 11.3 can now be formulated:
Stability Criterion The system i Figure 11.3 is stable, i.e. the synchronous machine will not fall out of step, if it after a disturbance is in the region of attraction of the equilibrium point [0 0]T . 1 The criterion above implies that the system must have decelerated so that = 0 before the point = 0 is reached, otherwise the generator will fall out of step. The discussion above has disclosed a number of important features of the system without that any equation has been solved. The solutions to eq. (11.7) are not easy to solve analytically. It is seen that eq. (11.7) is the same as the equation describing a mathematical pendulum, and the exact solutions are given by elliptic integrals. Usually a mathematical pendulum is studied by its linearised equations, but one cannot get all information regarding stability from the these equations. The linearised equations can however give other kinds of useful information, and we will come back to that in section 11.4. In the next section the behaviour of the stable and unstable solutions will be discussed.
11.1.2
Stable and Unstable Solutions
In order to verify the conclusions drawn from the qualitative discussion above simulated solutions to eq. (11.7), i.e. for the system in Figure 11.1, will be shown for a few cases. The disturbance considered here is a three phase to earth fault on one of the lines close to the generator. That the fault is close to the generator implies that the electric power during the fault is zero. The faulty line is disconnected be be the distance protections, see Appendix A, causing that the reactance Xe is changed when the fault is cleared. By varying the fault clearing time dierent types of solutions to the swing equation can be obtained.
1 It is clear that = 0 is the only possible stable equilibrium point. P () is periodic with the period 2, but if the system has moved so that is close to another stable equilibrium point, will be so large that it is outside the region of attraction of this equilibrium point.
114

Electrical Power Before and After the Fault 2
1.8
1.6
1.4 Electric Power (p.u.)
1.2
0.8
0.6
0.4
0.2
20
40
60
80 100 Rotor Angle (degrees)
120
140
160
180
Figure 11.4. The electric power before the fault, Pe = 1.86 sin , after the fault, Pe = 1.26 sin , and the mechanical power, Pm = 1, as function of the rotor angle, .
The system is thus described by the following equation 2H d2 = Pm0 Pe,max sin D 0 dt2 (11.9)
The electric power is as follows: Pre-fault Pe,max = 1.86 p.u., post fault Pe,max = 1.26 p.u. and during the fault Pe = 0. Power curves before and after fault are shown in Figure 11.4 where also the constant power from the turbine Pm = 1 p.u. is drawn. Before the fault = 32.5 , which corresponds to the left intersection between Pe = 1.86 sin and Pm = 1. During the fault the rotor will accelerate since Pe = 0 and consequently Pm Pe > 0. When the fault is cleared the electric power will follow the curve Pe = 1.24 sin . If the system is stable it will settle down to the point = 53.5 which is the left intersection between Pe = 1.24 sin and Pm = 1. When the rotor angle is between 53.5 and 180 - 53.5 = 126.5 , which corresponds to the right intersection between Pe = 1.24 sin and Pm = 1, then Pm Pe < 0 and the rotor will decelerate. Should the rotor move beyond this point, then Pm Pe > 0 and the rotor will accelerate and the system becomes unstable. It is often said that the generator falls out of step or loses synchronism.
115
The longer the duration of the fault, the more the rotor will be accelerated and the larger the rotor angle. If the duration of the fault is too long there is a risk that the rotor swing is so large that it passes the right equilibrium point and the synchronous machine loses synchronism. To verify and demonstrate this the system has been simulated for fault clearing times equal to 4, 6, 6.5, and 8 cycles. 2 Results from these simulations are shown in Figures 11.5 to 11.8 where both the swing curves and the phase portraits of the system are shown. The phase portrait, or phase trajectory, is the trajectory of the system in the ( )-plane during the transient. The damping D is set to 0.02. The value of the damping determines how fast a stable solution will converge to its equilibrium point, and this value has minor inuence on the rst swing directly after the fault. 3 With a fault clearing time of 4 cycles the system is clearly stable, see Figure 11.5. The maximal rotor angle is 90 which gives an ample margin to the critical value of 126.5 . The system will eventually converge to the point = 53.5 . If the fault clearing time is increased to 6 cycles, see Figure 11.6, the stability margin is much smaller. The rotor swings out to 120 and it is very close to the critical point where it will start to accelerate. If the fault clearing time is increased to 6.5 cycles the curves of Figure 11.7 are obtained, and here the system is unstable. The generator cannot be decelerated enough before it reaches the critical point, but it passes this point and is further accelerated and loses synchronism. This type of instability is called rst swing instability. The phase portrait in this case shows that rotor is close to be decelerated before reaching the critical point, but it never reaches a zero value of . The previous fault clearing time, i.e. 6 cycles, is called critical fault clearing time. The last case shows a simulation with the fault clearing time equal to 8 cycles, see Figure 11.8. In this case the rotor is accelerated for such a long time that it passes the critical point very soon. The deceleration period could hardly be seen in the swing curve. The phase portrait shows also that the deceleration just after fault clearing is clearly insucient to stabilise the system.
Fault clearing times are often given in cycles of the power frequency. The reason for this is that part of the clearing time comes from the time needed for the breakers to clear the fault. Since the fault current only can be intercepted at a zero crossing, it is natural to measure the breaker time in cycles. Modern high voltage breakers can break a fault current in 2 3 cycles. Should still faster breaker times be required, special breakers could be installed, but these are more expensive than standard breakers. Protections are also described in Appendix A. 3 This is the case for realistic values of D. If the damping is very large, and positive, this will increase the stability signicantly in the system. The value used here is typical in a system where no special equipment has been installed to increase the damping, e.g. Power System Stabilisers on the voltage regulators.
116

Swing Curve for Generator Against Infinite Bus 150
100 Rotor Angle (degrees) 50 0 0
0.5
1.5 Time(s)
2.5
Phase Portrait of Generator Against Infinite Bus
4 Relative Angular Velocity of Rotor(rad/s)
6 0 50 Rotor Angle (degrees) 100 150
Figure 11.5. Swing curve and phase portrait for generator swinging against an innite bus. Fault clearing time = 4 cycles. Critical fault clearing time = 6 cycles.

117
0.5
1.5 Time(s)
2.5
118

0.5
1.5 Time(s)
2.5
Figure 11.7. Swing curve and phase portrait for generator swinging against an innite bus. Fault clearing time = 6.5 cycles. Critical fault clearing time = 6 cycles.

119
0.5
1.5 Time(s)
2.5
120
11.2
Equal Area Criterion
It was concluded above that even for quite simple systems, like the one in Figure 11.1 it was hard to get analytical solutions to the system equations. But quite often the exact behaviour of the solutions is not of interest, but most important is to determine if the system is stable or not after a given contingency. By using the equal area criterion that will derived in this section, it is possible with rather simple calculations to investigate the stability of the system in Figure 11.1 for dierent disturbances. The equal area criterion can be extended to more complicated systems, so that preliminary stability analyses could be done very fast. The swing equation of the system in Figure 11.1 can be written as d2 0 (Pm Pe,max sin ) = dt2 2H or 0 d2 = Pa 2 dt 2H with the accelerating power Pa dened by Pa = Pm Pe,max sin (11.10)
(11.11)
(11.12)
As the interest here is focused on the dynamics during the rst swing, the damping can be neglected, i.e. D = 0. From the Figures 11.5 to 11.8 it can be concluded that a necessary, but not sucient, condition for stability is that there exists a moment in time tm during the swing such that (tm ) = 0. The corresponding angle is m . Some formal manipulations with eq. (11.11) together with the condition (tm ) = 0 will provide a stability criterion. If eq. (11.11) is multiplied with one gets
2 0 d Pa 2 = dt 2H
(11.13)
which can be written as 1d 2 dt d dt

2
d dt
0 2H
Pa d
(11.14)
that can be integrated to give d = dt 0 H
Pa d + C
i
(11.15)
with C being a constant of integration, which is 0, since = 0 when = i . (i is the pre-fault rotor angle.) Thus a necessary condition for stability is that there is an angle m such that 0 H
m i
Pa d = 0
(11.16)
11.2. Equal Area Criterion or
121
m i
Pa d = 0
(11.17)
The stability criterion can thus be formulated as: The system is stable if there exists an angle m such that the area below the accelerating power Pa in the P () diagram between i and m vanishes. In most cases of practical interest the dynamics could be divided into two dierent phases: A rst phase when the rotor is accelerated, and a second phase when it is decelerated. Assume that the rotor is accelerated up to = c and that it is decelerated when c m . Then two dierent areas in the P () plane can be dened as
c
Aa =
i
(Pm Pe ( ))d (Pe ( ) Pm )d
(11.18)
Ar =
c
(11.19)
Indices a and r denote here accelerating and decelerating (retarding) areas, respectively. The angle c is most conveniently chosen as the angle when the fault is cleared. The stability criterion can now be formulated as The system is stable if there exists an angle m such that the areas Aa and Ar are equal, i.e. Aa = Ar . This latter formulation have given rise to the name equal area criterion for this stability criterion. This is illustrated in Figure 11.9. The advantage with the equal area criterion is that the stability of a system can be investigated without any excessive computational eorts. But the price paid is that the time t was eliminated from the equations. Since the time does not appear explicitly in the equations, actions to improve the stability must formulated in the angle space, which usually is not so practical. It is of course of more interest to have these actions formulated in terms of times, e.g. critical fault clearing times. By the equal area criterion critical fault clearing angles could be formulated, and these must then be transformed to critical fault clearing times. Sometimes this can be done easily, sometimes it is not straightforward. In the second formulation of the equal area criterion, an angle c was introduced. This angle denes for which values that Pa should be calculated as positive and negative, respectively. The choice of c has no inuence on the result from the equal area criterion, which is obvious from the rst formulation. When applying the equal area criterion in a practical case it is natural to chose c as the angle when the fault is cleared. At this angle the
122

P max sin e,
P() e Ar P m0
Aa
Figure 11.9. Application of the equal area criterion after a disturbance.
topology of the system is often changed, and thus the expression for Pe (), resulting in dierent expressions for Pe in the computation of Aa and Ar . Dierent stability related problems can be solved by using the equal area criterion. One common application is the calculation of how fast a fault must be cleared to ensure stability of the system. The maximum fault clearing time for which the system remains stable is called critical fault clearing time, as mentioned above, and the corresponding angle is called critical fault clearing angle. The equal area criterion can also be used to calculate the maximum power that can be transmitted for a given fault scenario. A third application is to determine if a system is stable or not for given data and disturbances. It is strongly recommended that problems from the problem set are solved to get an understanding of the application of the equal area criterion.
11.3
Lyapunov Stability Criterion
The theory of the Russian mathematician A. Lyapunov (1857 -1918) oers powerful tools for stability analysis of dynamical systems. In this theory the following theorem, rst formulated in 1899, forms the basis.
11.4. Small Signal Analysis
123
Theorem. Consider the dynamical system x = f (x), where x0 is an equi librium point. If in a region around the point x0 there exists a real valued and continuously dierentiable function V such that 1. V(x) > 0 and V(x0 ) = 0, 2. V(x) = grad(V(x)) f (x) 0, then x0 is a stable equilibrium point. If V(x) < 0, then x0 is asymptotically stable. The existence of the function V is also a necessary condition for stability. A problem with Lyapunovs theorem is that it does not give any information on how to nd the function V. It is obvious from the theorem that V can be regarded as a form of energy, and the theorem says that if the system is stable, it will stay in the region around x = x0 where V is positive denite provided the energy introduced to the system by a disturbance is suciently small. For the system studied above the following function can be shown to fulll the requirements of Lyapunovs theorem V(, ) = with Vk = and Vp = Pm0 (0 ) + Pe,max (cos 0 cos ) (11.22) Vk can here be interpreted as kinetic energy and Vp as potential energy. It is straightforward to show that V = D 2 by using eq. (11.6), and consequently the second requirement of the above theorem is fullled if D 0. It is clear that V is positive denite in . In Figure 11.10 the potential energy, Vp is shown for a system with Pm0 = 0.3 and Pe,max = 1 and it seen that it is positive denite in a neighbourhood of the equilibrium point 0 = arcsin(Pm0 /Pe,max ). Thus if the energy injected into the system during the fault will be less than the potential energy at the local maximum of Vp , the system will be stable. It can be shown that the above Lyapunov function, V, gives the same stability criterion as the equal area criterion. H 2 + Pm0 (0 ) + Pe,max (cos 0 cos ) = Vk + Vp (11.20) 0 H 2 0 (11.21)
11.4
Small Signal Analysis
The equal area criterion oers a tool to make stability analyses of the system in Figure 11.1 by rather simple computations.. For this system it is also possible to do other types of analyses that do not require any more extensive
124

Potential energy part of Lyapunov function
1.8
1.6
1.4
1.2
Vp
0.8
0.6
0.4
0.2
0.2 100
50
50
100
150
200
theta (degrees)
Figure 11.10. The function Vp of eq. (11.22).
computations. One method often used to analyse non-linear system is to linearise the system and study this latter system. Linear analysis is usually not so computationally demanding, and there exists many dierent tools for linear analysis. If the non-linearities of the system are small, it is clear that the linear system is a good approximation of the non-linear one. But it can be shown that for non-linear systems in general, provided the system equations are locally dierentiable, the linear system would give valuable information about the non-linear system, e.g. regarding its stability. Small signal, or linear, analysis will now be applied to the system in Figure 11.1. The system in Figure 11.1 is described by x1 x2 with [x1 x2 ]T = [ ]T It is now assumed that the system is in an equilibrium point [x10 x20 ]T = [0 0]T (11.25) (11.24) = x2 0 (Pm0 Pe,max sin x1 Dx2 ) 2H (11.23)
11.4. Small Signal Analysis
125
K M
Figure 11.11. Mass and spring in analogy with the system in Figure 11.1.
The linearised system describes the behaviour of small oscillations around the equilibrium point eq. (11.25). A rst order Taylor expansion of eq. (11.23) about the equilibrium point gives x1 x2 = 0 1 0 0 Pe,max cos x10 D 2H 2H x1 x2 (11.26)
From eq. (11.26) the frequency of small oscillations about the equilibrium point can be calculated and they are given by the eigenvalues of the 2x2 matrix, the Jacobian matrix, in the right hand side of eq. (11.26). First it is assumed that there is no damping in the system, i.e. D = 0. The eigenvalues of eq. (11.26) are then purely imaginary, which corresponds to (undamped) oscillations if 0 0 /2, = 0 Pe,max cos 0 2H (11.27)
and the frequency of oscillation is given by p = 0 Pe,max cos 0 2H (11.28)
Eq. (11.28) is similar to the equation describing the motion for a mass connected to a wall with a spring according to Figure 11.11. The frequency of oscillation for this system is given by KM = K M (11.29)
By comparing eq. (11.28) with eq. (11.29) it is seen that the spring constant K in the mechanical system corresponds to Pe,max cos 0 , which sometimes is called the stiness of the electrical system or its synchronising power. Furthermore, M corresponds to 2H/0 . From eq. (11.27) it is also seen that if /2 < 0 < , the eigenvalues will be real, and one of them will be positive corresponding to an unstable mode. This was earlier also concluded by the qualitative discussion in section 11.1.1, and this is now veried by the mathematical analysis done here.
126
It was earlier pointed out that the damping is not signicant for large disturbances in the system. For large disturbances it is the ability of the system to stay within the region of attraction of the equilibrium point that is important. This is depending on the stiness of the system. The damping is more important for how the system settles down to the equilibrium point when it has survived the rst swing. The linearised could thus be used for studies of the damping in power systems. For the simple system studied here the eigenvalues with D = 0 become D = j 2 with p given by eq. (11.28) and D = 0 D 2H (11.31) (p )
2
D 2
(11.30)
As the damping in a power system usually is quite small, eq. (11.30) can be approximated by D jp (11.32) 2 The frequency of oscillations can also in this case with good accuracy be approximated by eq. (11.28), and it is thus determined basically by the stiness of the system and the inertia constant of the machine.
11.5
Methods to Improve System Stability
It is obvious that its desirable that a power system can withstand as many disturbances as possible without becoming unstable. However, it is realised that it is not possible to design a system that can cope with all conceivable contingencies, so one has to restrict the considered disturbances in the design of the system. Usually, one considers only the most frequent faults in the system. It is also important to consider the consequence of a fault or disturbance when designing a system. Generally it can be said that the more investments that can be made in a system, the more robust the system can be made. One is thus faced with a technical-economical optimisation problem that is very complex. Furthermore, there are a number of parameters that are both very important in the optimisation process and at the same time very hard to quantify. Such a parameter is the value that consumers put on uninterruptable supply of electric power of good quality. This optimisation process has gained a lot of interest during the last few years due to the liberalisation (de-regulation) of the power market taking place all over the world. Theoretically the optimum value of reliability of power supply is when the marginal increase in the value perceived by the customers is equal to the cost of the investment made to achieve this increase in reliability.
11.5. Methods to Improve System Stability
127
Practically it is of course impossible to deduce this optimum, and one has to judge by dierent methods if the consumers are satised with the reliability of the system, and if not if they are willing to pay more to improve it. To have some more clear and practical rules for system design a number of deterministic rules have been used over the years. One such rule is the so called (N 1)-rule, which says that a system should be able to supply customer loads with any component, generator, line, etc., out of service. The disconnection of the component should preceded by a fault and the transients triggered by this fault should be considered. A more conservative approach would be to apply an (N 2)-criterion, which would result in a more reliable but more expensive system. Modern methods apply probabilistic tools. In many of these approaches one considers both the probability of a fault and its consequences. By assigning a measure, i.e. a number, for this combination of fault probability and consequence a risk based security analysis can be made. Some of these aspects are also discussed in chapter 14. This section will be concluded by a discussion of dierent ways to improve the angular stability in power systems. The stabilising eect of the dierent methods can be veried by studying Figures 11.3 and 11.9. Which method to apply in a given situation depends on a large number of parameters and must be determined from case to case, usually after extensive studies. In some cases, one or several methods can be excluded, e.g. the construction of new lines since the needed permission cannot be granted from the authorities. The most common methods to increase system angular stability are: Increase of the inertia constant of the generators. This makes the rotors more dicult to accelerate in connection with faults and the risk for losing synchronism is reduced. In most cases this is a very expensive means and only in special cases it can be applied, e.g. by installing a ywheel on a small hydro unit. Increase of system voltage. This increases Pe,max and for given power Pm the stability margin is increased. Reduction of the transfer reactance Xe . This will also increase Pe,max as in the previous case. This can be achieved by constructing parallel lines, or by installing series capacitors on existing lines or new lines. By installing series capacitors the eective reactance of the line is reduced. This method has been used extensively over the years, e.g. in the Swedish system. Installation of fast protections and fast breakers. In this way the time with a fault connected can be reduced and thereby the time during which the generator rotors are accelerated. The ability for the system
128
11. Power Swings in a Simple System to decelerate the rotor swings is increased. Another way is to use automatic re-closure after the fault is cleared, see Appendix A. Implementation of fast valving in steam turbines. By fast control of the mechanical power during and after a fault, the acceleration of the rotors can be decreased. It cannot be applied to nuclear power plants by security reasons. The method has not been used to any larger extend, since it is claimed to impose large thermal and mechanical stresses on valves, turbines, etc. Installation of braking resistors. These are shunt resistors that are connected by breakers fast after a fault close to a critical generator. The electric load of the critical machine increases and the risk for losing synchronism is reduced. Stability control of controllable devices such as High Voltage Direct Current (HVDC), controllable series capacitors, controllable reactive shunt compensation (SVC), etc. These devices are usually to expensive to install just for stabilising the system, but when they are installed, the cost to use their controllability for stabilisation is usually marginal.
12
Oscillations in Multi-Machine Systems
This chapter contains a very brief overview of modelling of multi-machine power systems. Many of the issues mentioned here will be dealt with in the advanced course on power system analysis
12.1
Simplied Model
The simplied model considered in the previous section has of course a limited application to realistic systems. The analysis of such a simple model is primary motivated by the insight it gives, and by that the model can be used for simplied, preliminary calculations. In this section we will expand the model to a system with more than one machine. The same assumptions for modelling is done, i.e.: The damping is neglected. The mechanical power is constant. The synchronous machines are modelled as constant voltage sources behind the transient reactance. The phase angle for the voltage source coincides with the rotor angle. Loads are represented by impedances (admittances). A model based on these assumptions is usually referred to as the classical model for electro-mechanical oscillations. Since the loads are represented by impedances the load buses can be eliminated and included in the node-admittance matrix, and the power which ows out from the generators can be written as a function of the voltage magnitudes and phase angles of the dierent machines. Pei = Pei (E, ) (12.1)
where E is a vector of voltages, of the N generators, and is a vector of their rotor angles T = 1 2 . . . N (12.2) 129
130
12. Oscillations in Multi-Machine Systems
and analogously for E. The swing equation for machine i can now be written as i i 0 = (12.3) (Pmi Pei ) i 2Hi To get a compact form of this system of equations we introduce x= 1 1 . . . N N
T
(12.4)
and the N coupled equations in (12.5) can be written as x = f (x) (12.5)
where f is a vector valued function containing the swing equations for the dierent machines. This system of swing equations is not easily analysed in the general case. In most cases numerical methods must be used for studies of the solutions. It can also be done by studying the linearised system according to the method in ection 11.4. It is found that the system can be compared to a mechanical system with N masses connected with springs and with spring constants determined by the stinesses of the lines. The eigenvalues and eigenvalues of the Jacobian matrix give the swing modes of the system, which often give a very good picture of how the real system behaves at disturbances. Actual oscillations in the system can often with good accuracy be approximated with linear combinations of these swing modes. Another way to analyse the system is by solving the non-linear equations with numerical integration (time-simulations).
12.2
Detailed Model
A number of simplications were earlier introduced to facilitate the analysis of the equations by utilising well known mathematical methods. Unfortunately, the used model has some important shortcomings which makes it less appropriate for more accurate studies of the dynamics of the power system. The most serious shortcomings are: The load model is too simple. A large part of the loads consists of asynchronous machines and other drives and their characteristics do not follow that of an impedance at voltage oscillations. In principle, the loads consume as much power as the generators produce. Therefore it is equally motivated to develop realistic load models as to make detailed models of the synchronous machines. The dynamics of the synchronous machines is too simple. Here it is only determined by two states, i and i , while a realistic model needs more states for a satisfactory description. Dierent uxes in the
12.2. Detailed Model
131
synchronous machines will vary with time and dierent regulators, e.g. for voltage and frequency, will have inuence on the process. Saturation phenomena in generators can also have inuence. Modern controllable equipment like static compensators (SVC) and high-voltage direct current transmissions (HVDC), will in many cases have a great inuence on the stability of the system. Therefore these should be modelled in an adequate way. The shortcomings originate essentially from a simplied modelling of important components in the system. A more detailed description gives naturally a better description, but the overview and handiness of the problem are at the same time lost. More complicated models imply that more state variables must be introduced. In the classical model every generator is modelled by two states, i and i , while a more detailed model demands more states. To model exciters, voltage regulators and turbine regulations more states also need to be introduced. Components like SVC and HVDC, with associated control equipment generally also require dynamic modelling, which imply that new state variables must be introduced. For more realistic load models, the load buses cannot be eliminated, as was done in the classical model. This implies that the equations for loads will formulate a constraint that the state variables must fulll. The above outlined general model for electro-mechanical oscillations have the following properties, which are worthwhile to comment: It is dierential-algebraic. The classical model resulted in a pure differential equation. Dierential-algebraic systems are much more complicated to analyse than systems which can be modelled by a pure system of dierential equations. The system is non-linear. The system is non-time invariant, i.e. the time occurs explicitly in the equations. This depends on events in the system, e.g. disturbances and connections, but even on protective actions of dierent kind. For realistic power systems the number of state variables is very large. As an example, the Nordic system contains hundreds of generators and the total number of states at simulations can be up to several thousands. The used analysis methods are based, either on the linearised system or on timesimulations of the complete model.
132
12. Oscillations in Multi-Machine Systems
13
Voltage Stability
Voltage stability has been a limiting factor in many power systems during the last years. Due to a gradual increase in load the reactive power balance can be critical at peak load conditions. This chapter summarises the basic concepts and provides basic analyses for two simple test systems. A number of stability indicators are introduced an applied to one of the test systems. The importance of load charactestics and load dynamics are demonstrated.
N THE PREVIOUS chapters it was shown that an active power imbalance could give rise to rotor oscillations, which can cause instability if the disturbances are large enough. If all generators in a system rotate with the same electrical angular velocity (in average over a certain time period) the system is said to be synchronously stable. (Sometimes it is said to be in angle stability). This type of stability is strongly connected to the active power balance in the system. As has been seen from the static analysis the reactive power and the voltage are strongly coupled to each other. As a rule of thumb it can be said that if there is net production of reactive power in a node, the voltage is high in the node, while a decit of reactive power implies that the voltage is low. It could now be asked whether the interplay between Q and U can also cause instabilities similar to the angular instability discussed earlier. It turns out that this is the case. The underlying mechanisms for this voltage instability will be dealt with in this chapter together with some methods of analysis.
13.1
Mechanisms of Voltage Instability
In chapter 8 voltage instability was classied into slow and fast ones. Even if the physical mechanism is the same, i.e. a lack of reactive power, the power system contingencies and power system structures causing them are dierent. The nal break down, collapse, in the system could be quite fast also in a slow voltage instability, and the actual instability can sometimes be disguised, which makes this kind of instability troublesome to detect in a real system. Consequently it can be hard to get reliable indicators to activate preventive actions if the time needed is available. 133
134
13. Voltage Stability
13.1.1
Long Term Voltage Instability
The time scale of the long term, or slow, voltage instability could be tens of seconds to hours. It develops from a gradual lack of reactive power at a node or in a part of the system. Even if the voltage stability phenomenon is closely coupled to the reactive power balance the active power plays a very important role indirectly. Many loads have an inductive power factor that is not unity, and thus an increase in active power will cause an increase in reactive power consumption. Furthermore, the reactive losses in transmission lines are proportional to the square of the active power transmitted, so heavily loaded lines will inuence the reactive power balance signicantly. A very deceptive role is sometimes played by synchronous machines and other sources of reactive power. These devices can often supply more reactive power than their steady state rating for a limited time, sometimes up to several minutes, and during this period of time the reactive power balance and voltage could be kept at acceptable values. However, when protections in these devices act to avoid that the equipment is damaged, the reactive power output is drastically reduced with a very fast voltage collapse as a result. The contributing factors to a slow voltage collapse are thus: Load increase Load recovery after faults Reactive losses in line due to high power transfer (This could be caused by line outages due to contingencies.) Loss of reactive supply (Capacitors, synchronous machines, etc.) For the second bullet above, the control and dynamics of load tapchangers are very important. A voltage reduction in the high voltage system will propagate to the lower voltage levels. This will usually result in a reduction of the load powers. However, load tap-changers on the transformers will try to restore the voltage on the low voltage side and thereby increasing the load power. This has in most cases a de-stabilising eect on the voltage stability.
13.1.2
Short Term Voltage Instability
For the short term, or fast, voltage instability the actual voltage collapse comes very soon after the disturbance, just a few seconds or less. In these instabilities, loads or other devices with special reactive power characteristics are often involved. Such a load is induction motors that just after fault consume a lot of reactive power. Other loads that could be very detrimental in this context are some types of power electronics based equipment. As an
13.2. Simple Systems for Analysis of Voltage Stability
135
Infinite bus
Ze
Ul-f Load
UN0
Pl + jQl
Figure 13.1. Simple system for analysis of voltage stability.
f example could a line commutated HVDC converter be de-stabilising in weak system, and if not properly accounted for in the system design it could lead to system collapse for critical contingencies. t Induction motors aretroublesome also in the slow voltage instability. This is particularly the case in countries where the peak load contains a large fraction of air conditioning, which is driven by induction motors. If the voltage gets low, these could stall and hence implying a higher reactive power consumption.
13.2
Simple Systems for Analysis of Voltage Stability
The simple model in Figure 11.1 turned out to include all essential properties needed to get an understanding of electro-mechanical oscillations in a power system and of the physical mechanisms that govern the dynamic process. A corresponding simple system for voltage stability is shown in Figure 13.1, which shows how a load is supplied via a line. It shall be emphasised that a detailed analysis of voltage stability in a power system is considerable more complicated than what shall be shown here, but the fundamental mechanisms can be illustrated by the simple system in Figure 13.1. If it is assumed that the line is lossless, i.e. Ze = jXe the relations below hold Ul UN Pl = sin (13.1) Xe Ql = Ul UN cos Ul 2 Xe (13.2)
The angle is not of interest and can be eliminated by moving the non cosine term to the left hand side in (13.2), thereafter squaring both sides of the new equation and of (13.1). After adding the new equations one gets: U2 Pl + Ql + l Xe
2 2
Ul UN Xe
(13.3)
136
Nose Curve 1.2
0.8 Load Voltage (U ) Pl,max

l
0.6
0.4
0.2
0.2
0.4
0.6
0.8
1 1.2 Load Power (P )

l
1.4
1.6
1.8
Figure 13.2. Relation between active power and voltage magnitude in the load node. The source voltage UN is chosen so that nominal voltage gives nominal load power. (UN = 1.044, Xe = 0.3 p.u., and power factor of the load = 1.)
It is now assumed that the reactive load varies in a given way with Pl and Ul , e.g. so that the load has constant power factor: Ql =k Pl (13.4)
If the relation between Pl and Ul is plotted in a diagram, Figure 13.2 is obtained, where it is assumed that k in eq. (13.4) is zero, i.e. the power factor = 1. The type of diagram shown in Figure 13.2 is often called a nose curve by obvious reasons, or a PU-curve. From this gure three observations can be made concerning the number of (theoretically) possible operating points: 1. For Pl < Pl,max there are two solutions 2. There is exactly one solution for Pl = Pl,max .
13.2. Simple Systems for Analysis of Voltage Stability 3. There is no solution for Pl > Pl,max
137
These observations will now briey be discussed. For values on Pl less than Pl,max it was concluded that there are two solutions. Physically the upper solution corresponds to high voltage and low current, while the lower one corresponds to low voltage and high current. In the power system the upper solution is of interest, since it gives lower losses. Furthermore in most situations this solution has, which will be shown in section 13.3, other properties which makes it attractive. For Pl = Pl,max there is exactly one solution. The value on Pl,max depends on Xe and UN , but also on the power factor for the load. The more reactive power the load draws the less will Pl,max be. For loads larger than Pl,max there is obviously no solution. 1 For loads larger than Pl,max there is obviously no solution, and it is consequently clear that the system is unstable if a load greater than Pl,max , is connected. The analysis made above was quite simplied. It was assumed that the load power Pl could be varied as an independent parameter. However, for a realistic physical load the load power will be a function of the load voltage Ul , i.e. Pl = Pl (Ul ). This function is often referred to as the load (voltage) characteristics. The behaviour of the function Pl (Ul ) depends on the time scale considered. In a shorter time scale it is the physical processes involved in the load device that determine the characteristics, while in a longer time scale various control systems must also be taken into account. It turns out that in a the shorter time scale many loads could be modelled as Pl = Pl0 Ul Ul0
kp
(13.5)
for not too large variations around the operation point Ul = Ul0 , with kp typically 1 2. In the longer time scale control systems strive for keeping the load power constant for many load devices. Load modelling is very important in voltage stability analysis, but the detailed knowledge about load models in real systems is often very rudimentary. Much work is being done in the eld, but still more work is needed. Another simple system that could exhibit voltage unstable behaviour is shown in Figure 13.3. In this system there are two fairly strong systems, 1 and 2, interconnected by long transmission lines. Between the two areas 1 and 2 there is a small system connected, denoted by index m in Figure 13.3.
1 Mathematically the point Pl = Pl,max is called a bifurcation point with the active load as parameter. At this point the two solutions coalesce to one, which will disappear when Pl is further increased. It can be worth mentioning that it is not always that the maximum available load power is a bifurcation point, but this depends on the independent parameter in the system. In the studied case Pl is the parameter, but sometimes it can be the load impedance or another parameter, e.g. numbers of connected asynchronous machines. In these cases maximum available load power and the bifurcation point do not necessarily coincide. See also the discussion in section 13.3
138
~
U1
Um
Figure 13.3. Simple system equivalent for studying voltage instability.
Um U1 U2
Figure 13.4. Voltage phasors for the system in Figure 13.3.
For simplicity it is assumed that X1 = X2 , U1 = U2 , and that no load is consumed at the intermediate bus. The voltage Um will then be Um = U1 cos(/2) (13.6)
where is the angle between the voltages U1 and U2 , see Figure 13.4. If now the angle is increased due to an increase in the transmitted power or in a fault increasing the reactance between areas 1 and 2, the voltage Um will decrease. If there is no reactive power resources at this point it might lead to a voltage collapse in this area, that could in worst case inuence the whole system.
13.3
Analysis of Voltage Stability
The simple system of Figure 13.1 will now be analysed for a few cases. In the analysis of voltage stability a few indicators have been introduced to simplify
~
U2
X1
X2
13.3. Analysis of Voltage Stability
139
the stability analysis. These indicators are based in small signal properties of the system and are consequently derived from linearised models. As concluded from the analysis of synchronous stability, linear analysis is often very powerful, but most be used with care since it only captures the local behaviour of the system around an operating point. To determine the global behaviour the full non-linear model must be employed.
13.3.1
Stability Indicators
In this section voltage stability indicators of the system in Figure 13.1 will be discussed. Similar indicators with similar properties exist also for larger systems and these are often used in voltage stability studies of real systems. Voltage Sensitivity Factor The Voltage Sensitivity Factor VSF is one of the most used indicators or indices of voltage stability. It is dened for a given node i as V SFi = Ui Qi (13.7)
The physical interpretation of V SF is that it measures the change in voltage magnitude at a given node as a consequence of a reactive power injection in that node. Since all voltage control in a power system is based on that a reactive power injection results in voltage increase, it is clear that a stability criterion is V SFi > 0 (13.8) Furthermore, it is reasonable to require that an injection of reactive power in any node will not decrease the voltage in any node in the system. This property is often called voltage regularity. The V SF -indicator will be demonstrated on the system in Figure 13.1 later on. U/E-indicator Another indicator that is more closely related to the voltage control in a power system is given by U/E, where U is the voltage magnitude of a given bus and E is a voltage that is controllable, e.g. at the terminals of a synchronous machine. Also in this case the requirement of stability is given by U >0 (13.9) E Qg /Ql-criterion This criterion relates the reactive power consumed at the load(s) to the reactive power generated by the synchronous machines. In Figure 13.1 it
140
thus relates the reactive power injected at the innite bus, here called Qg , to the reactive power consumed by the load Ql . The stability criterion is Qg >0 Ql (13.10)
It can be shown that all the stability criteria above are equivalent in most realistic systems. Which one to use is then dependent on which system variables that are known and controllable in a given situation. The three criteria have also that in common that at the transition point to instability they tend to innity before going negative.
13.3.2
Analysis of Simple System
To gain a better understanding of the mechanisms governing voltage stability a few dierent cases will be analysed. The system in Figure 13.1 will be used in all cases. Inuence of System Parameters Here the inuence of a few system parameters will be investigated. The nose curve in Figure 13.2 was drawn for a lossless system, i.e. = 90 , and the power factor of the load being equal to unity. In Figure 13.5 the power factor has been varied, and in Figure 13.6 nose curves for a number of dierent values of the impedance angle is plotted. In all the cases the magnitude of the source voltage UN has been adjusted so that when the load is 1 p.u. the load voltage Ul = 1 p.u. From Figure 13.5 it can be concluded that the power factor has a signicant inuence on the nose curve, and thereby on the voltage stability. If the reactive consumption of the load could be compensated at the load point, this results in higher Pl,max but also in a better voltage regulation around the nominal operating point. The conclusion drawn from Figure 13.6 is that the impedance angle does not inuence the nose curve signicantly for realistic values of high voltage systems. A X/R ratio of 10 results in 84 and the deviation form the curve cooresponding to = 90 is insignicant. Thus, in the following it will be assumed that = 90 . Inuence of Load Characteristics The voltage characteristics of the loads have a very large impact on the voltage stability as will be shown in the following. It should be pointed out the the actual dynamics occurring in a real system during the course of events to be discussed here could be very complicated and that the description given here is simplied to emphasise the fundamental dynamics in the system. The very detailed interplay between load dynamics and voltage control might result that the actual course of events turn out to be dierent.

Nose Curve 1.2
141
1 3 0.8 1 Load Voltage (Ul) 2
0.6
0.4
0.2
0.5
1 Load Power (Pl)
1.5
Figure 13.5. Nose curves for dierent values of the power factor of the load of the system in Figure 13.1. Curve 1: Power factor = 1. Curve 2: Power factor = 0.95 (inductive). Curve 3: Power factor = 0.95 (capacitive)
Load Admittance Inuence on Load Power In Figure 13.7 the load characteristics of an impedance load, i.e. with kp = 2 in eq. (13.5) has been plotted together with two dierent nose curves. A very important property of a load is the load power controllability, i.e. the possibility to increase or decrease a load power by controlling suitable parameters. If it assumed that the physical load in Figure 13.7 is a heater, this control is achieved by changing the resistance (conductance) of the load. Assuming that the voltage is constant the load power can then be increased by decreasing the resistance of the load, i.e. to increase the conductance of the load. In a more general case this is usually expressed as an increase in load power is achieved by increasing the load admittance. Increasing the load admittance in a general case, when the load is more complicated than a simple resistor as in the heater case, should be interpreted as connecting more load devices of the same kind at the load node. Actually this has
142
Nose Curve 1.2
0.8 Load Voltage (U ) 3 0.6 2

l
0.4
0.2
0.5
1 Load Power (P )
l
1.5
Figure 13.6. Nose curves for dierent values of the impedance angle of Ze of the system in Figure 13.1. Curve 1: = 90 . Curve 2: = 80 . Curve 3: = 70
formed the basis for one denition of voltage stability: A system is said to be voltage stable if an increase in load admittance results in an increase in load power. Let us now revert to Figure 13.7. If the nose curve of the system is according to curve 1 it is seen that an increase in load admittance, dashed load curve, gives an increase in load power. Correspondingly a decrease in load admittance, dotted curve, for the operating point on nose curve 1, gives a decrease in load power. Thus the system is voltage stable according to the denition above for the undisturbed (nominal system). If a disturbance occurs increasing the line impedance drastically, the nose curve will be according to curve 2 in Figure 13.7. The new operating point will the intersection between the load curve and the nose curve resulting a lower load power and a lower voltage. If there is a controller trying to restore the load power to nominal it will increase the load admittance (dashed load curve), which will result in a voltage decrease and a load power decrease. This could be explained in the following way. The load admittance increase results in an increase

Nose Curve and Load Characteristics 1.2
143
1 1 0.8 Load Voltage (U )

l
0.6
0.4 2
0.2
0.2
0.4
0.6
0.8
1 1.2 Load Power (P )

l
1.4
1.6
1.8
Figure 13.7. Nose curve for undisturbed system, curve 1, and for a system after a contingency resulting in increase of Ze , curve 2. Load characteristics for impedance load at nominal operating point, solid curve. Load curve with 10% increase in load admittance, dashed curve. Load curve with 10% reduction in load curve, dotted curve.
of the load current. This current increase is however so large that the load voltage drops so much that the load power decreases. In the nominal operating point of nose curve 1, the voltage drop due to the current increase is however so small that the net eect is an increase in the load power. The situation is clearly unstable for nose curve 2, and the system will collapse. In Figure 13.7 the voltage after the disturbance, nose curve 2, is very low and in a real system protections would have been activated. It should also be noted that since Pl,max of nose curve 2 is much lower than the nominal power, the load power could never recover to the pre-fault value. Generally it can be seen that an increase in load admittance will give an increase in load power if dPld dPl > (13.11) dU dU where Pld refers to the load device, i.e. the load characteristics, and Pl to the load power of the nose curve.
144 VSF Criterion
The V SF criterion was introduced above and that will now be analysed for dierent load characteristics. This criterion is closely related to how the voltage control will act in the system. It is rather straightforward to calculate V SF for dierent load characteristics and operating points on the nose curve. From such an analysis the following can be concluded: For constant power loads, i.e. kp = 0 in eq. (13.5), operating points on the upper half of the nose curve are stable, i.e. V SF > 0, while those on the lower part are unstable or V SF < 0. For impedance loads V SF > 0 for all operating points, indicating voltage regularity. For 0 < kp < 2 operating points on the upper part of the nose curve are stable. In addition a part of the lower part is also stable ranging from Pl,max to lower Pl values. This interval increases with kp with boundaries according to the two previous conclusions. Since in a longer time scale most loads could be approximated by constant power loads, only permanent operation at the upper part of the nose curve is feasible. During transients operation on the lower part could be acceptable. Load Dynamics The discussion in this paragraph relates to discussion concerning load admittance, but deal with a somewhat dierent disturbance. The system in Figure 13.1 is considered again. It is now assumed that a voltage increase of 5% occurs in the sending end system, i.e. UN is increased by 5%. In Figure 13.8 the nose curves and load characteristics corresponding to an impedance load is depicted. The load is controlled in such a way that the load power should be constant, e.g. thermostat controlled heater. The voltage increase resulted in a load power increase and consequently the load power control orders a decrease in load admittance. This decrease in load admittance results in the dashed load characteristics in Figure 13.8, and the system can settle to a new operating point at the desired load power, i.e. 1 p.u., but at a slightly higher voltage. An analysis of the behaviour of the system if the sending end voltage is slightly decreased shows that the system will settle at 1 p.u. load power at a slightly lower voltage. Now assume that the system is operating on the lower part of the nose curve, see Figure 13.9. According to the V SF criterion this operating point is stable. The same disturbance as in the previous case is considered. From Figures 13.8 and 13.9 it can be concluded that both systems are stable according to the U/E criterion since an increase in the sending voltage

145
1.2
1 2 Load Voltage (U )
l
0.8 1 0.6
0.4
0.2
0.5 Load Power (P )

l
1.5
Figure 13.8. Nose curve for undisturbed system, curve 1, and for a system where the sending end voltage has been increased by 5%, curve 2. Solid load characteristics corresponds to load admittance before voltage increase. Dashed load characteristics corresponds to decreased load admittance to give load power = 1 p.u. for system after disturbance.
causes an increase in the load voltage. It is seen that voltage increase in Figures 13.9 implies a load power increase. The same load power controller as in the previous case would order a decrease in the load admittance to decrease the load power. Such a load admittance decrease is represented by the dashed load characteristics in Figures 13.9. But, the load admittance decrease results in a load power increase. If there are no restrictions on the action of the controller it will continue to decrease the load admittance until the point A of Figure 13.9 is reached. A similar analysis shows that a voltage decrease in the sending end system will result in that the load admittance will increase until (theoretically) the point (0,0) is reached. It is thus clear that the initial operating point of Figures 13.9 is unstable.
146
1.2 A 1
2 1
Load Voltage (U )
0.8
0.6
0.4
0.2
0.5 Load Power (P )

l
1.5
Figure 13.9. Nose curve for undisturbed system, curve 1, and for a system where the sending end voltage has been increased by 5%, curve 2. Solid load characteristics corresponds to load admittance before voltage increase. Dashed load characteristics corresponds to decreased load admittance.
QU curves Finally another method sometimes used to analyse voltage stability will be demonstrated. In the system of Figure 13.1 one can calculate the amount of reactive power that must be provided to the network for a given load power for dierent load voltages. If this is done for the system of Figure 13.1 assuming Pl = 1 p.u., a curve according to Figure 13.10 is obtained. The same data as in Figure 13.2 have been used. Since Figure 13.2 was calculated with unity power factor of the load, the load voltages corresponding to Q = 0 of Figure 13.10 should correspond to the voltages giving Pl = 1 p.u. in Figure 13.2. This is easily veried. Furthermore, the V SF criterion for stability says that an reactive power injection should result in a higher voltage. In the QU curve this corresponds to a positive derivative, which are for load voltages higher than 0.6 p.u.
147
QU Curve
0.3 0.2 0.1 Reactive Power Demand of Network 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.3
0.4
0.5
0.6
0.7 Load Voltage (U )

l
0.8
0.9
1.1
Figure 13.10. QU curve for the system of Figure 13.1. (UN = 1.044, Xe = 0.3 p.u., Pl = 1.0 p.u.)
148
14
Control of Electric Power Systems
This chapter gives an introduction to control of power systems. Both the faster distributed controls and slower centralised controls are described.
POWER SYSTEM is planned, designed and built to supply the consumers with electrical energy considering:
Economy Quality Supply security (Reliability) Environmental impact The last point has during the latest years become more and more important. This can be seen from the fact that most of the larger companies nowadays have special departments responsible for that the companies follow laws and regulations within the environment area. The three rst points above constitute the basis of the optimisation that the power companies must make regarding their investments and daily operation. It is evident that economical considerations must be regarded when deciding issues about quality and supply security for the system. The higher quality and supply security that is required, the more expensive the electrical energy becomes. Methods based on mathematical analysis and decision making theory (operations research) have been worked out to handle this optimisation problem, which contains several considerations/adjustments that are hard to formulate in a stringent mathematical form. As cheap and safe access to electrical energy is of great importance for almost all activities in a modern society, it is not surprising that the decisions makers in the power industry were pioneers to utilise advanced optimisation methods and other analysis tools for expansion and investment planning. In this compendium we will limit ourselves to discuss how a power system can be controlled to fulll the constrains regarding quality. A very short introduction is also given about the overall control in a power system, which besides its aim to fulll the quality even plays an important role in the supply security. 149
150
14. Control of Electric Power Systems Three important factors that decide the electrical quality are: Frequency variations Voltage variations Waveform of voltage and current
The two rst factors will be dealt with in this course, but the last point is getting greater importance because more nonlinear components are connected to the electrical system, e.g. power electronics. These give rise to harmonics, see Chapter 14 in the Static Analysis. It has been shown earlier that the frequency is a very good indicator on the active power balance in a power system. The frequency is constant when the same amount of electrical power is produced as consumed by the loads, including system losses. If this is not the case frequency changes will occur. It can also be noted that the frequency is the same in the whole system at steady state. This depends on that the active power easily can be transported in the system (that is why it has been possible to build the power systems of todays size). It is often said that the active power is a global quantity, unlike the reactive power which is a local quantity, since reactive power cannot be transmitted over longer (electrical) distances. (This is because of that X normally is much greater than R in a power system, at least for transmission and sub-transmission networks.) The frequency of the system is reduced when a load increase is not compensated for by a corresponding increase of the turbine power of the connected generators. The power decit decelerates the generator rotors and consequently the frequency is reduced. Frequency reductions also arise when production is lost, e.g. as a consequence of failures in the system which lead to that protections disconnect the failed equipment. Too large reductions of the frequency could lead to system collapse, since a lot of equipment in the power stations, e.g. power supply systems, do not tolerate too low frequencies. A load reduction in the system which is not compensated for by a reduction of turbine power leads to a frequency increase. The permitted stationary frequency deviation in interconnected power systems is at normal operation typically 0.2 % or 0.1 Hz. There are also sometimes requirements on how large the dierence between actual, physical time and the time corresponding to the integrated frequency from the system can be. This time dierence is normally not permitted to be larger than 10s. It is also important that the voltage deviations in the system is limited. This is of importance for the connected loads, but a good voltage prole is also essential for keeping the losses low and for utilising the reactive reserves to establish a secure operation of the system. Voltage control is, as been pointed out earlier, a more local control than the frequency control. If the
14.1. Control of Active Power and Frequency
151
Figure 14.1. Typical variation of the load in a power system. The upper curve shows the highest hourly average values over the year and the lower curve the lowest hourly average value.
voltage deviate from the set value in a node, the control action must be made in this or a nearby node.
14.1
Control of Active Power and Frequency
Since the stored energy in the system is relatively small, the constant of inertia for a typical generator is 4 s, the electrical energy must be produced in the same moment as it is consumed by the loads. Since the load varies, a certain power reserve must constantly be availably. In addition for the daily variation, see Figure 14.1, there are continuously spontaneous load variations up to 2 % of the total load during a minute. The generation reserves are generally divided into dierent groups according to their properties: spinning, supplementary and back-up, see Figure 14.2. The reserves are at operation planning divided into normal operation and disturbance reserve after the cause of the needed reserve. At the operation planning of the power system forecasts over the expected load are carried out continuously. Forecasted values can never exactly
152
MW
Total Instaleed Capacity
14. Control of Electric Power Systems
Unavailable Generation Reserves Generation = Loads + Losses
Operation Planning Disturbance Reserve Normal Operation Reserve
Operation Back-up
Supplementary Spinning
Figure 14.2. The dierent generation reserves in a power system.
coincide with real values, why a reserve which can adjust the dierence is required. Therefore frequency variations arise, which must be compensated for. Both power imbalances due to incorrect load forecasts and to occasional load variations are controlled with the normal operation reserve. Events which can lead to utilisation of the disturbance reserve are: generator trips or line trips. Disturbances of that kind can lead to frequency reductions and reserves must often be put into operation. In the operation the reserves are divided after needed time for activation, see Figure 14.2. To keep reserves is expensive and therefore it is of interest to minimise the needed eort of the reserves for maintaining the wanted reliability and security.
14.1.1
Spinning reserve
The spinning reserve is co-ordinated in the Nordic system and is above all located in hydro units with turbine control, see Figure 14.3. The turbine control is activated within some seconds if the frequency deviates from the normal and changes then the turbine setting for the regulating generators by changing the guide vane opening to the turbines. When the system frequency is changed the power demand of certain loads is also changed, specially for motors, in such a way that a frequency increase leads to increased power consumption and a frequency decrease gives lower power consumption. This frequency dependence of the load stabilises the frequency. Components for heat- and light production are fairly insensitive to frequency variations. The resulting change in active power, P , at a frequency change in the system is: P = Pg Pl (14.1) where: Pg = the total change in output power from the generators participating in the frequency control, Pl = the total frequency dependent
14.1. Control of Active Power and Frequency
153
U f Ufield Turbine P mech P mech
Uset Ugen
~
f fset f
Power System
Figure 14.3. Generator unit with voltage regulator and turbine regulator for frequency control.
load power change. The frequency dependence of the load can often be approximated as Pl = Df (14.2) where D is a constant. The value on D is typically 0 4 %/Hz, but large variations occur depending on the kind of load. The frequency control is performed so that, at least for small frequency deviations, the total change in output power of the frequency controlling generators are (in steady state) proportional to the frequency change, f . The relation between changes in active power and frequency deviations can then be written (14.3) Pg = Rf The constant R is determined by the setting of the turbine control of the frequency controlling generators, and is named the systems frequency regulation constant or frequency droop and is normally given in MW/Hz. The frequency regulation for a system is established so that the frequency does not fall outside given limits at certain dimensional disturbances which can occur, e.g. trip of the largest generator in the system. It shall be noted that the variation of frequency, computed from the above relations, are the steady state values which are obtained after the frequency controllers have responded. The transient frequency deviation during the process following a large disturbance can be considerably larger than the stationary one, up to some Hz. In Figure 14.4 it is shown how the frequency can vary after a disturbance in the system. As seen the maximal frequency deviation is considerably greater than the stationary one. The transient behaviour of the frequency is highly dependent on the
154
Figure 14.4. The transient frequency response after a production loss. The diagram is a recording from the Nordic system at a loss of a 1000 MW unit on the 24 November 1983.
characteristics of the turbine regulators of the generators, and in this respect the hydro turbines and steam turbines are totally dierent. Shortly it can be said that the transient frequency deviation is considerably larger if the frequency control is performed by hydro power. This is due to that the hydro turbine is a non-minimum-phase system.
14.1.2
Supplementary Reserves
Supplementary reserves must also be available to cope with fast changes in the generation plan at large forecast deviations or at permanent outages when the spinning reserve has been used and must be restored. From the remaining deviation of frequency and the knowledge about the system frequency droop it is straightforward for the system operators to compute how much new power production that must be connected to restore the system within the prescribed limits. This can either be done by connecting new production or by increasing the turbine power at generators in operation. In some systems this is handled manually, but in some countries, e.g. in continental Europe and USA, this is at certain power companies accomplished by special slow automatic controllers, so-called Automatic Generation Control (AGC). This regulation is considerably slower than frequency control and must be designed not to interact with the frequency control in an adverse way.
14.2. Control of Reactive Power and Voltage
155
14.1.3
Back-Up Reserves
The slow reserve shall be availably within about 2 h. It is needed for security reasons to restore the fast reserve when this has been used. The slow reserve can consist of thermal power which is kept at stand by. This implies that pressure and temperature is kept on a level which permits synchronisation and load recovery within 2 hours. Co-generation back units aggregate for which it is possible to increase the electrical power can also be used.
14.2
14.2.1
Control of Reactive Power and Voltage

Reactive Power Control
In steady state operation both active power balance and reactive power balance must be maintained. The reactive power generated by synchronous machines and capacitances must be equal to the reactive power of the loads plus the reactive transmission losses. If the active power balance is not kept, the frequency in the system will be inuenced, while an imbalance in reactive power will result in that the voltages in the system dier from the desired ones. If the power system is operated in the correct way, the voltage drops on the lines are usually small. The voltages in the nodes of the system will then almost be the same (at voltage prole). In this case the transmission system is eectively used, i.e. primarily for transmission of active power, and not for transmission of reactive power. As known from the Static Analysis the voltage in a system is strongly effected by the reactive power ow. Consequently the voltage can be controlled to desired values, by control of the reactive power. Increased production of reactive power gives higher voltage nearby the production source, while an increased consumption of reactive power gives lower voltage. Therefore it is of great interest to study which components and devices which can be used to regulate the reactive power in a power system. While the active power is entirely produced in the generators of the system, there are several sources of reactive power. In the other hand the reactive power cannot be transported over long distances in the system, since normally X >> R in a power system. Important producers of reactive power are: Overexcited synchronous machines Capacitor banks The capacitance of overhead lines and cables Important consumers of reactive power are:
156 Inductive static loads
Under-excited synchronous machines Induction motors Shunt reactors The inductance of overhead lines and cables Transformer inductances Line commutated static converters For some of these the reactive power is easy to control, while for others it is practically impossible. The reactive power of the synchronous machines is easily controlled by means of the excitation. Switching of shunt capacitors and reactors can also control the reactive power. If thyristors are used to switch capacitors and/or thyristors are used to control the current through shunt reactors, a fast and step-less control of the reactive power can be obtained. Such a device is called SVC (Static Var Compensator). As has been shown earlier it is most eective to compensate the reactive power as close as possible to the reactive load. There are certain high voltage taris to encourage large consumers, e.g. industries, and electrical distributions companies to compensate their loads in an eective way. These taris are generally designed so that the reactive power is only allowed to reach a certain percentage of the active power. If this percentage is exceeded, the consumer has to pay for the reactive power. The high voltage network is in that way primarily used for transmission of active power. The reactive looses of power lines and transformers depend on the size of the reactance. At overhead-transmission line the reactance can be slightly reduced by the use of multiple conductors. The only possibility to radically reduce the total reactance of a transmission line s to connect a series capacitor, see 8.1.4 in the Static Analysis.
14.2.2
Voltage Control
The following factors inuence primarily the voltages in a power system: Terminal voltages of synchronous machines Impedances of lines Transmitted reactive power Turns ratio of transformers
14.2. Control of Reactive Power and Voltage
157
U1
zk
U 2
U2
N1
N2
Figure 14.5. Transformer with variable turns-ratio (tap changer).
A suitable use of these leads to the desired voltage prole. The generators are often operated at constant voltage, by using an automatic voltage regulator (AVR). The output from this controls via the electric eld exciter the excitation of the machine so that the voltage is equal to the set value, see Figure 14.3. The voltage drop caused by the generator transformer is sometimes compensated totally or partly for, and the voltage can consequently be kept constant on the high voltage side of the transformer. Synchronous compensators are installed for voltage control. These are synchronous machines without turbine or mechanical load, which can produce and consume reactive power by controlling the excitation. Nowadays new installations of synchronous compensators are very rare. The impact of the impedances of the lines on the reactive power balance, and thereby the voltage, have been analysed in the Static Analysis. These are generally not used for control of the reactive power. Series capacitors are generally installed to increase the active transmission capacity of a line. From the static analysis it is also known that the reactive power transmitted, has a great impact on the voltage prole. Large reactive transmissions cause large voltage drops, why these should be avoided. Instead, the production of reactive power should be as close as possible to the reactive loads. This can be made by the excitation of the synchronous machines, which have been described above. However, there are often no synchronous machines close to the load, so the most cost-eective way is to use shunt capacitors which are switched according to the load variations. An SVC can be economically motivated if fast response or accuracy in the regulation is required. Shunt reactors must sometimes be installed to limit the voltages to reasonable levels. In networks which contain a lot of cables this is also necessary, since the reactive generation from these is much larger than from overhead lines. (C is larger and X is smaller.) An important method for controlling the voltage in power systems is by changing the turns ratio of a transformer. Certain transformers are
158
equipped with a number of taps on one of the windings. Voltage control can be obtained by switching between these taps, see Figure 14.5. If switching during operation can be made by means of tap changers, this possibility to voltage control is very eective and useful. Normally the taps are placed on the high voltage winding, the upper side, since then the lowest current needs to be switched. If N1 is the number of turns on the upper side and N2 is the number of turns on the lower side, the turns ratio of transformer is dened as = N1 N2 (14.4)
Then the relation between the voltage on the high voltage side, U1 , and on the low voltage side, U2 , at no load is U2 = U1 (14.5)
If the voltage decreases on the high voltage side, the voltage on the lower side can be kept constant by decreasing , i.e. by switching out a number of windings on the high voltage side. When the transformer is loaded eq. (14.5) is of course not correct, since the load current gives a voltage drop over the leakage reactance of the transformer, zk , but the same principle can still be applied at voltage control. Transformers with automatic tap changer control are often used for voltage control in distribution networks. The voltage at the consumers can therefore be kept fairly constant even though voltage variations occur at the high voltage network. Time constants in these regulators are typically some ten seconds. Sometimes the turns ratio cannot be changed during operation, but just manually when the transformer is o load. In this case one can only change the voltage level in large but not control the voltage variations in the network.
14.3
Supervisory Control of Electric Power Systems
The frequency- and voltage control described above are performed by local controllers, but overall and central control is also needed to secure safe and economical operation of the power system. This control is performed from control centres which generally are hierarchically organised. Often there is a national centre on the top, which supervises the national transmission system and co-ordinate the frequency control. In certain cases this centre can also co-ordinate the operation of dierent power companies, which in that way can optimise their operation. Next level contains a number of regional control centres which co-ordinate and supervise the regional networks and
14.3. Supervisory Control of Electric Power Systems
159
Normal State Restoration System Collapse Extremis State Alert State
Emergency State
= Planned (Intended) Changes = Unplanned Changes (Disturbances)

Figure 14.6. The dierent operating states of a power system.
control the generation within that region. Next level of operation centre has the task of controlling the operation for one or several power stations. The information ow between these dierent levels is very large and there is an extensive communication system parallel with the electrical power system to handle this. Computers are extensively used to perform the overall control of the power system, but also the operators make important decisions in the control process. The dierent states that a power system can be found in is often described according to Figure 14.6. Normal operation is the state which is desired for a system. Besides that all quantities in the system are within the limits given by the equipment and the system, there are margins so that the system can tolerate certain predetermined disturbances, e.g. loss of lines and generation, without endangering the system security. When the system is in alert operation all quantities are still within the allowed limits, but the margins which existed for normal operation have been lost as a consequence of some disturbance(s). In alert operation power can be supplied as usual, and the consumers do not notice that any disturbance has aected the system. In this state the goal for the control is to bring the system back to normal operation through connecting equipment to establish margins to cope with new disturbances. Emergency operation can occur if further disturbances occur or if a large disturbance strikes the system. Some quantities, e.g. voltages or power ows on the lines, are now outside the permitted limits. Still there is no shortage of power, but fast actions must be taken otherwise further equipment will be tripped by dierent protections. If it is not succeeded to bring the system to alert operation and gradually
160
14. Control of Electric Power Systems Operating Margins OK No No No Rating Limits OK OK No No Pgen = Pload + Ploss OK OK OK No
Normal State Alert State Emergency State Extremis State
Table 14.1. Characteristics of the dierent operating states of a power system.
to normal operation, there is a risk to end up in extremis operation. In this state there is not enough with power to supply the connected loads, and the frequency is decreasing. Load shedding is usually needed to save the system. First of all such load is shed, that easily can be connected again without any larger damages to the consumers, for example electrical heaters or domestic loads. As a last resort loads, which demand long time for restoration, e.g. process industry, are shed. If this defence is not successful the system, or often parts of the system, will collapse, i.e. be de-energised, which is a state which must be avoided by all means. If this state should be reached the system must as fast as possible be brought back, often part by part, to normal operation. This process is called system restoration. The transitions which are marked with unbroken lines in the Figure 14.6 are initiated and performed from dierent control centres. (The broken lines denote undesired disturbances.) The central control plays an important role to get the system into normal operation. In Table 14.1 the characteristics of the dierent operating states of a power system are summarized.
Appendix A
Phase-Shifting Transformers
In three-phase transformers it is possible to arrange the windings in such way that not only the voltage magnitudes between the primary and secondary sides are dierent but also the phase angles. In this appendix a number of windings arrangements are described that could accomplish such a phase shift.
OR SINGLE-PHASE TRANSFORMERS the relation between the primary and secondary side voltages is a real number, known as the turnsratio, during sinusoidal steady state conditions. In three phase systems the windings can be arranged in dierent ways, and in some of these congurations a phase-shift is also introduced between the primary and secondary side. Consider the two three-phase transformers in Figure A.1. The transformer at the top has both the primary and secondary windings connected in a Y, called a Y-Y or wye-wye connection and denoted Yy. Obviously there is only a transformation of the magnitude of the voltages between the two sides. In the transformer at the bottom however, the primary side is delta connected and the secondary side is connected in Y, often denoted Dy or y. In addition a phase shift is here introduced between the primary and secondary windings, in this case 30 . Obviously other multiples of 30 phase shifts can also be obtained. An obvious way to incorporate the phase-shift in the transformer model would be to allow for a complex turns-ratio, i.e. t = a ej , as explained in subsection 2.2.2. As also shown there, a phase shift of a transformer has an inuence on the power ow, predominantly the active power ow, if one or more parallel paths exist. It is also shown that the magnitude of the turns-ratio, i.e. a in the expression above, inuences mostly the reactive power ow. A phase-shifting transformer would thus be a tool to control the active power ows on parallel paths. However, the transformers introducing multiples of 30 phase shift would in most cases be to crude for a meaningful power ow control. In realistic systems phase-shifts up to maximum 10 20 or so would be very powerful. In addition it would be desirable to be able to control the phase-shift depending on the loading conditions in the network. Thus a transformer according to Figure A.2 where the Ur etc can be controlled, both in mag161
162
R
S T
A. Phase-Shifting Transformers
r s t
S t
n R
S T
r s t
s
T S
n
Figure A.1. Two dierent congurations of three-winding transformer.
nitude and in phase, would be of value. This is the so called regulating transformer. The regulating function of these transformers is achieved through a tap changer that can change the turns on one or more windings and thereby changing the turns-ratio. Depending on the design both the magnitude and the phase of the turns-ratio could be changed. It is obvious that a tap-changer in the transformer connection in Figure A.1 would only change the magnitude of the turns-ratio. (It is assumed that the tap-changers on the three phases are identical and operated identically.) Such arrangements could thus only be used for voltage control. This is one of the main methods to control the voltage in the distribution systems. In order to introduce a controllable phase-shift to control the active power ow more complicated winding congurations must be introduced. There are several ways of achieving this, and the solution chosen depends on a number of parameters, such as phase-shift range to be controlled, voltage dierence between the primary and secondary sides, power rating, etc. Phase-shifting transformers, or phase-shifters, are quite rare in comparison to non phase-shifting transformers 1 . However, phase-shifters play an im1
With non phase-shifting transformers are here meant transformers where a possible
163
DUr Ur Us Ut Ur+DUr=UR Us+DUs=US Ut+DUt=UT UT DUt Ut US
Figure A.2. Basic principle of a regulating transformer.
Ur
UR
Us DUs
portant role particularly in highly meshed systems where power transfers over long distances take place. Phase-shifters can then re-direct the power ows to circuits that are not so highly loaded. Another function can be to direct the power ow to high voltage lines from highly meshed networks at lower voltage levels. In Europe and North America several phase-shifters are installed for power ow control. Since the traditional tap-changers are mechanical devices, with a typical minimum time interval between subsequent switchings of several tens of seconds, these controllers can only be used for steady state power ow control. Tap-changers with thyristor switches are today available, and with these much faster switchings can be achieved. Thereby also dynamic power swings could be damped. For this purpose also completely power electronics based controllable devices are developed, so called FACTS devices (FACTS = Flexible AC Transmission Systems), which normally oer a much higher degree of controllability. To obtain a phase-shift, a voltage in quadrature to the phase voltage must be added. Therefore, phase-shifters are also called quadrature transformers or quadrature boosters. Figures A.3 and A.4 show two examples of how phase-shifters can be designed. In Figure A.3 the voltage U is inserted by a series transformer and the phase-shift is always accompanied by a certain change in magnitude. For the transformer in Figure A.4 an almost pure phase-shift can be obtained for small phase-shifts.
phase-shift does not inuence the power ow. Transformers connected Yd are commonly used as generator step-up transformers and as transformers feeding distribution networks. However, since no parallel paths exists, or the parallel transformers are identical, the phase-shift does not play any role for the power ow.
164
A. Phase-Shifting Transformers
Ur Us Ut
Ur+DUr Us+DUs Ut+DUt Booster DUt Ut
DUr Ur a Us DUs
Figure A.3. Phase-shift obtained by a series boost transformer.
Ur
DUr Ur+DUr aU r Ut Us DUs DUt
Us Ut
Us+DUs Ut+DUt
Figure A.4. Phase-shift obtained by a quadrature voltage.
Appendix B
Protections in Electric Power Systems
In this Appendix a brief summary of how protections are designed and how they function is given. The very important distance protections and their operating principles are discussed. Some special protections and system wide protections that are of relevance for power system stability is briey reviewed.
IFFERENT TYPES of protections are installed to protect the equipment in an electric power system. Their task is to disconnect failed or overloaded equipment or parts of the system to avoid unnecessary damages on equipment and personnel. The purpose is also to limit the impact of failures on the parts of the system that have not failed. Special types of protection are the system protections. Their task is to prevent collapse (black out) of the system or parts of the system. An intensive development of protections based on modern information technology is going on both regarding hardware and software. On the hardware side microprocessors have been used over a long time to implement dierent functions in the protections, and with the recent developments more and more complicated functions can be implemented in a reliable way. Powerful methods like signal processing, state estimation, and articial intelligence, are being integrated into the protections. In general the functions which earlier were handled with separate relays are increasingly being integrated with other functional units for control and supervision. Furthermore, more complicated criteria for activation of protections can be applied. The interested reader is referred to the literature for further information. The summary here is concentrated on general principles for protections.
B.1
Design of Protections
A protection for an electric power system comprises the following parts: Measurement device with current- and/or voltage transformers and other sensors measuring the relevant quantities. Relay which when certain conditions are fullled sends signals to a circuit breaker or another switching device. This relay was earlier a 165
166
B. Protections in Electric Power Systems separate unit, but can in modern protections be a part of a larger unit for protection, supervision and control. Circuit breakers which execute the given instruction(s) from the relay. Telecommunication system is mainly used at distance (line) protections to get a faster and more reliable performance. Power supply systems which shall secure the power supply to the protection system, even with faults in the system.
The requirements on a protection system are that they should be dependable, secure, selective, sensitive, and fast. Dependability means that the protection should react and do its action when a fault occurs for which it is designed to react for. To achieve desired dependability double or even triple sets of certain parts of the protection or of signal paths might be needed. Malfunctions can be divided into not occurring or unwanted operations. Normally none occurring operations are more serious malfunctions than unwanted ones. Security means that the protection should not react when no fault occurs or when a fault for which its not intended to react occurs. Selectivity implies that not more than necessary pieces of equipment and apparatuses are disconnected to isolate a fault. Sensitivity is needed to detect failures which cause small fault currents, e.g. high impedance faults. This implies that the risk for misoperations increases at small disturbances, e.g. at energisation of transformers, or at high load operation but normal operation. The protection should react fast to secure that damages on persons and equipment are prevented or limited. The protections are often classied according to the object that they protect. An example is shown in Figure B.1. If a failure occurs within an indicated area in Figure B.1 this area should be isolated from the rest of the network. Many of the protections which protect separate pieces of equipment or parts of a system with occupy a limited physical area are so called current dierential protections. These protections measure the dierence between two currents, which in normal operation should be equal, and the protection is activated if this deviation exceeds a predetermined value. Both dierences in amplitude and phase can trigger the relay. The principle for a current dierential protection is shown in Figure B.2.
B.2. Distance Protections
167
Generator Protection Distance Protection
Transformer Protection
Bus Bar Protection
Figure B.1. The dierent protection zones in a power systems.
B.2
B.2.1
Distance Protections
General Principles
So called distance protections are important protections concerning stability and dynamics in a power system. Their task is to disconnect faulted lines or cables. Since large parts of the power system consist physically of lines and these are exposed to dierent disturbances, e.g. lightning strokes, down falling trees etc., it is important that those faults can be isolated to minimise the impact on the rest of the system. The most common faults are ground (earth) faults, i.e. short circuits between two or more phases and ground (shunt faults). Also interruptions in the lines can occur (series faults). The operating principle of the distance protection is shown in Figure B.3. Current and voltage are measured in both ends of the line and from these an apparent impedance can be calculated: Z = U/I. In normal operation this impedance varies within a certain area (large and almost resistive values on Z), but if a fault occurs it will drastically change. The given value depends on where on the line the failure occurs, and from system parameters as line data and short circuit capacity, it can be calculated where the fault has occurred. For each distance protection there are several protection zones dened in the Z plane according to Figure B.4. A low value on Z implies that the fault is close to the measurement. From line data and short circuit capacity
168
B. Protections in Electric Power Systems
Protected Object
IA
IB
Relay Logics
Signal(s) to Breaker(s)
Figure B.2. Principles of a current dierential protection.
Instrument Transformers Breaker #2
Breaker #1
U, I Relay Logic Telecommunication
U, I Relay Logic
Figure B.3. The operating principle of a distance protection.
it can then be decided if the fault is in the protected line, within Zone 1, or not. If that is the case, a trip order is given to the breaker at the same station within some milliseconds, typically 10 ms, after Z has reached Zone 1. At the same time a trip order is given to the breaker in the other end of the line. This latter trip order is not needed for isolation of the fault, if the protection system in the other end works as it should, but this trip order (transfer trip) increases the security in the system. If the measured value on Z is in Zone 2 or 3, it implies that the fault is outside the actual line. This implies that neither breaker 1 nor 2 in Figure B.3 shall be opened. If the breakers, which according to the protection plane should isolate the fault, are not operated by some reason, other breakers which are further away from the fault must isolate it. These secondary
B.2. Distance Protections
169
Zone 3 Zone 2 Zone 1 Normal Value of Z
Figure B.4. Dierent zones in a distance protection.
breakers will be used rst after it is clear that the primary breakers have not isolated the fault. Therefore if, Z is in Zone 2, the breaker does not get the trip order until typically some hundred milliseconds have passed. To coordinate and tune the settings of the protections to give a fast, reliable, sensitive and selective protection system is a complicated and an important task in an electric power system. In modern protection systems dierent areas can be dened according to Figure B.4 with in principal arbitrary geometric shapes, which facilitates the work. A plan comprising the dierent areas of protections and time settings is usually called a selectivity plan. The work to establish a selectivity plan is often very time consuming because it should be appropriate for every feasible state of operation, i.e. for dierent numbers of generators and lines connected and also at dierent load levels. Often trade-os must be made to reach acceptable results.
B.2.2
Automatic Re-Closure
When a ground fault or a short circuit occurs an ionised plasma (arc) that carries the fault current is often formed. This arc remains ionised as long as a the current ows through it. If the fault current is extinguished it is usually sucient for the plasma to cool down during some hundred milliseconds to rebuild the isolation so that the line can be re-connected. This is used in many systems and the line is automatically re-energised after a given time period after fault clearing. This is the case in many systems, where failed 400 kV lines, or lines at higher voltage levels, are automatically re-closed after 400 ms disconnection. If the failure would still remain, the line is kept disconnected during a larger time period, typically 800 ms, before the next re-closure is done. If the failure would remain after the second re-closure
170
B. Protections in Electric Power Systems
attempt, more attempts to re-close will not be made. In this latter case the isolation of the line has probably been permanently damaged and it must be repaired before the line can be put in operation again. (This scenario is typical. Due to specic conditions in dierent countries, deviations can occur.) Automatic re-closure must be made with certain carefulness when it is made close to large thermal power plants, e.g. nuclear power plants. The connection can in this case cause large stresses on the generator shaft if it occurs at certain phase positions. The re-closure can in this case be done in the remote end of the line seen from the thermal power plant, and then a synchronised re-closure at the other end is made. (With synchronised closure is meant that the voltage over the breaker is zero when it is closed.) In that way the transients in the system are drastically reduced, but of course the re-closure takes longer time. In some systems all fault clearings on the high voltage grid are n made on all three phases. I other systems one phase clearing is used. This means that only the faulty phase(s) is (are) disconnected at fault clearing and re-closure.
B.3
Out of Step Protections
A synchronous machine which has fallen out of step, i.e. its angular velocity does not coincide with the angular velocity of the net, has lost the synchronism with the system, and the machine must be disconnected. In order to supply electrical power to the system it must be phased in to the system later on. During the time period when the synchronous machine falls out of step large current pulses will pass through the generator, and if these becomes too many they can damage the generator and limit its life time. Furthermore, vibrations that can jeopardise the generator can arise. To protect a synchronous machine which has fallen out of step the synchronous machine is equipped with a out of step rely. Also in this case an impedance is dened from voltage and current. If the operation state of the generator is in the critical impedance zone a clock is started. If the generator then comes into the critical zone repeatedly times, this is a criterion that the generator is falling out of step, and the protection gives a trip order to the generator breakers. Generators are also equipped with other protections against overload (stator current protection), over-excitation, etc.
B.4
System Protections
System protections are special types of protection, the primary task of which is not to isolate failed equipment, but to prevent that the total system or large parts of it collapse. System protections often use information from
B.4. System Protections
171
several dierent points in the system or quantities which can give a reliable diagnosis of the state of the system. These systems often works in a time scale which is considerably longer than the more device oriented protections which were considered earlier, typically several seconds or minutes. An example on a system protection is load shedding. This is used to avoid that the frequency in the system decreases under acceptable values if the generation capacity has dropped in the system. The load shedding then disconnects predetermined loads depending on how much and how fast the frequency is falling. Voltage collapse protection is another system protection, the task of which is to prevent voltage collapse in the system. Load shedding uses only the frequency as input signal, while the voltage collapse protection often uses several dierent quantities as input signals.

Modelling and Analysis of Electric Power Systems

Uploaded by

Copyright:

Available Formats

Modelling and Analysis of Electric Power Systems

Uploaded by

Document Information

Original Description:

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

Modelling and Analysis of Electric Power Systems

Uploaded by

Copyright:

Available Formats

What are some of the main components and models used in power flow analysis?

What are some of the main components and models used in power flow analysis?

What are two methods described for solving the power flow problem?

What are two methods described for solving the power flow problem?

Modelling and Analysis of Electric Power Systems

ETH Z rich u March 2006

4 Nodal Formulation of the Network Equations

Contents Problem Solvability . . . . . . . . . . . . . . . . . . . . . . . . 34 37 37 39 40 41 42 44 45 48 48 52 55 57 61 63 66 66 66 67 68 71

Power System Dynamics and Stability

v . 97 . 98 . 100 . 100 . 100

10 The 10.1 10.2 10.3

Contents 165 . 165 . 167 . 167 . 169 . 170 . 170

Power Flow Analysis

1.2. Fault Current Analysis

Fault Current Analysis

i+di Cdx u+du

Lines and Cables

2.1. Lines and Cables

with gkm = and

where the complex voltages are

Uk e jqk U jqk k ke k Ikm Ikm Uk e jqk U jqk k ke k Ikm

Um e jqm Um e jqm m Imk Imk m

Um e jqm Um e jqm m Imk m

Figure 2.4. In-phase transformer model

Figure 2.5. Equivalent model for in-phase transformer

1:akme jjkm p qp = qk + jkm Up = akmUk

Figure 2.6. Phase-shifting transformer with tkm = akm ejkm .

Figure 2.7. System for Example 2.6.

Unied Branch Model

Uk e jqk U jq k ke k k Ikm Ikm 1:tkm 1:tkm

Uq e jqk Uq e jqk q q tmk:1 tmk:1

Um e jqm Um e jqm m Imk Imk m

Figure 2.8. Transformer symmetrical model.

Uk e jqk Uk e jqk 1:tkm k 1:tkm Ikm k Ikm

Uq e jqq zkm zkm

t :1 Uqq jqq mk e tmk:1 q

Um e jqm Um e jqm m Imk Imk m

Figure 2.9. Unied branch model extended (-model).

Figure 2.10. A shunt connected to bus k.

load load Pk = Pk (Uk )

Qload = Qload (Uk ) k k

Figure 2.11. Model of a load connected to bus k. gen Ik k I gen k

Figure 2.12. Model of a generator connected to bus k.

Stator Current Heating Limit

Stator current heating limit

End region heating limit

Field Current Heating Limit

Stator End Region Heating Limit

2 Qkm = Uk (bkm + bsh ) + Uk Um bkm cos km Uk Um gkm sin km km

where the notation km = k m is introduced.

3. Active and Reactive Power Flows

Uk Um bkm cos km + Uk Um gkm sin km

3. Active and Reactive Power Flows

Phase-Shifting Transformer with akm = 1

2 Qkm = Uk bkm + Uk Um bkm cos(km + km )

3.4. Unied Power Flow Equations

Unied Power Flow Equations

3. Active and Reactive Power Flows

Ikm ; for k = 1, . . . , N mk Ik is the net current injection