Artículo Original

Evaluación del WRF como herramienta para determinar el rendimiento de un sistema fotovoltaico en condiciones de cielo despejado

Evaluation of the WRF as a tool to determine the performance of a photovoltaic system in clear sky conditions

Flavia Gutiérrez Muiña *
Centro de Gestión de la Información y Desarrollo de la Energía, CUBAENERGIA, Cuba
Alina Roig Rassi
Centro de Gestión de la Información y Desarrollo de la Energía, CUBAENERGIA, Cuba
Israel Borrajero Montejo
Instituto de Meteorología, Cuba
Bárbara Liz Miravet Sánchez
Centro de Gestión de la Información y Desarrollo de la Energía, CUBAENERGIA, Cuba

Evaluación del WRF como herramienta para determinar el rendimiento de un sistema fotovoltaico en condiciones de cielo despejado

Revista Cubana de Meteorología, vol. 28, núm. 4, e07, 2022

Instituto de Meteorología

Recepción: 01 Agosto 2022

Aprobación: 03 Agosto 2022

RESUMEN: En los sistemas fotovoltaicos se considera de vital importancia la evaluación de la eficiencia de la planta durante largos periodos con el objetivo de optimizar su rendimiento y maximizar su confiabilidad. La imposibilidad de adquirir en Cuba instrumentos medidores de la intensidad de la radiación solar para monitorear su desempeño, así como las limitantes para su mantenimiento y calibración, impone la necesidad de estimarla a través del empleo de modelos numéricos. En el presente trabajo se realiza la evaluación del modelo WRF como alternativa para determinar el rendimiento de una instalación fotovoltaica a partir de la estimación de la intensidad de la radiación solar en condiciones de cielo despejado. Se emplearon como valores de referencia datos de intensidad de la radiación solar medidos durante marzo de 2014 por una celda calibrada. Para realizar los ajustes entre la intensidad de la radiación solar en el plano horizontal (EICelda) y energía solar de onda corta acumulada en una hora en la superficie a cielo despejado (EIWRFCS), se emplearon 5 muestras con diferentes condiciones de cercanía entre las variables. Dentro de los resultados obtenidos vemos que, es posible emplear el modelo WRF como alternativa teniendo en cuenta que entre las horas de sol a analizar exista una diferencia ≤ 150 Wh/m2.

Palabras claves: Intensidad de la radiación solar a cielo despejado, rendimiento de sistemas fotovoltaicos.

ABSTRACT: In photovoltaic systems the evaluation of the efficiency of the plant during long periods is considered of vital importance in order to optimize their performance and maximize their reliability. The impossibility of acquiring in Cuba instruments to measure the intensity of solar radiation to monitor the performance photovoltaic systems, as well as the limitations for their maintenance and calibration imposes the need to estimate it through the use numerical models. In the present work, an assessment of the WRF model is carried out as an alternative to determine the performance of photovoltaic facilities, based on the estimation of the intensity of solar radiation under clear sky conditions. The solar radiation intensity data measured during the month of March 2014 by a calibrated cell is taken as a reference. 5 samples with different conditions of proximity between the variables were used to obtain fitted dependencies between the intensity of solar radiation in the horizontal plane (EICelda) and accumulated downwelling clear sky shortwave flux at bottom (EIWRFCS). It was obtained that it is possible to use the WRF model as an alternative, taking into account that there is a difference between the hours of sunshine to be analyzed difference ≤ 150 Wh/m2.

Keywords: Intensity of solar radiation in clear sky conditions, performance of photovoltaic systems.

INTRODUCCIÓN

El desarrollo industrial y tecnológico alcanzado a partir de la segunda mitad del siglo XX ha provocado un aumento en la explotación de los combustibles fósiles, situación que se refleja en dos aspectos fundamentales: el agotamiento, antes de lo previsto, de los yacimientos, y el incremento de la contaminación ambiental debido a la emisión de grandes cantidades de gases de efecto invernadero y de contaminantes a la atmósfera.

