# Transient Stability Enhancement of Multimachine Power System Using Robust and Novel Controller Based CSC-STATCOM

Authors: Sandeep Gupta; Ramesh Kumar Tripathi
Journal: Advances in Power Electronics
 (2015)

---

## Abstract

A current source converter (CSC) based static synchronous compensator (STATCOM) is a shunt flexible AC transmission system (FACTS) device, which has a vital role as a stability support for small and large transient instability in an interconnected power network. This paper investigates the impact of a novel and robust pole-shifting controller for CSC-STATCOM to improve the transient stability of the multimachine power system. The proposed algorithm utilizes CSC based STATCOM to supply reactive power to the test system to maintain the transient stability in the event of severe contingency. Firstly, modeling and pole-shifting controller design for CSC based STATCOM are stated. After that, we show the impact of the proposed method in the multimachine power system with different disturbances. Here, applicability of the proposed scheme is demonstrated through simulation in MATLAB and the simulation results show an improvement in the transient stability of multimachine power system with CSC-STATCOM. Also clearly shown, the robustness and effectiveness of CSC-STATCOM are better rather than other shunt FACTS devices (SVC and VSC-STATCOM) by comparing the results in this paper.


---

## Body

## 1. Introduction

The continuous developments of electrical loads due to the modification of the society structure result in today’s transmission structure to be faced close to their stability restrictions. So the renovation of urban and rural power network is more and more necessary. Due to governmental, financial, and green climate reasons, it is not always possible to construct new transmission lines to relieve the power system stability problem at the existing overloaded transmission lines. As a result, the utility industry is facing the challenge of efficient utilization of the existing AC transmission lines in power system networks. So transient stability, voltage regulation, damping oscillations, and so forth are the most important operating issues that electrical engineers are facing during power transfer at high levels.In above power quality problems, transient stability is one of the most important key factors during power transfer at high levels. According to the literature, transient stability of a power system is its ability to maintain synchronous operation of the machines when subjected to a large disturbance [1]. While the generator excitation system with PSS (power system stabilizer) can maintain excitation control and stability it is not adequate to sustain the stability of power system for large faults or overloading occurs near to generator terminals [2].So many researchers worked on this problem in finding the solution for a long time. These solutions are such as using wide-area measurement signals [3], phasor measurement unit [4], and flexible AC transmission system. In these solutions, one of the powerful methods for enhancing the transient stability is to use flexible AC transmission system (FACTS) devices [5–8]. Even though the prime objective of shunt FACTS devices (SVC, STATCOM) is to maintain bus voltage by absorbing (or injecting) reactive power, they are also competent of improving the system stability by diminishing (or increasing) the capability of power transfer when the machine angle decreases (increases), which is accomplished by operating the shunt FACTS devices in inductive (capacitive) mode.In many research papers [2, 9–11], the different types of these devices with different control techniques are used for improving transient stability. In between these FACTS devices, the STATCOM is valuable for enhancement power system dynamic stability and frequency stabilization due to the more rapid output response, lower harmonics, superior control stability and small size, and so forth [7, 12]. By their inverter configuration, basic type of STATCOM topology can be realized by either a current-source converter (CSC) or a voltage-source converter (VSC) [13–17]. But recent research confirms several merits of CSC based STATCOM over VSC based STATCOM [18, 19]. These advantages are high converter reliability, quick starting, and inherent short-circuit protection, and the output current of the converter is directly controlled and in low switching frequency this reduces the filtering requirements compared with the case of a VSC. Therefore CSC based STATCOM is very useful in power systems rather than VSC based STATCOM in many cases.Presently the most used techniques for controller design of FACTS devices are the Proportional Integration (PI), PID controller [20], pole-shifting controller, and linear quadratic regulator (LQR) [21]. But LQR and pole-shifting algorithms give quicker response in comparison to PI and PID algorithm [22]. LQR controller gain (K) can be calculated by solving the Riccati equation and K is also dependent on the two cost functions (Q,R). So Riccati equation solvers have some limitations, which relate to the input arguments. But pole shifting method does not face this type of any problem. So pole shifting method gives a better and robust performance in comparison to other methods.The main contribution of this paper is the application of proposed pole-shifting controller based CSC-STATCOM for improvement of power system stability in terms of transient stability by injecting (or absorbing) reactive power. In this paper, the proposed scheme is used in the multimachine power transmission system with dynamic loads under a grievous disturbance condition (three-phase fault or heavy loading) to enhancement of power system transient stability studies and to observe the impact of the CSC based STATCOM on electromechanical oscillations and transmission capacity. Furthermore, the obtained outcomes from the proposed algorithm based CSC-STATCOM are compared to the obtained outcomes from the other shunt FACTS devices (SVC and VSC-STATCOM) which are used in previous works [23, 24].The rest of the paper is prepared as follows. Section2 discusses the circuit modeling and proposed pole-shifting controller designing for CSC based STATCOM. A two-area power system is described with a CSC-STATCOM device in Section 3. Simulation results of the test system with and without CSC based STATCOM for severe contingency are shown in Section 4, to improve the transient stability of the multimachine power system. Comparison among different shunt FACTS devices (SVC, VSC-STATCOM, and CSC-STATCOM) is also described in Section 4. Finally, Section 5 concludes this paper.


