1

Estudio de Flujo de Carga (26 Noviembre 2019)

María Fernanda Valle Florez

Abstract- In electrical engineering, the study of load flow is diagrama de una sola línea, define la configuración del
an important tool that involves numerical analysis applied to a sistema y la ubicación y el tamaño de las cargas, la
power system. In the study of the flow of load usually a generación y el equipo. Está organizado en una lista de
simplified notation is used like the diagram unifiliar and the datos que define el modelo matemático para el componente
system by unit, and concentrates in several forms of the
electrical power AC (for example, voltages, angles of the
del sistema de alimentación de caché y cómo los
voltages, active power and reactive power ). This study analyzes componentes están conectados entre sí un ejemplo de esta
the power systems operating in the steady state. There are configuración se da en la (figura 1), a esta configuración
several programs that implement the study of power flow. esquemática a es a la que se denomina diagrama unifilar,
Index Terms—Flujo de Carga, Newton Raphson, Gauss- donde mediante este podemos representar el complejo
Seidel. sistema de potencia de una manera reducida y sencilla de
I. INTRODUCCIÓN

U no de los procedimientos más usados en el análisis

del sistema de potencia es el cálculo del flujo de
carga. La planificación, el diseño y la operación
de los sistemas de energía requieren cálculos para analizar
el rendimiento en estado estable del sistema de potencia
bajo varias condiciones de operación y para estudiar los
efectos que se generan por los cambios en la
configuración del equipo.
La mayoría de los estudios de carga baja para análisis de
sistemas se basan en modelos de red estáticos. En general, analizar.
los servicios públicos utilizan sistemas automáticos de
control de supervisión y adquisición de datos (SCADA).
Dichos sistemas se utilizan principalmente para operar Fig. 1 Diagrama de una línea típica de un sistema eléctrico industrial
para el estudio de carga.
herramientas para la optimización de la generación, el
control de línea de enlace, la distribución, las pérdidas y el
control de línea de enlace. Esta discusión se ocupa
I. DATOS DE ENTRADA
únicamente de los modelos de red estática y su análisis.
Los datos se dividirán en categorías de datos: datos del
Debido a que el problema del flujo de carga se refiere a sistema, datos de bus, datos de generador, datos de línea y
la operación equilibrada y en estado estable de los sistemas datos de formularios de respuesta.
de energía, se utiliza un modelo de secuencia positiva y  Datos del sistema
monofásica del sistema de energía. El flujo de carga es
también básico para diferentes tipos de de estudios como Por las facilidades que da trabajar con el sistema en
los son el estudio de corto circuito, estabilidad, arranque de p.u. para el desarrollo del sistema, se trabajara con esta
motor y de armónicos. herramienta en vez de utilizar los datos en voltios,
amperios y límites. El sistema por unidad permite
II. OBJETIVOS algoritmos de solución más fáciles y eficientes y también
simplifica el análisis del flujo de carga. La conversión de
A. Conocer el estudio de flujo de carga y sus los datos del sistema a una representación por unidad
diferentes aplicaciones. requiere la selección de una base kVA y un voltaje base.
B. El objetivo de un programa de flujo de potencia es La selección de la base kVA y el voltaje de la base
de analizar el desempeño de una red eléctrica en régimen especifican el impedimento de la base y la tensión.
permanente bajo diferentes condiciones de explotación. Se  Datos del bus
trata de la herramienta de análisis de base para el Los datos del bus describen el bus central y la carga y
planeamiento, el diseño y la explotación de cualquier red las derivaciones conectadas a ese bus. Los datos incluyen:
de energía eléctrica, ya sea de transporte o de distribución
1. Bus número.
de energía, pública o industrial.
2. Bus nombre.
C. El objetivo de este estudio es realizar una 3. Tipo de bus.
evaluación general del sistema donde se instalará la
compensación de reactivo. 4. Carga.
5. Shunt.
III. REPRESENTACION DEL SISTEMA
6. Por unidad de voltaje y ángulo.
La información del sistema, que se muestra en el
 2

