33

CAPTULO 3

SOFTWARE Y MODELACIN DE COMPONENTES

En estudios de protecciones los elementos modelados deben ser lo ms cercanos a la realidad, por tal motivo se debe poner especial atencin en el desarrollo de los datos para los elementos que van a estar involucrados en las simulaciones, adems de usar un software adecuado para dichos estudios. El software ATP es una de las mejores herramientas para tales propsitos, pero no existe la literatura suficiente como para, a partir de ella, saltarse los detalles de la modelacin de los elementos usados en este trabajo y simplemente ingresar los datos. Por tal motivo, este captulo tratar con detalle la preparacin e ingreso de datos, especialmente para los transformadores de potencia y corriente de la Subestacin Pomasqui, en los elementos que posee ATP.

3.1 SOFTWARE USADO

3.1.1 ATP (Alternative Transients Program) [11] [12] El ATP, que en espaol es Programa Alternativo para Transitorios

Electromangnticos, es la versin no comercial del EMTP (Electromagnetic Transients Program), aunque requiere una licencia. Fueron desarrollados para realizar estudios de transitorios electromagnticos como electromecnicos, dando mayor nfasis a los primeros. Con el tiempo se fueron agregando subrutinas que ampliaban su uso a modelaciones ms complejas como: modelos detallados de transformadores de potencia, lneas y cables, sistemas de control y ltimamente creacin de nuevas herramientas realizadas por el usuario mediante lenguaje de programacin MODELS, propio para ATP, as como incluir dentro del programa MODELS, programas externos tales como C++, MATLAB, FORTRAN, etc.

34

La escritura de los modelos en ATP se la realiza llenando religiosamente bien las tarjetas correspondientes a cada elemento del modelo, as como de los nodos que los interconectan, resultando una manera muy rgida y complicada. Por tal motivo, Kristian Hoidalen desarroll, y ha venido actualizando, el ATPDraw, que es un preprocesador grfico en el cual se realizan archivos CAD (dibujos) de cada elemento que posee ATP, interconectndolos y formando as una representacin grfica del sistema que se va a modelar. Los elementos poseen una ventana en la cual se ingresan los datos requeridos por el ATP; estos datos son orientados a su posicin correcta en las tarjetas del ATP, evitando as que se cometa errores, aunque cierto tipos de elementos como BCTRAN, LCC, XFORMER, etc., producen un archivo de salida de extensin PCH el cual es ledo por el ATP.

El ATP est compuesto por los siguientes archivos: *.atp.- Es el archivo propio de ATP, formado por tarjetas. *.lis. *.[Link] un archivo texto de salida. Sirve para observar el clculo de Es un archivo de salida y de entrada. Salida porque hace un parmetros y para leer el tipo de error cuando se los comete. tratamiento de los datos que se ingresan en el interfaz grfico, para obtener los datos que necesita el ATP. Entrada porque dichos datos necesarios son orientados a los lugares exactos en las tarjetas del ATP. *.pl4.- Es un archivo de salida que se usa para realizar los grficos de las seales seleccionadas, ya sea de voltaje, corriente, TACS o MODELS, que son las que ofrece el ATP. Hay que sealar que el ATP solamente ofrece las seales de salida en funcin del tiempo.

La forma en que ATP resuelve las ecuaciones que describen al modelo es usando la Regla Trapeziodal, por lo que hay que realizar el ajuste del paso del tiempo de acuerdo al estudio que se vaya a realizar.1

Para una visin ms amplia del ATP y ATPDraw consultar las referencias [11] y [12].

35

ATPDraw tambin posee un comando para insertar un postprocesador grfico, usualmente el PlotXY, el cual nos sirve para realizar las grficas de los archivos de salida PL4; el PlotXY adems cuenta con una herramienta que sirve para hacer un estudio, aunque no muy detallado, de los armnicos de dichas seales.

3.1.2 TOP (The Output Processor) Es un postprocesador para realizar grficos de alta calidad. Al igual que ATP, es un software libre creado por Electrotek Concepts, Inc. Es capaz de leer una gran variedad de archivos de salida para grficos, entre ellos archivos PL4. Adems posee una gran variedad de operaciones entre grficos, as como un algoritmo para realizar el filtrado de armnicas, filtro pasabajos, pasaaltos y pasabandas, etc. [13] En el ATPDraw se puede incluir el llamado al TOP por medio de la herramienta Command Edit. Este software nos servir para exportar los grficos de salida PL4, en archivos tipo Texto que pueden ser ledos por otros programas tales como Matlab, Mathcad o Excel; adems brinda archivos de salida en formato Comtrade que es til cuando se quiere comparar seales reales con las seales simuladas.

