Study of ShortCircuit Currents Around D？tmarovice Power Station
 Author: Ali Shehab Abdulwadood
 Organization: Ali Shehab Abdulwadood
 Publish: Transactions on Electrical and Electronic Materials Volume 15, Issue3, p117~124, 25 June 2014

ABSTRACT
The calculation of shortcircuit currents is important for power systems operation and restoration, and for determining the means to protect human lives and properties. In this paper, a part of a power system network, around the Dĕtmarovice power station in Czech Republic, was simulated by the well known program EMTPATPDraw (Electromagnetic Transients ProgramAlternative Transient Program), and shortcircuit currents and voltages were calculated at different points in the electric network and presented as a time function by the PlotXY program. Calculations were done just for phasetoground, and for the threephase shortcircuit at the Kunčice substation. The results were important for determining the characteristics of the equipment required to withstand or break the shortcircuit current; for this reason, the calculations were repeated using earthfault resistances only for the case of busbar KUN shown in
Figs. 5 and6 .

KEYWORD
Shortcircuit , Shortcircuit currents and voltages , Phasetoground , Threephase , EMTPATPDraw

1. INTRODUCTION
Normally, short circuit studies involve calculating a bolted threephase fault condition, which can be characterized as all three phases 'bolted' together to create a zero impedance connection. This establishes the 'worst case' (highest current) condition that results in maximum threephase thermal and mechanical stress in the system. From this calculation, other types of fault conditions can be approximated. This “worst case” condition should be used to specify proper interrupting rating requirements, selectively coordinate the system and provide component protection [4].
In a threephase system, a fault may equally involve all three phases. A bolted fault acts as if the three phases were connected together with links of zero impedance prior to the fault, i.e., the fault impedance itself is zero and the fault is limited by the system and machine impedances only. Such a fault is called a symmetrical threephase bolted fault, or a solid fault.
Generally, such faults give the maximum shortcircuit currents and form the basis of the calculations of shortcircuit duties on switching devices. Faults involving one, or more than one, phase and ground are called unsymmetrical faults. Unsymmetrical faults are more common than threephase faults, but under certain conditions, linetoground fault or double linetoground fault current values may exceed threephase symmetrical fault currents. Shortcircuit calculations are, thus, the primary study whenever a new power system is designed or an existing system is to be expanded or upgraded [7].
In a distribution system, several variables affect calculated bolted threephase shortcircuit currents. Appropriate variable values must be selected for the specific application analysis. Electrical installations almost always require protection against shortcircuits wherever there is an electrical discontinuity. This most often corresponds to points where there is a change in conductor crosssection. The shortcircuit current must be calculated at each stage of the installation to determine the characteristics of the equipment required to withstand or break the fault current [4].
Sources of a shortcircuit current that are normally considered include utility sources, local sources, synchronous motors, induction motors and alternate power sources [4].
The consequences of shortcircuits can vary depending on the type and duration of the fault, the installation point (as in a physical location) where the fault occurred and the shortcircuit power. Consequences include:
 At the fault location: the presence of electrical arcs, resulting in damage to insulation, welding of conductors, initiation of fires and endangerment to life.
 On the faulty circuit: electrodynamic forces, resulting in the deformation of the busbars, disconnection of cables, and excessive rise in temperature due to an increase in Joule losses, risking damage to insulation.
 On other circuits in the network or in nearby networks: voltage drop during the time required to clear the fault, ranging from a few milliseconds to a few hundred milliseconds.
 Dynamic instability and/or the loss of machine synchronization, etc.
Shortcircuit calculations should be done at all critical points in the system. These would include the service entrance, transfer switches, panel boards, load centers, motor control centers, disconnects and motor starters. The various types of shortcircuits that can occur during electrical installations are as follows [3]:
 Phasetoground (80% of faults).
 Phasetophase (15% of faults). This fault often degenerates into a 3phase fault.
 Threephase (only 5% of initial faults).
