Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Carousel with three slides shown at a time. Use the Previous and Next buttons to navigate three slides at a time, or the slide dot buttons at the end to jump three slides at a time. Vibrating Wind Turbine

Yuxi Ji, Wei He, … Meng Zhan

Amirhossein Sajadi, Rick Wallace Kenyon & Bri-Mathias Hodge

Farah Echiheb, Yasmine Ihedrane, … Madiha ElGhamrasni

Gurupraanesh Raman, Gururaghav Raman & Jimmy Chih-Hsien Peng

Weijia Yang, Per Norrlund, … Urban Lundin

Samyak Tamrakar, Michael Conrath & Stefan Kettemann

MA Hannan, Jamal Abd.Ali, …ZY Dong

Huihui Song, Xuewei Zhang, … Yanbin Qu

Libor Pekař, Mengjie Song, ... František Zezulka

Scientific Reports volume 13, Article number: 2476 (2023 ) Cite this article

The doubly-fed induction generator (DFIG) with virtual inertia control and reactive damping control gives a renewable energy generation system inertia and damping characteristics similar to those of a thermal power plant, and the parameters of the control strategy have a direct impact on the small-signal stability of the system. This paper firstly introduces the operating characteristics and control strategies of DFIG-based damping control and virtual inertia control, establishes a small-signal model of the control-based DFIG integrated interconnected system, and investigates the effects of virtual inertia and reactive damping values on the small-signal stability of the system; then, the maximum damping ratio of the interval oscillation mode in small disturbance analysis is taken as the optimization objective, and the control parameters are the optimization variables. An optimization method of inertia and damping parameters is established for improving the small disturbance stability of the system. The results show that the optimization procedure could improve the damping ratio of the interval oscillation mode while ensuring the system frequency. The effects of virtual inertia and reactive damping values on the small signal stability of the system are investigated, and an optimal allocation model and method for virtual inertia used to improve the small disturbance stability of the system is proposed.

New energy generation is gradually replacing fossil energy generation in the grid. Power electronics connect most new energy generation to the grid and disconnect during operation. Renewable energy generation has minimal inertia. Thus, ordinary renewable energy producing equipment cannot regulate power supply frequency caused by grid oscillations1. Power system inertia decreases. This threatens power system frequency stability2.

Virtual inertia control is to use the kinetic energy stored in the rotor of the wind turbine to increase the electromagnetic power output of the inverter3 and slow down the frequency drop4 in a short time . In terms of ways to improve the control method, The author of5 proposes a general virtual inertia control strategy that can regulate the inertia power provided by the wind turbine by modifying the corresponding virtual inertia control coefficients. Some researchers have proposed an online under-frequency load shedding (UFLS) strategy considering virtual inertia control of wind turbines6. It can use the FM potential of wind farms to their fullest and optimize under-frequency load shedding by using fuzzy techniques to cut loads more accurately and speed up frequency recovery. The way that virtual inertia is spread out in the system can also affect how stable the system is. Some researchers have pointed out that the optimal placement of virtual inertia can be studied from two perspectives: frequency stability and small-signal stability. Some researchers have used the performance metric of the parametric number of parameters to show the stability of the system7. The optimal inertia layout problem is turned into a non-convex problem, and a way to find a locally best solution is shown. However, because the effect of random fluctuations in measurements during modeling is not controllable and the damping coefficient is not a decision variable in the problem setup8, the situation becomes complicated when considering extensions to more detailed system models and specifications. The paper9 proposes an optimization procedure to answer the question of how changes in inertia affect the stability of the power system and where it is most appropriate to place the device providing the virtual inertia. But the article doesn't take into account how the inertia system interacts with other control loops, and it's too complicated to solve the nonlinear problem. For practical purposes, a simpler approach is needed.

However, the wind turbine may affect the small disturbance stability of the system after the introduction of virtual inertia feedback link. The oscillation mode of a doubly-fed wind turbine after the introduction of virtual inertia has been analyzed. This has the potential to deteriorate the small disturbance stability of the system and to reduce the damping ratio of the inter-regional oscillation modes. In the case of insufficient damping, the inertia change of the system after disturbance may lead to severe transient frequency oscillations10. Many studies in recent years have shown that reactive power control with wavelet transform also plays an important role in improving the system stability. In11, the effect of active and reactive power control on the synchronous stability of the system was analyzed and it was found that reactive power damping control (RPDC) helps to counteract the side effects of virtual inertia control(VIC) on the damping of electromechanical oscillations. In12 author explored the role of RPDC through root trajectory analysis and simulation studies. The importance of RPDC loop parameter selection was pointed out, and it was noted that the integration loop of RPDC may be involved in system oscillations13. Some researchers investigated the enhancement of RPDC combined with VIC on the damping of doubly-fed motor oscillations14. However, the effect on the stability of small disturbances when the virtual inertia is increased is not discussed. According to some researchers, the optimal position and configuration amount of virtual inertia need to consider frequency stability and small signal stability15. In this paper, the role of RPDC and VIC in power system stability is investigated, and a procedure is used to calculate the optimal values of control parameters under the premise of ensuring system stability.

This paper first investigates the effect of damping and inertia of the dominant eigenvalues of small signal model. The damping and inertia values depend on the state of the control link in a DFIG, thus demonstrating the interaction between the damping and inertia performances and the system state. Based on these findings, updated rules that can improve the damping ratio of the grid are provided. The main idea is to use linearized summation methods to handle the sensitivity of dynamic eigenvalues of a nonlinear system to achieve an optimal system state equilibrium. The influence of variations in the parameters involved in the oscillation mode damping ratio on the sensitivity of the damping ratio is analyzed parametrically. Finally, by solving the proposed optimal configuration model of inertia and damping parameters, the optimal configuration scheme is obtained. The feasibility of the proposed optimization method is verified by simulation.

