Sensitivity Analysis and Small-Signal Stability of Grid Following Converters

Análisis de sensibilidad y estabilidad de pequeña señal en convertidores seguidores de red

DOI 10.22430/22565337.2383 Logotecnologicas PDF Table Figure

Received: April 6, 2022
Accepted: June 6, 2022
Available: July 1, 2022

How to cite / Cómo citar
S. Sepúlveda-García; A. Garcés-Ruíz, “Sensitivity Analysis and Small-Signal Stability of Grid Following Converters,” TecnoLógicas, vol. 25, nro. 54, e2383, 2022. https://doi.org/10.22430/22565337.2383

Abstract

This article presents a small-signal study, based on a sensitivity analysis, of a voltage source converter (VSC) operating as a grid-following device. This type of operation is standard for integrating distributed energy resources directly into a primary grid without considering the control over specific variables (e.g., voltage and frequency), thus maximizing the amount of injected active power. The present study adopted the classical vector-oriented control based on the linearization method of the averaged model; hence, this analysis is limited to linear controls. The main objective of this study is to evaluate how each controller parameter influences the system’s stability. To conduct the sensitivity analysis, an averaged model in the dq reference frame of the VSC was employed to describe the system dynamics in the equivalent state-space model. Afterward, the stability was analyzed through a sensitivity analysis of the eigenvalues of the corresponding state matrix. The numerical results demonstrate that the main problems for the stability of VSCs operating as grid-following converters are threefold: a high value of the filter inductance, a non-ideal impedance that appears depending on the connection point, and poor coordination of the parameters of the Kp controllers. These results are compared to a simple bifurcation analysis of the state matrix, which consists of a diagram that describes the variation in eigenvalues when a parameter changes, thus proving their validity.

Keywords: Grid-following converters, Hopf bifurcation, Sensitivity analysis, Vector-oriented control, Voltage source converter.

Resumen

El presente artículo presenta un estudio de pequeña señal, basado en el análisis de sensibilidad, del convertidor Voltage Source Converter (VSC) operando como dispositivo de seguimiento de red. Este tipo de operación es común para la integración de recursos energéticos distribuidos directamente a la red primaria sin considerar el control sobre variables específicas como la tensión o la frecuencia, maximizando así la cantidad de potencia activa inyectada. El presente estudio adopta el clásico Control Vectorial Orientado basado en la linealización del modelo promediado, por tanto, este análisis es limitado a controles lineales. El objetivo principal de esta investigación consistió en evaluar cómo influye cada parámetro del controlador en la estabilidad del sistema. Con el fin de realizar el análisis de sensibilidad, se empleó un modelo medio en el marco de referencia dq del VSC para describir la dinámica del sistema en el espacio de estados equivalentes; posteriormente, se buscó la estabilidad mediante un análisis de sensibilidad de los valores propios de la matriz de estado correspondiente, cuyos resultados numéricos demostraron que los principales problemas para la estabilidad de los VSC que operan como convertidores de seguimiento son tres: un valor elevado de la inductancia del filtro, la presencia de una impedancia no ideal que aparece en función del punto de conexión, y una mala coordinación de los parámetros de los controladores Kp. Estos resultados se han comparado con un simple análisis de bifurcación de la matriz de estado, que consiste en un gráfico que describe el movimiento de los valores propios cuando cambia un parámetro, corroborando así su validez.

Palabras clave: Convertidores seguidores de red, Bifurcación tipo Hopf, Análisis de sensibilidad, Control vectorial orientado, Convertidor de fuente de tensión.

1. INTRODUCTION

Power electronic converters are key elements in the modern electrical grid, they are immerse from applications of High Voltage Direct Current applications (HVDC) [1] to Low Voltage Direct Current (LVDC) [2], [3] applications. Moreover, the penetration of Distributed Energy Resources (DERs) into the grids make a strongly dependence of this devices [4].

