Stability Domain Analysis and Enhancement of Squirrel Cage Induction Generator Wind Turbines in Weak Grids

: There are signiﬁcant concerns regarding the stability of increased wind power generation in weak power grids. This paper investigates and improves the stability of Wind Turbine Squirrel Cage Induction Generators (WT-SCIGs) with series compensation and weak interconnections to the power grid. Detailed time-domain and state-space modeling have revealed new bifurcations and oscillatory modes for a WT-SCIG connected radially to a weak grid through a series compensated line. The stability domain analyses are carried out by computing bifurcations in the system by analyzing eigenvalues of the linearized system. The analyses demonstrate for the ﬁrst time how the degree of compensation at which the Hopf bifurcation occurs depends on the X / R ratio of the line, operating slip of the induction generator, and voltage regulator parameters as well as the time delays in measurements. A new damping controller is proposed, which greatly improves the dynamic stability of the WT-SCIG and eliminates destructive Hopf bifurcations in weak grids for a wide range of series compensation. This allows for a much larger percentage of series compensation than what is usually possible, while avoiding instabilities, thereby maximizing the power transfer capability.


Introduction
Utilizing series capacitors in power systems is one of the most economical ways to improve the power flow through existing transmission lines [1]. Unfortunately, the level of series compensation is often limited in many systems by the existence of Hopf bifurcations giving rise to subsynchronous resonance (SSR), which has been known to cause severe instabilities in transmission systems [2][3][4].
Subsynchronous interactions in Doubly Fed Induction Generator (DFIG)-based wind farms are well documented. However, only a few articles [5][6][7][8] have been published on SSR in Wind Turbine-Squirrel Cage Induction Generators (WT-SCIGs). WT-SCIGs operate at a fixed speed mostly in grid-connected mode. Their main advantages are simplicity in construction and maintenance [9]. The generator requires reactive power for its operation, which results in a poor power factor and is sometimes corrected using shunt capacitor banks at the generator terminals [10] or a Static Var Compensator (SVC). In this paper, it is shown that an SVC can be modified with an improved controller to eliminate all critical Hopf bifurcations in a weak grid.
Since renewable energy sources, such as wind power, may be located in remote areas, their connections to the grid are usually through long feeders of high impedance, making the interconnection weak. This has raised significant concerns about the stability of renewable power generations in weak power grids. In reference [11], the stability of voltage source inverters is investigated in a weak grid characterized by a low Short Circuit Ratio (SCR), low inertia and high impedance. In [5], a detailed eigenvalue analysis is presented for a Type-1 wind turbine to study the influence of series compensation and 1.
Comprehensive analysis of bifurcations in weak grids and uncovering new modes, which have not been reported in previous research.

2.
Proposing a new voltage regulator design for the SVC to aid in damping SSR along with a design of a new damping controller. 3.
Validating the performance of the proposed controller through detailed time domain simulations when the system is subjected to frequency disturbances on the grid side, mimicking the action of speed governors, which has not been documented in the earlier research.

Dynamic Modeling of WT-SCIG and Controller
The single line diagram of the WT-SCIG system considered in this study is shown in Figure 1. The mathematical model of the induction generator and the network was developed based on [12] and is presented below. Following this, the mathematical models of the voltage regulator and the subsynchronous damping controller (SSDC) are presented.

SCIG Weak Grid
High X/R ratio Low SCR � Figure 1. Single line diagram of the WT-SCIG system with variations in grid voltage and frequency.

Dynamics of Induction Generator
The detailed dynamic model of the SCIG is presented in this section. The dynamics of the stator and rotor windings are described by the following equations: where X s = X ls + X m , X r = X lr + X m , s g is the slip of the induction generator, which is given by s g = 1 − ω r , and ω s is the synchronous speed in rad/s. The electromechanical dynamics of the wind turbine are described by the swing equation given as the following: where T e = X m (I qs I dr − I ds I qr ) is the electrical torque output of the wind turbine, T w is the torque input from the wind and H g is the inertia constant of the wind turbine in seconds. It is to be noted that Equations (1) and (2) are nonlinear.

