Ups - TTS208 Analisis de Cargas de Redes Electricas

“ANÁLISIS DE CARGABILIDAD DE REDES ELÉCTRICAS DE DISTRIBUCIÓN BASADO EN LA
MODIFICACIÓN DEL ALGORITMO DE GAUSS-SEIDEL”

UNIVERSIDAD POLITÉCNICA SALESIANA
SEDE QUITO
CARRERA:
INGENIERÍA ELÉCTRICA
Trabajo de titulación previo a la obtención del título de

INGENIERO ELÉCTRICO
TEMA:
“ANÁLISIS DE CARGABILIDAD DE REDES ELÉCTRICAS DE
DISTRIBUCIÓN BASADO EN LA MODIFICACIÓN DEL ALGORITMO DE
GAUSS-SEIDEL”
AUTOR:
HENRY SAÚL GANCINO BUSTILLOS
TUTOR:
ING. DIEGO FRANCISCO CARRIÓN GALARZA
Quito, noviembre 2020

Henry Saúl Gancino Bustillos
“ANÁLISIS DE CARGABILIDAD DE REDES ELÉCTRICAS DE DISTRIBUCIÓN

BASADO EN LA MODIFICACIÓN DEL ALGORITMO DE GAUSS-SEIDEL”
.
Universidad Politécnica Salesiana

Ingeniería Eléctrica
Breve reseña historia e información de contacto:

Henry Saúl Gancino Bustillos (Y'1996-J'07). Realizó sus
estudios secundarios en el colegio técnico Sucre, se graduó de
Bachiller Técnico en Electricidad. Egresado de la Carrera de
Ingeniería Eléctrica de la Universidad Politécnica Salesiana. Su
trabajo de investigación se basa en modificar el algoritmo de
Gauss-Seidel para que pueda ser utilizado en sistemas de
distribución y con ayuda de los resultados del flujo de potencia
encontrar el punto crítico en la línea en donde sobrepasado ese
punto va a existir una expansión térmica en los conductores.
hgancino@est.ups.edu.ec
Dirigido por:
Diego Francisco Carrión Galarza (Y'1981-M'12). Se graduó
como Ingeniero Eléctrico en la Universidad Politécnica Salesiana,
Ecuador en el año 2010 y en la actualidad está trabajando para
lograr su título de doctor en ingeniería de la Universidad
Politécnica Bolivariana – Medellín, Colombia. Es profesor e
investigador de la Universidad Politécnica Salesiana – Quito,
Ecuador. En la actualidad es miembro del grupo de investigación
GIREI (Grupo de Investigación en Redes Eléctricas Inteligentes
– Smart Grid Research Group).
dcarrion@ups.edu.ec
Todos los derechos reservados:

Queda prohibida, salvo excepción prevista en la ley, cualquier forma de reproducción,
distribución, comunicación pública y transformación de esta obra para fines comerciales,
sin contar con la autorización de los titulares de propiedad intelectual. La infracción de
los derechos mencionados puede ser constitutiva de delito contra la propiedad intelectual.
Se permite la libre difusión de este texto con fines académicos o investigativos por
cualquier medio, con la debida notificación a los autores.
DERECHOS RESERVADOS ©2020
Universidad Politécnica Salesiana
QUITO-ECUADOR
I
DECLARATORIA DE COAUTORÍA DEL DOCENTE TUTOR/A
Yo, Diego Francisco Carrión Galarza declaro que bajo mi dirección y asesoría fue
desarrollado el trabajo de titulación denominado “ANÁLISIS DE CARGABILIDAD DE
REDES ELECTRICAS DE DISTRIBUCION BASADO EN LA MODIFICACIÓN DEL
ALGORITMO DE GAUSS-SEIDEL” realizado por Henry Saúl Gancino Bustillos,
obteniendo un producto que cumple con todos los requisitos estipulados por la
Universidad Politécnica Salesiana para ser considerados como trabajo final de titulación.
Quito, noviembre 2020
………………………………….
Diego Francisco Carrión Galarza
Cédula de identidad: 1713703062
II
CESIÓN DE DERECHOS DE AUTOR
Yo, Henry Saúl Gancino Bustillos, con documento de identificación N° 1725114076,

manifiesto mi voluntad y cedo a la Universidad Politécnica Salesiana la titularidad sobre
los derechos patrimoniales en virtud de que soy autor del trabajo de grado/titulación
intitulado: “ANALISIS DE CARGABIIDAD DE REDES ELÉCTRICAS DE
DISTRIBUCIÓN BASADO EN LA MODIFICACIÓN DEL ALGORITMO DE GAUSS-
SEIDEL”, mismo que ha sido desarrollado para optar por el título de: Ingeniero Eléctrico,
en la Universidad Politécnica Salesiana, quedando la Universidad facultada para ejercer
plenamente los derechos cedidos anteriormente.
En aplicación a lo determinado en la Ley de Propiedad Intelectual, en mi condición de
autor me reservo los derechos morales de la obra antes citada. En concordancia, suscribo
este documento en el momento que hago entrega del trabajo final en formato digital a la
Biblioteca de la Universidad Politécnica Salesiana.
.............................................
Nombre: Henry Saúl Gancino Bustillos
Cédula: 1725114076
Fecha: noviembre 2020
III
DEDICATORIA
El presente trabajo de investigación se lo dedico a Dios quien nunca dejo de bendecirme

durante este largo proceso de estudio, con mucho cariño a mis padres, quienes han sido
un pilar fundamental en el transcurso de mi vida brindándome su apoyo en todo momento.
A mi hermano y a mi cuñada quienes vivieron de cerca todo el proceso para obtener este
magnífico logro.
IV
AGRADECIMIENTOS
Agradezco a Dios por bridarme su guía y sus bendiciones durante toda mi vida, a mis
padres que siempre creyeron en mí y me dieron todo el apoyo necesario para culminar mis
estudios universitarios.
Agradezco a la Universidad Politécnica Salesiana, que me abrió las puertas para forjarme
como profesional, a mis amigos con los que viví todo el proceso universitario.
Agradezco a la Empresa Eléctrica Quito, por darme la oportunidad de realizar mis

practicas preprofesionales y despertar mi interés por el área de distribución eléctrica.
Agradezco al Ingeniero Diego Carrión, quien es el tutor del presente artículo, por la
apertura para dirigir el trabajo desarrollado.
V
ÍNDICE GENERAL
1. Introducción ........................................................................................................... 2
2. Gauss-Seidel Para Flujos de Potencia. ................................................................... 4
2.1. Flujos de Potencia en Redes Radiales ................................................................ 5
2.2. Parámetros de Cargabilidad. .............................................................................. 6
2.3. Expansión térmica en los conductores ............................................................... 8
2.4. Comportamiento Térmico en los conductores. .................................................. 9
2.5. Pandeo en las Líneas ........................................................................................ 11
3. Estimación de Conductores en Distribución ........................................................ 12
3.1. Selección de los Conductores Eléctricos.......................................................... 12
3.2. Determinación de la Corriente por un Flujo de Carga. .................................... 13
3.3. Efecto del Exceso de Corriente y Temperatura en los Conductores ................ 14
4. Metodología Para Determinar Flujos en los Sistemas de 13 y 34 nodos. ............ 14
4.1. Algoritmo Para la Aplicación de Gauss-Seidel en Distribución. ..................... 15
4.2. Verificación de los Resultados por DIgSILENT Power Factory. .................... 16
4.3. Calibre de Conductores Seleccionados de Acuerdo al Nivel de Corriente. ..... 17
4.4. Límite de Cargabilidad de los Alimentadores de 13 y 34 Nodos. ................... 18
4.5. Expansión Térmica del Sistema Alimentador de 13 y 34 Nodos de la IEEE. . 19
5. Resultados ............................................................................................................ 20
5.1. Flujo de Potencia. ............................................................................................. 20
5.2. Flujos en DIgSILENT de los Sistemas de 13 y 34 Nodos. .............................. 20
5.3. Punto Crítico Antes de que Exista una Extensión Térmica. ............................ 21
6. Conclusiones ........................................................................................................ 22
7. Trabajos futuros ................................................................................................... 23
VI
8. Referencias ........................................................................................................... 23
8.1. Matriz de Estado de Arte.................................................................................. 27
8.2. Resumen de Indicadores .................................................................................. 30
VII
INDICE DE FIGURAS
Figura 1. Etapas de planificación ........................................................................................................ 4
Figura 2. Temática de la investigación.............................................................................................. 4
Figura 3. Modelo de impedancias en serie y en derivación. ....................................................... 4
Figura 4. Topología radial. ................................................................................................................... 6
Figura 5. Perdidas de una red subterránea de distribución. ......................................................... 7
Figura 6. Perdidas de una red aérea de distribución. ..................................................................... 7
Figura 7. Límite del nivel de sobrecarga SIL. ................................................................................. 7
Figura 8. Curva características de las líneas. ................................................................................... 8
Figura 9. Estado de la curva a diferentes temperaturas ................................................................ 9
Figura 10. Representación gráfica de las variables del balance térmico ............................... 11
Figura 11. Sistema de alimentación de 13 nodos de la IEEE ................................................... 15
Figura 12. Sistema de alimentación de 34 nodos de la IEEE. .................................................. 15
Figura 13. Simulación del flujo del alimentador de 13 nodos de la IEEE. ........................... 16
Figura 14. Simulación del flujo del alimentador de 34 nodos de la IEEE. ........................... 17
Figura 15. Corriente proporcional a temperatura en el sistema de 13 nodos de la IEEE. 18
Figura 16. Corriente proporcional a temperatura en el sistema de 34 nodos de la IEEE. 19
Figura 17. Punto crítico promedio en el sistema de 13 nodos de la IEEE. ........................... 19
Figura 18. Punto crítico promedio en el sistema de 34 nodos de la IEEE. ........................... 20
Figura 19. Resultados del flujo de potencia del alimentador de 13 nodos de la IEEE...... 20
Figura 20. Resultados del flujo de potencia del alimentador de 13 nodos de la IEEE...... 20
Figura 21. Flujo en DIgSILENT del alimentador de 13 nodos. ............................................. 20
Figura 22. Flujo en DIgSILENT del alimentador de 34 nodos. ............................................... 21
Figura 23. Punto crítico de la configuración 601. ........................................................................ 21

VIII
Figura 31. Resumen e indicador de la temática - Estado del arte............................................ 30
Figura 32. Indicador de formulación del problema - Estado del arte. .................................... 30
Figura 33. Indicador de solución - Estado del arte ...................................................................... 31
IX
INDICE DE TABLAS
Tabla 1. Variables del balance térmico. ........................................................................... 10
Tabla 2: Variables del balance térmico ............................................................................ 11
Tabla 3: Datos del alimentador de 34 nodos de la IEEE.................................................. 13
Tabla 4: Datos del alimentador de 13 nodos de la IEEE.................................................. 13
Tabla 5: Tipo de conductor eléctrico del sistema de 13 nodos de la IEEE. ..................... 17
Tabla 6: Tipo de conductor eléctrico del sistema de 34 nodos de la IEEE. ..................... 17
Tabla 7: Configuración de las líneas del sistema de 13 nodos de la IEEE. ..................... 18
Tabla 8: Configuración de las líneas del sistema de 13 nodos de la IEEE. ..................... 19
Tabla 9:Matriz del estado del arte. ................................................................................... 27
X
ANÁLISIS DE CARGABILIDAD DE REDES ELÉCTRICAS DE
DISTRIBUCIÓN BASADO EN LA MODIFICACIÓN DEL
ALGORITMO DE GAUSS-SEIDEL.
Resumen Abstract
En los sistemas alimentadores, una de las In the feeder systems, the main function is
funciones principales es la distribución de the distribution of energy, that is distribute
energía, misma que se distribuye a través through the electric conductors with a
de conductores eléctricos con una radial topology. The bad operation of this
topología radial. El mal funcionamiento de systems is produced by several reason, the
estos sistemas se produce por diferentes principal reason occur when the current
razones, una de las principales ocurre surpasses the overload factor in the line,
cuando la corriente supera el factor de resulting in the conductor surpasses its
sobrecarga en la línea, creando que el thermal limit, but this can be avoided by
conductor supere su límite térmico, sin analysing the potency flow of the
embargo, esto se puede evitar analizando distribution systems, so that exist reliable
el flujo de potencia de los sistemas de results in radial systems and of this form
distribución, para que existan resultados establish the critical current point in the
confiables en sistemas radiales y de esta conductor before thermal expansion
manera determinar el punto de corriente occur.
crítico en el conductor antes de que ocurra
una extensión térmica. For the mentioned, the purpose of this
article is to adapt of the algorithm of
Por lo mencionado, el siguiente articulo Gauss-Seidel in distribution systems, and
tiene como finalidad la adaptación del to decide a current critical point in the
algoritmo de Gauss-Seidel en sistemas de conductors of the 13 and 34 nodes feeder
distribución, y determinar un punto crítico systems of the IEEE.
de corriente en los conductores de los Through the results obtained, it will be
sistemas alimentadores de 13 y 34 nodos may to define at which current level a
del IEEE. thermal extension begins to exist in the
A través de los resultados obtenidos se conductors and the results of the potency
podrá definir a que nivel de corriente flow will be verified with a simulator
empieza a existir una extensión térmica en which will affirm the results obtained.
los conductores y se verificarán los
resultados del flujo de potencia con un
simulador el cual afirmará los resultados
obtenidos.
Palabras Clave: Máxima transferencia de Keywords: Maximum energy transfer,

