Diseño de Energia Eolica en La Guajira

ANLISIS DE GENERACIN COMPLEMENTARIA ENTRE
ENERGA HIDRULICA Y ELICA

CASO: GENERACIN ISAGEN PROYECTOS ELICOS EN LA
GUAJIRA COLOMBIANA

ARMANDO JOS EALO OTERO
Ingeniero Civil

Trabajo Dirigido de Grado como requisito parcial para optar al ttulo de Magister en
Ingeniera - Recursos Hidrulicos

Maestra en Ingeniera - Recursos Hidrulicos
Escuela de Geociencias y Medio Ambiente
Facultad de Minas

Medelln - Antioquia
2011

ANLISIS DE GENERACIN COMPLEMENTARIA ENTRE
ENERGA HIDRULICA Y ELICA
CASO: GENERACIN ISAGEN PROYECTOS ELICOS EN LA
GUAJIRA COLOMBIANA

ARMANDO JOS EALO OTERO
Ingeniero Civil

DIRECTOR:
Germn Poveda Jaramillo - I.C. MsC. PhD.

CO - DIRECTOR:
Oscar Mesa Snchez - I.C. MsC. PhD.

Maestra en Ingeniera - Recursos Hidrulicos
Escuela de Geociencias y Medio Ambiente
Facultad de Minas
Universidad Nacional de Colombia Sede Medelln
Medelln Antioquia
2011
Anlisis de generacin complementaria entre energa
hidrulica y elica. Caso: generacin ISAGEN proyectos
elicos en La Guajira colombiana

Armando Jos Ealo Otero
Ingeniero Civil

iii

RESUMEN

El anlisis de complementariedad de los recursos elico e hdrico parte de la justificacin
de la diversificacin de matriz energtica de Colombia, que respalde a la energa
hidrulica, la cual se ve afectada directamente por los eventos macroclimticos como el
ENSO y por el cambio climtico antropognico.

A partir de informacin diaria y mensual de series de velocidad del viento y de caudales,
se analiz el recurso viento y energticamente la complementariedad. Los datos
utilizados, a diferentes alturas de medicin, comprenden varias estaciones del IDEAM e
ISAGEN ubicadas en la Costa Caribe colombiana. Para el caso de La Guajira, el anlisis
de la estacin Puerto Bolvar arroj un ciclo diurno definido de la velocidad del viento con
un mximo de 8,7 m s
-1
entre las 15:00 y 16:00 horas y un mnimo de 4,5 m s
-1
. El ciclo
anual reporta mximos de 7,2 m s
-1
en el mes de marzo y en julio con 8,2 m s
-1
, cuando se
presenta el verano en el centro del pas, y mnimos de 5,0 m s
-1
durante el mes de
octubre. Se verific el aumento de la velocidad con la altura entre 20 m y 60 m, con un
coeficiente de Hellman de 0,14. Se determin la direccin Este como la predominante,
asociada al patrn de circulacin de los vientos alisios. Se analiz la correlacin a nivel
trimestral de los ndices del ENSO y la Oscilacin del Atlntico Norte con la velocidad.
Asociado al desplazamiento de la ZCIT, a nivel mensual se identific mayor
complementariedad del viento con caudales del norte y centro del pas, caso contrario que
con los ros del sur como el Pata, Cauca y Luisa. El anlisis de correlacin de la
velocidad del viento con los caudales afluentes a los Proyectos y Centrales
Hidroelctricas de ISAGEN reportan el mayor coeficiente con ro Sogamoso. Con la
velocidad del viento y las caractersticas del aerogenerador Nordex 60/1300 se calcul la
produccin energtica comparndola aisladamente con la generacin histrica de las
centrales de ISAGEN sin encontrar complementariedad en el ciclo promedio anual. La
energa media anual de cada Proyecto se estim a partir de un modelo bajo la RNO y se
calcul el tamao de parques elicos que igualen dichas las energas medias anuales,
requiriendo un parque de 140 MW para igualar Amoy y 1.040 MW para Sogamoso, lo
que indica aumento en capacidades instaladas. El anlisis integrado realizado con
Sogamoso el cual posee embalse regulador, sugiere un aumento en energa firme y
viabiliza financieramente si se opera integradamente con la Central, un Parque de 200
MW.


Ingeniero Civil

iv

AGRADECIMIENTOS

A ISAGEN S.A. E.S.P, por brindar la informacin y los datos necesarios para adelantar
este trabajo investigativo.

A los profesores Germn Poveda y Oscar Mesa por su asesora y tutora durante la
elaboracin y el fundamento terico de ste documento.

Al Ingeniero William ngel, por su orientacin tcnica en el desarrollo de las simulaciones
energticas.

A los ingenieros Luis Posada y Javier Mndez por permitirme participar en los proyectos
de parques elicos que actualmente desarrolla ISAGEN.

DEDICATORIA:

A Margoth y Armando, los pilares de mi ingeniera.


Ingeniero Civil

v

CONTENIDO

AGRADECIMIENTOS .................................................................................................................... iv
1 INTRODUCCIN ..................................................................................................................13
2 OBJETIVOS ..........................................................................................................................17
2.1 OBJETIVO GENERAL ...................................................................................................17
2.2 OBJETIVOS ESPECFICOS ..........................................................................................17
3 JUSTIFICACIN ...................................................................................................................18
4 ESTADO DEL ARTE .............................................................................................................24
4.1 Del diagnstico y caracterizacin del potencial elico de la zona de estudio...................24
4.2 Costos Parques Elicos Internacionales. Regulacin y Normatividad
Colombiana ...............................................................................................................................25
4.3 Del Protocolo de Kioto y los Mecanismos de Desarrollo Limpio - MDL ...........................28
4.4 Complementariedad del recurso Hdrico y Elico ...........................................................29
5 ANLISIS DEL RECURSO ELICO DE LA GUAJIRA COLOMBIANA ..................................36
5.1 Impacto del Evento ENSO (El Nio - Oscilacin del Sur) sobre los vientos de La
Guajira 47
5.1.1 Variacin del ciclo anual de la velocidad el viento ante el Evento ENSO .................49
5.1.2 Ciclo anual de la temperatura de la zona y su relacin la velocidad el viento
y el Evento ENSO .................................................................................................................51
5.1.3 Ciclo anual del ciclo diurno de la velocidad del viento .............................................53
6 COMPLEMENTARIEDAD ENTRE LA HIDROLOGA DE VARIAS REGIONES DE
COLOMBIA CON LA VELOCIDAD DEL VIENTO DE LA GUAJIRA. ..............................................56
7 MODELACIN Y ANLISIS DE GENERACIN DE ENERGA DE PARQUES
ELICOS GENRICOS EN LA GUAJIRA .....................................................................................61
7.1 Aerogeneradores Nordex N60/1.300 ..............................................................................61
8 ANLISIS DE GENERACIN COMPLEMENTARIA HIDRULICA - ELICA ........................65
8.1 CICLO ANUAL DE LAS SERIES DE CAUDALES MENSUALES DE LOS SITIOS
DE INTERS ............................................................................................................................66
8.2 ESCALA INTERANUAL: RELACIN DE LOS CAUDALES MEDIOS
MENSUALES CON EL FENMENO EL NIO-OSCILACIN DEL SUR (ENSO) ......................67
8.2.1 Anlisis serie de caudales del Proyecto Hidroelctrico Sogamoso ..........................67
8.2.2 Anlisis serie de caudales del Proyecto Hidroelctrico del ro Amoy .....................68
8.2.3 Anlisis serie de caudales de la Central Hidroelctrica Miel I ..................................68

Ingeniero Civil

vi

8.3 ANLISIS DE COMPLEMENTARIEDAD DE LAS SERIES HIDROLGICAS
CON LA SERIE DE VIENTOS LA GUAJIRA ..............................................................................69
8.4 ANLISIS DE GENERACIN DE ENERGA HISTRICA CENTRALES ISAGEN
VS. PARQUES ELICOS GENRICOS....................................................................................70
8.5 ANLISIS DE GENERACIN COMPLEMENTARIA PROYECTOS
HIDROELCTRICOS VS. PARQUES ELICOS .......................................................................72
8.5.1 Caractersticas del Modelo de Simulacin de Generacin de Energa.....................73
8.5.2 Funcionamiento del Modelo ...................................................................................74
8.5.3 Regla de Operacin del Embalse ...........................................................................74
8.5.4 Confiabilidad en la Meta .........................................................................................75
8.5.5 Representacin grfica de los resultados de la simulacin .....................................76
8.6 COMPARACIONES ENERGTICAS AISLADAS ...........................................................78
8.6.1 Proyecto Sogamoso Parque Elico Genrico 1.040 MW......................................78
8.6.2 Proyecto Amoy Parque Elico Genrico 140 MW ..............................................82
8.7 ANLISIS ENERGTICO INTEGRADO DE COMPLEMENTARIEDAD ..........................83
8.7.1 Anlisis Integrado P. H. Sogamoso Parques Genricos .......................................84
8.8 ANLISIS FINANCIERO DE LA GENERACIN INTEGRADA .......................................87
9 CONCLUSIONES ..................................................................................................................92
10 BIBLIOGRAFA Y FUENTES DE INFORMACIN..............................................................97


Ingeniero Civil

vii

LISTA DE FIGURAS

Figura 1 Distribucin de la Canasta Energtica Mundial, Energas Renovables al finalizar el
2008 (GW) ................................................................................................................. 13
Figura 2 Crecimiento anual de la capacidad instalada de Energa Elica. AWEA - Annual-
Wind-Report-2009. ..................................................................................................... 14
Figura 3 Distribuciones de densidad de energa elica. Para (a) de invierno boreal (DEF) y
(b) del verano boreal (JJA), por un perodo de ocho aos entre 2000 y 2007. La
escala de grises se utiliza para mostrar la topografa de la tierra. ............................... 15
Figura 4 Distribucin de la canasta de Generacin de Energa Elctrica en Colombia a
diciembre de 2008. UPME - Plan de Expansin de Referencia 2009-2023. ............... 18
Figura 5. Expansin de la Energa Elctrica Colombiana 2012 -2018. (ACOLGEN, 2009) ........... 18
Figura 6 Pronstico de produccin de gas. UPME - : Plan de Expansin de Referencia 2009-
2023. ......................................................................................................................... 20
Figura 7 Evolucin embalse agregado del SIN en Eventos ENSO - El Nio (XM, en lnea). ......... 21
Figura 8 Comportamiento del precio de Bolsa, Contratos y Mc. entre enero 2007 y octubre de
2009 (XM, 2009). ....................................................................................................... 21
Figura 9 Precios de oferta por combustible principal. (XM, 2009). ............................................... 22
Figura 10 Distribucin de potencia de energa elica en el Mar Caribe colombiano y en
particular de la zona de La Guajira Colombiana. Adaptada de Liu, Tang y Xie,
2008. ......................................................................................................................... 25
Figura 11 Distribucin de los costos de capital estimados de Parques elicos en Europa.
Adaptada de Blanco, (2008) ....................................................................................... 26
Figura 12 Complementariedad hidrulica - elica entre las series de caudales medios
mensuales del ro Nare y la velocidad promedia de vientos en la Alta Guajira
(Franco y Dyner, 2004). ............................................................................................. 30
Figura 13 Operacin aislada y conjunta de embalse de 0,2 del ro Guavio con parque
Jeprachi de igual capacidad instalada (Adaptada de Vergara et al, 2010). ................. 33
Figura 14 Operacin aislada y conjunta de embalse de 0,5 del ro Guavio con parque
Jeprachi de igual capacidad instalada (Adaptada de Vergara et al, 2010). ................. 33
Figura 15 Energa Firme para el ro Guavio como resultado de la Operacin Conjunta y
aislada (Adaptada de Vergara et al, 2010). ................................................................. 33
Figura 16 Ciclo diurno promedio de la velocidad del viento para la zona de estudio .................... 38
Figura 17 Estacin Manaure, comparacin ciclo diurno velocidad del viento ciclo diurno
temperatura ............................................................................................................... 39
Figura 18 Estacin Flor de la Sabana, comparacin ciclo diurno velocidad del viento ciclo
diurno temperatura ..................................................................................................... 39
Figura 19 Ciclo anual promedio de la velocidad del viento para la zona de estudio...................... 41
Figura 20 Variabilidad de las rosas de viento en la zona de estudio ............................................ 42
Figura 21 Estacin Puerto Bolvar, comparacin temperatura media mensual vs. velocidad
del viento mensual ..................................................................................................... 43

Ingeniero Civil

viii

Figura 22 Estacin Puerto Bolvar, comparacin precipitacin media mensual vs. velocidad
del viento mensual ..................................................................................................... 43
Figura 23 Estacin Puerto Bolvar, anlisis espectral de la velocidad promedia mensual del
viento ......................................................................................................................... 44
Figura 24 Estacin Puerto Bolvar, anlisis espectral de la velocidad promedia horaria del
viento ......................................................................................................................... 44
Figura 25 Estacin Flor de la Sabana - variacin de la velocidad del viento con respecto a la
altura ......................................................................................................................... 45
Figura 26 Estacin Manaure - variacin de la velocidad del viento con respecto a la altura ......... 45
Figura 27 Componentes bsicos de una turbina elica moderna, de eje horizontal con una
caja de cambios (Diseo por el Laboratorio Nacional de Energa Renovable -
NREL). Adaptada de SRREN IPCC, 2007. .............................................................. 46
Figura 28 Comparacin de la serie estandarizada de la velocidad del viento con el ndice
Ocenico del Nio (ONI) ............................................................................................ 47
Figura 29 Estacin Puerto Bolvar, comparacin de ciclos anuales de velocidad del viento,
aos El Nio, La Nia y normales............................................................................... 50
Figura 30 Estacin virtual Reanlisis, comparacin de ciclos anuales de velocidad del viento,
Figura 31 Estacin Ernesto Cortizzos, comparacin de ciclos anuales de velocidad del viento,
Figura 32 Estacin Puerto Bolvar, comparacin de la serie estandarizada de la temperatura
media mensual con el ndice Ocenico del Nio (ONI) ............................................... 52
Figura 33 Estacin Puerto Bolvar, coeficientes de correlacin para n meses de rezago de la
serie estandarizada de la temperatura media mensual con respecto al ndice
Ocenico del Nio (ONI) ............................................................................................ 52
Figura 34 Estacin Puerto Bolvar, comparacin de la serie de la temperatura con la
velocidad media mensual ........................................................................................... 52
Figura 35 Estacin Puerto Bolvar, ciclo anual del ciclo diurno de la velocidad del viento
(m/s). Promedio en aos normales. ............................................................................ 53
(m/s). Promedio en aos El Nio. ............................................................................... 54
(m/s). Promedio en aos La Nia. .............................................................................. 54
Figura 38 Localizacin de las estaciones sobre los ros seleccionados para estudio y
estacin Puerto Bolvar referencia para el anlisis de complementariedad............... 57
Figura 39 Anlisis de complementariedad regional entre ros del norte y noroeste colombiano
y la estacin Puerto Bolvar. ...................................................................................... 58
Figura 40 Anlisis de complementariedad regional entre ros del centro de Colombia y la
estacin Puerto Bolvar. ............................................................................................. 59
Figura 41 Anlisis de complementariedad regional entre ros del sur de Colombia y la
estacin Puerto Bolvar. ............................................................................................. 59
Figura 42 Curva de Potencia Terica Aerogenerador Nordex N60/1.300. (Nordex Energy
Catlogo N60/1300) ................................................................................................... 62
Figura 43 Comparacin de Tamao de algunos aerogeneradores tpicos. Crecimiento en el
tiempo. Adaptada de SRREN IPCC, 2007. .............................................................. 63
Figura 44 Ajuste a nivel mensual de velocidades del Viento vs. Potencia .................................... 64

Ingeniero Civil

ix

Figura 45 Comparacin de la produccin energtica de Parque Elico Simulado vs. Jeprachi
2004 2005. .............................................................................................................. 64
Figura 46 Comparacin de ciclo anual de produccin energtica promedio Parque simulado
20 MW vs. Jeprachi. .................................................................................................. 65
Figura 47 Localizacin de Sitios de Inters: Centrales, Proyectos y Estaciones de Velocidad
del Viento. .................................................................................................................. 66
Figura 48 Ciclo Anual de caudales medios mensuales para los sitios de inters .......................... 67
Figura 49 Complementariedad caudales y vientos, series mensuales multianuales
estandarizadas .......................................................................................................... 69
Figura 50 Comparacin entre la generacin mensual entre Miel I y un Parque Elico de 396
MW ............................................................................................................................ 70
Figura 51 Comparacin entre el ciclo anual de generacin mensual entre Miel I y un Parque
Elico de 396 MW ...................................................................................................... 71
Figura 52 Comparacin entre el ciclo anual de generacin mensual entre Jaguas y un Parque
Elico de 170 MW ...................................................................................................... 71
Figura 53 Comparacin entre el ciclo anual de generacin mensual entre San Carlos y un
Parque Elico de 1240 MW ........................................................................................ 72
Figura 54. Regla Normal de Operacin del Proyecto Sogamoso Operacin del embalse .......... 75
Figura 55. Proyecto Sogamoso Caudal medio mensual, caudal de diseo y Operacin del
embalse ..................................................................................................................... 76
Figura 56. Proyecto Sogamoso - Simulacin Energa Media y Firme Mensual ............................. 76
Figura 57. Proyecto Sogamoso Relacin Potencia Promedio Mensual vs. Capacidad
Instalada .................................................................................................................... 77
Figura 58. Proyecto Sogamoso Energa Media Mensual ............................................................. 77
Figura 59. Produccin Energtica Proyecto Sogamoso vs. Parque elico genrico que iguala
su Energa Media. ...................................................................................................... 78
Figura 60. Comparacin de Produccin Energtica de Parque 1.040 MW vs. Sogamoso 820
MW ............................................................................................................................ 79
Figura 61 Ciclo Anual de Produccin Energtica del Proyecto Sogamoso 820 MW: Normal, El
Nio y La Nia ........................................................................................................... 79
Figura 62 Ciclo Anual de Produccin Energtica del Parque 1040 MW: Normal, El Nio y La
Nia ........................................................................................................................... 80
Figura 63 Ciclo Anual de Produccin Energtica del Proyecto Sogamoso 820 MW y del
Parque 1.040 MW: Evento ENSO El Nio ............................................................... 80
Figura 64 Ciclo Anual de Produccin Energtica del Proyecto Sogamoso 820 MW y del
Parque 1040 MW: Evento ENSO La Nia ................................................................ 81
Figura 65 Diferencias de Produccin Energtica Parque 1.040 MW vs. Sogamoso 820 MW ....... 82
Figura 67 Produccin Energtica Proyecto Amoy vs. Parque elico genrico que iguala
energa Media ............................................................................................................ 82
Figura 67 Comparacin de Produccin Energtica de Parque 140 MW vs. Amoy 80 MW .......... 83
Figura 68 Simulacin de Generacin de Energa variacin meta o demanda vs. energa
generada ................................................................................................................... 85
Figura 69 Simulacin de Generacin de Energa variacin meta o demanda vs. energa
generada ................................................................................................................... 85
Figura 70 Variacin del precio de energa entre 2009 y 2023. Adaptada de UPME, 2009. ........... 89
Figura 71 Variacin VPN y TIR para los diferentes Parques Elicos Genricos ........................... 89

Ingeniero Civil

x

Figura 72 Proyecto Sogamoso: Variacin VPN y TIR, operacin integrada con beneficios. ......... 90
Figura 73 Proyecto Sogamoso: Variacin VPN y TIR, operacin integrada sin beneficios............ 90
Figura 74 Variacin VPN y TIR para los diferentes Parques Elicos Genricos con beneficios
de operacin integrada............................................................................................... 91


Ingeniero Civil

xi

LISTADO DE TABLAS

Tabla 1 Informacin de estaciones de registro de velocidad del viento utilizadas .......................... 36
Tabla 2 Matriz de correlaciones trimestrales entre el ndice ONI y la velocidad del viento en
la estacin Puerto Bolvar .............................................................................................. 48
Tabla 3 Matriz de correlaciones trimestrales entre el ndice SOI y la velocidad del viento en
la estacin Puerto Bolvar .............................................................................................. 48
Tabla 4 Matriz de correlaciones trimestrales entre el ndice ONI y la velocidad del viento
Nodo Reanlisis. ........................................................................................................... 49
Tabla 5 Matriz de correlaciones trimestrales entre el ndice ONI y la Temperatura en C .............. 52
Tabla 6 Informacin de estaciones de caudales utilizadas ............................................................ 56
Tabla 7 Matriz de correlaciones trimestrales caudales medios mensuales Sogamoso vs. ONI ...... 68
Tabla 8 Matriz de correlaciones trimestrales caudales medios mensuales Amoy vs. ONI ............ 68
Tabla 9 Matriz de Correlaciones Trimestrales Caudales medios mensuales Miel I vs. ONI ........... 68
Tabla 10 Matriz de correlaciones la velocidad del viento caudales: estacin Reanlisis,
estacin Puerto Bolvar vs. Caudales sitios de inters ................................................... 69
Tabla 11 Caractersticas del Proyecto Sogamoso para la simulacin integrada ............................ 84
Tabla 12 Resultados de la simulacin Integrada Sogamoso Parques Elicos Genricos ........... 86


Ingeniero Civil

xii

LISTA DE ABREBIATURAS

ENSO: El Nio Southern Oscillation - El Nio Oscilacin del Sur
ISAGEN S.A. E.S.P.: Empresa de servicios pblicos de carcter mixta, Generadora
y comercializadora de soluciones energticas
IDEAM: Instituto de Hidrologa, Meteorologa y Estudios Ambientales de Colombia
ONI: Oceanic Nio Index - ndice Ocenico Evento El Nio
SOI: Southern Oscillation Index - ndice Oscilacin del Sur
ZCIT: Zona de Convergencia Intertropical
NAO: North Atlantic Oscillation - Oscilacin del Atlntico Norte
IPCC: Intergovernmental Panel on Climate Change - Panel Intergubernamental de
Cambio Climtico
NCEP: National Centers for Environmental Prediction
NCAR: The National Center for Atmospheric Research
SIN: Sistema Interconectado Nacional
AWEA: American Wind Energy Association
VPN: Valor Presente Neto
TIR: Tasa Interna de Retorno
MDL: Mecanismos de Desarrollo Limpio
FAZNI: Fondo de Ayuda de Zonas no Interconectadas
UPME: Unidad de Planeacin Minero Energtica
PEN: Plan Energtico Nacional
EAS: Energas Alternativas
ANDI: Asociacin Nacional de Industriales
CREG: Comisin de Regulacin de Energa y Gas
ENFICC: Energa Firme para el Cargo por Confiabilidad
GEI: Gases de Efecto Invernadero
MAVDT: Ministerio de Ambiente Vivienda y Desarrollo Territorial
DEF: Trimestre diciembre enero febrero. Definicin similar para MAM, JJA, SON


Ingeniero Civil

13

1 INTRODUCCIN

Las energas renovables no convencionales
1
(ISAGEN S.A. E.S.P., 2005) usualmente
tienen significativamente menores impactos sociales y ambientales, pero cuando son
comparadas con la energa derivada de recursos convencionales; tales como los
combustibles fsiles, stas son sinnimo de riesgo financiero (Berry y Jaccard, 2000).
Estos menores impactos o beneficios ambientales pueden ser tanto locales, por los
menores aportes de emisiones contaminantes, como globales por la reduccin emisiones
de gases de efecto invernadero.

Segn Berry y Jaccard (2000), hay tres principales razones de la discrepancia entre los
beneficios sociales y econmicos de las energas renovables y su elevado costo
financiero, en relacin con las fuentes convencionales de generacin de tipo
contaminantes:

Algunas leyes proveen subsidios para las fuentes de generacin convencional.
Los costos totales de la contaminacin no estn incluidos en la evaluacin
financiera de las fuentes de generacin de electricidad convencional.
Las tecnologas renovables siempre han sido asociadas a los altos costos. Estos
costos relativos disminuirn en el tiempo con la comercializacin generalizada
debido a las economas de aprendizaje y de escala en la fabricacin de equipos.

Figura 1 Distribucin de la Canasta Energtica Mundial, Energas Renovables al finalizar
el 2009 (GW). Fuente: Renewables 2011 - Global Status Report, REN21. 2011.

1
Entindase en este documento como energas renovables no convencionales: Geotrmica, Solar Fotovoltaica, Elica,
Biomasa, Mareomotriz y la utilizacin de Hidrgeno.

Ingeniero Civil

14

En la Figura 1, se observa la distribucin actual las fuentes de generacin de energa
elctrica mundial. De las energas renovables, la hidroelctrica es la de mayor porcentaje
con solo el 3%, y la elica solo llega al 0,3% de la generacin mundial. Es clara entonces
la necesidad urgente de encontrar los mecanismos que tiendan a disminuir el consumo de
las fuentes de generacin fsil, de lo contrario ser prcticamente imposible detener el
cambio global.

Una de las principales fuentes de generacin de energa que est aumentando su
participacin en el mercado es la energa elica. Segn la AWEA, al finalizar el 2008, se
instalaron alrededor del mundo ms de 27.000 MW, lo que aument el total mundial en
capacidad instalada de energa elica a 120.000 MW (Ver Figura 2).

Figura 2 Crecimiento anual de la capacidad instalada de Energa Elica. Fuente: AWEA - Annual-
Wind-Report-2009.

Segn el trabajo de Joselin e Iniyan (2005), an queda un gran potencial aprovechable.
Se ha estimado que aproximadamente 10 millones de MW de energa estn
continuamente disponible en el viento de la tierra. Por su parte Pryor y Barthelmine
(2009) argumentan que la expansin de la capacidad instalada de energa elica est
preparada para desempear un papel clave en la mitigacin del cambio climtico.
Sin embargo, los parques elicos tienen asociados altos costos de instalacin si se
comparan con los de las de otras fuentes convencionales. Colombia no es la excepcin;
hoy en da comparar los costos asociados a la construccin y operacin de un parque
elico es mucho menos rentable que el de las centrales hidroelctricas de tamao
comparable. Sin embargo, en los ltimos aos dicho costo ha disminuido
significativamente y se espera que siga disminuyendo dado el continuo crecimiento del
mercado de la energa elica, tal y como se observa en la Figura 2.
A pesar que la regulacin y normatividad del Mercado de Energa en Colombia est
diseada para la construccin y operacin de centrales Hidroelctricas y Termoelctricas,

Ingeniero Civil

15

no se ha cuantificado claramente cul sera el beneficio ambiental de combinar la energa
elica e hidrulica, de tal manera que se desplazara la generacin trmica. Como se
presentar ms adelante, se ha demostrado la existencia de complementariedad de las
dos fuentes de generacin; la elica en La Guajira colombiana, y la hidrulica en la zona
central del pas.
Por otro lado, y analizando el rgimen de vientos a nivel global, en la Figura 3 se aprecia
en los mapas, construidos a partir de 8 aos de mediciones de redar (tcnica basada en
la utilizacin de dispersimetro), que la potencia promedio de los vientos del mundo, tiene
los mximos en Diciembre, Enero y Febrero correspondiente al invierno boreal (norte del
planeta), y en los meses de Junio, Julio y Agosto los mnimos para la misma zona (Liu,
et. al, 2008).

