22 Anti Islanding
22 Anti Islanding
22 Anti Islanding
RJ, Brazil
⎡ ⎛ 2 ⎞ ⎛ 2 ⎞⎤
sin θ sin⎜ θ − π ⎟ sin ⎜ θ + π ⎟⎥ ⎡va ⎤
⎡vd ⎤ 2 ⎢ ⎝ 3 ⎠ ⎝ 3 ⎠ ⎢ ⎥
⎢v ⎥ = ⎢ ⎥ ⋅ vb
⎢ ⎥ (1)
⎣ q ⎦ 3 ⎢cos θ cos⎛⎜ θ − 2 π ⎞⎟ cos⎛⎜ θ + 2 π ⎞⎟⎥ ⎢v ⎥
⎢⎣ ⎝ 3 ⎠ ⎝ 3 ⎠⎦ ⎣ c ⎦
⎥
Grid/Load
abc
PLL dq
Fig. 6. Diagram block of the current controller of VSI. Fig. 8. GEVS extensions to constant current-controlled three-phase inverter
systems.
Figure 7 shows the GEVS diagram block and Fig. 8 shows the
modification into inverter control because of GEVS inclusion.
Band-Pass
Filter Gain
Fig. 9. Three-phase diagram of the IBDG connected to the system through an
R+jωL distribution line feeding a RLLLCL parallel local load.
Or in dq coordinates:
BPF
Fig. 7. GEVS diagram block.
⎡iRa ⎤ ⎡va ⎤
⎢i ⎥ = 1 ⋅ ⎢v ⎥ The steady-state operation point x0 can be obtained through a
(15)
⎢ Rb ⎥ R ⎢ b ⎥ steady-state analysis like a load flow.
⎣⎢iRc ⎦⎥ ⎢⎣vc ⎦⎥
L
⎡ a⎤
v ⎡iLa ⎤ The linearized equations are then:
⎢v ⎥ = L ⋅ d ⎢i ⎥
(16)
⎢ b⎥ L
dt ⎢ ⎥
Lb
⎢⎣vc ⎥⎦ ⎢⎣iLc ⎥⎦ Δvd = RΔiNd + LpΔiNd − ω0 LΔiNq − LiNq0 ΔωPLL − E sin(δ 0 )Δδ PLL (27)
Δvq = RΔiNq + LpΔiNq + ω0 LΔiNd + LiNd 0 ΔωPLL − E cos(δ 0 )Δδ PLL (28)
⎡iCa ⎤ ⎡v a ⎤
⎢i ⎥ = C ⋅ d ⎢v ⎥ Δid = ΔiRd + ΔiLd + ΔiCd + ΔiNd (29)
(17)
⎢ Cb ⎥ L
dt ⎢ ⎥
b
Δiq = ΔiRq + ΔiLq + ΔiCq + ΔiNq (30)
⎣⎢iCc ⎦⎥ ⎣⎢vc ⎦⎥
ΔiRd = (1 / RL )Δvd (31)
Or in dq coordinates: ΔiRq = (1 / RL )Δvq (32)
Δvd = LL pΔiLd − ω0 LL ΔiLq − LL iLq 0 ΔωPLL (33)
⎡iRd ⎤ 1 ⎡vd ⎤ Δvq = LL pΔiLq + ω0 LL ΔiLd + LLiLd 0 Δω PLL (34)
⎢i ⎥ = ⋅⎢ ⎥ (18)
⎣ Rq ⎦ RL ⎣vq ⎦ ΔiCd = C L pΔvd − ω0CL Δvq − CL vq 0 ΔωPLL (35)
⎡vd ⎤ ⎡ LL p − ω PLL LL ⎤ ⎡iLd ⎤ ΔiCq = CL pΔvq + ω0CL Δvd + CL vd 0 ΔωPLL (36)
⎢v ⎥ = ⎢ ⋅⎢ ⎥ (19)
⎣ q ⎦ ⎣+ ω PLL LL LL p ⎥⎦ ⎣iLq ⎦ Δud = Ls pΔid (37)
⎡iCd ⎤ ⎡ CL p − ω PLLCL ⎤ ⎡vd ⎤ Δu q = Ls pΔiq (38)
⎢i ⎥ = ⎢ ⋅⎢ ⎥ (20)
⎣ Cq ⎦ ⎣+ ω PLLC CL p ⎥⎦ ⎣vq ⎦ pΔu d = k pi p (Δidref + Δidp − Δid ) + kii (Δidref + Δidp − Δid ) (39)
pΔuq = k pi p (Δiqref − Δiq ) + kii (Δiqref − Δiq ) (40)
The (vd,vq) versus (id, iq) DG terminal characteristics is given
pΔω PLL = k pPLL pΔvq + kiPLL Δvq (41)
by equation (5). Assuming the equivalent distribution system
bus as a perfectly balanced three-phase voltage source, stiff, pΔδ PLL = Δω PLL (42)
with constant angular frequency ω0 and magnitude E, then: Δidp = K v Δvd (43)
0.016
Model Validation 0.014
Δid (p.u.)
confronted against the average model of the inverter so we can 0.01
conclude about your validation.
0.008
0.004
The single DG system shown in Fig. 9 average model were
implemented in SimPowerSystems™/Matlab®. Equation (44), 0.002
which contains both Differential and Algebraic Equations
0
(DAE), was integrated in time domain using Matlab® ode15s, 2 2.02 2.04 2.06 2.08 2.1
Time (s)
a variable order solver based on the numerical differentiation
Fig. 11 Dynamic response of direct axis current id.
formulas (NDF). Because of the model’s algebraic equations,
matrix E is singular and ode15s can solve DAE with singular id Step Response
mass matrix. ref
0.02
Average Model
The parameters used in the single DG model are listed in Table 0.018 Small−Signal Model
1 and the response to a step change of 10% in idref is shown in 0.016
Fig. 10–Fig. 13.
0.014
0.012
Table 1. Parameters used in single DG simulations.
ΔP (p.u.)
0.005
0.8
Δvd (p.u.)
0
0.6
Δω (rad/s)
−0.005 0.4
−0.01 0.2
0
−0.015
−0.2
−0.02
2 2.02 2.04 2.06 2.08 2.1
Time (s) −0.4
Fig. 10. Dynamic response of direct axis voltage vd. 2 2.02 2.04 2.06 2.08 2.1
Time (s)
Fig. 13 Dynamic response of angular frequency ω.
One can assess the stability of the IBDG grid-connected Proc. 2007 IEEE Power Engineering Society General Meeting, pp.1-4,
24-28 Jun 2007.
operation because of the positive feedback control through the
[13] E. Virnik, “Stability Analysis of Positive Descriptor Systems” Linear
eigenanalysis of matrix pencil A,E in equation (44). The Algebra and its Applications, vol.429, no.10, pp.2640-2659, Nov. 2008.
system is asymptotically stable if all its finite generalized
eigenvalues have negative real part [13]. The generalized Acknowledgments
eigenvalues λ of the matrix pencil A, E are the roots of the
characteristic equation: The authors would like to thanks São Paulo State Research
Foundation (FAPESP) and Coordination for the Improvement
det (A − λE) = 0 (47)
of Higher Level Personnel (CAPES) for the support of this
work.
Conclusions
Tiago R. Ricciardi is M.Sc. graduate student with the Department of Electrical
Energy Systems of University of Campinas, Campinas, SP, Brazil. His research
The presented small-signal model is useful to analyze the
interests are distributed generation and power systems stability and control. (e-
impact of voltage positive-feedback anti-islanding schemes on mail: [email protected]).
grid-connected IBDG systems stability. The linearized set of Diogo Salles is Ph.D. graduate student with the Department of Electrical
equations can well represent the system against small Energy Systems of University of Campinas, Campinas, SP, Brazil. His research
perturbations, with the advantages of the state-space interests are distributed generation, distribution systems and power systems
protection, stability and control. (e-mail: [email protected]).
representation. This state-space representation allows the use of
a whole set of linear control techniques to direct stability Walmir Freitas is Associate Professor with the Department of Electrical
assessment, without the need of slow time domain simulations. Energy Systems of University of Campinas, Campinas, SP, Brazil. His research
interests include distributed generation and distribution systems. (e-mail:
[email protected]).
Through the modal analysis of the DAE representation of the
system one can conduct an extensive parametric sensitivity Xiaoyu Wang is Assistant Researcher with the Department of Electrical
study in order to assess the limit penetration of multiple IBDGs Engineering at Tsinghua University, Beijing, China. His research interests
in distribution systems. include distributed generation and power system stability. (e-mail:
[email protected]).
References
[1] N. Jenkins, R. Allan, P. Crossley, D. Kirschen e G. Strbac, “Embedded
Generation”, London: The Institution of Electrical Engineers (IEE), 2000,
292 p.
[2] A. M. Borbely e J. F. Kreider, "Distributed Generation: The Power
Paradigm for the New Millennium", Boca Raton, FL: CRC Press, 2001,
416 p.
[3] IEEE Standard for Interconnecting Distributed Resources with Electric
Power Systems, IEEE Standard 1547-2003, New York, Jun 2003.
[4] IEEE Recommended Practice for Utility Interface of Photovoltaic (PV)
Systems, IEEE Standard 929-2000, New York, Jan 2000.
[5] W. Bower e M. Ropp, “Evaluation of Islanding Detection Methods for
Utility-Interactive Inverters in Photovoltaic Systems” Sandia National
Laboratories, Albuquerque, NM, Tech. Rep. SAND 2002-3591, Nov.
2002.
[6] J. Stevens, R. Bonn, J. Ginn, S. Gonzalez e G. Kern, “Development and
Testing of an Approach to Anti-Islanding in Utility-Interconnected
Photovoltaic Systems” Sandia National Laboratories, Albuquerque, NM,
Tech. Rep. SAND 2000-1939, Aug. 2000.
[7] Z. Ye, R. Walling, L. Garces, R. Zhou, L. Li e T. Wang, "Study and
Development of Anti-Islanding Control for Grid-Connected Inverters"
National Renewable Energy Laboratory (NREL), Golden, CO, Tech.
Rep. NREL/SR-560-36243, May. 2004.
[8] X. Wang, “Investigation of Positive Feedback Anti-Islanding Scheme for
Inverter-Based Distributed Generation” Ph.D. thesis, Dept. Elect. and
Comp. Eng., Univ. Alberta, Edmonton, 2008.
[9] W. Freitas, “Métodos para Análise de Proteção Anti-ilhamento de
Geradores Distribuídos” Thesis, School of Elect. and Comp. Eng., Univ.
of Campinas, 2007.
[10] X. Wang, W. Freitas, W. Xu e V. Dinavahi, "Impact of DG Interface
Controls on the Sandia Frequency Shift Antiislanding Method" IEEE
Transactions on Energy Conversion, vol.22, no.3, pp.792-794, Sep. 2007.
[11] X. Wang e W. Freitas, "Impact of Positive-Feedback Anti-Islanding
Methods on Small-Signal Stability of Inverter-Based Distributed
Generation" IEEE Transactions on Energy Conversion, vol.23, no.3,
pp.923-931, Sep. 2008.
[12] X. Wang, W. Freitas, W. Xu e V. Dinavahi, "Impact of Interface Controls
on the Steady-State Stability of Inverter-Based Distributed Generators" in