EE354 - Adicional - Flujo de Potencia - 2017-II PDF

Universidad Nacional de Ingeniera
Facultad de Ingeniera Elctrica y Electrnica
Escuela Profesional de Ingeniera Elctrica
Curso: EE354N
Anlisis de Sistemas de Potencia II
Expositor : Manfred F. Bedriana Arons

Doctor en Ingeniera Elctrica
Ingeniero Electricista C.I.P. N 95644
2017
Escuela Profesional de Ingeniera Elctrica, FIEE-UNI
Curso EE354 Anlisis de Sistemas de Potencia II
Material Adicional
Flujo de Potencia
03/09/2017 [email protected] 2
Flujo de Potencia
Conceptos y Mtodo de Gauss-Seidel
A. Introduccin (1)
Para la evaluacin de los sistemas de energa elctrica comnmente se usan
los estudios de flujo de potencia.
Esta herramienta es usada en el planeamiento, operacin econmica y
control de un sistema de potencia.
Dentro de la evaluacin de los sistemas de energa elctrica, las empresas
elctricas requieren:
Calcular las tensiones en barra para valores de carga dados y con la
programacin de la generacin ya definida en el despacho econmico.
Calcular las prdidas activas y reactivas en las lneas de transmisin y en los
transformadores.
Calcular los flujos de potencia activa y reactiva en las lneas de transmisin y
en los transformadores.
A. Introduccin (2)
Estudios de flujo de potencia herramientas computacionales.
A. Introduccin (3)
Los algoritmos de flujo de potencia usan la matriz admitancia de barra YBUS
debido al alto grado de dispersidad de las matrices de red de sistemas reales.
Grado de dispersidad (GD) Porcentaje de elementos nulos de la matriz.
nmero de elementos nulos

GD 100%
nmero total de elementos
NB 2 ( NB 2 NR)
GD 100%
2
NB
Ejemplo: Sistema brasileo
(sur/sudeste/centro-oeste)
NB = 1663 barras
NR = 2349 ramas
16632 (1663 2 2349)

GD 100% 99.77%
16632
B. Evolucin del flujo de potencia (1)

Antes se usaban los analizadores de redes (network analyzers). Consistan de paneles
con componentes del sistema modelados a escala por medio de resistores,
inductancias y capacitores.
El primer mtodo prctico de solucin del problema de flujo de potencia a travs de
un computador digital Ward, Hale 1956 (mtodo de la matriz Ybus)
Ward, J. B., Hale, H. W., Digital Computer Solution of Power-Flow Problems," AIEE
Transactions, Part III Power Apparatus and Systems, Vol. 75, pp. 398-404, June 1956.
Mtodo de Newton aplicado al flujo de potencia (caractersticas de convergencia
excelentes) Van Ness 1959
J.E. Van Ness, "Iteration methods for digital load flow studies", Trans. AIEE (Power
Apparatus and Systems), Vol. 78A, p. 583, August1959.
Tcnicas de almacenamiento compacto (rapidez en los mtodos y exige menos
cantidad de memoria) Tinney, Walker 1967.
W. F. Tinney, and J .W. Walker, Direct solutions of sparse network equations by
optimally ordered triangular factorization, Proc. IEEE 55, 1801-1809, 1967.
Las empresas de electricidad consideran al Newton como mejor mtodo.
B. Evolucin del flujo de potencia (2)

Mtodo desacoplados rpidos (mtodos ms rpidos que el Newton convencional)
Stott, Alsac 1974 (mtodo X-B)
B. Stott, O. Alsac, Fast decoupled load flow, IEEE Transactions on Power Apparatus
and Systems, 1974.
Van Amerongen1989 (mtodo B-X)
R.A.M. van Amerongen, A general-purpose version of the fast decoupled load flow, IEEE
Transactions on Power Systems, 1989.
Monticelli, 1990 (hiptesis y fundamentos del mtodo B-X)
A.J. Monticelli, A. Garcia, O.R. Saavedra, Fast decoupled load flow: hypothesis,
derivations, and testing, IEEE Transactions on Power Systems, 1990.
Mtodo robustos (mtodos Newton modificado con deteccin de factibilidad)
Iwamoto 1981 (optimizacin de paso en coordenadas rectangulares)
S. Iwamoto, S. Tamura, A load flow calculation method for ill-conditioned power
systems, IEEE Transactions on Power Apparatus and Systems, vol.PAS-100, n.4, 1981.
Castro 2000 (optimizacin de paso en coordenadas polares)
Luciana M. C. Braz, Carlos A. Castro, Carlos A. E Murari, A critical evaluation of step
size optimization based load flow methods, IEEE Transactions on Power Systems, vol.15,
n.1, 2000.
C. Conceptos iniciales (1)