This section first explains the control model and its use including DFIG. Following that, the small-signal stability analysis methodology is presented.

A DFIG with virtual inertia control can provide a connected system with additional virtual inertia. Also, the damping characteristics of a synchronous machine can be simulated using reactive damping control. The control diagram is shown in Fig. 1.

DFIG grid-connected system with damping and virtual inertia control system.

A crucial component of DFIG grid-connected systems is multi-timescale control. The dynamic properties of the storage element (rotor) at this time scale may be studied using the rotor speed time scale. The analysis of the virtual inertia response using rotor dynamic analysis is suited for the rotor speed time scale, which can be utilized to analyze the dynamic properties of the energy storage element (rotor) at this time scale. The analysis of the virtual inertia response characteristics utilizing rotor dynamics is suited for the rotor speed time scale, which can be utilized to analyze the dynamic features of the energy storage element (rotor) at this time scale. As a result, in order to give a clearer and more specific solution to the issue of DFIG equivalent inertia assessment, this study A simplified model of the DFIG is built for the rotor speed time scale, and the following assumptions are established and the following assumptions are made.

Ignoring the dynamic changes in the magnetic chain and the stator resistance losses. Ignoring the variations in wind speed and mechanical torque.

Regulation of AC current control is ignored. The output current vectors of the RSC and GSC can track their reference values instantly

The net-side converter is less dynamic compared to the rotor-side converter, so the net-side converter dynamics are ignored.

Figure 2 illustrates a typical vector control strategy for a doubly-fed wind turbine based on end-voltage orientation. There are various forms and capacities of energy storage elements within the doubly-fed wind turbine, such as mechanical rotors, DC capacitors and AC inductors. For each energy storage element, conventional and transient control are set up for different purposes. In order to maintain a stable state of each type of energy storage element, conventional controls such as rotor speed control, DC voltage control and AC current control are designed respectively; in order to protect each type of energy storage element from its own stress, corresponding transient controls and hardware circuit protection are designed respectively, such as emergency pitch control/brake device, transient current command control, compensation control, etc. shown in grey dashed lines. The DFIG's output power varies with grid disturbance. According to storage capacity, each energy storage element's state changes. Depending on storage capacity, each energy storage element's state varies, driving each conventional/transient controller to act or process is multi-scale and sequential.

Typical control strategy and multi-scale division for doubly-fed wind turbines.

In this scenario, conventional control is employed for shallow faults and deeper faults. When the fault is shallow, conventional control is utilized; otherwise, transient control is used. Paper16 describes the reaction of a doubly-fed turbine in various modes. This work focuses on the control time scale of a doubly-fed wind turbine's rotor speed when the fault is shallow, i.e. when conventional control mode is applied. Modeling/analysis. Inertia dominates this scale's reaction. Speed, inertia, and reactive power management dominate at this size17. In addition, this paper only focuses on the occurrence of three symmetrical fault events in the grid, the negative sequence control under asymmetrical faults is not considered for the time being.

The doubly-fed turbine is continuously operated in the maximum power point tracking (MPPT) mode due to the control of the rotor-side inverter. The commonly used virtual inertia control is an inertia model based on the frequency feedback on the wind power generation in full power mode used to provide an inertial response to the wind turbine and adjust the wind power's active output18. This control model is shown in Fig. 3.

Virtual inertia control link block diagram.

The control link in Fig. 2 can be described by the following equation:

After the virtual inertia control loop is introduced, the active power input of the DFIG rotor-side inverter is given by:

where the introduced intermediate variable is \(\dot{x}_{\omega }\) ; Tω is the time constant; active power in the MPPT mode is Popt ; output power of the virtual inertia control link is PJ; virtual inertia control link scale factor is Kω, and it has a constant positive value.

To facilitate the analysis, the differential loop of the virtual inertia control link is simplified to a purely differential link19, and the active reference can be expressed as follows:

where the virtual inertia of a DFIG is provided by the virtual control link; the measured angular frequency fmeas = ωmeas/(2π); The equivalent inertia of the virtual inertia control link to the power system is

With the introduction of reactive power damping control, the reactive power generated by the control loop is given by20:

where Kd is the damping factor of the reactive damping control link; Qbase is the reactive power compensation generated by the system, and normally it has a value of zero.

Thus, doubly-fed motors with virtual inertia control and reactive damping control can provide virtual inertia and damping support to the grid. However, during electromechanical transients, the control loop may cause low-frequency oscillations due to the coupling of the doubly-fed motor with other synchronous generators.

In the power system analysis, it is considered that multiple DFIG wind turbines form a wind farm. The model is mainly used for small signal analysis of multi-generator power systems. In power systems, the synchronous generator (SG) is represented by a third-order model, where the excitation power unit is represented by a first-order model21. Assuming that all turbines in a wind farm have the same type and operating conditions, turbine parameters can be approximated to the calculated parameters of a single turbine22. Thus, the equivalent model of a wind farm can be described by the set of dynamic characteristic equations of a DFIG, while the doubly-fed generator model is shown in the Appendix. Its state variable x is x = [βref ωr ωgen θs isd isq Ed Eq udc x1 x2 x3 x4 x5 x6 x7 x8 δi ωi Efi E´ i ], where δi ωi Efi E´ i and are the first work angle, speed, excitation potential, and transient potential of the ith synchronous generator. The state variables that are often of interest when exploring how inertia control changes the coupling between the dynamic characteristics of a doubly-fed wind turbine and the system are as follows:

where δ1 and ω1 are the work angle and speed of the synchronous generator, respectively; s, θshaft and ωtur are the state variables of the mechanical loop of the doubly-fed wind turbine; the remaining variables denote the state variables of the induction motor and control link of the doubly-fed wind turbine. Near the equilibrium point, the differential equations of these state variables are linearized to get the characteristic matrix equation for the small-signal stability analysis of a power system with a wind farm that gets power from two different sources.