7. Bus base kV A continuación se mostrara un ejemplo de la conexión

El número de bus es normalmente el índice principal de y se mostrara gráficamente lo anterior mencionado sobre
la información sobre el bus. Por ejemplo, se utiliza para la conexión de un generador a un bus y sus
definir las conexiones de línea en los datos de línea y se especificaciones.
usará para obtener la salida de un bus durante la ejecución
del programa, normalmente el nombre es solo para usos
informativos. El tipo de bus es un código de entrada del
programa que permite al programa organizar
adecuadamente los buses para la solución de flujo de  Datos de línea
carga. 'Típicamente, los cuatro tipos son:
Los datos también se ingresan para la línea de caché en
el sistema. Aquí, el término "línea" se refiere a todos los
(1) Autobuses de carga elementos que conectan dos buses, incluidas las líneas de
(2) Autobuses de generador transmisión, los cables, los reactores en serie, los
condensadores en serie y los transformadores. Los
(3) Autobuses de giro
elementos de datos incluyen:
(4) Autobuses desconectados
1. Resistencia.
Los términos bus de "carga" y bus de "generador" no
2. Reactancia.
deben ser escritos demasiado literalmente. Un bus de
"carga" es cualquier bus que no tiene un generador. Un bus 3. carga de suceptancia (capacitancia shunt)
de "carga" no necesita tener carga; simplemente puede ser 4. clasificaciones de línea.
un punto de interconexión para dos o más líneas. Un bus 5. Línea en servicio o fuera de servicio en código.
"generador" también podría tener carga conectada. El bus
6. Línea conectada en Shunt.
de "oscilación" o "holgura" es un tipo especial de bus de
generador que necesita el proceso de solución.
Como se describe en 4.7, las líneas están representadas por
un modelo 7 con resistencia en serie y reactancia y la mitad
La carga de bus normalmente se ingresa en MW y Mvar
de la susceptibilidad de carga colocada en cada extremo de
a voltaje nominal. Normalmente, la carga se trata como un
la línea. La resistencia, la reactancia y la susceptancia
MVA constante, es decir, independiente del voltaje. En
generalmente se ingresan en unidades por unidad o
algunos casos, también se ingresará un componente
porcentaje, según la convención del programa.
constante de constancia o impedancia constante de la
carga, de modo que la carga sea una función de la tensión.
Se debe tener cuidado para garantizar que se utilice la  Datos del transformador
convención de señalización adecuada para distinguir los Se requieren datos adicionales para los
reactores de los condensadores según se define para el transformadores. Esto se puede ingresar como parte de los
programa de carga particular que se está utilizando. Los datos de línea o como una categoría de datos separada
datos del bus generalmente contienen una estimación dependiendo del programa de flujo de carga particular que
inicial del voltaje que se utilizará como punto de partida se esté utilizando. Este dato adicional usualmente incluye:
para la solución de flujo de carga.
1. Ajuste de tap en por unidad.
2. Angulo del Tap en grados
 Datos del generador
3. Máxima posición del tap
Los datos de Generador se ingresan para cada generador
4. Mínima posición del tap
en el sistema, incluido el generador de giro de sistema. Los
datos definen la salida de potencia del generador y cómo el 5. Rango de voltaje programado con tamaño de paso de
generador controla la tensión. Los elementos de datos toma o y voltaje programado fijo utilizando una
normalmente ingresados son: aproximación de toma continua.
1. Potencia Real de Salida en MW. Los últimos tres elementos de datos son necesarios solo
para los transformadores de cambio de toma de carga
2. Máxima potencia reactiva de salida en MVAR (es
(LTC) que varían automáticamente su configuración de
decir máximo limite reactivo de la maquina)
toma para controlar el voltaje en un lado del
3. Mínima potencia reactiva de Salida en MVAR (es transformador.
decir MINIMO limite reactivo de la maquina).
En la fig. 3 se da un breve ejemplo de cómo es la
4. Voltaje programado en por unidad. conexión de los elementos anteriormente mencionados y
5. Generador en servicio o fuera de servicio en código. además se dan sus respectivas especificaciones, aparte de
Otros elementos de datos que pueden incluirse son la esto podemos ver la representación de un sistema de
base de MVA del generador y la impedancia interna del potencia mediante el diagrama unifilar.
generador para uso en estudios de cortocircuito y
dinámicos. El uso de los elementos de datos anteriores para
determinar el voltaje del generador y la potencia reactiva.
 3

