Journal of Modern Power Systems and Clean Energy

ISSN 2196-5625 CN 32-1884/TK

网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

Torsional Oscillation Damping Analysis and Suppression Strategy for PMSG-based Wind Generation System  PDF

  • Yizhuo Ma
  • Graduate (Student Member, IEEE)
  • Jin Xu
  • Guojie Li (Senior Member, IEEE)
  • Keyou Wang (Member, IEEE)
Key Laboratory of Control of Power Transmission and Conversion Ministry of Education, Shanghai Jiao Tong University, Shanghai 200240, China

Updated:2025-03-26

DOI:10.35833/MPCE.2024.000219

  • Full Text
  • Figs & Tabs
  • References
  • Authors
  • About
CITE
OUTLINE

Abstract

External disturbances can induce torsional oscillation with weak damping in the shaft system of permanent magnet synchronous generators (PMSGs) based wind generation system, thereby inducing low-frequency oscillations. However, the influence of electromagnetic torque on the shaft system damping and corresponding parameter laws have been scarcely explored. We define the electrical damping coefficient as a quantitative measure for the influence of electromagnetic torque on the shaft system damping. The torsional oscillation damping characteristics of the shaft system under vector control are analyzed, and the transfer function for electromagnetic torque and speed is derived. Additionally, we elucidate the mechanism by which the electromagnetic torque influences the shaft system damping. Simultaneously, laws describing the influence of wind speed, system parameters, and control parameters on the torsional oscillation damping are analyzed. Accordingly, the optimal damping angle of the shaft system a torsional oscillation suppression strategy is proposed to compensate for with uncertainty in the parameters affecting damping. The studied system is modeled using MATLAB/Simulink, and the simulation results validate the effectiveness of the theoretical analysis and proposed torsional oscillation suppression strategy.

I. Introduction

MANY countries have recently established carbon neutrality targets. In line with these targets, the installed wind power capacity is projected to grow by 430 GW from 2022 to 2025, reflecting an upward trend in wind power adoption [

1], [2]. The permanent magnet synchronous generator (PMSG) based wind generation system has gained popularity owing to benefits such as gearbox-free operation and high energy conversion efficiency and performance [3], [4]. The integration of large-scale grid-connected wind power systems poses various challenges to power systems, including power quality, voltage and frequency control, security, and stability [5]. Therefore, the dynamic characteristics of PMSG-based wind generation systems must be examined along with their interactions with power grids [6]. Furthermore, analyzing and understanding the dynamic characteristics of the shaft system in PMSG-based wind generation systems can notably enhance the safe and stable operation of power systems [7].

Compared with the doubly-fed induction generator, the shaft system of PMSG is simpler and consists of only three main components: wind turbine (WT), low-speed transmission shaft, and generator. These components are directly connected without a gearbox through a low-speed transmission [

8], [9]. Owing to its multipole structure, the shaft system of PMSG is more flexible than that in a conventional power plant [10]. In addition, the PMSG has more pole pairs in addition to the flexible shaft system [11].

Although a single-mass model may explain the transient instability in a PMSG under drastic disturbances, a double-mass model is required to accurately represent the system dynamics [

12]. This is because the double-mass model accurately represents the flexibility and mechanical oscillations of the shaft system.

Under severe external disturbances such as short-circuit faults, the shaft system in a PMSG-based wind generation system experiences torsional oscillation, resulting in power oscillation in the grid-connected wind power system [

11]. Torsional oscillations not only induce fatigue in the transmission shaft, thereby reducing its service life, but also give rise to low-frequency oscillations, potentially compromising the safety and stability of the entire PMSG-based wind generation systems [12]. Hence, the torsional oscillations in PMSG-based wind generation systems must be investigated [13].

The torsional oscillation characteristics of PMSG-based wind generation systems are influenced by three factors: input mechanical torque, inherent torque of the shaft system, and input electromagnetic torque under external disturbances [

14]. The former two factors are determined by the intrinsic system characteristics, including the generator rotor [15]. The effect of the electromagnetic torque on the torsional oscillation and its potential for oscillation suppression have been analyzed in [16] and [17]. However, a research gap persists concerning the specific parameters and underlying laws that govern the influence of the electromagnetic torque on torsional oscillation characteristics. Therefore, the damping characteristics of PMSG transmission shafting should be unveiled, and its influence mechanism should be understood. This study provides insights into the torsional oscillation characteristics of PMSG-based wind generation systems. Bridging this research gap may have practical implications for designing and optimizing torsional oscillation damping controllers.

To enhance the stability of PMSG-based wind generation systems, various active damping control methods have been devised to mitigate torsional oscillation [

18]-[20]. In [18], a simple real power control strategy is introduced based on the rapid torque control for a direct-drive wind energy conversion system employing a PMSG. This study also outlines the parameter tuning procedure for the proposed control strategy. In [19], a sensorless active damping control strategy for direct-driven permanent magnet WT generators uses cascade observers to estimate the speeds of the WT and PMSG. In [20], the torsional oscillation damping is improved by adding a damping transfer function proportional to the speed difference between the WT and generator. Although previous studies have suggested methods to suppress torsional oscillations, to the best of our knowledge, few of them have thoroughly and quantitatively analyzed the mechanism of the influence of electromagnetic torque on the shaft system damping. In addition, an optimal damping angle compensation method has not been devised to improve the system stability.

