IEEE LATIN AMERICA TRANSACTIONS, VOL. 2, NO.

3, SEPTEMBER 2004

201

Modelado y Simulacin de Transformadores

de Corriente con el Programa EMTP/ATP
Jorge Alberto Aguilar Camarena, Member, IEEE

Resumen -- El presente trabajo describe el modelado y

simulacin de transformadores de corriente (TCs) utilizados
con relevadores de proteccin, con datos de pruebas de
laboratorio, de saturacin y relacin, determinando su
respuesta ante diferentes condiciones del sistema elctrico. El
anlisis del modelo se basa en la simulacin del TC, utilizando
el modelo de transformador saturable del ElectroMagnetic
Transient Program (EMTP/ATP) y la comparacin con
oscilografas de las pruebas de saturacin. Se concluye
enfatizando el estudio de los TCs ante fenmenos transitorios
en escenarios significativos para establecer la adecuada
seleccin, asegurando la operacin confiable de los relevadores
de proteccin.

En el presente trabajo se utiliz TCs, cdigo C, que

cubre a la mayora de los transformadores de corriente tipo
toroidal (bushing), con devanados uniformemente
distribuidos y perdidas despreciables, donde la clase de
precisin puede ser determinada utilizando un burden
estndar para determinar el voltaje de excitacin.
Esta relacin de voltaje se basa en el voltaje que el TC
puede soportar, a travs del burden, para un valor de 20
veces la corriente nominal simtrica, sin exceder 10% de
error en la relacin de transformacin. La tabla I muestra la
relacin de los TCs cdigo C, entre el burden estndar
secundario y el voltaje secundario.

Palabras clave EMTP/ATP, simulacin digital, transformadores de corriente, saturacin.

I. INTRODUCCIN

L EMTP/ATP es un programa para el modelado y

simulacin digital, ampliamente utilizado y reconocido
en el ambiente universitario. Para el ingeniero de
proteccin, la seleccin de la clase de precisin de los
transformadores de corriente (TCs) y el anlisis del mismo
son de suma importancia ante diversos disturbios, ya que
estn directamente relacionados con la precisin y
desempeo de los relevadores de proteccin.
El comportamiento de los TCs, ante el crecimiento de
los sistemas de energa, puede verse afectado por:
a)

Corrientes de falla asimtricas con componente

decreciente de CD.
b) Magnitudes elevadas de corrientes de falla,
combinadas con alto burden en secundario del TC.
La salida de corriente secundaria puede ser afectada de
manera drstica cuando el TC trabaja en la regin no lineal
de la caracterstica de excitacin [1].
De acuerdo con ANSI, la clase de precisin de los TCs
utilizados para proteccin est determinada por el estndar
IEEE C57.13-1993 [2]. La capacidad del TC y el rango de
voltaje esta definida por una letra. El cdigo de designacin
esta indicada por las letras C y T.

El autor est en el programa de Ingeniera Elctrica, dentro de la Facultad

de Ingeniera Mexicali, Universidad Autnoma de Baja California,
Mexicali B.C. Mxico. (email: [email protected])

TABLA I
TC TIPO C BURDEN Y VOLTAJE

Burden
1
2
4
8

Cdigo C
C100
C200
C400
C800

Voltaje
100 V
200 V
400 V
800 V

II. COMPORTAMIENTO DE LOS TCs

El comportamiento del transformador de corriente
puede ser estimado de acuerdo con el concepto rea volttiempo [1]. El voltaje visto en el secundario del TC est en
funcin de la corriente de falla iF y el total de burden ZB,
como;

dI
dt

>V @

(1)

donde N es el nmero de vueltas del secundario y es el

flujo del ncleo en webers (Wb). Integrando (1) para
determinar la densidad de flujo a travs de la seccin
transversal en el ncleo tenemos:

NI

NBA

v dt

(2)

De la expresin anterior, en la figura 1 el rea bajo la

curva de la seal de voltaje representara la densidad de
flujo de saturacin (Bs).

202

IEEE LATIN AMERICA TRANSACTIONS, VOL. 2, NO. 3, SEPTEMBER 2004

El transformador de corriente que se utiliz para la prueba es

cdigo C, tipo OE12, que tiene las siguientes caractersticas;

20
[kA ]
15
10