técnica de Gauss-Seidel y la técnica de Newton-Raphson.

Ambos resuelven las ecuaciones de bus en forma de
admitancia. Este sistema de ecuaciones es muy usado
V. METODOS DE SOLUCION DE FLUJO DE CARGA debido a la preparación de los datos y la facilidad con que
El desarrollo de hardware fue de gran necesidad para se puede formar y modificar la matriz de admitancia del
poder solucionar de una manera más eficiente, óptima y bus en los casos posteriores.
exacta las ecuaciones de red, ya que este procedimiento a En un estudio de flujo de carga, los parámetros principales
mano en un sistema de gran escala era casi imposible de son:
solucionar, además de ser un desarrollo muy lento. Ward y P = Potencia activa en la red de transmisión
Hale desarrollaron una carga exitosa (programa bajo
Q = Potencia reactiva en la red de transmisión
usando un procedimiento iterativo de Newton modificado
en 1956). La aplicación del algoritmo de iteración de | V | = Magnitud del voltaje del bus (voltaje a tierra)
Gauss-Seidel siguió poco después. La investigación en δ = Ángulo del voltaje del bus referido a una referencia
algoritmos continuó y el método de Newton-Raphson se común
introdujo en carly 1960's. En los años intermedios se ha La Tabla 1 resume las especificaciones eléctricas
realizado una investigación considerable para mejorar el estándar para los tres tipos de bus. Las clasificaciones "bus
rendimiento de estos algoritmos, haciéndolos más generador" y "bus de carga" no deben tomarse como
robustos, capaces de manejar componentes adicionales del absolutas. Habrá, por ejemplo, ocasiones en las que P y | V
sistema de energía y permitiendo cambios en la red mucho puedan especificar un bus de carga pura.
más grandes. TABLA. 1
ESPECIFICACIONES DEL BUS DE FLUJO DE CARGA
A. Formulación del problema.
El cálculo del flujo de carga es un problema de
solución de red. Los voltajes y corrientes están
relacionados por la siguiente ecuación:
[ I ] =[ Y ]∗[ V ] ( 2 )
Donde:
I: Vector de corrientes totales en secuencia positiva
en los nodos de la red (buses)
V: Vector de voltajes totales en secuencia positiva en
los nodos de la red (buses)
Y: Matriz de Admitancias de la red.
Se debe tener en cuenta que esta ecuación es una 1) Técnica iterativa de Gauss-Seidel
ecuación algebraica lineal con coeficientes Está basado en los métodos de resolución iterativos de
complejos. ecuaciones lineales de los matemáticos alemanes Carl
Friedich Gauss y Phillip Ludwig von Seidel. Es un método
Las condiciones de los terminales en cada bus se bastante sencillo de emplear y su metodología es la más
describen normalmente en términos de potencia simple de todos, sin embargo, requiere de un gran número
activa y reactiva (P y Q). La corriente de bus de iteraciones para converger a una solución.
podemos escribirlas de la siguiente manera para Para realizar el flujo de carga por este método debe
facilitar el cálculo de la misma: llevarse el sistema a por unidad y deben seguirse los
¿ siguientes pasos:
( Pi + Q i )
I i= ¿ (3) 1) Definir tipos de barras según la clasificación anterior
Vi y obtener todos los parámetros de la red (datos de
Combinando las ecuaciones (2) y (3) tenemos: impedancias).
¿ 2) Se asumen todos los voltajes iniciales de las barras
( Pi+ Q i )
¿ =[ Y ]∗[ V ] (4) iguales a 1p.u y con ángulo cero, excepto aquellas barras
Vi que sean PV donde se conoce |V| y en la Slack donde se
La ecuación 4 no es lineal y no puede resolverse conoce módulo y ángulo.
fácilmente mediante técnicas de matriz de forma cerrada. 3) Comienza a iterarse actualizando V para todas las
Debido a esto, las soluciones de flujo de carga se obtienen barras PQ, utilizando la relación:
mediante procedimientos que involucran técnicas
iterativas.
B. Algoritmos de solución iterativa
Muchas de estas son variaciones de una u otra de las
dos técnicas básicas que la industria utiliza hoy en día: la Donde 𝑃𝑁𝑒𝑡𝑎 = 𝑃𝑔𝑒𝑛 − 𝑃𝑑𝑒𝑚 y 𝑄𝑁𝑒𝑡𝑎 = 𝑄𝑔𝑒𝑛 −
 4

