Symmetrical and Unsymmetrical Fault Currents of A Wind Power Plant
Symmetrical and Unsymmetrical Fault Currents of A Wind Power Plant
Conference Paper
NREL/CP-5500-53463
December 2011
Contract No. DE-AC36-08GO28308
NOTICE
The submitted manuscript has been offered by an employee of the Alliance for Sustainable Energy, LLC
(Alliance), a contractor of the US Government under Contract No. DE-AC36-08GO28308. Accordingly, the US
Government and Alliance retain a nonexclusive royalty-free license to publish or reproduce the published form of
this contribution, or allow others to do so, for US Government purposes.
This report was prepared as an account of work sponsored by an agency of the United States government.
Neither the United States government nor any agency thereof, nor any of their employees, makes any warranty,
express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of
any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately
owned rights. Reference herein to any specific commercial product, process, or service by trade name,
trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation,
or favoring by the United States government or any agency thereof. The views and opinions of authors
expressed herein do not necessarily state or reflect those of the United States government or any agency thereof.
Available electronically at http://www.osti.gov/bridge
Available for a processing fee to U.S. Department of Energy
and its contractors, in paper, from:
U.S. Department of Energy
Office of Scientific and Technical Information
P.O. Box 62
Oak Ridge, TN 37831-0062
phone: 865.576.8401
fax: 865.576.5728
email: mailto:[email protected]
Available for sale to the public, in paper, from:
U.S. Department of Commerce
National Technical Information Service
5285 Port Royal Road
Springfield, VA 22161
phone: 800.553.6847
fax: 703.605.6900
email: [email protected]
online ordering: http://www.ntis.gov/help/ordermethods.aspx
Cover Photos: (left to right) PIX 16416, PIX 17423, PIX 16560, PIX 17613, PIX 17436, PIX 17721
Printed on paper containing at least 50% wastepaper, including 10% post consumer waste.
Member IEEE
[email protected]
M. Singh
Member IEEE
[email protected]
E. Muljadi
Fellow, IEEE
[email protected]
I. INTRODUCTION
2
3
As shown in Figure 4, the fault current is driven by the decaying flux trapped in the rotor winding as represented by the
right portion of equation (1). The larger the leakage inductances (), the smaller is the fault current amplitude. The fault
current dies out after the flux driving the fault current depleted to zero. Note that the DC and AC transient components of
the SCC flowing out of the stator windings induce fault currents in the rotor winding and vice versa until the magnetic
flux is depleted.
Min Isc-peak
Max Isc-peak
WTG
Max
ISC_PEAK
Min
ISC_PEAK
Type 1
2
Type 2
2V s
X
'
S
2V s
X
'
S
2V s
2Vs
X S'
' 2
S
+ (9 Rr' ) 2
Type 3
2
2Vs
X S'
1.1
IRATED
Type 4
1.1
IRATED
For a Type 1 WTG, the maximum SCC is based on the assumption that the DC offset is at the worst condition, and the
minimum SCC is calculated by assuming that the DC offset is
zero. For a Type 2 WTG, the maximum value is computed when
= 0 . The minimum value is computed when the slip
reaches 10% above synchronous speed. And for a Type 3 WTG,
the maximum value is computed when the crowbar shorts the
rotor winding and the minimum value is computed when the
power converter can follow the commanded current (i.e., in case
the fault occurs far away from the point of interconnection (POI),
the remaining terminal voltage is sufficiently high enough to let
the power converter operate normally and supply the commanded currents). Note that for a symmetrical fault, the actual fault
current for each phase is different from the other phases due to
the fact that the time of the fault occurs at a different phase angle
for different phases, thus affecting the DC offset. For a Type 4
WTG, the stator current can always be controlled because of the
nature of power converter, which is based on a current-controlled
voltage source converter.
A time domain simulation is performed in PSCAD, and the
steady-state calculations are performed using Mathcad and
Cyme software for a symmetrical fault. The results are tabulated in Table II. The calculated results from different software platforms are very close to the approximation listed in
Table I. Note, that only Type 1 and Type 4 are listed. The
Type 2 and Type 3 WTGs will respond differently because of
the existence of the external rotor resistance in a Type 2 WTG
and the activation of the crowbar circuit in a Type 3 WTG,
which will respond non-linearly to the fault. The SCC for a
Type 2 and Type 3 WTG, as indicated in Table I, will have
the size difference between the SCC of the Type 1 and Type 4
WTGs.
a)
Type 1 WTG.
b)
Type 2 WTG.
c)
Type 3 WTG.
d)
Type 4 WTG.
TABLE II
ISC_PEAK COMPARISON FOR DIFFERENT SOFTWARE PLATFORMS
WTG
Table I
PSCAD
Cyme
Mathcad
Type
Min
Max
3.4 p.u.
6.3
p.u.
5 p.u.
5.5 p.u.
3.8 p.u.
0 p.u.
1.1
p.u.
1.1 p.u.
1.1 p.u.
1.1 p.u.
In Figure 13, the SCC for a Type 1 WTG is shown both for
the three-phase currents and the corresponding sequence
components. The changes in positive sequence and the sudden appearance of the negative sequence are also shown. The
absence of the zero sequence current is a consequence of
winding connections.
Type 4
I2-
I2o
I6+
I6-
I6o
1.40
1.20
1.00
Other
Types
I2+
0.80
0.60
0.40
0.20
0.00
20
1.20
10
1.00
I3a
0.80
-10
I3b
0.60
-20
I3c
0.40
1.05
1.1
1.15
1.2
1.25
1.3
1.35
1.4
0.20
TIME (sec)
0.00
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
In Figure 14, the SCC for a Type 4 WTG shows the fault
currents in its sequence current components. At the POI
(Figure 14a), there exist both the zero sequence and the negative sequence currents because of the substation transformer
winding connection (YgYg) and collector system capacitances
respectively. As shown in Figure 14b, at the generator terminals however, there is a pad-mounted transformer (Yg)
that will block the zero sequence component, and the Type 4
WTG produces a positive sequence component (refer to neg.
sequence switch in the equivalent circuit shown Figure 12).
In Table III, the SCC at the POI is computed for a SLG
fault using different software platform. It is shown that the
SCC results match for all software for Type 1 WTG. A small
mismatch between the MathCAD and PSCAD results can be
seen in Table 3. This discrepancy exists because in
MathCAD, we remove the capacitances from the circuit.
TABLE III
COMPUTED SCC AT THE POINT OF INTERCONNECTION OF THE WPP
IN SEQUENCE COMPONENTS IN PER UNIT
FOR AN SLG FAULT WITH DIFFERENT SOFTWARE PLATFORMS
WTG
Type 1
Type 4
I_Seq
PSCAD
Cyme
Math-cad
1.42
1.58
1.67
0.58
0.57
0.55
0.63
0.57
0.65
1.4
1.44
1.35
1.0
0.95
1.0
0.25
0.0
0.0
TABLE IV
COMPUTED SCC AT THE GENERATOR TERMINALS
IN SEQUENCE COMPONENTS IN PER UNIT
FOR AN SLG FAULT WITH DIFFERENT SOFTWARE PLATFORMS
WTG
Type 1
Type 4
WTG
I_Seq
PSCAD
Cyme
Mathcad
0.0
0.0
0.0
0.95
0.86
0.95
I_Seq
PSCAD
Cyme
MathCAD
0.0
0.0
0.0
0.8
0.86
1.05
0.58
0.56
0.55
0.0
0.0
0.0
0.62
0.59
0.65
1.0
1.0
1.0
0.0
0.0
0.0
0.4
0.0
0.0
1.0
1.0
1.0
0.0
0.0
0.0
Type 1
Type 4
In Table VI, the SCC at the POI is presented for LLG fault.
It is shown that we have a very good match between PSCAD,
Cyme, and MathCAD calculations. In comparison to Table
V, we can see the presence of the zero sequence current in
LLG fault as expected in Table VI.
TABLE VI
COMPUTED SCC AT THE POINT OF INTERCONNECTION OF THE WPP
IN SEQUENCE COMPONENTS IN PER UNIT
FOR AN LLG FAULT WITH DIFFERENT SOFTWARE PLATFORMS
WTG
Type 1
Figure 15: An example of typical output panel for SCC calculation with
Cyme.
Type 4
I_Seq
PSCAD
Cyme
Mathcad
1.55
1.63
1.48
1.17
1.11
1.1
0.55
0.62
0.8
1.4
1.49
1.7
0.9
0.25
0.0
the corresponding sequence components are presented in Figure 17. As shown in Figure 16, it is difficult to discern the
type of faults that occur in the line. In comparison, from Figure 17, it is obvious there is a distinction between the fault
currents for the LL fault and the LLG fault.
IV. CONCLUSIONS
In this paper, the SCC contributions of different WTGs for
faults at the terminal of the generator were simulated using
time domain simulations and steady-state calculations. Two
power system commercial software platforms were used
(PSCADTM, and CymeTM), and a multipurpose mathematical
computer program (MathCADTM) is also used to compute the
SCC.
A simplified method to compute the SCC for a symmetrical fault is presented and it is tabulated in Table I. The SCC
results were tabulated in Table II, comparing the size of the
SCC at the POI for three different methods. Note that only
Type 1 and Type 4 WTGs are used because they represent the
maximum and minimum SCC contribution. The calculations
for Type 2 and Type 3 WTGs are affected by time-of-fault
occurrence and the action of the external rotor resistance control and the crowbar action, thus, the SCC contribution is usually lower than the Type 1 WTG, but it is higher than the
Type 4 WTG.
The unsymmetrical faults were simulated and tabulated.
As shown in the Table III through Table V, the unsymmetrical fault calculations from three different software packages
shows a good agreement for unsymmetrical fault calculations.
Each WPP is unique. Therefore, recommended practice
from local reliability organizations, manufacturers, transmission planners, wind plant developers, and local utilities should
be followed very closely when performing studies of WPP.
Main : Graphs
1.25
I SCC POI
1.00
0.75
0.50
0.25
0.00
-0.25
-0.50
-0.75
-1.00
sec
1.00
1.20
a)
1.40
1.60
1.80
2.00
1.80
2.00
1.25
I SCC POI
1.00
0.75
0.50
0.25
0.00
-0.25
-0.50
-0.75
-1.00
sec
1.00
1.20
1.40
1.60
V. ACKNOWLEDGMENT
Main : Graphs
0.300
I2+
I2-
I2o
0.250
0.200
0.150
VI. REFERENCES
0.100
[1]
0.050
0.000
x
1.00
1.20
a)
1.40
1.60
1.80
2.00
[2]
0.450
I2+
I2-
[3]
I2o
0.400
0.350
0.300
[4]
0.250
0.200
0.150
0.100
[5]
0.050
0.000
x
1.00
1.20
1.40
1.60
1.80
2.00
[6]
[7]
U.S. Department of Energy Energy Efficiency and Renewable Energy, 20% Wind Energy by 2030 Increasing Wind Energys Contribution to U.S. Electricity Supply, May, 2008.
J. Charles Smith, Michael R. Milligan, Edgar A. DeMeo and Brian
Parsons, "Utility wind Integration and operating impact state of the art,"
IEEE Trans. Power Systems, vol. 22, pp. 900-908, Aug. 2007.
IEEE PES Wind Plant Collector System Design Working Group, Wind
Power Plant Grounding, Overvoltage Protection, and Insulation Coordination, Proceedings of the 2009 IEEE Power and Energy Society
General Meeting.
Nader Samaan, Robert Zavadil, J. Charles Smith and Jose Conto,
Modeling of Wind Power Plants for Short Circuit Analysis in the
Transmission Network, in Proc. of IEEE/PES Transmission and Distribution Conference, Chicago, USA, April 2008.
J. Moren, S.W.H. de Haan, Short-Circuit Current of Wind Turbines
with Doubly Fed Induction Generator, IEEE Transactions on Energy
Conversion, Vol. 22, No. 1, March 2007.
Sulawa, Zabara, et al. Short circuit current of induction generators.
IEEE ISCAS 2007 proceedings.
E. Muljadi, V. Gevorgian, Short Circuit Modeling of a Wind Power
Plant, in Proc. 2011 IEEE Power Engineering Society General Meeting.
VII. BIOGRAPHIES
Vahan Gevorgian (M97) graduated from the Yerevan Polytechnic Institute (Armenia) in 1986. During
his studies he concentrated on electrical machines. His
thesis research dealt with doubly-fed induction generators for stand-alone power Systems. He obtained his
Ph.D. degree in Electrical Engineering Dept. from the
State Engineering University of Armenia in 1993. His
dissertation was devoted to a modeling of electrical
transients in large wind turbine generators.
Dr. Gevorgian is currently working at the National
Wind Technology Center (NWTC) of National Renewable Energy Laboratory (NREL) in Golden, Colorado, USA, as a research engineer. His current
interests include modeling and testing of various applications of small wind
turbine based power systems.
Mohit Singh (M2011) received his M.S. and Ph.D.
in Electrical Engineering from the University of Texas,
Austin in 2007 and 2011 respectively. His research is
focused on dynamic modeling of wind turbine generators.
Dr. Singh is currently working at the National Renewable Energy Laboratory (NREL) in Golden, Colorado, USA, as a post-doctoral researcher in transmission and grid integration of renewable energy. His
current interests include modeling and testing of various applications of wind
turbine generators and other renewable energy resources. He is a member of
the IEEE. He is involved in the activities of the IEEE Power and Energy
Society (PES).
Eduard Muljadi (M82-SM94-F10) received
his Ph. D. (in Electrical Engineering) from the University of Wisconsin, Madison. From 1988 to 1992,
he taught at California State University, Fresno, CA.
In June 1992, he joined the National Renewable
Energy Laboratory in Golden, Colorado. His current
research interests are in the fields of electric machines, power electronics, and power systems in
general with emphasis on renewable energy applications. He is member of
Eta Kappa Nu, Sigma Xi and a Fellow of the IEEE. He is involved in the
activities of the IEEE Industry Application Society (IAS), Power Electronics
Society, and Power and Energy Society (PES).
He is currently a member of various committees of the IAS, and a member
of Working Group on Renewable Technologies and Dynamic Performance
Wind Generation Task Force of the PES. He holds two patents in power
conversion for renewable energy.