**A Neural Network-Based Model Reference Control Architecture for Oscillation Damping in Interconnected Power System**

**Waqar Uddin 1, Nadia Zeb 2, Kamran Zeb 1,3, Muhammad Ishfaq 1, Imran Khan 1, Saif Ul Islam 1, Ayesha Tanoli 4, Aun Haider 4, Hee-Je Kim 1,\* and Gwan-Soo Park <sup>1</sup>**


Received: 12 July 2019; Accepted: 20 September 2019; Published: 24 September 2019

**Abstract:** In this paper, a model reference controller (MRC) based on a neural network (NN) is proposed for damping oscillations in electric power systems. Variation in reactive load, internal or external perturbation/faults, and asynchronization of the connected machine cause oscillations in power systems. If the oscillation is not damped properly, it will lead to a complete collapse of the power system. An MRC base unified power flow controller (UPFC) is proposed to mitigate the oscillations in 2-area, 4-machine interconnected power systems. The MRC controller is using the NN for training, as well as for plant identification. The proposed NN-based MRC controller is capable of damping power oscillations; hence, the system acquires a stable condition. The response of the proposed MRC is compared with the traditionally used proportional integral (PI) controller to validate its performance. The key performance indicator integral square error (ISE) and integral absolute error (IAE) of both controllers is calculated for single phase, two phase, and three phase faults. MATLAB/Simulink is used to implement and simulate the 2-area, 4-machine power system.

**Keywords:** power oscillations; UPFC; non-linear control; neural network; model reference control

#### **1. Introduction**

The electric power system is composed of various interconnected generating stations that are connected to load centers through the transmission and distribution systems. The interconnection of different generating stations increases the complexity of the power network. Due to the stochastic nature of the consumer profile, a mismatch occurs between the load demand and supply. The load variation causes a change in system frequency, voltage, active power, and mechanical power. The variation in these parameters causes fluctuation in the rotor speed. As a result, a low frequency (0.1 Hz to few hertz) oscillation in active power takes place, known as a power oscillation. These power oscillation either decays after a specified time or increases in magnitude until the system collapse (blackout) occurs [1,2].

In order to overcome the issues of stability due to power oscillations and control power flow in a system, a Flexible AC Transmission System (FACTS) has been proposed by scientists and researchers.

FACTS devices serve many functions in a power system, i.e., control of power flow, oscillation damping, voltage control, and improvement of transient stability. The Unified Power Flow Controller (UPFC) is an important member of the FACTS family. The UPFC is composed of two voltage source converters coupled through a DC link. The converters provide both series and shunt compensation to the transmission line, hence separating control of active and reactive power becomes possible.

The problem of Low-Frequency Oscillation (LFO) is tackled through FACTS devices in the modern power system. FACT devices and a UPFC based on FACT devices are used mostly for mitigating the LFOs. These devices are also used for enhancing the stability of the power system. The UPFC ensures the security and stability of the system and damps out the LFOs.

Different linear and nonlinear control techniques are applied to mitigate the LFOs [2]. A modified non-dominated sorting genetic algorithm-II (MNSGA-II) is proposed for control of the UPFC to damp out power oscillation in [3]. The proposed design of the UPFC controller is based on a multi-objective optimization problem. However, the controller response is slow in the case of a heavy loaded condition. In [4], the authors presented Variable Neighborhood Search (VNS) optimization algorithm for tackling the problem of oscillations. Although the design process includes high complexity, the proposed controller responds well during normal conditions. However, the author did not consider the heavy loaded or fault condition. An imperialist competition algorithm is used for damping power oscillation in [5]. The controller is designed for wide range of operating conditions, but the author considered a Single Machine Infinite Bus (SMIB) system. The artificial bee colony (ABC) algorithm [6] and Support Vector Regression (SVR) [7] are also proposed for design of a supplementary damping controller. However, the design is equipped with an additional Power System Stabilizer (PSS), along with the UPFC, which will increase system cost and complexity. The author in [8] proposed adaptive input–output feedback linearization control to achieve a non-linear stability estimator and 4th order linearized model. However, due to the parallel application of these procedures, the system faces instability during disturbances. The linear observer excludes the system nonlinearities and only approximate system parameters. The sub-transient state of the synchronous generator is neglected after the fault clearance, while the FACTS device is turned on and off by using an input–output controller. In [9], the researcher proposed a co-evolutionary cultural algorithm (CPCA) based on particle swarm optimization. The author proposed a PSS, along with UPFC, for oscillation damping. But the response of the controller is slow as it takes too much time in damping the power oscillations. In [10], a PSS is used to stabilize the LFOs. Also, in case of failure of the PSS, a supplementary control based on a Linear Quadratic Regulator (LQR) regulates the control input to FACTS devices. In [11], the author proposed a Fuzzy Neural Non-linear Proportional Integral controller for damping LFOs. The positive aspects of this controller are its adaptive nature in learning rate, online tuning, simple structure, and parameter identification. Inter-area oscillations can be efficiently damped through neural predictive control. It will also enhance the stabilization of the power system [12]. A Takagi-Sugeno controller produces gain over the large operating region. However, it cannot guarantee the stability of the system, the robustness of uncertain perturbation in the system, and parameter variations. Transient stability of a machine connected to an infinite bus system is enhanced by using non-linear coordination of generator excitation in [13].

