Improved Droop Control Strategy of Multiple Energy Storage Applications in an AC Microgrid Based on the State of Charge

Distributed energy storage technology is used to stabilize the frequency and voltage of the microgrid operating in islanded mode. However, due to the inconsistent state of charge (SoC) of the energy storage unit (ESU), the active power output of the ESU cannot be shared reasonably. On the basis of stabilizing voltage and frequency, this paper presents a power exponential function droop control (PEFDC) strategy considering the SoC. In this control strategy, the ESU is allowed to adjust the output power adaptively according to its own SoC level during discharge and reaches SoC equilibrium. Simulation models are built to compare the PEFDC strategy with conventional droop control (CDC) and power function droop control (PFDC) approaches. The simulation results illustrate the superiority of the proposed control strategy over the other two methods. Finally, the hardware-in-the-loop experiment is conducted to verify the effectiveness of the PEFDC strategy.


Introduction
The development and use of new energy sources have become trends for the present and the future, as the serious pollution caused by traditional energy sources has brought environmental problems to the fore.Distributed-based microgrids are attracting widespread interest, as they can better exploit the advantages of distributed power sources [1,2].Distributed energy storage technology is an effective measure to meet the needs of new energy consumption and enhance the stability of the microgrid.Multiple energy storage units (ESUs) can be connected to the AC microgrid via inverters to expand the capacity of the microgrid [3][4][5].Stabilizing the frequency of the AC microgrid is the basic task of various control strategies, and each ESU can be cut off or put in depending on the scale of the load, which improves the stability of the microgrid [6,7].
Due to the lack of support from the large power grid, the voltage and frequency of an AC microgrid in islanded mode tend to fluctuate and an effective control strategy is required [8][9][10].When faults and load variations occur, distributed ESUs will adjust their outputs to stabilize the system frequency so that the power quality is guaranteed.Most of the distributed ESUs are connected to the microgrid through converters for capacity expansion.The main control strategies for ESUs include PQ control, droop control, constant voltage control, etc. [11][12][13][14].Droop control has many advantages, such as easy implementation and a plug-and-play feature.Therefore, it is widely used in islanded microgrids [15][16][17].The conventional droop control strategy ignores the existence of line impedance, resulting in deviations in power output of the ESUs.In extreme cases, large circulating currents will occur between the ESUs, affecting the stable operation of the system [1,18].In addition, the power is not shared rationally among ESUs because the conventional droop control does not take the SoC of the batteries into account.In Ref. [19], the multiple of frequency and the differential of the DC-side voltage deviation are introduced to the droop strategy.The system inertia is increased and the dynamic performance of the system is improved.However, the stability of the system decreases significantly in the case of large load fluctuations.The multiple operating modes of the converter are integrated into a single functional loop in Ref. [20], which enhances the stability of the system and reduces transient fluctuations during mode switching.In Ref. [21], a power function droop control based on SoC is proposed, but the method is applied to DC microgrids.In Ref. [22], a simple relationship between droop coefficient and charge state is constructed by introducing SoC directly in the droop coefficient.However, the results show that the control effect needs to be improved.Similarly, a charge control strategy for energy storage based on SoC is investigated in Ref. [23].In Ref. [24], a charge-state-based AC microgrid power control strategy is studied, which still provides good control performance when the SoC difference between the ESUs is small.Nonetheless, the shortcoming is that the transmission of the SoC reference value requires communication among the ESUs, which is in conflict with the plug-and-play principle.Moreover, the stability analysis is missed.Another droop control method considering the SoC of the battery is proposed in Ref. [25], where the ESU adjusts its output power depending on its SoC.However, the method is only available for ESUs with an SoC of more than 30%.The method introduces only one exponential factor, the size of which directly affects the amplitude of frequency fluctuation, the control performance is deteriorated due to the lack of constraints from other factors.
On the basis of research efforts reported in studies, this paper proposes a power exponential function droop control (PEFDC) strategy for multiple-ESU applications in an islanded microgrid based on the SoC of each ESU.The proposed method constructs an improved power exponential function between the SoC of ESU and the droop coefficient.Compared with other methods, this strategy not only stabilizes the AC bus voltage and the frequency when the load changes, but also enables each ESU to adjust its power output dynamically depending on its own SoC.The ESU with higher SoC produces more power and the one with lower SoC produces less power.Finally, reasonable sharing of power and SoC balance are achieved.Moreover, this droop method can work at low SoC and works by measuring local information, so communication is not required between ESUs.
The remainder of the paper is organized as follows: Section 2 illustrates the typical topology of the AC microgrid and the basic principle of the droop control.Section 3 introduces the proposed power exponential function droop control.Section 4 analyzes the stability by the small signal method.Section 5 gives the simulation and hardware-in-theloop experiment results, several droop methods are compared and the design procedures of the factors are also discussed.Section 6 gives the conclusion.

