SimSEE

SimSEE

Introducción a la Operación Óptima de

Sistemas de Generación de Energía
Eléctrica

1/102
 SimSEE

Organización.

Curso Básico. (5 clases de 2 hrs).

Conceptos generales y conocimiento de la herramienta SimSEE.
Curso Intermedio. (10 clases de 2 hrs).
Fundamentos de la operación óptima y profundización en el uso
de la plataforma SimSEE.
Webinar. Caso Belice (1 hr).
Ejemplo de aplicación de SimSEE para el análisis de la expansión
de la generación de Belice.

2/102
 SimSEE
Calendario Curso Básico De 9:00 a 12:00 (Ecuador; GMT-5)
Calendario Curso Intermedio De 9:00 a 12:00 (Ecuador; GMT-5)
Webinar Caso Belice De 9:00 a 11:00 (Ecuador; GMT-5)

20 020
20 2

3/102
 SimSEE

Curso Básico.

4/102
 SimSEE

Curso básico.

1) El sistema de generación. Conceptos básicos. Despacho óptimo en

un sistema muy simple. Operación óptima de sistemas dinámicos.
Ruben Chaer 02/03/2020 9:00 - 11:00
2) Ejemplo de aplicación – Caso Uruguay. Gonzalo Casaravilla.
04/03/2020 9:00 - 11:00
3) Ejemplo con Demanda, generadores térmicos. Cálculo del Costo
Marginal y Beneficio Marginal. Ingresos por energía y
disponibilidad. Ximena Caporale 06/03/2020 9:00 - 11:00
4) Se agrega al ejemplo modelo de la Red con 2 nodos y Arcos.
Congestión. Se agregan al ejemplo centrales hidroeléctricas con
Embalse. - Política de Operación. Ruben Chaer 09/03/2020 9:00 -
11:00
5) Se agregan al ejemplo generación Eólica y Solar. Ruben Chaer
11/03/2020 9:00 - 11:00
5/102
 SimSEE

Aprobación del curso Básico.

Durante el desarrollo del curso, se irá construyendo un sistema ficticio con

Centrales Térmicas, Hidráulicas, Eólica y Solar para el abastecimiento de
una demanda dada.
Cada alumno deberá tener instalado SimSEE e ir replicando los ejercicios
realizados en clase.
En la última sesión del curso se suministrará un formulario con preguntas
multiple opción que, para ser contestadas, se deberá realizar simulaciones
sobre el sistema construido.
Para la aprobación del curso se requiere un 80% de respuestas correctas.
6/102
 SimSEE

Curso Básico. - Clase 1

1) El sistema de generación. Conceptos básicos.

2) Despacho óptimo en un sistema muy simple.

3) Operación óptima de sistemas dinámicos.

Instructor: Ruben Chaer 02/03/2020 9:00 - 11:00

7/102
 SimSEE

SimSEE
https://simsee.org

Plataforma de Simulación de Sistemas de Energía Eléctrica

●
100% OOP.
●
Sala de Juego.
●
Actores.
●
Parámetros Dinámicos.
●
Monitores.
●
Modelos Estocásticos.
●
Gratuito.
●
Código Abierto.

8/102
 SimSEE

Sistema de Energía Eléctrica

9/102
 SimSEE

Posibles usos de SimSEE.

Optimización del uso de los recursos de generación.

Costo Maginal y Beneficio de Sustitución.
Cálculo del beneficio económico de un proyecto.
Cálculo de precios de equilibrio entre Demanda y Generación.
Figuras de Riesgo.
Planificación de Inversiones.
Optimización de la Operación de Mediano y Corto Plazo.
Optimización de las compras de combustibles.
Manejo de la incertidumbre hidrológica.
Volatilidad del precio del petróleo sobre los costos de generación.
Cálculo de la Potencia Firme necesaria en el sistema.

10/102
 SimSEE

El Sistema

Generadores
Demandas
Red eléctrica
Interconexiones

11/102
 SimSEE

El Simulador

Una Herramienta

12/102
 SimSEE

¿Por qué y para qué?

13/102
 SimSEE

Horizonte de tiempo y Paso de tiempo

Paso 1 Paso k

(t)

Horizonte

t_inicial t_final

Horizonte de Tiempo: ventana de tiempo en la que vamos a

