Revista Web
Revista Web
Enero 2018
Edición No. 14
(ISSN 1390-5074)
Filtro activo de potencia para instalaciones bifásicas con capacidad de balance de corrientes
64
M. Fajardo, J. Viola, F. Quizhpi
Metodología para Identificar Áreas de Control de Voltaje en un sistema eléctrico de potencia aplicando simulación Monte Carlo
72
L. Bonilla, D. Echeverría, J. Cepeda
Consideraciones Conceptuales sobre Análisis de Estabilidad Transitoria en el Ajuste del Tiempo de Recierre Monofásico de la Línea de Transmisión
San Rafael - El Inga de 500 kV 80
L. Dután, J. Cepeda.
EFICIENCIA ENERGÉTICA
Estudios eléctricos de la interconexión temporal del Sistema Eléctrico Interconectado Petrolero Extendido (SEIP-E) de Petroamazonas EP al Sistema
Nacional Interconectado (S.N.I.) 92
R. Sanguil, M. Solís
Análisis Espacio Temporal del Impacto al Sistema de Distribución de Nuevos Tipos de Carga con Alto Consumo de Energía: Caso de Estudio en
CENTROSUR con la Cocina de Inducción 101
S. Zambrano, B. Jaramillo, J. Cárdenas, M. Mejía, A. Padilha, J. Melo
Operación Óptima del Sistema Eléctrico Aislado de la Isla San Cristóbal, con Generación Renovable y Vehículos Eléctricos
110
R. Barreto, PF. O. Gomis J. Heredia
TÉCNOLÓGICOS E INNOVACIÓN
Evaluación de la Distorsión Armónica debida a la Operación del Compensador Estático de Potencia Reactiva en la Subestación Chorrillos
122
B. Jadán, F. Ordoñez, W. Vargas, H. Arcos
Comité Editorial
La Revista Técnica “energía” de CENACE, ha sido tradicionalmente publicada con una periodicidad anual (en el 2018 se
ha planificado modificar la periodicidad a semestral), y su distribución no tiene costo alguno. Se permite la reproducción
total o parcial de esta revista citando la fuente.
Indización
Revista Técnica “energía” está indizada en las bases de datos de “EBSCO Information Services”. El acceso a la versión
electrónica se encuentra en: [Link]
Revista Técnica “energía” está indizada en el catálogo Latindex, Sistema Regional de Información en Línea para Revistas
Científicas de América Latina, el Caribe, España y Portugal. El acceso a la catalogación electrónica se encuentra en
[Link]
MISIÓN DE LA REVISTA TÉCNICA “energía”
La Revista Técnica “energía”, auspiciada por el Operador Nacional de Electricidad, CENACE, está orientada a la
publicación del conocimiento científico en el país, permitiendo la difusión del conocimiento generado en las diferentes áreas
de CENACE, así como en otras instituciones relacionadas con los sectores eléctrico y energético, colegios profesionales,
institutos de investigación e instituciones de educación superior, cuyos temas sean de interés dentro del ámbito eléctrico,
energético y empresarial.
Dentro de la Revista Técnica “energía” existen dos categorías y cada autor o autores deben indicar en la presentación de
su trabajo técnico, en qué categoría desean participar:
CATEGORÍA I
Comprende los trabajos relacionados con aplicaciones prácticas o de desarrollo basados en la ejecución de los procesos de
las diferentes instituciones del sector eléctrico y energético.
CATEGORÍA II
Comprender los artículos basados en trabajos científicos y académicos de investigación de Universidades e Institutos, así
como también trabajos académicos de tesis de grado para la obtención de títulos de tercer nivel o cuarto nivel (especialidad,
maestría o doctorado).
El temario referencial de la Revista Técnica “energía” es el siguiente:
EFICIENCIA ENERGÉTICA
• Calidad de energía
• Energías renovables: regulación, comercialización e incorporación al sistema eléctrico
• Uso eficiente de la energía eléctrica y control de la demanda
GESTIÓN EMPRESARIAL
• Gestión de riesgos
• Gestión de recursos humanos
• Sistemas de gestión integrados
• Responsabilidad social
• Gestión del conocimiento
EDITORIAL
• Sistemas de Almacenamiento y
Generación Eficientes,
• Prácticas Socio-ambientales en la
gestión energética.
Aplicación Práctica / Practical Issues
Abstract Resumen
This article proposes a Spatial Load Forecasting Este artículo propone un modelo de Proyección
(SLF) model applied to an electric distribution Espacial de la Demanda (PED) aplicado a un
system for the short, medium and long term sistema eléctrico de distribución para el corto,
horizon, which includes the estimation of the mediano y largo plazo, lo cual incluye la estimación
magnitude and geographic location of electricity de la magnitud y la ubicación geográfica de la
demand energy from new consumers. demanda de energía eléctrica de nuevos clientes.
Index terms− Spatial Load Forecasting, Spatial Palabras clave− Proyección Espacial de la
Analysis, Geographic Information Systems, Land- Demanda, Análisis Espacial, Sistemas de
Use based Simulation, Cellular Automata. Información Geográfico, Simulación del uso del
suelo, Autómata celular.
7
Edición No 14, Enero 2018
8
Zambrano et al. / Modelo de Simulación Jerárquico para la Proyección Espacio Temporal de la Demanda Eléctrica
9
Edición No 14, Enero 2018
3.1. Creación del Mapa de Densidad de La zona de estudio se dividió en dos partes: el área
Demanda del Año Base urbana de Cuenca y otra para el resto del área de
concesión [1]. La Tabla 1 detalla los factores para esta
Se estable un año base en el cual se determina segunda zona. Mientras que el número de factores
la ubicación geográfica de los clientes y luminarias. considerados para el área urbana de Cuenca es diez,
Para este estudio, el año base corresponde al 2015, ya que se cuenta con información más detallada
con ello se crea el mapa de densidad de demanda para y completa obtenida del Plan de Desarrollo y
cuatro grupos de consumo: Residencial, Comercial, Ordenamiento Territorial de Cuenca (PDOT).
Industrial y Otros.
Tabla 1: Factores Espaciales para el Análisis del Área de
Mediante estudios de investigación y Concesión sin la parte Urbana de Cuenca
caracterización de la carga [2], se determinaron las
principales variables de interés, tales como el factor Ítem Descripción
de carga, el factor de demanda y factor de conversión
Factores de proximidad
del consumo mensual de energía a un valor de
potencia máxima para cada tipo o clase de cliente; F1 Proximidad a cantones
además se obtiene los perfiles de carga unitarios. F2 Proximidad a autopistas
F3 Proximidad a calles
A partir del valor de demanda máxima del sistema,
F4 Proximidad a centros educativos
en el mes y hora del año base, se distribuye dicho
valor en función de la potencia medida en clientes con F5 Proximidad a servicios de salud
demanda, potencia estimada en la gran mayoría de F6 Proximidad a parroquias
clientes (conversión del consumo de energía mensual) Factores locales
y potencia de luminarias [2]. El resultado es un mapa F7 Pendiente de la superficie
temático, Fig. 4, que muestra la densidad de demanda
F8 Suelo restringido
(kW/km2), en este caso de los clientes residenciales.
Las zonas de color rojo representan las áreas de alta Factores de entorno
densidad de carga correspondiente a zonas urbanas de F9 Densidad de demanda en las inmediaciones
las ciudades, en cambio las zonas de color amarillo
representan las áreas de baja densidad de carga. 3.2.1 Clasificación
Influencia 0 a 10 m = valor 1
10 a 28 m = valor 2
28 a 54 m = valor 3
54 a 89 m = valor 4
89 a 140 m = valor 5
Figura 4: Mapa de Demanda del Año Base para Clientes
Residenciales
Distancia a la carretera
3.2. Análisis de los Factores Espaciales y (metros)
10
Zambrano et al. / Modelo de Simulación Jerárquico para la Proyección Espacio Temporal de la Demanda Eléctrica
11
Edición No 14, Enero 2018
(3)
12
Zambrano et al. / Modelo de Simulación Jerárquico para la Proyección Espacio Temporal de la Demanda Eléctrica
Existen varias técnicas de asignación, entre Figura 12: Resultado de la Asignación de Nuevos Clientes en el
ellas: Montecarlo [3], Autómata celular (AC) [4] Cantón Gualaceo perteneciente a la SE #15
y Multiagente [10]. En este trabajo se utilizó AC,
debido al menor error obtenido al tomar una muestra
de análisis.
13
Edición No 14, Enero 2018
(10)
4. RESULTADOS
14
Zambrano et al. / Modelo de Simulación Jerárquico para la Proyección Espacio Temporal de la Demanda Eléctrica
15
Edición No 14, Enero 2018
16
Artículo Académico / Academic Paper
Abstract Resumen
The application of the phasor technology offers La aplicación de la tecnología fasorial ofrece nuevas
new possibilities for monitoring, protection, posibilidades para la supervisión, protección,
analysis and control of Electric Power Systems análisis y control de Sistemas Eléctricos de
(EPS). In this project, an algorithm is implemented Potencia (SEP). En este proyecto se implementa
to localize a failure using phasor data, which un algoritmo para la localización de fallas usando
improves traditional methods. The implemented datos fasoriales, el cual permitiría mejorar los
algorithm distinguishes responses of several métodos tradicionales. El algoritmo implementado
linearized models against measurements mode distingue las respuestas de varios modelos
during the transient response after the failure. linealizados frente a las mediciones fasoriales ante
Based on this, a Maximum Likelihood A-posteriori un transitorio en el SEP. En base a esto, un detector
(MAP) detector is implemented, which uses de Máxima Probabilidad A-posteriori (MAP)
statistical inference from several hypotheses to es implementado, el cual se basa en la inferencia
determine the probability of a failure from the estadística de varias hipótesis que determinan la
measurements. The algorithm is implemented probabilidad de ocurrencia de falla. El algoritmo
in the DPL software module PowerFactory se implementa en el módulo DPL del software
DIgSILENT. The different experiments prove PowerFactory de DIgSILENT y se realizan
that the developed software is able to distinguish distintas pruebas que demuestran la capacidad de
the most likely scenario from a list of possible la herramienta de software de escoger la hipótesis
scenarios. This work uses simulated- scenarios más probable entre una lista de posibles hipótesis.
and PMU’s measurements. El trabajo presentado es enteramente realizado en
el campo de la simulación, y se emulan datos que se
obtendrían de PMUs.
Index terms− Phasor measurement, Localization Palabras clave− Medición Sincrofasorial, Localización
of failures, Maximum A-posteriori Probability de fallas de SEP, Maximum A-Posteriori Probability
(MAP) (MAP).
17
Edición No 14, Enero 2018
Debido a las condiciones climáticas, diferentes Para el propósito de este estudio, la topología del
agentes externos, y a la longitud de las líneas de sistema, incluyendo los parámetros de sus elementos,
transmisión (L/T), estos elementos son susceptibles son conocidos al igual que las condiciones de estado
a eventos de falla, siendo importante el estudio de estacionario del sistema. Adicionalmente, se puede
este tipo de fenómenos con el objeto de proteger a los establecer una lista completa de posibles eventos
diferentes equipos del Sistema Eléctrico de Potencia de falla de línea, y sus respectivas características de
(SEP) y evitar cortes de servicio [1]. estado estable post-falla a través de un análisis de
contingencia.
La rápida ubicación de un evento de falla ayuda
a mejorar los tiempos de restauración del sistema y La ubicación de una falla es determinada mediante
permite un mejor control en tiempo real de sus efectos. el establecimiento de un evento de falla dentro de
Los operadores de SEP utilizan equipos y sistemas la lista de fallas posibles. Para el efecto se emplean
automáticos de protección que incluyen relés que los datos (emulados) de PMUs reportados una vez
cuentan con algoritmos de localización de fallas en despejada la falla. Matemáticamente la solución al
base a cálculos de impedancias de L/T. Sobre la base problema de ubicación se basa en la modelación lineal
de estos cálculos, una amplia variedad de algoritmos del SEP y en la determinación del vector de medición
de detección de falla han sido desarrollados e de las PMUs.
implementados en varios sistemas eléctricos [2], [3],
[4], [5], [6], [7], [8]. El algoritmo de ubicación de falla se fundamenta
en un detector de Máxima Probabilidad A-posteriori
La ubicación de fallas en L/T realizada por los (MAP). El análisis del detector MAP muestra que
sistemas de protección es eficaz, aunque susceptible el evento de falla se puede identificar mediante la
a errores vinculados al ajuste y calibración de los comparación de los datos de medida de las PMU’s con
equipos y a tiempos de operación. Por otra parte, el las respuestas pronosticadas por el modelo del sistema
sistema de protección no puede proveer información para cada evento de falla [11].
en tiempo real sobre la variación de los parámetros de
la L/T y de las características de carga. 2.1. Planteamiento de Hipótesis
Las unidades de medición fasorial (PMU) Se consideran m hipótesis de los posibles eventos
actualmente permiten tomar mediciones necesarias de falla, la hipótesis se vincula a la ocurrencia de
para realizar el análisis de la operación de los sistemas un evento de falla en una línea de transmisión de la
eléctricos y de eventos de falla en tiempo real. Las red eléctrica y cada hipótesis tiene una probabilidad a
PMUs, son usadas en diferentes sistemas eléctricos priori donde
alrededor del mundo, ya que permiten obtener datos
análogos de voltaje y de corriente sincronizados en 2.2. Modelación del Sistema Eléctrico de Potencia
un GPS de cualquier parte del sistema eléctrico [9]. [12], [13]
18
Chávez et al. / Herramienta de Software para detección y ubicación de fallas en un SEP
(4) donde:
19
Edición No 14, Enero 2018
Las mediciones realizadas por las PMUs tienen el : Vector que contiene las respuestas promedio
siguiente modelo: pronosticadas por el modelo del sistema en el
horizonte de tiempo para el evento de falla i, para
(14) una condición inicial considerada después del
tiempo de compensación.
Donde p es la medición de los ángulos de rotor y
velocidades angulares de un conjunto de generadores (17)
obtenida de las mediciones de las PMUs, después de
un proceso de pre-estimación de estas variables. La : Representa su correspondiente propagación.
matriz tiene valores de 0 o 1 los cuales especifican las
variables de estado dinámico que se miden;
es una matriz de distribución de señal gaussiana de (18)
ruido blanco con media cero y una matriz de varianza
. : Representa el determinante de la matriz Q.
Una forma de decidir qué evento ha ocurrido es : Matriz de observabilidad por cada hipótesis
elegir la hipótesis con la máxima probabilidad, dados
los datos de las PMUs. La técnica de ubicación de
la falla que se presenta se basa en un detector de
máxima probabilidad a posteriori (MAP). El detector (19)
MAP, puede identificar un evento de falla mediante
la comparación de los datos medidos por la PMU con : Matriz con matrices de identidad ordenada de la
las respuestas pronosticadas por el modelo dinámico siguiente forma:
linealizado para cada evento de falla. El objetivo
del detector es decidir entre posibles hipótesis cual (20)
es la más probable. En este proyecto, las hipótesis
consideradas son:
La regla de decisión de la ecuación (16) elige
la hipótesis más probable entre . Se hacen
comparaciones secuenciales en parejas de hipótesis
para encontrar la hipótesis más probable dado el
Cada hipótesis se modela suponiendo que tiene un conjunto de mediciones.
valor conocido de probabilidad a priori. El detector
MAP emplea un enfoque bayesiano para elegir la Por ejemplo, considere tres eventos de falla o
hipótesis con mayor probabilidad a posteriori. En hipótesis: y si es más probable que ocurra
concreto, la selección del detector es la hipótesis que la hipótesis , de acuerdo con la regla de
si: decisión MAP, se comparará de nuevo con y la
más probable entre ellas es la hipótesis seleccionada .
(15) 3. IMPLEMENTACIÓN DEL ALGORITMO
20
Chávez et al. / Herramienta de Software para detección y ubicación de fallas en un SEP
21
Edición No 14, Enero 2018
22
Chávez et al. / Herramienta de Software para detección y ubicación de fallas en un SEP
23
Edición No 14, Enero 2018
24
Chávez et al. / Herramienta de Software para detección y ubicación de fallas en un SEP
25
Artículo Académico / Academic Paper
Abstract Resumen
Index terms− MVMO, PSO, Power System Palabras clave− MVMO, PSO, Estabilizadores
Stabilizer (PSS), Interconnected National System, de sistemas de potencia (PSS), Sistema Nacional
oscillation mode, damping, optimization, stability. Interconectado (SNI), modo de oscilación,
amortiguamiento, optimización, estabilidad.
26
Edición No 14, Enero 2018
27
Arcos et al. / Sintonización de los Estabilizadores de Sistemas de Potencia del Sistema Nacional Interconectado Ecuatoriano
Donde:
28
Edición No 14, Enero 2018
Tabla 4: Modos Oscilatorios del SNI 2018 con PSSs sin Calibrar utilizan unidades equivalentes para representar a
las centrales de generación. Esta definición consiste
Modos de oscilación SNI 2018 sin sintonizar PSSs
Frecuencia de Razón de en colocar una unidad que representa a la totalidad
Modo Amortiguamiento amortiguamiento de una central usando la facilidad disponible en
[Hz] [%] PowerFactory [3].
1259 - 1260 1,70749 2,1795
1265 - 1266 1,6755 2,6529 6. COMPARACIÓN DE RESULTADOS
1467 - 1468 1,0299 4,1412 La comparación de los resultados obtenidos con los
1699 – 1700 0,5990 2,9508 algoritmos MVMO y PSO se basa en tres parámetros:
incremento de amortiguamiento del modo oscilatorio
1711 – 1712 0,6555 4,1857
(análisis de pequeña señal), respuesta del sistema a
1713 - 1714 0,6555 4,1857 transitorios electromecánicos (análisis dinámico del
sistema) y tiempos de convergencia y ejecución.
La determinación de los factores de participación
para cada modo de oscilación permite identificar las 6.1. Análisis de Pequeña Señal
unidades de generación con mayor presencia en cada Este análisis permite determinar el incremento
modo. A manera de ejemplo se presenta el factor de del amortiguamiento de los modos oscilatorios del
participación asociado al modo 1265 – 1266 en la Fig. sistema. Posterior a la aplicación de ambos métodos
3, en donde se observa que las máquinas con mayor de optimización se realizaron sendos análisis de
contribución son Coca Codo Sinclair y Baba. pequeña señal obteniéndose los resultados que se
presentan en la Tabla 6.
Tabla 6: Resultados de Incremento de Amortiguamiento en el
SNI 2018 por Modo
29
Arcos et al. / Sintonización de los Estabilizadores de Sistemas de Potencia del Sistema Nacional Interconectado Ecuatoriano
30
Edición No 14, Enero 2018
31
Arcos et al. / Sintonización de los Estabilizadores de Sistemas de Potencia del Sistema Nacional Interconectado Ecuatoriano
32
Artículo Académico / Academic Paper
Abstract Resumen
This work analyzes the current demand for El trabajo analiza la demanda actual de la zona
the coverage area of Substation 57 Pomasqui de cobertura de la subestación 57 Pomasqui de
owned by the Empresa Electrica Quito (EEQ) propiedad de la Empresa Eléctrica Quito (EEQ), y
and it establishes the forecasted demand at establece la proyección de la demanda al año 2018
2018 by using a technique for microgrids. To do usando el método por microáreas. Para ello, el estudio
this, the study considers two major stages: 1) considera dos grandes etapas: 1) determinación de la
determination of the maximum demand at each demanda máxima de cada microárea, y 2) estimación
microgrid, and 2) estimation of the land use del área habilitada de construcción de cada zona,
based on the Plan de Uso y Ocupacion del Suelo basándose en el Plan de Uso y Ocupación del Suelo
(PUOS) promulgated by Plan Metropolitano de (PUOS) promulgado en el Plan Metropolitano de
Ordenamiento Territorial (PMOT) of the city of Ordenamiento Territorial (PMOT) de la ciudad de
Quito. The results obtained in these two stages Quito. Los resultados obtenidos en estas dos etapas
are used, in conjunction with the logistic curve, to se utilizan, en combinación con la curva logística,
make the forecasting of the demand at 2018. That para realizar la proyección de la demanda al año
is, it is based on the fact that the growth of the 2018. Es decir, se basa en el hecho que el crecimiento
electric demand follows the behavior of the logistic de la demanda eléctrica sigue el comportamiento
curve and it considers the results of the forecast de la curva logística y considera los resultados de
of the demand. This establishes a procedure that la proyección de la demanda. Esto establece un
allows to estimate the demand for the year 2018. procedimiento que permite estimar la demanda
Finally, with the results of the forecasted demand para el año 2018. Finalmente, con los resultados
for 2018, power flow studies are performed using de la demanda proyectada al año 2018, se realizan
the CYMDIST software to estimate the future estudios de flujos de potencia usando el software
operating conditions. This helps detecting possible CYMDIST con el fin de estimar las condiciones
problems and discussing alternative solutions in operativas futuras. Esto permite detectar posibles
the short and medium term distribution planning. problemas y comentar alternativas de solución en
la planificación de distribución a corto y mediano
plazo.
Index terms− Logistic Curve, Microgrids, Palabras clave− Curva logística, Microáreas,
Electricity Demand Forecasting, PUOS, PMOT Proyección de la demanda, PUOS, PMOT y
and Electrical Power Distribution System. Sistema Eléctrico de Distribución.
33
Edición No 14, Enero 2018
34
Jimenez et al. / Estudio por Microáreas de la Demanda Actual y Futura de la Subestación 57 Pomasqui
Desde este sitio, se extiende los ejes norte, sur, este 2.2. Curva Logística
y oeste hasta los límites del área de cobertura de la
EEQ [4]. Para el desarrollo de este trabajo se utiliza la
función logística o curva logística, muy utilizada en
Con el sistema de coordenadas planteado, se traza modelos de crecimiento poblacional, la cual refleja
una red de cuadrículas de un kilómetro cuadrado el comportamiento característico de la demanda
que ocupe la totalidad de la zona de influencia de eléctrica [5] [6]. En la Fig. 2 se muestra en detalle la
la S/E 57 Pomasqui. Se escoge un estudio uniforme curva logística.
y homogéneo con microáreas de un kilómetro
cuadrado debido a que la superficie a analizar es
extensa y la información a manejar elevada. Para el
presente estudio, en el área de cobertura de la S7E
57 Pomasqui se establecen 240 microáreas de un
kilómetro cuadrado.
35
Edición No 14, Enero 2018
Tabla 1: Datos para Simulación Demanda no Coincidente 2.4. Plan de Uso y Ocupación del Suelo
36
Jimenez et al. / Estudio por Microáreas de la Demanda Actual y Futura de la Subestación 57 Pomasqui
Una vez que se identifica la microárea saturada Finalmente, la demanda adicional por concepto
se debe determinar el factor de saturación, el cual de cocinas de inducción se obtiene al multiplicar la
resulta de la relación entre el área neta habilitada de la demanda individual de la cocina seleccionada (2,8
microárea saturada con la demanda eléctrica existente kW) por el número de cocinas, y por el factor de
en la misma, la expresión viene dada en [kVA/m²]. coincidencia. A continuación, se muestra un ejemplo
En la Tabla 5 se definen los factores de saturación de cálculo para la microárea NE0721 que cuenta con
para la demanda máxima no coincidente y demanda 1 558 usuarios residenciales.
máxima coincidente según el tipo de uso de suelo.
Ya que se espera que un 80% de los usuarios tenga
Tabla 5: Factores de Saturación por Tipo de Uso de Suelo cocina de inducción, se obtiene que el número de
usuarios reales de cocinas de inducción es 1 246.
Demanda máxima no coincidente Demanda máxima coincidente
De-
Factor de Demanda Factor de El cálculo del factor de coincidencia se lo realiza
manda
Uso vigente saturada
saturación Uso vigente saturada saturación
tomando en cuenta la demanda diversificada del
[kVA] [kVA/m²] [kVA] [kVA/m²] número de usuarios residenciales de cada microárea,
Residencial 2.301,00 5,8912E-03 Residencial 1.818,00 4,5924E-03 en este caso 1 246 usuarios, dividido para la demanda
Industrial 315 6,6790E-04 Industrial 287 6,0850E-04 máxima diversificada de un usuario, que es igual
Agrícola Agrícola
a 0,8 kW y 3,7 kW respectivamente [3] [6]. Para el
125 5,3820E-04 82 3,5300E-04
Residencial Residencial procedimiento detallado se usa (1):
Equipa- Equipa-
18 2,1675E-03 14 1,6859E-03
miento miento
Protección
15 2,6068E-03
Protección
10 1,7379E-03
(1)
Ecológica Ecológica
El PEC busca reemplazar el uso de Gas Licuado El factor de coincidencia para este caso es 0,2162.
de Petróleo (GLP) por electricidad para la cocción Para calcular la demanda por cocinas de inducción se
y calentamiento de agua en el sector residencial. El emplea (2).
ingreso de este tipo de cargas representa un incremento
adicional de la demanda eléctrica que debe ser (2)
incorporada a la proyección de cada microárea que
cuente con algún usuario residencial. Se espera que Donde:
el plan tenga una cobertura final superior al 80% del
total de los abonados residenciales [3] [6].
37
Edición No 14, Enero 2018
continuación, a modo de ejemplo se desarrolla el en el año 2035 (en un plazo de 21 años tomando como
proceso mencionado para el escenario de demanda no año cero el 2014), para las áreas establecidas como
coincidente usando la información correspondiente RNR, agrícola residencial, equipamiento y protección
de la Tabla 4 y 5. El mismo proceso se replica para el ecológica se define 2050 como año de saturación (en
escenario de demanda coincidente. un plazo de 36 años tomando como año cero el 2014).
38
Jimenez et al. / Estudio por Microáreas de la Demanda Actual y Futura de la Subestación 57 Pomasqui
El mismo procedimiento se aplica para el Tabla 7: Demanda Proyectada por Primario - Demanda No
resto de microáreas, considerando la demanda del Coincidente
factor de saturación acorde al tipo de uso de suelo
predominante. Demanda total
Demanda proyectada por fase
proyectada Factor de
Alimentador [kVA]
2018 potencia
Además, en base a la información analizada de la [kVA] A B C
proyección de la demanda saturada en cada microárea, 57A 9 519,74 3 337,10 3 166,39 3 016,24 0,947
se observan algunos casos en los cuales ésta supera 57B 8 890,18 2 861,72 2 923,59 3 104,86 0,967
a la demanda usada como referencia del factor de 57C 9 757,25 3 233,85 3 364,88 3 158,52 0,961
saturación. Este hecho imposibilita el cálculo directo 57D 4 149,10 1 389,83 1 320,60 1 438,67 0,963
de la demanda al año 2018, por lo que se desarrolla 57E 13 914,69 4 560,81 4 502,57 4 851,32 0,975
un proceso adicional al detallado anteriormente. Este 57F 7 078,85 2 118,83 2 676,91 2 283,11 0,969
proceso adicional se lo desagrega en [12]. 57G 11 048,22 3 365,69 3 737,49 3 945,04 0,947
39
Edición No 14, Enero 2018
40
Jimenez et al. / Estudio por Microáreas de la Demanda Actual y Futura de la Subestación 57 Pomasqui
Para el transformador T2 los primarios 57A y Tabla 12: Transformadores Sobrecargados Demanda Máxima
57E tienen cargabilidades de 81,63% y 117,60% No Coincidente
respectivamente, este último ya supera el 100%
de cargabilidad y probablemente presentará
inconvenientes en su operación normal. Los primarios
restantes, 57D y 57F presentan cargabilidades
menores al 60% lo cual indica un funcionamiento
adecuado. El transformador en el primario 57C es el mismo
que presenta una cargabilidad alta en la simulación
Los primarios que sobrepasan el 80% de para el año 2015 pero en este caso, supera el 30% de
cargabilidad asociado al transformador T1 son los su capacidad nominal, mientras los transformadores
alimentadores 57C con 81,97% y 57G con 93,83%, agrupados en el alimentador 57E tienen cargabilidades
respectivamente, lo cual indica que en unos pocos que sobrepasan a la capacidad nominal en un 16%.
años en el futuro a partir del 2018 estarán al límite
de su capacidad. Lo mismo sucede con el primario De igual manera, para el escenario de demanda
57B que con un 75,47% se encuentra en un rango coincidente se cuenta con tres transformadores
aceptable pero muy próximo a bordear los límites sobrecargados, los cuales se detallan en la Tabla 13.
establecidos. De los transformadores de distribución restantes, 2
presentan cargabilidades entre el 80% y 60%; y los
Para el escenario de demanda coincidente demás tienen cargabilidades inferiores al 60%.
mostrado en la Tabla 11 se observa una tendencia
similar, el primario 57E sobrepasa el 100% de Tabla 13: Transformadores Sobrecargados Demanda Máxima
cargabilidad y dependiendo de la clase del TC puede Coincidente
provocar la operación de su respectiva protección. Los
primarios 57D y 57F se mantienen con cargabilidades
bajas mientras que los alimentadores restantes (57A,
57B, 57C, 57G) se encuentran entre el rango de 60%
y 80%, aceptable para la operación normal de la
subestación. El análisis de cargabilidad de los transformadores
de distribución indica que la gestión del sistema de
Tabla 10: Cargabilidad Alimentadores S/E 57 Pomasqui Año
distribución en los próximos años debe dar prioridad
2018-Dem. No Coincidente
al cambio de los equipos que se muestran en la Tabla
12 y en la Tabla 13.
4. CONCLUSIONES Y RECOMENDACIONES
41
Edición No 14, Enero 2018
42
Jimenez et al. / Estudio por Microáreas de la Demanda Actual y Futura de la Subestación 57 Pomasqui
43
Artículo Académico / Academic Paper
[Link]@[Link]; [Link]@[Link]
2
Coorporacion Nacional de Electricidad CNEL EP, Unidad de Negocio Sucumbíos
[Link]@[Link]
Abstract Resumen
44
Edición No 14, Enero 2018
45
Guamán et al. / Estudio de Coordinación de Protecciones del Sistema de Subtransmisión de CNEL EP Sucumbíos
Máximas Mínimas
corrientes de C.C. corrientes de C.C.
BARRA 69 kV
3 1
1 [kA] 3 [kA]
[kA] [kA]
Jivino de CNEL 17.454 15.167 15.355 13.332
Sacha 7.737 5.885 6.559 5.005
Francisco de Orellana 3.509 3.251 2.93 2.745
Payamino 3.48 3.201 2.906 2.703 Figura 2: Respuesta del Esquema PUTT ante una Falla al 90%
con Respecto a la Subestación A
46
Edición No 14, Enero 2018
Cuando la falla ocurre al 50% de la línea como se Figura 6: Impedancia Calculada por el Relé 21 de la
Subestación A ante una Falla de Alta Impedancia
observa en la Fig. 5, ambos relés de distancia de la
subestaciones A y B opera en Zona 1, envían señal de La impedancia que calcula el relé 21 ubicado en la
disparo y disparan sus interruptores instantáneamente, subestación B en el diagrama R-X se puede observar
no se necesita confirmación de señal del extremo en la Fig. 7.
remoto.
47
Guamán et al. / Estudio de Coordinación de Protecciones del Sistema de Subtransmisión de CNEL EP Sucumbíos
Puede darse el caso que exista una falla de alta Ante fallas de alta impedancia y si la línea es muy
impedancia que provoque que los relés de distancia corta, el relé de distancia puede ver la falla en Zona 2,
de las dos subestaciones vean la falla en Zona 2 lo que provocaría que en el esquema de teleprotección
como se ilustra en la Fig. 8, lo que provoca que en PUTT nunca se envié la señal de disparo para que
el esquema PUTT no se envíe la señal de disparo al opere en Zona 2 acelerado, por lo que el relé de
extremo remoto y por lo tanto sus protecciones operen distancia operaría de manera temporizada. Por ello se
en tiempos de Zona 2 (entre 0.3s y 0.6s [5]), por lo puede concluir que el esquema POTT es mucho más
cual no se despeja la falla de manera instantánea. confiable que el esquema PUTT para fallas de alta
impedancia.
En el caso del sistema de subtransmisión de CNEL
EP Sucumbíos, debido a que las líneas donde se van
a implementar los relés de distancia con esquema de
teleprotección son relativamente cortas, se recomienda
utilizar el esquema POTT para detección de fallas de
alta impedancia, contando además como respaldo
con una protección de sobrecorriente direccional de
neutro.
En la Tabla 5 resume los ejemplos expuestos
anteriormente y la Tabla 6 muestra la comparativa
entre los dos esquemas de teleprotección.
48
Edición No 14, Enero 2018
a implementar los relés de distancia con esquema de • Para la coordinación, se considera que existe
teleprotección son relativamente cortas se recomienda pérdida de los canales de comunicación para
utilizar el esquema POTT para la detección de fallas el relé de distancia debido a que en la práctica
de alta impedancia con respaldo de 67 N para tener operaría de manera instantánea por el esquema
una mayor sensibilidad a este tipo de fallas. de teleprotección propuesto (POTT), por
ello actúa como protección de respaldo y
5. COORDINACIÓN DE PROTECCIONES para ciertos eventos de falla que se describe
a continuación actúa en diferentes zonas de
Antes de realizar el análisis correspondiente, es protección, esto tiene por objetivo verificar la
necesario aclarar cómo se establece la coordinación, correcta coordinación de las protecciones de
debido a que el sistema de subtransmisión posee respaldo.
ciertos tramos radiales y mallados.
• Se sugiere la instalación de relés de bajo
Para tramos radiales se plantea de la siguiente voltaje en terminales de carga de una línea.
manera: Estos relés de bajo voltaje no se coordinan
con las demás protecciones, pues su objetivo
• Protección principal: 50/51 de fase y neutro,
para todos los eventos de falla que ocurran no es la detección y despeje de falla, sino
dentro de la zona protegida. la desconexión completa de la línea en
• Protección de respaldo: 50/51 de fase y condiciones de post-falla.
neutro del elemento aguas arriba para La presencia de relés de bajo voltaje (27), se debe
todos los eventos de falla en caso de que la a que por ser un sistema radial, hay extremo fuerte
protección principal no opere en los casos que y débil en algunos puntos de la red. Se implementa
corresponda. el relé de bajo voltaje en el extremo débil debido
Para tramos mallados o anillos se plantea de la a que para cualquier tipo de falla la corriente de
siguiente manera: cortocircuito en el extremo débil es cero (menor a
la corriente nominal) y cualquier otro tipo de relé
• Protección principal: 87L, para todos los operaría de manera errónea en este extremo. Para
eventos de falla que ocurran dentro de la zona poder aislar el elemento en falla el relé de bajo voltaje
protegida. operaría luego de un tiempo prudencial en que la falla
• Protección principal: 21 (F/N) para todos los ha sido despejada desde el extremo fuerte. El relé 27
eventos de falla que ocurran dentro de la zona
opera mediante la detección de la caída de voltaje en
protegida con canales de comunicación.
el extremo débil (el punto donde éste se instala) como
• Protección de respaldo: 21F, 67 (F/N), para
todos los eventos de falla en caso de que la se observa en la Fig. 21.
protección principal no opere en los casos que En base a los criterios mencionados anteriormente
corresponda. y luego de realizar los estudios requeridos mediante
Es necesario aclarar que no se sugiere los relés PowerFactory se proponen los siguientes ajustes para
21N en las subestaciones Lago Agrio, Jivino de los relés de bajo voltaje (27) instalados en el sistema
CNEL y Jivino de TRANSELECTRIC, debido a de subtransmisión de CNEL EP Sucumbíos como se
que las líneas a proteger son relativamente cortas y puede observar en la Tabla 7.
para fallas con contacto a tierra mediante una alta
Tabla 7: Relés de Bajo Voltaje para el Sistema de
impedancia de falla, estos relés no operan debido al Subtransmisión de CNEL EP Sucumbíos
subalcance presentado.
Esto se debe a que la resistencia de falla al ser PROTECCIÓN DE BAJO VOLTAJE
la suma de la resistencia de arco más la resistencia Voltaje Tiempo
Línea S/E
con contacto a tierra puede llegar a ser muy (pu) (s)
elevada (experiencia del operador del sistema de Coca-Sacha
Francisco de
0.85 3
subtransmisión de Sucumbíos, máxima de 8 Ω, por Orellana
ser zona selvática y con elevados índices de incidencia Sacha-Jivino Sacha 0.85 3
de falla en ciertos tramos) provocando una variación Tarapoa-
muy significativa de la resistencia (en especial para Tarapoa 0.85 3
Shushufindi
líneas cortas) y el cambio del ángulo de fase, lo cual Shshufindi-
Shushufindi 0.85 3
provocaría que el relé no opere por subalcance, Jivino
las pruebas realizadas mediante simulaciones con Lago Agrio-
Lago Agrio 0.85 3
Jivino
impedancia de falla desde los 10 Ω a 15 Ω confirman
Celso
que los relés 21N no detectan la falla, comprobando
Castellano- Celso Castellanos 0.85 3
lo que la literatura técnica indica [8] y recomienda Lago Ario
para este tipo de líneas, establecer como protección Lumbaquí-
de respaldo un 67N para una mayor sensibilidad con Celso Lumbaquí 0.85 3
fallas de alta impedancia que puede darse el caso. Castellanos
49
Guamán et al. / Estudio de Coordinación de Protecciones del Sistema de Subtransmisión de CNEL EP Sucumbíos
50
Edición No 14, Enero 2018
PROTECCIÓN
TIPO DE Resistencia
Figura 15: Respuesta de la Protección de Distancia de Fase S/E BAHÍA 87L 211 671 67N1
FALLA de falla [Ω]
(21) Ubicado en la S/E Jivino de CNEL ante un Evento de
3F 0 SI SI SI NO
Falla Trifásico
Jivino de 2F 5 SI SI SI NO
Jivino de
Como se puede observar en la Fig. 16, el relé de TRANS
ELECTRIC CNEL 69
1F 10 SI NO SI SI
kV
sobrecorriente direccional de fase (67, protección de
2F-T 10 SI NO SI SI
respaldo) opera de manera temporizada, lo cual es
correcto debido a que la falla ocurre en la dirección Falla al 50% desde la S/E Jivino de CNEL
2F-T 20 SI NO SI SI
51
Guamán et al. / Estudio de Coordinación de Protecciones del Sistema de Subtransmisión de CNEL EP Sucumbíos
52
Edición No 14, Enero 2018
53
Aplicación Práctica / Practical Issues
Abstract Resumen
This work presents an analysis of the power Este trabajo presenta un análisis de la evaluación
system transient stability assessment using the de la estabilidad transitoria utilizando como
concept of Center of Inertia COI as reference. referencia el Centro de Inercia COI del sistema de
First, a synchronous machine model for transient potencia. En primer lugar, se presenta un modelo
stability studies has been presented. This model de máquina síncrona utilizado en estudios de
uses the dq0 reference frame based on Park estabilidad transitoria, cuyo modelo se basa en un
transform. Afterwards, the importance of marco de referencia dq0 utilizando la transformada
choosing a specific angle reference for defining de Park. Posteriormente, se analiza la importancia
any Park-based model is analyzed. This aspect de la elección de un ángulo de referencia como
has several important consequences that motivate un paso indispensable para definir cualquier
the selection of a general time-variable angle modelo basado en la transformada de Park. Esto
reference known as Center of Inertia COI. With tiene algunas consecuencias importantes, que
the purpose of generalizing the benefits of using motivarán la elección de un ángulo de referencia
COI for transient stability assessment, formal general variable en el tiempo conocido como
mathematical demonstrations are developed. In el Centro de Inercia COI. Con el propósito de
addition, simulation results of several study cases generalizar los beneficios del uso del COI para
analyzed at two test power systems are presented. realizar análisis de estabilidad transitoria, se
By this way, it is verified that COI gives a general desarrollan demostraciones matemáticas formales
reference frame properly designed for assessment y se presentan resultados de simulaciones de
of power system transient stability. varios casos de estudio realizados en dos sistemas
de prueba. De esta forma, se verifica que el uso
del COI brinda el marco referencial apropiado
para analizar estabilidad transitoria en sistemas
eléctricos de potencia.
Index terms− Center of Inertia, Rotor Angle, Palabras clave− Ángulo de Rotor, Centro de Inercia,
Synchronous Machine, and Transient Stability. Estabilidad Transitoria, Máquina Sincrónica.
54
Edición No 14, Enero 2018
55
Echeverría y Cepeda / Evaluación de Estabilidad Transitoria de Sistemas de Potencia utilizando el concepto de COI
1.1.1 El estator
El estator es la parte estacionaria, hecho en forma
cilíndrica, provisto en el interior de la superficie de
(1)
ranuras equidistantes y su respectivo devanado para
las tres fases (a, b y c), desplazadas en espacio entre
sí por 120 grados eléctricos [3].
Generalmente, los ángulos se miden en radianes
1.1.2 El rotor eléctricos o grados eléctricos. El ángulo cubierto por
El rotor es la parte móvil, que gira dentro del un par de polos es 2π rad o 360 grados eléctricos (ver
estator, el cual contiene el devanado de campo Fig. 2).
(excitación, alimentado con corriente continua) y
devanados de amortiguamiento formados por barras
conectadas a anillos en los extremos del rotor que
forman la denominada jaula de ardilla. La corriente
continua en el bobinado de campo produce un
campo magnético, que es estacionario con respecto
al rotor, pero que gira con respecto al estator con
una velocidad igual a la velocidad del rotor. De esta
Figura 2: La Relación entre Ángulos Mecánicos y
manera, voltajes alternos, desplazados por 120º en el
Eléctricos [3]
tiempo, son inducidos en los devanados del estator,
a una frecuencia dada por la velocidad del rotor y el Para establecer la relación entre grados eléctricos
número de pares de polos magnéticos [3]. usados en expresiones de onda de voltaje y corriente,
y grados geométricos que definen la posición del
1.1.3 Principio de Funcionamiento
estator, la ecuación se multiplica por 2πt, obteniendo:
Cuando el generador está conectado a la red, los
devanados del estator producen corrientes alternas (2)
con una pulsación igual a la de la red. Estas corrientes
generan a su vez un campo magnético giratorio (3)
que interactúa con el campo producido por la
corriente de excitación, dando como resultado un par Lo que resulta:
electromagnético (Te) que, en condiciones de estado
estacionario, es igual al par mecánico aplicado al rotor (4)
(Tm). Para obtener un par constante, tanto los campos
magnéticos del estator como los campos magnéticos 2.2. Modelo Electromecánico del Generador
del rotor deben girar con la misma velocidad (para ser Sincrónico: Ecuación de Oscilación
sincrónicos) [3].
La ecuación de oscilación (o ecuación de movimiento)
Si Tm no es igual a Te, la velocidad del rotor y que gobierna el movimiento del rotor de una máquina
por lo tanto la velocidad del campo electromagnético sincrónica se basa en un principio elemental de
del rotor es diferente de la velocidad del campo dinámica que establece que el par de aceleración es
electromagnético del estator. Por lo tanto, el producto del momento de inercia del rotor por su
corrientes parásitas se producen en las barras aceleración angular [4],[5].
de amortiguamiento, con una frecuencia igual a
la frecuencia de deslizamiento. Estas corrientes
inducidas producen un campo electromagnético que, (5)
según la ley de Lenz, se opone al cambio que las ha
producido. Por lo tanto, si la velocidad del rotor es
sobre-sincrónica, las corrientes inducidas causan un Donde J es el momento total de inercia de la masa
aumento del par electromagnético Te, y el rotor se del rotor (en kg-m2) y qm es el desplazamiento del
desacelera. De lo contrario, si la velocidad del rotor rotor con respecto al eje estacionario, en radianes
es sub-sincrónica, las corrientes inducidas causan la mecánicos (rad).
disminución del par electromagnético y el rotor se
Bajo la operación en estado estable del generador,
acelera [3].
Tm y Te son iguales y el par de aceleración Ta es cero.
El número de los polos magnéticos (P) determina la En este caso, no hay aceleración o desaceleración de
velocidad mecánica wm del rotor y la frecuencia f de las la masa del rotor y la velocidad constante que resulta
corrientes del estator. Denotando por p=P/2 el número es la velocidad sincrónica [4].
de pares de polos, y n la velocidad del rotor (rpm), la
Como se mide con respecto al eje de referencia
frecuencia (Hz) de la fuerza electromotriz inducida
estacionario sobre el estator, es una medición absoluta
viene dada por la ecuación , que en estado estacionario
del ángulo del rotor. En consecuencia, continuamente
es igual a la frecuencia del sistema de potencia [3].
se incrementa con el tiempo aun a velocidad
56
Edición No 14, Enero 2018
(7)
Donde V es la magnitud máxima de la onda y
es un ángulo eléctrico que se integra a partir la
y, frecuencia eléctrica base:
(8)
(13)
Se recuerda, de la dinámica elemental, que la
Dado que no hay una opción natural para escoger
potencia es igual al par por la velocidad angular y así,
t0 en un sistema que opera continuamente, 0 debe ser
al multiplicar la ecuación por , y reemplazando en
elegido arbitrariamente [6]. Obsérvese que, para una
se obtiene [4]:
forma de onda de tiempo dada como en , la elección
de 0 afecta el valor de .
(9) Los voltajes trifásicos, denotados va, vb y vc, una
vez aplicada la transformada de Park aparecen como:
57
Echeverría y Cepeda / Evaluación de Estabilidad Transitoria de Sistemas de Potencia utilizando el concepto de COI
orientación del eje de campo directo en las máquinas 4. MARCO DE REFERENCIA GENERAL
i y j, y este ángulo es independiente del ángulo de DADO POR EL COI
referencia particular que se elija. Además, dada una El Centro de Inercia (COI) ofrece un marco
medición de voltaje o corriente arbitraria (como en la de referencia físico que permite representar de la
ecuación (12)), el ángulo de fase también depende mejor manera posible el problema de estabilidad
de la elección de las coordenadas base (ángulo de transitoria. Además, la transformación de la ecuación
referencia), de manera que todos los ángulos relativos de oscilación en las coordenadas del centro de inercia
se conservan si el sistema de coordenadas base es no sólo ofrece un marco físico para la formulación,
cambiado. Dado que todas las cantidades del sistema sino que también elimina la energía asociada con el
dependen solamente de los ángulos relativos, la movimiento del centro inercial que no contribuye a
elección de las coordenadas base es verdaderamente la determinación de estabilidad. Por lo tanto, es muy
arbitraria [6]. recomendable referenciar la respuesta dinámica del
Analizando la Fig. 3, se revela un hecho importante sistema al COI para obtener un comportamiento de
con respecto al sistema de referencia de la máquina. estabilidad transitoria adecuado [7].
Está claro que todos los ángulos de las máquinas (di, En esta sección se presenta una demostración
para el i-ésimo generador) se calculan con respecto matemática formal que el Centro de Inercia COI ofrece
a las coordenadas de red. En otras palabras, los un marco de referencia general independientemente
ángulos del rotor se definen como desplazamientos de cuál es la referencia de ángulo utilizada por el
desde el ángulo de referencia arbitrario [6]. Por lo algoritmo de integración. Esta consideración se
tanto, la elección de un ángulo de referencia es un analizó basándose en el hecho de que DIgSILENT
paso indispensable para definir cualquier modelo PowerFactory utiliza el ángulo de rotor del generador
basado en la transformada de Park. Esto tiene algunas oscilante como referencia en el proceso de integración
consecuencias importantes, que motivarán la elección en el dominio del tiempo (como se explicó en la
de un ángulo de referencia variable en el tiempo Sección 3), mientras que la teoría tradicional define
conocido como el Centro de Inercia COI, que se una referencia síncrona arbitraria asociada a la
tratará más adelante. velocidad nominal síncrona wn.
58
Edición No 14, Enero 2018
(19)
Sea δ’COI System el ángulo de rotor equivalente del
COI del sistema calculado en esta referencia (COI’),
entonces:
(20)
(21)
Figura 4: Diagrama del Sistema WSCC – 9 barras, 3
Generadores
(22)
(23)
(24)
(25)
(26)
Figura 5: Diagrama del Sistema New England – 39 Barras, 10
El ángulo del rotor del i-ésimo generador referido Generadores
al COI, referido a su vez a COI’, se define como sigue:
5.1. Evaluación de Estabilidad Transitoria del
Sistema WSCC
(27)
Se desarrolló una herramienta computacional que
permite calcular el ángulo del Centro de Inercia COI
(28) para ser utilizado como ángulo de referencia en el
sistema de potencia. El software es implementado en
(29) el lenguaje de programación DPL de PowerFactory,
versión 15 [9], y permite analizar la seguridad
Comparando y , es fácil concluir que Entonces, operativa de un sistema de potencia a partir del análisis
se ha demostrado que el cálculo de de estabilidad transitoria basado en el concepto de
COI da una referencia general que es independiente Centro de Inercia.
de la referencia de ángulo utilizada por el algoritmo Para probar lo analizado en la Sección 3 y lo
de integración. demostrado en la Sección 4, sobre la importancia de
definir un ángulo de referencia general, a continuación,
5. APLICACIÓN Y RESULTADOS se presentan los resultados de las simulaciones de una
Para comprobar la efectividad de utilizar el marco falla trifásica en línea 7-5 del sistema WSCC, la cual
de referencia general denominado COI para el análisis es despejada a los 180 ms, y se analizan los resultados
de estabilidad transitoria, se realizan simulaciones en de los ángulos de rotor obtenidos en el Programa
dos sistemas de prueba con el propósito de evaluar PowerFactory de DIgSILENT.
su estabilidad transitoria: i) el sistema de potencia de Caso 1: Ángulo de referencia G1
prueba de la Western System Coordinating Council
(WSCC), de 9 barras y 3 generadores [8], ver Fig. 4; y, La Fig. 6 presenta los resultados de cada uno de
ii) el sistema de potencia de prueba de New England, los ángulos de rotor del sistema de 9 barras. En esta
de 39 barras y 10 generadores [7], ver Fig. 5. gráfica se aprecia que el ángulo de rotor del G1 es
cero, debido a que es la máquina de referencia. En
cambio, los generadores G2 y G3 presentan pérdida de
sincronismo, dado que sus ángulos crecen hasta más
de 180º, que es el valor establecido como referencia
para determinar pérdida de sincronismo.
59
Echeverría y Cepeda / Evaluación de Estabilidad Transitoria de Sistemas de Potencia utilizando el concepto de COI
Figura 7: Velocidad de las Máquinas, con G1 de Referencia Figura 9: Velocidad de las Máquinas, con G2 de Referencia
La Fig. 8 presenta los resultados de cada uno de La Fig. 10 presenta los resultados de los ángulos
los ángulos de rotor del sistema de 9 barras cuando de rotor referenciados al Centro de Inercia COI del
el G2 es la referencia. En este sentido, se aprecia que sistema de 9 barras. Cabe destacar que las mismas
el ángulo de rotor del G2 es cero, debido a que es la señales son obtenidas usando como datos de entrada
máquina de referencia. En cambio, el generador G1 ya sean los resultados del Caso 1 o del Caso 2.
aparentemente presenta pérdida de sincronismo, dado
que el valor de su ángulo crece hasta más de 180º.
Estos valores difieren de los obtenidos en el
Caso 1, y no permiten tener un claro análisis de
estabilidad transitoria del sistema de 9 barras a partir
de los resultados obtenidos del comportamiento
de los ángulos de rotor, razón por la cual dejan de
ser confiables para obtener una conclusión directa
de cuáles son las máquinas que realmente pierden
sincronismo.
La Fig.9 presenta el comportamiento de la
velocidad de las máquinas sincrónicas, donde
Figura 10: Ángulos de Rotor Referenciados al Centro de Inercia
se aprecia que el sistema se divide en dos áreas
coherentes, la primera establecida por el generador En la gráfica se aprecia claramente cuál es el
G1 y la segunda por los generadores G2 y G3. verdadero comportamiento de los ángulos del sistema
Es importante resaltar que el comportamiento de y los generadores G2 y G3 son los que realmente
60
Edición No 14, Enero 2018
61
Echeverría y Cepeda / Evaluación de Estabilidad Transitoria de Sistemas de Potencia utilizando el concepto de COI
62
Edición No 14, Enero 2018
63
Aplicación Práctica / Practical Issues
Abstract Resumen
The inclusion of high consumption electric loads La inclusión de cargas eléctricas de alto consumo
such as induction stove and electric heaters has como cocinas de inducción y calentadores eléctricos
made it necessary to install a second phase in ha hecho necesaria la instalación de una segunda
many homes. This brings with it the problem of fase en muchos domicilios. Esto trae acarreado
unbalances in the currents of each phase. The el problema de desbalances en las corrientes
present work proposes the use of two inverter de cada fase. El presente trabajo propone el uso
topologies acting as active two-phases power filters. de dos topologías de inversores actuando como
In addition to the reactive power compensation filtros activos de potencia bifásicos. Además de la
capability, these filters can balance active power capacidad de compensación de potencia reactiva
in each phase for 120-phase and 180-degree two- estos filtros pueden balancear la potencia activa en
phases systems. cada fase para sistemas bifásicos con 120° y con
180° entre fases.
Index terms− Back to Back, active power filter, Palabras clave− Back to Back, filtro activo de
two-phase systems, pulse width modulation, load potencia, sistemas bifásicos, modulación de
balance, load unbalance. ancho de pulso, cargas balanceadas, cargas
desbalanceadas.
64
Edición No 14, Enero 2018
65
Fajardo et al. / Filtro Activo de Potencia para Instalaciones Bifásicas con Capacidad de Balance de Corrientes
neutro. Esta topología, al igual que la BtB mostrada Si se utiliza este mismo principio para un sistema
en la Fig. 3 tiene capacidad para compensar tanto bifásico con desfase de 120º entre fases las expresiones
potencia reactiva como activa. Algunas de las ventajas de las moduladoras enviadas a cada rama son:
de esta topología con relación a la BtB de puentes H
son: menor dispositivos de conmutación y no tener (7)
estados de conmutación que produzcan cortocircuitos
en el bus DC, por lo que no es necesario dividir los
inductores para contrarrestar este efecto. (8)
La topología propuesta requiere un rediseño de la Sin embargo, a diferencia del sistema trifásico,
modulación de ancho de pulso enviada a cada rama para obtener voltajes alternos entre las líneas VAC y
del puente para poder obtener tensiones apropiadas a VBC lo único que se requiere es eliminar la componente
los distintos sistemas bifásicos. continua de los voltajes de VAN y VBN, por lo tanto, el
voltaje de la tercera fase que puede cumplir con este
requerimiento será:
(9)
(10)
Figura 4: Filtro Bifásico Propuesto
(11)
3.1. Modulación en la Tercera Rama de Neutro
En un puente trifásico operando en condiciones Se utiliza el mismo concepto del sistema trifásico
equilibradas las expresiones de las moduladoras para los sistemas bifásicos con desfase de 180º entre
enviadas a cada rama de puente fueron tomadas de fases y las expresiones para el voltaje de fase y de
[10] y son: línea son las siguientes.
(12)
(1)
(13)
(2)
(14)
(3)
(15)
Los voltajes de correspondientes de línea resultan: Adicionalmente con las ecuaciones (10), (11) y
(15) se puede determinar que se requiere al menos un
bus DC de 350V para obtener un voltaje de red de
(4) 175Vpp en la carga, cabe recalcar que esto sería con el
índice de modulación “m” en su máximo de 1 y que
además la tensión de 175Vpp es para estar justo por
(5)
encima del valor pico de la tensión estándar de línea.
La topología y la estrategia de modulación
(6) propuestas se prueban conectando cargas resistivas de
20W y un bus de 300VDC. La Fig. 5 muestra el valor
medio de las tensiones de salida cuando se modula un
En , (2) y (3) se observa que los voltajes de fase sistema bifásico con separación de 120°.
tienen una componente continua y una alterna, pero
en los voltajes de línea expresados en (4), (5) y (6)
la componente continua se elimina al restar una fase
de la otra.
66
Edición No 14, Enero 2018
En la Fig. 6 se puede observar que se produce un Figura 7: Circuito Equivalente del Filtro Activo Bifásico
voltaje sinusoidal en la carga de aproximadamente
150VDC. La topología fue simulada con modulaciones donde Ts es el periodo de muestreo e iAPF es la corriente
para redes de 120º de desfase como de 180º teniendo de compensación del inversor. Esta expresión es
un buen funcionamiento para ambos casos. idéntica para ambas fases. Si se reemplaza iAPF(k+1)
por la referencia de corriente deseada se puede
despejar el valor de tensión que el inversor debe
sintetizar
(17)
67
Fajardo et al. / Filtro Activo de Potencia para Instalaciones Bifásicas con Capacidad de Balance de Corrientes
68
Edición No 14, Enero 2018
69
Fajardo et al. / Filtro Activo de Potencia para Instalaciones Bifásicas con Capacidad de Balance de Corrientes
Figura 15: Corrientes Compensadas para las Fases “A” y “B” REFERENCIAS BIBLIOGRÁFICAS
[4] Leng S., Liu W., Chung I. Y., Cartes D., “Active
power filter for three-phase current harmonic
cancellation and reactive power compensation,”
in 2009 American Control Conference, pp. 2140–
2147, June 2009.
[6] Fajardo M., Viola J., Restrepo J., Quizhpi F., Aller
Figura 17: Corrientes Compensadas para las Fases “A” y “B” J., “Two-phase active power filter direct current
control with capacitor voltages estimation and
6. CONCLUSIONES balance,” in 2015 IEEE Workshop on Power
Electronics and Power Quality Applications
El filtro de potencia activo propuesto con una (PEPQA), pp. 1–6, June 2015.
topología para sistemas bifásico es capaz de balancear
potencia activa entre las fases además de cumplir con [7] Alcala J., Cardenas V., Miranda H., Perez-
el requerimiento de eliminación de armónicos de Ramirez J., Charre S., “A three-phase back-to-
corriente y mejora del factor de potencia, permitiendo back converter for reactive power compensation,
un equilibrio de potencia entre las fases. El control se current harmonic filtering and active power
basa en un algoritmo de control directo de corriente compensation,” in 2013 IEEE Energy Conversion
con un controlador PI para estabilizar el voltaje del Congress and Exposition, pp. 2371–2377, Sept
bus DC. El APF se simulo ́ para el funcionamiento en 2013.
70
Edición No 14, Enero 2018
[8] Tang C. Y., Chen Y. F., Chen Y. M., Chang Y. R., Julio César Viola.- Nació en La
“DC-link voltage control strategy for three-phase Paz (ER), Argentina, en 1975.
back-to-back active power conditioners,” IEEE Obtuvo el título de Ingeniero
Transactions on Industrial Electronics, vol. 62, Electrónico en el año 2000 en
pp. 6306–6316, Oct 2015. la Universidad Tecnológica
Nacional, Paraná (ER), Argentina,
[9] Viola J., Fajardo M., Aller J. M., Quizhpi F., y el título de Doctor en Ingeniería
Restrepo J., Ginart A., “Back-to-back active en el año 2008 en la Universidad
power filter for current balancing in two-phase Simón Bolívar, Caracas,
systems,” in 2017 Ninth Annual IEEE Green Venezuela. Desde el año 2012 ha estado en la Universidad
Technologies Conference (GreenTech), pp. 139– Politécnica Salesiana como Investigador del Proyecto
144, March 2017. Prometeo – SENESCYT y desde el 2016 está vinculado
como profesor titular. Sus intereses en investigación se
[10] Rashid M., “ELECTRONICA DE POTENCIA, centran en la electrónica de potencia para el control de
Circuitos, dispositivos y aplicaciones”. México. máquinas eléctricas y compensación de energía reactiva,
3a Edición. Prentice Hall, 2004. controles neuronales y difusos y energías alternativas,
áreas en las cuales posee más de 40 artículos científicos
[11] Gimenez M., Guzman V., Restrepo J., Aller arbitrados.
J., Bueno A., Viola J., Millan A., Cabello A.,
“PLATAFORMA: Development of an integrated
dynamic test system to determine power Flavio Alfredo Quizhpi.- Nació
electronics systems performance,” Revista de la en Cuenca, Ecuador, en 1969.
Facultad de Ingenieria - UCV, vol. 23, no. 3, pp. Recibió el Título de Ingeniero
91 – 102, 2008. Electrónico en la Universidad
[12] Viola J., Restrepo J., Quizhpi F., Gimenez M., Politécnica Salesiana en el 2003,
Aller J., Guzman V., Bueno A., “A flexible Licenciado en Ciencias de la
hardware platform for applications in power Educación en la Universidad
electronics research and education,” in Electrical Politécnica Salesiana en el
Power Energy Conference (EPEC), 2014 IEEE, 1999, Especialista en Educación
pp. 1–6, Nov. 2014. Superior en la Universidad del Azuay en el 2002.
Actualmente cursando el programa Doctoral en Ingeniería
Eléctrica en la Universidad Simón Bolívar de Venezuela.
Dedicado al área de investigación de los Sistemas
Marco Vinicio Fajardo.- Nació Eléctricos de Potencia, Confiabilidad, FACTs e Inversores
en Ambato-Ecuador en 1990. Multinivel
Recibió su título de Ingeniero
Electrónica de la Universidad
Politécnica Salesiana en 2013; de
Master en Ingeniería Biomédica
de la Universidad Politécnica
de Madrid, Madrid-España en
2015 Actualmente, se encuentra
cursando sus estudios de Doctorado en la Universidad
Simón Bolivar, Caraca-Venezuela y su campo de
investigación se encuentra relacionado con los Sistemas
de Potencia y Eficiencia Energética. Sus intereses de
investigación se centran en Sistemas Eléctricos de
Potencia, FACTs e Inversores Multinivel y eficiencia
energética.
71
Artículo Académico / Academic Paper
Abstract Resumen
This Document presents a Methodology to Este documento presenta una Metodología para
Identify the Areas of Voltage Control (ACV) in an Identificar las Áreas de Control de Voltaje (ACV)
electrical power system. The use of data mining en un sistema eléctrico de potencia. El uso de
and the concept of electric distance have been used minería de datos y el concepto de distancia eléctrica
in this work to identify the ACV in a power system han sido utilizados en este trabajo para identificar
for different probabilistic scenarios of system las ACV en un sistema de potencia para distintos
operation considering N-1 contingencies, which escenarios probabilísticos de operación del sistema
were formulated by Monte Carlo simulation. ACV considerando contingencias N-1, los cuales fueron
identification in power systems allows to identify formulados mediante simulación Monte Carlo.
reactive power resources of each area to establish La identificación de las ACV en un sistema de
operational strategies in order to improve the potencia permite identificar los recursos de
control of the voltage in the system. The proposed potencia reactiva de cada área y con ello establecer
methodology has been implemented to identify the estrategias operativas a fin de mejorar el control
ACV in the National Interconnected System of del voltaje en el sistema. La metodología planteada
Ecuador. ha sido implementada para identificar las ACV en
el Sistema Nacional Interconectado de Ecuador.
Index terms− Areas of Voltage Control, Electric Palabras clave− Áreas de Control de Voltaje,
Distance, Data mining, Monte Carlo simulation. Distancia Eléctrica, minería de datos, simulación
Monte Carlo.
72
Edición No 14, Enero 2018
73
Bonilla et al. / Metodología para Identificar Áreas de Control de Voltaje en un SEP Aplicando Simulación Monte Carlo
(1)
(5)
En ocasiones el método Monte Carlo ha sido
utilizado para obtener resultados estadísticos
2.3. Técnicas de Minería de Datos
(representados generalmente por su valor medio y
desviación estándar), sin embargo en este trabajo
El concepto de minería de datos ha sido creado
no se pretende usar la simulación Monte Carlo para
principalmente para el análisis de datos multivariante.
obtener información estadística. El objetivo del
Dentro de las técnicas utilizadas en este trabajo son
presente trabajo es usar esta técnica probabilística para
análisis de componentes principales y el método de
generar una gran cantidad de escenarios de operación
agrupación de datos Fuzzy C Meas (FCM).
mediante la variación aleatoria de la carga en función
de una distribución de probabilidad y la selección de
contingencias N-1 en líneas de transmisión, y de esta 2.3.1 Análisis de componentes principales
manera determinar la matriz de sensibilidad (para
cada escenario) a partir de la solución de flujos de El análisis de componentes principales es una
potencia. técnica de análisis multivariante, para reducir
las variables de un conjunto de datos, sin perder
la información proporcionada por las variables
2.2. Distancia Eléctrica
originales. Esta técnica permite transformar las
variables originales, por lo general correlacionadas
El concepto de distancia eléctrica involucra a
entre sí, en nuevas pocas variables (no correlacionadas)
la matriz dQ/dV, la cual forma parte de la matriz
conocidas como componentes principales (PCs). Las
74
Edición No 14, Enero 2018
componentes principales están ordenadas de tal forma de los grados de pertenencia de un dato específico en
que las primeras componentes retienen las máximas todos los clusters es igual a 1.
variaciones presentes en las variables originales.
Una generalización de las componentes principales FCM se basa en la minimización de la función J:
son las coordenadas principales. El cálculo de las
coordenadas principales de un conjunto de datos (X)
puede ser obtenida a partir de una matriz de distancias (8)
D (matriz cuadrada de dimensión n x n, en donde n
representa el número de elementos u observaciones
del conjunto de datos, que para el presente trabajo
En donde:
representa el número de barras o nodos del sistema
que son considerados en el análisis).
uij : es el grado de pertenencia o membrecía de la
A partir de una matriz de distancias se puede i-ésima observación al j-ésimo cluster.
construir la matriz de similitud Q (similar a la matriz
de covarianza) la cual representa la variabilidad n : es el número de datos.
entre elementos, y puede ser obtenida de la siguiente
manera [4] [6]: cl : es el número total de clusters.
La matriz de coordenadas principales, Yr, puede La función J se minimiza si cumple con las
ser obtenida a partir de los valores (Λr) y vectores ecuaciones (9) y (10):
propios (Vr) de Q, de la siguiente manera [4] [6]:
(9)
(7)
2.3.2 Método de agrupación de datos Fuzzy C Means En este trabajo se utilizó la técnica FCM para
agrupar las barras en ACV, sin embargo, dado que el
La gran cantidad de datos e información que se algoritmo de FCM es una herramienta de análisis de
dispone actualmente ha hecho necesario contar con datos, el mismo no considera la configuración natural
técnicas que permitan extraer información implícita de la red eléctrica, de hecho, los grupos de barras que
en enormes bases de datos. Para manejar una gran resultan de la aplicación de la FCM no corresponden
cantidad de datos surge las técnicas de agrupación de necesariamente a las áreas eléctricas reales, por lo que
datos, las cuales son técnicas analíticas diseñadas con es necesario añadir una rutina adicional al algoritmo
el objetivo de detectar patrones de comportamiento de agrupamiento que compruebe la conectividad de
consistentes entre diferentes variables de tal manera las barras.
que esta información pueda ser agrupada o clasificada
de acuerdo a algún comportamiento [7]. Dentro de En este sentido, con el objetivo de identificar la
las técnicas utilizadas para la agrupación de datos conectividad de las barras se utiliza un algoritmo que
se encuentra la técnica Fuzzy C Meas (FCM) cuyo emplea la teoría de Goderya [9]. A continuación se
algoritmo de clasificación se describe a continuación presenta un resumen del algoritmo.
[7] [8]:
Sea Ar(mxn) la matriz de incidencia de un área r
El algoritmo FCM asigna a cada dato, un valor del sistema de potencia, de dimensión n x m (m es el
de pertenencia (valores entre 0 y 1) a cada uno de los número de ramas o líneas de un área r y n es el número
clusters formados, entonces un valor pertenece a un de barras de la misma área r) entonces la matriz de
cluster específico si su valor de pertenencia es mayor conectividad Cr(nxn) puede ser definida como sigue:
al grado de pertenencia de los otros clusters. La suma
75
Bonilla et al. / Metodología para Identificar Áreas de Control de Voltaje en un SEP Aplicando Simulación Monte Carlo
(11)
(12)
(13)
76
Edición No 14, Enero 2018
5. CONCLUSIONES
77
Bonilla et al. / Metodología para Identificar Áreas de Control de Voltaje en un SEP Aplicando Simulación Monte Carlo
Considerando contingencias N-1, cuando se [3] [Link], J.C. Sabonnadiere, J.-[Link], and
formaron 6 Áreas de Control de voltaje, tanto en J.-P. Paul, ‘‘Structural analysis of the electrical
hidrología alta como en baja, estas se conformaron system: application to secondary voltage control
por las barras que se muestran en la Tabla 5. in France,’’ IEEE Trans Power Syst, vol4, no.2,
pp. 479 – 486, May 1989.
La barra de la central Sopladora, si bien
físicamente se encuentra cerca de la central Molino, [4] J. Cepeda, “Evaluación de la Vulnerabilidad del
eléctricamente está conectada a la subestación Sistema Eléctrico de Potencia en Tiempo Real
Milagro, es por ello que se justifica que en todos los usando Tecnología de Medición Sincrofasorial”,
escenarios simulados la barra de la central Sopladora Tesis de Grado para obtener el título de Doctor
y las barras de Milagro, San Bartolo y Machala en Ingeniería Eléctrica, Instituto de Energía
siempre se encontraban en una misma Área de Eléctrica, Universidad Nacional de San Juan, San
Control de Voltaje. Juan, Argentina, Diciembre 2013, ISBN: 978-
987-33-4328-5.
6. RECOMENDACIONES
[5]
[Link] and [Link], “PowerFactory
Application for Power System Analysis“, Agosto
Las Áreas de Control de Voltaje identificadas en
2014.
un sistema eléctrico de potencia permiten identificar
las reservas de potencia reactiva de cada Área y con
[6] D. Peña, Análisis de Datos Multivariantes,
ello se puede determinar las zonas con déficit de
Editorial McGraw-Hill, España, cap. 1 – 8, Enero
potencia reactiva que podrían causar problemas de
2002.
estabilidad de voltaje. Es por ello que se recomienda
para trabajos futuros diseñar una metodología
[7] J. Rojas Julio Chavarro, Ricardo Laverde
para determinar automáticamente las reservas de
“Técnicas de Lógica Difusa aplicadas a la minería
potencia reactiva en función de cada una de las ACV
de datos“, Universidad Tecnológica de Pereira.
identificadas en el S.N.I.
ISSN 0122-1701.
Considerando que la estabilidad de voltaje de
[8] K .Hammouda, and F. Karray, “A Comparative
un sistema eléctrico de potencia es un fenómeno
Study of Data Clustering Techniques”, University
esencialmente local, se recomienda utilizar las
of Waterloo, Ontario, Canada.
ACV identificadas en el S.N.I. para realizar estudios
eléctricos de estabilidad de voltaje. [9] F. Goderya, A. Metwally, and O. Mansour, “Fast
Detection and Identification of Islands in Power
Systems”, IEEE Transactions on Power Apparatus
and Systems, pp. 217 – 221, 1980.
78
Edición No 14, Enero 2018
79
Aplicación Práctica / Practical Issues
E-mail: [Link]@[Link]
2
Operador Nacional de Electricidad CENACE, Quito, Ecuador
E-mail: jcepeda@[Link]
Abstract Resumen
This work analyzes the implications of several Este trabajo analiza las implicaciones de varias
conceptual considerations, such as modeling consideraciones conceptuales, tales como
and integration characteristics, synchronous características de modelación e integración,
machine modeling and saturation modeling, tipo de modelación de la máquina sincrónica y
in the results of the transient stability analysis efecto de la modelación de la saturación, en los
facing the occurrence of single-phase faults resultados del análisis de estabilidad transitoria
followed by opening and reclosing at San ante la ocurrencia de fallas monófasicas seguidas
Rafael–El Inga 500 kV transmission line. For de apertura y recierre monofásico en la línea
this, the software DIgSILENT PowerFactory, de transmisión San Rafael–El Inga de 500 kV.
ATP and EMTP-RV are used to model and to Para esto, se usan los programas DIgSILENT
simulate the dynamic behavior of the system PowerFactory, ATP y EMTP-RV, para modelar
under different considerations. The results show y simular el comportamiento dinámico del
that the models of synchronous machines based sistema ante diferentes consideraciones. Los
on Park transform present greater problems resultados muestran que los modelos de máquinas
of stability under unbalanced perturbations or sincrónicas basados en la Transformada de Park
asymmetric conditions. Likewise, it has been presentan mayores problemas de estabilidad ante
also demonstrated how differences of modeling perturbaciones desbalanceadas o condiciones
and integration algorithms can lead to different asímetricas. Asimismo, se demuestra cómo
dynamic behaviors in the system, which should be las diferencias de modelación y algoritmos
adequately taken into account in the studies. de integración pueden provocar diferentes
comportamientos dinámicos en el sistema, lo cual
debe ser adecuadamente tomado en consideración
en los estudios.
Index terms− ATP, EMTP-RV, PowerFactory, Palabras clave− ATP, EMTP-RV, PowerFactory,
Park transform, phase domain, saturation, single- Transformada de Park, dominio de fase, saturación,
phase reclosing. recierre monofásico.
80
Edición No 14, Enero 2018
81
Dután y Cepeda / Consideraciones sobre Análisis de Estabilidad Transitoria en el ajuste del tiempo de recierre monofásico
(CCS) ante recierres monofásicos de la línea San un valor bajo de voltaje (Vind_A); el cual, debido a que
Rafael – El Inga de 500 kV, en un escenario crítico, el aire en el punto de falla está caliente e ionizado,
tomando en cuenta que la máquina sincrónica es puede mantener la circulación de la corriente de arco
representada con dos tipos de modelos matemáticos secundario. Sin embargo, después de un tiempo,
diferentes: i) modelo basado en Transformada de con la ayuda del reactor de neutro y factores como
Park; y, ii) modelo en el Dominio de Fase. Se modela la velocidad del viento y condiciones ambientales, la
y valida el AVR de acuerdo a registros reales y se corriente de arco secundario tiene una alta probabilidad
modela el PSS en base a documentos IEEE. Debido a de autoextinguirse[3]. El tiempo desde la apertura
la configuración de la conexión de la central CCS al de la fase hasta que el aire recupera totalmente su
Sistema Nacional Interconectado (SNI) del Ecuador, propiedad aislante puede ser estimado mediante (1)
se usa un modelo máquina–barra infinita, el cual se que está en función del nivel de voltaje de la línea de
adecúa apropiadamente para representar las diferentes transmisión [8].
configuraciones de generación, topología y demanda.
VL
Los resultados muestran que usando modelos de t = 10.5 + (1)
máquinas sincrónicas basados en Park implementado 34.5
en DIgSILENT y en el modelo SM59 de ATP se
requiere de tiempos muertos menores a 500 ms para donde t es el tiempo muerto en ciclos y VL es el voltaje
que las máquinas de CCS no pierdan la estabilidad, de la línea de transmisión en kV rms.
mientras que usando el modelo de máquina en el
Dominio de Fase implementado como modelo SM58 Usando para sistemas a 500 kV, el tiempo muerto
en ATP y el modelo de Park mejorado, implementado es de 26 ciclos si la frecuencia es 60 Hz, que equivale
en EMTP-RV, el tiempo muerto para el escenario a 430 ms, el cual es aproximado a 500 ms para este
bajo estudio puede ser superior a 1 s. Asimismo, nivel de voltaje.
se demuestra cómo las diferencias de modelación El origen del arco secundario es el acoplamiento
y algoritmos de integración pueden provocar capacitivo entre las fases sanas y la fase abierta; este
diferentes comportamientos dinámicos del modelo efecto puede ser neutralizado mediante el uso de un
del sistema, lo cual debe ser adecuadamente tomado arreglo de cuatro reactores (tres de fase y uno de
en consideración en los estudios. neutro) conectados a la línea de transmisión, de tal
manera que la corriente de arco secundario pueda,
2. MARCO TEÓRICO REFERENCIAL
mediante una sintonización perfecta del reactor de
2.1. Recierre Monofásico, Arco Secundario y neutro, llegar a ser cero [3]; en la práctica, éste valor
Estabilidad Transitoria no se consigue. En [9] se presenta una metodología
para la sintonización de reactores de neutro del
La corriente primaria de falla a través del aire sistema de EAV del Ecuador.
provoca su calentamiento y lo ioniza, perdiendo
su capacidad de aislamiento. Sobre la base de lo Para visualizar el efecto de una apertura y recierre
expuesto por [3], en la Fig. 1 se observa que una falla trifásico y de la apertura y recierre monofásico en la
fase-tierra (fase A) se la puede despejar mediante la estabilidad transitoria, es posible usar el criterio de
apertura de, únicamente, la fase fallada (interruptores áreas iguales. Si el esquema empleado para despejar
S1 y S2). la falla es el ATF-RTF, entonces, de acuerdo a la Fig.
2, el área A1 es mayor al área A2 lo que significa
que la máquina sincrónica pierde el sincronismo; en
VC este esquema, para evitar este evento, los tiempos
VB
B1 - B0
S1
VA de recierre trifásico deben ser muy pequeños y no
B1 - B0 siempre es posible aún con el uso de interruptores
B1 - B0
B0
ultrarrápidos. En cambio, en la Fig. 3 se observa que
B0 Vind_A
B0
en un esquema AMF-RMF el área A1 es menor que
b1 el área A2, lo que implica que la máquina no perderá
C
b1
b1 Corriente de sincronismo.
Arco Secundario
B
n’
S2
A bn
g
82
Cepeda, Echeverría y Argüello / Incorporación de Inteligencia en el S.N.T. Ecuatoriano Proyecto Yachak
83
Dután y Cepeda / Consideraciones sobre Análisis de Estabilidad Transitoria en el ajuste del tiempo de recierre monofásico
Voltaje de armadura
P.U.
CCS 4x205 MVA Corriente de armadura
Zeq
ING SRF AVR 1.6
Y Δ
S1_1 S1_2 PSS
LT_1
U1, U2, U3 y U4 1.4 3 1
Δ 8 km
118 km
Zeq A A
S2_1 F S2_2 1.2 2
LT_2 4x205 MVA
3x600 MVA
118 km 8 km AVR
Y Δ 1
PSS
U5, U6, U7 y U8
0.8
Zeq
0.6
Figura 4: Diagrama Unifilar del Sistema Bajo Estudio 0.4 1) Saturación-circuito-abierto (1092 A)
2) Saturación-corto-circuito (1062 A)
0.2 3) AGL (1049 A)
3.1. Escenario Crítico del Sistema Bajo Estudio Corriente de campo 1 P.U. = 1049 A
0 P.U.
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 Corriente
2.4
De acuerdo a la Fig. 4, la topología más crítica de campo
para analizar la estabilidad transitoria de la central Figura 5: Curva de Saturación de una Unidad de la Central
CCS resulta cuando se encuentra operando una sola CCS
LT entre las subestaciones San Rafael y El Inga y,
de acuerdo a estudios iniciales basados en flujos de Cada unidad de generación posee un AVR del tipo
potencia, transmitiendo 740 MW desde CCS; esto EXC9000 tal como se muestra en la Fig. 6 y un PSS
implica que, para el caso crítico de estabilidad, se del tipo IEEE PSS2B que se muestra en la Fig. 7.
debe considerar cuatro unidades conectadas. Bajo Vref
estas condiciones, el objetivo es analizar la estabilidad
Qt
transitoria de las cuatro unidades de la central CCS, 1
1+sTR
KIR +
+
Ufmax
representadas con diferentes modelos de máquina, Uf
Vt
ante diferentes tiempos muertos de recierre. 1
1+sTR
-
+
+
KR
1+sTA3
1+sTA4
1
1+sTS
PSS
Ufmin
El análisis se realiza mediante simulaciones Pt 1
1+sTR PSS2B
llevadas a cabo usando tres tipos de programas Según IEEE
Δω 421.5-2005
computacionales: PowerFactory de DIgSILENT, 1
1+sTR
cuyo modelo de máquina se basa en la Transformada
de Park [12], ATP el cual dispone de los modelos
Figura 6: Modelo Simplificado de AVR de cada Unidad de la
de máquina SM59 [5] (Transformada de Park) y
Central CCS
SM58 (modelo en el dominio de fase) y EMTP-RV
cuyo modelo de máquina se basa en una mejora al
UST max
algoritmo de Park. Ks1
1+sT1 1+sT3 1+sT5 PSS
1+sT2 1+sT4 1+sT6 UST min
En todos los casos, el estudio se fundamenta en un
ω N
modelo máquina–barra infinita; en PowerFactory el sTW1
1+sTW1
sTW2
1+sTW2
1
1+sTR
+
+
1+sT8
1+sT9 M
+
-
equivalente de red está ubicado en la barra de 500 kV
de la SE El Inga, mientras que en ATP y EMTP-RV, Ks3
84
Edición No 14, Enero 2018
3.3. Validación del Modelo de AVR según acuerdo al diagrama de la Fig. 4, el modelo se ajustó
Pruebas de Campo y Modelación en a las condiciones del sistema durante la prueba real
PowerFactory del AVR realizada en campo. La Fig. 9 presenta la
Durante la sintonización en campo de los PSS de la respuesta dinámica de la salida del AVR ante la prueba
central CCS se realizaron mediciones de la respuesta de pulso descrita tanto de los registros de campo como
del AVR de una unidad. La metodología consistió de las simulaciones realizadas en PowerFactory y
en inyectar un pulso de un 5% sobre el valor de ATP. Es posible apreciar la buena proximidad de las
referencia del AVR de una unidad en vacío, aplicado tres señales, lo que ha permitido validar los modelos
durante 1 s. Este evento de prueba se reprodujo en el tanto de PowerFactory como de ATP. Además, en
programa PowerFactory, lo que permitió realizar la caso de interés, un análisis más profundo respecto de
validación de la respuesta dinámica del AVR basado la validación del modelo del AVR de CCS puede ser
en la aceptable proximidad de su salida simulada encontrado en [16].
con los datos registrados en campo. En ATP, de
85
Dután y Cepeda / Consideraciones sobre Análisis de Estabilidad Transitoria en el ajuste del tiempo de recierre monofásico
4.1. Recierre Monofásico en ATP Modelo SM59 En la Fig. 12 se muestran oscilaciones amortiguadas
del ángulo mecánico. En la Fig. 13 se observa la
En las Fig. 10 y 11 se muestran el ángulo mecánico variación de la velocidad mecánica respecto a la
y la potencia de salida de una unidad de CCS, para un sincrónica. En la Fig. 14 se muestra el voltaje terminal,
tiempo muerto de 500 ms, donde se observa que las mientras que en la Fig. 15 se observa que la salida del
máquinas pierden estabilidad, mientras que para un AVR sigue, en oposición de fase, al voltaje terminal;
tiempo muerto de 300 ms se obtiene una estabilidad es decir, cuando el voltaje terminal disminuye la
oscilatoria. Sin embargo, de acuerdo a lo señalado salida del AVR se incrementa y viceversa. En las Fig.
en [6], éstos resultados pueden deberse al modelo de 16 y 17 se muestra la potencia activa y reactiva de
máquina basado en la Transformada de Park) que es una unidad de generación, mientras que de la Fig.
susceptible a perder estabilidad numérica si se toma 18 se observa la respuesta en oposición de fase del
en cuenta la curva de saturación y/o si el sistema PSS ante la variación de la potencia activa; es decir,
es no simétrico. Simulaciones adicionales fueron cuando la potencia de salida disminuye la salida del
realizadas tomando en cuenta lo sugerido en [7], para PSS se incrementa y viceversa. Además, de la Fig. 17
lo cual se disminuyó el paso de integración numérica se observa que durante la sincronización de la línea se
de la simulación obteniéndose resultados iguales a los producen cambios muy drásticos en la potencia activa
mostrados en las Fig. 10 y 11. de cada máquina sincrónica: la potencia cambia de
184 MW a 325 MW casi de manera instantánea; este
comportamiento debe ser analizado para determinar si
la máquina es capaz de soportar esta variación.
86
Edición No 14, Enero 2018
Figura 12: Ángulo Mecánico para Tiempos de Recierre de 500 Figura 16: Potencia Reactiva para Tiempos de Recierre de 500
ms y 1000 ms ms y 1000 ms
87
Dután y Cepeda / Consideraciones sobre Análisis de Estabilidad Transitoria en el ajuste del tiempo de recierre monofásico
DIgSILENT
infinita en donde el equivalente de red del SNI se
conecta en la barra de 500 kV de la SE El Inga. La
Fig. 20 presenta el diagrama unifilar del sistema en
estudio donde una de las líneas San Rafael–El Inga, a
500 kV, está fuera de servicio. Figura 22: Potencias Activa y Reactiva para Tiempo de
Recierre de 500 ms
XL_SRaf_RL8_500
XL_Inga_RL6_500
0.033
28.0
-0.0
0.0
0.0
0.441
0.000
0.000
-35.4
16.3
0.0
0.0
0.0
0.0
G_HPAS_CCS
491.4
0.98
-735.6
743.2
130.5
0.881
-37.6
SG
90.5
SNI
~
-735.6
366.8
0.431
0.433
739.9
0.871
0.887
-35.0
-93.6
15.7
41.5
13.8
1.00
0.034
0.033
28.9
28.0
-0.0
491.4
0.0
0.98
7.1
491.1
483.0
0.98
0.97
-0.0
6.9
XL_Inga_RL8_500
XL_SRaf_RL6_500
SanRafael
El_Inga
88
Edición No 14, Enero 2018
4.5. Discusión
En esta sección se han presentado una serie de
resultados de simulaciones referentes al mismo
fenómeno (análisis de tiempos de recierre monofásico
en la línea de transmisión San Rafael–El Inga de
500 kV cuando la línea en paralelo se encuentra
indisponible), considerando varios posibles aspectos
conceptuales de modelación y simulación, que se
ha demostrado, afectan los resultados del análisis.
Estas consideraciones conceptuales se resumen a
continuación:
• Tres programas de simulación diferentes:
Figura 24: Potencias Activa y Reactiva para Tiempo de ATP, EMTP-RV y PowerFactory. A este
Recierre de 500 ms, Sin Saturación
respecto, es importante mencionar que cada
software emplea diferente modelación y
Finalmente, el mismo análisis de recierre diferentes algoritmos de integración.
monofásico a los 500 ms de ocurrido el despeje
• Dos tipos de modelo de máquinas sincrónicas:
de la falla se realiza considerando la saturación
Modelo Basado en la Transformada de Park
exponencial en los generadores pero con una potencia
y Modelo en el dominio de fase.
de cortocircuito del sistema equivalente de 3200
MVA (un sistema ligeramente menos robusto). En • Dos consideraciones diferentes en la
este caso, la respuesta de la simulación indica que modelación del sistema equivalente (barra
la central CCS pierde sincronismo en la segunda infinita): sistema reducido o equivalente
oscilación, como se muestra en las Fig. 25 y 26. estático con potencia de cortocircuito.
• Inclusión o no del efecto de la saturación
magnética de la máquina sincrónica.
De los análisis realizados se ha logrado demostrar
lo siguiente:
• Las características propias de cada
programa de simulación (ATP, EMTP-RV
o PowerFactory), relacionadas tanto con la
forma de modelación como con los algoritmos
de integración, hacen que los resultados
sean diferentes. Así, las simulaciones en
ATP y EMTP-RV muestran un sistema con
déficit de torque sincronizante, en tanto
que las simulaciones en PowerFactory
muestran un sistema con falta de torque
de amortiguamiento. Estas características
Figura 25: Ángulo del Rotor para Tiempo de Recierre de 500
ms, 3200 MVA de Potencia de Cortocircuito del SNI
hacen que para ATP y EMTP-RV resulta
beneficioso no considerar la saturación de
las máquinas sincrónicas en tanto que para
PowerFactory esto resulta perjudicial.
• En el caso de PowerFactory, en donde se
ha modelado al SNI como un equivalente
estático, se aprecia claramente que la
conexión de CCS con un sistema más robusto
(con mayor potencia de cortocircuito) mejora
apreciablemente la condición de estabilidad.
• La elección del tipo de modelo de la máquina
sincrónica resulta crítico. A este respecto, se
ha verificado las afirmaciones hechas en [6]
respecto que el modelo basado en el marco
de referencia dq0 resulta inestable ante
perturbaciones desbalanceadas o condiciones
Figura 26: Potencias Activa y Reactiva para Tiempo de asimétricas y en [11] donde se indica que una
Recierre de 3200 MVA de Potencia de Cortocircuito del SNI mejora al algoritmo dq0 permite que este
89
Dután y Cepeda / Consideraciones sobre Análisis de Estabilidad Transitoria en el ajuste del tiempo de recierre monofásico
modelo funcione de una manera adecuada en [4] Ramnarayan Patel, T. S. Bhatti, and D. P. Kothari,
eventos donde la red está desbalanceada y se “Improvement of power system transient stability
toma en cuenta la saturación de la máquina. by coordinated operation of fast valving and
braking resistor - IEEE Xplore Document,” IEE
5. CONCLUSIONES Y RECOMENDACIONES Proc. - Gener. Transm. Distrib., vol. 150, no. 3,
pp. 311–316, May 2003.
En este artículo, se han puesto en evidencia
varias consideraciones conceptuales relacionadas
[5] H. W. Dommel, EMTP Theory Book. Microtran
con la condición de estabilidad transitoria respecto
Power System Analysis Corporation, 1996.
del ajuste del tiempo de recierre monofásico de la
línea de transmisión San Rafael–El Inga de 500
[6] E. Haginomori, T. Koshiduka, H. Ikeda, and J.
kV. En este sentido, se ha demostrado que aspectos
Arai, Power System Transient Analysis: Theory
básicos de modelación resultan críticos al momento
and Practice Using Simulation Programs (ATP-
de definir los tiempos críticos de recierre monofásico
EMTP). John Wiley & Sons, 2016.
del sistema. Estos aspectos son: el programa de
simulación (ATP, EMTP-RV o PowerFactory), el tipo
[7] IEEE, “Re-examination of Synchronous Machine
de modelo de la máquina sincrónica, la inclusión de
Modeling Techniques for Electromagnetic
saturación, entre otros posibles. Un aspecto a resaltar
Transient Simulations,” IEEE Transactions on
es que se ha confirmado el concepto expresado por [6]
Power Systems, vol. 22. pp. 1221–1230, Aug-
respecto de la inestabilidad del modelo basado en la
2007.
transformada de Park cuando se simulan condiciones
desbalanceadas.
[8] IEEE, “IEEE Guide for Automatic Reclosing
Sobre la base de lo mencionado, es oportuno of Circuit Breakers for AC Distribution and
destacar que un adecuado trabajo de modelación debe Transmission Lines,” IEEE Std C37104-2012
necesariamente estar acompañado de un proceso de Revis. IEEE Std C37104-2002, pp. 1–72, Jul.
validación de modelos. En este sentido, no es posible, 2012.
únicamente con resultados de simulación, definir
qué modelo deberá ser el apropiado para realizar las [9] W. Dután, L. Dután, “Sintonización de reactores
simulaciones, más aún cuando se pretende analizar de neutro en líneas de 500 kV del sistema
fenómenos de alta criticidad (como el recierre eléctrico ecuatoriano,” presentado en las XXVI
monofásico ante condiciones de alto estrés). A este Jornadas de Ingeniería Eléctrica y Electrónica -
respecto, mientras no se haya realizado el proceso de EPN, Quito, 2016.
validación de modelos, lo recomendable es optar por
aquel que brinde los resultados más conservativos a [10] E. W. Kimbark, Power System Stability. John
costa de una inevitable des-optimización operativa. Wiley & Sons, 1995.
Como trabajo futuro se ha planteado investigar [11] Ulas Karaagac, “Synchronous Machine Modeling
con mayor profundidad sobre técnicas de validación Precision and Efficiency in Electromagnetic
de modelos en base a mediciones, lo cual permitirá, Transients,” Ph.D. Thesis, École Polytechnique
en última instancia, definir apropiadamente cuál es el de Montréal, Montréal, 2011.
mejor marco de referencia para realizar estudios de
impacto crítico. [12] DIgSILENT PowerFactory, “Technical Reference
Documentation - Synchronous Machine.” 2011.
REFERENCIAS BIBLIOGRÁFICAS
[1] IEEE, “Single phase tripping and auto reclosing [13] P. M. Anderson and A. A. Fouad, Power System
of transmission lines-IEEE Committee Report,” Control and Stability, 2nd Ed. Wiley India Pvt.
IEEE Transactions on Power Delivery, vol. 7. Limited, 2008.
pp. 182–192, Jan-1992.
[14] P. Kundur, N. J. Balu, and M. G. Lauby, Power
[2] E. G. and A. C. Héctor J. Altuve, Normann system stability and control. McGraw-Hill, 1994.
Fischer, and Armando Guzmán, “Tutorial on
Single-Pole Tripping and Reclosing,” presented [15] L. Dubé, Users guide to MODELS in ATP, Apr-
at the 39th Annual Western Protective Relay 1996.
Conference, Spokane, Washington, 2012.
[16] J. Cepeda and P. Verdugo, “Automatic Voltage
[3] E. W. Kimbark, “Suppression of Ground-Fault Regulator Model Validation Based on Mean-
Arcs on Single-Pole-Switched EHV Lines by Variance Mapping Optimization and Field Tests,”
Shunt Reactors,” IEEE Trans. Power Appar. Rev. Téc. “energía” No 14, enero de 2018.
Syst., vol. 83, no. 3, pp. 285–290, Mar. 1964.
90
Edición No 14, Enero 2018
[17]
SIEMENS, “Adaptive Autoreclosure to Jaime Cepeda Campaña.-
Increase System Stability and Reduce Stress to Nació en Latacunga, Ecuador
Circuit Breakers,” presented at the The Nordic en 1981. Recibió el título
Workshop in Power System Protection and de Ingeniero Eléctrico en la
Control, Trondheim-Norway, 2017. Escuela Politécnica Nacional
en 2005, y el de Doctor en
Ingeniería Eléctrica en la
Universidad Nacional de San
Luis Manuel Dután Amay.-
Nació en la parroquia Guapán Juan en 2013. Entre 2005 y 2009
de la ciudad de Azogues, trabajó en Schlumberder y en el CONELEC. Colaboró
provincia de Cañar. Recibió el como investigador en el Instituto de Energía Eléctrica,
título de Ingeniero Eléctrico de Universidad Nacional de San Juan, Argentina y
la Escuela Politécnica Nacional en el Instituto de Sistemas Eléctricos de Potencia,
donde también ha cursado Universidad Duisburg-Essen, Alemania entre 2009
estudios de Maestría en y 2013. Su tesis de doctorado fue galardonada con
Ingeniería Eléctrica. Desde el el Primer Lugar del Premio “Domingo Faustino
2009 es colaborador de la Corporación Eléctrica del Sarmiento” a la mejor investigación realizada en la
Ecuador CELEC EP - Unidad de Negocio Transelectric provincia de San Juan, Argentina en 2014. Recibió el
en el departamento de Planificación y desde 2014 en premio Innovadores Ecuatorianos Menores de 35 años
la Subgerencia de Proyectos de Expansión – Área del MIT Technology Review en 2015. Actualmente
de Diseño de Subestaciones. Sus áreas de interés se se desempeña como Subgerente Nacional de
encuentran relacionadas con simulación y análisis de Investigación y Desarrollo del CENACE y como
transitorios electromagnéticos, equipo primario de Profesor a Tiempo Parcial en la Escuela Politécnica
subestaciones y coordinación de aislamiento de líneas Nacional. Sus áreas de interés incluyen la evaluación
de transmisión y subestaciones. de vulnerabilidad en tiempo real y el desarrollo de
Smart Grids.
91
Aplicación Práctica / Practical Issues
PETROAMAZONAS EP
E-mail: Ramiro_Sanguil@[Link]; Marcelo_Solis@[Link]
Abstract Resumen
Article 48 (Chapter No III) of “Reglamento para El Artículo 48, CAPÍTULO III del “Reglamento
el Libre Acceso a los Sistemas de Transmisión y para el Libre Acceso a los Sistemas de Transmisión
Distribución” published by ARCONEL (Agencia y Distribución” de la Agencia de Regulación y
de Regulación y Control de Electricidad), Control de Electricidad (ARCONEL) promueve
promotes the free access of third agents to existing el libre acceso de terceros agentes a la capacidad
or remaining transmission and distribution de transporte existente o remanente de los sistemas
capacity. In this context, the Oil Industry de transmisión y distribución. En este marco, la
Petroamazonas EP planed the expansion project empresa hidrocarburífera Petroamazonas EP
and the interconnection of its Electrical Power planifica el proyecto de expansión y la conexión
System (SEIP-E) to Ecuadorian Transmission de su Sistema Eléctrico Petrolero – Extendido
System (SNT) at 69 kV in the First Stage. (SEIP-E) al Sistema Nacional de Transmisión
(SNT) a nivel de 69 kV – Fase 0.
The principal results of electrical studies of the
interconnection between SNT and SEIP-E are En el presente documento se resumen los
presented in this document. Moreover, electrical principales resultados de los estudios eléctricos
studies in steady state and transient state are also en estado estable y dinámico de la interconexión.
included. In addition, this document presents the Además, se presenta la propuesta de interconexión
temporary interconnection proposal, which allows temporal que permitirá el intercambio de potencia/
the interchange of Power/Energy between the energía desde el SNT al SEIP-E.
SNT and SEIP-E.
92
Edición No 14, Enero 2018
Este proyecto de interconexión temporal permitirá, En general la interconexión del SEIP-E al SNT
desde el punto de vista técnico y económico, el presenta un buen desempeño técnico para los dos
desarrollo de la industria petrolera, abasteciendo sistemas eléctricos, tanto en condiciones normales
sus cargas con energía de tipo hidroeléctrico y de operación como bajo eventos de contingencia en
dinamizando la matriz productiva, mediante [2]: estado dinámico. En estado estable, la integración
de los sistemas permite mantener los voltajes en
1) Optimización del uso de recursos naturales valores admisibles de operación en las barras de las
no renovables para generación eléctrica subestaciones de los sistemas integrados, así como
(gas, petróleo). su factor de potencia en barra de interconexión,
de acuerdo con los nuevos límites sobre bandas de
2) Desplazamiento de generación a diésel, del variación de voltaje y factores de potencia del Sistema
cual el país es deficitario (reducir gastos). Nacional de Transmisión.
93
Sanguil y Solís./ Estudios Eléctricos de la Interconexión Temporal del SEIP-E de Petroamazonas EP al S.N.I. – fase 0
La interconexión se realizará por medio de cables 2.2. Modelo DIgSILENT Power Factory
de 69 kV, 244 mm2 doble circuito, con una longitud
aproximada de 4 km entre la S/E Shushufindi de 69 La propuesta de interconexión ha sido modelada
kV y la estructura de retención GV58. Para ello, un en el programa computacional Power Factory versión
circuito de la línea de transmisión Jivino – Shushufindi 2017 de DIgSILENT. Esta base de datos consolida el
de 230 kV de Transelectric EP será energizada a 69 SNT y SEIP-E, considerando: a) plan de expansión de
kV desde el patio de la S/E Jivino mediante una bahía las redes transmisión; b) escenarios de la demanda de
con esquema de barra principal y transferencia. los sistemas eléctricos; y, c) períodos hidrológicos del
Sistema Nacional Interconectado (S.N.I.).
Es importante destacar que, la estructura de
retención GV58 se encuentra aproximadamente Se modela el SEIP-E a nivel de 69kV formado
a 20 km de la S/E Jivino 230/69 kV de CELEC por: 12subestaciones de 69/13.8kV, 2 líneas de
EPTRANSELECTRIC. Las características de la subtransmisión enterradas simple circuito de 69 kV,
línea de interconexión entre el SEIP-E y el SNT se 10 líneas de subtransmisión aéreas simple circuito de
muestran en la Fig. 1. 69 kV y centrales térmicas de generación petróleo,
gas/petróleo, gas y diésel. La demanda del SEIP-E
En la Tabla 1 se presentan las características se caracteriza por su comportamiento prácticamente
eléctricas del cable de conexión del SEIP-E al plano en el tiempo, cuyo valor se encuentra alrededor
SNT, el cual es un cable con conductor de aluminio de 140 MW.
monofásico con aislamiento XLPE en configuración
enterrada triangular (trefoil).
3. METODOLOGÍA PARA ANÁLISIS DE
INTERCONEXIÓN
Tabla 1: Características Eléctricas del Cable de Interconexión
94
Edición No 14, Enero 2018
0,90
Factor de Potencia
0,85
4.1. Transferencia y Factor de Potencia 0,80
0,75
0,60
importadas por el SEIP-E desde el SNT, las cuales son Lluvioso Estiaje Lluvioso Estiaje Lluvioso Estiaje
medidas en la línea Nueva Jivino_69 kV - Estructura Año 2017 Año 2018 Año 2019
95
Sanguil y Solís./ Estudios Eléctricos de la Interconexión Temporal del SEIP-E de Petroamazonas EP al S.N.I. – fase 0
Voltajes (p.u.)
1,010
0,995
Dmax Dmed Dmin Dmax Dmed Dmin
0,990
Las contingencias 1 y 2 no afectan la integridad
0,980 de los sistemas debido a que, luego del despeje de las
0,970
0,976 fallas, los parámetros de frecuencia, voltaje y ángulo
de voltaje en barras de las subestaciones alcanzan una
0,960
Dmax Dmed Dmin Dmax Dmed Dmin nueva condición de equilibrio.
Lluvioso Estiaje
A continuación, para la contingencia 1, en la Fig.
8 se muestra como los parámetros de voltaje, ángulo
Figura 6: Voltajes en Barras del SEIP-E - Año 2017 CON
Interconexión
de voltaje frecuencia y derivada de la frecuencia en
barras de 69 kV del SNT. En la Fig. 9 se presenta la
respuesta de frecuencia de las unidades de generación
del SEIP-E, y en la Fig. 10 se muestra el ángulo de
rotor de las unidades de generación del SEIP-E. En
estas figuras se observa que las variables analizadas
alcanzan nuevas condiciones de equilibrio, la cuales
se encuentran dentro de los rangos aceptables de
operación poscontingencia.
96
Edición No 14, Enero 2018
Del análisis de las fig. 9 y 10 se observa que homologar Esquema de Alivio de Carga (EAC)
existe un error en la inicialización de las variables, instalado en el sistema petrolero, considerando:
este error es causado por la inicialización de los
modelos dinámicos del SVC Chinu y de las Centrales • Los ajustes de deslastre de carga del S.N.I.
Hidráulicas Chivor y Guavio, estos elementos
pertenecen al sistema colombiano. • La contingencia más severa para el SEIP-E
(pérdida de la interconexión del SEIP-E y
S.N.I.).
DIgSILENT
1.0404
0.8929
0.7453
0.5978
0.4502
0.3027
-0.2000 5.8399 11.880 17.920 23.960 [s] 30.000
B_Jivino_69: Voltage, Magnitude in p.u.
Jivino_69\B2: Voltage, Magnitude in p.u. • La frecuencia de operación de las unidades de
6.7560
5.2539
3.7517
generación del SEIP-E.
2.2496
0.7475
-0.7546
-0.2000 5.8399 11.880 17.920 23.960 [s] 30.000
B_Jivino_69: Voltage, Angle in deg
Jivino_69\B2: Voltage, Angle in deg • La máxima transferencia de energía desde el
60.154
60.107
60.059
S.N.I. hacia el SIEP-E.
60.012
59.964
59.917
0.6154
0.4692
Carga del SEIP-E, el cual es resultado del análisis de
las contingencias dinámicas en la interconexión del
0.3231
0.1769
0.0308
-0.1154
-0.2000 5.8399 11.880
B_Jivino_69: Derivative of El. Frequency in 1/s
Jivino_69\B2: Derivative of El. Frequency in 1/s
17.920 23.960 [s] 30.000
SEIP-E con el SNT.
Carga
Frecuencia Tiempo
Paso Campo total paso
[Hz] [s]
[MW]
Lago Agrio, CPF,
1 6,37 59,40 0,5
Sacha
CPF, Shushufindi,
2 6,92 59,20 0,6
ZPF
97
Sanguil y Solís./ Estudios Eléctricos de la Interconexión Temporal del SEIP-E de Petroamazonas EP al S.N.I. – fase 0
DIgSILENT
1,0393
0,8247
0,6100
0,3954
4,6503
SHUSHUFINDI_69: Voltage, Magnitude in p.u.
funciones gradiente de frecuencia ANSI 81R, bajo
-3,8448
-12,340 voltaje ANSI 27 y potencia inversa ANSI 32P.
-20,835
-29,330
-37,825
61,356
60,660 y ANSI 27, con los siguientes ajustes:
59,965
59,269
58,574
57,878
-0,1000 5,9199 11,940
SHUSHUFINDI_69: Electrical Frequency in Hz
17,960 23,980 [s] 30,000 • ANSI 81R:-3.20 Hz/s – 60 ms.
11,218
6,5044
1,7910
-2,9223 • ANSI 27: 0.8 p.u. – 500 ms.
-7,6356
-12,349
-0,1000 5,9199 11,940 17,960 23,980 [s] 30,000
DIgSILENT
1.00
0.90
de energía: 1) a nivel de 230 kV alimentado por los Francisco_de_Orellana_138\B1: Ángulo de Voltaje sin la desconexion del SEIP del SNI
Francisco_de_Orellana_138\B1: Ángulo de Voltaje con la desconexion del SEIP del SNI
San Rafael y Jivino (con una conexión fuerte Francisco_de_Orellana_138\B1: Derivada de la frecuencia con la desconexion del SEIP del SNI
11.078
el SEIP-E exporta potencia reactiva hacia SHS_JIV_69kV: Potencia Reactiva sin la desconexion del SEIP del SNI
SHS_JIV_69kV: Potencia Reactiva con la desconexion del SEIP del SNI
98
Edición No 14, Enero 2018
DIgSILENT
60.574
543.56
[Hz]
Imaginary Part [rad/s]
60.026
326.14
59.479
108.71
58.931
-1.00E+0.. -8.00E+0.. -6.00E+0.. -4.00E+0.. Real Part [1/s] -2.00E+0..
-108.71
58.384
-326.14
57.836
-0.1000 8.9199 17.940 26.960 35.980 [s] 45.000
CPF_69: Frecuencia con la desconexion del SEIP del SNI
SECOYA_69: Frecuencia con la desconexion del SEIP del SNI
ZPF_69: Frecuencia con la desconexion del SEIP del SNI
CPF_69: Frecuencia sin la desconexion del SEIP del SNI
SECOYA_69: Frecuencia sin la desconexion del SEIP del SNI
-543.56
ZPF_69: Frecuencia sin la desconexion del SEIP del SNI
Stable Eigenvalues
Unstable Eigenvalues
99
Sanguil y Solís./ Estudios Eléctricos de la Interconexión Temporal del SEIP-E de Petroamazonas EP al S.N.I. – fase 0
Por último, debido a que el SEIP-E no tiene [3] Agencia de Regulación y Control de Electricidad
participación en los factores de amortiguamiento (ARCONEL), “Regulación No. CONELEC
menores al 5%, no es necesario modificar los 004/01 Calidad del servicio eléctrico de
parámetros de ajuste de Reguladores de Voltaje distribución”. [Online] Available at: [Link]
(AVR) y Reguladores de Velocidad (Gobernadores) [Link]/regulaciones
de los generadores instalados en el SEIP-E.
[4] Agencia de Regulación y Control de Electricidad
(ARCONEL), “Regulación No. CONELEC
7. CONCLUSIONES Y RECOMENDACIONES
006/00 Procedimientos de Despacho y
Operación”. [Online] Available at: [Link]
Los resultados mostrados en este informe
[Link]/regulaciones
muestran que la alternativa de interconexión posee
un buen desempeño técnico en estado estable y
[5] SMARTPRO S.A., “Estudios eléctricos en estado
dinámico. Adicionalmente, es importante recalcar
estable de la interconexión del sistema eléctrico
que el proyecto de interconexión temporal permitirá
interconectado petrolero extendido (SEIP-E)
el desarrollo de la industria de explotación petrolera,
de Petroamazonas EP al sistema nacional
abasteciendo sus cargas con generación hidroeléctrica
interconectado (S.N.I.) – fase 0”, Consultoría de
de bajo costo, logrando mejorar su matiz energética y
Petroamazonas EP, Quito, Ecuador, Mayo, 2017.
dinamizar la matriz productiva.
[6] DIgSILENT Power Factory, “User’s Manual
Se demuestra que ante contingencias severas en el
Version 15”, DIgSILET GmbH, Gomaringen,
S.N.I. es necesaria la desconexión del SEIP-E. Esto
Germany, 2013
con la finalidad de evitar el colapso de los sistemas
interconectados. Para aislar el SEIP-E del SNI se
utilizaron las funciones de protección gradiente de
frecuencia y bajo voltaje, cuyos parámetros fueron
Ramiro Sanguil Andrade.-
analizados en los diferentes escenarios operativos,
Nació en Ambato en 1985.
hidrológicos y de demanda del SNI. Además, la
Recibió su título de Ingeniero
actualización del EAC del SEIP-E, permite mantener
Eléctrico de la Escuela
parámetros eléctricos dentro de los rangos admisibles
Politécnica Nacional en 2012.
de operación.
Desde 2015 trabaja en el
Departamento de Soluciones
Por último, es importante recomendar la
Energéticas de Petroamazonas
actualización o generación de nuevas regulaciones
EP. Actualmente cursa la
para análisis de interconexión de sistemas eléctricos.
Maestría en Redes Eléctricas Inteligentes en la
En estas regulaciones se deberá profundizar en el
Escuela Politécnica Nacional.
análisis dinámico, modos de oscilaciones, factores de
amortiguamiento, modelación de la carga, entre otros.
Marcelo Solís Sanchez.-
Nació en Ambato en 1984.
REFERENCIAS BIBLIOGRÁFICAS Obtuvo el título de Ingeniero
en Electrónica en la Escuela
[1] Agencia de Regulación y Control de Electricidad Politécnica del Ejército en
(ARCONEL), “Reglamento para el Libre 2009. Actualmente trabaja en
Acceso a los Sistemas de Transmisión y el Departamento de Soluciones
Distribución”. [Online] Available at: [Link] Energéticas de Petroamazonas
[Link] EP.
100
Aplicación Práctica / Practical Issues
E-mail: [Link]@[Link]
2
Universidad de Cuenca
E-mail: brian.jaramillo1@[Link]; jose.cardenas1007@[Link]
3
Universidad Estatal de São Paulo - UNESP
E-mail: maanmejia@[Link]
4
Universidad Federal ABC - UFABC
E-mail: padilha@[Link]; [Link]@[Link]
Abstract Resumen
Due to the incursion of new appliances with high- Debido a la incursión de nuevos electrodomésticos
energy consumption, the demand of the residential con alto consumo de energía, la demanda del
sector of an Electricity Distribution Utility can be sector residencial de una Empresa Eléctrica de
increased suddenly, affecting the operation of the Distribución puede ser incrementada de forma
electrical network and, consequently, changing repentina, afectando la operación de la red
eléctrica y, consecuentemente, cambiando los
the investment requirements for the expansion of
requerimientos de inversión para la expansión de
the electrical network. la red.
In this work was used a methodology composed En este trabajo se utilizó una metodología
of a spatial-temporal model and a method of compuesta por un modelo espacio temporal y un
impact analysis in the low voltage circuits of the método de análisis del impacto en los circuitos
electrical network. The model characterizes the de bajo voltaje de la red eléctrica. El modelo
intention of the inhabitants of the city to buy this caracteriza la intención de los habitantes de la
kind of appliances. On the other hand, the method ciudad en comprar esta clase de electrodomésticos.
determines the new values of operating, voltages Por otro lado, el método determina los nuevos
and the overload that can cause the supply of the valores de operación, voltajes y la sobrecarga que
demand for these appliances. puede provocar el suministro de la demanda para
estos electrodomésticos.
The results obtained provide valuable information
Los resultados obtenidos proporcionan información
on the current state of electrical installations of the valiosa del estado actual de las instalaciones de la
Electric Utility and the tendency of the impact that Empresa Eléctrica y la tendencia del impacto, que
can cause the supply of the demand of this type puede causar el suministro de demanda para este
of loads in the secondary circuits of the electrical tipo de cargas, en los circuitos secundarios del
network. sistema eléctrico.
Index terms− Distribution System Planning, Palabras clave− Planificación del Sistema de
Spatial-Temporal Estimation, Impact Studies, Distribución, Estimación Espacio Temporal,
Geographic Information Systems. Estudios de Impacto, Sistemas de Información
Geográfico.
En la actualidad, diversos tipos de cargas eléctricas El modelo espacio temporal está estructurado
como: termo tanques, aire acondicionado, duchas por dos módulos: uno espacial y uno temporal. En
eléctricas, cocinas de inducción, etc. son ofertados en [1] estos módulos fueron explicados, detallando los
el mercado con el objetivo de aumentar el confort en datos de entrada y su construcción en sistemas de
el hogar. La compra masiva de estos equipos puede información geográficos (GIS). De esta forma, en
afectar considerablemente la demanda de energía este artículo se mencionará los pasos principales del
eléctrica del sector residencial, y más aún, si el uso anterior modelo, para facilitar el entendimiento de la
de estos aparatos está siendo incentivado por políticas metodología seguida en este trabajo.
gubernamentales de un país [1].
2.1. Módulo Espacial
El gobierno de Ecuador creó un plan de incentivos
para disminuir el consumo de gas liquidificado El módulo espacial se fundamenta en la regresión
de petróleo (GLP). Una de las acciones de estos espacial GWR (Geographically Weighted Regression).
incentivos es el Programa de Eficiencia en la Cocción La regresión GWR permite modelar la distribución
(PEC) que motiva a la población el uso de cocinas espacial de la proporción de residencias con intención
de inducción eléctrica [2]. El PEC trae grandes de comprar cocinas de inducción para cada
desafíos a las Empresas de Distribución, desde el región de una área en estudio, con el fin de considerar:
punto de vista de la operación de la red eléctrica y los patrones geográficos, características socioeconómicas
proyectos futuros para su expansión. De esta forma, típicas y la influencia de las regiones vecinas en el
deben realizarse estudios detallados del impacto que proceso de estimación de [4]. De este modo, en
causará en la red, el suministro de la demanda de este módulo espacial se utiliza la expresión:
estas cocinas. Estos estudios permitirán determinar
las adecuaciones necesarias para satisfacer la nueva (1)
demanda.
Siendo: un conjunto de información socioeconómica
La compra de una cocina de inducción por parte representativa de la región es un parámetro de
de una familia depende de varios factores: precio de sensibilidad asociado a cada variable socioeconómica
la cocina, costo de instalación y de mantenimiento, utilizada y es la diferencia entre el valor estimado
etc. Por causa de estos factores, algunos habitantes por la regresión y el valor datos recolectados de
pueden estar poco interesados. Debido a esto, la proporción de residencias con intención de compra
se requiere un estudio que permita identificar la del producto Estos datos recolectados no
ubicación espacial de las residencias interesadas en consideran la influencia de la vecindad en la compra.
comprar esta clase de electrodoméstico y su demanda Esa influencia es considerada por el modelo a través
asociada. El resultado del pronóstico espacial de esta de los parámetros Los valores de pueden
nueva demanda, puede ser representado por medio de ser determinados por medio de sistemas GIS, por
mapas temáticos, para facilitar la visualización de la ejemplo, ArcGIS [5] o Terraview [6].
distribución espacial de la densidad de carga [3].
2.2. Módulo Temporal
En este artículo se utiliza una metodología
compuesta por un modelo espacio temporal y un Una regresión de distribución logística permite
método de análisis del impacto en las redes de bajo modelar el crecimiento a lo largo del tiempo de una
voltaje. El modelo espacio temporal determina el variable que está acotada por un valor finito [7].
número de domicilios que tienen la intención de De esta forma, el módulo temporal caracteriza la
compra de cocinas de inducción. Este modelo utiliza evolución temporal de la proporción de residencias
la información del censo demográfico, para realizar con intención de compra de la cocina de inducción
un estudio del potencial mercado, y así, estimar el por medio de la expresión:
porcentaje de residencias interesadas en la compra
del producto. Los resultados de este modelo, son
utilizados en el método de análisis del impacto en los (2)
circuitos secundarios, considerando como impacto, la
sobrecarga en las instalaciones eléctricas de la empresa
de distribución, el incremento de pérdidas técnicas Donde representa el número máximo de
y el incumplimiento de los niveles permisibles de residencias que pueden comprar la cocina; es
voltajes. Los resultados obtenidos proporcionan un la proporción de residencias con intención de compra
insumo importante para la planificación y operación en el primer año de estudio (k=0); t es el período
de la red de empresas de distribución de energía de estudio, número de años, y a es el coeficiente de
eléctrica. crecimiento anual de residencias con intención de
compra.
102
Zambrano et al. / Análisis Espacio Temporal del Impacto de Nuevos Tipos de Carga con Alto Consumo de Energía
2.3. Algoritmo para la Construcción del Modelo a causa del ingreso de nuevos servicios y a la gestión
comercial, es necesario realizar una fase de validación
El algoritmo del modelo espacio temporal puede de la información disponible, en base a inspecciones
ser construido en GIS, utilizando como información en campo y a depuración de la información en GIS,
de entrada los valores de para el primer año de forma de garantizar la conectividad eléctrica hasta
de estudio, el conjunto de variables socioeconómicas clientes y luminarias.
representativas X de la región , los parámetros
y a. La determinación de esos parámetros Adicionalmente, es necesario preparar una base de
se explican en [1]. El valor de puede ser parámetros eléctricos requeridos por la herramienta
obtenido a través de estudios de viabilidad económica. de flujo de potencia, esto es, parámetros físicos y
eléctricos de cables, conductores, transformadores y
Los pasos del modelo espacio temporal, se estructuras de soporte. Una característica de las redes
describen a continuación: de BV es el uso de transformador monofásico con
toma central, con ello, la herramienta de simulación
Paso 1. Determinar el crecimiento temporal debe ser capaz de modelarlo [8]. Este trabajo utilizó
por medio de (2). como herramienta de flujo de potencia el software
“Open Distribution System Simulator” (OpenDSS)
Paso 2. Calcular los valores de usando (1); [9].
Paso 3. Distribuir los valores de para cada La Fig. 1 muestra un ejemplo del resultado de
transformador dentro de cada región la construcción de los modelos de red, la imagen
superior corresponde a la información en GIS que
Después de haber calculado los valores de incluye una foto satelital y la imagen inferior al
para cada año del período de estudio, estos valores “plot” en el OpenDSS. El ejemplo corresponde a un
deben ser representados en mapas temáticos. Estos transformador de 75 kVA y su red asociada, donde se
mapas ayudan en la identificación de las instalaciones destaca que un punto de carga puede estar relacionado
eléctricas que precisan de un mayor análisis en la a uno o más clientes (medidores de energía).
evaluación del impacto del suministro de energía de
la nueva carga.
103
Edición No 14, Enero 2018
Perfil de carga diario del transformador MV/BV: para 3.3.2 Paso 2: Perfil de carga del transformador
ello se utiliza un perfil diario con una resolución o
intervalo de medición, generalmente, de 10 minutos. Este paso corresponde al barrido del perfil diario
de mediciones del transformador, para este trabajo se
Perfil de carga de la cocina de inducción: la Fig. 2 utilizó una resolución de 10 minutos.
presenta el perfil promedio de la cocina de inducción,
destacándose que el pico de demanda se encuentra al 3.3.3 Paso 3: Definición del Caso Base
mediodía [2].
En los pasos anteriores es necesario ajustar el
flujo de potencia a los datos de medición, a esto se
le denomina “Caso Base”, distribuyendo el valor
de la medición del transformador a cada una de las
cargas (clientes y luminarias). En el proceso de ajuste
se debe considerar las pérdidas en la red (líneas y
transformador).
104
Zambrano et al. / Análisis Espacio Temporal del Impacto de Nuevos Tipos de Carga con Alto Consumo de Energía
3.3.5 Paso 5: Simulaciones de flujos de carga en el Número de residencias donde hay más de 5
escenario “n” habitantes y por lo menos una persona permanente
en el hogar.
En cada escenario se realiza un barrido a todas
las cargas y se adiciona las potencias activas (P) y Estas variables fueron seleccionadas basándose en
reactivas (Q), conforme a los clientes seleccionados la experiencia de los planeadores. Adicionalmente,
en dicho escenario y el nivel de carga definido en el fue realizado un estudio de correlación estadística
paso 1. Posteriormente se simula un flujo de carga y usando el método de coeficiente de correlación de
se almacena los resultados similares al caso base. Spearman. Ese método fue escogido porque es más
adecuado para relaciones que son no lineales. Otros
En total para cada uno de los circuitos de BV se métodos pueden ser utilizados en función de los datos
simulan: 144 (hora: minuto del perfil diario) x [4 de entrada [1].
(Distribución de carga) + 7 escenarios] = 1584 flujos
de carga. Para que el software ArcGIS pueda realizar el
cálculo de GWR (expresión (1)), es necesario informar
4. RESULTADOS la estructura de vecindad entre las regiones del área
de estudio. De esta forma un archivo geo-referenciado
El área de estudio corresponde a la zona urbana de la ciudad Cuenca, fue cargado en el programa con
de la ciudad de Cuenca, Ecuador y el período las divisiones por sectores censitarios. Cada uno de
establecido para el análisis de crecimiento de carga es esos sectores representa una región donde se han
de 5 años. En esta sección se presentan los resultados recolectado los datos del censo.
obtenidos aplicando esta metodología. Las tensiones
se encuentran en 12,7 kV (MV) / 120-240 V (BV), 4.1.2 Datos para el módulo temporal
con un sistema multi-aterrado.
El módulo temporal requiere de los parámetros:
4.1. Resultados del Modelo Espacio Temporal
El modelo espacio temporal fue implementado en En este trabajo se asumió que el valor máximo de
el software ArcGIS 10.3 [5]. Este software contiene cocinas que podían ser instaladas, era igual
una librería llamada “Modeling Spatial Relationships” al total de residencias existentes en cada sector. El
de la caja de herramientas “Spatial Statistics Tools”, porcentaje de residencias en comprar las cocinas
que facilita el cálculo de los parámetros del módulo en el año inicial, , es determinado por un estudio de
temporal y de los valores del módulo espacial. viabilidad económica. El crecimiento de la instalación
de nuevas cocinas en cada sector, , fue calculado
4.1.1 Datos de entrada para el módulo espacial considerando el crecimiento natural de clientes
residenciales de la empresa de energía, y un registro
El módulo espacial requiere de las variables: histórico de instalación de cocinas en el sector.
. Adicionalmente, para el cálculo de este índice, fue
Los valores de representan el porcentaje considerado un factor que representa el interés de la
de residencias de cada región , que tienen condiciones población en comprar electrodomésticos, denominado
económica para comprar las cocinas, en el año k. Estos factor de innovación tecnológica.
valores no consideran la influencia de la vecindad, en
los resultados obtenidos. En el primer año k=0 (año 4.1.3 Crecimiento espacio temporal anual del
inicial), los fueron calculados a través de un porcentaje de domicilios con intención de comprar
estudio potencial de mercado. En los siguientes años cocinas de inducción
del período de estudio (k>0) el módulo temporal,
determinó estos valores, simulando el crecimiento de En las Fig. 4 y 5, se presentan los mapas
los resultados de . temáticos de la ciudad en estudio con los resultados
de la aplicación del modelo espacio temporal. Estos
105
Edición No 14, Enero 2018
106
Zambrano et al. / Análisis Espacio Temporal del Impacto de Nuevos Tipos de Carga con Alto Consumo de Energía
107
Edición No 14, Enero 2018
escenario 6 los transformadores alcanzan una carga Con respecto a los resultados del caso de estudio,
del 60%; y, para un escenario del 100% de cocinas bajo las condiciones actuales, el incremento de esta
la máxima carga promedio está alrededor del 147%. carga en los 5 años de predicción, es decir un 10%
de incursión de cocinas de inducción hasta el año 5,
Esta carga de alto consumo adiciona pérdidas no presenta impactos considerables en las redes de
técnicas a la red, para el caso Base, las pérdidas bajo voltaje de CENTROSUR, debido a que existe un
de energía son del 2,87% y durante el período nivel de reserva importante en los transformadores de
de predicción (escenario 6) no es significativo su distribución y en la capacidad de carga de las redes; sin
incremento, sino hasta el 100% de cocinas donde embargo en el caso de que exista algún acontecimiento
alcanza un valor de 3,5%. Este último valor equivale que modifique estas condiciones de incursión de
a un incremento de 0,63 puntos porcentuales en las cocinas, será necesario calibrar nuevamente el módulo
pérdidas totales del sistema (7,13%). temporal y analizar el rendimiento del sistema.
El proceso metodológico planteado, soportado en [8] Electric Power Research Institute EPRI, “Modeling
herramientas de planificación, es muy efectivo para Single-phase, Center-tapped Distribution
analizar los probables impactos en la red, producto de Transformers”, DSS Tech Note, Mar. 2009.
la incursión masiva de este tipo de cargas importantes.
El mismo puede ser utilizado para analizar el impacto [9] R. C. Dugan and T. E. Mcdermott, “An Open
en las redes de distribución en medio voltaje. Source Platform for Collaborating on Smart Grid
Adicionalmente, se puede analizar el impacto de Research,” in Power and Energy Society General
fuentes de energía distribuida y almacenamiento de Meeting, 2011 IEEE, pp. 1–7.
energía o su combinación con este tipo de cargas.
108
Zambrano et al. / Análisis Espacio Temporal del Impacto de Nuevos Tipos de Carga con Alto Consumo de Energía
109
Artículo Académico / Academic Paper
Abstract Resumen
This article proposes an algorithm for optimum Este artículo propone un algoritmo para
operation of isolated electrical systems, which operación óptima de sistemas eléctricos aislados,
minimizes the costs of producing electricity in a el cual minimiza los costos de producción de
mix of wind turbines and thermal generators. electricidad en un mix de aerogeneradores y
Intelligent recharging of electric vehicles and generados térmicos. Se incluye en el análisis la
energy storage systems is included in the analysis. recarga inteligente de vehículos eléctricos y de
To model the randomness of wind generation and sistemas de almacenamiento de energía. Para
electricity demand, a scenario tree is generated modelar la aleatoriedad de la generación eólica
by the Neural Gas method. The proposed y de la demanda de electricidad, es generado un
methodology is applied to a case study on the island árbol de escenarios por el método Neural Gas. La
of San Cristóbal in the Galapagos Archipelago metodología propuesta es aplicada a un caso de
as a contribution that encourages to reduce the estudio en la Isla San Cristóbal en el Archipiélago
high dependence of fossil fuels for the production de Galápagos como una contribución que incentive
of electricity and for the transportation of goods a reducir la alta dependencia de combustibles
and people and in the to reduce the economic and fósiles para la producción de electricidad y para
environmental problems of fuel transportation la transportación de bienes y personas, y en la
from the continent. búsqueda de reducir los problemas económicos y
ambientales del transporte de combustibles desde
el continente.
Index terms− Isolated power systems, renewable Palabras clave− Sistemas eléctricos aislados,
generation, electric vehicles, Galápagos Islands, generación renovable, vehículos eléctricos, Islas
stochastic program. Galápagos, programación estocástica.
110
Edición No 14, Enero 2018
111
Barreto et al. / Operación Óptima del Sistema Eléctrico Aislado de la Isla San Cristóbal