System Structure
A typical topology of an AC microgrid is shown in Figure 1, where n multiple energy ESUs are connected to the AC bus via bidirectional DC/ac converters.The ac grid consists of 'x' DC sources and 'y' ac sources connected through DC/AC converters and AC/AC converters respectively, while AC loads are connected to the AC bus directly.The microgrid has two operation modes, i.e., grid-connected mode and islanded mode, which can be switched via a static transfer switch (STS).In grid-connected mode, the frequency and voltage of the AC microgrid are supported by the large grid and the power is shared among the micro-sources by means of power tracking [26].Each micro-source operating in islanded mode supports the load independently and stabilizes the voltage and frequency of the microgrid to ensure the power quality when load variation occurs.
In islanded operation mode of the microgrid, the control strategy enables coordinated control of the charging and/or discharging of each ESU, which not only helps to stabilize the voltage and frequency, thus improving the power quality, but also allows each ESU to adjust its output active power in real time according to its own SoC.ESUs with high SoC deliver more power, while those with low SoC deliver less power.Eventually, power balance is achieved, namely, power is shared among the ESUs in a reasonable way, thus making the most of each ESU.In islanded operation mode of the microgrid, the control strategy enables coordinated control of the charging and/or discharging of each ESU, which not only helps to stabilize the voltage and frequency, thus improving the power quality, but also allows each ESU to adjust its output active power in real time according to its own SoC.ESUs with high SoC deliver more power, while those with low SoC deliver less power.Eventually, power balance is achieved, namely, power is shared among the ESUs in a reasonable way, thus making the most of each ESU.

Droop Control of AC microgrids
In AC microgrids, the use of simple conventional droop control in a distributed energy storage system allows the power to be shared proportionally according to the droop coefficients without communication [27].The droop characteristic curve of an AC microgrid is shown in Figure 2, where the frequency and voltage droop coefficients are kp and kq, respectively.In steady state, the AC microgrid operates at Point A, where the voltage and frequency are Vac,1 and fac,1, respectively.When the load changes, the ESU can adjust the output power dynamically according to the droop characteristics.For example, when the load increases, the output power of the converter will increase, and this leads to decreases in frequency and voltage.The operating state will move from Point A to Point B. Therefore, the active and reactive power outputs of the converter can be adjusted by varying the values of the droop factors on the basis of frequency and voltage stability.
The equations of conventional droop control can be obtained as where  and  are the actual frequency and voltage of the microgrid, respectively,

Droop Control of AC Microgrids
In AC microgrids, the use of simple conventional droop control in a distributed energy storage system allows the power to be shared proportionally according to the droop coefficients without communication [27].The droop characteristic curve of an AC microgrid is shown in Figure 2, where the frequency and voltage droop coefficients are k p and k q , respectively.In steady state, the AC microgrid operates at Point A, where the voltage and frequency are V ac,1 and f ac,1 , respectively.When the load changes, the ESU can adjust the output power dynamically according to the droop characteristics.For example, when the load increases, the output power of the converter will increase, and this leads to decreases in frequency and voltage.The operating state will move from Point A to Point B. Therefore, the active and reactive power outputs of the converter can be adjusted by varying the values of the droop factors on the basis of frequency and voltage stability.In islanded operation mode of the microgrid, the control strategy enables coordinated control of the charging and/or discharging of each ESU, which not only helps to stabilize the voltage and frequency, thus improving the power quality, but also allows each ESU to adjust its output active power in real time according to its own SoC.ESUs with high SoC deliver more power, while those with low SoC deliver less power.Eventually, power balance is achieved, namely, power is shared among the ESUs in a reasonable way, thus making the most of each ESU.

