Active Disturbance Rejection Control Scheme for Reducing Mutual Current and Harmonics in Multi-Parallel Grid-Connected Inverters

With the increasing penetration of renewable energy sources into modern power systems, parallel inverters with LCL filters are commonly employed in the grid interface, giving rise to potential resonance problems. Among the different resonances, interactive resonance is triggered by interaction among inverters when different current references are applied to parallel inverters. It may also feature a mutual current that circulates among inverters instead of flowing into the grid and introduces harmonics or instability in the control system. In this paper, active disturbance rejection control based on a reduced-order extended state observer (RESO) was proposed for parallel inverters. With the proposed scheme, the interaction between inverters is considered as an exogenous disturbance caused by other parallel inverters, estimated by the RESO, and rejected by the controller. In the results, the mutual current and interactive harmonics, calculated via fast Fourier transform, were reduced with the proposed control scheme. Thus, the lower total harmonic distortion of each current was achieved. Additionally, the robust stability and less model-dependent control design are the other additive advantages over the derivative filtered capacitor voltage feedforward-based active damping using PI control. The simulation and real-time experimental results of the conventional and proposed scheme, obtained using the hardware-in-loop, were presented to verify the theoretical analysis under the similar and different current reference cases.


Motivation and Incitement
Interest in distributed power-generation systems based on renewable energy sources (RESs) for sustainable development of the environment is growing. In this context, the pulse-width-modulated grid-connected inverter (GCI) has become the most widespread topology for delivering high-quality power from RESs to the grid. To attenuate the high-frequency switching harmonics generated by pulse-width modulation (PWM), the LCL filter has been increasingly adopted between the inverter and power grid [1].
The inherent resonance in the LCL filter can make the controller very complex and require a proper design with consideration of the grid impedance to avoid instability [2,3]. In renewable energy systems, it is common to connect multiple inverters in parallel with the grid to enhance the total generation capacity, leading to the resonance issue becoming more complex [4]. In such systems, internal, series, and interactive resonances have been observed [5]. Among them, interactive resonance, which is also known as parallel resonance, arises, owing to the mutual interaction between parallel GCIs operating under different current references [5][6][7].

Literature Review
Several studies have been conducted to mitigate the interactive resonance by introducing a virtual harmonic resistance [5], an active damper based on a high-bandwidth power converter [8], or capacitor current-based active damping [9][10][11]. The interactive resonance causes a mutual current-also known as an interactive current-which circulates among inverters instead of flowing into the grid. The relationship between the mutual current and the interactive resonance has been derived, and the adverse effects of the mutual current on the stability of parallel GCIs have been investigated [6,7,12]. However, in these studies, the harmonics introduced by the mutual current were not considered. In [13], the mutual current harmonics were investigated. The harmonics were divided into the low-and high-frequency components. However, the current-reference uncertainty in the RESs was not considered. In previous studies [14][15][16], the instability caused by the mutual current has been investigated under asynchronous PWM carriers or grid-impedance variation. However, the mutual current harmonics may increase the total harmonic distortion (THD) of inverter currents under operation with different references. The general LCL filter design guideline for preventing mutual current instability is provided in [17]. However, the dynamics of currents under different current references were not considered. Thus far, no effective solution for reducing the mutual current and harmonics under different current references in multi-parallel GCIs has been reported.
Active disturbance rejection control (ADRC), which is a disturbance observer-based control method, provides an active and effective way to handle complex/uncertain systems with minimal information about the system [18]. The ADRC estimates all the exogenous and endogenous disturbances, as a value called the lumped disturbance, in the system with the model-independent extended state observer (ESO), which is the main part of ADRC, and compensates them in the control law [19]. The reduction of the resonance harmonics for the LCL filter type GCI was investigated with ADRC in [20][21][22]].

