Active-Current Control of Large-Scale Wind Turbines for Power System Transient Stability Improvement Based on Perturbation Estimation Approach

This paper proposes an active-current control strategy for large-scale wind turbines (WTs) to improve the transient stability of power systems based on a perturbation estimation (PE) approach. The main idea of this control strategy is to mitigate the generator imbalance of mechanical and electrical powers by controlling the active-current of WTs. The effective mutual couplings of synchronous generators and WTs are identified using a Kron-reduction technique first. Then, the control object of each WT is assigned based on the identified mutual couplings. Finally, an individual controller is developed for each WT using a PE approach. In the control algorithm, a perturbation state (PS) is introduced for each WT to represent the comprehensive effect of the nonlinearities and parameter variations of the power system, and then it is estimated by a designed perturbation observer. The estimated PS is employed to compensate the actual perturbation, and to finally achieve the adaptive control design without requiring an accurate system model. The effectiveness of the proposed control approach on improving the system transient stability is validated in the modified IEEE 39-bus system.


Introduction
Transient stability, which reflects the ability of power systems in maintaining synchronism when subjected to a severe fault, is an important factor in the design and operation of power systems [1].A transient stability control strategy (TSCS) is essential to prevent a system from losing synchronism.The traditional TSCS includes generator tripping [2], load shedding [3,4], and excitation boosting control [5], etc.However, the improvement provided by the excitation boosting control is limited due to the restriction of maximum exciting current.The generator tripping and load shedding would bring about undesired economic loss.Considerable efforts are still devoted to finding effective and economical TSCSs to improve the system transient stability.
The penetration level of wind power in systems has increased significantly around the world during the past decades [6,7].Among the new installed wind turbines (WTs), the doubly fed induction generator (DFIG) is the most popular technology for high power conversion efficiency and low manufacturing cost [8,9].The DFIGs are connected to the system though power electronic converters and could provide flexible controllability.Although the integration of large-scale WTs could either benefit or deteriorate the system transient stability, depending on the actual situation, transient stability can be improved significantly by applying a supplementary control strategy on WTs [10,11].Reference [11] proposed a supplementary control strategy to modulate the active-power of WTs based on the terminal bus frequency for a predefined interval subsequent to a large disturbance.A current reference control strategy was proposed in Reference [12] to temporarily provide reactive current support during a terminal voltage dip, and this consequently improved the system transient stability.Reference [13] proposed a synthetic control to provide reactive-current support and to reduce the active-current for transient stability improvement.However, most existing TSCSs are designed for single WT.For a power system integrated with multiple large-scale WTs, how to coordinate these WTs to improve the transient stability has not been investigated.
The nonlinearities inherited in power systems play an important role in TSCS design.To cope with these nonlinearities, many nonlinear control approaches, such as controlled Lyapunov approach [14,15] and feedback linearization approach [16,17], have been implemented.However, those approaches require an accurate model of power systems, which is difficult to obtain in reality.A perturbation estimation (PE) approach provides an alternative to solve this problem [18].The main idea of PE is to treat the nonlinearities and external disturbances as a perturbation state (PS), and to then try to estimate it by a designed perturbation observer (PO) in real time.With the compensation of the estimated PS, the control design does not require accurate models of power systems, and it consequently provides good robustness to external disturbances.
In this study, a perturbation estimation-based control (PEC) is proposed to control the active-current of WTs for improving power system transient stability.The mutual couplings of generators and WTs are identified using a Kron-reduction approach, and then the control object of each WT is assigned based on the identified mutual couplings.For each WT, an individual controller is developed based on a PE approach to control its active-current for transient stability improvement without requiring accurate models of power systems.Finally, the developed control is tested in the modified IEEE 39-bus system, and the simulation results verify its effectiveness for improving the power system transient stability as well as its robustness to wind speed variation and communication delay in wide area measurement system (WAMS).
The rest of this paper is organized as follows.Section 2 analyzes the mutual couplings between generators and WTs.The design of PEC is presented in Section 3, and Section 4 gives the simulation results of the proposed control approach.Finally, Section 5 concludes the paper.

