FORMULACION DEL PROBLEMA DE FLUJO DE POTENCIA

En este capítulo se tratarán los métodos clásicos de flujo de carga, pensados

inicialmente para sistemas de potencia pero que pueden ser extendidos a
sistemas de distribución.

La solución del problema de flujos de potencia es básica para la mayoría de los

análisis que se realizan en sistemas eléctricos de potencia. Esto ha tenido como
consecuencia que el esfuerzo que se ha dedicado al desarrollo de métodos de
solución sea notable. Los métodos de flujo de carga se aplican en la planeación,
operación y control de sistemas eléctricos de potencia y distribución, donde los
flujos de potencia se resuelven para diferentes casos, por ejemplo, en la
evaluación de la seguridad, donde se requiere determinar la operación ante
contingencias, o como elemento auxiliar en estudios de reconfiguración de redes
de distribución, localización de capacitores en las mismas, o para evaluar las
condiciones iniciales en estudios de fallas, entre otros.
Si se tiene un sistema eléctrico de distribución en el cual se conocen las cargas en
las barras, los parámetros de la red, y se conocen dos de sus cuatro variables
características (V , θ , P y Q ), se puede determinar el estado de operación del
sistema.

Los diferentes elementos que conforman una red de distribución están diseñados
para operar bajo condiciones de equilibrio, y esta se logra de manera balanceada
debido a que la carga en cada una de las barras es igualmente balanceada. Es
esta 25 la razón por la cual el flujo de carga a través de los elementos puede ser
calculado mediante un equivalente monofásico. Como la red de distribución es
trifásica, se debe descomponer en tres redes de secuencia, denominadas de
secuencia positiva, negativa y cero. Pero solo la red de secuencia positiva
presenta efecto sobre el flujo de carga en condiciones balanceadas.

Los estudios de flujo de carga son utilizados en la planeación del sistema eléctrico,
para determinar las condiciones operativas optimas. El flujo de carga proporciona
información como son los voltajes nodales, flujos de potencia por las líneas,
perdidas entre otros.
Básicamente, el problema de flujos de potencia convencional puede definirse
como el cálculo de voltajes nodales y, posteriormente, el de flujos de potencia a
través de cada elemento de la red, para valores conocidos de generación y carga
nodales en MW y MVAr, en un instante de tiempo específico.
La solución del problema puede o no estar sujeta a restricciones de red, tales
como límites de generación de potencia activa y reactiva, magnitud de voltajes
nodales, así como flujos en elementos, entre otras.
Los métodos clásicos mas empleados para el análisis de flujo de carga son
desarrollados en base a la matriz de red Ybus

Los elementos que componen la matriz Ybus se pueden describir de la siguiente

forma:
El modelo matemático que se utiliza para solucionar el problema de flujo de carga
es expresado en función de los voltajes en las barras, estos voltajes se pueden
expresar 26 en componentes polares o rectangulares, pero usualmente se
expresan en forma polar.

La ecuación mediante la cual se relacionan las corrientes inyectadas en las barras

en función de los voltajes de las barras, y la matriz Ybus es la siguiente:

Si se tiene que Pi y Qi son las potencias activa y reactiva especificadas en la barra

i . Entonces la ecuación que relaciona estas variables con la tensión en la barra i ,
y la matriz Ybus es:

En forma polar se expresa de la siguiente forma:

El método de flujo de carga debe hallar el conjunto de soluciones (voltajes), que

hacen que la potencia calculada sea aproximadamente igual a la potencia
especificada en cada barra. El error de potencia en la barra i se calcula como la
diferencia entre la potencia inyectada por el sistema externo e interno, esta
diferencia se presenta en el flujo de carga en el cual la potencia del sistema
externo es fijada en un valor especifico, y la del sistema interno se determina
basándose en los voltajes en las barras y en los parámetros de la red, solo hasta
que se alcanza el voltaje de equilibrio, las dos 27 potencias coinciden y el error de
potencia será menor que un valor establecido (tolerancia). Los errores de potencia
se hallan de la siguiente forma:
Si en la barra no existe generación o carga, las potencias programadas equivalen
a cero.

ΔPi y ΔQi son los errores de potencia, cuando la potencia programada y la