Figura 3 Distribuciones de densidad de energa elica. Para (a) de invierno boreal (DEF) y (b)
verano boreal (JJA), por un perodo de ocho aos entre 2000 y 2007. La escala de grises se utiliza
para mostrar la topografa de la tierra. Tomado de Liu, et. al, 2008.
Una de las desventajas de la energa elica es su variabilidad en el tiempo, lo que
aparentemente disminuye su confiabilidad para el abastecimiento de las necesidades de
energa elctrica.

La generacin de energa poco a poco est evolucionando hacia las fuentes renovables, y
las empresas que no quieran quedarse rezagadas debern enlistarse en el aprendizaje
tecnolgico de las fuentes no convencionales para generar electricidad.


Ingeniero Civil

16

En Colombia, dentro de la caracterizacin y diagnstico general de la zona de vientos de
La Guajira colombiana, las Empresas Pblicas de Medelln (EPM) han jugado un papel
destacado, sobre todo por la experiencia de cerca de 10 aos en estudios y 6 en
operacin comercial de su parque elico piloto Jeprachi en la Guajira. No obstante, al
tratarse de un negocio estratgico, tal informacin no es pblica y los anlisis con relacin
a la velocidad del viento a la altura del aerogenerador y las bitcoras de operacin y
mantenimiento del Parque son estrictamente confidenciales.

Este documento estudia diferentes aspectos que se deben considerar previos a la
caracterizacin energtica del recurso elico teniendo en cuenta su variacin intradiaria e
intranual. Previo a lo anterior en el captulo 3, se presentan elementos a manera de
justificacin con los cuales se considera que diversificar la matriz de generacin de
energa en Colombia a partir de otras fuentes no convencionales de energa se requiere
con carcter prioritario, describiendo el estado actual de las fuentes de generacin y el
comportamiento del mercado desde su creacin.
Seguido a lo anterior, en el captulo 4 se presentan elementos de carcter terico que
respaldan el hecho de incluir en la matriz el potencial de la energa elica, algunos
aspectos tcnicos de la tecnologa y su complementariedad con el recurso hdrico.
En el captulo 5 analiza el recurso elico en La Guajira colombiana, presentando los ciclos
o variacin diaria, mensual y anual del recurso, su correlacin con datos del reanlisis
NCEP/NCAR y el impacto de las dos fases del sistema El Nio-Oscilacin del Sur (ENSO
por sus siglas en ingls, El Nio Southern Oscillation) sobre la magnitud de las
velocidades de los vientos en la atmsfera inferior. Se incluye tambin un anlisis de la
direccin del viento en varias estaciones de la zona norte de Colombia (la mayora sobre
el departamento de La Guajira) y adicionalmente se presenta con datos de campo cmo
vara la magnitud de la velocidad del viento con la altura.
En el captulo 6 se realiza un anlisis del comportamiento promedio del ciclo anual de los
caudales y velocidades del viento en tres regiones del pas: norte, centro y sur. En el
captulo 7 se realiza un anlisis similar con la informacin de los Proyectos y Centrales
Hidroelctricas de ISAGEN, con el objeto de identificar cul o cules se complementaran
mejor con parques elicos de diferentes capacidades instaladas ubicados en La Guajira.
En el captulo 8 se realiza un anlisis de generacin teniendo en cuenta, para parques
elicos genricos un aerogenerador en particular, igual a los del Parque Jeprachi y para
los proyectos hidroelctricos el desarrollo e implementacin de un modelo de operacin
de embalses mediante la Regla Normal de Operacin, y teniendo en cuenta un escenario
de comparacin de generaciones individuales y otro integrado. Los anlisis consideran el
componente econmico, con el fin de cuantificar los beneficios para un inversionista como
ISAGEN S.A. de la complementariedad del recurso hidro elico.


Ingeniero Civil

17

2 OBJETIVOS
2.1 OBJETIVO GENERAL

Caracterizar el recurso elico y analizar la complementariedad de la generacin hidrulica
del centro de Colombia con la generacin elica en La Guajira colombiana, y cuantificar
los beneficios econmicos producto del aprovechamiento de esta fuente de generacin de
energa no convencional.

2.2 OBJETIVOS ESPECFICOS

- Caracterizar la variabilidad del recurso elico en la zona de La Guajira colombiana en
distintas escalas de tiempo (desde la diurna hasta la interanual) a partir de informacin
de torres de medicin meteorolgicas, as como de modelos climticos.

- Analizar la complementariedad intra-anual entre los recursos de generacin elico e
hidrulico, en las diferentes regiones del pas.

- Analizar la posible complementariedad natural entre el recurso elico en La Guajira
colombiana y el recurso hdrico en la zona de las centrales y proyectos hidroelctricos
de ISAGEN S.A. E.S.P.

- Cuantificar aislada e integradamente con un proyecto hidroelctrico con embalse, la
complementariedad energtica de parques elicos de diferentes tamaos.

- A partir de estimaciones financieras, cuantificar la capacidad instalada (MW) necesaria
para que un parque elico sea rentable teniendo en cuenta la cuantificacin del
beneficio asociado a la complementariedad, y de manera general las restricciones
legales, de mercado y regulatorias de Colombia.


Ingeniero Civil

18

3 JUSTIFICACIN

Colombia posee un alto potencial en recurso elico disponible principalmente en La
Guajira. A la fecha solo es aprovechado por el Parque Piloto Jeprachi de propiedad de
las Empresas Pblicas de Medelln. Este parque de 19.5 MW.

La matriz de potencia o capacidad instalada energtica nacional es hdrica-trmica con
valores porcentuales aproximados del orden de 65% - 35% respectivamente, sin embargo
a la generacin de energa de estas capacidades instaladas, vara segn la poca del
ao, la influencia del Evento ENSO, entre otros factores. En la Figura 4 se observa la
distribucin de la capacidad instalada en MW para generacin energa elctrica.

Figura 4 Distribucin de la Capacidad Instalada de Generacin de Energa Elctrica en Colombia
(MW) a diciembre de 2008. Fuente: UPME - Plan de Expansin de Referencia 2009-2023.

Para los prximos aos se espera un reajuste de estos porcentajes. En la Figura 5 a y b,
se observa el alto porcentaje en capacidad instalada hidrulica, alrededor de un 66%, y
para el 2018 se espera que esta capacidad aumente a un 72%. Como es de comn
conocimiento, este tipo de generacin produce muy pocas emisiones de Gases de Efecto
Invernadero - GEI y ms an, si se compara con otros tipos de generacin a partir de
combustibles fsiles (trmica a carbn, trmicas a gas, generacin Diesel, etc.).

2008
a

2018
b
Figura 5. Expansin de la Energa Elctrica Colombiana 2012 -2018. (ACOLGEN, 2009)


Ingeniero Civil

19

El aumento de la capacidad instalada en su mayora hidrulica, obedece a las
obligaciones de energa firme obtenidas por los diferentes agentes mediante un proceso
de subasta de energa firme regulado por la Comisin de Regulacin de Energa y Gas
CREG. La Obligacin de Energa Firme (OEF) es un producto diseado para garantizar la
confiabilidad en el suministro de energa en el largo plazo a precios eficientes. En esta
subasta la mayor participacin correspondi a centrales hidroelctricas.

Sin embargo, los grandes embalses para proyectos hidroelctricos de alta capacidad
instalada ocasionan fuertes impactos ambientales y sociales negativos. A pesar de no
ser el objeto de este trabajo, el anlisis de la generacin de GEI producto del
almacenamiento de agua de los grandes embalses, segn informacin disponible en el
sitio Web de HydroQuebec, los grandes embalses de las hidroelctricas inducen la
descomposicin de una pequea fraccin de la biomasa inundada (bosques, materia
orgnica y tipos de suelo) y un aumento de la fauna acutica y la vegetacin en el
depsito. El resultado es una mayor concentracin gases de efecto invernadero (GEI),
despus del llenado del embalse, principalmente CO
2
(dixido de carbono) y una pequea
cantidad de CH
4
(metano). Sin embargo, estas emisiones son temporales y con un pico de
dos a cuatro aos despus de que el depsito est lleno. Durante la siguiente dcada, las
emisiones de CO
2
disminuyen gradualmente y vuelven a los niveles emitidos por los lagos
y los ros vecinos. Ver Tremblay, A., Schetagne, R. 2006; y Tremblay et. al. 2004.

Las restricciones ambientales cada vez son mayores y necesarias para la conservacin
de la biodiversidad nacional y la sobre-dependencia del agua como recurso primario para
la generacin de energa elctrica del pas, debern empujar a la diversificacin de la
canasta energtica. Segn el Informe de Sntesis del Grupo de Trabajo I del IPCC, los
resultados de algunos escenarios de la relacin directa del cambio climtico con el evento
El Nio Oscilacin del Sur - ENSO indican que si es probable que el calentamiento
mundial produzca una mayor variabilidad y extremos en las lluvias y un aumento del
riesgo de sequas y crecidas asociadas a las dos fases del ENSO (El Nio y La Nia), lo
que impacta directamente a la generacin hidroelctrica. Ver IPCC, Cambio climtico
2007: Resmenes de los Grupos de trabajo.

En la actualidad el respaldo y confiabilidad que tiene el Sistema Interconectado Nacional
- SIN, se soporta en la capacidad instalada trmica que a 2008 es del orden del 33% (ver
Figura 4), del cual el 25% corresponde a generacin con gas natural, el cual hace parte
de las fuentes de generacin de energa elctrica con combustible fsil y segn la Unidad
de Planeacin Minero Energtica - UPME, en los escenarios futuros de disponibilidad en
Colombia, el pronstico de la produccin de dicho combustible a mediano plazo muestra
un crecimiento hasta el ao 2011, momento a partir del cual, como consecuencia de la
declinacin de las reservas proyectadas declaradas, se iniciara una disminucin en la
produccin de los mismos.

Pese a que el sector de gas natural ha tenido importantes avances en los temas de
desarrollo de infraestructura y cobertura del servicio en Colombia, en el mediano plazo,
este desarrollo podra no ser suficiente para atender la creciente demanda de este
energtico. Ver Figura 6.

Ingeniero Civil

20

Figura 6 Pronstico de produccin de gas. Fuente: UPME: Plan de Expansin de Referencia 2009-
2023.

Conjugando los anteriores elementos y teniendo en cuenta la futura demanda de energa
elctrica del pas, as como la propuesta de interconexin con Centroamrica, la
disponibilidad del gas para generacin y los otros mercados demandantes del recurso
cada vez estar sujeta a mayores restricciones. Lo anterior sugiere la posibilidad de que a
que la confiabilidad podra buscarse a travs de la construccin de plantas termoelctricas
a carbn, dada la abundancia de este recurso en el pas, sin embargo, este es otro
combustible fsil.

Dado lo anterior, la dependencia del recurso hdrico; los posibles eventos El Nio fuertes y
duraderos, el calentamiento global, la deforestacin misma en las cuencas de los ros
colombianos, la incertidumbre con la produccin y precios del gas, reforzara la
vulnerabilidad anteriormente mencionada.

El evento de El Nio 19911992, impact la oferta de energa elctrica producindose un
racionamiento entre marzo de 1992 y abril de 1993. En la Figura 7 se muestra la variacin
del embalse agregado del sistema, el impacto de El Nio mencionado y la baja respuesta
del SIN ante el evento por indisponibilidad de las plantas trmicas.

Las consecuencias de la ocurrencia de El Nio 1991-1992 reflejaron la dependencia del
sistema, predominantemente hdrico con respecto a la variabilidad natural del recurso
hdrico, y por ende su vulnerabilidad y dependencia de la energa trmica de respaldo. La
experiencia ayud a plantear a largo plazo la implementacin de un plan de emergencia
con la puesta en operacin de 2.138 MW entre 1992 y 1994 con la puesta en marcha en
general de plantas trmicas, adicionalmente, oblig la reforma del sector elctrico en
bsqueda de la eficiencia en aspectos como precios, confiabilidad, calidad del servicio y
cobertura (XM, 2009). Se considera que esta implementacin se realiz sin la
planificacin adecuada y sin estudios detallados de optimizacin de la canasta energtica
y ms como respuesta de miedo ante racionamientos.

Ingeniero Civil

21

Figura 7 Evolucin embalse agregado del SIN en Eventos ENSO - El Nio. Tomada de XM,
Informe Especial Resolucin CREG 135/97, 2009

En la Figura 7 se observa el respaldo trmico de generacin durante el siguiente evento
de El Nio 1997-1998. En este caso no se present racionamiento de energa.

Es claro que al incluir un mayor porcentaje de generacin trmica, el precio de
generacin, transmisin y venta de energa aumenta. Un ejemplo es el presentado
durante los aos 2008 y 2009. La Figura 8 muestra el comportamiento del precio de bolsa
(valor de compra y venta intradiurnas del sistema interconectado Colombiano),
presentando un crecimiento anual de 56%, al pasar de un promedio en 2008 de 89,14
$/kWh a 139,57 $/kWh en diciembre de 2009. Sin embargo, el incremento en los precios
de contratos entre generadores y consumidores no regulados fue del 18%, con 88,78
$/kWh en 2008 y 104,74 $/kWh en 2009, pesos constantes de diciembre de 2009 (XM,
2009).

Figura 8. Comportamiento del precio de Bolsa y Contratos entre enero 2007 y octubre de 2009
(XM, 2009).

El comportamiento anterior puede observarse tambin en los precios de oferta por
combustible presentado en la Figura 9, donde el precio de oferta del agua est totalmente
ligado al precio de bolsa, sobre todo en los meses de agosto a diciembre. Segn la NOAA

Ingeniero Civil

22

el Evento ENSO 2009 2010 inici desde el mes de junio con un ndice ONI de 0,6
(www.cpc.ncep.noaa.gov).

Figura 9. Precios de oferta por combustible principal. (XM, 2009).

Lo expuesto anteriormente sugiere que si en los meses secos y extremadamente secos
(durante la ocurrencia de El Nio) el precio de la energa en bolsa aumenta, un recurso
capaz de generar ante dicha eventualidad como sera el elico obtendra mejores
beneficios financieros en esos perodos.

Cabe anotar que segn la regulacin actual del mercado elctrico colombiano las plantas
elicas operaran como centrales hidroelctricas a Filo de Agua, es decir generaran todo
el recurso disponible y entraran en la base de la generacin para ayudar a cubrir la
demanda diaria. Lo anterior sugiere que los Parques Elicos generaran constantemente,
siempre y cuando estn disponibles.

Como se mencion anteriormente, en energa generada, la proporcin es hidrulica 77%
aproximadamente y trmica 18%; al depender de los aportes hidrolgicos, las pocas de
sequa que se presentan durante evento ENSO- El Nio hacen indispensable contar con
plantas de generacin con energa firme, que remplacen la energa generada por
hidroelctricas, para atender la demanda. De no contar con estos recursos, los usuarios
tendran que ser racionados, con los correspondientes costos sobre la economa nacional
y el bienestar de la poblacin.

Uno de los principios subyacentes en un sistema de precios como el diseado para el
MEM en Colombia, es que este debe proporcionar la seal econmica de largo plazo para
la expansin de la capacidad instalada requerida por el pas. Asimismo, la evolucin y el
comportamiento de los precios deben reflejar el nivel de confiabilidad en el suministro que
est dispuesta a pagar la demanda nacional.

Sin embargo, la volatilidad de los precios en la bolsa, y presentado en la Figura No. 8, que
se explica en gran parte por el elevado componente hidrulico, la estacionalidad climtica
(siete meses de invierno y cinco meses de verano) y la aparicin peridica de El Nio,
puede constituir un riesgo considerable para aquellos generadores que deben disponer de
fuentes de financiacin de sus proyectos de generacin, si no se cuenta con mecanismos
que cubran estas eventualidades. Por estas consideraciones, se encontr indispensable

Ingeniero Civil

23

implementar un esquema de remuneracin que permita hacer viable la inversin en los
recursos de generacin necesarios atender la demanda de manera eficiente en
condiciones crticas de abastecimiento hdrico, a travs de la estabilizacin de los
ingresos del generador. Este esquema est incorporado en la legislacin colombiana
desde la Ley 143 de 1994, artculo 23.

Tras diez aos de aplicacin ininterrumpida del Cargo por Capacidad, la CREG dise un
nuevo esquema basado en un mecanismo de mercado denominado Cargo por
Confiabilidad, que opera desde el primero de diciembre de 2006. Este mecanismo
conserva lo esencial del esquema de liquidacin, facturacin y recaudo que garantiz, con
xito, durante los diez aos continuos el pago a los generadores del Cargo por
Capacidad.

Uno de los componentes esenciales del nuevo esquema es la existencia de
las Obligaciones de Energa Firme (OEF), que corresponden a un compromiso de los
generadores respaldado por activos de generacin capaces de producir energa firme
durante condiciones crticas de abastecimiento. Este nuevo esquema permite asegurar la
confiabilidad en el suministro de energa en el largo plazo a precios eficientes.

Para estos propsitos, se subastan entre los generadores las OEF que se requieren para
cubrir la demanda del Sistema. El generador al que se le asigna una OEF recibe una
remuneracin conocida y estable durante un plazo determinado, y se compromete a
entregar determinada cantidad de energa cuando el precio de bolsa supera un umbral
previamente establecido por la CREG y denominado Precio de Escasez. Dicha
remuneracin es liquidada y recaudada por el ASIC y pagada por los usuarios del SIN, a
travs de las tarifas que cobran los comercializadores. Las OEF del nuevo Cargo por
Confiabilidad establecen un vnculo jurdico entre la demanda del MEM y los generadores,
que permite, tanto a generadores como a usuarios del sistema, obtener los beneficios
derivados de un mecanismo estable en el largo plazo y que da seales e incentivos para
la inversin en nuevos recursos de generacin, garantizando de esta forma el suministro
de energa elctrica necesario para el crecimiento del pas.


Ingeniero Civil

24

4 ESTADO DEL ARTE

4.1 Del diagnstico y caracterizacin del potencial elico de la zona de estudio

La expansin de la capacidad instalada de energa elica puede jugar un papel
determinante en la mitigacin del cambio climtico. Sin embargo, segn Pryor y
Barthelmie (2009), la energa elica tambin es susceptible al cambio climtico global. En
el norte de Europa han identificado aumento en las velocidades del viento para los ltimos
aos del siglo XXI; sin embargo, los modelos usados an llevan consigo una gran
incertidumbre.

De ah la importancia de cuantificar, diagnosticar y disminuir la incertidumbre respecto a la
disponibilidad del recurso en la alta Guajira colombiana. Hoy se tiene conocimiento de las
caractersticas geogrficas e hidrolgicas de la zona. La Guajira se identifica con orografa
plana, seca, semidesrtica y con su territorio dividido en dos zonas con caractersticas
diferentes: la alta Guajira, ms al norte, semidesrtica y con poca vegetacin, y la baja
Guajira, menos rida y con alguna vegetacin, especialmente pastos para la ganadera y
con produccin agrcola en las estribaciones de la sierra Nevada de Santa Marta. Dentro
de la alta Guajira se destaca la presencia de la serrana La Macuira (Parque Nacional
Natural), ecosistema nico debido a la presencia de montaas y bosques que contrastan
con la geografa desrtica de la zona.

En cuanto al recurso elico, podra decirse que el avance ms importante respecto a su
caracterizacin ha sido la publicacin del Atlas de Viento y Energa Elica de Colombia,
UPME IDEAM, 2006. Dicho estudio se constituye en la herramienta bsica para la
estimacin preliminar del potencial y la identificacin en el territorio nacional de zonas que
mantienen velocidades de viento altas y persistentes a lo largo del ao. Se destacan
dentro de ellas, la pennsula de La Guajira, el centro del litoral Caribe en Bolvar, la Isla de
San Andrs y algunos sectores de Boyac (Isaac et al., 2008).

La Guajira posee un potencial energtico estimado entre 5 y 20 GW, los vientos son
predominantemente alisios
2
con una velocidad promedio de 9,8 m/s a 50 metros de altura
y baja turbulencia (Isaac, 2009).

La variacin intra - anual de la velocidad del viento (el ciclo anual) est relacionada entre
otras cosas, con los perodos de precipitacin y sequa directamente influenciadas por la
migracin meridional de la Zona de Convergencia Intertropical (ZCIT) durante el ciclo
anual, la cual acta sobre todo el territorio nacional (Poveda, 2004). Adicional a esto,
dada la posicin de la zona de estudio (Norte de Colombia) la ZCIT solo pasar una vez al
ao, indicando tambin el porqu del ciclo unimodal que se presenta en la velocidad del
viento y la precipitacin en la zona.

2
Circulan entre los trpicos, desde los 30-35 de latitud hacia el ecuador. Se dirigen desde las altas presiones
subtropicales, hacia las bajas presiones ecuatoriales. El movimiento de rotacin de la Tierra desva a los alisios hacia el
oeste, y por ello soplan del nordeste al sudoeste en el hemisferio norte y del sudeste hacia el noroeste en el hemisferio sur.

Ingeniero Civil

25

Para el ciclo diurno, la variabilidad del viento est asociada a la distribucin de la
temperatura, resultante del ciclo de insolacin superficial en las 24 horas del da
fuertemente marcada en la zona tropical (Poveda, 2004).

Retomando los resultados de Liu, Tang y Xie (2008), con relacin a la velocidad del viento
en la zona martima de La Guajira colombiana, en la Figura 10 se presentan las
densidades de energa elica a partir de datos de radar de ocho (8) aos de mediciones.
En esta figura pueden observarse las particularidades en la zona, marcadas como una
concentracin de mximas densidades de potencia elica en la Figura 3.

Figura 10 Distribucin de potencia de energa elica en el Mar Caribe colombiano y en particular de
la zona de La Guajira Colombiana. Adaptada de Liu, Tang y Xie, 2008.

4.2 Costos de referencia en Parques Elicos Internacionales. Regulacin y
Normatividad Colombiana

Segn Blanco (2008), hoy en da en la Unin Europea, los costos de generacin de un
parque elico en tierra estn entre 4,5 y 8,7 cent/kWh. Los Costos de produccin han
aumentado en ms del 20% en los ltimos tres (3) aos, principalmente debido a un
aumento del precio de algunas materias primas estratgicas en un momento en que la
demanda mundial se ha disparado. Sin embargo, la posicin competitiva de las
inversiones en energa elica versus otras tecnologas no ha sido alterada. En el largo
plazo, cabra esperar que los costos de produccin bajen; y que esto sea suficiente para
compensar el alto precio de los insumos, que en su mayora estn representados por los
aerogeneradores como puede apreciarse en la Figura 11. Depender en gran medida la
aplicacin de polticas correctas, como Investigacin y Desarrollo (I+D) para nuevos
materiales, Operacin y Mantenimiento (O & M) con dispositivos de control remoto, la
energa elica marina, las turbinas y subestructuras, la introduccin de tecnologas
avanzadas en el emplazamiento, las tcnicas de pronstico, el acceso a una financiacin
adecuada, y la estabilidad jurdica, como medida de largo plazo (Blanco, 2008).

Ingeniero Civil

26

Figura 11 Distribucin de los costos de capital estimados de Parques elicos en Europa. Adaptada
de Blanco (2008).

Esta situacin podra ser aplicable en Colombia considerando un escenario de instalacin
de grandes parques elicos prximos aos en Colombia, dado a que los aerogeneradores
seran importados desde Europa, en primera instancia, Estados Unidos, o tal vez el
mercado emergente de algunas plantas de fabricacin en Brasil o Argentina, como el caso
de IMPSA. Sin embargo se requerirn incentivos y polticas de fomento a la generacin
elica.

De ah la importancia de esta investigacin, la cual plantea para la energa elica,
escenarios en los cuales se consideren los beneficios ambientales, el futuro de los costos
de maquinarias y equipos para los parques elicos, as como tambin la regulacin,
normatividad y precios del mercado elctrico colombiano.

A la fecha, el sistema de interconexin de la zona norte del pas cuenta con una
estructura radial, vulnerable y en el que no se prevn proyectos en el corto plazo de lneas
de transmisin. Adems cabe destacar que precisamente la zona de mayor potencial
elico tiene una alta capacidad instalada de centrales trmicas y para la generacin elica
existiran desventajas asociadas con las rampas de entrega de potencia energtica,
producto de la variabilidad intradiurna de la velocidad del viento. Adicional a esto, no
existen en el mercado elctrico colombiano criterios especficos para reglamentar la
operacin de parques elicos (Isaac, 2009).

Como punto a favor, en el ao 2003, se cre el Comit 185 de Energa elica por parte
del ICONTEC. A la fecha ha desarrollado dos anteproyectos de NTC: Requisitos de
seguridad de los sistemas de generacin elicos, Ensayo Curva de Potencia de
Aerogeneradores, y algunas adaptaciones de las normas IEC 61400 (UPME ICONTEC,
2008). Estos proyectos apuntan a la implementacin de una regulacin competente dentro
del MEM de Colombia.

Isaac et al. (2008) resaltan la iniciativa importante del Plan Energtico Nacional 2006-
2025, en el cual se plantean varias estrategias relacionadas con la implementacin de las
energas renovables, principalmente enfocadas a la alimentacin de zonas rurales y no
interconectadas. Las sugerencias apuntan a la definicin de subsidios directos y la
reestructuracin de otros existentes, as como el fortalecimiento de los esquemas
institucionales para promocin y desarrollo.


Ingeniero Civil

27

Por su parte, la Unidad de Planeacin Minero Energtica UPME con el apoyo del Banco
Mundial y la firma consultora Estudios Energticos, bajo la motivacin de analizar algunas
polticas que permitan la promocin de las energas renovables, evaluaron un posible
escenario de expansin en generacin contemplando la participacin de energa elica,
realizaron diversos anlisis al comportamiento de la expansin propuesta considerando
algunos escenarios y analizaron el comportamiento de los resultados de la expansin
propuesta sobre un modelo financiero (UPME ANDI, 2008).

Como resultado principal, el estudio revel que en el sistema se hace necesario
desarrollar un cargo de sostenibilidad, debido a que para conseguir la TIR que espera un
inversionista con los mecanismos y costos actuales los ingresos no seran suficientes
para la implementacin de ste tipo de alternativas (UPME ANDI, 2008).