El estudio de flujo de potencia asume lo siguiente:
Operacin a condiciones balanceadas (cargas balanceadas y lneas transpuestas).
Los componentes son modelados con un equivalente monofsico.
Cada barra est asociada a 4 cantidades: inyeccin de potencia activa (P),
inyeccin de potencia reactiva (Q), mdulo de la tensin en barra (V), ngulo de
la tensin en barra ().
Notacin:
Potencia de inyeccin que
entra a la barra (+)
Potencia de inyeccin que sale
de la barra ()

Tipos de barra en la red
Barra de carga (PQ).
Se conoce la inyeccin de potencia activa (P) y reactiva (Q).
Barra de generacin (PV)
Se conoce la inyeccin de potencia activa (P) y el mdulo de la tensin (V).
Barra slack (barra swing o barra V)
Se conoce la magnitud de la tensin (V) y el ngulo ().
La barra slack puede representar a un gran sistema de potencia (barra infinita).
Representa la referencia angular del sistema (ngulo conocido).
La potencia de la barra slack asume las prdidas activas en la red. Comnmente
debe ser un generador con mayor capacidad de regulacin de frecuencia.

Tipos de barra en la red
Ejemplo sistema de 3 barras
Generacin Generacin Generacin barra 2

Generacin barra 1
Inyeccin de potencia positiva
Inyeccin de potencia positiva
Carga barra 1 Carga barra 2

Inyeccin de potencia negativa Inyeccin de potencia negativa
Nota:
Pueden conectarse cargas en Carga barra 3
las barras slack y PV. Inyeccin de potencia negativa
Estas barras no pierden su

condicin de tipo de barra a
pesar de la carga.

Que significa tener una referencia de mdulo y ngulo de tensin?
(consulta hecha en clase)
Todo sistema elctrico AC debe poseer una referencia para los mdulos y
ngulos de las tensiones.
En el caso de los sistemas de potencia se tiene:
1. La referencia de mdulos de tensiones de barra est implcitamente expresado al
definir en el nodo de tierra (punto equipotencial, valor de 0 Volts de tensin).
Nota: Al formar la matriz Ybus incluimos implcitamente el nodo de tierra, los trminos
de shunt yi0 estn definidos para Vi 0 y V0 = 0.
2. La referencia de ngulos la tensin lo da la barra slack, es decir, sea cualquier valor
del ngulo de tensin de la barra slack, este servir de referencia para los dems
ngulos de tensin en barra calculados en la red elctrica.
Barra ngulo de barra slack Barra i

slack (referencia de ngulos
de V del sistema)
Vslack slack Shunt: yi0
Nodo de tierra 0
0
(referencia de mdulos de V del sistema)

Que significa tener una referencia de mdulo y ngulo de tensin?
Ejemplo.
Cambio de ngulo de la barra slack
Barra Barra Barra Barra

Barra Barra
slack PQ slack PQ
Transf. PV PV
Transf.
1 2 3 1 2 3
Lnea Lnea
G G
. . . .
P12 = 105 MW P23 = 100 MW Carga P12 = 105 MW P23 = 100 MW Carga
Q12 = 50 MVAR Q23 = 20 MVAR Q12 = 50 MVAR Q23 = 20 MVAR
Referencia angular del sistema I: slack = 0 Referencia angular del sistema II: slack = 10
(equivalente a sumar 10 a todos los ngulos de las
barras del sistema I)
Los valores de flujo de potencia se mantienen
iguales (solo dependen de la diferencia angular)

Flujos de potencia por las redes elctricas
1
2
Nota:
En toda barra se cumple que la
suma algebraica de inyecciones
de potencia (total que entra y
3 total que sale) es igual a cero.
Ejemplo en barra 1:
+ Pg1 Pc1 P12 P13= 0

