# Optimization Control Strategy for Islanded Parallel Virtual Synchronous Generators

Authors: Hailin Hu; Fu Feng; Tao Wang; Xiaofeng Wan; Xiaohua Ding
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

Virtual Synchronous Generators (VSGs) can reduce the frequency and power oscillation in the grid. For parallel multi-VSGs in island microgrid, the differences of equivalent output impedance and line impedance affect the orderly sharing of reactive power and restraining of circulating current greatly. The power sharing and circulating current characteristics of VSGs are analyzed in this paper. To solve this problem, an optimization control method for parallel VSGs is proposed, which included the inner voltage and current loop and the outer reactive power loop. In the inner voltage and current loop, a virtual complex impedance including a resistive component and an inductive component is introduced. It reduces the coupling between active power and reactive power, and it reduces the disorderly distribution of power due to impedance differences. In the reactive power loop, the output voltage feedback and integrator links are adopted, which reduce the deviation of output voltage and restrain circulating current. The effects on the equivalent output impedance in the different control parameters and virtual complex impedance are analyzed, the effects on system stability in the different resistive components of virtual complex impedance are analyzed, and the proper parameters are selected. Simulation and experimental results show the correctness and validity of the proposed control method.


---

## Body

## 1. Introduction

By increasing the inertia and damping links in an active power control loop, the grid support capability in a weak grid of virtual synchronous generators (VSGs) is summarized in [1, 2]. When multiple VSGs are operated in parallel, there are differences in the output voltage, equivalent output impedance, and connection line impedance of the individual VSGs, which are caused by the different controller parameters, nonlinearities in the power devices, differences in filter parameters, and the distance from the grid connection point. These differences affect the accuracy of power distribution, increase the fundamental current circulation between the VSGs, and may cause overcurrent faults [3]. Therefore, it is of great value to study the power distribution and current circulation suppression between parallel VSGs. [4, 5].Focusing on the power distribution problem of parallel VSGs, in [6], an integral controller is added in the reactive power loop to eliminate the influence of connection line impedance on reactive power distribution. However, the coupling of active and reactive power is not considered. In [7], a feedforward angular frequency is added to the power frequency link to realize switching under different operation modes. A small-signal model and several constraints are used to optimize the dynamic and static characteristics of the parameters, but the implementation method is complicated. In [8], function terms for the power and droop coefficient in traditional droop control are introduced. This improved the dynamic regulation performance, but it did not eliminate the steady-state error caused by the connection line impedance. In [9], a new multiloop control strategy is described, the virtual reactance can be selected according to the external inductance, but the conturation principle of virtual impedance in parallel VSGs is not considered. In [10], a joint control method for the virtual inductor and external output inductor is proposed, which is used to improve the applicability of the inverter P-f and Q-V droop control in medium- and low-voltage networks. However, it is difficult to obtain the external output inductance accurately, and the method is not practical. In [11], a parallel VSGs control method based on virtual impedance is proposed. This approach improves the excitation link to achieve adequate control of the AC bus voltage, but the power coupling problem is not considered. In [12], a consistency-based parallel VSGs distributed control method is described. Although this enhances the system frequency and active power characteristics, it requires communication between the parallel VSGs.This paper describes a method that improves the traditional virtual synchronous power generation technology, based on analysis of the parallel operation characteristics of multiple VSGs. The optimization control strategy includes the outer reactive power loop and the inner voltage and current loop. In the inner voltage and current loop, a virtual complex impedance including a resistive component and an inductive component is introduced. It reduces the effects on decoupling between active power and reactive power, and it reduces the effects on accurate power sharing due to the impedance difference. In the reactive power loop, the output voltage feedback and integrator links are adopted, which will reduce the deviation of output voltage and restrain circulating current.


## 2. Operating Characteristics of Parallel VSGs

Figure1 shows the simplified structure of a traditional parallel VSGs system. For a single VSG, it connects to the AC bus through the connection line and grid-connected switch. Zline is the impedance of the connection line, and K is the grid-connected switch. Zload is the common load, L1 is the filter inductance, and R1 is the sum of the equivalent internal resistance of the filter inductance and the equivalent internal resistance of the power device. UC is the VSG output capacitor voltage, also known as the terminal voltage. I2 is the grid-side output current, and I1 is the converter-side current.Figure 1
Structure diagram of parallel VSGs.A simple model with two parallel VSGs is shown in Figure2. U1 and U2 are the voltage amplitudes of the VSG bridge arms’ midpoint, E is the AC bus voltage amplitude, Zo is the equivalent output impedance of the VSG Zo=Ro+jXo, and Zline is the connection line impedance Zline=Rline+jXline.Figure 2
Two VSG parallel models.The transmission impedanceZ is equal to the sum of the equivalent output impedance and the connection line impedance, φ is the transmission impedance phase angle, and θ is the phase angle difference between the bus voltage and the voltage of the bridge arm’s midpoint. Io and IL represent the output current of the VSG and load current.


### 2.1. Analysis of VSG Power Transmission Characteristics

