Small-Signal Modeling and Analysis of Grid-Connected Inverter with Power Differential Droop Control

The conventional voltage and frequency droop control strategy in grid-connected inverter suffers a major setback in the presence of disturbance by producing oscillations. Adding a power differential term in droop controller is an effective way to address such drawback. In this paper, grid-connected inverter’s small-signal models of the conventional droop control and the power differential droop control are established. The eigenvalues of the models are then determined by system matrix. The eigenvalues analysis is presented which helps in identifying the relationship between the system stability and controller parameters. It is concluded that the damping ratio of dominant low-frequency eigenvalues increased and the oscillation caused by the disturbance is suppressed when a power differential term is added to the droop control method. The MATLAB/Simulink models of grid-connected inverter with both control strategies are also established to validate the results of small-signal analysis.


Introduction
The concept of distributed power generation and distributed energy storage has come to gain more attention because of the several advantages it offers over centralized power generation and storage [1][2][3][4].Power electronic converters, typically the inverters, are used to interface these distributed generators (DGs) to the grid [5,6].They make the sources more flexible to control and operate.However, their physical negligible inertia also makes the system prone to oscillations when subjected to disturbances.
DGs in a microgrid are expected to operate in an autonomous manner.And their operation and dynamics are strongly contingent upon the connected sources and on the power regulation control of the power electronic converters.The conventional means of control in DG is carried out by mimicking conventional generators.Based on this, the droop control is presented in [7][8][9][10][11].A droop control scheme uses local power to detect changes in the system and adjust the operating points of generators accordingly.The scheme achieves accurate power sharing while maintaining regulation of the DG's voltage magnitude and frequency.However, the stability and performance of a control scheme are contingent on system parameters, their state, and the droop gain.Analysis of previous works has shown that the selection of the droop gain is crucial for the stability and performance of the DG.Empirical formulas are used in the selection of these gains [12].However, it is not sufficient to analyze the consequences of the droop gains on the stability of the DG and microgrid as a whole.To address this problem, small-signal modeling is used [11][12][13][14][15].
On the other hand, the droop controller cannot provide sufficient low-frequency damping to compensate for the lack of the physical inertia in DG [14,15].It still makes the system potentially susceptible to oscillation by disturbances.In conventional synchronous generators, power system stabilizer (PSS) is used to rectify the same oscillation problems in the output power.The PSS improves the damping of the electromechanical oscillations by acting through the generator excitation system to produce electrical torque in addition to the damping torque according to the speed deviation generated [16].Using the same idea, it can be stipulated that adding a power differential term to the droop control is a possible means of solving the oscillation problem associated with the inverter in DGs.And it has been found to be effective in [15].However, the stability and performance of the DG have not yet been fully analyzed and compared with The circuit of inverter and the controller.
the droop control after the addition of the power differential term.
In this paper, complete small-signal models for DG are established for both the droop control and the power differential droop control.Further, the eigenvalue distribution of each control strategy is plotted to compare the damping ratios.Stability analysis of the DG under both control strategies is conducted in order to know the relation between stability and control parameters.Finally, MATLAB/Simulink simulations are used to validate the small-signal model and analysis.
The rest of the paper is in the following order.In Section 2, the structure of DG is introduced.In Section 3, a complete small-signal model for DG with the droop control is presented.Based on this work, in Section 4, modeling of DG with power differential droop control is presented.Stability analysis and simulation results are provided in Section 5.And conclusions are drawn in Section 6.

Grid-Connected Inverter
DG commonly interfaces with microgrid through voltage source inverter (VSI).As shown in Figure 1, VSI is interfaced with the microgrid at bus V  .Modeling of the system is in two main sections: the VSI and the bus.The VSI is divided into two main parts: the power block and the control block.The three-phase inverter bridge, output LCL filter, and transmission line make up the power block, while the control block consists of the power controller, voltage controller, and current controller.The switching effect of the inverter can be ignored for cases of high frequency.Assuming an ideal voltage source for the inverter, the DC bus dynamics are also ignored.With the exception of the ignored components mentioned, state models are presented for all the above components in the next section.

Model with Droop Control
3.1.Power Controller. Figure 2 shows the block diagram of power controller, where V o , V o ,  o , and  o are --axis component of inverter output voltage V o and current  o .Here,  is calculated as three-phase active power and  as singlephase reactive power.These power components are passed through low-pass filters, shown as (1), to obtain the active power  and reactive power . c represents the cut-off frequency of low-pass filters.Hence, To mimic the inertia characteristics of conventional generators and provide a degree of negative feedback, an artificial droop characteristic is introduced in the inverter frequency  as in (2). is set according to the droop gains  and .And phase  is set by integrating the frequency as in (3), where  n is rated frequency. rate is the rated active power.Hence, Similarly, a droop characteristic is also introduced in the voltage magnitude as given in (4).Here, the voltages are under the - synchronous rotating frame with reference to the inverter's frequency.Therefore, the output voltage magnitude reference is represented at the -axis, and the -axis reference is set to zero.Therefore, Equations ( 1) to ( 4) are then rearranged and linearized to obtain the state-space form as shown in All matrices in (5) have been shown in the Appendix, the same as below.