En el estudio de Vergara et. al. 2010, se exploran varias alternativas enmarcadas dentro
de la normativa colombiana para evaluar la construccin de parques elicos dentro de una
viabilidad financiera. Dentro de estos se plantean la participacin activa en los MDL, la
reduccin del impuesto sobre la renta, lo que mejorara la rentabilidad de los proyectos;
exenciones de los cargos del SIN, como contribuciones ambientales o para el FAZNI
(Fondo de Ayuda para Zonas no Interconectadas) y sobre todo una reestructuracin
regulatoria para todas las fuentes que participan en el mercado. Este ltimo aspecto
impacta tanto a la generacin elica, que si se evala un parque elico de 300 MW con un
costo de instalacin de $1.800/kW e ingresos producto de una confiabilidad entre un 20%
y un 36% atraeran la inversin privada, volviendo viable; siempre y cuando no se incluya
el pago por confiabilidad y los cargos del SIN. El proyecto alcanzara una TIR del 14%.

Con relacin a la normatividad y regulacin en Colombia se destacan las siguientes:

- Ley 697 de 2001. Formula el establecimiento de lineamientos de las polticas,
estrategias e instrumentos para el fomento y promocin de las fuentes no
convencionales de energa, con prelacin en las ZNI.

- El Plan Nacional de Desarrollo 2006 2010, establece como una de sus metas
sustituir capacidad proveniente de combustibles fsiles con Energas Alternativas
(EAS) de 0,35 MW a 6,9 MW en las Zonas No Interconectadas ZNI.

- El PEN plantea el uso de las EAS como una alternativa vlida de abastecimiento
energtico que disminuye efectos del cambio climtico y evita la destruccin de
ecosistemas.

A pesar de los adelantos mundiales las EAS, en Colombia no han aumentado como se
esperaba debido a los altos costos de inversin y el bajo conocimiento de las tecnologas,
as como la competencia de otras fuentes de generacin de ms bajo costo. Las ZNI
pueden ser propicias para el desarrollo y aplicacin de las EAS debido a los altos costos
de operacin de las plantas que all suministran la energa elctrica, las cuales en su
mayora generan con combustibles fsiles.

El primer proyecto elico en Colombia, Jeprachi, se hizo viable gracias a factores como el
nivel de velocidad de los vientos en La Guajira, las econmicas alcanzadas por disminuir

Ingeniero Civil

28

costos de transporte por su cercana a Puerto Bolvar, las facilidades de transporte
terrestre, y la zona aduanera especial y el impacto mnimo sobre el medio ambiente
(Isaac, 2008).

Otros beneficios que promueven el desarrollo de proyectos de energas renovables en
Colombia son los Incentivos tributarios, explcitos en la Ley 778 de 2002 en los artculos
18 y 95:

- Artculo 18. Exencin de renta por quince aos:
- Generacin de electricidad con energa elica, biomasa y residuos agrcolas.
- Reduzcan emisiones y permitan obtener certificados dentro del esquema del
Mecanismo de Desarrollo Limpio.
- 50% de los ingresos provenientes de la venta de certificados de reduccin de
emisiones se destina a inversin social.

- Artculo 95. Exclusin de IVA:
- Importacin de equipos y maquinaria para proyectos que generen certificados de
reducciones de emisiones GEI.

Por ltimo, se resalta el actual proyecto de resolucin CREG 092 de 2011, por la cual se
define la metodologa para determinar la Energa Firme para el Cargo por Confiabilidad
ENFICC de plantas elicas. En este proyecto se presentan dos escenarios definidos por
Parques existentes y proyectos como tal. Para el primero el clculo considera la
informacin promedia mensual del viento, el ndice de indisponibilidad histrica forzada, el
factor de conversin y la capacidad efectiva neta del Parque. Para el caso de los
proyectos, si la informacin disponible de vientos es inferior a 10 aos en promedios
mensuales, la energa firme se calcula con un factor constante multiplicado con la
capacidad efectiva neta de la planta (CEN). Estos factores se estimaron a partir de la
informacin histrica de generacin del Parque Elico Jeprachi. En el captulo 8 se
retoman algunos aspectos mencionados anteriormente.
4.3 Del Protocolo de Kioto y los Mecanismos de Desarrollo Limpio - MDL

Uno de los resultados ms importantes de la cumbre sobre Medio Ambiente y Desarrollo
de 1992 fue la aprobacin de tres tratados internacionales jurdicamente vinculantes,
entre los cuales se encuentra la Convencin Marco de las Naciones Unidas sobre el
Cambio Climtico - CMNUCC. Hasta el momento, cerca de 186 Estados hacen parte de la
CMNUCC. Lo anterior hace que la CMNUCC sea uno de los acuerdos sobre desarrollo
sostenible con mayor respaldo poltico a nivel global (MAVDT, 2002).

El objetivo de la CMNUCC es estabilizar las concentraciones de gases efecto invernadero
en la atmsfera, a un nivel que impida efectos peligrosos en el sistema climtico.

El Protocolo de Kioto dentro de la Convencin Marco de Naciones Unidas sobre Cambio
Climtico, busca desarrollar el objetivo de la Convencin y fortalecer la respuesta de la
comunidad internacional al cambio climtico.


Ingeniero Civil

29

Establece para los pases incluidos en el Anexo I de la Convencin, obligaciones
cuantificadas de reduccin de emisiones de Gases de Efecto Invernadero - GEI. En su
conjunto promediado estos pases debern reducir en un 5,2% sus emisiones colectivas
de seis GEI, para el periodo comprendido entre los aos 2008 y 2012. Este compromiso
de reduccin de emisiones es diferenciado: la Unin Europea debe reducir sus emisiones
en un 8%, Estados Unidos 7%, Japn y Canad 6%, mientras Australia podr aumentar
un 8% (MAVDT, 2002).

Los pases en desarrollo como Colombia, que hacen parte de la CMNUCC pero no estn
incluidos en su Anexo I, es decir no cuentan an con una meta concreta y obligatoria de
reduccin de emisiones de Gases de Efecto Invernadero y se les permite incluso un
aumento de emisiones, en razn a que se les reconoce su derecho a avanzar en la
realizacin de su modelo de desarrollo econmico y social (MAVDT, 2002).

El Protocolo establece una serie de mecanismos que dan flexibilidad a las partes para el
cumplimiento de sus compromisos: Comercio de Emisiones, que permite que los pases
industrializados comercien entre s los derechos de emisin; la implementacin conjunta
entre las Partes del Anexo I, en virtud del cual los pases podrn intercambiar unidades de
reduccin de emisiones que resulten del desarrollo de proyectos encaminados a mitigar
las emisiones de GEI; y el Mecanismo de Desarrollo Limpio para implementar entre
pases del Anexo I y pases en desarrollo (MAVDT, 2002).

Colombia cuenta con un proceso reglamentado y estandarizado para la aprobacin
nacional de proyectos aplicables al MDL. Varios proyectos de generacin de energa
elctrica se han acogido a los mecanismos, entre los cuales se encuentran el Proyecto
Hidroelctrico del ro Amoy (80 MW a Filo de Agua) y el Parque Elico Jeprachi (19,5
MW instalados), el cual realiz venta de reducciones de emisiones por US $3,2 millones.

Sin embargo en el estudio de Vergara et. al. 2010, se concluye que los ingresos por
CO
2
generados dentro de los MDL no cambian significativamente la capacidad para
compensar los altos costos de inversin. Con informacin de 2007 analizada, los costos
de inversin incluso a un precio de CO
2
e de 50 dlares EE.UU., la energa elica no es
todava la opcin de menor costo. Dentro de este rango de ingresos, crditos de
carbono no afectan efectivamente la clasificacin de las opciones, lo que demuestra que
los MDL, no es suficiente para promover la energa alternativa con cero emisiones de
carbono, bajo las condiciones existentes en Colombia. Por lo tanto, las otras
opciones polticas son necesarias para facilitar el acceso al mercado de la energa elica.

4.4 Complementariedad del recurso Hdrico y Elico

Franco y Dyner (2004) describen que histricamente para Colombia, la complementacin
trmica permite el dimensionamiento econmico de los proyectos hidrulicos para las
condiciones operativas del ao hidrolgico medio, bajo la premisa de que en los aos
secos las centrales trmicas suplan la base de la demanda y las hidrulicas el pico. Por
stas y otras razones se emprendi, por parte del gobierno colombiano desde el ao
1995, la iniciativa de vincular capital privado y promover la expansin del sistema. Por
ejemplo, para el perodo de los 1993 2000, la alternativa de generacin trmica se

Ingeniero Civil

30

convirti en el eje de la expansin elctrica. En este periodo la composicin del parque
generador pas de un 80% de composicin hidrulica y 20% trmica, a un 68% y 32%
respectivamente. De tal forma que con la nueva composicin la capacidad instalada del
SIN afront el Evento ENSO - El Nio ocurrido en 1998, y si bien fue muy intenso, no
hubo racionamiento como el ocurrido en los aos 1991 -1992.

Esta estructura planteada, a pesar de ser confiable es fuertemente impactada por la
volatilidad del precio del gas, el carbn siendo adems sumamente contaminante. Por lo
anterior la importancia de incluir a la matriz energtica el componente de generacin a
partir de fuentes renovables no convencionales. Esto plantea una altsima vulnerabilidad
hacia el futuro en caso de establecerse un sistema de impuestos por produccin de CO
2
.

La complementariedad de la hidrologa con el recurso elico colombiano se ha estudiando
con anterioridad. En la Figura 12 se presenta un ejemplo de complementariedad, el caso
de la serie de caudales mensuales del ro Nare y el rgimen de vientos en La Guajira
colombiana.

Figura 12 Complementariedad hidrulica - elica entre las series de caudales medios mensuales
del ro Nare y la velocidad promedia de vientos en la Alta Guajira (Franco y Dyner, 2004).

La complementariedad hidrulica - elica puede entenderse como la oportunidad de
generar energa con una de las tecnologas, cuando se presenta escasez del recurso
energtico de la otra. Franco y Dyner, (2004) revelan, que para el caso de las Empresas
Pblicas de Medelln, sus futuros proyectos elicos y el Parque existente Jeprachi
producirn ms energa durante los perodos de disminucin de la generacin hidrulica,
en una estacin seca, evitando un mayor descenso en el nivel de sus embalses.

Sin embargo, los resultados de aprovechar la complementariedad entre los recursos
elico e hdrico, a pesar de presentar disminucin en la rentabilidad en un portafolio
conformado por las fuentes hdrica y elica, se logra reducir las variables de riesgo.
Mediante el uso de la tcnica de Dinmica de Sistemas se ha evidenciado que el valor de
riesgo de las compras de energa en bolsa se reduce de forma importante, especialmente
durante la fase clida del ENSO - El Nio, as EPM no se vera forzado a las costosas
compras en bolsa para cumplir sus contratos a largo plazo en momentos de sequa, que
afectan disminuyendo la rentabilidad de su negocio de generacin de energa (Franco y
Dyner, 2004).

Se debe aclarar que los agentes generadores de energa elctrica en el Mercado de
Energa Mayorista venden su mayor parte de energa en contratos a largo plazo, y los

Ingeniero Civil

31

excedentes generados van a la bolsa de energa (Mercado de corto plazo Spot). El
mercado colombiano tiene sus particularidades, la regulacin y la normatividad estn
ntimamente ligadas a las fuentes convencionales de energa y a las caractersticas
econmicas propias. Una de esas particularidades es la participacin estatal (pblica) en
ms del 90% de la transmisin de energa elctrica (Interconexin Elctrica S.A. - ISA) y
cerca del 35% de la generacin de energa del pas (ISAGEN y Empresas Pblicas de
Medelln - EEPPM).

La utilidad de operar en condiciones de complementariedad entre las tecnologas que
conforman un Portafolio Tecnolgico, ha sido objeto de estudio para aprovechamientos
hidrulicos - elicos en otras partes del mundo.

El trabajo de Belanger y Gagnon (2002), desarrollado en la provincia de Quebec Canad,
utiliz datos horarios reales de velocidad del viento entre 1990-1996 presenta la
simulacin de varias combinaciones de la complementariedad de los dos sistemas,
incluyendo las series de caudales y las demandas asociadas al sistema de generacin.
Los resultados no fueron alentadores, e indican que la energa elica requiere una
instalacin de respaldo adicional para compensar las fluctuaciones del recurso. Al parecer
el problema es de escala, en grandes redes y grandes tamaos de parques elicos,
puede que no sea necesario construir la capacidad adicional, ya que la confiabilidad
podra ser mantenida a partir de otras herramientas como la compra de servicios
auxiliares a otros generadores.

Esto tambin representa un costo oculto de la energa elica. En tal estudio tambin se
analiz el impacto del recurso elico en escenarios de caudales mximos y mnimos. En
Quebec, con ms de 30.000 MW de capacidad hidroelctrica, la instalacin de unos
pocos cientos de MW de energa elica no tendra un impacto significativo sobre los
almacenamientos de los embalses. Tal vez la instalacin de 2000 3000 de MW si
lograran un efecto.

Por su parte, Denault et al. (2009), presentan otro estudio realizado en Quebec con
informacin desde 1958 a 2003, utilizaron la simulacin de la operacin de un portafolio
basado en anlisis de riesgos, teniendo en cuenta tres escenarios que consideraban;
produccin total de energa hidrulica, produccin total energtica elica y combinacin
de produccin energtica. En los tres escenarios, la energa elica proporcion un efecto
sustancial de la diversificacin y una mejora del perfil de riesgo para cualquier penetracin
de la energa elica, de hasta un 30%. Sin embargo, este ltimo estudio no consider la
complementariedad estacional de los dos recursos naturales.

Para el caso de Brasil, Odilon et al. (2009) muestran que las polticas de incentivos para la
energa elica que se estn aplicando desde la ltima dcada en el pas, promovieron el
desarrollo de algunos de proyectos, principalmente a travs de las fuentes de energa
alternativas. Para el noreste de Brasil, el ro San Francisco es su principal recurso
generador de energa elctrica, dicho ro presenta disminucin significativa en sus
caudales en las pocas de verano, cuando los vientos aumentan en dicha zona. Por su
parte, en las ltimas dcadas la generacin de energa elica ha demostrado su
adecuacin a la escala de GW, necesarias para una contribucin efectiva a los sistemas
elctricos. Se ha demostrado con informacin de primera mano la complementariedad

Ingeniero Civil

32

estacional hidrologa - vientos en las reas del noreste de Brasil, y se ha analizado su
posible efecto sobre la viabilidad de la estabilizacin estacional de la oferta de energa en
la zona interconectada de dicho pas.

En el documento final del estudio del Banco Mundial, Wind Energy in Colombia A
Framework for Market Entry, se realiza un anlisis de complementariedad entre los
recursos elico e hdrico, indicando que la energa producto del viento podra contribuir al
sistema cundo este ms lo necesita, durante los perodos secos y en los picos de
demanda energtica. En dicho documento tambin se expone el hecho de que el Sistema
Interconectado Nacional podra ser afectado durante fuertes eventos El Nio dada a su
gran componente hdrico (Vergara et al, 2010).

Se utilizaron los datos existentes de generacin de energa a partir de Jeprachi (para el
perodo comprendido entre febrero 2004 - marzo 2009) y los registros de la velocidad del
viento de la estacin Puerto Bolvar se ampliaron para cubrir el perodo de 1985 a 2008
para evaluar la capacidad de generacin elica durante los perodos de sequa. El anlisis
considera cuatro ros con el desarrollo de energa hidroelctrica: Guavio, Nare, Cauca y
Magdalena. La sequa ms grave en estas cuencas corresponden con el perodo de El
Nio de abril 1991 a julio 1992, cuando ocurri racionamiento de energa, y desde abril
1997 a mayo 1998 (el segundo evento El Nio ms fuerte del siglo XX), cuando los
precios del mercado fueron muy altos, obligando a definir cambios regulatorios del
sector. Durante estos perodos la generacin estimada del viento estuvo muy por encima
del valor medio. Es decir, durante los perodos de extrema sequa relacionada con El
Nio, la energa elica en el norte de Colombia fue superior a la media ofreciendo cierta
compensacin al dficit de energa hidrulica. Lo anterior podra indicar que la
complementariedad de los recursos tambin funciona a la escala de tiempo interanual.

En el estudio de Vergara et al., (2010), la complementariedad tambin se explor
mediante el anlisis de la operacin conjunta de un sistema simple que consiste en una
operacin de parques elicos en conjunto con una planta de energa hidroelctrica de
tamao similar para cada uno de los ros estudiados y para una gama de tamaos de
embalses. Los resultados sugieren que la energa firme de la explotacin conjunta de
parques elicos y la energa hidroelctrica supera la operacin aislada de la planta de
energa hidroelctrica y del parque elico. Este resultado es vlido para una amplia gama
de tamaos posibles de embalses estudiados y de localizaciones. La fuerte
complementariedad que la operacin conjunta de parques elicos y centrales
hidroelctricas no ha sido reconocida, ni es remunerada por el sistema actual de
regulacin adoptado en Colombia.

En las Figura 13 y Figura 14 se observa que la operacin conjunta aumenta la capacidad
de almacenamiento del sistema complementario sobre todo cuando aumenta el tamao
del embalse.

En la Figura 15, se presenta que la energa firme que resulta de la operacin
conjunta del parque elico y la planta hidroelctrica hipottica es mayor que la suma de
las energas aisladas de los dos proyectos.


Ingeniero Civil

33

Figura 13 Operacin aislada y conjunta de embalse de 0,2 del ro Guavio con parque Jeprachi de
igual capacidad instalada
3
(Adaptada de Vergara et al, 2010).

Figura 14 Operacin aislada y conjunta de embalse de 0,5 del ro Guavio con parque Jeprachi de
igual capacidad instalada (Adaptada de Vergara et al, 2010).

Figura 15 Energa Firme para el ro Guavio como resultado de la Operacin Conjunta y aislada
(Adaptada de Vergara et al, 2010).

3
Volumen de embalse expresado como una fraccin del ingreso promedio de energa obtenida con Guavio

0.000
0.200
0.400
0.600
0.800
1.000
1.200
D
e
c
-
8
4
J
a
n
-
8
6
J
a
n
-
8
7
J
a
n
-
8
8
J
a
n
-
8
9
J
a
n
-
9
0
J
a
n
-
9
1
J
a
n
-
9
2
J
a
n
-
9
3
J
a
n
-
9
4
J
a
n
-
9
5
J
a
n
-
9
6
J
a
n
-
9
7
J
a
n
-
9
8
J
a
n
-
9
9
J
a
n
-
0
0
J
a
n
-
0
1
J
a
n
-
0
2
J
a
n
-
0
3
J
a
n
-
0
4
J
a
n
-
0
5
J
a
n
-
0
6
J
a
n
-
0
7
J
a
n
-
0
8
R
e
s
e
r
v
o
i
r

v
o
l
u
m
e
/
R
e
s
r
v
o
i
r

c
a
p
a
c
i
t
y
GUAVIO RIVER RESERVOIR OPERATION
Reservoir Size =.2
Isolated
Joint
0.000
0.200
0.400
0.600
0.800
1.000
1.200
D
e
c
-
8
4
J
a
n
-
8
6
J
a
n
-
8
7
J
a
n
-
8
8
J
a
n
-
8
9
J
a
n
-
9
0
J
a
n
-
9
1
J
a
n
-
9
2
J
a
n
-
9
3
J
a
n
-
9
4
J
a
n
-
9
5
J
a
n
-
9
6
J
a
n
-
9
7
J
a
n
-
9
8
J
a
n
-
9
9
J
a
n
-
0
0
J
a
n
-
0
1
J
a
n
-
0
2
J
a
n
-
0
3
J
a
n
-
0
4
J
a
n
-
0
5
J
a
n
-
0
6
J
a
n
-
0
7
J
a
n
-
0
8
R
e
s
e
r
v
o
i
r

v
o
l
u
m
e
/
R
e
s
r
v
o
i
r

c
a
p
a
c
i
t
y
GUAVIO RIVER RESERVOIR OPERATION
Reservoir Size =.5
Isolated
Joint
0.000
0.200
0.400
0.600
0.800
1.000
1.200
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
F
i
r
m

e
n
e
r
g
y
/
M
e
a
n

e
n
e
r
g
y
Reservoir size/Mean annual energy
FIRM ENERGY (Fraction of mean energy)
Guavio River
Guavio R.(isolated)
Jepirachi (isolated)
Guavio R. + Jepirachi
(isolated)
Guavio R. + Jepirachi
(joint)

Ingeniero Civil

34

4.5 Elementos tericos para el anlisis de informacin

Para el procesamiento y anlisis de la informacin se utilizaron conceptos bsicos de la
Estadstica y Anlisis Espectral de series de tiempo. Del libro de Montgomery y Runger,
1997, se tomaron los conceptos de: Desviacin Estndar, Error Estndar de la Media
definidos por:

Donde v
i
son las velocidades del viento a la escala considerada.

Donde es la desviacin estndar, v
i
son las velocidades del viento a la escala
considerada y es la media de los datos de viento considerados.

Para el anlisis espectral se aplic la FFT, una manera ms rpida de calcular la
Transformada Discreta de Fourier, la cual, bsicamente, descompone una seal en senos
y cosenos de diferentes frecuencias y con diferentes amplitudes. La transformada discreta
de Fourier (DFT) de una serie de datos (registro) se define como:

Con n = -N/2+1+ N/2, donde:
con
representa una
serie de N Datos en el tiempo, separados entre s , y
es el resultado de la DFT para

el ciclo asociado para la frecuencia dada por n/N (Meja et al, 1999).

Para encontrar la relacin entre las variables de inters se utiliz el coeficiente de
correlacin de Pearson, el cual es un ndice que mide la relacin lineal entre dos variables
aleatorias cuantitativas. A diferencia de la covarianza, la correlacin de Pearson es
independiente de la escala de medida de las variables.

En este caso cuando se presenta el anlisis para dos variables aleatorias x e y; el
coeficiente de correlacin de Pearson se simboliza con la letra
x,y
, siendo la expresin
que nos permite calcularlo:

Donde:
XY
es la covarianza de (X,Y)
X
es la desviacin tpica de la variable X
Y
es la desviacin tpica de la variable Y

Para los clculos presentados en el documento el coeficiente se calcul a partir de la
funcin de Microsoft Excel COEF.DE.CORREL(matriz1; matriz2) (tomado de ayuda

Ingeniero Civil

35

Microsoft Excel 2010). La sintaxis de la funcin COEF.DE.CORREL tiene los siguientes
argumentos (argumento: valor que proporciona informacin a una accin, un evento, un
mtodo, una propiedad, una funcin o un procedimiento.):

Matriz 1 Variable 1. Un rango de celdas de valores.
Matriz 2 Variable 2. Un segundo rango de celdas de valores

En algunos casos se presenta la correlacin de la informacin de dos variables
analizadas, mediante representacin grfica de regresin lineal de los datos de cada
variable; el software realiza el ajuste por mnimos cuadrados y cuantifica dicho ajuste
presentando la ecuacin lineal y el denominado coeficiente de determinacin (r
2
), el
coeficiente de correlacin ser entonces la raz cuadrada del r
2
.

Para cuantificar o establecer el grado de complementariedad se utilizar el coeficiente de
correlacin. El coeficiente ser positivo siempre y cuando las variables sean directamente
proporcionales (al aumentar la variable x aumenta la variable y) y negativo en caso de
comportamiento inversamente proporcional. La complementariedad estar asociada a
comportamiento inverso en este caso de las variables de caudal y viento.

Para los clculos de coeficientes de correlacin se consider un intervalo de confianza
mnimo del 95%.

Para las matrices de correlacin trimestral, se realizaron promedios trimestrales de las
variables analizadas correlacionndolas de manera cruzada y abarcando los cuatro (4)
trimestres del ao as:

Matriz de
correlaciones
trimestrales
genrica
Variable 1
DEF
Diciembre
Enero
Febrero
MAM
Marzo
Abril
Mayo
JJA
Junio
Julio
Agosto
SON
Septiembre
Octubre
Noviembre
Variable 2
DEF
MAM
JJA
SON


Ingeniero Civil

36

5 ANLISIS DEL RECURSO ELICO DE LA GUAJIRA COLOMBIANA

Como se mencion anteriormente, la Guajira colombiana cuenta con un alto potencial con
respecto a otras regiones del pas. Para el anlisis de la velocidad del viento se utiliz
informacin de varias estaciones climatolgicas de la zona. En principio se analiz el
comportamiento de la velocidad y direccin del viento. La informacin utilizada y
debidamente procesada incluye cuatro (4) estaciones climatolgicas del IDEAM y tres (3)
estaciones climatolgicas especiales de ISAGEN, ver Tabla 1.

Tabla 1 Informacin de estaciones de registro de velocidad del viento utilizadas
Estacin
Localizacin
Latitud
Longitud
Municipio -
Departamento
Perodo de
informacin
Cant. Unidad
Informacin
Faltante (%)
Variables Resolucin
Puerto
Bolvar
1211'N
7155' W
Uribia - La
Guajira
1986 - 2006 20 ao 21
Velocidad y
direccin del viento
- temperatura
horaria -
mensual
Almirante
Padilla
1132' N
7256' W
Riohacha - La
Guajira
1975 - 1984 9 ao 40
Velocidad y
direccin del viento
horaria
Parque
Elico
1212'36.60"N
721'20.10"W
Uribia - La
Guajira
2008 - 2010 18 mes 3
Velocidad y
direccin del viento
10 minutos
Nazaret
1211' N
7117' W
Uribia - La
Guajira
1974 - 1976 3 ao 50
Velocidad y
direccin del viento
horaria
La Mina
1107' N
7236' W
Hatonuevo - La
Guajira
1988 - 1991 4 ao 47
Velocidad y
direccin del viento
horaria
Virtual
NCEP/NCAR
124'60.00" N
725'0.00" W
Uribia - La
Guajira
1948 - 2009 61 ao 0
Velocidad del
viento
mensual
Flor de la
Sabana
1131'34.60" N
72 2'37.56" W
Maicao - La
Guajira
2009 4 mes 2
Velocidad y
direccin del viento
10 minutos
Manaure
1148'30.80" N
7223'18.90" W
Manaure - La
Guajira
2009 4 mes 2
Velocidad y
direccin del viento
10 minutos

Las estaciones del IDEAM capturan datos horarios de direccin y velocidad del viento en
superficie a una altura de 10 m y cuentan con veletas a la misma altura para registrar la
direccin del viento. Los datos de direccin del viento se encuentran discretizados en 8
sectores: N, NE, E, SE, S, SW, W, NW. Se har nfasis en la estacin Puerto Bolvar
dado que tiene el mayor perodo de registros y adicionalmente cuenta con datos de
precipitacin y temperatura a nivel horario para cada variable.

Las estaciones especiales de ISAGEN registran informacin de velocidad y direccin del
viento, as como temperatura, humedad relativa y presin atmosfrica a 20, 40 y 60
metros de altura, cada 10 minutos. La direccin del viento es registrada a una resolucin
de grados (azimut). En la Figura 16 se presenta el resumen de las estaciones, periodo de
registro, ubicacin, y dems informacin relevante.

