Abstract

Transient stability assessment (TSA) of the power system is a crucial issue with escalating demands and large operational constraints. Real-time TSA allows for deciding and monitoring of the relevant preventive/corrective control actions depending on the dynamic behavior of the system components. To assess this, coherency of generating machines is to be found. After determination of the coherent machines, any corrective or preventive action can be initiated by the system operator to maintain stability of the system during occurrence of any severe contingency. The Transient Severity Index (TSI) introduced in this paper has proven to be an interesting alternative for determining generator coherency. Furthermore, the numerical values of this index are employed to construct a supervised learning-based classifier and the ranking method with the help of system load and generation as input features. This framework employs the support vector machine (SVM) to perform the ranking of the generators based on severity and classify them into vulnerable and nonvulnerable machines. The results are validated on the IEEE 10-generator, 39-bus test (New England) system. It is observed that the proposed index and the supervised learning engine give satisfactory results and both are aligned with the published approaches.

1. Introduction

Modern power network is emerged from the vertically integrated structure to the deregulated unbundled structure. This transformation enables consumers and generators to participate in the competitive business environment. Issues of power quality, reliability, risk, and security become prominent in this environment. High power demand has made an essential criterion to maintain the power system secure under the stressed operating conditions. Transient stability (T/S) of the power system is an important denominator and defined by the CIGRE committee [1] as “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.” Stability studies can be subdivided into two categories: the first one is assessment and the second one is control. Ensuring T/S under any unforeseen contingency is a challenge for a system operator. This fact initiated the research in the direction of assessment methods along with preventive control methods for ensuring system stability.

In the past, various methods have been applied by researchers to assess system stability. Stability analysis-oriented research is mainly based on the evaluation of the violation of the limit of power system quantities and the behavior of the generator under contingency conditions. Some methods [25] proposed for stability assessment in the latest literature were based on risk or probabilistic, pattern recognition, and time domain simulation (TDS) techniques. TSA is directly concerned about the large disturbance. A risk-based methodology was proposed in [24]. A probabilistic small-disturbance stability analysis (PSSA) method with uncertainties was proposed by Preece and Milanović [2]. In this approach, the probability density functions (PDFs) were designed for the damping of the critical oscillatory electromechanical modes by modeling the stochastic variation of system uncertainties. However, this article did not give any analysis for large disturbances. A linear function-based technique was used in [3] to calculate a risk-based index for TSA. A new technique to identify the severity of the power substation under possible operating scenarios was proposed by Da Silva et al. [4]. The research work presented in [5] was based on the pattern discovery-based fuzzy classification. The pattern identification algorithm was modified by using centroid deviation analysis and predisturbance information of the training data set. This centroid deviation analysis-based technique is very difficult for practical implementation. The Lyapunov stability criteria were given in [6] for differential/algebraic models of power systems which include the effect of generator damping and nonlinear loads. Cheung and Chow [7] investigated the simulation of slow dynamics in two-time-scale power systems. In this work, an approach was proposed to obtain the slow dynamics by projecting the trajectory of the postfault system onto its slow manifold. Direct stability analysis of AC/DC power systems using a structure-preserving energy function (SPEF) was proposed in [8]. Bhat et al. [9] gave a concept of the single machine equivalent (SIME) T/S method combined with Simulink to model T/S phenomena and the Graphical User Interface Design Environment (GUIDE) software of MATLAB. A method based on the network decomposition algorithm was presented in [10] to implement simultaneous unbalances in power system T/S.

In the current decade, many researchers have proposed stability indices for assessment of stability of large power systems. Transient security indices determined in [11, 12] can be used to measure the synchronism state of a generating unit. A small signal stability index was proposed in [13] for power network dynamic assessment by employing TDS. The value of this index was calculated by the systems’ eigenvalues using dynamic simulation. Several indices were presented in [14, 15] for TSA. A practical and heuristic index based on T/S in large power systems was proposed in [14] for fast contingency ranking. A didactic procedure based on simulation for T/S of a multimachine power system was proposed in [16]. Some of the techniques presented in the literature were based on the element T/S, angle stability index, frequency stability index, and transient voltage stability index.

