Journal of Modern Power Systems and Clean Energy

ISSN 2196-5625 CN 32-1884/TK

网刊加载中。。。

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

确定继续浏览么?

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

Compound Compensation Control for Improving Low-voltage Ride-through Capability of Virtual Synchronous Generators  PDF

  • Zhiyuan Meng 1
  • Xiangyang Xing 1 (Member, IEEE)
  • Xiangjun Li 2 (Senior Member, IEEE)
  • Jiadong Sun 1
1. School of Control Science and Engineering, Shandong University, Shandong, China; 2. China Electric Power Research Institute, Beijing, China

Updated:2025-05-27

DOI:10.35833/MPCE.2024.000404

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

Abstract

The virtual synchronous generator (VSG), utilized as a control strategy for grid-forming inverters, is an effective method of providing inertia and voltage support to the grid. However, the VSG exhibits limited capabilities in low-voltage ride-through (LVRT) mode. Specifically,the slow response of the power loop poses challenges for VSG in grid voltage support and increases the risk of overcurrent, potentially violating present grid codes. This paper reveals the mechanism behind the delayed response speed of VSG control during the grid faults. On this basis, a compound compensation control strategy is proposed for improving the LVRT capability of the VSG, which incorporates adaptive frequency feedforward compensation (AFFC), direct power angle compensation (DPAC), internal potential compensation (IPC), and transient virtual impedance (TVI), effectively expediting the response speed and reducing transient current. Furthermore, the proposed control strategy ensures that the VSG operates smoothly back to its normal control state following the restoration from the grid faults. Subsequently, a large-signal model is developed to facilitate parameter design and stability analysis, which incorporates grid codes and TVI. Finally, the small-signal stability analysis and simulation and experimental results prove the correctness of the theoretical analysis and the effectiveness of the proposed control strategy.

I. Introduction

DISTRIBUTED generation systems (DGSs), consisting of renewable energy sources and grid-connected inverters, represent a booming technology for addressing the energy crisis problems [

1]. In recent years, the voltage source inverters (VSIs) are increasingly integrated into the power grid [2]. However, the electrical equipment often exhibits low inertia and weak damping characteristics, thereby posing a threat to the stable operation of the power system [3]-[5]. Thereupon, the virtual synchronous generator (VSG), as a form of grid-forming (GFM) inverter, is proposed to provide inertia and voltage support for the grid [6]-[8], thereby ensuring the stable operation of VSIs.

The VSG aims to emulate the characteristics of the synchronous machines in the VSIs. However, they exhibit limitations during the grid faults, as the inverters are not built to handle prolonged overcurrent situations, reflecting a comparatively weaker overcurrent capability than that of traditional synchronous machines [

9].

To equip the VSG with fault ride-through capabilities during the grid faults, various low-voltage ride-through (LVRT) strategies are developed to tackle challenges like overcurrent and rapid power tracking. Presently, two primary LVRT control strategies for VSG are prevalent. One LVRT control strategy transitions the VSG into a grid-following (GFL) inverter during the grid faults [

10], [11]. The other LVRT control strategy incorporates additional control loops specifically for compensating internal potential or power angle [12]-[15].

Although GFL inverters have advantages in current limiting and tracking speed [

16]-[18], mode switching causes the VSG to lose the voltage source characteristics. The use of phase-locked loops (PLLs) in GFL inverters might induce instability in systems with high penetration of renewable energy [19]. Additionally, to mitigate the effects of mode switching, a seamless switching control strategy is essential to maintain state variable consistency before and after the grid faults. This necessitates the simultaneous implementation of both GFL and GFM control modes, as detailed in [11], which significantly increases the computational demands on the controller.

The transient virtual impedance (TVI) control strategy is proposed in [

12], [13], which introduces virtual impedance (VI) to effectively reduce output voltage of the VSG, thereby facilitating current limitation. Importantly, the TVI is designed to gradually decrease to zero, minimizing any negative impact on the system.

Despite significant research into LVRT control and small-signal modeling of TVI, the practical application of TVI still encounters several challenges. Firstly, while TVI effectively limits current, it does not enhance the power response speed of VSG during the grid faults, which is in contradiction with the existing grid codes. For example, [

20] mandates that the reactive power adjustment time for VSG should not exceed 60 ms when a symmetrical fault occurs in the grid voltage. Secondly, grid codes dictate that the power output of VSG during the grid faults should be a function of the output voltage [20], [21]. This dynamic adjustment is frequently overlooked in both small-signal [22] and large-signal [12] models. To address these problems, TVI must be integrated with additional control strategies for effective fault ride-through capability, and the mathematical modeling of VSG must align with the stipulations of grid codes.

TVI can be regarded as an indirect form of internal potential compensation (IPC). References [

14] and [15] explore direct IPC techniques. In [14], the output voltage of the VSG is directly incorporated into the control loop as a compensation term, which effectively addresses overcurrent problem caused by the grid faults. However, [14] does not directly control reactive power but adjusts the output voltage instead, which fails to meet the reactive power support standards required by the grid codes. Reference [15] introduces a refined IPC formula that utilizes output voltage and phase angle, and proposes a power angle compensation strategy that reacts to changes in the power angle during the grid faults. This strategy successfully mitigates steady-state overcurrent problem and decreases the regulation time of the power loops. However, [15] only resolves the steady-state overcurrent and does not suppress the transient current generated during the grid faults and fault recovery. Furthermore, the formulas for compensation in [15] are based on the condition of maximum output current, which may lead to inaccurate compensation. Additionally, the LVRT strategy described in [15] necessitates extra voltage sensors, thereby increasing equipment costs.

In summary, the current research problems can be summarized as insufficient analysis of the impact of the grid faults on the dynamic performance of VSG, inadequate accuracy of internal potential and power angle compensation, and a lack of large-signal and small-signal modeling that accounts for grid codes.