𝑄𝑑em.
4) Cuando se trate de una barra PV para la cual no se
conoce Q, a esta hay que estimarla siguiendo la siguiente
ecuación.

Esta Q representa solamente una cifra aproximada para

los valores actuales de V. Una manera simple de ver la
ecuación anterior es que todas las admitancias de las
barras con las que existe conectividad a la barra PV
evaluada, serán multiplicadas por el voltaje en cada una de
esas barras según el valor más actualizado de ese voltaje. Del lado izquierdo de la ecuación se tiene el vector de
5) Una vez calculada Q aprox. se procede a calcular el correcciones y del lado derecho el vector de errores. La
voltaje (en modo fasorial) en esa barra PV, siguiendo la matriz J es referida como la matriz Jacobiana. Para
ecuación: empezar debe calcularse el vector de errores del sistema,
esto viene dado por la diferencia entre la potencia
calculada y el valor de potencia especificada, esto para la
cantidad de barras PQ y PV, por tanto puede decirse que
cada barra PQ, se tendrá un desajuste ΔP y un ΔQ
correspondiente a dicha barra, y para cada barra PV se
En el caso de una barra PV, la Qgen se asume como la
tendrá un ΔP correspondiente.
Q aproximada calculada anteriormente y Qdem dependerá
de si existe una carga en esa barra que demande Q. Hay Por otro lado, el vector de correcciones para voltajes y
que aclarar, que si bien las barras PV no suelen tener ángulos de fase se arma igualmente según el número de
cargas en su propia barra, es posible que eso ocurra. Es de barras PV y PQ. Para las barras PQ deben determinarse las
notar, que en el caso en que no exista demanda, la 𝑄𝑁𝑒𝑡𝑎 variables Δδ y Δ |V|, mientras que para las PV sólo debe
siempre será negativa debido al signo de Q aprox. determinarse Δδ. El orden de ambos vectores es como se
describe en la imagen.
6) Se repite el proceso para todas las barras PV tomando
siempre en cuenta que los valores de los voltajes en las 3) Se calculan los ΔP y ΔQ del vector de errores. Estos
ecuaciones deben ser los más recientes obtenidos vienen dados por 𝛥𝑃 = 𝑃𝑁𝑒𝑡𝑎 − 𝑃𝑐𝑎𝑙𝑐 y
7) Las barras PQ tienen P y Q conocidas, por lo que no ΔQ=𝑄𝑁𝑒𝑡𝑎 − 𝑄 𝑐𝑎𝑙𝑐 , donde las potencias netas vienen
es necesario calcular la Q aproximada, y dado que siempre dadas por 𝑃𝑁𝑒𝑡𝑎 = 𝑃𝑔𝑒𝑛 − 𝑃𝑑𝑒𝑚 y 𝑄𝑁𝑒𝑡𝑎 = 𝑄𝑔𝑒𝑛 − 𝑄𝑑𝑒𝑚
son barras modeladas como demanda, y si hubiere y donde por otro lado, las potencias calculadas serán:
generadores, su potencia activa y reactiva generada son
fijas. Tomando esto en cuenta se aplica la misma ecuación
de voltaje (4), en las cuales PNeta = Pgen − Pdem y QNeta
= Qgen − Qdem
8) Una vez que todos los voltajes y ángulos en todas las
barras hayan sido calculados, se habrá completado una
iteración. Para culminar el proceso debe cumplirse que
𝑉𝑖𝑘 − 𝑉𝑖𝑘−1 ≤ 𝜀. Donde ε será estipulado según Las ecuaciones anteriores provienen de las ecuaciones
convenga, tanto para la magnitud del voltaje como para su de flujo de potencia por las líneas. Las sumatorias
ángulo de fase (típicamente 0.001). anteriores pueden entenderse como la sumatoria de los
efectos de todas las líneas que confluyen en la barra i.
4) Se crea la matriz Jacobiana. Esta puede ser
2) Técnica iterativa de Newton Raphson
dividida en cuatro submatrices.
Este es el método más conocido y utilizado por su gran
precisión y convergencia para un gran número de
escenarios en pocas iteraciones.
Igualmente, es conveniente tener el sistema en por
unidad para realizar el flujo de carga. Los pasos para
realizar este método son:
1) Definir clasificación de barras del sistema, tener los
parámetros (impedancias, frecuencia) de la red y asignar
un número a cada barra.
2) Se define la ecuación de Newton-Raphson para SEP’s
con las variables a calcular, excluyendo para el estudio la
barra Slack En primer lugar, tenemos la submatriz J11, la cual
 5

