CD 10783
CD 10783
CD 10783
I
DECLARACIÓN DE AUTORÍA
Yo Vinicio Sebastián Ribadeneira Ramírez, declaro bajo juramento que el trabajo aquí
descrito es de mi autoría; que no ha sido previamente presentado para ningún grado o
calificación profesional; y, que he consultado las referencias bibliográficas que se incluyen
en este documento.
II
DEDICATORIA
Todos los actos nobles de mi vida se los dedico con todo el amor a mi hija Cayetana, mi
razón de existir, quien a su corta edad supo entender que debía finalizar un desafío
personal para cumplir una promesa. El sacrificio de años y muchas madrugadas se
plasman en este documento.
III
AGRADECIMIENTO
Por la gestión oportuna para solventar dificultades al Ing. MSc. Carlos Herrera
Por la información proporcionada a Ing. MSc. Antonio Ortiz EEASA. Ing. Luis Marcial
EEASA. Ing. Diego Andrade EEQ. Ing. Julio Esparza ELEPCO
IV
ÍNDICE DE CONTENIDO
AVAL.......................................................................................................................I
DECLARACIÓN DE AUTORÍA .............................................................................. II
DEDICATORIA ..................................................................................................... III
AGRADECIMIENTO ............................................................................................ IV
ÍNDICE DE CONTENIDO...................................................................................... V
RESUMEN .......................................................................................................... VII
ABSTRACT ........................................................................................................ VIII
1. INTRODUCCIÓN ............................................................................................... 1
1.1 Pregunta de investigación ......................................................................... 8
1.2 Objetivo General ....................................................................................... 8
1.3 Objetivos Específicos ................................................................................ 8
1.4 Hipótesis o alcance (de ser el caso) .......................................................... 9
1.5 Marco Teórico ......................................................................................... 10
1.5.1. Aspectos a ser considerados en Generación Distribuida. .................... 10
1.5.1.1. Función de la Generación Distribuida en la evolución de las redes ......
de distribución eléctrica. .................................................................. 12
1.5.1.2. Consecuencia directa de la Generación Distribuida en las redes .........
de distribución eléctrica. .................................................................. 14
1.5.2. Tipos de tecnologías aplicadas para Generación Distribuida. .............. 15
1.5.3. Funciones de probabilidad bivariada. ................................................... 18
2.METODOLOGÍA ............................................................................................... 25
2.1. Descripción matemática de una red eléctrica simple para comprobar la regla
de 2/3 ................................................................................................................... 25
2.2.Uso de OPENDSS para cálculo de pérdidas de potencia. .............................. 32
2.2.1.Descripción general..................................................................................... 32
2.2.2.Principales elementos de un SED a ser modelados en OPENDSS. ............ 35
2.2.3.Modelación en OPENDSS de la red de prueba de 123 nodos de la IEEE. .. 41
2.3.Uso del software Cyme en el modelamiento de redes de distribución de las......
empresas eléctricas en Ecuador ......................................................................... .45
2.4.Utilización de Matlab para implementación de la metodología. ...................... 55
2.4.1.Modelamiento de elementos eléctricos para OPENDSS. ............................ 55
2.4.2Interface de comunicación de Matlab con OPENDSS. ................................. 71
2.4.3.Implementación del método probabilístico bivariado. .................................. 71
2.4.4.Descripción del programa base. .................................................................. 75
V
2.4.5.Modificación del programa base para aproximación – ajuste de resultados ....
............................................................................................................................. 84
2.4.6.Programa base para análisis de energía. .................................................... 85
3.RESULTADOS Y DISCUSIÓN .......................................................................... 88
3.1.Resultados ..................................................................................................... 88
3.1.1Redes de prueba IEEE. ................................................................................ 88
3.1.2Redes de la Empresa Eléctrica Quito EEQ. ............................................... 122
3.1.3Redes de la Empresa Eléctrica Ambato EEASA. ....................................... 133
3.1.4Redes de la Empresa Eléctrica Cotopaxi ELEPCO. ................................... 143
3.2.Discusión ..................................................................................................... 154
3.2.1.Número de casos aleatorios necesarios. ................................................... 154
3.2.2.Método probabilístico bivariado versus trivariado. ..................................... 156
3.2.3.Análisis para la selección de la mejor opción en GD. ................................ 158
3.2.4.Comparación de resultados con CYME. .................................................... 164
3.2.5.Comparación del mejoramiento de red de prueba de 37 Nodos IEEE con ......
otro método. ....................................................................................................... 167
4.CONCLUSIONES ........................................................................................... 170
5.REFERENCIAS BIBLIOGRÁFICAS ................................................................ 174
ANEXO A. Ejemplos de aplicación de las funciones de probabilidad Rayleigh y ......
Weibull para cálculo de disponibilidad de energía eólica .................................... 179
ANEXO B. Análisis eléctricos factibles con OPENDSS ...................................... 182
ANEXO C. Comparación de resultados de CYME con IEEE para la red de prueba .
de 123 nodos. .................................................................................................... 183
ANEXO D. Estructura del programa base para determinar la mejor ubicación y ......
tamaño de Capacitores – Reactores en GD con la metodología probabilística ........
bivariada. ........................................................................................................... 193
ANEXO E. Configuración para conexiones de Generación Distribuida. ............. 198
VI
RESUMEN
La Generación Distribuida reduce las pérdidas de energía en las redes eléctricas y mejora
los perfiles de voltaje. Ya que resulta extremadamente complejo describir en un sistema de
ecuaciones las pérdidas de energía de una red de distribución sin ninguna simplificación.
Se utiliza el software OPENDSS para modelar la red y a través de Matlab se envía las
instrucciones de control y análisis de resultados. Para replicar la mayor parte de escenarios
posibles, se utiliza la probabilidad bivariada uniforme como una herramienta matemática
que reproduzca cada lugar y potencia de GD para seleccionar la mejor opción que reduzca
las pérdidas, considerando un análisis adicional con perfiles de energía. El empleo de
capacitores se lo realiza junto con GD para el mejoramiento de perfiles de voltaje. La
validación de resultados se lo hace comparando los niveles de precisión de CYME con
OPENDSS en la red de 123 nodos de la IEEE. La utilización de generación eólica y
fotovoltaica como recurso energético renovable es lo más utilizado en GD, como un primer
paso a la migración de redes inteligentes.
VII
ABSTRACT
Distributed generation reduces energy loss in power supplies and improves voltage profiles.
As it is extremely complex, describe a system of energy loss equations in a distribution
feeder without any simplification. It uses the OPENDSS software to shape the feeder and,
through Matlab, send control instructions and analysis of results. To replicate the best
possible scenario, it uses the bivariate probability uniform, as the mathematical tool that
reproduces each location and the power of DG to choose the best option that reduces loss,
considering additional analysis with energy profiles. The use of capacitors is done in
conjunction with GD to improve the voltage profiles. Validation of the results is done by
comparing CYME precision levels with OPENDSS in the 123 node IEEE feeder. The use of
wind and photovoltaic generation, as the renewable energy resource, is the most widely
used in GD, as the first step to migrate smart grids.
VIII
1. INTRODUCCIÓN
El constante aumento de la demanda eléctrica, hace que las redes de distribución operen
a ciertas horas a su máxima capacidad o sobrecargadas, lo que implica el aumento de las
pérdidas técnicas por efecto Joule. Para solucionar este problema, las empresas eléctricas
distribuidoras tradicionalmente optan por la ampliación o reforzamiento de las redes de
distribución del sistema eléctrico, lo que conlleva una mayor inversión económica en
infraestructura, para ampliar la capacidad de transmisión y distribución de energía y poder
satisfacer las necesidades de los clientes con mejores parámetros de calidad [1]. De
continuar dicha situación, se puede esperar que las necesidades energéticas de los
consumidores no puedan ser satisfechas oportunamente, siendo uno de los impedimentos
el tiempo de construcción de nuevas redes de distribución y el costo económico de las
mismas.
Otra alternativa puede ser la inserción de Generación Distribuida (GD) de baja potencia.
Las fuentes de GD, también conocida como generación in-situ, generación embebida,
generación descentralizada, generación dispersa o energía dispersa [1] consiste en la
conexión de pequeñas fuentes de energía a los lugares más próximos a las cargas. Entre
las principales características de la inserción de GD a sistemas de distribución eléctrica
reales consta:
• La GD puede ser del tipo renovable y no renovable con el beneficio adicional de que al
utilizar fuentes de energía renovables se disminuye la emisión de CO2.
• En el caso de una fuente energética por radiación solar como son los paneles solares,
la energía estará presente en un determinado periodo del día; es decir, que durante las
1
primeras horas de la mañana tendrá una baja producción energética e irá aumentando
hasta que la máxima producción de energía ocurre al mediodía, hora en la que existe
un pico de radiación solar, luego disminuye paulatinamente hasta un valor cero al
finalizar la tarde cuya gráfica Potencia vs. Tiempo es similar a una “Campana de Gauss”
[1], esto sin contar con la eventualidad que en ciertas ocasiones habrán nubes que se
interpongan entre los rayos solares y los paneles solares.
Ya que la GD es una opción para disminuir las pérdidas de potencia por efecto Joule; a
más de mejorar el perfil de voltaje y liberar la capacidad de la red. Para esto, se requiere
2
desarrollar una metodología práctica y efectiva. En publicaciones científicas han sido
desarrolladas varias metodologías para GD que van desde métodos analíticos con varias
simplificaciones [3] hasta cálculos más complejas como lo indicado en [4], [5] y [6].
Para determinar que metodología es la más conveniente para la reducción de pérdidas con
GD, se ha revisado publicaciones científicas referente a GD en las que consta las
metodologías numéricas como:
3
[36] utiliza BT como una variante a los métodos analíticos con varias simplificaciones
manifestada en [33]. La BT resuelve simultáneamente UOGD y ubicación óptima de fuentes
de potencia reactiva (capacitores) como lo indicado en [37]. [38] utiliza un modelo
estocástico continuo para determinar la ubicación óptima de generación distribuida. Esto
es resuelto con una combinación de AG, BT y Búsqueda Discreta (BD). La Optimización
por Enjambre de Partículas (OEP) es aplicado para resolver un modelo de UOGD con factor
de potencia considerando la potencia variable en cargas como lo indicado en [39]. El
mejoramiento de OEP con varios tipos de GD y potencia reactiva se indica en [40]. También
hay un modelo hibrido de AG y OEP planteado en [41]. La Optimización por Colonia de
Hormigas (OCH) para UOGD mediante un algoritmo heurístico es propuesta en [42];
mientras que la Colonia Artificial de Abejas (CAA) es otro método propuesto en [43]. La
Evolución Diferencial (ED) es utilizada para los cálculos de sensibilidad de voltaje en
ubicación y tamaño de GD como lo indica [44]. [45] utiliza otro tipo de algoritmo heurístico
para evaluar- la energía no suministrada vs. las pérdidas técnicas y determinar un lugar
aproximado de GD. Un método heurístico basado en aproximaciones de valores determina
la UOGD para minimizar los costos del sistema de confiabilidad es desarrollado en [46]. Un
desarrollo más complejo de la evolución diferencial se resuelve por un análisis de
sensibilidad de voltaje y análisis de sensibilidad de pérdidas mediante ecuaciones de flujo
de potencia con restricciones como lo indicado en [47]. Un método heurístico especial que
calcula las regiones de alta probabilidad para localización de GD es desarrollado en [48].
Finalmente hay métodos heurísticos interactivos para resolver UOGD en dos etapas que
son las técnicas de agrupamiento y búsqueda exhaustiva como lo indica en [49] o para
pequeñas redes de distribución como se presenta en [50].
Consideraciones generales:
Por lo tanto, del análisis a las metodologías presentadas en las citas bibliográficas
anteriores, se plantea que los métodos probabilísticos, como los utilizados en [13], [14],
[20], [35], [38] y [51] son los más factibles de ser aplicados para un mejoramiento de las
redes de distribución eléctrica por medio de GD por su facilidad de aplicación ya que
permiten:
• Que el modelo de carga sea diverso y con potencia variable en el tiempo acorde a la
topología de una red de distribución real.
• Que se realicen de forma simultánea análisis de las principales variables que forman
parte de una GD (posición, tamaño, número de GD y/o factor de potencia).
4
• Que el método aplicado determine la solución aproximada de acuerdo al objetivo de
estudio; así, por ejemplo, reducir las pérdidas de potencia, reducir costos de inversión
de GD, aumentar la capacidad de penetración de GD, aumentar el costo/ beneficio de
GD, mejorar los perfiles de voltaje, etc.
Las situaciones que debe considerar la metodología a ser aplicada consisten en: Analizar
el comportamiento del sistema eléctrico a ser intervenido a través de un flujo de potencia,
establecer si existe la necesidad y conveniencia de instalar GD en dicha red de acuerdo al
nivel de pérdidas de potencia y caídas de voltaje; además se debe considerar que la
mayoría de alimentadores primarios son: de gran longitud, con diversidad de tipos de
cargas y caídas de voltaje significativas en los tramos finales. Una vez analizada dichos
resultados y condiciones respectivamente, se establece el límite máximo de GD que podría
ser instalado de acuerdo a la capacidad de la red y cargas instaladas.
También se deben considerar las características propias del perfil generación eléctrica
según el tipo de fuente energética a ser utilizada.
5
Para poder determinar la mejor ubicación y capacidad de GD con recursos eólicos y
fotovoltaicos, se debe tomar en cuenta las siguientes características de este tipo de
energía: limitada capacidad de potencia, bajo factor de planta, la generación eólica y
fotovoltaica pueden producir mayor energía que la demanda en un intervalo de tiempo; y,
los lugares físicos para la instalación de equipos de generación son limitados por la
infraestructura y tamaño.
Con este trabajo, se pretende brindar una solución aproximada respecto a una solución
ideal o exacta para determinar el mejor lugar y capacidad de GD para la conexión a la red
de distribución en análisis, de acuerdo al perfil de carga de los transformadores de
distribución de la red, topología del sistema eléctrico y las características de las fuentes de
energía solar y eólica. Como beneficio para la red de distribución se obtendrá: la reducción
de pérdidas técnicas por efecto Joule, mejora de los perfiles de voltaje y aumento de la
capacidad de transmisión de energía.
Las simulaciones numéricas mediante la teoría de probabilidades son una buena técnica
para resolver problemas físicos y matemáticos de gran complejidad, siendo está opción
ampliamente utilizada en lugar de métodos analíticos de difícil solución.
En este trabajo por medio del uso de un método probabilístico, pretende determinar (con
resultados aproximados) la mejor ubicación y capacidad GD, con el objetivo de reducir las
pérdidas técnicas en la red de distribución. Las redes de distribución deben considerar
prioritariamente las siguientes variables dentro del análisis:
6
• Modelamiento de cargas concentradas o distribuidas a nivel de media tensión.
De las soluciones obtenidas para cada caso, se puede realizar un nuevo análisis de
simulación numérica en la “vecindad” de los mejores resultados para aumentar la precisión
del resultado final.
7
1.1 Pregunta de investigación
¿Cómo determinar la mejor ubicación y capacidad de Generación Distribuida para la
reducción de pérdidas técnicas por efecto Joule en un alimentador primario?
• Utilizar datos de configuración de la red a ser modelada; así como perfiles de carga y
disponibilidad del recurso energético de las fuentes de generación eólica o fotovoltaica
a ser implementada, todo esto para la programación de flujos de potencia en OpenDSS
e identificar las pérdidas del sistema sin GD.
8
1.4 Hipótesis o alcance
El alcance del presente trabajo es plantear una metodología probabilística para la
reducción de pérdidas en un sistema de distribución a nivel de alimentadores con la
inserción de GD, para lo cual el desarrollo incluye:
• Se planteará una metodología probabilística que aplique una simulación numérica para
la generación de casos aleatorios combinados de ubicación – capacidad para GD. Se
establecerá de forma práctica el número mínimo de casos aleatorios que deberán ser
generados con el suficiente nivel probabilístico de confianza para el análisis de cada
caso según el tamaño de la red.
• Se programará en Matlab una interface de comunicación con OpenDSS para que este
último realice el cálculo de flujos de potencia y determinación de pérdidas por efecto
Joule de cada uno de los casos combinados de GD, producto de la simulación numérica
obtenida de la distribución bivariada que mejor se ajuste a las características de la red
a ser analizada y mejorada.
9
ya que está tesis es una propuesta metodológica probabilística y por lo tanto
experimental.
• La inserción de GD en una red eléctrica causa que el flujo de potencia sea redefinido,
por lo que las protecciones eléctricas deben ser reajustadas a la nueva configuración
propuesta [61].
Según [62] los recursos distribuidos pueden emplearse según la oferta y demanda
energética del sistema de distribución (situación que no ocurre en un sistema de
transmisión). La fuente de generación puede estar instalada del lado del consumidor,
ocasionando que en ciertas condiciones se reduzca los niveles de demanda máxima,
incrementar la eficiencia energética de las industrias y la confiabilidad de la red de
distribución.
10
La evaluación del impacto de la inserción de GD en la operación de los sistemas de
distribución debe considerar aspectos como: operación en estado estable, coordinación de
protecciones, comportamiento dinámico del sistema y prestación de servicios auxiliares de
las unidades de generación. Estos impactos pueden ser analizados para comparar las
soluciones obtenidas en cada escenario antes y después de la introducción de GD.
Entre los desafíos que debe superar la integración de los recursos energéticos para una
GD eficiente consta [63]:
• Diseño de esquemas para una apropiada gestión de la demanda permitiendo que la red
se adapte a las necesidades de los clientes.
11
• Desarrollar estrategias de mercado que permitan la migración a las redes inteligentes
para la conexión de diferentes tipos de cargas en una integración perfecta a lo largo del
tiempo.
Figura 1. Estructura general de una instalación de Generación Distribuida en una red eléctrica
[64].
12
Figura 2. Configuración tradicional de un sistema de distribución eléctrica [65].
Figura 3. Inserción de generación distribuida en una red eléctrica para medio y bajo voltaje [65].
13
Figura 4. Inserción de generación distribuida en una red eléctrica para medio y bajo voltaje [65].
Las pérdidas de energía por efecto Joule son directamente proporcionales a la resistencia
del conductor y el cuadrado de la corriente. Al instalar GD cerca de las cargas de una red
eléctrica se reconfigura el flujo de potencia y por lo tanto se reduce la magnitud de la
corriente total que circula desde la principal fuente que es la subestación reductora, de tal
forma que las pérdidas por efecto Joule también se reducen. Las pérdidas en sistemas de
distribución se evalúan en función de un periodo de tiempo; es decir en MWh por año,
comenzando con las redes sin GD.
14
Posteriores análisis se realizan con diferentes configuraciones de los generadores de
acuerdo a: niveles de potencia disponible, ubicación en la red de distribución, número de
unidades, capacidad de control, factor de planta, etc. También se debe tomar en cuenta
los diferentes tipos cargas y demandas energéticas, las cuales se deben relacionar y
modelar matemáticamente mediante curvas tipo. Sin embargo, dependiendo de la
configuración de GD sujeta a prueba, se obtienen diferentes resultados de reducción de
pérdidas y en algunos casos aumenta según las pruebas publicadas en [66] como lo
muestra la Figura 5.
Según la Agencia Internacional de Energía [64] las fuentes de energía renovables son
definidos como los recursos que generalmente no están sujetos a agotamiento por su
abundancia en la naturaleza como: el calor o la luz solar, la fuerza del viento, la materia
orgánica (biomasa), caídas de agua, energía oceánica y calentamiento geotérmico. Las
tecnologías desarrolladas para estas fuentes energéticas son las unidades micro
hidráulicas, celdas fotovoltaicas, turbinas de viento, sistemas térmicos solares, baterías de
almacenamiento, las turbinas de gas de ciclo combinado y ciertos motores de combustión
interna que utilizan bio combustibles. Los rangos de potencia instalada de las tecnologías
de generación eléctrica se muestran en la Tabla 1.
15
Para la inserción de GD en redes de distribución eléctrica se debe considerar los costos de
inversión de la tecnología a ser aplicada versus la eficiencia de cada sistema y la
disponibilidad del recurso energético. De acuerdo a lo desarrollado en [68], la Figura 6,
muestra algunas de las tecnologías de generación eléctrica.
De los recursos energéticos renovables, la energía hidroeléctrica y eólica son los de mayor
utilización. Así, por ejemplo, existen países donde se instalan micro centrales de
generación hidroeléctricas a nivel de distribución para abastecer la demanda energética de
cargas especificas en ubicaciones lejanas; sin embargo, la hidroelectricidad está limitada
a la disponibilidad de dicho recurso. Ante ello varios países han optado por el desarrollo,
implementación y explotación de tecnologías renovables no convencionales como lo indica
la Figura 7.
16
Los países de la Organización para la Cooperación y Desarrollo Económico (OCDE) han
invertido más recursos en energías renovables no convencionales cuyos resultados
tienden a la energía eólica, solar y biocombustibles como lo indica la Figura 8.
17
Figura 10. Tasa de crecimiento anual mundial de generación solar fotovoltaica [71].
Figura 11. Expresión gráfica de una variable aleatoria en teoría de probabilidades [72] .
18
Figura 12. Expresión gráfica de la medibilidad de una variable aleatoria [72].
𝜒 𝜒1 𝜒2 𝜒3 …
𝑓(𝜒) Ρ1 Ρ2 Ρ3 …
Figura 13. Ejemplos gráficos de variable aleatoria discreta y continua con 𝑓(𝜒) [72].
Donde 𝑓(𝜒) se conoce como función de densidad o función de probabilidad; mientras que
19
Aplicando la definición de función de distribución:
La función de distribución será:
Figura 14. Ejemplo de función de distribución con variables aleatorias discretas [72].
Para el caso de una función de densidad con v.a continua dada por:
20
Figura 15. Ejemplo de función de distribución con variable aleatoria continua [72].
Existen una gran cantidad de funciones de distribución de tipo discreto o continuo. Los
nombres que poseen son por las situaciones en las que surgieron, el nombre del
descubridor del método o la persona que lo utilizo en alguna aplicación especial. El objetivo
principal de utilizar un determinado tipo de función de distribución, es que describa de forma
correcta el evento dentro del espacio de probabilidad (Ω, ℱ, Ρ) de acuerdo al espacio
muestral Ω obtenido. Entre las principales funciones que se menciona en [72] consta:
uniforme discreta, Bernoulli, binomial, geométrica, binomial negativa, hipergeométrica,
Poisson, uniforme continua, exponencial, gama, beta, Weibull, normal, ji–cuadrada, t, F.
1 𝓍 ∈ {𝓍1 , 𝓍2 , … . . , 𝓍𝑛 }
𝑠𝑖
𝑛𝑚 𝓎 ∈ {𝓎1 , 𝓎2 , … . . , 𝓎𝑚 }
𝑓(𝓍, 𝓎) =
{0 𝑒𝑛 𝑜𝑡𝑟𝑜 𝑐𝑎𝑠𝑜
1
𝑠𝑖 𝓍 ∈ {𝓍1 , 𝓍2 , … . . , 𝓍𝑛 }
𝑓𝑋 (𝓍) = { 𝑛
0 𝑒𝑛 𝑜𝑡𝑟𝑜 𝑐𝑎𝑠𝑜
1
𝑓𝑌 (𝓎) = {𝑚 𝑠𝑖 𝓎 ∈ {𝓎1 , 𝓎2 , … . . , 𝓎𝑚 }
0 𝑒𝑛 𝑜𝑡𝑟𝑜 𝑐𝑎𝑠𝑜
Donde X y Y son independientes si
𝑓𝑋.𝑌 (𝓍, 𝓎) = 𝑓𝑋 (𝓍) ∙ 𝑓𝑌 (𝓎)
21
Figura 17. Ejemplo de función de probabilidad uniforme bivariada discreta [73].
{0 𝑒𝑛 𝑜𝑡𝑟𝑜 𝑐𝑎𝑠𝑜
1
𝑓𝑋 (𝓍) = { (𝑏 − 𝑎) 𝑠𝑖 𝓍 ∈ {𝑎, 𝑏 }
0 𝑒𝑛 𝑜𝑡𝑟𝑜 𝑐𝑎𝑠𝑜
1
𝑓𝑌 (𝓎) = { (𝑑 − 𝑐) 𝑠𝑖 𝓎 ∈ {𝑐, 𝑑 }
0 𝑒𝑛 𝑜𝑡𝑟𝑜 𝑐𝑎𝑠𝑜
0 𝑠𝑖 𝓍 < 𝑎 𝑜 𝓎 < 𝑐
(𝑥−𝑎)(𝓎−𝑐) 𝑎≤𝓍<𝑏
(𝑏−𝑎)(𝑑−𝑐)
𝑠𝑖 𝑐≤𝓎<𝑑
(𝑥−𝑎) 𝑎≤𝓍<𝑏
𝑠𝑖
𝐹(𝓍, 𝓎) = (𝑏−𝑎) 𝓎≥𝑑
(𝓎−𝑐) 𝑐≤𝓎<𝑑
(𝑑−𝑐)
𝑠𝑖 𝑥≥𝑏
{1 𝑠𝑖 𝓍 ≥ 𝑏 , 𝓎 ≥ 𝑑
22
Figura 20. Ejemplo de función de distribución uniforme bivariada continua [73].
𝜎1 2 𝜌𝜎1 𝜎2
∑=( )
𝜌𝜎1 𝜎2 𝜎2 2
Donde:
1. ∑ es simétrica
−1 1 𝜎2 2 −𝜌𝜎1 𝜎2
3. ∑ = (1−𝜌2)𝜎 ( ).
2 2
1 𝜎2 −𝜌𝜎1 𝜎2 𝜎1 2
𝑎
4. ∑ es positiva definida en (𝑎, 𝑏) ∑ ( ) > 0, con (𝑎, 𝑏) ≠ (0,0)
𝑏
1 1 −1 𝑥 − 𝜇
1
𝑓(𝓍, 𝓎) = 𝑒𝑥𝑝 (− (𝑥 − 𝜇1 , 𝓎 − 𝜇2 ) ∑ (𝓎 − 𝜇 ))
2𝜋√𝑑𝑒𝑡 Σ 2 2
1 1 −1 𝑥
𝑓(𝓍, 𝓎) = 𝑒𝑥𝑝 (− (𝑥, 𝓎) ∑ (𝓎 ))
2𝜋√𝑑𝑒𝑡 Σ 2
Donde:
1 𝜌 −1 1 1 −𝜌
∑=( ) , ∑ = (1−𝜌2) ( )
𝜌 1 −𝜌 1
23
Figura 21. Ejemplo de función de densidad normal bivariada estándar [73].
Ante ello, una mejor opción es la representación gráfica en dos dimensiones para las
funciones 𝑓𝑋 (𝓍) y 𝑓𝑌 (𝓎) cómo se lo indica en la Figura 23.
Para el objetivo de ésta tesis de maestría, basta con representar gráficamente a las
variables X y Y que serán la posición y potencia instalada de GD, tal como se analizará en
la sección 2.4.3 de este trabajo.
24
2. METODOLOGÍA
Figura 24. Esquema de alimentador primario con carga uniformemente distribuida [75].
Donde una carga cualquiera está ubicada a una distancia “x”. Las pérdidas de potencia por
efecto Joule tiene una componente activa y otra reactiva que se lo representa como:
Como una restricción adicional, el generador distribuido a ser insertado solo inyecta
potencia activa por lo que la ecuación 2.1 quedaría:
𝐼1 2 cos ∅ . 𝑅 + 𝐼1 2 sin ∅ . 𝑋 = (𝐼 cos ∅)2 𝑅 − 2(𝐼 cos ∅ . 𝐼𝐺𝐷 )2 𝑅 + 𝐼𝐺𝐷 2 𝑅 + (𝐼 sin ∅)2 𝑋
25
Ya que cos ∅ = 1 y sin ∅ = 0 se simplifica a:
Por lo tanto, la variación de pérdidas por potencia ∆𝑃𝑒𝑟𝑃 se la indica en las ecuaciones (2.3)
y (2.4)
∆𝑃𝑒𝑟𝑃 = 𝐼 2 𝑅 − 𝐼1 2 𝑅 (2.3)
Para calcular las pérdidas totales, la ecuación 2.6 debe ser una función continua e
integrable; por lo tanto:
1
𝑃𝑒𝑟𝑃 = ∫ 𝑑𝑃𝑒𝑟𝑃
𝑥=0
1
𝑃𝑒𝑟𝑃 = (𝐼1 + 𝐼1 𝐼1 + 𝐼2 2 )𝑅
2
(2.7)
Con la inserción de un generador distribuido, como lo indica la Figura 25, se inyecta una
corriente 𝐼𝐺𝐷 y la ecuación 2.5 puede ser redefinida como:
𝐼1 − (𝐼1 − 𝐼2 )𝑥 𝑥1 ≤ 𝑥 ≤ 1
Por lo tanto, las nuevas pérdidas del sistema con un generador distribuido son:
𝑥1 1
′
𝑃𝑒𝑟𝑃 = 3 ∫ [𝐼1 − (𝐼1 − 𝐼2 )𝑥 − 𝐼𝐺𝐷 ]2 𝑅𝑑𝑥 + 3 ∫ [𝐼1 − (𝐼1 − 𝐼2 )𝑥 − 𝐼𝐺𝐷 ]2 𝑅𝑑𝑥
𝑥=0 𝑥=𝑥1
26
Figura 25. Esquema de alimentador primario con carga uniformemente distribuida y un GD [75].
𝐼𝐺𝐷 𝐼2
Siendo: 𝑐= 𝜆=
𝐼1 𝐼1
3𝑥1 𝑐
Δ𝑃𝑒𝑟𝑃 = [2 − 𝑥1 + 𝑥1 𝜆 − 𝑐] (2.11)
1+𝜆+𝜆2
27
Por lo expuesto en la Figura 26, para un valor ideal de 𝜆 = 0 existe un valor definido en 𝑐 y
𝑥1 para una máxima reducción de pérdidas. Sin embargo, al aumentar el valor de 𝜆 a ¼,
para valores de 𝑐 < 0.2 la reducción de pérdidas es menor y la máxima reducción de
pérdidas posible, se lo puede conseguir instalando el GD al final del alimentador primario
𝑥1 ≈ 1, pero matemáticamente la distancia óptima está en 𝑥1 > 1.
28
En el caso de introducir un segundo generador distribuido con una misma potencia al
primero, las ecuaciones 2.9 y 2.13 se pueden modificar de la siguiente manera:
𝑥1 𝑥2
′
𝑃𝑒𝑟𝑃 = 3 ∫ [𝐼1 − (𝐼1 − 𝐼2 )𝑥 − 2𝐼𝐺𝐷 ]2 𝑅𝑑𝑥 + 3 ∫ [𝐼1 − (𝐼1 − 𝐼2 )𝑥 − 𝐼𝐺𝐷 ]2 𝑅𝑑𝑥
𝑥=0 𝑥=𝑥1
1
Si la función objetivo es reducir las pérdidas, al derivar la ecuación 2.16 e igualarla a cero:
1 (2𝑖 − 1)𝑐
𝑥𝑖,𝑜𝑝𝑡 = −
1 − 𝜆 2(1 − 𝜆)
(2.17)
∑(2𝑖 − 1) = 𝑛2
𝑖=1
(2.19)
𝑛
𝑛(𝑛 + 1)
∑𝑖 =
2
𝑖=1
(2.20)
𝑛
𝑛(𝑛 + 1)(2𝑛 + 1)
∑ 𝑖2 =
6
𝑖=1
(2.21)
𝑛
1 𝑛
∑ =
1−𝜆 1−𝜆
𝑖=1
(2.22)
29
Al reemplazar las ecuaciones 2.19, 2.20, 2.21 y 2.22 en 2.18 se obtiene:
3𝛼𝑐 𝑐 2 𝑛(4𝑛2 − 1)
Δ𝑃𝑒𝑟𝑃,𝑜𝑝𝑡 = [𝑛 − 𝑐𝑛2 + ]
1−𝜆 12
(2.23)
Para determinar el factor de compensación 𝑐 adecuado de cada GD, que reduzca las
pérdidas Δ𝑃𝑒𝑟𝑃 a un nivel óptimo, se deriva la ecuación 2.23 e iguala a cero, obteniendo:
2
𝑐=
2𝑛 + 1
(2.24)
En caso de insertar un generador distribuido los resultados preliminares son:
2 2 8𝛼
𝑐= 𝑥1 = Δ𝑃𝑒𝑟𝑃,𝑜𝑝𝑡 =
3 3 (1 − 𝜆) 9 (1 − 𝜆)
De acuerdo a la Figura 24, el alimentador primario es de tipo radial, con una sola troncal y
cargas uniformemente distribuidas. Un alimentador primario real posee: más de una
troncal, circuitos primarios trifásicos y monofásicos, redes desbalanceadas, circuitos
secundarios, cargas distribuidas de forma no uniforme en potencia y ubicación, elementos
de protección y maniobra, capacitores y otros elementos que no fueron considerados en el
análisis anterior. Un ejemplo real es el indicado en la Figura 29.
En un alimentador real existe una topografía compleja por la ubicación y potencia instalada
de las cargas, esto hace que el parámetro 𝜆 no sea factible de calcular mediante funciones
algebraicas ya que, en cada bifurcación de los circuitos, las relaciones de corriente son
diferentes; además, el parámetro 𝜆 depende directamente de la uniformidad de la
distribución de cargas.
30
Figura 29. Diagrama del circuito del alimentador de prueba de 8500 nodos [76].
31
2.2. Uso de OPENDSS para cálculo de pérdidas de potencia.
Sin embargo, no todos los elementos son una relación matricial lineal simple, por lo que la
ecuación de corriente es compleja y está en función de ecuaciones diferenciales como lo
indica la siguiente expresión:
𝐼𝑇𝑒𝑟𝑚 (𝑡) = Ϝ(𝑉𝑇𝑒𝑟𝑚 , [𝐸𝑠𝑡𝑎𝑑𝑜], 𝑡)
La función F varía acorde al tipo de simulación realizada, de tal manera que los elementos
deben ser capaces de ser representados como matrices parciales, al resolver sistemas de
𝜕Ϝ
ecuaciones diferenciales mediante el operador .
𝜕𝑉
Las cargas y los generadores son representados por equivalentes Norton con una 𝑌𝑝𝑟𝑖𝑚 y
una corriente de compensación 𝐼𝑐𝑜𝑚𝑝 o inyección de corriente que compensa la porción no
lineal como lo indica la Figura 31.
32
Figura 31. Representación de equivalente Norton [55].
OPENDSS provee los algoritmos para construir las matrices nodales de admitancia que
formen una nueva matriz de admitancias de todo el sistema 𝑌𝑠𝑖𝑡𝑒𝑚𝑎 y así resolver el primer
sistema de ecuaciones [𝑰] = [𝒀][𝑽]. Luego se realiza el mismo procedimiento con un vector
𝐼𝑐𝑜𝑚𝑝 en cada elemento del sistema por lo que se debe resolver el sistema de ecuaciones
[𝑰𝒄𝒐𝒎𝒑 (𝑽)] = [𝒀𝒔𝒊𝒔𝒕𝒆𝒎𝒂 ] [𝑽]. Este procedimiento se repite interactivamente bajo la siguiente
condición [𝑽𝒏+𝟏 ] = [𝒀𝒔𝒊𝒔𝒕𝒆𝒎𝒂 ]−𝟏 𝐼𝑐𝑜𝑚𝑝 [𝑽𝒏 ] 𝒏 = 𝟎, 𝟏, 𝟐, … .. hasta que el valor de 𝒏 sea tal
que [∆𝑽] = [𝑽𝒏+𝟏 ] − [𝑽𝒏 ] < 𝟎. 𝟎𝟎𝟎𝟎𝟏 𝒑. 𝒖. como lo indica la Figura 32:
33
feeders”, cuyos resultados oficiales fueron verificados utilizando OPENDSS, en especial
con [76] donde el número de nodos hace imposible realizar un cálculo de forma manual.
34
Tabla 2. Principales comandos de funciones en OPENDSS.
Los comentarios son muy útiles al momento de redactar el script que es posible mediante
los caracteres // o ! al inicio de la línea de programación o /*……*/ para múltiples líneas;
mientras que el carácter ∼ o la palabra “more” para seguir editando la línea de
programación en otros reglones.
35
En el caso de las líneas eléctricas los parámetros son la resistencia y la inductancia. La
ecuación de impedancia según [77] es:
(2.24)
En OPENDSS los parámetros que modelen la impedancia de la línea eléctrica se los indica
en la Tabla 3.
36
Tabla 5. Listado de comandos para modelamiento geométrico en líneas de distribución en
OPENDSS.
37
Tabla 7. Listado de comandos para modelamiento de transformadores en OPENDSS.
Ya que la demanda energética está en función del tiempo, es necesario definir la curva de
demanda tipo para cada carga mediante el comando “loadshape” el cual corresponde a un
conjunto de números no enteros entre 0 a 1, que son multiplicados a la potencia instalada
de la carga. Con ello se puede modelar un flujo de potencia para un determinado tiempo,
38
ya sea: diario, mensual, anual, etc. Los parámetros que definen el “loadshape” son los
indicados en la Tabla 9.
39
Los bancos de capacitores poseen diferentes configuraciones de conexión a la red. Si es
en delta, el segundo terminal no existe y si es en estrella el segundo terminal debe ser
especificado si es aterrizado o no en tierra. Los comandos para modelar el “capacitor
object” de los indica en la Tabla 11.
40
2.2.3. Modelación en OPENDSS de la red de prueba de 123 nodos de la
IEEE.
• Modelamiento de la “Vsource”
41
• Modelamiento de los transformadores de distribución. Mediante la instrucción “redirect”
se ejecuta el script “Transformadores_123.dss” que debe estar en la misma carpeta de
Windows donde se ejecuta el script principal.
42
• Los capacitores son modelados en el archivo “Capacitores_123.dss”
• Los seccionadores son elementos de maniobra para la interconexión con otras redes o
elementos de la red y son modelados en el archivo “Seccionadores_123.dss”
43
• Para obtener respuestas gráficas de OPENDSS se conFigura el archivo
“123nodosIEEEcoord.dat” que contiene las coordenadas de los nodos de la siguiente
forma:
• Ya que cada script ejecuta ordenes línea por línea de forma descendente (como
cualquier software de programación en texto), este debe recibir instrucciones para la
resolución matemática del circuito y también las instrucciones para mostrar en forma de
texto y gráficos los resultados obtenidos.
44
2.3. Uso del software Cyme en el modelamiento de redes de
distribución de las empresas eléctricas en Ecuador.
Este software considera diferentes tipos de análisis como lo indicada en la primera parte
de Figura 49 y para propósitos de este trabajo de tesis se utilizará únicamente los flujos de
potencia desequilibrado. Para ello es importante reconocer el entorno de trabajo del Cyme
como lo indica la Figura 36.
45
Figura 37. Creación de perfil de usuario y ubicación de directorios en CYME.
46
(equipos), datos de la topología de cada red e información relevante del proyecto en
estudio. Las etapas son descritas en la Figura 40.
47
Figura 41. Inserción de información en bases de datos Access en CYME.
Figura 42. Alimentadores de la zona noreste de la EEQ cargada en la base de datos en CYME.
48
Figura 43. Alimentadores de la S/E Huachi de la EEASA cargada en la base de datos en CYME.
Lo correcto es activar (cargar) un solo alimentador a la vez del inventario de redes, para
realizar el análisis de flujo de potencia y obtención de resultados por medio de reportes
para extracción de datos. Sin embargo, cuando el modelamiento de los elementos de la
red es incorrecto, CYME lo reporta como inconsistencias de la red, mostrando resultados
erróneos o sin la ejecución del flujo de potencia u otro análisis como lo muestra las Figuras
44 y 45.
49
Figura 45. Reporte errores de modelamiento de un alimentador de la EEASA en CYME.
50
La presentación gráfica de resultados se la puede personalizar de acuerdo a las
preferencias visuales del usuario como lo indica la Figura 47, para este trabajo se ajustaron
los resultados de datos de red como lo indica la Figura 48
Figura 48. Presentación gráfica de datos en la red de prueba de 123 nodos de la IEEE en
CYME.
51
Figura 49. Configuración para ajustes de resultados de flujo de potencia en CYME.
Figura 50. Resultados gráficos de flujo de potencia en la red de prueba de 123 nodos de la IEEE
en CYME.
Entre los elementos más comunes que se identifican en una red de distribución eléctrica
consta: equivalente de fuente (representa la subestación), barra, transformador de
potencia, reguladores de voltaje, interruptores – seccionadores, nodos, líneas aéreas,
líneas subterráneas, transformadores de distribución, cargas concentradas y capacitores.
Las gráficas que representan los elementos eléctricos constan en la Figura 51.
52
software como OPENDSS. La obtención de ésta información es a través de los reportes
como lo indica la Figura 52.
53
Figura 54. Parámetros adicionales de líneas aéreas y subterráneas en CYME.
Cada hoja de reporte CYME debe ser copiada o exportada a Excel, como lo indica la Figura
55. El archivo .xlsx de cada tipo elemento se lo debe guardar como archivo individual .csv
(comma separated values) con la respectiva nomenclatura acorde al tipo de elemento como
se lo indica en la Figura 56.
Figura 56. Archivos .csv con los tipos de elementos eléctricos exportados desde CYME.
54
2.4. Utilización de Matlab para implementación de la
metodología.
b) Para cada línea que representa un elemento eléctrico de la red, se identifica los
caracteres importados y los datos separados por comas
55
e) Los resultados se copian, elimina comillas y pega en un scrip de OPENDSS y se lo
graba como “cargas_Alimentador.dss” como se lo muestra en la Figura 57.
56
Figura 58. Procesamiento de información de capacitores al lenguaje de programación de
OPENDSS.
a. Leer los archivos .csv y asignar a cada variable con una cadena alfanumérica “ { } ”
b. Cada fila representa una configuración geométrica de una línea aérea, así como las
características eléctricas y longitud de cada tramo, las cuales deben ser clasificadas en
variables con una nomenclatura fácilmente identificable para cada propósito de cálculo.
57
c. De acuerdo a la designación de fases monofásica, bifásica o trifásica de cada tramo, se
debe clasificar las variables designadas en el paso b y posterior unión en matrices
mayores, proceso igual al establecido en [77] para obtener la matriz primitiva y de
impedancias.
58
e. Con las matrices debidamente establecidas en los pasos c y d se realizan los cálculos
para la obtención de matrices primitivas de impedancia, primitivas de admitancia,
matrices de coeficiente de potencial eléctrico y finalmente de capacitancia.
59
g. Los resultados se copian, elimina comillas y espacios para pegar en un scrip de
OPENDSS y se lo graba como “lineasaereas_Alimentador.dss”. cómo se lo muestra en
la Figura 59.
a) Leer los archivos .csv y asignar a cada variable con una cadena alfanumérica “ { } ”
b) Identificar los caracteres importados para cada elemento separado por comas y
convertirlos a datos numéricos.
60
c) Se asigna a la variable “asig_subterraneo”, las variables del paso b de acuerdo a las
fases designadas de cada tramo y según la estructura básica de programación de
OPENDSS y se une en una sola cadena de texto alfanumérica para cada elemento.
61
Figura 60. Procesamiento de información de configuración de líneas subterráneas al lenguaje de
programación de OPENDSS.
a. Leer los archivos .csv y asignar a cada variable con una cadena alfanumérica “ { } ”
b. Identificar los caracteres importados para cada elemento separado por comas y
convertirlos a datos numéricos.
62
c. Transformar de caracteres a números únicos para cada nodo a fin de establecer
comparaciones posteriores para la asignación correcta de nodos a cada elemento.
63
d. De acuerdo a la ubicación de los seccionadores – interruptores y transformadores, para
una correcta designación de nomenclatura a los nodos de las líneas aéreas o
subterráneas, es necesario comparar numéricamente las designaciones de nodos de
cada elemento de tal forma que no se repitan y exista una adecuada conexión entre los
elementos y repliquen de forma exacta el modelamiento de la red.
64
e. Identificar todos los nodos a los que se conectan las líneas aéreas o subterráneas a
medio voltaje, para elaborar un listado de nodos (sin que se repita) con los que se puede
conectar uno o varios generadores distribuidos.
65
f. Se asigna a las variables “asig_top_aereo” y “asig_top_subterraneo”, las variables del
paso b y d de acuerdo a las fases designadas de cada tramo entre nodos y según la
estructura básica de programación de OPENDSS y se une en una sola cadena de texto
alfanumérica para cada elemento.
66
h. Se asigna a la variable “asig_InSc” y “asig_Trafo”, las variables del paso b y g de
acuerdo a los seccionadores – interruptores o transformadores conectados a cada tramo
mediante nodos y según la estructura básica de programación de OPENDSS y se une
en una sola cadena de texto alfanumérica para cada elemento.
67
i. Se asigna a la variable “asig_coord_InSc”, “asig_coord_Trafo”, “asig_coord_aereo” y
“asig_coord_subterraneo” las variables de coordenadas de cada nodo de los elementos
y según la estructura básica de programación de OPENDSS que se une en una sola
cadena de texto alfanumérica para cada elemento.
j. Los resultados se copian, elimina comillas, reemplazan espacios y pega en varios scrip
de OPENDSS y se los graba como “Secc_Inter_Alimentador.dss”,
“Transformadores_Alimentador.dss”, “Topologia_Alimentador.dss” y
“Coordenadas_Alimentador.dat” cómo se lo muestra en las Figuras 61, 62, 63 y 64.
68
Figura 62. Procesamiento de información de configuración de transformadores al lenguaje de
programación de OPENDSS.
k. Producir líneas de programación con las variables del paso e para la inserción de un
generador distribuido por cada caso de análisis desde Matlab hacia OPENDSS desde 3
instancias del programa base como lo indica la sección 2.4.4
69
l. Los resultados se copian, elimina comillas, reemplaza espacios, se crea nuevas comillas
y pega en archivos .txt para pegarlos posteriormente en las 3 instancias del scrip base
de Matlab “inicio_Alimentador_un_gen_paralelo.m” que se lo explica en la sección 2.4.4
literal e).
Figura 65. Procesamiento de información para inserción de generación distribuida con todos los
casos posibles desde Matlab a OPENDSS.
70
2.4.2. Interface de comunicación de Matlab con OPENDSS.
Matlab posee controladores llamados ActiveX que son pequeñas aplicaciones integradas
y permiten una comunicación con otro software como el COM Server de OPENDSS de la
Figura 34, la comunicación es bidireccional como la indica la Figura 66. La activación del
COM Server en el entorno Windows se lo realiza en el ambiente DOS a través de la
instrucción Regsvr32 OpenDSSEngine.dll en el directorio de instalación de OPENDSS. El
inicio de las comunicaciones se lo realiza desde Matlab mediante la creación de una función
como se indica a continuación:
71
Al comprender los conceptos planteadas en la sección 1.5.3, Matlab contiene comandos
de programación que replican de forma exacta las funciones de probabilidad deseadas sin
tener que plantear ecuaciones o funciones complejas como las planteadas a manera de
ejemplo en la sección 1.5.3. Se debe elegir la adecuada función de probabilidad que cubra
la mayor parte de soluciones de las variables posición y potencia instalada de la red de
distribución, dentro de un rango finito de valores a ser simulados.
La mejor opción es la función de probabilidad uniforme ya que todos los datos tienen la
misma probabilidad de ser generados únicamente en el rango indicado. Las otras funciones
aptas son: Normal desplazada, inversa Gauss, Burr, Weibull, Logística y Localización a
escala T. Las otras funciones no son aptas para una combinación bivariada ya que
concentran los valores de probabilidad en zonas extremas.
Para una simulación de n casos, con una combinación probabilística bivariada uniforme
discreta, todos los datos generados está dentro del rango finito designado; asi por ejemplo,
para una simulación de 1000 casos, con un rango finito de 0 a 100 para el eje X y de 0 a
10000 para el eje Y cuya gráfica es la indicada en la Figura 68.
72
Figura 68. Simulación de datos para una función probabilística bivariada uniforme discreta.
Los gráficos para otro tipo de combinaciones probabilísticas bivariadas se indica en las
Figuras 69 y 70.
Figura 69. Simulación de datos para diferentes funciones probabilísticas bivariada combinadas.
73
Figura 70. Simulación de datos para diferentes funciones probabilísticas bivariada combinadas.
El criterio para el análisis gráfico de la generación de números aleatorios para cada tipo de
distribución bivariada, primero se observa en las gráficas del extremo izquierdo el rango
finito para el eje X y Y que genera de manera automática Matlab, la concentración o
dispersión de datos para una determinada zona; así como la cantidad de datos atípicos. Lo
importante es observar la cantidad de datos que se encuentran dentro del rango finito
deseado y los que están fuera. Las segundas gráficas muestran las mismas distribuciones
que las del lado izquierdo y a la misma escala de la Figura 68. Aquí se observa la
homogeneidad o no de distribución de datos. La combinación de una función de distribución
con otra, para formar una función bivariada es inmensa; sin embargo, el objetivo de este
trabajo de tesis se requiere uniformidad = homogeneidad, para que todos los casos
posibles tengan la misma probabilidad de ser generados. Tomando en cuenta las
segundas, terceras y cuartas gráficas de las Figuras 69 y 70, la cuantificación de los datos
para 1000 casos aleatorios se los muestra en la Figura 71.
74
Figura 71. Ponderación de casos idóneos en una simulación de 1000 datos para combinaciones
bivariadas de funciones probabilísticas.
Por lo expuesto en la Figura 68, 69, 70 y 71, las mejores combinaciones para funciones de
distribución bivariada son la Uniforme y Weibull.
Este programa realiza las interacciones y comunicación con OPENDSS para el análisis de
resultados obtenidos y determinar la mejor opción de generación distribuida. Para ello es
importante mencionar las diferentes etapas como se indica a continuación:
75
b) Matlab inicia la primera interacción con OPENDSS con la programación indicada en la
sección 2.4.2 para evaluar el estado inicial de la red con la ejecución de un flujo de
potencia con redes desbalanceadas. De la resolución obtenida se hace la exportación
de datos en archivos formato “.csv” de voltajes, corrientes, potencias, cargas,
generadores, pérdidas, capacidad y sobrecargas. También se puede mostrar las
variables antes mencionadas en archivo texto para una visualización preliminar de
resultados, que para este caso lo necesario será voltajes y pérdidas. La presentación
gráfica es importante en el trazado de la red mediante coordenadas de cada nodo para
la presentación de potencia y pérdidas; así como los perfiles de voltaje de cada fase.
c) De los resultados del paso anterior se extraen la potencia activa y aparente máxima, el
factor de potencia y las pérdidas totales. Estos valores son la base para determinar la
potencia máxima de la GD y el porcentaje de pérdidas a ser reducido.
76
d) De la teoría y ejemplos expuestos en la sección 1.5.3 y las combinaciones bivariadas
escogidas en la sección 2.4.3, se genera casos aleatorios para cada función de
probabilidad donde los límites son: el número de nodos en el eje X y la potencia máxima
en el eje Y. Para la función de densidad Weibull es necesario especificar dichos límites
y determinar el número de casos efectivos para someterlos a simulación.
77
De acuerdo al número de casos generados, se representa gráficamente la densidad de
datos como lo indica la Figura 72. Por medio de la instrucción round, se aproxima a
números enteros, en especial para los nodos donde podría ser ubicado un GD.
Figura 72. Vectores aleatorios generados para una distribución probabilística bivariada uniforme
con 500 y 5000 casos en la red de 123 nodos IEEEE.
Los casos generados o también llamados vectores aleatorios pueden ser generados
simultáneamente para 2, 3 o, 4 GD respetando los límites de número de nodos y potencia
máxima admisible a la red.
Figura 73. Vectores aleatorios generados para una distribución probabilística bivariada uniforme,
con 500 casos para 4 GD en la red de 123 nodos IEEEE.
78
e) Para cada vector aleatorio, se debe hacer una simulación de flujo de potencia en
OPENDSS y los resultados de cada caso deben ser exportados para análisis de
reducción de pérdidas. Sin embargo, un gran número de casos representa varias
limitaciones por tiempo de ejecución para cada simulación y recurso computacional
disponible. La manera correcta de reducir tiempos es a través de la computación
paralela como lo indicado en [80], donde Matlab realiza de forma simultánea varias
interacciones con OPENDSS; es decir, se obtienen varios resultados de flujos de
potencia al mismo tiempo. El número de interacciones lo determina la cantidad de
núcleos que posee el procesador del CPU, cuyo fragmento de programación se lo
describe como:
79
En el comentario indicado como “primera programación, asig1_prog” se pega el texto
indicado en la primera parte de la Figura 65. Este texto es generado para la red sujeta al
análisis y contiene cada posible nodo donde se pueda conectar el GD. Los resultados del
caso generado son exportados a un archivo “.csv” y las pérdidas totales de cada caso son
almacenadas en una sola variable. El archivo generado posteriormente es eliminado para
evitar la innecesaria ocupación de espacio en el disco duro.
f) Con el registro de resultados de todos los casos, se selecciona el caso que reduce al
máximo las pérdidas en la red. La selección del caso incluye: el número del nodo
(ubicación), el nombre del nodo, la distancia hacia la subestación y la potencia del GD;
que deberán ser presentados en la ventana de comandos de Matlab. La asignación del
mejor caso debe ser graficada para determinar la localización del vector de GD.
80
g) Con la selección del mejor vector aleatorio, se realiza un nuevo flujo de potencia desde
Matlab hacia OPENDSS con el factor de potencia fp=1 para el GD. Se exportan los
resultados de las mismas variables que se indicaron en el literal b) con la mejor opción.
h) Con los resultados del literal f) para fp=1, se requiere un análisis del GD a diferentes
valores de fp, con el objetivo de seleccionar un mejor valor para la reducción de pérdidas
y verificar los perfiles de voltaje de la red. El rango de valores para el análisis fue
ingresado en el literal a), cuyo fragmento de programación se lo describe como:
81
i) Con el registro de los resultados del literal h) se realiza la selección de la mejor opción
del factor de potencia para la máxima reducción de pérdidas, similar al literal f) se realiza
la exportación de los resultados de las mismas variables de los literales b) y g). Con
estos valores se realiza una última simulación de flujos de potencia que serán los
resultados definitivos.
82
j) Finalmente se leen los valores de los resultados de los voltajes obtenidos en los literales
b), g) e i) para compararlos estadísticamente con el valor promedio y varianza. El
mejoramiento de estos valores se los verifica de forma numérica y gráfica.
83
2.4.5. Modificación del programa base para aproximación – ajuste de
resultados.
El análisis de alimentadores primarios puede ser una tarea compleja y extensa por la gran
cantidad de nodos, lo que haría suponer que, a más nodos, se requieren más vectores
aleatorios como lo indica la Figura 72. Con el objetivo de aumentar la confiabilidad de
resultados, para que la cantidad de vectores aleatorios sean suficientes para el análisis
muestral por probabilidad, se redondean los números que contiene cada vector aleatorio
(discretización) por medio de instrucción round(número,-cifras_aproximadas) y el cambio
se realiza al final del literal d) de la sección 2.4.4 como se indica a continuación:
Al discretizar los valores de ubicación y potencia para cada vector aleatorio las nuevas
graficas se muestra en la Figura 74.
Figura 74. Vectores aleatorios generados con discretización de 2 cifras para una distribución
probabilística bivariada uniforme con 500 casos en la red de 123 nodos IEEE.
Del vector aleatorio preseleccionado, con el mismo procedimiento del literal f) de la sección
anterior y que está representado gráficamente en la Figura 74, se debe hacer una nueva
84
interacción en la vecindad del vector aleatorio preseleccionado, siguiendo nuevamente el
procedimiento desde los literales d, e, f y g de la sección anterior, la nueva representación
gráfica se lo indica en la Figura 75.
Figura 75. Vectores aleatorios generados en la vecindad del vector preseleccionado para una
distribución probabilística bivariada uniforme con 500 casos en la red de 123 nodos
IEEE.
Figura 76. Vectores aleatorios generados en la vecindad del vector preseleccionado para una
distribución probabilística bivariada uniforme con 5000 casos en la red de 123 nodos
IEEEE.
85
las realizan en la cabecera de cada alimentador. También es necesario estimar el recurso
energético disponible para generación Eólica y Solar. El software HOMER (Hybrid
Optimization of Multiple Energy Resources) [81] permite la descarga de esta información a
través de los datos a nivel mundial que posee la NASA como de lo indica en la Figura 77.
Figura 77. Información de recurso energético eólico y solar para Ecuador [81].
b) Ingresar el número de nodo asignado como resultado obtenido en el ítem 2.4.5 para el
cálculo de la nueva potencia del GD.
c) Para el GD a ser insertado, se debe asignar los datos del recurso energético eólico o
solar obtenido. También se asigna la curva de perfil de carga de acuerdo al registro en
la cabecera del alimentador. Los datos de cada tipo son almacenados en un archivo
“.csv” como lo indica la Figura 78.
86
Figura 78. Archivos para asignación de perfiles de generación y carga.
d) Se eliminan otros archivos generados y la carpeta que contiene cada caso para evitar
ocupar espacio innecesario en el disco duro.
El resto del programa tiene el mismo procedimiento que el indicado en la sección 2.4.4.
87
3. RESULTADOS Y DISCUSIÓN
3.1. Resultados
Debido a la buena confiabilidad de los datos de las redes de prueba de la IEEE y que se
encuentran correctamente modelados en CYME, las pruebas se las realiza sobre las redes
de 13, 34, 37 y 123 nodos.
Tabla 13. Resultados importados a Matlab para red original de 13 nodos IEEE sin GD.
Variable Valor
Potencia activa del sistema: 3577.7000 [kW]
Potencia aparente del sistema: 3970.7611 [kVA]
Factor de potencia: 0.9010
Pérdidas por potencia activa: 110.5046 [kW]
El porcentaje de pérdidas: 3,0887%
Figura 79. Resultados directos de OPENDSS para la red original de 13 nodos IEEE sin GD.
Figura 80. Perfiles de voltaje para la red original de 13 nodos IEEE sin GD.
88
Al aplicar el programa base de la sección 2.4.4, a la red original de 13 nodos IEEE, se
obtiene lo indicado en la Tabla 14, la Figura 81 indica los resultados numéricos obtenidos
de la simulación en OPENDSS para factor de potencia fp=1 y fp=0.94 capacitivo. En la
Figura 82 se presentan los resultados gráficos del lugar y la potencia del generador
distribuido. Finalmente, en la Figura 83 indica los perfiles de voltaje para la red original de
13 nodos IEEE con GD con los factores de potencia 1 y 0.94 capacitivo.
Tabla 14. Resultados importados a Matlab para red original de 13 nodos IEEE con un GD.
Variable Valor
El mejor lugar del generador 1: nodo 671
A una distancia desde la S/E de: 4000.00 [ft]
A una distancia desde la S/E de: 1219.20 [m]
Con una potencia de: 3100.00 [kW]
Pérdida mínimizada a: 33.4220 [kW]
Porcentaje de disminución de pérdida: 69.76 %
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.9400
Pérdida mínimizada a: 23.8076 [kW]
Porcentaje de disminución de pérdida: 78.46 %
Con una ganacia adicional de: 8.70 %
Con respecto al generador solo con potencia activa
Figura 81. Resultados directos de OPENDSS para la red original de 13 nodos IEEE con un GD.
1 GD
Figura 82. Lugar y potencia de GD para red original de 13 nodos IEEE con un GD.
89
Factor de Potencia = 1 Factor de Potencia = 0.94
Figura 83. Perfiles de voltaje para la red original de 13 nodos IEEE con GD.
El análisis estadístico de las Figuras 80 y 83, para perfiles de voltaje de las fases ABC se
indica en la Figura 84 mediante diagramas de “caja de bigotes” y los valores promedio más
la desviación estándar que muestra la mejora numérica de los valores de voltaje:
Figura 84. Análisis estadístico de los perfiles de voltaje para la red original de 13 nodos IEEE.
90
nodo 650 y 632); por lo que se suprimen estos elementos para realizar una nueva
simulación y analizar el comportamiento de la red como se indica a continuación:
Tabla 15. Resultados importados a Matlab para red modificada de 13 nodos IEEE sin GD.
Variable Valor
Potencia activa del sistema: 3328.4000 [kW]
Potencia aparente del sistema: 4054.2010 [kVA]
Factor de potencia: 0.8210
Pérdidas por potencia activa: 126.3966 [kW]
El porcentaje de pérdidas: 3,7975%
Figura 85. Resultados directos de OPENDSS para la red modificada de 13 nodos IEEE sin GD.
Figura 86. Perfiles de voltaje para la red modificada de 13 nodos IEEE sin GD.
Tabla 16. Resultados importados a Matlab para red modificada de 13 nodos IEEE con un GD.
Variable Valor
El mejor lugar del generador 1: nodo 671
A una distancia desde la S/E de: 4000.00 [ft]
A una distancia desde la S/E de: 1219.20 [m]
Con una potencia de: 2900.00 [kW]
Pérdida mínimizada a: 54.9675 [kW]
Porcentaje de disminución de pérdida: 56.51 %
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.8400
Pérdida mínimizada a: 27.6475 [kW]
Porcentaje de disminución de pérdida: 78.13 %
Con una ganacia adicional de: 21.61 %
Con respecto al generador solo con potencia activa
91
Figura 87. Lugar y potencia de GD para red modificada de 13 nodos IEEE con un GD.
El análisis estadístico de las Figuras 86 y 88, para las fases ABC se indica en la Figura 89:
92
Figura 89. Análisis estadístico de los perfiles de voltaje para la red modificada de 13 nodos
IEEE.
Tabla 17. Resultados importados a Matlab para red modificada de 13 nodos IEEE con un GD y
mejoramiento con capacitor - inductor.
Variable Valor
El mejor lugar del Capacitor 1 Fase A: nodo 634
A una distancia desde la S/E de: 2510.00 [ft]
A una distancia desde la S/E de: 765.05 [m]
Con una reactancia de: -25.00 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 692
A una distancia desde la S/E de: 4010.00 [ft]
A una distancia desde la S/E de: 1222.25 [m]
Con una reactancia de: 144.00 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 671
A una distancia desde la S/E de: 4000.00 [ft]
A una distancia desde la S/E de: 1219.20 [m]
Con una reactancia de: -324.00 [kVAr]
Pérdidas por potencia activa: 24.5926 [kW]
Porcentaje de disminución de pérdida: 80.5433 %
Con una ganacia adicional de: 2.41 %
El porcentaje de pérdidas de potencia activa es: 0.7389 %
93
Figura 90. Lugar y reactancia de Capacitores para red modificada de 13 nodos IEEE con un GD
y FP.
Figura 91. Perfiles de voltaje para la red modificada de 13 nodos IEEE con GD y mejorada con
capacitores – inductores.
El análisis estadístico de las Figuras 86, 88 y 91, para las fases ABC se indica en la Figura
92:
94
Figura 92. Análisis estadístico de los perfiles de voltaje para la red modificada de 13 nodos IEEE
con mejora de capacitores.
Para el análisis de la red de prueba de 34 nodos IEEE se verifica que posee 2 reguladores
de voltaje entre los nodos 814 – 850 y 832 – 852 y capacitores en los nodos 844 y 848
según la Figura 93. El comportamiento de la red original 34 nodos IEEE se lo muestra en
la Figura 94.
Figura 94. Comportamiento de la red original de prueba 34 nodos IEEE sin GD.
Como se puede verificar en la Figura 93, hay incidencia directa de los reguladores de
voltaje y en menor medida de los capacitores; por lo tanto, se utilizará el mismo criterio que
en la red de 13 nodos IEEE. Los resultados de la red modificada 34 nodos IEEE se los
presenta en las tablas 18, 19, 20 y las Figuras 95, 96, 97, 98 y 99.
95
Tabla 18. Resultados importados a Matlab para red modificada 34 nodos IEEE sin GD.
Variable Valor
Potencia activa del sistema: 1542.6000 [kW]
Potencia aparente del sistema: 1737.9341 [kVA]
Factor de potencia: 0.8876
Pérdidas por potencia activa: 215.4289 [kW]
El porcentaje de pérdidas: 13,9653%
Figura 95. Perfiles de voltaje para la red modificada de 34 nodos IEEE sin GD.
Tabla 19. Resultados importados a Matlab para red modificada de 34 nodos IEEE con un GD.
Variable Valor
El mejor lugar del generador 1: nodo 834
A una distancia desde la S/E de: 183670.00 [ft]
A una distancia desde la S/E de: 55982.62 [m]
Con una potencia de: 1300.0000 [kW]
Pérdida mínimizada a: 98.0669 [kW]
Porcentaje de disminución de pérdida: 54.48 %
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.8500
Pérdida mínimizada a: 48.2561 [kW]
Porcentaje de disminución de pérdida: 77.60 %
Con una ganacia adicional de: 23.12 %
Con respecto al generador solo con potencia activa
Figura 96. Lugar y potencia de GD para red modificada de 34 nodos IEEE con un GD.
96
a) Factor de Potencia = 1 b) Factor de Potencia = -0.85
Figura 97. Perfiles de voltaje para la red modificada de 34 nodos IEEE con GD.
Tabla 20. Resultados importados a Matlab para red modificada de 34 nodos IEEE con un GD y
mejoramiento con capacitor - inductor.
Variable Valor
El mejor lugar del Capacitor 1 Fase A: nodo 820
A una distancia desde la S/E de: 153950.00 [ft]
A una distancia desde la S/E de: 46923.96 [m]
Con una reactancia de: -162.00 [kVAr]
El mejor lugar del Capacitor 1 Fase B: nodo 890
A una distancia desde la S/E de: 183500.00 [ft]
A una distancia desde la S/E de: 55930.80 [m]
Con una reactancia de: -60.00 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 890
A una distancia desde la S/E de: 183500.00 [ft]
A una distancia desde la S/E de: 55930.80 [m]
Con una reactancia de: -41.00 [kVAr]
Pérdidas por potencia activa: 46.8187 [kW]
Porcentaje de disminución de pérdida: 78.2672 %
Con una ganacia adicional de: 0.67 %
El porcentaje de pérdidas de potencia activa es: 3,0351%
Figura 98. Lugar y reactancia de capacitores para red modificada de 34 nodos IEEE con un GD
y FP.
97
Figura 99. Perfiles de voltaje para la red modificada de 34 nodos IEEE con GD y mejorada con
capacitores.
El análisis estadístico de las Figuras 95, 97 y 99 para las fases ABC se indica en la Figura
100:
Figura 100. Análisis estadístico de los perfiles de voltaje para la red modificada de 13 nodos IEEE
con mejora de capacitores.
Para el análisis de la red de prueba de 123 nodos IEEE se verifica que posee 3 reguladores
de voltaje entre los nodos 150 – 149, 160 – 67 y 25 –26 y capacitores en los nodos 83, 88,
98
90 y 92 según la Figura 101. El comportamiento de la red original 123 nodos IEEE se lo
muestra en la Figura 102.
Figura 102. Comportamiento de la red original de prueba 123 nodos IEEE sin GD.
Como se puede verificar en la Figura 102, hay incidencia directa de los reguladores de
voltaje y en menor medida de los capacitores; por lo tanto, se utilizará el mismo criterio que
en la red de 13 y 34 nodos IEEE. Los resultados de la red modificada 123 nodos IEEE se
los presenta en las tablas 21, 22, 23 y las Figuras 103, 104, 105, 106 y 107.
99
Tabla 21. Resultados importados a Matlab para red modificada 123 nodos IEEE sin GD.
Variable Valor
Potencia activa del sistema: 3442.2000 [kW]
Potencia aparente del sistema: 4005.1214 [kVA]
Factor de potencia: 0.8594
Pérdidas por potencia activa: 109.1186 [kW]
El porcentaje de pérdidas: 3,1700%
Figura 103. Perfiles de voltaje para la red modificada de 123 nodos IEEE sin GD.
Tabla 22. Resultados importados a Matlab para red modificada de 123 nodos IEEE con un GD.
Variable Valor
El mejor lugar del generador 1: nodo 67
A una distancia desde la S/E de: 3377.00 [ft]
A una distancia desde la S/E de: 1029.31 [m]
Con una potencia de: 2200.00 [kW]
Pérdida mínimizada a: 47.3908 [kW]
Porcentaje de disminución de pérdida: 56.57 %
Porcentaje pérdidas potencia activa: 1,3768%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.8600
Pérdida mínimizada a: 27.6191 [kW]
Porcentaje de disminución de pérdida: 74.69 %
Con una ganacia adicional de: 18.12 %
Con respecto al generador solo con potencia activa
Porcentaje pérdidas potencia activa: 0,8024%
Porcentaje de aporte de potencia de GD 63,9126%
100
Figura 104. Lugar y potencia de GD para red modificada de 123 nodos IEEE con un GD.
Tabla 23. Resultados importados a Matlab para red modificada de 123 nodos IEEE con un GD y
mejoramiento con capacitor - inductor.
Variable Valor
El mejor lugar del Capacitor 1 Fase A: nodo 18
A una distancia desde la S/E de: 2025.00 [ft]
A una distancia desde la S/E de: 617.22 [m]
Con una reactancia de: -401.00 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 97
A una distancia desde la S/E de: 3627.00 [ft]
A una distancia desde la S/E de: 1105.51 [m]
Con una reactancia de: 72.00 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 34
A una distancia desde la S/E de: 1350.00 [ft]
A una distancia desde la S/E de: 411.48 [m]
Con una reactancia de: -293.00 [kVAr]
Pérdidas por potencia activa: 26.4380 [kW]
Porcentaje de disminución de pérdida: 75.7713 %
Con una ganacia adicional de: 1.08 %
El porcentaje de pérdidas de potencia activa es: 0.7681 %
101
Figura 106. Lugar y reactancia de capacitores - inductor para red modificada de 123 nodos IEEE
con un GD y FP.
Figura 107. Perfiles de voltaje para la red modificada de 123 nodos IEEE con un GD y mejorada
con capacitores.
El análisis estadístico de las Figuras 103, 105 y 107 para las fases ABC se indica en la
Figura 108:
102
Figura 108. Análisis estadístico de los perfiles de voltaje para la red modificada de 123 nodos
IEEE con un GD y mejora con capacitores.
Aplicando el programa la sección 2.4.5 con 1 GD en la red modificada de 123 nodos IEEE
se obtienen los resultados con mayor precisión como se indica en las tablas 24, 25 y
Figuras 109, 110 y 111.
Tabla 24. Resultados importados a Matlab para red modificada de 123 nodos IEEE con un GD
por aproximación – ajuste.
Variable Valor
El mejor lugar del generador 1: nodo 67
A una distancia desde la S/E de: 3377.00 [ft]
A una distancia desde la S/E de: 1029.31 [m]
Con una potencia de: 2290.00 [kW]
Pérdida mínimizada a: 47.3069 [kW]
Porcentaje de disminución de pérdida: 56.65 %
Porcentaje pérdidas potencia activa: 1,3743%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.8700
Pérdida mínimizada a: 27.4017 [kW]
Porcentaje de disminución de pérdida: 74.89 %
Con una ganacia adicional de: 18.24 %
Con respecto al generador solo con potencia activa
Porcentaje pérdidas potencia activa: 0,7961%
Porcentaje de aporte de potencia de GD 66,5272%
a) Aproximación b) Ajuste
Figura 109. Lugar y potencia de GD para red modificada de 123 nodos IEEE con un por
aproximación – ajuste.
103
a) Factor de Potencia = 1 b) Factor de Potencia = -0.87
Figura 110. Perfiles de voltaje para la red modificada de 123 nodos IEEE con un GD por
aproximación – ajuste.
Tabla 25. Resultados importados a Matlab para red modificada de 123 nodos IEEE con un GD y
mejoramiento con capacitor - inductor por aproximación – ajuste.
Variable Valor
El mejor lugar del Capacitor 1 Fase A: nodo 35
A una distancia desde la S/E de: 2401.00 [ft]
A una distancia desde la S/E de: 731.82 [m]
Con una reactancia de: -381.00 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 101
A una distancia desde la S/E de: 3878.00 [ft]
A una distancia desde la S/E de: 1182.01 [m]
Con una reactancia de: 78.00 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 35
A una distancia desde la S/E de: 2401.00 [ft]
A una distancia desde la S/E de: 731.82 [m]
Con una reactancia de: -224.00 [kVAr]
Pérdidas por potencia activa: 25.2231 [kW]
Porcentaje de disminución de pérdida: 76.8847 %
Con una ganacia adicional de: 1.99 %
El porcentaje de pérdidas de potencia activa es: 0.7328 %
a) Aproximación b) Ajuste
Figura 111. Lugar y reactancia de capacitores e inductor para red modificada de 123 nodos IEEE
con un GD y FP por aproximación – ajuste.
104
a) Aproximación b) Ajuste
Figura 112. Perfiles de voltaje para la red modificada de 123 nodos IEEE con un GD y mejorada
con capacitores e inductor por aproximación – ajuste.
El análisis estadístico de las Figuras 110 y 112 para las fases ABC se indica en la Figura
113:
Figura 113. Análisis estadístico de los perfiles de voltaje para la red modificada de 123 nodos
IEEE con un GD y mejora con capacitores e inductor por aproximación – ajuste.
105
Al utilizar el programa base de las secciones 2.4.4 (directo), 2.4.5 (aproximación – ajuste)
y el Anexo D para 2 GD en la red modificada de 123 nodos IEEE, se obtienen los resultados
indicados en las tablas 26, 27 y Figuras 114, 115, 116, 117, 118 y 119.
Tabla 26. Resultados importados a Matlab para red modificada de 123 nodos IEEE con dos GD
de forma directa y por aproximación – ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del generador 1: nodo 60 nodo 67
A una distancia desde la S/E de: 3026.00 [ft] 3377.00 [ft]
A una distancia desde la S/E de: 922.32 [m] 1029.31 [m]
Con una potencia de: 2400.0000 [kW] 2290.0000 [kW]
Pérdida mínimizada a: 48.7033 [kW] 47.3069 [kW]
Porcentaje de disminución de pérdida: 55.37 % 56.65 %
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.8700 -0.8700
Pérdida mínimizada a: 20.6134 [kW] 18.1945 [kW]
Porcentaje de disminución de pérdida: 81.11 % 83.33 %
Con una ganacia adicional de: 18.13 % 18.34 %
Con respecto al generador 1 solo con potencia activa
El mejor lugar del generador 2: nodo 47 nodo 47
A una distancia desde la S/E de: 3351.00 [ft] 3351.00 [ft]
A una distancia desde la S/E de: 1021.38 [m] 1021.38 [m]
Con una potencia de: 800.0000 [kW] 880.0000 [kW]
Pérdida mínimizada a: 40.3971 [kW] 38.2119 [kW]
Porcentaje de disminución de pérdida: 62.98 % 64.98 %
Porcentaje pérdidas potencia activa: 1,1736% 1,1101%
Con el ajuste del valor del factor de potencia al G2...
Valor calculado es: -0.8400 -0.8600
Pérdida mínimizada a: 17.6674 [kW] 15.0090 [kW]
Porcentaje de disminución de pérdida: 83.81 % 86.25 %
Con una ganacia adicional de: 20.83 % 21.26 %
Con respecto al generador 2 solo con potencia activa
Porcentaje pérdidas potencia activa: 0,5133% 0,4360%
Porcentaje de aporte de potencia de GD 92,9638% 92,0923%
106
a) Primer generador b) Segundo Generador
Figura 115. Acercamiento al lugar y potencia de GD para red modificada de 123 nodos IEEE con
dos GD por aproximación – ajuste.
107
Tabla 27. Resultados importados a Matlab para red modificada de 123 nodos IEEE con dos GD y
mejoramiento con capacitor - inductor mediante aproximación y ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del Capacitor 1 Fase A: nodo 18 nodo 35
A una distancia desde la S/E de: 2025.00 [ft] 2401.00 [ft]
A una distancia desde la S/E de: 617.22 [m] 731.82 [m]
Con una reactancia de: -181.00 [kVAr] -166.00 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 62 nodo 60
A una distancia desde la S/E de: 3276.00 [ft] 3026.00 [ft]
A una distancia desde la S/E de: 998.52 [m] 922.32 [m]
Con una reactancia de: 196.00 [kVAr] 215.00 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 44 nodo 34
A una distancia desde la S/E de: 3101.00 [ft] 1350.00 [ft]
A una distancia desde la S/E de: 945.18 [m] 411.48 [m]
Con una reactancia de: -37.00 [kVAr] -20.48 [kVAr]
Pérdidas por potencia activa: 17.8031 [kW] 14.6931 [kW]
Porcentaje de disminución de pérdida: 83.6846 % 86.5347 %
Con una ganacia adicional de: -0.12 % -0.29 %
El porcentaje de pérdidas de potencia activa es: 0.5172 % 0.4269 %
108
El análisis estadístico de las Figuras 115, 117 y 119 para las fases ABC se indica en la
Figura 120:
Figura 120. Análisis estadístico de los perfiles de voltaje para la red modificada de 123 nodos
IEEE con dos GD y mejora con capacitores e inductor por aproximación – ajuste.
Tabla 28. Resultados importados a Matlab para red modificada de 123 nodos IEEE con tres GD
de forma directa y por aproximación – ajuste.
110
a) Primer generador b) Segundo Generador c) Tercer Generador
Figura 122. Acercamiento al lugar y potencia de GD para red modificada de 123 nodos IEEE con
tres GD por aproximación – ajuste.
111
Tabla 29. Resultados importados a Matlab para red modificada de 123 nodos IEEE con tres GD y
mejoramiento con capacitor - inductor mediante aproximación y ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del Capacitor 1 Fase A: nodo 110 nodo 113
A una distancia desde la S/E de: 5228.00 [ft] 5878.00 [ft]
A una distancia desde la S/E de: 1593.49 [m] 1791.61 [m]
Con una reactancia de: -60.0000 [kVAr] -26.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 101 nodo 57
A una distancia desde la S/E de: 3878.00 [ft] 2276.00 [ft]
A una distancia desde la S/E de: 1182.01 [m] 693.72 [m]
Con una reactancia de: 124.0000 [kVAr] 185.0000 [kVAr]
El mejor lugar del Inductor 1 Fase C: nodo 8 nodo 34
A una distancia desde la S/E de: 900.00 [ft] 1350.00 [ft]
A una distancia desde la S/E de: 274.32 [m] 411.48 [m]
Con una reactancia de: -38.0000 [kVAr] 20.4750 [kVAr]
Pérdidas por potencia activa: 12.3824 [kW] 12.8507 [kW]
Porcentaje de disminución de pérdida: 88.6523 % 88.2232 %
Con una ganacia adicional de: 0.20 % 0.0153 %
El porcentaje de pérdidas de potencia activa es: 0.3597 % 0.3733 %
112
El análisis estadístico de las Figuras 123, 124 y 126 para las fases ABC se indica en la
Figura 127:
Figura 127. Análisis estadístico de los perfiles de voltaje para la red modificada de 123 nodos
IEEE con dos GD y mejora con capacitores e inductor por aproximación – ajuste.
Finalmente se analiza las pérdidas de energía para la red de 123 nodos IEEE como lo
indica la sección 2.4.6. Se utiliza la ubicación de los generadores y capacitores – inductores
obtenidos en el resultado de las tablas 24, 25, 26, 27, 28 y 29 para 1GD, 2GD y 3GD
respectivamente. Los resultados se los indica en desde la Tabla 30 hasta la 36 y desde la
Figura 128 hasta la 142.
Tabla 30. Resultados importados a Matlab para red modificada 123 nodos IEEE con perfil de
energía sin GD.
Variable Valor
Energia activa consumida: 7289.7630 [MWh]
Energía aparente consumida: 8368.3804 [MVAh]
Pérdidas por potencia activa: 75.7680 [MWh]
El porcentaje de pérdidas: 1.0394%
113
Figura 128. Perfiles de voltaje para la red modificada de 123 nodos IEEE con perfil de energía y
sin GD.
Tabla 31. Resultados importados a Matlab para red modificada de 123 nodos IEEE con perfil de
energía y un GD.
Variable Valor directo
El mejor lugar del generador 1: nodo 67
A una distancia desde la S/E de: 3377.00 [ft]
A una distancia desde la S/E de: 1029.31 [m]
Con una potencia de: 1600.00 [kW]
Pérdida mínimizada a: 54.1643 [kWh]
Porcentaje de disminución de pérdida: 28.51 %
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.8800
Pérdida mínimizada a: 47.4953 [MWh]
Porcentaje de disminución de pérdida: 37.31 %
Con una ganacia adicional de: 8.80 %
Con respecto al generador solo con potencia activa
Total de pérdidas de energía por potencia: 0.6515 %
Figura 129. Ajuste de potencia de GD para red modificada de 123 nodos IEEE con energía y un
GD.
114
a) Factor de Potencia = 1 b) Factor de Potencia = -0.88
Figura 130. Perfiles de voltaje para la red modificada de 123 nodos IEEE con perfil de energía y
un GD.
Tabla 32. Resultados importados a Matlab para red modificada de 123 nodos IEEE con perfil de
energía, un GD y mejoramiento con capacitores.
Variable Valor directo
El mejor lugar del Capacitor 1 Fase A: nodo 18
A una distancia desde la S/E de: 2025.00 [ft]
A una distancia desde la S/E de: 617.22 [m]
Con una reactancia de: -114.00 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 97
A una distancia desde la S/E de: 3627.00 [ft]
A una distancia desde la S/E de: 1105.51 [m]
Con una reactancia de: -119.00 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 34
A una distancia desde la S/E de: 1350.00 [ft]
A una distancia desde la S/E de: 411.48 [m]
Con una reactancia de: -97.00 [kVAr]
Pérdidas por energía: 51.1312 [MWh]
Porcentaje de disminución de pérdida: 32.5161 %
Con una ganacia adicional solo con GD de: -4.79 %
Total de pérdidas de energía por potencia activa: 0.7014 %
Figura 131. Lugar y reactancia de capacitores para red modificada de 123 nodos IEEE con perfil
de energía, un GD y FP.
115
Figura 132. Perfiles de voltaje para la red modificada de 123 nodos IEEE con perfil de energía, un
GD y mejorada con capacitores.
El análisis estadístico de las Figuras 128, 130 y 132 para las fases ABC se indica en la
Figura 133:
Figura 133. Análisis estadístico de los perfiles de voltaje para la red modificada de 123 nodos
IEEE con perfiles de energía, un GD y mejora con capacitores.
116
Análisis de energía con 2GD
Tabla 33. Resultados importados a Matlab para red modificada de 123 nodos IEEE con perfil de
energía y dos GD.
Variable GD1 GD2
El mejor lugar del generador: nodo 67 nodo 47
A una distancia desde la S/E de: 3377.00 [ft] 3351.00 [ft]
A una distancia desde la S/E de: 1029.31 [m] 1021.38 [m]
Con una potencia de: 1600.00 [kW] 500.00 [kW]
Pérdida mínimizada a: 54.1643 [MWh] 58.0371 [MWh]
Porcentaje de disminución de pérdida: 28,5130% 23,4016%
Con el ajuste del valor del factor de potencia al generador...
Valor calculado es: -0.8800 -0.8500
Pérdida mínimizada a: 47.4953 [MWh] 41.9781 [MWh]
Porcentaje de disminución de pérdida: 37,3148% 44,5965%
Con una ganacia adicional de: 8,8019% 21,1950%
Con respecto al generador solo con potencia activa
Total de pérdidas de energía por potencia activa: 0,6515% 0,5758%
Figura 134. Ajuste de potencia de GD para red modificada de 123 nodos IEEE con energía y dos
GD.
117
Tabla 34. Resultados importados a Matlab para red modificada de 123 nodos IEEE con perfil de
energía, dos GD y mejoramiento con capacitores.
Variable Valor directo
El mejor lugar del Capacitor 1 Fase A: nodo 35
A una distancia desde la S/E de: 2401.00 [ft]
A una distancia desde la S/E de: 731.82 [m]
Con una reactancia de: -280.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 60
A una distancia desde la S/E de: 3026.00 [ft]
A una distancia desde la S/E de: 922.32 [m]
Con una reactancia de: -474.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 34
A una distancia desde la S/E de: 1350.00 [ft]
A una distancia desde la S/E de: 411.48 [m]
Con una reactancia de: -270.0000 [kVAr]
Pérdidas por energía: 152.5056 [MWh]
Porcentaje de disminución de pérdida: -101.2796 %
Con una ganacia adicional solo con GD de: -138.5896 %
Total de pérdidas de energía por potencia activa: 3.7113%
Figura 136. Lugar y reactancia de capacitores para red modificada de 123 nodos IEEE con perfil
de energía, dos GD y FP.
Figura 137. Perfiles de voltaje para la red modificada de 123 nodos IEEE con perfil de energía,
dos GD y mejorada con capacitores.
118
El análisis estadístico de las Figuras 128, 135 y 137 para las fases ABC se indica en la
Figura 138:
Figura 138. Análisis estadístico de los perfiles de voltaje para la red modificada de 123 nodos
IEEE con perfiles de energía, dos GD y mejora con capacitores.
Tabla 35. Resultados importados a Matlab para red modificada de 123 nodos IEEE con perfil de
energía y tres GD.
119
Figura 139. Ajuste de potencia de GD para red modificada de 123 nodos IEEE con energía y tres
GD.
Tabla 36. Resultados importados a Matlab para red modificada de 123 nodos IEEE con perfil de
energía, tres GD y mejoramiento con capacitores.
120
A una distancia desde la S/E de: 693.72 [m]
Con una reactancia de: -682.00 [kVAr]
El mejor lugar del Inductor 1 Fase C: nodo 34
A una distancia desde la S/E de: 1350.00 [ft]
A una distancia desde la S/E de: 411.48 [m]
Con una reactancia de: -568.00 [kVAr]
Pérdidas por potencia activa: 240.7809 [MWh]
Porcentaje de disminución de pérdida: 217.7870 %
Con una ganacia adicional de: 255.0970 %
Total de pérdidas de energía por potencia activa: 5,8595%
Figura 141. Lugar y reactancia de capacitores para red modificada de 123 nodos IEEE con perfil
de energía, tres GD y FP.
Figura 142. Perfiles de voltaje para la red modificada de 123 nodos IEEE con perfil de energía,
tres GD y mejorada con capacitores.
El análisis estadístico de las Figuras 128, 139 y 141 para las fases ABC se indica en la
Figura 143:
121
Figura 143. Análisis estadístico de los perfiles de voltaje para la red modificada de 123 nodos
IEEE con perfiles de energía, tres GD y mejora con capacitores.
Tabla 37. Resultados importados a Matlab para alimentador 49A EEQ sin GD.
Variable Valor
Potencia activa del sistema: 5915.4000 [kW]
Potencia aparente del sistema: 6078.1219 [kVA]
Factor de potencia: 0.9732
Pérdidas por potencia activa: 236.6355 [kW]
El porcentaje de pérdidas: 4,0003%
Figura 144. Perfiles de voltaje para alimentador 49A EEQ sin GD.
122
Tabla 38. Resultados importados a Matlab para alimentador 49A EEQ con un GD de forma
directa y por aproximación – ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del generador 1: nodo 1034261MA nodo 1031391MA
A una distancia desde la S/E de: 138278.60 [ft] 139248.70 [ft]
A una distancia desde la S/E de: 42147.32 [m] 42443.00 [m]
Con una potencia de: 1900.0000 [kW] 1890.0000 [kW]
Pérdida mínimizada a: 108.1166 [kW] 107.8529 [kW]
Porcentaje de disminución de pérdida: 54,3109% 54,4224%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.9800 -0.9800
Pérdida mínimizada a: 102.7964 [kW] 102.5600 [kW]
Porcentaje de disminución de pérdida: 56,5592% 56,6591%
Con una ganacia adicional de: 2,2483% 2,2367%
Con respecto al generador solo con potencia activa
El porcentaje de pérdidas de potencia activa es: 1,7378% 1,7338%
Porcentaje de aporte de potencia de GD 32,1196% 31,9505%
123
Tabla 39. Resultados importados a Matlab para alimentador 49A EEQ con un GD y mejoramiento
con capacitor – inductor de forma directa y por aproximación – ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del Capacitor 1 Fase A: nodo 45478MA nodo 681301MA
A una distancia desde la S/E de: 82539.10 [ft] 85957.60 [ft]
A una distancia desde la S/E de: 25157.92 [m] 26199.88 [m]
Con una reactancia de: -166.0000 [kVAr] -161.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 82602MA nodo 70628MA
A una distancia desde la S/E de: 88709.60 [ft] 58050.30 [ft]
A una distancia desde la S/E de: 27038.69 [m] 17693.73 [m]
Con una reactancia de: -120.0000 [kVAr] -167.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 165478MA nodo 114357MA
A una distancia desde la S/E de: 808.80 [ft] 35400.60 [ft]
A una distancia desde la S/E de: 246.52 [m] 10790.10 [m]
Con una reactancia de: 223.0000 [kVAr] 230.0000 [kVAr]
Pérdidas por potencia activa: 107.1703 [kW] 112.7295 [kW]
Porcentaje de disminución de pérdida: 54,7108% 52,3615%
Con una ganacia adicional de: -1,8484% -4,2975%
El porcentaje de pérdidas de potencia activa es: 1,8117% 1,9057%
124
El análisis estadístico de las Figuras 144, 146 y 148 para las fases ABC se indica en la
Figura 149:
Figura 149. Análisis estadístico de los perfiles de voltaje para alimentador 49A EEQ con un GD y
mejora con capacitores.
125
Tabla 40. Resultados importados a Matlab para alimentador 49A EEQ con dos GD de forma
directa y por aproximación – ajuste.
126
a) Primer generador b) Segundo Generador
Figura 152. Acercamiento al lugar y potencia de GD para alimentador 49A EEQ con dos GD por
aproximación – ajuste.
127
Tabla 41. Resultados importados a Matlab para alimentador 49A EEQ con dos GD y
mejoramiento con capacitor - inductor mediante aproximación y ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del Capacitor 1 Fase A: nodo 72451MA nodo 50420MA
A una distancia desde la S/E de: 219604.70 [ft] 78482.70 [ft]
A una distancia desde la S/E de: 66935.51 [m] 23921.53 [m]
Con una reactancia de: -123.0000 [kVAr] -201.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 559072MA nodo 28144MA
A una distancia desde la S/E de: 124865.40 [ft] 94466.10 [ft]
A una distancia desde la S/E de: 38058.97 [m] 28793.27 [m]
Con una reactancia de: -159.0000 [kVAr] -184.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 80335MA nodo 8297MA
A una distancia desde la S/E de: 83456.70 [ft] 146096.00 [ft]
A una distancia desde la S/E de: 25437.60 [m] 44530.06 [m]
Con una reactancia de: 59.0000 [kVAr] 48.0000 [kVAr]
Pérdidas por potencia activa: 100.9746 [kW] 96.4830 [kW]
Porcentaje de disminución de pérdida: 57,3291% 59,2272%
Con una ganacia adicional de: -3,5148% -3,1487%
El porcentaje de pérdidas de potencia activa es: 1,7070% 1,6310%
128
El análisis estadístico de las Figuras 144, 154 y 156 para las fases ABC se indica en la
Figura 157:
Figura 157. Análisis estadístico de los perfiles de voltaje para alimentador 49A EEQ con dos GD y
mejora con capacitores e inductor por aproximación – ajuste.
129
Análisis de energía para 1GD y 2GD
Tabla 42. Resultados importados a Matlab para alimentador 49A con perfil de energía sin GD.
Variable Valor
Energia activa consumida: 6012.9303 [MWh]
Energía aparente consumida: 6159.1426 [MVAh]
Pérdidas por potencia activa: 201.8696 [MWh]
El porcentaje de pérdidas: 3,3573%
Figura 159. Perfiles de voltaje para alimentador 49A con perfil de energía y sin GD.
Tabla 43. Resultados importados a Matlab para alimentador 49A con perfil de energía y GD.
Variable solo G1 GD1 GD2
El mejor lugar del generador 1: 1031391MA 1031391MA 412501MA
A una distancia desde la S/E de: 139248.70 [ft] 139248.70 [ft] 94700.20 [ft]
A una distancia desde la S/E de: 42443.00 [m] 42443.00 [m] 28864.62 [m]
Con una potencia de: 1500.00 [kW] 900.0000 [kW] 500.0000 [kW]
Pérdida mínimizada a: 139.9367 [MWh] 148.7454 [MWh] 160.4083 [MWh]
Porcentaje de disminución de pérdida: 30,6797% 26,3161% 20,5387%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.9900 -0.9900 1.0000
Pérdida mínimizada a: 139.7462 [MWh] 131.0369 [MWh] 130.4032 [MWh]
Porcentaje de disminución de pérdida: 30,7740% 35,0883% 35,4023%
Con una ganacia adicional de: 0,0944% 8,7722% 14,8636%
Con respecto al generador solo con potencia activa
Total de pérdidas de energía por potencia: 2,3241% 2,1793% 2,1687%
Figura 160. Ajuste de potencia de GD para alimentador 49A con energía y dos GD.
130
a) Un GD b) Dos GD
Figura 161. Perfiles de voltaje para alimentador 49A con perfil de energía con GD.
Tabla 44. Resultados importados a Matlab para alimentador 49A con perfil de energía, con GD y
mejoramiento con capacitores.
Variable con G1 con G1 y G2
El mejor lugar del Capacitor 1 Fase A: 681301MA 50420MA
A una distancia desde la S/E de: 85957.60 [ft] 78482.70 [ft]
A una distancia desde la S/E de: 26199.88 [m] 23921.53 [m]
Con una reactancia de: -253.0000 [kVAr] -96.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: 70628MA 28144MA
A una distancia desde la S/E de: 58050.30 [ft] 94466.10 [ft]
A una distancia desde la S/E de: 17693.73 [m] 28793.27 [m]
Con una reactancia de: -275.0000 [kVAr] -548.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: 114357MA 8297MA
A una distancia desde la S/E de: 35400.60 [ft] 146096.00 [ft]
A una distancia desde la S/E de: 10790.10 [m] 44530.06 [m]
Con una reactancia de: -186.0000 [kVAr] -559.0000 [kVAr]
Pérdidas por energía: 237.3900 [MWh] 1150.6617 [MWh]
Porcentaje de disminución de pérdida: -17,5957% -470,0025%
Con una ganacia adicional solo con GD de: -48,3697% -500,7765%
Total de pérdidas de energía por potencia activa: 3,9480% 19,1365%
a) Un GD b) Dos GD
Figura 162. Lugar y reactancia de capacitores para alimentador 49A con perfil de energía, con GD
y FP.
131
a) Un GD b) Dos GD
Figura 163. Perfiles de voltaje para alimentador 49A con perfil de energía, con GD y mejorada con
capacitores.
El análisis estadístico de las Figuras 159, 161 y 163 para las fases ABC se indica en la
Figura 164:
Figura 164. Análisis estadístico de los perfiles de voltaje para alimentador 49A con perfiles de
energía, un GD y mejora con capacitores.
132
3.1.3. Redes de la Empresa Eléctrica Ambato EEASA.
Tabla 45. Resultados importados a Matlab para alimentador Misahualli EEASA sin GD.
Variable Valor
Potencia activa del sistema: 1258.4000 [kW]
Potencia aparente del sistema: 1321.2942 [kVA]
Factor de potencia: 0.9524
Pérdidas por potencia activa: 42.8250 [kW]
El porcentaje de pérdidas: 3,4031%
Figura 165. Perfiles de voltaje para alimentador Misahualli EEASA sin GD.
Tabla 46. Resultados importados a Matlab para alimentador Misahualli EEASA con un GD de
forma directa y por aproximación – ajuste.
133
a) Directo b) Aproximación – Ajuste
Figura 166. Lugar y potencia de GD para alimentador Misahualli EEASA con un GD de forma
directa y por aproximación – ajuste.
Tabla 47. Resultados importados a Matlab para alimentador Misahualli EEASA con un GD y
mejoramiento con capacitor – inductor de forma directa y por aproximación – ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del Capacitor 1 Fase A: nodo MTA_S_21176 nodo MTA_S_151997
A una distancia desde la S/E de: 61958.00 [ft] 69580.00 [ft]
A una distancia desde la S/E de: 18884.80 [m] 21207.98 [m]
Con una reactancia de: -53.0000 [kVAr] -49.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo MTA_S_290835 nodo MTA_S_739761
A una distancia desde la S/E de: 225502.30 [ft] 228778.50 [ft]
A una distancia desde la S/E de: 68733.10 [m] 69731.69 [m]
Con una reactancia de: -30.0000 [kVAr] -32.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo MTA_S_466156 nodo MTA_S_21976
A una distancia desde la S/E de: 155496.10 [ft] 99746.70 [ft]
A una distancia desde la S/E de: 47395.21 [m] 30402.79 [m]
Con una reactancia de: 20.0000 [kVAr] 27.0000 [kVAr]
Pérdidas por potencia activa: 6.5578 [kW] 6.5464 [kW]
Porcentaje de disminución de pérdida: 84,6870% 84,7136%
Con una ganacia adicional de: 0,2370% 0,1847%
El porcentaje de pérdidas de potencia activa es: 0,5211% 0,5202%
134
a) Directo b) Aproximación – Ajuste
Figura 168. Lugar y reactancia de capacitores - inductor para alimentador Misahualli EEASA con
un GD y FP de forma directa y por aproximación – ajuste.
El análisis estadístico de las Figuras 165, 167 y 169 para las fases ABC se indica en la
Figura 170:
135
Figura 170. Análisis estadístico de los perfiles de voltaje para alimentador Misahualli EEASA con
un GD y mejora con capacitores.
Tabla 48. Resultados importados a Matlab para alimentador Misahualli EEASA con dos GD de
forma directa y por aproximación – ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del generador 1: MTA_S_9361 MTA_S_22995
A una distancia desde la S/E de: 112477.90 [ft] 125868.60 [ft]
A una distancia desde la S/E de: 34283.26 [m] 38364.75 [m]
Con una potencia de: 900.0000 [kW] 920.0000 [kW]
Pérdida mínimizada a: 10.8473 [kW] 10.4022 [kW]
Porcentaje de disminución de pérdida: 74,6706% 75,7100%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.9700 -0.9700
Pérdida mínimizada a: 6.2521 [kW] 6.7051 [kW]
Porcentaje de disminución de pérdida: 85,4008% 84,3430%
Con una ganacia adicional de: 10,7302% 8,6330%
Con respecto al generador 1 solo con potencia activa
El mejor lugar del generador 2: MTA_S_739753 MTA_S_69444
A una distancia desde la S/E de: 226685.10 [ft] 243823.90 [ft]
A una distancia desde la S/E de: 69093.62 [m] 74317.52 [m]
Con una potencia de: 200.0000 [kW] 80.0000 [kW]
Pérdida mínimizada a: 8.2319 [kW] 8.6651 [kW]
Porcentaje de disminución de pérdida: 80,7778% 79,7663%
Con el ajuste del valor del factor de potencia al G2...
Valor calculado es: -0.9800 1.0000
Pérdida mínimizada a: 6.1629 [kW] 6.7051 [kW]
Porcentaje de disminución de pérdida: 85,6091% 84,3430%
136
Con una ganacia adicional de: 4,8313% 4,5768%
Con respecto al generador 2 solo con potencia activa
El porcentaje de pérdidas de potencia activa es: 0,4897% 0,5328%
Con una ganacia adicional de: 10,7302% 8,6330%
Con respecto al generador 1 solo con potencia activa
El mejor lugar del generador 2: MTA_S_739753 MTA_S_69444
A una distancia desde la S/E de: 226685.10 [ft] 243823.90 [ft]
A una distancia desde la S/E de: 69093.62 [m] 74317.52 [m]
Con una potencia de: 200.0000 [kW] 80.0000 [kW]
Pérdida mínimizada a: 8.2319 [kW] 8.6651 [kW]
Porcentaje de disminución de pérdida: 80,7778% 79,7663%
Con el ajuste del valor del factor de potencia al G2...
Valor calculado es: -0.9800 1.0000
Pérdida mínimizada a: 6.1629 [kW] 6.7051 [kW]
Porcentaje de disminución de pérdida: 85,6091% 84,3430%
Con una ganacia adicional de: 4,8313% 4,5768%
Con respecto al generador 2 solo con potencia activa
El porcentaje de pérdidas de potencia activa es: 0,4897% 0,5328%
Porcentaje de aporte de potencia de GD 87,4126% 79,4660%
137
a) Factor de Potencia = 1 b) Factor de Potencia = -0.97 y 1
Figura 175. Perfiles de voltaje para alimentador Misahualli EEASA con dos GD por aproximación
– ajuste.
Tabla 49. Resultados importados a Matlab para alimentador Misahualli EEASA con dos GD y
mejoramiento con capacitor - inductor mediante aproximación y ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del Capacitor 1 Fase A: MTA_S_651146 MTA_S_19963
A una distancia desde la S/E de: 57657.70 [ft] 55069.00 [ft]
A una distancia desde la S/E de: 17574.07 [m] 16785.03 [m]
Con una reactancia de: -21.0000 [kVAr] -84.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: MTA_S_650848 MTA_S_23446
A una distancia desde la S/E de: 85763.20 [ft] 24874.40 [ft]
A una distancia desde la S/E de: 26140.62 [m] 7581.72 [m]
Con una reactancia de: -29.0000 [kVAr] -41.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: MTA_S_21976 MTA_S_72092
A una distancia desde la S/E de: 99746.70 [ft] 130890.80 [ft]
A una distancia desde la S/E de: 30402.79 [m] 39895.52 [m]
Con una reactancia de: 41.0000 [kVAr] 39.0000 [kVAr]
Pérdidas por potencia activa: 5.9968 [kW] 6.7130 [kW]
Porcentaje de disminución de pérdida: 85,9970% 84,3246%
Con una ganacia adicional de: 0,3879% -0,0184%
El porcentaje de pérdidas de potencia activa es: 0,4765% 0,5335%
138
a) Directo b) Aproximación – Ajuste
Figura 177. Perfiles de voltaje para alimentador Misahualli EEASA con dos GD y mejorada con
capacitores e inductor de forma directa y por aproximación – ajuste.
El análisis estadístico de las Figuras 165, 175 y 177 para las fases ABC se indica en la
Figura 178:
Figura 178. Análisis estadístico de los perfiles de voltaje para alimentador Misahualli EEASA con
dos GD y mejora con capacitores e inductor por aproximación – ajuste.
139
Figura 179. Vista espacial alimentador Misahualli EEASA con dos GD y Capacitores.
Resultados importados a Matlab para alimentador Misahualli EEASA con perfil de energía
sin GD.
Tabla 50. Resultados importados a Matlab para alimentador Misahualli EEASA con perfil de
energía sin GD.
Variable Valor
Energia activa consumida: 2642.4087 [MWh]
Energía aparente consumida: 2718.4042 [MVAh]
Pérdidas por potencia activa: 28.9212 [MWh]
El porcentaje de pérdidas: 1,0945%
Figura 180. Perfiles de voltaje para alimentador Misahualli EEASA con perfil de energía y sin GD.
140
Tabla 51. Resultados importados a Matlab para alimentador Misahualli EEASA con perfil de
energía y GD.
Variable solo G1 GD1 GD2
El mejor lugar del generador 1: MTA_S_156760 MTA_S_54636 MTA_S_651532
A una distancia desde la S/E de: 126318.50 [ft] 113938.70 [ft] 244330.00 [ft]
A una distancia desde la S/E de: 38501.88 [m] 34728.52 [m] 74471.78 [m]
Con una potencia de: 600.0000 [kW] 300.0000 [kW] 100.0000 [kW]
Pérdida mínimizada a: 17.9071 [MWh] 21.1071 [MWh] 23.1771 [MWh]
Porcentaje de disminución de pérdida: 38,0831% 27,0186% 19,8612%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: 1.0000 1.0000 0.9800
Pérdida mínimizada a: 17.9071 [MWh] 19.0105 [MWh] 18.8668 [MWh]
Porcentaje de disminución de pérdida: 38,0831% 34,2679% 34,7648%
Con una ganacia adicional de: 0,0000% 7,2494% 14,9036%
Con respecto al generador solo con potencia activa
Total de pérdidas de energía por potencia: 0,6777% 0,7194% 0,7140%
Figura 181. Ajuste de potencia de GD para alimentador Misahualli EEASA con energía y un GD.
a) Un GD b) Dos GD
Figura 182. Perfiles de voltaje para alimentador Misahualli EEASA con perfil de energía y un GD.
141
Tabla 52. Resultados importados a Matlab para alimentador Misahualli EEASA con perfil de
energía, un GD y mejoramiento con capacitores.
Variable con G1 con G1 y G2
El mejor lugar del Capacitor 1 Fase A: MTA_S_151997 MTA_S_38941
A una distancia desde la S/E de: 69580.00 [ft] 107547.00 [ft]
A una distancia desde la S/E de: 21207.98 [m] 32780.33 [m]
Con una reactancia de: -12.0000 [kVAr] -168.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: MTA_S_739761 MTA_S_20888
A una distancia desde la S/E de: 228778.50 [ft] 66750.80 [ft]
A una distancia desde la S/E de: 69731.69 [m] 20345.64 [m]
Con una reactancia de: -66.0000 [kVAr] -281.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: MTA_S_21976 MTA_S_70695
A una distancia desde la S/E de: 99746.70 [ft] 117212.00 [ft]
A una distancia desde la S/E de: 30402.79 [m] 35726.22 [m]
Con una reactancia de: -55.0000 [kVAr] -251.0000 [kVAr]
Pérdidas por energía: 32.2241 [MWh] 180.8105 [MWh]
Porcentaje de disminución de pérdida: -11,4203% -525,1833%
Con una ganacia adicional solo con GD de: -49,5035% -563,2664%
Total de pérdidas de energía por potencia activa: 1,2195% 6,8426%
a) Un GD b) Dos GD
Figura 183. Lugar y reactancia de capacitores para alimentador Misahualli EEASA con perfil de
energía, un GD y FP.
a) Un GD b) Dos GD
Figura 184. Perfiles de voltaje para alimentador Misahualli EEASA con perfil de energía, un GD y
mejorada con capacitores.
142
El análisis estadístico de las Figuras 180, 182 y 184 para las fases ABC se indica en la
Figura 185:
Figura 185. Análisis estadístico de los perfiles de voltaje para alimentador Misahualli EEASA con
perfiles de energía, un GD y mejora con capacitores.
Tabla 53. Resultados importados a Matlab para alimentador La Maná ELEPCO sin GD.
Variable Valor
Potencia activa del sistema: 2515.1000 [kW]
Potencia aparente del sistema: 2797.4744 [kVA]
Factor de potencia: 0.8991
Pérdidas por potencia activa: 293.2266 [kW]
El porcentaje de pérdidas: 11,6586%
143
Figura 186. Perfiles de voltaje para alimentador La Maná ELEPCO sin GD.
Tabla 54. Resultados importados a Matlab para alimentador La Maná ELEPCO con un GD de
forma directa y por aproximación – ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del generador 1: MTA_S_201366 MTA_S_50366
A una distancia desde la S/E de: 64661.30 [ft] 61752.70 [ft]
A una distancia desde la S/E de: 19708.76 [m] 18822.22 [m]
Con una potencia de: 1500.0000 [kW] 1510.0000 [kW]
Pérdida mínimizada a: 221.1288 [kW] 218.8660 [kW]
Porcentaje de disminución de pérdida: 24,5877% 25,3594%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.9500 -0.9400
Pérdida mínimizada a: 210.9903 [kW] 209.1703 [kW]
Porcentaje de disminución de pérdida: 28,0453% 28,6660%
Con una ganacia adicional de: 3,4576% 3,3066%
Con respecto al generador solo con potencia activa
El porcentaje de pérdidas de potencia activa es: 8,3889% 8,3166%
Porcentaje de aporte de potencia de GD 59,6398% 60,0374%
144
a) Factor de Potencia = 1 b) Factor de Potencia = -0.94
Figura 188. Perfiles de voltaje para alimentador La Maná ELEPCO con un GD por aproximación –
ajuste.
Tabla 55. Resultados importados a Matlab para alimentador La Maná ELEPCO con un GD y
mejoramiento con capacitor – inductor de forma directa y por aproximación – ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del Capacitor 1 Fase A: MTA_S_76981 MTA_S_77029
A una distancia desde la S/E de: 121486.80 [ft] 124453.80 [ft]
A una distancia desde la S/E de: 37029.18 [m] 37933.52 [m]
Con una reactancia de: -127.0000 [kVAr] -125.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: MTA_S_49394 MTA_S_255610
A una distancia desde la S/E de: 25528.10 [ft] 29159.50 [ft]
A una distancia desde la S/E de: 7780.96 [m] 8887.82 [m]
Con una reactancia de: -189.0000 [kVAr] -162.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: MTA_S_99874 MTA_S_255043
A una distancia desde la S/E de: 116808.50 [ft] 113850.00 [ft]
A una distancia desde la S/E de: 35603.23 [m] 34701.48 [m]
Con una reactancia de: -114.0000 [kVAr] -86.0000 [kVAr]
Pérdidas por potencia activa: 210.5980 [kW] 208.9108 [kW]
Porcentaje de disminución de pérdida: 28,1791% 28,7545%
Con una ganacia adicional de: 0,1338% 0,0885%
El porcentaje de pérdidas de potencia activa es: 8,3733% 8,3063%
145
a) Directo b) Aproximación – Ajuste
Figura 190. Perfiles de voltaje para alimentador La Maná ELEPCO con un GD y mejorada con
capacitores de forma directa y por aproximación – ajuste.
El análisis estadístico de las Figuras 186, 188 y 190 para las fases ABC se indica en la
Figura 191:
Figura 191. Análisis estadístico de los perfiles de voltaje para alimentador La Maná ELEPCO con
un GD y mejora con capacitores.
146
Figura 192. Vista espacial alimentador La Maná ELEPCO con un GD y capacitores.
Tabla 56. Resultados importados a Matlab para alimentador La Maná ELEPCO con dos GD de
forma directa y por aproximación – ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del generador 1: MTA_S_50372 MTA_S_50366
A una distancia desde la S/E de: 57229.80 [ft] 61752.70 [ft]
A una distancia desde la S/E de: 17443.64 [m] 18822.22 [m]
Con una potencia de: 1500.0000 [kW] 1500.0000 [kW]
Pérdida mínimizada a: 220.2089 [kW] 218.8592 [kW]
Porcentaje de disminución de pérdida: 24,9015% 25,3618%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.9200 -0.9200
Pérdida mínimizada a: 194.6915 [kW] 194.5906 [kW]
Porcentaje de disminución de pérdida: 33,6037% 33,6381%
Con una ganacia adicional de: 8,7023% 8,2764%
Con respecto al generador 1 solo con potencia activa
El mejor lugar del generador 2: MTS_S_16325 MTA_S_76981
A una distancia desde la S/E de: 121995.00 [ft] 121486.80 [ft]
A una distancia desde la S/E de: 37184.08 [m] 37029.18 [m]
Con una potencia de: 400.0000 [kW] 350.0000 [kW]
Pérdida mínimizada a: 207.5383 [kW] 207.3921 [kW]
Porcentaje de disminución de pérdida: 29,2226% 29,2724%
Con el ajuste del valor del factor de potencia al G2...
Valor calculado es: -0.9300 -0.9200
Pérdida mínimizada a: 192.2721 [kW] 192.4842 [kW]
Porcentaje de disminución de pérdida: 34,4288% 34,3565%
Con una ganacia adicional de: 5,2063% 5,0841%
Con respecto al generador 2 solo con potencia activa
El porcentaje de pérdidas de potencia activa es: 7,6447% 7,6531%
Porcentaje de aporte de potencia de GD 75,5437% 73,5557%
147
a) Directo b) Aproximación – Ajuste
Figura 193. Lugar y potencia de GD para alimentador La Maná ELEPCO con dos GD de forma
directa y por aproximación – ajuste.
148
a) Factor de Potencia = 1 b) Factor de Potencia = -0.87 y -0.86
Figura 196. Perfiles de voltaje para alimentador La Maná ELEPCO con dos GD por aproximación
– ajuste.
Tabla 57. Resultados importados a Matlab para alimentador La Maná ELEPCO con dos GD y
mejoramiento con capacitor - inductor mediante aproximación y ajuste.
Variable Valor directo Valor ajustado
El mejor lugar del Capacitor 1 Fase A: MTA_S_256075 MTA_S_77242
A una distancia desde la S/E de: 124455.70 [ft] 139294.40 [ft]
A una distancia desde la S/E de: 37934.10 [m] 42456.93 [m]
Con una reactancia de: -40.0000 [kVAr] -38.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: MTA_S_76804 MTA_S_50371
A una distancia desde la S/E de: 68146.10 [ft] 56448.30 [ft]
A una distancia desde la S/E de: 20770.93 [m] 17205.44 [m]
Con una reactancia de: 40.0000 [kVAr] 41.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: MTA_S_323820 MTA_S_99717
A una distancia desde la S/E de: 114595.20 [ft] 123440.70 [ft]
A una distancia desde la S/E de: 34928.62 [m] 37624.73 [m]
Con una reactancia de: -69.0000 [kVAr] -65.0000 [kVAr]
Pérdidas por potencia activa: 192.1725 [kW] 192.4348 [kW]
Porcentaje de disminución de pérdida: 34,4628% 34,3733%
Con una ganacia adicional de: 0,0340% 0,0168%
El porcentaje de pérdidas de potencia activa es: 7,6407% 7,6512%
149
a) Directo b) Aproximación – Ajuste
Figura 198. Perfiles de voltaje para alimentador La Maná ELEPCO con dos GD y mejorada con
capacitores e inductor de forma directa y por aproximación – ajuste.
El análisis estadístico de las Figuras 186, 196 y 198 para las fases ABC se indica en la
Figura 199:
Figura 199. Análisis estadístico de los perfiles de voltaje para alimentador La Maná ELEPCO con
dos GD y mejora con capacitores e inductor por aproximación – ajuste.
150
Figura 200. Vista espacial alimentador La Maná ELEPCO con dos GD y Capacitores.
Tabla 58. Resultados importados a Matlab para alimentador La Maná ELEPCO con perfil de
energía sin GD.
Variable Valor
Energia activa consumida: 6158.0533 [MWh]
Energía aparente consumida: 6440.3157 [MVAh]
Pérdidas por potencia activa: 275.2695 [MWh]
El porcentaje de pérdidas: 4,4701%
Figura 201. Perfiles de voltaje para alimentador La Maná ELEPCO con perfil de energía y sin GD.
151
Tabla 59. Resultados importados a Matlab para alimentador La Maná ELEPCO con perfil de
energía y un GD.
Variable solo G1 GD1 GD2
El mejor lugar del generador 1: MTA_S_50366 MTA_S_50366 MTA_S_76981
A una distancia desde la S/E de: 61752.70 [ft] 61752.70 [ft] 121486.80 [ft]
A una distancia desde la S/E de: 18822.22 [m] 18822.22 [m] 37029.18 [m]
Con una potencia de: 1300.0000 [kW] 1100.0000 [kW] 300.0000 [kW]
Pérdida mínimizada a: 233.6439 [MWh] 234.3955 [MWh] 263.0487 [MWh]
Porcentaje de disminución de pérdida: 15,1218% 14,8487% 4,4396%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.9800 -0.9700 -0.9800
Pérdida mínimizada a: 231.6241 [MWh] 219.5583 [MWh] 219.0022 [MWh]
Porcentaje de disminución de pérdida: 15,8555% 20,2388% 20,4408%
Con una ganacia adicional de: 0,7338% 5,3901% 16,0012%
Con respecto al generador solo con potencia activa
Total de pérdidas de energía por potencia: 3,7613% 3,5654% 3,5564%
Figura 202. Ajuste de potencia de GD para alimentador La Maná ELEPCO con energía y un GD.
a) Un GD b) Dos GD
Figura 203. Perfiles de voltaje para alimentador La Maná ELEPCO con perfil de energía y un GD.
152
Tabla 60. Resultados importados a Matlab para alimentador La Maná ELEPCO con perfil de
energía, un GD y mejoramiento con capacitores.
Variable con G1 con G1 y G2
El mejor lugar del Capacitor 1 Fase A: MTA_S_77029 MTA_S_77242
A una distancia desde la S/E de: 124453.80 [ft] 139294.40 [ft]
A una distancia desde la S/E de: 37933.52 [m] 42456.93 [m]
Con una reactancia de: -23.0000 [kVAr] -561.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: MTA_S_255610 MTA_S_50371
A una distancia desde la S/E de: 29159.50 [ft] 56448.30 [ft]
A una distancia desde la S/E de: 8887.82 [m] 17205.44 [m]
Con una reactancia de: -19.0000 [kVAr] -602.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: MTA_S_255043 MTA_S_99717
A una distancia desde la S/E de: 113613.30 [ft] 123440.70 [ft]
A una distancia desde la S/E de: 34629.33 [m] 37624.73 [m]
Con una reactancia de: -21.0000 [kVAr] -222.0000 [kVAr]
Pérdidas por energía: 230.7357 [MWh] 2657.0629 [MWh]
Porcentaje de disminución de pérdida: 16,1783% -865,2587%
Con una ganacia adicional solo con GD de: 0,3227% -881,1143%
Total de pérdidas de energía por potencia activa: 3,7469% 43,1478%
a) Un GD b) Dos GD
Figura 204. Lugar y reactancia de capacitores para alimentador La Maná ELEPCO con perfil de
energía, un GD y FP.
a) Un GD b) Dos GD
Figura 205. Perfiles de voltaje para alimentador La Maná ELEPCO con perfil de energía, un GD y
mejorada con capacitores.
153
El análisis estadístico de las Figuras 201, 203 y 205 para las fases ABC se indica en la
Figura 206:
Figura 206. Análisis estadístico de los perfiles de voltaje para alimentador La Maná ELEPCO con
perfiles de energía, un GD y mejora con capacitores.
3.2. Discusión
154
Tabla 61. Registro de simulaciones para red de prueba 13 Nodos IEEE.
Casos Tiempo [s] Potencia [kW] Val min % Val min fp opt % Nodo
generados 1fp
100 11,2953 3000,0000 55,10 56,4045 27,4792 78,2595 692
200 21,4087 3100,0000 55,41 56,1642 27,4807 78,2584 671
225 24,7137 2900,0000 54,97 56,5119 27,6475 78,1264 671
250 27,6288 2900,0000 54,97 56,5119 27,6475 78,1264 671
275 30,4328 2800,0000 55,00 56,4858 27,9801 77,8632 671
300 31,1516 3000,0000 55,10 56,4045 27,4792 78,2595 671
400 40,6941 2900,0000 54,97 56,5119 27,6475 78,1264 692
Tabla 64. Registro de simulaciones para red de prueba 123 Nodos IEEE.
Casos Tiempo [s] Potencia Val min % Val min fp opt % Nodo
generado [kW] 1fp
s
800 84,7954 2600,0000 48,4249 55,6217 28,0114 74,3294 67
900 98,1622 2100,0000 48,6106 55,4516 29,1432 73,2922 72
1000 113,6847 2600,0000 48,4249 55,6217 28,0114 74,3294 67
2000 238,8309 2300,0000 47,3090 56,6444 27,3874 74,9012 67
2500 315,5281 2300,0000 47,3090 56,6444 27,3874 74,9012 67
3000 366,3885 2300,0000 47,3090 56,6444 27,3874 74,9012 67
4000 520,5660 2300,0000 47,3090 56,6444 27,3874 74,9012 67
Como se observa en las tablas 61, 63 y 64 se establece que para un cálculo directo
(sección 2.4.4), por cada nodo de la red se deben generar al menos 20 simulaciones para
obtener un resultado confiable y que por cada segundo se realizan 10 simulaciones. El
método de aproximación y ajuste (sección 2.4.5) es de gran ayuda para aumentar la
155
precisión de resultados; así, por ejemplo, para un alimentador de 1.000 nodos se requieren
al menos 20.000 casos de simulación en un método de cálculo directo y en el cálculo por
aproximación y ajuste se requieren al menos 200 casos. Esto se consigue con la instrucción
de redondeo de números “round” en el que aproxima las cantidades de centenas a
unidades para posiciones de GD y de centenas a decenas para potencias de GD. De los
resultados obtenidos en los alimentadores de EEQ, EEASA y ELEPCO (tablas 38, 40, 46,
48, 54 y 56) la diferencia de reducción de pérdidas entre el método directo y por
aproximación – ajuste es mínima, lo que asegura que la cantidad de simulaciones fue la
adecuada para obtener los resultados esperados.
Tabla 65. Comparación de resultados importados a Matlab para red modificada de 123 nodos
IEEE con un GD.
156
a) Bivariada b) Trivariada
Figura 207. Comparación de perfiles de voltaje para la red modificada de 123 nodos IEEE con un
GD.
Tabla 66. Comparación de resultados importados a Matlab para red modificada de 123 nodos
IEEE con un GD y mejoramiento con capacitor - inductor.
Variable Bivariado Trivariado
El mejor lugar del Capacitor 1 Fase A: nodo 18 nodo 135
A una distancia desde la S/E de: 2025.00 [ft] 2025.00 [ft]
A una distancia desde la S/E de: 617.22 [m] 617.22 [m]
Con una reactancia de: -401.00 [kVAr] -394.0000 [kVAr]
El mejor lugar del Inductor 1 Fase B: nodo 97 nodo 81
A una distancia desde la S/E de: 3627.00 [ft] 5302.00 [ft]
A una distancia desde la S/E de: 1105.51 [m] 1616.05 [m]
Con una reactancia de: 72.00 [kVAr] 87.0000 [kVAr]
El mejor lugar del Capacitor 1 Fase C: nodo 34 nodo 50
A una distancia desde la S/E de: 1350.00 [ft] 3851.00 [ft]
A una distancia desde la S/E de: 411.48 [m] 1173.78 [m]
Con una reactancia de: -293.00 [kVAr] -183.0000 [kVAr]
Pérdidas por potencia activa: 26.4380 [kW] 26.1413 [kW]
Porcentaje de disminución de pérdida: 75.7713 % 76,0432%
Con una ganacia adicional de: 1.08 % 1,17%
El porcentaje de pérdidas de potencia activa es: 0.7681 % 0.7594 %
a) Bivariada b) Trivariada
Figura 208. Comparación de perfiles de voltaje para la red modificada de 123 nodos IEEE con un
GD y mejorada con capacitores – inductor.
157
El análisis estadístico de las Figuras 103, 207 y 208 para las fases ABC se indica en la
Figura 209:
Figura 209. Análisis estadístico de los perfiles de voltaje para la red modificada de 123 nodos
IEEE con un GD y mejora con capacitores.
Por lo observado en los resultados de las tablas 21, 24, 26 y 28 para la red de prueba de
123 nodos de la IEEE, el aumento de lugares para la inserción de GD reduce más las
pérdidas técnicas. Para un GD reducen en 1.3743%, con dos GD en 1.1101% y con tres
GD en 1.0554%, con mejoras del 56.65%, 64.98% y 66.71% respectivamente, además los
perfiles de voltaje en promedio y varianza mejoraron considerablemente según las Figuras
158
113, 120 y 127 para un, dos y tres GD respectivamente. La diferencia entre dos y tres GD
es poca, por lo que es práctico aplicar este método solo con uno y dos GD para el análisis
de alimentadores reales.
Como mejor opción para reducción de pérdidas, todos los resultados donde se ubican los
GD son trifásicos; por lo tanto, la variación del factor de potencia de la GD incide
directamente en el comportamiento de las 3 fases como lo indica las Figuras 105, 110, 116,
117, 123 y 124. En algunos casos una de las fases tiene un voltaje mayor al 1 por unidad,
situación poco común en alimentadores. Sin embargo, el ajuste del factor de potencia
ayuda en la reducción de las pérdidas técnicas como lo demuestra los resultados de las
tablas 113, 120 y 127 donde los valores a considerar son: 0.8024% para un GD, 0.4360%
con dos GD en y 0.3659% con tres GD, con mejoras del 74.89%, 83.81% y 88.86%
respectivamente. La variación del factor de potencia en unidades de generación y fuentes
de potencia es posible con conversores de electrónica de potencia AC – AC como lo
indicado en el Anexo E.
Para verificar la eficacia de la inserción de capacitores en la red, sin la variación del factor
de potencia (FP) en las nuevas unidades de GD, se aplica el mismo método en la red de
prueba de 123 IEEE con los resultados indicados en la Figura 210:
159
Como se puede verificar en el gráfico 210, la variación del factor de potencia de las
unidades de GD es más eficaz para el mejoramiento de la red. Por lo tanto, se utilizará la
variación de factor de potencia de las unidades de GD en los análisis de los alimentadores
de las distribuidoras.
La generación de casos aleatorios para cada capacitor, es únicamente para los lugares
factibles de conexión y depende de la topología de la red como lo indica las Figuras 106,
111, 118 y 125, donde se observa claramente los espacios verticales en blanco donde no
es factible la conexión de un capacitor monofásico. La variación de reducción de pérdidas
por inserción de capacitores es mínima como lo indica las Tablas 23, 25, 27 y 29.
Al comparar los resultados de ubicación y capacidad de cada capacitor y/o inductor de las
tablas 23, 25, 27 y 29 se observa una variación considerable entre el método directo y por
aproximación – ajuste a pesar de que se trata de una red de 123 nodos. Esto se debe a
que una mínima variación en la ubicación, potencia y factor de potencia de los GD, incide
de manera sensible en los resultados de los capacitores y/o inductores. La repetibilidad de
los resultados en capacitores y/o inductores tiene el mismo valor de confianza que con
generadores ya que se utiliza la misma metodología.
Los resultados de análisis de energía para capacitores indicados en las tablas 32, 34 y 36
tienen una diferente interpretación ya que los valores de potencia reactiva tienen un valor
superior en comparación a los resultados de las tablas 25, 27 y 29. Esto se debe a que el
factor de potencia de la red varía en función de la curva de perfil de demanda y en esta
metodología para el cálculo del tamaño del capacitor solo se está utilizando el valor
promedio de voltaje de la red más cercano al 1 p.u. con la menor varianza posible.
160
Esta metodología con capacitores es efectiva para análisis de potencia y no de energía;
por lo tanto, los valores de reactancia capacitiva (-) deberán ser los resultados de los
análisis de potencia y en caso que los resultados fueran con reactancia inductiva (+) es
preferible no instalar ningún reactor o inductor.
Aplicar estos criterios expuestos en los párrafos anteriores a los alimentadores reales, es
complejo y depende de: el recurso energético disponible, la longitud total de la red, número
de nodos, posición geográfica, la configuración de la red por topología, las pérdidas
técnicas, las caídas de voltaje, el desbalance entre fases y la carga instalada.
Para el análisis energético, los perfiles de voltaje de la Figura 159 son similares la Figura
144 pero a una escala diferente. Las potencias de operación promedio de los GD de la
Tabla 43 se reducen en comparación a la potencia instalada de los GD de la Tabla 38 y
40. La variación del factor de potencia de los GD es mínima entre las tablas 43, 38 y 40.
Hay un mejoramiento de perfiles de voltaje promedio de las Figuras 161 en relación a la
Figura 159. Los valores de potencia reactiva registrados en la Tabla 44 para dos GD están
fuera del rango de perfiles de voltaje permitido de acuerdo a la Figura 163b y por eso los
niveles de pérdidas de energía se incrementarían, siendo una mejor opción utilizar los
valores de potencia reactiva de las Tablas 39 y 41 solo para capacitores.
161
Análisis del alimentador Misahualli Empresa Eléctrica Ambato
Para el caso de la EEASA, la zona Oriental del cantón Tena es muy lejana y distante, con
una distancia de hasta 95 km desde la cabecera del alimentador hasta el punto más lejano,
opera a 13.8 kV y las caídas de voltaje llegan hasta en 0.92 p.u. en la fase A como lo
evidencia la Figura 165. Este es un alimentador ideal para aplicar GD. Al comparar los
resultados con un GD y dos GD, las pérdidas se reducen a 0.5265% para un GD (Tabla
46) y 0.5328% para dos GD (Tabla 48). La potencia instalada de un GD (Tabla 46) es del
65.95% y dos GD (Tabla 48) es 79.46% de la potencia instalada de la red. Al incorporar los
capacitores mejoran los perfiles de voltaje como lo indica la Figura 169 para un GD y la
Figura 177 para dos GD manteniendo los mismos niveles de pérdidas en 0,5211% (Tabla
47) y 0,5202% (Tabla 49) respectivamente. Las opciones de uno o dos GD son adecuadas
ya que están dentro de los rangos de tolerancia, aunque con dos GD mejora más los
perfiles de voltaje.
Para el análisis energético, los perfiles de voltaje de la Figura 165 son similares la Figura
180 pero a una escala diferente y dentro del rango de tolerancia por voltaje mínimo
permitido. Las potencias de operación promedio de los GD de la Tabla 51 se reducen en
comparación a la potencia instalada de los GD de la Tabla 46 y 48. La variación del factor
de potencia de los GD es hasta 0.03 entre las tablas 51, 46 y 48. Hay un mejoramiento de
perfiles de voltaje promedio de las Figuras 182 en relación a la Figura 180, siendo la fase
B la más favorecida para un GD. Los valores de potencia reactiva registrados en la Tabla
52 para dos GD están fuera del rango de perfiles de voltaje permitido de acuerdo a la Figura
184b y por eso los niveles de sobrevoltajes y pérdidas de energía se incrementarían
considerablemente, siendo una mejor opción utilizar los valores de potencia reactiva de las
Tablas 47 y 49 solo para capacitores.
162
dos GD manteniendo los mismos niveles de pérdidas en 8,3063% (Tabla 55) y 7,6512%
(Tabla 57) respectivamente. La opción de dos GD es más conveniente ya que posee
mejores perfiles de voltaje (Figura 198.b) en relación a la de un GD (Figura 190.b) aunque
todavía están fuera de los rangos de tolerancia. Con el incremento de puntos de GD se
reconfigura la red, disminuyen las pérdidas y los capacitores que mejoran los perfiles de
voltaje de la red. Para el caso de tres GD el mejoramiento de los perfiles de voltaje es
importante pero no tiene un cambio significativo en relación a dos GD, como se lo indica
en la Figura 211.
a) Dos GD b) Tres GD
Figura 211. Comparación de perfiles de voltaje para el alimentador La Maná con GD y mejorada
con capacitores.
Para el análisis energético, los perfiles de voltaje de la Figura 201 son similares la Figura
186 pero a una escala diferente y fuera del rango de tolerancia por voltaje mínimo permitido.
Las potencias de operación promedio de los GD de la Tabla 59 se reducen en comparación
a la potencia instalada de los GD de la Tabla 54 y 56. La variación del factor de potencia
de los GD es hasta 0.06 entre las tablas 59, 54 y 56. Hay un mejoramiento de perfiles de
voltaje promedio de las Figuras 203 en relación a la Figura 186, sin mayor diferencia entre
un GD y dos GD. Los valores de potencia reactiva registrados en la Tabla 60 para dos GD
están fuera del rango de perfiles de voltaje permitido de acuerdo a la Figura 205b y por eso
los niveles de sobrevoltajes, y pérdidas de energía se incrementarían considerablemente,
siendo una mejor opción utilizar los valores de potencia reactiva de la Tabla 57 solo para
capacitores.
163
3.2.4. Comparación de resultados con CYME.
A fin de que los resultados de esta metodología puedan ser verificados en el software que
utilizan las empresas de distribución en el Ecuador para modelación de redes de potencia,
en el Anexo C se presenta los niveles de precisión de cálculos entre CYME y la IEEE para
la red de prueba de 123 nodos, con resultados adecuados. Con la incorporación de GD a
la red de 123 nodos se obtiene los mismos niveles de precisión como se indica en las desde
la Tabla 67 hasta la Tabla 72.
Tabla 67. Resultado en CYME para red modificada de 123 nodos IEEE con un GD.
Tabla 68. Resultado en CYME para red modificada de 123 nodos IEEE con un GD y fp=–0.87.
164
Al comparar los resultados de las Tablas 67 y 68 con la Tabla 24 hay una diferencia de
precisión del 0.0162% y 0.0107% con fp=–0.87 por reducción de pérdidas.
Tabla 69. Resultado en CYME para red modificada de 123 nodos IEEE con dos GD.
Tabla 70. Resultado en CYME para red modificada de 123 nodos IEEE con dos GD y fp=–0.86.
Al comparar los resultados de las tablas 69 y 70 con la Tabla 26 hay una diferencia de
precisión del 0.0344% y 0.0161% con fp=–0.86 por reducción de pérdidas.
165
Tabla 71. Resultado en CYME para red modificada de 123 nodos IEEE con tres GD.
Tabla 72. Resultado en CYME para red modificada de 123 nodos IEEE con tres GD y fp=–0.85.
Al comparar los resultados de las Tablas 71 y 72 con la Tabla 28 hay una diferencia de
precisión del 0.0329% y 0.0086% con fp=–0.85 por reducción de pérdidas.
166
3.2.5. Comparación del mejoramiento de red de prueba de 37 Nodos IEEE
con otro método.
Los perfiles de voltaje presentados en [35] sin GD como lo indica la Figura 212.a no
corresponde con exactitud a los resultados presentados en [56] por lo que en la Figura
218.b se presenta el valor corregido.
El autor de [35] realiza el análisis de los perfiles de voltaje para 3 escenarios como lo indica
las Figuras 213 cuyo orden de eje horizontal es por nombre de nodo.
Figura 213. Perfiles de voltaje para red de 123 nodos IEEE con GD.
Como se observa en las Figuras 213 los voltajes están por encima de 1 p.u. y los nodos de
conexión de GD se los indica en la Tabla 73.
Tabla 73. Esquema de candidatos para los tipos, lugar y tamaño de GDs [35].
167
Al aplicar la metodología probabilística bivariada a la red de 37 nodos IEEE, con los mismos
criterios utilizados en la sección 3.1.1, se obtiene los resultados indicados en las Figuras
214, 215 y Tabla 74.
a) sin GD b) con GD
Figura 214. Perfiles de voltaje para red modificada de 37 nodos IEEE.
Figura 215. Análisis estadístico de los perfiles de voltaje para la red modificada de 37 nodos
IEEE, un GD y mejora con capacitores.
168
Tabla 74. Resultados para red modificada 37 nodos IEEE con GD.
Variable Valor
Potencia activa del sistema: 2326.2000 [kW]
Potencia aparente del sistema: 2696.9130 [kVA]
Factor de potencia: 0.8625
Pérdidas por potencia activa: 117.8369 [kW]
El porcentaje de pérdidas: 5,0656%
El mejor lugar del generador 1: nodo 702
A una distancia desde la S/E de: 2810.00 [ft]
A una distancia desde la S/E de: 856.49 [m]
Con una potencia de: 2200.00 [kW]
Pérdida mínimizada a: 45.3956 [kW]
Porcentaje de disminución de pérdida: 61,4759%
Porcentaje pérdidas potencia activa: 1,9515%
Con el ajuste del valor del factor de potencia al G1...
Valor calculado es: -0.8800
Pérdida mínimizada a: 26.7173 [kW]
Porcentaje de disminución de pérdida: 77,3269%
Con una ganacia adicional de: 15,8510%
Con respecto al generador solo con potencia activa
Porcentaje pérdidas potencia activa: 1,1485%
Aporte de potencia de GD 94,5748%
Al comparar los perfiles de voltaje de [35] con esta metodología para la red de 37 nodos
IEEE, se comprueba que esta metodología es mejor ya que los promedios de cada voltaje
son más cercanos a la unidad (Figura 215) y la varianza es la más cercana a cero. También
es mejor la presentación de perfiles de voltaje con un orden descendente de distancia en
el eje horizontal ya que de esa forma se puede apreciar las caídas de voltaje según la
distancia hacia la fuente, como se lo indica en la Figura 216.
a) sin GD b) con GD
Figura 216. Perfiles de voltaje para la red modificada 37 nodos IEEE con GD.
169
4. CONCLUSIONES
• Los resultados de los análisis de flujos de potencia son como una referencia inicial del
comportamiento de la red ante la inserción de GD, un análisis más detallado es al
establecer los valores de demanda máxima, mínima y media de cada alimentador para
compáralo con un análisis de energía y establecer el valor más adecuado de potencia
de GD.
• Aunque los resultados de análisis de flujos de potencia para GD son una referencia del
máximo valor posible de potencia de nueva generación a ser instalada en la red, los
resultados de potencia de GD por análisis de energía son hasta un 35% menor, siendo
un valor más cercano a la realidad de la operación de la red, tomando en cuenta que no
siempre toda la demanda del recurso energético renovable estará disponible, según la
región geográfica a ser analizada.
170
de las redes IEEE para demostrar la eficiencia de la inserción de GD y aplicar este
mismo criterio para el cálculo de GD en los alimentadores de la EEQ, EEASA y
ELEPCO.
• En los resultados y análisis de la red de 123 nodos IEEE (sección. 3.1.1 y 3.2.3) la
solución numérica para mejoramiento de perfiles de voltaje es con capacitores y
reactores; sin embargo, en alimentadores reales, no se emplea los reactores ya que la
característica de las cargas en la red es inductiva y se necesita inyectar reactivos
capacitivos. Por lo tanto, para respuestas con reactivos inductivos no deberán ser
tomadas en cuenta.
• A medida que se incrementa los puntos de GD, disminuye la potencia instalada en cada
punto, los promedios de voltaje se aproximan a 1 p.u y la solución teórica es un GD por
cada carga en cuyo caso sería un abastecimiento autónomo y la red de distribución
perdería su razón de ser.
• Cuando el recurso energético no sea suficiente para abastecer como GD, la distribuidora
tiene como opción la instalación de un transformador de potencia a través de una línea
de subtransmisión que eléctricamente reemplazaría al GD en posición y potencia y con
los equipos de electrónica de potencia podrá ajustar el factor de potencia de salida a la
red.
171
Asi por ejemplo, se lo demuestra en la Figura 211 donde no existe una variación
significativa entre dos GD y tres GD para el alimentador La Maná de ELEPCO
RECOMENDACIONES:
172
para el número de vectores aleatorios a ser analizados; por lo tanto, se deberá tomar en
cuenta todas las características de hardware para la aplicación en redes de distribución
de gran tamaño como las que fueron analizadas en esta tesis.
173
5. REFERENCIAS BIBLIOGRÁFICAS
[1] H. L. Willis y W. G. Scott, Distributed power generation: planning and evaluation. New York:
Marcel Dekker, 2000.
[2] G. Carpinelli, G. Celli, S. Mocci, F. Pilo, y A. Russo, «Optimisation of embedded generation
sizing and siting by using a double trade-off method», IEE Proc. - Gener. Transm. Distrib., vol.
152, n.o 4, p. 503, 2005, doi: 10.1049/ip-gtd:20045129.
[3] H. L. Willis, «Analytical methods and rules of thumb for modeling DG-distribution
interaction», 2000, vol. 3, pp. 1643-1644, doi: 10.1109/PESS.2000.868774.
[4] D. Q. Hung, N. Mithulananthan, y R. C. Bansal, «Analytical Expressions for DG Allocation in
Primary Distribution Networks», IEEE Trans. Energy Convers., vol. 25, n.o 3, pp. 814-820, sep.
2010, doi: 10.1109/TEC.2010.2044414.
[5] P. M. Costa y M. A. Matos, «Avoided losses on LV networks as a result of microgeneration»,
Electr. Power Syst. Res., vol. 79, n.o 4, pp. 629-634, abr. 2009, doi:
10.1016/j.epsr.2008.09.014.
[6] D. Q. Hung y N. Mithulananthan, «Multiple Distributed Generator Placement in Primary
Distribution Networks for Loss Reduction», IEEE Trans. Ind. Electron., vol. 60, n.o 4, pp. 1700-
1708, abr. 2013, doi: 10.1109/TIE.2011.2112316.
[7] N. S. Rau y Yih-Heui Wan, «Optimum location of resources in distributed planning», IEEE
Trans. Power Syst., vol. 9, n.o 4, pp. 2014-2020, nov. 1994, doi: 10.1109/59.331463.
[8] P. N. Vovos, G. P. Harrison, A. R. Wallace, y J. W. Bialek, «Optimal Power Flow as a Tool for
Fault Level-Constrained Network Capacity Analysis», IEEE Trans. Power Syst., vol. 20, n.o 2, pp.
734-741, may 2005, doi: 10.1109/TPWRS.2005.846070.
[9] A. Keane y M. O’Malley, «Optimal Allocation of Embedded Generation on Distribution
Networks», IEEE Trans. Power Syst., vol. 20, n.o 3, pp. 1640-1646, ago. 2005, doi:
10.1109/TPWRS.2005.852115.
[10] A. Keane y M. O’Malley, «Optimal Utilization of Distribution Networks for Energy Harvesting»,
IEEE Trans. Power Syst., vol. 22, n.o 1, pp. 467-475, feb. 2007, doi:
10.1109/TPWRS.2006.888959.
[11] M. F. AlHajri, M. R. AlRashidi, y M. E. El-Hawary, «Improved Sequential Quadratic
Programming Approach for Optimal Distribution Generation Deployments via Stability and
Sensitivity Analyses», Electr. Power Compon. Syst., vol. 38, n.o 14, pp. 1595-1614, dic. 2010,
doi: 10.1080/15325008.2010.492451.
[12] P. N. Vovos y J. W. Bialek, «Direct Incorporation of Fault Level Constraints in Optimal Power
Flow as a Tool for Network Capacity Analysis», IEEE Trans. Power Syst., vol. 20, n.o 4, pp.
2125-2134, nov. 2005, doi: 10.1109/TPWRS.2005.856975.
[13] Y. M. Atwa y E. F. El-Saadany, «Probabilistic approach for optimal allocation of wind-based
distributed generation in distribution systems», IET Renew. Power Gener., vol. 5, n.o 1, p. 79,
2011, doi: 10.1049/iet-rpg.2009.0011.
[14] Y. M. Atwa, E. F. El-Saadany, M. M. A. Salama, y R. Seethapathy, «Optimal Renewable
Resources Mix for Distribution System Energy Loss Minimization», IEEE Trans. Power Syst.,
vol. 25, n.o 1, pp. 360-370, feb. 2010, doi: 10.1109/TPWRS.2009.2030276.
[15] L. F. Ochoa, C. J. Dent, y G. P. Harrison, «Distribution Network Capacity Assessment: Variable
DG and Active Networks», p. 9.
[16] C. J. Dent, L. F. Ochoa, y G. P. Harrison, «Network Distributed Generation Capacity Analysis
Using OPF With Voltage Step Constraints», IEEE Trans. Power Syst., vol. 25, n.o 1, pp. 296-304,
feb. 2010, doi: 10.1109/TPWRS.2009.2030424.
174
[17] L. F. Ochoa y G. P. Harrison, «Minimizing Energy Losses: Optimal Accommodation and Smart
Operation of Renewable Distributed Generation», IEEE Trans. Power Syst., vol. 26, n.o 1, pp.
198-205, feb. 2011, doi: 10.1109/TPWRS.2010.2049036.
[18] S. Porkar, P. Poure, A. Abbaspour-Tehrani-fard, y S. Saadate, «Optimal allocation of
distributed generation using a two-stage multi-objective mixed-integer-nonlinear
programming», Eur. Trans. Electr. Power, vol. 21, n.o 1, pp. 1072-1087, ene. 2011, doi:
10.1002/etep.497.
[19] A. Kumar y W. Gao, «Optimal distributed generation location using mixed integer non-linear
programming in hybrid electricity markets», IET Gener. Transm. Distrib., vol. 4, n.o 2, p. 281,
2010, doi: 10.1049/iet-gtd.2009.0026.
[20] R. S. Al Abri, E. F. El-Saadany, y Y. M. Atwa, «Optimal Placement and Sizing Method to
Improve the Voltage Stability Margin in a Distribution System Using Distributed Generation»,
IEEE Trans. Power Syst., vol. 28, n.o 1, pp. 326-334, feb. 2013, doi:
10.1109/TPWRS.2012.2200049.
[21] N. Khalesi, N. Rezaei, y M.-R. Haghifam, «DG allocation with application of dynamic
programming for loss reduction and reliability improvement», Int. J. Electr. Power Energy
Syst., vol. 33, n.o 2, pp. 288-295, feb. 2011, doi: 10.1016/j.ijepes.2010.08.024.
[22] D. Zhu, R. P. Broadwater, K. Tam, R. Seguin, y H. Asgeirsson, «Impact of DG Placement on
Reliability and Efficiency With Time-Varying Loads», IEEE Trans. Power Syst., vol. 21, n.o 1, pp.
419-427, feb. 2006, doi: 10.1109/TPWRS.2005.860943.
[23] D. Singh, R. K. Misra, y D. Singh, «Effect of Load Models in Distributed Generation Planning»,
IEEE Trans. Power Syst., vol. 22, n.o 4, pp. 2204-2212, nov. 2007, doi:
10.1109/TPWRS.2007.907582.
[24] S. Kotamarty, S. Khushalani, y N. Schulz, «Impact of distributed generation on distribution
contingency analysis», Electr. Power Syst. Res., vol. 78, n.o 9, pp. 1537-1545, sep. 2008, doi:
10.1016/j.epsr.2008.01.020.
[25] H. Khan y M. A. Choudhry, «Implementation of Distributed Generation (IDG) algorithm for
performance enhancement of distribution feeder under extreme load growth», Int. J. Electr.
Power Energy Syst., vol. 32, n.o 9, pp. 985-997, nov. 2010, doi: 10.1016/j.ijepes.2010.02.006.
[26] J. O. Kim, S. K. Park, K. W. Park, y C. Singh, «ispersed Generation Planning Using Improved
Hereford Ranch Algorithm», p. 6.
[27] C. L. T. Borges y D. M. Falcão, «Optimal distributed generation allocation for reliability, losses,
and voltage improvement», Int. J. Electr. Power Energy Syst., vol. 28, n.o 6, pp. 413-420, jul.
2006, doi: 10.1016/j.ijepes.2006.02.003.
[28] T. N. Shukla, S. P. Singh, V. Srinivasarao, y K. B. Naik, «Optimal Sizing of Distributed
Generation Placed on Radial Distribution Systems», Electr. Power Compon. Syst., vol. 38, n.o 3,
pp. 260-274, ene. 2010, doi: 10.1080/15325000903273403.
[29] D. Singh, D. Singh, y K. S. Verma, «Multiobjective Optimization for DG Planning With Load
Models», IEEE Trans. Power Syst., vol. 24, n.o 1, pp. 427-436, feb. 2009, doi:
10.1109/TPWRS.2008.2009483.
[30] R. K. Singh y S. K. Goswami, «Optimum Siting and Sizing of Distributed Generations in Radial
and Networked Systems», Electr. Power Compon. Syst., vol. 37, n.o 2, pp. 127-145, ene. 2009,
doi: 10.1080/15325000802388633.
[31] R. K. Singh y S. K. Goswami, «Optimum allocation of distributed generations based on nodal
pricing for profit, loss reduction, and voltage improvement including voltage rise issue», Int. J.
Electr. Power Energy Syst., vol. 32, n.o 6, pp. 637-644, jul. 2010, doi:
10.1016/j.ijepes.2009.11.021.
[32] J.-H. Teng, Y.-H. Liu, C.-Y. Chen, y C.-F. Chen, «Value-based distributed generator placements
for service quality improvements», Int. J. Electr. Power Energy Syst., vol. 29, n.o 3, pp. 268-
274, mar. 2007, doi: 10.1016/j.ijepes.2006.07.008.
175
[33] M. Gandomkar, M. Vakilian, y M. Ehsan, «A Genetic–Based Tabu Search Algorithm for
Optimal DG Allocation in Distribution Networks», Electr. Power Compon. Syst., vol. 33, n.o 12,
pp. 1351-1362, dic. 2005, doi: 10.1080/15325000590964254.
[34] E. Haesen, J. Driesen, y R. Belmans, «Robust planning methodology for integration of
stochastic generators in distribution grids», IET Renew. Power Gener., vol. 1, n.o 1, p. 25,
2007, doi: 10.1049/iet-rpg:20060012.
[35] Z. Liu, F. Wen, y G. Ledwich, «Optimal Siting and Sizing of Distributed Generators in
Distribution Systems Considering Uncertainties», IEEE Trans. Power Deliv., vol. 26, n.o 4, pp.
2541-2551, oct. 2011, doi: 10.1109/TPWRD.2011.2165972.
[36] K. Nara, Y. Hayashi, K. Ikeda, y T. Ashizawa, «Application of tabu search to optimal placement
of distributed generators», 2001, vol. 2, pp. 918-923, doi: 10.1109/PESW.2001.916995.
[37] M. E. H. Golshan y S. Ali Arefifar, «Optimal allocation of distributed generation and reactive
sources considering tap positions of voltage regulators as control variables», Eur. Trans.
Electr. Power, vol. 17, n.o 3, pp. 219-239, may 2007, doi: 10.1002/etep.130.
[38] C. Novoa y T. Jin, «Reliability centered planning for distributed generation considering wind
power volatility», Electr. Power Syst. Res., vol. 81, n.o 8, pp. 1654-1661, ago. 2011, doi:
10.1016/j.epsr.2011.04.004.
[39] A. M. El-Zonkoly, «Optimal placement of multi-distributed generation units including different
load models using particle swarm optimisation», IET Gener. Transm. Distrib., vol. 5, n.o 7, p.
760, 2011, doi: 10.1049/iet-gtd.2010.0676.
[40] W. Prommee y W. Ongsakul, «Optimal multiple distributed generation placement in
microgrid system by improved reinitialized social structures particle swarm optimization»,
Eur. Trans. Electr. Power, vol. 21, n.o 1, pp. 489-504, ene. 2011, doi: 10.1002/etep.456.
[41] M. H. Moradi y M. Abedini, «A combination of genetic algorithm and particle swarm
optimization for optimal DG location and sizing in distribution systems», Int. J. Electr. Power
Energy Syst., vol. 34, n.o 1, pp. 66-74, ene. 2012, doi: 10.1016/j.ijepes.2011.08.023.
[42] L. Wang y C. Singh, «Reliability-Constrained Optimum Placement of Reclosers and Distributed
Generators in Distribution Networks Using an Ant Colony System Algorithm», IEEE Trans. Syst.
Man Cybern. Part C Appl. Rev., vol. 38, n.o 6, pp. 757-764, nov. 2008, doi:
10.1109/TSMCC.2008.2001573.
[43] F. S. Abu-Mouti y M. E. El-Hawary, «Optimal Distributed Generation Allocation and Sizing in
Distribution Systems via Artificial Bee Colony Algorithm», IEEE Trans. Power Deliv., vol. 26, n.o
4, pp. 2090-2101, oct. 2011, doi: 10.1109/TPWRD.2011.2158246.
[44] L. D. Arya, A. Koshti, y S. C. Choube, «Distributed generation planning using differential
evolution accounting voltage stability consideration», Int. J. Electr. Power Energy Syst., vol.
42, n.o 1, pp. 196-207, nov. 2012, doi: 10.1016/j.ijepes.2012.04.011.
[45] H. Hamedi y M. Gandomkar, «A straightforward approach to minimizing unsupplied energy
and power loss through DG placement and evaluating power quality in relation to load
variations over time», Int. J. Electr. Power Energy Syst., nov. 2011, doi:
10.1016/j.ijepes.2011.10.001.
[46] B. Banerjee y S. M. Islam, «Reliability based optimum location of distributed generation», Int.
J. Electr. Power Energy Syst., vol. 33, n.o 8, pp. 1470-1478, oct. 2011, doi:
10.1016/j.ijepes.2011.06.029.
[47] D. H. Popović, J. A. Greatbanks, M. Begović, y A. Pregelj, «Placement of distributed generators
and reclosers for distribution network security and reliability», Int. J. Electr. Power Energy
Syst., vol. 27, n.o 5-6, pp. 398-408, jun. 2005, doi: 10.1016/j.ijepes.2005.02.002.
[48] H. M. Khodr, M. R. Silva, Z. Vale, y C. Ramos, «A probabilistic methodology for distributed
generation location in isolated electrical service area», Electr. Power Syst. Res., vol. 80, n.o 4,
pp. 390-399, abr. 2010, doi: 10.1016/j.epsr.2009.10.001.
176
[49] F. Rotaru, G. Chicco, G. Grigoras, y G. Cartina, «Two-stage distributed generation optimal
sizing with clustering-based node selection», Int. J. Electr. Power Energy Syst., vol. 40, n.o 1,
pp. 120-129, sep. 2012, doi: 10.1016/j.ijepes.2012.02.012.
[50] G. N. Koutroumpezis y A. S. Safigianni, «Optimum allocation of the maximum possible
distributed generation penetration in a distribution network», Electr. Power Syst. Res., vol. 80,
n.o 12, pp. 1421-1427, dic. 2010, doi: 10.1016/j.epsr.2010.06.005.
[51] M. F. Shaaban, Y. M. Atwa, y E. F. El-Saadany, «DG allocation for benefit maximization in
distribution networks», IEEE Trans. Power Syst., vol. 28, n.o 2, pp. 639-649, may 2013, doi:
10.1109/TPWRS.2012.2213309.
[52] L. Rincón, Introducción a la probabilidad, 2da ed. México, DF: Departamento de Matemáticas
Facultad de Ciencias UNAM Circuito Exterior de CU 04510 México CDMX, 2014.
[53] A. Hadian, M.-R. Haghifam, J. Zohrevand, y E. Akhavan-Rezai, «Probabilistic approach for
renewable dg placement in distribution systems with uncertain and time varying loads», jul.
2009, pp. 1-8, doi: 10.1109/PES.2009.5275458.
[54] Y. G. Hegazy, M. M. A. Salama, y A. Y. Chikhani, «Adequacy assessment of distributed
generation systems using monte carlo simulation», IEEE Trans. Power Syst., vol. 18, n.o 1, pp.
48-52, feb. 2003, doi: 10.1109/TPWRS.2002.807044.
[55] R. C. Dugan, D. Montenegro, y A. Ballanti, «Reference Guide The Open Distribution System
Simulator (OpenDSS)», p. 199, 2018.
[56] «Resources | PES Test Feeder». http://sites.ieee.org/pes-testfeeders/resources/ (accedido
jul. 19, 2018).
[57] «CYME - Software, Análisis de sistemas de distribución».
http://www.cyme.com/es/software/cymdist/ (accedido jul. 08, 2019).
[58] R. H. Lasseter y P. Paigi, «Microgrid: A Conceptual Solution», IEEE 35th Annu. Power Electron.
Spec. Conf., vol. vol 6, n.o PESC 04, p. 6, 2004.
[59] M. Meiqin, L. Chang, y D. Ming, «Integration and Intelligent Control of Micro-Grids with Multi-
Energy Generations: A Review», Proc. IEEE Int. Conf. Sustain. Energy Technol., p. pp 777–780,
2008.
[60] S. Morozumi, «Micro-grid Demonstration Projects in Japan», IEEE Power Convers. Conf.-
Nagoya, n.o PCC 07, p. 8, 2007.
[61] A. Arulampalam, M. Barnes, A. Engler, A. Goodwin, y N. Jenkins, «Control of power electronic
interfaces in distributed generation microgrids», Int J Electron 919, pp. 503–523, 2004.
[62] D. Moskovitz, «Profits and Progress Through Distributed Resources», Regul. Assist. Proj., p.
35, 2000.
[63] J. A. P. Lopes, N. Hatziargyriou, J. Mutale, P. Djapic, y N. Jenkins, «Integrating distributed
generation into electric power systems: A review of drivers, challenges and opportunities»,
Electr. Power Syst. Res., p. 15, 2007.
[64] M. S. Mahmoud y F. M. Al-Sunni, Control and optimization of distributed generation systems.
Cham Heidelberg New York: Springer, 2015.
[65] M. Nigris, F. Bernardelli, y M. F. Coviello, «REDES INTELIGENTES DE ENERGÍA EN AMERICA
LATINA Y EL CARIBE: VIABILIDAD Y DESAFIOS», Chile, oct. 2010, p. 145.
[66] D. Bertini, R. Cicoria, J. S. de Assis, A. Silvestri, D. Falabretti, y V. Olivieri, «La penetrazione
della Generazione Distribuita nelle attuali reti elettriche italiane», Proy. Investig. Sobre El Sist.
Eléctrico Reti Attive Generazione Distrib. E Sist. Accumulo, vol. código 10000473, p. 121, 2009.
[67] T. Ackermann, «Distributed generation: a definition», Electr. Power Syst. Res., p. 10, 2001.
[68] LABEIN TECNALIA, «Guía Básica de la Generación Distribuida», Dir. Gen. Ind. Energ. Minas
Comunidad Madr., p. 65, 2005.
[69] INTERNATIONAL ENERGY AGENCY, «OVERVIEW OF RENEWABLES AND WASTE IN THE
WORLD», https://www.iea.org/statistics/renewables/, p. 12, 2019.
177
[70] I. Marti, «IEA WIND TCP ANNUAL REPORT», https://community.ieawind.org/publications/ar,
p. 164, 2017.
[71] INTERNATIONAL ENERGY AGENCY, «SNAPSHOT 2019 OF GLOBAL PHOTOVOLTAIC MARKET»,
Httpwwwiea-Pvpsorgindexphpid92, p. 10, 2019.
[72] L. Rincon, Introducción a la probabilidad, 3.a ed. Mexico D.F: Facultad de Ciencias UNAM,
2014.
[73] L. Rincon, Curso Intermedio de Probabilidad, 2.a ed. Mexico D.F: Facultad de Ciencias UNAM,
2007.
[74] B. E. V. Lizcano, «Evaluación de errores de datos sismológicos y sus efectos sobre los
parámetros y modelos derivados de estos para una región determinada cubierta por la red
del Observatorio Sismológico del SurOccidente», Universidad del Valle, Colombia, Santiago de
Cali, 2003.
[75] S. Ramírez Castaño, Redes de Distribución de Energía, 3.a ed. Manizales, Colombia:
Universidad Nacional de Colombia, 2006.
[76] R. C. Dugan y R. F. Arritt, «The IEEE 8500-Node Test Feeder», IEEE Power Energy Soc., p. 6,
2010.
[77] W. H. Kersting, Distribution System Modeling and Analysis, Third Edition, 3.a ed. New Mexico
State University: CRC Press, 2012.
[78] R. C. Dugan, D. Montenegro, y A. Ballanti, «Reference Guide The Open Distribution System
Simulator (OpenDSS)», p. 199, 2018.
[79] R. Baldick, E. Dept, y J. V. Staron, «IEEE 123 Node Test Feeder», IEEE Power Energy Soc., p. 40,
2014.
[80] «Requirements and Ports for MATLAB Parallel Server - MATLAB & Simulink - MathWorks
América Latina». https://la.mathworks.com/help/matlab-parallel-server/requirements-and-
ports-for-matlab-parallel-server.html (accedido abr. 11, 2020).
[81] «HOMER Pro - Microgrid Software for Designing Optimized Hybrid Microgrids».
https://www.homerenergy.com/products/pro/index.html (accedido jul. 08, 2019).
[82] S. Pfenninger y I. Staffell, «Long-term patterns of European PV output using 30 years of
validated hourly reanalysis and satellite data», Energy, vol. 114, pp. 1251-1265, nov. 2016,
doi: 10.1016/j.energy.2016.08.060.
[83] I. Staffell y S. Pfenninger, «Using bias-corrected reanalysis to simulate current and future
wind power output», Energy, vol. 114, pp. 1224-1239, nov. 2016, doi:
10.1016/j.energy.2016.08.068.
[84] S. Ross, Simulation, 5.a ed. Estados Unidos, California: Elsevier, 2013.
[85] M. Alonso, «Gestión Óptima de Potencia Reactiva en Sistemas Eléctricos con Generación
Eólica», Carlos III de Madrid, España, Leganés, 2010.
[86] M. A. Rosato, Small Wind Turbines for Electricity and Irrigation, 1.a ed., vol. Taylor&Francis
Group, 1 vols. New York, Estados Unidos: CRC Press, 2019.
[87] B. Zhao, C. Wang, y X. Zhang, Grid-Integrated and Standalone Photovoltaic Distributed
Generation Systems: Analysis, Design, and Control. Chichester, UK: John Wiley & Sons
Singapore Pte. Ltd, 2017.
178
ANEXO A. Ejemplos de aplicación de las funciones de probabilidad
Rayleigh y Weibull para cálculo de disponibilidad de energía eólica
Con una densidad del aire de 𝜌 = 1.22 (𝑘𝑔/𝑚3 ) y un área de cobertura de la turbina 𝑆 =
1 (𝑚2 ) se aplica la ec. (1.9) para intervalos de velocidad desde 0 a 14 (m/s) cuyos
resultados se los presenta en la Tabla A.1 y Figura A.1
Tabla A.1. Función de probabilidad de Rayleigh para una velocidad promedio de 6 (m/s) [86].
179
Figura A.1. Función de probabilidad de Rayleigh para una velocidad promedio de 6 (m/s) [82].
180
Tabla A.2. Función de probabilidad de Weibull para una velocidad promedio de 6 (m/s) [86].
Figura A.2. Función de probabilidad de Weibull para una velocidad promedio de 6 (m/s) [82].
181
ANEXO B. Análisis eléctricos factibles con OPENDSS
182
ANEXO C. Comparación de resultados de CYME con IEEE para la
red de prueba de 123 nodos.
183
184
185
186
187
188
189
190
191
192
ANEXO D. Estructura del programa base para determinar la mejor
ubicación y tamaño de Capacitores – Reactores en GD con la
metodología probabilística bivariada.
a) Se generan vectores aleatorios para posición y tamaño de cada capacitor en las fases
A, B y C para una distribución bivariada uniforme, Weibull y sus combinaciones, cuya
estructura es similar al literal d) de la sección 2.4.4 como se indica a continuación:
193
c) Ya que no es factible realizar conexiones de capacitores en todos los nodos de la red,
es necesario seleccionar los casos idóneos, la programación es personalizada para
cada alimentador con la sección “quinta programación, asig5_prog” y es generada de
la misma forma de los literales k) y l) de la sección 2.4.1, como se indica a continuación:
d) Con los casos idóneos de capacitores seleccionados, se debe hacer una simulación de
flujo de potencia para cada vector del capacitor de la fase A y el resultado del mejor
GD. Una estructura similar a la del literal e) de la sección 2.4.4 pero en esta interacción
se aumenta los registros de voltaje promedio y varianza como se indica a continuación:
194
195
e) Se selecciona la mejor opción de voltaje promedio más próxima al 1 p.u. para el
capacitor A sin que aumente las pérdidas obtenidas en el literal g) de la sección 2.4.4
con el fragmento de programación que se indica a continuación:
f) Se repite los literales d) y e) para los capacitores B y C. Con los resultados finales se
aplica la misma estructura de programación indicada en el literal i) y se aumenta los
resultados de los vectores con las mejores posiciones y potencias para capacitores. El
fragmento de programación se indica a continuación:
196
197
ANEXO E. Configuración para conexiones de Generación
Distribuida.
Figura E.1. Diseños de turbinas con 3 alabes e inserción en eje central [86].
198
almacenamiento de energía en bancos de baterías u otras aplicaciones (cuarta parte de la
Figura E.2).
199
200
MSc. Sebastián Ribadeneira Ramírez
ESCUELA POLITÉCNICA NACIONAL
201