The relevant part of the system characteristic equation can represent the coupling relationship between the DFIG and SG dynamic characteristics.

where A11, A22, and A33 describe the dynamic characteristics of the SG, DFIG, and the mechanical part of the wind turbine, respectively; the other variables indicate the characteristics of the coupling of the parts with other different dynamic components.

Equation (11) can also be expressed as a state equation consisting of the state variables of the synchronous generator SG, the DFIG of the doubly-fed turbine and the mechanical part of the DFIG.

where A12 and A13 relate to the system tidal equation; they contain non-zero elements regardless of whether the equivalent load has dynamic characteristics23.

In this section, the conventional small-signal DFIG model and the small-signal model with virtual inertial and reactive damping control loops are presented. The expressions of the eigenvalues and damping ratios of the system under the action of damped and inertial control loops are derived based on the simple model of the two devices. In addition, an attempt is made to reveal the influence of the inertia and damping control loops on the small-signal stability of the system using the eigenvalue resolution method24. Finally, the eigenvalue analysis method commonly used to analyze the stability of small signals and their parameters, including the equation of state, eigenvalues, and participation factors of the components corresponding to eigenvalues, is introduced.

In the conventional DFIG control loop section, the time constant of the control link is much smaller than that of the mechanical section. The time constants of transient processes dominated by rotor currents are much smaller than the mechanical time constants15. Therefore, the dynamic process defined by Eqs. (1)–(8) is much faster than the dynamic process expressed by Eqs. (9)–(12). Thus, the steady-state solution to the dynamic process represented by Eqs. (1)–(8) is given by:

where Xm is the mutual stator-rotor inductance; Xs is the stator reactance; Ut is the stator voltage; Iqr is the rotor q-axis current; Psref is the active power.

In the maximum wind power tracking mode, the speed control module can rapidly provide a reasonable power output Psref based on the current wind speed and wind turbine speed. The total output power of DFIG is balanced with the wind power. According to Eq. (9), when there is an outward disturbance, the control link of wind turbines can rapidly adjust Iqr to balance the system inertia, which makes the dynamic mechanical characteristics of the DFIG decoupled from the system.

According to Eqs. (A9)–(A12) in appendix, only the rotor side current and the corresponding control parameter changes affect the damping parameters. This indicates that the additional damping control method based on the rotor-side current loop exhibits decoupling characteristics for other system parameter changes.

After introducing the virtual inertia control link and reactive damping control, the expanded state variables are expressed as follows:

The previous section analyzes the coupling characteristics of variables based on time constants of oscillatory modes. On the time scale, the rotor-side voltage has a high response rate. It is relatively decoupled from the other state variables and thus does not affect the small disturbance stability of the system25. Since the time constants of the inertia and damping controls are close to those of the transient electromechanical oscillation process, the control links of virtual inertia and system damping are the main factors affecting the small disturbance stability of doubly-fed generators.

This paper studies only the small-signal stability for small disturbances of inertia and damping. Therefore, parameters unrelated to the small-signal stability can be omitted or simplified. Thus, only the critical system nodes containing system parameters are analyzed, while nodes that do not have damping or inertia are removed from the model, so the model may not maintain the original network. Next, a small disturbance model containing dynamic system parameters is constructed.

For a doubly-fed wind turbine and a synchronous generator in the system, the network between them can be approximated to the components connected to the network, forming the system shown in Fig. 4.

Diagram of two areas power system.

where PG and PD are the active power transmitted via the contact line through SG and DFIG. E'∠ δ is the instantaneous voltage of SG. UG is the voltage value at the contact. U is the voltage on the busbar of other equipment, set as the ideal condition. Z1 and Z2 are the line impedance parameters. δ and θ denote the phase angle difference from the power inputs SG and DFIG to the contact point voltage UG, respectively.

The DFIG is set as the equilibrium node at the access point with the grid, the generator adopts the classical model, and the mechanical power PG is constant. Then the equation of motion of the node can be written as:

where HG, PG are the inertia and electromagnetic power of the synchronous generator, respectively; D is the damping factor.

The equation of perturbation motion of the equilibrium node is obtained by finding the amount of small perturbation to the above equation. The conventional unitary oscillation equation for SG can be expressed by26:

where DG is the constant of inertia and damping; ΔPG is the active power variation; Δδ is the power angle variation of an SG; p is the differential operator. Let the active power P of the DFIG cause the phase angle change θ and the reactive power Q cause the system voltage amplitude change ΔUG. Then, each variable is set to its initial value, denoted by the subscript "0." ΔPG and ΔQG can be expressed by:

By substituting (17) and (18) into (15), the system equation is obtained as follows:

Let the time constant of a PLL be infinitesimal; then, the frequency obtained by the virtual inertia control loop in a DFIG is the system frequency. Assuming that the angular frequency variation is consistent in the system presented in Fig. 3, the active output power variation with virtual inertia can be obtained by:

where Δωc is the measured angular frequency value. Since the reactive output affects only the voltage amplitude U at the common connection point, the reactive damping control loop can be used to denote QG.

The value of the reactive power damping control output derived from Eq. (4) is given by:

Further, combining Eqs. (22) and (20), we have:

The power variation ΔPu at the connection point is expressed by:

By combining Eqs. (24)–(26), it can be obtained that:

where \(K_{1} = \left( {\frac{{U\sin \left( {\delta_{0} - \theta_{0} } \right)}}{{Z_{2} }} - \frac{{E^{\prime } \sin \theta_{0} }}{{Z_{1} }}} \right)\frac{{Z_{1} K_{d} }}{{2U_{G0} - E^{\prime } \cos \theta_{0} }}\) , \(a_{0} = K_{G} + K_{L}\) .

Further, by substituting (27) into (20) can be expressed as follows:

where \(H_{DFIG} = \frac{{K_{\omega } a_{0} z_{2} }}{{2\pi U_{G} U\cos \left( {\delta_{0} - \theta_{0} } \right)}} > 0\) , \(a_{1} = \frac{{E^{\prime}U_{G} U\cos \theta_{0} \cos \left( {\delta_{0} - \theta_{0} } \right)}}{{E^{\prime}z_{2} \cos \theta_{0} + Uz_{1} \cos \left( {\delta_{0} - \theta_{0} } \right)}}\) , \(K_{QD}^{\prime } = \left( {\frac{{U\sin \left( {\delta_{0} - \theta_{0} } \right)}}{{Z_{2} }} + \frac{{K_{L} }}{{K_{G} }}\frac{{E^{\prime } \sin \theta_{0} }}{{Z_{1} }}} \right)\frac{{Z_{1} K_{Q} }}{{2U_{G0} - E^{\prime } \cos \theta_{0} }}\) , \(D_{DFIG} = K_{G} K_{QD}^{\prime } /a\) . where HDFIG and DDFIG denote the equivalent inertia and damping coefficients provided by the control loop, respectively; HDFIG is defined by virtual inertia control; DDFIG is defined by reactive damping control. In short, the virtual inertia control enhances the equivalent inertia of the system.

Eigenvalue analysis is an efficient method to analyze the damping characteristics of a system and the oscillation modes present in the system based on eigenvalues of the state matrix of the system. Its mathematical basis is the Lyapunov linearization method. The differential equation consisting of state variables near the equilibrium point is linearized to obtain the eigenmatrix equation for small-signal stability analysis of a power system containing a doubly-fed wind farm27.

In This study, the effect of a wind turbine on the small disturbance stability of a system after the introduction of inertia is observed. Using the mode analysis method, different oscillation modes dominated by certain system components can be obtained. The participation factors and damping ratios of different state variables are derived according to the corresponding characteristic roots, which indicate the degree of participation of the state variables in different system components in the oscillation mode and the influence on the damping ratio. In addition, the degree of influence of a specific component on the small-signal stability of a system is expressed. The specific method and participation will be explained in detail in the next section.

When the wind turbine load changes, the output torque imbalance reduces the resonance frequency and may cause low frequency oscillations in the output power of the wind power system, thus affecting the safety and stability of the wind turbine and its system. Among the main oscillations caused are:

The DFIG unit internal shaft system oscillation model with an oscillation frequency of 2.18 and a relatively low damping ratio of 0.04.

The intra-regional oscillation containing units G3 and G4 with an oscillation frequency of 1.44 Hz and an oscillation damping ratio of 0.12.

The intra-regional oscillation model containing synchronous units G1 and G2. The oscillation frequency and damping ratio of this model are 1.25 and 0.06 respectively.

Inter-regional oscillation model dominated by the control link of the wind farm. In the control coefficient variation has a significant effect on the oscillation frequency.

For the system defined by (20), the dominant eigenvalues of the system and the corresponding damping ratios are respectively given by28:

where mi is the inertia parameter in a system, di is the damping parameter, and a is the system electrical quantity parameter.

According to Eq. (30), the increase in inertia m causes the damping ratio of a system to decrease, which can affect the power angle stability of the system. An appropriate increase in the damping coefficient can improve the damping ratio of the system. It mainly affects the inter-regional oscillation mode, which indicates that grid-connected wind turbines cause a decrease in the small-signal stability of the system when virtual inertia control is used. To ensure the small-signal stability of the system, reactive damping control needs to be used in conjunction with virtual inertia control to achieve an optimal damping ratio of the system.

This section analyzes the effect of the doubly-fed wind turbine on system stability under small disturbances from the perspective of the system, considering the addition of inertia and reactive damping controls and how the added control loop changes the coupling relationship between the dynamic characteristics of the doubly-fed wind turbine and the system.

In small disturbance stability analysis, the system state-space equation has often been used. Still, the analytical expression of the mode damping ratio considering system parameters can be difficult to obtain by the state-space analysis. To solve this problem, a solution algorithm based on the sensitivity of the damping ratio to the system parameters is proposed in this paper. Since the DFIG incidental control loop parameters mainly affect the interval oscillation damping ratio, the main goal of optimization is to maximize the damping ratio for the worst case of interval oscillation, ensuring optimal attenuation in terms of system oscillations. It is assumed that both damping and inertia can be adjusted within constraints.

The damping ratio ζ is that it indicates the nonlinearity in sensitivity to the decision variables. The goal of an optimization algorithm is to find a local optimal solution by iteratively solving a linearized approximation. Specifically, the damping and inertial sensitivities of the damping ratio are used to compute a new damping ratio and sensitivity to the continuous parameter update of the damping ratio linearization in the next iteration, which it is necessary to have information on two parameters. First, calculating the sensitivity requires finding system eigenvalues, which must be recalculated at each iteration. In this study, an approximation of system eigenvalues is used. The second concern relates to the iteration approach. In this study, an iterative approach is adopted, fixing the step size of each parameter. The values of damping ratios are compared after the first and second iterations. The iterative algorithm stops if the new value is worse than the previous value or if the difference between them is less than the predefined constraint value. Although this approach has been proven to be useful in improving the stability of simple networks, employing more powerful solution strategies for complex systems might be necessary.

The superscripts "0," "ν," and "ν + 1" correspond to the first, current, and next iterations, respectively. The difference between the current and next iterations is expressed by:

For the optimized cost function, the constraints| are expressed as follows.

It has been shown that the inertia m has different effects on the small disturbance stability and frequency stability. Increasing the value of m can improve the frequency stability of the system, but it can also deteriorate the small disturbance stability of the system. To ensure the system has good small disturbance stability and frequency stability at the same time, the value of m should be determined considering the constraints of the two types of stability simultaneously.