The output active power and reactive power of VSGii=1,2 can be expressed as follows:(1)Pi=UiEZicosθi−E2Zicosφi+UiEZisinθisinφi,Qi=UiEZicosθi−E2Zisinφi−UiEZisinθisinφi.According to (1), the output active power and reactive power are related to the voltage amplitude and phase angle difference and are affected by the transmission impedance. By setting the voltage and current control loop parameters, the equivalent output impedance of the VSG can be controlled, so that the equivalent output inductive reactance is much larger than the equivalent output resistive impedance, Xo≫Ro. The transmission impedance of different voltage levels has different impedance ratios. In low-voltage distribution networks, the resistive impedance is much larger than the inductive impedance, Xline≪Rline. In medium-voltage distribution networks, Xline≈Rline. This causes the transmission impedance to exhibit complex characteristics. To decouple the active power from the reactive power, the transmission impedance is set to be inductive using virtual impedance technology. The power transmission is approximately shown as follows:(2)Pi=UiEXoi+Xlineiθi,Qi=EUi−EXoi+Xlinei.From the abovementioned analysis, when the transmission impedance is inductive, the active power transmission is mainly affected by the voltage phase angle, whereas the reactive power transmission is affected by the voltage amplitude. [13] The active power and reactive power are approximately decoupled.The VSG power control equation can be described as follows:(3)Jdωdt=P∗ωn−Piωn+Dpωn−ωi,Ui=En+DqQ∗−Qi,where J is the virtual inertia, P∗ is the reference value of active power, Pi is the actual output active power, DP is the damping coefficient, ωn is the rated angular frequency, and ωi is the virtual angular frequency. Ui is the amplitude of VSG output voltage, En is the rated voltage amplitude, Dq is the reactive power droop coefficient, Q∗ is the reference value of reactive power, and Qi is the actual output reactive power. Combining equations (2) and (3), the VSG active power closed-loop control diagram is shown as Figure 3. The reactive power closed-loop control diagram is shown as Figure 4.Figure 3
Active power closed-loop control diagram of the VSG.Figure 4
Reactive power closed-loop control diagram of the VSG.There is an integral link in the active loop of the virtual synchronous control strategy; when the VSG runs in the steady state, the input of the integral link is 0.(4)Pi=P∗+Dpωn−ωi.According to (4), when the transmission impedance of the VSG is inductive, the active power distribution is related to the active power command value and the active-frequency droop coefficient, and active power distribution is not affected by transmission impedance.For two VSGs with the same rated capacity, as long as the active power command value and the active-frequency droop coefficient are equal, the active power can be divided equally. Therefore, the active power transmission is highly robust.Similarly, the VSG output reactive power equation can be obtained as follows:(5)Qi=En−E+Q∗DqE/Xi1+E/XiDq.According to (5), the reactive power distribution is affected by the droop coefficient, the reactive power reference value, and the transmission impedance. When the reactive power of two VSGs with equal capacity is evenly distributed, the droop coefficients, reactive power setting values, and transmission impedances of each VSG must be equal. When the actual system is running, the system interference, the nonlinearity of the power device, the parameter offset, and the difference in connection impedance make it difficult to meet the conditions for the orderly distribution of reactive power.


### 2.2. Analysis of Current Circulation Characteristics of the VSG

The output current characteristics of parallel VSGs are analyzed in the fundamental wave frequency domain. Equation (6) can be obtained from Figure 2.(6)U1∠θ1−Zline1Io1=E∠0,U2∠θ2−Zline2Io2=E∠0,Io1+Io2=IL=E∠0Zload.In the abovementioned equations, the transmission impedance is the sum of the equivalent output impedance and the connection line impedance,Zi=Zoi+Zlinei. If the equivalent output impedance and connection line impedance of the two VSGs are equal, the following equation can be obtained:(7)Io1=U1∠θ1−U2∠θ22Z1+E2Zload,Io2=U2∠θ2−U1∠θ12Z1+E2Zload.When the equivalent output impedance of the VSG is inductive and the connection line impedance is far less than the equivalent output impedance, the transmission impedance is approximately equal to the equivalent output impedance,Zi≈Xoi, and the following equation can be obtained:(8)Io1=U1∠θ1−U2∠θ22jXo1+E2Zload=IH1+IoL1,Io2=U2∠θ2−U1∠θ22jXo1+E2Zload=IH2+IoL2.According to equation (8), as the equivalent output impedance of the two VSGs is inductive and equal and the connection line impedance is much smaller than the equivalent output impedance, the output current of the VSG includes two components: the load current IoL and the circulation current IH. The load current is distributed equally, and the circulating current is determined by the output voltage amplitude difference, phase difference, and equivalent output impedance of the parallel VSGs. [13].The fundamental circulating current between two VSGs is defined as follows:(9)IH=Io1−Io22=U1∠θ1−U2∠θ22jXo1.If the virtual voltage amplitude of the VSGs is equal, there is only a phase difference, and the phase difference will cause active power circulation, which will transfer active current circulation from the advanced phase part to the delayed phase part in the parallel VSGs. If the virtual voltage phase of the VSG is equal, there is only an amplitude difference, and the reactive current circulation will flow from the high-voltage amplitude part to the low-voltage amplitude part. The output impedance of the VSG with high-voltage amplitude is inductive, and the output impedance of the VSG with low-voltage amplitude is capacitive. If both the amplitude and phase of the output voltages of the two VSGs are not equal, then the circulating current contains both active and reactive components [14].From the analysis of the current circulation characteristics, it can be seen that the current circulation of the parallel VSGs is suppressed by achieving the output voltage and transmission impedance of the two VSGs is equal. The suppression of the current circulation can also be realized by increasing the equivalent output impedance.


