Next Article in Journal
Animal Nutrition and Welfare in Sustainable Production Systems
Next Article in Special Issue
Optimal Power Flow in Wind–Photovoltaic Energy Regulation Systems Using a Modified Turbulent Water Flow-Based Optimization
Previous Article in Journal
The Path Constitution of Platform Evolution: An Organizational Momentum View
Previous Article in Special Issue
A Novel Fairness-Based Cost Model for Adopting Smart Charging at Fast Charging Stations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Small Signal Stability Analysis of a Microgrid in Grid-Connected Mode

Department of Electrical Engineering, College of Engineering, Jouf University, Sakaka 72388, Saudi Arabia
Sustainability 2022, 14(15), 9372; https://doi.org/10.3390/su14159372
Submission received: 4 July 2022 / Revised: 22 July 2022 / Accepted: 25 July 2022 / Published: 31 July 2022

Abstract

:
Microgrid stability issues are classified into three categories: transient, voltage, and small signal stability (SSS). Small variations in the load demand and small perturbations in the control system and line impedance parameters can cause instability, which can be avoided by performing an SSS analysis. This paper focuses on investigating the impact of line impedance and passive filter parameters on the stability of a MG in grid-connected mode. Therefore, a MG system was represented mathematically, before performing an SSS analysis that calculated the stability margin of the MG parameters. A sensitivity analysis was performed to determine those parameters highly participating in the SSS. The mathematical results were validated using the simulation results, which were obtained using MATLAB Simulink.

1. Introduction

A microgrid (MG) is a group of micro-sources and loads in a subsystem, which can work in an islanding mode or in a grid-connected mode. If the MG is in grid-connected mode, its frequency and voltage will be dominated by the grid, which makes stability less crucial than that when it is in the islanding mode, where the variation of load demands causes major disturbances to the frequency and voltage, which can lead to destabilization [1,2].
Stability issues in MGs are classified as small signal, transient, and voltage stability. Small signal stability (SSS) is related to the feedback gains of controllers, changes in power demands, and small perturbations in system parameters. A MG is transiently stable if it reaches the steady state condition after large disturbances such as faults, large load steps, and switching to the islanding mode [3].
The power output of renewable energies changes with time, which is considered a perturbation to the system. This power variation makes it challenging for MGs to produce reliable and stable power, especially when there is variation in the load. Therefore, it is crucial to study the impact of small power variations on MG stability, which are considered SSS [4].
SSS can be studied using the eigenvalue theorem, where the system is unstable if at least one of the eigenvalues is on the right half plane of an imaginary axis. However, knowing the eigenvalues is not possible without obtaining mathematical models of the whole system. Therefore, a MG system should be represented mathematically by linear equations, before forming the state space model that allows stability analysis [5,6].
Power demand variations and parameter perturbations are considered SSS issues that can affect the reliability of MGs. Therefore, it is crucial to provide accurate mathematical models for MGs to analyze their stability conditions with respect to small changes to MGs parameters and power demand variations [4]. Power systems behavior is represented by non-linear differential equations. Thereafter, the stability condition can be tested by using eigenvalue analysis, which can be compared with time domain simulations. The accuracy of the mathematical model representing a MG can be measured by the percentage of matching between the analytical results and time domain results. A system is considered stable if it is able to function again around the same operating point after a disturbance. For example, the voltage of a system is stable if the voltage at all buses is settled after any kind of disturbance. If the voltage becomes unstable, its behavior will be either falling or rising [6,7].
When a MG is in grid-connected mode it is dominated by the grid, which makes it less possible for the MG to be unstable. Stability is a big concern when a MG is in islanding mode, where stability is highly related to control gains, load fluctuations, and line impedance variation [8].
The authors of [9] studied the SSS of a MG in islanding mode, with two distributed generators. The parameters of current, voltage, and LC passive filter were varied in small steps to study the impact on stability. In [10], stability state was examined by applying a load step to a MG in the islanding mode. A mathematical model, simulation model, and experiment setup were provided, which showed the same results. The authors of [11] used SSS analysis to study the impact of droop control gain variation and load variation on the stability of a MG, with three distributed generators in islanding mode. Similarly, in [12] SSS analysis was performed on a MG with two distributed generators in islanding mode, with respect to droop control gains and load variation. The authors of [13] developed a small signal model of a MG in islanding mode, to analyze the stability and robustness of the MG, and then an adaptive control method was proposed to adjust the control gains in response to disturbances. In [14,15,16,17,18,19,20] the SSS of MGs in islanding mode was studied when controller gains were increased. In [18], the SSS of a MG in islanding mode was resolved, while proposing a method that simultaneously ensures stability and optimal power flow. Most research related to the SSS of MGs in islanding mode has focused on designing controllers to manage the variation of the load demands, while ensuring stability of the MG system [19].
The authors of [21] used an SSS analysis to examine the impact of droop control gains on the stability of a MG in grid-connected mode. SSS analysis was implemented in [22] to discover the stability margin of control gains in a MG in grid-connected mode. The authors of [23] used an SSS analysis to study the stability of a MG in grid-connected mode, with respect to load demand variations. The authors of [24] derived the SSS limits of a phase locked loop, responsible for synchronizing a wind farm with the grid. The authors of [25] provided a small signal model of a MG with dynamic loads and transmission lines, before investigating the parameter sensitivity and proposing methods to the enhance stability of MGs under different scenarios and configurations.
It is important to study any possible disturbance of MGs, especially in the islanding mode, because the inertia is low and MGs are normally associated with renewable energy sources. It is assumed that MGs are stable when they are in the grid-connected mode, as they will be dominated by the grid. Therefore, previous research has focused on the SSS analysis of MGs in islanding mode with respect to control gains and load demand variations. Only a few articles discussed SSS when the MG is in grid-connected mode, as the MG will be dominated by the grid and small disturbances to control gains and load demands do not normally cause instability. Furthermore, there are scarce published data on investigating the SSS of MGs in grid-connected mode, with respect to LCL filter and line impedance parameter variations. Therefore, it was decided in this article to investigate the SSS of a MG when it is in grid-connected mode.
The main contributions of this article are summarized as follows:
  • Deriving a state space model of a MG in grid-connected mode, which involves a LCL filter, line impedance, and control system.
  • Investigation of the impact of the LCL filter and line impedance parameters on the stability of a MG in grid-connected mode.
  • Sensitivity analysis of MG parameters was performed to investigate the most sensitive parameters that affect the SSS of a MG in grid-connected mode.
