DESPACHO ECONOMICO DESPACHO ECONOMICO
CAGS
OPERACIN OPTIMA DE SISTEMAS
DE POTENCIA
Conlleva:
Economa en la operacin.
Seguridad del sistema
Emisiones de plantas que operan con base en
CAGS
Emisiones de plantas que operan con base en
combustibles fsiles.
Vertimientos de agua en embalses.
Por ahora solo se analizar el primer aspecto,
conocido con el nombre de: Problema del Despacho
Econmico
El objetivo mas frecuente del despacho
econmico, es minimizar el costo total de
la generacin de potencia activa (MW) de
varias plantas generadoras, a la vez que
DESPACHO PTIMO
CAGS
varias plantas generadoras, a la vez que
se logra alimentar todas las cargas
adecuadamente.
Aun cuando las plantas hidrulicas,
pueden tener un costo de operacin
muy bajo, estn limitadas por la
disponibilidad de agua, cuando se
OPERACIN OPTIMA DE SISTEMAS
DE POTENCIA
CAGS
disponibilidad de agua, cuando se
contempla un periodo largo.
Por simplicidad se iniciar analizando
plantas trmicas.
Flujo de carga. Hay dos variables conocidas y
dos desconocidas por cada barra. Las
variables conocidas son: potencia activa y
reactiva para las barras PQ, la potencia activa
y la magnitud del voltaje en las barras PV y
OPERACIN OPTIMA DE
SISTEMAS DE POTENCIA
CAGS
y la magnitud del voltaje en las barras PV y
magnitud y ngulo de voltaje en la barra Slack.
Si se permite que las variables especificadas
previamente (potencia activa, potencia reactiva
y magnitud de voltaje) varen en un rango, se
tendran infinitas soluciones al flujo de carga.
Operacin econmica de generadores en un
barraje.
Costo de operacin de un generador.
El principal componente del costo de operacin
de un generador es el costo del combustible.
DESPACHO ECONOMICO
CAGS
de un generador es el costo del combustible.
El costo de mantenimiento es menor. El costo
de combustible es especialmente alto en
plantas trmicas y nucleares. En plantas
hidrulicas prcticamente no tiene costo, por
la facilidad de almacenar el agua.
DESPACHO ECONOMICO
La curva de entrada-salida de una unidad
trmica, se puede expresar en millones de
kilocaloras/hora o costo de operacin ($)/hora
versus salida en MW. Este costo se puede
determinar experimentalmente.
CAGS
CURVA Costo - MW
0
200
400
600
800
1000
1200
1400
1600
8
8
.
3
8
.
6
8
.
9
9
.
2
9
.
5
9
.
8
1
0
.
1
1
0
.
4
1
0
.
7
1
1
1
1
.
3
1
1
.
6
1
1
.
9
1
2
.
2
1
2
.
5
1
2
.
8
1
3
.
1
MW
C
o
s
t
o
DESPACHO ECONOMICO
La curva anterior, tiene un punto de
MW
min
, abajo del cual, puede no ser
econmicamente factible despachar la
mquina o puede ser tcnicamente
CAGS
mquina o puede ser tcnicamente
imposible su operacin.
Despacho Econmico
CURVA Costo - MW
0
200
400
600
800
1000
1200
1400
1600
C
o
s
t
o
CAGS
0
8
8
.
3
8
.
6
8
.
9
9
.
2
9
.
5
9
.
8
1
0
.
1
1
0
.
4
1
0
.
7
1
1
1
1
.
3
1
1
.
6
1
1
.
9
1
2
.
2
1
2
.
5
1
2
.
8
1
3
.
1
MW
MW moneda b P a IC
i Gi i i
/ ) ( + =
Gi
i
dP
dC
A la pendiente de la curva de costos, se le conoce
como: Costo incremental del combustible (IC)
hora moneda d P b P a C
i Gi i Gi i i
/
2
1
2
+ + =
Despacho ptimo en un barraje
Si a priori se conoce cuales generadores (P
Gi
) van a
ser despachados, para satisfacer una demanda P
D
en un barraje determinado
P P
Es una desigualdad estricta con el
fin de disponer de reserva rodante
CAGS
D Gi
P P
max ,
Donde P
Gi
corresponde a la potencia nominal del
generador i. Se cumple adems que:
k i P P P
Gi Gi Gi
,......, 2 , 1 ,
max , min ,
=
fin de disponer de reserva rodante
Despacho ptimo en un barraje
El costo de operacin es insensible a la carga reactiva de los
generadores. Es despreciable
La pregunta es: Cual es la forma optima por la cual se logra
atender la demanda P
D
de la carga, entre los generadores
conectados a la barra?
La respuesta a esto se obtiene, si minimizamos el costo operativo:
CAGS
=
=
k
i
Gi i
P C C
1
) (
Cumpliendo a su vez, con la desigualdad indicada previamente,
para cada generador. En tal caso:
0
1
=
=
k
i
D Gi
P P
Ecuacin de balance de
Potencia
Despacho ptimo en un barraje
Dado que es no lineal y adems C
j
es
independiente de , entonces, este es un
problema de programacin no lineal separable
) (
Gi i
P C
) ( i j P
Gi
Si nos olvidamos provisionalmente de la restriccin
relacionada con la desigualdad previa, el problema
CAGS
relacionada con la desigualdad previa, el problema
se puede resolver por una metodologa de
optimizacin conocida como Multiplicadores de
Lagrange. El Langragiano se define como:
= =
=
k
i
k
i
D Gi Gi i
P P P C
1 1
) (
L
Multiplicador de Lagrange
Despacho ptimo en un barraje
La optimizacin se logra cuando
C
P
Gi
0
= =
L
CAGS
k i
P
C
Gi
i
,...., 2 , 1 ; = =
Gi
i
P
C
Donde
es el costo incremental del i
th
generador
en $/MWh
=
= =
Gk
k
G G
P
C
P
C
P
C
.. ..........
2
2
1
1
Despacho ptimo en un barraje
De acuerdo a lo anterior, el punto de
cargabilidad optima de los generadores,
corresponde al punto en el cual, el costo
incremental de todos ellos es igual. A las
ecuaciones anteriores se les denomina,
CAGS
ecuaciones anteriores se les denomina,
Ecuaciones de Coordinacin. El conjunto de
ecuaciones a solucionar esta compuesto de n
derivadas para n generadores y la ecuacin de
demanda. En total son n+1 ecuaciones. Las
incgnitas son las n generaciones y el
multiplicador de Lagrange.
Despacho ptimo en un barraje
Ejemplo
Dos generadores. Ecuaciones de costos:
4 * 30 * 125 . 0
5 * 40 * 1 . 0
2
2
2 2
1
2
1 1
+ + =
+ + =
G G
G G
P P C
P P C
2 , 1 , 125 20 = i P
Gi
CAGS
Se requiere generar 150 MW. Posibilidades y costos
G1 G2 C
G1
C
G2
C
total
30 120 1295 5404 6699
50 100 2255 4254 6509
80 70 3845 2716.5 6561.5
100 50 5005 1816.5 6821.5
Despacho ptimo en un barraje
Ejemplo
Solucin por
multiplicadores de
Lagrange.
150
30 * 25 . 0
40 * 2 . 0
2 1
2
2
2
1
1
1
= +
= + =
= + =
G G
G
G
G
G
P P
P
P
C
P
P
C
CAGS
2 1 G G
Son tres ecuaciones, tres incgnitas.
Solucin: P
G1
= 61.11, P
G2
= 88.89,
= 52.22
Los despachos obtenidos, cumplen la restriccin:
2 , 1 , 125 20 = i P
Gi
Despacho ptimo en un barraje
Ejemplo
Costos alrededor del punto ptimo.
G
1
G
2
C
G1
C
G2
C
total
CAGS
60 90 2765 3714.5 6479.5
61.11 88.89 2822.8 3656.4 6479.2
62 88 2869.4 3612 6481.4
Despacho ptimo en un barraje
Calcular el despacho ptimo para diferentes escenarios de carga.
Cuando el despacho de uno de los generadores supera o esta por
debajo de sus valores limites, se genera, el valor mximo o
mnimo respectivamente. Lo anterior es valido siempre y cuando
la condicin determinante sea obligar a la generacin
Ejemplo. Para el caso anterior, calcular el despacho ptimo para
una generacin total de 40 MW.
CAGS
40
30 * 25 . 0
40 * 2 . 0
2 1
2
1
= +
= +
= +
G G
G
G
P P
P
P
Solucin: P
G1
= 0, P
G2
= 40, =40
Se redespacha como P
G1
= 20, P
G2
= 20
Despacho ptimo en un barraje
Tarea
Para el caso anterior, determinar el
despacho ptimo para un rango de carga
desde 40 hasta 250 MW, en incrementos
CAGS
desde 40 hasta 250 MW, en incrementos
de 5 MW. Suponer que siempre se debe
generar por lo menos el mnimo.
Despacho ptimo en un barraje
Supongamos dos
generadores que tienen
los siguientes costos:
100 * 30 * 125 . 0
120 * 40 * 1 . 0
1
2
2 2
1
2
1 1
+ + =
+ + =
P P C
P P C
G
G
Asimismo, alimentan una carga que tiene la siguiente
CAGS
Asimismo, alimentan una carga que tiene la siguiente
demanda diaria.
Demanda Horaria
00:00 06:00
06:00 18:00
18:00 23:59
0
50
100
150
200
250
00:00 04:48 09:36 14:24 19:12 00:00 04:48
Hora
M
W
P
mn
= 76
P
mx
= 220
Despacho ptimo en un barraje
Supongamos que el solo prender un generador,
ocasiona un gasto de 400. Es conveniente, desde un
punto de vista financiero, apagar uno de los
generadores, en los periodos de baja carga?
Obviamente ambos generadores deben estar en
CAGS
76
30 * 25 . 0
40 * 2 . 0
2 1
2
1
= +
= +
= +
G G
G
G
P P
P
P
Obviamente ambos generadores deben estar en
servicio, en el periodo de carga alta.
En el periodo de demanda mnima se tiene que:
56
20
2
1
=
=
G
G
P
P
Despacho ptimo en un barraje
De acuerdo a lo anterior, se tendran los dos
siguientes escenarios
G
1
G
2
C
G1
C
G2
Costo
Fijo
C
total
CAGS
0 76 0 3102 400 37624
20 56 960 2172 0 37584
Los 400 adicionales, solo se le asignan al escenario
en el cual G
1
= 0, ya que en este caso, es necesario
prender el generador para las siguientes 12 horas
con demanda total de 220 MW
Despacho ptimo teniendo en cuenta
prdidas
En este caso, las prdidas pueden afectar la forma en que
distribuimos la generacin entre varias plantas. El problema de
optimizacin es entonces minimizar:
=
=
k
i
Gi i
P C C
1
) (
CAGS
en todo momento, sujeto a satisfacer la demanda, incluyendo
prdidas de transmisin
) 1 ( 0
1
=
=
L
k
i
D Gi
P P P
k: nmero total de generadores
P
gi
: generacin de la i
esima
planta
P
D
: suma de las demandas en cada barra
P
L
: Prdidas totales del sistema de transmisin
Ecuacin de
balance de
Potencia
Despacho ptimo
Se plantea el Langragiano:
= =
=
k
i
k
i
L D Gi Gi i
P P P P C
1 1
) (
L
CAGS
Es posible demostrar que si se supone que el factor de
potencia de cada carga permanece constante, las
prdidas del sistema, se pueden expresar como una
funcin de la generacin de activa en cada planta.
) ,....., , , (
3 2 1 Gk G G G L L
P P P P P P =
Despacho ptimo
El despacho ptimo se logra cuando
k i
P
P
P
C
P
Gi
L
Gi
i
Gi
,...., 2 , 1 , 0 = =
L
Reorganizando
CAGS
Reorganizando
k i L
P
C
o
P
P
P
C
i
Gi
i
Gi
L
Gi
i
,...., 2 , 1 , ) 2 (
1
= =
=
|
|
\
|
|
|
\
|
=
Gi
L
i
P
P
L Donde
1
1
Factor de
Penalidad del
generador i
esimo
Despacho ptimo
El multiplicador Langragiano estar en
$/MWh, cuando el costo del combustible este
en $/hora. La ecuacin anterior, indica que el
costo mnimo de operacin se logra cuando el
costo incremental de operacin de cada planta
CAGS
costo incremental de operacin de cada planta
multiplicado por su factor de penalidad, sea
idntico para todos los generadores.
Las (k+1) variables (P
G1
, P
G2
,., P
Gk
, ) se
pueden calcular a partir de los k despachos
ptimos (2) y de la ecuacin de balance de
potencia (1).
Despacho ptimo
A la derivada parcial
Gi
L
P
P
Se le conoce con el nombre de Prdidas Incrementales
de Transmisin (ITL) del i generador. La ecuacin
CAGS
de Transmisin (ITL) del i
esimo
generador. La ecuacin
(2) tambin se puede escribir como
[ ] k i ITL IC
i i
,...., 2 , 1 ) ( 1 ) ( = =
A esta ecuacin se le conoce con el nombre de
Ecuacin Exacta de Coordinacin
Prdidas de Transmisin
Forma aproximada
para expresar las
prdidas como una
funcin de las
potencias activas de
Gn mn
k
m
k
n
Gm L
P B P P
= =
=
1 1
P
Gm
, P
Gn
= Potencia activa en el m-
CAGS
potencias activas de
los generadores:
Coeficientes-B
P
Gm
, P
Gn
= Potencia activa en el m-
avo o n-avo generador
B
mn
= Coeficiente de prdidas, los
cuales son constantes, bajo
determinadas condiciones de
generacin
Si P
G
,s estn en MW, B
mn
debe estar en MW
-1
. B
mn
= B
nm
Prdidas de Transmisin
La ecuacin de
prdidas se puede
escribir en forma
matricial, como:
=
=
k
ik
G
G
G
G
t
G L
B B B
B B B
B y
P
P
P
BP P P
2 22 21
12 11
2
1
M M
L
L
M
CAGS
kk k k Gk
B B B P
2 1
B es una matriz simtrica. Para un sistema de tres plantas se tendra:
1 3 31 3 2 23 2 1 21
2
3 33
2
2 22
2
1 11
2 2 2
G G G G G G G G G L
P P B P P B P P B P B P B P B P + + + + + =
Y
= =
Gn mn
k
m
k
n
Gm
Gi Gi
L
P B P
P P
P
1 1
Prdidas de Transmisin
=
=
+ +
k
j
Gj ij
Gi
L
Gi ij Gi
k
i m
m
Gi mn Gm Gn in
k
i n
n
Gi
Gi Gi
L
P B
P
P
P B P P B P P B P
P P
P
1
1 1
2
1 C
CAGS
i Gi i Gi i i
d P b P a C + + =
2
2
1
i Gi i
Gi
i
b P a
P
C
+ =
Recordando que
y
[ ] k i ITL IC
i i
,...., 2 , 1 ) ( 1 ) ( = =
( )
=
=
+ = +
= + +
k
i j
j
Gj ij i Gi ii i
k
j
Gj ij i Gi i
P B b P B a
P B b P a
1
1
2 2
2
Prdidas de Transmisin
k i
B
a
P B
b
P
ii
i
k
i j
j
Gj ij
i
Gi
....., , 2 , 1 ;
2
2 1
1
=
|
\
|
+
=
CAGS
La ecuacin anterior, junto con la ecuacin de balance de potencia, son la
base para calcular los despachos ptimos de los generadores. El
despacho del generador i es funcin del despacho de los generadores j. El
proceso es iterativo, iniciando en algn valor razonable. El proceso se
detiene cuando los despachos no se modifican entre iteracin e iteracin.
Despacho ptimo
1. Seleccionar un valor para:
2. Suponer P
Gi
= 0,
i=1,2,.,k
3. Solucionar (1) en forma
iterativa para los P
gi
,s
4. Calcular P
L
Gn mn
k k
Gm L
P B P P
=
k i
B
a
P B
b
P
ii
i
k
i j
j
Gj ij
i
Gi
....., , 2 , 1 ;
2
2 1
1
=
|
\
|
+
=
0
=
(1)
CAGS
4. Calcular P
L
5. Verificar si se satisface la
ecuacin de balance de
potencia
6. Si se satisface, detener
proceso, en caso
contrario, incrementar o
decrementar de acuerdo
a la siguiente expresin y
continuar con paso 3.
Gn mn
m n
Gm L
= = 1 1
<
=
L
k
i
D Gi
P P P
1
> |
\
|
< |
\
|
=
=
, 0
: , 0
1
1
L
k
i
D Gi
L
k
i
D Gi
P P P
P P P Si
Ejercicio 1
Para el sistema de dos barras de la figura, se transmiten 100 MW
de la planta 1 a la carga. Se tienen unas prdidas de 10 MW.
Calcular la generacin en cada planta y la potencia recibida por la
carga, cuando es de $25/MWh
El costo incremental es:
MWh P
P
C
G
G
/ $ 16 * 02 . 0
1
1
1
+ =
CAGS
MWh P
P
C
P
G
G
G
/ $ 20 * 04 . 0
2
2
2
1
+ =
GEN1
CARGA
GEN2
N2 N1
Ejercicio (cont)
P
G2
no tiene efectos en las prdidas, ya que su carga
esta en la misma barra.
B
22
= 0 y B
12
= 0 = B
21
2
P B P =
2
11
100 10 = B
CAGS
1 11 G L
P B P =
Entonces
1
11
11
001 . 0
100 10
=
=
MW B
B
20 2 2 04 . 0
16 2 2 02 . 0
1 21 2 22 2
2 12 1 11 1
= + +
= + +
G G G
G G G
P B P B P
P B P B P
G
1
G
2
= + +
=
k
j
Gj ij i Gi i
P B b P a
1
2
Ejercicio (Cont)
Al sustituir valores de B y
P
G1
= 128.57 MW
P
G2
= 125 MW
Las prdidas de transmisin sern:
CAGS
MW P
L
53 . 16 57 . 128 001 . 0
2
= =
La carga es:
P
D
= P
G1
+ P
G2
P
L
= 128.57 + 125 16.53 = 237.04
Ejercicio 2
Para el ejercicio anterior, calcular el ahorro debido a tener en cuenta
prdidas.
Ya se conoce el despacho teniendo en cuenta prdidas. Si no se
tienen en cuenta, se procede con la metodologa inicialmente vista:
20 * 04 . 0 16 * 02 . 0 + = + P P
CAGS
04 . 237 0001 . 0
20 * 04 . 0 16 * 02 . 0
2
1 2 1
2 1
+ = +
+ = +
G G G
G G
P P P
P P
Al solucionar el sistema anterior, da:
P
G1
= 275.18, P
G2
= 37.59
Ejercicio 2
El cambio en costo para el G
1
por pasar de la solucin anterior, a la nueva es:
hr
P P dP P
G G G G
/ $ 69 . 2937
16 01 . 0 ) 16 02 . 0 (
18 . 275
57 . 128
1
2
1
18 . 275
57 . 128
1 1
=
+ = +
CAGS
h
P P dP P
hr
G G G G
/ $ 43 . 2032
20 02 . 0 ) 20 04 . 0 (
/ $ 69 . 2937
59 . 37
125
2
2
2
59 . 37
125
2 2
=
+ = +
=
Ahorro: 2937.69 2032.43 = $905.26/hr
FORMULA PRDIDAS DE
TRANSMISIN
Mtodo preciso: Kron. Muy complicado
Presentaremos un mtodo simplificado
Sea el siguiente sistema:
CAGS
D
p
p
I
I
M
1
1
=
SM-17
SM-24
L36
N7
20 kV
L41
N33
20 kV
N10
20 kV
1
2
Id
Ip1
Id
El generador 1, provee
toda la potencia Id
necesaria.
Se define:
En forma similar para el generador 2
D
p
p
I
I
M
2
2
=
Factores de Factores de
Distribucin de Distribucin de
corriente corriente
FORMULA PRDIDAS DE
TRANSMISIN
Cuando ambos generadores estn en
funcionamiento
N33
Ig1
CAGS
SM-17
SM-24
L36
N7
20 kV
L41
N33
20 kV
N10
20 kV
1
2
Ig1
Ip
Id
Ig2
2 2 1 1
* *
G p G p p
I M I M I + =
FORMULA PRDIDAS DE
TRANSMISIN
Supuestos:
A) Todas las corrientes de carga tienen el mismo
ngulo de fase, con respecto a una referencia comn.
( )
i Di i i Di
I I =
CAGS
( )
i Di i i Di
I I =
i
Angulo del voltaje Angulo de la carga
De acuerdo a esto, tiene sentido que el ngulo vare
en un rango muy estrecho. Es decir, es casi constante
FORMULA PRDIDAS DE
TRANSMISIN
b) Se supone que la relacin X/R es igual para todas
las lneas y transformadores.
Dado que las relaciones X/R son iguales, las corrientes
por los elementos deben tener ngulos iguales a los de
las cargas.
CAGS
Con base en estas dos suposiciones se concluye
que I
P1
e I
D
tienen el mismo ngulo de fase y por lo
tanto los factores de distribucin de corriente (M
pi
)
son prcticamente nmeros reales.
FORMULA PRDIDAS DE
TRANSMISIN
Sea
2 2 2 1 1 1
= =
G G G G
I I y I I
Con base en las dos ecuaciones anteriores se tiene que
( ) ( )
2
2 2 2 1 1 1
2
2 2 2 1 1 1
2
sin sin cos cos
G P G P G P G P P
I M I M I M I M I + + + =
CAGS
( ) ( )
2 2 2 1 1 1 2 2 2 1 1 1
sin sin cos cos
G P G P G P G P P
I M I M I M I M I + + + =
) cos( 2
2 1 2 1 2 1
2
2
2
2
2
1
2
1
2
+ + =
G G P P G P G P P
I I M M I M I M I
Al expandir esta ecuacin y simplificar, tenemos
Si se tiene en cuenta que
2 2
2
2
1 1
1
1
cos 3 cos 3 V
P
I y
V
P
I
G
G
G
G
= =
FORMULA PRDIDAS DE
TRANSMISIN
Si R
p
es la resistencia de la rama p, las prdidas totales
del sistema sern:
=
P P L
R I P
2
3
Sustituyendo la ecuacin para I
P
se tiene que:
CAGS
( )
( )
+ =
p
P P
G
p
P P P
G G
p
P P
G
L
R M
V
P
R M M
V V
P P
R M
V
P
P
2
2
2
2
2
2
2
2
2 1
2 1 2 1
2 1 2 1
2
1
2
1
2
1
2
1
cos
cos cos
) cos( 2
cos
FORMULA PRDIDAS DE
TRANSMISIN
La ecuacin anterior es equivalente a
2
2 22 2 1 12
2
1 11
2
G G G G L
P B P P B P B P + + =
Donde
= R M B
2
1
CAGS
Donde
( )
( )
=
=
p
P P
p
P P P
p
P P
R M
V
B
R M M
V V
B
R M
V
B
2
2
2
2
2
2
22
2 1
2 1 2 1
2 1
12
2
1
2
1
2
1
11
cos
1
cos cos
) cos(
cos
1
22 12 11
, , B B B
A los coeficientes se les conoce como
coeficientes de prdidas o coeficientes B
FORMULA PRDIDAS DE
TRANSMISIN
Para el caso general de k plantas:
= =
=
k
m
k
n
Gn mn Gm L
P B P P
1 1
CAGS
Donde
=
p
P Pn Pm
n m n m
n m
mn
R M M
V V
B
cos cos
) cos(
Adems de los supuestos ya presentados, son
necesarios los siguientes supuestos adicionales, en
caso que se necesite que los coeficientes-B,
permanezcan constantes a diferentes despachos.
FORMULA PRDIDAS DE
TRANSMISIN
Todas las corrientes de carga, cambian en la misma proporcin con
CAGS
Todas las corrientes de carga, cambian en la misma proporcin con
respecto a la corriente total
Las magnitudes de voltaje permanecen constantes en todas las plantas
La relacin de potencia reactiva a la activa (o equivalentemente el
factor de potencia) permanece constante en todas las plantas
Los ngulos de los voltajes de fase en las plantas permanecen constantes.
EJEMPLO
La figura siguiente muestra
un sistema que consta de
dos plantas (1 y 2)
conectadas a las barras 1 y 2
respectivamente. La red
contiene 2 cargas y cuatro
lneas. La barra de referencia
I1
I2
a
1
c
b
0
2
4
Ref Bus
V=1/_0 pu
CAGS
lneas. La barra de referencia
tiene un voltaje 1 pu a cero
grados. Las corrientes e
impedancias son:
04 . 0 01 . 0 06 . 0 015 . 0
04 . 0 01 . 0 06 . 0 015 . 0
9 . 0 6 . 3 4 . 0 6 . 1
25 . 0 1 5 . 0 2
j Z j Z
j Z j Z
j I j I
j I j I
d b
c a
d b
c a
+ = + =
+ = + =
= =
= =
d
carga1
carga2
3
Ejemplo
Todas las corrientes
de carga, tienen una
relacin constante con
respecto a la corriente
total.
2174 . 0
15 . 1 6 . 4
25 . 0 1
7826 . 0
15 . 1 6 . 4
9 . 0 6 . 3
=
=
+
=
=
+
j
j
I I
I
j
j
I I
I
d c
c
d c
d
CAGS
total.
d c
7826 . 0 , 2174 . 0 , 2174 . 0 , 1
1 1 1 1
= = = =
d c b a
M M M M
7826 . 0 , 2174 . 0 , 7826 . 0 , 0
2 2 2 2
= = = =
d c b a
M M M M
Ejemplo
Generalmente se
debera correr un flujo
de carga. En este
ejemplo no es
necesario
o
o
9 . 4 051 . 1 09 . 0 048 . 1
) 06 . 0 015 . 0 )( 4 . 0 6 . 1 ( 1
05 . 6 066 . 1 1125 . 0 06 . 1
) 06 . 0 015 . 0 )( 5 . 0 2 ( 0 . 1
2
1
= + =
+ + =
= + =
+ + =
j
j j V
pu j
j j V
CAGS
Los ngulos de las corrientes de las plantas sern: ) , (
2 1 c b a
I I I I I + = =
( ) 1 0 cos cos
14
6 . 2
65 . 0
tan ; 14
2
5 . 0
tan
2 1
1
2
1
1
= =
=
|
\
|
= =
|
\
|
=
o
o o
Los factores de potencia de las plantas:
( ) 9393 . 0 14 05 . 6 cos
1
= + =
o o
pf
( ) 946 . 0 14 9 . 4 cos
2
= + =
o o
pf
Ejemplo
Los coeficientes de prdidas son:
02224 . 0
9393 . 0 * 066 . 1
7826 . 0 * 01 . 0 2174 . 0 * 01 . 0 2174 . 0 * 015 . 0 1 * 015 . 0
2 2
2 2 2 2
11
=
+ + +
= B
CAGS
00406 . 0
946 . 0 * 9393 . 0 * 051 . 1 * 066 . 1
7826 . 0 * 01 . 0 2174 . 0 * 01 . 0 015 . 0 * 7826 . 0 * ) 2174 . 0 (
01597 . 0
946 . 0 * 051 . 1
7826 . 0 * 01 . 0 2174 . 0 * 01 . 0 2174 . 0 * 01 . 0 7826 . 0 * 015 . 0
2 2
12
2 2
2 2 2 2
22
=
+ +
=
=
+ + +
=
B
B
Ejemplo
En ocasiones, la potencia de generacin est
en PU. En tal caso, los coeficientes-B, tambin
debern estar en pu. En base de 100 MVA:
02224 . 0
= B
CAGS
100
00406 . 0
100
01597 . 0
100
02224 . 0
12
22
11
=
=
=
B
B
B
FLUJO DE CARGA PTIMO
Planteamiento general. Se iniciar planteando una funcin objetivo
y restricciones de igualdad (ecuaciones de flujo de carga). Luego
se generalizar para incluir desigualdades (limites)
Minimizar la funcin objetivo escalar:
=
k
Gi i
P C C ) (
CAGS
=
i
Gi i
P C C ) (
Sujeto a las ecuaciones de flujo de carga. Barras PQ:
( )
( ) 0 sin
0 cos
1
1
= + +
= +
=
=
i j ij ij j
n
j
i i
i j ij ij j
n
j
i i
Y V V Q
Y V V P
FLUJO DE CARGA PTIMO
Barras PV
( ) 0 cos
1
= +
=
i j ij ij j
n
j
i i
Y V V P
En la barra I
th
P P P =
CAGS
Di Gi i
Di Gi i
Q Q Q
P P P
=
=
Donde P
Gi
y P
di
son las demandas en la barra i
Las ecuaciones para las barras PV y PQ se pueden expresar en
forma vectorial
FLUJO DE CARGA PTIMO
=
barra PV para P, en ecuacin
barra PQ para Q, en ecuacin
barra PQ para P, en ecuacin
) , ( y x f
= barra PQ para ,
barra PQ para ,
x
i
i
V
El vector de variables dependientes es:
CAGS
=
barra PV para ,
barra PQ para , x
i
i
El vector de variables dependientes es:
y el vector de variables independientes es:
=
p
u
V
P
Q
P
V
i
i
i
barra PV para ,
barra PV para ,
barra PQ para ,
barra PQ para ,
slack barra para ,
slack barra para ,
y
1
1
1
FLUJO DE CARGA PTIMO
El vector de variables independientes y se puede
particionar en un vector u de variables de control, las
cuales se modificarn hasta lograr un valor ptimo de
la funcin objetivo y un vector p, de parmetros fijos o
no controlables. Las variables controlables pueden ser
las magnitudes de voltaje en las barras PV y P sobre
CAGS
las magnitudes de voltaje en las barras PV y P
Gi
sobre
las barras con potencia de generacin.
FLUJO DE CARGA PTIMO
De acuerdo a lo
anterior, el problema
de optimizacin se
puede plantear como
( ) u x C
u
, min
CAGS
Sujeto a las
restricciones de
igualdad
( ) 0 , , = p u x f
FLUJO DE CARGA PTIMO
Para solucionar el problema de optimizacin,
se define la funcin del Langragiano como:
( ) ) , , ( ) , ( , , p u x f u x C p u x
T
= L
CAGS
( ) ) , , ( ) , ( , , p u x f u x C p u x = L
Donde es un vector de multiplicadores de
Lagrange de las mismas dimensiones que f(x,u,p)
Las condiciones necesarias para minimizar la funcin
de Lagrange sin restricciones son
FLUJO DE CARGA PTIMO
0
0
=
T
i
T
i
u
f
u
C
u
x
f
x
C
x
L
L
CAGS
( ) 0 , , = =
p u x f
L
Estas ecuaciones son no lineales y solo se
pueden resolver por medio de un proceso
iterativo. Un mtodo muy utilizado es el metodo
del gradiente o mtodo del menor paso
descendente
FLUJO DE CARGA OPTIMO
La tcnica para resolver las ecuaciones
anteriores, es ajustar el vector de control
u, para que pase una solucin valida
(conjunto de valores x que satisfacen las
ecuaciones de restriccin para unos
CAGS
ecuaciones de restriccin para unos
vectores u y p dados. Es decir la solucin
del flujo de carga) en la direccin del
paso descendente mas pequeo
(gradiente negativo), a un nuevo valor
con un valor de gradiente menor.
