𝑮𝟐 𝑮𝟏
N2 N1
TALLER
Juan Manuel Montoya Tovar
CEO World Energy Advisor
N3
[email protected] +51 949891256
Demanda
Introducción al Despacho Económico
con aplicación en GAMS
CONTENIDO
1. Tópicos de Despacho Económico
2. ¿Qué es un Modelo?
3. Flujo de Carga
4. ¿Qué es un flujo de carga óptimo?
BALANCE GENERACIÓN - DEMANDA
59.8 60 60.2 60.4
59.6
59.4 60.6
59.2 60.8
Frecuencia
1
¿Qué
centrales
DEMANDA GENERACIÓN deben operar
para cubrir la
Demanda?
Juan Manuel Montoya Tovar
CEO World Energy Advisor
RECURSOS PARA LA GENERACIÓN
ELÉCTRICA
G
N°orden D
E
Hidráulica
N E
E M
R
N° orden
Térmica
POOL
A
A
N
C N° orden
Eólica
I D
N° orden
O Solar A
N 2
¿Qué criterios tiene
para la operación?
Juan Manuel Montoya Tovar
CEO World Energy Advisor
RECURSOS PARA LA GENERACIÓN
ELÉCTRICA
US$/MWh
Diesel
Carbón
GN
Agua (río)
Agua (mar)
Sol
Viento
Vapor 3
¿Se puede Optimizar los
recursos (costos) para la
generación eléctrica?
Juan Manuel Montoya Tovar
CEO World Energy Advisor
CUESTIONAMIENTOS
¿Qué centrales deben
1 operar para cubrir la HIDRO TERMO EÓLICA SOLAR
Demanda?
¿Qué criterios tiene para EFICIE ECONO SEGURI
2 la operación?
NCIA MÍA DAD
A criterio se puede
¿Se puede Optimizar los seleccionar, pero el criterio
3 recursos (costos) para la
COMBU
STIBLE
AGUA
es subjetivo, no sería mejor
generación eléctrica? apoyarse en un Modelo?
Juan Manuel Montoya Tovar
CEO World Energy Advisor
¿QUÉ ES UN MODELO?
Los estudios realizados En la práctica, el Dadas las diferentes
sobre sistemas de sistema eléctrico se características entre
energía eléctrica se puede representar dichos fenómenos o
basan en la simulación todo tipo de eventos, previsión del despacho
de fenómenos reales fenómenos simultáneos de generación, no es
usando modelos que o de mayor horizonte posible definir un único
representan el mismo de tiempo, según la modelo capaz de dar
comportamiento que el conveniencia del una imagen real para
sistema real. estudio. todo el rango de
tiempos.
Juan Manuel Montoya Tovar
CEO World Energy Advisor
¿QUÉ ES UN MODELO?
Sistema Eléctrico
Es la representación
de la realidad en forma
simplificada para poder
analizarla
Juan Manuel Montoya Tovar
CEO World Energy Advisor
¿QUÉ ES UN MODELO?
Sistema Eléctrico
Es la representación
de la realidad en forma
simplificada para poder
analizarla
Es importante indicar que cuanto más cercano a la realidad se
desea modelar, más difícil es resolverlo, por ello es que la
mayoría de modelos tienden a ser de programación lineal
Juan Manuel Montoya Tovar
CEO World Energy Advisor
PROCESO
MODELO
Determinación de la estructura de un modelo y el hecho
de proporcionar los datos a dicho modelo
La tarea de
desarrollar un
modelo de un
sistema se puede
dividir en 2 sub tareas La estructura, fija las características que se deben de
considerar en el sistema eléctrico, mientras que, los
datos suministran los valores de dichas características
Juan Manuel Montoya Tovar
CEO World Energy Advisor
PROCESO
MODELO
Primero se debe tener un profundo conocimiento de las características de
funcionamiento y los parámetros que definen a cada uno de los equipos que
integran el sistema a simular
PROCESO
Caso contrario los resultados suelen carecer de validez por haber utilizado una
metodología y modelos que no son apropiados
Juan Manuel Montoya Tovar
CEO World Energy Advisor
TIPS
MODELO
Representa Selecciona
A la hora de diseñar el
cada componente del
Preguntate modelo hay que decidir
sistema eléctrico-
que parámetros es
energético de forma que se Analiza necesario incluir en el
puedan realizar análisis y
Qué expresiones matemáticas mismo para disponer de
estudios sobre dicho
describen las características de cada Cuáles, entre las distintas información suficiente
sistema
elemento del sistema que pueda ser representaciones o modelos, es el significativa y que otros
incluida como parte de la simulación que mejor describe el sistema que parámetros se pueden
deseas implementar obviar
Evalúa
Una vez que el modelo este desarrollado, el ajuste del mismo se realice mediante la
comparación de los resultados obtenidos con otros modelos similares. Es la única forma
de garantizar que el modelo esté bien diseñado
Juan Manuel Montoya Tovar
CEO World Energy Advisor
FLUJO DE CARGA
Sistema Eléctrico
¿En qué consiste?
𝑉1 𝜃1 Flujo de Carga
𝑉2 𝜃2
Consiste en calcular los valores de las tensiones en P,Q
𝑉3 𝜃3
módulo y ángulo, en todos los nudos de un sistema P,Q
eléctrico
𝑉7 𝜃7
P,Q P,Q
𝑉5 𝜃5
𝑉6 𝜃6 𝑉4 𝜃4
P,Q P,Q
A partir de estos valores, se determina el flujo de
potencia activa y reactiva en cada elemento del 4
sistema, y se obtiene las pérdidas en dichos elementos, ¿Se pudo Optimizar los
la potencia activa generada, entre otros. recursos (costos) para la
generación eléctrica?
Juan Manuel Montoya Tovar
CEO World Energy Advisor
FLUJO DE CARGA ÓPTIMO
La operación óptima del
sistema, en general,
comprende criterios de El objetivo principal en el
economía de operación, problema de suministro eléctrico
seguridad del sistema, o también llamado problema de
liberaciones óptimas de agua Despacho Económico, es reducir
al mínimo el costo total de la
en la generación
hidroeléctrica, etc. Todas estas
FOP generación (Costo de Producción
o Costo de Operación), al mismo
consideraciones pueden llevar tiempo que satisfacen la demanda
a requisitos en conflicto y de las cargas y las pérdidas en los
suele hacerse un compromiso enlaces de transmisión.
para la operación óptima del
sistema.
Juan Manuel Montoya Tovar
CEO World Energy Advisor
FLUJO DE CARGA ÓPTIMO
EN CONCLUSIÓN:
Tiene una función objetivo
El flujo de carga óptimo es aquel flujo
que obedece un objetivo y ante la
Sujeto a restricciones de igualdad y/o
búsqueda de éste fin, debe respetar las desigualdad
condiciones impuestas por el sistema.
Juan Manuel Montoya Tovar
CEO World Energy Advisor
MODELAMIENTO
1. Centrales Térmicas
2. Centrales Hidráulicas
3. Líneas de Transmisión
Flujo de Carga
Pérdidas en la transmisión
Modelamiento del Costo del Combustible
Costo ($)
𝑭𝑪 𝑷𝑮 = 𝑨 + 𝑩𝑷𝑮+C𝑷𝟐𝑮
𝑦 𝜀
De acuerdo a la base matemática del
𝐶𝐼𝑥 < 𝐶𝐼𝑦 despacho económico, en un sistema
𝑥 𝜀 puramente térmico, una unidad
𝜕𝐶𝑜𝑠𝑡𝑜
λ=
𝜕𝑃𝐺
= 𝑇𝑎𝑛𝜀 = 𝐶𝐼 participa en la formación del costo
𝑷𝑮(MW)
Pmín
marginal únicamente cuando opera
Pmáx
dentro de sus límites operativos (entre
Costo ($) su potencia mínima y máxima).
𝑦 𝜀 𝑭𝑪 𝑷𝑮 = 𝑨 + 𝑩𝑷𝑮
Cuando la unidad alcanza su límite
𝐶𝐼𝑥 = 𝐶𝐼𝑦 operativo máximo, se convierte en una
𝑥 𝜀
λ=
𝜕𝐶𝑜𝑠𝑡𝑜
= 𝑇𝑎𝑛𝜀 = 𝐶𝐼
unidad tomadora del costo marginal
𝜕𝑃𝐺
𝑷𝑮(MW)
resultante del sistema.
Pmín Pmáx
CI: Costo Incremental
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento del Embalse
Q [𝑚3 /𝑠]
caudal
Vmáx[M𝑚3 ]
S [𝑚3 /𝑠] Volumen máximo 𝑃ℎ ; Potencia de la Central hidráulica
vertimiento 𝜀: 𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙𝑎 ℎ𝑖𝑑𝑟𝑜
3.6: 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑑𝑒 𝑐𝑜𝑛𝑣𝑒𝑟𝑠𝑖ó𝑛
volumen mínimo
Vmín[M𝑚3 ]
𝜀[𝑀𝑊/𝑚3 ]
Central
Hidráulica
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento del Flujo de Carga
Basado en el método de Newton para encontrar potencias activas y potencias
reactivas y voltajes en cada nodo con el objetivo de minimizar el costo de
FOP AC operación
Flujo de carga que no toma en cuenta las pérdidas en las líneas y se centra en
determinar la potencia real. Está direccionado a proveedores que desean vender
FOP DC energía tan caro como sea posible y consumidores que deseen comprar energía tan
barato como sea posible
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento del Flujo de Carga
Barra 1 Barra 2
𝑹 + 𝒋𝑿
𝑉1 𝜃1 𝑉2 𝜃2
El flujo de carga de potencia activa y potencia reactiva es:
𝑉1 . 𝑉2 𝑉1
𝑃12 = 2 − cos(𝜃1 − 𝜃2 . 𝑅 + 𝑋. 𝑠𝑒𝑛(𝜃1 − 𝜃2 )
𝑅 + 𝑋2 𝑉2
𝑉1 . 𝑉2 𝑉1
𝑄12 = − cos(𝜃1 − 𝜃2 . 𝑋 − 𝑅. 𝑠𝑒𝑛(𝜃1 − 𝜃2 )
𝑅2 + 𝑋2 𝑉2
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento del Flujo de Carga en DC
Barra 1 Barra 2
𝑹 + 𝒋𝑿
𝑉1 𝜃1 𝑉2 𝜃2
Para un flujo de carga en DC, lo que se hace es considerar:
R ≈ 0
𝑉1 = 𝑉2 ≈ 1 𝑃𝑈
𝜃1 − 𝜃2 → 𝑚𝑢𝑦 𝑝𝑒𝑞𝑢𝑒ñ𝑜 Por ejemplo en el Modelo que
utiliza COES es un Modelo de flujo
de carga en DC con pérdidas
(están linealizadas)
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento del Flujo de Carga en DC
Haciendo: Haciendo: Haciendo:
R≈0 𝑉1 = 𝑉2 ≈ 1 𝑃𝑈 𝜃1 − 𝜃2 → 𝑚𝑢𝑦 𝑝𝑒𝑞𝑢𝑒ñ𝑜
𝑉1. 𝑉2 1 1
𝑃12 = 𝑠𝑒𝑛(𝜃1 − 𝜃2 ) 𝑃12 = 𝑠𝑒𝑛(𝜃1 − 𝜃2 ) 𝑃12 ≈ 𝜃 − 𝜃2
𝑋 𝑋 𝑋 1
𝑃12 ≈ 𝑌 𝜃1 − 𝜃2
𝑉1 .𝑉2 𝑉1 1
𝑄12 = − cos(𝜃1 − 𝜃2 𝑄12 = 1 − cos(𝜃1 − 𝜃2 𝑄12 ≈ 0
𝑋 𝑉2 𝑋
Flujo de potencia DC en cada línea
está definido en términos de dicha
expresión
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento del Flujo de Carga en DC
Proceso para el cálculo de un flujo de potencia en DC
Para empezar se parte Lo siguiente es formar
del cálculo de los la matriz admitancia y
componentes en Por tratar de ubicar la
Unidad como la dirección de los flujos
resistencia, reactancia,
potencia,voltajes, etc.
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento de las Pérdidas
en la Transmisión
Barra 1 Barra 2
𝑹 + 𝒋𝑿
𝑉1 𝜃1 𝑉2 𝜃2
Las pérdidas totales
𝑆 = ∆𝑉(𝐼)ҧ ∗ = 𝑃 + 𝑗𝑄…(1)
en la línea serán:
Donde: ∆𝑉 = 𝑉1 ∠𝜃1 − 𝑉2 ∠𝜃2 … (2)
En un flujo DC los voltajes 𝑉1 = 𝑉2 = 1 … (3)
en todas las barras es 1 PU:
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento de las Pérdidas
en la Transmisión en DC
y
Teniendo en cuenta la ecuación (2) y (3) su desarrollo sería el siguiente: (a,b)
b
∆𝑉 = cos𝜃1 + jsen𝜃1 − [cos𝜃2 + jsen𝜃2]
r
𝜃
a x
∆𝑉 = cos𝜃1 − 𝑐𝑜𝑠𝜃2 + 𝑗[sen𝜃1 − sen𝜃2 ]
𝑎 + 𝑗𝑏 = 𝑟∠𝜃 = 𝑟(𝑐𝑜𝑠𝜃 + 𝑗𝑠𝑒𝑛𝜃)
A B
La corriente viene dado por:
∗
∆𝑉
(𝐼)ҧ ∗ =
𝑅 + 𝑗𝑋
Considerando la ecuación (1) y las expresiones anteriores, las pérdidas totales son:
∗
∆𝑉 𝐴2 + 𝐵 2 𝐴 + 𝑗𝐵 . 𝐴 − 𝑗𝐵 = 𝐴2 + 𝐵 2
𝑆 = ∆𝑉. =
𝑅 + 𝑗𝑋 𝑅 + 𝑗𝑋
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento de las Pérdidas
en la Transmisión en DC
𝐴2 = 𝑐𝑜𝑠 2 𝜃1 + 𝑐𝑜𝑠 2 𝜃2 − 2cos𝜃1 cos𝜃2
𝐵 2 = 𝑠𝑒𝑛2 𝜃1 + 𝑠𝑒𝑛2 𝜃2 − 2sen𝜃1 sen𝜃2
𝐴2 + 𝐵 2 = 1 + 1 − 2cos𝜃1 cos𝜃2 − 2sen𝜃1 sen𝜃2
Por lo tanto las pérdidas totales: 𝑃𝑟𝑜𝑝𝑖𝑒𝑑𝑎𝑑:
2[1 − (cos𝜃1 cos𝜃2 + sen𝜃1 sen𝜃2)] cos 𝜃1 − 𝜃2 = cos𝜃1 cos𝜃2 + sen𝜃1 sen𝜃2
𝑆=
𝑅 + 𝑗𝑋
2[1 − (cos 𝜃1 − 𝜃2 ]
𝑆= 𝐴𝑟𝑡𝑖𝑓𝑖𝑐𝑖𝑜:
𝑅 + 𝑗𝑋
2[1 − (cos 𝜃1 − 𝜃2 ] Multiplicar a num. y den. por (𝑅 − 𝑗𝑋)
𝑆= . 𝑅 − 𝑗𝑋 … (4)
𝑅2 + 𝑋 2
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
Modelamiento de las Pérdidas
en la Transmisión en DC
De la ecuación (4), la pérdida activa y la pérdida reactiva viene dado por:
2[1 − cos 𝜃1 − 𝜃2 ] 2[1 − (cos 𝜃1 − 𝜃2 ]
𝑆= . 𝑅 + 𝑗 − .𝑋 𝐴𝑟𝑡𝑖𝑓𝑖𝑐𝑖𝑜:
𝑅2 + 𝑋2 𝑅2 + 𝑋 2
Multiplicar al denominador por 𝑅 2 + 𝑋 2
P: Pérdida Activa Q: Pérdida Reactiva
Aplicando un artificio sobre la Pérdida Activa: 𝑅
𝑋 𝑐𝑜𝑠𝜑 =
𝜑 𝑅 2 + 𝑋2
2[1 − cos 𝜃1 − 𝜃2 ] 𝑅
𝑃= . 𝑅
𝑅2 + 𝑋 2 𝑅2 + 𝑋 2
2[1 − cos 𝜃1 − 𝜃2 ]
𝑃= . 𝑐𝑜𝑠𝜑 𝑋≫𝑅
𝑅2 + 𝑋2 Las Pérdidas en cada línea
2[1 − cos 𝜃1 − 𝜃2 ] está definido en términos de
𝑃= . 𝑐𝑜𝑠𝜑 𝑃 = 𝑌. 2[1 − cos 𝜃1 − 𝜃2 ]𝑐𝑜𝑠𝜑 dicha expresión
𝑋
Centrales Centrales Flujo de Pérdidas en la
Térmicas Hidráulicas Carga transmisión
GAMS
1. Esquema general de la formulación
2. Estructura de un Modelo en GAMS
3. Ejemplo de un Modelo de un
Sistema de Despacho Hidro -
térmico con 6 barras
Esquema general de la Formulación
Esquema de Estructura de un Ejemplo de un
formulación Modelo en GAMS Modelo
Esquema general de la Formulación
Esquema de Estructura de un Ejemplo de un
formulación Modelo en GAMS Modelo
Esquema general de la Formulación
Esquema de Estructura de un Ejemplo de un
formulación Modelo en GAMS Modelo
Estructura de un Modelo en GAMS
Declaración de conjuntos (SETS)
Declaración de entrada (SCALARS, PARAMETERS, TABLE)
Resultados
Declaración de variables
Declaración de ecuaciones
Modelo
Resolución del problema de optimización
Esquema de Estructura de un Ejemplo de un
formulación Modelo en GAMS Modelo
Ejemplo de un Modelo de un Sistema de
Despacho Hidro térmico con 6 barras
H1 T1 T2 T3 Calcular el despacho óptimo de las centrales eléctricas y
luego conteste las siguientes interrogantes:
Barra 4 Barra 1 Barra 2 Barra 3
1. Valorice la energía producida por las generadoras eléctricas y
la energía consumida por la demanda, indique ¿la diferencia
Distribuidora de valorizaciones es igual a la diferencia entre energías
(pérdidas)?
2. En caso de estar en mantenimiento la central hidráulica
¿Varían mucho los costos marginales? y explique si
aumentara en 1MW la capacidad de la Térmica 2 ¿cuánto
sería el ahorro en el costo de operación?
3. Si la línea que enlaza las barras 3 y 6 se redujera su
capacidad a 130 MW ¿cuánto aumentaría el costo de
Barra 5 Barra 6 operación? ¿Qué sucede con los costos marginales cuando
ocurre la congestión?
Acería Minera
Esquema de Estructura de un Ejemplo de un
formulación Modelo en GAMS Modelo
Ejemplo de un Modelo de un Sistema de
Despacho Hidro térmico con 6 barras
Barra 4 Barra 5 Barra 6
Etapas
Cementera Distribuidora Mina
Unidades Térmicas 1 70.68 78.54 98.17
2 65.9 73.23 91.53
A B Pmáx Pmín 3
4
65.95
63.04
73.28
70.04
91.6
87.55
Líneas de Transmisión
1 200 16.669 80 0
5 63.44 70.49 88.11
2 200 9.333 120 0 6 68.72 76.36 95.44 resistencia reactancia Capacidad
Línea
3 200 12.727 150 0 7 68.34 75.94 94.92 (pu) (pu) (MW)
8 69.65 77.39 96.74 1-2 2.00057876 3.99967036 350
9 69.63 77.36 96.7
1-4 1.17560949 4.70608466 350
10 72.69 80.77 100.96
Central Hidráulica 11 72.34 80.38 100.48
2-3 0.77076766 3.84582299 350
12 75.26 83.62 104.53 3-6 1.92692897 9.61460651 350
Potencia Máxima 240 MW 13 73.83 82.03 102.54 4-5 1.00026702 1.99979046 350
Potencia Mínima 30 MW 14 71.56 79.51 99.39 5-6 1.00011275 2.99988055 350
Caudal de entrada 35 m3/s 15 71.39 79.32 99.15
16 71.71 79.68 99.6
Volumen inicial 250 MM3
17 71.44 79.37 99.22
Volumen máximo 450 MM3 18 70.91 78.79 98.48
Volumen mínimo 50 MM3 19 85.62 95.14 118.92
Volumen meta 125 MM3 20 89.2 99.11 123.88
Coeficiente de producción 3.846 21 90 100 125
Constante de conversión 22 86.7 96.34 120.42
3.6 23 77.22 85.8 107.25
caudal a volumen
24 72.42 80.46 100.58
Esquema de Estructura de un Ejemplo de un
formulación Modelo en GAMS Modelo
Muchas Gracias!
wea
world energy advisor
Juan Manuel Montoya Tovar
CEO World Energy Advisor
[email protected] +51 949891256