Droop Control of AC microgrids
In AC microgrids, the use of simple conventional droop control in a distributed energy storage system allows the power to be shared proportionally according to the droop coefficients without communication [27].The droop characteristic curve of an AC microgrid is shown in Figure 2, where the frequency and voltage droop coefficients are kp and kq, respectively.In steady state, the AC microgrid operates at Point A, where the voltage and frequency are Vac,1 and fac,1, respectively.When the load changes, the ESU can adjust the output power dynamically according to the droop characteristics.For example, when the load increases, the output power of the converter will increase, and this leads to decreases in frequency and voltage.The operating state will move from Point A to Point B. Therefore, the active and reactive power outputs of the converter can be adjusted by varying the values of the droop factors on the basis of frequency and voltage stability.
The equations of conventional droop control can be obtained as where  and  are the actual frequency and voltage of the microgrid, respectively,  and  are the active and reactive power outputs of the converter after low-pass The equations of conventional droop control can be obtained as where f ac and V ac are the actual frequency and voltage of the microgrid, respectively, p acl and q acl are the active and reactive power outputs of the converter after low-pass filtering, respectively, f * ac and V * ac are the reference values of frequency and voltage, respectively, P * ac and Q * ac are the nominal powers of the converter, respectively, while k p and k q are the droop control coefficients, respectively.

Improved Droop Control Strategy
To solve the problem of unreasonable active power sharing in conventional droop control, this paper proposes an improved droop control strategy, which takes the state of charge of the ESU into account, and no communication is required among multiple ESUs to achieve reasonable sharing of the active power.The proposed control strategy constructs a power exponential function of the droop factor with the SoC of the ESU.The improved droop expression is where m, n and p are the correction factors of the droop coefficient, and SoC i is the state of charge of the ith ESU.From Equation (3), the active power output of two paralleled ESUs during discharge is related to the SoC, and the ratio of the active power deviations is expressed as As can be seen from Equation ( 4), the output power of the ESU is positively related to the battery's SoC.The factors m, n and p mainly affect the rate of equilibrium and cause the system frequency fluctuation.The effects of these three factors are analyzed in detail below.

Mathematical Model of the System
A mathematical model is developed to analyze the effect of the factors m, n and p on system frequency and load power sharing performance.The SoC estimation of lithium-ion batteries usually uses the coulomb counting method [28]; the calculation formula is where SoC i0 is the initial charge state of the lithium battery.The output power of the ESU is where U bati and i bati are the output voltage and current of the lithium battery, respectively, C e is the capacity of the battery.So, the SoCs of the two batteries can be expressed as Since the frequency is a global variable, it is clear from Equation (3) that the power output of the converter only depends on the droop coefficient, so Equations ( 7) and ( 8) can be further deduced as where P load is the active power of the load.Mathematical simulations were carried out for Equations ( 9) and (10); the results are shown in Figure 3.As can be seen from the figure, the SoC of the two ESUs approach together gradually.
where  is the active power of the load.Mathematical simulations were carried out for Equations ( 9) and ( 10); the results are shown in Figure 3.As can be seen from the figure, the SoC of the two ESUs approach together gradually.

Proposed Converter Schemes
The detailed structure of the battery-dominated hybrid microgrid under islanded operation is shown in Figure 4. Energy storage batteries connect to the AC microgrid via DC/AC converters.LC filter is used at the output of each converter to improve the power quality.The load is connected directly to the AC bus.If the phase and amplitude of the voltage are known, the reference value of the voltage can be obtained by means of voltage composition:

Proposed Converter Schemes
The detailed structure of the battery-dominated hybrid microgrid under islanded operation is shown in Figure 4. Energy storage batteries connect to the AC microgrid via DC/AC converters.LC filter is used at the output of each converter to improve the power quality.The load is connected directly to the AC bus.
where  is the active power of the load.Mathematical simulations were carried out for Equations ( 9) and (10); the results are shown in Figure 3.As can be seen from the figure, the SoC of the two ESUs approach together gradually.

