(Created page with " ==Introducción== En los países desarrollados del mundo se ha venido cambiando el método de generación de energía, dejando atrás la generación con energías provenien...") |
|||
Line 1: | Line 1: | ||
+ | <!-- metadata commented in wiki content | ||
+ | ==Variabilidad del costo de incertidumbre en la operación de generadores eólicos Variability in uncertainty cost for Wind Energy Generation (WEG)== | ||
+ | '''Luis Avellaneda<sup>a</sup>, Juan Monsalve<sup>b</sup>, Sergio Rivera<sup>c</sup>''' | ||
+ | --> | ||
− | ==Introducción== | + | ==Resumen== |
+ | Frente al problema de la variabilidad de la velocidad del viento en la generación eólica y la dificultad que esto conlleva en el despacho económico de energía eléctrica, en este se presentan desarrollos formulaciones matemáticas (funciones), denominadas costos de incertidumbre(UCF). El término se utiliza en un modelo compuesto de dos partes, uno implica costos de penalización por sobrestimación y el otro debido a la subestimación de una potencia programada. Mediante formulación matemática analítica y simulaciones de Montecarlo, se ha encontrado el valor esperado de los costos de incertidumbre, considerando la estocasticidad del recurso primario (viento) mediante distribuciones de probabilidad. Sin embargo, no se pueden evaluar con exactitud los resultados desde un punto de vista analítico, sin tener en cuenta la variabilidad alrededor del valor promedio. En este artículo se calcula la función analítica de la varianza, entendida como medida de dispersión de los valores esperados de los costos de incertidumbre para el caso concreto de los operadores eólicos; asimismo, se valida esta formulación mediante simulaciones de Montecarlo. La función encontrada es usada para observar la dispersión de los datos en torno a su valor medio esperado, mediante la aplicación de la desviación estándar que es sencillo determinar a partir de la varianza obtenida. | ||
+ | |||
+ | ''Palabras Clave: Energía Eólica Controlable, Costos de Incertidumbre, Estudios Probabilísticos, Simulaciones de Monte Carlo'' | ||
+ | |||
+ | ==Abstract== | ||
+ | |||
+ | Faced with the problem of wind speed in generation and the difficulty that this entails in the economic future of electric power, mathematical formulations (functions), called uncertainty costs, have been developed. The term is used in a model composed of two parts, an implication of penalization costs for overestimation and the other due to the underestimation of a programmed power. Through analytical mathematics and Monte Carlo simulations, we have found the expected value of the costs of uncertainty, the stability of primary resources through the distribution of probability. However, the results can not be evaluated accurately from an analytical point of view, without taking into account the variability around the average value. In this article we calculate the analytical function of the variance, we understand as a measure of dispersion of the expected values of uncertainty costs for the specific case of wind operators; Likewise, this formulation has been validated through Monte Carlo simulations. The function found is used to observe the dispersion of the data around its expected average value, by applying the standard deviation. | ||
+ | |||
+ | '''keywords''' Controllable wind energy, uncertainty costs, probabilistic studies, Monte Carlo simulation. | ||
+ | |||
+ | ==1 Introducción== | ||
En los países desarrollados del mundo se ha venido cambiando el método de generación de energía, dejando atrás la generación con energías provenientes de fuentes fósiles y nucleares, y dando paso a la implementación de la generación con energías renovables. Esta transición se ha notado con la construcción de grades parques de generación de tipo eólico. Sin embargo debido a la estocasticidad del recurso primario, este tipo de generación se considera como no programable en el despacho económico óptimo <span id='citeF-1'></span>[[#cite-1|[1]]]. El desarrollo de modelos de despacho económico teniendo en cuenta la estocasticidad de las fuentes primarias se empieza a estudiar en artículos como el desarrollado en el 2008 por John Hetzer, David C. Yu y Kalu Bhattarai <span id='citeF-2'></span>[[#cite-2|[2]]] , en el que plantean la necesidad de incorporar sistemas de conversión de energía eólica (WECS) en el problema de despacho económico. Este documento introduce los términos de costos de sobreestimación y subestimación de la potencia eólica disponible, añadiéndolos como factores al despacho económico clásico; para ello caracterizan la variabilidad del viento mediante la función de probabilidad de densidad de Weibull. Más adelante en 2012 JunHua Zhao, Fushuan Wen, Zhao Yang Dong, Yusheng Xue, and Kit Po Wong, en un artículo publicado por la IEEE <span id='citeF-3'></span>[[#cite-3|[3]]],proponen una función de distribución de la potencia eólica a partir del comportamiento de la velocidad del viento según la distribución de Rayleigh, desarrollan así un algoritmo de optimización basándose en un método de decodificación de Enjambre de Partículas Discreto (PSO) para solucionar el problema de despacho. Los artículos anteriores incorporan métodos heurísticos a la solución del despacho en generación eólica; en el mismo sentido el siguiente articulo <span id='citeF-4'></span>[[#cite-4|[4]]] expone el método ABC- Weighted-Sum, aplicándolo para obtener una mejor configuración que relacione los costos de generación con las emisiones en la generación térmica a gas, implementando la generación eólica para observar la disminución total en las emisiones. | En los países desarrollados del mundo se ha venido cambiando el método de generación de energía, dejando atrás la generación con energías provenientes de fuentes fósiles y nucleares, y dando paso a la implementación de la generación con energías renovables. Esta transición se ha notado con la construcción de grades parques de generación de tipo eólico. Sin embargo debido a la estocasticidad del recurso primario, este tipo de generación se considera como no programable en el despacho económico óptimo <span id='citeF-1'></span>[[#cite-1|[1]]]. El desarrollo de modelos de despacho económico teniendo en cuenta la estocasticidad de las fuentes primarias se empieza a estudiar en artículos como el desarrollado en el 2008 por John Hetzer, David C. Yu y Kalu Bhattarai <span id='citeF-2'></span>[[#cite-2|[2]]] , en el que plantean la necesidad de incorporar sistemas de conversión de energía eólica (WECS) en el problema de despacho económico. Este documento introduce los términos de costos de sobreestimación y subestimación de la potencia eólica disponible, añadiéndolos como factores al despacho económico clásico; para ello caracterizan la variabilidad del viento mediante la función de probabilidad de densidad de Weibull. Más adelante en 2012 JunHua Zhao, Fushuan Wen, Zhao Yang Dong, Yusheng Xue, and Kit Po Wong, en un artículo publicado por la IEEE <span id='citeF-3'></span>[[#cite-3|[3]]],proponen una función de distribución de la potencia eólica a partir del comportamiento de la velocidad del viento según la distribución de Rayleigh, desarrollan así un algoritmo de optimización basándose en un método de decodificación de Enjambre de Partículas Discreto (PSO) para solucionar el problema de despacho. Los artículos anteriores incorporan métodos heurísticos a la solución del despacho en generación eólica; en el mismo sentido el siguiente articulo <span id='citeF-4'></span>[[#cite-4|[4]]] expone el método ABC- Weighted-Sum, aplicándolo para obtener una mejor configuración que relacione los costos de generación con las emisiones en la generación térmica a gas, implementando la generación eólica para observar la disminución total en las emisiones. | ||
Line 9: | Line 24: | ||
El articulo mencionado anteriormente valida la función analítica del valor esperado con simulaciones de Montecarlo obteniendo resultados con bajo error. La velocidad del viento utilizada como variable aleatoria y la cantidad de datos utilizados en este estudio, hacen que esta medida de tendencia central sea insuficiente para describir el conjunto de datos, por ello se hace necesario elaborar un análisis estadístico de la dispersión de estos datos respecto al valor medio evaluado. Una medida de dispersión utilizada en este artículo es la varianza, que describe la esperanza del cuadrado de la desviación de cada variable respecto a la media de los datos. En los capítulos siguientes se hará el planteamiento de la función analítica de varianza a partir de su definición expresada en términos del valor esperado que se obtuvo en <span id='citeF-1'></span>[[#cite-1|[1]]]. Esta función analítica se validará con la simulación de Montecarlo de un generador eólico, teniendo en cuenta como variable aleatoria la velocidad del viento modelada como una función de probabilidad ya usada en WEG, la función de Weibull. Se presenta también un estudio de optimización teniendo en cuenta que potencia programada puede tener tanto un menor costo de incertidumbre así como una menor variabilidad de los datos, dicha optimización se propone de forma estadística utilizando la desviación estándar calculada a partir de la función de varianza encontrada para finalmente obtener conclusiones del trabajo realizado. | El articulo mencionado anteriormente valida la función analítica del valor esperado con simulaciones de Montecarlo obteniendo resultados con bajo error. La velocidad del viento utilizada como variable aleatoria y la cantidad de datos utilizados en este estudio, hacen que esta medida de tendencia central sea insuficiente para describir el conjunto de datos, por ello se hace necesario elaborar un análisis estadístico de la dispersión de estos datos respecto al valor medio evaluado. Una medida de dispersión utilizada en este artículo es la varianza, que describe la esperanza del cuadrado de la desviación de cada variable respecto a la media de los datos. En los capítulos siguientes se hará el planteamiento de la función analítica de varianza a partir de su definición expresada en términos del valor esperado que se obtuvo en <span id='citeF-1'></span>[[#cite-1|[1]]]. Esta función analítica se validará con la simulación de Montecarlo de un generador eólico, teniendo en cuenta como variable aleatoria la velocidad del viento modelada como una función de probabilidad ya usada en WEG, la función de Weibull. Se presenta también un estudio de optimización teniendo en cuenta que potencia programada puede tener tanto un menor costo de incertidumbre así como una menor variabilidad de los datos, dicha optimización se propone de forma estadística utilizando la desviación estándar calculada a partir de la función de varianza encontrada para finalmente obtener conclusiones del trabajo realizado. | ||
− | == | + | ==2 Descripción matemática de los costos de incertidumbre para generación eólica== |
En los sistemas eólicos la variabilidad del viento determina la potencia disponible de un generador, por tanto, dada una potencia programada existe la posibilidad de que la potencia disponible sea mayor o menor. Esta incertidumbre de potencia tiene asociado un costo: Si se tiene mayor potencia disponible de la que se programó, se desperdicia el recurso primario disponible, pues no se evaluó con exactitud el valor real de la velocidad del viento y por otro lado si por razones meteorológicas la velocidad de viento es baja junto con una potencia obtenida menor a la programada, significa que se evaluó por encima de la capacidad real que tiene el generador, entonces la potencia faltante para alcanzar la programada se ve traducida en un costo. .<span id='citeF-1'></span>[[#cite-1|[1]]]<span id='citeF-5'></span>[[#cite-5|[5]]] | En los sistemas eólicos la variabilidad del viento determina la potencia disponible de un generador, por tanto, dada una potencia programada existe la posibilidad de que la potencia disponible sea mayor o menor. Esta incertidumbre de potencia tiene asociado un costo: Si se tiene mayor potencia disponible de la que se programó, se desperdicia el recurso primario disponible, pues no se evaluó con exactitud el valor real de la velocidad del viento y por otro lado si por razones meteorológicas la velocidad de viento es baja junto con una potencia obtenida menor a la programada, significa que se evaluó por encima de la capacidad real que tiene el generador, entonces la potencia faltante para alcanzar la programada se ve traducida en un costo. .<span id='citeF-1'></span>[[#cite-1|[1]]]<span id='citeF-5'></span>[[#cite-5|[5]]] | ||
Line 169: | Line 184: | ||
Este resultado será fundamental para el planteamiento y desarrollo de la varianza. <span id='citeF-1'></span>[[#cite-1|[1]]]<span id='citeF-5'></span>[[#cite-5|[5]]] | Este resultado será fundamental para el planteamiento y desarrollo de la varianza. <span id='citeF-1'></span>[[#cite-1|[1]]]<span id='citeF-5'></span>[[#cite-5|[5]]] | ||
− | == | + | ==3 Formulación de la varianza en los costos de incertidumbre== |
La varianza es una medida de dispersión definida como la esperanza del cuadrado de la desviación respecto a la media y puede ser definida de la siguiente forma: | La varianza es una medida de dispersión definida como la esperanza del cuadrado de la desviación respecto a la media y puede ser definida de la siguiente forma: | ||
Line 257: | Line 272: | ||
|} | |} | ||
− | === | + | ===3.1 Varianza de costos de penalización por subestimar=== |
Como se puede ver en la ecuación (13) la varianza se puede dividir en dos partes, de la siguiente forma: | Como se puede ver en la ecuación (13) la varianza se puede dividir en dos partes, de la siguiente forma: | ||
Line 324: | Line 339: | ||
|} | |} | ||
− | === | + | ===3.2 Varianza de costos de penalización por sobrestimar=== |
De forma similar al planteamiento y solución de la varianza de costos de penalización por subestimar, se obtiene: | De forma similar al planteamiento y solución de la varianza de costos de penalización por subestimar, se obtiene: | ||
Line 403: | Line 418: | ||
Por definición <span id='citeF-2'></span>[[#cite-2|[2]]] los costos de penalización por subestimar ocurren cuando los de sobrestimar no están, solo existe una posibilidad de que se encuentren y es en un solo punto (<math display="inline">W_{w,s,i}</math>), al desarrollar matemáticamente este punto el resultado es cero; es decir que E[XY] se puede considerar como la suma de los valores esperados, con esto ya se tienen todos los valores analíticos. | Por definición <span id='citeF-2'></span>[[#cite-2|[2]]] los costos de penalización por subestimar ocurren cuando los de sobrestimar no están, solo existe una posibilidad de que se encuentren y es en un solo punto (<math display="inline">W_{w,s,i}</math>), al desarrollar matemáticamente este punto el resultado es cero; es decir que E[XY] se puede considerar como la suma de los valores esperados, con esto ya se tienen todos los valores analíticos. | ||
− | == | + | ==4 Validación de la varianza de los costos de incertidumbre mediante simulación == |
Con el objetivo de validar la fórmula encontrada analíticamente se plantea la simulación de Montecarlo para un generador eólico utilizando los parámetros de la tabla 1, la simulación consiste en 100000 datos aleatorios de la velocidad del viento teniendo en cuenta la distribución de Weibull adaptada específicamente a la Raylegh, estos datos se encuentran acumulados alrededor del valor del viento nominal dado en la tabla de parámetros, 15 m/s. | Con el objetivo de validar la fórmula encontrada analíticamente se plantea la simulación de Montecarlo para un generador eólico utilizando los parámetros de la tabla 1, la simulación consiste en 100000 datos aleatorios de la velocidad del viento teniendo en cuenta la distribución de Weibull adaptada específicamente a la Raylegh, estos datos se encuentran acumulados alrededor del valor del viento nominal dado en la tabla de parámetros, 15 m/s. | ||
Line 507: | Line 522: | ||
|} | |} | ||
− | == | + | ==5 Conclusiones== |
Este articulo propone la formulación matemática de la varianza de los costos de incertidumbre que representa la medida de dispersión de los valores esperados de los UCF para el caso específico de los generadores eólicos. Su importancia radica en que muestra la variabilidad de los costos de incertidumbre al analizar distintos escenarios de la potencia programada, específicamente en sistemas de generación eólicos donde la fuente de energía primaria (velocidad del viento) tiene asociada una gran incertidumbre. | Este articulo propone la formulación matemática de la varianza de los costos de incertidumbre que representa la medida de dispersión de los valores esperados de los UCF para el caso específico de los generadores eólicos. Su importancia radica en que muestra la variabilidad de los costos de incertidumbre al analizar distintos escenarios de la potencia programada, específicamente en sistemas de generación eólicos donde la fuente de energía primaria (velocidad del viento) tiene asociada una gran incertidumbre. |
Frente al problema de la variabilidad de la velocidad del viento en la generación eólica y la dificultad que esto conlleva en el despacho económico de energía eléctrica, en este se presentan desarrollos formulaciones matemáticas (funciones), denominadas costos de incertidumbre(UCF). El término se utiliza en un modelo compuesto de dos partes, uno implica costos de penalización por sobrestimación y el otro debido a la subestimación de una potencia programada. Mediante formulación matemática analítica y simulaciones de Montecarlo, se ha encontrado el valor esperado de los costos de incertidumbre, considerando la estocasticidad del recurso primario (viento) mediante distribuciones de probabilidad. Sin embargo, no se pueden evaluar con exactitud los resultados desde un punto de vista analítico, sin tener en cuenta la variabilidad alrededor del valor promedio. En este artículo se calcula la función analítica de la varianza, entendida como medida de dispersión de los valores esperados de los costos de incertidumbre para el caso concreto de los operadores eólicos; asimismo, se valida esta formulación mediante simulaciones de Montecarlo. La función encontrada es usada para observar la dispersión de los datos en torno a su valor medio esperado, mediante la aplicación de la desviación estándar que es sencillo determinar a partir de la varianza obtenida.
Palabras Clave: Energía Eólica Controlable, Costos de Incertidumbre, Estudios Probabilísticos, Simulaciones de Monte Carlo
Faced with the problem of wind speed in generation and the difficulty that this entails in the economic future of electric power, mathematical formulations (functions), called uncertainty costs, have been developed. The term is used in a model composed of two parts, an implication of penalization costs for overestimation and the other due to the underestimation of a programmed power. Through analytical mathematics and Monte Carlo simulations, we have found the expected value of the costs of uncertainty, the stability of primary resources through the distribution of probability. However, the results can not be evaluated accurately from an analytical point of view, without taking into account the variability around the average value. In this article we calculate the analytical function of the variance, we understand as a measure of dispersion of the expected values of uncertainty costs for the specific case of wind operators; Likewise, this formulation has been validated through Monte Carlo simulations. The function found is used to observe the dispersion of the data around its expected average value, by applying the standard deviation.
keywords Controllable wind energy, uncertainty costs, probabilistic studies, Monte Carlo simulation.
En los países desarrollados del mundo se ha venido cambiando el método de generación de energía, dejando atrás la generación con energías provenientes de fuentes fósiles y nucleares, y dando paso a la implementación de la generación con energías renovables. Esta transición se ha notado con la construcción de grades parques de generación de tipo eólico. Sin embargo debido a la estocasticidad del recurso primario, este tipo de generación se considera como no programable en el despacho económico óptimo [1]. El desarrollo de modelos de despacho económico teniendo en cuenta la estocasticidad de las fuentes primarias se empieza a estudiar en artículos como el desarrollado en el 2008 por John Hetzer, David C. Yu y Kalu Bhattarai [2] , en el que plantean la necesidad de incorporar sistemas de conversión de energía eólica (WECS) en el problema de despacho económico. Este documento introduce los términos de costos de sobreestimación y subestimación de la potencia eólica disponible, añadiéndolos como factores al despacho económico clásico; para ello caracterizan la variabilidad del viento mediante la función de probabilidad de densidad de Weibull. Más adelante en 2012 JunHua Zhao, Fushuan Wen, Zhao Yang Dong, Yusheng Xue, and Kit Po Wong, en un artículo publicado por la IEEE [3],proponen una función de distribución de la potencia eólica a partir del comportamiento de la velocidad del viento según la distribución de Rayleigh, desarrollan así un algoritmo de optimización basándose en un método de decodificación de Enjambre de Partículas Discreto (PSO) para solucionar el problema de despacho. Los artículos anteriores incorporan métodos heurísticos a la solución del despacho en generación eólica; en el mismo sentido el siguiente articulo [4] expone el método ABC- Weighted-Sum, aplicándolo para obtener una mejor configuración que relacione los costos de generación con las emisiones en la generación térmica a gas, implementando la generación eólica para observar la disminución total en las emisiones.
Para hacer posible la implementación de (WEG) Juan Camilo Arévalo, Fabian Santos y Sergio Rivera, en [1] se introduce la formulación analítica del valor esperado de los costos de incertidumbre, los cuales determinan costos de penalización por subestimar o sobrestimar la potencia programada, lo que permite al operador incluir estos costos en el valor de operación del generador eólico.
El articulo mencionado anteriormente valida la función analítica del valor esperado con simulaciones de Montecarlo obteniendo resultados con bajo error. La velocidad del viento utilizada como variable aleatoria y la cantidad de datos utilizados en este estudio, hacen que esta medida de tendencia central sea insuficiente para describir el conjunto de datos, por ello se hace necesario elaborar un análisis estadístico de la dispersión de estos datos respecto al valor medio evaluado. Una medida de dispersión utilizada en este artículo es la varianza, que describe la esperanza del cuadrado de la desviación de cada variable respecto a la media de los datos. En los capítulos siguientes se hará el planteamiento de la función analítica de varianza a partir de su definición expresada en términos del valor esperado que se obtuvo en [1]. Esta función analítica se validará con la simulación de Montecarlo de un generador eólico, teniendo en cuenta como variable aleatoria la velocidad del viento modelada como una función de probabilidad ya usada en WEG, la función de Weibull. Se presenta también un estudio de optimización teniendo en cuenta que potencia programada puede tener tanto un menor costo de incertidumbre así como una menor variabilidad de los datos, dicha optimización se propone de forma estadística utilizando la desviación estándar calculada a partir de la función de varianza encontrada para finalmente obtener conclusiones del trabajo realizado.
En los sistemas eólicos la variabilidad del viento determina la potencia disponible de un generador, por tanto, dada una potencia programada existe la posibilidad de que la potencia disponible sea mayor o menor. Esta incertidumbre de potencia tiene asociado un costo: Si se tiene mayor potencia disponible de la que se programó, se desperdicia el recurso primario disponible, pues no se evaluó con exactitud el valor real de la velocidad del viento y por otro lado si por razones meteorológicas la velocidad de viento es baja junto con una potencia obtenida menor a la programada, significa que se evaluó por encima de la capacidad real que tiene el generador, entonces la potencia faltante para alcanzar la programada se ve traducida en un costo. .[1][5]
Como se mencionó, el costo de incertidumbre(UCF) puede tener dos valores, uno donde hay un excedente de potencia que no se despacho y otro en el que hace falta para cumplir con lo programado, a cada uno de estos escenarios se los define como costos por subestimar y costos por sobrestimar, respectivamente.
Siendo la potencia disponible y la potencia programada se definen los costos por subestimar y sobrestimar de la siguiente manera:
|
(1) |
Siendo un coeficiente de penalización por subestimar
|
(2) |
Siendo un coeficiente de penalización por subestimar
Finalmente, los costos de incertidumbre corresponden al valor esperado de las funciones de costos de penalización por subestimar y sobrestimar E(X)[5], para este caso aplicados a la generación eólica se tiene:
|
(3) |
Donde UFC es la función de costos de incertidumbre. Para cada caso es desarrollada una integral:
|
(4) |
|
(5) |
donde:
es la PDF de la potencia WEG del generador i,
teniendo en cuenta que la función de probabilidad sigue la distribución de Weibull se puede plantear finalmente la función de costos de incertidumbre UCF. Para ello es necesario definir los límites o condiciones en las que se cumplen los costos de penalización por subestimar o sobrestimar, [1] y [5] definen estos limites de la siguiente manera:
|
(6) |
|
(7) |
|
(8) |
Aplicando estas condiciones se obtienen las siguientes funciones de costos de penalización por subestimar y sobrestimar [1]:
|
(9) |
|
(10) |
Este resultado será fundamental para el planteamiento y desarrollo de la varianza. [1][5]
La varianza es una medida de dispersión definida como la esperanza del cuadrado de la desviación respecto a la media y puede ser definida de la siguiente forma:
|
(11) |
Siendo E[X] el valor esperado de una función de distribución de probabilidad X.
Específicamente en este caso, E[X] es una función de costos de penalización por subestimar y sobrestimar; es decir, se tendrán dos resultados de varianza, uno debido a los costos de penalización por subestimar y otro por sobrestimar:
|
(12) |
|
(13) |
Para hallar la varianza total se debe tener en cuenta que los costos de penalización por subestimar y sobrestimar son eventos dependientes, como consecuencia se utilizó la siguiente formula:
|
(14) |
Es decir:
|
(15) |
Como se puede ver en la ecuación (13) la varianza se puede dividir en dos partes, de la siguiente forma:
|
(16) |
se soluciona a partir de metodos de sustitución simple y se puede deducir de (11). El resultado de la varianza por subestimar es el siguiente:
|
|
(17) |
|
(18) |
|
(19) |
De forma similar al planteamiento y solución de la varianza de costos de penalización por subestimar, se obtiene:
|
(20) |
Al igual que en el caso anterior, a partir de sustituciones simples y resultado de (8) es posible encontrar la solución, obteniendo el siguiente resultado:
|
|
(21) |
|
(22) |
|
(23) |
Para la varianza total de los costos de incertidumbre como ya se mencionó se utilizó la covarianza, la cual es obtenida con la siguiente formula:
|
(24) |
Por definición [2] los costos de penalización por subestimar ocurren cuando los de sobrestimar no están, solo existe una posibilidad de que se encuentren y es en un solo punto (), al desarrollar matemáticamente este punto el resultado es cero; es decir que E[XY] se puede considerar como la suma de los valores esperados, con esto ya se tienen todos los valores analíticos.
Con el objetivo de validar la fórmula encontrada analíticamente se plantea la simulación de Montecarlo para un generador eólico utilizando los parámetros de la tabla 1, la simulación consiste en 100000 datos aleatorios de la velocidad del viento teniendo en cuenta la distribución de Weibull adaptada específicamente a la Raylegh, estos datos se encuentran acumulados alrededor del valor del viento nominal dado en la tabla de parámetros, 15 m/s.
Entradas | ||
Variable | Descripción y Unidad | Valor |
Potencia programada al generador eólico [MW] | 100 | |
Velocidad de inicio de generación [m/s] | 5 | |
Velocidad nominal del viento [m/s] | 15 | |
Velocidad de parada de generador [m/s] | 45 | |
Potencia nominal de salida [MW] | 150 | |
Constante lineal [MW/m/s] | 15 | |
Constante independiente [MW] | -75 | |
Parámetro de escala Raylegh [m/s] | 15.9577 | |
Costo de penalización por sub estimar al generador [$ /MW] | 300 | |
Costo de penalización por sobre estimar al generador [$ /MW] | 700 | |
N | Número de iteraciones | 100000 |
![]() |
Figure 1: Velocidad del viento, Potencia disponible y Potencia programada |
Luego, en la figura 2 se presenta la cantidad de escenarios en los cuales se subestima y se sobrestima la potencia programada, como se definió en (1) y (2).
![]() |
Figure 2: Histograma costos por subestimar y sobrestimar |
![]() |
Figure 3: Histograma UCF |
![]() |
Figure 4: UCF al variar la potencia programada |
Finalmente se procede a calcular la varianza mediante simulación de Montecarlo al igual que se simulo los UCF, obteniendo (), al compararlo con el valor obtenido analíticamente se obtuvo un error de 0.0043.
Además, se hace el mismo análisis que con los costos de incertidumbre variando la potencia programada observando que cuando este tiene una menor varianza significa una menor dispersión para la cantidad de escenarios.
![]() |
Figure 5: Varianza al cambiar la Potencia programada |
Este articulo propone la formulación matemática de la varianza de los costos de incertidumbre que representa la medida de dispersión de los valores esperados de los UCF para el caso específico de los generadores eólicos. Su importancia radica en que muestra la variabilidad de los costos de incertidumbre al analizar distintos escenarios de la potencia programada, específicamente en sistemas de generación eólicos donde la fuente de energía primaria (velocidad del viento) tiene asociada una gran incertidumbre.
Pensando en el desarrollo y la alta penetración que comienzan a tener las energías renovables, el desarrollo de UFC presenta una gran aplicación en la optimización del despacho económico, pues minimiza riesgos tanto económicos como técnicos por sobrestimar o subestimar la potencia programada. Asimismo, el desarrollo matemático que se muestra de la varianza de UFC, añade un valor de validez y confianza a la optimización del despacho de generadores eólicos(WEG).
Finalmente el proceso matemático para obtener la varianza de UFC en WEG planteado, sirve como punto de partida para el futuro desarrollo de funciones matemáticas para varianza de UFC de otros sistemas de generación. Estos son sistemas de energías renovables o sistemas presentes en microgrids de generación distribuida, en los cuales ya se desarrollo un valor esperado de los costos de incertidumbre, como son generación fotovoltaica (PVG), vehículos eléctricos conectados a la red (PEV) y pequeñas centrales hidroeléctricas (PCHs)[6][7].
[1] Arevalo, J. C. Santos, F. and Rivera S. R. (2018) "Uncertainty Cost Functions for Solar Photovoltaic Generation, Wind Energy Generation, and Plug-In Electric Vehicles: Mathematical expected Value and Verification by Monte Carlo Simulation", Volume x. Power and Energy Conversion 1-37
[2] Hetzer, J. Yu, D. C. and Bhattarai K. (2008) "An Economic Dispatch Model Incorporating Wind Power", Volume 23. IEEE TRANSACTIONS ON ENERGY CONVERSION 1-9
[3] Zhao, J. Fushuan, W. Zhao, Y. Yusheng, X. and Kit P. W. (2012) "Optimal Dispatch of Electric Vehicles and Wind Power Using Enhanced Particle Swarm Optimization", Volume 8. IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS 1-11
[4] Haddi, S. and Tarek, B. (2015) "Economic/Emission Dispatch including wind power using ABC-Weighted-Sum", Volume 3. International Journal of Engineering and Technical Research (IJETR) 1-6
[5] Arevalo, J. C. Santos, F. and Rivera S. R. (2017) "Application of Analytical Uncertainty Costs of Solar, Wind and Electric Vehicles in Optimal Power Dispatch", Volume 22. Revista ingeniería 324-346
[6] Molina, F. S. Perez, S. J. and Rivera S. R. (2017) "Uncertainty Cost Function Formulation in Small Hydropower Plants Inside a Microgrid", Volume 8. Ingenierías USBMed 1-9
[7] Yao, Y. Wenzhong, G. and James Momoh. (2015) "Economic Dispatch for Microgrid Containing Electric Vehicles via Probabilistic Modeling", Volume 3. NATIONAL RENEWABLE ENERGY LABORATORY (NREL) 1-8
Published on 16/09/19
Accepted on 11/09/19
Submitted on 18/01/19
Volume 35, Issue 3, 2019
DOI: 10.23967/j.rimni.2019.09.003
Licence: CC BY-NC-SA license