## 2.1. Analysis of VSG Power Transmission Characteristics

The output active power and reactive power of VSGii=1,2 can be expressed as follows:(1)Pi=UiEZicosθi−E2Zicosφi+UiEZisinθisinφi,Qi=UiEZicosθi−E2Zisinφi−UiEZisinθisinφi.According to (1), the output active power and reactive power are related to the voltage amplitude and phase angle difference and are affected by the transmission impedance. By setting the voltage and current control loop parameters, the equivalent output impedance of the VSG can be controlled, so that the equivalent output inductive reactance is much larger than the equivalent output resistive impedance, Xo≫Ro. The transmission impedance of different voltage levels has different impedance ratios. In low-voltage distribution networks, the resistive impedance is much larger than the inductive impedance, Xline≪Rline. In medium-voltage distribution networks, Xline≈Rline. This causes the transmission impedance to exhibit complex characteristics. To decouple the active power from the reactive power, the transmission impedance is set to be inductive using virtual impedance technology. The power transmission is approximately shown as follows:(2)Pi=UiEXoi+Xlineiθi,Qi=EUi−EXoi+Xlinei.From the abovementioned analysis, when the transmission impedance is inductive, the active power transmission is mainly affected by the voltage phase angle, whereas the reactive power transmission is affected by the voltage amplitude. [13] The active power and reactive power are approximately decoupled.The VSG power control equation can be described as follows:(3)Jdωdt=P∗ωn−Piωn+Dpωn−ωi,Ui=En+DqQ∗−Qi,where J is the virtual inertia, P∗ is the reference value of active power, Pi is the actual output active power, DP is the damping coefficient, ωn is the rated angular frequency, and ωi is the virtual angular frequency. Ui is the amplitude of VSG output voltage, En is the rated voltage amplitude, Dq is the reactive power droop coefficient, Q∗ is the reference value of reactive power, and Qi is the actual output reactive power. Combining equations (2) and (3), the VSG active power closed-loop control diagram is shown as Figure 3. The reactive power closed-loop control diagram is shown as Figure 4.Figure 3
Active power closed-loop control diagram of the VSG.Figure 4
Reactive power closed-loop control diagram of the VSG.There is an integral link in the active loop of the virtual synchronous control strategy; when the VSG runs in the steady state, the input of the integral link is 0.(4)Pi=P∗+Dpωn−ωi.According to (4), when the transmission impedance of the VSG is inductive, the active power distribution is related to the active power command value and the active-frequency droop coefficient, and active power distribution is not affected by transmission impedance.For two VSGs with the same rated capacity, as long as the active power command value and the active-frequency droop coefficient are equal, the active power can be divided equally. Therefore, the active power transmission is highly robust.Similarly, the VSG output reactive power equation can be obtained as follows:(5)Qi=En−E+Q∗DqE/Xi1+E/XiDq.According to (5), the reactive power distribution is affected by the droop coefficient, the reactive power reference value, and the transmission impedance. When the reactive power of two VSGs with equal capacity is evenly distributed, the droop coefficients, reactive power setting values, and transmission impedances of each VSG must be equal. When the actual system is running, the system interference, the nonlinearity of the power device, the parameter offset, and the difference in connection impedance make it difficult to meet the conditions for the orderly distribution of reactive power.


## 2.2. Analysis of Current Circulation Characteristics of the VSG

The output current characteristics of parallel VSGs are analyzed in the fundamental wave frequency domain. Equation (6) can be obtained from Figure 2.(6)U1∠θ1−Zline1Io1=E∠0,U2∠θ2−Zline2Io2=E∠0,Io1+Io2=IL=E∠0Zload.In the abovementioned equations, the transmission impedance is the sum of the equivalent output impedance and the connection line impedance,Zi=Zoi+Zlinei. If the equivalent output impedance and connection line impedance of the two VSGs are equal, the following equation can be obtained:(7)Io1=U1∠θ1−U2∠θ22Z1+E2Zload,Io2=U2∠θ2−U1∠θ12Z1+E2Zload.When the equivalent output impedance of the VSG is inductive and the connection line impedance is far less than the equivalent output impedance, the transmission impedance is approximately equal to the equivalent output impedance,Zi≈Xoi, and the following equation can be obtained:(8)Io1=U1∠θ1−U2∠θ22jXo1+E2Zload=IH1+IoL1,Io2=U2∠θ2−U1∠θ22jXo1+E2Zload=IH2+IoL2.According to equation (8), as the equivalent output impedance of the two VSGs is inductive and equal and the connection line impedance is much smaller than the equivalent output impedance, the output current of the VSG includes two components: the load current IoL and the circulation current IH. The load current is distributed equally, and the circulating current is determined by the output voltage amplitude difference, phase difference, and equivalent output impedance of the parallel VSGs. [13].The fundamental circulating current between two VSGs is defined as follows:(9)IH=Io1−Io22=U1∠θ1−U2∠θ22jXo1.If the virtual voltage amplitude of the VSGs is equal, there is only a phase difference, and the phase difference will cause active power circulation, which will transfer active current circulation from the advanced phase part to the delayed phase part in the parallel VSGs. If the virtual voltage phase of the VSG is equal, there is only an amplitude difference, and the reactive current circulation will flow from the high-voltage amplitude part to the low-voltage amplitude part. The output impedance of the VSG with high-voltage amplitude is inductive, and the output impedance of the VSG with low-voltage amplitude is capacitive. If both the amplitude and phase of the output voltages of the two VSGs are not equal, then the circulating current contains both active and reactive components [14].From the analysis of the current circulation characteristics, it can be seen that the current circulation of the parallel VSGs is suppressed by achieving the output voltage and transmission impedance of the two VSGs is equal. The suppression of the current circulation can also be realized by increasing the equivalent output impedance.


