CAP

SOLUCIO

5.1
problema

ITULO

N AL PROBLEMA DE FLUJO DE POTENCIA

Definicion

del

Si en un sistema de potencia se especican los valores de carga nodal, se conocen los para
metros de la red y ademas se conocen dos de sus cuatro variables caractersticas (V, , P,
Q), es posible determinar el estado de equilibrio nodal y global del sistema. Esta tarea se
denomina ujo de potencia.
Los elementos de la red de transmision son disenados para condiciones de operacion
en equilibrio, ademas la operacion de los mismos es efectuada de manera balanceada debido
a que la carga nodal es equilibrada, por lo tanto, el ujo de potencia a traves de los elementos
presenta condiciones de balanceo, esto es, no existe circulacion de corriente por neutro ni tierra,
en estas condiciones, si la red trifasica es descompuesta en tres redes de secuencia,
denominadas de secuencia positiva,negativa y cero, solo la red de secuencia positiva tiene
efecto en la operacion.
Por lo anterior, al efectuar analisis en estado estacionario en redes de alta potencia, la red
se modela por su equivalente monof sico de secuencia positiva.
Los estudios de Flujos de Potencia se aplican en:

Planeacion y diseno del Sistema Electrico.

 Determinacion de las condiciones operativas. El ujo de Potencia proporciona informacio

n de suma importancia acerca de la operacion de las redes electricas como es:
Voltajes nodales, lujos de potencia por las lneas, perdidas, etc.
Las matrices de redes empleadas en los diferentes estudios de sistemas electricos son:

La matriz admitancia de barra Ybus

La matriz impedancia de barra Zbus

Los metodos clasicos mas comunmente empleados de ujo de carga son desarrollados
con base en la matriz Ybus . Los de corto circuito estan basados en la matriz Zbus .
Los elementos de la matriz Ybus son descritos de la siguiente manera:
YiJ = |YiJ | iJ = |YiJ |CosiJ + J |YiJ |SeniJ
YiJ = GiJ + J BiJ
El modelo matem tico del sistema para resolver el problema de ujo de carga se expresa en
funcion de voltajes nodales, los cuales podran ser dados en componentes polares o
rectangulares, siendo la mas usada en los escritos tecnicos la representacion polar.
Vi = |Vi | i = |Vi |(Cosi + JSeni )
La ecuacion que relaciona la corriente inyectada en el nodo en funcion de voltajes
nodales y la matriz Ybus es la siguiente:
n
I = Y V + Y V + + Y V
=
i
i1 1
i2 2
in n

n=1

V
in n

Sean Pi y Qi las potencias real y reactiva especicadas en el nodo i. La ecuacion

que relaciona estas variables con el voltaje nodal y la matriz Ybus es:
P J Q = V n
i
i
V
n

in

n=1

Expresado en forma polar es:

P JQ =n
i
i

|Y

n=1

V V | +
in n i in
n
i

El error de potencia en el nodo (i) se calcula como la diferencia entre la potencia inyectada
por el sistema externo e interno, esta diferencia se presenta en el proceso del ujo de carga en el
cual la potencia del sistema externo es establecida en un valor, la del sistema interno se calcula
con base en los voltajes nodales y parametros de la red, as que solo hasta cuando se
alcance el voltaje de equilibrio- las dos potencias coinciden y el error de potencia sera
menor que un valor establecido, calculado as:

P ) P
P = (P
=P
i
Gi
Di
i calc
i prog
P

i calc

Q ) Q
Q = (Q
=Q
i
Gi
Di
i calc
i prog
Q

i calc

Si no existe generacion o carga sus terminos correspondientes equivalen a cero.

Pi y Qi : Errores que ocurren durante el desarrollo del ujo de potencia. Cuando la
potencia
programada y la calculada se aproximan entonces: |Pi | y |Qi |

Las ecuaciones de potencia activa y reactiva inyectada pueden ser expresadas en diferente
forma:

PGi

P i1
P i2

P ij

PDi
Pi programada

Pi calculada
1

QG i

Q i1
Q i2

Q ij

QD i
Qi programada

Qi calculada

Figura 5.1: Potencia activa y reactiva inyectada

1. Primera
P +JQ = V
k

P + J Q = (e + J f n )
k
k
k
k

m=1

Pk + J Qk =

(e

k
B
)
km

n
k

(G

km

J B )(e J f )
km m
m

f B
Jf G
Je
+ J f )(e G
k m km
m km
m km
m

m=1

Pk =

e f B
e e G
+f e B
+f f G
k m km
k m km
k m km
k m km

m=1

Qk =

n
m=1

km

m=1

f f B
e e B
e f G
f e G
k m km
k m km
k m km
k m km

2. Segunda
P +JQ = V
k

Pk + J Qk =
Pk + J Qk =

m=1
n

Qk =

km

m=1

Vk eJk Ykm e

m=1
n |
m=1

km

Jm

Vm e

V Y
V |
k km m k
m
km

|V Y
V |Cos(
)
k km m
k
m
km

Pk =

n
k

m=1

|V Y
V |Sen(
)
k km m
k
m
km

3. Tercera
P +JQ = V
k

Pk + J Qk =
Pk + J Qk =

n
k

km

m=1

Vk eJk Ykm e

m=1

km

Vm e

Jm

JB )
V V (Cos
+ JSen
)(G
k m
km
km
km
km

m=1
n

Pk = Vk

Vm (Gkm Coskm + Bkm Senkm )

m=1

Q =V n
k
k

m=1

5.2
Tareas
formulacion

B
V (G
Sen
Cos
)
m km
km
km
km

basicas

en

la

Se denen dos tareas basicas en la formulacion del problema de ujo e carga, la primera
se reere al planteamiento de modelo matematico y la segunda a la seleccion del metodo para
la solucion de dicho modelo.
En la formulacion del modelo matematico las ecuaciones nodales se plantean en funcion
de las cuatro variables nodales y los parametros del sistema. Estas ecuaciones siguen la
siguiente extructura:
F (P , Q , |V i|, i) =
Ni Ni
0
Por cada nodo son expresadas dos ecuaciones como se presenta a seguir:
P =
i

N
n=1

|Y

V V |Cos( )
in i n
i
n
in

Q = N |Y V V |Sen( )
i
in i n
i
n
in
n=1

El modelo a solucionar es escrito de la siguiente forma:

Pi = (P iEspecicado

PCalculado ) = 0;
i

Qi = (QEspecicado
i

Q
Calculado

)=0

En la seleccion del metodo numerico se debera tener en cuenta el tipo de modelo a

resolver (lineal o nolineal), el tipo de aplicacion que se desee resolver, entre las que se cuentan
(operacion, planeacion, diseno), la calidad de respuesta esperada y los requerimientos de
memoria.

5.3

Resena historica

. Ecuaciones que describen las relaciones nodales son:

IN = YN VN
VN = ZN IN
Para resolver este sistema de ecuaciones es necesario utilizar el computador, sin embargo
este solo se aplica a partir de la decada de los 60.
Inicialmente para estudiar las redes electricas fue empleado el analizador de redes, con
este se reproduca la red a pequena escala en el laboratorio, la restriccion es que solo era
posible reproducir sistemas pequenos.
El desarrollo del computador permite que las ecuaciones de red y sus metodos de soluci
on, adquieran gran importancia.
- Inicialmente fueron desarrollados los metodos directos: IN = YN VN
Problemas: Lento y en muchos casos no convergen
- Luego aparecio: VN = ZN IN
Con este se mejoro la convergencia.
Problemas: Exigencia de memoria y lento
- Posteriormente se desenvuelve el Gauss.
Requiere poca memoria.
Problemas: Lento y diverge facilmente por malos condicionamientos.
- Mas adelante aparece el metodo Newton-Raphson, mostrando poderosas propiedades
de convergencia, pero con excesivos requerimientos de memoria, haciendolo no
competitivo en ese momento, debido al manejo de matrices de gran tamano con
computadores lentos y de poca capacidad.

Tinney y Walker presentaron la solucion a dicho problema, con la aplicacion de las

ma- trices dispersas, al estudio de las redes electricas.
- Posteriormente aparece el ujo de carga desacoplado,desarrollado por B. Stott, mejorando
en rapidez y requerimiento de memoria. Mas adelante aparece el ujo de carga
linealizado y el fcdc

5.4

Metodos de solucion

Los mas estudiados por los especialistas del tema son:

Iterativos ( Lineales y no lineales)

Gauss
Guss - Seidel

Acoplado

No-Lineales (Linealizados)

:
.Metodo de Newton Raphson
de Newton
Desacoplados : .Desacoplado
.Desacoplado rapido
de Newton
.Ward y Hale

N oLineales

5.5

Acoplados :
.Metodo de la No linealidad de segundo orden
Desacoplados : .Metodo de la No linealidad de segundo orden desacoplado

Condiciones y Definicion del problema

El sistema se estudia en condiciones estacionarias. La red opera de manera balanceada.

Cantidades o variables asociadas a cada nodo

G , , QQ
G
PP

|V D|, D

V ariables nodales

De las seis variables nodales mencionadas anteriormente, dos son conocidas, estas son la
demanda de potencia activa y reactiva: PD y QD
As las ecuaciones nodales se plantean en funcion de las siguientes cuatro variables nodales

Potencia inyectada (S S )
G
D
.activa

.reactiva
P G , QG
|V |, V oltaje
.magnitud

.angulo

5.6

Clasificacion de las variables

- Variables incontrolables: PD y QD
- Variables de control: PG y QG
- Variables de estado: |V | y

5.7

Lmites practicos de las variables

De estado
|Vi |min |Vi | |Vi |max
Lmite mnimo por estabilidad y el maximo por
aislamiento.
Capacidad de transmision | | max | |
i

De control
P
P
Gmin

5.8

P
Gi

;Q

Gmax

Gmin

Q
Gi

Gmax

Tipos de Nodos (Clasificacion de los Nodos)

Nodo Flotante, Slack, de holgura o de referencia (V, )

E ste nodo es el encargado de establecer el balance global del sistema, es decir,

asume los sobrantes (faltantes) de carga de los otros nodos, mas las perdidas del
sistema. Es
planteado de la siguiente manera:
n

i=1

P n P = P
G
D
p
erdidas
i=1

Generalmente se selecciona el nodo de mayor generacion. Este nodo generador

debera
tener la suciente cantidad de reactivos para el control de la magnitud del voltaje.
Se supone conocida la potencia generada en todos los nodos generadores, excepto en el
nodo denominado slack, que corresponde al nodo generador de mayor capacidad; adema
s las demandas son conocidas.
n

Potencia generada conocida PkG =

PG
i=1
i = slack

Potencia demandada PkD =

PD

i=1

Al reescribir las ecuaciones anteriores, el balance global del sistema cumple con la siguiente
ecuacion:
P
+P P =P
Gslack
kG
kD
p
erdidas
De la ecuacion anterior, el nodo slack debe generar una potencia para garantizar el
balance global del sistema:
P
= P + P + P
Gslack
kG
kD
p
edidas
Por lo tanto la potencia generada en el nodo slack debera ser la suciente para cubrir
la
diferencia entre (PkG + PkD ) mas el valor de las perdidas. Las perdidas son una
variable
no conocidas, por lo tanto, la variable PGslack tendra que ser libre para que se
pueda cumplir la ecuacion anterior.
Al ser la magnitud de voltaje una variable especicada, el QGslack sera una variable libre
(no especicada). El nodo slack entonces se encargara de establecer el balance global
del sistema, tanto en potencia activa como en potencia reactiva, cubriendo los faltantes
(sobrantes) que presente el sistema.
En este nodo las variable son clasicadas as:
|V |, Variables especicadas
P, Q Variables desconocidas.
El angulo se asume a 0o y sera referencia para los demas nodos del sistema.

Nodo de Voltaje controlado (PV).

Generalmente se asumen los nodos con generacion representativa.
Las variable son clasicadas as
: P, V Variables
especicadas Q, Variables
desconocidas
El lmite de reactivos podra ser considerado, dependiendo de la aplicacion.
Q

Q
imin
i
Q
imax

A este tipo de barra pertenecen aquellos nodos que tengan elementos con capacidad para
controlar la magnitud del voltaje.
Barras candidatas a ser denominadas de voltaje controlado:
Generadores
Compensadores sincronicos
Compensadores estaticos activos (controlados por tiristores)

 Transformadores con combinadores automaticos de Taps bajo carga, siempre y

cuando su operacion este en el rango de control.

 Nodo de Carga (P,Q)

Corresponde a los nodos de carga o nodos donde la generacion es muy pequena. Esto
es, representa las barras donde la carga predomina sobre la generacion. Ejemplo al
mismo barraje se conecta carga y generacion, siendo este ultimo poco representativo
frente a la magnitud de la carga. Tambien representa las barras de paso.
Las variable son clasicadas as
:
P, Q Variables especicadas
V, Variables desconocidas
Una Forma alternativa de denir el nodo de carga considerando restriccion de voltaje
es la siguiente:
Barras de carga (P-Q) donde la magnitud de voltaje se requiere que este dentro de unos
lmites
cantidad prejados: Vmin V Vmax . Este tipo de barra permite establecer la
de potencia reactiva que se necesita inyectar en la red a n de mantener una operacio
n adecuada.
En el problema de ujo de potencia este tipo de nodo es operado de la siguiente manera:
Cuando uno de los lmites de voltaje barra es sobrepasado, el nodo cambia de tipo a
ser (P - V),Q queda entonces como una incognita. Resuelto el ujo se calcula la
diferen Q
cia (Q
), que indica la cantidad de reactivos (L-C) que
especif icado
calculado
deberan ser
instalados.

5.9

Metodo de Gauss y Gauss Seidel utilizado en el

pro- blema de Flujo de Potencia

Tuvieron gran acogida en los anos 60 por su facilidad en la implementacion (mnimos

requer- imientos de memoria y facilidad en su programacion).
Ecuaciones planteadas:
P J Q =n V
k
k
Y
V
km m

m=1

Pk J Qk = V Ykk Vk + V
k

Ykm Vm

k
m=1
m= k

P JQ
k

= Vk +
Ykk

m= k
m=1

V
km m

De la ecuacion anterior, se despeja el voltaje del nodo k:

Vk = 1
Y kk

Pk

Qk

m=1
m= k

Ykm Vm

Ademas en los nodos tipo voltaje controlado se calcula:

Q =
k

V
km m

; Siendo (Qk = QGk QDk )

I
mag

m=1

Q =V n
k
k

m=1

(G

km

B Cos
Sen
)
km
km
km

En estos nodos se debera satisfacer que:

Q
Gk min
Gk
Q
Gk max
El crtero de parada:
i
k
max.|v k+1
v i | Pare

Por facilidad se acostumbra suponer el vector de voltajes iniciales en los nodos del sistema
de la siguiente manera:
- En barra (P-Q): 1 0o
- En barra (P-V): |Vespecif icado| 0o
| 0
- En barra (V-): |V
o
slack

En la barra de referencia se conoce la magnitud y el angulo del voltaje, por lo tanto se

tienen n-1 incognitas, as que en el calculo de los voltajes nodales no se tendra encuenta la
ecuacion de este nodo.

ALGORITMO DE GAUSS SEIDEL

Parametros de la red (Lneas y

transformadores)

Informacion Nodal
v = 0 (contador de iteraciones)
Condiciones iniciales (Voltajes)

v = v + de
1 nodos)
k = 1 (contador

si

Es k el nodo slack

k=k+1
Calcular
Es k nodo Generador

si
no

Q = I n (V
k
m
Y
V )
km m

no

k=k+1

Calcular

m=1

Q Q

Q
Y1

Vv =

kk

kmin

P jQ
k

n
k

Ykm Vm

m=1 (m=k)

kmax

kmax

Lmite violado
=Q
o Q =
k
kmin
k
V

Se calcularon voltajes en todos los nodos ?

Normalizar voltajes en nodos Generadores

l=0
l = l +(Especicado)
1

V v = V lv
o
l
l
l
l nodos de Generaci
on

v+1

v+1

v
= |V k V k |
|
|

Fin del proceso (resultados complementarios)

Chequeo de convergencia M

ax V k

|
|V v|

Ejercicio: Dado el siguiente sistema de energa electrica de 3 nodos, determinar los

valores de voltajes nodales, ujos de potencia por las lneas y perdidas. Emplear para su
solucion el metodo de Gauss.
2

Figura 5.2: Red electrica de 3 nodos

Lnea

p.u.

Y /2
p.u.
1 2 0.01
1 3 0.03
2 3 0.04

p.u.

0.06 0.005
0.20 0.01
0.25 0.015

= 100 MVA

base

Datos nodales
N odo

1
2
3

PG
QG
PD
QD
V oltaje
(MW ) (MVAR) (M W ) (MVAR) (p.u.)
150
30

80
15

Voltaje nominal del sistema: 115 KV

60
20
45

30
10
15

1.02
1.01

5.10

Metodo de Newton Raphson en la solucion del

pro- blema de Flujo de Potencia

El metodo de Newton Raphson es del tipo iterativo que linealiza las funciones en el punto de
trabajo. Metodo recomendado en la solucion de ecuaciones algebraicas no lineales, como
es el caso de las ecuaciones nodales de potencia inyectada que describen las redes el
ectricas en estado estacionario.

5.10.1
atico

Modelo

matem

Las ecuaciones nodales describen el comportamiento del sistema en estado estacionario. En

cada nodo se plantea una ecuacion como se describe a seguir:
1

S i1

SGi

S i2

S ij

SDi
Sistema externo

Sistema interno

Figura 5.3: Planteamiento de la ecuacion de potencia nodal inyectada

SN i = PN i + J QN i = SGi SDi (Sistema externo)
SN i = V i

V (Sistema interno)

im m

m=1

En estado estacionario la potencia inyectada por el sistema externo debera ser igual a la
inyectada por el sistema interno, se debe cumplir entonces que:
S
Ni

= V eJ i
i

m=1

eJim V eJ m ( Formula que establece el balance nodal)

im

Con base en el tipo de nodo, las ecuaciones son planteadas como ecuaciones del subsistema 1
y subsistema 2. Las ecuaciones del subsistema 1 establecen el modelo matematico del
problema de ujo de carga y son escritas como sigue:
Pi = P esp P cal =
i
i
(PG
i

PD ) Vi
i

n
m=1

Vm (Gim Cosim + Bim Senim ) = 0

n
cal
G
i
Qi = Qesp
Qi = (Q i QD ) Vim=1 Vm (Gim Senim Bim Cosim ) = 0
i

5.10.2

Metodo de solucion de Newton

Este metodo se basa en la expansion en series de taylor y es descrito as:

F (X k + X k ) = F (X k ) + F (X k )X k + ..... = 0
Se desprecian potencias superiores a 1 y se despeja el k-esimo incremento de la variable
de estado X.
X k = F (X k ) F (X k )
1

X k+1 = X k + X k
Si se tiene un modelo de n ecuaciones con n incognitas, el problema linealizado toma la
siguiente forma:

F1 (x)
F (x)

F (x)
n

F1

X1
F
2

21
X

F1

X2

F2
X2

F1
X
Fn
2

Xn

.
.

. .

X
n
F
1

X
Fn
2

X1

n
X

n
n

La notacion del conjunto de ecuaciones anteriores es: F (x) = [ J ] X

n

ik
Observaciones a la ecuacion de potencia nodal: Pi + J Qi = Vi
Y kV
k=1
Las potencias netas inyectadas dependen de su propia tension,
as como tambien

de
las tensiones de los otros nodos con los cuales dicho nodo esta interconectado y de
los
parametros de la red.

Si los nodos i y k no estan interconectados no aparece termino Yik en la matriz Ybus por
lo tanto el termino (i k) no incide en la sumatoria de la potencia nodal.
En los sistemas de transmision R
X, ademas la diferencia angular es pequena.
Las ecuaciones son sinosoidales y de segundo orden.
Estructura de la matriz Jacobiana:

X
.

Dado el siguiente sistema de tres nodos, escribir el modelo matematico linealizado del
prob- lema para resolverlo por el metodo de Newton Raphson.

Figura 5.4: Red electrica de 3 nodos

P1
P

P1
1
P
2

P
Q31 =
Q2
Q

P1
2
P
22
P
23
Q
21
Q
22
Q
3
2

Q1
1
Q
2
1

P1
V1
P
V 2

P2
3
P
3

Q1
V1
Q
V 2

Q2
3
Q
3

P1
V2
P

P2
V3
P
V 3

V22
P
V23
Q

V21
Q
V22
Q
V 3
2

Q2
V3
Q
V 3

2
V3
V1
V2
3

Como se observa de la gura anterior, las ecuaciones nodales estan descritas en funcio
n de (P, Q, V y ). El modelo matematico del problema de ujo de carga toma la
siguiente estructura matricial.

NPQ
+

NPV
V

NPQ

=
Q

Los errores de potencia (P y Q) son calculados como la diferencia entre la potencia

inyectada por el sistema externo (especicado) e interno (calculado) de la siguiente manera:
P = P

especif icado

calculado
Q = Q

especif icado

calculado
Q
El modelo de n ecuaciones linealizadas para el problema de ujo de carga es escrito de la
siguiente manera:

V
V

Las ecuaciones podran ser planteadas en forma polar o rectangular, dependiendo de

como se describa la variable de estado (voltaje).
Vp = ep + J fp ; Vp = |Vp |eJP
Ecuaciones del Modelo de flujo de Carga de Newton:
En el sistema de energa electrica podran ser planteadas 2.n ecuaciones nodales (n de
poten- cia activa y n de potencia reactiva), que conforman el modelo matematico y que se
denominan sistema de ecuaciones 1 y 2.
Las ecuaciones del sistema 1 son aquellas que conforman el modelo del ujo de carga, el
resto de ecuaciones (2.n - ecuaciones del Sistema 1) se denominan ecuaciones del sistema 2.
Las ecuaciones que conforman el sistema 1 dependen del tipo de nodo, as: En
nodos generadores se plantea ecuacion de potencia activa y en el nodo de carga se plantea
tanto la ecuacion de potencia activa como de reactiva.
Con las ecuaciones del sistema 2 se podran determinar resultados complementarios, como
son la potencia activa y reactiva generada en el nodo Slack y la potencia reactiva generada en
los nodos generadores.
Terminos de la matriz jacobiana:
Esta matriz esta compuesta por terminos diagonales y no diagonales. Los terminos
diago- nales escritos en forma polar son:
Pk Pk Qk Qk
, V , , V
k

Los terminos fuera de la diagonal en forma polar, son descritos as:

Pk Pk Qk Qk
, V , , V
La descripcion polar es la mas frecuentemente utilizada en la literatura especializada.
El ujo de carga es la etapa inicial de otros problemas, como es el caso del estimador de
estado, ujo de carga desacoplado, entre otros y en la mayora de estos se parte del hecho
de que el modelo se plantea en componentes polares.
sica planteada en componentes polares:
Ecuacion b
n Y
P + JQ = V
eJkm V eJm
k
k
k
km
m
eJk
m=1
Deduccion de los terminos de la diagonal: (

Qk

Pk + J Qk = V 2 Ykk e
e

Jk

kk

+ Vk

n
m=1
m= k

Jkm Jm Jm
Ykm e k
V
e

Pk

Qk =Jk
J Vk
+J
e

Ykm

km

Jm

Vme

m
m

J
eJe k e eJ
kk kmJ m )
Y km
Q
+J
+
J
k =
Vm
(Vk
J

V2
Ykk

m=
1

!
k

Pk +

J Qk

eJ +k J
= kk Q )
J
V +J
2 Yk(Pk
k

k
k

) +J
= + Qk )
J
J J (P
V 2 Bk
(G kk
kk

Separando en parte
real e imaginaria se
obtiene:
Pk

= Qk

V2
Bk
k

Qk = + Pk
k
V
2
G
kk

Deduccion de los
P la
terminos de
Q
k y
diagonal:
( k VV
k k
V ) Vk
k

Jk

= Yke
2 k J
+k
kk
V
e
+

J Jm

Ye

k km
m
Vm

m=1
m= k

=
nY
J
2V e e
e e
Q
J mk
V+ Jk 2 Yk Jkk 1= mJkm
Vk
m
m
Vk + Vkk =
k
Vm
k
Pk
V V

=
n
2 e e e eJ
V
Q
J
+ Jk YkkJkk Y Jkm
Vk
km
V + Vkk
m Vm
k
m
=1

Pk
V V
k

Pk
V V
k

=
+ J Qk

V
Q2 e
+ Jk YkJkk
Vk
Vk +
k
Pk

=
) + J Qk

+
QV
J
2
k
+ J (G B Pk
Vk kk
k
V
k
k

Partiendo en parte real e

imaginaria se tiene:
Pk
= + Pk
V V V
k
2
G
kk

Q
k

V = V 2 B

+Q
Vkk

kk

fuera de la diagonal: ( Pk m

Qk

eJkm
n Y
P
V ekm
Jm
k

+m=1

J
Q
k

=
V
k

e
J

= J Vk J Jm
e
e
eJk Ykm
km
Vm

= + J fk
)
J (e J fm )
P
J )(Gkm B
km m
+(ek
J

Q
k

m
Partiendo en parte real e imaginaria se
obtiene:
Pk = e
f +e
f = Hkm
m mf m m
m e B e G k G f B
k k k
k k k
k
m
m
m
m

Qk

= ek

e k +f
e
f = Jkm
m ek
m
m fm fm
k
k
G
B
B
G
k
k
k
m m m

k
P
Deduccion de los terminos
k
Q
fuera
de
la
diagonal:
(
VVm y

k m
V )
V
m

Pk
Vm

m
= J
e
QV
k k
J
+ k
J V Ykme
km

J
e
m

V
B

L
k
k

=
2

J VQk = Y J V Jm
k k
Vm V V
k
e m me
Jk m
e

Partiendo en parte real e imaginaria

se obtiene:

e f B +
Gkm
k
k m km
V= ek emP
fk em Bkm + fk fm Gkm =
Nkm
V

V= fk em Gkm fk fm Bkm
ek em Bkm ek fm Gkm =

Lkm
Terminos
de la
diagonal:

H
k
k

Q
k

Nkk = Pk + Gkk V 2
J =P G V2
kk

kk

Terminos fuera de la diagonal:

Hkm = Lkm = ek em Bkm ek Fm Gkm + Fk em Gkm Fk Fm Bkm
N = J = e e G e f B + f e B + f f G
km

km

k m

km

k m

km

k m

km

k m

km

Ejercicio: Dada la siguiente red de 3 nodos y 3 lneas, una demanda de 720 MW y

410
MVARS y una generacion de 800 MW, solucionarlo por el metodo de Newton Raphson .
G1

D1

L2

L3

D2

D3
G3

Figura 5.5: Red electrica de 3 nodos

Voltaje de operacion = 230 KV
Pbase = 100 MVA
Datos de nodo

N odo
N odo
N odo
Datos de lneas y transformadores
Nodo A
1
1
2

METODO

GAUSS SEIDEL

VOLTAJES NODALES

1
2
3

0.82530
1.02

ANGULOS [Grados]

1
2
3

-14.08927
-0.45146

GENERACIO
N[MW]
POTENCIA
ACTIVA

1
2
3

0.97214
400.17283

POTENCIA REACTIVA[Mvar]

1
2
3

0.91405
292.22033

DEMANDA
POTENCIA ACTIVA [MW]

1
2
3

300
220

FLUJOS DE POTENCIA
POTENCIA ACTIVA [MW]

1-2
1-3
2-3

8.92751
-179.09535

POTENCIA REACTIVA [Mvar]

1-2
1-3
2-3

-9.78492
-95.814

FLUJOS DE POTENCIA
POTENCIA ACTIVA [MW]

2-1
3-1
3-2

-8.91736
-189.09019

POTENCIA REACTIVA [Mvar]

2-1
3-1
3-2

-7.02733
149.4276

PERDIDAS EN LAS LINEAS

POTENCIA ACTIVA [MW]

1-2
1-3
2-3

0.01016
9.99484

POTENCIA REACTIVA [Mvar]

1-2
1-3
2-3

-16.99226
53.61364


ALGORITMO PARA LA SOLUCIO N DEL PROBLEMA
DE FLUJO DE POTENCIA POR NEWTON RAPHSON

Parametros de la red (Lneas y

transformadores)

Informacion Nodal
Condiciones iniciales (Voltajes), para

metros
k = 0 (contador de iteraciones)

. Calcular: P (i = 1

Q (i = 1,nd)

,nd)
i

. Comprobar que max |Pi |, |Qi | Solucion

H N
J L

P
Q

V /V

Resolver el modelo linealizado

=

V /V

H N
J L

P
Q

. Actualizar las variables

V k+1 = V k (1 +

k+1

V k

= +

k = k + 1

Ejemplo: Dado un sistema de 4 nodos y 4 lneas, como se muestra en la gura 5.6,

solu- cionarlo por el metodo de Newton Raphson.

Figura 5.6: Red electrica de 4 nodos

Datos de lneas:
Lnea
L
L1
L2
L3

Nodo inicial y nal

Y /2
sh en

Serie en p.u.

Datos de barras:
Barra
1
2
3
4
Resultados:
Inicio del Subproblema 1

8.9852 44.8360i 3.8156 + 19.0781i 5.1696 + 25.8478i

0
3.8156 + 19.0781i 8.9852 44.8360i

0
5.1696 + 25.8478i

Y
5.1696 + 25.8478i

bus

Condiciones iniciales: Voltajes incognitas = 1.0 p.u.

Angulos incognitas = 0o
Primera iteracion
Errores de potencia [P |Q] =
Cumple tolerancia ?

J acobiano =

1.5966 1.9395 2.2129 0.4465 0.8345

[NO]

26 .3648 15.4209

3.0842
41.7857
45.4429 5.2730
9.0886

Resolviendo para y V , y recalculando y V tenemos:

Segunda iteracion
Errores de potencia [P |Q] =
Cumple tolerancia ?

0.0323 0.0645 0.0359 0.0342 0.0620

[NO]

44.3749

J
acobiano
26.1256 =15.1217 41.2473

10.3562

4.1296 2.1828

Resolviendo para y V , y recalculando y V tenemos:

Tercera iteracion:
Errores de potencia [P |Q] =

0.0339 0.1455 0.0424 0.0414 0.1646 .10

3

Cumple tolerancia ?

[SI]

J acobiano

44.3270

26.1026 15.0938 41.1963

10.33721

Por lo tanto, la solucion nal para las variables de estado V y son

|V |
Anguloso

Resultados con base en ecuaciones del subproblema 2

Para el calculo se aplica una potencia base 100 MVA
Potencias activas inyectadas:
P
P12
P3
P

=
=
=
=

136.7948 Mw
169.9966 Mw
199.9854 Mw
237.9958 Mw

Potencias Reactivas inyectadas:

Q1
Q2
Q3
Q

=
=
=
=

83.4977 MVAR
105.3459 MVAR
123.9235 MVAR
131.8393 MVAR

Flujos de Potencia:

P
P12
P13
P24
P34
P21
P31
P42
43

V1
= 1.0000
=
0

4.1183 2.1655

5.10.3

Transformadores reguladores

Los transformadores reguladores se pueden usar para controlar los ujos de potencia activa y
reactiva en un circuito. El modelo de este transformador para estudios de ujos de carga fue
explicado en 3.3.12. (transformador en fase - control de potencia reactiva)
En la siguiente gura se presenta una representacion detallada de este transformador.
j

1: t

t Vi

Vi
-

VJ

Figura 5.7: Representacion del transformador regulador

El modelo equivalente del transformador mostrado anteriormente es:
ty

Vi

t(t-1) y

Vj

(1-t)y
-

Figura 5.8: Modelo del transformador regulador

t es representado por un real, ya que no existe desplazamiento angular entre el primario y
secundario, por lo tanto solo existe cambio en la magnitud.
En el problema de ujo de potencia los taps podran ser includos de manera manual
o automatica.
Manejo manual de Taps en transformadores
En el caso manual la variable t no es includa explicitamente en la formulacion del
Newton Raphson. En este despues de cada iteracion los Taps deben ser ajustados usando
una tecnica de desplazamiento asi:

tnuevo

viej
o

reg

=t
+ ci (V
iiteracion)
i
i

V ). (t siempre va atras de V en una

i

= Voltaje que se
V quiere mantener con el Tap.

c = [ 0.5
1.0] valores tipicos empleados.
i

En cada iteracion se calcula un nuevo valor de t el cual se

reemplaza en el modelo del transformador para posteriormente
incluirlo en la matriz Ybus , la que posteriormente se emplea en el c
alculo de las potencias inyectadas y de la matriz jacobiana.
Manejo del Taps automatico
transformador bajo carga

del

El taps podra ser manejado automaticamente en el ujo de carga,

para esto en el nodo en el que se controla el voltaje con los taps del
transformador se efectua un cambio de variable, voltaje nodal se
cambia por la variable taps. En la siguiente gura se presenta un
sistema de cinco nodos que contiene 2 transformadores reguladores.

Figura 5.9: Red electrica de 5 nodos con

transformadores con taps
El modelo matematico del ujo de carga includos los dos
transformadores de la gura anterior se presenta a seguir.
P1
P
P2
P3
P4
Q5
Q2
3

Los terminos de la matriz Jacobiana se calculan asi:

Pi Qi Pi Qi
t2 t , t2 t , t3 t , t3 t
2

5.10.4

Solucion del flujo de potencia mediante el metodo de

Newton
Raphson Desacoplado

Caracterstica propia de cualquier sistema de transmision de potencia es la fuerte interrelacio

n que existe entre:
P
Q |V |
Se tomara provecho de esta caracterstica para obtener un modelo de ujo de carga
que emplee menos memoria y tiempo de computo y que entregue resultados lo
sucientemente precisos.
En esta forma de trabajo, resuelve separadamente o desacoplalos problemas (P )
y |V |). Este desacople no incluye independencia absoluta de los subproblemas.
(Q
Ecuaciones basicas y elementos de la Matriz Funcional
Caractersticas de los sistemas de transmision.

Alta relacion reactancia/resistencia

B

ik

Susceptancia

>G

ik

C onductancia

La diferencia angular ik entre nodos adyacentes es pequena

Senik

ik Cosik

1.0

Considerando estas expresiones en el calculo de los elementos de la matriz Jacobiana en

coordenadas polares y utilizando las aproximaciones anteriores se tiene:

P
Q

H N
J L

V /V

Caractersticas de la submatriz N
Las conductancias Gik estan asociadas a valores de coseno.

 Las susceptancias Bik estan asociadas a valores de seno.

Se obtiene entonces valores relativamente pequenos.

En la submatriz J.
Ocurre lo mismo que en la submatriz N

En la submatriz H y L
El elemento Bik (susceptancia) esta acompanada de valores coseno
Se obtienen entonces valores relativamente grandes
Por lo tanto el peso asociado de las submatrices H y L es mayor que el de las submatrices
J y N, lo cual permite despreciar el efecto de estas u
ltimas.

P
QN

H O
O L

V /V

Se obtienen entonces dos ecuaciones desacopladas de un menor orden (casi la mitad de las
utilizadas en los metodos acoplados).
Caractersticas de las matrices H y L:

Tiene la misma forma de la matriz admitancia nodal

Son matrices simetricas en posicion no en valor
Son funcion de los valores de voltaje por lo cual estan cambiando constantemente
Tamano de la matriz H (NPQ + NPV)
Tamano de la matriz L (NPQ)


ALGORITMO PARA LA SOLUCIO N DEL PROBLEMA DE
FLUJO DE POTENCIA DESACOPLADO DE NEWTON

Parametros de la red (Lneas y

transformadores)

Informacion Nodal
Condiciones iniciales (Voltajes), para

metros
k = 0 (contador de iteraciones (P-))
k 1 = 0 (contador de iteraciones (Q-V))
2

. Calcular P = P P P (V, Y )
N

Comprobar

P
< Comprobar

Q <
N

PN = H k k

k+1

= + k

k =k +1
1

. Calcular QN = QG QD QN (V Y )

Comprobar

Q < Comprobar
N

Solucion

P <
N

N = L
QV

k
Vk

V k+1 = V k + V
k

k =k +1
2

Ejemplo: Resolver el sistema de dos barras mostrado en la gura 5.10 por el metodo
ujo de carga desacoplado de Newton.

k Nodo k: tipo V,

i km

m Nodo m: tipo P,Q

Z km = R km + JXkm

Vk k

i mk
Vm m

y/2 = Jb sh

y/2 = J b sh

km

km

Figura 5.10: Red electrica de 2 nodos

Datos del sistema
Barra 1:

(V-)

Entonces: V1 = 1.0 p.u 1 = 0.0 : incognitas : P1 , Q1

Barra 2: (P-Q)
Entonces: P2 = 0.30 p.u. Q2 = 0.07 p.u. incognitas : V2 , 2
Tolerancia de convergencia: = 0.002 p.u.
Parametros del sistema (lnea de transmision):
R1

= 0.2 p.u. , X1

= 1.0 p.u. , y/2 = 0.02 p.u.

(Pbase = 100 MW)

Expresiones requeridas en el calculo
P2 = 0, 1923V 2 0, 1923V2Cos2 + 0, 9615V2Sen2
Q = 0, 9415V 2 0, 1923V Sen2 0, 9615V Cos
2
2
2
2
2
P = 0, 30 P
2
2
Q = 0, 07 Q
2
2
H = Q + 0, 9415V 2
22

L22 =

2
1

2
2

[Q2 + 0, 9415V V]

Proceso iterativo de solucion

Condiciones iniciales
i) Kp = Kq = 1; p = q = 0; 2 = 0o ; V o = 1.0; = 0.003

Primera iteracion (P-)

ii) P2 (V 2o , 2o ) = 0, 1923(V2 o )2 0, 1923V2 o Cos2o + 0, 9615V
2 o Sen
2o
P2 (V 2o , 2o ) = 0, 1923(1)2 0, 1923(1)(1) + 0, 9615(1)(0) = 0
2 o ,2 o ) =
P2 (V2 o , 2 o ) = 0, 30 0, 0 P2 (V
0, 30
iii) |P2 (V 2o , o2)| = 0, 30 < = 0, 003 [NO]
2 o )2
iv) H22 (V2o , 2o ) = Q2 (V2 o , 2o ) + 0, 9415(V
o
o 2
o
o
o
o
o
Q2 (V2 , 2 ) = 0, 9415(V
0,
)Sen
0,
Cos

2
2
2 9615V
2
2
)
1923(V
2
o
o
Q (V2 , 2 ) = 0, 9415 0, 9615 = 0, 02
H22 (V2o , 2o ) = (0, 02) + 0, 9415 = 0, 9615
2 (V o , o ) = H2 2(V o , o ) o
P
2 o =
= 0, 31201

0,30
2
2
2
2
2
0,9615
v) 21 = 2o + 2o = 0 0, 31201 2o 0, 3121
vi) p = 0 + 1 = 1
vii) Kq = 1 (x)