This paper is organized as follows: Section 2 presents the MG system and its parameters. Section 3 presents the state space of the MG system. Section 4 reports the simulation model that was modelled in MATLAB, while Section 5 shows the analysis and results of the MG model.

2. Proposed Controller Design

A MG supporting a three-phase grid with the proposed control method is represented in Figure 1, and the system parameters are detailed in Table 1. The figure shows that the MG is in grid-connected mode. The input is a DC voltage that is converted to AC through the inverter. However, the output voltage of the inverter is mostly a square wave, which requires a filter to convert the output square wave into a sine wave. The LCL filter is mainly used to interface the MG with the grid and to reduce the ripples produced by the high switching frequency of the inverter.
A three-phase voltage and current should be transformed into a two-phase orthogonal rotating reference frame before entering a control system. The control system is working in the rotating reference frame (dq) and is then transferred back to the three-phase (abc) before entering the VSI, to appear at the input of the LCL filter. The reference active and reactive power are considered as inputs to the control system, where the inverter output power is controlled by controlling the inverter output current, while the inverter output voltage should be fixed. The active power reference is set to 500 W, while the reactive power is set to zero, in order to make the power factor of the MG be at unity. Therefore, the total reactive power demand should be covered by the grid.
The droop controller is used to generate the reference currents based on the input reference values of the power. Therefore, the droop controller is responsible for meeting any loads connected to the MG. The current controller is a PI controller, which is responsible for comparing between the measured inverter output currents and the reference currents produced by the droop controller. The reference currents are adjusted by the PCC voltage and power references. The output of the PI current controller is a voltage that is fed to the PWM generator to be compared with the carrier signal, to generate the appropriate switching sequence and to allow controlling the voltage at the output of the inverter.
A phase locked loop (PLL) is used to decrease the frequency deviation between the grid and MG. Frequency deviation can severely affect the power sharing accuracy of a MG in grid-connected mode. The PLL is responsible for tracking the phase angle of the grid with the phase angle of the voltage at the PCC. Integrating the output frequency of the PLL provides the reference phase angle that is required for the transformation blocks.

3. State Space Model

To study the small signal stability of the system shown in Figure 1, the system is divided into five submodules. Each submodule is linearized and represented in the state space form. The general form of state space models is provided in [5,8]. The A matrix, which is the state matrix, of the whole system is formed by the combination of all of the submodules on a common reference frame. The stability was tested by analyzing the eigenvalues resulting from the combined system A matrix.

3.1. Phase-Locked Loop