Proposed Converter Schemes
The detailed structure of the battery-dominated hybrid microgrid under islanded operation is shown in Figure 4. Energy storage batteries connect to the AC microgrid via DC/AC converters.LC filter is used at the output of each converter to improve the power quality.The load is connected directly to the AC bus.If the phase and amplitude of the voltage are known, the reference value of the voltage can be obtained by means of voltage composition: If the phase and amplitude of the voltage are known, the reference value of the voltage can be obtained by means of voltage composition: Figure 5 shows the process of reference voltage generation, and the key is the phase angle and amplitude of the voltage, which can be obtained by Equations ( 2) and (3).They are known as frequency droop and voltage droop.
Electronics 2021, 10, x FOR PEER REVIEW 6 of 16 Figure 5 shows the process of reference voltage generation, and the key is the phase angle and amplitude of the voltage, which can be obtained by Equations ( 2) and (3).They are known as frequency droop and voltage droop.After obtaining the reference voltage, the d-axis and q-axis components under the stationary coordinate can be obtained by Park transformation.The d-axis voltage is coupled with the q-axis voltage.In order to facilitate the design of the controller and to achieve independent control, a front decoupling should be performed before the current controller.It should be noted that a large virtual impedance is added before obtaining the voltage reference values for the d and q axes to share the reactive power.The approximate decoupling of the converter output power is achieved by remodeling the equivalent output impedance without changing the control parameters.Figure 6 shows the structure of the voltage and current double closed-loop control, where the d-axis voltage is compared with the d-axis reference voltage and the error is sent to a PI controller to obtain the d-axis reference current.Then, the reference current is sent through the current inner-loop controller and the output is added to the d-axis voltage to obtain a reference control voltage, which is used as a reference signal for PWM.

Small Signal Stability Analysis
From Equation (3), the improved droop expressions for the two paralleled ESUs are obtained as The small signal perturbations of  1 and  2 are After obtaining the reference voltage, the d-axis and q-axis components under the stationary coordinate can be obtained by Park transformation.The d-axis voltage is coupled with the q-axis voltage.In order to facilitate the design of the controller and to achieve independent control, a front decoupling should be performed before the current controller.It should be noted that a large virtual impedance is added before obtaining the voltage reference values for the d and q axes to share the reactive power.The approximate decoupling of the converter output power is achieved by remodeling the equivalent output impedance without changing the control parameters.Figure 6 shows the structure of the voltage and current double closed-loop control, where the d-axis voltage is compared with the d-axis reference voltage and the error is sent to a PI controller to obtain the daxis reference current.Then, the reference current is sent through the current inner-loop controller and the output is added to the d-axis voltage to obtain a reference control voltage, which is used as a reference signal for PWM.
Electronics 2021, 10, x FOR PEER REVIEW 6 of 16 Figure 5 shows the process of reference voltage generation, and the key is the phase angle and amplitude of the voltage, which can be obtained by Equations ( 2) and (3).They are known as frequency droop and voltage droop.After obtaining the reference voltage, the d-axis and q-axis components under the stationary coordinate can be obtained by Park transformation.The d-axis voltage is coupled with the q-axis voltage.In order to facilitate the design of the controller and to achieve independent control, a front decoupling should be performed before the current controller.It should be noted that a large virtual impedance is added before obtaining the voltage reference values for the d and q axes to share the reactive power.The approximate decoupling of the converter output power is achieved by remodeling the equivalent output impedance without changing the control parameters.Figure 6 shows the structure of the voltage and current double closed-loop control, where the d-axis voltage is compared with the d-axis reference voltage and the error is sent to a PI controller to obtain the d-axis reference current.Then, the reference current is sent through the current inner-loop controller and the output is added to the d-axis voltage to obtain a reference control voltage, which is used as a reference signal for PWM.

Small Signal Stability Analysis
From Equation (3), the improved droop expressions for the two paralleled ESUs are obtained as The small signal perturbations of  and  are