As stated, generator coherency has a considerable application in power system operation and control. The sole concept of coherent groups of generators in a power system is the similar behavior of TDS responses of generators, when they are subjected to a perturbation [17]. This is called coherency. In this consideration, generators that have the similar postdisturbance rotor angle or speed variation characteristics are called coherent and are placed in the similar coherent group. Hence, it is required to monitor and examine the relationship or similarity of rotor angle deviation to find the coherent nature of the power system components. The generators’ dynamic response under disturbances can also be recognized by the deviations in phase angles of voltage or/and current phasors of the system.

Several methods have been introduced in the literature for identification of coherent behavior of the generators and classification of them, according to their similar behavior (TDS characteristics). After identification of the coherent groups, control strategies can be applied on them. In order to initiate any preventive control under stressed condition, it is desirable to find the coherency between generators. In general, coherency classification techniques can be divided mainly into two types. The techniques that are placed in the first type are based on model reduction and required computation of the eigenvalues and eigenvectors of the power system [18]. For example, a synchronic modal equivalencing (SME) was proposed [18] for structure-preserving dynamic equivalencing of large power system models. Techniques that are placed in the second type are based on disturbances and use TDS to find coherent groups of generators. For example, several studies have used the rotor trajectory index (RTI) [12], Fourier spectrum [19], or fast Fourier dominant interarea mode [20], principal component analysis [21], independent component analysis [22], hierarchical clustering methods [23, 24], fuzzy c-medoids algorithm [25, 26], wavelet [27], and Hilbert–Huang transform [28] to find coherent generators. Several techniques proposed in literature cannot always be reliable for all applications, such as system protection schemes or remedial action schemes, dynamic equivalencing, and controlled islanding. This maloperation is found due to the degree of coherency between generators or nongenerating buses. This may change for different disturbances at different power system operating conditions.

Generally, real-time TSA is an approach to find the fast and accurate prediction of the system stability status (either stable or unstable) in real time by considering the future behavior of the generator under the disturbed operating condition. In the literature, various topologies were proposed to forecast the system stability states by using different types of artificial intelligence-based techniques [2932]. Hashiesh et al. proposed a supervised learning technique for T/S state prediction by the artificial neural network (ANN) [30]. The SVM-based techniques [31] were used for power system contingency classification and ranking in [29, 32]. The real-time TSA method using the SVM was also proposed in [33]. The prediction of the real-time transient state of the power can be done successfully by the above-discussed methods, but these methods are used only as a classifier. The stability status of the individual generator for any contingency cannot be determined by these methods. For the control action like generator rescheduling, knowledge of the individual generator state (either stable or unstable) is mandatory. To overcome these deficiencies, we propose an SVM-based method to predict the status of the individual generator and overall system by using the coherent characteristics of the generators.

In this paper, the T/S status of the power system is analyzed by simulating the power system under different operating conditions, on different fault locations, and for different fault clearing times. The status of system stability is identified by using a new index, Transient Severity Index (TSI). This index is based on rotor angle deviation of the generators in the postfault state. A supervised learning engine based on the LS-SVM is proposed to rank the generators according to the values of TSIs. The potential features for this supervised learning engine are extracted from TDS of the power system under dynamic operating conditions. These potential features are real, reactive power at generating and load buses.

The remaining part of this paper is organized as follows: Section 2 describes power system modeling, Section 3 contains the mathematical framework of the proposed TSI along with the model of the LS-SVM. In Section 4, simulation and results are presented, while Section 5 highlights major conclusions drawn from the work and also presents a solid milieu for future work.

2. Power System Modeling

The TDS method is used to solve the differential equations involved in power system stability analysis. TDS of the power system is required for DSA to find whether the system can maintain synchronism during disturbance or not. This decision is taken by monitoring the movement of trajectories of rotor angles during a perturbation period. The swing equation shows the transient behavior of the system. If the trajectories of the rotor angle of either a single generator or a group of generators are found to have continuous increment without limit with reference to remaining machines, then the system is unstable. Another phenomenon is if rotor angles of all working system generators remain bounded within their respective permissible limits, then the system is stable [17, 34]:where is the transfer admittance between the and generators, is the moment of inertia, is the electrical power output, is the damping coefficient, is the mechanical power input, is the rotor speed deviation, is the rotor angle of the machine, and is the rotor angle of the machine. The deviation between generator rotor angles with reference to time can be found by using the set of swing equations (1)–(3):

