Study of Impulse Voltage Distribution in Transformer Windings

Download as pdf or txt
Download as pdf or txt
You are on page 1of 5

2010 IEEE International Conference on Power and Energy (PECon2010), Nov 29 - Dec 1, 2010, Kuala Lumpur, Malaysia

Study of Impulse Voltage Distribution in


Transformer Windings
Hendri Masdi *, Norman Mariun *, Azah Mohamed **and Noor Izzri Abdul Wahab*
* Department of Electrical and Electronic Engineering, Faculty of Engineering,
University Putra Malaysia 43400 Serdang, Selangor Darul Ehsan, Malaysia.
** Department of Electrical, Electronic and Systems Engineering, Faculty of Engineering,
Universiti Kebangsaan Malaysia, 43600 UKM, Bangi, Selangor, Malaysia.
AbstractIn this paper, analysis impulse voltage distribution
in transformer winding, and winding of the voltage impulse
response with a standard full wave and clipped wave especially
at
the
ends
of
the
entanglement
presented.
Compared with experimental data EMTPDC simulation, the
result of t <1.0 s voltage difference occurs. For t > 2.0 s was
obtained near the same voltage large, between experimental
and simulation results. So we can say that circuit simulation
can
mimic
the
real
situation
EMTDC/PSCAD.
Construction of turns and the thickness of insulation (Cs and
Cg) transformer, poorer have considerable effect on the
distribution center of experiencing transient voltage
transformer(transitional-symptoms). Winding section closest
to the high voltage high voltage gradient than most parts of
another entanglement, tests were performed in a standard
form according to the current surge of IEC-71 8/20s.
Important information from this study is that the simulation
gives a more pessimistic outcome (EMTDC/PSCAD
calculations show a larger stress values compared with
experimental). Although it is difficult to imitate the behavior
of entanglement in the initial transient, but the real situation
is not worse. Thus, this simulation can be used as EMTPDC
experimental substitution want to know the response to the
current surge of turns.
KeywordsEMTDC/PSCAD; transformer model; Transient
overvoltage.

II.

I.
INTRODUCTION
Analysis impulse voltage distribution in transformer
winding, and winding of the voltage impulse response
with a standard full wave and clipped wave especially at
the ends of the entanglement. For transient study in
electric power system can nowadays be done by
simulation with the help of software packages dedicated
for that purpose, such as EMTDC/PSCAD and others.
Studies that involve substation require a good skill and
knowledge in how to model more or less correctly the
equipment included in substations model. There has been
a lot of discussion regarding transient model of
substations main equipment, such as: power transformer,
Surge Arrester and others.
So far, the transient model of a power transformer
consists of a circuit of capacitance that represent the
capacitance of the winding insulation to ground, and the
inter-turn capacitance. To model properly a power
transformer that could reproduce the phenomenon close to

978-1-4244-8946-6/10/$26.00 2010 IEEE

real situation is very important, especially if one wishes to


study the stresses borne by the winding during occurrence
of transient overvoltage [1,2,3].
To realize that goal, we experiments and simulation
were conducted on a distribution power transformer
whose capacity was: 100kVA, 20.000V/400V. To
experimentally simulate a lightning current Surge Current
Generator from Schaffner Type NSG 857 was utilized,
which can produce surge current up to 7,5 kA of 8x20
microsecond.
In the experiment, current of 5 kA was injected to the
transformer winding, and then the over voltage between
the transformer HV-terminal and ground and also between
the transformer conductor layer to ground were measured.
From the measurement, a waveform and voltage
distribution along the winding can then be obtained [6].
Further, simulations were conducted, by applying the
same current waveform. The simulation was performed
using the EMTDC/PSCAD; and the transformer constants
included in the simulation have been derived from the
physical dimension of the transformer winding [7].

379

TRANSFORMER CONSTRUCTION AND ITS


CONSTANTS

Schematic drawing of the transformer investigated in


the study can be seen in figure 1, where the HV-winding
was constructed of 9 layers. The transformer was a three
phase winding. In the experimental simulation, the surge
current was injected to the tap 1 and the tap 9 connected
to ground. These taps represent the HV-terminal and the
neutral point of the transformer, respectively.

Connector
(a)
RSTN

(a)

In the simulation the following assumptions have also


been adopted:

The HV-winding under test was considered as