indica la variación de potencia activa respecto al ángulo en calculadas en la barra i.

cada barra. Se varían las potencias (de barras PQ y PV) 6) Una vez obtenidos todos los valores del Jacobiano, se
respecto a ángulos (de barras PQ y PV). ejecuta la ecuación de Newton Raphson (2) y se obtiene el
vector de correcciones con Δδ y ΔV/|V|. Una vez obtenido,
se calculan los nuevos δ y V :

Con estos nuevos valores se tienen los nuevos

Es importante destacar que en el caso de J12 se varían voltajes con módulos y ángulos, concluyendo así una
iteración.
las potencias (de barras PQ y PV) respecto a los voltajes
desconocidos de las barras PQ, es decir, los voltajes en la 7) Se repetirá el proceso con los nuevos valores de
derivada parcial corresponden a los voltajes de las barras voltaje hasta que se llegue a la convergencia deseada.
PQ.
VI. ANALISIS DE FLUJO DE CARGA
Se realizan estudios de flujo alto para verificar el
funcionamiento de un sistema existente en condiciones
normales o de interrupción, para ver si el sistema existente
es capaz de suministrando cargas adicionales planificadas,
o para verificar y comparar nuevas alternativas para
Para el caso de la submatriz J21, se varía la potencia adiciones al sistema para suministrar nueva carga o
reactiva de las barras PQ respecto a los ángulos de las mejorar el rendimiento del sistema.
barras PV y PQ. Generalmente, el ingeniero de estudio tiene un
conjunto predefinido de criterios que el sistema debe
cumplir. Éstos incluyen:
A. Criterios de voltaje
Los criterios de voltaje generalmente se dividen en un
rango de voltaje aceptable para condiciones normales y un
rango más amplio de voltaje aceptable en condiciones de
corte.
Para el caso de la submatriz J22 la variación viene dada
B. Los flujos en las líneas y los transformadores deben
entre la potencia reactiva de las barras PQ y los voltajes de
estar dentro de los valores térmicos definidos.
dichas barras.
Estos criterios térmicos para líneas y transformadores
5) Se calculan los elementos de la matriz Jacobiana
también pueden tener tal división, lo que permite una
utilizando las siguientes aproximaciones de las derivadas.
capacidad de sobrecarga temporal debido a la constante de
Es conveniente dividir los elementos diagonales de las
tiempo térmica del equipo o la capacidad adicional de
matrices de los no diagonales, en este sentido tenemos las
enfriamiento forzado de los transformadores.
siguientes consideraciones.
C. Las salidas reactivas del generador deben estar
dentro de los límites definidos por el generador
Elementos Diagonales:
Un estudio normalmente comienza con la preparación
de casos básicos para representar los diferentes modos de
operación del sistema o planta. La condición de operación
normalmente elegida es la carga máxima. Cuando la carga
máxima se produce en diferentes momentos en diferentes
partes del sistema, pueden ser necesarios varios casos base.
Los casos base deben representar condiciones de operación
realistas.
Elementos No-Diagonales: Los casos base se analizan para determinar si los
voltajes y los retardos están dentro de los rangos
aceptables. Si se observan problemas de voltaje o
sobrecarga, se pueden hacer cambios en el sistema a los
datos de flujo de carga y el caso se resuelve en caso de que
los cambios sean efectivos para solucionar el problema.
Para remediar los problemas de baja tensión, los cambios
posibles incluyen:
Donde Pik y Qik se refieren a los flujos de potencia 1. Cambio en las posiciones de tomas del
entre las barras i y k y Pi y Qi son las potencias transformador.
 6