The MG has to be synchronized with the main grid, as a lack of synchronization may cause large transient currents to occur at the point of common coupling (PCC). To ensure the synchronization process, the grid voltage amplitude, phase angle, and frequency should be accurately matched. This process is performed by the PLL.
The PLL in Figure 2 is used to track the MG frequency with the grid frequency and provide the transformation angle. The dynamics of the PLL are represented in (1)–(3). The small signal linearized state space model of the PLL is represented in (4) and (5).
It is assumed that Φ P L L = v o q d t to simplify the generation of the state-space model.
˙ P L L = k p , P L L v o q + k i , P L L Φ P L L
Φ ˙ P L L = v o q
ω P L L = k p , P L L v o q + k i , P L L Φ P L L
[ Δ ˙ P L L Δ Φ ˙ P L L ] = A P L L [ Δ P L L Δ Φ P L L ] + B P L L [ Δ i l d q Δ v o d q Δ i o d q ]
[ Δ ω P L L ] = [ C P L L ] [ Δ P L L Δ Φ P L L ] + D P L L [ Δ i l d q Δ v o d q Δ i o d q ]
where
A P L L = [ 0 k i , P L L 0 0 ] ,   B P L L = [ 0 0 0 k p , P L L 0 0 0 0 0 1 0 0 ] C P L L = [ 0 k i , P L L ] ,   D P L L = [ 0 0 0 k p , P L L 0 0 ]

3.2. Droop Controller

Figure 3 shows the power control method, where the reference frequency is compared with the output frequency of the PLL and the error signal is passed through the droop frequency gain and added to the commanded active power. On the other side, the d-component of the inverter output voltage at the PCC is compared with the reference voltage, and the error signal is passed through the droop voltage gain and added to the commanded reactive power. The outputs of the droop controller are the current reference values. The dynamics of the droop controller are represented in (6)–(9). The small signal linearized state space model of the droop controller is represented in (10) and (11).
P ˙ = P * + m ( ω * ω P L L )
Q ˙ = Q * + n ( E * v o d )  
i l d * = 2 3 v o d P ˙ = 2 3 v o d ( P * + m ( ω * ω P L L ) )  
i l q * = 2 3 v o d Q ˙ = 2 3 v o d ( Q * + n ( E * v o d ) )  
[ Δ P ˙ Δ Q ˙ ] = A D [ Δ P Δ Q ] + B D 1 [ Δ i l d q Δ v o d q Δ i o d q ] + B D 2 [ Δ ω P L L ]  
[ Δ i l d * Δ i l q * ] = C D [ Δ P Δ Q ] + D D 1 [ Δ i l d q Δ v o d q Δ i o d q ] + D D 2 [ Δ ω P L L ]  
where
A D = [ 0 0 0 0 ] ,   B D 1 = [ 0 0 0 0 0 0 0 0 n 0 0 0 ] B D 2 = [ m 0 ] ,   C D = [ 0 0 0 0 ] D D 1 = [ 0 0 2 3 v o d 2 ( P * + m ( ω * ω P L L ) ) 0 0 0 0 0 2 3 v o d 2 ( Q * + n ( E * ) ) 0 0 0 ] D D 2 = [ 2 m 3 v o d 0 ]

3.3. Current Controller

Two PI current controllers are shown in Figure 4, and they are used to provide the output voltage of the inverter. The controllers compare the measured currents at the coupling inductor and the reference currents, which are the outputs of the previous submodule. The dynamics of the current controller are represented in (12)–(15). The small signal linearized state space model of the current controller is represented in (16) and (17).
γ ˙ d = i l d * i l d
γ ˙ q = i l q * i l q
v i d * = ω n L c i l q + k p c ( i l d * i l d ) + k i c γ d
v i q * = ω n L c i l d + k p c ( i l q * i l q ) + k i c γ q
[ Δ γ ˙ d Δ γ ˙ q ] = A C [ Δ γ d Δ γ q ] + B C 1 [ Δ i l d * Δ i l q * ] + B C 2 [ Δ i l d q Δ v o d q Δ i o d q ]
[ Δ v i d * Δ v i q * ] = C C [ Δ γ d Δ γ q ] + D C 1 [ Δ i l d * Δ i l q * ] + D C 2 [ Δ i l d q Δ v o d q Δ i o d q ]
where
A C = [ 0 0 0 0 ] ,   B C 1 = [ 1 0 0 1 ] B C 2 = [ 1 0 0 0 0 0 0 1 0 0 0 0 ] ,   C C = [ k i c 0 0 k i c ] D C 1 = [ k p c 0 0 k p c ] ,   D C 2 = [ k p c ω n L c 0 0 0 0 ω n L c k p c 0 0 0 0 ]

3.4. LCL Filter