Stability Analysis 4.1. Small Signal Stability Analysis
From Equation (3), the improved droop expressions for the two paralleled ESUs are obtained as The small signal perturbations of f ac1 and f ac2 are where the capital letters represent steady state quantities and the symbols with 'ˆ' mean small signal perturbation quantities.Small signal perturbation is applied to Equations ( 7) and ( 8), the results in frequency domain are expressed as The power from the ESU needs to be filtered by a low-pass filter: pacli = H l paci (18) where H l is the transfer function of a second-order low-pass filter: where ξ is the damping coefficient and ω 0 is the cut-off frequency.After eliminating pac1 and pac2 by substituting Equation ( 18) into Equations ( 16) and ( 17), we obtain After substituting ( 20) into ( 14) and ( 21) into ( 15), ŜoC 1 and ŜoC 2 are eliminated: (23) To simplify Equations ( 22) and ( 23), set g im1 and g im2 as Thus, Equations ( 22) and ( 23) can be expressed as In a steady microgrid, the power is balanced and the frequency is the same everywhere: Apply small signal perturbation to ( 28) and (29), then solve simultaneously with ( 22) and ( 23), the characteristic equation of the system can be obtained as where The stability of the system with PEFDC can be analyzed through the characteristic equation.ESU 1 and ESU 2 have the same rated voltage and the same capacity as 500 V and 20 Ah, respectively.The initial droop factor k p is 0.000001.The cut-off frequency and damping factor of the low-pass filter are 157 rad/s and 0.707, respectively.Figure 7a shows the change of eigenvalues when m = 3, n = 5, p = 5, SoC 1 = 0.9 and SoC 2 changes from 0.8 to 0.1 with an interval of 0.1.As SoC decreases, the main eigenvalue moves to the left away from the imaginary axis on the real axis and all eigenvalues are on the left half-plane, so the system is stable.The relationship between SoC 1 and SoC 2 is equivalent and will not be analyzed further.Figure 7b shows the change of eigenvalues when m = 3, n = 5, SoC 1 = 0.9, SoC 2 = 0.8 and p varies from 1 to 200 with an interval of 1.It can be observed that the eigenvalues are always on the left half-plane as p increases.When p is 163, the system starts to appear with a pair of conjugate complex eigenvalues (−156.8± 7.715j) at which the system will start to oscillate, so it is not appropriate to take a value of p greater than 163.As m = 3, p = 5,SoC 1 = 0.9 and SoC 2 SoC 2 = 0.8, the change of eigenvalues when n varies from 30 to 1 with an interval of 1 is analyzed, which is shown in Figure 7c.As n decreases, the dominant eigenvalues move away from the imaginary axis, so the rapidity of the system is enhanced.Because all the eigenvalues are on the real axis on the left half-plane, the system is smooth and free of overshoot.In Figure 7d, as n = 5, p = 5,SoC 1 = 0.9 and SoC 2 = 0.8, the system is stable as all eigenvalues are on the negative half-axis of the real axis as the factor m increases from 1 to 30. where The stability of the system with PEFDC can be analyzed through the characteristic equation.ESU 1 and ESU 2 have the same rated voltage and the same capacity as 500 V and 20 Ah, respectively.The initial droop factor  is 0.000001.The cut-off frequency and damping factor of the low-pass filter are 157 rad/s and 0.707, respectively.Figure 7a shows the change of eigenvalues when m = 3, n = 5, p = 5,  = 0.9 and  changes from 0.8 to 0.1 with an interval of 0.1.As SoC decreases, the main eigenvalue moves to the left away from the imaginary axis on the real axis and all eigenvalues are on the left halfplane, so the system is stable.The relationship between  and  is equivalent and will not be analyzed further.Figure 7b shows the change of eigenvalues when m = 3, n = 5,  = 0.9,  = 0.8 and p varies from 1 to 200 with an interval of 1.It can be observed that the eigenvalues are always on the left half-plane as p increases.When p is 163, the system starts to appear with a pair of conjugate complex eigenvalues (−156.8± 7.715j) at which the system will start to oscillate, so it is not appropriate to take a value of p greater than 163.As m = 3, p = 5,  = 0.9 and  = 0.8, the change of eigenvalues when n varies from 30 to 1 with an interval of 1 is analyzed, which is shown in Figure 7c.
As n decreases, the dominant eigenvalues move away from the imaginary axis, so the rapidity of the system is enhanced.Because all the eigenvalues are on the real axis on the left half-plane, the system is smooth and free of overshoot.In Figure 7d, as n = 5, p = 5,  = 0.9 and  = 0.8, the system is stable as all eigenvalues are on the negative half-axis of the real axis as the factor m increases from 1 to 30.

System Frequency Deviation
In this paper, the sharing of power on the distributed energy storage side is analyzed and the conventional virtual impedance droop control method is used to solve the reactive power sharing problem.In islanded operation, the AC microgrid supported by the ESUs will incur frequency deviations, as can be seen from Equation (3), which is rewritten as As can be seen from Equation (31), the frequency deviation of the AC microgrid in

System Frequency Deviation
In this paper, the sharing of power on the distributed energy storage side is analyzed and the conventional virtual impedance droop control method is used to solve the reactive power sharing problem.In islanded operation, the AC microgrid supported by the ESUs will incur frequency deviations, as can be seen from Equation (3), which is rewritten as As can be seen from Equation (31), the frequency deviation of the AC microgrid in islanded operation is related to the initial droop coefficient and the factors m, n and p.Therefore, in order to ensure that the system frequency deviation is within a reasonable range, an appropriate value for the initial droop coefficient and the factors m, n and p need to be specified.It can be seen from Equation (3) that the droop coefficient of the converter satisfies the following relation: In order to ensure stable operation of the system within a reasonable frequency range, we need to ensure that the product of the maximum value of the droop coefficient and the difference between the output power and the rated power is less than the maximum allowed deviation of the AC microgrid system frequency, i.e., satisfying Equation (33).

