DOI QR코드

DOI QR Code

Performance Improvement Strategy for Parallel-operated Virtual Synchronous Generators in Microgrids

  • Zhang, Hui (College of Automation and Information Engineering, Xi'an University of Technology) ;
  • Zhang, Ruixue (College of Automation and Information Engineering, Xi'an University of Technology) ;
  • Sun, Kai (State Key Lab. of Power Systems, Department of Electrical Engineering, Tsinghua University) ;
  • Feng, Wei (State Key Lab. of Power Systems, Department of Electrical Engineering, Tsinghua University)
  • Received : 2018.07.07
  • Accepted : 2018.12.03
  • Published : 2019.03.20

Abstract

The concept of virtual synchronous generators (VSGs) is a valuable means for improving the frequency stability of microgrids (MGs). However, a great virtual inertia in a VSG's controller may cause power oscillation, thereby deteriorating system stability. In this study, a small-signal model of an MG with two paralleled VSGs is established, and a control strategy for maintaining a constant inertial time with an increasing active-frequency droop coefficient (m) is proposed on the basis of a root locus analysis. The power oscillation is suppressed by adjusting virtual synchronous reactance, damping coefficient, and load frequency coefficient under the same inertial time constant. In addition, the dynamic load distribution is sensitive to the controller parameters, especially under the parallel operation of VSGs with different capacities. Therefore, an active power increment method is introduced to improve the precision of active power sharing in dynamic response. Simulation and experimental is used to verify the theoretical analysis findings.

Keywords

I. INTRODUCTION

Recently, the solar power as renewable energy source, have been connected with utility grids by using a current-controlled inverter, which result in an increasingly serious effect on utility grids. Consequently, researchers [1]-[6] have proposed a virtual synchronous generator (VSG) technology, which make the inverters have the same characteristics of synchronous generator. The VSG is a new and interesting means of large-scale distributed generation (DG) in connecting with utility grids.

A conventional droop controller, which simulates the droop characteristic of a synchronous generator, is often adopted by inverters to share the active and reactive power in islanded mode [7-10]. However, the system’s frequency stability cannot be guaranteed because the equivalent inertia is often small. Thus, a VSG control strategy for islanding mode is presented by introducing the rotor motion equation of a synchronous generator, thereby enhancing the frequency stability of a common alternating current (AC) bus. An improved control strategy of synchronverters is proposed [11], [12]. The inverter and synchronous generator are mathematically equivalent, but no current loop exists in the controller of synchronverters. Thus, the power quality of a grid-connected current is easily affected. The current loop is introduced, and the control strategy for parallel-operated VSGs is investigated to enhance the power quality [13]. However, the design methodology for parallel-operated VSGs with different capacities is not considered. A synchronous generator emulation control strategy for a voltage source converter (VSC)-HVDC station [14] is designed, in which the desired frequency support to a low-inertia grid is provided by the VSC station with the proposed control strategy. The system stability cannot be guaranteed when numerous, varying loads are connected to a weak power system. Thus, the effects of suppressing frequency fluctuations by using a VSG have been investigated [15], [16]. However, few studies have focused on stability analysis. A VSG with adaptive inertia is considered [17], [18] such that great inertia is used during acceleration, and small inertia is utilized to increase the deceleration effect. The introduction of a VSG’s virtual inertia improves the frequency stability of an AC bus. However, power oscillation is not considered when a VSG is connected to the common AC bus [19]. A small-signal model for a VSG is established [20]-[22], and the stability in grid-connected and islanded modes is analyzed. However, the parallel-operated mode has not been evaluated.

A great virtual inertia effectively suppresses AC bus frequency fluctuation under frequent load disturbance and increases power oscillation when the VSG is connected to the common AC bus. Hence, a conventional droop control method is presented [23] to suppress power oscillation by introducing differential negative feedback. However, this method seriously affects the stability of a common AC bus. A power system stabilizer (PSS) is adopted [24], [25] to enhance the transient stability of a standalone microgrid (MG) system, which can also suppress power oscillation. The introduction of PSS parameters adds complexity to the design of controller parameters. Virtual admittance is adopted [26] to reduce the control runaway capability when the measured current is affected by harmonic, transient, and measured noises, which suppresses the power oscillation. However, the selection of virtual admittance parameter is not designed. Therefore, a control strategy is investigated to suppress power oscillation without influencing the system’s frequency stability. Moreover, the VSG’s output frequency slowly fluctuates with the occurrence of load disturbances, which is different from conventional droop control, and the dynamic response of the active power allocation of paralleled VSGs is inaccurate when the frequency variations of two VSGs are different. The power oscillation is suppressed by optimizing the parameters of the virtual impedance of VSG in parallel-operated VSGs with same capacities [29]. But the parallel operation of VSG with different capacities will increase power oscillation and return circuit current. Therefore, other parameters J, D,and Kω, of the parallel VSG with same capacity are analyzed, and two methods to suppress power oscillation are proposed in this paper. In addition, the parameter design method for parallel-operated VSGs with different capacities is investigated.

An improved control strategy of VSG is proposed in Section II. The causes of power oscillation due to large virtual inertia are analyzed in Section III, two conditions are proposed: 1) increasing the difference coefficient, the damping coefficient, and the virtual inertia with the same proportion and 2) increasing the virtual synchronous reactance in an appropriate range. The parameter design method for parallel-operated VSGs with different capacities is investigated in Section IV. The effectiveness of the proposed method is verified through simulation and experiment in Sections V and VI, respectively, and the conclusions are presented in Section VII.