Para los anlisis regionales se utiliz informacin mensual de un punto en la malla de los
resultados del Proyecto del Reanlisis Climtico de National Center for Enviromental
Prediction - NCEP/National Center for Atmospheric Research-NCAR,
(www.cdc.noaa.gov/data/gridded/data.narr.html) en la coordenada 124'60.00" N - 72
5'0.00" W, ubicado cerca de la estacin Puerto Bolvar.

Adicionalmente, se utiliz informacin pblica y de libre acceso de energa generada por
el Parque Jeprachi disponible en la base de datos Nen de XM (www.xm.com.co). Este

Ingeniero Civil

37

Parque cuenta con 15 aerogeneradores de 1,3 MW con torres a una altura de 60 m
medidos hasta el eje de la mquina. Para obtener la informacin de velocidad horaria del
viento, se utilizaron los valores de la curva caracterstica de la mquinas,
(Aerogeneradores Nordex-N60/1.300 kW), se dibuj la curva de potencia y se identificaron
las ecuaciones de ajuste a los datos tanto para la funcin original que relaciona la
velocidad del viento (ms
-1
) con la Potencia (MW), es decir en MW vs. ms
-1
.


Ingeniero Civil

38

0,0
2,0
4,0
6,0
8,0
10,0
12,0
0
7
:
0
0
0
9
:
0
0
1
1
:
0
0
1
3
:
0
0
1
5
:
0
0
1
7
:
0
0
1
9
:
0
0
2
1
:
0
0
2
3
:
0
0
0
1
:
0
0
0
3
:
0
0
0
5
:
0
0
V
e
l
o
c
i
d
a
d
(
m
/
s
)
Hora
Estacin Parque Elico
Ciclo Diurno Promedio - Velocidad del Viento
Parque Elico
0,0
2,0
4,0
6,0
8,0
10,0
12,0
0
7
:
0
0
0
9
:
0
0
1
1
:
0
0
1
3
:
0
0
1
5
:
0
0
1
7
:
0
0
1
9
:
0
0
2
1
:
0
0
2
3
:
0
0
0
1
:
0
0
0
3
:
0
0
0
5
:
0
0
V
e
l
o
c
i
d
a
d

(
m
/
s
)
Hora
Estacin Puerto Bolivar
0,0
2,0
4,0
6,0
8,0
10,0
12,0
0
7
:
0
0
0
9
:
0
0
1
1
:
0
0
1
3
:
0
0
1
5
:
0
0
1
7
:
0
0
1
9
:
0
0
2
1
:
0
0
2
3
:
0
0
0
1
:
0
0
0
3
:
0
0
0
5
:
0
0
V
e
l
o
c
i
d
a
d

d
e
l

V
i
e
n
t
o

(
m
/
s
)
Hora
Estacin La Mina
0,0
2,0
4,0
6,0
8,0
10,0
12,0
0
7
:
0
0
0
9
:
0
0
1
1
:
0
0
1
3
:
0
0
1
5
:
0
0
1
7
:
0
0
1
9
:
0
0
2
1
:
0
0
2
3
:
0
0
0
1
:
0
0
0
3
:
0
0
0
5
:
0
0
V
e
l
o
c
i
d
a
d

d
e
l

V
i
e
n
t
o

(
m
/
s
)
Hora
Estacin Nazareth
Figura 16 Ciclo diurno promedio de la velocidad del viento para la zona de estudio

Estacin Manaure

Estaciones IDEAM. Anemmetros y
Veletas a 10 metros.

Estaciones Especiales ISAGEN S.A.
Anemmetros y Veletas a 20, 40 y 60
metros. Medicin Temperatura.

Parque Elico Jeprachi. Conversin de
Velocidades a partir de Potencias.

Estacin Virtual Datos Reanlisis

Estacin Flor de la Sabana

0,0
2,0
4,0
6,0
8,0
10,0
12,0
0
7
:
0
0
0
9
:
0
0
1
1
:
0
0
1
3
:
0
0
1
5
:
0
0
1
7
:
0
0
1
9
:
0
0
2
1
:
0
0
2
3
:
0
0
0
1
:
0
0
0
3
:
0
0
0
5
:
0
0
V
e
l
o
c
i
d
a
d
(
m
/
s
)
Hora
Estacin Almirante Padilla
Puerto Bolivar 20 aos
Nazareth 3 aos
Apto. AlmirantePadilla 9 aos
La Mina 4 aos
Intercor 2 aos
Flor de la Sabana 4 meses
Manaure 4 meses

Ingeniero Civil

39

Del procesamiento de la informacin de direccin y magnitud de la velocidad viento en la
zona del Departamento de la Guajira, se identific la variacin intradiaria (Ciclo Diurno), el
comportamiento Mensual Intranual (Ciclo Anual) y la Rosa de Vientos.

La variabilidad intradiurna de la velocidad del viento est asociada con la temperatura del
aire, resultante del ciclo de insolacin superficial en las 24 horas del da, fuertemente
marcado en la zona tropical (Poveda, 2004). Los valores ms bajos de la velocidad del
viento ocurren en horas la madrugada entre las 03:00 y 05:00 justo antes de salir el sol, y
los mximos al final de la tarde antes de que la temperatura de la superficie comience a
debilitarse (entre las 15:00 y 16:00 horas), tal como se observa en las estaciones Flor de
la Sabana y Manaure, las cuales registran la temperatura de la zona. Ver Figura 16,
Figura 17 y Figura 18.

Figura 17 Comparacin ciclo diurno velocidad del viento ciclo diurno temperatura en la estacin
Manaure.

Figura 18 Comparacin ciclo diurno velocidad del viento ciclo diurno temperatura en la estacin
Flor de la Sabana.

Por otro lado, la variacin intranual de la velocidad del viento (ciclo anual) presentada en
la Figura 19 est relacionado entre otros, con los perodos de precipitacin y sequa
directamente influenciadas por la Zona de Convergencia Intertropical - ZCIT, la cual acta
20
21
22
23
24
25
26
27
28
29
30
0,0
2,0
4,0
6,0
8,0
10,0
12,0
14,0
0
7
:
0
0
0
9
:
0
0
1
1
:
0
0
1
3
:
0
0
1
5
:
0
0
1
7
:
0
0
1
9
:
0
0
2
1
:
0
0
2
3
:
0
0
0
1
:
0
0
0
3
:
0
0
0
5
:
0
0
T
e
m
p
e
r
a
t
u
r
a

(
C
)
V
e
l
o
c
i
d
a
d

(
m
/
s
)
Hora
Estacin Manaure
Relacin Ciclos Diurnos Promedios - Velocidad del Viento y Temperatura
60m 40m
20m T
20
21
22
23
24
25
26
27
28
29
30
0,0
2,0
4,0
6,0
8,0
10,0
12,0
14,0
0
7
:
0
0
0
9
:
0
0
1
1
:
0
0
1
3
:
0
0
1
5
:
0
0
1
7
:
0
0
1
9
:
0
0
2
1
:
0
0
2
3
:
0
0
0
1
:
0
0
0
3
:
0
0
0
5
:
0
0
T
e
m
p
e
r
a
t
u
r
a

(
C
)
V
e
l
o
c
i
d
a
d
(
m
/
s
)
Hora
Relacin Ciclos Diurnos Promedios - Velocidad del Viento y Temperatura
60m 40m
20m T

Ingeniero Civil

40

sobre todo el territorio nacional (Poveda, 2004). Adicional a esto, dada la posicin de la
zona de estudio (norte de Colombia) la ZCIT solo pasar una vez al ao, lo cual se refleja
en un ciclo anual unimodal en la velocidad del viento y la precipitacin promedios
mensuales en la zona. De igual manera, existe una relacin entre el aumento de la
velocidad del viento, producto del aumento de la temperatura y, aunque la variabilidad sea
mayor en la velocidad, las dos variables presentan mximos justo antes de empezar el
perodo lluvioso, este comportamiento se aprecia en los meses de junio, julio y agosto.
Ver Figura 21 y Figura 22.
En la Figura 19 se observa que las mayores velocidades se presentan en el sitio
denominado Cabo de la Vela, corroborado con la estacin Puerto Bolvar y la estacin
Virtual estimada a partir del Reanlisis NCEP/NCAR. La diferencia en magnitud de
velocidades se debe a que las otras estaciones estn afectadas por el entorno de su
emplazamiento: para el caso de la estacin Nazaret la disminucin en la magnitud de la
velocidad podra deberse a la barrera orogrfica que impone la Serrana La Macuira. Por
otro lado, las estaciones Almirante Padilla y La Mina, se ven afectadas, la primera por
estar ubicada dentro del Aeropuerto (zona urbana) y la segunda dentro de la Mina del
Cerrejn (campo INTERCOR). Sin embargo, es clara la semejanza en el comportamiento
del ciclo anual de los vientos en las distintas zonas de la Guajira colombiana.
Es necesario conocer el comportamiento del ciclo anual de la velocidad, con el objeto de
verificar si existe alguna correlacin entre los perodos de verano en el centro y norte del
Pas. Considerando el desplazamiento de la ZCIT por el territorio nacional, se observa
que los perodos secos del centro de pas, diciembre - enero febrero y junio julio
agosto (donde se encuentran la mayora de los embalses de generacin del SIN), podran
ser respaldados por los mayores vientos en La Guajira, dado que se tiene una poca
comn de verano correspondiente al trimestre junio julio agosto.
En relacin a la direccin de los vientos en la Guajira, se estimaron las Rosas de Viento
para todas las estaciones de medicin. Se encontr un altsimo porcentaje de la direccin
de la velocidad del viento dentro del Cuadrante NE E, correspondiente a los 45 90 de
Azimut. Esto indica la predominancia de la componente zonal del viento, claramente
explicada por la presencia de los vientos Alisios del este. El anlisis no revel cambios
significativos intranuales en la direccin del viento, sin embargo, en la Figura 20 se
aprecia que a medida que nos trasladamos hacia la zona sur de la Guajira la componente
meridional toma mayor importancia.


Ingeniero Civil

41

Figura 19 Ciclo anual promedio de la velocidad del viento para la zona de estudio. Las barras de
error representan el error estndar de estimacin de la media mensual respectiva.
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
8,0
9,0
10,0
Jun Jul Ago Sep Oct Nov Dic Ene Feb Mar Abr May
V
e
l
o
c
i
d
a
d

(
m
/
s
)
Mes
Estacin Almirante Padilla
Ciclo Anual Promedio - Velocidad del Viento
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
8,0
9,0
10,0
V
e
l
o
c
i
d
a
d

(
m
/
s
)
Mes
Estacin La Mina
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
8,0
9,0
10,0
V
e
l
o
c
i
d
a
d

(
m
/
s
)
Mes
Estacin Virtual Reanlisis
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
8,0
9,0
10,0
V
e
l
o
c
i
d
a
d
(
m
/
s
)
Mes
Estacin Nazareth
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
8,0
9,0
10,0
V
e
l
o
c
i
d
a
d

(
m
/
s
)
Mes
Estacin Puerto Bolivar

Ingeniero Civil

42

Figura 20 Variabilidad de las rosas de viento en la zona de estudio

0%
5%
10%
15%
20%
25%
30%
N
NE
E
SE
S
SW
W
NW
Estacin La Mina
Rosa de Vientos Multianual
0.1-1.5 1.6-3.3 3.4-5.4
5.5-7.9 >=8
0%
5%
10%
15%
20%
N
NE
E
SE
S
SW
W
NW
Estacin Apto Almirante Padilla
0.1-1.5 m/s 1.6-3.3 m/s 3.4-5.4 m/s
5.5-7.9 m/s >=8.0 m/s
0%
5%
10%
15%
20%
25%
30%
N
NE
E
SE
S
SW
W
NW
Estacin Nazareth
0.1-1.5 m/s 1.6-3.3 m/s 3.4-5.4 m/s
5.5-7.9 m/s >=8 m/s
0%
15%
30%
45%
60%
75%
N
NE
E
SE
S
SW
W
NW
Estacin Puerto Bolvar
>= 7 m/s 3.6 -7 m/s 0.1-3.5 m/s

Ingeniero Civil

43

Figura 21 Ciclo anual promedio de la temperatura media mensual y de la velocidad del viento en la
estacin Puerto Bolvar.

Figura 22 . Ciclo anual promedio de la temperatura media mensual y de la velocidad del viento en
la estacin Puerto Bolvar.