Mutual Couplings between Generators and WTs
A power system, consisting of n generators and m WTs, is considered to analyze the effective mutual couplings of generators and WTs as depicted in Figure 1.
For analysis purposes, several assumptions are made as follows: • System loads are modeled as constant impedance loads and can therefore be absorbed into the node admittance matrix.

•
A subtransient model, with subtransient potential in series with subtransient impedance, is considered for synchronous generators.
It should be noted that these assumptions are only for mutual coupling analysis and control strategy design.The designed control strategy is tested in the detailed model of power systems, where the detailed generator model is equipped with an automatic voltage regular (AVR) and governor.The composite load model is also considered.It should be noted that these assumptions are only for mutual coupling analysis and control strategy design.The designed control strategy is tested in the detailed model of power systems, where the detailed generator model is equipped with an automatic voltage regular (AVR) and governor.The composite load model is also considered.With the generator subtransient impedances and load impedances included in the node admittance matrix, we obtain where E G = [E G1 , . . ., E Gn ] T and I G = [I G1 , . . ., I Gn ] T are the generator potential and current vectors, respectively.Y GG , Y Gw , Y wG , and Y ww are the partitioned matrices of the reduced system admittance matrix with all nodes other than the generator internal nodes and WT terminal nodes eliminated, respectively.V w = [V w1 , . . ., V wm ] T is the voltage vector of WT terminal nodes.I w = [I w1 , . . ., I wm ] T is the WT injection current vector.Kron-reduction is used to determine the effective equivalent admittance of generators and WTs by eliminating the voltage vector V w , and we obtain The ith row, jth column element of Y GG , Y GG (i,j), is the Kron-reduced equivalent admittance between generator i and generator j, determining the effective mutual coupling of them.The ith row, jth column elements of Y Gw , Y Gw (i,j), is the Kron-reduced equivalent admittance between generator i and WT j, determining the effective mutual coupling of them.
From (2), the generator current I G and, further, the generator electrical power P G = Re(E G ,I G ) can be impacted by controlling WT injection I w .Certainly, these impacts are different due to the different mutual couplings between generators and WTs, which can be quantified by the associated Kron-reduced equivalent admittance.The larger the equivalent admittance is, the stronger the mutual coupling is.

Control Object Assignation
Following a severe fault, instability occurs when one or some generators separate from the rest irrecoverably due to the imbalance of their mechanical and electrical power.Mitigating the generator imbalance power by controlling WT injection currents will bring about a transient stability improvement.However, the effects of controlling WT injection currents on mitigating the generator imbalance power are different due to the different effective mutual couplings of generators and WTs, which can be identified by Y Gw .The stronger the effective mutual couplings are, the more prominent the effects are.
Based on the analysis in Section 2, a control object assignation approach is proposed from the control effect point of view.System generators are clustered into non-lapped generator groups based on the mutual couplings of generators and WTs, and each cluster is controlled by a WT.The control object assignation approach consists of 6 steps, as follows: 1.
Establish the node admittance matrix of power systems Y sys .2.
Employ the Kron-reduction approach on Y sys to obtain Y' Gw .

3.
Transform Y' Gw to Y mn by setting the elements of Y' Gw , whose amplitudes are less than a specified value κ, as zeros to neglect the weak mutual couplings between generators and WTs. 4.
Denote the ith row of Y mn as Y i , representing the effective mutual couplings of generator i and each WT.

5.
Generator i is clustered into generator group j (Denoted as GG j ), where j is calculated by Assign GG j as the control object of WT j.
It should be noted that some generators may not be clustered into any group due to their weak mutual couplings with each WT.These generators are not assigned as the control object of any WT, because they can hardly be controlled by any WT.

