Adaptive Virtual Impedance Droop Control Based on Consensus Control of Reactive Current

: It is difﬁcult to achieve accurate distribution of reactive power based on conventional droop control due to the line impedance mismatch in an islanded microgrid. An adaptive virtual impendence method based on consensus control of reactive current is proposed in this paper. A distributed control structure without the central controller has been established. In this structure, each distributed generation unit (DG) is an independent agent, one-way communication is used between the adjacent DGs, and the reactive power sharing is equivalent to a problem of reactive power current consensus. It has been proven that the system is asymptotically stable under the proposed control strategy. When the adjacent DG’s reactive power is not proportionally distributed, the current weight error term will generate a virtual impedance correction term through the proportional-integral controller based on the reactive current consensus control strategy, thus introducing adaptive virtual impedance to eliminate mismatches in output impedance between DGs. Reactive power auto-proportional distribution can be achieved without knowing the line impedance. At the same time, the power control loop is simpliﬁed and the virtual impedance compensation angle is employed to compensate the decreased reference voltage magnitude and varied phase angle due to the introduction of the virtual impedance, so the stability of the system can be improved. Finally, the correctness and effectiveness of the proposed strategy are veriﬁed by modeling analysis and microgrid simulations.


Introduction
The increasing demand for renewable energy-based power generation systems consisting of inverters necessitates the design of sophisticated inverter control systems [1][2][3]. Apart from ensuring the quality of power output, the inverter control system should be able to coordinate the operation of multiple distributed generators (DGs) and distribute power to loads according to the capacity of each power source [4]. Research on energy management, reliability of component, coordination, and optimization of renewable energy in hybrid microgrids has been conducted [5][6][7][8][9]. The droop control system, based on the droop characteristic of synchronous generators, offers the benefits of sparse communication and 'plug and play'. Therefore, it is widely used in microgrid control [10][11][12]. While the P-f droop control enables accurate sharing of active power during real-world microgrid operation, the Q-V droop control exhibits poor performance with regards to reactive power sharing because it is susceptible to the influence of line impedance.
where 0 E and 0 ω are no-load voltage and frequency, respectively, i m and i n are the active and reactive droop coefficients of inverter I, respectively.
In [26], the virtual impedance is realized by adjusting the parameters of the double-loop PI controller. The output voltage is given by [26]: where ( ) G s is voltage gain, and ( ) o Z s is the equivalent output impedance which can be calculated using Equation (8) It can be seen from Equation (8) that the equivalent output impedance is not only related to the filter parameters L and C, but also influenced by the PI controller. Using the design parameters stated in [26] (see Table 1), the equivalent output impedance is inductive, and the influence of the line impedance can be ignored. The output impedance Bode diagram is shown in Figure 2.      It can be seen by analyzing Figure 2, while the output impedance of the inverter can be inductive at 50 Hz, the output impedance can only be about 0.01 mΩ, which limits the range of the line impedance parameters. That is, only when a relatively small line impedance is chosen, the dominant status of the output impedance can be ensured. Moreover, changing the double-loop PI controller parameters can only affect small adjustments to the output impedance, and if consideration is solely given to the design of the output impedance when choosing PI control parameters, the system stability margin will inevitably be affected negatively, resulting in reduced stability.
In this study, a new virtual impedance , V i Z using direct series connection is designed. Thus, the system output impedance is given by: where i Z is the equivalent output impedance of the micro-source, The static virtual impedance always ensures the system impedance is mainly inductive. The adaptive virtual impedance automatically adjusts the virtual impedance magnitude according to the equivalent output impedance of the adjacent DG. Hence, accurate reactive power sharing is achieved by eliminating the output impedance mismatch. This method does not rely on the adjustment of the double-loop control parameters of the system and, therefore, allows for a greater range of impedance adjustment compared with the previous method.