2. Incremento en el voltaje del horario del generador Generador No. Unidades S (MVA) Pbase S (pu)
3. Adición de condensadores de derivación Salto Grande 7 150 100 10.5
4. Reconfiguración del sistema para cambiar la carga a Baygorria 3 36 100 1.08
líneas menos cargadas Palmar 3 111 100 3.33
5. Desconexión de reactores en derivación. Terra 2 40 100 0.8
6. Adición de líneas o transformadores. Montevideo E 3 4º 62.5 100 3.1875
Para remediar las cargas de líneas pesadas o de 5º 100
transformadores, se aplican la mayoría de los mismos 6º 156.25
remedios. En general, los dos primeros de los remedios Montevideo L 2 142.75 100
anteriores no ayudarán a cargas pesadas debido a grandes condiciones del caso base y se hacen adiciones, si es
flujos de potencia real (vatios). El flujo de potencia real necesario, para que el sistema cumpla con los criterios de
está determinado por la separación de fase del bus de rendimiento. Luego se estudian las condiciones de
suministro que conduce al ángulo de fase del bus de carga, interrupción y nuevamente se requieren cambios en el
con magnitudes de voltaje que tienen un efecto secundario. sistema. Los estudios de los sistemas futuros varían en
Sin embargo, el flujo de potencia reactiva se determina cuanto a que usualmente hay más formas alternativas de
principalmente por la magnitud del voltaje con la potencia resolver los problemas encontrados. El flujo alto es la
reactiva que va desde el bus de mayor voltaje al bus de herramienta que permite comparar las alternativas en
menor voltaje. El flujo de energía real y reactiva, que está términos de su efectividad en condiciones normales y de
influenciado principalmente por diferentes restricciones. contingencia. Junto con otros estudios, así como datos de
costo y confiabilidad, los resultados conducen a la
La toma nominal del transformador puede cambiar la
selección de la mejor alternativa.
relación relativa de la tensión en el bus primario y
secundario y, por lo tanto, pueden cambiar el flujo de VII. EJEMPLO DE ESTUDIO DE FLUJO DE CARGA CON
potencia reactiva, mientras que el flujo de potencia real no SOFTWARE
se ve afectado en gran medida por un cambio en la
El ejemplo planteado consiste en una red simplificada del
posición de toma.
sistema eléctrico uruguayo. Se detalla el esquema unifilar en
Cuando los voltajes y flujos de la caja base están en el la figura siguiente:
rango deseado, el sistema debe examinarse para verificar el
funcionamiento en condiciones anormales (análisis de
contingencia). Estas condiciones incluyen:
* Pérdida de una línea de transmisión o cable.
* Pérdida de un transformador.
* Pérdida de un generador. Condiciones anormales de
suministro.
El análisis del flujo de carga se usa para diseñar un
sistema que tenga un buen perfil de voltaje y cargas de
línea aceptables durante el funcionamiento normal y que
continuará funcionando de manera aceptable cuando una o
más líneas se vuelvan inoperantes debido al daño de la
línea, los rayos, el fallo de los transformadores, etc.
Además de los beneficios descritos anteriormente, un
estudio de la potencia reactiva puede ayudar para hallar el
mínimo en las ramas, puede llevar a pérdidas de línea
reducidas y una mejor distribución de voltaje. La
reducción en la demanda de KVA debido a la corrección
del factor de potencia puede llevar a menores facturas de
servicios públicos para una planta industrial. El tamaño y
la ubicación de los condensadores de corrección del factor
de potencia y la configuración de los voltajes programados Fig. 4 Representación del sistema de potencia del ejercicio propuesto
del generador y las posiciones de la toma del
Para el modelado de los distintos elementos de la red se
transformador se pueden estudiar con los flujos de carga.
tuvieron las siguientes consideraciones:
Los estudios de carga ahora también pueden
proporcionar datos para el control automático de la carga y
la demanda, si es necesario, además de esto también se Generadores: Para el tratamiento de los generadores de
utiliza para verificar los efectos del crecimiento futuro de energía eléctrica se tomaron de forma de tener un
la carga y la efectividad de las adiciones planificadas. generador equivalente con la suma de las potencias
Estos estudios se realizan de la misma manera que los aparentes de todas las máquinas en cada usina generadora.
estudios del sistema actual. Las cargas futuras se Para el modelo equivalente del generador se tomó el
determinan y se introducen en el nodal. Se estudian las mismo como una inyección de potencia activa fija P, y de
 7