Controlled System Model
The models of WT and generator used for control design are described in detail in this subsection.The detailed model of DFIG is given in Figure 2, which consists of a wind turbine, a transmission system, a wound rotor induction generator, and a back-to-back converter.In addition, a crowbar system is used to protect the converter from over-current during the fault, and consequently, this achieves the continuous operation of the DFIG during the transient.It should be noted that some generators may not be clustered into any group due to their weak mutual couplings with each WT.These generators are not assigned as the control object of any WT, because they can hardly be controlled by any WT.

Controlled System Model
The models of WT and generator used for control design are described in detail in this subsection.The detailed model of DFIG is given in Figure 2, which consists of a wind turbine, a transmission system, a wound rotor induction generator, and a back-to-back converter.In addition, a crowbar system is used to protect the converter from over-current during the fault, and consequently, this achieves the continuous operation of the DFIG during the transient.The mechanical power Pt captured by a wind turbine from the wind can be expressed as where ρ, r, and A are the air density, rotor radio, and rotor-swept area of a wind turbine, respectively; vwind, Cp, λ, ωt, and β are the wind speed, power coefficient, tip-speed ratio, turbine speed, and pitch angle, respectively.A pitch controller is employed to control the pitch angle β for preventing a turbine from over-speed, as shown in Figure 3.The mechanical power P t captured by a wind turbine from the wind can be expressed as where ρ, r, and A are the air density, rotor radio, and rotor-swept area of a wind turbine, respectively; v wind , C p , λ, ω t , and β are the wind speed, power coefficient, tip-speed ratio, turbine speed, and pitch angle, respectively.A pitch controller is employed to control the pitch angle β for preventing a turbine from over-speed, as shown in Figure 3.
where ρ, r, and A are the air density, rotor radio, and rotor-swept area of a wind turbine, respectively; vwind, Cp, λ, ωt, and β are the wind speed, power coefficient, tip-speed ratio, turbine speed, and pitch angle, respectively.A pitch controller is employed to control the pitch angle β for preventing a turbine from over-speed, as shown in Figure 3.A transmission system, consisting of a low-speed drive shaft, gearbox, and a high-speed drive shaft, transmits the mechanical power Pm to the wound induction generator, which can be modeled as a two-mass shaft model in (5).A transmission system, consisting of a low-speed drive shaft, gearbox, and a high-speed drive shaft, transmits the mechanical power P m to the wound induction generator, which can be modeled as a two-mass shaft model in (5).
where H t , D t and P t are the inertia time constant, damping constant, and mechanical torque of a wind turbine, respectively; H g , D g , ω g , and P g are the inertia time constant, damping constant, speed, and mechanical torque of a DFIG, respectively; K s , D s , θ s , and ω s are the shaft stiffness, damping constant, torsional twist, and synchronous speed, respectively.
The wound induction generator model in the d-q reference frame can be expressed as follows: where ω slip = ω s − ω g is the slip of DFIG.V sd , V sq and V rd , V rq represent the d-axis and q-axis components of the stator voltage and rotor voltage, respectively.i sd , i sq and i rd , i rq represent the d-axis and q-axis components of the stator current and rotor current, respectively.ψ sd , ψ sq and ψ rd , ψ rq represent the d-axis and q-axis components of the stator flux and rotor flux, respectively.L s , L r , and L m represent the stator inductor, rotor inductor, and the mutual inductor of stator and rotor, respectively.R s and R r represent the stator resistance and rotor resistance, respectively.A back-to-back converter consists of a gird side converter (GSC) and a rotor side converter (RSC), where the GSC maintains the DC voltage of the capacitor and regulates the reactive power of the GSC exchanged with the grid, and the RSC performs the maximum power point tracking (MPPT) control and regulates the stator reactive power injected into the grid.The standard vector control technique is applied on both the RSC and GSC, and the d-axis of the synchronously rotating frame is oriented along the stator voltage V w as shown in Formula (7), where the stator resistance R s is neglected for being much smaller than its reactance L s .
Therefore, the injection current of DFIG j can be expressed Energies 2018, 11, 1995 where the d-axis stator current i sd,j and q-axis stator current i sq,j can be controlled by the rotor active-current i rd,j and rotor reactive-current i rq,j respectively as shown in (9).
i sd,j = L m,j i rd,j /L s,j i sq,j = − V wj + L s,j L m,j i rq,j /L s,j ω s (9) where L s,j and L m,j represent the stator inductor and mutual induction of stator and rotor of DFIG j, represectively.V wj is the amplitude of the terminal voltage V wj .
A typical converter controller is shown in Figure 4, which consists of a fast time-scale current controller, a medium time-scale DC voltage controller, and a slow time-scale speed controller.The proposed PEC is integrated in the RSC, and the auxiliary signal u j , derived from the PEC, is added on the active-current reference I rdref,j to modulate the active-current of DFIG for improving power system transient stability.In terms of transient stability, the time-scale of interest (approximate 0.1-1 s) is much slower than the time-scale of the current controller (approximate 0.01-0.1 s), as shown in Figure 4 and demonstrated in [19].Therefore, the dynamics of the current controller can be omitted, and the d-axis current Ird can be considered as instantaneously tracking their reference values, as follows: The swing Formula models, which link the rotor acceleration to the imbalance between the mechanical power supplied by the governor and the electrical power, are considered to describe the dynamics of generators in power systems, as follows: where Mi is the inertial constant.δi and ωi are the rotor angle and rotor speed of generator i, respectively.Pmi and PGi are the mechanical and electrical power of generator i, respectively.
PGi can be expressed as In terms of transient stability, the time-scale of interest (approximate 0.1-1 s) is much slower than the time-scale of the current controller (approximate 0.01-0.1 s), as shown in Figure 4 and demonstrated in [19].Therefore, the dynamics of the current controller can be omitted, and the d-axis current I rd can be considered as instantaneously tracking their reference values, as follows: The swing Formula models, which link the rotor acceleration to the imbalance between the mechanical power supplied by the governor and the electrical power, are considered to describe the dynamics of generators in power systems, as follows: where M i is the inertial constant.δ i and ω i are the rotor angle and rotor speed of generator i, respectively.P mi and P Gi are the mechanical and electrical power of generator i, respectively.
P Gi can be expressed as Substituting Formula ( 8) and ( 9) into (12), we obtain where E Gi and δ Gi are the amplitude and phase of E Gi , respectively.V wj and θ wj are the amplitude and phase of V wj , respectively.Y Gw (i,j) and α ij are the amplitude and phase of Y Gw (i,j), respectively.From Formula (13), P Gi consists of two terms: P GiG and P Giw , which represent the exchanged power with the other generators and WTs, respectively.It is clear that the active-power of generator i can be regulated by controlling either the injection active-current or reactive-current of WTs.Considering the fact that the transient stability is high relative to the active-power, only the active-current of WTs is considered in PEC design.
From ( 11)-( 13), the model consisting of n j generators within GG j and WT j can be expressed as I drre f ,n j + u n j − c Qn j I sq,n j y j = h j (δ 1 , ω 1 , . . ., δ n j , ω n j ) where n j is the total number of generators in GG j .y j = h j (δ 1 , ω 1 , . . .,δ nj , ω nj ) is the output variable as illustrated in Formula (15).