Constraints (32) and (33) limit the total amount of inertia and damping of the system in the optimization process, respectively; meanwhile, constraints (34) and (35) denote the boundaries of damping and inertia, respectively. Further, step size constraints (36) and (37) limit the variations of parameters in each iteration. Namely, the optimization algorithm cannot converge in a nonlinear domain without such a constraint. Then, the sensitivity optimization is given by:

Equation (38) calculates a new damping ratio for the oscillatory mode using the damping ratio cumulative for the control parameters k and m sensitivity. The superscript v denotes the number of iterations and will not be described separately in this paper.

The optimization goal is to maximize the objective function under defined constraints. The termination condition is set to check whether the damping result obtained in two consecutive iterations is less than an arbitrary constant defined according to the accuracy requirement. In addition, if the number of iterations reaches the upper limit, the calculation is also terminated.

Based on the sensitivity of the damping ratio to the variable parameters, the algorithm used to solve the optimization model (31)–(37) is proposed, and its flowchart is presented in Fig. 5. This algorithm includes three main steps constraint setting, iterative solution determination, and termination condition definition, which are described in detail in the following.

Diagram of two areas power system.

It is necessary to determine conditions that satisfy the minimum inertia required by the system and the maximum inertia provided by the virtual inertia control link according to constraints (32)–(37). Then, the number of iterations is set to zero, and the initial values k0 and m0 are assigned to T(0); T(k) denotes the parameter of the mth wind farm at the kth iteration.

In this phase, ζ is first solved for with according to Eq. (28). As explained previously, the algorithm is used to determine the direction of optimization through sensitivity. In the (k + 1)st iteration, after the mth wind farm ζ needs to be determined whether it is less than zero, and if it is less than zero then it is set to zero. After the NDFIG wind farms have been obtained, it is judged whether the constraints are satisfied and if so, ζ(v) is updated to ζ(k + 1); if not, the kth result is retained.

Since the optimization objective is to maximize the damping ratio of the interval oscillation mode, then the termination condition can be set to be less than any constant chosen according to the accuracy requirements. In addition, if the number of iterations reaches the upper limit, the calculation is also terminated.

In this paper, the simulation system shown in Fig. 6 is built in DIgSILENT/PowerFactory. The diagram shows that the system is divided into two areas connected by transmission lines, with the left side as the sending end and the right side as the receiving end. Gen1 and Gen2 are synchronous generators on the sending side, and Gen2 is the reference motor. On the other hand, there are two wind farms, Wind Farm 1 and Wind Farm 2, on the receiving side. These farms have turbines fed from both sides and two synchronous generators, Gen3 and Gen4, which have a power output of 900 MW on the transmitter side. A study of how wind power with virtual inertia control affects the system performance in an interconnected grid found that when wind farms with virtual inertia control are on the transmitter side, they can make the oscillation pattern between regions unstable. According to the paper29, when the wind farm with virtual inertia control is at the sending end, the oscillation pattern between regions will be unstable. The theoretical derivation of the paper shows that the suppression of system oscillations is better when the virtual inertia is assigned far from the center of inertia or in the low inertia region. In the Appendix, the generator and excitation parameters, as well as the initial control parameters of the control system discussed in this paper, are given.

Simulation diagram of a four-generator dual-zone system.

The frequency stability of the power system is related to the unit inertia30. To ensure frequency stability, the minimum total inertia of the system should satisfy Eq. (38).

where Ht is the overall inertia of the system; HSG,n is the time inertia constant of the nth SG; Hω,m is the equivalent inertia of the mth doubly-fed wind turbine; SSG and SDFIG are the equivalent inertia values of the synchronous generator and doubly-fed wind turbine, respectively; Hω,m is the equivalent inertia of the mth doubly-fed turbine; SSG and SDFIG are the nominal capacities of the SG and DFIG, respectively; NSG and NDFIG are the nominal capacities of the synchronous generator and doubly-fed turbine, respectively; Hmin indicates the minimum inertia value for stable operation of the system, and it is expressed as follows :

This example's investigation starts with a sensitivity analysis to determine the impact of virtual inertia. The study first checks, through sensitivity analysis, the influence of the virtual inertia distribution on the interval oscillation pattern, and then validates the correctness of the suggested model and method. The study of this case first validates, via sensitivity analysis, the influence of virtual inertia allocation on the interval oscillation mode, and then checks the efficacy of the suggested model and method.

Let I1 and I2 be the sensitivity values of interval oscillation damping ratio ζ to inertia as m1 in area 1 and inertia as m2 in area 2. Combining pattern frequencies and participation factors, the interval oscillation modes of the system can be obtained. The sensitivity analysis enables a more intuitive insight into the effect of the virtual inertia m and the damping parameter d on the damping ratio ζ of the interval oscillation mode. There are 2 main factors that influence the distribution of virtual inertia. One is the electrical distance of the virtual inertia from the load and the other is the capacity of the wind farm where the virtual inertia control link is added. In this paper, we first consider the influence of inertia and the electrical distance from the load center on the sensitivity results. I1 and I2 are the sensitivities of the interval oscillation damping ratios for wind farm 1 and wind farm 2 respectively; the damping parameters are taken to be the same value. The capacities of wind farms 1 and 2 are the same, and L1 and L2 are the lengths of lines 1 and 2. Combining the mode frequency with the participation factor, the ratio of the inter-zone oscillation mode of the system can be obtained as follows.

As can be seen from Table 1, when two wind farms have the same capacity, keeping L2 length at 10 km, I1: I2 increases as L1 grows. Adding virtual inertia to wind farm 1 improves the inter-area oscillation mode damping ratio. Since the lengths of L1 and L2 measure the electrical distance from the wind farm to the load node, it can be seen that in this example, assigning the virtual inertia to a wind farm that is further away from the load electrical distance is more likely to The damping ratio of the grid interval oscillation mode is improved. When the line length L1 is greater than L2, the damping ratios of the interval oscillation modes of the two wind farms are related to I1 and I2 as follows: where the solid line indicates the damping ratio represent I1. The dashed line shows the damping ratio represent I2.

