Next Article in Journal
Additive Manufacturing and Performance of E-Type Transformer Core
Next Article in Special Issue
Sizing Procedure for System Hybridization Based on Experimental Source Modeling in Grid Application
Previous Article in Journal
Inhomogeneities and Cell-to-Cell Variations in Lithium-Ion Batteries, a Review
Previous Article in Special Issue
Comparative Study of Hysteresis Controller, Resonant Controller and Direct Torque Control of Five-Phase IM under Open-Phase Fault Operation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparative Analysis of a Power System Stability with Virtual Inertia

Department of Electrical Power Engineering, Czech Technical University, 166 27 Prague, Czech Republic
*
Author to whom correspondence should be addressed.
Energies 2021, 14(11), 3277; https://doi.org/10.3390/en14113277
Submission received: 6 April 2021 / Revised: 21 May 2021 / Accepted: 1 June 2021 / Published: 3 June 2021

Abstract

:
The paper investigates the stability of a power system with synchronverters. A synchronverter is a control strategy for voltage source converters that introduces virtual inertia by mimicking synchronous machines. The authors picked a commonly known IEEE 9 bus and IEEE 39 bus test case systems for the test case studies. The paper presents the power system’s modal analysis with Voltage Source Converters (VSCs) controlled as synchronverters, vector control, or Rate of Change of Frequency-based Virtual Synchronous Generator, thus comparing different approaches to VSC control. The first case study compares selected control algorithms, the IEEE 9 bus system, with one VSC in the paper. The results demonstrate the benefits of synchronverters over other control strategies. The system with synchronverters has a higher minimal damping ratio, which is proven to be the case by numerical simulations. In the second case study, the effects of virtual inertia placement were investigated. The computations showed that placement is indeed important, however, the control strategy is as important. Besides, the system with synchronverters exhibits better stability characteristics. The paper demonstrates that the application of synchronverters is feasible and can meet the demand for algorithms that bring the benefits of virtual inertia.

1. Introduction

Power systems have dramatically changed in the last decade. The penetration of Renewable Energy Sources (RESs) combined with plans to phase out nuclear power plants in many countries shapes the industry. Those developments brought new dilemmas to power system operators and planners. A major one is how to deal with low grid inertia. There is no doubt that inertia plays a significant role in short-term frequency control. Since it naturally stores energy, which is used during a spike in a power imbalance in the gird. The importance of inertia for frequency regulation is demonstrated in Figure 1.
Therefore, many people in the industry are trying to solve the upcoming problem with control of a low inertia grid. Some researchers proposed to change the requirements on Rate of Change of Frequency (ROCOF) as a simple solution. This approach was adopted in Ireland [1]. There are other proposals, one that stands out the most is an introduction of virtual inertia. For instance, Hydro-Québec started to worry about low inertia in 2003 and set up requirements for power plant operators [2]. Similar ideas were proposed across the Atlantic by ENTSO-E in 2017 [3]. However, this document does not set mandatory requirements but instead introduces National Power System operators’ guidelines in Europe. Furthermore, it outlines some examples of Canada that already require frequency response by wind power plants. Reference [2] probably preceded others in setting up virtual inertia guidelines.
The most cost-effective solution to a low inertia grid is to make RES power plants behave like synchronous generators. The fundamental idea is to control VSC so it will respond to a disturbance like a synchronous generator. Intensive research of that topic yielded many different control strategies. The authors primarily divide them into two categories grid following and grid forming. The grid-following control for synchronization requires a stiff grid reference. Vector control is typical grid-following control that requires Phase Locked Loop (PLL) for synchronization with the grid [5]. Different proposed control strategies emulate inertia response, however, some still do it through PLL-provided measurement of frequency [6,7]. The authors call such control strategies based on frequency measurement Rate of Change of Frequency-based Virtual Synchronous Generators or in abbreviation ROCOF VSG. Some of such proposed topologies [8] utilize simple power-based calculations. Others build on vector control and introduce voltage reference signal that emulate inertia response [9,10]. The other control approach is to emulate synchronous generators and make a VSC essentially a grid-forming device. Synchronous generator based topologies are synchronverters [11,12] and Virtual Synchronous Machine (VISMA) [13,14]. The advantage that emulation of inertia brings that it becomes a tunable parameter in the system. The rotating mass of the rotor gives the inertia of the generator. The energy storage limits virtual inertia on the other hand that Renewable Energy Source has. Usually, it is a capacitor or a battery.
Existing research papers mainly focused on analyzing the stability of a synchronverter connected to simple networks without comparing them to other common control strategies. For instance an interesting studies [15,16] on the application of synchronverters to damp oscillations grid did not discuss other possible control algorithms that can be an alternative to synchronverters. Additionally, in these papers, the dynamics of the DC bus is not considered. Furthermore, paper [17] investigates the interaction between grid following and grid forming devices. However, the test case system is focused on the stability of power converters, not on the power system dynamics with several generators. Thus, the main contribution of the paper to the field of power system stability is the following. The research focuses on the dynamics and stability of a power system with multiple generators. The paper compares different VSC control strategies that represent the main categories outlined by the authors. The chosen control strategies are vector control, ROCOF VSG [9] and synchronverters [11]. The paper presents the comparative analysis of the investigated control topologies in IEEE 9 bus and IEEE 39 bus benchmark systems. For the comparison, the authors conducted modal analysis and numerical simulations to verify the conclusions. In the first case, the authors analyze the IEEE 9 bus system’s stability with only one Photovoltaic Power Plant. In the second case where the test case system is large, authors investigate the effects of placement of VSCs with different control strategies. The conclusion emerged from the placement analysis that the system’s performance depends not only on the placement of virtual inertia but on the control strategy that implements it.

2. Mathematical Models of the Devices

The section provides mathematical models of the devices connected to the power system. Furthermore, it explains some assumptions that have been made in the modeling. Additionally, it briefly describes how the devices can be connected to the network.

2.1. Synchrnous Generator: Two-Axis Model