isolated from the other two windings and from
the LV-winding.
The transformer was assumed as air-core
transformer. This assumption was applied since
the influence of iron only significant for the case
of switching impulse.

In the study, the simulation was conducted using two


type of circuit model; the first model included the
transformers capacitances only, while in the second all
transformers constants were taken into account. Figure 2
and figure 3 illustrate the transformer models 1 or 2 used
in the simulation by EMTDC/PSCAD and figure 4 is the
output from the system simulation.

(b)

Figure 1: Construction of the transformer


a. Transformer under study
b. Inner construction of the transformer

Arrestor Module:
Two electrical connections for the arrester
to go to module, double click on the icon
VHigh
High

Steep
Front
N Arrester

Neut
VNeut

These transformer constants were calculated using the


principle of transformer modeling as proposed in the
reference [6].

Ae
2 6 7 .2 [o h m ]

C g 1 ,1 [n F ]

C g 1 ,2 [n F ]

C g 1 ,3 [n F ]

C g 1 ,4 [n F ]

C g 1 ,5 [n F ]

C g 1 ,6 [n F ]

C g 1 ,7 [n F ]

C g 1 ,8 [n F ]

C g 1 ,9 [n F ]

Etank1

the LV connection and the Neutral are


combined since the LV service wire
and the neutral wire are twisted
together and are strongly coupled

TABLE I

L1- 9

III.

Capacitance (pF)
Cs

8,5.10-3
8,7.10-3
9,0.10-3
9,2.10-3
9,5.10-3
9,7.10-3
1,0.10-2
1,02.10-2

Cg
3,4.10-4
3,5.10-4
3,7.10-4
3,8.10-4
4,0.10-4
4,2.10-4
4,4.10-4
4,6.10-4
4,8.10-4

Inductan
ce
(mH)
1.1
1.1
1.1
1.1
1.1
1.1
1.1
1.1

Bx

B
-

TIME

+
Ae

Bx

Surge

Simple Lightning Surge:


I = 50 * [ EXP(-50e-6 * t) - EXP(-1.2e-6 * t) ]

(a)

TRANSFORMER CONSTANTS
Conductor
Layers
(from-to)
L1- 1
L1- 2
L1- 3
L1- 4
L1- 5
L1- 6
L1- 7
L1- 8

2 6 9 .1 [o h m ]

6 8 0 .4 [o h m ]

C s 1 ,1 [n F ]

C s 1 ,2 [n F ]

C s 1 ,3 [n F ]

C s 1 ,4 [n F ]

C s 1 ,5 [n F ]

C s 1 ,6 [n F ]

C s 1 ,7 [n F ]

Etank

C s 1 ,8 [n F ]

C s 1 ,9 [n F ]

The transformers constants have been calculated from


the physical dimension. Each winding layer corresponds
with one transformer winding segment of the transformer
model, where uniformly distributed, so each layer can be
represented by a lumped circuit component [4,5,6]. The
inner-turn and winding-to-ground capacitance, inductance
and resistance of transformer obtained from calculation
based on physical dimension are shown in the table 1.

Resistance
(Ohm)
4
4
4
4
4
4
4
4

CIRCUIT MODEL FOR EMTDC/PSCAD

(b)

In the numerical simulation, we applied the subroutine


JMARTI SETUP where the transformer was considered
as distributed and frequency-dependent parameters. The
surge current used in the simulation was the type 15 of
EMTDC/PSCAD.

Figure 2: Model 1, consisting of transformers capacitances only.


(a) Simulation Model 1 and (b)Schematics Model 1

380

0.0

Arrestor Module:
Two electrical connections for the arrester
to go to module, double click on the icon

Vltg: Tank1 to Line

-1.0k

VHigh
Cs1,5 [nF]

Cs1,4 [nF]

Cs1,3 [nF]

Cs1,2 [nF]

L1,6 [H]

L1,5 [H]

L1,4 [H]

L1,3 [H]

L1,2 [H]

Cs1,1 [nF]
L1,1 [H]

-2.0k

L1,1 [H]
R1,1 [ohm] L1,1[H]

Etank

R..,n [ohm]

R1,6 [ohm]

R1,5 [ohm]

R1,4 [ohm]

R1,3 [ohm]

R1,2 [ohm]

R1,1 [ohm]

R1,1 [ohm]

Steep
Front
N Arrester

2 6 9 .1 [o h m ]

Cs1,6 [nF]

L.., n [H]

6 8 0 .4 [o h m ]

Cs. .,n [nF]

Voltage (kV)

High

Neut

-3.0k

-4.0k

VNeut
Ae Bx
2 6 7 .2 [o h m ]

C g 1 ,1 [n F ]

C g 1 ,2 [n F ]

C g 1 ,3 [n F ]

C g 1 ,4 [n F ]

C g 1 ,5 [n F ]

C g 1 ,6 [n F ]

C g .., n [n F ]

the LVconnection and the Neutral are


combined since the LVservice wire
and the neutral wire are twisted
together and are strongly coupled

-5.0k

TIME

+
Ae Bx

Surge
-6.0k

Simple Lightning Surge:


Etank1

-7.0k

I = 50 * [ EXP(-50e-6 * t) - EXP(-1.2e-6 * t) ]

0.000

0.001m 0.002m 0.003m 0.004m 0.005m 0.006m 0.007m 0.008m 0.009m 0.010m

(a)

(b)
Figure 4: Output Simulation EMTDC/PSCAD Model 1, consisting of
transformers capacitances only.
(a) Output Cs Simulation Model 1 and
(b) Output Cg Simulation Model 1

(b)
IV.

Figure 3: Model 2, including all parameters R, L, C

EXPERIMENTAL SETUP

The experimental measurement was carried out with


help of the setup shown in figure 5. The surge current
was generated from the generator Shaffner NSG 587,
which can produce impulse current up to 5 kA peak
value, with different waveform. In this simulation, we
applied current waveform of 8/20s.

(a) Simulation Model 2 and (b) Schematics Model 2

Where:
Cg1; Cg2...: Capacitance between conductor layers
to ground
Cs1, Cs2...: Capacitance between layers
L11, L12, : Inductance of each conductor layer
R11, R12, : Resistance of each conductor layer

LV Winding

Computer

GPIB
Card

Generator
Impuls

1
2
3
4
5
6

Oscilloscope

Probe HV

7
8
9

HV Winding

Figure5: Experimental Setup

(a)

381

The experiments were done by injecting current to the


upper terminal of the transformer, and the lower terminal
was grounded through a metal strip of copper. The
voltages measured at the winding taps were recorded by
PC, which was connected to the oscilloscope through
GPIB card. During the test, the LV-winding was shortcircuited and connected to ground, as well.

V.

RESULTS

Results of experimental measurement shown in the


figure 6 indicate the behavior of the transformer, when
surge current is injected to the HV terminal (at layer 1).
The voltage measured at the conductor layers (L2, L5,
L8) oscillate during the first 4 s, where the highest
voltage was at the layer L1 (the layer near to the
transformers HV terminal, the curve not presented in the
figure 5). After 6s the voltage at the layer tends to be
stabilised.

Figure6: Voltages at the layers L2, L5, L8 during 6 s after current


injection to the winding.

Figure 7: Comparison between experimental measurement and


numerical simulation.

The figure 7 shows comparison between experimental


measurement and numerical simulation. As can be seen,
the discrepancy between the curves is important at t 1

382

s. The curve obtained from the experimental


measurement decreased toward zero.
After t 1 s, the simulation and experimental
measurement curves were in good concordance that
means the response of the transformer can be adequately
simulated. The numerical simulation shown in the figure
5 was obtained by using the model 2. Simulation by using
the model 1 was not presented in this report, it gave
however similar characteristic.
Results of the study show that the numerical
simulation could not follow the real phenomena at the
beginning of transient period (t < 1s), where at this time
period the capacitances have important influence during
transient.
As the capacitance at each layer applied in the EMTP
simulation has been obtained through approximate
calculation and also assumed as lumped component, the
response at higher frequencies would not follow exactly
the phenomena. This might be caused by stray
capacitances of the layers to ground, which play an
important role for higher frequencies. Experimental
curves of the layers L1 up to L9, show a sharp decrease at
t 1 s, which might eventually be due to reflection at
ground terminal. Yet, this assumption is not very clear,
and needs to be further investigated. The simulation curve
shows in fact same decrease at that instance, but not very
important.
From the comparison between experimental and
numerical simulation it can be concluded that numerical
calculation
can
simulate
adequately
transient
characteristic of a transformer after some microseconds
of surge current flowing through the winding. To
simulate more or less precisely at t 0-1 s proves to be
difficult, since the estimation of the transformer
capacitances could not be done precisely for higher
frequency component of the surge currents.
However, from this study, it has been observed that the
experimental measurements have given results, where the
simulation curve is always higher than the measured
ones. It means the simulation tends to result in more
pessimistic values in overvoltage estimation.
Plot of voltage distribution at t 6 s, both of
experimental and numerical simulation (model 1 and 2) is
presented in the figure 8. As already known from theory,
the voltage gradient at the upper turns (nearest to the HV
terminal) will be the most important. Nearly 60% of the
voltage is born by layer 1 and layer 2.
Again, the
simulation curve indicates more severe situation
compared to the real condition.