## 2. Mathematical Modeling of Pole-Shifting Controller Based CSC-STATCOM

### 2.1. CSC Based STATCOM Model

In this section, to verify the response of the STATCOM on dynamic performance, the mathematical modeling and control strategy of a CSC based STATCOM are needed to be presented. So in the designing of controller for CSC based STATCOM, the state space equations from the CSC-STATCOM circuit must be introduced. To minimize the complexity of mathematical calculation, the theory ofdq transformation of currents has been applied in this circuit, which makes the d and q components independent parameters. Figure 1 shows the circuit diagram of a typical CSC based STATCOM.Figure 1
The representation of CSC based STATCOM.iSR, iSS, iST: line current; vCR, vCS, vCT: voltages across the filter capacitors; vR, vS, vT: line voltages; Idc: dc-side current; Rdc: converter switching and conduction losses; Ldc: smoothing inductor; C: filter capacitance; L: inductance of the line reactor; R: resistance of the line reactor.The basic mathematical equations of the CSC based STATCOM have been derived in the literature [19]. Therefore, only brief details of the primary equations for CSC-STATCOM are given here for the readers’ convenience. Based on the equivalent circuit of CSC-STATCOM as shown in Figure 1, the differential equations for the system can be achieved, which are derived in the abc frame and then transformed into the synchronous dq frame using dq transformation method [25]:
(1)ddtIdc=-RdcLdcIdc-32LdcMdVd-32LdcMqVq,(2)ddtId=-RLId+ωIq-1LEdn+1LVd,(3)ddtIq=-ωId-RLIq+1LVq,(4)ddtVd=-1CId+ωVq+1CMdIdc,(5)ddtVq=-1CIq-ωVd+1CMqIdc.In above differential equationsMd and Mq are the two input variables. Two output variables are Idc and Iq. Here, ω is the rotation frequency of the system and this is equal to the nominal frequency of the system voltage. Id and Iq are the d-axis and q-axis components of the line current. Md and Mq are d-axis and q-axis components of the system modulating signal, respectively. Ed and Eq are direct and quadrature axis of system voltage. Here Eq is taken as a zero. Vd and Vq are the d-axis and q-axis components of the voltage across filter capacitor, respectively.Equation (1) shows that controller for CSC based STATCOM has nonlinearity characteristic. So this nonlinear property can be removed by accurately modeling of CSC based STATCOM. From (1) to (5), it can be seen that nonlinear property in the CSC-STATCOM model is due to the part of Idc. This nonlinear property is removed with the help of active power balance equation. Here, it has been assumed that the power loss in the switches and resistance Rdc is ignored and the turns ratio of the shunt transformer isn : 1. Here it is considered that generated power from the AC-source side (Pac) is equal to absorbed power by DC-source side (Pdc) and power loss at this time is negligible. So this can be written in equation form as
(6)Pac=Pdc⟹-3Ed2nId=LdcIdcddtIdc+RdcIdc2.After using power balance equation and mathematical calculation, nonlinear characteristic is removed from (1). Finally the equation is obtained below:
(7)ddtIdc2=-2RdcLdcIdc2-3EdLdcnId.In (7) state variable (Idc) is replaced by the state variable (Idc2), to make the dynamic equation linear. Finally the resulting better dynamic and robust model of the STATCOM in matrix form can be derived as
(8)ddtidc2idiqvcdvcq=-2RdcLdc3EdLdcn0000-RLωo1L00-ωo-RL01L0-1C00ωo00-1C-ωo0*idc2idiqvcdvcq+0000001C001C*IidIiq+0-1L000*Ed.Above modeling of CSC based STATCOM is written in the form of modern control methods, that is, state-space representation. For state-space modeling of the system, next section has been considered.


### 2.2. Pole-Shifting Controller Design

