# Improving Power System Stability with Gramian Matrix-Based Optimal Setting of a Single Series FACTS Device: Feasibility Study in Vietnamese Power System

Authors: Le Van Dai; Doan Duc Tung; Tran Le Thang Dong; Le Cao Quyen
Journal: Complexity
 (2017)

---

## Abstract

The Vietnamese power system has experienced instabilities due to the effect of increase in peak load demand or contingency grid faults; hence, using flexible alternating-current transmission systems (FACTS) devices is a best choice for improving the stability margins. Among the FACTS devices, the thyristor-controlled series capacitor (TCSC) is a series connected FACTS device widely used in power systems. However, in practice, its influence and ability depend on setting. For solving the problem, this paper proposes a relevant method for optimal setting of a single TCSC for the purpose of damping the power system oscillations. This proposed method is developed from the combination between the energy method and Hankel-norm approximation approach based on the controllability Gramian matrix considering the Lyapunov equation to search for a number of feasible locations on the small-signal stability analysis. The transient stability analysis is used to compare and determine appropriate settings through various simulation cases. The effectiveness of the proposed method is confirmed by the simulation results based on the power system simulation engineering (PSS/E) and MATLAB programs. The obtained results show that the proposed method can apply to immediately solve the difficulties encountering in the Vietnamese power system.


---

## Body

## 1. Introduction

Motivation. When the social, political, and technological aspects develop, the demand of electric power grows rapidly, resulting in the increase in the scale and complexity of power systems. Some characteristics, such as long transmission distances over weak grids, highly variable generation patterns, and heavy load, tend to increase the wide-area electromechanical oscillations. These oscillations threaten the secure operation of the power systems and if not controlled efficiently can lead to generator outages, line tripping, and even large-scale blackouts.During over several decades, a number of electrical power systems in the world have been faced with the serious power system blackouts, such as Indian on July 30 and 31, 2012 [1]; Dubai on June 9, 2005; Malaysia on Jan. 13, 2005; Kuwait on Nov. 01, 2004; Australia on Aug. 14, 2004; Shanghai on Aug. 27, 2003; Sweden/Denmark on Sept. 23, 2003; Finland on Aug., 1992; Western France on Jan. 12, 1987; Florida on Jan. 12, 1987; Florida, USA, on May 17, 1985; Belgium on Aug. 4, 1982 [2–5]. Typically, the event is on August 10, 1996, western blackout of Western Systems Coordinating Council (WSCC) interconnection, caused by the negative damping of the 0.25 Hz western interarea mode [6]. As a result, leading to the split of the system into four large islands, over 7.5 million customers experienced outages ranging from a few minutes to nine hours, and total load loss is 30,500 MW because of the poor oscillations damping. Figures 1(a) and 1(b) show the measured and simulated power swings of this event. The measured data show that the standard planning models could be unreliable predictors of oscillatory behavior.Figure 1
The event of the small-signal instability during the blackout WSCC system on Aug. 10, 1996: (a) the observed California Oregon Interconnections power (Dittmer control center); (b) the simulated California Oregon Interconnections power (initial base case).
(a) (b)In recent years, the economic tempo in the Vietnam has been developed rapidly; the total load has increased continuously. Accordingly, the Vietnamese power system has been faced with the serious power system blackouts, such as on Dec. 27, 2006; July 20, 2007; Apr. 09, 2007; and the latest event on May 22, 2013. All of the technical problems that the Vietnamese power system identified after these events had already been reported [7, 8].The power system oscillations occur in the power systems because of the contingencies, such as the grid faults and sudden load changes; the dampening of these oscillations is necessary for a secure system operation. If the controlled systems react quickly against faults, the power system stability will enhance significantly. The advanced power electronics has led to a new design called flexible alternating-current transmission systems (FACTS) by Electrical Power Research Institute (EPRI). The FACTS devices make more use of the exiting capacities in the power system and enhance the power system stability. For example, the parameters in the power system are controlled and the load flow is modified to preclude the overload of transmission lines after the grid faults. The FACTS devices are widely used to improve the efficiency of power system operation. However, the benefits derived from FACTS controllers, such as the small-signal stability and transient stability that depend greatly on their optimal placement in the power systems [9, 10]. Therefore, looking for the optimal placement of FACTS devices in the large-scale power systems is an interesting research topic.State of the Science. The damping of electromechanical oscillations in the power systems is a matter to belong with power engineers and become one of the most interesting research topics [11, 12]. The PSSs are a profound influence solution on improving the power system stability. This solution is based on the excitation system and the generator speed deviation signals, in which the speed deviation is used as the supplementary control signals in order to provide to the automatic voltage regulators. The fast progress in the field of power electronics had contributed to the development of the electrical power industry based on the controllable utilization of FACTS [13]. This device can dampen both the local- and interarea oscillations [10], whereas PSSs could only dampen the local-area oscillations [14]. Thus, FACTS is the best choice to enhance the stability margin of the existing power systems [9, 15] and especially to TCSC. This device is considered as one of the most effective FACTS and is widely used to regulate the power flow in transmission lines, dampen the power oscillations, mitigate subsynchronous resonance (SSR), enhance power system stability, and so forth by changing the reactance of transmission line [16]. There is hard sledding for the researchers regarding how to determine optimal location of FACTS.The methods for solving location problems can be classified into two categories: (i) analytical techniques and (ii) heuristic optimization approaches [17–19]. Among heuristics, Particle Swarm Optimization (PSO) has some advantages as less control parameters, no require preconditions (continuity or differentiability of objective functions), and slow computational burden, which make it popular to solve optimal location and design problems of FACTS [20–24]. However, PSO had drawback as slow convergence in search stage and limit of local search ability; furthermore, the algorithm will not work out the problems of scattering and optimization [15]. Genetic algorithm (GA) has been proposed to solve optimal location and design problems of FACTS [9, 15, 19, 25]. This method has been applied to obtain promising results. However, downside of GA is that the requested run time is very long when studying the large-scale systems [26]. In [24], optimal location of TCSC is determined via PSO whose objective is to maximize small-signal stability. In [9, 23], a mathematical objective function is used to derive objective function for GA and PSO by considering the load-ability in the system security margins, respectively. In [27], the authors used a method based on hybrid between PSO and GA for optimal allocation of TCSC to enhance voltage stability and reduce power system losses. The brainstorm optimization algorithm (BSOA) is new heuristic optimization algorithm proposed by the authors in [28] to seek the optimal location of TCSC for enhancing the voltage profile. The obtained results are better than PSO and GA.For the analytical approaches, the modal controllability index has been developed by the authors in [29] to find suitable location of TCSC for dampening interarea model of oscillations. However, the authors consider the simulation with or without FACTS placed in the power system to calculate the maximum controllability index values corresponding to critical mode. Particularly, their main interest is in the input signals that do not know what is occurring at output of FACTS through observability index values. In [30], Vaidya and Rajderkar used the sensitivity-based method and the authors in [31] applied real power flow performance index sensitivity to determine the optimal location of TCSC to enhance the power system stability. In [32], the eigenvalue sensitivity method is utilized to find optimal location of controllable series capacitors for dampening power system location. The eigenvalue sensitivity values are calculated based on the modal controllability and observability indices of series reactance modulation. In [33], the suitable location of series compensators is determined based on trajectory sensitivity analysis. The objective is to maximize the benefit of series compensators in order to improve the generator rotor angle stability. The energy method based on Gramian matrices is another technique developed by the authors in [34, 35] to determine the optimal setting of TCSC and static VAR compensator (SVC) obtained with promising result. However, the difficult of this method is when calculating the too large-scale power systems.It is observed that most of existing methods in the previous literatures have been proposed recently for the location of FACTS. These methods have several drawbacks; firstly, the computation of critical modes may be questionable in case of large-scale power system since they may not be unique. Moreover, the computation of them also depends on the local or interarea modes. Secondly, the computation of participation factors is only based on the state variables and neglects the input-output behavior. Thirdly, it just focuses on analyzing the small-scale power systems. Therefore, in order to overcome these drawbacks, this paper is a continuation of [34] and combines the Hankel-norm approximation method to determine the best location for installing TCSC with objective for damping power system oscillation of the practical power system, to wit Vietnamese power system. It indicates that using the Gramian matrix-based method to calculate the complexity and large-scale power systems takes a lot of time since the system state matrix is very large. Therefore, the Hankel-norm approximation method [36] is proposed to solve such problem. The selection of the input signal for TCSC controller is an important consideration for seeking the optimal location to dampen the interarea oscillations, in which the line reactive and active power, line current, and bus voltage are all good selections. In this study, the active power in transmission line is considered as an effective input signal for TCSC controller. The contingency cases are considered based on the active power perturbation signals in the transmission lines that were selected on the basis of the real power line flow performance index (PI) introduced in [37]. In addition, in order to determine the parameters of the TCSC controller, the optimal method in [38] is also considered in this paper.Contribution. In this paper, a relevant method for determining the optimal placement of TCSC controller is proposed to enhance the stability of large-scale power systems. This proposed method is developed from the energy approach based on the controllability Gramian matrix of the linearized system. The multimachine power system with TCSC controller is expressed in the form of a differential algebraic equation (DAE) model. The controllability Gramian matrix is obtained from the unique positive definite solutions of the Lyapunov equation and it depends on the control input-output and state matrixes of system. The optimal location is selected based on the maximum total Gramian energy calculated from the contingency outage cases (disturbance in lines) on the small-signal stability analysis that means the control input must insure the smallest control energy. In line with this purpose, the Hankel-norm approximation method is applied to reduce the number of state variables when dealing with large-scale power systems.The main new contributions of this paper are summarized as follows:(i)
To develop a relevant method to determine the optimal location of FACTS on the small-signal stability analysis(ii)
To propose an association between proposed method and the Hankel-norm approximation method to limit the time calculation, so that the proposed method can be easily implemented for complexity and large-scale power systems.The remainder of this paper is organized as follows: Section2 addresses the principle of characteristics of the Vietnamese power system and the differential algebraic equation (DAE) model of the power system with and without TCSC controller. The proposed method for optimal location of TCSC controller based on the Gramian matrices and the Hankel-norm approximation method are introduced in Section 3. The case studies and conclusions are given in Sections 4 and 5, respectively. Finally, the algorithm for the transformation matrix can be found in the Appendix.


## 2. Theoretical Analysis

### 2.1. Characteristics of the Vietnamese Power System

The Vietnamese 500 kV power system operated in 1994. The length of transmission line of about 1,500 km was connected from the HoaBinh to PhuLam power stations with total rating of 2,700 MVA. It reached the length of about 5690 km and total capacity of all the 500 kV substations of 22,800 MVA in 2015. Several areas of 500 kV transmission lines are compensated by shunt reactors of about 70% and series capacitors of about 60%. The total generation capacity of the Vietnamese power system could reach 60,000 MW by the end of 2020 [39]. The result of load flow calculation on the 500 kV power system is shown in Figure 2.Figure 2
The result of load flow calculation on the Vietnamese 500/220 kV power system 2020.


### 2.2. DAE Model of Power System

The methodological method for dynamic modeling of generalm-machine and n-bus power system has been described in [40] and is applied for this work. In this model, each synchronous generator is represented by two-axis flux decay dynamic model along with IEEE type I. The differential algebraic equation (DAE) model of the power system without TCSC controller can be expressed as follows:(1)xsys=fxsys,ysys,usysxsys0=x00=gxsys,ysys,usysysys0=y0in which xsys, ysys, and usys are, respectively, the state, algebraic, and input vectors and are defined as(2)xsys=δi,ωi,Eqi′,Edi′,Efdi,VRi,RFiT,ysys=Vj,θj,Idi,IqiT,usys=TMi,VREFiT,i=1,…,m;j=1,…,n,whereT is the transpose operator,δ is the rotor angle of generator,ω is the speed of generator,V is the voltage magnitude of bus,θ is the power angle of bus,Id are the d-axis components of the current of generator,Iq are the q-axis components of the current of generator,VR is the input amplifier voltage of the excitation of generator,RF is the stabilizer feedback variable of the excitation of generator,TM is the electrical power of generator,Eq′ is the q-axis component of the internal voltage of generator,Ed′ is the d-axis component of the internal voltage of generator,Efd is the d-axis component of the field voltages of generator,VREF is the reference voltage of generator.Next, the linearized model is given as [40](3)ΔX˙sys00=Asys′Bsys1′Bsys2′Csys1′Dsys11′Dsys12′Csys2′Dsys21′Dsys22′ΔXsysΔysys1Δysys2+Esys00ΔUsys.It can be identified as(4)Asys=Asys′,Bsys=Bsys1′Bsys2′,Csys=Csys1′Csys2′,Dsys=Dsys11′Dsys12′Dsys21′Dsys22′,(5)Esys=Esys100.Equation (3) can be changed to another form as follows [40]:(6)Δx˙sys=AsysΔxsys+BsysΔysys+EsysΔusys,0=CsysΔxsys+DsysΔysys,where(7)Δxsys=Δδi,Δωi,ΔEqi′,ΔEdi′,ΔEfdi,ΔVRi,ΔRFiT,Δysys=Δysys1,Δysys1=θ1V1⋯Vm∣θ2⋯θnVm+1⋯VnT,Δusys=ΔTMi,ΔVREFiT.