Dynamics of Transmission Line, Shunt and Series Capacitors
The dynamics of the series capacitor compensated transmission line connecting the terminals of the SCIG to the grid is given by the following equation: 1 ω s X t 0 0 X t d dt Shunt capacitors are installed at the terminals of the SCIG to supply the necessary reactive power to maintain the terminal voltage of the generator at an acceptable level. Equation (4) determines the dynamics of voltage across the shunt capacitor. Depending on whether the shunt compensation is fixed or variable, the susceptance of the shunt capacitor B sh = 1/X sh is constant or is the output of the voltage regulator of the SVC, which is described in more detail in Section 2.3.
Voltage across the series capacitor is described by the following differential equation:

Design of Voltage Regulator
The block diagram of the voltage regulator for the SVC is given in Figure 2. The terminal voltage of the generator is sensed through a transducer, which has a time constant of T m1 seconds and is compared with the reference voltage. The voltage error is processed by a PI controller, having proportional and integral gains K p1 and K i1 , respectively. The output of the PI controller gives the reactive power reference (Q re f ) for the shunt capacitor. This is compared with the actual value of the reactive power injected by the shunt capacitor (Q sh = V 2 m B sh ) at the generator terminals, which is measured by a transducer with time constant T m2 seconds. The reactive power error is processed by a second PI controller with proportional and integral gains K p2 and K i2 , respectively. The output of the second PI controller gives the susceptance B sh of the SVC. The proportional gain of the second PI controller enhances the stability domain (up to 95% series compensation) of the system with measurement transducers with higher values of time constants, which is shown in Section 3.3.3. The presence of a second PI controller is helpful in tuning the response of the system if there are limits imposed on the gains of the first PI controller. There is an additional signal V sup , which is a supplementary signal added to the first summing junction of the voltage regulator, which is used to control critical bifurcations in the system due to SSR. This is explained in more detail in the following section. The dynamics of the measurement blocks of the voltage regulator are described by the following equations: Equations (8)-(11) define the dynamics of the two PI controllers.

Design of the SSDC
The design of the subsynchronous damping controller (SSDC) consists of a washout block in series with a compensator, which provides a phase lead. The structure of the SSDC is very similar to that of a power system stabilizer used in synchronous generators to eliminate negative damping introduced by high gain exciters [13,14]. The signals X sig in Figure 3 considered in this paper are the generator slip (s g ) and the real power at the generator terminals (P e ). The dynamics of the SSDC are given by the following equations: 1 + 1 + 1 +

Results
Eigenvalue analysis and time domain simulations were performed on the system under study and the results are presented in the sections that follow. The parameters of the induction generator, network and wind turbine can be found in Appendix A.

Stability Domain Analysis of WT-SCIG with Fixed Shunt Capacitor Compensation
Bifurcations are identified in the system under study by computing the eigenvalues of the linearized mathematical model. The critical bifurcations encountered in this study are the Hopf bifurcation and the saddle-focus [15] bifurcation, which result in a change in the stability of the equilibrium point. In addition to this, node-focus bifurcations were also observed in the power system, which may or may not cause a change in the stability of the equilibrium point and is, therefore, non-critical [16]. Stability domain of the system with respect to series compensation is defined by the Hopf bifurcation point at which SSR occurs.

Eigenvalue Trajectories for Different Generator Slips
Bifurcations in the system with fixed shunt capacitor compensation at the SCIG terminals are analyzed, considering the degree of series compensation K c = X c /X t as the bifurcation parameter. Eigenvalues of the system with an X/R ratio of 30 are calculated and plotted on the complex plane when K c is varied from 0.005 to 0.98 for two different generator slips s g = 1% and s g = 3%, which correspond to low and high value wind speeds, respectively, as shown in Figure 4. The variations in torque input from the wind are assumed to be slow and hence are neglected in this analysis (∆T w ≈ 0). This greatly simplifies the computations, as the operating slip s g does not have to be calculated for each level of series compensation. Eigenvalues are presented for the system for 60% series compensation in Table 1 for a slip of 1%. It can be observed from the table that the system consists of 5 pairs of complex conjugate eigenvalues and 1 real eigenvalue. Among the complex conjugate eigenvalues, there are 3 pairs, which have frequencies above 60 Hz (Supersynchronous modes λ 1 , λ 2 , λ 3 ), and 2 pairs have frequencies below 60 Hz (Subsynchronous modes λ 4 , λ 5 ). The complex eigenvalues corresponding to λ 4 has a positive real part resulting in an unstable system making it the critical eigenvalue pair. The participation factors of the critical eigenvalue pair are presented in Table 2.
It was observed from the eigenvalue analysis that the supersynchronous modes are stable for the entire range of series compensation considered for s g of 1% and 3%. However, the damping of the subsynchronous mode decreases as K c increases, and when K c is between 0.575 and 0.58, a Hopf bifurcation occurs due to SSR when s g is 1%. However, at a higher slip s g of 3%, the system is unstable for K c = 0.005 due to the presence of a real eigenvalue on the right half of the complex plane, resulting in the system being unstable. The equilibrium is called a 'saddle-focus', characterized by three leading eigenvalues given by λ 1,2 = −8.8894 ± j16.957 and λ 3 = 3.5009. The participation factors for these leading eigenvalues are presented in Table 3.    As the compensation increases, the unstable eigenvalue moves to the left half of the complex plane, and when K c > 0.315, λ 3 mode crosses over from the right half to the left half of the complex plane, resulting in a bifurcation, making the system stable. The system remains stable as the percentage series compensation (%K c = K c × 100%) increases to 60% (K c = 0.60) beyond which a Hopf bifurcation results, making the system unstable due to SSR. The critical percentage series compensation at which instability due to Hopf bifurcation is observed is slightly larger for a higher operating slip of the induction generator.

Eigenvalue Trajectories for Different Grid Strengths
Increasing the X/R ratio of the line decreases the Short Circuit Ratio (SCR) at the point of interconnection (POI), as the Thevenin impedance of the system up to the POI increases as the reactance of the line increases. This results in a weak interconnection to the grid. The eigenvalues of the system for a generator slip of 1% are calculated for X/R ratios of 15, 20 and 25 and are plotted in Figure 5. Supersynchronous modes are stable for the entire range of series compensation considered but the subsynchronous mode becomes less damped as the degree of series compensation is increased. The non-oscillatory mode is stable for the entire range of compensation. SSR occurs as the percentage series compensation approaches 71.8%, 65.6% and 61.2% for X/R ratios of 15, 20 and 25, respectively, which shows that as the grid becomes weaker, there is a higher possibility of SSR at lower values of series compensation. This is further evident in Figure 6, where the critical value of series compensation is plotted as a function of X/R ratio of the network, which gives the stability boundary with respect to series compensation.

Stability Domain Analysis of WT-SCIG with SVC at the Stator Terminals
An SVC equipped with a voltage controller is used at the generator terminal in place of a fixed capacitor to supply the needed reactive power to maintain the voltage at the bus at V re f . The voltage controller introduces 4 new states in the analysis: V m and Q m , which correspond to measured signals of terminal voltage and reactive power, and W 1 and W 2 , which correspond to the intermediate states of the PI controllers. The parameters of the SVC voltage regulator for this system under consideration (unless otherwise specified) are presented in Table 4.
1.05 pu 1 ms 1 ms 0.1 0 40 40 Eigenvalues of the system are plotted on the complex plane as K c is varied from 0.001 to 0.95 for different operating conditions. It should be noted that the equilibrium point corresponding to the slip of the generator does not change as the series compensation is varied because the generator terminal voltage is held constant by the SVC.