## 3. Optimization Control Strategy of Parallel VSGs

### 3.1. Introduction of Virtual Impedance

From the analysis in the previous section, the active power and the reactive power are decoupled by setting the transmission impedance to be purely resistive or purely inductive, and this is the premise of independent control of active power and reactive power. Increasing the equivalent output impedance of the virtual synchronous generator can reduce the current circulation caused by the difference between the equivalent output impedance and the connection impedance. To change the transmission impedance characteristics of the VSG, the output current of the grid side is used as the feedback, and a virtual complex impedance consisting of resistance and inductive components is introduced. The negative resistive virtual impedance is used to reduce the resistive component of transmission impedance, and the transmission impedance is approximately purely inductive. It can reduce the coupling of active power and reactive power transmission. The positive inductive virtual impedance is introduced, the transmission impedance increases, and the current circulation caused by differences in the equivalent output impedance and connection line impedance is reduced. After introducing the virtual complex impedance, the voltage compensation is shown as follows:(10)uvd=Rvdi2d−ωLvdi2q,uvq=Rvqi2q+ωLvqi2d,where Rvq and Rvd are the virtual impedances in the dq rotation coordinates. To reduce the resistive component of the equivalent output impedance, the virtual impedance Rv takes a negative value, and Lvq and Lvd are the virtual inductive resistances in the dq coordinates [13].To increase the control precision and control the current circulation, the voltage-type VSG structure is adopted, and the terminal voltage control loop is added after the power control loop. The structure of the inner voltage and current control loop with virtual complex impedance is shown as Figure5.Figure 5
Voltage and current loop structure after adding virtual impedance.In Figure5, kuP and Kui are the proportional (P) and integral (I) gains of the voltage-loop PI controller, kiP is the proportional gain of the current inner loop, kPwm is the converter gain, Zvs is the virtual complex impedance (Zvs=Rv+Lvs), u∗ is the voltage-loop voltage command value, uv∗ is the voltage-loop reference value after introducing the virtual impedance, u is the output voltage, i1 is the converter-side current, and i2 is the grid-side output current. A mathematical model of the output voltage can be established according to Figure 5. When the virtual impedance is not introduced, it can be written as follows:(11)CsUs=U∗s−Uskup+kuis−I2s+UsCskipkpwm−Us−L1s+R1I2sL1s+R1.Equation (11) can be transformed into the following equation:(12)Us=kipkpwmkups+kuiU∗sL1Cs3+kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui,−L1s2+kipkpwm+R1sI2sL1Cs3+kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui.It is equivalent to the following equation:(13)Us=GsU∗s−ZosI2,where Gx is the equivalent closed-loop voltage gain and Z0s is the equivalent output impedance. These terms can be written as follows:(14)Gs=kipkpwmkups+kuiL1Cs3kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui,Zos=L1s2+kipkpwm+R1sL1Cs3+kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui.After adding the virtual impedance, the following equation can be obtained:(15)Uv∗s=U∗s−ZvsI2s.Substituting (16) into (13) gives(16)Us=GsU∗−ZvsI2−ZosI2,=GsU∗−GsZvs+ZosI2,where GsZvs+Z0s is the equivalent output impedance after adding the virtual complex impedance and Z0′s is shown as the following equation:(17)Z0′s=GsZvs+Zos.According to equation (15) and the parameters in Table 1, the equivalent output impedance can be illustrated using a Bode diagram without the added virtual impedance. Figure 6 is the equivalent output impedance Bode diagram for Kui  = 0, 100, 500, and 1000. It can be seen that when Kui = 0, the low-frequency equivalent output impedance is resistive, and it is the complex impedance characteristic at the fundamental frequency. When Kui=500, the equivalent output impedance is inductive at the fundamental frequency. The larger the value of Kui, the wider the frequency band that belongs to the inductive output impedance. After the introduction of the virtual inductive impedance Lv, the equivalent output impedance is inductive at low and fundamental frequencies. The fundamental output impedance is approximately equal to the virtual inductive impedance Lv, which weakens the effect of changes in the filter inductance parameters on the output impedance. Therefore, the equivalent output impedance at the fundamental frequency can be designed by introducing a virtual complex impedance to match the impedance of the parallel VSGs, and the aim of improving power distribution accuracy and current circulation suppression can be achieved.Table 1
VSG parameters.
Parameter/unitValueL1/mh5R1/Ω0.1C/uf10Zline/Ω0.5 + 0.06jLv/mh1kup0.3kui600kip2kpwm250Rv/Ω−0.4Figure 6
Equivalent output impedance Bode diagram with no virtual impedance.When selecting the virtual negative impedanceRv, the power decoupling and system stability requirements must be comprehensively considered. After introduction of the virtual impedance, the equivalent output impedance of the VSG at the fundamental frequency is approximately equal to the virtual impedance, so the transmission impedance at the fundamental frequency is given as follows:(18)Zs=Zvs+Zlines=Rline+Rv+sLline+Lv.The formulaRv=−Rline is established, and active power and reactive power transmission can achieve decoupling.The resistive component of transmission impedance has a damping effect on nonfundamental frequency disturbances. In the resistive component of transmission impedance setting, both the power decoupling and the nonfundamental frequency stability should be taken into consideration. The equivalent output impedance of the VSG after introducing the virtual impedance can be expressed as follows:(19)Zs=Zo′s+Zlines,=GsZvs+Zos+Zlines.The relationship between the output harmonic currentids and the voltage disturbance uds is shown as follows:(20)idsuds=1Zs=Ys.Therefore, the stability of the system can be analyzed by identifying the distribution of the characteristic roots of the characteristic equation. Combining (17), (18), and (20), the following equation can be obtained:(21)Ys=TsLlines+RlineTs+kipkpwmkups+kuiLvs+Rv+L1s2+kipkpwm+R1s,Ts=L1Cs3+kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui.WhenRline=0.5Ω and Lline=0.4mh, according to the parameters in Table 1, the root locus of Ys with respect to Rv is shown in Figure 7. The pole S1 determines the stability of YS. As Rv increases, S1 moves to the right of the complex plane and reaches the origin when Rv=−0.5Ω. Therefore, |Rv| should not be too large, it should be kept at a certain difference with Rline, so as S1 has a certain distance from the imaginary axis, and the system has a certain stability margin.Figure 7
Root locus ofYs for increasing Rv.