Relacin
Clase de Precisin
Burden
Resistencia secundario

5
0
-5

1200/5
C200
2
0.0025 /vuelta

-10
-15
0 .0 0 5 1

0 .0 4 3 1

0 .0 2 4 1

(f ile s im u la c io n _ TC _ id e a l1 . p l4 ; x -v a r t ) c : I F B

0 .0 6 2 0
-

c :IF 1B -

0 .0 8 1 0
c :IF B

[s ] 0 .1 0 0 0

--c : I F 1 B -

Fig. 1. Voltaje en burden para corriente de falla asimtrica.

Las mediciones para la obtencin de la curva de saturacin

V-I rms se realizaron de acuerdo con el estndar IEEE
C57.13.1993 utilizando la mxima relacin de transformacin
del TC y arrojando los resultados que se muestran en la tabla II.

Considerando mxima asimetra del voltaje en (2),

podemos establecer que

NI

NB s A

TABLA II
VALORES DE PRUEBA CURVA DE EXCITACIN

Ve (voltios)
20
40
80
140
160
180
200
220
240
280

i ( e  Rt /L  cos Z t ) dt

B F

Z B i F (1  e  Rt / L )  sen Z t
R

(3)

donde:
Bs = Densidad [Wb/m2] de flujo de saturacin
N = Nmero de vueltas del TC
A = rea [m2] del ncleo del TC
iF = Magnitud de la corriente [A] secundaria de falla
ZB = Burden [ohm] en terminales del secundario del TC
L/R = Constante de tiempo [s] del circuito primario fallado
Utilizando el valor mximo de la expresin (3),
tenemos:

NB s A Z

Z B iF  1
R

(4)

Considerando que el rango de voltaje en el TC puede

soportar 20 veces su valor a travs del burden de acuerdo
con [2], podemos establecer que para prevenir la saturacin
la expresin (4) quedar como:

Ie (amperes)
0.01
0.02
0.03
0.04
0.05
0.06
0.08
0.14
1
10

Ze (ohmios)
2000
2000
2666.66667
3500
3200
3000
2500
1571.42857
240
28

La impedancia de la rama de magnetizacin (Ze) de la

tabla II se obtuvo por la divisin directa de cada valor
medido del voltaje de excitacin (Ve) por su valor
correspondiente de corriente de excitacin (Ie).
La curva de saturacin del TC bajo prueba con los datos de
medicin de tabla II queda como muestra la figura 2, donde se
puede apreciar la caracterstica no lineal de la impedancia de
magnetizacin tomando valores desde 2000 para una
corriente de 0.01 amperes hasta un mximo de 3500 para 0.04
amperes. Este ltimo es el punto de mxima permeabilidad
localizada tangencialmente a 45 de la curva de saturacin.
Curva de Saturacin

20

X
iF
 1

(5)

1000

Ve (V)

100

III. CARACTERSTICA DEL TC EN PRUEBA

10

Las caractersticas de excitacin que determinan la clase

de precisin de los TCs, pueden ser obtenidas de curvas V-I
rms de fabricantes y por pruebas de saturacin en campo o
laboratorio, en este ltimo caso del presente trabajo.

0,01

0,1

1
Ie (A)

Fig. 2. Curva caracterstica de saturacin del TC.

10

CAMARENA : MODELING AND SIMULATION OF CURRENT

Como parte de la misma prueba se grabaron registros

oscilogrficos de la corriente de excitacin (Ie) para
observar la forma de onda y posteriormente poder
compararla con el modelo de simulacin, como se
presenta en la figura 3 para un voltaje de excitacin de
240 Volts.

203

LP = Inductancia [:] en devanado primario

VRP = Nivel de voltaje [V] devanado primario
RS = Resistencia [:] en devanado secundario
LS = Inductancia [:] en devanado secundario
VRS = Nivel de voltaje [V] devanado secundario
RMS = Caracterstica de saturacin en valores RMS
Para la prueba de simulacin con el ATP se utilizaron la
caracterstica de saturacin de figura 2, los datos
caractersticos del TC y los datos de tabla II, quedando la
informacin de entrada para la simulacin para saturacin
de acuerdo con la tabla III.

TABLA III
PARMETROS DE PRUEBA TC