Contribution and Paper Organization
In this paper, a mutual current and resonance harmonics reduction scheme for multi-parallel GCIs with the reduced-order ESO (RESO)-based ADRC, irrespective of the current reference, was proposed. In addition, the common and mutual current expressions were derived with the proposed scheme. The interactive resonance among inverters was treated as an exogenous disturbance, estimated by RESO, and compensated by the proposed scheme. The damping of interactive and common resonances with the proposed method is explained by the frequency response method in the z-domain under different numbers of parallel inverters. The effective resonance damping resulted in greatly reduced mutual current and harmonics in the parallel grid-connected inverters. Furthermore, the stability of the system and disturbance rejection capabilities of the system were analyzed and compared with the conventional PI control using the derivative filter capacitor voltage feedforward-based active damping. To verify the performance of the proposed scheme, several experiments and harmonic analyses were performed for a system composed of two 3kVA inverters in parallel, to consider the worst case, using the real-time hardware-in-loop (HIL) simulation under similar and different reference cases.
The rest of the paper is organized as follows. The mathematical model of the multi-parallel GCI system with LCL filters is derived in Section 2. In Section 3, each element of the overall system with the proposed control scheme is explained, and the design guidelines for the RESO are provided. The common and mutual current expressions in relation to the current reference with the proposed Energies 2019, 12, 4363 3 of 21 method are derived in Section 4. In Section 5, the resonance damping, stability analysis, and disturbance rejection capabilities of the proposed scheme are presented and compared with those of the conventional derivative filtered capacitor voltage feedforward-based active damping. A simulation and real-time experimental verification for validating the analysis performed in the preceding sections under different current-reference conditions are presented in Section 6. Conclusions are drawn in Section 7.

Mathematical Modeling of Multi-Parallel GCIs
The operation of multiple (n) parallel three-phase GCIs with LCL filters is shown in Figure 1. Each inverter is provided with a constant DC voltage V dci (i = 1, . . . , n). Z 1i and Z 2i represent the inverter-and grid-side inductor impedances, respectively, Z 3i represents the filter capacitance impedances, and Z g represents the inductive grid impedance. The inverter-side currents and voltages are represented by i i and v i , respectively. e g represents the ideal and balanced grid voltage, and v pcc represents the voltage at the point of common coupling (PCC). The inverter parameters are assumed to be identical, i.e.,

Mathematical Modeling of Multi-Parallel GCIs
The operation of multiple (n) parallel three-phase GCIs with LCL filters is shown in Figure 1. Each inverter is provided with a constant DC voltage Vdci (i = 1,…, n). Z1i and Z2i represent the inverter-and grid-side inductor impedances, respectively, Z3i represents the filter capacitance impedances, and Zg represents the inductive grid impedance. The inverter-side currents and voltages are represented by i i and v i , respectively. eg represents the ideal and balanced grid voltage, and vpcc represents the voltage at the point of common coupling (PCC). The inverter parameters are assumed to be identical, i.e., Each current i i makes (n − 1) minor contributions to the other j inverters and one major contribution to the grid. The minor contribution is defined as mutual currents i mj (j = 1, . . . , n and j i), which circulate among other inverters. The major contribution is defined as a common current i c. The multi-parallel GCIs is a multi-input-multi-output system [23]. The relationship between the output currents i i and input voltages v i for multi-parallel GCIs can be derived as where Y(s) is the n × n admittance matrix of the system. In Equation (2), the impact of each inverter's own voltage v i (s) on its own current i i (s) is represented by diagonal elements Y ii (s), and the impact of other inverter voltages v j (s) on the current i i (s) is represented by non-diagonal elements Y ij (s). For an n parallel inverter system, an index i corresponds to the inverter under consideration, and an index j corresponds to each of the remaining (n − 1) inverter modules. All the diagonal elements are the same when we assume identical parameters for each inverter. Similarly, all non-diagonal elements are the same. By using the superposition principle, the diagonal and non-diagonal elements can be derived as [6,7,12] Y ii (s) = n−1 n · G plant (s) where G plant (s) and G coupling (s) represent the LCL filter transfer function and grid impedance coupling transfer function, respectively.
Here, ω res and ω res1 are defined as the interactive and common resonance frequencies, respectively.
From Equation (4), it can be easily seen that G plant (s) and G coupling (s) become identical when L g is zero, and accordingly, all the non-diagonal elements Y ij (s) become zero. Consequently, no interaction among inverters is observed [7]. However, when L g is not zero, each inverter interacts with not only a grid but also the other inverters. In this case, the current i i (s) can be expressed as where i ms (s) represents the sum of (n − 1) contributions of current i i (s) as a mutual current. From Equation (2), i i (s) can be expressed as From Equations (3), (6) and (7), the currents i ms (s) and i c (s) are derived as i ms (s) = 1 n G plant (s) where ω g represents the grid frequency. All the parameters of the multi-parallel GCIs are presented in Table 1.

