Smart Grid and Renewable Energy
Vol. 4  No. 1 (2013) , Article ID: 28122 , 8 pages DOI:10.4236/sgre.2013.41005

Identification of the Optimum Locations of Power System Stabilizers in a Multimachine Power System Using Second Order Eigenvalue Sensitivity Analysis

Abdul Mahabuba1, Abdullah Khan2

1College of Engineering and Applied Sciences, Al Ghurair University, Dubai, UAE; 2Department of Electrical and Electronics Engineering, B.S. Abdur Rahman University of Science and Technology, Chennai, India.

Email: mifarz@yahoo.com

Received April 8th, 2012; revised November 8th, 2012; accepted November 15th, 2012

Keywords: Power System Stabilizer (PSS); Small Signal Stability; First Order Eigenvalue Sensitivity (FOES); Second Order Eigenvalue Sensitivity (SOES); Participation Factor (PF)

ABSTRACT

This research paper deals with the identification of the best location of the Power System Stabilizers (PSS) and also the tuning of PSS parameters in order to improve the overall dynamic stability of multi machine power systems. The location of PSS is determined by identifying the critical modes and their corresponding first and second order eigenvalue sensitivities. In this formulation, sensitivity analysis of a particular mode can be performed with only its eigenvalues and their left and right eigenvectors. The simplicity and efficiency of this approach sharply contrast to the complexity of the traditional approach, where all eigenvalues and eigenvectors are required at the same time. The effectiveness of this method in selecting the optimum location for placement of PSSs is compared with the participation factor method. The proposed sensitivity theory used to identify the best PSS location in a five machine, eight bus El-Metwally and Malik System to increase the damping of both local and inter area modes for various operating conditions.

1. Introduction

Power systems are highly complex systems that contain non-linear and time varying elements. Many power systems face the problem of troublesome dynamic oscillations in the range of 0.2 to 2.5 Hz associated with some poorly damped swing modes. Power system stabilizers (PSSs) are commonly used to damp these oscillations and increase damping of swing modes [1]. In the application of PSS to increase the damping of a certain troublesome mode in a multi machine power system the very first step is to determine the best PSS location(s). For a local mode the job of selecting PSS location(s) is not difficult because only few machines are involved in this local oscillation and there are only a few choices. But for the inter area mode large number of machines may be involved in the oscillation. This makes the PSS location selecting problem very complicated [2]. The most commonly used approach for stability siting was the Eigen vector method proposed by DeMello et al. [3]. This method too may fail in certain circumstances and may lead to undesirable stabilizer location as found out in various studies. Hiyama [4] presented an approach using the concept of coherent groups. This method suffers from the major disadvantage that generators within one coherent group for a large disturbance do not necessarily remain in the same group under small-disturbance conditions.

Zhang and Coonick [5] proposed a new method based on the method of inequalities for the coordinated synthesis of PSSs parameters in multi-machine power system in order to enhance overall system small signal stability. Antonio et al. [6] presented a method that simultaneously tune multiple power system damping controllers using Genetic Algorithms. Cai and Elrich [7] suggested the simultaneous coordinated tuning of the series FACTS Power Oscillation Damping controller in multi-machine power system. Doi and Abe [8] developed a new coordinated synthesis method by combining eigenvalue sensitivity analysis and linear programming applied to this method is used to synthesize the coordination of power system stabilizers in a new multi machine system.

In this paper, the first order and second order eigenvalue sensitivity theory presented in [9-12] used for selection of optimum locations of PSS and the same is compared with the participation factor method. This method is used to identify the best PSS locations in a five-machine, eight bus system (El-Metwally and Malik System) to increase the damping of an inter area mode for three different operating conditions. The dynamic response of rotor angle deviation of the system confirms that the prediction of the best PSS location by the proposed method is correct and accurate. Section 2 discusses about first order and second order eigenvalue sensitivity analysis. Section 3 presents the dynamic stability model of multimachine power systems. Section 4 discusses the proposed method of optimum locations of PSSs. Section 5 provides the simulation results and Section 6 gives the conclusion.

2. Eigenvalue Sensitivity

The linearized model is as shown in Equation (1)

(1)

The eigenvalues of a matrix are given by non-trivial solutions to the equation :

(2)

where A is an n x n matrix and is a nx1 vector.

For any eigenvalue, the n column vector, which satisfies Equation (2), is called the right eigenvector of A associated with the eigenvalue.