Para la estacin Puerto Bolvar se estimaron las anteriores periodicidades del viento a dos
escalas de tiempo consideradas; ciclo diario y ciclo anual, mediante el uso de la
Transformada Rpida de Fourier FFT, definida como:
1
0
2
1
) (
N
k
N
ink
i n
e f
N
t f
t

con n = -N/2+1+ N/2, Donde:
i i
ft f con t i t
i
A = 1 , , 1 , 0 = N i , representa una
serie de N datos en el tiempo, separados entre s t A , y
n
f
es el resultado de la DFT
(Transformada Discreta de Fourier) para el periodo asociado con la frecuencia dada por
n/N t A (Meja et al., 1999).
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
8,0
9,0
20
22
24
26
28
30
32
V
e
l
o
c
i
d
a
d

d
e
l

V
i
e
n
t
o

(
m
/
s
)
T
e
m
p
e
r
a
t
u
r
a

(
C
)
Mes
Ciclos Promedios Anuales: Velocidad del Viento y Temperatura
Temperatura
Velocidad
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
8,0
9,0
0,0
20,0
40,0
60,0
80,0
100,0
120,0
V
e
l
o
c
i
d
a
d

d
e
l

V
i
e
n
t
o

(
m
/
s
)
P
r
e
c
i
p
i
t
a
c
i
n

(
m
m
)
Mes
Ciclos Promedios Anuales: Velocidad del Viento y Precipitacin
Precipitacin
Velocidad

Ingeniero Civil

44

En la Figura 23 y Figura 24 se observa picos espectrales importantes en los perodos de
12 y 24 horas, as como para 6 y 12 meses representando respectivamente el semiciclo y
ciclo anual de la velocidad del viento.

Figura 23 Anlisis espectral de la velocidad promedia mensual del viento en la estacin Puerto
Bolvar.

Figura 24 Anlisis espectral de la velocidad promedia horaria del viento en la estacin Puerto
Bolvar.

Adicionalmente, se estudi la variacin de la velocidad del viento con respecto a la altura.
Para las estaciones de Flor de la Sabana y Manaure se elaboraron perfiles verticales de
velocidad del viento. La velocidad vara con la altura, siguiendo aproximadamente una
ecuacin conocida como ley exponencial de Hellmann, de la forma:
o
|
|
.
|
\
|
=
0
0
H
H
v v
,

donde:
0
1
2
3
4
5
6
7
8
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
D
e
n
s
i
d
a
d

e
s
p
e
c
t
r
a
l
Periodicidad (meses)
Anlsis Espectral de Velocidad del Viento Promedio Mensual
0
5
10
15
20
25
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100
D
e
n
s
i
d
a
d

e
s
p
e
c
t
r
a
l
Periodicidad (Horas)
Anlisis Espectral de Velocidad del Viento Promedio Horaria

Ingeniero Civil

45

v
: es la velocidad (m s
-1
) a la altura
H
(m)
0
v
: es la velocidad (m s
-1
) a una altura
0
H de referencia, normalmente referida a una
altura de 10 m.
o
: es el coeficiente de friccin o exponente de Hellmann (el cual vara con la topografa
del terreno).
La relacin funcional entre la velocidad con la altura se presentan en la Figura 25 y Figura
26, para los valores promedios, mximos y mnimos diarios.

Figura 25 - Variacin de la velocidad del viento con respecto a la altura en la Estacin Flor de la
Sabana.

Figura 26 - Variacin de la velocidad del viento con respecto a la altura en la estacin Manaure.

En las Figura 25 y Figura 26 se observa que la velocidad del viento vara segn la ley
exponencial de Hellmann corroborando el comportamiento logartmico de la variacin del

0
10
20
30
40
50
60
70
0,0 2,0 4,0 6,0 8,0 10,0 12,0 14,0
A
l
t
u
r
a

s
o
b
r
e

l
a

s
u
p
e
r
f
i
c
i
e

(
m
)
Velocidad del viento (m/s)
Variacin de la Velocidad del Viento con la altura
Prom Diario
Max Diario
Min Diario

0
10
20
30
40
50
60
70
0,0 2,0 4,0 6,0 8,0 10,0 12,0 14,0
A
l
t
u
r
a

s
o
b
r
e

l
a

s
u
p
e
r
f
i
c
i
e

(
m
)
Velocidad del Viento (m/s)
Estacin Manaure
Variacin de la Velocidad del Viento con la altura
Promedio Diario
Mx Diario
Mn Diario

Ingeniero Civil

46

viento en la capa lmite atmosfrica. Para la estacin Manaure se calcul un coeficiente
o
= 0,19, y para Flor de la Sabana, el valor de
o
corresponde a 0,20.
La importancia de analizar los perfiles de velocidad radica en que las palas (aspas) de los
aerogeneradores estarn sujetas a diferentes cargas durante la operacin debido a que
normalmente el rea abarcada por las palas comprende un crculo con dimetros entre los
30 m y los 100 m, es decir cuando la pala se encuentre en la parte superior recibir mayor
carga de viento que las dos palas restantes.
Este estudio no pretende analizar las diferentes cargas de viento a las cuales est
sometido un aerogenerador, solo se presenta como una referencia relevante, y los
factores calculados solo servirn para cuantificar la velocidad del viento a la altura
requerida. En la Figura 27 se presentan las caractersticas tpicas de un aerogenerador
convencional de eje horizontal y tres aspas, se presenta la direccin del viento y la
inclinacin producto del aumento de la velocidad con altura mencionada anteriormente.

Figura 27 Componentes bsicos de una turbina elica moderna, de eje horizontal con una caja de
cambios (Diseo por el Laboratorio Nacional de Energa Renovable - NREL). Adaptada de SRREN IPCC, 2007.
Otro punto importante es que a partir del perfil de velocidades obtenido es posible
identificar el lmite de la capa lmite atmosfrica en un punto de la vertical en el que no hay
cambios significativos de velocidad con la altura. Ese punto indica posiblemente la altura
mxima a la que pueden instalarse aerogeneradores, dado que para alturas mayores no
aumentar la produccin energtica. Para los anlisis realizados en este trabajo se
utilizar el aerogenerador de eje horizontal y operacin a barlovento Nordex N60/1300
con una altura de buje de 60 m.

Ingeniero Civil

47

5.1 Variabilidad Interanual de la velocidad del viento en La Guajira: Impacto del
Evento ENSO (El Nio - Oscilacin del Sur).

A escala Interanual se requiere de informacin de estaciones con un perodo de registro
con al menos 20 aos, la cual es representativa dado que incluye varios Eventos ENSO.
Por ello, solo se utiliz la estacin Puerto Bolvar, que cuenta con cerca de 26 aos de
informacin horaria disponible.
En la Figura 23 para el anlisis espectral mensual, se observa un pico para los 55 meses
(aproximadamente 4 aos), muy probablemente explicado por el ENSO, y otro pico de
menor potencia entre los 25-30 meses, posiblemente asociado a la Oscilacin Cuasi-
Bienal, la cual no es objeto de ste anlisis (Poveda, 2004).
Para examinar la influencia del ENSO sobre el ciclo anual de los vientos de la Guajira
Colombiana se cuantifica el grado de dependencia lineal existente entre dichas variables,
mediante la estimacin de correlaciones para diferentes rezagos. Para el anlisis se
estandariz la serie de velocidad del viento y se compararon con las anomalas
estandarizadas del Oceanic Nio Index ONI, como se muestra en la Figura 28.

Figura 28 Evolucin simultnea de las series mensuales estandarizadas de la velocidad del viento
en la estacin Puerto Bolvar con el ndice Ocenico del Nio (ONI)
4

El coeficiente de correlacin simultneo (rezago 0) entre las dos series es de 0,26 con un
intervalo de confianza del 95%. Se calcularon los coeficientes hasta 5 meses de rezago
del ndice con respecto a la velocidad del viento sin aportar mejores resultados. La poca
correlacin podra estar relacionada con la estandarizacin de la serie, la cual remueve el
ciclo anual de la variable.

4
Los valores de las anomalas se extrajeron del sitio Web
www.cpc.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml, consultada en diciembre de 2010.

-4
-3
-2
-1
0
1
2
3
1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005
Relacin entre la serie estandarizada de velocidad del viento y el ONI
ONI
Velocidad Mensual

Ingeniero Civil

48

Con base en los resultados anteriores, se procedi a calcular los coeficientes de
correlacin trimestrales y mensuales de las series de los ndices ONI y SOI con la serie no
estandarizada de las velocidades del viento (se consideraron aos hidrolgicos, es decir
desde junio del ao i-1 a mayo del ao i). Los resultados de las correlaciones trimestrales
cruzadas se presentan en las Tabla 2, Tabla 3 y Tabla 4.
Tabla 2 Matriz de correlaciones trimestrales entre el ndice ONI y la velocidad del viento en la
estacin Puerto Bolvar

Tabla 3 Matriz de correlaciones trimestrales entre el ndice SOI y la velocidad del viento en la
estacin Puerto Bolvar

De la Tabla 2 y Tabla 3 se observa que las mayores correlaciones se presentan para el
ndice ONI en el trimestre SON con las velocidades en los trimestres JJA y SON.
Comportamiento igual presenta el ndice del SOI. Los signos de las correlaciones indican
su asociacin inversa. El ONI representa el aumento o disminucin de la temperatura
superficial del mar y de acuerdo con este ndice, en condiciones El Nio (La Nia), el ONI
debe ser igual superior (igual o inferior) a medio grado Celsius de temperatura por
encima o por debajo de la media (anomala). Por lo tanto se considera en este caso, si las
variables son directamente proporcionales, que al aumentar el ndice ONI se presentarn
aumentos en la velocidad del viento y viceversa.
Contrario al ONI, el ndice de Oscilacin del Sur SOI es de tipo atmosfrico. Se calcula
como las diferencias normalizadas entre la presin atmosfrica medida en Tahit
(Polinesia Francesa, representa el sector central del Pacfico tropical) y Darwin (norte de
Australia, representa el sector occidental). Es un indicador de las fluctuaciones de la masa
atmosfrica entre los sectores centro-oriental y occidental del pacfico. Cuando este ndice
es positivo, es posible establecer la ocurrencia de un evento La Nia (o enfriamiento del
ocano), pero cuando el ndice es negativo durante varios meses consecutivos, es posible
que se est dando un evento El Nio. Por lo anterior cuando se presenta una disminucin
en el SOI se presentarn valores mayores de velocidad del viento, siempre y cuando
estas dos variables sean inversamente proporcionales.
Para el ONI, los resultados de las correlaciones son del orden de 0,47 con una
confiabilidad del 96%. Para los trimestres en los cuales el ciclo anual presenta las
mayores velocidades del viento, es decir simultneamente en JJA-JJA y SON-SON y
DEF MAM JJA SON
DEF 0,24 0,19 -0,02 -0,03
MAM 0,17 0,17 0,24 0,31
JJA 0,06 0,08 0,43 0,47
SON 0,14 0,15 0,50 0,47
ONI
DEF MAM JJA SON
DEF -0,35 -0,29 -0,35 -0,15
MAM -0,28 -0,22 -0,06 -0,43
JJA -0,23 -0,27 -0,45 -0,45
SON -0,24 -0,32 -0,54 -0,58
SOI

Ingeniero Civil

49

antecediendo un trimestre el ndice a los viento (JJA-SON). En concordancia, el SOI se
present una correlacin mayor 0,58 con una confiabilidad del 95%.
Los datos del Reanlisis (nodo ubicado a 20 km de la estacin Puerto Bolvar)
presentaron el mismo comportamiento; adems de tener una correlacin de 0,66 con los
datos de vientos en la estacin de Puerto Bolvar, se encontraron resultados similares con
relacin a la influencia del ENSO, como se muestra en la Tabla 4. Las mayores
correlaciones se tienen para el ONI en JAA asociado a velocidades del viento en JJA y
SON-JJA.
Tabla 4 Matriz de correlaciones trimestrales entre el ndice ONI y la velocidad del viento Nodo
Reanlisis.

5.1.1 Variacin del ciclo anual de la velocidad del viento con el ENSO

Para este anlisis se utilizaron los datos a nivel mensual de la estacin Puerto Bolvar, la
estacin virtual del Reanlisis y, como apoyo, los datos de la estacin del aeropuerto
Ernesto Cortizzos de la ciudad de Barranquilla. El anlisis se realiz promediando las
velocidades del viento durante eventos de El Nio y La Nia, definidos segn la
Administracin Nacional Ocenica y Atmosfrica (NOAA), de Estados Unidos
(www.cdc.noaa.gov/data/gridded/data.narr.html) y comparndolos con las condiciones de
aos normales.
En la Figura 29 se observa el impacto positivo de El Nio en la mayora de los meses, con
valores de velocidad del viento por encima de las condiciones normales. Caso contrario
se presenta durante La Nia, particularmente entre agosto y marzo del ao hidrolgico.

DEF MAM JJA SON
DEF -0,10 0,34 - 0,04 - 0,20
MAM -0,10 0,43 - 0,19 0,30
JJA 0,05 0,21 - 0,46 0,33
SON 0,06 0,12 - 0,51 0,33
ONI
4
5
6
7
8
9
10
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
V
e
l
o
c
i
d
a
d
(
m
/
s
)
Mes
Nia Nio Normal

Ingeniero Civil

50

Figura 29 Estacin Puerto Bolvar, comparacin de ciclos anuales de velocidad del viento, aos El
Nio, La Nia y normales.

Para el nodo virtual del Reanlisis se presenta un comportamiento similar para la segunda
mitad del ao, mayores velocidades del viento, sin embargo, el primer semestre es
totalmente contrario, entre enero y junio las velocidades del viento aumentan durante La
Nia y disminuyen durante El Nio.

Figura 30 Estacin virtual Reanlisis, comparacin de ciclos anuales de velocidad del viento, aos
El Nio, La Nia y normales.
Cabe anotar que la estacin virtual tiene mayor perodo de registros, desde enero 1960
hasta abril de 2009. Este perodo comprende un mayor nmero de ocurrencias de ambas
fases del evento ENSO; sin embargo, los valores encontrados indican que la ocurrencia
de El Nio no solo disminuira las velocidades del viento en los meses de enero a mayo,
sino que en estos es cuando se presentan disminuciones de los niveles de los embalses
en el centro del pas.
Resultados similares se presentan en el trabajo Ruz y Bernal (2009), donde a partir de
informacin del Reanlisis NCEP/NCAR e ICOADS identificaron la variacin del ciclo
anual de los vientos, en los trimestres DEF y MAM la velocidad del viento disminuye
(aumentan las frecuencias de vientos dbiles) y mayor durante La Nia (aumentan las
frecuencias de vientos fuertes). En cambio, en los trimestres JJA y SON la velocidad del
viento se intensifica durante El Nio y disminuye durante La Nia, esto se debe a que el
chorro tropical superficial del Caribe, el cual domina la variacin espacio temporal de los
vientos de la zona, obedece al forzamiento conjunto de las anomalas clidas en el
Pacfico oriental tropical y de las anomalas de presin en el alto del Atlntico Norte
(teleconexiones).
Con el fin de verificar el comportamiento de la velocidad del viento entre los meses de
enero y mayo ante la ocurrencia del Evento ENSO y dado a que las estaciones de Puerto
Bolvar y Reanlisis arrojan resultados inversos (es decir, velocidades del viento menores
y mayores durante El Nio en los meses de enero a mayo, Reanlisis y Puerto Bolvar
respectivamente), se recurri a la estacin del aeropuerto Ernesto Cortizzos, la cual se
encuentra ubicada en la Costa Caribe colombiana y cuenta con un perodo de registros
4
5
6
7
8
9
10
V
e
l
o
c
i
d
a
d

(
m
/
s
)
Mes
Nia Normal Nio

Ingeniero Civil

51

que cubre varios Eventos ENSO. En este caso, a pesar de existir una notoria diferencia
en el ciclo anual, presentando los mayores valores de velocidad entre enero y mayo, y un
marcado comportamiento unimodal, el efecto del evento ENSO se refleja de forma similar
que para el caso de los datos de la estacin virtual del Reanlisis NCEP/NCAR.

Figura 31 Estacin Ernesto Cortizzos, comparacin de ciclos anuales de velocidad del viento, aos
El Nio, La Nia y normales.

Particularmente en la Figura 31 se observa que durante La Nia, la velocidad del viento
aumenta entre enero y junio en un 4% con relacin a los valores promedios. Utilizando los
resultados arrojados por la estacin Ernesto Cortizzos, Reanlisis y Puerto Bolvar se
infiere que entre los meses de julio y diciembre en promedio, la velocidad aumenta
durante la ocurrencia de El Nio; con los datos de reanlisis se estim que en promedio la
velocidad aumenta un 8% con relacin a los valores medios de velocidad, incluso en los
meses de agosto, septiembre y noviembre el aumento es superior al 12%.

5.1.2 Ciclo anual de la temperatura superficial del aire y su relacin con la
velocidad el viento y el ENSO
Un anlisis adicional realizado a partir de la informacin estandarizada de la temperatura
media mensual disponible (en C), indic una mayor correlacin para entre un rezago de 2
a 4 meses de la temperatura con respecto al ndice ONI y del orden de 0,49 con una
confiabilidad del 95%. Ver Figura 32 y Figura 33.

2
3
4
5
6
7
V
e
l
o
c
i
d
a
d

(
m
/
s
)
Mes
ENSO - La Nia ENSO - El Nio Aos Normales
-3
-2
-1
0
1
2
3
1
9
8
8
1
9
8
8
1
9
8
8
1
9
8
9
1
9
8
9
1
9
9
0
1
9
9
0
1
9
9
0
1
9
9
1
1
9
9
1
1
9
9
2
1
9
9
2
1
9
9
3
1
9
9
3
1
9
9
3
1
9
9
4
1
9
9
4
1
9
9
5
1
9
9
5
1
9
9
5
1
9
9
6
1
9
9
6
1
9
9
7
1
9
9
7
1
9
9
8
1
9
9
8
1
9
9
8
1
9
9
9
1
9
9
9
2
0
0
0
2
0
0
0
2
0
0
0
2
0
0
1
2
0
0
1
2
0
0
2
2
0
0
2
2
0
0
3
2
0
0
3
2
0
0
3
2
0
0
4
2
0
0
4
2
0
0
5
2
0
0
5
2
0
0
5
Relacin entre la serie estandarizada de la Temperatura Promedio Mensual y el ONI
Temp E
ONI

Ingeniero Civil

52

Figura 32 Estacin Puerto Bolvar, comparacin de la serie estandarizada de la temperatura media
mensual con el ndice Ocenico del Nio (ONI)

Figura 33 Estacin Puerto Bolvar, coeficientes de correlacin para n meses de rezago de la serie
estandarizada de la temperatura media mensual con respecto al ndice Ocenico del Nio (ONI)

Analizando la velocidad del viento con relacin a la temperatura, se observ un
comportamiento de rezago similar, dado a que las dos variables se encuentran
estrechamente ligadas. Se realiz una correlacin entre la series mensual de velocidad y
temperatura encontrando un valor de -0,51 (confiabilidad del 95%) rezagada la velocidad
4 meses con respecto a la temperatura. Ver Figura 34.

Figura 34 Estacin Puerto Bolvar, comparacin de la serie de la temperatura con la velocidad
media mensual

De igual forma como se presenta en la Tabla 5, las correlaciones trimestrales marcaron
nuevamente al fenmeno El Nio asociado a la temperatura en su mximo DEF - DEF
(ONI - velocidad) y antecedindose un trimestre (MAM- DEF) al igual que SON-SON y
antecedindose el ndice JJA SON.
Tabla 5 Matriz de correlaciones trimestrales entre el ndice ONI y la Temperatura en C

Con relacin al comportamiento de la velocidad del viento en la zona de La Guajira, se
puede concluir que existe un ciclo diurno y anual promedio definido, el cual es de utilidad
a la hora de realizar clculos de produccin energtica de Parques Elicos en la zona, en
especial para el ciclo anual se identifican los perodos secos en los cuales la velocidad del
0
0,1
0,2
0,3
0,4
0,5
0,6
R
e
z
-
4
R
e
z
-
3
R
e
z
-
2
R
e
z
-
1
R
e
z
0
R
e
z
1
R
e
z
2
R
e
z
3
R
e
z
4
R
e
z
5
R
e
z
6
C
o
e
f
i
c
i
e
n
t
e

d
e

C
o
r
r
e
l
a
c
i
n
(
9
5
%

c
o
n
f
i
a
b
i
l
i
d
a
d
)
Correlacin entre Temp Estandarizada - ONI
23
24
25
26
27
28
29
30
31
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
8,0
9,0
10,0
1
9
8
8
1
9
8
8
1
9
8
9
1
9
8
9
1
9
9
0
1
9
9
0
1
9
9
1
1
9
9
1
1
9
9
2
1
9
9
2
1
9
9
3
1
9
9
3
1
9
9
4
1
9
9
4
1
9
9
5
1
9
9
5
1
9
9
6
1
9
9
6
1
9
9
7
1
9
9
7
1
9
9
8
1
9
9
8
1
9
9
9
1
9
9
9
2
0
0
0
2
0
0
0
2
0
0
1
2
0
0
1
2
0
0
2
2
0
0
2
2
0
0
3
2
0
0
3
2
0
0
4
2
0
0
4
2
0
0
5
2
0
0
5
T
e
m
p
e
r
a
t
u
r
a

(
C
)
V
e
l
o
c
i
d
a
d

d
e
l
V
i
e
n
t
o
(
m
/
s
)
Estacin Puerto Bolvar - Relacin entre la Velocidad del Viento y la Temperatura
Velocidad (m/s)
Temperatura (C)
DEF MAM JJA SON
DEF 0,57 0,24 0,28 0,14
MAM 0,49 0,15 0,27 0,22
JJA 0,12 -0,11 0,06 0,30
SON 0,05 -0,21 -0,01 0,33
ONI
Temperatura (C)

Ingeniero Civil

53

viento aumenta, los perodos en los cuales se podra presentar la complementariedad con
el recurso hdrico y el impacto positivo de El Nio en la velocidad del viento, en los meses
que se presentan los mnimos hidrolgicos. Sin embargo las mayores correlaciones no se
alcanzan en el trimestre DEF, tiempo en el cual El Nio a partir de su mximo desarrollo
ha disminuido los embalses histricamente, en algunos casos a nivel crtico.
5.1.3 Ciclo anual del ciclo diurno de la velocidad del viento

A partir de los datos horarios de la estacin Puerto Bolvar se construyeron las figuras que
se presentan a continuacin. En la Figura 35 se presenta la combinacin de los valores
ms altos al finalizar la tarde y dicha magnitud intensificada en los meses de mxima
velocidad (junio - agosto). Similares esquemas se realizaron para el Evento ENSO El Nio
y La Nia.

Figura 35 Ciclo anual del ciclo diurno de la velocidad del viento (m/s). Promedio en aos normales
en la estacin Puerto Bolvar.

Los valores de la velocidad del viento para los aos influenciados por el ENSO - El Nio a
las 16:00 horas del mes de julio son mayores que para la misma hora del mismo mes de
los aos La Nia. Sin embargo, nuevamente se refleja que para todos los meses del ao
durante la ocurrencia de El Nio la velocidad ser mayor. Ver Figura 36 y Figura 37.

Para los aos La Nia, en promedio, los mximos diarios se presentan en horas de la
noche entre las 21:00 y 23:00. Tambin se observan las menores velocidades del viento
registradas cuando ms impacta el evento (en la segunda mitad del ao) y en horas de la
maana, como se observa en la Figura 37.
La importancia de conocer la variacin del ciclo diurno de la velocidad del viento durante
las dos fases del ENSO, radica en los cambios de disponibilidad energtica de un parque
elico en el sitio, si bien es claro, en condiciones normales los mximos se presentarn a
mitad de la tarde aproximadamente a las 15:00 horas, durante La Nia, las mayores

Jun
Jul
Ago
Sep
Oct
Nov
Dic
Ene
Feb
Mar
Abr
May
Estacin Puerto Bollvar
Ciclio Anual del Cliclo Diurno de la Velocidad del Viento (m/s)
Aos Normales
10,0-12,0
8,0-10,0
6,0-8,0
4,0-6,0
2,0-4,0
0,0-2,0
V
e
l
o
c
i
d
a
d

d
e
l

v
i
e
n
t
o

(
m
/
s
)


Ingeniero Civil

54

velocidades se presentarn en horas de la noche cuando aumenta la demanda horaria de
energa elctrica en el pas. Este anlisis deber ampliarse con la informacin de otras
estaciones de la zona, y de igual manera los estudios de complementariedad energtica a
nivel horario, que permitan simular escenarios de oferta energtica al SIN.
A pesar de ser importante el aspecto anterior, el anlisis energtico de
complementariedad se realizar a nivel mensual considerando los ciclos promedios
anuales de series de velocidad del viento y series de caudales.

Figura 36 Estacin Puerto Bolvar, ciclo anual del ciclo diurno de la velocidad del viento (m/s).
Promedio en aos El Nio.

Figura 37 Estacin Puerto Bolvar, ciclo anual del ciclo diurno de la velocidad del viento (m/s).
Promedio en aos La Nia.

Jun
Jul
Ago
Sep
Oct
Nov
Dic
Ene
Feb
Mar
Abr
May
Aos Nio
10,0-12,0
8,0-10,0
6,0-8,0
4,0-6,0
2,0-4,0
0,0-2,0
Jun
Jul
Ago
Sep
Oct
Nov
Dic
Ene
Feb
Mar
Abr
May
Aos Nia
10,0-12,0
8,0-10,0
6,0-8,0
4,0-6,0
2,0-4,0
0,0-2,0
V
e
l
o
c
i
d
a
d

d
e
l

v
i
e
n
t
o

(
m
/
s
)

V
e
l
o
c
i
d
a
d

d
e
l

v
i
e
n
t
o

(
m
/
s
)


Ingeniero Civil

55


Ingeniero Civil

56

6 COMPLEMENTARIEDAD ENTRE LA HIDROLOGA DE VARIAS REGIONES DE
COLOMBIA CON LA VELOCIDAD DEL VIENTO DE LA GUAJIRA.

En este captulo se analiza el comportamiento de la hidrologa de varias regiones del pas
con relacin a la velocidad del viento en la zona de La Guajira colombiana. El objetivo
particular en este caso es identificar, con informacin de estaciones hidromtricas
localizadas en diferentes zonas del pas, las que complementan mejor su hidrologa con la
velocidad del viento, e independientemente que se puedan construir proyectos
hidroelctricos en dichos ros. Lo anterior teniendo en cuenta que las centrales
hidroelctricas y los futuros proyectos de la empresa ISAGEN S.A. E.S.P. tambin estn
localizados en diferentes regiones del pas.

Como es de conocimiento general ISAGEN ha instalado torres de medicin de la
velocidad del viento y dems variables meteorolgicas para el desarrollo de proyectos de
instalacin de Parques Elicos en varios sitios de La Guajira.

Dada la variabilidad de la hidrologa de Colombia, se realiz un anlisis con informacin
de caudales medios mensuales de varias regiones del pas contrastndola con la serie de
velocidad del viento de la estacin Puerto Bolvar localizada en La Guajira colombiana. La
informacin utilizada comprende las series de caudales promedios mensuales de los ros
Don Diego, Sin, Peralonso, Catatumbo y Pamplonita localizados en la zona norte y
nororiente del pas, los ros Sucio, Sumapaz, Atrato y Luisa localizados en la zona central
del pas, y los ros Pata, Guaitar, Cauca, La Plata y San Pedro, localizados en el sur del
territorio colombiano. En la Figura 38 se presenta la localizacin de las estaciones sobre
los ros mencionados, y en la Tabla 6 se presenta la informacin detallada de dichas
estaciones.
Tabla 6 Informacin de estaciones de caudales utilizadas
Estacin Corriente Localizacin
Municipio -
Departamento
Perodo de
informacin
Cant. Variable
Promedi
o (m
3
s
-1
)
Resolucin
Puente
Carretera
Don Diego
1115'N
7342' W
Santa Marta
Magdalena
1965 - 2004 40 Caudal 39,24 mensual
Villalosada La Plata
219 N'
7557' W
La Plata
Huila
1964 - 2003 40 Caudal 40,51 mensual
Agroyaco Guaitara
134' N
7727' W
Los Andes
Nario
1965 - 2002 38 Caudal 79,05 mensual
Bellavista Atrato
634' N
7634' W
Bojay
Choc
1965 - 2002 38 Caudal 2420,8 mensual
Quince
Letras
Catatumbo
830'N
7315' W
Teorama
N. de Santander
1964 - 2003 40 Caudal 181,31 mensual
Pavo Real Luisa
413' N
7512' W
Valle de San Juan
Tolima
1959 - 2003 45 Caudal 5,06 mensual
La Don
Juana
Pamplonita
742' N
7236' W
Chinacota
Santander
1965 - 2003 39 Caudal 7,08 mensual
Paletara Negro
210' N
7628'W
Purac
Cauca
1965 - 2005 41 Caudal 2,00 mensual
Puente
Guascas
Pata
136' N
7726' W
Policarpa
Nario
1965 - 2002 38 Caudal 225,89 mensual
Cornejo Peralonso
755' N
7239' W
Cayetano
Santander
1969 - 2003 35 Caudal 16,36 mensual
El Anil Sucio
652' N
7608' W
Uramita
Antioquia
1965 - 2002 38 Caudal 35,89 mensual
Itarca San Pedro
130' N
7528' W
La Montaita
Caquet
1964 - 2001 38 Caudal 73,25 mensual
Montera Sin
846' N
7552' W
Montera
Crdoba
1963 - 2005 43 Caudal 376,56 mensual
Lomitas Cauca
207' N
7627' W
Purac
Cauca
1970 - 2004 34 Caudal 4,33 mensual

Ingeniero Civil

57

Figura 38 Localizacin de las estaciones sobre los ros seleccionados para estudio y estacin
Puerto Bolvar referencia para el anlisis de complementariedad.
0
20
40
60
80
100
120
140
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro San Pedro - Ciclo Anual
0
1
2
3
4
5
6
7
8
9
10
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Luisa - Ciclo Anual
0
10
20
30
40
50
60
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Sucio - Ciclo Anual
0
50
100
150
200
250
300
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Catatumbo - Ciclo Anual
0
2
4
6
8
10
12
14
16
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Pamplonita - Ciclo Anual
0
5
10
15
20
25
30
35
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Peralonso - Ciclo Anual
0
10
20
30
40
50
60
70
80
90
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Don Diego - Ciclo Anual
0
100
200
300
400
500
600
700
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Sin - Ciclo Anual
0
500
1000
1500
2000
2500
3000
3500
4000
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Atrato - Ciclo Anual
0
2
4
6
8
10
12
14
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Cauca - Ciclo Anual
0
50
100
150
200
250
300
350
400
450
500
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro Pata - Ciclo Anual
0
10
20
30
40
50
60
70
80
C
a
u
d
a
l
(
m
3
/
s
)
Mes
Ro La Plata - Ciclo Anual
Convenciones:
Condicin Normal
ENSO - El Nio
ENSO - La Nia

Ingeniero Civil

58

Con el objeto de hacer comparables los caudales de las diferentes estaciones, con los
datos de velocidad del viento, fue necesario estandarizar las series a partir de sus valores
promedios y desviaciones estndar. El anlisis regional realizado grfico y con
coeficientes de correlacin inversos calculados con las series estandarizadas, present
los siguientes comportamientos:

Para el caso de los ros ubicados al norte del pas, y especficamente los localizados en el
noroccidente colombiano en la zona de los Santanderes, es decir los ros Pamplonita y
Peralonso, se presenta una complementariedad asociada a los coeficientes de correlacin
inversos ms altos. Para los ros Catatumbo y Sin, la complementariedad solo se
presenta en algunos meses del ao, entre septiembre y noviembre; en los meses de junio
a agosto, se presentan condiciones similares, es decir altas velocidades del viento y
caudales altos. Sin embargo, cuando se presentan vientos mnimos, es decir entre
septiembre y noviembre, los valores de caudales de estos ltimos son menores que los
ubicados sobre los Santanderes. Ver Figura 39.

Figura 39 Anlisis de complementariedad regional entre ros del norte y noroeste colombiano y la
estacin Puerto Bolvar.

Para el caso de los ros del centro del pas, los ros Sucio, Atrato y Luisa (localizados en
las cordilleras central y occidental), en los meses de septiembre, octubre y noviembre
cuando se presenta las menores velocidades del viento durante el ao, se presentan los
caudales ms altos a nivel mensual, y por ende las mayores correlaciones inversas. Ver
Figura 40.

El ro Sumapz localizado en el centro del pas en la cordillera oriental, no presenta
comportamientos de complementariedad.

-2,5
-2,0
-1,5
-1,0
-0,5
0,0
0,5
1,0
1,5
2,0
2,5
Pto Bolvar Sin Catatumbo Peralonso Pamplonita
Correlacin Pto. Bolvar
Catatumbo 0,42
Sin -0,19
Peralonso -0,70
Pamplonita -0,56

Ingeniero Civil

59

Figura 40 Anlisis de complementariedad regional entre ros del centro de Colombia y la estacin
Puerto Bolvar.
Para el caso de los ros en el sur de Colombia, se observ la menor complementariedad
regional, Los ros Cauca, San Pedro y La Plata son los de menor correlacin con la serie
de vientos de la estacin Puerto Bolvar y por ende los de menor complementariedad,
presentan altos caudales al igual que las velocidades del viento. Sin embargo, se resalta
que los ros Guatara y Pata se comportan intranualmente de manera similar que la
velocidad el viento, mejor an, si se analiza con un rezago la serie de caudales un (1) mes
la correlacin inversa podra ser mucho ms alta. Ver Figura 41.

Figura 41 Anlisis de complementariedad regional entre ros del sur de Colombia y la estacin
Puerto Bolvar.
-2,5
-2,0
-1,5
-1,0
-0,5
0,0
0,5
1,0
1,5
2,0
2,5
Pto Bolvar Atrato Sumapaz Sucio Luisa
-2,5
-2,0
-1,5
-1,0
-0,5
0,0
0,5
1,0
1,5
2,0
2,5
Pto Bolvar Cauca Pata La Plata San Pedro Guaitara
Cauca 0,54
San Pedro 0,33
Pata -0,14
Guaitara 0,37
La Plata 0,40
Atrato -0,48
Sumapaz 0,03
Sucio -0,65
Luisa -0,34

Ingeniero Civil

60

El anlisis regional anterior refleja la mayor complementariedad entre los recursos para
los ros de la zona nororiente del pas, sin embargo, en el centro del pas algunas
centrales hidroelctricas existentes o futuros proyectos, podran complementarse con
parques elicos en La Guajira.

Teniendo en cuenta los resultados encontrados en dichas zonas y analizando la
localizacin de las Centrales de Generacin en operacin y los Proyectos Hidroelctricos
de ISAGEN, se realizaron los anlisis de complementariedad que se presentan en el
captulo 7.


Ingeniero Civil

61

7 SIMULACIN Y ANLISIS DE GENERACIN DE ENERGA DE PARQUES
ELICOS GENRICOS EN LA GUAJIRA

En este captulo se presenta el desarrollo de la modelacin de distintos parques elicos
de capacidades instaladas asociadas a las centrales hidroelctricas de ISAGEN,
analizando la energa generada por las centrales existentes y por dichos parques, a partir
de la informacin de vientos de la estacin Puerto Bolvar y la ayuda de la estacin virtual
del reanlisis, e informacin de generacin de energa disponible en el sitio web de XM. El
anlisis permitir determinar la existencia de correlacin en la generacin elica e hdrica
de ISAGEN, el cual se presenta en el captulo 8.

La simulacin de la produccin energtica de un parque elico genrico comprende
muchos elementos como el tipo de aerogenerador o turbina, la clase de vientos segn la
Norma IEC 61400 I (INTERNATIONAL STANDARD IEC 61400-1, Third edition, Wind
turbines Part 1: Design requirements, 2005), la distribucin de los aerogeneradores en
la zona del parque y sobre todo la serie de velocidades del viento a resolucin horaria.

Para la simulacin energtica se seleccion la turbina Nordex 1300, la seleccin de este
aerogenerador se bas en tener un referente comparativo en la zona, dado a que los
aerogeneradores del Parque Elico Jeprachi tienen este tipo de turbinas. Adicionalmente
durante la seleccin de dicha mquina EEPPM consider otras mquinas existentes a la
fecha de construccin y las consideraciones de la ubicacin del parque, terreno desrtico
y altas temperaturas, entre otras. (Revista EPM Vol. 15, 2004).

Cuando se cuenta con la informacin suficiente, para los clculos de energa se utilizan
software que permiten considerar los diferentes equipos, alturas, dimetros y curvas de
potencia de los aerogeneradores. Los softwares como el WindPro (desarrollado por EMD
International A/S) y WASP (The Wind Atlas Analysis and Application Program,
desarrollado por National Laboratory for Sustainable Energy de Dinamarca) son algunos
de los usados para las estimaciones energticas, sin embargo, requieren licencia e
informacin detallada de un parque elico como la distribucin de mquinas en el rea,
entre otras mencionadas anteriormente. Este anlisis utiliz informacin horaria para los
clculos energticos junto con la curva de mquina de la N60, incluyendo un factor del
20% por indisponibilidad por operacin y mantenimiento, prdidas por estelas generadas
entre aerogeneradores, prdidas asociadas por variacin en la temperatura en la zona y
prdidas elctricas.

Otra consideracin que no tiene en cuenta la estimacin energtica cuando no se realiza
con la ayuda de un software es la variacin de la velocidad del viento por la altura dentro
del rea cubierta por el aerogenerador.
7.1 Aerogeneradores Nordex N60/1.300
Estas mquinas diseadas por la compaa Nordex de Alemania tienen un dimetro de
rotor de 60 metros e igual altura de bujes, es decir, las palas son de 30 metros de largo.
En la Figura 42 se presenta su respectiva curva de potencia. (Revista EPM Vol 15, 2004).


Ingeniero Civil

62

Figura 42 Curva de Potencia Terica Aerogenerador Nordex N60/1.300. (Nordex Energy Catlogo N60/1300)

La potencia mxima que puede alcanzar la turbina es de 1.300 kW (1,3 MW), cual se
logra con una velocidad ptima de 15 ms
-1
, la velocidad mnima para generacin debe ser
de 4 ms
-1
y la velocidad mxima a la cual el equipo an genera la mxima potencia es del
orden de 25 ms
-1
.

Como se mencion anteriormente ests maquinas fueron las seleccionadas dentro del
estudio del Parque Jeprachi hace alrededor de 10 aos, desde entonces la industria de
los fabricantes de aerogeneradores ha evolucionado significativamente. Algunos
fabricantes como GENERAL ELECTRIC hoy en da ofrecen aerogeneradores de 2,5 MW
y hasta 3,0 MW para turbinas offshore (parques elicos en el mar) y la compaa VESTAS
con turbinas que alcanzan los 3,0 MW para parques en tierra. Esto indica cierto criterio
conservador para la estimacin energtica de los parques elicos supuestos en este
anlisis, no solo por las capacidades de las mquinas de hoy en da, sino tambin a las
alturas que se estn diseando (lo que aumenta la velocidad del viento) y las mejoras en
los sistemas de caja multiplicadora y en otros elementos del aerogenerador que aumentan
su eficiencia.

En la Figura 43 se observan la evolucin y tendencias respecto al tamao de mquinas y
sus respectivos tamaos de palas.
Curva Terica N60/1300
Velocidad del Viento Vs. Potencia
0
200
400
600
800
1000
1200
1400
1600
0 3 6 9 12 15 18 21 24 27
P
o
t
e
n
c
i
a

(
k
W
)

Ingeniero Civil

63

Figura 43 Comparacin de Tamao de algunos aerogeneradores tpicos. Crecimiento en el tiempo.
Adaptada de SRREN IPCC, 2007.
Para el clculo energtico se utiliz la serie de velocidades del viento de cada hora de la
estacin Parque Elico acotando la produccin entre 4 y 25 ms
-1
. Luego se utiliz el
promedio de las potencias generadas en cada intervalo de tiempo y se calcul la energa
media mensual agregada.

Sin embargo, como solo se tiene un ao de informacin a nivel horario de la estacin, las
series velocidades del viento se alargaron en el tiempo con la ayuda de la informacin
horaria de la estacin Puerto Bolvar, extendiendo as, la serie mensual reconstruida
desde 1986 al 2005. Luego, con la informacin energtica mensual se realiz un ajuste
para extender la serie de energa generada desde enero de 1960 a diciembre de 2005
con la ayuda de los datos mensuales de velocidad del viento Estacin Virtual del
Reanlisis, generando una regresin lineal entre la velocidad y la energa a nivel mensual,
como se observa en el ajuste de la Figura 44. La correlacin entre la serie Puerto Bolvar
y el Reanlisis es de 0,66, con el 95% de confiabilidad.

Cabe anotar que las velocidades del viento de las estaciones de Puerto Bolvar y
Reanlisis son a la altura de 10 m y que para realizar la reconstruccin de los datos
faltantes y extrapolaciones se utiliz un coeficiente de Hellmann a la altura de 60 m e igual
a 0,14; y el cul se calcul a partir del perfil vertical de velocidad del viento.

Por ltimo, la serie de energas generadas para un (1) aerogenerador se escala
dependiendo el nmero de mquinas a simular, es decir dependiendo de la capacidad
instalada de los parques supuestos para el anlisis.


Ingeniero Civil

64

Figura 44 Ajuste a nivel mensual de velocidades del Viento vs. Potencia en la estacin Puerto
Bolvar. Correlacin

En la Figura 45 y a manera de comparacin, se presenta la produccin energtica
mensual de un Parque Genrico como el propuesto en este anlisis de 19,5 MW
emplazado en la zona de Puerto Bolvar comparada con el Parque Elico Jeprachi el cual
dispone de datos de generacin media mensual obtenidos del portal Nen de XM
(www.xm.com.co). Para la comparacin con el Parque Jeprachi se tom el perodo
comn entre series comprendido entre el 2004 y 2005.

La produccin real sigue igual comportamiento que el Parque simulado alcanzando
valores cercanos a nivel mensual. Se aclara que la informacin disponible en el sistema
Nen no incluye la disponibilidad de aerogeneradores que producen la energa mensual
presentada, por lo que se desconoce el nmero de mquinas asociadas a dicha
produccin.

Figura 45 Comparacin de la produccin energtica de Parque Elico Simulado vs. Jeprachi 2004
2005.
0
1
2
3
4
5
6
7
8
9
10
0 2 4 6 8 10 12 14
P
u
e
r
t
o

B
o
l
v
a
r
Reanlisis
Correlacin Reanlisis - Puerto Bolvar Velocidad del Viento Vs. Potencia
y = 110.55x - 372.19
R
2
= 0.9807
-
200
400
600
800
1,000
1,200
1,400
- 2 4 6 8 10 12 14 16
Velocidad (m/s)
P
o
t
e
n
c
i
a

(
k
W
)
0
1
2
3
4
5
6
7
8
9
10
m
a
r
-
0
4
a
b
r
-
0
4
m
a
y
-
0
4
j
u
n
-
0
4
j
u
l
-
0
4
a
g
o
-
0
4
s
e
p
-
0
4
o
c
t
-
0
4
n
o
v
-
0
4
d
i
c
-
0
4
e
n
e
-
0
5
f
e
b
-
0
5
m
a
r
-
0
5
a
b
r
-
0
5
m
a
y
-
0
5
j
u
n
-
0
5
j
u
l
-
0
5
a
g
o
-
0
5
s
e
p
-
0
5
o
c
t
-
0
5
n
o
v
-
0
5
d
i
c
-
0
5
G
W
h
-
m
e
s
Mes - ao
Comparacin de Produccion Energtica de Parques Elicos: 2004 -2005
Parque 20 MW
Jeprrachi

Ingeniero Civil

65

De igual forma en la Figura 46 se presenta el ciclo anual de produccin energtica de un
Parque de 19,5 MW emplazado en Puerto Bolvar teniendo en cuenta las consideraciones
planteadas anteriormente y comparndolas con el Parque Elico Jeprachi el cual posee
igual potencia y tipo de mquina, y con el objeto de comparar el comportamiento de
ambos Parques.

Figura 46 Comparacin de ciclo anual de produccin energtica promedio Parque simulado 20 MW
vs. Jeprachi.
A pesar de mantener similar comportamiento, la diferencia en magnitud de produccin
media mensual radica en que el factor de planta calculado es diferente al real; mientras
que el Parque genrico simulado obtiene una produccin energtica media anual de 81,9
GWh con un factor de planta de 0,46, Jeprachi solo alcanza en su promedio histrico
52,6 GWh con un factor de planta de 0,32.
8 ANLISIS DE GENERACIN COMPLEMENTARIA HIDRULICA - ELICA

Retomando lo expuesto en el captulo anterior, y a partir de la modelacin de una (1)
turbina elica se simularn parques elicos de diferentes capacidades instaladas
asociadas a los proyectos y centrales hidroelctricas de ISAGEN con el objetivo de
identificar la complementariedad de las dos fuentes y cuantificar las capacidades elicas
requeridas para generar igual energa media de las hidroelctricas. Para lo anterior fue
necesario realizar un anlisis previo de correlacin de series hidrolgicas y de velocidad
del viento y revisar el comportamiento hidrolgico ante Eventos ENSO para identificar los
comportamientos de dichas variables ante El Nio y La Nia.

La informacin de caudales utilizada en este anlisis comprende las series de caudales
medios mensuales totales afluentes a las Centrales Hidroelctricas y Proyectos de
Generacin de ISAGEN de los ros Sogamoso (484,3 m
3
s
-1
), Amoy (18,2 m
3
s
-1
), La Miel
(83,5 m
3
s
-1
), Calderas (8,5 m
3
s
-1
), Jaguas (37,2 m
3
s
-1
) y los caudales artificiales afluentes
a la central hidroelctrica San Carlos (136,8 m
3
s
-1
). Se anota que los ros Sogamoso y
Amoy hacen parte de los proyectos hidroelctricos que obedecen a los mismos nombres.
0
1
2
3
4
5
6
7
8
9
10
G
W
h

-
m
e
s
Mes
Comparacin de Produccion Energtica de Parques Elicos
Jeprrachi
Parque 20 MW

Ingeniero Civil

66

La Miel, Calderas, San Carlos y Jaguas son centrales hidroelctricas en operacin de
ISAGEN. En la Figura 47 se observa la localizacin de los sitios y estaciones analizadas
en este estudio.

Figura 47 Localizacin de Sitios de Inters: Centrales, Proyectos y Estaciones de Velocidad del
Viento.

8.1 CICLO ANUAL DE LAS SERIES DE CAUDALES MENSUALES DE LOS SITIOS
DE INTERS

La Figura 48 presenta el ciclo anual de las series de caudales medios mensuales
afluentes a las Centrales Hidroelctricas y Proyectos de ISAGEN. Adicionalmente se
presenta el ciclo anual de velocidad del viento en las estaciones Puerto Bolvar y
Reanlisis NCEP/NCAR.

Para las series de caudales se observa el comportamiento bimodal del ciclo anual
promedio.

Proyecto Hidroelctrico
Sogamoso
Amoy
Estacin Reanlisis
Central
Hidroelctrica Miel I
Central Hidroelctrica San
Carlos
Central Hidroelctrica Jaguas
Central Hidroelctrica Calderas

Ingeniero Civil

67

Figura 48 Ciclo Anual de caudales y vientos medios mensuales para los sitios de inters

De la figura anterior se observa que tanto las centrales hidroelctricas existentes como los
proyectos, en los meses de julio a septiembre se presentan mnimos en la hidrologa,
meses en los cuales la velocidad del viento es mayor. Similarmente para los meses entre
mediados de diciembre y marzo.
8.2 ESCALA INTERANUAL: RELACIN DE LOS CAUDALES MEDIOS MENSUALES
CON EL FENMENO EL NIO-OSCILACIN DEL SUR (ENSO)

Para las series de caudales medios mensuales se identific el ciclo anual y su correlacin
con los eventos ENSO. Las correlaciones se realizaron con la serie completa de cada sitio
sin estandarizar para evitar la remocin del ciclo anual y con la serie de valores del
Oceanic Nio Index - ONI.

Adicionalmente se realizaron anlisis de correlacin trimestral para cada serie de
caudales con el ONI. Para todas las series de caudales se estimaron las correlaciones,
con resultados que indican altos valores con relacin al ndice ONI. La serie que mayor
correlacin present con el ENSO corresponde al ro La Miel. A continuacin, a manera
de ejemplo, se presentan algunos de los resultados encontrados para La Miel, Sogamoso
y Amoy.
8.2.1 Anlisis serie de caudales del Proyecto Hidroelctrico Sogamoso

Coeficiente Correlacin con el ndice ONI: - 0,24 Rezago 3 meses de disminucin de
caudales en relacin con el ndice ONI. Las correlaciones trimestrales se observan en la
Sogamoso 1960 - 2003
Amoy 1974 - 2006
Estacin Puerto
Bolvar 1986-2007
Estacin Reanlisis
Central Hidroelctrica Miel I
1963 - 2008
Central Hidroelctrica San Carlos
1984 - 2008
Central Hidroelctrica Jaguas
1955-2008
Central Hidroelctrica Calderas
1991- 2007
0
100
200
300
400
500
600
700
800
900
C
a
u
d
a
l
(
m
3
/
s
)
Ciclo Anual Promedio - Caudales Sogamoso
0
5
10
15
20
25
30
C
a
u
d
a
l
(
m
3
/
s
)
Ciclo Anual Promedio - Caudales Amoy
0
2
4
6
8
10
12
14
16
C
a
u
d
a
l (
m
3
/ s
)
Ciclo Anual Promedio - Caudales Calderas
0
20
40
60
80
100
120
140
C
a
u
d
a
l (
m
3
/ s
)
Ciclo Anual Promedio - Caudales Miel I
0
10
20
30
40
50
60
C
a
u
d
a
l (
m
3
/ s
)
Ciclo Anual Promedio - Caudales Jaguas
0
10
20
30
40
50
60
C
a
u
d
a
l ( m
3
/ s
)
Ciclo Anual Promedio - Caudales San Carlos

Ingeniero Civil

68

Tabla 7; es importante destacar que el trimestre DEF es el ms impactado por el evento
ENSO.

Tabla 7 Matriz de correlaciones trimestrales caudales medios mensuales Sogamoso vs. ONI

8.2.2 Anlisis serie de caudales del Proyecto Hidroelctrico del ro Amoy

Coeficiente Correlacin con el ndice ONI: - 0,37 Rezago 3 meses de disminucin de
caudales en relacin con el ndice ONI. Las correlaciones trimestrales se observan en la
Tabla 8; es importante destacar que al igual que la serie Sogamoso el trimestre DEF es el
ms impactado por el Evento ENSO.

Tabla 8 Matriz de correlaciones trimestrales caudales medios mensuales Amoy vs. ONI

8.2.3 Anlisis serie de caudales de la Central Hidroelctrica Miel I

Coeficiente Correlacin con el ndice ONI: - 0,39, sin rezago. Las correlaciones
trimestrales se observan en la Tabla 9; es importante destacar que al igual que la serie
Sogamoso y la serie Amoy, el trimestre DEF es el ms impactado por el Evento ENSO.
Adicional a este, en el trimestre JJA tambin se presentan correlaciones que se
consideran altas.

Tabla 9 Matriz de Correlaciones Trimestrales Caudales medios mensuales Miel I vs. ONI

Se resalta que a pesar que la serie natural afluente a la Central Hidroelctrica San Carlos
se correlaciona altamente con el ndice ONI, a esta central le ingresan otros aportes
modificados por las polticas de operacin del embalse de la Central Hidroelctrica
Guatap, por lo cual en el anlisis energtico no se utiliz esta central (la de mayor
capacidad instalada del pas) para productos comparativos con un parque elico genrico.

DEF MAM JJA SON
DEF -0,44 -0,26 -0,02 -0,02
MAM -0,32 -0,15 -0,22 -0,21
JJA -0,10 0,06 -0,26 -0,31
SON -0,02 0,17 -0,19 -0,24
Caudales Medios Mensuales Sogamoso (m
3
/s)
ONI
DEF MAM JJA SON
DEF -0,59 -0,39 -0,27 -0,31
MAM -0,49 -0,31 -0,29 -0,42
JJA -0,21 -0,16 -0,22 -0,38
SON -0,10 -0,06 -0,15 -0,35
ONI
Caudales Medios Mensuales Amoy (m
3
/s)
DEF MAM JJA SON
DEF -0,50 -0,34 0,03 -0,09
MAM -0,41 -0,38 -0,26 -0,27
JJA -0,15 -0,25 -0,44 -0,41
SON -0,06 -0,15 -0,44 -0,40
ONI
Caudales Medios Mensuales Miel I (m
3
/s)

Ingeniero Civil

69

8.3 ANLISIS DE COMPLEMENTARIEDAD DE LAS SERIES HIDROLGICAS CON
LA SERIE DE VIENTOS LA GUAJIRA

Luego de realizar el anlisis por separado de las series de caudales en los sitios de
inters y las series de viento en La Guajira, se procedi a comparar cada serie de
caudales con la serie de velocidad mensual del viento de la estacin Puerto Bolvar y la
estacin Virtual Reanlisis. Para hacer la comparacin, se estandarizaron las series
mencionadas. En la Figura 49 se observa el comportamiento de las series analizadas. En
general es clara la complementariedad de los caudales en el interior del pas y los vientos
en la Guajira colombiana.
La serie mensual del ro Amoy presenta un comportamiento diferente a las dems series
de caudales, presentando valores de caudales altos cuando se presenta el verano en las
dems series de caudales y especficamente en los meses de junio y julio. Esto tambin
se ve reflejado en el coeficiente de correlacin igual a cero (0) de dicha serie con la serie
de vientos de la estacin de Puerto Bolvar.

Figura 49 Complementariedad caudales y vientos, series mensuales multianuales estandarizadas

Las series de caudales medios mensuales presentaron los siguientes coeficientes de
correlacin con las de velocidad del viento. Ver Tabla 10.
Tabla 10 Matriz de correlaciones la velocidad del viento caudales: estacin Reanlisis, estacin
Puerto Bolvar vs. Caudales sitios de inters
Central - Proyecto Hidroelctrica (o)
Puerto Bolvar
vs. Caudales
Reanlisis vs.
Caudales
Central Hidroelctrica San Carlos - 1240 MW (1987 2005) -0,32 -0,52
Central Hidroelctrica Jaguas - 170 MW (1987 2005) -0,32 -0,47
Central Hidroelctrica Miel I - 396 MW (1987 2005) -0,43 -0,28
Central Hidroelctrica Calderas - 19.5 MW (1991 2001) -0,24 -0,41
Proyecto Hidroelctrico Amoy - 80 MW (1987 2005) 0,00 -0,12
Proyecto Hidroelctrico Sogamoso - 820 MW (1987 2003) -0,38 -0,55
-3
-2
-1
0
1
2
3
Calderas Amoy Miel I Sogamoso Jaguas San Carlos Pto Bolvar

Ingeniero Civil

70

Es de resaltar que la serie del Reanlisis presenta mayores correlaciones con los
caudales, a excepcin de la serie de Miel I. No existen datos disponibles para calificar
este comportamiento; sin embargo, se puede observar en la Tabla 10 que la serie del
Reanlisis disminuye su correlacin con los caudales mensuales de norte a sur, en
relacin con los sitios de inters; es decir, Sogamoso que es el sitio ms al norte tiene la
mayor correlacin -0,55 y Amoy -0,12, la que se encuentra ubicada ms al sur de los
sitios de inters. Estos resultados son coherentes con el anlisis presentado para las
diferentes regiones del pas.
8.4 ANLISIS DE GENERACIN DE ENERGA HISTRICA CENTRALES ISAGEN
VS. PARQUES ELICOS GENRICOS.

Considerando la informacin disponible en la base de datos de Nen de XM, agencia
administradora del mercado de energa en Colombia, se revis la complementariedad de
la generacin histrica de las plantas hidroelctricas de ISAGEN comparndolas con
parques elicos genricos de igual capacidad instalada.

A pesar de existir correlacin simultnea inversa (sin rezago) entre la serie de caudales
afluentes a la Central Miel I y la serie de vientos, las energas medias mensuales
expresadas en kWh no reportan complementariedad, esto sucede porque la poltica
operacin del embalse est asociada al comportamiento del sistema, el cual se opera
para satisfacer una demanda con una oferta hdrica, que a la fecha no ofrece
complementariedad con otro recurso y por lo tanto la operacin comercial de ISAGEN .
Ver Figura 50.

Igual comportamiento se observa en el ciclo anual de generacin media mensual del
Parque Elico genrico y la central Miel I, ver Figura 51.

Figura 50 Comparacin entre la generacin mensual entre Miel I y un Parque Elico de 396 MW

0
50
100
150
200
250
300
G
e
n
e
r
a
c
i
n

G
W
h
-
m
e
s
Mes - ao
Miel I
Parque Elico 396 MW

Ingeniero Civil

71

Figura 51 Comparacin entre el ciclo anual de generacin mensual entre Miel I y un Parque Elico
de 396 MW

Para el caso de las centrales San Carlos y Jaguas, comparndolas con parques de
similares capacidades instaladas se observa complementariedad en los primeros meses
del ao y en la mitad del ao con un mximo del aporte elico en los meses de julio y
agosto. Es importante aclarar que la operacin de las plantas obedece a criterios de
comercializacin de la generacin o despacho de las plantas por parte del agente,
considerando las condiciones del mercado de energa, criterios y condiciones no
contempladas en este anlisis.

El comportamiento de las series de generacin media mensual se presenta en las Figura
52 y Figura 53.

Figura 52 Comparacin entre el ciclo anual de generacin mensual entre Jaguas y un Parque
Elico de 170 MW

0
200
400
600
800
1000
1200
1400
1600
1800
0
20
40
60
80
100
120
140
160
180
200
E
n
e
r
g
a

M
e
d
i
a

A
n
u
a
l

(
G
W
h
)
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
)
Parque 396 MW
Miel I
E. Media Miel I
E. Media Parque
0
100
200
300
400
500
600
700
800
0
10
20
30
40
50
60
70
80
90
E
n
e
r
g
a