El concepto de fuente de energía alternativa surge hace más de 40 años, como aquella que puede sustituir a las energías convencionales, ya sea debido al menor efecto contaminante y/o por la posibilidad de renovarse (Stolik, 2019), constituyendo una parte importante de la solución hacia un desarrollo sostenible. El conjunto de las energías renovables comparte la misma fuente, el Sol (Morradán, 2012), del que se abastece directamente la energía fotovoltaica, la cual se presenta como alternativa ventajosa desde el punto de vista económico y ambiental.

La producción de energía eléctrica por medio de sistemas fotovoltaicos se caracteriza por ser fluctuante ya que depende directamente de los cambios en la intensidad de la radiación solar. En estos se considera de vital importancia la evaluación de su eficiencia durante largos períodos con el objetivo de optimizar su rendimiento y maximizar su confiabilidad. Los estudios que se llevan a cabo con este fin se deben realizar bajo cielo despejado; condición ideal para el funcionamiento del generador fotovoltaico y de los equipos de medición de intensidad de la radiación (como piranómetros y celdas fotovoltaicas), los cuales son calibrados bajo este requerimiento (Guinand & Quintero, 2020).

La problemática de la carencia de datos suficientes de radiación y el uso de equipamientos costosos y específicos ha provocado que se hayan realizado numerosos intentos para estimarla a partir de formulaciones sencillas con base empírica o teórica (Sayago, Ovando, & Willington, 2011). Entre los métodos de estimación se encuentran la utilización de imágenes de satélites y el empleo de modelos numéricos; la selección de uno u otro depende del horizonte de pronóstico y la resolución espacial y temporal necesarias.

En Cuba, el uso de modelos numéricos tiene particular importancia ya que las posibilidades de adquisición de instrumentos medidores de la intensidad de la radiación solar, su mantenimiento y calibración son limitadas. Al cierre de mayo del 2019, el país contaba con 65 parques fotovoltaicos instalados al Sistema Electroenergético Nacional (SEN) en funcionamiento y 15 en proceso de terminación. Muy pocos cuentan con instrumentos con la calidad adecuada para la medición de la radiación global que permitan monitorear el rendimiento del parque.

El presente trabajo tiene como objetivo principal evaluar el empleo del modelo WRF como alternativa para determinar el rendimiento de una instalación fotovoltaica a partir de la estimación de la intensidad de la radiación solar en condiciones de cielo despejado.

MATERIALES Y MÉTODOS

Área de estudio

Se seleccionó como área de estudio el polígono integral para el monitoreo y evaluación de paneles fotovoltaicos, ubicado en la azotea del Centro de Gestión de la Información y Desarrollo de la Energía (CUBAENERGÍA). A partir de contar con el sistema de monitoreo se pudo determinar el coeficiente característico del sistema fotovoltaico como:

C o e f c a r a c t = E n e r g í a e l é c t r i c a g e n e r a d a ( k W h ) E n e r g í a s o l a r i n c i d e n t e k W h m 2 A r e a d e s u p e r f i c i e c a p t a d o r a ( m 2 ) (1)

Modelo numérico WRF - ARW

Para la realización de este estudio se empleó también el modelo numérico WRF en su versión 3.9.1.1. con el núcleo dinámico ARW. Para las corridas se seleccionaron tres dominios anidados (Figura 1): un dominio externo con resolución de 9 km; un segundo dominio con resolución de 3 km y un tercer dominio de 1 km. El centro de todos los dominios coincide con el emplazamiento del polígono integral.

Dominios de modelación
Figura 1
Dominios de modelación

Los esquemas de parametrización utilizados se muestran en la Tabla 1.

Tabla 1
Esquemas de parametrización
Esquemas de parametrización

Para realizar las comparaciones entre las mediciones reales y estimadas por el modelo se utilizó la intensidad de la radiación solar incidente en el panel acumulada para cada hora, denominada EICelda; así como las variables ACSWDNB (Accumulated Downwelling Shortwave Flux At Bottom) y ACSWDNBC (Accumulated Downwelling Clear Sky Shortwave Flux At Bottom), obtenidas a través del modelo y nombradas EIWRF y EIWRFCS respectivamente.

Los datos recopilados pertenecen al mes de marzo de 2014 en un horario comprendido entre las 8:00 y las 19:00 horas.

Procesamiento de datos

Una vez que se tienen los datos de las tres variables para cada día se compararon los valores de intensidad de la radiación solar acumulados en una hora teniendo en cuenta las siguientes condiciones:

Para realizar los ajustes entre EICelda y EIWRFCS se emplearon datas seleccionadas a partir de diferencias existentes entre ambas variables (Tabla 2) y que cumplieran las condiciones previamente mencionadas.

Tabla 2
Datas seleccionadas para determinar ajustes entre EICelda y EIWRFCS
Datas seleccionadas para determinar ajustes entre EICelda y EIWRFCS

A partir de estos datos se determinó para qué diferencia se obtiene el mejor ajuste para estimar la intensidad de la radiación solar acumulada en una hora bajo condiciones de cielo despejado a partir del empleo del WRF-ARW.

Evaluación de los resultados

Para evaluar la calidad de los resultados obtenidos se realizó un análisis cuantitativo de los mismos a través de funciones estadísticas sencillas descritas por (Galán, Arce, Koch, & Lara, 2015). Los estadígrafos que se utilizaron fueron: el error medio (MBE) (Ecuación 2), el error porcentual medio (Ecuación 3), el error cuadrático medio (RMSE) (Ecuación 4) y el índice de correlación de Pearson (rP) (Ecuación 5).

M B E = 1 n ( E I W R F C S i - E I C e l d a i ) (2)

M P E % = 1 n ( E I W R F C S i - E I C e l d a i ) R C e l d a i 100 (3)

R M S E = 1 n ( E I W R F C S i - E I C e l d a i ) 2 1 2 (4)

r P = n E I W R F C S i E I C e l d a i - E I W R F C S i E I C e l d a i n E I W R F C S i 2 - E I W R F C S i 2 n E I C e l d a i 2 - E I C e l d a i 2 (5)

El resultado óptimo es aquel que consiga que los errores estadísticos tomen valor nulo y el coeficiente de correlación rP tome el valor de la unidad. Mientras más cerca estén los valores a los óptimos mayor será la exactitud del modelo empleado.

RESULTADOS Y DISCUSIÓN

Comparación de los valores de intensidad de la radiación solar

Los valores horarios acumulados de intensidad de la radiación solar estimados por el modelo (EIWRF y EIWRFCS) se comportan de manera similar a los medidos por la celda fotovoltaica los días 1ro (Figura 2 (a) ) y 8 de marzo, días típicos de cielo despejado; mientras que en el resto de los días la variable EIWRF es completamente diferente a las otras dos (Figura 2 (b) ), debido a la influencia de la nubosidad. Se puede apreciar también que el comportamiento de la variable EIWRFCS presenta la misma tendencia que el de la medida por la celda fotovoltaica, aunque el modelo presenta valores menores durante la mañana y mayores durante la tarde.

Comparación entre las variables EICelda, EIWRF y EIWRFCS para los días 1ro (a) y 28 de marzo (b)
Figura 2
Comparación entre las variables EICelda, EIWRF y EIWRFCS para los días 1ro (a) y 28 de marzo (b)