Proposed Mutual Current Reduction Scheme
The overall system with the proposed RESO-based ADRC current control scheme for the multi-parallel three-phase GCIs is shown in Figure 2, where i d i and i q i represent the d and q components of i i , respectively, and v d gi and v q gi represent the d and q components of v pcc , respectively. The current i i , which is influenced by all the other j inverters voltages v j , is chosen as a feedback to the current control.
The v pcc is sensed only for the current and voltage synchronization around ω g with the phase-locked loop (PLL). The current references i d ire f and i q ire f are generated independently by the active power P iref and reactive power Q iref references applied to each inverter, respectively. The six-pulse IGBT gate signals s i1 to s i6 are also generated separately for each inverter by applying space-vector modulation (SVM) to the control signals u abc i provided by each RESO-based ADRC control block after performing inverse dq transformation on control signals u d i and u RESO-based ADRC current control Overall system with the proposed reduced-order extended state observer (RESO)-based active disturbance rejection control (ADRC) current control scheme.

Proposed Mutual Current Reduction Scheme
The overall system with the proposed RESO-based ADRC current control scheme for the multiparallel three-phase GCIs is shown in Figure 2  The detailed diagram with the proposed control scheme from Figure 2 for any inverter i in the z-domain is shown in Figure 3. In this paper, only i d i (z) regulation is explained, for simplicity, and a similar configuration for i q i (z) regulation is adopted. The i d i (z) is compared with i d ire f , followed by proportional control with the gain K p selected according to the required performance of the control. The estimated lumped disturbance z d 2i calculated by the RESO is subtracted to compensate for the exogenous disturbance, i.e., the effect of the sum of the other j inverters voltages. In the result, the control signal u d i (z) is generated after dividing by the turning parameter b, which is equal to 1/L 1 . The z −1 represents an inherent delay caused by computation and the PWM update, and the inverter is simplified as a linear amplifier with gain K PWM . Finally, the v d i (z) is generated. Similarly, v d j (z) is generated by i d j (z) current regulation of any other j inverters. The control plant, which comprises multi-parallel GCIs in the discrete time domain is represented as a block-diagram according to Equation  (9). With the consideration of the zero-order-hold (ZOH) effect in Figure 3, the discrete representations of G d plant (s) and G d coupling (s) are given as [24] Other j# of inverter voltages = The estimated lumped disturbance z 2i d calculated by the RESO is subtracted to compensate for the exogenous disturbance, i.e., the effect of the sum of the other j inverters voltages. In the result, the control signal u i d (z) is generated after dividing by the turning parameter b, which is equal to 1/L1.