Reactive Power Sharing Based on Reactive Current Control
In [35], it is noted that in the typical "power-voltage-current" closed-loop droop control, the power control loop increases the complexity of the system and affects its dynamic characteristics due to the power calculations required. A method to realize droop control by controlling active and reactive currents was proposed in [36]. The key merit of this control method is that the system is incapable of yielding infinitely impulse current when short circuit occurs ( 0 i E = ). This paper applies this method to the design of adaptive virtual impedance.
Active and reactive currents ai I and ri I of equivalent impedance Z can be calculated as: It can be seen by analyzing Figure 2, while the output impedance of the inverter can be inductive at 50 Hz, the output impedance can only be about 0.01 mΩ, which limits the range of the line impedance parameters. That is, only when a relatively small line impedance is chosen, the dominant status of the output impedance can be ensured. Moreover, changing the double-loop PI controller parameters can only affect small adjustments to the output impedance, and if consideration is solely given to the design of the output impedance when choosing PI control parameters, the system stability margin will inevitably be affected negatively, resulting in reduced stability.
In this study, a new virtual impedance Z V,i using direct series connection is designed. Thus, the system output impedance is given by: where Z i is the equivalent output impedance of the micro-source, Z line,i is the line impedance and Z V,i is the virtual impedance consisting of the static virtual impedance Z * V,i and the adaptive virtual impedance Z adp,i .The static virtual impedance always ensures the system impedance is mainly inductive. The adaptive virtual impedance automatically adjusts the virtual impedance magnitude according to the equivalent output impedance of the adjacent DG. Hence, accurate reactive power sharing is achieved by eliminating the output impedance mismatch. This method does not rely on the adjustment of the double-loop control parameters of the system and, therefore, allows for a greater range of impedance adjustment compared with the previous method.

Reactive Power Sharing Based on Reactive Current Control
In [35], it is noted that in the typical "power-voltage-current" closed-loop droop control, the power control loop increases the complexity of the system and affects its dynamic characteristics due to the power calculations required. A method to realize droop control by controlling active and reactive currents was proposed in [36]. The key merit of this control method is that the system is incapable of yielding infinitely impulse current when short circuit occurs (E i = 0). This paper applies this method to the design of adaptive virtual impedance.
Active and reactive currents I ai and I ri of equivalent impedance Z can be calculated as: Considering that (3) and (4) can be rewritten as: The active current and reactive current droop characteristics are defined as: In order to satisfy the proportional distribution of the load reactive power, the droop coefficient and reactive power should be designed to be inversely proportional [30]: When using reactive current control, the following equation should be satisfied: Substituting (13) into Equation (15), the reactive current flow from the DGi to load is obtained as: To satisfy Equation (17), the following equation can be derived from Equation (18): As it can be seen in Equation (19) that in order to achieve accurate reactive power distribution, k ni should be proportional to the line reactance Xi. Considering Equation (17), the line reactances must be designed to be inversely proportional to the reactive current, respectively: where X i ∼ = Z i as previously described.

Microgrid Communication Topology
If each DG in the microgrid is regarded as a node, and the communication links between the DGs are equated to edges, the topological model of the entire microgrid can be represented by a digraph (directed graph). Usually, a digraph G refers to an ordered triplet (V(G), E(G), A(G)), where V(G) = (v 1 , v 2 , · · · v N ) is a set of non-empty nodes, E(G) ⊂ V(G) × V(G) is an edge set that edges of points with communication associated and A(G) is a correlation function which is typically an associated adjacency matrix A(G) = [a ij ] ∈ N×N . In this paper, the digraph is assumed to be time-invariant, i.e., A(G) is a constant. The edge from node j to node i is denoted as (v j , v i ).
, then a ij = 1 and node j is called a neighbor of node i. The set of all nodes that satisfy the above condition is denoted as Many properties of the digraph are manifested in the above matrix. In this paper, adjacent nodes in the DG network are treated as follows: Node i can obtain information from node j but node j cannot obtain information from node i, with both the indegree and outdegree of each node set to 1, as shown in Figure 3. Therefore, the communication between the DGs is regarded as a one-way distributed loose communication. Moreover, the network topology in Figure 3 contains a directional spanning tree. When a communication line is faulty, the system still maintains high reliability and the sparse network structure reduces the communication cost.
The Laplacian matrix is defined as properties of the digraph are manifested in the above matrix. In this paper, adjacent nodes in the DG network are treated as follows: Node i can obtain information from node j but node j cannot obtain information from node i, with both the indegree and outdegree of each node set to 1, as shown in Figure 3. Therefore, the communication between the DGs is regarded as a one-way distributed loose communication. Moreover, the network topology in Figure 3 contains a directional spanning tree. When a communication line is faulty, the system still maintains high reliability and the sparse network structure reduces the communication cost.

