Performance and security enhancement using generalized optimal unified power flow controller under contingency conditions and renewable energy penetrations

In this paper, a novel flexible AC transmission system (FACTS) device named generalized optimal unified power flow controller (GOUPFC) is introduced to control the power flows in multi transmission lines and to regulate the voltages and angles at the load buses. The detailed power injection modeling of GOUPFC is presented in this paper. The optimal location of GOUPFC is determined based on line collapse proximity indicator (LCPI). A multi-objective function is framed in terms of average voltage deviation (AVDI), real power loss (Ploss) and average line collapse proximity indicator (LCPIavg) to test the effectiveness of the proposed device. The simulation studies are performed on standard IEEE 57-bus test system under single line contingency and considering various renewable energy source (RES) penetrations. The control parameters of GOUPFC are optimized by using whale optimization (WO-BAT) algorithm, by hybridizing WOA and BAT algorithms, and the superiority of WO-BAT is observed in minimizing the proposed objective function and enhancing the voltage profile.

loadability with OUPFC device and is presented in [22]. The economic operation is performed in power systems to analyze the effect of OUPFC under weather changing conditions as well as variable loading conditions [23]. Transmission system security is enhanced under single line contingency conditions with OUPFC device under different renewable energy generations and is proposed in [24].
In light of the reviewed works, we claim that the following are the major contributions of this paper.
1. For the first time, a novel and advanced FACTS device, namely GOUPFC and its steady-state modeling, is introduced. 2. Besides, the optimal location and sizing of GOUPFC are determined using a novel hybrid approach WO-BAT. 3. In order to improve the exploitation features of WOA, a predefined search space for locations is defined using line collapse proximity indicator (LCPI) and its convergence features are improved using BAT algorithm. 4. The impact of GOUPFC is analyzed on IEEE 57-bus system considering line contingencies and renewable energy uncertainty.

Methods
It is possible to extend the voltage and power flow control beyond what can achievable with the OUPFC by using a new configuration is called as generalized optimal unified power flow controller (GOUPFC). To the best of our knowledge, no research work has been developed in the area of GOUPFC power injection modeling and its application to achieve the optimal operation of the power system. Firstly, the power injection modeling of GOUPFC device is presented by considering switching losses, and then, the optimal location is determined based on LCPI index. The control parameters of GOUPFC device are optimized with WO-BAT, WOA and BAT algorithms. Besides, the performance of GOUPFC is investigated under single line contingencies and with support of various renewable penetrations for standard IEEE-57 test system to minimize the objective function formulated in terms of AVDI, P loss and LCPI avg and to enhance voltage profile. This paper is organized in the following sections. Section 3 describes the operating principle of GOUPFC. Section 4 presents the mathematical modeling of GOUPFC. Section 5 presents the overall objective function formulation. Section 6 presents WO-BAT algorithm description, and Sect. 7 presents the determination of GOUPFC optimal location and the analytical results evaluated from different case studies.

Operating principle of GOUPFC
The transmission system security in the multi-lines can be enhanced by including another OUPFC in different line is not a cost-effective solution instead going for the new device named GOUPFC which is the best choice. The single line diagram of GOUPFC is shown in Fig. 1. It consists of one exciting transformer with four windings such as one primary winding and three secondary windings named as secondary, tertiary and quaternary and two triple winding injecting transformers each with one primary and two secondaries named as secondary and tertiary. The PST in one line is coupled to the two secondary windings of the exciting and injecting transformers, and the PST in the other line is connected to the quaternary and secondary windings of the exciting and injecting transformer, respectively. One of the three converters is connected in shunt with a bus, and the remaining two converters are connected in series with two different transmission lines through the tertiary winding of the exciting and injecting transformers. The overall configuration of GOUPFC can able to control the total six control parameters which includes the controllable voltage magnitude and phase angles at bus-í and independent real and reactive power flows in the two lines. Figure 1 shows the two transmission lines i-j and i-k is connected with GOUPFC at buses j and k, respectively.
The series converter and PST will inject a controllable voltage magnitude and phase in series into the lines through injecting transformers. The shunt converter can able to (1) supply or absorb the real power needed by the series converter; (2) exchange the controllable reactive power with line; and (3) regulate the dc link voltage.