In the projectile transient phenomena, the inertial center of the system is taken as the reference frame for the calculations. The machines’ rotor angles with respect to center of inertia [35] are used to detect whether the system is stable or not. For an machine system with the inertia constant and rotor angle of the machine, the inertial center is determined as per (4).

3. Proposed Methodology

In this work, a numerical integration technique (Runge–Kutta method) is used to solve the swing equation. The transient security status is determined by monitoring the swing in rotor angle trajectories and deviation in rotor angle with respect to the constraint for transient instability. The constraint is given as follows:where is the maximum allowable value of the relative rotor angle for secure operation. During the case of disturbance, if the relative rotor angle violates in a time interval , then the system is considered as insecure , else it is considered as secure . For this work, the maximum allowable value of the relative rotor angle for secure operation is taken as [36].

The ranking of the generating machines for every insecure contingency needs to be discovered with less computational burden and time. Therefore, the rotor angle trajectory-based severity index, called the Transient Severity Index (TSI), is proposed in this paper. The TSI is determined from TDS, and the process of determination is shown in Figure 1. The TSI can be used to assess the stability of the power system and to rank the criticality of the generators. Hence, it gives the synchronizing condition for a given hard contingency. In this work, we present a method to rank the machines by employing the TSI. Equation (6) shows the mathematical formulation of the TSI:where is the value of rotor angle deviation in degrees, is the trajectory constant, and is the time at which is measured. for all the generators is noted after of fault clearing time (FCT). A classification boundary is shown in Table 1. The numerical value of the TSI is an indicator of the unstable or stable system state of the power system.

The generators can be divided into three groups: the first is highly critical (most advanced) generators (CGs), the second is advanced generators, and the third is noncritical (least advanced) generators (NGs). Ranking of generators can be found by their corresponding generator TSI value. Furthermore, this ranking can be used in the process of preventive and corrective control action.

3.1. Sensitivity Analysis of K

In this section, sensitivity analysis of the trajectory constant (K) is presented. On the basis of the values of K, the TSI can take any three values as per (8)–(10).

Case 1:

Case 2:

Case 3:

The variation of and TSI with respect to is observed, and the results of this analysis are exhibited in Figure 2. It shows TSI trajectories for nine different values of . varies from 120° to 360° in the step of 30°. By inspecting the variation of the TSI values with respect to for the constant value of , that higher value of reduces the span of the index. As the value of decreases, the span of the index increases. It is a basic characteristic to design a scalar index. Hence, the trajectory at K = 120° (i.e., ) exhibits the largest span to classification and is used for further analysis.

3.2. Design of the Supervised Learning Engine

The values of the rotor angle are viewed as a multidimensional vector (a list of index values), and a classification boundary (hyperplane) is required to classify the generators in the abovementioned three coherent groups. Design of the SVM has two aspects: the first is to reduce the number of features by any feature selection method and the second is to find the largest margin (hyperplane) between both the states of the system (stable and unstable). For handling these aspects, particle swarm optimization (PSO) [3739] is employed in this work.

The TSI can be used to assess the stability of the power system and rank the generator. A supervised learning engine based on the LS-SVM is proposed to predict the coherent group and rank the generators according to the values of the TSI as depicted in Figure 3. In LS-SVMs, the input data are mapped in a linear form with the high-dimensional feature space with the help of kernel functions. To solve the classification problem, least-square loss function is used in this paper.

The input features vector for this engine consisting of real and reactive power at generator buses ( and ) and load buses ( and ) is defined by the following equation:where , in which and represent the total no. of generator buses and load buses, respectively.