It can therefore be concluded that in this example system, the allocation of virtual inertia to wind farms with a long electrical distance from the load or a large capacity results in a better damping ratio for the inter-area oscillation modes and is more conducive to small disturbance stability of the system.

This section demonstrates the influence of wind farm capacity on the sensitivity results, proving that the main factor influencing the values of the wind turbine control parameters d and m is the wind farm capacity. Only one interval oscillation mode exists in this example system. The results of the sensitivity analysis are shown in Table 2.

When two wind farms have the same capacity, I1 and I2 increase with the growth of transmission lines L1 and L2, which indicates that adding virtual inertia to wind farm 1 will improve the inter-area oscillation mode damping ratio. In this case, the virtual inertia is assigned to wind farms that are further away from the load electrically, which improves the grid's inter-area oscillation mode damping ratio. In this paper, an optimal turbine configuration location in the system is selected directly to exclude the influence of the configuration of virtual inertia devices of the system on the small disturbances analysis.

Under the premise of guaranteeing the maximum virtual inertia of system capacity provided by the equipment in the system, as constrained by Eq. (39), the total capacity of the two wind farms is kept constant. The capacity ratio of the two wind farms is continuously changed; namely, the capacity of wind farm 2 is reduced, while the capacity of wind farm 1 is increased. In addition, the sensitivity ratio I1: I2 increases with the capacity share of wind farm 1. Thus, in this case, the degree of influence of virtual inertia on system stability under small disturbances is positively related to the capacity share of wind farm 1. Hence, virtual inertia to wind farms with a larger capacity can increase the inter-area oscillation damping rate.

Next we analyze the sensitivity of the damping ζ ratio to the fan damping control parameter K. The sensitivity equation is given by (21). Using the capacity parameters of the DFIG unit in Table 1 Case 4, the values of the sensitivity for each of the two regions are found, as shown in Fig. 7.

Comparison of sensitivity of oscillation mode damping ratio to damping coefficient.

According to Fig. 8, the sensitivity value of the oscillatory mode damping ratio to the damping factor of the system depends mainly on the magnitude of the inertia value. It indicates that the inertia value influences the configuration of the DFIG damping parameters in the system and is negatively related to the virtual inertia provided by DFIG. Also, Fig. 7 shows that the initial sensitivity of a region with a small capacity is considerable. Still, the sensitivity decreases fast with the virtual inertia, and the sensitivity of both areas tends to be the same. This indicates that the increasing effect of inertia on the damping ratio reaches a threshold value when the inertia increases to a particular matter.

Comparison of sensitivity of oscillation mode damping ratio to damping coefficient.

Next, we utilize the characteristics listed in Table 1 for Case 4 to calculate how sensitive the damping ratio is to the magnitude of inertia. Figure 8 displays the outcomes using the derivative of the damping ratio with respect to the damping factor d as the vertical axis and the inertia control factor m as the horizontal axis.

The solid line in Fig. 9 shows the sensitivity of the damping ratio to the virtual inertia in region 1, and the dashed line indicates the sensitivity of the damping ratio to the virtual inertia in region 2. According to the results in Fig. 8, the sensitivity values in both regions increase with inertia. However, the increase in the sensitivity value is more pronounced when increasing the virtual inertia in region 1. It can therefore be seen that in this example, when the damping coefficients are the same in both regions, the degree of influence of the virtual control coefficient on the stability of the system with small disturbances is positively related to the capacity share of wind farm 1, i.e. assigning the virtual inertia to a wind farm with a larger capacity in this example will improve the mode damping ratio.

Changing trends of inter-area oscillation.

According to the previous conclusions, the sensitivity of the damping ratio to the damping factor d depends mainly on the inertia value in a specific region of a system. In this paper, we selected the data of case 4 (Table 1) and the optimization algorithm (37) mentioned earlier is employed to determine an optimal damping ratio of the system and find the values of inertia and damping coefficient corresponding to optimal damping ratio. In the considered case, the minimum inertia of the system required for frequency stabilization is set to Hmin = 5.12 s, and the following constraint can be obtained based on Eq. (40):

Using the data of Case 4 (Table 1) converted according to the maximum virtual inertia output power provided by different capacities, it is set that m1 = 0.15, m2 = 0.25, d1 = 0.92, and d2 = 0.92. In the algorithm presented in Section III, the optimization step and termination condition are set to u = 0.1 and Δζ ≤ 1 × 10−8, respectively. The initial values and constraints are brought into the program for calculation and the optimization ends when 10 calculations are run and ζ reaches its maximum value.

The surface shown in Fig. 10 illustrates the relationship between the damping ratio of the interval oscillation mode and the control parameters m and d. The scatter trajectory in the figure represents the optimization iteration process. The results show that the critical eigenvalues of the optimized interval oscillation mode are shifted to the left in the real part and to the horizontal axis in the imaginary part, and the small disturbance stability of the system is improved. The correctness of the proposed algorithm is verified.

Comparison of eigenvalues after damping ratio optimization.

In the following, we will derive the relationship between the state variable participation coefficient of the system control loop and the damping ratio of the inter-area oscillations from the change in the eigenvalues derived in Fig. 9, where we derive the change in the participation coefficient. The corresponding characteristic roots of the inter-area oscillation pattern in Fig. 10 are selected and ranked according to the magnitude of the state variable participation coefficients. As the system parameters change over the course of the iteration, the state variable participation factors change accordingly. The results obtained are shown in Fig. 11.

Participation factors for state variables.

As shown in Fig. 11 the state variables of the control parameters dominate the oscillation mode. The virtual inertia control state variables xω1 and xω2 increase gradually with the parameters' values. The reactive damping control state variables participation factors xD1 and xD2 have the same changing. This is consistent with the previous conclusion that the inertia parameter assignment dominates the damping ratio. As the parameters increase, the damping ratio of the system increases, indicating that the system stability under small disturbances is improved.