The switches of the inverter are assumed to be ideal and they produce no losses, which makes the inverter voltage the same as the demanded voltage v i = v i * . The dynamics of the LCL filter with the feeder are represented in (18)–(23). The PLL frequency, which is the inverter frequency, is considered to be the common reference frame. Therefore, all of the equations are represented in the inverter reference frame, except that the voltage at the grid bus needs to be converted according to the transformation matrix, as in (24) and (25). The grid voltage is transformed to the common reference frame according to (26), where the transformation angle δ represents the phase angle difference between the grid voltage and the inverter voltage. Finally, the equations are linearized and represented in the state-space form. The small signal linearized state space model of the LCL filter is represented in (27).
i ˙ l d = R c L c i l d + ω P L L i l q + 1 L c v i d 1 L c v o d
i ˙ l q = R c L c i l q ω P L L i l d + 1 L c v i q 1 L c v o q
v ˙ o d = ω P L L v o q + 1 C f i l d 1 C f i o d
v ˙ o q = ω P L L v o d + 1 C f i l q 1 C f i o q
i ˙ o d = ( R g + R f ) L g + L f i o d + ω P L L i o q + 1 L g + L f ( v o d v g d )
i ˙ o q = ( R g + R f ) L g + L f i o q ω P L L i o d + 1 L g + L f ( v o q v g q )
f D Q , G = T g l o b a l l o c a l f d q , g
T g l o b a l l o c a l = [ cos ( δ ) sin ( δ ) sin ( δ ) cos ( δ ) ] , δ = ( ω ω P L L )
v g d = V g D cos ( δ 1 ) V g Q sin ( δ 1 ) v g q = V g D sin ( δ 1 ) + V g Q cos ( δ 1 )
where f D Q ,   G is the global frequency, f d q , g is the local frequency, and V g , D Q is the grid bus voltage in the global reference frame.
[ Δ i ˙ l d q Δ v ˙ o d q Δ i ˙ o d q ] = A L C L [ Δ i l d q Δ v o d q Δ i o d q ] + B L C L 1 [ Δ v i d Δ v i q ] + B L C L 2 [ Δ δ 1 ] + B L C L 3 [ Δ ω P L L ]
where
A L C L = [ R c L c ω P L L 1 L c 0 0 0 ω P L L R c L c 0 1 L c 0 0 1 C f 0 0 ω P L L 1 C f 0 0 1 C f ω P L L 0 0 1 C f 0 0 1 L g + L f 0 ( R g + R f ) L g + L f ω P L L 0 0 0 1 L g + L f ω P L L ( R g + R f ) L g + L f ]
B L C L 1 = [ 1 L c 0 0 1 L c 0 0 0 0 0 0 0 0 ] ,   B L C L 2 = [ 0 0 0 0 V g D sin ( δ 1 ) L g + L f V g D cos ( δ 1 ) L g + L f ] ,   B L C L 3 = [ i l q i l d v o q v o d i o q i o d ]

3.5. Bus

The grid angular frequency is represented by ω g , and the inverter angular frequency is ω P L L . The phase difference between the grid and the inverter is δ 1 , which is used to transform the grid voltage from the global frame to the local frame. This angle varies consistently with the variation of the inverter frequency, as described in (28) and (29). The small signal linearized state space model of the bus is represented in (30).
δ ˙ 1 = ω g ω P L L
δ ˙ 1 = ω g k p , P L L v o q k i , P L L Φ P L L
[ Δ δ ˙ 1 ] = A B U S [ Δ δ 1 ] + B B U S [ Δ ω c o m ]
where
A B U S = [ 0 ] ,   B B U S = [ 1 ]

3.6. Full System