Multi-parallel GCIs
The z -1 represents an inherent delay caused by computation and the PWM update, and the inverter is simplified as a linear amplifier with gain KPWM.
current regulation of any other j inverters. The control plant, which comprises multi-parallel GCIs in the discrete time domain is represented as a block-diagram according to Equation (9). With the consideration of the zero-order-hold (ZOH) effect in Figure 3, With the proposed scheme, the interaction influence of other inverters on each current can be treated as an exogenous disturbance, whereas the coupling effect between the dq variables is treated as an endogenous disturbance.
In general, the state-space model of a first-order time-delayed system with an extended state can be represented as Overall system with the proposed RESO-based ADRC current control scheme. Proposed RESO-based first-order ADRC current control scheme for multi-parallel GCIs.
With the proposed scheme, the interaction influence of other inverters on each current can be treated as an exogenous disturbance, whereas the coupling effect between the dq variables is treated as an endogenous disturbance.
In general, the state-space model of a first-order time-delayed system with an extended state can be represented as .
where x 1 (t) and x 2 (t) represent the system states and y(t) is the output of the system. b is the control gain parameter, u c (t − T s ) is the control input with a time delay of sampling time T s , and f (t) represents the lumped disturbance in the system. Generally, the second-order ESO, which is the core part of first-order ADRC, is constructed as where z 1 (t) and z 2 (t) are the estimated values of x 1 (t) and x 2 (t), respectively, and β 1 β 2 is the gain matrix of the ESO. Assuming that x 1 (t) is known through measurement, Equation (13) can be reduced to an estimation of z 2 (t): .
Using Equation (12), we can rewrite Equation (14) as where the observer gain β 2 is selected as an observer bandwidth ω o . After applying the Laplace transform to Equation (15) and separating the two terms u c (s) and y(s), the output-to-estimation transfer function G zy (s) and the control-to-estimation transfer function G zu (s) are derived as G zy (s) = Transforming Equation (16) into the discrete time domain with consideration of the zero-order-hold (ZOH) effect and then replacing y(z) with i d i (z), u c (z) with u d i (z), and z 2 (z) with z d 2i (z) yields the following G zi (z) and G zu (z) in Figure 3:

Modeling of the Proposed Scheme
In this section, the overall dynamics of the system with the proposed control scheme are modeled. From Figure 3, z d 2i (z) and u d i (z) can be expressed as By substituting Equations (17) and (18) into Equation (19), we obtain where G 1 (z) and G 2 (z) are defined as follows: With a well-tuned RESO, the observer output z d 2i (z) closely tracks the lumped disturbances in the multi-parallel GCIs. The proposed scheme can actively compensate for the effect of lumped disturbances in real time. As indicated by Equation (21), the proposed scheme does not rely on detailed knowledge of the system. The expression for v d i (z) from Figure 3 based on Equation (20) are derived as follows: Similarly, the expression for v d j (z) can be derived as From Figure 3, the discrete time-domain representation of i d c (z) can be obtained as By substituting Equations (22) and (23) into Equation (24), i d c (z) can be expressed as By applying the current-separation scheme to Figure 3 for considering i d c (z) only, we obtain By substituting Equation (26) into Equation (25), the relationship between i d c (z) and all the current references can be obtained: Similarly, the discrete time-domain representation of i d ms (z) from Figure 3 can be obtained as By substituting Equations (22) and (23) into Equation (28), i d ms (z) can be expressed as The current i d ms (z) arises only when dissimilar currents exist among inverters; therefore, we obtain n j=1, j i By substituting Equation (30) into Equation (29), the following expression for i d ms (z) is obtained: It can be observed from Equation (31) that when i d ire f = i d jre f , the mutual current dynamics can be ignored.
By applying the dq transformation and discretizing Equation (6), the expression for the current i d i (z) can be obtained: Substituting Equations (27) and (31) into Equation (32) gives the complete expression of i d i (z): The expression in Equation (33) can be simplified by separating the i d i (z)/i d ire f (z) and i d i (z)/ i d jre f (z) characteristics and introducing T(z) and R(z): where G mut (z) and G com (z) are the two closed-loop transfer functions representing the i d ms (z) and i d c (z), respectively, from Equation (34) and are expressed as The performance of the multi-parallel GCIs can be evaluated using Equation (34), which includes i ms (z) and i c (z) dynamics simultaneously.