energía, limite térmico, cargabilidad, flujo thermal limit, chargeability, potency flow,
de potencia, dilatación en el conductor. dilatation in the conductor.
1
1. Introducción
Las redes eléctricas de distribución, son método. La topología de las redes de

las que se encargan de transportar energía distribución es radial por lo que se aplica
eléctrica de medio voltaje para su un flujo de potencia de nodo a nodo
comercialización, ya sea su entrega al tomando en consideración las
cliente en medio o bajo voltaje. Dichas impedancias de las líneas, los voltajes y
redes parten desde una subestación las cargas que circulan por estas. El
reductora y su topología es radial y/o en algoritmo de Gauss-Seidel se basa en una
anillo abierto por si existe alguna serie de iteraciones hasta llegar a un error
contingencia. Estas deben ser capaces de muy mínimo, prácticamente este es el
soportar la carga que circula sin presentar problema de usar Gauss-Seidel en
ningún problema, cabe indicar que las distribución ya que los datos obtenidos no
redes de distribución pueden estar sobre son muy confiables pero algunos autores
la superficie (sistema aéreo) o debajo de realizan una cierta modificación al
la superficie (sistema soterrado) y es método de Gauss-Seidel para que sea
fundamental determinar el límite de aplicable en sistemas eléctricos de
cargabilidad que pueden soportar los distribución [2]–[4].
conductores eléctricos.
En los artículos [5], [6]consideran un
Se tomo ejemplos de sistemas de nodo principal y le dan el nombre de
distribución de la IEEE para estudiar los Slack mientras que todos los demás nodos
efectos causados cuando la corriente que son considerados nodos PQ, si bien se
circula por el conductor sobrepasa el sabe que los flujos en transmisión son en
nivel para el que fue diseñado, conocidos sistemas eléctricos con una topología en
también como sistemas alimentadores. Es malla y en distribución su topología es
necesario identificar que este artículo radial al momento de tomar o ir
trata únicamente del análisis del sistema considerando solo dos nodos a la vez para
eléctrico a nivel de distribución, tomando la realización del flujo de potencia los
en consideración la parte de resultados van a ser válidos.
comercialización eléctrica ya que se la
referencia como puntos de carga. El En [1],[7]–[9] muestran que Gauss-Seidel
algoritmo de Gauss-Seidel para flujos de si es aplicable en distribución o en una red
potencia generalmente se los utiliza en radial, por lo cual el objetivo de utilizar
transmisión debido a que toma en cuenta este algoritmo es encontrar el punto de
la parte capacitiva que se producen en las demanda critico que pueden soportar las
líneas debido al alto voltaje, a la distancia líneas de distribución antes de que exista
de la línea y a la altura en la que se una extensión o dilatación térmica en los
encuentra, pero en este artículo conductores, con el fin de visualizar el
académico se aplicará Gauss-Seidel en calibre ideal de conductor para que no
distribución, tomando en cuenta el exista tal efecto, el cual causa que el
artículo[1], el autor realiza flujos de conductor sufra una dilatación
potencia en distribución por el método de volumétrica y con oscilaciones de viento,
Gauss-Seidel, pero debido a las distancias al no existir la debida distancia entre las
y a las magnitudes de voltaje que son fases, estas rompan la distancia de
inferiores a las magnitudes de voltaje que seguridad causando cortocircuitos o
existe en transmisión no se obtienen arranques del sistema, por lo cual al
resultados confiables por lo que se realizar el análisis se va a prever este tipo
realizan ciertas adecuaciones en el de contingencias o accidentes en las
algoritmo o en la forma de aplicar el líneas. Además, una vez que ya exista una
2
dilatación volumétrica en los este algoritmo el problema es que tiene un
conductores, se analizará comparando los exceso de memoria por lo que se
resultados obtenidos por el flujo con el almacenan datos de cada iteración. Por lo
(Surge Impedance Loading o Impedancia que se toma en consideración que el
natural de la línea) SIL nominal de la problema de los resultados es por la
línea, si estos necesitan ser cambiados o topología y por los valores de impedancia
solo bastará con una adecuación o ajuste que existen, después de analizar[11]–
en el conductor. Este efecto se visualiza [13], en estos documentos se puede
en la parte de transmisión como en la de determinar que se usan otros algoritmos o
distribución por lo cual se va a tomar métodos para realizar un análisis en
datos o estudios realizados en transmisión distribución. La forma en la que se
y referirlos a la parte de distribución. modifica el flujo o la forma en la que se
le aplica se detallará en los siguientes
Al momento que se habla de dilatación puntos del artículo, además cuando se
térmica también se habla del pandeo que habla de un análisis de flujo en
existe en las líneas debido al aumento de distribución por lo general realizan el
temperatura en los conductores causados análisis de nodo a nodo tomando en
por niveles de sobre carga o porque las cuenta todo el tramo, sus características y
condiciones de los conductores eléctricos se realiza un análisis de carga trifásico.
no son los adecuados para soportar esa
cantidad de carga, por lo cual, con el Al analizar los resultados del flujo de
algoritmo de Gauss-Seidel se obtendrán potencia por Gauss-Seidel se los
resultados que ayudaran a determinar compararán con los resultados obtenidos
hasta qué punto de carga pueden soportar en el software DIgSILENT POWER
las redes eléctricas de distribución antes FACTORY con el método de (Newton-
de que exista una expansión térmica. Se Raphson, clásicas ecuaciones de
verificarán los datos obtenidos por un potencia) el cual se encargará de
software en el cual se pueden determinar comparar los resultados obtenidos por el
flujos de potencia en líneas de flujo realizado en Matlab. Con los
distribución. Una herramienta resultados confirmados se realizará la
fundamental para la comprensión de las prevención de que exista una extensión
redes eléctricas de distribución son los térmica en los conductores, por lo cual se
flujos de potencia monofásicos en los va a encontrar un punto de demanda
cuales toman en cuenta un sistema crítico, con lo que se va a empezar a
balanceado para las tres fases, estos son determinar la extensión térmica que se va
diseñados específicamente para el a producir en el conductor eléctrico o el
análisis de carga en transmisión pero se nivel de corriente que va a pasar antes de
los utilizan también en distribución que exista la dilatación en el cable, lo cual
debido a la facilidad de cálculo por lo que va a evitar que en los cables de
se asume que son sistemas balanceados, distribución exista un pandeo en las líneas
para que el almacenamiento de memoria y esto pueda ocasionar algún conflicto
del flujo no sea tan extendido y exista una con las otras líneas, va a prevenir el
convergencia más rápida en los contacto con las otras líneas u objetos
resultados [5], [6], [10]. cercanos, esto también ayudara a que la
La aplicación en distribución de los flujos vida útil del conductor sea más larga.
de potencia monofásicos no da buenos
resultados de convergencia, la decisión de
analizar las líneas por Gauss-Seidel es por
la creación de la matriz de impedancia lo
cual dan muestras más rápido, pero en
3
Flujos de Potencia 2. Gauss-Seidel Para Flujos de
en Distribución
Potencia.
Modificación del método de Gauss -Seidel
para aplicarlo en distribución (resultados
mas confiables) El método de Gauss-Seidel se lo utiliza
Comparación de los resultados con los de
DIgSILENT (Verificar convergencia de los
para determinar los voltajes en las barras
resultados) o en los nodos y sus ángulos, además de
Estimación de la corriente y el calibre de encontrar estos dos parámetros también
conductor de cada tramo de los sistemas
de 13 y 34 nodos se puede encontrar a través de una serie
Determinación del punto crítico en el de iteraciones la potencia reactiva y
que empieza a darse la dilatación en los
conductores eléctricos activa que circula entre dos nodos, el
Evaluación de la relación del nivel de número de iteraciones corresponden al
temperatura-corriente desde el punto
que se empieza a dar la extensión
nivel de tolerancia que se le dará
térmica en los conductores dependiendo de la utilización de las
Figura 1. Etapas de planificación líneas. Algo que ayuda a entender los
flujos de potencia es la relación que existe
En este artículo se propone un modelo de
entre 𝑃, 𝑄, |𝑉| y 𝛿, en muchos casos se
flujo de potencia modificado para
desprecia la resistencia serie con la
analizar modelos radiales de la IEEE, con
finalidad de facilitar la resolución del
los resultados obtenidos se va a
flujo, tomando en consideración que en la
determinar el límite de cargabilidad de las
mayoría de los sistemas su resistencia es
líneas de distribución en los cuales se
muy baja por lo cual no afecta los
determinará a que punto de corriente
resultados.
empieza a existir una extensión térmica
en los conductores por lo que se podrá
llegar a la conclusión hasta que nivel de
corriente puede circular por los
conductores sin que el aumento de
temperatura en los mismos produzca una
dilatación, ya que si existe la dilatación
puede causar problemas con el sistema
debido al pandeo que se produce por el
aumento de longitud en la línea.
Figura 3. Modelo de impedancias en serie
y en derivación.
𝑆1 = 𝑉1 ∗ 𝐼1∗ (1)
A continuación, se detallará el modelo

general para flujos de potencia en
sistemas eléctricos de dos nodos. Las
potencias que son inyectadas al nodo 1 se
la conoce como 𝑆1 y se la determina por
(1), la intensidad de corriente está
compuesta en los casos de las líneas de
transmisión por la sumatoria de la
corriente que pasa por 𝑌𝑠ℎ y 𝑍𝑠𝑒𝑟𝑖𝑒 , se las
puede identificar en la figura 3 y para los
sistemas de distribución la corriente es la
que circula solo por 𝑍𝑠𝑒𝑟𝑖𝑒 , Gauss-Seidel
Figura 2. Temática de la investigación se lo utiliza para N barras o nodos, donde
4
se calcula el voltaje en una barra K y esta
𝑛
es expresada en función de la potencia (8)
𝑄𝑖 = ∑|𝑉𝑖 𝑌𝑖𝑘 𝑉𝑘 | sin(𝛿𝑘 − 𝛿𝑖 + 𝛾𝑖𝑘 )
activa y reactiva, suponiendo cargas 𝑘=1
conectadas a las barras o estimando los
voltajes en las otras barras, dando un En el algoritmo de Gauss-Seidel se
inicio plano a las que no tienen datos el encuentra el método iterativo y los pasos
cual es el encargado de dar valores de para empezar a realizar este método,
volteje de 1 + 0𝑗 𝑝𝑢 a los nodos en también se da el tipo de barras o de nodo
donde no están especificados sus valores que se denomina antes de realizar el
de voltaje[14], [15]. En (2) determina el algoritmo lo cual se va a detallar a
modelo iterativo del voltaje. continuación.
𝑛 (2) En el método de Gauss-Seidel se define

1 𝑆𝑘∗ como primer paso una barra Slack, la cual
𝑉𝑘 = − ∑ 𝑌𝑘𝑖 𝑉𝑖
𝑌𝑘𝑘 𝑉𝑘∗ es la barra de referencia de todo el sistema
𝑖=1
( 𝑖≠𝑗 ) desde donde va a empezar el proceso
Donde: iterativo, como segundo paso se da un
𝑌𝑘𝑖 es la impedancia en p.u. inicio plano lo que significa dar un valor
𝑆𝑘∗ es la potencia aparente en p.u.
de 1.0 + 𝑗0 a los voltajes desconocidos,
𝑉𝑖 es el voltaje en p.u.
el tercer paso es definir las barras PQ a las
𝐼1∗ es la corriente en p.u.
que se les conoce también como barras de
𝑌11 𝑌12 cargas y se las resuelve para encontrar las
(3)
𝑌𝐵𝑈𝑆 = [ ] componentes reales e imaginarias de los
𝑌21 𝑌22
voltajes y como último paso se deben
𝑍11 𝑍12 (4) determinar las barras PV la cual se
𝑍𝐵𝑈𝑆 ≈ (𝑌𝐵𝑈𝑆 )−1 = [ ]
𝑍21 𝑍22 resuelve con (2) para determinar el
voltaje 𝑉𝑘 .
𝐼𝐵𝑈𝑆 = 𝑌𝐵𝑈𝑆 ∗ 𝑉𝐵𝑈𝑆 (5)
2.1. Flujos de Potencia en Redes
𝑉𝐵𝑈𝑆 = 𝑍𝐵𝑈𝑆 ∗ 𝐼𝐵𝑈𝑆 (6) Radiales
En (3) se observa la matriz de admitancias Las redes eléctricas de distribución al