calculada son muy similares entre si, se cumple que:
En este punto se pueden aceptar los voltajes correspondientes como las variables
que representan el estado real del sistema. Las potencias activas y reactiva
pueden ser expresadas de la siguiente forma:

GENERALIDADES DE LOS SISTEMAS DE DISTRIBUCION

Las redes de distribución presentan características muy particulares, que los
diferencian de los sistemas de transmisión o de potencia. Entre estas se
distinguen:

♦ Topologías radiales.
♦ Múltiples conexiones (monofásicas, bifásicas, entre otras).
♦ Cargas de distinta naturaleza.
♦ Líneas de resistencia comparables a la reactancia.
♦ Líneas sin transposiciones.

Los sistemas de distribución son generalmente radiales, esto es, el flujo de

potencia nace de una sola barra. Esta barra principal se reconoce como la
subestación (Nodo Slack) que alimenta al resto de la red, ver figura 2.2. En la
subestación se reduce el voltaje del nivel de alta tensión (A.T.), al de media
tensión (M.T.).

Generalmente se utiliza para el control de tensión en el lado de M.T. un

transformador con regulador de tensión (Taps) o en su defecto un banco de
condensadores. El cambiador automático de derivaciones en transformadores de
potencia AT/MT de sistemas de distribución permite realizar el cambio de
derivaciones con carga conectada (Load Tap Changer, LTC, o Cambiador de
Derivaciones Bajo Carga, CDBC). La tensión en barras de M.T. de la subestación
oscila normalmente entre 1 y 05.1 en p.u, por efecto del compensador por caída
de línea que posee el sistema de control de CDBC. En las horas de pico fluye
mayor corriente por las líneas, lo que provoca una mayor caída de tensión en las
mismas. Este problema se atenúa en parte con el compensador de caída de línea
que eleva la tensión en barras de la subestación de potencia en estas horas. La
distribución se efectúa luego en el nivel de media tensión (M.T.) o en baja tensión
(B.T.), dependiendo del tipo de usuario a alimentar. Los usuarios residenciales o
comerciales se alimentan en B.T. Los industriales se alimentan en M.T. o en B.T.,
según los requerimientos particulares de cada uno de ellos.
En los sistemas de distribución se pueden encontrar diferentes tipos de
conexiones: trifásicas, bifásicas o monofásicas. Si bien es cierto en M.T.
predominan las redes trifásicas es frecuente encontrar cargas bifásicas,
especialmente en zonas rurales. Sin embargo, en B.T. es donde se encuentran las
más variadas conexiones, consecuencia de una mayoría de cargas residenciales
de naturaleza monofásicas. Los desequilibrios que se generan en B.T. tratan de
amortiguarse repartiendo equitativamente las cargas en las tres fases. En este
trabajo se emplearán sistemas que operan con cargas equilibradas.

Otro aspecto importante de los sistemas de distribución es la presencia de cargas

de distinta naturaleza. Los tipos de carga que comúnmente se encuentran son:
residenciales, comerciales, industriales y agro-industriales. Cada una de estas se
caracteriza por poseer un factor de potencia particular y un determinado
comportamiento frente a las variaciones de tensión, y temperatura.

Contrariamente a lo que sucede en sistemas de transmisión, en distribución la

resistencia de las líneas es comparable a su reactancia. Generalmente la razón R
X tiene un amplio rango de variación, pudiendo llegar a ser bastante menor que 1.

Los sistemas de distribución son esencialmente radiales y no presentan nodos PV.

Este tipo de sistemas eléctricos se compone básicamente de un nodo fuente
(Slack) y nodos PQ.

En los sistemas de distribución no existen transposiciones. La causa es que en

estos sistemas las líneas son cortas (menos de 50 Km.) [15]. Esto permite que las
caídas de tensión debido a los acoplamientos entre las fases sean
desequilibradas. Por esta causa, el modelamiento más acertado de las líneas es a
través de una matriz simétrica llena, de 3× 3 . Sin embargo, en múltiples
aplicaciones es posible hacer un modelamiento monofásico equivalente como el
que se presenta.

CLASIFICACION DE LOS METODOS DE FLUJO DE CARGA