Feature selection aims to find the most important information from a given set of features. As this task can be seen as an optimization problem, the combinatorial growth of the possible solutions may be inviable for an exhaustive search. In this paper, PSO is used for feature selection and SVM is used for maximum margin classifier design [40]. A mean of standard deviation-based objective function is optimized by using PSO. In this approach, PSO creates different data packets with random feature selection and obtains the standard deviation of each data packet and maximizes the mean of the standard deviation of different data packets. The final feature vector is that which has the maximum mean or large Euclidean distance of features. The vector of potential input features for training of the SVM can be found from the following equation:

The target for the SVM output is the generator severity rank and coherent group of all the generators at a given time instant, which is defined as target vector as shown in the following equation:

Authors simulated 1,000 different cases; out of these, 85% are used for training and validation and the remaining 15% are used for testing purpose.

4. Simulation Results

The proposed strategy is tested on the global standard benchmark IEEE 10-generator and 39-bus (New England) power network as shown in Figure 4. The system has 46 transmission lines, 39 buses, 10 generators, and 12 transformers. Bus 39 is considered as the slack bus during simulations. All the TDS works are carried out using MATLAB [41], MATPOWER [42], and PSAT Version 2.1.10 [43] with Intel CoreTM i3-3110M platform with 2.40 GHz and 6 GB RAM. The symmetrical fault with ground is initiated at random time from 1.1 s to 1.2 s and cleared at 1.3 s with different random perturbations of the real and reactive loads on all the buses and real and reactive power generation at the generator buses. System total load is varied from 95% to 110% of the base case total load. Fault is applied on all types of buses. To demonstrate the effectiveness of the proposed method, six typical severe contingencies have been taken into consideration as described in Table 2 for studies. The proposed procedure is utilized to identify coherency information for generation rescheduling to initiate preventive control. The experimental procedure and its evaluation methodology are as follows:Step 1: For any insecure contingency, after the observation of rotor angle trajectories, three coherent groups are made to classify all the generators according to the value of the TSI as in Table 1.Step 2: All the critical generators (most advanced generators), also considered in a group of critical machines, form the generation decreasing group.Step 3: All noncritical machines or machines of the least advanced generator group form the generation increasing group.Step 4: After the coherency and stability assessment, initiate the preventive control as suggested by the proposed method.

Table 3 shows that the results (rank and coherent group) obtained from the LS-SVM are aligned with RTI [12]. For the justification of these results, TDS is also performed, and the results of six cases are presented below.

Case I is shown for 97.5% random load variation, and three-phase fault at bus 13 is initiated at 1.1272 s and cleared at 1.3 s by opening the breakers to isolate the line 13-14. The state of T/S for all generators is determined by observing rotor angle swings in the time interval. As per Table 3, TSI value of the generators G1 and G2 is violate stability criteria (i.e., TSI > 1), so these generators are considered as unstable generators. Figure 5 shows the relative rotor angle of all the generators with respect to COI. The rotor swing of the generators G2 (at bus 31) and G3 (at bus 32) violates the threshold value, causing them to lose synchronism. Both the generators are going out of step with respect to the other generators, making the whole system transiently unstable.

Case II is shown for 110% random load variation, and three-phase fault at bus 7 is initiated at 1.1498 s and cleared at 1.3 s by opening the breakers to isolate the line 6-7. As per Table 3, TSI values of the generators G4 and G5 violate stability criteria (i.e., TSI > 1), so these generators are considered as unstable generators. Classification results for G3 obtained from the RTI and TSI are conflicting. Figure 6 shows the relative rotor angle of all the generators with respect to COI. The rotor swing of the generators G4 (at bus 33) and G5 (at bus 34) rises above the threshold value, causing them to lose synchronism. Both the generators are going out of step with respect to the other generators, making the whole system transiently unstable. However, the generator G3 (at bus 32) is also going out of the stability equilibrium, but after FCT, it regains the stable state after a few cycles. TDS shown in Figure 6 indicates that this generator exhibits self-healing characteristics and does not pose any threat to system stability.

