Parameter Calculation of a Small Wind Turbine Model

Cálculo de parámetros del modelo de una turbina eólica

(160 w)

Carlos A. Ramírez Andrés Julián Saavedra-Montes Ph.D.

Departamento de Energía Eléctrica y Automática
Medellín, Colombia
caramireg@unal.edu.co, ajsaaved@unal.edu.co

Resumen— Un método para calcular los parámetros del coeficientes de una función prestablecida. En [5] y [6] se
modelo de una turbina eólica, de tres palas, eje horizontal y baja presenta una aproximación de la función del Cp a un
potencia es propuesto en este artículo. El método destaca un polinomio. Otra alternativa para calcular la curva son los
experimento en el cual se registran la velocidad del viento, la métodos estocásticos presentados en [6–8]. El más conocido de
velocidad angular del aerogenerador y la potencia eléctrica
entregada a la carga. Un polinomio de quinto orden que
estos es el método de bines, el cual consiste en un
representa la relación entre el coeficiente de potencia y la procedimiento de reducción de datos agrupándolos y
relación de velocidades es obtenido con los datos registrados. El promediándolos en intervalos denominados bines. Métodos de
método propuesto no requiere el uso de un túnel de viento y identificación y optimización para determinar los coeficientes
puede ser llevado a cabo en el lugar de instalación del de la función Cp(Ȝ) han sido presentados en [9–11].
aerogenerador. En este artículo se presenta un método experimental para
calcular los parámetros de una turbina eólica de baja potencia.
Abstract— A method to estimate the model parameters of a El artículo esta organizado de la siguiente forma: en la Sección
small wind turbine is proposed in this paper. The method
II se presenta el modelo estático de la turbina eólica de baja
highlights an experiment in which the wind speed, the generator
angular speed and the real power consuming by the load are potencia. El proceso para caracterizar el modelo de una turbina
registered. A fifth order polynomial to represent the ratio es presentado la Sección III y finalmente en la Sección IV se
between the power coefficient and the tip speed ratio is obtained presentan las conclusiones.
from the registered data. The proposed method does not require
a wind tunnel and can be carried out in the location of the wind II. MODELO ESTÁTICO DE UNA TURBINA EÓLICA DE BAJA
generator. POTENCIA
La potencia disponible en el viento Pv, calculada para un
Index Terms — Wind energy, power coefficient, bins, wind
área específica A, está dada por la Ec. 1 [12]:
turbine.
Palabras clave — Energía eólica, coeficiente de potencia, bines, ܲ௩ ൌ ͲǤͷǤ ߩǤ ‫ܣ‬Ǥ ‫ ݒ‬ଷ (1)
turbina de viento.

I. INTRODUCCIÓN Donde la potencia está dada en W, ȡ es la densidad del aire

en kg/m3, v es la velocidad del viento en m/s y A está dada en
Actualmente se quiere disminuir la dependencia de los m2. Las turbinas eólicas solo pueden extraer una porción de la
combustibles fósiles para proveer energía eléctrica, dado que potencia presente en el viento, esta porción está representada
son una de las causas del calentamiento global y sus por el coeficiente de potencia calculado con la Ec. 2:
consecuencias [1]. Una de las estrategias para acompañar el
uso de los combustibles fósiles son las fuentes alternativas de ‫ܥ‬௣ ൌ ܲ௧ Ȁܲ௩ (2)
energía eléctrica. Esta estrategia ha impulsado el desarrollo de
los sistemas eólicos en los últimos años [2].
Donde Pt es la potencia que es aprovechada por la turbina y
Para el análisis y control de los sistemas de generación
está en W. El coeficiente de potencia del rotor de una turbina
eólica es necesario contar con modelos que simulen su
eólica varía en función del tip speed ratio o Ȝ, el cual está
comportamiento. El modelo estático de una turbina eólica
definido como la relación entre la velocidad de barrido de la
posee tres parámetros: la densidad del aire, el área de barrido y
turbina y la velocidad del viento. Esta relación se muestra en la
la función del coeficiente de potencia Cp(Ȝ), la cual es conocida
Ec. 3:
como la curva de Cp. Algunos métodos para hallar la curva son
presentados en la literatura. En [3] y [4] se presenta un método
ߣ ൌ ߱Ǥ ܴȀ‫ݒ‬ (3)
basado en un filtro de Kalman, el cual encuentra el valor de los

