Next Article in Journal
Analysis of the Effect of the Variable Charging Current Control Method on Cycle Life of Li-ion Batteries
Next Article in Special Issue
Control for Three-Phase LCL-Filter PWM Rectifier with BESS-Oriented Application
Previous Article in Journal
A Voltage-Based Approach for Series High Impedance Fault Detection and Location in Distribution Systems Using Smart Meters
Previous Article in Special Issue
Nonlinear Control of a Doubly Fed Generator Supplied by a Current Source Inverter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Control Strategy for Smooth Power Tracking of a Grid-Connected Virtual Synchronous Generator Based on Linear Active Disturbance Rejection Control

1
School of Electric Power Engineering, Nanjing Institute of Technology, Nanjing 211167, China
2
School of Energy and Environment, Southeast University, Nanjing 211189, China
3
School of Electrical and New Energy, China Three Gorges University, Yichang 443002, China
*
Author to whom correspondence should be addressed.
Energies 2019, 12(15), 3024; https://doi.org/10.3390/en12153024
Submission received: 15 June 2019 / Revised: 26 July 2019 / Accepted: 5 August 2019 / Published: 6 August 2019
(This article belongs to the Special Issue Control Strategies for Power Conversion Systems)

Abstract

:
The power quality of new energy resources has received tremendous attention recently. The control approach for the inverter, an interface between the new energy resources, and the infinite bus system is of vital importance. For the virtual synchronous generator (VSG), one of the research hotspots in the inverter control field, there are some challenges remaining to be dealt with. First is the contradiction between the rapid response and overshoot of active power output if VSG is connected to the grid. Secondly, the active power is deeply influenced by the fluctuation of gird frequency and this may bring power oscillation to VSG in weak grids. In this article, an active power controller for power tracking of grid-connected VSG is designed based on linear active disturbance rejection control (LADRC) by compensating for the lumped disturbance in a feedforward fashion. The parameters of the controller are analyzed and tuned in the frequency domain to acquire a desirable control performance. Moreover, the robustness of the control system is also considered. Simulation results show that the designed control system can transmit active power to the grid in a timely manner with no overshoot, as demanded. Additionally, it can output active power steadily according to the power reference without using a phase-locked loop (PLL) when the grid frequency has different features of fluctuation. In addition, the simulation results demonstrate that the improved VSG has strong robustness to the model parameter perturbation and mismatch.

1. Introduction

With the depletion of traditional energy sources and scalable applications of wind power and solar power, the penetration of electronic power equipment has increased rapidly in the grid. Traditionally, droop control strategy [1], active-reactive power (PQ) control strategy [2,3], and voltage-frequency (Vf) control strategy [3,4] can be chosen for inverter control, where droop control is capable of both the off-grid and on-grid modes. PQ control and Vf control are usually coordinated with each other to operate in both island and grid-connected modes for inverters. However, the characteristics of low inertia and damping of the three control strategies bring increasing uncertainty to the power grid operation.
Previous studies [5,6] proposed the concepts of a virtual synchronous generator (VSG) and virtual synchronous machine (VISMA) separately as alternatives, which emulate the active power-frequency regulating characteristics of a synchronous generator (SG) and synchronous machine (SM). The fluctuation of grid frequency is suppressed and the power quality is improved to a certain extent. The stability of the operation of the power system is also enhanced. VSG is now a research hotspot in microgrids.
VSG does not work at the rated operation point while working in island mode, and its frequency and voltage drop according to the droop characteristic curves of synchronous generators. If connected to the infinite bus system, VSG embodies the rotor inertia characteristic of SG, which causes the overshoot of active power, and thus the following transition process. The contradiction between the adjusting time and overshoot of active power has received extensive attention in academia, such that several different settling schemes have been put forward. One study [7] proposed a control strategy for the frequency and power of VSG based on Bang-Bang control. The angular frequency of VSG and its rate of change are considered to adjust the inertia momentarily, and the fluctuation is reduced. The second order differential value of output angular frequency and the differential value of virtual mechanical power were introduced in previous work [8] to adjust the inertia actively. The method of adjusting the inertia dynamically can also be found in other studies [9,10,11]. The proposed strategies in previous work [7,8,9,10,11] can indeed reduce the overshoot, however, the noise of the output power and frequency would be amplified when taking derivatives and that may affect the dynamic performance of VSG. Another study [12] proposed an improved VSG algorithm that is raised to increase stability and decrease the overshoot of power by adding an angle frequency differential term. In another study [13], a lead-lag item is introduced into the control structure and the pole-zero positions of closed-loop transfer function are changed by adding a lead item to the original inertia block. The damping ratio increases and the response accelerates in the presence of the added zero point in the left half of the S plane. If the grid frequency falls to a lower value, the output current of conventional VSG will reach a high value. This may cause the deviation between the power reference and the output power, and may even result in the action of overcurrent protection [13]. Additionally, the characteristics of conventional VSG will be deteriorated if the inverter has weak ties to the grid, and it is likely to result in output power oscillation [14]. Therefore, to meet more requirements in different scenarios, it is necessary to propose a supplementary or coordinated control strategy for VSG, which can handle output power steadily according to the power reference generated by the dispatching center and have the offset of power eliminated. To suppress the deviation of active power with consideration for the frequency fluctuation of the grid, one study [15] proposed a control strategy that adds a feedforward block for the deviation of frequency, where the conventional damping control in VSG is removed and a derivative term is added into the transfer function from the active power reference to the output active power, thus the steady-state error of power is eliminated when the grid frequency has a step change. Considering the electricity quantity of the direct current (DC) side, one study [16] proposed a coordinated control strategy that can switch between the conventional droop mode and the constant power control mode based on VSG. However, there still exists overshoot in the adjustment period. In another study [17], the grid frequency is tracked and that is set as the reference of the inverter instantaneously, and that can eliminate the deviation of power. However, there are problems that remain to be solved in phase-locked loop (PLL), such as the slow response speed and the complexity of designing parameters for PLL. Hence, different types of improved PLL were proposed in previous work [18,19,20]. The existing research studies focus on different aspects separately, such as reducing the active power overshoot or suppressing the deviation of active power when the grid frequency fluctuates. It is observed that a control strategy for VSG has not been put forward until now that can both track the power reference smoothly without overshoot and realize the purpose of the offset-free regulation of power when the disturbance exists. Hence, a smooth power tracking strategy for grid-connected VSG is proposed in this article accordingly.
Active disturbance rejection control (ADRC) was proposed by Prof. Han, which does not rely on the precise model of the controlled object [21,22]. The unmodeled parts and the outside disturbance in the sum are seen as the lumped disturbance, which can be observed by the extended state observer (ESO) and compensated by a feedforward controller. Compared with PID control, disturbance suppression is more advanced for ADRC. However, a large number of parameters need to be tuned in ADRC and it is difficult to demonstrate the stability of the control system by classical control theories. Prof. Gao proposed the linear active disturbance rejection control (LADRC). The controller is linearized and parameterized in previous studies [23,24] that promote the application of ADRC in engineering.
On the basis of the above introduction, a smooth power tracking approach for grid-connected VSG is proposed based on LADRC, which has three advantages, as follows:
  • The proposed control strategy based on LADRC can have VSG output active power to the grid as demanded without overshoot in time compared with the conventional VSG.
  • This control strategy enables VSG to transmit power to the grid stably in accordance with the power reference without PLL and it is insensitive to the gird frequency fluctuation.
  • The unmodeled parts can be seen as constituent parts of lumped disturbance and the improved VSG has strong robustness to the parameter perturbation and model-plant mismatch.
This article is organized as follows. Section 2 introduces the basic principles of VSG briefly. In Section 3, the drawbacks of transmitting active power to the grid are analyzed by deriving the power transfer function based on the conventional VSG control strategy. Furthermore, the controlling approach for active power of grid-connected VSG is stated based on LADRC, which illustrates the details of modeling. The method for parameter tuning of the designed controller is obtained through frequency domain analyses in Section 4, which contributes to parameter tuning and the performance of the controller. Simulation results are shown in Section 5 to validate the proposed control strategy for power control of grid-connected VSG and Section 6 gives the conclusion of this article with future work proposals.

2. Basic Principles of VSG

To emulate the primary frequency modulation characteristic of SG, we assume that the number of pole pairs is one, and the classical equations of the second-order VSG for inverters can be described as
P m = P ref + 1 K f ( ω n ω )
E = E 0 + 1 K q ( Q ref Q e )
{ J d ω d t = P m P e ω D ( ω ω n ) θ = ω d t
where P m , P ref , P e , Q ref , and Q e are the mechanical power, the active power reference, the output power, the reactive power reference, and the output reactive power, respectively; ω and ω n are denoted as the rotor angular frequency and the rated angular frequency, respectively; E and E 0 indicate the voltage reference and the rated peak value of the three-phase voltage, respectively; K f and K q are represented as the active power-frequency droop coefficient and the reactive power-voltage droop coefficient, respectively; J is the virtual inertia, D is the damping coefficient, and θ is the electrical angle. The topology and control structure of the grid-connected VSG are shown in Figure 1.
As shown in Figure 1, U dc is the constant voltage source containing wind power, solar power, and energy storage system. Here, L 1 , L 2 , and C are the filter inductance, the line inductance, and the filter capacitance, respectively; I CA , I CB , and I CC are the three-phase currents of the capacitance; E A , E B , and E C , and U A , U B , and U C are introduced to represent the output three-phase voltages of the inverter and the three-phase voltages of the capacitor separately; I Cd and I Cq , and U d and U q are the d-axis and q-axis current of the capacitance, and d-axis and q-axis voltage of the capacitance, respectively; Z L represents the local loads and the power is denoted as P Local . The breaker is controlled to achieve the switch between the on-grid and off-grid modes of VSG. Besides, the method for pre-synchronization of VSG refers to previous work [25], which proposed the approach of orienting from the q-axis component of grid voltage. The control block is shown in Figure 2, where U is the peak value of the three-phase voltage of VSG, U g is the three-phase voltage RMS of the grid, and φ is the initial phase angle of U g . Additionally, to realize better performance of VSG, a proportional-integral controller is introduced into the reactive power control loop when the VSG is connected to the grid. The control block is shown in Figure 2 [17,26,27] and the coefficients of proportionality and integration are 1 / K q and K iq , respectively.

3. Control Strategies for VSG

3.1. Performance Analysis of Conventional Grid-Connected VSG

The equivalent topology of the gird-connected VSG and the grid can be illustrated as Figure 3.
As shown in Figure 3, the VSG is equivalent to a SG and it is connected to the grid by three-phase transmission lines, where Z is the line impedance and Z = R + j ω L 2 = R + j X . Here, δ is the power angle of VSG and U 0 is modeled as the three-phase voltage RMS of VSG. In general, the line is inductive owing to the presence of transformers and filter inductance, thus the apparent power transmitted to the grid by VSG can be given by the following equation [28]:
S = 3 U I = 3 U U 0 δ U g 0 jX = P e + j Q e = 3 U 0 U g X sin δ + j   ( 3 U 0 U g X cos δ 3 U g 2 X )
where δ is generally seen as a small value, sin δ δ . Then, P e in Equation (4) can be linearized as [12]:
P e = 3 U 0 U g δ X
Now, we let ω = ω n when calculating the torque deviation, and we also have the switching of electronic devices, the fine adjustment process of voltage, and the current double closed loop omitted. Combining Equations (1)–(3) with (5), the equivalent active power control block diagram of the grid-connected VSG can be obtained as Figure 4.
Here, ω g is the real-time angular grid frequency and C = 3 U 0 U g X . Letting Δ ω ng = ω n ω g , therefore, the output active power can be modeled as:
P e ( s ) = C K f K f J ω n s 2 + ( K f D ω n + 1 ) s + C K f P ref ( s ) + C ( K f J ω n s + K f D ω n + 1 ) K f J ω n s 2 + ( K f D ω n + 1 ) s + C K f Δ ω ng ( s )
We suppose that the grid frequency does not fluctuate, i.e., ω n = ω g . Then, Equation (6) can be transformed into:
P e ( s ) = C K f K f J ω n s 2 + ( K f D ω n + 1 ) s + C K f P ref ( s )
The coefficients of the characteristic equation in Equation (7) are all greater than zero, and the system is stable. Further, the damping ratio of Equation (7) can be obtained as:
ξ 1 = D 2 ω n J C + 1 2 K f 1 J ω n C
It is clear that ξ 1 increases, the overshoot decreases, and the system responds more slowly with the increase of D . If J decreases, the overshoot will increase, and thus the response speed will accelerate. There exists a contradiction to be solved between the rapid response and the overshoot of output power.
The final value theorem [29] can be applied to Equation (7), and the steady-state power output to the grid can be calculated as P e = P ref . However, in reality, the grid frequency is always fluctuating, and it is difficult for VSG to operate stably at the rated frequency. Therefore,
P e ( s ) Δ ω ng ( s ) = C ( K f J ω n s + K f D ω n + 1 ) K f J ω n s 2 + ( K f D ω n + 1 ) s + C K f
Supposing Δ ω ng is the unit step signal, the steady-state power deviation can be deduced as
Δ P e = D ω n + 1 K f
It can be perceived from Equation (10) that there is a negative correlation between K f and Δ P e , and a positive correlation between D and Δ P e separately. However, under the condition that ω n ω g , the presence of the difference between P e and P ref cannot be eliminated if K f and D are changed to reduce the deviation. The overload coefficient of the inverter shall not exceed 1.5 in practice generally, and if Δ ω ng is too large, the actual output power of the VSG may exceed the maximum permissible power of the inverter. As one example, the extreme operation of the VSG may trigger its overcurrent protection, which will cause it to detach from the grid; as another example, if VSG is not disconnected from the grid, the long period of overloading can shorten its life. Hence, a strategy that enables VSG output power according to the reference is needed if ω g fluctuates or in some other extreme circumstances.

3.2. Control Strategy for Power of Grid-Connected VSG Based on LADRC

3.2.1. Principles of the Second-Order LADRC

It can be observed that the controlled object (6) is a second-order model, thus we design a second-order LADRC power controller to solve the above problems accordingly. The second-order LADRC controller consists of a linear extended state observer (LESO) and linear state error feedback (LSEF). The advantage of LADRC is that LESO can regard the lumped disturbance as one of the variable states of the system to be observed. Its basic control structure is shown in Figure 5 [30], where r is the reference and y is the output of the controlled object. The estimation for the output of the controlled object, the derivative value of the output, and the lumped disturbance are denoted as z 1 , z 2 , and z 3 , respectively. Additionally, b is the gain of the controlled object, u is the compensated control value, e is the difference between y and z 1 ; l 1 , l 2 , and l 3 are the different error feedback gain factors in LESO.
Here, f is denoted as the lumped disturbance, and the time-domain expression of the controlled object can be written as:
y ¨ = f ( y , y ˙ , w , t ) + b u
The model of LESO can be designed as:
{ e = y z 1 z ˙ 1 = z 2 + l 1 e z ˙ 2 = z 3 + l 2 e + b u z ˙ 3 = l 3 e
Noticing Equation (11), the expression of the controlled object can be written in state space form as:
[ x ˙ 1 x ˙ 2 x ˙ 3 ] = [ 0 1 0 0 0 1 0 0 0 ] [ x 1 x 2 x 3 ] + [ 0 b 0 ] u + [ 0 0 1 ] f
y = [ 1 0 0 ] [ y y ˙ f ]
Combining Equation (12) with Equations (13) and (14), the expression of the three-order LESO can be deduced as:
[ z ˙ 1 z ˙ 2 z ˙ 3 ] = [ l 1 1 0 l 2 0 1 l 3 0 0 ] [ z 1 z 2 z 3 ] + [ 0 b 0 l 1 l 2 l 3 ] [ u y ]
Then, taking the Laplace transformation of Equation (15), we have:
Z ( s ) = 1 s 3 + l 1 s 2 + l 2 s + l 3 [ b s l 1 s 2 + l 2 s + l 3 b s 2 + b l 1 s l 2 s 2 + l 3 s b l 3 l 3 s 2 ] [ U ( s ) Y ( s ) ]
Based on the estimated states, LSEF can be designed as:
u = K p ( r z 1 ) K d z 2 z 3 b
where K p and K d are the gains of LSEF. According to Equations (16) and (17), it can be seen that:
U ( s ) = 1 b [ K p R ( s ) C ( s ) Y ( s ) H ( s ) C ( s ) ]
where
C ( s ) = s 3 + l 1 s 2 + l 2 s + l 3 s 3 + ( l 1 + K d ) s 2 + ( l 1 K d + l 2 + K p ) s
H ( s ) = ( K p l 1 + k d l 2 + l 3 ) s 2 + ( K p l 2 + k d l 3 ) s + K P l 3 s 3 + l 1 s 2 + l 2 s + l 3
Combining Equation (11) with Equations (18)–(20), the simplified control structure can be drawn in Figure 6.
As shown in Figure 6, the control structure of active power does not become complicated, and it can be deduced as the classical feedback control structure where the disturbance is considered. Based on the estimations of output and disturbance, LSEF can compensate for the controlled variable in feedback, which will quickly make the output track the reference without offset. If the observations are reliable, then z 1 y , z 2 y ˙ , z 3 f . According to Equations (12) and (17), the differential equation can be derived as:
y ¨ = K d y ˙ K p y + K p r
The closed-loop transfer function of the system can be deduced by the classical control theory as:
Y ( s ) R ( s ) = K p s 2 + K d s + K p
A previous study [24] proposed a parameter tuning method for LADRC based on the bandwidths, which can track the reference without overshoot in a timely manner. Therefore, l 1 , l 2 , and l 3 , and K p and K d are relevant to the bandwidth of LESO and that of LSEF, which are noted as ω o and ω c , respectively. The tuning details are given by:
{ l 1 = 3 ω o ,   l 2 = 3 ω o 2 ,   l 3 = ω o 3 K d = 2 ω c ,   K p = ω c 2

3.2.2. Control Strategy for Active Power of VSG Based on LADRC

According to Equation (11), Equation (6) can be derived in the time domain as:
d 2 P e d t 2 = ( K f D ω n + 1 ) K f J ω n d P e d t C P e J ω n + C d Δ ω ng d t + C ( K f D ω n + 1 ) K f J ω n Δ ω ng + C J ω n P ref = f + C J ω n P ref = f + b u
Thus, the equivalent control structure of power based on LADRC can be obtained according to Section 3.1 as Figure 7.
Here, P u is the compensated control value and the transfer function of power can be stated as:
P e ( s ) = N ( s ) Q ( s ) P ref ( s ) + M ( s ) Q ( s ) [ ω n ( s ) ω g ( s ) ] = N ( s ) Q ( s ) P ref ( s ) + M ( s ) Q ( s ) Δ ω ng ( s )
where
N ( s ) = C K f ω c 2 ( s 3 + 3 ω o s 2 + 3 ω o 2 s + ω o 3 )
M ( s ) = b C K f J ω n s 4 + b C [ K f J ω n ( 3 ω o + 2 ω c ) + K f D ω n + 1 ] s 3 + b C [ K f J ω n ( 6 ω o ω c + 3 ω o 2 + ω c 2 ) + ( K f D ω n + 1 ) ( 3 ω o + 2 ω c ) ] s 2 + b C ( K f D ω n + 1 ) ( 6 ω o ω c + 3 ω o 2 + ω c 2 ) s
Q ( s ) = q 1 s 5 + q 2 s 4 + q 3 s 3 + q 4 s 2 + q 5 s + q 6
In Equation (28), q 1 q 6 are the constant values obtained by deducing the transfer function in Figure 7. It is obvious that the transient process of power is influenced by two factors directly, which are the power reference and Δ ω ng , respectively. Besides, according to Equations (27) and (28), the steady-state power output has no tracking offset if Δ ω ng is a step change value. It indicates the ability to suppress the disturbance under the LADRC control strategy for VSG. In addition, it can be further seen in Equation (24) that Δ ω ng is not the only part of the lumped disturbance. The influence of the difference between power reference and the output active power and its differential value, the influence of the difference of angular frequency reference, and the grid angular frequency and its differential value are all taken as the components of the lumped disturbance. Furthermore, some other factors that are unmodeled, which bring deviation to the active power, are all regarded as the components of lump disturbance and compensated. If the lumped disturbance can be observed and compensated precisely, an integrator series structure can be obtained. Therefore, the output can track the reference without overshoot and deviation, even though the plant is not accurate enough. This further shows the advantage that LADRC does not rely on the precise model of the controlled object.

4. Parameters Tuning for Controller

4.1. Nominal Performance Concerning

The stability of the control system is the premise for LADRC to achieve great performance in the power control of the grid-connected VSG. The stability margin of LADRC should be determined and the controller parameters need to be tuned first. According to previous work [24], there are three parameters remaining to be determined: b , ω o , and ω c . Now, we combine Equation (24) with Table 1, and b can be calculated as:
b = C J ω n = 4597
In order to clarify the relationship between the parameters when the system is stable, the characteristic equation of the closed loop transfer function can be obtained as:
q 1 s 5 + q 2 s 4 + q 3 s 3 + q 4 s 2 + q 5 s + q 6 = 0
where
q 1 = b K f J ω n
q 2 = b K f J ω n ( 3 ω o + 2 ω c ) + b ( K f D ω n + 1 )
q 3 = b K f J ω n ( 6 ω o ω c + ω c 2 + 3 ω o 2 ) + b ( K f D ω n + 1 ) ( 3 ω o + 2 ω c ) + b C K f
q 4 = b ( K f D ω n + 1 ) ( 6 ω o ω c + ω c 2 + 3 ω o 2 ) + C K f ( ω o 3 + 3 ω o ω c 2 + 6 ω c ω o 2 + 3 b ω o + 2 b ω c )
q 5 = C K f ( 6 b ω o ω c + 3 b ω o 2 + b ω c 2 + 3 ω o 2 ω c 2 + 2 ω c ω o 3 )
q 6 = C K f ω c 2 ω o 3
The Lienard–Chipard criterion can be applied to judge whether the system is stable or not. The necessary and sufficient condition for stability is that the coefficients of the characteristic equation and each odd (or even) order Hurwitz determinant are greater than zero separately. It is noticeable that the coefficients of the characteristic equation are all greater than zero. Therefore, the necessary and sufficient condition for the stability of the system can be simplified as [31]:
Δ 5 = ( q 2 q 3 q 1 q 4 ) ( q 4 q 5 q 3 q 6 ) ( q 2 q 5 q 1 q 6 ) 2 > 0
Thus, preliminarily, the relationship between ω o and ω c can be determined when the system is stable, according to the constraint in Equation (37). However, the Lienard–Chipard criterion cannot reflect the degree of stability. To obtain the reliable parameters of LADRC further, it is necessary to analyze the influence of different parameters on the performance of the LADRC controller. Combined with the above analyses, the initial values are chosen as ω o 0 = 300 and ω c 0 = 300 . The influence on the performance of the controller will be discussed with the changes of ω o and ω c .

4.1.1. Parameter ω c Increases When ω o Remains Unchanged

The parameter ω c increases from 0.25 ω c 0 to 4 ω c 0 gradually. The pole-zero map of the transfer function from P ref to P e and the Bode diagram of the transfer function from Δ ω ng to P e are given in Figure 8 and Figure 9, respectively.
As we can see in Figure 8, the closed-loop poles move away from the imaginary axis gradually and the response of the output accelerates with the increase of ω c . However, the angle between the dominant poles and the negative real axis increases first and then decreases, which can be perceived as the damping of the system experiencing the process of first decreasing and then increasing. When ω c is determined as the smallest value of 0.25 ω c 0 , the system has the largest damping value, which means that a larger value of ω c may cause overshoot of the power output. In addition, it can be perceived that ω c does not affect the position of the zero point. According to Appendix A, if the value of ω c is too large, the high-frequency measurement noise would be amplified more obviously, and this will deteriorate the control performance. In Figure 9, with the increase of the frequency from the band of 1 rad / s to 100 rad / s , the ability of anti-disturbance will be reduced. Besides, if ω c increases, the gain of Δ ω ng will decrease, and the output is less sensitive to the low frequency disturbance. This is to say, if the disturbance rejecting ability needs to be enhanced, we can increase ω c on the premise that the measurement noise is bearable for the system.

4.1.2. Parameter ω o Increases When ω c Remains Unchanged

The parameter ω o increases from 0.25 ω o 0 to 4 ω o 0 gradually. The pole-zero map of the transfer function from P ref to P e and the Bode diagram of the transfer function from Δ ω ng to P e are given in Figure 10 and Figure 11, respectively.
It can be seen from Figure 10 that the closed-loop poles move away from the imaginary axis gradually with the increase of ω o , which is similar to Figure 8. The increase of ω o can also accelerate the response speed of the system, which can track the power reference more rapidly. When ω o is small, the closed-loop zero points and dominant poles form dipoles and the equivalent order of the system can be reduced to some extent. However, with the increase of ω o , the dipoles no longer exist. Besides, the ability of rejecting disturbance increases when ω o increases, as shown in Figure 11.
At the same time, as introduced in Appendix A, the measurement noise would be amplified more evidently with the increase of ω o . Moreover, depending on the combinations of Figure 8 and Figure 10, a phenomenon can be seen showing that ω o has a deeper influence on the positions of dominant poles and the dynamic characteristics of the system when the two parameters have the same change. When the frequency of disturbance is 1 rad/s, the range of the gain of Δ ω ng , which is a part of the lumped disturbance, is from 23.2 dB to 51.8 dB and from 15.6 dB to 59.4 dB in Figure 9 and Figure 11, respectively. When the values of the two parameters have the same change, the change of ω o has a stronger influence on the ability to suppress the disturbance than the change of ω c .
Additionally, b s is the set value of b in the controller. In this article, the value of b can be calculated out. Supposing the condition that the gain of the model is hard to obtain, we can draw the pole-zero map of N ( s ) Q ( s ) and the Bode diagram of M ( s ) Q ( s ) as Figure 12 and Figure 13, respectively. Therefore, the influences of different values of b s on the control performance can be analyzed, letting ω o 0 = ω c 0 = 300 and b s change from 0.05 b to 3 b .
It can be seen in Figure 12 that when b s is smaller than 0.2 b , the system will collapse. Further, if b s is determined as a value greater than approximately 0.2 b , the system will be stable. Moreover, it can be seen in Figure 13 that when b s is a relatively large value, the gain of disturbance will increase. In other words, there is a tradeoff between the stability of the system and the ability to suppress the disturbance when the gain of the controlled object is hard to obtain. We need to tune b s from a large value to a small value for stability consideration. When the system is stable, we can decrease the value of b s concerning disturbance-rejection.
Above all, the tuning method for the parameters can be summarized as follows:
  • The appropriate initial values for ω o and ω c should not be too large, as stated in the constraints given before. If the dynamic response speed is slow or the capability needs to be enhanced, ω o and ω c can be increased until the overshoot exists or the influence of high-frequency noise is obvious.
  • Owing to ω o having a more pronounced influence on the dynamic performance of the system, we can reduce ω o first until the overshoot disappears or the influence of noise is controllable. Then, we can increase ω c until the overshoot exists or the influence of the noise is obvious. After this, the tuning sequences can be exchanged until the ideal performance requirement is met.
  • If the gain of the controlled object is not clear, the reasonable initial value for b s should not be too small for stability consideration. If the system is unstable, we can increase b s . Besides, if the capability of rejecting the disturbance needs to be enhanced, we can then decrease b s .

4.2. Robust Performance Concerning

As far as the robustness of the control system is concerned, we can deduce the sensitivity function S ( s ) and the complementary sensitivity function T ( s ) of the system as [32,33]:
S ( s ) = Q ( s ) N ( s ) Q ( s )
T ( s ) = N ( s ) Q ( s )
The necessary and sufficient condition for the uncertain system’s robust stability is:
ω > 0 ,   T ( j ω ) < 1 | Δ ( j ω ) |
For application, the maximum magnitudes of S ( s ) and T ( s ) are denoted as:
M S = max ω > 0 | S ( j ω ) |
M T = max ω > 0 | T ( j ω ) |
As obtained in previous work [34], the range of M S is from 1.2 to 2.0 and the corresponding range of M T is from 1 to 1.5. It is known that a smaller M S can make the output of the system more insensitive to the uncertainties of the controlled plant, which can strengthen the system’s stability. If a model is not precise enough, the reasonable values of M S and M T are recommended as smaller values in a previous study [35]. Based on the above discussions, to guarantee robust operation, M S is determined during the range from 1.1 to 1.2 and M T is determined during the range from 1 to 1.5. On the basis of above analysis, we let ω o = 420 and ω c = 70 , so that M S = 1 . 14 and M T = 1 . The robustness can be guaranteed and the robust performance of the system will be further testified in Section 5 under the condition that the parameter of the controlled object perturbates.

5. Case Study

For clarification of this control approach, here we can include the operating procedures of the improved VSG based on the proposed control strategy in this article (called “LADRC-VSG”). Figure 14 shows the flow chart of the operation of LADRC-VSG.
  • The VSG starts.
  • The VSG runs under island-mode and the breaker maintains “off” status. Meanwhile, the controller detects the grid-connected power reference. If the grid-connected power reference does not exist, continue with Step 2, otherwise the pre-synchronization starts and we proceed to the next step.
  • The system judges whether pre-synchronization is completed or not. If it is uncompleted, the procedure continues pre-synchronization, or else it moves to Step 4.
  • The breaker is closed and the inverter runs under the mode of LADRC-VSG. Meanwhile, the system detects whether the grid-connected power reference is still there or not. If it is, the procedure of Step 4 continues, or else it goes back to Step 2.
Simulations are carried out in the MATLAB/Simulink environment to verify the proposed control strategy. The grid is set as an infinite bus power system and the internal impedance of the grid is zero. The dynamic responses are studied separately based on the conventional VSG, the VSG based on setting the grid frequency measured by PLL as the frequency reference of VSG [36] (called “PLL-VSG”), and LADRC-VSG. The initial phase angle of the grid voltage is set as 15 . Here, P ref is 0 kW, as is the case in island mode, and Q ref is 0 kvar. The VSG loads 20 kW at first, P ref changes to 40 kW at 0.5 s, and it has a step change to 60 kW at 2 s. The lower limit of the control value of LADRC-VSG is set as 25 kW and the total simulation time is 3.5 s.
The simulation cases are divided into two parts, which include the nominal plant in Section 5.1 and the parametric perturbating plant in Section 5.2. Thus, the ability to suppress the disturbance and the robustness of the system will be verified separately. As observed, the cases studied in previous work [15,16] are based on the assumption that the form of grid frequency fluctuation is a step change. Further, the universality of the analysis needs to be considered, thus the cases based on the step change, ramp change, and sinewave form change of grid frequency are studied in Section 5.1, respectively. Moreover, the performance of the mismatching model is studied. The condition of a step change of grid frequency is introduced in detail, and other conditions are illustrated briefly.

5.1. Performance Based on the Nominal Model

5.1.1. Step Change of Grid Frequency

The grid frequency falls to 49.9 Hz at 2.5 s and is restored to normal at 3 s. The frequency of LADRC-VSG and the grid frequency measured by the PLL are depicted in Figure 15, where the PLL only works during the period of pre-synchronization to acquire the phase angle of the grid voltage before the VSG is connected to the grid. The method for pre-synchronization refers to a previous study [25], in which a PI controller is applied and the maximum frequency reaches the value of 50.7 Hz during the period. Additionally, the PLL withdraws after the switch of the breaker is closed for conventional the VSG and LADRC-VSG. From 0 s to 0.5 s, the inverter runs under the droop mode of VSG and the frequency of VSG is lower than the rated frequency. Pre-synchronization control starts at 0.5 s and the frequency of the VSG is identical to the grid frequency at about 1.16 s, thus VSG is connected to the grid and outputs power according to the power reference from then onwards.
The difference between the conventional VSG and PLL-VSG is that the PLL-VSG adopts the grid frequency as the frequency reference. Therefore, when the frequency does not fluctuate, the conventional VSG shares almost the same waveform with PLL-VSG (0–2.5 s). Figure 16 shows the output power of VSG under different control strategies. From 0 s to 0.5 s, the output power has the same value of 20 kW as the load power. Between 1.16 s and 2 s, the grid frequency is normal and the grid frequency is equal to the VSG reference frequency. Under the condition that the grid frequency is normal (between1.16 s and 2 s), we define
σ P e = P en ( max ) P ref P ref × 100 %
Under the circumstance that the grid frequency fluctuates (between 2.5 s and 3 s), we denote that
d P e = P ed ( max ) P ref
where σ P e is the percent overshoot of P e and P en ( max ) is the maximum output value of power when the grid frequency is normal (from 1.16 s to 2 s); d P e is the maximum deviation of P e , and P ed ( max ) is the maximum output value of power when the grid frequency fluctuates (from 2.5 s to 3 s). In Figure 16, we can see that the power overshoot of the conventional VSG and PLL-VSG still exists when connected to the grid with the damping ratio of 0.926. The P en ( max ) of the conventional VSG and PLL-VSG is 43.9 kW and the σ P e is 9.7%, whereas obviously LADRC-VSG can track the power reference rapidly with no overshoot, and the σ P e of LADRC-VSG is equal to zero. The adjustment time of conventional VSG and PLL-VSG is 0.19 s, which certifies the excellent performance of tracking the power reference of the designed controller. The grid frequency drops by 0.1 Hz from 2.5 s to 3 s, according to Equation (10), and the steady-state power output of the conventional VSG should be 79.74 kW at this time theoretically, which is in agreement with the simulation results. By contrast, the maximum output power of LADRC-VSG is 65.74 kW and there is no oscillation in the dynamic regulation process. As for PLL-VSG, the maximum output power is 69.5 kW. Thus, the d P e values of LADRC-VSG, conventional VSG, and PLL-VSG are 5.74 kW, 25.13 kW, and 9.5 kW from 2.5 s to 3 s, respectively. In the process of grid frequency recovery, the regulating time of LADRC-VSG is less than 0.1 s, which is lower than that of 0.22 s for conventional VSG and PLL-VSG. Further, it can be seen in Figure 17, under the condition that the active power is controlled precisely, that the oscillation intensity of the reactive power of LADRC-VSG is lower than conventional VSG and PLL-VSG. All of the above results show that LADRC is reliable in power control of the grid-connected VSG.
Figure 18, Figure 19, and Figure 20 show the output current waveforms of LADRC-VSG, conventional VSG, and PLL-VSG, respectively. The output current of VSG under LADRC control strategy has no overshoot and the current is relatively smooth. From 2.5 s to 3 s, the steady-state current of LADRC-VSG is 126.4 A and the current is almost not affected by the fluctuation of grid frequency. At the same time, the steady-state error of current of PLL-VSG can also be eliminated, whereas the maximum output current of the conventional VSG is up to 176 A and the steady-state error of the current is nearly 49.6 A. The goals of disturbance suppression and fast tracking for the power reference without overshoot are achieved.

5.1.2. Ramp Change of Grid Frequency

The grid frequency has a ramp change from 50 Hz at 2 s to 49.9 Hz at 2.5 s, and the grid frequency remains at 49.9 Hz after this point. The output active power waveforms of LADRC-VSG, conventional VSG, and PLL-VSG are shown in Figure 21. From 2.5 s to 3 s, the d P e values of LADRC-VSG, conventional VSG, and PLL-VSG are 0.6 kW, 19.8 kW, and 0.6 kW, respectively. The dynamic error of the power of PLL-VSG and LADRC-VSG is suppressed to a low value and the active power shows no sudden change when the grid frequency fluctuates. At this time, compared with PLL-VSG, the advantage of LADRC-VSG is that it can track the power reference without rapid overshoot.

5.1.3. Sinewave form Change of Grid Frequency

The grid frequency has a sinewave form change with the amplitude of 0.05 Hz at 2 s, which returns to normal at 2.5 s. The cycle of the sinewave form fluctuation is 0.5 s. The performances of LADRC-VSG, conventional VSG, and PLL-VSG can be seen in Figure 22.
The varying angular frequency of disturbance between 2.5 s and 3 s is 2 × 2 × π = 12.566   rad / s and the output active power of LADRC-VSG has a slight fluctuation. This phenomenon can be explained as the capability to suppress disturbance being reduced because the varying frequency of disturbance increases, and this is consistent with the analysis in Section 4, even though the d P e of LADRC-VSG is 1.2 kW, which is less than the 10.2 kW value of the conventional VSG and 1.9 kW value of PLL-VSG under the same condition. With different forms of disturbance of frequency, the deviation of power can all be suppressed to an ideal low value, even as low as zero, which reveals that LADRC-VSG is not sensitive to the change of the form of disturbance. The simulations in Section 5.1.1, Section 5.1.2 and Section 5.1.3 convince us that the designed LADRC-VSG controller can reject the external disturbance reliably and it can achieve steady VSG active output power according to the command.

5.2. Control Performance Based on the Parametric Perturbating Model

Section 5.1 reveals the effectiveness of the control method under the condition that the gird frequency has different forms of fluctuation. In real situations, however, the parameters of the filters, the transmission lines, and the inverters may be perturbed for various reasons, such as equipment aging, failure, and other factors. Here, we take the perturbation of the parameter of L 2 of the line inductance as an example to illustrate the robust stability of the system.
Here, L 2 changes from 0.404 mH to 0.1 mH and the damping ratio ξ 1 decreases from 0.926 to 0.461. For robustness testing, b s is not changed in the controller, other simulating conditions are as poor as the conditions in Section 5.1.3, and the grid frequency has a sinewave form change. The amplitude of change is 0.05 Hz and the cycle of the fluctuation is 0.5 s. The grid frequency recovers at 3 s and the simulation results are presented in Figure 23, Figure 24, Figure 25, Figure 26, Figure 27 and Figure 28.
As we can see from Figure 23, with the perturbation of L 2 , the P en ( max ) values of the conventional VSG and PLL-VSG increases to 54.87 kW between 1.16 s and 2 s, whereas the LADRC-VSG can continue transmitting active power to the grid without overshoot. The σ P e values of LADRC-VSG, conventional VSG, and PLL-VSG are 0, 37.17%, and 37.17%, respectively. Additionally, the d P e values of LADRC-VSG, conventional VSG, and PLL-VSG are 1.44 kW, 11.6 kW, and 4.9 kW, respectively. The adjustment time of LADRC-VSG is 0.12 s, which is much shorter than that the 0.8 s value for conventional VSG and PLL-VSG. LADRC-VSG can output power stably in the case of grid frequency fluctuations, which illustrates the advantage of LADRC, namely that it does not rely on a precise model of the controlled plant. The output reactive power waveforms of LADRC-VSG, conventional VSG, and PLL-VSG are shown in Figure 24. It can be seen that the oscillation of active power will cause instability of the reactive power of conventional VSG and PLL-VSG owing to the presence of coupling between active power and reactive power. Moreover, even though the three control strategies all adopt the method of PI control for reactive power, the better performance of active power of LADRC-VSG causes much less oscillation of reactive power. In addition, the corresponding output angular frequencies of the inverter under different control strategies are shown in Figure 25, where we can see that the angular frequency of LADRC-VSG has better performance. Compared with the performance of the conventional VSG and PLL-VSG, the model mismatching parts in LADRC-VSG can be seen as disturbance, and that can be observed by LESO and compensated by LSEF without delay; thus, the strong robustness can be guaranteed. The estimates for P e and d P e / d t and the lumped disturbance are shown in Figure 26, Figure 27 and Figure 28, respectively, where we can see that the observation for P e is of high-precision and the observed disturbance is unitless, as the disturbance is caused by a gathering of many factors in Figure 28. The simulations carried out can illustrate the stability and robustness of the system further, which testify to the validity of the proposed control strategy. The comparative analyses of the controllers’ index are given in Table 2 and Table 3.
As seen in the analyses above, LADRC-VSG can track the power reference without overshoot by compensating for the lumped disturbance, even though the model is mismatched. At the same time, it is noteworthy that the steady-state accuracy of the active power is reduced a little when the controlled object is mismatched, which can be perceived when comparing Figure 22 and Figure 23. As shown in the amplified picture in Figure 22, the steady-state output power is between 59.8 kW and 60.2 kW, whereas under the same condition (the grid frequency does not fluctuate between 1.16 s and 2.5 s), it can be observed that the steady-state output power is between 59.5 kW and 60.5 kW in the amplified picture in Figure 23. This phenomenon reveals that we require a relatively more exact value for the gain of the controlled object to achieve better performance of the controller when the precise model of the controlled object has many uncertainties, so that the capability to reject disturbance can be further enhanced.

6. Conclusions

A control strategy for power tracking of a grid-connected VSG based on LADRC is proposed in this article. On the basis of deducing the power equation of the conventional grid-connected VSG, this article clarifies that VSG has some shortcomings in the power regulation period, such as the contradiction between the response speed of the power output and output overshoot. Additionally, the overcurrent of VSG may cause its dispatch from the grid. Hence, the active power control strategy based on LADRC is proposed. The lumped disturbance is compensated, including the simplified part of the model, the influence of the difference of power reference and the output active power and its differential value, and the difference of the angular frequency reference and the grid angular frequency and its differential value. Then, by deriving the transfer function of the second-order LADRC, the equivalent structure of the grid-connected VSG based on LADRC strategy is obtained. After that, the influences of different parameters of LADRC on the dynamic characteristics of the system are analyzed through frequency domain analysis methods and the robust stability of the control system is analyzed further as well. The simulation results show that LADRC has great performance in power control of the grid-connected VSG, which can improve the power quality of the grid-connected inverters. When the grid frequency is normal, LADRC-VSG can output power according to the power reference quickly and without overshoot. Moreover, LADRC-VSG can maintain robust stability when the parameters perturbate and the analyses show that the LADRC-VSG has a greater performance index. The future vision will be focused on the application of LADRC in the pre-synchronization period of VSG and load-sharing control for VSGs in the microgrid.

Author Contributions

Conceptualization, Y.Z. and J.Z.; formal analysis Y.Z., J.Z., and X.D.; data curation, J.Z., and X.D.; writing—original draft preparation, Y.Z., P.G., and P.Z.; writing—review and editing, J.Z., X.D., and X.Z.

Funding

This research was funded by National Natural Science Foundation of China (No. 51607083).

Acknowledgments

All the authors are grateful for the support from National Natural Science Foundation of China and the first author thanks his family and friends for their encouragements over the years.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

If the measurement noise is considered, the control structure can be drawn as Figure A1.
Figure A1. Control structure for the power of the VSG based on LADRC when the measurement noise is considered.
Figure A1. Control structure for the power of the VSG based on LADRC when the measurement noise is considered.
Energies 12 03024 g0a1
When ω c changes from 0.25 ω c 0 to 4 ω c 0 , and ω o maintains the value of ω o 0 as unchanged, the Bode diagram of the transfer function from the measurement noise ( P noise ) to the output power ( P e ) is drawn as in Figure A2. When ω o increases from 0.25 ω o 0 to 4 ω o 0 and ω c maintains the value of ω c 0 as unchanged, the Bode diagram of the transfer function from P noise to the P e is drawn as in Figure A3. It can be seen in the two figures that the amplification of the high-frequency noise in the output becomes more obvious with the increase of ω c and ω o . Thus, when the system is influenced by the measurement noise or the performance of the system is challenged, we must decrease ω c , ω o , or both.
Figure A2. Bode diagram of P e ( s ) P noise ( s ) when ω c increases.
Figure A2. Bode diagram of P e ( s ) P noise ( s ) when ω c increases.
Energies 12 03024 g0a2
Figure A3. Bode diagram of P e ( s ) P noise ( s ) when ω o increases.
Figure A3. Bode diagram of P e ( s ) P noise ( s ) when ω o increases.
Energies 12 03024 g0a3

References

  1. Engler, A.; Soultanis, N. Droop control in LV-grids. In Proceedings of the 2005 International Conference on Future Power Systems, Amsterdam, The Netherlands, 16–18 November 2005; p. 6. [Google Scholar] [CrossRef]
  2. Chen, M.R.; Wang, H.; Zeng, G.Q.; Dai, Y.X.; Bi, D.Q. Optimal PQ Control of Grid-Connected Inverters in a Microgrid Based on Adaptive Population Extremal Optimization. Energies 2018, 11, 2107. [Google Scholar] [CrossRef]
  3. Adhikari, S.; Li, F. Coordinated Vf and PQ control of solar photovoltaic generators with MPPT and battery storage in microgrids. IEEE Trans. Smart Grid 2014, 5, 1270–1281. [Google Scholar] [CrossRef]
  4. Lopes, J.P.; Moreira, C.L.; Madureira, A.G.; Resende, F.O.; Wu, X.A.; Jayawarna, N.A.; Zhang, Y.A. Control strategies for microgrids emergency operation. In Proceedings of the 2005 International Conference on Future Power Systems, Amsterdam, The Netherlands, 16–18 November 2005; pp. 57–78. [Google Scholar] [CrossRef]
  5. Beck, H.P.; Hesse, R. Virtual synchronous machine. In Proceedings of the 2007 9th International Conference on Electrical Power Quality and Utilisation, Barcelona, Spain, 9–11 October 2007; pp. 1–6. [Google Scholar] [CrossRef]
  6. Visscher, K.; De Haan, S.W.H. Virtual synchronous machines (VSG’s) for frequency stabilisation in future grids with a significant share of decentralized generation. In Proceedings of the CIRED Seminar 2008: SmartGrids for Distribution, Frankfurt, Germany, 23–24 June 2008; pp. 1–4. [Google Scholar] [CrossRef]
  7. Alipoor, J.; Miura, Y.; Ise, T. Power system stabilization using virtual synchronous generator with alternating moment of inertia. IEEE J. Emerg. Sel. Top. Power Electron. 2015, 3, 451–458. [Google Scholar] [CrossRef]
  8. Vaddipalli, R.K.N.; Lopes, L.A.; Rathore, A.K. Virtual Synchronous Generator with Variable Inertia Emulation Via Power Tracking Algorithm. In Proceedings of the 2018 2nd International Conference on Power, Energy and Environment: Towards Smart Technology (ICEPE), Shillong, India, 1–2 June 2018; pp. 1–6. [Google Scholar] [CrossRef]
  9. Fan, W.; Yan, X.; Hua, T. Adaptive parameter control strategy of VSG for improving system transient stability. In Proceedings of the 2017 IEEE 3rd International Future Energy Electronics Conference and ECCE Asia (IFEEC 2017-ECCE Asia), Kaohsiung, Taiwan, 3–7 June 2017; pp. 2053–2058. [Google Scholar] [CrossRef]
  10. Li, D.; Zhu, Q.; Lin, S.; Bian, X.Y. A self-adaptive inertia and damping combination control of VSG to support frequency stability. IEEE Trans. Energy Convers. 2017, 32, 397–398. [Google Scholar] [CrossRef]
  11. Rahman, F.S.; Kerdphol, T.; Watanabe, M.; Mitani, Y. Optimization of virtual inertia considering system frequency protection scheme. Electr. Power Syst. Res. 2019, 170, 294–302. [Google Scholar] [CrossRef]
  12. Haizhen, X.; Xing, Z.; Fang, L.; Fubin, M.; Rongliang, S.; Hua, N. An improved virtual synchronous generator algorithm for system stability enhancement. In Proceedings of the 2015 IEEE 2nd International Future Energy Electronics Conference (IFEEC), Taipei, Taiwan, 1–4 November 2015; pp. 1–6. [Google Scholar] [CrossRef]
  13. Xu, H.; Zhang, X.; Liu, F.; Mao, F.; Shi, R.; Yu, C.; Yu, Y. Virtual Synchronous Generator Control Strategy Based on Lead-lag Link Virtual Inertia. Proc. CSEE 2017, 37, 1918–1926. [Google Scholar] [CrossRef]
  14. Li, X.; Hu, Y.; Shao, Y.; Chen, G. Mechanism analysis and suppression strategies of power oscillation for virtual synchronous generator. In Proceedings of the IECON 2017—43rd Annual Conference of the IEEE Industrial Electronics Society, Beijing, China, 5–8 November 2017; pp. 4955–4960. [Google Scholar] [CrossRef]
  15. Li, M.; Wang, Y.; Zhou, H.; Hu, W. A phase feedforward based virtual synchronous generator control scheme. In Proceedings of the 2018 IEEE Applied Power Electronics Conference and Exposition (APEC), San Antonio, TX, USA, 4–8 March 2018; pp. 3314–3318. [Google Scholar] [CrossRef]
  16. Yan, X.; Jia, J.; Li, Z.; Ge, Y. Power Control and Smooth Mode Switchover for Grid-connected Virtual Synchronous Generators. Autom. Electr. Power Syst. 2018, 42, 91–99. [Google Scholar] [CrossRef]
  17. Shintai, T.; Miura, Y.; Ise, T. Oscillation damping of a distributed generator using a virtual synchronous generator. IEEE Trans. Power Deliv. 2014, 29, 668–676. [Google Scholar] [CrossRef]
  18. Khan, S.; Bletterie, B.; Anta, A.; Gawlik, W. On Small Signal Frequency Stability under Virtual Inertia and the Role of PLLs. Energies 2018, 11, 2372. [Google Scholar] [CrossRef]
  19. He, X.; Geng, H.; Yang, G. A generalized design framework of notch filter based frequency-locked loop for three-phase grid voltage. IEEE Trans. Ind. Electron. 2018, 65, 7072–7084. [Google Scholar] [CrossRef]
  20. Du, H.; Sun, Q.; Cheng, Q.; Ma, D.; Wang, X. An Adaptive Frequency Phase-Locked Loop Based on a Third Order Generalized Integral. Energies 2019, 12, 309. [Google Scholar] [CrossRef]
  21. Han, J.-Q. Auto-disturbances-rejection controller and its’ applications. Control Decis. 1998, 1, 19–23. [Google Scholar] [CrossRef]
  22. Han, J. From PID to active disturbance rejection control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  23. Gao, Z. Active disturbance rejection control: A paradigm shift in feedback control system design. In Proceedings of the 2006 American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; p. 7. [Google Scholar] [CrossRef]
  24. Gao, Z. Scaling and bandwidth-parameterization based controller tuning. In Proceedings of the 2003 American Control Conference, Denver, CO, USA, 4–6 June 2003; pp. 4989–4996. [Google Scholar] [CrossRef]
  25. Lv, Z.; Sheng, W.; Zhong, Q.; Liu, H.; Zeng, Z.; Yang, L.; Liu, L. Virtual Synchronous Generator and Its Applications in Micro-grid. Proc. CSEE 2014, 34, 2591–2603. [Google Scholar] [CrossRef]
  26. Shintai, T.; Miura, Y.; Ise, T. Reactive power control for load sharing with virtual synchronous generator control. In Proceedings of the 7th International Power Electronics and Motion Control Conference, Harbin, China, 2–5 June 2012; pp. 846–853. [Google Scholar] [CrossRef]
  27. Zheng, T.; Chen, L.; Guo, Y.; Mei, S. Comprehensive control strategy of virtual synchronous generator under unbalanced voltage conditions. IET Gener. Transm. Distrib. 2017, 12, 1621–1630. [Google Scholar] [CrossRef]
  28. Hu, W.; Wu, Z.; Sun, C.; Song, Y.; Yuan, K. Modeling and parameter setting method for grid-connected inverter of energy storage system based on VSG. Electr. Power Autom. Equip. 2018, 13–23. [Google Scholar] [CrossRef]
  29. Chen, J.; Lundberg, K.H.; Davison, D.E.; Bernstein, D.S. The final value theorem revisited-infinite limits and irrational functions. IEEE Control Syst. Mag. 2007, 27, 97–99. [Google Scholar] [CrossRef]
  30. Zhang, F.; Chen, X.; Zhang, Y.; Jiang, H. Research on Power Transmission Capability of VSC-HVDC Based on Second Order LADRC. Energy Procedia 2019, 158, 2592–2598. [Google Scholar] [CrossRef]
  31. Yuan, D.; Ma, X.; Zeng, Q.; Qiu, X. Research on frequency-band characteristics and parameters configuration of linear active disturbance rejection control for second-order systems. Control Theory Appl. 2013, 30, 1630–1640. [Google Scholar] [CrossRef]
  32. Skogestad, S.; Postlethwaite, I. Multivariable Feedback Control: Analysis and Design, 2nd ed.; Wiley: New York, NY, USA, 2007; p. 22. [Google Scholar]
  33. Wu, M.; He, Y.; Yu, J. Robust Control Theory; Higher Education Press: Beijing, China, 2010; pp. 24–32. [Google Scholar]
  34. Åström, K.J.; Panagopoulos, H.; Hägglund, T. Design of PI controllers based on non-convex optimization. Automatica 1998, 34, 585–601. [Google Scholar] [CrossRef]
  35. Sun, L.; Li, D.; Gao, Z.; Yang, Z.; Zhao, S. Combined feedforward and model-assisted active disturbance rejection control for non-minimum phase system. ISA Trans. 2016, 64, 24–33. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, F. Research on Microgrid Inverter Control Strategy Based on Virtual Synchronous Generator. Ph.D. Thesis, Hefei University of Technology, Hefei, China, April 2015. [Google Scholar]
Figure 1. Topology and general control structure of grid-connected inverter.
Figure 1. Topology and general control structure of grid-connected inverter.
Energies 12 03024 g001
Figure 2. Overall schematic control diagram of grid-connected virtual synchronous generator (VSG).
Figure 2. Overall schematic control diagram of grid-connected virtual synchronous generator (VSG).
Energies 12 03024 g002
Figure 3. Equivalent diagram of grid-connected VSG.
Figure 3. Equivalent diagram of grid-connected VSG.
Energies 12 03024 g003
Figure 4. Power control diagram of grid-connected VSG.
Figure 4. Power control diagram of grid-connected VSG.
Energies 12 03024 g004
Figure 5. Structure of second-order linear active disturbance rejection control (LADRC).
Figure 5. Structure of second-order linear active disturbance rejection control (LADRC).
Energies 12 03024 g005
Figure 6. Simplified structure of LADRC.
Figure 6. Simplified structure of LADRC.
Energies 12 03024 g006
Figure 7. Control structure for power of VSG based on LADRC.
Figure 7. Control structure for power of VSG based on LADRC.
Energies 12 03024 g007
Figure 8. Pole-zero map of N ( s ) Q ( s ) with changing ω c (the “ × ” represents the pole points and “ ” represents zero points; the same below).
Figure 8. Pole-zero map of N ( s ) Q ( s ) with changing ω c (the “ × ” represents the pole points and “ ” represents zero points; the same below).
Energies 12 03024 g008
Figure 9. Bode diagram of M ( s ) Q ( s ) with changing ω c (different colors represent different gains; the same below).
Figure 9. Bode diagram of M ( s ) Q ( s ) with changing ω c (different colors represent different gains; the same below).
Energies 12 03024 g009
Figure 10. Pole-zero map of N ( s ) Q ( s ) with changing ω o .
Figure 10. Pole-zero map of N ( s ) Q ( s ) with changing ω o .
Energies 12 03024 g010
Figure 11. Bode diagram of M ( s ) Q ( s ) with changing ω o .
Figure 11. Bode diagram of M ( s ) Q ( s ) with changing ω o .
Energies 12 03024 g011
Figure 12. Pole-zero map of N ( s ) Q ( s ) when b s changes.
Figure 12. Pole-zero map of N ( s ) Q ( s ) when b s changes.
Energies 12 03024 g012
Figure 13. Bode diagram of M ( s ) Q ( s ) when b s changes.
Figure 13. Bode diagram of M ( s ) Q ( s ) when b s changes.
Energies 12 03024 g013
Figure 14. The procedure of the operation of linear active disturbance rejection control-virtual synchronous generator (LADRC-VSG).
Figure 14. The procedure of the operation of linear active disturbance rejection control-virtual synchronous generator (LADRC-VSG).
Energies 12 03024 g014
Figure 15. Frequencies of LADRC-VSG and grid measured by phase-locked loop (PLL).
Figure 15. Frequencies of LADRC-VSG and grid measured by phase-locked loop (PLL).
Energies 12 03024 g015
Figure 16. Output active power of VSG under different control strategies.
Figure 16. Output active power of VSG under different control strategies.
Energies 12 03024 g016
Figure 17. Corresponding output reactive power of VSG under different control strategies.
Figure 17. Corresponding output reactive power of VSG under different control strategies.
Energies 12 03024 g017
Figure 18. Output current (Phase A) of LADRC-VSG.
Figure 18. Output current (Phase A) of LADRC-VSG.
Energies 12 03024 g018
Figure 19. Output current (Phase A) of conventional VSG.
Figure 19. Output current (Phase A) of conventional VSG.
Energies 12 03024 g019
Figure 20. Output current (Phase A) of PLL-VSG.
Figure 20. Output current (Phase A) of PLL-VSG.
Energies 12 03024 g020
Figure 21. Output power of VSG under different control strategies.
Figure 21. Output power of VSG under different control strategies.
Energies 12 03024 g021
Figure 22. Output power of VSG under different control strategies.
Figure 22. Output power of VSG under different control strategies.
Energies 12 03024 g022
Figure 23. Output power of VSG under different control strategies.
Figure 23. Output power of VSG under different control strategies.
Energies 12 03024 g023
Figure 24. Corresponding output reactive power of VSG under different control strategies.
Figure 24. Corresponding output reactive power of VSG under different control strategies.
Energies 12 03024 g024
Figure 25. Output angular frequencies of inverter under different control strategies.
Figure 25. Output angular frequencies of inverter under different control strategies.
Energies 12 03024 g025
Figure 26. Output power observed by linear extended state observer (LESO).
Figure 26. Output power observed by linear extended state observer (LESO).
Energies 12 03024 g026
Figure 27. Derivative value of active power observed by LESO.
Figure 27. Derivative value of active power observed by LESO.
Energies 12 03024 g027
Figure 28. The lumped disturbance observed by LESO.
Figure 28. The lumped disturbance observed by LESO.
Energies 12 03024 g028
Table 1. Key parameters of the system.
Table 1. Key parameters of the system.
ParameterSymbolValueParameterSymbolValue
U g V 220 U dc V 800
D N m rad 1 100 J kg m 2 0.8
K q 3330 L 1 mH 0.6
K iq 0.005 L 2 mH 0.404
K f 0.0628 C 1 μ F 1500
R Ω 0.1 ω n rad 1 s 314.16
Table 2. Comparison of d P e (2.5–3 s, one decimal place reserved).
Table 2. Comparison of d P e (2.5–3 s, one decimal place reserved).
Type of DisturbanceLADRC-VSG (kW)VSG (kW)PLL-VSG (kW)
Step change of frequency5.725.19.5
Ramp change of frequency0.619.80.6
Sinewave form change of frequency1.210.21.9
Model mismatch; Sinewave form change of frequency1.411.64.9
Table 3. Comparison of σ P e (1.16–2 s, one decimal place reserved).
Table 3. Comparison of σ P e (1.16–2 s, one decimal place reserved).
Type of ModelLADRC-VSG (%)VSG (%)PLL-VSG (%)
Normal model09.79.7
Mismatch model037.237.2

Share and Cite

MDPI and ACS Style

Zhang, Y.; Zhu, J.; Dong, X.; Zhao, P.; Ge, P.; Zhang, X. A Control Strategy for Smooth Power Tracking of a Grid-Connected Virtual Synchronous Generator Based on Linear Active Disturbance Rejection Control. Energies 2019, 12, 3024. https://doi.org/10.3390/en12153024

AMA Style

Zhang Y, Zhu J, Dong X, Zhao P, Ge P, Zhang X. A Control Strategy for Smooth Power Tracking of a Grid-Connected Virtual Synchronous Generator Based on Linear Active Disturbance Rejection Control. Energies. 2019; 12(15):3024. https://doi.org/10.3390/en12153024

Chicago/Turabian Style

Zhang, Yaya, Jianzhong Zhu, Xueyu Dong, Pinchao Zhao, Peng Ge, and Xiaolian Zhang. 2019. "A Control Strategy for Smooth Power Tracking of a Grid-Connected Virtual Synchronous Generator Based on Linear Active Disturbance Rejection Control" Energies 12, no. 15: 3024. https://doi.org/10.3390/en12153024

APA Style

Zhang, Y., Zhu, J., Dong, X., Zhao, P., Ge, P., & Zhang, X. (2019). A Control Strategy for Smooth Power Tracking of a Grid-Connected Virtual Synchronous Generator Based on Linear Active Disturbance Rejection Control. Energies, 12(15), 3024. https://doi.org/10.3390/en12153024

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop