Transient stability analysis of DFIG-based wind farm-integrated power system considering gearbox ratio and reactive power control

Nowadays, integration of large-scale wind farms (WFs) into power systems is experiencing rapid growth. As this rapid integration can affect transient stability significantly, employing doubly fed induction generator (DFIG)-based wind turbines, which have shown better behavior regarding system stability, has attracted much attention. This research contributes to the literature by investigating the transient stability of the power system with increasing penetration of DFIG-based WFs. In the proposed framework, the current-balance form is utilized for the network equations, and in this way, transient stability is performed using time-domain simulation. According to the simulation results, when the rate of wind power generation exceeds 0.7 per-unit, the increasing trend of the critical clearing time (CCT) is reversed and the CCT decreases greatly with the increased wind power penetration. In addition, the reactive power compensation by DFIG, the gearbox ratio, the power system strength, and DFIG parameters are comprehensively investigated as effective parameters on transient stability. Since the rated rotor speed of DFIG significantly impacts the electrical torque and machine currents, the reduction of the rated rotor speed due to the change of the gearbox ratio has been investigated as one of the effective factors to improve the transient stability. The simulation results demonstrate the effectiveness of the proposed approach in improving power system transient stability.


Introduction
Electricity generation based on fossil fuel power stations can be considered as one of the largest sources of emissions.This challenge promoted the importance of renewable units integration as they can create a net-zero-emissions energy system [1].Emerging transmission and distribution systems based on renewable energies have many environmental and economic advantages compared to traditional systems due to the intermittent nature of their production sources [2][3][4].Also, due to the coordinated operation of many infrastructures, energy hub systems based on renewable energies increase energy efficiency in addition to reducing emissions [5][6][7][8][9].Meanwhile, wind turbines have penetrated more rapidly compared to other renewable resources [10][11][12].According to the latest report of the Global Wind Energy Council, the global wind industry had its second-best year in 2021, with almost 94 GW of capacity added globally, trailing behind the 2020's record growth by only 1.8%.This brings total installed wind capacity to 837 GW, a growth of 12.4% compared to 2020 [13].
Transient stability, which is defined system's ability to maintain synchronism in a disturbance, has been assessed by considering synchronous generators (SGs) in different studies in the literature.In a disturbance state in the system where the wind power penetration (WPP) is not considerable, wind farms (WFs) are disconnected from the power system when the fault occurs.However, due to the increasing the WPP in the power systems, the study of their behavior and interaction with other power generation equipment and loads seems important [14,15].Analyzing the performance of the power system when dealing with dynamic and transient phenomena is one of the most significant issues in this field.As a result, this topic has been noticed by researchers and various studies have examined the effects of the high penetration of doubly fed induction generator (DFIG)-based WFs on the dynamic behavior of the power system and transient stability in recent decades.
In recent years, transient stability analysis of WFs integrated power systems has received increasing attention in terms of power system topologies, wind power generation technologies, fault types, and location as well as innovative models, methods, and approaches [16][17][18][19][20][21][22][23].In this way, innovative frameworks follow the transient stability enhancement of a wind-penetrated power system through the power loss index method and the self-sorting analytical approach [16], Gauss-Hermite integral-based multi-point estimation method [17], two-stage model [18], supervisory adaptive predictive control [19], hybrid adaptive control of flywheel energy storage units [20], aggregated model [21], datadriven self-tuning additional sliding mode controller [22], and practical unit commitment framework constrained by the requirement of converter-driven stability [23].
Furthermore, from the existing wind power generation technologies, DFIG can be mentioned as a very attractive technology that has the most application in the industrial field.Accordingly, in the last decade, the transient stability analysis in power systems with the high penetration of DFIG-based WFs has become an interesting field of research [24][25][26][27][28][29][30][31][32][33][34][35].As the first performed research in this field, in [24], the simultaneous implicit (SI) based on time-domain (T-D) simulation has been used and a suitable model has been provided for transient stability assessment (TSA).Also, the researchers in [25] have also examined the effect of neglecting stator transients in the DFIG model.Comparison of the integration a common SG with a DFIG in power system when a 3-phase fault occurs, as one of the first works performed in this field, has been considered [26].Also, a comprehensive study has been performed in [27] that approves that DFIGbased WF integration improves the system response to small disturbances.
The current state of the art in this field focused on innovative models and approaches to improve transient stability in a WF equipped with DFIGs.Reference [28] uses a 9-bus system, integrated with a DFIG WFs, to analyze its behavior during disturbances, and [29] utilizes a 6-bus test system to analyze stability by stator-flux-oriented vector control approach for both stator-and rotor-side converter (RSC).A nonlinear control and unified power flow controller for transient stability improvement with transient energy function (TEF) also are proposed in [30].Semi-aggregation method through a nonlinear integral backstepping control [31] and active diode bridge fault current limiter topology by limiting the torque transients, currents peak, drop, active and reactive powers, and terminal voltage [32], also can be utilized.Further, a coordination approach to damp out the oscillations in a DFIG-integrated power system aiming at transient stability improvement using a power oscillation damper for static var compensator and power system stabilizer is developed in [33].Reference [34] employs DFIG by TEF technique and super-twisting differentiator in the multi-machine power system to improve transient stability.TEF is also employed by [35] to detect real-time transient stability.
Furthermore, the authors in [36] investigate the transient stability of the DFIG-integrated power system using the phase-locked loop and establish the grid-connected system model.Reference [37] analyzes transient stability in the hydro-quebec system containing a DFIG-based WFs through the multilevel static synchronous compensator.A control topology that integrates the robust performance of fractional order sliding mode control with the high charging/discharging rate and low maintenance requirements of super-capacitors employed in a DFIG-based wind energy system aiming at transient stability improvement [38].Reference [39] develops a noncontrolled fault current limiter model with reactive power support in a DFIG-integrated power system to enhance transient stability based on a nonsuperconducting reactor located on the rotor side.A new RSC control method is also proposed in [40], which utilizes DFIGs to improve stability and low-voltage ride-through (LVRT).
Considering various situations of the integration of DFIG WFs into the power systems, the reviewed studies have rarely been successful in compressive analyses of the transient behavior of the power systems in the presence of WFs, as the presented results in the above-mentioned literature are limited.The DFIG-based WF is very sensitive during the grid faults.In the event of fault, the terminal voltage of DFIG decreases expressively from the rated voltage, and the RSC is highly effected due to injects of high fault current.Therefore, the performance of the DFIG under these conditions can considerably affect system stability [41].In this way, it is very important to investigate a suitable method to enhance the stability and LVRT capability during the fault [42,43].Therefore, as the first work in this paper, the impact of Q DFIG and P DFIG has been investigated independently on improving the transient stability of the power system.According to the simulation results, when the rate of wind power generation exceeds 0.7 per-unit, the increasing trend of the critical clearing time (CCT) is reversed and the CCT decreases greatly with the increased WPP.Also, the generation of reactive power by DFIG has a desirable impact on transient stability.However, considering how DFIG supplies reactive power, the increased current in the rotor can be seen, leading to an increase in the converter rating.Thus, one of the solutions proposed in this paper to prevent the increase in the converter rating can be the use of external reactive power compensation devices in the DFIG bus.Furthermore, as the rated rotor speed in DFIG can be chosen arbitrarily by adjusting the gearbox ratio [44] and has a significant impact on the electrical torque and machine currents, the rated rotor speed is also analyzed as one of the factors affecting the transient stability improvement.Based on this, the reduction of the rated rotor speed due to the change of the gearbox ratio is investigated as one of the effective factors to improve the transient stability.In addition, the power system strength and DFIG parameters are comprehensively investigated as effective parameters on transient stability.
According to the above-mentioned points, the contributions of this paper can be summarized as follows: • Using the current-balance form for the network equations and development of the grid-connected DFIG equivalent circuit for the dynamic studies of the power systems.• Analysis of the effect of both P DFIG (by changing the capacity of the WF equivalent machine) and Q DFIG on transient stability improvement.• Analysis of the effect of change in the rated rotor speed of DFIG by adjusting the gearbox ratio on transient stability improvement.• Investigating the effect of the power system strength and DFIG parameters, specially DFIG stator resistance on power system transient stability.
The remaining parts of the paper are organized as follows: Section 2 provides power system components modeling and describes DFIG-based WF-integrated power system dynamic modeling.The simulation results are presented and discussed in Sect.3, and finally, conclusions are drawn in Sect. 4.