II. DESIGN OF A SELF-SYNCHRONIZED VSC

Fig. 1 shows a typical AC MG, in which the DGs and storage elements (SEs) are connected in parallel via the DC bus and then connected to the AC buses via self-synchronized VSCs. The VSG is constituted by DGs, SEs, and self-synchronized VSCs. For convenience, VSGs and self-synchronized VSCs are not distinguished in this study. VSGs often work in parallel to increase the capacity of an MG. Therefore, the analysis and discussion in this study focus on parallel-operated VSGs.

E1PWAX_2019_v19n2_580_f0001.png 이미지

Fig. 1. Typical configuration of a VSG-based MG.

Fig. 2 shows the control block diagram of VSGs. The VSG control strategy consists of active power frequency, virtual synchronous reactance, and voltage/current control loops. uabc is the output voltage of a VSG, ugabc denotes the voltage of the AC bus, iabc indicates the output current, and iLabc represents the inductor current. The VSG contains a parallel presynchronization unit based on PLL to connect the AC bus.

E1PWAX_2019_v19n2_580_f0019.png 이미지

Fig. 2. Control block diagram of the VSG parallel system.

A. Active Power Frequency Control

The active power frequency control of a VSG is reflected in the equations of prime-mover regulation and rotor motion [19]. The rotor motion is expressed as