The combined submodules are represented in (31) where there are 13 states and five inputs. The A s y s represents the state matrix of the full model, and B s y s represents the input matrix to the system.
[ Δ x ˙ s y s ] = A s y s [ Δ x s y s ] + B s y s [ Δ u s y s ]
u = [ ω P L L P * Q * V g D V g Q ] T x = [ P L L Φ P L L P Q γ d γ q i l d i l q v o d v o q i o d i o q δ 1 ] T
A s y s = A P L L 0 0 B P L L 0 B D 2 C P L L A D 0 B D 1 + B D 2 D P L L 0 B C 1 D D 2 C P L L B C 1 C D A C B C 2 + B C 1 D D 1 + B C 1 D D 2 D P L L 0 B L C L 3 C P L L + B L C L 1 D C 1 D D 2 C P L L B L C L 1 D C 1 C D B L C L 1 C C A L C L + B L C L 1 D C 2 + B L C L 1 D C 1 D D 1 + B L C L 1 D C 1 D D 2 D P L L + B L C L 3 D P L L B L C L 2 B B U S C P L L 0 0 B B U S D P L L A B U S 13 × 13
A s y s = 0 K i , P L L 0 0 0 0 0 0 0 K p , P L L 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 K i , P L L · m 0 0 0 0 0 0 0 K i , P L L · m 0 0 0 0 0 0 0 0 0 0 0 n 0 0 0 0 0 2 · K i , P L L · m 3 · v o d 0 0 0 0 1 0 2 · P * m ω P L L ω * 3 · v o d 2 2 · K i , P L L · m 3 · v o d 0 0 0 0 0 0 0 0 0 0 1 2 · Q * + E * · n 3 · v o d 2 0 0 0 0 0 K i , p L L · i l q 2 · K i , p L L · K p c · m 3 · L c · v o d 0 0 K i c L c 0 K p c L c R c L c ω P L L ω n 1 L c 2 · K p c · P * m ω P L L ω * 3 · L c · v o d 2 K p , p L L · i l q 2 · K p c · K p , p L L · m 3 · L c · v o d 0 0 0 0 K i , P L L · i l d 0 0 0 K i c L c ω n ω P L L K p c L c R c L c 2 · K p c · Q * + E * · n 3 · L c · v o d 2 K p , P L L · i l d 1 L c 0 0 0 0 K i , P L L · v o q 0 0 0 0 1 C f 0 0 ω p L L + ( K p , L L · v o q ) 1 C f 0 0 0 K i , P L L · v o d 0 0 0 0 0 1 C f ω P L L K p , P L L · v o d 0 1 C f 0 0 K i , P L L · i o q 0 0 0 0 0 0 1 L g + L f K p , P L L · i o q R g + R f L g + L f ω P L L V g D s i n δ 1 L g + L f 0 K i , P L L · i o d 0 0 0 0 0 0 0 1 L g + L f ( K p , P L L · i o d ) ω P L L R g + R f L g + L f V g D c o s δ 1 L g + L f 0 K i , P L L 0 0 0 0 0 0 0 K p , P L L 0 0 0 13 × 13

4. Simulation Model

The MG system was modelled in MATLAB Simulink, as shown in Figure 5, and the model parameters are reported in Table 1. The solver type in MATLAB was set to ode45x with a time step of 1 µs. The total simulation time was eighteen minutes and it was performed using a personal computer with an Intel® Core™ i7-8750H CPU @ 2.2 GHz processor and memory (RAM) of 16.00 GB.

5. Results and Discussion

Figure 6 shows the active and reactive power output of the MG where they match the reference values of 500 W and 0 VAR, respectively. Small perturbations in the parameter values of any system may change the system’s behavior significantly, causing instability. The effect of parameter variation on a specific system behavior can be determined by calculating the parameter sensitivity. The parameter sensitivity analysis gives information about parameters that can cause instability with only small changes in their values. Therefore, parameter values can be designed carefully to avoid instability. The parameter sensitivity is calculated by
S e n p j λ i = | λ i ( p + Δ p j ) λ i ( p ) | Δ p j
where λ i ( p ) is the ith eigenvalue of the system before perturbation, λ i ( p + Δ p j ) is the ith eigenvalue of the system after perturbation, and Δ p j is the difference between the original parameter value and the perturbed parameter value. The parameter sensitivity analysis is illustrated in Figure 7, which shows that the most sensitive parameters in the MG system were Lg and Lf. Therefore, it was decided to study the impact of changing their values on the stability condition.
The eigenvalues of the combined small signal model can be related to the system parameters and states. This relationship can help us to optimize the system dynamic response and stability. The trace of the eigenvalues when changing the grid side inductance of the LCL filter is shown in Figure 8. The inductance value was varied in steps of 0.3 mH in the time domain simulation performed by MATLAB. The figure shows that the value of Lf that destabilized the system was 6.2 mH. In order to assess the accuracy of the mathematical model, it was decided to simulate the system with different values of Lf and measure the d-component of the voltage at the PCC. It was found in the simulation that the system became unstable at 6.2 mH of Lf, as shown in Figure 9. The figure also shows that when reducing the Lf value by 0.3 mH, the system was stable, while the system was unstable when it was increased by 0.3 mH.
When varying the transmission line inductance in the mathematical model in steps of 0.3 mH, it was found that the system became unstable at 6.1 mH, as shown in Figure 10. This value is validated by the simulation results shown in Figure 11. The figure shows that when applying 6.1 mH and above, the system is unstable, with the appearance of high voltage values with oscillations. When reducing the inductance value by 0.3 mH and when applying the original setting of 3.44 µH, the system was stable, and the voltage stayed around 120 V, without any high oscillations.
The inverter-side inductance of the LCL filter did not affect the stability of the system, as shown in Figure 12. This mathematical result is validated by the simulation results shown in Figure 13. The figure shows that the system was stable, despite the significant increase in the inductance value. These results are confirmed by the parameter sensitivity analysis shown in Figure 7, which shows that the MG system was not highly sensitive to changes in the inverter-side inductance.
The parameter sensitivity analysis showed that the poorly damped modes were related to the LCL filter parameters and the feeder inductance. Consequently, the damping ratio of these modes can be increased by increasing the series damping resistances of these sensitive parameters. However, passive damping has some drawbacks, such as increasing the power loss in the filter and weakening the ability of the LCL filter to attenuate the harmonics generated by the switching in the inverter. Furthermore, it is not possible to increase the damping of the feeder because, in practice, this represents the distance between the MG and grid. Therefore, active damping should be designed to avoid decreasing the system’s efficiency.