In networks with isolated nodes, or ineffectively earthed, the connection of the phase with the earth dose not called singlephase shortcircuit, but qualitatively a ground fault.
In this study, only phasetoground and threephase shortcircuits are calculated to investigate shortcircuit currents and voltages.
From Figure 1a, the characteristic equations for a singlephase short circuit can be written directly as follows [12]:
In a real singlephase short circuit, the voltage at the damaged phase does not decrease to zero. This type of short circuit can occur only in networks with earthed nodes of transformers or generators. For shortcircuits, phase A current can be written as follows:
where
E is the phase voltage.
Z_{1}, Z_{2} and Z_{0} are the positive, negative and zero impedances, respectively.
The characteristic equations for a threephase short circuit is as follows (Fig. 1(b)):
For the phase currents, assuming a symmetrical short:
where
Switching from the state before the fault to the new state is accompanied by a transient. In singlephase shortcircuits (Fig. 2(a)), the transient is caused by the balance capacities in phases B and C, where its energy shifts to inductance, which will be transformed into the energy of the coil’s magnetic field and vice versa. This results in an oscillating circuit, whose natural frequency depends on the values of capacities and inductance and oscillates the phase A voltage from the phase value U_{A} to, or near, zero and U_{B} and U_{C} to the values U'_{B} and U'_{C} calculated by equations 5 and 6 below:
The diagram of Fig. 2(a) can be replaced by Fig. 2(b), and the frequency of the circuit is calculated by equation 7 below:
The results, obtained by ATPDraw [2], which is a graphical, mousedriven preprocessor to the ATP version of the Electromagnetic Transients Program (EMTP), are drawn by PlotXY program [11].
2. THE STUDIED CASE
The studied case (Fig. 3) corresponds to a part of a power network around a branch from the Dĕtmarovice power station (EDE)  block No. 2 in Czech Republic. Block No. 2 (EDE2) has a generator with a threephase voltampere rating of about SG=235 MVA and a voltage of about 15.75 kV and it supplies to the power network 110 kV a power of about 160 MW via the block transformer (EDET). A part of this power, about 100 MW, is consumed in the Kunčice substation (KUN), and 60 MW is fed to the supergrid 220 kV of Lístkovec (LIS2) via the distribution transformer (LIST). Both parts of the power are transferred via the Vratimov substation (VRA), considered to be under no load. All the substations are connected by AlFe110 kV transmission lines. A single line, No. 692 (two conductors in a bundle 2×670 mm^{2}, length 20.4 km), is connected between the Dĕtmarovice power station (EDE) and the Vratimov substation (VRA); a double line No. 605606 (240 mm^{2}, length 5.168 km) between the Vratimov substation (VRA) and the Kunčice substation (KUN); and a double line No. 641642 (450 mm^{2}, length 8.107 km) between the Vratimov substation (VRA) and the Lístkovec substation (LIS).
3. CONSTRUCTION OF THE MODELS
The power system has two supply stations, the block of Dĕtmarovice (EDE), which contains a synchronous generator (EDE2) and a block transformer (EDET), and the supergrid 220 kV of Lístkovec (LIS2), in addition to two transformers, three transmission lines and one load of Kunčice (KUN). ATPDraw has sufficient models to represent this power system, but some components are needed to calculate shortcircuit currents and voltages, such as measuring probes and a combination of switches with resistances.
3.1 Supply stations
The synchronous generator (SG=235 MVA/15.75 kV) of the Dĕtmarovice supply stations (EDE) is represented by the ATPDraw model SM59. Parameters of the generator are listed in Table 1 below. The voltage magnitude is calculated as in equation 8 [14] and represents the steadystate voltage at the terminals of the machine.
The supergrid 220 kV of Lístkovec (LIS2) is represented as an ideal source of voltage (AC3ph type 14) with amplitude given by equation 9:
3.2 Transformers
The power system contains two transformers, the block transformer (EDET), type EC 41 M, manufactured by Škoda PlzeňCzech Republic, and the distribution transformer (LIST), manufactured by Škoda Plzeň. Both transformers are modeled by the ATPDraw model BCTRAN. Table 2 describes the parameters of the transformers.
The other values are as follows: three phases, two windings, a 3legged stacked core, test frequency of 50 Hz and DYn1 connection for EDET and YNyn0 connection for LIST.
3.3 Transmission lines
The 110 kV transmission lines, made of AlFe, are modeled by the ATPDraw model Lines/Cables (LCC). The parameters and location of the conductors on the towers for each line are listed below:
 The single line No. 692 has two conductors in a bundle: cross sectional area of 2×670 mm2, length of 20.4 km, an overhead line, 3phase, type PI, with skin effect, Rho (ground resistivity) = 120 Ωm, and initial frequency = 50 Hz. The parameters and location of the conductors on the tower are given in Table 3.
 Double lines No. 605606 are made of AlFe: cross sectional area of 240 mm2, length 5.168 km, overhead lines 3phase, type PI, with skin effect, Rho (ground resistivity)=100 Ωm, and initial frequency = 50 Hz. The parameters and the locations of the conductors on the tower are given in Table 4.
 Double lines No. 641642 are made of AlFe: cross sectional area of 450 mm2, length 8.107 km, overhead line, 6phase, type PI, with skin effect, Rho (ground resistivity)=100 Ωm, and initial frequency = 50 Hz. The parameters and the location of the conductors on the tower are given in Table 5.