realizar el estudio del sistema.
Paso de Tiempo: intervalo seleccionado para realizar la
integración de las ecuaciones del Sistema.
14/102
 SimSEE

Demanda - Potencia y Energía.

•

o
ing
ado

Dom
Sáb
Pico de la noche
Pico del mediodía
Valle de la madrugada

Serie diezminutal. Fuente ADME.

15/102
 SimSEE

Variación diezminutal

60

40 dP10m-dPH
dP10m

20
MW

0
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1

-20

-40

-60

16/102
 SimSEE

Reservas de corto plazo

Reserva Operativa. (2.1% de la Demanda; < 10 minutos).

Reserva 10 minutos. (2.1% de la Demanda; 10 a 20 minutos).
Reserva Fría. (3% de la Demanda; 20 minutos).
Seguimiento de la Demanda Horaria.

17/102
 SimSEE

Estacionalidad de la demanda.

Demanda diaria - Uruguay 2013

40
35
30
25
20
GWh

15
10
5
0
2013-01

2013-04

2013-07

2013-10

2013-12
18/102
 SimSEE

Evolución de largo plazo.

Demanda de energía eléctrica de Uruguay a nivel de generación.
12000 15,0%

10000 10,0%

[%] Crecimiento inter anual

8000 5,0%

6000 0,0%
GWh

4000 -5,0%
Crisis económica
2000 -10,0%
Restricciones
0 -15,0%

2005

2010
1965

1970

1975

1980

1985

1990

1995

2000

2015

2020
Fuente: ADME
19/102
 SimSEE

Generadores

En SimSEE, los Generadores son los Actores que producen energía.

Según su tecnología se agrupan en: Térmicos, Hidráulicos,
Eólicos, Solares, etc.
Los Generadores térmicos están formados por una o más unidades
turbina-generador (o moto-generador). Las turbinas pueden ser
del tipo turbo-vapor, donde se produce la expansión del vapor
generado en una caldera, o turbinas aero-derivativas como la de
los aviones tipo jet, en que la combustión se produce en la misma
turbina expandiéndose los gases de la combustión en su interior.
Las centrales térmicas más comunes utilizan como combustible
Carbón, FuelOil pesado, Gas Oil, Gas Natural, bio‑ combustibles.

20/102
 SimSEE

Sistema de transporte

Conectando a los Generadores y a las Demandas, tenemos el

Sistema de Transporte que en detalle está compuesto por el
conjunto de líneas, barras, transformadores etc. que hacen posible
el transporte de la energía eléctrica desde los generadores hasta
las demandas en condiciones de potencia, tensión y frecuencia
requeridos.
En SimSEE, bastará con definir barras “virtuales” de conexión
donde se resumen las grandes áreas de carga y “arcos” entre las
barras representando las capacidades de transporte y sus pérdidas
entre dichas áreas.

21/102
 SimSEE

Despacho y restricciones de Nodo

Resolver el DESPACHO significa decidir en cada instante de tiempo

qué máquinas serán las que generen y en qué cantidad, así como los
flujos por las interconexiones y arcos para cubrir las diferentes
demandas. El problema de despacho consiste en suministrar las
demandas al menor costo posible.
En la simulación del sistema, el despacho se resuelve en cada
intervalo de tiempo o paso de simulación.
En cada nodo de un sistema se debe cumplir en todo instante el
balance de Potencia. Es decir, la suma de las potencias inyectadas al
nodo debe ser CERO.

22/102
 SimSEE

Despacho y restricciones de Nodo

23/102
 SimSEE
Despacho y restricciones de Nodo

Podemos plantear el problema de despacho como:

mín ∑ c k (PG k )dt

(k )
@ ∑ PG k −∑ PD h =0
k k
Donde:
PGk es la potencia inyectada por el generador k
PDh es la potencia retirada por la demanda h en el nodo en
cuestión
ck(PGk) es el costo de generación del generador k expresado en
USD/h cuando entrega una potencia PGk y
dt es la duración del paso de tiempo en horas.

24/102
 SimSEE

Despacho y restricciones de Nodo

Entonces, el problema es minimizar el costo de generación en el

paso de tiempo, cumpliendo con las restricciones de nodo.

Este planteo es una simplificación del problema de despacho