[3]

Figure 8: Comparison of voltage distribution obtained from


experimental measurement and simulation, along the transformer
winding at t 6 s.

VI.

CONCLUSION

Transformer model used for transient condition is


normally represented by its capacitance to ground, the
value of which is so for given based on approximated
value depending on the rated voltage and power capacity
of transformer [6,8,9]. This very simple model can
sufficiently represent the transformer when rough
prediction on transient overvoltage due to lightning surge
current is needed [7,11].
The study here presented shows that the more complete
and detail model is sufficiently good for transient
simulation. Somehow, the model could still not be quite
precise when it concerns the first microseconds of the
transient period. On the other hand, after some
microseconds the concordance between the simulation
and measurement curves is good.
The information obtained from this study which can
justify the estimation of transient overvoltage by
numerical simulation, is the fact that it resulted in more
pessimistic overvoltage values compared to the real
condition [10].
From the study, it has been discovered that transformer
capacitances, especially those of conduction layer to
ground, had the most important role in numerical
simulation. Our study shows that there is only small
discrepancy of voltage values obtained by using the
model 1 and model 2.

REFERENCES

[1]
[2]

A. Marti, J.R and S. Chimklai, Simplified ThreePhaseTransformer Model for Electromagnetic Transient Studies,
IEEE Trans. On Power Delivery, Vol. 10, N0. 3, July, 1993.
B.V.Sumangala and G.R.Nagabhushana, Analysis Of Surge
Voltage Distribution In A Model Transformer For Different Types
Of Surges With Turn Resolution And Its Validation. TENCON
2006. 2006 IEEE Region 10 Conference Digital Object Identifier:
10.1109/TENCON.2006.344190 Publication Year: 2006 , Page(s):
1 4.

383

Mobammad S. Naderi, and M. Vakilian A Method for Modeling


Power Transformer Tap Winding for Transient Studies,
Conference Record of the 2004 IEEE International Symposium on
Electrical Insulation, Indianapolis, IN USA, 19-22 September
2004
[4] Morcheds, A.S., et al., A High Frequency Transformer Model for
the EMTP, IEEE Trans. On Power Delivery, Vol. 8, No. 3, July,
1993.
[5] Rodrigo H., and H.Q.S. Dang, Behaviour of Transformer
Winding Under Surge Voltages, Departement of Electrical
Engineering, Royal Melbourne Institute of Tecnology, Australia,
1999.
[6] Xiaoyan Zhu, Huaying Dong, Guishu Liang and Changhu Ji, A
new hybrid model of transformer windings under very fast
transient overvoltages, Electrical Machines and Systems, 2008.
ICEMS 2008. International Conference on Publication Year: 2008
, Page(s): 4296 - 4301
[7] Francis De Leon, Adam Semlyen Complete Transformer Model
for Electromagnetic Transients, IEEE Transactions on Power
Delivery, Vol.9, No. 1, January 1994, pp 231-239.
[8] Pritinda, Chowdhuri, Electrical Transient in Power System,
Research Studies Press Ltd, John Willey & Sons Inc, England,
1996.
[9] Greenwood, A, Electrical Transient in Power Systems, second
edition, John Willey & Sons, Canada, 1991.
[10] Electromagnetic Transient Program (EMTP) Application
Guide, Westinghose Electric Corporaton Power Planning and
Operation Program Electritical System Division, California, 1986.
[11] Jouni Pylvninen, Kirsi Nousiainen, and Pekka Verho, Studies to
Utilise Calculated Condition Information and Measurements for
Transformer Condition Assessment, International Review of
Electrical Engineering, August 2009. Vol. 4. n. 4, pp. 684-689.

You might also like