We analyze the torsional oscillation damping characteristics under vector control using an electromagnetic damping torque method. The electrical damping coefficient is defined as a quantitative representation of the effect of the electromagnetic torque on the shaft system damping. Then, the transfer function for electromagnetic torque and speed is derived, and the influence of the electromagnetic torque on the shaft system damping is characterized. Further, the impact of wind speed, system parameters, and control parameters on torsional oscillation damping are examined in detail. Based on these findings, we propose a torsional oscillation suppression strategy for the active power control section of a machine-side converter (MSC). This strategy compensates for the damping angle of the shaft system and ensures the maximum damping.

II. Modeling and Control of PMSG-based Wind Generation System

The topology of a grid-connected PMSG-based wind generation system typically includes the shaft system, PMSG, back-to-back full-power converter, transformers, and the control sections of the grid-side converter (GSC) and MSC, as shown in Fig. 1, where MPPT is short for maximum power point tracking; PI is short for proportional and integral; PCC is short for point of common coupling; and PLL is short for phase-locked loop. This section presents a mathematical model of the considered system using the per-unit system.

Fig. 1  Typical topology of PMSG-based wind generation system.

A. PMSG Model

The PMSG is controlled in the d-q rotating coordinates with the d-axis aligned with magnetic flux linkage of PMSG rotor ψf. The stator voltage of PMSG is given by [

21]:

usd=-Rsisd-Ldωebdisddt+ωgLqisqusq=-Rsisq-Ldωebdisqdt-ωgLdisd+ωgψf (1)

where usd and usq are the d- and q-axis stator terminal voltages, respectively; isd and isq are the d- and q-axis stator currents, respectively; Rs is the stator resistance of PMSG; ωeb is the base value of stator angular frequency; ωg is the angular frequency of PMSG rotor; and Ld and Lq are the d- and q-axis self-inductances of PMSG stator, respectively.

The speed of a megawatt PMSG is relatively low, and most PMSG is nonsalient surface-mounted (Ld=Lq). Therefore, the electromagnetic torque Te can be formulated as [

20]:

Te=ψfisq (2)

B. Double-mass Model

The double-mass model can be expressed as [

16]:

2Htdωtdt=Tm-Tsh2Hgdωgdt=Tsh-Tedθshdt=ωeb(ωt-ωg)Tsh=Kshθsh+Dsh(ωt-ωg) (3)

where Ht and Hg are the inertial time constants of WT and PMSG, respectively; ωt is the WT speed of generator rotor; θsh is the torsion angle of WT relative to the generator rotor; Ksh is the stiffness coefficient of shaft system; Dsh is the damping coefficient of shaft system; and Tm, Te, and Tsh are the mechanical, electromagnetic, and shaft system torques, respectively.

C. MSC Model

The active power generated by PMSG Pe is expressed as:

Pe=ψfisqωg (4)

From (4), Te and Pe can be accurately controlled by varying isq. The active power control loop of MSC is shown in Fig. 2.

Fig. 2  Active power control loop of MSC.

In Fig. 2, the superscript * represents the reference value; Kp1 and Ki1 are the proportional and integral coefficients for the power outer loop of MSC, respectively; and Kp2 and Ki2 are the proportional and integral coefficients of the current inner loop on MSC, respectively. Pe* can be expressed as:

Pe*=koptωg3    0<v<vr1              vrv<vmax (5)

where kopt is the MPPT curve coefficient; and v, vr, and vmax are the actual, rated, and maximum wind speeds, respectively.

Based on Fig. 2, the variable transfer of the active power control loop can be formulated as [

22]:

isq*=(Pe*-Pe)Kp1+Ki1susq=-(isq*-isq)Kp2+Ki2sisq=-usqLqs/ωeb+Rs (6)

III. Torsional Oscillation Damping Characteristics

This section presents the effect of electromagnetic torque Te on the torsional oscillation damping characteristics during MPPT. We derive the incremental transfer function of the electromagnetic torque and rotational speed difference, calculate the quantitative expression of electrical damping based on the electrical damping of synchronous generators (SGs), and analyze the effect of electrical damping on shaft system damping. Our findings provide a foundation for analyzing the influence of the PMSG parameters on the shaft system damping characteristics.

A. Electrical Damping Characteristics of PMSG

Using electromagnetic damping analysis, the linearized electromagnetic torque of PMSG ΔTe can be expressed as:

ΔTe=DeΔωΔ+KeΔθsh (7)

where De and Ke are the electrical damping and synchronization coefficients of the PMSG, respectively; and ΔωΔ is the difference between Δωg and Δωt, i.e., ΔωΔ=Δωt-Δωg.

We omit self-damping and the generator rotor. Using (3) and (7), the transfer function for torsional angle θsh and mechanical torque Tm is given by:

ΔθshΔTm=ωeb2Hts2+Dsh2Ht+Dsh2Hg-De2Hgs+ωebKsh2Ht+Ksh2Hg-Ke2Hg (8)

The damping attenuation factor for torsional oscillation ξ can be determined using (8) as:

ξ=Dsh(Ht+Hg)-DeHt4HtHgωosc (9)

where ωosc is the natural oscillation angular frequency of the shaft system, and ωosc=ωeb[Ksh/(2Ht)+(Ksh-Ke)/(2Hg)].

The electrical damping coefficient De has a negative correlation with the torsional oscillation damping coefficient ξ. When De<0, the phase difference between ΔTe and ΔωΔ lies in (90°,270°), resulting in an increase in ξ. This in turn leads to a positive damping of the shaft system, enhancing its stability. When De>0, the phase difference between ΔTe and ΔωΔ falls in (0,90°), causing ξ to decrease. Consequently, the shaft system experiences negative damping and thus instability.

B. Shaft System Damping Characteristics of PMSG

To obtain the small-signal output power of the PMSG, the MPPT curve can be approximated by a linear function around the steady-state operating point as:

ΔPe*=3koptωg02Δωg (10)

where ΔPe* is the linearized increment of the reference generator output power; Δωg is the linearized increment of the generator speed; and ωg0 is the generator speed at the stable running point.

The linearized increment of reference electromagnetic torque on the MPPT curve can be formulated as:

ΔTe*=2koptωg0Δωg (11)

By combining (2) and (4)-(6), the electromagnetic torque can be expressed as:

Te=ψfisq=ψfkoptn1n2ωg3-n1n2TeωgRs+Lqs/ωeb+n2 (12)

where n1=Kp1+Ki1/s; and n2=Kp2+Ki2/s.

Linearizing both sides of the equation simultaneously allows to formulate the incremental transfer function for the linearized electromagnetic torque ΔTe and ΔωΔ as (13).

ΔTeΔωΔ=-ψfHtHt+Hg2koptωg02[Kp1Kp2s2+Lqs3/ωeb+(Rs+Kp2+ωg0ψfKp1Kp2)s2+(Kp1Ki2+Ki1Kp2)s+Ki1Ki2][Ki2+ωg0ψf(Kp1Ki2+Ki1Kp2)]s+ωg0ψfKi1Ki2 (13)

Further, considering (13), we have (14).

m=2koptωg02HtHt+HgG(s)=-ψfKp1Kp2s2+(Kp1Ki2+Ki1Kp2)s+Lqs3/ωeb+(Rs+Kp2+ωg0ψfKp1Kp2)s2+     Ki1Ki2[Ki2+ωg0ψf(Kp1Ki2+Ki1Kp2)]s+ωg0ψfKi1Ki2 (14)

The derivation process is detailed in Supplementary Material A.

The linearized electromagnetic torque can be expressed as ΔTe=mG(s)ΔωΔ.

According to (7) and (13), the real part of the incremental transfer function for ΔTe and ΔωΔ corresponds to the electrical damping coefficient De. Moreover, as indicated by (9), De is the sole controllable parameter influencing shaft system damping. Hence, investigating the impact of system parameters on De is essential for enhancing the damping.

IV. Influence of System Parameters on Torsional Oscillation Damping Characteristics

In this section, we examine the correlation between De and various system parameters: ωg0, Rs, Kp1, Ki1, Kp2, Ki2, and v. Specifically, we investigate the effects of wind speed v, the PMSG parameters, and control parameters on De and the oscillation characteristics of shaft system. Accordingly, we establish a theoretical basis for the proposed torsional oscillation suppression strategy outlined in Section V.

The main parameters of the PMSG-based wind generation system provided by China Wind Power Group Limited are listed in Table I. In addition, the control parameters are tuned using Simulink to achieve the optimal performance.

TABLE Ⅰ  Main Parameters of PMSG-based Wind Generation System
ParameterValue
Base power value Sb 2 MW
Base value of AC phase voltage Ub 575 V
Base value of DC-link voltage Udcb 1150 V
Base value of stator angular frequency ωeb 377 rad/s
Base value of PMSG rotor speed ωmb 7.85 rad/s
Base value of grid angular frequency ωb 377 rad/s
Rated frequency of PMSG fn 60 Hz
Pole pair of PMSG np 48
Magnetic flux linkage of PMSG rotor ψf 1.18842 p.u.
d-axis self-inductance of PMSG stator Ld 0.5131 p.u.
q-axis self-inductance of PMSG stator Lq 0.5131 p.u.
Stator resistance of PMSG Rs 0.0001 p.u.
Inertia time constant of WT Ht 6.69 s
Inertia time constant of PMSG Hg 1 s
Damping coefficient of shaft system Dsh 1 p.u.
Stiffness coefficient of shaft system Ksh 1.6 p.u.
Proportional coefficient of power outer loop on MSC Kp1 1 p.u.
Integral coefficient of power outer loop on MSC Ki1 20 s-1
Proportional coefficient of current inner loop on MSC Kp2 1 p.u.
Integral coefficient of current inner loop on MSC Ki2 10 s-1
Capacitance of DC-link voltage C 0.6232 p.u.
Proportional coefficient of voltage outer loop on GSC Kp3 5 p.u.
Integral coefficient of voltage outer loop on GSC Ki3 300 s-1
Proportional coefficient of current inner loop on GSC Kp4 0.83 p.u.
Integral coefficient of current inner loop on GSC Ki4 5 s-1
Grid inductance Lcg 0.55 p.u.
Grid resistance Rcg 0.006 p.u.
Lowpass filter for Pe ωP 100 rad/s
PLL filter ωf 100 rad/s
Proportional coefficient of PLL Kp,pll 4.1 p.u
Integral coefficient of PLL Ki,pll 200 s-1

A. Influence of Wind Speed and Stator Resistance on De

Figure 3 shows the Bode diagram of ΔTe/ΔωΔ with various ωg0 ranging from 0.6 to 1.0 p.u., where θosc is the phase angle of transfer function G(s) at ωsoc. The phase lag of ΔTe with respect to ΔωΔ falls between 90° and 180° near ωosc.

Fig. 3  Bode diagram of ΔTe/ΔωΔ with various ωg0.

Hence, De is negative, resulting in positive damping provided by ΔTe for the shaft system, thus improving the system stability. The De values corresponding to various ωg0 are presented in Table II. De decreases as ωg0 increases, indicating an enhanced effect of the electromagnetic torque on the shaft system damping.

TABLE Ⅱ  De Values Corresponding to Various ωg0
ωg0 (p.u.)De
0.6 -0.248
0.7 -0.334
0.8 -0.428
0.9 -0.527
1.0 -0.626

The variation range of the stator resistance should be moderated when investigating the impact of temperature changes on Rs. The Bode diagram of ΔTe/ΔωΔ and the corresponding De are shown in Fig. 4 and Table III, respectively, with various Rs. The curves in the diagram are very similar, suggesting that changes in Rs have a small influence on De. Thus, the effect of ΔTe on the shaft system damping is negligible and does not affect the stability of the shaft system under MPPT.

Fig. 4  Bode diagram of ΔTe/ΔωΔ  with various Rs.

TABLE Ⅲ  De Values Corresponding to Various Rs
Rs (p.u.)De (p.u.)
0.0001 -0.6260
0.0003 -0.6260
0.0005 -0.6260
0.0007 -0.6260
0.0009 -0.6259

B. Influence of Control System Parameters on De

The active power control loop comprises two PI controllers and four control parameters: Kp1, Ki1, Kp2, and Ki2. We set Kp1=1 p.u., Ki1=20 s-1, Kp2=1 p.u., and Ki2 =10 s-1 as default control parameters. Subsequently, we select parameter within the ranges of 1-3 p.u. for Kp1, 5-25 s-1 for Ki1, 1-3 p.u. for Kp2, and 10-30 s-1 for Ki2.

Using (13), the Bode diagram of ΔTe/ΔωΔ with various Kp1 ranging from 1 to 3 p.u. is shown in Fig. 5. The phase lag of ΔTe with respect to ΔωΔ falls between 90° and 180° near ωosc for Kp1 of 1.0, 1.5, 2.0, 2.5, and 3.0 p.u.. In such cases, De is negative, indicating that Te contributes to positive damping of the shaft system, thus enhancing its stability. When Kp1=3.0 p.u., the phase lag of Te with respect to ΔωΔ is between 0° and 90° near ωosc. In this case, the negative value of De suggests that Te still contributes to positive damping of the shaft system, but its stability deteriorates. Table IV lists the De values corresponding to various Kp1. As the PMSG operates under the MPPT, De transitions from negative to positive as Kp1 increases. This indicates a weakening effect of the electromagnetic torque on the shaft system damping, which is unfavorable to stability.

Fig. 5  Bode diagram of ΔTe/ΔωΔ with various Kp1.

TABLE Ⅳ  De Values Corresponding to Various Kp1
Kp1 (p.u.)De  (p.u.)
1.0 -0.63
1.5 -0.51
2.0 -0.30
2.5 0.01
3.0 0.31

We also use (13) to obtain the Bode diagram of ΔTe/ΔωΔ with various Ki1 ranging from 5 to 25 s-1 in steady states. As shown in Fig. 6, for Ki1 of 5, 10, 15, 20, and 25 s-1, the phase lag of ΔTe with respect to ΔωΔ falls between 90° and 180° near ωosc. In such cases, De is negative, indicating that Te contributes to the positive damping of shaft systems, thereby enhancing its stability. The De values corresponding to various Ki1 are listed in Table V. When the PMSG operates under the MPPT, De  decreases as Ki1 increases, thereby augmenting the damping effect of the electromagnetic torque on the shaft system and improving its stability.

Fig. 6  Bode diagram of ΔTe/ΔωΔ with Ki1.

TABLE V  De Values Corresponding to Various Ki1
Ki1 (s-1)De (p.u.)
5 -0.36
10 -0.63
15 -0.96
20 -1.37
25 -1.91

The analysis regarding Kp2 and Ki2 and their impacts on the shaft stability is analogous to the analysis conducted for Kp1 and Ki1, respectively, as given in Figs. 7 and 8 and Tables VI and VII, which is thus not discussed further in this subsection.

Fig. 7  Bode diagram of ΔTe /ΔωΔ  with Kp2 .

Fig. 8  Bode diagram of ΔTe/ΔωΔ with various Ki2 .

TABLE VI  De Values Corresponding to Various Kp2
Kp2 (p.u.)De  (p.u.)
1.0 -0.63
1.5 -0.41
2.0 -0.15
2.5 0.10
3.0 0.33
TABLE VII  De Values Corresponding to Various Ki2 
Ki2 (s-1)De  (p.u.)
10 -0.41
15 -0.63
20 -0.88
25 -1.20
30 -1.61

V. Torsional Oscillation Suppression Strategy

We present a system comprising MPPT, MSC, and GSC control strategies. The MSC aims to regulate the output active power Pe and the d-axis stator current.

An additional damping controller is introduced into the active power control loop to mitigate torsional oscillations in a PMSG-based wind generation system.

A. Proposed Torsional Oscillation Suppression Strategy

To enhance the effect of the electromagnetic torque on the shaft system damping, De should be reduced, as discussed in Section IV. Therefore, we introduce an additional damping controller into the active power control loop. Figure 9 shows the linearized model of the transfer relationship between ΔTm and ΔωΔ, which is obtained using (1) and (13). The dotted line in Fig. 9 represents the addition of damping controller H(s). The transitive relation between ΔTe and ΔωΔ can be derived as:

ΔTe=(m+H(s))G(s)ΔωΔ=(m+H(jω))G(jω)ΔωΔ=mG(jω)ejθΔωΔ+H(jω)G(jω)ejθ+φΔωΔ (15)

Fig. 9  Linearized model of transfer relationship between  ΔTm and ΔωΔ.

where θ is the phase angle of transfer function G(s); and φ is the phase angle of transfer function H(s).

With the additional damping controller in the active power control loop, De near ωosc can be expressed as:

De=De0+DeHDe0=mRe(|G(jωosc)|ejθosc)DeH=Re(|H(jωosc)||G(jωosc)|ejφosc) (16)

where De0 is the electrical damping coefficient without the additional damping controller; DeH is the additional electrical damping coefficient with the additional damping controller; and φosc is the phase angle of transfer function H(s) at ωosc.

From (16) and the expression of G(s), the system reaches stability, with both the amplitude and phase of G(s) remaining stable once the operating conditions, control parameters, and system parameters are determined. Therefore, based on G(s), we can adjust the amplitude and phase of |H(jωosc)| to compensate for phase φosc by modifying additional damping controller H(s) and the internal parameters. This adjustment aims to decrease De, enhance the damping effect of Te on the shaft system, and improve the stability of the shaft system.

The additional damping controller H(s) comprises two primary components: bandpass filter controller Hf(s) and phase compensation controller Hi(s). The bandpass filter controller Hf(s) employs a second-order bandpass filter expressed as:

Hf(s)=2ξfωnss2+2ξfωns+ωn2 (17)

where ξf is the damping ratio of Hf(s); and ωn is the center angular frequency of Hf(s).

We set the damping ratio of Hf(s) to be ξf=0.15. In addition, the center angular frequency of the bandpass filter, ωn=12.748 rad/s, corresponds to the inherent oscillation angular frequency of the shaft system.

The phase compensation controller Hi(s) enhances the controller performance, and its transfer function is given by:

Hi(s)=Ktod1+sT11+sT22 (18)

where Ktod is the damping controller gain; and T1 and T2 are the lead and lag correction time constants, respectively.

The additional damping controller H(s) is formulated as:

H(s)=Hf(s)Hi(s)=Ktod1+sT11+sT222ξfωnss2+2ξfωns+ωn2 (19)

Figure 10 shows a model of a PMSG-based wind farm (WF) connected to a four-machine two-area system. In this model, the WF comprises two PMSG-based wind generation systems with different parameters, which are equivalent to two 350 MW WFs, denoted as WF1 and WF2, respectively [

23]. The parameters for the AC grid and SGs are shown in Supplementary Material B and have some modifications. The parameters for the PMSG-based wind generation system of WF1 are listed in Table I, whereas those for WF2 are provided in Supplementary Material C.

Fig. 10  Model of PMSG-based WF connected to four-machine two-area system.

This model is employed to evaluate the torsional oscillation damping characteristics of PMSG-based wind generation systems with various parameters. We also investigate the independence of the torsional oscillation suppression strategies when multiple PMSG-based wind generation systems with different parameters are connected to power grid. Additionally, the coupling characteristics of the damping controller and torsional oscillation are examined.

Increasing Ktod enhances the torsional oscillation suppression ability of the damping controller but introduces a new oscillation mode. Excessive Ktod values can lead to positive real roots, thereby impacting the system stability. Therefore, the impact of Ktod on the system stability should be analyzed. Supplementary Material D presents a detailed linearization model of a PMSG-based WF connected to the four-machine two-area system.

The main oscillations are categorized into eight modes (modes 1-8): ① local oscillation mode dominated by the SGs in area 1; ② local oscillation mode dominated by the SGs in area 2; ③ torsional oscillation mode dominated by shaft system of WF1; ④ oscillation mode dominated by the PLL; ⑤ interarea oscillation mode dominated by all the SGs; ⑥ oscillation mode dominated by the proposed additional damping controller in WF1; ⑦ torsional oscillation mode dominated by shaft system of WF2; and ⑧ oscillation mode dominated by the proposed additional damping controller in WF2.

First, we vary Ktod1 from 0 to 24 with an increment step of 0.08. The root loci of oscillation modes are shown in Fig. 11(a). As shown in Fig. 11(b), when Ktod1 increases, modes 3 and 6 initially converge and then gradually diverge. The absolute value for the damping ratio of mode 3 first increases and then decreases, and when Ktod1 reaches 0.48, it reaches the maximum value. Additionally, the variation in Ktod1 has no effect on modes 7 and 8. Concurrently, we increase Ktod2 from 0 to 24 with an increment step of 0.08. As shown in Fig. 12, with an increase in Ktod2, the root loci of modes 7 and 8 follow the same trend as those of modes 3 and 6 when Ktod1 increases, respectively. Furthermore, the changes in the damping characteristics of shaft system for WF2 do not affect those for WF1. Therefore, when dealing with multiple WTs connected to AC power systems, a shaft system damping controller should be designed considering the characteristics of each PMSG-based wind generation system.

Fig. 11  Root loci and damping ratio of oscillation modes 1-8 with various Ktod1. (a) Root loci. (b) Damping ratio.

Fig. 12  Root loci and damping ratio of oscillation modes 1-8 with various Ktod2. (a) Root loci. (b) Damping ratio.

B. Implementation of Proposed Torsional Oscillation Suppression Strategy

The proposed oscillation suppression strategy is summarized in the following steps.

Step 1:   determine the system operating and control parameters such as wind speed, inertial time constants, stator resistance, Kp1, Ki1, Kp2, and Ki2.

Step 2:   calculate the torsional oscillation frequency.

Step 3:   obtain the Bode diagram of ΔTe/ΔωΔ using (13).

Step 4:   calculate the phase angle of the Bode diagram at the torsional oscillation frequency and obtain the compensation angle.

Step 5:   calculate the parameters of the phase compensation and bandpass filter controllers based on the compensation angle.

Step 6:   obtain the optimal damper controller gain Ktod based on eigenvalue analysis.

Based on the parameters listed in Table I, when the system reaches steady state, the phase-frequency characteristics of G(s) are depicted in Fig. 13 with various ωg0 ranging from 0.7 to 1.1 p.u.. Figure 13 shows that θosc ranges from -132° to -122° at ωosc of 12.748 rad/s. For instance, taking ωg0=1.1 p.u. as a reference, the phase angle at ωosc in the shaft system stabilizes at approximately -122° during steady-state operation.

Fig. 13  Phase-frequency characteristic of G(s) with various ωg0.

Based on (9) and (13), when ΔTe and ΔωΔ are in antiphase, De reaches its minimum value, and the damping provided by ΔTe is maximized, considerably benefiting the shafting stability. Hence, H(s) is used to adjust the phase angle of G(s) at ωosc to -180°. The phase angle that Hi(s) must compensate is φosc=-130°-θosc=-53°. Time constants T1 and T2 are computed as 0.56 and 0.2 s, respectively. The phase-frequency characteristic of G(s) after incorporating H(s) is shown in Fig. 14. By comparing Figs. 13 and 14, the phase angle of G(s) at ωosc has been adjusted to -180°, maximizing the damping effect of ΔTe.

Fig. 14  Phase-frequency characteristic of G(s) with various ωg0 after incorporating H(s).

VI. Case Study

We conduct simulations using the MATLAB/Simulink platform to evaluate the effectiveness of the proposed torsional oscillation suppression strategy. The simulations are executed on a personal computer equipped with a 3.3 GHz processor and 16 GB of RAM to implement the model depicted in Fig. 1. The converter is modeled as a digital system with a switching frequency of 2 kHz and sampling time of 5×10-6 s.

The parameters of the PMSG-based wind generation system listed in Table I are obtained from China Wind Power Group Limited based on an actual system. In addition, the control parameters are fine-tuned using Simulink to achieve optimal performance.

A. Evaluation of Proposed Torsional Oscillation Suppression Strategy

Throughout the simulation, the wind speed is maintained constant at 12.1 m/s, which corresponds to the rated wind speed. The time-domain simulation incorporates a three-phase short-circuit fault on node B12 of the studied system at t=50 s. The fault persists for 0.2 s before being cleared at t=50.2 s, after which the system returns to steady state.

The time-domain simulation is aimed to investigate two scenarios with and without applying the proposed torsional oscillation suppression strategy, i.e., scenarios 1 and 2, respectively.

The simulation results of studied system in scenarios 1 and 2 are shown in Fig. 15. The application of the double-mass model and proposed torsional oscillation suppression strategy results in a stabilization of ωg  after the short-circuit fault is cleared. In addition, the system attains transient stability, with Pg, Udc, and voltage at PCC Uo returning to steady-state levels, accompanied by oscillations following the disturbance. Implementing the proposed torsional oscillation suppression strategy mitigates transient instability and reduces the likelihood of low-frequency oscillations during the short-circuit fault. Consequently, the oscillations in ωg, Udc, Uo, and Pg are attenuated, enhancing the dynamic and steady-state operation of the studied system and improving the power grid stability compared with the results without the additional damping controller.

Fig. 15  Simulation results of studied system in scenarios 1 and 2. (a) ωg. (b) Pg. (c) Udc. (d) Uo.

Then, the fundamental parameters affecting the damping of the shaft system, including wind speed, system parameters, and control parameters, are analyzed. The simulation curves of ωg and Pg with different wind speeds under MPPT are shown in Fig. 16. Figure 17 shows the simulation curves of ωg and Pg with various Rs under MPPT. Additionally, Fig. 18 illustrates the simulation curves of ωg and Pg with various Kp1 and Ki1 under MPPT. The results for Kp2 and Ki2 are ignored because they are similar to those for Kp1 and Ki1, respectively.

Fig. 16  Simulation curves of ωg  and Pg  with different wind speeds under MPPT. (a) ωg with v=1.0 p.u.. (b) ωg with v=0.8 p.u.. (c) ωg with v=0.6 p.u.. (d) Pg with v=0.8 p.u.. (e) Pg  with v=1.0 p.u.. (f) Pg  with v=0.6 p.u..

Fig. 17  Simulation curves of ωg  and Pg with various Rs under MPPT. (a) ωg. (b) Pg.

Fig. 18  Simulation curves of ωg and Pg with various Kp1 and Ki1. (a) ωg with various Kp1. (b) Pg with various Kp1. (c) ωg with various Ki1. (d) Pg with various Ki1.

Figures 16-18 show that the oscillation amplitudes of ωg  and Pg are smaller and the oscillation attenuation is faster with increase of Ki1 and Ki2, indicating the shaft system damping improves. However, the oscillation amplitude of ωg  and Pg is larger and the oscillation attenuation is slower with increase of Kp1 and Kp2, indicating that the shaft system damping deteriorates. Additionally, the simulation results are nearly insensitive to variations in stator resistance Rs within a reasonable range influenced by temperature.

B. Comparison of Proposed Torsional Oscillation Suppression Strategy with Other Strategies

To evaluate the performance of the proposed torsional oscillation suppression strategy, we compared it with common strategies, including speed-feedback-based damping and active damping based on torque estimation [

19].

We consider two cases that illustrate the effectiveness of the proposed torsional oscillation suppression strategy. In case 1, wind speed v remains constant at 12.1 m/s, as in Section VI-A. A three-phase short-circuit fault occurs at node B12 of the system and t=50 s in the time-domain simulation. The fault persists for 0.2 s before being cleared at t=50.2 s, after which the system returns to steady state. In case 2, the initial wind speed is 12.1 m/s and abruptly drops to 7.3 m/s at t=50 s.

Figure 19 shows the transient response curves of ωg, Pg, Udc, Uo, Te, and θsh in response to the three-phase short-circuit fault using the proposed torsional oscillation suppression strategy, active damping strategy considering torque estimation, and damping strategy based on speed feedback [

3], [19].

Fig. 19  Transient response curves of ωg, Pg, Udc, Uo, Te, and θsh in response to three-phase short-circuit fault using different control strategies. (a) ωg. (b) Pg. (c) Udc. (d) Uo. (e) Te. (f) θsh.

The proposed torsional oscillation suppression strategy and other strategies improve the shaft system damping and suppress the oscillatory instability of the PMSG-based wind generation system during the three-phase short-circuit fault. The proposed torsional oscillation suppression strategy is more effective in suppressing torsional oscillations than the comparison strategies. Additionally, the proposed torsional oscillation suppression strategy attenuates the oscillations of ωg, Pg, Udc, Te, and θsh faster than other strategies. These results indicate that the proposed torsional oscillation suppression strategy outperforms other strategies in suppressing torsional oscillations.

Figure 20 shows that following a drop in wind speed, Pg, Te, and ωg initially decrease, and then stabilize after oscillation attenuation.

Fig. 20  Transient response curves of ωg, Pg, Udc, Uo, Te, and θsh following a drop in wind speed. (a) ωg. (b) Pg. (c) Udc. (d) Uo. (e) Te. (f) θsh.

Likewise, Udc, θsh, and Uo oscillate before reaching steady state. The proposed torsional oscillation suppression strategy demonstrates a superior performance in attenuating the oscillations of ωg, θsh, and Pg more rapidly compared with other strategies. This superiority is evident in Fig. 20(f), where the oscillation of θsh is promptly suppressed. The simulation results validate the superiority of the proposed torsional oscillation suppression strategy in suppressing torsional oscillations.

C. Verification of Effect of Ktod on System Stability

In this subsection, we analyze the influence of Ktod on the system stability. Figure 21 shows the transient response curves of ωg, Pg, Udc, Te, and θsh with various Ktod. Figure 21(a) and (e) shows that ωg and θsh decrease while the shaft system damping increases with increasing Ktod, reaching its peak when Ktod=0.48. When Ktod increases to 2, the shaft system damping decreases, and the decay rates of ωg and θsh decline.

Fig. 21  Transient response curves of ωg, Pg, Udc, Te, and θsh with various Ktod. (a) ωg. (b) Pg. (c) Udc. (d) Te. (e) θsh.

In Fig. 21(b)-(d), an increase in Ktod results in an increase in the active power reference value, leading to higher oscillation amplitudes in Pg, Udc, and Te. Therefore, the proper selection of Ktod is crucial for effective oscillation suppression. If Ktod is excessively small, the suppression effect becomes prominent. By contrast, if Ktod is very large, the suppression effect diminishes, increasing the instantaneous values of Pg, Udc, and Te, and potentially causing system instability.

VII. Conclusion

We analyze the torsional oscillation damping characteristics based on the electromagnetic damping torque method. Subsequently, the transfer function for electromagnetic torque and speed is derived, and the influence of electromagnetic torque on the shaft system damping is explained.

The electrical damping coefficient De is the only controllable parameter that affects the shaft system damping and is negatively correlated with damping coefficient ξ of the torsional oscillation. Additionally, the real part of the incremental transfer function of ΔTe  and ΔωΔ is equivalent to electrical damping coefficient De.

The influence of wind speed and system control parameters on torsional oscillation damping is determined. We find that v, Ki1, and Ki2 are positively correlated with the shaft system damping, while Kp1 and Kp2 are negatively correlated with shaft system damping, and a change in stator resistance Rs within a reasonable range has a negligible effect on damping. Therefore, to ensure the stable operation of a PMSG-based wind generation system, the high-speed operation is recommended. In addition, while ensuring an active power tracking rate and stability in the high-frequency range, we recommend minimizing Kp1 and Kp2, as well as increasing Ki1 and Ki2.

We propose an oscillation suppression strategy on the active power control loop to compensate for the damping angle between ΔTe and ΔωΔ and thus improve the stability of the shaft system based on the previous analysis. The studied system is modeled using MATLAB/Simulink, and the simulation results show that the proposed torsional oscillation suppression strategy outperforms other strategies in eliminating torsional oscillations. The simulation results verifies the effectiveness of our theoretical analysis and control strategy.

In future work, we will investigate an adaptive compensation method for the damping angle of a PMSG-based wind generation system considering the frequency offset.

References

1

H. Liu and C. Liu, “Frequency regulation of VSC-MTDC system with offshore wind farms,” Journal of Modern Power Systems and Clean Energy, vol. 12, no. 1, pp. 275-286, Jan. 2024. [Baidu Scholar] 

2

B. Sahu and B. P. Padhy, “Design of power system stabilizer for DFIG-based wind energy integrated power systems under combined influence of PLL and virtual inertia controller,” Journal of Modern Power Systems and Clean Energy, vol. 12, no. 2, pp. 524-534, Mar. 2024. [Baidu Scholar] 

3

A. Beiki and M. Rahimi, “An efficient sensorless approach for energy conversion enhancement and damping response improvement in permanent magnet synchronous generator (PMSG) based wind turbines,” International Transactions on Electrical Energy Systems, vol. 29, no. 1, pp. 18-35, Jan. 2019. [Baidu Scholar] 

4

K. A. Singh, A. Chaudhary, and K. Chaudhary, “Three-phase AC-DC converter for direct-drive PMSG-based wind energy conversion system,” Journal of Modern Power Systems and Clean Energy, vol. 11, no. 2, pp. 589-598, Mar. 2023. [Baidu Scholar] 

5

J. Tan and Y. Zhang, “Coordinated control strategy of a battery energy storage system to support a wind power plant providing multi-timescale frequency ancillary services,” IEEE Transactions on Sustainable Energy, vol. 8, no. 3, pp. 1140-1153, Jul. 2017. [Baidu Scholar] 

6

C. Zhang, Z. Li, Q. Gao et al., “Damping effects on torsional oscillation of DFIG drive-chain using different control strategies,” Proceedings of the CSEE, vol. 33, no. 27, pp. 135-144, Sept. 2013. [Baidu Scholar] 

7

M. Nasiri and R. Mohammadi, “Peak current limitation for grid side inverter by limited active power in PMSG-based wind turbines during different grid faults,” IEEE Transactions on Sustainable Energy, vol. 8, no. 1, pp. 3-12, Jan. 2017. [Baidu Scholar] 

8

H. Geng, D.-W. Xu, B. Wu et al., “Control and Stability analysis for the permanent magnetic synchronous generator based direct driven variable speed wind energy conversion system,” Proceedings of the CSEE, vol. 29, no. 33, pp. 68-75, Nov. 2009. [Baidu Scholar] 

9

G. Wang, Q. Shi, L. Fu et al., “Mechanism analysis of torsional vibration for directly-driven wind turbine with permanent magnet synchronous generator shaft system with virtual inertia control,” Electric Machines and Control, vol. 18, no. 8, pp. 8-16, Aug. 2014. [Baidu Scholar] 

10

L. Chen, X. Du, B. Hu et al., “Drivetrain oscillation analysis of grid forming type-IV wind turbine,” IEEE Transactions on Energy Conversion, vol. 37, no. 4, pp. 2321-2337, Dec. 2022. [Baidu Scholar] 

11

Z. Hao and Y. Yu, “Torque analysis for the influence of DFIG-based wind farm on damping of power system,” Transactions of China Electrotechnical Society, vol. 26, no. 5, pp. 152-158, Oct. 2011. [Baidu Scholar] 

12

A. J. G. Westlake, J. R. Bumby, and E. Spooner, “Damping the power-angle oscillations of a permanent-magnet synchronous generator with particular reference to wind turbine applications,” IEE Proceedings-Electric Power Applications, vol. 143, no. 3, pp. 269-280, May 1996. [Baidu Scholar] 

13

A. H. Rajaei, M. Mohamadian, S. M. Dehghan et al., “PMSG-based variable speed wind energy conversion system using Vienna rectifier,” European Transactions on Electrical Power, vol. 21, no. 1, pp. 954-972, Jan. 2011. [Baidu Scholar] 

14

J. Yao, X. Zeng, and J. Li, “Shaft oscillation characteristics of grid-connected doubly-fed induction generator-based wind power generation system,” Transactions of China Electrotechnical Society, vol. 32, no. 6, pp. 123-135, Jun. 2017. [Baidu Scholar] 

15

L. Wang, L. Cheng, Y. Sun et al., “Damping control of drive-train system of DFIG to compensate phase lag characteristics of electromagnetic torque-generator speed closed loop,” Power System Technology, vol. 38, no. 12, pp. 3333-3340, Mar. 2014. [Baidu Scholar] 

16

Z. Zhang, X. Zhao, L. Fu et al., “Stability and dynamic analysis of the pmsg-based WECS with torsional oscillation and power oscillation damping capabilities,” IEEE Transactions on Sustainable Energy, vol. 13, no. 4, pp. 2196-2210, Oct. 2022. [Baidu Scholar] 

17

S. S. H. Yazdi, Y. Shokri-Kalandaragh, and M. Bagheri, “Power system stability improvement considering drive train oscillations of virtual synchronous generator-regulated type-4 wind turbines,” IET Renewable Power Generation, vol. 17, no. 3, pp. 579-603, Feb. 2023. [Baidu Scholar] 

18

O. Alizadeh and A. Yazdani, “A strategy for real power control in a direct-drive PMSG-based wind energy conversion system,” IEEE Transactions on Power Delivery, vol. 28, no. 3, pp. 1297-1305, Jul. 2013. [Baidu Scholar] 

19

M. Rahimi and A. Beiki, “Efficient modification of the control system in PMSG-based wind turbine for improvement of the wind turbine dynamic response and suppression of torsional oscillations,” International Transactions on Electrical Energy Systems, vol. 28, no. 8, pp. 16, Aug. 2018. [Baidu Scholar] 

20

O. Alizadeh and A. Yazdani, “A control strategy for power regulation in a direct-drive WECS with flexible drive-train,” IEEE Transactions on Sustainable Energy, vol. 5, no. 4, pp. 1156-1165, Oct. 2014. [Baidu Scholar] 

21

W. Du, Y. Wang, Y. Wang et al., “Analytical examination of oscillatory stability of a grid-connected PMSG wind farm based on the block diagram model,” IEEE Transactions on Power Systems, vol. 36, no. 6, pp. 5670-5683, Nov. 2021. [Baidu Scholar] 

22

F. Zhou, “Drive-train torsional vibration suppression of large scale PMSG-based WECS,” Protection and Control of Modern Power Systems, vol. 7, no. 1, pp. 1-13, Dec. 2022. [Baidu Scholar] 

23

B. Shao, S. Zhao, B. Gao et al., “Adequacy of the single-generator equivalent model for stability analysis in wind farms with VSC-HVDC systems,” IEEE Transactions on Energy Conversion, vol. 36, no. 2, pp. 907-918, Jun. 2021. [Baidu Scholar]