Prdidas elctricas por las redes
Ejemplo: lnea de transmisin
Prdidas de transmisin

Balance de potencia (funcin de la barra slack)
Ejemplo sistema de 3 barras (barra 2 con carga)
Pprdidas1
Pprdidasi
Pprdidas2 Pprdidas3

Diagrama de admitancias + fuentes de inyeccin de corriente
(modelo lineal)
Ejemplo sistema de 3 barras (barra 2 con carga)
Carga o generacin modelada

1 2 como inyecciones de corriente
3
Nota:
Las variables del sistema se Susceptancias de las lneas +
resuelven a travs de: compensadores shunt
I = YV
(relacin lineal entre I y V)

Diagrama de admitancias + fuentes de inyeccin de potencia
(modelo no lineal) = Problema de flujo de potencia
Ejemplo sistema de 3 barras (generacin y carga)
Elementos contenidos en la
matriz Ybus
Carga o generacin modelada

como inyecciones de potencia
Carga (-)
Generador (+)
D. Ecuacin de flujo de potencia (1)

KCL (Ley de corrientes de Kirchhoff )
Barras vecinas a barra i
Admitancias en p.u.
Inyeccin de corriente en la barra i:
D. Ecuacin de flujo de potencia (2)

Ecuaciones de inyeccin de potencia
Relacin no lineal entre P y V
La formulacin matemtica anterior del problema de flujo de potencia

resulta en un conjunto de ecuaciones algebraicas no lineales que deben ser
resueltas por mtodos iterativos:
Mtodo Gauss-Seidel
Mtodo de Newton
E. Mtodo de Gauss-Seidel
Resuelve ecuaciones algebraicas no lineales.
Mtodo de desplazamientos sucesivos.
Mtodo iterativo:
Toma una funcin y la ordena de la siguiente forma x = g (x)
{pueden haber diversos arreglos de este tipo}
Hacer una estimacin inicial de la variable x : x[0] valor inicial
Encontrando una mejora iterativa de x[k], que es x[k+1] = g (x[k])
Se alcanza una solucin cuando la diferencia entre dos iteraciones
es menor que una exactitud especificado : | x[k+1] - x[k] |
Factores de aceleracin
Se puede mejorar la velocidad de convergencia : > 1
Modificacin del paso: la mejora se encuentra como
x[k+1] = x[k] + ( g( x[k] ) x[k] )
F. Ejemplo de mtodo de Gauss-Seidel (1)
Encontrar las races de la ecuacin: f (x) = x3 6x2 + 9x 4 = 0
Paso 1: Despejar la ecuacin g (x) de la siguiente forma :
9 x x3 6 x 2 4
1 6 4
x x 3 x 2 g ( x)
9 9 9
Paso 2: A partir de un valor inicial x[0] = 2, se llevan acabo varias iteraciones.
1 6 4
x[1] g ( x[ 0] 2) (2) 3 (2) 2 2.2222
9 9 9
1 6 4
x[ 2] g ( x[1] 2.2222) (2.2222) 3 (2.2222) 2 2.5173
9 9 9
1 6 4
x[3] g ( x[ 2] 2.5173) (2.5173)3 (2.5173) 2 2.8966
9 9 9
x[ 4] 3.3376
x[5] 3.7398
x[ 6] 3.9568
x[ 7 ] 3.9988
x[8] 4.0000 Convergencia a una solucin

con un factor de aceleracin de 1.25
Partiendo con un valor inicial de x[0] = 2
1 6 4
g ( x[ 0 ] 2) (2)3 (2) 2 2.2222
9 9 9
x[1] 2 1.25[2.2222 2] 2.2778 Aumento de velocidad
1 6 4
g ( x[1] 2.2778) (2.2778)3 (2.2778) 2 2.5902
9 9 9
x[ 2] 2.2778 1.25[2.5902 2.2778] 2.6683

Comparacin
Sin factor de aceleracin Con factor de aceleracin
Disminucin del nmero de

iteraciones
G. Gauss-Seidel para sistemas de ecuaciones (1)