The two-axis model or 4th order model is used in the power system model described in the book [18]. This model closely enough approximates the generator’s behavior even though it neglects the dynamics of one damping winding in d and one in q axis. Furthermore, the stator transients are neglected [19]. This introduces some errors, however, the results are more conservative. Thus it is not a problem for stability studies. For the ith generator, the model can be stated as follows.
T d 0 , i d E q , i d t = E q , i ( x d , i x d , i ) i d , i + E f d , i
T q 0 , i d E d , i d t = E d , i + ( x q , i x q , i ) i q , i
d δ i d t = ω 0 ( ω i ω s )
J i d δ i d t = M m e c h , i E q , i i q , i E d , i i d , i ( x d , i x d , i ) i d , i i q , i T F W
Stator algebraic equations:
E q , i = R s , i i q , i + x d , i i d , i + v q , i
E d , i = R s , i i d , i x q , i i q , i + v d , i
E q , i , E d , i transient voltage in d q axis
E f d , i field voltage
i q , i , i d , i current in d q axis
J i inertia constant
w i generator angular velocity
δ i load angle
w 0 nominal angular velocity
w s synchronous angular velocity
T F W damping torque to emulate the effect of damping windings

2.2. VSC Electric Part Model

Hereafter the circuit equations of the Voltage Source Converter are stated. Figure 2 shows the converter’s line diagram with an L-type filter. A further assumption is to neglect DC side switching losses since they are not relevant for grid dynamics/interaction. Thus it simplifies the dynamics of the DC side of the VSC.
m d v d c 2 = L f d i d , i d t ω s L f i q , i + R f i d , i + v d , i
m q v d c 2 = L f d i q , i d t + ω s L f i d , i + R f i q , i + v q , i
C d c d v d c d t = p p v p v s c v d c
p v s c = ( m d v d c i d + m q v d c i q ) / 2
L f , i filter inductance
R f , i filter resistance
m q , i , m d , i modulation vector in d q
C d c DC capacitor
v d c DC voltage
i q , i , i d , i injected current in d q axis
p p v power produced by PV panels
p v s c power injected to the grid
v d , i , v q , i voltage in PCC
Figure 2. VSC schematic diagram.
Figure 2. VSC schematic diagram.
Energies 14 03277 g002

2.3. Synchronverter Model in d q Reference Frame

The paper [11] introduced the notion of a synchronverter, a converter that mimics a synchronous generator. Figure 3 shows the block diagram of the synchronverter. This section provides the mathematical model in d q frame of a synchronverter connected to the grid.
The main idea behind the synchronverter control strategy is that it has, to some degree, mimic the behavior of a synchronous machine [20]. Hence, the mathematical model of synchronverters in the d q frame should be similar to the synchronous generator equations. Assuming that the virtual machine has only three winding in d q frame:
v d = R d i d + d ψ d dt ω ψ q
v q = R q i q + d ψ q dt + ω ψ d
v f = R f i f + d ψ f dt
The relations between fluxes and currents are following:
ψ d ψ q ψ f = L d d 0 L f d 0 L q q 0 L d f 0 L f f · i d i q i f
where L f d is the magnetic coupling inductance between filed circuit and d stator axis. Furthermore, in synchronverters, the assumption is that there is only one-way magnetic coupling i.e., from rotor t stator. Therefore, L d f = 0 and L d d = L f , L f being leakage inductance in a real machine or it is filter inductance, in the case of the synchronverter. Additionally, the derivative of the field current is assumed to be 0. With the assumption that R d = R q = R s R s is filter resistance Equations (11)–(13) can be simplified as follows
v d = R s i d + ω r L f i q L f d i d dt
v q = R s i q ω r L f i d L f d i q dt + ω r ψ f
where ψ f = L f d i f
Clearly (15) and (16) resemble the equations of the VSC electrical part (7) and (8). Thus modulations signals become m d = 0 and
m q = ω r ψ f
The swing equation emulating the mechanical part:
J d ω r d t = p m e c h ω r p s y n c h ω r D p ( ω r ω n )
where power injected by synchronverter p s y n c h in a synchronverter rotating frame is calculated considering assumptions (10) and (17). Thus,
p s y n c h = ω r ψ f v d c 2 i q r
Virtual mechanical power input— p m e c h is set to control the DC bus voltage [21,22]. Thus, the equations for the controller are:
p m e c h = K p , p ( v d c v d c , r e f ) + p r e g , p i
d p r e g , p i d t = K i , p ( v d c v d c , r e f )
where K p , p and K i , p that are proportional and integral constants respectively
The virtual load angle of the synchronverter is computed similarly to a synchronous machine:
d δ d t = ω r ω s
As Figure 3 shows that in PCC the synchronverter can control reactive power injection and/or voltage. Assuming that the synchronverter is set to be P-V node then K i , q = 0
d ψ d t = K i , v ( v r e f v p c c ) + K i , q ( q r e f q )
where v p c c is the voltage amplitude in PCC and K i , v is integral constant of voltage regulator.

2.4. Vector Control of VSC

Voltage Source Converters are usually equipped with vector control. Vector control block diagram is shown in Figure 4. This control strategy allows decoupling the control of active and reactive power injections using Park transformation (A1). In power system transient stability studies, the effects of some parts are neglected, for example, losses in the converter. The mathematical model is provided here in addition to the equations of the converter. Here the calculation of modulation of voltage vector is provided. It should be noted that the control objective in the authors’ models is the voltage in PCC, not the reactive power injection. Furthermore, the chosen model utilizes a PI regulator for DC voltage control and an Integral (I) regulator for PCC voltage. Additionally, for the voltage shift’s compensation due to filter resistance and reactance, the feedforward shift is used for modulation [23]. Synchronization with the grid is provided by synchronous reference frame PLL (SRF-PLL). A thorough derivation of equations is not necessary in this case since it is a commonly known control strategy.
m d = v d ω s L f i q + K p , i d ( i d , r e f i d ) + M d
m q = v q + ω s L f i q + K p , i q ( i q , r e f i q ) + M q
d M d d t = K i , i d ( i d , r e f i d )
d M q d t = K i , i q ( i q , r e f i q )
d p r e g , p i d t = K i , d c ( v d c v d c , r e f )
d i d , r e f d t = K i , p ( p r e g , p i + K p , d c ( v d c v d c , r e f ) p v s c )
d i q , r e f d t = K i , v ( v r e f v p c c )
d I p l l d t = k i , p l l ( v q c v q r e f )
d ω p l l d t = I p l l + k p , p l l ( v q c v q r e f )
d θ p l l d t = ω p l l ω s
Figure 4. Vector control block diagram.
Figure 4. Vector control block diagram.
Energies 14 03277 g004
Note, currents regulator constants are the same, thus K p , i d = K p , i q and K i , i d = K i , i q .

