Programación de la Generación

2. Despacho uninodal
Mayo 2020
 Restricción de
igualdad

Despacho económico uninodal

En un sistema uninodal, las
unidades generadoras y la
carga están conectadas a
un nodo. Se desprecian las
perdidas de transmisión.

El despacho económico es un N
Funci
problema de optimización, min Ft   Fi
ón
donde se minimiza el costo i 1
objeti
s.a.
total de operación del sistema, N
vo
sujeto a la ecuación de balance Pc   Pi Restricción
i 1 de balance
de potencia:
 Solución del despacho económico (DE)
por el método de Lagrange
Función de Lagrange La función de Lagrange
considera la función objetivo
L  FT   mas la restricción de balance
Donde: de potencia multiplicado por
N
N
el operador de Lagrange :
FT   Fi
i 1
;   Pc  P  0
i 1
i

N N
L   Fi   ( Pc   Pi )
i 1 i 1 La solución del DE calcula: ; Pi
Derivada de la función igual a cero y el costo operativo de cada unidad
dL Condición para dFi
0   0
dPi el mínimo dPi
N
dFi
  0 Ec. coordinación Pc   Pi
dPi i 1
Despacho económico uninodal
En un sistema uninodal, las
unidades generadoras y la
carga están conectadas a un
nodo.
Se desprecian las perdidas de
transmisión.

N
Funci
El despacho económico es un min FT   Fi ón
problema de optimización, s.a.
i 1
objeti
donde se minimiza el costo N
vo
total de operación del sistema, P P
i 1
i c
Restricción
de balance
sujeto a la ecuación de balance
de potencia: Pi min  Pi  Pi max Restricción
de capacidad
Condiciones de Kunh Tucker

Complementan las condiciones del Lagrangiano para incluir las

restricciones de desigualdad como términos adicionales.

d Fi
 ; Pimin  Pi  Pimax
d Pi
d Fi
; Pi  Pimax
d Pi

d Fi
; Pi  P i
min
d Pi
 Ejemplo1: DE uninodal
Un sistema de generación suministra a una
demanda de 52 MW, determine el despacho
económico:
 Costo Incremental
 Potencia de salida de cada unidad
 Costo de operación de cada unidad
y del sistema

Unidad H [Mbtu/h] Pmin [MW] Pmax [MW]

1 H1 = 0.234 P12 +2.112 P1 +112.8 12.63 21.05
2 H2 = 0.0022 P22 +8.71 P2 +69.91 11.29 19.50
3 H3 = 0.1032 P32 +6.119 P3 +79.61 12.16 20.23

Costo del combustible = 1.76 $/MBtu

Función costo F1 = 0.412 P12 + 3.72 P1 + 198

F ($/h) = 1.76 ($/MBtu)*H (MBtu/h) F2 = 0.004 P22 + 15.33 P2 + 123
F3 = 0.182 P32 + 10.77 P3 + 140
Ejemplo1(Cont.)
La solución del DE calcula:  , Pi dF1/dP1 = 0.824 P1 + 3.72 = λ
y el costo operativo de cada unidad dF2/dP2 = 0.008 P2 + 15.33 = λ
dFi dF3/dP3 = 0.364 P3 + 10.77 = λ
  0 P1 + P2 + P3 = 52
dPi
N La solución del sistema de ecuaciones:

P P
i 1
i c
λ ($/MWh) =
P1 (MW) =
15,52
14,33
Aplicar condiciones
de Kunh Tucker
P2 (MW) = 24,59 P2max = 19.5
Solución: P3 (MW) = 13,08

F ($/h) = 1.76 ($/MBtu)*H (MBtu/h)

Potencia de salida
λ ($/MWh) = 16,80 Costo operativo de las unidades
Luego: P1 (MW) = 15,89 F1 ($/h) = 361,5
F1 = 0.412 P12 + 3.72 P1 + 198 P2 (MW) = 19,5 F2 ($/h) = 423,4
F2 = 0.004 P22 + 15.33 P2 + 123 P3 (MW) = 16,61 F3 ($/h) = 369,1
F3 = 0.182 P32 + 10.77 P3 + 140
Costo operativo del sistema
FT ($/h) = 1154,1
Ejemplo1(Cont.) Gráficas Costo incremental
dF1/dP1 = 0.824 P1 + 3.72 = λ
dF2/dP2 = 0.008 P2 + 15.33 = λ
dF3/dP3 = 0.364 P3 + 10.77 = λ
P1 + P2 + P3 = 52

La solución del sistema de ecuaciones:

λ ( $/ MWh) = 15,52
P 1 ( MW) = 14,33
P2 ( MW) = 24,59
P3 ( MW) = 13,08

Aplicar condiciones
de Kunh Tucker
P2max = 19.5 MW

La solución final del DE

Potencia de salida Costo operativo de las unidades
P1 (MW) = 15,89 F1 ($/h) = 361,5
P2 (MW) = 19,5 F2 ($/h) = 423,4
P3 (MW) = 16,61 F3 ($/h) = 369,1

λ ($/MWh) = Costo operativo

16,80 del sistema
FT ($/h) = 1154,1
 Método del gradiente
Es un método iterativo de búsqueda de la Expandir el primer miembro de la
solución. ecuación en serie de Taylor de primer
 Iniciar con un valor de λ, orden, alrededor del punto de operación,
 Por extrapolación obtener un mejor valor de despreciar los términos de orden mayor,
λ.
 Continuar hasta que ΔPi sea menor que la  k ΔP (k) ΔP  k 
Δ    k
 df    
(k) ng
 dPi 
exactitud especificada.
   
i 1  d 

 dλ 
De la ecuación de costo incremental, obtener
Pi que es función de λ. ΔP  k 
  bi Δ  k
  k    k   Δ k 
d Fi ng
1
   2 c i Pi  bi Pi 
d Pi 2 ci 
i 1 2c i

De la ecuación de balance de potencia, obtener

Pi que es función de λ.
 k   bi
Ng

  bi  k
 PC   Pi(k)
Ng Ng
Pi  ΔP
f ( )   Pi    PC 2 ci i 1
i 1 i 1 2c i
Diagrama de Flujo - Método del gradiente
Inicio

LEER:
Pc, Fi, PL
λ˳, ΔP = 0.1

λ  k 1  λ  k   Δλ  k   k   bi
Pi 
2 ci

ΔP  k  Ng
Δ k   ng ΔP  k   PC   Pi(k)
1
 2c
i 1
i 1

SI
NO
ΔP < 0.1

Imprimir:

λ, Pi, Fi, FT

FIN
 Ejemplo2: DE (Método del Gradiente)
Un sistema de generación suministra a una
demanda de 52 MW, determine el despacho
económico por el método del gradiente:
 Costo Incremental
 Potencia de salida de cada unidad
 Costo de operación de cada unidad
y del sistema

dF1/dP1 = 0.824 P1 + 3.72 = λ

dF2/dP2 = 0.008 P2 + 15.33 = λ
F1 = 0.412 P1 + 3.72 P1 + 198; [Px=21.05 MW]
2
dF3/dP3 = 0.364 P3 + 10.77 = λ
F2 = 0.004 P2 + 15.33 P2 + 123; [Px=19.50 MW]
2
16  3.72
P1(1)   14.90
F3 = 0.182 P32 + 10.77 P3 + 140; [Px=20.23 MW] 0.824
1ra. Iteración: 16  15.33
Demanda = 52 MW P2(1)   83.75
0.008
λ0 = 16
16  10.77
ΔP = 0.1 P3(1)   14.37
0.364
 Ejm2: Método del Gradiente (Cont.)
16  3.72 (1)  (0)  Δ(0)
P1
(0)
  14.90
0.824
16  15.33 (1)  16  0.473  15.527
P2(0)   83.75
0.008 15.527  3.72
16  10.77 P1(1)   14.329
P3(0)   14.37 0.824
0.364 15.527  15.33
3 P2(1)   24.62  P2(1)  P2max  19.5
ΔP (0)  PC   Pi(0)  52 - (14.90  83.75  14.37) 0.008
i 1 15.527  10.77
P3(1)   13.07
ΔP (0)
 61.02 0.364
3
ΔP (1)
 PC   Pi(1)  52 - (14.329  19.5  13.07)
ΔP (0) i 1
Δ  (0)
1 ΔP (1)  5.10
 2c
i
5.1
 61.02 Δ(1)   1.288
Δ  (0)
 0.473 1 1
1

1

1 0
0.824 0.008 0.364 0.824 0.364
 Ejm2: Método del Gradiente (Cont.)
Costo operativo de cada unidad
( 2 )  (1)  Δ(1)
( 2 )  15.527  1.288  16.815 F1 = 0.412 P12 + 3.72 P1 + 198
F1 =361.13 $/h
F2 = 0.004 P22 + 15.33 P2 + 123
P2(1)  P2max  19.5MW
F2 =423.46 $/h
16.815  3.72
P1(2)   15.89 MW F3 = 0.182 P32 + 10.77 P3 + 140
0.824 F3 =369.10 $/h
16.815  10.77 Costo total de operación
P3(2)   16.61MW
0.364
FT = F1+ F2+ F3
3
FT = 361.13 +423.46+369.1
ΔP (2)
 PC   Pi(1)  52 - (15.89  19.5  16.61)
i 1 FT = 1153.7 $/h
ΔP (2)  0 Converge!!