Voltage and Current Controller.
The controllers for both the voltage and the current are shown in Figure 3.Both controllers employ standard proportional-integral (PI) regulators with decoupling and feedforward control loops.The voltage controller generates the reference vector for the current controller.Set outputs of integrator of PI regulators are   ,   (voltage controller) and   ,   (current controller).Then, the dynamics of the voltage and current controller can be given as where  pv and  iv are the proportional and integral gains of voltage controller. is the feedforward gain. pc and  ic are the proportional and integral gains of current controller.
The aforesaid equations can be linearized and the following state-space model of voltage and current controller unit can be obtained: ] . (7)

Output LCL Filter and Line.
Assuming that the inverter bridge losses are negligible, then The output filter and line model can be represented by the following equations: where  f and  c are equivalent resistances of  f and  c .And  g and  g are equivalent resistance and inductance of line.
It can be seen that bus voltage V b is the input of this model, and there is a phase difference ( 2 ) between V b and inverter's output voltage V o .Hence, V b needs to be converted to the inverter's reference frame.The equations for transformation are shown as where  bus is the amplitude of V b .Linearizing the above equations, the state-space model of the output filter and line is given by 3.4.Bus Model.Assume that the amplitude and frequency of the bus are constant.Then,  2 changes according to the inverter's frequency.For simplicity, the inverter's reference frame is set as common reference frame.The dynamics of  2 can then be given as And the state-space model of the bus can be obtained as

Complete Model of the Inverter.
A complete model of the grid-connected inverter is achieved by combining all submodels (5), ( 7), (10), and ( 12), and shown in

Model with Power Differential Droop Control
To improve inverter damping of low-frequency oscillation modes, by reference to PSS's operational principle, the droop characteristics are modified as follows [15]: Since modification only affects the power controller (5) and bus (12) models, the previous power controller models are updated as in ( 15) so as to construct a new system model with the above modification.Hence, The bus model with power differential term is modified as follows: The system model with power differential term droop control is represented as

Stability Analysis and Simulation
To evaluate and compare the performances of both control schemes, according to Figure 1, simulation models have been built via MATLAB/Simulink.System parameters are given in Table 1, and initial conditions of the system are given in Table 2.These steady-state operating point conditions were  obtained from a MATLAB/Simulink time-step simulation of the system.According to the established small-signal models, all eigenvalues calculations are based on Table 2.The eigenvalues of high and intermediate frequency modes have weak influence on the stability of the system [14].Therefore, in this paper, point of focus of analysis is on the low-frequency modes.

Active Power Droop Gain 𝑚.
Figure 4 shows the dominant low-frequency eigenvalues as a function of the active power droop gain  (1 − 5 ≤  ≤ 1 − 3), the reactive power droop gain  = 5 − 5, and the power differential gain  d =  d = 8−6.It shows that as  is increased, the complexconjugate pair of the dominant eigenvalues of the two control strategies move towards unstable region making the system more oscillatory.When  = 8−4, the dominant eigenvalues of the droop control enter the right half plane, eventually leading to instability.Contrary to that, the eigenvalues of the power differential droop control stay on the left half plane.These trajectories show that, under the same condition, the range of  increases when the power differential term is added to the conventional droop control.Figure 4 also shows that the complex-conjugate eigenvalues of the power differential droop control have a larger absolute real value than the droop control, which means that power differential term increases system damping and contributes to system stability.
In order to verify the correctness of the above smallsignal stability analysis, an effective method that selects some parameters which lead to system instability or oscillation is carried out.Then, observe whether the response of simulation model is consistent with the small-signal analysis.
Figure 5 shows output power response for 20% step change in load, when  = 8−4 for the two control strategies.It can be seen that, after the disturbance, the output power with the droop control starts to oscillate and eventually results in instability, while output power with the power differential droop control is steady.
Figure 6 shows the same 20% step change in load when  = 4−4 for the two control strategies.From the calculation result of eigenvalues, when  = 4−4, the complex-conjugate eigenvalues of two control strategies are −6.9 ± 52.2 (droop control) and −27.7 ± 47.4 (power differential droop control) with damping ratios 0.13 (droop control) and 0.5 (power differential droop control).From Figure 6, it can be seen that the active power of the power differential droop control has less oscillation than the droop control.In Figure 6(a), it can be seen that the oscillation period is 0.12 s, which is the same oscillation period given by complex-conjugate eigenvalues.These results validate the small-signal model and the analysis.