Primera iteracion (Q-V)

x) Q2 (V2o , 2 1 ) = 0, 9415(V2 o )2 0, 1923V
2 o Sen
2 1 0, 9615V
2 o Cos
2 1
Q2 (V2o , 2 1 ) = 0, 9415 0, 1923Sen(0, 31201) 0, 9615Cos(0, 31201)
o 1
Q2 (V2 , 2 ) = 0, 08545
2 o , 2 1 ) =
Q2 (V2o , 21 ) = 0, 07 0, 08545 Q2 (V
0, 01545
xi) |Q2 (V 2o , 21 )| = 0, 01545 < = 0, 003 [N O]
xii) L22 (V o , 1 ) = V 21o [Q2 (V2o , 12) + 0, 9415(V o )2 ] 2
2
2
22
o
L (V2 , 21 )

= 0, 08545 + 0, 9415 = 1, 026950,01545

2
o
1
Q (V2 , 2 ) = L22 (V2 o , 2 1 )V2 o V2 o = 1,02695= 0, 01505
V o = 0, 01505
2

xiii) V21 =
= V2 o + V2 o = 1, 0 0, 01505 = 0, 98495
V 1 0, 98495
2

xiv) q = 0 + 1 = 1
xv) Kp = 1

Segunda iteracion (P-)

P2 (V 21 , 21 ) = 0, 1923(0,
98495)2 0, 1923(0, 98495)Cos(0, 31201)+
0.9615(0, 98495)Sen(0, 3120)
P2 (V 21 , 21 ) = 0, 28442
0, 1923V Cos + 0, 9615V
(V 1 , (1 )0,=28442)
0, 1923(VP
1 )22 (V
1
2 1 ,2 1 ) = 1 0, 01558
P2 (V2 1 , 2 1 ) =ii)P
0,230
2
2 2
2
2
2
2
1
1
iii) |P (V 2 , 2 )| = 0, 01558 < = 0, 003 [NO]
2 1 )2
iv) H22 (V21 , 21 ) = Q2 (V2 1 , 21 ) + 0, 9415(V
1 )2
1
1 Cos21
2
2
2
Q2 (V 1 , 1 ) =2 0, 9415(V
0,2 1923V 0, 9615V
Q2 (V21 , 2 1 ) = 0, 98495[0, 9415(0, 98495)0, 1923Sen(0, 31201)0, 9615Cos(0, 31201)]
Q2 (V21 , 2 1 ) = 0, 07021
2 = 0, 84316
0,01558
H22 (V21 , 21 ) = 0, 07021 + 0, 9415(0, 98495)
2
P2 (V2 1 , 2 1 ) = H22 (V2 1 , 21 )2 1 1 = 0,84316= 0, 01848
v) 22 = 21 + 21 = 0, 31201 0, 01848 = 0, 33049
2 = 0, 33049 (rad)
2