### 3.2. Reactive Power Loop Improvement

From the analysis in the previous section, it is known that the transmission of active power and reactive power is decoupled, when the transmission impedance of the VSG is inductive. Due to the integration link in the active power loop, the active power allocation of the parallel VSGs is related to the active power reference value and the active power droop coefficient, and it is not affected by the transmission impedance. The reactive power distribution is affected not only by the droop coefficient but also by the reactive power reference value and the transmission impedance. At the same time, the load voltage fluctuates because of the load fluctuation and the droop characteristics of the reactive power loop. The fluctuations of load voltage can be restrained, and the output voltage can be stabilized within a certain range by introducing a negative feedback linkKeE∗−E in the reactive power loop (here, Ke is the feedback coefficient, E∗ is the reference value of the load voltage, and E is the actual value of the load voltage). The integral link is introduced into the power loop, and it achieves the decoupling of reactive power and transmission impedance.The improved reactive power control loop structure is shown in Figure8. When the system is stable, the input of the integration link is 0. The expression of power and voltage is shown as follows:(22)DqiQi−Q∗=KeE∗−E.Figure 8
Improved virtual synchronous reactive power-voltage control block diagram.From equation (22), it can be determined that the reactive power distribution is related to the droop coefficient, reactive power reference value, and voltage feedback coefficient.It is independent of the transmission impedance. As long as the voltage feedback coefficients are equal and the reactive power reference value and droop coefficient are set in inverse proportion to the rated capacity, the ordered distribution of reactive power can be achieved [15, 16].


### 3.3. Structure of Optimization Control Strategy

A structure diagram of the inner voltage and current control loop is shown in Figure9. The terminal voltage is controlled in the synchronous rotating coordinate system.Figure 9
Structure of the inner voltage and current loop.The terminal voltage loop adopts the PI regulator, the current loop adopts the P regulator, and the decoupling compensation is added to the controller to eliminate the cross coupling of the dq-axis voltage and current, and it improves the system dynamic performance. In the figure, the subscripts “d” and “q” indicate the synchronous rotation dq coordinates, Ucdq and ildq are the terminal voltage and converter-side current in the synchronous rotation dq coordinate system, and ω is the virtual electrical angular frequency of the VSG [17, 18].The improved control strategy of parallel VSGs is illustrated in Figure10. In the figure, L1 is the output filter inductor, R1 is the sum of the internal resistance of the filter inductor and the internal resistance of the power device, C is the filter capacitor, and Lline and Rline are the connection line impedances. i1abc is the inductor current on the converter side, i2abc is the load current, and ucabc is the voltage of capacitor.Figure 10
Improved parallel virtual synchronization control strategy.The instantaneous powerP and Q calculation formulas are shown as equation (23), and the reference value U∗ of the midpoint voltage of the bridge arm is obtained by the improved power loop.(23)P=1.5ucdi2d+ucqi2q,Q=1.5ucdi2d−ucdi2q,where Ucd and Ucq are the d and q coordinate components of the filter capacitor voltage Uc, while i2d and i2q are the d and q coordinate components of the load current i2.Ignoring the role of the filter capacitorC, the relationship among the bridge midpoint voltage, terminal voltage, and output current of the converter side is given by equation (24). The relationships between voltage and current in the dq coordinate system are given by equations (25) and (26) [11].(24)ucabc=uabc−R1i1abc−L1di1abcdt,(25)ucducq=uduq−Y−1i1di1q,(26)Y=1R12+X12R1X1−X1R1,uduq=EcosφEsinφ,where i1d and i1q are the output currents of the converter side in the dq coordinate system, Ud and Uq are the dq-axis components obtained by dq decomposition of the midpoint voltage of the bridge arm using the terminal voltage vector orientation, Ucd and Ucq are the dq-axis components of the terminal voltage, Y is the filter impedance matrix, and X1 is the inductive reactance X1=ωL1. The phase angle φ represents the phase angle difference between the bridge midpoint voltage vector and the terminal voltage vector, which is equal to the integral of the difference between the virtual rotor angular velocity ω and the terminal voltage electrical angular velocity ωc of the virtual synchronous control [18].The reference value of the terminal voltage under the synchronous rotation coordinates is shown as equation (25). The terminal voltage is sent to inner the voltage current control loop to get the modulation voltage, the control signals of power components are obtained through sinusoidal pulse width modulation, and power equipartitioning and circulating current suppression of the parallel VSGs are realized.