nodales diseñada para un flujo de igual que las redes de transmisión tienen
potencia de dos barras, con lo cual se una capacidad máxima de transferencia
puede determinar la corriente que circula, de energía desde un punto a otro, este
a través de (5). Con la cual al despejar el límite de transferencia toma en cuenta las
volteje, esta se invierte y quedaría de la variaciones de voltaje y estabilidad en las
siguiente forma (6). líneas de una longitud larga, en cambio en
las líneas que tienen una longitud corta
En (3) y (4) se observa que las ecuaciones solo se toma en cuenta la corriente
son lineales y esta está acorde al modelo máxima que puede soportar el conductor
de flujo que se va implementar ya que eléctrico y la potencia que puede fluir por
esta modelación de la red eléctrica es el conductor. Los flujos de potencia en
lineal pero los datos que se conocen son distribución van a presentar condiciones
los de las potencias por lo cual se diferentes a los flujos en transmisión, se
determina por (7) y (8). encuentran características en las redes de
distribución principalmente su topología,
diferentes conexiones, toda clase de
𝑛
(7) cargas, líneas resistencias iguales a sus
𝑃𝑖 = ∑|𝑉𝑖 𝑌𝑖𝑘 𝑉𝑘 | cos(𝛿𝑘 − 𝛿𝑖 + 𝛾𝑖𝑘 )
𝑘=1
5
reactancias, entre otras características
[16]. En donde la potencia activa va a ir de un
Existen planteamientos de flujos de nodo hacia otro tomando en cuenta para
potencia para un tipo de topología en la dirección del flujo el ángulo de mayor
especial (radial o árbol) al cual no ser voltaje y el de menor voltaje, en cambio
mallada realiza adaptaciones en los la potencia reactiva toma en cuenta las
métodos tradicionales de flujos de magnitudes de voltaje para su dirección,
potencia, para determinar resultados que además en la potencia activa su dirección
converjan a una solución correcta. El siempre tiene un solo sentido mientras
articulo [9] describe un método iterativo que en la potencia reactiva su dirección va
por tramos de la red para solucionar los a depender de la inyección de reactivos a
problemas de flujos tradicionales en la línea o al tramo. Como se ha
topologías que no sea la de malla, mencionado antes que la topología en
simplifica el procedimiento del flujo distribución es distinta a la de
realizando particiones en el sistema. Para transmisión, los flujos en los sistemas
encontrar una adecuada solución para eléctricos de distribución parte desde un
todo el sistema, se realiza un número de nodo principal, a este nodo principal se le
particiones, después de determina un conoce como la fuente de inyección de
orden para cada partición energía por lo que se entiende que es una
subestación reductora.
650
2.2. Parámetros de Cargabilidad.

645 633 634
El nivel y la distribución de la
El análisis de carga en los sistemas
684 675
eléctricos de distribución sirve o son
671 692
utilizados para encontrar el voltaje que
existe en cada nodo de la red, la potencia
activa y reactiva que circulan por todas
las líneas, pérdidas en los
652 680
transformadores y una gran variedad de
Figura 4. Topología radial. demanda (cargas resistivas, inductivas y
Según [8], [9] describen límites de carga capacitivas).
en los cuales se usa un sistema general de
dos nodos o dos barras de los cuales parte También existen ventajas en los sistemas
así la investigación de aplicar Gauss- eléctricos de distribución lo cual una de
Seidel a un sistema de distribución, ellas es que los transformadores son los
aplicando separaciones y realizando alimentadores de cada carga, se considera
flujos en dos nodos. El problema con la una ventaja porque las cargas estarán
utilización del método tradicional de alimentas en un misma barra sin importa
Gauss-Seidel es la excesiva memoria que si es la carga es monofásica, bifásica o
usa para encontrar la convergencia de sus trifásica, la cargas en distribución como
valores por lo cual la solución a tal ya se mencionó anteriormente no todas
problema, será determinar un flujo por las cargas son trifásicas por lo que en
cada tramo del sistema y con los estos sistemas se utiliza una distribución
resultados obtenidos realizar el siguiente de carga para que el sistema sea lo más
flujo y continuar asi hasta llegar al último balanceado posible, el balanceo de carga
tramo del sistema, esto se lo ayuda a que la corriente que circula por
realizará para evitar el exceso de memoria cada conductor sea la adecuada y circule
en el método. la misma corriente por cada conductor
eléctrico. Para determinar la cargabilidad
6
que pueden soportar las líneas se Cuando existen valores menores a los de
determina el SIL [1], [16],[17]. La la demanda del punto de intersección que
capacidad de potencia activa como se puede ver en la figura 7, la línea
reactiva que circulan en una línea deben registra valores capacitivos en cambio
estar entre los parámetros de operación cuando los valores son más altos su
del sistema. Estos parámetros están en comportamiento es inductivo. Los
función de la distancia de la línea, el valores de potencia que circulan por las
voltaje de cada nodo y el calibre del líneas deben cumplir con las condiciones
conductor. Para encontrar el flujo de del SIL.
carga ideal se usará el flujo de potencia.
Las impedancias que existen en las redes
eléctricas de distribución son resistivas e
inductivas, en estos sistema la parte
capacitiva de la impedancia es muy baja
por lo cual no se la toma en
consideración, en estas solo existe la
parte capacitiva en las redes que son
soterradas debido al efecto que se da entre
la fase y la tierra del conductor, pero en
las redes eléctricas de distribución aéreas
no se da el efecto capacitivo debido a la
altura que se encuentran, a su longitud y Figura 7. Límite del nivel de sobrecarga
al nivel de voltaje[17]. SIL.
Para determinar el valor de potencia ideal
que pueden circular por las líneas (SIL),
la que se determina por (9), se necesita
conocer algunos parámetros como la
impedancia en serie, en derivación y la
impedancia característica de la línea. El
SIL se va a utilizar en los sistemas de
distribución debido a que se necesita
conocer el límite térmico de los
conductores de los sistemas del IEEE
Figura 5. Perdidas de una red para evaluar a que punto va a existir una
subterránea de distribución. dilatación en los conductores. Cuando el
valor de voltaje en las líneas es alto su
potencia de igual manera cumpliendo
proporcionalmente al nivel de voltaje, el
SIL sirve para determinar el número y
calibre de conductores adecuados para la
potencia que van a circular por estos, al
aumentar el número de conductores su
Pérdidas Resistivas Pérdidas Inductivas
impedancia característica va a disminuir
y se la determina a través de (10), la
impedancia característica de un tramo de
la línea es aquella que es casi siempre
Figura 6. Perdidas de una red aérea de resistiva con un factor de potencia muy
distribución. cerca de uno.
7
𝑧𝑐 = √(𝑧𝑠 )(𝑧𝑑 ) (9)
∗
(𝑉𝐿𝐿 )2 (10)
𝑆𝐼𝐿 = [ ]
𝑧𝑐
Donde:
𝑧𝑐 Impedancia Característica [Ω]
𝑧𝑠 Impedancia Serie [Ω/l]
𝑧𝑑 Impedancia en derivación [Ωl]
2.3. Expansión térmica en los

Figura 8. Curva características de las
conductores
líneas.
Al momento que fluye corriente por un Los conductores que se utilizan en las
conductor eléctrico de un punto a otro redes aéreas de distribución son cables
punto se va a generar un cambio de que no contienen un aislamiento o son
temperatura (Aumento de temperatura) popularmente conocidos como
en el conductor eléctrico mejor conocido conductores desnudos, se va a tomar en
como efecto Joule, para que exista una consideración la temperatura del
diferencia de temperatura además de la ambiente, el calor que se produce cuando
corriente también tiene mucha influencia circula la corriente eléctrica y la
los factores climáticos. La temperatura radiación, todas estas variables se le
del ambiente y la potencia eléctrica que puede aplicar en la ecuación de balance
fluye por los conductores no son las de calor, se la puede referir a distribución
mismas a toda hora del día, estas son ya que generalmente se la usa en
variables estocásticas por lo que en el transmisión. Los cálculos fueron
conductor eléctrico la temperatura realizados con la metodología que
cambia a cada instante, también se toma propone [17]. En este documento se
en cuenta que las líneas eléctricas de un calcula el balance de calor y la
extremo a otro no son totalmente rectas temperatura que puede soportar un
existe una curva en las líneas eléctricas a conductor, según estos artículos [18]–
la cual se la denomina flecha, esa sería la [20] realizan otros métodos para calcular
curva principal de las redes eléctricas, el balance de calor los cuales son muy
estas se producen por el peso del parecidos al artículo mencionado
conductor y aunque los conductores estén anteriormente.
ajustados correctamente va a existir la
curva (flecha) como se observa en la Un ciclo que se da en el conductor es
figura 8, cuando se habla de que las cuando presenta una expansión en la
variables que afectan a la temperatura del longitud del mismo esto causa el aumento
conductor son estocásticas aparecen de la flecha en el conductor por lo cual
diferentes ciclos térmicos los cuales esto puede romper la distancia de
pueden ocasionar efectos en el conductor disrupción entre los conductores de
eléctrico. diferente fase o con cualquier estructura
aterrizada y modificando la estructura del
conductor, existe una relación entre la
curva característica del conductor y la
8
temperatura esto se puede observar en la 2.4. Comportamiento Térmico en
figura 8 y 9. los conductores.
Esto se puede determinar por la ecuación

de balance de calor en las cuales se
encuentran descritas anteriormente,
existe un estado estable o un estado
estacionario el cual se determina o se
controla la corriente que debe pasar por el
conductor según las características
propias del mismo, según el método
seleccionado no toma en cuenta algunos
problemas o factores en el conductor, por
ejemplo no se toma en cuenta el
Figura 9. Estado de la curva a diferentes calentamiento magnético en el conductor
temperaturas (calefacción magnética) lo cual debido al
aumento de calor causa que los
Otro ciclo térmico es cuando por el conductores eléctricos aumenten su
conductor atraviesa ciclos repetitivos de oposición al paso de corriente al menos
corriente por un transcurso de tiempo 2% y el otro efeto que no se considera es
prolongado ocurren daños en la el efecto corona pero este solo afecta
composición física del conductor por lo cuando existen disrupciones eléctricas en
que su vida útil disminuye. las uniones o terminales de los
conductores, partiendo de estos
Al determinar la cargabilidad en las redes parámetros se calculara el equilibrio
de distribución o transmisión se debe térmico en los conductores con (11).
calcular la temperatura del conductor Todos los componentes térmicos se
eléctrico (Tc) cuando existen variaciones describen de mejor manera en la tabla 1.
de corriente que circulan por el 1
[𝐼 2 𝑅
(𝑇𝑐) + 𝑞 − 𝑞 − 𝑞 ] = 0 (11)
𝑠 𝑐 𝑟
conductor, también es indispensable 𝑚𝐶𝑝
conocer la distancia que aumenta en la Donde:
línea cuando la temperatura aumenta, se 𝑚𝐶𝑝 es la capacidad térmica del
determina la distancia también por las conductor, es la relación del calor y masa
características propias del conductor. del conductor [𝐽/𝑚°𝐶]
También se debe tomar en cuenta si el 𝑅(𝑇𝑐) es la relación de temperatura y la
terreno es plano o existe inclinaciones en resistencia del conductor, sus unidades
el terreno esto también influye en la son [Ω/𝑚]
distancia que crecerá la flecha otro factor 𝐼es la intensidad De electrones que
que se debe de considerar para la circula por un conductor[A]
distancia que va a tomar la flecha es el 𝑞𝑠 es la suma de calor por radiación en
tiempo en el que ha estado operando el [W/m]
conductor ya que esto influye en las 𝑞𝑐 es la disminución de calor por
características propias del conductor lo convención en [W/m]
que significa más oposición al paso de
𝑞𝑟 es la disminución de calor por
corriente, más temperatura y mayor radiación en [W/m]
tensión de ruptura[21], [22].
9
Tabla 1. Variables del balance térmico. cuenta variables como como la altitud, la
Unidad de azmitud, el área, la absortividad y
Variable Descripción medida
principalmente la radiación solar que
Θ Ángulos de incidencia de [°] grados incide en el conductor, algunas de estas
los rayos del sol variables no son totalmente
independientes, algunas dependen del
𝐻𝑐 Altitud del sol [°] grados
sistema horario y en el lugar donde se
realiza el estudio. Las pérdidas de calor
𝑍1 Azmitud de la línea [°] grados que se producen en el conductor son por
convección 𝑞𝑐 y por radiación 𝑞𝑟 , las
pérdidas causadas por convención pueden
α Absortividad solar []
ser naturales o forzadas estas están
relacionadas con el enfriamiento natural
A’ Área proyectada por [𝑚2 /𝑚] del aire y se define a las pérdidas de calor
unidad de longitud
natural como aquel en donde la velocidad
𝑄𝑆 Radiación solar total [𝑊/𝑚2 ]
del viento es igual a cero en cambio la
forzada es cuando existen velocidades de
viento mayores a cero. Se calculan la
𝑝𝑆 Densidad del aire [𝐾𝑔/𝑚3] temperatura y la corriente que se
producen en la elongación del conductor
D Diámetro del conductor [𝑚𝑚] para encontrar una a proximidad en la
resistencia lineal debido a la temperatura
y se las puede calcular con (12)
𝑇ℎ𝑖𝑔ℎ Temperatura máxima [°𝐶]
𝑏0 𝑎0 𝑇0 + 𝑏0 𝑎 𝑡 (12)
𝑇𝑐 (𝑡) = − + 𝑒 0
𝑇𝑎 Temperatura ambiente [°𝐶] 𝑎0 𝑎0
Los coeficientes 𝑎0 y 𝑏0 se los calcula con

𝑇𝑓𝑖𝑙𝑚 Temperatura del aire [°𝐶]
(13) y (14)
𝛿𝐼 2 (13)
𝑢𝑓 [𝑃𝑎. 𝑆] 𝑎0 =
Viscosidad del aire
𝑚𝐶𝑝
𝐾𝑓 Conductividad térmica [𝑊/𝑚°𝐶] 1 (14)

𝑏0 = [𝛽𝐼 2 + 𝑞𝑠 − 𝑞𝑐 − 𝑞𝑟 ]
del aire 𝑚𝐶𝑝
𝑉𝑤 Velocidad del viento [𝑚/𝑠]