Resonance Damping and Stability Analysis
The multi-parallel GCIs system is stable if and only if the currents i ms (z) and i c (z) are stable. From Equation (35), the loop-gain expressions T c (z) and T m (z) for the common and mutual current stability, respectively, are obtained as and The resonance damping was analyzed via the frequency-response method using Bode plots with different values of n. The stability analysis was performed on the Bode plot by using the Nyquist stability criterion in the frequency range above 0 dB to find the phase margin (PM) of the system above the −180 • phase line. Additionally, the gain margin (GM) of the system is measured to determine how much increase in the control bandwidth is possible. The parameters of the proposed control scheme are given in Table 2. The parameter b in the proposed scheme plays an important role in achieving a larger observer bandwidth to improve the disturbance estimation performance of the RESO. A higher value of ω o can be obtained when a higher value of b is adopted. Generally, ω o is limited by the sampling frequency f s [20]. The Bode plot of the loop gain of the mutual current with the proposed scheme is compared with that of conventional derivative filtered capacitor voltage-based active damping using PI control [7] in Figure 4. The corresponding loop-gain expressions for the conventional scheme are provided in the Appendix A (Equations (A1) and (A2)). The figure indicates that the interactive resonance in the mutual current is independent of the number of parallel inverters. Hence, the frequency response of the mutual current remains the same under different values of n. A small PM of 4 • is observed with the conventional scheme owing to the ineffective resonance damping. Additionally, the resonance peak at f res amplifies the resonance, which introduces the mutual current in the parallel GCIs. Moreover, a further increase in the bandwidth of the conventional scheme may violate the system stability. However, the proposed scheme effectively damps the interactive resonance and achieves a larger PM of 69 • and a similar crossover frequency as the conventional scheme. With the proposed scheme, no resonance peak is observed above 0 dB at f res , which results in a significantly reduced mutual current in the parallel GCIs.  The frequency response of the loop gain of the common current with the conventional scheme under different values of n is shown in Figure 5. The ω res1 is a variable frequency depending on the number of parallel inverters at the PCC; hence, the frequency-response characteristics of the common current are varied according to n. As indicated by the figure, the conventional scheme cannot achieve enough resonance damping, and its disturbance rejection capabilities are highly affected with an increase in the number of inverters. Additionally, only a small PM in the range of 11°-23° is achieved with the conventional scheme, which limits the bandwidth of the control. In contrast, the proposed scheme effectively damps the common resonance regardless of n and maintains a higher stability of the system, disturbance-rejection capabilities, and the crossover frequency of the control, as shown in Figure 6. Figure 7a,b show the pole-zero maps of the closed-loop transfer functions Gcom_PI(z) and Gmut_PI(z) with the conventional scheme for two parallel GCIs under kp variation. As shown in Figure 7, the conventional scheme introduces a pole near the unit circle in the common and mutual current, which causes resonance in the system and reduces the stability margin of the system. However, in the polezero map of the proposed scheme shown in Figure 8a,b, we observe that the common and interactive resonances are compensated via pole-zero cancelation. Additionally, the common-and mutualcurrent poles are well inside the unit circle. Hence, the entire system is commonly and mutually stable. The frequency response of the loop gain of the common current with the conventional scheme under different values of n is shown in Figure 5. The ω res1 is a variable frequency depending on the number of parallel inverters at the PCC; hence, the frequency-response characteristics of the common current are varied according to n. As indicated by the figure, the conventional scheme cannot achieve enough resonance damping, and its disturbance rejection capabilities are highly affected with an increase in the number of inverters. Additionally, only a small PM in the range of 11 • -23 • is achieved with the conventional scheme, which limits the bandwidth of the control. In contrast, the proposed scheme effectively damps the common resonance regardless of n and maintains a higher stability of the system, disturbance-rejection capabilities, and the crossover frequency of the control, as shown in Figure 6.    Figure 7a,b show the pole-zero maps of the closed-loop transfer functions G com_PI (z) and G mut_PI (z) with the conventional scheme for two parallel GCIs under k p variation. As shown in Figure 7, the conventional scheme introduces a pole near the unit circle in the common and mutual current, which causes resonance in the system and reduces the stability margin of the system. However, in the pole-zero map of the proposed scheme shown in Figure 8a,b, we observe that the common and interactive resonances are compensated via pole-zero cancelation. Additionally, the common-and mutual-current poles are well inside the unit circle. Hence, the entire system is commonly and mutually stable.

Performance Verification
The proposed scheme was verified by simulation and real-time experimental implementation via the HIL technique with a 1 µ s step size for the two-parallel GCIs, where the control algorithm was implemented using real-time digital control circuitry.
The performance of the proposed scheme was investigated under two cases: Case I represents the different current references, i.e., to the q-axis current control of inverters #1 and #2 remain zero in both cases. Figure 9 shows the simulation response of the conventional scheme illustrated in [7]. During Case I, the mutual current is generated in the inverter #2 current, as shown in Figure 9b, owing to the ineffective interactive resonance damping. Consequently, more resonance harmonics are observed in