3.1.3 MATHCAD 14.0 Es un software para matemticas pensada para clculos tcnicos de ingeniera. Mathcad proporciona todas las capacidades resolutivas, la funcionalidad y la solidez necesarias para el clculo, la manipulacin de datos y el diseo industrial [14]. Una de las ventajas que presenta es al momento de escribir las ecuaciones, las cuales se editan en forma convencional, es decir, tal como vienen escritas en el papel; adems, se puede colocar texto en cualquier parte de la hoja de trabajo, sin que afecte a los procesos matemticos. Por otra parte, trae un men de ayuda con muchos ejemplos realizados, entre ellos un algoritmo para la Transformada Discreta de Fourier.

36

Para el caso especfico de este trabajo, Mathcad es usado para la graficacin de corrientes versus corrientes; dichas corrientes son tomadas desde los archivos de salida del TOP, cuya formato es Texto, y que son ledos por Mathcad. Esto es de gran utilidad, sobretodo para que los fenmenos obtenidos en las seales de salida de ATP sean visualizados de mejor manera, y por lo tanto resulte ms didctico el trabajo.

3.2 MODELACIN DEL AUTOTRANSFORMADOR DE POTENCIA DE LA SUBESTACIN POMASQUI (AT POMASQUI)

Es una de las partes ms importantes a modelar, en especial su rama de saturacin e histresis, las cuales cumplen un papel crucial en la aparicin de las corrientes Inrush. Por este motivo se va a dar un tratamiento detallado de su obtencin e insercin en el modelo de transformador usado en ATP, a partir de los datos de prueba de vaco, como de cortocircuito.

3.2.1 DATOS DE PLACA DEL AT POMASQUI Los datos de placa, tomados del reporte de pruebas Transelectric S.A., estn resumidos en la tabla siguiente: proporcionados por

RegistrodeControl

SIEMENS
EQUIPO:

HOJADEACOMPAAMIENTO DELOSENSAYOSFINALES
TIPO: FRECUENCIA: 60Hz 25Z27735

MLUN8154 AUTOTRANSFORMADOR3 TensinNominal: Norma: PCQ: Potencia: ANSIC57.12.90 180/240/300MVA 230/13813.8kV CLIENTE: TRANSELECTRIC

Tabla 3.1. Datos de placa del autotransformador de la S/E Pomasqui

37

En el Anexo 1 se encuentran los datos de pruebas de vaco y de cortocircuito del autotransformador de Pomasqui, tomados del reporte de pruebas [17]. Para la modelacin del transformador de potencia, el ATP tiene tres modelos de transformador trifsicos, a los cuales se puede agregar una rama externa para simular la rama de magnetizacin y, por consiguiente, los efectos de saturacin, histresis, flujo residual y flujo en estado estable. Dichos modelos son: Saturable 3 phase, XFORMER y BCTRAN. Se trabajar con el modelo BCTRAN.

3.2.2 MODELO BCTRAN La subrutina BCTRAN es un preprocesador de ATP el cual da como resultado la modelacin matricial del transformador de potencia a partir de sus valores nominales, as como de los datos de pruebas de vaco y cortocircuito. Tambin presenta la factibilidad de conectar una rama externa que modele la saturacin del ncleo. Una de las particularidades de este modelo de ATP es que, mediante un archivo de salida tipo .lis, presenta las matrices [A] [R], o las matrices [R] [L]; adems los valores de resistencia e inductancia de magnetizacin. Los datos necesarios que sern ingresado al modelo son: los voltajes primarios, secundarios y terciarios con sus respectivas potencias, tipo de conexin, tipo de ncleo, datos de pruebas de vaco y cortocircuito

3.2.3 RAMA DE MAGNETIZACIN: SATURACIN Y LAZO DE HISTRESIS Las pruebas de vaco se realizan colocando una fuente de voltaje variable en los terminales de menor voltaje y, empezando desde un valor igual a un porcentaje inferior al cien por ciento del voltaje nominal, se aumenta hasta un porcentaje mayor al cien por ciento; adems se toman los datos de las corrientes de cada fase y se obtiene un promedio. Para la rama de saturacin, el ATP tiene una herramienta llamada SATURATION, que sirve para obtener valores de corriente y flujo necesarios para modelar externamente la saturacin del ncleo del transformador. En esta librera se deben ingresar los valores de corriente rms (pu)

38

