Memoria Técnica Inel-MT-2003-2023

Estimación de curva de saturación de

transformadores de potencia

1 Propósito
La presente memoria técnica detalla la metodología para estimar la curva de
de potencia
ERNC

saturación de transformadores de transformadores de potencia utilizando el software

ATP/ATPDraw.
de un proyecto
transformadores

2 Alcance
Estimar la curva de saturación de un transformador de potencia de 100 MVA,
13.8/220 kV, YNd11, a ser implementada en un modelo computacional que considera
reactiva

sus características eléctricas y ferromagnéticas.

de saturación de
la compensación

3 Procedimiento

A. Revisión de la información del proyecto

curva

La subestación de potencia de referencia cuenta con dos transformadores de potencia

Estimación de de

para la central hidroeléctrica ubicadas de manera adyacente a la casa de máquinas

Determinación

conectadas cada una a su bahía respectiva en 220 kV, tal como se muestra en la
Figura 1.

Leyenda:
Instalaciones existentes
Instalaciones proyectadas
Ubicación de la malla a tierra
Instalaciones proyectadas

Figura 1 - Diagrama unifilar de la subestación de potencia

La Tabla 1 se muestra los datos básicos del transformador, donde es preciso señalar
que no se cuenta con la información de los ensayos de cortocircuito o de circuito
abierto, las cuales son necesarias para la representación de la curva de saturación
del transformador.

Inel – Escuela Técnica de Ingeniería

Programa de Especialización en el software ATP-EMT → [Link]
Elaborado por Marlon Lujan y Miguel Bejarano
 Tabla 1 – Características eléctricas del transformador de potencia.
Identificación Valor

Tipo Trifásico
Potencia ONAN 100 MVA
Tensión nominal 13.8/220 kV
Frecuencia 60 Hz
Grupo de conexión Ynd11
Impedancia de cortocircuito 13.0 %

B. Estimación de la curva de saturación con dos pendientes

Se construye la curva de saturación preliminar con dos pendientes, correspondientes

a la curva saturada y no saturada, donde el punto de inflexión es el “knee point” el
cual se representa con un valor comúnmente utilizado de 1.2 p.u. [1] y la relación
entre las pendientes se utiliza un valor típico de 500 [2].
Se calcula tanto la inductancia lineal y no lineal de acuerdo al siguiente desarrollo:

𝑆𝑛 𝑃𝑘 𝑉𝑠 2
𝑰𝒃𝒔 = = 4183.7 𝐴 ; 𝑹𝒌 = = 0.0229Ω ; 𝒁𝒃𝒔 = = 1.9Ω
√3 × 𝑉𝑠 3 × 𝐼𝑏𝑠 2 𝑆𝑛
𝑅𝑘

de potencia
𝑹% = × 100% = 1.2% ; 𝑿% = √𝑍𝑐𝑐 2 − 𝑅%2 = 12.9445%
𝑍𝑏𝑠
𝑋%

un proyecto ERNC
𝑿𝒑𝒖 = = 0.129445𝑝𝑢
100%

reactivadedetransformadores
𝑋𝑝𝑢
𝑳𝒔𝒂𝒕 = × 𝑍𝑏𝑝 = 0.0006539𝐻
2×𝜋×𝑓
𝑳𝒏𝒐𝒔𝒂𝒕 = 500 × 𝐿 = 0.32695067𝐻
Donde:

curva de saturación
𝑆𝑛 : Potencia nominal del transformador.
𝑉𝑝 : Tensión del lado secundario.
𝐼𝑏𝑠 : Corriente base en el lado secundario compensación
𝑅𝑘 : Resistencia total referido al secundario
𝑃𝑘 : Pérdidas en el cobre
de lade

𝑍𝑏𝑝 : Impedancia base en el secundario

Estimación

𝑅%: Resistencia total referido al secundario en porcentaje.

Determinación

𝑋𝑝𝑢 : Reactancia referida al secundario en p.u.