Y los coeficientes 𝛿 y 𝛽 se los calcula con
(15) y (16)
(𝑅𝐻 − 𝑅𝐿 ) (15)
𝛿=
Según los artículos [18], [23] habla (𝑇𝐻 − 𝑇𝐿 )
acerca de las ganancias y las pérdidas de
calor, las ganancias de calor se dan 𝛽 = 𝑅𝐿 + 𝛿𝑇𝐿 (16)
cuando aumenta la temperatura del
conductor debido a la radiación solar y al Donde:
paso de corriente por el conductor. Las 𝑅𝐻 𝑦 𝑅𝐿 es la resistencia del conductor
ganancias de calor que se producen por la medida a diferente temperatura.
radiación solar que absorbe el conductor 𝑇𝐻 𝑦 𝑇𝐿 es diferente valor de temperatura
se la conoce como 𝑞𝑠 y para calcularla se del conductor
necesita saber el ángulo de radiación
hacia el conductor, para encontrar la
ganancia de calor se deben tener en
10
Y (17) es para calcular la corriente que Tabla 2: Variables del balance térmico
circula por el conductor cuando existe un Unidad de
Variable Descripción medida
cambio de temperatura.
(17)
𝑎0 𝑚𝐶𝑝
𝐼=√
𝛿 ⧍𝑦 Desnivel entre apoyos [𝑚]
2.5. Pandeo en las Líneas 𝑆 Separación entre [𝑚]

apoyos
El pandeo en las redes eléctricas se da

𝑇𝑒𝑛 Tensión longitudinal [𝐾𝑔]
cuando existe una extensión a lo largo del
conductor, por lo cual esto se puede
conocer con la temperatura y la cantidad 𝐻 Tensión horizontal [𝐾𝑔]
de corriente que circula por el conductor 𝐿 Longitud del [𝑚]
eléctrico, estos dos parámetros se conductor
encuentran en el literal anterior. La
configuración física del conductor se ve 𝑊 Peso por unidad de [𝐾𝑔/𝑚]
longitud del conductor
expuesta a cambios debido al aumento de
temperatura en el conductor, va ocasionar
una variación en la longitud del cable por
lo que esta variación de longitud tal vez 𝑇 Temperatura de [°𝐶]
conductor
no puede cumplir con las distancias
mínimas de seguridad, para determinar 𝑂 Punto más bajo del []
esta variable de longitud también se debe conductor
saber cuál es el voltaje de disrupción entre
conductores de diferente fase y la 𝑆𝐴 Distancia del apoyo A [𝑚]
al punto O
estructura en la que se encuentra sujeto,
esto se va a determinar a través de (18), 𝑆𝐵 Distancia del apoyo B [𝑚]
todas las variables mencionadas al punto O
anteriormente se describen mejoren la
tabla 2, en donde se utilizan funciones 𝐿𝐴 Longitud del [𝑚]
conductor desde apoyo
hiperbólicas para encontrar todos los A al punto O
aspectos necesarios según los siguientes
artículos [18], [24].
𝐿𝐵 Longitud del [𝑚]
conductor desde apoyo
B al punto O
𝑦𝐴 Altura desde el punto [𝑚]

de apoyo de A hasta al
punto O
𝑌𝐵 Altura desde el punto [𝑚]

de apoyo de B hasta al
punto O
𝐷 Flecha [𝑚]
𝑋𝐹 Distancia horizontal [𝑚]

Figura 10. Representación gráfica de las
variables del balance térmico
11
Para determinar la extensión térmica del Como último paso para determina la
conductor se realiza un modelo variación de longitud que se produjo en el
geométrico en un punto cualquiera del conductor, se realiza cálculos de la
conductor para evaluar la elongación longitud desde el punto A al punto O y
horizontal del mismo más conocido como desde el punto B al punto O se los calcula
(x) y se lo calcula por (18). a través de las siguientes ecuaciones.
𝐻 𝑊𝑆𝐴
ℎ(𝑥) = ℎ𝐴 − [cosh ( ) (18) 𝐿 = 𝐿𝐴 + 𝐿𝐵 (26)
𝑊 𝐻
𝑊(𝑆𝐴 − 𝑥
+ cosh ( )]
𝐻 𝐻 𝑊𝑆𝐴 (27)
𝐿𝐴 = sinh ( )
𝑊 𝐻
Donde:
𝑆𝐴 es el apoyo que se encuentra en el nivel 𝐻 𝑊𝑆𝐵 (28)
mas bajo, y se lo determina por el 𝐿𝐵 = sinh ( )
𝑊 𝐻
siguiente criterio.
3. Estimación de Conductores en
𝑆 𝑆𝑖 ⧍𝑦 ≥ 0 (19)
𝑆𝐴 = { 𝑃𝐵 } Distribución
𝐴 − 𝑆𝑃𝐵 𝑆𝑖 ⧍𝑦 < 0
La estimación de los conductores

Y 𝑆𝑃𝐵 se determina mediante (20) eléctricos en las líneas de distribución se
⧍𝑦
los realiza a través de la corriente que va
𝑆𝐴 𝐻 2 (20) a circular por los conductores, en los
𝑆𝑃𝐵 = − sinh−1 [ ]
2 𝑊 𝐻 𝑊𝑆
sinh−1 ( 𝐴 ) casos donde el sistemas es aéreo el
𝑊 𝐻2
conductor recomendado es un cable
Una vez que se calcula el modelo desnudo que se tiende por todo el sistema,
geométrico que toma el conductor al por lo general el cable que se tiende por
aumentar la temperatura en el mismo se el sistema es de un mismo calibre para
procede a calcular la flecha con (23) que evitar los empalmes o puntos de conexión
se produce en el conductor eléctrico y aun sabiendo que no en todas las partes
para encontrar el valor de la flecha es del sistema se necesita el mismo calibre
necesario determinar la pendiente 𝑚 que del conductor, en los casos donde se
se forma con (21) y la distancia horizontal requiere que el sistema sea soterrado se
entre el puto O y D con (22). necesita conductores con un tipo de
aislamiento especial que soporte otros
⧍𝑦 (21) tipos de parámetros para que funcione
𝑚=
𝑆 adecuadamente el sistema[25].
𝐻 ⧍𝑦 (22)
𝑥𝑓 =
̅̅̅ sinh−1 ( )
𝑊 𝑆 3.1. Selección de los Conductores
Eléctricos.
𝐷 = 𝑦𝑟 (𝑥
̅̅̅)
𝑓 − 𝑦𝑐 (𝑥
̅̅̅)
𝑓 (23)
Para determinar el calibre de los
En donde los valores 𝑦𝑟 y 𝑦𝑐 se los conductores que se van a colocar en el
determina de la siguiente manera. sistema, como primer paso de debe
calcular la corriente que va a circular por
𝐻 ̅̅̅
𝑊𝑆 el conductor con los datos del sistema a la
𝑦𝑟 (𝑥̅𝑓 ) = [cosh (
𝐵
) − 1] − 𝑚(𝑆𝐵 (24)
𝑊 𝐻 cual se la va a denominar 𝐼𝑠 , que será la
− 𝑥𝑓 )
corriente del sistema, después se elige la
𝐻 ̅̅̅
𝑊𝑆
corriente nominal la cual se la denomina
𝐵 (25)
𝑦𝑐 (𝑥̅𝑓 ) =
𝑊
[cosh (
𝐻
) − 1] como 𝐼𝑛 y por ultimo según tablas ya
12
establecidas por los fabricante de desnudo que se debe colocar. Se debe
conductores se elige una 𝐼𝑐 la cual debe destacar que se coloca en la mayoría de
ser mayor o igual a la 𝐼𝑛 [26],[27]. sistemas un conductor del mismo calibre
para casi todos los tramos por lo cual se
Con estos parámetros establecidos se ha decidido realizar el estudio de
debe cumplir (1). dilatación térmica en el primer tramo de
cada sistema (13 y 34 nodos) porque en el
𝐼𝑆 =< 𝐼𝑛 =< 𝐼𝑐 (1) análisis de este, los resultados de
corriente de los dos sistemas de
El conductor seleccionado debe ser distribución se visualizaran en las tablas
desnudo y tiene que ser capaz de soportar 3 y 4, el tramo principal es el que va a
la corriente que vaya a pasar por este con soportar la mayor corriente debido a la
un cierto margen de tolerancia. sumatoria de corrientes de todas las
ramificaciones del sistema y como es el
3.2. Determinación de la Corriente mismo calibre de conductor para todo el
por un Flujo de Carga. sistema el tramo principal será en donde
se genera la mayor temperatura en todo el
Al realizar un flujo de carga en sistemas conductor[11],[31]–[34].
eléctricos los resultados que dan son de
voltaje y potencia por lo cual a través de Tabla 3: Alimentador de 34 nodos.
estos parámetros se puede determinar la Descripción Valor Unidad
corriente que circula de los diferentes Nodo A 800 -----
tramos del sistema, en (2) se puede Nodo B 802 -----
despejar el valor de la corriente[28]–[30]. V1 1 p.u.
V2 0.2622+j0.3729 p.u
𝑆1 = 𝑉1 ∗ 𝐼1∗ (2) Z serie 1.0259+j1.0239 p.u.
S 0.3098-j0.0543 p.u
Partiendo de (2) se puede despejar el I 1.1387+j0.3729 p.u.
valor de la corriente y así poder Tabla 4: Alimentador de 13 nodos.

Descripción Valor Unidad
determinar el calibre del conductor como
Nodo A 850 -----
se ve en (3) y (4).
Nodo B 832 -----
𝑆1∗
V1 1 p.u.
(3)
𝐼1 = ∗ V2 0.8317+j0.0908 p.u
𝑉1 Z serie 0.431+j1.2128 p.u.
S 0.2516-j0.7389 p.u
𝑆1∗ (4) I 0.1827+j0.1649 p.u.
= 𝑉1 𝑌𝑠ℎ + (𝑉1 − 𝑉2 )𝑌𝑠𝑒𝑟
𝑉1∗
En la tabla 3 y 4 se puede observar los
Como en los sistemas radiales resultados obtenidos en los tramos
seleccionados no existe un efecto principales de los alimentadores de 13 y
capacitivo en las líneas debido a 34 nodos de la IEEE obtenidos mediante
diferentes factores por lo cual las la modificación del método de Gauss-
impedancias en derivación no existen la Seidel como se describe en las tablas los
ecuación quedaría de la siguiente manera. resultados del voltaje y de la potencia son
necesarios para determinar la corriente
𝐼1 = (𝑉1 − 𝑉2 )𝑌𝑠𝑒𝑟 (5) que pasa de un nodo hacia otro. Con el
resultado de la corriente ya se puede
Con (5) se determinaría la corriente que determinar el calibre de los conductores
pasa de un nodo a otro y con esta ya se idóneos que deben ser colocados en los
puede determinar el calibre de conductor sistemas[11],[35]–[38].
13
temperatura se debe de considerar
𝑆𝐵3𝛷 (6) fundamentalmente la resistividad por lo
𝐼𝐵 = que esta indica si es un buen conductor
√3𝑈𝐵𝐿𝑖𝑛𝑒𝑎
por ejemplo la ρ del cobre es de
𝐼𝑟𝑒𝑎𝑙= 𝐼𝐵 ∗ 𝐼 (7) 1.7𝑋10−8 Ohm − m. Con (8) se puede
determinar la resistencia cuando la
Donde: temperatura aumenta.
𝑆𝐵3𝛷 es la potencia real en VA.
𝑈𝐵𝐿𝑖𝑛𝑒𝑎 es el voltaje base en V. 𝑅= 𝑅𝑂 ∗ (1 + 𝛼𝑡) (8)
𝐼 es la corriente calculada por el flujo en
p.u. Donde:
𝐼𝐵 es la corriente base en A. 𝑅𝑂 es la resistencia a una temperatura
inicial
Para obtener el valor real de la corriente 𝛼 es el coeficiente de temperatura en la
se va a utilizar (6), (7) y así determinar un resistencia.
conductor que cumpla con las
características requeridas. Si se necesita calcular la resistividad del
conductor a una temperatura diferente a
3.3. Efecto del Exceso de Corriente 20°C se utiliza (9) la cual sirve para
y Temperatura en los determinar a cualquier temperatura a
Conductores diferencia de (8) que sirve solo cuando la
temperatura es a 20°C.
La resistencia al paso de corriente que
tienen los conductores eléctricos depende 𝜌𝛳 = 𝜌20°𝐶 + 𝜌20° ∗ 𝛼(𝑇 − 20°) (9)
de la temperatura a la que se encuentre el
conductor ya que la resistividad (r) va a Donde:
cambiar dependiendo a la temperatura 𝜌𝛳 es resistividad a la temperatura
que este se encuentre. La variación de requerida.
temperatura es proporcional a la corriente 𝛼 es coeficiente de temperatura.
que pasa por el conductor por lo cual si la 𝜌20°𝐶 es resistividad a 20°C.
corriente aumenta la temperatura 𝑇 es la temperatura a la que se desea
aumenta y va a existir un cambio en la calcular.
resistencia (R) la cual va a limitar la
corriente que va a pasar por dicho 4. Metodología Para Determinar
conductor. De forma común la Flujos en los Sistemas de 13 y 34
temperatura que suele estar dada para la nodos.
resistividad es de 20° C y se denomina
(α) la variación de impedancia por
Gauss-Seidel es uno de los método para
temperatura y corriente[19],[39], [40].
determinar flujos de potencia en líneas de
transmisión, pero en este artículo se
Para evitar daños o fallas en el sistema se
realizó la modificación del método y se lo
recomienda que a mayor corriente debe
uso en líneas de distribución acoplando el
ser mayor la sección del conductor
método según [1], [16],[41] con estos
expresada en 𝑚𝑚2 , además el documentos se puede determinar un
calentamiento (perdidas) o la elevación resultado efectivo en líneas de
de temperatura se producen por que la distribución realizando flujos de potencia
estructura física de los conductores no es entre dos nodos ya que así la
perfecta por lo tanto también causa caídas convergencia de los datos es más acertada
de voltaje en los terminales del y se eliminaría el conflicto con la
conductor. Para controlar la elevación de modificación en el método de Gauss-
14
Seidel en líneas de distribución por los parámetros que se describió para el
resultados no exactos o por la excesiva sistema de 13 nodos de la IEEE, con los
memoria que con lleva usar el método en cuales se va a determinar el calibre del
distribución. conductor y se va a realizar el estudio
desde que corriente va a empezar una
650
extensión térmica en los conductores.
645 632 633 634 Antes de realizar el análisis de la

extensión térmica es necesario
representar los valores que dio el flujo de
684 671 692 675 potencia en sus valores reales y así
encontrar cual es la sección de conductor
adecuado para soportar esa corriente sin
que exista la dilatación en el conductor y
652 680
determinar desde que valor de corriente
Figura 11. Sistema de alimentación de incrementa la temperatura en el
13 nodos de la IEEE conductor y empieza la extensión térmica
en el cable.
El estudio de flujos entre dos nodos se los
realizo en los sistemas de distribución de
13 y 34 nodos de la IEEE. En el sistema 4.1. Algoritmo Para la Aplicación
de 13 nodos, figura 11, se realizaron 12 de Gauss-Seidel en
flujos de potencia de nodo a nodo, cada Distribución.
flujo consiguiente que se realizo fue con
los datos ya calculados del anterior para Con la finalidad de determinar el límite de
que el resultado sea más confiable y no cargabilidad en las líneas de los sistemas
haya exceso de memoria. Para obtener alimentadores de 13 y 34 nodos de la
resultados de voltaje y potencia con lo IEEE, se implementó el algoritmo de
que ya se puede realizar el cálculo para flujo de potencia modificado para dos
encontrar la corriente y así determinar el nodos, mismo que considera condiciones
calibre del conductor. de un flujo de potencia normal claro que
848 evitando o excluyendo la parte capacitiva
822 846 del sistema ya que en sistemas de
652
820 844
distribución no existe la parte capacitiva,
resistencia en derivación o porque su
864 842
818
valor es muy bajo para ser considerada en
802
806
808
812
814 850
816
824
826
858
834
860 836 840
el flujo de potencia.
800
862
832
888 890
810
838
Algoritmo del Flujo de Potencia
para el Sistema de 13 y 34 Nodos de
852
la IEEE.
828 830 854 856
Figura 12. Sistema de alimentación de Paso 1: Definir los datos de los

34 nodos de la IEEE. alimentadores de 13 y 34 nodos de
la IEEE.
En el sistema de 34 nodos de la IEEE, 𝑉𝑖+1 :Voltaje del primer nodo.
figura 12, se realizaron 33 flujos de 𝑉𝑗 :Voltaje del segundo nodo.
potencia debido a que los flujos se los 𝑍𝑠𝑒𝑟𝑖𝑒 : Impedancia serie entre los nodos.
realiza de cada uno de los tramos del 𝑃𝑖 : Potencia Activa de cada nodo.
sistema, para determinar los mismos 𝑄𝑖 : Potencia Reactiva de cada nodo.
15
𝑆𝑘∗ : Potencia aparente conjugada. aplicándole una modificación en el
algoritmo para que no exista un excesivo
Paso 2: Definir las matrices de almacenamiento de datos y los resultados
admitancias para dos nodos. puedan llegar a una convergencia más
𝑌 𝑌 rápida. Estos resultados se los comparará
𝑌𝐵𝑈𝑆 = [ 11 12 ]
𝑌21 𝑌22 con los resultados que se obtengan en la
simulación, con el software DIgSILENT
Paso 3: Realización del flujo. PowerFactory se verificará o validará la
j=2 modificación del algoritmo.
for i=1:n
𝑛
𝑆𝑘∗ 650
− ∑ 𝑌𝑘𝑖 𝑉𝑖
𝑉𝑘∗
𝑖=1 646 645 633 634
𝑖≠𝑗 300 ft 500 ft 500 ft XFM-1
conf 603 conf 603 conf 602 S=160+110i
end 632
2000 ft
𝑛 conf 601
1 𝑆𝑘∗
𝑉𝑘 = − ∑ 𝑌𝑘𝑖 𝑉𝑖
𝑌𝑘𝑘 𝑉𝑘∗
684 692
300 ft switch
𝑖=1 conf 604
( 𝑖≠𝑗 ) 611 671
675
300 ft 2000 ft
500 ft
conf 605 conf 601 conf 606
Paso 4: Cálculo de la Potencia
Activa y Reactiva en cada nodo. 652 680
𝑛 800 ft 1000 ft
conf 607 conf 601
𝑃𝑖 = ∑|𝑉𝑖 𝑌𝑖𝑘 𝑉𝑘 | cos(𝛿𝑘 − 𝛿𝑖 + 𝛾𝑖𝑘 ) Figura 13. Simulación del flujo del
𝑘=1 alimentador de 13 nodos de la IEEE.
𝑛
𝑄𝑖 = ∑|𝑉𝑖 𝑌𝑖𝑘 𝑉𝑘 | sin(𝛿𝑘 − 𝛿𝑖 + 𝛾𝑖𝑘 ) En la simulación del alimentador de 13

𝑘=1 nodos de la IEEE se puede comparar los
resultados obtenidos en el flujo de
potencia realizado en Matlab con los
Paso 5: Cálculo de la Corriente en resultados obtenidos en Power Factory y
p.u. llegar a determinar si los resultados en
Matlab son iguales o muy parecidos a los
𝐼1 = (𝑉1 − 𝑉2 )𝑌𝑠𝑒𝑟 resultados obtenidos en el simulador, los
cuales.
Paso 6: Determinación de la
corriente a su valor real. Con los resultados verificados ya se
puede empezar a determinar a qué nivel
𝑆𝐵3𝛷 de corriente empieza a existir una
𝐼𝐵 =
√3𝑈𝐵𝐿𝑖𝑛𝑒𝑎 dilatación en el cable y así poder
determinar otros parámetros que se dan
𝐼𝑟𝑒𝑎𝑙= 𝐼𝐵 ∗ 𝐼 por la dilatación en el cable como la
longitud que aumenta y cuanto crece la
Paso 7: Fin. flecha. La comprobación de los
resultados da paso a determinar la
4.2. Verificación de los Resultados corriente que circula por cada tramo del
por DIgSILENT Power sistema, con estos datos se puede
Factory. determinar que calibre de conductor se
debería colocar en cada tramo del sistema
Los flujos de potencia que se realizaron para así obtener una selección adecuada
en los sistemas de 13 y 34 nodo de la de conductores.
IEEE con la metodología de Gauss-Seidel
16
848
530 ft
conf 301
Tabla 5: Tipo de conductor eléctrico del
822
13740 ft
846
3640 ft
sistema de 13 nodos de la IEEE.
conf 302 conf 301
652 Node Node Corrie Confi Fase Neutro
820
49150 ft
844
1350 ft
A B nte(A) g. ACS ACSR
conf 301
conf 302
842
R
818 864
1710 ft 280 ft
conf 302
1620 ft
conf 302
conf 301 632 645 133.2 603 1/0 1/0
808
802 814 860 836
32230 ft 850 824 840
2580 ft 29730 ft 858 2020 ft 2680 ft
conf 300 10 ft V=0.0640 + 0.6125i 860 ft
conf 300 conf 300 4900 ft conf 301 conf 301
conf 301 S=0.0032-0.0011i conf 301
conf 301 834
806 826 5830 ft
812 816
800
1730 ft
conf 300
37500 ft 310 ft
3030 ft
conf 303
conf 301 632 633 208.10 602 4/0 4/0
300 conf 301 862
832 280 ft
10 ft 888 890 conf 301
conf 301 0 ft 10560 ft
810 XFM-1 conf 300
5804 ft
conf 303 838 633 634 0 XFM
4860 ft
conf 304 -1
852
36830 ft
830 conf 301
828
20440 ft 645 646 133.2 603 1/0 1/0
conf 301
840 ft
301 854 856
520 ft 23330 ft
conf 301 conf 303
650 632 319.4 601 4/0 4/0

Figura 14. Simulación del flujo del
alimentador de 34 nodos de la IEEE. 684 652 125.1 607 1/0 1/0
632 671 200.4 601 4/0 4/0

En el sistema de 13 y 34 nodos de la IEEE
Se aprecia la realización de los flujos de
potencia por tramos o flujos de potencia 671 684 100.04 604 1/0 1/0
de dos nodos por lo que en el sistema de

13 nodos se realizaron 12 flujos de 671 680 200.08 601 4/0 4/0
potencia en los cuales se representaron las
cargas de acuerdo a los datos que 671 692 0 Switc
h
especifica el sistema. En el sistema de 34
nodos se realizaron 33 flujos de potencia 684 611 145.34 605 1/0 1/0
en el cual se colocaron los datos de
acuerdo a los datos especificados en el 692 675 0 606 1/0 1/0
sistema.
4.3. Calibre de Conductores En la selección de conductores en el

Seleccionados de Acuerdo al sistema de 13 y 34 nodos de la IEEE se
Nivel de Corriente. toma en cuenta los tramos principales y
los secundarios con los cuales se puede
La selección del calibre del conductor se determinar que calibre es el apropiado
lo realiza con base a datos que se (tablas 5 y 6) y recomendar un mismo
obtuvieron en Matlab y se verificaron en calibre para diferentes tramos y así tener
DIgSILENT Power Factory con lo cual una acertada selección de conductores en
después de validar los datos se determinó el sistema en los nodos donde no existe
la corriente que pasa por cada tramo del una carga el paso de corriente es cero pero
sistema, esto se realizó en el sistema de por recomendaciones de completar el
13 y 34 nodos dando a continuación los sistema se dio un calibre de conductor
calibres que se especifican en las tablas 5 basado en características generales del
y 6 correspondientemente. sistema.
Tabla 6: Tipo de conductor eléctrico del
sistema de 34 nodos de la IEEE.
Nod Nod Corriente( Confi Fase Neutr
eA eB A) g. ACS o
R ACS
R
800 802 133.2 300 1/0 1/0
17
Nod Nod Corriente( Confi Fase Neutr 4.4. Límite de Cargabilidad de los
eA eB A) g. ACS o
R ACS Alimentadores de 13 y 34
R Nodos.
802 806 133.2 300 1/0 1/0
806 808 133.2 300 1/0 1/0 Desde la estructura física del cable
808 810 10.34 303 #4 #4 6/1 eléctrico la potencia transferida máxima
6/1 depende de la relación de temperatura con
808 812 133.2 300 1/0 1/0
la corriente a la cual se la conoce como
812 814 133.2 300 1/0 1/0 límite térmico, en las líneas de
814 850 133.2 301 #2 #2 6/1
transmisión no se permite la circulación
6/1 de energía cercana al límite térmico en la
816 818 100.12 302 #4 #4 6/1
6/1
línea, en cambio en sistemas de
816 824 100.12 301 #2 #2 6/1 distribución en muchas ocasiones estas
6/1 operar al borde del límite térmico de la
818 820 100.12 302 #4 #4 6/1
6/1 línea según (IEEE Std. 738-2002). En
820 822 100.12 301 #2 #2 6/1 este apartado se aplicará una metodología
6/1
824 826 100.12 302 #4 #4 6/1 para determinar el límite térmico en los
6/1 sistemas alimentadores de 13 y 34 nodos,
824 828 133.2 302 #4 #4 6/1
6/1 esto se detalla en la tabla 7 y 8.
828 830 99.4 303 #4 #4 6/1
6/1
830 854 99.4 301 #2 #2 6/1 Tabla 7: Configuración del sistema de
6/1 13 nodos.
832 858 99.4 301 #2 #2 6/1 Config. Fase Neutro
6/1
832 888 99.4 XFM-
ACSR ACSR
1
834 860 80.34 301 #2 #2 6/1
6/1 601 556,500 26/7 4/0 6/1
834 842 80.34 301 #2 #2 6/1
6/1 602 4/0 6/1 4/0 6/1
836 840 80.34 301 #2 #2 6/1
6/1
836 862 80.34 301 #2 #2 6/1 603 1/0 1/0
6/1
842 844 80.34 301 #2 #2 6/1 604 1/0 1/0
6/1
844 846 80.34 301 #2 #2 6/1 605 1/0 1/0
6/1
846 848 80.34 301 #2 #2 6/1
6/1
850 816 133.2 301 #2 #2 6/1
6/1
852 832 80.34 301 #2 #2 6/1
6/1
854 856 80.34 303 #4 #4 6/1
6/1
854 852 80.34 301 #2 #2 6/1
6/1
858 864 15.88 302 #4 #4 6/1
6/1
858 834 80.34 301 #2 #2 6/1
6/1
860 836 80.34 301 #2 #2 6/1
6/1
862 838 10.14 304 #2 #2 6/1
6/1
888 890 10.24 300 1/0 1/0
Figura 15. Corriente proporcional a

temperatura en el sistema de 13 nodos de
la IEEE.
18
Se encontró una relación entre el nivel de nivel de temperatura va a existir una
corriente y la temperatura, en la figura expansión térmica en los conductores del
(13) se puede apreciar que si la corriente sistema de 34 nodos.
aumenta en el conductor la temperatura
también sufre un aumento proporcional 4.5. Expansión Térmica del
después del análisis en el aumento de Sistema Alimentador de 13 y 34
temperatura se verificara a que Nodos de la IEEE.
temperatura empieza la dilatación térmica
en los conductores. Según (IEEE Std. 738-2002) sobre la
capacidad térmica en los conductores, el
Tabla 8: Configuración del sistema de factor de sobrecorriente en los
34 nodos. conductores en sistemas de distribución
Config. Fase Neutro pueden operar a un 200% de su capacidad
sin que haya problemas o exista una
ACSR ACSR dilatación en el conductor, con lo cual se
va a determinar la dilatación térmica
300 1/0 1/0 haciendo referencia a la temperatura
cuando se duplica el valor de la corriente
301 #2 6/1 #2 6/1 nominal que puede pasar por el
conductor. En la figura (15) y (16) se
302 #4 6/1 #4 6/1 puede observar hasta qué punto de
corriente puede pasar por el conductor
303 #2 6/1 #2 6/1 antes y después de que exista una
extensión térmica en el mismo, se tomó
en consideración todos los parámetros
para determinar el límite de corriente que
puede pasar por los sistemas de 13 y 34
nodos de la IEEE.
Figura 16. Corriente proporcional a

temperatura en el sistema de 34 nodos de
la IEEE.
En la figura (14) se puede apreciar la Figura 17. Punto crítico promedio en el
relación entre corriente y temperatura en sistema de 13 nodos de la IEEE.
los conductores del sistema de 34 nodos
de la IEEE, en esta relación se determina
la proporcionalidad del aumento de
corriente y el aumento de temperatura en
el conductor y así poder determinar a qué
19
848
V=0.3419 + 0.2501i
S=0.3779+0.4759i
822 846
V=0.1558 + 0.1439i V=-0.5681 + 0.2159i
S=0.2794-1.4919i S=0.0794+0.2078i
652
820 844
V=1.0788 + 1.0678i V=0.5792 + 0.7891i
S=0.0587-0.2906 S=0.0991-0.0730i
842
818 864 V=0.1987 + 0.2710i
V=0.8572 + 0.2144i V=0.2858 + 0.1724i S=0.1577-0.3811i
S=0.0206-0.1626i S=0.2351-0.6169i
808
802 814 860 836 840
V=3.6426 + 3.0198i 850 824 858
V=-0.2622 + 0.3729i V=-0.9916 + 0.0058i V=0.1841 + 0.5629i V=-0.6015 + 0.1655i V=0.7029 + 0.5141i
S=0.0963-0.1101i V=-0.0308 + 0.0056i V=0.0640 + 0.6125i V=0.5195 + 0.3434i
S=0.3098-0.0543i S=-0.0126+0.5364i S=0.0132+0.1746i S=-0.0106+0.2267i S=0.3400+0.5231i
S=-2.8246-1.3463i S=0.0032-0.0011i S=0.0017-0.0025i
834
806 826 V=-0.6289 + 0.2593i
812 816
V=0.5372 + 0.5679i V=1.6165 + 0.8578i
800 S=0.3444+0.2562i
V=0.6503 + 0.1908i V=0.1600 + 0.1367i
S=0.4713-0.3964
S=0.0920+0.1107i
S=1.3408-1.5240i S=0.0107+0.0244i 862
832 V=0.2313 + 0.2794i
V=0.0043 + 0.0056i 888 890 S=0.1318+0.3864i
S=0.0243+0.8719i V=1.2959 + 2.5654i V=0.7256 + 0.7131i
810 S=0.0045+5.6548e-05i S=2.5630-1.3007i
V=0.7413 + 0.1077i
S=0.4587-4.7093i 838
V=2.5427 + 1.8804i
S=0.1994-0.3030i
852
V=-0.9935 + 0.0047i
830 S=0.0041+0.0128i
V=0.8747 + 1.3407i
828 S=0.0661-0.0792i
V= 2.8221 + 2.1118i
S=0.5221-0.2214i 854 856
V=1.8288 + 1.3504i V=
S=0.9769-1.1851i S=0.0931-2.2715i
Figura 20. Resultados del flujo de

Figura 18. Punto crítico promedio en el potencia del alimentador de 13 nodos de
sistema de 34 nodos de la IEEE. la IEEE.
5. Resultados 5.2. Flujos en DIgSILENT de los
Sistemas de 13 y 34 Nodos.
La metodología propuesta será sometida
a dos casos de aplicación considerando En las figuras 19 y 20, se podrán
que cada caso será enfocado a las normas visualizar los resultados obtenidos por el
IEEE Std. 738-2002, es decir al límite software DIgSILENT Power Factory en
térmico de dos sistemas alimentadores de el cual se pueden verificar que los
13 y 34 nodos respectivamente con una anteriores resultados si fueron los
topología radial. correctos con lo que ya da paso a realizar
el cálculo de la corriente que circula por
5.1. Flujo de Potencia.
cada tramo de los sistemas.
En las figuras 17 y 18, se podrá visualizar 650
el resultado del flujo de potencia en los

dos sistemas alimentadores en donde se 646
V=0.1381 + 0.1399i
645
V= -1.8794 + 0.3330i
633
V=0.0997 + 0.1011i
634
V=1.1000 + 2.0000i
S=1.2786+1.0426i
detalla el voltaje, la potencia reactiva y S=0.9219+0.6148i
632
S=1.5590+1.9857i S=0.1212-0.4073i
activa en cada nodo. Con estos resultados V=-0.8317+j0.0908

S=-0.2516-j0,7389
se podrá determinar la corriente que pasa
684 692
por cada tramo del conductor. V=-0.1242 + 0.0329i
S=1.7483+1.5191i
V=0
S=0
611 675
671
650 V=0.1949 + 0.1977i V=-0.6417 + 0.0531i
V=0.1577 + 0.1250i
S=0.3984+0.2354i S=0.7999+1.7163 S=-0.9496+-2.3076
646 645 633 634

V=0.1381 + 0.1399i V= -1.8794 + 0.3330i V=0.0997 + 0.1011i V=1.1000 + 2.0000i
S=0.9219+0.6148i S=1.2786+1.0426i S=0.1212-0.4073i 652 680
S=1.5590+1.9857i
632 V=0.6394 + 0.2441i V=0.1693 + 0.4977i
V=-0.8317+j0.0908 S=0.1867+0.0198i S=0.4888+1.1231i
S=-0.2516-j0,7389
684 692
Figura 21. Flujo en DIgSILENT del
V=-0.1242 + 0.0329i V=0
S=1.7483+1.5191i S=0 alimentador de 13 nodos.
611 675
671
V=0.1949 + 0.1977i V=-0.6417 + 0.0531i
V=0.1577 + 0.1250i
S=0.3984+0.2354i S=0.7999+1.7163 S=-0.9496+-2.3076
652 680
V=0.6394 + 0.2441i V=0.1693 + 0.4977i
S=0.1867+0.0198i S=0.4888+1.1231i
Figura 19. Resultados del flujo de potencia

del alimentador de 13 nodos de la IEEE.
20
848
V=0.3419 + 0.2501i
S=0.3779+0.4759i
822 846
V=0.1558 + 0.1439i V=-0.5681 + 0.2159i
S=0.2794-1.4919i S=0.0794+0.2078i
652
820 844
V=1.0788 + 1.0678i V=0.5792 + 0.7891i
S=0.0587-0.2906 S=0.0991-0.0730i
842
818 864 V=0.1987 + 0.2710i
V=0.8572 + 0.2144i V=0.2858 + 0.1724i S=0.1577-0.3811i
S=0.0206-0.1626i S=0.2351-0.6169i
808
802 814 860 836 840
V=3.6426 + 3.0198i 850 824 858
V=-0.2725 + 0.3829i V=-0.9916 + 0.0058i V=0.1841 + 0.5629i V=-0.6015 + 0.1655i V=0.7029 + 0.5141i
S=0.0963-0.1101i V=-0.0308 + 0.0056i V=0.0640 + 0.6125i V=0.5195 + 0.3434i
S=0.3188-0.0541i S=-0.0126+0.5364i S=0.0132+0.1746i S=-0.0106+0.2267i S=0.3400+0.5231i
S=-2.8246-1.3463i S=0.0032-0.0011i S=0.0017-0.0025i
834
806 826 V=-0.6289 + 0.2593i
812 816
V=0.5372 + 0.5679i V=1.6165 + 0.8578i S=0.0920+0.1107i
V=0.6503 + 0.1908i V=0.1600 + 0.1367i
S=0.3444+0.2562i S=0.4713-0.3964
S=1.3408-1.5240i S=0.0107+0.0244i 862
832 V=0.2313 + 0.2794i
V=0.0043 + 0.0056i 888 890 S=0.1318+0.3864i
S=0.0243+0.8719i V=1.2959 + 2.5654i V=0.7256 + 0.7131i
810 S=0.0045+5.6548e-05i S=2.5630-1.3007i
V=0.7413 + 0.1077i
S=0.4587-4.7093i 838
V=2.5427 + 1.8804i
S=0.1994-0.3030i
852
V=-0.9935 + 0.0047i
830 S=0.0041+0.0128i
V=0.8747 + 1.3407i
828 S=0.0661-0.0792i
V= 2.8221 + 2.1118i
S=0.5221-0.2214i 854 856
V=1.8288 + 1.3504i V=
S=0.9769-1.1851i S=0.0931-2.2715i
Figura 22. Flujo en DIgSILENT del Figura 24. Punto crítico de la

alimentador de 34 nodos. configuración 602.
5.3. Punto Crítico Antes de que

Exista una Extensión Térmica.
En las figuras 21 a 28, se podrá visualizar

que nivel de corriente es la que tiene que
pasar por los conductores de los sistemas
alimentadores de 13 y 34 nodos de la
IEEE para que el incremento de
temperatura en los mismos sea lo
suficiente alto para que supere su límite
térmico o en otras palabras que empiece a Figura 25. Punto crítico de la
sufrir una expansión térmica en los configuración 603.
conductores de los sistemas.
Figura 26. Punto crítico de la

Figura 23. Punto crítico de la
configuración 605.
configuración 601.
21
Figura 27. Punto crítico de la Figura 30. Punto crítico de la
configuración 300. configuración 303.
6. Conclusiones
El uso del algoritmo de Gauss-Seidel para

flujos de potencia es considerado
inexacto cuando se lo aplica en sistemas
de distribución pero cuando a este
algoritmo se le realiza modificaciones en
la forma de aplicarlo en el sistema es
decir, el algoritmo se aplica de dos nodos
en dos nodos hasta cubrir todo el sistema,
los resultados que da el algoritmo son
aceptables, ya que al realizarlo solo en
dos nodos el exceso de memoria que
Figura 28. Punto crítico de la causa cuando se lo realiza de la forma
configuración 301. tradicional, ya no ocurre aplicando la
forma anteriormente mencionada por lo
cual los resultados convergen más rápido
y dan una solución confiable.
La utilización de un simulador para

determinar los flujos de potencia en los
sistemas alimentadores es importante
para la comparación de resultados y
verificar si los resultados fueron los
correctos, los resultados generados por el
simulador dieron iguales o muy cercanos
a los resultados obtenidos por la
modificación del algoritmo de Gauss-
Seidel aplicado en sistemas de
Figura 29. Punto crítico de la distribución, por lo tanto se llega a la
configuración 302. determinación que la modificación
propuesta dan resultados válidos y sin
exceso de memoria.
22
Los resultados obtenidos por el flujo de 388.
potencia ayudan a determinar la corriente [2] X. M. Huang, Y. J. Zhang, and H.
nominal que pasa por cada tramo de los C. Huang, “Automatic reactive
sistemas por lo que se concluye el calibre power control in distribution
ideal de los conductores que deben ser network based on Feeder Power
colocados en los diferentes tramos de los Factor Assessment,” Asia-Pacific
sistemas de distribución dejando como Power Energy Eng. Conf.
resultados la adecuada selección de APPEEC, vol. 2015-March, no.
conductores eléctricos. March, pp. 2–5, 2014, doi:
10.1109/APPEEC.2014.7066071.
En las redes eléctricas de distribución el [3] P. S. Bhowmik, S. P. Bose, D. V.
factor de sobrecarga que es permitido en Rajan, and S. Deb, “Power flow
estas es de 2 por lo cual se evaluó con analysis of power system using
corrientes más altas que superaban el power perturbation method,”
valor de sobrecarga de los conductores en PEAM 2011 - Proc. 2011 IEEE
los sistemas de distribución y se Power Eng. Autom. Conf., vol. 3,
determinó el punto crítico en el cual la no. 2, pp. 380–384, 2011, doi:
temperatura en el conductor incrementa 10.1109/PEAM.2011.6135117.
tanto que empieza a sufrir cambios la [4] D. P. Chassin, P. R. Armstrong, D.
forma física del conductor, empieza a G. Chavarr??a-Miranda, and R. T.
darse una expansión térmica en el Guttromson, “Gauss-Seidel
conductor. El nivel de corriente tiene que accelerated: Implementing flow
superar al menos dos veces el nivel para solvers on field programmable
el que están diseñados los conductores gate arrays,” 2006 IEEE Power
para que ocurra una extensión térmica en Eng. Soc. Gen. Meet. PES, 2006,
ellos. doi: 10.1109/pes.2006.1709227.
[5] B. A. Carré, “Solution of Load-
7. Trabajos futuros Flow Problems by Partitioning
El estudio del análisis de cargabilidad en Systems into Trees,” IEEE Trans.
las líneas de distribución por la Power Appar. Syst., vol. PAS-87,
modificación del algoritmo de Gauss- no. 11, pp. 1931–1938, 1968, doi:
Seidel, da la apertura de realizar estudios 10.1109/TPAS.1968.292015.
relacionados con la aplicación de otros [6] H. Sato, “Computation of power
flujos de potencia que se los utiliza system loadability limits,” Proc.
únicamente en transmisión para ver si a IEEE Power Eng. Soc. Transm.
estos realizándoles un ajuste en la Distrib. Conf., vol. 3, no. ASIA
estructura o en la forma en la que se PACIFIC, pp. 1707–1711, 2002,
aplican se puede llegar a dar resultados doi: 10.1109/tdc.2002.1177711.
confiables en distribución. [7] S. Li, X. Han, and M. Wang, “A
New Method for Reference
8. Referencias Network Considering Factor,” no.
51407106, pp. 2014–2017, 2016.
[1] F. T. Toledo and M. A. V. León, [8] V. V. Reddy, G. Yesuratnam, and
“Computational algorithm for the M. S. Kalavathi, “Impact of
analysis of loadability in voltage and power factor change
distribution systems,” Proc. 2018 on primary distribution feeder
IEEE 38th Cent. Am. Panama power loss in radial and loop type
Conv. CONCAPAN 2018, pp. 2–7, of feeders,” Proc. - ICETEEEM
2018, doi: 2012, Int. Conf. Emerg. Trends
10.1109/CONCAPAN.2018.8596 Electr. Eng. Energy Manag., pp.
23
70–76, 2012, doi: [16] J. A. Treece, “Bootstrap Gauss-
10.1109/ICETEEEM.2012.64944 Seidel load flow,” Proc. Inst.
46. Electr. Eng., vol. 116, no. 5, p.
[9] K. Peng et al., “Gauss-Newton 866, 1969, doi:
hybrid power flow algorithm for 10.1049/piee.1969.0161.
AC-DC distribution system,” [17] B. S. Mohan, C. B.
China Int. Conf. Electr. Distrib. Shankaralingappa, and R. Prakash,
CICED, vol. 2016-Septe, no. “Fractional frequency
Ciced, pp. 10–13, 2016, doi: transmission system to enhance
10.1109/CICED.2016.7576280. line loadability in EHVAC and
[10] M. T. U. Naik, “Solving UHVAC long transmission line,”
Eigenvalues of Symmetric,” no. 2014 Int. Conf. Adv. Electron.
Iccmc, pp. 922–925, 2017. Comput. Commun. ICAECC 2014,
[11] C. Marcelo and M. Montecinos, pp. 0–5, 2015, doi:
“PARA SISTEMAS DE 10.1109/ICAECC.2014.7002477.
DISTRIBUCION,” Pontificia [18] S. Milena and T. Gutiérrez,
Universidad Católica de Chile, “Comportamiento de conductores
1989. eléctricos usados en líneas de
[12] G. . Gilbert, D.E.Bouchard, and A. transmisión ante esfuerzos
Y. Chikhani, “Comparison of electromecánicos y térmicos
load,” R. Mil. Coll. Canada combinados,” issue addressed this
Kingston, Ontario, K7K 5LO, vol. thesis is Behav. Electr. Conduct.
1, no. 1, pp. 850–853, 1998. used Transm. lines, vol. 2011, no.
[13] S. Mohagheghi, G. K. Stefopoulos, 5, p. 70, 2011.
and M. Parniani, “Reactive [19] L. Coria, “Analisis De Flujos De
compensation techniques for Carga,” vol. 1, no. 1, pp. 117–202,
increasing loadability of long 2015.
primary distribution lines,” Proc. [20] H. Zhang, H. Wang, and S. Zheng,
37th Annu. North Am. Power “Analysis of the thermal expansion
Symp. 2005, vol. 2005, pp. 233– effect on measurement accuracy of
237, 2005, doi: Rogowski coils,” Diangong Jishu
10.1109/NAPS.2005.1560530. Xuebao/Transactions China
[14] F. Romero et al., “Development of Electrotech. Soc., vol. 22, no. 5,
a genetic algorithm for evaluating pp. 18–23, 2007.
the performance of overhead [21] J. G. Michopoulos, M. Young, and
power distribution lines and A. Iliopoulos, “A multiphysics
proposing solutions to improve theory for the static contact of
distribution line safety,” IET Conf. deformable conductors with fractal
Publ., vol. 2013, no. 615 CP, pp. rough surfaces,” IEEE Trans.
10–13, 2013, doi: Plasma Sci., vol. 43, no. 5, pp.
10.1049/cp.2013.0982. 1597–1610, 2015, doi:
[15] W. T. Huang and S. T. Chen, 10.1109/TPS.2015.2416980.
“Sequential three-phase power [22] S. D. Lee, K. Y. Shin, H. J. Song,
flow calculation for radial D. Il Lee, and B. U. Min, “The sag
distribution systems via three- and fatigue properties of
phase Z-bus distribution factor,” STACIR/AW as a high
2009 4th Int. Conf. Innov. Comput. temperature, low sag conductor,”
Inf. Control. ICICIC 2009, no. 1, 2006 Int. Conf. Power Syst.
pp. 647–651, 2009, doi: Technol. POWERCON2006, vol.
10.1109/ICICIC.2009.329. 1, no. 1, pp. 1–7, 2006, doi:
24
10.1109/ICPST.2006.321531. Power Transmission Lines,” IEEE
[23] J. Luo, L. Zhang, and Y. Liu, Lat. Am. Trans., vol. (in press), no.
“Study on the thermal expansion X, pp. 1–8, 2019.
of EHV cable line with large size [30] D. Carrión, E. Inga, J. W.
conductor in tunnel,” 2010 IEEE Gonzalez, and R. Hincapié,
PES Transm. Distrib. Conf. Expo. “Optimal geographical placement
Smart Solut. a Chang. World, vol. of phasor measurement units based
1, no. 1, pp. 1–5, 2010, doi: on clustering techniques,” Proc. -
10.1109/TDC.2010.5484463. 2016 51st Int. Univ. Power Eng.
[24] J. Kwon and K. W. Hedman, Conf. UPEC 2016, vol. 2017-
“Transmission expansion planning Janua, pp. 1–6, 2016, doi:
model considering conductor 10.1109/UPEC.2016.8114003.
thermal dynamics and high [31] D. Carrion, J. W. Gonzalez, I. A.
temperature low sag conductors,” Isaac, G. J. Lopez, and H. A.
IET Gener. Transm. Distrib., vol. Cardona, “Load Characterization
9, no. 15, pp. 2311–2318, 2015, Based on Voltage and Current
doi: 10.1049/iet-gtd.2015.0257. Phasorial Measurements in Micro-
[25] S. Lukasz and R. Waldernar, “The Grids,” Proc. - 2017 Int. Conf. Inf.
International Conference on Syst. Comput. Sci. INCISCOS
Advanced Power System 2017, vol. 2017-Novem, no. 2, pp.
Automation and Protection 1–6, 2018, doi:
Temperature Dependent Dynamic 10.1109/INCISCOS.2017.23.
Loadability Control for [32] D. Carrión, J. W. González, I. A.
Transmission Lines,” pp. 1981– Issac, and G. J. López, “Optimal
1985, 2011. fault location in transmission lines
[26] R. Adapa, M. Shirmohamadi, and using hybrid method,” 2017 IEEE
R. Gray, “SLiM (Sagging Line PES Innov. Smart Grid Technol.
Mitigator) - Mitigating excess sag Conf. - Lat. Am. ISGT Lat. Am.
in overhead transmission lines,” 2017, vol. 2017-Janua, no. 3, pp.
Proc. Inaug. IEEE PES 2005 Conf. 1–6, 2017, doi: 10.1109/ISGT-
Expo. Africa, vol. 2005, no. July, LA.2017.8126757.
pp. 168–174, 2005, doi: [33] D. Carrion, A. Ayo, and J. W.
10.1109/pesafr.2005.1611808. Gonzalez, “Under Frequency Load
[27] D. López and D. Carrión, “Método Disconnection Scheme Based on
Heurístico de Ubicación Óptima Improvement to Semi-Adaptative
de Centros de Transformación y Model,” Proc. - 2019 Int. Conf.
Enrutamiento de Redes Eléctricas Inf. Syst. Comput. Sci. INCISCOS
de Distribución,” Universidad 2019, vol. 1, no. 3, pp. 361–365,
Politécnica Salesiana, Quito, 2016. 2019, doi:
[28] F. Mosquera and D. Carrión, 10.1109/INCISCOS49368.2019.0
“Localización óptima de plantas 0063.
virtuales de generación en sistemas [34] J. Pilatásig Lasluisa and D.
eléctricos de potencia basados en Carrión, “Resiliencia de Sistemas
flujos óptimos de potencia,” I+D Eléctricos de Potencia mediante la
Tecnológico, vol. 16, no. 2, pp. 1– Conmutación de Líneas de
12, 2020, doi: Transmisión – Estado del arte,”
10.33412/idt.v16.2.2827. I+D Tecnológico, vol. 16, no. 2,
[29] S. Pinzón, D. Carrión, and E. Inga, pp. 1–10, 2020, doi:
“Optimal Transmission Switching 10.33412/idt.v16.2.2834.
Considering N-1 Contingencies on [35] J. Palacios and D. Carrión, “Estado
25
del arte de la planeación de
expansión de sistemas de
transmisión,” I+D Tecnológico,
vol. 16, no. 2, p. 8, 2020, doi:
10.33412/idt.v16.2.2835.
[36] P. Escudero and D. Carrión,
“Expansión de un sistema de
transmisión mediante LOPF-AC
Expansion of a transmission
system using LOPF-AC,” Inge
Cuc, vol. 14, no. 2, pp. 116–125,
2018, doi:
10.17981/ingecuc.14.2.2018.11.
[37] A. Sanchez and D. Carrion,
“Modeling of the behavior power
flow on transmission lines based
on voltage and current
synchronopasors,” IEEE Lat. Am.
Trans., vol. 16, no. 4, pp. 1142–
1149, 2018, doi:
10.1109/TLA.2018.8362149.
[38] D. Santiana and D. Carrión,
“Simulación de contingencias
basados en flujos dinámicos de
potencia para sistemas eléctricos
de transmisión.,” Universidad
Politécnica Salesiana, 2019.
[39] J. J. P. Calupiña and D. Carrión,
“Planeación óptima de flujos de
potencia reactiva en sistemas
eléctricos de transmisión basado
en optimización no lineal entera
mixta,” Universidad Politécnica
Salesiana, 2019.
[40] David Alberto Pérez Cruz and D.
Carrión, “Ubicación Óptima De
Dispositivos Facts En Redes
Eléctricas De Transmisión
Mediante Búsqueda Inteligente,”
Universidad Politécnica Salesiana,
2017.
[41] P. A. Achig and D. Carrión,
“Compensación óptima de
potencia reactiva en sistemas
eléctricos de transmisión mediante
búsqueda armónica,” Universidad
Politécnica Salesiana, 2019.
26
5
4
3
2
1
ITEM
AÑO DE PUBLICACIÓN
2014
1969
2017
2018
2005
CUDA
Ms.Teja
TÍTULO DEL ARTÍCULO
Bootstrap Gauss-Seidel load flow

DATOS
Loadability in Distribution Systems

Measurement Precision of Rogowski Coils
Computational Algorithm for the Analysis of

Analysis of the Thermal Expansion Effect on
Implementing the Gauss Seidel Algorithm for
Automatic reactive power control in distribution

Solving Eigenvalues of Symmetric Matrices with
network based on Feeder Power Factor Assessment

8.1. Matriz de Estado de Arte
3
1
1
17
16
CITAS
REDES DE DISTRIBUCION
GAUSS-SEIDEL
FLUJOS DE POTENCI EN REDES RADIALES

TEMÁTICA
PARAMETROS DE CARGABILIDAD
EXPANSION DE LOS CONDUCTORES
ESTIMACION DE CONDUCTORES
SELECCIÓN DE CONDUCTORES
DETERMINACION DE LA CORRIENTE POR MEDIO DE UN

FLUJO
FORMULACIÓN
DEL PROBLEMA
EFECTO DEL EXCESO DE CORRIENTE EN LOS CONDUCTORES
METODOLOGIA
ALGORITMO APLICADO EN DISTRIBUCION
REALIZACION DE GAUSS-SEIDEL EN REDES DE

DISTRIBUCION
CONVERGENCIA DE LOS RESULTADOS EN DIgSILENT
Tabla 9:Matriz del estado del arte.
CALIBRE DEL CONDCUTOR SEGÚN EL NIVEL DE CORRINTE

PROBLEMA
LIMITE DE CARGABILIDAD EN LAS LINEAS

RESTRICCIONES DEL
CONDUCTORES ESPECIFIOS
EXPANSION TERMICA EN LOS CONDUCTORES
ALGORITMO EN MATLAB
FLUJO EN DIgSILENT
DETERMINACION DEL PUNTO CRITICO DE CORRIENTE

MODIFICACIÓN DEL ALGORITMO DE GAUSS-SEIDEL
RELACION TEMPERATURA-CORRIENTE
DETERMINAR EL LIMITE TERMICO EN LOS CONDUCTORES

PROPUESTAS PARA
GAUSS-SEIDEL MODIFICADO
RESOLVER EL PROBLEMA
                   
                     
                    
                 
                    
RESULTADOS DE DIGSILENTEN POWER FACTORY
APLICACIÓN EN MATLAB
VALORES DE VOLTAJE, POTENCIA Y CORRIENTE
VERIFICACION DE CONVERGENCIA DE RESULTADOS

ANÁLISIS DE CARGABILIDAD DE REDES ELÉCTRICAS DE DISTRIBUCIÓN BASADO EN LA
LIMITE DE CORRIENTE ANTES DE LA EXPANSION TERMICA

SOLUCIÓN
PROPUESTA
LIMITE DE CORRIENTE Y TEMPERATURA

27
EFECTOS EN LAS REDES ELECTRICAS
DIFERENTES ALGORITMOS
       
       
       
      
       
FORMULACIÓN RESTRICCIONES DEL PROPUESTAS PARA SOLUCIÓN
DATOS TEMÁTICA  DEL PROBLEMA  PROBLEMA  RESOLVER EL PROBLEMA  PROPUESTA 




RESULTADOS DE DIGSILENTEN POWER FACTORY
LIMITE DE CORRIENTE ANTES DE LA EXPANSION

DETERMINACION DE LA CORRIENTE POR MEDIO
VALORES DE VOLTAJE, POTENCIA Y CORRIENTE


CALIBRE DEL CONDCUTOR SEGÚN EL NIVEL DE

EFECTOS EN LAS REDES ELECTRICAS 

EFECTO DEL EXCESO DE CORRIENTE EN LOS
EXPANSION DE LOS CONDUCTORES 
DETERMINAR EL LIMITE TERMICO EN LOS

DETERMINACION DEL PUNTO CRITICO DE

CONVERGENCIA DE LOS RESULTADOS EN



ESTIMACION DE CONDUCTORES 
VERIFICACION DE CONVERGENCIA DE
SELECCIÓN DE CONDUCTORES 
DIFERENTES ALGORITMOS 


FLUJO EN DIgSILENT 
ALGORITMO EN MATLAB










METODOLOGIA
CONDUCTORES
CONDUCTORES
DISTRIBUCION
GAUSS-SEIDEL
RESULTADOS
DE UN FLUJO
CORRIENTE
CORRINTE
DIgSILENT
TERMICA
CITAS
ITEM
6 2006
Gauss-Seidel Accelerated: Implementing Flow
Solvers on Field Programmable Gate Arrays
15                      
7 2006
Power Flow Analysis of Power System Using Power
Perturbation method
12                        
8 1968
Solution of Load-Flow Problems by
Partitioning Systems into Trees
12                        
9 2002 Computation of Power System Loadability Limits 12                          
A New Method for Reference Network Considering
10 2016 Contingent Events Based on Line Outage 3                        
Distribution Factor
Sequential Three-Phase Power Flow Calculation for
11 2009 Radial Distribution Systems via Three-Phase Z-Bus 3                            
Distribution Factor
Impact ofVoltage and Power Factor Change on

12 2012 Primary Distribution Feeder Power Loss in Radial 4                       
and Loop Type ofFeeders
DEVELOPMENT OF A GENETIC ALGORITHM

FOR EVALUATING THE
13 2013
PERFORMANCE OF OVERHEAD POWER
DISTRIBUTION LINES AND PROPOSING
0                         
SOLUTIONS TO IMPROVE DISTRIBUTION
LINE SAFETY
14 2016
Gauss-Newton Hybrid Power FlowAlgorithm for
AC-DC Distribution System
1                        
15 1989
FLUJO DE POTENCIA TRIFASICO PARA
SISTEMAS DE DISTRIBUCION
19                   
28
FORMULACIÓN RESTRICCIONES DEL PROPUESTAS PARA SOLUCIÓN
DATOS TEMÁTICA  DEL PROBLEMA  PROBLEMA  RESOLVER EL PROBLEMA  PROPUESTA 





EFECTOS EN LAS REDES ELECTRICAS 

EFECTO DEL EXCESO DE CORRIENTE EN LOS
CALIBRE DEL CONDCUTOR SEGÚN EL NIVEL

EXPANSION DE LOS CONDUCTORES 
DETERMINAR EL LIMITE TERMICO EN LOS

DETERMINACION DEL PUNTO CRITICO DE

CONVERGENCIA DE LOS RESULTADOS EN
DETERMINACION DE LA CORRIENTE POR

ESTIMACION DE CONDUCTORES 
RESULTADOS DE DIGSILENTEN POWER

SELECCIÓN DE CONDUCTORES 
VERIFICACION DE CONVERGENCIA DE
LIMITE DE CORRIENTE ANTES DE LA

VALORES DE VOLTAJE, POTENCIA Y

DIFERENTES ALGORITMOS 



FLUJO EN DIgSILENT 

ALGORITMO EN MATLAB




EXPANSION TERMICA



MEDIO DE UN FLUJO



METODOLOGIA
CONDUCTORES
CONDUCTORES
DISTRIBUCION
GAUSS-SEIDEL
DE CORRINTE
RESULTADOS
CORRIENTE
CORRIENTE
DIgSILENT
FACTORY
CITAS
ITEM
A COMPARISON OF LOAD FLOW ANALYSIS

16 1998 USING DISTFLOW, GAUSS-SEIDEL, AND 42                          
OPTIMAL LOAD FLOW ALGOFUTHMS
Reactive Compensation Techniques for Increasing

17 2005
Loadability of Long Primary Distribution Lines
2                     
18 2015 ANALISIS E FLUJOS DE CARGA 1                           

Comportamiento de conductores eléctricos usados
19 2011 en líneas de transmisión ante esfuerzos 8                      
electromecánicos y térmicos combinados
The Sag and Fatigue Pn)perties of STACIR/AW as a
20 2006 High 13                          
Temperature, Low Sag Conductor
                    