Eigenvalues for Low and High Wind Speeds
Eigenvalues of the system with an SVC at the generator terminal are plotted in Figure 7 for wind speeds of 14 and 18 m/s. When the wind speed is 18 m/s, the operating slip of the SCIG is 3.051%. For a small value of compensation K c = 0.001, 7 complex conjugate pairs of eigenvalues and one real eigenvalue all on the left half of the complex plane are observed. Among the complex conjugate pairs of eigenvalues, 4 of them are of low frequency: λ 1 = −974.04 ± j14.893, λ 2 = −12.633 ± j17.108, λ 3 = −0.4911 ± j7.3045 and λ 4 = −18.911 ± j376.28, which follow different pathways on the complex plane as K c is varied. As the series compensation is increased, λ 1 starts moving toward the right half of the complex plane, and its frequency of oscillation decreases as well. When K c changes from 0.074 to 0.075, a node-focus bifurcation results from λ 1 splitting into 2 real eigenvalues: λ 11 = −975.15, λ 12 = −972.57, which then move in opposite directions as K c is increased. However, these eigenvalues remain in the left half of the complex plane as K c is increased to 0.95. As K c changes from 0.697 to 0.698, a second node-focus bifurcation takes place where λ 3 = −7.3577 ± j0.51245 transforms into 2 real eigenvalues: λ 31 = −7.6109, λ 32 = −7.136, and move in opposite directions. An increase in K c also causes λ 4 to move toward the right half of the complex plane; a Hopf bifurcation occurs in the system when K c changes from 0.742 to 0.743 as λ 4 , which is now −0.0118 ± j97.353, changes to 0.02075 ± j97.166, and the system becomes unstable due to SSR. When the eigenvalues of the system for lower wind speed of 14 m/s are analyzed, 6 complex conjugate eigenvalue pairs (stable foci) and 3 stable nodes are observed. Among the 6 complex conjugate pairs, 3 are low frequency modes: λ 1 = −4.2745 ± j9.4714, λ 2 = −8.2775 ± j14.712 and λ 3 = −18.917 ± j376.28. The 3 stable nodes are λ 4 = −1021.3, λ 5 = −925.66 and λ 6 = −46.858. As K c is increased, λ 1 moves further into the left half of the complex plane and the frequency of oscillation decreases; at K c = 0.559, λ 1 = −9.4476 ± j0.52247, which is a stable focus, changes into 2 stable nodes: λ 11 = −9.6715 and λ 12 = −9.2474, which is the result of a node-focus bifurcation. Eigenvalue λ 2 moves further into the left half of the complex plane as K c increases. When K c changes from 0.703 to 0.704, λ 3 = −0.0259 ± j103.85 crosses over to the right half of the complex plane and changes to λ 2 = 0.0057 ± j103.66 as a result of a Hopf bifurcation due to SSR. When K c increases, λ 4 and λ 5 move toward each other; when K c changes from 0.938 to 0.939, λ 4 = −979.87 and λ 5 = −968.33 collide and transform into a pair of complex-conjugate eigenvalues λ 45 = −974.14 ± j4.0936 as the result of a node-focus bifurcation. It should be noted that the system is already unstable, due to the earlier Hopf bifurcation. Additionally, λ 6 starts moving along the real axis toward the right as K c increases, but it still remains in the left half of the complex plane.

Eigenvalues for Different Grid Strengths
Eigenvalues are computed for various levels of series compensation for 3 different grid strengths, which are indicated by the X/R ratio of the line and are traced on the complex plane as given in Figure 8. The system is analyzed for a lower wind speed of 14 m/s because it is observed from the previous section that subsynchronous oscillations occur at low values of K c when the power generated is low. The system with X/R = 10 has a node-focus bifurcation as K c changes from 0.563 to 0.564 as the stable focus λ = −9.5344 ± j0.4132 transforms into 2 stable nodes each moving in the opposite direction as K c is increased. When K c changes from 0.857 to 0.858, the complex conjugate eigenvalue pair λ cr = −0.0528 ± j73.976 crosses over to the right half of the complex plane and becomes λ cr = 0.0069 ± j73.81 as a result of a Hopf bifurcation, which gives rise to unstable oscillations at a subsynchronous frequency. Similar behavior is encountered in the system for higher X/R ratios but at lower values of K c . The node-focus bifurcation occurs when K c changes from 0.558 to 0.559 and from 0.557 to 0.558 for X/R ratios of 20 and 30, respectively. The system undergoes a Hopf bifurcation when K c is close to 0.704 and 0.623 withλ cr = 0.0057 ± j103.66 and λ cr = 0.0085 ± j120.02 for X/R ratios of 20 and 30, respectively. This shows that as the strength of the grid decreases, the percentage series compensation at which a Hopf bifurcation occurs becomes smaller, making it crucial to implement a damping controller in order to avoid this phenomenon. The stability boundary with respect to series compensation and the X/R ratio of the network is shown in Figure 9, which confirms that the stability domain with respect to series compensation becomes smaller for weaker networks. However, the system with SVC has a slightly larger stability domain, compared to the system with fixed shunt capacitors at the stator terminals of the SCIG. (%)

Unstable operating region
Stable operating region Figure 9. Stability boundary of the system using SVC with respect to the X/R ratio of network.

Eigenvalues for Different Measurement Time Delays
Time delays in measurement of voltage magnitude and reactive power are modeled using a first order transfer function. Eigenvalues of the system are traced when K c is varied from 0.001 to 0.95 and are plotted on the complex plane for 3 different measurement time delays in the voltage regulator which can be found in Figure 10. The system is analyzed for higher wind speed V w = 18 m/s because certain unique bifurcations occur in this system at this wind speed. When the measurement time delay is 1 ms, the system undergoes a node-focus bifurcation at K c = 0.075 and K c = 0.704 and a Hopf bifurcation due to SSR at K c = 0.743. With a measurement delay of 10 ms, the system undergoes a Hopf bifurcation at K c = 0.679 and a node-focus bifurcation at K c = 0.691 as a stable focus transforms into 2 stable nodes but the system is already unstable at K c = 0.679. As the measurement time delay is increased to 100 ms, for K c = 0.001, there are 2 pairs of complex conjugate eigenvalues λ 1 = 1.4651 ± j19.368 and λ 2 = 0.2911 ± j9.7375 on the right half of the complex plane, which results in the system being unstable. When K c increases, these two pairs move toward the left half of the complex plane and when K c = 0.038, λ 1 = 1.1153 ± j19.34 and λ 2 = −0.0075 ± j10.014. Though the real part of λ 2 is negative, the system still is unstable, as λ 1 has a positive real part. When K c approaches 0.172, λ 1 changes to −0.0068 ± j19.208, which means all the eigenvalues of the system have negative real parts, indicating a stable system. Since the system changes from an unstable state to a stable state because of the movement of a complex conjugate pair of eigenvalues to the left half of the complex plane from the right, the system undergoes a reverse Hopf bifurcation.
The reason for this phenomenon can be explained as follows: the system is unstable for small values of series compensation, due to increased time delays in measurement. When the series compensation is increased, the system becomes more and more stable because increasing the series capacitive reactance decreases the effective line reactance, which decreases the time constant of the electrical circuit, thus compensating for measurement delays. Hence, the system can be stabilized by increasing the series compensation above a certain threshold when instability due to large measurement delays is experienced. The upper limit of K c is determined by the Hopf bifurcation, which occurs due to SSR at K c = 0.683 where λ cr = 0.0223 ± j108.57. Finally, a node-focus bifurcation takes place at K c = 0.865, where a stable focus changes into 2 stable nodes but the system is already unstable due to SSR.

Stability Domain Analysis of WT-SCIG with SVC and SSDC
In this section, the performance of the designed damping controllers based on slip signal (SSDC-A) and real power (SSDC-B) is analyzed through computing eigenvalue trajectories as K c is varied from 0.001 to 0.95 for various X/R ratios, wind speeds and measurement time delays. The controller gains of the SVC and SSDC are tuned for the various scenarios considered to obtain a stable system characterized by all eigenvalues on the left half of the complex plane up to a series compensation of 95%. The parameters of the network and wind speeds are the same as those considered in Section 3.2.

