Libro Generacion

Sistemas de Generación de
Energı́a
Por
Alejandro Garcés Ruiz.
Programa de Ingenierı́a Eléctrica

Universidad Tecnológica de Pereira
2
Universidad Tecnológica de Pereira.

3
Introducción
L OS sistemas de generación de energı́a desempeñan un papel cada
vez mas importante en el desarrollo económico de las naciones, por
tanto, el análisis de estos sistemas debe tener en cuenta el impacto
económico, social y ambiental. No obstante, para entender este impacto,
es necesario conocer de forma detallada el componente técnico implı́cito
en la operación y el planeamiento.
El propósito fundamental del presente libro, es introducir conceptos

relativos a los sistemas de generación de energı́a orientado a estudiantes
de últimos semestres de ingenierı́a. Para ello, es indispensable que el
estudiante esté familiarizado con conceptos de análisis de sistemas de
potencia y matemáticas básicas.
Una de las caracterı́sticas novedosas frente a la literatura clásica,

es el enfoque local; los diferente capı́tulos están enmarcado en el
contexto colombiano, en cuanto se estudian principalmente los sistemas
de generación termoeléctrica e hidroeléctrica, teniendo en cuenta dos
partes fundamentales en las que se divide el libro: aspectos técnicos de
los sistemas de generación y operación económica de los mismos.
En la primera parte se describe de forma general el funcionamiento

de los sistema de generación de energı́a más comunes en Colombia,
iniciando, en el capı́tulo 1 con una descripción detallada del sistema
eléctrico colombiano, teniendo en cuenta aspectos técnicos como la
capacidad instalada y las princiaples centrales de generación y aspectos
de tipo regulatorio que permiten entender el funcionamiento del sistema.
Posteriormente se presentan capı́tulos relacionados con el funcionamiento
de los sistemas térmicos e hidráulicos, haciendo especial enfásis en
estos últimos por su importancia estratégica para el paı́s. Se presenta
igualmente, un capı́tulo sobre hidrologı́a (capı́tulo 4) el cual permite
entender algunos aspectos relacionados con el carácter estocástico de ésta
variable y su repercusión con otras variables como la capacidad instalada
y la confiabilidad. El capı́tulo 5, presenta de forma general, las principales
caracterı́sticas de las turbinas en centrales hidráulicas mientras que en el
capı́tulo 6 se muestran algunos aspectos de la operación del sistema de
4
generación desde el punto de vista eléctrico.
En la segunda parte se presenta la operación economica de los sistemas

de generación de energı́a. Para ello se parte del capı́tulo 7 el cual es una
introducción a los conceptos de optimización, el objetivo fundamental
de éste capı́tulo es introducir conceptos fundamentales como función
objetivo, variables de descisión o restricciones, los cuales son básicos
para el modelamiento matemático de todo tipo de problemas no solo
en el campo de la ingenierı́a eléctrica. Posteriormente, se aplican estos
conceptos a los sistemas puramente térmicos y los sistemas hidrotérmicos
(capı́tulos 8 y 9). El capı́tulo 10 presenta las prinicipales caracterı́sticas
del mercado electrico colombiano desde el punto de vista matemático y
regulatorio. Finalmente, se presentan algunos clásicos de los sistemas
como la simulación de costos y el planeamiento de la expansión de la
generación.
Una de las caracterı́sticas únicas de este texto, es el uso de una

herramienta computacional altamente aplicada en la ingenierı́a y la
investigación de operaciones: el lenguaje GAMS. El uso de esta
herramienta en el campo de la investigación de operaciones es altamente
conocido, pero no existı́a un texto que aplicara directamente estos
conocimientos a la ingenierı́a eléctrica.
Mucho del material aquı́ presentado esta sustentado por trabajos

de investigación previos realizados en el grupo de investigación en
planeamiento de sistemas eléctricos y la maestrı́a en ingenierı́a eléctrica
de la Universidad Tecnológica de Pereira.
Contenido
I Aspectos técnicos de los sistemas de generación 1
1 El sistema eléctrico colombiano. 3
1.0 Estructura del sector eléctrico . . . . . . . . . . . . . . . . . . 5
1.1.1 Regulación y planeamiento . . . . . . . . . . . . . . . 6
1.1.2 Operación y administración . . . . . . . . . . . . . . . 9
1.1.3 Órganos consultores y asesores . . . . . . . . . . . . . 10
1.1.4 Fondos de apoyo . . . . . . . . . . . . . . . . . . . . . 11
1.1 El sistema colombiano . . . . . . . . . . . . . . . . . . . . . . 12
1.2.5 Centrales hidroeléctricas . . . . . . . . . . . . . . . . . 12
1.2.6 Plantas térmicas . . . . . . . . . . . . . . . . . . . . . . 14
1.2.7 Energı́a eólica . . . . . . . . . . . . . . . . . . . . . . . 15
1.2 Transacciones internacionales . . . . . . . . . . . . . . . . . . 15
1.3 La generación de energı́a como variable macroeconomica . . 15
5
6 CONTENIDO
2 Plantas Térmicas. 19
2.0 Funcionamiento de una central térmica . . . . . . . . . . . . 21
2.1 Partes constructivas . . . . . . . . . . . . . . . . . . . . . . . . 22
2.2 Principales ciclos termodinámicos . . . . . . . . . . . . . . . 24
2.3.8 Ciclo de Carnot . . . . . . . . . . . . . . . . . . . . . . 25
2.3.9 Ciclo Rankine . . . . . . . . . . . . . . . . . . . . . . . 27
2.3 Plantas de ciclo combinado. . . . . . . . . . . . . . . . . . . . 29
3 Centrales Hidroeléctricas. 33
3.0 Clasificación de las centrales hidroeléctricas . . . . . . . . . . 36
3.1.10 Central de pasada (a filo de agua) . . . . . . . . . . . 36
3.1.11 Central de embalse . . . . . . . . . . . . . . . . . . . . 42
3.1.12 Clasificación . . . . . . . . . . . . . . . . . . . . . . . 44
3.1 Aspectos técnicos de las centrales de embalse . . . . . . . . . 45
4 Hidrologı́a básica 51
4.0 Reservorios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
4.1.13 Algunas definciones . . . . . . . . . . . . . . . . . . . 52
4.1 Estudio de variables promedio . . . . . . . . . . . . . . . . . 54
4.2 Estudio de variables extremas . . . . . . . . . . . . . . . . . . 57
5 Turbinas hidráulicas 65
5.0 Caracterı́sticas generales . . . . . . . . . . . . . . . . . . . . . 66

CONTENIDO 7
5.1 Tipos de turbinas hidráulicas . . . . . . . . . . . . . . . . . . 67
5.2.14 Turbina pelton . . . . . . . . . . . . . . . . . . . . . . 68
5.2.15 Turbina francis . . . . . . . . . . . . . . . . . . . . . . 72
5.2.16 Turbina kaplan . . . . . . . . . . . . . . . . . . . . . . 75
5.2 Elección de la turbina en una central hidráulica . . . . . . . . 77
5.3 Control de velocidad (Gobernador) . . . . . . . . . . . . . . . 79
6 Aspectos operativos. 83
6.0 Curva de capacidad de un generador. . . . . . . . . . . . . . 84
6.1 Control de tensión. . . . . . . . . . . . . . . . . . . . . . . . . 88
6.2 Control automático de la generación (AGC) . . . . . . . . . . 90
6.3.17 Control primario de frecuencia . . . . . . . . . . . . . 92
6.3.18 Control secundario . . . . . . . . . . . . . . . . . . . 94
6.3.19 Control Terciario . . . . . . . . . . . . . . . . . . . . . 101
6.3.20 Reserva rodante . . . . . . . . . . . . . . . . . . . . . 102
6.3 Lecturas complementarias . . . . . . . . . . . . . . . . . . . . 103
II Aspectos económicos de los sistemas de generación 105
7 Introducción a la Optimización 107
7.0 Multiplicadores de Lagrange . . . . . . . . . . . . . . . . . . 109
7.1.21 Interpretación económica . . . . . . . . . . . . . . . . 111

8 CONTENIDO
7.1 Programación dinámica . . . . . . . . . . . . . . . . . . . . . 112
7.2.22 Algoritmo de solución . . . . . . . . . . . . . . . . . . 113
7.2 Programación lineal . . . . . . . . . . . . . . . . . . . . . . . . 116
7.3 Algoritmos de búsqueda aleatoria . . . . . . . . . . . . . . . . 118
8 Operación de sistemas térmicos. 121
8.0 Principales caracterı́sticas de las plantas térmicas . . . . . . . 122
8.1.23 Restricciones . . . . . . . . . . . . . . . . . . . . . . . 123
8.1.24 Costos de operación . . . . . . . . . . . . . . . . . . . 124
8.1 Despacho económico de plantas térmicas . . . . . . . . . . . 127
8.2.25 Restricciones. . . . . . . . . . . . . . . . . . . . . . . . 132
8.2.26 Red de transmisión . . . . . . . . . . . . . . . . . . . 137
8.2.27 Pérdidas de la red (Modelo DC) . . . . . . . . . . . . 141
8.2.28 Pérdidas de la red (Modelo AC) . . . . . . . . . . . . 144
8.2 Arranque y parada de unidades térmicas. . . . . . . . . . . . 150
8.3 Generación de un equivalente térmico. . . . . . . . . . . . . . 153
9 Operación de sistemas hidrotérmicos. 159
9.0 Modelo matemático lineal . . . . . . . . . . . . . . . . . . . . 164
9.1 Modelo no lineal de mediano plazo . . . . . . . . . . . . . . . 166
9.2.29 Costo del agua. . . . . . . . . . . . . . . . . . . . . . . 173
9.2.30 Múltiples embalses . . . . . . . . . . . . . . . . . . . 175

CONTENIDO 9
9.2 Modelo no lineal de corto plazo . . . . . . . . . . . . . . . . . 176
10 El mercado colombiano. 183
10.0 Modelo de mercado . . . . . . . . . . . . . . . . . . . . . . . 184
10.1.31 Concepto de competencia . . . . . . . . . . . . . . . 186
10.1.32 Bolsa de energı́a . . . . . . . . . . . . . . . . . . . . . 188
10.1.33 Pérdidas de potencia . . . . . . . . . . . . . . . . . . 192
10.1.34 Restricciones y reconciliaciones . . . . . . . . . . . . 193
10.1.35 Contratos bilaterales . . . . . . . . . . . . . . . . . . 199
10.1 Poder de mercado . . . . . . . . . . . . . . . . . . . . . . . . 202
10.2 Estrategias de los agentes . . . . . . . . . . . . . . . . . . . . 204
11 Simulación probabilı́stica de costos 207
11.0 Demanda . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 208
11.1 Generación . . . . . . . . . . . . . . . . . . . . . . . . . . . . 210
11.2.36 Salidas forzadas . . . . . . . . . . . . . . . . . . . . . 212
11.2.37 Confiabilidad . . . . . . . . . . . . . . . . . . . . . . 216
12 Planeamiento de la generación 219
12.0 Modelo estático . . . . . . . . . . . . . . . . . . . . . . . . . . 220
12.1.38 Función objetivo . . . . . . . . . . . . . . . . . . . . 221
12.1.39 Restricciones . . . . . . . . . . . . . . . . . . . . . . 222

10 CONTENIDO
A Métodos de Aforo. 227
A.0 Llenado de deposito . . . . . . . . . . . . . . . . . . . . . . . 227
A.1 Método del flotador . . . . . . . . . . . . . . . . . . . . . . . 228
A.2 Método del vertedero . . . . . . . . . . . . . . . . . . . . . . 231
B Formulación de modelos en Gams. 235
B.0 programación básica . . . . . . . . . . . . . . . . . . . . . . . 236
B.1 Solvers y modelos disponibles. . . . . . . . . . . . . . . . . . 239
B.2 manejo de ciclos y estructuras condicionales . . . . . . . . . 241
B.3.40 estructura condicional . . . . . . . . . . . . . . . . . . 241
B.3.41 estructura for . . . . . . . . . . . . . . . . . . . . . . . 242
B.3.42 Los operadores ord y card . . . . . . . . . . . . . . . 243
B.3.43 estructura loop . . . . . . . . . . . . . . . . . . . . . . 244
B.3 Manejo de archivos . . . . . . . . . . . . . . . . . . . . . . . . 245
B.4.44 Archivos gdx . . . . . . . . . . . . . . . . . . . . . . . 246

Parte I
Aspectos técnicos de los sistemas

de generación
1
El sistema eléctrico colombiano.
1
Contenido
Caracterı́sticas
Regulación y planeamiento
Operación y administración
Asesorı́a
Fondos de apoyo
3
4 CAPÍTULO 1. EL SISTEMA ELÉCTRICO COLOMBIANO.
L A estructura para el suministro de la energı́a eléctrica fue el

resultado de un prolongado proceso de intervención estatal, que se inició
prácticamente en 1928 con la expedición de la Ley 113 que declaró de
utilidad pública el aprovechamiento de la fuerza hidráulica. Desde
entonces funcionó de manera centralizada hasta las reformas efectuadas
en 1994. Durante el viejo esquema, las compañı́as estatales mantenı́an un
poder monopólico sobre un área determinada e integradas verticalmente,
prestaban los servicios de generación, transmisión y distribución. Este tipo
de monopolio sobre un área especı́fica, se debió al desarrollo regional que
presentaba el paı́s.
La primera central en
Colombia se construyó en
Bogotá en 1890 (8 años
Más tarde el sistema eléctrico colombiano se interconectó, y fue ası́
después que Tomas Alba
Edison instalara la
como nació ISA -Interconexión Eléctrica S.A-, permitiendo el intercambio
iluminación en la ciudad
de Nueva York)
de energı́a entre los sistemas regionales, con el fin de lograr el mejor
aprovechamiento de la capacidad energética de todo el sistema. ISA se
encargaba de la coordinación del suministro de electricidad, siguiendo
procesos de optimización, en donde se minimizaban los costos del sistema,
del planeamiento de la expansión del sistema de generación y transmisión
y, si era necesario, de la construcción y operación de las nuevas centrales
de generación.
Durante los años ochenta, el Sector Eléctrico Colombiano entró en

crisis, al igual que en la mayorı́a de paı́ses de América Latina. Esta
situación se debió especialmente al subsidio de tarifas y a la politización
de las empresas estatales, lo cuál generó un deterioro en el desempeño
de este sector. Al mismo tiempo, se desarrollaron grandes proyectos de
generación, con sobre-costos y atrasos considerables, lo que llevó a que
finalmente el sector se convirtiera en una gran carga para el Estado.
Por otro lado, en todo el mundo comenzó a ponerse en duda la eficacia

de los monopolios estatales para prestación de los servicios públicos,
iniciándose grandes reformas en algunos paı́ses tales como el Reino
Unido, Noruega y Chile.
Entre los principales cambios se destacan:

5
• Introducir competencia en el sector eléctrico
• Permitir la inversión privada, llegando al punto de privatizar las

compañı́as estatales.
• Eliminar la integración vertical, separando los negocios de

transmisión, distribución y generación.
• Dejar al estado solamente el papel de ente regulador
Ante los hechos anteriormente mencionados, a principios de los

años noventas se vio la necesidad en Colombia de modernizar el sector
eléctrico, abriéndolo a la participación privada, y siguiendo un esquema
similar a los paı́ses pioneros en este desarrollo, en especial el Reino
Unido. Esta reestructuración se realizó con las leyes 142 (Ley de
Servicios Públicos) y 143 (Ley Eléctrica) de 1994, las cuales definieron el
marco regulatorio para establecer las condiciones que permitieran que su
desarrollo estuviese determinado bajo la sana competencia. Estas leyes
crearon el Mercado Mayorista de Energı́a Eléctrica. La reglamentación de
este mercado fue desarrollada por la Comisión de Regulación de Energı́a
y Gas -CREG. Para este propósito, la Comisión se asesoró de consultores
nacionales e internacionales y con apoyo de las empresas del mismo
sector, promulgó las reglamentaciones básicas y puso en funcionamiento
el nuevo esquema a partir del 20 de julio de 1995.
1.1 Estructura del sector eléctrico
Aunque en teorı́a un esquema de mercado aumenta la eficiencia del

sector eléctrico, es función del estado definir reglas orientadas a fomentar
la competencia, eficiencia y seguridad energética del paı́s; ası́ mismo,
debe establecer los lineamientos en cuanto al planeamiento, operación
y administración del sistema. Para esto, se cuenta con una serie de
organismos de carácter estatal con funciones claramente definidas (ver
gráfico 1.1):
1.1.1
Los órganos de regulación y planeamiento determinan las politicas del
Regulación y
sector a mediano y largo plazo teniendo en cuenta aspectos económicos
planeamiento
como la competencia entre agentes del mercado y aspectos técnicos como
la calidad y seguridad en el suministro.
CREG: comisión de regulación de energı́a y gas. La conforman el ministro

de minas y energı́a, el ministro de hacienda y crédito público, el
director del departamento nacional de planeación, cinco (5) expertos
en asuntos energéticos nombrados por el presidente para un periodo
de cuatro (4) años y por el superintendente de servicios públicos
domiciliarios el cual asiste con voz pero sin voto.
La función de esta comisión consiste en regular la prestación
de los servicios públicos domiciliarios de energı́a eléctrica y gas
combustible, con el objeto de lograr que estos servicios se presten al
mayor número posible de ciudadanos, al menor costo posible para
los usuarios y con una adecuada remuneración a los prestadores del
servicio, garantizando calidad, cobertura y expansión.
UPME: unidad de planeamiento minero-energética Tiene por
objetivo planear en forma integral y permanente el desarrollo
y aprovechamiento de los recursos energéticos y mineros. Para
lograr el cumplimiento de estos objetivos, desarrolla, entre otras, las
siguientes actividades 1 :
• Establecer los requerimientos minero-energéticos de la
población y los agentes económicos del paı́s.
• Establecer la manera de satisfacer dichos requerimientos
teniendo en cuenta los recursos minero-energéticos existentes.
• Elaborar y actualizar el Plan Energético Nacional y el Plan de
Expansión del sector eléctrico.
• Desarrollar análisis económicos de las principales variables
sectoriales y evaluar el comportamiento e incidencia del sector
minero energético en la economı́a del paı́s.
1
Decreto Nº 255 enero 28 de 2004
7
Politicas sectoriales
Departamento Ministerio de Presidencia
MME
Planeación Hacienda
Nombra
S S PD
Asiste Control
CREG
(Regulación)
Asesora
Asesora
UPME CAPT CNO CND

Operación
Agentes
GENERACIÓN Despacha
Planes
Indicativos IS A
AS IC
INVERS IONIS TAS TRANS MIS IÓN (MercadoMEM )
PÚBLICOS Y
PRIVADOS DIS TRIBUCIÓN
US UARIOS
COMERCIALIZACIÓN
Vigila
Figura 1.1: Sistema eléctrico colombiano
• Evaluar la conveniencia económica y social del desarrollo de

fuentes y usos energéticos no convencionales.
• Evaluar la rentabilidad económica y social de las exportaciones
de los recursos mineros y energéticos.
• Prestar servicios técnicos de planeación y asesorı́a.
• Organizar, operar y mantener la base única de información
estadı́stica oficial del sector minero-energético.
• Establecer los indicadores de evaluación del sector
minero-energético, con el fin de elaborar informes que
cuantifiquen su gestión.
• Conceptuar sobre la viabilidad técnica y financiera de los

proyectos presentados, para ser financiados por el fondo
de apoyo financiero para la energización de las zonas no
interconectadas, FAZNI, y fondos especiales para energización
rural.
Las entidades del sector minero energético tienen la obligación

de hacer llegar toda la información necesaria para el adecuado
cumplimiento de las funciones de la unidad de planeación minero
energética.
SSPD: Superintendencia de servicios públicos domiciliarios Es un

organismo de carácter técnico. Creado por la Constitución de 1991
para ejercer el control, la inspección y la vigilancia de las entidades
prestadoras de servicios públicos domiciliarios. Desarrolla las
siguientes actividades:
• Publica las evaluaciones de gestión realizadas a los prestadores

y proporciona la información pertinente a quien la solicite.
• Da conceptos a las Comisiones de Regulación y Ministerios
que ası́ lo requieran en relación con los servicios públicos
domiciliarios.
• Certifica que la estratificación ha sido correcta, cuando se
trate de otorgar subsidios con los recursos nacionales y a
exigencia de la nación, para ello se basa en los resultados de
las estratificaciones enviadas por los municipios y distritos del
paı́s, en sus áreas urbana, centros poblados y rural.
• Establece los sistemas de información y contabilidad que deben
aplicar los prestadores de servicios públicos domiciliarios.
Define la información que las empresas deben proporcionar
sin costo al público y señala los valores que deben pagar las
personas por la información especial que pidan a las empresas
de servicios públicos, sino hay acuerdo entre el solicitante y la
empresa.
• Atiende los recursos de apelación que en subsidio interpongan
suscriptores y usuarios, una vez se haya resuelto el recurso de
reposición ante el prestador del servicio.
9
• Resuelve las apelaciones contra lo decidido por los Personeros

Municipales, por impugnaciones contra la elección de vocales
de control. Resuelve en segunda instancia los recursos
de reposición que interpongan los usuarios, en materia de
estratificación.
1.1.2 Operación
El planeamiento de corto plazo, denominado planeamiento de la y
operación, es realizado por una serie de entidades independientes, a saber: administración
CND: Centro nacional de despacho. Es la dependencia encargada de

la planeación, supervisión y control de la operación integrada de
los recursos de generación, interconexión y transmisión del sistema
interconectado nacional. Está igualmente encargado de dar las
instrucciones a los Centros Regionales de Despacho para coordinar
las maniobras de las instalaciones con el fin de tener una operación
segura, confiable y ceñida al reglamento de operación y a todos los
acuerdos del Consejo Nacional de Operación.
CRD: Centros Regionales de Despacho. Son centros de supervisión y
control de la operación de las redes, subestaciones y centrales de
generación localizadas en una misma región, cuya función es la
de coordinar la operación y maniobras de esas instalaciones, con
sujeción, en lo pertinente, a las instrucciones impartidas por el CND,
en desarrollo de las previsiones contenidas en el reglamento de
operación, con el fin de asegurar una operación segura y confiable
del sistema interconectado.
ASIC: Administrador del sistema de intercambios comerciales.
Dependencia, encargada del registro de los contratos de energı́a a
largo plazo; de la liquidación, facturación, cobro y pago del valor
de los actos o contratos de energı́a en la bolsa por generadores
y comercializadores; del mantenimiento de los sistemas de
información y programas de computación requeridos; y del
cumplimiento de las tareas necesarias para el funcionamiento
adecuado del Sistema de Intercambios Comerciales (SIC).
LAC: Liquidador y administrador de cuentas del sistema de transmisión

nacional - STN. Dependencia de Interconexión Eléctrica S.A. E.S.P
- ISA, que participa en la administración del MEM, encargada
de liquidar y facturar los cargos de uso de las redes del Sistema
Interconectado Nacional que le sean asignadas, de determinar el
ingreso regulado a los transportadores y de administrar las cuentas
que por concepto del uso de las redes se causen a los agentes del
mercado mayorista.
1.1.3 Órganos Éstos son los encargados de proponer polı́ticas de corto, mediano
consultores y y largo plazo para la operación del sistema y asesorar a los demás
asesores estamentos en la toma de decisiones:
CNO: Consejo Nacional de Operación. Organismo que tiene como

función principal acordar los aspectos técnicos para garantizar que la
operación integrada del sistema interconectado nacional sea segura,
confiable y económica, y ser el órgano ejecutor del reglamento de
operación.
Las decisiones del Consejo Nacional de Operación pueden ser
recurridas ante la Comisión de Regulación de Energı́a y Gas.
CAC: Comité asesor de comercialización Asiste a la CREG en los aspectos

de los comercializadores en el mercado de energı́a mayorista (MEM).
CAPT: Comité asesor de planeamiento de la transmisión. Asesora

a la UPME en la compatibilización de criterios, estrategias y
metodologı́as para el plan de expansión del sistema de transmisión
nacional. El CND no hace parte del CAPT pero asiste a sus reuniones
y le asesora en los aspectos de su competencia.
11
1.1.4 Fondos de
Dadas las caracterı́sticas de desarrollo desigual de las diferentes apoyo
regiones del paı́s, las cuales, impiden una competencia perfecta,
especialmente en zonas de poco desarrollo económico, se crean una serie
de fondos de apoyo:
FAZNI: Fondo de apoyo financiero para la energización de las zonas no

interconectadas. Es una cuenta especial del Ministerio de Minas
y Energı́a (MME) cuya función es financiar obras de desarrollo
energético en los departamentos contemplados en el artı́culo 309 de
la constitución (antiguas intendencias 2 y comisarias 3 ) ası́ mismo,
los departamentos de Chocó, Meta y Caquetá. La administración
de los recursos la realiza el Comité de Administración del fondo
(CAFAZNI) conformado por el director del departamento nacional
de planeamiento, el ministro de minas y energı́a y el ministro de
hacienda.
FOES: Fondo de energı́a social. Su objetivo es destinar hasta $40 por

cada kW/h de energı́a eléctrica que consuman usuarios ubicados en
zonas de difı́cil gestión y menor desarrollo. Se recauda mediante los
recursos de las rentas por congestión calculadas por el ASIC y de las
exportaciones de energı́a.
FAER: Fondo de apoyo financiero para la energización de zonas

rurales interconectadas. Financia proyectos de inversión para
la construcción e instalación de nueva infraestructura eléctrica,
que permita ampliar la cobertura de energı́a en las zonas rurales
interconectadas. El ASIC recauda $1 por cada kW/h despachado en
la bolsa del mercado de energı́a mayorista.
IPSE: Instituto de planificación y promoción de soluciones energéticas

para las zonas no interconectadas. Tiene como función ejecutar los
lineamientos y polı́ticas del ministerio de minas y energı́a en las
2
Arauca, Casanare, Putumayo, el Archipielago de San Andres, Providencia y Santa
Catalina
3
Amazonas, Guaviare, Guainia, Vaupes y Vichada
zonas no interconectadas; adelantar investigaciones y estudios sobre

las necesidades de estas zonas y la viabilidad técnica y financiera de
los proyectos.
1.2 El sistema colombiano
Colombia presenta una capacidad instalada de aproximadamente

13GW de los cuales el 35% es de origen termoelectrico mientras el restante
65% corresponde a la generación hidraulica4 . Las fuentes de energia
alternativa como la generación eólica presenta aún un porcentaje muy
bajo con respecto a la demanda del paı́s. No obstante, estas fuentes
de energı́a presentan una importancia estrategı́ca para el desarrollo del
sistema de generación, especialmente con la entrada en vigencia de nuevas
regulaciones ambientales que limiten el uso la generación térmica por
su emisión de gases de efecto invernadero y de grandes centrales de
generación por su alto impacto ambiental y social.
1.2.5 Centrales El mayor porcentaje de generación en Colombia corresponde a

hidroeléctricas centrales hidroeléctricas, esto reduce los costos de operación del sistema
pero lo hace vulnerable frente a condiciones de sequı́a, es allı́ donde
las plantas térmicas juegan un papel importante. La hidroelectricidad
tiene como caracterı́stica económica fundamental, los altos costos durante
la construcción que se ven compensados por los costos de operación.
En general una inversión en generación hidráulica presenta periodos de
retorno altos, los organismos de planeamiento deben entonces garantizar
reglas claras a largo plazo para motivar la inversión en este sector.
Algunas centrales están inmersas en cadenas hidrológicas como

4
No obstante, éste procentaje puede variar según la hidrologı́a anual pues la
capacidad de generación de la mayorı́a de emblases depende directamente de la cantidad
de agua almacenada, igualmente, las unidades de generación térmica pueden salir
temporalmente del sistema debido a mantenimientos o fallas inesperadas
13
la mostrada en la figura 1, estos sistemas permiten un mayor

aprovechamiento del recurso hı́drico pero requiere un mayor nivel de
análisis en el planeamiento de la operación.
Río Nare
El Peñol
1168,9Mm3
Guatapé
560 MW
Río Guatapé
Rebosadero
Río Biscocho
Rio Concepción
S an Lorenzo
J aguas 149,9 Mm3
170
Playas Rebosadero
56,1 Mm3
Playas
Río Calderas
200 MW Río Nare
Rio S an Carlos
Punchina
Otros afluentes 50,1 Mm3
Calderas
0,2 Mm3
S an Carlos
Calderas 1240MW
Río S amanáNorte
Figura 1.2: Candena hidráulica Nare - Gautapé - San Carlos
En la tabla 1.1 se muestran un resumen de las principales centrales

hidroeléctricas en Colombia.
Las centras hidráulicas presentan un acople temporal en la operación,

esto significa que una decisión operativa presente afecta los costos
operativos futuros. Esto se debe a la capacidad de almacenamiento de
agua de las grandes centrales de generación. Este tema será analizado en
Tabla 1.1: Principales centrales hidráulicas en Colombia