Those DERs usually operate in DC, and hence requires a converter, named Voltage Source Converter (VSC), to integrate the renewable resource into the primary grid [5]. Although this is a relatively new technology in Colombia, grid operators have an increasing interest in studying the impact of these devices on the power distribution grid. There is a consensus in the necessity of diversify the generation matrix, and increase the penetration of DERs into the grid which naturally lead to an increase in the presence of VSCs [6]. The international community identified that the stability analysis of inverter-based systems is a fundamental challenge for the massive integration of renewable energy in the near future, as demonstrated in the most recent document of the IEEE-PES Task Force on Microgrid Stability Definitions, Analysis, and Modeling [7].

Power electronic converters introduce new dynamic phenomena compared to conventional power systems dominated by the dynamic behavior of synchronous machines. The issues related to the stability of these new grids are described in [7], [8], where it is clear that the control of the converter is a fundamental element in the study of small-signal stability [9].

The classical vector-oriented control employed for electrical machines has been adapted to be applicable in power electronic converters and nowadays is one of the most commercial [10], [11]. This control has several loops that interact together, generating complex dynamics, however, in order to simplify the analysis, some studies make an approximation by analyzing the average model of the VSC, which only takes into account the dynamic response and not the complete model which includes the switch losses [9].

Several studies have been conducted to analyze the stability of the voltage source converter, such as the unstable phenomenon in the three-phase case [12]. Furthermore, it has been reported that one of the main problems affecting the small-signal stability of the VSC system is the influence of a non-ideal grid in the system [13]–[15]. This study has different approaches for the analysis of the stability. In [16], [17] it was analyzed directly in the Nonlinear equivalent; however, this approach is not suitable for large systems. The passivity theory is used in [18] in order to study the small-signal stability of the VSC. A recent approach made use of the Graph theory together with the state space [19], this approach demonstrates to be more robust and faster than the studies previously described.

However, it is still preferred to use linearization methods due to its simplicity to calculate the instability [20]. Recently there has also been an interest in stability analysis of the single-phase converter due to its wide applications in electrical systems [14], [21].

Nevertheless, sensitivity analysis is rarely discussed for these converters, despite its usefulness in determining the elements and parameters of the system that cause the most problems to the stability in a numerical way [22]. Therefore, this article aims to fill this gap, showing a sensitivity analysis for the VSC and determining how the parameters affect the system's stability. The stability analysis of inverter-based power systems is classified into control-system stability and power supply balance stability [23] The former is divided into the electric machine and converter stability, whereas the latter is divided into voltage and frequency stability. The study presented here devises the converter stability, although it also concerns voltage stability on the dc-link, both from the small-signal perspective. As mentioned before this approach is limited to controllers based on average models, other controllers such as the Sliding Mode control are more robust, however some industrial applications still make use of PIs controllers in cascade with the simplified average models due to its ease implementation [24].

The rest of the paper is organized as follows: First, the converter model is described mathematically in Section II. Next, in Section III, the mathematical model is linearized around the operating point. Finally, in Section IV, a sensitivity analysis of the system's elements is performed. Conclusions are presented in Section V, followed by a brief list of references. This article was previously introduced at the international conference X-Sicel.

2. METHODOLOGICAL ASPECTS

2.1 Model of the VSC

The employed model of the VSC is a simplified model with an ideal grid and switching losses close to zero, this approximation is known as average model, and provides a fair estimation of the stability of the VSC; for more robust or detailed applications it is necessary to make use of different type of controllers. The model is continuous with a modulation index mi ∈ [-1,1]. The control to be studied is the conventional vector-oriented control, the most common approach in commercial converters. The circuit, together with its control, is described in Figure 1 [25].

Voltage source converter (VSC) with vector-oriented control (VOC)
Figure 1. Voltage source converter (VSC) with vector-oriented control (VOC). Source: Created by the authors.

2.1.1 Model of the grid

Kirchhoff’s second law is applied to the loop comprised of the converter output, the filter, and the point of common coupling (PCC), as given in (1) and (2).

                   (1)

                   (2)

All variables are evident from the figure above. Note that the system is described in the dq frame because the variables are easier to control by PI-type controllers in a stationary frame of reference than in a rotating one. Similarly, the relationship for the non-ideal grid is given in (3) and (4).

                   (3)

                   (4)

From the set of equations described above, it is possible to demonstrate that the tensions in the point of common coupling are (5) and (6), respectively.

                   (5)

                   (6)