## 3.1. Introduction of Virtual Impedance

From the analysis in the previous section, the active power and the reactive power are decoupled by setting the transmission impedance to be purely resistive or purely inductive, and this is the premise of independent control of active power and reactive power. Increasing the equivalent output impedance of the virtual synchronous generator can reduce the current circulation caused by the difference between the equivalent output impedance and the connection impedance. To change the transmission impedance characteristics of the VSG, the output current of the grid side is used as the feedback, and a virtual complex impedance consisting of resistance and inductive components is introduced. The negative resistive virtual impedance is used to reduce the resistive component of transmission impedance, and the transmission impedance is approximately purely inductive. It can reduce the coupling of active power and reactive power transmission. The positive inductive virtual impedance is introduced, the transmission impedance increases, and the current circulation caused by differences in the equivalent output impedance and connection line impedance is reduced. After introducing the virtual complex impedance, the voltage compensation is shown as follows:(10)uvd=Rvdi2d−ωLvdi2q,uvq=Rvqi2q+ωLvqi2d,where Rvq and Rvd are the virtual impedances in the dq rotation coordinates. To reduce the resistive component of the equivalent output impedance, the virtual impedance Rv takes a negative value, and Lvq and Lvd are the virtual inductive resistances in the dq coordinates [13].To increase the control precision and control the current circulation, the voltage-type VSG structure is adopted, and the terminal voltage control loop is added after the power control loop. The structure of the inner voltage and current control loop with virtual complex impedance is shown as Figure5.Figure 5
Voltage and current loop structure after adding virtual impedance.In Figure5, kuP and Kui are the proportional (P) and integral (I) gains of the voltage-loop PI controller, kiP is the proportional gain of the current inner loop, kPwm is the converter gain, Zvs is the virtual complex impedance (Zvs=Rv+Lvs), u∗ is the voltage-loop voltage command value, uv∗ is the voltage-loop reference value after introducing the virtual impedance, u is the output voltage, i1 is the converter-side current, and i2 is the grid-side output current. A mathematical model of the output voltage can be established according to Figure 5. When the virtual impedance is not introduced, it can be written as follows:(11)CsUs=U∗s−Uskup+kuis−I2s+UsCskipkpwm−Us−L1s+R1I2sL1s+R1.Equation (11) can be transformed into the following equation:(12)Us=kipkpwmkups+kuiU∗sL1Cs3+kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui,−L1s2+kipkpwm+R1sI2sL1Cs3+kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui.It is equivalent to the following equation:(13)Us=GsU∗s−ZosI2,where Gx is the equivalent closed-loop voltage gain and Z0s is the equivalent output impedance. These terms can be written as follows:(14)Gs=kipkpwmkups+kuiL1Cs3kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui,Zos=L1s2+kipkpwm+R1sL1Cs3+kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui.After adding the virtual impedance, the following equation can be obtained:(15)Uv∗s=U∗s−ZvsI2s.Substituting (16) into (13) gives(16)Us=GsU∗−ZvsI2−ZosI2,=GsU∗−GsZvs+ZosI2,where GsZvs+Z0s is the equivalent output impedance after adding the virtual complex impedance and Z0′s is shown as the following equation:(17)Z0′s=GsZvs+Zos.According to equation (15) and the parameters in Table 1, the equivalent output impedance can be illustrated using a Bode diagram without the added virtual impedance. Figure 6 is the equivalent output impedance Bode diagram for Kui  = 0, 100, 500, and 1000. It can be seen that when Kui = 0, the low-frequency equivalent output impedance is resistive, and it is the complex impedance characteristic at the fundamental frequency. When Kui=500, the equivalent output impedance is inductive at the fundamental frequency. The larger the value of Kui, the wider the frequency band that belongs to the inductive output impedance. After the introduction of the virtual inductive impedance Lv, the equivalent output impedance is inductive at low and fundamental frequencies. The fundamental output impedance is approximately equal to the virtual inductive impedance Lv, which weakens the effect of changes in the filter inductance parameters on the output impedance. Therefore, the equivalent output impedance at the fundamental frequency can be designed by introducing a virtual complex impedance to match the impedance of the parallel VSGs, and the aim of improving power distribution accuracy and current circulation suppression can be achieved.Table 1
VSG parameters.
Parameter/unitValueL1/mh5R1/Ω0.1C/uf10Zline/Ω0.5 + 0.06jLv/mh1kup0.3kui600kip2kpwm250Rv/Ω−0.4Figure 6
Equivalent output impedance Bode diagram with no virtual impedance.When selecting the virtual negative impedanceRv, the power decoupling and system stability requirements must be comprehensively considered. After introduction of the virtual impedance, the equivalent output impedance of the VSG at the fundamental frequency is approximately equal to the virtual impedance, so the transmission impedance at the fundamental frequency is given as follows:(18)Zs=Zvs+Zlines=Rline+Rv+sLline+Lv.The formulaRv=−Rline is established, and active power and reactive power transmission can achieve decoupling.The resistive component of transmission impedance has a damping effect on nonfundamental frequency disturbances. In the resistive component of transmission impedance setting, both the power decoupling and the nonfundamental frequency stability should be taken into consideration. The equivalent output impedance of the VSG after introducing the virtual impedance can be expressed as follows:(19)Zs=Zo′s+Zlines,=GsZvs+Zos+Zlines.The relationship between the output harmonic currentids and the voltage disturbance uds is shown as follows:(20)idsuds=1Zs=Ys.Therefore, the stability of the system can be analyzed by identifying the distribution of the characteristic roots of the characteristic equation. Combining (17), (18), and (20), the following equation can be obtained:(21)Ys=TsLlines+RlineTs+kipkpwmkups+kuiLvs+Rv+L1s2+kipkpwm+R1s,Ts=L1Cs3+kipkpwm+R1Cs2+kipkpwmkup+1s+kipkpwmkui.WhenRline=0.5Ω and Lline=0.4mh, according to the parameters in Table 1, the root locus of Ys with respect to Rv is shown in Figure 7. The pole S1 determines the stability of YS. As Rv increases, S1 moves to the right of the complex plane and reaches the origin when Rv=−0.5Ω. Therefore, |Rv| should not be too large, it should be kept at a certain difference with Rline, so as S1 has a certain distance from the imaginary axis, and the system has a certain stability margin.Figure 7
Root locus ofYs for increasing Rv.