Case III is shown for 95% random load variation, and three-phase fault at bus 10 is initiated at 1.1468 s and cleared at 1.3 s by opening the breakers to isolate the line 10-11. As per Table 3, TSI value of the generator G3 is more then stability criteria (i.e., TSI > 1), so this generator is considered as an unstable generator. Figure 7 shows the relative rotor angle of all the generators with respect to COI. The rotor swing of the generator G3 (at bus 32) violates the threshold value, causing it to lose synchronism. This generator is going out of step with the rest of the generators, making the whole system transiently unstable. However, the generator G2 (at bus 31) is also going out of step, but after FCT, it regains its stable state after a few cycles. TDS shown in Figure 7 indicates that this generator exhibits self-healing characteristics and does not pose any threat to system stability.

Case IV is shown for 105% random load variation, and three-phase fault at bus 14 is initiated at 1.1917 s and cleared at 1.3 s by opening the breakers to isolate the line 14-15. The state of T/S for all generators is determined by observing rotor angle swings in the time interval. As per Table 3, TSI values of the generators G4 and G5 are more than stability criteria (i.e., TSI > 1), so these generators are considered as unstable generators. Figure 8 shows the relative rotor angle of all the generators with respect to COI. The rotor swing of the generators G4 (at bus 33) and G5 (at bus 34) violates the threshold value, causing them to lose synchronism. Both the generators are going out of step with respect to the other generators, making the whole system transiently unstable.

Case V is shown for the base case load, and three-phase fault at bus 26 is initiated at 1.1610 s and cleared at 1.3 s by opening the breakers to isolate the line 26-27. As per Table 3, TSI value of the generator G9 is more than stability criteria (i.e., TSI > 1), so this generator is considered as an unstable generator. Figure 9 shows the relative rotor angle of all the generators with respect to COI. The rotor swing of the generator G9 (at bus 38) violates the threshold value, causing it to lose synchronism. This generator is going out of step with the other generators, making the whole system transiently unstable.

Case VI is shown for 105% random load variation, and three-phase fault at generator no. 5, at bus 34 is initiated at 1.1052 s and cleared at 1.3 s. The state of T/S for all generators is determined by observing rotor angle swings in the time interval. As per Table 3, TSI value of the generator G5 is more then stability criteria (i.e., TSI > 1), so this generator is considered as an unstable generator. Figure 10 shows the relative rotor angle of all the generators with respect to COI. The rotor swing of the generator G5 (at bus 34) violates the threshold value, causing it to lose synchronism. This generator is going out of step with the rest of the generators, making the whole system transiently unstable. To make the system stable for this case, the proposed method suggests rescheduling of the generating units.

Results shown in Table 3 are validated through nonlinear TDS as shown in Figures 510. It is observed that the results of TDS and the supervised learning architecture are aligned with each other. This comparison establishes the efficacy of the proposed approach.

Figure 11 shows the value of standard deviation for different values of the trajectory constant. It has been observed that the standard deviation is highest for K = 120°, which indicates that the Euclidean distance of the samples is higher for K = 120° as compared to the others. Hence, these cases are easily separable by a linear discriminant function.

5. Conclusion

Power system stability is a matter of concern in the present deregulated scenario. Fast and accurate stability assessment methods are more meaningful in the deregulated and dynamic framework of the power system. The problem of stability assessment has been addressed in this paper. Following are the major conclusions of this work:(i)The index TSI has been proposed in this paper. This index can serve two objectives: the first is to classify the generators in coherent groups and the second is to rank them according to the severity level.(ii)A supervised learning engine has been presented. This engine employs the LS-SVM method to rank and classify the generators. A swarm intelligence-based well-known technique PSO has been applied to perform feature selection and SVM to design maximum margin classification. It has been observed that the results of the proposed supervised engine are aligned with the numerical values of TSIs obtained using other approaches.(iii)The results of this supervised engine have been validated through nonlinear TDS. It has been observed that the results of the proposed index and supervised engine are aligned with each other. Furthermore, the TDS results verify the effectiveness of this approach.

The application of different feature selection methods and different topologies of neural architectures lies in the scope of the future work.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge the support and encouragement provided by the authorities of the Malaviya National Institute of Technology, Jaipur, and Swami Keshvanand Institute of Technology, Jaipur.