Calculo de Potencia de Aerogenerador PDF
Calculo de Potencia de Aerogenerador PDF
Calculo de Potencia de Aerogenerador PDF
ESCUELA DE INGENIERIA
Profesor Supervisor:
RODRIGO ESCOBAR MORAGAS
cuando lo necesité.
i
AGRADECIMIENTO
Quiero agradecer a mi familia por apoyarme en todo momento y brindarme todas las
condiciones necesarias para poder desarrollarme como ingeniero y persona.
A mi profesor guía, Rodrigo Escobar por ser un cable a tierra para poder finalizar de
buena manera esta tesis.
A las empresas AngloAmerican S.A y ABB S.A que con su visión emprendedora
invirtieron en parte de este proyecto lo que logró que este estudio pudiese ser llevado a
cabo.
A la gente que me apoyó con sus conocimientos técnicos y profesionales para poder
obtener las competencias necesarias para completar este desafío a Mario González,
Leonardo Ramírez, Adolfo Escobar, Nibaldo Candia, Jaime Zamora, al departamento de
Mecánica y Eléctrica por atender mis consultas y a todos los auxiliares y secretarias del
Departamento de Mecánica por su buena disposición.
ii
ÍNDICE GENERAL:
AGRADECIMIENTO .......................................................................................................ii
ABSTRACT .......................................................................................................................ii
1 INTRODUCCIÓN ......................................................................................................... 1
1.3 Tecnologías eólicas: Aerogeneradores de eje vertical versus eje horizontal ........... 5
2 OBJETIVOS ................................................................................................................... 7
4.1.- Métodos utilizados con los datos de viento obtenidos en terreno ....................... 43
iii
5 Resultados y Discusión ................................................................................................. 79
5.1.- Análisis de las curvas de viento obtenidas a partir de los datos en terreno ......... 79
iv
ÍNDICE DE TABLAS
Tabla 6: Eficiencias de turbina Giga Star para estimación del VAWT prototipo. .......... 39
Tabla 10: Resultado de los filtros aplicado para un período de configuración. ............... 81
Tabla 11: Tabla resumen de la cantidad de errores para el período de medición. ........... 82
Tabla 13: Tabla de resultados con las características mensuales más importantes. ........ 89
Tabla 14: Comparación de los factores del coeficiente de cizalle ( en distintos sitios 90
Tabla 15: Resultados de coeficientes de forma y escala según meses del estudio .......... 94
Tabla 18: Coeficientes que minimizan la MC de la estimación de velocidad por PLP. 100
Tabla 19: Estimación desde cada anemómetro a los otros. ........................................... 102
v
Tabla 22: Tabla resumen de la energía exportada e importada...................................... 120
Tabla 23: Multiplicaciones posibles con la combinación de las cajas epicíclicas. ........ 139
vi
ÍNDICE DE FIGURAS
Figura 2: Frecuencia de velocidad del viento en Chagres con Curva Weibull. ............... 10
Figura 13: Esquema General del drive controlador provisto por ABB. .......................... 25
Figura 17: Comportamiento de la frecuencia de los vientos en Chagres para un día. ..... 84
Figura 19: Resultado de las turbulencias para las velocidades del perfil......................... 87
Figura 20: Proyección de las velocidades del viento promedio según LLP y PLP. ........ 92
Figura 21: Perfil de velocidad diario para el mes de Agosto del 2011. ........................... 93
Figura 22: Comparación de curvas de probabilidad Weibull para distintos terrenos. ..... 95
vii
Figura 23: Velocidad del viento medida por las distintas sondas. ................................... 96
Figura 27: Curva revoluciones vs tiempo cada 1 minuto para el 16/11/2011................ 107
Figura 28: Rpm en comparada con las velocidades diarias en 10 minutos.................... 107
Figura 29: Frecuencia de las Revoluciones por minuto del aerogenerador. .................. 108
Figura 30: Velocidad de rotor en aerogenerador respecto del viento sin carga. ............ 109
Figura 31: Tip Speed Ratio vs velocidad del viento en vacío. ....................................... 110
Figura 32: Primera curva de torque programada en el software DriveAP. .................... 111
Figura 33: Interfaz del programa adaptativo con el programa de la curva inicial. ........ 112
Figura 34: Análisis diario de parámetros más importantes cada 10 minutos................. 114
Figura 37: Valor del Coeficiente de potencia según el TSR. ......................................... 119
Figura 40: Representación de los desplazamientos sufridos por el brazo inferior......... 133
Figura 41: Representación del factor de seguridad para el análisis del brazo inferior .. 134
viii
Figura 46: Interfaz gráfica del programa DriveWindow. .............................................. 144
Figura 48: Velocidad del viento promediada cada diez minutos ................................... 145
ix
LISTA DE ACRÓNIMOS
x
RESUMEN
A continuación se analizó la velocidad angular del rotor. Con esto se pudo ver los
límites de las revoluciones para calcular la relación de transmisión según el motor que se
encuentre disponible. Luego se clasificó el aerogenerador según los principios físicos
que ocurren alrededor de sus álabes, sustentación o arrastre, lo que llevará a estimar su
eficiencia.
i
ABSTRACT
During the last three years, Chile has been developing a huge market related to wind
energy resources. Many sites have already been tested and the strong wind potential of
our country is being projected in the construction of new wind parks.
The VAWT were introduced in the market of small wind generators and were driven by
the fact that they are easily built, they have low decibel emission and they work
regardless of the wind direction. Lately, in the search for more efficiency and scale of
economies, bigger prototypes of this technology have been developed. Due to this, it is
necessary to look for a standardized methodology in order to obtain the power curve of
the VAWT. This methodology should be based on the advances in the characterization
made in the HAWT.
The first part of this methodology consists in studying the specific place where the
turbine is going to be located. All the aerodynamic factors should be measured with the
geographic aspects, which enable a complete characterization of the wind turbine.
Furthermore, it will be necessary to estimate a representative value of the wind speed
profile, the density, and the swept area, in order to be able to calculate the power of the
wind current.
Afterwards, the behavior of the angular velocity of the rotor should be evaluated. This
will allow the calculation of the transmission according to the available motor. With this
result, it will be possible to classify our wind turbine according to the physics principles
which are involved between its blades.
Finally, all the electric parameters delivered by the turbine will be measured according
to the wind speed. The curves of power versus the wind speed and the power coefficient
versus the TSR will be developed. After this is done, there will be a comparison between
the energy production of the power curve and the energy measured by the sensors.
All of these procedures will be applied in a VAWT turbine, to show the magnitude and
the shape of the results of this method.
ii
1
1 INTRODUCCIÓN
Desde el descubrimiento de la teoría del pico de Hubbert1, los científicos del planeta han
estado intentando desarrollar nuevas fuentes energéticas que puedan garantizar de forma
económica la posible escasez de fuentes no renovables. Esto, unido a la teoría del
calentamiento global, hace necesaria la búsqueda de fuentes alternativas de energía y de
diversificación de las ya existentes.
1
http://en.wikipedia.org/wiki/Peak_oil
2
http://www.bp.com/sectiongenericarticle800.do?categoryId=9037191&contentId=7068648
2
económica de las turbinas de eje horizontal para algunos nichos, éstas son las turbinas de
eje vertical, las cuales actualmente se encuentran mayormente desarrolladas en bajas
potencias y para ser instaladas en lugares urbanos.
Chile, por su diferenciada geografía de larga costa, tiene lugares con alto potencial de
velocidad y densidad de viento, lo que lo hace un país atractivo para invertir en energías
renovables, especialmente en energía eólica. Es por esto que esta tesis se ha enfocado en
aportar a un posible desarrollo de torres eólicas que podrían ser construidas en Chile,
ahorrando los enormes costos que tiene el traslado de las torres, debido a su gran
volumen y delgada geometría.
En Chile existe una extensa costa donde es muy probable encontrar buenos vientos para
la generación eólica, lo que se puede comprobar con el creciente número de estaciones
de medición eólica que se han instalado durante los últimos años. Las estaciones
meteorológicas para evaluación de proyectos eólicos son fácilmente identificables con
respecto al resto debido a la altura del mástil que sostiene los anemómetros de medición.
Los mejores vientos se encuentran en la zona sur del país, en la Región de Magallanes,
sin embargo en esa región no se cuenta con grandes consumos eléctricos. Este sistema
interconectado solo abastece al 0,54% de la energía consumida en Chile, por lo que estos
vientos no pueden ser aprovechados en los grandes consumos de la minería y las
ciudades ubicadas mayormente al norte del país, por lo que sólo el parque eólico más
3
CNE; CDEC-SIC; CDEC-SING (2011). Generación Bruta Sistemas SING (1999-2011)
SIC (1996-2011) Recuperado el 10 de Noviembre del 2011, de www.cne.cl.
4
pequeño está ubicado en esta región. El resto de los parques se encuentran en la costa de
la Cuarta Región, la cual cuenta con velocidades de viento considerables y está ubicada
más cerca de consumos mayores, como es la ciudad de Santiago.
Para poder estimar de manera acertada los vientos en la costa de Chile, la Universidad
de Chile ha publicado un mapa eólico, basado en distintos datos de las zonas de
medición y usando una extrapolación para resultados entre estos puntos4. Si bien esta
aproximación puede resultar irreal en la práctica, pues todo proyecto eólico necesita de
mediciones in situ de por lo menos un año para reducir la incertidumbre del proyecto, da
un primer indicador global de dónde se encontrarían los mejores vientos del país.
En figura 1 se aprecia una vista general del mapa eólico desarrollado por la Universidad
de Chile, enfocado en los vientos de la costa, que muestra cómo los vientos al extremo
sur son de mayor magnitud que los vientos costeros del litoral central y norte.
4
CERE-UMAG
5
2. La velocidad tangencial del extremo de los álabes es menor, por lo que genera baja
intensidad de ruido, menos vibraciones, levanta poco polvo y tiene menor impacto en la
fauna silvestre.
3. Son más fáciles de instalar, requieren de menor mantenimiento y las cajas reductoras
y otros equipos eléctricos pueden ser colocados en la base de la torre.
1. Hasta el momento se han logrado menores coeficientes de potencia que los HAWT
(Cp=0.45)
2. Debido al método constructivo, es más fácil que los HAWT puedan llegar a mayor
altura, que es donde se encuentran las mayores velocidades de viento.
3. En general, necesitan más material por metro cuadrado de superficie cubierta, lo que
aumenta los costos de la turbina.6
5
http://energy.gov/oe/technology-development/smart-grid.
6
Ropatec AG (2003). HAWT versus VAWT Small VAWTs find a clear niche,
Renewable Energy Focus, 1471 0846/03, 44-46.
6
La tecnología más utilizada para extraer energía mecánica o eléctrica del viento son las
turbinas de eje horizontal. Estas turbinas se caracterizan por generar una cierta potencia
con un determinado costo de inversión, el cual debe estar siendo constantemente
comparado con los de otras tecnologías de energía, tanto eólicas como de otras fuentes.
Para realizar la comparación entre HAWT y VAWT es que en esta tesis se desarrollará
un método para obtener la potencia y la eficiencia de una turbina de eje vertical.
La necesidad de desarrollar un método para los VAWT viene del crecimiento paulatino
que han tenido los aerogeneradores verticales de pequeña escala para usos urbanos.
Estos pequeños aerogeneradores han ido aumentando su tamaño para buscar mayores
eficiencias y economías de escala, llegando a haber algunas turbinas comerciales de gran
escala para esta tecnología (30 kW)7 diseñadas para ser utilizados dentro o en las afueras
de una ciudad. Este crecimiento de envergadura provoca un incremento en las
inversiones en esta tecnología, aumentando el interés de crear un procedimiento
estandarizado para la caracterización de estas turbinas.
El método desarrollado para obtener la curva de potencia de la turbina debe ser acorde a
las normas internacionales de medición de viento8 y por ende debe respetar los
intervalos de muestreo, los análisis recomendados para esos desarrollos y otras
características fundamentales de esas normas; de esta manera se tendrá más exactitud a
la hora de evaluar un proyecto singular o un parque eólico y se podrá evaluar la turbina
en los terrenos ya censados mediante las normas de medición existentes.
7
www.ropatec.com
8
No existe una norma específica para censar el viento pero si normas en donde se recomiendan ciertos
procedimientos tales como la IEC 614001.
7
2 OBJETIVOS
En el siguiente capítulo se muestra una descripción del terreno utilizado para hacer el
estudio, tomando en cuenta las principales características de la materia prima disponible
que es el viento; luego se describe la turbina que será analizada por la metodología; y se
detallan todos los sistemas e instrumentos necesarios para la rápida y correcta obtención
de la curva de potencia, que es el objetivo principal de esta tesis.
Luego se describirán los componentes y las características del prototipo utilizado para
esta tesis, se dedicará especial atención a aquellas partes que afecten directamente a los
cálculos de la potencia, tales como: superficies alares, pesos, geometría de estructuras,
componentes de transmisión y roces producidos en la interacción de las distintas piezas,
entre otros.
- Esta variación de velocidad que ocurre entre la mañana y la tarde ayuda a las labores
de mantenimiento o de adaptación que se deben hacer para montar o hacer cambios de
configuraciones en el aerogenerador, permitiendo un fácil acceso para realizar trabajos.
La tabla 2 y la tabla 3 muestran un resumen con todas las principales características del
terreno durante el período en que se hicieron las pruebas. Es importante tener en cuenta
que durante el proyecto fueron aumentando los sistemas de monitoreo y, por lo tanto, en
los primeros meses se cuenta con menor información que al final.
Este análisis inicial tomó en cuenta varios factores que caracterizan al viento. Es más
completo que los análisis meteorológicos que se hacen en la mayoría de los terrenos,
pero para validar la turbina se necesitará ser aún más específicos en cuanto a las
particularidades implícitas en el sitio para poder validar la turbina, considerando la
mayoría de las condiciones que los terrenos puedan tener.
2500.0
Frecuencia de Velocidad de viento al año
Numero de datos [1 min]
2000.0
1500.0
Frec6m
1000.0 Frec12m
Frec18m
500.0
W12m
0.0
10.5
11.5
12.5
13.5
14.5
15.5
6.5
0.5
1.5
2.5
3.5
4.5
5.5
7.5
8.5
9.5
4.00
3.50 de altura
3.00
2.50
2.00
[m/s]
1.50
PROMEDIO A 12M
1.00
0.50
0.00
Es importante destacar que las normas utilizadas para obtener las curvas de potencia de
los aerogeneradores de eje horizontal indican, que para obtener la curva de potencia de
manera correcta es de importancia que el terreno en el cual se haga la medición sea del
tipo planicie, sin embargo esta aseveración es contrastada por los recientes estudios, en
donde se ha demostrado que se pueden llegar a excelentes resultados para las curvas de
potencia aplicando los mismos procedimientos obtenidos de la norma para terrenos con
pendientes entre 3º y 12º en zonas cercanas a las mediciones9. Estas afirmaciones son de
gran utilidad en esta situación, pues las laderas presentes en el terreno no desviarán de
mayor manera los resultados obtenidos mediante el procedimiento.
El aerogenerador construido es de eje vertical, lo que quiere decir que el eje que
transmite la potencia está ubicado en dirección perpendicular al flujo del viento y, por
tanto, los álabes captadores de la energía cinética de viento también lo están. La empresa
a cargo de este proyecto comenzó la construcción del aerogenerador pensando en
obtener una potencia de 100kW. Por lo tanto las piezas que lo forma están diseñadas
para esta potencia.
9
U. Bunse, H. Mellinghoff, O. Haack (2006). Uncertainty of annual energy production
for a specific turbine model based on a set of IEC 61400-12 measurements. Recuperado
de German Wind Energy Insitute GmbH (DEWI): http://www.dewi.de
10
http://www.wired.com/science/discoveries/news/2006/10/71908 Recuperado el 20 de Diciembre
13
Primero se muestra la figura 5 que marca los componentes que serán descritos a
continuación:
La torre es una estructura compuesta en su mayor parte por perfiles cilíndricos soldados.
Fue pensada tanto para su fácil transporte, como para que resistiera los esfuerzos
generados por la velocidad del viento. La torre puede ser acoplada fácilmente a un
14
camión y la base no excede los 3 m de longitud para que pueda ser transportada a través
de las carreteras de Chile.
Tiene un perfil cuadrado que le da estabilidad y se ubica sobre una base inferior en la
que va apoyada. La idea principal de la torre es que sea capaz de resistir el peso del
rotor, que es de 7 toneladas aproximadamente, dependiendo de la configuración de rotor
que se utilice, además de soportar las cargas variables ejercidas por el viento. Para esto
existe una base cuadrada de 12 m que va acoplada por debajo de la torre, para así poder
resistir de mejor manera los torques que pudiesen producir un volcamiento.
Como se vio en la Introducción, una de las ventajas de los equipos de eje vertical es que
ellos pueden aprovechar la energía del viento que viene en todas las direcciones. Es por
este motivo que la torre es completamente simétrica y será instalada en lugares que
presenten grandes velocidades de viento sin importar su dirección; una alta velocidad del
viento acelerará la toma de mediciones y con esto podrá obtener una mejor y más pronta
caracterización, lo que lleva a tomar mejores conclusiones sobre su funcionamiento.
El rotor es el componente que está fijo a la torre, pero que permite el libre giro de los
patagios alrededor del eje cardánico central. Tiene dos estructuras de anclaje firme; una
es la unión con la torre mediante soldadura y la segunda es la que mantiene suspendidos
a los brazos inferior y superior mediante unos pasadores. Todo esto va montado sobre un
rodamiento que debe estar diseñado para resistir el peso del rotor completo.
La conexión entre el rotor y los patagios esta provista por los brazos, los cuales se
descomponen en brazo inferior y brazo superior, y están pensados para soportar dos
patagios de 800mm, o bien uno de 1600mm. Estos brazos son telescópicos por lo que
tienen la posibilidad de estirarse u encogerse, barriendo de esa manera distintas áreas y
pudiendo variar su rango de operación estimado, el cual está entre los 10.2 y los 12.2 m
de diámetro.
A continuación la figura 7 muestra una vista de los brazos superiores unidos con el rotor
antes de ser montados.
En la tabla 4 se muestran los principales pesos de los brazos, los patagios y el rotor
vistos desde sus valores unitarios hasta los pesos correspondientes a la configuración
final analizada para esta tesis.
Una de las grandes características de los aerogeneradores de eje vertical es que el eje de
salida del rotor es vertical, por lo tanto se puede crear un sistema de transmisión,
mediante cardanes, flanges y cajas multiplicadoras, que lleve la potencia mecánica del
rotor a la base de la torre, haciendo que el generador eléctrico y el variador de
frecuencia se ubiquen en un lugar cercano a la base de la torre, disminuyendo el peso en
la parte superior de la torre y, por ende, parte considerable de los esfuerzos que ello
genera.
y al aumentar el tamaño del brazo la velocidad angular debe disminuir para mantener las
proporciones de la velocidad tangencial.
Para dimensionar qué cajas pueden ser parte de la estructura del aerogenerador, se
contempló tanto la potencia de la máquina a la cual servían, como el ancho del eje con el
cual transmiten la potencia.
iv) Generador
El motor eléctrico fue provisto durante 6 meses por la empresa ABB, como aporte al
desarrollo del proyecto y fue el primer motor de imanes permanentes en Chile en tener
una capacidad nominal de 100 kW. Luego el generador fue comprado por Poliex para
seguir con las investigaciones.
Torque: 2010 Nm
Frecuencia: 50 Hz
b) Dimensiones y Características
Carcasa: largo: 1,408m; ancho: 0,59m; alto: 0,852m; diámetro brida: 0,66m.
v) Sistema de Control
11
Especificaciones y dimensiones obtenidas de ABB.
24
El drive fue seleccionado mediante el programa DriveSize provisto por ABB, el cual
recomienda el variador de frecuencia que mejor se comporte según las distintas
características del motor y la carga que recibirá. En este caso la carga se tomó como
similar a la que produce una bomba hidráulica, es decir carga cuadrática.
a) Datos de Placa
Potencia nominal: 90 kW
b) Dimensiones y Características
Figura 13: Esquema General del drive controlador provisto por ABB.12
El generador y drive serán instalados dentro de la torre, cercanos al nivel de suelo, a
diferencia de los aerogeneradores de eje horizontal, lo que disminuirá las dificultades
asociadas a la instalación y mantención de los equipos. En un primer piso se instalará el
drive, luego en un segundo piso el generador eléctrico. Hay que notar que ambos se
encuentran separados a una distancia prudente permitiendo la ventilación de los equipos.
La estructura mostrada en la figura será dispuesta dentro de una cabina que los aislará
del exterior.
Es importante destacar que el drive cuenta con un sistema de medición el cual es capaz
de medir las variables más importantes dentro de la generación eléctrica, como son los
voltajes, las corrientes, el torque hecho sobre el eje, las temperaturas, el voltaje del bus
de corriente continua, y otros parámetros que se mencionarán más adelante.
Para la medición de las variables correctas, la IEC ha definido que hay que medir en la
salida del sistema de corriente eléctrica antes del transformador13. En este caso inyectó a
380 V L-L por lo que no habrá un transformador elevador antes de conectarnos a la red
de AngloAmerican.
12
Planos entregado por ABB S.A.
13
International Electrotechnical Commission (2005). Power performance measurements of electricity
producing wind turbines. (61400-12-1 Ed1). Holland: Author.
26
Para poder medir y registrar todas las variables que influyen en la generación eléctrica
del aerogenerador, se debió instalar una serie de instrumentos que sean capaces de medir
y registrar de manera confiable los parámetros que ayudarán a obtener la curva de
potencia del aerogenerador. En esta sección se incluyen recomendaciones importantes
para que la instalación del sistema funcione correctamente. Estas sugerencias son parte
de la experiencia de la instalación del sistema en la turbina descrita en la sección
anterior.
aerogenerador
Para esta situación la utilización de 3 anemómetros fue muy importante para medir bien
el elevado factor de coeficiente de cizalle del sitio, este coeficiente indica la variación
que presenta la velocidad del viento respecto de la altura. En este caso el coeficiente es
producido debido a los distintos obstáculos que existen en la geografía cercana al
emplazamiento, la cual hace variar considerablemente la velocidad del viento con la
altura, en comparación con los terrenos costeros generalmente utilizados en Chile para la
generación eoloeléctrica.
27
Los anemómetros son de copas que es el tipo de anemómetro recomendado por la IEC, y
fueron comprados nuevos para que estuvieran calibrados. Al estar todos a poca
diferencia de altura, da la posibilidad de contrastar las mediciones de unos con otros para
ver si se ha mantenido la calibración durante su vida útil, lo cual es una gran ventaja.
El mástil que sostiene los anemómetros para el caso de las turbinas de eje vertical se
debe encontrar a una distancia de 1.5 D, en donde para este caso se definió √ ,
donde A es el área que representa el plano vertical que atraviesa al rotor por su
ecuador.14 En este caso la distancia al mástil debería ser 10.2 m.
14
International Electrotechnical Commission (2005). Power performance measurements
of electricity producing wind turbines. (61400-12-1 Ed1). Holanda: Autor.
28
para así comprobar la relación de multiplicación de las cajas, ver rápidamente cuándo
una de estas cajas se frene o se cizallen los pernos de la unión al cardán. El prototipo a
analizar contempla dos cajas multiplicadoras, por ende se usarán dos sensores de rpm.
Con ambos datos más las mediciones de las revoluciones hechas por el motor eléctrico
se pueden obtener datos de todas las velocidades a las que está sometido el cardán.
Los dos sensores instalados dan la capacidad de prever cualquier pérdida de transmisión
que ocurra en el eje cardán, y ayuda a contrastar la veracidad de los datos. El sensor de
medición de las revoluciones es muy simple y a la vez muy preciso; consiste de un
sensor magnético que envía ondas de pulso cuando el giro hace pasar un metal a
milímetros del sensor.
Estos sensores de medición deben almacenar cada registro con su hora de ocurrencia
tomada desde el reloj central del computador. Con esto se puede obtener la curva de
revoluciones por minuto sobre la velocidad del viento que tiene cada una de las
configuraciones del aerogenerador. La curva de revoluciones es la primera fase antes de
obtener la curva de potencia de éste, y para ser obtenida, la base de datos de la velocidad
angular debe ser unida a la base de datos de la velocidad del viento.
Se monitoreó la temperatura de las cajas para verificar que no estén ocurriendo pérdidas
excesivas de calor por roce en los contactos que ellas tienen.
29
Esta categoría considera: Voltajes (de salida del motor y de salida hacia la red eléctrica),
amperajes, potencias (activa, reactiva y aparente), frecuencia, armónicas, energía
generada, torque eléctrico, corriente en el bus DC y factores de potencia, por esto se
cuenta con dos sistemas independientes.
El primer sistema es el provisto por el drive, el cual tiene su propio analizador de redes,
de alta velocidad de registro para almacenar variables, pues puede tomar datos cada
milésimas de segundo, pero tiene poca capacidad para monitorear distintas variables,
pues solamente acepta hasta 6 variables distintas al mismo tiempo. Es fundamental que
el drive pueda accionar al motor con un control de torque programable (DTC), para así
poder optimizar la generación eléctrica.
Todos los instrumentos instalados en la localidad van conectados según sus distintos
sistemas de transmisión de información a un computador central de bajo consumo, que
almacena los datos y es capaz de ser accedido remotamente para poder obtenerlos desde
cualquier punto en que se tenga conexión a internet.
Para el correcto control del aerogenerador se instalaron cuatro cámaras en los puntos
más importantes de la estructura. Estas cámaras pueden ser vistas a través del mismo
acceso remoto que se utiliza para rescatar la información, y entrega a distancia datos
sobre el estado de los distintos componentes del aerogenerador. Las cámaras están
ubicadas en las siguientes posiciones: exterior con vista al aerogenerador completo, en el
rotor, en la segunda caja multiplicadora y en la entrada del generador eléctrico.
31
Los sensores son muy sensibles a la corriente eléctrica que los alimenta y a la corriente
que alimenta al computador que recibe la información de los datos crudos, por lo que se
recomienda crear un sistema eléctrico confiable y que genere corriente continua y
corriente alterna de buena calidad para que así los registros de los sensores no se vean
afectados por distorsiones. Si se conecta el computador y los sensores a la misma red
que está conectada al drive, lo más probable es que muchos datos se pierdan.
Para que todo el sistema sea accesible desde internet, se debe tener una conexión
confiable y de buena velocidad, dependiendo de la cantidad de cámaras que se ubique.
Si se quieren transmitir los datos en tiempo real, se recomienda utilizar servicio de
internet a través de telefonía de la compañía que tenga mejor señal en la zona, que sea
3G o superior en velocidad y colocar la antena lo más alejada posible de la radiación del
drive.
a) Router 3G D-Link.
32
15
La precisión de cada instrumento debe verse en los documentos IEC y Measnet asociados a esta tesis.
33
teórica que tendrá la turbina eólica. Diseños menores llevarán a encontrar posibles fallas
en el sistema, y diseños mayores llevarán a pérdidas de eficiencia en roce y peso de
sobra.
Para definir el límite teórico primero debió introducir el significado del coeficiente de
potencia ( ) el que indica la fracción de potencia que el aerogenerador es capaz de
extraer del viento y está directamente relacionado con la eficiencia que presenta la
turbina eólica. Con el resultado de este coeficiente se pudo ver cuál será la potencia
disponible en la turbina a partir de una cierta potencia en la corriente de viento. La
relación está dada por la siguiente ecuación:
[1]
Como es imposible que la turbina eólica extraiga toda la potencia del viento definió en
primer lugar . Luego, Albert Betz consideró que:
- El flujo de aire es incompresible debido a que la velocidad del viento para este tipo de
proceso ocurre a velocidades Match inferiores a 0.3 ( ).
- A través de la turbina eólica ocurre una expansión del aire que ingresa al aerogenerador
disminuyendo así su velocidad a medida que el aire la atraviesa.
34
Sin embargo, la velocidad final no puede ser nula debido a que debe haber circulación
constante del aire para que este se mantenga ingresando al sistema. Por lo que existe un
cociente entre la velocidad de entrada y la velocidad de salida que maximiza la potencia
entregada por el flujo de aire.
En la figura 15 mostramos un esquema de la expansión del flujo que cruza a través del
aerogenerador.
Luego utilizando la diferencia de energía entre los fluidos de entrada y salida se obtuvo
una ecuación de la potencia extraída [2]. Esta expresión fue igualada con la potencia
obtenida del flujo de aire considerando la Segunda Ley de Newton obteniendo una
fuerza que al ser multiplicada por la velocidad entrega la potencia en función de la
velocidad de entrada y de salida [3]. Estas ecuaciones igualadas quedan de la siguiente
forma:
( [2]
35
( [3]
ecuación [1] obtiene una expresión de la potencia útil en base a la velocidad de entrada y
salida, pero realizando el siguiente cambio de variable , queda solamente la
potencia en función de la velocidad de entrada.
( ( [4]
Esta expresión se deriva con respecto a y se iguala a cero para obtener su máximo,
con eso se obtiene la relación entre las velocidades de entrada y de salida, tal que se
produzca la mayor cantidad de potencia disponible teórica para la turbina.
( ) [5]
Luego de la demostración sobre la potencia teórica máxima para toda máquina que
convierta energía eólica en potencia mecánica, se procede a entrar más en detalle sobre
el coeficiente de potencia obtenido en aerogeneradores de eje vertical. Para esto se
diferencian los tipos de aerogeneradores verticales y se muestran las eficiencias
calculadas teórica y experimentalmente que se han podido alcanzar en la literatura de
este tipo de tecnología para los diferentes modelos de turbina.
Los aerogeneradores Savonius solamente funcionan mediante las fuerzas de arrastre, por
ende sólo son capaces de girar debido a la diferencia de presión entre ambos álabes, por
lo que sus perfiles alares están optimizados para que las diferencias de arrastre entre el
borde de ataque y el de salida produzcan la mayor diferencia de presión en todo el giro,
generando una velocidad angular y un torque en el eje principal. Estas máquinas son
fáciles de construir y son utilizadas principalmente como bombas de agua, pues tienen
partida automática y tienden a ser máquinas de elevado torque, pero son poco eficientes
alcanzando coeficientes de potencia experimentales cercanos de 0.0516
La tabla 5 muestra las potencias esperadas según los coeficientes de potencias obtenidos
de los aerogeneradores de la literatura según el área de superficie que debiese tener el
aerogenerador descrito en la sección anterior.
16
National Institute for Rural Engineering (2000). Output Characteristics of Savonius
Windmill Used as a Power Source Recuperado el 5 de Julio del 2011, de
http://www.nkk.affrc.go.jp/eng/
17
Kjellin, J. et al. (2010). Power coefficient measurement on a 12 kW straight bladed
vertical axis wind Turbine, Renewable Energy, 36 (2011) 3050-3053. Uppsala, Sweden.
37
Destacamos que se debiese tener una eficiencia del 65% para que el aerogenerador
construido por Poliex Ltda alcance los 100kW a una velocidad nominal de 12 m/s.
Para que la turbina sea comercialmente rentable, lo mejor es compararla con turbinas ya
existentes en el mercado, pero las turbinas de eje vertical se encuentran bien
desarrolladas para bajas potencias, aproximadamente entre los 1KW hasta los 10 KW en
donde son utilizadas en zonas urbanas, principalmente en la techumbre de casas y
edificios, por lo tanto resulta difícil compararlas con una equivalente en el mercado.
poder mejorar sus eficiencias económicas. En este nicho puede ver que los modelos
basados en el arrastre se quedan atrás en comparación con los modelos basados en
sustentación. Por ese motivo para áreas de superficie perpendicular superiores a los 35
existe en su mayoría turbinas del tipo Darrieus en el mercado. El área barrida por la
corriente es cercana a los 144 lo que obliga a comparar el prototipo puesto a prueba
en esta tesis con este tipo de tecnología, que es también, la más eficiente en cuanto al
factor de potencia.
18
Ropatec AG (2011). Features Ropatec´s small wind turbine vertical axis Recuperado el 10 de
Octubre del 2011, de http://www.ropatec.com
39
Tabla 6: Eficiencias de turbina Giga Star para estimación del VAWT prototipo.
Para el procedimiento tuvo en cuenta que los aerogeneradores que serán sometidos a este
análisis suelen ser prototipos en desarrollo, por lo que el procedimiento está estructurado
para potenciar la optimización de los equipos, tomando en cuenta que van cambiando
componentes para lograr mejores resultados. El procedimiento debe cerrar ciclos de
medición de manera rápida para minimizar: el tiempo de toma de decisiones, la
obtención de resultados y el tiempo que demora realizar los cambios en la estructura
física, que, en algunos casos, suelen ser cambios de larga duración de implementación y
otras veces más cortos. Estos tiempos se mueven aproximadamente entre tres días y tres
semanas dependiendo también del tamaño del aerogenerador que esté a prueba y de las
42
máquinas disponibles por los desarrolladores. Pero la medición también debe contemplar
ciclos de medición más largos para así tener análisis que contemplen una mayor muestra
de datos y sean más confiables. Debido a esto es imprescindible tener un registro lo más
completo y catalogable posible, para así tener la posibilidad de realizar cambios de
configuraciones rápidos sin perder información de cada configuración y cada cambio
realizado; con esto se podrá volver a analizar los resultados pasados, según los nuevos
resultados que se van obteniendo.
El análisis mensual ayuda a tomar en cuenta las variables generales del terreno y a
aquéllas que deben ser analizadas con un mayor número de datos. Es fundamental que
los análisis diarios hechos se vean reflejados en los cálculos mensuales como los
promedios de velocidad según la hora del día. Mientras más se vean reflejados estos
43
resultados diarios en los análisis mensuales, significa que el análisis diario muestra una
mejor aproximación como ciclo cerrado de datos.
Luego será necesario integrar todas estas variables medidas para obtener la curva de
potencia que describe al aerogenerador, es decir en base a este pensamiento se obtendrá
una curva de potencia por día que sería el menor ciclo completo que ha definido y esta
curva que podrá ser comparable con la de otros días para verificar el correcto
funcionamiento de los sistemas de medición para, finalmente, unir los distintos datos
diarios para formar finalmente una curva que este respaldada por una gran cantidad de
datos en el análisis por configuración.
Otra variable que afecta la potencia es la densidad del viento, que a su vez depende de la
temperatura y de la presión del sitio. Es por eso que se analizó el rango de temperaturas
y presiones que hay en el terreno para poder determinar qué tanto variará la densidad de
acuerdo a las distintas horas del día y épocas del año. De esta manera, hacia el final de la
sección, se obtiene una estimación para la potencia del viento de acuerdo a las
características obtenidas del terreno.
También serán mostrados efectos macro sobre las características del viento tales como el
comportamiento de la velocidad a lo largo del día y año, lo cual indicará en qué
momentos es preferible realizar cambios en las configuraciones del aerogenerador y en
qué momentos del día debe estar listo para que se tomen las pruebas de potencia; así
como los componentes de coeficiente de cizalle y de rugosidad del terreno, entre otros,
que diferencian a este terreno de otros que ya han sido estudiados.
El resultado de todas estas variables del viento entregará como resultado qué tan apto es
el terreno para lograr caracterizar a la turbina eólica, sin olvidar la importancia de otras
características logísticas y financieras que son estrictamente necesarias para mirar este
desarrollo como una posible inversión. Estas características son la cercanía a Santiago,
para la obtención de repuestos, realización de mantenciones, traslados más rápidos y
económicos, y la ubicación de una planta que tenga una red trifásica disponible para
poder integrar la energía producida a esta red.
Un análisis diario de las variables que se miden con los sensores instalados en el
aerogenerador es fundamental para tener una idea del comportamiento diario a través de
un periodo más extenso y para detectar posibles fallas en los sensores o en la estructura
45
del aerogenerador que retrase la recolección de datos fidedignos para ser procesados. En
el Anexo E se presenta un gráfico en dónde uno de los anemómetros presenta errores,
con esto se pretende mostrar cómo se distingue inmediatamente un tipo de error de
medición de un anemómetro cuando se le compara con los otros anemómetros que
tuvieron resultados dentro de lo esperado. Es posible percatarse de estos errores
mediante la realización y análisis de los gráficos de los parámetros monitoreados,
especialmente diseñados para este fin, y de los filtros de los datos que serán explicados
posteriormente.
En este caso es importante destacar que más que obtener las características del viento
necesarias para determinar la energía producida por el terreno al colocar una cierta
turbina ya caracterizada (procedimiento realizado por empresas que evalúan terrenos
para prospección eólica), se necesita sensar el terreno para caracterizar la turbina por lo
que es importante contar con más datos y análisis de los que normalmente son pedidos,
como: el índice y frecuencia de turbulencia del viento, la obtención de la velocidad
representativa del perfil de vientos, promedios de las velocidades en distintos intervalos
de tiempo del día, entre otras características, los que hacen que el análisis del terreno sea
más completo que los análisis que usualmente se hacen para caracterizar los
emplazamientos con potencial eólico.
En este caso, para realizar una correcta manipulación de los datos diarios, se construyó
una Macros de Excel programada tanto con Visual Basic como con código fuente de
Excel, capaz de mostrar los parámetros que describen de mejor manera al viento de
forma automática para las distintas bases de datos.
Los anemómetros tienden muchas veces a cometer errores de medición que estropearían
un análisis de datos que no incluyera una limpieza de los datos almacenados en las bases
de datos. Algunos de estos datos mal registrados pueden ser descubiertos mediante
46
simple inspección, por ejemplo, ver un valor mayor que 20 m/s de velocidad de viento
cuando en los instantes anterior y posterior son bastante menores, mientras que otros son
difíciles de captar como, por ejemplo, que de una hora a otra el promedio de viento haya
aumentado en más de 5 m/s. Debido a la relevancia y la frecuencia de estos sucesos
existen estudios que indican cómo se deben limpiar estos datos erróneos mediante una
serie de filtros que deben ser aplicados a los datos.
Existen pruebas tanto para las variables velocidad del viento, dirección del viento,
temperatura, presión barométrica y velocidad vertical del viento, siendo las primeras dos
las más importantes y que consideran una mayor cantidad de filtros que el resto de las
variables.19
Es importante destacar que al validar datos eólicos se deben observar los datos con una
frecuencia de 1 Hz y se registran en intervalos de 10 minutos. Los valores que se deben
registrar son: promedio de velocidad del viento, promedio de dirección del viento,
velocidad máxima ocurrida durante el intervalo de tiempo, velocidad mínima registrada
durante el intervalo de tiempo y desviación estándar de la velocidad del viento. Para este
proyecto la frecuencia máxima de muestreo del anemómetro y datalogger es 0,0166 Hz
que representa intervalos de un minuto, por lo que se realiza el registro de datos con la
máxima frecuencia disponible por el anemómetro. En el caso de los anemómetros
dispuestos por la empresa no se distingue la capacidad de escanear con registrar, por lo
que no era posible tener bases de datos tomadas cada 10 minutos con desviaciones
estándar de muestreos cada un minuto, sino que se debió registrar todo a la frecuencia
máxima lo que daba como resultado grandes bases de datos.
Debido a esta diferencia fundamental, los filtros que se utilizan están divididos para los
datos crudos o datos registrados cada un minuto y los datos promedios, que son los datos
filtrados, según los promedios de diez minutos realizados para cada estudio de
configuración del aerogenerador. Para los primeros, no se contó con una medición de
19
Pucci V., Bahamondes F. y Escobar R. (2010). Informe Ingeniería DICTUC: Buenas
prácticas en prospección Eólica y Solar. Disponible en Dictuc S.A. www.dictuc.cl
47
desviación estándar por lo que se consideró que los datos cada 1 minuto tienen menor
cantidad de filtros que los de 10 minutos.
Las tablas 7,8 y 9 muestran los valores críticos utilizados por los distintos criterios para
esta tesis. Si se quieren ver las tablas originales para ser usadas en validación de datos
generales se puede ver la norma o el informe del Dictuc citado como referencia de esta
información.
1) Pruebas de Rango
Son las pruebas que indican el rango máximo que pudiesen tener ciertas variables de
velocidad y dirección. Son las que se usan más comúnmente y son más accesibles
mediante una inspección visual del archivo. Es importante que la persona que realice el
chequeo sea consciente de que estos criterios están diseñados para la prospección eólica
de terrenos que cumplen con las características necesarias para las turbinas de eje
horizontal a ciertas alturas de medición, por lo tanto algunos serán cambiados.
Pruebas de Rango
Parámetro de muestra Criterio de Validación
Velocidad del viento: horizontal
• Promedio offset < Promedio < 25m/s
• Desviación estándar 0 < DE < 3m/s
• Ráfaga máxima offset < Max < 30m/s
Dirección del viento
• Promedio 0° < Promedio < 360°
• Desviación estándar 3° < DE <75°
• Ráfaga máxima 0° < Max < 360°
Para el caso de esta tesis el terreno difiere bastante de lo que sería un sitio óptimo para la
instalación de turbinas eólicas de eje horizontal, por lo que se varió algunos de los
48
rangos sugeridos por la norma y adaptar otros que ésta no considera. Las principales
diferencias entre las características del viento en Chagres y el viento que se usa para la
otra tecnología de turbinas es que en este caso se está midiendo desde una baja altura y
con una cantidad considerable de accidentes geográficos y arbustos alrededor, lo que
hace que haya una gran diferencia de velocidad entre las distintas alturas que se están
midiendo y que el viento presente mayor frecuencia en distintas direcciones que en otros
terrenos. Esto por un lado obliga, por ejemplo, a aumentar la diferencia de medición
permitida entre un anemómetro y otro, y también a considerar un rango más amplio de
diferencias de dirección entre una veleta y las siguientes.
2) Pruebas de relacionales
Estas pruebas llegan a sus resultados mediante la creación de restricciones que comparan
los valores de un anemómetro con otro. La cantidad de pruebas que se pueden hacer
dependen en gran medida de la cantidad de sensores que se tenga, por ejemplo, si en una
torre se utiliza solamente un anemómetro, éste no se podrá comparar con otro para ver si
la medición se encuentra dentro de un cierto rango esperable de acuerdo al terreno en el
que se está midiendo y a la experiencia anterior.
Las pruebas relacionales que define la norma están estructuradas de acuerdo a una cierta
configuración de la torre de medición (mediciones a 20, 40 y 60 m) que es la forma más
estándar y básica de medir los terrenos, pero como, en la actualidad han habido varios
cambios, y las torres están siendo estructuras de mayor altura, se necesita una medición
distinta, por lo que estos filtros deben ser solamente tomados como referencia y
adaptados para las características de la torre a la cual se aplicarán.
En la tabla 8 se muestran los filtros que se han elegido para este proyecto según sus
características. Los filtros originales pueden ser vistos en el documento citado
anteriormente.
49
Pruebas Relacionales
*Parámetro de muestra Criterio de validación
Velocidad de viento: horizontal
• Ráfaga máxima vs. Promedio Ráfaga máxima ≤ 2,5*Promedio
• Δ Promedio 18m/12m ** ≤ 2,0m/s
• Δ Diario máximo 18m/12m ≤ 5m/s
• Δ Promedio 18m/6m ≤ 4.2m/s
• Δ Diario máximo 18m/6m ≤ 7,5m/s
• Δ Máximo ≤ 2,0m/s
Dirección del viento
• Δ Promedio 18m/12m ≤ 90°
3) Pruebas de Tendencias
Son las pruebas que se basan en las tendencias pasadas que ha tenido un mismo sensor,
es decir, vigilan que no ocurran cambios bruscos o inesperados en las variables cuando
se cuente con un registro que ya lleva una cierta cantidad de tiempo tomando datos. De
estas pruebas se tomará solamente una, la cual indica que dentro de una hora de
promedio de viento no puede existir un valor promedio de 10 minutos que la sobrepase
por más de 5 m/s.
Pruebas de Tendencia
Parámetro de muestra Criterio de validación
Promedio de la Velocidad del Viento (todos los tipos de sensor)
• Cambio de 1 hora < 5.0 m/s
50
Una vez que se eliminan los valores entregados por estos filtros se obtendrán datos
confiables que pueden empezar a ser procesados para obtener resultados sobre la
caracterización del viento y la potencia de la máquina.
Además, los filtros deben ser realizados con contadores para que entreguen un resumen
de frecuencia final del número de datos filtrados para cada categoría y según cada filtro.
Los datos de viento son tomados en períodos de 60 segundos. Esta medición tiene una
frecuencia de registro 10 veces mayor que los estudios de prospección de terrenos
eólicos. Esto debido a que el viento presenta fluctuaciones bastante rápidas con respecto
al tiempo, lo que hace que los gráficos cada 1 minuto puedan parecer demasiado
irregulares para tener una apreciación del comportamiento del aerogenerador cuando se
toman intervalos diarios de medición. Por este motivo, se realizan promedios de las
mediciones cada diez minutos lo que muestra curvas más suaves con menores peaks de
velocidad para cada sensor. La idea del intervalo de diez minutos es trabajar con tiempos
que cumplan con los requerimientos de visualización y con los tiempos que recomiendan
las normas IEC, tanto para estimación de producción de potencia como para prospección
eólica. Hacerlo de esta manera facilitará la comparación de esta turbina con las
actualmente desarrolladas, y la curva obtenida por la turbina será fácilmente integrable
al procedimiento y resultados de todos los terrenos que han sido caracterizados según la
norma.
Con los datos promediados según los intervalos recomendados se construye un gráfico
de velocidad de viento de los distintos anemómetros respecto del tiempo, el cual es muy
efectivo para indicar cualquier desperfecto que pueda ocurrir en los anemómetros y la
causa de ese desperfecto. Con estos gráficos, durante el desarrollo de la tesis, se
descubrió que uno de los defectos más frecuentes se manifestaba cuando la carga de
alimentación del computador recolector de datos era poco regulada, por lo que hay que
51
Para estimar que tan bueno es un terreno en cuanto su potencial eólico, es fundamental
conocer la frecuencia de velocidad del viento del terreno, es decir, cuántas horas del día,
mes o año contiene a cada una de las velocidades del viento que presenta la distribución
del perfil de vientos. Esta frecuencia se debe obtener tanto para los datos registrados
cada un minuto, como para los registrados cada diez minutos, pues para la estimación de
la energía se utilizarán ambas frecuencias.
20
Mahyoub H.Al Buhairi (2006). A statistical analysis of wind speed data and
assessment of wind energy potential in Taiz-Yemen, Physics Department Taiz
Univeristy. Republic of Yemen. Recuperado de Assiut University:
http://www.aun.edu.eg
52
De la misma forma en que se obtienen las frecuencias de las velocidades del viento, el
programa diseñado registra la frecuencia en las direcciones del viento. Este registro es
tomado como precaución, debido a la importancia que tiene para sistemas eólicos de eje
horizontal, el cual no es simétrico respecto del eje vertical y, por lo tanto, debe
orientarse según la dirección del viento. Es por esta diferencia que la máquina puede ser
instalada en lugares de alta velocidad de viento, pero que pueden ser desechados por
tener una frecuencia de dirección muy desequilibrada entre los distintos puntos
cardinales. Estas zonas tienden a ser locaciones con cerros tipo valles o laderas, en
donde la dirección de las corrientes de viento dependen más de las condiciones
geográficas del lugar que de los cambios de presión producidos por la temperatura en las
diferentes partes de la atmósfera. Para este caso no es trascendental tener las frecuencias
de direcciones, ya que el aerogenerador funciona de igual forma con velocidades que
vienen desde distintas direcciones, simplificando en gran medida los procedimientos
para obtener la energía generada, sin embargo igual fueron calculadas por la importancia
que tienen en la caracterización del terreno.
53
Toda corriente de viento presenta turbulencias durante su viaje, debido a que las
presiones atmosféricas y las superficies terrestres interactúan de forma irregular,
haciendo que la velocidad del viento varíe notablemente con respecto a pequeñas
fracciones del tiempo. Esta turbulencia afecta el grosor de la capa límite en los perfiles
alares que funcionan según el principio de sustentación afectando su eficiencia, siendo
importante que en cualquier caracterización de una turbina queden claramente
especificados cuáles eran los índices de turbulencia del terreno en el que se hizo la
medición eólica.
Para cada medición realizada por los anemómetros se debe calcular la desviación
estándar y la velocidad promedio del período, para luego estimar la intensidad de
turbulencia (TI), que está dada por la siguiente fórmula:
, donde:
21
Ammoint (2011), Measurement equipment for the wind industry Recuperado el 27 de
Agosto del 2011, de http://www.ammonit.com
54
Esta turbulencia debe ser calculada mediante un escaneo de frecuencia 1 Hz por parte
del anemómetro. Y en un período de 10 minutos.
En el archivo de datos en donde se realizan todos estos cálculos, es muy importante que
para cada columna de valores de las velocidades del viento a las distintas alturas y en
distintos promedios de estas velocidades se obtenga el valor máximo y el valor mínimo
de los valores promedios registrados, los cuales son una primera aproximación de cómo
fue el día, para determinar si conviene analizar las revoluciones por minutos y la
potencia para ese día; además muestra si hubo defectos en las mediciones aquel día,
siendo sobre todo la velocidad máxima un filtro que se debe imponer a los datos.
Los análisis mensuales, que con el tiempo finalmente llegan a ser informes anuales, son
importantes porque permiten evaluar económicamente un proyecto eólico. Esto
posibilitará comparar la turbina caracterizada con otras turbinas en el largo plazo y así
poder determinar qué turbina se comportaría de manera más eficiente de acuerdo a las
características de viento de un cierto terreno, especialmente cuando la turbina a
desarrollar cuenta con una geometría bastante distinta, la cual puede tener una curva de
potencia ideal para diferentes características de viento del terreno.
Es importante que se realice este análisis mensual para todo el período que toma la
caracterización de la turbina. Recordar que el objetivo de este análisis de largo plazo es
poder tener una apreciación general tanto del terreno, en todo el tiempo, como una
55
apreciación general de los vientos durante el período en que la turbina estuvo sometida a
estudio.
Los análisis mensuales se realizan mediante un programa hecho en Excel Visual Basic
en donde, ingresando los datos mediante una planilla Excel, éstos son ordenados,
modificados y graficados en la medida que muestren la mayor cantidad de características
del perfil de viento del terreno a utilizar. Estas características se pueden resumir en las
siguientes:
Son tablas que contienen un resumen con la información más importante del terreno,
como el período de la medición, los anemómetros y dataloggers usados, los promedios
de velocidad a distintas alturas, velocidades máximas, los parámetros para las curvas
Weibull de ese período, los coeficientes de cizalle y rugosidad de terreno para el
período.
Son gráficos que muestran la ocurrencia de los registros más importantes. Con la
frecuencia del viento se puede estimar la potencia que produciría una turbina; con la
dirección, se puede ver de dónde viene el viento y por determinar cómo orientar un
parque; y, con la turbulencia, es posible estimar cuál es la calidad de ese viento lo que
modificará en alguna medida la curva de potencia que se usará en conjunto con la
frecuencia para estimar el AEP (Annual Energy Production) del terreno.
Estos estudios tienen directa relación con el Log Law Profile y el Power Law Profile
utilizados para el cálculo estimativo de extrapolación de datos de velocidad. La
extrapolación de los datos está basada en determinar el valor de un término que permita
proyectar la velocidad medida a otras alturas que sean importantes de conocer. Mientras
mayor sea el término significa que mayor es la diferencia de las velocidades con la
altura. Se debe caracterizar el viento del terreno con ambos factores, porque éstas son las
formas más utilizadas y así se puede comparar mejor con otros sitios ya certificados.
(
( )
(
[6]
( )
Dónde: (
Este cálculo solamente se puede realizar cuando se tienen dos o más anemómetros bajo
la misma línea vertical, y se realiza para toda la combinación de alturas disponibles. En
esta evaluación, al tener 3 anemómetros, se encontró tres parejas distintas de pares, por
lo que se obtuvo tres valores los cuales promediados entregarán el cizalle específico del
terreno. Esto debe ser realizado cada mes para así llevar un registro preciso de la
variación del cizalle de acuerdo a las épocas del año.
22
Measnet (2009). Measnet Procedure: Power Performance measurement procedure,
Version 5, December Recuperado el 7 de Marzo del 2011, de http://www.measnet.com
57
más sensible en magnitud que el coeficiente de cizalle y usa una proyección logarítmica
en vez de exponencial.
( (
(
(
( )
(
(
Dónde: (
Esta curva es una de las más importantes; entrega información sobre la velocidad del
viento en metros por segundo [ y la dirección en grados sexagesimales para
cada una de las horas del día de todo el mes.
Cuando se realiza el grafico utilizando la velocidad del viento, la curva permite conocer
durante qué horas del día existe mayor probabilidad de encontrar vientos altos y eso
puede llevar a definir horarios de trabajo y horarios de pruebas eléctricas para la
construcción y desarrollo de la turbina.
f) Obtención de los parámetros de la curva Weibull que más se asemeje a los datos:
Esto permitirá comparar el comportamiento del viento en Chagres con otros terrenos
que ya hayan tenido una evaluación eólica.
Según la norma IEC61400-12 las velocidades de viento necesarias para poder hacer una
adecuada medición deberían estar en el siguiente rango: desde velocidades 1 m/s
menores que la velocidad en la que se comienza a generar potencia (cut-in), hasta
velocidades 1.5 veces mayores que el 85% de la potencia nominal de la turbina. Para el
rango inferior claramente no hay problema, pues en los resultados se presentaron
momentos en que existe velocidad de viento menor a 0.5 m/s, lo que sin lugar a dudas es
menor que lo necesario.
menor de los casos. Estos datos deben estar presentes un período de tiempo mínimo de
30 minutos.23
Una vez realizado el estudio del terreno para estimar si tiene potencial para caracterizar
una turbina, se debe hacer una estimación adecuada de la cantidad de potencia que trae
la corriente de viento del terreno. De esta manera se podrá estimar la eficiencia de
conversión de energía eólica en eléctrica del aerogenerador.
Para estimar la energía producida por una corriente de viento lo mejor es estimar un
perfil de velocidad del viento para todo el rango del área barrida por el aerogenerador.
Para obtener este perfil de manera precisa habría que colocar una cantidad ilimitada de
anemómetros cubriendo el rango de alturas anteriormente descrito, pero como eso no es
posible, para esta tesis, se colocaron 3 anemómetros a distancias de 6 m entre ellos y con
esas mediciones se pretende realizar una proyección de las velocidades intermedias que
ocurren en el perfil de vientos.
23
International Electrotechnical Commission (2005). Power performance measurements
of electricity producing wind turbines. (61400-12-1 Ed1). Holanda: Autor.
24
Ammoint (2011): Measurement equipment for the wind industry at
www.ammonit.com.
60
Las turbinas de eje vertical comenzaron su desarrollo con modelos pequeños y ahora han
ido aumentando su tamaño hasta buscar un límite con el de las horizontales. Esta menor
envergadura hace que sea más económica la obtención de un perfil del viento bien
definido.
Para estimar la velocidad a alturas no censadas del perfil, se tienen tres alternativas de
interpolación; la primera, es usar una interpolación lineal identificando el coeficiente de
aumento de velocidad real que ocurre a mayores velocidades; luego es posible usar las
61
Para el desarrollo de esta tesis se cuenta con una turbina de baja altura, el extremo
inferior del álabe está a 6 m sobre el nivel del suelo, mientras que el extremo superior se
ubica a 18 m sobre el nivel del suelo, esto hace que sea accesible económicamente
colocar anemómetros a las alturas más importantes del perfil de viento, extremo inferior,
extremo superior y al centro del perfil. Además se debe tomar en cuenta la altura del
perfil que debe ser aproximado, en este caso con 12 m que es una distancia pequeña, en
donde se pueden obtener buenas estimaciones. Debido a que la interpolación es a poca
distancia, y con varios datos indica que una buena aproximación lineal sería ideal para
estimar la velocidad del perfil completo, pero esto debe ser corroborado con un análisis
ECM de los puntos.
Antes de analizar los datos se deben descartar los datos de velocidad que tengan un valor
menor al de cut-in de la turbina, pues éstos pueden distorsionar los resultados de las
velocidades que realmente se quiere estimar. Si no se conoce el valor de cut-in se debe
estimar uno de la literatura, experiencia y práctica.
1) Estimación lineal
La estimación lineal busca evaluar cada una de las velocidades en que no se tiene
anemómetro a partir de una o más de las velocidades en las que sí se tienen
anemómetros para obtener un perfil completo de velocidades.
25
J.F. Manwell; J.G. McGowan y A.L. Rogers (2002). Wind Energy Explained Theory,
Design and Application. University of Massachusetts, Amherst, USA.
62
A medida que se tiene una mayor velocidad relativa del perfil, la diferencia de
velocidades entre anemómetros ubicados a diferentes alturas va aumentando. Se debe
encontrar la tasa promedio de esa variación respecto de la velocidad desde cada uno de
los anemómetros hacia los otros. La tasa está dada por la siguiente expresión:
( ( (
Esta tasa ( se puede encontrar resolviendo una regresión lineal entre ambos datos, o
bien encontrar el valor que tenga el menor error cuadrático medio de la diferencia del
valor estimado con el valor real.
{∑ ( ( ))}
( )
Esta regla se usa encarecidamente para determinar tanto el perfil de velocidad de viento
de un sitio, como para determinar la velocidad de viento que se tendrá a una determinada
altura fuera del rango medido, es decir, es una herramienta que se usa para interpolación
y extrapolación. Para este caso solamente interpolará, puesto que debido a las bajas
alturas máximas que se tienen en el proyecto es económicamente viable medir la
velocidad del viento para todo el rango que recorren las aspas.
63
Para obtener las estimaciones de la velocidad del viento a distintas alturas primero se
debe estimar el coeficiente de rugosidad ( ). Para esto se debe clasificar el terreno
según sus características topológicas, su vegetación y otros obstáculos que pudiesen
impedir el flujo laminar del viento. El coeficiente de rugosidad determina las
variaciones de la velocidad del viento según la altura de medición mediante una relación
logarítmica. Se calcula principalmente mediante dos formas.
Una forma de obtener el valor del coeficiente de rugosidad es mediante el uso de la tabla
mostrada en el anexo G, esta forma no es muy utilizada para estimar el valor, debido a
que está demasiado sujeta a la interpretación personal sobre las características del
terreno, lo que la hace imprecisa. Pero sí es bueno tenerla en cuenta como una referencia
para verificar los resultados.
Dónde:
(
64
De esta manera se deben compararon los resultados con los obtenidos por la tabla del
anexo y los obtenidos mediante el cálculo anterior.
Una vez con el coeficiente definido se puede reemplazar en [11] para obtener la
velocidad del viento que se estima desde los anemómetros extremos al anemómetro
central.
Para obtener el valor de según datos obtenidos, se debió utilizar la igualdad superior
despejada que queda de la siguiente forma:
( )
( (
( )
Una vez hemos estimado la velocidad a la altura del anemómetro central se resuelve la
misma minimización del error que para el caso lineal utilizando la ecuación [9] donde se
obtienen los mejores valores de para estimar la velocidad del anemómetro intermedio.
( ( ( )
En este caso no existen tablas para el valor de alpha ( , por lo tanto se debe calcular
primero su valor colocando más de un anemómetro en el lugar de medición. El número
de valores que se obtendrá de la velocidad intermedia a estimar es igual al número de
anemómetros con el que se mida el perfil.
Teniendo en cuenta esta fórmula es que se realizó la misma optimización que para el
caso anterior. Pero en este caso se debe encontrar aquel valor de que minimiza la
media del error. De esta manera podremos tener un método para comparar cómo se
comporta la predicción realizada con esta estimación con las otras predicciones hechas
con anterioridad. La fórmula utilizada para estimar es la siguiente:
( ( ) ( ( )
( (
Una vez determinados los coeficientes de cada función que minimizan el error
cuadrático medio del error, se comparó la combinación de todos los anemómetros que
llevan al menor error de la estimación del anemómetro central. Para esto, se minimizó
cada una de las estimaciones, ponderando en la fracción la distancia total a la que se
encuentran del anemómetro ( . Esto se resuelve de la siguiente manera:
{ [ ∑ ( )]}
26
J.F. Manwell; J.G. McGowan y A.L. Rogers (2002). Wind Energy Explained Theory,
Design and Application. University of Massachusetts, Amherst, USA: [Wiley]
66
Este resultado entrega un vector que contiene los valores de error cuadrático medio
de las mejores aproximaciones calculadas mediante el método anterior. Finalmente, el
menor valor del vector mostrará cuál de todas las aproximaciones estima de mejor
manera el perfil y, ésta será la utilizada para estimar todo su rango.
Esta curva de perfil tendrá como entrada todos los datos de anemómetros que se
encuentren en él, y su integral dividida por la altura tendrá como resultado la velocidad
representativa del perfil.
Para crear una función que caracterice el perfil es primordial definir aquellos
coeficientes que valoricen las aproximaciones de los anemómetros de acuerdo a la altura
que se quiere conocer. Estos coeficientes se calcularán de acuerdo al inverso de la
distancia del anemómetro con el punto que se quiera medir: por ejemplo, si se quiere
medir a los 17m, entonces el anemómetro de menor distancia, y por ende de mayor
relevancia, es el de 18 m.
Luego de esto se debe calcular la curva que permita estimar de mejor manera la
velocidad en ese lugar con cada uno de los anemómetros. Para eso se usa el método de
aproximación calculado anteriormente y se determinan los coeficientes de aproximación
óptimos desde todas las distancias. Con esos puntos fue posible armar curvas que
describan a los coeficientes para todas las alturas de perfil.
67
Combinando ambos parámetros se obtuvo una ecuación que estima todas las velocidades
en el perfil. Luego la integral de esta curva en su recorrido entregará el valor de la
velocidad equivalente.
Una vez definida la velocidad del perfil es muy importante saber cuál es la cantidad de
energía y potencia que tiene el viento que está pasando a través de la máquina, valor que
entrega el límite superior de la energía extraíble desde la turbina. Lo anterior junto con
los resultados de la potencia empírica, señala la eficiencia de la máquina que se está
analizando.
Para poder conocer la potencia del viento es importante conocer de manera precisa la
densidad del aire del lugar y ver la variación de ésta según las distintas alturas y épocas
del año.
Posteriormente se deben observar las diferencias de la densidad del aire de acuerdo a las
distintas épocas del año y horas del día, dependiendo de las horas en las que se genere
68
potencia que son las relevantes. Estos valores se deben comparar con los publicados por
la el documento confeccionado por Measnet27.
curva de potencia adaptando los resultados según una fórmula presentada. Es necesario
señalar que la norma IEC61400-12 señala un menor rango de densidades en donde
. Pero como el último documento mencionado fue escrito en una
Finalmente, una vez definidos todos estos parámetros, será posible obtener de manera
precisa la potencia en la corriente de viento que pasa a través del aerogenerador.
27
Measnet: Power performance measurement procedure. Version 5 December – 2009.
69
Dado que el proceso de obtención de una curva de potencia está directamente ligado al
proceso de optimización del aerogenerador, los resultados de curvas de revoluciones
ayudarán a comparar los comportamientos obtenidos en cada configuración para ver cuál
equipo es más eficiente.
Para obtener los resultados de las revoluciones del prototipo hay que basarse
principalmente en dos métodos. Uno es de rápida inspección, que se lleva a cabo en los
análisis diarios y que indica la performance del aerogenerador y el correcto
funcionamiento de todos los sensores de medición asociados a la curva de revoluciones.
Este análisis debe ser automático e integrado a las bases de datos y a los programas que
analizan en primera instancia la información. El segundo análisis basa su procedimiento
en la norma de curvas de potencia para aerogeneradores de eje horizontal y entrega
resultados en su mismo formato, siendo una primera aproximación a la obtención de la
curva de potencia al describir la velocidad angular del eje de transmisión.
Con los datos necesarios para obtener la curva de revoluciones del aerogenerador, se
realizó una adaptación a las bases de datos según los programas realizados para el
70
Este análisis fue realizado tomando en consideración el método descrito por la norma
IEC 61400-12 sobre “Wind turbine power performance testing”, aplicado también por la
institución Measnet en su documento de mejora de la norma actual, y es considerado
71
como un buen método para aproximar curvas de potencia.28 Por lo que se continúa
usando en la actualidad. En la norma, el método no se aplica para determinar una curva
de revoluciones ni con carga ni en vacío, pero para efecto de este procedimiento es
importante agregarlo debido a que nos permite comprender y optimizar el tipo de
aerogenerador que se está caracterizando, lo que afectará muchas decisiones futuras
como la estimación inicial de la curva de torque eléctrico a integrar al sistema.
El método consiste en tomar los datos empíricos mediante sensores que guarden la
información con el tiempo real de la medición coordinándose entre ellos para luego
formar bases de datos que contengan en sus columnas toda la información necesaria para
obtener los resultados. Luego esos datos deben ser filtrados según su clase para no
procesar ningún error de parte de los sensores (anemómetros, tacómetros, analizador de
potencia). Estos filtros deben ser tanto programados como visuales, y su fin es buscar
anomalías en los datos. La frecuencia de muestreo debe ser la mayor posible que
permitan los sensores, en general se recomienda una frecuencia de 1 Hz con un registro
de 1 y 10 minutos para este caso.
Una vez determinada la base de datos correcta, se puede dividir todos los datos en
muestras representativas de la velocidad del viento ordenándolas de forma ascendente.
Se consideran representativos intervalos de 0.5 m/s, que son los mismos intervalos que
se realizan en los estudios de prospección de terrenos eólicos. Una vez que se tienen
separados los datos según la velocidad del viento, se calcula el valor que tiene el menor
error cuadrático medio con respecto a los datos, para fijarlo como el valor de las
revoluciones para esa velocidad representativa de viento.
28
Power coefficient measurement on a 12 kW straight bladed vertical axis wind Turbine
J. Kjellin*, F. Bülow, S. Eriksson, P. Deglaire, M. Leijon, H. Bernhoff
72
Primero se define el término Tip Speed Ratio (TSR) como el cuociente entre la
velocidad tangencial en un instante de tiempo del álabe y la velocidad del viento del
terreno para ese mismo período. Éste término está directamente relacionado con el
coeficiente de potencia de la turbina, y por lo tanto con su eficiencia de transformación
de la energía de la corriente eólica que circula a través de ella.
Para cada geometría de turbina existe un valor de TSR tal, que el coeficiente de potencia
es máximo para todo el rango de las velocidades de viento del terreno. Esto implica que
ese valor es el valor en que la turbina opera con la mayor eficiencia y éste depende en
gran medida de los principios físicos con los que opera la turbina, de su diseño y del
número de álabes, especialmente si aplican las fuerzas de sustentación. Por ejemplo,
para las turbinas Savonius los TSR óptimos son similares a 0.3 y los TSR máximos están
en valores cercanos a los 1.229, debido a que este tipo de turbina funciona por el empuje
del aire, es difícil que los álabes superen la velocidad del viento que los impulsa. Esto
hace que los Savonius sea máquinas de baja revolución y alto torque. Para los Darreius
el valor óptimo del TSR ronda los 3 y los TSR máximo en la práctica tienen un valor
aproximado de 4.330, presentando velocidades de revoluciones por minuto bastante
mayores que los Savoniuns. Por este motivo los Darreius son máquinas de alta
revolución y bajo torque.
Para darse cuenta qué tipo de aerogenerador se está caracterizando, se debe obtener la
curva de revoluciones según el método de los contenedores explicado anteriormente para
el aerogenerador en vacío, es decir, conectado al motor eléctrico pero sin carga resistiva
que lo detenga. De esta curva se debe desprender el comportamiento del TSR para luego
29
Ragheb M. (2011). Optimal Tip Speed Ratio. Recuperado el 4 de Mayo del 2011 de la
University of illinois at urbana-champaign http://www.illinois.edu
30
J. Kjellin et al. (2010). Power coefficient measurement on a 12 kW straight bladed
vertical axis wind Turbine, Renewable Energy, 36 (2011) 3050-3053. Uppsala, Sweden.
73
compararlo con los resultados de cada uno de los grandes diseños de turbinas verticales
y poder clasificarla como tal.
La base para optimizar el aerogenerador es conocer la magnitud del torque que se debe
asignar al aerogenerador según la cantidad de viento que exista en ese momento, pero
como el analizador del drive no viene con un anemómetro integrado, se utiliza el valor
de las revoluciones por minuto para definir el torque resistivo del motor. En la sección
anterior se describió el método de dependencia entre las revoluciones y la velocidad del
viento con un pequeño retraso debido a la inercia del sistema que por lo general es
considerada como despreciable, debido a que incluso en intervalos de promedio cada 1
minuto es muy difícil de apreciar. Es importante notar que la relación entre las
revoluciones y la velocidad del viento depende fuertemente de la curva de torque
ingresada por lo que esta optimización es un proceso iterativo y, para cada curva de
torque programada, la curva de revoluciones versus velocidad de viento debe ser
recalculada.
Para programar la curva de torque en este proyecto se utilizó programas suministrados
por la empresa proveedora del motor y el drive ABB S.A. Ambos componentes son
controlados mediante los programas DriveWindows y DriveAP, que son capaces de
manipular el sistema DTC del variador de frecuencia y por ende, programar una curva
74
de torque cualquiera en el drive, para así alcanzar el punto de MPPT del aerogenerador.
En el Anexo H se puede ver la interfaz de los programas.
El proceso comienza con la obtención de la primera curva de torque. Para llegar a esta
curva se puede realizar el siguiente procedimiento: El método se utiliza en situaciones de
vientos menores a 1 m/s y consiste en accionar el generador eléctrico como motor.
Como éste está conectado a un variador de frecuencia, es posible controlar la velocidad
de giro del motor, lo que hace que el proceso sea seguro. Una vez accionado, se hizo
girar al motor en el mismo sentido de su diseño, midiendo todas las variables mecánicas
y eléctricas de ese proceso, hasta llegar a recorrer todo el rango de velocidades angulares
obtenidas con la prueba en vacío. Los datos de velocidad angular con la potencia
consumida por el motor para hacerlo girar a esa velocidad, entregan una curva de torque
inicial que tendrá la misma forma que el torque que el viento aplica al aerogenerador
cuando está en vacío y por ende, será poco eficiente debido a que su TSR es el máximo
( que está alejado del óptimo. La curva obtenida tendrá la misma forma que la del
viento, pero será menor a ésta manteniendo la seguridad en el proceso.
El procedimiento entregará una primera curva que servirá para ver la tendencia de la
curva, pudiendo ser cuadrática, cúbica, etc, y para probar el funcionamiento de todos los
componentes eléctricos, especialmente la velocidad de respuesta y magnitud del control
de torque lo que indicará qué tan suave trabajará la turbina eólica.
Para los aerogeneradores que funcionen bajo el principio de sustentación, es muy
importante definir una velocidad llamada velocidad de cut-in o de corte de entrada. Esto
se hace debido a que algunos tipos de aerogeneradores, como el Darreius, tienen
dificultades para salir de la inercia inicial31 por lo que es preferible tenerlo sin torque
eléctrico hasta que entra en su velocidad de sustentación de diseño que depende de la
velocidad del viento y de las revoluciones por minuto a las que está girando. Las
31
Vásquez Leiva M. (2007). Diseño de una turbina eólica de eje vertical con rotor
mixto. (Tesis de maestría no publicada). Universidad Austral de Chile, Valdivia, Chile.
75
velocidades de cut-in son prácticamente iguales para todas las turbinas tanto verticales
como horizontales, siendo de 3 m/s para las turbinas de pequeña escala y de 4 m/s para
las de gran escala. Para programar una curva de cut-in se debe tener en cuenta la curva
de revoluciones obtenida de la metodología anterior. Luego se observan las revoluciones
del motor a los 3 o 4 m/s y, finalmente, se programa en el drive para que se encienda el
control de torque cuando las revoluciones alcancen la magnitud de la curva de rpm. Para
el valor de cut-out se debe cortar el suministro de torque eléctrico de parte del drive o
bien hacer que éste decaiga de manera de mantener la potencia máxima, pero saliéndose
del TSR óptimo que tiene la curva de torque del aerogenerador. La forma de programar
el valor de cut-in y de cut-out puede variar según la marca y modelo del drive que se
utilice para programar la curva y de los parámetros que sea capaz de medir.
Una vez solucionado el problema del control de torque inicial se comenzó el análisis de
los resultados de los parámetros eléctricos que influyen en la generación eléctrica y con
ello obtener las curvas características de la configuración del aerogenerador.
comparar las frecuencias de los datos de viento y las revoluciones con las de potencia.
Estos gráficos debiesen tener un comportamiento similar lo cual se debe verificar.
Con este método se logró obtener la dependencia de la velocidad del viento con la
potencia. Se separaron las mediciones en intervalos de velocidad de viento de 0,5 m/s
según las mediciones de 1 y 10 minutos32. Estos valores serán graficados y se entregarán
las tablas correspondientes con los valores para la velocidad de viento, las revoluciones
por minuto, la potencia, el coeficiente de potencia, el TSR, las turbulencias y el
coeficiente de cizalle, que son los factores más influyentes en las variaciones de la curva
de potencia.
Según el IEC los resultados se deben presentar en dos tablas denominadas Tabla A y
Tabla B. La diferencia entre ambas es que la primera cuenta con todos los datos
obtenidos excepto los eliminados por mantenimientos de la turbina, efecto de capas de
hielo en los álabes y dirección fuera del rango de operación definido, mientras que la
Tabla B se deben excluir además todos los datos en donde la turbina entre en régimen de
cut-off.
Para los casos en que los datos de velocidad del viento no alcancen la velocidad de cut-
off, se puede extrapolar linealmente la curva, tal que para vientos de mayor intensidad se
mantiene la potencia obtenida según el último punto de medición de la curva, la que será
plana al final. Este supuesto será considerado como válido siempre que se logre llegar a
32
Measnet (2009): Power Quality Measurement Procedure, Version 4.
77
Una vez que se obtuvo a la curva de potencia que rige a la turbina, es posible obtener la
energía que la turbina generará según las condiciones de viento específicas de un
terreno. Finalmente, para las turbinas eólicas, la energía generada es el flujo principal de
ingresos de un parque eólico.
Los procedimientos antes realizados permiten tener todas las bases para realizar este
cálculo. Se tomó la distribución de viento perteneciente a un período significativo de
tiempo en donde se haya mantenido una configuración y curva de torque específicos
para una turbina. Esa distribución se multiplicó por la curva de torque obtenida mediante
el procedimiento descrito. Ambos resultados se encuentran en intervalos de igual
magnitud por lo que solamente se debe multiplicar la potencia por el tiempo que estuvo
actuando para obtener la energía en kW-h generada por las condiciones. Ese valor debe
78
ser comparado con la energía medida por el analizador de corrientes trifásicas, y de esa
manera, se puede observar si la metodología se aproxima a la inyección real de energía
que hubo en el sistema durante el período de análisis.
5 Resultados y Discusión
5.1.- Análisis de las curvas de viento obtenidas a partir de los datos en terreno
Para lograr esto, los cálculos deben ser desarrollados iterativamente con bastante
frecuencia y, como la unión de bases de datos tiene un formato similar cada día, se creó
un código de programación en Visual Basic y Macros de Excel que es capaz de realizar
los cálculos y la gráfica de las variables exigidas por el método, usando los datos de
entrada correspondientes. Este programa fue diseñado, en el marco de esta tesis, para ser
adaptado a otras bases de datos eólicas con formatos de presentación distintos. Pero
solamente la obtención del perfil representativo de viento debe ser realizada en algún
programa que tenga presente análisis numérico matemático.
80
Para llevar un control preciso del estado de viento y del correcto funcionamiento de los
anemómetros, se revisó prácticamente todos los días el estado de las curvas de viento
mediante la metodología expuesta. Esto permitió descubrir que los meses del año en que
se encuentran mayores promedios de viento son: Noviembre, Diciembre, Enero y
Febrero, meses en los que se podían medir velocidades máximas de 9.5 m/s a 12 m de
altura y promediadas cada 10 minutos. Los otros meses eran caracterizados por algunos
días de velocidad un poco menor y otros de velocidad muy baja, entre 5 y 6 m/s
promediadas cada 10 minutos. Esto hace que se necesite una mayor cantidad de días
para hacer las pruebas y conseguir una muestra significativa de vientos de alta
velocidad. Estas velocidades no permitieron obtener la caracterización completa de una
curva de potencia estándar, pues, por lo general, la velocidad de viento nominal de las
turbinas se encuentra cerca de los 12-13 m/s para los VAWT.
Por el motivo antes mencionado, este estudio requiere realizar un análisis diario, ya que
de esa manera se tiene un registro de datos muy adaptables para hacer los requeridos. Se
pueden juntar los días de acuerdo a las configuraciones o simplemente comparar días
diferentes para una misma configuración. Además un día es un período de tiempo
aceptable para no tener demasiados archivos dentro de la caracterización de una
configuración.
En el siguiente apartado se verán los distintos análisis que se realizaron durante las
mediciones diarias y las mediciones por configuración, durante cada período de
medición se mostrará la utilidad de cada resultado en el desarrollo de optimización del
equipo y de la caracterización del sitio de la instalación.
Los filtros impuestos a los datos de viento programados en los archivos fueron bastante
útiles en entregarnos información sobre cuál estaba siendo la calidad de los registros por
parte de los anemómetros. Se programan contadores de frecuencia que resumían la
81
Resumen
Altura(m) Errores % Total
v6m 177 2.1
v12m 149 1.7
v18m 108 1.3
De este resultado se puede observar que un 2% de los datos presenta errores según los
filtros definidos, valor que parece razonable. Los datos con anomalías son eliminados de
la base de datos para no incluir errores al cálculo de los resultados posteriores.
8
Velocidad de viento y Rpm
7 Velocidad vs Tiempo
6
5
4
3 V6m
2
V12m
1
V18m
0
1:01:05 PM
2:01:05 PM
3:01:05 PM
4:01:05 PM
5:01:05 PM
6:01:05 PM
7:01:05 PM
8:01:05 PM
9:02:06 PM
1:00:04 AM
2:00:04 AM
3:00:04 AM
4:00:04 AM
5:01:05 AM
6:01:05 AM
7:01:05 AM
8:01:05 AM
9:01:05 AM
12:00:04 AM
10:01:05 AM
11:01:05 AM
12:01:05 PM
10:02:06 PM
11:02:06 PM
Tiempo
Se puede apreciar que la distancia entre ellas va aumentando cuando el viento relativo
del perfil aumenta, es decir, a bajas velocidades de viento las curvas están muy pegadas
pero a altas velocidades estas tienden a separarse más una de otras.
12
Frecuencia de la velocidad del viento [1 min]
10
Frecuencia [hr]
6 Frec6m(hr)
4 Frec12m(hr)
Frec18m(hr)
2
0
13.5
10.5
11.5
12.5
14.5
15.5
4.5
0.5
1.5
2.5
3.5
5.5
6.5
7.5
8.5
9.5
En el gráfico se puede observar que existe una gran cantidad de horas con viento entre 0
y 1 m/s. Esto ocurre debido a la baja altura de las mediciones y a la gran diferencia de
velocidad existente entre los vientos de la tarde con respecto a los vientos en la mañana,
periodo donde se acumulan muchos datos de velocidad cero. Este factor de diferencias
entre velocidades de viento en la tarde y en la mañana es determinante para el cálculo de
la potencia obtenida, ya que este terreno entregará una menor cantidad de tiempo de
vientos altos que aquellos terrenos que presentan un perfil de viento constante. Esto hará
que se necesite una mayor cantidad de días para poder completar el tiempo necesario
85
para que la curva de potencia reúna los datos suficientes para ser fidedigna, lo que
retrasará los tiempos del proyecto en cuanto a la actualización de sus configuraciones.
20
18 Frecuencia de dirección
16
14
Frecuencia [%]
12
10 Frec6m(%)
8
Frec12m(%)
6
4 Frec18m(%)
2
0
100
120
140
160
180
200
220
240
260
280
300
320
340
0
20
40
60
80
Cada terreno posee flujos más o menos turbulentos de acuerdo a la zona en donde se
realice la medición, que puede ser el mar o zonas con mayor volumen obstáculos; y a
factores climáticos que puedan ocurrir en el lugar, como la nubosidad variable y
tendencia a remolinos entre otras. La turbulencia afecta en gran medida a la
generación33. Por lo que es importante el estudio de turbulencias en un terreno para
determinar la incertidumbre sobre la curva de potencia obtenida y el cómo se comportará
ésta al ser llevada a otros terrenos con características de turbulencias distintas.
33
U. Bunse, H. Mellinghoff, O. Haack (2006). Uncertainty of annual energy production
for a specific turbine model based on a set of IEC 61400-12 measurements, German
Wind Energy Insitute GmbH (DEWI).
87
1.8
1.6 Gráfico de Turbulencias
1.4
Turbulencias
1.2 Turbulencia 6m
1
Turbulencia 12m
0.8
0.6 Turbulencia 18m
0.4 IEC Category A
0.2 IEC Category C
0
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
31
Velocidad del viento [m/s]
Figura 19: Resultado de las turbulencias para las velocidades del perfil.
Este gráfico muestra la cantidad de turbulencias promedio en intervalos de medida 0,5
m/s de largo. El gráfico sigue las tendencias normales mostradas en los estudios
eólicos34, a velocidades bajas las turbulencias tienden a ser muy altas debido a que los
cambios de velocidad cero a velocidades superiores a cero tienen una alta desviación
estándar. Esto se da en parte por la inercia del anemómetro a velocidades pequeñas,
luego el valor comienza a disminuir terminando por debajo de los estándares
recomendados por la IEC, esto señala que el terreno a utilizar no presenta grandes
turbulencias por lo que no hay una distorsión en los datos, sólo la normal producida en
la mayoría de los terrenos.
En el anexo J, se puede ver una comparación entre las curvas de turbulencia obtenidas a
las distintas alturas de Chagres con otras curvas obtenidas en un lugar costero típico del
litoral de la Cuarta Región de Chile llamado LV2, para mostrar una comparación de lo
similar que son las turbulencias en ambos sitios.
34
Measnet (2009). Measnet Procedure: Power Performance measurement procedure,
Version 5, December.
88
La velocidad equivalente máxima alcanzada en el período fue de 9.2 m/s para el perfil
en promedios de 10 minutos. Esta velocidad no alcanza a ser la velocidad nominal
estimada de la turbina entre 12 y 13 m/s, por lo tanto en este terreno no se logrará
modelar el perfil completo del aerogenerador.
Se recomienda hacer estos análisis con los datos mensuales, para reconocer
características macro, tanto anuales como mensuales del terreno. Estas características
macro describirán diversas acciones que pueden ayudar a la logística o resolución de
problemas durante la caracterización del aerogenerador, como la posibilidad de mejorar
el rango de velocidades de viento, y tiempos aptos para trabajar en cambios de
configuraciones. Cada una de estas ventajas serán analizadas según cada uno de los
puntos a continuación.
89
La tabla 13 resume las características más importantes del terreno según lo requerido
para caracterizar la turbina. Está hecha para ser lo primero que se entrega como informe
de resultados sobre el sitio.
Tabla 13: Tabla de resultados con las características mensuales más importantes.
REPORTE VIENTOS
ANTECEDENTES
PERÍODO 1 de Octubre al 31 de Octubre
LOGGER Em50
INTERVALO DE MUESTREO 1 minuto
NUMERO DE REGISTROS 44581
% DE RECUPERACIÓN DE DATOS 94
PROMEDIO A 6M [m/s] 2.6
PROMEDIO A 12M [m/s] 3.2
PROMEDIO A 18M [m/s] 3.8
VALOR MAXIMO 6M [m/s] 7.8
VALOR MAXIMO 12M [m/s] 9.2
VALOR MAXIMO 18M [m/s] 10.7
DESVIACIÓN ESTANDAR 0.59
WEIBULL K 13.01
WEIBULL C 0.29
Coeficiente de Rugosidad (m) 0.91
Alpha ( ) 0.35
90
Al igual que para los análisis diarios y de configuración se obtienen los mismos gráficos
Tabla 14: Comparación de los factores del coeficiente de cizalle ( en distintos sitios
Chagres 2011
Febrero Junio Agosto
CR Alpha ( CR Alpha ( CR Alpha (
6m-12m 0.948 0.424 1.100 0.563 0.912 0.388
12m-18m 1.082 0.434 1.065 0.421 1.025 0.391
6m-18m 0.990 0.427 1.097 0.510 0.942 0.390
Promedio 1.007 0.428 1.087 0.498 0.959 0.390
LV2 2011
Junio Agosto
CR Alpha ( CR Alpha (
10m-20m 0.0162 0.0824 0.0099 0.0753
20m-30m 0.0017 0.0559 0.0027 0.0593
10m-30m 0.0087 0.0726 0.0066 0.0694
Promedio 0.0089 0.0703 0.0064 0.0680
91
Luego es apreciable que existen diferencias en los factores según las épocas del año, esto
está en directa relación con las diferencias en velocidad que encontradas durante el año
para cada terreno, pero las diferencias en un mismo sitio no son tan grandes por lo que
se podría estimar un rango del coeficiente de cizalle para todo el año.
Finalmente, se tiene que las diferencias entre ambos terrenos son notables. Los factores
de rugosidad y cizalle de Chagres son bastante mayores que los coeficientes de el otro
terreno. Esto señala que diferencia de velocidad varía notablemente con la altura y es
fácil de comprobar observando las curvas de velocidad respecto del transcurso de un día,
en donde en Chagres se puede apreciar una diferencia notable a pesar de la poca altura
que diferencia a los anemómetros
Con ambos resultados anteriores para cada mes, se usaron los métodos de extrapolación
explicados anteriormente para estimar las velocidades que se tendrán a velocidades
superiores de las medidas por los anemómetros. Estas extrapolaciones están diseñadas,
según la literatura, para pronosticar solamente a 1/3 más de altura de la medida en ese
caso, en este caso se podría llegar a 24 m de altura pero, para efectos de comparación
con otros estudios eólicos aproximó a 100 m de altura dejando abierta la posibilidad de
aumento de errores en la estimación.
En la figura 20 se muestra una gráfica que contempla las proyecciones de velocidad del
viento según la altura usando ambas reglas:
92
5
Velocidad de viento promedio [m/s]
4.5
Proyección de Velocidad de viento
4
3.5
3
2.5
Vel LLP [m/s]
2
Vel PLP [m/s]
1.5
1
0.5
0
2 10 20 30 40 50 60 70 80 90 100 110
Altura sobre el suelo [m]
Figura 20: Proyección de las velocidades del viento promedio según LLP y PLP.
Al comparar estas curvas con las de otros estudios se aprecia que éstas tienen una
pendiente bastante pronunciada. Esto se debe a los altos valores que obtenidos tanto para
el coeficiente de rugosidad como para el coeficiente de cizallamiento. Esta curva es muy
importante porque señala que, si realizadas las mediciones, los vientos obtenidos son
muy pequeños, un aumento en la altura de la torre provocará significativos aumentos en
la velocidad del viento, entregando la opción en caso de no recorrer por completo el
rango de velocidades requerido, de poder aumentar la altura de la torre para llegar a
cubrir esos rangos asumiendo el costo asociado a levantarla.
razón de la magnitud de la diferencia entre los vientos en la noche, mañana y tarde del
terreno y, de esa manera, poder planificar los tiempos para los cambios de configuración
y para la toma de mediciones. En la figura 21 se expone un gráfico que muestra las
variaciones del perfil de velocidad según las horas de día.
3.5
Perfil de velocidad diario
3
Velocidad del viento [m/s]
2.5
2
Vel 6m [m/s]
1.5
Vel 12m [m/s]
1 Vel 18m [m/s]
0.5
0
9
0
1
2
3
4
5
6
7
8
10
11
12
13
14
15
16
17
18
19
20
21
22
23
Horas del día [hr]
Figura 21: Perfil de velocidad diario para el mes de Agosto del 2011.
Aquí se aprecia fácilmente una diferencia de 3 veces en el promedio entre los vientos de
la mañana en comparación con lo vientos vespertinos. Esta irregularidad es bastante
peculiar de este terreno, y ayuda en el sentido de poder realizar mantenciones al equipo
pero disminuye la exposición de horas a vientos de magnitud elevada.
Tabla 15: Resultados de coeficientes de forma y escala según meses del estudio
PERÍO Ener Febrer Marzo Abril Mayo Junio Julio Agost Septiemb Octubr
DO o '11 o '11 '11 '11 '11 '11 '11 o '11 re '11 e '11
WEIB
ULL K 0.31 0.30 0.27 0.24 0.28 0.23 0.22 0.24 0.32 0.29
WEIB
ULL C 15.80 13.30 10.63 5.65 7.28 3.44 3.28 4.39 13.72 13.01
En estos resultados los coeficientes obtenidos son bastante distintos a los coeficientes
típicos de un terreno eólico. En un terreno tipo utilizado para prospección eólica se
tienen parámetros del orden de K=1.2 y C=5.8 lo que le da una forma de cima de baja
altura cercana al x=4 y una caída lenta que significa un aumento de las frecuencias
superiores.
En el caso de Chagres existe un valor alto para el parámetro de forma lo que hace que la
cima de la curva tienda a ubicarse cercana a un promedio de 15 lo que sería excelente
para un terreno eólico, pero el parámetro de escala, al ser tan pequeño, hace que la cima
se mueva a valores muy cercanos de los 0 m/s y que la curva vaya cayendo lentamente
hacia al final permitiendo la existencia de velocidades superiores.
Estos resultados reflejan que el terreno presenta muchos tiempos con velocidades
pequeñas pero con presencia de velocidades altas lo que lo hace bueno para el desarrollo
de una turbina.
0.16
Comparación curvas Weibull
0.14
0.12 Weibull Ch Enero '11
(15.8 y 0.31) h=12m
0.1
Weibull Ch Julio '11 (3.28
0.08 y 0.22) h=12m
0.06 Weibull X Enero'08 (1.701
y 6.274) h=40m
0.04
Weibull X Julio '07 (1.881
0.02 y 7.08) h=40m
0
10.5
11.5
12.5
13.5
14.5
0.5
1.5
2.5
3.5
4.5
5.5
6.5
7.5
8.5
9.5
Se puede mencionar que el terreno Weibull ideal tendría los siguientes parámetros: K=
infinito y C = Velocidad nominal de la turbina.
Para el sitio las velocidades de operación nominal son realmente difíciles de obtener
debido a la baja altura en que se están tomando las mediciones y se encuentra el rotor de
la turbina. Se logra llegar velocidades de 9.2 m/s a altura de rotor en mediciones de
intervalos de 10 minutos con las cuales se puede obtener una primera aproximación
bastante completa de la curva de potencia. Por lo tanto, en el terreno no es posible
realizar la caracterización completa de la turbina a la altura de trabajo que se definió.
En la figura 23 se visualiza cómo se comportan las sondas en un día que presentó una
variada cantidad de velocidades de viento para este terreno:
96
Figura 23: Velocidad del viento medida por las distintas sondas.
En la figura se aprecia que el perfil de las velocidades tiene una correlación muy alta, es
decir todos los perfiles tienen alzas y bajas de velocidad al mismo tiempo siempre
manteniendo una diferencia constante entre sus valores que depende de la velocidad del
viento que se tenga.
1.2
Factor lineal en Velocidad a 12m
1
Factor lineal [m]
0.8
0.6
0.4 FL
0.2
10.5
11.5
12.5
1
9
1.5
10
11
12
13
0.5
2.5
3.5
4.5
5.5
6.5
7.5
8.5
9.5
Velocidad del viento [m/s]
1.6
Coeficiente de Cizalle vs Velocidad del Viento
Coeficiente de cizalle[-]
1.4
1.2
1
0.8
0.6
Alpha 12m
0.4
0.2
0
10.5
11.5
12.5
0.5
1.5
2.5
3.5
4.5
5.5
6.5
7.5
8.5
9.5
12
Coeficiente de Rugosidad en Velocidad 12m
Coeficiente de Rugosidad [m]
10
CR
4
10.5
11.5
12.5
1.5
1
11
0.5
2.5
3.5
4.5
5.5
6.5
7.5
8.5
9.5
10
12
13
Velocidad del viento [m/s]
Se decide tomar en cuenta para el análisis solamente los datos de velocidad mayores
para los que el coeficiente de cada aproximación se estabiliza.
1) Estimación lineal
Se Necesita comprobar qué tan buena es la estimación lineal para predecir la velocidad
central (12 m) desde cada una de las extremas (6 y 18 m). Para eso se tomó 33.000
datos, y se les aplicó el filtro de velocidades inferiores a 3 m/s que entrega el análisis
anterior y se realiza los procedimientos expuestos en la metodología.
99
Los resultados de ambos cálculos se pueden ver en la tabla 16 que caracteriza la fórmula
antes mostrada.
ESTIMACIÓN LINEAL
V6 → V12 V18 → V12 (V6+V18)/2 → V12 V6 , V18 → V12
Coeficiente Coeficiente Coeficiente Coeficiente
1.314 0.825 1.021 1.313 , 0.824
MECM MECM MECM MECM
0.412 0.667 0.422 0.373
Se puede apreciar que las estimaciones llegan a un resultado muy similar siendo ambas
por separado la que estima de mejor manera la velocidad a los 12 m teniendo un error
cuadrático medio de 0.373.
Estos resultados deben ser comparados con los anteriores para poder definir cuál de estas
dos técnicas aproxima de mejor manera la estimación para esta configuración del
aerogenerador, tomando en consideración de lo que se quiere realizar es una
interpolación en una corta distancia.
Para eso se comparó los resultados obtenidos en las distintas aproximaciones y ver cuál
de ellas presenta menores niveles de media cuadrática del error de la aproximación.
Según los resultados se pudo ver que todas las estimaciones se asemejan igual a la
estimación de velocidades a 12 m con un ECM de 0.373. Esta semejanza es tanto para
los datos de 6 a 12 m como para los datos de 18 a 12 m y es mejor en su media
cuadrática, por lo que indiscutiblemente se eligió la estimación lineal por su simplicidad
101
para programar. Esta similitud es debido a la corta distancia que tienen los álabes y a la
gran cantidad de anemómetros para esa distancia.
Este resultado es una base para aproximar una velocidad del perfil de viento que
dependa de la velocidad medida en los extremos y de la velocidad medida a la altura de
rotor. Reemplazando los resultados obtenidos en la curva señalada en la metodología
obtuvo la siguiente velocidad equivalente para el estudio.
| |
Luego se calculó los mejores estimadores de cada anemómetro para todo el perfil de
viento. Para eso se debió de sumar los resultados de aproximación hacia los 12 m con los
resultados de aproximación a 6 y 18 m. La tabla 19 muestra los resultados:
102
Con estos puntos se pueden obtener curvas que describan el comportamiento para todo
el perfil. El resultado de la aproximación entrega curvas lineales con sobre el 0.98
( *
( *
( ( * *
(
( | |
( * ( * ( ) ( ) ( ) ( )
[ ]
| |
Esta curva es la mejor estimación de la velocidad del viento para cada una de las alturas
que tiene el perfil y debe ser integrada y dividida en la altura del perfil que en este caso
es 12 m para poder obtener la velocidad representativa del perfil.
La integral multiplicada por el ancho (a) de la turbina representa el caudal (Q) de aire
que circula por la sección de la turbina. Esta dependencia se puede ver en la siguiente
fórmula:
∫ ( [23]
Al resolver y evaluar la integral [23] para todo el recorrido de la altura del perfil de
viento obtenemos la siguiente expresión:
∫ (
Ésta es el área bajo la curva del perfil entre los 6 y 18 m que es donde circula la
corriente que pasa a través del aerogenerador, y debe ser dividida por la altura para
obtener la velocidad representativa del perfil. Quedando representada por la siguiente
ecuación:
[25]
Para este caso se tienen los registros de una estación meteorológica instalada a 100 m del
aerogenerador la cual tiene 5 años de datos de presión y temperatura para el
emplazamiento, que son datos medidos durante bastante tiempo y cercanos al lugar de
trabajo.
De estos datos se obtiene que la temperatura y la presión tienen los perfiles mostrados en
el anexo L, tanto para invierno como verano; por lo tanto es posible obtener que la
presión se mantiene aproximadamente en 964 milibares para la hora en que existe una
mayor velocidad del viento, las 16:00 horas, mientras que la temperatura en invierno
tiene un valor aproximado de 15 º C y para el verano un valor aproximado de 21 ºC.
En el anexo M se muestra una tabla en donde se calcula la densidad del terreno mediante
la ley general de los gases y entrega valores de para el invierno y de
para el verano, los que son valores bastante parecidos entre ellos.
Estos valores de densidad están dentro de los rangos establecidos para realizar la
correcta obtención de la curva de potencia según las recomendaciones de Measnet 35, las
cuales indican que el valor de la densidad debe mantenerse dentro de un rango de
.
35
Measnet (2009). Measnet Procedure: Power Performance measurement procedure,
Version 5, December.
105
En este caso, se proyecta el plano rotor sobre cualquier plano que tenga contenido el eje
cardánico, entregará un rectángulo el cual va variando sus dimensiones según las
distintas configuraciones que se le den al aerogenerador. Ya se vio que este prototipo es
adaptable, por lo que el diámetro y el largo de patagios pueden variar cambiando el área
de superficie del viento.
A continuación están las fórmulas con las que se obtiene la energía y potencia del
viento, las que dependen de los parámetros ya encontrados y descritos para este caso:
( ( (
Una vez definidos ambos y obtenidos los resultados de la potencia eléctrica generada por
la máquina se obtendrá su eficiencia eólica.
eje del generador. Existen tres puntos de medición de las revoluciones por minuto: el
primer lugar corresponde al rotor superior; el segundo, es la salida de la caja
multiplicadora colocada en el rotor; y el tercero, es la medición realizada por el drive
ACS 800. Esta medición solamente se realiza cuando la máquina está encendida y por
ende ejerciendo el torque de magnetización de los imanes.
En esta parte es fundamental la medición de revoluciones del eje central el cual debe ser
agregado a los análisis diarios y por configuración realizados anteriormente. Los análisis
mensuales no se utilizarán mucho en esta sección, pues por lo general dentro de un mes
se prueba más de una configuración del equipo.
El ánimo para esta sección surgió, primero, para conocer la cantidad de revoluciones que
habrán en el motor y así mantenerlo lo más cercano a su velocidad nominal, pero sin
arriesgar a sobrepasar su velocidad angular máxima; y, segundo, para clasificar a qué
tipo pertenece el aerogenerador según los principios aerodinámicos de funcionamiento
de los álabes.
Los primeros resultados son sin el generador conectado al eje transmisor, puesto que con
esto medirán las capacidades de giro del rotor respecto al viento sin tener en cuenta los
roces asociados al giro del motor que, multiplicados en torque por las cajas
multiplicadoras, se hacen más significativos.
14
Velocidad, rpm y Potencia
12 Curva de Rpm [1min]
10
8
6 RPM Rotor [rpm]
4
V6m [m/s]
2
V12m[m/s]
0
1:49:10 PM
2:58:10 PM
4:07:10 PM
5:16:10 PM
6:25:10 PM
7:34:09 PM
8:44:08 PM
9:53:08 PM
11:31:10 AM
1:09:08 AM
2:19:08 AM
3:28:10 AM
4:37:09 AM
5:46:10 AM
6:55:10 AM
8:04:10 AM
9:13:10 AM
12:00:08 AM
10:22:10 AM
12:40:10 PM
11:02:08 PM
V18m[m/s]
12
Velocidad vs Tiempo (Promedio 10 minutos)
Velocidad del viento [m/s]
10
8
6
V6m
4
2 V12m
0 V18m
1:00:00 AM
8:00:00 PM
2:00:00 AM
3:00:00 AM
4:00:00 AM
5:00:00 AM
6:00:00 AM
7:00:00 AM
8:00:00 AM
9:00:00 AM
1:00:00 PM
2:00:00 PM
3:00:00 PM
4:00:00 PM
5:00:00 PM
6:00:00 PM
7:00:00 PM
9:00:00 PM
12:00:00 AM
12:00:00 PM
10:00:00 PM
11:00:00 PM
10:00:00 AM
11:00:00 AM
Rpm
gráfico deberán ser comparados con los resultados obtenidos de toda la configuración y
así poder comprobar el correcto análisis de los datos.
Este gráfico es una inspección rápida para ver el funcionamiento del rotor y de la
transmisión mecánica del aerogenerador y así poder descubrir de manera rápida
cualquier anomalía en las cajas cardanes y/o crucetas del sistema de transmisión.
Además se puede ver la rápida reacción que tiene el rotor respecto de los cambios de
velocidad de viento, pues ambas curvas se comportan de manera muy similar sin
presentar mayor cantidad de retrasos.
La figura 29 muestra las revoluciones por minuto obtenidas para la prueba del
aerogenerador en vacío.
4.5
Frecuencia de rpm cada 1 minuto
4
3.5
Frecuencia de rpm [1min]
2.5
2
FrecRpm
1.5
0.5
0
11.5
10.5
12.5
13.5
14.5
15.5
16.5
17.5
0.5
1.5
2.5
3.5
4.5
5.5
6.5
7.5
8.5
9.5
que la velocidad nominal. Es decir puede haber pérdida de eficiencia pero no hay riesgo
de sobrepasar las revoluciones máximas del generador.
En la figura 30 se muestran los resultados obtenidos para el TSR del sistema, cuando el
motor eléctrico está desconectado de la transmisión mecánica, lo que permite disminuir
el roce del giro del rotor y por ende aumentar el nivel de revoluciones medidas al
máximo obtenible para cada velocidad.
14.0
12.0
10.0
8.0
6.0 RPM Rotor [rpm] 1min
4.0 RPM Rotor [rpm] 10min
2.0
0.0
10.3
3.2
0.3
1.2
2.2
4.2
5.2
6.2
7.2
8.3
9.3
Figura 30: Velocidad de rotor en aerogenerador respecto del viento sin carga.
Es posible apreciar que la velocidad del viento es directamente proporcional a las
revoluciones del equipo y que la curva de un minuto se comporta muy similar a la de 10
minutos en un principio para luego separarse debido a que la de 1 minuto aumenta el
rango de datos utilizado.
En la figura 31 se muestra el gráfico del TSR según la velocidad equivalente del viento
para el registro en intervalos de 1 y 10 minutos.
110
1.2
TSR para 1 y 10 minutos
1.0
0.8
TSR [-]
0.6
TSR 10min
0.4
TSR 1min
0.2
0.0
10.3
10.8
3.7
0.3
0.7
1.2
1.7
2.2
2.7
3.2
4.2
4.7
5.2
5.7
6.2
6.7
7.2
7.7
8.3
8.8
9.3
9.7
Velocidad Equivalente de viento [m/s]
La primera curva de torque que se suministró al drive fue una curva que se obtuvo a
partir de la prueba inicial que consiste utilizar el generador como motor y ver el torque
que se necesita para llegar a las 20 RPM, el cual es un valor un poco mayor al obtenido
en la curva de revoluciones de los resultados anteriores. Se utilizó 20 rpm porque esa es
la velocidad nominal del motor vista desde el rotor con las cajas multiplicadoras
elegidas. Se definió un cut-in de 10% de la velocidad nominal del motor. El obtener esta
curva en un lugar sin viento es un buen punto de partida para comenzar a ver cómo se
debiese comportar la curva de torque del generador.
Se obtuvieron los puntos de torque y revoluciones por minuto para distintas potencias
suministradas en el motor eléctrico y la unión de los puntos formó una curva que puede
ser descrita mediante una curva de tendencia cuadrática la cual se muestra a
continuación:
( [28]
La forma de esta curva corrobora que el aerogenerador tiende a funcionar como una
bomba o un ventilador en donde la curva de torque es cuadrática según la velocidad del
fluido que se esté impulsando. Luego en la figura 32 se muestra la gráfica la función
descrita anteriormente:
Torque vs rpm
25
Torque [%deTorqueNominal]
20
15
10
Torque
5
0
0 10 20 30 40 50 60 70 80 90 100
Revoluciones del motor [%deRpmNominal]
Se aprecia que para las 20 revoluciones por minuto del rotor se alcanza un torque del
23% del torque nominal, lo que corresponde a 462 Nm. Ambos datos indican que se
necesitó una potencia de 11 kW para hacer girar el generador a esa velocidad. Esta
potencia representa la potencia que ejerce la fuerza de roce para cada velocidad de rotor.
Figura 33: Interfaz del programa adaptativo con el programa de la curva inicial.
Este primer programa trabaja controlando la función de Direct Torque Control del drive
incluye la programación de la velocidad cut-in, en este caso 10 % de la velocidad
nominal y la curva de torque respecto de las revoluciones por minuto como función
cuadrática, que en este caso representa a la curva cuadrática del gráfico anterior.
Los resultados eléctricos de esta curva de torque serán mostrados en la sección siguiente
sobre la metodología de obtención de curva de potencia.
113
Estas curvas son el dato más importante que caracteriza a una turbina eólica, pues con
estas curvas pudo predecir la energía que generará según las distintas características de
viento que tenga un terreno, o sea, la información entregada por las mediciones de
frecuencia, turbulencias, velocidad y dirección, entre otras, de un terreno más la curva de
potencia son las que hacen posible la evaluación de un proyecto eólico. Es
responsabilidad de cada propietario entregar esta curva la cual debe ser obtenida
mediante los procedimientos usados en la norma de generación de potencia eólica.
Este análisis consiste en caracterizar las variables respecto del tiempo para verificar
errores y estimar cómo son las variaciones durante el período de registro de datos en
plazos de registro de 1 y 10 minutos.
Las siguientes curvas de potencia son obtenidas a partir de la configuración 18, la cual
consiste en la torre con 12 patagios de cuerda 0.8mts, 2 patagios de cuerda 1.6 m y de
largo 12 m, dispuestos de manera equidistante unos de otros a un diámetro de 12.2 m.
A continuación se presenta la figura 34 que presenta uno de los análisis diarios que
deben hacerse utilizando todos los parámetros necesarios para el cálculo de la potencia
del aerogenerador.
114
12
Potencia vs Tiempo
10
Velocidad, rpm y Potencia
6 V6m [m/s]
V12m [m/s]
4
V18m [m/s]
2
Rpm [rpm]
0 Pot [kW]
1:21:37 PM
8:01:37 PM
4:02:03 PM
6:44:08 PM
2:45:10 PM
9:25:10 PM
5:26:04 PM
1:28:10 PM
8:09:08 PM
6:40:05 AM
2:41:07 AM
9:21:58 AM
5:23:07 AM
1:24:09 AM
8:04:09 AM
4:06:02 AM
6:48:09 AM
12:00:04 AM
10:46:03 AM
12:07:08 AM
10:43:06 PM
-2 12:03:07 PM
Los datos de la calidad de la corriente eléctrica no son analizados todavía, hasta que
podamos obtener la curva de potencia cercana al máximo.
contenedores
Al aplicar nuevamente este método a los datos incluyendo los valores de la potencia
obtenida para formular una tabla que especifique todas las características señaladas en la
metodología, esta tabla servirá para optimizar la potencia del aerogenerador. Luego con
los datos contenidos en la tabla podremos graficar una curva que muestre la potencia
115
generada en función de la velocidad equivalente del viento, que será el resultado final de
la tesis.
La primera curva de potencia probada fue la curva inicial y está en base a la curva de
torque obtenida mediante la primera inicialización. Esta curva, como era de esperar,
produjo resultados muy bajos de potencia debido al poco torque que presentaba, por lo
que no será mostrada como resultado.
Luego se decidió aumentar el valor de la curva de torque a un valor que produjera más
potencia, pero que mantuviera la tendencia cuadrática de ésta. La segunda curva
programada es la mostrada a continuación:
El resultado entregado por esta curva se muestra en la tabla 20 la cual contempla todos
los datos importantes para la potencia clasificados según intervalos de 0.5 m/s de
velocidad de viento. Este resultado es el obtenido para la configuración descrita en la
primera parte de la tesis y posibilita la futura optimización del prototipo según la
metodología provista.
Los resultados se muestran para para intervalos de medición de 10 minutos (Tabla 20) y
1 minuto (Tabla 21), de esa manera podremos comprobar la precisión de cada uno de
ellos para predecir la energía generada durante el intervalo de medición.
En este caso, solamente presentamos los resultados finales según la Tabla B, pues nunca
podremos llegamos a las velocidades de cut -off debido a las condiciones del terreno.
116
Standard Standard
Hub Height Wind
RPM Rotor [rpm] Potencia [kW] Turbulence deviation of Wind Shear deviation of
Bin-no. [-] Speed [m/s] Cp [-] 10min Lamda [-] 10min No of data sets
10min 10min Uncertainty [-] Turbulence Exponent [-] Wind Shear
10min
Uncertainty [-] Exponent [-]
Hub Height
RPM Rotor [rpm] Potencia Lamda [-]
Bin-no. [-] Wind Speed Cp [-] 10min No of data sets
1min [kW]1min 10min
[m/s] 1min
Se puede apreciar la tabla cada 10 minutos contiene toda la información sobre la curva
de potencia mientras que la tabla de 1 minuto le faltan algunos cálculos debido a faltas
en la capacidad de registro de los anemómetros. Por lo que se tabula sólo la información
posible de obtener.
La potencia registrada presenta como parte de ella el consumo del drive lo que reduce la
potencia final obtenida desde la turbina. Se aprecia que el coeficiente lambda, para cada
uno de los intervalos de velocidad de viento, se mantiene con un pequeño aumento a
medida que crece la velocidad del viento.
Es apreciable que en la medición de mayor velocidad de viento tiene menores datos de
frecuencia por lo que presenta más errores en la medición, por este motivo es muy
importante que se lleguen a velocidades por sobre la velocidad nominal.
118
Para poder ver más claramente la curva de potencia de la turbina es que de esta tabla se
muestran las figuras 35 y 36 en donde se grafica la potencia versus la velocidad del
viento para 10 y 1 minuto respectivamente.
1.200
1.000
Curva de Potencia [10min]
0.800
0.600
0.400
Potencia
0.200
Potencia [kW] 10min
0.000
0.2
1.2
2.2
3.3
4.2
5.2
6.3
7.2
8.2
-0.200 9.1
-0.400
-0.600
-0.800
Velocidad a la altura del eje [m/s]
0.40
0.20 Potencia [kW]1min
0.00
0.26
1.19
2.21
3.24
4.22
5.23
6.24
7.24
8.23
9.21
10.16
-0.20
-0.40
-0.60
-0.80
Velocidad a la altura del eje [m/s]
aerogenerador el que debe alimentar este consumo mientras está girando. Este consumo
es significativo para la potencia que está generando.
Es posible ver que existe poca diferencia entre el análisis de un minuto sobre el análisis
cada 10 minutos el primero llega hasta los 10.1 m/s mientras que el segundo lo hace sólo
hasta los 9.1 m/s lo que es menor que la velocidad nominal que debería tener el
aerogenerador. Los mejores resultados están presentes en el análisis de 10 minutos
debido a que se llega a una menor velocidad de viento para los mismos datos de potencia
producida.
Es importante obtener los resultados para el coeficiente de potencia y así darse cuenta de
la eficiencia que está teniendo la turbina. Para esto se tomaron en cuenta los valores del
coeficiente de potencia y TSR de la tabla anterior, fueron ordenados en función del TSR
y graficados para ver cómo se comporta el coeficiente.
0.025
0.020
0.015
0.010 Cp [-] 10min
0.005
0.000
0.11 0.16 0.20 0.23 0.26 0.27 0.28 0.30 0.32 0.32 0.33 0.33
Tip Speed Ratio
Debido a que la curva en todo el rango medido es creciente, aún no se ha sido capaz de
encontrar el valor del TSR que genere el máximo de potencia, por lo que para las
siguientes iteraciones de curva de potencia se deberán programar curvas de torque
menores para poder conseguir un mayor rango de TSR.
Los analizadores de redes trifásicas que fueron instalados cuentan con medidores de
energía integrados los que constan de variables acumulativas. Tomando el valor inicial y
el final se puede verificar la energía extraída y la consumida por el aerogenerador.
6. Análisis y Conclusiones
el período de medición y que haya menor pérdida de datos debido a los filtros aplicados
a los datos.
económico censar la velocidad del viento en todo su perfil. Además, como estas turbinas
están diseñadas para ser ubicadas en sitios urbanos, en donde hay mayor turbulencia y
coeficiente de cizalle, el método de la norma actual presentaría errores mayores para
estimar la velocidad representativa.
importante contar con una norma que certifique la igualdad de estos intervalos, tanto
para sensar los terrenos como para obtener la curva de potencia, que mejore la
estimación de la energía y que ayude a comparar la magnitud de las distintas curvas
mostradas de los aerogeneradores del mercado.
Al comparar la eficiencia del equipo con la de otros modelos similares se podrá dar
cuenta que es bastante menor. Esto puede ser atribuido a la geometría de los álabes y al
hecho de que varios de los componentes de transmisión y generación de energía eléctrica
se encuentran sobredimensionados con respecto a las dimensiones del rotor. Esto
disminuiría el desempeño del aerogenerador, pues la potencia mecánica disponible se
consumiría en accionar aquellos mecanismos sobredimensionados como lo son
principalmente el rotor, las cajas multiplicadoras y el generador eléctrico. Lo que daría
resultados menores a lo esperado.
Se esperaba que el diseño del álabe fuera impulsado por el efecto sustentación que es el
que ha demostrado mejor eficiencia pero, al parecer, la hendidura que tiene el perfil
patagio provoca una caída en pérdida aerodinámica continuada (“stall”) durante los
flujos de trabajo que no le permiten entrar en sustentación manteniéndolo
constantemente en bajo régimen de arrastre. Esto es fácil de notar debido a que
difícilmente se superan los valores de TSR = 1.1, incluso para las pruebas realizadas en
vacío.
126
El método propuesto en esta tesis describe bien la energía generada por la turbina, que es
el objetivo primordial de la tesis, obteniendo errores del 3.1% para el método de los
intervalos aplicado según el período más corto de medición de datos (1 minuto). Sin
embargo sería bueno probar el método en una turbina comercial y comparar los
resultados obtenidos con los suministrados por la empresa proveedora de la turbina.
BIBLIOGRAFÍA
Ropatec AG (2003). HAWT versus VAWT Small VAWTs find a clear niche,
Renewable Energy Focus, 1471 0846/03, 44-46.
Ropatec AG (2011). Features Ropatec´s small wind turbine vertical axis Recuperado el
10 de Octubre del 2011, de http://www.ropatec.com
Mahyoub H.Al Buhairi (2006). A statistical analysis of wind speed data and assessment
of wind energy potential in Taiz-Yemen, Physics Department Taiz Univeristy. Republic
of Yemen. Recuperado de Assiut University: http://www.aun.edu.eg
J.F. Manwell; J.G. McGowan y A.L. Rogers (2002). Wind Energy Explained Theory,
Design and Application. University of Massachusetts, Amherst, USA: [Wiley]
128
Ragheb M. (2011). Optimal Tip Speed Ratio. Recuperado el 4 de Mayo del 2011 de la
University of illinois at urbana-champaign http://www.illinois.edu
Vásquez Leiva M. (2007). Diseño de una turbina eólica de eje vertical con rotor mixto.
(Tesis de maestría no publicada). Universidad Austral de Chile, Valdivia, Chile.
ANEXOS
130
Esfuerzos admisibles
Las constantes permitidas para esta estructura fueron calculadas de la siguiente manera
Para esta simulación se asignó una carga de magnitud igual al peso soportado por el
brazo. Para determinarla se distribuyó de igual manera la carga soportada por el brazo
superior e inferior, la que corresponde al peso de dos patagios. Por lo tanto la carga
correspondiente a cada brazo (superior e inferior) corresponde a 1859,99 N (flecha
morada, figura 39). El empotramiento indicado en la figura por flechas verdes, fue fijado
en el extremo sujeto al rotor y un segundo empotramiento simulando el tensor de apoyo.
Por último se representa la gravedad con una flecha color rojo, y se estimó con un valor
de 9,81 m/s2.
132
Resultado de la Modelación
Figura 41: Representación del factor de seguridad para el análisis del brazo inferior
135
Con la intención de obtener mayores revoluciones por minuto se efectuó durante las
pruebas de validación, los cambios de cajas multiplicadores, mejorando la relación de
estas y midiendo su resultado posteriormente. A continuación se muestran fotografías de
cada una de ellas y se detallan.
Corona Fija, Porta satélites con traslación, satélites en rotación, sol en rotación.
137
Diámetro exterior: 88 mm
Diámetro interior: 66 mm
Corona Fija, Porta satélites con traslación, satélites en rotación, sol en rotación.
139
Easy Vertical HE
Wind Speed Bin (m/s) Swept Area [m2] VAWT Power [VA] Wind Power [VA] Efficiency
4 3.42 27 134 20.1
5 3.42 53 261 20.1
6 3.42 91 451 20.1
7 3.42 144 716 20.1
8 3.42 217 1068 20.3
9 3.42 318 1521 20.9
10 3.42 447 2086 21.4
11 3.42 610 2777 22.0
12 3.42 911 3605 25.3
13 3.42 1155 4583 25.2
14 3.42 1200 5725 21.0
15 3.42 1200 7041 17.0
16 3.42 1200 8545 14.0
17 3.42 1200 10250 11.7
18 3.42 1200 12167 9.9
Simply vertical
Wind Speed Bin (m/s) Swept Area [m2] VAWT Power [VA] Wind Power [VA] Efficiency
4 6.6 64 258 24.8
5 6.6 127 503 25.2
6 6.6 224 870 25.8
7 6.6 362 1381 26.2
8 6.6 551 2061 26.7
9 6.6 799 2935 27.2
10 6.6 1117 4026 27.7
11 6.6 1497 5359 27.9
12 6.6 1957 6957 28.1
13 6.6 2506 8845 28.3
14 6.6 3152 11047 28.5
15 6.6 3600 13588 26.5
16 6.6 3600 16490 21.8
17 6.6 3600 19780 18.2
18 6.6 3600 23480 15.3
141
Maxi Vertical
Wind Speed Bin (m/s) Swept Area [m2] VAWT Power [VA] Wind Power [VA] Efficiency
4 11.75 112 459 24.4
5 11.75 224 896 25.0
6 11.75 394 1548 25.4
7 11.75 638 2458 26.0
8 11.75 970 3670 26.4
9 11.75 1407 5225 26.9
10 11.75 1965 7168 27.4
11 11.75 2634 9540 27.6
12 11.75 3445 12385 27.8
13 11.75 4410 15747 28.0
14 11.75 5547 19668 28.2
15 11.75 6977 24190 28.8
16 11.75 8776 29358 29.9
17 11.75 10000 35214 28.4
18 11.75 10000 41801 23.9
Big Star
Wind Speed Bin (m/s) Swept Area [m2] VAWT Power [VA] Wind Power [VA] Efficiency
4 34.4 241 1343 17.9
5 34.4 500 2623 19.1
6 34.4 1005 4533 22.2
7 34.4 1734 7198 24.1
8 34.4 2796 10744 26.0
9 34.4 4276 15297 28.0
10 34.4 6068 20984 28.9
11 34.4 8346 27930 29.9
12 34.4 11360 36260 31.3
13 34.4 14894 46102 32.3
14 34.4 19166 57580 33.3
15 34.4 19900 70821 28.1
16 34.4 19900 85950 23.2
17 34.4 19900 103094 19.3
18 34.4 19900 122379 16.3
142
A continuación se expone un gráfico que señala la diferencia existente ente tomar los
datos como promedios de 1 minuto a graficarlos como promedios de 10 minutos, para el
caso del segundo el gráfico disminuye notablemente la magnitud de los valores. Las
mediciones fueron realizadas para el mismo intervalo de tiempo.
10.0
9.0
Velocidad del viento (Promedio 1 minuto)
8.0
Velocidad [m/s]
7.0
6.0
5.0
4.0
3.0 V6m
2.0
1.0 V12m
0.0
V18m
11:00:00 AM
1:00:00 PM
2:00:00 PM
3:00:00 PM
4:00:00 PM
5:00:00 PM
6:00:00 PM
7:00:00 PM
8:00:00 PM
9:00:00 PM
1:00:00 AM
2:00:00 AM
3:00:00 AM
4:00:00 AM
5:00:00 AM
6:00:00 AM
7:00:00 AM
8:00:00 AM
9:00:00 AM
12:00:00 AM
10:00:00 AM
12:00:00 PM
10:00:00 PM
11:00:00 PM
Tiempo [1min]
8
7
6
5
4
3 V6m
2
1 V12m
0 V18m
4:00:00 PM
1:00:00 PM
2:00:00 PM
3:00:00 PM
5:00:00 PM
6:00:00 PM
7:00:00 PM
8:00:00 PM
9:00:00 PM
1:00:00 AM
2:00:00 AM
3:00:00 AM
4:00:00 AM
5:00:00 AM
6:00:00 AM
7:00:00 AM
8:00:00 AM
9:00:00 AM
11:00:00 AM
12:00:00 PM
10:00:00 PM
11:00:00 PM
12:00:00 AM
10:00:00 AM
3
Turbulencias de Chagres vs LV2
Intensidad de turbulencias [-]
10.5
11.5
12.5
13.5
14.5
15.5
7.5
0.5
1.5
2.5
3.5
4.5
5.5
6.5
8.5
9.5
Turb CH 18m
Hora
Día 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 Prom Max
1 10 10 9 9 9 9 9 9 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 9.9 10.5
2 10 10 10 10 10 10 10 10 10 11 11 13 14 15 16 16 16 15 13 12 11 10 10 11 11.8 16.4
3 10 8 9 8 7 6 7 7 7 9 11 14 16 18 19 19 18 17 15 14 12 11 11 9 11.8 18.9
4 7 6 6 5 6 6 8 7 7 8 8 9 10 10 10 10 11 11 10 8 7 6 6 7 7.9 10.7
5 6 6 6 5 4 4 3 4 3 6 9 12 14 16 16 17 16 14 12 11 10 9 7 6 9.0 16.6
6 6 4 5 4 3 3 2 2 3 6 9 12 13 14 15 16 16 14 13 11 10 9 8 7 8.5 15.8
7 6 6 5 6 5 5 5 5 5 8 11 13 15 17 19 19 19 18 14 12 11 9 8 8 10.3 18.8
8 8 6 8 7 7 7 6 5 7 9 12 15 18 20 21 21 20 19 17 14 13 10 9 10 12.1 21.1
9 9 8 8 9 7 7 6 8 8 9 12 15 17 19 20 21 21 19 17 14 12 11 10 9 12.3 21.5
10 9 8 8 7 6 6 5 5 5 8 12 15 16 17 16 16 15 14 13 13 12 12 11 11 10.9 16.9
11 10 10 10 10 10 9 10 9 9 10 11 11 11 12 12 14 13 12 11 11 10 8 7 7 10.3 13.8
12 7 7 6 5 6 5 5 4 5 7 10 13 15 17 18 19 19 17 15 12 10 7 7 6 10.1 19.1
13 6 7 7 7 6 5 5 5 5 7 9 11 13 13 14 13 12 12 11 10 9 9 9 9 8.9 13.5
14 9 8 8 8 8 8 8 8 8 9 9 10 10 11 11 11 11 11 11 10 10 9 9 9 9.4 11.3
15 8 8 8 9 8 8 9 9 9 10 11 11 10 9 9 8 8 8 8 8 8 8 8 8 8.8 11.3
16 8 9 9 10 10 10 10 10 10 9 9 10 12 13 13 13 13 12 11 10 9 8 8 7 10.2 13.4
17 6 6 5 6 6 5 4 3 5 7 10 12 13 14 14 14 13 12 11 9 9 9 9 8 8.7 14.0
18 7 7 8 8 8 8 8 7 8 8 9 10 11 11 11 10 9 9 8 8 8 8 7 7 8.4 11.4
19 7 7 7 7 7 7 6 6 6 7 9 10 11 12 12 12 12 11 9 9 8 8 7 6 8.5 12.2
20 5 5 4 4 4 4 2 2 4 8 10 11 13 13 14 14 14 13 11 9 9 8 7 5 8.0 13.8
21 4 4 3 3 2 2 2 1 4 6 9 11 14 15 16 16 15 14 13 12 11 8 6 6 8.2 15.9
22 6 6 5 6 5 5 3 4 5 7 10 13 15 15 15 16 15 13 12 10 9 8 7 7 9.0 15.8
23 6 6 6 6 6 6 6 6 6 7 7 8 9 10 10 11 12 11 9 7 5 5 5 4 7.3 11.8
24 5 4 4 4 2 2 3 1 4 6 10 13 15 18 18 17 16 15 13 10 9 8 7 7 8.8 18.1
25 6 5 4 3 2 1 1 1 3 7 10 13 15 17 18 18 17 15 13 11 10 9 8 8 8.9 18.0
26 8 8 8 8 8 8 8 7 8 10 11 12 13 13 12 12 12 12 11 11 10 10 10 10 10.0 13.4
27 9 9 9 9 9 9 9 9 9 10 10 11 12 12 13 13 13 12 11 9 9 9 9 9 10.1 13.1
28 9 6 6 5 5 4 4 4 5 7 9 11 13 15 16 17 15 14 14 13 12 10 8 8 9.6 16.9
29 7 6 6 5 6 5 6 5 8 11 15 17 18 20 22 23 23 22 18 16 14 11 12 10 12.8 22.8
30 9 8 10 9 8 7 7 8 9 13 16 19 22 23 24 24 23 21 19 18 16 15 12 10 14.6 23.8
31 9 7 7 6 6 5 4 4 7 11 13 15 17 18 18 18 17 16 15 14 13 12 11 10 11.2 17.6
NOTAS1.-
: CASILLAS EN BLANCO = > DATO INVALIDADO POR CALIBRACION O FALLA PROMEDIO DEL MES 9.9
Hora
Día 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 Prom Max
1 962 963 963 963 963 963 963 963 964 964 965 965 965 965 965 964 964 964 965 965 965 964 964 964 964 965
2 964 964 964 964 964 964 964 965 965 966 967 967 968 968 968 967 966 965 964 964 964 964 964 964 965 968
3 963 963 963 963 963 962 963 963 962 963 964 964 964 964 964 963 962 961 960 960 960 959 960 960 962 964
4 960 959 959 959 959 959 960 959 962 962 962 963 964 964 964 964 965 964 963 962 962 962 962 962 962 965
5 963 963 963 962 962 962 961 960 961 963 964 964 965 964 964 963 962 962 961 961 961 961 961 961 962 965
6 961 961 961 961 961 961 962 961 961 963 965 965 965 965 965 964 964 964 964 964 964 963 963 964 963 965
7 964 963 963 962 962 961 963 963 963 964 965 965 966 965 965 964 963 963 962 962 962 961 961 961 963 966
8 962 961 961 961 961 960 960 960 961 961 962 963 965 964 964 963 963 962 962 962 962 962 962 962 962 965
9 962 962 962 962 961 960 961 961 962 963 964 964 965 965 965 964 964 963 962 962 962 962 962 962 962 965
10 962 962 962 962 962 962 962 962 962 963 965 966 965 965 964 963 963 963 963 963 963 963 963 963 963 966
11 962 962 962 962 961 962 962 963 963 964 965 965 966 967 966 967 966 965 965 965 965 965 965 965 964 967
12 965 965 964 964 963 963 963 963 963 964 965 965 966 965 965 964 963 962 962 961 961 961 961 961 963 966
13 960 961 961 960 960 960 960 960 960 961 962 963 964 964 964 963 963 963 963 963 964 964 964 964 962 964
14 964 964 964 964 964 965 965 965 965 966 966 966 966 966 966 965 965 965 964 964 963 963 963 963 965 966
15 963 962 962 962 962 962 962 962 962 962 962 962 962 962 961 961 960 960 960 960 961 961 960 960 961 963
16 960 959 959 959 959 958 958 959 959 960 961 962 963 965 964 964 964 963 962 963 963 963 963 962 961 965
17 962 962 962 962 962 961 962 962 962 963 964 965 965 965 964 964 963 962 962 963 963 963 963 963 963 965
18 963 963 963 963 963 964 964 964 965 965 965 966 966 966 965 965 964 964 963 963 964 964 964 964 964 966
19 964 964 964 964 965 965 965 965 966 966 967 968 968 968 967 968 968 967 967 967 968 968 968 968 966 968
20 967 967 967 967 966 966 966 966 966 967 968 969 969 969 968 968 967 967 967 967 967 967 967 967 967 969
21 966 966 966 966 966 965 966 966 966 966 968 969 969 968 967 966 965 964 964 963 964 964 963 962 966 969
22 962 962 961 960 960 959 959 959 960 960 962 962 961 962 962 962 961 960 960 961 961 960 960 960 961 962
23 961 962 962 961 962 962 962 962 964 964 965 965 965 965 965 965 963 962 963 962 963 963 963 962 963 965
24 962 962 961 961 961 960 960 960 961 962 963 964 964 963 963 962 962 962 962 962 962 962 962 962 962 964
25 962 962 962 961 961 961 961 961 961 962 964 964 964 964 963 963 963 962 961 961 961 961 961 961 962 964
26 962 962 962 961 962 962 962 963 964 965 966 966 966 964 965 965 964 964 964 964 964 964 964 964 964 966
27 964 964 963 963 963 964 964 964 964 965 966 966 967 967 967 967 967 967 967 967 967 968 968 969 966 969
28 969 969 969 969 969 969 968 969 969 970 971 971 970 970 969 968 967 967 966 966 966 965 965 965 968 971
29 964 964 964 963 963 963 963 964 965 966 966 967 967 966 966 965 965 963 962 962 962 962 962 962 964 967
30 961 961 961 961 960 960 960 961 961 962 963 964 964 963 962 961 961 960 960 960 961 961 960 960 961 964
31 960 960 960 959 959 959 960 960 961 962 963 964 964 963 963 963 962 961 961 961 961 961 961 961 961 964
NOTAS1.-
: CASILLAS EN BLANCO = > DATO INVALIDADO POR CALIBRACION O FALLA PROMEDIO DEL MES 963
148
Hora
Día 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 Prom Max
1 14 14 14 13 12 11 11 11 11 12 13 14 15 16 17 17 18 18 17 16 15 14 14 13 14.2 17.7
2 13 12 11 11 10 10 9 9 12 14 16 16 17 17 18 18 18 18 17 16 15 14 14 14 14.2 18.3
3 13 13 13 13 13 13 13 13 13 14 14 15 16 17 18 18 18 18 18 18 17 16 14 14 15.0 18.3
4 13 12 12 11 11 10 9 9 13 16 18 19 21 22 24 24 24 24 22 21 20 18 17 16 16.9 24.4
5 15 14 14 13 12 11 11 12 13 16 18 19 20 22 22 23 22 22 21 19 17 15 14 14 16.7 22.7
6 13 13 12 12 12 12 12 11 12 14 15 16 17 18 20 21 22 22 21 19 17 16 14 14 15.6 21.6
7 13 12 12 12 12 12 12 12 13 14 15 16 17 18 19 20 21 21 20 19 17 17 16 15 15.5 20.8
8 14 13 12 12 12 11 10 11 13 15 17 19 21 22 23 24 24 23 22 20 18 17 16 15 16.8 23.8
9 14 14 13 12 12 9 9 11 13 16 19 21 23 24 26 26 26 25 24 23 21 19 18 18 18.2 26.5
10 17 15 13 12 11 11 10 11 16 20 22 24 26 26 27 27 25 24 24 23 21 19 18 17 19.1 26.7
11 17 16 14 13 12 12 12 11 15 18 19 19 19 18 19 19 19 19 18 17 17 16 15 14 16.1 19.5
12 14 14 13 12 13 12 12 12 13 14 16 18 19 21 23 24 23 21 20 18 17 16 15 16.5 23.8
13 14 14 13 13 12 12 11 11 12 15 17 18 19 21 22 24 24 24 23 22 20 18 17 16 17.2 24.1
14 15 15 14 13 13 12 11 11 14 15 16 18 20 22 22 21 21 20 19 18 17 16 15 14 16.4 21.9
15 14 14 14 14 14 14 14 14 14 14 15 16 16 17 17 18 18 18 17 16 15 15 14 14 15.2 17.8
16 13 13 13 13 13 13 13 13 13 13 14 16 17 18 19 19 18 17 17 16 15 15 15 14 15.1 18.8
17 14 14 12 12 12 13 14 14 15 16 17 20 21 22 23 24 24 23 22 20 19 18 17 16 17.6 24.0
18 16 15 15 14 13 12 11 12 14 17 19 21 23 24 24 23 23 23 23 21 19 17 16 15 17.9 23.7
19 15 14 13 13 12 12 11 11 13 15 17 20 22 23 25 25 24 23 23 21 20 18 17 16 17.6 25.3
20 15 14 14 13 12 12 10 11 14 16 18 20 22 22 23 23 22 22 22 21 19 17 16 15 17.2 23.0
21 15 12 11 11 10 9 9 9 14 18 20 22 23 24 25 26 25 24 23 22 20 18 17 15 17.6 25.5
22 15 14 13 12 11 10 9 10 14 17 19 22 23 24 25 26 25 24 23 22 20 18 17 16 17.9 25.6
23 15 15 14 13 11 10 10 11 15 18 20 22 23 24 24 24 24 24 23 22 20 18 17 16 17.9 24.5
24 15 15 14 13 13 12 11 11 12 13 15 16 18 19 21 22 22 22 21 20 18 17 16 16 16.4 22.2
25 15 15 13 12 11 10 10 10 13 17 19 21 23 23 23 23 22 22 21 20 18 17 16 15 17.0 23.3
26 15 14 13 13 11 11 11 12 12 12 14 15 18 19 20 21 21 20 20 19 17 15 14 13 15.4 20.9
27 12 12 12 11 11 11 11 12 12 12 14 16 20 21 21 21 21 20 19 18 16 15 14 14 15.3 21.4
28 14 13 12 12 11 11 9 10 13 16 19 21 23 25 26 26 26 26 25 24 22 21 19 18 18.4 26.4
29 18 17 15 14 13 12 12 12 16 20 24 24 25 26 26 25 25 24 23 22 21 19 18 18 19.5 26.2
30 17 16 15 13 12 11 11 11 13 17 19 20 22 24 24 24 25 24 24 22 21 20 18 18 18.4 24.5
31 17 16 16 13 12 12 11 11 15 17 20 21 22 22 23 23 23 22 21 20 18 17 16 15 17.7 23.5
NOTAS1.-
: CASILLAS EN BLANCO = > DATO INVALIDADO POR CALIBRACION O FALLA PROMEDIO DEL MES 16.8
Hora
Día 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 Prom Max
1 962 962 961 961 961 962 962 962 963 964 964 964 965 965 964 964 964 962 960 960 962 961 961 962 962 965
2 963 962 962 962 961 961 960 962 964 965 965 965 965 965 965 965 965 965 964 964 964 964 964 964 963 965
3 963 963 963 963 963 963 963 964 964 965 965 965 965 965 965 964 964 963 963 963 962 962 962 962 964 965
4 961 961 961 961 961 961 961 961 962 963 963 963 963 963 963 962 962 962 962 961 961 961 961 961 962 963
5 961 961 961 961 960 960 961 961 963 964 964 964 964 964 964 964 963 963 963 962 962 962 963 963 962 964
6 963 962 962 962 962 962 962 963 964 964 965 965 965 965 964 965 964 964 963 963 962 962 962 962 963 965
7 962 962 961 961 961 961 962 962 962 963 964 964 963 963 963 963 963 963 962 962 962 962 962 962 962 964
8 962 961 961 961 961 961 962 962 964 965 965 965 965 965 965 965 965 964 964 963 963 963 963 963 963 965
9 963 962 962 962 961 961 961 962 963 965 965 965 965 965 964 964 964 964 963 963 962 962 962 962 963 965
10 962 961 960 959 959 958 959 959 960 962 962 962 962 962 961 961 961 960 960 959 959 959 959 959 960 962
11 959 958 958 957 956 957 957 958 959 960 960 960 960 960 961 961 961 960 960 960 960 960 960 961 959 961
12 961 961 960 960 960 961 961 961 962 963 963 964 963 963 963 963 963 962 962 962 962 962 962 962 964
13 962 962 962 962 962 962 962 962 963 964 965 965 964 964 964 964 964 964 963 962 962 962 962 962 963 965
14 962 962 961 961 961 960 961 961 962 963 963 963 963 963 962 962 962 962 962 961 961 962 962 962 962 963
15 962 962 961 961 961 961 962 962 962 963 963 964 964 964 964 964 963 963 963 962 962 962 962 962 962 964
16 962 962 962 962 962 962 962 962 963 963 964 964 965 965 964 964 964 963 963 962 962 962 962 962 963 965
17 962 962 961 961 961 962 962 963 963 964 964 965 965 965 964 964 964 964 963 962 962 962 962 962 963 965
18 962 962 962 962 962 962 962 962 964 965 965 965 965 965 965 965 964 964 963 963 962 962 963 963 963 965
19 962 962 962 961 961 962 961 962 963 964 964 964 964 963 963 963 963 962 962 961 961 960 960 960 962 964
20 960 960 959 959 959 959 959 960 962 963 964 964 964 964 964 964 964 964 963 963 962 963 963 963 962 964
21 963 962 962 961 961 961 961 961 963 964 965 964 964 964 964 964 964 963 963 962 962 962 962 962 963 965
22 962 961 961 960 960 960 960 961 962 964 964 964 964 963 963 962 963 963 962 962 962 962 962 962 962 964
23 962 961 961 961 960 960 960 961 962 963 964 964 963 963 963 963 963 963 962 962 962 962 962 962 962 964
24 962 962 961 961 961 961 961 962 963 964 964 964 964 964 963 963 963 963 962 962 962 962 962 962 962 964
25 962 962 961 961 960 960 961 961 962 963 963 963 963 963 963 962 962 962 961 961 960 960 961 960 962 963
26 960 960 959 959 959 959 959 960 960 961 961 961 962 962 962 962 962 962 960 961 961 960 960 960 960 962
27 960 961 961 960 960 961 961 962 962 963 964 965 965 964 964 964 962 963 962 962 962 962 962 962 962 965
28 961 961 961 961 961 961 961 961 963 964 964 965 964 964 964 964 964 964 963 963 963 963 963 963 963 965
29 963 963 962 961 961 960 961 961 962 964 964 964 964 963 963 963 963 963 962 961 961 961 961 961 962 964
30 961 961 961 960 960 960 960 960 961 962 963 963 963 963 963 961 962 962 961 961 960 960 961 961 961 963
31 960 960 960 959 959 959 959 959 961 962 962 962 962 962 962 962 961 961 960 960 960 960 960 960 960 962
NOTAS1.-
: CASILLAS EN BLANCO = > DATO INVALIDADO POR CALIBRACION O FALLA PROMEDIO DEL MES 962
149