Considerando un sistema de n ecuaciones
Ordenando cada ecuacin para una de las variables
G. Gauss-Seidel para sistemas de ecuaciones (2)

Pasos
Se asume una solucin aproximada para las variables independientes.
Se busca el resultado en una nueva solucin aproximada
En el mtodo de Gauss-Seidel, los valores actualizados de las variables

calculadas en las ecuaciones anteriores se utilizan inmediatamente en
la solucin de las ecuaciones siguientes.
H. Ecuaciones de flujo de potencia (1)

La ecuacin inicial de inyeccin de potencia
En la forma de Gauss-Seidel

Reescribiendo las ecuaciones de inyeccin potencia para hallar P y Q
Inyeccin de potencia activa y reactiva

Barras de carga (PQ): P y Q permanecen fijas.
Barras de generacin (PV): P permanecen fija.
Las corrientes y las potencias se expresan como inyectndose a las barras:

Para generacin la inyeccin es positiva.
Para cargas la inyeccin es negativa.

El conjunto de ecuaciones son:
La potencia especificada Psch y Qsch es igual a la diferencia entre generacin y

carga de potencia activa y reactiva, respectivamente. Son trminos constantes
en todo el proceso de convergencia.

Diferencia entre potencia especificada e inyeccin de potencia
Potencia especificada en la barra k

(generacin menos carga )
Pksch = Pgk Pck
Potencia total transmitida por las

lneas desde la barra k
(potencia de inyeccin en barra k)

Ordenando las ecuaciones en trminos de los elementos de Ybus
I. Solucin del mtodo de Gauss-Seidel

Caractersticas del sistema
Ya que ambas componentes (V y ) son especificadas para la barra slack
hay 2(n 1) ecuaciones que deben resolverse en forma iterativa, esto porque
el nmero total de variables V y es igual a 2(n 1).
Para las barras de carga (PQ), la potencia activa y reactiva son conocidas:
El mdulo de la tensin y ngulo deben ser estimados.
Los ngulos en general estn muy prximos, por lo que un valor inicial de
0 es el apropiado.
Para las barras de generacin (PV), la potencia activa y tensin son

conocidas
La potencia activa es programada.
La potencia reactiva es calculado sobre el valor estimado de tensin.
Las tensiones complejas se hallan como:
donde V sch es el valor de tensin especificado para la barra PV.

J. Ejemplo sistema 3 barras (1)

Sea el diagrama unifilar de impedancias por unidad (Sbase = 100 MVA).
Usando el mtodo de Gauss-Seidel, determine el valor de la tensin fasorial

en las barras 2 y 3, con una precisin de 2 decimales.

Se obtiene el diagrama unifilar de admitancia por unidad.
En las barras de carga son calculados los valores especificados de potencia

por unidad:

Valores iniciales:
Se aplica el mtodo de Gauss-Seidel, llevando a cabo 7 iteraciones:

Proceso de Gauss-Seidel converge al cabo de 7 iteraciones:
Despus de calculadas las tensiones en las barras (mdulo y ngulo), se

calcula la potencia de inyeccin de la barra slack
K. Flujo de potencia y prdidas en la lnea (1)

Despus de calcular los mdulos y ngulos de las tensiones en barras se
calculan los flujos de potencia y las prdidas en las ramas de la red.
Las lneas y los transformadores son ramas de la red elctrica.
La direccin del flujo de corriente se define como mostrado en la Figura para
un elemento del sistema de potencia (ejemplo: lnea de longitud mediana)
Sij Sji
El flujo de potencia se define para cada extremo de la rama, similar al flujo de

corrientes.
Ejemplo: El flujo de potencia Sij sale de la barra i y fluye hacia la barra j.
K. Flujo de potencia y prdidas en la lnea (2)

El flujo de potencia y la corriente son definidos como:
Prdidas de potencia: Sprdidas ij = Sij + Sji
Nota:
Sij Sji Como se ve en la figura, las

prdidas elctricas entre las
barras i,j involucran la suma
de prdidas elctricas de los
elementos: yi0, yij e yj0.
L. Ejemplo sistema 3 barras

Conocidas las tensiones se calculan los flujos de potencia en las lneas.
Barra i
Pij en MW
(Pprdidas, ij en MW)
Qij en MVAR (Qprdidas, ij en MVAR)
Las prdidas son calculadas como:
L. Ejemplo sistema 3 barras (MATLAB) (1)

