FLUJO DE POTENCIAS.

1. Introducción.

En análisis de flujo de potencias (o flujo de cargas) constituye la herramienta principal para

estudiar un sistema eléctrico de transporte o distribución de energía eléctrica en régimen
permanente y poder tomar decisiones relativas a:

• Explotación del sistema (simulación de acciones posibles sobre el sistema existente).

• Planificación del sistema (estable futuros).

El problema se puede formular así: Tomando como datos las potencias demandadas por los
consumidores y las suministradas por los generadores, determinar las tensiones en cada nudo
del sistema, en régimen permanente y equilibrado, así como los flujos de potencia activa y
reactiva a través de los distintos elementos del sistema (líneas, transformadores, reactancias y
condensadores).

La resolución del flujo de cargas tiene distintas etapas:

• La primera etapa consiste en la determinación de las tensiones complejas en todos los

nudos del sistema. Dado que los datos de cargas vienen generalmente proporcionados
en términos de potencia y no como impedancia, no es posible aplicar los métodos de
cálculo basados en corrientes de nudos y tensiones de mallas. Por otra parte, los
generadores son fuentes de potencia, no fuentes de tensión o de corriente. Todo ello
da lugar a que la formulación del problema conste de un sistema de ecuaciones no
lineales.
• La segunda etapa consiste en la realización de unos cálculos rutinarios, a partir de los
resultados de la primera etapa, que conducen a la determinación de los flujos de
potencia activa y reactiva, pérdidas nen los elementos, etc.

2. Modelo de red.

Como paso previo al planteamiento y resolución de problemas, se presenta el modelo de red

empleado para este tipo de estudios, o más bien, los modelos empleados para representar los
distintos elementos del sistema.

Consideremos el sistema de la figura. Se trata de una red mallada, habitual en los sistemas de
transporte.

Esquema unifilar de un sistema de energía eléctrica

1
 • Se asume que las potencias de las cargas 𝑆⃗D i son conocidas y su valor constante.
• En la mayoría de los casos se suponen ideales los transformadores de las centrales
generadoras.

Así, las potencias 𝑆⃗G i consideradas son las suministradas a los nudos desde los lados de
alta de los transformadores correspondientes. Debido a ello, en los esquemas unifilares
solo se representa el símbolo del generador, en lugar del conjunto generador-
transformador.

• Para los transformadores no se suelen considerar los desfases introducido por el tipo de
conexión, lo que equivales a considerar un índice horario h igual a cero. Se desprecian
pérdidas en el hierro y la corriente de magnetización.

Según se ha visto, para los transformadores con relación de transformación igual a la

nominal en el sistema por unidad, el circuito equivalente queda reducido a una
impedancia serie conectada entre los nudos a los que se conecta el primario y el
secundario. Si además despreciamos las pérdidas en el cobre, la impedancia se reduce
a una reactancia.

• Las líneas de interconexión se representan por su modelo en régimen permanente

concretado en el circuito equivalente en π.

En las líneas de transporte es habitual despreciar la resistencia serie. La conductancia se

desprecia en todos los casos.

En la siguiente tabla se muestran de forma resumida los modelos usados para cada elemento de
interconexión (transformador, línea) entre los nudos. Para cada uno de ellos se muestra su
representación gráfica, las ecuaciones que reflejan su comportamiento en relación con las
variables eléctricas en sus extremos y el circuito equivalente en π. Es importante destacar que
las ecuaciones tienen la forma: I = YU, donde aparecen los parámetros admitancia.

2
En la figura anterior, representativa del esquema unifilar de un sistema eléctrico, se entenderá
que se está considerando el sistema a nivel de transporte. En esa situación, los generadores
serán las centrales que generan la electricidad, las cargas son los grandes consumidores
industriales y/o las subestaciones que interconectan con la red de subtransporte.

Se podría aplicar la misma técnica a sistemas eléctricos con otros niveles de tensión. Así, en el
nivel de subtransporte, los nudos en los que se representa un generador serán aquellos en los
que se interconecta con la red de transporte mediante subestaciones, las cargas serán los
consumidores industriales conectados a este nivel de tensión y/o las subestaciones de
interconexión con la red de distribución.