Consensus Control
According to the consensus control strategy, in order to control the reactive current, the ni ri k I term of each DG should be assigned the same value as shown in Equation (17). Thus, the reactive power distribution based on consensus control of the reactive current can be transformed into the synchronization problem of the first-order linear multi-agent system [37]. Linearization of Equation (17) can be formulated as: According to the consensus control theory, ri I u is an adjoint control term. If the term is not 0, there is a reactive power deviation between the local and its neighbors: where aij is an adjacency matrix element that reflects the connection of DGs. Then the matrix form of the consensus control system can be expressed as: where the above variables are defined as . The stability of the proposed control system is proven below. In this paper a simple quadratic Lyapunov function is used to show the stability of the system. However, [38] indicates that for some applications using a non-quadratic Lyapunov function can also lead to a better performance.
Constructing Lyapunov energy function as:

Consensus Control
According to the consensus control strategy, in order to control the reactive current, the k ni I ri term of each DG should be assigned the same value as shown in Equation (17). Thus, the reactive power distribution based on consensus control of the reactive current can be transformed into the synchronization problem of the first-order linear multi-agent system [37]. Linearization of Equation (17) can be formulated as: k n1 .
According to the consensus control theory, u I ri is an adjoint control term. If the term is not 0, there is a reactive power deviation between the local and its neighbors: u I ri = −C nI r e k ni I ri (22) where C nI r is the coupling gain and e k ni I ri is the reactive current weight error between local and adjacent DG units: where a ij is an adjacency matrix element that reflects the connection of DGs. Then the matrix form of the consensus control system can be expressed as: k n . I r = u I r (24) u I r = −C nI r e k n I r (25) e k n I r = Lk n I r (26) where the above variables are defined as k n I rN ] T , u I r = [u I r1 , u I r2 , · · · u I rN ] T , e k n I r = [e k n1 I r1 , e k n1 I r1 , · · · e k nN I rN ] T , k n I r = [k n1 I r1 , k n2 I r2 , · · · k nN I rN ] T . The stability of the proposed control system is proven below. In this paper a simple quadratic Lyapunov function is used to show the stability of the system. However, [38] indicates that for some applications using a non-quadratic Lyapunov function can also lead to a better performance.
Constructing Lyapunov energy function as: where P = P T , P is a positive definite matrix. The following equation can be obtained from Equation (27): . V k n I r = e T k n I r P . e k n I r = e T k n I r PLk n . I r (28) Combining the equations above, the following equation can be obtained: . V k n I r = e T k n I r P . e k n I r = e T k n I r PLu I r (29) Defining H ≡ L, Substituting Equation (25) into (29), and the following equation is described as: This paper considers the configuration of four different rated power DG units. Letting P = diag{1, 1, 1, 1}, then the Laplacian of the communication topology is: Then: The matrix PH + H T P is positive definite, so . V k n I r < 0, using Lyapunov's second law, the system is determined to be asymptotically stable.