Fig. 3. Registro oscilogrfico de Ie en prueba de saturacin.

IV. MODELO DE SIMULACIN

Los transformadores de corriente, como se ha
comentado, estn sujetos a fenmenos de saturacin ante
fallas de corto circuito transitorias, por lo que se utiliz el
programa EMTP/ATP [3] y su interfase grfica ATPDraw
[4] para simular y reproducir los datos de oscilografas de
pruebas.
Aunque existen algunas guas para el modelado de
transformadores de corriente [5-6], se utiliz como
referencia el documento presentado en [7]. El
transformador de corriente es representado utilizando el
elemento de Transformador Saturable de acuerdo con la
figura 4

Devanado 1
Bajo voltaje

Datos
IO
FO
RMAG
RP
LP
VRP
RS
LS
VRS
RMS

Valores
0.1131
0.7503
3500
0.6
0
240
0
1E-7
1
1

As, de acuerdo con la figura 5, el circuito resultante en

el ambiente grfico ATPDraw es:

Devanado 2
Alto voltaje

Fig. 5. Modelo de simulacin para saturacin de TC.

Fig. 4. Modelo de Transformador Saturable del ATP.

donde:
I0 = Corriente [A] en la rama de magnetizacin
F0 = Flujo [Wb] en la rama de magnetizacin
RMAG = Resistencia [:] en la rama de magnetizacin
RP = Resistencia [:] en devanado primario

Para las pruebas de saturacin y relacin del

transformador de corriente, se utiliz la prestacin
$PARAMETER para la definicin de variables, en
conjunto con POCKET CALCULATOR VARIES
PARAMETERS (PCVP), para manejar en un solo
modelo varias simulaciones de acuerdo con las pruebas
de laboratorio.
El circuito equivalente para la simulacin de relacin
de transformacin (RTC) qued representado por la
figura 6:

204

IEEE LATIN AMERICA TRANSACTIONS, VOL. 2, NO. 3, SEPTEMBER 2004

Con los resultados anteriores, se procedi a realizar

simulaciones con datos reales del sistema elctrico de
energa de distribucin al que estaran conectados los TCs,
para observar su comportamiento ante fallas de acuerdo con
el diagrama trifsico siguiente, figura 9:
X/R = 21

Fig.6. Modelo de simulacin para RTC.

Circuito

13.8 KV

V. RESULTADOS DE LA SIMULACIN
Para la prueba de la caracterstica de saturacin, se
obtuvieron los resultados de la figura 7, con un voltaje de
excitacin de 240 V rms:

1200/5

Fig. 9. Sistema de distribucin simulado.

3

Ie = 1.0028 [A-RMS]

Los datos del sistema elctrico utilizado para la

simulacin se detallan en la tabla IV.

Ie [A]

TABLA IV
PARMETROS DEL SISTEMA ELCTRICO SIMULADO

-1

TC

-2

-3
16.67

22.71

28.75

34.79

40.83

46.88

(f ile Curv a_Saturacion_TC_UABC_RMS.pl4; x-v ar t) c:

-SEC1

52.92
t: IERMS

600/5

58.96 [ms] 65.00

Fig. 7. Corriente de excitacin (Ie): forma de onda y valor RMS

De acuerdo con las figuras 3 y 7, las formas de onda y

valores de la corriente de excitacin son muy similares.
Para la corriente de excitacin (Ie ) RMS, la medicin de
laboratorio de acuerdo a la prueba es Ie = 1 amper, y para
la simulacin con el ATP toma el valor Ie = 1.0028
amperes.
Los resultados de las pruebas de relacin con el modelo
de simulacin del ATP se muestran en la figura 8, grfica
de factor de correccin de relacin.

Burden
[:]
1.62+0.784

PCC
[MVA]
193

X/R

Tensin
[kV]

21

13.8

En funcin de lo anterior se simul falla trifsica que

sera el caso menos probable a la salida del sistema de
energa con un instante de ocurrencia 0 para inicio del
corto circuito obtenindose una respuesta con menor
componente de DC para el transitorio. El registro de la
simulacin de los comportamientos obtenidos de las
corrientes se muestra en la figura 10.
10
[kA]
5

Factor Correccin Relacin (RCF)

-5

1.01
-10

RCF

1
0.99

Isecundario TC

-15