En cuanto a los valores acumulados en 24 horas, se encontró una variabilidad significativa en las estimaciones del modelo con respecto a la variable EIWRF, lo cual puede estar relacionado con las parametrizaciones escogidas que resultaron adecuadas para estimar los días de cielo despejado y presentaron grandes inexactitudes para los días con generación de nubes. La variable EIWRF se comporta de manera similar a EICelda solamente los días 8, 9 y 14 de marzo; mientras que EIWRFCS muestra similitud los días 1ro, 4, 5, 6, 7, 10, 11, 17, 24, 28, 29 y 30. En la Figura 3 se puede observar como los acumulados diarios de EIWRFCS sobrestiman a los medidos por la celda y presentan una distribución suave. Por su parte, los valores de EIWRF son subestimados y poseen una distribución más abrupta.

Marcha diaria de EICelda, EIWRF y EIWRFCS
Figura 3
Marcha diaria de EICelda, EIWRF y EIWRFCS

Por tanto, para determinar el desempeño de un sistema fotovoltaico a partir del modelo WRF- ARW se debe utilizar la variable EIWRFCS donde no hay influencia de la nubosidad.

Ajustes entre EICelda y EIWRFCS

Los resultados de los ajustes se muestran en la Figura 4, donde se aprecia como a medida que se disminuye el valor de la diferencia entre EICelda y EIWRFCS, ambas variables se acercan a la línea de tendencia y la pendiente de la ecuación se aproxima a 1, es decir, el valor estimado por el WRF (EIWRFCS) se aproxima al valor medido por la celda fotovoltaica (EICelda). Sin embargo, cuando disminuimos la diferencia, decrece la cantidad de datos que cumplen esa condición y la muestra no sería representativa en comparación con la población.

Ajustes entre EICelda y EIWRFCS para diferencias de 60 Wh/m2 (a), 75 Wh/m2 (b), 100 Wh/m2 (c), 150 Wh/m2 (d) y para todas las horas de sol (e)
Figura 4
Ajustes entre EICelda y EIWRFCS para diferencias de 60 Wh/m2 (a), 75 Wh/m2 (b), 100 Wh/m2 (c), 150 Wh/m2 (d) y para todas las horas de sol (e)

En la Tabla 3 se muestran los valores de los estadígrafos calculados. A medida que aumentan las diferencias seleccionadas, nótese que el error medio (MBE), el cuadrático medio (RMSE) y el porcentual medio (MPE%) aumentan, alejándose de cero (valor óptimo), lo que significa que la EIWRFCS se aleja de la EICelda. Por otra parte, a medida que aumenta la diferencia disminuye el índice de correlación de Pearson (rP), es decir, disminuye la correlación que existe entre EIWRFCS y EICelda; aunque de manera general los valores de rP son aceptables.

Tabla 3
Valores de los estadígrafos utilizados en el estudio.
Valores de los estadígrafos utilizados en el estudio.

Para determinar el ajuste más adecuado es necesario buscar un equilibrio entre la dispersión de los valores con respecto a la línea de tendencia, el porciento que la muestra representa del total y el valor de los estadígrafos aplicados.

A partir de los resultados mostrados en la Tabla 3 se observa que, aunque el error cuadrático medio más pequeño corresponde a la data 1, donde las diferencias entre las variables son menores o igual a 60, la muestra que cumple estas condiciones también es pequeña, pues solo contempla el 25.6 % de todas las horas empleadas en el estudio. Sin embargo, la data 4, donde las diferencias entre las variables son menores o igual a 150, a pesar de que su error cuadrático medio duplica a la data 1, contempla más del 50% de todas las horas empleadas en el estudio y tiene valores de estadígrafos y correlación de Pearson que se consideran aceptables para la evaluación realizada. Por esta razón se considera como ajuste más adecuado la data 4, que viene dada por la expresión:

E I W R F C S = 0.9125 E I C e l d a + 57.331

Determinación del rendimiento del sistema fotovoltaico

Para determinar el rendimiento de un sistema fotovoltaico se verificó cuan cercanos son los valores del coeficiente característico calculados según la ecuación 1 empleando como energía solar incidente los valores de EIWRFCS y EICelda, teniendo en cuenta las cinco datas definidas. Por último, se aplicaron los estadígrafos MBE y RMSE para verificar las diferencias entre ambas variables. Los resultados obtenidos se muestran en la Tabla 4.