El estudio se realizará considerando el sistema eléctrico trifásico en régimen permanente y

equilibrado, lo que permitirá hacer un análisis por fase.

Se define la potencia compleja inyectada en un nudo como:

𝑆⃗ i = 𝑆⃗G i - 𝑆⃗D i (1)

La potencia 𝑆⃗ i es lo que queda de potencia suministrada por el generador i, 𝑆⃗Gi, una vez sustraída
la potencia demandada en el mismo nudo i.

Se define la corriente inyectada en un nudo como:

𝐼⃗ i = 𝐼⃗G i - 𝐼⃗D i (2)

Siendo 𝐼⃗Gi la corriente suministrada por la fase a del generador del nudo i, 𝐼⃗D i es la corriente
demandada por la fase a de la carga conectada al nudo i y, por tanto, 𝐼⃗ i es la corriente de la
fase a inyectada en el sistema desde el nudo i.

3. Matriz de admitancias de un sistema eléctrico.

Consideremos la siguiente figura que representa el esquema unifilar de un sistema eléctrico de

4 nudos, que tiene generador y carga en cada nudo.

En la siguiente figura se muestra el circuito equivalente por fase del sistema de la figura anterior,
donde las líneas se han representado mediante su equivalente en π

3
 ✓ Aplicando el método de análisis de las tensiones de nudos,

✓ tomando el neutro como nudo de referencia, y

✓ teniendo en cuenta las ecuaciones que rigen el funcionamiento de los distintos

elementos de interconexión del sistema,

podemos escribir la siguiente expresión matricial:

(3)

donde:

4
Cada término 𝑌⃗⃗ii (i = 1, 2, 3,4) se denomina admitancia propia del nudo i y es igual a la suma
algebraica de todas las admitancias unidas directamente a ese nudo.

Cada término 𝑌⃗⃗ik (i, k = 1, 2, 3,4) se denomina admitancia mutua entre los nudos i y k, y es igual
a la suma (con signo negativo) de todas las admitancias conectadas directamente entre esos
nudos. Además, siempre que no haya transformadores desfasadores se cumple: 𝑌⃗⃗ik = 𝑌⃗⃗ki.

La expresión matricial anterior se puede escribir de forma compacta de la siguiente manera:

𝐼⃗i = ∑𝑛𝑘=1 𝑌
⃗⃗ ik 𝑈
⃗⃗K i = 1, 2, 3, …, n (4)

O en forma de matriz:

IBUS = YBUS UBUS

YBUS se conoce como matriz de admitancia de nudos. La dimensión de esa matriz es n x n, siendo
n el número de nudos independientes de la red. El número total de nudos siempre es n + 1,
incluyendo en nudo de referencia (neutro).

4. El problema del flujo de potencias.

Como ya se ha dicho, el problema del flujo de cargas o de potencias no es otro que la

determinación de las tensiones en cada uno de los nudos de un sistema trifásico de energía
eléctrica en régimen permanente y equilibrado, Una vez conocidas las tensiones en los nudos se
puede determinar el flujo de potencia activa y de potencia reactiva en las líneas y
transformadores, así como las pérdidas del sistema.

En la siguiente figura se muestra una fase de un nudo genérico i con su fuente de potencia
equivalente, suministrando al nudo una corriente 𝐼⃗i , una potencia aparente compleja 𝑆⃗ i, siendo
⃗⃗i la tensión en el mismo.
𝑈

La potencia compleja entrante en el nudo i desde la fuente (inyectada) será:

𝑆⃗ i = 𝑈
⃗⃗i 𝐼⃗i * i = 1, 2, …, n

donde n es el número de nodos del sistema, sin incluir el neutro.

Teniendo en cuenta la expresión (4) y sustituyendo en la anterior, tendremos:

(5)

Utilizando la siguiente notación:

5
Obtenemos el siguiente sistema no lineal de n ecuaciones complejas.

(6)

Que se puede descomponer en 2n ecuaciones reales siguientes:

(7.1)

(7.2)