Reactive Power Gain 𝑛.
Figure 7 shows the dominant low-frequency eigenvalues as a function of the reactive power droop gain  (1 − 5 ≤  ≤ 7.5 − 4), the active power droop gain  = 8 − 5, and the power differential gain  d =  d = 8 − 6.It can be seen that as  is increased, the complex-conjugate dominant eigenvalues of the two control strategies move away from the imaginary axis, and the real eigenvalues move towards the right half plane, eventually leading the system to instability.It can also be seen that the power differential droop control has larger damping than the droop control at low-frequency oscillation modes.
Figure 8 shows the output power response for 20% step change in load when  = 5 − 4 for the two control strategies.Again, from the calculation result of eigenvalues, when  = 5 − 4, the eigenvalues of the two control strategies are both in the left half plane, so both systems are stable.The complexconjugate eigenvalues of the droop control are −25.4 ± 24.7; its given oscillation period of 0.25 s and damping ratio of 0.72 correspond to rapid decay of oscillation.The complexconjugate eigenvalues of the power differential droop control are −124 ± 54, and the real eigenvalue is −16.7,so the real eigenvalues play a dominant role; there is no low-frequency oscillation mode.It can be observed that the results in Figure 8 closely match the above analysis.
The simulation results in Figures 5, 6, and 8 are consistent with the stability and dynamic characteristics given by the small-signal model, demonstrating the validity of the smallsignal modeling and analysis.Considering the need for stability and suitable damping, this paper takes  = 8 − 5 and  = 5 − 5.    9 shows the dominant low-frequency eigenvalues as a function of the power differential gains

Analysis for Power Differential Gains
It can be seen that increasing  d ,  d causes the complexconjugate dominant eigenvalues of the power differential droop to move in a spiral manner.At the initial stage of the increment, the eigenvalues move away from the right half plane while approaching the real axis.However, further increment causes the eigenvalues to move away from the real axis and close to the right half plane, eventually moving into the right half plane.Therefore, selecting small values for  d ,  d causes the system to have small damping.Selecting large values also causes the system to have more oscillations in case of a disturbance, eventually leading to instability.
Figure 10 shows the output power response for the power differential droop control for a 20% step in load when  d ,  d take values of 1 − 7, 4 − 6, and 2.75 − 5 under initial conditions.It can be seen that the output power in Figure 10(a) has more oscillations than that in Figure 10(b), which is because when  d ,  d = 1 − 7, the system damping is smaller than when  d ,  d = 4 − 6.In Figure 10(c), it can be seen that the system becomes unstable after the load step.This however is expected because the value of  d ,  d in this case moves the dominant complex-conjugate eigenvalues into the right half plane.Obviously, 4 − 6 is an appropriate value of  d and  d , when  = 8 − 5 and  = 5 − 5.

Conclusion
The grid-connected inverter with the droop control is potentially susceptible to oscillation due to its poor damping of lowfrequency mode.A solution is adding power differential term in droop controller.In this paper, the small-signal models of grid-connected inverter with both the droop control and the power differential droop control are established.The eigenvalues of these models were analyzed.And the following conclusions are obtained: the damping increased after adding power differential term, and the system oscillation is suppressed.The range of active power droop gain  is increased compared with the droop control under stable conditions.The system with the power differential droop control has poor damping both on the small and on the large power differential gains  d ,  d .And the system will be unstable with larger  d ,  d .These results were verified by simulation model via MATLAB/Simulink.It was observed that the simulation system dynamics closely match the conclusion of analysis.

A. Matrices in Section 3
A.2. Voltage Controller.One has the following: Mathematical Problems in Engineering And, for current controller, A.3.Output LCL Filter and Line.One has the following: ] ,

B. Matrices in Section 4
One has the following:

Figure 4 :
Figure 4: Dominant locus of eigenvalue with  increasing.

Figure 7 :
Figure 7: Dominant locus of eigenvalue with  increasing.

DD+ B V1 D Pk dt 0 BB
3 d  c  o −3 d  c  o −3 d  c  o −3 d  c  o 0 0 3 c  o 3 c  o 3 c  o 3 c  o 0 0  c  o − c  o − c  o  c  o P dt = [0 0 −3 d  c  o −3 d  c  o −3 d  c  o −3 d  c  o ] , Pk dt = [ 0 0 − d  c  o  d  c  o  d  c  o − d  c  o C1 D V1 C Pk dt B C1 C V 0 B C1 D V1 D Pk dt + B C1 D V2 + B C2 0 LCL1 D C1 D V1 C Pk dt + B LCL3 C P dt B LCL1 D C1 C V B LCL1 C C A LCL + B LCL1 (D C1 D V2 + D C2 ) +B LCL1 D C1 D V1 D Pk dt + B LCL3 D P dt B LCL2