Mathematical modeling of GOUPFC
A GOUPFC can be represented by three controllable voltage source converters, and two-phase shifting transformers connected in two transmission lines is shown in Fig. 2. The series converter controllable voltages and PST voltages are given by, where r 1 , r 2 and k 1 , k 2 are the per unit voltage magnitudes of two series converters and two PSTs, respectively;γ 1 and γ 2 are the phase angles of series converters; σ 1 and σ 2 are the PST phase angles; and V i is the bus-i voltage. The voltages and phase angles are operating in the limits specified as follows.
GOUPFC having three voltage source converters among those one is shunt converter, and the remaining two are series converters. Shunt converter is placed at bus-I, and the two series converters are placed in the lines i-j and i-k, respectively. In addition to the voltage source converters, GOUPFC employing two separate PSTs which are incorporated in the two separate lines where the two series converters are placed. The total voltages injected in the two lines are the phasor sum of the voltages obtained from the series converters and PSTs, and these are given by, The GOUPFC is modeled in series-connected and shunt-connected voltage source models. The voltages in the two transmission lines behind the line reactance can be written mathematically as,

Series-connected voltage source model
The series-connected voltage source model of GOUPFC with three buses i, j, and k is shown in Fig. 2a, and the corresponding Norton's equivalent current source model is shown in Fig. 2b. Let the currents flowing in the two lines i-j and i-k are I se 1 and I se 2 , respectively, are given by where X se 1 and X se 2 are the line reactances of the two lines and their corresponding susceptances, respectively, given by, B se 1 = 1/X se 1 and B se 2 = 1/X se 2 .
The independent complex power injections of GOUPFC at buses i, j and k are given as, Let θ ij = θ i − θ j ; θ ik = θ i − θ k . By using Euler's and trigonometric identities, the real and reactive power injections at the buses i, j and k are calculated as follows.
The equivalent power injection modeling for the series-connected voltage source model is shown in Fig. 3a. The amount of complex power supplied by the combination of series converter and PST in the individual lines is derived as follows The real and reactive powers supplied by the two converters in the two transmission lines are given by,

Shunt-connected voltage source model
The equivalent circuit for the shunt-connected voltage source model is shown in Fig. 3b. This model provides the equivalent power injections at the GOUPFC shunt bus to the two series branches through the converter and PST combination. The voltage at the sending end is controlled by the reactive power injection at the GOUPFC shunt converter. The amount of real power supplied by the shunt converter is equal to the real power consumed by the two series converters. Therefore, the real power inserted at the shunt converter is given by: Let assume a constant reactive power injection Q sh at bus-i and the apparent power injection at the shunt bus-i is given as

Final GOUPFC modeling
The final PIM of GOUPFC is achieved by summing up the equations obtained in series-and shunt-connected models. The corresponding equivalent circuit representing GOUPFC power injection is shown in Fig. 4. The resultant real and reactive power equations at the GOUPFC buses are given as

Overall objective function (OOF) formulation
An OPF-based multi-objective function is formulated to minimize the overall objective function (OOF) is given as, where w 1 , w 2 and w 3 are the weighing factors considered as one. f 1 , f 2 and f 3 are the AVDI, P loss and average LCPI, respectively.

Average voltage deviation index (AVDI)
The voltage profile of the test system subjected to contingency condition is evaluated in terms of stability index called AVDI. It is given as [23], where |V i,ref | is the reference voltage at bus-i, |V i | is the voltage at bus-i under contingencies and nb is the number of load buses. In specific, AVDI is used to explain the average voltage deviation of the entire network in comparison with the reference bus voltage (i.e., maximum of all bus voltages). Under each contingency, the voltage profile may change, and thus, there is a change in AVDI. The more AVDI means, more voltage imbalance in the network.