## 3.2. Reactive Power Loop Improvement

From the analysis in the previous section, it is known that the transmission of active power and reactive power is decoupled, when the transmission impedance of the VSG is inductive. Due to the integration link in the active power loop, the active power allocation of the parallel VSGs is related to the active power reference value and the active power droop coefficient, and it is not affected by the transmission impedance. The reactive power distribution is affected not only by the droop coefficient but also by the reactive power reference value and the transmission impedance. At the same time, the load voltage fluctuates because of the load fluctuation and the droop characteristics of the reactive power loop. The fluctuations of load voltage can be restrained, and the output voltage can be stabilized within a certain range by introducing a negative feedback linkKeE∗−E in the reactive power loop (here, Ke is the feedback coefficient, E∗ is the reference value of the load voltage, and E is the actual value of the load voltage). The integral link is introduced into the power loop, and it achieves the decoupling of reactive power and transmission impedance.The improved reactive power control loop structure is shown in Figure8. When the system is stable, the input of the integration link is 0. The expression of power and voltage is shown as follows:(22)DqiQi−Q∗=KeE∗−E.Figure 8
Improved virtual synchronous reactive power-voltage control block diagram.From equation (22), it can be determined that the reactive power distribution is related to the droop coefficient, reactive power reference value, and voltage feedback coefficient.It is independent of the transmission impedance. As long as the voltage feedback coefficients are equal and the reactive power reference value and droop coefficient are set in inverse proportion to the rated capacity, the ordered distribution of reactive power can be achieved [15, 16].


## 3.3. Structure of Optimization Control Strategy

A structure diagram of the inner voltage and current control loop is shown in Figure9. The terminal voltage is controlled in the synchronous rotating coordinate system.Figure 9
Structure of the inner voltage and current loop.The terminal voltage loop adopts the PI regulator, the current loop adopts the P regulator, and the decoupling compensation is added to the controller to eliminate the cross coupling of the dq-axis voltage and current, and it improves the system dynamic performance. In the figure, the subscripts “d” and “q” indicate the synchronous rotation dq coordinates, Ucdq and ildq are the terminal voltage and converter-side current in the synchronous rotation dq coordinate system, and ω is the virtual electrical angular frequency of the VSG [17, 18].The improved control strategy of parallel VSGs is illustrated in Figure10. In the figure, L1 is the output filter inductor, R1 is the sum of the internal resistance of the filter inductor and the internal resistance of the power device, C is the filter capacitor, and Lline and Rline are the connection line impedances. i1abc is the inductor current on the converter side, i2abc is the load current, and ucabc is the voltage of capacitor.Figure 10
Improved parallel virtual synchronization control strategy.The instantaneous powerP and Q calculation formulas are shown as equation (23), and the reference value U∗ of the midpoint voltage of the bridge arm is obtained by the improved power loop.(23)P=1.5ucdi2d+ucqi2q,Q=1.5ucdi2d−ucdi2q,where Ucd and Ucq are the d and q coordinate components of the filter capacitor voltage Uc, while i2d and i2q are the d and q coordinate components of the load current i2.Ignoring the role of the filter capacitorC, the relationship among the bridge midpoint voltage, terminal voltage, and output current of the converter side is given by equation (24). The relationships between voltage and current in the dq coordinate system are given by equations (25) and (26) [11].(24)ucabc=uabc−R1i1abc−L1di1abcdt,(25)ucducq=uduq−Y−1i1di1q,(26)Y=1R12+X12R1X1−X1R1,uduq=EcosφEsinφ,where i1d and i1q are the output currents of the converter side in the dq coordinate system, Ud and Uq are the dq-axis components obtained by dq decomposition of the midpoint voltage of the bridge arm using the terminal voltage vector orientation, Ucd and Ucq are the dq-axis components of the terminal voltage, Y is the filter impedance matrix, and X1 is the inductive reactance X1=ωL1. The phase angle φ represents the phase angle difference between the bridge midpoint voltage vector and the terminal voltage vector, which is equal to the integral of the difference between the virtual rotor angular velocity ω and the terminal voltage electrical angular velocity ωc of the virtual synchronous control [18].The reference value of the terminal voltage under the synchronous rotation coordinates is shown as equation (25). The terminal voltage is sent to inner the voltage current control loop to get the modulation voltage, the control signals of power components are obtained through sinusoidal pulse width modulation, and power equipartitioning and circulating current suppression of the parallel VSGs are realized.