𝐿𝑠𝑎𝑡 : Inductancia Saturada
𝐿𝑛𝑜𝑠𝑎𝑡 : Inductancia no Saturada

De acuerdo a las inductancias calculadas se tiene cuenta los puntos ya conocidos en

la curva de la corriente pico vs flujo magnetizante ambos en valores por unidad:

𝐼𝑝𝑝1 = 0; 𝛷𝑝1 = 0; 𝐼𝑝1 = 0; 𝛷1 = 0; 𝛷𝑝2 = 1.2𝑝𝑢

Inel – Escuela Técnica de Ingeniería
Programa de Especialización en el software ATP-EMTP→ [Link]
Elaborado por Marlon Lujan y Miguel Bejarano
 Luego los puntos restantes son calculados de acuerdo a la Figura 2:

Φp2 ∗ √2 Vp
Φ2 = × = 35.86605 Wb
2πf √3
Φ2 − Φ1
Ipp2 = = 0.026221pu
Lnosat
Ipp3 = Ipp2 ∗ 10 = 0.26221pu
(Φp2 − Φp1 ) ∗ (Ipp3 − Ipp2 )
de potencia

Φp3 = + Φp2 = 1.2216pu

ERNC

(Ipp2 − Ipp1 ) ∗ 500

de saturación de
la compensación
curva de un proyecto
transformadores
reactiva

Figura 2 – Curva de saturación del tipo dos pendientes

Estimación de de
Determinación

C. Estimación de la curva de saturación con ajuste polinómico

Se aplica un ajuste de la curva, asintóticamente a la curva del tipo doble pendiente,

para representar de manera óptima el flujo magnetizante de acuerdo a la ecuación
de la corriente que recorre la reactancia magnetizante iterando el valor de la variable
𝛷0 para conseguir la mayor aproximación asintótica.
Para el ajuste de la curva, calculamos el valor de la inductancia saturada y no
saturada de acuerdo a las pendientes de la curva del tipo doble pendiente como se
muestra:
(Φp2 − Φp1 )
lM = m1 = = 45.765708pu
(Ipp2 − Ipp1 )
(Φp3 − Φp2 )
lsat = m2 = = 0.0915314pu
(Ipp3 − Ipp3 )
Además, para el ajuste tomaremos en cuenta esta ecuación [3]:
ΦM ΦM ksat
iMX = (1 + | | ) ( p. u. )
lM Φ0

Donde:
𝑖𝑀𝑋 : Corriente (p.u.) a treves de la reactancia magnetizante.
𝛷𝑀 : Flujo de magnetización (p.u.)
𝑙𝑀 : Reactancia lineal (p.u.)
Inel – Escuela Técnica de Ingeniería
Programa de Especialización en el software ATP-EMT → [Link]
Elaborado por Marlon Lujan y Miguel Bejarano
 𝑙𝑠𝑎𝑡 : Reactancia Saturada (p.u.)
𝛷0 : Parámetro de ajuste de la característica polinomial (p.u.)
𝑘𝑠𝑎𝑡 : Exponente de saturación (grado del polinomio)
Se calcula la corriente que recorre la reactancia magnetizante con un valor de Φ0 =
1 p.u. se obtiene los resultados tabulados en la Tabla 2 y graficados en la Figura 3
para la iteración 1, como el resultado grafico no se aproxima asintóticamente a la
parte saturada, entonces se realiza una iteración 2 escogiendo Φ0 = 1.02 p.u. y se
obtiene la curva polinomial – Iteración 2 mostrado en la Figura 3 con los resultados
tabulado en la Tabla 2.

Figura 3 – Curva de saturación del tipo polinomial con 2 iteraciones.

de potencia
Tabla 2 – Corriente a través de la reactancia magnetizante

Ipico [pu]

un proyecto ERNC
Φpico [pu]

reactivadedetransformadores
Iteración 1 Iteración 2

0.03 0.000656 0.000656