vi) p = 1 + 1 = 2
vii) kq = 1 (x)

Segunda iteracion (Q-V)

x) Q2 (V21 , 2 1 ) = 0, 9415(V2 1 )2 0, 1923V
2 1 Sen
2 2 0, 9615V
2 1 Cos
2 2
Q2 (V 1 , 1 ) = (0, 98495) [0, 9415(0, 98495) 0, 1923Sen(0, 33049) 0, 9615Cos(0,
2
2
33049)]
2
1
Q (V2 , 2 1 ) = 0, 07906
Q2 (V21 , 21 ) = 0, 07 (0, 07906) Q2 (V2 1 , 21 ) = 0, 00906
xi) Q2 (V21 , 21 )| = 0, 00906 < = 0, 003 [NO]
xii) L22 (V 1 , 2 ) = 211 [Q2 (V21 , 2 2) + 0, 9415(V 1 )2 ] 2
2
2
V
1
21 22
L22(V , ) = 0,98495
[0, 0706 + 0, 9415(0, 98495)2] = 1, 00760
0,00906
Q2 (V21
, 22 ) = L22 (V2 1 , 22 )V2 1 V2 1 = 1,0076
= 0, 00899
V 1 = 0, 00899
2

xiii) V22 =
= V2 1 + V2 1 = 0, 98495 0, 00899 = 0, 97596
V 2 0.97596
2