### 2.3. DAE Model of Power System with TCSC Controller

#### 2.3.1. TCSC Controller

The main role of TCSC is to control fast the active power flow, increase the power transfer on transmission line, and enhance the stability of the power system. The basic structure TCSC consists of a fixed series capacitor bank C in parallel with a thyristor-controlled reactor (TCR), as shown in Figure3(a). It can control the continuous power flow on the alternating-current (AC) line with a variable series capacitive reactance. This series reactance is adjusted through variation of firing angle α; the effective reactance Xtcsc depends on three regions: (i) inductive region, which starts increasing from XL/XC value to infinity (XL(α)=XC); (ii) capacitive region, which starts increasing from infinity to capacitive reactance XC; and (iii) resonance region, which occurs between these two regions. In order to avoid the resonance region, the steady state limits of the firing angle are chosen to be 90° ≤ α ≤ 180° with the resonant point αr = 128.4°. In this paper, for investigating power system stability after a fault, TCSC should operate in capacitive region and the limit of the steady state to the firing angle is chosen to be αr≤α≤180°, as shown in Figure 3(c).Figure 3
The TCSC controller: (a) structure; (b) equivalent; (c) reactance versus firing angle characteristic cure.
(a) (b) (c)The variableXtcsc can be obtained by using some of the control strategies and the feedback signal of the TCSC controller. This feedback signal can be the signal of active power, reactive power, current of transmission line, or transmission angle [41]. In this paper, the main objective is to dampen the power oscillations. The active power signal in transmission line is chosen [42] and the control strategy for TCSC is shown in Figure 4 [34].Figure 4
The transfer function mode of TCSC controller.The new equivalent impedance of the line, when placed TCSC, as shown in Figure3(b), can be obtained as [43](8)Xnew_line=Xline+Xtcscin which the relationship between the firing angle α and the impedance of the TCSC at fundamental frequency can be derived as follows [44]:(9)Xtcsc=-XC+K12σ+sin⁡2σ-K2cos2⁡σϖtan⁡ϖσ-tan⁡σ,σ=π-α,ϖ=XCXL,XLC=XCXLXC-XL,K1=XC+XLCπ,K2=4XLC2πXL,where σ is the conduction angle, ϖ is the TCSC ratio. XL=ωL is the reactance of the inductor, and XL=1/ωC is the fixed capacitive impedance.


#### 2.3.2. The Unification of TCSC Controller in Power System

The TCSC has been instated on the transmission line between basesk and t of an n-bus power system. The injected power flow into buses k and t is given ib the following equations [45].At busk,(10)Pk=ΔGktVk2-VkVtΔGktcos⁡θk-θt+ΔBktsin⁡θk-θt,Qk=-ΔBktVk2-VkVtΔGktsin⁡θk-θt-ΔBktcos⁡θk-θt.Similarly, at bust,(11)Pt=Vt2ΔGkt-VkVtΔGtkcos⁡θk-θt+ΔBtksin⁡θk-θtQt=-Vt2ΔBkt+VkVtΔGktsin⁡θk-θt-ΔBtkcos⁡θk-θt,where Pk, Qk, Pt, and Qt are the active and reactive powers injected at buses k and t, respectively. Also, Vk, θk, Vt, and θt are voltage magnitudes and phase angles of buses k and t, respectively. ΔG and ΔB depend on TCSC reactance and are given as(12)ΔGtk=XtcscRlineXtcsc-2XlineRline2+Xline2Rline2+Xtcsc-Xline2,ΔBkt=XtcscRline2-Xline2+XtcscXlineRline2+Xline2Rline2+Xtcsc-Xline2,where Rline and Xline represent the resistance and reactance of the line, respectively, and Xtcsc is the optimal value of TCSC.The linearized model of TCSC is given as follows:(13)ΔX˙tcsc=AtcscΔXtcsc+BtcscΔθk,ΔVk,Δθt,ΔVtT,ΔPk,ΔQk,ΔPt,ΔQtT=CtcscΔXtcsc+DtcscΔθk,ΔVk,Δθt,ΔVtT.Incorporating (3) and (13), the DAE model of the multimachine power system with TCSC controller can be described as follows:(14)ΔX˙sysΔX˙tcsc00=Asys′A1tcscBsys1′Bsys2′A2tcscAtcscB1tcsc′B2tcsc′Csys1′C1tcscDsys11′Dsys12′Csys2′C2tcscD1tcsc′D2tcsc′ΔXsysΔXtcscΔysys1Δysys2+Esys-tcsc000ΔUsys-tcsc.It can be identified as(15)Anew=Asys′A1tcscA2tcscAtcsc,Bnew=Bsys1′Bsys2′Btcsc′Btcsc′,Cnew=Csys1′C1tcscCsys2′C2tcsc,Dnew=Dsys11′Dsys12′D1tcsc′D2tcsc′,Enew=Esys-tcsc,0,0,0T.Therefore, (14) can be changed to another form as follows:(16)Δx˙new=AnewΔx+BnewΔy+EnewΔu,0=CnewΔx+DnewΔy.


## 2.1. Characteristics of the Vietnamese Power System

The Vietnamese 500 kV power system operated in 1994. The length of transmission line of about 1,500 km was connected from the HoaBinh to PhuLam power stations with total rating of 2,700 MVA. It reached the length of about 5690 km and total capacity of all the 500 kV substations of 22,800 MVA in 2015. Several areas of 500 kV transmission lines are compensated by shunt reactors of about 70% and series capacitors of about 60%. The total generation capacity of the Vietnamese power system could reach 60,000 MW by the end of 2020 [39]. The result of load flow calculation on the 500 kV power system is shown in Figure 2.Figure 2
The result of load flow calculation on the Vietnamese 500/220 kV power system 2020.


## 2.2. DAE Model of Power System

The methodological method for dynamic modeling of generalm-machine and n-bus power system has been described in [40] and is applied for this work. In this model, each synchronous generator is represented by two-axis flux decay dynamic model along with IEEE type I. The differential algebraic equation (DAE) model of the power system without TCSC controller can be expressed as follows:(1)xsys=fxsys,ysys,usysxsys0=x00=gxsys,ysys,usysysys0=y0in which xsys, ysys, and usys are, respectively, the state, algebraic, and input vectors and are defined as(2)xsys=δi,ωi,Eqi′,Edi′,Efdi,VRi,RFiT,ysys=Vj,θj,Idi,IqiT,usys=TMi,VREFiT,i=1,…,m;j=1,…,n,whereT is the transpose operator,δ is the rotor angle of generator,ω is the speed of generator,V is the voltage magnitude of bus,θ is the power angle of bus,Id are the d-axis components of the current of generator,Iq are the q-axis components of the current of generator,VR is the input amplifier voltage of the excitation of generator,RF is the stabilizer feedback variable of the excitation of generator,TM is the electrical power of generator,Eq′ is the q-axis component of the internal voltage of generator,Ed′ is the d-axis component of the internal voltage of generator,Efd is the d-axis component of the field voltages of generator,VREF is the reference voltage of generator.Next, the linearized model is given as [40](3)ΔX˙sys00=Asys′Bsys1′Bsys2′Csys1′Dsys11′Dsys12′Csys2′Dsys21′Dsys22′ΔXsysΔysys1Δysys2+Esys00ΔUsys.It can be identified as(4)Asys=Asys′,Bsys=Bsys1′Bsys2′,Csys=Csys1′Csys2′,Dsys=Dsys11′Dsys12′Dsys21′Dsys22′,(5)Esys=Esys100.Equation (3) can be changed to another form as follows [40]:(6)Δx˙sys=AsysΔxsys+BsysΔysys+EsysΔusys,0=CsysΔxsys+DsysΔysys,where(7)Δxsys=Δδi,Δωi,ΔEqi′,ΔEdi′,ΔEfdi,ΔVRi,ΔRFiT,Δysys=Δysys1,Δysys1=θ1V1⋯Vm∣θ2⋯θnVm+1⋯VnT,Δusys=ΔTMi,ΔVREFiT.


## 2.3. DAE Model of Power System with TCSC Controller

### 2.3.1. TCSC Controller

The main role of TCSC is to control fast the active power flow, increase the power transfer on transmission line, and enhance the stability of the power system. The basic structure TCSC consists of a fixed series capacitor bank C in parallel with a thyristor-controlled reactor (TCR), as shown in Figure3(a). It can control the continuous power flow on the alternating-current (AC) line with a variable series capacitive reactance. This series reactance is adjusted through variation of firing angle α; the effective reactance Xtcsc depends on three regions: (i) inductive region, which starts increasing from XL/XC value to infinity (XL(α)=XC); (ii) capacitive region, which starts increasing from infinity to capacitive reactance XC; and (iii) resonance region, which occurs between these two regions. In order to avoid the resonance region, the steady state limits of the firing angle are chosen to be 90° ≤ α ≤ 180° with the resonant point αr = 128.4°. In this paper, for investigating power system stability after a fault, TCSC should operate in capacitive region and the limit of the steady state to the firing angle is chosen to be αr≤α≤180°, as shown in Figure 3(c).Figure 3
The TCSC controller: (a) structure; (b) equivalent; (c) reactance versus firing angle characteristic cure.
(a) (b) (c)The variableXtcsc can be obtained by using some of the control strategies and the feedback signal of the TCSC controller. This feedback signal can be the signal of active power, reactive power, current of transmission line, or transmission angle [41]. In this paper, the main objective is to dampen the power oscillations. The active power signal in transmission line is chosen [42] and the control strategy for TCSC is shown in Figure 4 [34].Figure 4
The transfer function mode of TCSC controller.The new equivalent impedance of the line, when placed TCSC, as shown in Figure3(b), can be obtained as [43](8)Xnew_line=Xline+Xtcscin which the relationship between the firing angle α and the impedance of the TCSC at fundamental frequency can be derived as follows [44]:(9)Xtcsc=-XC+K12σ+sin⁡2σ-K2cos2⁡σϖtan⁡ϖσ-tan⁡σ,σ=π-α,ϖ=XCXL,XLC=XCXLXC-XL,K1=XC+XLCπ,K2=4XLC2πXL,where σ is the conduction angle, ϖ is the TCSC ratio. XL=ωL is the reactance of the inductor, and XL=1/ωC is the fixed capacitive impedance.


### 2.3.2. The Unification of TCSC Controller in Power System

The TCSC has been instated on the transmission line between basesk and t of an n-bus power system. The injected power flow into buses k and t is given ib the following equations [45].At busk,(10)Pk=ΔGktVk2-VkVtΔGktcos⁡θk-θt+ΔBktsin⁡θk-θt,Qk=-ΔBktVk2-VkVtΔGktsin⁡θk-θt-ΔBktcos⁡θk-θt.Similarly, at bust,(11)Pt=Vt2ΔGkt-VkVtΔGtkcos⁡θk-θt+ΔBtksin⁡θk-θtQt=-Vt2ΔBkt+VkVtΔGktsin⁡θk-θt-ΔBtkcos⁡θk-θt,where Pk, Qk, Pt, and Qt are the active and reactive powers injected at buses k and t, respectively. Also, Vk, θk, Vt, and θt are voltage magnitudes and phase angles of buses k and t, respectively. ΔG and ΔB depend on TCSC reactance and are given as(12)ΔGtk=XtcscRlineXtcsc-2XlineRline2+Xline2Rline2+Xtcsc-Xline2,ΔBkt=XtcscRline2-Xline2+XtcscXlineRline2+Xline2Rline2+Xtcsc-Xline2,where Rline and Xline represent the resistance and reactance of the line, respectively, and Xtcsc is the optimal value of TCSC.The linearized model of TCSC is given as follows:(13)ΔX˙tcsc=AtcscΔXtcsc+BtcscΔθk,ΔVk,Δθt,ΔVtT,ΔPk,ΔQk,ΔPt,ΔQtT=CtcscΔXtcsc+DtcscΔθk,ΔVk,Δθt,ΔVtT.Incorporating (3) and (13), the DAE model of the multimachine power system with TCSC controller can be described as follows:(14)ΔX˙sysΔX˙tcsc00=Asys′A1tcscBsys1′Bsys2′A2tcscAtcscB1tcsc′B2tcsc′Csys1′C1tcscDsys11′Dsys12′Csys2′C2tcscD1tcsc′D2tcsc′ΔXsysΔXtcscΔysys1Δysys2+Esys-tcsc000ΔUsys-tcsc.It can be identified as(15)Anew=Asys′A1tcscA2tcscAtcsc,Bnew=Bsys1′Bsys2′Btcsc′Btcsc′,Cnew=Csys1′C1tcscCsys2′C2tcsc,Dnew=Dsys11′Dsys12′D1tcsc′D2tcsc′,Enew=Esys-tcsc,0,0,0T.Therefore, (14) can be changed to another form as follows:(16)Δx˙new=AnewΔx+BnewΔy+EnewΔu,0=CnewΔx+DnewΔy.


## 2.3.1. TCSC Controller