The inter-area oscillation mode can be reflected in the dynamic response of the active power of the transmission line31. The amplitude of the oscillations reflects the small disturbance stability of the system. Due to virtual inertia control and reactive damping control, a doubly-fed generator is like to the synchronous generator in terms of inertia and damping characteristics. Namely, it can support the active system power for a short time after system power fluctuations. Assume that a three-phase short-circuit fault occurs at node B8 at t = 1 s, and the fault is addressed after 0.1 s. Then, power fluctuations at the contact line before and after optimization are as shown in Fig. 12.

Power across the contact line after a fault.

The curves in these plots show that the system amplitude and oscillation times become significantly larger without the addition of reactive control and virtual inertial control to the DFIG. The parameters of the virtual inertia control and reactive power damping control are set separately and compared with the parameters derived in “Parameter analysis and Simulation implementation” section that include both controls. When no auxiliary control measures are introduced into the system, the system oscillates and tends to destabilize due to an inadequate damping ratio after a three-phase fault has occurred. The system stabilizes only after several cycles when other stabilization measures have taken effect. The WTGs, with the introduction of virtual inertia in the region, reduce the active output during this period, reducing the amount of active power mismatch in the system during the transient process. Of the four cases, the values of the two control parameters derived using the method of sensitivity in this paper give the best results in terms of oscillation amplitude and oscillation time after fault occurrence. It also allows the reactive and active power of the doubly-fed wind turbine to maintain its rated output value to the maximum extent possible in the event of a fault. Simulation results show that the optimization method mentioned in this paper performs best in terms of oscillation amplitude and oscillation time of the transmission line, and the active power response curve decays to steady state the fastest, proving the effectiveness of the proposed optimization method.

In this paper, a virtual inertia optimization allocation method, which considers the small disturbance stability of the grid, is proposed. In addition, an optimal solution algorithm is developed considering the damping ratio sensitivity. The effectiveness of the proposed method is verified by time-domain simulations. This paper demonstrates that new techniques, such as virtual inertia control in systems with high proportional turbine penetration, can allow a system to adjust inertia and damping levels. The results indicate that the system damping ratio decreases with the inertia coefficient, thus affecting the system power angle stability. Sensitivity to parameters using damping ratios, justifying configurations with larger inertia and damping in wind farms with high capacity. Finally, the proposed selection method of optimal parameters for virtual inertia and damping can be used to maximize the small disturbance stability of a system under constraints. The arithmetic examples verify the effectiveness of the proposed selection method. The results show that the properly selected parameters can significantly improve the small disturbance stability while keeping the system's power angle stable.

Due to the limitation of the existing analysis methods, the effect of virtual inertia and its distribution characteristics on the small disturbance stability of a power system cannot be clearly explained. Therefore, this paper analyzes the influence of the virtual inertia distribution on the small disturbance stability of the system from an empirical point of view and provides a reference for the selection of capacity and parameters of the virtual inertia control. Based on the proposed optimization model and method, constraints and condition parameters can be extended. However, a quantitative description of the mechanism of the influence of damping and inertia control loops on the stability of small disturbances in DFIG containing power systems deserves further study.

The datasets generated and/or analysed during the current study are available in the [NAME] repository, [PERSISTENT WEB LINK TO DATASETS]. Reprints and permissions information is available at www.nature.com/reprints. The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials.

Tang, M., Li, Z., Li, G., Qie Z. & Yu, Y. “Power oscillation characteristics of large power systems under random disturbances. In: “2021 IEEE 5th Conference on Energy Internet and Energy System Integration (EI2), Taiyuan, China, pp. 2866–2870 (2021).

Kerdphol, T., Rahman, F. S., Mitani, Y., Watanabe, M. & Küfeoǧlu, S. K. "Robust virtual inertia control of an islanded microgrid considering high penetration of renewable energy. IEEE Access 6, 625–636 (2018).

Ndirangu, J., Nderu, J. & Irungu, G. “Impact of variation of virtual inertia and virtual damping on frequency and power angle responses of virtual inertia emulation controlled converter. In: 2022 IEEE PES/IAS PowerAfrica,Kigali, Rwanda, 1–4 (2022).

Tuo, M. & Li, X. Optimal allocation of virtual inertia devices for enhancing frequency stability in low-inertia power systems. In: 2021 North American Power Symposium (NAPS). 1–6 (2021).

Peng, Ke., Wei, Z., Chen, J. & Li, H. Hierarchical virtual inertia control of DC distribution system for plug-and-play electric vehicle integration. Int. J. Electr. Power Energy Syst. https://doi.org/10.1016/j.ijepes.2021.106769,128,(106769) (2021).

Wu, X., Dai, J., Tang, Y. & Xue, F. Adaptive under-frequency load shedding scheme for power systems with high wind power penetration considering operating regions. IET Gener. Transm. Distrib. 16, 4400–4416 (2022).

Untea, A. E. et al. Broiler meat fatty acids composition, lipid metabolism, and oxidative stability parameters as affected by cranberry leaves and walnut meal supplemented diets. Sci. Rep. 12, 21618 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Song, D. et al. Coordinated optimization on energy capture and torque fluctuation of wind turbines via variable weight NMPC with fuzzy regulator. Appl. Energy 312, 118821 (2022).

Gründel, L., Alipaz-Dicke, A., Trinh, M., Storms, S., Brecher, C. & Bitterolf, D. Frequency-based identification routine for the inertial parameters of an industrial robot. In: 2022 26th International Conference on Methods and Models in Automation and Robotics (MMAR). 352–357 (2022).

Zhang, W. et al. Online estimation of inertia and frequency regulation capability based on small perturbation. Energy Rep. 8(5), 650–656 (2022).

Liu, R. et al. An improved virtual inertia control strategy for low voltage AC microgrids with hybrid energy storage systems. Energies 15, 442. https://doi.org/10.3390/en15020442 (2022).

Fan, L., Yin, H. & Miao, Z. On active/reactive power modulation of DFIG-based wind generation for interarea oscillation damping. IEEE Trans. Energy Convers. 26(2), 513–521 (2011).

Yin, H. et al. Interaction mechanism analysis and additional damping control for hybrid multi-infeed HVDC system. IEEE Trans. Power Deliv. 37(5), 3904–3916 (2022).

Xu, X. et al. The Effect of DFIG on Electromechanical Oscillation with Combined Reactive Power Damping Control and Virtual Inertial Control. In 2019 IEEE 8th International Conference on Advanced Power System Automation and Protection (APAP) IEEE, 2019.

Xiong, L. et al. Specified time consensus control of ESSs for frequency support with DFIG wind turbines. Int. J. Electr. Power Energy Syst. 135, 107520 (2022).

Chang Yuanzhu, Hu. et al. Fault current analysis of type-3 WTs considering sequential switching of internal control and protection circuits in multi time scales during LVRT. IEEE Trans. Power Syst. 33(6), 6894–6903 (2018).

Tang, W. et al. Modeling of DFIG-based wind turbine for power system transient response analysis in rotor speed control timescale. IEEE Trans. Power Syst. 33(6), 6795–6805 (2018).

Jia, Y., Huang, T., Li, Y. & Ma, R. Parameter setting strategy for the controller of the DFIG wind turbine considering the small-signal stability of power grids. IEEE Access 8, 31287–31294 (2020).

Liao, K., Xu, Y., He, Z. & Dong, Z. Y. Second-order sliding mode based P-Q coordinated modulation of DFIGs against interarea oscillations. IEEE Trans. Power Syst. 32(6), 4978–4980 (2017).

Liao, K. et al. A sliding mode based damping control of DFIG for interarea power oscillations. IEEE Trans. Sustain. Energy 8(1), 258–267 (2017).

Chaban, A., Perzyński, T., Popenda, A., Figura, R. & Levoniuk, V. Mathematical modeling of transient processes in the susceptible motion transmission in a ship propulsion system containing a shaft synchronous generator. Energies 15(9), 3266 (2022).

Suvorov, A., Askarov, A., Kievets, A. & Rudnik, V. A comprehensive assessment of the state-of-the-art virtual synchronous generator models. Electr. Power Syst. Res. 209, 108054 (2022).

Jianzu, Hu., Lei, Yu., Chi, Y. & Tian, X. Analysis on the inertia and the damping characteristics of DFIG under multiple working conditions based on the grid-forming control. Energy Rep. 8(8), 591–604 (2022).

Jia, Y., Li, Y. & Li, W. Study on effect of DFIG on small signal stability of power system with high share wind power. China Int. Electr. Energy Conf. (CIEEC) 2017, 787–791 (2017).

Liu. W., Li. W., Wang. P., Li. D., Li. Z. & Xu, G. Transient Process Analysis of Multi-Physical Parameters of Turbine Generator During Out of Phase Synchronization Fault. IEEE Trans. Energy Convers. 37(3), 2058-2068(2022).

Ma, J. et al. Research on the impact of DFIG virtual inertia control on power system small-signal stability considering the phase-locked loop. IEEE Trans. Power Syst. 32(3), 2094–2105 (2017).

Xu, X., Liang, D., Wang, Z., Li, D., Zhang, A. & Xing, H. The effect of DFIG on electromechanical oscillation with combined reactive power damping control and virtual inertial control. In 2019 IEEE 8th International Conference on Advanced Power System Automation and Protection (APAP). 1284–1288 (2019).

Liu, C., Diao, S., Cai, G., Zhang, Y. & Yu, F. Construction of detailed energy function of double-fed induction generator for small signal stability analysis of power systems with a high penetration of double-fed induction generator. IET Gener. Transm. Distrib. 16, 4761–4774 (2022).

Wang, M. et al. Impedance modeling and stability analysis of DFIG wind farm with LCC-HVDC transmission. IEEE J. Emerg. Sel. Top. Circuits Syst. 12(1), 7–19 (2022).

Pulgar-Painemal, H., Wang, Y. & Silva-Saravia, H. On inertia distribution, inter-area oscillations and location of electronically-interfaced resources. IEEE Trans. Power Syst. 33(1), 995–1003 (2018).

Akhbari, A., Rahimi, M. & Khooban, M. H. Direct current grid-based doubly-fed induction generator wind turbines: Real-time control and stability analysis. IET Power Electron. 15, 1158–1173 (2022).

The authors would like to express their sincere gratitude to the experts on their unlimited support in the validation of the framework and suggestions as well as providing some data. The authors want to thank the editor and anonymous reviewers for their valuable comments and suggestions to improve the version of this article. This work was supported in part by the Natural Science Foundation of Xinjiang Autonomous Region under Grant (2020D01C068), in part by the National Natural Science Foundation of China under Grant (52167016).

College of Electrical Engineering, Xinjiang University, Urumqi, 830047, CO, People’s Republic of China

Liu Qi, Wu jiahui & Wang Haiyun

CGN New Energy Investment (Shenzhen) Co., Ltd., Xinjiang Branch, Xinjiang, 830011, China

Zhang Hua & Yang Jian

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

L.Q. and W.J. wrote the main manuscript text and W.H.，Z.H. and Y.J. prepared ALL figures and tables. All authors reviewed the manuscript.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Qi, L., jiahui, W., Haiyun, W. et al. Effect of DFIG control parameters on small signal stability in power systems. Sci Rep 13, 2476 (2023). https://doi.org/10.1038/s41598-023-29278-5

DOI: https://doi.org/10.1038/s41598-023-29278-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Windmill Generator For Home Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.