978Ǧ1Ǧ4614Ǧ6/12/$31.00 ©2012 IEEE

 Donde Ȧ es la velocidad angular del eje de la turbina dada
en rad/s, y R es el radio de la turbina dado en m. La relación . TABLA I. DATOS DE PLACA DEL AEROGENERADOR.
Cp(ë) para turbinas eólicas de tres palas fijas, es decir donde el Parámetro Valor
Peso 6,00 Kg
ángulo de ataque del viento sobre las palas es igual a cero, se Velocidad de arranque 2,70 m/s
presenta en la Ec. 4 y la Ec. 5 [13]. Además la Fig. 1, la cual es Máxima velocidad del viento 49,17 m/s
presentada en la referencia [14], muestra la relación Cp(ë) para Potencia Nominal 160 W a 12,50 m/s
turbinas horizontales de 2, 3 y 20 palas. Radio 0,58 M

§c · 
C p (O ) c1 ¨ 2  c3 ¸ e Oi  c5 O
c4
A. Cálculo de la densidad del aire
© Oi ¹
(4)
El valor de la densidad del aire depende de la temperatura y

 0.035
la presión del lugar donde se registran los datos. Para calcular

Oi O
1 1 un valor promedio de la densidad del aire se instaló una
(5)
estación meteorológica que registra datos de temperatura,
humedad, presión atmosférica, velocidad del viento y dirección
Donde c1, c2, c3, c4 y c5 son coeficientes númericos con del viento cada 48 s. La estación almacena cada 5 min el
valores para cada rotor. promedio de los datos registrados cada 48 s.
El rango de medida para la velocidad del viento es de 0 -
44.4 m/s. El valor de la densidad del aire es calculado con la
Ec. 7, la cual es propuesta en la norma IEC 61400 [7], donde B
es la presión atmosférica dada en kg/m.s2, T es la temperatura
absoluta del aire en K y Ro es la constante del aire seco 287,05
J/kg/K. Los datos de presión y temperatura corresponden a la
información almacenada por la estación entre el 11 de abril del
2011 y el 13 de Julio del 2012. Estos datos se encuentran
disponibles en un sitio web [15]. El valor promedio de la
densidad es ȡ = 0,9838 Kg/m3:

ߩ ൌ ‫ܤ‬Ȁሺܴை ܶሻ (7)

B. Calculo del área de barrido de la turbina

