1. Introduction
The transition to power systems involving more power electronics-based converters represents a challenge today. An increase in the penetration of inverter-interfaced renewable energy sources (such as wind and sun) may lead to instabilities within the power system [
1,
2]. For this reason, accurate converter models, taking into consideration both the converter hardware parameters and its control algorithm, will be necessary. Moreover, some grid codes are now requiring dynamic models on different degree of details for simulation and analysis at the power system level [
3,
4].
Several techniques are presently available for this purpose. The most popular are frequency domain analysis using impedance models and the eigenvalue analysis using state-space models [
5]. Regarding the impedance models, various approaches are available in the literature [
6,
7,
8,
9,
10,
11], as well as experimental procedures [
12] to obtain the black-box frequency response of a converter without knowing its internal parameters. Regarding the system-level analysis, the Generalized Nyquist Criterion can be used [
13] if the individual equivalent impedance models are available. However, impedance models fail to provide an immediate comprehension of the poles in the system and how they are influenced by the system parameters. Moreover, by observing the Bode plot of the equivalent impedance of a converter, the effect of some poles can be canceled by corresponding zeros of the system.
On the contrary, the state-space approach has the advantage of clearly identifying the modes of the system under study and their correlation both to its states (e.g., controller, filter) using participation factors [
14] and to the system parameters using sensitivity analysis. Moreover, the state-space model can provide an immediate feedback of the different poles time constants, without incurring the effects of zero-pole cancellation that can be present in the impedance model representation.
Traditionally, compared to synchronous generators, grid-connected inverters have been working with much smaller time constants. These time constants are related to both the physical hardware of the converters (i.e., grid side filters) and the digital controllers. However, to obtain a higher integration of renewable energy sources, inverters may be required to provide ancillary services to the grid, in order to maintain its stability and quality (frequency, voltage and harmonic content). These services typically operate with much slower time constants, in the order of magnitude of the line frequency or less. For these reasons, a multi time scale analysis of the system may be necessary. Also, simplified models may be desirable, focusing only on certain phenomena, to reduce the computational burden of otherwise extensive simulations to assess the stability and the behavior of the system.
However, state-space modeling is a very difficult task when dealing with complex systems internally interconnected. An increase in the number of state variables, inputs and outputs of the system leads to complicated matrix expressions. Moreover, another disadvantage is its low flexibility: if components are added, removed or modified, the modeling process must be repeated from the beginning. This disadvantage is particularly evident when modeling converters equipped with advanced cascaded or more emerging controllers. Such controllers may be implemented to provide additional features beyond the traditional current or voltage control of grid-connected converters and are an enabling technology to integrate more renewable energy sources in the future power systems.
An example of outer advanced controller are the virtual synchronous generators (VSGs), which can provide virtual inertia, harmonic compensation and reactive support during voltage dips [
15,
16,
17,
18,
19,
20,
21,
22].
The literature proposes some examples of state-space modeling of such controllers [
18], but following a monolithic approach, i.e., deriving the complete model at once. This kind of approach is feasible, but it lacks of flexibility, modularity and scalability (i.e., if something changes, the model must be derived again). Moreover, the number of state variables involved dramatically increases the order of the obtained system. This complexity leads to a larger modeling effort and it is prone to errors. Therefore, troubleshooting cannot be easily done during the derivation of the model.
A valid methodology to tackle such issues is the Component Connection Method (CCM). This method solves the lack of modularity and simplifies the modeling procedure, reducing both the analytical effort and the probability of making mistakes during the derivation, by separately modeling the single components of the system and then linking them using sparse interconnection matrices, as shown in
Figure 1.
To provide a global and practical perspective on this technique, this paper reviews the requirements of power electronics modeling in power systems and the most recent applications of CCM to this matter. Moreover, it contributes with a generalized step-by-step Component Connection Method (CCM)-based modeling approach of the individual converter unit, considering both the hardware and the control part. Compared to the classical CCM theory and existing applications, purely algebraic blocks are introduced in this paper and it is described how to include them in the CCM in a completely general way. Such blocks have no physical meaning in power systems applications, but they simplify the modeling of power electronics controllers (e.g., calculating the current references starting from the power setpoints and the grid voltage).
This paper is organized as follows. In
Section 2 the state of the art of CCM applied to power electronics-based power systems is reviewed, along with its most recent applications.
Section 3 presents a step-by step procedure of CCM-based modeling. A practical particular example of such method by modeling the Simplified Virtual Synchronous Compensator (S-VSC) [
23], is given in
Section 4.
Section 5 presents some analyses, which can be performed thanks to the state-space model and the experimental validation of the presented case study. Finally, the conclusions are provided in
Section 6.
2. State-Space Modeling Literature
With the spread of distributed generation and microgrids based on power electronics converters, accurate and efficient modeling techniques have become necessary. Such methods must guarantee flexibility (i.e., being easily applicable to different systems), scalability (i.e., can model more than one converter and can merge multiple grid sub-portions) and various detail level, depending on the specific analysis needed.
For this purpose, various levels of complexity and detail have been proposed during the years. At first, the focus was on low frequency subsynchronous oscillations, neglecting the faster behavior of current/voltage controllers [
24]. In such cases, the inverter was considered as an ideal voltage source, controlled in amplitude and frequency by active and reactive droop controllers. However, such level of simplification is not always accurate in predicting instabilities or poorly damped resonances occurring at higher frequencies (innermost current or voltage control bandwidth and grid-side filters resonance frequency). Therefore, more accurate models have been proposed [
25,
26].
Given these requirements and summarizing the benefits of state-space modeling briefly mentioned in
Section 1:
The modes of the system under study are evident in their frequency and damping. Therefore, poorly damped oscillations can be easily identified;
The interaction among the modes and the states can be analyzed by means of participation factors. The states associated with the most critical modes are then identified;
A reduced order model of the system can be obtained, by neglecting the dynamic behavior of some of its parts. This allows faster simulations without a decrease in the quality of the results.
CCM is a valid modeling approach. In fact, the main advantage of CCM is the possibility of decomposing an articulated system into its fundamental blocks, which can then be individually modeled, and the possibility of reconfiguring the system by simply adding new blocks or just modifying the interconnection matrices. Such advantage is especially beneficial when modeling a single converter unit. This allows an easy an quick derivation of the state-space model of the converter, even when a single or more control layers (e.g., current controller) or hardware component (e.g., grid side filter) are modified in a later stage of the design.
CCM has been traditionally used at power system level to analyze multi-node networks with several conventional generating units or loads starting from [
27]. Later, it has been specifically used to model the early grid applications of power electronics: static var compensators and HVDC transmission systems [
28,
29]. This method has gained renewed interest in the last decade, when it was applied to modern power electronics-based power systems. First, in 2014 CCM was applied to the specific modeling of wind farms and their connection to the main grid [
8,
30,
31]. The same research evolved then in CCM modeling applied to multi converter power system, being able to include highly detailed inverter models (i.e., including accurate control loops and digital delays models) in wider inverter-based power systems.
CCM is currently used to analyze the stability and the interaction among conventional (i.e., phase locked loop (PLL)-based) and new controllers (i.e., Synchronverters) for renewable generation [
32,
33,
34]. In these works, the CCM is the first stage of a more advanced state-space
-analysis [
34] to consider the uncertainties of the modeled plants.
CCM is still also currently adopted in power system studies, to include inverter-based grid nodes into larger scale studies (e.g., Subsynchronous torsional interactions [
35]). In such cases, its flexibility is exploited in order to simplify the modeling procedure neglecting aspects at frequency out of the range of interest (e.g., digital control delay). The history of CCM is schematically summarized in
Figure 2.
3. Generalized CCM Modeling Approach
This section considers a step-by-step procedure to model a grid-connected three-phase two level inverter, interfaced to the grid by means of an LCL filter and is depicted in
Figure 3. The DC supply of the inverter is considered to be an ideal voltage source. Thanks to the modular approach guaranteed by the CCM, a more accurate model of the DC source, or even the model of a more complex DC side, e.g., a DC microgrid [
36,
37], can be easily included. The grid is modeled as a Thévenin equivalent circuit.
The CCM consists of the modular analysis of the system under modeling. First, the single components are identified and linearized around the chosen operating point. Then, the connections between the components are defined by means of sparse interlinking matrices. Finally, the state-space model of the overall system is obtained by linear algebra calculation.
This process is performed by the following general step-by-step procedure described below:
Define the system inputs and outputs . Examples of inputs are external disturbances such as references variations (e.g., active and reactive power) and grid faults (i.e., frequency, voltage amplitude or phase variations). The outputs of the system can be chosen arbitrarily, depending on the aspect under study. Examples are the power (or current) injected by the inverter and the measured grid frequency using PLLs or VSGs.
Identify the single components. Two kinds of components can be used: dynamic (with state variables) and algebraic (without state variables). In general, n dynamic components and m algebraic components are identified and linearized.
Each
k-th linearized dynamic component can be written in the following state-space form:
where
represents the state variables vector of the block,
the inputs of the component and
its outputs. The matrices
,
,
and
are the component state-space matrices. Examples of the dynamic components are LCL filters and integral regulators.
Algebraic components are defined to simplify the modeling procedure. In a traditional state-space representation, only dynamic components are defined. However, this distinction can be very useful when it comes to the modeling of a controller. Many algebraic calculations, such as divisions, lead to complicated linearized expressions and it can be convenient to simplify the process by defining such algebraic components. They do not have state variables, but only inputs
and outputs
related as:
Droop controllers, sums and divisions to calculate the current reference from the active and reactive power references are examples of algebraic blocks.
The vector of the state variables of the global system
is the aggregate of the
state variable vectors of the single components. The list of the components inputs
U and the list of the components outputs
Y can be defined in an equivalent way as follows:
An aggregated model is defined by composing the individual system matrices as follows:
where:
In (5), it must be noted that the and matrices must be extended by adding null elements to compensate for the algebraic blocks. In particular, a number of null columns equal to the total number of inputs of all the algebraic blocks must be attached to the matrix. Also several null rows equal to the sum of the algebraic block outputs must be attached to the matrix.
Define the connection matrices, to connect the single components.
These matrices
T are sparse and connect the inputs and outputs of the system as follows:
The global linearized system state-space model can be obtained as:
where:
when the system structure is modified, U and Y must be modified according to the new blocks, as well as the connection matrices.
As mentioned in the introduction, the CCM is useful when performing a multi time scale analysis of a system. The reason behind this kind of analysis is that the most recent controllers, enabling power electronic converters provide ancillary services to the grid (e.g., VSGs) typically work on different time scales. The inner control loops (voltage, current) are operating at the time scale of the switching frequency. On the other hand, the outer level controllers, managing for example the virtual inertia, operate with time constants in the order of magnitude of seconds. To perform a simplified analysis with a traditional state-space modeling approach, the system has to be completely remodeled, meaning time loss and possibility of introducing errors. The CCM solves this problem: the user can easily exclude some levels of the system by substituting them with purely algebraic blocks. This way the computational burden for simulations is reduced, while preserving the correctness of the analysis at the time scale of interest. This simplification affects only the state variables of the
k-th excluded block, but does not affect its inputs and outputs. Therefore, a dynamic block can be transformed into an algebraic block and the global system is easily derived again, without any need to change the connection matrices. The poles relative to the neglected blocks will not be present in the global system anymore. The pole map of the system and a participation factors analysis, which will be described in
Section 5, are a good starting point to decide which blocks can be neglected to perform which study, since they give a clear and straightforward view of the time constant of each pole and which blocks are related to them.
Moreover, the CCM is flexible and solves also the issues related to the comparison or improvements of the inverter control and hardware. In fact, CCM decompose an articulated system into its fundamental blocks, which can then be individually modeled, and allows the reconfiguration of the system by simply adding new blocks or just modifying the interconnection matrices. Such advantage is especially beneficial when modeling a single converter unit. The benefits are especially evident in comparison to the traditional monolithic modeling approach, which has to be repeated again when any part of the system is modified. This flexibility of CCM allows an easy and quick derivation of the new state-space model of the converter, even when a single or more control layers (e.g., current controller) or hardware component (e.g., grid side filter) are modified in a later stage of the design.
A practical example of the CCM superiority can be done by comparing it with the traditional state-space modeling. The literature [
18] gives examples of VSG modeling with traditional techniques. In this example the complete state-space model of a grid-connected converter controlled with a cascaded VSG is described. However, the equivalent model is obtained directly, with a monolithic approach. The result is a 19 × 19 state variable matrix, with no easy identification of which part of this matrix corresponds to which block of the converter. With the CCM, as demonstrated by the approach of (5), this is not the case. Each functional block of the system under study corresponds to a well defined matrix, which can be modified in a later stage. For example, in case a different current control strategy is implemented in a later stage of the design, its inputs and outputs will not probably change, but only its internal structure. Several controllers can be therefore compared, as well as different approximations of the delay of the digital controller.
4. Modeling of the S-VSC
In this section, an example of the application of the CCM is given. The system under study is the Simplified Virtual Synchronous Compensator (S-VSC). This model has been chosen as a representative of a VSG-based controller providing ancillary services to the grid. However, CCM can be applied to model any other controller available in the literature, by being a general modeling approach.
The system consists of the inverter hardware part, its control (current controller and S-VSC) and the connection to the equivalent grid, as shown in
Figure 4. The system is modeled in per unit with
base voltage,
base power and
base angular speed.
More details on this model are available in the literature [
23].
4.1. Defining System Inputs and Outputs
First, the inputs and the outputs of the system are defined. The system inputs are defined to model the possible electrical perturbations that can influence the system. In particular, step variations in the active and reactive power references (, ), grid frequency , grid voltage amplitude and grid voltage phase angle variations. This choice allows full testing of the system under not-rated operating conditions.
The system outputs
are arbitrarily selected, depending on the quantities of interest. In this example, they have been chosen as follows:
where
and
are the active and reactive powers injected by the inverter.
4.2. Component Identification
The components are identified as shown in
Figure 4.
dynamic blocks and
algebraic blocks are defined as follows:
Dynamic components: LCL filter, Inverter control loops (PI and delay model), S-VSC electrical part (stator and damper), S-VSC power loops (mechanical part and excitation control);
Algebraic components: Power reference calculation (power to current) and grid perturbation model.
4.3. Dynamic Components Definition
4.3.1. LCL Filter
The inverter is connected to the grid by means of an LCL filter, as depicted in
Figure 5. The filter is modeled in the (
d,
q) frame synchronous to the S-VSC virtual rotor position
, rotating at
.
The equations modeling this component (in per unit) are as follows:
This block has the following state variables
, inputs
and outputs
:
The output can be useful in the case of a multi-inverter plant where the voltage amplitude control at the plant point of common coupling (PCC) is required.
By linearizing (
10), the necessary state-space matrices
are obtained.
4.3.2. Current Controller
The current controller considered in this model has a sampling time of
and it is based on a proportional integral (PI) regulator implemented in the (
d,
q) synchronous reference. The proportional and integral gains of the PI regulator are
and
. The total delay
of the digital controller and the modulation is modeled with a first order Padé approximation. A higher order approximation can easily be considered by modifying this block. The controller complete block diagram is shown in
Figure 6.
The following state variables
are identified (the integrator
and the delay model
), as well as the inputs
and outputs
of this block:
The state-space matrices of this block are .
4.3.3. S-VSC Electromagnetic Equations
The electrical and magnetic equations of both virtual stator and rotor of the S-SVC [
23] are (in per unit):
where
,
,
and
are the virtual flux linkages of the machine (
d-axis,
q-axis, damper winding and excitation, respectively);
and
are the machine virtual currents;
and
are the virtual stator resistance and inductance;
and
are the virtual damper parameters, tuned as described in [
38].
Therefore, the block state variables
, inputs
and outputs
are defined as follows:
The state-space matrices
of this block are obtained from the linearization of (
13).
4.3.4. S-VSC Power Loops Equations
The active and reactive power control equations of the S-SVC [
23] are (in per unit):
where
and
are the virtual active and reactive power of the S-VSC. The S-VSC is always operated at zero reference power in order to obtain better stability and damping [
23].
H is the inertia constant of the machine,
is the virtual load angle of the machine (
q-axis to grid voltage vector).
is the peak value of the voltage across the filter capacitor
and
is the gain of the excitation control.
This block state variables
, inputs
and outputs
are as follows:
The state-space matrices
of the power loops block are again obtained by the linearization of (
15).
4.4. Algebraic Components Definition
4.4.1. Power Reference Calculation
In this algebraic block the following calculation is performed to obtain the current references
from the power references
and the grid voltage
:
where:
Are the sums of the external references and the S-VSC compensation power references.
Due to the division involved in (
17), a dedicated block is justified to simplify the modeling process.
There are no state variables, while the inputs
and outputs
of the block are:
Being an algebraic block, only the matrix is obtained.
4.4.2. Grid Perturbations
This algebraic block generates the grid voltage variations according to the external disturbances , given as system inputs.
The grid voltage vector is defined as:
Any phase displacement
in the grid results in the variation of the load angle
. As seen in the previous S-VSC power loops section, the grid frequency
variation leads to a load angle variation and it is, therefore, included in (
15). The vector diagram of the grid and S-VSC is shown in
Figure 7, as well as the angle variation due to a phase jump in the grid.
This block is defined by the following inputs
and outputs
:
And from the matrix
, obtained by linearizing (
20).
4.5. Aggregated Model
Once all the dynamic and algebraic blocks have been defined, the aggregated model can be obtained. First the aggregated state variables
, inputs
U and outputs
Y are defined:
Then, the aggregated model matrices are:
Note that the matrices
and
are extended with several null elements equal to the number of inputs (8 in total) and outputs (4 in total) of the two algebraic blocks. The detailed expressions of the matrices in (
23) are available in the
Appendix A.
4.6. Connection Matrices
Now the connection matrices can be obtained by linking the inputs
U and the outputs
Y of the aggregated system together with the global inputs
and outputs
. The matrix
is zero, since there is no direct feed-through from the system inputs to the system outputs.The detailed expressions of the connection matrices are available in the
Appendix A.
4.7. Derivation of the Global State Space Model
Finally, by applying (
8), the system matrices
,
,
,
are obtained. It must be noted that such model could have also been derived in a monolithic way, by studying the converter as a whole. However, by applying CCM, the model can be handled more easily and allows a quick and straightforward modification, in case some parts of the controller or of the physical system are changed. A practical and straightforward example of this advantage is the modification of the current feedback of the current controller. Usually, for grid-connected converters interfaced with LCL filters to the grid, two types of current feedbacks are adopted: grid side current feedback and converter side current feedback [
39]. In the proposed example, the converter side current feedback was adopted. Therefore, the outputs of the LCL filter block
features the converter side current. However, if the current feedback is shifted to the grid side current, only the
vector and the relative matrices
and
must be modified. Each other block is not altered by this modification. Also the interconnection matrices are kept constant, being the two feedbacks equivalent from a signal routing point of view. A second example of the reconfiguration by adding new blocks is represented by the addition of an external droop controller. If a primary frequency regulation is required to the plant, then an extra proportional frequency controller is added, in order to generate the necessary active power references for the control. With the adopted CCM, the existing model is not modified, but an extra algebraic block is added before the reference calculation block. This extra block will receive as inputs the nominal and the actual grid frequency and generate the active droop power reference, which is then summed to the external power references.
5. Analysis and Validation
5.1. Analysis of the Derived State-Space Model
The poles
of the system can be obtained numerically by solving:
In
Figure 8, the poles of the system (also listed in
Table 1) are depicted in the complex plane. The parameters of the state-space model have been chosen according to the experimental setup and are listed in
Table 2. The classical observations about the stability and damping of the system can be done. Moreover, the poles of the system can visually be distinguished as follows:
The high frequency poles (1–4) are related to the LCL filter (in the range of 2 kHz);
the equivalent poles of the digital control (5–8, 11, 12);
the poles (9, 10) of the virtual stator of the S-VSC and the grid (50 Hz);
the low frequency poles (13–16) of the electromechanical part of the S-VSC (excitation control, damper winding and swing equation) in the range of 1 Hz.
To clearly define which phenomena are related to the single poles, the analysis of the participation factors [
14] can be performed. The participation factor
of the
k-th state to the
i-th mode is defined as:
where
is the
k-th value of the right eigenvector
and
is the
k-th value of the left eigenvector
.
From (
26), the states and the modes of the system can be correlated as the results in
Table 1. In this way, the poles can be associated with the physical quantities that influence them.
From the point of view of a multi-timescale analysis, the poles whose time constants are too small or too large can be excluded to simplify the model. For example, in the considered case study, the model can be simplified to analyze the low frequency behavior of the system (e.g., integration in a wider network to study the inertial behavior of the system during a frequency drop), by neglecting the high frequency poles related to the LCL filter physical parameters and from the current controller, which can be approximated with a direct feed-through with no dynamics, being orders of magnitude faster than the electromechanical dynamics of the S-VSC. The blocks and their interconnections can be rearranged to obtain a simplified model. This simplified model is useful to study the low frequency interaction within a more complex network.
Moreover, the state-space model allows the study of the behavior of the system under a change of parameters. A practical example is the variation of the short circuit ratio (SCR) in the connection to the grid, representing the inverse of the magnitude of the grid impedance in per unit. This analysis is especially important for wind power applications connected to very weak grids. In fact, due to the often isolated geographical positions, long cables or overhead lines connections are necessary and SCR ranging down to 1.5 or even less have been reported [
40].
In
Figure 9 the SCR of the system is swept from 10 to 1.5. The S-SVC damper and excitation control parameters are tuned accordingly. As it can be seen from the pole map, the poles relative to the LCL filter (1–4) are altered in their frequency, as the grid side inductance is modified. The S-VSC stator poles (9, 10) are also affected, being the stator in an equivalent series connection with the grid inductance. The effect is a lower damping of such poles that can be however easily compensated for by increasing the virtual stator resistance
. Finally, the low frequency mechanical poles (13–15) do not lose damping, thanks to the tuning algorithm described in [
38], but only change their natural frequency. The excitation control pole (16) is not modified, as it is tuned to always obtain the desired time constant.
As previously mentioned, the state-space model can be used for a wide range of analyses. Analyses of interest are, for example, the pole map for different operating points (i.e., operation under load) and the pole variation due to the influence of parametric uncertainty of the physical components (i.e., grid inductance, filter parameters) on the tuning of both the controllers and the VSG.
5.2. Experimental Validation
The obtained state-space model can then be verified by a comparison with either a simulation of the complete system (e.g., Simulink, PLECS, PSCAD and other simulation tools can be used) or experimentally.
In this case, the S-VSC has been implemented on a dSPACE DS 1007, controlling a 15 kVA three-phase inverter connected to a grid emulator, as shown in
Figure 10. The control runs at
kHz, which also corresponds to the switching frequency
of the converter.
The parameters of both the experimental setup and the S-VSC are listed in
Table 2.
The following four tests have been performed to validate the state-space model using all the five available inputs , , , and :
Step active power reference variation ;
Step reactive power reference variation ;
Step drop of Grid frequency ;
Grid voltage drop with phase jump .
In the first test (
Figure 11) the active power fed to the grid and the frequency of the S-VSC are compared with the outputs of the state-space model. The active power reference rises from 0.2 pu to 0.3 pu at
s. This choice allows a more general test of the state-space model. In fact, as it can be seen in (
A5), a large part of the matrix
depends on the initial active and reactive power references
and
. Therefore, a non zero initial power operating point has been chosen, better to verify this part of the model. As it can be seen from
Figure 11 both the S-VSC virtual speed
and the injected active power are similar to the results from the state-space model. The state space-model, as it has been described, does not take into account the non-ideal behavior of the measurement process (e.g., noise, delays), which on the other hand may affect the experimental setup.
In the second test (
Figure 12), the reactive power injected by the inverter is compared, when the reactive power reference rises from 0.1 pu to 0.2 pu. As with Test 1, the initial non-zero operating point has been chosen better to test the model. Also in this case, the obtained model is accurate enough to simulate the step transient in the reactive power reference. Again, the state-space model is not including the non-ideal behavior of the experimental setup. This non-ideal behavior can be neglected when dealing with longer timescales, such as in Test 2. In fact, the delays and the noise act with much shorter time constants, which do not affect the behavior during these transients (operating in the time scale of seconds).
Test 3 (
Figure 13) deals with the input relative to the grid frequency variation
. A step frequency variation of −0.2 Hz has been applied by the grid emulator. The S-VSC tries to compensate the frequency drop by injecting active power into the grid, providing a virtual inertial effect. The state-space model well predicts both the S-VSC speed
profile (
Figure 13a) and the amount of active power injected (
Figure 13b). The modeling of the frequency variations, as obtained here, can be useful in system level studies, when analyzing frequency regulation and inertial support. This kind of analysis also helps studying the inertial effect resulting from a different tuning of the S-VSC (the key parameters are in this case the inertia constant
H and the virtual damper parameters
and
, which are directly related to the electromechanical damping
of the virtual machine [
38]).
Finally, Test 4 (
Figure 14) validates the behavior of the system in case of a realistic grid fault. During the fault, the voltage of the grid drops to 95% of its nominal value and a phase jump of
is applied. For such tests, the quantities of interest are the reactive power injected into the grid (
Figure 14a) to provide reactive support and the virtual excitation flux of the S-VSC
(
Figure 14b) to evaluate the time constant of the excitation control. As expected, the state-space model is well predicting both the shape and the amplitude of such quantities, proving itself as a good analysis tool to correctly tune the S-VSC to guarantee the desired reactive response during grid faults.
6. Conclusions
In this paper, a general step-by-step procedure to easily obtain the state-space model of a grid-connected converter equipped with advanced cascaded controllers using CCM is presented. This mathematical method is suitable to model cascaded controllers with interconnected signals thanks to its modular approach. Moreover, it is easily possible to modify or reconfigure the system, without remodeling it.
This approach can be used to derive state-space models of grid-tied inverters equipped with advanced cascaded controllers or emerging control techniques, such as VSGs in order to study and simulate the stability of a power system with a higher integration of renewable energy sources. In addition, state-space models are favorable in case of multi-timescale analysis of the system: the eigenvalues of the system can be obtained, as well as their correlation to the system state variables. Therefore, simplified models can be obtained, neglecting high or low frequency phenomena which are out of the scope of the analysis.
The single models obtained following this procedure can represent a useful tool for both power electronics designers and grid operators. The first can obtain preliminary information on the weak spots of the system they are designing and test their behavior under different operating conditions. From the grid side, such models would be provided by the manufacturers with different level of details [
3,
4] to be integrated in a more system-level analysis.
This step-by-step approach has been verified with a practical example of the modeling of a VSG model (S-VSC), which is then validated experimentally on a 15 kVA grid-connected inverter.
Author Contributions
Conceptualization, F.M., E.C., R.B., T.D. and F.B.; methodology, F.M. and E.C.; software, F.M.; validation, F.M.; formal analysis, F.M. and E.C.; investigation, F.M. and E.C.; resources, S.M., R.B. and F.B.; data curation, F.M.; writing–original draft preparation, F.M.; writing–review and editing, F.M., R.B., T.D. and F.B.; visualization, F.M., S.M., R.B., T.D. and F.B.; supervision, R.B., T.D. and F.B.; project administration, S.M., R.B., T.D. and F.B. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Conflicts of Interest
The authors declare no conflict of interest.
Abbreviations
The following abbreviations are used in this manuscript:
CCM | Component Connection Method |
VSG | Virtual Synchronous Generator |
HVDC | High Voltage Direct Current |
S-VSC | Simplified Virtual Synchronous Compensator |
PLL | Phase Locked Loop |
PI | Proportional Integral |
Appendix A
The subscript 0 indicates the value around which the linearization is performed.
Current controller matrices:
S-VSC electromagnetic part matrices:
S-VSC power loops matrices:
Power to current reference calculation:
Grid perturbation matrix:
Connection matrices. Only the non-zero elements are given:
References
- Li, C. Unstable Operation of Photovoltaic Inverter from Field Experiences. IEEE Trans. Power Deliv. 2018, 33, 1013–1015. [Google Scholar] [CrossRef]
- Mollerstedt, E.; Bernhardsson, B. Out of control because of harmonics-an analysis of the harmonic response of an inverter locomotive. IEEE Control Syst. Mag. 2000, 20, 70–81. [Google Scholar] [CrossRef]
- Terna Grid Code. Allegato A17—CENTRALI EOLICHE—Condizioni Generali di Connessione Alle Reti AT Sistemi di Protezione Regolazione e Controllo; Terna Grid Code Attachment A17; TERNA: Roma, Italy, 2018. [Google Scholar]
- Terna Grid Code. Allegato A68—CENTRALI FOTOVOLTAICHE—Condizioni Generali di Connessione Alle Reti AT Sistemi di Protezione Regolazione e Controllo; Terna Grid Code—Attachment A68; TERNA: Roma, Italy, 2018. [Google Scholar]
- Wang, X.; Blaabjerg, F. Harmonic Stability in Power Electronic-Based Power Systems: Concept, Modeling, and Analysis. IEEE Trans. Smart Grid 2019, 10, 2858–2870. [Google Scholar] [CrossRef] [Green Version]
- Sun, J. Small-Signal Methods for AC Distributed Power Systems—A Review. IEEE Trans. Power Electron. 2009, 24, 2545–2554. [Google Scholar] [CrossRef]
- Sun, J. Impedance-Based Stability Criterion for Grid-Connected Inverters. IEEE Trans. Power Electron. 2011, 26, 3075–3078. [Google Scholar] [CrossRef]
- Wang, X.; Blaabjerg, F.; Wu, W. Modeling and Analysis of Harmonic Stability in an AC Power-Electronics-Based Power System. IEEE Trans. Power Electron. 2014, 29, 6421–6432. [Google Scholar] [CrossRef]
- Wang, X.; Harnefors, L.; Blaabjerg, F. Unified Impedance Model of Grid-Connected Voltage-Source Converters. IEEE Trans. Power Electron. 2018, 33, 1775–1787. [Google Scholar] [CrossRef] [Green Version]
- Wen, B.; Boroyevich, D.; Burgos, R.; Mattavelli, P.; Shen, Z. Analysis of D-Q Small-Signal Impedance of Grid-Tied Inverters. IEEE Trans. Power Electron. 2016, 31, 675–687. [Google Scholar] [CrossRef]
- Liu, H.; Xie, X.; Liu, W. An Oscillatory Stability Criterion Based on the Unified dq-Frame Impedance Network Model for Power Systems With High-Penetration Renewables. IEEE Trans. Power Syst. 2018, 33, 3472–3485. [Google Scholar] [CrossRef]
- Francis, G.; Burgos, R.; Boroyevich, D.; Wang, F.; Karimi, K. An algorithm and implementation system for measuring impedance in the D-Q domain. In Proceedings of the 2011 IEEE Energy Conversion Congress and Exposition, Phoenix, AZ, USA, 17–22 September 2011; pp. 3221–3228. [Google Scholar] [CrossRef] [Green Version]
- Belkhayat, M. Stability Criteria for AC Power Systems with Regulated Loads. Ph.D. Thesis, Purdue University, West Lafayette, IN, USA, 1997. [Google Scholar]
- Kundur, P. Power System Stability and Control; McGraw-Hill Education: New York, NY, USA, 1994. [Google Scholar]
- Fang, J.; Li, H.; Tang, Y.; Blaabjerg, F. On the Inertia of Future More-Electronics Power Systems. IEEE J. Emerg. Sel. Top. Power Electron. 2019, 7, 2130–2146. [Google Scholar] [CrossRef]
- Beck, H.; 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]
- Zhong, Q.; Weiss, G. Synchronverters: Inverters That Mimic Synchronous Generators. IEEE Trans. Ind. Electron. 2011, 58, 1259–1267. [Google Scholar] [CrossRef]
- D’Arco, S.; Suul, J.A.; Fosso, O.B. A Virtual Synchronous Machine implementation for distributed control of power converters in SmartGrids. Electr. Power Syst. Res. 2015, 122, 180–197. [Google Scholar] [CrossRef]
- Driesen, J.; Visscher, K. Virtual synchronous generators. In Proceedings of the 2008 IEEE Power and Energy Society General Meeting—Conversion and Delivery of Electrical Energy in the 21st Century, Pittsburgh, PA, USA, 20–24 July 2008; pp. 1–3. [Google Scholar] [CrossRef]
- Van Wesenbeeck, M.P.N.; De Haan, S.W.H.; Varela, P.; Visscher, K. Grid tied converter with virtual kinetic storage. In Proceedings of the 2009 IEEE Bucharest PowerTech, Bucharest, Romania, 28 June–2 July 2009; pp. 1–7. [Google Scholar] [CrossRef]
- Bevrani, H.; Ise, T.; Miura, Y. Virtual synchronous generators: A survey and new perspectives. Int. J. Electr. Power Energy Syst. 2014, 54, 244–254. [Google Scholar] [CrossRef]
- Rodríguez, P.; Citro, C.; Candela, J.I.; Rocabert, J.; Luna, A. Flexible Grid Connection and Islanding of SPC-Based PV Power Converters. IEEE Trans. Ind. Appl. 2018, 54, 2690–2702. [Google Scholar] [CrossRef]
- Mandrile, F.; Carpaneto, E.; Bojoi, R. Grid-Tied Inverter with Simplified Virtual Synchronous Compensator for Grid Services and Grid Support. In Proceedings of the 2019 IEEE Energy Conversion Congress and Exposition (ECCE), Baltimore, MD, USA, 29 September–3 October 2019; IEEE: Baltimore, MD, USA, 2019; pp. 4317–4323. [Google Scholar] [CrossRef]
- Coelho, E.A.A.; Cortizo, P.C.; Garcia, P.F.D. Small-signal stability for parallel-connected inverters in stand-alone AC supply systems. IEEE Trans. Ind. Appl. 2002, 38, 533–542. [Google Scholar] [CrossRef]
- Agorreta, J.L.; Borrega, M.; López, J.; Marroyo, L. Modeling and Control of N-Paralleled Grid-Connected Inverters With LCL Filter Coupled Due to Grid Impedance in PV Plants. IEEE Trans. Power Electron. 2011, 26, 770–785. [Google Scholar] [CrossRef]
- Bottrell, N.; Prodanovic, M.; Green, T.C. Dynamic Stability of a Microgrid With an Active Load. IEEE Trans. Power Electron. 2013, 28, 5107–5119. [Google Scholar] [CrossRef] [Green Version]
- Uudrill, J.M. Dynamic Stability Calculations for an Arbitrary Number of Interconnected Synchronous Machines. IEEE Trans. Power Appar. Syst. 1968, PAS-87, 835–844. [Google Scholar] [CrossRef]
- Gaba, G.; Lefebvre, S.; Mukhedkar, D. Comparative analysis and study of the dynamic stability of AC/DC systems. IEEE Trans. Power Syst. 1988, 3, 978–985. [Google Scholar] [CrossRef]
- Arabi, S.; Rogers, G.; Wong, D.; Kundur, P.; Lauby, M. Small signal stability program analysis of SVC and HVDC in AC power systems. IEEE Trans. Power Syst. 1991, 6, 1147–1153. [Google Scholar] [CrossRef]
- Hou, P.; Ebrahimzadeh, E.; Wang, X.; Blaabjerg, F.; Fang, J.; Wang, Y. Harmonic stability analysis of offshore wind farm with component connection method. In Proceedings of the IECON 2017—43rd Annual Conference of the IEEE Industrial Electronics Society, Beijing, China, 29 October–1 November 2017; IEEE: Beijing, China, 2017; pp. 4926–4932. [Google Scholar] [CrossRef]
- Wang, Y. Stability Assessment of Inverter-Fed Power Systems. Ph.D. Thesis, Aalborg University, Aalborg, Denmark, 2017. [Google Scholar]
- Rosso, R.; Engelken, S.; Liserre, M. Analysis of the Behavior of Synchronverters Operating in Parallel by Means of Component Connection Method (CCM). In Proceedings of the 2018 IEEE Energy Conversion Congress and Exposition (ECCE), Portland, OR, USA, 23–27 September 2018; IEEE: Portland, OR, USA, 2018; pp. 2228–2235. [Google Scholar] [CrossRef]
- Rosso, R.; Engelken, S.; Liserre, M. Analysis of the Parallel Operation Between Synchronverters and PLL-Based Converters. In Proceedings of the 2019 IEEE Energy Conversion Congress and Exposition (ECCE), Baltimore, MD, USA, 29 September–3 October 2019; IEEE: Baltimore, MD, USA, 2019; pp. 2583–2590. [Google Scholar] [CrossRef]
- 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]
- Dimitrovski, R.; Mehlmann, G.; Luther, M.; Ndreko, M.; Winter, W.; Trunk, K.; Rachinger, S. Comprehensive computer program for SSTI analysis of AC/DC power systems based on the component connection method. Int. J. Electr. Power Energy Syst. 2019, 111, 508–516. [Google Scholar] [CrossRef]
- Dragičević, T.; Lu, X.; Vasquez, J.C.; Guerrero, J.M. DC Microgrids—Part I: A Review of Control Strategies and Stabilization Techniques. IEEE Trans. Power Electron. 2016, 31, 4876–4891. [Google Scholar] [CrossRef] [Green Version]
- Dragičević, T.; Lu, X.; Vasquez, J.C.; Guerrero, J.M. DC Microgrids—Part II: A Review of Power Architectures, Applications, and Standardization Issues. IEEE Trans. Power Electron. 2016, 31, 3528–3549. [Google Scholar] [CrossRef] [Green Version]
- Mandrile, F.; Carpaneto, E.; Bojoi, R. VSG Simplified Damper Winding: Design Guidelines. In Proceedings of the IECON 2019—45th Annual Conference of the IEEE Industrial Electronics Society, Lisbon, Portugal, 14–17 October 2019; IEEE: Lisbon, Portugal, 2019; pp. 3962–3967. [Google Scholar] [CrossRef]
- Teodorescu, R.; Liserre, M.; Rodriguez, P. Grid Converters for Photovoltaic and Wind Power Systems; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
- Kjaer, P.C.; Gupta, M.; Martinez, A.; Saylors, S. Comparison of Experiences with Wind Power Plants with Low SCR. In Proceedings of the IEEE PES General Meeting, Denver, CO, USA, 26 March 2015. Technical Panel. [Google Scholar]
Figure 1.
Modular state-space modeling of a grid connected converter including controllers.
Figure 1.
Modular state-space modeling of a grid connected converter including controllers.
Figure 2.
Applications and trends of CCM in time.
Figure 2.
Applications and trends of CCM in time.
Figure 3.
Diagram of the system to be modeled.
Figure 3.
Diagram of the system to be modeled.
Figure 4.
Block diagram of the modeled S-VSC [
23] control with the physical system.
Figure 4.
Block diagram of the modeled S-VSC [
23] control with the physical system.
Figure 5.
LCL model in the stationary (,) frame.
Figure 5.
LCL model in the stationary (,) frame.
Figure 6.
Block diagram of the PI regulator and equivalent delay (PI + delay) of the digital control and modulation.
Figure 6.
Block diagram of the PI regulator and equivalent delay (PI + delay) of the digital control and modulation.
Figure 7.
Vector diagram of the grid voltage and the (d,q) S-VSC rotating reference frames. The angle is defined as the load angle variation after a fault event.
Figure 7.
Vector diagram of the grid voltage and the (d,q) S-VSC rotating reference frames. The angle is defined as the load angle variation after a fault event.
Figure 8.
Pole map of the S-SVC. The poles are numbered according to
Table 1.
Figure 8.
Pole map of the S-SVC. The poles are numbered according to
Table 1.
Figure 9.
Pole map of the S-SVC when the SCR changes from 10 to 1.5. The arrows indicate a decrease in the SCR. The poles are numbered according to
Table 1. From left to right: (
a) Complete pole map; (
b) Magnifications of the lower frequency poles.
Figure 9.
Pole map of the S-SVC when the SCR changes from 10 to 1.5. The arrows indicate a decrease in the SCR. The poles are numbered according to
Table 1. From left to right: (
a) Complete pole map; (
b) Magnifications of the lower frequency poles.
Figure 10.
Experimental setup where the controller is based on dSPACE.
Figure 10.
Experimental setup where the controller is based on dSPACE.
Figure 11.
Test 1: Step in active power reference from 0.2 pu to 0.3 pu. From left to right: (a) S-VSC virtual rotor speed (Hz) from experimental test (Exp) and state-space model (SS); (b) Active power injected from the inverter (pu) from experimental test (Exp), filtered active power (Filt) and output of the state-space model (SS).
Figure 11.
Test 1: Step in active power reference from 0.2 pu to 0.3 pu. From left to right: (a) S-VSC virtual rotor speed (Hz) from experimental test (Exp) and state-space model (SS); (b) Active power injected from the inverter (pu) from experimental test (Exp), filtered active power (Filt) and output of the state-space model (SS).
Figure 12.
Test 2: Step in reactive power reference from 0.1 pu to 0.2 pu. Reactive power injected from the inverter (pu) from experimental test (exp), filtered reactive power (Filt) and output of the state-space model (SS).
Figure 12.
Test 2: Step in reactive power reference from 0.1 pu to 0.2 pu. Reactive power injected from the inverter (pu) from experimental test (exp), filtered reactive power (Filt) and output of the state-space model (SS).
Figure 13.
Test 3: Grid frequency step drop from 50 Hz to 49.8 Hz. From left to right: (a) S-VSC virtual rotor speed (Hz) from experimental test (Exp) and state-space model (SS); (b) Active power injected into the grid (pu) from experimental test (Exp), filtered active power (Filt) and output of the state-space model (SS).
Figure 13.
Test 3: Grid frequency step drop from 50 Hz to 49.8 Hz. From left to right: (a) S-VSC virtual rotor speed (Hz) from experimental test (Exp) and state-space model (SS); (b) Active power injected into the grid (pu) from experimental test (Exp), filtered active power (Filt) and output of the state-space model (SS).
Figure 14.
Test 4: Grid voltage drop % with phase jump . From left to right: (a) Reactive power injected into the grid (pu) from experimental test (Exp), filtered reactive power (Filt) and output of the state-space model (SS); (b) S-VSC virtual excitation flux from the experimental test (Exp) and state-space model (SS).
Figure 14.
Test 4: Grid voltage drop % with phase jump . From left to right: (a) Reactive power injected into the grid (pu) from experimental test (Exp), filtered reactive power (Filt) and output of the state-space model (SS); (b) S-VSC virtual excitation flux from the experimental test (Exp) and state-space model (SS).
Table 1.
Poles and participation factors. Only significant factors are displayed (larger than 0.1).
Table 1.
Poles and participation factors. Only significant factors are displayed (larger than 0.1).
Poles | 1–2 | 3–4 | 5–6 | 7–8 | 9–10 | 11–12 | 13–14 | 15 | 16 |
(Hz) | 2084.58 | 1988.07 | 1513.86 | 272.92 | 50.34 | 31.85 | 1.38 | 1.35 | 0.16 |
(ms) | 0.35 | 0.35 | 0.11 | 0.61 | 21.87 | 5.17 | 167.39 | 117.67 | 999.67 |
| 0.216 | 0.229 | ∼1 | 0.959 | 0.145 | 0.966 | 0.691 | – | – |
| 0.2065 | 0.2079 | 0.1532 | 0.2639 | – | – | – | – | – |
| 0.2065 | 0.2079 | 0.1532 | 0.2609 | – | – | – | – | – |
| – | – | – | – | 0.2740 | – | – | – | – |
| – | – | – | – | 0.2810 | – | – | – | – |
| 0.2471 | 0.2475 | – | – | – | – | – | – | – |
| 0.2471 | 0.2474 | – | – | – | – | – | – | – |
| – | – | – | – | – | 0.5750 | – | – | – |
| – | – | – | – | – | 0.5759 | – | – | – |
| – | – | 0.6859 | 0.1876 | – | – | – | – | – |
| – | – | 0.6863 | 0.1875 | – | – | – | – | – |
| – | – | – | 0.4200 | 0.2282 | – | – | – | – |
| – | – | – | 0.4178 | 0.2199 | – | – | – | – |
| – | – | – | – | – | – | 2.2608 | 5.5357 | – |
| – | – | – | – | – | – | 1.6308 | −2.2624 | – |
| – | – | – | – | – | – | 1.6362 | −2.2699 | – |
| – | – | – | – | – | – | – | – | 1.0002 |
Table 2.
Inverter and S-VSC parameters for the state-space model and of the experimental setup.
Table 2.
Inverter and S-VSC parameters for the state-space model and of the experimental setup.
Parameter | Value | Parameter | Value |
---|
| V | | 15 kVA |
| 2 mH | | 5 F |
| 1 mH | | 3 mH |
| 3.77 | | 710.6 /s |
| 0.1 pu | | 0.02 pu |
H | 4 s | | 0.22 |
| 1.048 pu | | 0.278 s |
| 10 kHz | | 700 V |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).