Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Cited by Google 
Similars in
SciELO 
Similars in Google 
Permalink
1. INTRODUCCIÓN
En Colombia, “las fuentes disponibles de información de recurso solar indican que el país cuenta con una irradiación promedio de 4,5 kWh/m2/d (UPME, IDEAM 2005), la cual supera el promedio mundial de 3,9 kWh/m2/d y está muy por encima del promedio recibido en Alemania (3,0 kWh/m2/d) (ARTICSUN, SF) país que hace mayor uso de la energía solar fotovoltaica”5. Así mismo, la región caribe no es la excepción, presenta amplios sectores de brillo solar medio diario en la Guajira y en el norte del Cesar, con promedios entre 7 y 9 horas de sol al día (hSd); en gran parte del centro y el resto del norte de la región, valores entre 6 y 7 (hSd); mientras que al sur de esta región está entre 4 y 6 (hSd)1), (11.
En virtud de lo anterior, la presente investigación estuvo enfocada en hacer un análisis completo de la zona que permitiera un aprovechamiento de la energía solar fotovoltaica, a través del uso de paneles solares para el suministro de electricidad en la operación de la planta de tratamiento de agua potable, las instalaciones administrativas y de laboratorio. Con la motivación de desarrollar un proyecto amigable con el medio ambiente, en cuanto a la disminución de emisiones asociadas al consumo de combustible; autosostenible en el tiempo dado que no requiere un pago de factura mensual por parte del acueducto para el suministro de energía, y que favorezca la continuidad del servicio de agua potable a los usuarios.
En cuanto al área de estudio, se seleccionó el municipio González que se encuentra ubicado al sur del departamento del Cesar, limita por el norte, oriente y occidente con el departamento de Norte de Santander y por el sur con el municipio de Río de Oro (Cesar). Sus coordenadas geográficas son 8°23’ latitud norte y 73°19’ longitud oeste de Greenwich. La cabecera municipal está a 1.240 metros sobre el nivel del mar, la temperatura varía entre los 19° y 22° grados centígrados, el clima es templado, ligeramente húmedo en la parte alta y semihúmedo en la parte baja6.
Con la finalidad de estimar el potencial de radiación en la zona de estudio, se obtuvo y analizó la información secundaria, procedente de la estación climática Aeropuerto Aguas Claras del municipio de Ocaña, tomando como referencia la distancia más próxima y que su latitud de 8,315°N es similar con respecto a la ubicación de la obra a 8,393°N, donde se tuvieron en cuenta datos de Tmáx: temperatura máxima, [°C] y Tmín: temperatura mínima [°C] de los años: 1985, 1996, 2001, 2011, y 2012. De igual manera se tomó información del aplicativo web disponible en la página POWER conjunto de datos solares de la NASA. Posteriormente, se estimó el potencial de radiación a través del modelo matemático propuesto por Bristow-Campbell, teniendo en cuenta la disponibilidad de información, y el buen ajuste presentado en los estudios publicados por el Dr. Becquer Frauberth Camayo, en Perú, tituladosDesarrollo del modelo Bristow Campbell para estimar la radiación solar global de la región de Junin”,2yValidación y aplicación del modelo Bristow Campbell para estimar la radiación solar global de la región de Junin3, en los que se discutieron las variables y condiciones de uso para el modelo trabajado. Así mismo, resalta su buen ajuste para el cálculo de la radiación; por lo que es posible afirmar que este tipo de modelos empíricos proveen una buena aproximación al dato real de la zona que se desee estimar. Cabe destacar que las condiciones de la zona de estudio del autor mencionado, comparte características similares a las de la presente investigación, en lo referente a la disponibilidad de la información.
Finalmente, partiendo de la ausencia de líneas de transmisión de energía eléctrica cercana a la ubicación del proyecto (zona no interconectada - ZNI), se planteó un dimensionamiento del sistema energético, de una instalación solar fotovoltaica autónoma, a partir del potencial de radiación de la zona, y los consumos estimados para suplir el requerimiento energético en la operación de la planta de tratamiento de agua potable, las instalaciones administrativas y de laboratorio del acueducto municipal. De acuerdo con la literatura se propone la siguiente ruta:6), (7), (10
Dimensionamiento del subsistema de generación: Se selecciona el tipo de panel solar y se determina el número de paneles necesarios.
Dimensionamiento subsistema de acumulación: Se realizará la selección del tipo de acumulador considerando la capacidad nominal diaria de la batería y la profundidad de descarga.
Dimensionamiento del regulador-inversor: Se determinó la potencia del regulador inversor en relación con la distribución de energía a través de los elementos que lo conforman.
Es de resaltar que cada vez son más los proyectos estructurados con criterio sostenible como el caso del estudio denominadoMarco de referencia para la planificación de generación distribuida en zonas no interconectadas, teniendo en cuenta los efectos ambientales, tales como el calentamiento global y las emisiones de gases8.
2. METODOLOGÍA
La Fig. 1 muestra el diseño metodológico llevado a cabo en la presente investigación, clasificada de tipo aplicada, descriptiva y experimental cualitativa. La cual se dividió principalmente en dos partes, la primera enfocada a determinar el potencial de radiación de la zona y la segunda proponer un dimensionamiento del sistema fotovoltaico.
PARTE I: POTENCIAL DE RADIACIÓN DE LA ZONA
La radiación solar es la energía procedente del sol en forma de ondas electromagnéticas. La irradiancia es la densidad de potencia incidente en una superficie9. La metodología que se muestra a continuación tomó como referencia el estudio desarrollado por el Dr. Becquer Frauberth Camayo, que es un modelo de tipo empírico propuesto por Bristow-Campbell (1984), que permite el cálculo de la radiación solar global, tal y como se describe en la Ecuación (1) y está en función de la temperatura y de la irradiación solar extraterrestre, definida como la fracción de la energía total incidente en la parte superior de la atmósfera.
(1)
Donde:
H g :Radiación solar global,[Wh/m 2 ] H e :Irradiación solar extraterrestre,[Wh/m 2 ] T max :temperatura máxima, [°C] T min :temperatura mínima, [°C] a: coeficientes b: coeficientes c: coeficientes
Para resolver el modelo planteado en la Ecuación (1) se deben determinar el valor de los coeficientes y el valor de la irradiación solar extraterrestre.
1.1 Determinación de los coeficientes para la zona de estudio: Inicialmente se ajusta el valor del coeficiente “a” que representa el máximo valor característico de cada área de estudio y depende de la latitud, y “b y c” que determinan el efecto del incremento del ΔT, sobre la temperatura de la atmósfera. Para estimar el coeficiente “a” se analizaron datos diarios de la Estación de Aguas Claras Aeropuerto Ocaña de Tmáx y Tmín; así como también los datos de Radiación solar Global [KWH/m2/día] de la página POWER NASA, para los años 1985, 1996, 2001, 2011, y 2012. El cual corresponde a un único valor a lo largo del periodo de estudio.
Los otros coeficientes “b” y “c” se calcularon diariamente con la Ecuación (2) y la Ecuación (3), respectivamente.
(2)
(3)
Donde:
Latitud (φ), [grados].
1.2 Estimación de la radiación solar extraterrestre (H eh ): Luego de tener el coeficiente “a” de la zona, se determina la radiación solar extraterrestre a partir del desarrollo de la Ecuación (4)
(4)
Donde:
I SC Constante solar0,082MJm -2 día -1 (E 0 ) = Factor de corrección de la excentricidad de la órbita terrestre (w sh ) = Ángulo horario para una superficie horizontal (δ) = Declinación solar en radianes
Como se puede evidenciar se desconocen los parámetros relacionados con el factor de corrección de la excentricidad de la órbita terrestre, el ángulo horario para una superficie horizontal y la declinación solar en radianes. Por tanto, se procede de la siguiente manera:
1.2.1 Determinación del factor de corrección de la excentricidad de la órbita terrestre: Se determinó a partir del desarrollo de la Ecuación (5)
(5)
Donde:
d n corresponde al número del día Juliano del año (1 ≤d n ≤ 365), Variando desde 1 para primero de enero, hasta 365 para 31 de diciembre.
1.2.2 Determinación del ángulo horario para una superficie horizontal: Se determinó a partir del desarrollo de la Ecuación (6)
(6)
Donde:
(w sh ) = Ángulo horario para una superficie horizontal
1.2.3 Determinación de la declinación solar en radianes: Se determinó a partir del desarrollo de la Ecuación (7)
(7)
Donde:
(δ) = Declinación solar en radianes
Una vez aplicado el modelo matemático y determinando el potencial de radiación en la zona de estudio, se realizó el dimensionamiento del sistema solar fotovoltaico4 7.
PARTE II: DIMENSIONAMIENTO DEL SISTEMA FOTOVOLTAICO
Para el dimensionamiento del sistema solar fotovoltaico se tuvo en cuenta el potencial de radiación disponible en la zona (enunciado en la PARTE I), el requerimiento energético en la operación de la planta de tratamiento de agua potable, las instalaciones administrativas y el laboratorio del acueducto municipal y el dimensionamiento del subsistema de generación, de acumulación y del regulador inversor12.
2.1 Estimación del consumo: Se determinaron las cargas y consumos de todos los elementos eléctricos de la planta de tratamiento de agua potable, las instalaciones administrativas y de laboratorio del acueducto municipal, teniendo en cuenta que este será operado por un sistema híbrido, hidráulico y solar. En todo caso cabe resaltar que el presente estudio se enfoca únicamente en el sistema solar fotovoltaico.
2.2 Dimensionamiento del sistema solar fotovoltaico: Se considera un montaje autónomo, provisto de un sistema de acumulación, que posibilita alimentar la carga en aquellos momentos que el sistema fotovoltaico no se encuentra operativo. A continuación, se muestran cada una de las partes que lo conforman:
2.2.1 Dimensionamiento del sistema de generación: Es formado por un grupo de paneles montados sobre una estructura que lo sujeta, buscando el máximo de radiación solar incidente. Para calcular el número de módulos necesarios se considera la Ecuación (8) y la Ecuación (9).
Escenario 1: El número total de módulos fotovoltaicos que se deben instalar se calcula a partir del consumo diario
(8)
Donde:
L md : Consumo medio diario mensual (wh/día) P MPP,STC : Potencia pico del módulo en condiciones STC (w p ) G d : Radiación solar global. También indica la HPS P G PR: Factor global de funcionamiento del generador
Escenario 2: Se estima la potencia mínima (P min) para cubrir las necesidades del consumo mencionado anteriormente y se calculan cuántos paneles cubren este requerimiento.
(9)
Donde
P min: Potencia mínima E D : Consumo diario G GEM : Irradiación estándar (1kw/m2) G dm (α,β): Irradiación sobre el generador P R : Rendimiento energético de la instalación
2.2.2 Dimensionamiento del sistema de acumulación: Sirve para almacenar energía y dar una tensión estable, aunque los paneles no estén recibiendo radiación. La carga eléctrica de una batería es el número de Ah que es capaz de almacenar9.
La (Energía diaria) se convierte en Ah que es la unidad de capacidad y se denomina consumo diario (Load daily), para lo anterior se divide la energía en el voltaje. De acuerdo con la Ecuación (10)
(10)
Donde:
E D : Energía diaria V: Voltaje
La capacidad de la batería se calcula según la Ecuación (11)
(11)
Donde:
C:capacidad del acumulador Ah A:autonomía del sistema en días. 6 horas, es decir 0,25 días en el que la batería dará energía y no recibirá. L D : Consumo diario medio de carga en Ah P Dmax : Profundidad de descarga máxima.PD= 70% η rd : Rendimiento energético del regulador + batería. Consumo regulador y batería un 20% de su potencia η rd 80% η inv : Rendimiento energético del inversor. Consumo inversor 15% de su potencia η inv 85%
2.2.3 Dimensionamiento del regulador: Controla las sobrecargas y las descargas en las baterías, lo que acorta su vida útil, evita la descarga de las baterías a través de los paneles durante la noche, cuando el voltaje de salida del panel FV es nulo.
La tensión nominal será la misma que la de la batería. Se eligió en función de la máxima corriente que puede controlar (ítem a y b).
De los paneles al regulador: La corriente es la de cortocircuito (ICC) del campo solar.
De la batería hacia las cargas. Se calcula la corriente de paso, a partir de la potencia de las cargas.
Adicionalmente, soporta al mismo tiempo una sobrecarga del 25% de las corrientes anteriores.
2.2.4 Dimensionamiento del inversor: Transforma la energía eléctrica continua en energía eléctrica alterna. Determinada a partir de la suma de potencia de todos los equipos que funcionan simultáneamente sobre el voltaje.
3. RESULTADOS
PARTE I: ESTIMACIÓN DEL POTENCIAL DE RADIACIÓN DE LA ZONA
El coeficiente “a” se determinó asociado a la pendiente de la línea recta, teniendo en cuenta el planteamiento que se evidencia a continuación, para el cual, se estableció un valor diario, para cada uno de los años seleccionados.
Ecuación de la línea recta:y = m*x
(12)
Donde:
: Corresponde al valor de “y” a: Valor de la pendiente
: Valor de “x”
A modo de ejemplo la Fig. 2 permite evidenciar la secuencia de cálculo para los primeros días del año 1985. Se procesaron un total de 1825 datos (columna de dn) correspondientes a los días analizados durante los 5 años del periodo de estudio. La columna de resultados se calcula una vez haya definido el coeficiente.
Con respecto al análisis global la Fig. 3 muestra un valor estimado para “a” de 0,7858.
Para ampliar el análisis y definir el valor del coeficiente se realizó un histograma de frecuencia tal y como se evidencia en la Fig. 4, donde el valor más próximo se situó en el rango comprendido entre [0,77 y 0,88]. Para instaurar el valor exacto se calculó el PERCENTIL.50 el cual da una aproximación más exacta al valor deseado.
Fuente: Los autores.
Fig. 4. Análisis de frecuencia de datos para la determinación de la constante “a”
Siendo como resultado más cercano, un valor para “a” de 0,79. Posteriormente se determinaron los coeficientes “b” y “c” los cuales estuvieron entre los rangos [1, 11-1, 62 y 0,02 -0,06], respectivamente.
Una vez estimado el valor de las contantes se procedió a determinar la radiación para los años 2016, 2017 y 2018, tal y como se evidencia en la Fig. 5, se realizó el análisis para cada día durante los tres años, estableciendo que los meses con menor radiación fueron enero y diciembre. Distribuidos en los primeros 31 días del año y los últimos 60 días. De igual manera, es posible apreciar que con respecto al año 2016 la radiación de 2018 presentó un aumento aproximado de 0,5 KWH/m2 día. A partir de los resultados obtenidos se establecen 5 horas de Sol Pico (HSP).
PARTE II: DIMENSIONAMIENTO DEL SISTEMA FOTOVOLTAICO
En relación con los consumos, cabe resaltar que el acueducto debe operar las 24 horas durante los 7 días de la semana de forma continua. La Tabla I muestra la caracterización de todos los consumos establecidos en la planta.
Se estimó un consumo diario de 17,31 KWh/día, para suplir el requerimiento energético de la planta de tratamiento de agua potable, las instalaciones administrativas y de laboratorio del acueducto municipal. En virtud de las diferentes facilidades eléctricas planeadas para instalar en la PTAP. La Fig. 6 muestra el perfil de carga para cada elemento.
De acuerdo con los voltajes, cargas y horarios de operación de cada elemento y en concordancia con las fuentes de energía seleccionada hidráulica y fotovoltaica se tiene lo siguiente:
Sistema fotovoltaico:
Turbina y/o fotovoltaico:
Cabe resaltar que para realizar el dimensionamiento del sistema solar fotovoltaico se definieron las especificaciones técnicas del sistema, teniendo en cuenta que la presente investigación constituye un soporte técnico teórico en relación con la selección del sistema que fue implementado en el acueducto municipal de González.
Dimensionamiento del sistema de generación: Se consideraron paneles solares de la marca TRINA 335 W, referencia TSM-PD14 constituido por 72 células módulo multicristalino. La Tabla II muestra las especificaciones dadas por el fabricante.
Para determinar el número de paneles solares se relacionan las ecuaciones (8), (9) y la Tabla II.
Escenario 1
De acuerdo con los cálculos realizados se requieren 12 módulos solares.
Escenario 2
Poniendo 12 paneles solares de 335W se obtiene una potencia de 4,02 kW, cubriendo las necesidades de la planta.
Dimensionamiento del sistema de acumulación: La batería seleccionada corresponde a la referencia LPG SERIES-GEL BATTERY, LPG12-200 [12V200AH].
La selección de la tensión se basa en la potencia requerida en los equipos, 2764,25 W (Tabla I) y en la Tabla III se presentan los voltajes recomendados para cada caso9.
Se considera una potencia de 48 V.
Se dimensiona la capacidad del banco de baterías, el cual debe poseer un Voltaje de 48V, con una autonomía de 6 horas diarias al 70% de descarga diaria, en ausencia de ambas fuentes de energía (hidroeléctrico y fotovoltaico).
La capacidad de la batería se calcula:
Por lo tanto, para suplir los 48 V y la capacidad de la batería, se seleccionaron 4 baterías de 12 V y 200Ah.
Dimensionamiento del regulador: El modelo seleccionado corresponde a un VS70, en relación con las características eléctricas, cuenta con una corriente máxima de salida de 70 A.9
De los paneles al regulador: De las especificaciones del panel solar seleccionado: ICC = 9,35 A. Y teniendo en cuenta que se instalaron 6 filas de 2 paneles, la ICCT = 9,35*6 = 56,1 A
De la batería hacia las cargas (Tabla I). Se calcula la corriente de paso, a partir de la potencia de las cargas.
Adicionalmente soporta al mismo tiempo una sobrecarga del 25% de las corrientes anteriores. La Tabla IV resume las condiciones de operación del regulador.
Se necesita un regulador de 48V, con una Ientra > 70A y una Isal > 58,57A
Dimensionamiento del inversor: Se seleccionó un inversor Studer Innotec referencia XTM 2600- 48 Xtender, el cual cumple con los requerimientos del sistema.
A La Fig. 7 muestra el diagrama preliminar de la configuración del sistema solar fotovoltaico
6. DISCUSIÓN Y CONCLUSIONES
El modelo matemático seleccionado corresponde al propuesto por Bristow-Campbell, dado que es el que presenta mejor ajuste con las variables disponibles en el área de estudio, corroborado con los resultados publicados por el Dr. Becquer Frauberth Camayo, en Perú. Donde las condiciones de la zona de estudio son similares a las presentadas en este trabajo.
Se estimó el potencial de radiación solar diaria de 5 KWh/m2/día utilizando el modelo propuesto por Bristow-Campbell, tomando como variables la temperatura máxima, la temperatura mínima, la radiación solar, latitud, día, meses, y año; proporcionado por la estación climática Aeropuerto Aguas Claras, del municipio de Ocaña, considerada la más cercana al área de estudio, así como la página del Proyecto POWER de la NASA.
Para el sistema solar fotovoltaico, el cual opera como una fuente de energía de respaldo, con una autonomía de operación de 6 horas al día, se realizó un dimensionamiento preliminar, el cual consta de:12 paneles solares de la marca TRINA 335 W, referencia TSM-PD14, constituido por 72 células módulo policristalino; 4 baterías de 12 V - 200Ah, cuya configuración permitió obtener 48V; un regulador VS70, con una corriente máxima de salida de 70 A; y un inversor Studer Innotec referencia XTM 2600-48 Xtender.