Real power loss (P loss )
The real power loss is given mathematically as, where nb is no. of buses; nl is no. of lines; k is the line number; r k and I k are the resistance of line k and current through it, respectively; Y ii and θ ij are the shunt admittance at bus i and its angle, respectively; V i and δ i are the voltage magnitude and its angle at bus-i, respectively; V j and δ j are the voltage and its angle at busj, respectively.

Line collapse proximity index (LCPI)
The error in evaluating the voltage stability due to neglecting the line charging reactance and the magnitude and relative directions of active and reactive power can be overcome by considering an index called LCPI [25]. LCPI is given as, where A and B are network parameters magnitude; α and β are the phases of A and B, respectively; and V i is the sending end voltage.P j and Q j are the real and reactive power flow; for a line, LCPI should always less than 1 failing leads to instability.

Constraints
The OOF expressed in Eq. (25) is subjected to various equality and inequality constraints [21] as given below.

Equality constraints
The active and reactive power balance equations are the equality constraints that can be expressed, for all the buses expect FACTS incident buses, as The active and reactive power balance equations at the FACTS incident buses are given as,

Inequality constraints
The inequality constraints considered for optimization problem are given as follows.

WO-BAT algorithm
In this paper, the control parameters of GOUPFC device are optimized with the hybrid algorithm WO-BAT. The detailed description about WO-BAT is described as follows.

Whale optimization algorithm
The exploitation and exploration are the two stages involved in WOA [26]. In exploitation stage, surrounding prey and spiral updation of the position vector are modeled. This method is called bubble net attacking (BNA) method [27]. Random searching of prey is performed in exploitation stage. After identifying the position of prey, humpback whales surround them. Initially, the optimal location of prey in the search space is not defined. Therefore, this algorithm considers the present solution which is the optimal prey and the other search delegates (agents) will drifted toward the best search delegates. This can be expressed mathematically as, In the exploration stage, to compel the search delegates to move far away from the local whale, the coefficient vector − → C 1 is used to generate random values less or greater than one. The location of delegates is identified based on this random selection rather than the best search delegate. This gives the global solution by overcoming the location solution. This can be given mathematically, where −−→ P rand is the random position vector.

Bat algorithm
Bats are the creatures with echo-location abilities. They create a loud sound pulse and receive the echo from the neighboring objects. They can guess the position of the neighboring object by using the time delay of echo sound. They measure the shape and the direction of object by sound pulse comparative amplitude analysis received at the ear. They investigate and simplify the data collected and figured out an image in brain to identify the neighboring object. The concept of bat algorithm (BA) and its mathematical modeling is provided in [28].
To determine the location of prey, bats fly randomly in the search space with a velocity v i . Later, they change their positions (x i ) with constant frequency (f min ), different wavelengths (β) and loudness A 0 . The new solutions for position and velocity are given mathematically as, where β is a random number obtained by uniform distribution in the range [0, 1]. x* is the present global best; the minimum and maximum frequency limits (f min and f max ) are taken 0 and 100, respectively.

WO-BAT algorithm
The major limitation of WOA is its convergence speed in obtaining global solution. However, this can be overcome by embedding BAT algorithm partially to WOA to increase the exploration. In this hybrid algorithm, the global position is updated based on condition technique, which means if the present solution is better than the old solution, then the old one can be replaced with the present solution. The detailed flowchart for this hybrid algorithm is shown in Fig. 5.