3.4 Loads
The load in the Kunčice substation (KUN) can be easily modeled by a standard component RLC3 as a resistive load, connected to a star with a grounded node. Load consumed is about 100 MW, with 110 kV. The resistance value will be [14]:
Finally, the scheme of the Dĕtmarovice power station  block No. 2 shown in Fig. 3, will be modeled by ATPDraw as shown in Fig. 4 below.
4. RESULTS AND DISCUSSION
In this study, we investigate the shortcircuit currents and voltages at the bus of the Kunčice substation (KUN), as well as at some points such as the secondary side of the Dĕtmarovice transformer (EDET) and at the secondary side of the Lístkovec transformer (LIST). The shortcircuit calculations will be done only for phasetoground and for threephase at the Kunčice substation (KUN).
For each analysis in this study, 200 statistical simulations were performed in total. A switching simulation was represented by a statistical switch, modeled at ATP as SW_STAT, with a Gaussian probability distribution and described by a standard deviation σ (set to 1 μs). For all simulations, the starting time of the fault was set to 0.0378 s, which is the time the voltage of phase A crosses the maximum.
4.1 Phasetoground shortcircuit at the Kun？ice substation
Assuming that the fault occurs at phase A, the following figures illustrate the imbalance situations that could occur at the different points in the network.
4.1.1 Busbar KUN
Figures 5 and 6 show the currents of the three phases of the load of Kunčice and the voltages at this busbar. The current and the voltage at the damaged phase A (○) decreased to zero, while those at the faultless phases B (□) and C (△) were slightly affected. The amplitude of the current at all phases (Fig. 5), in the case of steady state, reached 752 A and because the fault, the amplitude of the current at phases B and C increased to almost 811.5 A and 861.4 A, respectively (increased about 7.9% and 14.5%, respectively). In Fig. 6, the amplitude of the voltage reached 88.7 kV at steady state and because the fault increased to 95.7 kV (7.89%) at phase B and 101.65 kV (14.6%) at phase C.
4.1.2 At the secondary side of the D？tmarovice transformer EDET
The amplitudes of the three phase steady state currents that flowed from the secondary side of the transformer EDET (Fig. 7) are 1,164.1 A at phase A (○), 1,207.6 A at B (□) and 1,184 A at C (△).because the fault, the amplitude of the phase A current increased to a maximum 3,794.2 A (225.9%). This disturbance appears in the curves of voltages (Fig. 8), whose the amplitudes at the three phase steady state were almost 89.85 kV, because the fault decreased to 51.1 kV at phase A (decreased 43.1%).
4.1.3. At the secondary side of the Listkovec transformer list
The amplitudes of the three phase steady state currents that flowed from LIST (Fig. 9) are 441.23 A at phase A (○), 451.59 A at B (□) and 418.36 A at C (△).because the fault the amplitude of the phases A current increased to a maximum of 8,225.3 A (1,764.2%). The amplitudes of the three phase steady state voltages (Figure 10) are almost 89.1 kV because the fault decreased to 41.31 kV (53.6%) at phase A.
4.2 Threephase shortcircuit at the Kun？ice substation
Based on the above explanations, the following imbalance situations occur:
4.2.1 Busbar KUN
Figures 11 and 12 show the currents and voltages at the three phases. The current and the voltage at the damaged phases decreased to almost zero. The amplitude of the current at all phases (Fig. 11) at steady state reached 753 A and the amplitude of the voltage reached 89 kV (Fig. 12).
4.2.2. At the secondary side of the D？tmarovice transformer EDET
The amplitudes of the steady state currents at all phases (Fig. 13) are 1,164.1 A at phase A (○), 1207.6 A at B (□) and 1,184.1 A at C (△).because the fault, the amplitudes of the currents increased to a maximum 4,219.5 A at phase A (262.5%), 7,135.6 A at B (490.9%) and 6,979.4 A at C (489.4%). This disturbance appears in the curves of voltages (Fig. 14), whose amplitudes of the three phase steady state are almost 90 kV, because the fault decreased to almost 40 kV (55.6%) in all phases.
The transient disturbance at the start of the fault showed that the transient voltage at phase A reached the value of 131.16 kV, i.e. increased by about 45.7% in comparison with that at the steady state.
4.2.3 At the secondary side of the Listkovec transformer LIST
The amplitudes of the three phase steady state currents that flowed from the LIST (Fig. 15) are 443.71 A at phase A (○), 451.71 A at B (□) and 429.37 A at C (△).because the fault, the amplitude of the current increased to a maximum 10,077 A (2,171.1%) at phase A, 16,019 A (3,446.3%) at B and 16,090 A (3,647.4%) at C. The amplitudes of the three phase steady state voltages (Fig. 16) are almost 89.2 kV because the fault decreased to almost 33.5 kV (62.4%) in all phases.
At the start of the fault, the transient voltage reached the value of 112.7 kV (26.3%) at phase A.
5. NOTES
A fault such as indirect connection of one phase with the ground in networks with an isolated node or in compensated networks (node is connected via arcsuppression coil) is called an earthfault (or ground fault). This fault may occur in normal operation and currents will leak from phase to earth. Earthleakage currents have small values and are capacitive [12], while shortcircuit currents usually reach multiple values higher than those at steady state operation and are inductive.
In this study, we assumed that the earthfault resistance was zero, although it has a decisive influence on the Yvoltage (U0) value and thereby on earthfault protections. For example, the functions of some earthfault protections, such as GSC12 [12], are blocked by the value of (U0); protection will not react if the value does not reached 30~40% of the nominal value.
Shortcircuit calculations are possible using earthfault resistances if the power system is correctly modeled and the study aims to find the correct configuration for the network to avoid faults. In ATPDraw, the earthfault resistance is simply modeled by the component "Resistor", which will be onepole or threepole grounded. For example, with earthfault resistance of 10 Ω for the scheme corresponding to Figs. 5 and 6, the curves of the currents and voltages will differ as shown below. The amplitude of the current at phase A reached 538.92 A, not zero (Fig. 5), and the amplitude of the voltage reached 63.6 kV.
6. CONCLUSIONS
This study demonstrated the imbalance situations due to short circuits at different locations in power systems. Shortcircuits may occur in well designed power systems and produce large decaying transient currents, in most cases much above the system load currents. Such a case can be studied to identify the worst possible cases in a power system by modeling using a suitable program, such as EMTPATPDraw, and simulating the faults at different locations in the network.
The high oscillations on the voltage curves due to a fault originate from the oscillating circuit described previously in Fig. 2.
The worst case in this simulation (the fault started when the voltage of phase A crosses the maximum, t=0.0378 s), which need to be protected against, was the rise of a short circuit current at the secondary side of the Listkovec transformer LIST for both phasetoground and threephase faults. Referring to Fig. 9 and Fig. 15, the amplitude of the current at phase A increased to a maximum 8,225.3 A (1,764.2%) because of the phasetoground fault and increased to 16,090 A (3,647.4%) at phase C because of the threephase fault. These values decreased when the voltage was not the maximum.
The worst case was the threephase fault, which is called a symmetrical threephase bolted fault. Knowing this, the shortcircuit to the faulty section of the electrical system should be limited by applying appropriate switching devices that can operate under shortcircuit conditions without damage and isolate only the faulty section.

