Abstract

The virtual impedance can suppress the low-frequency oscillation effectively caused by the power coupling effect for the grid-forming direct-drive permanent magnet synchronous generator (D-PMSG). However, there are few studies on the influence of virtual impedance on the subsynchronous oscillation caused by the interaction between the grid-forming D-PMSG and the weak grid. In this paper, the mechanism of virtual impedance improving the stability of virtual synchronous generator (VSG) control for the converter is revealed from the perspective of power coupling, and the influence of virtual impedance on the system stability in subsynchronous frequency range was analysed based on the sequence impedance model. The analysis results show that the virtual inductance cannot effectively improve the oscillation suppression effect, but the virtual resistance can greatly improve the damping characteristics of the system which has a good oscillation suppression effect. However, excessive increase in virtual resistance will cause power coupling oscillation, even subsynchronous oscillation, which shows limitations in subsynchronous oscillation suppression. To solve the problem, a dynamic voltage compensation control method based on active disturbance rejection control (ADRC) was proposed, and the parameters are designed by bandwidth method. The method can suppress the subsynchronous oscillation effectively in a certain frequency range under different grid strength and disturbance conditions. And it shows strong adaptability and good system robustness in weak grid. Finally, the effectiveness of the proposed oscillation suppression strategy was verified by time domain simulation.

1. Introduction

Under the background of continuous growth of global wind power installed capacity, to make up for the lack of system inertia caused by the introduction of high proportion power electronic equipment, the grid-forming VSG control technology is applied to the grid-side converter of D-PMSG, which provides inertia support to the power system by mimicking the transient characteristics of synchronous generator [1–3]. Under the condition of high proportion power electronic equipment, the interaction between D-PMSG and the weak grid is prone to subsynchronous oscillation, which is particularly prominent in remote areas [4, 5]. The virtual impedance can suppress the low-frequency oscillation effectively caused by the power coupling effect for the grid-forming D-PMSG. However, there are few studies on the influence of virtual impedance on the subsynchronous oscillation caused by the interaction between the grid-forming D-PMSG and the weak grid [6, 7]. The paper will focus on this issue to carry out research.

The sequence impedance analysis method is widely used in the stability analysis of new energy grid connection because of its easy engineering measurement and simple form [8–10]. The mechanism and characteristics of the subsynchronous oscillation in weak grid were analysed according to distribution characteristics of the impedance model for D-PMSG in the frequency domain [11, 12]. The equivalent capacitive characteristics of D-PMSG interact with the inductive characteristics of the weak grid and generate electrical resonance. Then, the subsynchronous oscillation occurs with resonance divergence under the equivalent negative damping impact of phase-locked loop [13]. To avoid the negative damping impact and provide inertia support to the power system, the grid-forming VSG control was proposed to apply in the grid-side converter of D-PMSG [14, 15].

The rotor motion state equation is introduced in VSG control which is based on the droop control to mimic the characteristics of the synchronous generator (SG) [16]. Compared with traditional vector control, VSG control could provide inertia support to the system and avoid the negative damping effect caused by the phase-locked loop [17]. Hence, VSG control is widely used in new energy grid-connected operation under the background of the high proportion of power electronic equipment. However, due to the influence of its own control structure, the problem of power coupling oscillation will be generated inevitably [18]. To solve the problem of low-frequency oscillation caused by power coupling, supplementary damping controllers have been introduced into the power loop [19]. This method can improve the damping characteristics and the stability of the system, but it cannot change the power coupling effect of VSG. Considering the power decoupling, the virtual impedance was introduced to eliminate the coupling effect, which can suppress the low-frequency oscillation effectively [20]. In References [21, 22], the virtual impedance was used to improve the transient stability in the microgrid. Reference [23] uses virtual impedance to change the output impedance of the HVDC system and improve the transmission capacity of the system. In addition, the virtual impedance could be used to improve the quality of output current [24]. The virtual impedance can suppress the reactive circulation by keeping the symmetry of the line impedance in the multi-VSG system and improve the stability of the system [25]. In Reference [26], the virtual impedance is combined with inertia to improve the response speed of system frequency. To achieve accurate voltage control, voltage-current double closed loop was often introduced into VSG. But the subsynchronous oscillation will occur easily when interacting with the weak grid. It can be known from Table 1 that the virtual impedance could suppress the power coupling oscillation effectively within several hertz. However, there are limitations in the suppression of subsynchronous oscillation from several hertz to 50 Hz. To solve the problem, a dynamic voltage compensation control method based on ADRC was proposed. This method could suppress the subsynchronous oscillation in a certain frequency range and show strong robustness. Finally, the effectiveness of the proposed oscillation suppression strategy was verified by time domain simulation.

Compared to the previous research, the main contributions of this article are as follows: First, the mechanism of virtual impedance improving the stability of grid-forming converter is revealed from the perspective of power coupling. Second, the analysis results based on the sequence impedance model provide a theoretical basis for the study of virtual impedance to improve the stability in the subsynchronous range. Third, the ADRC controller parameters were designed by the bandwidth method, which could suppress the subsynchronous oscillation in a certain frequency range and show strong robustness. And it could provide a reference for the study of oscillation suppression methods.

This paper is constructed as follows: The topology and control scheme of grid-forming D-PMSG is introduced in Section 2. In Section 3, the influence of virtual impedance on system stability is analysed. And the dynamic voltage compensation control based on ADRC is designed in Section 4. In Section 5, the effectiveness of the proposed oscillation mitigation strategy is verified by simulation experiments. Finally, Section 6 concludes this article.

2. Topology and Control Scheme of Grid-Forming D-PMSG

2.1. Mathematical Model of PMSG

The stator voltage equation in three-phase stationary reference frames can be expressed as formula (1) by the Kirchhoff circuit theory and electromagnetic induction law. where is stator winding resistance; , , and are stator terminal voltage, respectively; , , and are stator current, respectively; and , , and are three-phase winding flux linkages, respectively.

Stator flux equation can be written as follows: where is rotor flux linkage, is equivalent excitation inductance, and is leakage inductance.

The electromagnetic torque equation can be obtained as follows: where is the number of pole pairs, is stator flux vector, and is stator current vector.

2.2. The Basic Principle of VSG Control

The rotor-side converter of grid-forming D-PMSG adopts voltage-current double closed loop control strategy. The grid-side converter adopts grid-forming VSG control strategy, which includes power outer loop and voltage-current inner loop control. Its main topology and control structure are shown in Figure 1.

The rotor motion state equation is introduced in VSG control which is based on the droop control in active power loop to mimic the transient characteristics of the SG. It could participate in system frequency regulation and provide inertial support to improve the frequency stability of the power system. The excitation voltage equation is also applied in the reactive power loop to regulate the bridge arm potential of the grid-side converter, and the terminal voltage can be accurately controlled by voltage and current double closed loop.

3. Sequence Impedance Modeling and Stability Analysis of Grid-Side Converter-Based VSG

3.1. Mechanism Analysis of Virtual Impedance to Improve the Stability of Grid-Side Converter

The output power transmission diagram of the distributed power is shown in Figure 2. The voltage is the public bus voltage. is the output voltage of distributed power. and are the equivalent impedance of the output line and virtual impedance. and are power angle differences and impedance angle.

According to Figure 2, the mathematical model of power transmission can be written as

It can be seen from formula (4) that there is a strong coupling relationship between active power and reactive power due to the impedance angle. This coupling characteristic will affect the stability of the system and even generate power coupling oscillation. To realize the decoupling of active and reactive power, it is necessary to keep the line mainly inductive. It means that the impedance angle θ needs to be set to .

The power angle is usually quite small; then, sin . Therefore, formula (4) can be simplified to

From formula (5), it can be seen that the active power is related to the power angle, and the reactive power is related to the voltage, and both of them are decoupled. The traditional droop control and VSG control are based on formula (5). In practice, the output line impedance is difficult to guarantee in pure inductive but contains resistive components. The virtual impedance can be used to keep mainly inductive, as shown in Figure 2. Therefore, the mechanism of virtual impedance improving the grid-forming converter stability is based on power decoupling.

3.2. Sequence Impedance Modeling of Grid-Side Converter

The sequence impedance model of grid-side converter was established based on harmonic linearization theory [27]. The output voltage, output current, and filter inductor current of phase A are expressed as (A.1) and (A.2).

3.2.1. The Model of Power Loop

In the axis reference frames, the output power can be expressed as where , , , and are the system output voltage and output current in axis reference frames, respectively. The frequency domain expressions are given in (A.3) and (A.4).

Using the convolution theorem in the frequency domain and ignoring the double frequency component, the expressions of active power and reactive power in the frequency domain can be obtained as follows: where denotes complex conjugation.

According to the motion state equation of rotor, the expression of VSG output phase angle is expressed as

Substituting (7) into (9) and neglecting quadratic term of small signal, the frequency domain expression of VSG output phase angle can be obtained as where , , and ; is the voltage fundamental phase, and is the initial phase.

Similarly, according to excitation voltage equation, the frequency domain expression of bridge arm potential amplitude can be obtained as follows: where .

3.2.2. The Model of Voltage-Current Double Closed Loop

Frequency domain expression of voltage out loop reference value is given as where and are voltage drop of virtual impedance, and it can be expressed in frequency domain as follows:

According to (14), converting the three-phase output voltage signal to two phase dq axis reference frame. The adopts VSG output phase angle. Frequency domain expressions of the trigonometric function and the Park transformation matrix are shown in (A.5)–(A.11).

Applying the Laplace transform to (14), the frequency domain expression of output voltage in dq axis reference frames can be obtained as

Similarly, the frequency domain expression of the inductor current can be obtained (see (A.12) and (A.13)).

3.2.3. Positive and Negative Sequence Impedance Expressions

Modulation signal expression in dq axis reference frames can be expressed as where

The bridge arm potential expression of A phase is given as follows: where is the bridge arm potential, is the DC side voltage, is the modulation coefficient, and is the modulation signal of phase A after the inverse Park transform for (17).

The current frequency response equivalent circuit is shown in Figure 3, and the current frequency response expression can be expressed as follows:

The sequence impedance in frequency domain can be expressed as

Neglecting coupling effect, the positive and negative sequence impedance expressions can be obtained as where

3.3. Stability Analysis for VSG in Weak Grid

The stability characteristics of grid-side converter-based VSG are shown in Figure 4. It can be seen from the phase-frequency curve that the negative sequence impedance of the grid-side converter presents inductive within 100 Hz, while the positive sequence impedance presents capacitive characteristics as a whole. Especially in the range of 20 Hz-44 Hz, the positive sequence impedance phase is lower than -90o which shows the characteristics of negative damping and capacitance. It can be seen from the amplitude-frequency curve that there is an intersection. It is between the positive sequence impedance curve and the weak grid impedance curve at 30 Hz, and the system loses stability when the phase difference exceeds 180o. The equivalent capacitance of the grid-side converter interacts with the inductive characteristics of weak grid and generates electrical resonance, which is easy to diverge and produce subsynchronous oscillation under the impact of negative damping. Therefore, the grid-forming D-PMSG has the risk of instability under weak grid conditions.

3.4. Influence of Virtual Impedance on Grid-Side Converter Stability

It can be known from the literature analysis that the virtual impedance can suppress the low-frequency oscillation effectively caused by power coupling within several hertz, but there are few studies focused on the suppression of subsynchronous oscillation from several hertz to 50 Hz. The sequence impedance model of the grid-side converter was established, and the influence of the virtual impedance on the system stability in the subsynchronous frequency range was analysed.

3.4.1. Influence of on Grid-Side Converter Stability

The influence of on grid-side converter stability is shown in Figure 5. The resonance frequency increases with which increases from the amplitude-frequency curve. Due to the weakening of the system capacitance, the intersection point with the grid impedance shifts slightly to the right, which is beneficial to stay away from the negative damping interval. But the overall change is relatively small. The phase-frequency curve shows that the negative damping interval of the system expands obviously with the increase of within 100 Hz. The negative damping interval is expanded obviously, and the system stability is weakened. It will further aggravate the risk of subsynchronous oscillation. Above 100 Hz, the positive and negative sequence impedances show transient inductive and negative damping characteristics, and the curves gradually overlap. The system near 1000 Hz is capacitive again due to the impact of the output filter capacitor. Therefore, the subsynchronous oscillation is not suppressed by increasing the virtual inductance.

3.4.2. Influence of on Grid-Side Converter Stability

The influence of on the grid-side converter stability is shown in Figure 6. The positive and negative sequence impedance amplitude gains decrease with the increase of the virtual resistance from the amplitude-frequency curve, but the overall change is relatively small. The phase-frequency curve shows that the positive sequence impedance within 100 Hz is capacitive and decreases with the increase of , while the negative sequence impedance is inductive and decreases with the increase of . Especially, the negative damping interval of positive sequence impedance reduced significantly with the increase of from the initial 20 Hz-44 Hz frequency range. It reduces the risk of subsynchronous oscillation greatly in grid-forming D-PMSG and improves system stability. The positive and negative sequence impedance curves above 100 Hz gradually coincide and are inductive. The system near 1000 Hz is capacitive again due to the influence of the output filter capacitor. Therefore, increasing the virtual resistance could reduce the risk of subsynchronous oscillation in grid-forming D-PMSG.

3.5. Influence of Virtual Resistance on Power Coupling

Although increasing the virtual resistance is beneficial to suppress the system oscillation for grid-forming D-PMSG in weak grid, it will deepen the coupling degree between active power and reactive power and even further affect the system stability.

According to formula (4), a small disturbance model corresponding to the steady state can be written as where and and are the bridge arm potential and power angle at stable operating point.

To analyse power coupling degree quantitatively, the relative gain matrix analysis method (RGA) is used [28]. According to Reference [20], the relative gain matrix can be listed as where represent the coupling coefficient .

The relationship between coupling coefficient and virtual resistance is shown in Figure 7. When the power angle and virtual inductance remain constant, the coupling coefficient between active and reactive power increases as the virtual resistance increases. Therefore, increasing the virtual resistance will deepen the power coupling and even cause power coupling oscillation.

From above analysis, it can be known that increasing virtual inductance could weaken the system capacitance to a certain extent and shift the resonant frequency to the right, but the shift amplitude is limited. And the negative damping interval expands significantly with the virtual inductance increasing, which further aggravates the risk of subsynchronous oscillation, while increasing the virtual resistance can narrow the negative damping interval significantly and reduce the risk of subsynchronous oscillation. However, increasing the virtual resistance will deepen the coupling degree between active power and reactive power and even increase the risk of power coupling oscillation. There are limitations on the subsynchronous oscillation suppression by increasing the virtual resistance. Hence, a dynamic voltage compensation control method based on ADRC was proposed in this paper.

4. The Design of Dynamic Voltage Compensation Control Based on ADRC

4.1. The Basic Principles of ADRC

Active disturbance rejection control has strong robustness and high control precision, which is an effective way to solve the problem of nonlinear systems [29, 30]. Taking the second-order control system as an example, the ADRC mainly includes three parts: nonlinear differential tracker, extended state observer, and nonlinear state error feedback. The basic control structure of ADRC is shown in Figure 8.

Second-order system can be expressed as follows: where and are the state variables of the controlled object; is an unknown function, which contains the lumped disturbance and the uncertain part of the system; is the system control input; is the system control output; and is the feedback coefficient. Let , .

The ADRC controller is designed according to the following steps: (1)Arranging the transition process with the set value as inputwhere is the input of the ADRC system, and are the output of tracking differentiator, is the velocity factor, is the filtering factor, is the simulation step length, and is the fastest feedback function. (2)The extended state observer (ESO) uses output and input to track the system state variables and and sets the estimated value of the state variable to and . The disturbance variable is extended to a new state variable , which is estimated to be where

, , and are a set of parameters (3)State error feedback law is presented as follows:where is a set of parameters (4)Disturbance compensation process is presented as follows:where is the compensation factor

4.2. The Parameter Design of ADRC

The error equation (33) can be obtained by subtracting (26) from (29).

Take,,, and. Using bandwidth method [30] and applying the Laplace transform to equation (33), the transfer function between lumped disturbance error and disturbance can be expressed as where , , and ; is the bandwidth of extended state observer.

It can be seen from Figure 9 that the amplitude gain curve shifts to the right with the increase of and the disturbance attenuation interval expands in the lower frequency band. The amplitude gain is required to be negative at 50 Hz, and the lumped disturbance error converges to stable quickly in 0 to 50 Hz. When the bandwidth is larger, the accuracy of the extended state observer is higher, and the influence of noise is also greater. Considering comprehensively, the bandwidth of the extended state observer is set as 400 rad/s; as the compensation factor is mainly used to regulate the compensation voltage which is set as 14.3.

The dynamic voltage compensation control structure based on ADRC is shown in Figure 10. The inductor current is used as control variable, and the compensation voltage is used as control input. To suppress the subsynchronous oscillation and improve the robustness of the system, the dynamic compensation voltage is generated continuously by observing the harmonic component of the inductor current. The dynamic compensation voltage without changing the virtual resistance parameter was used to offset the harmonic voltage drop in the virtual impedance part and suppress the subsynchronous oscillation well in weak grid. This method solves the problem of power coupling oscillation caused by increasing the virtual resistance.

5. Simulation Experiment Results

To verify the effectiveness of the proposed control strategy, a model of 3 MW grid-forming D-PMSG is established on MATLAB/Simulink software platform. The simulation parameters are shown in Table 2.

5.1. Simulation Results under Different Grid Strength Conditions

The grid impedance was switched from normal value to 0.17 mH at 0.8 s; then, was used to simulate the weak grid. The oscillation of output current occurred at 30 Hz, while the oscillation of active and reactive power occurred at 20 Hz, as shown in Figure 11. The dynamic voltage compensation control was started at 1.8 s, and the grid inductance was switched back to normal at 3 s.

It can be seen from Figure 12(a) that after the dynamic voltage compensation control was started at 1.8 s, the active power and reactive power of the system were attenuated to stable quickly from the maximum oscillation amplitude of 0.33 MW and 0.38 Mvar. And the output current was also stabilized through several cycles of adjustment. The grid inductance was switched to 0.22 mH to simulate a weaker grid at 0.8 s. As shown in Figure 12(b), the oscillation amplitude of the system became larger than before. But after the dynamic voltage compensation control started at 1.8 s, the system oscillation gradually disappeared. Although the oscillation amplitude and frequency changed with the weakening of grid strength, the system still converged to stable through a short adjustment. Therefore, the dynamic voltage compensation control based on ADRC can suppress the subsynchronous oscillation in a certain grid strength range, which shows strong adaptability in weak grid.

5.2. Simulation Results of Terminal Voltage Step Disturbance

The subsynchronous oscillation in grid-forming D-PMSG occurred when the grid inductance was switched from normal to 0.22 mH at 0.8 s. The dynamic voltage compensation control was started, and 3.5% rated voltage step disturbance was applied at 1.8 s. The disturbance was removed at 2.7 s; then, the grid inductance was switched back to normal at 3 s. It can be seen from Figure 13 that the system output power and current oscillation were still suppressed effectively under voltage disturbance conditions. The reactive power drops down 0.77 Mvar which is affected by the droop coefficient and still converged to stable rapidly and shows strong anti-interference ability.

5.3. Simulation Results of Active Power Command Step Disturbance

The grid inductance was switched from normal to 0.22 mH at 0.8 s, and then, the subsynchronous oscillation of the system occurred. The dynamic voltage compensation control was started at 1.8 s, and 1 MW active power command step disturbance was applied. The disturbance was removed at 2.7 s, and then, the grid inductance was switched back to normal at 3 s. It can be seen from Figure 14 that the output power and current oscillation of the system converged to stable by short adjustment, which is only a few cycles of adjustment time more than that without disturbance and shows good system robustness.

5.4. Simulation Results of Increasing Virtual Resistance

The grid inductance was switched to 0.22 mH at 0.8 s, and then, the subsynchronous oscillation of the system occurred. The dynamic voltage compensation control was started, and 0.024 Ω virtual resistance was applied at 1.8 s. The virtual resistance continuously increased to 0.05 Ω at 2.5 s. The resistance was decreased to 0.024 Ω, and the grid inductance was switched back to normal at 3.0 s. It can be seen from Figure 15 that the system oscillation has a faster convergence speed and better stability than the case without virtual resistance during 1.8 to 2.5 s. The power curve diverged to form a new oscillation with the virtual resistance increasing to 0.05 Ω after 2.5 s. Therefore, increasing the virtual resistance within a certain range is beneficial to improve the stability of grid-forming D-PMSG. However, increasing the virtual resistance continuously will deepen the power coupling effect and form a new power coupling oscillation, which is not conducive to the stability of the system. The oscillation suppression method on ADRC and virtual resistance is compared in Table 3.

The above simulation results are consistent with the theoretical analysis. There are limitations in solving the oscillation suppression problem by simply increasing the virtual resistance. The dynamic voltage compensation control based on ADRC can effectively mitigate the subsynchronous oscillation for grid-forming D-PMSG in weak grid and avoid the problems of power coupling oscillation caused by increasing virtual resistance.

From the simulation analysis, it can be seen that the virtual resistance can improve the oscillation suppression effect, but the continuous increase will deepen the power coupling and form a new power coupling oscillation. Meanwhile, this method needs to change the virtual resistance parameter frequently, which shows weak adaptive ability and poor robustness. The dynamic voltage compensation control can suppress the subsynchronous oscillation in a certain grid strength range without changing the virtual resistance parameter. However, the harmonic current oscillation amplitude increases, and the dynamic compensation control accuracy decreases in the face of a weaker grid. Combining the virtual resistance and dynamic voltage compensation control can not only further improve the oscillation suppression effect but also maintain sufficient control accuracy in the face of a weaker grid and reduce the pressure of ADRC parameter design greatly. Therefore, how to determine the virtual resistance boundary range of power coupling oscillation and combine ADRC with virtual resistance to improve the effect of subsynchronous oscillation suppression is the direction of further research.

6. Conclusion

In this paper, the influence of the virtual impedance on the subsynchronous oscillation caused by the interaction between the grid-forming D-PMSG and the weak grid was analysed, and a dynamic voltage compensation method based on ADRC was proposed. The following conclusions are obtained by modeling and simulation: (1)This paper revealed the mechanism of virtual impedance for improving the stability of grid-forming converter from the perspective of power coupling, and the influence of the virtual impedance on subsynchronous oscillation was analysed by the sequence impedance model(2)The analysis results show that increasing the virtual inductance will enlarge the negative damping interval, which is not conducive to system stability. And increasing the virtual resistance can greatly reduce the negative damping interval and improve system stability. However, increasing virtual resistance continuously will deepen the power coupling degree and cause a new power oscillation. Therefore, there are limitations in suppressing the subsynchronous oscillation for grid-forming D-PMSG in weak grid by changing the virtual resistance(3)To solve the problem, a dynamic voltage compensation control method based on ADRC was proposed. The ADRC controller parameters were designed by the bandwidth method, which can suppress the subsynchronous oscillation in a certain frequency range and show strong adaptability and good system robustness in weak grid

Appendix

The output voltage, output current, and filter inductor current of phase A are expressed as follows:

The frequency domain expression of (A.1) can be written as where , , , , , , , , and .

Frequency domain expressions of system output voltage and current in axis reference frames are expressed as follows:

According to the trigonometric function relationship, (A.5) can be obtained.

Applying the Laplace transform to (A.5), the trigonometric function can be written as

The expression of the Park transformation matrix is expressed as follows:

The Park transformation matrix is linearized as follows: where and are perturbation component and fundamental component, and the matrix expressions are given as follows:

Inductance current frequency domain expression can be obtained as

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 51967016), Natural Science Foundation of Inner Mongolia Autonomous Region (no. 2022LHQN05001), Inner Mongolia Autonomous Region Key R&D and Achievement Transformation Program (no. 2023YFHH0077), and Basic Scientific Research Expenses Program of Universities directly under Inner Mongolia Autonomous Region (nos. JY20220094 and JY20230066).