Datos: Montar archivo de entrada
DataBus3.m
Diagrama unifilar de impedancias

por unidad (Sbase = 100 MVA)
Archivo de entrada: DataBus3.m

disp(' This program saves bus data and line data in the MATLAB workspace')
disp(' ')
basemva = 100;
disp(' 1. busdata --> workspace')

% Bus Bus Voltage Angle ---Load---- -------Generator----- Static Mvar
% No code Mag. Degree MW Mvar MW Mvar Qmin Qmax +Qc/-Ql
busdata=[1 1 1.05 0.0 0.0 0.0 0.0 0.0 -999 999 0
2 0 1.0 0.0 256.6 110.2 0.0 0.0 0 0 0
3 0 1.0 0.0 138.6 45.2 0.0 0.0 0 0 0];
disp(' 2. linedata --> workspace')

% Line code
% Bus bus R X 1/2 B = 1 for lines
% nl nr p.u. p.u. p.u. > 1 or < 1 tr. tap at bus nl
linedata=[1 2 0.02 0.04 0.0 1
1 3 0.01 0.03 0.0 1
2 3 0.0125 0.025 0.0 1];
L. Ejemplo sistema 3 barras (MATLAB) (2)

Simular flujo de carga
Resolver mediante el mtodo de Gauss-Seidel en MATLAB. Usar los sgtes programas:
lfybus: forma la matriz admitancia de barra
lfgauss: ejecuta el mtodo de Gauss-Seidel para flujo de potencia
Mostrar resultados
Usar los sgtes programas:
busout: reporte de resultados de barras
lineflow: reporte de resultados de flujos de potencia
Barra i
Pij en MW
(Pprdidas, ij en MW)
Qij en MVAR
(Qprdidas, ij en MVAR)
M. Ejemplo MATLAB (1)

En la Figura se ve el famoso sistema
IEEE 30 barras, el cual proviene del
American Electric Power Service Co.
Datos
Resolver mediante el mtodo de Gauss-

Seidel en MATLAB:
usar programas lfybus, lfgauss
Mostrar resultados: busout, lineflow
M. Ejemplo MATLAB (2)

