Calculation of Power System Critical Loading ... - Semantic Scholar

Calculation of Power System Critical Loading Conditions Ian A. Hiskens

Yuri V. Makarov

Department of Electrical and Computer Engineering The University of Newcastle, Callaghan, NSW, 2308, Australia

Summary The paper considers the problem of nding points of maximum loadability which are closest (in a local

sense) to the power system operating point. This optimization problem leads to a set of equations which describe such critical points. Not all solutions of this set of equations are critical points. The paper therefore explores the nature and characteristics of solutions. A two stage algorithm is proposed for solving the critical point problem. The rst stage is simply to nd a point of maximum loadability which lies is a speci ed direction. The second stage uses a continuation method to move from that initial point to the desired critical point. The proposed algorithm is tested on an eight bus power system example.

1

INTRODUCTION

The usual approach to the assessment of power system security, in a quasi-static sense, is to assume some particular loading pattern, then determine the corresponding point of maximum loadability [1, 2, 3]. Nose curve ideas, i.e., P-V and/or Q-V curves, follow from this approach. However such a method of assessment can result in optimistic predictions of security. The system may be much more sensitive to load changes for some dierent loading pattern. This paper presents a method for determining the most critical loading condition, and the associated security margin. The power ow is central to the calculation of quasistatic security margins. It describes the steady state relationships between system voltages and parameters such as loads, voltage setpoints and network impedances. That description takes the form of a set of n nonlinear algebraic equations y0 + f(x) = 0 (1) where y0 2 Rn is the vector of speci ed independent parameters such as active and reactiven powers of loads and generators or xed voltages, x 2 R is the state, consisting of nodal voltages. The vector function f(x) de nes the sum of power ows or currents into each bus from the rest of the network. The power ow equations f(x) generally de ne a mapping from state space to a subset of parameter space, i.e., ?f : Rn ! L where L  Rn. The region L de nes the set of parameters for which power ow solutions exist. As parameters y vary, the solutions of (1) will also move in state space. Parameters y may move to a point where two solutions coalesce, with further variation of y resulting in the disappearance of that solution. Behaviour of that form is referred to as a saddle node bifurcation. It follows from the Implicit Function Theorem that at such bifurcation points det Dx f = det J(x) = 0 (2) i.e., J(x) is the Jacobian matrix of f(x). Further, we see that the region L must be bounded by surfaces of points satisfying (2). We shall de ne the boundary as  = f(x; y) : x; y 2 Rn; y + f(x) = 0; det J j(x;y) = 0g (3) The projection of  onto state space and parameter space will be referred to as x ; y respectively.

To ensure adequate security of a power system, it is important that the operating point be prevented from moving too close to . That is, operating points should always be (at least) some speci ed distance from . As parameters y correspond to physical quantities that can be measured and controlled, it is useful to consider this distance in terms of parameter space, i.e., d(y) = ky ? y0 k (4) where y0 is the operating point, and y 2 y , i.e., y is a point on the solution boundary. The shortest (or critical) distance miny2y d(y) gives a measure of power system security in the most dangerous direction of loading. In addition, the critical vector (y0 ? y) de nes the optimal way of controlling the power system to maximise security. Its largest components indicate parameters which contribute most to the security conditions [4, 5]. In this paper we address the issue of robustly nding the minimum distance to . This question has been investigated before, see for example [4, 5]. We are proposing a continuation approach to nding the points on  which are closest (in a local sense) to the operating point. We shall refer to these points as critical points. The paper is organised as follows. Section 2 establishes the mathematical description of critical points. Properties of the solutions of the critical point problem are discussed in Section 3. Section 4 proposes an algorithm for nding the critical points. An eight bus example is considered in Section 5. Conclusions are given in Section 6.

2

FORMULATION OF THE PROBLEM

One of the rst things to consider in the formulation of the optimization problem (5) min d(y) y 2y is that not all parameters are free to vary. Some parameters, such as power injected at buses which have no load or generation, must always be xed. We therefore adopt the following power ow formulation. Let m equations in (1) contain xed values of parameters y2 = y20 = const. The other (n ? m) parameters y1 are free to vary. Then the system (1) can be rewritten as y1 + f1 (x) = 0 (6) 0 y2 + f2 (x) = 0 (7)