Design of PEC
The transient stability control object is to restore the generator rotor speeds back to the synchronous speed ω s , consequently making generators synchronize again.Therefore, in this paper, the output variable y j is selected as the average value of rotor speeds of generators within GG j , as follows: To employ a PE approach, the linearization model of a power system should be established by calculating the derivative of the variable y j until the control input signal u wj appears For Formula ( 16), a PS ψ j is defined in Formula ( 17) to represent the comprehensive perturbation of the nonlinearities and parameter variation.
where b 0j is a positive constant.Substituting Formula ( 16) into (15), we obtain the following: where λ is the derivative of ψ j .Instead of calculating ψ j according to its definition in Formula ( 18), a PO is designed for Formula (17) to estimate ψ j in real time.
From Formula ( 18) and ( 19), the error dynamics of a PO can be obtained by The observer coefficients should be set properly so as to place the pole of Equation ( 20) at the desired location in the open left-half complex plane.By choosing a positive constant as the pole of Equation ( 21), the observer coefficients can be easily obtained.
where λ j is the chosen pole.
To achieve quick and accurate estimation of PS, the pole λ j of ( 21) should be set large enough.However, an overlarge λ j will increase the sensitivity to the noise in signaling.Therefore, the selection of λ j must make a trade-off between the tracking performance and the noise tolerance.
With the estimated PS, the final control law is given in ( 22), where the estimated PS ψ j is used to reject the comprehensive disturbance.
where k j is a positive constant.The larger the value of k j is, the better performance PEC provides.However, an overlarge k j will exploit more kinetic energy stored in the rotating mass, which aggravates the fluctuation of rotor speed and may violate the specified speed range of [0.7, 1.25] p.u.Therefore, the selection of k j should make a trade-off between the control performance and the rotor speed fluctuation.
To prevent a PEC from activating in cases of small rotor speed deviation due to the small disturbance occurred in the power system, a dead-band block is designed as follows: where ε j is a small positive constant, which is set by the power system operator.
Energies 2018, 11, 1995 9 of 15 Moreover, before adding the output of PEC on the active-current reference, the control output of a PEC should pass through a limiter to prevent overmodulation.A PEC exploits the available kinetic energy of a WT to improve system transient stability.However, the available kinetic energy depends on the turbine speed, which varies with the wind speed under a MPPT control mode.The overmodulation of a PEC may cause the turbine speed to violate the permissible range, consequently causing WT tripping.In this paper, the output control of PEC u j is limited within the range of [−0.3, 0.3], based on the simulation results.How to select appropriate configurations of limits under different wind speed conditions requires further research, which is not the focus of this paper.
To clearly illustrate the structure of the proposed PEC, a block diagram is given in Figure 5.
where εj is a small positive constant, which is set by the power system operator.Moreover, before adding the output of PEC on the active-current reference, the control output of a PEC should pass through a limiter to prevent overmodulation.A PEC exploits the available kinetic energy of a WT to improve system transient stability.However, the available kinetic energy depends on the turbine speed, which varies with the wind speed under a MPPT control mode.The overmodulation of a PEC may cause the turbine speed to violate the permissible range, consequently causing WT tripping.In this paper, the output control of PEC uj is limited within the range of [−0.3, 0.3], based on the simulation results.How to select appropriate configurations of limits under different wind speed conditions requires further research, which is not the focus of this paper.
To clearly illustrate the structure of the proposed PEC, a block diagram is given in Figure 5.