\(\left\{\begin{array}{l} J \frac{\mathrm{d} \omega}{\mathrm{d} t}=T_{m}-T_{e}-D \Delta \omega=\frac{P_{m}}{\omega}-\frac{P_{e}}{\omega}-D\left(\omega-\omega_{N}\right) \\ \frac{\mathrm{d} \delta}{\mathrm{d} t}=\omega-\omega_{N} \end{array}\right.,\)       (1)

\(P_{m}=P_{r e f}+K_{\omega}\left(\omega_{r e f}-\omega\right),\)       (2)

where ω is the rotor angular velocity (rad/s), ωref represents the reference rotor angular velocity (rad/s), ωN refers to the rated rotor angular velocity (rad/s), Tm stands for the mechanical torque (N·m), Te indicates the electromagnetic torque (N·m), Pref denotes the reference power (W), Pm signifies the mechanical power (W), Pe is the electromagnetic power (W), D represents the damping coefficient, J denotes the virtual inertia (kg·m2), Kω indicates the load frequency coefficient, and δ refers to the power angle (rad).

In the case of a small excursion of angular speed, ω can be represented by ωN. Hence, Equation 1 can be changed to

\(\left\{\begin{array}{l} J \frac{\mathrm{d} \omega}{\mathrm{d} t} \approx T_{m}-T_{e}-D \Delta \omega=\frac{P_{m}}{\omega_{N}}-\frac{P_{e}}{\omega_{N}}-D\left(\omega-\omega_{N}\right) \\ \frac{\mathrm{d} \delta}{\mathrm{d} t}=\omega-\omega_{N} \end{array}\right..\)       (3)

Therefore, the relationship between the frequency and output active power of a VSG by using Equations 2 and 3 can be written as

\(\frac{\omega_{r e f}-\omega}{P_{r e f}-P}=-\frac{1}{J \omega_{N} s+\left(D \omega_{N}+K_{\omega}\right)}=-\frac{m}{\tau s+1},\)       (4)

where

\(\left\{\begin{array}{l} \tau=\frac{J \omega_{N}}{D \omega_{N}+K_{\omega}} \\ m=\frac{1}{D \omega_{N}+K_{\omega}} \end{array}\right..\)

B. Virtual Synchronous Reactance

The virtual synchronous reactance of a VSG is formed by synchronous generator's basic equation, as shown in Equation 5, which simulates the electrical part of the generator:

\(\left\{\begin{array}{l} u_{d}^{*}=-R i_{d}+L \frac{d i_{d}}{d t}+\omega L i_{q}+E \\ u_{q}^{*}=-R i_{q}+L \frac{d i_{q}}{d t}-\omega L i_{d} \end{array}\right.,\)       (5)

where ud* and uq* are the output voltage reference values in the dq frame; id and iq represent the direct and quadrature output currents, respectively; and L and R denote the virtual synchronous reactance (H) and resistance (Ω), respectively. However, R is set to zero so that the VSG's equivalent output impedance indicates induction.

A VSG has a reactive power voltage droop characteristic due to the existence of virtual synchronous reactance in the controller. Therefore, VSGs can realize reactive power sharing. On the basis of power flow calculation and by ignoring the virtual reactive power loss via virtual synchronous reactance, the relationship between voltage amplitude U and reactive power Q can be calculated as

\(U=E-\frac{\omega L}{E} Q.\)       (6)

The VSG is suitable for hierarchical distribution networks due to virtual synchronous reactance L, and the output voltage and reactive power of VSG have a droop relationship. Besides, the output reactive power can reduce AC bus voltage drop by removing the reactive power voltage droop control loop from the traditional droop control. Thus, the voltage quality can be improved compared with a conventional droop controller.

III. SUPPRESSION OF POWER OSCILLATION FOR PARALLEL-OPERATED VSGS

The small-signal model is established to analyze the power oscillation of two paralleled VSGs based on MGs, as shown in Fig. 3. E1∠δ1 and E2∠δ2 are the output voltages of VSG1 and VSG2, respectively, and U∠0 represents the AC bus voltage. As the line impedance is ignored in the analysis, Z1 and Z2 stand for the equivalent output impedances of VSG1 and VSG2, respectively, and the load is ZL.

E1PWAX_2019_v19n2_580_f0020.png 이미지

Fig. 3. Equivalent circuit of two-VSG systems.

The active power and angular frequency with small disturbance are calculated as

\(\left\{\begin{array}{l} \omega=\omega_{s}+\Delta \omega \\ P_{e}=P_{e s}+\Delta P, \end{array}\right.\)       (7)

where ωs and Pe are at the steady-state operating point and Δω indicate the small disturbances of the VSG output frequency and ΔP denotes the small disturbances of the VSG active power. Therefore, the state equation can be derived from Equations 2, 3, and 4 as follows:

\(\left\{\begin{array}{l} \frac{d \Delta \delta_{1}}{d t}=\Delta \omega_{1} \\ \frac{d \Delta \delta_{2}}{d t}=\Delta \omega_{2} \\ \frac{d \Delta \omega_{1}}{d t}=-\left(\frac{K_{\omega 1}}{J_{1} \omega_{N}}+\frac{D_{1}}{J_{1}}\right) \Delta \omega_{1}-\frac{\Delta P_{e 1}}{J_{1} \omega_{N}} \\ \frac{d \Delta \omega_{2}}{d t}=-\left(\frac{K_{\omega 2}}{J_{2} \omega_{N}}+\frac{D_{2}}{J_{2}}\right) \Delta \omega_{2}-\frac{\Delta P_{e 2}}{J_{2} \omega_{N}}, \end{array}\right.\)       (8)

where ∆δi refers to the power angle increment of each VSG. Then, the virtual electromagnetic power of the inverter can be calculated as

\(\left\{\begin{array}{l} \Delta P_{e 1}=S_{E 1} \Delta \delta_{12} \\ \Delta P_{e 2}=S_{E 2} \Delta \delta_{12}, \end{array}\right.\)       (9)

where

\(\left\{\begin{array}{l} S_{E I}=E_{1} E_{2}\left(-\left|G_{12}\right| \sin \delta_{12}+\left|B_{12}\right| \cos \delta_{12}\right) \\ S_{E 2}=E_{1} E_{2}\left(-\left|G_{12}\right| \sin \delta_{12}-\left|B_{12}\right| \cos \delta_{12}\right), \end{array}\right.\)       (10)

δ12 is the phase difference between VSG1 and VSG2, δ12 = δ1 − δ2, and G12 and B12 denote the conductance and susceptance between VSG1 and VSG2, respectively.

Therefore, the small-signal state-space model of two paralleled VSGs can be derived from Equations 7, 8, and 9 as

\(\left(\begin{array}{c} \frac{d \Delta \delta_{12}}{d t} \\ \frac{d \Delta \omega_{1}}{d t} \\ \frac{d \Delta \omega_{2}}{d t} \end{array}\right)=\left(\begin{array}{ccc} 0 & 1 & -1 \\ -\frac{S_{E 1}}{J_{1} \omega_{N}} & -\left(\frac{K_{\omega 1}}{J_{1} \omega_{N}}+\frac{D_{1}}{J_{1}}\right) & 0 \\ -\frac{S_{E 2}}{J_{2} \omega_{N}} & 0 & -\left(\frac{K_{\omega 2}}{J_{2} \omega_{N}}+\frac{D_{2}}{J_{2}}\right) \end{array}\right)\left(\begin{array}{c} \Delta \delta_{12} \\ \Delta \omega_{1} \\ \Delta \omega_{2} \end{array}\right).\)       (11)

Therefore, the characteristic equation can be obtained as

\(s^{3}+A s^{2}+B s+C=0,\)       (12)

where

\(\begin{array}{l} A=\frac{K_{\omega 1}}{J_{1} \omega_{\mathrm{N}}}+\frac{D_{1}}{J_{1}}+\frac{K_{\omega 2}}{J_{2} \omega_{\mathrm{N}}}+\frac{D_{2}}{J_{2}} \\ B=\frac{K_{\omega 1} K_{\omega 2}}{J_{1} J_{2} \omega_{N}^{2}}+\frac{D_{2} K_{\omega 1}}{J_{1} J_{2} \omega_{N}}+\frac{D_{1} K_{\omega 2}}{J_{1} J_{2} \omega_{N}}+\frac{D_{1} D_{2}}{J_{1} J_{2}}+\frac{S_{E 1}}{J_{1} \omega_{N}}-\frac{S_{E 2}}{J_{2} \omega_{N}} \\ C=\frac{1}{J_{1} J_{2} \omega_{N}^{2}}\left(K_{\omega 2} S_{E 1}-K_{\omega 1} S_{E 2}+D_{2} \omega_{N} S_{E 1}-D_{1} \omega_{N} S_{E 2}\right). \end{array}\)

Suppose that the same control parameters are adopted by VSG1 and VSG2. Then, the variables in Equation 12 can be simplified as follows:

\(\begin{array}{l} A^{\prime}=2\left(\frac{K_{\omega}}{J \omega_{0}}+\frac{D}{J}\right) \\ B^{\prime}=\frac{K_{\omega}^{2}}{J^{2} \omega_{0}^{2}}+\frac{2 D K_{\omega}}{J^{2} \omega_{0}}+\frac{D^{2}}{J^{2}}+\frac{1}{J \omega_{0}}\left(S_{\mathrm{E} 1}-S_{\mathrm{E} 2}\right) \\ C^{\prime}=\frac{1}{J^{2} \omega_{0}^{2}}\left[\left(K_{\mathrm{o}}+D \omega_{0}\right)\left(S_{\mathrm{E} 1}-S_{\mathrm{E} 2}\right)\right]. \end{array}\)

On the basis of Equations 4 and 12, the characteristic equation of the small-signal model can be reorganized as follows:

\(s^{3}+\frac{2}{\tau} s^{2}+\left[\frac{1}{\tau^{2}}+\frac{m}{\tau}\left(S_{E 1}-S_{E 2}\right)\right] s+\frac{m}{\tau^{2}}\left(S_{E 1}-S_{E 2}\right)=0.\)       (13)

As shown in Fig. 4, the equivalent impedance between two VSGs can be calculated as

\(Z_{12}=\frac{Z_{1} Z_{2}+Z_{1} Z_{L}+Z_{2} Z_{L}}{Z_{L}}=R_{12}+j X_{12},\)       (14)

where R12 and X12 are the resistive and inductive components of impedance, respectively.

E1PWAX_2019_v19n2_580_f0005.png 이미지

Fig. 4. Poles and zeros of: (a) D = 3, Kω = 1,000, L = 1 mH, J = 0.1 to 5. (b) D = 3, J = 1, L = 0.1 to 50 mH, Kω = 1,000. (c) J = 1, L = 1 mH, Kω = 1,000, D = 0 to 50. (d) D = 3, J = 1, L = 1 mH, Kω = 10 to 3,000. (e) L = 1 mH, m = 0.000001 to 0.01.

The load capacity and control parameters of two VSGs are identical. Therefore, hypothetically, δ12 ≈ 0, and load impedance ZL is set to 1 Ω. Furthermore, the relationship between SE1 and SE2 can be derived from Equation 10 to obtain

\(S_{E 1}-S_{E 2}=\frac{4 E_{1} E_{2}}{\omega^{3} L^{3}+4 \omega L}.\)       (15)

Then, the item of SE1 − SE2 in the characteristic equation can be replaced.

Fig. 4 presents the root locus analysis of the model. As shown in Fig. 4(a), the dominant pole becomes close to the imaginary axis, and the system becomes unstable with the increase in virtual inertia J. However, the dominant pole moves toward the real axis with large virtual synchronous reactance L, which indicates that the system oscillation is suppressed, as shown in Fig. 4(b). Similarly, the system becomes stable with the increase in D and Kω, as shown in Figs. 4(c) and (d), respectively. However, the dominant pole moves away from the real axis when the droop coefficient m increases and the system inertia time constant is unchanged, which indicates that the system damping is deteriorated. Moreover, when increasing inertial time constant τ, the dominant pole will move toward the imaginary axis, and the stability margin of the system is decreased, as shown in Fig. 4(e).

Therefore, two control strategies to suppress power oscillation are presented on the basis of the analysis:

(1) The inertial time should be kept constant with the increase of active frequency droop coefficient m, that is, D, J, and Kω should be enhanced with the same scale.

(2) Increasing virtual synchronous reactance L.

IV. PARALLEL-OPERATED VSGS WITH DIFFERENT CAPACITIES

Fig. 2 shows the principle of the proposed control strategy. Generated real power P and reactive power Q are [27], [28] as follows:

\(P_{i}=\frac{E_{i} U}{Z_{i}} \sin \delta_{i},\)       (16)

\(Q_{i}=\frac{E_{i} U}{Z_{i}} (\cos \delta_{i} - U),\)       (17)

where i = 1, 2.

In this study, E1 = E2 = E, and the line impedance is not considered. Therefore, Zi = ωLi, and the output active power of two VSGs is shown as follows

\(P_{i}=\frac{E U}{\omega L_{i}} \sin \delta_{i}.\)       (18)

With the increase in resistive load, the voltage drop caused by virtual impedance can be ignored compared with that caused by the reactive power on virtual synchronous reactance. The active power can be calculated as

\(P_{i}+\Delta P_{i}=\frac{E U}{\omega L_{i}} \sin \left(\delta_{i}+\Delta \delta_{i}\right).\)       (19)

Generally, δ is less than 30° for a synchronous generator based on Equations 16, 18, and 19. In this case, sin δ ≈ δ and cos δ ≈ 1. Thus, the dynamic active power of two VSGs can be calculated as

\(\Delta P_{i}=\frac{E U}{\omega L_{i}} \sin \left(\delta_{i}+\Delta \delta_{i}\right)-\frac{E U}{\omega L_{i}} \sin \left(\delta_{i}\right)=\frac{E U}{\omega L_{i}} \sin \Delta \delta_{i}\)       (20)

where ∆Pi (i = 1, 2) represent the dynamic active power for VSG1 and VSG2, respectively, and ∆δi (i = 1, 2) are the power angle variations. The load capacity ratio of VSG1 and VSG2 is set to C1: C2.

Sin ∆δ can be simplified to ∆δ because the power angle change is often relatively small. Therefore, we can rearrange Equation 20 as

\(\Delta P_{i}=\frac{E U}{\omega L_{i}} \Delta \delta_{i}.\)       (21)

The active power disturbance capacity ratio must remain constant when the load disturbance occurs to guarantee the active power sharing in the dynamic response, which is expressed as follows:

\(\frac{\Delta P_{1}}{\Delta P_{2}}=\frac{C_{1}}{C_{2}}.\)       (22)

Assuming that the resistive load increases at t0 and the time span in which the frequency of the common AC bus changes to steady-state is ∆t, the power angle increment of two VSGs can be calculated as

\(\Delta \delta_{i}=\int_{t_{0}}^{t_{0}+\Delta t} \Delta \omega_{i} d t.\)       (23)

Considering that Pref = 0, Equation 24 can be derived by combining Equations 4, 21, 22, and 23, which is expressed as follows:

\(\frac{\Delta \omega_{1}(s)}{\Delta \omega_{2}(s)}=\frac{E_{1} X_{2}}{E_{2} X_{1}} \frac{J_{1} \omega_{N} s+\left(D_{1} \omega_{N}+K_{\omega 1}\right)}{J_{2} \omega_{N} s+\left(D_{2} \omega_{N}+K_{\omega 2}\right)}=\frac{C_{1}}{C_{2}}.\)       (24)

The active power disturbance capacity ratio of VSGs in the dynamic response remains constant. The control parameters of two VSGs are determined on the basis of the following:

\(\frac{L_{2}}{L_{1}}=\frac{J_{1}}{J_{2}}=\frac{D_{1}}{D_{2}}=\frac{K_{\omega 1}}{K_{\omega 2}}=\frac{C_{1}}{C_{2}}.\)       (25)

On the basis of Equations 6 and 17, the output reactive power of two VSGs when the inductive load increases can be calculated as

\(Q_{i}=\frac{E\left(U-\frac{\omega L_{i}}{E} \Delta Q_{i}\right)}{\omega L_{i}}\left[\cos \delta-\left(U-\frac{\omega L_{i}}{E} \Delta Q_{i}\right)\right],\)       (26)

where ∆Qi (i = 1, 2) are the dynamic reactive power values for VSG1 and VSG2, respectively.

Therefore, dynamic reactive power ∆Qi based on Equations 17 and 26 can be calculated as

\(\Delta Q_{i}=(2 U-\cos \delta-1) \frac{E}{\omega L_{i}}.\)       (27)

Considering that the reactive power disturbance capacity ratio must be kept constant when load disturbance occurs, we have

\(\frac{\Delta Q_{1}}{\Delta Q_{2}}=\frac{C_{1}}{C_{2}}.\)       (28)

The reactive power disturbance capacity ratio of VSGs in the dynamic response remains constant based on Equations 27 and 28, and the control parameters of two VSGs can be calculated as

\(\frac{L_{2}}{L_{1}}=\frac{C_{1}}{C_{2}}.\)       (29)

On the basis of Equations 25 and 29 designed parameters of VSGs with different capacities, the dynamic response of active and reactive power allocated by capacity can be improved.

V. SIMULATION RESULT

To verify the accuracy of the proposed control and improved parameter design methodology, a simulation model of two paralleled VSGs is established in MATLAB/Simulink. Table I lists the parameters of the controller and power stage.

TABLE I PARAMETERS OF SIMULATIONS

E1PWAX_2019_v19n2_580_t0001.png 이미지

A. Simulation for Pre-synchronization

The phase of the VSG must be the same as that of the AC bus before the VSG is connected to the AC bus, and the pre-synchronization for parallel-operated VSGs must be determined. Fig. 5 shows the simulated voltage waveforms of a VSG and AC bus. The voltage phase of the VSG can be synchronized to the AC bus by using a pre-synchronization method.

E1PWAX_2019_v19n2_580_f0007.png 이미지

Fig. 5. Simulation waveform for VSG output voltage and AC bus voltage.

B. Simulation for Power Oscillation Suppression

As shown in Fig. 6, VSG1 normally operates with a resistive load of approximately 400 W in the beginning, and VSG2 is switched at 1.0 s. Active power oscillation remarkably increases with the increase in virtual inertia J when the VSG is connected to the AC bus by keeping all other parameters unchanged. This condition is consistent with the small-signal analysis, Fig. 4(a) shows that large virtual inertia affects the stability of the VSG parallel system.

E1PWAX_2019_v19n2_580_f0008.png 이미지

Fig. 6. Simulation waveforms for active power with different J values.

Fig. 7 shows the same scenario is used to determine the power damping influence of droop coefficient m. The overshoot and oscillation of the output active power are effectively suppressed with a large damping coefficient and virtual inertia when the damping coefficient, virtual inertia, and load frequency coefficient are accordingly adjusted with the same time constant τ. This finding is consistent with the results of smallsignal analysis, as shown in Fig. 4(e). System damping is enhanced, and the oscillation is suppressed with large values of J, D, and Kω.

E1PWAX_2019_v19n2_580_f0009.png 이미지

Fig. 7. Simulation waveforms for active power with different m values.

As shown in Fig. 8, the time of frequency curve recovery to steady state is unchanged with the increase of J, D, and Kω, which is approximately 1.0 s. The stability of frequency is not affected.

E1PWAX_2019_v19n2_580_f0010.png 이미지

Fig. 8. Simulation waveforms for AC bus frequency with different parameters.

The same scenario is used to evaluate the power damping influence of virtual synchronous reactance L, as shown in Fig. 9. Active power fluctuation is smoothed and suppressed with large values of L. This condition is consistent with the results of small-signal analysis. As shown in Fig. 4(b), system damping is enhanced, and the oscillation is suppressed with large values of L.

E1PWAX_2019_v19n2_580_f0011.png 이미지

Fig. 9. Simulation waveforms for active power with different L values.

C. Simulation for Parallel-operated VSGs with Different Capacities

Simulation: Parallel-operated VSGs with resistive loads operate, resistive and inductive loads are added at 1.5 s, and the load capacity ratio is set to C1:C2 = 5:3. Table II lists the simulation parameters.

TABLE II PARAMETERS OF SIMULATIONS

E1PWAX_2019_v19n2_580_t0002.png 이미지

As shown in Fig. 10(a), the active power sharing performance is improved after the optimization of parameters. The system is moved to stable operation point within 0.1 s, whereas the original system takes approximately 2.5 s. Moreover, the dynamic reactive power sharing performance is improved, as exhibited in Fig. 10(b). As presented in Fig. 10(c), the relative frequency difference between VSG1 and VSG2 is not zero due to the non-optimized parameters in 1–3 s.

E1PWAX_2019_v19n2_580_f0012.png 이미지

Fig. 10. Simulation waveforms of the VSG parallel system with different load capacities. (a) Active power of parallel-operated VSGs. (b) Reactive power of parallel-operated VSGs. (c) Frequency of parallel-operated VSGs.

VI. EXPERIMENTAL RESULTS

In order to verify the effectiveness of the control strategy, an experimental platform was developed, which includes two VSGs, an LC filter, a DSP control board, and a resistive load, as shown in Fig. 11. Table III lists the parameters of the power stage.

E1PWAX_2019_v19n2_580_f0013.png 이미지

Fig. 11. Experimental platform of parallel-operated VSGs. (a) Schematic of the experimental setup. (b) Photo of the experimental platform.

TABLE III PARAMETERS OF EXPERIMENTS

E1PWAX_2019_v19n2_580_t0003.png 이미지

A. Contrast Experiment for Traditional Droop and VSG Controls

The experimental waveforms of the VSG and AC bus is shown as Fig. 12. The voltage of the VSG can be synchronized to the AC bus through pre-synchronization.

E1PWAX_2019_v19n2_580_f0014.png 이미지

Fig. 12. Experimental waveform of the output voltage of VSG and AC bus voltage.

B. Power Oscillation Suppression for Parallel-operated VSGs

Figs. 13 and 14 present the operation waveforms of VSGs with different values of J, D, and Kω. These waveforms are obtained by using a D/A converter. As shown in Figs. 13(a), (b), and (c), a large active power oscillation occurs in the parallel system when VSG2 is connected to the AC bus with small J, D, and Kω values. However, the output active power of VSGs becomes smooth with the increase of J, D, and Kω. The inertia time is kept constant. Figs. 14(a), (b), and (c) show the frequency waveforms of the AC bus. The parameters are consistent with Figs. 13(a), (b), and (c), in which the active load increases and active power decreases. In addition, the time of frequency curve recovery to steady state is unchanged with the increase of J, D, and Kω at the same rate, which is approximately 1 s. Hence, the frequency stability of the AC bus is unchanged. Thus, the increase of J, D, and Kω at the same rate can decrease the power oscillation and overshoot, and improve the parallel system's stability.

E1PWAX_2019_v19n2_580_f0015.png 이미지

Fig. 13. Output active power experimental waveforms of VSG parallel system. (a) D = 2, J = 0.5, Kω = 100, L = 0.006. (b) D = 3, J = 0.75, Kω = 150, L = 0.006. (c) D = 4, J = 1, Kω = 200, L = 0.006.

E1PWAX_2019_v19n2_580_f0016.png 이미지

Fig. 14. Frequency waveforms of the AC bus .(a) D = 2, J = 0.5, Kω = 100, L = 0.006. (b) D = 3, J = 0.75, Kω = 150, L = 0.006. (c) D = 4, J = 1, Kω = 200, L = 0.006.

Fig. 15 exhibits the operation waveforms of VSGs with different values of L. A large power oscillation occurs in the parallel system when VSG2 is connected to the AC bus with a small virtual synchronous reactance. However, the output active power of VSGs becomes smooth with the increase of virtual synchronous reactance. Thus, the increase of L can reduce the overshoot, suppress the power oscillation, and improve the parallel system's stability. The increase of L is effective for suppressing power oscillations.

E1PWAX_2019_v19n2_580_f0017.png 이미지

Fig. 15. Output active power experimental waveforms of VSG parallel system. (a) D = 3, J = 0.6, Kω = 150, L = 0.005. (b) D = 3, J = 0.6, Kω = 150, L = 0.008. (c) D = 3, J = 0.6, Kω = 150, L = 0.01.

C. Parallel-operated VSGs with Different Capacities

The design methodology of parallel-operated VSGs with different capacities is verified by experiment. Table IV lists the controller parameters.

TABLE IV OPTIMIZED SYSTEM PARAMETERS FOR THE EXPERIMENT

E1PWAX_2019_v19n2_580_t0004.png 이미지

Fig. 16(a) presents the experimental waveforms of load step disturbance without optimized parameters. The current amplitude of VSG1 slowly increases. However, the VSG2 current slowly decreases, and the output currents of the two VSGs’ dynamic current distributions are not proportional to the load capacity ratio.

E1PWAX_2019_v19n2_580_f0018.png 이미지

Fig. 16. Experimental waveforms of VSG parallel system with different load capacities. (a) Parameters are not optimized. (b) Parameters are optimized.

Fig. 16(b) shows the experimental waveforms of load step disturbance with optimized parameters. The dynamic currents of the two VSGs meet the load capacity ratio when a disturbance occurs, which is consistent with the simulation results. The experimental results verify the effectiveness of the designed parameters for parallel-operated VSGs with different capacities.

VII. CONCLUSION

This study investigates the mechanism of power oscillation in parallel-operated VSGs on the basis of a small signal model, and an improved control strategy is presented without altering inertia time constant to suppress power oscillation: 1) The inertial time should be kept constant with the increase of active frequency droop coefficient m, that is, D, J, and Kω should be enhanced with the same scale; 2) Increasing virtual synchronous reactance L. The simulation and experiment results verify the effectiveness of the proposed control strategy. The frequency stability of the common AC bus is improved, and power oscillations are suppressed. Moreover, the dynamic power sharing precision is effectively enhanced by using an improved parameter design method for parallel-operated VSGs with different capacities.

References

  1. H.-P. Beck and R. Hesse, "Virtual synchronous machine," in Proceedings of the 9th International Conference on Electrical Power Quality and Utilisation, pp. 1-6, 2007.
  2. F. Gao and M. R. Iravani, “A control strategy for a distributed generation unit in grid-connected and autonomous modes of operation,” IEEE Trans. Power Del., Vol. 23, No. 2, pp. 850-859, Apr. 2008. https://doi.org/10.1109/TPWRD.2007.915950
  3. J. Alipoor, Y. Miura, and T. Ise, "Distributed generation grid integration using virtual synchronous generator with adoptive virtual inertia," 2013 IEEE Energy Conversion Congress and Exposition, pp. 4546-4552, 2013.
  4. T. Shintai, Y. Miura, and T. Ise, “Oscillation damping of a distributed generator using a virtual synchronous generator,” IEEE Trans. Power Del., Vol. 29, No. 2, pp. 668-676, Apr. 2014. https://doi.org/10.1109/TPWRD.2013.2281359
  5. J. Alipoor, Y. Miura, and T. Ise, “Power system stabilization using virtual synchronous generator with alternating moment of inertia,” IEEE J. Emerg. Sel. Topics Power Electron., Vol. 3, No. 2, pp. 451-458, Jun. 2015. https://doi.org/10.1109/JESTPE.2014.2362530
  6. Q. C. Zhong and G. Weiss, "Static synchronous generators for distributed generation and renewable energy," Power Systems Conference and Exposition, 2009. PSCE '09. IEEE/PES, pp. 1-6, 2009.
  7. M. John, P. A. Mendoza-Araya, and G. Venkataramanan, "Small signal impedance measurement in droop controlled AC microgrids," 2014 IEEE Energy Conversion Congress and Exposition (ECCE), pp. 702-709, 2014.
  8. J. C. Vasquez, J. M. Guerrero, M. Savaghebi, and R. Teodorescu, "Modeling, analysis, and design of stationary reference frame droop controlled parallel three-phase voltage source inverters," Power Electronics and ECCE Asia (ICPE & ECCE), 2011 IEEE 8th International Conference on, pp. 272-279, 2011.
  9. J. C. Vasquez, J. M. Guerrero, M. Savaghebi, J. Eloy-Garcia, and R. Teodorescu, “Modeling, analysis, and design of stationary-reference-frame droop-controlled parallel threephase voltage source inverters,” IEEE Trans. Ind. Electron., Vol. 60, No. 4, pp. 1271-1280, Apr. 2013. https://doi.org/10.1109/TIE.2012.2194951
  10. T. Wu, Z. Liu, J. Liu, S. Wang, and Z. You, “A unified virtual power decoupling method for droop-controlled parallel inverters in microgrids,” IEEE Trans. Power Electron., Vol. 31, No. 8, pp. 5587-5603, Aug, 2016. https://doi.org/10.1109/TPEL.2015.2497972
  11. Q. C. Zhong and G. Weiss, “Synchronverters: Inverters that mimic synchronous generators,” IEEE Trans. Ind. Electron., Vol. 58, No. 4, pp. 1259-1267, Apr. 2011. https://doi.org/10.1109/TIE.2010.2048839
  12. Q. C. Zhong, P. L. Nguyen, Z. Ma, and W. Sheng, “Selfsynchronized synchronverters: inverters without a dedicated synchronization unit,” IEEE Trans. Power Electron., Vol. 29, No. 2, pp. 617-630, Feb. 2014. https://doi.org/10.1109/TPEL.2013.2258684
  13. C. Cheng, Z. Zeng, H. Yang, and R. Zhao, "Wireless parallel control of three-phase inverters based on virtual synchronous generator theory," Electrical Machines and Systems (ICEMS), 2013 International Conference on, pp. 162-166, 2013.
  14. M. Guan, W. Pan, J. Zhang, Q. Hao, J. Cheng, and X. Zheng, “Synchronous generator emulation control strategy for voltage source converter (VSC) stations,” IEEE Trans. Power Syst, Vol. 30, No. 6, pp. 3093-3101, Nov. 2015. https://doi.org/10.1109/TPWRS.2014.2384498
  15. Y. Hirase, O. Noro, K. Sugimoto, K. Sakimoto, Y. Shindo, and T. Ise, "Effects of suppressing frequency fluctuations by parallel operation of virtual synchronous generator in microgrids," 2015 IEEE Energy Conversion Congress and Exposition (ECCE), pp. 3694-3701, 2015.
  16. Y. Hirase, O. Noro, K. Sugimoto, K. Sakimoto, Y. Shindo, and T. I. Osaka, "Effects and analysis of suppressing frequency fluctuations in microgrids using virtual synchronous generator control," Industrial Electronics Society, IECON 2015 - 41st Annual Conference of the IEEE, pp. 43-49, 2015.
  17. J. Alipoor, Y. Miura, and T. Ise, "Distributed generation grid integration using virtual synchronous generator with adoptive virtual inertia," 2013 IEEE Energy Conversion Congress and Exposition, pp. 4546-4552, 2013.
  18. J. Alipoor, Y. Miura, and T. Ise, “Power system stabilization using virtual synchronous generator with alternating moment of inertia,” IEEE J. Emerg. Sel. Topics Power Electron., Vol. 3, No. 2, pp. 451-458, Jun. 2015. https://doi.org/10.1109/JESTPE.2014.2362530
  19. Y. Du, J. Su, M. Mao, and X. Yang, "Autonomous controller based on synchronous generator dq0 model for micro grid inverters," IEEE ECCE Asia, pp. 2645-2649, 2011.
  20. J. Liu, Y. Miura, and T. Ise, “Comparison of dynamic characteristics between virtual synchronous generator and droop control in inverter-based distributed generators,” IEEE Trans. Power Electron., Vol. 31, No. 5, pp. 3600-3611, May 2016. https://doi.org/10.1109/TPEL.2015.2465852
  21. J. Liu, Y. Miura, and T. Ise, "Dynamic characteristics and stability comparisons between virtual synchronous generator and droop control in inverter-based distributed generators," 2014 International Power Electronics Conference (IPECHiroshima 2014 - ECCE ASIA), pp. 1536-1543, 2014.
  22. H. Wu, X. Ruan, D. Yang, X. Chen, W. Zhao, Z. Lv, and Q.-C. Zhong, “Small-signal modeling and parameters design for virtual synchronous generators,” IEEE Trans. Ind. Electron., Vol. 63, No. 7, pp. 4292-4303, Jul. 2016. https://doi.org/10.1109/TIE.2016.2543181
  23. J. M. Guerrero, L. G. de Vicuna, J. Matas, M. Castilla, and J. Miret, “A wireless controller to enhance dynamic performance of parallel inverters in distributed generation systems,” IEEE Trans. Power Electron., Vol. 19, No. 5, pp. 1205-1213, Sep. 2004. https://doi.org/10.1109/TPEL.2004.833451
  24. R. M. Santos Filho, P. F. Seixas, P. C. Cortizo, G. Gateau, and E. A. A. Coelho, "Power system stabilizer for comunicationless parallel connected inverters," 2010 IEEE International Symposium on Industrial Electronics, pp. 1004-1009, 2010.
  25. W. A. Habri and Y. L. A. Magid, "Power system stabilizer for power distribution control of parallel inverters in a grid -Connected micro-grid system," Innovative Smart Grid Technologies-Middle East (ISGT Middle East), 2011 IEEE PES Conference on, pp. 1-8, 2011.
  26. P. Rodriguez, C. Citro, I. Candela, J. Rocabert, and P. Rodriguez, "Flexible grid connection and islanding of SPCbased PV power converters," IEEE Trans. Ind. Appl. Vol. 54, No. 3, pp. 2690-2702, May/Jun. 2018. https://doi.org/10.1109/TIA.2018.2800683
  27. T. Wildi, Electrical Machines, Drives and Power Systems, 6th ed. EnglewoodCliffs, NJ, USA: Prentice-Hall, 2005.B.
  28. B. Singh, R. Saha, A. Chandra, and K. Al-Haddad, "Static synchronous compensators (STATCOM): A review," IET Power Electron., Vol. 2, No. 4, pp. 297-324, Jul, 2009. https://doi.org/10.1049/iet-pel.2008.0034
  29. Y. Wei, H. Zhang, Q. Song, and K. Sun, "Control strategy for parallel-operated virtual synchronous generators," 2016 International Power Electronics and Motion Control Conference (IPEMC-ECCE Asia 2016), pp. 2015-2021, 2016.