M
e
d
i
a

A
n
u
a
l

(
G
W
h
)
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
)
Parque 170 MW
Jaguas
E. Media Jaguas
E. Media Parque

Ingeniero Civil

72

Figura 53 Comparacin entre el ciclo anual de generacin mensual entre San Carlos y un Parque
Elico de 1240 MW

En las Figura 51, Figura 52 y Figura 53 se observa que la energa media anual de los
parques elicos genricos es un poco mayor que la de las centrales hidroelctricas; esto
se debe a que los factores de planta (o capacidad) histricos reflejan las disponibilidad del
recurso, la capacidad mxima de la planta sin incluir disponibilidad y criterios de operacin
comercial. Para el caso de la central Miel I corresponde a 0,4, Jaguas a 0,47 y San Carlos
a 0,52; estos ndices se calcularon con la informacin disponible en el portal del
administrador del Mercado de Energa, XM.

No se realiza un anlisis similar para la central Calderas dado a sus perodos
interrumpidos de generacin debido a eventos externos a su operacin comercial.
8.5 ANLISIS DE GENERACIN COMPLEMENTARIA PROYECTOS
HIDROELCTRICOS VS. PARQUES ELICOS

En este numeral se analiza la produccin energtica de un proyecto a filo de agua y uno
con embalse de ISAGEN con relacin a parques elicos que igualen su energa media
simulados de manera aislada y con el objeto de comparar los ciclos anuales de las
generaciones en condiciones normales y durante las dos fases del ENSO, y la posibilidad
de complementariedad entre ambas.

Para las simulaciones de energa hidrulica de proyectos hidroelctricos se elabor un
modelo de simulacin de energa basado en la Regla Normal de Operacin a nivel
mensual, el cual considera la operacin aislada de los proyectos del Sistema
Interconectado Nacional. Cabe anotar que la operacin futura de las centrales est ligada
a muchos otros elementos como son las polticas internas de operacin de ISAGEN, los
contratos de venta de energa a largo plazo, la variabilidad de los precios de energa y las
compras y ventas de energa en el mercado spot. A este nivel, el anlisis no considera
estos elementos mencionados.
0
1000
2000
3000
4000
5000
6000
0
100
200
300
400
500
600
E
n
e
r
g
a

M
e
d
i
a

A
n
u
a
l

(
G
W
h
)
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
)
Parque 1240 MW
San Carlos
E. Media San Carlos
E. Media Parque

Ingeniero Civil

73

8.5.1 Caractersticas del Modelo de Simulacin de Generacin de Energa

A continuacin se describen los criterios, supuestos y la metodologa del modelo de
simulacin energtica utilizado; se definen los parmetros iniciales o de entrada para el
modelo (para cada proyecto), el criterio o regla de operacin del embalse, las
restricciones asociadas y los resultados o datos de salida del modelo.

Las caractersticas generales del modelo de simulacin energtica utilizado son las
siguientes:

- Tipo de Operacin: Se consider la operacin aislada y bajo la utilizacin de la
Regla Normal de Operacin (RNO). Se aplica para proyectos con embalse regulador
y para el caso de Amoy; dada su condicin de Filo de Agua.
- Resolucin Temporal y series hidrolgicas: La resolucin temporal del modelo es
mensual dado a que se quiere comprobar la complementariedad observada entre los
ciclos anuales de los recursos elico e hdrico. El mismo intervalo temporal de la
variacin y regulacin del embalse tambin se consider mensual.
- Regla de Operacin: La operacin del embalse se basa en la Regla Normal de
Operacin (RNO), en funcin de los caudales disponibles por almacenamiento o por
aportes del ro, la cual permite maximizar la produccin de energa firme (confiable)
del proyecto y minimizar los vertimientos, a partir del cumplimiento de una meta con
un 95% de confiabilidad.
- Variables de la simulacin: El modelo representa las principales variables
dinmicas del proceso de generacin, tales como:
- Relacin Volumen del embalse cota de operacin del mismo. El modelo permite
obtener el volumen asociado (en Hm
3
) a cada nivel de operacin a travs de la curva
de calibracin del embalse (Cota Volumen). Para las evaluaciones realizadas se
consider una relacin lineal entre la cota de operacin, el volumen til y salto bruto
(y neto) asociado.
- Se consideran las prdidas en energa por friccin sobre el salto bruto
aprovechable en cada una de las alternativas. Las prdidas friccionales fueron
calculadas teniendo en cuenta las longitudes y caractersticas principales de las
conducciones: rea, permetro, tipo de material de revestimiento, etc. Las prdidas
locales se consideraron como igual al 10% de las prdidas por friccin longitudinales.
- Restricciones. Las principales restricciones que considera el modelo de simulacin
son:
Caudal Ecolgico: El modelo considera una disminucin en el caudal afluente por
el efecto del caudal ecolgico (se puede considerar una variacin mensual), sin
embargo para el anlisis energtico del Proyecto Hidroelctrico Sogamoso no se
consider el caudal ecolgico por ser tener descarga a pie de Presa.
Capacidad mxima de operacin: La mxima capacidad de operacin de la planta
est asociada a la operacin de la misma bajo los condicionantes de caudal de

Ingeniero Civil

74

diseo y mxima cota de operacin, es decir, cuando la potencia reportada es
igual a potencia instalada de la central.
Serie Hidrolgica de caudales promedios mensuales: la serie de datos hidrolgicos
de entrada al modelo corresponde a los caudales promedios mensuales afluentes
al embalse.
Eficiencias: El modelo considera eficiencias constantes para los equipos
electromecnicos as: Generador: 98%, Transformador: 99.5% y Turbinas: 92%.
Es importante considerar que en la operacin real de la planta la eficiencia de las
turbinas depende del nivel de operacin del embalse y que en este anlisis no se
consider este efecto sobre la generacin energtica.
8.5.2 Funcionamiento del Modelo

El modelo de simulacin parte de los parmetros bsicos anteriormente mencionados,
como son: la serie hidrolgica de caudales medios diarios, el caudal de diseo
seleccionado, las eficiencias de los equipos electromecnicos, la altura de la presa (salto
bruto), entre otros.
Con respecto al clculo de las prdidas de energa por friccin, se considera que la
geometra de la conduccin es de tipo circular y se adopta un n de Manning por tramo de
tipo de material de revestimiento, a partir de estos, se calcula el coeficientes total de
prdidas hidrulicas K el cual considera el total de las prdidas de energa (prdidas por
friccin y prdidas locales).
La expresin matemtica utilizada es: H
f
= KQ
2
; donde Q es el caudal de diseo
considerado en m
3
s
-1
y H
f
son las prdidas hidrulicas (locales y por friccin).
Teniendo en cuenta los parmetros bsicos el modelo calcula la variacin de los saltos
netos asociados a los niveles de operacin del embalse y su potencia (MW) y la energa
(GWh-mes) asociada para cada mes agregando las energas medias diarias.
Al finalizar la modelacin y realizando la agregacin mensual se obtiene la Energa Media
Anual en unidades de GWh - ao.
8.5.3 Regla de Operacin del Embalse

La regla utilizada es la Regla Normal de Operacin (RNO) para embalses, la cual se
utiliza para maximizar la produccin de energa firme y minimizar los vertimientos.
El funcionamiento de la regla consiste en atender para cada perodo simulado una meta
de caudal (m
3
s
-1
) preestablecida.
Las simulaciones realizadas deben cumplir con el 95% de confiabilidad para todo el
perodo considerado y la meta en otras palabras corresponder a la energa firme de la
simulacin.

Ingeniero Civil

75

Para cada perodo simulado, se parte de un volumen o caudal necesario para generar la
meta y se compara con las disponibilidades de agua definidas en la RNO del embalse. La
disponibilidad corresponde al volumen almacenado de agua ms el caudal afluente en un
perodo establecido de tiempo, menos el caudal ecolgico considerado. S la
disponibilidad es menor que la meta, se presentar una falla y queda registrado para su
posterior conteo y para el clculo de confiabilidad. Si la disponibilidad es mayor o igual
que la meta, se cumple la meta y se genera esa cantidad de agua y almacenando los
excedentes sin sobrepasar la capacidad establecida del embalse.
Si se presenta una disponibilidad mayor que la meta y el embalse est a la mxima
capacidad de almacenamiento, se genera cumpliendo la meta, pero como el embalse est
en su nivel mximo, los excedentes permiten generar energa secundaria manteniendo a
capacidad total el embalse al final del perodo y minimizando los vertimientos. En este
caso la generacin estar restringida por el caudal de diseo de la planta. Si la
disponibilidad del recurso excede el volumen til del embalse y la capacidad mxima de
turbinacin, se genera la meta de energa, se genera toda la energa secundaria posible
para ese perodo de tiempo y adicionalmente los excedentes se vierten.
De esta manera se presentaran vertimiento en la planta cuando se presenten la
combinacin de las restricciones de volumen mximo del embalse y la turbinacin mxima
asociada al caudal de diseo. En la Figura 54 se muestra un ejemplo de la RNO para el
Proyecto Sogamoso.

Figura 54. Regla Normal de Operacin del Proyecto Sogamoso Operacin del embalse.

Para el caso del Proyecto Amoy la operacin a filo de agua, solo permite entregar la energa
producto de del caudal afluente sin la posibilidad de regulacin.
8.5.4 Confiabilidad en la Meta

Para calcular la energa generada, el modelo realiza los clculos a partir de una meta
mensual en volumen preestablecida y, de acuerdo a la confiabilidad calculada, se define
RNO - Proyecto Sogamoso
0
200
400
600
800
1000
1200
1400
1600
1800
2000
0 800 1600 2400 3200 4000 4800 5600 6400 7200 8000 8800
Volumen Disponible Hm
3
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l
(
G
W
h
-
m
e
s
)

Vertimientos

Meta
+
Ener
ga
Secu
ndari
a
Meta
(Energa firme)
Falla


Ingeniero Civil

76

aumentar o disminuir la meta con el fin de alcanzar la confiabilidad deseada. El proceso
es de tipo iterativo para evaluar la alternativa deseada, el cual permite ajustar la meta
para la confiabilidad media del 95%.
8.5.5 Representacin grfica de los resultados de la simulacin