State error feedback control
Perturbation observer

Test System
To test the performance of the developed PEC, a modified IEEE 39-bus system integrated with four DFIG-based wind farms (WFs) of rated capacity 600 MW was established on a DIgSILENT platform (version 2018 SP1, DIgSILENT/PowerFactory Company, Berlin, Germany), as depicted in Figure 6.The four WFs operated at different wind speed conditions, as shown in Table 1, and the total wind power accounted for 18.9% of the overall power generation.In practice, each WF contained a large number of small capacity WTs.Establishing the detailed model of every WF would result in a heavy simulation burden and increase the simulation time.Therefore, the aggregation technique, which is commonly used in the published papers [11], was taken by modeling a WF as a single WT of equivalent capacities.The configurations of the WT were taken from Formula [20].The subtransient models were considered for generators.All generators other than G 1 , which is an equivalent external generator, were equipped with AVR type IEEE-ET1 and speed governor type IEEE-G1.The composite load models consisting of 50% constant load and 50% induction motor load were considered.The subtransient models were considered for generators.All generators other than G1, which is an equivalent external generator, were equipped with AVR type IEEE-ET1 and speed governor type IEEE-G1.The composite load models consisting of 50% constant load and 50% induction motor load were considered.
The modified IEEE 39-bus system.

Configuration of PEC
The element amplitudes of Y ′ Gw were calculated as shown in Formula ( 24) to analyze the mutual couplings of generators and WTs.It was clear, for each WT, that there only existed a small number of strong coupled generators.By setting κ = 10 −2 , generators were clustered into four generator groups {GG1}, {GG2}, {GG3}, and {GG4} consisting of {G9}, {G8, G10}, {G2, G3}, and {G4, G5, G6, G7}, respectively, and these groups were assigned as the control object of each WT.
After control object assignation, an individual controller was developed for each WT as described in Section 3.1.There were only 3 parameters that required configuration for the PEC: the observer pole λj and the positive constants b0j and kj.These parameters were selected by trial and