Furthermore, it is possible to demonstrate that the dynamics of the currents in the ac side in the dq frame are (7) and (8).

                   (7)

                   (8)

The set of equations describing the system's dynamics is completed with the voltage dynamics at the DC side (9). This voltage can be obtained through the power balance as follows.

                   (9)

The vector-oriented control has a structure with different modules named phase-locked-loop (PLL), inner-loop, and outer-loop. Each of these modules is briefly described below.

2.1.2 Phase-Locked-Loop (PLL)

The phase-locked-loop is in charge of coordinating the converter with the pre-existing network; the block diagram is described in Figure 2.

Phase Locked Loop (PLL)
Figure 2. Phase Locked Loop (PLL). Source: Created by the authors.

By defining an auxiliary variable γpll. The model can be mathematically rewritten at the PI controller's output [26], as shown in (10) and (11).

                   (10)

                   (11)

2.1.3 Inner Loop

The inner loop controls the currents of the converter on the AC side; the block diagram is shown in Figure 3.

Inner Loop
Figure 3. Inner Loop. Source: Created by the authors.

Similarly, a set of auxiliary variables γd and γq are defined to describe the model mathematically as (12)-(15).

                   (12)

                   (13)

                   (14)

                   (15)

The variables with a bar, such as x ̅ are references. Equations (5) and (6) are replaced into (10) and (11), so it is possible to obtain (16) and (17).

                   (16)

                   (17)

That represents the control variables of the system.

2.1.4 Outer loop

Finally, the outer loop is in charge of controlling the DC voltage of the system; the block diagram is described in Figure 4.

Outer Loop
Figure 4. Outer Loop. Source: Created by the authors.

Which can be described as (18)-(19).

                   (18)

                   (19)

Since the unitary power factor is assumed, there is no control over reactive power.

2.2 State-Space Model

Equations (7) and (19) are collected. Observe that the state variables are x=[id, iq, vdc, γd, γq, γv ]T, the control variables y=[vds, vqs ]T, and the input variables u=[vdg, vdc, Idc, iq*]T. Note that the variables θ' and γpll are not considered since, in this representation, there are no control and system state variables in the PLL model and vice versa; therefore, the frequency is supposed to be ideal. These equations follow the form described in (20) and (21).

                   (20)

                   (21)

The aim is to obtain a system of linear equations like (22), which may be linearized around the equilibrium point, using the expressions described in (23) and (24) [27].

                   (22)

                   (23)

                   (24)

The system's initial conditions can be obtained by considering stable operation (i.e., the dynamics are zero; therefore, the derivatives are equal to zero), these values are given in Table 1.

Table 1. Initial conditions
Source: Created by the authors.
Initial conditions

The results of the linearization around the operating point are given by (25), where the elements J3 of the third row are given by (26)-(31).

                   (25)

                   (26)

                   (27)

                   (28)

                   (29)

                   (30)

                   (31)

2.2.1 Sensitivity Analysis

The sensitive analysis aims to determine the influence of each element of the state matrix on the system's eigenvalues. That is to say, to evaluate how a change in the value of one component could affect the position of the eigenvalues. The participation factor pi expresses this relation which indicates the rate of change of a specific eigenvalue against an element of the state matrix as shown in (32).

                   (32)

Where akm corresponds to the entry in row k, with column m; φi corresponds to the eigenvectors on the left and ϕi to the eigenvectors on the right. Notice that (25) exposes the variation of the eigenvalue λi concerning the element akm; likewise, it is possible to observe that the final result is a three-dimensional matrix, known as a tensor.

It is difficult to visualize the results in a tensor. Therefore, the influence of a specific parameter (e.g., the inductance), immersed in different elements of the state matrix, is obtained by applying the chain rule in (32) , resulting in (33).

                   (33)

Where x corresponds to the parameter whose influence on the oscillation mode i is to be known, finally, it is necessary to clarify that the numerical value obtained from the participation factors does not have much relevance. However, their relative value does. It is taken to measure how much a variation affects or not one or another oscillation mode.

3. RESULTS AND DISCUSSION

The stability analysis was performed in a realistic system with system parameters and constants shown in Tables (2) and (3), respectively.