Therefore we have,

(3)

Similarly, the n-row vector which satisfies

(4)

which is called the left eigenvector associated with the eigenvector.

Consider the equation,

(5)

Sensitivityof eigenvalues with respect to elements of,

(6)

Premultiplying the Equation (6) by & noting that & results,

All elements of for all and 1 for

(7)

Differentiate Equation (7) with respect to parameter gives Equation (8):

(8)

Differentiating Equation (8) with respect to another PSS parameter, then simplify yield the second order eigenvalue sensitivity [14,15] as below:

(9)

The Equation (9) requires knowledge of all eigenvalues and eigenvectors, which is very large problem to solve when the system matrix is very large. To find the second order eigenvalue sensitivities without having to know all the eigenvalues and eigenvectors, the second order eigenvalue sensitivities can be written in terms of two first order eigenvector sensitivities [4] as Equation (10):

(10)

If the eigenvector sensitivities and are known, then the above equation can easily be computed.

3. Dynamic Stability Model of Multi Machine Power System

In stability analysis of a multi-machine system, modeling of all the machines in a more detailed manner is exceedingly complex in view of the large number of synchronous machines to be simulated. Therefore simplifying assumptions are usually made in modeling the system. In this studies two axis model is used for all machines in the test system.

The linearized state equations for the two axis model under the assumption are given by Anderson [13] as follows:

(11)

where the state variables for the machine are

-direct axis component of voltage behind transient reactance

-quadrature axis component of voltage behind transient reactance

-angular velocity of rotor

-rotor angle in radians The state variables of the machine including the exciter and PSS are given by

(12)

A five machine, eight bus El-Metwally and Malik test system (Figure 1) is chosen to investigate the optimum location of the PSS and the tuning of the PSS parameters. The system data is taken from [1].

4. Proposed Method for the Optimum Location of PSSs

The step-wise procedure for identifying the best location of PSS as detailed below.

Step 1: Linearize the system at a chosen operating condition and from the linearized system matrix, identify the rotor swing modes including the nature of modes whether local or inter area modes using eigen sensitivity analysis.

Step 2: Determine the critical rotor swing modes which have damping ratio less than a specified value, say

Figure 1. Test system.

0.45.

Step 3: Steps 1 and 2 are repeated for each chosen operating condition.

Step 4: The critical swing modes of all the operating conditions identified are arranged in the ascending order of the damping ratio. To each critical swing mode associate an operating condition; this has the least damping ratio.

Step 5: The best location for damping the first critical swing mode in the list is obtained using the eigen value sensitivity and participation factors.

Step 6: Step 5 is repeated for all other critical swing modes.

The following three operating conditions are chosen for analysis of the test system:

1) Full load operating condition.

2) Operating condition 1) with load increased by 25%.

3) Full load operating condition with one of the two tie lines between buses 6 and 7 removed.

From Table 1, it is seen that there are four critical swing modes. Mode shapes are obtained by plotting the speed component of the right Eigen vectors of the swing modes as in Figure 2.

From the mode shapes, the first mode is the inter machine oscillation local to area 1, with G2 & G3 swinging

Table 1. Eigenvalues—Operating condition (a).

Figure 2. Mode shapes of swing modes.

against G5. The second mode is the inter machine oscillation local to area 1, with G2 swinging against G3 & G5. The third mode is the inter machine oscillation local to area 2, with G1 swinging against G4. The fourth mode is an inter-area mode, with generators G2, G3, G5 of area 1 swinging against generators G1 and G4 of area 2. The critical modes with their damping ratio, nature of the modes, first order eigen sensitivity value for the three operating conditions are shown in Tables 2-4 respectively.

The critical swing modes identified in all the three operating conditions are arranged in the ascending order of the damping ratio. The operating conditions and related details of the four critical rotor swing modes are given in Table 5.

From Table 5, for the first critical mode, machine 2 is having the highest first order sensitivity (Table 2) and PSS is located at machine 2. Similarly for the second critical mode, machine 2 is having the highest first order sensitivity (Table 2) and PSS is located at machine 2. For the third critical mode, machine 4 is having the highest first order sensitivity (Table 2) and PSS is located at machine 4. For the fourth critical mode, machine 2 is having the highest first order sensitivity (Table 3) and PSS is located at machine 2.

Table 2. Modes and first order eigenvalue sensitivities— operating condition (a).

Table 3. Modes and first order eigenvalue sensitivities— operating condition (b).

Table 4. Modes and first order eigenvalue sensitivities— operating condition (c).

The same procedure outlined in first order is repeated for second order eigenvalue sensitivities to determine the optimum location of PSS. According to second order eigenvalue sensitivity, the results reveal that there are two oscillatory rotor swing modes for each of the three operating conditions which are shown in Tables 6-8 for the operating conditions 1), 2) and 3) respectively.

Table 5. Ranking—Critical swing modes (FOES).

Table 6. Modes and second order eigenvalue sensitivities— operating condition (a).

Table 7. Modes and second order eigenvalue sensitivities— operating condition (b).

The critical swing modes identified in all the three operating conditions are arranged in the ascending order of the damping ratio. The operating conditions and other related details of the two critical rotor swing modes are given in Table 9.

From Table 9, for the first critical mode, machine 4 is having the highest second order sensitivity (Table 7) and PSS is located at machine 4. Similarly for the second critical mode, machine 2 is having the highest second order sensitivity (Table 6) and PSS is located at machine 2. Hence PSSs are to be located at machines 2 and 4 for two critical rotor swing modes.

The tuning of PSS parameters using Conventional technique has been given below. The gain and time constants of CPSS for the five machines for three operating condition 1), 2), 3) are shown in Table 10.

The effectiveness of the proposed method of location and tuning are investigated by carrying out simulation of the linear model of the system for three operating conditions with disturbances.

Table 8. Modes and second order eigenvalue sensitivities— operating condition (c).

Table 9. Ranking—Critical swing modes (SOES).

Table 10. Tuned parameters of PSS.

5. Simulation Results

The application of the FOES and SOES is demonstrated by applying it to the five-machine, eight bus system. The time domain simulations of rotor angle deviation of the system are carried out without PSS and with PSSs located on machines 2 & 4 after tuning the PSS parameters. The corresponding system response of initial perturbation of 5˚ in rotor angle at machine 2, 3, 4 and 5 with respect to machine 1 for operating condition 1) and 2) are shown in Figures 3(a)-(d) to 4(a)-(d) respectively.

Similar type of responses was occurring for the operating condition 3) also. The effectiveness of the first and second order eigenvalue sensitivities methods is compared with the participation factor method in selecting the optimum location for placement of PSSs [16]. Figures 5(a)-(c) are the comparison of the system responses for operating conditions 1), 2) and 3) under 5˚ rotor angle perturbation.

6. Conclusions

In this paper, an eigen-sensitivity of the system matrix is developed using a dominant eigenvalue and the corresponding eigenvectors. The proposed sensitivity can provide valuable information in enhancing power system stability identifying the siting of compensation device PSS. The results of the method indicate that machines 2 and 4 are the optimum locations for installing PSSs.

The eigenvalue sensitivity analysis is used to find out optimum location of PSS for different operating conditions and gives the best damping of all the critical modes of the system. The work in this paper is concerned with the location Power system stabilizers in multimachine power systems and tuning its parameters. This work has got great practical relevance for the modern power systems, where the ever increasing demand for the electrical power and the integration of large utilities have been regularly causing the system to operate at the threshold of its stability limits, which calls for adequate precautionary measures to be taken to ensure rotor angle stability of the synchronous machines. An efficient computer-programming package for analyzing the dynamic stability of power systems has been developed in MATLAB. Using the first and second order sensitivity the

(a)(b)(c)(d)

Figure 3. (a) System response of -op.condition (a) with an initial perturbation of 5˚; (b) System response of -op.condition (a) with an initial perturbation of 5˚; (c) System Response of -op.condition (a) with an initial perturbation of 5˚; (d) System response of -op.condition (a) with an initial perturbation of 5˚.

(a)(b)(c)(d)

Figure 4. (a) System response of -op.condition (b) with an initial perturbation of 5˚; (b) System response of -op.condition (b) with an initial perturbation of 5˚; (c) System Response of -op.condition (b) with an initial perturbation of 5˚; (d) System response of -op.condition (b) with an initial perturbation of 5˚.

(a)(b)(c)

Figure 5. (a) System Response of -op.condition (a) with an initial perturbation of 5˚; (b) System Response of - op.condition (b) with an initial perturbation of 5˚; (c) System Response of -op.condition (c) with an initial perturbation of 5˚.