general que veremos más adelante, en que el costo a minimizar es
el de generación en el horizonte de tiempo (conjunto de pasos) y
en el que además de las restricciones de nodo se deben de
verificar un conjunto mayor de restricciones como las impuestas
por los embalses, las líneas de transmisión etc.

25/102
 SimSEE

Sistema de transporte

Arco A
Arco B Nodo 2
Nodo 1

D2
G1 D1 G2

26/102
 SimSEE

●
Racionamiento
●
Falla
●
Déficit
●
Demanda insatisfecha
●
Costo de Falla

27/102
 SimSEE

Déficit o Falla

Cuando no es posible cubrir la demanda, no hay forma de satisfacer

una o más de las ecuaciones de balances de potencias. En ese caso se
produce un DÉFICIT o FALLA en el suministro de energía y habrá
demanda insatisfecha.

En la simulación del sistema, se agregan unos generadores ficticios

que llamamos MÁQUINAS DE FALLA, de forma que se puedan
satisfacer las ecuaciones de balance de potencia en todo instante.

Se puede saber en qué momentos se produce déficit y en qué

cantidad, observando el despacho de estas MÁQUINAS DE FALLA.

28/102
 SimSEE

Escalones de Falla

Se suele agregar más de una MÁQUINA DE FALLA,

asignándoles diferentes costos variables de operación.
En SimSEE, para cada Demanda se deben definir los escalones de
falla especificando su profundidad y costo.

29/102
 SimSEE

El DESPACHO del sistema en un instante dado es el estado del

sistema en lo que se refiere a cuáles máquinas se encuentran
prendidas y apagadas y qué potencia producen, junto con la
información de las potencias por las interconexiones que existan y
la potencia de la demanda.

El objetivo último de un sistema de generación de energía

eléctrica es el de suministrar la demanda al menor costo
posible. Cuando hablamos de costo de suministro estamos
hablando de costo en el largo plazo.

30/102
 SimSEE

Orden de mérito y Costo Marginal.

Sólo Costos Variables.

31/102
 SimSEE

Costo Futuro
+∞
CF p ( x k )=
ahora
⟨ ∑
k= k ahora
q
(k− k ahora )
ce k
⟩ r

El Costo Futuro (CF) es el Valor Esperado de la integral en el

tiempo, desde ahora hasta el infinito, del costo de cada Etapa
(ce).
El costo de etapa (ce) es la suma de los costos de compra de
combustibles más los costos de no suministro de la demanda más
costos de importaciones y menos los ingresos que se obtengan por
la exportaciones de energía hacia otros sistemas.
El factor q tiene en cuenta la tasa de actualización para el cálculo
del Valor Presente del costo de la operación futura.
El Valor esperado es sobre el conjunto de posibles realizaciones
de los procesos estocásticos involucrados (roturas de máquinas,
lluvia, viento, sol, etc.).
32/102
 SimSEE

Objetivo del OED = min(< CF >)

El Organismo Encargado del Despacho tiene como objetivo, en

todo momento, minimizar el valor esperado del CF.
Cuando en el sistema hay “reservorios”, la decisión de usar los
recursos almacenados en un instante de tiempo afectan los costos
futuros, dado que no dispondremos en el futuro de esos recursos
utilizados.
Cuando no hay reservorios, las decisiones tomadas en un tiempo
dado no afectan las posibilidades de tomar decisiones en el futuro.
En estas circunstancias, minimizar CF, es simplemente minimizar
el costo en cada etapa (o paso de tiempo).

33/102
 SimSEE

Ejemplo – Despacho de un sistema sencillo (sin reservorios).

Dos centrales:

G1: 80MW y 110USD/MWh

G2: 110MW y 200USD/MWh
Costo de Falla: 400USD/MWh.

34/102
 SimSEE
Ejemplo- Demanda 168 horas

35/102
 SimSEE

Curva de Oferta ( orden de mérito )

450

400
Costo Marginal [USD/MWh]

350

300

250

200 G1
G2
150 F

100

50

0
0 20 40 60 80 100 120 140 160 180 200 220 240

DEMANDA [MW]

36/102
 SimSEE

Ejemplo - Despacho horario resultante.

37/102
 SimSEE

Bandas Horarias, Bloques, Postes o Patamares

En la simulación de un SEE, en cada paso de tiempo hay que resolver

con qué generadores se suministra la demanda.
El paso de tiempo se debe hacer los suficientemente pequeño como
para poder suponer que la potencia de cada Actor es constante en el
paso de tiempo.
Para respetar el Balance Instantáneo de Potencia se debe cumplir que
la potencia media de cada Actor (Demandas y Generadores) sea una
buena representación de las potencias durante el paso de tiempo.
Cuanto más chico sea el paso de tiempo, mejor será la representación
del sistema y menores serán los errores. Por supuesto que a costa de
un mayor esfuerzo de cálculo.

38/102
Ejemplo, si el paso de tiempo fuese semanal. SimSEE

Demanda media = 155.74 MW

39/102
 SimSEE

Curva de Oferta ( orden de mérito )

450

400
Costo Marginal [USD/MWh]

350

300

250

200 G1
G2
150 F

100

50

0
0 20 40 60 80 100 120 140 160 180 200 220 240

DEMANDA [MW]

40/102
 SimSEE

Ejemplo, con Paso Semanal.

La energía de G1 + G2 es suficiente para abastecer la energía de la

semana, … pero falta POTENCIA.
Estamos ante un problema de POTENCIA y no de ENERGIA.
41/102
 SimSEE

Paso de tiempo dividido en POSTES o Bandas Horarias

como simplificación para mantener pasos de tiempo
grandes, respetando las restricciones de POTENCIA.

Se introduce el concepto de BANDA HORARIA, BLOQUE,

POSTE o PATAMAR como una subdivisión del Paso de Tiempo.
La idea es desordenar el tiempo dentro del Paso de Tiempo,
reagrupando las horas de acuerdo al requerimiento de potencia.
Se definen BANDAS HORARIAS, de diferente duración, en las
que el nivel de Demanda está bien representado.
Se debe especificar un POSTE “De Punta” y otro “De Valle” lo
suficientemente angostos como para representar adecuadamente
los requerimientos máximos y mínimos de potencia.

42/102
 SimSEE

Monótona de Carga y Postes Horarios

Poste 1 Poste 2 Poste 3

43/102
 SimSEE

Resolución por Postes

Resolvemos el despacho para cada poste considerando la potencia

media de la demanda en el poste, en las ecuaciones de balance de
potencia de los nodos.
Para el cálculo de costos tenemos en cuenta la duración de cada
poste.
Con esta definición en tres bandas dentro de un paso de tiempo, la
solución será prender las máquinas G1 y G2 durante todo el paso
de tiempo.
G1 estará “a pleno” durante todos los POSTES.
G2 estará “a pleno” solo en el POSTE 1 (punta).
En el POSTE 1, la suma de las potencias de G1 y G2 no es
suficiente y por tanto es necesario despachar Falla.

44/102
 SimSEE

Monótona de Carga y Postes Horarios

Poste 1 Poste 2 Poste 3

45/102
 SimSEE

Resolución por Postes

Desventajas:
Es una aproximación y como tal introduce errores o pérdida de
detalle en lo que estamos simulando.
No es posible tener en cuenta detalles como el costo de arranque y
parada de las centrales o la frecuencia con que éstas admiten ser
prendidas, dado que dentro del paso de tiempo se ha perdido la
secuencia horaria.
En SimSEE, se construyeron modelos de centrales que admiten
prendido y apagado por POSTE o prendido apagado por PASO
DE TIEMPO.

46/102
 SimSEE

Modelo de Central Térmica

El rendimiento de una central de generación térmica está determinado

por los rendimientos de sus turbinas, del generador eléctrico, del tipo
de combustible utilizado, de la temperatura del aire de alimentación,
etc.
Resumiendo todos los factores y conocido el costo del combustible,
representaremos el costo de generación como una función
dependiente de la potencia generada.

47/102
 SimSEE

Modelo de Generador Térmico Básico.

USD/h

cv inc

P máx MW

NPostes
C= ∑ P j cv inc Δ t j
j=1

48/102
 SimSEE

Modelo de Generador con Mínimo Técnico.

USD/h
cv inc

C mín

P mín P máx MW

NPostes
C= ∑ A j ( C mín +(P j −P mín )cv inc ) Δ t j
j=1