Considering the revised system description (6),(7), it is shown in [6] that solutions of (5) are a subset of the solutions of the constrained optimization problem ext ky10 + f1 (x)k2 (8) x y20 + f2 (x) = 0 (9) where `ext' denotes extrema of the cost function (8). The optimization is subject to nonlinear constraints (9). The cost function (8) de nes the square of the distance between the points y1 and y10 , with both points belonging to the constraint hyperplane y2 = y20 = const. Using a Lagrangian multiplier approach, it is further shown in [6] that solutions of (8),(9) are given by the set of equations

?s + y10 + f1(x) = 0

y20 + f2(x) J1t(x)s + J2t(x)

= 0 = 0 The last equation in (10) can be rewritten as

(10)

J t (x)s = 0 (11) where J t = [J1t J2t ] and s = [st t ]t. If s 6= 0, the Jacobian matrix J(x) is singular, and the vector s is a left eigenvector corresponding to a zero eigenvalue. Therefore, considering the optimization problem (8),(9), and the condition (11) for s 6= 0, we can conclude that critical points, i.e., points on  that are minimal distance (locally) from the operating point y0 , satisfy the system (10).

Therefore, because y2 = y20 , the compontent ?s of the vector ?s is the distance vector y ? y0 . Further, it is known that the vector s is normal to the singular hypersurface y [5]. De ne the intersection of y and the y2 = y20 hyperplane as y1 ;y20 . It follows that s is normal to the hypersurface y1 ;y20 . So the vector y ? y0 is orthogonal to y1 ;y20 . This con rms that nontrivial solutions of (10) correspond to local minima, maxima or saddles of the distance from the point y0 to y1 ;y20 . The minimum over the distances associated with all the nontrivial solution points characterizes the global \level" of power system security.

3.3

A Graphical Illustration

A graphical illustration of a nontrivial solution point of (10) is given by Figure 1. The solution point must lie somewhere on the intersection of the singular margin y and the constraint hyperplane y2 = y20 = const. The left eigenvector s at the nontrivial solution point y is perpendicular to the singular boundary, and its component s coincides with the vector from the singular point y to the operating point y0 . The component  (the Lagrange multiplier vector) of s is orthogonal to the constraint hyperplane. y2 = const det J(x)= 0

-sλ -s

-λ

3 SOLUTIONS OF THE CRITICAL POINT PROBLEM

yo2

The system (10) can be considered as an extended power ow problem where the usual power ow equations (the rst two equations in (10)) are supplemented by the singularity condition (11). State space, i.e., the space of unknown variables, is now extended to include the additional variables s . So the variables are (x; s) 2 R2n.

3.1

Trivial and Nontrivial Solutions

There are two kinds of solutions to the critical point problem (10):  Trivial solutions corresponding to the condition s = 0. Those solutions are actually the solutions of the usual power ow problem (1). They are global minima (zeros) of the distance function (4). All trivial solutions coincide in parameter space.  Nontrivial solutions conforming to the condition s 6= 0. Those solutions belong to the singular margin  given by (3). Solution of (10) can result in either trivial or nontrivial solutions, depending on initial estimates of the variables and the numerical solution technique used for solving the problem. A technique which produces nontrivial solutions is proposed in Section 4.

3.2

Distance and the Left Eigenvector

Consider nontrivial solutions of (10). It can be seen from the rst equation of (10), and (6) that at critical points, y1 ? y10 = y1 ? s + f1 (x) = ?s (12)

y

yo

y2

Figure 1: Graphical illustration of a nontrivial solution point.

3.4

A 3 Bus Example

Consider the 3 bus system of Figure 2. The values P1 , P2 are considered as free parameters y1 in (6). The voltage magnitudes V1 = 1pu and V2 = 1pu are taken as xed parameters y20 in (7). The optimization problem (8) therefore transforms to (13) ext [(P1 ? 0:4)2 + (P2 + 0:6)2] 1 ;2 Figure 3 shows the plane of state variables 1 ; 2, and Figure 4 the plane of free parameters P1 ; P2, with parameters V1 ; V2 xed. Singular margins x (dashed lines) and contours of the cost function (13) are plotted on the plane 1 ; 2 in Figure 3. Solutions of the optimization problem (13) are also shown. Points A1,A2 are minima. They correspond to trivial solutions of the critical point problem (10). A1 is the `normal' operating point. Points marked B and C are nontrivial solutions of (10), with B1


B4 being saddle points of (13), and C1,C2 maxima. All nontrivial solutions of (10) lie on the singular margin x .

P1 = 0.4

1

2

0.3+j1

V 1= 1 δ 1 = var

All solutions of (13) except A2 are shown in Figure 4. (A2 coincides with A1, so is not marked.) This gure also shows vectors from the operating point A1 to the solutions of (10). Each of these vectors is normal to the singular margin y . (The apparent absence of orthogonality of the vectors with respect to the singular margin in Figure 4 is caused by a dierence in the horizontal and vertical scales of the gure.) Figure 3 shows clearly that B1