El área de barrido es el área que está delimitada por las
Fig. 1. Comportamiento del coeficiente de potencia Cp para turbinas eólicas
palas de la turbina, en el caso de una turbina de eje horizontal y
con diferentes números de palas.
tres palas, el área de barrido corresponde al área de un círculo
de radio R y está dada por la Ec. 8, por lo tanto A = 1,075 m2.
El modelo estático de una turbina eólica de tres palas se
obtiene a partir de las Ec. 1 y 2, y se presenta en Ec. 6:
‫ ܣ‬ൌ ߨǤ ܴଶ (8)
ܲ௧ ൌ ͲǤͷǤ ߩǤ ‫ܣ‬Ǥ ‫ ݒ‬ଷ Ǥ ‫ܥ‬௣ ሺߣሻ (6)
C. Procedimiento para determinar la función Cp (Ȝ)
Los parámetros que se deben obtener para modelar una Para determinar la función Cp(Ȝ) es necesario registrar
turbina eólica son la densidad del aire ȡ, el area de barrido del datos durante un intervalo de tiempo y posteriormente calcular
rotor A y la función Cp (Ȝ) con sus correspondientes Cp y Ȝ para cada conjunto individual de datos hasta lograr
coeficientes. describir la función de Cp(Ȝ). Para calcular el Cp se utiliza la
Ec. 2, y es necesario medir o calcular la potencia de la turbina
III. CARACTERIZACIÓN DEL MODELO DE UNA TURBINA EÓLICA
y la potencia disponible en el viento. Para calcular Ȝ se utiliza
DE BAJA POTENCIA
la Ec. 3, y es necesario medir la velocidad del viento y la
La turbina eólica de baja potencia que se ha seleccionado velocidad angular del rotor.
para ser caracterizada está localizada en la sede del Instituto En la Fig. 2 se presenta el montaje propuesto para registrar
Tecnológico Metropolitano ITM de la Ciudad Medellín, latitud la señal de voltaje AC de salida del generador vAC, la señal de
6.245 N y longitud 75.550 W. La turbina eólica hace parte de voltaje de salida del anemómetro va y las señales de voltaje DC
un aerogenerador Air Brezze 40 cuyos datos de placa se vDC y corriente DC iDC en la carga electrónica. Las señales son
muestran en la tabla I. registradas con un osciloscopio digital. El osciloscopio está
A continuación se describe la estimación de los tres programado para enviar los valores de las siguientes señales: la
parámetros del modelo estático de la turbina eólica; la densidad frecuencia de la señal voltaje ac de salida del generador fvAC, la
del aire ȡ, el area de barrido del rotor A y la función Cp (Ȝ). cual se utiliza para calcular la velocidad angular del rotor Ȧ. El
osciloscopio también entrega la frecuencia de la señal de
voltaaje de salida del anemómeetro fva, la cuual se utiliza para consstante, corrien
nte constante, resistencia coonstante y pottencia
calcuular la velociddad del vientoo, y el voltaje promedio
p VDCC y la consstante. Los ranngos de operaación son: 0.1 – 120 V, 240 0 A, y
corriiente promediio IDC en la caarga, los cuales se utilizan para 12000 W. La carrga es prograamada en el modo resisttencia
calcuular la potenccia de salida ded la turbina. Cabe notar quue la consstante durantee la prueba. Ell valor de la reesistencia de carga
c
eficiiencia del geneerador sincrón nico de imanes permanentess y la RL se
s selecciona para extraer el e máximo dee la potencia de la
eficiiencia del recttificador se asumen iguales a uno, con essto la turbina frente a diferentes
d veloocidades de viiento. La seleccción
potencia eléctrica medida en la carga se aprooxima a la poteencia del valor
v de la ressistencia es preesentada en el ítem E.
de saalida de la turbbina y se calcuula como el prroducto del vooltaje L
Los datos en nviados al coomputador coorresponden a las
y la corriente prommedio de la caarga electrónicca. señaales fva, fvAC, IDC y VDC, con c esta infoormación se deben
d
calcular las variables v, Ȧ, Pv y Pt, paara esto se deben d
estabblecer las relaciones
r enntre señales medidas poor el
osciloscopio y lass variables requueridas para oobtener Cp (Ȝ)..
D. C
Caracterizacióón de las relacciones entre laas señales meddidas
y las variabless de la funciónn Cp (Ȝ).
Para medir la velocidad del vientto fue neceesario
P
caraacterizar el anemómetro quee viene integrrado en la estaación
meteeorológica. ElE anemómetroo funciona coon un encodeer, el
cuall genera una señal rectangular de voltajje cuya frecuuencia
esta dada por la velocidad
v de giro
g del anemóómetro, la cuaal está
dadaa por la veelocidad del viento. Parra caracterizaar el
anemmómetro se usó u un motor DC para emular e la rotaación
Fig. 2. Montaaje para registrar datos y determinaar la función Cp(ëë) provvocada por el viento en el sensor. Para medir la señ ñal de
Para
P la ubicaación de la estación metteorológica y del voltaaje de salida del
d anemómettro se usó el osciloscopio
o d
digital
aerogenerador se siguieron las recomendaciiones dadas por p la Tekttronix TDS 1000C-EDU,
1 simultáneameente se registtró la
norm
ma IEC 641000 [7], según laa cual la estacción meteoroló ógica veloocidad de viennto mostrada por la interffaz de la estaación
debee colocarse a una distanciaa del aerogeneerador entre 2 y 4 meteeorológica. La Fig. 4 pressenta el monttaje realizado para
vecees el diámetrro del rotor,, y entre esstas distanciaas se caraacterizar el aneemómetro.
recomienda una ddistancia de 2,52 veces el diámetro
d del rotor.
r L
La velocidad del motor DC D fue variadda a través de d su
Adem más, la estacción meteoroológica debe localizarse en e el voltaaje de alimenttación y conseecuentemente varió el registtro de
sentiido contrario de la direccción de viennto más espeerada veloocidad de viennto medido por la estación meteorológica
m a. Los
duraante la prueba.. En nuestro caso
c la distanccia tomada es de 3 datoos fueron regisstrados con un n intervalo de 2 min despuéés del
metrros en direccióón noroccidentte, ver Fig 3. cammbio del voltajje para que ell sistema lograara condiciones de
estaddo estacionariio. En la Fig. 5 se muestrann las interfacees del
osciloscopio, izqu uierda, y la estación meteoorológica, derrecha,
regisstradas en un n mismo instaante de tiemppo. La interfaz del
osciloscopio mueestra la señal de voltaje dell anemómetroo y la
frecuuencia de la misma
m y la innterfaz de la estación
e muesttra la
veloocidad del viennto equivalentte a la velociddad de rotacióón del
motoor. Para cadda velocidad del motor se registraro on la
frecuuencia de la señal
s de salidaa del anemómetro y la veloccidad
de viento
v equivalente. Los datoos registradoss se muestran en la
Fig. 6 y su relacióón es ajustada por la Ec. 9.

‫ ݒ‬ൌ ͲǤ͸͹ͻͺ݂‫ݒ‬௔ െ ͲǤͳ͵Ͳ͹ (9)

Fig. 3. Montajje del aerogeneraddor y la estación meteorológica.

Las
L señales sson registradaas con un ossciloscopio diigital
Tekttronix TDS30034C. El osciiloscopio está programado para
regisstrar las señalees durante un tiempo de meedida tm = 4s y para
enviar los datos ddurante un tiem mpo de envío te = 1s. Adem más la
frecuuencia de mueestreo es 2,5 kHz.
k Por lo tannto el oscilosccopio
envía al computaddor paquetes de d datos cada 5 s. Los paqu uetes
de datos contienenn a fvAC, fva, VDC
D e IDC.
L carga electtrónica DC programable BK
La B Precision 8514
puedde ser program mada para op perar en cuattro modos: vooltaje Fig. 4. Montaje de prrueba del anemóm
metro.
 de immanes perman nentes, con las cuales se cáálculo la resisttencia
a la cual se entregga la máxima potencia.
p
L
La prueba para construirr las curvas característicaas de
poteencia versus corriente y de d voltaje veersus corrientee del
geneerador consissten en increementar graduualmente la carga c
coneectada al geneerador para unau velocidad constante dee giro
del rotor
r del geneerador, mientrras se registraan la corrientee y el
voltaaje en la carrga, este prooceso se repite para diferrentes
veloocidades de girro del rotor deel generador.
P realizar la prueba se reetiró el conjunnto de las palas y el
Para
bujee del aerogenerador y se acoopló un motor DC para emuular el
Fig. 5. Interfaaces del osciloscoopio y la estación meteorológica
m giro producido por p el vientoo, para mantener la veloccidad
consstante se usó ele driver del motor
m DC, y para
p incremenntar la
cargga conectada al generadorr se usó una carga electró ónica
proggramable DC en modo corrriente constannte, con la cuual se
increementó la corrriente demanddada por la carrga en pasos de d 0,1
A. Con la info formación reccolectada de las pruebaas se
consstruyeron las curvas
c de las Fig.8
F y Fig. 9..

F
Fig. 6. Característiica del anemómettro

La
L Ec. 10 es uutilizada para calcular
c la veloocidad del rotor de Fig. 7. Señales de
d voltaje de salid
da del anemómettro y el generadorr
la tu
urbina, esta exppresión muesttra la relación que existe enttre la
velocidad angularr de un generaador y la frecuuencia de su señal E la Fig. 8 se
En s puede ver comoc la potenncia que entreega el
de voltaje de salidda. P es el núm
mero de polos del generadorr. geneerador para caada velocidad angular tienee un punto mááximo
de potencia,
p este punto
p se denoomina máximoo punto de pottencia
߱ ൌ ͶǤ ߨǤ ݂‫ݒ‬஺஼ Ȁܲ (10) mppp, en este puntto se puede caalcular la resisstencia equivaalente
Rmppp que se encuuentra con la corriente
c y el voltaje en el punto
p
E la Fig. 7 see presenta la seeñal va en el canal uno y la señal
En mppp. En la Fig 8 y en la Fig.9 F se trazaaron las líneaas de
vAC en el canal dos.
d Adicionaalmente la intterfaz muestraa las resisstencia constaante corresponndientes a los valores
v RL míínimo
frecuuencias de las señales, fva = 8,618 Hz y fvvAC = 52,45 Hzz. y RL máximo, entre e las que se encuentraan las resisteencias
Pv es calculadaa con la Ec.1,, el valor de velocidad de viento Rmppp para las diferentes
d velo
ocidades de girog del rotor.. Los
calcuulado con la Ec.
E 9 y el vallor estimado de d la densidadd del valoores de Rmpp se muestran en n la tabla II.
aire. Pt es aproxim
mada a la poteencia eléctricaa consumida porp la E valor de la resistencia dee carga RL, se calcula talquee con
El
ga, la cual es calculada con IDC, VDC y la Ec.
carg E 11: cuallquier valor ded velocidad de d giro de laa turbina se puedap
extraaer la máximaa potencia del generador, ccomo se pued de ver
ܲ௧ ؆ ܲ݁ ൌ ‫ܫ‬஽஼ Ǥ ܸ஽஼஼ (11) en la Fig.8 este valor
v de RL debe
d estar en el intervalo ded RL
míniimo y RL mááximo, por essto se tomó el e promedio de d los
L aproximacción de la Ec.. 11 es necesaria ya que no
La n se valoores de Rmpp p con el fin dee estar lo máás cerca posibble de
cuennta con un medidor
m de toorque para caalcular la poteencia todaas las resistenccias de Rmpp y con esto RL=1,269
= Ÿ
mecánica de laa turbina. Esta aproximaación incluye la
eficiiencia del generador y del rectificador enn el cálculo de Cp. F. E
El método de Bines.
B
E. Calculo
C de la rresistencia de carga RL P procesar la informació
Para ón obtenida dee las medidas de la
veloocidad del viiento, la velocidad angullar y la pottencia
Para
P programmar la carga electrónica DC D fue neceesario eléctrica con el fin
n de obtener la
l función del Cp (Ȝ) se utiliizó el
calcuular el valor de la resisten ncia de cargaa, este valor debe métoodo de bines que
q se presentta en el estándar IEC 61400--12-1
garaantizar la extraacción del máxximo de la pottencia de la turrbina [7], [8]. En este método
m los daatos se agrupaan en función de la
frentte a diferenntes velocidaddes de vientto, para estoo se veloocidad de vieento en interv valos de 0,5 m/s y para cada
consstruyeron las curvas caraacterísticas dee potencia veersus
corriiente y de volltaje versus coorriente de gennerador sincróónico
intervalo se encuentra el valor promedio de cada variable. Para este motivo se seleccionó el polinomio de quinto orden
calcular Cp el método de bines usa las Ec. 12, 13,14 y 15. presentado en la Ec. 16. La tabla III presenta los coeficientes
estimados para el polinomio de quinto orden y también incluye
el valor de la densidad del aire en el lugar donde se registraron
las medidas y el área de barrido de la turbina.