Method of Realizing Virtual Impedance
From Figure 4, the following equation can be obtained: where Z V,i = R V,i + jω i L V,i is the equivalent virtual impedance, k L and k R are the proportional gains (which are used to adjust the value of the adaptive virtual impedance), L * V,i represents the static virtual inductance and the reactive current error term e k ni I ri is processed by the proportional integral controller (PI) to yield the virtual impedance correction term δI r,i .
When the reactive power output from DGs are out of proportion, the current error term will be generated according to Equation (23). After the error term is processed by the proportional integral controller, the virtual impedance correction term will be generated and the adaptive virtual impedance of the system will be adjusted accordingly. In this paper, Section 3 indicates that the magnitude of the system impedance should be adjusted so that it is inversely proportional to the reactive load shared.
where , , , is the equivalent virtual impedance, L k and R k are the proportional gains When the reactive power output from DGs are out of proportion, the current error term will be generated according to Equation (23). After the error term is processed by the proportional integral As previously discussed, since the equivalent output impedance Z o (s) of the micro-source mentioned in [26] is only 0.01 mΩ, the adjustment range is limited. In this paper, the static virtual impedance Z * V (s) is connected in series and Z * V (s) >> Z o (s). In this case, the equivalent virtual impedance of the micro-source becomes the dominant parameter. Thus, by adjusting the value of the adaptive virtual impedance Z adp (s), the magnitude of the output impedance can be changed to achieve proportional sharing of reactive power.
After the introduction of the virtual impedance, the equivalent circuit of the DG is shown in Figure 5. controller, the virtual impedance correction term will be generated and the adaptive virtual impedance of the system will be adjusted accordingly. In this paper, Section 3 indicates that the magnitude of the system impedance should be adjusted so that it is inversely proportional to the reactive load shared. As previously discussed, since the equivalent output impedance After the introduction of the virtual impedance, the equivalent circuit of the DG is shown in Figure 5.  Table 2 and the Bode plot of the micro-source equivalent output impedance is shown in Figure 6.   Taking the virtual impedance Z V (s) = Z adp (s) + Z * V (s) to be inductive, the expression for the micro-source equivalent virtual impedance is shown in Equation (35), where the gain and phase angle of the transfer function G(s) at the fundamental frequency are 1 and 0, respectively. The droop control parameters are shown in Table 2 and the Bode plot of the micro-source equivalent output impedance is shown in Figure 6.
It can be determined by analyzing Figure 6 that the design method proposed in this paper ensures that the virtual impedance is inductive when the fundamental frequency is 50 Hz. Meanwhile, the magnitude of the virtual impedance reaches 0.252 Ω. Thus, a larger impedance adjustment range is available, compared with that offered by the design method proposed in [26]. parameters are shown in Table 2 and the Bode plot of the micro-source equivalent output impedance is shown in Figure 6.   Figure 6. Bode diagram of output impendence. Figure 6. Bode diagram of output impendence.

Reference Voltage Phase Shift Compensation Using Virtual Impedance Angle
The voltage drop V vi_d and V vi_q caused by the introduction of virtual impedance in synchronous reference frame can be expressed as (36) and (37) where I di and I qi are active and reactive currents respectively. The new voltage reference is calculated by the following equation where V dref and V qref are the initial reference voltage signals, so the change in control loop voltage in the dq-coordinate system as a result of introducing virtual impedance is depicted in Figure 7: Energies 2018, 11, x FOR PEER REVIEW 9 of 16 It can be determined by analyzing Figure 6 that the design method proposed in this paper ensures that the virtual impedance is inductive when the fundamental frequency is 50 Hz. Meanwhile, the magnitude of the virtual impedance reaches 0.252 Ω. Thus, a larger impedance adjustment range is available, compared with that offered by the design method proposed in [26].