Tabla 4
Valores de los estadígrafos para comparar los valores del coeficiente característico determinados a partir de EIWRFCS y EICelda.
Valores de los estadígrafos para comparar los valores del coeficiente característico determinados a partir de EIWRFCS y EICelda.

El menor valor de error medio (MBE) se obtiene para la data 4, que contempla más del 50% de todas las horas empleadas en el estudio; mientras que la data 1 muestra el menor valor del error cuadrático medio (RMSE), pero esta solo contempla el 25.6 % de todas las horas empleadas en el estudio.

CONCLUSIONES

La estimación de la intensidad de la radiación solar mediante el empleo del modelo WRF, con fines fotovoltaicos, tiene particular importancia en Cuba ya que las posibilidades de adquisición de instrumentos medidores, su mantenimiento y calibración es limitada.

Con este trabajo se ha demostrado la capacidad del modelo WRF como alternativa para evaluar el rendimiento de un sistema fotovoltaico a partir de la estimación de la intensidad de la radiación solar, para lo cual es necesario: la utilización de la variable Accumulated Downwelling Clear Sky Shortwave Flux At Bottom (ACSWDNBC, denominada en el trabajo EIWRFCS) para disminuir los efectos de la nubosidad, y que entre los valores medidos y estimados exista una diferencia menor o igual a 150 Wh/m2.

REFERENCIAS

Teixeira-Branco, V., Alonso-Suárez, R., de Almeida, E., Porrini, C., Gutiérrez, A., & Cazes, G. (2018). Evaluación del pronóstico de irradiación solar diaria en Uruguay utilizando el modelo WRF. In Anales del VII Congresso Brasileiro de Energia Solar, Gramado.

Cardona, M., & Piliougine, M. (2017). Medida de la irradiancia y temperatura. España: Máster Oficial en Tecnología de los sistemas de Energía Solar Fotovoltaica. Universidad Internacional de Andalucía.

Galán, R., Arce, A., Koch, C., & Lara, P. (2015). Modelo de cuantificación del potencial fotovoltaico de España. España: Colección Innovación y Conocimiento.

Quintero Santana, S., & Guinand Rodríguez, K. (2020). Protocolo y procedimiento aplicado a instrumentación piranómetro para medición de radiación solar.

Morradán, J. (2012). Piranómetro fotovoltaico. Pamplona, España: Tesis presentada para optar por el título académico de Ingeniero Técnico Industrial Eléctrico. Escuela Técnica Superior de Ingenieros Industriales y de Telecomunicación. Universidad Pública de Navarra.

Sayago, S., Bocco, M., Ovando, G., & Willington, E. A. (2011). Radiación solar horaria: modelos de estimación a partir de variables meteorológicas básicas. Avances en Energías Renovables y Medio Ambiente,15.

Stolik, D. (2019). Energía fotovoltaica para Cuba. La Habana, Cuba: Editorial CUBASOLAR.

Notas de autor

Contribución de los autores: Flavia Gutiérrez Muiña: Conceptualización, Conservación de datos, Análisis formal, Investigación, Metodología, Visualización, Redacción - borrador inicial. Alina Roig Rassi: Conceptualización, Investigación, Metodología, Administración de proyecto, Recursos, Redacción - revisión y edición. Israel Borrajero Montejo: Metodología, Recursos, Supervisión, Validación, Visualización. Bárbara Liz Miravet Sánchez: Redacción - revisión y edición.

*Autor para correspondencia: Flavia Gutiérrez Muiña. E-mail: flavia@cubaenergia.cu

Declaración de intereses

Conflicto de intereses: Los autores declaran que no existen conflicto de intereses
HTML generado a partir de XML-JATS4R por