In our paper, we present the neural network (NN)-based model reference controller (MRC) control scheme with a UPFC for mitigation and damping LFOs that is caused by the fault, disturbances, or due to loss of synchronism in connected machines. The MRC controller uses a NN that has the capability to learn and store information of non-linearities, using the auto-regression model instead of tuning of input parameters. The MRC is the integration of the (1) reference model, (2) NN plant, (3) NN controller, and (4) plant. The incorporation of all these four parts promises the robustness, intelligence, and adaptation of parameters. The MRC has an inherently adaptive mechanism, which will help in tracking of variation of process variables. The ability of a NN to learn and model complex and nonlinear relationships plays a vital role in control of a complex problem with uncertainties. After initial weight adjustment, it acquires the ability to predict uncertainties. These features of an

NN-based MRC controller will tackle the uncertainties and will have a better response than other conventional control schemes. The weight of the NN is adjusted with the Back Propagation (BP) algorithm. Eventually, it will damp out the oscillation quickly and protect the system from complete collapse. Moreover, the performance of proposed controller is compared with the PI controller for the same system to validate our claim.

The paper is organized as follows: Section 2 shows the general system model, while Section 3 discusses the mathematical modelling of the system. The supplementary control and results are discussed in Sections 4 and 5, respectively. Section 6 concludes this paper.

#### **2. System Model**

A UPFC installed, 2-area, 4-generator connected power system is depicted in Figure 1. This model is used to study LFOs of 0.1 to few Hertz.

**Figure 1.** One-line diagram of power system with installed unified power flow controller (UPFC).

A DC capacitor voltage regulator is used to damp out the oscillation produces due to the occurrence of disturbances in power, but it has a negative impact on the generated transients [14] because the required amount of power is transmitted to the series inverter for injection to the transmission line through DC link. Similarly, when variations in magnitude of DC link voltage occurs, it affects power transfer capability of the shunt inverter. Hence, the active power is injected into the system with the help of series inverter [15].

The synchronous generator connected to the system is used to generate power that is delivered through the transmission lines. The installation point of the UPFC is between bus 8 and bus 9.

The DC link capacitor is used to make the back to back connection of Voltage Sources Converters (VSCs) possible and transfer the power from the shunt converter to the series converter [16]. The main purpose of Series VSC (VSC-SE) is to control the injection of active power and reactive power to the system, while Shunt VSC (VSC-SH) is used to adjust active power flow to series converter and maintain DC link voltage. The connection of the shunt and series VSC is depicted in Figure 1. The connection of both the converters is such that it diminishes the total power injection to the system.

$$P\_{\text{S}\text{e}} + P\_{\text{S}\text{h}} = 0.\tag{1}$$