6. Conclusions

SSS has been discussed in this article for a MG in grid-connected mode, to improve the control design and also the stability of the system. A controller proposal, to control a MG in grid-connected mode, was introduced in this paper. The whole system of the MG, including the control system, was divided into submodules, and the small signal model for each submodule was produced in detail. The small signal models of the submodules were combined together on the inverter reference frame, to generate the A matrix of the full system, which was used to analyze the stability of the full system. A sensitivity analysis of the MG parameters was performed, before applying small variations to the most sensitive parameters, to the test stability. The aim of this dynamic modelling was to design the MG parameters that ensure stability. It was found that increasing the inductance destabilizes the MG system. When increasing the grid side inductance of the LCL filter beyond the value of 6.2 mH, the system became unstable; and when increasing the transmission line inductance beyond the value of 6.1 mH, the system also became unstable. The small signal model mathematical results were validated by comparison with the simulation results produced using MATLAB Simulink. Identifying stability margins is a crucial aspect of increasing grid resilience and of avoiding instability, which causes power outages and high costs, resulting from damages to the grid infrastructure. Furthermore, power outages may affect the economy indirectly by reducing sales and productivity. In future work, a solar power plant will be added, alongside electrical loads that represent the real loads of consumers, before studying the SSS, based on variations in the power source and electrical loads.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author would like to acknowledge Jouf University for supporting this research project.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Wang, Y.; Phung, T.; Ravishankar, J. Impact of battery storage on micro-grid transient performance. In Proceedings of the 2014 IEEE Region 10 Conference (TENCON 2014), Bangkok, Thailand, 22–25 October 2014. [Google Scholar]
  2. Tran, Q.; Sanseverino, E.; Zizzo, G.; Silvestre, M.; Nguyen, T.; Tran, Q. Real-time minimization power losses by driven primary regulation in islanded microgrids. Energies 2020, 13, 451. [Google Scholar] [CrossRef] [Green Version]
  3. Majumder, R. Some aspects of stability in Microgrids. IEEE Trans. Power Syst. 2013, 28, 3243–3252. [Google Scholar] [CrossRef]
  4. Zeng, Z.; Yang, H.; Zhao, R. Study on small signal stability of microgrids: A review and a new approach. Renew. Sust. Energy Rev. 2011, 15, 4818–4828. [Google Scholar] [CrossRef]
  5. Rasheduzzaman, M. Small Signal Modeling and Analysis of Microgrid Systems. Ph.D. Dissertation, Department of Electrical and Computer Engineering, Missouri University of Science and Technology, Rolla, MO, USA, 2015. [Google Scholar]
  6. Gibbard, M.J.; Pourbeik, P.; Vowles, D.J. Small-Signal Stability, Control and Dynamic Performance of Power Systems; University of Adelaide Press: Adelaide, Australia, 2015; ISBN 978-1-925261-03-5. [Google Scholar]
  7. Wang, Y.; Lu, Z.; Min, Y.; Wang, Z. Small signal analysis of microgrid with multiple micro sources based on reduced order model in islanding operation. In Proceedings of the 2011 IEEE Power and Energy Society General Meeting, Detroit, MI, USA, 24–28 July 2011. [Google Scholar]
  8. Abbasi, A.K.; Mustafa, M.W. Mathematical model and stability analysis of inverter-based distributed generator. HPC 2013, 2013, 195038. [Google Scholar] [CrossRef]
  9. Zhang, Y.; Jiang, Z.; Yu, X. Small-signal modeling and analysis of parallel-connected voltage source inverters. In Proceedings of the 2009 IEEE 6th International Power Electronics and Motion Control Conference, Wuhan, China, 17–20 May 2009. [Google Scholar]
  10. Rasheduzzaman, M.; Mueller, J.; Kimball, J. Small-signal modeling of a three-phase isolated inverter with both voltage and frequency droop control. In Proceedings of the 2014 IEEE Applied Power Electronics Conference and Exposition, Fort Worth, TX, USA, 16–20 March 2014. [Google Scholar]
  11. Adelpour, M.; Hamzeh, M.; Sheshyekani, K. Comprehensive small-signal stability analysis of islanded synchronous generator-based microgrids. SEGN 2021, 26, 100444. [Google Scholar] [CrossRef]
  12. Paspatis, A.G.; Konstantopoulos, G.C.; Dedeoglu, S. Control design and small-signal stability analysis of inverter-based microgrids with inherent current limitation under extreme load conditions. Electr. Power Syst. Res. 2021, 193, 106929. [Google Scholar] [CrossRef]
  13. Qiuye, S.; Bingyu, W.; Xiaomeng, F.; Shiyan, H. Small-signal stability and robustness analysis for microgrids under time-constrained dos attacks and a mitigation adaptive secondary control method. Sci. China Inf. Sci. 2022, 65, 162202. [Google Scholar] [CrossRef]
  14. Mohamed, Y.; El-Saadany, E. Adaptive decentralized droop controller to preserve power sharing stability of parallel inverters in distributed generation microgrids. IEEE Trans. Power Electron. 2008, 23, 2806–2816. [Google Scholar] [CrossRef]
  15. Raju, E.S.N.; Jain, T. Small signal modelling and stability analysis of AC islanded ac microgrid with inverter interfaced DGs. In Proceedings of the 2014 International Conference on Smart Electric Grid (ISEG), Guntur, India, 19–20 September 2014. [Google Scholar]
  16. Hadisupadmo, S.; Hadiputro, A.; Widyotriatmo, A. A small signal state space model of inverter-based microgrid control on single phase AC power network. Internetworking Indones. J. 2016, 8, 71–76. [Google Scholar]
  17. Katiraei, F.; Iravani, M.R.; Lehn, P.W. Small-signal dynamic model of a micro-grid including conventional and electronically interfaced distributed resources. IET Gener. Transm. Distrib. 2007, 1, 369–378. [Google Scholar] [CrossRef] [Green Version]
  18. Deepak, P.; Ramtin, M.; Tuncay, A.; Ali, D. Small signal stability constrained optimal power flow for inverter dominant autonomous microgrids. IEEE Trans. Ind. Electron. 2021, 69, 7318–7328. [Google Scholar] [CrossRef]
  19. Aram, M.; Ali, G.; Mohammad, B. Microgrid small-signal stability analysis considering dynamic load model. IET Renew. Power Gener. 2021, 15, 2799–2813. [Google Scholar] [CrossRef]
  20. Han, Y.; Qian, Z.; Shao, D.; Lin, Q.; Yin, S.; Chen, M.; Ai, Q. Small signal stability analysis of microgrid with multiple parallel inverters. In Proceedings of the IOP Conference Series: Earth and Environmental Science, Zhuhai, China, 15–17 January 2021. [Google Scholar]
  21. Chen, X.; Zhang, C.; Huang, Q.; Ofori-Oduro, M. Small-signal modeling and analysis of grid-connected inverter with power differential droop control. HPC 2016, 2016, 3965945. [Google Scholar] [CrossRef]
  22. Zhang, Y.; Liu, H.; Song, P.; Wu, L. Small-signal stability analysis of multiple grid-connected virtual synchronous generators. In Proceedings of the 2018 China International Conference on Electricity Distribution (CICED), Tianjin, China, 17–19 September 2018. [Google Scholar]
  23. Huang, L.; Wu, C.; Zhou, D.; Blaabjerg, F. Comparison of three small-signal stability analysis methods for grid-following inverter. In Proceedings of the 2021 International Aegean Conference on Electrical Machines and Power Electronics (ACEMP) & 2021 International Conference on Optimization of Electrical and Electronic Equipment (OPTIM), Brasov, Romania, 2–3 September 2021. [Google Scholar]
  24. Wenjuan, D.; Wenkai, D.; Haifeng, W. Small-signa stability limit of a grid-connected pmsg wind farm dominated by the dynamics of plls. IEEE Trans. Power Syst. 2020, 35, 2093–2107. [Google Scholar] [CrossRef]
  25. Uros, M.; Ognjen, S.; Petros, A.; Evangelos, V.; Duncan, C.; Gabriela, H. Understanding small-signal stability of low-inertia systems. IEEE Trans. Power Syst. 2021, 36, 3997–4017. [Google Scholar] [CrossRef]
  26. Reznik, A.; Godoy, M.; Al-Durra, A.; Muyeen, S. LCL filter design and performance analysis for grid-interconnected systems. IEEE Trans. Indust. App. 2014, 50, 1225–1232. [Google Scholar] [CrossRef]