TABLA II. VALORES DE RMPP

Velocidad
IMpp A VMpp V RMpp Ÿ
angular rpm
579 6,997 11,590 1,656

528 6,498 9,930 1,528

486 5,997 9,060 1,511

453 4,997 8,600 1,721

350 4,997 6,200 1,241

Fig. 8. Curvas de potencia versus corriente 325 4,497 4,830 1,074

300 3,997 4,600 1,151

280 2,997 5,200 1,735

260 3,496 3,900 1,116

240 2,996 3,600 1,201

220 2,497 2,500 1,001

200 2,247 2,100 0,935

180 1,747 1,900 1,088

160 1,377 1,650 1,198

140 1,247 1,100 0,882

Fig. 9. Curvas de voltaje versus corriente

ே
‫ݒ‬௜ ൌ ሺͳȀܰ௜ ሻ σ௝ୀଵ
೔
‫ݒ‬௡ǡ௜ǡ௝ (12)
ே೔
߱௜ ൌ ሺͳȀܰ௜ ሻ σ௝ୀଵ ߱௡ǡ௜ǡ௝ (13)
ே೔
ܲ௜ ൌ ሺͳȀܰ௜ ሻ σ௝ୀଵ ܲ௡ǡ௜ǡ௝ (14)
ଷ
‫ܥ‬௣ǡ௜ ൌ ሺܲ௜ ȀͲǡͷǤ ߩ‫ݒܣ‬௜ ሻ (15)

Donde vi es la velocidad del viento normalizada y promedio

en el bin i, vn,i,j es la velocidad del viento normalizada del
conjunto de datos j en el bin i, Ȧi es la velocidad del viento Fig.10 Curvas ajustadas de Cp versus Ȝ
normalizada y promedio en el bin i, Ȧn,i,j es la velocidad del
viento normalizada del conjunto de datos j en el bin i Pi es la ‫ܥ‬௣ ൌ ܿଵ Ǥ ߣହ ൅ ܿଶ Ǥ ߣସ ൅ ܿଷ Ǥ ߣଷ ൅ ܿସ Ǥ ߣଶ ൅ ܿହ Ǥ ߣ ൅ ܿ଺ (16)
potencia eléctrica normalizada y promediada en el bin i, Pn,i,j
es la potencia eléctrica normalizada del conjunto de datos j en TABLA III. PARÁMETROS DEL MODELO DE UNA TURBINA EÓLICA DE BAJA POTENCIA
el bin j, Cp,i es el coeficiente de potencia en el bin i, Ni es el Parámetro Valor
número de datos en el bin i. c1 0,000104
Los datos tomados el 19 y el 28 de Julio del 2012 fueron c2 -0,002167
c3 0,014742
procesados con el método de bines y a partir de los datos c4 -0,033909
obtenidos de Cp y Ȝ se construyó la Fig.10, en está figura se c5 0,027660
ajustaron tres funciones con el fin de describir la función de Cp c6 -0,000512
y Ȝ. Las funciones utilizadas son polinomios de orden tres, ȡ 0,9838 kg/m3
cuatro y cinco. Los polinomios de orden tres y cuatro presentan A 1,075 m2
valores negativos haciendo que la función presente valores
negativos de Cp, lo cual contradice su definición, Ec. 2. Por
 IV. CONCLUSIONES [6] K. Mok, “Identification of the power coefficient of
Un método experimental para calcular los parámetros del wind turbines,” Power Engineering Society General
modelo de una turbina eólica, de tres palas, eje horizontal y Meeting, 2005. IEEE. pp. 2078–2082 Vol. 2, 2005.
baja potencia (160W) fue propuesto en este artículo. El método
es dependiente del viento disponible en el que se encuentra [7] I. E. C. Standard, “61400-12, International
instalado el aerogenerador. Un polinomio de quinto orden fue Electrotechnical Commission (IEC)(2002),” Wind
obtenido para representar la función Cp (Ȝ). La densidad del Turbine Generator Systems Part, vol. 12.
aire fue calculada con valores promedio de la temperatura y
presión atmosférica registradas en el lugar de ubicación del [8] A. Llombart, S. J. Watson, D. Llombart, and J. M.
aerogenerador. El área de barrido fue calculada con los datos Fandos, “Power Curve Characterization I: improving
del fabricante de la turbina. El método propuesto no requiere el the bin method,” in Proceedings of the International
uso de un túnel de viento y puede ser llevado a cabo en el lugar Conference of Renewable Energies and Power Quality
de instalación del aerogenerador. (ICREPQ), 2005.