xiv) q = 1 + 1 = 2
xv) Kp = 1 (ii)

Tercera iteracion (P )

P2 (V 22 , 22 ) = 0,
29887
P2 (V 2 , 2 ) = 0, 30 (0, 29887) = 0, 00113
2

iii) |P2 (V 22 , 22 )| = 0, 00113 < =

0, 003 [SI]
El subproblema
ir al paso
(viii)
0,
ii) P P(V 2convergio
, 2 ) = 0, 1923(V
1923V 2 Cos2 + 0, 9615V
2 )2

2
2
2
2
2
2
viii) k = 0 (ix)
p
ix) El subproblema QV aun no converge, pues Kq = 1, entonces ir a (x)

Tercera iteracion (Q-V)

x) Q2 (V22 , 2 2 ) = 0, 9415(V2 2 )2 0, 1923V
2 2 Sen
22
2 2
0, 9615Cos
2
2
Q2 (V2 , 2 ) = 0, 07074
2 2 ,2 2 ) =
Q2 = (V22 , 22 ) = 0, 07 0, 07074 Q2 (V
0, 00074
xi) |Q2 (V 22 , 22 )| = 0, 00074 < = 0, 003 [SI]
Entonces el subproblema QV convergio (xvi)
xvi) K = 0 (xvii)
q
xvii) El subproblema P tambien convergio, ya que Kp =
0 xviii) Convergencia
lograda:
2 = 2 = 0, 33049 (rad.)
V = V 2 = 0, 97596 (p.u.)
2