The input control signal to UPFC are *mSe*, *mSh*, δ*Se*, and δ*Sh*. *mSe* and *mSh* is used to represent the amplitude modulation index of series and shunt compensation, respectively. However, δ*Se* and δ*Sh* are used for the representation of series and shunt compensation phase angles [2,17]. In the case of a series converter, δ*Se* and *mSe* are the important factors that need to be determined. Because δ*Se* is used to control real power flow in to the high voltage direct current (HVDC) transmission line, and *mSe* is responsible for voltage control. Similarly, in the case of the shunt converter, the reactive power of the AC terminal can be controlled through *mSh*, and δ*Sh* is used to control DC link voltage (*Vdc*). VSC-SE injects voltage to the AC transmission line through a series connected transformer. This series injected voltage is used to tackle different problems, i.e., phase shift, series compensation, voltage control, and combination of these parameters. VSC-SH is used to take real power from AC terminal and transfer it to DC terminal. This power is then transferred by VSC-SE to the system for compensation. These features make the UPFC a flexible FACTS device that has the properties of all other FACTS devices. The block diagram of the UPFC general control is depicted in Figure 2. *V*<sup>1</sup> and *V*<sup>2</sup> are bus voltages to which VSC-SH and VSC-SE are connected, respectively, while "*i*" is line current.

**Figure 2.** UPFC control.

#### **3. Mathematical Modelling of Power System**

The model used in this paper is composed of a UPFC installed, 4-machine, two-area power system. The UPFC is composed of a VSC − SH, VSC − SE, and DC link.

#### *Mathematical Model*

The non-linear model of the UPFC is presented below from Equation (2) to Equation (3). This model is used for the analysis of stability of the power system. The parameters that affect the dynamics of the system are considered in this model, and resistance and transients of transformer is neglected. The modelling equations given below are transformed to dq − axis by using Park's transformation [17,18]. Series VSC and shunt VSC are modelled as:

$$
\begin{bmatrix} V\_{ds\epsilon} \\ V\_{qs\epsilon} \end{bmatrix} = \begin{bmatrix} 0 & -X\_{\mathsf{s}\epsilon} \\ X\_{\mathsf{s}\epsilon} & 0 \end{bmatrix} \begin{bmatrix} I\_{ds\epsilon} \\ I\_{qs\epsilon} \end{bmatrix} + \frac{m\_{\mathsf{s}\epsilon}V\_{d\epsilon}}{2} \ast \begin{bmatrix} \cos\delta\_{\mathsf{s}\epsilon} \\ \sin\delta\_{\mathsf{s}\epsilon} \end{bmatrix} \tag{2}
$$

*Energies* **2019**, *12*, 3653

$$
\begin{bmatrix} V\_{dsh} \\ V\_{qsh} \end{bmatrix} = \begin{bmatrix} 0 & -X\_{sl} \\ X\_{sh} & 0 \end{bmatrix} \begin{bmatrix} I\_{dsh} \\ I\_{qsh} \end{bmatrix} + \frac{m\_{sl}V\_{dcl}}{2} \ast \begin{bmatrix} \cos \delta\_{sl} \\ \sin \delta\_{sl} \end{bmatrix}.\tag{3}
$$

The subscript *sh*' and *se* is used for shunt and series converter parameters, respectively. *Xse*, *Xsh* are equivalent reactance of series and shunt transformer. *Vdse*, *Vqse* are the dq voltage component of series branch, while *Idse*, *Iqse* are *dq*-component of series converter current. The dynamic of DC voltage at DC link is represented as:

$$\frac{d\upsilon\_{\rm dc}}{dt} = \frac{3m\_{\rm sh}}{4\mathcal{C}\_{\rm dc}} \left[ \begin{array}{cc} \cos\delta\_{\rm sh} & \sin\delta\_{\rm sh} \end{array} \right] \left[ \begin{array}{cc} I\_{\rm shh} \\ I\_{\rm qsh} \end{array} \right] + \frac{3m\_{\rm sc}}{4\mathcal{C}\_{\rm dc}} \left[ \begin{array}{cc} \cos\delta\_{\rm sc} & \sin\delta\_{\rm sc} \end{array} \right] \left[ \begin{array}{cc} I\_{\rm dsc} \\ I\_{\rm qsc} \end{array} \right] . \tag{4}$$

*Vdc* shows the DC link voltage, and *Cdc* is the capacitances of the DC link capacitor. The model of the UPFC connected power system is presented by:

$$
\dot{\delta} = a\rho\_0 \Delta \omega\_\prime \tag{5}
$$

$$
\dot{\omega} = \frac{P\_{\text{mech}} - P\_{\text{elcc}} - D \Delta \omega}{M},
\tag{6}
$$

$$\dot{E}\_{fd} = -\frac{E\_{f-d}}{T\_A} + \frac{K\_A}{T\_A}(V\_{s0} - V\_s)\_\prime \tag{7}$$