Table 2. System Parameters
Source: Created by the authors.
Parameter Value
Capacitance (C) 2200 μF
Filter inductance (Lf) 55 mF
Filter Resistance (Rf) 50 mΩ
Source inductance (Lg) 10 mH
Source Resistance (Rg) 20 mΩ
Source voltage (vdg) 380 V
DC voltage (vdc) 1000 V
DC current (Idc) 2 A

Table 3. Control Parameters
Source: Created by the authors.
Parameter Value
kpi 29.33
kpv 1.1729
kii 7106
kiv 312.66

The model was implemented in OpenModelica and was compared with the linearization around the operation point, proving a good approximation, as depicted in Figure 5. Note that it is not strictly necessary to use per unit variables.

Comparison between the statespace model and the complete model
Figure 5. Comparison between the statespace model and the complete model. Source: Created by the authors.

The system's stability can be analyzed by positioning the eigenvalues of the state matrix as mentioned in (34).

                   (34)

For the base case, they take the value described in (35).

                   (35)

This equation assumes that the system operates in steady-state; however, we want to answer the question: how do the system parameters affect the system's stability?

This question can be solved in the first instance by applying (25) where the generic case (i.e., a three-dimensional matrix) is obtained; however, visualizing the results is somewhat complex. With the help of (26), it is possible to get more visible results. Note that each numerical value in Table 4 and Table 5 indicates the rate of change of the eigenvalue concerning the system parameter.

Table 4. Participation factors (real part) system parameters
Source: Created by the authors.
Parameter Lf C
p1 3796.68 -4171.90
p2 3796.68 -4171.90
p3 1057.86 2245.72
p4 1057.86 2245.72
p5 4854.54 -1.14 × 10-27
p6 4854.54 -1.14 × 10-27

Table 5. Participation factors (real part) control parameters
Source: Created by the authors.
Parameter kpi kpv kii kiv
p1 -6.83 -91.424 -0.0035 0.371
p2 -6.83 -91.424 -0.0035 0.371
p3 -3.65 56.75 0.0035 -0.37
p4 -3.65 56.75 0.0035 -0.37
p5 -9.09 3.92 × 10-30 -1.911 × 10-19 -3.49 × 10-32
p6 -9.09 3.92 × 10-30 -1.911 × 10-19 -3.49 × 10-32

As can be deduced, the participation factors, like the eigenvalues, come in conjugate pairs. As mentioned, their numerical value does not have that much impact; what matters are the associated relative values (as well as their sign). It is possible to deduce, for example, that a slight change in the inductance value will noticeably affect the positioning of all the eigenvalues (especially λ1,2,5,6), reducing the stability margins. While in the case of an increase in the capacitance value, the stability margin of the eigenvalues λ1,2 is increased at the cost of a reduction in the margins of the eigenvalues λ4,5.

Contrary to the case of the inductance, an increase in the constant kpi will increase the stability margins of the system; it is also noticeable how the real part of the eigenvalues λ5,6 is only influenced by the system's inductance and the current controller's proportional action.

It is also remarkable how the increase of the other constants can improve and worsen the positioning of some eigenvalues.

Similarly, the imaginary part of the eigenvalues is essential since they affect the so-called damping ratios described in (36). A value lower than 10% of these damping ratios affects the system, and control processes must be developed to increase them.

                   (36)

Hence it is worth asking, how do changes in controller parameters and constants affect the imaginary part? The answer to this question may be obtained from Tables 6 and 7.

Table 6. Participation factors (imaginary part) system parameters
Source: Created by the authors.
Parameter Lf C
p1 -690.37 -96230.6
p2 690.37 96230.6
p3 -777.58 22583.16
p4 777.58 -22583.16
p5 16.01 3.44×10-27
p6 -16.01 -3.44×10-27

Table 7. Participation factors (imaginary part) control parameters
Source: Created by the authors.
Parameter kpi kpv kii kiv
p1 0.26 149.74 0.018 0.115
p2 -0.26 −149.74 -0.018 -0.115
p3 −3.14 −72.74 0.013 0.113
p4 3.14 72.74 -0.013 -0.113
p5 −8.92 3.86×10-30 0.033 -2.20×10-32
p6 8.92 -3.86×10-30 -0.033 2.20×10-32