y los de voltaje rms (pu) obtenidos de las pruebas de vaco, pero referentes a una sola bobina, la cual se supondr representativa de las otras dos siempre y cuando los bobinados sean hechos con conductores de la misma dimensin e igual nmero de vueltas. Por esta razn hay que efectuar un tratamiento a los valores de corriente obtenidos en la prueba de vaco, pues al ser realizada en el terciario, que tiene una conexin delta, las corrientes de fase no son iguales a las corrientes de lnea; los valores de voltaje se mantienen.

If =

Il 3

(3.1)

V f = Vl

(3.2)

Basados en las pruebas de vaco realizadas al transformador de Pomasqui, y usando las ecuaciones (3.1) y (3.2) se obtienen los siguientes valores de voltajes y corrientes por fase:
Irms(A) 0,629312 0,773649 0,912213 1,03923 1,14893 1,28749 2,03805 4,68808 10,5771 20,1264 Vrms(kV) 6,805 8,282 9,653 11,033 12,412 13,105 13,804 14,476 15,183 15,856

Tabla 3.2. Resultados obtenidos a partir de datos de pruebas de vaco del AT Pomasqui

La librera SATURATION exige que las corrientes y voltajes estn en pu, en sus propias bases. Del reporte de pruebas se observa que el voltaje base es 13.8 kV y la potencia base es 180 MVA. Entonces, como necesitamos solo para una bobina,

39

la potencia base dividimos para tres; con estos resultados calculamos la corriente base. VB = 13.8 kV SB = 60 MVA

I B ( A) =

S B 1000 (kVA) 13.8 (kV )

I B = 4347.83 A
De esta manera se obtienen los valores de voltaje rms (pu) y corriente rms (pu), que se muestran en la tabla 3.3.
Irms(pu) 1,45E04 1,78E04 2,10E04 2,39E04 2,64E04 2,96E04 4,69E04 1,08E03 2,43E03 4,63E03 Vrms(pu) 0,49311594 0,60014493 0,69949275 0,79949275 0,89942029 0,94963768 1,00028986 1,04898551 1,10021739 1,14898551

Tabla 3.3. Corriente y voltaje rms (pu) de pruebas de vaco para el AT Pomasqui

Los valores de la tabla 3.3 estn referidos al terciario. En este trabajo se va modelar la rama de magnetizacin conectada en el primario, por lo que dichos datos debern cambiar de base. Este clculo lo puede realizar el Saturation, para ello simplemente debemos cambiar los valores de voltaje base y potencia base en las tarjetas para dichos datos. Entonces se tiene: VB = 132.79 kV SB = 100 MVA

40

Figura 3.1. Ingreso de datos en SATURATION para obtener la curva de saturacin referida al primario del AT Pomasqui

Obtenindose los siguientes resultados:

Figura 3.2. Archivo de salida tipo punch para la curva de saturacin del AT Pomasqui

Para modelar el lazo de histresis, ATP posee una herramienta llamada HYSTERESIS cuyos datos de entrada son los penltimos valores de corriente y flujo [11], (donde se superponen los grficas de flujo e histresis), del archivo de salida punch de la rutina SATURATION, figuras 3.2 y 3.3.

41

Lazo de Histrresis
110
3

500

Flujo (Wb-T)

10

10

500

Curva de Saturacin Histresis + Histresis 110

3

Corriente (A)

Figura 3.3. Lazo de histresis para el AT Pomasqui

Cabe sealar que HYSTERESIS solamente presenta la parte inferior del lazo de histresis, ATP internamente genera el lazo superior [11].

Figura 3.4. Datos de salida de HYSTERESIS para el AT Pomasqui

42

Figura 3.5. Lazo inferior de histresis calculado por ATP para el AT Pomasqui

Estos datos se ingresan en los elementos externos no lineales de ATP; se eligi trabajar con el elemento L(i) Type 96 ya que es una herramienta que nos permite, adems, tener un control del flujo inicial y el flujo residual, y que nos sern de gran utilidad. En serie con cada elemento se coloca una resistencia de igual valor a la del bobinado primario, esto para tener en cuenta la atenuacin de la onda de intensidad en la simulacin del transitorio de conexin del transformador [16]. Para estudios de energizacin, el acoplamiento magntico entre los circuitos y la caracterstica no lineal de la rama de magnetizacin son las ms importantes2.

3.3 MODELACIN CORRIENTE

DE

LOS

TRANSFORMADORES

DE

El Transformador de Corriente se modela usando Saturable 1 phase, el cual nos brinda la posibilidad de ingresar los datos de prueba de los TCs que se realizan en fbrica. Este elemento presenta un modelo igual al de la figura 3.7, con N=1, es decir, con un solo circuito.

