Journal of Modern Power Systems and Clean Energy

ISSN 2196-5625 CN 32-1884/TK

网刊加载中。。。

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

确定继续浏览么?

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

Impedance Modeling and Stability Factor Assessment of Grid-connected Converters Based on Linear Active Disturbance Rejection Control  PDF

  • Hailiang Xu (Member, IEEE)
  • Fei Nie
  • Zhongxing Wang
  • Shinan Wang
  • Jiabing Hu (Senior Member, IEEE)
China University of Petroleum (East China), Qingdao, China; Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China

Updated:2021-11-23

DOI:10.35833/MPCE.2021.000280

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

Abstract

With the increase of converter-based renewable energy generation connected into the power grid, the interaction between renewable energy and grid impedance has introduced lots of new issues, among which the sub- and super-synchronous oscillation phenomenon makes a big concern. The linear active disturbance rejection control (LADRC) is a potential way to improve the damping characteristics of the grid-connected system, but the key factors and influencing mechanism on system stability are unknown. This paper establishes the equivalent impedance and coupling admittance models of a typical three-phase grid-connected converter. Then, the influence of the key factors such as the bandwidth of the LADRC and grid impedance on the stability and frequency coupling effect is assessed in detail. Finally, the theoretical analysis results are verified by simulations and experiments.

I. Introduction

NOWADAYS, more and more renewable energy generations are being connected into the power grid, where the grid-connected converters play an irreplaceable role between the renewable energy and the power grid [

1]. It should be noted that large-scale renewable energy generations such as wind and photovoltaic farms are usually located in remote areas. Hence, the long distance transmission and multistage transformers make the grid impedance in-negligible in stability analysis [2], [3]. Moreover, with the increase of grid impedance, the grid short-circuit ratio (SCR) decreases gradually, which will impair the current quality of grid-connected renewable energy [4].

The interaction between the renewable energy generation and weak grid may lead to serious stability issues. For instance, in Hami wind farm located in Xinjiang Uygur Autonomous Region of China, a sub-synchronous oscillation accident caused by the interaction between the wind converters and the weak grid occurred in 2015 [

5]. Note that such oscillation phenomenon has a common feature, i.e., a pair of interdependent sub- and super-synchronous oscillation components will appear synchronously. It is usually called the frequency coupling (FC) [6], or mirror frequency (MF) effect [7].

It is essential to point out that this kind of oscillation is quite different from that in the conventional single input and single output (SISO) system [

8], [9]. Due to the asymmetrical control links such as phase-locked loop (PLL) [10], and the considerable grid impedance [11], the current response will contain a disturbance component with its frequency at 2f1-fp or fp-2f1 in addition to the original frequency at fp. Thus, the system becomes a single input and double output (SIDO) system, rather than the SISO one [12].

Consequently, the original stability criterion for the SISO system is no longer applicable. Instead, a generalized Nyquist criterion (GNC) can be utilized [

13], [14]. Due to the complex analysis of the GNC [15], the parameters of the system cannot be conveniently modified according to the stability margin. In [7], an improved sequence impedance modeling method is proposed, which can greatly simplify the modeling process. The method is also valid when the frequency coupling effect is considered [16], [17]. In the dq frame, a method to transform the impedance model from a multi-input multi-output (MIMO) system into an SISO one is proposed in [9]. This transformation can greatly simplify the stability analysis [18].

Various methods have been presented to suppress the oscillation components [

8], [19], [20]. In [19], an asymmetrical current control method is proposed to make up for the asymmetrical characteristics introduced by PLL. In [8], a flexible symmetrical control structure is designed by introducing the concept of complex phase angle. Theoretical analysis and experimental results demonstrate that the symmetrical PLL can significantly suppress the frequency coupling effect. Nevertheless, this method inevitably needs to alter the internal structure of PLL [21], and makes the stability analysis process more complicated.

As an alternative, the linear active disturbance rejection control (LADRC) is applied into the PLL in [

22], which can suppress the effect of FC caused by the grid impedance and unsymmetrical control loops. The original concept of the LADRC can be traced back to the ADRC in [23]. And LADRC is proposed in [24] by linearizing the controller and extended state observer (ESO), which makes the algorithm simple and easy to implement. Nowadays, the LADRC has been widely used in engineering because of its unique disturbance-rejection characteristic. Especially, in the grid-connected converter system, LADRC shows its superior control performance. In [25], the LADRC is introduced into the DC-link voltage controller to replace the conventional proportional-integral (PI) regulator for variable speed hydro-electric plants, and the transient performance of the LADRC-based system with a stiff grid is satisfactory. In [26] and [27], the LADRC is applied to the DC-link voltage control loop, which enhances the robustness of the grid-connected system. In addition, an improved LESO is proposed to form a new generalized integrator-ESO (GI-ESO) in PLL [28]. This new type of LADRC has an excellent performance in suppressing the background harmonics of the grid voltage.

Although the LADRC has been applied in converter-based grid-connected systems, the stability issue has not yet been sufficiently investigated, especially when the grid impedance is considerable. In addition, the LADRC is a potential way to improve the damping characteristics of grid-connected systems, but the key factors and the influencing mechanism on stability are unknown.

To address this issue, this paper presents a comprehensive stability assessment of LADRC-based converter control system with weak grid connection. And the single-parameter LADRC is adopted in the grid-connected system, considering that it cannot only effectively reduce the number of parameters to be tuned, but also has a ideal phase margin, compared with the traditional PI control [

29], [30]. As a specific contribution, this paper establishes the equivalent impedance and coupling admittance models of a typical three-phase pulse width modulation (PWM) converter under weak grid condition, where the LADRC is utilized into the DC-link voltage control. Furthermore, the system stability is analyzed according to frequency characteristic curves and Nyquist curves. Moreover, the influence of the key factors on system stability and frequency coupling effect is assessed in detail. Simulations and experiments are finally carried out to validate the correctness of the theoretical analysis.

II. LADRC-based DC-link Voltage Control of Grid-connected Converter

A. Configuration of Three-phase Grid-connected Converter

Figure 1 shows the main circuit and control structure of a three-phase grid-connected converter. Note that the LADRC is implemented in the DC-link voltage control loop, while the d- and q-axis current control and PLL still use the traditional PI regulator. In Fig. 1, ud and uq are the dq-axis components of PCC voltage; idref and iqref are the reference currents of dq-axis components, respectively; mi is the PWM signal; md and mq are the dq-axis components of mi; Q1-6 is the output signal of space vector pulse width modulation (SVPWM); uPCC is the point of common coupling (PCC) voltage; ig is the grid-connected current; Udc and Udcref are the DC-link voltage and its reference, respectively; θPLL is the output angle of PLL; HPLL and Hi are the PI controller of PLL and inner current loop, respectively; Cdc is the DC-link capacitor; Rload is the equivalent resistance of the DC-link load; and Lf and Lg are the equivalent inductances of filter and the power grid, respectively. For simplification, the grid impedance is expressed as Zg=jω0Lg, and ω0 is the fundamental frequency.

Fig. 1 Main circuit and control structure of three-phase grid-connected converter with LADRC.

As for the DC-link circuit, it can be given as:

CdcdUdcdt=idc-iload (1)

where idc and iload are the total DC-link current and the equivalent load current, respectively.

In order to obtain the parameters in LADRC, the relation between Udc and idref needs to be determined firstly. Note that Udc is a dynamic quantity, while idref is a steady-state one. Thus, it is necessary to make the following assumptions: ① the power losses on the converters and the internal resistor of the inductive filter can be neglected, i.e., PACPDC; ② the DC-link voltage and grid-connected current are well controlled, i.e., UdcUdcref and ididref [

31].

PAC=32U1id32U1idrefPDC=UdcidcUdcrefidcidc3U12Udcrefidref (2)

where U1 is the fundamental magnitude of grid voltage.

Based on (1) and (2), the relation between idref and Udc can be expressed by:

dUdcdt=3U12CdcUdcrefidref-UdcCdcRload (3)

It is notable that this is a transitional expression that can be used to derive the control parameters of LADRC.

B. Design of Second-order LADRC

Suppose that a second-order LADRC is applied in the DC-link voltage loop. Then, its model can be represented by:

y¨=-a1y˙-a2y+w+(b-b0)uftotal+b0u (4)

where u, y, and w are the input, output and unknown external disturbances, respectively; a1 and a2 are the parameters of the unknown system; b is the unknown input gain, while b0 is a known value; and ftotal is the total disturbance, containing both the internal and external uncertainties.

Assume that y=Udc and u=idref. By taking the derivation of (3), the expression of b0 can be obtained, as show in (5).

d2Udcdt2=-3U12RloadCdc2Udcrefb0idref+Udc(CdcRload)2ftotal (5)

Introduce state variables x1, x2, and extended state variable x3, with x1=y, x2=y˙, and x3=ftotal. Then, the LESO can be organized as:

z˙1=z2+k1(y-z1)z˙2=z3+k2(y-z1)+b0uz˙3=k3(y-z1) (6)

where z1, z2, and z3 are the estimations of x1, x2, and ftotal, respectively; and k1, k2 and k3 are the observer gains, respectively.

If the estimation error between z3 and ftotal can be neglected, i.e., z3x3=ftotal, the original nonlinear control system becomes a linear integrated series control system. Therefore, the linear state error feedback (LESF) law can be designed as:

u=1b0[l2(r-z1)-l1z2-z3] (7)

where r is the input reference.

Finally, the expression of LESO (6) and the expression of LSEF law (7) make up the whole LADRC system.

Note that there are five parameters to be confirmed, i.e., l1, l2, k1, k2, and k3. In [

29], it reveals that the single-parameter structure enables the system to obtain a fast step-free response without overshoot, and it is superior to the traditional PI controller in suppressing low-frequency external and high-frequency internal disturbances. Hence, a parameter ωL called bandwidth of the LADRC is introduced. Consequently, l1, l2, k1, k2, and k3 can then be expressed by ωL as:

l1=2ωLl2=ωL2k1=3ωLk2=3ωL2k3=ωL3 (8)

It is worth noting that ωL corresponds to the control bandwidth of the traditional PI controller, and is different from the controller bandwidth ωc, which is usually 3-10 times the observer bandwidth ωo in the traditional LADRC. In other words, ωL=ωc=ωo. Consequently, this design can reduce the number of control parameters, which further simplifies the tuning process of parameters. In addition, it can also make the system have a satisfactory phase margin [

30].

C. Analysis of Tracking and Anti-interference Characteristics

Taking the design of d-axis current controller as an example, and considering the signal-sampling delay of the inner current loop, the relationship between idref and Udc is shown in Fig. 2, where kpi and kii are the control parameters of inner current loop, respectively.

Fig. 2 Control structure of inner current loop.

According to Fig. 2, the transfer function between idref and Udc can be deduced as:

Gp(s)=3U1Rload2Udcref(RloadCdcs+1)(4Tss+1) (9)

where Ts is the sampling time.

Figure 3 shows the DC-link voltage control loop, where the LADRC is adopted. Assume U(s), Y(s), and R(s) are the Laplace transforms of u, y, and r, respectively. For a clear description, define C1(s) and C(s) to be the respective feedforward and feedback controllers, with their transfer functions shown in (10). Consequently, Fig. 3 can be simplified to the standard structure of LADRC, as shown in Fig. 4.

Fig. 3 DC-link voltage control loop.

Fig. 4 Simplified control structure of LADRC.

C1(s)=ωL2(s+ωL)3ωL3(10s2+5ωLs+ωL2)C(s)=ωL3(10s2+5ωLs+ωL2)b0s(s2+5ωLs+10ωL2)=5ωLb0+ωL2b0s+10b0sCLADRC_PID(s)1s2+5ωLs+10ωL2CLADRC_LPF(s) (10)

where CLADRC_PID(s) and CLADRC_LPF(s) are the equivalent proportional integral differential (PID) and second-order filter parts in C(s), respectively.

The frequency response characteristics of the open-loop transfer function are shown in Fig. 5(a). It can be observed that the two parts of the LADRC, i.e., CLADRC_PID(s) and CLADRC_LPF(s), are indivisible. It means that only when the two parts get combined, we can obtain the satisfactory bandwidth and phase margin.

Fig. 5 Frequency characteristic curves of L(s). (a) LLADRC(s)=C(s)Gp(s), LLADRC_PID(s) = CLADRC_PID(s)Gp(s),   and   LLADRC_LPF(s) = CLADRC_LPF(s)Gp(s). (b)LLADRC(s)=C(s)Gp(s),  LPI(s)=HPI(s)Gp(s),  and  LLPF(s)=HLPF(s)Gp(s).

To make it more clear, Fig. 5(b) compares bode diagrams of L(s). Note that the bandwidths of these three controllers are set to be the same. With the same control bandwidth, the system phase margin with LADRC is much higher than that of the PI or LPF, implying that the system stability can be evidently enhanced when the LADRC is adopted.

Then, according to Fig. 4, the close-loop transfer function of the LADRC-based system can be derived as:

y=C1(s)C(s)GP(s)1+C(s)GP(s)r+GP(s)1+C(s)GP(s)ftotal (11)

Figure 6 depicts the frequency characteristic curves with different ωL. As shown in Fig. 6(a), when ωL increases from 300 rad/s to 700 rad/s, the control bandwidth of the tracking signal gets enlarged. However, the attenuation ability of the high frequency signal will be weakened. To ensure enough tracking and anti-interference performance, ωL should be 5-10 times of the bandwidth of traditional PI controller [

25].

Fig. 6 Frequency characteristic curves with different ωL. (a) Tracking characteristic. (b) Anti-interference characteristic.

Figure 6(b) shows the anti-interference characteristics of the LADRC as ωL increases. It is notable that with the increase of ωL, the anti-interference characteristic on the middle and low frequency band will be influenced, while the high frequency band is not effected significantly. However, because of the high frequency roll-off, the high-frequency internal disturbance will be suppressed as well.

III. Impedance Modeling of Three-phase Grid-connected Converter with LADRC

A. Equivalent Impedance Modeling of System with LADRC

When the grid impedance Lg is considered, the coupling current disturbance Iωp-2ω0 is superimposed with the original one, which intensifies the frequency coupling effect of the system. This disturbance transfer process can be represented as the signal flow diagram shown in Fig. 7.

Fig. 7 Disturbance transfer process considering Lg.

In Fig. 7, Ypp(s) and Ypn(s) represent the equivalent admittance from current disturbance Iωp and coupling current disturbance Iωp-2ω0 to voltage disturbance Uωp, and Ypp(s) and Ypn(s) can be deduced as (12) and (13).

Ypp(s)=M1+H2,ωpI1*+H2,ωp-2ω0I1-[(s-jω0)Cdc+1/Rload]/(3kPWM)+H4,ωpI1*+H4,ωp-2ω0I1-H3,ωpI1*-H3,ωp-2ω0I1-kPWMUdcrefH2,ωp-kPWMUdcrefH4,ωp-kPWMM1-1+kPWMUdcrefH3,ωp-kPWMUdcrefH2,ωp-2ω0-(s-j2ω0)Lf-kPWMUdcrefH4,ωp-2ω0-kPWMM1*kPWMUdcrefH3,ωp-2ω0M1*+H1,ωpI1*+H1,ωp-2ω0I1M1+H2,ωpI1*+H2,ωp-2ω0I1-[(s-jω0)Cdc+1/Rload]/(3kPWM)+H4,ωpI1*+H4,ωp-2ω0I1-kPWMUdcrefH1,ωp-sLf-kPWMUdcrefH2,ωp-kPWMUdcrefH4,ωp-kPWMM1-kPWMUdcrefH1,ωp-2ω0-kPWMUdcrefH2,ωp-2ω0-(s-j2ω0)Lf-kPWMUdcrefH4,ωp-2ω0-kPWMM1* (12)
Ypn(s)=M1*+H1,ωpI1*+H1,ωp-2ω0I1-[(s-jω0)Cdc+1/Rload]/(3kPWM)+H4,ωpI1*+H4,ωp-2ω0I1-H3,ωpI1*-H3,ωp-2ω0I1-kPWMUdcrefH1,ωp-sLf-kPWMUdcrefH4,ωp-kPWMM1-1+kPWMUdcrefH3,ωp-kPWMUdcrefH1,ωp-2ω0-kPWMUdcrefH4,ωp-2ω0-kPWMM1*kPWMUdcrefH3,ωp-2ω0M1*+H1,ωpI1*+H1,ωp-2ω0I1M1+H2,ωpI1*+H2,ωp-2ω0I1-[(s-jω0)Cdc+1/Rload]/(3kPWM)+H4,ωpI1*+H4,ωp-2ω0I1-kPWMUdcrefH1,ωp-sLf-kPWMUdcrefH2,ωp-kPWMUdcrefH4,ωp-kPWMM1-kPWMUdcrefH1,ωp-2ω0-kPWMUdcrefH2,ωp-2ω0-(s-j2ω0)Lf-kPWMUdcrefH4,ωp-2ω0-kPWMM1* (13)

where Hi,ωp and Hi,ωp-2ω0 (i=1, 2, 3, 4) are the influence coefficients of each disturbance component on the PWM waveforms, respectively; M1=(U1-sLfI1)/(kPWMUdcref) is the fundamental component of PWM waveform; U1 and I1 are the fundamental voltage and current, respectively; Uωp-2ω0 is the coupling voltage disturbance; ‘*’ represents a conjugation of the variable.

Note that Ynn(s) and Ynp(s) can be deduced by Ypp(s) and Ypn(s) [

32], [33], i.e.,

Ynn(s)=Ypp*(s-j2ω0)Ynp(s)=Ypn*(s-j2ω0) (14)

Therefore, the equivalent impedance of the grid-side converter considering FC under the weak grid condition can be written as:

Zaceq(s)=UωpIωp=1Ypp(s)-Ypn(s)Ynp(s)Zg(s-j2ω0)1+Ynn(s)Zg(s-j2ω0)Yc(s) (15)

where Yc(s) is the equivalent admittance of reactive frequency-coupling effect in Fig. 7.

B. Equivalent Impedance Modeling of Control Links

The disturbance component of DC-link voltage in frequency domain can be expressed as:

Udc,ωp-ω0=3kPWMMωpI1*+Mωp-2ω0I1+M1*Iωp+M1Iωp-2ω0(s-jω0)Cdc+1/Rload (16)
Map(f)=Mωp=H1,ωpsIωp+H2,ωpsIωp-2ω0+H3,ωpsUωp+H4,ωpsUdc,ωpf=±fpMωp-2ω0=H1,ωp-2ω0sIωp+H2,ωp-2ω0sIωp-2ω0+H3,ωp-2ω0sUωp+H4,ωp-2ω0sUdc,ωpf=±(fp-2f1) (17)

where Mωp and Mωp-2ω0 are the disturbances of PWM waveform at frequency of fp and fp-2f1, respectively.

As a result, a disturbance component with frequency fp-f0 will be contained in the reference of inner current loop or the DC-link output of the voltage loop controller, i.e.,

Idref,ωp-ω0(f)=Udc,ωp-ω0(f)C(sjω0)f=±fp (18)

In addition, θPLL will be affected by the existence of voltage disturbance component. Further, the d- and q-axis components of the grid-connected current will be influenced by the angle disturbance in the output of PLL. Considering the coordinate transformation links, the disturbance of output angle ΔθPLL [

16] can be expressed as:

ΔθPLL(f)=jTPLL(s)Uωpf=±(fp-f1) (19)
TPLL(s)=HPLL(s)1+U1HPLL(s) (20)

where HPLL(s)=(kpp+kip/s)/s is the forward channel transfer function of PLL, and kpp and kip are the control parameters of PLL.

Then, the d- and q-axis components of the grid-connected current affected by the disturbance of output angle can be deduced by:

Id(f)=Iωp+Iωp-2ω0I1(sin φi1)ΔθPLLf=±(fp-f1)Iq(f)=jIωp±jIωp-2ω0±I1(cos φi1)ΔθPLLf=±(fp-f1) (21)

where φi1 is the angle of fundamental current.

According to Fig. 1, the d- and q-axis components of the PWM modulated signal can be expressed as:

md(s)=[(Udcref-Udc)C(s)-id]Hi(s)-ω0Lfiqmq(s)=(iqref-iq)Hi(s)+ω0Lfid (22)

where Hi(s)=kpi+kii/s is the transfer function of current PI controller.

And the PWM waveform of phase a can be expressed as:

ma(s)=cos(θPLL(s))mds-sin(θPLL(s))mq(s) (23)

By considering (16)-(23), the coefficients in (12) and (13) can be solved and obtained as:

H1,ωp(s)=Hi(s-jω0)-jω0LfH1,ωp-2ω0(s)=0 (24)
H2,ωp(s)=0H2,ωp-2ω0(s)=Hi(s+jω0)+jω0Lf (25)
H3,ωp(s)=j12TPLL(s-jω0)Hi(s-jω0)C(s-jω0)+TPLL(s-jω0){[jω0Lf-Hi(s-jω0)]I1+M1}H3,ωp-2ω0(s)=j12TPLL(s+jω0)Hi(s+jω0)C(s+jω0)+TPLL(s+jω0){[jω0Lf+Hi(s+jω0)]I1*-M1*} (26)
H4,ωp(s)=12C(s-jω0)Hi(s-jω0)H4,ωp-2ω0(s)=12C(s+jω0)Hi(s+jω0) (27)

IV. Key Factor Assessment on System Stability

According to the impedance-based stability criterion, if the ratio between the grid impedance Zg and the inverter output impedance Zeq satisfies the Nyquist criterion, the system would be stable [

32]. This section will evaluate the influence of LADRC on the system stability.

A. Stability Mechanism Analysis

Figure 8 displays the bode diagrams of Zeq(s), Zg(s), and Ycouple(s) based on PI controller with different SCRs of the system. The damping characteristics of the system get deteriorated when the SCR decreases. As a consequence, the magnitude-frequency characteristic curves of Zeq(s) and Zg(s) intersect at several points, which are not conducive to the stability of the system. For instance, when SCR=2, the two curves intersect at four frequencies, i.e., 72 Hz, 80 Hz, 171 Hz, and 175 Hz, with the phase differences being 13°, 227°, 215°, and 95°, respectively. Note that, if the phase difference is more than 180°, the system will fall into an oscillation risk area.

Fig. 8 Frequency characteristic curves with PI controller (SCR = 2, 4, and 8). (a) Zeq(s) and Zg(s). (b) Ycouple(s).

In addition, as demonstrated in Fig. 8(b), the decrease of SCR increases the magnitude-frequency characteristic curve of Ycouple(s). It indicates that the intensity of the FC effect is negatively correlated with the SCR. Figure 9 presents the Nyquist curves with PI controller. With the decrease of SCR, the corresponding Nyquist curves gradually encircles (-1, j0). And when SCR=2, the system becomes unstable.

Fig. 9 Nyquist curves with PI controller. (a) SCR=2,  4,  and 8. (b) fv=20,  30 and 40 Hz.

In [

10], it has revealed that reducing the bandwidth of the control loop can enhance the stability and depress the FC effect of the system under weak grid, which is verified by the Nyquist curve shown in Fig. 9(b). Furthermore, the stability and FC degree of the system with low SCR are compared when the traditional PI controller, the bandwidth reduced PI (BR-PI) controller, and LADRC work, and comparison results are shown in Fig. 10.

Fig. 10 Frequency characteristic curves with various kinds of control strategy. (a) Zeq(s) and Zg(s). (b) Ycouple(s).

From Fig. 10, it can be observed that when LADRC is introduced, both the magnitude-frequency characteristic curves of Zeq(s) and Ycouple(s) show better performance than the other two controllers. Though reducing the bandwidth of PI controller can improve the stability of the system and reduce the degree of frequency coupling, it will worsen the dynamic performance of the system. Figure 10 also indicates that after adopting LADRC, the system stability can be enhanced with weak grid by improving the damping characteristics of the system.

B. Effect of Key Factors on Stability

In order to study the influence of LADRC parameters and grid impedance on the stability of the system, the frequency characteristic curves and Nyquist curves of Zeq(s) and Ycouple(s) are drawn when the values of the ωL and Lg are changed.

Figures 11 and 12 show the frequency characteristic curves with different Lg and ωL, while Fig. 13(a) and (b) presents the corresponding Nyquist curves. As shown in Fig. 11(a), there is almost no intersection between Zeq(s) and Zg(s). Even if at the narrow intersection area (when the SCR is further reduced to 1.5), the phase difference is still less than 180°. Similar to Fig. 8(b), it can be observed from Fig. 11(b) that the magnitude of the coupling admittance increases a little as the SCR decreases. However, the Nyquist curves do not surround (-1, j0) in these cases, as shown in Fig. 13(a). It indicates that after the introduction of the LADRC, the system can be stable with a fairly low SCR.

Fig. 11 Frequency characteristic curves with LADRC when SCR=1.5, 1.75, and 2. (a) Zeq(s) and Zg(s). (b) Ycouple(s).

Fig. 12 Frequency characteristic curves with LADRC when ωL = 100, 500, 700, and 1000 rad/s. (a) Zeq(s) and Zg(s). (b) Ycouple(s).

Fig. 13 Nyquist curves with LADRC. (a) SCR=1.5,  1.75,  and  2. (b) ωL=100,  500,  700, and  1000 Hz.

As shown in Fig. 12(a), the magnitude-frequency characteristic curve of Zeq(s) at low frequency band gets large as ωL increases, which means that the damping of the system can be enhanced.

Meanwhile, as shown in Fig. 12(b), the magnitude of the coupling admittance decreases when ωL increases, indicating that the FC effect is significantly suppressed. Furthermore, although the Nyquist curves in Fig. 13(b) gradually approximate (-1, j0) with the increase of ωL, the adjustable range is much larger than that of the traditional PI controller. Besides, compared with the results shown in Figs. 8 and 9, the grid-connected system can remain stable and have a satisfactory controllable bandwidth range with the LADRC under weak grid.

In summary, it can be concluded that the LADRC is able to enhance the stability and robustness of the system under weak grid. Furthermore, it can effectively suppress the FC effect caused by the grid impedance and asymmetrical control links. But the bandwidth of LADRC should not be too large, otherwise the contrary results may be obtained.

V. Simulation and Experimental Results

In this section, simulation and experiments are carried out to verify the stability analysis of the grid-connected converter system with LADRC. A 2 kW converter prototype is built in the lab, as shown in Fig. 14. Simulation and experimental parameters are shown in Table I and Table II.

Fig. 14 Hardware of experimental platform.

TABLE I Simulation System Parameters
ParameterValue
Udcref 650 V
U1 311 V
Lg 6.3, 3.2, 1.6 mH (SCR = 2, 4, 8)
Lf 3.5 mH
Cdc 4400 μF
Rload 20 Ω
kpp 0.367
kip 21.036
kpv 0.50, 0.75, 1.007
kiv 28.78, 64.76, 115.15
kpi 4.003
kii 2289
ωL 100, 300, 500, 700 rad/s
TABLE II Expertimental System Parameters
ParameterValue
Udcref 440 V
U1 179.59 V
fswitch, fsample 20, 10 kHz
Lg 18, 21, 24 mH (SCR=2.5,  2.1, and  1.9)
Lf 1.5 mH
Cdc 1230 μF
Rload 100 Ω
kpp 0.26
kip 2.3
kpv 0.04
kiv 22
kpi 27
kii 15
ωL 100, 300, 500, 700 rad/s

A. Stability Validation of LADRC

The simulation results with different control strategies are shown in Fig. 15(a). It demonstrates that when SCR decreases from 8 to 2, the oscillation of grid-connected current and DC-link voltage gradually intensifies. When SCR is reduced to 2, a pair of sub- and super-synchronous oscillation appears at 20 Hz and 80 Hz, respectively. This result corresponds to the theoretical analysis result shown in Fig. 8. Moreover, it can be observed from Fig. 15(b) that when LADRC is introduced, the sub- and super-synchronous oscillation can be suppressed well when SCR=2, and the total harmonic distortion (THD) value also reduces from 18.23% to 4.94% evidently.

Fig. 15 Simulation results with different control strategies. (a) Udc, Ig, and THD with PI controller. (b) Udc, Ig, and THD with LADRC. (c) Magnitude with PI controller. (d) Magnitude with LADRC.

Figure 16 gives the experimental results with different control strategies, it can be observed that the LADRC has a better performance than traditional PI controller in suppressing the sub- and super-synchronous oscillation in the grid-connected current. Furthermore, the system can remain stable when SCR is reduced to 1.9.

Fig. 16 Experimental results with different control strategies. (a) Udc and Ia with PI controller when SCR= 2.5. (b) Fast Fourier transform (FFT) result of Ia with PI controller when SCR= 2.5. (c) Udc and Ia with LADRC controller when SCR= 2.5. (d) FFT result of Ia with LADRC controller when SCR= 2.5. (e) Udc and Ia with PI controller when SCR= 2.1. (f) FFT result of Ia with PI controller when SCR= 2.1. (g) Udc and Ia with LADRC controller when SCR= 1.9. (h) FFT result of Ia with LADRC controller when SCR= 1.9.

B. Parameter Robustness and Dynamic Performance of LADRC

As mentioned above, changing the bandwidth of LADRC will have a certain impact on the system stability. As shown in Fig. 17, with the increase of ωL, the stability of the system will be enhanced to some extent. However, when ωL increases to 700 rad/s, high frequency ripples will emerge, which is consistent with the analysis in Section II. The ability to attenuate the high frequency disturbance will weaken with a too high bandwidth of ωL. Nevertheless, ωL can still change in a large range to obtain better robustness.

Fig. 17 Experimental results with different ωL. (a) ωL=100  rad/s. (b) ωL= 300  rad/s. (c) ωL=500  rad/s. (d) ωL=700  rad/s.

To test the dynamic performance of the system after the introduction of the LADRC, the value of the Udcref is suddenly changed from 440 V to 450 V when SCR=2.1. Figure 18 shows the experimental results with different ωL when Udcref increases from 440 V to 450 V. A bigger ωL can decrease the response time and there is no overshoot.

Fig. 18 Experimental results with different ωL when Udcref increases from 440 V to 450 V. (a) ωL=100  rad/s. (b) ωL=300  rad/s. (c) ωL=500  rad/s. (d) ωL=700  rad/s.

Figure 19 shows experimental results with different ωL when Iqref increases from 0 to 3 A. Similarly, the system remains stable when ωL changes in a wide range. It is worth noting that the transient response process is shortened obviously, and there is no overshoot, either. This is consistent with the situation when Udcref changes.

Fig. 19 Experimental results with different ωL when Iqref increases from 0 to 3 A. (a) ωL=100  rad/s. (b) ωL=300  rad/s. (c) ωL=500  rad/s. (d) ωL=700  rad/s.

In summary, the system stability is greatly improved by LADRC. Although the bandwidth of LADRC can be adjusted in a large range, attention should be paid to avoid the influence of system stability due to more ripple.

VI. Conclusion

This paper establishes the equivalent impedance and the coupling admittance models considering the influence of LADRC, inner current loop controller, PLL, and the grid impedance. Based on these models, the influences of several key factors are studied. And the specific conclusions can be summarized as follows.

1) Based on the analysis of the established impedance model, the damping characteristic of system at low SCR is significantly enhanced with the introduction of LADRC, and the FC effect can be greatly suppressed.

2) Only ωL needs to be adjusted during the tuning process of single-parameter structure LADRC. Moreover, the system can keep robust when ωL is changed in a wide range.

3) Though ωL can be tuned in a wide range, it is necessary to set the ωL selection at 2-5 times of the bandwidth of the traditional PI controller to guarantee a good dynamic performance. This can be used as a guide in the tuning of LADRC parameters.

4) Considering the satisfactory dynamic and static performance of the LADRC, it can also be used in other control loops of the grid-connected system, e.g., the PLL and the current control. Similarly, the stability issue needs to be evaluated carefully when adopting the LADRC.

References

1

M. N. Naz, S. Imtiaz, M. K. L. Bhatti et al., “Dynamic stability improvement of decentralized wind farms by effective distribution static compensator,” Journal of Modern Power Systems and Clean Energy, vol. 9, no. 3, pp. 516-525, May 2021. [Baidu Scholar

2

S. Song, Z. Wei, Y. Lin et al., “Impedance modeling and stability analysis of PV grid-connected inverter systems considering frequency coupling,” CSEE Journal of Power and Energy Systems, vol. 6, no. 2, pp. 279-290, Jun. 2020. [Baidu Scholar

3

X. Zhang, D. Xia, Z. Fu et al., “An improved feedforward control method considering PLL dynamics to improve weak grid stability of grid-connected inverters,” IEEE Transactions on Industry Applications, vol. 54, no. 5, pp. 5143-5151, Sept. 2018. [Baidu Scholar

4

E. Marrazi, G. Yang, and P. Weinreich-Jensen, “Allocation of synchronous condensers for restoration of system short-circuit power,” Journal of Modern Power Systems and Clean Energy, vol. 6, no. 1, pp. 17-26, Jan. 2018. [Baidu Scholar

5

X. Xie, Y. Zhan, J. Shair et al., “Identifying the source of subsynchronous control interaction via wide-area monitoring of sub/super-synchronous power flows,” IEEE Transactions on Power Delivery, vol. 35, no. 5, pp. 2177-2185, Oct. 2020. [Baidu Scholar

6

Y. Liao and X. Wang, “Stationary-frame complex-valued frequency-domain modeling of three-phase power converters,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 8, no. 2, pp. 1922-1933, Jun. 2020. [Baidu Scholar

7

A. Rygg, M. Molinas, C. Zhang et al., “A modified sequence-domain impedance definition and its equivalence to the dq-domain impedance definition for the stability analysis of AC power electronic systems,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 4, no. 4, pp. 1383-1396, Dec. 2016. [Baidu Scholar

8

D. Yang, X. Wang, F. Liu et al., “Symmetrical PLL for SISO impedance modeling and enhanced stability in weak grids,” IEEE Transactions on Power Electronics, vol. 35, no. 2, pp. 1473-1483, Feb. 2020. [Baidu Scholar

9

C. Zhang, X. Cai, A. Rygg et al., “Sequence domain SISO equivalent models of a grid-tied voltage source converter system for small-signal stability analysis,” IEEE Transactions on Energy Conversion, vol. 33, no. 2, pp. 741-749, Jun. 2018. [Baidu Scholar

10

Y. Xu, H. Nian, T. Wang et al., “Frequency coupling characteristic modeling and stability analysis of doubly fed induction generator,” IEEE Transactions on Energy Conversion, vol. 33, no. 3, pp. 1475-1486, Sept. 2018. [Baidu Scholar

11

L. Jia, X. Ruan, W. Zhao et al., “An adaptive active damper for improving the stability of grid-connected inverters under weak grid,” IEEE Transactions on Power Electronics, vol. 33, no. 11, pp. 9561-9574, Nov. 2018. [Baidu Scholar

12

W. Wu, L. Zhou, Y. Chen et al., “Sequence-impedance-based stability comparison between VSGs and traditional grid-connected inverters,” IEEE Transactions on Power Electronics, vol. 34, no. 1, pp. 46-52, Jan. 2019. [Baidu Scholar

13

Y. Li, Z. Shuai, X. Liu et al., “Stability analysis and location optimization method for multiconverter power systems based on nodal admittance matrix,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 9, no. 1, pp. 529-538, Feb. 2021. [Baidu Scholar

14

R. Wang, Q. Sun, D. Qin et al., “Steady-state stability assessment of AC-busbar plug-in electric vehicle charging station with photovoltaic,” Journal of Modern Power Systems and Clean Energy, vol. 8, no. 5, pp. 884-894, Sept. 2020. [Baidu Scholar

15

A. Emami-Naeini and R. L. Kosut, “The generalized Nyquist criterion and robustness margins with applications,” in Proceedings of 2012 IEEE 51st IEEE Conference on Decision and Control (CDC), Maui, USA, Feb. 2012, pp. 226-231. [Baidu Scholar

16

M. Cespedes and J. Sun, “Impedance modeling and analysis of grid-connected voltage-source converters,” IEEE Transactions on Power Electronics, vol. 29, no. 3, pp. 1254-1261, Mar. 2014. [Baidu Scholar

17

M. K. Bakhshizadeh, X. Wang, F. Blaabjerg et al., “Couplings in phase domain impedance modeling of grid-connected converters,” IEEE Transactions on Power Electronics, vol. 31, no. 10, pp. 6792-6796, Oct. 2016. [Baidu Scholar

18

J. Samanes, A. Urtasun, E. L. Barrios et al., “Control design and stability analysis of power converters: the MIMO generalized bode criterion,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 8, no. 2, pp. 1880-1893, Jun. 2020. [Baidu Scholar

19

X. Zhang, S. Fu, W. Chen et al., “A symmetrical control method for grid-connected converters to suppress the frequency coupling under weak grid conditions,” IEEE Transactions on Power Electronics, vol. 35, no. 12, pp. 13488-13499, Dec. 2020. [Baidu Scholar

20

B. Hu, H. Nian, M. Li et al., “Impedance characteristic analysis and reshaping method of DFIG system based on DPC without PLL,” IEEE Transactions on Industrial Electronics, vol. 68, no. 10, pp. 9767-9777, Oct. 2021. [Baidu Scholar

21

D. Zhu, S. Zhou, X. Zou et al., “Small-signal disturbance compensation control for LCL-type grid-connected converter in weak grid,” IEEE Transactions on Industry Applications, vol. 56, no. 3, pp. 2852-2861, May 2020. [Baidu Scholar

22

Z. Xie, Y. Chen, W. Wu et al., “Admittance modeling and stability analysis of grid-connected inverter with LADRC-PLL,” IEEE Transactions on Industrial Electronics, vol. 68, no. 12, pp. 12272-12284, Dec. 2021. [Baidu Scholar

23

J. Han, “From PID to active disturbance rejection control,” IEEE Transactions on Industrial Electronics, vol. 56, no. 3, pp. 900-906, Mar. 2009. [Baidu Scholar

24

Z. Gao, “Scaling and bandwidth-parameterization based controller tuning,” in Proceedings of the 2003 American Control Conference, Denver, USA, Nov. 2003, pp. 4989-4996. [Baidu Scholar

25

B. Guo, S. Bacha, M. Alamir et al., “A robust LESO-based DC-link voltage controller for variable speed hydro-electric plants,” in Proceedings of 2019 IEEE International Conference on Industrial Technology (ICIT), Melbourne, Australia, Jul. 2019, pp. 361-366. [Baidu Scholar

26

X. Zhou, Q. Liu, Y. Ma et al., “DC-Link voltage research of photovoltaic grid-connected inverter using improved active disturbance rejection control,” IEEE Access, vol. 9, pp. 9884-9894, Jan. 2021. [Baidu Scholar

27

X. Zhou, Q. Liu, Y. Ma et al., “A modified LADRC-based DC-link voltage controller for photovoltaic grid-connected inverters,” Electronics, vol. 10, no. 8, pp. 1-17, Apr. 2021. [Baidu Scholar

28

B. Guo, S. Bacha, M. Alamir et al., “Generalized integrator-extended state observer with applications to grid-connected converters in the presence of disturbances,” IEEE Transactions on Control Systems Technology, vol. 29, no. 2, pp. 744-755, Mar. 2021. [Baidu Scholar

29

G. Tian and Z. Gao, “Frequency response analysis of active disturbance rejection based control system,” in Proceedings of 2007 IEEE International Conference on Control Applications, Singapore, Nov. 2007, pp. 1595-1599. [Baidu Scholar

30

Q. Zheng and Z. Gao, “Active disturbance rejection control: between the formulation in time and the understanding in frequency,” Control Theory and Technology, vol. 14, no. 3, pp. 250-259, Aug. 2016. [Baidu Scholar

31

M. Merai, M. W. Naouar, I. Slama-Belkhodja et al., “An improved DC-link voltage control strategy for grid connected converters,” IEEE Transactions on Power Electronics, vol. 33, no. 4, pp. 3575-3582, Apr. 2018. [Baidu Scholar

32

Y. Liu, X. Zhou, H. Yu et al., “Sequence impedance modeling and stability assessment for load converters in weak grids,” IEEE Transactions on Industrial Electronics, vol. 68, no. 5, pp. 4056-4067, May 2021. [Baidu Scholar

33

Y. Liu, X. Zhou, Y. Chen et al., “Sequence impedance modeling and stability analysis for load converters with inertial support,” IEEE Transactions on Power Electronics, vol. 35, no. 12, pp. 13031-13041, Dec. 2020. [Baidu Scholar