Simulation Parameters
To verify the effectiveness of the proposed improved droop control strategy, a parallel model of two energy storage batteries was built in the MATLAB/Simulink environment to simulate the islanded operation of a multi-storage AC microgrid.The converter carrier frequency is 6000 Hz, the output frequency of the converters is 50 Hz and the simulation step of the system is 30 µs.Other simulation parameters are shown in Table 1.0.9

Sharing of Load Power
As shown in Table 1, all the parameters in the system are the same, the difference is the state of charge of the energy storage battery, and if required, take n = 5, p = 7 and m = 3.Three methods will be compared: the conventional droop control (CDC), the power function droop control (PFDC) and the power exponential function droop control (PEFDC).Figure 8a shows the comparison of the active power output of the ESU (ESU) where the blue, green and red lines represent the active power output of the ESUs with the application of CDC, PFDC and PEFDC, respectively.As can be seen from the figure, the CDC strategy allows the active power to be shared evenly, but cannot adjust the output active power according to the actual remaining power of the ESU.The green lines (PFDC) and red lines (PEFDC) show that ESUs can output active power according to their SoC, and the red lines represent the better control performance.Figure 8b shows the SoC difference between the two ESUs.The SoC difference between ESUs cannot be equalized and even increases when applying CDC.The PFDC can reduce the difference of SoC, but the SoC equalization process is slow.Finally, the PEFDC strategy achieves fast SoC equalization and the ∆SoC approaches 0 when t = 2000 s.
red lines (PEFDC) show that ESUs can output active power according to their SoC, and the red lines represent the better control performance.Figure 8b shows the SoC difference between the two ESUs.The SoC difference between ESUs cannot be equalized and even increases when applying CDC.The PFDC can reduce the difference of SoC, but the SoC equalization process is slow.Finally, the PEFDC strategy achieves fast SoC equalization and the  approaches 0 when t=2000s.

Comparisons
Table 2 compares the three droop control methods applied to the AC microgrid.It can be seen that the proposed power exponential droop control method has superiority over the other methods.The control method can be used in AC microgrids with multiple energy storages.

Role and Design of Factors m, n and p
From Section 4.1, it is known that the system is stable over the range of  ∈ [1,30],  ∈ [1,30] and  ∈ [1,163] during the variation of SOC.To further narrow the range of