Configuration of PEC
The element amplitudes of Y Gw were calculated as shown in Formula (24) to analyze the mutual couplings of generators and WTs.It was clear, for each WT, that there only existed a small number of strong coupled generators.By setting κ = 10 −2 , generators were clustered into four generator groups {GG 1 }, {GG 2 }, {GG 3 }, and {GG 4 } consisting of {G 9 }, {G 8 , G 10 }, {G 2 , G 3 }, and {G 4 , G 5 , G 6 , G 7 }, respectively, and these groups were assigned as the control object of each WT.
After control object assignation, an individual controller was developed for each WT as described in Section 3.1.There were only 3 parameters that required configuration for the PEC: the observer pole λ j and the positive constants b 0j and k j .These parameters were selected by trial and error, and all controllers took the same configuration: λ j = 10, b 0j = 1, k j = 10.Therefore, the observer coefficients l j1 and l j2 were set as 20 and 100, respectively. (24)

General Performance of PEC
The performance of the PEC was evaluated by applying a three-phase-to-ground fault at the middle of the transmission line between buses 15 and 16 at 1 s and then isolating the fault line at 1.12 s.The rotor angles and rotor speeds of all generators without and with the PEC were compared, as shown in Figure 7, which clearly illustrates the effectiveness of the proposed PEC in improving the transient stability of power systems.

General Performance of PEC
The performance of the PEC was evaluated by applying a three-phase-to-ground fault at the middle of the transmission line between buses 15 and 16 at 1 s and then isolating the fault line at 1.12 s.The rotor angles and rotor speeds of all generators without and with the PEC were compared, as shown in Figure 7, which clearly illustrates the effectiveness of the proposed PEC in improving the transient stability of power systems.The dynamics of all WTs without and with the PEC are compared in Figure 8.The active-powers of WTs were regulated by controlling the active-current of WTs with the PEC to improve the system transient stability, as shown in Figure 8a,b.It should be noted that the crowbar of DFIG W4 was activated at 1.01 s to protect converter from over-current, and then it was removed at 1.12 s after the fault clearance, consequently achieving the continuous operation of the WT.The developed PEC was activated at 1.13 s after the removal of the crowbar to control the active-current of the WT.The PEC exploited the kinetic energy stored in the rotating mass by changing the turbine speed.However, due to the large inertia of the rotating mass, a small speed fluctuation occured, which did not threaten the stable operation of the WT (shown in Figure 8d).The dynamics of terminal voltage are compared in The dynamics of all WTs without and with the PEC are compared in Figure 8.The active-powers of WTs were regulated by controlling the active-current of WTs with the PEC to improve the system transient stability, as shown in Figure 8a,b.It should be noted that the crowbar of DFIG W 4 was activated at 1.01 s to protect converter from over-current, and then it was removed at 1.12 s after the fault clearance, consequently achieving the continuous operation of the WT.The developed PEC was activated at 1.13 s after the removal of the crowbar to control the active-current of the WT.The PEC exploited the kinetic energy stored in the rotating mass by changing the turbine speed.However, due to the large inertia of the rotating mass, a small speed fluctuation occured, which did not threaten the stable operation of the WT (shown in Figure 8d).The dynamics of terminal voltage are compared in Figure 8e,f, which illustrates that the developed PEC does not violate the specified voltage limits ([0.7-1.2]p.u.) during or following a severe fault.The critical clearing time (CCT) is a widely used metric in evaluating the system transient stability.The CCTs with and without the PEC under different fault conditions are compared in Table 2.It can be seen that the transient stability margin with the PEC increases effectively compared to that without the PEC, which verifies the effectiveness of the proposed PEC.Since the developed PEC takes wide-area measurements as its control input signal, the communication delay in the WAMS had an adverse effect on the control performance of the PEC.In reality, the communication delay in the WAMS using optical fiber was usually within the range of 80-150 ms [21].The dynamics of generator rotor angle and rotor speed with 100 ms and 200 ms communication delays are shown in Figure 12, respectively, which are similar to the dynamics with no communication delay, as shown in Figure 7.