2
2

Observacion:
Si fuera efectuada una iteracion adicional, los valores obtenidos de 2 y V2 son los
siguientes:
= 0, 33186 (rad.)
2

V2 = 0, 97522 (p.u.)

5.10.5
apido

Newton Raphson Desacoplado r

Ademas de las simplicaciones hechas al Newton Raphson desacoplado, se contemplan algunas

simplicaciones adicionales con el n de mejorar el metodo en lo que se reere a requerimiento
de memoria y tiempo de computo sin desmejorar la calidad de los resultados.
Las
simplicaciones al modelo matematico son hechas con base en las caractersticas que
presentan los sistemas de transmision, tal como se describe a continuacion:

Alta relacion reactancia/resistencia por lo tanto Bik > Gik

La diferencia
angular entre nodos adyacentes es pequena, por lo tanto:
Sen( ) 0
Cos( ) 1.0
k

 De experiencia obtenida en el estudio de sistemas de potencia, utilizando valores en por

unidad, se sabe que
Qk << Bkk V 2

En sistemas bien condicionados, se puede asumir que los valores de las tensiones de los
nodos del sistema son proximas a 1.
Si se aplican estas simplicaciones a las submatrices H y L se obtiene:
Hkk = Qk V 2 Bkk
k
H = V V (G Sen B Cos )
km
k m
km
km
km
km
L =Q V B
2

kk

kk

= V V (G

km

Sen

km

km

Cos

km

km

Sea V una matriz diagonal cuyos elementos no nulos son las magnitudes de las tensiones de
barra del sistema, se obtiene entonces:
H =VH
L=VL
Quedando los terminos de la matriz jacobiana de la siguiente forma:
= Vm (Gkm SenH
Bkm Coskm )
km

Vk Bkk
H kk = Qk
Vk
= G Sen B Cos
km
km
km
B km L km
=
k
Q
kk
2
L
km

kk

As las ecuaciones del metodo de Newton desacoplado asume la forma de:

P = H
V
Q
V

= L V

Ahora introducimos las siguientes aproximaciones:

Coskm 1
|Bkm| |Gkm Senkm |
|Bkk V 2 | |Qk |

La tercera aproximacion es verdadera, ya que las reactancias shunt (cargas, reactores,

capaci- tores y shunt de lneas) son mucho mayores que las reactancias serie de las lneas y
transfor- madores. Con las aproximaciones, los terminos anteriores toman la siguiente forma:
H

kk

V B
Bk kk

km

B km

H
L

L
kk
kk
Considerando adicionalmente que V V 1 se tendra:
m
k
H B
L B
Siendo la estructura de las matrices B y B semejantes a la de la matriz Ybus , con las
siguientes diferencias:

En B
En B

no aparece la lnea y la columna de la barra V

no aparecen las lneas y columnas de las barras V y PV

Las ecuaciones del metodo desacoplado rapido son las siguientes:

P = B
V
V

= B V

Donde las submatrices B y B son constantes y dependen solamente de los parametros

de la red.
Finalmente el metodo desacoplado rapido presenta un mejor desempeno cuando son
despre1
km y tenemos que:
km
ciadas las resistencias serie en la formacion de B , entonces b =
X
B
= 1

B km =
kk

B
kk

Xkm
mk

1
X
km

= Bkm k
= BkkB

Donde Bkm y Bkk son elementos de la matriz de susceptancias B y Xkm es la reactancia

serie de una lnea de transmision.

45


ALGORITMO PARA LA SOLUCIO N DEL PROBLEMA
DE FLUJO DE POTENCIA RAPIDO DE NEWTON
Parametros de la red (Lneas y
transformadores)
Informacion Nodal

Condiciones iniciales (Voltajes), para

metros
k1 = 0 (contador de iteraciones (P-))
k2 = 0 (contador de iteraciones (Q-V))
Conforman las matrices [B ] y [B ]
. Calcular P = P P P (V, Y )
N

Comprobar

P
< Comprobar

Q <
N

N=
P
k
V

[B ]k

k+1 = k + k

k =k +1
1

. Calcular QN = QG QD QN (V Y )

Comprobar

Q < Comprobar
N

Solucion

P <
N

N
Q
k
V

= [B ]V
k

V k+1 = V k + V
k

k =k +1
2

46

Ejemplo: Resolver el mismo problema del ejemplo anterior, empleando el metodo de

New- ton Raphson desacoplado rapido
Expresiones necesarias
P2 = 0, 1923V 2 0, 1923V2Cos2 + 0, 9615V2Sen2
2
Q = 0, 9415V 2 0, 1923V Sen 0, 9615V Cos
2
2
2
2
2
P = 0, 30 P
2

Q = 0, 07 Q
2

1
12

B =
=1
B = B = 0, 9415

22

Comentarios a la matriz jacobiana del flujo de carga desacoplado rapido

- En la formacion de B no participa la lnea y la columna de la barra (V )
- En la formacion de B no participan las lneas y columnas de las barras (V ) y (PV)
- Dimension de B : NPQ + NPV
- Dimension de B : NPQ

Condiciones iniciales
i) Kp = Kq = 1; p = q = 0; 2o = 0, 0; V2 0 = 1, 0; = 0, 003

Primera iteracion (P )
o )2
o Sen o 2 = 0, 02
2
2
ii) P2 (V o , o ) = 20, 1923(V
0,21923V
Coso 2+ 0, 9615V
0 = o
30
P2 (V o , o ) = 0,
0, 30
2
2
P (V 2o ,

2
o2)

iii) |
| < ? [NO]
0, 30

2 (V2o ,o
2)
P
o
0,30 =
iv)
= B o
11
o
=
V
2
1
o
o
v) 2 =
2 + 2 = 0 0, 30 = 0, 30
1 = 0, 30

vi) p = 0 + 1 = 1
vii) K = 1 ir al paso (x)
q

Primera iteracion (Q-V)

x) Q2 (V20 , 2 1 ) = 0, 9515(V2 o )2 0, 1923V
2 o Sen
2o
2 o Cos
2 o = 0, 07977
0, 9615V
Q2 (V20 , 21 ) = 0, 07 0, 07977 = 0, 00977
xi) |Q2 (V 20 , 21 )| < ?
0, 00977 < 0, 003 [NO]

01038
0 1
o
0,00977 = 0,V
xii) Q2 (Vo2 ,2 ) = B V o
V2
1(0,9415)
2
2

47

xiii) V21 = V2 o + V2 o = 1, 0 0, 01038 = 0, 98962

V 1 = 0, 98962
2

xiv) q = 0 + 1 = 1
xv) kp = 1 (ii)

Segunda iteracion (P )
2 1 + 20, 9615V 1 Sen
2 1
ii) P2 (V 1 , 1 )2 = 20, 1923(V 1 )2 2 0, 1923V 1 Cos
P2 (V 21 , 21 ) = 0, 27467
P2 (V2 1 , 2 1 ) = 0, 30 (0, 2746) = 0, 02533
iii) |P2 (V 21 , 12 ) < ?
0, 02533 < 0, 0033 [NO]
0, 02560

1 1
1
0,02533 =
iv) P2 (V12 ,2 ) = B 1
1(0,98962)
=
2
2
V
v) 22 = 21 + 21 = 0, 30 0, 02560 = 0, 32560
2 = 0, 32560

vi) p = 1 + 1 = 2 vii)
kq = 1 ira(x)