V. AGRADECIMIENTOS [9] Z. Xu, Q. Hu, and M. Ehsani, “Estimation of Effective

Wind Speed for Fixed-Speed Wind Turbines Based on
Este trabajo fue soportado por el Departamento Frequency Domain Data Fusion,” IEEE Transactions
Administrativo de Ciencia, Tecnología e Innovación on Sustainable Energy, vol. 3, no. 1, pp. 57–64, Jan.
COLCIENCIAS, a través de la beca-pasantía Jóvenes 2012.
Investigadores de 2011 y de la beca doctoral 095-2005.
También fue soportado por la Universidad Nacional de
[10] A. Monroy, “Real-time identification of wind turbine
Colombia a través del grupo de investigación GAUNAL bajo
rotor power coefficient,” Power, vol. 3, no. 1, pp.
los proyectos IDENT_EXGEN y SMART_ALEN. Los autores
3690–3695, 2006.
agradecen el apoyo del Laboratorio de Electrónica Industrial y
de Potencia del Instituto Tecnológico Metropolitano ITM.
[11] A. Science, “Modeling and Optimization of Power
VI. REFERENCIAS Coefficient Using 2 K Factorial Methodology,”
Energy, pp. 2–7, 2010.
[1] L. A. Barroso, H. Rudrick, F. Sensfuss, and P.
Linares, “The Green Effect,” Power and Energy [12] O. Anaya-Lara, N. Jenkins, J. Ekanayake, P.
Magazine, IEEE, vol. 8, no. 5, pp. 22–35, 2010. Cartxright, and M. Hughes, Wind energy generatio:
Modelling and control, Fist. 2009, p. 289.
[2] R. Moreno, G. Strbac, F. Porrua, S. Mocarquer, and B.
Bezerra, “Making Room for the Boom,” Power and [13] R. I. Ovando, J. Aguayo, and M. Cotorogea,
Energy Magazine, IEEE, vol. 8, no. 5, pp. 36–46, “Emulation of a Low Power Wind Turbine with a DC
2010. motor in Matlab/Simulink,” Power Electronics
Specialists Conference, 2007. PESC 2007. IEEE. pp.
[3] H. Zhihong, Z. Yuan, and X. Chang, “State estimation 859–864, 2007.
for wind turbine system based on Kalman filter,”
Systems and Control in Aerospace and Astronautics, [14] M. E. Cruz, V. M. Guzmán Arguis, M. I. Giménez, H.
2008. ISSCAA 2008. 2nd International Symposium on. Fernández, and J. Restrepo, “Emulador eólico para
pp. 1–3, 2008. aerogeneradores de baja potencia,” Universidad,
Ciencia y Tecnología, vol. 13, pp. 329–338, 2009.
[4] D. Bourlis and J. a. M. Bleijs, “A wind speed
estimation method using adaptive Kalman filtering for [15] Instituto Tegnológico Metropolitano, “Datos estación
a variable speed stall regulated wind turbine,” 2010 meteorológica.” [Online]. Available:
IEEE 11th International Conference on Probabilistic http://galeon.com/metereologicomede/Medellin_folder
Methods Applied to Power Systems, pp. 89–94, Jun. /ITM.html. [Accessed: 14-Aug-2012].
2010.

[5] G. Son, S. Member, H. Lee, and J. Park, “Estimation

of Wind Turbine Rotor Power Coefficient Using RMP
Model,” Energy, 2009.