Si observamos las ecuaciones (7) vemos que a cada nudo i están asociadas 2 ecuaciones y cuatro
magnitudes eléctricas: el módulo de la tensión Ui, el argumento δi, la potencia activa inyectada
Pi, y la potencia reactiva inyectada Qi.

El planteamiento preciso del problema se basa en consideraciones de índole práctica, pero

también en consideraciones matemáticas. Se desarrollan a continuación estas consideraciones.

4.1. Tipos de nudos.

Los nudos del sistema eléctrico se clasifican en tres categorías: Nudos degeneración o PV, nudos
de carga o PQ y nudo oscilante o slack.

En todos los nudos se supone conocida la potencia demandada: S⃗⃗ espDi = P espDi + jQ espDi

La notación “esp” hace referencia a datos “especificados”.

Nudos de generación o nudos PV.

En un nudo de este tipo, la potencia activa PGi y el módulo de la tensión Ui son conocidos, esto
es, se consideran datos.

En consecuencia, son conocidos P espi = P espGi – P espDi y U espi

Son incógnitas Q i y δ i .

Nudos de carga o nudo PQ.

Son los nudos en los que se conoce el consumo de tanto de potencia activa como de potencia
reactiva y en los que o bien no hay generación o si la hay esta es fija y se toma como dato
(pequeños generadores sin regulador de tensión).

Por tanto, son datos P espi = P espGi – P espDi y Q espi = Q espGi – Q espDi

6
Son incógnitas Ui y δi.

En la realidad el 80-90% son nudos PQ.

Cabe destacar que en algún nudo de carga puede haber conectados condensadores, el nudo
será identificado como nudo PQ si los condensadores suministran potencia reactiva fija, pero se
considerará nudo PV si son utilizados para mantener una tensión especificada.

Nudo oscilante o nudo slack.

Si de todos los nudos del sistema, una parte fueran PV y el resto PQ, resultaría que todos los
valores Pi serían datos del problema, lo cual es incompatible con la exigencia de balance de
potencia activa.

En concreto, el balance de potencia activa viene dado por la expresión:

(8)

Las pérdidas de activa, Ppérd no pueden saberse de antemano, por lo que no todas las Pi pueden
ser datos, al menos una debe dejarse como incógnita.

Como todas las potencias PDi son datos, en necesario que quede algún generador con la PGi sin
especificar. Al nudo que se le asigna esa característica se le identifica como nudo 1 y se le
denomina nudo oscilante o slack. En el nudo slack son datos U1esp y δ1esp = 0 ( se toma como
ángulo de referencia) e incógnitas P1 y Q1.

En la siguiente se indican los tipos de nudos:

Tipos de nudos

4.2. Formulación precisa.

Se enumeran los nudos de tal forma que el nudo 1 es el nudo oscilante, los nudos 2, 3, …, m son
los nudos PV y los nudos m +1, …, n son los nudos PQ. La formulación del problema es entonces
la siguiente:

⃗⃗1, (P2, U2), …, (Pm, Um), 𝑆⃗m+1, …, 𝑆⃗n

Datos: 𝑈

Hallar: 𝑆⃗1, (Q2, δ2), …, (Qm, δm), 𝑈

⃗⃗m+1, …, 𝑈
⃗⃗n

Se trata de resolver un sistema de 2n ecuaciones algebraicas reales no lineales con 2n incógnitas

reales. La obtención de la solución tiene dos etapas diferentes.

Primera etapa:

7
Consiste en realizar un proceso iterativo para resolver las ecuaciones definidas por (7.1) y (7.2)
que relacionan las incógnitas.

En el proceso iterativo no se considera las dos ecuaciones al nudo 1, que proporcionan de forma
explícita las incógnitas (P1 y Q1) en función de las tensiones de los nudos. Por tanto, debe
resolverse un sistema con 2n-2 ecuaciones, implícitas no lineales, para determinar 2n-2
incógnitas.

Del mismo modo, las ecuaciones definidas por (7.2) correspondientes a los nudos PV (i = 2, 3, …,
m) proporcionan de forma explícita las potencias reactivas que son incógnitas en esos nudos.

En la resolución del proceso iterativo debe tenerse en cuenta una serie de restricciones
impuestas por las limitaciones físicas de los equipos y por las condiciones de funcionamiento
adecuado del sistema. Estas son:

• Los módulos de las tensiones en los nudos han de cumplir la condición:

Uimin ≤ Ui ≤ Uimáx i = 1, 2, 3, …, n

Los equipos están diseñados para funcionar con una tensión nominal, siendo posible
oscilaciones dentro de un margen de tolerancia. En el nudo 1 y los nudos PV (2, 3, …, m)
el módulo de la tensión es dato, y se expresará dentro de ese margen de tolerancia. En
los nudos PQ (m+1, …, n), donde los módulos de la tensión son incógnitas, se
comprobará que se cumple la restricción una vez resuelto el proceso iterativo,
concluyéndose que si no se cumple el funcionamiento del sistema no sería correcto y
habría que introducir cambios en las variables consideradas datos de partida.

• Los argumentos de las tensiones en los nudos han de cumplir la condición:

│δi – δk│≤ │δi – δk│máx i, k = 1, 2, 3, …, n

Esta condición equivale a decir que las potencia activas, a través de líneas (Pik), no debe
superar un valor máximo impuesto. Si no se cumple el funcionamiento del sistema no
sería correcto y habría que introducir cambios en las variables consideradas datos de
partida

• Debido a las limitaciones físicas de los generadores, los valores de PGi están sometidos a
las restricciones:

PGimin ≤ PGi ≤ PGimáx i = 1, 2, 3, …, n

En los nudos PV (2, 3, …, m) las potencias activas son dato, deben estar dentro de ese
margen de tolerancia, teniendo en cuenta además el balance de potencia activa:

(P activas demandadas = P activas generadas + pérdidas)

Pero este ajuste no podrá hacer hasta que este resulto el problema, es decir,
determinada P1 y las pérdidas.

• Debido a las limitaciones físicas de los generadores, los valores de Q Gi están sometidos
a las restricciones:

Q Gimin ≤ Q Gi ≤Q Gimáx i = 1, 2, 3, …, n

8
 En los nudos PV (2, 3, …, m) las potencias reactivas son incógnitas, debe tenerse en
cuenta esta restricción en el proceso iterativo de resolución. En concreto, en cada
iteración debe comprobarse si se cumple o no la restricción.

Segunda etapa:

Consiste en calcular lo siguiente:

Interconexión entre dos nudos

1. Las potencias en el nudo slack (nudo 1) mediante la ecuación correspondiente al

nudo 1:

(9)

2. El flujo de potencias en la líneas teniendo en cuenta su modelo en π mediante.

(10)

3. Las pérdidas totales del sistema, que pueden calcularse mediante la suma de las
inyecciones de todos los nudos, esto es:

(11)

5. Método de Gauss-Seidel.

Se trata de la aplicación a sistemas no líneas de una técnica de resolución iterativa de sistemas

lineales. Apenas se utiliza en la actualidad dado que presenta de problemas de convergencia y
necesita muchas iteraciones para llegar a la solución aceptable.

Caso I. Consideremos en primer lugar el problema del flujo de potencias para un sistema en que
no hay focos PV, esto es, todo los nudos son PQ, excepto uno que es el nudo oscilante.

⃗⃗1, 𝑆⃗2, 𝑆⃗3, …, 𝑆⃗n

Datos: 𝑈

Hallar: 𝑆⃗1, 𝑈
⃗⃗2, 𝑈
⃗⃗3, …, 𝑈
⃗⃗n

Hacemos uso de la ecuación (5), y tendremos:

(12)

9
(13)

La ecuación (12) proporciona el valor de la potencia compleja 𝑆⃗1, una vez sean conocidas todas
las tensiones de los nudos 𝑈 ⃗⃗1, 𝑈
⃗⃗2, …, 𝑈
⃗⃗n. Como 𝑈
⃗⃗1 es un dato, solamente hay que hallar 𝑈
⃗⃗2, 𝑈
⃗⃗3,
…, 𝑈⃗⃗n. Así tenemos n-1 incógnitas complejas, que se hayan al resolver n-1 ecuaciones complejas
dadas por (13).