Eigenvalue Trajectories with SSDC for Low and High Wind Speeds
Eigenvalue trajectories of the system for an X/R ratio of 20 with SSDC-A are computed for wind speeds V w of 14 and 18 m/s and are presented in Figure 11 on the left. The subsynchronous mode starts moving toward the right half of the complex plane but even for a compensation of K c = 0.95, the corresponding eigenvalue pairs are λ sub = −3.2013 ± j79.178 and λ sub = −8.2223 ± j72.016 for wind speeds of 14 and 18 m/s, respectively, which have negative real parts. The eigenvalue closest to the imaginary axis is λ = −0.00276, which corresponds to V w = 18 m/s. Similarly, for the system with SSDC-B, eigenvalue trajectories are computed for wind speeds of 14 and 18 m/s and are plotted in Figure 11 to the right. It is observed that the damping of the subsynchronous mode starts reducing as the amount of series compensation K c is increased from 0.001 up to 0.793; beyond that, the damping starts increasing as the compensation is increased. A similar phenomenon is observed for V w = 18 m/s, except that the damping of the subsynchronous mode decreases for K c up to 0.725; beyond this, the damping increases. However, all the eigenvalues of the system for both low and high wind speeds lie on the left half of the complex plane, which shows that the Hopf bifurcation is eliminated for the entire range of series compensation considered. The parameters of the SVC voltage regulator and SSDC-A and B are listed in Table 5. Eigenvalues of the system with SSDC-A and SSDC-B are traced on the complex plane for varying grid strengths, which is a function of the X/R ratio. It was shown in Section 3.2.2 that as the strength of the grid decreases, the percentage series compensation at which the system becomes unstable due to SSR becomes smaller. In Figure 12, it is shown that with SSDC-A and SSDC-B, the eigenvalues of the system are all found to lie on the left half of the complex plane, thus eliminating the Hopf bifurcation due to SSR. In SSDC-B, with controller gains K p1 , K p2 = 0.1 and K i1 , K i2 = 50 and a stabilizer gain K SS of 0.3, all the modes are found to be well damped for X/R ratios to about 15. Beyond that, the stabilizer gain is increased to 0.7 and all the modes are found to be well damped for X/R = {10, 20, 30}, and those are shown in Figure 12b. The values of control gains of the SVC and SSDC are presented in Table 6. Re-tuning of the stabilizer and PI controllers in the voltage regulator are not needed with SSDC-A for all the X/R ratios considered here.  Eigenvalues of the system with SSDC for 3 values of measurement delays (1, 10 and 100 ms) are plotted on the complex plane and are presented in Figure 13. The parameters of SSDC-A and -B and that of the SVC with SSDC-A are same as those presented in Table 5.
For the system with SSDC-A, all the modes are damped and no re-tuning of the PI controls of the voltage controller is needed. However, with SSDC-B, the controller gains of the SVC have to be re-tuned in order to obtain a stable system up to 95% series compensation for each of the time delays considered here; those values are presented in Table 7. Hence, it is shown that both SSDC-A and SSDC-B are able to mitigate subsynchronous oscillations as well as the unstable low frequency modes, which appear for a time delay of 100 ms in the system without the damping controller.

