JORNADAS TCNICAS ISA - 2012
METODOLOGA PARA EL ANLISIS DE PERTURBACIONES EN SISTEMAS DE
POTENCIA MEDIANTE ATPDRAW CASO DE PUENTE ABIERTO EN LA LNEA
ESMERALDA YUMBO 2 A 230kV
Diego Alejandro Tejada Arango
Sadul Segundo Urbaez Gonzalez
Especialista Planeacin Operacin
Direccin Planeacin de la Operacin
XM-Colombia
dtejada@[Link]
Especialista Operacin Red
Direccin Gestin de la Operacin
ISA-Colombia
surbaez@[Link]
Categora: Sistemas de control, proteccin y telecomunicaciones
RESUMEN
Este documento presenta el anlisis, la metodologa, simulacin y experiencia de anlisis posoperativo de
Interconexin Elctrica S.A. (ISA) asociada a una falla en un circuito compensado serie por un puente
abierto que ocasion indisponibilidad del circuito fallado por casi 24 horas. Este artculo hace especial
nfasis en la localizacin de la falla y metodologa de simulacin para la caracterizacin desde el punto de
vista elctrico del evento.
PALABRAS CLAVES.
Interconexin Elctrica S.A. (ISA), Sistema de transmisin Nacional (STN), ATPDraw, fallas longitudinales,
recierre, registrador de fallas, compensacin serie, Sistema de Informacin de Descargas (SID), anlisis de
fallas.
INTRODUCCIN
Las redes de transmisin de un sistema de
potencia, estn constantemente expuestos fallas
ocasionadas
por
descargas
atmosfricas,
vegetacin, acercamientos con circuitos de otros
niveles de tensin, atentados mal intencionados,
entre otros. Por tanto, las fallas ms comunes en
las lneas de transmisin son los cortos circuitos a
tierra y entre fases, como lo son fallas
monofsicas, trifsicas, bifsicas a tierra y
bifsicas aisladas. Sin embargo, las redes de
transmisin son cada vez ms viejas y empiezan a
ser ms comunes fallas del tipo longitudinal, como
conductores rotos (tradicionalmente ms comunes
en niveles de distribucin). En la Tabla 1 se
presenta una estadstica de las fallas en circuitos
de transmisin presentadas en el perodo 2001 2010 en la red de ISA.
Causa
Descarga Atmosfrica
AMI (Actos Mal Intencionados)
Quema o Fuego Bajo la Lnea
Dao en Conductor - Bajante
Acercamiento Conductores
Fallas
1,781
445
59
22
25
Causa
Alta Vegetacin
Cada de rbol
No Establecida
Otros
Total
Fallas
81
10
1,257
96
3,776
Porcentaje
2.1%
0.3%
33.3%
2.5%
-
Fuente: Elaboracin Propia
Tabla 1. Parmetros para los equivalentes.
Como se aprecia en la Tabla 1, las fallas
asociadas a ruptura de conductores y daos en
bajantes han representado un promedio de 2 fallas
al ao en la red de transmisin de ISA. Por esta
razn el anlisis, caracterizacin y simulacin de
este tipo de fallas es de especial inters y para ello
en este artculo se describe y analiza el evento del
26/12/2010 en el circuito Esmeralda Yumbo 2 a
230kV, el cual tuvo como causa de la falla, un
puente abierto en la fase C de la torre 321.
Porcentaje
47.2%
11.8%
1.6%
0.6%
0.7%
1
JORNADAS TCNICAS ISA - 2012
1. DESCRIPCIN DEL EVENTO
El 26/12/2010 a las 18:33 se present el disparo
del circuito Esmeralda - Yumbo 2 a 230 kV con 130
MW ante falla en el circuito y operacin de sus
protecciones principales, ver Figura 1. La falla se
localiz en la torre 321 a 16.51 km desde Yumbo
por un puente abierto en la fase C del circuito, en
la Figura 2 se presenta el registro fotogrfico de la
falla mediante la inspeccin area realizada por
ISA. La lnea tiene una compensacin serie en la
subestacin San Marcos 230kV la cual qued por
by-pass a raz de la falla en el circuito.
Fuente: Inspeccin de la lnea realizada por ISA
Figura 2. Fotografa de la causa de la falla
2. ANALISIS DEL EVENTO ESMERALDA
YUMBO 2
Como se mencion en numerales anteriores la
falla se produjo por un puente abierto en la fase C
de la torre 321. En la torre fallada la fase C est
en la parte superior y la fase B en la intermedia
(ver Figura 2), al abrirse el puente de la fase C, se
hizo contacto con la fase B y produjo la falla entre
las fases del circuito. Se produjo el disparo de
ambos extremos por recierre en falla en Esmeralda
y falla en tiempo de reclamo en Yumbo. La
compensacin serie se fue a by-pass (cierre del
interruptor 2CS10), luego se realiz intento de
cierre desde Yumbo el cual fue exitoso, luego se
realiz intento de cierre manual desde Esmeralda
el cual no fue exitoso por falla en las fases B y C,
cuando se realiz este cierre desde Esmeralda, en
Yumbo se detect la falla slo en la fase B y la
baha realiz recierre monofsico exitoso en
Yumbo. Ambos cierres manuales se realizaron
con la compensacin serie en by-pass. Debido a
este rechazo en la energizacin, se procedi a
revisar la localizacin de la falla e inspeccin del
circuito.
Fuente: ISA-Sistema de Transmisin Nacional
Figura 1. Diagrama unifilar
2.1 RELES INSTALADOS PARA PROTECCION DEL
CIRCUITO
Para la proteccin del circuito se tienen
instalados en ambos extremos de la lnea, rels
PXLN y LFDC ambos del fabricante Alsthom. El
2
JORNADAS TCNICAS ISA - 2012
rel PXLN es un rel con principio de operacin
distancia y como particularidad para la proteccin
de esta lnea su zona 1 est temporizada 120ms,
con el fin de dar prioridad al rel LFDC, el es de
un principio de operacin por cantidades
superimpuestas (V, I).
2.2 ANLISIS DE LOS REGISTROS DE FALLA EN
ESMERALDA
En la baha BL Esmeralda 230 kV - Yumbo 2 se
detect Falla en la fase C en el rel LFDC (no se
detect falla en la fase B, las tensin de fase no
disminuy y la corriente no aument lo suficiente
como para que los deltas generaran disparo de la
fase B) el rel PXLN present omisin de
actuacin en la primera falla lo cual es normal
porque la zona 1 est temporizada 120ms. La falla
fue aclarada por el interruptor L250, con apertura
de la fase C a los 77 ms y con recierre a los 802
ms (tiempo muerto: 725 ms, una vez realizado el
recierre se present nuevamente falla en la fase C
a los 805 ms, aclarada por el interruptor L250, con
apertura de las fases ABC a los 872 ms por orden
del rel LFDC y del PXLN por su funcin de
recierre en falla, ver Figura 3.
Fuente: Registrador de Fallas en Esmeralda
Figura 4. Anlisis de la corriente de prefalla
Fuente: Registrador de Fallas en Esmeralda
Fuente: Registrador de Fallas en Esmeralda
Figura 3. Oscilografa del evento en Esmeralda
Figura 5. Anlisis de la tensin durante el ciclo de
recierre monofsico
De este registro es importante destacar dos
aspectos importantes para detectar la apertura de
la fase C en algn punto de la lnea. Primero, la
corriente de prefalla de la fase C es capacitiva, lo
que indica que ya no haba flujo de potencia por
esta fase en el circuito, ver Figura 4. Segundo, la
tensin de la fase C durante el ciclo de recierre
monofsico es igual al de la fase B, lo cual indica
la unin de ambas fases, ver Figura 5.
2.3 ANLISIS DE LOS REGISTROS DE FALLA EN
YUMBO
En la baha BL Yumbo 230 kV - Esmeralda 2 se
detect falla en las fases B y C en el rel LFDC, el
rel PXLN present omisin de actuacin en la
primera falla lo cual es normal porque la zona 1
est temporizada 120ms. La falla fue aclarada por
el interruptor L190, con apertura de las fases ABC
a los 70 ms y con recierre trifsico Barra Viva
Lnea Muerta a los 663 ms (tiempo muerto: 593
ms). Una vez realizado el recierre se present
nuevamente falla en la fase B a los 805 ms
(despus del recierre del otro extremo) el rel
LFDC dispar nuevamente, la segunda falla fue
aclarada por el interruptor L190, con apertura de
las fases ABC a los 867 ms, ver Figura 6.
JORNADAS TCNICAS ISA - 2012
100.4km con el rel PXLN, en Yumbo el rel PXLN
no gener registro para el evento y la localizacin
con el RDF INDACTIC de Yumbo que marc una
localizacin entre 29.9km y 42km desde esta
subestacin utilizando el loop fase B a
tierra. Como se puede apreciar con estas
localizaciones se tiene una localizacin errada de
la falla porque la suma de ellas corresponde a la
longitud real de la lnea. Esta situacin ocasion
que la inspeccin de la lnea se realizara en un
sector diferente al real de la falla, lo cual gener
extra costos en la inspeccin y mayor
indisponibilidad del circuito.
Fuente: Registrador de Fallas en Yumbo
Figura 6. Oscilografa del evento en Yumbo
Al igual que para el extremo de Esmeralda, para el
extremo de Yumbo tambin vale la pena destacar
de los registros de la falla la corriente de prefalla
en la fase C, la cual tambin indicada apertura del
conductor antes de la falla, ver Figura 4.
Fuente: Registrador de Fallas en Yumbo
Figura 7. Anlisis de la corriente de prefalla
2.4 ANLISIS DEL TIPO DE FALLA
Se analiz la falla desde ambos extremos,
mostrando que por la caracterstica de la falla, sta
fue bifsica en un extremo y monofsica en el
otro. Ahora bien, una pregunta que surge en el
momento del anlisis es si la falla fue bifsica a
tierra o bifsica aislada. Con el fin de entender un
poco mejor el tipo de falla se procedi a realizar
una sincronizacin de los registros de ambos
extremos y mediante el anlisis fasorial, en el
momento de la falla, poder determinar el tipo de
falla presentada. En la Figura 8, se muestra como
los fasores de corriente de las fases B en ambos
extremos de la lnea suman el valor negativo del
fasor de corriente de la fase C en el extremo donde
la falla se present de manera bifsica. Esto nos
lleva a lanzar la hiptesis de que la falla fue
bifsica aislada de tierra, esta hiptesis ser
validad mediante la simulacin del evento en el
ATPDraw en el numeral 4 del presente artculo.
Fuente: Elaboracin Propia
Figura 8. Anlisis fasorial de las corrientes de falla
A pesar de ser una falla sin resistencia de falla, por
sus caractersticas elctricas, produce un alto error
en el clculo de la localizacin si empleamos
mtodos de localizacin basados en la impedancia
de la lnea.
2.5 ANLISIS DE LA LOCALIZACIN DE LA
FALLA
La longitud total del circuito es de 193.3 km, la
localizacin de la falla desde Esmeralda fue de
4
JORNADAS TCNICAS ISA - 2012
Por esta razn se procedi se procedi a validar la
localizacin de la falla con el registro de la
corriente capacitiva desde el extremo de
Yumbo. Este mtodo de localizacin est basado
en que para una lnea de transmisin en estado
estable a frecuencia industrial la magnitud de la
corriente capacitiva de la lnea es directamente
proporcional a la longitud de la lnea. Para el caso
de una fase abierta la corriente capacitiva
corresponder a los kilmetros desde la
subestacin hasta el punto donde se encuentra
abierta la fase.
Para nuestro caso particular del circuito Esmeralda
Yumbo 2, al energizar desde Yumbo el registro
mostr que las fases sanas (A y B) tena una
corriente de 88.9Amp capacitivos y que la fase
afectada (C) tena una corriente de 8.74Amp
capacitivos, ver Figura 9. Haciendo una regla de
tres (si para 193.3km son 88.9km, para 8.74Amp
sern 20km) se determin que la localizacin de la
falla era alrededor de 20km desde Yumbo. El CTE
procedi a realizar la inspeccin en la nueva
localizacin y se determin que la falla estaba en
la torre 321 a 16.51km desde Yumbo (error de
3.49km en la localizacin mediante este mtodo).
simulacin en ATPDraw. Para el modelado de
cualquier sistema elctrico existen procedimientos
dentro del proceso de generacin de un caso de
simulacin y que en todos los casos en que se
quiera reproducir un evento deben ser aplicados
de una u otra forma. Por esta razn, se proponen
los siguientes puntos importantes para todo
proceso de modelado de fallas en sistemas de
potencias:
1. Se debe visualizar el punto donde ocurri la
falla y a partir de eso determinar que partes de
todo el sistema se modelaran de forma
detallada y cules de manera sencilla.
2. Obtener los datos necesarios para la
representacin en ATPDraw de los dispositivos
que sern modelados. En caso de no
conseguir estos datos, utilizar datos tpicos o
datos de un dispositivo similar.
3. Por medio de un Software de simulacin de
sistemas de potencia (para el caso de ISA
Power Factory DIgSILENT) obtener los datos
de los equivalentes de Thevenin (R0, X0, R1 y
X1) de las zonas del sistema que no se
modelaran en detalle.
4. Al tener todos los datos necesarios, se
representa en el software ATPDraw cada uno
de los dispositivos (Lneas, Transformadores,
Reactores, Condensadores, etc.) que se
encuentran dentro de la zona donde ocurri la
falla y son relevantes en el anlisis del evento.
5. Al tener modelado el sistema (sin la falla), se
debe sintonizar el modelo con los valores de
prefalla (Tensiones y Corrientes) entregados
por los registros de los Registradores de falla y
rels.
6. Al tener sintonizado el sistema seleccionamos
el punto real donde ocurri la falla en el
modelo, en caso de ser en una lnea de
transmisin, utilizar la distancia entregada por
el localizador de fallas como una primera
aproximacin.
7. Por medio de los registros se debe visualizar el
tipo de falla que ocurri en el evento. En caso
de no estar muy claro en el registro, se debe
realizar
una
sensibilizacin,
aplicando
diferentes fallas al modelo hasta encontrar la
que represente las formas de ondas del
Fuente: Elaboracin Propia
Figura 9. Magnitudes de la corriente capacitiva
3. MODELO DE SIMULACION EN ATPDRAW
Con el anlisis del evento presentado en los
numerales anteriores, se procedi a corroborar el
comportamiento de la falla mediante una
5
JORNADAS TCNICAS ISA - 2012
registro durante la falla (la falla se ubicar de
acuerdo al tem anterior).
8. Al tener el tipo de falla, por medio de los
interruptores
del
sistema
modelado
representaremos la secuencia de apertura y
cierre de cada uno de los interruptores que
intervienen en el evento.
9. Luego se simula en el programa y se observa
que los aportes de corrientes sean
consecuentes con los aportes vistos en el
registro, en caso contrario, mover el punto de
falla.
10. Dentro del proceso de sintonizacin de la falla,
existen ciertos parmetros que pueden ser
variados, entre los cuales estn la ubicacin
de la falla, impedancia de la falla y los valores
de secuencia cero.
evento. Los parmetros obtenidos se presentan
en la Tabla 2.
Parmetro
Icc 3f [kA]
Icc 1f [kA]
R0 [ohm]
X0 [ohm]
R1 [ohm]
X1 [ohm]
Vprefalla [kV]
S/E Esmeralda
13.48
13.15
1.306
9.863
0.817
9.279
131
S/E Yumbo
10.39
11.48
0.979
8.048
0.902
12.044
130
Fuente: Elaboracin Propia
Tabla 2. Parmetros para los equivalentes.
3.2 LNEA DE TRANSMISIN
La lnea de transmisin Esmeralda Yumbo 2 a
230kV est en una estructura de circuito doble con
un diseo de torre Normalizada Clase II, la
configuracin es vertical, con un conductor de fase
AAAC 1200 61/0 y un cable de guarda:
ALUMOWELD AW 7 NUM9 7/0. En la Tabla 3 se
presentan los parmetros de la lnea y en la Figura
11 la configuracin del modelo creado en
ATPDraw para la simulacin.
De acuerdo a la metodologa expuesta, a
continuacin se presentan los parmetros de
modelado para el caso del evento de la lnea
Esmeralda Yumbo 2.
3.1 EQUIVALENTES DE CORTO CIRCUITO
Para los equivalentes de corto circuito se
emplearon los elementos LINESY_3 (para modelo
de las redes de secuencia) y ACSOURCE (para la
fuente de tensin trifsica), ver Figura 10.
Parmetro
Longitud [km]
Ao de entrada
R0 [ohm/km]
X0 [ohm/km]
B0 [S/km]
R1 [ohm/km]
X1 [ohm/km]
B1 [S/km]
Lnea Esmeralda Yumbo 2
193.3
1971
0.4013
1.3275
2.1592
0.0547
0.4618
3.4880
Fuente: XM Parmetros Tcnicos (PARATEC)
Tabla 3. Parmetros de la lnea de transmisin.
Fuente: Elaboracin Propia
Figura 10. Elementos usados en los equivalentes
Los parmetros de secuencia (R0, X0, R1, X1)
empleados para la simulacin se obtuvieron a
partir de los niveles de corto circuito trifsico y
monofsico de cada subestacin. La tensin de
prefalla lnea-neutro se obtuvo de los registros del
6
JORNADAS TCNICAS ISA - 2012
Fuente: Elaboracin Propia
Figura 13. Modelo realizado en ATPDraw
Fuente: Elaboracin Propia
4. SIMULACIONES DEL EVENTO CON
ATPDRAW
Con los elementos del numeral 3 y la informacin
de los tiempos de operacin reales del evento del
numeral 2, se procedieron a realizar las
simulaciones con el ATPDraw. Para confirmar la
hiptesis de una falla bifsica a aislada de tierra,
se procedi a realizar este tipo de falla en el
modelo y a validar los valores de simulacin
respecto a los valores reales de los registros de
falla. A continuacin en la Figura 14, Figura 15,
Figura 16 y Figura 17, se presenta una
comparacin de los resultados y algunos valores
de especial inters para realizar una comparacin,
la primera grfica siempre ser asociada al registro
real y la segunda al resultado de la simulacin.
Figura 11. Modelo de lnea usado en ATPDraw
3.3 COMPENSACIN SERIE
Para la simulacin de la compensacin serie de la
lnea, se emple un modelo simplificado
condensador ms interruptor de by-pass, ver
Figura 12.
Fuente: Elaboracin Propia
Figura 12. Modelo empleado para la compensacin
serie
Este condensador serie, compensa la lnea en un
50% de la lnea, tiene una capacidad de 3x22Mvar,
tensin nominal de 31.5kV y tiene para su
proteccin un spark gap, un varistor y un
interruptor. Para esta simulacin slo se model el
interruptor debido a que producto de la energa
producida en el varistor hizo que la compensacin
rpidamente se fuera a by-pass (cierre del
interruptor) y la compensacin no tuviera ms
efecto en la falla.
3.4 MODELO PARA LA SIMULACION
Con todos los parmetros de entrada descritos en
los numerales anteriores y los datos obtenidos del
anlisis de la falla, se cre el modelo de simulacin
en ATPDraw, presentado en la Figura 13:
Fuente: Elaboracin Propia
Figura 14. Tensiones en el extremo de Esmeralda
7
JORNADAS TCNICAS ISA - 2012
Fuente: Elaboracin Propia
Fuente: Elaboracin Propia
Figura 17. Corrientes en el extremo de Yumbo
Figura 15. Tensiones en el extremo de Yumbo
En la Figura 18, se presenta una comparacin de
las corrientes de falla punto a punto entre la
simulacin y el registro real de falla.
Fuente: Elaboracin Propia
Figura 18. Comparacin de las corrientes de falla
De acuerdo con los resultados anteriores, una
comparacin entre los resultados entregados por la
simulacin y los valores reales dados por el
registrador de fallas, se puede concluir que
definitivamente la falla que ocurri durante el
evento fue bifsica aislada de tierra, debido a que
representa con fidelidad las formas de ondas y
valores de tensin y corriente en la simulacin
realizada en ATPDraw.
Fuente: Elaboracin Propia
Figura 16. Corrientes en el extremo de Esmeralda
5. CONCLUSIONES Y RECOMENDACIONES
Las fallas asociadas a apertura de puentes o
ruptura de conductores representan un gran
JORNADAS TCNICAS ISA - 2012
reto para los ingenieros de operacin, debido a
que para este tipo de fallas los mtodos de
localizacin por impedancia pueden generar
errores que afectan la disponibilidad de la
lnea, en este artculo adems de presentar
una caracterizacin de la falla ocurrida en el
circuito Esmeralda Yumbo 2, se realiz
nfasis en la localizacin de falla mediante
corriente capacitiva, el cual para el caso de
ISA ha generado inspecciones en la lnea con
este tipo de fallas ms efectivas que los
mtodos de localizacin por impedancia.
20, NO. 3, JULY 2005, pp. 2038-2044
[6] ISA, Registro De Falla De La Subestacin
Esmeralda. 26/12/2010.
[7] ISA, Registro De Falla De La Subestacin
Yumbo. 26/12/2010.
[8] ISA, Metodologa para el mejoramiento del
proceso anlisis de perturbaciones mediante la
incorporacin de herramientas de simulacin
de sistemas de potencia para el equipo de
anlisis operativo de ISA. 2007.
En los sistemas de potencia ocurren diferente
tipo de fallas a diario, sin embargo fallas como
la presentada en este artculo son poco
comunes. Las redes de transmisin cada vez
ms viejas y con mayores flujos de potencia,
hacen ms propenso los circuitos a este tipo
de falla. Este hecho hace que la simulacin de
ese tipo de fallas para su caracterizacin sea
de vital importancia para la operacin del
sistema de potencia y en este artculo hemos
planteado
una
metodologa
para
la
caracterizacin de la falla mediante el
programa de simulacin ATPDraw con
resultados satisfactorios.
Sadul Segundo Urbaez Gonzalez
La caracterizacin de la falla mediante el
ATPDraw comprob que se trataba de una
falla bifsica aislada, no simtrica, lo cual
representa todo un reto para las protecciones
elctricas. Por esta razn se recomienda
incluir este tipo fallas dentro de las pruebas a
rels de proteccin.
Diego Alejandro Tejada Arango
Naci en Fonseca, Colombia. Se gradu de
Ingeniera Elctrica en la Universidad Industrial de
Santander (UIS) - Bucaramanga, Colombia, en
1994, y termin sus estudios de maestra en
Potencia Elctrica en La Universidad Industrial de
Santander- 1998. Empleado de Interconexin
Electrica S.A E.S.P desde 1997, donde ha
desempeado los cargos de Analista de Operacin
y actualmente desempea el cargo de Especialista
de Operacin de la Direccin Gestin de la
operacin de la Subgerencia Tcnica de
Interconexin Elctrica S.A E.S.P
Naci en Medelln, Colombia. Se gradu de
Ingeniera Elctrica en la Universidad Nacional de
Colombia sede Medelln y como Especialista en
Transmisin y Distribucin de Energa Elctrica en
la Universidad Pontificia Bolivariana.
Ha trabajado en la Universidad de Antioquia como
profesor en simulacin de sistemas de potencia, en
Ingeniera Especializada S.A. como ingeniero
consultor, Analista Operativo en Interconexin
Elctrica S.A y actualmente se desempea como.
Especialista de Planeacin en XM. S.A.
REFERENCIAS BIBLIOGRFICAS.
[1] DOMMEL H. W., "Emtp Rule Book", Can/Am
EMTP User Group, 1995
[2] DOMMEL H. W., "Electromagnetic Transients
Program Manual (EMTP Theory Book)",
Portland, OR: Bonneville Power Administration,
1986
[3] FURNAS,
"Transitrios
Elctricos
Coordenao De Isolamento ", 1987
[4] HIDALEN H. K., "Graphical Preprocessor To
ATP", 2008
[5] MARTNEZ J. A., GUSTAVSEN B., DURBAK
D., "Parameter Determination For Modeling
System TransientsPart I: Overhead Lines",
IEEE Transactions On Power Delivery, VOL.
9