En términos de variables reales serían 2n – 2 incógnitas, que se hayan al resolver 2n – 2

ecuaciones.

Las ecuaciones complejas dadas por (13) se pueden escribir de la siguiente manera:

(14)

⃗⃗i, y separando el término de 𝑌

Dividendo por 𝑈 ⃗⃗ii queda:

(15)

Así se tienen n -1 ecuaciones complejas, implícitas no lineales, de donde se pueden determinar

⃗⃗i de la forma matricial:
n -1 incógnitas 𝑈

x = F (x)

donde x es el vector formado por las variables incógnitas.

El método de Gauss-Seidel permite resolver un sistema de ecuaciones de esta característica

mediante un proceso iterativo. El proceso se describe en los siguientes puntos:

• Para empezar las iteraciones se toma un conjunto de valores de tensiones de partida.

⃗⃗i = 1 ∠0º, siendo i = 1, 2, …,n
Típicamente 𝑈

• Se calculan los valores de las tensiones mediante el proceso iterativo definido por la
expresión:

(16)

10
 En las incógnitas aparece el superíndice que hace referencia al número de iteración. El
valor r = 0 se reserva para la iteración con los valores iniciales o de partida. En la iteración
r + 1 se calculan los valores de las incógnitas 𝑈 ⃗⃗i(r + 1) a partir de los valores
correspondientes en la iteración r.

• El proceso continúa hasta que se cumpla una condición de parada, como puede ser que
se satisfaga la desigualdad siguiente:

(17)

Siendo ε un valor prefijado.

Es decir, se para el proceso de iteración cuando la diferencia entre las tensiones de dos
iteraciones consecutivas, en valor absoluto, sea menor que la cantidad ε para todos los
nudos (típicamente ε = 10-4).

Al resolver estas ecuaciones se procede en el orden de los nudos 2, 3, … Entonces cuando se va

a calcular⃗⃗⃗⃗
𝑈3(r + 1) ya antes se ha calculado 𝑈 ⃗⃗2(r + 1). Como presumiblemente 𝑈 ⃗⃗2(r + 1) es una mejor
estimación que 𝑈 ⃗⃗2(r), parece lógico usar 𝑈 ⃗⃗2(r + 1) en lugar de 𝑈 ⃗⃗2(r). Del mismo modo, cuando
vayamos a calcular 𝑈 ⃗⃗4(r + 1), pueden usarse 𝑈
⃗⃗2(r + 1) y 𝑈⃗⃗3(r + 1).

Este razonamiento permite modificar el método denominado Gauss, definido por la expresión
(16), que pasa a llamarse Gauss-Seidel, definido por la expresión:

(18)

Caso II. Consideremos ahora el problema general en el que existen nudos PV y nudos PQ. En este
caso, tenemos.

⃗⃗1, (P2, U2), …, (Pm, Um), 𝑆⃗m +1, 𝑆⃗m + 2, …, 𝑆⃗n

Datos: 𝑈

Hallar: 𝑆⃗1, (Q2, δ2), …, (Qm, δm), 𝑈

⃗⃗m+1, 𝑈
⃗⃗m+2, …, 𝑈
⃗⃗n

Para resolver el problema hay que hacer algunos cambios en lo visto en el caso anterior. La
expresión para las iteraciones se escribe de forma más adecuada así:

(19)

11
Es la misma expresión que la (18), pero ahora aparece Qi(r + 1) en lugar de Qi.

Para los nudos de carga o PQ (i = m+1, …, n) resulta que Pi y Qi son datos y el proceso iterativo
será el mismo que se ha descrito en el Caso I, por lo que se puede ignorar el superíndice aplicado
a Qi.

Para los nudos de generación o PV (i = 2, …, m) las potencias reactivas Qi son incógnitas, siendo
funciones explícitas de la tensiones de nudos y, por tanto, pueden ser calculadas después de
cada iteración r. Así, de las ecuaciones (14) dadas por:

Se obtiene:

(20)

Teniendo en cuenta, como se ha dicho en el Caso I, que las ecuaciones se resuelven siguiendo el
orden de los nudos 2, 3, …, es posible usar los valores de las tensiones más actualizados,
entonces la expresión para la determinación de las Qi en cada iteración será:

(21)

y estos serán los valores que se introducen en (19).

Como en los nudos PV el módulo de la tensión es un dato, lo que se calcula con (19) es el
argumento de la tensión, esto es:

(22)

y las tensiones correspondientes (módulos y argumentos) serán:

(23)

12
En los nudos PV, en cada iteración hay que comprobar si Qi está dentro de los límites prácticos
definidos por Q Gimin ≤ Q Gi ≤Q Gimáx i = 1, 2, 3, …, n. Solo entonces tiene sentido determinar δi(m+1)
quedando perfectamente definido 𝑈 ⃗⃗i(r+1). Si no se cumple la restricción, se hace que se cumpla
imponiendo un valor para Qi deforma que el nudo se trata como si fuera PQ pasando a ser
incógnita la tensión(módulo y argumento) en ese nudo.

6. Método de Newton-Raphson.

6.1. Concepto matemático.

Consideremos una función de una sola variable f (x), Su representación gráfica puede ser, por
ejemplo así:

Si hacemos f(x) = 0 tendremos una ecuación cuya solución exacta nos daría las raíces de dicha
ecuación, o lo que es lo mismo, el punto de corte de la función con el eje de abscisas (eje x).

Podemos hacer un desarrollo en serie de Taylor de la función f(x) en el entorno de un punto

cualquiera, por ejemplo, x(0), que puede considerarse una aproximación o estimación de la
solución exacta.

La ecuación tendrá entonces la forma siguiente:

(24)

𝑓´(𝑥 (0) ) 𝑓´´(𝑥 (0) )

𝑓(𝑥) = 𝑓 (𝑥 (0) ) + (𝑥 − 𝑥 (0) ) + (𝑥 − 𝑥 (0) )2 + ⋯ = 0
1! 2!
Haciendo un cambio de variable definido por Δx(o) = x – x(0), que representa la diferencia entre la
solución exacta y la estimación dada por x(0), la ecuación anterior se puede expresar así:

(25)

𝑓´(𝑥 (0) ) (0) 𝑓´´(𝑥 (0) )

𝑓(𝑥) = 𝑓 (𝑥 (0) ) + 𝛥𝑥 + (𝛥𝑥 (0) )2 + ⋯ = 0
1! 2!
Si consideramos todos los términos del desarrollo en serie, la resolución de la ecuación
proporciona la solución exacta.

Ahora bien, si consideramos solo los dos primeros términos del desarrollo en serie de Taylor,
tendremos la ecuación siguiente:

13
(26) 𝑓(𝑥 (0) ) + 𝑓´(𝑥 (0) )(𝑥 − 𝑥 (0) ) = 0

O bien:

(27) 𝑓 (𝑥 (0) ) + 𝑓´(𝑥 (0) )𝛥𝑥 (0) = 0

Ambas expresiones son equivalentes, en la ecuación (26) la incógnita es x, en la ecuación (27) la

incógnita es 𝛥𝑥 (0) , estando relacionadas ambas incógnitas por la relación Δx(o) = x – x(0). La
ecuación, en cualquiera de sus formas, no proporcionará la solución exacta ya que solo hemos
considerados unos pocos elementos del desarrollo en serie de Taylor (como se ha dicho, de
considerarse el desarrollo de Taylor completo, esto es, todos sus términos, si se obtendría la
solución exacta).

En cualquier caso, se considere la ecuación de una u otra forma, el primer término de la igualdad
es la recta tangente a la curva f(x) en el punto x(0), por lo que la ecuación correspondiente que
resulta al igualar a cero permite determinar el punto de corte de la recta con el eje de abscisas
(eje x), determinándose así una nueva raíz, x(1), que es una mejor aproximación a la solución
exacta tal como se muestra en la siguiente figura:

Siendo:

(28)

x(1) = x(0) + Δx(0)

Ordenando las ecuaciones (26) y (27), se pueden escribir así:

(29)

−𝑓´(𝑥 (0) )(𝑥 − 𝑥 (0) ) = 𝑓(𝑥 (0) )

(30)