[Fig. 1.] Singlephase and threephase shortcircuits. (a) singlephase and (b) threephase.

[]

[]

[]

[]

[Fig. 2.] Network diagram and directions of energy flow. (a) network diagram and (b) equivalent circuit.

[]

[]

[]

[Fig. 3.] Scheme of the D？tmarovice power station  block No. 2.

[Table 1.] The parameters of the generator.

[]

[]

[Table. 2.] Parameters of BCTRAN.

[Table 3.] The parameters of the transmission line No. 692.

[Table. 4.] The parameters of transmission line No. 605606.

[Table 5.] The cable input data.

[]

[Fig. 4.] The ATPDraw model of the network.

[Fig. 5.] The currents to the load KUN.

[Fig. 6.] The voltages at the busbar KUN.

[Fig. 7.] The currents at EDET.

[Fig. 8.] The voltages at EDET.

[Fig. 9.] The voltages at EDET.

[Fig. 10.] The voltages at EDET.

[Fig. 11.] The currents to the load KUN.

[Fig. 12.] The voltages to the load KUN.

[Fig. 13.] The currents at EDET.

[Fig. 14.] The voltages at EDET.

[Fig. 15.] The currents at LIST.

[Fig. 16.] The voltages at LIST.

[Fig. 17.] The currents to the load KUN  (□) with and (○) without earthfault resistance.

[Fig. 18.] The voltages at the busbar KUN  (□) with and (○) without earthfault resistance.