Reference Voltage Phase Shift Compensation Using Virtual Impedance Angle
The voltage drop where di I and qi I are active and reactive currents respectively. The new voltage reference is calculated by the following equation where dref V and qref V are the initial reference voltage signals, so the change in control loop voltage in the dq-coordinate system as a result of introducing virtual impedance is depicted in Figure 7: After developing Equation (40) v θ can be calculated using Equation (41) from reference [39].  Figure 7 one can calculate the virtual angle added by the virtual impendence, which is between V dqnewref and V dqref . Equation (40) describes the relation between V dqnewref and V dqref .
The change in the angle at different virtual impedances is shown in Figure 8. It can be determined from Figures 7 and 8 that the introduction of the virtual impedance changes not only the amplitude of the reference voltage, but also its phase angle. In order to compensate for the deviation in the reference voltage phase angle caused by the adaptive virtual impedance, a virtual impedance angle is proposed, as shown in Figure 9: Figure 9. Voltage reference generator block.
The phase shift due to the introduction of the virtual impedance is eliminated by introducing v θ into the voltage reference loop. Note that i ω is the angular frequency of DG, obtained by the active power control loop, while ref V is provided directly. As discussed in Section 4.2, the virtual impedance correction term is obtained through consensus control of the reactive current to eliminate the output impedance mismatch and, thus, achieve proportional sharing of reactive power. There is no need to obtain reactive power information throughout the whole process. Therefore, this paper has simplified the power control by removing the reactive power control loop and keeping the original active power control loop. Reference [35] points out that, after this streamlining, the voltage droop characteristic is still maintained. The overall control topology is depicted in Figure 10. It can be determined from Figures 7 and 8 that the introduction of the virtual impedance changes not only the amplitude of the reference voltage, but also its phase angle. In order to compensate for the deviation in the reference voltage phase angle caused by the adaptive virtual impedance, a virtual impedance angle is proposed, as shown in Figure 9: It can be determined from Figures 7 and 8 that the introduction of the virtual impedance changes not only the amplitude of the reference voltage, but also its phase angle. In order to compensate for the deviation in the reference voltage phase angle caused by the adaptive virtual impedance, a virtual impedance angle is proposed, as shown in Figure 9: Figure 9. Voltage reference generator block.
The phase shift due to the introduction of the virtual impedance is eliminated by introducing v θ into the voltage reference loop. Note that i ω is the angular frequency of DG, obtained by the active power control loop, while ref V is provided directly. As discussed in Section 4.2, the virtual impedance correction term is obtained through consensus control of the reactive current to eliminate the output impedance mismatch and, thus, achieve proportional sharing of reactive power. There is no need to obtain reactive power information throughout the whole process. Therefore, this paper has simplified the power control by removing the reactive power control loop and keeping the original active power control loop. Reference [35] points out that, after this streamlining, the voltage droop characteristic is still maintained. The overall control topology is depicted in Figure 10. The phase shift due to the introduction of the virtual impedance is eliminated by introducing θ v into the voltage reference loop. Note that ω i is the angular frequency of DG, obtained by the active power control loop, while V ref is provided directly.
As discussed in Section 4.2, the virtual impedance correction term is obtained through consensus control of the reactive current to eliminate the output impedance mismatch and, thus, achieve proportional sharing of reactive power. There is no need to obtain reactive power information throughout the whole process. Therefore, this paper has simplified the power control by removing the reactive power control loop and keeping the original active power control loop. Reference [35] points out that, after this streamlining, the voltage droop characteristic is still maintained. The overall control topology is depicted in Figure 10.
control of the reactive current to eliminate the output impedance mismatch and, thus, achieve proportional sharing of reactive power. There is no need to obtain reactive power information throughout the whole process. Therefore, this paper has simplified the power control by removing the reactive power control loop and keeping the original active power control loop. Reference [35] points out that, after this streamlining, the voltage droop characteristic is still maintained. The overall control topology is depicted in Figure 10. . Block diagram of the entire system control scheme. Figure 10. Block diagram of the entire system control scheme.

Simulation Results
The microgrid shown in Figure 11 was modelled in MATLAB R2014b/Simulink (MATLAB R2014b, MathWorks, Natick, MA, USA) to verify the effectiveness of the proposed control method. The microgrid contains four DGs and the line impedance between the DGs is set to resistive-inductive. The specific DG parameters, voltage and current double loop controller parameters, load parameters, and line impedance parameters are shown in Table 3.

Simulation Results
The microgrid shown in Figure 11 was modelled in MATLAB R2014b/Simulink (MATLAB R2014b, MathWorks, Natick, MA, USA) to verify the effectiveness of the proposed control method. The microgrid contains four DGs and the line impedance between the DGs is set to resistiveinductive. The specific DG parameters, voltage and current double loop controller parameters, load parameters, and line impedance parameters are shown in Table 3.     The simulation was designed to verify the effectiveness of the proposed control strategy with respect to reactive power sharing and voltage control. Four inverters with rated capacities in the ratio of 1.5:3:4:4 were used in the simulation. Between 0 and 1 s of the simulation period, each system operates with a local load. After 1 s, the loads were added. Figure 12 shows the simulation results of the conventional droop control strategy and the control strategy based on consensus control of the reactive current. It can be found that between 1 and 2 s the active power is proportionally shared under both control strategies (as shown in Figure 12a,b). However, due to the mismatch of the line impedance, the reactive power sharing of the conventional droop control strategy is relatively poor (as shown in Figure 12c). In contrast, the proposed control strategy allows the equivalent output impedance of each micro-source to be effectively matched, thereby enabling the proportional sharing of reactive power in the ratio of 1.5:3:4:4 (as shown in Figure 12d). The system frequency will eventually converge and the frequency waveform obtained by the control strategy proposed in this paper is smoother (as shown in Figure 12e-f). In the simulation, the load switching is performed in 1 s. During the switching process, the virtual impedance will automatically change with the load distribution (as shown in Figure 12g-h). It can be seen, finally, that the virtual reactance dominates the system, which satisfies the P-f, Q-V droop control requirements.  Figure 12d). The system frequency will eventually converge and the frequency waveform obtained by the control strategy proposed in this paper is smoother (as shown in Figure 12e-f). In the simulation, the load switching is performed in 1 s. During the switching process, the virtual impedance will automatically change with the load distribution (as shown in Figure 12g-h). It can be seen, finally, that the virtual reactance dominates the system, which satisfies the P-f, Q-V droop control requirements. To verify the effectiveness of introducing the virtual impedance compensation angle θv, the output power of the first inverter is taken as an example, as shown in Figure 13. Since the virtual impedance compensation angle eliminates the phase angle difference of the reference voltage, the power waveform is smoother and more stable compared with the situation when the compensation angle is not introduced.  To verify the effectiveness of introducing the virtual impedance compensation angle θ v , the output power of the first inverter is taken as an example, as shown in Figure 13. Since the virtual impedance compensation angle eliminates the phase angle difference of the reference voltage, the power waveform is smoother and more stable compared with the situation when the compensation angle is not introduced. To verify the effectiveness of introducing the virtual impedance compensation angle θv, the output power of the first inverter is taken as an example, as shown in Figure 13. Since the virtual impedance compensation angle eliminates the phase angle difference of the reference voltage, the power waveform is smoother and more stable compared with the situation when the compensation angle is not introduced. The voltage sag at the microgrid AC bus is shown in Figure 14. As the conventional droop control is deviation-based, voltage sags may occur, especially when a virtual impedance is added. Compared with the control method based on direct introduction of virtual impedance, the bus voltage sag caused by the adaptive virtual impedance control method proposed in this paper is low.  The voltage sag at the microgrid AC bus is shown in Figure 14. As the conventional droop control is deviation-based, voltage sags may occur, especially when a virtual impedance is added. Compared with the control method based on direct introduction of virtual impedance, the bus voltage sag caused by the adaptive virtual impedance control method proposed in this paper is low. To verify the effectiveness of introducing the virtual impedance compensation angle θv, the output power of the first inverter is taken as an example, as shown in Figure 13. Since the virtual impedance compensation angle eliminates the phase angle difference of the reference voltage, the power waveform is smoother and more stable compared with the situation when the compensation angle is not introduced. The voltage sag at the microgrid AC bus is shown in Figure 14. As the conventional droop control is deviation-based, voltage sags may occur, especially when a virtual impedance is added. Compared with the control method based on direct introduction of virtual impedance, the bus voltage sag caused by the adaptive virtual impedance control method proposed in this paper is low.  A simulation was carried out to verify the reliability of the system under the reactive current consensus control strategy when communication failure occurs. In the simulation, the communication link between DGs 1 and 2 fails. Despite the failed communication link, the system is able to maintain the ratio of 1.5:3:4:4 in the sharing of reactive power while active power is controlled (as shown in Figure 15a,b). Therefore, when one communication link fails, the proposed control strategy still outperforms the conventional droop control method with regards to reactive power sharing. Even when communication fails, the adaptive virtual impedance based on reactive current consistency proposed in this paper still automatically changes with the load distribution during the load switching process (as shown in Figure 15c,d) and, ultimately, virtual reactance remains a dominant role in the system.
(as shown in Figure 15a,b). Therefore, when one communication link fails, the proposed control strategy still outperforms the conventional droop control method with regards to reactive power sharing. Even when communication fails, the adaptive virtual impedance based on reactive current consistency proposed in this paper still automatically changes with the load distribution during the load switching process (as shown in Figure 15c,d) and, ultimately, virtual reactance remains a dominant role in the system.