Comparisons
Table 2 compares the three droop control methods applied to the AC microgrid.It can be seen that the proposed power exponential droop control method has superiority over the other methods.The control method can be used in AC microgrids with multiple energy storages.From Section 4.1, it is known that the system is stable over the range of n ∈ [1, 30], m ∈ [1, 30] and p ∈ [1, 163] during the variation of SoC.To further narrow the range of values taken, the following should be done.Given by (3), a droop coefficient function Z(SoC) is defined as The effect of the parameters m, n, p and SoC on Z(SoC) is shown in Figure 9. Figure 9a is a surface plot of the droop coefficient for m = 3, n = 5, SoC ∈ [0.1, 0.95] and p ∈ [0, 10].It can be seen from the plot that when SoC varies and p is in the range of [0,9], the variation of Z(SoC) is smoother, which is beneficial to the design of the parameters, so p can be taken in the range of [0,9].Figure 9b is a surface plot of Z(SoC) for m = 3, p = 5, SoC ∈ [0.1, 0.95] and n ∈ [0, 10].As can be seen from the graph, the surface plot is relatively smooth, so n can be taken in the range of [1,10].Similarly, Figure 9c is a surface plot of Z(SoC) for n = 5, p = 5, SoC ∈ [0.1, 0.95] and m ∈ [0, 10]; the surface is smooth, as m varies and is suitable for the design of the parameters.
Figure 10 shows the variation of the ESU output power and the fluctuation of system frequency when the factor n is taken to different values, while m = 3 and p = 5.It can be seen from Figure 10a,b that the factor n has little influence on the power output from the ESUs, but mainly affects the frequency of the system.The larger n is, the smaller the frequency fluctuation of the system.Therefore, the system frequency can be stabilized by adjusting the value of n.
[0,10].It can be seen from the plot that when SOC varies and p is in the range of [0,9], the variation of Z(SOC) is smoother, which is beneficial to the design of the parameters, so p can be taken in the range of [0,9].Figure 9b is a surface plot of Z(SOC) for m = 3, p = 5,  ∈ [0.1,0.95] and  ∈ [0,10].As can be seen from the graph, the surface plot is relatively smooth, so n can be taken in the range of [1,10].Similarly, Figure 9c is a surface plot of Z(SOC) for n = 5, p = 5,  ∈ [0.1,0.95] and  ∈ [0,10]; the surface is smooth, as m varies and is suitable for the design of the parameters.Figure 10 shows the variation of the ESU output power and the fluctuation of system frequency when the factor n is taken to different values, while m = 3 and p = 5.It can be seen from Figure 10a,b that the factor n has little influence on the power output from the ESUs, but mainly affects the frequency of the system.The larger n is, the smaller the frequency fluctuation of the system.Therefore, the system frequency can be stabilized by adjusting the value of n.  Figure 11 illustrates the output power of ESU and the system frequency fluctuations when the factor p is taken differently, while m = 3 and n = 5.As can be seen from Figure 11a,b, an increase in p mainly affects the output power of the ESUs.The larger p is, the more significant the distinction between the output powers of the ESUs.Therefore, the system frequency can be stabilized and good control performance is achieved by selecting a reasonable combination of factors n and p.The ESU with high SoC outputs a large power while the other outputs a small power, and they converge gradually.Figure 11 illustrates the output power of ESU and the system frequency fluctuations when the factor p is taken differently, while m = 3 and n = 5.As can be seen from Figure 11a,b, an increase in p mainly affects the output power of the ESUs.The larger p is, the more significant the distinction between the output powers of the ESUs.Therefore, the system frequency can be stabilized and good control performance is achieved by selecting a reasonable combination of factors n and p.The ESU with high SoC outputs a large power while the other outputs a small power, and they converge gradually.
Figure 11 illustrates the output power of ESU and the system frequency fluctuations when the factor p is taken differently, while m = 3 and n = 5.As can be seen from Figure 11a,b, an increase in p mainly affects the output power of the ESUs.The larger p is, the more significant the distinction between the output powers of the ESUs.Therefore, the system frequency can be stabilized and good control performance is achieved by selecting a reasonable combination of factors n and p.The ESU with high SoC outputs a large power while the other outputs a small power, and they converge gradually.

Hardware-in-the-Loop Experiment
To verify the validity of the proposed control strategy, the hardware-in-the-loop experiment is conducted to verify the validity of the algorithm.The experimental environment is shown in the Figure 12.Two main aspects are considered: (a) the correctness of the strategy; (b) the effectiveness of the strategy in the event of sudden changes in load.The experimental parameters are shown in Table 3.

Hardware-in-the-Loop Experiment
To verify the validity of the proposed control strategy, the hardware-in-the-loop experiment is conducted to verify the validity of the algorithm.The experimental environment is shown in the Figure 12.Two main aspects are considered: (a) the correctness of the strategy; (b) the effectiveness of the strategy in the event of sudden changes in load.The experimental parameters are shown in Table 3.At 1600 s, the active load is increased by 5 kW and the active powers from the ESUs are shown in Figure 13.It can be seen that the ESU with a higher charge state outputs more power, while the other outputs less power, and they tend to be the same as the SoC approaches.When the load increases suddenly, the strategy still works.At 1600 s, the active load is increased by 5 kW and the active powers from the ESUs are shown in Figure 13.It can be seen that the ESU with a higher charge state outputs more power, while the other outputs less power, and they tend to be the same as the SoC approaches.When the load increases suddenly, the strategy still works.At 1600 s, the active load is increased by 5 kW and the active powers from the ESUs are shown in Figure 13.It can be seen that the ESU with a higher charge state outputs more power, while the other outputs less power, and they tend to be the same as the SoC approaches.When the load increases suddenly, the strategy still works.Figure 14 shows the waveforms of the bus voltage before and after the increase in load.It is obvious that the voltages are stable.Figure 15 shows the current variation of the converter, from which it can be seen that the converter currents increase rapidly when the load increases and reach stability quickly.Figure 14 shows the waveforms of the bus voltage before and after the increase in load.It is obvious that the voltages are stable.Figure 15 shows the current variation of the converter, from which it can be seen that the converter currents increase rapidly when the load increases and reach stability quickly.