Robustness to the Communication Delay
Since the developed PEC takes wide-area measurements as its control input signal, the communication delay in the WAMS had an adverse effect on the control performance of the PEC.In reality, the communication delay in the WAMS using optical fiber was usually within the range of 80-150 ms [21].The dynamics of generator rotor angle and rotor speed with 100 ms and 200 ms communication delays are shown in Figure 12, respectively, which are similar to the dynamics with no communication delay, as shown in Figure 7.The CCTs with different communication delays are given in Table 2.The increase of communication delay has little impact on the performance of the PEC.Therefore, the proposed PEC is robust to the possible and extreme delays in the communication network.

Conclusions
This paper proposes an active-current control for large-scale WTs to improve power system transient stability.The control object of each WT is assigned based on the identified effective mutual couplings of generators and WTs using a Kron-reduction approach, and an individual controller is designed for each WT using a PE technique to control its active-current for improving the system transient stability.The design of the PEC does not require an accurate model of power systems.The simulation results, carried in the modified IEEE 39-bus system, verify the effectiveness of the proposed PEC on the system transient stability improvement under different fault conditions as well as its robustness to wind speed variation and communication delay in WAMS.The CCTs with different communication delays are given in Table 2.The increase of communication delay has little impact on the performance of the PEC.Therefore, the proposed PEC is robust to the possible and extreme delays in the communication network.

Figure 5 .
Figure 5.The block diagram of PEC.

Figure 5 .
Figure 5.The block diagram of PEC.

Figure 7 .
Figure 7. Dynamics of the generator (a) rotor angle and (c) rotor speed without PEC; Dynamics of the generator (b) rotor angle and (d) rotor speed with PEC.

Figure 7 .
Figure 7. Dynamics of the generator (a) rotor angle and (c) rotor speed without PEC; Dynamics of the generator (b) rotor angle and (d) rotor speed with PEC.

Energies 2018 ,Figure 8 .
Figure 8e,f, which illustrates that the developed PEC does not violate the specified voltage limits ([0.7-1.2]p.u.) during or following a severe fault.

Figure 10 .Figure 10 .Figure 11 .
Figure 10.Dynamics of the generator (a) rotor angle and (b) rotor speed without PEC; Dynamics of the generator (c) rotor angle and (d) rotor speed with PEC.

Figure 12 .
Figure 12.Dynamics of generator (a) rotor angle; (c) rotor speed with 100 ms delay; dynamics of generator (b) rotor angle; (d) rotor speed with 200 ms delay.

Figure 12 .
Figure 12.Dynamics of generator (a) rotor angle; (c) rotor speed with 100 ms delay; dynamics of generator (b) rotor angle; (d) rotor speed with 200 ms delay.

Table 1 .
Wind speeds and active-power of different WFs.