2.5. ROCOF VSG

The third examined strategy is somewhat between vector control and synchronverters. It also provides an inertia response to disturbance. However, the inertia control response is built on top of vector control. There are different topologies that emulate inertia response by using frequency measurement. Authors name such control approach a Rate of Change of Frequency Virtual Synchronous Generator or in abbreviation ROCOF VSG. This particular topology was proposed in the paper [9]. Figure 5 shows the block diagram of the examined control strategy. Fundamentally the mathematical model of this topology is similar to vector control with only one difference. The DC bus voltage reference consists of two signals. The first is the voltage bus reference, and the second is the virtual inertia response signal that modifies the DC bus reference. Thus to build mathematical model of a VSC control as ROCOF VSG in the Equations (25)–(33) DC bus voltage reference signal is going to be:
v d c , r e f = v d c , n o m i n a l + V J , r e f + D p ( ω p l l ω n )
Figure 5. ROCOF VSG (Rate of Change of Frequency Virtual Synchrnous Generator) block diagram.
Figure 5. ROCOF VSG (Rate of Change of Frequency Virtual Synchrnous Generator) block diagram.
Energies 14 03277 g005
The control signal of virtual inertia response emulation is as follows:
d V J , r e f d t = ( H p ω ˙ p l l V J , r e f ) / T j
where H p emulate inertia and D p damping constants. The time delay T j simulates the gradual power injection from the kinetic energy of a synchronous generator. Thus, ROCOF VSG provides a power response based on frequency change. It should be noted that this control strategy is not grid forming in itself. It does support frequency but still uses PLL to track the frequency of the grid. Therefore, ROCOF VSG is very sensitive to the ability of PLL to track frequency correctly without oscillations.

2.6. Network Model

In transient stability studies, conventionally, the dynamics of lines are neglected. Since in comparison to time constants of machines, the propagation in lines happens almost instantly. Therefore the network equations become algebraic equations. The general form of a power system mathematical model has the following form:
d x d t = f ( x , y , t )
0 = g ( x , y , t )
where x is the state vector, which depending on model, may contain variables of current injections, excitation controllers, turbine governors etc. Additionally, y vector usually corresponds to voltages.
The general form of network algebraic equations that describe power flow is following
Y · V I ( x , V ) = 0
where I ( x , V ) respects variable impedance loads. However, for simplifications often the constant impedance loads are used in stability studies [18,24,25]. Thus, equations can be simplified
Y · V = I ( x )
Figure 6 shows the relation between different reference frames. In the authors’ implementation, all devices are modeled in their own rotating frame for more modular code, and the current injections are transformed using Park transformation Equation (A1) (Appendix A). The transformation between different d q rotating frames can be made directly by simplifying transformation matrices. Let’s examine ith generator with angular velocity w r , i . The network has synchronous velocity ω s , sometimes it is calculated as a fictitious center of inertia reference or fixed to one of the generators. Consequently, machine electrical angle is θ r = ω r d t or θ r = ω s t + δ . As Figure 6 shows to transform a vector from machine d q reference frame to two network d q reference frame. Firstly, the Inverse of the park transformation Equation (A1) where angles are machine electrical angle must be applied. Then from a b c natural frame back to network d q frame, the park transformation Equation (A1) with grid electrical angle is applied. It can be written more strictly in matrix form:
f g r i d = T d q , g r i d · T d q , i 1 f r o t o r
Assuming that zero component can be neglected, i.e., the symmetrical network:
f d , i g r i d f q , i g r i d = cos ( δ i ) sin ( δ i ) sin ( δ i ) cos ( δ i ) · f d , i r o t o r f q , i r o t o r

3. Stability Comparison

In this section, the authors compare a power system’s stability with a photovoltaic power plant controlled as a synchronverter, or vector control strategy, or ROCOF VSG. Additionally, a brief description of the methodology that was used for comparison is provided. For the case studies, the authors picked a commonly know IEEE 9 bus and IEEE 39 bus test systems [26]. IEEE 9 bus is complex enough but still can be easily understood. Therefore, for the IEEE 9 bus test system, the authors chose to carry out modal analysis and numerical simulations. The modal analysis provides an insight into system dynamics, and transient stability verifies whether the system can remain stable after big disturbances such as faults. For the case study of VSCs in the IEEE 39 bus test system, authors using modal analysis compare the quality of control for different configurations of the network. This test system allows demonstration of the effects of placement of virtual inertia since it is much larger. The mathematical models of the devices were implemented in Wolfram Mathematica [27]. Wolfram Mathematica combines powerful symbolic capability with all possible numerical tools for solving algebraic differential equations.

3.1. Modal Analysis

For small-signal stability analysis the general form of power system mathematical model (36) and (37) is linearized. Hence, the linear form of the model is as follows:
Δ x ˙ = A · Δ x + B · Δ v
where Δ x is the state vector of the power system, A is the system matrix corresponding to devices constants and B is the input matrix that correlates with node voltages. Algebraic equations describe the power flow in the network.
Δ i = Y N · Δ v
where Y N network admittance matrix and Δ i is the vector of current injections
The current injections of the devices are the following:
Δ i = C D · Δ x + D D · Δ v
where C D and D D is matrix corresponding to individual devices. Thus, algebraic equations can be simplified:
Y N · Δ v = C D · Δ x + D D · Δ v
Then the solution for the vector of voltages is following
Δ v = ( Y N D D ) 1 C D · Δ x
Hence, the system of differential-algebraic equations become the system of only differential equations:
Δ x ˙ = A · Δ x + B ( Y N D D ) 1 C D · Δ x
Consequently, A s y s = A + B ( Y N D D ) 1 C D is new system matrix with reduced algebraic equation. Thus, the small-signal stability can be analyzed by finding the eigenvalues of A s y s . Furthermore, the damping ratio ζ i of individual eigenvalues might be evaluated using the following formula:
ζ i = σ i σ i 2 + ω i 2
where the assumed form of the ith eigenvalue is λ i = σ i + j ω i . The eigenvalues of the system matrix obtained by previous computations are then used for comparing the stability of examined strategies. The Lyapunov indirect method allows gaining insights into the system dynamics without directly solving differential equations. The ability of a linear system to return to the steady-state after disturbance is defined by m a x ( R e ( λ ) ) . Furthermore, the damping ratio of eigenvalues predicts whether the oscillations provided by the eigenvalue persist only for a short time after disturbance. Fundamentally, it is a relation between exponential decay of the amplitude and the oscillations given by the complex part of the eigenvalue. Moreover, as part of the modal analysis, the authors calculated the participation factors for individual modes. A more thorough explanation can be found in Sauer’s book [18], or Kundur’s book [24]. The participation matrix is calculated as follows:
p k , i = | w k , i | | v k , i | k = 1 n | w k , i | | v k , i |
where w k , i and v k , i are components of left and right eigenvectors corresponding to ith eigenvalue. Inherently the sum of all participation factors for individual modes is equal to unity.

3.2. Case Study: IEEE 9 Bus

Western Electricity Coordinating Council three-machine system (Figure 7) is a classic test case system for stability studies. The model includes three machines, nine buses, and three loads. The authors chose to use constant impedance loads in the models. Furthermore, the machines are equipped with IEEE type 1 exciter and Steam turbine governor model (can be found in Sauer’s book [18]). Parameters of the system can be found in Appendix B.
In this test case, the PS with only one VSC connected to node 6 was analyzed. The PV power production was 0.6 pu. Firstly, the power system mathematical model was put together using equations in Section 2. After that, the equations were linearized, and the reduced system matrix was obtained. Computed eigenvalues of the system with synchronverters, vector control and ROCOF VSG are shown in Table 1, Table 2 and Table 3 correspondingly. Please note, only variables with a participation factor greater than 0.2 are displayed. Furthermore, the participation factor matrices are visualized in Figure A1, Figure A2 and Figure A3 (Appendix D). The eigenvalues with the smallest damping ratios ξ are highlighted in bold font in the tables. For the system with synchronverters, the damping ratio ξ m i n = 0.0298 , which was the best result of all examined strategies. Vector control had the lowest ξ m i n = 0.0062 damping ratio. The application of ROCOF VSG improved it, resulting in ξ m i n = 0.0112 . Figure 8 visualizes some of the results for demonstration of eigenvalues placements. In the case of synchronverters, most eigenvalues were closer to the x-axis, which means that the damping of these frequencies was higher than in the case of vector control or ROCOF VSG. Furthermore, the dashed lines (black for synchronverters, red for ROCOF VSG, and blue for vector control) show the lines corresponding to the smallest damping ratio. Additionally, in the case of synchronverters, the m a x ( R e ( λ ) ) was greater than for other control strategies. The maximum real part of the eigenvalues is a good indicator of the quality of control that provides information about the system’s ability to return to a steady-state after disturbance. For the ROCOF VCG and vector control, the eigenvalue with m a x ( R e ( λ ) ) was introduced by PLL.
The most severe disturbance that occurs in PS is a three-phase fault. Therefore, the authors ran a simulation with the following scenario. The three-phase fault occurred on line 8–9. After the fault was cleared, the faulted line was disconnected. Clearing times t c l = 0.1 s and t c l = 0.2 s were considered. The simulation results for the load angle of the third machine are presented in Figure 9 and Figure 10. For the t c l = 0.1 s, as Figure 9 shows vector control and ROCOF VSG had higher load angle deviation from steady-state value than synchronverters. Besides, the ROCOF VSG expectedly was slightly better than vector control. For t c l = 0.2 s, the system with vector control was unstable. ROCOF VSG allowed higher load angle deviation and had prolonged oscillations than synchronverters. It was expected because as Table 3 shows, ROCOF VSG had lower damping and greater m a x ( R e ( λ ) ) . Therefore, the overall dynamic stability of the power system was increased by applying synchronverter control. The synchronverter increased the system’s total inertia, thus generators sped up longer due to higher system inertia, which resulted in lower load angle deviation from steady-state value. Since vector control did not provide additional inertia the generator loses synchronism for short-circuit with t c l = 0.2 s. ROCOF VSG improved the system’s performance compared to vector control, but the resulting dynamics of the system had higher load angle deviations and oscillations than the system with synchronverters.

3.3. Case Study: IEEE 39 Bus