DFIG components modeling
In power system stability analysis, the low-frequency dynamics over a time frame of 0.1-10 s after a fault are of interest [45].Thus, fast dynamics that are damped out very quickly are ignored.Accordingly, in this paper, to analyze the transient stability, a two-axis model along with a static excitation system has been employed for SG [46].The components of DFIG are aerodynamic system, drive-train system, induction generator, and power electronics converters.In this way, all parts of the DFIG are modeled, separately.
The wind turbine aerodynamic system has been usually represented by an algebraic model, and thus, for stability studies, its mechanical dynamics are ignored, Eq. ( 1) [47]. where Also, in this paper, the blade angle pitch control has not been considered and the wind speed has been assumed to be constant, as in other works [11,12,27,35,48].At any given wind speed, there is an optimal wind turbine speed, in which the power extracted from the wind turbine is maximum (λ opt , C p, max ).Accordingly, by combining Eqs. ( 1) and ( 2), the equations of the wind turbine power and torque can be rewritten as Eqs.( 3) and (4) [11]. where The curve of speed-optimal torque (power) for DFIG is shown in Fig. 1.The rated wind speed is 15 m/s.According to Fig. 1 there are two different operating regimes for DFIG.In the sub-rated regime, the rated wind speed value is above the wind speed and the DFIG must follow the curve of maximum torque.However, in the rated regime, the rated wind speed value is below the wind speed and the speed and torque of the turbine must be limited to their rated values.The rated rotor speed can be chosen arbitrarily by adjusting the gearbox ratio, and the rated rotor speed has a significant impact on the electrical torque and machine currents, but due to the small amplitude, it has little impact on the rotor voltage [47].
Choosing higher rated rotor speed results in lower electrical torque because the torque is defined by Eq. ( 5): T e, rated = P ag, rated /ω r, rated (5) Therefore, in terms of the converter rating, the higher value of rated rotor speed is better due to create less the amplitude of the machine currents.Thus, in the present paper, the rated rotor speed of DFIG has been considered equal to 1.2 per-unit.Also, according to experimental studies and power system simulations, a two-mass model is preferred for the drive-train system of the wind turbine when studying the stability [11,47].Thus, the two-mass drive-train model is used for DFIG in this paper.
Normally, the stator transients in asynchronous generators and SGs have not been considered for dynamic analysis of the power system.Therefore, the currents and voltages are represented as phasors in the transient stability simulation programs [49].Accordingly, the stator transients are also ignored in this paper, and thus, the dynamic equations of the induction generator are expressed as follows (Eqs.(6-10)) [50].
V dr (9) where T e = e d i ds + e q i qs , Also, the electrical control system of DFIG consists of two parts: the RSC and the grid side converter (GSC).The objective of the RSC is to limit the rotor speed in the rated operating regime and to maximize the power extracted of the wind turbine in a sub-rated operating regime.On the other hand, the objective of the GSC is to always control the power factor and to keep the dc-link voltage constant [11,47].In terms of speed, the inner-loop current controllers of the power electronic converters are so fast compared with laboratory results, and theoretical considerations have shown that the speed operation of is much faster than the electromechanical transients studied in the stability analysis of the power system, and thus, their operation does not have a significant effect on the power system transient stability studies.In this way, a new reference current can be generated in 10 ms or

Sub-synchronous
Rated regime Sub-rated regime (1) (2) (3) Fig. 1 Typical reference power (electromagnetic torque) based on the rotor speed less, and it should be noted that 10 ms is a normal time step in power system dynamic simulations [12,47,51].As a result, change in the current controllers can be considered instantaneous.
Accordingly, considering that the GSC controls are done instantaneously, the AC current of the GSC is such that the active power injected into the grid matches that of the RSC at the unity power factor.Therefore, the active power transfer of the rotor is made completely by the GSC, and the DC-link voltage remains constant [12].Accordingly, it is assumed in this paper that the converters are ideal (the GSC is modeled as the current source and the RSC is modeled as the voltage source) and the DC-link voltage is constant.
According to what is said above, the model used for DFIG has been based on two basic assumptions: (1) electromagnetic transients in the stator and the branch between the inverters and the grid can be ignored.This assumption has no significant impact on electromechanical oscillations and helps to solve stator equations together with grid equations, and (2) it is assumed that the current control can be considered instantaneous, while the dynamics in the current control loops are ignored.The time constant of the current control is about one cycle (10 ms), which is quite faster than the time scale of transient stability studies (in the order of seconds).Therefore, the fast dynamics of such a control system have no impact on electromechanical oscillations.

The grid-connected DFIG modeling
When the way of connecting the machine to the grid is expressed, the network equations are extracted, and dynamic modeling is completed.It should be noted that the transient stability simulation is generally simplified using the differential-algebraic equation (DAE) model with network equations Fig. 2 The vector diagram of the machine (SG and DFIG) and grid reference frame in the current-balance form because the network-admittance matrix should be refactored only when a disturbance is occurring [46].In the same direction, Fig. 2 shows the vector diagram of how the synchronous reference frames of machines (the reference dq-frame) and the grid (the reference DQ-frame) are connected.Figure 2a shows the terminal voltage of the machine (SG) that is expressed in terms of the reference DQ-frame (Eq.( 11)) [46].
Accordingly, the DAE model of the network equations in the current-balance form can be expressed as Eqs.( 12) and (13) for SG [46].Similarly, Fig. 2b shows the connection between dqreference frame of DFIG and DQ-reference frame of the power system.Thus, the terminal voltage equation of DFIG is as follows (Eq.( 14)): By using the classical vector control method for DFIG modeling [11], the DFIG stator flux is aligned with the d-axis.Therefore, V s = V ds + j V qs , V qs = |V s |, V ds = 0, andδ = V s .In this way and according to the points expressed, the DFIG equivalent circuit in the grid-connected condition is presented in Fig. 3.
Therefore, the network equations in the current-balance form are obtained using the equivalent circuit developed for DFIG (Eqs.(15,16)).