0.98

Iprimario TC

0.97
0

50

100

150

Porciento corriente primaria

Fig. 8. Factor de correccin de relacin.

200

-20
0.08

0.10

0.12

0.16

0.14

(f ile simulacion_TC_ideal1.pl4; x-v ar t) c:S1

factors:
1
240
offsets:
0
0

c:P2C 1
0

0.18

0.20

0.22 [s] 0.24

Fig. 10. Corrientes primaria (Ip) y secundaria (Is) TC 1200/5 C200.

CAMARENA : MODELING AND SIMULATION OF CURRENT

En la figura 10 se observa que la condicin de

saturacin crtica permanece por cerca de 100 ms (6 ciclos),
tiempo en el cual se pueden dar sobre alcances por ejemplo
de relevadores de sobrecorriente instantneos y
temporizados teniendo una respuesta con un retardo de
tiempo mayor que el ajustado por la disminucin de
corriente. Lo mismo para relevadores conectados
diferencialmente, ya que en la comparacin de las
corrientes secundarias tendran distintos valores.
Con el modelo anterior se pueden simular fallas de
mayor impacto factibles de presentarse en los sistemas
elctricos de distribucin. Con las magnitudes obtenidas se
pueden grabar archivos de acuerdo con el estndar IEEE
C37.111-1991 Common Format for Transient Data
Exchange (COMTRADE) para verificar los ajustes de los
relevadores de proteccin que intervengan en los casos de
inters y estimar el burden que podran soportar los TCs
ante eventos transitorios de acuerdo con el criterio de la
expresin (6).
VI. CONCLUSIONES
Se ha analizado la respuesta de los transformadores de
corriente tipo bushing probados en el laboratorio de
Ingeniera Elctrica de la UABC, con el apoyo del
programa Alternative Transient Program (ATP).
Los resultados reproducen las caractersticas
presentadas en las oscilografas de registros de las pruebas
de saturacin y permiten simular los modelos de TCs con
escenarios significativos para comprobar los ajustes
ptimos de las protecciones asociadas a los sistemas de
energa.
VII. BIBLIOGRAFA
[1]
[2]
[3]
[4]
[5]

[6]

[7]

Stanley E. Zocholl, Jeff Roberts, Gabriel Benmouyal, Selecting CTs

to Optimize Relay Performance, WPRC, Oct. 1996.
IEEE Std .C57.13-1993, IEEE standard Requirements for
Instruments Transformers, IEEE, 1993.
ATP/EMTP CanAm Users Group, ATP Rule Book; 1995.
H. K. Hoidalen, ATPDraw Users Manual, Sintef, Versin 3.5,
Windows, Octubre, 2002.
Working Group C-5 of IEEE PSRC, Mathematical Models for
Current, Voltage, and Coupling Capacitor Voltage Transformers,
IEEE Trans. Power Delivery, vol.15, n1, Jan. 2000.
S. E. Zocholl, D. W. Smaha, SCS. Current Transformer Concepts.
Presentado a la 19th Western Protective Relay Conference, Spokane,
Washingtong, Oct. 20-22, 1992.
R. Folkers, SEL, Determine Current Transformer Suitability Using
EMTP Models, Presentado a la 26th Western Protective Relay
Conference, Spokane, Oct 26-28, 1999.

VIII. BIOGRAFA
Jorge Alberto Aguilar Camarena. Obtuvo la
Licenciatura como Ingeniero Mecnico Electricista,
con especialidad en Electricidad, y el grado de Maestro
en Ciencias en Ingeniera de Sistemas por la
Universidad Autnoma de Baja California (UABC).
Actualmente cursa estudios dentro del Programa

205

Doctoral en Ingeniera. Desde 1992 se desempea como Maestro en la

carrera de Ingeniera Elctrica en la UABC, en las reas de Protecciones y
Sistemas de Distribucin de Energa. Trabaj en la Comisin Federal de
Electricidad, de 1987 a 2003, en los procesos de Transmisin y
Distribucin, especialidad de Protecciones. Es miembro IEEE y auditor
lder en ISO 9001:2000 dentro del sistema de calidad de la Facultad de
Ingeniera. Con inters profesional en el anlisis de sistemas elctricos de
Proteccin, Potencia y Calidad de Energa.