Conclusions
In this paper, reactive power sharing within an isolated microgrid using virtual impedance has been investigated.
The microgrid system adopts a distributed control structure without a central controller, in which each DG acts as an independent agent and adjacent DGs engage in one-way communication. The reactive power sharing is achieved through consensus control of the reactive current. In addition, the microgrid system proves to be asymptotically stable under this control strategy.
The virtual impedance correction term is obtained by the consensus control strategy of reactive current. Together with the static virtual impedance, the adaptive virtual impedance is introduced to eliminate the output impedance mismatch between adjacent micro-sources, thereby achieving proportional reactive power sharing without knowing specific line impedance parameters. Moreover, it has been proved that the proposed method can perform well even when a communication link fails.
The virtual impedance compensation angle θv is introduced to compensate for phase angle offset in the micro-source reference voltage caused by the introduction of virtual impedance, thus, the

Conclusions
In this paper, reactive power sharing within an isolated microgrid using virtual impedance has been investigated.
The microgrid system adopts a distributed control structure without a central controller, in which each DG acts as an independent agent and adjacent DGs engage in one-way communication. The reactive power sharing is achieved through consensus control of the reactive current. In addition, the microgrid system proves to be asymptotically stable under this control strategy.
The virtual impedance correction term is obtained by the consensus control strategy of reactive current. Together with the static virtual impedance, the adaptive virtual impedance is introduced to eliminate the output impedance mismatch between adjacent micro-sources, thereby achieving proportional reactive power sharing without knowing specific line impedance parameters. Moreover, it has been proved that the proposed method can perform well even when a communication link fails.
The virtual impedance compensation angle θ v is introduced to compensate for phase angle offset in the micro-source reference voltage caused by the introduction of virtual impedance, thus, the stability of the system can be enhanced. The power control loop is simplified without affecting the original droop characteristics, thereby reducing system complexity and limiting the voltage sag.
Author Contributions: The paper was a collaborative effort among the authors. Z.L. and Q.W. carried out relevant theoretical research, performed the simulation, analyzed the data, and wrote the paper. Y.Z., J.Z. and E.M. provided critical comments.