Figure 1. Three-phase voltage source inverter connected to a grid via a control method and a LCL filter.
Figure 1. Three-phase voltage source inverter connected to a grid via a control method and a LCL filter.
Sustainability 14 09372 g001
Figure 2. PLL.
Figure 2. PLL.
Sustainability 14 09372 g002
Figure 3. Droop controller.
Figure 3. Droop controller.
Sustainability 14 09372 g003
Figure 4. Current controller.
Figure 4. Current controller.
Sustainability 14 09372 g004
Figure 5. Screen capture of a MG system in grid-connected mode modelled in MATLAB Simulink.
Figure 5. Screen capture of a MG system in grid-connected mode modelled in MATLAB Simulink.
Sustainability 14 09372 g005
Figure 6. Output power of the MG.
Figure 6. Output power of the MG.
Sustainability 14 09372 g006
Figure 7. Parameter sensitivity analysis of the MG system.
Figure 7. Parameter sensitivity analysis of the MG system.
Sustainability 14 09372 g007
Figure 8. Trace of the eigenvalues of the linearized mathematical model under the variation of Lf in the range of 0.16 mH and 0.01 H.
Figure 8. Trace of the eigenvalues of the linearized mathematical model under the variation of Lf in the range of 0.16 mH and 0.01 H.
Sustainability 14 09372 g008
Figure 9. Simulation results of the voltage at the PCC when varying Lf.
Figure 9. Simulation results of the voltage at the PCC when varying Lf.
Sustainability 14 09372 g009
Figure 10. Trace of the eigenvalues of the linearized mathematical model under the variation of Lg in the range of 3.44 µH and 0.01 H.
Figure 10. Trace of the eigenvalues of the linearized mathematical model under the variation of Lg in the range of 3.44 µH and 0.01 H.
Sustainability 14 09372 g010
Figure 11. Simulation results of the voltage at the PCC when varying Lg.
Figure 11. Simulation results of the voltage at the PCC when varying Lg.
Sustainability 14 09372 g011
Figure 12. Trace of the eigenvalues of the linearized mathematical model under the variation of Lc in the range of 4.77 mH and 0.01 H.
Figure 12. Trace of the eigenvalues of the linearized mathematical model under the variation of Lc in the range of 4.77 mH and 0.01 H.
Sustainability 14 09372 g012
Figure 13. Simulation results of the voltage at the PCC when varying Lc.
Figure 13. Simulation results of the voltage at the PCC when varying Lc.
Sustainability 14 09372 g013
Table 1. Parameters of the grid-tied inverter system [26].
Table 1. Parameters of the grid-tied inverter system [26].
ParameterSymbolValue
Inverter-side parasitic resistance R c 0.25   Ω
Grid-side parasitic resistance of the LCL filter R f 0.25   Ω
Transmission line parasitic resistance R g 0.25   Ω
Inverter-side inductance of the LCL filter L c 4.77   mH
Grid-side inductance of the LCL filter L f 0.16   mH
Transmission line inductance L g 3.44   μ H
Filter capacitance C f 14.7   μ F
Proportional gain of the current controller k p c 14.24
Integral gain of the current controller k i c 2.11 × 10 4
Frequency droop gainm 1 × 10 6
Voltage droop gain n 1 × 10 6
Inverter input DC voltage v d c 180   V
Switching frequency f s 8   kHz
Grid voltage magnitude v g 120   V
Fundamental frequency f n 50   Hz
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alnuman, H. Small Signal Stability Analysis of a Microgrid in Grid-Connected Mode. Sustainability 2022, 14, 9372. https://doi.org/10.3390/su14159372

AMA Style

Alnuman H. Small Signal Stability Analysis of a Microgrid in Grid-Connected Mode. Sustainability. 2022; 14(15):9372. https://doi.org/10.3390/su14159372

Chicago/Turabian Style

Alnuman, Hammad. 2022. "Small Signal Stability Analysis of a Microgrid in Grid-Connected Mode" Sustainability 14, no. 15: 9372. https://doi.org/10.3390/su14159372

APA Style

Alnuman, H. (2022). Small Signal Stability Analysis of a Microgrid in Grid-Connected Mode. Sustainability, 14(15), 9372. https://doi.org/10.3390/su14159372

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