Simulation results and analyses
In this paper, a single machine infinite bus (SMIB) system has been employed for TSA.To investigate the effect of WPP on power system transient stability, a WF is added to the base power system in the point of common connection (PCC) bus.
It should be noted that as with other work done in this area [11,12,14,35], the WF is modeled with aggregated wind turbines as an equivalent single machine.In this way and according to Sect. 2, the dynamic modeling of the power system can be represented using a set of DAEs.where u is the system's input and x and y, respectively, represent the state and algebraic variables.As a result, the function f represents the power system's differential equations, while the function g represents its algebraic equations.The equations should be solved step-by-step in the time domain for investigating transient stability using a T-D simulation, making the SI method one of the best approaches [12].In this approach, nonlinear differential equations are converted into nonlinear algebraic equations using trapezoidal numerical integration at each time-step, and the set of equations is then solved by the Newton-Raphson technique.It should be noted that the requirement for numerical simulation of a power system in the transient state is to calculate the initial conditions of the system variables that are obtained by load flow.Accordingly, employing the equivalent circuit of the DFIG and the power system, T-D simulation has been performed by MATLAB and the transient stability status is analyzed comprehensively.Figure 4 shows the power system under study that has been used in Ref. [12].The parameters and base case operation of the power system are presented in Appendix.A.

