Stability Analysis of Multi-Inverters Grid-Connected System Using Improved Current Source Impedance Ratio Criterion.
Recently, Electric vehicles (EVs) are becoming an ecologically attractive alternative to gasoline driven cars. Large-scale integration of the EVs can serve as a distributed and mobile energy source in the electricity market, which construct a new form of energy feedback system , . EVs power station is an intelligent power station with good application prospects. When the energy generated by power systems is insufficient, the station can deliver energy from the local EVs storage to grid , .
Conventional energy feedback system connected to the grid with a single inverter in Point of Common Coupling (PCC), which required storage battery centrally stored and caused a higher grid-connected inverter power capacity. In the EVs power station, multi-points connection to grid is more convenient for the large number of inverters. And the local use of the EVs storage battery can also reduce the power capacity of the inverter.
Usually, the effect of mutual coupling is small in a strong grid even with a large number of inverters at multi-points connection to grid. However, in the EVs power station, the energy feedback system is connected to the grid through medium or low voltage transformer causing a weak grid which would bring stability problem due to the leakage and line inductors. In order to solve this problem, some papers have been published. In  a simplified equivalent method has been proposed based on the small signal model of multi-inverters parallel system. In  equivalent impedance model of the multi-inverters parallel system has been established, and the authors concluded that the resonance peaks of multi-inverters parallel system changed with the number of parallel inverters. The problem of power quality in distributed generation system has been analysed in  caused by series-parallel resonance of a large number of inverters. Furthermore, the literature  pointed that multi-inverters grid-connected system would generate grid voltage harmonic distortion and grid-connected current oscillation in multi-inverters grid-connected distributed generation system.
In this paper, a simplified equivalent model of the multi-inverters grid-connected system is established based on the circuit network theory. The control performance and stability of the system are analysed by improving traditional voltage source impedance criterion to current source impedance ratio criterion. Finally, the simulation results are given to verify the effectiveness of multi-inverter model and the improved current source impedance criterion.
II. MODEL OF MULTI-INVERTERS GRID-CONNECTED SYSTEM
In the EVs power station, the input of the energy feedback system is storage battery and the output is connected to the three-phase grid with output filters. In order to simplify the analysis, the model of multi-inverters grid-connected system can be equivalent as Fig. 1 .
In Fig. 1, LCL filters are used in the system due to their effective switching harmonic attenuation. [Z.sub.11], [Z.sub.21], .... [Z.sub.n1] are inverter side inductor impedances; [Z.sub.12], [Z.sub.22], .... [Z.sub.n2] are grid side inductor impedances; [Z.sub.13], [Z.sub.23], .... [Z.sub.n3] are capacitor impedances; [Z.sub.g] is the grid impedance; [V.sub.g] is the grid voltage, [i.sub.01] ... [i.sub.0n] is the inverter side inductor current. Matrix form of the transfer function between the grid side inductor current [i.sub.11] ... [i.sub.1n] and the output voltage of inverter [v.sub.01] ... [v.sub.0n] can be presented as follows:
[mathematical expression not reproducible]. (1)
G(s) is the coupling function among multi-inverters grid-connected system in (1). Due to the coupling of parameters in multi-inverters system, it is difficult to accurately obtain the expression of G(s). Considering the maturity of production and manufacturing process, assume the output characteristics of inverters are consistent. Therefore all diagonal elements of G(s) are equal [mathematical expression not reproducible]. G(s) can be expressed as a symmetric matrix, and (1) can be rewritten as follows:
[mathematical expression not reproducible]. (2)
According to (2), the relationship between the n-th grid side inductor current [i.sub.1n] and the n-th output voltage of inverter [v.sub.0n] is given in
[G.sub.m] = [G.sub.11][v.sub.0n] + [G.sub.12] ([v.sub.01] + ... + [v.sub.0n-1]). (3)
Considering the output characteristics of multi-inverters are voltage source, the influence of nonlinear modulation can be ignored. The transfer function [G.sub.m] (s) between [i.sub.1n] and [v.sub.0n] can be presented as follows
[G.sub.m](s) = [i.sub.ln]/[v.sub.0n] [approximately equal to] [G.sub.11] + (n - 1) [G.sub.12], (4)
where n [not equal to] 1. In (4), [G.sub.11] reflects the relationship between the output voltage and the output current of single inverter. [G.sub.12] reflects the relationship between the output voltage of other inverters and the output current of the single inverter. According to superposition theorem of the circuit network, [G.sub.11], [G.sub.12] and [G.sub.m] (s) can be obtained respectively.
1. Assuming the output voltage of the inverters is zero, and [mathematical expression not reproducible], Fig. 1 can be equivalent to Fig. 2, and [G.sub.11] can be obtained
[mathematical expression not reproducible], (5)
2. Assuming the output voltage of inverters is zero, Fig. 1 can be equivalent to the Fig. 3, and [G.sub.12] can be obtained
[mathematical expression not reproducible], (6)
3. Based on (5) and (6), [G.sub.m] can be obtained
As shown in (7), the difference between equivalent model of the multiple grid-connected inverters and the single grid-connected inverter is n x [Z.sub.g]. So the multi-inverters can be seen as a single inverter with n time's grid impedance.
III. IMPROVED CURRENT SOURCE IMPEDANCE RATIO CRITERION FOR STABILITY ANALYSIS
Traditional stability analysis usually uses loop gain analysis method based on classical control theory. And closed-loop transfer function of converters is established according to small signal model. When the number of the converters is large, the interaction among individually designed power modules may cause the instability of the whole system. And the orders of the closed-loop transfer function could increase exponentially with the numbers of converters rising, which causes loop gain analysis method become more complex.
As a result, defining impedance specification is an effective way to analyse stability for this system. This method realizes the decoupling of the power supply and the load. Only the input and output impedance characteristics of source and load need to be considered. Traditional voltage source impedance ratio criterion is mainly based on (8). [Z.sub.s] is the impedance of power supply, [Z.sub.i] is the load impedance, [V.sub.s] is the voltage of power supply, [V.sub.i] is the voltage of the load. The impedance ratio [T.sub.m] = [Z.sub.s] / [Z.sub.i] is regarded as loop gain to analyse system stability
[V.sub.i] = [V.sub.s] x [Z.sub.i]/[Z.sub.s] + [Z.sub.i] = [V.sub.s] x 1/1 + [Z.sub.s]/[Z.sub.i]. (8)
However, the traditional impedance ratio criterion is mainly suitable for the applications that characteristics of the power supply and load are voltage source , . The grid-connected system only can be regarded as current source in the view of the grid, and the grid-connected current must maintain the same frequency and phase with the grid to ensure the stability of the system. The traditional impedance criterion is no longer applicable. In this paper an improved current source impedance ratio criterion has been proposed.
Based on the equivalent transformation of circuit networks as shown in Fig. 4, the source is grid-connected inverters and the load is the grid. Therefore, improved current source impedance criterion can be obtained as (9). [I.sub.s] is the equivalent current source. [I.sub.L] is the equivalent load current
[mathematical expression not reproducible], (9)
Compared with (8) and (9), the impedance ratio based on current source mode is [T.sub.m] = [Z.sub.g] / [Z.sub.s]. Considering the disturbance of grid, the output current can be written as follows
[mathematical expression not reproducible]. (10)
Equation (10) is the same with the small signal model of grid-connected inverter. Therefore, the impedance ratio [T.sub.m] = [Z.sub.g] / [Z.sub.s] can be regarded as stability criterion of the grid-connected system.
The stability can be analysed by Nyquist criterion. In order to overcome the conservatism and forbidden region of the improved current source criterion, the following two points are mainly considered.
1. When [absolute value of [T.sub.m]] [less than or equal to] 1, the system can keep stability in full frequency range.
2. When [absolute value of [T.sub.m]] > 1, the stability analysis should combine with the phase margin. Due to the existence of grid impedance and other reactive elements, resonance and harmonic distortion is easy to appear in AC system, which is different from DC system.
IV. STABILITY ANALYSIS OF MULTI-INVERTERS GRID-CONNECTED SYSTEM
A. Impedance Characteristics of LCL Grid-connected System
1) Equivalent Grid Impedance
Due to the nonlinear elements of power loads, it is difficult to obtain accurate model of the grid. Usually it is approximately equivalent to resistance [R.sub.g] and inductance [L.sub.g] in series. The resistance-inductance equivalent model consists with actual measurement of the grid in high frequency range. And the characteristic in low frequency range is mainly decided by transformer leakage reactance and line impedance of middle and low voltage network to PCC. Therefore, the value of grid impedance [Z.sub.g] can be estimated according to the impedance voltage ratio of the transformer. Equivalent model of the grid is given as
[Z.sub.g] = [R.sub.g] + s[L.sub.g]. (11)
2) Output Impedance of LCL Grid-Connected Inverters
In order to eliminate the influence of resonance caused by LCL grid-connected filters, passive damping and active damping are widely used in practice. Active damping with virtual resistance algorithm can effectively suppress resonant peak caused by LCL filters and improve system efficiency , which is used to analyse the output impedance of the grid-connected inverters in this paper.
The dual-loop control diagram of grid-connected inverter with virtual impedance is shown in Fig. 5.
The outer loop is grid-connected current control and the inner loop is capacitor current control. And the open-loop transfer function of the system is given as
[G.sub.open] = [K.sub.pwm] ([K.sub.p]s + [K.sub.i])/[s.sup.2]([s.sup.2][L.sub.1][L.sub.2]C + s[K.sub.c][L.sub.2]C + [L.sub.1] + [L.sub.2]), (12)
where [L.sub.1] is the inverter side inductor, [L.sub.2] is the grid side inductor, C is the parallel capacitor, [K.sub.pwm] is the system gain, [K.sub.p] and [K.sub.i] are parameters of PI controller [G.sub.0](s). [K.sub.c] is the virtual resistance feedback gain. [i.sub.1] is the grid-connected current, [i.sub.r] is given current of grid-connected current, [i.sub.c] is capacitor current, [v.sub.c] is capacitor voltage.
According to Fig. 5, (13) can be obtained as follows:
[mathematical expression not reproducible]. (13)
According to (12) and (13), the output impedance [Z.sub.o] can be obtained
[mathematical expression not reproducible]. (14)
When internal resistance of capacitor and inductor are ignored, (14) can be equivalent to (15)
[mathematical expression not reproducible]. (15)
B. Stability Analysis
According to parameters in Table I, frequency characteristics of [Z.sub.o] and [Z.sub.g] are shown in Fig. 6 in three cases (n = 1, 20, 50).
The performance of the multi-inverters grid-connected system is analysed based on the current source impedance criterion. [Z.sub.o] can be regard as the impedance of power supply in Section III. The impedance ratio [T.sub.m] is given by
[mathematical expression not reproducible]. (16)
The amplitude-frequency characteristics of [T.sub.m] in three cases (n = 1, 20, 50) are shown in Fig. 7. According to impedance characteristics shown in (16), the characteristic in high frequency range is mainly decided by n[L.sub.g] / [L.sub.2], while the coupling point between multi-inverters output impedance and grid impedance often appears in middle and low frequency range. Therefore, stable operating characteristics of the system should consider the phase margin.
It can be seen from Fig. 7 that the stability and harmonic quality meets requirement when n = 1 and [absolute value of [T.sub.m]] is less than 1 in the range of full frequency. When n = 20, impedance crossover frequency is about 800 Hz and the phase margin is about 90[degrees], also the system has good stability. When n = 50, impedance crossover frequency is about 300 Hz, the phase margin is about 80[degrees]. With the increasing of the number of grid-connected inverters, the impedance crossover frequency and phase margin will be decreased regularly. Furthermore, the low-frequency oscillation would be leads to excessive harmonics distortion and resonance.
Nyquist curves in three cases (n = 1, 20, 50) are shown in Fig. 8. The phase margin of multi-inverters grid-connected system is decreased with the increasing of the number of grid-connected inverters. When it is decreased to a certain degree, the system will be unstable due to resonance.
V. SIMULATION RESULTS
Simulation results are present to verify the theoretical analysis. Simulation parameters are shown in Table I.
In order to verify the accuracy of multi-inverters grid-connected equivalent model, two different simulating ways are implemented. And grid-connected current steps at 0.045 s to analyse dynamic characteristics of the model. [I.sub.oeq], [V.sub.seq] are the output current and voltage of the equivalent model proposed in Section II. [I.sub.0], [V.sub.s] are the output current and voltage of multi--paralleled inverters.
When the impedance [Z'.sub.g] of multi-inverters grid-connected system is equal to n time's grid impedance, it can be seen from Fig. 9 that the equivalence keeps correspondence with these two cases.
As shown in Fig. 10, when the grid impedance is not equal to n time's grid impedance, the system does not meet the proposed equivalent relationship and output characteristics are quite different in these two cases.
Also simulation results are present to verify the improved current source impedance ratio criterion. Simulation waveform of a single LCL grid-connected inverter (n = 1) with virtual resistance algorithm is shown in Fig. 11. The power quality of output voltage and grid-connected current is higher.
Simulation waveforms of grid-connected inverters (n = 20) are shown in Fig. 12. The system stability is good and the Total Harmonic Distortion (THD) of grid-connected current is 0.69 %.
Simulation waveforms of grid-connected inverters (n = 50) are shown in Fig. 13. The system is basically stable, but the THD of grid-connected current is 6.31 %. THD exceeds grid connection requirements seriously, especially the low order harmonic current. This situation mainly caused by the lower impedance crossover frequency. It is difficult to suppress by adjusting control bandwidth or algorithm.
It can be seen that simulation results are consistent with the analysis in Section IV. When the number of paralleled-inverters is excessive, the power quality of grid-connected current would be poor. Therefore, the number of the inverters in EVs power station needs to be limited when multi-inverters are connected to the grid at the same time.
In this paper, multi-inverters grid-connected system for EVs Power Station is studied. The equivalent model of the multi-inverters grid-connected system is obtained through the equivalent analysis with circuit network theory. Improved current source impedance is proposed to analysis stability of multi-inverters grid-connected system. And the resonance problem of multi-inverters grid-connected system is also analysed. The stability analysis of multi-inverters grid-connected system can provide theoretical and engineering guidance for construction of EVs power station system.
 B Sun, T. Dragicevic, F. D. Freijedo, J. C. Vasquez, J. M. Guerrero, "A control algorithm for electric vehicle fast charging stations equipped with flywheel energy storage systems", IEEE Trans. Power Electronics, vol. 31, no. 4, pp. 6674-6685, 2016. DOI: 10.1109/TPEL.2015.2500962.
 Z. Xu, W. Su, Z. Hu, Y. Song, H. Zhang, "A hierarchical framework for coordinated charging of plug-in electric vehicles", IEEE Trans. on Smart Grid, vol. 7, no. 1, pp. 428-438, 2016. DOI: 10.1109/TSG.2014.2387436.
 S. Rivera, Bin Wu, "Electric vehicle charging station with an energy storage stage for split-DC bus voltage balancing", IEEE Trans. on Power Electronics, vol. 3, no. 3, pp. 2376-2386, 2017. DOI: 10.1109/TPEL.2016.2568039.
 M. R. Sarker, H. Pandzic, M. A. Ortega-Vazquez, "Optimal operation and services scheduling for an electric vehicle battery swapping station", IEEE Trans. Power Systems, vol. 30, no. 2, pp. 901-910, 2015. DOI: 10.1109/TPWRS.2014.2331560.
 J. H. R. Enslin, P. J. M. Heskes, "Harmonic interaction between a large number of distributed power inverters and the distribution network", IEEE Trans. Power Electronics, vol. 19, no. 6, pp. 1586-1593, 2004. DOI: 10.1109/TPEL.2004.836615.
 Jinwei He, Yun Wei Li, D. Bosnjak, B. Harris, "Investigation and active damping of multiple resonances in a parallel-inverter based microgrid", IEEE Trans. Power Electronics, vol. 20, no. 1, pp. 234-246, 2013. DOI: 10.1109/TPEL.2012.2195032.
 Zhiyong Chen, An Luo, Yandong Chen, Mingshen Li, "Resonance features of multi-paralleled grid-connected inverters and its damping method", in Int. Power Electronics and Application Conf. and Exposition, 2014, pp. 120-125. DOI: 10.1109/PEAC.2014.7037840.
 A. Kotsopoulos, P. J. M. Heskes, M. J. Jansen, "Zero-crossing distortion in grid-connected PV inverters", IEEE Trans. Industrial Electronics, vol. 52, no. 2, pp. 558-565, 2005. DOI: 10.1109/TIE.2005.844254.
 J. L. Agorreta, M. Borrega, J. Lopez, L. Marroyo, "Modeling and control of N-paralleled grid-connected inverters with LCL filter coupled due to grid impedance in PV plants", IEEE Trans. on Power Electronics, vol. 26, no. 3, pp. 770-785, 2011. DOI: 10.1109/TPEL.2010.2095429.
 Byungcho Choi, Dongsoo Kim, Donggyu Lee, Seungwon Choi, Jian Sun, "Analysis of input filter interactionsin switching power converters", IEEE Trans. on Power Electronics, vol. 22, no. 2, pp. 452-460, 2007. DOI: 10.1109/TPEL.2006.889925.
 Xiaogang Feng, Jinjun Liu, Fred Lee, "Impedance specifications for stable DC distributed power systems", IEEE Trans. on Power Electronics, vol. 17, no. 2, pp. 157-162, 2002. DOI: 10.1109/63.988825.
 P. A. Dahono, "A control method to damp oscillation in the input LC filter", in 33rd Annual IEEE Power Electronics Specialists Conf., 2002, pp. 1630-1635. DOI: 10.1109/PSEC.2002.1023044.
Hu Guozhen (1,2), Mao Ling (1), Liu Jun (1), Zhang Lei (1), Fang Zhi Jian (2)
(1) School of Electrical and Electronic Information Engineering, Hubei Polytechnic University, Huangshi, China, 435003
(2) School of Electrical Engineering, Wuhan University, Wuhan, China, 430072
Manuscript received 27 May, 2018; accepted 16 December, 2018.
This research was funded by a grant from China Postdoctoral Science Foundation (168154) and Natural Science Foundation of Hubei Province (2017CFB195, D20164501).
Caption: Fig. 1. Equivalent model of multi-inverters grid-connected system.
Caption: Fig. 2. Equivalent model for calculating [G.sub.11].
Caption: Fig. 3. Equivalent model for calculating [G.sub.12].
Caption: Fig. 4. Equivalent model of improved current source impedance criterion.
Caption: Fig. 5. Dual-loop control diagram of grid-connected inverter.
Caption: Fig. 6. Frequency characteristics of Zo and Zg.
Caption: Fig. 7. Frequency characteristics of impedance ratio [T.sub.m].
Caption: Fig. 8. Nyquist curves of impedance ratio [T.sub.m].
Caption: Fig. 9. Simulation waveforms ([Z'.sub.g] = n[Z.sub.g]).
Caption: Fig. 10. Simulation waveforms ([Z'.sub.g] [not equal to] n[Z.sub.g]).
Caption: Fig. 11. Simulation waveforms of LCL grid-connected system (n = 1).
Caption: Fig. 12. Simulation waveforms of LCL grid-connected system (n = 20): (a) Grid-connected voltage and current waveforms; (b) THD of grid-connected current.
Caption: Fig. 13. Simulation waveforms of LCL grid-connected system (n = 50): a) Grid-connected voltage and current waveforms; b) THD of grid-connected current.
TABLE I. SIMULATION PARAMETERS. Parameters Value DC voltage 390 V [L.sub.1] = [L.sub.2] 0.5 mH C 5 uF [K.sub.p] 9 [K.sub.i] 12500 [L.sub.g] 0.3 mH
|Printer friendly Cite/link Email Feedback|
|Author:||Guozhen, Hu; Ling, Mao; Jun, Liu; Lei, Zhang; Jian, Fang Zhi|
|Publication:||Elektronika ir Elektrotechnika|
|Date:||Apr 1, 2019|
|Previous Article:||Q-Learning Based Failure Detection and Self-Recovery Algorithm for Multi-Robot Domains.|
|Next Article:||A Novel Self-Powered Smart Current Sensor for Power Equipment.|