The pole-shifting technique is one of the basic control methods which are employed in feedback control system theory. Theoretically pole-shifting technique is to set the preferred pole position and to move the pole position of the system to that preferred pole position, to get the desired system outcomes [26]. Here poles of system are shifted because the position of the poles related directly to the eigenvalues of the system, which control the dynamic characteristics of the system outcomes. But for this method, the system must be controllable. In the dynamic modeling of systems, state-space equations involve three types of variables: state variables (x) and input (u) and output (y) variables with disturbance (e). So comparing (8) with the standard state-space representation, that is,
(9)x˙=Ax+Bu+Fe,y=Cx.We get the system matrices as(10)x=Idc2idiqVdVqT;u=IidIiqT;e=Ed;y=Idc2IqT,A=-2RdcLdc-3EdLdc0000-RLω1L00-ω-RL01L0-1c00ω00-1c-ω0;B=0000001c001c;C=1000010000T;F=0-1L000.In the above equations (9) five system states, two control inputs, and two control outputs are presented, where x is the state vector, u is the input vector, A is the basis matrix, B is the input matrix, and e is disturbance input.If the controller is set as(11)u=-K*x+J*yref+N*e,then the state equation of closed loop can be written as(12)x˙=A-B*K*x+J*yref+B*N*e+F*e.Here, for steady state condition:(13)x˙=0.Then,(14)y=C*A-B*K-1*B*J*yref+B*N*e+F*e,
where J=C*-A-B*K-1*B-1 and N=C*-A+B*K-1*B-1*(C*-A+B*K-1*F); these constant values are found out from a mathematical calculation for tracking the reference output value (yref) by the system output value (y). Here K is the state-feedback gain matrix. The gain matrix K is designed in such a way that (15) is satisfied with the desired poles:
(15)sI-(A-BK)=s-P1s-P2⋯s-Pn,
where P1,P2,…,Pn are the desired pole locations. Equation (15) is the desired characteristic polynomial equation. The values of P1,P2,…,Pn are selected such that the system becomes stable and all closed-loop eigenvalues are located in the left half of the complex-plane. The final configuration of the proposed pole-shifting controller based CSC-STATCOM is shown in Figure 2.Figure 2
Control structure of pole-shifting controller based CSC-STATCOM.All the required parameters for the structure of pole-shifting controller based CSC-STATCOM are given in the appendix.


## 2.1. CSC Based STATCOM Model

In this section, to verify the response of the STATCOM on dynamic performance, the mathematical modeling and control strategy of a CSC based STATCOM are needed to be presented. So in the designing of controller for CSC based STATCOM, the state space equations from the CSC-STATCOM circuit must be introduced. To minimize the complexity of mathematical calculation, the theory ofdq transformation of currents has been applied in this circuit, which makes the d and q components independent parameters. Figure 1 shows the circuit diagram of a typical CSC based STATCOM.Figure 1
The representation of CSC based STATCOM.iSR, iSS, iST: line current; vCR, vCS, vCT: voltages across the filter capacitors; vR, vS, vT: line voltages; Idc: dc-side current; Rdc: converter switching and conduction losses; Ldc: smoothing inductor; C: filter capacitance; L: inductance of the line reactor; R: resistance of the line reactor.The basic mathematical equations of the CSC based STATCOM have been derived in the literature [19]. Therefore, only brief details of the primary equations for CSC-STATCOM are given here for the readers’ convenience. Based on the equivalent circuit of CSC-STATCOM as shown in Figure 1, the differential equations for the system can be achieved, which are derived in the abc frame and then transformed into the synchronous dq frame using dq transformation method [25]:
(1)ddtIdc=-RdcLdcIdc-32LdcMdVd-32LdcMqVq,(2)ddtId=-RLId+ωIq-1LEdn+1LVd,(3)ddtIq=-ωId-RLIq+1LVq,(4)ddtVd=-1CId+ωVq+1CMdIdc,(5)ddtVq=-1CIq-ωVd+1CMqIdc.In above differential equationsMd and Mq are the two input variables. Two output variables are Idc and Iq. Here, ω is the rotation frequency of the system and this is equal to the nominal frequency of the system voltage. Id and Iq are the d-axis and q-axis components of the line current. Md and Mq are d-axis and q-axis components of the system modulating signal, respectively. Ed and Eq are direct and quadrature axis of system voltage. Here Eq is taken as a zero. Vd and Vq are the d-axis and q-axis components of the voltage across filter capacitor, respectively.Equation (1) shows that controller for CSC based STATCOM has nonlinearity characteristic. So this nonlinear property can be removed by accurately modeling of CSC based STATCOM. From (1) to (5), it can be seen that nonlinear property in the CSC-STATCOM model is due to the part of Idc. This nonlinear property is removed with the help of active power balance equation. Here, it has been assumed that the power loss in the switches and resistance Rdc is ignored and the turns ratio of the shunt transformer isn : 1. Here it is considered that generated power from the AC-source side (Pac) is equal to absorbed power by DC-source side (Pdc) and power loss at this time is negligible. So this can be written in equation form as
(6)Pac=Pdc⟹-3Ed2nId=LdcIdcddtIdc+RdcIdc2.After using power balance equation and mathematical calculation, nonlinear characteristic is removed from (1). Finally the equation is obtained below:
(7)ddtIdc2=-2RdcLdcIdc2-3EdLdcnId.In (7) state variable (Idc) is replaced by the state variable (Idc2), to make the dynamic equation linear. Finally the resulting better dynamic and robust model of the STATCOM in matrix form can be derived as
(8)ddtidc2idiqvcdvcq=-2RdcLdc3EdLdcn0000-RLωo1L00-ωo-RL01L0-1C00ωo00-1C-ωo0*idc2idiqvcdvcq+0000001C001C*IidIiq+0-1L000*Ed.Above modeling of CSC based STATCOM is written in the form of modern control methods, that is, state-space representation. For state-space modeling of the system, next section has been considered.