DFIG WFs integration analyses
It is assumed that because of the fault (a three-phase short circuit in the line X e2 ), the PCC bus voltage is decreased by 0.2 per-unit.In this paper, to analyze more accurately the impact of increasing WPP on the transient stability status, 3 cases have been studied.Thus, as the WPP in the power system increases, the capacity of the WF equivalent machine also changes and increases.For each case, the impact of increasing WPP on the power system transient stability is analyzed.The rate of wind power generation is assumed to be randomly variable.The rate of WPP is the rate of power generated by the wind relative to the total power delivered to the infinite bus in the base case.To perform TSA and T-D simulation, it is assumed that the time response of the system dynamics is simulated to be 20 s in the steady-state operating conditions and then the fault occurs.In this section, the fault clearing time for all cases is considered constantly to be 0.1 s.In this way, the simulations are performed and P DFIG and the bus voltage of DFIG, and the rotor angle and bus voltage of SG, are monitored.
For Case A, Figs. 5 and 6 show the simulation results.According to Fig. 6a, as long as the WPP in the power system is limited to 16.67%, the grid status is stable.However, when the WPP is increased to 19.44%, the grid experiences a multiswing rotor angle instability.Subsequently, the condition of the grid remains unstable with a higher increase (23.61% and 27.78%) in the WPP.
For Case B, simulation has been also performed in a similar way, and P DFIG and the bus voltage of DFIG, as well as the bus voltage and rotor angle of SG have been monitored, the results of which are shown in Figs.7 and 8.In this case, when the WPP increases to 33.33%, the condition of the grid is stable.However, when the WPP is increased to 38  47.22%, the grid experiences a multi-swing rotor angle instability.According to Fig. 8a, for the rate of WPP of 55.56%, the curve of the rotor angle is separately presented.It can be observed that the system will be unstable (multi-swing) in this state due to the rotor angle increasing.For Case C, the T-D simulation has been performed similarly by setting the fault duration of 0.1 s and the results  Fig. 8 The curve of the SG for Case B a the rotor angle and b the bus voltage obtained are illustrated in Figs. 9 and 10.According to Fig. 10a, when the WPP increases to a certain threshold (66.67%), the transient stability status is suitable.However, when WPP continues to increase, the grid experiences transient instability and severe voltage drop.
To analyze the results obtained, the CCT is obtained employing T-D simulation, as given in Table 1.According to Table 1, two points can be understood.First, for each case, when the WPP is increased, the CCT also is increased in comparison with the base case condition.Thus, the condition of transient stability initially improves.By continuing to increase the WPP up to a pre-specified threshold value, the CCT is still increasing.Then, the CCT is decreased with the continuation of the increase in WPP.Accordingly, the condition of transient stability initially improves with increasing the WPP, but then the situation subsequently deteriorates.The important point here is that in each of three cases, when the rate of the power generated by the WF is limited to about 0.7 per-unit (a certain threshold), the CCT increases with the increased WPP in the system.However, when the rate of wind power generation exceeds 0.7 per-unit, the trend is reversed and the CCT decreases greatly with the increased WPP.Second, due to the lower inertia of a DFIG compared to a SG, the higher the wind power generation in the grid, the higher the grid oscillations caused by disturbance.According to the last row of Table 1, when the rate of the power generated by the WF is limited to about 0.7 per-unit, first-swing transient instability occurs due to the insufficient synchronizing torque in the power system.However, when the rate of wind power generation exceeds 0.7 per unit, multi-swing transient instability occurs due to the insufficient damping torque.It should be noted that the rated rotor speed of DFIG is considered to be 1.2 per-unit, and it is assumed that DFIG has no reactive power exchange with the grid.

Reactive power change analysis
In this section, the effect of Q DFIG has been investigated on the transient stability status.Since voltage recovery after fault clearing by a DFIG is weaker than a SG, the value of Q DFIG is a vital issue in terms of TSA [11].To do this, in all cases, the rate of reactive power generated by DFIG is increased to 0.9 per-unit.It is assumed that the short circuit fault is the same as in Sect.3.1.Accordingly, CCT has been calculated based on TSA and T-D simulation.The simulation results obtained in this section are presented in Fig. 11.According to Fig. 11, for all cases, CCT increases when Q DFIG is increased; hence, the transient stability status has been improved.Also, as in the previous subsection, when the WPP increases up to a certain threshold value, the CCT also increases, but then the situation deteriorates.As Q DFIG increases, the stator and rotor currents increase.Accordingly, considering the rates of Q DFIG in three different states of -0.5, 0, and 0.5 per-unit, the curves of the rotor and stator current in the steady-state for the rated and sub-rated regimes are shown in Fig. 12.As Fig. 12 demonstrates, when Q DFIG of 0.5 per-unit, the stator current, especially the rotor current, increases significantly in both regimes.On the other hand, to decrease the power electronic converter rating, it is usually assumed that DFIG has no reactive power exchange with the grid [27,44].Thus, if the aim is to increase Q DFIG to improve its transient stability, the steady-state rotor current of DFIG and thus the converter rating must be increased.Generally, the reactive power can be injected into the grid by DFIG in two ways [44]: Scenario A: It is assumed that no portion of the reactive power is supplied by the rotor.In this way, the GSC operates at unity power factor, and all reactive power is transferred from the stator to the grid.As stated in Sect.2, the equivalent circuit presented for DFIG in this paper has been obtained Fig. 12 The magnitude of DFIG a rotor and b stator currents in the steady state accordingly.Scenario B: In this state, a part of the DFIG reactive power is supplied by the GSC, and therefore, the GSC does not operate at unity power factor.
In Scenario A, it is assumed that all the reactive power generated by DFIG that is equal to 0.9 per-unit is supplied by the stator, and the steady-state rotor current is considered.Here also, the conditions of the fault in the power system are the same as in the previous section.Figure 13 shows the rotor current of DFIG for all cases.As can be seen, the steady-state rotor current has increased between 21 and 33% for all three cases.Additionally, the lower the capacity of the WF, the greater the increase in rotor current due to the increase in the generated reactive power.Thus, the RSC rating should be increased in this state.On the other hand, in this state, all the generated reactive power is supplied by the DFIG stator and is injected into the grid.Therefore, the GSC operates at Reactive power = 0.9 Pu Igsc = 0.083 Pu Igsc = 0.076 Pu Igsc = 0.087 Pu Fig. 14 The steady-state GSC current for all cases with maximum power capture of the WF unity power factor.In Fig. 14, the GSC current for all cases is shown.Considering Fig. 14, the GSC current is very low in this state.Therefore, regarding the converter rating, this state of supplying the reactive power by DFIG does not increase the rate of the GSC.
In Scenario B, the GSC does not operate at unity power factor and supplies a part of the reactive power generated by DFIG in the operating conditions of reactive power control.Accordingly, the reactive power of the GSC is no longer equal to zero.
Therefore, in this state, the GSC current changes as Eq. ( 19) and is placed in the equivalent circuit of DFIG: Thus, it is assumed in this section that the reactive power of 0.3 per-unit is supplied by the GSC and the reactive power of 0.6 per-unit is supplied by the stator, totally constituting the reactive power of 0.9 per-unit generated by DFIG.Figures 15  and 16 show the RSC and GSC current.In this state, the RSC current increases by about 13% and includes a less percentage of increase than Scenario A, while the GSC current includes a more percentage of increase than in Scenario A and increases by 26.8% in Case A. Reactive power = 0.9 Pu Igsc = 0.252 Pu Igsc = 0.229 Pu Igsc = 0.268 Pu Fig. 16 The steady-state GSC current for all cases with maximum power capture of WF Reference [42] is one of the works conducted in this area.In [42], the impact of Q DFIG on the transient stability status has been investigated and it only shows that the generation of reactive power by DFIG has a desirable impact on transient stability.In this subsection, according to what is stated above, the impact of Q DFIG on the transient stability status improvement was comprehensively investigated.Based on the results of the simulations, the generation of reactive power by DFIG has a desirable impact on transient stability.However, considering how DFIG supplies reactive power in both scenarios, the increased current in the rotor can be seen, leading to an increase in the converter rating.Thus, one of the solutions proposed in this paper to prevent the increase in the converter rating can be the use of external reactive power compensation devices in the DFIG bus.