Porqu en el bloque de barras del formato de datos Data_IEEE30.m se debe
llenar valores de tensin en las barras de carga (por ejemplo en 1 p.u.) y no
valores en 0, a pesar que se supone que el problema de flujo de potencia no fija
la tensin de estas barras?
Formato de barras del archivo Data_IEEE30.m
disp(' 1. busdata --> workspace')
% IEEE 30-BUS TEST SYSTEM (American Electric Power)
% Bus Bus Voltage Angle ---Load---- -------Generator----- Static Mvar
% No code Mag. Degree MW Mvar MW Mvar Qmin Qmax +Qc/-Ql
busdata=[1 1 1.06 0.0 0.0 0.0 0.0 0.0 0 0 0
2 2 1.043 0.0 21.70 12.7 40.0 0.0 -40 50 0
3 0 1.0 0.0 2.4 1.2 0.0 0.0 0 0 0
4 0 1.06 0.0 7.6 1.6 0.0 0.0 0 0 0
5 2 1.01 0.0 94.2 19.0 0.0 0.0 -40 40 0
6 0 1.0 0.0 0.0 0.0 0.0 0.0 0 0 0
7 0 1.0 0.0 22.8 10.9 0.0 0.0 0 0 0
8 2 1.01 0.0 30.0 30.0 0.0 0.0 -30 40 0
9 0 1.0 0.0 0.0 0.0 0.0 0.0 0 0 0
Los valores que se llenan de mdulo y ngulo de las tensiones de las barras PQ son
usados para establecer el punto inicial. El mtodo de Gauss-Seidel, as como el
mtodo de Newton, necesitan un buen punto inicial para que el proceso se acerque a
zona de convergencia de la solucin deseada. En caso contrario el proceso diverge o
entra a otra zona de convergencia.
Flujo de Potencia
Mtodo de Newton-Raphson
y Desacoplados
N. Mtodo de Newton-Raphson (1)

Convergencia cuadrtica
Mtodo matemtico superior al de Gauss-Seidel.
Ms eficiente para redes grandes

El nmero de iteraciones requerido para la solucin es independiente del
tamao del sistema.
Las ecuaciones de Newton-Raphson son repartidas en forma natural al

sistema de potencia.
Resolver el mdulo y ngulo de las tensiones, teniendo especificadas las
inyecciones de potencia activa y reactiva.
La calidad de la solucin del mtodo de Newton-Rapshon depende de la

proximidad del punto inicial a la solucin.

Es un mtodo de aproximacin sucesiva que usa la serie de Taylor
Considere la funcin : f (x) = c, donde x no se conoce.
Sea una estimacin inicial x[0], tambin tenemos que x[0] es una pequea
desviacin que va hacia la solucin correcta
f (x[0] + x[0] ) = c Mtodo optimista, porque supone que

con x[0] se encuentra la solucin
Realizando la expansin de Taylor alrededor del punto x[0]

Asumiendo que el error, x[0] , es pequeo, se desprecia los trminos de mayor
orden, resultando as
donde
Reordenado la ecuacin se tiene:
Pendiente de la recta tangente a la

curva f(x) en el punto x[0]
O. Ejemplo mtodo de Newton
Recta tangente a la curva

f(x) en el punto x[0] = 6
A partir de un valor inicial x[0] = 6, se llevan acabo 3 iteraciones.

Menor nmero de iteraciones comparado a Gauss-Seidel (7 iteraciones).
P. Ecuaciones de Flujo de Potencia

Ley de Kirchhoff para inyeccin de corrientes
ij ngulo de la
admitancia Yij
Inyeccin de potencia activa y reactiva
Sustituyendo cada Ii
Dividiendo la potencia en su parte real y reactiva
Q. Matriz Jacobiana (1)

Poniendo las ecuaciones en su forma iterativa
Formando las ecuaciones matriciales del sistemas de ecuaciones
Forma general de la ecuacin para encontrar una solucin
La ecuacin iterativa
Matriz Jacobiana
Q. Matriz Jacobiana (2)

La matriz Jacobiana contiene las primeras derivadas del conjunto de
funciones matriciales P y Q respecto a las variables V y .
R. Proceso iterativo del mtodo de Newton

Mismatches o residuos de potencia
Definido como la diferencia entre la potencia especificada Psch y la calculada Pi
(llamada tambin potencia de inyeccin).
La potencia especificada Psch y Qsch es igual a la diferencia entre generacin y

carga de potencia activa y reactiva, respectivamente. Son trminos constantes en
todo el proceso de convergencia.
Estimamos los desvos de los mdulos y ngulos de tensiones de barra.
Actualizamos las estimaciones de mdulos y ngulos de tensiones de barra.
S. Tipos de Barra y formacin de la matriz Jacobiana (1)

Barra slack o barra swing
Una barra con generador debe ser seleccionada y definida como la referencia
angular
El mdulo y ngulo de la tensin son conocidos para esta barra.
El ngulo se selecciona arbitrariamente a 0.
Esta barra no es incluida en la formacin de la matriz Jacobiana.
Barra generador (PV)
El mdulo de la tensin y la potencia activa de inyeccin son conocidos.
El ngulo de la tensin y la potencia reactiva se deben calcular.
Esta barra esta incluida en la parte de potencia activa de la matriz Jacobiana.
Barra de carga (PQ)
Se conocen la inyeccin de potencia activa y reactiva.
Esta barra esta totalmente incluida en la matriz Jacobiana (parte de potencia activa
y reactiva).
S. Tipos de Barra y formacin de la matriz Jacobiana (2)

Ejemplo de matriz Jacobiana
Barra Datos Incgnitas
1
Generacin Generacin 2
3
Matriz Jacobiana
P2 P2 P2

P2 2 3 V3
2
P P3 P3 P3

3 3 V3
3
Q3 2 V3
Q 3 Q3 Q3
2 3 V3
P2 P2 P2

2 3 V3
2
P P3 P3
3
V3
3
3
2 V3
Q3 Q3 Q3
2 3 V3
T. Pasos del mtodo de Newton-Raphson (1)

1. Establecer el punto flat start (llamado de mnima energa)
Para barras de carga, se establecen valores iguales a las tensiones de la barra
slack o 100
Para barras de generador, se igualan los ngulos con la barra slack o 0
2. Calculo de los residuos de potencia (mismatches de potencia)

Para cada barra de carga, calcular las inyecciones P y Q usando las tensiones
estimadas del sistema.
Para barras de generacin, calcular la inyeccin P.
Obtener las potencias especificadas (constantes en todo el proceso)
Obtener la variacin de potencia, P y Q
3. Formando la matriz jacobiana

Usando las derivadas parciales de P y Q respecto a los mdulos y ngulos de
tensiones en barra.
T. Pasos del mtodo de Newton-Raphson (2)

4. Encontrar la solucin matricial (eligiendo a o b)
a. Usar la inversa de la matriz Jacobiana y multiplicar por la variacin de
potencia.
b. Realizar la eliminacin gaussiana sobre la matriz Jacobiana con el vector b

igual a los residuos de potencia
Calculamos y V.
5. Encontrar nuevas estimaciones para los mdulos y ngulos de tensin.
6. Repetir el proceso hasta que los residuos de potencia (mismatches) sean
menor que una tolerancia especificada.
U. Ejemplo MATLAB
Resolver mediante el mtodo de

Newton en MATLAB:
usar programas lfybus, lfnewton
Que diferencia encuentra entre el

mtodo Gauss-Seidel y el mtodo de
Newton?
V. Flujo de potencia desacoplado rpido (1)

Para lneas y transformadores que tienen una relacin X /R alto:
Los cambios de la potencia activa P
Son menos sensibles a cambios en el mdulo de la tensin V
Son ms sensibles a cambios en el ngulo de fase
Los cambios en la potencia reactiva Q

Son menos sensibles a cambios in el ngulo de fase,
Son ms sensibles a cambios en la mdulo de la tensin, V
Las submatrices jacobianas JQ y JPV tienden a ser mas pequeas en

comparacin a los valores de JP y JQV
Las submatrices de la matriz Jacobiana JQ y JPV se pueden aproximar

iguales a cero .
B. Stott, O. Alsac, Fast decoupled load flow, IEEE Transactions on Power
Apparatus and Systems, 1974.

Elementos JPV
Elementos JQ

La ecuacin matricial es separada en dos ecuaciones desacopladas
Requiere mucho menos tiempo para resolver el sistema desacoplado en

comparacin al mtodo extendido de Newton-Raphson.
Las submatrices JP y JQV pueden simplificarse aun ms para eliminar la

necesidad de recalcular las submatrices durante cada iteracin.
Algunos trminos en cada elemento son relativamente pequeos y se pueden
despreciar.
Esto resulta en una matriz Jacobiana con trminos constantes.
W. Trmino de la diagonal JP del Jacobiano
X. Trmino fuera de la diagonal JP del Jacobiano
Y. Trmino de la diagonal JQV del Jacobiano
Z. Trmino fuera de la diagonal JQV del Jacobiano
AB. Mtodo de Flujo de potencia desacoplado rpido

Ecuacin de variacin de potencia de cada termino de las matrices JP y JQV
Formado por los valores de

la parte imaginaria de la
matriz Ybus (Y = G + jB),
pero estructuras diferentes!
B diferente a B
Ecuacin matricial para las matrices JP y JQV
La matriz Ybus tiene una

dimensin 2n x 2n (n nmero
de barras).
La matriz B tiene una
dimensin menor que n x n
AC. Ejemplo mtodo desacoplado rpido (1)

Encontrar la solucin de flujo de potencia usando el flujo de potencia
desacoplado rpido.
Las admitancias del sistema son:

La matriz B ser igual a:
Primera iteracin:

Primera iteracin
Proceso iterativo
AD. Ejemplo MATLAB

Resolver mediante el mtodo de

Desacoplado Rpido en MATLAB:
usar programas lfybus, lfdecouple
Que diferencia encuentra entre el

mtodo Newton y el mtodo
desacoplado rpido?
AE. Comparacin entre mtodo Newton y desacoplado

Mtodo de Newton Mtodo Desacoplado Rpido
Recta tangente a la curva

g(x) en el punto x(0)
Prximo de
la solucin
Siguientes rectas siguen siendo tangentes a la curva g(x) Todas las recta tiene la misma pendiente
diferentes matrices usadas en todo el proceso de misma matriz usada en todo el proceso de
convergencia (matriz Jacobiana se actualiza en c/iteracin) convergencia (matriz B o B)
AF. Transformadores y control de tensin (1)

Transformador considerando una impedancia de dispersin ykm
Skm Smk
Modelo del transformador

Efecto al cambiar el tap
Posicin nominal del tap
Se obtiene el siguiente circuito

Posiciones del tap
a = 1 p.u. (valor nominal)
a < 1 p.u.
Nota:
Efecto capacitivo en barra k
Efecto inductivo en barra m
a > 1 p.u.
Nota:
Efecto inductivo en barra k
Efecto capacitivo en barra m

EE354 - Adicional - Flujo de Potencia - 2017-II PDF

Cargado por

Copyright:

Formatos disponibles

EE354 - Adicional - Flujo de Potencia - 2017-II PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

EE354 - Adicional - Flujo de Potencia - 2017-II PDF

Cargado por

Copyright:

Formatos disponibles

Universidad Nacional de Ingeniera

Facultad de Ingeniera Elctrica y Electrnica

Escuela Profesional de Ingeniera Elctrica

Expositor : Manfred F. Bedriana Arons

nmero de elementos nulos

16632 (1663 2 2349)

B. Evolucin del flujo de potencia (1)

B. Evolucin del flujo de potencia (2)

C. Conceptos iniciales (1)

C. Conceptos iniciales (2)

C. Conceptos iniciales (3)

Generacin Generacin Generacin barra 2

Carga barra 1 Carga barra 2

Estas barras no pierden su

C. Conceptos iniciales (4)

Barra ngulo de barra slack Barra i

C. Conceptos iniciales (5)

Barra Barra Barra Barra

Q12 = 50 MVAR Q23 = 20 MVAR Q12 = 50 MVAR Q23 = 20 MVAR

C. Conceptos iniciales (6)

C. Conceptos iniciales (7)

C. Conceptos iniciales (8)

C. Conceptos iniciales (9)

Carga o generacin modelada

C. Conceptos iniciales (10)

Carga o generacin modelada

D. Ecuacin de flujo de potencia (1)

Barras vecinas a barra i

Inyeccin de corriente en la barra i:

D. Ecuacin de flujo de potencia (2)

Relacin no lineal entre P y V

La formulacin matemtica anterior del problema de flujo de potencia

x[k+1] = x[k] + ( g( x[k] ) x[k] )

F. Ejemplo de mtodo de Gauss-Seidel (1)

Encontrar las races de la ecuacin: f (x) = x3 6x2 + 9x 4 = 0

Paso 1: Despejar la ecuacin g (x) de la siguiente forma :

F. Ejemplo de mtodo de Gauss-Seidel (2)

Paso 2: A partir de un valor inicial x[0] = 2, se llevan acabo varias iteraciones.

F. Ejemplo de mtodo de Gauss-Seidel (3)

F. Ejemplo de mtodo de Gauss-Seidel (4)

Encontrar las races de la ecuacin: f (x) = x3 6x2 + 9x 4 = 0

Partiendo con un valor inicial de x[0] = 2

F. Ejemplo de mtodo de Gauss-Seidel (5)

Disminucin del nmero de

G. Gauss-Seidel para sistemas de ecuaciones (1)

Ordenando cada ecuacin para una de las variables

G. Gauss-Seidel para sistemas de ecuaciones (2)

Se busca el resultado en una nueva solucin aproximada

En el mtodo de Gauss-Seidel, los valores actualizados de las variables

H. Ecuaciones de flujo de potencia (1)

H. Ecuaciones de flujo de potencia (2)

Inyeccin de potencia activa y reactiva

Las corrientes y las potencias se expresan como inyectndose a las barras:

H. Ecuaciones de flujo de potencia (3)

La potencia especificada Psch y Qsch es igual a la diferencia entre generacin y

H. Ecuaciones de flujo de potencia (4)

Potencia especificada en la barra k

Pksch = Pgk Pck

Potencia total transmitida por las