Hence, this paper introduces a compound compensation control strategy for VSG to address the outlined problems, which is aimed to improve LVRT performance of VSG. The main contributions of the proposed control strategy in this paper are as follows:

1) A virtual power angle is defined to derive accurate internal potential and direct power angle compensation (DPAC) without additional measurement equipment. The analysis via the transfer function of power loop shows that the grid faults extend the regulation time and increase overcurrent durations. Consequently, a compound compensation control strategy is proposed to enable rapid tracking of the power reference during the grid faults.

2) The effects of TVI control parameters on the dynamic performance and stability of VSG are examined through large-signal model.

3) A small-signal model is developed, which considers the reactive power support characteristics required by grid codes. The feasibility of the developed small-signal model is demonstrated through eigenvalue matrix analysis.

The rest of this paper is organized as follows. Section II derives the principle of VSG. Section III proposes the composite compensation control strategy. Sections IV and V show the simulation validation and experimental validation. Section VI concludes this paper.

II. Principle of VSG

A. VSG Control Structure

The active power control loop (APCL) is designed to generate the internal potential phase [

23], which is given as:

ω=ωn+1Js+DpPref-Pωndθdt=ωdδgdt=ω-ωn=Δω (1)

where J is the virtual moment of inertia; Dp is the damping coefficient; ω and ωn are the output angular frequency of the VSG and the rated angular frequency of the grid, respectively; P is the active power delivered to the point of common coupling (PCC); Pref is the active power reference; θ is internal potential phase; and δg is the power angle.

The reactive power control loop (RPCL) is used to generate the internal potential signal, which is expressed as:

KdMdt=Qref-Q+Dq(UCn-VCm) (2)

where Dq is the voltage-droop coefficient; VCm and UCn are the PCC voltage magnitude and rated internal potential magnitude, respectively; Q is the reactive power delivered to the PCC; Qref is the reactive power reference; M is the variation of the internal potential relative to UCn; and K is the inertia coefficient.

Moreover, as shown in Fig. 1, the internal potential Edref=M+UCn, and Eqref=0, where Eqref is the q-axis reference voltage of VSG. The voltage-droop coefficient Dq affects the reactive power reference during the grid faults. Thus, it affects the reactive power output. Therefore, we set Dq=0 in the further analysis. Considering the RPCL is governed by an integrator, it is feasible to achieve the reactive power output target.

Fig. 1  Conventional VSG control structure.

The power delivered to the PCC can be derived as:

P=1.5(vCdisd+vCqisq)Q=1.5(vCqisd-vCdisq) (3)

where vCd and vCq are the PCC voltages on the synchronous dq coordinates, respectively; and isd and isq are the inverter-side currents on the synchronous dq coordinates, respectively.

Then, the PCC voltage magnitude VCm can be expressed as:

VCm=vCd2+vCq2 (4)

The VI, consisting of virtual resistance Rv and virtual inductance Lv, is designed to generate reference PCC voltages vCdref and vCqref [

24]. The voltage and current control loops (VCCLs) are used to produce modulated signal eabc* for VSI. Subsequently, pulse width modulation (PWM) pulses are generated by sinusoidal PWM (SPWM).

Figure 1 shows the conventional VSG control structure, where Cf is the filter capacitor; Lg is the line inductance; isabc is the inverter-side current; vCabc is the PCC voltage; and vgabc is the grid voltage.

B. Reactive Power Support Characteristics

The magnitude of reactive current injected into the grid depends on the severity of PCC voltage sags [

20], [25], which is described as:

Iqref=0                                  VCmUCn0.9-kIN0.9-VCmUCn    0.2VCmUCn<0.9-1.05IN                      VCmUCn<0.2 (5)

where Iqref is the reactive current reference; k is the current gain and is set to be 1.5; and IN is the nominal grid current magnitude.

Then, the active current injected into the grid is given as:

Idref=IN2-Iqref2    IN2-Iqref200                       IN2-Iqref2<0 (6)

where Idref is the active current reference.

Thereupon, the power reference of the VSG can be derived as:

Pref=1.5VCmIdrefQref=-1.5VCmIqref (7)

A power-oriented reference generation method for VSG during LVRT is presented in (7). When the VI is applied, (7) remains unchanged, as mentioned in [

26]. Additionally, to ensure that the current during the grid faults is consistent with (5) and (6), a current-oriented power reference generation can be used, which is discussed in detail in Supplementary Material A.

When Dq=0, it means that the RPCL loses droop control. Thus, the RPCL can be expressed as:

KdMdt=Qref0+(0.9UCn-VCm)1.5VCmkINUCnDqe-Q (8)

where Qref0 is the reactive power reference under the normal condition, which is set to be 0; and Dqe is the equivalent damping coefficient.

The RPCL expression for a conventional VSG is expressed as:

KdMdt=Qref0+(UCn-VCm)Dq-Q    Dq=QmaxkqUCn (9)

where kq[0.05,0.10] [

6], [27]; and Qmax is the maximum reactive power reference.

Compared with (9), (8) shows that during the grid faults, the damping coefficient becomes a function of the PCC voltage. Under the normal condition, the VSG outputs little reactive power and the damping coefficient has a minimal effect. Therefore, even though Dq is set to be 0 in this paper, the VSG still possesses an equivalent damping coefficient Dqe.

To further demonstrate the rationale for setting Dq to be 0, simulation validation is conducted. The conventional VSG in (9) allows the VSG to support the PCC voltage during the grid faults. However, this causes the VSG to output excessive reactive power, resulting in overcurrent and inaccurate reactive power control, as illustrated in Fig. 2, where Qmax=10 kvar, kq=0.1, and Vg and Vgn are the amplitude of grid voltage and its normal value, respectively.