Los métodos de Flujo de Carga en sistemas de distribución se pueden clasificar en

tres grupos, a saber:

Métodos de Flujo de Carga convencionales.

Estas son metodologías diseñadas en forma general, que pueden ser aplicadas
tanto en sistemas eléctricos de potencia, como en sistemas de distribución. Entre
los métodos de Flujo de carga convencionales se encuentran:

♦ Método de Gauss-seidel:
Este método presenta la desventaja de que no siempre converge a una solución, o
de que a veces converge muy lentamente, especialmente cuando la matriz YBUS
no es diagonalmente dominante.

♦ Método de Newton Raphson:

Este es un método muy confiable, más rápido desde el punto de vista
computacional y más económico en los requisitos de almacenamiento de
información. Se basa en la expansión en series de Taylor para una función de
dos o más variables para resolver el problema de flujos de potencia.

♦ Método de Newton Raphson Desacoplado Rápido:

Para sistemas de transmisión de potencia de gran escala, representa una
alternativa para mejorar la eficiencia computacional y reducir los
requerimientos de memoria
Cabe anotar que la metodología de Newton Raphson Desacoplado Rápido fue
diseñada básicamente para ser aplicada en sistemas eléctricos de potencia, sin
embargo, se puede modificar para adaptarse a las características de los sistemas
de distribución.

Métodos de Flujo de Carga de barrido.

Existen métodos que aprovechan la característica radial de los sistemas de

distribución. Sus algoritmos están basados en un barrido atrás/adelante sobre un
sistema escalonado, denominados métodos de barrido. Entre los métodos de
Flujo de Carga de barrido se encuentran:

♦ Método de Suma de Corrientes:

Este método hace uso de manera separada de la primera y segunda ley de
Kirchhoff y de la relación iV.

♦ Método de Suma de Potencias:

El método de Suma de Potencias se basa en resolver para cada rama del
sistema la ecuación.

♦ Método de Newton Modificado:

Está basado en que la matriz jacobiana presenta la forma T UDU , donde U es
una matriz triangular superior constante, y depende de la topología de la red,
y D es una matriz diagonal.

Otros métodos de Flujo de Carga.

Existen otras metodologías para el cálculo de Flujo de Carga en sistemas de
distribución, que pueden presentar características similares a las metodologías
convencionales y de las de barrido, deben estar contenidas en una categoría
diferente:

♦ Método de Racijic-Bose:
Es una versión modificada del método de Newton Raphson Desacoplado
Rápido.
33

♦ Método de La Cadena:
Este método se basa en el hecho de que las redes radiales después de ser
ordenadas, la inversa de la matriz Jacobiana correspondiente, no genera
elementos de relleno, asi que el proceso de inversión puede efectuarse en
función de la disposición topológica de la red.

♦ Método de Rotación de ejes:

Se basa en un flujo de carga usando Newton Raphson Desacoplado Rápido
(NRDR), pero con algunas pequeñas modificaciones. El NRDR no converge
cuando se utilizan tablas de datos de sistemas eléctricos de distribución, por la
relación R
X que se presenta en estos sistemas.

♦ Método Basado en la Formulación de la YBUS :

Esta metodología fue propuesta por Young-Hyun Moon, Byoung-Kon Choi, SeHo
Kim, Bok-Nam Ha y Jung-Ho Lee en un artículo de la IEEE en Mayo de
1999.
 Clasificación de los métodos de flujo de carga

METODOS TRADICIONALES DE FLUJO DE POTENCIA

Los métodos tradicionales de flujo de carga como Gauss-Seidel, Newton Raphsón,

entre otros, son metodologías diseñadas especialmente para ser aplicadas en
sistemas eléctricos de potencia, pero mediante algunas adaptaciones pueden ser
usados en sistemas de distribución. En el desarrollo de este capítulo se
presentarán las ecuaciones que manejan los métodos de flujo de carga
convencionales y se podrá observar el algoritmo mediante el cual se efectúa el
cálculo del flujo de carga.

MÉTODO DE NEWTON RAPHSON

El método de Newton Raphson es un método diseñado especialmente para ser

aplicado en sistemas eléctricos de potencia, pero puede ser usado en sistemas de
distribución.