0.1 0.002185 0.002185
0.4 0.008740 0.008740
0.6 0.013116 0.013115
0.8 0.018095 0.017937

curva de saturación
1 0.043701 0.038086
1.05 0.070640 0.058382
1.1 0.124438 0.098636 compensación
1.15 0.229596 0.177051
1.2 0.430200 0.326384
de lade

1.25 0.803596 0.604103

Estimación

1.3 1.482369 1.108722

Determinación

1.31 1.672241 1.249855

1.33 2.123589 1.585326

D. Elaboración del modelo computacional

Para la simulación se realiza una energización en vacío como se muestra en la Figura

4, además se considera la curva de magnetización del núcleo representada por
inductancias no lineales del tipo 98 en el software ATP/ATPDraw [4], clic aquí para
descargar.

Inel – Escuela Técnica de Ingeniería

Programa de Especialización en el software ATP-EMTP→ [Link]
Elaborado por Marlon Lujan y Miguel Bejarano
 V
CURRN IRUSH BCT
I Y

100%

I
Figura 4 – Modelo computacional – Energización en vacío.

E. Ajuste de la corriente inrush de acuerdo a la potencia del

de potencia
ERNC

transformador
de un proyecto

El ajuste final se realiza en el valor pico de la corriente inrush de acuerdo a la curva

transformadores

presentada en la Figura 5, se toma de referencia la potencia del transformador igual

a 100 MVA lo cual corresponde a la relación de la corriente nominal y la corriente
inrush entre los valores 3.25 y 4.75, estos valores nos ayudaran a delimitar la
reactiva

corriente Inrush pico para nuestro ajuste final.

curva
Estimación de de
Determinación de saturación de
la compensación

Figura 5 – Relación Inrush/Inominal de un transformador de potencia.

A partir de la Tabla 2 se implementa la curva de magnetización en las inductancias

no lineales como se muestra en la Figura 6.

Figura 6 – Curva de magnetización del núcleo del transformado de potencia.

Inel – Escuela Técnica de Ingeniería
Programa de Especialización en el software ATP-EMT → [Link]
Elaborado por Marlon Lujan y Miguel Bejarano
 La Figura 7 muestra la corriente inrush, producto de la energización del
transformador, en el lado primario del transformador de potencia, con una corriente
pico de fase a tierra de 1219.7 A, lo cual corresponde a 4.65 veces la corriente
nominal del lado primario.

Figura 7 – Corriente inrush del transformador.

4 Conclusiones
Se estableció y validó el procedimiento para la estimación de curva de saturación de
un transformador de potencia de 100 MVA. Dicho procedimiento es válido para

de potencia
diferentes tamaños de transformador.

un proyecto ERNC
En futuros trabajos se mostrará cómo usar las subrutinas SATURA e HYSDAT para

reactivadedetransformadores
emplearlas cuando se conocen los ensayos de circuito abierto y cómo corregir errores
de convergencia en dichas subrutinas cuando las características no son
monótonamente crecientes.

5 Referencias bibliográficas

curva de saturación
[1] A. N. Washington Luiz, «Transformer Modeling for Transient Studies,» 1994.
compensación
[2] C. P. Jansen, J. G. Slootweg y R. A. de Groot, Calculation of Transformer Inrush
Currents occurring during the Energizing of the Public Grid after a major Black Out,
de lade

2005 IEEE Russia Power Tech, 2005, pp. 1-7.

Estimación

[3] DIgSILENT, «DIgSILENT PowerFactory User Manual 2023,» 2023.

Determinación

[4] Can/Am EMTP User Group, ATP Rule Book-Alternative Transient Program, Portland,
Oregon, 1995.

[5] G. Ziegler, «Siemens Digital Differential Protection».

Inel – Escuela Técnica de Ingeniería

Programa de Especialización en el software ATP-EMTP→ [Link]
Elaborado por Marlon Lujan y Miguel Bejarano