Performance Verification
The proposed scheme was verified by simulation and real-time experimental implementation via the HIL technique with a 1 µ s step size for the two-parallel GCIs, where the control algorithm was implemented using real-time digital control circuitry.
The performance of the proposed scheme was investigated under two cases: Case I represents the different current references, i.e., to the q-axis current control of inverters #1 and #2 remain zero in both cases. Figure 9 shows the simulation response of the conventional scheme illustrated in [7]. During Case I, the mutual current is generated in the inverter #2 current, as shown in Figure 9b, owing to the ineffective interactive resonance damping. Consequently, more resonance harmonics are observed in

Performance Verification
The proposed scheme was verified by simulation and real-time experimental implementation via the HIL technique with a 1 µs step size for the two-parallel GCIs, where the control algorithm was implemented using real-time digital control circuitry.
The performance of the proposed scheme was investigated under two cases: Case I represents the different current references, i.e., i d 1re f = 5 A and i d 2re f = 0, and Case II represents the same current references, i.e., i d 1re f = 5 A and i d 2re f = 5 A to each inverter. The current references i q 1re f and i q 2re f to the q-axis current control of inverters #1 and #2 remain zero in both cases. Figure 9 shows the simulation response of the conventional scheme illustrated in [7]. During Case I, the mutual current is generated in the inverter #2 current, as shown in Figure 9b, owing to the In Figure 10, the simulation responses of the proposed scheme are shown. The proposed scheme effectively damps the common and interactive resonance in Cases I and II. Consequently, no mutual current is observed in the inverter #2 current of Figure 10b. Additionally, no interactive resonance harmonics are observed in i1 and ic of Figure 10a,c for both Cases I and II. The smooth tracking response of the common current shown in Figure 11c during a step change at 20 ms explains the larger stability margin, as described in Section 5. The real-time experimental results of the conventional scheme for Cases I and II are shown in Figure 11. With the adoption of Case I, the mutual current is activated and circulates in inverter #2, In Figure 10, the simulation responses of the proposed scheme are shown. The proposed scheme effectively damps the common and interactive resonance in Cases I and II. Consequently, no mutual current is observed in the inverter #2 current of Figure 10b. Additionally, no interactive resonance harmonics are observed in i 1 and i c of Figure 10a,c for both Cases I and II. The smooth tracking response of the common current shown in Figure 11c during a step change at 20 ms explains the larger stability margin, as described in Section 5. In Figure 10, the simulation responses of the proposed scheme are shown. The proposed scheme effectively damps the common and interactive resonance in Cases I and II. Consequently, no mutual current is observed in the inverter #2 current of Figure 10b. Additionally, no interactive resonance harmonics are observed in i1 and ic of Figure 10a,c for both Cases I and II. The smooth tracking response of the common current shown in Figure 11c during a step change at 20 ms explains the larger stability margin, as described in Section 5. The real-time experimental results of the conventional scheme for Cases I and II are shown in Figure 11. With the adoption of Case I, the mutual current is activated and circulates in inverter #2,   Figure 12. No mutual current is flowing in i2 when Case I is adopted. In the results, no harmonic oscillation in the steady state is observed in the currents i1, i2 and ic, regardless of the current reference. Additionally, the oscillations in ic during the step change (from Case I to Case II) are removed, which explains the larger stability margin obtained with the proposed scheme, as shown in Figure 12b. Furthermore, the proposed control scheme achieves robust performance with minimum information about the system required without adopting grid voltage feedforward. This theory can also be applied to a system with an uncertain gird impedance or imbalance vpcc.
Case  Figure 11. Experimental results for Cases I and II under the conventional scheme presented in [7]. Figure 11. Experimental results for Cases I and II under the conventional scheme presented in [7].
The real-time experimental results of the conventional scheme for Cases I and II are shown in Figure 11. With the adoption of Case I, the mutual current is activated and circulates in inverter #2, as indicated by i 2 . The mutual current introduces harmonics in i 1 , i 2 and i c during Case I. During Case II, the mutual current is zero, and interactive resonance harmonics do not exist. Additionally, during the step change in inverter #2 from Case I to II, oscillations in the common current are observed, as shown in Figure 11b, which depicts the simulation performance shown in Figure 10.
The real-time experimental performance of the proposed scheme during Cases I and II is shown in Figure 12. No mutual current is flowing in i 2 when Case I is adopted. In the results, no harmonic oscillation in the steady state is observed in the currents i 1 , i 2 and i c , regardless of the current reference. Additionally, the oscillations in i c during the step change (from Case I to Case II) are removed, which explains the larger stability margin obtained with the proposed scheme, as shown in Figure 12b. Furthermore, the proposed control scheme achieves robust performance with minimum information about the system required without adopting grid voltage feedforward. This theory can also be applied to a system with an uncertain gird impedance or imbalance v pcc .  The harmonic analysis using the fast Fourier transform (FFT) for the conventional scheme on i1, i2 and ic is shown in Figure 13a,b for Cases I and II, respectively. The resonance is excited in the common current around fres and fres1, which may result in higher THD of i1 and ic during Case I, as shown in Figure 13a. The THD of i1, i2 and ic is slightly improved during Case II, as shown in Figure  13b.
The harmonic analysis using the FFT for the proposed scheme on currents i1, i2 and ic is shown in Figure 14a,b for each case. The harmonics around fres and fres1 are effectively damped for both Case I and Case II. Additionally, the THD of i1, i2 and ic is greatly improved. During Case I, the THD values for i1 and ic measured at the frequency ωg are calculated as 5.02% and 0.65%, respectively. Similarly, the THD values of i1, i2 and ic during Case II are calculated as 4.84%, 4.84% and 0.20%, respectively. The harmonic analysis using the fast Fourier transform (FFT) for the conventional scheme on i 1 , i 2 and i c is shown in Figure 13a,b for Cases I and II, respectively. The resonance is excited in the common current around f res and f res1 , which may result in higher THD of i 1 and i c during Case I, as shown in Figure 13a. The THD of i 1 , i 2 and i c is slightly improved during Case II, as shown in Figure 13b.
The harmonic analysis using the FFT for the proposed scheme on currents i 1 , i 2 and i c is shown in Figure 14a,b for each case. The harmonics around f res and f res1 are effectively damped for both Case I and Case II. Additionally, the THD of i 1 , i 2 and i c is greatly improved. During Case I, the THD values for i 1 and i c measured at the frequency ω g are calculated as 5.02% and 0.65%, respectively. Similarly, the THD values of i 1 , i 2 and i c during Case II are calculated as 4.84%, 4.84% and 0.20%, respectively.

