MODELADO Y SIMULACIÓN DE SISTEMAS FOTOVOLTAICOS INTEGRADOS
EN REDES DE DISTRIBUCIÓN – CASO DE ESTUDIO SISTEMA
FOTOVOLTÁICO DE LA UNIVERSIDAD AUTÓNOMA DE OCCIDENTE
PEDRO ALEJANDRO MACANA TASCON
UNIVERSIDAD AUTÓNOMA DE OCCIDENTE
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE ENERGÉTICA Y MECÁNICA
PROGRAMA DE INGENIERÍA ELÉCTRICA
SANTIAGO DE CALI
2018
MODELADO Y SIMULACIÓN DE SISTEMAS FOTOVOLTAICOS INTEGRADOS
EN REDES DE DISTRIBUCIÓN – CASO DE ESTUDIO SISTEMA
FOTOVOLTÁICO DE LA UNIVERSIDAD AUTÓNOMA DE OCCIDENTE
PEDRO ALEJANDRO MACANA TASCON
2116235
Proyecto de Grado para optar al título de
Ingeniero Electricista
Director
PAUL ANDRÉS MANRIQUE CASTILLO
Ingeniero Electricista
UNIVERSIDAD AUTÓNOMA DE OCCIDENTE
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE ENERGÉTICA Y MECÁNICA
PROGRAMA DE INGENIERÍA ELÉCTRICA
SANTIAGO DE CALI
2018
Nota de aceptación:
Aprobado por el Comité de Grado en
cumplimiento de los requisitos exigidos
por la Universidad Autónoma de
Occidente para optar al título de
Ingeniero Electricista
RICARDO MORENO
Jurado
ROSAURADEL PILAR CASTRILLÓN
Jurado
Santiago de Cali, 05 febrero de 2018
3
CONTENIDO
pág.
RESUMEN 13
INTRODUCCIÓN 14
1. ASPECTOS GENERALES 17
1.1 SISTEMA ELÉCTRICO COLOMBIANO 17
1.2 ENERGÍA SOLAR FOTOVOLTAICA 21
1.2.1 Tipos de celdas fotovoltaicas 21
1.2.2 Principio de funcionamiento de las celdas fotovoltaicas 22
1.2.3 Parámetros de funcionamiento. 23
1.2.4 Efectos de la irradiancia y la temperatura 24
1.2.5 Modulo solar fotovoltaico. 25
1.2.6 Sistema solar fotovoltaico 25
1.2.7 Clasificación de los sistemas fotovoltaicos 26
[Link] Aplicaciones autónomas 26
[Link] Aplicaciones conectadas a la red 26
1.2.8 Energía solar fotovoltaica en el mundo. 27
1.2.9 Potencial de energía solar en Colombia 29
1.2.10 Autogeneración y cogeneración 30
[Link] Aplicaciones de autogeneración 32
1.2.11 Building Attached/Adapted/Applied Photovoltaics (BAPV) 32
1.2.12 Situación actual BAPV en el Valle del Cauca 33
2. CARACTERIZACIÓN 36
2.1 PROPUESTA METODOLOGICA PARA EL MODELADO Y
SIMULACION DE SISTEMAS FOTOVOLTAICOS 36
2.2 SISTEMA FOTOVOLTAICO 150 KWP (FASE 1) 37
2.1.1 Ubicación 38
2.1.2 Objeto 38
2.1.3 Características módulos fotovoltaicos marca Panasonic 38
[Link] Panasonic VBHN240 38
[Link] Panasonic VBHN230 39
2.1.4 Tipos de esquemas - Sistema fotovoltaico 150 kW 40
[Link] Parqueadero ala norte 40
[Link]. Parqueadero ala sur. 41
[Link]. Azotea edificio administrativo 42
2.1.5 Protecciones 43
2.1.6. Radiación 44
2.1.7. Rendimiento 45
2.3 SISTEMA FOTOVOLTAICO 250 KWP (FASE 2) 45
2.2.1. Ubicación 46
2.2.2 Objeto 46
2.3.3 Características módulo fotovoltaico marca Hareon 47
[Link] Hareon HR-265W 47
2.3.4 Tipos de esquemas - Sistema fotovoltaico 250 kW 48
[Link] Edificio - Aulas 48
[Link] Edificio bienestar universitario y cafetería 48
5
2.2.5 Protecciones 49
2.2.6 Rendimiento 49
2.4 CARACTERIZACIÓN DEMANDA GENERAL - RED ELÉCTRICA
UAO 51
2.3.1 Enero 52
2.3.2 Febrero 54
2.3.4 Marzo. 55
2.3.4 Abril 56
2.3.5. Mayo 57
2.3.6 Junio 59
2.3.7 Julio 60
3. DESCRIPCIÓN HERRAMIENTA ETAP PARA LA SIMULACIÓN DE
SISTEMAS FOTOVOLTAICOS 63
3.1 SISTEMA FOTOVOLTAICO EN ETAP 64
3.1.1 Módulo fotovoltaico (PV) 64
3.1.2 PV Panel page 66
3.1.3 PV Array page 69
3.1.4 Physical Page 74
3.1.5 Inverter page 74
3.2 ANÁLISIS DE FLUJO DE POTENCIA 77
4. VERIFICACIÓN DEL MODELO DE ETAP 78
4.1 INTERPRETACIÓN DE RESULTADOS UTILIZANDO TECNICAS
ESTADISTICAS 83
4.1.1 CORRELACIÓN 83
6
[Link] Potencia medida Vs irradiancia 83
[Link] Potencia simulada Vs irradiancia 84
[Link] Potencia medida Vs Potencia simulada 85
4.1.2 REGRESIÓN LINEAL SIMPLE 86
[Link] Ajuste lineal: Potencia medida Vs Potencia simulada 87
4.1.3 MAPE 88
[Link] Caso 1: Valores obtenidos mediante la simulación 89
[Link] Caso 2: Valores obtenidos mediante el ajuste de la simulación 90
5. INTEGRACIÓN DEL SISTEMA FOTOVOLTAICO CONECTADO A
RED 91
5.1 IDENTIFICACIÓN DE LA TOPOLOGÍA DE RED 91
5.2 MODELAMIENTO DE LA RED 92
5.3 RESULTADOS Y ANALISIS DE LA SIMULACIÓN EN ETAP 96
5.3.1 Periodo académico – ES 97
5.3.2 Periodo académico – FDS 99
5.3.3 Periodo de vacaciones – ES 101
5.3.4 Periodo de vacaciones – FDS 102
5.4 AULAS 1: GENERACION VS DEMANDA 104
6. CONCLUSIONES 108
BIBLIOGRAFIA 109
7
LISTA DE FIGURAS
pág.
Figura 1. Capacidad de generación SIN por tipo de combustible 2015 18
Figura 2. Generación térmica – Emisiones de CO2 20
Figura 3. Procesos ocurridos en una celda solar fotovoltaica 22
Figura 4. Efectos sobre la corriente y la tensión en un módulo 24
Figura 5. Capacidad y adiciones de energía solar FV 28
Figura 6. Capacidad y adiciones de energía solar FV, 10 países
líderes (2015) 28
Figura 7. Potencial de energía solar, algunas ciudades capitales 30
Figura 8. SFV Azotea edificio central – UAO 34
Figura 9. SFV Estacionamientos – Alkosto 34
Figura 10. SFV Azotea – Mercamio 35
Figura 11. Metodología para el modelado de SFV conectados en
redes de distribución 36
Figura 12. Vista superior del campus universitario 38
Figura 13. Diagrama unifilar SFV ala norte 41
Figura 14. Diagrama unifilar SFV ala sur 42
Figura 15. Diagrama unifilar SFV azotea 43
Figura 16. Disposición caja combinadora 43
Figura 17. Irradiación diaria promedio 44
Figura 18. Generación diaria promedio mensual 45
Figura 19. Emplazamiento sistema fotovoltaico fase II 46
Figura 20. Diagrama unifilar SFV Edificio - aulas 48
8
Figura 21. Diagrama unifilar SFV edificio bienestar universitario y
cafetería. 49
Figura 22. Generación diaria promedio SFV fase II – Segundo
semestre 2017 50
Figura 23. Diagrama unifilar UAO 52
Figura 24. Demanda diaria promedio – Enero 52
Figura 25. Demanda diaria promedio – Febrero 54
Figura 26. Demanda diaria promedio – Marzo 55
Figura 27. Demanda diaria promedio – Abril 57
Figura 28. Demanda diaria promedio – Mayo 58
Figura 29. Demanda diaria promedio – Junio 59
Figura 30. Demanda diaria promedio – Julio 60
Figura 31. Perfil de demanda mensual – Primer semestre 2017 62
Figura 32. Arreglo fotovoltaico 64
Figura 33. Curva I Vs V de un módulo fotovoltaico 65
Figura 34. Especificaciones de módulo fotovoltaico 66
Figura 35. Diferentes curvas de I-V 68
Figura 36. Librería de módulos fotovoltaicos 69
Figura 37. Especificaciones del módulo fotovoltaico 69
Figura 38. Conexión serie y paralelo 70
Figura 39. Calculadora de irradiación en ETAP 72
Figura 40. Physical page 74
Figura 41. Invertir page – PV Array editor 75
Figura 42. Editor de elemento inversor 76
Figura 43. Librería de conductores en ETAP 77
9
Figura 44. Modelo del sistema fotovoltaico en ETAP versión 12.6 78
Figura 45. Curva de generación diaria promedio anual 80
Figura 46. Generación de energía mensual (2016) 81
Figura 47. Flujograma para verificar el modelo de ETAP 82
Figura 48. Correlación Potencia medida Vs Irradiancia 84
Figura 49. Correlación Potencia simulada Vs Irradiancia 85
Figura 50. Correlación Potencia medida Vs Potencia simulada 86
Figura 51. Correlación Potencia medida Vs Potencia simulada
ajustada 88
Figura 52. Modelo red UAO – sistema fotovoltaico 96
Figura 53. Generación – Demanda Periodo Académico ES 98
Figura 54. Perfiles de tensión – Periodo Académico ES 98
Figura 55. Generación – Demanda Periodo Académico FDS 100
Figura 56. Perfiles de tensión – P. Académico FDS 100
Figura 57. Generación – Demanda Periodo Vacaciones ES 102
Figura 58. Perfiles de tensión – P. Vacaciones ES 102
Figura 59. Generación – Demanda Periodo Vacaciones FDS 103
Figura 60. Aulas 1: Generación Vs Demanda – Día hábil 105
Figura 61. Aulas 1: Generación Vs Demanda – Día sábado 106
Figura 62. Aulas 1: Generación Vs Demanda – Día domingo 107
10
LISTA DE TABLAS
pág.
Tabla 1. Capacidad efectiva neta del SIN a Julio de 2016 19
Tabla 2. Factores característicos del sistema eléctrico colombiano 19
Tabla 3. Registro de proyectos por tipo de tecnología 29
Tabla 4. Especificaciones eléctricas módulo VBHN240 39
Tabla 5. Especificaciones mecánicas módulo VBHN240 39
Tabla 6. Especificaciones eléctricas módulo VBHN230 40
Tabla 7. Especificaciones mecánicas módulo VBHN230 40
Tabla 8. Especificaciones eléctricas módulo HR-265W 47
Tabla 9. Especificaciones mecánicas módulo HR-265W 47
Tabla 10. Generación total de energía SFV fase II – Segundo semestre
2017 50
Tabla 11. Promedios de consumo por jornada - Enero 53
Tabla 12. Promedios de consumo por jornada - Febrero 54
Tabla 13. Promedios de consumo por jornada - Marzo 56
Tabla 14. Promedios de consumo por jornada - Abril 57
Tabla 15. Promedios de consumo por jornada - Mayo 58
Tabla 16. Promedios de consumo por jornada – Junio 59
Tabla 17. Promedios de consumo por jornada – Julio 61
Tabla 18. Promedios horarios de generación de energía por mes 79
Tabla 19. Total de energía generada mensualmente por hora 81
Tabla 20. Calculo de errores porcentual por mes (2016) 82
Tabla 21. Resultados caso 1 89
11
Tabla 22. Resultados caso 2 90
Tabla 23. Resultado filtrando registros según irradiancia 90
Tabla 24. Parámetros de transformadores de la red 91
Tabla 25. Temperatura horaria en la ciudad de Cali 92
Tabla 26. Irradiancia horaria promedio por mes 93
Tabla 27. Irradiancia promedio 93
Tabla 28. Demanda de potencia por hora Periodo Académico ES 94
Tabla 29. Demanda de potencia por hora Periodo Académico FDS 94
Tabla 30. Demanda de potencia por hora P. Vacaciones ES 95
Tabla 31. Demanda de potencia por hora P. Vacaciones FDS 95
Tabla 32. Resultados periodo académico ES 97
Tabla 33. Resultados periodo académico FDS 99
Tabla 34. Resultados periodo de vacaciones ES 101
Tabla 35. Resultados periodo de vacaciones FDS 103
Tabla 36. Producción de energía simulada para el SFV del edificio de
aulas 1 105
12
RESUMEN
Las energías renovables han alcanzado un elevado grado de desarrollo
tecnológico que ha permitido su rápida inserción en variados esquemas de
electrificación: desde los sistemas energéticos aislados, hasta los conectados a
redes eléctricas convencionales, pasando por esquemas de generación eléctrica
que complementan su operación con otras fuentes energéticas renovables y no
renovables. Esta variedad de aplicaciones pone de manifiesto la imperiosa
necesidad de evaluar la pertinencia de la integración de las energías renovables
dentro de esquemas tradicionales de generación convencional. En ese orden de
ideas, la finalidad principal de este trabajo, es la elaboración de un método que
permita analizar la integración de la energía solar fotovoltaica en redes de
distribución, a partir del uso de software de simulación de sistemas eléctricos de
potencia ETAP. Para ello, se escogió la red eléctrica de la Universidad Autónoma
de Occidente, que se constituye como la institución educativa con el sistema
fotovoltaico de mayor potencia instalada a nivel nacional, el cual es descrito en la
primera parte de este trabajo. Después de esto, se plantea un método, que
permite determinar la validez del modelo construido en ETAP, mediante la
herramienta estadística Minitab. Y por último, se procede a obtener los flujos de
potencia para diferentes condiciones operacionales del sistema.
Palabas claves: energía solar fotovoltaica, redes de distribución, software ETAP,
esquemas de electrificación.
13
INTRODUCCIÓN
Actualmente en Colombia los sistemas de distribución operan con el viejo modelo:
La generación sigue a la carga, que fueron implementados desde el siglo anterior,
donde la generación de energía eléctrica se realizaba en zonas apartadas de los
centros de consumo. En nuestro país, la mayor parte de la energía eléctrica
producida, se ha conseguido gracias a las grandes centrales hidroeléctricas
ubicadas en la zona centro del país, y plantas térmicas instaladas regularmente en
la costa; estas últimas, operan con combustibles fósiles, como el carbón, el diésel,
el gas, entre otros. Según el último boletín de la UPME, a mediados del año 2016,
la capacidad efectiva de generación para el recurso hídrico, era del 70% 1, siendo
la región andina y la pacífica, las zonas en donde se aprovecha más este recurso.
Alineados con los cambios tecnológicos que se han presentado alrededor del
mundo, Colombia adopto el 13 de Mayo de 2014 la ley 1715 que promueve e
incentiva la penetración de fuentes de energía renovable no convencionales
(FERNC) en el país y el uso eficiente de la energía eléctrica, y con ella hasta el
mes de Agosto del año 2016, se tiene el registro de una capacidad efectiva con
FERNC de 2.077,23 MW 2, la cual representa el 25,65% del total de proyectos de
generación registrados por la UPME a esa fecha. Estas condiciones de cambio
que se avecinan, sugieren de una alta preparación en los sistemas de distribución,
que si bien, seguirán realizando la misma labor, ahora tendrán que estar en
condiciones de permitir un flujo bidireccional, que provengan de distintas fuentes
de generación.
A raíz de los incentivos tributarios de la ley mencionada anteriormente, muchos
usuarios del sistema, en su mayoría comerciales e industriales, han decidido
apostar por sistemas de generación propia, que ayuden de alguna manera, a
solucionar o a respaldar, el suministro de energía eléctrica en sus instalaciones.
Esto marca una tendencia hacia la implementación de fuentes de energía propia
(autogeneración) en nuestro país, lo cual seguirá creciendo a medida que los
usuarios evidencien los ahorros económicos que se puedan obtener con la
integración de fuentes de energía renovable, que para el caso particular de este
trabajo, será la energía solar fotovoltaica.
1
UNIDAD DE PLANEACIÓN DE MINAS Y ENERGÍA (UPME). Ministerio de Minas y energías. [en
línea]. En: Boletín estadístico de minas y energía 2012 – 2016. Octubre de 2016. p. 18.
[consultado: 30 de marzo de 2017]. Disponible en Internet:
[Link]
2
UNIDAD DE PLANEACIÓN DE MINAS Y ENERGÍA (UPME). Ministerio de Minas y energías. [en
línea]. En: Registro de proyectos de generación Agosto 2016. 31 agosto 2017. p. 4. [consultado: 30
de marzo de 2017]. Disponible en Internet:
[Link]
pdf
14
Normalmente cuando se quiere implementar un sistema fotovoltaico en una
edificación, por ejemplo, se suele recurrir a un software como PVSyst o PVWatt
para estimar la generación de energía eléctrica que tiene una instalación durante
un año. Pero como su objetivo es el de estimar producción de energía, se deja
pasar por alto ciertos aspectos que son importantes para la operación del sistema,
como por ejemplo, el comportamiento de la red durante una jornada diurna con la
integración de un sistema fotovoltaico, conocer la proporción con la cual disminuye
la energía entregada por la red, gracias a la inyección de potencia de la instalación
fotovoltaica. Es precisamente esto, lo que se busca resolver en este trabajo
mediante la obtención de los flujos de potencia en la red eléctrica general de la
Universidad Autónoma de Occidente bajo distintas condiciones de operación.
En razón de los expuesto anteriormente, la estructura de este documento está
organizada en cinco capítulos que atienden al objetivo fundamental de analizar por
medio de la herramienta computacional de análisis de sistemas de potencia ETAP,
el comportamiento de los flujos de potencia derivados de la integración de un
sistema fotovoltaico en un red de distribución, que permitan determinar la
viabilidad de este tipo de instalaciones.
El Capítulo uno corresponde al marco introductorio de la investigación, en donde
se describe de manera general como está constituida actualmente la canasta
energética de Colombia, y los aspectos fundamentales de la energía solar
fotovoltaica. Además, se realiza una revisión acerca de algunos de los sistemas
fotovoltaicos de mayor capacidad instalados en la región del Valle del Cauca.
El Capítulo dos presenta la caracterización del sistema fotovoltaico de 150 kWp
que fue puesto en marcha a mediados del 2015, la segunda fase que entro en
operación en julio de 2017, y los perfiles de demanda que tuvo la institución
durante los primeros siete meses del 2017. Estos representan los requerimientos
necesarios para el modelado de la red eléctrica del caso bajo estudio.
El Capítulo tres corresponde a una descripción general de la herramienta
computacional ETAP, que permite conocer los requerimientos para el modelado
de sistemas de generación fotovoltaica.
En el Capítulo 4 se verifica el modelo construido en ETAP del sistema fotovoltaico
de 150 kWp, mediante un análisis de correlación, un ajuste lineal y el cálculo del
error porcentual absoluto medio (MAPE), que implican variables como: potencia
medida, potencia simulada e irradiancia.
15
Finalmente, en el Capítulo cinco se presentan los resultados obtenidos de las
simulaciones de flujo de potencia, integrando los sistemas fotovoltaicos de la fase I
(150 kWp) y II (250 kWp), con la red eléctrica de la Universidad Autónoma de
Occidente, bajo distintas condiciones operativas.
16
1. ASPECTOS GENERALES
1.1 SISTEMA ELÉCTRICO COLOMBIANO
Ana María Sandoval en su monografía del sector de electricidad y gas
colombiano: Condiciones actuales y retos futuros, señala que desde principios
del siglo XX, las redes eléctricas de Colombia, tenían como objetivo proporcionar
de alumbrado público a los municipios, que por lo regular eran sistemas
independientes, que eran alimentados mediante grandes centrales de generación
ubicadas en zonas rurales apartadas de los centros de consumo. A mediados del
siglo, bajo la intención de intercambiar energía entre regiones, ciudades y hasta
países, se crea la empresa ISA (Interconexión Eléctrica S.A.) que se encargaría
de construir y operar las líneas de interconexión. 3
El panorama nacional no ha sido muy distinto con respecto a los sistemas de
electrificación mundial, y es que la interconexión entre regiones del país, creó una
cadena productiva que va desde la generación, pasa por la transmisión y
distribución, para culminar con la comercialización. Estos cuatro dominios han
permitido en buena medida, mejorar la calidad de vida de los habitantes del
territorio colombiano. Aquí se evidencia, cómo después de pasado más de un
siglo, se sigue manejando un esquema convencional que maneja el modelo de la
generación sigue a la carga. A pesar de esto, el sistema eléctrico nacional ha
mantenido el objetivo de seguirse expandiendo, en busca de atender las
necesidades de nuevos usuarios, ofreciendo el servicio de la energía eléctrica, con
una buena calidad de potencia, de suministro y de servicio, e incluso, está
intentando cambiar poco a poco, el paradigma que se tiene, mediante los
proyectos de ley como la 1715 del 2014.
El país bajo sus condiciones climáticas, ha logrado establecer dos fuentes
principales de generación eléctrica. La que mayor predominación tiene, es la
hídrica, la cual representa el 70% de la capacidad efectiva neta, para mediados
del año 2016, y la segunda, es a base de combustibles fósiles (Plantas terminas),
como el gas y el carbón, entre otros, las cuales representan un 18% de la
capacidad instalada, sin considerar las térmicas menores. 4 Estos dos recursos se
han repartido en dos zonas del país, el primero, se consigue muy fácilmente en la
región centro del país, y debido a la ausencia del recurso hídrico en la región
caribe, se optó por la generación de energía eléctrica mediante las centrales
3
SANDOVAL, Ana María. Monografías del sector de electricidad y gas colombiano: Condiciones
actuales y retos futuros [en línea]. Archivos de economía. 2004. [consultado: 4 de abril de 2017].
Disponible en Internet: [Link]
4
UPME, Boletín estadístico de minas y energía 2012 – 2016, Op. cit., p.18
17
térmicas en esa zona del país. En la siguiente figura, se presenta la distribución de
las plantas de generación, en el SIN.
Figura 1. Capacidad de generación SIN por tipo de combustible 2015
Fuente: UPME. Capacidad efectiva de energía eléctrica SIN año 2015 [imagen].
Boletín estadístico de minas y energía 2012 – 2016. Octubre de 2016. P 14.
[Consultado: marzo 30 de 2017]. Disponible en Internet:
[Link]
18
Para el mes de julio de 2016, el sistema interconectado nacional tenía una
capacidad efectiva neta cerca de 16.540 MW 5. En la tabla 1, se muestra la
capacidad instalada del SIN, por tipo de tecnología.
Tabla 1. Capacidad efectiva neta del SIN a Julio de 2016
Centrales Capacidad Instalada (MW)
Hidráulica 11.512,95
Carbón 1.369,50
Gas 1.697,95
Eólica 18,42
Otros 1.940,80
Total 16.539,62
En 2015, Colombia generó 66,5 TWh de energía eléctrica 6, de los cuales el 33%
corresponde a generación mediante el recurso de los combustibles fósiles. En la
siguiente tabla se muestra desde el 2008 al 2016, el porcentaje de penetración de
generación térmica en el SIN, junto con el promedio de las emisiones de CO2
requeridas para entregar un kWh al sistema.
Tabla 2. Factores característicos del sistema eléctrico colombiano
Generación Promedio anual
Año
térmica (%) (gCO2/kWh)
2008 15% 95
2009 27% 164
2010 28% 185
2011 17% 102
2012 20% 125
2013 28% 192
2014 30% 200
2015 33% 221
2016 28% 190
La siguiente gráfica permite observar con mayor claridad, la evolución mensual de
la generación térmica, y las emisiones de CO2 durante los últimos nueve años.
5
Ibíd., p. 17.
6
Ibíd., p. 18.
19
Figura 2. Generación térmica – Emisiones de CO2
Fuente: Seguimiento a la huella de carbono de la generación de energía eléctrica
de las plantas despachadas centralmente en Colombia durante el año 2016. En:
XM – FILIAL DE ISA, p.2
De la anterior figura, se puede extraer la correlación directa que hay entre la
generación térmica del SIN, y las emisiones de CO2. El último y mayor incremento
que tuvo el promedio de emisiones de CO2, fue a finales del 2015 a raíz de los
cambios climáticos producidos por el fenómeno del Niño en nuestro país. Debido a
la relación que existe entre estas dos variables, surge la imperiosa necesidad de
tener una diversificación en la canasta energética del país. Para contribuir a ello,
se aprobó el plan de expansión de generación y transmisión para el periodo 2016
– 2030, con el que se busca adicionar durante los 15 años siguientes, 5.362 MW,
de los cuales el 37,7% se espera que sean provenientes de fuentes de energía
renovable no convencional7. Esto podría ser un buen motivante para los
departamentos de Antioquia, Meta, Caquetá, Valle del cauca y Santander, que
según el IDEAM en un informe presentado en noviembre de 2016, son las
regiones que más aportan en emisiones de efecto invernadero en el país 8, para
que empiecen a apostar por proyectos donde se involucren nuevas tecnologías de
generación, como lo son la solar, eólica o biomasa. Algunas empresas ya han
7
Colombia instalará 2.025 MW de energías renovables [en línea]. En: Revista Eólica y del Vehículo
Eléctrico. 24 de marzo de 2017. [Consultado: abril 1 de 2017]. Disponible en internet:
[Link]
8
IDEAM y PNUD presentan Inventario Nacional de Emisiones de Gases de Efecto Invernadero [en
línea]. Programa de las Naciones Unidas para el Desarrollo (UNPD). 3 de noviembre de 2016.
[Consultado: abril 1 de 2017]. Disponible en internet:
[Link]
[Link]
20
comenzado a implementar esquemas de autogeneración, que operan de manera
descentralizada. Estas soluciones, en su mayoría son sistemas de generación
fotovoltaica instalados normalmente en la terraza o en los techos de las
edificaciones. Por tanto, a continuación se abordara el concepto de energía solar
fotovoltaica.
1.2 ENERGÍA SOLAR FOTOVOLTAICA
La palabra fotovoltaica proviene de la palabra photo que representa luz, y voltaico
que hace referencia a electricidad, por tanto, la energía fotovoltaica aprovecha la
radiación solar para producir energía eléctrica. El principio se basa en la absorción
de la radiación solar por parte de un material semiconductor, que constituye las
denominadas células fotovoltaicas, provocando un desplazamiento de cargas en
su interior que originan la generación de una corriente DC. 9 Esto se produce a raíz
del efecto fotovoltaico, que fue descubierto por Becquerel en 1839, pero solo fue
hasta 1954 que pudo ser aprovechada por los laboratorios de Bell, quienes
desarrollaron la primer célula fotovoltaica funcional con un semiconductor de silicio
de unión pn, cuya eficiencia era del 6%. 10
1.2.1 Tipos de celdas fotovoltaicas. El material semiconductor más empleado
por las empresas fabricadoras de celdas fotovoltaicas es el silicio, un elemento
que puede encontrarse habitualmente en la arena de cuarzo (SIO2), cuya
disponibilidad es ilimitada pues es una materia prima inagotable. Dependiendo del
proceso de fabricación, las celdas fotovoltaicas se pueden clasificar en:
• Silicio Monocristalino (Si-m): El silicio se refina hasta obtener un bloque de
cristales muy puros, este se corta en láminas y es la base de las celdas. Si bien
estas celdas son de buena eficiencia, entre el 15 y 18%, son bastantes costosas.
En este caso el silicio está compuesto de cristales orientados de la misma forma.
La red cristalina es uniforme en todo el material y tiene muy pocas
imperfecciones. 11
9
FUNDACIÓN DE LA ENERGÍA DE LA COMUNIDAD DE Madrid (FENERCOM). Guía básica de la
generación distribuida. [en línea]. fenercom p. 22. [Consultado: 3 de abril de 2017]. Disponible en
Internet: [Link]
[Link]
10
GRUPO NAP. Energía Solar Fotovoltaica [En linea]. Colegio Oficial de Ingenieros de
Telecomunicación. Madrid: Ibergraphi, 2002. P 20. [consultado: 2 de junio de 2017]. Disponible en
internet: [Link]
11
UPME. Costos indicativos de generación eléctrica en Colombia. [En línea]. Integral: Ingeniería de
consulta. Abril de 2005. 250 p. [Consultado: 5 de junio de 2017]. Disponible en internet:
[Link]
[Link]
21
• Silicio Policristalino (Si-P): Este tipo es menos puro, con menos costo y menor
eficiencia, entre el 10 y el 16%, no está formado por cristales orientados de la
misma forma. El proceso de cristalización no es tan cuidadoso y la red cristalina
no es uniforme en todo el material, pudiéndose apreciar las zonas donde los
cristales tienen la misma orientación. 12
• Silicio Amorfo (Si-A): Las celdas fotovoltaicas más recientes emplean silicio
amorfo y se reconocen como celdas de película delgada. Su menor costo a pesar
de su baja eficiencia que esta alrededor de 7 a 10%, lo hace atractivo. En el silicio
amorfo no hay red cristalina y se obtiene un rendimiento inferior a los
anteriormente nombrados, sin embargo, posee la ventada, además de su bajo
costo, de ser un material muy absorbente por lo que basta una fina capa para
captar la luz solar 13.
1.2.2 Principio de funcionamiento de las celdas fotovoltaicas. La UPME en el
documento Costos indicativos de generación eléctrica en Colombia,
mencionan que: cuando la luz del sol incide sobre el material, los fotones generan
pares electrón-hueco tanto en el semiconductor tipo p como en el n. Los
electrones entonces librados de sus enlaces, se comportan ahora como electrones
libres, difundiéndose por el material. Igualmente ocurre con los huecos14. En la
siguiente figura se presenta un esquema simplificado del funcionamiento de una
celda solar fotovoltaica.
Figura 3. Procesos ocurridos en una celda solar fotovoltaica.
Fuente: UPME. Procesos ocurridos en una celda solar fotovoltaica [imagen]. Costos
indicativos de generación eléctrica en Colombia. P 251. [Consultado: marzo 5 de junio de
2017]. Disponible en Internet:
tp://[Link]/siel/documentos/documentacion/generacion/costos_indicativos_generacion_e
[Link]
12
Ibíd., p. 250.
13
Ibíd., p. 250.
14
Ibíd., p. 251.
22
Los electrones que se encuentran en exceso en la capa n del semiconductor
migran hacia la p y los huecos de la capa p del semiconductor, hacia la n. De esta
manera, se establece un campo eléctrico fijo en la unión pn. Esta región es el
responsable de establecer la corriente eléctrica a partir de los electrones. El
potencial asociado a la unión pn depende del material que conforma el arreglo de
la celda, la concentración de impurezas y de los portadores intrínsecos. Los
portadores de carga minoritarios en cada alcanzan la unión, en donde el campo
eléctrico existente se encarga de desplazarlos a la región donde estos electrones
y huecos son portadores de carga mayoritaria. La región p se carga positivamente
mientras la capa n negativamente, dando lugar a una fuerza electromotriz, que
permite la circulación de la corriente por el circuito exterior. De esta manera se
consigue el efecto fotovoltaico. Sin embargo, no toda la energía del sol puede
convertirse en energía eléctrica, pues parte de ella, se refleja en la celda o pasa a
través de ella y en algunos casos los electrones se recombinan antes de llegar a la
región p o n. Adicionalmente, cuando la radiación solar es muy débil, no es
suficiente para liberar los electrones, pero cuando es muy intensa, solo parte de
esa energía se usa para liberar los electrones, el resto se convierte en calor 15.
1.2.3 Parámetros de funcionamiento. Las celdas fotovoltaicas presentan unos
parámetros de funcionamiento que dependen de la configuración, y las
condiciones de operación de la misma. A continuación, se presentan estas
características 16:
• Intensidad de cortocircuito (Isc): Corresponde a la intensidad de corriente
eléctrica que se obtiene de la célula cuando, tras haber cortocircuitado sus
terminales, la tensión en bornes es nula. Es la máxima corriente obtenible. Valores
típicos de 10 a 40 mA por cm^2 de celda.
• Tensión en circuito abierto (Voc): Es la tensión en bornes máxima que puede
obtenerse de la célula, cuando no hay conectado ningún consumo y la intensidad
es nula. En dispositivos de Si de tipo medio se sitúa alrededor de los 0,6V,
mientras que los de GaAs es alrededor de 1V.
• Punto de máxima potencia (Imppt y Vmppt): La potencia eléctrica extraída de
la celda corresponde al producto del voltaje y la corriente. Esta potencia es cero
tanto en circuito abierto como en cortocircuito de la celda. Existe un punto de
máximo valor de la potencia denominado potencia pico, el cual corresponde al
valor que entrega la potencia bajo unas condiciones estándar de la radiación
incidente, entendida como el producto de la corriente y la tensión pico.
15
Ibíd., p. 252.
16
MANRIQUE, Paul. Celdas y módulos [Diapositivas]. Curso de sistemas fotovoltaicos conectados
a red. Universidad Autónoma de Occidente. 2016. p. 16, 17, 18.
23
• Factor de forma (FF): Es el cociente entre la máxima potencia obtenible sobre el
producto de los valores extremos de la celda (Isc y Voc). Básicamente el factor de
forma es una medida de la desviación de la rectangularidad, siendo un indicador
de la calidad de la celda.
• Eficiencia: La eficiencia de la conversión energética de una celda se define
como el cociente entre la potencia máxima que puede ser suministrada a una
carga y la potencia de la radiación solar incidente.
1.2.4 Efectos de la irradiancia y la temperatura. La mayoría de los módulos
fotovoltaicos son fabricados bajo unas condiciones de prueba estándar (STC por
sus siglas en ingles), que son: irradiancia de 1.000W/m^2, la incidencia de la luz
es perpendicular al módulo, la distribución espectral es de 1,5 A.M., y la
temperatura nominal de operación de celda (NOCT por sus siglas en inglés) es de
25°C. Estas condiciones no se van a tener en la práctica y es por esta razón, que
tanto la corriente y la tensión del módulo, se verán afectadas, por lo cual el
rendimiento se verá disminuido. Para la corriente generada por la celda, es
prácticamente proporcional a la radiación incidente sobre ella, siendo despreciable
el efecto de la radiación sobre el incremento o decremento de la tensión de la
celda. En la siguiente figura, se aprecia dicho comportamiento. Mientras que el
incremento o decremento de la temperatura afecta sensiblemente la cantidad de
tensión generada en la célula, siendo menor a incrementos de la temperatura y
viceversa 17. En la siguiente figura, se aprecian los comportamientos de la corriente
(Izquierda) y la tensión (Derecha).
Figura 4. Efectos sobre la corriente y la tensión en un módulo
Fuente: Manrique, Paul. Celdas y módulos [Diapositivas]. Curso de sistemas
fotovoltaicos conectados a red. Universidad Autónoma de Occidente. 2016. p. 20
17
MANRIQUE, Paul. Celdas y módulos [Diapositivas]. Curso de sistemas fotovoltaicos conectados
a red. Universidad Autónoma de Occidente. 2016. p. 20
24
1.2.5 Modulo solar fotovoltaico. Un módulo fotovoltaico es un conjunto de
celdas interconectadas eléctricamente y protegidas contra la intemperie,
generalmente, tienen una cubierta frontal de vidrio templado y un marco de
aluminio que facilita su transporte e instalación. Debido a que no tienen partes
móviles, los módulos son muy confiables y duraderos. Algunos fabricantes los
garantizan hasta por 25 años contra defectos de fábrica y reducción de
rendimiento, además requieren de mínimo mantenimiento.
1.2.6 Sistema solar fotovoltaico. Para alcanzar los niveles de voltaje y potencia
requeridos para suplir alguna necesidad de consumo, se construyen los sistemas
fotovoltaicos, combinando módulos en serie y/o paralelo para conformar un
arreglo. Una de las principales ventajas de los sistemas fotovoltaicos, es la
flexibilidad con la que se pueden concebir los diseños y posterior instalación de los
mismos, debido a su modularidad. Es decir, se pueden establecer configuraciones
de modulo, dependiendo de los parámetros de potencia que tenga el inversor, al
cual se va a conectar el arreglo de módulos.
La curva corriente-voltaje (IV) del arreglo se asemeja a la curva de los módulos
individualmente, con la única salvedad, que los parámetros estarán dados al
número de módulos que se tengan conectados tanto en serie como en paralelo.
Para el caso de la tensión, solo se verá modificada, cuando los módulos se
conectan en serie. Por tanto, el voltaje del sistema fotovoltaico estará dado por la
siguiente expresión
Vmppt SFV = �Vmppt del modulo� ∗ (numero de modulos en serie) [1]
Para el caso de la intensidad de corriente, su valor estará establecido por el
número de módulos que se conectes en paralelo. Por tanto, la corriente nominal
del sistema fotovoltaico se establece de la siguiente manera:
Imppt SFV = �𝐼𝐼mppt del modulo� ∗ (numero de modulos en paralelo) [2]
Por último, la potencia nominal del arreglo, dependerá de la suma de la potencia
nominal de cada módulo. Entonces:
Pmppt SFV = �𝑃𝑃mppt del modulo� ∗ (numero de modulos del arreglo) [3]
El diseño de un sistema fotovoltaico depende esencialmente de los requerimientos
del consumidor. Por lo regular, los sistemas fotovoltaicos se constituyen por el
25
modulo, controladores de carga, batería, inversor y las cargas, que pueden ser AC
o DC, dependiendo del tipo de aplicación que se le vaya a dar. Por tal razón, a
continuación se presentan algunas aplicaciones en las cuales se emplea un
sistema fotovoltaico.
1.2.7 Clasificación de los sistemas fotovoltaicos. Los sistemas fotovoltaicos
son un conjunto de componentes mecánicos, eléctricos y electrónicos que
concurren para captar y transformar la energía solar disponible, en energía
eléctrica. La clasificación de estos sistemas, se puede realizar en función de la
aplicación. Entre ellas se distingues dos, que son las aplicaciones autónomas y la
otra conectado a red. Siendo esta ultima la que se tratará en el trabajo.
[Link] Aplicaciones autónomas. Como su nombre lo indica, son instalaciones
que producen electricidad sin tener ningún tipo de conexión con la red de uso.
Algunas de las aplicaciones pueden ser: 18
• Telecomunicaciones: Telefonía rural, vía radio.
• Electrificación de zonas rurales y aisladas: Estas instalaciones, que se pueden
realizar en cualquier lugar, están pensadas para países y regiones en desarrollo y
todas aquellas zonas en que no existe acceso a la red eléctrica.
• Señalización: Se aplica a señales de tráfico luminosas, formadas por diodos
LED, alimentados por un módulo fotovoltaico y una batería.
• Bombeo de agua: Estas instalaciones están pensadas para lugares tales como
granjas, ranchos, etc. Se pueden realizar en cualquier lugar. Su uso puede ser
tanto para agua potable como para riego.
[Link] Aplicaciones conectadas a la red. En este caso, los sistemas
fotovoltaicos son conectados a las redes de distribución, y normalmente en
nuestro país, este tipo de soluciones están dirigidas a cubrir cierta parte de la
demanda en horas del día, y en la jornada de la noche, la red es quien cubre toda
la demanda. Las dos aplicaciones más distinguidas son 19:
18
DIAZ, Tomas. Instalaciones solares fotovoltaicas. (técnico en instalaciones eléctricas y
automáticas). [En línea]. España: McGraw-Hill / interamericana. 2010. p. 10. [Consultado 20 de
mayo de 2017]. Disponible en Internet:
[Link]
19
Ibíd., p. 11.
26
• Granjas solares: recintos en los que se concentra un número determinado de
instalaciones fotovoltaicas de diferentes propietarios con el fin de vender la
energía eléctrica producida a la compañía eléctrica con la cual se haya establecido
el contrato. En Colombia, la primera central fotovoltaica empezara operaciones a
final de este año, la cual tiene una potencia instalada cercana a los 10 MW, y
estará instalada en el municipio de Yumbo en el departamento del Valle del
Cauca.
• Edificios fotovoltaicos: Esta es una de las aplicaciones que en los últimos años
se ha venido desarrollando. Actualmente se ha empleado el termino BIPV
(fotovoltaica integrada a edificaciones) que consiste en combinar un elemento
constructivo del edificio, con células fotovoltaicas para generar energía, es decir,
los productos tienen doble funcionalidad. Hay otro concepto que se tratara más
adelante y corresponde al BAPV (fotovoltaica adaptada a edificaciones), el cual
consiste en instalar módulos fotovoltaicos sobre la terraza de los edificios para el
aprovechamiento de la energía solar. Este tipo de generación de energía eléctrica,
es el más adecuado para las zonas urbanas, pues no provoca efectos ambientales
adversos.
Los sistemas fotovoltaicos es una de las tecnologías que se utilizan como
generación distribuida, que para el Unión de América de distribución de potencia
(DPCA) se define como: “cualquier tecnología de generación a pequeña escala
que proporciona electricidad en puntos más cercanos al consumidor que la
generación centralizada y que se puede conectar directamente al consumidor o a
la red de transporte o distribución”20. Por tanto, a continuación, se realiza una
breve descripción de las posibles topologías con las cuales se conectan a la red,
estos sistemas de generación.
1.2.8 Energía solar fotovoltaica en el mundo. Esta fuente de energía renovable
no convencional ha sido la tecnología con la mayor pendiente de crecimiento con
respecto las otras, y es que tan solo en el año 2015, según REN21 (Red de
políticas en energía renovable para el siglo 21) se adicionaron 50 GW a la
capacidad de generación a nivel mundial. 21 En la siguiente figura, se puede
apreciar la cantidad de adiciones que se han tenido año tras año, desde el 2005.
20
FUNDACIÓN DE LA ENERGÍA DE LA COMUNIDAD DE Madrid (FENERCOM). Guía básica de
la generación distribuida. [En línea]. Fenercom p. 9. [consultado: 3 de abril de 2017]. Disponible en
Internet: [Link]
[Link]
21
REN21. Energía renovable 2016: Reporte de la situación mundial. [en línea]. ren21. p. 22.
[Consultado: 10 de abril de 2017]. Disponible en Internet: [Link]
renewables/global-status-report/
27
Figura 5. Capacidad y adiciones de energía solar FV
Fuente: REN21. Energía renovable 2016: Reporte de la situación mundial.
[Imagen]. REN21. p. 22. [Consultado: 10 de abril de 2017]. Disponible en Internet:
[Link]
Los cuatro países que mayor participación tienen a nivel mundial en energía solar
fotovoltaica son China, Alemania, Japón y Estados Unidos, siendo el primero, que
mayor aumento tuvo, en su matriz energética con la tecnología solar fotovoltaica
en el año 2015. Para ese año, la capacidad efectiva de energía solar fotovoltaica
era de 277 GW, y China era participe, con cerca del 16%. Japón fue el segundo
país con mayor aumento en esta tecnología, adicionando para ese año 11 GW.
Figura 6. Capacidad y adiciones de energía solar FV, 10 países líderes (2015)
Fuente: REN21. Energía renovable 2016: Reporte de la situación mundial.
[Imagen]. ren21. p. 22. [Consultado: 10 de abril de 2017]. Disponible en Internet:
[Link]
28
La energía solar fotovoltaica en Colombia, se podría decir, que es una tecnología
emergente, y que apenas, en los últimos años ha empezado a implementarse.
Muestra de ello, son los proyectos que se registraron en la UPME en el año de
2016. Resaltado se puede observar, que energía solar es la segunda tecnología
con mayor número de proyectos registrados, pero está en la posición número 4, si
de capacidad se hablara, siendo superada por la hídrica, térmica, y la eólica. Esto
podría deberse, a que los proyectos de generación fotovoltaica, están destinados
para autoconsumo, que bien podrían estar conectado en sistemas de distribución.
Tabla 3. Registro de proyectos por tipo de tecnología
Fuente: UNIDAD DE PLANEACIÓN DE MINAS Y ENERGÍA (UPME). Ministerio
de Minas y energías. Registro de proyectos de generación Agosto 2016. [Imagen].
31 agosto 2017. p. 4. [Consultado: 30 de marzo de 2017].
1.2.9 Potencial de energía solar en Colombia. En Colombia, se tiene
actualmente un Atlas de Radiación Solar de Colombia 22, que contiene 13 mapas,
donde se representa la distribución espacial del potencial energético, uno por cada
mes, y uno adicional en donde se consigna el promedio anual. En los mapas se
encuentra el valor promedio diario de la radiación. Este atlas es fruto de un trabajo
en conjunto entre el Instituto de Hidrología, Meteorología y Estudios Ambientales
(IDEAM) y la Unidad de Planeación Minero Energética (UPME), quienes desde el
siglo anterior, han elaborado este trabajo, hasta la última versión que fue
publicada en el año de 2014.
Como se puede apreciar en la gráfica siguiente, se consigna el potencial
energético de algunas ciudades del país, dentro de esas, se encuentran las 4
ciudades principales, donde se podría pensar, que hay mayor probabilidad de
integrar la tecnología solar fotovoltaica a edificaciones, pues en estas ciudades, es
donde se tiene mayor desarrollo económico, social, tecnológico, etc.
22
UPME; IDEAM. Mapas de radiación solar global sobre una superficie plana. [en línea]. Ministerio
de Minas y Energía, 2014. [Consultado: 4 de abril de 2017]. Disponible en Internet:
[Link]
29
Figura 7. Potencial de energía solar, algunas ciudades capitales
Promedio de iradiación anual
6,0
kWh/m^2/día
5,0
4,0
3,0
2,0
1,0
0,0
Gracias al potencial de energía solar que se tiene en el territorio colombiano,
surge la necesidad de tener una estrategia que permita la inclusión de esquemas
de autogeneración, que contribuyan de manera directa a la reducción del impacto
negativo sobre el medio ambiente. Un ejemplo de ello, surgió en la ciudad de San
Francisco a comienzos de este año, que mediante el uso de la legislación, entro
en vigor una ley promovida por Scott Wiener, supervisor y miembro de comité
urbanístico y de medio ambiente de la ciudad de San Francisco, que busca exigir
a todos los edificios nuevos de más de diez pisos, destinar al menos el 15% del
área del techo, en la instalación de sistemas fotovoltaicos23. De esta manera, se
logra que las nuevas construcciones sean capaces de producir cierta parte de la
energía consumida. Estas aplicaciones de los sistemas fotovoltaicos, se les
denomina BAPV, pero antes proseguir con este concepto, es importante mirar la
regulación colombiana, quienes establecieron dos conceptos que permite clasificar
la generación local, que son: la autogeneración y cogeneración.
1.2.10 Autogeneración y cogeneración. Los sectores industrial y comercial en
Colombia, son los primeros en incursionarse en este tipo de tipo de soluciones,
aunque son pocas. Además, no se ha utilizado el poder legislativo para
implementar políticas que permitan consolidar una unión firme entre los sectores
de la construcción y la energía, que permitan implementar soluciones energéticas
como se hace en San Francisco. En el país, se comenzó a tener ciertas
aproximaciones con el concepto de generación distribuida en el artículo 5 de la
23
ALVAREZ, Raúl. San Francisco será la primera ciudad de EEUU en exigir paneles solares en
edificios nuevos [en línea]. Xataka. 20 de abril de 2016. [consultado: 7 de abril de 2017]. Disponible
en Internet: [Link]
paneles-solares-en-edificios-nuevos
30
Ley 1715 de 2014, en donde se definen los términos de autogeneración y
cogeneración. El primero lo definen como: “aquella actividad realizada por
personas naturales o jurídicas que producen energía eléctrica principalmente, para
atender sus propias necesidades. En el evento en que se generen excedentes de
energía eléctrica a partir de tal actividad, estos podrán entregarse a la red, en los
términos que establezca la Comisión de Regulación de Energía y Gas (CREG)
para tal fin”24. El segundo se define como: “producción combinada de energía
eléctrica y energía térmica que hace parte integrante de una actividad
productiva” 25
En un estudio realizado por un consorcio entre las empresas Hart y
Regeneración para la UPME, determinaron que la capacidad instalada mediante
el concepto de autogeneración para el 2014, fue de 1.193,1 MW, repartidos entre
el sector industrial con 234 MW, el sector comercial con 4,1 MW, y el sector
petróleo con 955 MW. Estas plantas de generación, se basan en tres fuentes, que
son: gas natural (56%), carbón (15%) e hidráulica (29%) 26. Para el caso de la
cogeneración en el mismo año, fue de 691,7 MW, siendo solo los sectores de
industrial y petróleo, los participantes. El industrial con una capacidad instalada de
596,7 MW y el sector petróleo con 95 MW. Para la cogeneración, se emplearon las
siguientes fuentes primarias: biomasa (38%), gas natural (40%), carbón (18%), y
crudo (4%). El 23% de las soluciones energéticas de autogeneración están en el
departamento de Bolívar, y lo sigue Antioquia con el 20%. Mientras que las plantas
de cogeneración, están instaladas mayoritariamente en el departamento del Valle
del Cauca, con cerca del 55% 27, debido a los ingenios azucareros de esta región,
que han visto en la cogeneración, una buena rentabilidad.
La UPME considera la autogeneración a pequeña y gran escala. Según la
resolución 281 de 2015, se considera autogeneración a pequeña escala, cuando
la capacidad instalada del sistema de generación, no superara 1 MW. Además, en
el artículo [Link].4.6 del decreto 348 del 1 marzo de 2017, menciona los
siguientes parámetros para ser considerado autogenerador a pequeña escala:
24
COLOMBIA. CONGRESO DE LA REPUBLICA. Ley 1715 del 13 mayo de 2014. Por medio de la
cual se regula la integración de las energías renovables no convencionales al sistema energético
nacional. [En línea]. upme. p 3. [Consultado: 20 de abril de 2017]. Disponible en internet:
[Link]
25
Ibid. Pag 3
26
UPME. Capacidad instalada de autogeneración [en línea]. Consorcio HART-RE. [Consultado: 21
de abril de 2017]. Disponible en internet:
[Link]
[Link]
27
UPME. Capacidad instalada de cogeneración 2014 [en línea]. Consorcio HART-RE. [Consultado:
21 de abril de 2017]. Disponible en internet:
[Link]
[Link]
31
• La potencia instalada debe ser igual o inferior al límite máximo determinado
por la UPME para la autogeneración a pequeña escala (1MW).
• La energía eléctrica producida por la persona natural o jurídica se entrega
para su propio consumo, sin necesidad de utilizar activos de uso del sistema de
transmisión regional (STR) y/o sistemas de distribución local (SDL).
• La cantidad de energía sobrante o excedente podrá ser cualquier porcentaje
del valor de su consumo propio.
• Los activos de generación pueden ser de propiedad de la persona natural o
jurídica o de terceros y la operación de dichos activos puede ser desarrollada por
los propietarios o terceros.
[Link] Aplicaciones de autogeneración. Cuando se emplea el término de
generación distribuida, se pueden identificar cuatro tipos de aplicaciones
básicas28, que son:
• Generación para carga base: Sistema de generación principal y continua,
interconectado a la red de distribución para operaciones de compra y venta de
energía.
• Generación para carga en picos: Sistema alterno de respaldo para reducir el
consumo desde la red de distribución en los periodos de mayor demanda.
• Generación aislada: Sistema de generación para poblaciones totalmente
aisladas del sistema de interconexión eléctrica.
• Generación para soporte de la red de distribución: Sistema de respaldo
empleado especialmente en empresas de alto consumo energético solo en ciertos
periodos del año, o para casos en los que se requiere elevar los niveles de
confiabilidad en el suministro eléctrico.
1.2.11 Building Attached/Adapted/Applied Photovoltaics (BAPV). El término
BAPV, hace referencia a la energía solar fotovoltaica adaptada a las edificaciones,
que ha empezado a tomar fuerza en los últimos tiempos en nuestro país, pues
empresas nacionales e internacionales han apostado por esquemas de
28
MANTILLA, Análisis del esquema de generación distribuida como una opción para el sistema
eléctrico colombiano. Op, Cit. 100.
32
autogeneración con algunos usuarios del sector industrial y comercial. De acuerdo
con Martínez, el concepto como tal, “consiste en añadir sistemas fotovoltaicos en
los edificios, que en la mayoría de los casos se instalan sobre una estructura fija
anclada en la terraza o techo de los edificios, y por lo general, afectan de manera
negativa la estética del edificio, pero con el beneficio de tener una generación de
energía eléctrica en sitio” 29.
En la literatura, se encuentran dos términos que vinculan la tecnología solar
fotovoltaica con el sector de la construcción: el primero es el BAPV, y el segundo
es el BIPV, que se refiere a Fotovoltaica Integrada a Edificaciones. Jelle 30
menciona que este último se diferencia por ser un elemento convencional de
construcción, pero que genera energía. Es decir, son elementos que tienen doble
funcionalidad, y ejemplo de ello, se puede ver reflejado en tejas, claraboyas,
fachadas, ventanales, entre otras. Mientras que para el primer término, el sistema
fotovoltaico se considera como un complemento adicional del edificio, por lo cual
no se considera como un elemento directamente relacionado con los aspectos
funcionales de las estructuras.
1.2.12 Situación actual BAPV en el Valle del Cauca. A continuación, se
muestran los 3 sistemas fotovoltaicos más grandes de la región del Valle del
Cauca:
• Universidad Autónoma de Occidente - Cali:
Esta universidad de la ciudad de Cali, tiene el sistema solar fotovoltaico de mayor
potencia instalado en una institución educativa. El proyecto se realizó bajo la
coordinación del grupo de investigación de energías GIEN, en alianza con la
empresa de energías del pacifico EPSA E.S.P. La instalación fotovoltaica está
distribuida en tres partes, en la azotea del edificio central, y los parqueaderos de
ala sur y ala norte, que sumados, tienen una potencia 150 kWp. Actualmente se
empezó la expansión de este sistema, con 250 kWp adicionales, siendo el sistema
solar fotovoltaico con mayor capacidad en una institución educativa, con 400 kWp.
29
MARTINEZ, Alberto. Energía solar fotovoltaica integrada en la edificación: modelizado y análisis
del efecto del sombreado en la captación de irradiación. [en línea]. Doctorado en innovación en
ingeniería de producto y procesos industriales. La Rioja. Universidad de la Rioja. Facultad: Escuela
Técnica Superior de Ingeniería Industrial. Departamento de Ingeniería Mecánica, 2016. p. 34.
[consultado: 12 de abril de 2017]. Disponible en Internet:
[Link]
30
JELLE, Bjorn. Building Integrated Photovoltaics: A Concise Description of the Current State of the
Art and Possible Research Pathways. [en línea]. MDPI, Energies. 21 de septiembre de 2016. p. 2.
[consultado: 9 de abril de 2017]. Disponible en internet: [Link]
vor
33
Figura 8. SFV Azotea edificio central – UAO.
Fuente: UNIVERSIDAD AUTONOMA DE OCCIDENTE (UAO). [Imagen]. El
sistema solar fotovoltaico de mayor potencia instalado en una institución educativa
en Colombia. Noticias. Santiago de Cali, viernes 8 de mayo de 2015. [Consultado:
11 de abril de 2017]. Disponible en Internet:
[Link]
• Alkosto Híper ahorro - Cali
En el 2013, Alkosto inauguró una planta de generación fotovoltaica en los techos
de los estacionamientos de sus visitantes. La instalación cuenta con 498 módulos
fotovoltaicos, que cubren un área de 750 metros cuadrados y una potencia
instalada de 114 kWp que se interconecta a la red, y además de ello, instaló
sistemas de aires acondicionados de alta eficiencia VRF, iluminación LED, sistemas
de seguridad CCTV, que representan sin lugar a duda, un gran aporte en cuanto a
la eficiencia energética. Este sistema fotovoltaico genera en promedio 430
kWh/día.
Figura 9. SFV Estacionamientos – Alkosto.
Fuente: EMPRESAS Y [Link]. [Imagen]. Alkosto, utiliza la generación
fotovoltaica en Cali. Noticias, 14 de junio de 2013. [Consultado: 11 de abril de
2017]. Disponible en Internet: [Link]
34
• Mercamio – Cali
Este sistema fotovoltaico se inauguró en el 2016, y está instalado al sur de la
ciudad en el supermercado Mercamio. Cuenta con 408 módulos de 260W, lo que
representa una potencia instalada de 106 kWp. El proyecto fue montado por la
empresa Tecmac ingeniería.
Figura 10. SFV Azotea – Mercamio.
Fuente: TECMAC INGENIERÍA. Proyectos [Imagen]. Mercamio. [Consultado: 11
de abril de 2017]. Disponible en Internet:
[Link]
35
2. CARACTERIZACIÓN
2.1 PROPUESTA METODOLOGICA PARA EL MODELADO Y SIMULACION
DE SISTEMAS FOTOVOLTAICOS
El proceso metodológico que se puede observar en la siguiente figura, plantea el
desarrollo de cuatro etapas que abarcan: el levantamiento de la información, el
modelado de la misma, validación del modelo y evaluación de la red.
Figura 11. Metodología para el modelado de SFV conectados en redes de
distribución
La primera etapa de la propuesta consiste en el levantamiento de la topología de
la red eléctrica, que para este caso será la red de la Universidad Autónoma de
Occidente. Esto incluye la caracterización del sistema fotovoltaico instalado en la
institución, y el comportamiento de las cargas del cliente, la cual será desarrollada
en este capítulo.
En la segunda etapa que corresponde al modelado, se busca plasmar en la
herramienta de ETAP, todos los datos conseguidos de la red eléctrica, que se
compone de generadores fotovoltaicos, transformadores, inversores, cargas, entre
otros elementos más, que permitirán tener una representación acercada al sistema
real.
36
El proceso de decisión que se observa en la figura 11, hace referencia al recurso
solar, que si bien puede haber un registro de medición de temperatura como lo fue
en el caso bajo estudio de la Universidad Autónoma de Occidente, debido al
piranometro instalado cuando entro en operación el sistema fotovoltaico de 150
kWp en el año de 2015. Este trabajo se desarrolló con los registros del año 2016.
Pero si no se tienen estos registros, la otra opción la brinda la herramienta, pues
ésta cuenta con una calculadora de irradiancia, que se presentara más adelante
en el capítulo 3 (figura 39), la cual brinda al usuario un valor instantáneo de
irradiancia (W/m^2) según sea la hora y la ubicación del sistema fotovoltaico.
La tercera etapa, que está consignada en el capítulo cuatro de este trabajo,
pretende validar que tan aproximado se encuentra en modelo construido en la
herramienta ETAP, con el sistema real instalado en la institución, y para ello fue
importante tener la caracterización de los sistemas fotovoltaicos. La validación del
mismo, se basó en análisis de correlación, ajustes lineales de las variables reales
y pronosticadas, que permitieron al final, determinar el error porcentual absoluto
medio de estas variables bajo análisis.
Y por último en la etapa cuatro, se busca establecer condiciones típicas de
operación del sistema eléctrico, que más se aproximen al modelo real. Por ello, a
la hora de realizar la caracterización en este capítulo, se tuvieron en cuenta por
separado, los días hábiles y los días festivos. Además de ello, es importante hacer
la salvedad, que debido a limitaciones de la herramienta, no es posible realizar un
análisis durante periodo largos, lo cual conllevo a realizar las simulaciones de la
red en estado estable. Más adelante el autor dará conocimiento acerca de las
consideraciones que se tomaron para representar de mejor manera el
comportamiento de la red, durante un semestre en la universidad.
2.2 SISTEMA FOTOVOLTAICO 150 KWP (FASE 1)
La siguiente caracterización del sistema fotovoltaico de la fase I y fase II, se
realizará con base a la norma internacional IEC 60364-7-712, que fue elaborara
por el comité técnico 64 de IEC: Instalaciones eléctricas y protección contra los
choques eléctricos. Esta norma se basa en las instalaciones eléctricas en edificios
y específicamente la parte 7-712 involucra las reglas para las instalaciones y
emplazamientos especiales de sistemas con alimentación solar fotovoltaica.
La universidad Autónoma de occidente es pionera a nivel nacional, en cuanto a
generación fotovoltaica. Pues desde el año 2015, se instaló el primer sistema
fotovoltaico de 150 kWp en la institución, convirtiéndose en la planta de
generación fotovoltaica más grande instalada en una institución educativa
37
2.1.1 Ubicación. La localización de este sistema fotovoltaico es en la ciudad de
Cali, Colombia: 3°21’12.4” Latitud Norte y 76°31’21.2 Longitud Oeste. Las
coordenadas están expresadas para esta zona en grados-minutos-segundos. Para
un mayor entendimiento, estas también se pueden expresar en grados decimales:
3.53449° Norte, -76.522554° Oeste. La siguiente imagen es la vista superior del
sitio de emplazamiento de los tres subsistemas que componen el sistema
fotovoltaico de 150 kW de la Universidad Autónoma de Occidente de Cali.
Figura 12. Vista superior del campus universitario
Fuente: Google Maps. [Imagen]. Vista superior Universidad Autónoma de
Occidente. [Consultado el 18 de Agosto de 2017].
2.1.2 Objeto. La Universidad Autónoma de Occidente desde hace unos pocos
años, una de las metas a largo plazo, es la de tener un campus sostenible, y bajo
esa idea y con ayuda de la empresa EPSA E.S.P, el sistema fotovoltaico se
convirtió en una realidad. Capaz de generar un promedio mensual alrededor de los
600 kWh/mes, valor que vendría siendo aproximadamente el 5% de la demanda
mensual de la universidad.
2.1.3 Características módulos fotovoltaicos marca Panasonic.
[Link] Panasonic VBHN240. En las tablas siguientes, se consignan las
especificaciones eléctricas y mecánicas, de los módulos de referencia VBH240,
utilizados para el sistema.
38
Tabla 4. Especificaciones eléctricas módulo VBHN240
Parámetro Valor
Pmax 240W
Vpm 43,7 V
Ipm 5,51 A
Voc 52,4 V
Isc 5,85 A
Coef. Temp (Pmax) -0,3%/°C
Coef. Temp (Voc) -0,126 V/°C
Coef. Temp (Isc) 1,76 mA/°C
NOCT 48,3 °C
CEC PTC Raiting 223,5 W
Eficiencia de celda 21,60%
Eficiencia de módulo 19%
Watts/Ft^2 17,7 W
Voltaje max del sistema 600 VDC
Capacidad del fusible serie 15 A
Tolerancia de garantía (-/+) 0% / +10 %
Tabla 5. Especificaciones mecánicas módulo VBHN240
Parámetro Valor
Diodos de bypass 3
Área del módulo 1.26 m^2
Masa 15 kg
Dimensiones (LxWxH) 1580x798x35 mm
Longitud del cable 1030/880 mm
Tamaño del cable/tipo No. 12 AWG/PV cable
Tipo de conector Multi-Contact tipo IV (MC4)
Carga estática del viento y nieve 2.400 Pa.
[Link] Panasonic VBHN230. En las tablas siguientes, se consignan las
especificaciones eléctricas y mecánicas, de los módulos de referencia VBH230,
utilizados para el sistema.
39
Tabla 6. Especificaciones eléctricas módulo VBHN230
Parámetro Valor
Pmax 230W
Vpm 42,3 V
Ipm 5,45 A
Voc 51,2 V
Isc 5,83 A
Coef. Temp (Pmax) -0,3%/°C
Coef. Temp (Voc) -0,126 V/°C
Coef. Temp (Isc) 1,76 mA/°C
NOCT 48,3 °C
Eficiencia de celda 21 %
Eficiencia de módulo 17,9 %
Voltaje max del sistema 600 V
Capacidad del fusible serie 15 A
Tabla 7. Especificaciones mecánicas módulo VBHN230
Parámetro Valor
Diodos de bypass 3
Área del módulo 1.28 m^2
Masa 15 kg
Dimensiones (LxWxH) 1580x812x35 mm
Los módulos fueron instalados en diferentes configuraciones dependiendo de la
disponibilidad de área, y la capacidad de los inversores. En la siguiente sección,
se precisará los esquemas de conexión que fueron empleados para la integración
del sistema fotovoltaico.
2.1.4 Tipos de esquemas - Sistema fotovoltaico 150 kW
[Link] Parqueadero ala norte. Esta zona tiene 188 módulos fotovoltaicos de
240Wp de la marca PANASONIC, tecnología HIT, referencia VBHN240.
Adicionalmente, estos módulos fotovoltaicos se conectarán a dos inversores SMA
STP20000TL-US cada uno de potencia nominal de 20 kW. Esta zona de parqueo
tendrá una potencia instalada de 45.12 kWp. Los 188 módulos se dividirán en dos
secciones, cada una de 94 módulos. Cada sección estará compuesta por cinco
cadenas de 10 módulos conectados en serie y cuatro cadenas de 11 módulos
conectados en serie de acuerdo con el siguiente diagrama esquemático:
40
Figura 13. Diagrama unifilar SFV ala norte
Fuente: Grupo de Investigación de Energías (GIEN - UAO). Distribución del
sistema solar fotovoltaico del campus Valle de Líli Universidad Autónoma de
Occidente.
[Link]. Parqueadero ala sur. Esta zona cuenta con 225 módulos fotovoltaicos
divididos en dos secciones: una de 110 módulos de 240Wp de la marca
PANASONIC, tecnología HIT, referencia VBHN240 y otra de 115 módulos de
230Wp del mismo fabricante y tecnología, referencia VBHN230. Adicionalmente,
estos módulos fotovoltaicos están conectados a dos inversores SMA
STP24000TL-US cada uno de potencia nominal de 24 kW. Esta zona de parqueo
tendrá una potencia instalada de 52.85 kWp y tendrá el siguiente diagrama
esquemático:
41
Figura 14. Diagrama unifilar SFV ala sur
Fuente: Grupo de Investigación de Energías (GIEN - UAO). Distribución del
sistema solar fotovoltaico del campus Valle de Líli Universidad Autónoma de
Occidente.
[Link]. Azotea edificio administrativo. Esta zona tiene dispuesto 230 módulos
fotovoltaicos de 230Wp de la marca PANASONIC, tecnología HIT, referencia
VBHN230. Adicionalmente, estos módulos fotovoltaicos están conectados a dos
inversores SMA STP24000TL-US cada uno de potencia nominal de 24 kW. Esta
zona de azotea cuenta con una potencia instalada de 52.9 kWp y posee el
siguiente diagrama esquemático:
42
Figura 15. Diagrama unifilar SFV azotea
Fuente: Grupo de Investigación de Energías (GIEN - UAO). Distribución del
sistema solar fotovoltaico del campus Valle de Líli Universidad Autónoma de
Occidente.
2.1.5 Protecciones. En las figuras de la sección anterior, se puede observar que
para cada una de las cadenas que componen los arreglos de módulos, se tiene
una protección contra sobre corrientes mediante un fusible de 15A en el lado DC,
el cual se encuentra ubicado en las cajas combinadoras. Para tener mayor
claridad, a continuación se presenta una figura permite conocer interior de una
caja combinadora. Para este sistema fotovoltaico se emplearon 12.
Figura 16. Disposición caja combinadora
43
En la figura anterior se pueden observar los fusibles que protegen a cada uno de
las cadenas de módulos, ubicados en la parte superior izquierda e inferior
derecha. En la parte superior derecha se encuentra ubicado el dispositivo de
protección contra sobretensiones (DPS) en el lado de DC. En el lado AC, existe
protección contra sobre intensidades; en total son seis, tres de 30A y tres de 40A.
Además, cada uno de los módulos cuenta con tres diodos de by-pass con el fin de
permitir el flujo de corriente por las celdas sin ningún problema, cuando se tiene
algún sombreado sobre el elemento. También es de importancia mencionar que
estos dos tipos de módulos, están bajo la conformidad de la norma IEC61215.
Según el voltaje máximo de los módulos es de 600V, razón por la cual, cada una
de las cadenas del sistema fotovoltaico no supera ni los 500V.
2.1.6. Radiación. En la ubicación de este sistema fotovoltaico, hay un gran
potencial del recurso solar, pues gracias al piranómetro que se instaló en la azotea
de la universidad, se cuenta con un registro de datos, que permiten evaluar dicho
recurso. Para este trabajo en particular, se analizaron los datos para el año 2016,
que se extrajeron mediante el portal de Sunny, que es la marca de los inversores.
A continuación se mostrará la radiación promedio diaria que se tuvo para cada uno
de los 12 meses del año.
Figura 17. Irradiación diaria promedio
Iradiación diaria promedio
6,00
kWh/m^2 / dia
5,00
4,00
3,00
2,00
1,00
-
De la figura anterior se puede afirmar que existe una enorme potencial de energía
por explotar. Fruto de ello, es la expansión de este sistema fotovoltaico, que para
mediados del 2017, incremento su potencia instalada a 400 kW. Según los datos
obtenidos, la radiación promedio diaria en el año de 2016, es de 4,97 kWh/m^2.
44
2.1.7. Rendimiento. En esta sección se presentara el rendimiento que tuvo el
sistema fotovoltaico para el año 2016. En la siguiente figura, se observará la
energía diaria promedio que entrego el sistema fotovoltaico de 150 kWp en los
doce meses del año en mención.
Figura 18. Generación diaria promedio mensual
Generación diaria promedio de energía
700
600
500
kWh/dia
400
300
200
100
-
Las variables generación y radiación tienen una relación proporcional, razón por la
cual se evidencia la disminución de energía para el último mes del año. Para ese
mismo mes, la radiación diaria promedio fue de 3,92 kW/m^2. Mientras que para el
mes de marzo, donde se tuvo una mayor generación de energía, la radiación
promedio diaria fue de 5,16 kW/m^2. La producción de energía en el mes de
marzo, fue de 18,833 MW/mes y la producción total del año, fue de 208,56
MWh/año.
2.3 SISTEMA FOTOVOLTAICO 250 KWP (FASE 2)
A principios de julio de 2017, entro en operación la segunda fase del proyecto
fotovoltaico de la Universidad Autónoma de Occidente. Con la expansión, el
sistema fotovoltaico sumo un total de 1.546 módulos, quedando el sistema
fotovoltaico con una capacidad instalada de 402 kWp31.
31
CELSIA. Celsia pone en operación segunda fase de techos solares en Universidad Autónoma de
Occidente [En línea]. Sala de prensa, 30 de junio de 2017. [consultado: 20 julio de 2017].
Disponible en internet: [Link]
autonoma
45
2.2.1. Ubicación. La expansión del sistema fotovoltaico se compone de 6
subsistemas, uno por cada bloque de aulas, uno en el edificio de bienestar
universitario y por último, el subsistema de del edificio de la cafetería. Se podría
considerar la misma localización del sistema de 150 kWp que se mostró
anteriormente, pues las coordenadas para la zona más lejana que es la cafería,
solo varia en unas cuantas milésimas (3.53282° N, -76.520729° O). En la siguiente
imagen se muestra el sitio de emplazamiento de los subsistemas.
Figura 19. Emplazamiento sistema fotovoltaico fase II
Fuente: CELSIA. Celsia pone en operación segunda fase de techos solares en
Universidad Autónoma de Occidente [Imagen]. Sala de prensa, 30 de junio de
2017. [Consultado: 20 julio de 2017]. Disponible en internet:
[Link]
2.2.2 Objeto. La institución en vista de los beneficios ambientales y económicos
que se obtuvieron de la primera experiencia, se decidió junto con su aliado EPSA,
aumentar su capacidad instalada, que pasó de atender el 5% al 14,8% de la
demanda de energía del campus; con lo que se evitará la emisión de 221.8
toneladas de CO2. 32
32
Ibíd., párr. 2. Disponible en internet: [Link]
universidad-autonoma
46
2.3.3 Características módulo fotovoltaico marca Hareon
[Link] Hareon HR-265W. En las tablas siguientes, se consignan las
especificaciones eléctricas y mecánicas, de los módulos de referencia HR-265W,
utilizados para el sistema nuevo.
Tabla 8. Especificaciones eléctricas módulo HR-265W
Parámetro Valor
Pmax 265W
Vpm 30,71 V
Ipm 8,63 A
Voc 37,81 V
Isc 9,24 A
Coef. Temp (Pmax) -0,44%/°C
Coef. Temp (Voc) -0,32%/°C
Coef. Temp (Isc) 0,055%/°C
NOCT 45°C ± 3 °C
Temperatura de operación - 40°C ~ +85°C
Eficiencia de módulo 16,33%
Voltaje max del sistema 1000 VDC
Capacidad del fusible serie 15 A
Tolerancia de garantía (-/+) 0 ~ +5 W
Tabla 9. Especificaciones mecánicas módulo HR-265W
Parámetro Valor
Tipo de celda Poli-cristalina
Número de celdas 60 (6x10)
Masa 19 kg
Dimensiones (LxWxH) 1636x992x35 mm
Tamaño del cable 4mm^2
No. de diodos By-pass 3
Para la expansión del sistema fotovoltaico, se utilizó solo una referencia de
modulo, lo que brinda cierta simetría a la hora de realizar la configuración de los
esquemas de conexión. En el siguiente apartado, se muestra de manera detallada,
la configuración de los subsistemas.
47
2.3.4 Tipos de esquemas - Sistema fotovoltaico 250 kW
[Link] Edificio - Aulas. En la institución actualmente existen 4 bloques en donde
se disponen aulas de clases. Cada uno de estos bloques cuenta con 144 módulos
fotovoltaicos de 265 W marca Hareon de referencia HR-265W conectados a dos
inversores Fronius Symo con una potencia nominal de 17,5 kW. La potencia
instalada para cada edificio es de 38,16 kWp, y la configuración de cada
subsistema está compuesta por cuatro cajas de combinación, dos para cada
inversor. Cada una de las cajas de combinación, tiene conectado 36 módulos
fotovoltaicos dispuestos en dos cadenas, cada una de ellas con 18 módulos en
serie. A continuación se puede observar de mejor manera la conexión para cada
subsistema.
Figura 20. Diagrama unifilar SFV Edificio - aulas
Fuente: Empresa De Energía Del Pacífico S.A. E.S.P. (EPSA)
[Link] Edificio bienestar universitario y cafetería. Al extremo del campus se
sitúan dos subsistemas, uno para el edificio de bienestar universitario, y el otro
para el edificio donde se encuentra la cafetería. El techo de este último, soporta
180 módulos, y el edificio de bienestar, soporta 188, para un total de 368 módulos
fotovoltaicos, obteniendo así, una potencia instalada de 97,52 kWp, repartidos en
4 inversores, dos para cada edificio de marca Fronius Symo de 24 kW. A cada
inversor se conectan dos cajas de combinación, cada una permite construir un
arreglo de 46 módulos dispuestos en dos cadenas, cada una con 23 módulos
conectados en serie.
48
Figura 21. Diagrama unifilar SFV edificio bienestar universitario y cafetería.
Fuente: Empresa De Energía Del Pacífico S.A. E.S.P. (EPSA)
2.2.5 Protecciones. En los esquemas que se presentaron justamente en la
sección anterior, se puede observar que para cada una de las cadenas que
componen los arreglos, se tiene una protección contra sobre corrientes mediante
un fusible de 15A en el lado DC, el cual se encuentra ubicado en las cajas
combinadoras. En esa misma caja de combinación se tiene un interruptor de 25
ADC a 1.000 VDC; esto aplica para todas las cajas combinadoras. El lado DC
también cuenta con protección contra sobretensiones, y para ello, se ubicó un
DPS en cada extremo, es decir, justo después de salir de la caja de combinación,
justo antes de entrar al inversor, debido a la extensa distancia que había entre
estos dos. Por otro lado, en la parte AC, para cada uno de los inversores Symo de
17,5 kW (Aulas), se ubicó un interruptor de 30A por fase, mientras que para los
inversores de 24 kW (Cafetería y Bienestar Universitario), se tienen interruptores
de 40 A por Fase.
Según el voltaje máximo permitido por los módulos en una cadena, es de
1.000VDC, razón por la cual, cada una de las cadenas del sistema fotovoltaico
está por debajo de los 750VDC. También es importante mencionar, que al igual
que los módulos Panasonic, estos también están bajo conformidad de la norma
IEC61215.
2.2.6 Rendimiento. Para el rendimiento de la planta de generación de la fase I, se
presentó la producción de energía promedio diaria que se tuvo para los doce
meses del año 2016. En esta ocasión, se debe tener en cuenta que la instalación
fotovoltaica entro en operación en el mes de julio de 2017. Por tal razón, la
información que se presentara a continuación, corresponde al promedio de
producción diaria de energía que tuvo cada subsistema de la fase II en el segundo
semestre del año 2017.
49
Figura 22. Generación diaria promedio SFV fase II – Segundo semestre 2017
Fuente: Empresa De Energía Del Pacífico S.A. E.S.P. (EPSA)
El arreglo fotovoltaico que mayor energía aporta al día, es aquel que se encuentra
ubicado en los techos del edificio de bienestar universitario (B. U.) aportando un
promedio semestral cercano a los 200 kWh/día. Los subsistemas de las aulas uno,
dos y tres, generan al día, alrededor de 150 kWh/día. Mientras que el arreglo que
menos energía aporta diariamente es el que está ubicado en aulas 3, pues cuenta
con un promedio semestral diario de 132 kWh/dia. Adicionalmente en la siguiente
tabla, se consigna la suma total de energía generada para cada uno de los meses
ya mencionados.
Tabla 10. Generación total de energía SFV fase II – Segundo semestre 2017
Generación de energía mensual - kWh/mes
SFV\Mes Julio Agosto Septiembre Octubre Noviembre Diciembre
Aulas 1 5.017,17 4.986,24 5.019,80 4.370,51 4.433,89 4.445,70
Aulas 2 4.903,25 4.901,28 4.992,91 4.386,86 4.490,81 4.485,56
Aulas 3 3.672,23 3.306,13 4.311,56 4.398,19 4.518,78 4.134,59
Aulas 4 4.968,22 4.861,72 5.017,25 4.396,51 4.527,56 4.505,58
B. U. 6.404,78 6.208,57 6.450,63 5.672,24 5.771,24 5.827,30
Cafeteria 6.331,44 6.094,10 6.295,23 5.551,71 5.578,84 5.561,94
Total 31.297,09 30.358,04 32.087,38 28.776,02 29.321,12 28.960,67
50
La mayor producción de energía se obtuvo en el mes de septiembre, en donde la
institución tuvo una generación de energía igual a 32.087 kWh/mes, mientras que
la más baja, se obtuvo para el mes de octubre con una generación de 28.776
kWh/mes, es decir, 3.311 kWh menos que para el mes de septiembre. Al finalizar
el año 2017, la generación de energía fue de 180.800 kWh/semestre.
2.4 CARACTERIZACIÓN DEMANDA GENERAL - RED ELÉCTRICA UAO
En este apartado, se presentará una caracterización general de la demanda de la
Universidad Autónoma de Occidente. La red eléctrica de la institución, está
constituida por tres subestaciones, cada una de ellas, alimentada por
transformadores de 500, 1000 y 800 kVA. El transformador 1 de 500 kVA,
empleado en la subestación uno, alimenta las cargas del edificio central, los
laboratorios ubicados en los sótanos uno y dos; El transformador dos de 1000
kVA, empleado para la subestación dos, contiene las cargas de las aulas, la
cafetería, bienestar universitario, la planta de tratamiento de aguas residuales
(PTAR), y la planta de tratamiento de agua potable (PTAP), siendo estas, las
cargas de mayor representación conectadas a dicho nodo; para terminar, el
transformador tres se encarga de alimentar la carga conectada a la subestación
tres, la cual representa a las unidades enfriadoras o chillers.
Este trabajo no tendrá como prioridad, detallar las cargas que se conectan aguas
debajo de las subestaciones, pues se considerara una sola carga general, que
represente cada uno de los circuitos ramales, que se conectan a dichos nodos.
Estas cargas se determinaron mediante los datos registrados por los medidores
Scheneider Electric instalados en cada una de las subestaciones de la
universidad. La obtención de la información fue posible gracias a una plataforma
de dicho fabricante, que tiene como nombre PowerMonitor. La información
extraída corresponde a los primeros siete meses del año 2017, con la cual se
pudo constatar la diferencia que se tiene con la demanda de los días entre
semana, y los días festivos (sábado y domingo). A continuación, se presenta
mensualmente el promedio diario de dichos comportamientos.
51
Figura 23. Diagrama unifilar UAO
2.3.1 Enero. En el primer mes del año, se tuvo un comportamiento esperado, pues
fue el tercer consumo de energía más bajo del semestre. Este comportamiento
radica en el periodo vacacional que se tuvo en las dos primeras semanas. Cuando
se iniciaron clases, la demanda comenzó a estabilizarse. En la siguiente figura se
presenta el perfil diario promedio de demanda para cada una de las
subestaciones, teniendo en cuenta tanto los días hábiles, como los festivos.
Figura 24. Demanda diaria promedio – Enero
52
La demanda entre semana de la institución difiere totalmente de los consumos de
los fines de semana. Por tal razón, en las gráficas de cada mes, se visualizarán
curvas con el nombre ES y FDS. El primero significa Entre semana y el segundo
corresponde a Fines de semana; la palabra Sub hace referencia a la subestación,
que están indicadas con el número uno, dos y tres. Para este caso en particular,
se obtuvo una mayor demanda de energía en la subestación tres en las horas del
mediodía, y es de esperarse, pues en esas horas es cuando los chillers realizan
un mayor esfuerzo por mantener una temperatura adecuada en cada una de las
zonas de la universidad. Normalmente, en horas del mediodía, la subestación dos
es aquella que exige mayor energía, debido a que de allí se alimentan las cargas
de la cafetería. La demanda diaria promedio para las subestaciones uno, dos y
tres, en un día ES para el mes de enero, es de 156.58, 143.97, 134.46 kWh para
la Sub1, Sub2 y Sub3, respectivamente. Mientras que para un día FDS es 65.58,
83.51 y 17.17 kWh, en ese mismo orden. Adicionalmente, también es importante
tener una diferenciación de los consumos dependiendo de la jornada durante el
día. Por ello, la siguiente tabla permite presentar los promedios horarios de
consumo para las tres jornadas que se manejan en la institución. Estas se
clasifican en: Jornada diurna, que va desde las 6:00 hasta las 17:00; la jornada
nocturna está comprendida desde las 18:00 hasta las 21:00; y la jornada de la
madrugada se tiene desde las 22:00 hasta las 05:00 horas. En esta ocasión, solo
se analizará los promedios del consumo total.
Tabla 11. Promedios de consumo por jornada - Enero
SUB 1 SUB 2 SUB 3
TOTAL
JORNADA [kVA] [kVA] [kVA]
ES FDS ES FDS ES FDS ES FDS
Diurna
195,76+52,95i 61,49+20,69i 176,77+58,63i 90,96+24,7i 214,01+54,6i 26,01-25,4i 586,54+166,18i 178,46+19,99i
(6:00 - 17:00)
Nocturna
204,68+48,62i 77,82+14,78i 165,05+43,13i 78,13+17,16i 184,22+42,49i -0,15-36,94i 553,95+134,24i 155,8-5i
(18:00 - 21:00)
Madrugada
73,59+16,34i 65,59+15,32i 84,24+18,9i 75,03+16,12i -0,15-36,99i -0,15-37,1i 157,68-1,75i 140,47-5,66i
(22:00 - 05:00)
Según el promedio total de las tres subestaciones, la jornada en la que mayor
consumo de energía se tiene es para la diurna, con 609,62 kVA en un día entre
semana (ES). El consumo promedio para un día festivo para la misma jornada,
representa el 29,4% del consumo promedio entre semana. Mientras que para la
jornada nocturna, el consumo promedio de un día ES fue de 569,98 kVA y para un
día FDS es tan solo de 155,88 kVA, es decir, el 27% de un día entre semana. Para
la última jornada que corresponde a la madrugada, la variación no es tan
representativa, pues el promedio que se obtuvo para la misma fue de 157,68 kVA
para ES y 140,58 kVA para FDS, teniendo solamente un desfase de
aproximadamente 17 kVA.
53
2.3.2 Febrero. Ya en este mes, los consumos de energía se comenzaron a
normalizar debido al inicio de clases del primer semestre del año. De hecho, el
pico de consumo se consigue por la demanda de la subestación dos después del
mediodía para los días entre semana. Incluso, para el promedio de los días FDS
de este mes, se obtuvo un consumo cercano a los 191 kW en horas del almuerzo,
que representa un 62% aproximadamente del consumo máximo que se tuvo en el
promedio de días ES de la subestación dos. La demanda diaria promedio en el
mes de febrero para las subestaciones uno, dos y tres, en un día ES, es de
157.25, 196.56, 107.46 kWh, respectivamente. Mientras que para un día FDS es
75.33, 119.16 y 33.8 kWh, respectivamente.
Figura 25. Demanda diaria promedio – Febrero
Los comportamientos para cada jornada, se muestran a continuación.
Tabla 12. Promedios de consumo por jornada - Febrero
SUB 1 SUB 2 SUB 3
TOTAL
JORNADA [kVA] [kVA] [kVA]
ES FDS ES FDS ES FDS ES FDS
Diurna
194,44+66,18i 78,04+32,16i 250,65+96,71i 152,22+55,72i 167,5+38,08i 67,09-11,23i 612,59+200,97i 297,35+76,65i
(6:00 - 17:00)
Nocturna
212,13+59,74i 81,27+21,68i 240,55+78,66i 89,33+24,62i 147,65+23,45i 1,21-37,22i 600,33+161,85i 171,81+9,08i
(18:00 - 21:00)
Madrugada
75,53+19,42i 67,47+18,04i 94,04+24,66i 84,49+22,99i 1,11-37,31i 1,03-37,4i 170,68+6,77i 152,99+3,63i
(22:00 - 05:00)
Al igual que el mes anterior, el comportamiento se sigue manteniendo, pues la
jornada que mayor consumo de energía tuvo, es la diurna, con 644,71 kVA en un
54
día ES. El consumo promedio para un día festivo para la misma jornada,
representó en este caso, el 47,63% del consumo promedio entre semana, es decir
18.23 puntos más que el mes anterior. Para la jornada nocturna, el consumo
promedio de un día ES fue de 621,76 kVA y para un día FDS es tan solo de
172,04 kVA, es decir, el 27,66% de un día entre semana. Para la última jornada
que corresponde a la madrugada, la variación tampoco fue representativa, pues
solo están separados por 17,78 kVA, siendo el día ES de mayor consumo.
2.3.4 Marzo. Para el tercer mes del año, ya se han nivelado los consumos de la
universidad, razón por la cual, los perfiles de demanda diarios empiezan a tener
cierta simetría entre sí. El consumo promedio diario del mes de marzo, fue el más
alto de los siete meses bajo análisis, siendo de 534.9 kWh/día. Además, al igual
que el caso anterior (febrero), el pico de demanda de la subestación dos, se
obtiene en horas del mediodía. Y la demanda diaria promedio en el mes de marzo
para las subestaciones uno, dos y tres, en un día ES, es de 196.08, 227.27,
111.55 kWh, respectivamente. Mientras que para un día FDS es 72.72, 115.9 y
26.50 kWh, respectivamente. A continuación, se presenta el comportamiento diario
promedio que se obtuvo para este mes.
Figura 26. Demanda diaria promedio – Marzo
Los comportamientos para cada jornada, se muestran a continuación
55
Tabla 13. Promedios de consumo por jornada - Marzo
SUB 1 SUB 2 SUB 3
TOTAL
JORNADA [kVA] [kVA] [kVA]
ES FDS ES FDS ES FDS ES FDS
Diurna
253,48+66,9i 77,04+32,65i 287,56+107,95i 145,38+52,98i 179,3+50,15i 52,2-17,47i 720,34+225i 274,62+68,16i
(6:00 - 17:00)
Nocturna
262,61+56,41i 76,34+18,38i 291,16+95,25i 89,34+24,52i 129,47+18,27i 1,37-36,88i 683,24+169,93i 167,05+6,02i
(18:00 - 21:00)
Madrugada
77,66+19,44i 64,99+18,52i 104,88+28,37i 86,24+23,19i 1,58-36,92i 1,46-36,98i 184,12+10,89i 152,69+4,73i
(22:00 - 05:00)
El mes de marzo es aquel donde se consiguieron los promedios más altos para las
jornadas diurnas y nocturnas. Para la primera se obtuvo un consumo de 720,69
kVA en un día entre semana (ES). El consumo promedio para un día festivo en la
misma jornada, representa el 39,2% del consumo promedio entre semana.
Mientras que para la jornada nocturna, el consumo promedio de un día ES fue de
704,05 kVA y para un día FDS es tan solo de 167,15 kVA, es decir, el 23,7% de un
día entre semana. Para la última jornada que corresponde a la madrugada, un día
ES tuvo un promedio de 184,44 kVA y 152,76 kVA para un día FDS, teniendo
solamente un desfase de aproximadamente 31,68 kVA.
2.3.4 Abril. Como era de esperarse, para el mes de abril, también se obtuvo un
comportamiento muy parecido a la demanda diaria promedio de energía de los
últimos dos meses presentados anteriormente (febrero y marzo). En este caso, se
obtuvo un pico máximo en la subestación dos de 322.36 kW a las 3:00 P.m.
mientras que para la subestación uno, alcanzo su máxima demanda a las 5:00
P.m. de 292.28 kW, para el promedio diario ES. Además, el consumo de energía
de las unidades enfriadoras que están conectadas a la subestación tres, se
aumentó en un 2,3% con respecto al mes de marzo, pasando de 111.55 kWh/día a
114.13 kWh/día. Por último, la demanda diaria promedio en el mes de abril para
las subestaciones uno, dos y tres, en un día ES, es de 179.29, 214.14 y 114.13
kW/h, respectivamente. Mientras que para un día FDS es 58.18, 105.72 y 25.81
kW/h, respectivamente.
56
Figura 27. Demanda diaria promedio – Abril
Los comportamientos para cada jornada, se consignan en la siguiente tabla.
Tabla 14. Promedios de consumo por jornada - Abril
SUB 1 SUB 2 SUB 3
TOTAL
JORNADA [kVA] [kVA] [kVA]
ES FDS ES FDS ES FDS ES FDS
Diurna
227,75+61,26i 57,09+27,98i 270,43+102,11i 128,66+43,65i 178,85+50,53i 48,14-21,01i 677,03+213,9i 233,89+50,62i
(6:00 - 17:00)
Nocturna
248,94+53,16i 64,8+15,16i 272,55+91,53i 85,72+21,41i 123,54+16,21i 7,78-35,3i 645,03+160,9i 158,3+1,27i
(18:00 - 21:00)
Madrugada
72,36+17,27i 57,26+16,3i 100,52+27,26i 81,32+20,95i 1,51-36,72i 1,43-36,9i 174,39+7,81i 140,01+0,35i
(22:00 - 05:00)
En este mes, el consumo promedio de energía para un día ES de la jornada diurna
fue de 710 kVA. El consumo promedio para un día festivo de la misma jornada,
representó en este caso, el 33,7% del consumo promedio entre semana. Para la
jornada nocturna, el consumo promedio de un día ES fue de 664,79 kVA y para un
día FDS es tan solo de 158,3 kVA, es decir, el 23,81% de un día entre semana.
Para la última jornada que corresponde a la madrugada, la variación fue por
aproximadamente 34 kVA, siendo el día ES de mayor consumo.
2.3.5. Mayo. En el quinto mes del año, y ulimo mes del periodo academico, se
obtuvo un comportamiento que no sucedió en ninguno de los meses anteriores, y
es la paridad que tuvieron las curvas referidas a ES Sub1 y ES Sub2. Dicho
comportamiento se produjo a raiz del solapamiento que comienza a existir para
57
estas fechas finales en los laboratorios de la universidad. El ultimo mes del
semestre es cuando los estudiantes recurren a realizar sus trabajos finales en las
salas o laboratorios de la institucion.
Figura 28. Demanda diaria promedio – Mayo
Los comportamientos para cada jornada, se consignan en la siguiente tabla.
Tabla 15. Promedios de consumo por jornada - Mayo
SUB 1 SUB 2 SUB 3
TOTAL
JORNADA [kVA] [kVA] [kVA]
ES FDS ES FDS ES FDS ES FDS
Diurna
245,32+60,95i 82,7+31,55i 244,44+93,7i 127,74+45,89i 183,62+56,17i 58,57-16,97i 673,38+210,82i 269,01+60,47i
(6:00 - 17:00)
Nocturna
255,02+53,21i 69,45+14,84i 280,6+72,59i 98,78+16,35i 147,55+26,1i 7,51-36,82i 683,17+151,9i 175,74-5,63i
(18:00 - 21:00)
Madrugada
77,51+18,73i 62,22+16,99i 108,81+13,98i 94,23+15,23i 2,35-37,04i 1,44-36,81i 188,67-4,33i 157,89-4,59i
(22:00 - 05:00)
En el mes de mayo se obtuvo un consumo promedio de 705,61 kVA en un día
entre semana, para un día festivo de la misma jornada, se tuvo un consumo
promedio de 275,72 kVA, estando cerca de 430 kVA por debajo del día ES. En la
jornada nocturna, el consumo promedio de un día ES fue de 699,85 kVA y para un
día FDS es de 175,83 kVA, es decir, el 25,12% de un día entre semana. Para la
última jornada de la madrugada, un día ES tuvo un promedio de 188,71 kVA y
157,95 kVA para un día FDS, teniendo solamente un desfase de
aproximadamente 30,76 kVA.
58
2.3.6 Junio. Para el final del semestre, el comportamiento de la demanda se
redujo notablemente a raíz de la finalización del periodo académico en la
institución. Este comportamiento no se había presentado en ninguno de los meses
anteriores, normalmente siempre hubo cierta paridad entre la subestación uno y
dos en cuanto al consumo de energía eléctrica en horas del mediodía. El pico que
se puede observar en la subestación dos en horas de la tarde – noche, pueden ser
debidas a la iluminación. En la subestación tres, también disminuyó el promedio
diario en un 17% con respecto al promedio de consumo del mes de marzo para los
días ES. La demanda diaria promedio en el mes de junio para las subestaciones
uno, dos y tres, en un día ES, es de 150.48, 98.34, 95.63 kW/h, respectivamente.
Mientras que para un día FDS es 63.29, 62.28 y 23.64 kW/h, respectivamente.
Figura 29. Demanda diaria promedio – Junio
Los comportamientos para cada jornada, se consignan en la siguiente tabla.
Tabla 16. Promedios de consumo por jornada – Junio
SUB 1 SUB 2 SUB 3
TOTAL
JORNADA [kVA] [kVA] [kVA]
ES FDS ES FDS ES FDS ES FDS
Diurna
187+57,35i 58,51+27,59i 83,08+16,83i 45,49+7,78i 155,28+33,51i 47,6-18,46i 425,36+107,69i 151,6+16,91i
(6:00 - 17:00)
Nocturna
200,54+53,07i 79,64+19,73i 158,82+1,41i 83,07-16,29i 103,27+4,77i 1,31-36,74i 462,63+59,25i 164,02-33,3i
(18:00 - 21:00)
Madrugada
70,6+17,28i 62,28+16,62i 90,55-22,45i 77,07-25,46i 2,33-36,44i 1,49-36,84i 163,48-41,61i 140,84-45,68i
(22:00 - 05:00)
En este mes se obtuvo el consumo promedio más bajo de los primeros siete
meses del 2017, y fue de 438,78 kVA para un día entre semana; para un día
59
festivo de la misma jornada, se tuvo un consumo promedio de 152,54 kVA. En la
jornada nocturna, el consumo promedio de un día ES fue de 466,4 kVA y para un
día FDS es de 167,36 kVA, es decir, el 35,88% de un día entre semana. Para la
última jornada de la madrugada, un día ES tuvo un consumo promedio de 168,69
kVA y 148,06 kVA para un día FDS, teniendo solamente un desfase de
aproximadamente 20,62 kVA.
2.3.7 Julio. El último mes del periodo de vacaciones, fue una constante del mes
junio para cada una de las subestaciones. Para esta fecha, entro en operación la
fase II del proyecto fotovoltaico de la institución. Razón por la cual, se observa que
para la curva FDS Sub 2, en el rango de 10:00 a.m. a 2 p.m. se entregan
excedentes a la red, debido a que la demanda es inferior a la producción de
energía en ese momento.
Figura 30. Demanda diaria promedio – Julio
Los comportamientos para cada jornada, se consignan en la siguiente tabla.
60
Tabla 17. Promedios de consumo por jornada – Julio
SUB 1 SUB 2 SUB 3
TOTAL
JORNADA [kVA] [kVA] [kVA]
ES FDS ES FDS ES FDS ES FDS
Diurna
194,63+57i 50,64+26,57i 123,95+31,74i 21,75-3,64i 186,2+54,44i 56,33-14,68i 504,78+143,18i 128,72+8,25i
(6:00 - 17:00)
Nocturna
207,46+53,86i 66,54+15,57i 192,52+1,96i 85,39-18,05i 123,26+11,25i 4,93-34,74i 523,24+67,07i 156,86-37,22i
(18:00 - 21:00)
Madrugada
69,73+16,9i 58,74+16,82i 92,75-18,16i 82,02-22,55i 5,11-34,79i 5,08-34,82i 167,59-36,05i 145,84-40,55i
(22:00 - 05:00)
Para este último mes bajo estudio, el consumo promedio de energía para un día
ES de la jornada diurna fue de 524,69 kVA. El consumo promedio para un día
festivo de la misma jornada, representó en este caso, el 24,58% del consumo
promedio entre semana. Para la jornada nocturna, el consumo promedio de un día
ES fue de 527,52 kVA y para un día FDS fue de 161,21 kVA, es decir, el 30,56%
de un día entre semana. Para la última jornada que corresponde a la madrugada,
la variación fue por aproximadamente 20,05 kVA, siendo el día ES de mayor
consumo.
Finalmente, lo que se buscaba con el análisis de la información anteriormente
mostrada, era validar la similitud de los consumos de energía que tiene la
universidad durante un periodo académico, tanto en un día hábil de la semana,
como para los días festivos. Además, se quería tener en cuenta los
comportamientos que se tienen durante el periodo de vacaciones, momento donde
la red eléctrica realiza un menor esfuerzo para atender la demanda de la
institución. Para darle validez a esto, se hace pertinente presentar los perfiles
promedios de demanda diarios que se tuvieron para cada uno de los meses del
bajo análisis del año 2017.
61
Figura 31. Perfil de demanda mensual – Primer semestre 2017
La conclusión que se obtuvo con respecto a estos comportamientos de consumo
de energía eléctrica en la universidad, corresponde a la uniformidad del consumo
que se tiene tanto en el periodo académico, como en un periodo vacacional. Para
realizar un análisis del comportamiento de la red frente, no hay que olvidar que la
curva de generación de energía que se tiene con una instalación fotovoltaica está
dada en un rango de 6:00 Am. A 6:00 P.m. Lo que permitirá en el capítulo 5,
realizar varias condiciones de operación, bajo distintos comportamientos de
consumo, que representaran a los días hábiles y los días festivos tanto del periodo
académico, como del periodo de vacaciones.
62
3. DESCRIPCIÓN HERRAMIENTA ETAP PARA LA SIMULACIÓN DE
SISTEMAS FOTOVOLTAICOS
La importancia de las herramientas de simulación cobró importancia a medida que
los sistemas de transmisión y distribución expandían su rango de operación. En
vista a ello, surgió la necesidad por parte de las empresas de energía (operadores
de red), a tener una óptima preparación antes y durante la concepción de las
redes, e incluso a la hora de planear la expansión de los sistemas, estas
herramientas juegan un papel importante, pues permiten conocer en detalle las
nuevas condiciones de operación a las cuales estaría sometido el sistema
eléctrico de potencia. Los sectores industrial y comercial, no han sido ajenos a
esta realidad, pues estos también poseen su sistema eléctrico, con la única
diferencia de ser redes de menor escala.
Actualmente el país está pasando por una época de transición, que mediante
elementos legislativos, han permitido a empresas en su mayoría del sector
industrial y comercial, incluir en su red eléctrica sistemas de generación de
energía que permitan tener cierta independencia de la red y además, que dicha
fuente de generación sea un respaldo que permita generar ahorros. En Colombia,
según la ley existen dos maneras de tener generación propia, una es la
cogeneración y la otra es la autogeneración. La primera ha sido explotada
mayormente por ingenios azucareros que vieron en esta, una opción rentable que
permite mediante un proceso térmico, generar energía. Mientras que la segunda,
ha despertado gran curiosidad en el territorio, pues ha habido empresas
interesadas en explotar el recurso de sol, y por medio de la tecnología solar
fotovoltaica han logrado tener procesos de autoconsumo, que permiten por un
lado mejorar sus facturas de servicio, y por el otro, contribuir al mejoramiento del
medio ambiente con soluciones sostenibles.
Dada la proliferación del uso de sistemas fotovoltaicos en el país, se hace
pertinente la utilización de una herramienta que permita determinar la producción
de energía de un sistema solar fotovoltaico teniendo en cuenta la variabilidad del
recurso, y su comportamiento conectado a la red de distribución desde el punto de
vista de simulación de sistemas eléctricos de potencia. Bajo este contexto, El
análisis de flujos de potencia, es el estudio que se pone en manifiesto en este
trabajo. Este permite evaluar el comportamiento de los componentes del sistema
eléctrico de la Universidad Autónoma de Occidente, donde se podrán determinar
los voltajes en los nodos, la corriente, la potencia activa y reactiva de la red, bajo
ciertas condiciones operativas del mismo. El análisis de flujo de potencia, también
permite:
63
• Determinar el efecto sobre la red eléctrica de la incorporación de nuevos
circuitos de carga y generación, que den lugar una modificación de las condiciones
de operación del sistema.
• Planificar acciones de contingencias ante perdidas temporales de la red
principal.
• Establecer condiciones óptimas de operación del sistema eléctrico.
• Establecer puntos de compensación reactiva.
El estudio se realizará en el software de simulación ETAP, que cuenta con una
capacidad de detalle en los elementos de diseño, que lo hace diferenciador a otro
software. ETAP permite simular el comportamiento de fuentes de energía
renovable como la solar y la eólica, teniendo en cuenta la variabilidad del recurso
que se tiene a lo largo del día al emplear este tipo de tecnología. La posibilidad de
analizar el comportamiento de la integración de sistemas de energía solar
fotovoltaica con la red, es la finalidad de este trabajo, y para ello, es importante
empezar explicando la manera cómo funciona este tipo de tecnología en el
software.
3.1 SISTEMA FOTOVOLTAICO EN ETAP
3.1.1 Módulo fotovoltaico (PV). El módulo fotovoltaico ha comenzado a tomar
importancia en la red eléctrica, por la masificación que ha tenido en los últimos
años. Este tipo de tecnología convierte la radiación del sol en corriente continua
mediante la utilización de semiconductores y luego es convertida por corriente
alterna a través de los inversores. El elemento PV Array en ETAP, se utiliza para
representar los paneles individuales conectados tanto en serie como en paralelo.
A continuación se muestra el orden que se tiene antes de consolidar un arreglo
fotovoltaico, para obtener unos parámetros deseados en un sistema fotovoltaico.
Figura 32. Arreglo fotovoltaico
Fuente: ETAP Help 12.6. Photovoltaic (PV) Array. [Imagen]
64
Un generador fotovoltaico en ETAP, se puede constituir por un arreglo de módulos
que puede estar conectado tanto en serie como en paralelo. La versión 12.6 del
software cuenta con una librería corta, pero ETAP permite al usuario digitar los
parámetros más importantes de un módulo fotovoltaico, como lo son la potencia, la
tensión de circuito abierto, la corriente de cortocircuito, el voltaje, la corriente en el
punto de máxima potencia, la eficiencia, la temperatura de celda bajo condiciones
STC (NOCT) y hasta los coeficientes de variación por efectos de la temperatura
ambiente. Con estos datos de placa, el software construye las curvas corriente Vs
voltaje para cada caso en particular. En la siguiente figura, se puede apreciar una
curva corriente Vs voltaje de un módulo fotovoltaico. En ella se puede apreciar los
cuatro parámetros principales del funcionamiento de un módulo, en la parte
superior izquierda se encuentra la corriente de cortocircuito, que para este caso es
alrededor de 2,7 A, y ocurre cuando los terminales del módulo, están
cortocircuitados entre sí. El último punto de la parte inferior derecha hace
referencia a la tensión de circuito abierto (21 V) y ocurre cuando no se tiene
ninguna carga conectada a sus terminales, por lo tanto no hay circulación de
corriente. Por último se tienen los puntos correspondientes a Imp y Vmp, los
cuales permiten obtener el punto operación de máxima potencia. Es decir, es el
momento en donde el modulo fotovoltaico, entrega su mayor capacidad.
Figura 33. Curva I Vs V de un módulo fotovoltaico
Fuente: Software ETAP. Help 12.6. [Imagen] PV panel page.
La curva IV está determinada por los valores de voltaje y corriente, los cuales son
afectados de manera directa por la temperatura de las celdas, que dependiendo
las condiciones de radiación, y temperatura ambiente, su valor será modificado.
En el caso del voltaje, es inversamente proporcional con respecto a la temperatura
de celda, mientras que para la corriente, su valor aumentara a medida que la
65
temperatura de celda sea mayor. En la siguiente ecuación, se evidencia cómo la
temperatura de celda es afectada por los parámetros S (irradiancia) y Tamb
(temperatura ambiente), que varían durante el día, debidos a la aleatoriedad del
recurso y a las condiciones climatológicas.
𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁−20
𝑇𝑇𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 = 𝑇𝑇𝐴𝐴𝐴𝐴𝐴𝐴 + 80
∗ 𝑆𝑆 [4]
A continuación, se presentaran los recuadros que aparecen en ETAP, para que
mediante la interfaz, el usuario pueda consignar cada uno de los parámetros del
elemento fotovoltaico.
3.1.2 PV Panel page. Las especificaciones eléctricas de un módulo fotovoltaico
son definidas en el editor PV Panel, como se muestra a continuación:
Figura 34. Especificaciones de módulo fotovoltaico
Fuente: Software ETAP Help 12.6. Elemento PV array [Imagen]. PV Array
editor/PV panel page.
• Power: en esta casilla se ingresa la potencia nominal en Watts del módulo FV.
También se tiene la posibilidad de consultar los módulos de la librería que posee
el software. Los parámetros de estos módulos no son posibles modificarlos, pues
corresponden a los datos por defecto de los fabricantes.
66
• Tol. P: corresponde a la tolerancia del módulo fotovoltaico. Normalmente este
dato es dado por el fabricante, pero ETAP, no lo utiliza para ninguno de sus
cálculos, pues solo es utilizado con fines informativos.
• Vmp: como se ha venido mencionando, este dato de voltaje hace referencia a
la tensión cuando se tiene el punto de máxima potencia. El Vmp ingresado para un
solo módulo.
• Voc: corresponde a la tensión de circuito abierto del módulo.
• %Eff: este parámetro es calculado por ETAP, y representa la eficiencia del
módulo en porcentaje. Esta dado por la siguiente ecuación:
𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃
%𝐸𝐸𝐸𝐸𝐸𝐸 = 𝑊𝑊 [5]
𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴[𝑚𝑚2 ]∗𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼 𝑏𝑏𝑏𝑏𝑏𝑏𝑏𝑏 [ 2 ]
𝑚𝑚
• Imp: corresponde a la corriente en el punto de máxima potencia y se ingresa
en amperios.
• Isc: corresponde a la corriente de cortocircuito del módulo individual.
• %Fill Factor: muestra el factor de relleno del módulo FV calculado en
porcentaje. Se define como el área rectangular en una curva I Vs V. Este factor
debe ser mayor que 0.7, lo cual indica que el módulo está en buenas condiciones.
Se calcula mediante la siguiente ecuación:
𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃
FF = (𝑉𝑉𝑉𝑉𝑉𝑉)∗(𝐼𝐼𝐼𝐼𝐼𝐼) [6]
Coeficientes de ajuste de rendimiento: la temperatura afecta el rendimiento de
los módulos fotovoltaicos. La magnitud del Voc es inversamente proporcional a la
temperatura de celda. En la mayoría de las células de silicio cristalino la variación
del Voc con respecto a la temperatura es de -0,5%/°C. Las más eficientes tienen
una variación del -0,35%/°C. Para el caso de las células del silicio amorfo, su
variación esta entre -0,2 y -0,3 %/°C. Mientras que para el caso de la corriente,
esta aumenta proporcionalmente con respecto a la temperatura de celda, debidos
al aumento en el número de portadores generados térmicamente en la célula. Sin
embargo, este efecto es ligero, aproximadamente 0,065%/°C, para las células de
silicio cristalino y 0,09%/°C para las células de silicio amorfo. La siguiente figura
muestra las curvas I Vs V que se pueden tener para una célula de silicio cristalino
a diversas temperaturas. La corriente cambia con la irradiancia, y el voltaje con la
temperatura.
67
Figura 35. Diferentes curvas de I-V
Fuente: ETAP Help 12.6. PV panel page. [Imagen]. Parámetros de afectación de
curvas I-V.
• Alpha Isc: en este campo se ingresa el coeficiente de ajuste para la
corriente de cortocircuito.
• Beta Voc: en este campo se ingresa el coeficiente de ajuste para el voltaje
de circuito abierto.
• Delta Voc: se ingresa el coeficiente de ajuste para el voltaje de circuito
abierto. Pero para este caso, este coeficiente se utiliza para calcular el voltaje de
circuito abierto basado en diferentes niveles de irradiancia a la radiación base.
En el recuadro base, se registran los datos de temperatura e irradiancia en
condiciones STC, y NOCT.
• Temp: Se ingresa la temperatura en grados Celsius (°C) que es utilizada
por el fabricante para determinar la potencia nominal del módulo fotovoltaico. La
base predeterminada para la temperatura es de 25°C (Condiciones STC).
• Irrad: se introduce la irradiancia base en W/m^2, con la cual el fabricante
determina la potencia nominal del módulo fotovoltaico. Según las condiciones
STC, este valor corresponde a 1000 W/m^2.
• NOCT: Corresponde a la temperatura nominal de operación de la celda, y
se ingresa en °C.
En el botón Library, se puede acceder a una ventana que contiene una serie de
fabricantes de módulos fotovoltaicos. A continuación se muestra la interfaz que se
despliega.
68
Figura 36. Librería de módulos fotovoltaicos
Fuente: Software ETAP. Help 12.6. Elemento PV array [Imagen]. PV Array
editor/PV panel page/Library quick pick.
3.1.3 PV Array page. Las especificaciones eléctricas de un arreglo de paneles
fotovoltaicos son definidas en el editor PV Array, como se muestra a continuación:
Figura 37. Especificaciones del módulo fotovoltaico
Fuente: Software ETAP Help 12.6. Elemento PV Array. [Imagen]. PV Array
editor/PV Array.
69
PV Panel
• Watt/Panel: se refiere a la potencia de cada módulo fotovoltaico. Este campo
no puede ser modificado en esta ventana.
• # in series: en este campo se define el número de módulos que estarán
conectados en serie. Los módulos en serie, determinan el voltaje que tendrá el
arreglo.
• # in paralelo: aquí se establece la cantidad de módulos que se instalaran en
paralelo. Los módulos en paralelo, determinan la corriente que tendrá el arreglo.
A continuación, en la figura 16 se explica de mejor manera, las variantes que se
tienen en el arreglo, cuando se tiene una configuración serie o paralelo.
Figura 38. Conexión serie y paralelo
Fuente: Software ETAP Help 12.6. Elemento PV Array [Imagen]. PV Array page –
DC PV Array Editor.
PV Array (Total)
• # of Panels: aquí se determina la potencia total que se tendrá con el arreglo de
módulos conectados tanto en serie como en paralelo.
• Volts DC: en este campo, el software realiza el cálculo de la potencia nominal
que se tendrá con el número de módulos conectados en serie.
• kW DC: corresponde a la potencia total DC en kW y se calcula en base a la
cantidad de módulos fotovoltaicos que componen el arreglo.
70
• Amps DC: corresponde a la corriente nominal que se tendría en el arreglo, y
es determinada por el número de paneles que se conecten en paralelo.
• Generation Category: son las etiquetas que se dan a diferentes formas de
generación. Sus nombres pueden ser modificados, como lo es para este caso, en
donde se simulará por horas del día. De esta manera, se logra conseguir las
curvas de generación de un sistema fotovoltaico.
• Irradiance: este campo muestra la irradiancia en el módulo fotovoltaico y está
dado por W/m^2 y puede ser definida por el usuario, o puede ser obtenida
mediante la calculadora de irradiación del programa. La potencia del arreglo (MPP
kW), está determinada en función de este valor de irradiancia, y la temperatura
ambiente (Ta).
• Ta: en este espacio se define la temperatura ambiente en °C, a la cual está
sometido el arreglo fotovoltaico.
• Tc: corresponde a la temperatura de celda, la cual es calculada mediante la
ecuación [Link] valor de temperatura tiene una relación inversa con la eficiencia y
la potencia de entrega de los módulos en el arreglo.
• MPP kW: ETAP calcula la potencia máxima que entregara el arreglo de
módulos fotovoltaicos, teniendo en cuenta los tres últimos parámetros.
Irradiance Calculator
Esta calculadora de irradiancia permite obtener información de una ubicación
específica, que es ingresada por el usuario junto con la fecha y la hora. La
calculadora determinará la irradiancia teórica en W/m^2, a nivel del mar. A
continuación se muestra la ventana que se despliega al oprimir el botón
Irradiance Calc, y se describe cada uno de los campos que contiene.
71
Figura 39. Calculadora de irradiación en ETAP
Fuente: Software ETAP Help 12.6. Elemento PV Array [Imagen]. PV Array page –
DC PV Array Editor / irradiance calc.
Location information
• Latitude: se ingresa la latitud en grados. El norte del ecuador, se define como la
dirección positiva.
• Longitude: se ingresa la longitud en grados. El oeste del primer meridiano, se
define como la dirección positiva.
• Time zone: se introduce el tiempo universal coordinado (UTC), para las
coordenadas ingresadas.
• Local time: en este campo se selecciona la hora del día, para la cual se requiere
el valor de irradiancia correspondiente.
• Date: aquí se establece la fecha para la cual se requiere el cálculo de
irradiancia.
72
• Boton Calculate: al dar click en él, ETAP determina la posición solar y la
irradiancia para ese instante de tiempo, con los datos de entrada que se
mencionaron anteriormente.
Calculation
• Declination: corresponde al alguno aparente del sol que esta al norte o al sur
del plano ecuatorial de la tierra.
• Equation of Time: La ecuación del tiempo es la diferencia entre el tiempo solar
aparente y el tiempo solar medio medido en un instante dado en el mismo punto
de la Tierra. En cualquier momento dado esa diferencia es la misma en todas
partes.
• Solar altitude: es el ángulo de elevación del sol. Es decir, el ángulo entre la
dirección del centro geométrico del disco aparente del sol y el horizonte
• Solar Azimuth: Para un observador, el azimut solar es el ángulo medido en
sentido horario desde el norte hasta el plano vertical formado por el sol y la
ubicación del observador.
• Solar time: es el cálculo del paso del tiempo basado en la posición del sol en el
cielo. La unidad fundamental del tiempo solar es el día. Cuando el sol es visible,
un observador en cualquier longitud puede medir la posición del sol en el cielo y
calcular su ángulo de la hora, que se interpreta como hora local para ese
observador.
• Sunrise: representa el instante en que el borde superior del sol aparece por
encima del horizonte en el oriente.
• Sunset: corresponde a la desaparición diaria del sol debajo del horizonte en el
occidente como resultado de la rotación de la tierra.
• Air mass: la masa de aire es la cantidad de luz solar absorbida a través de la
atmosfera. Esta trayectoria se compara generalmente con una trayectoria vertical
directamente al nivel del mar, que se designa como AM = 1 AM.
73
• Irradiance: en este campo se consigna la irradiance calculada por ETAP, según
los datos de entrada que son ingresados en este despliegue.
3.1.4 Physical Page. La herramienta también tiene una ventana en donde se
consignan las medidas físicas del módulo fotovoltaico, como el largo, el ancho, la
profundidad y el peso. A continuación se muestra una imagen para dar mayor
claridad al despliegue.
Figura 40. Physical page
Fuente: Software ETAP Help 12.6. Elemento PV Array [Imagen]. PV Array /
Physical page – DC PV Array editor.
Hasta aquí, se han presentado las ventanas más importantes del elemento PV
Array, sin tener en cuenta el convertidor DC/AC. ETAP tiene otro elemento
adicional, el cual tiene integrado a él, un inversor, y el despliegue de la ventana
que contiene sus parámetros, se muestra a continuación.
3.1.5 Inverter page. En esta pestaña del elemento, el usuario puede visualizar las
especificaciones eléctricas del inversor que está asociado al módulo fotovoltaico.
Sí el usuario lo que desea es ingresar o modificar parámetros del inversor, debe
dar click en el botón Inverter Editor. A continuación se muestra el despliegue de
la ventana Inverter, que es contenida por el elemento PV Array.
74
Figura 41. Invertir page – PV Array editor
Fuente: Software ETAP Help 12.6. Elemento PV Array [Imagen]. PV Array editor /
Inverter page.
PV Array – Total Rated
En este conjunto, se muestran los parámetros de voltaje, corriente y potencia
totales, que se tendrán dependiendo del número de módulos fotovoltaicos
empleados en el sistema, y la configuración del diseño, que se establecen en la
pestaña PV Array.
Inverter
En esta parte, se consigna la información básica, con respecto a la potencia
DC/AC del elemento inversor. A continuación se describe cada campo.
• ID: este campo corresponde a un nombre que identifica el inversor que se está
utilizando para un arreglo de módulos fotovoltaicos.
• DC: representa la potencia nominal en kW del lado de corriente continua del
inversor. Junto a este campo de la misma fila, se encuentran otros 3. Estos
corresponden al nivel de tensión DC que maneja el equipo, la FLA (full load
75
current), corresponde a la corriente DC a plena carga en amperios. Por último, se
tiene el campo de eficiencia en porcentaje del inversor.
• AC: en la fila de corriente alterna, se tiene el primer campo correspondiente a
la potencia nominal de salida del inversor en kVA, seguido de la tensión en kV, y al
igual que la fila DC, tiene un campo para la corriente en amperios, y por último, se
tiene el campo para asignar el factor de potencia del equipo en porcentaje.
En la figura 19 se mostró la edición desde la ventana del elemento conjunto,
módulo-inversor, y sus parámetros son generales. Pero si se desea tener un
mayor nivel de detalle, se debe ingresar al Rating del elemento, pulsando el botón
Inverter Editor. La ventana que se despliega, se muestra a continuación.
Figura 42. Editor de elemento inversor
Fuente: Software ETAP Help 12.6. Elemento Inverter [Imagen]. Inverter editor /
Rating.
En la pestaña Rating se puede especificar el rango de tensión DC para su
operación en el conjunto DC Rating, la eficiencia del inversor se puede ingresar
según la cargabilidad del mismo. También se puede la contribución del
cortocircuito en el sistema de distribución, además de poder determinar el tipo de
76
conexión que se tiene a tierra. En las siguientes pestañas del elemento, se permite
configurar categorías de generación, establecer fuentes de armónicos, entre otros.
Volviendo a la pestaña inverter de la ventana del elemento PV Array, en la parte
inferior derecha, se encuentra ubicado un botón Cable library, este corresponde
al conductor que va desde la salida del arreglo fotovoltaico, hasta la entrada del
inversor. ETAP tiene una librería en donde se pueden encontrar diversos tipos de
conductores, a los cuales se les ingresa la longitud del cable que se va a utilizar.
En la siguiente figura, se puede apreciar la ventana que se despliega al ingresar a
la librería de conductores.
Figura 43. Librería de conductores en ETAP
Fuente: Software ETAP Help 12.6. Elemento PV Array [Imagen]. PV Array /
Inverter / Cable library. [Consultado: 20 de abril de 2017].
3.2 ANÁLISIS DE FLUJO DE POTENCIA
El propósito de una simulación es básicamente entender el comportamiento de un
sistema en particular, y poder evaluar diferentes estrategias para el
funcionamiento del mismo. Es esto, lo que se busca realizar con el sistema
eléctrico de la Universidad Autónoma de Occidente, mediante el módulo Load
Flow de ETAP. El módulo para el análisis de flujo de potencia en el software
calcula las tensiones de las barras, los factores de potencia de las derivaciones,
las corrientes y los flujos de potencia en un sistema eléctrico ya sea radial o anillo.
Al igual que la mayoría de los software de simulación de sistemas eléctricos de
potencia, ETAP permite seleccionar entre varios métodos de simulación, según
sea la preferencia del usuario. La herramienta tiene los siguientes métodos como
opción: Gauss-Seidel, desacoplado rápido, Newton Raphson desacoplado y por
último, el método de Newton Raphson. Este último, será el que se utilizara para
este estudio.
77
4. VERIFICACIÓN DEL MODELO DE ETAP
En este apartado, solo se enfatizara en el modelado de los sistemas fotovoltaicos
que están ubicados en el edificio central y los parqueaderos del ala sur y norte,
que juntos, suman una potencia de 150 kWp. La validación se realizará mediante
una regresión lineal que permite obtener cuál es la correlación qué existe entre
dos variables. En este caso, se busca determinar la correlación entre la potencia
real medida en el año del 2016, registrada por los equipos de medición del
sistema, y la potencia simulada en la herramienta ETAP. Además, después de
establecidos los valores de estas dos variables, se pretende también
correlacionarlas con la variable irradiación, que es medida en el sitio de la planta.
Para determinar la variable Potencia Simulada, se empleó en la herramienta
ETAP, el elemento modulo-inversor. En este se puede ingresar la cantidad de
módulos en serie y paralelo que constituyen el arreglo del sistema fotovoltaico.
Para el caso del sistema de la universidad, cuenta con 2 convertidores DC/AC
(inversores) de 20 kW, y 4 con una potencia de 24kW. Cada uno cuenta con 2
entradas mppt. Por tanto, se establecieron 6 arreglos de módulos. Es decir, 6
elementos modulo-inversor en la herramienta. A continuación se presenta el
modelado de estos elementos en ETAP.
Figura 44. Modelo del sistema fotovoltaico en ETAP versión 12.6
Las simulaciones se realizaron con base a la irradiación registrada el año anterior
(2016) por el piranometro de la universidad. Normalmente, la curva diaria de un
generador fotovoltaico, se establece entre las 7:00 A.m. y a las 6:00 P.m. Pero
para este trabajo se tuvo que establecer una consideración con respecto a este
78
rango, pues la herramienta cuenta con 10 categorías de generación, cada una de
ellas representará una hora del día. En este caso, no se tendrán en cuenta las dos
horas del día, que se ubican en los extremos del rango, es decir, las 7:00 A.m. y
las 6 P.m. Según las mediciones en campo, la energía entregada a la red por el
sistema fotovoltaico en esas dos horas, no superan ni siquiera los 7 kWh, que
representa cerca del 1,5% de la generación diaria. Por tanto, el rango que se
estableció para el análisis, inicia desde las 8:00 A.m. hasta las 5:00 P.m. Además,
las simulaciones realizadas para un día, equivalen a 10, una por hora. Lo que
representaría en el año, 3.650 simulaciones de flujo de potencia. Por lo cual, la
elección de un rango de 10 horas, permite tener mayor practicidad para el
estudio.
Después de simulados todos los días del año, bajo las condiciones anteriormente
descritas, se pudo determinar los promedios mensuales de generación de energía.
Estos se consignan en la siguiente tabla:
Tabla 18. Promedios horarios de generación de energía por mes
Promedios de generación diaria [kWh] Promedio
Hora Anual por
Enero Febrero Marzo Abril Mayo Junio Julio Agosto Septiembre Octubre Noviembre Diciembre hora
08:00 10,660 12,909 14,548 13,869 19,076 18,216 15,394 15,807 19,650 21,937 18,555 13,675 16,191
09:00 31,512 35,469 35,469 35,079 41,434 43,382 38,211 39,745 44,549 45,889 46,058 28,336 38,761
10:00 58,649 59,964 64,060 57,952 67,816 68,764 67,909 64,847 69,198 70,196 66,863 45,699 63,493
11:00 79,674 82,901 79,439 73,816 81,789 82,721 80,949 78,733 84,632 86,543 79,040 62,613 79,404
12:00 92,647 88,922 91,912 88,029 86,722 89,742 82,690 78,805 88,763 89,465 79,248 64,600 85,129
Medición
13:00 91,218 89,905 96,336 86,493 84,191 84,527 91,000 75,003 87,300 86,197 78,082 63,140 84,449
14:00 80,717 83,971 87,490 81,977 80,137 74,146 80,832 71,416 79,119 77,825 67,914 54,997 76,712
15:00 67,195 65,158 73,667 59,923 60,597 61,031 62,050 63,993 64,437 59,245 50,553 40,499 60,696
16:00 47,284 49,180 50,240 42,694 36,784 43,948 44,310 43,153 41,617 35,280 35,108 24,098 41,141
17:00 23,662 27,555 25,512 20,626 18,127 22,163 26,985 22,063 18,823 14,203 14,557 11,672 20,496
08:00 0,000 0,000 0,000 0,001 0,001 0,001 0,001 0,001 0,638 1,144 0,001 0,002 0,149
09:00 15,952 23,061 16,013 24,826 33,708 31,627 22,389 36,652 37,095 38,408 41,954 17,841 28,294
10:00 58,872 60,945 57,293 56,897 65,772 69,833 69,926 72,815 79,733 75,258 75,900 51,933 66,265
11:00 89,774 89,311 91,584 78,690 89,936 90,167 89,963 101,556 101,000 97,742 92,667 85,261 91,471
12:00 108,355 102,552 104,667 96,655 102,152 104,567 92,370 105,852 111,033 104,161 95,600 91,000 101,580
Simulación
13:00 109,581 101,934 113,292 100,345 101,387 98,733 102,074 99,704 109,767 99,923 95,620 88,931 101,774
14:00 103,290 102,414 103,083 93,527 98,226 91,633 98,778 96,889 104,933 100,226 85,606 80,565 96,598
15:00 88,710 81,424 89,708 77,138 77,355 72,400 75,407 88,815 86,100 74,774 62,761 57,522 77,676
16:00 60,389 61,552 62,375 44,559 46,282 55,967 51,720 64,815 51,801 37,989 40,506 29,289 50,604
17:00 20,222 28,292 21,643 14,650 10,393 16,053 23,557 19,364 18,816 7,743 8,855 3,314 16,075
Con los promedios mensuales de generación de energía, se obtuvieron los
promedios anuales, y con estos, se logró realizar un cruce entre la potencia
medida, y la potencia simulada. Esto permitió determinar de manera visual, la
potencia que fue entregada por el sistema fotovoltaico en el año 2016, desde el
punto de vista de simulación, y la medición en campo.
79
A la hora de realizar las simulaciones se tiene que tener presente, que en algunos
de los arreglos que se establecieron en el modelo de ETAP, diferían del número
real de módulos que componen los arreglos en la instalación real. El primer caso,
se tiene para los arreglos del parqueadero. Éste tiene dos inversores, y a cada
uno de ellos, llegan dos circuitos desde dos cajas de combinación. Estas cuentan
con dos arreglos: uno es de cinco cadenas y diez módulos en serie, mientras que
la otra son cuatro cadenas y once módulos en serie. En ETAP, se pudo establecer
un arreglo general para cada inversor de nueve cadenas, de diez módulos en
serie. Por tanto, se agregaron al modelo real, seis módulos fotovoltaicos por cada
arreglo, que representan cerca de 1,5 kW. Para el subsistema del parqueadero
del ala sur, en ETAP se pudo establecer un arreglo de 110 módulos, dispuestos en
diez cadenas y once módulos en serie. Mientras que sitio, se encuentran 115
módulos, 60 (6x10) llegan a una caja de combinación y 55 (5x11) llegan a la otra,
para un total de 115 módulos, que se conectan a uno de los dos inversores de 24
kW ubicados en esta área; cerca de 2,35 kW se despreciaron del modelo del
subsistema del parqueadero de ala sur. Para el último caso que contempla los
arreglos de la azotea del edificio central, también se despreciaron diez módulos en
el modelo, pues el arreglo que se pudo establecer fue de once cadenas y diez
módulos en serie para cada inversor. A continuación se presenta el
comportamiento de un día promedio para las dos variables bajo análisis.
Figura 45. Curva de generación diaria promedio anual
También se considera importante presentar la información de manera mensual, es
decir, la sumatoria de energía producida en cada hora del rango establecido, para
los doce meses del año 2016. Para ello, se muestra las siguientes tablas.
80
Tabla 19. Total de energía generada mensualmente por hora
Generación mensuales - kWh/mes Generación
Hora
Enero Febrero Marzo Abril Mayo Junio Julio Agosto Septiembre Octubre Noviembre Diciembre Anual
08:00 330,5 374,4 347,8 416,1 591,3 546,5 477,2 490,0 589,5 680,1 556,7 423,9 5823,9
09:00 976,9 1028,6 822,3 1052,4 1284,5 1301,5 1184,5 1232,1 1336,5 1422,6 1381,8 878,4 13902,0
10:00 1818,1 1739,0 1512,6 1738,6 2102,3 2062,9 2105,2 2010,3 2075,9 2176,1 2005,9 1416,7 22763,5
11:00 2469,9 2404,1 1917,5 2214,5 2535,5 2481,6 2509,4 2440,7 2539,0 2682,8 2371,2 1941,0 28507,3
12:00 2872,1 2578,7 2222,3 2640,9 2688,4 2692,3 2563,4 2443,0 2662,9 2773,4 2377,4 2002,6 30517,3
Medición
13:00 2827,8 2607,2 2339,0 2594,8 2609,9 2535,8 2821,0 2325,1 2619,0 2672,1 2342,5 1957,3 30251,5
14:00 2502,2 2435,2 2062,1 2459,3 2484,3 2224,4 2505,8 2213,9 2373,6 2412,6 2037,4 1704,9 27415,6
15:00 2083,0 1889,6 1780,4 1797,7 1878,5 1830,9 1923,5 1983,8 1933,1 1836,6 1516,6 1255,5 21709,3
16:00 1465,8 1426,2 1205,1 1280,8 1140,3 1318,4 1373,6 1337,7 1248,5 1093,7 1053,3 747,0 14690,6
17:00 733,5 799,1 606,7 618,8 562,0 664,9 836,5 684,0 564,7 440,3 436,7 350,1 7297,3
08:00 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 19,1 35,5 0,0 0,0 54,8
09:00 494,5 668,8 384,3 720,0 1044,9 948,8 638,5 989,6 1112,8 1190,6 1258,6 410,3 9861,9
10:00 1825,0 1767,4 1375,0 1650,0 2038,9 2095,0 1960,0 1966,0 2392,0 2333,0 2277,0 1194,5 22873,9
11:00 2783,0 2590,0 2198,0 2282,0 2788,0 2705,0 2498,0 2742,0 3030,0 3030,0 2780,0 1961,0 31387,0
12:00 3359,0 2974,0 2512,0 2803,0 3166,7 3137,0 2614,0 2858,0 3331,0 3229,0 2868,0 2093,0 34944,7
Simulación
13:00 3397,0 2956,1 2719,0 2910,0 3143,0 2962,0 2872,0 2692,0 3293,0 3097,6 2868,6 2045,4 34955,7
14:00 3202,0 2970,0 2474,0 2712,3 3045,0 2749,0 2772,0 2616,0 3148,0 3107,0 2568,2 1853,0 33216,5
15:00 2750,0 2361,3 2153,0 2237,0 2398,0 2172,0 2135,0 2398,0 2583,0 2318,0 1882,8 1323,0 26711,2
16:00 1872,1 1785,0 1497,0 1292,2 1434,7 1679,0 1431,5 1750,0 1554,0 1177,7 1215,2 673,6 17362,0
17:00 626,9 820,5 519,4 424,9 322,2 481,6 671,0 522,8 564,5 240,0 265,7 76,2 5535,7
Con estas dos tablas en donde se consignan los datos obtenidos mediante los
medidores del sistema fotovoltaico de 150 kWp de la institución, y los datos
obtenidos mediante las simulaciones en la herramienta ETAP, se pudo realizar la
siguiente grafica que se muestra a continuación.
Figura 46. Generación de energía mensual (2016)
Se podría afirmar, que los valores de generación obtenidos mediante las
simulaciones, tienen una alta simetría con los valores medidos, y prueba de esto,
81
se puede observar en la tabla 19, donde se presenta el error porcentual que se
obtuvo para cada mes del año 2016. En ella se observa que el mes de septiembre
es donde se obtiene el mayor error porcentual, siendo este de 17,19 %, seguido
de enero con un error igual a 12,33%; estos meses fueron en donde mayor
radiación mensual promedio se tuvo en el año, de 5,56 kWh/m^2 y 5,2 kWh/m^2,
respectivamente.
Tabla 20. Calculo de errores porcentual por mes (2016)
Enero Febrero Marzo Abril Mayo Junio Julio Agosto Septiembre Octubre Noviembre Diciembre Total
Medición 18,08 17,28 14,82 16,81 17,88 17,66 18,30 17,16 17,94 18,19 16,08 12,68 202,88
Simulación 20,31 18,89 15,83 17,03 19,38 18,93 17,59 18,53 21,03 19,76 17,98 11,63 216,90
Error porcentual 12,33% 9,32% 6,86% 1,29% 8,42% 7,19% 3,87% 8,01% 17,19% 8,62% 11,85% 8,26% 6,91%
En este punto, los datos obtenidos mediante las simulaciones en la herramienta
debían ser evaluados, y para ello, el autor estableció un método con el cual podría
ser validado dicho modelo. Este se muestra a continuación en la figura 48.
Figura 47. Flujograma para verificar el modelo de ETAP
82
4.1 INTERPRETACIÓN DE RESULTADOS UTILIZANDO TECNICAS
ESTADISTICAS
4.1.1 CORRELACIÓN
En primer lugar, se pretende utilizar un diagrama de dispersión para estudiar la
relación existente entre dos variables cuantitativas; en este punto, se tienen las
variables Potencia medida Vs Irradiancia, Potencia simulada Vs Irradiancia, y
Potencia medida Vs Potencia simulada, bajo interés, siendo la última
correlación, la de mayor interés. Un análisis de correlación intenta medir la fuerza
con la cual dos variables se relacionan entre sí 33. El análisis permite determinar la
magnitud y la dirección de dicha relación; la primera indica la fuerza de la relación
que existe entre las variables y la segunda, se refiere a si ésta es positiva o
negativa. Estas dos características se expresan mediante el coeficiente de
correlación (rxy ), que varía entre -1 y 1, siendo (rxy ) = −1 una relación lineal
perfecta negativa, (rxy ) = 0 una relación inexistente entre las variables, y (rxy ) = 1
una relación lineal perfecta positiva. Cuando las relaciones no son ni -1, ni 1, pero
varían dentro de ese rango, se les denomina relaciones imperfectas, y se
clasifican de la siguiente manera 34:
Si (𝑟𝑟𝑥𝑥𝑥𝑥 ) = 0 ---------------- Nula
Si 0 < (𝑟𝑟𝑥𝑥𝑥𝑥 ) < 0,2 ------- Muy Baja
Si 0,2 ≤ (𝑟𝑟𝑥𝑥𝑥𝑥 ) < 0,4 ----- Baja
Si 0,4 ≤ (rxy ) ≤ 0,6 ----- Moderada
Si 0,6 ≤ (𝑟𝑟𝑥𝑥𝑥𝑥 ) < 0,8 ----- Alta
Si 0,8 ≤ (𝑟𝑟𝑥𝑥𝑥𝑥 ) < 1 -------- Muy alta
Si (𝑟𝑟𝑥𝑥𝑥𝑥 ) = ±1 -------------- Perfecta
Los resultados que se muestran a continuación, fueron determinados a partir de la
herramienta para estadística, Minitab.
[Link] Potencia medida Vs irradiancia. En este primer análisis, se tenían datos
reales medidos en campo. El primero corresponde a la variable independiente,
33
WALPOLE, MYERS. Probabilidad y estadística para ingeniería y ciencias. [en línea]. novena
edición. México: Pearson. p. 431. [Consultado 28 de julio de 2017]. Disponible en Internet:
[Link]
[Link]
34
SIERRA, Enrique. Pruebas de correlaciones de Pearson y Spearman. [Diapositivas] Scribd.
[Consultado: 29 de julio de 2017]. Disponible en internet:
[Link]
enrique-sierra
83
que esta consignada en el eje de las abscisas como irradiancia. En el eje de las
ordenadas, se encuentra la variable potencia, que hace referencia a la potencia
entregada por el sistema fotovoltaico de 150 kWp. En la siguiente grafica se
muestra el diagrama de dispersión obtenido mediante el análisis de correlación.
Figura 48. Correlación Potencia medida Vs Irradiancia
Diagrama de dispersión
Potencia Medida [kW] = 6,269 + 103,3 Irradiancia [kW/m^2]
140
S 8,28892
120
R-cuad. 92,5%
R-cuad.(ajustado) 92,5%
Potencia Medida [kW]
100
80
60
40
20
0,0 0,2 0,4 0,6 0,8 1,0 1,2
Irradiancia [kW/m^2]
Para esta primera correlación, hay alrededor de seis puntos atípicos en el
diagrama, es decir, son puntos que no caen o no tiene cercanía sobre la recta.
Pero como se puede observar en el diagrama de dispersión, el grado de relación
entre estas dos variables es alto, pues el coeficiente de correlación es de 92,5%,
además de ser una relación positiva. Es importante mencionar, que la función que
aparece encima del grafico de dispersión, hace referencia al modelo matemático
que mejor representa los puntos que rodean la recta del ajuste lineal.
[Link] Potencia simulada Vs irradiancia. Ahora, se procede a realizar la
relación entre una variable experimental y una variable medida en campo. La
primera como su nombre lo indica, es la potencia simulada, y la segunda hace
referencia a la irradiancia, siendo la misma variable que se analizó en la anterior
correlación. El diagrama de dispersión para este segundo caso, se presenta a
continuación:
84
Figura 49. Correlación Potencia simulada Vs Irradiancia
Diagrama de dispersión
Potencia Simulada [kW] = - 10,45 + 148,4 Irradiancia [kW/m^2]
175
150 S 8,42448
R-cuad. 96,1%
Potencia Simulada [kW]
125 R-cuad.(ajustado) 96,1%
100
75
50
25
0,0 0,2 0,4 0,6 0,8 1,0 1,2
Irradiancia [kW/m^2]
En esta segunda correlación, se refleja un comportamiento esperado. Pues desde
que se realizaron las simulaciones de flujo de potencia en la herramienta ETAP, se
pudo comprobar, que los arreglos de módulos fotovoltaicos no entregaban
potencia cuando el valor de irradiancia era en promedio, inferior a los 0,250
kW/m^2. Esto puede deberse al algoritmo interno que tiene el elemento en ETAP,
que impide el aporte de potencia a la red. Normalmente estos casos se presentan
a eso de las ocho de la mañana, y las pérdidas del inversor son igual a la potencia
generada, por lo cual la entrega de potencia es nula. Esto se ve reflejado de igual
manera en la correlación, pues en la gráfica anterior, se observa como para los
valores de irradiancia menores a 0,25 kW/m^2 aproximadamente, la potencia
simulada es igual a cero. En el momento que se supera esta franja, se observa
como los puntos empiezan a alinearse con respecto a la línea recta, y como era de
esperarse, la dispersión de los datos se reduce rotundamente, debido al modelo
utilizado por la herramienta ETAP en las simulaciones. El coeficiente de
correlación para este caso es de 96,1%, es decir, se mejora en 3,6 puntos.
[Link] Potencia medida Vs Potencia simulada. Para este último caso de
correlación, se busca determinar la relación entre las dos variables de mayor
interés en este trabajo. Con estas se pretende encontrar un modelo matemático,
que permita conocer el comportamiento de la potencia medida, con los valores de
potencia simulados en ETAP. En el siguiente gráfico, se muestra el
comportamiento de las variables bajo análisis.
85
Figura 50. Correlación Potencia medida Vs Potencia simulada
Diagrama de dispersión
Potencia Medida [kW] = 14,88 + 0,6751 Potencia Simulada [kW]
120
S 9,35189
100 R-cuad. 90,5%
R-cuad.(ajustado) 90,5%
Potencia Medida [kW]
80
60
40
20
0 20 40 60 80 100 120 140
Potencia Simulada [kW]
En esta oportunidad, se replica la misma situación que se presentó en el caso
anterior. Cuando se tenía un valor de irradiancia menor a 0,25 kW/m^2, la entrega
de potencia del sistema fotovoltaico, fue cero. Mientras que para valores de
irradiancia mayor a ese límite, la potencia aportada a la red, se encontraba cerca
de los 50 kW. Justamente esto, es una consecuencia para que en el diagrama de
dispersión, los puntos de la correlación tengan una fragmentación que ocurre
cuando la potencia simulada es aproximadamente de 20 kW. En el límite inferior
se encuentran en su mayoría, los puntos simulados en la primera hora del rango
seleccionado, es decir, las ocho de la mañana. Normalmente, a esa hora del día la
energía aportada por el sol, es insuficiente para que el generador fotovoltaico en la
herramienta ETAP, pueda aportar potencia. Al superar ese intervalo de la primera
hora del rango, los valores de Potencia simulada, comienzan a tener cierta
proporcionalidad con respecto a la variable Potencia medida, ya sea que se
encuentren esos valores por encima o por debajo de la referencia, que es la
ecuación de la recta determinada por la regresión. Por ello, a partir del límite
superior de los 25 kW, la relación entre las variables comienza a ser mucho más
fuerte, tanto así, que la correlación toma un valor de 90,5%.
4.1.2 REGRESIÓN LINEAL SIMPLE
Con un análisis de regresión lineal simple se consiguen dos cosas: con la primera
se busca describir la relación lineal que existe entre dos variables cuantitativas
mediante la ecuación de la recta que mejor se ajuste a los datos de interés; la
segunda consiste en utilizar esa ecuación que se determinó, para realizar una
predicción de los valores de una variable (dependiente), usando la información
aportada por la otra variable (independiente). Esta ecuación mencionada, se le
86
conoce como ecuación de regresión lineal simple, y está representada de la
siguiente manera:
𝑌𝑌𝑖𝑖′ = 𝐵𝐵 ∙ 𝑋𝑋𝑖𝑖 + 𝐴𝐴
Dónde: 𝑌𝑌𝑖𝑖′ : valor estimado de Y para el valor 𝑋𝑋𝑖𝑖 .
A: constante de regresión lineal.
B: Coeficiente de regresión lineal.
La pendiente de la recta, está determinada por el coeficiente de regresión lineal, y
para ello, hay varias maneras de interpretar dicho valor. Por tanto:
• Si B > 0 la relación lineal es positiva y el valor absoluto de B representa el
número de unidades que tiende a aumentar la variable Y, por cada unidad que
aumenta la variable X.
• Si B < 0, la relación lineal es negativa y el valor absoluto de B representa el
número de unidades que tiende a disminuir la variable Y por cada unidad que
aumenta la variable X.
• Si B=0, la ecuación de regresión lineal no es el modelo más adecuado para
describir la relación entre las variables involucradas.
La constante de regresión lineal indica el valor correspondiente a la variable
dependiente cuando la variable independiente asume un valor igual a cero. Es
decir, es el valor con el que se corta el eje de las ordenadas.
[Link] Ajuste lineal: Potencia medida Vs Potencia simulada. Para este
apartado, se tendrá en cuenta solamente el último análisis de correlación que se
realizó en la sección anterior, aquel que involucra las variables potencia medida y
potencia simulada. La intención de realizar un ajuste lineal a los datos, es poder
realizar un mejor acercamiento entre los valores medidos en campo, y los valores
simulados en ETAP. Y de esta manera poder predecir con un menor margen de
error, la potencia entregada por el sistema fotovoltaico en la herramienta
computacional ETAP.
87
Los datos obtenidos de las simulaciones fueron ingresados al programa MiniTab,
con el cual se realizaron las correlaciones, que permitió determinar el modelo de
regresión que más se ajustaba a cada caso. Para este, la ecuación con la que se
lograba un mejor ajuste, era la siguiente:
Pot. Medida = 14,88 + 0,6721 * Pot. Simulada
Recalculando los valores de potencia simulada con el ajuste del modelo, se
espera obtener un mayor acercamiento a los valores obtenidos en campo. Esto se
corroborará con la siguiente grafica que se muestra a continuación.
Figura 51. Correlación Potencia medida Vs Potencia simulada ajustada
Diagrama de dispersión
Potencia Medida [kW] = 0,0019 + 1,000 Potencia simulada ajustada [kW]
120
S 9,35189
100 R-cuad. 90,5%
R-cuad.(ajustado) 90,5%
Potencia Medida [kW]
80
60
40
20
20 40 60 80 100 120
Potencia simulada ajustada [kW]
El coeficiente de correlación para este caso ajustado, es el mismo que para la
correlación que involucra las variables Potencia medida y Potencia simulada,
que fue de 90,5%, pues lo único que se logró con el ajuste, fue acotar los datos
para que se tuviera una mejor proporcionalidad con respecto a la variable
Potencia medida. Por ello, la modificación que se puede evidenciar en la figura
anterior, es el punto de partida del eje de las abscisas. Éste paso de ser 0 kW a 15
kW. En la siguiente sección, se verifica mediante el error porcentual absoluto
medio (MAPE), la mejoría que se tuvo realizando el ajuste lineal.
4.1.3 MAPE
El Error Porcentual Absoluto Medio (MAPE o Mean Absolute Percentage Error) es
un indicador de desempeño que mide el tamaño del error (absoluto) en términos
porcentuales. El hecho que se estime una magnitud del error porcentual lo hace
88
un indicador frecuentemente utilizado por los encargados de elaborar pronósticos
debido a su fácil interpretación 35. La fórmula utilizada para el cálculo del MAPE es
la siguiente:
𝑛𝑛
100% 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑖𝑖 − 𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃ó𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑖𝑖
𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀 = �� �
𝑛𝑛 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑖𝑖
𝑖𝑖=1
Este indicador se utilizará en tres oportunidades para el caso Potencia medida Vs
Potencia simulada: la primera, corresponde a las condiciones reales de
simulación, es decir, sin tener en cuenta el modelo de regresión lineal; en la
segunda, se tendrá en cuenta los nuevos valores de potencia obtenidos mediante
el ajuste lineal (Predicción); y para la tercera se realizará una consideración, que
más adelante se describirá. Por lo pronto, se desarrollaran las dos primeras.
[Link] Caso 1: Valores obtenidos mediante la simulación. Para este caso, se
compararon uno a uno los datos de medición y simulación de potencia para todo el
año del 2016, con excepción de 19 días en donde no se tuvo registros de
irradiancia. Estos días fueron: del 26/marzo al 1/abril, del 29 julio al 3 de agosto,
del 8 al 12 de diciembre, y el 27 de diciembre. Por parte de la variable Potencia
medida, hubo que omitir 7 registros, pues pudo deberse a problemas de
comunicación, que impidieron la obtención de la medida. En total fueron 3.414
registros de la variable Potencia medida que fueron contrastados con la variable
Potencia Simulada. A continuación se muestra un breve resumen de los
resultados obtenidos.
Tabla 21. Resultados caso 1
Sumatoria del error
Muestra MAPE
porcentual absoluto
3414 1.419,88 41,6%
El error absoluto promedio para los registros donde la potencia simulada arrojaba
un valor de cero, fue de alrededor 15 kW. Fueron cerca de 706 registros que
presentaban dicho comportamiento, lo que representa aproximadamente el 22%
de la totalidad de los datos. Esto da respuesta al porqué de un MAPE de 41,6%,
hay muchos registros donde la potencia entregada por el sistema fotovoltaico en
ETAP, es nulo. Por ello, se plantea el siguiente caso.
35
GEO Tutoriales. Error porcentual absoluto medio (MAPE) en un pronóstico de demanda [en
línea]. Blogs Gestión de Operaciones. 26 de Enero de 2015, párr. 1. [Consultado: 30 de julio de
2017]. Disponible en Internet: [Link]
porcentual-absoluto-medio-mape-en-un-pronostico-de-demanda/
89
[Link] Caso 2: Valores obtenidos mediante el ajuste de la simulación. Como
se mencionó anteriormente, este caso toma mayor importancia a partir del error
porcentual absoluto medio que se tuvo en el primer caso. Pero gracias a la
regresión lineal realizada, se pudo tener un ajuste para esos valores de cero, pues
con el modelo matemático, se sobrescribe un valor de aproximadamente 15 kW. A
continuación se presenta el resumen de los resultados.
Tabla 22. Resultados caso 2
Sumatoria del error
Muestra MAPE
porcentual absoluto
3414 976,89 28,62%
Con la modificación de los registros, el error se redujo en unos 12,98 puntos. Pero
a pesar de disminuir el MAPE, el problema que se tiene con la entrega de
potencia nula, solo se corrigió en la mayoría de los registros, pero existen otros en
los cuales el error absoluto porcentual, supera el 2.500%. Esto ocurre debido a
que el modelo matemático determinado, no discrimina entre rangos de irradiancia
que permiten según sea el caso, una mayor o menor entrega de potencia. Es
decir, sí para el MAPE del caso uno, el error porcentual absoluto era de 100%, es
porque el dato de Potencia simulada correspondía a un valor de cero, y para el
dato de Potencia medida se supone a modo explicativo, un valor de 1 kW. Por
tanto, cuando se realiza el ajuste, el valor que tomaría la variable Potencia
simulada ajustada sería de 14,88 kW. Por tanto el error porcentual absoluto para
este dato en específico corresponde a 1.388%, es decir, 13,88 veces la Potencia
medida.
Para concluir este capítulo, es importante mencionar que el modelo se ajusta con
lo que se desea en esta investigación, siempre y cuando se tenga una condición
mínima de irradiancia de 250W/m^2. En la sección [Link] se presentó como los
valores de potencia se hacían cero cuando la irradiancia no superaba dicho
umbral. Por tanto, cuando se supera ese umbral, el modelo comienza a tener un
menor error. Esto se puede evidenciar, filtrando los datos según la irradiancia;
para este caso, se excluirán 511 registros, los contienen una irradiancia menor a
100W/m^2, y error porcentual se redujo en aproximadamente en trece puntos
porcentuales.
Tabla 23. Resultado filtrando registros según irradiancia
Sumatoria del error
Muestra MAPE
porcentual absoluto
2903 453,92 15,64 %
90
5. INTEGRACIÓN DEL SISTEMA FOTOVOLTAICO CONECTADO A RED
Este capítulo tiene como finalidad, realizar el modelamiento y la simulación del
sistema eléctrico de la universidad autónoma de occidente en el software de
sistemas de potencia ETAP, que permitirá evaluar ciertas condiciones típicas de
operación, mediante un análisis de flujo de potencia.
5.1 IDENTIFICACIÓN DE LA TOPOLOGÍA DE RED
Como ya se mencionó en la sección 2.3, la red eléctrica de la institución es de tipo
radial, y la demanda instalada, se aproxima a los 1000 kVA. Está compuesta por
tres nodos principales con un voltaje nominal de 480 V, los cuales se conectan
mediante tres transformadores a la red de Emcali de 13,2 kV. El primero de ellos,
tiene una potencia nominal de 800 kVA, el segundo es de 1000 kVA Y la potencia
del tercero es de 500 kVA. Los tres con una relación de tensión de 13,2/0,48 kV.
Para el caso de las cargas, se modelaron tres cargas, que representarán la
demanda total en cada una de las subestaciones. Adicional a ello, los generadores
fotovoltaicos de 150 kWp y 250 kWp caracterizados en el capítulo 2, se conectan a
la subestación 1 y 2, respectivamente. Cada subestación cuenta con bancos de
compensación capacitiva, de 100 kVar, 120 kVar y 300 kVar para las
subestaciones 1, 2 y 3, respectivamente, pero como se va a presentar más
adelante, estos bancos no tendrán importancia en el análisis, puesto que el aporte
de energía que realizan los sistemas fotovoltaicos durante el día, permite tener un
óptimo nivel de tensión en cada uno de los nodos.
Para los parámetros de los transformadores, se optó por seleccionar unos valores
típicos, que por defecto el software calcula. Según la ventana de ayudas de ETAP,
el cálculo de la impedancia típica y la relación X/R provienen de dos fuentes: Ansi
C57.12.10 y del Power System Handbook by Beeman. Por tanto, a continuación
se muestran los parámetros de los transformadores que alimentan las cargas de
cada una de las subestaciones. Los valores para la secuencia positiva y cero, son
iguales.
Tabla 24. Parámetros de transformadores de la red
Elemento Capacidad Voltaje R (%) X (%) Zcc (%) X/R
Transformador 1 800 kVA 13,2/0,48 kV 1,77 5,471 5,75 3,09
Transformador 2 1000 kVA 13,2/0,48 kV 0,979 5,666 5,75 5,79
Transformador 3 500 kVA 13,2/0,48 kV 0,979 5,666 5,75 5,79
Finalmente, todas las cargas de la red fueron clasificadas como trifásicas
91
balanceadas de tipo potencia constante, cuya capacidad oscilada entre los 244
kVA para la carga uno, 285 kVA para la carga dos y 205 kVA para la carga tres.
En la siguiente sección, se explicará de manera más detallada, el comportamiento
de estas bajo las simulaciones.
5.2 MODELAMIENTO DE LA RED
La finalidad que se tiene con este apartado, es mostrar cómo fue la elaboración
del modelo de la red, y cuáles fueron las condiciones de operación que se
establecieron para la evaluación de la red con la integración del sistema
fotovoltaico.
En primer lugar, para el modelado del sistema fotovoltaico de la fase I y la fase II,
se establecieron una serie de condiciones con las cuales se realizaran las
simulaciones. Una de ellas es establecer un rango de operación del sistema
fotovoltaico, que estaría limitado desde las 8 de la mañana, a 5 de la tarde. Estos
diez intervalos permitirán conocer el comportamiento de un día bajo unas
condiciones operativas; la segunda corresponde a la irradiancia que se va a
utilizar para el rango anteriormente mencionado. Para ello, se extrajo del portal del
inversor Symo, la irradiancia por hora de todo el año del 2016, y de esta
información se obtuvo los promedios horarios de cada una de las horas que
compone dicho rango; por último, se debía establecer para este rango, los valores
de la temperatura ambiente en sitio. Para esto, se recurrió a un pronóstico de
temperatura de un día al azar en la ciudad de Cali (12 de septiembre de 2017),
que permitirá realizar las simulaciones. En la tabla que se muestra a continuación,
se presentan los valores de temperatura que se tomaron para el rango de 8:00
a.m. a 5 p.m.
Tabla 25. Temperatura horaria en la ciudad de Cali
Temperatura
Hora
(°C)
08:00 26,7
09:00 28,5
10:00 30,0
11:00 30,7
12:00 31,0
13:00 30,7
14:00 30,1
15:00 29,3
16:00 28,4
17:00 27,2
92
Vale la pena aclarar, que se debe redondear los datos para ser ingresados a la
herramienta, pues éste no acepta números reales. Pero como era de esperarse, la
temperatura en Cali, siempre se mantiene dentro de ese promedio. Para el caso
de la irradiancia, a continuación se muestra una tabla en donde se puede observar
su comportamiento durante el año anterior en la institución universitaria. Los
valores se muestran en kW/m^2.
Tabla 26. Irradiancia horaria promedio por mes
Mes
Enero Febrero Marzo Abril Mayo Junio Julio Agosto Septiembre Octubre Noviembre Diciembre
Hora
07:00 0,0024 0,0025 0,0033 0,0073 0,0134 0,0122 0,0096 0,0095 0,0168 0,0170 0,0187 0,0101
08:00 0,0561 0,0677 0,0796 0,0834 0,1101 0,1031 0,0857 0,0988 0,1266 0,1282 0,1070 0,0997
09:00 0,2089 0,2379 0,2186 0,2385 0,2812 0,2877 0,2466 0,3000 0,3161 0,3152 0,3338 0,2147
10:00 0,4495 0,4511 0,4369 0,4292 0,4858 0,5001 0,5010 0,5317 0,5684 0,5378 0,5539 0,4017
11:00 0,6652 0,6568 0,6782 0,5975 0,6567 0,6488 0,6398 0,7435 0,7538 0,7289 0,6902 0,6189
12:00 0,8223 0,7829 0,8044 0,7367 0,7492 0,7756 0,6841 0,7976 0,8534 0,7950 0,7213 0,6970
13:00 0,8305 0,7993 0,8860 0,7712 0,7578 0,7337 0,7649 0,7544 0,8449 0,7614 0,7424 0,6654
14:00 0,7685 0,7820 0,7812 0,7127 0,7264 0,6672 0,7248 0,7227 0,7920 0,7309 0,6409 0,6101
15:00 0,6374 0,6057 0,6559 0,5698 0,5558 0,5186 0,5442 0,6412 0,6210 0,5398 0,4859 0,4324
16:00 0,4377 0,4552 0,4693 0,3778 0,3611 0,4019 0,3886 0,4638 0,4111 0,3311 0,3334 0,2878
17:00 0,2443 0,2666 0,2439 0,1938 0,1757 0,2031 0,2394 0,2305 0,2083 0,1478 0,1612 0,1345
18:00 0,0680 0,0829 0,0836 0,0614 0,0482 0,0607 0,0856 0,0750 0,0470 0,0298 0,0294 0,0369
19:00 0,0025 0,0048 0,0036 0,0026 0,0013 0,0026 0,0045 0,0031 0,0005 0,0000 0,0000 0,0004
Con estos datos se obtuvo los promedios para cada una de las horas del día bajo
estudio. Por tanto, a continuación se presentan los valores de irradiancia utilizados
para estimar la producción de energía durante cada hora del día.
Tabla 27. Irradiancia promedio
Irradiancia
Hora
[W/m^2]
08:00 95
09:00 267
10:00 487
11:00 673
12:00 768
13:00 776
14:00 722
15:00 567
16:00 393
17:00 204
Ya establecidos unos parámetros de entrada con los cuales se puede determinar
la entrega de potencia hora tras hora, se debe conocer también, los
comportamientos de demanda que tiene la universidad. Para ello, en el Capítulo
93
dos se realiza una caracterización de los consumos, con estos, el autor propone
cuatro condiciones que permitirán evaluar modos de operación de la red. Los
comportamientos están involucrados directamente con el periodo del año, por tal
razón, las condiciones de operación estarán dadas por los siguientes nombres:
• Periodo académico – días entre semana (ES): Para este primer caso, se tendrá
en cuenta el promedio de consumo horario para los días hábiles de los primeros 5
meses del año 2017, desde enero hasta mayo. A continuación se presenta el
comportamiento diario promedio de la carga durante un periodo académico.
Tabla 28. Demanda de potencia por hora Periodo Académico ES
PERIODO ACADEMICO - DIAS ENTRE SEMANA
SUB 1 SUB 2 SUB 3
Hora
kW kVar kW kVar kW kVar
08:00 209,41 62,16 212,87 76,96 168,66 39,90
09:00 245,30 64,54 259,23 97,35 191,11 52,71
10:00 232,91 64,57 266,21 101,64 196,31 56,44
11:00 228,52 64,42 280,66 107,88 198,94 58,99
12:00 218,40 64,81 282,47 106,76 202,78 60,45
13:00 203,22 65,74 250,73 100,93 199,96 57,43
14:00 220,83 64,45 266,09 102,81 203,81 61,12
15:00 257,13 64,65 297,48 113,17 207,34 62,54
16:00 272,41 63,55 287,12 107,49 204,59 62,07
17:00 274,82 63,04 269,37 97,20 196,56 58,56
Promedio 236,30 64,19 267,22 101,22 197,01 57,02
• Periodo académico – fines de semana (FDS): En este segundo caso, se tienen
los mismos meses que para el caso anterior, con la única diferencia de que se
tendrán, solo los consumos de los días festivos de un periodo académico. Es
decir, de enero a mayo. La siguiente tabla muestra el comportamiento diario
promedio para este caso.
Tabla 29. Demanda de potencia por hora Periodo Académico FDS
PERIODO ACADEMICO - FINES DE SEMANA
SUB 1 SUB 2 SUB 3
Hora
kW kVar kW kVar kW kVar
08:00 91,55 29,55 122,00 39,67 78,00 -5,06
09:00 84,05 30,34 138,65 49,57 67,50 -10,66
10:00 73,43 30,16 148,36 51,48 66,93 -10,97
11:00 63,56 31,62 151,28 54,03 68,41 -10,14
12:00 51,27 31,21 150,26 55,36 63,69 -14,27
13:00 50,23 30,46 147,65 54,85 63,90 -13,26
14:00 47,19 28,95 144,16 55,54 52,02 -19,04
15:00 62,42 30,34 128,15 44,28 50,94 -19,18
16:00 72,41 28,92 116,22 39,56 44,25 -21,85
17:00 83,77 28,89 101,62 31,30 34,78 -26,47
Promedio 67,99 30,04 134,84 47,57 59,04 -15,09
94
• Periodo vacacional – días entre semana (ES): Para este tercer escenario que
se propone, se considera los consumos que comprenden los meses de junio y
julio, que representan las vacaciones de verano en la universidad. En la tabla 22,
se registran los comportamientos promedio de consumo diario para este caso.
Tabla 30. Demanda de potencia por hora P. Vacaciones ES
PERIODO DE VACACIONES - DIAS ENTRE SEMANA
SUB 1 SUB 2 SUB 3
Hora
kW kVar kW kVar kW kVar
08:00 175,06 56,15 117,78 22,99 156,36 33,80
09:00 210,48 59,64 112,62 41,99 172,14 43,50
10:00 198,51 59,95 95,29 35,28 192,23 56,85
11:00 184,87 59,49 69,31 31,70 190,12 56,60
12:00 187,75 59,63 75,39 32,34 189,23 57,22
13:00 170,54 58,65 50,02 23,64 184,15 52,96
14:00 184,37 58,40 84,16 27,92 190,84 58,44
15:00 214,49 58,42 117,19 34,20 187,54 56,27
16:00 227,61 59,21 141,75 32,06 182,96 53,39
17:00 242,12 56,58 166,89 19,56 163,31 36,47
Promedio 199,58 58,61 103,04 30,17 180,89 50,55
• Periodo vacacional – fines de semana (FDS): Por último, se tienen los fines de
semana de los meses junio y julio. Y estos tuvieron los comportamientos que se
muestran a continuación.
Tabla 31. Demanda de potencia por hora P. Vacaciones FDS
PERIODO DE VACACIONES - FINES DE SEMANA
SUB 1 SUB 2 SUB 3
Hora
kW kVar kW kVar kW kVar
08:00 81,13 28,41 82,54 6,82 68,36 -9,14
09:00 59,00 27,78 43,23 12,64 71,04 -7,62
10:00 49,10 27,68 19,45 11,79 67,31 -11,91
11:00 37,00 29,96 -7,36 10,63 62,02 -14,72
12:00 24,80 29,42 -15,20 9,62 62,99 -13,44
13:00 26,17 27,64 -11,40 9,10 55,98 -14,27
14:00 26,79 25,33 -6,55 9,82 53,90 -15,66
15:00 46,53 26,46 15,47 -2,38 50,09 -16,70
16:00 62,36 26,73 41,45 -3,72 50,12 -17,43
17:00 82,23 28,70 65,38 -10,49 49,99 -15,80
Promedio 49,51 27,81 22,70 5,38 59,18 -13,67
Se consideró que debido a la poca variación que se tiene con los consumos de la
universidad durante los días de un periodo académico y los días de un periodo
vacacional, los cuatro escenarios que se presentaron anteriormente, permitirán
analizar la interacción que tiene el sistema fotovoltaico con la red eléctrica de la
universidad, teniendo en cuenta ciertas condiciones de operación. En la siguiente
sección, se pretende mostrar los resultados obtenidos.
95
5.3 RESULTADOS Y ANALISIS DE LA SIMULACIÓN EN ETAP
En esta sección tiene como finalidad presentar los resultados obtenidos de las
simulaciones para el análisis de flujo de potencia bajo los cuatro escenarios diarios
de operación que se mencionaron anteriormente. La creación del modelo con el
cual se trabajara fue realizada en la versión 16 de ETAP. Éste, se presenta a
continuación:
Figura 52. Modelo red UAO – sistema fotovoltaico
Como se puede ver en la imagen, los arreglos de módulos ubicados en el lado
izquierdo del modelo, corresponde al sistema fotovoltaico de 150 kWp instalado en
el año del 2015. Mientras que la fase II del sistema, está modelado en la parte
derecha de la figura. Estos sistemas no fueron modelados de la misma manera
que se utilizó para la validación del modelo en el Capítulo cuatro (elemento
módulo-inversor). Para este caso, los módulos se modelaron por separado de los
inversores, y de esta manera se lograba tener el número de módulos exactos en
cada uno de los arreglos de los subsistemas. Esto implica algo, y es que para
llegar al flujo de potencia AC, primero se debe simular un flujo de potencia DC,
obteniendo así, la potencia activa entregada por cada arreglo. Esto se realizó en
diez ocasiones, debidas a las diez horas que se tiene en el rango seleccionado.
En este punto es importante mencionar, que el software considera al inversor
como una fuente de tensión. Razón por la cual, fue posible omitir seis inversores
en los subsistemas de la fase II, uno por cada dos arreglos de módulos. En el
capítulo dos, se mostró que cada edificio está compuesto por dos inversores, pero
en este caso y debido a las condiciones de simulación, solo fue necesario utilizar
un elemento para cada subsistema, con la única condición de que ese inversor
96
representara la potencia nominal de dos. Adicional a ello, la potencia entregada
diariamente por el generador fotovoltaico no variará en ninguno de los cuatro
escenarios planteados, pues desde un principio se estableció una irradiancia de
un día típico, basado en los promedios horarios del año 2016. Caso contrario
ocurrirá con las cargas que si tendrán variación para cada uno de los ejercicios,
pues estas serán las encargadas de recrear las distintas condiciones de operación
que tiene la institución durante un día cotidiano. A continuación, se muestran los
resultados obtenidos de para cada uno de los escenarios establecidos.
5.3.1 Periodo académico – ES. El resultado obtenido para el día típico que se
determinó para estas condiciones, fue el siguiente:
Tabla 32. Resultados periodo académico ES
Generación Demanda [kVA]
Hora
Red [kVA] SFV 150 kWp SFV 250 kWp S/E - 1 S/E - 2 S/E - 3
08:00 545,7+195,3i 18,08 29,52 209,4+56,9i 213+80,7i 168,6+48,8i
09:00 578,2+230,5i 45,80 74,85 245,3+66,6i 259,2+98,2i 191,1+55,3i
10:00 492,8+228,7i 78,08 128,28 233+63,3i 266,1+100,8i 196,2+56,8i
11:00 438,4+232,6i 103,91 170,60 228,5+62,1i 280,6+106,3i 199+57,6i
12:00 401,1+230,9i 116,22 191,68 218,3+59,3i 282,4+107i 202,7+58,7i
13:00 348,3+213i 117,22 193,62 203,2+55,2i 250,6+94,9i 200+57,9i
14:00 404,6+225,9i 109,96 181,54 220,9+60i 266,1+100,8i 203,9+59i
15:00 530,8+251,6i 89,24 146,40 257,1+69,8i 297,4+112,6i 207,2+60i
16:00 597,3+253,3i 64,70 106,08 272,5+74i 287+108,7i 204,5+59,2i
17:00 649,6+246,8i 36,06 58,92 274,8+74,7i 269,3+102i 196,6+56,9i
Total diario 4986,8+2308,6i 2060,77 7004,5+2224,1i
De la información anterior, se puede deducir que la energía requerida de la red se
disminuye cuando hay una mayor irradiancia, es decir en horas del mediodía. Para
la 1:00 P.m. la red de distribución entrego 348 kW, siendo ésta, la hora del rango
donde el sistema eléctrico de la universidad requirió menos energía por parte de la
red. Según este caso, para esa misma hora, el sistema fotovoltaico inyecto la
mayor potencia al sistema eléctrico, siendo capaz de cubrir el 47,5% del consumo
total de la institución, mientras que para las 8:00 A.m. el sistema fotovoltaico,
entrego la menor de potencia, tan solo suministro el 8% de la demanda total para
esa hora. Al terminar el rango, de 8:00A.m. a 5:00 P.m., la potencia suministrada
por el sistema fotovoltaico, cubre el 29,41% de la demanda total de dicho rango.
Claro está, que ese porcentaje se verá mucho más reducido, si se compara con la
demanda total de un día, es decir, incluyendo las horas que están por fuera del
rango de generación fotovoltaica. En la figura que se mostrara a continuación, se
podrá ver el comportamiento de la potencia activa generada y la demanda activa
dentro del rango bajo estudio.
97
Figura 53. Generación – Demanda Periodo Académico ES
En la figura anterior, se puede apreciar más claramente el comportamiento de las
variables importantes que componen el sistema. Se puede observar el efecto que
tiene la integración del sistema fotovoltaico con la red eléctrica de la Universidad
Autónoma de Occidente, el cual es capaz de reducir el suministro de energía de la
red en horas del mediodía. En la siguiente gráfica, se muestra los perfiles de
tensión que se consiguieron bajo estas condiciones operativas. La tensión en la
subestación uno, fue la que más vario durante el rango de 98,96% (12:00 hrs) a
98,28% (17:00 hrs), y su aumento se produjo cuando el sistema fotovoltaico
aportaba la mayor cantidad de energía a la red que fue de 117,22 kW.
Figura 54. Perfiles de tensión – Periodo Académico ES
98
Ahora, es el turno de analizar el comportamiento de los días festivos para el
mismo periodo académico.
5.3.2 Periodo académico – FDS. El resultado obtenido para el día típico que se
determinó para estas condiciones, fue el siguiente:
Tabla 33. Resultados periodo académico FDS
Generación Demanda [kVA]
Hora
Red [kVA] SFV 150 kWp SFV 250 kWp S/E - 1 S/E - 2 S/E - 3
08:00 244,5+65,4i 18,08 29,52 91,6+40,5i 122+43,1i 78-19,9i
09:00 170,4+69,8i 45,80 74,85 84+37,1i 138,6+48,9i 67,5-17,2i
10:00 84,4+68,3i 78,08 128,28 73,4+32,4i 148,3+52,3i 67-17,1i
11:00 12,3+64,7i 103,91 170,60 63,6+28,1i 151,2+53,4i 68,4-17,5i
12:00 -38,4+60,4i 116,22 191,68 51,3+22,7i 150,2+53i 63,7-16,3i
13:00 -44,6+59i 117,22 193,62 50,2+22,2i 147,6+52,1i 63,9-16,3i
14:00 -44,3+59,3i 109,96 181,54 47,2+20,8i 144,1+50,9i 52-13,3i
15:00 8,4+60,2i 89,24 146,40 62,4+27,6i 128,1+45,2i 51-13i
16:00 63,5+62,1i 64,70 106,08 72,4+32i 116,2+41i 44,2-11,3i
17:00 125,8+64,7i 36,06 58,92 83,8+37i 101,6+35,9i 34,8-8,89i
Total diario 582+633,9i 2060,77 2618,3+625,41i
En la figura anterior se puede observar el comportamiento horario para la
condición operativa de este día típico, que corresponde a un día festivo
comprendido entre los meses enero a mayo. Para este escenario, se tiene algo
esperado, pues se tienen algunos valores negativos en la columna Red y se
registraron para las 12:00, 1:00, y 2:00 P.m. Esto indica que durante esas tres
horas, la generación fotovoltaica fue suficiente para atender la demanda de la
institución, e incluso, fue capaz de entregar excedentes a la red. La energía
aportada durante el rango, represento el 79% de la demanda. Las horas cuando la
institución requirió de una mayor entrega de energía por parte de la red, fue para
las dos primeras horas del rango, y la ultimo, de 244, 170 y 126 kW para las 8:00
A.m., 9:00 A.m. y 5:00 P.m., respectivamente. Para este caso, la energía aportada
por el sistema fotovoltaico durante el rango de 10%, representa el 79% de la
demanda obtenida para ese mismo rango. A continuación, se presentan en una
gráfica del comportamiento de la demanda y la generación para esta condición.
99
Figura 55. Generación – Demanda Periodo Académico FDS
Los niveles de demanda para esta condición, se redujo cerca del 60%. Razón por
la cual el sistema fotovoltaico es capaz de proporcionar energía suficiente para
alimentar la carga existente, por algunas horas. Por tanto, en un fin de semana
bajo esta condición, la universidad tiene un ahorro de 4.385,9 kW con respecto a
un día hábil del periodo académico bajo el rango establecido. En este caso, los
perfiles de tensión tuvieron un comportamiento similar al caso anterior, pues la
tensión en la subestación uno fluctúa debido al aporte de potencia que tiene el
sistema fotovoltaico de 150 kWp, pues este en horas del mediodía aporta 23,4
veces la demanda de la carga conectada a la subestación mencionada
anteriormente. Mientras que la tensión en la subestación dos se mantiene
constante en horas del mediodía debido a que en ese momento, la generación
representaba solo 1,31 veces la demanda de esa subestación. Incluso, en horas
de la tarde tiene a incrementar la tensión, pues la demanda empieza a descender.
Figura 56. Perfiles de tensión – P. Académico FDS
100
5.3.3 Periodo de vacaciones – ES. El resultado obtenido para el día típico que
se determinó para estas condiciones, fue el siguiente:
Tabla 34. Resultados periodo de vacaciones ES
Generación Demanda [kVA]
Hora
Red [kVA] SFV 150 kWp SFV 250 kWp S/E - 1 S/E - 2 S/E - 3
08:00 403+135i 18,08 29,52 175+51,4i 117,7+34,5i 156,3+43,7i
09:00 377+148,9i 45,80 74,85 210,6+61,8i 112,6+33i 172,2+48,1i
10:00 283+145i 78,08 128,28 198,6+58,3i 95,3+27,9i 192,3+53,7i
11:00 174,1+132,2i 103,91 170,60 184,8+54,3i 69,3+20,3i 190,1+53,1i
12:00 149,6+134,6i 116,22 191,68 187,8+55,2i 75,4+22,1i 189,2+52,9i
13:00 98,9+120,6i 117,22 193,62 170,5+50,1i 50+14,6i 184,2+51,5i
14:00 172,6+136,5i 109,96 181,54 184,4+54,2i 84,2+24,6i 190,8+53,3i
15:00 287,4+154,8i 89,24 146,40 214,6+63i 117,1+34,3i 187,6+52,4i
16:00 384,4+165,7i 64,70 106,08 227,5+66,8i 141,7+41,5i 182,9+51,1i
17:00 480+173,9i 36,06 58,92 242,1+71,1i 166,9+48,9i 163,4+45,6i
Total diario 2810+1447,2i 2060,77 4835,1+1393,3i
Para los resultados obtenidos de esta simulación, hay que realizar una
observación con respecto a la demanda de energía. Los datos de consumo que se
extrajeron de los medidores de las subestaciones, correspondían al periodo
vacacional del año 2017, periodo en el cual iniciaba operaciones la segunda fase
del sistema fotovoltaico, que está conectado a la subestación. Por esta razón, en
la columna S/E – 2, se puede observar que para las horas del mediodía, cuando
hay más brillo solar, la demanda que registra el medidor se empieza a disminuir
debido a la inyección de potencia del sistema fotovoltaico de 250 kWp. Por tal
motivo, en los últimos dos casos, no se tendrán resultados tan satisfactorios. Pero
pueden representar una guía para próximos estudios que se deseen realizar.
El comportamiento del sistema fotovoltaico es el mismo que los casos anteriores,
lo único que está variando es el comportamiento de las cargas. Para este
escenario, el punto máximo se tuvo a las 5:00 P.m. con una demanda de potencia
igual a 572 kW. La hora en la cual se entregó mayor potencia fue a la 1:00 P.m.
con 310 kW, que representaba el 77% de la demanda que había para ese
momento. Este día típico que se estableció para este caso, terminado las diez
horas, la red tuvo que entregar 2.799 kWh, lo cual representa 2.168 kWh menos
que para el caso uno, donde se tenían los días entre semana de un periodo
académico. De este escenario también se puede mencionar que al finalizar el
rango, la generación fotovoltaica logro cubrir el 43% de la demanda total. La
siguiente gráfica permite presentar el comportamiento de la demanda y la
generación que se tuvo para esta condición.
101
Figura 57. Generación – Demanda Periodo Vacaciones ES
Al igual que los casos anteriores, la tensión en la subestación tres, a la cual se
conectan los chillers, solo tiene una leve desviación, pero se mantiene dentro del
rango del 99,5%. Mientras que las tensiones en las subestaciones uno y dos,
varían durante el mediodía con un margen de aproximadamente 0,5%. Las
tensiones en estos dos barrajes empiezan a tener una disminución a partir de las
14:00 hrs, pues es cuando empieza a incrementar la demanda, llegando incluso a
caer por debajo del 98,5% a las 17:00 hrs.
Figura 58. Perfiles de tensión – P. Vacaciones ES
5.3.4 Periodo de vacaciones – FDS. El resultado obtenido para el día típico que
se determinó para estas condiciones, fue el siguiente:
102
Tabla 35. Resultados periodo de vacaciones FDS
Generación Demanda [kVA]
Hora
Red [kVA] SFV 150 kWp SFV 250 kWp S/E - 1 S/E - 2 S/E - 3
08:00 170,4+47i 18,08 29,52 81,1+45,6i 68,1+16,1i 68,4-15,8i
09:00 53,4+27,6i 45,80 74,85 59+33,1i 43,2+10,2i 71-16,4i
10:00 -68,3+17,8i 78,08 128,28 49,1+27,6i 19,5+4,61i 67,3-15,5i
11:00 -179+7,26i 103,91 170,60 37+20,8i -7,35-1,74i 62-14,3i
12:00 -230,5-0,665i 116,22 191,68 24,8+13,9i -15,2-3,6i 63-14,5i
13:00 -235,1+2,52i 117,22 193,62 26,2+14,7i -11,4-2,7i 56-12,9i
14:00 -213+3,95i 109,96 181,54 26,8+15i -6,56-1,55i 53,9-12,5i
15:00 -120,8+19,6i 89,24 146,40 46,5+26,1i 15,5+3,66i 50,1-11,6i
16:00 -15,4+33,8i 64,70 106,08 62,3+35i 41,5+9,82i 50,1-11,6i
17:00 103,3+50,8i 36,06 58,92 82,2+46,2i 65,4+15,5i 50-11,66i
Total diario -735+209,665i 2060,77 1299,49+191,54i
En este último escenario, se tuvo el comportamiento esperado, pero no en las
mismas proporciones. Es decir, la potencia entregada por parte de las
simulaciones, está entregando una potencia adicional a la suministrada realmente
por la planta, pues como se mencionó anteriormente, los datos de consumo
corresponden al mes donde entro en operación la fase dos del proyecto
fotovoltaico. Según los resultados, el sistema fotovoltaico permite entregar
excedentes desde las 11:00 A.m. hasta las 4:00 P.m. y al finalizar el rango, el
balance indica que la energía inyectada a la red en ese rango de diez horas es de
alrededor 735 kWh. Para este último gráfico, el comportamiento de la generación y
demanda se presentan a continuación.
Figura 59. Generación – Demanda Periodo Vacaciones FDS
103
Ya en este último caso, se puede observar la carencia de aceptación que tienen
estos datos. La demanda se ve afectada con la variación de la generación en la
subestación dos. Además de ello, se puede inferir de los resultados, que la
generación fotovoltaica es para este día típico, 1,5 veces la demanda de energía
tenida para dicho rango. Normalmente, después de entrar en operación la
segunda fase, en un día festivo, la universidad entrega excedentes por el orden de
los 400 kWh/día, mientras que para el caso simulado, la entrega de energía fue de
750 kWh/Rango, es decir, se está teniendo casi el doble de la entrega de energía.
Por tanto, este escenario no resulto ser muy provechoso, por los datos de
consumo que se tuvieron. Sin embargo, en la siguiente grafica se puede apreciar
cómo el comportamiento de las tensiones mantuvo la misma proporción con
respecto a los casos anteriores, con la única diferencia de llegar casi al 100,2%.
Esto debido a que el nodo de la subestación dos, se comportó gran parte del
rango como un nodo generador, mientras que para la subestación uno, la
demanda estuvo por debajo del promedio.
Figura 60. Perfiles de tensión – P. Vacaciones FDS
Para culminar con este trabajo, el autor desarrollará la siguiente sección con el fin
de determinar la conveniencia de utilizar almacenamiento en alguno de los
subsistemas instalados en la institución, para ello, se seleccionó el edificio de
Aulas 1. Esto permitirá determinar la necesidad de realizar un trabajo futuro
entorno a ello.
5.4 AULAS 1: GENERACION VS DEMANDA
Esta sección se desarrollara un comparativo entre la generación y la demanda en
el edificio de aulas uno, lo que permitirá determinar la conveniencia de contar con
almacenamiento en el sitio. El ejercicio se realizara con datos de febrero, marzo,
104
abril y mayo del 2017, para el caso de la demanda, los cuales corresponden a los
meses de un periodo académico en la institución. Mientras que para el caso de la
generación solo se pudo conseguir información horaria de una semana, que
corresponde a la segunda semana del mes de septiembre de 2017, que va del
cuatro al diez del mes en mención. El análisis se realizará con tres posibles
escenarios, en donde se tendrá la demanda promedio horaria de los días hábiles,
la demanda horaria de día un sábado y un día domingo. Pueda que los datos que
se utilizaran para el análisis no sean congruentes en cuanto a las fechas en las
cuales se registraron las variables, pero servirán observar el comportamiento de
un día cotidiano en el sitio. Para comenzar, se muestra el comportamiento de un
día hábil, pero antes de ello, en la tabla siguiente se consigna el comportamiento
de la potencia entregada por el subsistema de aulas 1 en las simulaciones de
ETAP.
Tabla 36. Producción de energía simulada para el SFV del edificio de aulas 1
Hora 8:00 9:00 10:00 11:00 12:00 13:00 14:00 15:00 16:00 17:00
kWh 4,5 11,4 19,59 26,03 29,25 29,54 27,69 22,32 16,19 8,98
Figura 61. Aulas 1: Generación Vs Demanda – Día hábil
Como se puede observar en este primer caso, la generación medida del
subsistema fotovoltaico instalado en la terraza del edificio de aulas 1, en ningún
momento del día supera la demanda de energía en ese edificio. La curva de color
rojo que se aprecia en la gráfica corresponde al comportamiento plasmado en la
tabla 36, gracias a las simulaciones obtenidas de la herramienta ETAP. Dicha
generación simulada llega a tener tan solo 1 kWh de diferencia con respecto a la
demanda a las 13:00 hrs, siendo en ese momento, la producción de energía igual
a 29,54 kWh, mientras que la demanda llega a 30,54 kWh. Al terminar el día, la
generación medida del subsistema fotovoltaico es de 162,24 kWh/día, lo que
105
representa el 29,32% de la energía consumida, pues para este caso, el día tuvo
un consumo de 553,25 kWh. Ahora, se procede con el caso del día sábado y para
ello, se muestra el comportamiento en la figura siguiente.
Figura 62. Aulas 1: Generación Vs Demanda – Día sábado
Si bien el consumo de energía se disminuyó en gran proporción, pasando de
553,25 kWh/día a 311,88 kWh/día, el comportamiento de la generación, no tuvo
mayores cambios. Pues en este perfil de consumo se tiene un pico al mediodía de
aproximadamente 35 kWh, mientras que la generación fue de 18,48 kWh y 29,25
kWh para la medida y la simulada, respectivamente. Se debe tener en cuenta que
la curva de la generación simulada será la misma en los tres casos, pues es la
misma que se utilizó en los casos planteados para un día típico de la sección 5.3.
Al culminar el día la producción de energía fue de 135,96 kWh/día, lo que
representa el 43,6% del consumo para el día sábado escogido para este ejercicio.
Y por último, se tiene el comportamiento de un día festivo.
106
Figura 63. Aulas 1: Generación Vs Demanda – Día domingo
Sin duda alguna, el mejor escenario es el que se tiene en la figura anterior, pues la
generación de energía llega a ser en horas del mediodía cerca de seis veces la
demanda, la cual se mantiene constante durante todo el día. Dicho
comportamiento de generación corresponde al día diez de septiembre del 2017.
Pero lo que permite ver esta gráfica, es que siempre y cuando se encuentre dentro
del rango de utilidad de un sistema fotovoltaico, su generación superará el
consumo en el edificio de aulas uno, y de seguro ocurrirá lo mismo en aulas dos,
tres y cuatro. El consumo para este día festivo bajo estudio, fue de solo 57,55
kWh/día, que comparado con un día hábil de clases, se redujo en casi 500
kWh/día, lo cual permite tener excedentes en el día del orden de los 118,12
kWh/día, puesto que para este caso, se produjeron 175,68 kWh, lo que permite
concluir, que la utilización de fuentes de almacenamiento solo serían conveniente
para los días festivos en donde se entregan excedentes a la red.
107
6. CONCLUSIONES
Los sistemas de distribución pueden beneficiarse en gran medida de la presencia
de fuentes de generación distribuida. Prueba de ello, se consigue con la
investigación realizada acerca de los perfiles de demanda diarios que tiene la
Universidad Autónoma de Occidente. Con estos, se pudo comprobar el alivio que
tuvo la red debido a la entrada en operación del sistema fotovoltaico de la fase II,
llegando a generar en promedio, 980 kWh/día, cuando inicialmente, el sistema de
la fase I, llegaba a generar en promedio unos 600 kWh/día, para un total cercano a
los 1.6 MWh generados diariamente. Tanto así, que en los días festivos, la
institución llegó a generar en el mes de agosto de 2017, excedentes a la red de
aproximadamente 400kWh/día. También se pudo comprobar la uniformidad de los
perfiles diarios de demanda que se tienen durante periodos específicos del año.
Razón por la cual, en este trabajo se establecen cuatro condiciones de operación
que permitían crear unos escenarios de simulación.
Con los resultados obtenidos durante la evaluación del modelo desarrollado en la
herramienta computacional ETAP, es claro que existe una fuerte correlación entre
las variables de irradiancia y potencia, tanto medida como simulada. Además, de
esta última se puede determinar, que la herramienta tiene una alta confiabilidad
para valores de irradiancia superiores a los 250 W/m^2, a partir de ese límite, la
potencia entregada por el sistema fotovoltaico modelado, se empieza a acercar
mucho más a la referencia que se tenía, la cual corresponde a la potencia medida
que entregó el sistema en el año de 2016. El método que se estableció con el
ajuste lineal, permitió pasar de un error porcentual absoluto medio de 41,6% a
28,6%, logrando reducir 13 puntos el error absoluto medio. Además, se puede
considerar un punto adicional para la herramienta, pues a pesar que ésta no fue
diseñada para realizar estimaciones de energía durante un largo periodo de
tiempo, en el capítulo cuatro queda evidenciado que es posible realizar un estudio
de este tipo, a pesar de no ser una herramienta práctica para este tipo de
utilidades.
El modelado y simulación para el caso de estudio propuesto de la Universidad
Autónoma de Occidente, se constituye en el principal aporte de este trabajo, dado
que integra fuentes de generación no convencionales a un sistema de distribución
típico que se tiene normalmente en usuarios industriales o comerciales. El cual
permite conocer el comportamiento de los flujos de potencia en una red, frente a
diferentes condiciones operativas que se tienen a lo largo de un día, con perfiles
de consumo diferentes. Los resultados que se obtuvieron para los casos bajo
estudio, permitieron corroborar mediante una herramienta computacional para la
simulación de sistemas de potencia, que efectivamente, si se están entregando
108
excedentes a la red los días festivos. Además, de estimar el porcentaje de energía
que se estaría cubriendo durante unos días típicos de consumo en la universidad.
BIBLIOGRAFIA
ALVAREZ, Raúl. San Francisco será la primera ciudad de EEUU en exigir paneles
solares en edificios nuevos [en línea]. Xataka. 20 de abril de 2016. 6 párr.
[Consultado: 7 de abril de 2017]. Disponible en Internet:
[Link]
exigir-paneles-solares-en-edificios-nuevos
CELSIA. Celsia pone en operación segunda fase de techos solares en Universidad
Autónoma de Occidente [en línea]. Sala de prensa, 30 de junio de 2017.
[Consultado: 20 julio de 2017]. Disponible en internet: [Link]
de-prensa/celsia-techos-solares-universidad-autonoma
COLOMBIA. CONGRESO DE LA REPUBLICA. Ley 1715 (13, mayo, 2014). Por
medio de la cual se regula la integración de las energías renovables no
convencionales al sistema energético nacional. [En línea]. upme. p 3. [Consultado:
20 de abril de 2017]. Disponible en internet:
[Link]
Colombia instalará 2.025 MW de energías renovables [en línea]. En: Revista
Eólica y del Vehículo Eléctrico. [Consultado: abril 1 de 2017]. 24 de marzo de
2017. Disponible en internet: [Link]
instalara-2-025-mw-de-energias-renovables/
FUNDACIÓN DE LA ENERGÍA DE LA COMUNIDAD DE Madrid (FENERCOM).
Guía básica de la generación distribuida. [En línea]. fenercom.65 p. [consultado: 2
de abril de 2017]. Disponible en Internet:
[Link]
[Link]
GEO Tutoriales. Error porcentual absoluto medio (MAPE) en un pronóstico de
demanda [en línea]. Blogs Gestión de Operaciones. 26 de Enero de 2015, párr. 1.
[consultado: 30 de julio de 2017]. Disponible en Internet:
[Link]
absoluto-medio-mape-en-un-pronostico-de-demanda/
109
GRUPO NAP. Energía Solar Fotovoltaica. [En línea]. Colegio Oficial de Ingenieros
de Telecomunicación. Madrid: Ibergraphi, 2002. P 20. [Consultado: 2 de junio de
2017]. Disponible en internet:
[Link]
IDEAM y PNUD presentan Inventario Nacional de Emisiones de Gases de Efecto
Invernadero. [En línea]. Programa de las Naciones Unidas para el Desarrollo
(UNPD). [Consultado: abril 1 de 2017]. Disponible en internet:
[Link]
ideam-y-pnud-presentan-inventario-nacional-de-emisiones-de-gases-de-efecto-
[Link]
JELLE, Bjorn. Building Integrated Photovoltaics: A Concise Description of the
Current State of the Art and Possible Research Pathways. [En línea]. MDPI,
Energies. 21 de septiembre de 2016. 30 p. [Consultado: 9 de abril de 2017].
Disponible en Internet: [Link]
MANTILLA, Juan Miguel, et. Al. Análisis del esquema de generación distribuida
como una opción para el sistema eléctrico colombiano [En línea]. En: Revista
facultad de ingeniería Universidad Antioquia. No. 44: junio de 2008. p. 97-110.
[Consultado: 3 de abril de 2017]. Disponible en Internet:
[Link]
MARTINEZ, Alberto. Energía solar fotovoltaica integrada en la edificación:
modelizado y análisis del efecto del sombreado en la captación de irradiación. [En
línea]. Doctorado en innovación en ingeniería de producto y procesos industriales.
La Rioja. Universidad de la Rioja. Facultad: Escuela Técnica Superior de
Ingeniería Industrial. Departamento de Ingeniería Mecánica, 2016. 317 p.
[Consultado: 12 de abril de 2017]. Disponible en Internet:
[Link]
REN21. Energía renovable 2016: Reporte de la situación mundial. [En línea].
ren21 32 p. [Consultado: 10 de abril de 2017]. Disponible en Internet:
[Link]
SANDOVAL, Ana María. Monografías del sector de electricidad y gas colombiano:
Condiciones actuales y retos futuros [en línea]. Archivos de economía. 2004. 121
p. [Consultado: 4 de abril de 2017]. Disponible en Internet:
[Link]
110
SIERRA, Enrique. Pruebas de correlaciones de Pearson y Spearman.
[Diapositivas] Scribd. [Consultado: 29 de julio de 2017]. Disponible en internet:
[Link]
spearman-dr-enrique-sierra
UNIDAD DE PLANEACIÓN DE MINAS Y ENERGÍA (UPME). Ministerio de Minas
y energías. [En línea]. Boletín estadístico de minas y energía 2012 – 2016.
Octubre de 2016. 201 p. [consultado: 30 de marzo de 2017]. Disponible en
Internet: [Link]
----.----. Registro de proyectos de generación Agosto 2016. [En línea]. [Link]
31 agosto 2017. 17 p. [Consultado: 30 de marzo de 2017]. Disponible en Internet:
[Link]
acion_Ago2016.pdf
----.----. Plan energético nacional Colombia: Ideario energético 2050. (Enero,
2015). [En línea]. upme p. 85. [Consultado: 31 de marzo de 2017]. Disponible en
Internet: [Link]
----.----. Capacidad instalada de autogeneración [en línea]. Consorcio HART-RE.
[Consultado: 21 de abril de 2017]. Disponible en internet:
[Link]
[Link]
UNIDAD DE PLANEACIÓN DE MINAS Y ENERGÍA (UPME). Capacidad instalada
de cogeneración 2014. [En línea]. Consorcio HART-RE. [Consultado: 21 de abril
de 2017]. Disponible en internet:
[Link]
[Link]
----.----. Costos indicativos de generación eléctrica en Colombia. [en línea].
Integral: Ingeniería de consulta. Abril de 2005. 250 p. [consultado: 5 de junio de
2017]. Disponible en internet:
[Link]
os_generacion_ee.pdf
WALPOLE, MYERS. Probabilidad y estadística para ingeniería y ciencias. [En
línea]. Novena edición. México: Pearson. p. 431. [Consultado 28 de julio de 2017].
Disponible en Internet:
111
[Link]
[Link]
112