Rated rotor speed change analysis
The rated rotor speed can be chosen arbitrarily by adjusting the gearbox ratio, and the rated rotor speed has a significant impact on the electrical torque and machine currents, but due to the small amplitude, it has little impact on the rotor voltage.As mentioned earlier, in terms of the converter rating, higher value of rated rotor speed is better because it will create less the amplitude of the machine currents [44].Accordingly, in the all of previous sections, the rated rotor speed of DFIG is considered to be 1.Fig. 17 The stator and rotor currents and the electromagnetic torque with different rated rotor speed be noted that choosing higher rated rotor speed results in lower electrical torque.This is an issue that in the previous studies have not addressed it.Figure 17 shows the impact of changing the rated rotor speed by performing the steady-state simulation for DFIG.
As can be seen in Fig. 17, by decreasing the rated rotor speed to 1 per-unit, the currents of DFIG increase in the steady state and the condition deteriorates with respect to the converter rating.However, due to the impact of the rated rotor speed on the electromagnetic torque, the impact of changing the rated rotor speed of the DFIG on power system transient stability is studied in this section, which is considered the main contribution of this paper.For this purpose, the rated rotor speed is reduced from 1.2 per-unit to 1 per-unit by changing the gearbox ratio for all three cases and TSA is performed.Also, it is assumed that the fault of the power system is the same as in Sect.3.1, and the time response of the system dynamics is simulated to be 20 s in the steadystate operating conditions and then the fault occurs.The fault In the final part, the purpose is to investigate the effect of the DFIG parameters on transient stability by increasing the WPP.In the previous works, the effect of DFIG inertia as an effective parameter on transient stability has been investigated.It should be noted that DFIG inertia plays a crucial role in power system transient stability.Higher inertia provides enhanced system stability by contributing to the inertial response during disturbances and mitigate voltage and frequency deviations [44].However, according to the simulations performed, the effect of DFIG stator resistance on power system transient stability is analyzed in this paper.Accordingly, to perform the simulations in this section, all of the DFIG parameters are kept constant, and we only set the stator resistance of the DFIG in 3 different steps equal to 0.005, 0.01 and 0.05 .In this way, TSA is done for all cases and the CCT is calculated.Based on the simulation results in Table 2, the CCT is increased with increasing the value of the DFIG stator resistance; hence, the transient stability

Conclusion
A comprehensive analysis of the influence of the effective parameters on transient stability of the power system is covered in this paper.At first, this work has been carried out by considering the penetration of the DFIG-based WFs, the reactive power compensation by DFIG, and as well as the gearbox ratio as important factors affecting transient stability.The simulation results indicate that when the rate of the power generated by the DFIG WF is limited to about 0.7 per-unit, the CCT increases with the increased WPP in the system.However, when the rate of wind power generation exceeds 0.7 per-unit, the trend is reversed and the CCT decreases greatly with the increased WPP.Also, the simulation results show that reactive power compensation by DFIG has a desirable impact on transient stability.However, according to how reactive power is provided by DFIG, rotor current increases, leading to an increase in the converter rating.Thus, one of the solutions proposed in this paper to prevent the increase in the converter rating can be the usage of external reactive power compensation devices in the DFIG bus.Also, according to the simulation results in all studied cases, the transient stability power system condition is improved by decreasing the rated rotor speed from 1.2 to 1 per-unit.Finally, other parameters considered for TSA are the power system strength and the stator resistance of the DFIG.As the simulation results indicate, the CCT is increased with increasing the power system strength and the value of the DFIG stator resistance; hence, the transient stability condition improves.Future work is under way by the authors to further enhance the proposed approach.The authors intend to analyze the proposed approach in a multi-machine system and investigate its challenges in practical applications.

Fig. 3
Fig. 3 DFIG equivalent circuit in the grid-connected condition for dynamic studies

Fig. 4
Fig.4 The schematic of the SMIB system

Fig. 5
Fig. 5 The curve of the DFIG for Case A a the active power and b the terminal voltage

Fig. 7
Fig. 7 The curve of the DFIG for Case B a the active power and b the bus voltage

Fig. 11
Fig.11The impact of the increase of reactive power generated by DFIG on CCT

Table 1
The curve of the SG for Case C: a the rotor angle and b the bus voltage 2 per-unit.On the other hand, it should

Table 2
Comparison of the CCT with increasing the DFIG stator resistance

Table 3
The parameters of the SMIB system understudyThe SG and excitation system The SMIB system The DFIG wind farm