−𝑓´(𝑥 (0) )𝛥𝑥 (0)= 𝑓 (𝑥 (0) )

El proceso descrito se puede repetir, ahora haciendo el desarrollo en serie de Taylor en el

entorno de x(1), que es una mejor aproximación a la solución exacta. Considerando solo los dos
primeros términos del desarrollo en serie y considerando Δx(1) = x – x(1), obtendríamos las
siguientes ecuaciones equivalentes entre sí:

(31)

−𝑓´(𝑥 (1) )(𝑥 − 𝑥 (1) ) = 𝑓(𝑥 (1) )

14
(32)

−𝑓´(𝑥 (1) )𝛥𝑥 (1)= 𝑓 (𝑥 (1) )

Ambas representan la ecuación de la recta tangente a la curva f(x), que al igualarla a cero,
permitiría obtener la raíz x(2), que no es la solución exacta pero si es una mejor aproximación a
la misma:

(33) x(2) = x(1) + Δx(1)

El procedimiento iterativo descrito permite obtener la mejor aproximación a la solución exacta

de una ecuación.

Tal procedimiento puede generalizarse para sistemas de n ecuaciones con variables x1, x2, …, xn,
si bien no es posible una representación gráfica al igual que la realizada para una ecuación de
una variable.

Sea el sistema de ecuaciones algebraicas no lineales:

Supóngase que x1(o), x2(o), … xn(o) es una estimación inicial de la solución y que Δx1(o), Δx2(o), …,
Δxn(o) son las desviaciones de las variables que sumadas a las soluciones iniciales estimadas dan
la solución exacta. Como es evidente se cumple:

Cada una de las ecuaciones del sistema anterior puede desarrollarse en serie de Taylor en el
entorno del punto correspondiente a la solución estimada inicial, es decir, en el entorno del
punto x1(o), x2(o), … xn(o), teniendo:

(34)

Siendo (∂fi / ∂xj)(0), …, (∂fi / ∂xn)(o) son las derivadas parciales de fi respecto de x1, x2, …, xn ,
evaluadas en x1(o), x2(o), … xn(o).

15
Nótese que solo estamos considerando las derivadas de primer orden, es decir, despreciamos
las derivadas de orden superior, lo que implica que la solución del sistema no será la solución
exacta, solo una mejor aproximación.

Las ecuaciones con el desarrollo de Taylor se pueden escribir de forma matricial así:

(35)

Que en notación compacta se puede poner:

(36)

Donde:

es el vector de variables transpuesto.

es el vector de funciones transpuesto.

Y J es la matriz Jacobiana o simplemente Jacobiano del sistema de ecuaciones, que se obtiene

derivando las funciones f con respecto a x y particularizando para x(o). La expresión representa
un sistema de acuaciones lineales donde las incógnitas son Δx1(o), Δx2(o), …, Δxn(o).

Se trata de un sistema de acuaciones algebraicas lineales que se puede resolver mediante

triangulación o retrosustitución.

Sumando los valores obtenidos a las estimaciones iniciales se obtiene:

xi(1) = xi(0) + Δxi(o) i = 1, 2, 3, …, n

Como hemos dicho, solo consideramos los primeros elementos del desarrollo de Taylor, así que
x1(1), x2(1), …, xn(1) no serán las soluciones exactas, serán una mejor solución que x1(0), x2(0), …, xn(0).
Al identificarse con el superíndice (1) se entiende que es el resultado de la primera iteración. Si
se repite el proceso, a partir de la última estimación de la solución, se tiene, en general, el
algoritmo iterativo expresado por:

Donde Δxi(r) es la solución que se obtiene de una ecuación semejante a la (35) o la (36), de forma
compacta, cambiando el superíndice (0) por (r), esto es:

16
(37)

Las iteraciones terminan cuando se cumpla:

(38)

Siendo ε un valor prefijado.

6.2. Aplicación al problema de análisis del flujo de potencias.

En el caso general de que existan nudos PV y nudos PQ, el vector de incógnitas será:

(39)

Las ecuaciones al resolver en un problema de flujo de cargas son las ecuaciones (7.1) y (7.2) ya
vistas:

excluyendo las correspondientes al nudo oscilante y las correspondientes a la potencia reactiva

de los nudos PV, pero expresadas de la siguiente manera:

(40)

donde los denominados residuos de potencias son:

(41)

En este planteamiento tenemos un sistema de 2n – m – 1 ecuaciones reales en las 2n – m – 1

incógnitas que formas el vector x.

El sistema de ecuaciones dado por (37) se puede escribir así:

(42)

17
O de forma compacta:

(43)

La utilización de ΔUi /Ui facilita simplifica los términos del Jacobiano. La matriz Jacobiana está
compuesta por cuatro submatrices J1, J2, J3 y J4. Las derivadas parciales calculadas a partir de
(41) se reflejan en la siguiente tabla:

Observando las ecuaciones de las tabla se deduce una importante propiedad de la matriz
jacobiana:

En cada submatriz, los términos fuera de la diagonal principal (i ≠ j) contiene el valor del elemento
serie del circuito equivalente en π de la línea que interconecta los nudos. Si no existe conexión
directa entre los dos nudos, el elemento correspondiente de la matriz YBUS es nulo y, por tanto,
el término correspondiente de J también es cero.

El procedimiento o proceso para resolver el `problema del flujo de potencias por el método de
Newton-Raphson, se concreta en los siguientes pasos:

• PASO 1. La primera iteración se realiza a partir de una estimación inicial de la solución,

típicamente 𝛿⃗i(0) = 0, 𝑈
⃗⃗i(0) = 1 (perfil plano)

• PASO 2. Calcular el valor de las funciones [ΔP(x)│ΔQ(x)]T para x(r).

• PASO 3. Calcular la matriz Jacobiana J(r) para x(r) utilizando la Tabla 6-V.
• PASO 4. Calcular funciones [Δδ│ΔU/U]T del sistema de ecuaciones algebraicas lineales
dado por (42) o (43).
• PASO 5. Actualizar variables:

18
 (44)

• Parar el proceso iterativo cuando se cumpla:

(45)

Una vez terminado el proceso iterativo descrito serán conocidas las tensiones complejas
(módulos y argumentos) en todos los nudos. Conocidos tales valores se sustituirán en las
ecuaciones (9), (10) y (11), que proporcionarán los valores de las incógnitas correspondientes
al nudo oscilante P1 y Q1, los flujos de potencias en las líneas 𝑆⃗ik y las pérdidas.

Del mismo modo, la sustitución de las tensiones complejas en las ecuaciones dadas por (7.2)
correspondientes a los nudos PV (i = 2, 3, …, m) nos proporciona el valor de las potencias
reactivas incógnitas en dichos nudos. Ahora bien, el dispositivo que regula la tensión en los
nudos PV tiene una capacidad limitada para inyectar o absorber potencia reactiva, si se alcanza
algunos de los límites Qimin o Qimáx la tensión en el nudo regulado no puede mantenerse en el
valor Uiesp. Esto obligaría a empezar de nuevo el problema cambiando las tensiones especificadas
en los nudos PV, lo cual no parece práctico. Es más adecuado determinar la potencia reactiva en
cada nudo PV después de cada iteración, introduciendo los valores de las incógnitas calculados
en esa iteración en las ecuaciones correspondientes:

(46)

Si el valor calculado Qi en alguno de los nudos excede los límites, se hace que cumpla el límite
asignando un valor concreto coincidente con el valor del límite rebasado, esto es, Qiesp = Qimin o
Qiesp = Qimáx, con lo cual este nudo pasará a ser considerado nudo PQ, siendo entonces incógnita
la tensión correspondiente (módulo y argumento). Como consecuencia de pasar de ser nudo PV
a ser nudo PQ, es necesario reintroducir ΔQi en el vector de funciones e ΔUi en el vector de
incógnitas, cambiando entonces la matriz Jacobiana. Si en la iteración posterior (r) resulta que
Ui(r) > Uiesp con Qiesp = Qimáx o bien Ui(r) < Uiesp con Qiesp = Qimin , el nudo pasará a ser considerado
de nuevo PV. El programa de cálculo avisa si la solución se alcanza manteniéndose la violación
de límites.

19