Calidad de La Energia
Calidad de La Energia
Calidad de La Energia
Director
GABRIEL ORDOEZ PLATA
Ingeniero Electricista, PhD.
Codirector
VICTOR BARRERA NEZ
Ingeniero Electricista, MIE (C)
Agradecimientos
Los autores expresan sus agradecimientos:
A Gabriel Ordez Plata, Doctor Ingeniero Industrial en Ingeniera Elctrica, Profesor Titular de la
Universidad de Santander y Director del trabajo de grado, por brindarnos su apreciable
colaboracin e invaluables conocimientos, experiencia y sobretodo por su calidad humana.
A Hermann Ral Vargas Torres, Doctor Ingeniero Industrial en Ingeniera Elctrica, Profesor Titular
de la Universidad de Santander, por sus valiosas ideas en la automatizacin de la herramienta
SimRF.
A Jorge Andrs Cormane, Magster en Ingeniera Elctrica (C), por brindarnos su trabajo como
referencia, sus conocimientos de programacin en el campo de la interaccin entre las plataformas
informticas ATP y Matlab y por ser uno de los orientadores de Ronald.
A Juan Jos Mora, Ph.D (C) por su colaboracin en la revisin constante del proyecto y el
suministro de informacin para la elaboracin de la herramienta SimRF.
Este trabajo de grado muestra el estudio realizado para poder caracterizar dos circuitos de
distribucin de media tensin en ATPDraw y el proceso que sigue una herramienta computacional
desarrollada en Matlab, para simular y analizar de manera automtica distintas fallas de
cortocircuito en cada uno de los sistemas elctricos. Los circuitos modelados fueron provistos por
la empresa Comercializadora de Energa S.A. (CODENSA). La herramienta informtica fue basada
en la plataforma computacional ATP para realizar la simulacin de las fallas de cortocircuito y en
Matlab para efectuar el anlisis de resultados de la simulacin y sistematizacin del proceso de
simulacin. El proceso que lleva a cabo la aplicacin informtica se resume en el modelado del
sistema en ATPDraw, el copiado del modelo original con los cambios en la localizacin, tipo y
resistencia de falla, la simulacin de estas copias, la conversin de formato de los resultados de las
simulaciones y finalmente el anlisis de estos resultados en Matlab. El estudio de calidad se basa
en el anlisis de los resultados de las simulaciones, este anlisis consiste en el clculo de ciertos
parmetros elctricos a partir de las seales de tensin y de corriente medidas en la subestacin
en la cabecera de cada circuito. Estos parmetros se conocen como descriptores y algunos son el
tiempo de duracin y frecuencia del transitorio en la seal de tensin, variacin en la magnitud de
la corriente eficaz entre el estado prefalla y postfalla, entre otros.
1
Trabajo de investigacin
2
Facultad de Ingenieras Fsico-Mecnicas, Ingeniera Elctrica, Director: PhD. Gabriel Ordez Plata
vii
KEY WORDS: Power systems, electric distribution, electric power quality, Alternative Transients
Program, EMTP, ATPDraw, programming on Matlab.
This thesis shows the study for the characterization of two medium voltage distribution networks on
ATPDraw and the process a software developed in Matlab follows in order to automatically simulate
and analyze short-circuit faults on both electric systems. The modeled circuits were provided by
Comercializadora de Energa S.A. (CODENSA). The software developed was based on ATP
platform for simulations and on Matlab to analyze simulation results and to systematize the
simulation process. The process the software follows can be sum up by the modeling of the
distribution network on ATPDraw, the copy of the original model with changes of fault location, type
and resistance, the simulation of these copies, the format conversion of the simulation results and
finally the analysis of these simulation results on Matlab. The study of power quality is based on the
analysis of the simulation results, this analysis consists on the calculation of certain electric
parameters of the voltage and current signals measured on the main substation of the circuit. These
parameters are known as descriptors and some of them are the duration and frequency of voltage
signal transients, variations on the magnitude of RMS current between values before and after the
fault takes place, among others.
3
Research
4
Facultad de Ingenieras Fsico-Mecnicas, Ingeniera Elctrica, Director: PhD. Gabriel Ordez Plata
viii
TABLA DE CONTENIDO
1. INTRODUCCIN ............................................................................................. 6
2. CALIDAD DE LA ENERGA ELCTRICA ........................................................ 8
2.1. TRANSITORIOS ELECTROMAGNTICOS ........................................... 10
2.2. DESCRIPTORES.................................................................................... 12
2.2.1. Caracterizacin de los huecos de tensin ....................................... 12
2.2.2. Caracterizacin de las seales transitorias...................................... 15
2.3. FALLAS DE CORTOCIRCUITO ............................................................. 17
3. MODELADO DE ELEMENTOS DE CIRCUITOS DE DISTRIBUCIN........... 19
3.1. FUENTE.................................................................................................. 19
3.2. ALIMENTADORES O LNEAS DE DISTRIBUCIN ............................... 20
3.2.1. Parmetros concentrados ................................................................ 20
3.2.2. Parmetros distribuidos constantes ................................................. 21
3.2.3. Parmetros distribuidos dependientes de la frecuencia................... 21
3.3. CARGAS................................................................................................. 23
3.4. PROTECCIONES ................................................................................... 24
3.5. FALLAS DE CORTOCIRCUITO ............................................................. 24
4. SELECCIN DE MODELOS.......................................................................... 25
4.1. FUENTE.................................................................................................. 25
4.2. ALIMENTADORES O LNEAS DE DISTRIBUCIN ............................... 27
4.3. CARGAS................................................................................................. 31
4.4. PROTECCIONES ................................................................................... 33
4.5. FALLAS DE CORTOCIRCUITO ............................................................. 34
5. SISTEMATIZACIN DE LA SIMULACIN .................................................... 36
5.1. SELECCIN DE PLATAFORMAS INFORMTICAS.............................. 36
5.2. INFORMACIN GENERAL DE LAS PLATAFORMAS INFORMTICAS
SELECCIONADAS ............................................................................................ 37
5.3. INTERACCIN ENTRE PLATAFORMAS INFORMTICAS................... 38
5.3.1. Modelado de circuitos en ATPDraw y generacin del cdigo fuente 39
5.3.2. Copia de los ejecutables y estructura del rbol de carpetas............ 41
5.3.3. Modificacin del cdigo fuente......................................................... 42
5.3.4. Archivo de Instrucciones Batch_SimRF.bat y copia de la plantilla... 44
5.3.5. Archivo de instrucciones conv_RF.bat y conversin a archivos .mat
46
5.4. PROCEDIMIENTO DETALLADO PARA EJECUTAR LA HERRAMIENTA
INFORMTICA .................................................................................................. 47
6. SIMPLIFICACIN Y CARACTERIZACIN DE LOS CIRCUITOS REALES .. 50
7. RESULTADOS............................................................................................... 61
7.1. CIRCUITO SF-15 .................................................................................... 63
7.1.1. Circuito SF-15 con protecciones ...................................................... 63
7.1.2. Circuito SF-15 sin protecciones ....................................................... 69
7.2. CIRCUITO UM-33 ................................................................................... 72
7.2.1. Circuito UM-33 con protecciones ..................................................... 73
ix
LISTA DE FIGURAS
LISTA DE TABLAS
LISTA DE ANEXOS
1. INTRODUCCIN
Conociendo ahora lo que se quiere llevar a cabo y con el fin de centrar al lector en
el tema, se tratarn los siguientes tpicos en los captulos subsiguientes.
7
5
Norma Tcnica Colombiana NTC 5000
9
8
Comisin Reguladora de Energa y Gas, CREG
9
IEEE 1159 de 1995
11
Los transitorios oscilatorios son el cambio repentino de una seal a una frecuencia
diferente a la de suministro. Lo anterior se presenta en condicin de estado
estacionario de la tensin, la corriente, o ambas, donde los valores instantneos
cambian de polaridad rpidamente. Estos cambios se caracterizan por su
frecuencia, duracin y magnitud.
El transitorio tipo impulso por otra parte, se define como el cambio sbito a una
frecuencia distinta a la de suministro. Esto se presenta en condicin de estado
estacionario para la tensin, la corriente o ambos, y es unidireccional en polaridad.
Se caracterizan normalmente por los tiempos de subida y cada de la seal. Los
transitorios tipo impulso normalmente son causados por descargas atmosfricas y
pueden provocar transitorios oscilatorios por induccin electromagntica.
2.2. DESCRIPTORES
Los descriptores anteriores pueden presentarse en tres posibles zonas del hueco
de tensin:
Para los huecos de tensin no monofsicos las formas de onda son similares a las
de la figura 3.
Duracin de cada de tensin
trifsica
Magnitud de
cada de
tensin trifsica
Fase A
Fase B
Fase C
10
Greenwood Allan, Electrical Transients in Power Systems. Electric Power Engineering Department,
Troy, New York. Second edition.
11
Karimi, Masoud; Mokthari, Hossein; Iravani, Reza. Wavelet based on line disturbance detection for power
quality applications. IEEE Transactions on Power Delivery, vol. 17,No. 1, January 2002.
16
Informacin de la firma del sistema: Para cada fase se tiene un vector fila
que contiene la relacin corriente-tiempo (I/t) para los intervalos en que
est presente la falla, adems de esto, los tiempos entre disparos de las
protecciones, cuando estas actan para aislar una falla. Tambin contiene
informacin como la posicin de la muestra y tiempo de disparo.
Las fallas que comnmente ocurren en los sistemas de energa elctrica son las
de cortocircuito.
Existen varias formas de que ocurran este tipo de fallas, de las cuales se
enumeran para un sistema trifsico la de fase a tierra, fase a fase, fase-fase a
tierra, todas las fases y todas las fases a tierra.
Fase A tierra
Fase B tierra
Fase C tierra
Fase A fase B
Fase B fase C
Fase C fase A
Fase A fase B tierra
18
En este mtodo, cada tipo de falla tiene un circuito equivalente con el fin de
realizar el anlisis del sistema de transmisin de energa elctrica en el cual
ocurri la falla. Si el lector desea profundizar ms en este tema puede referirse al
escrito de Grainger/Stevenson 12 .
12
Grainger/Stevenson, Anlisis de Sistemas de Potencia. McGraw Hill.
19
Fuente
Alimentadores o lneas de distribucin
Cargas
Protecciones
Fallas de cortocircuito
3.1. FUENTE
El modelo considerado fue brindado por el software ATPDraw, el cual simula una
fuente trifsica de tensin o de corriente, aterrizada, con forma de onda sinusoidal
de nica frecuencia, en la cual se especifican los valores de amplitud, frecuencia,
desfase, tiempo de comienzo y tiempo de finalizacin, para representar este
elemento.
20
13
Orlando Hevia, Comparacin de los modelos de lnea del ATP
22
Para los dos modelos anteriores cabe decir que presentan el mismo inconveniente
que aquellos modelos de parmetros distribuidos constantes con respecto al
tiempo de paso, donde se exige que este sea menor que el tiempo de
propagacin.
14
Taku Noda, Akihiro Ametani. User Instructions of Noda Setup in ATP. Julio 1998.
23
3.3. CARGAS
3.4. PROTECCIONES
4. SELECCIN DE MODELOS
4.1. FUENTE
Poseer una potencia fija: Esto se debe a que en la ocurrencia de una falla
de cortocircuito, el comportamiento de la fuente ser el incremento de la
magnitud de corriente y el decrecimiento de la magnitud de la tensin.
15
CODENSA, Comercializadora de Energa S.A. E.S.P.
26
Los resultados obtenidos por medio de esta prueba muestran que el modelo de la
fuente, garantiza cadas de tensin en el evento de una falla de cortocircuito. Para
mostrar esto, en la figura 5 se muestra la tensin eficaz durante un periodo de
1seg, para el caso de una falla trifsica a los 250ms y se despej a los 750ms en
la ubicacin mostrada en la figura 4.
27
Tensin 20
Eficaz (V) *10 3
16
12
0
0.0 0.2 0.4 0.6 0.8 [s] 1.0
(file PRUEBAFUENTE.pl4; x-var t) t: XX0031 t: XX0039 t: XX0047
t (ms)
Figura 5. Tensin eficaz en terminales de la fuente cuando ocurre una falla.
Para los circuitos de media tensin estudiados, se encontr que para la tensin de
11 kV el nivel de falla es de 200MVA 16 . El nivel de tensin fue extrado de los
datos de los transformadores de las subestaciones principales como se ve en el
Anexo 1.
16
Math H.J. Bollen. Understanding Power Quality Problems, Voltage Sags and Interruptions. IEEE Press. 1960
28
Haciendo entonces una revisin del documento de Hevia, se descartan por este
mismo problema los modelos de parmetros distribuidos constantes, lo cual deja
el modelo Pi y el modelo de T. Noda.
17
Orlando Hevia, Comparacin de los modelos de lnea de ATP
18
Comercializadora de Energa S.A. E.S.P.
29
Viendo esto, solo queda usar el modelo Pi, el cual presenta tiempos de simulacin
razonables de 15 s. Entonces todas las lneas de los circuitos UM-33 y SF-15
fueron modeladas utilizando el equivalente Pi de la rutina LCC (Line Cable
Constants) trifsica del ATP.
El bloque LCC de ATP permite ingresar los datos de la lnea, como lo son, la
distancia entre fases, datos de los cables y otros datos como la resistencia de
tierra. Esto se muestra en la figura 7.
log(| Z |) log(| Z |)
2.3 2.3
Modelo Implementado
Modelo Pi Exacto
a. b.
Figura 8. Barrido de frecuencia para modelo Pi.
a. Secuencia cero, b. Secuencia positiva
2.2 log(| Z |) 2.2 log(| Z |)
Modelo Implementado
Modelo Pi Exacto
1.2 1.2
Modelo Implementado
0.1 0.1
Modelo Pi Exacto
log(freq) log(freq)
-0.9 -0.9
1.8 2.6 3.5 4.4 1.8 2.6 3.5 4.4
a. b.
Figura 9. Barrido de frecuencia para modelo J.Mart.
a. Secuencia cero, b. Secuencia positiva
31
Esto no ocurre con el modelo J.Mart (figura 9), sustentando esto que el modelo de
Mart es ms apropiado para altas frecuencias que el modelo Pi. Sin embargo, al
pensar en el tiempo necesario para calcular una solucin en los circuitos de media
tensin estudiados (2hr por falla en J.Mart vs. 15s por falla en modelo Pi), se
escoge el modelo Pi.
4.3. CARGAS
As que se decidi tomar los transformadores como cargas trifsicas que trabajan
a plena carga durante el transcurso de la simulacin.
Para el caso del circuito SF-15, en donde se compensa la carga en los nodos
N006 y N051, el clculo de sta compensacin se lleva a cabo por medio de un
coeficiente encontrado en la tabla 2 19 .
Factor Factor K
de
potencia Factor de potencia deseado
inicial 0.86 0.87 0.88 0.89 0.9 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1
0.66 0.545 0.572 0.599 0.626 0.654 0.683 0.712 0.743 0.775 0.810 0.847 0.888 0.935 0.996 1.138
0.67 0.515 0.541 0.568 0.596 0.624 0.652 0.682 0.713 0.745 0.779 0.816 0.857 0.905 0.966 1.108
0.68 0.485 0.512 0.539 0.566 0.594 0.623 0.652 0.683 0.715 0.750 0.787 0.828 0.875 0.936 1.078
0.69 0.456 0.482 0.509 0.537 0.565 0.593 0.623 0.654 0.686 0.720 0.757 0.798 0.846 0.907 1.049
0.70 0.427 0.453 0.480 0.508 0.536 0.565 0.594 0.625 0.657 0.692 0.729 0.770 0.817 0.878 1.020
0.71 0.398 0.425 0.452 0.480 0.508 0.536 0.566 0.597 0.629 0.663 0.700 0.741 0.789 0.849 0.992
0.72 0.370 0.397 0.424 0.452 0.480 0.508 0.538 0.569 0.601 0.635 0.672 0.713 0.761 0.821 0.964
0.73 0.343 0.370 0.396 0.424 0.452 0.481 0.510 0.541 0.573 0.608 0.645 0.686 0.733 0.794 0.936
0.74 0.316 0.342 0.369 0.397 0.425 0.453 0.483 0.514 0.546 0.580 0.617 0.658 0.706 0.766 0.909
0.75 0.289 0.315 0.342 0.370 0.398 0.426 0.456 0.487 0.519 0.553 0.590 0.631 0.679 0.739 0.882
0.76 0.262 0.288 0.315 0.343 0.371 0.400 0.429 0.460 0.492 0.526 0.563 0.605 0.652 0.713 0.855
0.77 0.235 0.262 0.289 0.316 0.344 0.373 0.403 0.433 0.466 0.500 0.537 0.578 0.626 0.686 0.829
0.78 0.209 0.236 0.263 0.290 0.318 0.347 0.376 0.407 0.439 0.474 0.511 0.552 0.599 0.660 0.802
0.79 0.183 0.209 0.236 0.264 0.292 0.320 0.350 0.381 0.413 0.447 0.484 0.525 0.573 0.634 0.776
0.80 0.157 0.183 0.210 0.238 0.266 0.294 0.324 0.355 0.387 0.421 0.458 0.499 0.547 0.608 0.750
0.81 0.131 0.157 0.184 0.212 0.240 0.268 0.298 0.329 0.361 0.395 0.432 0.473 0.521 0.581 0.724
0.82 0.105 0.131 0.158 0.186 0.214 0.242 0.272 0.303 0.335 0.369 0.406 0.447 0.495 0.556 0.698
0.83 0.079 0.105 0.132 0.160 0.188 0.216 0.246 0.277 0.309 0.343 0.380 0.421 0.469 0.530 0.672
0.84 0.053 0.079 0.106 0.134 0.162 0.190 0.220 0.251 0.283 0.317 0.354 0.395 0.443 0.503 0.646
0.85 0.026 0.053 0.080 0.107 0.135 0.164 0.194 0.225 0.257 0.291 0.328 0.369 0.417 0.477 0.620
0.86 --- 0.027 0.054 0.081 0.109 0.138 0.167 0.198 0.230 0.265 0.302 0.343 0.390 0.451 0.593
0.87 --- --- 0.027 0.054 0.082 0.111 0.141 0.172 0.204 0.238 0.275 0.316 0.364 0.424 0.567
0.88 --- --- --- 0.027 0.055 0.084 0.114 0.145 0.177 0.211 0.248 0.289 0.337 0.397 0.540
0.89 --- --- --- --- 0.028 0.057 0.086 0.117 0.149 0.184 0.221 0.262 0.309 0.370 0.512
0.90 --- --- --- --- --- 0.029 0.058 0.089 0.121 0.156 0.193 0.234 0.281 0.342 0.484
0.91 --- --- --- --- --- --- 0.030 0.060 0.093 0.127 0.164 0.205 0.253 0.313 0.456
0.92 --- --- --- --- --- --- --- 0.031 0.063 0.097 0.134 0.175 0.223 0.284 0.426
0.93 --- --- --- --- --- --- --- --- 0.032 0.067 0.104 0.145 0.192 0.253 0.395
0.94 --- --- --- --- --- --- --- --- --- 0.034 0.071 0.112 0.160 0.220 0.363
0.95 --- --- --- --- --- --- --- --- --- --- 0.037 0.078 0.126 0.186 0.329
0.96 --- --- --- --- --- --- --- --- --- --- --- 0.041 0.089 0.149 0.292
0.97 --- --- --- --- --- --- --- --- --- --- --- --- 0.048 0.108 0.251
0.98 --- --- --- --- --- --- --- --- --- --- --- --- --- 0.061 0.203
0.99 --- --- --- --- --- --- --- --- --- --- --- --- --- --- 0.142
19
Comisin Nacional para el Ahorro de Energa, CONAE, Mxico. Compensacin del factor de potencia.
www.conae.gob.mx/wb/CONAE/CONA_438_compensacion_del_fac
33
Para las cargas a compensar en ambos nodos del circuito SF-15 se asumi un
factor de potencia de 0,85 en atraso, compensado a 0,9 en atraso, lo cual implica
una constante K=0,135. El clculo de las potencias se ve en la tabla 3.
Reactancia
Potencia aparente Potencia activa Potencia reactiva de la
capacitiva de la
NODO de la carga de la carga compensacin
compensacin
S 3 (kVA) P 1 (kW) Q 1 (kVAR)
Xc ()
N006 300 85 11,47 3776,8
N051 475 134,583 18,17 2384,3
Tabla 3. Clculos de compensacin con K=0,135
4.4. PROTECCIONES
Constantes Valores
A 0,047
B 0,183
P 0,02
Tabla 4. Constantes de la curva moderadamente inversa.
20
Westinghouse Electric Company. Applied Protective Relaying. 1976
35
Como se puede ver en la figura 10, las fases y tierra tienen un interruptor
individual que se puede controlar a partir de tiempos de cerrado y apertura.
Tambin se encuentra una resistencia de falla por fase y una a tierra. Estas
resistencias toman valores diferentes dependiendo del tipo de falla, por ejemplo en
caso de una falla de la fase A a tierra, se tiene que el interruptor que se debe
cerrar es el de la fase A en el tiempo al cual ocurre la falla (0,02 s), la resistencia
de falla RA debe tomar el valor de la resistencia de falla Rf, el interruptor de tierra
estar cerrado desde el comienzo de la simulacin y la resistencia de tierra RG
deber tener un valor muy pequeo (0,001 ).
5. SISTEMATIZACIN DE LA SIMULACIN
Las fallas en los circuitos de distribucin pueden ser simuladas con la ayuda de
sistemas computacionales tales como ATP y Matlab. Estos sistemas contribuyen a
agilizar el proceso de anlisis de fallas y a responder algunas preguntas referentes
a estudios de calidad, que absolverlas mediante el uso de otros sistemas, no solo
sera difcil y riesgoso, sino que implicara altos costos econmicos.
Las conclusiones a las que llegaron los autores del mencionado documento, se
describen en la tabla 5.
21
F.Jurado, N.Acero, J.Carpio, M.Castro. Using various computer tools in electrical transient studies.
37
De igual manera, fue necesario tener en cuenta el documento cuyos autores son
Nan Zhang y Mladen Kezunovic, en razn a que all se describe una herramienta
informtica basada en la interactividad entre Matlab y ATP para el anlisis de
fallas en sistemas elctricos 22 .
22
N.Zhang, M.Kezunovic. Implementing and advanced simulation tool for comprehensive fault analysis.
38
Con respecto al ATP, conviene resaltar que existen numerosas versiones para
diversos tipos de usuarios. Para este estudio se encontr compatible la versin del
ATPDraw 3.5 Mingw32, la cual viene provista de cuatro ejecutables para su
instalacin (en el CD de instalacin se encuentran en la ruta
/Instaladores/GNU/Mingw32/Easyinstall). Estos archivos clasificados en su orden
de ejecucin para una correcta instalacin, son los siguientes:
InstATP111.exe InstATP114.exe
InstATP111_add.exe
InstATP111_lib.exe
ATPLnchUpdate.exe
El primer archivo ejecutable a instalar puede ser una de las dos versiones, siendo
la 114 ms reciente que la 111. El ltimo archivo ejecutable es muy importante,
dado que actualiza la aplicacin ATPLauncher, la cual es base de la etapa de
simulacin de la herramienta creada en este trabajo de grado.
Estas copias de la plantilla sern luego simuladas una a una utilizando el archivo
runATP.exe por medio de otro archivo llamado Batch_SimRF.bat. Este archivo .bat
instruye al computador el orden de las simulaciones y es creado a su vez,
mediante un algoritmo programado en Matlab.
Una vez se instale el ATPDraw, el usuario deber seleccionar los modelos de los
elementos del sistema que ms se adecuen al estudio a realizar y deber disear
el circuito en la aplicacin. Luego deber almacenar este circuito grfico (archivo
con extensin .adp) y realizar una simulacin cualquiera. Esto da como resultado
la generacin de un archivo de extensin .atp en la carpeta WORK de ATP, el cual
ser usado como plantilla o cdigo fuente.
ATP MATLAB
Modelar circuitos en Copia de los ejecutables y estructura
ATPDraw (seccin 5.3.1.) del rbol de carpetas (seccin 5.3.2.)
Figura 11. Proceso automtico para la simulacin de fallas utilizando ATP y MATLAB
Al elaborar el circuito grfico se debe tener en cuenta que para ramales trifsicos y
bifsicos los nodos deben ser nombrados siguiendo el formato NXXX, en el cual
XXX representa desde el nodo 000 en bornes de la fuente, hasta el 999; para
ramales monofsicos la convencin estar dada por NXXXF en el cual XXX
representa desde el nodo 000 en bornes de la fuente, hasta el 999 y F representa
la fase donde est conectado el ramal monofsico (A, B o C).
41
Conviene precisar que si el usuario utiliza bloques LCC 23 para representar las
lneas de transmisin o de distribucin de energa, al realizar la simulacin, estos
bloques generarn archivos .atp en la carpeta PROJECT/LCC. Estos archivos
debern ser dejados all, ya que sern referenciados en el archivo de extensin
.atp.
23
Lines Cables Constants Rutine, ATP
42
Una vez es generado el archivo modelo *.adp, este crea el cdigo fuente *.atp al
ser simulado; en la figura 13 se presenta la parte del cdigo fuente que debe ser
modificada mediante la escritura de las resistencias de falla de cada fase y de
tierra.
SwA com, SwB com, SwC com y com Swg y los tiempos de conmutacin
entre la barra y el interruptor, de cada una de las fases y tierra dependiendo el tipo
de falla.
Las lneas de texto que hacen referencia a los tiempos de conmutacin deben
respetar la convencin NXXXF donde XXX representa la barra donde ocurre la
falla y F la fase donde est conectado cada interruptor del elemento fallador,
adems se incluye la conexin a tierra por medio del interruptor Swg.
Para realizar las fallas en las fases, los tiempos que van en la columna <Tclose>
(tiempo de cierre) deben ser menores que los tiempos que van en la columna
<Top/Tde> (tiempo de apertura); para las fallas a tierra se recomienda introducir
un tiempo de conmutacin en la columna <Tclose> en -1 para lograr que el
interruptor a tierra actue correctamente.
En la figura 15 se ilustra una falla bifsica entre las fases A y B en la barra N001
que ocurre a los 0,02 s del tiempo total de simulacin.
runATP archivo1.atp
runATP archivo2.atp
runATP archivo3.atp
del for *.dbg *.tmp *.lis *.bin
move *.pl4 archivosPL4/Falla_XX
Figura 16. Lista de nombres de los archivos con extensin .atp a ser copiados.
En esta etapa se genera un gran nmero de copias del cdigo fuente del mismo
sistema de distribucin elctrica. Estas copias son dirigidas a la carpeta Falla_XX
correspondiente. Para cada una de estas carpetas se generan archivos, cuya
cantidad depende del tipo de falla, nmero de nodos del sistema y total de
resistencias de falla.
pl42mat *.pl4
move *.mat archivosMAT/Falla_XX
Se seleccionan cada una de las opciones de la interfaz que son las siguientes:
Anallisis de
PROCESAMIENTO DE DATOS resultados
DESCRIPTORES
EN MATLAB
El segundo criterio se tom por falta de informacin con respecto a las cargas
particulares. Al no conocer el nmero de usuarios conectados a cada
transformador, a qu fases, ni el rgimen de carga de cada transformador, fue
necesario suponer la conexin de los transformadores a plena carga, resultando
esto en la concentracin de los usuarios y en la eliminacin consecuente de las
redes de baja tensin.
Nodo ATP Nodos Inicial - Final Tramo Long. (m) Transfo. Carga (kVA) Cable
N041-N043 170067 170071 A 251 ACSR 2/0 Desnudo
N043-N044 P PROTECCION
N044-N045 170071 5920581 A 111 ACSR 2/0 Desnudo
N045-N046 5920581 6850487 A 102 6850487 150 ACSR 2/0 Desnudo
N045-N047 5920581 170091 A 141 ACSR 2/0 Desnudo
N047-N048 170091 5920641 A 248 5920641 45 ACSR 2/0 Desnudo
N047-N049 170091 5920607 A 66 ACSR 2/0 Desnudo
N049-N050 5920607 5920611 A 63 5920611 10 ACSR 2/0 Desnudo
N049-N051 5920607 170105 A 38 170105 150 ACSR 2/0 Desnudo
N043-N052 170071 170110 A 96 ACSR 2/0 Desnudo
N052-N053 170110 7505094 A 66 7505094 45 ACSR 2/0 Desnudo
N052-N054 170110 170124 A 126 ACSR 2/0 Desnudo
N054-N055 P PROTECCION
N055-N056 170124 7556594 A 179 ACSR 1 Desnudo
N056-N057 7556594 7556316 A 31 7556316 45 ACSR 2/0 Desnudo
N056-N058 7556594 6850502 A 65 6850502 225 ACSR 1 Desnudo
N054-N059 170124 170158 A 180 ACSR 2/0 Desnudo
N059-N060 170158 170162 A 74 170162 113 ACSR 4/0 Desnudo
N059-N061 P PROTECCION
N061-N062 170158 5920655 A 9 ACSR 2/0 Desnudo
N062-N063 5920655 5920669 A 321 ACSR 2 Desnudo
N063-N064 5920669 5920672 A 11 ACSR 2/0 Desnudo
N064-N065 5920672 5920686 A 23 ACSR 2/0 Desnudo
N064-N066 5920672 5920690 A 235 ACSR 2 Desnudo
N066-N067 5920690 7559127 A 65 7559127 75 ACSR 2/0 Desnudo
N066-N068 5920690 5920701 A 58 5920701 (2) 150 ACSR 2/0 Desnudo
N068-N069 5920701 7389715 A 114 7389715 45 ACSR 2 Desnudo
N059-N070 170158 170177 A 74 ACSR 2/0 Desnudo
N070-N071 P PROTECCION
N071-N072 170177 170181 A 19 ACSR 2/0 Desnudo
N072-N073 170181 170215 A 1400 170215 30 ACSR 1 Desnudo
N073-N074 170215 170234 A 663 170234 30 ACSR 1 Desnudo
N074-N075 170234 5920781 A 81 ACSR 2/0 Desnudo
N075-N076 5920781 170249 A 249 170249 25 ACSR 1 Desnudo
N070-N077 170177 170272 A 460 170272 45 ACSR 2/0 Desnudo
N077-N078 170272 170287 A 77 ACSR 2/0 Desnudo
N078-N079 170287 5022664 A 89 Al 2/0 Desnudo
N079-N080 5022664 5022678 A 108 5022678 75 ACSR 2 Desnudo
N079-N081 5022664 170291 A 132 ACSR 2/0 Desnudo
N081-N082 P PROTECCION
N082-N083 170291 5920806 A 16 ACSR 2/0 Desnudo
N083-N084 5920806 5002540 A 132 5002540 150 Al 2/0 Desnudo
N084-N085 5002540 5920837 A 358 5920837 75 ACSR 2/0 Desnudo
N081-N086 170291 5022681 A 24 Al 2/0 Desnudo
N086-N087 5022681 5920841 A 109 ACSR 4/0 Desnudo
N087-N088 5920841 7360301 A 26 7360301 75 ACSR 2/0 Desnudo
N086-N089 5022681 6850581 A 102 ACSR 2/0 Desnudo
55
Nodo ATP Nodos Inicial - Final Tramo Long. (m) Transfo. Carga (kVA) Cable
N089-N090 P PROTECCION
N090-N091 6850581 7360346 A 103 7360346 75 Al 1 Desnudo
N091-N092 7360346 7360329 A 92 7360329 113 ACSR 4 Desnudo
N092-N093 7360329 6850564 A 610 6850564 10 ACSR 4 Desnudo
N089-N094 6850581 5920868 A 88 ACSR 2/0 Desnudo
N094-N095 5920868 5920885 A 156 ACSR 2/0 Desnudo
N095-N096 5920885 5788531 S 83 5788531 300 Cu 2/0 Triplex
N094-N097 5920868 170311 A 44 170311 113 ACSR 2/0 Desnudo
N097-N098 170311 170325 A 175 ACSR 2/0 Desnudo
N098-N099 170325 170330 S 14 Cu 2/0 Triplex
N099-N100 170330 5920899 A 69 ACSR 2/0 Desnudo
N100-N101 5920899 5920901 A 31 5920901 113 ACSR 2/0 Desnudo
N100-N102 5920899 5920914 A 218 ACSR 2/0 Desnudo
N102-N103 5920914 5920928 A 87 5920928 150 ACSR 2/0 Desnudo
N102-N104 5920914 170344 A 516 170344 10 ACSR 2/0 Desnudo
N104-N105 170344 170363 A 436 170363 15 ACSR 2 Desnudo
N105-N106 170363 170378 A 224 ACSR 2 Desnudo
N106-N107 170378 170382 A 689 170382 30 ACSR 4 Desnudo
N106-N108 170378 6825290 A 30 ACSR 2 Desnudo
N108-N109 P PROTECCION
N109-N110 6825290 170397 A 56 ACSR 2/0 Desnudo
N110-N111 P PROTECCION
N111-N112 170397 6954351 A 259 6954351 113 ACSR 2 Desnudo
N112-N113 6954351 170598 A 130 ACSR 2 Desnudo
N113-N114 170598 170602 A 328 170602 15 Al 4 Desnudo
N114-N115 170602 170617 A 131 Al 4 Desnudo
N115-N116 P PROTECCIN
N116-N117 170617 170621 A 66 ACSR 2/0 Desnudo
N117-N118 170621 170636 A 548 170636 30 Al 1 Desnudo
N118-N119 170636 494755 A 914 494755 15 Al 1 Desnudo
N115-N120 170617 170655 A 507 170655 15 Al 4 Desnudo
N120-N121 170655 170660 A 401 Al 4 Desnudo
N121-N122 170660 170674 A 72 170674 15 Al 1 Desnudo
N121-N123 170660 170689 A 99 Al 4 Desnudo
N123-N124 170689 170693 A 10 ACSR 2/0 Desnudo
N124-N125 170693 170708 A 835 170708 15 Al 1 Desnudo
N123-N126 170689 170712 A 577 170712 25 Al 4 Desnudo
N126-N127 170712 170731 A 892 170731 15 Al 4 Desnudo
N127-N128 170731 643892 A 201 Al 4 Desnudo
N128-N129 643892 5921022 A 748 5921022 15 Al 4 Desnudo
N128-N130 643892 170751 A 445 Al 4 Desnudo
N130-N131 P PROTECCION
N131-N132 170751 170770 A 677 170770 15 Al 4 Desnudo
N132-N133 170770 170784 A 664 170784 15 Al 4 Desnudo
N130-N134 170751 170799 A 369 170799 15 ACSR 2/0 Desnudo
N134-N135 170799 170803 A 270 170803 15 Al 4 Desnudo
N135-N136 170803 6825315 A 178 6825315 15 Al 4/0 Desnudo
56
Nodo ATP Nodos Inicial - Final Tramo Long. (m) Transfo. Carga (kVA) Cable
N135-N137 170803 170818 A 448 170818 15 Al 4 Desnudo
N137-N138 170818 170822 A 387 170822 75 Al 4 Desnudo
N138-N139 170822 5921040 A 352 Al 4 Desnudo
N139-N140 5921040 5921053 A 311 5921053 25 Al 4 Desnudo
N139-N141 5921040 5921084 A 1093 5921084 25 Al 4 Desnudo
N141-N142 5921084 5921098 A 190 ACSR 2/0 Desnudo
N142-N143 P PROTECCION
N143-N144 5921098 5921113 A 346 5921113 25 Al 4 Desnudo
N144-N145 5921113 7298125 A 440 7298125 15 Al 2/0 Desnudo
N142-N146 5921098 5921100 A 58 Al 4 Desnudo
N146-N147 P PROTECCION
N147-N148 5921100 5827660 A 797 5827660 15 ACSR 2/0 Desnudo
N148-N149 5827660 5827642 A 523 5827642 15 Al 4 Desnudo
N149-N150 5827642 5827625 A 804 5827625 25 Al 4 Desnudo
N150-N151 5827625 5827409 A 897 5827409 15 Al 4 Desnudo
N151-N152 5827409 5921161 A 434 ACSR 2/0 Desnudo
N152-N153 5921161 5827366 A 1090 5827366 10 Al 4 Desnudo
N152-N154 5921161 5827383 A 533 5827383 30 Al 4 Desnudo
N110-N155 170397 5920962 A 518 5920962 25 ACSR 2 Desnudo
N155-N156 5920962 170416 A 115 ACSR 2 Desnudo
N156-N157 170416 5920976 A 321 ACSR 2/0 Desnudo
N157-N158 5920976 5920980 A 386 ACSR 2 Desnudo
N158-N159 5920980 170435 A 210 170435 30 ACSR 2/0 Desnudo
N159-N160 P PROTECCION
N160-N161 170435 170545 A 279 ACSR 2/0 Desnudo
N161-N162 170545 170550 A 320 170550 10 ACSR 2 Desnudo
N162-N163 170550 170564 A 395 170564 25 ACSR 2/0 Desnudo
N163-N164 170564 170579 A 439 170579 25 ACSR 2 Desnudo
N159-N165 170435 643911 A 309 ACSR 2/0 Desnudo
N165-N166 643911 170440 A 427 ACSR 2 Desnudo
N166-N167 170440 170469 A 154 170469 75 ACSR 2/0 Desnudo
N167-N168 170469 170488 A 286 ACSR 2/0 Desnudo
N168-N169 170488 170492 A 433 170492 10 ACSR 2 Desnudo
N169-N170 P PROTECCION
N170-N171 170492 170856 A 179 ACSR 2/0 Desnudo
N171-N172 170856 170841 A 861 170841 15 ACSR 2 Desnudo
N172-N173 170841 5984317 A 433 ACSR 2/0 Desnudo
N173-N174 5984317 5984321 A 263 ACSR 2 Desnudo
N174-N175 5984321 5984334 A 202 ACSR 2/0 Desnudo
N175-N176 5984334 5984348 A 472 5984348 (2) 10 ACSR 2 Desnudo
N169-N177 170492 170507 A 497 ACSR 2/0 Desnudo
N177-N178 170507 170511 A 480 170511 10 ACSR 2/0 Desnudo
N177-N179 170507 6954316 A 293 6954316 15 ACSR 2/0 Desnudo
N179-N180 P PROTECCION
N180-N181 6954316 5921189 A 15 ACSR 2/0 Desnudo
N181-N182 5921189 5984274 A 409 ACSR 2 Desnudo
N182-N183 5984274 6954273 A 421 6954273 15 ACSR 2 Desnudo
57
Nodo ATP Nodos Inicial - Final Tramo Long. (m) Transfo. Carga (kVA) Cable
N183-N184 6954273 7296567 A 65 7296567 30 ACSR 2 Desnudo
N182-N185 5984274 5984288 A 456 5984288 15 ACSR 2 Desnudo
N185-N186 5984288 6954302 A 314 6954302 30 ACSR 2 Desnudo
N179-N187 6954316 6954975 A 928 6954975 15 ACSR 2 Desnudo
N187-N188 6954975 170861 A 743 170861 15 ACSR 2/0 Desnudo
N188-N189 170861 6954196 A 1561 6954196 15 ACSR 2 Desnudo
N189-N190 6954196 170913 A 865 170913 25 ACSR 2 Desnudo
N189-N191 6954196 170928 A 594 170928 75 ACSR 2 Desnudo
N191-N192 170928 170932 A 478 ACSR 2 Desnudo
N192-N193 P PROTECCION
N193-N194 170932 170947 A 585 170947 5 ACSR 2 Desnudo
N194-N195 170947 170951 A 536 170951 25 ACSR 2/0 Desnudo
N195-N196 170951 7466991 A 297 7466991 15 ACSR 2/0 Desnudo
N192-N197 170932 170966 A 280 ACSR 2 Desnudo
N197-N198 170966 170971 A 135 170971 30 ACSR 2/0 Desnudo
N198-N199 P PROTECCION
N199-N200 170971 171059 A 315 ACSR 2 Desnudo
N200-N201 171059 5984212 A 537 5984212 15 ACSR 2/0 Desnudo
N201-N202 5984212 171063 A 76 ACSR 2/0 Desnudo
N199-N203 170971 170985 A 529 ACSR 2 Desnudo
N203-N204 170985 116081 A 353 ACSR 2 Desnudo
N204-N205 116081 5800607 A 194 ACSR 2 Desnudo
N203-N206 170985 170990 A 290 170990 25 ACSR 2 Desnudo
N206-N207 170990 171006 A 103 ACSR 2/0 Desnudo
N207-N208 171006 171025 A 1226 ACSR 2 Desnudo
N208-N209 171025 171030 A 6 ACSR 2/0 Desnudo
N209-N210 171030 171044 A 454 171044 5 ACSR 2 Desnudo
Tabla 7. Listado descriptivo del circuito UM-33
Nodo ATP Nodos Inicial - Final Tramo Long. (m) Transfo. Carga (kVA) Cable
N017-N018 401822 401837 S 41 401837 400 Cu 2/0 Tripolar
N018-N019 P PROTECCION
N019-N020 401837 401861 S 116 401861 400 Cu 4/0 Triplex
N014-N015 401837 401803 S 112 Cu 4/0 Triplex
N014-N016 401803 401818 S 32 401818 30 Cu 2/0 Triplex
N018-N021 401837 5011295 S 225 Cu 2/0 Triplex
N021-N022 5011295 401909 S 20 401909 225 Cu 2/0 Tripolar
N022-N023 401909 403002 S 81 Cu 2/0 Tripolar
N023-N024 403002 476017 S 213 Cu 4/0 Triplex
N024-N025 476017 476002 S 89 Cu 2/0 Triplex
N025-N026 476002 401932 S 86 Cu 2/0 Tripolar
N026-N027 401932 401913 S 40 401913 500 Cu 4/0 Triplex
N027-N028 401913 401928 S 23 Cu 4/0 Triplex
N028-N029 401928 401947 S 50 Cu 2/0 Tripolar
N029-N030 P PROTECCION
N030-N031 401947 401951 S 12 401951 500 Cu 2/0 Tripolar
N030-N031 401947 401951 S 12 Cu 2/0 Tripolar
N029-N032 401947 475996 S 173 Cu 2/0 Tripolar
N032-N033 475996 28965 S 94 28965 300 Cu 4/0 Triplex
N033-N034 P PROTECCION
N034-N035 28965 401966 S 48 Cu 4/0 Triplex
N035-N036 401966 1479336 S 211 Cu 2/0 Tripolar
N036-N037 1479336 85819 S 38 85819 315 Cu 2/0 Triplex
N037-N038 85819 5011599 S 110 Cu 2/0 Tripolar
N038-N039 5011599 401990 S 7 Cu 4/0 Triplex
N039-N040 401990 402006 S 10 402006 225 Cu 2/0 Tripolar
N039-N040 401990 402006 S 10 Cu 4/0 Triplex
N039-N041 401990 402011 S 18 Cu 2/0 Tripolar
N033-N042 28965 28923 S 166 28923 160 Cu 2/0 Tripolar
N042-N043 P PROTECCION
N043-N044 28923 5089171 S 26 5089171 300 Cu 2/0 Triplex
N042-N045 28923 28939 S 5 Cu 2/0 Triplex
N045-N046 28939 28902 S 47 Cu 4/0 Triplex
N046-N047 28902 28853 S 66 Cu 4/0 Triplex
N047-N048 28853 28895 S 14 28895 300 Cu 2/0 Triplex
N048-N049 28895 28874 S 153 28874 400 Cu 2/0 Triplex
N049-N050 28874 6491933 S 17 Cu 2/0 Tripolar
N050-N051 6491933 28869 S 19 28869 475 Cu 4/0 Triplex
N051-N052 28869 5018046 S 70 Cu 4/0 Triplex
N052-N053 5018046 5018063 S 103 Cu 2/0 Triplex
N053-N054 5018063 28827 S 22 28827 600 Cu 2/0 Tripolar
N054-N055 28827 28757 S 11 Cu 2/0 Tripolar
N055-N056 28827 28806 S 109 Cu 2/0 Tripolar
N056-N057 28806 28811 S 23 28811 75 Cu 2/0 Triplex
N057-N058 P PROTECCION
N058-N059 28811 402996 S 165 402996 75 Cu 2/0 Triplex
N059-N060 402996 7561468 S 197 Cu 2/0 Triplex
59
Nodo ATP Nodos Inicial - Final Tramo Long. (m) Transfo. Carga (kVA) Cable
N059-N061 402996 7574795 A 82 7574795 75 ACSR 2 Desnudo
N059-N062 402996 6108198 S 284 6108198 50 Cu 2/0 Triplex
N062-N063 6108198 87471 S 211 Cu 4/0 Triplex
N057-N064 28811 28778 S 55 Cu 2/0 Triplex
N064-N065 28778 28783 S 30 Cu 2/0 Tripolar
N065-N066 28783 28799 S 14 28799 150 Cu 2/0 Monopolar
N066-N067 28799 5088068 S 119 5088068 225 Cu 2/0 Monopolar
N067-N068 5088068 5011875 S 76 Cu 2/0 Monopolar
N068-N069 5011875 28741 S 197 Cu 2/0 Tripolar
N069-N070 28741 28832 S 9 Cu 2/0 Triplex
N070-N071 P PROTECCION
N071-N072 28832 28721 S 200 28721 150 Cu 2/0 Tripolar
N072-N073 28721 28666 S 116 28666 300 Cu 2/0 Tripolar
N073-N074 28666 6873685 S 103 6873685 225 Cu 2/0 Triplex
N074-N075 6873685 28902 S 40 28902 300 Cu 2/0 Triplex
Tabla 8. Listado descriptivo del circuito SF-15
Otro aspecto que cabe resaltar con respecto a la caracterizacin de los circuitos,
es el de la ubicacin de las protecciones. Dado que esta informacin no fue
suministrada por la empresa, las nicas protecciones que se sabe que existen en
la localizacin considerada, son los reconectadores del circuito UM-33.
La ubicacin de los fusibles fue designada para el circuito UM-33 con el criterio de
aislar alimentadores secundarios, es decir, que se desprendan del principal y que
tengan dos o ms cargas conectadas.
controlar el factor de potencia de tal forma que para los suscriptores o usuarios
no residenciales, y de los residenciales conectados a un nivel de tensin superior
a uno (1), el factor de potencia deber ser igual o superior a 0,9. 24 . Por lo tanto se
puede afirmar que el factor de potencia de cada usuario ser de alrededor de 0,9,
por lo que se toma este valor para realizar los clculos subsecuentes.
24
Empresa Electrificadota de Santander S.A. E.S.P. Norma tcnica de ESSA.
61
7. RESULTADOS
El motivo de esta divisin radica en que los descriptores que se calculan, pueden
clasificarse en aquellos que necesitan de simulaciones con alta frecuencia de
muestreo y corta duracin en tiempo de la simulacin, y en aquellos que necesitan
larga duracin del tiempo de simulacin y baja frecuencia de muestreo.
cada nodo del circuito (del N001 al N074). Las resistencias de falla
utilizadas son las seis mencionadas en la seccin 4.5. (0,5 , 5 , 10 , 20
, 30 y 40 ).
7000
[A]
5250
3500
1750
-1750
-3500
-5250
-7000
0.00 0.04 0.08 0.12 0.16 [s] 0.20
(file CodensaSF15conProtecciones.pl4; x-var t) c:X0102A-X0025A c:X0102C-X0025C c:X0102B-X0025B
Figura 20. Corrientes de fase frente a falla trifsica a tierra en N073 de SF-15. (Fase A Rojo,
Fase B Verde, Fase C Azul).
Este comportamiento sucede dado que la falla ocurre aguas abajo de un fusible.
As que el tiempo que perdura la corriente de falla es consecuencia del tiempo que
demora el fusible en aislar la falla. La tensin se comporta como se ve en la figura
21.
65
20
[kV]
15
10
-5
-10
-15
-20
0.00 0.04 0.08 0.12 0.16 [s] 0.20
(file CodensaSF15conProtecciones.pl4; x-var t) v:X0025A v:X0025B v:X0025C
Figura 21. Tensiones de lnea ante una falla trifsica a tierra en N073 de SF-15.(Fase A Rojo,
Fase B Verde, Fase C Azul).
8000
[A]
6000
4000
2000
-2000
-4000
-6000
-8000
0.0 0.1 0.2 0.3 0.4 0.5 0.6 [s] 0.7
(file CodensaSF15conProtecciones.pl4; x-var t) c:X0102A-X0025A a.
20
[kV]
15
10
-5
-10
-15
-20
0.0 0.1 0.2 0.3 0.4 0.5 0.6 [s] 0.7
(file CodensaSF15conProtecciones.pl4; x-var t) v:X0025A b.
67
8000
[A]
6000
4000
2000
-2000
-4000
-6000
-8000
0.0 0.1 0.2 0.3 0.4 0.5 0.6 [s] 0.7
(file CodensaSF15conProtecciones.pl4; x-var t) c:X0102B-X0025B c.
25.00
[kV]
18.75
12.50
6.25
0.00
-6.25
-12.50
-18.75
-25.00
0.0 0.1 0.2 0.3 0.4 0.5 0.6 [s] 0.7
(file CodensaSF15conProtecciones.pl4; x-var t) v:X0025B d.
68
8000
[A]
6000
4000
2000
-2000
-4000
-6000
-8000
0.0 0.1 0.2 0.3 0.4 0.5 0.6 [s] 0.7
(file CodensaSF15conProtecciones.pl4; x-var t) c:X0102C-X0025C e.
25.00
[kV]
18.75
12.50
6.25
0.00
-6.25
-12.50
-18.75
-25.00
0.0 0.1 0.2 0.3 0.4 0.5 0.6 [s] 0.7
(file CodensaSF15conProtecciones.pl4; x-var t) v:X0025C f.
69
Figura 22. Seales de tensin y corriente al ocurrir una falla trifsica en N063
a. Corriente de fase A ante falla en N063, b. Tensin de lnea ante falla en N063, c. Corriente de
fase B ante falla en N063, d. Tensin de lnea ante falla en N063, e. Corriente de fase A ante falla
en N063, f. Tensin de lnea ante falla en N063
El comportamiento del resto del circuito frente a los diferentes tipos de falla con
sus resistencias de falla, vara en funcin de la magnitud de la corriente. Esta
variacin se presenta principalmente en el tiempo de duracin de la corriente de
falla, siendo este periodo corto cuando la magnitud de la corriente de falla es
grande y siendo este periodo largo cuando la magnitud de la corriente es pequea.
Ya que la corriente depende de la resistencia de falla y del tipo de falla, se podra
decir que el tiempo de accin de menor duracin ocurre cuando existe una falla
trifsica con resistencia de falla 0,5 en el nodo N030 y la mayor duracin en
caso de presentarse una falla bifsica con una resistencia de falla de 40 en el
nodo elctricamente ms alejado de la fuente.
a.
b.
71
c.
Figura 23. Ondas de tensin falla tipo 1 en N030 con resistencia de falla de 0.5
a. Onda de tensin de fase A, b. Onda de tensin de fase B, c. Onda de tensin de fase C. .
a.
b.
72
c.
Figura 23. Ondas de corriente falla tipo 1 en N030 con resistencia de falla de 0.5
a. Onda de corriente de fase A, b. Onda de corriente de fase B, c. Onda de tensin de fase C.
Se puede notar en esta figura que existen tres intervalos en los que la corriente de
falla circula por la red. El primer intervalo mostrado es desconectado por el
reconectador actuando con su curva rpida, esto ocurre tambin para el segundo
74
intervalo. El tercer intervalo se calcula con la curva lenta del reconectador, as que
al ocurrir la falla aguas abajo del fusible, este ltimo intervalo corresponde al
tiempo que demora el fusible en desconectar la falla. Si la falla acontece aguas
arriba del fusible, el comportamiento ser el mismo, con la diferencia de que el
tiempo del ltimo intervalo ser mayor. Tambin se puede ver en la figura 25, que
la tensin sufre una cada muy leve, casi inapreciable sin realizar un detalle de la
onda. Esto muestra la lejana de la fuente a la que ocurre la falla.
1600
[A]
1200
800
400
-400
-800
-1200
-1600
0.0 0.3 0.6 0.9 1.2 [s] 1.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) c:NF0A -NF01A a.
75
20
[kV]
15
10
-5
-10
-15
-20
0.0 0.3 0.6 0.9 1.2 [s] 1.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) v:NF01A b.
2000
[A]
1500
1000
500
-500
-1000
-1500
-2000
0.0 0.3 0.6 0.9 1.2 [s] 1.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) c:NF0B -NF01B c.
76
20
[kV]
15
10
-5
-10
-15
-20
0.0 0.3 0.6 0.9 1.2 [s] 1.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) v:NF01B d.
2000
[A]
1500
1000
500
-500
-1000
-1500
-2000
0.0 0.3 0.6 0.9 1.2 [s] 1.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) c:NF0C -NF01C e.
77
20
[kV]
15
10
-5
-10
-15
-20
0.0 0.3 0.6 0.9 1.2 [s] 1.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) v:NF01C f.
Figura 25. Seales de tensin y corriente al ocurrir una falla trifsica en N196
a. Corriente de fase A ante falla en N196, b. Tensin fase A ante falla en N196, c. Corriente de
fase B ante falla en N196, d. Tensin fase B ante falla en N196, e. Corriente de fase C ante falla en
N196, f. Tensin fase C ante falla en N196.
Las divisiones que se desprenden del alimentador principal desde el nodo N001 al
N108 se comportan de similar manera que el ejemplo mostrado por las figuras 20
y 21 para el circuito SF-15, ya que su esquema de proteccin es rel-fusible. Las
fallas que ocurren en este mismo intervalo de nodos, pero sobre el alimentador
principal, muestran el mismo comportamiento que el de la figura 22 del circuito SF-
15.
Las diferencias entre uno y otro circuito se aprecian al analizar tanto las
caractersticas de las cadas de tensin, como las de las seales de corriente
(diferencias de magnitud entre estado prefalla y postfalla, tiempo de actuacin de
las protecciones, etc.).
78
Se puede ver en la figura 26, las ondas de corriente y tensin al ocurrir una falla
bifsica de las fases B y C en el nodo N116, la cual presenta una resistencia de
falla de 0,5 .
700
[A]
525
350
175
-175
-350
-525
-700
0.0 0.1 0.2 0.3 0.4 [s] 0.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) c:NF0A -NF01A a.
79
20
[kV]
15
10
-5
-10
-15
-20
0.0 0.1 0.2 0.3 0.4 [s] 0.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) v:NF01A b.
3000
[A]
2000
1000
-1000
-2000
-3000
0.0 0.1 0.2 0.3 0.4 [s] 0.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) c:NF0B -NF01B c.
80
20
[kV]
15
10
-5
-10
-15
-20
0.0 0.1 0.2 0.3 0.4 [s] 0.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) v:NF01B d.
3000
[A]
2000
1000
-1000
-2000
-3000
0.0 0.1 0.2 0.3 0.4 [s] 0.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) c:NF0C -NF01C e.
81
20
[kV]
15
10
-5
-10
-15
-20
0.0 0.1 0.2 0.3 0.4 [s] 0.5
(file CodensaUM33conProteccionesUltimate.pl4; x-var t) v:NF01C f.
Figura 26. Ondas de corriente y tensin para falla bifsica fases B-C en N116
a. Corriente fase A, b. Tensin fase A, c. Corriente fase B, d. Tensin fase B, e. Corriente fase C,
f. Tensin fase C.
a.
b.
c.
Figura 27. Variacin de la tensin eficaz en cada una de las fases entre prefalla y falla, de un
evento tipo 3 en las barras N001 a N210 con 6 valores de resistencia de falla. a. Variacin de la
tensin eficaz en la fase A, b. Variacin de la tensin eficaz en la fase B, c. Variacin de la tensin
eficaz en la fase C.
84
a.
b.
c.
Figura 28. Variacin de la tensin eficaz de lnea entre prefalla y falla, de una falla tipo 3 en las
barras N001 a N210 con 6 valores de resistencia de falla. a. Variacin de la tensin eficaz entre las
lneas AB, b. Variacin de la tensin eficaz entre las lneas BC, c. Variacin de la tensin eficaz
entre las lneas CA.
85
a.
b.
c.
Figura 29. Variacin de la corriente eficaz en cada una de las fases entre prefalla y falla, de una
falla tipo 3 en las barras N001 a N210 con 6 valores de resistencia de falla. a. Variacin de la
corriente eficaz en la fase A, b. Variacin de corriente eficaz en la fse B, c. Variacin de la
corriente eficaz en la fase C.
86
a.
b.
c.
Figura 30. Variacin de la corriente eficaz de lnea entre prefalla y falla, de una falla tipo 3 en las
barras N001 a N210 con 6 valores de resistencia de falla. a. Variacin de la corriente eficaz entre
las lneas AB, b. Variacin de la corriente eficaz entre las lneas BC, c. Variacin de la corriente
eficaz entre las lneas CA.
87
a.
b.
c.
Figura 31. Variacin de la potencia aparente en cada una de las fases entre prefalla y falla, de una
falla tipo 3 en las barras N001 a N210 con 6 valores de resistencia de falla. a. Variacin de la
potencia aparente en la fase A, b. Variacin de la potencia aparente en la fase B, c. Variacin de la
potencia aparente en la fase C.
88
a.
b.
c.
Figura 32. Variacin de la potencia aparente de lnea entre prefalla y falla, de una falla tipo 3 en
las barras N001 a N210 con 6 valores de resistencia de falla a. Variacin de la potencia aparente
entre las lneas AB, b. Variacin de la potencia aparente entre las lneas BC, c. Variacin de la
potencia aparente entre las lneas CA.
89
a.
b.
c.
Figura 33. Reactancia de falla por fase, de una falla tipo 3 en las barras N001 a N210 con 6
valores de resistencia de falla. a. Reactancia de falla en la fase A, b. Reactancia de falla en la fase
B, c. Reactancia de falla en la fase C.
90
b.
c.
Figura 34. Reactancia de falla de lnea, de una falla tipo 3 en las barras N001 a N050 con 6
valores de resistencia de falla. a. Reactancia de falla entre las lneas AB, b. Reactancia de falla
entre las lneas BC, c. Reactancia de falla entre las lneas CA.
91
a.
b.
c.
Figura 35. Variacin del factor de potencia en cada una de las fases entre prefalla y falla, de una
falla tipo 3 en las barras N001 a N210 con 6 valores de resistencia de falla. a. Variacin del factor
de potencia en la fase A, b. Variacin del factor de potencia en la fase B, c. Variacin del factor de
potencia en la fase C.
92
a.
b.
c.
Figura 36. Variacin del factor de potencia de lnea entre prefalla y falla, de una falla tipo 3 en las
barras N001 a N210 con 6 valores de resistencia de falla a. Variacin del factor de potencia entre
las lneas AB, b. Variacin del factor de potencia entre las lneas BC, c. Variacin del factor de
potencia entre las lneas CA.
93
a.
b.
c.
Figura 37. Frecuencia del transitorio de falla, de una falla tipo 3 en las barras N001 a N210 con 6
valores de resistencia de falla. a. Frecuencia del transitorio de falla fase A, b. Frecuencia del
transitorio de falla fase B, c. Frecuencia del transitorio de falla fase C.
CONCLUSIONES
A partir del trabajo desarrollado durante este proyecto de grado, los autores han
considerado las siguientes como las conclusiones ms relevantes.
Para los circuitos SF-15 (Ricaurte) y UM-33 (Pasquilla), se encontr que los
descriptores de mejor comportamiento son aquellos relacionados a la
tensin, corriente, reactancia de falla y potencia aparente, ya que presentan
una caracteristica regular en funcin de la distancia elctrica desde la
fuente al punto de falla.
El lector debe conocer que este trabajo de grado es un punto intermedio entre la
teora del anlisis de la calidad de la energa elctrica y la prctica de este anlisis
a la localizacin de fallas en sistemas de distribucin. Existen diversas
extensiones y trabajos que pueden ser basados en lo ya realizado por los autores
de este proyecto. A continuacin algunas recomendaciones y trabajos futuros.
BIBLIOGRAFA
[9], [13] P. Hevia, Orlando. Comparacin de los modelos de lnea del ATP
[en lnea]. Argentina. 8 p. Disponible en web:
http://iitree.ing.unlp.edu.ar/estudios/Archivos/modlin.pdf
[10] Taku Noda, Akihiro Ametani. User Instructions of Noda Setup in ATP.
Julio 1998.
Gabor Furst, Models Primer, for first time MODELS users, Primera
Edicin.
Paulo F. Ribeiro. Wavelet transform and advanced tool for analyzing non-
stationary harmonic distortions in power system. In Proc of the IEEE
International Conference on Harmonics in Power Systems, Bologna, Italy,
September 1994.
ANEXOS
REL DE SOBRECORRIENTE
MODEL RELE
INPUT XXA
XXB
XXC
DATA NCLOS
TCLOS
IPK
A
B
P
DIAL
VAR ar[1..200], -- Vector que recoge las muestras de corrientes (FASE A).
br[1..200], -- Vector que recoge las muestras de corrientes (FASE B).
c[1..200], -- Vector que recoge las muestras de corrientes (FASE C).
YYA, -- Salida Booleana para el suiche que despejar la Falla (FASE A).
YYB, -- Salida Booleana para el suiche que despejar la Falla (FASE B).
YYC, -- Salida Booleana para el suiche que despejar la Falla (FASE C).
IRMSA, -- Valor RMS de la corriente (FASE A).
IRMSB, -- Valor RMS de la corriente (FASE B).
IRMSC, -- Valor RMS de la corriente (FASE C).
PP, KKA, JJA, -- Ancho de Ventana RMS; Para realizar corrimiento
KKB,JJB,KKC,JJC, -- Contadores para Inicio de la lectura de la Irms
AUXA, FLGDA, -- Variables para medir Irms (FASE A).
AUXB, FLGDB, -- Variables para medir Irms (FASE B).
AUXC, FLGDC, -- Variables para medir Irms (FASE C).
FLAG1, -- Bandera que indica primer deteccin de sobrecorriente FALLA L-T.(Clase 1)
FLAG2, -- Bandera que indica primer deteccin de sobrecorriente FALLA L-L-T y
otras.(Clase 2)
FLGRT, -- Bandera con el fin de dar tiempo para reconocer una falla L-L y no confundir con
L-T.
FLAGD1, -- Bandera para el disparo en el cruce por cero
FLGD2A, -- Bandera para el disparo en el cruce por cero Fase A
FLGD2B, -- Bandera para el disparo en el cruce por cero Fase B
FLGD2C, -- Bandera para el disparo en el cruce por cero Fase C
FLGCL, -- Bandera que indica que ya se calculo los tiempo de accion del Recloser
IITT1, -- Valor instantneo de la Corriente de Falla L-T cuando debe actuar la proteccin.
IT2A, -- Valor instantneo de la Corriente de Falla L-L FASE A cuando debe actuar la proteccin
IT2B, -- Valor instantneo de la Corriente de Falla L-L FASE B cuando debe actuar la proteccin
IT2C, -- Valor instantneo de la Corriente de Falla L-L FASE C cuando debe actuar la proteccin
CONR1, -- Contador Nmero de acciones Rpidas. Inicialmente se le asignar NCLOS
CONR2, -- Contador Nmero de acciones Rpidas. Inicialmente se le asignar NCLOS
TTF -- Instante de Tiempo en el que se detecta la falla.
MM, -- Mayor de los mtiplos
MMA, -- Mltiplo de la Ipk Fase A
MMB, -- Mltiplo de la Ipk Fase B
MMC, -- Mltiplo de la Ipk Fase C
TEMRA1,TEMRA2, -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 1 y 2.
TEMRC1,TEMRC2, -- Temporizador Tiempo de CIERRE, curva RAPIDA, Falla L-T (1) y Falla L-L (2)
TSETRA, -- T1:Tiempo RAPIDO para ABRIR proteccin. Calculado con la Curva del recloser
TSETRC, -- T2:Tiempo RAPIDO para CERRAR proteccin, equivale a
OUTPUT YYA, YYB, YYC
EXEC
comment INICIALIZACION DE VARIABLES ------- endcomment
IF T=0 THEN
ar[1..200]:=0 -- Vector que recoge las muestras de corrientes (FASE A).
br[1..200]:=0 -- Vector que recoge las muestras de corrientes (FASE B).
c[1..200]:=0 -- Vector que recoge las muestras de corrientes (FASE C).
YYA:=BOOL(1) -- Salida Booleana para el suiche que despejar la Falla (FASE A).
YYB:=BOOL(1) -- Salida Booleana para el suiche que despejar la Falla (FASE B).
YYC:=BOOL(1) -- Salida Booleana para el suiche que despejar la Falla (FASE C).
IRMSA:=0 -- Valor RMS de la corriente (FASE A).
N
TTF:=T
ENDIF
ENDIF
IF FLAG2=1 THEN
FLAG1:=0
ENDIF
comment *****************************************************
******* CALCULO DE LOS TIEMPOS DE ACCION DEL RELE **********
************************************************** endcomment
IF (FLAG1=1) OR (FLAG2=1)
AND FLGCL=0 THEN
IF (T>=TTF+1.5*PP) AND (T<=TTF+1.5*PP+2*timestep) THEN
IF IRMSA>=IPK THEN
MMA:=IRMSA/IPK
ENDIF
IF IRMSB>=IPK THEN
MMB:=IRMSB/IPK
ENDIF
IF IRMSC>=IPK THEN
MMC:=IRMSC/IPK
ENDIF
MM:=max(MMA,MMB,MMC)
TSETRA:=DIAL*(B+ A/(MM**P-1)) -- T1:Tiempo RAPIDO para ABRIR proteccin
FLGCL:=1
ENDIF
ENDIF
IF (TEMRA1>TSETRA-6*PP AND TEMRA1<TSETRA-6*PP+10*timestep) THEN
IF IRMSA>=IPK THEN
MMA:=IRMSA/IPK
ENDIF
IF IRMSB>=IPK THEN
MMB:=IRMSB/IPK
ENDIF
IF IRMSC>=IPK THEN
MMC:=IRMSC/IPK
ENDIF
MM:=max(MMA,MMB,MMC)
TSETRA:=DIAL*(B+ A/(MM**P-1)) -- T1:Tiempo RAPIDO para ABRIR proteccin
ENDIF
comment **************************************************
*********** FALLAS L - T *******************
*********************************************** endcomment
IF FLAG1=1 -- Falla L-T
AND FLGCL=1 -- Y ya se clcularon los Tiempos de apertura
AND CONR1>0 THEN -- Y Contador Fast N1 > 0
IF (IRMSA<IPK)
AND (IRMSB<IPK)
AND (IRMSC<IPK)
AND YYA=1 AND YYB=1 AND YYC=1 THEN
FLAG1:=0 -- Bandera que indica primer deteccin de sobrecorriente FALLA L-T.(Clase 1)
FLGRT:=0 -- Bandera para un instante en que FLAG1 se hace uno y no debera.
FLAGD1:=0 -- Bandera para el disparo en el cruce por cero
IITT1:=0 -- Valor instantneo de la Corriente de Falla L-T
TTF:=0 -- Instante de tiempo en el que se detecta la falla.
MM:=0 -- Mayor de los mtiplos
MMA:=0 -- Mltiplo de la Ipk Fase A
MMB:=0 -- Mltiplo de la Ipk Fase B
MMC:=0 -- Mltiplo de la Ipk Fase C
TEMRC1:=0 -- Temporizador Tiempo de CIERRE, curva RAPIDA, clases 1.
TEMRA1:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 1.
ENDIF
IF TEMRA1<=TSETRA THEN -- Si Temp.Rap to Abrir N1 < T1
TEMRA1:=TEMRA1+timestep -- Siga temporizando
ENDIF
IF TEMRA1>TSETRA THEN -- Si Temp Rap to Abrir N1> T1
IF IRMSA>=IPK THEN
IF FLAGD1=0 THEN
IITT1:=XXA
FLAGD1:=1
ENDIF
Q
IF FLGD2A=0 THEN
IT2A:=XXA
FLGD2A:=1
ENDIF
IF FLGD2A=1 AND XXA*IT2A<0 THEN
YYA:=BOOL(-1)
ENDIF
IF FLGD2B=0 THEN
IT2B:=XXB
FLGD2B:=1
ENDIF
IF FLGD2B=1 AND XXB*IT2B<0 THEN
YYB:=BOOL(-1)
ENDIF
IF FLGD2C=0 THEN
IT2C:=XXC
FLGD2C:=1
ENDIF
IF FLGD2C=1 AND XXC*IT2C<0 THEN
YYC:=BOOL(-1)
ENDIF
ENDIF
IF (TEMRC2<=TSETRC) -- Si Temp Rap to Cerrar N2 < T2
AND (YYA=0 AND YYB=0 AND YYC=0) -- Y las tres protecciones estan abiertas
THEN
TEMRC2:=TEMRC2+timestep -- Siga temporizando
ENDIF
IF TEMRC2>TSETRC -- Si Temp Rap to Cerrar N1 > T2
AND (YYA=0 AND YYB=0 AND YYC=0) THEN -- Y la proteccion esta abierta
CONR2:=CONR2-1 -- Decrementar el contador de Rapidos
IF CONR2>0 THEN
YYA:=BOOL(1) -- Terminar la Temp y Cierra.
YYB:=BOOL(1)
YYC:=BOOL(1)
ENDIF
TEMRA2:=0 -- Reset Temp Rap to Abrir
TEMRC2:=0 -- Reset Temp Rap to Cerrar
FLGD2A:=0
FLGD2B:=0
FLGD2C:=0
IT2A:=0
IT2B:=0
IT2C:=0
ENDIF
ENDIF
ENDEXEC
ENDMODEL
RECONECTADOR TRIFSICO
MODEL REC3F
INPUT XXA
XXB
XXC
DATA NFAST
NDEL
TFAST
TDEL
IPK
Af
Bf
Pf
DIALf
Ad
Bd
Pd
DIALd
VAR a[1..200], -- Vector que recoge las muestras de corrientes (FASE A).
b[1..200], -- Vector que recoge las muestras de corrientes (FASE B).
S
IF (T>0)AND(T<=PP) THEN
b[JJB]:=XXB**2*timestep
AUXB:=AUXB+b[JJB]
JJB:=JJB+1
ENDIF
IF T>PP THEN
IF FLGDB=0 THEN
IRMSB:=sqrt(AUXB/PP)
FLGDB:=1
ENDIF
IF FLGDB=1 THEN
IRMSB:=( ( IRMSB**2*PP ) + timestep*XXB**2 - b[1] ) / PP
IF IRMSB<=0 THEN
IRMSB:=0
ENDIF
IRMSB:=SQRT(IRMSB)
b[JJB-1]:=XXB**2*timestep
FOR KKB:=1 to JJB+2 DO
b[KKB]:=b[KKB+1]
ENDFOR
ENDIF
ENDIF
comment ----------------------------------------------------------
------- CORRIENTE RMS (FASE C) ----------
------------------------------------------------------- endcomment
IF (T>0)AND(T<=PP) THEN
c[JJC]:=XXC**2*timestep
AUXC:=AUXC+c[JJC]
JJC:=JJC+1
ENDIF
IF T>PP THEN
IF FLGDC=0 THEN
IRMSC:=sqrt(AUXC/PP)
FLGDC:=1
ENDIF
IF FLGDC=1 THEN
IRMSC:= ( ( IRMSC**2*PP ) + timestep*XXC**2 - c[1] ) / PP
IF IRMSC<=0 THEN
IRMSC:=0
ENDIF
IRMSC:=SQRT(IRMSC)
c[JJC-1]:=XXC**2*timestep
FOR KKC:=1 to JJC+2 DO
c[KKC]:=c[KKC+1]
ENDFOR
ENDIF
ENDIF
comment *****************************************************
********* INSTANTE EN EL QUE SE DETECTA LA FALLA ***********
************************************************** endcomment
IF FLAG1=0
AND IRMSA>=IPK AND IRMSB<=IPK AND IRMSC<=IPK
OR IRMSB>=IPK AND IRMSA<=IPK AND IRMSC<=IPK
OR IRMSC>=IPK AND IRMSA<=IPK AND IRMSB<=IPK THEN
IF FLGRT=0 THEN
TTF:=T
FLGRT:=1
ENDIF
IF (FLGRT=1) AND (T>=TTF+2*PP/3) THEN
FLAG1:=1
ENDIF
ENDIF
IF FLAG2=0
AND IRMSA>=IPK AND IRMSB>=IPK
OR IRMSB>=IPK AND IRMSC>=IPK
OR IRMSC>=IPK AND IRMSA>=IPK THEN
FLAG2:=1
IF TTF=0 THEN
TTF:=T
ENDIF
V
ENDIF
IF FLAG2=1 THEN
FLAG1:=0
ENDIF
comment *****************************************************
******* CALCULO DE LOS TIEMPOS DE ACCION DEL RELE **********
************************************************** endcomment
IF (FLAG1=1) OR (FLAG2=1)
AND FLGCL=0 THEN
IF (T>=TTF+1.5*PP) AND (T<=TTF+1.5*PP+2*timestep) THEN
IF IRMSA>=IPK THEN
MMA:=IRMSA/IPK
ENDIF
IF IRMSB>=IPK THEN
MMB:=IRMSB/IPK
ENDIF
IF IRMSC>=IPK THEN
MMC:=IRMSC/IPK
ENDIF
MM:=max(MMA,MMB,MMC)
TSETRA:=DIALf*(Bf+ Af/(MM**Pf-1)) -- T1:Tiempo RAPIDO para ABRIR proteccin
TSETLA:=DIALd*(Bd+ Ad/(MM**Pd-1)) -- T3:Tiempo LENTO para ABRIR proteccin
FLGCL:=1
ENDIF
IF (TEMRA1>TSETRA-6*PP AND TEMRA1<TSETRA-6*PP+10*timestep)
OR (TEMLA1>TSETLA-6*PP AND TEMLA1<TSETLA-6*PP+10*timestep) THEN
IF IRMSA>=IPK THEN
MMA:=IRMSA/IPK
ENDIF
IF IRMSB>=IPK THEN
MMB:=IRMSB/IPK
ENDIF
IF IRMSC>=IPK THEN
MMC:=IRMSC/IPK
ENDIF
MM:=max(MMA,MMB,MMC)
TSETRA:=DIALf*(Bf+ Af/(MM**Pf-1)) -- T1:Tiempo RAPIDO para ABRIR proteccin
TSETLA:=DIALd*(Bd+ Ad/(MM**Pd-1)) -- T3:Tiempo LENTO para ABRIR proteccin
ENDIF
IF (TEMRA2>TSETRA-6*PP AND TEMRA2<TSETRA-6*PP+10*timestep)
OR (TEMLA2>TSETLA-6*PP AND TEMLA2<TSETLA-6*PP+10*timestep) THEN
IF IRMSA>=IPK THEN
MMA:=IRMSA/IPK
ENDIF
IF IRMSB>=IPK THEN
MMB:=IRMSB/IPK
ENDIF
IF IRMSC>=IPK THEN
MMC:=IRMSC/IPK
ENDIF
MM:=max(MMA,MMB,MMC)
TSETRA:=DIALf*(Bf+ Af/(MM**Pf-1)) -- T1:Tiempo RAPIDO para ABRIR proteccin
TSETLA:=DIALd*(Bd+ Ad/(MM**Pd-1)) -- T3:Tiempo LENTO para ABRIR proteccin
ENDIF
ENDIF
comment **************************************************
*********** FALLAS L - T *******************
*********************************************** endcomment
IF FLAG1=1 -- Falla L-T
AND FLGCL=1 -- Y ya se clcularon los Tiempos de apertura
AND CONR1>0 THEN -- Y Contador Fast N1 > 0
IF (IRMSA<IPK)
AND (IRMSB<IPK)
AND (IRMSC<IPK)
AND YYA=1 AND YYB=1 AND YYC=1 THEN
FLAG1:=0 -- Bandera que indica primer deteccin de sobrecorriente FALLA L-T.(Clase 1)
FLGRT:=0 -- Bandera para un instante en que FLAG1 se hace uno y no debera.
FLAGD1:=0 -- Bandera para el disparo en el cruce por cero
IITT1:=0 -- Valor instantneo de la Corriente de Falla L-T
TTF:=0 -- Instante de tiempo en el que se detecta la falla.
W
FLAG2:=0 -- Bandera que indica primer deteccin de sobrecorriente FALLA L-L-T y otras.(Clase 2)
FLGRT:=0 -- Bandera para un instante en que FLAG1 se hace uno y no debera.
FLGD2A:=0 -- Bandera para el disparo en el cruce por cero Fase A
FLGD2B:=0 -- Bandera para el disparo en el cruce por cero Fase B
FLGD2C:=0 -- Bandera para el disparo en el cruce por cero Fase C
IT2A:=0 -- Valor instantneo de la Corriente de Falla L-L FASE A
IT2B:=0 -- Valor instantneo de la Corriente de Falla L-L FASE B
IT2C:=0 -- Valor instantneo de la Corriente de Falla L-L FASE
TTF:=0 -- Instante de tiempo en el que se detecta la falla.
MM:=0 -- Mayor de los mtiplos
MMA:=0 -- Mltiplo de la Ipk Fase A
MMB:=0 -- Mltiplo de la Ipk Fase B
MMC:=0 -- Mltiplo de la Ipk Fase C
TEMRC2:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 2.
TEMLC2:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 2.
TEMRA2:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 2.
TEMLA2:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 2.
ENDIF
IF TEMRA2<=TSETRA THEN -- Si Temp.Rap to Abrir N2 < T1
TEMRA2:=TEMRA2+timestep -- Siga temporizando
ENDIF
IF TEMRA2>TSETRA THEN -- Si Temp Rap to Abrir N2 > T1
IF FLGD2A=0 THEN
IT2A:=XXA
FLGD2A:=1
ENDIF
IF FLGD2A=1 AND XXA*IT2A<0 THEN
YYA:=BOOL(-1)
ENDIF
IF FLGD2B=0 THEN
IT2B:=XXB
FLGD2B:=1
ENDIF
IF FLGD2B=1 AND XXB*IT2B<0 THEN
YYB:=BOOL(-1)
ENDIF
IF FLGD2C=0 THEN
IT2C:=XXC
FLGD2C:=1
ENDIF
IF FLGD2C=1 AND XXC*IT2C<0 THEN
YYC:=BOOL(-1)
ENDIF
ENDIF
IF (TEMRC2<=TSETRC) -- Si Temp Rap to Cerrar N2 < T2
AND (YYA=0 AND YYB=0 AND YYC=0) -- Y las tres protecciones estan abiertas
THEN
TEMRC2:=TEMRC2+timestep -- Siga temporizando
ENDIF
IF TEMRC2>TSETRC -- Si Temp Rap to Cerrar N1 > T2
AND (YYA=0 AND YYB=0 AND YYC=0) THEN -- Y la proteccion esta abierta
CONR2:=CONR2-1 -- Decrementar el contador de Rapidos
YYA:=BOOL(1) -- Terminar la Temp y Cierra.
YYB:=BOOL(1)
YYC:=BOOL(1)
TEMRA2:=0 -- Reset Temp Rap to Abrir
TEMRC2:=0 -- Reset Temp Rap to Cerrar
FLGD2A:=0
FLGD2B:=0
FLGD2C:=0
IT2A:=0
IT2B:=0
IT2C:=0
ENDIF
ENDIF
comment ---------- DERROTERO ACTUACION LENTA --------------- endcomment
IF CONR2=0 -- Falla L - L - T
AND FLGCL=1 -- Y ya se clcularon los Tiempos de apertura
AND CONL2>0 THEN -- Y Contador Fast N2 > 0
IF (IRMSA<IPK)
Z
AND (IRMSB<IPK)
AND (IRMSC<IPK)
AND YYA=1 AND YYB=1 AND YYC=1 THEN
FLAG2:=0 -- Bandera que indica primer deteccin de sobrecorriente FALLA L-L-T y otras.(Clase 2)
FLGRT:=0 -- Bandera para un instante en que FLAG1 se hace uno y no debera.
FLGD2A:=0 -- Bandera para el disparo en el cruce por cero Fase A
FLGD2B:=0 -- Bandera para el disparo en el cruce por cero Fase B
FLGD2C:=0 -- Bandera para el disparo en el cruce por cero Fase C
IT2A:=0 -- Valor instantneo de la Corriente de Falla L-L FASE A
IT2B:=0 -- Valor instantneo de la Corriente de Falla L-L FASE B
IT2C:=0 -- Valor instantneo de la Corriente de Falla L-L FASE
TTF:=0 -- Instante de tiempo en el que se detecta la falla.
MM:=0 -- Mayor de los mtiplos
MMA:=0 -- Mltiplo de la Ipk Fase A
MMB:=0 -- Mltiplo de la Ipk Fase B
MMC:=0 -- Mltiplo de la Ipk Fase C
TEMRC2:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 2.
TEMLC2:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 2.
TEMRA2:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 2.
TEMLA2:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA, clases 2.
ENDIF
IF TEMLA2<=TSETLA THEN -- Si Temp.Rap to Abrir N2 < T1
TEMLA2:=TEMLA2+timestep -- Siga temporizando
ENDIF
IF TEMLA2>TSETLA THEN -- Si Temp Rap to Abrir N2 > T1
IF FLGD2A=0 THEN
IT2A:=XXA
FLGD2A:=1
ENDIF
IF FLGD2A=1 AND XXA*IT2A<0 THEN
YYA:=BOOL(-1)
ENDIF
IF FLGD2B=0 THEN
IT2B:=XXB
FLGD2B:=1
ENDIF
IF FLGD2B=1 AND XXB*IT2B<0 THEN
YYB:=BOOL(-1)
ENDIF
IF FLGD2C=0 THEN
IT2C:=XXC
FLGD2C:=1
ENDIF
IF FLGD2C=1 AND XXC*IT2C<0 THEN
YYC:=BOOL(-1)
ENDIF
ENDIF
IF (TEMLC2<=TSETLC) -- Si Temp Rap to Cerrar N2 < T2
AND (YYA=0 AND YYB=0 AND YYC=0) -- Y las tres protecciones estan abiertas
THEN
TEMLC2:=TEMLC2+timestep -- Siga temporizando
ENDIF
IF TEMLC2>TSETLC -- Si Temp Rap to Cerrar N1 > T2
AND (YYA=0 AND YYB=0 AND YYC=0) THEN -- Y la proteccion esta abierta
CONL2:=CONL2-1 -- Decrementar el contador de Rapidos
IF CONL2>0 THEN
YYA:=BOOL(1) -- Terminar la Temp y Cierra si no se ha terminado .
YYB:=BOOL(1) -- la operacin Lenta.
YYC:=BOOL(1)
ENDIF
TEMLA2:=0 -- Reset Temp Rap to Abrir
TEMLC2:=0 -- Reset Temp Rap to Cerrar
FLGD2A:=0
FLGD2B:=0
FLGD2C:=0
IT2A:=0
IT2B:=0
IT2C:=0
ENDIF
ENDIF
AA
ENDEXEC
ENDMODEL
FUSIBLE
MODEL FUSE
INPUT XXIN
DATA IPK
A
B
P
DIAL
VAR af[1..200], -- Vector que recoge las muestras de corrientes.
YY, -- Salida Booleana para el suiche que despejar la Falla.
TT, -- Instante de tiempo en el que se Detecta la falla.
IRMS, -- Valor RMS de la corriente.
PP, KK, JJ, -- Ancho de Ventana RMS; Para realizar corrimiento ;Contador para Inicio Irms
AUX, FLGD, -- Variables para medir Irms
FLAG, -- Bandera que indica primer deteccin de sobrecorriente.
TEMRA, -- Temporizador Tiempo de APERTURA.
TSET, -- T1:Tiempo para ABRIR proteccin
MM, -- Multiplo de la corriente de Arranque.
FLGCL, -- Bandera que indica que ya se calculo los tiempo de accion del Fusible
FLAGD1, -- Bandera para apertura en el cruce por cero
IT1 -- Valor de la Corriente en el momento de la orden de apertura
ddd
OUTPUT YY
EXEC
comment ***********************************************************
******** INICIALIZACION DE VARIABLES *******************
******************************************************** endcomment
IF T=0 THEN
af[1..200]:=0 -- Vector que recoge las muestras de corrientes.
YY:=BOOL(1) -- Salida Booleana para el suiche que despejar la Falla.
IRMS:=0 -- Valor RMS de la corriente.
PP:=1/120 -- Ancho de Ventana RMS
KK:=1 -- Para hacer corrimiento de Componentes del vector "a".
JJ:=2 -- Contador para inicio de la IRMS
AUX:=0 -- Variable auxiliar para medir la IRMS
FLGD:=0 -- Bandera Auxiliar para medir la IRMS
FLAG:=0 -- Bandera que indica primer deteccin de sobrecorriente.
TEMRA:= 0 -- Temporizador Tiempo de APERTURA, curva RAPIDA
TSET:=1 -- Tiempo de apertura. Se calcula de la Curva del Fuse.
FLAGD1:=0 -- Bandera para el disparo en el cruce por cero
IT1:=0 -- Valor de la Corriente en el instante de la orden de disparo
MM:=1 -- Mltiplo de la corriente de Arranque
TT:=0 -- Instante de tiempo en el que se detecta la falla.
FLGCL:=0 -- Bandera que indica que ya se calculo los tiempo de accion del Fusible
ddd:=0
ENDIF
comment ***********************************************************
********* LECTURA DE LA CORRIENTE RMS *******************
******************************************************** endcomment
IF (T>0)AND(T<=PP) THEN
af[JJ]:=XXIN**2*timestep
AUX:=AUX+af[JJ]
JJ:=JJ+1
ENDIF
IF T>PP THEN
IF FLGD=0 THEN
IRMS:=sqrt(AUX/PP)
FLGD:=1
ENDIF
IF FLGD=1 THEN
IRMS:= ( ( IRMS**2*PP ) + timestep*XXIN**2 - af[1] ) / PP
IF IRMS<=0 THEN
IRMS:=0
ENDIF
IRMS:=SQRT(IRMS)
BB
af[JJ-1]:=XXIN**2*timestep
FOR KK:=1 to JJ+2 DO
af[KK]:=af[KK+1]
ENDFOR
ENDIF
ENDIF
comment ***********************************************************
********* PRIMER INSTATANTE DE LA FALLA *******************
******************************************************** endcomment
IF FLAG=0 AND IRMS>=IPK THEN
FLAG:=1
TT:=T
ENDIF
IF FLAG=1 AND IRMS<IPK THEN
FLAG:=0
TT:=0
ENDIF
IF FLAG=1 AND FLGCL=0 THEN
IF T>=TT+2*PP AND T<=TT+2*PP+2*timestep THEN
MM:=IRMS/IPK
TSET:=DIAL*(B+ A/(MM**P-1))
FLGCL:=1
ENDIF
ENDIF
comment ***********************************************************
********* ACTUACION DEL FUSIBLE *******************
******************************************************** endcomment
IF IRMS<IPK AND YY=1 THEN
FLAG:=0 -- Bandera que indica primer deteccin de sobrecorriente.
TEMRA:=0 -- Temporizador Tiempo de APERTURA, curva RAPIDA
FLAGD1:=0 -- Bandera para el disparo en el cruce por cero
IT1:=0 -- Valor de la Corriente en el instante de la orden de disparo
MM:=1 -- Mltiplo de la corriente de Arranque
TT:=0 -- Instante de tiempo en el que se detecta la falla.
ENDIF
IF FLAG=1 -- Inicio de la Accion del Rel
AND FLGCL=1 THEN -- Y ya se calcul el tiempo de accin del Fuse
IF IRMS>=IPK -- Si hay Sobrecorriente
AND DDD=0
AND TEMRA<=TSET THEN -- Si Temp.Rap to Abrir < T1
YY:=BOOL(1) -- Siga temporizando
TEMRA:=TEMRA+timestep
ENDIF
IF TEMRA>TSET THEN -- Si Temp Rap to Abrir > T1
IF FLAGD1=0 THEN
IT1:=XXIN
FLAGD1:=1
ENDIF
IF FLAGD1=1 AND XXIN*IT1<0 AND DDD=0 THEN
YY:=BOOL(-1)
DDD:=1
ENDIF
ENDIF
ENDIF
ENDEXEC
ENDMODEL
CC