## 2.2. Pole-Shifting Controller Design

The pole-shifting technique is one of the basic control methods which are employed in feedback control system theory. Theoretically pole-shifting technique is to set the preferred pole position and to move the pole position of the system to that preferred pole position, to get the desired system outcomes [26]. Here poles of system are shifted because the position of the poles related directly to the eigenvalues of the system, which control the dynamic characteristics of the system outcomes. But for this method, the system must be controllable. In the dynamic modeling of systems, state-space equations involve three types of variables: state variables (x) and input (u) and output (y) variables with disturbance (e). So comparing (8) with the standard state-space representation, that is,
(9)x˙=Ax+Bu+Fe,y=Cx.We get the system matrices as(10)x=Idc2idiqVdVqT;u=IidIiqT;e=Ed;y=Idc2IqT,A=-2RdcLdc-3EdLdc0000-RLω1L00-ω-RL01L0-1c00ω00-1c-ω0;B=0000001c001c;C=1000010000T;F=0-1L000.In the above equations (9) five system states, two control inputs, and two control outputs are presented, where x is the state vector, u is the input vector, A is the basis matrix, B is the input matrix, and e is disturbance input.If the controller is set as(11)u=-K*x+J*yref+N*e,then the state equation of closed loop can be written as(12)x˙=A-B*K*x+J*yref+B*N*e+F*e.Here, for steady state condition:(13)x˙=0.Then,(14)y=C*A-B*K-1*B*J*yref+B*N*e+F*e,
where J=C*-A-B*K-1*B-1 and N=C*-A+B*K-1*B-1*(C*-A+B*K-1*F); these constant values are found out from a mathematical calculation for tracking the reference output value (yref) by the system output value (y). Here K is the state-feedback gain matrix. The gain matrix K is designed in such a way that (15) is satisfied with the desired poles:
(15)sI-(A-BK)=s-P1s-P2⋯s-Pn,
where P1,P2,…,Pn are the desired pole locations. Equation (15) is the desired characteristic polynomial equation. The values of P1,P2,…,Pn are selected such that the system becomes stable and all closed-loop eigenvalues are located in the left half of the complex-plane. The final configuration of the proposed pole-shifting controller based CSC-STATCOM is shown in Figure 2.Figure 2
Control structure of pole-shifting controller based CSC-STATCOM.All the required parameters for the structure of pole-shifting controller based CSC-STATCOM are given in the appendix.


## 3. Two-Area Power System with CSC-STATCOM Facts Device

In this section, consider a two-area power system with a CSC-STATCOM at busb is connected through a long transmission system, where CSC-STATCOM is used as a shunt current source device. Figure 3 shows this type of representation. The dynamic model of the machine, with a CSC-STATCOM, can be written in the differential algebraic equation form as follows:
(16)δ˙=ω,(17)ω˙=1MPm-Peo-Pecsc.Figure 3
A single line diagram of two-area power system with CSC-STATCOM.Hereω is the rotor speed, δ is the rotor angle, Pm is the mechanical input power of generator, the output electrical power without CSC-STATCOM is represented by Peo, and M is the moment of inertia of the rotor. Equation (17) is also called swing equation. The additional factor of the output electrical power of generator from a CSC-STATCOM is Pecsc in the swing equation.Here for calculation ofPecsc, we assume the CSC-STATCOM works in capacitive mode. Then the injected current from CSC-STATCOM to test system can be written as
(18)Icsc=Icsc∠θbo-90°,
where θbo is the voltage angle at bus b in absentia of CSC-STATCOM. In Figure 3, the magnitude (Ed) and angle (θb) of voltage at bus b can be computed as
(19)θb=tan-1X2E1sinδX2E1cos⁡δ+X1E2,(20)Ed=X2E1cos⁡(δ-θbo)+X1E2cos⁡θboX1+X2+X1X2X1+X2Icsc.From (20), it can be said that the voltage magnitude of bus b (Ed) depends on the STATCOM current (Icsc). In (17), the electrical output power Pecsc of machine due to a CSC-STATCOM can be expressed as
(21)Pecsc=E1EdX1sinδ-θb.Finally, using (20) and (21) the total electrical output (Pe) of machine with CSC-STATCOM can be written as
(22)Pe=Peo+Pecsc⟹Pe=Peo+X1X2E1X1+X2X1Icscsinδ-θb.All above equations are represented for the capacitive mode of CSC-STATCOM. So if the inductive mode of operation is needed in the system thenIcsc is replaced by -Icsc in (18), (20), and (22). With the help of (17), the power-angle curve of the test system can be drawn for stability analysis as shown in Figure 4.Figure 4
Power-angle characteristic of the test system with a CSC-STATCOM.The power-angle (P-δ) curve of the test system without a CSC-STATCOM is represented by curve A (also called prefault condition) in Figure 4. Here the mechanical input is Pm, electrical output is Pe, and initial angle is δ0. When a fault occurs, Pe suddenly decreases and the operation shifts from point a to point b at curve B and thus the machine starts accelerating from point b to point c, where accelerating power   [Pa=(Pm-Pe)]>0. At fault clearing, Pe suddenly increases and the area a-b-c-d-a represents the accelerating area Aa as defined in (23). If the CSC-STATCOM operates in a capacitive mode (at fault clearing), Pe increases to point e at curve C (also called postfault condition). At this time Pa is negative. Thus the machine starts decelerating but its angle continues to increase from point e to the point f until it reaches a maximum allowable value δm at point f, for system stability. The area e-f-g-d-e represents the decelerating area Ad as defined in (23). From previous literature [1], equal area criterion for stability of the system can be written as
(23)∫δ0δcPm-Pefdδ=∫δcδmPep-Pmdδ⟹Aa=Ad.This equation is generated from Figure4, where δc is critical clearing angle. Pep is an electrical output for postfault condition. Pef is an electrical output during fault condition. From Figure 4 we clearly see, for capacitive mode operation (Icsc>0) of CSC-STATCOM device, the P-δ curve is not only uplifted but also displaced toward right and that endues more decelerating area and hence higher stability limit. After that Section 4 is considered, to verify the impact of the proposed controller based CSC-STATCOM modeling in the analysis of the multimachine power system stability in terms of transient stability studies.


## 4. Simulation Results

### 4.1. Multimachine Power System under Study

In this section, two-area four-machine power system is considered for the test system study. For this type of test system, a 500 kV transmission system with four hydraulic power plants (G1, G2, G3, and G4) connected through 500-km long transmission line is taken as shown in Figure5. Here combination of machine G1 and machine G2 represents the area 1 and combination of machine G3 and machine G4 represents the area 2. A rating of first (G1) and second (G2) power generation plant is 13.8 kv/1000 MVA. The electrical outputs of the third (G3) and fourth (G4) power plant are 3000 MVA. One 6000 MW large resistive load is connected at the bus B3 and another 1000 MW resistive load is connected at the bus B1 as shown in Figure 5. In this figure, 828 MW power is flowing through bus B2, 944 MW and 877 MW power are flowing through buses B4 and B5, respectively, and 2487 MW and 2643 MW power are flowing through buses B6 and B7, respectively. These values are based on power flow calculation. To improve the transient stability of the test-system after disturbances (faults or heavy loading), a pole-shifting controller based CSC-STATCOM is connected at the mid-point of transmission line (at bus B2). In order to previous work [27], shunt FACTS devices give maximum efficiency when connected at the mid-point of transmission line. Here each hydraulic generating unit is assembled with a turbine-governor set, excitation system, and synchronous generator. But the hydraulic generating unit is not explained in this paper. This is explained in [1]. All the data required for this test system model are given in the appendix.Figure 5
The single line diagram of the test-system model for transient stability study of multimachine power system.Now a severe contingency will be applied to the test system and to observe the impact of the pole-shifting controller based CSC-STATCOM for maintaining the system stability through MATLAB/SIMULINK.


### 4.2. Case I—Short-Circuit Fault

In this case, it is considered that a three-phase fault is occurring at near bus B1 att=0.1 s and is cleared at 0.277 s. The impact of system with and without CSC based STATCOM to this disturbance is shown in Figures 6, 7, 8, 9, and 10. Here simulation is carried out for 9 s to observe the nature of transients. It is clear that the system without CSC-STATCOM is unstable upon the clearance of the fault from Figures 6 to 8. But this system with pole-shifting controller based CSC-STATCOM is restored and stable upon the clearance of the fault and settled down after about 3 to 4 s in Figures 7 to 10. Synchronism in between four machine system is also maintained in these figures. In Figure 7,  δ1, δ2, δ3, and δ4 are the rotor angles of machines G1, G2, G3, and G4, respectively. Figure 7 represents the variation of rotor angles differences of the multimachine system. In Figure 8, speed of machines G1, G2, G3, and G4 is represented by ω1, ω2, ω3, and ω4, respectively. Here the critical clearing time (CCT) of fault is also found out for the test system stability by simulation. The fault time is increased in simulation to find the critical stability margin; thus CCT is obtained. CCT is defined as the maximal fault duration for which the system remains transiently stable [1]. CCT of the fault for the system with and without CSC-STATCOM is 325 ms and 276 ms, respectively, as shown in Table 1. So it shows that CCT of fault is also increased due to the impact of pole-shifting controller based CSC-STATCOM. Clearly, waveforms show that proposed topology is more effective and robust performance than that of the without pole-shifting controller based CSC-STATCOM in terms of settling time, CCT, and transient stability of the test-system.Table 1
CCT of disturbances for the system stability with different topologies.
Case number
System with different topologies
Critical clearing time (CCT)
1
Without CSC-STATCOM
276 ms
With CSC-STATCOM
325 ms
2
Without CSC-STATCOM
465 ms
With CSC-STATCOM
764 msSystem response without CSC-STATCOM for a three-phase fault (case I).
(a)
Positive sequence voltages at different buses B1, B2, and B3
(b)
Power flow at bus B2Variation of rotor angles differences of the test system for case I.
(a)
(b)
(c)Speed differences variation of the test system for case I.
(a)
Speed difference variation of machines G1 and G4
(b)
Speed difference variation of machines G3 and G4Test system response with pole-shifting controller based CSC-STATCOM for a three-phase fault (case I).
(a)
Positive sequence voltages at different buses B1, B2, and B3
(b)
Power flow at bus B2For case 1 (a) terminal voltagesVt1, Vt2, Vt3, and Vt4 of machine generators G1, G2, G3, and G4 with CSC-STATCOM, respectively, (b) reactive power inject (or absorb) by pole-shifting controller based CSC-STATCOM.
(a)
(b)