Study on the Thermal Expansion of EHV Cable Line
21 2010
with Large Size Conductor in Tunnel
1       
Transmission expansion planning model considering
22 2015 conductor thermal dynamics and high temperature 12                        
low sag conductors
                 

Temperature Dependent Dynamic Loadability
23 2011
Control for Transmission Lines
1        
SLiM (Sagging Line Mitigator) – Mitigating Excess
24 2005
Sag in Overhead Transmission Lines
3                        
Fractional Frequency Transmission System to
25 2015 Enhance Line Loadability in EHVAC & UHVAC 2                        
Long Transmission Line
Transmission management for congested power
26 2016 system: A review of concepts, technical challenges 3                       
and development of a new methodology
 CANTIDAD: 14 12 17 12 15 14 15 6 2 3 17 3 15 24 12 15 6 1 3 2 10 1 1 1 19 5 1 1 1 3 1
29
8.2. Resumen de Indicadores
TEMÁTICA
REDES DE
DISTRIBUCION
14
12
10
COMPORTAMIENTO 8 GAUSS-SEIDEL
TERMICO Y PANDEO 6
4
2
0
EXPANSION DE LOS FLUJOS DE POTENCI

CONDUCTORES EN REDES RADIALES
PARAMETROS DE
CARGABILIDAD
Figura 31. Resumen e indicador de la temática - Estado del arte.
FORMULACIÓN DEL PROBLEMA

10
0
ESTIMACION DE SELECCIÓN DE DETERMINACION DE EFECTO DEL EXCESO PANDEO EN LAS
CONDUCTORES CONDUCTORES LA CORRIENTE POR DE CORRIENTE EN LINEAS
MEDIO DE UN FLUJO LOS CONDUCTORES
Figura 32. Indicador de formulación del problema - Estado del arte.
30
SOLUCIÓN DADA MEDIANTE
14
12
10
0
APLICACIÓN EN VALORES DE VERIFICACION DE LIMITE DE LIMITE DE EFECTOS EN LAS
MATLAB VOLTAJE, CONVERGENCIA CORRIENTE CORRIENTE Y REDES
POTENCIA Y DE RESULTADOS ANTES DE LA TEMPERATURA ELECTRICAS
CORRIENTE EXPANSION
TERMICA
Figura 33. Indicador de solución - Estado del arte
31

Ups - TTS208 Analisis de Cargas de Redes Electricas

Cargado por

Copyright:

Formatos disponibles

Ups - TTS208 Analisis de Cargas de Redes Electricas

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ups - TTS208 Analisis de Cargas de Redes Electricas

Cargado por

Copyright:

Formatos disponibles

“ANÁLISIS DE CARGABILIDAD DE REDES ELÉCTRICAS DE DISTRIBUCIÓN BASADO EN LA