$$\dot{E}\_{\rm sq}^{\dot{\prime}} = \frac{E\_{f-d} - \left(\mathbf{x}\_d - \mathbf{x}\_d^{\prime}\right)\dot{\mathbf{t}}\_d - E\_{\rm sq}^{\prime}}{T\_{d0}^{\prime}}.\tag{8}$$

*Pmech* and *Pelec* represent mechanical and electrical power, respectively, where *M* and *D* represent inertia and damping coefficient, respectively. ω<sup>0</sup> represents synchronous speed. δ and ω are rotor angle and rotor speed, respectively, while Δω = (ω − 1). *Ef d*, *E sq*, and *V* are field, internal, and the terminal voltage of generator. *T <sup>d</sup>*<sup>0</sup> represents the time constant of the open circuit. *TA* and *KA* are time constant and exciter gain, respectively, while *Vs*<sup>0</sup> represents reference voltage. *xd* is the *d-axis* reactance of generator, and *x <sup>d</sup>* is transient reactance. Furthermore, the different terms used in Equation (5) to Equation (8) are given as:

$$\begin{array}{c} P\_{\text{clxc}} = V\_{s-q}I\_{s-q} + V\_{s-d}I\_{s-d} \\ V\_s = \sqrt{V\_d^2 + V\_{q'}^2} \\ V\_d = X\_d I\_q \quad V\_q = E\_q' - X\_d' I\_d \\ I\_d = I\_{Ed} + I\_{Rd} \; I\_q = I\_{Eq} + I\_{Rq} \end{array} \tag{9}$$

Similarly, the line current for VSC-SE and VSC-SH can be obtained as [17–19]:

$$I\_{\rm sdd} = -\frac{X\_{\rm sh}}{X\_d \Sigma} E\_q' + \frac{X\_{\rm shd} m\_{\rm sh} \sin \delta\_{\rm sh} V\_{\rm dc}}{2X\_d \Sigma} - \frac{X\_{\rm td}}{X\_d \Sigma} \left(\frac{m\_{\rm sc} \sin \delta\_{\rm sc} V\_{\rm dc}}{2} + V\_b \cos \delta\right) \tag{10}$$

$$I\_{\rm seq} = -\frac{X\_{\rm shq}m\_{\rm sh}\sin\delta\_{\rm sh}V\_{\rm dc}}{2X\_{\rm q}\Sigma} + \frac{X\_{\rm solq}}{X\_{\rm q}\Sigma} \left(\frac{m\_{\rm ak}\sin\delta\_{\rm ak}V\_{\rm dc}}{2} + V\_{\rm b}\sin\delta\right) \tag{11}$$

$$I\_{\rm shd} = \frac{X\_{\rm sc}}{X\_d \Sigma} E\_q' - \frac{X\_{\rm scd} m\_{\rm shl} \sin \delta\_{\rm sh} V\_{\rm dc}}{2X\_d \Sigma} + \frac{X\_{\rm shd}}{X\_d \Sigma} \left(\frac{m\_{\rm shl} \sin \delta\_{\rm sh} V\_{\rm dc}}{2} + V\_b \cos \delta\right) \tag{12}$$

$$I\_{shq} = -\frac{X\_{shq}m\_{sh}\cos\delta\_{sh}V\_{dc}}{2X\_{q}\Sigma} + \frac{X\_{tq}}{X\_{q}\Sigma} \left(\frac{m\_{sc}\cos\delta\_{sc}V\_{dc}}{2} + V\_{b}\sin\delta\right) \tag{13}$$

where as *Xd*Σ, *Xtd*, *Xq*Σ, and *Xtq* are defined in [18].

#### **4. Supplementary Control Design**

Supplementary control schemes with the MRC are introduced for mitigating LFOs in the power system. The control design is used in the UPFC internal control structure to enhance the stability of the power system. As the faults occur on the transmission system, the speed of a different generator connected to the power system deviates and become out of synchronism. This speed deviation (Δω) is treated as an input to the supplementary controller, as depicted in Figure 3. The controllers mitigate the LFOs by providing suitable input to the UPFC, hence improving the power quality and stabilizing the power system. The design of the controllers is based on model reference neural control, and it uses a nonlinear autoregressive moving average model. The MRC is the integration of the (1) reference model, (2) NN plant, (3) NN controller, and (4) plant. The different parts of controller will be discussed in detail in a subsection of this section.

**Figure 3.** Control structure for oscillation damping. MRC = model reference controller; VSC = Voltage Sources Converters; PQ = Active Power, Reactive Power; PWM = Pulse Width Modulation; PLL = Phase Lock Loop.

#### *4.1. Neural Network (NN)*

The NN, inspired by the biological system, is composed of interconnected building blocks known as neurons. The interconnection weight quantifies the strength of interconnection. Because of these interconnected weights, the NN acquires the ability to learn and model complex relationships. Hence, it becomes more intelligent in predicting the uncertainties. Due to the learning capability of NN, it has a wide range of applications in automation and optimization [20], image processing [21,22], speech recognition [23,24], control [25,26], modelling, and time series prediction [27,28].

The architecture of NN is composed of multiple input layer and hidden layers, as depicted in Figure 4. This architecture is known as "*N-1*" feedforward NNs. Whereas '*N'* represents the number of hidden neurons and it has only one output node. The design of the architecture of NN depends on the input weights '*zij*' and output weights 'γ*i*'. The formulation of the NN through the nonlinear function of hidden layer activation function 'ψ*i*' is presented as:

$$
\hat{\varrho}(k+1) = \sum\_{i=1}^{N} \gamma\_i \psi\_i \big( \sum\_{j=1}^{n+m} z\_{ij} X\_j \big) \tag{14}
$$

where *X* = [*u*1, *u*2, ... ., *um*+*n*] is the input vector to NN and 'ψ*i*' is the activation function of the '*i th*' hidden layer. The neuron in the hidden layer uses the nonlinear logistic function given by Equation (15):

$$\psi\_i(\sum\_{j=1}^{n+m} z\_{ij} X\_j) = \frac{1}{1 + \exp\{\sum\_{j=1}^{n+m} z\_{ij} X\_j\}}.\tag{15}$$

The equation can be rearranged as:

$$
\hat{y}(k+1) = \boldsymbol{\theta}^T \boldsymbol{\psi},\tag{16}
$$

whereas ψ (ψ1,ψ2, ... .,ψ*N*) *<sup>T</sup>*, <sup>θ</sup>*<sup>T</sup>* (γ1, <sup>γ</sup>2, ... ., <sup>γ</sup>*N*). The gradient descent optimization technique is used to tune and update the output weight only. This will, in turn, reduce the computational complexity of the NN. The update equation of the NN is given as:

$$
\gamma\_{k+1} = \gamma\_k a \big( (y\_r \times \mathfrak{H}) \times (y\_r \times \mathfrak{H})' \big). \tag{17}
$$

**Figure 4.** Neural network (NN) architecture.

#### *4.2. Model Reference Controller*

The NN controller is used in the MRC for tracking the reference signal, as well as for plant identification. The MRC uses two NN controllers and minimal online computation. One NN is used to train the controller, and the other is used for plant identification. The data set for training the NN controller is given in Appendix A. The architecture of the MRC is depicted in Figure 5. The controller produces a controlled output signal *u*(*t*) that acts as input to the plant. This control signal brings the output signal *yp*(*t*) closer to the reference signal *yr*(*t*).

According to [29], the different inputs and their respective outputs define the characteristics of the real-time process. The previous information of input and output is used to introduce system dynamics. Equation (18) is sampled at a regular interval of time:

$$y\_p(t) = f\_p\left[y\_p(t-1), \dots, y\_p(t-n\_a), u(t-k), \dots u(t-k-n\_b)\right] \tag{18}$$

where *k* ≥ 1 shows the time delay in number of samples, *fp*[(*y*, *u*)] is processed by non-linearity, and *na*, *nb* + 1 represents the process output and input, respectively. Equation (18) gets into the shape of the transfer function below by linearizing the process about the operating point.

$$\frac{y\_p(t)}{u(t)} = q^{-k} \frac{B(q^{-1})}{A(q^{-1})} = q^{-k} \frac{b\_0 + b\_1 q^{-1} + \dots + b\_{n\_b} q^{-n\_b}}{1 + a\_1 q^{-1} + \dots + a\_{n\_b} q^{-n\_a}}.\tag{19}$$

**Figure 5.** NN-based MRC.

In Equation (19), *yp*(*t*) and *u*(*t*) represents the deviation of process from its operating points. The parameters *ai* and *bi* are dependent on the operating condition, while *q*−<sup>1</sup> represents the backward shift operator. The equation for the generalized linear controller is given in Equation (20)

$$u(t) = \frac{1}{F(q^{-1})} \Big( H(q^{-1})\mathbb{R}(t) - G(q^{-1})y\_p(t) \Big),\tag{20}$$

where '*R*' is set point, while *F*, *G*, and *H* are the polynomial defined below:

$$\begin{cases} F(q^{-1}) = 1 + f\_1 q^{-1} + \dots + f\_{n\_f} q^{-n\_f} \\ G(q^{-1}) = g\_0 + g\_1 q^{-1} + \dots + g\_{n\_g} q^{-n\_g} \\ H(q^{-1}) = h\_0 + h\_1 q^{-1} + \dots + h\_{nh} q^{-nh} \end{cases} \tag{21}$$

The closed-loop transfer function of a linearized system is given by:

$$\frac{y\_p(t)}{R(t)} = q^{-k} \frac{B(q^{-1})H(q^{-1})}{A(q^{-1})F(q^{-1}) + q^{-k}B(q^{-1})G(q^{-1})}.\tag{22}$$

Equation (20) and the control structure *nf* , *ng* and *nh* of Equation (21) is used to design a linear controller model to track the reference point. The polynomial B is neglected to make Equation (22) equal to the characteristic equation of the system. Hence, B is replaced with another term, i.e., *B* - *q*−<sup>1</sup> = *<sup>b</sup>*0*B*+- *q*−<sup>1</sup> . Equation (22) will transform into the following equation:

$$\frac{y\_p(t)}{R(t)} = q^{-k} \frac{b\_0 H(q^{-1})}{A(q^{-1})F(q^{-1}) + q^{-1}b\_0 G(q^{-1})}$$

$$\frac{y\_p(t)}{R(t)} = q^{-k\_m} \frac{s\_0 + s\_1 q^{-1} + \dots + s\_{n\_s} q^{-n\_s}}{1 + t\_1 q^{-1} + \dots + t\_{n\_l} q^{-nt}}.\tag{23}$$

*S q*−<sup>1</sup> and *T q*−<sup>1</sup> are the polynomial that is used to set closed-loop response to a desired value. *km* is the time delay of the reference model. In order to enable the cancellation of *<sup>B</sup>*+- *q*−<sup>1</sup> , *F* - *q*−<sup>1</sup> is replaced with *B*+- *q*−<sup>1</sup> *F*1 - *q*−<sup>1</sup> , i.e., *F* - *q*−<sup>1</sup> = *B*+- *q*−<sup>1</sup> *F*1 - *q*−<sup>1</sup> . The minimal order of the unique controller solution is: 

$$m\_{\mathcal{g}} = n\_a - 1, \ n\_f = n\_b + k - 1, \ k\_m = k, \text{ while } H(q^{-1}) = \frac{1}{b\_0} S(q^{-1}) \tag{24}$$

and

$$m\_t \le n\_a + k - 1.\tag{25}$$

Hence, the structure of process *na* and *k* is used to limit the order of denominator of closed-looped transfer function *nt*. The order of *nt* should fulfill the condition given in Equation (25). The lower order of denominator is possible by making unused *ti* coefficient equal to zero. The following Equation (26) is used to introduce an observer into the controller design.

$$\mathcal{C}\{q^{-1}\} = 1 + c\_1 q^{-1} + c\_2 q^{-2} + c\_{\text{fl}\_c} q^{-n\_c},\tag{26}$$

while

$$m\_{\mathcal{L}} \le n\_a + k - 1 - n\_t. \tag{27}$$

Then, Equation (25) is updated as:

$$\mathcal{H}\{q^{-1}\} = \frac{1}{b\_0} \mathcal{V}\{q^{-1}\} \mathcal{C}\{q^{-1}\}.\tag{28}$$

Although due to nonlinearity, the parameter of controller will change, but the controller designed structure in Equation (20) will not change for any operating point of the plant in Equation (18) despite setting the observer polynomial *C* - *q*−<sup>1</sup> = 1 and *ti* coefficient to zero. As in Equation (20), a nonlinear controller design with the same input and output can smoothly interpolate between different parameters of linear controllers as the operating points of the plant changes. The MRC in Equation (29) is proposed by considering Equation (20) in difference equation form and linearizing a nonlinear system into the form of Equation (20)

$$u(t) = f\_l[u(t-1), \dots u(t-n\_f), r\_f(t), y\_p(t), \dots, y\_p(t-n\_\mathcal{g})].\tag{29}$$

In Equation (29), the NN will be used to realize *fr*[■], and Equation (24) will define the order of *ng* and *nf* , while Equation (27) defines *nc* and *rf*(*t*) = *S* - *q*−<sup>1</sup> *C* - *q*−<sup>1</sup> *r*(*t*).

#### Controller Adaption with BP

The gradient descent algorithm with time delay is implemented to update the parameter of NN. The output of NN with a multilayer of a neuron is given as:

$$u(t) = f\_0(\sum\_j z\_j^2 \{ f\_h(\sum\_j z\_{ij}^1 x\_i(t) + b\_j^1) \} + b^2). \tag{30}$$

*xi* is the input to the NN, and *b* and *z* are biases and weights of the NN, respectively. *f*<sup>0</sup> represents the nonlinear function of the input layer, and *fh* is function of the hidden layer. The cost function that needs to be minimized is given as:

$$J(t) = \frac{1}{2}e(tj)^2 = \frac{1}{2}(y\_r(t) - y\_p(t))^2. \tag{31}$$

*yr* and *yp* are the reference and measured output of the system, respectively. The weight parameter '*z*' is used in opposite direction in cost function gradient to minimize the error, given as:

$$\mathbf{z}(t) = \mathbf{z}(t-1) - \eta \frac{\partial I(t)}{\partial \mathbf{z}}.\tag{32}$$

The unknown 'η' represents the learning rate. The BP algorithm is used to calculate the gradient algorithm of Equation (32) Hence, Equation (18) and <sup>∂</sup>*J*(*t*) <sup>∂</sup>*yp*(*t*) <sup>=</sup> <sup>−</sup>*e*(*t*) gives:

$$\frac{\partial f(t)}{\partial z} = -\epsilon(t) \frac{\partial y\_p(t)}{\partial u(t-k)} \frac{\partial u(t-k)}{\partial z}.\tag{33}$$

In Equation (33), the first derivative is called process Jacobian. The second NN is used to estimate it as a process model, i.e., the NN model given by:

$$y^m(t) = f\_m[\bullet] = f\_0\{\sum\_j z\_j^{2m} \{ f\_h(\sum\_i z\_{ij}^{1m} \mathbf{x}\_i^m(t) + b\_j^{1m}) \} + b^{2m} \}.\tag{34}$$

The delayed model inputs/outputs are represented by *x<sup>m</sup> <sup>i</sup>* (*t*), as shown in Equation (18) in the form of *u*(*t* − *k*). Equation (35) shows one of the estimations of Jacobian.

$$\frac{\partial y\_p(t)}{\partial u(t-k)} \approx \frac{\partial y^m(t)}{\partial u(t-k)} = f\_0' \sum\_j z\_j^{2m} f\_{h\_j}' z\_{1j}^{1m} \,. \tag{35}$$

The process sign is used to approximate Jacobian if the model is not available. Equation (36) shows the update rule for the different parameter of the controller.

$$z(t) = z(t-1) + \eta e(t) \text{sgn} \left(\frac{\partial y\_p(t)}{\partial u(t-k)}\right) \frac{\partial u(t-k)}{\partial z}.\tag{36}$$

This equation is used to update the rule to get better dynamics of the system.

#### **5. Results and Discussion**

The performance of the designed controller is evaluated and compared with the PI controller under symmetrical and asymmetrical faults. The faults are applied for 50 ms, i.e., the fault is initiated is at *t* = 200 ms and cleared at *t* = 250 ms. Moreover, the fault is applied on the mid-point of the lower line of the double circuit line between buses *B*<sup>1</sup> and *B*2, as depicted in Figure 1. The proposed controller has the capability to mitigate the LFOs and enhance the stability of the power system. The parameters used for comparison of the proposed control structure are: (1) load bus, voltage, (2) speed deviation (Δ*w*), and (3) performance indices of controllers.

#### *5.1. Asymmetrical Faults*

Both asymmetrical and symmetrical faults are applied to the system for 50 ms (at *t* = 200 ms) to create oscillations in the system. The fault creates a voltage dip of 10% during the 1-Φ fault, while a 25% dip is created during the 2-Φ to ground fault. After the clearance of the 1-Φ fault, the voltage profile is restored at *t* = 2 s due to the adaptive and predictive nature of MRC, while the response with conventional PI is oscillatory and oscillation is damped out after *t* = 12 s, as shown in Figures 6 and 7, respectively. Similarly, the speed deviation shows high amplitude oscillation in the case of the PI controller, while the MRC has the capability to damp out the oscillation in comparatively less time. Figure 8 shows that the MRC has damped out the oscillation at *t* = 2 s, and speed deviation (Δ*w* ≈ 0) is approximately zero, while in comparison with MRC, the PI controller deviates from zero till *t* = 6 s, as depicted in Figure 8.

**Figure 8.** Speed deviation during 1-Φ to ground fault.

Similarly, during the 2-Φ to ground fault, the MRC has a better response in mitigation of oscillation as compared to the PI. As depicted in Figure 9, the MRC has smoothed the voltage response abruptly after fault clearance at *t* = 2.5 s, while the PI controller failed to show the robustness. The oscillation is completely damped out at *t* = 12 s, as shown in Figure 7. Similarly, speed deviation becomes zero at *t* = 2 s in the case of the MRC, while the PI vanishes the oscillations at *t* = 8 s.

**Figure 9.** Load bus voltage 2-Φ to ground fault.

#### *5.2. Symmetrical Faults*

In the case of an asymmetrical fault, the 3-Φ to ground fault is applied to evaluate the performance of the proposed controller, as depicted in Figure 10. The MRC has a better approach toward the mitigation of oscillations occurring after the voltage dip, as shown in Figure 11. The oscillation is completely damped out at *t* = 2.5 s in the case of the MRC, while the PI has low amplitude oscillations till *t* = 8 s. Similarly, speed deviation is subjected to larger amplitude high frequency oscillation and higher settling time using the PI, while the proposed controller settles the response to zero at *t* = 3 s.

**Figure 10.** Load bus voltage 3-Φ to ground fault.

**Figure 11.** Speed deviation during 3-Φ to ground fault.

Besides other performance evaluation parameter, Table 1 shows the performance indices of both controllers during all three types of faults. The lower value of integral square error (ISE) and integral absolute error (IAE) shows better response of controller. The MRC controller has both a lower value of ISE error and IAE.

**Table 1.** Performance indices. ISE = integral square error; IAE = integral absolute error.


#### **6. Conclusions**

The mitigation of power oscillations is a demanding task to enhance the stability of the power system. An intelligent, optimized, and robust design of the UPFC plays a vital role in damping of these oscillations. With these abilities, the MRC controller has shown tremendous performance in mitigating power oscillations, enhancing system stability, and improving the capability of the control system. Moreover, our control scheme shows robustness and intelligence during system perturbation. In order to show the effectiveness of proposed control scheme, it was compared with a traditional PI controller. The results clearly show the superiority of the NN-based MRC over the traditional PI controller, as it is better in oscillation damping. The NN-based MRC diminishes the amplitude of power oscillation, produced due to fault.

However, the MRC-based on a NN is complex algorithm. Due to its complexity, it requires considerable amount of time for calculation and implementation. Due to the integration of the NN controller in the MRC frame, it is associated with several issues and complications. The complications include selecting the input and output of the controller, as well as training of the NN for its weight adjustment.

Moreover, soon the proposed control will be extended to more than a 2-area system with a greater number of generation stations. Additionally, a new control scheme will be designed, and its results will be compared with both the PI and the MRC.

**Author Contributions:** W.U., N.Z., and K.Z. propose the main idea of the paper. W.U., M.I., S.U.I., A.T. and A.H have implemented the mathematical derivations, simulation verifications, and analyses. The paper is written by W.U. and is revised by N.Z., I.K., A.H., and H.-J.K. All the authors were involved in preparing the final version of this manuscript. Besides, this whole work is supervised and funded by H.-J.K. and G.-S.P.

**Funding:** This research was supported by Basic Research Laboratory through the National Research Foundations of Korea funded by the Ministry of Science, ICT and Future Planning (NRF-2015R1A4A1041584).

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Appendix A**


**Table A1.** Parametric Value of NN.

#### **References**


© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

*Article*