Manual de ATP Draw, pg 213

43

Figura 3.6. Modelacin del transformador por Saturable de ATP [11]

De acuerdo a lo dicho en 2.2.2 y en [9], se van a tener las siguientes simplificaciones: El secundario del TC se modela en el circuito 1, (winding 1 de la figura 3.6) En el circuito 2, (winding 2 de la figura 3.6), la resistencia R2 es igual a cero, la inductancia L2 es ms grande que cero, por ejemplo 10E-6. En el circuito 1, L1 es igual a cero, debido a que el secundario del TC, segn la norma ANSI, es de clase C y la reactancia de fuga es muy pequea. El valor de R1 es igual a la resistencia del secundario del TC. Adems, resistencias separadas se pueden agregar para modelar los conductores de conexin entre el TC y el burden y para modelar la carga que representa el rel diferencial. Debido a que RMAG es muy grande, se colocar el valor 0 ya que en ATP el numero cero indica una cantidad infinita. Para incluir la saturacin en el modelo, se usar la librera Saturation de ATP para obtener los puntos de corriente vs. flujo que se incluirn va comando INCLUDE, a partir de los datos de prueba.

44

Tambin se va a despreciar la histresis. Los datos de la saturacin se los calcula de igual manera que para el transformador de potencia, usando la subrutina SATURATION. Adems, hay que tener en cuenta la razn de transformacin del TC para el ingreso de Vrp y Vrs. Se proceder de la siguiente manera:

I nom prim TC I nom sec TC

Vrp Vrs

(3.3)

As, por ejemplo, para el TC colocado en el primario del transformador de Pomasqui, cuya relacin es 800/5, se tiene:

Vrp 800 160 = = Vrs 5 1

En el Anexo 1 se encuentran los datos de los Transformadores de Corriente para proteccin, los datos de las curvas de saturacin y el burden de cada uno de ellos, para el clculo de su precisin.

3.3.1 CONEXIN DE LOS TCs El modelo de rel diferencial a desarrollarse emular a los rels de ltima tecnologa, basados en microprocesador. Con este tipo de rel no es necesario realizar las conexiones en delta de los transformadores de corriente; esta conexin sirve para, por un lado corregir los desfasamientos de ngulos y las correcciones de la magnitud de la corriente debido a la presencia de una conexin delta en el transformador de potencia, y por otro para eliminar las componentes de secuencia cero, las cuales quedaban encerradas el circuito delta. En los rels modernos, dicha compensacin de ngulo, de fase y eliminacin de secuencia cero se realizan mediante algoritmos matemticos, as como tambin para la discriminacin de fallas internas, externas, energizaciones y sobreexcitacin. Debido a lo anterior, los secundarios de los TCs van a tener la conexin Y. En los circuitos secundaros se adicionarn unas resistencias que modelan la resistencia

45

de los conductores que unen los transformadores de corriente con el rel diferencial, adems una carga en Y que modela al rel mismo, y por ltimo unas resistencias que modelan las conexiones al neutro comn de los transformadores de corriente primarios y secundarios, y del rel diferencial.

3.4 ELEMENTOS ADICIONALES DEL SISTEMA A MODELAR

A ms de los elementos anteriores, se va a modelar el Sistema Nacional Interconectado, S.N.I., mediante circuitos equivalentes de cortocircuito, tanto de Colombia realizado en la S/E Jamondino, como del S.N.I. en la barra a 230 kV en la S/E Santa Rosa. Tambin se incluir en el modelo las lneas de transmisin Jamondino Pomasqui, Santa Rosa Pomasqui y la carga del autotransformador de Pomasqui.

3.4.1 EQUIVALENTES DE COLOMBIA Y DEL S.N.I. Para estos elementos se opt por usar el software DigSilent Power Factory 13.2, pues es el ms adecuado para realizar este trabajo. Los resultados que ofreci este programa son los siguientes valores de resistencia e inductancia equivalentes: a) Colombia: Subestacin Jamondino, barra de 230 kV
R sec 0 (/m) L sec 0 (/m) R sec + (/m) L sec + (/m) 7.553 71.893 6.269 70.882

Tabla 3.4. Impedancia equivalente de Colombia en S/E Jamondino

46

b)

Sistema Nacional Interconectado: Subestacin Santa Rosa, barra de

230 kV

R sec 0 (/m) L sec 0 (/m) R sec + (/m) L sec + (/m)