Results and discussions
The performance of GOUPFC device is analyzed under single line contingencies, and an OPF-based objective function is solved by using WO-BAT, WOA and BAT algorithms under the different RES penetrations. The test system data are taken from [29]. In comparison with earlier works [23,24], this work is focused on the development of mathematical modeling for the novel GOUPFC and its location and parameters optimization using novel hybrid approach WO-BAT in this work. The type of RES considered in this work is photovoltaic system, and its mathematical modeling is adapted from Ref. [23]. The simulations are carries under five different cases:

Optimal location of GOUPFC
Under base case conditions (i.e., no GOUPFC device and no RES), LCPI is calculated for all the lines of IEEE-57 test system and arranged in descending order as shown in Table 1. The top 5 ranked locations are only shown here after excluding the lines incident to generator and tap setting transformer. From Table 1, the two lines 38-49 (#78) and 37-38 (#50) are chosen for the placement of GOUPFC device. The shunt converter is placed at the bus#38 which is the common bus to the two line and the two series converters are placed in the two lines.

IEEE-57 bus description
IEEE-57 bus test system has 7 PV buses (i.e., 1, 2, 3, 6, 8, 9 and 12) and 50 PQ buses interconnected by 80 transmission lines. It has real and reactive power loads, respectively,  Fig. 6. From this figure, under no contingency conditions, GOUPFC along with 50% RES penetration level gives minimum OOF value for WO-BAT algorithm. In other words, as RES penetration increases the net effective loading on the system decreases. Thus, the OOF is decreased. However, by having GOPUF controls, it is further decreases. The optimal controls derived with WO-BAT are caused to improve the system performance than basic WOA and BAT, significantly.
In case2, the system is suffering with single line contingency and there is no support from GOUPFC device and RES generation. Under these conditions, the performance indices and OOF for no contingency and one line contingency are presented in Table 2. The system has real power loss and LCPI value with no contingency 27.8638 MW and 0.081, respectively, and the objective function value 27.9448. A contingency analysis is performed on the system without GOUPFC and RES, and the indices for the top ranked line 54-55 (#70) are as follows: P loss = 29.7828, ALCPI = 0.0860 and OOF = 29.8688. The individual objective functions and OOF values for the remaining lines are presented in Table 2.
In case 3, the integration of GOUPFC without RES is analysed under one line contingencies. The parameters of the GOUPFC device are optimized WO-BAT, WOA and BAT algorithms, and the results are demonstrated in Table 3. When there is no outage the real power loss is decreased to 22.9129 MW which is 17.77% decrement with reference to case 1. Under single line 54-55 (#70) contingency also GOUPFC has shown its superiority in terms of losses and LCPI value. The indices in one line contingency are P loss = 25.1792, ALSI = 0.0467, ALCPI = 0.0747 and OOF = 25.2539. Here, the real power loss and OOF have been reduced to 34.18% and 51.83%, respectively, with reference to case 1. The LCPI and OOF values for other lines are presented in Table 3. The  performance of GOUPFC under one line contingency at without RES for three algorithms is shown in Fig. 7.
In case 4, GOUPFC and 30% power penetration are together considered for the contingency. Here, the RES installed capacity is 30% of the real power load, i.e., 1250.8 × 0.3 = 375.24 MW. Therefore, the remaining load on the CE sources is 875.56 MW. Due to the integration of GOUPFC and RES into the system, the line flows and hence transmission losses are reduced under without and with line contingency conditions. The performance indices are evaluated for the both the conditions, and the losses are reduced by 37.57% and 34.18%, respectively. The optimized OOF values by WO-BAT algorithm are 17.4563 and 19.6639, respectively, for without and with line (54-55) contingency conditions. The individual objective functions and OOF for the remaining line contingencies are presented in Table 4. The performance of GOUPFC under single line contingency with0.3penetrations is shown in Fig. 8.
In case 5, the RES generation is enhanced to 50% and the support of GOUPFC device is still present. Here the RES installed capacity is 50% of the real power load, i.e., 1250.8 × 0.5 = 625.4 MW. Hence the remaining load on the CE sources is to be 625.4 MW. The control variables of GOUPFC device are optimized with three algorithms, and the results of WO-BAT are demonstrated here due to its superiority. Under no contingency condition, the performance indices with WO-BAT algorithm are P loss = 14.1166, ALCPI = 0.0566 and the OOF = 14.2489. Here, the power loss and OOF are diminished by 49.34% and 96.90%, respectively, with respect to case 1. Under contingency conditions, the line 54-55 (#70) with performance indices as P loss = 15.2089, ALCPI = 0.0506 and the OOF = 15.3261. Here also the combination of GOUPFC and 50% RES penetration level reduces the real power losses and OOF by 48.93% and 95.64%, respectively, with respect to case 1. The P loss , LCPI avg values and the corresponding OOF for the remaining line outages are presented in Table 5. The performance of GOUPFC under single line contingency with0.5penetrationsis shown in Fig. 9.  The control parameters of GOUPFC device are UPFC voltages ((r 1upfc (p.u), r 2upfc (p.u)) and UPFC phase angles (γ 1upfc (radians), γ 2upfc (radians)), PST ratios (k 1pst (p.u), k 2pst (p.u)) and PST angles (σ 1pst (radians), σ 2pst (radians)) for IEEE-57 bus system is optimized with WO-BAT algorithm and the corresponding parameters in case 5 is presented in Table 6. Furthermore, the comparison of convergence time for the three algorithms is given in Table 7. From this table, it can be observed that the hybrid WO-BAT algorithm is converged much faster than the individual WOA and BAT algorithms in all the case studies. Figure 10 shows the voltage profile at all the buses of IEEE-57 bus under with and without contingency conditions for base case, only with GOUPFC device and due to the presence of both GOUPFC and 50% RES support. From this figure, it is clear that the voltage profile has been improved much with GOUPFC than base case and GOUPFC with RES support than only GOUPFC device present in the system. From the above, it is concluded that the proposed GOUPFC device can improve the voltage stability not only under normal conditions but also under single line contingencies.
As the methodology is implemented on Standard IEEE 57-bus test system, and a significant improvement is observed in all the cases, we believe that this methodology also works on other networks. However, the computational efficiency of any algorithm may not quantifiable only on limited search space with small test systems. Thus, there is a need for evaluating the performance of proposed methodology on larger test systems, which we can consider as our future research work.

Conclusion
The detailed power injection modeling of the proposed FACTS device GOUPFC is presented in this paper. The optimal location of GOUPFC device is determined based on LCPI value. The performance indices and OOF value have been evaluated for different    case studies by considering GOUPFC device at different RES penetrations. The parameters of GOUPFC are optimized by using WO-BAT, WOA and BAT algorithms. With the combination of GOUPFC and increased RES penetrations the AVDI, losses, LCPI avg and OOF values are considerably reduced under no contingency conditions is observed in case 1. Under single line contingency conditions, without RES and GOUPFC, there is an increment in losses, LCPI index and reduced voltage profile is observed in case 2. All the performance indices and OOF value are significantly reduced with optimal placement of GOUPFC and additional support of RES penetrations 0%, 30% and 50% in the cases 3, 4 and 5, respectively. In addition, the voltage profile is also improved for the case GOUPFC with 50% RES level when compared to the other cases. From the simulation results, the hybrid algorithm WO-BAT is performed better than the individual algorithms such as WOA-BAT in minimizing AVDI, P loss , LCPI and OOF value. With the optimal location of GOUPFC, the overall system voltage profile is enhanced in IEEE-57 bus test system. Further, the impact the proposed device can be extended for largescale systems in the real-time scenario by incorporating two or more GOUPFC devices at optimal locations which is considered as the future scope of the present work.

FACTS Flexible AC transmission systems OUPFC
Optimal unified power flow controller GOUPFC Generalized optimal unified power flow controller WO-BAT Whale optimization-BAT LCPI Line collapse point indicator RES Renewable energy sources