Conclusions
An RESO-based ADRC control scheme is proposed for damping the interactive and common resonances present in LCL filter-type three-phase multi-parallel GCIs regardless of the current reference. The proposed scheme treats the interaction among inverters as an exogenous disturbance caused by other inverters in parallel, estimates it with the RESO, and compensates disturbances in

Conclusions
An RESO-based ADRC control scheme is proposed for damping the interactive and common resonances present in LCL filter-type three-phase multi-parallel GCIs regardless of the current reference. The proposed scheme treats the interaction among inverters as an exogenous disturbance caused by other inverters in parallel, estimates it with the RESO, and compensates disturbances in

Conclusions
An RESO-based ADRC control scheme is proposed for damping the interactive and common resonances present in LCL filter-type three-phase multi-parallel GCIs regardless of the current reference. The proposed scheme treats the interaction among inverters as an exogenous disturbance caused by other inverters in parallel, estimates it with the RESO, and compensates disturbances in the control law of ADRC. The larger stability margins for the common and mutual current, the preserved disturbance-rejection capabilities under different values of n, and the less model-dependent control design compared to conventional control scheme are the potential benefits of the proposed scheme. The performance of the proposed control scheme was verified and compared with that of the conventional control scheme presented in [7]. A harmonic analysis using the FFT indicates the low THD performance of the proposed method regardless of the current-reference conditions in the system with two parallel GCIs. Thus, the proposed scheme can be extended to a system with an uncertain grid impedance or grid voltage imbalance.