The main role of TCSC is to control fast the active power flow, increase the power transfer on transmission line, and enhance the stability of the power system. The basic structure TCSC consists of a fixed series capacitor bank C in parallel with a thyristor-controlled reactor (TCR), as shown in Figure3(a). It can control the continuous power flow on the alternating-current (AC) line with a variable series capacitive reactance. This series reactance is adjusted through variation of firing angle α; the effective reactance Xtcsc depends on three regions: (i) inductive region, which starts increasing from XL/XC value to infinity (XL(α)=XC); (ii) capacitive region, which starts increasing from infinity to capacitive reactance XC; and (iii) resonance region, which occurs between these two regions. In order to avoid the resonance region, the steady state limits of the firing angle are chosen to be 90° ≤ α ≤ 180° with the resonant point αr = 128.4°. In this paper, for investigating power system stability after a fault, TCSC should operate in capacitive region and the limit of the steady state to the firing angle is chosen to be αr≤α≤180°, as shown in Figure 3(c).Figure 3
The TCSC controller: (a) structure; (b) equivalent; (c) reactance versus firing angle characteristic cure.
(a) (b) (c)The variableXtcsc can be obtained by using some of the control strategies and the feedback signal of the TCSC controller. This feedback signal can be the signal of active power, reactive power, current of transmission line, or transmission angle [41]. In this paper, the main objective is to dampen the power oscillations. The active power signal in transmission line is chosen [42] and the control strategy for TCSC is shown in Figure 4 [34].Figure 4
The transfer function mode of TCSC controller.The new equivalent impedance of the line, when placed TCSC, as shown in Figure3(b), can be obtained as [43](8)Xnew_line=Xline+Xtcscin which the relationship between the firing angle α and the impedance of the TCSC at fundamental frequency can be derived as follows [44]:(9)Xtcsc=-XC+K12σ+sin⁡2σ-K2cos2⁡σϖtan⁡ϖσ-tan⁡σ,σ=π-α,ϖ=XCXL,XLC=XCXLXC-XL,K1=XC+XLCπ,K2=4XLC2πXL,where σ is the conduction angle, ϖ is the TCSC ratio. XL=ωL is the reactance of the inductor, and XL=1/ωC is the fixed capacitive impedance.


## 2.3.2. The Unification of TCSC Controller in Power System

The TCSC has been instated on the transmission line between basesk and t of an n-bus power system. The injected power flow into buses k and t is given ib the following equations [45].At busk,(10)Pk=ΔGktVk2-VkVtΔGktcos⁡θk-θt+ΔBktsin⁡θk-θt,Qk=-ΔBktVk2-VkVtΔGktsin⁡θk-θt-ΔBktcos⁡θk-θt.Similarly, at bust,(11)Pt=Vt2ΔGkt-VkVtΔGtkcos⁡θk-θt+ΔBtksin⁡θk-θtQt=-Vt2ΔBkt+VkVtΔGktsin⁡θk-θt-ΔBtkcos⁡θk-θt,where Pk, Qk, Pt, and Qt are the active and reactive powers injected at buses k and t, respectively. Also, Vk, θk, Vt, and θt are voltage magnitudes and phase angles of buses k and t, respectively. ΔG and ΔB depend on TCSC reactance and are given as(12)ΔGtk=XtcscRlineXtcsc-2XlineRline2+Xline2Rline2+Xtcsc-Xline2,ΔBkt=XtcscRline2-Xline2+XtcscXlineRline2+Xline2Rline2+Xtcsc-Xline2,where Rline and Xline represent the resistance and reactance of the line, respectively, and Xtcsc is the optimal value of TCSC.The linearized model of TCSC is given as follows:(13)ΔX˙tcsc=AtcscΔXtcsc+BtcscΔθk,ΔVk,Δθt,ΔVtT,ΔPk,ΔQk,ΔPt,ΔQtT=CtcscΔXtcsc+DtcscΔθk,ΔVk,Δθt,ΔVtT.Incorporating (3) and (13), the DAE model of the multimachine power system with TCSC controller can be described as follows:(14)ΔX˙sysΔX˙tcsc00=Asys′A1tcscBsys1′Bsys2′A2tcscAtcscB1tcsc′B2tcsc′Csys1′C1tcscDsys11′Dsys12′Csys2′C2tcscD1tcsc′D2tcsc′ΔXsysΔXtcscΔysys1Δysys2+Esys-tcsc000ΔUsys-tcsc.It can be identified as(15)Anew=Asys′A1tcscA2tcscAtcsc,Bnew=Bsys1′Bsys2′Btcsc′Btcsc′,Cnew=Csys1′C1tcscCsys2′C2tcsc,Dnew=Dsys11′Dsys12′D1tcsc′D2tcsc′,Enew=Esys-tcsc,0,0,0T.Therefore, (14) can be changed to another form as follows:(16)Δx˙new=AnewΔx+BnewΔy+EnewΔu,0=CnewΔx+DnewΔy.


## 3. The Proposed Approach

### 3.1. Gramian-Based Controllability and Observability

The system often has two properties, controllability and observability, which play an important role regarding the determination of the optimal location of TCSC in the power system. From that, the input and output variables need to be used in order to observe and control the system. Therefore, (16) can be redescribed in a state-space form as follows:(17)x˙t=Axt+But,yt=Cxt+Dut,wherex(t) is the state vector,y(t) is the output vector,u(t) is the input vector,A is the state matrix,B is the control matrix,C is the output matrix,D is the feed-forward matrix.


#### 3.1.1. Controllability

Definition 1.
System (17) is controllable over the interval [t0,t1], if any states (x0,x1∈R), there exists the input u(t): [t0,t1]→R, that drives the system from x0(t0) to x1(t1); this is true: the system is completely controllable despite the initial time and state.Every actuator in the power systems is the energy that can be limited, such that controllability matrix has been used for the purpose of dealing with the amount of input energy. This input energy is required to reach a given state from the origin. The property of controllability of the system can be described in a quantitative manner by the controllability function. Typically, it is defined for dynamic system as the minimum input energy that necessary drives the system from statex(-T)=0 to state x(0)=x0 and can be given as [46](18)Lcx0=minimizeu∈L2-T,0∫-T0ut2dt,subject  to:x˙t=Axt+But,x-T=0,x0=x0.It can be proven that the above system (17) has the transient controllability function, which is given as follows:(19)Lcx0=12x0TGc-1Tx0.Remark 2.
For the real dynamic system,Lc(x0) is large and the state x0 is difficult to reach. Thus, the system is uncontrollable.Remark 3.
MatricesA and B of the above system (17) are controllable in the time range (t0,t1), if and only if the controllability matrix Δ defined as (20) has full rank n, where n is the number of states. Notice that the controllability matrix has dimension n×nm, where m is the dimension of the input vector u.(20)Δ=BABA2B⋯An-1B.Remark 4.
MatricesA and B of the above system (17) are controllable if and only if the controllability Gramian matrix Gc(T) on horizon T defined as (21) has full rank n and is positive definite.(21)GcT=∫0Texp⁡AtBBTexp⁡ATtdt.


#### 3.1.2. Observability

The system states are the internal variables, which are hard to directly measure but the outputs can be measured easily at the same time. The observability property plays an important role in the analysis of optimal location.Definition 5.
System (17) is observable over the interval [t0,t1], if any states x0,x1∈R, there exists the output y(t): t0≤t≤t1<∞, assuming that the input u(t) is known; this is true: the system is observable despite the initial time and state.The observability property of the system can be characterized in a quantitative manner by the observability function. It is defined for dynamic system as output energy generated by the statex(0)=x0 (when u=0) and is given as [46](22)Lox0=12∫-T0yt2dt,subject  to:x0=x0,u≡0.It can be proven that system (17) has the transient observability function, which is given as follows:(23)Lox0=12x0TGo-1x0.Remark 6.
For the real dynamic system,Lo(x0) measures the effect of the initial condition on the output; if Lo(x0) is small, the effect of x0 in the output is confined, such that it is difficult to reach state x0. Thus, the system is unobservable.Remark 7.
MatricesA and C of the above system (17) are observable in the interval [t0,t1], if and only if the observability matrix ∇ defined as (24) has full rank n, where n is the number of states. Notice that the observability matrix has dimension np×n, where p is the dimension of the output vector y.(24)∇=CCACA2⋮CAn-1.Remark 8.
MatricesA and C of the above system (17) are observable if and only if the observability Gramian matrix Go(T) on horizon T defined as (25) has full rank n′ and is positive definite.(25)GoT=∫0Texp⁡ATtCTCexp⁡Atdt.If system (17) is asymptotically stability around the origin, the controllability and observability functions can be given by (26)Lc=minimizeu∈L2-∞,0x0=x0x-∞=012∫-∞0ut2dt,Lo=12∫0∞yt2dt,x0=x0,u≡0.Formal remarks for controllability and observability are white and black. In reality, some states are very calamitous to control or have little effect on the outputs. The degree of observability and controllability can be evaluated by the sizes ofGc and Go with infinite time horizon:(27)Gc0,∞=limT→∞∫0Texp⁡AtBBTexp⁡ATtdt,G00,∞=limT→∞∫0Texp⁡ATtCTCexp⁡Atdt.Notice that Gc and Go satisfy the following differential expression:(28)dGc0,tdt=AGc0,t+Gc0,tAT+BBT,dGo0,tdt=ATGo0,t+Go0,tA+CTC.It is hard to directly calculate the Gramian matrices from expression (27) because they consist of an exponential matrix and an integral. If all the eigenvalues of A have strictly negative real parts, as T→∞, there exists an easier way to calculate these matrices by solving the Lyapunov expression (29) below [47], which can be solved using MATLAB via gram. The obtained solutions are the unique positive definite.(29)AGc+GcAT+BBT=0,ATGo+GoA+CTC=0.Obviously, expression (29) shows that (i) the property of observability Gramian matrix depends on the output matrix A and state matrix C; accordingly, the output energy can be affected by properly choosing this matrix. (ii) The property of controllability Gramian matrix depends on the control matrix A and state matrix B; accordingly, the control energy can be affected by properly choosing this matrix. (iii) When the system is only detectable, the Gramian matrices will be only some nonnegative definite matrices.In particular, we can pose the problem seeking the minimum input energy (i.e., the control energy) that must derive system (17) from the initial state x0 to a final state xtf at time t=tf and this energy is defined by [34](30)J=exp⁡Atfx0-xtfTGc-exp⁡AtfGcexp⁡ATtfTexp⁡Atfx0-xtf.Clearly, expression (30) shows that minimizing J is the same as maximizing the Gramian matrix Gc. Furthermore, maximizing the transmitted total energy (potential and kinetic energy) from the actuators to the structure for a given input energy, it can be obtained as follows(31)E=∫0∞Ept+Ektdt=∫0∞xtxTtdt,where Ek is the kinetic energy and Ep is the potential energy. In this paper, we suppose that the input is a Dirac impulse that excites all frequencies and xt=exp⁡AtB is the impulse matrix. Hence, the obtained energy is the weightedtrace of controllability Gramian [34, 48]:(32)E=traceGc=∑i=1nσi,where σi is the ith Hankel singular value. The decomposition of this value denotes spatial energy decomposition contained in the impulse response. The sum of the singular values denotes the energy total in which each of the singular value represents the energy in a particular direction. The optimal location is determined based on the maximization of the total energy of the system at the fault clearing time. Therefore, (32) is a basis for selecting optimal control input.


### 3.2. Hankel-Norm Method

In order to compute easily, we focus on analyzing the infinite horizon Gramian. The Gramian matrices are computed by using expression (29) rather than solving expressions (21) and (25). However, it may be more difficult for analyzing the Vietnamese power system because that is a large-scale network. With such networks, the number of state variables could be big, leading to lost time for computing. In reality, determining the optimal location, we just need the number of major state variables that play an important role in order to analyze the linear system on the small-signal stability. Therefore, in order to solve this problem, the order reduction method is applied based on the computation of controllability and observability Gramian matrices introduced in [36].


#### 3.2.1. Balanced Realization

Considering system (17), in order to make effect for this study, the so-called similarity transformation M must need to modify the state variables and matrices of the system, but the output and input behaviors remain unchanged, such that the controllability and observability Gramian matrices satisfy the following:(33)Gc=Go=σ100⋯00σ20⋯0⋮⋮⋮⋱⋮000⋯σnwithσ1≥σ2≥⋯≥σn≥0.The so-called similarity transformation transforms the controllability and observability Gramian matrices in the following way:(34)G-c=MGcMT,G-o=M-1TGcM-1.However, the result turns out to be the invariant output and input behaviors, that is,MGcMT=M-1TGcM-1 (i.e., it relates to the Gramian of the original system by M).Expression (34) constructs a state-space transformation that is considered with the form f(t)=Mx(t); substituting this state-space transformation into (17), the new system can be obtained:(35)f˙t=A-ft+B-ut,yt=C-ft+D-ut,where A-=MAM-1; B-=MB; C-=CM-1; D-=D. Therefore, system (35) is a balanced form according to the given nonsingular transformation matrix M to change the system in such a way that the controllability and observability Gramian matrices satisfy the following condition:(36)G-c=G-o=diag⁡σ1,σ1,…,σn,σ1≥σ1≥⋯≥σn,in whichσ1,σ2,…,σn are real and positive numbers called the Hankel singular values and the transformation matrix is determined by using the algorithm [47] as shown in the Appendix.