As aforementioned, eigenvalues come in conjugate pairs, indicating that a change in each parameter will move the eigenvalues closer to (or further away from) the real axis depending on the case. For example, an increase in kpi is accompanied by an increase in the damping ratios of the system. Whereas for kii an increase moves the eigenvalues away from the real axis, implying that the damping ratios decrease significantly.

These results described above can be easily verified through a simple bifurcation analysis, where (34) is used, and a given parameter is varied, plotting the path of the eigenvalues as shown in Figure 6 and Figure 7.

Bifurcation of the system for L 50150 mH
Figure 6. Bifurcation of the system for L 50150 mH . Source: Created by the authors.

Bifurcation of the system for kii [5-100]×103
Figure 7. Bifurcation of the system for kii [5-100]×103. Source: Created by the authors.

As predicted from the analysis of the participation factors, an increase in the inductance value reduces the stability margin until a Hopf bifurcation occurs. Moreover, increasing the integral action of current control increases damping ratios. Consequently, the eigenvalues λ5,6 moves away from the real axis.

Finally, it is interesting to know what happens when a non-ideal network is in the system. In this scenario, the participation factors described in Table 8 are presented.

Table 8. Participation factors (real part) system parameters
Source: Created by the authors.
Parameter Lf Lg
p1,2 3947.41±288.29 475.8±55.07
p3,4 1180.67±686.03 161.97±84.83
p5,6 4856.19±509.26 -3.27×10-30±6.66×10-31

As can be seen, these indicate that with the presence of the non-ideal inductance, the stability margins of the first four eigenvalues are reduced. In addition, it can be noticed how it slightly changes the participation factors of the inductance L., becoming more "sensitive" to changes, which happens with the other parameters of the network and control. Therefore, the presence of a non-ideal network significantly affects the system's stability, reducing the stability margins and making the system more sensitive to changes in other parameters.

4. CONCLUSIONS

In this paper, a sensitivity analysis is used to determine numerically how system parameters and controller constants influence the stability of a VSC operating as a grid feeder. As remarked, this analysis can be helpful since it can easily detect the parameters that cause stability problems, the proposed approach is based in the average model neglecting the losses in the grid and the switch losses. In the same way, these variables are suitable to control to improve the positions of some eigenvalues, increase margins of stability, damping ratios, or simply design the converter to maintain stability when they are fixed variables As may be deduced from the results, an increase in the proportional action of the current controller increases the stability margins of all the eigenvalues, which is helpful since it can be coordinated with a high inductance (or capacitance) value for filtering, which notoriously these margins. From the perspective of the control, it is desired to have a high value for the integral action since it ensures zero steady-state error; nevertheless, as can be deduced from the analysis, a high value of this parameter reduce significantly the damping ratios of the system, hence, a study of the control and design of the VSC must be done with a corresponding study of stability. For future works, it is recommended to include more detailed models of the power source and if necessary, the incorporation of DC-DC converter in the model. In addition, it is important to analyze the influence of the R/X ratio in the case of microgrids; furthermore, the sensitivity analysis can be used in a general way in the study of stability of any electrical system that could be represented as its equivalent state-space model, for this reason it is recommended to make a sensitivity analysis of a modern controllers such as the Non-Linears widely proposed in the technical literature.

5. ACKNOWLEDGMENTS

This paper is a result of project Integra2023, code 111085271060, contract 80740-774-2020, funded by the Ministry of Science, Technology, and Innovation (Minciencias) and developed by the ICE3 Research Group at Universidad Tecnológica de Pereira (UTP).

CONFLICTS OF INTEREST

The authors have no conflicts of interest to declare. All co-authors have seen and agree with the manuscript's contents, and there is no additional financial interest to report, different from Minciencias and Universidad Tecnológica de Pereira. We certify that the submission is original work and is not under review at any other publication.

AUTHOR CONTRIBUTIONS

Simón Sepúlveda-García: Methodology, Implementation, Writing.

Alejandro Garcés-Ruíz: Conceptualization, Methodology, Writing, Validation.

6. REFERENCES