Este es un método muy confiable en su convergencia, más rápido desde el punto

de
vista computacional y más económico en los requisitos de almacenamiento de
información.
Se basa en la expansión en series de Taylor para una función de dos o más
variables
para resolver el problema de flujos de potencia.

Método de solución de Newton Raphson

El sistema de ecuaciones linealizado se escribe en forma completa y en
su forma compacta.
Al vector ~D se le llamo el vector de desajustes, también llamado vector de
residuos
por algunos autores. Este vector representa la diferencia entre los términos
independientes de cada ecuación, y el valor de dichos términos en función de las
incógnitas. Al vector ~C se le denomina vector de correcciones, contiene los
valores
que hay que agregar a las incógnitas de la k-esima iteración para mejorar
(corregir)
el valor anterior, en función del cual se calcularon dichos valores.

La formulación del método de Newton Raphson es directa, en esencia el problema

de
flujos consiste en calcular los voltajes nodales de la red, tomando en cuenta una
serie de restricciones que, en su expresión más simple, consiste en inyecciones
de
potencias conocidas. Dichas inyecciones constituyen las variables, mientras que
las funciones evaluadas en los valores de las incógnitas obtenidas en la
iteración k-esima, son las expresiones de las potencias.

Los elementos de dicho vector de desajustes serán igual a:

Donde las expresiones que definen a Pi y a Qi son las siguientes:

Por otro lado el vector de correcciones está compuesto por ΔVi y Δδ i .

Se puede así observar la formulación general del problema de flujos de potencia

en el método de Newton Raphson, es decir, se puede observar en términos de las
variables del problema de flujos de potencia.
Donde se muestran detalladamente los renglones que corresponden a la i-esima
barra, en el vector de desajustes, y su interacción con la m-esima barra, en el
vector de correcciones. Los elementos de la matriz Jacobiana muestran los
elementos correspondientes a dicha interacción.

Se observa que en el caso de los barras PQ , se asignan ambos elementos en el

vector de desajustes, ya que se conocen las inyecciones de potencia real y
reactiva. En las barras PQ , son incógnitas la magnitud y el angulo de voltaje, por
esta razón ambos aparecerán en el vector de correcciones, para este tipo de
barra. De esta forma resultan dos ecuaciones para cada barra PQ .
El número de ecuaciones que conforman el modelo matemático de flujos de
potencia
en el Newton Raphson será: 2n=2(n-1) . Para el nodo compensador no existe la
necesidad de escribir ecuación, ya que no se conocen las inyecciones de potencia
y
además la magnitud y el Angulo de voltaje en esta barra no son una incógnita.

La formulación anterior se conoce como formulación polar, debido a que las

variables
se expresan en forma polar. Existe otra formulación, denominada formulación
rectangular, sin embargo, no es tan popular como la formulación polar, debido a
que
esta es más eficiente en general.

A continuación, se desarrollan las expresiones correspondientes a los elementos

del
vector de desajustes y de la matriz Jacobiana.
Inicialmente se define el formato polar de voltajes y admitancias:

Las expresiones de las cantidades que forman el vector de desajustes se

definieron
previamente, las cuales al ser combinadas dan
como resultado las expresiones finales.
Se debe notar que el término Vi se introdujo dentro de la suma, debido a que el
índice de esta es k , y por lo tanto no se produce ninguna alteración realmente en
la
expresión.
Para desarrollar la expresión de la matriz Jacobiana, se definen las variables
matriciales del modelo como se indica:

La expresión matricial anterior implica las siguientes definiciones:

Teniendo presente lo anterior, se obtiene:

Elementos de [H]:

Elementos de [N]:

Elementos de [J]:

Elementos de [L]:
El proceso iterativo asociado se puede representar por la ecuación matricial. Que
muestra la ecuación de Newton Raphson en la iteración l-esima.

Conclusión
En este trabajo se presentó un estudio comparativo sobre diferentes metodologías
de cálculo de flujo de potencia para sistemas de distribución. El cual permite
observar el comportamiento de los diferentes métodos y así poder determinar ante
qué condiciones del sistema es más adecuado aplicar una u otra metodología.