#### 3.2.2. Order Reduction Realization

After obtaining the new system in a balanced form as shown in (35), the order reduction can be performed. The first step is to choose the order of k based on the purposes of reduction. The author in [36] has considered general situation that the Hankel singular values satisfy the condition, that is,(37)σ1≥σ2≥⋯≥σk≥σk+1=σk+2=⋯=σk+r>σk+r+1≥⋯≥σn>0.In this study, the system reduction is performed by eliminating all state variables corresponding to the Hankel singular values that are smaller than 10−3. Therefore, the system can be converted to the following order reduction form:(38)f˙1t=A-11f1t+A-12f2t+B-1utf˙2t=A-21f1t+A-22f2t+B-2ut;yt=C-1f1t+C-2f2t+D-ut,G-c=G-o=Σ-,where Σ-=[(diag⁡σ1,…,σk,σk+r+1,…,σn,)2-σk+12Ψr]=(Σ-12-σk+12Ψr) is the nonsingular and diagonal matrix. The vectors f1(t) and f2(t) are obtained by decomposing the state variable f(t) (i.e., f(t)=f1(t)f2(t)T) and have dimensions n-r and r, respectively. Ψr is the identity matrix and has dimension r. The matrices A-11,A-12;A-21,A-222×2, B-1;B-21×2, and C-1,C-22×1 are, respectively, partitioned from matrices A-, B-, and C- in balanced system (35) corresponding to the partitioned Gramian.If the number of inputs is more than that of outputs (the so-called all-pass dilation of the system), then,A-, B-, C-, and D- are replaced by A-T, B-T, C-T, and D-T, respectively [49]. Thus, the components for the so-called all-pass dilation of the system are defined as follows:(39)A~=Σ--1σk+12A-11T+Σ-1A-11Σ-1-σk+1C-1TUB-1T,B~=Σ--1Σ-1B-1+σk+1C-1TU,C~=C-1Σ-1+σk+1UB-1T,D~=D--σk+1U;Σ-=Σ-12-σk+12Ψr,where U is a unitary matrix and is solution of the following expression:(40)B-2+C-2TU=0.According to the conclusion in [49], the all-pass dilation of the system might be unstable. The final step is to compute the stable part of the all-pass dilation of the systems A~′, B~′, C~′, and D~′ based on the Hankel-norm method. As we consider a system that the number of inputs is more than the number of outputs, matrices A~′, B~′, C~′, and D~′ can be replaced by A~′T, B~′T, C~′T, and D~′T based on the following error of approximation:(41)Ts-Tks∞=σk+1,where T(s) and Tk(s) are the transfer functions of the original and the kth-order reduced systems.Clearly, expression (41) shows that the largest error for all frequencies is equal to the Hankel singular value k+1. It allows us to select a suitable value of k. Therefore, the matrices of the reduced system are obtained as follows:(42)Ak=A~′,Bk=B~′,Ck=C~′,Dk=D~′.


### 3.3. Algorithm

The main objective is to dampen the active power and generator angle oscillations in the power system when the active power perturbation occurs in the transmission line. The proposed algorithm is used to determine the optimal location of TCSC based on the Gramian critical energy, and the flowchart of the algorithm is shown in Figure5. The principles and details have been presented in Sections 3.1 and 3.2 and below is the interpretation.Figure 5
Flowchart of the proposed algorithm.First. Start.Step 1.
Place the TCSC controller in the lineith (i=1,2,3,…,n).Step 2.
Run steady state.Step 3.
Generate the active power perturbation signal in the transmission linejth (j=1,2,3,…,m) in which m is large enough. For each contingency of the active power perturbation j, the observed output υj and control input ξj were analyzed.Step 4.
Compute the matricesA, B, and C corresponding to the placement i and the active power perturbation signal j.Step 5.
Perform the order reduction for system based on the Hankel-norm method.Step 6.
Estimate the stable condition based on the state matrixAkij to exclude the unstable cases. If the condition is satisfied, Step 6 will be performed, and, vice versa, Step 1 will be iterated.Step 7.
Compute the controllability Gramian matrix of the new system (after performing the order reduction) corresponding to the active power perturbation signaljth.Step 8.
Iterate the steps from2 to 7 for computing a set of the active power perturbations (j=1,2,3,…,m).Step 9.
Compute the energy for each placement corresponding to a set of the active power perturbations(j=1,2,3,…,m).Step 10.
Iterate the steps from1 to 7 to calculate all the locations (i=1,2,3,…,n).Step 11.
Compare the maximum total energy to evaluate the optimal placement for TCSC controller.Finally. End.


## 3.1. Gramian-Based Controllability and Observability

The system often has two properties, controllability and observability, which play an important role regarding the determination of the optimal location of TCSC in the power system. From that, the input and output variables need to be used in order to observe and control the system. Therefore, (16) can be redescribed in a state-space form as follows:(17)x˙t=Axt+But,yt=Cxt+Dut,wherex(t) is the state vector,y(t) is the output vector,u(t) is the input vector,A is the state matrix,B is the control matrix,C is the output matrix,D is the feed-forward matrix.


### 3.1.1. Controllability

Definition 1.
System (17) is controllable over the interval [t0,t1], if any states (x0,x1∈R), there exists the input u(t): [t0,t1]→R, that drives the system from x0(t0) to x1(t1); this is true: the system is completely controllable despite the initial time and state.Every actuator in the power systems is the energy that can be limited, such that controllability matrix has been used for the purpose of dealing with the amount of input energy. This input energy is required to reach a given state from the origin. The property of controllability of the system can be described in a quantitative manner by the controllability function. Typically, it is defined for dynamic system as the minimum input energy that necessary drives the system from statex(-T)=0 to state x(0)=x0 and can be given as [46](18)Lcx0=minimizeu∈L2-T,0∫-T0ut2dt,subject  to:x˙t=Axt+But,x-T=0,x0=x0.It can be proven that the above system (17) has the transient controllability function, which is given as follows:(19)Lcx0=12x0TGc-1Tx0.Remark 2.
For the real dynamic system,Lc(x0) is large and the state x0 is difficult to reach. Thus, the system is uncontrollable.Remark 3.
MatricesA and B of the above system (17) are controllable in the time range (t0,t1), if and only if the controllability matrix Δ defined as (20) has full rank n, where n is the number of states. Notice that the controllability matrix has dimension n×nm, where m is the dimension of the input vector u.(20)Δ=BABA2B⋯An-1B.Remark 4.
MatricesA and B of the above system (17) are controllable if and only if the controllability Gramian matrix Gc(T) on horizon T defined as (21) has full rank n and is positive definite.(21)GcT=∫0Texp⁡AtBBTexp⁡ATtdt.


### 3.1.2. Observability

The system states are the internal variables, which are hard to directly measure but the outputs can be measured easily at the same time. The observability property plays an important role in the analysis of optimal location.Definition 5.
System (17) is observable over the interval [t0,t1], if any states x0,x1∈R, there exists the output y(t): t0≤t≤t1<∞, assuming that the input u(t) is known; this is true: the system is observable despite the initial time and state.The observability property of the system can be characterized in a quantitative manner by the observability function. It is defined for dynamic system as output energy generated by the statex(0)=x0 (when u=0) and is given as [46](22)Lox0=12∫-T0yt2dt,subject  to:x0=x0,u≡0.It can be proven that system (17) has the transient observability function, which is given as follows:(23)Lox0=12x0TGo-1x0.Remark 6.
For the real dynamic system,Lo(x0) measures the effect of the initial condition on the output; if Lo(x0) is small, the effect of x0 in the output is confined, such that it is difficult to reach state x0. Thus, the system is unobservable.Remark 7.
MatricesA and C of the above system (17) are observable in the interval [t0,t1], if and only if the observability matrix ∇ defined as (24) has full rank n, where n is the number of states. Notice that the observability matrix has dimension np×n, where p is the dimension of the output vector y.(24)∇=CCACA2⋮CAn-1.Remark 8.
MatricesA and C of the above system (17) are observable if and only if the observability Gramian matrix Go(T) on horizon T defined as (25) has full rank n′ and is positive definite.(25)GoT=∫0Texp⁡ATtCTCexp⁡Atdt.If system (17) is asymptotically stability around the origin, the controllability and observability functions can be given by (26)Lc=minimizeu∈L2-∞,0x0=x0x-∞=012∫-∞0ut2dt,Lo=12∫0∞yt2dt,x0=x0,u≡0.Formal remarks for controllability and observability are white and black. In reality, some states are very calamitous to control or have little effect on the outputs. The degree of observability and controllability can be evaluated by the sizes ofGc and Go with infinite time horizon:(27)Gc0,∞=limT→∞∫0Texp⁡AtBBTexp⁡ATtdt,G00,∞=limT→∞∫0Texp⁡ATtCTCexp⁡Atdt.Notice that Gc and Go satisfy the following differential expression:(28)dGc0,tdt=AGc0,t+Gc0,tAT+BBT,dGo0,tdt=ATGo0,t+Go0,tA+CTC.It is hard to directly calculate the Gramian matrices from expression (27) because they consist of an exponential matrix and an integral. If all the eigenvalues of A have strictly negative real parts, as T→∞, there exists an easier way to calculate these matrices by solving the Lyapunov expression (29) below [47], which can be solved using MATLAB via gram. The obtained solutions are the unique positive definite.(29)AGc+GcAT+BBT=0,ATGo+GoA+CTC=0.Obviously, expression (29) shows that (i) the property of observability Gramian matrix depends on the output matrix A and state matrix C; accordingly, the output energy can be affected by properly choosing this matrix. (ii) The property of controllability Gramian matrix depends on the control matrix A and state matrix B; accordingly, the control energy can be affected by properly choosing this matrix. (iii) When the system is only detectable, the Gramian matrices will be only some nonnegative definite matrices.In particular, we can pose the problem seeking the minimum input energy (i.e., the control energy) that must derive system (17) from the initial state x0 to a final state xtf at time t=tf and this energy is defined by [34](30)J=exp⁡Atfx0-xtfTGc-exp⁡AtfGcexp⁡ATtfTexp⁡Atfx0-xtf.Clearly, expression (30) shows that minimizing J is the same as maximizing the Gramian matrix Gc. Furthermore, maximizing the transmitted total energy (potential and kinetic energy) from the actuators to the structure for a given input energy, it can be obtained as follows(31)E=∫0∞Ept+Ektdt=∫0∞xtxTtdt,where Ek is the kinetic energy and Ep is the potential energy. In this paper, we suppose that the input is a Dirac impulse that excites all frequencies and xt=exp⁡AtB is the impulse matrix. Hence, the obtained energy is the weightedtrace of controllability Gramian [34, 48]:(32)E=traceGc=∑i=1nσi,where σi is the ith Hankel singular value. The decomposition of this value denotes spatial energy decomposition contained in the impulse response. The sum of the singular values denotes the energy total in which each of the singular value represents the energy in a particular direction. The optimal location is determined based on the maximization of the total energy of the system at the fault clearing time. Therefore, (32) is a basis for selecting optimal control input.


## 3.1.1. Controllability

Definition 1.
System (17) is controllable over the interval [t0,t1], if any states (x0,x1∈R), there exists the input u(t): [t0,t1]→R, that drives the system from x0(t0) to x1(t1); this is true: the system is completely controllable despite the initial time and state.Every actuator in the power systems is the energy that can be limited, such that controllability matrix has been used for the purpose of dealing with the amount of input energy. This input energy is required to reach a given state from the origin. The property of controllability of the system can be described in a quantitative manner by the controllability function. Typically, it is defined for dynamic system as the minimum input energy that necessary drives the system from statex(-T)=0 to state x(0)=x0 and can be given as [46](18)Lcx0=minimizeu∈L2-T,0∫-T0ut2dt,subject  to:x˙t=Axt+But,x-T=0,x0=x0.It can be proven that the above system (17) has the transient controllability function, which is given as follows:(19)Lcx0=12x0TGc-1Tx0.Remark 2.
For the real dynamic system,Lc(x0) is large and the state x0 is difficult to reach. Thus, the system is uncontrollable.Remark 3.
MatricesA and B of the above system (17) are controllable in the time range (t0,t1), if and only if the controllability matrix Δ defined as (20) has full rank n, where n is the number of states. Notice that the controllability matrix has dimension n×nm, where m is the dimension of the input vector u.(20)Δ=BABA2B⋯An-1B.Remark 4.
MatricesA and B of the above system (17) are controllable if and only if the controllability Gramian matrix Gc(T) on horizon T defined as (21) has full rank n and is positive definite.(21)GcT=∫0Texp⁡AtBBTexp⁡ATtdt.


## 3.1.2. Observability