A j =P j > P mín
49/102
 SimSEE

Operación óptima de un Sistema Dinámico.

50/102
 SimSEE

Modelado y Simulación.

51/102
 SimSEE

Estado de un Sistema Dinámico

X = Vector de información que capta
todo lo relevante del pasado para
calcular el futuro si se conocen las
entradas de aquí en adelante.

x1
x
X (t )= 2
...
xn []
52/102
 SimSEE

Entradas de Control y Entradas No Controlables

r(t)
X(t)
u(t)

r(t) : Entradas que no podemos controlar. Por ej.: Lluvias.

u(t) : Entradas sobre las que podemos actuar para guiar el sistema
por donde nos convenga (entradas de control). Por ej.: Potencia
despachada en cada generador.

53/102
 SimSEE

Modelo del Sistema Dinámico.

x k +1=f (x k , u k , r k , k )
y k =g( x k , u k , r k , k )

x, es el vector de estado del sistema.

y, es el vector de variables observadas o salidas del sistema.
u, es un vector controlable de entradas del sistema.
r, es un vector no-controlable de entradas del sistema.
k, es número de paso de tiempo.

54/102
 SimSEE

Operador del Sistema.

El OPERADOR se modela como un bloque de control que,

observando el estado y las entradas del sistema, genera los valores
del vector de control u(t).

55/102
 SimSEE

Política de Operación.

u k = p(x k , r k , k )

56/102
 SimSEE

Costo Futuro partiendo de un estado dado en un paso de

tiempo dado y aplicando una Política de Operación dada.

+∞
CF p (x k )=
ahora
⟨ ∑
k =k ahora
q
(k−k ahora)
ce (x k ,r k ,u k , k )
⟩ {r k }

57/102
 SimSEE

Ejemplo: Sistema con una variable de estado

V k+1=V k +(q aporte−q turbinado−q vertido )Δ t

0≤V k +1≤V máx
P=h ρ g η q turbinado
58/102
 SimSEE

Ejemplo: Sistema sin variable de estado

Motor Gen

P= η⋅PCI⋅q

59/102
 SimSEE

Vinculación temporal de las decisiones.

El uso de recursos almacenados (agua) en

el presente produce un incremento en los
costos de operación del futuro. La
postergación del uso de un recurso
almacenado produce un incremento en
los costos del presente.

La Política Óptima es la que equilibra la

afectación de costos entre presente y
futuro.

Este es el resultado de Hamilton-Jacobi.

60/102
 SimSEE

Restricción dinámica.

Si no hubiera restricciones
para el traslado de recursos
en el tiempo el CMO sería
el mismo para todas las
horas.

61/102
 SimSEE

Representación de la
incertidumbre.

62/102
 SimSEE

Fuentes
de aleatoriedad

• Demanda eléctrica.
• Caudales de aportes hídricos.
• Velocidad del viento.
• Radiación solar.
• Precio de los mercados regionales.
• Precios de los combustibles.
• Disponibilidad de combustibles.
• Roturas fortuitas.

63/102
 SimSEE

Disponiblidad de
Generadores, líneas de
transmisión, etc.
Rotura ->

Disponible Indisponible

<- Reparación

●
Si no representamos el estado de Disponible/Indisponible, con pasos de tiempo
de simulación pequeños, el efecto de la Inercia del proceso de
Rotura/Reparación es subestimado.
●
Cada unidad representada agrega una variable de estado booleana.

64/102
 SimSEE

Modelos CEGH.

• Conservan histogramas de amplitudes.

• Conservan correlaciones.

Mundo Real

TNL (lentes)

Mundo
Gaussiano

65/102
 SimSEE

Modelos CEGH.

Espacio Gaussiano:
Sistema lineal multi-variable
alimentado por vector de
Ruido Blanco Gaussiano. NLT
h=nr −1 h=m−1
X k +1= A h X k −h + B h R k −h NLT
∑ ∑
h=0 h=0
NLT
Espacio Real
Permite definir
NLT
reducciones del estado.
Permite la incorporación
NLT
de información de
Pronósticos.
NLT

66/102
 SimSEE

Modelos CEGH.

Asimilación de
Pronósticos

67/102
 SimSEE
Operador Sin Pronósticos.

68/102
 SimSEE
Operador Con Pronósticos.