Time Domain Simulation Results
The detailed non-linear model of the system with and without the designed SSDC is simulated and the system variables are plotted. The disturbance applied in the system is a three-phase fault at 1 s on the grid side, which causes a temporary voltage reduction of 15%. The fault is cleared in 6 cycles. In addition to this, a frequency disturbance is simulated using a mathematical function given by the following: where τ is the damping term and ω n is the frequency of the disturbance. Here, τ = 2.5 and ω n = 2 rad/s which is of the order of frequency involving electromechanical oscillations in the grid. This mimics the frequency changes in the grid due to the response of speed governors to disturbances in the grid. An X/R ratio of 20 is considered with a percentage series compensation of 71%, which is slightly above the Hopf bifurcation point due to SSR in the system without SSDC at a wind speed of 14 m/s. Parameters of the SVC voltage regulator for the system without SSDC can be found in Table 4, whereas, for the system with SSDC, the parameters of the voltage regulator and SSDC are listed in Table 5.
The terminal voltage of the generator has growing oscillations due to SSR, which is seen in Figure 14. If unchecked, this will result in dangerous values of overvoltages, which is harmful to the system. With the SSDC, the oscillations are damped, which can be observed in Figure 14a,b.
The traces of real and reactive power delivered to the grid when the system is subject to disturbances are shown in Figures 15 and 16, respectively. Without the SSDC, real power delivered to the grid starts oscillating between negative and positive values past 2 s (Figure 15a), implying that power flows back and forth between the SCIG and the grid. The same is observed in Figure 16 with the reactive power delivered to the grid. However, in Figures 15 and 16, the oscillations due to SSR are damped by the proposed controllers.   Oscillations in the stator and rotor currents result in growing oscillations in the electrical torque produced, which causes a mismatch in the mechanical torque input from the wind and the electrical torque output of the generator, which is manifested as oscillations of increasing magnitude in the speed of the rotor. However, the system with SSDC returns to the same operating point after a few cycles, which is observed from the dynamic torque-speed characteristics shown in Figures 17 and 18.

Conclusions
An SCIG-based wind turbine with fixed shunt capacitors at the stator terminals connected to the grid through a series compensated line becomes unstable because of a Hopf bifurcation due to SSR for both low and high wind speeds. The percentage series compensation at which the Hopf bifurcation occurs in the system is greater for higher wind speeds. This is because as the wind speed attains higher values, the magnitude of slip is higher, resulting in a less negative rotor resistance; hence, the system requires a higher percentage of series compensation to become unstable. There is a possibility of nonoscillatory instability for small values of percentage series compensation at higher wind speeds due to the presence of a saddle-focus equilibrium point, resulting in an unstable node. This makes the operation of the system impossible without a minimum level of series compensation when the wind speed is high. An SVC used in the place of a fixed capacitor is able to eliminate the unstable node at small values of series compensation at high wind speeds. In addition to a Hopf bifurcation, a node-focus bifurcation is also introduced into the system, which may not affect the stability of the system. In addition, the percentage series compensation at which SSR occurs is increased with the SVC. The Hopf bifurcation point, which defines the stability domain with respect to series compensation, is affected by the strength of the grid in systems with fixed capacitors as well as an SVC equipped with a voltage regulator. As the strength of the grid decreases with an increase in the X/R ratio, the stability domain with respect to the series compensation becomes smaller. This is due to the fact that as the X/R ratio of the line increases, for a given reactance, the resistance of the line becomes smaller. This causes the net resistance of the system to become negative due to the 'Induction Generator Effect' of SSR at lower values of percentage series compensation. However, the Hopf bifurcation point can be shifted in the system with SVC by proper adjustment of controller parameters in the voltage regulator. With higher time delays, there is a possibility of unstable low frequency oscillations, which are eliminated by higher levels of series compensation. Hence, stable operation of a system with higher measurement delays might not be possible without a minimum level of series compensation. A supplementary damping controller is needed for stable operation of the system for high levels of series compensation. A damping controller based on the generator slip signal as well as the real power at the terminals is able to successfully damp subsynchronous oscillations when tuned properly. This paper highlights new bifurcations in a WT-SCIG with a weak interconnection to the grid, considering a series compensated line, using a fixed capacitor and an SVC. The existence of a saddle-focus equilibrium point is shown in the system with fixed shunt capacitor compensation. An improved design of a voltage regulator for the SVC with a supplementary controller to mitigate subsynchronous oscillations is also presented. Eigenvalue trajectories are presented for the first time for variations in series compensation considering different wind speeds, grid strengths and measurement delays. It is demonstrated that the designed controller is successful in damping subsynchronous oscillations, which allows high levels of series compensation, enabling a higher amount of power transfer to the grid.

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

Abbreviations
The following abbreviations are used in this manuscript: