Next Article in Journal
Experimental Study on Sensitivity of Porosity to Pressure and Particle Size in Loose Coal Media
Next Article in Special Issue
Optimization of Solar Energy System for the Electric Vehicle at University Campus in Dhaka, Bangladesh
Previous Article in Journal
Resilience Assessment and Its Enhancement in Tackling Adverse Impact of Ice Disasters for Power Transmission Systems
Previous Article in Special Issue
Resilient Multiscale Coordination Control against Adversarial Nodes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lyapunov-Based Large Signal Stability Assessment for VSG Controlled Inverter-Interfaced Distributed Generators

Shanghai Jiao Tong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Energies 2018, 11(9), 2273; https://doi.org/10.3390/en11092273
Submission received: 9 August 2018 / Revised: 22 August 2018 / Accepted: 27 August 2018 / Published: 29 August 2018
(This article belongs to the Special Issue Distribution System Operation and Control)

Abstract

:
Inverter-interfaced distributed generators (IIDGs) have been widely applied due to their control flexibility. The stability problems of IIDGs under large signal disturbances, such as large load variations and feeder faults, will cause serious impacts on the system. The virtual synchronous generator (VSG) control is an effective scheme for IIDGs to increase transient stability. However, the existing linearized stability models of IIDGs are limited to small disturbances. Hence, this paper proposes a Lyapunov approach based on non-linearized models to assess the large signal stability of VSG-IIDG. The electrostatic machine model is introduced to establish the equivalent nonlinear model. On the basis of Popov’s theory, a Lyapunov function is derived to calculate the transient stability domain. The stability mechanism is revealed by depicting the stability domain using the locus of the angle and the angular frequency. Large signal stability of the VSG-IIDG is quantified according to the boundary of the stability domain. Effects and sensitivity analysis of the key parameters including the cable impedance, the load power, and the virtual inertia on the stability of the VSG-IIDG are also presented. The simulations are performed in PSCAD/EMTDC and the results demonstrate the proposed large signal stability assessment method.

1. Introduction

Inverter-interfaced distributed generators (IIDGs) feature in flexible control and quick response uses in different application scenarios. However, unlike conventional synchronous generators, inverters lack inertia [1]. High penetration of IIDGs may result in poor voltage and frequency response, and even instability with large-scale oscillation, asynchronism, and voltage collapse under large disturbances [2]. The virtual synchronous generator (VSG) control scheme is an effective solution to the assigned problem. By controlling the switching pattern of the inverter, the VSG emulates conventional synchronous generators. The VSG means to provide virtual inertia and additional damping that can reduce frequency deviations during disturbances [3,4]. The VSG controlled IIDG (VSG-IIDG) is becoming one of the most prospective renewable sources due to the outstanding features [5].
The stability issues of VSG-IIDGs are categorized into two types: small signal stability and large signal stability. Small signal stability of IIDG has been relatively perfected using the concept of synchronous generators in references [6,7]. The study of small signal stability uses the linearized model. Nyquist or Routh-Hurwitz stability criterion, eigenvalue analysis [8,9,10,11], and transfer function [12,13] are the regular tools to assess the small signal stability. However, the analysis of the small signal linearized model is only valid around the stable operating point yet not accurate under large disturbances. The behavior of VSG-IIDGs under large disturbances should be explored with nonlinear models of large signal stability. Large signal nonlinear analysis has wider application range and higher reliability. Besides, the stability domain can be obtained via the study of large signal stability. A large signal stable system is small signal stable yet the opposite is not true [14].
There are few studies presently on the large signal stability of VSG-IIDG due to various nonlinear factors [14]. Lyapunov-based methods are commonly used in the analysis of large signal stability. The main advantage of the methods is that they have a distinct computational performance. The stability assessment using Lyapunov theory is based on comparing the transient energy at fault clearing time with the critical energy without developing numerical computer models of the system [15]. Control models based on Lyapunov theory were proposed in references [16,17] for the integration of distributed generators into the distribution network. A Lyapunov function was established in [18] to prove the convergence of the proposed VSG control strategy. However, the transient stability mechanism and performance of VSG-IIDGs under large signal disturbances were not provided. Reference [19] investigated how transient energy of VSG-IIDG would be stored and released during disturbances using a Lyapunov function. Nevertheless, these studies address little on the assessment of large signal stability. The stability domain has not been depicted clearly. A valid Lyapunov function should be established to determine the stability domain. The mechanism of large signal stability, as well as the effect of parameters on the stability domain, needs to be further explored. Besides, further research is needed to establish the mathematical model of the VSG-IIDG for large signal nonlinear study. Reference [20] adopted the electrostatic machine model to establish the equivalent circuit of the droop-controlled IIDG. The Lyapunov function in [20] was not applicable to the stability assessment of the VSG-IIDG, but it provided reference values for the large signal stability study of the VSG-IIDG.
This paper focuses on large signal stability assessment of the VSG-IIDG. The contributions of the paper are threefold: (1) A nonlinear mathematical model of the VSG-IIDG is established by applying the equivalent model of an electrostatic machine. This nonlinear model combining both the electrical parts and control signals can be an analytical tool for the study of large signal stability (2). Based on Popov’s theory, a Lyapunov function is derived and the stability domain of IIDG is determined. This Lyapunov-based method has a distinct computational advantage. The stability assessment is based on comparing the transient energy of the postfault system with the critical energy without developing numerical computer models of the system (3). The large signal stability mechanism of VSG-IIDGs is revealed by analyzing the boundary of the stability domain. The boundary of the stability domain quantifies the magnitude of the deviation that the system can tolerate. The area of the stability domain reflects large signal stability in an intuitive way. The effect and sensitivity analysis of parameters on the stability domain are presented. The contributing factors of large signal stability are analyzed.
The sections of the paper are organized as follows: Section 2 analyzes the typical VSG control scheme of IIDG and establishes the equivalent electrostatic machine model. The nonlinear state matrix of the VSG-IIDG system is derived in Section 3. In Section 4, a Lyapunov function and the critical stability energy are figured out. The stability domain is defined accordingly. Simulation is conducted in Section 5 and the effect of parameters on the stability domain is analyzed.

2. Typical VSG Control Scheme and Its Equivalent Model for IIDG

2.1. The VSG Control Scheme of IIDG

The control scheme of the VSG-IIDG is shown in Figure 1. By controlling the switching pattern of the inverter, the VSG-IIDG has the dynamic properties of the synchronous generator (Figure 2). Given the fact that the capacity of the IIDG is much smaller than that of the host grid, the point of common coupling of VSG-IIDG is equivalent to an infinite bus (similar to single-machine-infinite-bus power system) and is regarded as a constant voltage source. In order to stress the key findings, the fluctuation of the renewable energy is ignored and the voltage at the DC-link capacitance is supposed to be constant. Besides, in order to simplify the analysis, the local load is represented by resistors and inductances.
For conventional synchronous generators, the inertia and damping play a significant role in terms of stability. However, the inner angle frequency of IIDGs is determined by the control scheme. The P-ω controller implements frequency adjustment in the VSG control [21]. The control function of the P-ω controller is expressed as Equation (1):
{ 2 H d ω d t = P r e f P o k Δ ω Δ ω = ω ω r e f
The P-ω controller imitates the behavior of a conventional synchronous generator and keeps the VSG-IIDG track the frequency of the host network.

2.2. VSG-IIDG Equivalent Electrostatic Machine Model

The stability assessment for VSG-IIDG depends on the nonlinear models. Andrade proposed the idea of modeling the inverter as an electrostatic machine [22]. The concept of electrostatic machine establishes a direct relationship between the DC and AC side of the inverter. This model allows the further introduction of traditional Lyapunov function and facilitates the analysis.
As shown in Figure 3, the IIDG is modeled as an electrostatic machine which is supplied by the direct voltage Udc. The electrostatic machine produces an electric field that induces alternating charges in the armature circuit. “Self” and “mutual” capacitances are included. The rotating reference frame (DQ0) is applied to link dc-side and ac-side. The electric field magnitude and the speed of the rotation parameters are related with Udc and the synthesized voltage reference vector g . Briefly, the changes in VSG-IIDG are represented by the changes in charges. Further derivation in [22] takes IIDG as a current source i i with equivalent resistance R e q and admittance X e q relevant to parameters of the inverter:
{ i i = ω C d c   d q 2 U d c j ω C d c d d 2 U d c R e q = 2 d q d d ω C d c ( d q 2 + d d 2 ) X e q = d q 2 d d 2 ω C d c ( d q 2 + d d 2 )
We substitute this electrostatic machine model into the IIDG system. The equivalent circuit of the IIDG system is shown in Figure 4. This equivalent circuit combines both the electrical part and control signals. It is used as the model for further study on large signal stability.

3. The Nonlinear Mathematical Model of VSG-IIDG

The equivalent circuit of IIDG system can be described as follows, by using Equation (1):
{ 2 H d θ i g   2 d t + k d θ i g d t P r e f + P o = 0 P o = α sin ( θ i g + θ α ) + β
where, α , θ α and β are related to the parameters of the equivalent circuit in Figure 4. The specific expression is presented in the Appendix A. Equation (3) can be seen as the expression of Figure 4 in mathematical form. Integrating the electrical parts with control signals, Equation (3) offers a mathematical model for stability study.
Next, the equilibrium point θ i g s is calculated from (5) when the system runs in a zero-deviation state:
θ i g   s = arcsin ( P r e f β α ) θ α
The state equation is normally based on equilibrium points so as to characterize the motion state of the system and facilitate the analysis of deviation. Define state variable x:
x = [ x 1 x 2 ] = [ θ i g θ i g s Δ ω ]
Then the mathematical model of the IIDG system is transferred from the equilibrium point to origin point. Hence the state equation in matrix formulation is:
{ x ˙ = A x + b f ( σ ) σ = c T x f ( σ ) = α 2 H   [ sin ( σ + θ i g s + θ α ) sin ( θ i g s + θ α ) ]
where, A = [ 0 0 1 k / 2 H ] , b = [ 0 1 ] , c = [ 1   0 ] .
Equation (6) is the expression describing the relationship between the input and the state of the IIDG system. The order of state matrix A is one, hence, matrix A is a singular matrix. It has two characteristic roots. One is zero and the other is k / 2 H . The control parameter H and k are greater than zero, so k / 2 H is in the open left half-plane. In order to facilitate the analysis, Equation (6) can be transformed into the following form through a non-singular transformation by reducing the order of x:
{ x ˜ ˙ = k 2 H   x ˜ f ( σ ) ε ˙ = f ( σ ) σ = 2 H k x ˜ + 2 H k ε
where, x ˜ = [ 0 1 ] x = x 2 = Δ ω .
The VSG-IIDG can be seen as a nonlinear system shown in Figure 5 where the control signals are tuned through both the nonlinear governor and the linear governor. The nonlinear part corresponds to the nonlinear function f ( σ ) in Equation (6).
G ( s ) = c T ( s I A )   1 b = 1 s ( s + k 2 H )
In terms of the linear part, the transfer function G ( s ) is obtained as in Equation (8). The nonlinear function f ( σ ) and the transfer function G ( s ) are both important in the study of large signal stability.

4. Lyapunov Function Construction and the Stability Domain Determination

In the following sections, the direct method of Lyapunov is applied to study the large signal stability of the VSG-IIDG. The approach uses a Lyapunov function V(x) to estimate the stability domain of the postfault system. The precondition is that this Lyapunov function satisfies the Popov’s theory on stability. The stability domain is defined by an inequality of V(x) < M, where M is a constant representing the critical stability energy. This method quantifies the extent of the deviation that the system can tolerate and features a remarkable computational advantage by simply comparing V with the critical energy rather than conventional step-by-step methods [15].

4.1. Popov stability Criterion

According to Popov’s Theory, in the system as (7), there exists a finite region of asymptotic stability if a real number q 0 exists that:
Re ( 1 + q j ω ) G ( j w ) 0  
for all ω > 0.
Figure 6 shows a graph drawing G ( j w ) on the complex plane. The crux of the Popov condition is to find a straight line of slope 1/q passing through the origin such that the imaginary part of G ( j w ) lies categorically under it with points of tangency permitted. In Figure 6, l 1 , l 2 , l 3 and l 4 satisfy the Popov condition, whereas l 5 and l 6 do not. There can be many straight lines meeting the condition (the lines between l 1 and l 4 ) and the figure only shows four of them. The different choices of the straight line, or to be exact, the values of q are related to the conservativeness of the stability analysis.
Apply the Popov condition to the system whose state equation is shown as (6):
q k   2 H 1 ( k 2 H ) 2 + ω 2 0
The inequality holds when q 2 H / k , and we choose q = 2 H / k + k 2 / 4 H 2 as indicated in [23]. Then we can use the q as a parameter to construct a Lyapunov function.

4.2. Lyapunov Function

This paper applies Kalman’s algorithm [24] to construct the Lyapunov function. The algorithm uses the q as a parameter. Kalman has proved that the Lyapunov function is positive definite with the semi-negative definite derivative. The Lyapunov function so derived satisfies the Popov condition.
(a)
Define the function W ( ω ) as:
W ( ω ) = Re ( 1 + q j ω ) G ( j ω ) ψ ( j w ) ψ ( j w ) = k 3   8 H 3
where:
ψ ( s ) = det ( s I + k 2 H   ) = s + k 2 H
(b)
Factorize W ( ω ) as: W ( ω ) = θ ( j ω ) θ ( j ω ) . Hence θ ( j ω ) is:
θ ( j ω ) = k 3   8 H 3
(c)
Define the leading coefficient of the polynomial in ascending powers when j ω = z in θ ( j ω ) a r . Then define a vector u with its components being the coefficients of the polynomial r ψ ( z ) θ ( z ) in ascending powers. Here the vector u reduces to a scalar in the VSG-IIDG model:
u = k 3   8 H 3
(d)
Figure out the symmetric positive definite matrix B ˜ by solving the Lyapunov matrix equation:
k 2 H   B ˜ B ˜ k 2 H = u u T
Here the Lyapunov matrix equation reduces to the scalar equation so that:
B ˜ = k 2   8 H 2
(e)
Get the Lyapunov function:
V ( x ˜ , ε , σ ) = x ˜ T B ˜ x ˜ + 1 2 q ε 2 + q 0 σ f ( σ )   d σ
To be more intuitive and convenient, the Lyapunov function needs to be transformed in terms of x. Hence the Lyapunov function of the VSG-IIDG is as follows:
V ( x 1 , x 2 ) = k 4 H   x 1 2 + ( H k + k 2 8 H 2 ) x 2 2 + x 1 x 2 + ( α k + α k 2 8 H 3 ) ( cos ( θ i g s + θ α ) cos ( x 1 + θ i g s + θ α ) x 1 sin ( θ i g s + θ α ) )

4.3. Lyapunov-Based Large Signal Stability Domain

The stability domain is defined by an inequality of form V(x) < M, where M is a constant representing the critical stability energy. With a Lyapunov function like (18), reference [25] uses Lagrange multipliers to define the critical stability energy M. Hence, the critical stability energy M of stability region is as Equation (19) in the system (6):
M = [ π 2 ( θ i g   s + θ α ) ] 2 c T b 1 c + q 0 π 2 ( θ i g s + θ α ) f ( σ ) d σ = [ π 2 ( θ i g s + θ α ) ] 2 2 ( k 4 32 H 4 + 4 H k 3 ) + ( α k + α k 2 8 H 3 ) [ 2 cos ( θ i g s + θ α ) sin ( θ i g s + θ α ) ( π 2 ( θ i g s + θ α ) ) ]
The inequality to find stability domain can be expressed as:
V ( x 1 , x 2 ) < M
Inequality (20) gives an effective way of determining the stability domain. If the operating points of postfault system satisfy inequality (20), it indicates that the VSG-IIDG is large-signal stable. The VSG-IIDG system can recuperate and return to the equilibrium stability where x 1 = x 2 = 0 . The complexity of the stability issue of VSG-IIDGs lies in two aspects. First, the VSG control is realized by power electronic devices, hence the features of inverters should be considered. Besides, the VSG emulates the behavior of a synchronous generator. The stability problems similar to conventional synchronous generators will definitely be involved [26]. Inequality (20) takes both of them into consideration. The impact of control parameters has been reflected in the coefficients. The stability domain can be depicted by solving the equation V ( x 1 , x 2 ) = M .
Table 1 presents the steps of determining whether a postfault IIDG system is large signal stable based on stability domain. First, the electrical and control parameters is given to establish the IIDG model. Then, the output power is obtained as Equation (27) and the equilibrium point is figured out according to Equation (4). The equilibrium point is used to derive the state equation and obtain b and c in Equation (6). Therefore, the Lyapunov function can be established as Equation (18), and the critical stability energy M of the stability region is calculated according to Equation (19). Then compare the value of the Lyapunov function with the critical M. Smaller function value indicates that the IIDG system is large signal stable.

5. Study Cases

5.1. Parameters of the Simulation Case

Simulations are performed in PSCAD/EMTDC to evaluate the large signal stability of the VSG-IIDG system. The topology and control scheme of the VSG-IIDG system is shown in Figure 1 and Figure 2. Parameters of the VSG-IIDG and the interconnected network are listed in Table 2. A solid three-phase line-to-ground fault of happens in the middle of the transmission line at 3 s.
According to Equation (4), the equilibrium point of VSG-IIDG is:
θ i g   s = 0.6756
Substitute the above data for Equation (6), and the state equation of the IIDG system yields as:
{ d x 1   d t = x 2 d x 2 d t = 0.133 f ( x 1 ) f ( x 1 ) = 1.2346 [ sin ( x 1 + 1.5442 ) sin ( 1.5442 ) ]

5.2. Large Signal Stability Analysis of the VSG-IIDG

The metallic short-circuit fault happens in the middle of the cable at the time of 3 s and is cleared after variable time delays. Once the fault is cleared, Δ ω and θ i g deviate from the equilibrium point at various degrees (different values of x 1 and x 2 ). It means operating points with different values of the Lyapunov function V ( x 1 , x 2 ) .
Figure 7 shows the contour plots of the Lyapunov function of operating points with different Δ ω and θ i g . Deviating from the equilibrium point at various degrees, the values of the Lyapunov function are diverse. For the points inside the stability domain, the value of the Lyapunov function is under the critical value M, which indicates that the system is asymptotic stability when operated in these points. Figure 8 shows that the stability domain is in a conical type. Whichever disturbance that has angular frequency and angle limited in the taper can recuperate the system.
Figure 8 shows the operational trajectories of five stable points (p1–p5) whose Lyapunov function value is less than the critical value M. The abscissa axis of time starts after the fault is cleared. These curves display the deviation of angular frequency and angle from a disturbance condition. It can be observed that after the fault is cleared, these operational trajectories tend to the equilibrium points where Δ ω = θ i g = 0 . Even though these operational trajectories deviate from the equilibrium point at various degrees, as time pass by, the system finally runs in a zero-deviation state.
When it comes to the unstable operating points, it is quite different. Take an unstable point m and a stable operating point n as example. The point m comes when the fault is cleared at 3.15 s and the deviation Δ ω = 11.5 , θ i g = 0.56 . While the point n comes when the fault is cleared at 3.1 s and the deviation Δ ω = 0.75 , θ i g = 0.1 . As shown in Figure 7, the Lyapunov function values of point m and n are on the different sides of the contour plane of the critical stability energy M. The Lyapunov function value of point n is smaller than the critical stability energy M, which indicated it is inside the stability domain. However, the Lyapunov function value of point m is greater than the critical stability energy M and it is an unstable point.
Figure 9 shows the comparison of simulations of the unstable point m and the stable operating point n. Due to different fault clearing time of these two points, the extent of disturbance varies (different values of Δ ω and θ i g ). As Figure 9a shows, the operational trajectory of the stable point n experiences damped oscillations before reaching the equilibrium stability where Δ ω = θ i g = 0 . The angular frequency gets slowly back to the reference value and the angle difference between the VSG-IIDG and the network reduces to zero. However, the trajectory of the unstable point m is divergent and of extreme volatility. Even when the fault is cleared, the angular frequency can’t keep in synchronism. The system can’t recuperate and tends to instability. The simulations shown in Figure 9 verify the analysis results.

5.3. The Impact of Parameters on the Stability Domain

The cable impedance, load power, and virtual inertia have a significant influence on the stability of VSG-IIDG [27,28]. The stability domain with different parameters will be explained in the next sections. The stability domain is determined by the boundary of the equation V(x) = M. The boundary quantifies the extent of disturbances that the VSG-IIDG system can endure. In the simulation, only one parameter is changed and the others are kept the same in different scenarios.
Figure 10a depicts the impact of the impedance of cable Z L ( R L + j X L ) on the stability domain. As impedance decreases, the stability domain expands in the ratio of equality. The increase in impedance means the connection between the IIDG and the network is less. When the supporting function of the network is weakened, the IIDG system is apt to instability.
Figure 10b shows the change of the stability domain when load power P D decreases. When load power decreases, the stability domain tends to increase. This is due to the fact that a larger load power adds more burden to the system. This is similar to the power angle stability of conventional synchronous generators, where lower power causes relative increase of the acceleration area and out-of-step situations will be rare.
The effect of different virtual inertia H on the stability domain is shown in Figure 10c. As the virtual inertia increases, the stability domain is prone to be smaller. What’s more, the shrink doesn’t happen in the ratio of equality. The shape of boundary changes and a peak appears.
Figure 10c shows that the rising virtual inertia has a negative effect on the large signal stability of the VSG-IIDG. In large signal stability of conventional synchronous generators, the change of inertia does not contribute to the change of the stability boundary [29]. However, it’s not the case in the VSG algorithm and Figure 10c witnesses a marked distortion of the stability domain as virtual inertia increases.
It should be noted that the rising virtual inertia will also slow down the response speed of the IIDG system since it can be seen as an integration constant. Hence the rising virtual inertia leads to a decrease not only in the distance between an operating point and the stability boundary but also the speed running from this point to the boundary. When the distance and speed are reduced simultaneously, it’s not sure whether the time duration will be shortened or not. That means, the critical cleaning time of the fault may not decrease even with a large virtual inertia (a smaller stability domain), since the response speed of the IIDG system is slow and there is still enough time to diagnose and handle the fault.
Table 3 presents the stability of the VSG-IIDG with different values of load power, cable impedance, and virtual inertia when the clearing time is 0.1 s. As shown in the table, when active power of local load increases, the IIDG system tends to instability. When load power and control parameter H remain unchanged, the IIDG system is unstable with larger cable impedance. Also, when the inertia is too large, the stability of the system goes from stable to unstable. The results show in Table 3 indicate that the IIDG system is apt to be unstable with larger load power and cable impedance or smaller virtual inertia.

5.4. Sensitivity Analysis of Parameters

Investigation of the large signal stability by ad hoc variations of the parameters is challenging, especially when several parameters act at the same time. Sensitivity analysis is helpful in identify which parameter should be modified in an easier way. Such studies are of high importance, considering the assessment of large signal stability in different scenarios with large expected variations in grid configurations, operating conditions, and system parameters.
This paper draws on the experience of transient stability analysis in conventional power system [30]. The sensitivity analysis of the load power, cable impedance, and virtual inertia on the stability domain area is performed. The definition of the sensitivity k s e n is in partial differential equation as:
k s e n = S d y
where, S d is the area of stability domain, y is the parameter to be analyzed (load power P D , cable impedance Z L , and virtual inertia H). The results are shown in Table 4, respectively.
It can be seen that the stability domain is mainly sensitive to load power, cable impedance, and virtual inertia. The sensitivity does not change a lot when different values of load power and cable impedance are adopted. When virtual inertia is kept at a small value, the area of stability domain changes slowly as inertia varies. It should be pointed out that only the virtual inertia can be modified to improve the stability of a control strategy. The load power and cable impedance cannot be influenced by the control. However, the load power and cable impedance can be selected within the proper range during the design of the system.

6. Conclusions

This paper assesses the large signal stability of the VSG-IIDG and derives the boundary of the stability domain:
(1)
The nonlinear mathematical model of the VSG-IIDG is established. The equivalent model of electrostatic machine is applied. Both the electrical parts and control signals are taken into account. This nonlinear model can be an analytical tool for the study of large signal stability.
(2)
A Lyapunov function is derived based on Popov’s theory to determine the stability of IIDG. By comparing the transient energy of the post-fault system with the critical energy, this Lyapunov-based method has a distinct computational advantage.
(3)
The stability domain is depicted and the large signal stability mechanism of VSG-IIDGs is revealed. The stability domain quantifies the magnitude of the deviation that the system can tolerate. The impacts and sensitivity analysis of parameters on the stability domain are presented. The results indicate that large disturbances may lead to instability of the VSG-IIDG with deviation and oscillation of angular frequency. The VSG-IIDG tends to instability with larger load power and cable impedance or smaller virtual inertia.
The study contributes to the further studies of large signal stability. The results have great significance for the design, operation, and planning of the VSG control scheme. What’s more, the analysis of transient process and mechanism of the VSG is helpful in further research on protection and fault analysis.

Author Contributions

Writing—original draft, data curation, formal analysis, M.L.; project administration, methodology, writing—review & editing, W.H.; writing—review & editing, validation, M.Y.; resources, supervision, N.T.

Funding

This research was funded by [National Key Research and Development Program of China] grant number [2017YFB0903202], [Shanghai Sailing Program] grant number [17YF1410200], [National Natural Science Foundation of China] grant number [51807117].

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

θ i the power angle of the IIDG
θ i g the power angle difference between θ i and θ g
θ i g s the equilibrium point
θ g the power angle of the common coupling point of VSG-IIDG
ωvirtual angular frequency
ω r e f angular frequency reference
Δ ω the angular frequency difference between ω and ω r e f
C a c , X C , Y C filtering capacitance of VSG-IIDG and the relative reactance and admittance
C d c capacitance of the DC link of VSG-IIDG
Hvirtual inertia of the VSG-IIDG
L L , R L , X L , Y L , Z L inductor, resistance, reactance, admittance, and impedance of the cable
L f , X f , Y f filtering inductor of VSG-IIDG and the relative reactance and admittance
Mcritical stability energy
P D , Q D , R D , X D , Y D active power, reactive power, equivalent resistance, reactance, and admittance of the local load
Pooutput active power measurement of the IIDG
P r e f active power reference of the IIDG
R e q , X e q , Y e q the equivalent resistance, reactance, and admittance of the equivalent IIDG current source
S d the area of stability domain
U d c voltage at the DC link of VSG-IIDG
U g voltage at the point of common coupling
d d , d q average values of voltage modulation for the duty cycle in d and q axis
g the synthesized voltage reference vector in VSG control
kthe damping factor of VSG-IIDG
k s e n the sensitivity
qPopov coefficient

Appendix A

In Figure 4, i i and i g are the equivalent current source of the IIDG and the DN. θ i g is the difference of the angle:
{ i i = i i θ i i g = U g Y L = i g θ g θ i g   = θ i θ g
Y f and Y C are the admittances of LC filter. Y L is the admittance of transmission line. Y D and Y e q are the equivalent admittance of load and IIDG. The expressions are:
{ Y f = ( j X f ) 1 Y C = ( j X C ) 1 Y L = ( R L + j X L ) 1 Y D = ( R D + j X D ) 1 Y e q = ( R e q + j X e q ) 1
The output voltage of IIDG V can be expressed as:
V = Y m i g + Y n i i  
where: { Y m = Y f + Y e q Y f Y L + Y L Y e q + Y f Y C + Y C Y e q + Y f Y D + Y D Y e q + Y f Y e q Y n = Y f Y f Y L + Y L Y e q + Y f Y C + Y C Y e q + Y f Y D + Y D Y e q + Y f Y e q .
Hence, the output power of the VSG-IIDG P o in Figure 5 can be expressed as:
P o = α sin ( θ i g   + θ α ) + β
where: { α = Y m * Y n Y L * Y m Y n * Y L * + Y n i i i g θ α = π 2 ( Y m * Y n Y L * Y m Y n * Y L * + Y n ) * β = Re [ ( Y m 2 Y L * + Y m 2 Y D * Y m ) i g 2 + Y n 2 ( Y L * + Y D * ) i i 2 ] .

References

  1. Shuai, Z.; Sun, Y.; Shen, Z.J.; Tian, W. Microgrid stability: Classification and a review. Renew. Sustain. Energy Rev. 2016, 58, 167–179. [Google Scholar] [CrossRef]
  2. Soni, N.; Doolla, S.; Chandorkar, M.C. Improvement of Transient Response in Microgrids Using Virtual Inertia. IEEE Trans. Power Deliv. 2013, 28, 1830–1838. [Google Scholar] [CrossRef]
  3. Liu, J.; Miura, Y.; Bevrani, H. Enhanced Virtual Synchronous Generator Control for Parallel Inverters in Microgrids. IEEE Trans. Smart Grid 2017, 8, 2268–2277. [Google Scholar] [CrossRef]
  4. Chen, L.; Wang, R.; Zheng, T. Optimal Control of Transient Response of Virtual Synchronous Generator Based on Adaptive Parameter Adjustment. Proc. CSEE 2016, 11, 5724–5731. [Google Scholar]
  5. Zhang, X.; Zhu, D.; Xu, H. Review of Virtual Synchronous Generator Technology in Distributed Generation. J. Power Supply 2012, 10, 1–6. [Google Scholar]
  6. Mehrasa, M.; Godina, R.; Pouresmaeil, E.; Vechiu, I.; Rodríguez, R.L.; Catalão, J.P. Synchronous active proportional resonant-based control technique for high penetration of distributed generation units into power grids. Presented at the IEEE Pes Innovative Smart Grid Technologies Conference Europe, Torino, Italy, 26–29 September 2017; pp. 1–6. [Google Scholar]
  7. Pouresmaeil, E.; Mehrasa, M.; Godina, R. Double synchronous controller for integration of large-scale renewable energy sources into a low-inertia power grid. Presented at the IEEE Pes Innovative Smart Grid Technologies Conference Europe, Torino, Italy, 26–29 September 2017; pp. 1–6. [Google Scholar]
  8. D’Arco, S.; Suul, J.A.; Fosso, O.B. A virtual synchronous machine implementation for distributed control of power converters in Smart Grids. Electr. Power Syst. Res. 2015, 122, 180–197. [Google Scholar] [CrossRef]
  9. Meng, J.; Wang, Y.; Shi, X. Control strategy and parameter analysis of distributed inverter based on virtual synchronous generator. Trans. China Electrotech. Soc. 2014, 29, 1–10. [Google Scholar]
  10. Du, Y.; Su, J.; Zhang, Z. A mode adaptive FM control method for microgrid. China J. Electr. Eng. 2013, 33, 67–75. [Google Scholar]
  11. Mehrasa, M.; Pouresmaeil, E.; Zabihi, S.; Vechiu, I.; Catalão, J.P. A multi-loop control technique for the stable operation of modular multilevel converters in HVDC transmission systems. Int. J. Electr. Power Energy Syst. 2018, 96, 194–207. [Google Scholar] [CrossRef]
  12. Cheng, C.; Yang, H.; Zeng, Z. Adaptive control method of rotor inertia for virtual synchronous generator. Autom. Electr. Power Syst. 2015, 19, 82–89. [Google Scholar]
  13. Huang, L.; Xin, H.; Huang, W. Quantitative analysis method of frequency response characteristics of power system with virtual inertia. Autom. Electr. Power Syst. 2018, 42, 31–38. [Google Scholar]
  14. Zhu, S.; Liu, K.; Qin, L. A review of transient stability analysis for power electronic systems. Proc. Chin. Acad. Electr. Eng. 2017, 37, 3948–3962. [Google Scholar]
  15. Kabalan, M.; Singh, P.; Niebur, D. Large Signal Lyapunov-Based Stability Studies in Microgrids: A Review. IEEE Trans. Smart Grid 2017, 8, 2287–2295. [Google Scholar] [CrossRef]
  16. Mehrasa, M.; Adabi, M.E.; Pouresmaeil, E. Direct Lyapunov control (DLC) technique for distributed generation (DG) technology. Electr. Eng. 2014, 96, 309–321. [Google Scholar] [CrossRef]
  17. Mehrasa, M.; Pouresmaeil, E.; Catalao, J.P.S. Direct Lyapunov Control Technique for the Stable Operation of Multilevel Converter-Based Distributed Generation in Power Grid. Emerg. Sel. Top. Power Electr. 2014, 2, 931–941. [Google Scholar] [CrossRef]
  18. Liu, Y.; Chen, J.; Hou, X. Dynamic frequency stabilization control strategy for microgrid based on the adaptive virtual inertial system. Autom. Electr. Power Syst. 2018, 42, 75–82. [Google Scholar] [CrossRef]
  19. Alipoor, J.; Miura, Y.; Ise, T. Power system stabilization using a virtual synchronous generator with alternating moment of inertia. IEEE J. Emerg. Sel. Top. Power Electr. 2015, 3, 451–458. [Google Scholar] [CrossRef]
  20. Andrade, F.; Kampouropoulos, K.; Romeral Martínez, J.L.; Vasquez Quintero, J.C.; Guerrero, J.M. Study of large-signal stability of an inverter-based generator using a Lyapunov function. Presented at the 40th Annual Conference of the IEEE Industrial Electronics Society, Dallas, TX, USA, 29 October–1 November 2014; pp. 1840–1846. [Google Scholar]
  21. Gao, F.; Iravani, M.R. A Control Strategy for a Distributed Generation Unit in Grid-Connected and Autonomous Modes of Operation. IEEE Trans. Power Deliv. 2008, 23, 850–859. [Google Scholar]
  22. Andrade, F.A.; Romeral, L.; Cusido, J. New model of a converter-based generator using electrostatic synchronous machine concept. IEEE Trans. Energy Convers. 2014, 29, 344–353. [Google Scholar]
  23. Pai, M.A.; Mohan, M.A.; Rao, J.G. Power System Transient Stability: Regions Using Popov’s Method. IEEE Trans. Power Appl. Syst. 1970, PAS-89, 788–794. [Google Scholar] [CrossRef]
  24. Kalman, R.E. Liapunov functions for the problem of Lur’e solving (8b): In automatic control. Proc. Natl. Acad. Sci. USA 1963, 49, 201–205. [Google Scholar] [CrossRef] [PubMed]
  25. Walker, J.A.; McClamroch, N.H. Finite regions of attraction for problem of Lur’e. Intern. J. Control 1967, 6, 331–336. [Google Scholar] [CrossRef]
  26. Wang, X.; Jiang, H.; Liu, H. Summary of the research on virtual synchronous generator grid connection stability. North China Power Technol. 2017, 9, 14–21. [Google Scholar]
  27. Renedo, J.; GarcíA-Cerrada, A.; Rouco, L. Active power control strategies for transient stability enhancement of AC/DC grids with VSC-HVDC multi-terminal systems. IEEE Trans. Power Syst. 2016, 31, 4595–4604. [Google Scholar] [CrossRef]
  28. Hu, T. A nonlinear-system approach to analysis and design of power-electronic converters with saturation and bilinear terms. IEEE Trans. Power Electr. 2011, 26, 399–410. [Google Scholar] [CrossRef]
  29. He, Y.; Wen, Z.; Wang, F. Analysis of Power System; The Second Volume; Huazhong University of Science and Technology Press: Wuhan, China, 1985. [Google Scholar]
  30. Fouad, A.A.; Vittal, V. Power System Transient Stability Analysis Using the Transient Energy Function Method; Prentice Hall: Upper Saddle River, NJ, USA, 1992. [Google Scholar]
Figure 1. The control scheme of VSG-IIDG.
Figure 1. The control scheme of VSG-IIDG.
Energies 11 02273 g001
Figure 2. Synchronous generator model.
Figure 2. Synchronous generator model.
Energies 11 02273 g002
Figure 3. Representation of the electrostatic machine.
Figure 3. Representation of the electrostatic machine.
Energies 11 02273 g003
Figure 4. The equivalent circuit of the IIDG system.
Figure 4. The equivalent circuit of the IIDG system.
Energies 11 02273 g004
Figure 5. Nonlinear VSG-IIDG model.
Figure 5. Nonlinear VSG-IIDG model.
Energies 11 02273 g005
Figure 6. G ( j w ) satisfying the Popov condition on the complex plane.
Figure 6. G ( j w ) satisfying the Popov condition on the complex plane.
Energies 11 02273 g006
Figure 7. Contour plots of V function.
Figure 7. Contour plots of V function.
Energies 11 02273 g007
Figure 8. Operational trajectories of stable points.
Figure 8. Operational trajectories of stable points.
Energies 11 02273 g008
Figure 9. Comparison of simulations when fault clearing time are different. (a) Operational trajectories; (b) Angle difference; (c) Angular frequency of the IIDG.
Figure 9. Comparison of simulations when fault clearing time are different. (a) Operational trajectories; (b) Angle difference; (c) Angular frequency of the IIDG.
Energies 11 02273 g009
Figure 10. The stability domain. (a) Under different impedance of cable; (b) Under different load power; (c) Under different virtual inertia.
Figure 10. The stability domain. (a) Under different impedance of cable; (b) Under different load power; (c) Under different virtual inertia.
Energies 11 02273 g010
Table 1. Rules and steps of the stability domain determination.
Table 1. Rules and steps of the stability domain determination.
Step 1Input the parameters: voltage, current, control parameters as H, k, etc.
Step 2Obtain the expression of output power as Equation (27)
Step 3Calculate the equilibrium point as Equation (4)
Step 4Establish Lyapunov function V( x 1 , x 2 ) as Equation (18)
Step 5Figure out the critical stability energy M as Equation (19)
Step 6If V( x 1 , x 2 ) < M, the system is large signal stable
Table 2. Parameters of the VSG-IIDG and control system.
Table 2. Parameters of the VSG-IIDG and control system.
ParametersValue
DC voltage U d c 1 kV
DC capacitance C d c 100 μF
Filtering capacitance C a c 400 μF
Filtering impedance L f 1 mH
AC voltage U g 311 V
Resistance of cable R L 0.2 Ω
Reactance of cable L L 1 mH
Power of load P D + j Q D 300 kW + j100 kVar
Virtual inertia H0.15
Damping factor k0.01
Table 3. Stability of the VSG-IIDG with different P D , Z L , and H.
Table 3. Stability of the VSG-IIDG with different P D , Z L , and H.
Load Power P D ( MW ) Cable Impedance Z L ( Ω ) Virtual Inertia HStability
0.30.2 + j0.3140.15stable
0.60.2 + j0.3140.15stable
10.2 + j0.3140.15unstable
0.30.5 + j0.7850.15stable
0.31 + j1.570.15unstable
0.30.2 + j0.3140.5stable
0.30.2 + j0.3142unstable
Table 4. Sensitivity of the stability domain area to P D , Z L , and H.
Table 4. Sensitivity of the stability domain area to P D , Z L , and H.
P D (MW)0.10.30.50.70.9
Sensitivity−4.991−4.981−4.701−4.687−4.295
Z L (Ω)0.2 + j0.3140.4 + j0.6280.6 + j0.9420.8 + j1.2561 + j1.57
Sensitivity−5.661−5.508−5.138−4.037−4.011
H0.10.511.52
Sensitivity−1.980−2.256−3.007−4.751−5.124

Share and Cite

MDPI and ACS Style

Li, M.; Huang, W.; Tai, N.; Yu, M. Lyapunov-Based Large Signal Stability Assessment for VSG Controlled Inverter-Interfaced Distributed Generators. Energies 2018, 11, 2273. https://doi.org/10.3390/en11092273

AMA Style

Li M, Huang W, Tai N, Yu M. Lyapunov-Based Large Signal Stability Assessment for VSG Controlled Inverter-Interfaced Distributed Generators. Energies. 2018; 11(9):2273. https://doi.org/10.3390/en11092273

Chicago/Turabian Style

Li, Meiyi, Wentao Huang, Nengling Tai, and Moduo Yu. 2018. "Lyapunov-Based Large Signal Stability Assessment for VSG Controlled Inverter-Interfaced Distributed Generators" Energies 11, no. 9: 2273. https://doi.org/10.3390/en11092273

APA Style

Li, M., Huang, W., Tai, N., & Yu, M. (2018). Lyapunov-Based Large Signal Stability Assessment for VSG Controlled Inverter-Interfaced Distributed Generators. Energies, 11(9), 2273. https://doi.org/10.3390/en11092273

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