Central MW Embalse Ubicación
3
San Carlos 1240 Punchina (50, 1M m − 554m) Antioquia
Guavio 1150 Guavio (760.9M m3 − 110m)
Chivor 1000 Esmeralda (591M m3 − 768m)
Calderas Calderas (0, 2M m3 ) Antioquia
Guatapé 560 El Peñol (1168, 9M m3 ) Antioquia
Betania 540 Huila
3
Playas 200 Playas (56, 1M m Antioquia
Jaguas 170 San Lorenzo(149, 9M m3 ) Antioquia
forma detallada en capı́tulos posteriores.
1.2.6 Plantas La existencia en el paı́s de yacimientos de petroleo, gas natural y

térmicas carbón permite la generación de energı́a eléctrica por medio de centrales
térmicas; la generación térmica es totalmente complementaria a las
centrales hidraulicas ya que presenta bajo costo de instalación pero altos
costos de producción. En Colombia, estas centrales están ubicadas
especialmente en la costa atlántica debido a que esta zona no cuenta con
cadenas montañosas que permitan la construcción de un embalse para
generar energı́a.
Las principales centrales de generación termoeléctrica son presentadas

en la tabla 1.2
Tabla 1.2: Principales centrales termoeléctricas en Colombia

Central MW Tipo Ubicación
Tebsa 750 Barranquilla
TermoSierra 473 2 Turbinas de gas Antioquia
1 Turbina de vapor
Termo-Flores 400 Turbina de gas Barranquilla
Turbina de vapor
15
1.2.7 Energı́a
El parque eólico Jepirachi es un proyecto de generación eólica eólica
desarrollado por EPM y ubicado en el municipio de Uribia (Guajira), zona
perteneciente al resguardo indı́gena Wayuu.
El parque está conformado por 15 aerogeneradores Nordex N 60/250 +

1300kW interconectados mediante una red subterránea a 13, 8kV con
una capacidad instalada total de 19, 5M W . Cada aerogenerador está
compuesto por una hélice de 3 álabes con eje horizontal ubicado sobre
una torre de 60m. Los álabes forman un rotor de 60m de diámetro. El
generador de 4a6 polos está conectado directamente al eje de la turbina
con un peso aproximado de 5, 5 toneladas. El parque es una experiencia
piloto que hace parte del Programa general de investigaciones, proyectos y
actividades asociadas para el desarrollo de la energı́a eólica en Colombia,
con el cual se pretende adquirir conocimientos sobre esta energı́a, verificar
su desempeño y realizar la adaptación tecnológica a las caracterı́sticas
particulares del paı́s.
1.3 Transacciones internacionales
Colombia se encuentra interconectada con Ecuador y Venezuela como

se muestra en la tabla 1.3. Existen planes para una interconexión con
Panamá lo cual le representarı́a al paı́s la posibilidad de integrarse
con Centroamerica. En cuanto a las interconexiones existentes, la
más importante se presenta con Ecuador la cual hasta mayo de
2005 ha representado unas exportaciones de 3487, 58GW h equivalentes
a U SD$280 310 000 e importaciones de 114, 06GW h que equivalen a
U SD$3 410 000. Ası́ mismo se presentaron rentas por congestión de
U SD$158 380 000.
Tabla 1.3: Enlaces internacionales

Enlace Tensión (kV) Capacidad (MW)
Venezuela Cuestecitas Cuatricentenario 220 150
San Mateo El Corozo 230 150
La Fria Zulia 115 30
Ecuador Panamericana Tulcan 138 35
Jamondino Pomasqui 230 250
1.4 La generación de energı́a como variable

macroeconomica
La generación de energı́a eléctrica presenta una fuerte relación con el

crecimiento de la economı́a, prueba de ello es la forma del crecimiento de
la demanda con respecto al crecimiento del producto interno bruto (PIB)
como se muestra en la figura 1.3.
Figura 1.3: Crecimiento del PIB y de la demanda del paı́s

17
Es por ello, que el éxito en la predicción de la demanda de

energı́a eléctrica depende ostensiblemente del pronostico del crecimiento
económicos del paı́s. Ası́ mismo, otras variables energéticas como el
costo internacional del combustible, el crecimiento del gas domiciliario
entre otros, afecta significativamente el consumo de energı́a eléctrica.
Finalmente, la demanda de potencia pico, parece seguir una tendencia
lineal creciente por lo cual es posible predecirla por medio de una
regresión lineal simple.
Lecturas complementarias
1. Escuela de ingenierı́a de Antioquia. http : //f luidos.eia.edu.co/

obrashidraulicas/articulos/centraleshidroelectricasdecol/
2. Sistema Colombiano http : //www.upme.gov.co/siel/
3. ISA El sector eléctrico colombiano: origenes, evolución y retos (1882-1999)

Interconexión eléctrica SA, 1999.
4. Proyecto de generación eólica. http : //www.eeppm.com/epmcom/

contenido/acercade/inf raestructura/generacion/Jepirachi/index.htm.
Plantas Térmicas.
2
Contenido
Funcionamiento.
Principales ciclos.
Planta de vapor
Planta de ciclo combinado.
19
20 CAPÍTULO 2. PLANTAS TÉRMICAS.
U NA central térmica produce energı́a eléctrica a partir de la combustión

de carbón, fuel-oil o gas en una caldera diseñada al efecto. 1
Las plantas térmicas son de importancia estratégica para un sistema

como el colombiano, dada la necesidad de disponer de electricidad de
manera constante e independientemente de las variaciones del clima.
Mientras las hidroeléctricas necesitan de 4 a 7 años para su construcción,
una planta térmica se puede poner en operación en 1 o 2 años. Ası́ mismo,
la ubicación de la misma es menos restringida debido a la masificación de
los sistemas de gasoductos.
Existen varios tipos de plantas térmicas, pero todas se componen de

tres elementos básicos:
• Un elemento que transforma energı́a quı́mica (combustión o

quemado de combustible.)
• Un elemento que transforma energı́a mecánica (turbina o motor).
• Un elemento que transforma energı́a eléctrica (generador o

alternador).
La energı́a quı́mica se puede obtener de tres formas básicas:
• Un motor de combustión interna: se denomina ası́ porque dicha

combustión se realiza en el mismo motor y no en un elemento
independiente. Estos motores aprovechan la expansión de los gases
producidos por la combustión del diesel o búnker en la cámara de
un cilindro.
• Turbina de vapor: funciona al quemar el combustible en una caldera,

generando vapor que impulsa la turbina y genera movimiento a la
máquina sı́ncrona.
1
Las plantas nucleares pueden ser consideradas como plantas térmicas solo que en
este caso la caldera es remplazada por un reactor.
21
• Turbina de gas: en ésta, el fluido que produce el movimiento está

constituido por los gases de la combustión en cámaras especiales,
que elevados a temperatura y presión, mueven los álabes de la
turbina y la hacen girar velozmente. El movimiento giratorio del
eje de la turbina se trasmite al rotor de un generador que es el que se
encarga de producir la electricidad.
Las plantas del primer tipo presentan una limitada capacidad de

generación en potencia y disponibilidad horaria. Esta tecnologı́a es
principlamente utilizada en zonas no interconectadas.
Este capitulo se dedica a las plantas del segundo y tercer tipo por ser
las más comunes en Colombia.
2.1 Funcionamiento de una central térmica
El funcionamiento de todas las centrales térmicas es semejante. El

combustible se almacena en depósitos adyacentes, desde donde se
suministra a la caldera; una vez allı́, los quemadores provocan la
combustión, en un sistema convencional, el fluido de trabajo (agua)
recibe calor de la combustión en un proceso que idealmente se produce
a temperatura constante. El vapor producido entra con una presión
elevada en la turbina la cual hace girar al alternador. La turbina consta
de tres cuerpos: de alta, media y baja presión, unidos por un mismo
eje. En la primera parte (alta presión) los álabes se presentan en mayor
cantidad con un tamaño reducido. La turbina de media presión posee
una menor cantidad de álabes de mayor tamaño que los anteriores. Por
último, a baja presión se tienen álabes aún más grandes y en menor
cantidad. El objetivo de esta triple disposición es aprovechar al máximo
la fuerza del vapor, ya que este pierde presión progresivamente, por
lo cual los álabes de la turbina se hacen de mayor tamaño cuando se
pasa de un cuerpo a otro. Este vapor, antes de entrar en la turbina, es
cuidadosamente deshumidificado pues de no hacerlo las pequeñı́simas
gotas de agua en suspensión que serı́an lanzadas a gran velocidad contra
los álabes, actuarı́an como si fueran proyectiles y erosionarı́an las paletas

hasta dejarlas inservibles.
A la salida de la turbina, el vapor con una presión reducida, es enviado

a los condensadores. Allı́ es enfriado hasta que el agua llegue a estado
lı́quido. El fluido es nuevamente conducido a los tubos que tapizan las
paredes de la caldera continuando el ciclo.
Para minimizar los efectos de la combustión de carbón sobre el medio

ambiente, la central posee una chimenea de gran altura 2 , que dispersa
los contaminantes en las capas altas de la atmósfera. Igualmente existen
precipitadores que retienen buena parte de los gases contaminantes en su
interior.
2.2 Partes constructivas
En la figura 2.1 se muestra un esquema de una central térmica

convencional.
1. Entrada de combustible: El combustible (normalmente carbón) es

transportado desde la mina hasta el triturador en donde es cortado
en piezas de aproximadamente 15cm de diámetro. El carbón
triturado es llevado a la caldera (o a tanques de almacenamiento) a
través de una banda transportadora. Antes de ingresar a la caldera el
carbón es pulverizado para mejorar su aprovechamiento calorı́fero.
2. Caldera: El carbón pulverizado es quemado transfiriendo energı́a

calorı́fera al fluido de trabajo (agua) produciendo vapor a alta
presión.
3. Chimenea: Evacúa los gases de emisión a la atmósfera de forma que

se garantice una dispersión suficiente en la atmósfera.
2
Existen chimeneas de más de 300 metros
23
3
2
5
4 1
6 6
Figura 2.1: Esquema de una planta térmica
4. Equipo de reducción de emisiones: Antes que los gases de

emisión sean dispersos en la atmósfera se debe reducir su nivel
contaminante. Por medio de precipitadores electrostáticos se
eliminan las partı́culas de censa de los gases de combustión. El oxido
de azufre es eliminado por medio de una reacción quı́mica que lo
transforma en yeso.
5. Torre de refrigeración: Su función es enfriar el agua del circuito

de refrigeración. El aire recorre el interior de la torre en sentido
ascendente, enfriando el agua que cae en sentido contrario, en forma
de gotas de forma análoga a una ducha.
6. Bomba del circuito de refrigeración: Asegura la circulación del agua

de refrigeración entre el condensador y la fuente frı́a (Rı́o).
7. Condensador:El vapor que a cedido parte de su energı́a a la turbina
es dirigido al condensador, donde pasa de nuevo al estado de agua
lı́quida antes de incorporarse de nuevo al ciclo. Para condensar el
vapor se utiliza una fuente frı́a en un circuito cerrado por la torre de
refrigeración.
8. Turbina de vapor: El vapor producido por la caldera mueve los
álabes haciendo girar la turbina. De esta forma, la energı́a contenida
en el vapor se transforma en energı́a mecánica de rotación. Esta
energı́a es transferida al generador que la transforma en energı́a
eléctrica la cual por medio de un transformador elevador es llevada
a niveles de transmisión.
2.3 Principales ciclos termodinámicos
Las plantas térmicas se clasifican principalmente por el combustible

utilizado y los ciclos termodinámicos que utilizan, en general, una unidad
térmica puede utilizar varios ciclos termodinámicos prácticos. En este
campo, los esfuerzos en la investigación, se encaminan al aumento de
la eficiencia de la planta, con el fin de aprovechar más la capacidad
de generación de cada combustible. Para ello se puede crear ciclos
extremadamente complejos, que van mas allá del alcance de este texto;
un tema igualmente importante es el impacto ambiental de los diferentes
tipos de tecnologı́as, un ciclo mas eficiente es en general más limpio. A
pesar de la diversidad de tecnologı́as existentes, todas ellas se basan en
el ciclo termodinámico básico el cual está conformado por cuatro etapas
El fluido ingresa a la caldera como agua a presión alta y temperatura

baja, el calor se transfiere al fluido a una presión aproximadamente
constante. A la entrada de la turbina, el fluido (vapor) presenta una
temperatura y presión alta efectuando trabajo a medida que el fluido se
25
expande al pasar de la presión de la caldera a la del condensador. En la

etapa de condensación, el agua regresa a su estado liquido a una presión
casi constante transfiriendo calor al exterior del sistema. Finalmente, el
agua es extraı́da del condensador por medio de una bomba que alimenta
la caldera para iniciar el ciclo nuevamente.
Caldera
Turbina
4
Generador
2
Bomba Condensador
Figura 2.2: Esquema del ciclo termodinámico
2.3.8 Ciclo de
El ciclo de Carnot es un proceso cı́clico reversible que utiliza un gas Carnot
ideal y que consta de dos procesos isotérmicos3 y dos adiabáticos4 como
se muestra en la figura 2.3:
• Expansión isotérmica (AB): A temperatura TH el gas absorbe calor.
• Expansión adiabática (BC): El gas pasa de la temperatura TH a TL .

3
a temperatura constante
4
aislados térmicamente
• Compresión isotérmica (CD). A temperatura TL el gas libera parte

del calor.
• Compresión adiabática (DA): El sistema pasa nuevamente de la
temperatura TL a TH .
T P
A
TH A B PH
B TH
TL D TL
D C
C
S V
Figura 2.3: El ciclo de Carnot
Teóricamente, no puede existir un ciclo termodinámico más eficiente

que el ciclo de Carnot, no obstante, en la practica la eficiencia de un planta
real es inferior debido a que no hay procesos totalmente reversibles5 ,
siempre existirán pérdidas por fricción.
La eficiencia del ciclo puede ser expresada en términos de la

Nicolas Léonard Sadi
temperatura máxima TH y mı́nima TL , ambas en Kelvin mediante la
Carnot (1796-1832) en 1824
publico un trabajo
expresión:
denominado ’Reflexions
sur la puissance motrice
du feu et sur les machines
propres á développer cette
TL
puissance’ en donde se
mostraba el ciclo que lleva ηCnt = 1− · 100% (2.1)
su nombre. La principal TH
motivación de Carnot para
estudiar los principios de
la termodinámica fue el
atraso en que se En un ciclo real, además de la eficiencia del cliclo ηCnt existe una
encontraba Francia con
respecto al Reino Unido eficiencia ηB asociada al dispositivo de calentamiento y el combustible
tras las guerras
napoleonicas utilizado. En general la eficiencia total ηg está dada por el producto de
las dos:
5
Se presentan cambios en la entropı́a
27
ηg = ηCnt · ηB (2.2)
Existen técnicas para aumentar la eficiencia del ciclo en detrimento de

la eficiencia ηB por lo cual es importante tomar como referencia ηg . Una
mejor medida para determinar el rendimiento de un ciclo real es la relación
de eficiencia la cual indica qué tanto se aleja un ciclo real del ciclo teórico:
ηg(real)
ηEf = (2.3)
ηg(teorico)
2.3.9 Ciclo
En la figura 2.4 se muestra el ciclo de Rankine simple en donde la Rankine
numeración de cada etapa corresponde a la gráfica 2.2.
PB
T
TB
3
PA
TB
2
TA
1 4
Figura 2.4: Ciclo de Rankine

la eficiencia de este ciclo puede ser determinada en función de las

temperaturas medias (TA y TB ):

TA
ηCRk = 1− · 100% (2.4)
TB
El rendimiento del ciclo Rankine puede ser aumentado de varias

formas:
Aumentar el vacı́o en el condensador: Mientras más baja sea la presión

de salida de la turbina menor es la temperatura de saturación del
vapor en el condensador TA y por ende la eficiencia del ciclo. La
única limitación frente a esta mejora son los costos adicionales en las
torres de enfriamiento.
Aumentar la temperatura TB : Con esto la temperatura media TB también
aumenta. La principal limitación frente a esta mejora son los
materiales de construcción (ya que por encima de 565o C es necesario
introducir aceros austenitı́cos).
Aumentar la presión inicial de vapor: De esta forma se obtiene una
mayor temperatura de saturación incrementando la temperatura
media.
Existen algunas variantes del ciclo Rankine para aumentar la eficiencia:
Precalentamiento: Consiste en utilizar vapor para precalentar el agua a

la entrada de la caldera. Un proceso de precalentamiento permite
aumentar la eficiencia del ciclo al aproximar la temperatura media
TB al valor de la temperatura máxima TB . Pueden ser utilizadas
varias etapas de precalentamiento, no obstante, a medida que
aumenta el número de precalentadores disminuya la eficiencia
adicional por lo que convencionalmente se utilizan pocas etapas.
Los ciclos con precalentamiento no solo incrementan la eficiencia
sino que racionalizan las caracteristicas necesarias de la turbina y
disminuye los problemas de humedad excesiva en algunas partes
del ciclo.
29
Recalentador : En este ciclo dos turbinas operan en serie. La primera

recibe el vapor de la caldera a presión elevada; a la salida de la
primera turbina, el vapor es reingresado a la caldera y es recalentado
para pasar a través de la segunda turbina (de baja presión). Las
etapas de recalentamiento disminuyen la humedad a la salida de la
turbina aunque puede presentar desventajas desde el punto de vista
de la confiabilidad.
Regeneración: En este ciclo el fluido de trabajo es calentado antes de

ingresar al condensador por una parte del vapor que sale de la
caldera, esto puede aumentar la eficiencia del ciclo en la medida en
que establesca una adecuada relación costo beneficio.
2.4 Plantas de ciclo combinado.
Las plantas de ciclo combinado presentan una eficiencia superior a las

demás plantas térmicas6 . El ciclo combinado requiere la co-existencia de
dos ciclos termodinámicos en un mismo sistema. En uno el fluido de
trabajo es el agua mientras que en el otro es un gas combustible, esto
supone la utilización de una turbina de vapor y una turbina de gas. Los
gases de escape a altas temperaturas de la turbina de gas son utilizados
por el ciclo de vapor para generar energı́a.
1. Estación de regulación y medida (ERM): El gas natural llega a la

central por medio de un gasoducto, y pasa a traves de la estación de
regulación y medida donde se eliminan las eventuales impurezas del
gas y se determina la cantidad de gas a utilizar y su poder calorı́fico,
igualmente se regula la presión de trabajo.
2. Casa de filtro: El aire aspirado debe ser sometido a un proceso de

limpiado para eliminar las partı́culas que transporta.
6
Se puede obtener un rendimiento 50% superior al rendimiento de una central de ciclo
único.
4
7 3
2
1
Figura 2.5: Esquema de un ciclo combinado
3. Turbina de gas: El compresor de la turbina aspira aire y lo comprime

para inyectarlo en la cámara de combustión donde participa en la
combustión del gas natural. Los gases producidos se encuentran
a alta presión (30 bar) y temperatura (1300ºC) y cede parte de su
energı́a a la turbina que a su vez la transmite a la máquina sı́ncrona.
4. Caldera de recuperación: Los gases de escape de la turbina

transmiten su energı́a al agua que recorre los tubos situados en la
caldera, transformandola en vapor el cual es dirigido a la turbina de
vapor.
5. Turbina de vapor: El vapor producido por la caldera mueve los

álabes haciendo girar la turbina.
6. condensador: El vapor que ha cedido su energı́a a la turbina es

dirigido al condensador, donde pasa de nuevo al estado liquido para
incoroporarse nuevamente al ciclo. Para condensar se utiliza el agua
31
de un rio ( o el mar) o un circuito cerrado mediante una torre de

refrigeración.
7. Torre de refrigeración: El agua de refrigeración es bombeada a lo

mas alto de la torre para posteriormente ser enfriada por el aire que
circula en sentido ascendente.
1. R.W. Haywood. Análisis termodinámico de plantas eléctricas. Editorial

limusa.1986.
2. Frederick T. Morse. Centrales eléctricas. Editorial Continental
3. VY. Rizhkin Centrales electricas (2 tomos). Editorial Mir.
4. W, Desmond E. Advanced thermodynamics for engineers Editorial John

Wiley and Sons, 1997
5. V.N, Hendrick C, Zeransky. Basic engineering thermodynamics

Editorial McGrae-Hill.
Centrales Hidroeléctricas.
3
Contenido
Caracterı́sticas básicas.
Clasificación
Tipos de presas.
Obras civiles
33
34 CAPÍTULO 3. CENTRALES HIDROELÉCTRICAS.
U NA central hidroeléctrica es aquella en la que la energı́a potencial

del agua almacenada en un embalse se transforma en energı́a cinética
para posteriormente transformarse en energı́a eléctrica por medio de un
sistema turbina-generador.
Normalmente, las centrales hidroeléctricas se construyen en los cauces

de los rı́os, creando un embalse para retener el agua que es conducida por
medio de una tuberı́a de presión hasta la casa de máquinas. La capacidad
de generación de energı́a de una central depende principalmente del
caudal turbinado (Q), de la caı́da neta Hn y de la eficiencia del grupo
turbina generador (η).
La potencia que puede ser transformada a partir de un diferencial de

masa en la toma de la presa esta dada por:
Los antiguos griegos

aprovechaban ya la dE dm
energı́a del agua, P = = g · Hn · (3.1)
utilizaban ruedas
hidráulicas para moler
dt dt
trigo, sin embargo, la
posibilidad de emplear
esclavos y animales
retrasó su aplicación
Esta expresión puede ser llevada a una forma puramente algebraica
generalizada.
teniendo en cuenta la densidad del agua ρ y la eficiencia de la turbina y el
generador η como se muestra en la siguiente ecuación:
P = g · ρ · η · Q · Hn (3.2)
Q es el caudal por la tuberı́a de presión y g representa la aceleración de

la gravedad.
No siempre las estructuras monumentales aseguran un capacidad de

generación importante ya que ésta depende tanto del caudal como la
altura.
En muchos embalses la altura Hn permanece casi constante por lo

cual la potencia generada toma una forma lineal dada por un factor de
35
turbinamiento ξ
ξ = g · ρ · η · Hn (3.3)
En otros casos se opta por un polinomio (normalmente de grado 4)

para representar la relación no lineal entre
Entre las ventajas de las centrales hidroeléctricas se destacan:
• No requiere combustible por lo cual los costos directos de

explotación pueden considerarse como cero.
• Es una forma de generación limpia.
• Las presas necesarias para el procesos de generación, pueden

producir otros beneficios como la protección frente a inundaciones,
riego de cultivos, navegación y turismo.
• Larga vida útil.
• El sistema utilizado es eficiente y seguro.
Entre las posibles desventajas se destacas:
• Los costos de inversión son extremadamente altos.
• El lugar de emplazamiento suele estar alejado de los grandes centros

de consumo, lo cual hace necesario la construcción de extensas lineas
de transmisión.
• Por lo general la construcción tarda un considerable intervalo de

tiempo.
• La disponibilidad de la energı́a puede fluctuar según la estación del

año.
• El carácter estocástico de la hidrologı́a produce un riesgo de

indisponibilidad que no está presente en las plantas térmicas.
• Aunque el proceso de generación es limpio se produce gran impacto

ambiental por la construcción de las represas.
3.1 Clasificación de las centrales

hidroeléctricas
La primera central
hidroelectrica se contruyo
en 1880 en Cragside
En Colombia se presentan básicamente dos tipos de centrales
Northumberland (Gran
Bretaña)
hidraulicas: las centrales de pasada y las centrales de embalse.
Las primeras corresponden las centrales a filo de agua, normalmente

pequeñas centrales hidroeléctricas ubicadas en los sistemas de
sub-transmisión y en las zonas no interconectadas. Estas centrales
presentan un menor impacto ambiental aunque su capacidad de
generación es limitada.
Las centrales de embalse corresponden a las grandes represas las cuales

constituyen aproximadamente el 65% de la capacidad instalada del paı́s.
3.1.10 Central En este tipo de central no existe una acumulación apreciable de agua
de pasada (a filo por lo cual la capacidad de la planta depende del caudal disponible en
de agua) el rı́o. Como no es posible almacenar agua, la generación debe seguir
las fluctuaciones del caudal disponible, por tal razón este tipo de central,
requiere un caudal suficientemente constante para asegurar a lo largo del
año una potencia determinada.
De acuerdo con la potencia instalada, la Organización Latinoamericana

de Energı́a (OLADE) ha clasificado este tipo de centrales de acuerdo a la
potencia y caı́da en metros de la siguiente forma:
37
En Colombia, las centrales de pasada o PCH1 en derivación presentan

una configuración similar a la que muestra en la figura 3.1
Nivel del agua
Azud
Cámara de carga
Canal de baja presión
Bocatoma Bloque de apoyo
Tu
Hc
be
ría
Bloque de anclaje
de
pr
es
ió
n
Casa de máquinas
Figura 3.1: Esquema de una pequeña central hidroeléctrica
Tabla 3.1: Clasificación de acuerdo a la potencia

Potencia [kW] Tipo
0 - 50 Microcentral
50 - 500 Minicentral
500 - 5000 Pequeña central
Tabla 3.2: Clasificación de acuerdo a la caı́da en metros

Baja Media Alta >
Microcentral 0 - 15 15 - 50 50
Minicentral 0 - 20 20 - 100 100
PCH 0 - 25 25 - 130 130
A continuación se describen cada una de las partes constructivas de

una PCH.
1
Pequeñas centrales hidroeléctricas
Toma
Convencionalmente se construye una pequeña presa o azud2 , como se

muestra en la figura 3.2, la cual produce una elevación del nivel de agua
en el rı́o para permitir la captación del caudal requerido.
Canal de baja presió
Vertedero
Rejilla Compuerta de lavado
Canal de desfogue
o
a mpe a d Dentellón
Z
Azud
Figura 3.2: Toma convensional
Por lo general este tipo de presa esta diseñada para que el agua vierta
por encima de ellas mediante aliviaderos de coronación (vertederos).
El azud debe estar acompañado de de un zampeado el cual disminuye

la velocidad del agua a fin de disminuir la erosión; igualmente, se
construye un dentellon para disminuir la presión del agua y mejorar el
anclaje de la estructura.
La toma se regula por medio de compuertas ubicadas aguas arriba

2
Este tipo de presa es también conocida como presa de vertedero
39
(cierre de seguridad) y aguas abajo (cierre de regulación) del azud.
Una rejilla debe ser ubicada a la entrada de la toma para evitar la

entrada de elementos que arrastra el rı́o; dichas rejillas deben diseñarse
con una pequeña pendiente para permitir su limpieza periódica.
Canal de baja presión.
El canal de baja presión es utilizado para conducir el agua desde la toma

hasta la cámara de carga; Generalmente, la conducción es realizada por
medio de un canal abierto revestido interiormente para evitar filtraciones
en el terreno. La velocidad a los largo del canal se mantiene constante
y depende de los parámetros de diseño del canal como se muestra en la
ecuación de Chezy.3
γ 12 √
u= · R·S (3.4)
k
donde:
γ: el peso especifico del agua.
k: coeficiente de fricción.
R: radio hidráulico dado por el cociente entre el área y el perı́metro en

contacto con el agua.
S: pendiente longitudinal del canal.
Con el fin de evitar fenómenos de sedimentación o erosión la velocidad

en el canal debe estar entre 0.7m/s y 2m/s.
Una relación mas utilizada es la de Manning:

3
ver anexo ??
Q 1 2 √
= · R3 · S (3.5)
A n
en donde n es el coeficiente de rugosidad de Manning; algunos valores

promedio de este coeficiente se muestran en la tabla 3.3
Tabla 3.3: Valores promedio de n

Material n
Concreto terminado 0,012
Concreto rugoso 0,014
Ladrillo 0,016
Tierra 0,025
Grava 0,029
Desarenador
La función del desarenador es sedimentar las partı́culas de material solido

suspendidas en el agua. El no contar con un desarenador puede traer las
siguientes consecuencias:
• Disminución de la capacidad real del tanque de presión.
• Desgaste en la tuberı́a de presión por efecto de las partı́culas solidas

(Disminución de la vida útil).
• Disminución de la sección del canal por efecto de la sedimentación.
Sistema de presión
Cuando el salto es superior a 15m es necesario construir una cámara de

carga conectada a la tuberı́a forzada o tuberı́a de presión. La cámara
de carga o tanque de presión debe tener la capacidad suficiente para la
partida o parada brusca de las turbinas.
41
Un fenómeno a tener en cuenta en el diseño del sistema de presión es

el golpe de ariete. Este consiste en una variación de presión en la tuberı́a
forzada por encima o por debajo de la presión normal. Se produce por
fluctuaciones en el caudal, por ejemplo en el cierre de las válvulas que
regulan el caudal de entrada a la turbina. En este caso la velocidad en la
parte inferior de la tuberı́a disminuye hasta llegar al reposo, esto ocasiona
un aumento en la presión del liquido que se propaga como una onda en
dirección de la cámara de carga (golpe de ariete positivo).
Un fenómeno opuesto ocurre cuando el sistema requiere aumentar

la potencia generada; en este caso, el regulador aumenta la cantidad de
caudal en la turbina, ocasionando una depresión brusca en la tuberı́a
(golpe de ariete negativo).
En ambos casos, la elasticidad de las paredes de la tuberı́a debe ser

llevada en cuenta al igual que la compresibilidad del agua. Este fenómeno
de elasticidad es el causante de la oscilación en el sistema, la fricción
reduce esta oscilación en un movimiento amortiguado.
Una forma de reducir los efectos del golpe de ariete es mediante un

pozo vertical denominado almenara o chimenea de equilibrio el cual se
ubica entre la cámara de presión y la casa de máquinas; este elemento,
expuesto a la presión atmosférica por la parte superior, permite equilibrar
las presiones mediante variaciones en el nivel del liquido. Se requiere una
almenara cuando se cumple la siguiente relación:
L
≤5 (3.6)
HB
donde:
L : longitud de la tuberı́a forzada.
HB : Caı́da bruta.
En general, cuando el tiempo de cierre del dispositivo de control es

inferior a 3s se requiere la instalación de una chimenea de equilibrio.
3.1.11 Central Las centrales de embalse acumulan una considerable cantidad del
de embalse lı́quido ”aguas arriba” en una zona denominada reservorio. Los embalses
permiten regular la cantidad de agua que es turbinada por lo cual el
volumen del embalse depende de la potencia a generar y de las afluencias.
La capacidad de almacenar agua permite producir energı́a aun en épocas
de sequı́a o estiaje. Este tipo de centrales exigen una mayor inversión
de capital aunque los kilovatios-hora generados tienen un costo menor.
A pesar de ser una tecnologı́a limpia, puede ocasionar importantes
impactos ambientales y sociales que deben ser estudiados en la etapa de
planeamiento para establecer una adecuada relación costo beneficio.
En general los diferentes tipos de embalses presentan un volumen útil

y un volumen muerto.
Altura del golpe de la ola
Lamina de agua
Vertedero
Volumen útil
Toma
Altura por sedimentos

(Volumen muerto)
Figura 3.3: Esquema de una presa de gravedad

43
El volumen útil se refiere a la cantidad de agua que puede ser utilizada

para el proceso de generación de energı́a, este se calcula de acuerdo al
estudio hidrológico.
El volumen muerto es el espacio destinado a la acumulación de

sedimentos, este se calcula de acuerdo a la vida útil planeada de la
siguiente forma:
Vm = Vs · Tvu (3.7)
en donde:
Vm : Volumen muerto.
Vs : Volumen de sedimentos medios anuales del rı́o (m3 /año)
Tvu : Vida útil esperada en años.
Normalmente el volumen muerto es de aproximadamente el 10% del

volumen útil.
Presas
En el caso de una PCH la presa (denominada azud) tienen como función

única elevar el nivel del agua para permitir la toma. En los grandes
emplazamientos, las presas tienen en general una doble finalidad: obtener
una elevación del nivel del agua para obtener un salto mayor y crear
un deposito para almacenar agua y de esta forma regular la utilización
de la misma. Para esto se requiere una fuerte cimentación denominada
fundamento o apoyos laterales denominados estribos. La parte mas alta
de la presa (corona) puede ser utilizada para construir vı́as de acceso a la
presa. Las superficies frontales de la presa son denominados paramentos
y son construidos normalmente de concreto o arcilla. La clasificación
depende principalmente del tipo de material y la forma del paramento.
3.1.12
Los embalses pueden ser clasificados según el tipo de material con el
Clasificación
cual es construido:
• Presas de hormigón.
– De gravedad.
– De boveda.
– De contrafuerte.
• Presas de elementos sin trabar (escollera).
– De piedra.
– De tierra.
• Presas de mixtas.
Las presas de gravedad son construidas en hormigón presentando

una forma triangular. Son las más duraderas y requieren menor
mantenimiento. Su estabilidad radica en su propio peso, por tal razón
la altura de este tipo de presas está limitada por la resistencia del terreno
(figura 3.4).
Al igual que la anterior, las presas de bóveda o arco son construidas

en hormigón pero a diferencia de las primeras, este tipo de presas utilizan
menor cantidad de hormigón ya que su estabilidad no radica en su propio
peso sino en la posibilidad de distribuir la carga a las paredes de la
montaña mediante una curvatura con concavidad hacia el embalse. Por
esta caracterı́stica, las presas de arco deben ser construidas en cañones
estrechos y rocosos que puedan soportar la presión del agua; esto dificulta
localizar emplazamientos adecuados.
Las presas de contrafuertes (o aligeradas) están constituidas por una

pared que soporta el agua y una serie de contrafuertes o pilares que
soportan la presión de la misma transfiriéndola a las fundaciones. Al
45
Figura 3.4: Ejemplo de una presa de gravedad
utilizar elementos estructurales mas complejos (vigas, columnas, etc)

utiliza menos hormigón.
Finalmente, las presas de elementos sin trabar pueden ser de piedra,

cuando mas del 50% del material de la estructura es piedra o de
arcilla(tierra) en caso de usar materiales con granulometria pequeña. El
ancho de estas estructuras es de 4 a 7 veces mayor a su altura pues debe
evitar las filtraciones de agua las cuales son inversamente proporcionales a
la distancia horizontal. La presa puede se homogenea en caso de utilizar el
mismo material o heterogenea en caso de tener una mezcla de materiales
compactados con granulometrias variadas. Las presas de elementos sin
trabar no soportan crecidas por lo cual los estudios del caudal de crecida
son especialmente importantes en el planeamiento de la construcción.
Figura 3.5: Presa de arco : Represa Hover (USA)
3.2 Aspectos técnicos de las centrales de

embalse
Utilizando las curvas del nivel del plano topográfico es posible obtener
una relación entre la altura de la presa y el volumen de agua almacenado,
esta relación puede ser representada como un polinomio de cuarto orden
como se muestra en la ecuación 3.8.
V = a0 + a1 · H + a2 · H 2 + a3 · H 3 + a4 · H 4 (3.8)
47
Figura 3.6: Presa de tierra : Subestación Calima
Esta relación puede ser llevada a una forma lineal dado el orden de
magnitud de los coeficientes superiores, no obstante, la linealización de
la función supone una aproximación importante que puede afectar la
exactitud de los cálculos en algunas aplicaciones.
La casa de máquinas de una central generadora está compuesta por

un conjunto de unidades generadoras con caracterı́sticas semejantes. Las
principales caracterı́sticas son:
• Potencia efectiva: Es la máxima potencia activa que puede ser

generada en régimen permanente cuando la unidad entra en
operación.
• Turbinamiento efectivo: Es definido como la cantidad de agua

turbinada que sometida a la caı́da efectiva genera el valor de la
potencia efectiva.
• Altura de caı́da efectiva: es definida como la menor caı́da liquida

sobre la cual la unidad puede generar su potencia efectiva.
• Factor de utilización: se define como la relación entre la potencia

promedio histórica y la capacidad especificada durante el diseño.
Algunas plantas presentan un bajo factor de utilización, esto
significa que a pesar de tener un gran potencial de generación en
la práctica genera menor potencia debido a las afluencias limitadas,
es por ello importante realizar un estudio hidrológico antes de la
construcción de la central.
En la fase de diseño de la presa son definidos los valores de potencia

nominal, altura y caı́da nominal de cada central. Una vez terminada la
fase de construcción de la central y puesta en operación se determinan los
demás valores. Un concepto altamente utilizado en el análisis de sistemas
hidroeléctricos es la submotorización de la casa de maquinas, la cual está
relacionada con una baja capacidad de turbinamiento comparado con las
afluencias históricas del reservorio. Las usinas submotorizadas pueden
causar un efecto denominado estrangulamiento de la cascada en donde
una baja capacidad de turbinamiento obliga al vertimiento de excedentes
del agua que llega a la turbina.
Para cada altura la central es capaz de producir una potencia a costa de

un turbinamiento maximo (turbinamiento efectivo).
La producción de potencia activa presenta dos comportamientos

distintos en donde el limite entre estos dos comportamientos es la caı́da
liquida efectiva. La operación para alturas inferiores a la caı́da liquida
efectiva esta limitada por la turbina mientras que para alturas superiores
la limitación la impone la maquina sı́ncrona. Una forma simplicada de
determinar el turbinamiento efectivo máximo es mediante la ecuación:
α
HL
Qen = Qef · (3.9)
Hef
en donde alfa toma diferentes valores según el caso:
• α = 0.5 cuando HL ≤ Hef y se trata de una turbina pelton o francis.
• α = 0.2 cuando HL ≤ Hef y se trata de una turbina kaplan.

49
• α = −1 cuando HL ≥ Hef
La potencia máxima se puede calcular como:
β
HL
P = Pef · (3.10)
Hef
Con β igual a 1.5, 1.2 y 0 respectivamente para los tres casos

considerados.
En cuanto a la eficiencia se tiene que para las grandes centrales

hidroelectricas el rendimiento está entre 0, 88 y 0, 94.
1. Ortiz Flórez, Ramiro.Pequeñas centrales hidroeléctricas. Editorial

McGrawHill. Bogotá 2001.
2. http://thales.cica.es/rd/Recursos/rd99/ed99-0226-01/paginaprincipal.html
3. G. Zoppeti. Centrales hidroeléctricas. Editorial Gustavo Gili.
4. Eletrobrás. Diretrizes para projetos de PCH. Manual de eletrobras

disponible en internet en www.eletrobras.com.br
Hidrologı́a básica
4
Contenido
Reservorios
Caudal asegurado
Caudal de crecida
51
52 CAPÍTULO 4. HIDROLOGÍA BÁSICA
La hidrologı́a es una ciencia clave en el estudio de los sistemas de

generación. Ésta permite determinar la capacidad real de la central y
planear la operación de la misma teniendo en cuenta la estacionalidad de
las afluencias. En este capı́tulo se estudiaran los principales conceptos en
la materia, aplicados a las centrales de generación hidráulica.
4.1 Reservorios
Conforme a la capacidad de regulación, las plantas son clasificadas

como de acumulación y de compensación. Los reservorios de
acumulación tienen gran capacidad de almacenar energı́a en forma de
agua. Por tal razón son responsables de la regulación del caudal del rı́o en
cada ciclo hidrológico. Las centrales con reservorios de acumulación son
denominadoas de reservorio. Los reservorios de compensación presentan
poca capacidad de regulación por lo que las plantas de compensación son
denominadas a filo de agua.
En una central a filo de agua el volumen máximo y mı́nimo son iguales

por lo que se espera un volumen útil constante.
4.1.13 Algunas Antes de iniciar el estudio detallado del caudal como una variable
definciones estocástica importante en la operación de los sistemas hı́dricos, es
necesario definir algunos conceptos:
Caudal de crecida
Corresponde al caudal máximo registrado (o esperado) en un

aprovechamiento; este valor es de especial interés para el diseño de
las obras civiles (azud, canal de conducción, vertederos, etc) pues éstas
deben estar en capacidad estructural de soportar estos caudales de
crecida.
53
La caracterı́stica estacional de las afluencias permite determinar un

Periodo de retorno asociado con el cual se diseñan las obras de desvió y
captación. Según la OLADE1 los periodos de retorno recomendado para
la captación son los siguientes:
• Microcentrales: 20 a 25 años
• Minicentrales: 50 a 100 años
• Pequeñas centrales: 100 a 150 años
Caudal mı́nimo
Debe ser definido a partir de los estudios ambientales por lo que es

denominado caudal ecológico. Este tiene especial importancia en los
sistemas que requieren obras de desviación o retención como túneles,
canales de baja presión o embalses los cuales limitan el caudal disponible
aguas abajo del emplazamiento.
Caudal asegurado
La estocasticidad en las afluencias lleva con sigo una variabilidad de la

potencia disponible en la unidad hidráulica, por tal razón, es importante
definir el caudal mı́nimo que se puede asegurar con una probabilidad alta
(normalmente entre el 85% al 95%).
Mientras que el caudal de crecida determina el diseño de las

estructuras civiles, el caudal asegurado permite el diseño de los elementos
electromecanicos y el cálculo de la potencia disponible final.
En la gran mayorı́a de pequeñas centrales hidroeléctricas, el caudal

asegurado es el que finalmente será conducido por la tuberı́a de presión
hacia la casa de máquinas, por lo que la potencia eléctrica resultante será
constante.
1
Organización Latinoamericana de Energı́a
En otros tipos de centrales, la potencia de salida puede ser tan variable

como las afluencias por lo que el caudal asegurado permite determinar el
factor de utilización de la misma. El caudal asegurado se calcula mediante
un estudio de variables promedio como se muestra en la siguiente sección
4.2 Estudio de variables promedio
Como se mostró en la ecuación 3.2 la potencia disponible en un

aprovechamiento depende principalmente de la altura neta y el caudal
disponible; no siempre se cuenta con mediciones históricas de caudal
por lo cual es necesario realizarlas manualmente, a este procedimiento se
denomina aforo2 .
La naturaleza estocástica de las afluencias hace que la potencia

disponible sea igualmente aleatoria. Para analizar esta caracterı́stica se
utiliza una gráfica de caudal vs tiempo denominada hidrograma. A partir
de los datos consignados en el hidrograma es posible construir una curva
de permanencia o duración con el cual se puede calcular la potencia
asegurada.
Al igual que la curva de duración de carga en el caso de la demanda

de energı́a, la curva de permanencia determina la cantidad de tiempo en
el cual se puede asegurar un caudal determinado o superior; si este de
tiempo se divide por el periodo total de las mediciones entonces se tiene
una probabilidad asociada.
A partir de los datos de caudal del año de mayor estiaje3 se construye

una tabla de frecuencia acumulada ordenando los intervalos de clase de
mayor a menor. El número de marcas de clase NC puede ser determinado
a partir del número total de datos N mediante la siguiente relación:
2
En el anexo A se presentan los principales métodos de aforo para para proyectos de
pequeña escala.
3
sequia
55
NC = 1 + 3, 3 · ln (N ) (4.1)
El ancho del intervalo de clase se calcula en base a los caudales

extremos (Qmax y Qmin ) y al número de intervalos de clase:
Qmax − Qmin
ΔQ = (4.2)
Nc − 1
Se determina la frecuencia de cada intervalo como el número de veces

que que cada dato está en el intervalo divido por el número total de
datos. Finalmente la curva puede ser aproximada a una forma continua
por ejemplo a una función exponencial.
En la figura 4.2 el punto de 95% indica que se asegura un caudal igual

o superior a Q con una probabilidad del 95%.
95 %
Figura 4.1: Curva de permanencia y duración
Este valor tiene especial interés para el estudio de factibilidad de una

central pues indica el caudal asegurado de la misma. Es posible ajustar los
puntos de la curva a una forma continua mediante un método de mı́nimos
cuadrados.
No siempre es posible tener datos de caudales en el punto de

aprovechamiento por lo cual es necesario hacer utilizar la ley de
continuidad que define la siguiente relación:
Aa
Q= · Qm (4.3)
Am
en donde
Aa : Área en el punto de aprovechamiento.
Am : Área en el punto de medición.
Qm : Caudal en el punto de medición.
En algunas ocasiones no se cuenta con datos confiables de caudales en

el punto de emplazamiento, esta situación se presenta cuando la relación
entre las áreas del punto de medición y aprovechamiento es muy alta (≥ 4)
o a que los datos son muy pocos en el tiempo ≤ 1 año.
Ejemplo 4.1
Construir la curva de permanencia para los siguientes datos de caudales:
El número de intervalos de clase se calcula como:
Nc = 1 + 3.3 · ln(12) (4.4)
se tienen 9 intervalos de clase. El ancho de cada intervalo de clase es

0.2625
La frecuencia relativa se calcula como la frecuencia por cada intervalo

sobre el número total de muestras N . La tabla 4.2 muestra el
procedimiento:
57
Tabla 4.1: Datos de caudales

Mes Caudal Q[m3 /s]
Enero 4.15
Febrero 3.58
Marzo 3.87
Abril 3.45
Mayo 3.21
Junio 2.05
Julio 2.35
Agosto 2.38
Septiembre 3.51
Octubre 3.84
Noviembre 3.59
Diciembre 3.97
En la gráfica 4.2 se muestra el resultado obtenido usando una

aproximación exponencial.
Q[m3 /s] = 4.8128 · e−0.0069·P [%] (4.5)
Con esta aproximación exponencial es posible determinar el caudal

asegurado con una probabilidad del 95% el cual es la base para calcular
la potencia disponible. En este caso particulas el caudal asegurado es de
2, 5008m3 /s
La curva de duración de caudales permite además determinar algunas

variables importantes en la central como la potencia asegurada la cual
puede ser calculada a partir de la ecuación 3.2. En el caso de una central
a filo de agua en derivación, ésta es la potencia nominal de la central,
mientras que en una central de embalse, es posible diseñar la turbina y
el generador para una potencia nominal mayor, teniendo en cuenta que
ésta potencia no se puede alcanzar el 100% del tiempo de operación. La
potencia obtenida por medio del caudal asegurado constituye entonces
una variable de suma importancia para el análisis economico de un
proyecto de generación.
Tabla 4.2: Determinación de la frecuencia acumulada.
Rango de clase Frecuencia Frecuencia Frecuencia

Marca de
Xi Xf f relativa fR acumulada fR clase Xm
4,15 3,89 2 16,67 16,67 4,02
3,89 3,63 2 16,67 33,33 3,76
3,63 3,36 4 33,33 66,67 3,49
3,36 3,10 1 8,33 75,00 3,23
3,10 2,84 0 0,00 75,00 2,97
2,84 2,58 0 0,00 75,00 2,71
2,58 2,31 2 16,67 91,67 2,44
2,31 2,05 1 8,33 100,00 2,18
4.3 Estudio de variables extremas
Existe una clara diferencia entre las variables medias y las variables
extremas pues mientras que las primeras tienden a una distribución
normal (debido al teorema del lı́mite central) las variables extremas
pueden tener distribuciones diversas. El el caso de las afluencias en
una central hidroeléctrica, existen menor cantidad de datos históricos
asociados a caudales de crecida por lo que corre mayor riesgo en el
análisis.
La probabilidad que una variable aleatoria (caudal) tome un valor

igual o inferior X está dada por una función de distribución de
probabilidad:
P (x ≤ X) = F (x) (4.6)
Si se tienen n observaciones independientes, la probabilidad Φ(x) de

que todos los valores resulten inferiores a X disminuye a una distribución
dada por el producto de los sucesos.
59
4.5
3.5
2.5
20 30 40 50 60 70 80 90 100
Figura 4.2: Ejemplo curva de masas

n
Φ(x) = P (xk ≤ X) = (F (x))n (4.7)
k=1
Ası́ que la probabilidad de que el mayor valor de x sea mayor a X

después de n observaciones es:
Φ(x) = 1 − (1 − F (x))n (4.8)
De lo anterior se tiene las siguientes observaciones:
• Se puede conocer la distribución Φ(x) a partir de F (x)
• Ésta función depende del número de observaciones.
La probabilidad p de ocurrencia de de un fenómeno puede ser asociado

con una frecuencia por lo cual el periodo de retorno T corresponde al
inverso multiplicativo:
1
p = P (x ≥ X) = (4.9)
T
En el caso de los caudales de crecida, el periodo de retorno T indica el

tiempo en el que se espera un caudal igual o superior al caudal de crecida
(Qmax ) con un riesgo asociado p. Las estructuras de contención como la
presa deben ser diseñadas para soportar este caudal por lo cual el estudio
de las variables extremas es clave para la seguridad del embalse.
Los valores extremos son independientes entre sı́, es necesario entonces

asignar un valor de probabilidad empı́rica a cada dato. Para ello se ordenan
los n registros históricos de forma decreciente sin importar la fecha de
ocurrencia, a cada uno de estos registros se le asigna una probabilidad en
función del orden i en la lista ordenada. Existen varias probabilidades
empı́rica históricamente utilizadas como se muestra en la tabla 4.3:
Tabla 4.3: Probabilidades empı́ricas

Nombre Expresión
California (1923) i/n
Hazen (1930) (2 · i − 1)/(2 · n)
Weibull(1939) i/(n + 1)
Gringorten(1963) (i − 0, 44)/(n + 0, 12)
Por ejemplo la probabilidad tipo Weibull asigna un periodo de retorno

n + 1 años al mayor de los datos históricos y n+1
n
al menor. Los métodos
más recomendables son los tres últimos aunque no se ha dicho la última
palabra en este tema en la literatura internacional.
Una vez asignada la probabilidad empı́rica, es necesario hacer un

ajuste de los datos a una distribución de probabilidad, para ello se pueden
usar pruebas como la de Kolmogorof-Smirnof o la de ψ 2 (ver anexo ??).
No obstante una forma simple y relativamente precisa es la regresión por
mı́nimos cuadrados.
Los métodos de mı́nimos cuadrados permiten ajustar un conjunto de

puntos a una lı́nea recta de la forma:
61
y =a+b·x (4.10)
En el caso de las distribuciones de probabilidad es posible realizar un

ajuste a una curva de la curvas dadas en la table 4.4
Tabla 4.4: Distribuciones de probabilidad ajustables

Nombre Expresión
−y
Gumbel Φ(y) = e−e con y = α · (x − μ)
Exponencial Φ(y) = e−y
y
Fréchet Φ(y) = e−e con x = α · (ln(x) − μ)
Para utilizar el método convensional de mı́nimos cuadrados es

necesario llevar la expresión a una lı́nea recta; en el caso de la función
de Gumbel se utiliza un doble logaritmo de la función Φ(x). Los valores
Φi = Φ(xi ) se asignan de acuerdo a la probabilidad empı́rica seleccionada
de tal forma que pi = 1 − Φi .
El ajuste encontrado permite realizar una extrapolación para

determinar el caudal de crecida para periodos de retorno más altos, esto
minimiza el riesgo de que un crecida presente un caudal superior al valor
diseñado para las diferentes estructuras.
La metodologı́a planteada puede ser utilizada para el cálculo de dos

caudales de crecida: el requerido para las obras de desvı́o y el de las
estructuras civiles propiamente.
A partir de la ecuación 4.8 es posible calcular el riesgo asociado Φ para

un periodo de retorno T y una vida útil n. Los valores del tiempo de
recurrencia son normalmente del orden 100 años para pequeñas centrales
hidroeléctricas con azud. Para centrales de mayor capacidad el tiempo de
recuerrencia se calcula de acuerdo al tipo de presa:
• Presas en concreto: 500 años.
• Presas en tierra: 1000 años.

Ejemplo 4.2
Realizar el análisis de valores extremos para el conjunto de datos presentado

en la table 4.5.
Tabla 4.5: Valores de caudal de crecida

Fecha Caudal de crecida
2005 15,2
2004 17,3
2003 19,5
2002 8,0
2001 9,3
2000 11,2
1999 13,0
1998 14,0
1997 16,8
1996 17,0
1995 15,7
1994 17,3
Inicialmente es necesario ordenar los términos de mayor a menor y

asignarles una probabilidad empı́rica como se muestra en la tabla 4.6 (En
este caso se usará la ecuación de Weibull).
Los valores empiricos se pueden ajustar a una curva dada por la

distribución de Gumbel ası́:

1
y = −ln ln (4.11)
1 − pi
El ajuste lineal da como resultado α = 0, 2706 y μ = 12, 6641. En la

gráfica 4 se muestra el ajuste obtenido.
No obstante como se muestra en la figura, es posible obtener dos

ajustes diferentes: el de los caudales como función de las probabilidades
63
Tabla 4.6: Probabilidades asignadas a los datos de la tabla 4.5

Caudal (x) i Probabilidad p
19,5 1 0,0769
17,3 2 0,1538
17,3 3 0,2308
17,0 4 0,3077
16,8 5 0,3846
15,7 6 0,4615
15,2 7 0,5385
14,0 8 0,6154
13,0 9 0,6923
11,2 10 0,7692
9,3 11 0,8462
8,0 12 0,9231
y el de éstas en función de los primeros. Para obtener un resultado mas

ajustado a la realidad es posible usar un ajuste ortogonal (ver anexo ??) el
cual minimiza el error entre la curva teórica y los datos reales teniendo
en cuenta un error presente en ambas variables, en otras palabras, se
minimiza la distancia medida perpendicularmente entre cada punto y la
recta de ajuste.
Finalmente, es posible obtener el caudal de crecida mediante el tiempo

de retorno el cual se puede calcular según el riesgo a que se está dispuesto.
Para las obras de desvı́o se puede tomar un periodo de retorno T = 150
años con una duración de 1 año en la construcción. En este caso el riesgo
de que el caudal sea mayor al establecido es de 0, 67%.
Para las obras definitivas se puede considerar un periodo de retorno

de T = 500 años y una vida útil de n = 50 años por lo que el riesgo es:
n
1
Φ=1− 1− = 9, 5% (4.12)
T
El caudal de crecida para esta situación se puede encontrar a partir

de las ecuaciones 4.11 y 4.10. En este caso es de Q = 35, 63. Como se
0.9
Ajustedelaprobabilidad
alcaudal
0.8
Probabilidaddeocurrencia
0.7
0.6
Ajustedelcaudal
alaprobabilidad
0.5
0.4
0.3
0.2
0.1
0
8 10 12 14 16 18 20
Caudaldecrecida
Figura 4.3: Gráfica de probabilidad vs Caudal de crecida
observa es un valor mucho mayor a los datos históricos lo que significa

una aproximación muy conservadora.
Maggio, Guillermo Eduardo . Análisis estadı́stico de variables extremas

Buenos Aires. 2003. Aplicaciones en hidrologı́a. (www.noldor.com.ar).
Turbinas hidráulicas
5
Contenido
Caracteristicas generales
Turbinas Pélton
Turbinas Francis
Turbinas Kaplan
65
66 CAPÍTULO 5. TURBINAS HIDRÁULICAS
L ATurbina es un dispositivo capaz de transformar la energı́a hidráulica

en energı́a mecánica. Su importancia en un sistema de generación de
energı́a radica en su relación con la estabilidad del sistema y el control
de la potencia activa, ası́ mismo, los costos asociados a cada tipo de
tecnologı́a tienen un impacto importante en el análisis económico de la
planta. Aunque la turbina es un dispositivo puramente mecánico, ésta
tiene un fuerte impacto en variables eléctricas que tienen que ver con la
eficiencia y la estabilidad de la central. Ası́ mismo, condiciones como las
salidas por mantenimiento de éste elemento pueden afectar la operación
económica del sistema y la confiabilidad del mismo.
En el estudio de las turbinas hidráulicas no se tienen en cuenta los

cambios de temperatura y se analizan en estado estable. Las turbinas
hidráulicas más utilizadas son las tipo Pelton, Francis y Kaplan. Cada una
de éstas presenta comportamientos distintos que deben ser analizados.
En este capı́tulo se presentan las principales caracterı́sticas estos tipos de
turbinas y algunos criterios básicos para la selección de las mismas.
5.1 Caracterı́sticas generales
Las máquinas hidráulicas al igual que las máquinas eléctricas

presentan una parte fija (distribuidor) y una parte móvil (rodete). Los
elementos constitutivos principales son los siguientes:
• distribuidor : su misión es dirigir el agua hacia el rodete. Convierte

la energı́a de presión en energı́a cinetica.
• rodete: está provisto de álabes móviles que permiten convertir

la energı́a cinética y de presión del agua, en energı́a cinética
rotacionales en el eje de la turbina para transferirla al generador
sı́ncrono.
• difusor: también denominado tubo de aspiración, permite el escape

del agua disminuyendo la presión hasta el vacı́o. Si la turbina no
posee tubo de aspiración, se denomina de escape libre.
67
5.2 Tipos de turbinas hidráulicas
De acuerdo a la presión del agua en el rodete se presenta la siguiente

clasificación:
• Turbina de acción (pelton): en este tipo de turbina el agua sale

del distribuidor a presión atmosférica por lo cual la totalidad de
la energı́a potencial gravitacional del agua se transforma en energı́a
cinética.
• Turbina de reacción (francis, kaplan): la energı́a potencial del agua

se transforma en parte en energı́a cinética y en parte en energı́a
de presión, por tal razón, el agua sale del distribuidor a una cierta
presión que disminuye a medida que el agua atraviesa el rodete
llegando a una presión de vacı́o en el difusor.
De igual forma, es posible presentar una clasificación de acuerdo con

la dirección del agua:
• Axiales (kaplan, hélice, bulbo): el agua ingresa en forma paralela al

eje de la turbina.
• Radiales (francis): el agua ingresa perpendicular al eje de la turbina.
• Tangenciales (pelton): el agua ingresa en forma lateral o tangencial

contra los álabes.
En general una turbina va acoplada a un alternador que ha de generar

electricidad a una determinada frecuencia (60Hz), por lo que su velocidad
debe ser tal que, conjugada con el número de pares de polos, produzca
dicha frecuencia. Para mantener la velocidad constante se cuenta con
un sistema de regulación de velocidad. En la medida que la carga en
el sistema aumenta la velocidad disminuye. En el caso de las turbinas
hidráulicas la velocidad es función del flujo de agua por lo cual el control

de frecuencia va ligado al caudal útil.
En ambos casos el mecanismo de control de frecuencia es llamado

gobernador. Éste toma medidas de la velocidad y hace ajustes en
la entrada de potencia mecánica para retornar a la velocidad sı́ncrona
después de una variación en la carga eléctrica del sistema.
5.2.14 Turbina Clasificada dentro de las turbinas de acción y de chorro libre

pelton tangencial, es la más utilizada para grandes saltos (Entre 60 y 1500 m). En
la figura 5.1 se muestra una turbina pélton de 6 inyectores y eje vertical.
Los álabes o casolétas reciben el chorro del inyector e impulsa al rodete

el cual puede ser de eje vertical u horizontal. Las turbinas pelton de eje
horizontal disponen de menos de tres inyectores pues de contar con más,
serı́a imposible evitar que el agua cayera sobre la rueda misma a la salida
de las casoletas.
La turbina Pelton debe su

En la turbina pelton al igual que en todas las turbinas de acción,
nombre al norteamericano
Lester Allan Pelton (1829-
la energı́a potencial gravitacional se transforma en energı́a cinética por
1908), quien buscando oro
en California, concibió la
medio del inyector(distribuidor). Utilizando el teorema de Bernoulli se
idea de una rueda con
cucharas periféricas que
obtiene la velocidad en la salida del distribuidor:
aprovechara la energı́a
cinética de un chorro de
agua proveniente de una
tuberı́a a presión (tobera),
incidiendo de forma
tangencial sobre la misma.
Ensayó diversas formas de
vc2 Pc v2 Pd
álabes hasta alcanzar una
Hc + + = Hd + d + (5.1)
2·g 2·g
primera patente de la
rueda en 1880 aunque γ γ
como tal, la turbina fue
patentada en 1889.
En la cámara de carga (c) la presión manométrica es cero al igual que

en la salida del distribuidor (d), igualmente, la velocidad en la entrada
de la tuberı́a forzada se asume despreciable, por tal razón, la ecuación de
Bernoulli toma la siguiente forma:
69
Figura 5.1: Esquema de una turbina pelton de eje vertical
vd2
HN = (5.2)
2·g
en donde HN corresponde al salto neto, luego de restar la pérdidas en

la tuberı́a de presión1 .
La velocidad vd se ve disminuida a la salida del distribuidor debido a

las pérdidas del mismo; la entrada del agua en la cazota (álabe) presenta
un valor de velocidad vea tangencial al rodete mientras que a la salida del
1
Las pérdidas en la tuberı́a de presión pueden aproximarse al 4% de HN para las PCHs
mismo presenta una velocidad vsa de igual magnitud a un ángulo α con

respecto a la dirección tangencial. La velocidad relativa del agua vr dada
por la diferencia entre la velocidad del agua v y la velocidad tangencial
del rodete u, representa el grado de compromiso entre la velocidad y el
torque generado. Un valor elevado de w puede hacer que el agua choque
de forma brusca con los álabes mientras que un valor muy bajo produce
que el agua se deslice suavemente por el mismo sin transmitirle cantidad
de movimiento al eje.
La componente horizontal de la fuerza sobre el álabe puede ser

determinada en función de la velocidad relativa del agua dado que ésta
es igual en magnitud a la entrada y la salida del álabe (ver figura 5.2):
RPM
Vs
a
u
Vr
V ea Vr
Figura 5.2: Esquema del álabe de una turbina pelton
F = ρ · Q · (vr − vr · cos(α)) (5.3)
donde
vr = vea − u (5.4)
La potencia que se transfiere al rodete es:
P = u · ρ · Q · (vr − vr · cos(α)) (5.5)

71
De esta ecuación se deduce que la máxima potencia se obtiene con un

ángulo α = 180o y una velocidad tangencial del rodete de u = v/2. En
la practica los álabes se construyen en forma de cuchara doble como se
muestra en la figura 5.3 para evitar empujes no balanceados en la flecha.
Figura 5.3: Esquema de los álabes de una turbina pelton
Para diseñar el diámetro D del rodete y el numero de polos de la

máquina sı́ncrona se deduce la siguiente expresión:
π
u= · D · RP M (5.6)
60
El caudal de salida en función del diámetro de la boquilla d se puede

obtener como:
π 2
Q= · d · vd (5.7)
4
1 d 1
En la practica 200
< D
< 7
con lo cual se puede establecer el diseño.
La eficiencia de la turbina pelton puede ser aproximada mediante la

expresión:
Cv2
η= · (1 + ζ) (5.8)
2
en done:
Cv : Coeficiente de la boquilla (representa la pérdida de velocidad del

agua en el inyector V = Cv · vd ).
ζ: representa la diferencia entre las velocidades relativas de entrada y
salida del agua.
Las turbinas Pelton presentan de 1 a 6 inyectores. La turbina de un

inyector ofrece ventajas en cuanto a economı́a y simplicidad además que
permite ubicarla en un eje horizontal; Las turbinas de varios inyectores se
construyen con eje vertical y permite mayor velocidad en la turbina.
5.2.15 Turbina Las turbinas francis presentan admisión radial y centrı́peta con tubo de
francis aspiración para evacuar el agua. Convensionalmente se construyen para
un rendimiento máximo y se clasifican en rápidas, normales y lentas; la
diferencia entre cada tipo radica en la forma y tamaño del rodete:
• Rodetes lentos: se usan para grandes saltos, su velocidad en RPM es

menor y su diametro es mayor con respecto al tubo de aspiración.
Los álabes aumentan su espesor para evitar remolinos.
• Rodetes normales: el agua ingresa radialmente y sale de forma axial.
La velocidad es mayor al tipo anterior y el diametro del rodete es
levemente mayor al del tubo de aspiración.
• Rodete rápido: en este tipo de rodete el cambio de dirección del agua
se efectúa mas bruscamente. Existe una distancia grande entre el
rodete y el distrbuidor, ası́ mismo el diametro es menor al del tubo
de aspiración.
Un elemento importante en las turbinas de reacción es la cámara

espiral, ésta tiene como función dirigir el agua en el distribuidor de
forma que la velocidad de entrada permanezca constante y con magnitud
relativamente baja.
73
Figura 5.4: Esquema de una turbina francis
Por otro lado el distribuidor tiene como función dirigir el agua a los
álabes del rodete, igualmente permite el control del caudal modificando
el angulo de entrada del agua. El regulador de velocidad actúa sobre las
directrices móviles. En general el tubo de aspiración ocasiona pérdidas de
energı́a hidraulica.
Finalmente el tubo de aspiración juega un papel importante en el

desempeño de las turbinas de reacción, en el caso de las turbinas francis de
rodete lento, el tubo de aspiración permite crear una presión de vacı́o en
la salida del rodete aumentado la eficiencia de la turbina. En las turbinas
francis de rodete rápido, al igual que en las kaplan y las turbinas hélice, la
función del tubo de aspiración es simplemente extraer el agua.
En la gráfica 5.6 se muestra una turbina francis y el diagrama de

presión asociado a ésta. Los puntos marcados con números son los
siguientes:
Figura 5.5: Cámara espiral en una turbina de reacción
1. Entrada del agua

2. Entrada al rodete
3. Salida del rodete
4. Salida del tubo de aspiración
Aplicando el teorema de Bernoulli en el punto 3 se tiene

P3 v2
H = Hs + Hef + + 3 + ht + hr + hd (5.9)
γ 2·g
en donde ht ,hr y hd son respectivamente las pérdidas en tuberı́a de

presión, en el rodete de la turbina y en el distribuidor de la misma. La
eficiencia puede ser entonces calculada como función de la energı́a que se
transfiere finalmente a la turbina :

Hef Hs ht + hd + hr P3 v32
η= =1− + − + (5.10)
H H H γ·H 2·g·H
75
1
1
Camara
Espiral
Presión
Rodete
Tube
ri
P re s a 2 2
ion
3
3
Tubo
Aspiración
4
4
Figura 5.6: Diagrama de presión de la turbina
donde, cuando P2 es negativa la eficiencia aumenta lo cual se logra

mediante el tubo de aspiración al impedir que la presión atmosférica
ingrese a la cámara.
5.2.16 Turbina
Las turbinas kaplan (figura 5.7) son utilizadas para saltos bajos y kaplan
caudales grandes. Es altamente efieiciente y controlable pero su costo es
superior al costo de las turbinas pélton y francis.
Los álabes móviles del rodete permiten una alta eficiencia aun en casos
de cargas parciales. Entre sus caracterı́sticas se destacan:
• Dimensiones relativamente bajas.
• Alta velocidad de rotación.

Figura 5.7: Esquema de una turbina kaplan
• Alta curva de eficiencia.
• Gran capacidad de sobrecarga.
• Buena capacidad frente a cargas parciales.
El rodete presenta pocos álabes orientados radialmente con la

capacidad de modificar su posición. La curvatura de estos álabes permiten
77
pocas pérdidas de carga.
La turbina kaplan fue

El rango de aplicación de las turbinas kaplan se ha incrementado. Ası́ desarrollada al rededor de
1915 por el ingeniero
mismo, algunos proyectos de generación hidroeléctrica descartados desde austriaco Viktor Kaplan
(1876-1934). La turbina es
el punto de vista económico o ambiental se hacen viables al usar este tipo similar al sistema de
propulsión de un barco.
de turbina.
5.3 Elección de la turbina en una central

hidráulica
Generalmente la elección de la turbina en un proyecto de generación

hidroeléctrico esta definido por la altura y el caudal de diseño. Las
turbinas mas utilizadas para grandes caı́das son las Pelton y las Francis,
a continuación se presenta un cuadro comparativo entre estas dos
tecnologı́as:
Tabla 5.1: Caracterı́sticas de las turbinas pelton y fracis

Pelton Francis
Robusta Menor peso
Menor erosión en los álabes Mayor rendimiento máximo
Reparaciones más sencillas Mayor aprovechamiento
Facilidad en la regulación Alternador más económico
Mejor rendimiento a cargas parciales Dimensión en planta reducida
Para alturas medianas la elección es normalmente entre turbinas

kaplan y francis, en la tabla 5.2 se muestran sus caracterı́sticas:
Tabla 5.2: Caracterı́sticas de las turbinas fracis y kaplan

Francis kaplan
Menor costo Mejor rendimiento a cargas parciales
Menor peligro de cavitación Mejor rendimiento a alturas variables
Una forma simple de seleccionar la turbina es mediante la gráfica 5.8

la cual relaciona caudal contra altura lı́quida. En general las turbinas
de acción como la pélton son utilizadas para caı́das grandes y caudales
relativamente bajos mientras que las turbinas de reacción como la kaplan
son usadas para caudales grandes y caı́das bajas. Las turbinas francis
tienen una amplia gama de aplicaciones por lo cual las curvas se cruzan
en algunos puntos.
Altura (m)
1000
500
300 Pelton
200
100 Francis
50
30
20
10
Banki Kaplan
5
3
2
3
Caudal (m /s )
0,2 0,5 1 2 3 5 10 20 30 50 100
Figura 5.8: Grafico de selección de turbinas
Otro elemento importante es el rendimiento. En una unidad

hidroeléctrica el proceso de transformación de energı́a cinética en energı́a
eléctrica esta sujeto a un rendimiento de las maquinas presentes en el
proceso. La función de rendimiento de una central en particular puede ser
representado como función de la caı́da liquida disponible de la potencia a
obtener. Con estos datos se puede crear una curva de rendimiento la cual
79
debido a su forma esta es llamada curva colina. Ası́ mismo, existen curvas
como la de la figura 5.9 la cual muestra el rendimiento de diferentes tipos
de turbinas frente a cargas parciales.
100
Kaplan
80
Pelton
Eficiencia (%)
60
FrancisRápida
40
FrancisLenta
20
Caudal (% Qn) 20 60 100 140
Figura 5.9: Eficiencia de diferentes turbinas frente a cargas parciales
Las curvas de rendimiento de los diferentes tipos de turbinas pueden

ser encontradas en los manuales entregados por los fabricantes. No
obstante, existen modelos generales desarrollados para las principales
turbinas hidráulicas.
Como se mostró anteriormente, el caudal utilizado en la central

puede ser dividido en varias turbinas, las cuales por facilidad de
manteniemiento se asumen de las mismas caracterı́sticas. En esta situación
es recomendable construir curvas de rendimiento para 1, 2 o mas turbinas
de forma independiente.
5.4 Control de velocidad (Gobernador)
El gobernador de la turbina es un equipo de control que permite

ajustar la salida de potencia mecánica para mantener un equilibrio entre
la potencia demandada y la potencia generada y de este forma mantener
la frecuencia del sistema en 60 Hz. Cuando decrece la demanda, el exceso
de potencia activa generada acelera la turbina haciendo que la frecuencia
aumente, para compensar este efecto, el gobernador de velocidad debe
reducir la admisión de caudal y de esta forma balancear el sistema.
Regulador
Aguja
S istema Inyector
de
Control
Álabe
Deflector
Rodete
Figura 5.10: Ejemplo del regulador de una turbina Pelton
En el caso de las turbinas Peltón el control es realizado por medio de

la aguja del inyector la cual regula la cantidad de caudal que choca con
81
los álabes. Igualmente, en este tipo de turbinas existe un dispositivo

adicional de control denominado deflector el cual modifica la dirección
del chorro disminuyendo el caudal efectivo sobre los álabes del rodete
sin modificar el caudal de salida del inyector; de esta forma se evitan los
efectos adversos producidos por el fenómeno del golpe de ariete sobre la
tuberı́a de presión.
En las turbinas francis el control es realizado por medio de los álabes

del estator.
1. RETScreen International. Small hydro project analysis.

www.retscreen.net.
2. Kjolle Arne. Mechanical Equipment. Norwegian University of

Science and Technology.
3. http://www.microhydropower.net/turbines.html
4. Streeter, Victor L and Wylie Benjamin. Mecánica de los fluidos.

McGrawHill.
5. Dı́ez, Pedro Fernández. Turbinas Hidraulicas. Departamente de

ingenierı́a eléctrica y energética. Universidad de Cantabria.
6. Niño, Jose Roberto; Duarte Carlos Arturo. Hidráulica de tuberı́as y

máquinas hidráulicas. Universidad Nacional de Colombia. Facultad
de Ingenierı́a. 2006.
Aspectos operativos.
6
Contenido
Curva de capacidad de la máquina
sincrona
Control de tensión
Control automático de la generación
83
84 CAPÍTULO 6. ASPECTOS OPERATIVOS.
L OSsistemas eléctricos deben ser operados de forma coordinada para

garantizar criterios de calidad y seguridad. En este capitulo se muestra
de forma general aspectos relacionados con el control de tensión y
de frecuencia realizao por cada una de las unidades de generación.
Inicialmente se define la curva de capacidad de las máquinas sı́ncronas la
cual es importante para entender conceptos como el costo de oportunidad
de la potencia reactiva. Posteriormente se presentan algunos conceptos
generales sobre el control de tensión en el sistema. Finalmente se muestran
algunas generalidades sobre el control de frecuencia.
6.1 Curva de capacidad de un generador.
La capacidad de generación de una máquina sı́ncrona está limitada por

la fuente de generación primaria entre otras restricciones de tipo técnico a
saber:
• Limite de potencia activa del primomotor.
• Corriente máxima soportada por el devanado de inducido.
• Corriente máxima soportada por la bobina de excitación.
• Lı́mite de estabilidad estática.
• Potencia activa mı́nima.
Estas restricciones influyen en la capacidad de generar potencia activa

y reactiva; la curva de capacidad describe estas restricciones permitiendo
un modelo más riguroso de la máquina sı́ncrona.
En su forma más simple, un generador de rotor liso es representado

como se muestra en la figura ?? de donde se obtiene la siguiente relación:
V = Ef − j · X · I (6.1)
85
jX
V
Ef I
Figura 6.1: Modelo simplificado de la máquina sı́ncrona
en donde:
V : Voltaje en terminales de la máquina (Puede ser tomado como 1pu)

Ef : Voltaje inducido, es proporcional a la corriente de campo por lo cual
una corriente de campo máxima corresponde a un Efmax
X : Inductancia de la máquina de eje directo.
I : corriente de lı́nea.
la potencia activa suministrada por la máquina (P +j ·Q = V ·I ∗ ) puede

ser reescrita de acuerdo con la ecuación 6.1 tomando como referencia el
voltaje V :
V2 Ef · V
−j· = −j · −P +j·Q (6.2)
X X
La dirección de −P es tomada en sentido contrario para permitir el

diagrama en el primer cuadrante. En la figura 6.2 se muestra la curva por
capacidad resultante en donde el área sombreada representa los puntos en
los cuales la máquina puede operar normalmente.
Las lı́neas verticales (P max y P min) representan los lı́mites de potencia

activa, el primero por efecto de la capacidad máxima de generación
Lím
ite
de
a rm
ad
ura
Lí
I
m
V
ite
de
X
ca
Ef V /
m
po
ilida d
V 2/X
b
E s ta
Pmin
Pmax
Figura 6.2: Curva de capacidad de un generador
(turbina) y el segundo por la mı́nima generación posible en el caso de las

plantas térmicas.
la curva de capacidad es construida para un valor determinado de V

2
(normalmente 1pu), por lo cual la distancia VX es conocida.
El lı́mite de armadura esta dado por la máxima corriente que puede

soportar los devanados del estator; esta limitante se transforma en un
circulo con centro en el origen y radio V · Imax .
La máxima corriente en el devanado induce una tensión máxima

·V
Efmax , esta restricción se representa por un circulo de radio Efmax
X
y
V 2
centro en el eje Y a una distancia X bajo el eje X.
El lı́mite de estabilidad estática es representado por un ángulo máximo

δ entre Ef y V , en teorı́a este ángulo es de 90o , sin embargo en la practica
87
se toman valores menores. En general existe un lı́mite inferior dado por

el calentamiento en los extremos de la armadura. Este lı́mite se determina
en forma experimental y corresponde aproximadamente a una lı́nea recta
desde el punto de factor de potencia 0.95 en adelanto sobre la curva lı́mite
de la armadura hasta el punto en que absorbe el 60% de la potencia
reactiva sin potencia activa.
La curva de capacidad corresponde a la superposición de todos los

lı́mites por lo cual no siempre presenta la misma forma, por ejemplo, en
algunos sistemas el lı́mite de potencia máxima puede ser muy superior
por lo cual no entrarı́a a acotar la región factible.
Una curva simplificada puede ser planteada en función de restricciones

lineales, no obstante la curva completa es la mejor representación para la
máquina.
Qmin ≤ Q ≤ Qmax (6.3)
Pmin ≤ P ≤ Pmax (6.4)
La curva de capacidad en una máquina sincrona puede ser utilizada

para entender el costo asociado a la potencia reactiva (costo de
oportunidad); aunque Q no tiene un costo directo de generación, si tiene
un costo asociado a la cantidad de generación activa que desplaza cada
M V Ar generado cuando se encuentra en uno de los lı́mites en la curva.
Una central trabajando en un punto en el interior de la curva está en

capacidad de generar 1 MW adicional de potencia activa requerida por
el sistema, si el primo-motor (turbina) lo permite, no obstante, cuando se
encuentra trabajando en un punto lı́mite de la curva1 es imposible generar
potencia activa adicional, -aunque la turbina lo permita- pues la potencia
reactiva entra como restricción.
El costo de oportunidad de la potencia reactiva es entonces el valor

1
Corriente máxima soportada por el devanado de inducido o Corriente máxima
soportada por la bobina de excitación.
equivalente en potencia activa desplazada por cada unidad de potencia

reactiva usada o en reserva.
6.2 Control de tensión.
Los generadores que alimentan la demanda del SIN están obligadas a

participar en el control de voltaje, por medio de la generación o absorción
de potencia reactiva de acuerdo con la curva de capacidad declarada en los
formatos de capacidad. La tensión puede ser controlada de forma directa
por medio de la inyección de reactivos, éstos pueden ser proporcionados
por diferentes elementos como generadores, condensadores sı́ncronos y
dispositivos basados en electrónica de potencia entre otros. No todos los
generadores tienen un sistema de control control de tensión o regulador
automático de tensión (AVR2 ), los sistemas de cogeneración por ejemplo
solo pueden mantener constante un factor de potencia por lo que la
tensión se hace variable; solo las centrales con capacidad de generación
superior a 1 MVA cuentan con sistemas AVR. En particular, el control de
tensión se realiza de acuerdo a tres niveles:
Control primario : es el control que se realiza localmente en las

subestaciones y plantas de generación en forma automática o
manual, sobre los recursos de compensación reactiva. El control
automático, se ejerce individualmente sobre cada uno de los equipos.
Control secundario : es el control automático que se realiza sobre el nivel
primario de control de un área de control, manteniendo el voltaje de
uno o varios nodos piloto dentro de una banda especificada por el
nivel terciario de control.
Control terciario es el control automático o manual que realiza el CND
con el objeto de coordinar la operación de las áreas de control, para
garantizar en lo posible la operación segura y eficiente del STN. El
nivel terciario de control determina los voltajes a ser fijados en los
nodos pilotos de cada una de las áreas de control.
2
Automatic voltage regulator
89
El control primario se realiza de forma local, por tal razón los

responsables de este control son los operadores de los equipos de
compensación reactiva. No obstante, algunos elementos del sistema
pueden ser controlados de forma centalizada desde el CND. El control
secundario y terciario debe ser efectuado de forma coordinada, por tal
razón el responsable directo es el operador de la red (CND).
En un nivel de control por encima del control terciario de tensión se

encuentran mecanismos de detección de de proximidad a un colapso de
tensión. En estos casos de niveles bajos de tensión puede ser necesario un
deslastre de carga3 .
Un colapso de tensión se presenta cuando hay disponibilidad de

potencia reactiva para abastecer los requerimientos de la demanda. En
este caso la tensión del sistema disminuye hasta el punto en el cual el
sistema no es controlable.
Niveles muy bajos en la tensión pueden ocasionar dificultades técnicas

en los diferentes elementos, por ejemplo en los transformadores y lı́neas
de transmisión, niveles de tensión bajos exigen un aumento en la corriente
para necesaria para transportar una potencia determinada, este aumento
en la corriente puede ocasionar a su vez una violación de los lı́mites de la
lı́nea o el transformador a tal punto de sacarlo de operación, ası́ mismo,
un aumento en la corrientes produce un aumento en las pérdidas técnicas.
En el caso de las cargas motrices un disminución en la tensión puede
ocasionar desconexión del elemento o mal funcionamiento del mismo.
Por el contrario, un aumento en las tensiones nodales puede ocasionar

la saturación de los transformadores, ası́ como problemas de aislamientos
en los diferentes elementos. Los equipos electrónicos y de control son
especialmente sensibles a los voltajes altos.
3
El deslastre de carga es la desconexión automática de carga ante eventos en el sistema
interconectado nacional o en parte de éste. Se puede efectuar como último recurso tanto
en el control de tensión como de frecuencia
6.3 Control automático de la generación

(AGC)
La demanda no es una variable estática, por el contrario, ésta

cambia a lo largo del dı́a de forma aleatoria; aunque el promedio de
demanda en una hora determinada es una variable predecible, el balance
carga-generación varia constantemente, por tanto existen diferencias entre
la generación real y la demanda real. El error en el pronóstico de demanda
oscila entre el 3 y el 5%.
Igualmente, los cambios en la topologı́a del sistema causados por

contingencias pueden modificar la demanda total. Esta variación en
la demanda provoca un desequilibrio entre la potencia generada y la
potencia demandada, el cual debe ser compensado para evitar efectos no
deseados sobre la frecuencia y los niveles de tensión.
Demanda Generación
60 Hz
59 Hz 61 Hz
Figura 6.3: Analogı́a de la relación frecuencia vs potencia activa
En general la relación entre la frecuencia y la potencia activa del

sistema se puede explicar por medio de una analogı́a relativa a una
balanza en donde se representa los recursos de generación vs las
necesidades del sistema como se muestra en la figura 6.3. Cuando la
generación total del sistema es capaz de abastecer totalmente la demanda
91
del sistema, incluidas las pérdidas de potencia activa, la frecuencia

permanece constante en el valor deseado de 60Hz. No obstante, frente a
una variación en la demanda o la generación, la frecuencia puede cambiar
a un punto cercano a la frecuencia nominal. Este rango de variación
permitida es denominado banda muerta.
Cuando la generación total del sistema es superior a la demanda del

mismo, la frecuencia tiende a aumentar debido a que la máquina sı́ncrona
se acelera al tener menor resistencia al giro. Igualmente, cuando la
demanda aumenta, la maquina sı́ncrona tiende a frenarse por lo que la
frecuencia disminuye.
Un frecuencia por fuera del valor nominal (60Hz) es una situación no

deseada debido a:
• La velocidad de los motores depende estrechamente de la frecuencia
• Algunos dispositivos de control pueden depender de la frecuencia
• La turbina puede presentar fallas a velocidades diferentes del rango

para el cual fueron diseñadas.
• El sistema puede perder estabilidad
Es relativamente simple definir el comportamiento de la frecuencia

frente a una variación de la carga, ésta se puede asimilar como una
balanza (ver figura ??) en donde la frecuencia aumenta en la medida
que se presente una mayor potencia generada con respecto a la potencia
demandada por el sistema. Este comportamiento es fácil de entender en
un sistema aislado, no obstante, en un sistema interconectado la situación
se hace mas compleja. En esta situación el sistema es separado en
sub-areas de control en donde si un área presenta déficit de potencia, otra
área puede entrar a dar asistencia para mantener la frecuencia constante.
Esto introduce un mayor grado de complejidad al problema, pues dada
una desviación de la frecuencia del sistema en estado permanente que a
su vez es causada por una desviación en el balance de potencias generadas
y demandas es necesario definir si ésta es causada por la propia área o por
otras áreas de control.
Para ello se define el Intercambio total de la red dado como la

suma algebraica de todos los intercambios de potencia entre dos áreas
adyacentes. Es evidente que el intercambio total de la red difiere entre
lo planeado y lo real.
6.3.17 Control El control primario de la velocidad es de carácter local y se refiere al

primario de gobernador de velocidad:
frecuencia
2 4
ωm = · ωe = · π · f (6.5)
p p
en donde ωm es la frecuencia mecánica de la turbina mientras que ωe es

la frecuencia eléctrica.
El gobernador de velocidad en una central hidráulica o térmica tiene

conceptualmente un funcionamiento similar. La mayor parte de éstos
sistemas son de carácter hidráulico-mecánico o electro-hidráulico. El
modelo del sistema puede ser descrito de la siguiente forma4 :

KG 1
ΔXE = ΔPc − · Δf (6.6)
1 + TG · S R
A su vez, el modelo de la turbina es descrito por medio de la ecuación
KT
ΔPM = · ΔXE (6.7)
1 + TT · S
Reemplazando estas dos ecuaciones se tiene lo siguiente

KT · KG 1
ΔPM = · ΔPc − · Δf (6.8)
(1 + TG ·) · (1 + TT · S) R
4
Ver referencia [2]
93
Usando el teorema del valor final
lim f (t) lim S · F (S)

= (6.9)
t→∞ S→0
y definiendo adecuadamente la constantes del problema, se tiene:
1
ΔPM = ΔPc − · Δf (6.10)
R
Esta relación está definida en el dominio del tiempo y en estado

permanente. Para el caso no controlado (ΔPc ), se le permite al sistema
modificar la frecuencia de la siguiente forma
1
ΔPM = − · Δf (6.11)
R
donde R es llamada regulación de velocidad.
Ejemplo 6.1
En la figura 6.4se muestra un sistema no controlado directamente, esto

significa que si se presenta una variación en la potencia demanda ΔPM la
frecuencia cambia inmediatamente a un valor Δf . Mostrar la respuesta en
frecuencia si en este caso las constantes de tiempo de la turbina y el generador
son despreciables en comparación con las constantes de tiempo del sistema.
La respuesta en frecuencia frente a una variación de 0.1pu en el instante

de tiempo t = 2s se presenta en la figura 6.5.
Usando el resultado obtenido a partir del teorema del valor final, y

considerando la constante de tiempo del sistema, se calcula para una
variación en potencia ΔPM = 0.01pu un cambio en resultante en la
Figura 6.4: Diagrama de bloques de un sistema no controlado
frecuencia de Δf = −0.0235. El lector puede comprobar igualmente el

resultado del sistema sin tener en cuenta el modelo de la turbina y el
control de velocidad (los cuales presentan una constante de tiempo muy
pequeña en comparación con la constante de tiempo del sistema).
6.3.18 Control El gobernador básicamente puede controlar un sistema aislado

secundario modificando la potencia activa frente a una variación en la frecuencia.
En un sistema interconectado, los cambios en la frecuencia deben ser
asumidos de forma coordianda por las diferentes centrales del sistema
agrupadas en áreas de control. A esto se denomina control secundario
de frecuencia.
En el caso del sistema colombiano, la generación fuera de mérito

necesaria para suplir los servicios de regulación secundaria de frecuencia,
es comercialmente, obligación de los generadores en proporción a la
potencia despachada.
A medida que la demanda global del sistema sufre variaciones con

respecto a su valor nominal se presentan desviaciones en las variables de
estado. El sistema de control debe detectar dichas variaciones y tomar
95
0.005
0.01
0.015
0.02
0.025
0.03
0.035
0 1 2 3 4 5 6 7 8 9 10
Figura 6.5: Respuesta en frecuencia del sistema sin cambio en la potencia

mecánica
en tiempo real las medidas correctivas que eliminen en el menor tiempo

posible estas desviaciones.
Las variables de estado son continuamente monitoreadas, algunas,

como la magnitud de la tensión, son monitoreadas de forma directa
mientras que otras, como el ángulo de la tensión, presentan una mayor
dificultad, razón por la cual deben ser determinadas indirectamente5 .
Las mediciones obtenidas son comparadas con las referencias deseadas
obteniendo un conjunto de señales de error.
El equipo de control debe responder automáticamente a estas señales

de error y ajustar la salida de potencia para mantener la frecuencia
constante (60Hz) como se muestra en la figura 6.6.
5
Por ejemplo a través de la potencia activa que es función de los ángulos en las
fv
fn =60
frecuencia (Hz)
P Potencia (MW)
Figura 6.6: Regulación de frecuencia
Como se mostró en la sección anterior, la pendiente decreciente que

relaciona la frecuencia con la potencia activa es denominada regulación
de velocidad (ecuación 6.11). Para normalizar esta relación, se define la
regulación unitaria de velocidad:
(fv − f ) /fn
Ru = (6.12)
P/Sn
en donde:
fv : frecuencia en vacı́o
f : frecuencia actual
fn : frecuencia nominal (60Hz)
P : potencia actual
Sn : potencia nominal
tensiones nodales
97
Ası́ que en cada unidad de generación se tiene
Snk Δf
ΔPk = − · (6.13)
Ruk fn
En un sistema interconectado, el control automático de la generación

se refiere especı́ficamente al control secundario de frecuencia. El
sistema es dividido en áreas de control encargadas de dar respuesta
a las fluctuaciones de la demanda propia. Desde luego puede existir
intercambios entre áreas del sistema, sin embargo cada área debe estar
en capacidad de absorber sus propios cambios de carga (figura 6.7).
Figura 6.7: Zonas de control en el sistema Colombiano
En la operación real, las áreas de control son encargadas de asegurar

también los intercambios de potencia entre áreas. Una diferencia entre el
flujo real de potencia y el flujo programado indica un error o desviación
de la generación en cada área. Éste error denominado ECA (error de
control de área) representa la discrepancia entre la generación las fuentes
de generación y la demanda total dentro del área de control.

ECA = Pkm − Pkm(esp) − (10 · β) · (f − fn ) (6.14)
donde
Pkm : Flujo de potencia actual por la lı́nea km
Pkm(esp) : Flujo de potencia programada en la lı́nea km
β : Constante de sesgo de frecuencia
f : frecuencia actual del sistema
fn : frecuencia nominal del sistema 60Hz
La constante de sesgo de frecuencia β es un estimado de la respuesta

en frecuencia del área de control. El valor de β es una aproximación de la
caracterı́stica de respuesta en frecuencia (CRF) la cual, como su nombre lo
indica, representa la respuesta del sistema a un cambio en la frecuencia.
La CRF es la respuesta combinada de todas las unidades de generación y
las cargas motrices frente a una variación en la frecuencia. La respuesta en
frecuencia depende principalmente de:
• Los parámetros de control de los diferentes gobernadores de

velocidad
• La condición del sistema de potencia cuando una desviación de

frecuencia ocurre.
• La respuesta en frecuencia de las cargas conectadas.
La CRF varia con respecto a las condiciones de operación, esto significa

que una pérdida de generación determinada puede producir respuestas en
frecuencia diferentes si ocurre en dos puntos de la curva de carga distintos.
En la gráfica ?? se muestra un ejemplo real de la comportamiento de la
frecuencia frente a la salida de una unidad generadora.
En el interior de cada área los reguladores de velocidad deben actuar

de forma coordinada y determinar los cambios a realizar frente a una
variación en la demanda ΔP .
99
60.06
60.03
60
Frecuencia [Hz]
59.97
59.94 Banda muerta propuesta
59.91
El AGC Retorna el sistema a 60Hz.
59.88
Respuesta de los gobernadores
59.85
F.R =220MW/0.2Hz de velocidad.
59.82 =1100 MW/Hz
59.79
f = 59.8 Hzy t = 18:45:16 Tiempo[min ]
59.76
18:51:50
18:44:38
18:47:31
18:46:05
18:48:58
18:50:24
18:41:46
18:43:12
Figura 6.8: Evolución de la frecuencia ante la pérdida de 220 MW de

Guavio.
ΔP = ΔPk (6.15)
donde el cambio en la potencia generada ΔPk esta relacionado con la

frecuencia por medio de la relación
Snk Δf
ΔPk = − · (6.16)
Ruk fn
el cambio total es entonces

Δf Snk
ΔP = − (6.17)
fn R uk
El cambio en la frecuencia total es común a todo el sistema

Δf −ΔP
= S (6.18)
fn nk
R uk
De donde se obtiene el cambio en cada cental
Snk
R uk
ΔPk = S (6.19)
nk
R uk
Ejemplo 6.2
Cuatro unidades de generación operan en paralelo a una frecuencia de 60Hz

para abastecer una demanda de 1300M W . Los datos nominales y los despachos
de estas unidades son mostrados en la tabla 6
Tabla 6.1: Valores nominales para el ejemplo

unidad S[M V A] nominal P [M W ] despachada R%
A 300 150 3
B 500 400 4
C 500 350 3
D 600 400 5
Súbitamente, el generador A sale de operación, determinar los cambios en la

frecuencia del sistema y las variaciones en la potencia activa generada por las
demás unidades antes de que ocurra cualquier acción de control suplementario.
El cambio en la potencia activa es ΔP = 150, con este valor y los datos

nominales de cada generador, se calcula la variación en la frecuencia.
Δf −150
= 500 500 600 = −3.644 × 10−3 (6.20)
fN 0.04
+ 0.03
+ 0.05
101
Esto significa una variación de −0.2186Hz en la frecuencia, el signo

negativo indica una disminución de la misma pues el efecto de la salida de
un generador es equivalente al incremento en la demanda. Esta variación
es superior a la máxima permitida dada por la banda muerta, por tal razón
es necesario realizar cambios en la generación. Usando la relación 6.16 se
tiene:
• ΔPB = 45.55M W por lo que la unidad B debe aumentar su

generación a PB = 445.55M W
• ΔPC = 60.73M W por lo que la unidad C debe aumentar su

• ΔPD = 43.72M W por lo que la unidad D debe aumentar su

6.3.19 Control
Se conoce con el nombre de regulación terciaria de frecuencia la Terciario
re-localización de los recursos en forma económica una vez se se hayan
finalizado las acciones de control automático. La re-localización de
recursos secundarios se realiza con el fin de mantener la frecuencia
dentro de la banda muerta y hacer frente a un desequilibrio inesperado
y transitorio entre la generación y la demanda. Ésta reubicación de
recursos de generación puede ser costosa y económicamente inviable, por
tal razón, una vez superada la etapa preliminar de los controles primario
y secundario, es necesario hacer ajustes más finos teniendo en cuenta
criterios de costos.
Tı́picamente los recursos de regulación terciaria se obtienen después de

la actuación de la regulación secundaria de frecuencia, es decir después de
10 minutos hasta 1 hora.
Finalmente, ante desbalances mayores por efecto de contingencias,

existen mecanismos adicionales a saber:
EDAC Esquema de desconexión automática de carga: Es el conjunto de

protecciones que actúan coordinadamente bajo un evento de pérdida
de generación de magnitud en el sistema previniendo un colapso
total mediante la desconexión de carga en cantidades especı́ficas y en
tiempos particulares. También es denominado eyección automática
de carga.
RAG Rechazo automático de generación: Se refiere a las acciones
de control en los generadores para disminuir la potencia activa
generada frente a una disminución brusca y considerable en la
demanda y de esta forma mantener la frecuencia en 60Hz.
6.3.20 Reserva Un generador puede modificar su potencia activa generada frente

rodante a una modificación de la frecuencia solo si tiene una capacidad de
generación adicional. Esta capacidad de generación es generalmente
denominada reserva rodante. La reserva rodante es la diferencia entre la
salida actual de potencia con respecto a la máxima capacidad disponible.
Un generador puede tener un valor grande de reserva rodante y sin
embargo no presentar una buena respuesta frente a una variación en la
demanda. Esto se debe a la velocidad de respuesta de la unidad, por tal
razón gran parte de la reserva rodante se le asigna a las plantas hidráulicas.
En general las unidades hidráulicas pueden pueden responder mas rápido
y con mayor disponibilidad de potencia; una central hidráulica puede
entregar 100 MW adicionales en menos de 10 segundos. Por el contrario,
las plantas térmicas con ciclo de vapor tienen una repuesta lenta, por
ejemplo una central de 500 MW funcionando a 300 MW puede generar
solo 20 MW adicionales frente a un requerimiento del sistema en un
tiempo de 10 segundos.
La frecuencia como toda variable real tiene un rango de error tolerable,

este error está asociado principalmente a los sistemas de control mecánicos
que impiden eliminar totalmente el error. Normalmente, la frecuencia en
el sistema puede mantenerse en valores de 60Hz ± 0, 02Hz, a este margen
de error se le denomina Banda muerta.
La reserva rodante se refiere a la capacidad disponible que puede

103
ser utilizada inmediatamente para suplir potencia en el evento que otra

unidad falle o ocurra un requerimiento inesperado de potencia. Las
turbinas generadoras para tales unidades generalmente necesitan estar en
funcionamiento dado que la reserva de otra unidad disponible necesita un
tiempo prudencial para llevar a su plena potencia y sincronizar con la red.
Los métodos para modelar la reserva rodante varı́an dependiendo de las
unidades que gobernarán.
Existen dos criterios determı́nisticas básicos para calcular la reserva

rodante en un sistema de generación:
Porcentaje de reserva : Según este criterio una reserva de entre el 15%

y el 35% es suficiente para un sistema. El porcentaje de reserva es
calculado como:
CG
R[%] = − 1 · 100% (6.21)
Dm
en donde CG es la capacidad de generación instalada y Dm es la
máxima demanda esperada.
Pérdida de la unidad más grande: Según este criterio la reserva del
sistema debe ser mayor o igual a la capacidad de la mayor unidad
del sistema.
6.4 Lecturas complementarias
1. Rodas, Dario Eliécer; D̈iagrama de capacidad de un alternador ¨ En:

Revista Scientia et technica. Año IX Nº 22 Oct. 2003.
2. Giraldo Buitrago, Didier; Tabares Gomez, Ivan; Teorı́a de control.
Universidad Tecnológica de Pereira. Facultad de Ingenierı́a
Eléctrica. 1997.
3. John J. Grainger; William D. Stevenson Análisis de sistemas de potencia.
McGraw-Hill. 1996.
4. Olle, Elgerd. Electric energy systems theory: an introduction

McGarwHill.
5. EPRI. Interconnected Power System Dynamics Tutorial.

Parte II
Aspectos económicos de los

sistemas de generación
105
Introducción a la Optimización
7
Contenido
Investigación operativa.
Multiplicadores de Lagrange.
Programación Dinámica.
Programación Lineal
107
108 CAPÍTULO 7. INTRODUCCIÓN A LA OPTIMIZACIÓN
L Ainvestigación operativa es una disciplina moderna caracterizada por

la utilización de teorı́as, métodos y algoritmos matemáticos aplicados a la
toma de decisiones.
El uso de técnicas de investigación operativa permite resolver los

problemas inherentes a la ingenierı́a eléctrica teniendo en cuenta los
aspectos técnicos y económicos del sistema. Para aplicar las metodologı́as
de solución, es necesario plantear el problema como un modelo
matemático de optimización, caracterizado por una función objetivo y un
conjunto de restricciones como se muestra a continuación:
min f (x) (7.1)

sujeto a
g(x) − a = 0 (7.2)
Durante la segunda guerra
mundial, un grupo de h(x) − b <0 (7.3)
investigadores militares
comenzaron a utilizar
técnicas matemáticas para
la ubicación óptima de
radares, polı́ticas de Como su nombre lo indica, la función objetivo f (x) representa el fin
maniobrabilidad de los
barcos y submarinos entre deseable en el proceso de toma de decisiones; esta función corresponde, en
otros. Este es el inicio
formal de la investigación la mayorı́a de los casos, a una variable económica a optimizar, por ejemplo
de operaciones
minimizar los costos o maximizar las utilidades.
En estos años posteriores a

la Segunda Guerra
Mundial, en Estados
Las restricciones corresponden a un conjunto de ecuaciones g(x) e
Unidos se encontró que la
eficaz coordinación de
inecuaciones h(x) que la solución óptima debe cumplir. Normalmente,
todas las energı́as y
recursos de la nación era las restricciones representan las variables de tipo técnico a tener en cuenta
un problema de tal
complejidad, que su en el sistema, por ejemplo las leyes de Kirchhoff, los lı́mites de tensión y
resolución pasaba
necesariamente por los potencia reactiva, las ecuaciones de flujo de potencia entre otros.
modelos de optimización.
Se ha estimado, de una
manera general, que si un
paı́s subdesarrollado
Tanto las restricciones como la función objetivo tienen asociadas un
utilizase los métodos de
optimización en sus
conjunto de variables independientes x que pueden ser modificadas
procesos, su producto
interior bruto (PIB) en el sistema para optimizar la función objetivo cumpliendo con las
aumentarı́a entre un 10 y
un 15% en tan sólo un año restricciones, éstas son denominadas variables de decisión.
109
Si tanto la función objetivo como las restricciones son ecuaciones

lineales y las variables de decisión toman valores continuos se dice que el
problema es de programación lineal (PL). De igual forma, si las variables
solo pueden tomar valores enteros el modelo se hace más complejo y debe
ser solucionado por técnicas de programación entera (PE) y programación
entera mixta (PEM). Finalmente, si la función objetivo o alguna de las
restricciones presenta una no linealidad, el problema es de programación
no lineal (PNL) o programación no lineal entera (PNLE) de acuerdo a la
caracterı́stica continua o discreta de las variables de decisión.
Para problemas de optimización no-lineal con restricción de igualdad

exclusivamente, la técnica más utilizada es la función Lagrangiana o
multiplicadores de Lagrange; esta metodologı́a además de dar solución
a este tipo de problemas, permite una interpretación geométrica y
económica que por su importancia en los sistemas de generación será
tratada inicialmente. Otra técnica importante es la programación
dinámica, utilizada en la solución de algunos problemas no lineales y
enteros de caracterı́sticas especiales como se verá mas adelante.
Los problemas de PL suelen ser solucionados por dos técnicas

fundamentales: el algoritmo simplex y el método de puntos interiores.
7.1 Multiplicadores de Lagrange
La forma más simple y comúnmente utilizada para encontrar puntos

extremos (máximos o mı́nimos) a una función f (x) es determinar aquellos
puntos x en que f (x) = 0. No obstante, en muchos problemas prácticos, el
modelo matemático no es totalmente definido por una función f (x) pues
debe cumplir además con una restricción g(x) = a. En estos casos debe ser
utilizado el método de los multiplicadores de Lagrange:
∇f (x) = λ · ∇g(x) (7.4)
Este método es esencialmente una extensión natural al método usual

en donde se considera una nueva función (Lagrangeana) dada por:
L(x) = f (x) + λ · (g(x) − a) (7.5)
en donde es claro que minL = minf , ası́ basta encontrar los valores de
x (variables de desición) y λ (multiplicador de lagrange)1 .
La función f (x) determina el grado de optimalidad de las variables x

mientras que la función g(x) determina la factibilidad. Una desviación de
En los siglos XVII y XVIII,
grandes matemáticos g(x) con respecto a la contante a no es aceptada como solución aún cuando
la función f (x) sea mı́nima, en este sentido, la función g(x) − a = 0 actúa
como Newton, Leibnitz,
Bernouilli y, sobre todo,
Lagrange, que tanto
habı́an contribuido al como una restricción al problema.
desarrollo del cálculo
infinitesimal, se ocuparon
de obtener máximos y
mı́nimos condicionados Un problema de optimización con una restricción de igualdad puede
de determinadas
funciones. ser representado como se muestra a continuación:
min f (x) (7.6)

sujeto a
g(x) − a = 0 (7.7)
Una función lagrangeana consiste en adicionar la restricción a la

función objetivo penalizando la infactibilidad por medio de una valor λ
como se muestra en la ecuación 7.8
min L(x) = f (x) + λ · (g(x) − a) (7.8)
Para minimizar esta nueva función se recurre a las llamadas

condiciones de primer orden:
1
La condición del método es una condición necesaria pero no suficiente
111
∂L ∂f ∂g
= +λ· =0 (7.9)
∂xi ∂xi ∂xi
∂L
= g(x) − a = 0 (7.10)
∂λ
Sea x el vector de valores x que cumple con estas condiciones. este

vector corresponde a una solución óptima del sistema.
7.1.21
Los multiplicadores de Lagrange tienen una importante interpretación Interpretación
económica que debe ser estudiada. económica
La solución óptima x será distinta si los recursos existentes en la

restricción a cambian, por lo cual x = x(a); reemplazando en la función
lagrangeana se obtiene:
L(x(a)) = f (x(a)) + λ(a) · (g(x(a)) − a) (7.11)
Derivando respecto a a se tiene:

∂L(x(a)) ∂x ∂x ∂λ
= ∇f · + λ · ∇g · −1 + · (g − a) (7.12)
∂a ∂a ∂a ∂a
Puesto que x(a) es una solución debe cumplir las condiciones de

primer orden (7.9 y 7.10) de donde se obtiene:
∂L(x(a))
= −λ(a) (7.13)
∂a
El multiplicador de Lagrange puede ser expresado como la derivada

de la función objetivo; lo anterior es causado por la equivalencia entre la
función objetivo y el lagrangeano en el punto x.
∂f (x(a))
λ(a) = − (7.14)
∂a
El multiplicador de Lagrange puede interpretarse como un valor de

sensibilidad de la función objetivo ante un relajamiento de la restricción,
es decir, si se permite a la restricción g(x) tomar un valor a + Δa la
función objetivo podrá tomar un valor menor en aproximadamente λ · Δa
unidades. Comúnmente se denomina precio sombra al valor de λ.
7.2 Programación dinámica
Múltiples problemas en ingenierı́a presentan un comportamiento

dinámico en el tiempo; en este tipo de problemas la solución es dada por
decisiones secuenciales en etapas sucesivas con el fin de minimizar el costo
total de dichas decisiones. La programación dinámica permite resolver
estos problemas.
En los problemas dinámicos, cada etapa se valora no solo por el costo

actual de tomar una decisión sino también por el costo futuro que se
origina a partir de ella.
Un problema dinámico puede ser representado por un grafo orientado

El objetivo es encontrar el camino óptimo en función de los costos C

para ir desde el punto inicial A hasta el punto F . El procedimiento para
encontrar la solución puede realizarse hacia adelante (forward) o hacia atrás
(backward) usando el principio de optimalidad de Bellman2 :
La decisión en un estado depende del estado en el cual está,

no de como se llegó a él pues toda la información sobre el
pasado se resume en el estado mismo. Todo subconjunto de
2
R. Bellman: Dynamic Programming, Princeton University Press, 1957.
113
C
C BC
C CF
B
C AB C DC F
C CE
A C AD
D
C EF
C DE
Figura 7.1: Grafico orientado
una solución óptima es a su vez una solución óptima para un

problema parcial.
En el caso mostrado en la figura, si se sabe que el camino óptimo para

ir de A a F pasa por B, entonces el problema parcial de ir de B a F está
contenida en la solución óptima global.
7.2.22
Sea k una etapa arbitrarı́a en el grafo que representa el problema, xk Algoritmo de
son las variables de decisión asociadas a dicha etapa y fk (xk ) el valor solución
acumulado de la función objetivo desde el estado k hasta el estado final
N . Se tiene que:
fk (xk ) = min {Cxk + fk+1 (xk+1 )} (7.15)
El algoritmo base consiste en agrupar soluciones a cada uno de los

problemas parciales hasta conseguir una solución al problema global.
Ejemplo 7.1
Encontrar la ruta mas económica para ir desde el estado A hasta el estado J

dados el siguiente grafo:
4 12
B E H
15 3
7 10 8
5
A C G
J
3
11 5
6 3
D 4 8
F I
Figura 7.2: Ejemplo de programación dinámica
En este caso se puede optar por soluciones constructivas3 que en este

caso consiste en elegir los caminos mas económicos sin tener una visión
global del problema. La solución constructiva haciendo un recorrido hacia
adelante es ACEHJ con un costos de 30 mientras que la solución haciendo
un recorrido en sentido contrario serı́a ADF GHJ con un coste de 24. La
solución utilizando un algoritmo basado en el principio de optimalidad
de Bellman se presenta a continuación.
Se toma en cada etapa el costo mı́nimo teniendo en cuenta el estado y

la historia del mismo:
3
Este tipo de algoritmos son comúnmente denominados algoritmos golosos y
corresponden solo a una primera aproximación a la solución del problema
115
Etapa 1:
Estado Coste Minimo Origen

B 7 A
C 5 A
D 6 A
Etapa 2:

E 11 B
G 13 F
F 10 D
Etapa 3:

H 21 G
I 17 G
Etapa 4:

J 22 I
La solución óptima tiene un coste de 21 mediante una ruta

que puede ser determinada por medio de un recorrido en sentido
inverso por las diferentes etapas JIGF DA. Como conclusión, los
algoritmos constructivos o ”miopes” presentan un grado de optimalidad
visiblemente inferior a las soluciones encontradas por medio de la
programación dinámica.
En ingenierı́a eléctrica, las técnicas de programación dinámica son

ampliamente utilizadas en problemas como:
1. Planeamiento de la expansión de la generación.

2. Programación de plantas térmicas.
3. Coordinación hidrotérmica.
7.3 Programación lineal
Aunque la mayorı́a de problemas en ingenierı́a son de carácter

no-lineal, frecuentemente son linealizados al rededor de un punto de
operación. Para resolver este tipo de problemas es necesario recurrir a
las teorı́as de programación lineal.
Un problema de programación lineal, al igual que todos los problemas

de optimización, es representado por una función objetivo y un conjunto
de restricciones, solo que en este caso una y otra es de carácter lineal como
se muestra en la ecuación 7.16
min cT · x (7.16)
sujeto a
A·x<b (7.17)
El matemático fránces Jean

Baptiste-Joseph Fourier
(1768-1830) fue el primero
c es un vector columna cuyos coeficientes representan los costos
en intuir, aunque de forma
imprecisa, los métodos de
asociados a cada una de las variables de decisión xi , la matriz A y el
lo que actualmente
llamamos programación
vector de recursos b representan las restricciones del sistema para formar
lineal y la potencialidad
que de ellos se deriva
un polı́tope convexo en el cual la solución óptima se presenta en un vértice.
Existen dos metodologı́as básicas para dar solución a este

tipo de problemas: el método simplex y los métodos de puntos
interiores. Igualmente existen metodologı́as enfocadas a problemas
con caracterı́sticas especiales como es el caso de flujo en redes.
”The tremendous power of

the simplex method is a
constant surprise to me”
El método simplex se basa en un principio según el cual la solución
G.B. Dantzig óptima a un problema de optimización lineal se encuentra localizada en
un vértice. Este método parte de un vértice dentro de la región factible del
117
problema y recorre vértices vecinos de mejor calidad (dada por la función

objetivo) hasta llegar a la solución óptima global.
El método de puntos interiores plantea una nueva función que tiene en

cuenta tanto a la función objetivo como a la desviación de las restricciones,
de forma análoga a la función lagrangeana. La metodologı́a parte de un
punto interior a las restricciones de desigualdad del sistema y recorre la
región factible en busca de la solución óptima. En el caso particular en que
el número de variables de decisión sean 3 o menos, el problema puede ser
representado de forma gráfica permitiendo una solución intuitiva. El método simplex fue
desarrollado en 1951 por
G.B. Dantzig en el United
States Bureau of Standards
SEAC COMPUTER,
ayudándose de varios
modelos por computador
Ejemplo 7.2 de la firma IBM.
Encontrar la solución óptima al siguiente problema de programación lineal.
max z = x + y (7.18)
sujeto a
x−2·y <0 (7.19)
x<4 (7.20)
y<4 (7.21)
x>0 (7.22)
y>0 (7.23)
(7.24)
El conjunto de restricciones puede ser representado como se muestra

en la figura 7.3
La región sombreada representa el conjunto de puntos factibles, es

decir, los puntos que cumplen con el conjunto de restricciones. La forma
de encontrar la solución óptima consiste en proyectar las curvas de nivel
y
y< 4
4
C
ur
z= 8
va
x> 0
de
ni
1) x< 4
ve
te
,
(1
en
l
di
ra
G
x - 2 y <0
0 x
y> 0 4
Figura 7.3: Solución gráfica al problema de PL
de z en el plano xy. La función z = x + y puede crecer indefinidamente

en la dirección (1, 1) mientras se cumplan las restricciones, es por ello
que la solución óptima se encuentra en un punto extremo del espacio de
soluciones (en este caso el punto (4, 4) con una función objetivo z = 8).
7.4 Algoritmos de búsqueda aleatoria
Además de las técnicas exactas como la programación lineal y los

multiplicadores de lagrange, existen técnicas que buscan aproximarse a
la solución optima de un problema de optimización mediante un conjunto
de pruebas aleatorias. Algunas de estas técnicas denominadas heurı́sticas
pueden incluir elementos de inteligencia artificial como los algoritmos
119
evolutivos y el recocido simulado. No obstante, para problemas de

relativa poca complejidad, es posible utilizar un algoritmo de búsqueda
estrictamente aleatoria, este algoritmo puede solucionar un problema del
tipo:
min F = {f (x) con x ∈ S} (7.25)
en donde S es el espacio de soluciones factibles, en otras palabras,

el conjunto de soluciones x que cumplen con las restricciones. Cabe
anotar que este tipo de metodologı́a no exige ningún tipo de caracterı́stica
especial en la función objetivo o el espacio de soluciones tal como
linealidad, convexidad, continuidad o posibilidad de hallar una derivada.
Es por esto, que estos algoritmos son altamente utilizados en diferentes
problemas de ingenierı́a.
El algoritmo básico genera de forma aleatoria e iterativa, un conjunto

de soluciones factibles x las cuales son evaluadas mediante la función
objetivo f para determinar la solución mı́nima. Desde luego, la solución
encontrada será tan solo una aproximación a la solución real, pero en
algunos casos, este tipo de solución aproximada es suficiente para la
aplicación correspondiente.
Sea ω = dim (S) la dimensión del espacio de soluciones, y τ la

dimensión de un conjunto s lo suficientemente cercano a la solución
óptima x por lo cual un valor x ∈ s es aceptado como solución
aproximada. Desde luego, el conjunto s no es conocido, pero es posible
definir la dimensión τ del mismo, la cual es un tipo de tolerancia del
algoritmo.
La probabilidad de que una solución x este en el espacio solución s es:
τ
p∗ = (7.26)
ω
Como cada evento es independiente, después de n iteraciones la

probabilidad de encontrar un x ∈ s es dada por:

N
P =1− (1 − p∗ ) (7.27)
k=1
por lo cual el número de iteraciones necesarias para obtener de forma

aleatoria una solución en un conjunto que contiene a la solución óptima
con una tolerancia y una probabilidad P es dado por:
ln (1 − P )
n= (7.28)
ln (1 − τ /ω)
Desde luego, este numero de iteraciones se convierte en una cota

superior la cual puede ser disminuida adicionando elementos especı́ficos
del problema particular, como la dirección del gradiente, o elementos de
inteligencia artificial como el concepto de selección natural. (este tipo de
metodologı́as va mas allá de los alcances de este libro).
Finalmente, se puede demostrar un mejor desempeño cuando

se utilizan direcciones aleatorias en cambio de puntos aleatorios,
especialmente cuando el espacio de soluciones tiende a ser convexo.
1. James C. Spall. Introduction to stochastic search and optimization.

Wiley-Interscience 2003.
2. http://departamentos.unican.es/macc/personal/profesores/castillo/Libro.htm
3. http://www.optimization-online.org/
4. http://valle.fciencias.unam.mx/intermat/ArticuloLag/articuloPDFb.pdf
5. http://w3.cnice.mec.es/eos/MaterialesEducativos/mem2003/programacion/
Operación de sistemas térmicos.
8
Contenido
Despacho económico
Enganche de unidades térmicas
Programación dinámica
121
122 CAPÍTULO 8. OPERACIÓN DE SISTEMAS TÉRMICOS.
A UNQUE el sistema colombiano presenta un alto porcentaje de

generación hidráulica, el estudio de los sistemas puramente térmicos
permite entender de forma inicial la aplicación de los modelos de
optimización a los sistemas de generación de energı́a eléctrica. Ası́ mismo,
la creciente importancia de los sistemas de cogeneración
Existen múltiples metodologı́as usadas para analizar la operación

de un sistema eléctrico de potencia. Algunas, como el flujo de carga,
permiten simular la operación de la red a partir de datos como la potencia
demandada, los niveles de tensión en los generadores y las potencias
generadas por los mismos (con excepción del nodo slack).
El problema de flujo de carga parte de dos datos en las centrales de

generación a saber: potencia activa generada y voltaje en terminales de
la máquina sı́ncrona. En este capı́tulo se mostrará cómo determinar la
potencia generada por parte de las centrales térmicas teniendo en cuenta
aspectos técnicos y económicos, ası́ mismo, se mostrará la relación con el
flujo de carga y las pérdidas del sistema. Se mostrarán dos problemas
clásicos: el despacho económico y la programación del arranque y parada
de las centrales térmicas.
8.1 Principales caracterı́sticas de las plantas

térmicas
La potencia eléctrica generada en una central térmica puede ser

determinada como función de las entradas de combustible, esta función
tiene en cuenta la eficiencia de la planta y el tipo de combustible utilizado.
Puesto que los costos variables de operación de la central están asociados
principalmente al combustible, es posible obtener una función de costo
frente a la potencia eléctrica generada.
Esta función puede ser llevada a una forma cuadrática o lineal por
tramos, en cualquier caso es sujeta de ser involucrada en un proceso de
optimización como parte de la función objetivo.
123
Estas funciones son aplicables tanto a las grandes centrales de

generación como a los pequeños cogeneradores.
El término cogeneración se aplica usualmente a los sistemas de

producción de energı́a y vapor para un proceso industrial. En este tipo
de sistemas, el excedente de energı́a puede ser inyectado a la red por lo
cual debe ser tenido en cuenta en los diferentes análisis.
En general, las decisiones operativas tomadas en los sistemas

exclusivamente térmicos durante una etapa de la programación no tienen
mayores consecuencias en los perı́odos subsiguientes. Sin embargo,
cabe señalar que esta afirmación pierde validez ante ciertos tipos de
contrato de compra de combustible que pueden restringir el consumo
del mismo, ası́ como cuando se considera de forma mas detallada
las restricciones de arranque y parada de las centrales ası́ como las
restricciones medioambientales de emisiones de gases contaminantes.
A diferencia de las plantas hidráulicas, los costos de operación de

una unidad termoeléctrica no dependen de la energı́a entregada por
otra unidad, esta caracterı́stica facilita el modelamiento matemático
permitiendo separar variables y procesos.
8.1.23
Las principales restricciones de las plantas térmicas se presentan a Restricciones
continuación:
Lı́mites de operación (Pmin < P < Pmax ): la potencia máxima

y mı́nima de un central térmica están dados por la capacidad del
primomotor y del lı́mite de estabilidad respectivamente. La potencia
mı́nima es del 10 al 30% para las plantas a gas y del 20 al 50% para las
plantas a carbón.
Mı́nimo tiempo de operación: una vez la planta es puesta en

marcha, ésta no puede apagarse repentinamente pues requiere un
tiempo que asegure temperatura constante en toda la unidad generadora,
particularmente en la turbina. De esta forma se controla la fatiga del
material, manteniendo los gradientes de temperatura dentro de los lı́mites
técnicos.
Mı́nimo tiempo de apagado: al igual que en la operación, se requiere

un tiempo que asegure temperatura constante cuando la unidad es puesta
fuera de servicio.
Limitaciones en unidades de planta: cuando se tienen múltiples

unidades generadoras en una misma planta pueden presentarse
restricciones de funcionamiento simultaneo de las mismas por
condiciones de mantenimiento o al compartir partes del proceso.
Rampas de encendido y apagado: debido a los gradientes de

temperatura, las unidades generadoras deben tomar o dejar carga eléctrica
de forma gradual, esta caracterı́stica puede ser modelada por medio de
rampas de encendido y apagado.
Limitaciones de combustible: la disponibilidad del combustible

puede variar según los tipos de contrato, las restricciones de transporte
del mismo, entre otros. Este hecho debe ser llevado en cuenta en el modelo
matemático.
8.1.24 Costos de Los costos totales de generación se clasifican en tres grandes grupos:
operación
• Costos asociados al capital deinversión.
• Costos de combustible.
• Costos de operación y mantenimiento.
Ası́ mismo, se pueden presentar dos categorı́as básicas: costos fijos y

costos variables.
Los costos fijos son independientes de la generación de la planta

e incluyen depreciación, intereses sobre el capital de inversión,
mantenimiento, empleados, entre otros. Los costos variables por
125
los costos asociados a la operación de las unidades de generación

térmica son los siguientes:
Costo de partida: Éste está asociado a la cantidad de combustible

necesario para llevar la caldera a la temperatura y presión adecuados para
operar la turbina; este costo depende del tiempo que ha estado apagada la
planta como se muestra en la figura:
C($)
C0
Figura 8.1: Costo de partida de una planta térmica
Esta relación muestra como el mayor costo posible se produce cuando

la planta arranca en frı́o mientras que disminuye a medida que el tiempo
en que está apagada disminuye. Esto puede ser expresado por la siguiente
relación:
C = C0 · (1 − e−φ·t ) (8.1)
en donde:
C0 : Costo de partida en frio
phi : Razón de enfriamiento de la unidad.
Costo de apagado: Este costo es asociado principalmente a la

operación y mantenimiento de la planta por lo cual es independiente del
tiempo y toma un valor fijo.
Costo del combustible: Como se describió anteriormente, el costo

de operación de una unidad térmica depende en gran medida del costo
del combustible, este costo esta asociado a la compra de combustible que
normalmente se realiza por contratos de mediano o largo plazo para evitar
que las fluctuaciones del mercado de combustible impacte al mercado
eléctrico.
Potencia[BTU]
Válvula 3
Válvula 1 Válvula 2
Potencia[MW]
Figura 8.2: Relación entre la potencia de entrada y salida de la planta

térmica
En la figura ?? se muestra la relación entre la potencia de entrada y

salida de una central térmica, este relación permite determinar la cantidad
de combustible necesario para generar determinada cantidad de potencia
activa. Teniendo en cuenta el costo del combustible, es posible llevar ésta
relación a un función de costo vs potencia generada la cual es utilizada
en los diferentes análisis económicos. La curva presenta una serie de
discontinuidades que representan los puntos de apertura o cierre de
las diferentes válvulas de admisión. No obstante, es posible hacer una
aproximación a una forma continua representada por un polinomio.
Normalmente las funciones de costo son representadas en formas

cuadráticas o lineales por tramos; la aproximación más comúnmente
utilizada es cuadrática:
127
ak
fk (Pk ) = · Pk2 + bk · Pk + ck (8.2)
2
en donde ak , bk son constantes asociadas a los costos variables (costos

de generación) de planta k; ck representa los costos fijos de la misma.
Como se dijo anteriormente, los costos de generación de la planta solo
dependen de la potencia generada por ella misma, por esta razón se dice
que el problema de despacho económico es independiente en el tiempo
(estático).
8.2 Despacho económico de plantas térmicas
El problema de despacho económico consiste en determinar la

cantidad de potencia generada por parte de cada una de las unidades
térmicas. Una primera aproximación consiste en eliminar las restricciones
de la red, esto implica que cada una de las unidades de generación
ası́ como la demanda están conectadas al mismo nodo. Un modelo ası́
planteado puede ser útil en el caso de tener unidades de generación
con caracterı́sticas distintas en una misma planta, igualmente como
aproximación inicial al problema en sistemas interconectados.
Un sistema con NT unidades térmicas conectadas como se muestra

en la figura 8.3 alimenta una demanda de potencia PD , las centrales
presentan funciones de costo como se describe en la ecuación 8.2. El
problema consiste en determinar la potencia generada Pk por cada una
de las unidades de tal forma que el costo operativo sea mı́nimo y que se
satisfaga la demanda.
El problema matemático puede ser planteado como se muestra a

continuación:
Función objetivo:
NT
M in Costo = fk (Pk ) (8.3)
k=1
Unidad 1 Unidad 2 Unidad NT
...........
PD
Figura 8.3: Esquema de despacho uninodal
Restricción de balance de potencia
NT
P k = PD (8.4)
k=1
La solución a este problema de programación no lineal con una

restricción de igualdad se obtiene mediante la teorı́a de multiplicadores
de Lagrange:
∂
λ= fk (Pk ) (8.5)
∂Pk
NT
Pk = P D (8.6)
k=1
Puesto que las funciones de costo son normalmente cuadráticas, el

sistema resultante es lineal por lo cual la solución es única.
λ k = ak · P k + b k (8.7)
En conclusión, para garantizar un despacho económico, es necesario

que el valor de λ sea el mismo en todas las unidades de generación y
129
que se satisfaga la demanda. El multiplicador λ es denominado costo

incremental del combustible y representa la razón de cambio del costo frente
a un aumento o disminución de una unidad de potencia del sistema.
Si un conjunto de generadores presenta costos incrementales diferentes,
es posible disminuir una unidad de carga del generador con mayor λ
y entregarla a la de menor valor, el ahorro con esta operación estará
dado por la diferencia de los costos incrementales; al nuevo punto de
operación, corresponde un nuevo conjunto de multiplicadores de lagrange
más cercanos entre sı́, de tal forma que se puede repetir la operación hasta
el punto en que todas las plantas presenten el mismo costo incremental,
en este punto, cualquier movimiento producirá un sobre-costo.
Ejemplo 8.1
Se tienen dos centrales generadoras con los parámetros mostrados en la tabla

8. Se desea alimentar una carga de 90M W .
Tabla 8.1: Datos de generación.

Nodo a b c
1 108 × 10−4 11,669 5
2 176 × 10−4 10,333 5
Para dar solución a este problema se deben cumplir dos condiciones:

igualdad en los costos incrementales de combustible y balance de potencia
generada y demandada.
λ = 108 × 10−4 · PG1 + 11, 669 = 176 × 10−4 · PG2 + 10, 333 (8.8)
PG1 + PG2 = 90 (8.9)
este sistema de ecuaciones tiene solución única:

PG1 = 8, 73M W (8.10)
PG2 = 81, 27M W (8.11)
El costo incremental de las dos plantas es el siguiente
$
λ = 11, 763284 (8.12)
MWh
Mientras que el costo total de la operación es de $1010, 1672 Los costos

fijos de la planta c no están directamente involucrados en el proceso
de optimización; como puede esperarse estos costos siempre serán los
mismos (aún cuando la planta no genere).
Si la demanda es modificada a un valor de 91M W , la operación

más económica (resultado del proceso de optimización) tiene un costo
de $1021, 9305 lo cual corresponde a un aumento de $11, 76329 (costo
incremental).
Ejemplo 8.2
Dos centrales generadoras deben abastecer una demanda variable para un

periodo de 24 horas. La demanda se muestra en la figura 8.4. Las funciones
de costo de cada una de las centrales son mostradas a continuación:
f1 = 3, 6 × 10−6 · P13 + 9, 8 · P1 + 8 (8.13)
f2 = 11, 2 · P2 + 10 (8.14)
131
600 MW
300 MW
24 h
Figura 8.4: Curva de duración de carga
El problema en este caso consiste en minimizar el costo total de

la energı́a despachada a lo largo de un horizonte de tiempo. Las
condiciones de optimalidad son las mismas ya que el problema del
despacho económico es independiente del tiempo. No obstante, en este
caso los costos incrementales resultantes no son de tipo lineal.
La demanda puede ser representada por la expresión
300
PD (t) = 600 − ·t (8.15)
24
La depedendencia de la demanda con el tiempo implica que tanto los

costos incrementales como las potencias generadas son tambien funciones
del tiempo. Aplicando los criterios de optimalidad se tiene:
3 · (3, 6 × 10−6 ) · (P1 (t))2 + 9, 8 = 11, 2 (8.16)
de donde se obtiene:

360 t ≤ 19, 2
P1 (t) = (8.17)
600 − 12, 5 · t t ≥ 19, 2

240 − 12, 5 · t t ≤ 19, 2
P2 (t) = (8.18)
0 t ≥ 19, 2
La razón de la discontinuidad en la respuesta, es la generación negativa

asociada a la planta 2 en caso de ser despachada con el criterio de
los multiplicadores de lagrange. Desde luego, una generación negativa
en este generador (lo cual significa la máquina sı́ncrona trabajando
como motor) no implica un ahorro en términos del costo, al contrario,
presentarı́a gastos constantes. El punto en el cual se presenta este
funcionamiento anormal de la máquina es en t = 19, 2 por lo cual, de
allı́ en adelante, la generación debe ser despachada por la primera central.
El costo incremental del combustible del sistema es el siguiente:

11, 2 t ≤ 19, 2
λ(t) = −6 2 (8.19)
10, 8 × 10 · (600 − 12, 5 · t) + 9, 8 t ≥ 19, 2
En la siguiente sección se mostrara de forma más detallada las

implicaciones de las restricciones sobre los costos incrementales.
8.2.25
Las centrales térmicas presentan restricciones con respecto a la
Restricciones.
capacidad de generación, las cuales deben ser modeladas en el despacho
económico. Un lı́mite de potencia máxima (Pmax ) esta dado por la
capacidad de la fuente primaria de energı́a. Ası́ mismo, existe un limite
inferior Pmin por debajo del cual la planta no está en capacidad de generar
de forma confiable. El modelo matemático resultante presenta la siguiente
forma:
Función objetivo
NT
M in Costo = fk (Pk ) (8.20)
k=1
Restricción de balance de potencia

133
NT
Pk = P D (8.21)
k=1
Restricciones de capacidad de generación
Pk(min) ≤ Pk ≤ Pk(max) (8.22)
La solución de este tipo de problemas exige el uso de herramientas

más sofisticadas de programación no lineal (ver anexo B), no obstante una
primera aproximación consiste en obviar las inecuaciones y considerar el
sistema con restricciones de igualdad; si el resultado no viola ninguna
restricción de capacidad el problema ha sido solucionado. En caso
contrario, se busca la restricción violada y se obliga a estar en el margen
correspondiente, de esta forma, la variable toma un valor fijo por lo cual
deja de participar en el problema de optimización, posteriormente se da
solución nuevamente al problema resultante, el proceso se repite de forma
anidada hasta que no se presente ninguna violación en las restricciones de
desigualdad.
Como ya se ha mencionado, los costos incrementales del combustible

representan el valor en que incurre el sistema al aumentar en 1M W la
demanda del mismo; bajo esta definición, cuando una restricción limita la
generación de una unidad, el costo incremental del sistema corresponde a
los costos incrementales de las centrales restantes como se muestra en la
figura 8.6.
Como ya se ha mencionado, los costos incrementales del combustible

representan el valor en que incurre el sistema al aumentar en 1M W la
demanda del mismo; bajo esta definición, cuando una restricción limita la
generación de una unidad, el costo incremental del sistema corresponde a
los costos incrementales de las centrales restantes como se muestra en la
figura 8.6.
Costos Incrementales
Costo incremental
del sistema Pmax
T1 T2 T3
Potencia
Figura 8.5: Costos incrementales de 3 plantas termicas con la tercera planta

en el lı́mite superior
Ejemplo 8.3
Cuatro unidades operan en conjunto para cubrir una demanda total de

800M W , encontrar la salida requerida en cada planta para que el despacho sea
óptimo teniendo en cuenta que la generación máxima en cada una de las plantas
es de 400M W y la generación mı́nima es de 10M W . Los costos incrementales de
combustible en M$W para cuatro unidades son los siguientes:
∂f1
λ1 = = 0, 0120 · P1 + 9 (8.23)
∂P1
∂f2
λ2 = = 0, 0096 · P2 + 6 (8.24)
∂P2
∂f3
λ3 = = 0, 0080 · P3 + 8 (8.25)
∂P3
135
Costos Incrementales
Costo incremental
del sistema Pmax
T1 T2 T3
Potencia
Figura 8.6: Costos incrementales de 3 plantas térmicas con la tercera planta

en el lı́mite superior
∂f4
λ4 = = 0, 0068 · P4 + 10 (8.26)
∂P4
Inicialmente, se plantea el sistema sin considerar las restricciones de

capacidad de cada planta, el sistema de 4 ecuaciones resultante es el
siguiente:
⎛ ⎞ ⎛ ⎞ ⎛ ⎞
0.0120 −0.0096 0.0000 0.0000 P1 −3
⎜ 0.0120 0.0000 −0.0080 0.0000 ⎟ ⎜ P2 ⎟ ⎜ −1 ⎟
⎜ ⎟·⎜ ⎟=⎜ ⎟
⎝ 0.0120 0.0000 0.0000 −0.0068 ⎠ ⎝ P3 ⎠ ⎝ 1 ⎠ (8.27)
1.0000 1.0000 1.0000 1.0000 P4 800
la solución encontrada es:

⎛ ⎞
92.4
⎜ 428.0 ⎟
P =⎜ ⎟
⎝ 263.6 ⎠ (8.28)
16.0
como se puede observar, P2 supera la máxima capacidad de la dicha

unidad, por tal razón se fija el valor en 400 y se repite el proceso pero esta
vez con las tres variables restantes:
⎛ ⎞ ⎛ ⎞ ⎛ ⎞
0.0120 −0.0080 0.0000 P1 8−9
⎝ 0.0120 0.0000 −0.0068 ⎠ · ⎝ P3 ⎠ = ⎝ 10 − 9 ⎠ (8.29)
1.0000 1.0000 1.0000 P4 800 − 400
la nueva solución es:
⎛ ⎞
98.96550
⎜ 400.0000 ⎟
P =⎜ ⎟
⎝ 273.4483 ⎠ (8.30)
27.5862
esta solución cumple con todas las restricciones y presenta un costo

total de $6882, 586 (asumiendo costos fijos de cero)1
1
Este problema puede ser solucionado utilizando la herramienta Gams como se
muestra en el anexo B
137
8.2.26 Red de
La red de transmisión también puede ser considerada como restricción transmisión
adicional, una forma simple de lograrlo es mediante un flujo DC.
P km
k
m
Figura 8.7: flujo de potencia en una lı́nea de transmisión
El flujo de potencia activa en una lı́nea del sistema que conecta los
nodos k y m esta dado por la ecuación 8.31
Pkm = gkm · Vk2 − bkm · Vk · Vm · sen (θkm ) − gkm · Vk · Vm · cos (θkm ) (8.31)
en donde
rkm
gkm = 2
(8.32)
rkm + x2km
−xkm
bkm = 2
(8.33)
rkm + x2km
En sistemas de potencia es posible considerar las siguientes

aproximaciones:
−1
bkm ≈ (8.34)
xkm
Vk = Vm ≈ 1 (8.35)
sen (θkm ) ≈ θkm (8.36)
cos (θkm ) ≈ 1 (8.37)
Bajo estas consideraciones la ecuación de flujo de potencia puede ser

re-escrita como:
1
Pkm = · (θk − θm) (8.38)
xkm
Esta ecuación lineal puede ser llevada en cuenta en el despacho

económico adicionando una restricción de igualdad por cada lı́nea
(ec.8.38). Ası́ mismo, la capacidad de cada lı́nea puede ser modelada
como:
|Pkm | < Pkm(max) (8.39)
Las restricciones de valor absoluto pueden ser fácilmente

transformadas en inecuaciones lineales:
− Pkm(max) < Pkm < Pkm(max) (8.40)
Las técnicas manuales antes mencionadas pueden ser utilizadas para

dar solución a este modelo, sin embargo, el uso del computador facilita
los cálculos en esta materia;
Ejemplo 8.4
El sistema mostrado en la figura 8.8 presenta dos plantas hidráulicas y dos

plantas térmicas. La generación de las plantas hidráulicas ha sido fijada en
139
30M W en cada una.2 Los datos de las plantas térmicas están dados en la tabla 8,
la base fijada es de 100M W . Cada unidad de generación térmica puede generar
entre 0 y 80M W Se desea encontrar un despacho económico considerando la red,
los datos de la misma están dados en la tabla 8
Térmica
Térmica
1
2
1 40 MW 2
60 MW
5
3 4
100 MW
Hidráulica Hidráulica
1 2
Figura 8.8: Sistema de 5 nodos
Tabla 8.2: Parámetros de la red (figura 8.8)

N1 N2 X [pu] Pmax[pu]
1 5 0.15 0.40
2 5 0.1 0.18
3 5 0.13 0.40
4 5 0.12 0.40
El modelo utilizando Gams se muestra a continuación:
2
El despacho de un sistema hidrotérmico será tratado en el capı́tulo siguiente.

$title Despacho economico de p l a n t a s t e r m i c a s
$ontext
Despacho de p l a n t a s t e r m i c a s
te ni en do en cuenta l a s r e s t r i c c i o n e s de l a red
$offtext
s e t s pt P l a n t a s t e r m i c a s / 1 * 2 /
n nodos d e l s i s t e m a / 1 * 5 /
l l i n e a s de t r a n s m i s i o n / 1 * 4 / ;
t a b l e DT( pt , * ) Datos de l a s p l a n t a s t e r m i c a s
a2 b Max Min Nodo
1 108 1166.9 0.8 0 1
2 176 1033.3 0.8 0 2;
parameter Carga ( n ) c a r g a e l s i s t e m a / 1 0.4 ,

2 0.6 ,
3 −0.3 ,
4 −0.3 ,
5 1.0/;
table Sistema ( l , * ) Datos d e l s i s t e m a de t r a n s m i s i o n
N1 N2 X Pmax
1 1 5 0.15 0.40
2 2 5 0.1 0.18
3 3 5 0.13 0.40
4 4 5 0.12 0.40;
variables
P ( pt ) P o t e n c i a generada por cada una de l a s p l a n t a s t e r m i c a s
pf ( l ) F l u j o de P o t e n c i a por l a s l i n e a s
Ang ( n ) Angulos nodales
Coste Costo de l a p o t e n c i a a c t i v a ;
P . l o ( pt ) = DT( pt , ’ Min ’ ) ;
P . up ( pt ) = DT( pt , ’ Max ’ ) ;
pf . l o ( l ) = − Sistema ( l , ’ Pmax ’ ) ;
pf . up ( l ) = Sistema ( l , ’ Pmax ’ ) ;
equation
fo funcion o b j e t i v o
rf ( l ) r e s t r i c c i o n de f l u j o de p o t e n c i a en l a s l i n e a s
pn ( n ) p o t e n c i a nodal
slack angulo en e l nodo s l a c k ;
fo . . Coste = e = sum ( pt , DT( pt , ’ a2 ’ ) * P ( pt ) * P ( pt ) / 2 + DT( pt , ’ b ’ ) * P ( pt ) ) ;
rf ( l ) . . pf ( l ) = e = ( sum ( n , ( Ang ( n ) $ ( ord ( n)= Sistema ( l , ’ N1 ’ ) ) ) )
− sum ( n , ( Ang ( n ) $ ( ord ( n)= Sistema ( l , ’ N2 ’ ) ) ) ) )
/Sistema ( l , ’ X ’ ) ;
pn ( n ) . . sum ( l , pf ( l ) $ ( ord ( n)= Sistema ( l , ’ N1’ ) ) ) −
sum ( l , pf ( l ) $ ( ord ( n)= Sistema ( l , ’ N2 ’ ) ) ) = e=
sum ( pt ,−P ( pt ) $ ( ord ( n)=DT( pt , ’ Nodo ’ ) ) ) + Carga ( n ) ;
s l a c k . . Ang ( ’ 1 ’ ) = e = 0 ;
model DespachoEconomico / a l l / ;
solve DespachoEconomico minimizing Coste using nlp ;
display Ang . l , P . l

141
La solución teniendo en cuenta la red es de 62M W en la planta

térmica 1 y de 78M W en la planta térmica 2. La solución sin tener en
cuenta la red es de 60M W y 80M W respectivamente. Para llegar a esta
solución se pueden realizar una secuencia de cálculos simples añadiendo
progresivamente las restricciones:
• Inicialmente se realiza el despacho sin tener en cuenta ninguna

restricción de desigualdad; en este caso la solución es P1 =
39, 72M W y P2 = 100.28M W . Los costos incrementales del
combustible son iguales (12.098$/M W (Los costos incrementales
no deben ser expresados en por unidad para una adecuada
interpretación).
• La planta térmica 2 viola su restricción de capacidad por lo cual se

debe fijar en este valor, la planta térmica 1 aumenta su generación a
60M W . Los costos incrementales están dados por la planta 1 la cual
esta en libertad de generar potencia adicional. λ = 12, 32$/M W
• Finalmente al ubicar la generación calculada en el sistema se

producen violaciones en la cargabilidad de las lı́neas, éstas implican
una nueva modificación en la generación a los valores arrojados
por GAMS de 62M W y 78M W con un costo incremental de λ =
12.33$/M W .
8.2.27 Pérdidas
Las pérdidas pueden ser incluidas al modelo utilizando un flujo de de la red
carga AC como restricción adicional. Sin embargo, existe una forma inicial (Modelo DC)
de considerarlas de forma parcial. Bajo el esquema de nodo único, la
principal restricción corresponde al balance de potencia; en este caso,
dicha ecuación debe ser modificada con las pérdidas de potencia activa:
NT
Pk = P D + PL (8.41)
k=1
en donde PL corresponde a las pérdidas totales del sistema. Bajo esta

consideración, las condiciones de optimalidad son las siguientes:

∂fk 1
λ= · (8.42)
∂Pk 1 − ∂PL
∂Pk
Nuevamente la condición básica corresponde al cumplimiento de

igualdad entre de los costos incrementales y el balance de potencia.
Aparece un nuevo término denominado factor de penalización por pérdidas
(ξk ) el cual modifica los costos incrementales del combustible, dicho factor
es definido como:
1
ξk = (8.43)
1 − ∂PL
∂Pk
El factor de penalización puede ser obtenido mediante las restricciones

de flujo de carga, no obstante, el uso de coeficientes de costo (B) permiten
un desarrollo aproximado como se describe a continuación:
El flujo de potencia dado por la ecuación 8.31 puede ser simplificado

teniendo en cuenta una expansión hasta el segundo término de la función
cos(θkm ) y voltajes 1pu:
2
θkm
cos(θkm ) ≈ 1 − (8.44)
2

θ2
Pkm ≈ gkm − gkm · 1 − km − bkm · θkm (8.45)
2
Las pérdidas de potencia activa en la lı́nea km están dadas por la suma

algebraica entre el flujo de potencia Pkm y Pmk :
Pperd(km) = Pkm + Pmk (8.46)

143
2
Pperd(km) = gkm · θkm (8.47)
Pperd(km) = gkm · (xkm · Pkm )2 (8.48)
Fácilmente se puede concluir que las pérdidas totales son función de

la potencia nodal pues ésta determina los ángulos en cada nodo (θkm =
f (PN )).
2
PL = gkm · θkm (8.49)
Una forma mas detallada de determinar las pérdidas como una función
de la potencia generada consiste en calcular los coeficientes de pérdidas
utilizando un flujo de carga AC. En este caso es necesario aplicar una
transformación invariante en potencia. Los coeficientes de pérdidas varı́an
con el nivel de tensión por lo cual deben ser recalculados en cada iteración.
Ejemplo 8.5
Para el ejemplo anterior (figura 8.8) asumir la resistencia equivalente igual

1
a x .
10 km
Las pérdidas como función de la potencia generada por las dos plantas
térmicas se expresa por la siguiente ecuación:
3 1
PL = · (PT 1 − 0, 4)2 + · (PT 2 − 0, 6)2 (8.50)
202 101
Para un despacho a nodo único (tabla 8 y demanda de 120M W )

los factores de penalización por pérdidas son respectivamente 0.99992 y
1.00804, con estos valores, los costos incrementales del combustible son
12, 0971 y 12, 1953, por esta razón es necesario una generación capaz de
asegurar que éstos costos incrementales sean iguales, se requiere entonces

dar solución al siguiente sistema de dos ecuaciones con dos incógnitas3 .
108 · PT 1 + 1166, 9 176 · PT 2 + 1033.3

6 = 2 (8.51)
1 − 202 · (PT 1 − 0, 4) 1 − 101 · (PT 2 − 0, 6)
PT 1 + PT 2 = 1, 4 + PL (8.52)
La solución a este sistema de ecuaciones es PT 1 = 0, 432775 y PT 2 =

0, 979506
8.2.28 Pérdidas Es posible considerar las pérdidas de la red usando un modelo más
de la red exacto, para esto es necesario hacer uso de los conceptos de análisis de
(Modelo AC) sistemas de potencia.
Para un punto de operación determinado, dado por las variables de

estado que en el caso de los sistemas de potencia son los voltajes nodales,
la potencia nodal puede ser calculada como:
Pk = PTk − PDk (8.53)
en donde PTk es la potencia generada por la planta térmica conectada

al nodo k mientras que PDk es la potencia demandada en dicho nodo, la
cual es un dato gracias a la función de predicción de la demanda.
Desde el punto de vista del sistema la potencia nodal puede ser

expresada en términos del voltaje y la corriente respectiva Sk = Vk · I∗k :
Pk + j · Qk = Vk · (cos(θk ) + j · sen(θk )) · I∗k (8.54)

3
En un sistema con una mayor cantidad de generadores este sistema de ecuaciones
debe ser solucionado usando un método numérico
145
en donde I∗k es la corriente nodal inyectada al sistema expresada en

forma compleja. De esta relación se obtiene

cos(θk ) sen(θk ) sen(θk ) cos(θk )
Ik = Pk · + Qk · + j · Pk · − Qk ·
Vk Vk Vk Vk
(8.55)
A su vez, las pérdidas de potencia del sistema son expresadas como

la suma de las potencias nodales, las cuales se determinan a partir de la
matriz de impedancia nodal ZBU S como:
SL = [I] · [ZBU S ] · [I]∗ (8.56)
La simetrı́a de la matriz ZBU S asegura la siguiente relación:
PL = [Real(I)] · [RBU S ] · [Real(I)] + [Imag(I)] · [RBU S ] · [Imag(I)] (8.57)
usando el resultado obtenido en la ecuación 8.55 se obtiene:
N N
cos(θk ) sen(θk ) cos(θm ) sen(θm )
PL = rkm · Pk · + Qk · · Pm · + Qm ·
k=1 m=1
Vk Vk Vm Vm

sen(θk ) cos(θk ) sen(θm ) cos(θm )
+ Pk · − Qk · · Pm · − Qm · (8.58)
Vk Vk Vm Vm
en donde rkm es un término general de la matriz RBU S mientras que

Pk es la potencia nodal (PGk − PDk ). Desde luego, esta ecuación incluye
los nodos de carga, los cuales se cancelan en el momento de derivar con
respecto a cada una de las potencias generadas, igualmente, la derivada
de las pérdidas con respecto a cada una de la potencias generadas es
equivalente a derivar con respecto a la potencia nodal correspondiente.

Finalmente, usando algunas identidades trigonométricas y la simetrı́a
caracteristica del problema, es posible mostrar la expresión en una forma
más simplificada:
N N N
N

rkm · cos(θkm ) rkm · sen(θkm ) · Qm
PL = · Pk · P m − 2 · · Pk
k=1 m=1
Vk · Vm k=1 m=1
Vk · Vm
N N
rkm · cos(θkm )
+ · Qk · Qm (8.59)
k=1 m=1
V k · Vm
se definen las siguientes constantes:
rkm · cos(θkm )
βkm = (8.60)
Vk · Vm
N
rkm · sen(θkm ) · Qm
βk0 = (8.61)
m=1
Vk · Vm
N N
rkm · cos(θkm )
β00 = · Qk · Qm (8.62)
m=1
V k · Vm
k=1
por lo que las perdidas pueden ser expresadas como:
N N N
PL = β00 − 2 · βk0 · Pk + βkm · Pk · Pm (8.63)
k=1 k=1 m=1
esta ecuación permite calcular el factor de penalización por pérdidas

de una forma más simple y exacta; no obstante, a cada modificación del
despacho, las constantes γ cambian igualmente, por lo que es necesario
re-calcular en flujo de carga en un algoritmo iterativo:
147
1. Inicialmente se calcula el despacho sin tener en cuenta las pérdidas

del sistema, de esta forma se obtiene un conjunto de potencias a
generar por parte de cada uno de los nodos P V .
2. A partir de los resultados del despacho se calcula un flujo de

carga por el método de Newton-Raphson4 ; con los nuevos datos de
tensiones y potencias reactivas es posible actualizar los valores de
beta
3. Se calcula el despacho teniendo en cuenta las pérdidas, se asumen los

factores beta constantes, por lo que se trata de un sistema no lineal
de ecuaciones (tal y como se trató en el modelo DC)
4. Si los costos incrementales del combustible sufren una variación por

encima de la tolerancia, se regresa al paso 2.
Ejemplo 8.6
Dar solución al ejemplo anterior teniendo en cuenta el modelo AC, asumir que
todas las cargas tienen una demanda de potencia reactiva de 10M V Ar y que las
tensiones en los nodos generadores son de 1pu. Los efectos capacitivos de las lineas
de transmisión son de Y /2 = 0, 0005
Los resultados del flujo de carga (partiendo de un despacho inicial sin

pérdidas) después de 5 iteraciones se muestran en la tabla 8
Tabla 8.3: resultados del flujo de carga

nodo V [pu] A[] Q[pu]
1 1,0000000 0,0000000 0,0403680
2 1,0000000 -0,5919422 0,0684400
3 0,9913941 -1,7081456 0,1000000
4
Ver referencia [1]
Las constantes de pérdidas β son las siguientes
⎛ ⎞
0, 0090636 −0, 0046879 −0, 0022059
βkm = ⎝ −0, 0046879 0, 0065632 −0, 0009455 ⎠ (8.64)
−0, 0022059 −0, 0009455 0, 0015902
⎛ ⎞
0, 32636588
βk = ⎝ 0, 37974879 · 10−5 ⎠ (8.65)
0, 39164764
β00 = 6, 626276315643167 · 10−5 (8.66)
La ecuación de pérdidas toma la siguiente forma
PL = 0, 734 · 10−3 + 0, 906 · 10−2 · PT21 + 0, 133 · 10−3 · PT 1

− 0, 938 · 10−2 · PT 1 · PT 2 − 0, 338 · 10−2 · PT 2 + 0, 656 · 10−2 · PT22 (8.67)
para la condición de operación descrita las pérdidas de potencia activa

son de 0, 0010791pu.
El sistema de ecuaciones resultante es el siguiente:
1080 · PT 1 + 11669
λ= (8.68)
0, 999867 − 0, 01812 · PT 1 + 0, 938 · 10−2 · PT 2
1760 · PT 2 + 10333
λ= (8.69)
1, 00338 + 0, 938 · 10−2 · PT 1 − 0, 01312 · PT 2
PT 1 + PT 2 = 1, 4 + PL (8.70)
La solución a este sistema de ecuaciones da como resultado un nuevo

despacho económico el cual modifica el punto de operación de la red, por
149
lo tanto, es necesario re-calcular las variables de estado (voltajes) y los

factores beta. En este caso la potencia generada por las plantas es de 0, 4276
y 0, 9758 respectivamente.
Los resultados obtenidos por el modelo AC son en la mayorı́a de los

casos muy similares a los obtenidos por medio del modelo DC razón por la
cual en muchos casos se prefiere el primero. Se puede comprobar como en
las siguientes iteraciones, los cambios en los factores β son despreciables
por lo cual es posible mantenerlos constantes.
Ejemplo 8.7
Un sistema de generación está compuesto por dos unidades térmicas. La

generación óptima de la planta 1 es 370M W . Calcular el costo incremental del
combustible para cada planta (teniendo en cuenta el factor de penalización por
pérdidas) y la potencia generada por la planta 2, si los costos del combustible de
dos unidades y las pérdidas están dadas por las siguiente ecuaciones:
F1 = 7, 74 · P1 + 0, 00107 · P12 (8.71)
F2 = 7, 72 · P2 + 0, 00072 · P22 (8.72)
PL = 0, 5 × 10−3 · P12 + 0, 2 × 10−3 · P22 (8.73)
Las condiciones de optimalidad son las siguientes:
7, 74 + 2 · 0, 00107 · P1 7, 72 + 2 · 0, 00072 · P2
= (8.74)
1 − 2 · 0, 5 × 10 · P1
−3 1 − 2 · 0, 2 × 10−3 · P2
de donde se obtiene directamente P2 = 849, 136M W . Usando la

ecuación de pérdidas PL se puede obtener la demanda real de potencia
y la información restante del problema:
• Pérdidas 212, 65M W (Igual al 17% de la potencia generada)
• Demanda de potencia 1006, 486M W
• Generación de potencia: 1219, 136M W
• Costo de la operación de la planta 1: 3010, 283
• Costo de operación de la planta 2: 7074, 47
• Costo total de la operación: 10084, 753
Un modelo que no considere las pérdidas llegarı́a al siguiente resultado
• Demanda de potencia 1006, 486M W
• Potencia generada por la unidad 1: 399, 256M W
• Potencia generada por la unidad 2: 607, 230M W
No obstante, en el momento de poner este despacho en el sistema,

el nodo slack (en este caso el nodo 2) tendrı́a que asumir una potencia
adicional de 153M W dada por la ecuación de pérdidas.
8.3 Arranque y parada de unidades térmicas.
El problema del despacho de plantas térmicas toma una caracterı́stica

dinámica debido a los costos de encendido y apagado asociado a la
operación de las mismas. Este problema denominado Unit commitment,
puede ser solucionado mediante la técnica de programación dinámica.
Como se mostró en secciones anteriores, el costo de arranque una planta
térmica depende del tiempo de apagado de la misma, no obstante, la
caracterı́stica asintótica de la función exponencial permite considerar un
costo constante igual al costo de arranque en frı́o cuando se discretiza el
151
tiempo en periodos extensos. Los costos de apagado pueden igualmente

ser considerados constantes.
La metodologı́a de programación dinámica requiere un árbol de

decisión sobre el cual aplicar el algoritmo de optimización. En dicho árbol,
cada fila corresponde a una combinación de unidades térmicas mientras
que cada columna representa un periodo de tiempo. Las conexiones entre
los nodos representan los costos de arranque o parada mientras que el
valor en cada nodo se calcula utilizando la técnica de los multiplicadores
de Lagrange en un despacho económico con las plantas activas en el
estado correspondiente. Existen condiciones infactibles en el gráfico, por
ejemplo, salidas por mantenimiento, despacho por debajo del lı́mite de
alguna planta activa, minimo tiempo de apagado o máximo tiempo de
encendido. En estos casos se debe eliminar la rama correspondiente o
asociarle un costo infinito.
El problema de enganche de plantas térmicas puede ser solucionado

igualmente por técnicas heurı́sticas o algoritmos combinatoriales tales
como las técnicas evolutivas, los algoritmos genéticos entre otros.
Ejemplo 8.8
Determinar el cronograma de arranque y parada para las unidades térmicas

descritas en la tabla 8 y la curva de carga mostrada en la figura 8.9.
Para las dos plantas térmicas se asume un valor de φ = 4h−1 y de Co =

400 en la ecuación ??. Esto significa que el costo de arranque es constante
(Co ) pues en cada periodo de tiempo la función ha alcanzado el estado
estacionario. Las dos plantas están en operación en el instante inicial. La
potencia mı́nima generada es de 20M W y la máxima es de 150M W . La
planta T 2 debe entrar en mantenimiento en la hora 22 por lo cual el último
estado es T 1. En la tabla 8.4 se muestra el despacho de potencia para cada
nivel de demanda y los tres posibles estados de operación (T 1,T 2,T 1−T 2).
120
Demanda [MW]
110
90
80 80
0 4 12 18 22 24
Tiempo [h]
Figura 8.9: Curva de carga discretizada
Tabla 8.4: Diferentes condiciones de despacho

Demanda T1 T2 T1 y T2
[MW] PA PB Costo PA PB Costo PA PB Costo
80 80 0 978 0 80 893 2,535 77,7648 893
90 90 0 1104 0 90 1011 8,732 81,2676 1010
110 110 0 1359 0 110 1253 21,127 88,8730 1247
120 120 0 1488 0 120 1377 27,324 92,6760 1366
Los periodos de demanda de 80M W y 90M W deben ser generados de

forma independiente por una de las dos plantas pues el despacho conjunto
implica una generación inferior a 20M W en la planta T 1. El gráfico 8.10
muestra el problema para ser resuelto mediante programación dinámica.
Los costos de parada de cada planta se asumen despreciables.
Los costos de cada nodo resultan de multiplicar la potencia despachada

en cada estado por el tiempo en horas del mismo. La solución se muestra
en la tabla 8.5. La ruta óptima (ACF HL) tiene un costo de $27482
T1 153
3912 10872 6624 5952 1956
B D I L
G
T1-T2
0 A 9976 5464
E J
T2
3572 10024 6066 5508

C F H K
Figura 8.10: Diagrama para solucionar el problema de unit commintment
8.4 Generación de un equivalente térmico.
Un sistema con múltiples plantas térmicas puede ser llevado a un

sistema equivalente utilizando un estrategia de optimización. Para esto
se realiza el proceso de despacho económico sin red utilizando todos los
generadores presentes en el sistema desde un lı́mite inferior igual a la
suma de los lı́mites inferiores de cada una de las centrales y un lı́mite
superior igual a la suma de los lı́mites superiores del sistema.
Este procedimiento se realiza en pasos discretos hasta obtener una

curva que relaciona la potencia demandada y los costos totales. Esta
curva puede ser llevada a una forma cuadrática por medio de un
ajuste de mı́nimos cuadrados. Los parámetros de ajuste de la térmica
corresponden a la planta térmica equivalente. No siempre la curva
resultante corresponde a una forma cuadrática, razón por la cual el
equivalente no siempre es una buena aproximación, ası́ mismo, la
discretización de la demanda puede influir en la exactitud del equivalente.
Las pérdidas del sistema pueden ser incluidas en el equivalente.
Es posible desarrollar un equivalente térmico resolviendo de forma

Tabla 8.5: Rutas óptimas para el ejemplo propuesto

Nodo Costo Origen
B 3912 A
C 3572 A
D 14784 B
E 13948 C
F 13596 C
G 20572 E
H 19662 F
I 26014 H
J 25526 H
K 25170 H
L 27482 J
analı́tica el sistema de ecuaciones resultante del proceso de optimización.

Este tipo de solución es practica en el caso de tener pocas centrales
térmicas, pero puede ser mas complejo en un sistema real. La solución
analı́tica es siempre preferible por su exactitud.
Las condiciones de optimalidad en un sistema de NT plantas térmicas

son las siguientes:
λ = ak · P k + b k (8.75)
NT
Pk = P D (8.76)
k=1
de donde se obtiene:
Np
1 1 bk
λ= N p 1
· PD + N p 1
· (8.77)
k=1 ak k=1 ak k=1
ak
definiendo las constantes equivalente aD y bD se tiene:

155

1
aD = N p 1
(8.78)
k=1 ak
Np
bk
b D = aD · (8.79)
k=1
ak
finalmente, la potencia a generar por parte de cada unidad térmica es

la siguiente:
λ − bk
Pk = (8.80)
ak
El uso de una planta térmica equivalente puede facilitar

procedimientos más complejos como el caso del despecho hidrotérmico
como se mostrará en el siguiente capitulo. El tiempo computacional
en el momento de realizar el equivalente se ve compensado con la
simplificación que se obtiene al llevar un sistema de múltiples plantas
térmicas a un sistema equivalente con una sola central.
Ejemplo 8.9
Generar un equivalente térmico para el sistema del ejemplo N 1
Se trata de un sistema con solo dos unidades de generación por lo que

la solución analı́tica es practica, para esto, se resuelve usa la ecuacion 8.77.
λ = 0, 0067 · PD + 11, 1609 (8.81)
las potencias generadas son:

P1 = 0, 6197 · PD − 47, 0423 (8.82)
P2 = 0, 3803 · PD + 47, 0423 (8.83)
Sin embargo, las restricciones de capacidad de la central modifican los

resultados; por ejemplo, por debajo de P1 = 75, 9091M W la generación
se hace negativa, por tal razón, la totalidad de la generación debe ser
abastecida por la central 2 :

0 PD ≤ 75, 9091
P1 = (8.84)
0, 6197 · PD − 47, 0423 PD ≥ 75, 9091

PD PD ≤ 75, 9091
P2 = (8.85)
0, 3803 · PD + 47, 0423 PD ≥ 75, 9091
la función de costo incremental toma una forma lineal por intervalos:

0, 0176 · PD + 10, 3330 PD ≤ 75, 9091
λ= (8.86)
0, 0067 · PD + 11, 1609 PD ≥ 75, 9091
Es posible aproximar estas funciones a una única lı́nea recta mediante

un proceso de mı́nimos cuadradados en un intervalo determinado, para
ello se deriva la siguiente función de error con respecto a aeq y beq :
Pc Pm
2
e= (aeq · PD + beq − λ(PD )) dPD + (aeq · PD + beq − λ(PD ))2 dPD
0 Pc
(8.87)
Para el caso (Pc = 75, 9091) con una potencia máxima de Pm = 300M W
los valores equivalente son aeq = 0, 0084 y beq = 10, 7948.
157
1. Allen J. Wood; Bruce F. Wollenberg. Power generation, operation, and

control John Wiley & Sons.
2. John J. Grainger; William D. Stevenson. Análisis de sistemas de

potencia. McGraw-Hill.
3. Mohamed E El-Hawary. Electrical power systems (design and analysis).

IEEE PRESS Power systems engineering series.
4. Dhillon Kothari. Power System Optimization. Prentice Hall of India.

2004.
Operación de sistemas
9
hidrotérmicos.
Contenido
Modelo lineal de largo plazo
Programación dinámica para problemas

de mediano plazo
Técnica exacta para problemas de corto

plazo
159
160 CAPÍTULO 9. OPERACIÓN DE SISTEMAS HIDROTÉRMICOS.
Humedo Normal
Utilizar
Embalses
S eco Deficit
Decisión
Humedo Vertimiento
No
Utilizar
Embalses
S eco Normal
Caudales
Hoy Consecuencia
Futuros
Figura 9.1: Dependencia temporal del problema
U NA forma de energı́a renovable abundante en la naturaleza es la

energı́a hidráulica resultante del aprovechamiento de las caı́das de agua.
En Colombia, es posible aprovechar este recurso gracias a las riqueza
hidrológica y la topografı́a montañosa. Sin embargo la generación térmica
es necesaria para disponer de energı́a eléctrica en tiempos de sequı́a.
El despacho hidrotérmico determina el cronograma de generación

hidráulica y térmica, de tal forma que se minimice el costo operativo a
lo largo de un horizonte de planeamiento. La necesidad de optimización
se debe a las caracterı́sticas del sistema de generación. La generación
hidráulica presenta un costo operativo bajo y una dependencia con
la hidrologı́a -variable estocastica- la cual introduce un riesgo de
racionamiento asociado con las temporadas de sequı́a. De otro lado, la
generación térmica presenta un riesgo inferior con un costo operativo alto.
A diferencia de los sistemas puramente térmicos, los sistemas

hidrotérmicos presentan caracterı́sticas dinámicas, es decir, una decisión
161
operativa presente tiene influencia en el costo de la operación futura del

sistema. Esto conlleva a que el planeamiento de corto plazo requiera los
resultados del mediano y éste a su vez los del largo plazo.
En la figura 9.1 se muestra un esquema del acoplamiento temporal

del problema. Se evidencian dos comportamientos crı́ticos: utilizar los
embalses al máximo puede llevar a un déficit de energı́a y racionamiento
en caso de seguir una temporada seca. De otro lado no utilizar
los embalses puede aumentar innecesariamente los costos operativos
presentes y producir vertimientos en la operación futura (lo cual implica
una pérdida del recurso hı́drico).
En la gráficas 9.2 se muestra las funciones de costo presente y futuro

de un sistema hidrotérmico. La función de costo presente determina
los costos de la operación con respecto al volumen almacenado, cuanto
mayor cantidad de agua se almacene, será menor la cantidad de agua
turbinada para la generación de energı́a, esto significa que la mayor parte
de la demanda debe ser suministrada por parte de la plantas térmicas
las cuales tienen un costo de operación alto, es por esto que la curva es
decreciente. De otra parte, si la cantidad de agua desembalsada es menor
se obtendrá un mayor volumen acumulado para los periodos siguientes,
esto significa que la operación futura puede ser abastecida por más
generación hidráulica desplazando las plantas térmicas y disminuyendo
por ende los costos operativos. El uso optimo del agua se obtiene cuando
están equilibrados los valores inmediatos y futuros. Es por esto que el
despacho hidrotérmico debe realizarse en un horizonte de planeamiento
cordinado a largo, mediano y corto plazo.
La estocasticidad de las afluencias puede ser modelada de varias

formas, a saber:
• Mediante una adecuada predicción: Usando herramientas de

predicción basadas en la inteligencia artificial como las redes
neurodifusas o mediante modelos ARIMA 1 .
• Mediante programación estocástica (una separación en escenarios

con una probabilidad determinada.)
1
Autoregresive Integrated Moving Average
Costo
FCF+FCP
FCP FCF
Optimo Embalse
Figura 9.2: Funciones de costo futuro y presente
• Mediante programación dinámica dual estocástica2 .
• Mediante modelos de simulación como la simulación de montecarlo.
En algunas ocasiones solo se requiere conocer el comportamiento

óptimo del sistema frente a un escenario crı́tico en cuyo caso el problema
se hace deterministico. Igualmente, la separación en escenarios permite
considerar de forma simplificada el comportamiento aleatorio de las
afluencias. Normalmente se consideran tres escenarios: optimista, medio
y pesimista; cada escenario en particular es de tipo deterministico.
La primera opción hace determinı́stico el problema de despacho

hidrotérmico, no obstante, la predicción de los caudales es cada vez más
difı́cil en especial para estudios a largo plazo. Un elemento importante
a tener en cuenta es el grado de detalle esperado en el problema y la
coordinación entre los diferentes tipos de planeamiento. La solución del
problema de largo plazo (normalmente considerado de forma estocástica
y sin considerar la red) alimentan el problema de mediano plazo y este a
su vez al problema de corto plazo. El despacho hidrotérmico constituye el
denominado planeamiento energético, en Colombia este planeamiento es de
tipo indicativo. Para problemas de planeamiento eléctrico (Planeamiento
2
La programación dinámica dual corresponde a la aplicación iterativa de una técnica
denominada corte de Benders. Para más información ver referencia [2]
163
de la expansión de la transmisión) el resultado del despacho hidrotérmico

es una variable de entrada, como en este caso la red de transmisión debe
ser sujeta a modificaciones, el despacho hidrotérmico debe realizarse a
nodo único. Para problemas de simulación de la operación, el despacho
debe realizarse considerando la red ya sea mediante un flujo DC o por
medio de un modelo de transportes.
En los sistemas térmicos las unidades presentan un costo directo, esto

significa que su costo depende solo de su nivel de operación y no de
las demás unidades como ocurre con las centrales hidráulicas. En los
sistemas hidrotérmicos, aunque el agua no presenta un coste directo, se
puede calcular un costo de oportunidad el cual está dado por la cantidad
de generación térmica desplazada.
Los modelos de despacho para corto y largo plazo difieren en el

detalle del modelo matemático utilizado y en la información de entrada
al problema. El despacho de largo plazo puede ser realizado por medio
de un modelo lineal el cual asume constante la caı́da H por lo cual se
establece una relación proporcional entre el caudal y la potencia generada.
Esta aproximación es normalmente válida en el caso de centrales con
turbinas pélton, en donde la casa de máquinas se encuentra a distancia
considerable de la toma y las variaciones en el nivel del embalse no
afectan significativamente la caı́da total. El uso de las metodologı́as de
programación dinámica permite tener en cuenta la relación no lineal entre
el caudal y la potencia generada por lo cual esta metodologı́a puede ser
usada para estudios de mediano y corto plazo. No obstante, la técnica de
multiplicadores de Lagrange puede ser usada con éxito en los problemas
de corto plazo en donde el volumen disponible es dado por el resultado
del despacho de largo plazo.
Las plantas térmicas presentan una función de costo cuadrático como

se mostró en el capı́tulo anterior, no obstante, es posible realizar una
simplificación lineal del problema como se describe en la siguiente sección.
9.1 Modelo matemático lineal
El despacho hidrotérmico determinı́stico puede ser modelado como un

problema de programación lineal si se considera el costo del combustible
como una función lineal al igual que la relación entre caudal y potencia
generada. A continuación se muestra de forma detallada el modelo
matemático.
Función objetivo:
N p NT
M in Cj · PT (jt) (9.1)
t=1 j=1
Conjunto de restricciones:
PH(it) = τi · Qit (9.2)
NT NH
PT (jt) + PH(it) = PD(t) (9.3)
j=1 i=1
Vi,t+1 = Vi,t + T · (Ait − Qit − Sit ) (9.4)
Vitmin ≤ Vi,t ≤ Vitmax (9.5)
Qitmin ≤ Qit ≤ Qitmax (9.6)
PT (jt)min ≤ PT (jt) ≤ PT (jt)max (9.7)

165
donde
• Np : número de periodos a condiderar.
• NT : número total de plantas térmicas.
• NH : número total de centrales hidroeléctricas.
• PD(t) : Demanda total en el periodo t.
• PH(it) : generación hidráulica para la planta i en el periodo p.
• PT (jt) : generación térmica para la planta j en periodo p.
• Vi,t : volumen del embalse i en periodo t.
• Qit : caudal turbinado en embalse i en periodo t.
• Sip : vertimientos asociados al embalse i en el periodo t.
• Ait : afluencias en el embalse i para el periodo t.
• Cj : costo unitario de la generación térmica en la planta j.
• τi : factor de turbinamiento en la unidad hidroeléctrica i.
• T : periodo de tiempo considerado.
Este modelo puede ser resuelto usando una técnica de programación

lineal como simplex o puntos interiores.
Se invita al lector a realizar la implementación en Gams del modelo

propuesto.
St
At
Vt
Qt
Figura 9.3: Variables hidráulicas
9.2 Modelo no lineal de mediano plazo
Como se mostró en el capı́tulo anterior la función de costo asociada a

cada central térmica presenta un forma cuadrática y no lineal como supone
el modelo presentado de despacho hidrotérmico. Una metodologı́a
relativamente eficiente para dar solución al problema de despacho
hidrotérmico es la programación dinámica. El problema de optimización
en este caso consiste en minimizar en cada periodo el costo presente (CP)
mas la función de costo futuro (FCF).
Como se mostró en capı́tulos anteriores la energı́a potencial obtenida

mediante el almacenamiento de agua en el reservorio se transforma en
energı́a eléctrica por medio del sistema turbina-generador. El agua ingresa
a través de la bocatoma recorriendo la tuberı́a de presión hasta la casa de
máquinas en donde el proceso de transformación es realizado, finalmente,
el caudal es extraı́do por medio de un canal de desfogue.
El agua represada puede ser retirada sin generar energı́a por medio del
vertedero de excesos el cual permite controlar el nivel del embalse.
El volumen de agua almacenada en el reservorio en el periodo t es

representado por la variable Vt . Según su uso se clasifica en volumen
muerto, útil o de seguridad.
167
El volumen muerto no puede ser aprovechado en la generación de

energı́a. El volumen de seguridad es el volumen existente entre el nivel
máximo operativo y el nivel máximo de aprovechamiento. Por encima
del nivel de aprovechamiento, la presa puede colapsar lo cual generarı́a
vertimientos.
Los reservorios de acuerdo a su capacidad de regulación del caudal son

clasificados como de acumulación y de compensación. Los reservorios
de acumulación tienen gran capacidad de almacenar energı́a en forma
potencial (agua), por lo cual son responsables de la regulación del rı́o,
contrario a lo que ocurre en los reservorios de compensación. Las
represas con reservorios de acumulación son generalmente denominadas
de reservorio mientras que las con reservorio de compensación son
denominadas a filo de agua.
Las centrales a filo de agua son restadas directamente de la demanda

por lo que no son incluidas en el proceso de optimización.
Las cascadas hidrologicas tambien deben ser representadas en el

modelo de despacho multiembalse. en la figura 1 se muestra una cadena
hidrológica tı́pica.
Un dato importante en el análisis de un reservorio son las afluencias

como se mostró en capı́tulos anteriores. El volumen afluente se clasifica
en dos tipos: natural e incremental. Las afluencias naturales constituyen
las que vienen del afluente sin tener en cuenta las demás usinas. Las
afluencias incrementales corresponden al acumulado de los vertimientos
y los turbinamientos de las centrales ubicadas aguas arriba en la cadena
hidrológica
El volumen almacenado puede colocarse en función de la caida liquida

de la central por medio de un polinomio. como se mostró en el capitulo 4.
En el caso de centrales con turbinas de reacción se puede obtener

un polinomio similar que relaciona la altura del canal de salida vs el
turbinamiento, Esta relación debe ser llevada en cuenta en el proceso de
optimización.
Es posible obtener una simplificación en el proceso linealizando estos

polinomios, el modelo resultante es una aproximación. Ası́, la relación

entre el caudal turbinado y la potencia generada por cada unidad
hidráulica es expresada mediante una función lineal como se muestra en
la ecuación 9.8:

τo + τ · PH si PH > 0
Q= (9.8)
0 si PH = 0
No obstante, la metodologı́a de programación dinámica permite la

inclusión de funciones no lineales sin aumentar significativamente la
complejidad del algoritmo. Una vez determinado el modelamiento no
lineal del problema es posible dar solución al mismo. Para ello es necesario
construir el grafo que representa la discretización del embalse y darle
solución por medio del principio de optimalidad de Bellman.
Inicialmente se determina el caudal requerido de acuerdo a la

discretización elegida para cada una de los enlaces mediante la ecuación
9.9:
(Vt − Vt−1 )
Q = At − (9.9)
T
Con este valor de caudal es posible determinar la potencia activa

generada por la central hidráulica mediante la ecuación 9.8. Si la
potencia activa es superior a la demanda del sistema entonces se producen
vertimientos. Esto no significa necesariamente que la solución definitiva
presente vertimientos, lo cual ocurrirı́a solo si la ruta óptima pasa por
la rama correspondiente. Una vez determinada la generación hidráulica
es posible determinar la generación térmica correspondiente simplemente
restando de la demanda:
PT (t) = PD(t) − PH(t) (9.10)
Los costos asociados a la rama son determinados con base en la

ecuación de costos de la planta térmica.
169
Ejemplo 9.1
Encontrar el despacho óptimo para un tiempo de 24h discretizado en periodos

de 8h en un sistema conformado por una unidad hidráulica y una unidad térmica.
El volumen inicial es de 2 · 106 m3 . (Discretizar el embalse en pasos de 1 · 106 m3 ).
Se espera que el volumen al final del periodo de planeamiento sea como mı́nimo
igual al inicial. Las caracterı́sticas del sistema son las siguientes:
Planta hidráulica : El embalse equivalente presenta un volumen mı́nimo

de 1 · 106 m3 y un volumen máximo de 3 · 106 m3 , la relación caudal
potencia generada esta dada por: Q = 2 + 0.1 · PH .
Planta térmica : La planta térmica presenta una caracterı́stica de costos

1
f = 100 · PT2 + 4 · PT con una potencia máxima de 500M W .
Demanda para cada periodo : Se consideraran tres periodos con

demandas de 300M W ,300M W y 800M W .
Afluencias : Las afluencias en los tres periodos considerados son de

35m3 /s, 30m3 /s y 40m3 /s
El costo por racionamiento es de $15/M W
La figura 9.4 muestra el grafo para el problema descrito
El procedimiento para calcular los costos de cada rama se muestra

a continuación, en este caso el estado H puede ser descartado ya que
la restricción principal es que se tenga al menos el mismo volumen de
partida. Cualquier solución que termine en H será mayor o igual a la
solución óptima que termina en I.
Transición AB :
3·106 −2·106
• S + Q = 35 − 8·3600
= 0, 277m3 /s
B E H
A C F I
D G
Figura 9.4: Grafo para el problema de despacho hidrotérmico
• Este valor está por debajo del mı́nimo necesario para generar
(2m3 /s) por lo cual la totalidad de la generación debe ser
asumida por la planta térmica:
• S = 0.277m3 /s
• PH = 0
• PT = 300M W
• Costo = $2100
Transición AC : Vi = Vf
• S + Q = 35m3 /s
• PH(disp) = 330M W
• PH = 300M W
• S = 3m3 /s
• PT = 0
• Costo = $0
171
Transición AD :
1·106 −2·106
• S + Q = 35 − 8·3600
= 69.722m3 /s
• PH(disp) = 677M W
• PH = 300M W
• PT = 0
• S = 37, 72m3 /s
• Costo = $0
Transición BE : Vi = Vf
• S + Q = 30m3 /s
• PH = 280M W
• PT = 20
• Costo = $84
Transición CE :
3·106 −2·106
• S + Q = 30 − 8·3600
= −4.722
• No hay capacidad de generación hidráulica para pasar de C a
E. Se podrı́a pasar a un estado intermedio entre E y F pero éste
no fue considerado en la discretización planteada (éste error es
el costo de la discretización en este problema).
• Costo = ∞
Transición DE :
3·106 −1·106
• S + Q = 30 − 8·3600
= −39.44
• Las afluencias no son suficientes para esta transición
• Costo = ∞
Transición BF :
2·106 −3·106
• S + Q = 30 − 8·3600
= 64.722m3 /s
• PH(disp) = 627.22M W
• PH = 300M W
• S = 32.722m3 /s
• PT = 0M W
• Costo = $0
Transición CF : Equivale a la transición BE Costo = $84
Transición DF : Equivale a la transición CE Costo = ∞
Transición BG :
1·106 −3·106
• S + Q = 30 − 8·3600
= 99.44m3 /s
• PH(disp) = 974M W Supera la demanda
• PH = 300M W
• S = 942.44m3 /s
• PT = 0M W
• Costo = $0
Transición CG : Equivale a la transición BF Costo = $0
Transición DG : Equivale a la transición BE Costo = $84
Transición EI :
2·106 −3·106
• S + Q = 40 − 8·3600
= 74.722m3 /s
• PH(disp) = 727.22M W
• PH = 727.22M W
• PT = 72.78M W
• Costo = $344
Transición FI :
• S + Q = 40m3 /s
• PH = 380M W
• PT = 420M W
• Costo = $3444
Transición GI :
173
2·106 −1·106
• S + Q = 40 − 8·3600
= 5.27m3 /s
• PH = 32.78M W
• PT n = 767.22M W Supera la capacidad de la planta
• PH = 500M W
• PR = 267.22M W (Racionamiento)
• Costo = $4012.8 + $4500
Los casos de costo infinito son rutas no factibles que el método descarta
de forma automática. Una vez definidos los costos de las diferentes
rutas se procede a usar programación dinámica para dar solución al
problema. En este caso la ruta óptima es ABEI cuyo costo es de $2528.
Los vertimientos son minimizados de forma implı́cita por lo que éstos
no requieren ser penalizados. La caracterı́stica dinámica del sistema se
hace evidente en la ruta ACGI la cual presenta costos iguales a cero en los
primeros periodos y sin embargo, es la más costosa además de presentar
racionamiento.
9.2.29 Costo del

El costo incremental del agua corresponde al multiplicador de agua.
Lagrange (o a la variable dual, en el caso del modelo lineal) asociado
a la restricción de caudal. Este costo incremental indica el cambio
en la función objetivo frente a una variación unitaria del caudal. En
múltiples análisis, este costo indica el valor cuantitativo del agua (costo de
oportunidad) y puede ser determinado para especificar la oferta de cada
central hidráulica.
Dada la relación entre caudal y potencia hidráulica es posible plantear

la siguiente aproximación:
ΔQt ∂Q
= (9.11)
ΔPH(t) ∂PH
Igualmente, de la ecuación de balance de potencia (ec 9.10) se obtiene:

ΔPT (t) = −ΔPH(t) (9.12)
esta relación indica que cada MW generado por parte de la central

hidráulica desplaza un MW de generación térmica. Finalmente, de la
función objetivo se tiene:
Np
∂f
ΔF = · ΔPT (t) (9.13)
t=1
∂PT (t)
reemplazando
Np
∂f /∂PT (t)
ΔF = − · ΔQ (9.14)
t=2
∂Q/∂PH(t)
en donde es posible definir el costo de oportunidad del agua para cada

uno de los periodos:
∂f /∂PT (t)
υt = − (9.15)
∂Q/∂PH(t)
el signo negativo indica que cada unidad de caudal adicional genera

una disminución en los costos totales de operación, igual a υt . Este valor
puede ser usado dentro de un algoritmo de optimización para definir
cual de los periodos deben tener mayor nivel de turbinamientos. Desde
luego las restricciones de desigualdad pueden modificar estos costos de
oportunidad, pero en general, es un buen indice de sensibilidad para
definir el valor del agua ya que ésta, como combustible, tiene un costo
igual a cero, ası́ que el costo del agua está asociado a la cantidad de
generación térmica que desplaza.
El valor del costo de oportunidad del agua es lo que al despacho

termico son los multiplicadores de lagrange, por tal razón, se constituye en
un indicador de sensibilidad de tipo económico en el análisis del mercado.
175
9.2.30 Múltiples
Hasta el momento se ha mostrado el proceso de optimización para embalses
un sistema simple conformado por una planta hidráulica y una planta
térmica. En el caso de tener varias plantas térmicas se puede recurrir a
un procedimiento de despacho económico como se mostró en el capı́tulo
anterior. Igualmente, un sistema de varios generadores térmicos puede
ser llevado a un generador térmico equivalente (ver sección 8).
En el caso de tener varias centrales hidroeléctricas es posible realizar un

embalse equivalente si los periodos de lluvia y estiaje en cada uno de los
embalses es similar. En caso contrario es posible utilizar la aproximación
de la programación dinámica optimizando cada una de las centrales de
forma independiente dejando las demás constantes. Por ejemplo, para dos
centrales hidráulicas de generación cuyas soluciones están dadas por X1 y
X2 y la función de costos es F (X1 , X2 ) se tiene el siguiente procedimiento,
similar al método de aproximaciones sucesivas de Gauss:
Inicialmente se encuentra una solución factible X2 (no necesariamente

óptima) y se da solución al problema F teniendo X2 constante. La
metodologı́a de solución puede ser programación dinámica convencional
pues la solución X2 se resta directamente de la demanda del sistema.
La solución encontrada por esta metodologı́a alimenta nuevamente el
problema de optimización pero esta vez el valor de X1 es constante y es el
valor de X2 la variable de decisión.
Uno de elementos fundamentales de este método, es la maldición

de la dimencionalidad pues cuanto mayor es el numero de estados
(discretizaciones del embalse) más exacta será la metodologı́a, pero el
esfuerzo computacional se hace extremadamente alto.
La red puede ser igualmente modelada mediante un flujo DC o por

medio de un modelo de transportes. De esta forma se pueden considerar
las restricciones de cargabilidad de las lı́neas de transmisión.
9.3 Modelo no lineal de corto plazo
El problema de corto plazo es alimentado con los resultados de los

problemas de mediano y largo plazo; en este caso el volumen disponible
es un dato caracterı́stico del problema por lo cual la principal variable de
interés es el despacho de la generación hidráulica y térmica.
Las plantas térmicas son representadas nuevamente por una forma

cuadrática dada por una función f (PT ). En el caso de existir
varias unidades térmicas éstas pueden ser modeladas por una térmica
equivalente como se mostró en el capitulo anterior o re-plantear el proceso
de optimización para considerar cada central.
a
F (PT ) = · P 2 + b · PT + c (9.16)
2 T
De forma análoga, las plantas hidráulicas pueden ser representadas

por una forma cuadrática entre el caudal turbinado y la potencia generada:
α
Q(PH ) = · PH2 + β · PH + γ (9.17)
2
El problema consiste en encontrar una ruta óptima para la generación

hidráulica y térmica tal que se cumpla el balance de potencia y se utilice
el volumen disponible. El carácter dinámico del problema hace que la
solución a encontrar sea una función y no un valor constante como en
el caso del despacho de plantas térmicas, esto significa determinar las
funciones PT (t),PH (t), Qt y F (t). El modelo matemático se presenta a
continuación para un periodo de tiempo T y un volumen disponible Vm :
T
min F (t) · dt (9.18)
0
restricción de balance de potencia para cada uno de los puntos en el

177
intervalo
PH (t) + PT (t) = PD (t) ∀t ∈ [0T ] (9.19)
restricción de volumen total utilizado
T
Q(t)dt = Vm (9.20)
0
La restricción de balance de potencia puede ser introducida en una

función aumentada teniendo en cuenta el multiplicador de Lagrange para
cada uno de los puntos t, en este caso, λ es una función del tiempo
pues la restricción se debe cumplir para cada uno de los puntos t. Con
respecto a la restricción de volumen, el multiplicador de lagrange es único
y representa el coeficiente de conversión del caudal. La función toma
entonces la siguiente forma:
T
min J = {F (t) + λ(t) · (PD (t) − PH (t) − PT (t)) + y · Q(t)} dt (9.21)
0
derivando con respecto a cada una de las variables del problema se

obtienen las condiciones de optimalidad
∂
F (t) − λ(t) = 0 ∀t ∈ [0 − T ] (9.22)
∂PT
∂
y· Q(t) − λ(t) = 0 ∀t ∈ [0 − T ] (9.23)
∂PH
PH (t) + PT (t) = PD (t) ∀t ∈ [0 − T ] (9.24)
T
Q(t)dt = Vm (9.25)
0
El sistema de ecuaciones resultante puede ser solucionado mediante

un método numérico aunque en algunos casos es posible darle solución
de forma analı́tica.
Ejemplo 9.2
Mostrar el despacho óptimo de un sistema conformado por una unidad

hidráulica y una unidad térmica para un periodo de 24h. La curva de carga es
mostrada en la figura 9.5, el volumen total es de 357500m3
800 MW
450 MW
300 MW
8h 16 h 24 h
Figura 9.5: Curva de carga para el problema de despacho
Las funciones de costo y turbinamiento son respectivamente:
F = 2 · PT + 0, 001 · PT2 (9.26)
Q = 50 · PH + 0, 01 · PH2 (9.27)
179
Usando las condiciones de optimalidad resultantes se obtiene:
λ(t) = 2 + 0, 002 · PT = y · (50 + 0, 01 · PH ) (9.28)
Teniendo en cuenta el balance de potencia, es posible expresar la

potencia hidraulica como función de coeficiente de conversión de caudal:
y · (50 + 0, 01 · PH ) − 2
PH + = PD (t) (9.29)
0, 002
Finalmente, teniendo en cuenta la restricción de volumen (ecuación

9.25) se obtiene el valor de y:
8 16 24
Qt dt + Qt dt + Qt dt = 357500 (9.30)
0 8 16
−3971 + 50000 · y + 250000 · y 2

− 600 · = 357500 (9.31)
(1 + 10 · y)2
esta ecuación presenta dos soluciones de las cuales se elige la solución

positiva: y = 0, 0446. Con este valor es posible determinar la generación
hidráulica y térmica en cada intervalo:
⎧ ⎫
⎨ 128, 3711 M W 0 ≤ t ≤ 8 ⎬
PH = 474.2094 M W 8 ≤ t ≤ 16 (9.32)
⎩ ⎭
232.1226 M W 16 ≤ t ≤ 24
⎧ ⎫
⎨ 171.6289 M W 0 ≤ t ≤ 8 ⎬
PH = 325.7906 M W 8 ≤ t ≤ 16 (9.33)
⎩ ⎭
217.8774 M W 16 ≤ t ≤ 24
Ejemplo 9.3
Un sistema de generación está compuesto por una central hidráulica y una

central térmica. La central térmica presenta una relación de costos de F = 1, 5·PT2
en donde PT es la potencia generada en M W . La central hidráulica tiene un
volumen disponible de 10800m3 y una relación caudal potencia generada de Q =
1, 8 · PH en donde PH está dado en M W y Q en m3 /s. Calcular el despacho
óptimo para un periodo de 24 horas si la demanda del sistema es la mostrada en la
figura 9.6
600 MW
300 MW
24 h
Figura 9.6: Demanda variable del sistema
La demanda del sistema puede ser representada por la ecuación:
300
PD = 600 − ·t (9.34)
24
La metodologı́a de solución es igual al ejemplo anterior salvo que

ésta vez la demanda presenta una forma continua. Las condiciones de
optimalidad siguen siendo las mismas:
λ(t) = 3 · PT (9.35)
181
λ(t) = 1, 8 · y (9.36)
24
300 1, 8
1, 8 · 600 − ·t− dt = 10800 (9.37)
0 24 3
de donde el valor de y es 333, 33 Con este valor es posible determinar

la generación en la central hidráulica y térmica:
PT = 200M W
PH = 400 − 12, 5 · t
Costo : 1 440.000
1. Mohamed E El-Hawary. Electrical power systems (design and analysis).

IEEE PRESS Power systems engineering series.

10
El mercado colombiano.
Contenido
El mercado eléctrico
Bolsa vs contratos
Restricciones de la red
Poder del mercado.
183
184 CAPÍTULO 10. EL MERCADO COLOMBIANO.
E Lsistema eléctrico Colombiano es operado mediante un enfoque de

mercado, por tal razón es importante tanto para las empresas generadoras
como para los comercializadores, entender el funcionamiento del mercado
eléctrico. En este capı́tulo se presenta una descripción general del
funcionamiento del mercado y de los procesos de despacho y re-despacho
del mismo, usando los conceptos de optimización presentados en
capı́tulos previos.
10.1 Modelo de mercado
La competencia es un elemento fundamental en el desarrollo de los

mercados, ésta permite la optimización implı́cita del mismo al aumentar
la calidad del servicio y reducir el precio. Una competencia perfecta
se presenta cuando ningún agente o empresa es capaz de influir sobre
el precio, para ello, generalmente se requiere el cumplimiento de las
siguientes condiciones:
• un gran número de competidores.
• perfecta información: Todos los competidores saben del precio de

sus rivales y todos los compradores saben de las ofertas de los
productores.
• la empresa puede vender la cantidad que desea de su producto.
El sistema eléctrico constituye un mercado atı́pico en la medida en que

se aleja totalmente de estas condiciones:
• el número de competidores es limitado y relativamente constante
• depende de variables aleatorias como la hidrologı́a
• la cantidad de generación depende de la demanda y esta es a su vez

inelástica
185
• imposibilidad de almacenamiento del producto

• la discontinuidad del servicio es intolerable
A diferencia de los demás mercados, el mercado eléctrico no sufre

variaciones considerables de demanda frente a modificaciones en el
precio, esto se conoce como inelasticidad de la demanda.
Dada estas circunstancias parece imposible obtener competencia en

un sistema eléctrico, no obstante, es posible crear un mercado eléctrico
mediante un axioma simple: la energı́a eléctrica como producto puede separarse
comercialmente de su transporte como servicio. Esto se logra mediante dos
esquemas de mercado: el mercado spot y los contratos bilaterales.
El mercado spot se refiere a la bolsa de energı́a, la cual define el

precio nodal frente a las ofertas y la demanda horaria. El mercado
spot es altamente volátil y puede ocasionar riesgos importantes para los
diferentes agentes del mercado. Existen varios modelos de mercado spot,
en el caso Colombiano la bolsa de energı́a presenta un tipo de despacho a
nodo único basado en el concepto de orden de méritos que será estudiado
mas adelante.
Los contratos bilaterales constituyen un mecanismo de compra y venta

de energı́a a largo plazo, por tal razón constituye un mecanismo de
cubrimiento frente al riesgo. Los contratos pueden ser fı́sicos o financieros.
En los primeros, el contrato obliga a la empresa generadora a abastecer
con sus recursos propios el valor de la energı́a vendida, mientras que en
los contratos financieros, la energı́a vendida por una empresa puede ser
suministrada por otra mediante una segunda compra-venta de energı́a
que puede ser realizada por medio de un contrato o mediante la bolsa
de energı́a.
Un mercado eléctrico puede tener una combinación de los dos modelos

- bolsa y contratos-, como es el caso colombiano, en el cual se tiene la figura
de la bolsa de energı́a y los contratos bilaterales de tipo financiero.
No debe confundirse la libre competencia con la ausencia de regulación

o intervención gubernamental. El fomento de la competitividad entre
empresas requiere frecuentemente una normativa reguladora minuciosa.
En el caso del sistema colombiano la regulación está a cargo de la CREG

mientras que el planeamiento de la transmisión lo realiza la UPME.
Puesto que el mercado supone separar el mercado eléctrico del servicio

de transporte de energı́a, el planeamiento se realiza en dos frentes:
planeamiento energético y planeamiento eléctrico.
El planeamiento energético se refiere básicamente al uso de los recursos

de generación, por ejemplo, mediante el despacho hidrotérmico. No
obstante, en la medida que se garantice la libre competencia de los
agentes generadores, el despacho hidrotérmico presenta un planeamiento
puramente indicativo.
De otro lado, el servicio de transporte de energı́a eléctrica, debe ser

centralizado por parte de la UPME que define el planeamiento eléctrico,
éste a diferencia del anterior es obligatorio debido al monopolio natural
de la transmisión. El planeamiento eléctrico determina las nuevas
interconexiones a realizarse en el sistema, ası́ como la entrada en operación
de nuevas centrales y subestaciones.
10.1.31
Colombia, al igual que la mayorı́a de los paı́ses, decidió modificar su
Concepto de
estructura regulatoria a un modelo de mercado. La tendencia básica en
competencia
los procesos de reestructuración es promover la competencia de agentes
privados. En la mayorı́a de los casos, la competencia ha reemplazado los
procedimientos de planeamiento de la expansión basados en optimización
centralizada.
Estos cambios son posibles gracias a los avances en la tecnologı́a de

los sistemas de generación. Históricamente, estos han sido un monopolio
natural que impedı́an la competencia. No obstante, gracias al aumento
en la eficiencia de las plantas térmicas y a la reducción en los precios
del gas, las plantas térmicas pueden entrar en competencia con las
plantas hidráulicas en algunas estaciones del año. Igualmente, en muchas
circunstancias es mas eficiente la construcción de nuevas centrales de baja
capacidad a la repotenciación de unidades existentes, este hecho aumenta
el número de participantes en el mercado lo cual va en beneficio de la
187
competencia.
Además de Colombia, paı́ses como Chile, Reino Unido, Nueva

Zelanda y regiones de los Estados Unidos como California, han sometido
sus respectivos sistemas a procesos de privatización o reestructuración.
Existe una clara diferenciación entre estos dos conceptos: un sistema
puede ser reestructurado en el sentido de modificar las reglas con que
se compra y se vende la electricidad promoviendo la competencia; este
proceso se puede realizar manteniendo o no el control estatal sobre los
activos del sistema.
En la teorı́a económica se tiene el concepto de competencia perfecta

en donde cada uno de los productores por sus caracterı́sticas no puede
influir por medio de sus decisiones sobre el precio de mercado. Esto
significa que cada productor presenta un producto al mercado con iguales
caracterı́sticas de:
• Tiempo y oportunidad
• Posición geográfica
• Calidad del producto
Igualmente se asume que cada empresa tiene aproximadamente el

mismo tamaño y poder en el mercado. En general, cada uno de
los generadores producen el mismo producto en cuanto a calidad, sin
embargo, la posición de los generadores puede influir en el mercado en
la medida en que generadores aislados aumentan las pérdidas técnicas
y las congestiones en las lı́neas (Cargabilidad). De igual forma, con los
problemas ambientales suscitados por la contaminación atmosférica, es
posible ver preferencia por las fuentes de energı́a limpia. Ası́ mismo, la
transmisión se constituye en un monopolio natural en donde elementos
técnicos como la estabilidad y la confiabilidad pueden favorecer ciertas
unidades de generación. El mercado eléctrico es entonces desregulado en
cuanto a la generación pero intervenido en la transmisión.
El negocio de la generación es un oligopolio en la medida en que

un pequeño número de empresas pueden influir de forma significativa
en el mercado mientras que otros deben someterse a las condiciones

establecidas por los primeros.
Para separar cada uno de los agentes del mercado, se crea el concepto
de bolsa de energı́a la cual determina los despachos óptimos desde el
punto de vista de los precios mientras que las restricciones de transmisión
son consideradas en una etapa de optimización posterior.
La generación de energı́a puede verse como un producto mientras que

el control de tensión y frecuencia, ası́ como los requisitos de cargabilidad,
seguridad y confiabilidad se presentan como servicios auxiliares.
10.1.32 Bolsa de Un dı́a antes de la operación, los generadores reportan al sistema

energı́a la generación disponible para cada hora del dı́a y el precio de oferta
correspondiente. El operador del mercado realiza un despacho económico
para cada hora de acuerdo al modelo mostrado a continuación:
N
M in f = λj · Pjt (10.1)
j=1
sa
N
Pjt = Dt (10.2)
j=1
Pjt(min) ≤ Pjt ≤ Pjt(max) (10.3)
en donde
Pjt : Potencia generada por la planta j en el intervalo t
λj : Precio de oferta del generador j
Dt : Demanda para el periodo t
Pjt(min) , Pjt(max) : Lı́mites operativos de la planta j (Disponibilidad)

189
Como se puede observar este modelo de despacho no tiene en cuenta

la red (nodo único). Las cadenas hidráulicas, por su condición de
acoplamiento, deben ofertar de forma conjunta.
Este problema de programación lineal puede ser fácilmente

solucionado por los métodos tradicionales como el simplex o puntos
interiores. Al igual que en el caso del despacho clásico de unidades
térmicas, en este modelo, los multiplicadores de Lagrange presentan
una interesante interpretación económica, pues el precio de bolsa (spot)
es calculado como la derivada de la función de costo con respecto a
demanda:
∂f
πt = λ∗t = (10.4)
∂Dt
En este caso, el precio de bolsa πt = λ∗t es un costo incremental dado por

el precio de oferta de la última unidad despachada j ∗ . Este concepto de
costo incremental o costo marginal permite determinar el precio de bolsa
mediante un conjunto de reglas simples o despacho ideal.
El despacho ideal determina la generación óptima por parte de cada

una de las unidades de generación, mediante el concepto de orden de
méritos. Las plantas declaradas como inflexibles -aquellas en las que
Pjt(min) = Pjt(max) - son despachadas inicialmente; luego, las centrales son
despachadas desde la mas económica a la mas costosa hasta completar
la demanda teniendo en cuenta la capacidad máxima de generación de
cada una. La última planta no inflexible despachada determina el precio
de bolsa, si ésta planta presenta una inflexibilidad -normalmente por el
lı́mite inferior- se declara inflexible en ese valor y se realiza nuevamente
el despacho; el resultado obtenido por esta metodologı́a es el mismo del
modelo matemático antes presentado.
En un mercado basado en bolsa cada planta recibe como ingresos el

precio de reventa el cual está dado por el producto entre su generación y
el precio de bolsa. De esta forma la ganancia neta de cada planta Rj está
dada por la diferencia entre el precio de reventa y sus costos operativos
como se muestra en la ecuación 10.5.
24
Rj = (πt − cj ) · Pjt (10.5)
t=1
Existen algunas reglas básicas en el mercado Colombiano las cuales son

definidas por la CREG:
• Las cadenas hidráulicas ofertan de forma conjunta, de no hacerlo el

CND asume el precio de oferta más bajo de la cadena.
• En caso de no tener información sobre la oferta, el CND asume los

datos del dı́a anterior.
• Cuando el nivel del embalse está por debajo de un mı́nimo operativo

superior el embalse es intervenido
• En caso de racionamiento, el precio de bolsa es definido de acuerdo

al costo de racionamiento asociado a esa hora.
• Los generadores que posean plantas o unidades de generación

conectadas al STN, con capacidad mayor o igual a 20M W están
obligados a ofertar para el Despacho Central.
• Los generadores que posean plantas menores, con capacidad mayor

o igual a 10M W y menor a 20M W , pueden optar por participar en
la oferta para el Despacho Central.
• Los autogeneradores producen energı́a eléctrica exclusivamente

para atender sus propias necesidades. Por lo tanto, no usan la red
pública para fines distintos al de obtener respaldo del Sistema.
• Igualmente, los cogeneradores pueden vender sus excedentes de

energı́a al sistema en cuyo caso se les obliga a mantener un factor
de potencia de 0, 9
En el caso del mercado colombiano, el modelo de bolsa es utilizado

en dos puntos distintos de la operación, como se muestra en la figura
10.1. El pre-despacho ideal determina la generación mas económica
para abastecer la demanda proyectada del sistema, se realiza un dı́a
191
antes de la operación por lo que sus datos de entrada son proyecciones

de la demanda y disponibilidad esperada. La red impone algunas
restricciones adicionales que obliga a realizar un despacho programado
levemente distinto al despacho ideal, este es el que se aplica propiamente
en el sistema, no obstante, se pueden presentar restricciones adicionales
causadas por contingencias o servicio de control terciario de AGC, esto
implica un redespacho en la operación final. Posterior a la operación, se
realiza el proceso de liquidación, el cual determina las desviaciones y las
reconciliaciones como se mostrará más adelante; en esta etapa se utiliza
nuevamente el modelo de bolsa de energı́a para determinar el precio de
bolsa solo que esta vez se utilizan los datos reales de la operación como
potencia demanda real y la disponibilidad comercial.
Pre-Despacho Despacho Undía antes

Ideal Programado de la operación
Día de la
Redespacho
operación
Despacho Posterior a la
Reconciliaciones
Ideal Operación
Restricciones
Nodoúnico de la red.
Figura 10.1: Modelo simplificado del despacho económico

10.1.33 Pérdidas En la gráfica 10.2 se muestra un sistema de sub-transmisión conectado

de potencia al sistema de transmisión nacional (STN) en donde su frontera comercial
es el nodo 102. El generador G1 puede abastecer parte de la demanda
del sistema pero no esta conectado directamente al STN (Generación
embebida), en estas situaciones existe una reglamentación para la
asignación de las pérdidas de potencia activa.
100 101
S TN
102
S istema
106
103 107
104 105
G1
Figura 10.2: Ejemplo de un sistema de subtransmisión
En el caso del sistema eléctrico Colombiano, existen medidores

localizados en diferentes puntos de la red que determinan los flujos
potencia desde y hacia la misma. Dichos medidores registran el valor
de MWh con dos cifras decimales teniendo en cuenta la dirección de
los flujos de potencia para ası́ determinar si se trata de importación o
exportación de energı́a. Las pérdidas se calculan como la diferencia entre
las exportaciones e importaciones totales al sistema de transmisión.
193
Cada comercializador asume en proporción a su demanda, una parte

de las pérdidas de energı́a; Ası́ mismo, una vez evaluada la distribución
de las pérdidas, es posible calcular la demanda comercial de cada
comercializador como la suma de su demanda propia y su participación
en las pérdidas. Igualmente, el comercializador debe asumir las pérdidas
en los niveles de tensión menores a los del STN. 1 .
En el caso de existir generación embebida en el sistema, el

comercializador asume las pérdidas asociadas a su generación, no
obstante, cuando la generación excede las necesidades del comercializador
es el generador el que asume las pérdidas necesarias para colocar su
excedente en fronteras del STN.
10.1.34
El modelo de bolsa de energı́a permite un despacho óptimo desde Restricciones y
el punto de vista de los costos teniendo en cuenta la capacidad de reconciliaciones
generación de las plantas; no obstante, al intentar realizar este despacho en
el sistema real ocurren violaciones (infactibilidades) dadas principalmente
por la capacidad de transporte de la red. Estas violaciones se denominan
restricciones y exigen que el despacho sea modificado.
El modelo de despacho con restricciones de cargabilidad genera como

resultado el despacho programado, el cual en una versión simplificada se
puede modelar como:
N
M in f = λj · Pjt (10.6)
j=1
sa
N
Pjt = Dt (10.7)
j=1
Pjt(min) ≤ Pjt ≤ Pjt(max) (10.8)
−Pijt(max) ≤ Pijt ≤ Pijt(max) (10.9)
1
Ver resolución CREG ...
en donde Pijt representa el flujo de potencia por la lı́nea ij en el

periodo t. Desde luego, este modelo puede ser complementado con las
restricciones de capacidad de potencia reactiva, lı́mites de tensión entre
otros.
El resultado del despacho programado puede diferir del despacho

ideal lo cual obliga al sistema a incurrir en costos adicionales. La
liquidación de estos costos se denomina reconciliación y puede presentarse
de dos formas a saber:
• Reconciliación positiva.
• Reconciliación negativa.
La reconciliación positiva se le paga a los generadores que son

despachados con un valor superior al del despacho ideal. Esta energı́a
complementaria es pagada a un precio de reconciliación dado por la
siguiente expresión:
Rec(+)j = min (λj , λmax ) (10.10)
en donde:
λj : Precio de oferta del generador j.
λmax : Tope definido por la CREG.
Un concepto igualmente importante son las desviaciones, cada recurso

de generación que se desvı́e del despacho programado horario, por
fuera de la franja de tolerancia del 5% verá afectadas sus transacciones
comerciales ası́:
El dinero que horariamente se determine en la bolsa de energı́a

por penalizaciones, se asigna a los comercializadores a prorrata de su
participación en la demanda total.
195
Tabla 10.1: Oferta de las unidades generadoras

Generador MW $/M W h
Térmica 1 80 50
Térmica 2 120 90
Hidráulica 1 90 35
Hidráulica 2 150 40
La cantidad de energı́a a generar por concepto de reconciliación

positiva debe ser compensada por otras plantas despachadas en orden
de merito, las cuales deben disminuir su generación real. Estos cambios
se denominan reconciliaciones negativas y se pagan a un valor dado
por el promedio entre el precio de bolsa y la oferta del generador
correspondiente:
λj + π t
Rec(−)j = (10.11)
2
Ejemplo 10.1
Determinar el pre-despacho ideal, el despacho programado y los precios de

reconciliación para el sistema mostrado en la figura 10.3 si los precios de oferta
están dados en la tabla 10.1. Las lineas pueden transportar máximo 50M W . Se
asume el costo máximo para la reconciliación positiva de $60 por MWh.
El despacho ideal es el siguiente:
Hidráulica 1 : 90 MW
Este despacho ideal puede ser encontrado mediante un sencillo

ordenamiento por méritos como se mostró anteriormente o con un modelo
en Gams como se muestra a continuación:
Térmica
Térmica
1
2
1 40 MW 2
60 MW
5
3 4
100 MW
Hidráulica Hidráulica
1 2
Figura 10.3: Ejemplo de un sistema de potencia

$title Ejemplo de despacho i d e a l
s e t s k generadores /H1 , H2 , T1 , T2 / ;
t a b l e Datos ( k , * ) Datos de o f e r t a
MW Costo
H1 90 35
H2 150 40
T1 80 50
T2 120 90;
parameter D demanda d e l s i s t e m a / 2 0 0 / ;
variable f funcion o b j e t i v o ;
positive variable
P(k ) p o t e n c i a generada ;
P . up ( k ) = Datos ( k , ’MW’ ) ;
equation
re restriccion ;
fo . . f = e = sum ( k , P ( k ) * Datos ( k , ’ Costo ’ ) ) ;
re . . sum ( k , P ( k ) ) = e = D;
Model Despachoideal / A l l / ;
s o l v e Despachoideal usi n l p minimizing f ;
Execute Unload ’ I d e a l . gdx ’ , P , f ;

Un dato interesante en el modelo usando Gams es el costo marginal

asociado a la restricción re el cual define el precio de bolsa. En un modelo
mas simple el precio de bolsa lo define la ultima planta no inflexible
197
despachada, en este caso la hidráulica 2 con $40. Este despacho satisface

la demanda de 200 MW y tiene en cuenta las restricciones de las diferentes
unidades de generación, no obstante, al ubicar este despacho en la red
se nota que las lı́neas 3 − 5, 4 − 5 y 5 − 2 sufrirán una sobrecarga. Para
compensar este hecho se utiliza el siguiente razonamiento:
Los generadores hidráulicos tienen el costo más bajo del sistema por lo
cual deben abastecer la mayor parte de la demanda posible, la capacidad
de las lı́neas limitan esta generación a 50 MW cada una. Los 100 MW
restantes deben ser generador por las plantas térmicas, en donde la planta
térmica 1 es mas eficiente por lo cual debe ser despachada con 80M W . Los
20 MW restantes deben ser generados por la unidad 2. Ası́ el despacho
programado es el siguiente:
Térmica 1 : 80 MW
Térmica 2 : 20 MW
Este nuevo despacho genera costos de reconciliación los cuales se

describen a continuación:
La unidad Hidráulica 1 deja de generar 40 MW los cuales se cancelan

a precio de reconciliación negativa ( 35+40
2
= $37, 5). De la misma forma
la planta hidráulica 2 presenta una reconciliación negativa de 60 MW a
un precio de 40+402
= $40. De otro lado, las plantas térmicas presentan
reconciliación positiva de 80 MW y 20 MW a un costo de $50 por MWh
para la planta térmica 1 y $60 por MWh para la planta 2.
Este problema puede ser solucionado usando Gams como se muestra

a continuación.

$ontext
Ejemplo de despacho programado
$offtext
s e t s k generadores /H1 , H2 , T1 , T2/

n nodos /N1 * N5/
l lineas /L1 * L4 / ;
t a b l e Datos ( k , * ) Datos de o f e r t a
MW Costo Nodo
H1 90 35 3
H2 150 40 4
T1 80 50 1
T2 120 90 2;
t a b l e L i n e a s ( l , * ) Datos de conexion
N1 N2 MW
L1 1 5 50
L2 2 5 50
L3 3 5 50
L4 4 5 50;
parameter Demanda ( n ) demanda en cada nodo /N1 4 0 ,
N2 6 0 ,
N3 0 ,
N4 0 ,
N5 1 0 0 / ;
variable
f funcion o b j e t i v o
Pf ( l ) F l u j o s de p o t e n c i a
Pn ( n ) P o t e n c i a nodal ;
Pf . up ( l ) = L i n e a s ( l , ’MW’ ) ;
Pf . l o ( l )= − L i n e a s ( l , ’MW’ ) ;
positive Variable
P(k ) P o t e n c i a generada ;
P . up ( k ) = Datos ( k , ’MW’ ) ;
equation
b(n ) b a l a n c e de p o t e n c i a nodal
u(n ) f l u j o s de p o t e n c i a ;
fo . . f = e = sum ( k , P ( k ) * Datos ( k , ’ Costo ’ ) ) ;
u ( n ) . . Pn ( n ) = e = sum ( l , Pf ( l ) $ ( L i n e a s ( l , ’ N1’ ) = ord ( n ) ) ) −
sum ( l , Pf ( l ) $ ( L i n e a s ( l , ’ N2’ ) = ord ( n ) ) ) ;
b ( n ) . . Pn ( n ) = e = sum ( k , P ( k ) $ ( Datos ( k , ’ Nodo ’ ) = ord ( n ) ) ) − Demanda ( n ) ;
Model Despachoprogramado / A l l / ;
s o l v e Despachoprogramado u s i n l p minimizing f ;
Execute Unload ’ Despacho . gdx ’ , Pf , P ;

199
10.1.35
Los contratos bilaterales de energı́a son intrumentos para el Contratos
cubrimiento del riesgo asociado con el precio de bolsa2 . Solo es posible bilaterales
realizar contratos entre los generadores y comercializadores del sistema,
para ello se debe presentar la solicitud de registro ante el ASIC con 5
dı́as de anticipación a la fecha en la cual entra en Operación Comercial
el contrato. Además, éste debe contener reglas o procedimientos claros
para determinar hora a hora, las cantidades de energı́a exigibles bajo el
contrato, y el precio respectivo. Se pueden presentar múltiples tipos de
contratos entre los agentes de un sistema, no obstante, existen algunos
contratos estándar a saber:
Pague lo contratado: El comercializador se compromete a pagar toda la

energı́a contratada a una tarifa pactada, independientemente que
ésta energı́a se consuma. Si el comercializador contrató una cantidad
mayor a su demanda efectiva, la diferencia se vende en bolsa. Este es
el único caso en que el comercializador vende energı́a en bolsa. En
la gráfica 10.4
Vende energía en bolsa

Compra energía en bolsa
MWh Contrato
Energía comprada por el contrato
Tiempo
Figura 10.4: Contrato tipo pague lo contratado
2
En especial en sistemas hidrotérmicos en donde la hidrologı́a puede ocasionar
volatilidad en el precio de bolsa.
Pague lo demandado: También denominado pague lo consumido. A

diferencia del anterior, este tipo de contrato define a priori la
cantidad de energı́a contratada, ésta se conoce al calcular la demanda
real del comercializador. El generador se compromete a suministrar
la totalidad de la demanda requerida por el comercializador (no
necesariamente de forma fı́sica). La figura 10.5 muestra este modelo.
MWh
Energía comprada por el contrato = Demanda
Tiempo
Figura 10.5: Contrato tipo pague lo demandado
Pague lo demandado con tope: Tiene la misma filosofı́a que el contrato

pague lo demandado solo que en este caso se pueden definir
restricciones en el mismo. El tope corresponde a una potencia
máxima por encima de la cual el comercializador se puede exponer
a la bolsa. La figura 10.6 muestra con detalle este tipo de contrato.
Todos los contratos son liquidados centralmente por el LAC. En

caso de que un comercializador suscriba varios contratos, estos deben
ser despachados teniendo en cuenta que los contratos tipo pague lo
contratado tienen prioridad sobre los pague lo demandado, es decir,
primero se despachan todos los contratos tipo pague lo contratado,
posteriormente se despachan los contratos tipo pague lo demandado
(con y sin tope), por orden de méritos. Si hay un empate entre dos
contratos, se distribuye la energı́a en forma proporcional al valor (MWh)
del contrato. En el sistema interconectado nacional (SIN) no existe la
figura de contratos fı́sicos, siendo todos ellos de carácter financiero, en
otras palabras, la energı́a despachada en contratos no incide de manera
directa en el despacho programado.
201
Energía comprada en bolsa
Tope
MWh
Energía comprada por el contrato
Tiempo
Figura 10.6: Contrato tipo pague lo demandado
El uso de contratos es un medio de cubrimiento frente al riesgo

en los mercados eléctricos. En caso de presentarse precios bajos en la
bolsa, los generadores obtienen ganancias aún sin ser despachados; los
comercializadores están dispuestos a pagar este sobre-costo para cubrirse
frente al riesgo que representan las alzas de precios. En la tabla 10 se
muestra de forma simplificada este aspecto.
Los contratos bilaterales presentan un tratamiento distinto cuando

se trata de usuarios regulados o usuarios no regulados. Las grandes
cargas del sistema (superiores a 2M W ) se consideran como usuarios
no regulados, éstos pueden comprar energı́a por medio de contratos
bilaterales directos o compra en bolsa. Los usuarios regulados
corresponden a la mayorı́a de cargas del sistema en especial las
residenciales, las cuales como es de esperarse, presentan una demanda
individual inferior a 2M W . Los contratos con destino al mercado regulado
se compran por medio de licitación publica en donde entra nuevamente en
juego en concepto de eficiencia.
Tabla 10.2: Caracterı́sticas del contrato tipo pague lo contratado

Agente Ventajas Desventajas
Comprador
Puede obtener mejores precios Riesgo de exposición a la bolsa
Preve el flujo de caja Pérdidas en caso de una
Ganancias en caso que el precio mala predicción de la demanda
en bolsa sea superior al contrato
Vendedor
Riesgo de exposición a la bolsa

Garantiza ingresos fijos en caso de no ser despachado
Ejemplo 10.2
Mostrar el proceso de liquidación para los siguientes contratos pactados para

un comercializador
• Un contrato tipo pague lo demandado (PLD1) con tope de 40M W

con el generador GT2 por valor de (62$/M W h).
• Un contrato de 30M W tipo pague lo contratado (PLC) por valor de

(30$/M W h) con el generador GT1.
• Un contrato tipo pague lo demandando (PLD2) con tope de 10M W

con el generador GH1 por valor de (50$/M W h).
En la tabla 10 se muestra el proceso de liquidación y los costos

correspondientes para el comercializador si el precio de bolsa es de
45$/M W h y la demanda real del comercializador es de 100M W h.
Es común que el precio de bolsa sea inferior al precio pactado en los

contratos, pero éste es un sobre-costo que algunos comercializadores están
dispuestos a pagar (quizás a falta de análisis más sofisticados) para evitar
el riesgo de exposición a la volatilidad de la bolsa.
203
Tabla 10.3: Caracterı́sticas del contrato tipo pague lo demandado

Agente Ventajas Desventajas
Comprador
Nunca vende exceso en bolsa Imposibilidad de tranzar

Si es sin tope, cubre en bolsa en situación
totalmente el riesgo de la bolsa de precios bajos
Vendedor
En caso de no ser despachado y

presentarse bajos precios en bolsa Expuesto a la demanda
obtiene ganancias del comprador
Tabla 10.4: Liquidación de contratos.

Contrato Valor Potencia Total
CPLC 30 30 900
CPLD1 50 10 500
CPLD2 62 40 2480
Bolsa 45 20 900
10.2 Poder de mercado
Se define poder de mercado como la habilidad de afectar

unilateralmente el precio de la energı́a. Solo se presenta competencia
perfecta cuando todos los agentes tienen asociado un poder de mercado
igual. No obstante, el mercado eléctrico no cumple las condiciones
de competencia perfecta debido a las restricciones técnicas y a la
heterogeneidad en la capacidad de generación de las diferentes centrales.
Existen algunas métricas que permiten definir el poder en el mercado

de determinada firma, una de las más utilizadas es el ı́ndice de Lerner.
Este ı́ndice es también llamado de precio-costo marginal y está dado por
la siguiente relación:
π t − λk
Lk = (10.12)
πt
en donde
πt : Precio de bolsa
λk : Costo marginal del productor k.
Desde luego las centrales mas eficientes presentan un costo marginal

menor lo que significa que tienen mayor poder sobre el mercado.
Los agentes con mayor poder tienen una mayor probabilidad de ser
despachados por lo que ejercen un mayor control sobre los precios. No
obstante el costo final de la energı́a comprada puede depender en gran
medida de las restricciones de la red, las cuales producen reconciliaciones.
10.3 Estrategias de los agentes
El mercado eléctrico es en realidad la superposición de varios

mercados entre los que se destacan la generación de energı́a eléctrica,
el servicio de AGC y control de tensión, la generación de potencia
reactiva y el aseguramiento del margen de seguridad. Esta caracterı́stica
permite a los diferentes agentes tener estrategias de oferta distintas y
complementarias, a saber:
• Algunos generadores -especialmente las centrales hidráulicas-

buscan ser despachadas por la bolsa en el esquema de orden
de méritos, en este sentido, su oferta busca ser cercana al costo
marginal de su operación para asegurar una porción del despacho.
Igualmente, un costo de oferta bajo permite menores egrésos en caso
de una reconciliación negativa.
205
• Otros generadores -en especial los térmicos- buscan ofertar a precios

altos para inducir un precio de bolsa igualmente alto, en caso de que
el nivel de los embalses sea insuficiente para abastecer la demanda.
Ası́ mismo, estos generadores pueden aumentar sus ingresos por
concepto de reconciliaciones positivas en caso de no ser despachados
por orden de méritos.
• Los precios de los contratos influyen en la bolsa de energı́a al
igual que ésta influye en los primeros. La bolsa de energı́a regula
los costos de los contratos pues una diferencia muy alta entre
éstos producirı́a que los comercializadores se volcaran en un solo
tipo de compra (bolsa o contratos). Igualmente, los costos de los
contratos dan señales de más largo plazo para que los generadores
realicen sus ofertas. Es de esperarse que los generadores busquen
ser despachados, reduciendo su oferta, para abastecer su propios
contratos. En caso de nos ser despachados, su ganancia está
garantizada en la medida en que el precio de bolsa sea menor al de
los contratos.
1. Unidad de planeamiento minero energético UPME. Visión del

Mercado Eléctrico Colombiano Bogotá, 2004.
2. Juan C. Morales Ruiz. Silvia Elena Cossio Mesa. Modelo de despacho
económico del sistema eléctrico Colombiano. X Simposio de especialistas
en planeamiento de la operación y expansión. Floranopolis Brasil
2006.
3. Introducción a los conceptos de equilibrio en economı́a. Universidad
Nacional de Colombia. Facultas de ciencias Económicas. Bogotá
2002.
11
Simulación probabilı́stica de costos
Contenido
Simulación probabilistica
Salidas forzadas
Índices de confiabilidad
Simulación de montecarlo.
207
208 CAPÍTULO 11. SIMULACIÓN PROBABILÍSTICA DE COSTOS
L simulación es una importante herramienta en el análisis de los

A
sistemas de generación de energı́a. Al igual que en el flujo de
carga en donde se pueden determinar las variables de estado del
sistema bajo determinadas condiciones reales o hipotéticas, la simulación
probabilı́stica permite encontrar el comportamiento del sistema ante
diferentes condiciones estocásticas.
Como en todos los problemas de simulación el modelo empleado

puede ser tan simple o complejo como se desee. En este caso se usan las
siguientes consideraciones:
• Los costos de operación de las unidades son lineales. (oferta en

bolsa)
• Las plantas son despachadas por orden de méritos.
• La demanda es modelada por una curva de duración de carga.
Modelos más complejos pueden considerar las polı́ticas regulatorias,

contratos y modelos hidrológicos.
11.1 Demanda
La demanda del sistema puede ser modelada de forma cronológica

y detallada por medio de una curva de carga o de forma monótona por
medio de una curva de duración de carga. En la figura 11.1 se muestra una
curva de duración, en donde la demanda es graficada en el eje Y mientras
la duración en horas es registrada en el eje X;
Aunque la curva de carga puede presentar datos importantes sobre

la operación del sistema, es la curva de duración de carga (CDC) la que
permite un modelo probabilı́stico de la demanda. La CDC constituye una
re-oganización cronológica de la curva de carga, en este sentido, puede
determinar la duración en horas de un determinado nivel de demanda.
209
1000 MW
800 MW
600 MW
400 MW
5 12 21 24
5h 7h 9h 3h
Figura 11.1: Curva de duración de carga
Cada punto Xk corresponde al tiempo en el que la demanda del

sistema es menor o igual al valor Yk . Es posible definir una función de
Y en X dada como
t = fCDC (P ) (11.1)
donde el tiempo t está dado en horas mientras que la potencia

demanda P está dada en M W . Normalizando esta relación es
posible interpretar esta función como una distribución acumulada de
probablidad:
T − fCDC (Pk )
P (P ≤ Pk ) = (11.2)
T
donde T es el periodo de tiempo que representa la curva1 .
El modelamiento probabilistico de la demanda mediante la curva de

duración de carga asume un patrón de comportamiento constante, esto
1
La curva de duración de carga se puede construir de forma horaria para un periodo
de 24h, un mes o un año.
supone una aproximación a la operación real ya que este comportamiento

cambia entre dı́as normales y festivos al igual que de un mes a otro.
11.2 Generación
Las plantas son ordenadas ascendentemente (orden de méritos) pues

el sistema debe garantizar criterios de eficiencia. De esta forma, se puede
encontrar una solución única para cada estado en la simulación pues
cuando la generación es superior a la demanda (superavit de energı́a), ésta
es abastecida desde la mas eficiente a la mas costosa. Si no se consideran
las salidas forzadas, la cantidad de energı́a asociada a cada unidad se
puede determinar superponiendo la capacidad de generación sobre la
curva de duración de carga, como se muestra en la figura 11.2. Las salidas
forzadas pueden ser modeladas mediante un indicador denominado tasa
de salidas forzadas o FOR por sus siglas en inglés (Forced Outage Rate).
Éste ı́ndice se determina mediante datos estadı́sticos de los tiempos
promedio en que una central presenta indisponibilidad. Una forma más
exacta de modelar esta indisponibilidad es mediante el ajuste a una
distribución de probabilidad.
1000 MW
800 MW
Generador B 600 MW
400 MW
Generador A
5 12 21 24
5h 7h 9h 3h
Figura 11.2: Modelo de la generación

211
Tabla 11.1: Unidades de generación disponibles

Unidad FOR S [M W ] C [$/M W h]
A 0 600 17,6
B 0,2 400 20,7
C 0,1 200 24,0
En algunos casos el número de estados hace necesario el uso de técnicas

computacionales como la simulación de Montecarlo2 .
Ejemplo 11.1
Una demanda dada por la curva 11.1 debe ser abastecida por tres unidades de
generación caracterizadas por la tabla 11.1, mostrar el proceso de simulación de
costos sin considerar las salidas forzadas
Los generadores son ordenados por méritos desde el más económico

al más costoso, de esta manera se asegura un despacho económico. La
generación de energı́a por parte de cada unidad, sin considerar las salidas
forzadas, está dada por el área bajo la curva superpuesta de la demanda y
la capacidad de generación como se muestra en la figura 11.2:
Para el generador A
SA
EA = fCDC (P ) · dP = 13800M W h (11.3)
0
Para el generador B
SA +SB
EB = fCDC (P ) · dP = 3400M W h (11.4)
SA
2
Ver referecia [3]
Para el generador C
SA +SB +SC
EC = fCDC (P ) · dP = 0M W h (11.5)
SA +SB
11.2.36 Salidas Cuando se considera las salidas de unidades de generación, el

forzadas procedimiento tiene en cuenta cada uno de los estados con su
correspondiente probabilidad 3 .
Cada evento debe ser analizado de forma independiente, en el caso

de N centrales de generación el número de estados es 2N dado que cada
central puede estar en una de dos situaciones operativas: encendido o
apagado. Para una curva de demanda discretizada en T periodos, el
número total de estados ET es:
ET = T · 2N (11.6)
La probabilidad asociada a cada estado puede ser calculada con base

en el ı́ndice de salida forzada de cada planta.
Ejemplo 11.2
Simular el sistema presentado en el ejemplo anterior teniendo en cuenta las

salidas forzadas de las centrales de generación.
La demanda debe ser suministrada teniendo en cuenta criterios

económicos, para ello se ordenan las centrales por orden de méritos;
posteriormente, es necesario definir las posibles condiciones de operación
teniendo en cuenta la tasa de salidas forzadas de cada unidad (FOR). En
3
Esta aproximación no siempre es correcta ya que asume tasa de eventos constantes.
213
este caso, la unidad de generación A tienen una tasa de salida de 0% lo

cual significa disponibilidad total. Los posibles estados de generación son
mostrados en la tabla 11.
Tabla 11.2: Posibles estados de operación

Estado A B C Probabilidad
1 0 0 0 0
2 0 0 1 0
3 0 1 0 0
4 0 1 1 0
5 1 0 0 (0, 2) · (0, 1) = 0, 02
6 1 0 1 (0, 2) · (1 − 0, 1) = 0, 18
7 1 1 0 (1 − 0, 2) · (0, 1) = 0, 08
8 1 1 1 (1 − 0, 2) · (1 − 0, 1) = 0, 72
En este problema en particular, los primeros 4 estados son descartados

pues la unidad generadora A esta disponible todo el tiempo. Los demás
estados tienen un tiempo dado por la probabilidad del mismo. A
continuación se muestra el procedimiento para cada uno de éstos estados:
Estado 5 : Solo el generador A está operando (ver figura 11.3).
1000 MW
800 MW
ENS 600 MW
400 MW
Generador A
5 12 21 24
5h 7h 9h 3h
Figura 11.3: Estado 5

• EA = 13800M W h · (0, 02) = 276M W h

• EB = 0
• EC = 0
• EN S = 3400M W h · (0, 02) = 68M W h (Energı́a no servida)
Estado 6 : Conectados el generador A y el generador C (ver figura 11.4).
1000 MW
ENS 800 MW
C 600 MW
400 MW
Generador A
5 12 21 24
5h 7h 9h 3h
Figura 11.4: Estado 6
• EA = 13800M W h · (0, 18) = 2484M W h

• EB = 0
• EC = 2400M W h · (0, 18) = 432M W h
• EN S = 1000M W h · (0, 18) = 180M W h
Estado 7 Conectados el generador A y el generador B (ver figura 11.5).
• EA = 13800M W h · (0, 08) = 1104M W h

• EB = 3400M W h · (0, 08) = 272M W h
• EC = 0
• EN S = 0
Estado 8 Todos los generadores funcionando (ver figura 11.5).

215
1000 MW
800 MW
Generador B 600 MW
400 MW
Generador A
5 12 21 24
5h 7h 9h 3h
Figura 11.5: Estados 7 y 8
• EA = 13800M W h · (0, 72) = 9936M W h

• EB = 3400M W h · (0, 72) = 2448M W h
• EC = 0
• EN S = 0
Total : Total de energı́a suministrada por cada unidad de generación

• EA = 13800M W h
• EB = 2720M W h
• EC = 432M W h
• EN S = 248M W h
• T otal = EA + EB + EC + EN S = 17200M W h (Total energı́a
demandada)
Costos : Una vez determinada la energı́a total generada por cada unidad
es posible conocer los costos totales de la operación:
• CA = $242880
• CB = $56304
• CC = $10368
• T otal = $309552 (Mas el costo de racionamiento)
11.2.37
La confiabilidad es la habilidad del sistema para cumplir su función.
Confiabilidad
En el caso de un sistema de generación de energı́a, su función es la de
suministrar potencia activa y reactiva ademas de mantener un perfil de
tensión y frecuencia. La confiabilidad está ı́ntimamente ligada con la
seguridad del sistema, esta última determinada mediante un análisis de
contingencias.
Existen múltiples ı́ndices de confiabilidad para el sistema de

generación, a continuación se presentan algunos de los más difundidos
en la literatura especializada:
LOLP : (Loss of load expectation): Es una medida de la probabilidad que

la demanda del sistema exceda la capacidad instalada durante un
periodo de tiempo.
EPNS : (Expected power not supplied): Es el valor esperado de la energı́a
no servida en el sistema.
Un elemento importante en el análisis de los sistemas de generación es

la redundancia en los elementos con el fin de garantizar la confiabilidad.
Aunque la capacidad instalada en generación supere la demanda, esto no
garantiza abastecimiento total de la misma. Por tal razón es importante
tener en cuenta estos ı́ndices de confiabilidad a la hora de planificar y
operar el sistema.
Ejemplo 11.3
Calcular los principales indices de confiabilidad para el ejemplo anterior.
• El total de energı́a no servida es EN S = 248M W h el cual

corresponde a un 1, 44% de la demanda total
217
• El tiempo en horas en que no se tiene disponibilidad suficiente para

abastecer la totalidad de la demanda está dado por los estados 5 y 6:
– Estado 5 : (5h + 7h) · (0, 02) = 0, 24h

– Estado 6: (5h) · (0, 18) = 0, 9h
– Total : 1, 14h
1,14
– LOLP = 24
= 0, 0475
1. Brian Manhire.EPRI monographs on simulation of electric power

production: Book 6: The probabilistic simulation workbook.
2. Zapata Carlos Julio, Confiabilidad en sistemas eléctricos. Universidad

Tecnológica de Pereira. Maestrı́a en Ingenierı́a eléctrica. 2005.
3. Zapata Carlos Julio, Garcés Lina Paola, Gomez Oscar. Modelamiento

de componentes del sistema compuesto generación-transmisión para
estudios de confiabilidad. Revista Scientia et Technica. Año X. Nº 25.
Agosto 2004.

5. Sullivan, Robert Lee. Power system planning. McGraw Hill.

12
Planeamiento de la generación
Contenido
Modelo desterminı́stico
219
220 CAPÍTULO 12. PLANEAMIENTO DE LA GENERACIÓN
E crecimiento de la demanda de energı́a eléctrica requiere un

L
incremento igual en la capacidad de generación del sistema. Para
determinar este crecimiento de forma óptima se recurre a los planes de
expansión de la generación . La selección de esquemas de generación,
especialmente en las plantas térmicas y los sistemas de co-generación,
requiere tener en cuenta los costos de inversión ası́ como la operación y
confiabilidad del sistema.
Tan importante como las metodologı́as de solución a los diferentes

problemas en ingenierı́a, es el adecuado planteamiento de los mismos.
Por tal razón en este capı́tulo se presentan los conceptos básicos sobre el
planeamiento de la generación y su solución mediante un software (en
este caso Gams). La forma en que estos programas dan solución este tipo
de problemas va mas allá de los objetivos de este libro 1
12.1 Modelo estático
El problema consiste en determinar la configuración óptima de

generadores térmicos para abastecer una demanda proyectada de energı́a.
Este problema presenta variables estocásticas, al igual que la casi la
totalidad de problemas en planeamiento; no obstante, se el planteamiento
corresponde a un modelo determinı́stico, en donde la demanda futura es
conocida.
Las variables de decisión son enteras pues corresponden a las unidades

térmicas a instalar. Ası́ Xk es el número de unidades del tipo k a instalarse.
El tipo de unidades factibles de ser seleccionadas esta determinado a priori
por los datos de fabricantes; estos datos corresponden a la capacidad de
generación, indisponibilidad proyectada, costos (tanto de inversión como
de operación), entre otros. Este modelo es denominado estático ya que
la totalidad de la inversión es realizada en el primer periodo aún cuando
pueda ser necesario solo en los últimos periodos.
1
El lector puede consultar las referencias [2] y [3] , entre los métodos más utilizados
en la solución de problemas de programación lineal entera mixta esta Branch and Bound
221
12.1.38 Función
La función objetivo consiste en minimizar la suma de los costos de objetivo
inversión, costos del combustible, operación y mantenimiento. Cada uno
de estos costos son desagregados a continuación:
Costos de inversión (Cik ): Corresponde a los costos de compra del activo

ası́ como su instalación en la central. Si se tiene el costo de los
generadores en el puerto (fuera del paı́s), el costo total debe ser
ajustado por un factor2
N
Ck · Xk (12.1)
k=1
Costos del combustible (Cck ): Son proporcionales a la potencia

promedio generada por la unidad en cada intervalo de tiempo
t. Dependen del tipo de tecnologı́a utilizada.
N
NT Cck · Pkt (12.2)
t=1 k=1
Costos de operación y mantenimiento: Incluyen dos términos: unos fijos

Cfk los cuales dependen de la unidad seleccionada y otros variables
Cvk que dependen de la potencia a generar por dicha unidad.
N NT N
Cfk · Xk + Cvk · Pkt (12.3)
k=1 t=1 k=1
2
Este factor de aproximadamente 2 según la referencia [4]
12.1.39
El sistema se asume despachado a nodo único, por tal razón las
Restricciones
principales restricciones corresponden al balance de potencia nodal:
Capacidad de generación: Cada uno de los generadores seleccionados

tienen una capacidad máxima Mk :
Pkt < Mk · Xk (12.4)
Balance de potencia: La potencia generada debe ser mayor que la

demanda en cada periodo (Dt ).
N
Pkt > Dt (12.5)
k=1
A este modelo se pueden agregar, entre otras restricciones adicionales,

las siguientes:
• Presupuesto máximo de inversión.

• Parámetros de confiabilidad.
• Costos de indisponibilidad.
• Numero máximo de unidades.
Ejemplo 12.1
Se desea realizar la expansión de un pequeño sistema térmico, la demanda

adicional proyecta para los siguientes tres periodos es (20 MW, 30 MW y 35MW);
los parámetros de las posibles unidades se dan en la tabla 12.1.
Se desea hacer un planeamiento de la expansión utilizando máximo

dos unidades del mismo tipo. El modelo escrito en Gams es el siguiente:
223
Tabla 12.1: Parámetros de las unidades de generación disponibles

Unidad Ci Cc Cf Cv M
G1 100 5.3 1.2 3.2 10
G2 180 5.2 1.3 3.1 15
G3 160 6 1.5 3.2 10
G4 130 5.8 1.3 3.2 10

$ t i t l e Planeamiento de l a expansion de l a g e n e r a c i o n
$ontext
Modelo d e t e r m i n i s t i c o para e l
de planeamiento de l a g e n e r a c i o n
s i n t e n e r en cuenta l a
indisponibilidad
$offtext
sets
k generadores / G1 * G4 /
t periodos / t1 * t3 / ;
parameters
Ci ( k ) c o s t o f i j o de ó i n v e r s i n
/ G1 100
G2 180
G3 160
G4 130 /
Cc ( k ) c o s t o s de c om b us ti b le
/ G1 5.3
G2 5.2
G3 6.0
G4 5.8 /
Cf ( k ) c o s t o s f i j o s de o p e r a c i o n
/ G1 1.2
G2 1.3
G3 1.5
G4 1.3 /
Cv ( k ) c o s t o s v a r i a b l e s de o p e r a c i o n
/ G1 3.2
G2 3.1
G3 3.2
G4 3.2 /
M( k ) capacidad de cada unidad
/ G1 10
G2 15
G3 10
G4 10 /
D( t ) demanda para un e s c e n a r i o [MW]
/ t1 20
t2 30
t3 35 /
variable
f f u n c i o n de c o s t o t o t a l
positive variables
fci f u n c i o n de c o s t o s de i n v e r s i o n
fcc f u n c i o n de c o s t o s de c om b u s ti b le
fcf f u n c i o n de c o s t o s f i j o s de o p e r a c i o n
fcv f u n c i o n de c o s t o s v a r i a b l e s
p ( k , t ) p o t e n c i a generada en cada periodo ;
integer variable
x(k ) numero de unidades de cada t i p o ;
x . up ( k ) = 2 ;
equations
qci equacion de c o s t o s de i n v e r s i o n
qcc equacion de c o s t o s de c om bu s ti bl e
qcf equacion de c o s t o s f i j o s de o p e r a c i o n
qcv equacion de c o s t o s v a r i a b l e s
qd ( t ) b a l a n c e de p o t e n c i a
qm( k , t ) capacidad de cada p l a n t a ;
fo . . f = e = f c i + fcc + f c f + fcv ;
qci . . f c i = e = sum ( k , Ci ( k ) * x ( k ) ) ;
qcc . . f c c = e = sum ( k , sum ( t , Cc ( k ) * P ( k , t ) ) ) ;
qcf . . f c f = e = sum ( k , Cf ( k ) * x ( k ) ) ;
qcv . . f c v = e = sum ( k , sum ( t , Cv ( k ) * P ( k , t ) ) ) ;
qd ( t ) . . D( t )= l = sum ( k , P ( k , t ) ) ;
qm( k , t ) . . P ( k , t ) = l = M( k ) * x ( k ) ;
model g e n e r a c i o n / a l l / ;
s o l v e g e n e r a c i o n minimizing f using mip ;

display p . l , x . l

La solución encontrada consiste en dos unidades del tipo G1 y una

unidad del tipo G2. Las potencias a generar pueden ser determinadas por
orden de méritos (Cv + Cc). Una forma poco eficiente de dar solución a
este problema consiste en verificar todas las configuraciones posibles de
generación para abastecer la demanda y evaluar sus costos, como es de
esperarse esta metodologı́a puede resultar inaplicable desde el punto de
vista computacional. Dado que las variables de decisión son discretas, el
problema no puede ser solucionado por técnicas como los multiplicadores
de lagrange. Debe ser usada una metodologı́a de programación lineal
entera mixta. (PLEM).
225
1. H.M.Khodr et all. A Linear Progranning Methodolory for the

optimization of Electric Power-Generation Schemes. En IEEE Trans.
on Power Systems. Vol 17 Nº3. Aug. 2002.
2. Gas Turbina World 1996 Handbook. Pequot.
3. Carlos Julio Zapata. Confiabilidad de sistemas eléctricos. Maestrı́a

en ingenerı́a eléctrica.
Métodos de Aforo.
A
E L registro de los caudales históricos constituyen la base para
determinar la generación de potencia asegurada en un aprovechamiento
hidroeléctrico. No obstante, en algunas ocasiones no se cuenta con esta
información por lo cual es necesario tomarla de forma manual. En este
anexo se presentan las principales formas de aforo para aprovechamientos
de pequeña y mediana escala.
A.1 Llenado de deposito
En arroyos con menos de 20m3 /s es posible desviar el caudal (o parte

de él) hacia un depósito de volumen conocido (V ), el procedimiento
consiste en cronometrar el tiempo (T ) que tarda en llenarse, el caudal
puede ser calculado como:
227
228 ANEXOS A. MÉTODOS DE AFORO.
Figura A.1: Método de llenado
V
Q= (A.1)
T
A.2 Método del flotador
Puede producir errores al rededor del 20%, consiste utilizar un flotador

(que puede ser una botella a medio llenar) para medir la velocidad
media del rio; se utiliza en los casos en que el nivel del agua tiene una
profundidad para permitir un recorrido sin obstrucciones del flotador y
presenta un velocidad relativamente baja. La longitud del recorrido (x)
debe ser superior a 10m.
El flotador debe ser lanzado antes de la zona demarcada y cronometrar

su recorrido en la parte media del lecho. Si en algún momento el flotador
toca el lecho del rı́o la medida debe ser descartada. Después de tomar
varias mediciones se puede obtener la velocidad superficial (Vsup )como:
x
Vsup = (A.2)
T
229
m
10
x=
b
y1
Figura A.2: Método del flotador
La velocidad promedio, con la cual se calcula el caudal promedio, debe

ser modificada por un factor K como se describe a continuación.
Se calcula el área (S)correspondiente a los trapecios resultantes de la

medición de las distintas profundidades, y los perı́metros Pi de cada uno.
(yi + yi−1 ) · b
Si = n
(A.3)
2
!
2
b b
Pi = yi−1 + yi + + + (yi − yi−1 )2 (A.4)
n n
Si la diferencia entre la relación área-perı́metro máximo y mı́nimo es

menor a 0.1 se tiene:

Si Si
− ≤ 0.1 (A.5)
Pi M ax Pi M in

1 Si
K = A · ln · +B (A.6)
n Pi
en donde A y B dependen del material del canal cono se muestra en la

tabal A.1:
Tabla A.1: Constantes A y B

Material A B
Canal de barro 0.0905 0.782
Canal pedregoso 0.0362 0.847
Canal de concreto 0.0150 0.893
El caudal promedio viene dado por:
Qprom = Vsup · K · ST (A.7)
en donde ST es el área total:
n
ST = Si (A.8)
i=1
Si la relación dada por la inecuación A.5 no se cumple el valor del

caudal promedio debe ser calculado como:
n
Vsup Si
Qprom = · A · ln +B · Si (A.9)
n i=1
Pi
Una forma más exacta de obtener la velocidad del rı́o consiste en

utilizar un elemento que permite medir la velocidad media del rı́o por
medio de señales acústicas proporcionales a la velocidad con que fluye el
231
agua en el rı́o. Este elemento es llamado molinete. Su forma estructural

permite tomar medidas a varias profundidades del rı́o obteniéndose
medidas más exactas (Superiores al 3% de error. Al igual que el método
del flotador, se requiere conocer el area transverzal del rı́o.
Molinete
Peso
Figura A.3: Esquema de un molinete
A.3 Método del vertedero
En arroyos con caudales superiores a 25 litros por segundo es posible

construir una pequeña presa en madera que permita una acertada
medición del caudal. Aplicando el teorema de Bernulli en los puntos (1) y
(2) de la gráfica A.5 se tiene
P1 v12 P2 v22
+ y1 + = + y2 + (A.10)
γ 2·g γ 2·g
si se aproxima v1 = 0 y referencia en el punto 2 se obtiene un forma

general de la velocidad para cualquier valor de y entre los puntos (2) y (3).
Figura A.4: Método del vertedero
"
v2 = 2·g·y (A.11)
El caudal puede ser calculado como:

Q= v · dA (A.12)
En el caso de un vertedero rectangular:
" H
1
Q= v · L · dy = 2 · g · L y 2 · dy (A.13)
0
233
(1)
H (2)
(3)
Figura A.5: Esquema del vertedero
De donde se obtiene la ecuación de Francis:
2" 3
Q= 2·g ·L·H2 (A.14)
3
Una aproximación mas real del fenómeno debe tener en cuenta las
constracciones en los bordes del vertedero por lo cual la equación de
Francis toma la siguiente forma:
3
Q = 1.84 · L · H 2 (A.15)
El valor de L debe tomarse de acuerdo a la siguiente restricción:
1 2
b≤L≤ b (A.16)
2 3
En el caso en que el rı́o presenta un ancho reducido es necesario un

vertedero triangular en donde se utiliza la ecuación de Thomsom:
Q = 1.4H 5/2 (A.17)

o
90
Figura A.6: Esquema del vertedero triangular

Formulación de modelos en Gams.
B
G AMSes un programa de computador ampliamente utilizado para dar
solución a complejos problemas de optimización matemática. Su nombre
corresponde a las siglas de General Algebraic Modeling System; entre sus
principales ventajas se destacan:
• Escalabilidad: El código fuente que representa el modelo es

independiente del tamaño del problema, esto significa que un mismo
problema de 10 o 100 variables presenta el mismo código fuente
aumentando solo la cantidad de datos de entrada y salida.
• Uso de solvers: Gams invoca de forma eficiente diferentes solvers

de tipo comercial (minos, cplex etc), sin modificar el código
correspondiente, de esta forma el usuario solo debe modelar el
problema y elegir el solver mas adecuado para su solución1 .
1
Esta caracterı́stica implica una ventaja en el modelamiento del problema, sin
embargo, en algunas situaciones es posible encontrar mejores resultados desarrollando
235
236 ANEXOS B. FORMULACIÓN DE MODELOS EN GAMS.
• Naturalidad en la sintaxis: la estructura de programación de Gams

está diseñada especı́ficamente para problemas de optimización por
lo cual es relativamente simple generar el modelo.
• Portabilidad: el programa está disponible en versiones de prueba

para diferentes sistemas operativos entre ellos Unix y Windows.
En esta sección se mostrarán las principales caracterı́sticas del

programa ası́ como una breve descripción de la sintaxis. Se utilizará una
metodologı́a basada en ejemplos aplicados al campo de los sistemas de
generación de energı́a.
B.1 programación básica
Se parte de un problema simple de despacho económico de plantas

térmicas: Es necesario despachar 0.95M W en un grupo de plantas
térmicas con costo de combustible dados en forma cuadrática, los datos
de cada generador son mostrados en la tabla B.1:
Tabla B.1: Datos de generación.

a
Generador 2
b Max Min
1 108 1166.9 2 0.5
2 176 1033.3 1.5 0.375
El problema uninodal presenta una función objetivo no lineal y un

conjunto de restricciones lineales de igualdad y desigualdad. La forma
más simple de dar solución a este problema es mediante la técnica de
multiplicadores de Lagrange, sin embargo, el uso del computador permite
un modelo más general en el que es posible modificar las diferentes
restricciones y determinar la sensibilidad de la función objetivo frente a
los parámetros.
un solver particular al problema (por ejemplo usando algoritmos evolutivos
especializados)
237
Inicialmente se nombra el modelo para lo cual se utiliza el comando

$title.

$ t i t l e Despacho economico uninodal .
$ontext
Problema de despacho economico de
p l a n t a s t e r m i c a s b a j o un esquema de nodo
unico .
$offtext

La estructura $ontext y $of f text permite generar comentarios que

facilitan la identificación del modelo cuando se almacena en el archivo. Es
uso de comentarios de este tipo es especialmente importante en los casos
en que se disponga de varios modelos con caracterı́sticas similares pues
permite su posterior identificación y administración.
El comando set es utilizado para indexar conjuntos de variables de

caracterı́sticas similares, en este caso el número de generadores; En general
la instrucción set permite designar los diferentes sub-ı́ndices del modelo.
la sintaxis en los comandos básicos consta de tres partes: inicialmente se
escribe el nombre del comando (en este caso set) y el nombre de la variable
o ecuación seguido por un comentario que indica la función de la misma,
finalmente, se escribe el valor correspondiente. es necesario ; para finalizar
cada orden.

s e t k i n d i c e de l a s p l a n t a s t e r m i c a s / T1 * T2 / ;

En el caso del comando set el valor debe ir encerrado entre dos rayas
oblicuas //; dicho valor puede ser un conjunto de nombres (por ejemplo
/termocentro, termocartagena/) o una serie como en este caso (/T 1 ∗ T 2/);
el uso de este tipo de sintaxis facilita escribir los modelos, por ejemplo,
en caso de tener cinco generadores el dato correspondiente puede ser
descrito como /T 1, T 2, T 3, T 4, T 5/ o simplemente mediante la expresión
simplificada /T 1 ∗ T 5/.
Los datos almacenados tienen carácter de texto por lo que no se pueden

realizar operaciones matemáticas directas con estos indices aún cuando
sean numéricos (en este caso, cada valor numérico es también almacenado
como un caractér)
Los datos tabulados se almacenan mediante el comando table, su

sintaxis es análoga al comando set ya descrito, en este caso el signo ∗
es utilizado para evitar declarar otro set con los valores de las columnas
(/a2, b, P max, P min/).

t a b l e gen ( k , * ) datos de l o s generadores
a2 b Pmax Pmin
T1 108 1166.9 2 0.5
T2 176 1033.3 1.5 0.375 ;

El comando parameter funciona de forma similar a table, se utiliza

principalmente para almacenar vectores (/1, 5, 108/) o datos únicos como
en este caso /0.95/.

parameter
d demanda d e l s i s t e m a / 0 . 9 5 / ;

Las variables son declaradas mediante variable cuando son irrestrictas

o positive variable cuando solo toman valores mayores a cero; éstas
corresponden a las variables de decisión del problema. puede tratarse
de una variable individual (x) o un conjunto de variables indexadas
(x(k),x(k, j)). los lı́mites superior e inferior (si existen) ası́ como el valor
inicial pueden ser intuitivamente declarados, para ilustrarlo se utilizará
una variable P g que representa la potencia generada por cada una de las
unidades térmicas. Los valores mı́nimo y máximo se toman de la tabla
gen para hacer el modelo más general y facilitar la entrada de datos. A
manera de ejemplo se define una condición inicial dada por la potencia
media generada por parte de cada unidad.

variables
costo c o s t o t o t a l de l a o p e r a c i o n ;
positive variables
Pg ( k ) p o t e n c i a generada por cada t e r m i c a ;
Pg . l o ( k ) = gen ( k , ’ Pmin ’ ) ;
Pg . up ( k ) = gen ( k , ’ Pmax ’ ) ;
Pg . l ( k ) = ( Pg . l o ( k ) + Pg . up ( k ) ) / 2 ;

Solo es necesario inicializar las variables en los casos en los cuales la

convergencia puede ser problemática, esto se presenta especialmente en
los problemas de programación no lineal y no lineal entera;
Las ecuaciones son declaradas en dos partes, inicialmente se define

un nombre y el comentario correspondiente; posteriormente se define
239
la ecuación precedida por dos puntos (..), las igualdades se declaran

mediante (= e =) mientras que las desigualdades son declaradas mediante
los sı́mbolos (= l =) o (= g =) si se trata de menor o mayor
respectivamente.

equations
re r e s t r i c c i o n de igualdad ;
fo . . c o s t o = e = sum ( k , ( gen ( k , ’ a2 ’ ) * Pg ( k ) / 2 + gen ( k , ’ b ’ ) ) * Pg ( k ) ) ;
re . . sum ( k , Pg ( k ) ) = e = d ;

Los comandos set, variable, equation, parameter permite su

declaración en singular o en plural según las necesidades.
(sets, variables, equations, parameters) Igualmente, el lenguaje no
distingue entre mayúsculas o minúsculas por los que es importante
seguir una nomenclatura propia para evitar confusiones.
El modelo debe ser nombrado (mediante el comando model) y

solucionado según el caso (minimizing o maximizing). La sentencia /all/
indica que el modelo utiliza todas las restricciones dadas. En caso de
tener varios modelos simplemente se enumeran las restricciones para cada
modelo.
Para dar solución al problema se usa el comando solve acompañado del

nombre del modelo y el tipo de optimización. Los tipos de optimización
más utilizados son los siguientes:
• lp : Programación lineal.
• nlp: Programación no lineal.
• mip: Programación lineal entera.
Aunque el programa reporta los resultados, estos se pueden definir

explı́citamente mediante display.

solve DespachoEconomico minimizing c o s t o using nlp ;
d i s p l a y Pg . l , c o s t o . l ;

B.2 Solvers y modelos disponibles.
GAM S es un programa de modelamiento matemático. Esto significa

que la solución al modelo debe ser encontrada por parte de otros
programas denominados solvers.
La interfaz gráfica del programa permite elegir entre varios solvers

para dar solución a un modelo en particular. Para elegir un solver se
hace clic en la pestaña < f ile > < options >. Posteriormente se elige
la pestaña solvers. Aparecerá una tabla con los solver disponibles y los
tipos de problemas que pueden resolver.
El comando option permite, entre otras cosas, definir el solver a utilizar;

para ello se establece el tipo de optimización y el solver correspondiente,
por ejemplo, si se desea utiliza el programa minos para un problema de
programación no lineal se procede de la siguiente forma:

o p t i o n nlp = minos55

Otra utilidad importante en el entorno grafico de GAMS es la librerı́a

de modelos. Ésta cuenta con una colección de problemas clásicos
modelados y solucionados directamente. Para acceder a la librerı́a de
modelos se hace clic en la pestaña < f ile > < M odel library >
Open gams model library.
Entre los modelos disponibles se encuentran:
• Despacho económico de plantas térmicas incluyendo las pérdidas de

transmisión. (DISPATCH).
• Modelo turco de planeamiento de la generación. (TURKPOW).
• Diseño optimo de turbogeneradores. (BEARING).
• Enganche de unidades térmicas. (FUEL).
• Despacho hidrotérmico. (HYDRO).

241
B.3 manejo de ciclos y estructuras

condicionales
En esta sección se hablará sobre el manejo de dos tipos básicos de

estructuras: las estructuras condicionales y los ciclos iterativos.
Los ciclos iterativos solo requieren ser utilizados en los casos en que no
sea posible usar el paralelismo implı́cito de las declaraciones en GAMS.
Por ejemplo, si se tiene un conjunto de restricciones con caracterı́sticas
similares, solo basta definirla con un indice (set) apropiado.
B.3.40
Las estructuras condicionales permiten tomar decisiones basado en estructura
la evaluación de alguna condición. De otra parte, los ciclos iterativos condicional
permiten solucionar cientos o miles de problemas de optimización que
presenten una estructura similar o que se interrelacionen temporalmente
por medio de alguna función conocida.
La forma más simple de definir un condicional es mediante $. Si por

ejemplo en un conjunto k se desea que solo las primeras 10 variables x(k)
presenten un lı́mite superior de 8, el modelo debe ser descrito ası́:

x . up ( k ) $ ( ord ( k ) < 1 0 ) = 8 ;

Una forma más elaborada de definir un codicinal es mediante la

estructura if . Ésta se utiliza cuando se tienen varias instrucciones que
deben ser ejecutadas cuando se cumple una condición, por ejemplo, para
el caso anterior:

loop ( k ,
i f ( ord ( k) <10 ,
x . up ( k ) = 8 ;
);
);

En este caso se utilizó una estructura iterativa para habilitar el

comando if . Ambos casos cumplen la misma función, no obstante, el
primero presenta mayor simplicidad por lo cual es la más aconsejable.
B.3.41 los ciclos iterativos pueden ser definidos por un comando f or. Para
estructura for mostrar el uso de este comanda se utilizará el ejemplo del despacho de
las plantas térmicas, pero esta vez se desea encontrar el despacho óptimo
para las 24 horas del dı́a de tal forma que la demanda aumente en 0.01 en
cada periodo.

$ T i t l e Despacho economico uninodal .
$OnText
Problema de despacho economico de
p l a n t a s t e r m i c a s b a j o un esquema de nodo
unico para l a s 2 4 horas d e l dia
$OffText
s e t s k i n d i c e de l a s p l a n t a s t e r m i c a s / T1 * T2 / ;
t a b l e gen ( k , * ) datos de l o s generadores

a2 b Pmax Pmin
T1 108 1166.9 2 0.5
T2 176 1033.3 1.5 0.375 ;
parameter d demanda del sistema / 0 . 9 5 / ;
variables
costo c o s t o t o t a l de l a o p e r a c i o n ;
positive variables
Pg ( k ) p o t e n c i a generada por cada t e r m i c a ;
Pg . l o ( k ) = gen ( k , ’ Pmin ’ ) ;
Pg . up ( k ) = gen ( k , ’ Pmax ’ ) ;
equations
re r e s t r i c c i o n de igualdad ;
fo . . c o s t o = e = sum ( k , ( gen ( k , ’ a2 ’ ) * Pg ( k ) / 2 + gen ( k , ’ b ’ ) ) * Pg ( k ) ) ;
re . . sum ( k , Pg ( k ) ) = e = d ;
s c a l a r t contador de i t e r a c i o n e s ;
for ( t = 1 to 2 4 ,
d = d + 0.01;
solve DespachoEconomico minimizing c o s t o using nlp ;
d i s p l a y Pg . l , c o s t o . l ;
);

243
Como se puede observar, la estructura es intuitiva y el modelo base es

el mismo; la diferencia radica en el momento de dar solución al problema.
Se define un escalar t el cual hará las veces de contador de iteraciones,
en cada iteración se aumenta la demanda en 0.01M W , para luego dar
solución al problema y finalmente mostrar los resultados.
B.3.42 Los
Los vectores definidos con el comando set son normalmente definidos operadores ord
en forma de texto aún cuando se trate de números naturales, esto significa y card
que en la definición

set i l o s numeros n a t u r a l e s / 1 * 1 0 / ;

la variable i toma los valores desde el uno hasta el diez pero en forma
de caracter; en otras palabras no es posible definir ningún condicional
de la forma $(i < 5). Para hacer esto es necesario usar el comando ord.
Este comando indica la posición en el orden de definición del vector. Por
ejemplo en la siguiente definición

set i p l a n t a s h i d r a u l i c a s / s a n c a r l o s , guavio , calima / ;

el comando ord(i) dará como resultado el número 2. Ası́, es posible

definir el condicional $(ord(i) < 5).
En el momento de definir una estructura condicional puede ser

importante determina el número de elementos de un vector particular,
para ello se utiliza el comando card. Por ejemplo, en el problema de
despacho hidrotérmico puede ser definida la condición final de volumen
de la siguiente forma:

V( t ) $ ( ord ( t ) = card ( t ) ) = 0 ;

B.3.43 Una forma igualmente simple de definir un proceso iterativo es

estructura loop mediante la estructura loop.
A diferencia de la estructura f or, ésta requiere de un vector definido

con antelación. Para ilustrar el uso de esta estructura se presenta un
problema de despacho hidrotérmico con una única central hidráulica
en donde el volumen en el periodo siguiente depende del mismo en el
periodo anterior y de las afluencias, turbinamientos y vertimientos. En
este caso, la restricción puede ser representada como:

loop ( t ,
V( t + 1 ) = V( t ) + N * (A( t ) − Q( t ) − S ( t ) ) ;
);

Se requiere que la variable tiempo sea definida inicialmente mediante

el comando set. No es necesario indicar el inicio o la condición final el del
ciclo pues ésta esta implı́cita en la definición del conjunto. No obstante,
es posible adicionar un condicional a la estructura, por ejemplo, para
asegurar una coherencia estructural en la metodologı́a es necesario indicar
que no se incluye el último periodo:

loop ( t $ ( ord ( t )< c a r t ( t ) −1) ,
V( t + 1 ) = V( t ) + N * (A( t ) − Q( t ) − S ( t ) ) ;
);

Otro ejemplo del uso de esta estructura es el programa para encontrar

la suma de los primeros 100 números naturales:

set k contador de i t e r a c i o n e s / 1 * 1 0 0 /
s c a l a r N suma de todos l o s numeros ;
N = 0;
loop ( k ,
N = N + ord ( k ) ;
);
display N

245
B.4 Manejo de archivos
En problemas de gran tamaño puede ser útil separar modelos y datos

en archivos independientes, para esto se utiliza el comando $include como
se muestra a continuación:

$ i n c l u d e MODELO1

M ODELO1 representa un modelo realizado anteriormente el cual

es llamado desde el archivo actual. Mediante esta metodologı́a es
posible separar al modelo del archivo que contiene los datos entrada del
problema. Ası́ mismo, existen problemas en donde varios modelos de
optimización trabajan coordinadamente para dar una solución global, por
ejemplo en un algoritmo constructivo de despacho hidrotérmico pueden
ser utilizados modelos de despacho de plantas térmicas y flujo de carga.
Otra herramienta importante es el manejo de archivos de salida. Es

posible almacenar los resultados de un modelo de optimización en un
archivo de texto independiente. A continuación se muestra la declaración
de un archivo llamado /Resultados.txt/:

f i l e Out / R e s u l t a d o s . t x t / ;

Este archivo se almacena en el mismo directorio del modelo (gamsdir)

y puede ser visualizado por cualquier editor de texto. Una vez declarado
el archivo, es posible almacenar datos invocando el alias definido (Out).

put Out ;
Out . nd = 0 ;

La secuencia de comandos Out.nd permite definir el número de

posiciones decimales de cada uno de los resultados obtenidos. En este
caso se redondea a los valores enteros.
Para almacenar los resultados de una variable x arbitraria se usa el

comando put. Igualmente, es posible almacenar texto. La lı́nea inclinada
/ permite saltar a la siguiente lı́nea en el archivo, de no usar esta lı́nea
los resultados se almacenan por filas dificultando la lectura y análisis del

archivo.

put ” Resultado : ” x . l //;

B.4.44 Archivos Una forma relativamente simple de mostrar los resultados o los datos
gdx del sistema es mediante el uso de archivos gdx (GAMS Data Exchange).
Este tipo de archivo funciona de forma análoga al comando display:

Execute Unload ’ R e s u l t a d o s . gdx ’ , x ;

En donde Resultados.gdx es el nombre del archivo en donde se

almacenarán los resultados de las variables x e y. Este archivo es
almacenado en la carpeta gamsdir por lo que solo basta abrirlo (desde la
interfaz gráfica de Gams) para ver de forma ordenada los resultados.
Índice
ARIMA, 161 Reconciliaciones, 193

Reservorio, 42
Bellman, 168 Restricciones, 108, 132, 180
Competencia, 186 Simulación de Montecarlo, 162,
Confiabilidad, 216 207, 211
Despacho hidrotérmico, 160 Trubinas, 68
Desregularización, 187
Economı́a, 16, 184
Francis, 72
Función objetivo, 108
Heurı́sticas, 119
Investigación operativa, 108
Jepirachi, 15
Kaplan, 75
Karush-Kuhn-Tucker, 112
Multiplicador de Lagrange, 173
OLADE, 36
Oligopolio, 188
Optimización no lineal, 109
PCH, 37
Pelton, 68
Plan de expansión, 220
Programación lineal, 116
247

Libro Generacion

Cargado por

Copyright:

Formatos disponibles

Libro Generacion

Cargado por

Información del documento

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

Libro Generacion

Cargado por

Copyright:

Formatos disponibles

Sistemas de Generación de

Programa de Ingenierı́a Eléctrica

Universidad Tecnológica de Pereira.

El propósito fundamental del presente libro, es introducir conceptos

Una de las caracterı́sticas novedosas frente a la literatura clásica,

En la primera parte se describe de forma general el funcionamiento

generación desde el punto de vista eléctrico.

En la segunda parte se presenta la operación economica de los sistemas

Una de las caracterı́sticas únicas de este texto, es el uso de una

Mucho del material aquı́ presentado esta sustentado por trabajos

I Aspectos técnicos de los sistemas de generación 1

1 El sistema eléctrico colombiano. 3

1.0 Estructura del sector eléctrico . . . . . . . . . . . . . . . . . . 5

1.1.1 Regulación y planeamiento . . . . . . . . . . . . . . . 6

1.1.2 Operación y administración . . . . . . . . . . . . . . . 9

1.1.3 Órganos consultores y asesores . . . . . . . . . . . . . 10

1.1.4 Fondos de apoyo . . . . . . . . . . . . . . . . . . . . . 11

1.1 El sistema colombiano . . . . . . . . . . . . . . . . . . . . . . 12

1.2.5 Centrales hidroeléctricas . . . . . . . . . . . . . . . . . 12

1.2.6 Plantas térmicas . . . . . . . . . . . . . . . . . . . . . . 14

1.2.7 Energı́a eólica . . . . . . . . . . . . . . . . . . . . . . . 15

1.2 Transacciones internacionales . . . . . . . . . . . . . . . . . . 15

1.3 La generación de energı́a como variable macroeconomica . . 15

2.0 Funcionamiento de una central térmica . . . . . . . . . . . . 21

2.1 Partes constructivas . . . . . . . . . . . . . . . . . . . . . . . . 22

2.2 Principales ciclos termodinámicos . . . . . . . . . . . . . . . 24

2.3.8 Ciclo de Carnot . . . . . . . . . . . . . . . . . . . . . . 25

2.3.9 Ciclo Rankine . . . . . . . . . . . . . . . . . . . . . . . 27

2.3 Plantas de ciclo combinado. . . . . . . . . . . . . . . . . . . . 29

3.0 Clasiﬁcación de las centrales hidroeléctricas . . . . . . . . . . 36

3.1.10 Central de pasada (a ﬁlo de agua) . . . . . . . . . . . 36

3.1.11 Central de embalse . . . . . . . . . . . . . . . . . . . . 42

3.1 Aspectos técnicos de las centrales de embalse . . . . . . . . . 45

4.1.13 Algunas deﬁnciones . . . . . . . . . . . . . . . . . . . 52

4.1 Estudio de variables promedio . . . . . . . . . . . . . . . . . 54

4.2 Estudio de variables extremas . . . . . . . . . . . . . . . . . . 57

5.0 Caracterı́sticas generales . . . . . . . . . . . . . . . . . . . . . 66

5.1 Tipos de turbinas hidráulicas . . . . . . . . . . . . . . . . . . 67

5.2.14 Turbina pelton . . . . . . . . . . . . . . . . . . . . . . 68

5.2.15 Turbina francis . . . . . . . . . . . . . . . . . . . . . . 72

5.2.16 Turbina kaplan . . . . . . . . . . . . . . . . . . . . . . 75

5.2 Elección de la turbina en una central hidráulica . . . . . . . . 77

5.3 Control de velocidad (Gobernador) . . . . . . . . . . . . . . . 79

6.0 Curva de capacidad de un generador. . . . . . . . . . . . . . 84

6.1 Control de tensión. . . . . . . . . . . . . . . . . . . . . . . . . 88

6.2 Control automático de la generación (AGC) . . . . . . . . . . 90

6.3.17 Control primario de frecuencia . . . . . . . . . . . . . 92

6.3.18 Control secundario . . . . . . . . . . . . . . . . . . . 94

6.3.19 Control Terciario . . . . . . . . . . . . . . . . . . . . . 101

6.3.20 Reserva rodante . . . . . . . . . . . . . . . . . . . . . 102

6.3 Lecturas complementarias . . . . . . . . . . . . . . . . . . . . 103

II Aspectos económicos de los sistemas de generación 105

7 Introducción a la Optimización 107

7.0 Multiplicadores de Lagrange . . . . . . . . . . . . . . . . . . 109

7.1.21 Interpretación económica . . . . . . . . . . . . . . . . 111

7.1 Programación dinámica . . . . . . . . . . . . . . . . . . . . . 112

7.2.22 Algoritmo de solución . . . . . . . . . . . . . . . . . . 113

7.2 Programación lineal . . . . . . . . . . . . . . . . . . . . . . . . 116

7.3 Algoritmos de búsqueda aleatoria . . . . . . . . . . . . . . . . 118

8 Operación de sistemas térmicos. 121