The system states are the internal variables, which are hard to directly measure but the outputs can be measured easily at the same time. The observability property plays an important role in the analysis of optimal location.Definition 5.
System (17) is observable over the interval [t0,t1], if any states x0,x1∈R, there exists the output y(t): t0≤t≤t1<∞, assuming that the input u(t) is known; this is true: the system is observable despite the initial time and state.The observability property of the system can be characterized in a quantitative manner by the observability function. It is defined for dynamic system as output energy generated by the statex(0)=x0 (when u=0) and is given as [46](22)Lox0=12∫-T0yt2dt,subject  to:x0=x0,u≡0.It can be proven that system (17) has the transient observability function, which is given as follows:(23)Lox0=12x0TGo-1x0.Remark 6.
For the real dynamic system,Lo(x0) measures the effect of the initial condition on the output; if Lo(x0) is small, the effect of x0 in the output is confined, such that it is difficult to reach state x0. Thus, the system is unobservable.Remark 7.
MatricesA and C of the above system (17) are observable in the interval [t0,t1], if and only if the observability matrix ∇ defined as (24) has full rank n, where n is the number of states. Notice that the observability matrix has dimension np×n, where p is the dimension of the output vector y.(24)∇=CCACA2⋮CAn-1.Remark 8.
MatricesA and C of the above system (17) are observable if and only if the observability Gramian matrix Go(T) on horizon T defined as (25) has full rank n′ and is positive definite.(25)GoT=∫0Texp⁡ATtCTCexp⁡Atdt.If system (17) is asymptotically stability around the origin, the controllability and observability functions can be given by (26)Lc=minimizeu∈L2-∞,0x0=x0x-∞=012∫-∞0ut2dt,Lo=12∫0∞yt2dt,x0=x0,u≡0.Formal remarks for controllability and observability are white and black. In reality, some states are very calamitous to control or have little effect on the outputs. The degree of observability and controllability can be evaluated by the sizes ofGc and Go with infinite time horizon:(27)Gc0,∞=limT→∞∫0Texp⁡AtBBTexp⁡ATtdt,G00,∞=limT→∞∫0Texp⁡ATtCTCexp⁡Atdt.Notice that Gc and Go satisfy the following differential expression:(28)dGc0,tdt=AGc0,t+Gc0,tAT+BBT,dGo0,tdt=ATGo0,t+Go0,tA+CTC.It is hard to directly calculate the Gramian matrices from expression (27) because they consist of an exponential matrix and an integral. If all the eigenvalues of A have strictly negative real parts, as T→∞, there exists an easier way to calculate these matrices by solving the Lyapunov expression (29) below [47], which can be solved using MATLAB via gram. The obtained solutions are the unique positive definite.(29)AGc+GcAT+BBT=0,ATGo+GoA+CTC=0.Obviously, expression (29) shows that (i) the property of observability Gramian matrix depends on the output matrix A and state matrix C; accordingly, the output energy can be affected by properly choosing this matrix. (ii) The property of controllability Gramian matrix depends on the control matrix A and state matrix B; accordingly, the control energy can be affected by properly choosing this matrix. (iii) When the system is only detectable, the Gramian matrices will be only some nonnegative definite matrices.In particular, we can pose the problem seeking the minimum input energy (i.e., the control energy) that must derive system (17) from the initial state x0 to a final state xtf at time t=tf and this energy is defined by [34](30)J=exp⁡Atfx0-xtfTGc-exp⁡AtfGcexp⁡ATtfTexp⁡Atfx0-xtf.Clearly, expression (30) shows that minimizing J is the same as maximizing the Gramian matrix Gc. Furthermore, maximizing the transmitted total energy (potential and kinetic energy) from the actuators to the structure for a given input energy, it can be obtained as follows(31)E=∫0∞Ept+Ektdt=∫0∞xtxTtdt,where Ek is the kinetic energy and Ep is the potential energy. In this paper, we suppose that the input is a Dirac impulse that excites all frequencies and xt=exp⁡AtB is the impulse matrix. Hence, the obtained energy is the weightedtrace of controllability Gramian [34, 48]:(32)E=traceGc=∑i=1nσi,where σi is the ith Hankel singular value. The decomposition of this value denotes spatial energy decomposition contained in the impulse response. The sum of the singular values denotes the energy total in which each of the singular value represents the energy in a particular direction. The optimal location is determined based on the maximization of the total energy of the system at the fault clearing time. Therefore, (32) is a basis for selecting optimal control input.


## 3.2. Hankel-Norm Method

In order to compute easily, we focus on analyzing the infinite horizon Gramian. The Gramian matrices are computed by using expression (29) rather than solving expressions (21) and (25). However, it may be more difficult for analyzing the Vietnamese power system because that is a large-scale network. With such networks, the number of state variables could be big, leading to lost time for computing. In reality, determining the optimal location, we just need the number of major state variables that play an important role in order to analyze the linear system on the small-signal stability. Therefore, in order to solve this problem, the order reduction method is applied based on the computation of controllability and observability Gramian matrices introduced in [36].


### 3.2.1. Balanced Realization

Considering system (17), in order to make effect for this study, the so-called similarity transformation M must need to modify the state variables and matrices of the system, but the output and input behaviors remain unchanged, such that the controllability and observability Gramian matrices satisfy the following:(33)Gc=Go=σ100⋯00σ20⋯0⋮⋮⋮⋱⋮000⋯σnwithσ1≥σ2≥⋯≥σn≥0.The so-called similarity transformation transforms the controllability and observability Gramian matrices in the following way:(34)G-c=MGcMT,G-o=M-1TGcM-1.However, the result turns out to be the invariant output and input behaviors, that is,MGcMT=M-1TGcM-1 (i.e., it relates to the Gramian of the original system by M).Expression (34) constructs a state-space transformation that is considered with the form f(t)=Mx(t); substituting this state-space transformation into (17), the new system can be obtained:(35)f˙t=A-ft+B-ut,yt=C-ft+D-ut,where A-=MAM-1; B-=MB; C-=CM-1; D-=D. Therefore, system (35) is a balanced form according to the given nonsingular transformation matrix M to change the system in such a way that the controllability and observability Gramian matrices satisfy the following condition:(36)G-c=G-o=diag⁡σ1,σ1,…,σn,σ1≥σ1≥⋯≥σn,in whichσ1,σ2,…,σn are real and positive numbers called the Hankel singular values and the transformation matrix is determined by using the algorithm [47] as shown in the Appendix.


### 3.2.2. Order Reduction Realization

After obtaining the new system in a balanced form as shown in (35), the order reduction can be performed. The first step is to choose the order of k based on the purposes of reduction. The author in [36] has considered general situation that the Hankel singular values satisfy the condition, that is,(37)σ1≥σ2≥⋯≥σk≥σk+1=σk+2=⋯=σk+r>σk+r+1≥⋯≥σn>0.In this study, the system reduction is performed by eliminating all state variables corresponding to the Hankel singular values that are smaller than 10−3. Therefore, the system can be converted to the following order reduction form:(38)f˙1t=A-11f1t+A-12f2t+B-1utf˙2t=A-21f1t+A-22f2t+B-2ut;yt=C-1f1t+C-2f2t+D-ut,G-c=G-o=Σ-,where Σ-=[(diag⁡σ1,…,σk,σk+r+1,…,σn,)2-σk+12Ψr]=(Σ-12-σk+12Ψr) is the nonsingular and diagonal matrix. The vectors f1(t) and f2(t) are obtained by decomposing the state variable f(t) (i.e., f(t)=f1(t)f2(t)T) and have dimensions n-r and r, respectively. Ψr is the identity matrix and has dimension r. The matrices A-11,A-12;A-21,A-222×2, B-1;B-21×2, and C-1,C-22×1 are, respectively, partitioned from matrices A-, B-, and C- in balanced system (35) corresponding to the partitioned Gramian.If the number of inputs is more than that of outputs (the so-called all-pass dilation of the system), then,A-, B-, C-, and D- are replaced by A-T, B-T, C-T, and D-T, respectively [49]. Thus, the components for the so-called all-pass dilation of the system are defined as follows:(39)A~=Σ--1σk+12A-11T+Σ-1A-11Σ-1-σk+1C-1TUB-1T,B~=Σ--1Σ-1B-1+σk+1C-1TU,C~=C-1Σ-1+σk+1UB-1T,D~=D--σk+1U;Σ-=Σ-12-σk+12Ψr,where U is a unitary matrix and is solution of the following expression:(40)B-2+C-2TU=0.According to the conclusion in [49], the all-pass dilation of the system might be unstable. The final step is to compute the stable part of the all-pass dilation of the systems A~′, B~′, C~′, and D~′ based on the Hankel-norm method. As we consider a system that the number of inputs is more than the number of outputs, matrices A~′, B~′, C~′, and D~′ can be replaced by A~′T, B~′T, C~′T, and D~′T based on the following error of approximation:(41)Ts-Tks∞=σk+1,where T(s) and Tk(s) are the transfer functions of the original and the kth-order reduced systems.Clearly, expression (41) shows that the largest error for all frequencies is equal to the Hankel singular value k+1. It allows us to select a suitable value of k. Therefore, the matrices of the reduced system are obtained as follows:(42)Ak=A~′,Bk=B~′,Ck=C~′,Dk=D~′.


## 3.2.1. Balanced Realization

Considering system (17), in order to make effect for this study, the so-called similarity transformation M must need to modify the state variables and matrices of the system, but the output and input behaviors remain unchanged, such that the controllability and observability Gramian matrices satisfy the following:(33)Gc=Go=σ100⋯00σ20⋯0⋮⋮⋮⋱⋮000⋯σnwithσ1≥σ2≥⋯≥σn≥0.The so-called similarity transformation transforms the controllability and observability Gramian matrices in the following way:(34)G-c=MGcMT,G-o=M-1TGcM-1.However, the result turns out to be the invariant output and input behaviors, that is,MGcMT=M-1TGcM-1 (i.e., it relates to the Gramian of the original system by M).Expression (34) constructs a state-space transformation that is considered with the form f(t)=Mx(t); substituting this state-space transformation into (17), the new system can be obtained:(35)f˙t=A-ft+B-ut,yt=C-ft+D-ut,where A-=MAM-1; B-=MB; C-=CM-1; D-=D. Therefore, system (35) is a balanced form according to the given nonsingular transformation matrix M to change the system in such a way that the controllability and observability Gramian matrices satisfy the following condition:(36)G-c=G-o=diag⁡σ1,σ1,…,σn,σ1≥σ1≥⋯≥σn,in whichσ1,σ2,…,σn are real and positive numbers called the Hankel singular values and the transformation matrix is determined by using the algorithm [47] as shown in the Appendix.


## 3.2.2. Order Reduction Realization

After obtaining the new system in a balanced form as shown in (35), the order reduction can be performed. The first step is to choose the order of k based on the purposes of reduction. The author in [36] has considered general situation that the Hankel singular values satisfy the condition, that is,(37)σ1≥σ2≥⋯≥σk≥σk+1=σk+2=⋯=σk+r>σk+r+1≥⋯≥σn>0.In this study, the system reduction is performed by eliminating all state variables corresponding to the Hankel singular values that are smaller than 10−3. Therefore, the system can be converted to the following order reduction form:(38)f˙1t=A-11f1t+A-12f2t+B-1utf˙2t=A-21f1t+A-22f2t+B-2ut;yt=C-1f1t+C-2f2t+D-ut,G-c=G-o=Σ-,where Σ-=[(diag⁡σ1,…,σk,σk+r+1,…,σn,)2-σk+12Ψr]=(Σ-12-σk+12Ψr) is the nonsingular and diagonal matrix. The vectors f1(t) and f2(t) are obtained by decomposing the state variable f(t) (i.e., f(t)=f1(t)f2(t)T) and have dimensions n-r and r, respectively. Ψr is the identity matrix and has dimension r. The matrices A-11,A-12;A-21,A-222×2, B-1;B-21×2, and C-1,C-22×1 are, respectively, partitioned from matrices A-, B-, and C- in balanced system (35) corresponding to the partitioned Gramian.If the number of inputs is more than that of outputs (the so-called all-pass dilation of the system), then,A-, B-, C-, and D- are replaced by A-T, B-T, C-T, and D-T, respectively [49]. Thus, the components for the so-called all-pass dilation of the system are defined as follows:(39)A~=Σ--1σk+12A-11T+Σ-1A-11Σ-1-σk+1C-1TUB-1T,B~=Σ--1Σ-1B-1+σk+1C-1TU,C~=C-1Σ-1+σk+1UB-1T,D~=D--σk+1U;Σ-=Σ-12-σk+12Ψr,where U is a unitary matrix and is solution of the following expression:(40)B-2+C-2TU=0.According to the conclusion in [49], the all-pass dilation of the system might be unstable. The final step is to compute the stable part of the all-pass dilation of the systems A~′, B~′, C~′, and D~′ based on the Hankel-norm method. As we consider a system that the number of inputs is more than the number of outputs, matrices A~′, B~′, C~′, and D~′ can be replaced by A~′T, B~′T, C~′T, and D~′T based on the following error of approximation:(41)Ts-Tks∞=σk+1,where T(s) and Tk(s) are the transfer functions of the original and the kth-order reduced systems.Clearly, expression (41) shows that the largest error for all frequencies is equal to the Hankel singular value k+1. It allows us to select a suitable value of k. Therefore, the matrices of the reduced system are obtained as follows:(42)Ak=A~′,Bk=B~′,Ck=C~′,Dk=D~′.


## 3.3. Algorithm