## 4. Simulation and Experiment

The model of two parallel VSGs is constructed in the Matlab/Simulink environment to verify the effectiveness of improved control strategy. The main parameters are listed in Table2. The total simulation time of the system is 2 seconds. From the start to 0.6 s, the first VSG runs alone, the other VSG is connected at 0.6 s, and the first VSG exits at 1.2 s.Table 2
System simulation parameters.
Parameter/unitValueUdc/V600U/V380J, Dp0.05, 15Dq480P∗/kW9Q∗/kvar6Lv/mh1P0/kW4Q0Q0/kvar3L1/mh, R1/Ω5, 0.1L2/mh, R2/Ω4.8, 0.1Zline1/Ω0.45 + 0.06jZline2/Ω0.55 + 0.07jRv/Ω−0.2Figure11 shows the simulation results of the voltage, current, and circulating current when two parallel VSGs are running with the traditional control strategy and the improved control strategy.Figure 11
Output voltage, current, and circulating current of two parallel VSGs: (a) traditional virtual synchronization control strategy; (b) improved virtual synchronization control strategy.
(a)(b)Figure12 shows the output power results of two VSGs with optimization control strategy.Figure 12
Output active power and reactive power of the parallel virtual synchronous control strategy: (a) traditional virtual synchronization control strategy; (b) improved virtual synchronization control strategy.
(a)(b)Due to the influence of impedance differences, the circulating current of parallel VSGs with traditional control strategy is relatively large, and the reactive power equipartitioning cannot be achieved.The dynamic current dividing response time of the improved control strategy is shorter, the amplitude of the circulating current is smaller in the steady state, and the effect of reactive power distribution is significantly improved.The experimental platform of two parallel VSGs is shown in Figure13, and the main parameters of the experiment are presented in Table 3. In the parallel VSGs experiment, the public load is 3 kW. During the experiment, an inverter works alone at the beginning, and another VSG is connected after a while.Figure 13
Pictures of the test platform.Table 3
Experimental parameters.
Parameter/unitValueUdc/V600kup0.07kui5Dqi190Zline1/Ω0.4 + 0.06jZline2/Ω0.5 + 0.07jL1/mh, R1/Ω5, 0.1L2/mh, R2/Ω4.8, 0.1Lv/mh1U/V380J, Dp0.02, 6kip2Ke0.2Q0/kvar1P0/kW3P∗/kW4Q∗/kvar1Rv/Ω−0.2Figure14 shows the output current and circulating current waveforms of two parallel VSGs with the traditional control strategy and the improved control strategy. Compared with the traditional control strategy, the improved control strategy can respond more quickly, and the circulating current is smaller.Figure 14
Output current and circulating current of two parallel VSGs: (a) traditional virtual synchronization control strategy; (b) improved virtual synchronization control strategy.
(a)(b)


## 5. Conclusions

In this paper, the problems of power equipartitioning and circulating current suppression in the parallel operation mode of VSGs were studied. The conclusions and requirements for future work are as follows.(1)
When multiple VSGs are connected in parallel, the resistance and inductance characteristics of the transmission impedance cause the coupling between reactive power and active power. The reactive power distribution is affected by the droop coefficient, the reactive power reference value, and the transmission impedance. The current circulation of the parallel VSGs is affected by the deviation of output voltage and the value of transmission impedance.(2)
The virtual complex impedance including a resistive component and an inductive component is introduced into the inner voltage and current control loop. It reduces the coupling between active power and reactive power, and it reduces the disorderly distribution of power due to impedance differences. The load voltage control loop and integral part are introduced in the reactive power link of the VSG, which reduce the deviation of output voltage and restrain circulating current.(3)
According to the comparative analysis of the influence of control parameters and virtual complex impedance on the equivalent output impedance, the control parameters and virtual complex impedance are selected.



---

*Source: 8971213-2020-10-19.xml*