### 4.3. Case II—Large Loading

For heavy loading case, one large series load centre (200 MW/6000 Mvar) is connected at near bus B1 (i.e., at near machine G2) in Figure5. This large loading is occurring only at time period 0.1 s to 0.466 s. Due to this disturbance, the simulation results of test system with and without CSC-STATCOM are shown in Figures 11 to 15. Clearly, the system is unstable in the absence of the pole-shifting controller based CSC-STATCOM device due to this disturbance in Figures 11 to 13 and 15. But system with pole-shifting controller based CSC-STATCOM is continuing stable condition in Figures 12 to 15. In Figure 12,  δ1, δ2, δ3, and δ4 are the rotor angles of machines G1, G2, G3, and G4, respectively. Figure 12 represents the variation of rotor angles differences of the multimachine system. In Figure 13, speeds of machines G1, G2, G3, and G4 are represented by ω1, ω2, ω3, and ω4, respectively. System voltages at different buses B1, B2, and B3 with proposed scheme are shown in Figure 14. Here CCT for the system with and without CSC-STATCOM is 764 ms and 465 ms, respectively, which are provided in Table 1. Clearly shown, CCT for test system is better due to the impact of pole-shifting controller based CSC-STATCOM. Hence the performance of the proposed scheme is still satisfactory in this case.Test-system response without CSC-STATCOM with a heavy loading (case II).
(a)
Positive sequence voltages at different buses B1, B2, and B3
(b)
Power flow at bus B2Variation of rotor angles differences of the test system in case II.
(a)
(b)
(c)Speed differences variation of the multimachine system in case II.
(a)
Speed difference variation of machines G1 and G4
(b)
Speed difference variation of machines G3 and G4Test system response with pole-shifting controller based CSC-STATCOM for a heavy loading (case II).
(a)
Positive sequence voltages at different buses B1, B2, and B3
(b)
Power flow at bus B2Speed versus rotor angle graph for (a) machine G1, (b) machine G2.
(a)
(b)


### 4.4. Case III—A Comparative Study

In order to show the robustness performance of the proposed scheme for improving the transient stability of the system, after this obtained outcomes from the proposed pole-shifting controller based CSC-STATCOM are compared to the obtained outcomes from the other shunt FACTS devices (SVC and VSC based STATCOM) which have been used in some most cited papers [23, 24]. In this case, it has been assumed that the test system is same for all shunt FACTS devices and all shunt FACTS devices have the same rating (200 Mvar). After this, we check the impact of these shunt FACTS devices on the test system with three-phase fault condition. The obtained simulation results with these shunt FACTS devices are compared in Figure 16. In this figure system with SVC is reached in unstable condition. Here three-phase fault duration is 0.1 s to 0.285 s. CCT for the system with different shunt FACTS devices is shown in Table 2. If heavy loading condition is applied during 0.1 s to 0.486 s in this case, then effect of this disturbance is shown in Figure 17. In all figures of Figures 16 and 17, multimachine power system with proposed topology based CSC-STATCOM is continuing stable condition. In these figures, the oscillations are damped more quickly and settled down after 3 to 4 sec by proposed topology based CSC-STATCOM in comparison to other devices. In this comparative study, transient stability margin is also increased by increasing the CCT with the help of CSC-STATCOM device from Table 2. Clearly, waveforms show that CSC-STATCOM is more effective and robust performance than that of other shunt FACTS devices (SVC and VSC-STATCOM) in terms of oscillation damping, settling time, CCT, and transient stability of the multimachine power system in Figures 16 and 17.Table 2
CCT of disturbances for the test system stability with different shunt FACTS devices.
S. number
Type of disturbances
FACTS devices
Critical clearing time (CCT)
1
Fault condition
With SVC
285 ms
With VSC-STATCOM
288 ms
With CSC-STATCOM
325 ms
2
Heavy loading condition
With SVC
486 ms
With VSC-STATCOM
489 ms
With CSC-STATCOM
764 msSystem response with different shunt FACTS devices (SVC, CSC-STATCOM, and VSC-STATCOM) for fault condition (case III).
(a)
Speed difference variation of machines G1 and G4
(b)
Speed difference variation of machines G3 and G4
(c)
Variation of rotor angles difference of machines G1 and G4
(d)
Variation of rotor angles difference of machines G3 and G4Figure 17
System response with different shunt FACTS devices (SVC, CSC-STATCOM, and VSC-STATCOM) in a heavy loading condition (case III).