Conclusion
In this paper, an improved droop control strategy of an AC microgrid with multienergy storage is proposed, and a power exponential function between the droop coefficient and SoC is constructed.Then, the small signal perturbation method is used to analyze the stability of the system.In addition, conventional droop and power function droop

Conclusion
In this paper, an improved droop control strategy of an AC microgrid with multienergy storage is proposed, and a power exponential function between the droop coefficient and SoC is constructed.Then, the small signal perturbation method is used to analyze the stability of the system.In addition, conventional droop and power function droop

Conclusions
In this paper, an improved droop control strategy of an AC microgrid with multienergy storage is proposed, and a power exponential function between the droop coefficient and SoC is constructed.Then, the small signal perturbation method is used to analyze the stability of the system.In addition, conventional droop and power function droop methods are also built for comparison.The role of the factors m, n and p and their effects on system stability are also analyzed.Finally, the simulation and HIL experiment are conducted.The results lead to the following conclusions.
(1) In an AC microgrid operating in islanded mode, the proposed control strategy allows the ESU to adjust its output adaptively based on its own SoC.The ESU with high SoC outputs a high power, and the one with low SoC outputs a low power.They output the same power at the same SoC.
(2) Small-signal analysis is used to analyze the stability of the system and to determine the ranges of factors m, n and p.The system frequency can be stabilized and the SoCs can be equalized quickly by adjusting the factors.
(3) The proposed control strategy does not require communication between ESUs.Compared to CDC and PFDC, the PEFDC has a faster SoC equalization rate and remains effective when the ESUs have small SoCs.
(4) The AC microgrid model with multiple EMUs is tested through the hardwarein-the-loop experiment, and the results verify the effectiveness of the proposed power exponential function droop control strategy.

Figure 1 .
Figure 1.A typical topology of an AC microgrid.

Figure 1 .
Figure 1.A typical topology of an AC microgrid.

Figure 1 .
Figure 1.A typical topology of an AC microgrid.

Figure 2 .
Figure 2. Droop characteristic curves of Q-V and P-f.

Figure 4 .
Figure 4. Structure of the battery-dominated hybrid microgrid under islanded operation.

Figure 4 .
Figure 4. Structure of the battery-dominated hybrid microgrid under islanded operation.

Figure 4 .
Figure 4. Structure of the battery-dominated hybrid microgrid under islanded operation.

Figure 5 .
Figure 5. Schematic diagram of the reference voltage generation.

Figure 6 .
Figure 6.Structure of the voltage and current double closed-loop control.

Figure 5 .
Figure 5. Schematic diagram of the reference voltage generation.

Figure 5 .
Figure 5. Schematic diagram of the reference voltage generation.

Figure 6 .
Figure 6.Structure of the voltage and current double closed-loop control.

Figure 6 .
Figure 6.Structure of the voltage and current double closed-loop control.

Figure 7 .
Figure 7. Distribution of eigenvalues: (a) changes with SoC; (b) changes with p; (c) changes with n; (d) changes with m.

Figure 7 .
Figure 7. Distribution of eigenvalues: (a) changes with SoC; (b) changes with p; (c) changes with n; (d) changes with m.

Figure 8 .
Figure 8. Comparisons of the three methods: (a) active power outputs; (b) SoC differences.

Figure 8 .
Figure 8. Comparisons of the three methods: (a) active power outputs; (b) SoC differences.

Figure 10 .
Figure 10.Influence of factors n on active powers and system frequency as m = 3 and p = 5: (a) influence of factor n on power allocation; (b) influence of factor n on frequency fluctuations.

Figure 10 .
Figure 10.Influence of factors n on active powers and system frequency as m = 3 and p = 5: (a) influence of factor n on power allocation; (b) influence of factor n on frequency fluctuations.

Figure 11 .
Figure 11.Influence of factors p on active power and system frequency as m = 3 and n = 5: (a) influence of factor p on power allocation; (b) influence of factor n on frequency fluctuations.

Figure 11 .
Figure 11.Influence of factors p on active power and system frequency as m = 3 and n = 5: (a) influence of factor p on power allocation; (b) influence of factor n on frequency fluctuations.

Figure 14 .Figure 15 .
Figure 14.The voltage variation of the bus.

Figure 15 .
Figure 15.The current variation of the converter.

Table 2 .
Comparisons of droop control methods.

Table 2 .
Comparisons of droop control methods.