0.923 14.298 2.403 22.183

Tabla 3.5. Impedancia Equivalente del SNI en S/E Santa Rosa

3.4.2 MODELACIN DE EQUIVALENTES DE COLOMBIA Y DEL S.N.I. Para la modelacin en ATP se usarn fuentes de voltaje trifsicas balanceadas de 230 kV rms cada una. Al ingresar el valor de la fuente hay que considerar que ATP trabaja con valores pico y por fase, considerando adems una conexin Y de la fuente. Entonces el valor a ingresar ser:

Vp f n =

Vrms (V ) 2 3 230 2 = 187794 (V) 3

Vp f n =

Fuente trifsica 230 kV rms R0 R+ L0 L+

S/E Jamondino 230 kV rms

Equivalente Colombia

Figura 3.7. Conexin del equivalente de Colombia en S/E Jamondino

47
S/E Santa Rosa 230 kV rms

Fuente trifsica 230 kV rms R0 R+ L0 L+

Equivalente Ecuador

Figura 3.9. Conexin del Equivalente del S.N.I. en S/E Santa Rosa

En ATP las resistencias e inductancias se modelarn usando el elemento LINESY_3, que considera acoplamiento de secuencias.

3.4.3 MODELACIN DE LNEAS DE TRANSMISIN Las lneas de transmisin Jamondino Pomasqui y Santa Rosa Pomasqui operan a 230 kV, con longitudes de 212.2 km y 46 km, respectivamente. Ambas son trifsicas, doble circuito, cada una con dos hilos de guarda. No se tomar en cuenta los aterramientos por segmentos.

Figura 3.9. Ingreso de Datos de Lneas de Transmisin en ATP

48

Los datos elctricos y fsicos de los conductores, y la configuracin de las torres fueron proporcionados por Transelectric S.A. El tipo de conductores de fase como de guardia son los mismos para las dos lneas, variando solamente en la distancia. Se usar la subrutina LINE CONSTANTS de ATP, la cual da como salida un modelo de lnea de parmetros distribuidos a partir de sus datos de diseo. Para una comprensin del uso de esta subrutina, se recomienda observar las referencias [11], [22] y [26]. Los datos de entrada de la geometra de las torres usadas y los datos fsicos de los conductores, en el interfaz grfico ATPDraw, son los mostrados en la figura 3.10.

3.4.4 MODELACIN DE LA CARGA DEL TRANSFORMADOR DE POMASQUI La simulacin de falla externa, en la barra de 138 kV de la S/E Pomasqui, se realizar a la carga actual a la que est alimentando el transformador. El valor, inicialmente en la forma de potencia activa y reactiva tomada de datos de flujo realizado por Transelectric S.A., se pasar a corriente. En estudios de protecciones se suele afectar la corriente de carga con un porcentaje, usualmente del 20%, para tomar en cuenta futuros incrementos en la carga. Luego de esto, y con el valor del voltaje secundario del transformador, se obtendr la impedancia de la carga.
DATOSDEENTRADADECARGA RESULTADOSPARALACARGA

Pactiva(MW)
150.00

Qreactiva(MVar) [Link](kV)
80.00 138.00

Corriente(A)
709.997+378.665i

Impedancia()
123.555+65.896i

[Link].

Tabla 3.6. Carga conectada al autotransformador de Pomasqui3

Datos proporcionados por Transelectric S.A. Divisin de Estudios y Plan de Expansin.

49

El elemento de ATP que modela la carga es el RLCY3.

3.4.5 MODELACIN DE OTROS COMPONENTES ADICIONALES Existen otros componentes adicionales muy tiles para la simulacin de energizaciones y fallas, tales como interruptores y ramas para conexiones de fallas monofsicas, fase-fase, fase-fase-tierra y trifsicas. Se explicar brevemente los interruptores y la falla monofsica.

a. INTERRUPTORES De la teora presentada en 2.1.3, se observa que el tiempo de cierre del interruptor que energiza al transformador es determinante para la aparicin de corrientes Inrush. Por tal motivo, se debe seleccionar un tiempo en la onda de voltaje de cualquiera de las tres fases, en la barra de Pomasqui, cuando realice el cruce por cero; este tiempo ser ingresado como tiempo de cierre. En ATP existe un tipo de interruptor controlado por tiempo denominado SWIT_3XT, que se usar para tal propsito. b. RAMA DE CONEXIN DE FALLA En esencia, es una resistencia de falla en serie con un interruptor controlado por tiempo, que se conecta a un nodo mediante un elemento de ATP, denominado SPLITER.