## 4.1. Multimachine Power System under Study

In this section, two-area four-machine power system is considered for the test system study. For this type of test system, a 500 kV transmission system with four hydraulic power plants (G1, G2, G3, and G4) connected through 500-km long transmission line is taken as shown in Figure5. Here combination of machine G1 and machine G2 represents the area 1 and combination of machine G3 and machine G4 represents the area 2. A rating of first (G1) and second (G2) power generation plant is 13.8 kv/1000 MVA. The electrical outputs of the third (G3) and fourth (G4) power plant are 3000 MVA. One 6000 MW large resistive load is connected at the bus B3 and another 1000 MW resistive load is connected at the bus B1 as shown in Figure 5. In this figure, 828 MW power is flowing through bus B2, 944 MW and 877 MW power are flowing through buses B4 and B5, respectively, and 2487 MW and 2643 MW power are flowing through buses B6 and B7, respectively. These values are based on power flow calculation. To improve the transient stability of the test-system after disturbances (faults or heavy loading), a pole-shifting controller based CSC-STATCOM is connected at the mid-point of transmission line (at bus B2). In order to previous work [27], shunt FACTS devices give maximum efficiency when connected at the mid-point of transmission line. Here each hydraulic generating unit is assembled with a turbine-governor set, excitation system, and synchronous generator. But the hydraulic generating unit is not explained in this paper. This is explained in [1]. All the data required for this test system model are given in the appendix.Figure 5
The single line diagram of the test-system model for transient stability study of multimachine power system.Now a severe contingency will be applied to the test system and to observe the impact of the pole-shifting controller based CSC-STATCOM for maintaining the system stability through MATLAB/SIMULINK.


## 4.2. Case I—Short-Circuit Fault

In this case, it is considered that a three-phase fault is occurring at near bus B1 att=0.1 s and is cleared at 0.277 s. The impact of system with and without CSC based STATCOM to this disturbance is shown in Figures 6, 7, 8, 9, and 10. Here simulation is carried out for 9 s to observe the nature of transients. It is clear that the system without CSC-STATCOM is unstable upon the clearance of the fault from Figures 6 to 8. But this system with pole-shifting controller based CSC-STATCOM is restored and stable upon the clearance of the fault and settled down after about 3 to 4 s in Figures 7 to 10. Synchronism in between four machine system is also maintained in these figures. In Figure 7,  δ1, δ2, δ3, and δ4 are the rotor angles of machines G1, G2, G3, and G4, respectively. Figure 7 represents the variation of rotor angles differences of the multimachine system. In Figure 8, speed of machines G1, G2, G3, and G4 is represented by ω1, ω2, ω3, and ω4, respectively. Here the critical clearing time (CCT) of fault is also found out for the test system stability by simulation. The fault time is increased in simulation to find the critical stability margin; thus CCT is obtained. CCT is defined as the maximal fault duration for which the system remains transiently stable [1]. CCT of the fault for the system with and without CSC-STATCOM is 325 ms and 276 ms, respectively, as shown in Table 1. So it shows that CCT of fault is also increased due to the impact of pole-shifting controller based CSC-STATCOM. Clearly, waveforms show that proposed topology is more effective and robust performance than that of the without pole-shifting controller based CSC-STATCOM in terms of settling time, CCT, and transient stability of the test-system.Table 1
CCT of disturbances for the system stability with different topologies.
Case number
System with different topologies
Critical clearing time (CCT)
1
Without CSC-STATCOM
276 ms
With CSC-STATCOM
325 ms
2
Without CSC-STATCOM
465 ms
With CSC-STATCOM
764 msSystem response without CSC-STATCOM for a three-phase fault (case I).
(a)
Positive sequence voltages at different buses B1, B2, and B3
(b)
Power flow at bus B2Variation of rotor angles differences of the test system for case I.
(a)
(b)
(c)Speed differences variation of the test system for case I.
(a)
Speed difference variation of machines G1 and G4
(b)
Speed difference variation of machines G3 and G4Test system response with pole-shifting controller based CSC-STATCOM for a three-phase fault (case I).
(a)
Positive sequence voltages at different buses B1, B2, and B3
(b)
Power flow at bus B2For case 1 (a) terminal voltagesVt1, Vt2, Vt3, and Vt4 of machine generators G1, G2, G3, and G4 with CSC-STATCOM, respectively, (b) reactive power inject (or absorb) by pole-shifting controller based CSC-STATCOM.
(a)
(b)


## 4.3. Case II—Large Loading