Como resultados del modelo se presentan los valores de caudales mensuales afluentes,
turbinados y vertidos, adicionalmente los valores de energa generada diariamente y la
agregacin para obtener la energa media mensual. En la Figura 55, Figura 56, Figura 57
y Figura 58 se presentan los resultados para la simulacin realizada del Proyecto
Sogamoso.

Figura 55. Proyecto Sogamoso Caudal medio mensual, caudal de diseo y Operacin del embalse

Figura 56. Proyecto Sogamoso - Simulacin Energa Media y Firme Mensual

Variacin de Afluencias y Volumen del Embalse
0
500
1000
1500
2000
2500
3000
3500
4000
e
n
e
-
6
0
e
n
e
-
6
2
e
n
e
-
6
4
e
n
e
-
6
6
e
n
e
-
6
8
e
n
e
-
7
0
e
n
e
-
7
2
e
n
e
-
7
4
e
n
e
-
7
6
e
n
e
-
7
8
e
n
e
-
8
0
e
n
e
-
8
2
e
n
e
-
8
4
e
n
e
-
8
6
e
n
e
-
8
8
e
n
e
-
9
0
e
n
e
-
9
2
e
n
e
-
9
4
e
n
e
-
9
6
e
n
e
-
9
8
e
n
e
-
0
0
e
n
e
-
0
2
Fecha
A
f
l
u
e
n
c
i
a
s

-

V
o
l
.

E
m
b

(
H
m
3

)
Vol. Embalse Qdiseo Afluencias
Energa Media Mensual Generada (GWh-mes)
0
100
200
300
400
500
600
700
e
n
e
-
6
0
e
n
e
-
6
2
e
n
e
-
6
4
e
n
e
-
6
6
e
n
e
-
6
8
e
n
e
-
7
0
e
n
e
-
7
2
e
n
e
-
7
4
e
n
e
-
7
6
e
n
e
-
7
8
e
n
e
-
8
0
e
n
e
-
8
2
e
n
e
-
8
4
e
n
e
-
8
6
e
n
e
-
8
8
e
n
e
-
9
0
e
n
e
-
9
2
e
n
e
-
9
4
e
n
e
-
9
6
e
n
e
-
9
8
e
n
e
-
0
0
e
n
e
-
0
2
Fecha
E
n
e
r
g
a

(
G
W
h
-
m
e
s
)
Energa (GWh-mes) Efirme (GWh-mes)

Ingeniero Civil

77

Figura 57. Proyecto Sogamoso Relacin Potencia Promedio Mensual vs. Capacidad Instalada

Figura 58. Proyecto Sogamoso Energa Media Mensual

Relacin Potencia Promedio Mensual
Capacidad Instalada (MW)
0
100
200
300
400
500
600
700
800
900
Mes
P
o
t
e
n
c
i
a

(
M
W
)
Cap. Inst. Potencia
Energa Promedio Mensual
0
50
100
150
200
250
300
350
400
450
500
550
Mes
E
n
e
r
g
a

(
G
W
h
)

Ingeniero Civil

78

8.6 COMPARACIONES ENERGTICAS AISLADAS

El primer ejercicio de comparacin consisti en tratar de igualar cada Proyecto
Hidroelctrico con un parque elico que genere la misma energa media multianual, es
decir, encontrar la capacidad instalada que iguale la produccin energtica de cada
proyecto.

8.6.1 Proyecto Sogamoso Parque Elico Genrico 1.032 MW

Los resultados de las simulaciones energticas se presentan comparativamente en la
Figura 59.

Figura 59. Produccin Energtica Proyecto Sogamoso vs. Parque elico genrico que iguala su
Energa Media.

Segn el clculo realizado utilizando los supuestos en este trabajo, la produccin media
del proyecto hidroelctrico Sogamoso es de 4.635 GWh-ao con una potencia instalada
de 820 MW, igual energa en el perodo comn se alcanzara con un parque elico de
1.032 MW (26% de aumento de capacidad instalada, obtenida con 794 turbinas
N60/1.300kW). El rea necesaria para un parque elico de las proporciones descritas es
de cerca de 26 km
2
, si se considera la separacin que debe existir entre filas de
aerogeneradores para evitar la prdida de eficiencia por el fenmeno de estelas. Es de
anotar el tamao de este Parque Elico, el cual se presenta como una comparacin
indicativa y para el cul en Colombia no existe regulacin ni normatividad existente que
permita su entrada al mercado, ni mucho lo menos la inviabilidad financiera de un Parque
de esta proporcin.

El ciclo anual de las energas medias mensuales (GWh-mes) se presenta en la Figura 60,
donde se observa la complementariedad de los dos recursos.
Produccin Energtica Parque Elico Vs. Sogamoso
0
100
200
300
400
500
600
700
E
n
e
-
7
0
E
n
e
-
7
1
E
n
e
-
7
2
E
n
e
-
7
3
E
n
e
-
7
4
E
n
e
-
7
5
E
n
e
-
7
6
E
n
e
-
7
7
E
n
e
-
7
8
E
n
e
-
7
9
E
n
e
-
8
0
E
n
e
-
8
1
E
n
e
-
8
2
E
n
e
-
8
3
E
n
e
-
8
4
E
n
e
-
8
5
E
n
e
-
8
6
E
n
e
-
8
7
E
n
e
-
8
8
E
n
e
-
8
9
E
n
e
-
9
0
E
n
e
-
9
1
E
n
e
-
9
2
E
n
e
-
9
3
E
n
e
-
9
4
E
n
e
-
9
5
E
n
e
-
9
6
E
n
e
-
9
7
E
n
e
-
9
8
E
n
e
-
9
9
E
n
e
-
0
0
Mes
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
-
a
o
)
Parque Sogamoso

Ingeniero Civil

79

Figura 60. Comparacin de Produccin Energtica de Parque 1.032 MW vs. Sogamoso 820 MW

Si se considera los efectos del ENSO sobre las producciones energticas se observan los
comportamientos mostrados en las Figura 61 y Figura 62. Se anota que la para calcular la
produccin energtica para cualquiera de las dos fases del Evento ENSO, se
consideraron los ndices ONI histricos. La agregacin entonces se realiz a nivel
mensual multianual.

Figura 61 Ciclo Anual de Produccin Energtica del Proyecto Sogamoso 820 MW: Normal, El Nio
y La Nia

0
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0
50
100
150
200
250
300
350
400
450
500
E
n
e
r
g
a

M
e
d
i
a

A
n
u
a
l

(
G
W
h
)
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
)
Parque 1,032 MW Sogamoso E Media
Promedio Multianual Produccin Energtica Mensual
Sogamoso, Sogamoso ENSO - El Nio y Sogamoso ENSO - La Nia
0
50
100
150
200
250
300
350
400
450
500
550
600
Mes
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
-
m
e
s
)
Sogamoso ENSO-El Nio ENSO-La Nia

Ingeniero Civil

80

Figura 62 Ciclo Anual de Produccin Energtica del Parque 1032 MW: Normal, El Nio y La Nia

De la Figura 62 se destaca el hecho de que los eventos La Nia no afectan los primeros
meses del ao, enero, febrero, marzo y abril. En la Figura 61 es claro el aumento en la
generacin de Sogamoso durante La Nia.

Figura 63 Ciclo Anual de Produccin Energtica del Proyecto Sogamoso 820 MW y del Parque
1.032 MW: Evento ENSO El Nio

Adicionalmente se observa el incremento en la Energa Media Anual del Parque Elico con
respecto al Proyecto Sogamoso.
Parque, Parque ENSO - El Nio y Parque ENSO - La Nia
0
50
100
150
200
250
300
350
400
450
500
550
Mes
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
-
m
e
s
)
Parque ENSO-El Nio ENSO-La Nia
0
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0
50
100
150
200
250
300
350
400
450
500
550
E
n
e
r
g
a

M
e
d
i
a

A
n
u
a
l

(
G
W
h
)
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
)
Parque 1,032 MW Sogamoso
E Media Sogamoso E Media Parque 1,032

Ingeniero Civil

81

Figura 64 Ciclo Anual de Produccin Energtica del Proyecto Sogamoso 820 MW y del Parque
1040 MW: Evento ENSO La Nia

En la Figura 63 y Figura 64 se presentan las comparaciones a nivel mensual para cada
fase del ENSO. Mientras la produccin energtica de Sogamoso en los aos El Nio y La
Nia con relacin a los aos normales se mantiene, la produccin energtica del parque
disminuye si se compara la generacin durante El Nio y La Nia. Lo anterior obedece a
que el Proyecto Sogamoso tiene la capacidad de regular su recurso afluente, diferente al
Parque genrico de 1.032 MW que opera como si fuera una planta a filo de agua, dado a
que no tiene como almacenar o regular el viento.

Con base en lo anterior se estim que durante El Nio, la disminucin de la generacin
por parte del Proyecto Sogamoso ser del 5%, porcentaje igual al aumento en produccin
energtica del Parque 1.032 MW y durante La Nia el Proyecto Sogamoso aumentar su
produccin del orden del 11,2%, y el Parque 1.032 MW disminuir su produccin en un
1,1%.

La Figura 65 presenta las diferencias mensuales de las producciones energticas del
Parque 1.032 MW con respecto al Proyecto Sogamoso. Se observa que durante El Nio
hay aumentos en las diferencias de generacin a favor del Parque para los meses de
Febrero, Marzo, Abril, Junio, Julio y Diciembre y durante La Nia el Proyecto Sogamoso
genera ms que el Parque y especficamente en los meses de Abril, Mayo, Agosto,
Septiembre, Octubre y Noviembre.
0
1000
2000
3000
4000
5000
6000
0
50
100
150
200
250
300
350
400
450
500
550
600
E
n
e
r
g
a

M
e
d
i
a

A
n
u
a
l

(
G
W
h
)
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
)
Parque 1,032 MW Sogamoso
E Media Sogamoso E Media Parque 1,032

Ingeniero Civil

82

Figura 65 Diferencias de Produccin Energtica Parque 1.040 MW vs. Sogamoso 820 MW

Los anteriores anlisis indican que el Parque Elico de 1.040 MW apoyara y ayudara a
afirmar la produccin energtica del proyecto Sogamoso en los meses crticos de la
ocurrencia de El Nio. Se debe aclarar que las condiciones regulatorias y de mercado de
un Parque Elico de este tamao an no estn establecidas y que lo presentando sera
en un escenario donde se brinden oportunidades de participacin (despacho) similares a
la tecnologa elica actual.

8.6.2 Proyecto Amoy Parque Elico Genrico 140 MW

Se realiza anlisis con un proyecto sin embalse regulador como el caso de Amoy. Los
resultados de las simulaciones energticas se presentan comparativamente en la Figura
66.

Figura 66 Produccin Energtica Proyecto Amoy vs. Parque elico genrico que iguala energa
Media

Segn el clculo realizado, la produccin media del Proyecto Hidroelctrico Amoy es de
590 GWh-ao con una potencia instalada de 80 MW, igual energa en el perodo comn
se alcanzara con un parque elico de 140 MW, lo que indica un 75% de aumento en
capacidad instalada, correspondiente a 108 turbinas N60/1300kW. El factor de planta de
Diferencias Parque - Sogamoso. Normal Vs. Evento ENSO - La Nia
-400
-300
-200
-100
0
100
200
Mes
D
i
f
e
r
e
n
c
i
a

E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
-
m
e
s
)
Normal ENSO-La Nia ENSO-El Nio
Produccin Energtica Parque Elico Vs. Amoy
0
10
20
30
40
50
60
70
80
90
100
e
n
e
-
7
4
e
n
e
-
7
5
e
n
e
-
7
6
e
n
e
-
7
7
e
n
e
-
7
8
e
n
e
-
7
9
e
n
e
-
8
0
e
n
e
-
8
1
e
n
e
-
8
2
e
n
e
-
8
3
e
n
e
-
8
4
e
n
e
-
8
5
e
n
e
-
8
6
e
n
e
-
8
7
e
n
e
-
8
8
e
n
e
-
8
9
e
n
e
-
9
0
e
n
e
-
9
1
e
n
e
-
9
2
e
n
e
-
9
3
e
n
e
-
9
4
e
n
e
-
9
5
e
n
e
-
9
6
e
n
e
-
9
7
e
n
e
-
9
8
e
n
e
-
9
9
e
n
e
-
0
0
Mes
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
-
a
o
)
Parque Amoy

Ingeniero Civil

83

Amoy, la buena regulacin natural de su caudal y la menor complementariedad con
respecto a la planta elica hacen que se requiera cerca del doble de su capacidad
instalada para obtener una energa media comparable.

El rea necesaria para un parque elico de las proporciones descritas es del orden de 3,5
km
2
.

Figura 67 Comparacin de Produccin Energtica de Parque 140 MW vs. Amoy 80 MW

Considerando lo anterior, se estim que durante El Nio la disminucin de la generacin
por parte del Proyecto ser del 5%, mientras que el aumento en produccin energtica del
Parque de 140 MW es del 3,8% y durante La Nia el Proyecto aumentar su produccin
del orden del 9,3 %, y el Parque 140 MW disminuir su produccin solo en un 1,1%. Ver
Figura 67.

8.7 ANLISIS ENERGTICO INTEGRADO DE COMPLEMENTARIEDAD

Este anlisis considera una operacin coordinada entre proyectos y supone que la
energa generada por el Parque Elico genrico entra sobre la base para meta de
generacin de la simulacin de los proyectos hidroelctricos, es decir que supone que
cuando se est generando en un Parque Elico se est almacenando agua en el
embalse, para el caso de los proyectos hidroelctricos con regulacin.

Para este anlisis integrado se consider que la energa media mensual generada por los
parques elicos se convierte en caudales mensuales utilizando el factor de conversin
medio de la central hidroelctrica considerada. Luego estos caudales entran en la base de
la afluencia de la Central Hidroelctrica disminuyendo o facilitando el cumplimiento de la
meta mensual, a exigir por parte de la RNO del modelo de simulacin sobre los caudales
propios de la Central.

0
100
200
300
400
500
600
700
0
10
20
30
40
50
60
70
E
n
e
r
g
a

M
e
d
i
a

(
G
W
h
)
E
n
e
r
g
a

M
e
d
i
a

M
e
n
s
u
a
l

(
G
W
h
)
Parque 140 MW Amoy E. Media

Ingeniero Civil

84

8.7.1 Anlisis Integrado P. H. Sogamoso Parques Elicos Genricos

Por lo anterior, se tom el proyecto Sogamoso y se incluy la generacin de Parques
Elicos Genricos aumentando su capacidad instalada considerando 20 MW, 50 MW, 100
MW, 200 MW y 400 MW, con el fin de identificar si hay aumentos o disminuciones en la
energa media del conjunto, y verificar si la energa firme anual integrada es mayor que la
suma de las energas firmes anuales aisladas.

Tabla 11 Caractersticas del Proyecto Sogamoso para la simulacin integrada
Caractersticas Proyecto Hidroelctrico Sogamoso
Volumen til Embalse 2.801
Volumen Inicial del Embalse 1.400
Factor de Conversin promedio (MW/m
3
/s) 1,22
Eficiencia Turbina (%) 92,0%
Eficiencia Generador (%) 98,0%
Eficiencia Transformador (%) 99,5%
Salto Bruto (m) 150,5
Perdidas (m) 9,0
Salto Neto (m) 141,5
Caudal de diseo (m
3
/s) 660
Capacidad Instalada (MW) 821,0
Capacidad mxima de Turbinacin (Hm
3
) 1.733,5

Las simulaciones se realizaron maximizando la energa firme, mediante la variacin de la
meta en caudal impuesta al modelo, con el objeto de realizar las comparaciones con las
diferentes capacidades de Parques Elicos.

En la Figura 68 se observa que la mayor energa media se alcanza para las metas
menores con valores por debajo del caudal medio del ro, esto aumenta la posibilidad de
almacenar en el embalse, mantener el nivel arriba mejorando la eficiencia y produccin
energtica. La mayor energa generada es secundaria, aunque esto aumenta los
vertimientos de la Central, sin embargo, la energa firme, para una confiabilidad del 95%,
se consigue aumentando la exigencia en meta sobre la Central y disminuyendo muy poco
la energa media, lo cual es de poca importancia ante la ganancia en energa firme.


Ingeniero Civil

85

Figura 68 Simulacin de Generacin de Energa del sistema Integrado variacin meta o demanda
vs. energa generada

En la Figura 69 se presenta la variacin de la energa media y firme integrada para
diferentes tamaos de parques Elicos. La simulacin se realiz maximizando energa
firme.

Figura 69 Simulacin de Generacin de Energa variacin meta o demanda vs. energa generada
(MW)

En la operacin integrada, a medida que se aumenta la capacidad instalada de parques
elicos genricos la energa firme aumenta significativamente, mientras que la energa
media es muy similar a la suma de las energas medias de la operacin aislada, y para
tamaos mayores a 400 MW la energa media integrada es menor; lo anterior se debe a
que para los mayores tamaos de capacidad instalada de parques, aumentan los
vertimientos del embalse del Proyecto Sogamoso y disminuye la energa secundaria.

1500
2000
2500
3000
3500
4000
4500
5000
200 220 240 260 280 300 320 340 360 380 400 413
Meta Q (m3/s)
E
n
e
r
g
a

A
n
u
a
l

(
G
W
h
-
a
o
)
EMEDIA
E FIRME
0
1000
2000
3000
4000
5000
6000
7000
0 20 MW 50 MW 100 MW 200 MW 400 MW
E
n
e
r
g
a

G
W
h
-
a
o
Tamao de Parques Elicos Genricos
Firme Integrado (Parque +Sogamoso)
Firme Parque + Firme Sogamoso
Media Integrado (Parque +Sogamoso)
Media Parque + Media Sogamoso

Ingeniero Civil

86

Como se mencion anteriormente, la energa firme integrada es mayor que la suma de las
energas firmes aisladas. Los valores de la simulacin integrada se observan en la Tabla
12.

Tabla 12 Resultados de la simulacin Integrada Sogamoso Parques Elicos Genricos
Parmetro 0 MW 20 MW 50 MW 100 MW 200 MW 400 MW
Meta de mxima Energa Firme (m
3
/s) 413 425 442 465 512 605
Meta de mxima Energa Firme (Hm
3
) 1.084,8 1.116,3 1.160,9 1.221,4 1.344,8 1.589,1
Confiabilidad (%) 95% 95% 95% 95% 95% 95%
Salto Neto Promedio (m) 140,2 140,0 139,8 140,0 140,2 140,5
Descarga promedio (m
3
/s) 454 463 476 496 537 612
Energa media mensual generada (GWh-mes) 386 393 403 420 455 518
Vertimientos promedios mensuales (m
3
/s) 29,45 29,83 30,80 33,05 38,32 54,29
Energa media anual (GWh) 4.635 4.712 4.838 5.046 5.456 6.218
Factor de Planta Medio 0,64 0,66 0,67 0,70 0,76 0,86
Energa Firme (GWh-ao) 3.564 3.631 3.751 3.977 4.425 5.244
Factor de Planta Firme 0,50 0,50 0,52 0,55 0,62 0,73
Caudal promedio mensual afluente (m
3
/s) 484,33 484,33 484,33 484,33 484,33 484,33
Afluencias Promedios (Hm
3
) 1.272,12 1.295,53 1.333,00 1.392,32 1.512,52 1.752,92
Capacidad Instalada Elica (MW) 20 50 100 200 400
Nmero de Aerogeneradores 15 39 77 154 308
Energa Media Parques Genricos (GWh-ao) 81,9 212,9 420,3 840,5 1681,0
Ene. Media Parque + Ene. Media Sog. (GWh-ao) 4.716,89 4.847,89 5.055,29 5.475,51 6.316,03
Diferencia Energa Media (GWh-ao) 4,86 9,71 9,68 19,62 97,75
Energa Firme Parques Genricos (GWh-ao) 36,19 94,09 185,77 371,53 743,07
Delta Energa Firme (GWh-ao) 67,13 187,54 413,34 861,39 1.680,48
Firme Parque + Firme Sogamoso (GWh-ao) 3.599,7 3.657,6 3.749,3 3.935,1 4.306,6
Diferencia de Energa Firme (GWh-ao) 30,9 93,5 227,6 489,9 937,4

En la Tabla 12 se presentan las energas firmes de los Parques Elicos, las cuales fueron
calculadas con las curvas de duracin de las series de energas medias mensuales, y
tomando la asociada al 95%. Se resalta que a la fecha en Colombia no existe
reglamentacin aplicable a los Parques Elicos con relacin al clculo de energa firme,
como la Energa Firme para el Cargo por Confiabilidad ENFICC.

El trabajo de Vergara et al. (2010) incluye una propuesta de clculo de este parmetro
para plantas elicas, basada en la generacin de energa del primer ao de operacin, y
ajustada anualmente considerando el perodo base inicial, este planteamiento considera
la poca informacin de velocidad del viento disponible y procura defender la importancia
de asignarle algn ingreso constante que ayude a viabilizar las plantas de generacin con
esta tecnologa.

Teniendo en cuenta lo expuesto anteriormente, y a manera de ejercicio, se utiliz el
software HIDENFICC desarrollado por la CREG para Centrales Hidroelctricas, y a partir
de la informacin de un Parque Genrico de 20 MW se clculo una Energa Firme para el
Cargo por Confiabilidad - ENFICC de dicho Parque Elico. Las consideraciones tienen en

Ingeniero Civil

87

cuenta la velocidad del viento promedia mensual como ingreso, las restricciones mximas
y mnimas de operacin de las turbinas y un factor de conversin calculado a partir de la
informacin a nivel horario utilizado en este informe.

Los resultados arrojaron una ENFICC del 95% PSS (es decir el nivel de energa que
Central Hidroelctrica puede generar con una probabilidad del 95%) igual a 20 GWhao,
la cual es inferior al valor calculado con curva de duracin de energas medias mensuales
e igual a 36,2 GWh-ao.

Como se observa en la Tabla 12, la energa firme aumenta significativamente para
parques mayores o iguales a 50 MW, lo que en primera medida indicara ingresos
adicionales que integralmente mejoraran las condiciones financieras de un Parque Elico
operando aisladamente.
8.8 ANLISIS FINANCIERO DE LA GENERACIN INTEGRADA

Con la informacin de generacin anual media y firme integrada, y a partir de las
generaciones aisladas del Proyecto Sogamoso y los diferentes Parques Elicos, se
cuantific el beneficio financiero a partir de una evaluacin con Valor presente Neto - VPN
y Tasa Interna de Retorno TIR.

Para lo anterior se utiliz la siguiente informacin:

Capacidad Instalada Parques Genricos: 20 MW, 50 MW, 100 MW, 200 MW y 400
MW.

Capacidad Instalada Proyecto Sogamoso: 820 MW

Costo de instalacin Parque Elico: 1.400 USD/kW (Vergara et. al, 2010).

Costo de instalacin Sogamoso: 1.707 USD/kW
5

Costo AOM Parques Elicos: 25 USD/kW instalados.

Otros costos: 2 USD/kW

Tasa de descuento para evaluacin financiera: 10%

Precio Energa Firme: 42 USD/MWh

Precio Energa Secundaria: 34 USD/MWh (80% del precio firme).

5
Informacin extrada de la presentacin ISAGEN dentro del Foro - LAWEA: "Elementos para el desarrollo de la energa
Elica en Colombia", julio de 2009.


Ingeniero Civil

88

Vida til Parques Elicos: 20 aos.

Vida til Sogamoso: 50 aos.

La inversin se considera puntual en el ao 0 de evaluacin.

No se consideran depreciaciones ni de obras civiles ni de equipos.

Como proyecto aplicable a MDL, se consider la venta de Certificados de
Reduccin de Emisiones de gases de efecto invernadero con un factor de emisin
de 0,2849 kg CO
2
e/kWh (Resolucin Ministerio de Minas y Energa No. 180947 de
junio de 2010), y una Tarifa de venta de 18 USD/tonCO
2
e (Vergara et al., 2010).
Con la informacin anterior, se tiene que el costo de inversin de los Parques Elicos va
desde 36 Millones USD para un parque de 20 MW hasta 720 Millones USD para un
parque de 400 MW. Para el caso del Proyecto Hidroelctrico Sogamoso se tiene una
inversin de 1.400 Millones USD.

Para los costos por Administracin, Operacin y Mantenimiento (AOM) de parques
elicos, en Vergara et al., (2010) se presenta para los EEUU su variacin en los ltimos
aos pasando de 35 USD/MWh en la dcada de los aos 80, a 20 USD/MWh en los aos
90, hasta 9 USD/MWh al da de hoy.

Para el caso del Parque Elico Jeprachi, el costo AOM igual al 2% de su inversin (28
Millones USD), lo que implica un costo por AOM de 28 USD/kW instalados. Teniendo en
cuenta lo anterior se opt por un valor intermedio del orden 25 USD/kW.
6

Para el precio de energa firme se utiliz la informacin del Plan de Expansin de la
UPME en el perodo 2009 2023 el cual sugiere que los precios del mercado se
comportarn como se presenta en la Figura 70. Para el escenario de demanda media se
calcul un valor promedio de cerca de 42 USD/MWh.

6
Foro Contexto y Normalizacin de la energa elica y solar en Colombia, UPME ICONTEC. Bogot, Noviembre 2008.
Experiencias en el Desarrollo de Proyectos de Energa Elica Parque Elico Jeprachi EPM.


Ingeniero Civil

89

Figura 70 Variacin del precio de energa entre 2009 y 2023. Adaptada de UPME, 2009.

A partir de los resultados de la simulacin descritos en la Tabla 12 y las consideraciones
anteriores, en la Figura 71 se observa el comportamiento del VPN y TIR para los Parques
Elicos Genricos, en los cuales no existen ingresos que superen los costos de inversin,
operacin y mantenimiento en el tiempo, lo cual tambin se ve reflejado en una TIR menor
que la tasa de descuento de referencia para la evaluacin considerada es 10%.

Figura 71 Variacin VPN y TIR para los diferentes Parques Elicos Genricos

Para el caso del Proyecto Sogamoso, los resultados de la operacin aislada e integrada
se presentan en la Figura 72, obtenindose un beneficio hasta de 2% en la TIR con un
parque de 400 MW.
0
15
30
45
60
75
90
105
120
135
e
n
e
-
0
9
j u
l -
0
9
e
n
e
-
1
0
j u
l -
1
0
e
n
e
-
1
1
j u
l -
1
1
e
n
e
-
1
2
j u
l -
1
2
e
n
e
-
1
3
j u
l -
1
3
e
n
e
-
1
4
j u
l -
1
4
e
n
e
-
1
5
j u
l -
1
5
e
n
e
-
1
6
j u
l -
1
6
e
n
e
-
1
7
j u
l -
1
7
e
n
e
-
1
8
j u
l -
1
8
e
n
e
-
1
9
j u
l -
1
9
e
n
e
-
2
0
j u
l -
2
0
e
n
e
-
2
1
j u
l -
2
1
e
n
e
-
2
2
j u
l -
2
2
e
n
e
-
2
3
j u
l -
2
3
U
S
$
/
M
W
h
DALTA DMEDIA DBAJA
Parque Elico - Variacin VPN y TIR: Operacin Aislada
$ -300.000.000
$ -250.000.000
$ -200.000.000
$ -150.000.000
$ -100.000.000
$ -50.000.000
$ 0
20 50 100 200 400
V
P
N