B4 are saddle points of the optimization problem (13). From Figure 4, we see that they satisfy (5) locally. On the other hand, C1,C2 are local maxima of (13). They are also points that locally satisfy maxy2y d(y).

P2 = -0.6 V 2= 1 δ 2 = var

0.2+j1

0.6+j2

V 3= 1 δ 3= 0

3



Depending on initial guesses of variables, and the solution technique, any of the points identi ed in Figures 3 and 4 could be obtained as a solution of (10). As seen in the example though, only some of those points are of interest to us. Thus, the problem is to nd saddle points on the boundary of the security region. This problem is quite dierent to usual optimization problems where minima or maxima are desired. An appropriate technique is discussed in Section 4.

Figure 2: Simple 3 bus power system. 3 B3

2

B4

C1 B2

Delta 2, rad

1

0

4 AN ALGORITHM FOR FINDING CRITICAL POINTS

A1 B1

-1

C2

A2

-2

-3 -3

-2

-1

0 Delta 1, rad

1

2

3

4.1

Figure 3: The 1 , 2 plane for the 3 bus example. B1 C1

A1

-1

-1.5

P2

B2

B4

-2 B3 -2.5

C2

-3

-3.5 -2

-1.5

-1

-0.5

0

0.5

1

1.5

2

Stage 1: Obtaining a Good Initial Estimate

The rst stage of the critical point algorithm must produce a good estimate of state variables x; s;  in the vicinity of the critical point. To achieve this, it is necessary to have some idea of the direction in parameter space from the operating point to the desired critical point. In practice this requirement does not restrict the usefulness of the method, as power system operators and planners will usually have a good idea of the way in which parameters of their system, such as loads, vary. Let the estimated loading direction be
y1 . Recall y2 = y20 . The initial estimate of the critical point can be taken as the point on the solution boundary  in the direction
y1 from the operating point. That point is given by

0

-0.5

The set of equations (10) that describe critical points also have solutions that are not of interest, e.g., trivial solutions where s = 0. Therefore, an algorithm for nding critical points must consist of two parts, (1) a way of obtaining a good estimate of the unknown state variables x; s;  in the vicinity of the critical point, and (2) a numerical technique that will converge reliably from that initial estimate to the critical point. Such an algorithm is described in this section.

2.5

P1