69/102
 SimSEE

Facilidad de integración de PRONÓSTICO.

h=nr −1 h=m−1
X k +1 = ∑ A h X k −h +S k + F k ∑ B h R k−h
h=0 h=0

Los sesgos (s) cambian la Guía del probabilidad 50% y los factores
de atenuación (f) regulan la inyección de ruido, permitiendo ir
desde un Pronóstico Deteminístico (ruido nulo) a la desaparición
del pronóstico (ruido histórico).

f 1 , k 0 ... 0
s1 , k

[]
S k = ...
sn , k [
F k = 0 f 2 , k ... 0
0 ... 0 f n , k ]
70/102
 SimSEE

Estado del proceso estocástico.

h=nr −1 h=m−1
X k +1 = ∑ A h X k −h + ∑ B h R k−h
h=0 h=0

La información del pasado, necesaria para calcular la salida del

sistema lineal en el espacio gaussiano, pasa a formar parte de El
Estado del Sistema.

X k , X k−1 , ... , X k−(n −1)

r

71/102
 SimSEE

Maldición de Belman.

)
,kk
,r
k
x
p(
=
k
u

Dim(u)×Dim( x)×Dim(r)×Dim(k )
72/102
 SimSEE

Reducción del Espacio de Estado en modelos CEGH

x1
x2 X
R
t
z= R( X )

A
u=PO z ( z , r , t )

X= A( z , w)
t
W es el ruido que permite poblar
el volumen de X que mapea
z en el z dado.

z1 73/102
 SimSEE

Tratamiento de un pronóstico en el modelado CEGH.

Gausianización.

Espacio REAL

P10
P50
P90

Espacio Gaussiano
P10
P50
x1
x2 P90

74/102
 SimSEE

Tratamiento del pronóstico en el espacio Gaussiano con

reducción.

P10 z= M R X
P50
x2
x1 P90 u=PO z ( z , r , t )
t R X= M A (t ) z+ B A (t) w

A
P10

P50
t
P90

z
z1
75/102
 SimSEE
http://vates.adme.com.uy

Pronos
Info. Pronos Pronos Eólica
Tiempo Hidro. Demanda Solar
Real

Pronos Despacho
VATES 168h
Tiempo
Continuo

Actualización cada una hora

WEB
PÚBLICA
ADME

76/102
 SimSEE

Optimización de la Operación.

El Jardín de las Delicias. EL BOSCO 1450-1516

77/102
 SimSEE

Política de uso de los recursos

¿Cuándo hay que elegir?

¿Qué es más económico?
Costo presente vs. Costo futuro.
Almacenes de energía.
Costos directo de tomar una decisión.
Costo por afectación del futuro.
Incertidumbre.
Valor Esperado y Riesgo.

Grabado de la serie “Los Caprichos”

Francisco de Goya, 1746 - 1828

78/102
 SimSEE

Sistema Dinámico, Operador y Política de Operación.

u k = p(x k , r k , k )

79/102
 SimSEE

Estado y
•

Poste Horario

Los Postes son un desorden del

tiempo.
Carece de sentido hablar de
estado por POSTE HORARIO.
El Estado será siempre por Paso
de Tiempo y nunca por Poste.

80/102
 SimSEE

Entradas no-controlables r(t)

•

Determinísticas y/o Aleatorias.

SIN ESTADO. Si son procesos aleatorios con memoria,
debemos identificarlos y representar su estado como parte del
Estado del Sistema.

81/102
 SimSEE

Costo Futuro y Costo de Etapa

•

∞
j−k
CF( x ,U k ,R k ,k)= ∑ q ⋅CE ( x j ,u j ,r j , j )
j=k
U k ={ uk ,uk+1 ,..... }
R k ={r k ,r k+1 ,..... }

x j+1=f ( x j ,u j ,r j , j ) Ecuación de transición o

Restricción dinámica

82/102
 SimSEE

Costo Futuro y Costo de Etapa

•

… recursión de Bellman

'
CF( x , U k , R k , k )=CE ( x , uk , r k , k ) +q⋅CF( x ,U k +1 , R k +1 , k +1 )
'
x = f ( x , uk , r k , k )

83/102
 SimSEE

Tasa de descuento
•

DurPaso/ DurAño
1
q= ( )
1+α
α=Tasas de descuento anual.
por ejemplo α=0.12 ; ( 12 % )

….. mmmhhhh ???

84/102
 SimSEE

Causalidad
•

Las decisiones del PRESENTE

pueden afectar el FUTURO
y No a la inversa???

85/102
 SimSEE

Minimizar el Costo Futuro

'
CF ( x , k )= ⟨mín { CE ( x ,u k ,r k , k ) +q⋅CF ( x , k+1 ) }⟩r
k
uk

86/102
 SimSEE

Minimizar el Costo Futuro

'
CF ( x , k )= ⟨mín { CE ( x ,u k ,r k , k ) +q⋅CF ( x , k+1 ) }⟩r
k
uk
'
x =f ( x , uk ,r k , k )
U k = { uk , uk +1 , . .. . . }= { uk , U k +1 }
R k = {r k , r k +1 , .. . . . }= { r k , R k +1 }

87/102
 SimSEE

Programación Dinámica Estocástica.

•

rk
xk xk 1

CE ''  
CF x '' , k  1

CE ' CF  x , k  1
'

t
uk

88/102
 SimSEE

Evolución del Estado

•

X k +1=f ( X k ,u k ,r k ,k )
Esta ecuación captura “la dinámica del sistema”.

NO-LINEAL y VARIANTE EN (t)

89/102
 SimSEE

Resolución Iterativa

90/102
 SimSEE

Valor del STOCK

Si pensamos que cada componente del estado x

representa un stock de un recurso (por ejemplo agua
embalsada), las derivadas de CF respecto de cada
variable pueden interpretarse como menos el valor que
le asignamos a una unidad de stock de esa variable.
Generalmente aumentar el stock de un recurso
disminuirá el CF por lo que estas derivadas son
negativas.

91/102
 SimSEE

Maldición de la dimensionalidad.

1
xk x1k 1

xk2
x 2k 1

t
kT  k  1T

N1 N 2 .... N nx
92/102
 SimSEE

Ejemplo CF(Bonete, EH)

93/102
 SimSEE

Ejemplo CF(Bonete, EH)

CF(1/1/2007-31/7/2009)

3 EH

1
70 72.4423 74.0951 75.4323 76.5859 77.6157 78.5546 79.423 80.2349 81
cota Bonete

94/102
 SimSEE

Ejemplo Valor del Agua

(Bonete, EH)

Valor del agua de Bonete (opt-sin sorteos) 1/1/2007

1200

1000

800
USD/MWh

600 Serie1

400

200

0
13
1

11

15

17

19

21

23

25

27

29

31

33

35

37

39

41

43

45

47

49
estado

95/102
 SimSEE

Optimización Con-Sorteos vs. Sin-Sorteos.

Valor del Agua de Bonete vs. máquinas del Sistema.

2500

2000

ConSorteos
SinSorteos
5ta
1500 6ta
CTR
PTI
SB
TGAA
1000 F1
F2
F3
F4

500

0
21

23

39

41
1

11

13

15

17

19

25

27

29

31

33

35

37

43

45

47

49
96/102
 SimSEE
Ej. Valor del Agua Bonete.
Estado hidrológico Super Seco (H1).

97/102
 SimSEE
Ej. Valor del Agua Bonete.
Estado hidrológico Super Húmedo ( H5 ).

98/102
 SimSEE

Pérdida de salto efectivo por caudal

erogado

Aproximación en SimSEE:
 
dh(QE) = caQE* QE + cbQE* QE^2

99/102
 SimSEE

Encadenamiento de centrales
hidroeléctricas.

a32 CH3
CH1 a12

a34

CH2 a24

CH4

100/102
 SimSEE

Control del vertimiento y erogado en

SimSEE.
SimSEE permite:
●
Fijar el vertimiento máximo en función del volumen
almacenado.
●
Especificar restricciones de caudal erogado mínimo.
●
Especificar restricciones de caudal vertido mínimo.

Erogado mínimo
Puede ser necesario para garantizar un nivel mínimo aguas abajo.
Navegabilidad – Toma de aguas.
Por control de crecidas – Protección de la presa.
Para mantener condiciones ambientales del río.

101/102
 SimSEE
FIN!!!. …

Muchas gracias por vuestra atención!.

102/102