New England 10 machine system (Figure 11) is a well-known benchmark test system. This test case provides a fairly complex and close to the real power system model. In the model, two VSC were considered, and the effects of control strategy on the system’s stability were investigated. Parameters of the system can be found in Appendix C.
The primary focus of the case study was to analyze the effect of placement of VSCs in different nodes. The following scenarios were considered. Two PV power plants connected through VSCs were placed in different nodes, excluding nodes connected to generators. Each power plant produced 0.8 pu. The following control quality criteria were computed: as in the previous case i.e., minimal damping and m a x ( R e ( λ ) . Both allowed us to gain insights into the dynamics of the system without solving the differential equations numerically. Figure 12, Figure 13 and Figure 14 show that by using synchronverters in most cases the minimal damping ratio was improved. It should be noted that in some node pairs, the results are missing since this configurations had unstable eigenvalues, thus were excluded from calculation. The application of ROCOF VSG also improved damping but overall to a smaller degree than synchronverters.
Figure 15, Figure 16 and Figure 17 show the results of computation of m a x ( R e ( λ ) ) for the IEEE 39 bus system with different control strategies. This quality of control index shows that the synchronverter improved the system’s ability to return to steady-state. Furthermore, it also shows that indeed the placement of virtual inertia was important and different configurations of the PS could result in substantial differences in the index value. On the contrary, for vector control, the calculated values of the m a x ( R e ( λ ) ) were the same for all considered configurations. It should be noted that for ROCOF VSG, the results were similar to vector control. That is explained by PLL primarily dictating the results for these control strategies. To summarize, the quality of control not only depended on the placement of virtual inertia [29] but on the control algorithm that it implemented. The results presented in this paper show that the quality of control could be affected by the control strategy for virtual inertia as much as placement.

4. Conclusions

The paper presents a thorough comparison of power system stability with different control strategies of VSCs. In the paper, three different control strategies were compared using modal analysis and numerical simulations. The first case study results proved that the synchronverter provides better damping and lower m a x ( R e ( λ ) ) , which shows that it can reach steady-state faster. Furthermore, numerical simulations (Figure 9) demonstrated that Power System with synchronverters performs better during three-phase faults than PS with vector control or ROCOF VSG. ROCOF VSG also improves transient stability comparing to vector control and allows the system to withstand longer clearing times, but the synchronverter does it better, as results provided in Figure 10 show. In the second case study, the effects of placement were analyzed in IEEE 39 bus benchmark system. The results of the analysis yielded that for vector control and VSG ROCOF, the configuration of the network does not change m a x ( R e ( λ ) ) . However, for synchronverters, this parameter is very changed with the placement of VSCs. Furthermore, the presented results that synchronverters are better than other control strategies in almost every possible configuration in terms of damping and m a x ( R e ( λ ) ) . The overall conclusion is that the control strategy that provides virtual inertia is as important as the placement of VSCs.
To summarize, the presented cases demonstrate that the power system with synchronverters performs better than the one with vector control or ROCOF VSG. From the point of modal analysis, the stability is improved by the introduction of synchronverters. Furthermore, transient stability performance is indeed enhanced by the application of a synchronverter. In the authors’ opinion, synchronverters are a compelling control strategy that introduces virtual inertia to the system and meets the modern power system’s requirements. Additionally, virtual inertia becomes a tunable parameter that gives more flexibility to power system operators.

Author Contributions

Conceptualization, L.V. and Z.M.; methodology, L.V.; software, L.V.; validation, Z.M.; formal analysis, L.V. and Z.M.; investigation, L.V.; resources, Z.M.; writing—original draft preparation, L.V.; writing—review and editing, L.V.; visualization, L.V.; supervision, Z.M.; project administration, Z.M.; funding acquisition, Z.M. Both authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Czech Technical University in Prague grant number SGS20/ 165/OHK3/3T/13.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ACAlternating Current
DCDirect Current
VSCVoltage scource converter
VSGVirtual Synchronous Generator
PCCPoint of Common Coupling
PFParticipation Factor
PLLPhase Locked Loop
PSPower System
PVPhotovoltaic
ROCOFRate of Change of Frequency
SRFSynchronous Reference Frame
RESRenewable Energy Source
IEEEInstitute of Electrical and Electronics Engineers

Appendix A. Park Transformation

T d q = 2 3 cos ( θ ) cos ( θ 2 π 3 ) cos ( θ + 2 π 3 ) sin ( θ ) sin ( θ 2 π 3 ) sin ( θ + 2 π 3 ) 1 2 1 2 1 2

Appendix B. IEEE 9 Bus: System Parameters

Machine parameters: H 1 = 23.64 , H 2 = 6.40 , H 3 = 1.50 , x d , 1 = 0.146 , x d , 2 = 0.8958 , x d , 3 = 1.3125 , x d , 1 = 0.0608 , x d , 2 = 0.1198 , x d , 3 = 0.1813 , x q , 1 = 0.0969 , x q , 2 = 0.8645 , x q , 3 = 1.2578 , x q , 1 = 0.0969 , x q , 2 = 0.1969 , x q , 3 = 0.25 , T d o , 1 = 8.96 , T d o , 2 = 6.0 , T d o , 3 = 5.89 , T q o , 1 = 0.31 , T q o , 2 = 0.535 , T q o , 3 = 0.6
VSC parameters: S n = 100 MVA, L f = 9.737 mH, R f = 147.402 µOhm, C d c = 722 mF, V d c = 40 kV.
Synchronverter parameters: J = 0.2 , D p = 0.15 , K i , v = 0.15 , K p , p = 57.32 · 10 4 , K i , p = 859.87 · 10 4
Vector Control: K i , v = 0.64 , K p , i d = 38.2 · 10 4 , K i , i d = 0.1433 , K i , p = 1.43 , K p , d c = 3.82 , K i , d c = 795.78 · 10 4 , K i , p l l = 7.03619 · 10 6 , K p , p l l = 63.33 · 10 4
ROCOF VSG: D p = 0.1 , T j = 0.2 , H p = 5 . All other parameters similar to vector control.

Appendix C. IEEE 39 Bus: System Parameters

VSC parameters: S n = 100 MVA, L f = 9.737 mH, R f = 147.402 µOhm, C d c = 722 mF, V d c = 40 kV.
Synchronverter parameters: J = 0.25 , D p = 0.25 , K i , v = 0.15 , K p , p = 57.32 · 10 4 , K i , p = 859.87 · 10 4
Vector Control: K i , v = 0.64 , K p , i d = 38.2 · 10 4 , K i , i d = 0.1433 , K i , p = 1.43 , K p , d c = 3.82 , K i , d c = 477 · 10 4 , K i , p l l = 7.03619 · 10 6 , K p , p l l = 63.33 · 10 4
ROCOF VSG: D p = 0.15 , T j = 0.3 , H p = 5 . All other parameters similar to vector control.

Appendix D. Figures

Figure A1. IEEE 9 bus: Matrix of Participation factors for the system with syncrhonverter.
Figure A1. IEEE 9 bus: Matrix of Participation factors for the system with syncrhonverter.
Energies 14 03277 g0a1
Figure A2. IEEE 9 bus: Matrix of Participation factors for the system with vector control.
Figure A2. IEEE 9 bus: Matrix of Participation factors for the system with vector control.
Energies 14 03277 g0a2
Figure A3. IEEE 9 bus: Matrix of Participation factors for the system with ROCOF VSG.
Figure A3. IEEE 9 bus: Matrix of Participation factors for the system with ROCOF VSG.
Energies 14 03277 g0a3

References

  1. EirGrid/SONI. RoCoF Modification Proposal—TSOs’ Recommendations; EirGrid/SONI: Dublin, Ireland, 2012. [Google Scholar]
  2. Hydro-Québec. Technical Requirements for the Connection of Generation Facilities to the Hydro-Québec Transmission System; Hydro-Québec: Montréal, QC, Canada, 2003. [Google Scholar]
  3. ENTSO-E. Need for Synthetic Inertia (SI) for Frequency Regulation; ENTSO-E: Brussels, Belgium, 2017. [Google Scholar]
  4. Kerdphol, T.; Rahman, F.S.; Watanabe, M.; Mitani, Y.; Turschner, D.; Beck, H.P. Enhanced Virtual Inertia Control Based on Derivative Technique to Emulate Simultaneous Inertia and Damping Properties for Microgrid Frequency Regulation. IEEE Access. 2019, 7, 14422–14433. [Google Scholar] [CrossRef]
  5. Milano, F.; Dörfler, F.; Hug, G.; Hill, D.J.; Verbič, G. Foundations and Challenges of Low-Inertia Systems (Invited Paper) (Invited Paper). In Proceedings of the 5th 2018 Power Systems Computation Conference (PSCC), Dublin, Ireland, 11–15 June 2018; pp. 1–25. [Google Scholar] [CrossRef] [Green Version]
  6. Sakimoto, K.; Miura, Y.; Ise, T. Stabilization of a power system with a distributed generator by a Virtual Synchronous Generator function. In Proceedings of the 8th International Conference on Power Electronics-ECCE Asia, Jeju, Korea, 30 May–3 June 2011; pp. 1498–1505. [Google Scholar]
  7. Wu, H.; Ruan, X.; Yang, D.; Chen, X.; Zhao, W.; Lv, Z.; Zhong, Q. Small-Signal Modeling and Parameters Design for Virtual Synchronous Generators. IEEE Trans. Ind. Electron. 2016, 63, 4292–4303. [Google Scholar] [CrossRef]
  8. Karapanos, V.; de Haan, S.; Zwetsloot, K. Real time simulation of a power system with VSG hardware in the loop. In Proceedings of the IECON 2011—37th Annual Conference of the IEEE Industrial Electronics Society, Melbourne, VIC, Australia, 7–10 November 2011; pp. 3748–3754. [Google Scholar] [CrossRef]
  9. Saeedian, M.; Pournazarian, B.; Seyedalipour, S.S.; Eskandari, B.; Pouresmaeil, E. Emulating Rotational Inertia of Synchronous Machines by a New Control Technique in Grid-Interactive Converters. Sustainability 2020, 12, 5346. [Google Scholar] [CrossRef]
  10. Fang, J.; Lin, P.; Li, H.; Yang, Y.; Tang, Y. An Improved Virtual Inertia Control for Three-Phase Voltage Source Converters Connected to a Weak Grid. IEEE Trans. Power Electron. 2019, 34, 8660–8670. [Google Scholar] [CrossRef] [Green Version]
  11. Zhong, Q.; Weiss, G. Synchronverters: Inverters That Mimic Synchronous Generators. IEEE Trans. Ind. Electron. 2011, 58, 1259–1267. [Google Scholar] [CrossRef]
  12. Zhong, Q.; Nguyen, P.; Ma, Z.; Sheng, W. Self-Synchronized Synchronverters: Inverters without a Dedicated Synchronization Unit. IEEE Trans. Power Electron. 2014, 29, 617–630. [Google Scholar] [CrossRef]
  13. Mo, O.; D’Arco, S.; Suul, J.A. Evaluation of Virtual Synchronous Machines with Dynamic or Quasi-Stationary Machine Models. IEEE Trans. Ind. Electron. 2017, 64, 5952–5962. [Google Scholar] [CrossRef] [Green Version]
  14. Beck, H.-P.; Hesse, R. Virtual Synchronous Machine. In Proceedings of the 9th International Conference on Electrical Power Quality and Utilisation, Barcelona, Spain, 9–11 October 2007. [Google Scholar]
  15. Van Emmerik, E.L.; França, B.W.; Aredes, M. A synchronverter to damp electromechanical oscillations in the Brazilian transmission grid. In Proceedings of the 2015 IEEE 24th International Symposium on Industrial Electronics (ISIE), Buzios, Brazil, 3–5 June 2015; pp. 221–226. [Google Scholar] [CrossRef]
  16. Brown, E.; Weiss, G. Using synchronverters for power grid stabilization. In Proceedings of the 2014 IEEE 28th Convention of Electrical & Electronics Engineers in Israel (IEEEI), Eilat, Israel, 3–5 December 2014; pp. 1–5. [Google Scholar] [CrossRef]
  17. Rosso, R.; Engelken, S.; Liserre, M. Robust Stability Investigation of the Interactions Among Grid-Forming and Grid-Following Converters. IEEE J. Emerg. Sel. Top. Power Electron. 2020, 8, 991–1003. [Google Scholar] [CrossRef]
  18. Sauer, P.W.; Pai, M.A. Power System Dynamics and Stability; Prentice Hall: Hoboken, NJ, USA, 1998. [Google Scholar]
  19. IEEE Guide for Synchronous Generator Modeling Practices and Applications in Power System Stability Analyses. IEEE Std 2003, 1110–2002. [CrossRef]
  20. Zhong, Q. Four-quadrant operation of AC machines powered by inverters that mimic synchronous generators. In Proceedings of the 5th IET International Conference on Power Electronics, Machines and Drives (PEMD 2010), Brighton, UK, 19–21 April 2010; pp. 1–6. [Google Scholar] [CrossRef]
  21. Vetoshkin, L.; Müller, Z. A comparative study of synchronverter stability. In Proceedings of the 2020 21st International Scientific Conference on Electric Power Engineering (EPE), Prague, Czech Republic, 19–21 October 2020. [Google Scholar]
  22. Vetoshkin, L.; Müller, Z. A supervisory MPC for synchronverter. In Proceedings of the 2020 21st International Scientific Conference on Electric Power Engineering (EPE), Prague, Czech Republic, 19–21 October 2020. [Google Scholar]
  23. Papangelis, L.; Debry, M.; Prevost, T.; Panciatici, P.; Van Cutsem, Y. Stability of a Voltage Source Converter Subject to Decrease of Short-Circuit Capacity: A Case Study. In Proceedings of the 2018 Power Systems Computation Conference (PSCC), Dublin, Ireland, 11–15 June 2018; pp. 1–7. [Google Scholar] [CrossRef] [Green Version]
  24. Kundur, P. Power System Stability and Control; McGraw-Hill: New York, NY, USA, 1994. [Google Scholar]
  25. Anderson, P.; Fouad, A. Power System Control and Stability; IEEE Press: Piscataway, NJ, USA, 2003. [Google Scholar]
  26. Dembart, B.; Erisman, A.; Cate, E.; Epton, M.; Dommel, H. Power System Dynamic Analysis: Phase I. Final Report. 1977. Available online: https://www.osti.gov/biblio/7296152 (accessed on 2 June 2021).
  27. Wolfram Research, Inc. Mathematica; Version 12.2; Wolfram Research, Inc.: Champaign, IL, USA, 2020. [Google Scholar]
  28. Hiskens, I. IEEE PES Task Force on Benchmark Systems for Stability Controls. 2013. Available online: http://www.sel.eesc.usp.br/ieee/IEEE39/New_England_Reduced_Model_(39_bus_system)_MATLAB_study_report.pdf (accessed on 10 September 2020).
  29. Poolla, B.K.; Bolognani, S.; Dörfler, F. Optimal Placement of Virtual Inertia in Power Grids. IEEE Trans. Autom. Control 2017, 62, 6209–6220. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Time intervals of frequency response during a contingency. Reprinted from [4].
Figure 1. Time intervals of frequency response during a contingency. Reprinted from [4].
Energies 14 03277 g001
Figure 3. Synchronverter block diagram.
Figure 3. Synchronverter block diagram.
Energies 14 03277 g003
Figure 6. Conversions between reference frames.
Figure 6. Conversions between reference frames.
Energies 14 03277 g006
Figure 7. WECC three-machine, nine bus system, Reprinted from [18].
Figure 7. WECC three-machine, nine bus system, Reprinted from [18].
Energies 14 03277 g007
Figure 8. Eigenvalues of the IEEE 9 bus system with VSC connected to Node 6. Black dashed line shows ξ m i n = 0.0298 of the system with synchronverters. Red dashed line corresponds to ξ m i n = 0.0112 for ROCOF VSG control. Blue dashed line corresponds to ξ m i n = 0.0062 for vector control.
Figure 8. Eigenvalues of the IEEE 9 bus system with VSC connected to Node 6. Black dashed line shows ξ m i n = 0.0298 of the system with synchronverters. Red dashed line corresponds to ξ m i n = 0.0112 for ROCOF VSG control. Blue dashed line corresponds to ξ m i n = 0.0062 for vector control.
Energies 14 03277 g008
Figure 9. IEEE 9 bus: Load angle of 3 generator during 3 phase fault for t c l = 0.1 s.
Figure 9. IEEE 9 bus: Load angle of 3 generator during 3 phase fault for t c l = 0.1 s.
Energies 14 03277 g009
Figure 10. IEEE 9 bus: Load angle of 3 generator during 3 phase fault for t c l = 0.2 s.
Figure 10. IEEE 9 bus: Load angle of 3 generator during 3 phase fault for t c l = 0.2 s.
Energies 14 03277 g010
Figure 11. New England 10-machine, 39 bus system, Adapted from [28].
Figure 11. New England 10-machine, 39 bus system, Adapted from [28].
Energies 14 03277 g011
Figure 12. IEEE 39 bus: Damping ratios for VSCs with vector control connected to different nodes.
Figure 12. IEEE 39 bus: Damping ratios for VSCs with vector control connected to different nodes.
Energies 14 03277 g012
Figure 13. IEEE 39 bus: Damping ratios for VSCs with ROCOF VSG control connected to different nodes.
Figure 13. IEEE 39 bus: Damping ratios for VSCs with ROCOF VSG control connected to different nodes.
Energies 14 03277 g013
Figure 14. IEEE 39 bus: Damping ratios for synchronverters connected to different nodes.
Figure 14. IEEE 39 bus: Damping ratios for synchronverters connected to different nodes.
Energies 14 03277 g014
Figure 15. IEEE 39 bus: m a x ( R e ( λ ) ) for VSCs with vector control connected to different nodes.
Figure 15. IEEE 39 bus: m a x ( R e ( λ ) ) for VSCs with vector control connected to different nodes.
Energies 14 03277 g015
Figure 16. IEEE 39 bus: m a x ( R e ( λ ) ) for VSCs with ROCOF VSG control connected to different nodes.
Figure 16. IEEE 39 bus: m a x ( R e ( λ ) ) for VSCs with ROCOF VSG control connected to different nodes.
Energies 14 03277 g016
Figure 17. IEEE 39 bus: m a x ( R e ( λ ) ) for synchronverters connected to different nodes.
Figure 17. IEEE 39 bus: m a x ( R e ( λ ) ) for synchronverters connected to different nodes.
Energies 14 03277 g017
Table 1. The eigenvalues and participation factors of the system with synchronverters.
Table 1. The eigenvalues and participation factors of the system with synchronverters.
EigenvaluesVariablePF
λ 1 , 2 = 72.53 ± 1265.85 j i q , s , V d c 0.5 , 0.2
λ 3 , 4 = 170.41 ± 305.73 j δ s , i d , s , ω s 0.32 , 0.27 , 0.22
λ 5 , 6 = 0.64 ± 12.86 j ω 3 , δ 3 0.38 , 0.39
λ 7 , 8 = 5.5 ± 7.95 j E f d , 2 , V r , 2 0.4 , 0.4
λ 9 , 10 = 5.34 ± 7.93 j E f d , 3 , V r , 3 0.38 , 0.39
λ 11 , 12 = 5.25 ± 7.86 j E f d , 1 , V r , 1 0.39 , 0.41
λ 13 , 14 = 0.25 ± 8.56 j ω 2 , δ 2 0.29 , 0.37
λ 15 , 16 = 1.15 ± 6.69 j p r e g , p i , V d c 0.37 , 0.32
λ 17 = 5.76 P s v , 2 , P s v , 3 0.27 , 0.38
λ 18 , 19 = 5.26 ± 0.45 j E d , 2 , E d , 3 , P s v , 3 0.2 , 0.24 , 0.25
λ 20 = 5.17 P s v , 1 0.79
λ 21 = 3.82 E d , 2 , E d , 2 0.39 , 0.44
λ 22 = 3.23 E d , 1 1
λ 23 = 2.37 M m e c h , 1 , M m e c h , 2 0.44 , 0.45
λ 24 = 2.32 M m e c h , 3 0.56
λ 25 , 26 = 1.12 ± 0.82 j ω 1 0.23
λ 27 , 28 = 0.39 ± 0.97 j E q , 1 0.2
λ 29 , 30 = 0.44 ± 0.75 j E q , 1 , R f , 1 0.25 , 0.2
λ 31 , 32 = 0.43 ± 0.49 j E q , 3 , R f , 3 0.26 , 0.31
λ 33 = 0.09 ψ s 0.99
Table 2. The eigenvalues and participation factors of the system with vector control.
Table 2. The eigenvalues and participation factors of the system with vector control.
EigenvaluesVariablePF
λ 1 , 2 = 34.72 ± 738.89 j i d , v s c , V d c 0.49 , 0.5
λ 3 = 74.59 i q , v s c 0.92
λ 4 , 5 = 0.17 ± 27.3 j ω p l l , θ p l l 0.49 , 0.49
λ 6 , 7 = 0.61 ± 12.84 j ω 3 , δ 3 0.38 , 0.39
λ 8 , 9 = 5.5 ± 7.95 j E f d , 2 , V r , 2 0.39 , 0.4
λ 10 , 11 = 5.33 ± 7.92 j E f d , 3 , V r , 3 0.37 , 0.38
λ 12 , 13 = 5.24 ± 7.84 j E f d , 1 , V r , 1 0.38 , 0.39
λ 14 , 15 = 0.12 ± 8.48 j ω 2 , δ 2 0.31 , 0.4
λ 16 = 5.76 P s v , 2 , P s v , 3 0.26 , 0.36
λ 17 , 18 = 5.26 ± 0.45 j E d , 2 , E d , 3 , P s v , 3 0.21 , 0.23 , 0.25
λ 19 = 5.16 P s v , 1 0.76
λ 20 = 4.28 M q 0.53
λ 21 = 3.23 E d , 1 1
λ 22 = 3.04 E d , 2 , E d , 3 , M q 0.25 , 0.28 , 0.3
λ 23 = 2.37 M m e c h , 1 , M m e c h , 2 0.41 , 0.46
λ 24 = 2.31 M m e c h , 3 0.59
λ 25 , 26 = 1.03 ± 0.95 j ω 1 0.29
λ 27 , 28 = 0.42 ± 1.15 j E q , 1 0.2
λ 29 , 30 = 0.15 ± 1.17 j i d , r e f , M d 0.49 , 0.48
λ 31 , 32 = 0.44 ± 0.75 j E q , 1 , R f , 1 0.21 , 0.28
λ 33 , 34 = 0.43 ± 0.49 j E q , 3 , R f , 3 0.24 , 0.29
λ 35 = 0.07 i q , r e f 0.97
λ 36 = 0.02 p r e g , p i 1
λ 37 = 0.02 I p l l 1
Table 3. The eigenvalues and participation factors of the system with ROCOF VSG.
Table 3. The eigenvalues and participation factors of the system with ROCOF VSG.
EigenvaluesVariablePF
λ 1 , 2 = 8.27 ± 739.62 j i d , v s c , V d c 0.49 , 0.5
λ 3 , 4 = 1.25 ± 26.25 j ω p l l , θ p l l 0.44 , 0.45
λ 5 , 6 = 11.41 ± 13.03 j i q , v s c , M q 0.44 , 0.4
λ 7 , 8 = 0.67 ± 12.91 j ω 3 , δ 3 0.39 , 0.4
λ 9 , 10 = 5.49 ± 7.95 j E f d , 2 , V r , 2 0.38 , 0.38
λ 11 , 12 = 5.32 ± 7.91 j E f d , 3 , V r , 3 0.32 , 0.33
λ 13 , 14 = 5.24 ± 7.82 j E f d , 1 , V r , 1 0.32 , 0.33
λ 15 , 16 = 0.31 ± 8.56 j ω 2 , δ 2 0.31 , 0.39
λ 17 = 5.85 j P s v , 2 , P s v , 3 0.3 , 0.32
λ 18 , 19 = 5.26 ± 0.44 j E d , 2 , E d , 3 , P s v , 3 0.21 , 0.22 , 0.23
λ 20 = 5.17 P s v , 1 0.81
λ 21 = 4.99 V J , r e f 0.98
λ 22 = 3.29 E d , 2 , E d , 3 0.37 , 0.39
λ 23 = 3.23 E d , 1 1
λ 24 = 2.37 M m e c h , 1 , M m e c h , 2 0.4 , 0.48
λ 25 = 2.31 M m e c h , 3 0.6
λ 26 , 27 = 1.01 ± 0.96 j ω 1 0.28
λ 28 , 29 = 0.44 ± 1.17 j E q , 1 0.2
λ 30 , 31 = 0.02 ± 1.17 j i d , r e f , M d 0.45 , 0.45
λ 32 , 33 = 0.44 ± 0.75 j E q , 1 , R f , 1 0.22 , 0.28
λ 34 , 35 = 0.43 ± 0.49 j E q , 3 , R f , 3 0.24 , 0.29
λ 36 = 0.07 i q , r e f 0.97
λ 37 = 0.02 p r e g , p i 1
λ 38 = 0.02 I p l l 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vetoshkin, L.; Müller, Z. A Comparative Analysis of a Power System Stability with Virtual Inertia. Energies 2021, 14, 3277. https://doi.org/10.3390/en14113277

AMA Style

Vetoshkin L, Müller Z. A Comparative Analysis of a Power System Stability with Virtual Inertia. Energies. 2021; 14(11):3277. https://doi.org/10.3390/en14113277

Chicago/Turabian Style

Vetoshkin, Lavr, and Zdeněk Müller. 2021. "A Comparative Analysis of a Power System Stability with Virtual Inertia" Energies 14, no. 11: 3277. https://doi.org/10.3390/en14113277

APA Style

Vetoshkin, L., & Müller, Z. (2021). A Comparative Analysis of a Power System Stability with Virtual Inertia. Energies, 14(11), 3277. https://doi.org/10.3390/en14113277

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