The main objective is to dampen the active power and generator angle oscillations in the power system when the active power perturbation occurs in the transmission line. The proposed algorithm is used to determine the optimal location of TCSC based on the Gramian critical energy, and the flowchart of the algorithm is shown in Figure5. The principles and details have been presented in Sections 3.1 and 3.2 and below is the interpretation.Figure 5
Flowchart of the proposed algorithm.First. Start.Step 1.
Place the TCSC controller in the lineith (i=1,2,3,…,n).Step 2.
Run steady state.Step 3.
Generate the active power perturbation signal in the transmission linejth (j=1,2,3,…,m) in which m is large enough. For each contingency of the active power perturbation j, the observed output υj and control input ξj were analyzed.Step 4.
Compute the matricesA, B, and C corresponding to the placement i and the active power perturbation signal j.Step 5.
Perform the order reduction for system based on the Hankel-norm method.Step 6.
Estimate the stable condition based on the state matrixAkij to exclude the unstable cases. If the condition is satisfied, Step 6 will be performed, and, vice versa, Step 1 will be iterated.Step 7.
Compute the controllability Gramian matrix of the new system (after performing the order reduction) corresponding to the active power perturbation signaljth.Step 8.
Iterate the steps from2 to 7 for computing a set of the active power perturbations (j=1,2,3,…,m).Step 9.
Compute the energy for each placement corresponding to a set of the active power perturbations(j=1,2,3,…,m).Step 10.
Iterate the steps from1 to 7 to calculate all the locations (i=1,2,3,…,n).Step 11.
Compare the maximum total energy to evaluate the optimal placement for TCSC controller.Finally. End.


## 4. Case Study

The Vietnamese 500/220 kV transmission system is used in the study to illustrate the effectiveness of the proposed method. This system consists of 29 substations of 500 kV, 162 substations of 220 kV, 16 double-circuit lines and 20 single-circuit lines of 500 kV, 205 double-circuit lines and 67 single-circuit lines of 220 kV, and 179 generator units. The generated total power is about 42,179 MW; the peak load demand is about 40,703 MW. The power system simulation engineering (PSS/E) program is used to analyze the transient and small-signal stability and MATLAB is used to calculate Gramian matrices and perform order reduction. All dynamic models such as generators, excitation systems, transmission lines, and loads are modeled by using PSS/E (from dynamic model library) [50], and all dynamic parameters are taken from [2] and the Vietnamese national load dispatch center (NLDC). The single-line diagram and the result of load flow calculation of the 500 kV voltage level are given in Figure 2.The optimal placement of TCSC in the Vietnamese network is determined based on (i) a combination of the controllability Gramian critical energy with the order reduction method on the small-signal stability analysis applied to search for several feasible locations; (ii) the transient stability analysis performed to compare and determine optimal locations through various simulation cases.In addition, the existing series compensation ratio on the Vietnamese 500 kV power system is considered to calculate and listed out in Table1.Table 1
The existing compensation ratio on the 500 kV system.
Case number The line between buses Compensation ratio (%) On circuit (1) NhoQuan-Hatinh (NQ-HT) 57.7 1 (2) NQ-HT 54.4 2 (3) VungAng-Danang (VA-DN) 64 1 (4) VA-DN 64 2 (5) Danang-Docsoi (DN-DS) 63 — (6) Docsoi-Pleiku (DS-Plei) 58 — (7) Thanhmy-Pleiku (TM-Plei) 52 — (8) Pleiku-Daknong (Plei-DakN) 76 — (9) Daknong-Caubong (DakN-CB) 61 — (10) Pleiku-Dilinh (Plei-DL) 70 — (11) Dilinh-Tandinh (DL-TD) 59 — (12) Pleiku-Caubong (Plei-CB) 70 —


### 4.1. Determine Feasible Locations

The correctness of the proposed method for the optimal location of TCSC is verified on the small-signal analysis; the dynamic model of TCSC is chosen as shown in Figure4 [34]. We suppose the compensation capacity of TCSC is set to the range [0.75, 1.2] of the uncompensated line. The active power perturbation signal in the line is considered as the input signal for the TCSC controller [51]. The optimal method in [38] was used to determine the parameters for the TCSC controller, and these parameters are listed in Table 2 [34].Table 2
Parameters of the TCSC controller.
Parameter Value Parameter Value T 1 0.1 sec K 0.75 T 2 0.1 sec X m a x 1.2Xline p.u. T 3 0.4 sec X m i n 0.25Xline p.u. T w 10 sec T t c s c 0.015 secThe studied contingency cases of active power perturbation were selected from single-line outage cases on the basis of the real power flow performance index (PI) introduced in [36], as follows:(43)PI=∑i=1nlυi2ϑPliPlimax2ϑ,where Pli and Plimax are the real power and rated power flow of the line i, respectively, ϑ is the exponent, and υi is a the real nonnegative weighting factor for line i. In this study, the value of the exponent was selected as 2.0 and 1.0 is the value of υi.Observing from (43) that all the power flows in the lines are within their limits, PI is small and, vice versa, it reaches out a high value when there are overloads. The studied contingency cases are listed in Table 3 by using (43); this corresponds to the active power perturbation signal in the line of the Vietnamese network.Table 3
Studied contingency cases.
Case number Active power perturbation signal in the line between buses (1) NhoQuan-HoaBinh (NQ-HB) (2) VietTri-Pitoong (VT-PT) (3) ThangLong-PhoNoi (TL-PN) (4) NhoQuan-SonLa (NQ-SL) (5) Pitoong-HoaBinh (PT-HB) (6) QuangNinh-HiepHoa (QN-HH) (7) HiepHoa-VietTri (HH-VT) (8) HiepHoa-Pitoong (HH-PT) (9) Pitoong-SonLa (PT-SL) (10) QuangNinh-ThangLong (QN-TL) (11) QuangNinh-PhoNoi (QN-PN) (12) QuangNinh-MongDuong (QN-MD) (13) ThuongTin-Nho Quan (TT-NQ) (14) ThuongTin-PhoNoi (TT-PN) (15) VungAng-HaTinh (VA-HT) (16) DaNang-ThanhMy (DN-TM) (17) SongMay-TanUyen (SM-TY) (18) SongMay-TanDinh (SM-TD) (19) NhaBe-PhuLam (NB-PL) (20) TanDinh-CauBong (TD-CB) (21) PhuLam-MyTho (PL-MT) (22) PhuLam-DucHoa (PL-DH) (23) DucHoa-MyTho (DH-MT) (24) MyTho-DuyenHai (MT-DH) (25) OMon-MyTho (OM-MT) (26) VinhTan-SongMay (VT-SM) (27) PhuMy-NhaBe (PM-NB) (28) CauBong-PhuLam (CB-PL) (29) NhaBe-MyTho (NB-MT) (30) CauBong-DucHoa (CB-DH) (31) SongMay-PhuMy (SM-PM)The Vietnamese 500/220 kV power system created the state matrix (A) having dimension [2371 × 2371] that corresponds with 179 generator units. By performing the equivalence of 179 generator units into 45 equivalent generator units, the state matrix now has dimension [515 × 515]. A TCSC controller has three states as the transducer, washout, and lead/lag added in the system. Therefore, the state matrix dimension of the equivalent power system model is [518 × 518].Once again, the state matrix of the system is reduced the dimension by applying the balanced order reduction method to eliminate the singular values that are smaller than 10−3. Figure 6 shows the eigenvalue of the system state matrix before and after using balanced reduction. Therefore, the state matrix (A11) of the new system has the dimension [45 × 45].Figure 6
Eigenvalue distribution before and after using balanced reduction.Figure7 shows the frequency response of the system before and after using balanced reduction. Obviously, the input-output behaviors of the system in both cases at the bandwidth range from 10−∞ to 105 are the same with the same control signal. The bandwidth has a range from 105 to 10∞; the frequency response of the reduced system is flat (a straight line) since its order is smaller than the original system. Therefore, we can conclude that two systems are equivalent in terms of the input-output behaviors.Figure 7
Frequency response of the original and reduced system.Figure8 shows the Hankel singular values after performing the balanced realization reduction technique. It can be observed from this figure that the 45 Hankel singular values have been finally retained in comparison with the 2374 state variables of the original system (including state variables of TCSC controller). This acquired singular values are smaller than the chosen threshold 10−3 (choice of this study). The eigenvalues distribution of the system before and after using balanced reduction is plotted in Figure 8. Therefore, the reduced system has 45 state variables, which can fully meet effects of the network to calculate the Gramian matrices.Figure 8
Distribution of Hankel singular values.The controllability Gramian energy indices were determined based on the proposed algorithm as shown in Section3.3 and the active power perturbation signals in the line as listed in Table 3 when the transient horizon is equal to infinity (T=∞) are given in Table 4, in which we only give several feasible locations having the high total energy value. From this table, it can be seen that the line PLei-CB is the suitable location for TCSC controller because it has the maximum total energy value; it means that the energy needed to drive the controllable state variables is smaller than other cases. Therefore, the line between buses PLei and CB is considered as the best location to install the TCSC controller.Table 4
Energy values according to several feasible locations of TCSC.
Case number Contingency cases (from Table3) (Trace ofGC) × 105 TCSC is placed on the line between buses None NQ-HT DN-VA PLei-TM PLei-DS DN-DS PLei-DL PLei-CB PLei-DakN DakN-CB DL-TD The Northern region (  1) (1) NQ-HB 0.0530 0.1354 0.1189 0.0780 0.0516 0.0614 0.0693 0.0780 0.0626 0.0290 0.0381 (2) VT-PT 0.0020 0.0142 0.0132 0.0020 0.0067 0.0080 0.0103 0.0118 0.0089 0.0055 0.0064 (3) TL-PN 0.0012 0.0947 0.0873 0.0374 0.0368 0.0483 0.0506 0.0526 0.0450 0.0215 0.0254 (4) NQ-SL 0.0242 0.0923 0.0840 0.0520 0.0354 0.0429 0.0458 0.0502 0.0418 0.0185 0.0242 (5) PT-HB 0.0034 0.0606 0.0575 0.0120 0.0240 0.0294 0.0307 0.0326 0.0282 0.0122 0.0156 (6) QN-HH 0.0010 0.0108 0.0099 0.0011 0.0044 0.0043 0.0048 0.0059 0.0047 0.0029 0.0032 (7) HH-VT 0.0012 0.0095 0.0097 0.0025 0.0044 0.0056 0.0061 0.0059 0.0054 0.0030 0.0031 (8) HH-PT 0.0032 0.0129 0.0122 0.0040 0.0061 0.0074 0.0092 0.0101 0.0079 0.0060 0.0055 (9) PT-SL 0.0010 0.0137 0.0117 0.0064 0.0059 0.0031 0.0074 0.0109 0.0071 0.0043 0.0057 (10) QN-TL 0.0029 0.0613 0.0564 0.0135 0.0242 0.0307 0.0337 0.0364 0.0301 0.0156 0.0176 (11) QN-PN 0.0032 0.0957 0.0887 0.0129 0.0372 0.0489 0.0514 0.0530 0.0458 0.0223 0.0253 (12) QN-MD 0.0216 0.0361 0.0343 0.0160 0.0140 0.0197 0.0194 0.0188 0.0171 0.0082 0.0090 (13) TT-NQ 0.0054 0.1523 0.1368 0.0850 0.0600 0.0766 0.0864 0.0941 0.0763 0.0382 0.0450 (14) TT-PN 0.0152 0.1173 0.1072 0.0630 0.0460 0.0598 0.0646 0.0685 0.0574 0.0279 0.0330 Total energy of (  1) 0.1385 0.9067 0.8278 0.3858 0.3569 0.4461 0.4897 0.5288 0.4383 0.2151 0.2571  The Southern region (  2) (15) VA-HT 0.0182 0.2777 0.2887 0.1660 0.1227 0.1616 0.1759 0.1900 0.1546 0.0780 0.0930 (16) DN-TM 0.0287 0.2015 0.4708 0.4982 0.5948 0.1894 0.3267 0.3003 0.2733 0.1642 0.1462 (17) SM-TY 0.0019 0.0037 0.0035 0.0002 0.0015 0.0022 0.0656 0.0125 0.0132 0.0296 0.0603 (18) SM-TD 0.0085 0.0453 0.0646 0.0399 0.0412 0.0416 0.3961 0.1184 0.0281 0.1382 0.3472 (19) NP-PL 0.0043 0.0683 0.0829 0.0492 0.0502 0.0460 0.1330 0.3648 0.1638 0.3823 0.0413 (20) TD-CB 0.0039 0.0268 0.0380 0.0359 0.0248 0.0232 0.1742 0.1055 0.5285 0.1123 0.1103 (21) PL-MT 0.0029 0.0158 0.0253 0.0229 0.0167 0.0188 0.0607 0.1077 0.0451 0.1217 0.0498 (22) PL-DH 0.0013 0.0079 0.0062 0.0013 0.0030 0.0021 0.0381 0.0175 0.0092 0.0255 0.0262 (23) DH-MT 0.0048 0.0180 0.0273 0.0095 0.0177 0.0189 0.0477 0.1083 0.0492 0.1166 0.0413 (24) MT-DHa 0.0018 0.0077 0.0174 0.0136 0.0117 0.0133 0.0425 0.0612 0.0269 0.0659 0.0312 (25) OM-MT 0.0023 0.0176 0.0256 0.0271 0.0166 0.0177 0.0672 0.0930 0.0419 0.1001 0.0498 (26) VT-SM 0.0082 0.0120 0.0289 0.0124 0.0195 0.0198 0.0969 0.0826 0.0343 0.0803 0.0911 (27) PM-NB 0.0011 0.0370 0.0341 0.0340 0.0184 0.0111 0.1370 0.1166 0.0620 0.0931 0.0798 (28) CB-PL 0.0043 0.0815 0.1142 0.0150 0.0734 0.0774 0.0584 0.5821 0.2359 0.6319 0.1125 (29) NB-MT 0.0034 0.0078 0.0114 0.0100 0.0077 0.0089 0.0774 0.0357 0.0123 0.0391 0.0502 (30) CB-DH 0.0150 0.0346 0.0499 0.0380 0.0330 0.0357 0.0347 0.2696 0.1048 0.3000 0.0600 (31) SM-PM 0.0112 0.0413 0.0424 0.0121 0.0245 0.0202 0.2338 0.0395 0.0303 0.0146 0.1432 Total energy of (  2) 0.1218 0.9043 1.3312 0.9853 1.0773 0.7079 2.1659 2.6053 1.8134 2.4934 1.5332  Total energy of (  1) + (  2) 0.2604 1.8111 2.1590 1.3711 1.4342 1.1540 2.6556 3.1341 2.2517 2.7085 1.7903


### 4.2. Retest Transient Stability

To verify the effectiveness of the proposed method, several dynamic cases are analyzed based on the transient stability to compare the suitable location (the line between buses PLei and CB) and other feasible ones (as shown in Table4). The optimal location of the TCSC controller is determined based on the rotor angle oscillation damping of the generator units having the large output power and the active power damping of the transmission line. In addition, the power distribution on 500/220 kV system at maximum and/or minimum load is considered. For maximum load, the transferred power on the 500 kV system from the North and Central to the South is very high. Thus, when the three-phase fault occurs on the lines such as Vungang (VA)-Danang (DN), and PLeiku (PLei)-DiLinh (DL), the system lost serious stability. The three-phase fault scenarios are considered in this study as listed in Table 5 and can be summarized briefly as follows.Table 5
Scenarios on three-phase fault for test of transient stability.
Case number Fault is nearly bus Line outage (1) 500 kV DN 500 kV DN-VA (2) 500 kV Plei 500 kV PLei-DL (3) 500 kV CB 500 kV TD-CB (4) 500 kV NQ 500 kV NQ-TTScenario Number 1. A solid three-phase fault occurs at 1 sec on the line 500 kV DN-VA of circuit 1, close to the bus 500 kV DN, and is cleared after 0.15 sec by tripping the faulted line.Scenario Number 2. A solid three-phase fault occurs at 1 sec on the line 500 kV PLei-DL, close to the bus 500 kV Plei, and is cleared after 0.15 sec by tripping the faulted line.Scenario Number 3. A solid three-phase fault occurs at 1 sec on the line 500 kV TD-CB, close to the bus 500 kV TD, and is cleared after 0.15 sec by tripping the faulted line.Scenario Number 4. A solid three-phase fault occurs at 1 sec on the line 500 kV TT-NQ, close to the bus 500 kV NQ, and is cleared after 0.15 sec by tripping the faulted line.Case 1.
The simulation was done on scenario number 1 based on the relative angle oscillations of generator, supposing that the system has been operating at maximum load in order to compare the difference locations of TCSC (that are on the lines PLei-CB and PLei-TM) and without TCSC. Figure9 plots the relative angle oscillations of generators of Yaly, PhuMy-3, VinhTan, and QuangNinh. Figures 9(a)–9(c) show the TCSC placed in the line PLei-CB, the angle oscillations of generators Yaly, PhuMy-3, and VinhTan are damped out in about 6 sec compared to that in the line PLei-TM and without TCSC controller. In particular, it can be observed from Figure 9(d) that the oscillations of generator QuangNinh are damped out faster than the case when the TCSC is placed in the line PLei-CB.Figure 9
The relative angle oscillations of generators: (a) Yaly; (b) PhuMy-3; (c) VinhTan; (d) QuangNinh.
(a) (b) (c) (d)Case 2.
The simulation of this case was done on scenario number 2 based on the active power of the line and the relative angle oscillations of generators to compare the suitable location of TCSC and other feasible locations as NQ-HT and DN-VA. In this case, the system has been operating at maximum load. Observing from Figure10, TCSC is placed in the line PLei-CB; the transient response is significantly improved compared with that in the two another locations. This could easily be construed from Table 4: the total energy value when TCSC is placed in the line PLei-CB is the larger.Figure 10
The transient response: (a) the relative angle oscillations of generator VinhTan; (b) the active power oscillations in the line CauBong-DucHoa; (c) the active power oscillations in the line Pleilu-ThanhMy.
(a) (b) (c)Case 3.
From Table4, it can be seen that the other feasible locations (DakN-CB and PLei-DL) have the second and third highest total energy values, respectively. The third simulation was carried out for scenario number 3 to compare these feasible locations with the suitable location (PLei-CB). In this case, the system is considered as operating at maximum load. Figure 11 plots the system response as the relative angle oscillations of generators (VinhTan and DuyenHai) and the active power oscillations in the lines (HaTinh-VungAng and PLeilu-ThanhMy). It can be observed from Figures 11(a)–11(d) that with TCSC placed in the line between PLei and CB the oscillations are damped faster compared with two other instances. Evidently, the effect is obtained in terms of oscillation damping from other two feasible locations to be the same because they have the same value of total energy value. Figure 11(e) shows the response of TCSC controller; it can be observed in this figure that the TCSC has injected the reactance to the grid during the period of fault and that the best influence is obtained in terms of reactive power flow output of TCSC when it is placed in the line PLei-CB.Figure 11
The transient response: (a) the relative angle oscillations of generator VinhTan (b); the relative angle oscillations of generator DuyenHai; (c) the active power oscillations in the line HaTinh-VungAng; (d) the active power oscillations in the line Pleiku-ThanhMy; (e) the response of TCSC controller.
(a) (b) (c) (d) (e)Case 4.
Observing from Table4, the highest total energy value was obtained on the electric power supply system in the Northern region corresponding to the TCSC controller placed in the line NQ-HT. Contrarily, in this controller placed in the line Plei-CB, the highest total energy value was obtained on the electric power supply system in the Southern region. This simulation case was done on scenario number 4 to compare two locations having the highest total energy value at two regions. In this case, we suppose that the system has been operating at maximum load. Figure 12 plots the relative angle oscillations of generators QuangNinh and VinhTan. Evidently, it can be observed from this figure that TCSC is placed in line PLei-CB; oscillations are damped faster compared to placement in the line NQ-HT despite the far fault location (in the line NQ-TT) compared with the TCSC placed at the Southern region (in the line PLei-CB).Figure 12
The relative angle oscillations of generator: (a) QuangNinh; (b) VinhTan.
(a) (b)Case 5.
Also, for comparison of the other locations, the simulation for this case was done on scenario number 2, supposing that the system has been operating at minimum load. Figure13 plots the transient response of the system. For maximum load, the transferred power on the 500 kV system from the North region to the South region is very high. Consequently, the amplitude deviation of the generator rotor angle and active power oscillations is very large. However, for minimum load, this deviation is negligible but it is still accordant with the proposed method. It could be seen from Figure 13 that the best effect is obtained in terms of oscillation damping of the generator rotor angle of SonLa and the active power in the lines of NhoQuan-SonLa and HaTinh-VungAng.Figure 13
The transient response: (a) the relative angle oscillations of generator SonLa; (b) the active power oscillations in the line NhoQuan-SonLa; (c) the active power oscillations in the line HaTinh-VungAng.
(a) (b) (c)


## 4.1. Determine Feasible Locations