For heavy loading case, one large series load centre (200 MW/6000 Mvar) is connected at near bus B1 (i.e., at near machine G2) in Figure5. This large loading is occurring only at time period 0.1 s to 0.466 s. Due to this disturbance, the simulation results of test system with and without CSC-STATCOM are shown in Figures 11 to 15. Clearly, the system is unstable in the absence of the pole-shifting controller based CSC-STATCOM device due to this disturbance in Figures 11 to 13 and 15. But system with pole-shifting controller based CSC-STATCOM is continuing stable condition in Figures 12 to 15. In Figure 12,  δ1, δ2, δ3, and δ4 are the rotor angles of machines G1, G2, G3, and G4, respectively. Figure 12 represents the variation of rotor angles differences of the multimachine system. In Figure 13, speeds of machines G1, G2, G3, and G4 are represented by ω1, ω2, ω3, and ω4, respectively. System voltages at different buses B1, B2, and B3 with proposed scheme are shown in Figure 14. Here CCT for the system with and without CSC-STATCOM is 764 ms and 465 ms, respectively, which are provided in Table 1. Clearly shown, CCT for test system is better due to the impact of pole-shifting controller based CSC-STATCOM. Hence the performance of the proposed scheme is still satisfactory in this case.Test-system response without CSC-STATCOM with a heavy loading (case II).
(a)
Positive sequence voltages at different buses B1, B2, and B3
(b)
Power flow at bus B2Variation of rotor angles differences of the test system in case II.
(a)
(b)
(c)Speed differences variation of the multimachine system in case II.
(a)
Speed difference variation of machines G1 and G4
(b)
Speed difference variation of machines G3 and G4Test system response with pole-shifting controller based CSC-STATCOM for a heavy loading (case II).
(a)
Positive sequence voltages at different buses B1, B2, and B3
(b)
Power flow at bus B2Speed versus rotor angle graph for (a) machine G1, (b) machine G2.
(a)
(b)


## 4.4. Case III—A Comparative Study

In order to show the robustness performance of the proposed scheme for improving the transient stability of the system, after this obtained outcomes from the proposed pole-shifting controller based CSC-STATCOM are compared to the obtained outcomes from the other shunt FACTS devices (SVC and VSC based STATCOM) which have been used in some most cited papers [23, 24]. In this case, it has been assumed that the test system is same for all shunt FACTS devices and all shunt FACTS devices have the same rating (200 Mvar). After this, we check the impact of these shunt FACTS devices on the test system with three-phase fault condition. The obtained simulation results with these shunt FACTS devices are compared in Figure 16. In this figure system with SVC is reached in unstable condition. Here three-phase fault duration is 0.1 s to 0.285 s. CCT for the system with different shunt FACTS devices is shown in Table 2. If heavy loading condition is applied during 0.1 s to 0.486 s in this case, then effect of this disturbance is shown in Figure 17. In all figures of Figures 16 and 17, multimachine power system with proposed topology based CSC-STATCOM is continuing stable condition. In these figures, the oscillations are damped more quickly and settled down after 3 to 4 sec by proposed topology based CSC-STATCOM in comparison to other devices. In this comparative study, transient stability margin is also increased by increasing the CCT with the help of CSC-STATCOM device from Table 2. Clearly, waveforms show that CSC-STATCOM is more effective and robust performance than that of other shunt FACTS devices (SVC and VSC-STATCOM) in terms of oscillation damping, settling time, CCT, and transient stability of the multimachine power system in Figures 16 and 17.Table 2
CCT of disturbances for the test system stability with different shunt FACTS devices.
S. number
Type of disturbances
FACTS devices
Critical clearing time (CCT)
1
Fault condition
With SVC
285 ms
With VSC-STATCOM
288 ms
With CSC-STATCOM
325 ms
2
Heavy loading condition
With SVC
486 ms
With VSC-STATCOM
489 ms
With CSC-STATCOM
764 msSystem response with different shunt FACTS devices (SVC, CSC-STATCOM, and VSC-STATCOM) for fault condition (case III).
(a)
Speed difference variation of machines G1 and G4
(b)
Speed difference variation of machines G3 and G4
(c)
Variation of rotor angles difference of machines G1 and G4
(d)
Variation of rotor angles difference of machines G3 and G4Figure 17
System response with different shunt FACTS devices (SVC, CSC-STATCOM, and VSC-STATCOM) in a heavy loading condition (case III).


## 5. Conclusions

In this paper, the dynamic modeling of CSC based STATCOM is studied and proposed pole-shifting controller for the best input-output response of CSC-STATCOM is presented in order to enhance the system transient stability of the multimachine power system with the different disturbances. The novelty in proposed approach lies in the fact that transient stability of a two-area four-machine power system is improved and critical clearing time of the disturbance is also increased. In this work, the proposed scheme is verified from MATLAB package. It has been also observed that pole-shifting controller based CSC-STATCOM can be more reliable and very effective rather than other shunt FACTS devices (SVC and VSC-STATCOM) in terms of oscillation damping, CCT, and transient stability of a multimachine power system. So, finally it can be said that CSC based STATCOM can be regarded as an alternative FACTS device to that of other shunt FACTS devices.



---

*Source: 626731-2015-02-04.xml*