Figure 4: The P1, P2 plane for the 3 bus example. Figure 4 shows the sections of the singular margin y plotted in the plane of free parameters P1 ; P2. The solution space L has a number of `layers', with each layer restricted by a section of the singular margin. These layers re ect the non-uniqueness of solutions of the power

ow equations. Each layer eectively generates a pair of power ow solutions. For example, there is a single layer at the point A1. Consequently, there are two distinct solutions of the power ow problem, A1,A2. They are shown in Figure 3.


y1 + y10 + f1 (x) = 0 (14) 0 y2 + f2 (x) = 0 (15) t J (x)s = 0 (16) st s = 1 (17) where is the loading parameter in the speci ed direction
y1 , and k
y1k = 1. An alternative formulation of (16),(17) uses the right eigenvector to achieve the singularity condition, rather than the left eigenvector s . Many techniques have been proposed for solving this problem, for example [1, 2, 3]. In some cases direct methods have been used, whilst others have applied continuation methods [7] to obtain the solution.

4.2

Stage 2: Motion Along the Singular Margin

Stage one of the algorithm provided us with a point on  in the vicinity of the desired critical point. Let that

point be x ;  ; s = [st t]t . We now wish to move from that point to the critical point. Consider the equations ( 
y1 + s ) ? s + y10 + f1(x) = 0 y20 + f2(x) = 0 J t(x)s = 0

(18) (19) (20)

The initial point x ;  ; s is a solution of (18)-(20) when  = 1. But when  = 0, the problem is exactly that of (10). So the critical point is a solution when  = 0. Therefore, as  is varied from 1 to 0, the solution of (18)-(20) is distorted from the initial point given by stage one, to the critical point. Notice that because of (20), all points along that path lie on . Also, (19) ensures that the path lies on the y2 = y20 hyperplane. In solving this continuation problem, it is helpful to scale s so that ksk =  at the initial point  = 1. This scaled s will still satisfy (20). It is shown in [6] that as  is varied from 1 to 0, the distance ky1 ? y10 k always reduces. This is an important property as it ensures that the stage two algorithm will never converge to maxima of (10). Further, the algorithm will always nd a loading condition which is more critical, i.e., which corresponds to a lower security margin, than the direction speci ed at stage one. Many numerical techniques exist for solving this continuation problem [3, 6, 7]. The solution method of [6] was used for the example given in Section 5.

5

AN 8 BUS EXAMPLE

The algorithm of Section 4 was tested on the eight bus example shown in Figure 5. Operating point values of generation and load parameters are given in Table 1. The system contains three generators with xed terminal voltage, and four nonzero loads. Bus 6 is the slack bus. The nominal voltage of all buses is 220kV. Line parameters and the operating point voltage pro le are given in [6]. 5 G G1

L

4

5

1

L

L5

1

7

G

4

8

3 L3

L7

6

4

2 G6 G

2

Figure 5: Eight bus test power system. Table 1 Bus parameters for the 8-bus system

Bus no. p 1 2 3 4 5 6 7 8

Generation, voltage Active power MW 32.0 600.0 -

Fixed voltage kV 220.0 220.0 220.0 -

Load Active power MW 16.0 64.0 256.0 1020.0 -

Reactive power Mvar 10.0 40.0 160.0 640.0 -

Table 2 Critical Points

I-st critical point Bus p

Voltage kV Re/Im 126.7 -27.9 123.3 182.2 50.3 181.7 2.9 219.9 24.4 173.7 220.0 0.0 133.0 -19.2 95.3 91.2

1 2 3 4 5 6 7 8

Vector y

? y0

MW Mvar or kV2 -48.5 0.0 0.0 0.0 0.0 0.0 199.8 0.0 210.4 0.0 0.0 0.0 0.0 0.0 0.0 0.0

II-nd critical point

Eigenvector MW Mvar or kV2 48.5 94.1 -238.5 -1.3 -217.7 8.2 -199.8 -0.8 -210.4 11.8 -1263.0 0.0 32.3 90.3 -188.0 109.6

Voltage kV Re/Im 81.6 -74.4 219.1 20.4 199.3 45.5 206.0 77.2 184.5 37.9 220.0 0.0 133.9 -42.7 171.0 1.0

Vector y

? y0

MW Mvar or kV2 -297.8 0.0 0.0 0.0 0.0 0.0 25.3 0.0 22.5 0.0 0.0 0.0 0.0 0.0 0.0 0.0

Eigenvector MW Mvar or kV2 297.8 232.0 -23.9 -0.5 -23.5 9.1 -25.3 -0.5 -22.5 12.5 -1826.3 0.0 64.6 115.8 -3.6 53.1

Table 3 Initial directions and convergence of the method

Initial loading directions

Convergence results

Experiment


P1 MW


P4 MW


P5 MW

Number of iterations

Solution

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

300.0 0.0 0.0 -300.0 0.0 0.0 300.0 300.0 0.0 300.0 300.0 0.0 -300.0 -300.0 0.0 300.0 -300.0 300.0 300.0 300.0 -300.0 -300.0 -300.0 -48.5 -297.8

0.0 300.0 0.0 0.0 -300.0 0.0 300.0 0.0 300.0 -300.0 0.0 300.0 300.0 0.0 -300.0 300.0 300.0 -300.0 300.0 -300.0 300.0 -300.0 -300.0 199.8 25.3

0.0 0.0 300.0 0.0 0.0 -300.0 0.0 300.0 300.0 0.0 -300.0 -300.0 0.0 300.0 300.0 300.0 300.0 300.0 -300.0 -300.0 -300.0 300.0 -300.0 210.4 22.5

24 4 4 4 7 6 4 7 8 6 4 22 5 4 10 0 0

I I I II s.p. s.p. I I I s.p. s.p. s.p. II II s.p. I I I s.p. s.p. II II II I II

The aim of the example was to determine the closest points on the power ow solution boundary  (the critical points), if the real power injections at buses 1, 4, and 5 were free parameters, i.e., allowed to vary from their operating point values. Two critical points were obtained using the algorithm of Section 4. Details of these points are given in Table 2. The length of the vector y ? y0 is 294.1MW for the rst critical point, and 299.7MW for the second critical point. In both cases the angle between y ? y0 (columns 3, 6 in Table 2) and s, which is formed from the elements of the left eigenvector s (columns 4, 7 in Table 2) that correspond to free parameters, is equal to 180deg. In obtaining these solutions, a number of dierent loading directions
y1 were used in the rst stage of the critical point algorithm. These loading directions are given in Table 3 (columns 2 to 4). The loading directions gave dierent (stage one) points on the solution boundary . Each of these points was used as the starting point for the second stage of the critical point algorithm. Table 3 (columns 5, 6) shows convergence results for the second stage when the solution technique of [6] was used. Column 6 indicates which of the critical point was converged to, or whether a singular point (s.p.) was encountered. Figure 6 shows trajectories of the second stage solution

process. Each trajectory starts from the point on  obtained from stage one for the dierent loading directions. The labels of these starting points correspond to the loading directions given in Table 3.

14 22

13

17 9

4 I

500

3

21

8

0

P5

2 16

II

7 12

23

0

6

18

5

-500

1

20 11

10

-400

600 400

-200

200

0

0

200

-200

400

-400

600 800

-600 P4

-800

P1

Figure 6: Trajectories of the stage two iterative process in the space of free parameters.

7

700

600 550

8

500 450

load ows of power systems (Fundamentals of the theory and computational methods), Publishing House of

350 Solution II

[2]

300 Solution I 250

2

3

4

5

6 Iterations

7

8

9

10

11

Figure 7: Distance changes during the stage two iterative process. The following observations were made about the stage two solution process:  The iterative solution process always moved along the solution space boundary , i.e., one eigenvalue of the power ow Jacobian J(x) was always zero.  The distance d(y) = ky ? y0 k steadily decreased as the iterative process moved from the initial point on  given by stage one, to the nal point (either of the critical points, or a singular point). This behaviour is shown in Figure 7.  In most cases solutions were obtained after 4-8 iterations. More iterations were required (10-24), and singular points were encountered, when inappropriate initial loading directions
y1 were chosen.  Convergence to trivial solutions or maxima never occured.

6

REFERENCES

[1] A.M. Kontorovich and A.V. Krukov, Stability limit

Singular points

400

200 1

ACKNOWLEDGEMENT

This work was sponsored in part by an Australian Electricity Supply Industry Research Board project grant \Voltage Collapse Analysis and Control".

650

Distance, MW

system. Points which (locally) provide this minimum distance satisfy a constrained optimization problem. The optimization problem leads to a set of equations which describe such critical points. Not all solutions of this set of equations are critical points however. Trivial solutions correspond to solutions of the usual power ow problem. Other nontrivial solutions describe extrema of the optimization problem that are not of interest. Care must therefore be taken to ensure that algorithms for nding critical points do in fact nd the correct type of points. A two stage algorithm can be used to nd critical points. Because there may be many critical points, it is necessary to provide an estimate of the direction in parameter space of the desired critical point. The rst stage of the algorithm nds a point on the solution space boundary which lies in that speci ed direction. The second stage uses a continuation method to move along the boundary from that initial point to the desired critical point. The distance from the operating point to the boundary point always decreases along the continuation path. The proposed algorithm converges reliably to desired critical points under normal conditions. However, if a very bad estimate of the direction of the critical point is used, singularity of the Jacobian of the critical point equations may occur.

CONCLUSIONS

The minimum distance from an operating point to the power ow solution space boundary, i.e., to points of maximum loadability, gives a measure of the security of a power

[3]

[4]

[5]

[6]

[7]

the Irkutsk University, Irkutsk, 1985 (in Russian). C.A. Ca~nizares and F.L. Alvarado, \Computational experience with the point of collapse method on very large AC/DC systems", Proc. NSF/ECC Workshop on Bulk Power System Voltage Phenomena II, Deep Creek Lake, MD, August 1991; published by ECC Inc., Fairfax, Virginia. I.A. Hiskens and R.J. Davy, \A technique for exploring the power ow solution space boundary", Technical Report EE9327, Department of Electrical and Computer Engineering, The University of Newcastle, Australia, May 1994. A.M. Kontorovich, A.V. Krukov, M.K. Lukina, Y.V. Makarov, V.E. Saktoev and R.G. Khulukshinov, Methods of stability indices computations for complicated power systems, Publishing House of the Irkutsk University, Irkutsk, 1988 (in Russian). I. Dobson and L. Lu, \New methods for computing a closest saddle node bifurcations and worst case load power margin for voltage collapse", IEEE Transactions on Power Systems, Vol. 8, No. 3, August 1993, pp. 905-913. Y.V. Makarov and I.A. Hiskens, \A continuation method approach to nding the closest saddle node bifurcation point", Proc. NSF/ECC Workshop on Bulk Power System Voltage Phenomena III, Davos, Switzerland, August 1994; published by ECC Inc., Fairfax, Virginia. R. Seydel, From Equilibrium to Chaos, Elsevier Science Publishing Co., New York, 1988.