Segunda iteracion (Q-V)

x) Q2 (V21 , 2 2 ) = 0, 9415(V1 2 )2 0, 1923V
2 1 Sen
2 2 0, 9615V
2 1 Cos
2 2
Q2 (V21 , 2 2 ) = 0, 08140
Q2 (V2o , 21 ) = 0, 07 0, 0814 = 0, 0114
xi) |Q2 (V 21 , 22 )| < ? [NO]
1 2
xii) Q2 (V12 ,2 ) = B V 2 V 2 = 0,01140 = 0, 01224
2
2
2
9415(0,98962)
2
V
2
V = 0, 01224 0,
xiii) V22 = V2 1 + V2 1 = 0, 98962 0, 01224 = 0, 97738
V 2 = 0, 97738
2

xiv) q = 1 + 1 = 2
xv) k = 1 ir al paso (ii)
p

Tercera iteracion (P )
ii) P2 (V 22 , 22 ) = 0, 1923(V2 2 )2 0, 1923V2 2 Cos2 2 + 0, 9615V
2 2 Sen
2 2
P2 (V 22 , 22 ) = 0, 29448
P2 (V2 2 , 2 2 ) = 0, 30 (0, 29448) = 0, 00502
iii) |P2 (V 22 , 22 ) < ? [NO]
0, 00514

2 2
2
0,00502 =
iv) P2 (V22 ,2 ) = B 2
1(0,97738)
2
=
2
V

48

v) 23 = 22 + 22 = 0, 32560 0, 00514 = 0, 33073

3 = 0, 33073
2

vi) p = 2 + 1 = 3
vii) k = 1 ir a (x)
q

Tercera iteracion (Q-V)

2 )2
2
3
2 Cos23
2
2
x) Q2 (V 2 , 3 ) = 20, 9415(V
0,2 1923V Sen 2 0, 9615V
Q2 (V22 , 23 ) = 0, 0716
2 2 , 23 ) | = 0, 07
Q
= 0016
|Q2 (V
0,
0, 0016
ya0716
que 0,
< 0, 003
2 (V 2 , 3 ) < ? Si,
2

El subproblema QV convergio y no se requiere de efectuar la iteraci

on. xvi) KQ = 0 e ir al paso (xvii)
xvii) El subproblema P aun no converge ya que Kp = 1 (ii)

Cuarta Iteracion(P ):
ii) P2 (V 22 , 23 ) = 0, 1923(V2 2 )2 0, 1923V2 2 Cos2 3 + 0, 9615V
2 2 Sen
2 3
2
22 23
P
29923
0,
0,
P(V2 (V, 2, )3 =
)=
30 (0, 29923) = 0, 00077
2
2
P (V 22 ,

2
32 )

iii)
| < ? Si, ya que 0, 00077 < 0, 003
viii) |K = 0 (ix)
p

ix) El subproblema QV tambien convergio, ya queKQ =

0 xviii) La solucion
obtenida
2 = 3 = 0, 33073
V2 = V2 = 0, 9774

Observacion:
El error es reducido despues de una iteracion adicional, obteniendoce:
= 0, 33193
2

V2 = 0, 97518

5.10.6

Metodo de Ward y Hale

Se presenta por su valor historico y se desarrolla a nivel informativo.

Concepto basico: Las potencias nodales son muy dependientes de los valores de tension
y angulo de fase de su propio nodo.
Lo anterior unido a las limitaciones en la capacidad de memoria de los computadores de la
epoca (decada de los anos 50), permitio plantear un metodo en el cual solo se considerar
an los elementos diagonales de la matriz Jacobiana.

49


ALGORITMO PARA LA SOLUCIO N DEL PROBLEMA
DE FLUJO DE POTENCIA WARD Y HALE
Parametros de la red (Lneas y
transformadores)
Informacion Nodal

Condiciones iniciales (Voltajes), para

metros
k = 0 (contador
de iteraciones)

. k=k+1
. Calcular los errores de potencia activa y reactiva

.Chequeo de convergencia P Ni < y QNi < Solucion

. Formar la matriz Jacobiana (4 elementos por nodo)

Pi

i
Q
i

Pi
Vi
Q
V i
i

. Encontrar valores incrementales del voltaje

. Actualizar los voltajes

Ecuaciones:El modelo se puede plantear tanto en coordenadas polares como en rectangulares

tal como se describe a continuacion.
Rectangulares:

P
QN i

k+1

= P
QN i /f
/fi P
QN i /e
/ei
Ni

Ni

i
f
e

Ni

Polares:

Ni

PN i

50

PN i /i PN i /Vi

= Q
/ N i
Ni
/Vi
i

k+1

V
i

Los terminos del jacobiano (derivadas parciales) ya fueron determinados al estudiar los m
etodos acoplados de Newton Raphson.
La matriz jacobiana tendra una estructura
diagonal cuyos elementos son (Hkk ), (Lkk ),(Jkk ), (Nkk ).

51

En cada nodo se plantean dos ecuaciones, y cada uno de estos presenta dos inco
gnitas
(Voltaje en magnitud y angulo o voltaje en parte real e
imaginario).
En general para los n nodos se tendran n sistemas y cada uno de estos conformado por
(P Q
i
i
)

5.10.7 Pequena Modificacion en el Metodo Desacoplado R

apido
En algunas aplicaciones de ujo de carga los elementos shunt equivalentes pueden asumir valores de admitancia anormalmente elevados, lo que puede llevar a dicultades en el m
etodo desacoplado rapido. La dicultad aparece debido a que la aproximacion:

|Bkk V 2 | |Qk |

ya no es mas verdadera.
La alteracion ocurre solamente en los elementos de la diagonal de la matriz B.

Analisis de la formacion Original

En la formulacion
original se tiene que:

kk

= Bkk +Qk

B =
kk

mk

bkm b sh +
k

Entonces:

mk

bsh
km

B kk = Lkk =

bkm bsh
k +
k

Analisis Modicado
Q =V
V (G Sen
k

k m

k m km

Q = B V2 +V
k

L kk

bsh
km

km

kk

= Bkk +Qk

km

V (G
k

Cos
km

km

Sen
km k

km

Cos
km

km

Reemplazando las ecuaciones anteriores

Lkk = 2Bkk +
1

Vm (Gkm Senkm Bkm Coskm )

V
k

mk

52

Considerando solamente lneas de transmision (o las simplicaciones correspondientes

en el caso de transformadores desfasadores y puros) se tiene que:
Y = Y = g jb = G + jB
km

km

km

km

km

km

Incorporando los valores de Gkm y Bkm y el valor de Bkk se tendra lo siguiente:

53

Lkk = 2

bkm + bsh +

bsh + 1
Vm
(gkm Sen
km + bkm
km ) V
Cos
m
k

bsh +2
kk =
m
'

bsh bm
+
bmkk

1Cos
+
gkm
km Sen
km

k
Si se consideran las siguientes
( simplicaciones (
( 1
Coskm ( (
b
Vk

k
m ((
y Vm

| (

Que equivaleka emplear las

aproximaciones:

Coskm 1
|Bkm| |Gkm Senkm |
V k Vm 1
Entonces se reducen a:
los
terminos
de L

km
Lkk = bsh b
2

sh
b + k
m

kk =
Lkk
2 =
bsh

bsh bkm

mk

La unica diferencia entre la

version original y la modicada
es el que en esta ultima la
suma de las susceptancias shunt
aparecen multiplicadas por 2.
Es importante observar que la
aproximacion obtenida considera
los = 0 y V = 1,0 p.u.

54