The correctness of the proposed method for the optimal location of TCSC is verified on the small-signal analysis; the dynamic model of TCSC is chosen as shown in Figure4 [34]. We suppose the compensation capacity of TCSC is set to the range [0.75, 1.2] of the uncompensated line. The active power perturbation signal in the line is considered as the input signal for the TCSC controller [51]. The optimal method in [38] was used to determine the parameters for the TCSC controller, and these parameters are listed in Table 2 [34].Table 2
Parameters of the TCSC controller.
Parameter Value Parameter Value T 1 0.1 sec K 0.75 T 2 0.1 sec X m a x 1.2Xline p.u. T 3 0.4 sec X m i n 0.25Xline p.u. T w 10 sec T t c s c 0.015 secThe studied contingency cases of active power perturbation were selected from single-line outage cases on the basis of the real power flow performance index (PI) introduced in [36], as follows:(43)PI=∑i=1nlυi2ϑPliPlimax2ϑ,where Pli and Plimax are the real power and rated power flow of the line i, respectively, ϑ is the exponent, and υi is a the real nonnegative weighting factor for line i. In this study, the value of the exponent was selected as 2.0 and 1.0 is the value of υi.Observing from (43) that all the power flows in the lines are within their limits, PI is small and, vice versa, it reaches out a high value when there are overloads. The studied contingency cases are listed in Table 3 by using (43); this corresponds to the active power perturbation signal in the line of the Vietnamese network.Table 3
Studied contingency cases.
Case number Active power perturbation signal in the line between buses (1) NhoQuan-HoaBinh (NQ-HB) (2) VietTri-Pitoong (VT-PT) (3) ThangLong-PhoNoi (TL-PN) (4) NhoQuan-SonLa (NQ-SL) (5) Pitoong-HoaBinh (PT-HB) (6) QuangNinh-HiepHoa (QN-HH) (7) HiepHoa-VietTri (HH-VT) (8) HiepHoa-Pitoong (HH-PT) (9) Pitoong-SonLa (PT-SL) (10) QuangNinh-ThangLong (QN-TL) (11) QuangNinh-PhoNoi (QN-PN) (12) QuangNinh-MongDuong (QN-MD) (13) ThuongTin-Nho Quan (TT-NQ) (14) ThuongTin-PhoNoi (TT-PN) (15) VungAng-HaTinh (VA-HT) (16) DaNang-ThanhMy (DN-TM) (17) SongMay-TanUyen (SM-TY) (18) SongMay-TanDinh (SM-TD) (19) NhaBe-PhuLam (NB-PL) (20) TanDinh-CauBong (TD-CB) (21) PhuLam-MyTho (PL-MT) (22) PhuLam-DucHoa (PL-DH) (23) DucHoa-MyTho (DH-MT) (24) MyTho-DuyenHai (MT-DH) (25) OMon-MyTho (OM-MT) (26) VinhTan-SongMay (VT-SM) (27) PhuMy-NhaBe (PM-NB) (28) CauBong-PhuLam (CB-PL) (29) NhaBe-MyTho (NB-MT) (30) CauBong-DucHoa (CB-DH) (31) SongMay-PhuMy (SM-PM)The Vietnamese 500/220 kV power system created the state matrix (A) having dimension [2371 × 2371] that corresponds with 179 generator units. By performing the equivalence of 179 generator units into 45 equivalent generator units, the state matrix now has dimension [515 × 515]. A TCSC controller has three states as the transducer, washout, and lead/lag added in the system. Therefore, the state matrix dimension of the equivalent power system model is [518 × 518].Once again, the state matrix of the system is reduced the dimension by applying the balanced order reduction method to eliminate the singular values that are smaller than 10−3. Figure 6 shows the eigenvalue of the system state matrix before and after using balanced reduction. Therefore, the state matrix (A11) of the new system has the dimension [45 × 45].Figure 6
Eigenvalue distribution before and after using balanced reduction.Figure7 shows the frequency response of the system before and after using balanced reduction. Obviously, the input-output behaviors of the system in both cases at the bandwidth range from 10−∞ to 105 are the same with the same control signal. The bandwidth has a range from 105 to 10∞; the frequency response of the reduced system is flat (a straight line) since its order is smaller than the original system. Therefore, we can conclude that two systems are equivalent in terms of the input-output behaviors.Figure 7
Frequency response of the original and reduced system.Figure8 shows the Hankel singular values after performing the balanced realization reduction technique. It can be observed from this figure that the 45 Hankel singular values have been finally retained in comparison with the 2374 state variables of the original system (including state variables of TCSC controller). This acquired singular values are smaller than the chosen threshold 10−3 (choice of this study). The eigenvalues distribution of the system before and after using balanced reduction is plotted in Figure 8. Therefore, the reduced system has 45 state variables, which can fully meet effects of the network to calculate the Gramian matrices.Figure 8
Distribution of Hankel singular values.The controllability Gramian energy indices were determined based on the proposed algorithm as shown in Section3.3 and the active power perturbation signals in the line as listed in Table 3 when the transient horizon is equal to infinity (T=∞) are given in Table 4, in which we only give several feasible locations having the high total energy value. From this table, it can be seen that the line PLei-CB is the suitable location for TCSC controller because it has the maximum total energy value; it means that the energy needed to drive the controllable state variables is smaller than other cases. Therefore, the line between buses PLei and CB is considered as the best location to install the TCSC controller.Table 4
Energy values according to several feasible locations of TCSC.
Case number Contingency cases (from Table3) (Trace ofGC) × 105 TCSC is placed on the line between buses None NQ-HT DN-VA PLei-TM PLei-DS DN-DS PLei-DL PLei-CB PLei-DakN DakN-CB DL-TD The Northern region (  1) (1) NQ-HB 0.0530 0.1354 0.1189 0.0780 0.0516 0.0614 0.0693 0.0780 0.0626 0.0290 0.0381 (2) VT-PT 0.0020 0.0142 0.0132 0.0020 0.0067 0.0080 0.0103 0.0118 0.0089 0.0055 0.0064 (3) TL-PN 0.0012 0.0947 0.0873 0.0374 0.0368 0.0483 0.0506 0.0526 0.0450 0.0215 0.0254 (4) NQ-SL 0.0242 0.0923 0.0840 0.0520 0.0354 0.0429 0.0458 0.0502 0.0418 0.0185 0.0242 (5) PT-HB 0.0034 0.0606 0.0575 0.0120 0.0240 0.0294 0.0307 0.0326 0.0282 0.0122 0.0156 (6) QN-HH 0.0010 0.0108 0.0099 0.0011 0.0044 0.0043 0.0048 0.0059 0.0047 0.0029 0.0032 (7) HH-VT 0.0012 0.0095 0.0097 0.0025 0.0044 0.0056 0.0061 0.0059 0.0054 0.0030 0.0031 (8) HH-PT 0.0032 0.0129 0.0122 0.0040 0.0061 0.0074 0.0092 0.0101 0.0079 0.0060 0.0055 (9) PT-SL 0.0010 0.0137 0.0117 0.0064 0.0059 0.0031 0.0074 0.0109 0.0071 0.0043 0.0057 (10) QN-TL 0.0029 0.0613 0.0564 0.0135 0.0242 0.0307 0.0337 0.0364 0.0301 0.0156 0.0176 (11) QN-PN 0.0032 0.0957 0.0887 0.0129 0.0372 0.0489 0.0514 0.0530 0.0458 0.0223 0.0253 (12) QN-MD 0.0216 0.0361 0.0343 0.0160 0.0140 0.0197 0.0194 0.0188 0.0171 0.0082 0.0090 (13) TT-NQ 0.0054 0.1523 0.1368 0.0850 0.0600 0.0766 0.0864 0.0941 0.0763 0.0382 0.0450 (14) TT-PN 0.0152 0.1173 0.1072 0.0630 0.0460 0.0598 0.0646 0.0685 0.0574 0.0279 0.0330 Total energy of (  1) 0.1385 0.9067 0.8278 0.3858 0.3569 0.4461 0.4897 0.5288 0.4383 0.2151 0.2571  The Southern region (  2) (15) VA-HT 0.0182 0.2777 0.2887 0.1660 0.1227 0.1616 0.1759 0.1900 0.1546 0.0780 0.0930 (16) DN-TM 0.0287 0.2015 0.4708 0.4982 0.5948 0.1894 0.3267 0.3003 0.2733 0.1642 0.1462 (17) SM-TY 0.0019 0.0037 0.0035 0.0002 0.0015 0.0022 0.0656 0.0125 0.0132 0.0296 0.0603 (18) SM-TD 0.0085 0.0453 0.0646 0.0399 0.0412 0.0416 0.3961 0.1184 0.0281 0.1382 0.3472 (19) NP-PL 0.0043 0.0683 0.0829 0.0492 0.0502 0.0460 0.1330 0.3648 0.1638 0.3823 0.0413 (20) TD-CB 0.0039 0.0268 0.0380 0.0359 0.0248 0.0232 0.1742 0.1055 0.5285 0.1123 0.1103 (21) PL-MT 0.0029 0.0158 0.0253 0.0229 0.0167 0.0188 0.0607 0.1077 0.0451 0.1217 0.0498 (22) PL-DH 0.0013 0.0079 0.0062 0.0013 0.0030 0.0021 0.0381 0.0175 0.0092 0.0255 0.0262 (23) DH-MT 0.0048 0.0180 0.0273 0.0095 0.0177 0.0189 0.0477 0.1083 0.0492 0.1166 0.0413 (24) MT-DHa 0.0018 0.0077 0.0174 0.0136 0.0117 0.0133 0.0425 0.0612 0.0269 0.0659 0.0312 (25) OM-MT 0.0023 0.0176 0.0256 0.0271 0.0166 0.0177 0.0672 0.0930 0.0419 0.1001 0.0498 (26) VT-SM 0.0082 0.0120 0.0289 0.0124 0.0195 0.0198 0.0969 0.0826 0.0343 0.0803 0.0911 (27) PM-NB 0.0011 0.0370 0.0341 0.0340 0.0184 0.0111 0.1370 0.1166 0.0620 0.0931 0.0798 (28) CB-PL 0.0043 0.0815 0.1142 0.0150 0.0734 0.0774 0.0584 0.5821 0.2359 0.6319 0.1125 (29) NB-MT 0.0034 0.0078 0.0114 0.0100 0.0077 0.0089 0.0774 0.0357 0.0123 0.0391 0.0502 (30) CB-DH 0.0150 0.0346 0.0499 0.0380 0.0330 0.0357 0.0347 0.2696 0.1048 0.3000 0.0600 (31) SM-PM 0.0112 0.0413 0.0424 0.0121 0.0245 0.0202 0.2338 0.0395 0.0303 0.0146 0.1432 Total energy of (  2) 0.1218 0.9043 1.3312 0.9853 1.0773 0.7079 2.1659 2.6053 1.8134 2.4934 1.5332  Total energy of (  1) + (  2) 0.2604 1.8111 2.1590 1.3711 1.4342 1.1540 2.6556 3.1341 2.2517 2.7085 1.7903


## 4.2. Retest Transient Stability

To verify the effectiveness of the proposed method, several dynamic cases are analyzed based on the transient stability to compare the suitable location (the line between buses PLei and CB) and other feasible ones (as shown in Table4). The optimal location of the TCSC controller is determined based on the rotor angle oscillation damping of the generator units having the large output power and the active power damping of the transmission line. In addition, the power distribution on 500/220 kV system at maximum and/or minimum load is considered. For maximum load, the transferred power on the 500 kV system from the North and Central to the South is very high. Thus, when the three-phase fault occurs on the lines such as Vungang (VA)-Danang (DN), and PLeiku (PLei)-DiLinh (DL), the system lost serious stability. The three-phase fault scenarios are considered in this study as listed in Table 5 and can be summarized briefly as follows.Table 5
Scenarios on three-phase fault for test of transient stability.
Case number Fault is nearly bus Line outage (1) 500 kV DN 500 kV DN-VA (2) 500 kV Plei 500 kV PLei-DL (3) 500 kV CB 500 kV TD-CB (4) 500 kV NQ 500 kV NQ-TTScenario Number 1. A solid three-phase fault occurs at 1 sec on the line 500 kV DN-VA of circuit 1, close to the bus 500 kV DN, and is cleared after 0.15 sec by tripping the faulted line.Scenario Number 2. A solid three-phase fault occurs at 1 sec on the line 500 kV PLei-DL, close to the bus 500 kV Plei, and is cleared after 0.15 sec by tripping the faulted line.Scenario Number 3. A solid three-phase fault occurs at 1 sec on the line 500 kV TD-CB, close to the bus 500 kV TD, and is cleared after 0.15 sec by tripping the faulted line.Scenario Number 4. A solid three-phase fault occurs at 1 sec on the line 500 kV TT-NQ, close to the bus 500 kV NQ, and is cleared after 0.15 sec by tripping the faulted line.Case 1.
The simulation was done on scenario number 1 based on the relative angle oscillations of generator, supposing that the system has been operating at maximum load in order to compare the difference locations of TCSC (that are on the lines PLei-CB and PLei-TM) and without TCSC. Figure9 plots the relative angle oscillations of generators of Yaly, PhuMy-3, VinhTan, and QuangNinh. Figures 9(a)–9(c) show the TCSC placed in the line PLei-CB, the angle oscillations of generators Yaly, PhuMy-3, and VinhTan are damped out in about 6 sec compared to that in the line PLei-TM and without TCSC controller. In particular, it can be observed from Figure 9(d) that the oscillations of generator QuangNinh are damped out faster than the case when the TCSC is placed in the line PLei-CB.Figure 9
The relative angle oscillations of generators: (a) Yaly; (b) PhuMy-3; (c) VinhTan; (d) QuangNinh.
(a) (b) (c) (d)Case 2.
The simulation of this case was done on scenario number 2 based on the active power of the line and the relative angle oscillations of generators to compare the suitable location of TCSC and other feasible locations as NQ-HT and DN-VA. In this case, the system has been operating at maximum load. Observing from Figure10, TCSC is placed in the line PLei-CB; the transient response is significantly improved compared with that in the two another locations. This could easily be construed from Table 4: the total energy value when TCSC is placed in the line PLei-CB is the larger.Figure 10
The transient response: (a) the relative angle oscillations of generator VinhTan; (b) the active power oscillations in the line CauBong-DucHoa; (c) the active power oscillations in the line Pleilu-ThanhMy.
(a) (b) (c)Case 3.
From Table4, it can be seen that the other feasible locations (DakN-CB and PLei-DL) have the second and third highest total energy values, respectively. The third simulation was carried out for scenario number 3 to compare these feasible locations with the suitable location (PLei-CB). In this case, the system is considered as operating at maximum load. Figure 11 plots the system response as the relative angle oscillations of generators (VinhTan and DuyenHai) and the active power oscillations in the lines (HaTinh-VungAng and PLeilu-ThanhMy). It can be observed from Figures 11(a)–11(d) that with TCSC placed in the line between PLei and CB the oscillations are damped faster compared with two other instances. Evidently, the effect is obtained in terms of oscillation damping from other two feasible locations to be the same because they have the same value of total energy value. Figure 11(e) shows the response of TCSC controller; it can be observed in this figure that the TCSC has injected the reactance to the grid during the period of fault and that the best influence is obtained in terms of reactive power flow output of TCSC when it is placed in the line PLei-CB.Figure 11
The transient response: (a) the relative angle oscillations of generator VinhTan (b); the relative angle oscillations of generator DuyenHai; (c) the active power oscillations in the line HaTinh-VungAng; (d) the active power oscillations in the line Pleiku-ThanhMy; (e) the response of TCSC controller.
(a) (b) (c) (d) (e)Case 4.
Observing from Table4, the highest total energy value was obtained on the electric power supply system in the Northern region corresponding to the TCSC controller placed in the line NQ-HT. Contrarily, in this controller placed in the line Plei-CB, the highest total energy value was obtained on the electric power supply system in the Southern region. This simulation case was done on scenario number 4 to compare two locations having the highest total energy value at two regions. In this case, we suppose that the system has been operating at maximum load. Figure 12 plots the relative angle oscillations of generators QuangNinh and VinhTan. Evidently, it can be observed from this figure that TCSC is placed in line PLei-CB; oscillations are damped faster compared to placement in the line NQ-HT despite the far fault location (in the line NQ-TT) compared with the TCSC placed at the Southern region (in the line PLei-CB).Figure 12
The relative angle oscillations of generator: (a) QuangNinh; (b) VinhTan.
(a) (b)Case 5.
Also, for comparison of the other locations, the simulation for this case was done on scenario number 2, supposing that the system has been operating at minimum load. Figure13 plots the transient response of the system. For maximum load, the transferred power on the 500 kV system from the North region to the South region is very high. Consequently, the amplitude deviation of the generator rotor angle and active power oscillations is very large. However, for minimum load, this deviation is negligible but it is still accordant with the proposed method. It could be seen from Figure 13 that the best effect is obtained in terms of oscillation damping of the generator rotor angle of SonLa and the active power in the lines of NhoQuan-SonLa and HaTinh-VungAng.Figure 13
The transient response: (a) the relative angle oscillations of generator SonLa; (b) the active power oscillations in the line NhoQuan-SonLa; (c) the active power oscillations in the line HaTinh-VungAng.
(a) (b) (c)


## 5. Conclusions

In this paper, a relevant stochastic method for the optimal placement of TCSC controller has been presented to enhance the rotor angle stability and dampen the power system oscillations in the multimachine systems. This proposed method is developed from the energy approach based on the controllability and observability Gramian matrices of the linearized multimachine systems. The optimal placement depends on the trace indices of the Gramian matrices that have been calculated on the active power perturbation in the line of the network (for this study, the applied network is the Vietnamese 500/220 kV power system).The optimal placement for TCSC controller is determined based on the Gramian critical energy values that have been calculated on the small-signal stability analysis. However, the acquired results showed that the power system could operate perfectly under the influence of the transient conditions.The time-domain simulation results on the transient stability analysis show that the rotor angle oscillations of generator and the power oscillations in the line are significantly dampened when the TCSC is placed in the line between PLeiku and CauBong, which has the maximum total energy value.The Gramian-based reduction method has been also introduced for the purpose of reducing the calculation time of the Gramian critical energy when dealing with the large-scale power systems.



---

*Source: 3014510-2017-01-29.xml*