MODIFICACIÓN DEL ALGORITMO DE GAUSS-SEIDEL”

Trabajo de titulación previo a la obtención del título de

Quito, noviembre 2020

“ANÁLISIS DE CARGABILIDAD DE REDES ELÉCTRICAS DE DISTRIBUCIÓN

Universidad Politécnica Salesiana

Breve reseña historia e información de contacto:

Todos los derechos reservados:

Quito, noviembre 2020

Yo, Henry Saúl Gancino Bustillos, con documento de identificación N° 1725114076,

El presente trabajo de investigación se lo dedico a Dios quien nunca dejo de bendecirme

Agradezco a la Empresa Eléctrica Quito, por darme la oportunidad de realizar mis

2. Gauss-Seidel Para Flujos de Potencia. ................................................................... 4

2.1. Flujos de Potencia en Redes Radiales ................................................................ 5

2.2. Parámetros de Cargabilidad. .............................................................................. 6

2.3. Expansión térmica en los conductores ............................................................... 8

2.4. Comportamiento Térmico en los conductores. .................................................. 9

2.5. Pandeo en las Líneas ........................................................................................ 11

3. Estimación de Conductores en Distribución ........................................................ 12

3.1. Selección de los Conductores Eléctricos.......................................................... 12

3.2. Determinación de la Corriente por un Flujo de Carga. .................................... 13

3.3. Efecto del Exceso de Corriente y Temperatura en los Conductores ................ 14

4. Metodología Para Determinar Flujos en los Sistemas de 13 y 34 nodos. ............ 14

4.1. Algoritmo Para la Aplicación de Gauss-Seidel en Distribución. ..................... 15

4.2. Verificación de los Resultados por DIgSILENT Power Factory. .................... 16

4.3. Calibre de Conductores Seleccionados de Acuerdo al Nivel de Corriente. ..... 17

4.4. Límite de Cargabilidad de los Alimentadores de 13 y 34 Nodos. ................... 18

4.5. Expansión Térmica del Sistema Alimentador de 13 y 34 Nodos de la IEEE. . 19

5.1. Flujo de Potencia. ............................................................................................. 20

5.2. Flujos en DIgSILENT de los Sistemas de 13 y 34 Nodos. .............................. 20

7. Trabajos futuros ................................................................................................... 23

8.1. Matriz de Estado de Arte.................................................................................. 27

8.2. Resumen de Indicadores .................................................................................. 30

Figura 1. Etapas de planificación ........................................................................................................ 4

Figura 2. Temática de la investigación.............................................................................................. 4

Figura 3. Modelo de impedancias en serie y en derivación. ....................................................... 4

Figura 4. Topología radial. ................................................................................................................... 6

Figura 5. Perdidas de una red subterránea de distribución. ......................................................... 7

Figura 6. Perdidas de una red aérea de distribución. ..................................................................... 7

Figura 7. Límite del nivel de sobrecarga SIL. ................................................................................. 7

Figura 8. Curva características de las líneas. ................................................................................... 8

Figura 9. Estado de la curva a diferentes temperaturas ................................................................ 9

Figura 11. Sistema de alimentación de 13 nodos de la IEEE ................................................... 15

Figura 12. Sistema de alimentación de 34 nodos de la IEEE. .................................................. 15

Figura 15. Corriente proporcional a temperatura en el sistema de 13 nodos de la IEEE. 18

Figura 16. Corriente proporcional a temperatura en el sistema de 34 nodos de la IEEE. 19

Figura 17. Punto crítico promedio en el sistema de 13 nodos de la IEEE. ........................... 19

Figura 18. Punto crítico promedio en el sistema de 34 nodos de la IEEE. ........................... 20

Figura 21. Flujo en DIgSILENT del alimentador de 13 nodos. ............................................. 20

Figura 22. Flujo en DIgSILENT del alimentador de 34 nodos. ............................................... 21

Figura 23. Punto crítico de la configuración 601. ........................................................................ 21

Figura 25. Punto crítico de la configuración 603. ........................................................................ 21

Figura 26. Punto crítico de la configuración 605. ........................................................................ 21

Figura 27. Punto crítico de la configuración 300. ........................................................................ 22

Figura 28. Punto crítico de la configuración 301. ........................................................................ 22

Figura 29. Punto crítico de la configuración 302. ........................................................................ 22

Figura 30. Punto crítico de la configuración 303. ........................................................................ 22

Figura 31. Resumen e indicador de la temática - Estado del arte............................................ 30

Figura 33. Indicador de solución - Estado del arte ...................................................................... 31