Fig. 2  Waveforms of inverter-side current and output power with conventional VSG control (Dq=Qmax/(kqUCn)).

Furthermore, [

26] indicates that the voltage mode characteristics of a GFM inverter can still be maintained by setting Dq to be 0 during a grid fault with the power references determined by the grid codes. If it is necessary to maintain the reactive power droop characteristic of the VSG under the normal condition, the strategy in [26] can be used. The detailed process is shown in Supplementary Material B.

C. Performance of VSG During Grid Faults

The grid faults are categorized into symmetrical and asymmetrical faults, with symmetrical faults posing the greatest severity. Therefore, this paper primarily focuses on symmetrical faults.

The equivalent impedance from the internal potential signal to the PCC can be expressed as:

Zeq=Req2+Xeq2=Rv2+(ωnLv)2=ZvθI=arctan ωnLfRv (10)

where Lf is the filter inductance; Req, Xeq, and Zeq are the equivalent resistance, reactance, and impedance, respectively; θI is the impedance angle; and Zv is the VI.

Lv is set equal to Lf in this paper, implying that no extra equivalent inductance is added. Consequently, Xeq=ωnLv=Xv, and Zeq equals the VI, i.e., Zv.

The circuit model of VSG in Fig. 1 is simplified as Fig. 3. The active power P and reactive power Q are calculated as:

Is=12Erefmδ-VCm0°ZvθIP=Re3VCM0°2I¯s=32VCmErefmZvsin(δ+θI)-VCm2RvZv2Q=Im3VCM0°2I¯s=32VCmErefmZvcos(δ-θI)-VCm2XvZv2 (11)

Fig. 3  Circuit model of VSG.

where bold font and superscript ()¯ represent the vector and conjugate vector, respectively; Erefm is the internal potential output by the VSG; Is is the vector of inverter-side current; δ is the phase difference between the internal potential and PCC voltage; and Re() and Im() are the real and imaginary parts, respectively.

It is worth noting that δ in (11) represents a virtual power angle, rather than the actual power angle. The calculation of the actual power angle requires data on grid voltage and grid impedance [

23], [24]. Acquiring this information typically demands additional measurement devices, which poses practical challenges and increases the costs. Therefore, this paper derives the expression of the virtual power angle by using VI and PCC voltage, which eliminates the need for extra hardware, allowing the calculations to be performed entirely through the software.

At the steady state, we set P=Pref and Q=Qref. To preserve the external characteristics of the circuit, the VI is set as Rv<<Xv. Thus, (11) can be rewritten as:

P32ErefmVCmXvsin δQ32ErefmVCmcosδ-VCm2Xv (12)
δarctan2PrefXv2QrefXv+3VCm2Erefm2XvQref+3VCm23VCmcosδ (13)

This adjustment enables the derivation of Erefm and δ by means of the power references and the PCC voltage.

At the steady state, the small-signal model of the VSG is given as [

22]:

P^32EsVCmcosδsXvδ^+32VCmsinδsXvE^Q^-32EsVCmsinδsXvδ^+32VCmcosδsXvE^ (14)

where the superscript  ^ represents the small perturbation; and Es and δs are the value of the internal potential and the power angle, which are equal to Erefm and δ calculated by (13), respectively.

Neglecting the coupling between APCL and RPCL, the loop gains are obtained by combining (1), (2), and (14), which can be given as:

Tp(s)=32EsVCmcosδsXv1(Js+Dp)ωn1sTq(s)=32VCmcosδsXv1Ks (15)

where Tp(s) and Tq(s) are the loop gains of APCL and RPCL, respectively.

The bode diagram of the power control loops under different PCC voltages is shown in Fig. 4. Figure 4 illustrates that the cutoff frequencies ωcp of the APCL and ωcq of the RPCL decrease as the PCC voltage sags, signifying a slowdown in the response speed of the power control loops during the grid faults.

Fig. 4  Bode diagram of power control loops under different PCC voltages. (a) Bode diagram of APCL. (b) Bode diagram of RPCL.

III. Compound Compensation Control Strategy

The compound compensation control strategy proposed in this paper comprises adaptive frequency feedforward compensation (AFFC), DPAC, IPC, and TVI. The proposed control strategy for LVRT is shown in Fig. 5.

Fig. 5  Proposed control strategy for LVRT.

A. AFFC

The proposed control strategy employs AFFC to increase the cutoff frequency of the APCL during the grid faults, thereby reducing the adjustment time of the APCL. Figure 6 presents the small-signal control block diagram of AFFC.

Fig. 6  Small-signal control block diagram of AFFC.

The adaptive compensation coefficient kA is given as:

kA=UCn2-EsVCmEsVCm (16)

The gain loop of AFFC TAFFC(s) is derived as:

TAFFC(s)=3UCn2cosδs2s(Js+Dp)ωnXv3UCn22s(Js+Dp)ωnXv (17)

where δs is approximately 0 at the steady state by (4)-(7).

AFFC raises the cutoff frequency of the APCL during the grid faults to the same level as that before the grid fault, thereby reducing the adjustment time.

B. DPAC

The vector diagram of the system before and during the grid faults is shown in Fig. 7, where the subscripts 0 and 1 represent the physical quantity before and during the grid faults, respectively; δg is the phase difference between the internal potential and grid voltage; vgq is the axis with a lag of vg0=90°; δCd0 is the phase angle of PCC voltage during grid fault on the synchronous d0q0 coordinates; δc2 is the change in angle of the internal potential Eref1 relative to E0; and vC1d0 and vC1q0 are the d-axis and q-axis components of PCC voltage vC1 on the synchronous d0q0 coordinates, respectively.

Fig. 7  Vector diagram of system before and during grid fault.

Figure 7 shows that the introduction of a compensation angle δc2 upon the detection of a grid fault enables the VSG to swiftly achieve a steady state. δc2 is calculated as:

δc2=δ1-δCd0=δ1-arctanvC1q0vC1d0 (18)

where δ1 is calculated by (13).

In order to derive δCd0 in the d0q0 coordinates before the grid faults, Fig. 5 depicts an integrator positioned at post-switch S1. Upon the detection of the PCC voltage sags, S1 transitions from 0 to 1, preserving the phase angle signal before the grid fault.

In (18), performing an arctan calculation on the PCC voltage might appear to simulate an open-loop PLL. However, this does not transform the VSG in a GFL mode. At the steady state, both δc1 and δc2 are constants, with their time derivatives equal to zero. Therefore, they have no impact on the output angular frequency ω. Therefore, ω is governed by the APCL both before and during the grid faults. This feature is the characteristic of a GFM inverter.

C. IPC

The results in Fig. 4 indicate that PCC voltage sags lead to a deceleration in the response speed of the RPCL, which extends the time needed for the internal potential signal to stabilize, causing prolonged periods of overcurrent. To address this problem, the IPC is introduced.

Following the detection of the PCC voltage sags, (13) is used to compute the internal potential signal. Then, S4 in Fig. 5 is switched from 0 to 1 to implement the IPC. The IPC acts as adaptive feedforward compensation, significantly reducing the regulatory impact on RPCL.

D. TVI

The TVI in Fig. 5 is adopted to reduce the transient current. Studies investigate the impact of the TVI on the GFM inverters through small-signal analysis [

12], [22]. However, small-signal models fail to capture the nonlinear behavior of the inverters during the grid faults [22]. Therefore, this paper examines the effects of the TVI on VSG in the LVRT mode using large-signal model. To highlight the significance of the TVI, a comprehensive large-signal model that incorporates only the TVI and grid codes are developed.

The HPF is expanded to define the state variable xVI, as shown in Fig. 8 and (19).

Fig. 8  Diagram of TVI control.

dxVIdt=-1TIxVI+1TIIsat (19)

where Isat=max (Im-Ith, 0), Im is the current magnitude, and Ith is the current threshold for starting TVI; and TI is the time constant of the high-pass filter (HPF).

Therefore, Rvt and Xvt are expressed as:

Rvt=dxVIdtTIkRXvt=Rvtσ (20)

where kR is the gain for TVI; Rvt and Xvt are the transient virtual resistance and inductance, respectively; and σ=Xvt/Rvt.

The gain kR is designed to limit the current during a bolted fault [

13], which is shown as:

UCn=Im,maxRvs2+Xvs2kR-b+b2-4ac2a (21)

where Im,max is the maximum current magnitude, which is set to be 1.5IN, i.e., Im,max=30 A; a, b, and c are shown in (SA1) in Supplementary Material A; and Rvs=Rv+Rv and Xvs=Xv+Xvt are the equivalent virtual resistance and equivalent virtual reactance, respectively.

Generally, σ should be greater than 3 to ensure that TVI exhibits inductive characteristics [

22]. The lower bound of kR is determined by setting σ=6. Therefore, kR is calculated as 0.2 Ω/A, and in subsequent analyses, kR might not be recalculated to simplify the analysis under various σ.

Generally, the bandwidth of the power loop is set much lower than that of the VCCL [

26]. Consequently, only the power loops and the TVI are considered during the establishment of the large-signal model of the VSG. The detailed large-signal model of the LCL circuit is derived in (SA2) in Supplementary Material A.

Then, the large-signal model of the VSG is formulated by combining (1), (2), (19), and (S2). The “ode45” command in MATLAB [

24] is used to solve the large-signal model with various TVI parameters during the grid faults. The current magnitude Im with various TVI parameters during the grid faults is illustrated in Fig. 9, where the grid voltage sags to 0.5Vgn at 2 s.

Fig. 9  Current magnitude with various TVI parameters during grid faults. (a) With various σ and without TVI. (b) With various TI.

Figure 9(a) reveals that without TVI, the VSG experiences significant transient current. As σ increases, there is a notable reduction in both the magnitude and oscillation duration of the transient current. However, an excessively high σ leads to a loss of the stability. Then, Fig. 9(b) illustrates that a moderate increase in TI can effectively mitigate the transient current and shorten the regulation time.

E. Description for Control Flow

Figure 10 presents the flowchart of the proposed control strategy. Upon detecting the PCC voltage sags (VCm0.9UCn), the VSG switches to the LVRT mode. Therefore, S2 and S3 are closed and S1 and S4 are switched from 0 to 1. This triggers the activation of AFFC, DPAC, and IPC, with TVI automatically engaging in response to the current levels. Following the fault recovery, the VSG is switched to the voltage recovery mode. Initially, S4 is switched from 0 to 1 to mitigate the overcurrent. A timer is then initiated to manage the transient process of the VSG and lessen the negative impacts of sudden removal of the compensation terms. The condition for the end of the voltage recovery mode is expressed as:

Fig. 10  Flow chart of proposed control strategy.

1) TTset, where T and Tset are the operation duration and preset time of the timer, respectively.

2) |Pref-P|Pth and |Qref-Q|Qth, where Pth and Qth are the active power and reactive power threshold, respectively.

After the voltage recovery mode ends, S1 is switched from 1 to 0. Yet, the power angle compensation δc2-δc1 continues because S2 and S3 remain engaged. As S1 is switched to 0, the input of the AFFC drops to 0. Thereupon, the power angle compensation value is used as feedback for the AFFC, facilitating a smooth exit from the power angle compensation. The specific process can be described as follows.

If δc1>δc2>0 and kc|δc2-δc1|>δth, we should close S5, where δth is the power angle threshold, and kc is the power angle compensation gain. Since the feedback quantity δc1-δc2 is positive, δc1 gradually decreases under the integrator in AFFC. Once |δc2-δc1|δth, we should open S2, S3, and S5, which means to end the power angle compensation smooth exit (PACSE) mode. The analytical approach for other conditions mirrors the above.

The flag for each switch is detailed in Supplementary Material B Table SBI.

Additionally, the proposed control strategy spans the entire LVRT process, rather than being limited to a specific segment. Consequently, the transition between these modes in Fig. 10 necessitates the switches and corresponding algorithms, inevitably increasing the complexity of the proposed control strategy. Despite this complexity, further sections present small-signal stability analysis and simulation results to demonstrate its feasibility and practicability.

IV. Simulation Validation

In this section, a small-signal model of the VSG is established to analyze the impact of the proposed control strategy on the stability of the VSG and to detail the design of control parameters. Based on this, the proposed control strategy is validated through simulations during both symmetrical and asymmetrical grid faults, and is compared with existing strategies.

A. Small-signal Stability Analysis and Design of Control Parameters

To ensure the generality of the results, the small-signal model is confined to the conditions where the PCC voltage sags within the range of [0.2, 0.9]UCn. By applying local linearization techniques, the small-signal modeling process is detailed in (SA3)-(SA10) in Supplementary Material A.

The selected analysis condition involves a grid voltage sag to 0.5Vgn. The trajectory of the eigenvalue variations of the state matrix formed by (SA3)-(SA10) in Supplementary Material A under different conditions is shown in Fig. 11.

Fig. 11  Trajectory of eigenvalue variations of state matrix under different conditions. (a) With different control strategies. (b) With proposed control strategy (σ increases from 6 to 30 with an interval of 2). (c) With proposed control strategy (TI increases from 0.01 to 1.00 with an interval of 0.05).

Specifically, the simulation results from large-signal model in Fig. 9 show that the VSG becomes unstable when σ is set to be 20 during the grid faults. Conversely, Fig. 11 demonstrates that even σ increases to 28, the VSG remains stable. As previously mentioned, the small-signal model cannot capture the nonlinear behavior induced by large disturbances. Therefore, the obtained stability boundaries may not be applicable. Nonetheless, the eigenvalue trajectories still reflect the dynamic performance of the VSG under different conditions, providing a basis for parameter design. Thereupon, σ is finally designed as 10 to balance the stability and VSG performance.

kR and σ are determined by the previous analysis. Then, this subsection continues to elaborate on the design of the remaining control parameters.

In deriving the power angle and IPC formulas, i.e., (13), it is assumed that Rv is much smaller than Xv with Xv=Xs=0.94 Ω. Therefore, Rv is designed as 0.02 Ω.

The current that the inverter can withstand should not exceed 1.5 times its rated current [

15]. This implies that the overcurrent suppression threshold Ith must be lower than this threshold, which is specified here as 24 A.

The peak fault current emerges approximately half a cycle post-fault and decays to steady state within approximately 10 cycles. Therefore, the time constant TI is set as half a grid voltage cycle, i.e., 0.01 s, which is used to suppress the maximum fault current. Upon the fault clearance, the preset time of the timer Tset is set to be 0.3 s, leaving sufficient time for the exit of the power angle compensation term.

Under the normal conditions, VSG power fluctuations are less than 0.05Pn [

28], where Pn is the rated power of 10 kW. Therefore, Pth and Qth are set to be 500 W and 500 var, respectively. Then, by using the power angle formula presented in (13), when the output power changes by 0.05Pn, (22) is derived to calculate the power angle variation, which is denoted as Δδ. In the PACSE mode, the calculated Δδ is approximately 0.003 rad when VCmUCn. Therefore, δth should be less than 0.003 rad. In addition, δth is designed to be 0.001 rad in this study.

Δδarctan2(1.05Pref)Xv3VCm2-arctan2PrefXv3VCm2 (22)

When initiating the PACSE mode, the VSG already reaches the steady state and δc20. Thus, the simplification of PACSE mode is shown in Fig. 12, where δc1,ini is the initial value of the integrator. Therefore, kc determines the duration of the PACSE mode. Considering the smoothness of this mode, kc is designed as 5.

Fig. 12  Simplification of PACSE mode.

The parameters J, Dp, and K are designed according to the guidelines provided in [

6] and [27]. Thus, they are not detailed in this paper.

B. Simulation Results During Symmetrical Grid Fault

To demonstrate the effectiveness of the proposed control strategy, a simulation verification is carried out. The simulation parameters are detailed in Table II.

TABLE II  Simulation Parameters
ParameterSymbolValue
Normal grid voltage magnitude Vgn 311 V
Nominal grid current amplitude IN 20 A
Rated angular frequency ωn 314 rad/s
Rated internal potential magnitude UCn 311 V
DC-link voltage Vdc 700 V
Filter inductance Lf 3 mH
Filter capacitance Cf 20 μF
Line inductance Lg 6 mH
Virtual moment of inertia J 0.06 kg·m2
Damping coefficient Dp 5 N·m·s
Inertia coefficient K 7 A·s

During the period when the grid voltage sags to 0.5Vgn, Fig. 13 illustrates the waveforms of PCC voltage, inverter-side current, and output power with the conventional VSG control. The results reveal that conventional VSG control results in considerable transient current and prolonged adjustment periods for both active power and reactive power during the symmetrical grid faults and voltage recovery mode, which complicates the codes for LVRT.

Fig. 13  Waveforms of PCC voltage, inverter-side current, and output power with conventional VSG control (Vg=0.5Vgn).

In contrast, Fig. 14 presents the simulation results with the proposed control strategy. The results illustrate the successful mitigation of transient current and swift achievement of target values for output power during the symmetrical grid faults. In addition, the PCC voltage can also quickly enter the steady state.

Fig. 14  Waveforms of PCC voltage, inverter-side current, output power, and switches with proposed control strategy (Vg=0.5Vgn).

Additionally, Fig. 14 shows that the timer stops at 1.3 s, initiating the PACSE mode, which ends by approximately 1.4 s.

The minimal impact on active power and current from PACSE mode proves the effectiveness of the proposed control strategy. Also, each switch operates accurately by the flowchart shown in Fig. 10.

The waveforms of IPC and power angle compensation are shown in Fig. 15, where Ecom is the output of IPC.

Fig. 15  Waveforms of IPC and power angle compensation (Vg=0.5Vgn).

In Fig. 15, Ecom is coincident with Erefm, indicating that the IPC can accurately calculate the PCC voltage and compensate it to the RPCL, which effectively reduces the adjustment time. The DPAC visibly contributes to power angle compensation as |δc2|>>|δc1| during the symmetrical grid faults. In the PACSE mode, δc gradually decreases to the threshold δth.

Figure 16 provides waveforms of PCC voltage, inverter-side current, output power, and switch with the proposed control strategy (Vg=0.2Vgn). The proposed control strategy demonstrates remarkable performance in suppressing transient currents and achieving fast power tracking.

Fig. 16  Waveforms of PCC voltage, inverter-side current, output power, and switch with proposed control strategy (Vg=0.2Vgn).

Furthermore, the calculation of the power angle compensation δc2 is based on the PCC voltage. The simulation results shown in Figs. 14-16 indicate that the proposed control strategy can quickly stabilize the PCC voltage, thereby ensuring the calculation of δc2.

The simulation results in Figs. 14-16 demonstrate the effectiveness of the proposed control strategy, which directly compensates the power angle and internal potential to reduce the regulatory effect of the power loop. Subsequently, Figs. 17 and 18 feature simulation results under various grid impedances and grid voltage sags. The proposed control strategy exhibits robust adaptability, ensuring that both output power and current comply with the grid codes for LVRT.

Fig. 17  Current and power waveforms under different grid impedances during LVRT. (a) Lg=9 mH. (b) Lg=15 mH. (c) Lg=20 mH.

Fig. 18  Current and power waveforms under different grid voltage sags during LVRT. (a) Vg=0.8Vgn. (b) Vg=0.4Vgn. (c) Vg=0.1Vgn.

The waveforms in Figs. 14-18 show a large overshoot and slow recovery process. The main reasons are as follows: ① since the APCL continues to operate in the voltage recovery mode, its inertia characteristic leads to power overshoot and prolonged recovery time; and ② the power angle compensation term remains even after the end of the voltage recovery mode and cannot be cleared until at least the time Tset elapses.

Considering that changing the control structure or parameters of the APCL requires additional switches or algorithms, which would further increase the control complexity, the recovery performance is improved by adjusting the preset time Tset. The simulation results are illustrated in Supplementary Material C Fig. SC1. It is evident that reducing Tset allows the power angle compensation term to exit quickly, reducing both overshoot and recovery time of active power.

It should be noted that the transient synchronization stability is also an important index for the VSG. After adopting the power-oriented reference generation method described in (5)-(7), as the PCC voltage magnitude VCm decreases, the active power reference Pref also decreases. This contributes to improving the transient stability of the VSG [

24]. Subsequently, to further demonstrate that the power-oriented reference generation method can ensure the transient synchronous stability of the VSG, the phase portraits of the VSG during different symmetrical grid faults under different grid inductance are plotted through simulation, as shown in Fig. 19, where point A is the initial operating point. It should be noted that the transient stability problems caused by the grid faults often occur under weak grid conditions. Thereupon, a larger grid-side inductance is selected, resulting in a lower short-circuit ratio (SCR) [24].

Fig. 19  Phase portraits of VSG during different symmetrical grid faults under different grid inductance. (a) Lg=20 mH and SCR=2.3. (b) Lg=30 mH and SCR=1.5.

In Fig. 19, even when the grid voltage sags to 0.1Vgn under weak grid conditions, the operating point starting from point A can still reach a steady state, i.e., Δω=0 and δ remains constant. This indicates that the proposed control strategy can effectively ensure the transient synchronous stability of VSG.

C. Comparisons with Different LVRT Control Strategies

The control strategies of the VSG during LVRT using internal potential and phase angle compensation are proposed in [

14], [15]. The innovation of this paper is to propose a more precise compensation formula, as shown in (13) and (18). To demonstrate the innovation, comparative simulations are depicted in Fig. 20, where subscripts 1, 2, and 3 are the control strategies in [14], [15], and this paper, respectively. Moreover, to ensure a fair comparison, the parameters of the power control loops and VCCL in the simulations for [14] and [15] are kept consistent with those used in this paper.

Fig. 20  Waveforms of current magnitude, output power, and RPCL output voltage with different LVRT control strategies.

These parameters are well tuned and the parameters involved in the strategies proposed in [

14] and [15] are also adjusted by the recommended strategies.

In [

14], compensation is only done when Ecom=VCm. However, due to the reactive power requirements during the grid faults and the effects of VI, VCm does not match Edref at the steady state, which is also validated in Fig. 15. Therefore, the RPCL still requires a longer time to reach the steady state.

In [

15], a control strategy involving both IPC and power angle compensation is employed.

The circuit and compensation formula are shown in Fig. 21 and (SC1) in Supplementary Material C, respectively, where igabc is the grid-side current, L2 is the grid-side inductance of the LCL filter, and vLabc is the output voltage of the LCL filter.

Fig. 21  Circuit in [

15].

The IPC in [

15] primarily serves to limit the current rather than directly estimate Edref in the steady state. As a result, the control strategy in [15] does not significantly enhance the response speed of the RPCL compared with the control strategy in [14]. When calculating δcd0, it is crucial to use the capacitor voltage components of the dq coordinates before the grid faults, as illustrated in Fig. 7. However, [15] does not change the coordinates during the grid faults, leading to inaccurate power angle compensation. This prevents the active power from swiftly tracking the setpoint. Additionally, [15] requires an additional set of voltage sensors to measure the output voltage of the LCL filter, which increases the equipment costs and overall complexity of the system setup.

In contrast to earlier control strategies shown in Fig. 20, the internal potential and power angle compensation formulas proposed in this paper offer enhanced accuracy without the need for additional voltage sensors. This refinement allows the output power to swiftly track the setpoint during the grid faults and maintain robust dynamic characteristics in the voltage recovery mode. Furthermore, the incorporation of TVI effectively reduces the transient current.

D. Simulation Results During Asymmetrical Grid Faults

The proposed control strategy also tackles the LVRT challenges presented by asymmetrical grid faults through the implementation of positive- and negative-sequence separation controls.

V. Experimental Validation

To prove the effectiveness of the proposed control strategy, experimental validation is carried out. The experimental parameters are shown in Table III. The experimental hardware platform is shown in Fig. 23. The grid impedance is composed of series inductors. In addition, the grid simulator is used to achieve the grid faults.

TABLE III  Experimental Parameters
ParameterSymbolValues
Normal grid voltage magnitude Vgn 100 V
Rated internal potential magnitude UCn 100 V
DC-link voltage Vdc 250 V
Virtual moment of inertia J 0.005 kg·m2
Damping coefficient Dp 0.76 N·m·s
Inertia coefficient K 3.3 A·s
Current threshold Ith 12 A
Nominal grid current amplitude IN 10 A
Power thresholds Pth, Qth 100 W, 100 var

Fig. 22  Waveforms of inverter-side current and output power during asymmetrical grid faults.

Fig. 23  Experimental hardware platform.

Figure 24 presents the experimental results with conventional VSG control during grid fault (Vg=0.2Vgn).

Fig. 24  Experimental results with conventional VSG control during grid faults (Vg=0.2Vgn). (a) Grid voltage and inverter-side current. (b) PCC power.

The VSG exhibits slow response within the power loop, coupled with considerable transient current. Conversely, Fig. 25 illustrates the experimental results with proposed control strategy during grid faults (Vg=0.2Vgn). The proposed control strategy enables swift stabilization of the VSG and effective transient current suppression.

Fig. 25  Experimental results with proposed control strategy during grid faults (Vg=0.2Vgn). (a) Grid voltage and inverter-side current. (b) PCC power, IPC, and power angle compensation.

The above results validate the theoretical derivation and confirm the effectiveness of the proposed control strategy for improving the LVRT performance of VSG.

VI. Conclusion

This paper examines the variations in power angle and PCC voltage of VSG during the grid faults and introduces a compound compensation control strategy. The conclusions of this paper are as follows:

1) The grid faults cause a reduction in the cutoff frequency of the power loop, thus leading to extended periods of power regulation and overcurrent.

2) The application of the internal potential and power angle compensation, which are derived from changes in the PCC voltage, active and reactive power references, and power angle, allows for rapid compensation of the power loop to its value in the steady state. This ensures that the VSG swiftly tracks the power references during grid faults.

3) Eigenvalue analysis and simulation studies confirm that the proposed control strategy effectively mitigates transient oscillations in the VSG during LVRT. The proposed control strategy demonstrates robust adaptability under different conditions. Furthermore, the large-signal model analysis reveals that TVI is instrumental in suppressing transient disturbances.

References

1

J. Fang, H. Deng, and S. M. Goetz, “Grid impedance estimation through grid-forming power converters,” IEEE Transactions on Power Electronics, vol. 36, no. 2, pp. 2094-2104, Feb. 2021. [Baidu Scholar] 

2

R. Zhang, C. Zhang, X. Xing et al., “Clustered voltage control and modulation scheme with switching loss reduction for high power hybrid cascaded converter,” IEEE Transactions on Industrial Electronics, vol. 71, no. 5, pp. 4911-4921, May 2024. [Baidu Scholar] 

3

C. Yang, L. Huang, H. Xin et al., “Placing grid-forming converters to enhance small signal stability of PLL-integrated power systems,” IEEE Transactions on Power Systems, vol. 36, no. 4, pp. 3563-3573, Jul. 2021. [Baidu Scholar] 

4

Y. Li, C. Lin, J. Hu et al., “PLL synchronization stability of grid-connected VSCs under asymmetric AC faults,” IEEE Transactions on Energy Conversion, vol. 37, no. 4, pp. 2438-2448, Dec. 2022. [Baidu Scholar] 

5

N. B. Lai, A. Tarrasó, G. N. Baltas et al., “External inertia emulation controller for grid-following power converter,” IEEE Transactions on Industry Applications, vol. 57, no. 6, pp. 6568-6576, Nov. 2021. [Baidu Scholar] 

6

H. Wu, X. Ruan, D. Yang et al., “Small-signal modeling and parameters design for virtual synchronous generators,” IEEE Transactions on Industrial Electronics, vol. 63, no. 7, pp. 4292-4303, Jul. 2016. [Baidu Scholar] 

7

B. Pournazarian, R. Sangrody, M. Lehtonen et al., “Simultaneous optimization of virtual synchronous generators parameters and virtual impedances in islanded microgrids,” IEEE Transactions on Smart Grid, vol. 13, no. 6, pp. 4202-4217, Nov. 2022. [Baidu Scholar] 

8

N. Mohammed, M. H. Ravanji, W. Zhou et al., “Online grid impedance estimation-based adaptive control of virtual synchronous generators considering strong and weak grid conditions,” IEEE Transactions on Sustainable Energy, vol. 14, no. 1, pp. 673-687, Jan. 2023. [Baidu Scholar] 

9

A. Pal, D. Pal, and B. K. Panigrahi, “A current saturation strategy for enhancing the low voltage ride-through capability of grid-forming inverters,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 70, no. 3, pp. 1009-1013, Mar. 2023. [Baidu Scholar] 

10

K. O. Oureilidis and C. S. Demoulias, “A fault clearing method in converter-dominated microgrids with conventional protection means,” IEEE Transactions on Power Electronics, vol. 31, no. 6, pp. 4628-4640, Jun. 2016. [Baidu Scholar] 

11

K. Shi, W. Song, P. Xu et al., “Low-voltage ride-through control strategy for a virtual synchronous generator based on smooth switching,” IEEE Access, vol. 6, pp. 2703-2711, Dec. 2017. [Baidu Scholar] 

12

T. Qoria, H. Wu, X. Wang et al., “Variable virtual impedance-based overcurrent protection for grid-forming inverters: small-signal, large-signal analysis and improvement,” IEEE Transactions on Smart Grid, vol. 14, no. 5, pp. 3324-3336, Sept. 2023. [Baidu Scholar] 

13

A. D. Paquette and D. M. Divan, “Virtual impedance current limiting for inverters in microgrids with synchronous generators,” IEEE Transactions on Industry Applications, vol. 51, no. 2, pp. 1630-1638, Mar. 2015. [Baidu Scholar] 

14

K. Koiwa, K. Inoo, T. Zanma et al., “Virtual voltage control of VSG for overcurrent suppression under symmetrical and asymmetrical voltage dips,” IEEE Transactions on Industrial Electronics, vol. 69, no. 11, pp. 11177-11186, Nov. 2022. [Baidu Scholar] 

15

R. Luo, H. Gao, F. Peng et al., “An adaptive low-voltage ride-through method for virtual synchronous generators during asymmetrical grid faults,” IEEE Transactions on Sustainable Energy, vol. 15, no. 3, pp. 1589-1600, Jul. 2024. [Baidu Scholar] 

16

Y. Jiang, C. Qin, X. Xing et al., “A hybrid passivity-based control strategy for three-level T-type inverter in LVRT operation,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 8, no. 4, pp. 4009-4024, Dec. 2020. [Baidu Scholar] 

17

A. Calle-Prado, S. Alepuz, J. Bordonau et al., “Model predictive current control of grid-connected neutral-point-clamped converters to meet low-voltage ride-through requirements,” IEEE Transactions on Industrial Electronics, vol. 62, no. 3, pp. 1503-1514, Mar. 2015. [Baidu Scholar] 

18

D. Shin, K. J. Lee, J. P. Lee et al., “Implementation of fault ride-through techniques of grid-connected inverter for distributed energy resources with adaptive low-pass notch PLL,” IEEE Transactions on Power Electronics, vol. 30, no. 5, pp. 2859-2871, May 2015. [Baidu Scholar] 

19

W. Wu, Y. Chen, L. Zhou et al., “Sequence impedance modeling and stability comparative analysis of voltage-controlled VSGs and current-controlled VSGs,” IEEE Transactions on Industrial Electronics, vol. 66, no. 8, pp. 6460-6472, Aug. 2019. [Baidu Scholar] 

20

Technical Requirements for Power Conversion System of Electrochemical Energy Storage System, Chinese National Standard GB/T 34120-2023, Dec. 2023. [Baidu Scholar] 

21

IEEE Standard for Interconnection and Interoperability of Inverter-Based Resources (IBRs) Interconnecting with Associated Transmission Electric Power Systems, IEEE Std 2800-2022, Apr. 2022. [Baidu Scholar] 

22

H. Wu and X. Wang, “Small-signal modeling and controller parameters tuning of grid-forming VSCs with adaptive virtual impedance-based current limitation,” IEEE Transactions on Power Electronics, vol. 37, no. 6, pp. 7185-7199, Jun. 2022. [Baidu Scholar] 

23

X. Xiong, C. Wu, B. Hu et al., “Transient damping method for improving the synchronization stability of virtual synchronous generators,” IEEE Transactions on Power Electronics, vol. 36, no. 7, pp. 7820-7831, Jul. 2021. [Baidu Scholar] 

24

X. Xiong, C. Wu, and F. Blaabjerg, “Effects of virtual resistance on transient stability of virtual synchronous generators under grid voltage sag,” IEEE Transactions on Industrial Electronics, vol. 69, no. 5, pp. 4754-4764, May 2022. [Baidu Scholar] 

25

Y. Yang, F. Blaabjerg, and Z. Zou, “Benchmarking of grid fault modes in single-phase grid-connected photovoltaic systems,” IEEE Transactions on Industry Applications, vol. 49, no. 5, pp. 2167-2176, Sept. 2013. [Baidu Scholar] 

26

M. G. Taul, X. Wang, P. Davari et al., “Current limiting control with enhanced dynamics of grid-forming converters during fault conditions,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 8, no. 2, pp. 1062-1073, Jun. 2020. [Baidu Scholar] 

27

M. Chen, D. Zhou, and F. Blaabjerg, “Enhanced transient angle stability control of grid-forming converter based on virtual synchronous generator,” IEEE Transactions on Industrial Electronics, vol. 69, no. 9, pp. 9133-9144, Sept. 2022. [Baidu Scholar] 

28

T. Qoria, Q. Cossart, X. Guillaud et al. (2019, Dec.). WP3-control and operation of a grid with 100% converter-based devices. [Online]. Availbale: https://www.researchgate.net/publication/338421895MIGRATE_D36_Requirement_guidelines_for_operating_a_grid_with_100_power_electronic_devices [Baidu Scholar]