reactiva Q variable dentro de un rango acotado [Qmin,

Qmax] controlada por un valor de consigna de tensión en
las barras a las que se conecta el generador.
Cargas: Para las cargas a considerar se decidió
tomar el peor caso. Para ello se ingresaron las
potencias activas y reactivas del pico de invierno
pasado. Esta información fue brindada por el
Despacho Nacional de Cargas. Se tomó potencia base
100 MVA.

Red de 150 kV

Líneas: Al tomar un sistema eléctrico simplificado, en

algunos casos se tomaron los diferentes tramos de líneas,
de forma de unirlos en una línea sola entre barras.

Transformadores: Se tomaron los siguientes datos

brindados por el Despacho Nacional de Cargas.

Red de 500 kV
 8

Se realizó una primera corrida en la cual no se colocaron

reguladores obteniéndose los siguientes resultados:

A) Entrada de datos:
 9

B) Resultados:
Registrándose los siguientes resultados:
 10

C) Modificaciones
Para elevar la tensión de las barras Nueva Palmira y
Conchillas se decidió colocar un regulador de tensión en
San Javier. Se obtuvieron los siguientes resultados:

1) Archivo de entrada de datos:

 11

2) Resultado:
 12

Observación: Nótese en los resultados la elevación de la

tensión en las barras de Nueva Palmira y Conchillas.

IX. EJERCICIO DE NEWTON RAPHSON

Dado un sistema de 4 nodos y 4 líneas, como se muestra en
la Fig 5, solucionarlo por el método de Newton Raphson.
 13

Fig.5. Representación del sistema de potencia de ejemplo.

Datos de líneas:
Resolviendo para Δδ y ΔV, y recalculando δ y V
tenemos:

Datos de barras: 3) Tercera iteración

 Errores de potencia
[ΔP|ΔQ] = [−0.0339 −0.1455 −0.0424 −0.0414
−0.1646]*10−3
 Cumple tolerancia: [SI]

Resultados:

Resolviendo para Δδ y ΔV, y recalculando δ y V

tenemos:

Inicio del subproblema 1

Condiciones iniciales:
1. Voltajes incógnitas = 1.0 p.u.
2. Ángulos incógnitas = 0º Resultados:
1) Primera iteración Para el cálculo se aplica una potencia base 100
 Errores de potencia MVA
[ΔP|ΔQ] = [−1.5966 −1.9395 2.2129 −0.4465 −0.8345]  P1 = 136.7948 Mw
 Cumple tolerancia: [NO]  P2 = −169.9966 Mw
 P3 = −199.9854 Mw
 P4 = 237.9958 Mw
Para el cálculo se aplica una potencia base 100
MVA
 Q1 = 83.4977 MVAR
Resolviendo para Δδ y ΔV, y recalculando δ y V
tenemos:  Q2 = −105.3459 MVAR
 Q3 = −123.9235 MVAR
 Q4 = 131.8393 MVAR
Flujos de Potencia:
 P12 = 38.6883 Mw Q12 = 22.2969
Mvar
2) Segunda iteración  P13 = 98.1065 Mw Q13 = 61.2008
 Errores de potencia Mvar
[ΔP|ΔQ] = [−0.0323 −0.0645 0.0359 −0.0342 −0.0620]  P24 = -131.5350 Mw Q24 = -74.1109
 Cumple tolerancia: [NO] Mvar
 P34 = -102.9101 Mw Q34 = -60.3650
Mvar
 P21 = -38.4616 Mw Q21 = -31.2349
Mvar
 P31 = -97.0753 Mw Q31 = -63.5585
Mvar
 14

 P42 = 133.2504 Mw Q42 = 74.9166

Mvar
 P43 = 104.7453 Mw Q43 = 56.9227
Mvar

VIII. CONCLUSIONES
1) Debería ser evidente una herramienta que predice el
rendimiento de sus sistemas eléctricos en diversas
condiciones de operación que realmente pueden
encontrarse antes de que ocurran estas condiciones es
de valor, una herramienta de este tipo se puede usar
para evitar costosos daños, equipos dañados y
posiblemente incluso la pérdida de vidas. La solución
de flujo de carga proporciona un medio para estudiar
sistemas en condiciones reales o hipotéticas. Los
resultados de la solución deben evaluarse y analizarse
con respecto a la operación presente y futura óptima.
Esto conduce a un diagnóstico del sistema tal como
existe. El análisis que muestra el camino para mejorar
la operación proporciona una base significativa para la
planificación futura del sistema.
2) En adición al análisis de flujo de potencia, llamado en
ocasiones el caso base, muchas implementaciones
software ejecutan otros tipos de análisis, tal como el
análisis de falla de cortocircuito, estudios de
estabilidad (transitorio y estado estable), unit
commitment y despacho económico. Particularmente,
algunos programas usan programación lineal para
encontrar el flujo de carga óptimo, es decir, las
condiciones que dan el costo más bajo por kilovatio-
hora entregado, además de las grandes utilidades ya
mencionadas en este artículo podemos contar con
estas otras implementaciones para las cuales nos
pueden servir el flujo de carga, lo cual lo hace un gran
método para la solución de los sistemas de potencia.
3) Una de las importancias más grandes que tiene el
análisis de flujo de carga se puede dar en su utilización
para predecir el funcionamiento de nuestro sistema
para la conexión de nuevas cargas, ya que gracias a
esto podemos diseñar el sistema para un futuro y no
solo para el presente.

REFERENCIAS

[1] IEEE Std 141-1986,Recomended Practice for Electric Power

Distribution for Industrial Plants (ANSI).
[2] BROWN,H,E,Solution of Large Networks by Matrix
Methods;New York;NY:John Wiley and Sons,1976
[3] “Análisis de Sistemas de Potencia”, notas de clase de
API1,Departamento de Ingeniería Eléctrica, Escuela Colombiana
de Ingeniería Julio Garavito,2019
[4] Alfredo C.,Claudio O. (2002, Marzo 10). Resolución del
problema del Flujo de Cargas [Online]. Available:
https://iie.fing.edu.uy/~jesus/proy_flucar.pdf
[5] Carlos H. (2016, Dec 10). Métodos de Flujo de Carga [Online].
Available:https://www.researchgate.net/publication/311796467_
Metodos_de_Flujo_de_Carga
[6] Wikipedia (2019, Feb 15). Flujo de Potencia [Online]. Available:
https://es.wikipedia.org/wiki/Flujo_de_potenci
15