machines creating the critical swing modes are identified, which can be proposed as the optimum locations for placing the power system stabilizers.

REFERENCES

  1. P. Lakshmi and M. A. Khan, “Design of a Robust Power System Stabilizer Using Fuzzy Logic for a Multi-Machine Power System,” International Journal on Electric Power Systems Research, Vol. 47, No. 4, 1998, pp. 39-46.
  2. E. Z. Zhout, O. P. Malik and G. S. Hope, “Theory and Method for Selection of Power System Stabilizer Location,” IEEE Transactions on Energy Conversion, Vol. 6, No. 1, 1991, pp. 170-176. doi:10.1109/60.73804
  3. De Mello, et al., “Coordinated Application of Stabilizers in Multi-Machine Power Systems,” IEEE Transactions on Power Apparatus and Systems, Vol. PAS-99, No. 3, 1980, pp. 227-230.
  4. T. Hiyama, “Coherency Based Identification of Optimum Site for Stabilizer Applications,” IEE Proceedings C on Generation, Transmission and Distribution, Vol. 130, No. 2, 1983, pp. 71-74.
  5. P. Zhang and A. H. Coonick, “Coordinated Synthesis of PSS Parameters in Multi-Machine Power Systems Using the Method of Inequalities Applied to Genetic Algorithms,” IEEE Transactions on Power Systems, Vol. 15, No. 3, 2000, pp. 811-816. doi:10.1109/59.867178
  6. A. L. B. Do Bomfim, G. N. Taranto and D. M. Falcao, “Simultaneous Tuning of Power System Damping Controllers Using Genetic Algorithms,” IEEE Transactions on Power Systems, Vol. 15, No. 1, 2000, pp. 163-169. doi:10.1109/59.852116
  7. L. J. Cai and L. Elrich, “Simultaneous Coordinated Tuning of PSS and FACTS Controller for Damping Power System Oscillations in Multi-Machine Systems,” 2003 IEEE Bologna Power Tech Conference Proceedings, Vol. 2, 23-26 June 2003, pp. 1025-1030.
  8. A. Doi and A. Abe, “Coordinated Synthesis of Power System Stabilizers in a Multi-Machine Power System,” IEEE Transactions on Power Apparatus and Systems, Vol. 103, No. 6, 1984, pp. 1473-1479. doi:10.1109/TPAS.1984.318486
  9. H.-K. Nam, Y.-K. Kim, K.-S. Shim and K. Y. Lee, “A New Eigen-Sensitivity Theory of Augmented Matrix & Its Applications to Power System Stability Analysis,” IEEE Transactions on Power Systems, Vol. 15, No. 1, 2000, pp. 363-369. doi:10.1109/59.852145
  10. P. Kundur, “Power System Stability and Control,” McGraw Hill Inc., New York, 1994.
  11. H. M. Zein El-Din and R. T. H. Alden, “Second Order Eigenvalue Sensitivites Applied to Power System Dynamics,” IEEE Transactions on Power Apparatus and Systems, Vol. PAS-96, No. 6, 1977, pp. 1928-1936. doi:10.1109/T-PAS.1977.32528
  12. J. V. Milanovie, “Identification of Electromechanical Modes and Placement of PSSs Using Relative Gain Array,” IEEE Transactions on Power Systems, Vol. 19, No. 1, 2004, pp. 410-417. doi:10.1109/TPWRS.2003.821454
  13. P. M. Anderson and A. A. Fouad, “Power System Control and Stability,” Iowa State University Press, Iowa, 1997.
  14. C. Y. Chung, K. W. Wang, C. T. Tse, X. Y. Bian and A. K. David, “Probabilistic Eigenvalue Sensitivity Analysis and PSS Design in Multimachine Systems,” IEEE Transactions on Power Systems, Vol. 18, No. 4, 2003, pp. 1439-1445. doi:10.1109/TPWRS.2003.818709
  15. C. Y. Chung, K. W. Wang, C. T. Tse and R. Niu, “Power System Stabilizer Design by Probabilistic Sensitivity Indexes,” IEEE Transactions on Power Systems, Vol. 17, No. 3, 2002, pp. 688-693. doi:10.1109/TPWRS.2002.800914
  16. Y. Y. H Su and C. L. Chen, “Identification of Optimum Location for Stabilizer Applications Using Participation Factors,” IEE Proceedings, Vol. 134, No. 3, 1987, pp. 237-245.