(
U
S
D
)
2,0%
2,5%
3,0%
3,5%
4,0%
Capacidades Instaladas (MW)
T
I
R

(
%
)
VPN TIR

Ingeniero Civil

90

Figura 72 Proyecto Sogamoso: Variacin VPN y TIR, operacin integrada con beneficios.

Teniendo en cuenta la diferencia de energas medias y firmes obtenidas de la operacin
integrada, los ingresos adicionales obtenidos se incluyeron en los ingresos anuales de los
Parques Genricos y se le restaron al Proyecto Sogamoso con el fin de asignarle la
ganancia adicional de la operacin integrada solo a los Parques Genricos . Los
resultados se presentan en la Figura 73 y Figura 74.

Figura 73 Proyecto Sogamoso: Variacin VPN y TIR, operacin integrada sin beneficios

En la Figura 74 se observa que para los parques de 100 MW, 200 MW y 400 MW se
presentan VPN mayores que cero y TIR por encima de la tasa de descuento. Para las
consideraciones financieras supuestas, para tamaos menores la energa firme adicional
no alcanza a cubrir los costos asociados aunque la energa secundaria no vare
significativamente. Adicionalmente para parques mayores a 400 MW se prev que aunque
haya un aumento de energa firme la energa secundaria perdida disminuir los beneficios
obtenidos en la operacin integrada.
Proyecto Sogamoso - Variacin VPN y TIR: Operacin aislada e Integrada
con Parques Elicos Genricos
$ 0
$ 50.000.000
$ 100.000.000
$ 150.000.000
$ 200.000.000
$ 250.000.000
$ 300.000.000
$ 350.000.000
$ 400.000.000
0
2
0
5
0
1
0
0
2
0
0
4
0
0
V
P
N

(
U
S
D
)
10,0%
11,0%
12,0%
13,0%
14,0%
15,0%
T
I
R

(
%
) VPN TIR
Proyecto Sogamoso - Variacin VPN y TIR: Operacin aislada e Integrada
con Parques Elicos Genricos
$ 0
$ 20.000.000
$ 40.000.000
$ 60.000.000
$ 80.000.000
$ 100.000.000
$ 120.000.000
$ 140.000.000
0
2
0
5
0
1
0
0
2
0
0
4
0
0
V
P
N

(
U
S
D
)
11,0%
11,2%
11,4%
11,6%
11,8%
12,0%
T
I
R

(
%
)
VPN TIR

Ingeniero Civil

91

Figura 74 Variacin VPN y TIR para los diferentes Parques Elicos Genricos con beneficios de
operacin integrada.

Es de aclarar que existen desventajas de la operacin conjunta en cuanto a la variacin
itraanual de los precios del mercado, los cuales tienen generalmente una alta correlacin
al comportamiento del ciclo hidrolgico nacional. Al operar conjuntamente la energa
mensual del sistema integrado se comporta al sistema hdrico convencional, caso
contrario cuando se operan los Parques Elicos por separado, los cuales entregan su
mayor energa mensual cundo el precio est ms alto, durante las pocas de sequa en
ciclo hidrolgico nacional. Los clculos indican una correlacin de 0,3 entre las energas
medias mensuales de un Parque de 20 MW con la serie de precios tomados del Plan de
Expansin de referencia, mientras que para la serie de energas medias mensuales de la
operacin integrada con igual capacidad instalada elica es de -0,5.
Parque Elico - Variacin VPN y TIR: Operacin Integrada
$ -10.000.000
$ -5.000.000
$ 0
$ 5.000.000
$ 10.000.000
$ 15.000.000
$ 20.000.000
20 50 100 200 400
V
P
N

(
U
S
D
)
7,0%
7,5%
8,0%
8,5%
9,0%
9,5%
10,0%
10,5%
11,0%
T
I
R

(
%
)
VPN TIR

Ingeniero Civil

92

9 CONCLUSIONES Y RECOMENDACIONES

La Guajira Colombiana presenta un ciclo diurno claramente definido con los mximos
entre las 15:00 y 16:00 horas, y los mnimos en horas de la madrugada (03:00 05:00),
justo antes de que el sol empiece a salir. Lo anterior se debe a la variacin de la
temperatura horaria debido a la insolacin asociada al clima tropical de esta regin.
El ciclo anual de la velocidad del viento est relacionado entre otros, con los perodos de
precipitacin y sequa directamente influenciadas por la Zona de Convergencia
Intertropical - ZCIT, la cual acta sobre todo el territorio nacional. Adicional a esto, dada la
posicin de la zona de estudio (norte de Colombia) la ZCIT solo pasar una vez al ao, lo
cual se refleja en un ciclo anual unimodal en la velocidad del viento y la precipitacin
promedios mensuales en la zona. De igual manera, existe una relacin entre el aumento
de la velocidad del viento, producto del aumento de la temperatura y, aunque la
variabilidad sea mayor en la velocidad, las dos variables presentan mximos justo antes
de empezar el perodo lluvioso, este comportamiento se aprecia en los meses de junio,
julio y agosto.
La importancia de la Oscilacin del Atlntico Norte impacta tanto en el anlisis de
correlacin mensual, como en el trimestral para la velocidad del viento y para la
temperatura, se observa la influencia del ndice en el trimestre MAM con el mismo
trimestre de la variable y antecedindola un trimestre, es decir MAM-JJA. Se debern
efectuar con informacin de ms estaciones, adicional a la estacin Puerto Bolvar con un
amplio perodo de registros para verificar si el comportamiento influye fuertemente en la
Zona.
Los datos virtuales del Reanlisis (NCEP/NCAR) se ajustan en buena medida a la
informacin real recopilada en superficie. A futuro se podran hacer correlaciones de los
datos del Reanlisis con los ndices del ENSO y la NAO y anlisis espectral de la variable,
previendo resultados importantes para la extrapolacin de series que permitan la
evaluacin en interanual de Parques Elicos.
En relacin a la direccin de los vientos en la Guajira, se encontr un altsimo porcentaje
de la direccin de la velocidad del viento dentro del Cuadrante NE E, correspondiente a
los 45 90 de Azimut, esto indica la predominancia de la componente zonal del viento,
claramente explicada por la presencia de los vientos Alisios del este. El anlisis no revel
cambios significativos intranuales en la direccin del viento, sin embargo, se encontr
que a medida que el anlisis se desplaza hacia la zona sur de la Guajira la componente
meridional toma mayor importancia. La importancia de tener una direccin predominante
radica en la eficiencia energtica de los parques elicos a instalar en la zona, dado a que
los cambios de direccin disminuyen la generacin y aumentan las prdidas por estelas.
El comportamiento de la velocidad vertical del viento se ajusta al comportamiento
exponencial descrito por la teora del Hellmann. Los coeficientes calculados entre 0,19 y
0,2 utilizaron para escalar las series de viento a la altura 60 metros, necesaria para
calcular las energas requeridas.


Ingeniero Civil

93

Las series de caudales medios mensuales de los Proyectos y Centrales Hidroelctricas
analizadas son impactados por los Eventos ENSO, El Nio y La Nia. Esto se identific
con la ayuda de coeficiente de correlacin con la serie del ndice ONI y con anlisis
espectral.

Con la ayuda del anlisis espectral se visualizaron picos para aproximadamente 4 aos, lo
que podra indicar impacto en la zona del Fenmeno de El Nio y por lo tanto aumento de
velocidades. Tanto la estacin Puerto Bolvar como la estacin virtual del Reanlisis
presentaron dichos aumentos en la poca de mejores vientos dentro del ciclo anual de la
zona JJA y corroborados por las correlaciones trimestrales.

Las mayores correlaciones trimestrales se observan cuando la velocidad es mayor, en el
trimestre JJA y para similar trimestre del ndice ONI (valores de coeficientes entre 0,43 y
0,50), a pesar que en este perodo los embalses del SIN an cuentan con agua disponible
para generacin y que la correlacin ms relevante para el anlisis de complementariedad
obedecera al trimestre DEF, en el trimestre JJA se presentan valores altos en el precio
de la energa. El anterior comportamiento se verific en las estaciones de Puerto Bolvar y
Reanlisis.

Adicional al punto anterior, se verific a partir de la informacin de vientos de las
estaciones de Puerto Bolvar, Reanlisis y Ernesto Cortizzos, que la velocidad del viento
aumenta durante El Nio entre los meses de agosto a noviembre (aumento hasta en 12%)
y disminuye durante los meses de enero a mayo; resultados contrarios a los encontrados
cuando se presenta La Nia, donde la velocidad del viento aumenta entre enero y junio en
un 4% con relacin a los valores promedios. Este aspecto impacta sobre la
complementariedad de los recursos ante la ocurrencia de El Nio, dado a que en el primer
trimestre del ao calendario, cuando an est presente el Evento y los embalses tengan
menor agua disponible para generacin, la velocidad del viento es menor.

La importancia de conocer la variacin del ciclo diurno de la velocidad del viento durante
las dos fases del ENSO, radica en los cambios de disponibilidad energtica de un parque
elico en el sitio, si bien es claro, en condiciones normales los mximos se presentarn a
mitad de la tarde aproximadamente a las 15:00 horas, durante La Nia, las mayores
velocidades se presentarn en horas de la noche cuando aumenta la demanda horaria de
energa elctrica en el pas. Lo anterior podra estudiarse dentro de un contexto de
anlisis de complementariedad intradiaria de los recursos en cuestin y desde el punto de
vista de generacin el comportamiento de la demanda diaria del SIN y la respectiva oferta
presentada por la generacin diaria integrada.
Las correlaciones inversas entre las series de viento de la estacin de puerto Bolvar,
elegida por su amplio perodo de registro y las series de caudales de distintas regiones
del pas, disminuyen a medida que el sitio de inters se encuentra ubicado ms hacia el
sur, lo cual indica una asociacin directa con la ZCIT. La mayor complementariedad entre
los recursos se presenta para los ros de la zona nororiente del pas, sin embargo, en el
centro del pas algunas centrales hidroelctricas existentes o futuros proyectos, podran
complementarse de manera adecuada con parques elicos en La Guajira.


Ingeniero Civil

94

El clculo de la produccin energtica a partir de la serie de vientos a nivel horario y
extendido con informacin mensual de estaciones como Puerto Bolvar y el nodo del
Reanlisis NCEP / NCAR provee informacin que se ajusta al Parque Elico Jeprachi
(utilizando el mismo aerogenerador Nordex N60/1.300), existente en la zona
considerndose en primera medida una aproximacin aceptable, a pesar que no
considera muchos elementos que utilizados por los Software para modelacin de Parques
Elicos como el WindPro y el WASP. Cabe anotar que la informacin disponible en la
base de datos de Nen de XM, con relacin a la produccin energtica de Jeprachi, no
incluye el nmero de mquinas asociadas a dicha produccin; y que los aerogeneradores
modernos son de mayores capacidades instaladas y ms eficientes.

El anlisis de las series histricas de produccin de energa de las centrales
hidroelctricas de ISAGEN, no presenta una correlacin directa con la produccin
energtica de Parques Elicos genricos de igual capacidad instalada, incluso en el caso
de la Central Miel I, la ubicada ms al sur de las Centrales analizadas, registra
prcticamente un comportamiento similar. Lo anterior podra obedecer a que el mercado
opera bajo una constante asociado con las centrales trmicas, las cuales operan cuando
se les necesite, independientemente de la poca del ao; y el comportamiento de la
hidrologa, es decir un mercado que no incluye en grandes cantidades otro recurso
complementario.

El anlisis de complementariedad entre la series de vientos de Puerto Bolvar y las series
de caudales afluentes a las centrales y proyectos hidroelctricos de ISAGEN, a nivel
mensual, report que: para la serie del ro Amoy se presenta un comportamiento
diferente a las dems series de caudales, presentando valores de caudales altos cuando
se presenta el verano en las dems series de caudales y especficamente en los meses
de Junio y Julio. Esto tambin se ve reflejado en el coeficiente de correlacin igual a cero
(0) de dicha serie con la serie de vientos de la estacin de Puerto Bolvar.
La serie del Reanlisis presenta mayores correlaciones con los caudales, a excepcin de
la serie de Miel I. No existen datos disponibles para calificar este comportamiento; sin
embargo, la serie del Reanlisis disminuye su correlacin con los caudales mensuales de
norte a sur, en relacin con los sitios de inters; es decir, Sogamoso que es el sitio ms al
norte tiene la mayor correlacin -0,55 y Amoy -0,12, la que se encuentra ubicada ms al
sur de los sitios de inters. Lo anterior verifica el comportamiento encontrado en el
anlisis de complementariedad regional.

Los factores de planta (o capacidad) histricos de las centrales hidroelctricas de ISAGEN
consideran la capacidad mxima de la planta sin incluir disponibilidad y criterios de
operacin comercial. Para el caso de la central Miel I 0,4, Jaguas 0,47 y San Carlos 0,52;
estos ndices se calcularon con la informacin disponible en el portal del administrador del
Mercado de Energa; XM. Lo anterior indica que los factores de planta de operacin son
menores a los proyectados en las simulaciones energticas, lo que se debe tener en
cuenta a la hora de realizar las cuantificaciones energticas y financieras cuando se
realizan estudios de factibilidad de tanto de Proyectos Hidroelctricos como Parques
Elicos. Caso similar se refleja en las producciones de Jeprachi y el Parque Elico
localizado en Puerto Bolvar.

Ingeniero Civil

95

El anlisis aislado de Parques Elicos genricos que producen la misma energa media
de los proyectos hidroelctricos Sogamoso y Amoy refleja complementariedad en los
momentos dnde los perodos de verano afectan la hidrologa afluente, incluso para el
caso de Sogamoso que tiene embalse regulador. Sin embargo para alcanzar igual energa
media se requieren Parques con un 26% ms capacidad instalada. Para el caso de un
parque que iguale la produccin energtica media anual de Sogamoso se requiere un
rea de 26 km
2
, si se considera la separacin que debe existir entre filas de
aerogeneradores para evitar la prdida de eficiencia por el fenmeno de estelas.

Mientras la produccin energtica de Sogamoso en los aos El Nio y La Nia con
relacin a los aos normales se mantiene, la produccin energtica del parque disminuye
si se compara la generacin durante El Nio y La Nia. Lo anterior obedece a que el
Proyecto Sogamoso tiene la capacidad de regular su recurso afluente, diferente al Parque
genrico de 1.032 MW que opera como si fuera una planta a filo de agua, dado a que no
tiene como almacenar o regular el viento. La disminucin de la generacin por parte del
Proyecto ser del 5%, porcentaje igual al aumento en produccin energtica del Parque
1.032 MW y durante La Nia el Proyecto aumentar su produccin del orden del 11,2%, y
el Parque 1.032 MW disminuir su produccin en un 1,1%. Lo anterior indica que el
Parque Elico de 1.032 MW ayudara a la produccin energtica del proyecto Sogamoso
en los meses crticos de la ocurrencia de El Nio.

La simulacin integrada de parques elicos genricos de diferentes tamaos: 20 MW, 50
MW, 100 MW, 200 MW y 400 MW con un proyecto con embalse regulador como el caso
de Sogamoso permitira obtener una mayor energa firme dado a que la generacin elica
entra en la base y baja la meta exigida al proyecto aumentando dicha energa, sin
embargo disminuye su energa secundaria y an ms para Parques Elicos de mayor
tamao. Lo anterior reflejado en beneficios financieros ayuda a los Parques Elicos dado
que los ingresos adicionales por energa firme, normalmente valoradas a un mayor precio,
incrementan los ingresos al Parque facilitando su viabilidad.

A pesar de que las evaluaciones financieras no consideran los detalles de depreciacin,
perodo de construccin, costos de estudios, costos detallados de operacin dentro del
SIN y dems, a grandes rasgos se observa que la complementariedad aporta beneficios
que bajo una simulacin con tamaos representativos de Parques Elicos aumenta la
viabilidad facilitando la futura inclusin de dichos Parques dentro del Mercado. Lo anterior
sugiere que la regulacin y normatividad existente deber evolucionar hacia la
participacin de las fuentes renovables no convencionales dentro de la canasta
energtica.

Se encontr que para los parques de 100 MW, 200 MW y 400 MW se presentan VPN
mayores que cero y TIR por encima de la tasa de descuento. Para las consideraciones
financieras supuestas, para tamaos menores la energa firme adicional no alcanza a
cubrir los costos asociados, aunque la energa secundaria no vare significativamente.
Adicionalmente para parques mayores a 400 MW se prev que aunque haya un aumento
de energa firme la energa secundaria perdida disminuir los beneficios obtenidos en la
operacin integrada.

Ingeniero Civil

96

Existen desventajas de la operacin conjunta en cuanto a la variacin itraanual de los
precios del mercado, los cuales se asocian completamente al comportamiento del ciclo
hidrolgico nacional, dado el alto componente de 66% de capacidad hidrulica del SIN. Al
operar conjuntamente la energa mensual del sistema integrado se comporta al sistema
hdrico convencional, caso contrario cuando se operan los Parques Elicos por separado,
los cuales entregan su mayor energa mensual cundo el precio est ms alto, durante las
pocas de sequa en ciclo hidrolgico nacional. Los clculos indican una correlacin de
0,3 entre las energas medias mensuales de un parque elico de 20 MW con la serie de
precios tomados del Plan de Expansin de referencia, mientras que para la serie de
energas medias mensuales de la operacin integrada con igual capacidad instalada
elica es de -0,5.

Para trabajos futuros es recomendable incluir dentro de los anlisis la operacin dentro
del Mercado de Energa, variando las capacidades instaladas como se present en este
informe. Se considera de inters que se tengan en cuenta todos los supuestos para que la
simulacin de parques elicos arroje resultados ms precisos y acordes con la operacin
y mantenimiento ejemplificado con el Parque Elico Jeprachi.

An cundo en el documento se verific que las generaciones de los parques elicos se
veran beneficiados con el aumento de precios en pocas de sequa y que en operacin
integrada, a pesar que se gana la energa firme se pierde la oportunidad de tener mayor
energa cuando los precios estn ms altos. Por lo anterior la evaluacin financiera
deber realizarse a nivel mensual para analizar este comportamiento.

Se deber profundizar en el comportamiento de la operacin integrada en pocas con
influencia del Evento ENSO e incluir dentro de los parmetros de la RNO criterios que
permitan incluir condiciones ms ajustadas al mercado energtico colombiano.

A manera de ejercicio se debern hacer anlisis a nivel horario dentro del SIN
considerando los picos de la energa elica y diferentes capacidades instaladas que
impacten a la matriz para verificar comportamiento asociados que permitan establecer o
recomendar estructuras de despacho, tanto de Centrales Hidroelctricas como Parques
Elicos.

Ingeniero Civil

97

10 BIBLIOGRAFA Y FUENTES DE INFORMACIN

Atlas de Viento y Energa Elica de Colombia. Unidad de Planeacin Minero Energtica
UPME, Instituto de Hidrologa, Meteorologa y Estudios Ambientales de Colombia
IDEAM. Ministerio de Minas y Energa. Bogot, 2006.

Autogeneracin en Colombia. 4. Encuentro de uso Racional de Energa UPME ANDI.
Bogot, Diciembre 2008.

Informe Anual de Viento 2009 - AWEA. http://www.awea.org/publications/reports/AWEA-
Annual-Wind-Report-2009.pdf. Consultada el 2 de Octubre de 2010.

BELANGER, Camille; GAGNON, Luc. Adding wind energy to hydropower. Energy. Les
Consultants LBCD inc. Rene-Levesque W. Energy Policy 30. Canad 2002.

BERRY, Trent; JACCARD, Mark. The renewable portfolio standard: design considerations
and an implementation survey. Energy Policy 29 (2001) 263-277. Julio del 2000.

BLANCO, Mara Isabel. The economics of wind energy. Department of Economics,
University of Alcala. Madrid, Spain. 2008.

CHAPRA Steven, RAYMOND P. Canale. Mtodos Numricos Para Ingenieros. Mxico:
McGraw-Hill, 2003.

CREG Comisin de Regulacin de Energa y Gas. Proyecto de resolucin 092 de julio
de 2011 - Por la cual se define la metodologa para determinar la energa firme para
Plantas Elicas. Ministerio de Minas y Energa. Bogot, 2011.

DENAULT, Michel; DUPUIS, Debbie and COUTURE, Sbastien. Complementarity of
hydro and wind power: Improving the risk profile of energy inflows. HEC Montreal
(Management Sciences) and GERAD, Risk Management, National Bank of Canada.
Energy Policy. Canada, 2009.

MEJA, F., MESA, O., POVEDA, G., VLEZ, J., HOYOS, C., MANTILLA R., BARCO, J.,
CUARTAS, A., MONTOYA M., y BOTERO, B. Distribucin Espacial y Ciclos Anual y
Semianual de la Precipitacin en Colombia. DYNA, No. 127, 7-24, Agosto, 1999.

Foro Contexto y Normalizacin de la energa elica y solar en Colombia, UPME
ICONTEC. Bogot, Noviembre 2008.

FRANCO, Juan Felipe; DYNER, Isaac. Evaluacin de un Portafolio de Generacin de
Electricidad, utilizando Dinmica de Sistemas. Energtica 32. Diciembre 2004.

INTERNATIONAL STANDARD IEC 61400-1, Third edition, Wind turbines Part 1: Design
requirements, 2005.


Ingeniero Civil

98

ISAAC, Idi Amin. Energa Elica en Colombia: Implicaciones para la Transmisin y la
Operacin. Foro de Energa Elica. Medelln. Julio 2009.

ISAAC, Idi Amin; BIECHL, Helmuth, GONZALEZ Jorge W. La Energa Elica en Alemania:
Experiencias a tener en cuenta para el caso Colombiano. Revista Investigaciones
Aplicadas. Universidad Pontificia Bolivariana. Medelln 2008.

ISAGEN. Fuentes no Convencionales de Generacin de Electricidad. ISAGEN S.A. E.S.P.
2005.

JOSELIN, Herbert G.M., INIYAN S., A review of wind energy technologies. Agosto de
2005.

LIU, Timothy; TANG, Wenqing and XIE, Xiaosu. Wind Power Distribution Over the Ocean.
Geophysical Research Letters. Vol. 3. Julio de 2008.

Ministerio de Ambiente, Vivienda y Desarrollo Territorial - MAVDT; Departamento Nacional
de Planeacin. Lineamientos de Poltica de Cambio Climtico Resumen Ejecutivo. Junio
2002.

MONTGOMERY Douglas, RUNGER George. Probabilidad y Estadstica Aplicadas a la
Ingeniera. Mxico: McGraw-Hill, 1997.

ODILON, A. C: do Amarante; SCHULTZ, Daro J.; BITTENCOURT, Rogiro M. y ROCHA,
Nelson A. Wind/Hydro complementary seasonal regimes in Brazil. 2001 DEWI Magazine
19.

POVEDA Germn, Gil Marta Mara y Quiceno Natalia. El Ciclo Anual de la Hidrologa de
Colombia en Relacin con el ENSO y la NAO. Boletn del Instituto Francs de Estudios
Andinos. 1998.

POVEDA, Germn. La Hidroclimatologa de Colombia: Una sntesis desde la Escala Inter-
decadal hasta la Escala Diurna. Revista de la Academia Colombiana de Ciencias, Vol. 28
(107), 201-222. , 2004.

PRYOR, S.C.; BARTHELMIE, R.J. Climate change impacts on wind energy: A review.
Atmospheric Science Program, College of Arts and Sciences, Bloomington 2009.

UPME. Las Energas Alternativas y la Cogeneracin y Autogeneracin en Colombia. 4.
Encuentro de Uso Racional de Energa UPME ANDI. Bogot, Diciembre 2008.

VERGARA, Walter, DEEB, Alejandro, TOBA, Natsuko, CRAMTON, Peter, LEINO, Irene.
Wind Energy in Colombia A Framework for Market Entry. Banco Mundial, 2010.

XM. Informe Especial Resolucin CREG 135/97 - Noviembre 2009. Disponible en
http://www.xm.com.co/Pages/HistoricodeInformesEspeciales.aspx. Consultada el 25 de
febrero de 2010.

Ingeniero Civil

99

UPME Plan de Expansin de Referencia 2009-2023. www.upme.gov.co

Cambio Climtico 2007 Informe de Sntesis. Grupo Intergubernamental de Expertos
sobre el Cambio Climtico, 2008.

ENSO Cycle: Recent Evolution, Current Status and Predictions. Update prepared by
Climate Prediction Center / NCEP. www.cpc.ncep.noaa.gov/products/expert_assessment/.
Consultada permanentemente.

REN21., Renewables 2011 - Global Status Report. REN21, Secretariat Paris, 2011.
Informe descargado en lnea, sitio WEB: www.ren21.net. Consultado el 24 de julio de
2011.

RESOLUCIN No. 180947 DE 2010 del Ministerio de Minas y Energa. Diario Oficial No.
47.734 de 8 de junio de 2010.

Revista Empresas Pblicas de Medelln. Volumen 15. Enero Junio de 2004.

Special Report on Renewable Energy Sources and Climate Change Mitigation Wind
Energy, Final Release. IPCC, 2007.

Sitio WEB HydroQuebec consultado el 20 de julio de 2011:
http://www.hydroquebec.com/sustainable-development/documentation/ges.html.

Sitio WEB EMD International A/S. http://www.emd.dk/WindPRO/Frontpage, consultado el
28 de julio de 2011.

The Wind Atlas Analysis and Application Program, desarrollado por National Laboratory
for Sustainable Energy de Dinamarca, http://www.risoe.dtu.dk/da/WAsP.aspx?sc_lang=en
consultada en julio de 2011.

TREMBLAY, A., SCHETAGNE, R. 2006. The Relationship Between Water Quality and
Greenhouse Gas Emissions in Reservoirs. International Journal on Hydropower & Dams,
13(1):103-107.

TREMBLAY, A., L. VARFALVY, C. ROEHM and M. GARNEAU. 2004. The issue of
greenhouse gases from hydroelectric reservoirs: from boreal to tropical regions. In:
Proceedings of the United Nations Symposium on Hydropower and Sustainable
Development. Beijing, China, October, 2004.

Diseño de Energia Eolica en La Guajira

Cargado por

Copyright:

Formatos disponibles

Diseño de Energia Eolica en La Guajira

Cargado por

Información del documento

Título original

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

Diseño de Energia Eolica en La Guajira

Cargado por

Copyright:

Formatos disponibles

ANLISIS DE GENERACIN COMPLEMENTARIA ENTRE

ENERGA HIDRULICA Y ELICA

es el resultado de la DFT para

También podría gustarte