Repositorio Institucional de Documentos

GRADO EN GEOLOGÍA AUTOR: PABLO VALDEARCOS PRUSÉN DIrectora: Ana Rosa Soria de Miguel Posibilidades del aprovechamiento geotérmico en las inmediaciones de Alfaro (La Rioja) Trabajo de fin de grado de geología (universidad de Zaragoza) Zaragoza, 8 de septiembre de 2014 Al tribunal de evaluación de los TFG 2013/2014 del Grado en Geología La Dra. Ana Rosa Soria de Miguel, directora del trabajo Fin de Grado presentado por Pablo Valdearcos Prusén y titulado “Posibilidades de aprovechamiento Geotérmico en las inmediaciones de Alfaro (La Rioja, España)” considera que dicho trabajo se ajusta a la normativa específica establecida a tal fin por el Departamento de Ciencias de la Tierra y por la Facultad de Ciencias de la Universidad de Zaragoza, por lo que doy mi aprobación para la exposición y defensa pública de dicho trabajo por parte de su autor. Atentamente Fdo: Ana Rosa Soria de Miguel 1 ÍNDICE 1 Resumen/Abstract 2 1. INTRODUCCIÓN 3 1.1. Objetivos y planteamiento del trabajo 3 1.2. Situación geográfica y geológica 3 1.3. Hidrogeología 5 1.4. Metodología 7 2. LA ENERGÍA GEOTÉRMICA 8 2.1. Energía geotérmica de muy baja temperatura 8 2.2. Captación de la energía geotérmica 9 3. ESTUDIOS PREVIOS 10 3.1. Necesidades de climatización de la vivienda 11 3.2. Estudio del terreno 11 3.3. Tipo de relleno del sondeo 12 3.4. Estudio del TRT (test de respuesta térmica) 13 4. CÁLCULO DEL APROVECHAMIENTO GEOTÉRMICO 19 4.1. Cálculo del aprovechamiento geotérmico de La Fonda 20 5. RENTABILIDAD DE LA INSTALACIÓN GEOTÉRMICA 25 6. CONCLUSIONS 27 7. REFERENCIAS BIBLIOGRAFICAS 28 8. ANEXOS 30 I. Resultados del TRT entre las 24 horas y las 26 horas y 40 mins 31 II. Relación de costes año a año entre la Energía geotérmica y la utilización de gas 34 2 Resumen La geotermia de muy baja temperatura es una energía que puede utilizarse en gran parte de la superficie terrestre, siendo una de las fuentes más baratas y eficientes para climatizar edificios. En este trabajo se muestra un estudio realizado en la localidad de Alfaro con el fin de determinar el rendimiento geotérmico del mismo y el tipo de sistema o instalación geotérmica más conveniente para el Hotel-restaurante (La Fonda) de esta localidad. La demanda energética del edificio se ha realizado mediante el análisis de los anteriores consumos energéticos de gas del edificio. El estudio geológico realizado mediante un sondeo inicial muestra que la litología del terreno es dominantemente arcillosa y con presencia de un acuífero a los 40m de profundidad. Tras el rellenado del pozo inicial se realiza un TRT o test de respuesta térmica, con el fin de analizar diversas propiedades del terreno tales como la conductividad, difusividad y resistencia térmica. Los datos obtenidos de todos estos estudios se introducen en un programa de ordenador (EED, Earth Energy Designer), juntos con otros datos y propiedades de la instalación, tales como el tipo de sondeo, tipo de fluido caloportador,... etc. A partir de todos estos parámetros, la instalación que requiere el edificio de La Fonda estará compuesta por 11 pozos de tipo U-simple de 103 metros de profundidad como mínimo y con 6 metros de espaciado entre cada uno. Finalmente se realiza el estudio de rentabilidad de la instalación geotérmica, amortizándose entre los 5 y 8 años de vida de la misma. Abstract The low temperature geothermal energy can be used quite a lot lengthwise the Earth surface, being one of the most economical energies to conditionate buildings. This study was carried out in the city of Alfaro to determine the geothermal efficiency of this area and the most advisable geothermal facility for the Hotel-Restaurant of the city (La Fonda). The energetical demand of the building has been determined by the analisis of previous gas consumptions. The geological study was carried out by an initial sounding. It shows that the main lithology in the land is chiefly clayey and that there is an aquifer 40m deep. After filling the initial sounding, a thermal respounse test (TRT) is carried out to determine different properties of the land, such as conductivity, difusivity or thermal resistance. The results obtained of these studies, as well as the information and properties of the facility, such as the type of sounding or the type of facility fluid, are inserted into a computer program (Earth Energy Designer, EED). Attending to all this information, the facility requiered in this case will be composed by 11 simple-U soundings, separated by 6 metres. This soundings will be at least 103metres deep. Eventually, the geothermal instalation rentability will be studied, being paid off through out the next 5-8 years. 3 1.- Introducción 1.1 Objetivos y planteamiento del trabajo El objetivo principal de este trabajo es estudiar las posibilidades del uso de la geotermia de muy baja temperatura o de baja entalpía para la climatización de edificios y uso de agua caliente. Este estudio se realiza mediante el análisis de una serie de parámetros de diferente naturaleza: hidrogeológicos, geológicos (litológicos), propiedades físicas del terreno tales como la conductividad térmica, clima,…etc. de la zona. A todos estos estudios, hay que añadirle las necesidades específicas de la vivienda que se va a climatizar, de tal manera que estas influirán de forma importante en el estudio, por lo que será una propiedad esencial a tener en cuenta. El trabajo de fin de grado que aquí se plantea se corresponde con un estudio de geología aplicada de carácter marcadamente interdisciplinar sobre el uso de la geotermia de muy baja temperatura para la climatización y uso de agua caliente sanitaria en una vivienda de la localidad de Alfaro (Logroño); concretamente, para el estudio de un hotel-restaurante der la localidad (La Fonda). Estos estudios pasan, en primer lugar, por conocer las necesidades de climatización o la demanda energética de la vivienda en concreto. Una vez conocidas estas necesidades, se procede a establecer las características del terreno en el que se va a realizar la Geotermia para determinar la capacidad de disipación de calor del terreno. A partir de este dato, se obtendrá la conductividad y la variación térmica de temperaturas que permite dicho terreno. La capacidad de disipación del terreno depende, principalmente, de la litología y de la presencia o no de acuífero; el cual, también condiciona el tipo de relleno en las perforaciones que deben realizarse. Todos estos datos, junto con las características climáticas de la zona en la que se ubica la vivienda a climatizar, propiedades de los materiales de instalación, etc., nos permitirán determinar la entidad de la instalación geotérmica que se precisa y si ésta es o no rentable. 1.2 Contexto geográfico y geológico La localidad de Alfaro se sitúa en la parte más oriental de La Rioja Baja. Se dispone a orillas del río Ebro, concretamente a la altura de la desembocadura del río Alhama. El área de estudio se sitúa dentro del Mapa geológico España a escala 1:50000, hoja nº 244 (Alfaro). En esta zona se pueden diferenciar toda una serie de materiales terciarios y cuaternarios, desde el Stampiense (Oligoceno), hasta las terrazas fluviales más modernas del Holoceno 4 (Figs.1.1 y 1.2). Estas terrazas fluviales se asientan sobre los materiales plegados de las sierras sur-pirenaicas navarras (situadas al norte y noreste) y la sierra de Cameros (situada al sur y suroeste). En Alfaro estas terrazas se disponen discordantemente sobre las Formaciones Lerín y Alfaro, que se identifican en el flanco sur del Anticlinal de Arguedas, una de las estructuras de la cadena sur-pirenaica Navarra (Fig. 1.3). Figura 1.1.- Cartografía geológica de la localidad de Alfaro y sus alrededores (Beroiz y Solé, 1977). Figura 1.2.- Leyenda del mapa geológico de Alfaro y sus alrededores (Beroiz y Solé, 1977). Las unidades presentes en los alrededores de la localidad de Alfaro son: - Formación Lerín: Comprende el Chatiense y el Aquitaniense inferior (Oligoceno superior-Mioceno inferior). Esta formación aparece a ambos lados del anticlinal de Arguedas con potencias que varían a ambos flancos del anticlinal. En el flanco sur, sobre el que se asientan las terrazas fluviales presentes en Alfaro, la potencia es de alrededor de 400 m. (Fig. 1.3). Litológicamente, esta formación tiene dos tramos litológicos distintos; uno más yesífero y otro más arcilloso. En la zona de estudio, en el flanco sur del anticlinal de T Ba-Bc C33-11 Formación Lerín 5 Arguedas, solo se aprecia el tramo más yesífero. En esta formación además, destacan algunos niveles de halita y otros más finos de areniscas y calizas intercalados con los yesos (Beroiz y Solé, 1977). El color o tono más característico de esta formación es el pardoamarillento. - Formación Alfaro: S dispone sobre las últimas capas de yeso de la Fm. Lerín. Abarca, desde el Aquitaniense inferior hasta el Vindoboniense, por lo que abarca casi todo el Mioceno. Esta formación, sin embargo, no aparece en el flanco norte del anticlinal de Arguedas, donde cambia de facies a la Fm. Tudela. Es complicado precisar la potencia de esta formación, debido a que no tenemos en ninguna zona el techo de la misma; no obstante, en la localidad de Alfaro se han medido hasta 50 m de potencia. (Fig. 1.3). Litológicamente, se compone de arcillas y limos que le aportan el tono rojizo a esta formación. Además aparecen intercalados niveles de areniscas con elevada proporción de granos de cuarzo y muy poco cementadas (Beroiz y Solé, 1977; Casas y Román, 1999). Figura 1.3.- Corte geológico N-S del flanco meridional del anticlinal de Arguedas (Beroiz y Solé, 1977) Los materiales cuaternarios de esta zona están representados por los sistemas de terrazas de los ríos Ebro y Alhama. Dentro de todo el sistema de terrazas del Ebro diferenciamos en los alrededores de Alfaro tres de ellas, situadas cronológicamente una a mediados del Pleistoceno, otra a finales y otra a comienzos del Holoceno (Beroiz y Solé, 1977). A su vez, también se distinguen las llanuras aluviales de los ríos ya mencionados y un sistema de glacis en la parte más oriental del pueblo, en la carretera de Castejón. 1.3 Hidrogeología En este apartado se expone la hidrogeología de la zona de estudio, a partir de los datos extraídos de diversos trabajos y estudios formalizados por la confederación hidrográfica del Ebro (CHE, 1995 a, b y c). La localidad de Alfaro se encuentra, hidrogeológicamente, dentro del dominio aluvial del Ebro y más concretamente dentro de la unidad hidrogeológica nº 405: “aluvial del Ebro: Lodosa-Tudela”, en su zona media (Fig. 1.4). En esta unidad hidrogeológica identifican las formaciones aluviales del río Ebro comprendidas en el tramo entre las localidades de Lodosa y Tudela así como los aluviales de sus afluentes Ega, Arga, Aragón, Cidacos y Alhama. Estos, junto con el Ebro, son los ríos que se relacionan con la unidad. En esta unidad se encuentra un único acuífero de carácter 6 libre, el aluvial del Ebro, constituido por las formaciones aluviales actuales y las terrazas bajas conectadas hidráulicamente con los cauces principales (terrazas de hasta 10 m). Las terrazas más altas, sin embargo, permanecen generalmente desconectadas, dando lugar a pequeños acuíferos colgados. Todos estos depósitos, descansan sobre las Fms. Alfaro y Lerín, constituidos principalmente por arcillas, yesos y margas, considerándolos como yacente de baja permeabilidad. Figura 1.4.- Cartografía hidrogeológica de la cuenca del río Ebro. En el recuadro está representada la unidad hidrogeológica en la que se encuentra la localidad de Alfaro (CHE, 1995c). Este acuífero varía de espesor con 10-20 m a la altura de Logroño, y llegando en la desembocadura del Alhama con el Ebro hasta 40 m, con espesores saturados comprendidos generalmente entre 7 y 25 m en la localidad de Alfaro (Fig. 1.5). No obstante, a medida que se asciende por el río Alhama, este espesor se va reduciendo hasta los 10-15 m de potencia. La geometría de la unidad es de perfil fusiforme, típica de acuíferos aluviales. La recarga de la unidad se realiza por infiltración del agua de lluvia y retornos de riego, además del almacenamiento en riberas en épocas de avenida, aportes de barrancos laterales y aportes subterráneos del aluvial aguas arriba de la unidad. Figura 1.5.- Cartografía de la U.H. Lodosa-Tudela. El círculo señalan la zona de estudio (CHE, 1995c). 7 Los recursos de esta unidad hidrogeológica se han estimado entre 135 y 206 hm3 al año. 1.4 Metodología A la hora de realizar un estudio de climatización en una vivienda o en cualquier recinto, deben seguirse una serie de rigurosos pasos con el fin de determinar, con la mayor fiabilidad posible, si es rentable una instalación de geotermia de baja temperatura en esa vivienda o recinto específico. En primer lugar es necesario conocer las necesidades de climatización de la vivienda. Este estudio, dependiendo de las características del edificio a climatizar, se puede hacer de forma manual, mediante un programa informático (Cype) o bien, a partir del gasto histórico del edificio, cuando se disponga de este dato. A continuación se realiza un sondeo inicial de 100 m de profundidad y se analiza y estudia la litología presente en la zona, así como la presencia o ausencia de acuíferos. Tras la realización de dicho estudio se procede al relleno del pozo con bentonita, si se carece de acuíferos someros o con gravas en presencia de acuíferos subsuperficiales. Después se realiza el Tests de Respuesta Térmica (TRT) en el sondeo inicial, cuyo objetivo es reproducir el comportamiento geotérmico de los diferentes horizontes litológicos frente a las inyecciones de calor en el terreno. Además según la normativa VDI-4640, en proyectos de más de 30 Kw, es obligatoria la realización de un test de respuesta térmico. Este Test, consiste en hacer circular entre 48 y 72 horas un fluido caloportador a través del pozo intercambiador de calor y monitorear la temperatura de entrada y salida del sistema. (Mattson et al., 2007). A partir de este test conoceremos una serie de parámetros como la potencia real que nos da el terreno, la conductividad térmica, la difusividad térmica y la resistencia térmica. Una vez conocidas las necesidades de la vivienda o potencia energética que precisa la vivienda, la litología, los datos físicos del terreno y la presencia o ausencia de acuíferos se procede a estudiar la posibilidad de la geotermia en el terreno. Este análisis se hace a partir de un programa informático denominado EED (Earth Energy designer), a partir del cual conoceremos las características de la estación geotérmica que cubra las necesidades del edificio (número de pozos que debemos hacer, profundidad de los mismos y la distancia o espaciado óptimo entre cada uno de ellos) (Blomberg et al., 2008). 14 partida previos (Fig. 3.4) que hacen referencia tanto a datos propios del sondeo (profundidad y diámetro), datos del tanque del que parte el agua y en el que es recogida, información del fluido y de las sondas, etc. Figura 3.4.- Conjunto de datos esenciales de partida para la realización de un TRT Una vez establecidos estos datos de partida ya se puede realizar es test. El TRT registra estos datos cada 60 minutos a lo largo de 75,78 horas. En Fig. 3.5 y Anexo 1 se muestran los registros de la maquina, en 11 columnas, que registran los siguientes parámetros: - Nº registro: muestra el número de registro que se está grabando - Tiempo: muestra el tiempo en segundos en el cual se ha tomado cada dato. Los datos se recogen en intervalos de 60 segundos. - Sonda 1: muestra la temperatura en el intercambiador donde actúan las resistencias. Esta temperatura no se tiene en cuenta para el cálculo, es meramente informativa. - Sonda 2: muestra la temperatura del fluido caloportador a la entrada del pozo. - Sonda 3: muestra la temperatura del fluido caloportador a la salida del pozo. - Sonda 4: muestra la temperatura exterior, temperatura ambiente. Esta temperatura se recoge solamente para comprobar que el ensayo no varía con la temperatura exterior. - Caudal: muestra el caudal que está pasando por el circuito. Este permanece relativamente constante, aún así se mide para tener dicho valor controlado. - Salto térmico: Diferencia entre las temperaturas de entrada y salida del pozo (S2 y S3). - Potencia calorífica: calcula la potencia disipada en el pozo. Lo calcula a partir del salto térmico y el caudal 15 - Temperatura media del fluido: calcula la media entre la entrada y la salida del pozo. - Ln t: calcula el logaritmo neperiano del tiempo. t (s) S1 (ºC) S2 (ºC) S3 (ºC) S4 (ºC) Caudal (l/min) Salto Térmico (S2-S3) Potencia Calor. (W) Tª media fluido (ºC) ln (t) 84000 25 26 23,4 13,1 18,7 2,60 3,38719333 24,70 11,34 84060 24,9 26 23,4 13,1 18,7 2,60 3,38719333 24,70 11,34 84120 24,8 26,1 23,4 13,1 18,8 2,70 3,53628 24,75 11,34 84180 24,9 26 23,4 13,1 18,7 2,60 3,38719333 24,70 11,34 84240 24,5 26,1 23,5 13,2 18,7 2,60 3,38719333 24,80 11,34 84300 25 26,1 23,4 13,2 18,8 2,70 3,53628 24,75 11,34 84360 25 26,1 23,4 13,2 18,7 2,70 3,51747 24,75 11,34 84420 24,9 26 23,4 13,2 18,7 2,60 3,38719333 24,70 11,34 84480 24,7 26 23,5 13,1 18,7 2,50 3,25691667 24,75 11,34 84540 24,4 26 23,5 13,2 18,7 2,50 3,25691667 24,75 11,34 84600 25 26,1 23,5 13,2 18,8 2,60 3,40530667 24,80 11,35 84660 25,2 26,1 23,5 13,2 18,7 2,60 3,38719333 24,80 11,35 84720 24,7 26,2 23,5 13,2 18,7 2,70 3,51747 24,85 11,35 84780 24,4 26 23,5 13,2 18,8 2,50 3,27433333 24,75 11,35 84840 25 26 23,5 13,2 18,7 2,50 3,25691667 24,75 11,35 84900 24,9 26,1 23,5 13,2 18,7 2,60 3,38719333 24,80 11,35 84960 24,6 26,1 23,5 13,3 18,7 2,60 3,38719333 24,80 11,35 85020 24,5 26,1 23,5 13,4 18,8 2,60 3,40530667 24,80 11,35 85080 24,6 26,2 23,5 13,4 18,7 2,70 3,51747 24,85 11,35 85140 25,3 26,1 23,5 13,5 18,7 2,60 3,38719333 24,80 11,35 85200 24,8 26,1 23,5 13,5 18,6 2,60 3,36908 24,80 11,35 Figura 3.5.- Datos registrados en el sondeo inicial de la fonda entre las 23,3 y las 23,6 horas desde el comienzo de la toma de datos De estos 10 tipos de datos (no se incluye el nº de registro), 6 de ellos son suministrados por el TRT a lo largo de las 75,78 horas, pero el resto son obtenidos a partir de una serie de cálculos. 16 Los datos propios del TRT son: el tiempo, la temperatura en el intercambiador donde actúan las resistencias (S1), las temperaturas de entrada y salida del pozo (S2 y S3), la temperatura ambiente (S4) y el caudal (el cual permanece relativamente constante). Los datos derivados y/o calculados con los datos registrados en el TRT son: el salto térmico, la potencia calorífica, la temperatura media del fluido en el pozo y el “ln” del tiempo. El salto térmico: El salto térmico es la diferencia entre la temperatura del agua a la entrada del pozo y la temperatura del agua en la salida del mismo  S3-S2. El sl salto térmico medio obtenido a lo largo del TRT realizado en La fonda es de entre 2,5 y 2,7 ºC La potencia calorífica: La potencia calorífica o potencia real del terreno, es la potencia o energía necesaria para aumentar la temperatura del fluido en una unidad de tiempo determinada. Este dato se calcula mediante una simple ecuación: P = Ce x Salto térmico x Caudal Donde P es la potencia calorífica del terreno, Ce es el calor específico del agua (valor constante), el salto térmico es obtenido a partir del TRT, y el caudal (a pesar de que es constante, se controlan sus medidas). Así pues la potencia calorífica media que se obtiene a lo largo del TRT realizado en La fonda varía entre 3,25 y 3,5 W. El “ln” del tiempo: Lo calcula directamente, obteniendo valores que varían de 11,34 a 11,35. La temperatura media del fluido en el pozo: Es la media aritmética del salto térmico de pozo: ((S2+S3)/2). Este dato dependerá en gran medida de la litología y sobretodo de la presencia de agua en el pozo. En el caso que nos ocupa varía entre 24.7 y 24.8. A partir de todos estos datos ya se pueden deducir diferentes propiedades del terreno como la conductividad térmica, la difusividad térmica y la resistencia térmica. Los cálculos y fórmulas utilizadas para la obtención de estos datos proceden del análisis y estudio del trabajo de Mattson et al. (2008). La conductividad térmica La conductividad térmica de un terreno es una propiedad física de los materiales que mide la capacidad de conducción de calor. En otras palabras, la conductividad térmica es también la capacidad de una sustancia de transferir la energía cinética de sus moléculas a otras moléculas adyacentes (Gálvez, 2013). 17 La conductividad térmica del terreno en el que se dispone La Fonda, la conocemos al realizar una gráfica donde se muestra el “ln” del tiempo y la temperatura media del fluido. En la gráfica resultante (Fig. 3.6) observamos la recta que marca la correlación entre ambos datos. La pendiente de esta recta nos dará el valor de la conductividad térmica del pozo. Figura 3.6.- Gráfica de correlación entre el ln del tiempo y la temperatura media del fluido a partir de los datos del TRT. Así, la pendiente de la línea resultante se ajusta a la fórmula: y = 1,433x + 8,402 A partir de esta fórmula se obtiene que la pendiente de la recta, que se corresponde con la conductividad térmica del terreno de La Fonda, es de 1,9291 W/mK Las temperaturas medias, necesarias para el cálculo de la conductividad térmica, entre las aguas de entrada y salida de un pozo, dependen de la litología presente y de la presencia de agua de acuíferos o no. Así pues, las conductividades térmicas resultantes de cada litología, en contacto o no con un acuífero, vienen representadas en la figura 3.7 (Pérez, 2010). Litología Humedad Arcillas Arcillas con abundantes niveles de areniscas Arcillas con yesos Limos con fragmentos de roca metamórfica Secas 0,5 1,5 1,2 0,8 En contacto con un acuífero 1,5-2,3 2,7-3 1,5-2,3 Figura 3.7.- Conductividades térmicas resultantes dependiendo del tipo de litología y la presencia o ausencia de acuíferos. Datos expresados en W/mK. 18 En nuestro caso, la litología que más se repite en el terreno es el de la arcilla, en ocasiones con areniscas (en el tramo superior) y otras veces con yesos (en el tramo inferior). Así pues teniendo en cuenta las variaciones del nivel freático, y la litología presente en esta zona, es lógico que el valor de la conductividad térmica en el terreno sea de alrededor de 2 W/mK (coincidiendo el valor obtenido con la media propuesta para esta litología en presencia de acuíferos). Difusividad térmica Es un índice que expresa la velocidad de cambio y flujo de temperaturas en un material hasta que alcanza el equilibrio térmico (David, 2009) Así pues la difusividad térmica de los materiales del terreno sobre el que se asienta La Fonda se calcula mediante la siguiente ecuación: α = λ / p*ce Donde α es la difusividad térmica, λ es la conductividad térmica (1,9291 W/mK), p es la densidad del agua (1000 kg/m3) y ce es el calor específico del agua (4180 J/KgºC). Así pues aplicando esta fórmula y teniendo en cuenta el valor de la conductividad térmica obtenido anteriormente se calcula la difusividad térmica del terreno de la Fonda, que va a ser de 4,615E-7 m2/s Resistencia térmica Es la capacidad de los materiales para oponerse al paso del calor. Es la magnitud inversa de la conductividad térmica (Gálvez, 2013) La resistencia térmica del terreno de La Fonda, se calcula mediante una compleja fórmula (Fig. 3.8): Figura 3.8.- Fórmula de la resistividad térmica de un terreno Esta fórmula depende del salto térmico (Tf-T0), la conductividad térmica del terreno (λ), el radio del pozo (rb), la inyección de calor constante en el terreno (qc), el tiempo (t), la 19 difusividad térmica (α), la constante de Euler (γ) y por la resistencia térmica (Rb )(Mattson et al., 2008) Así, mediante la aplicación de esta fórmula en nuestro terreno, obtenemos el valor de la resistencia térmica para los materiales que lo forman, que es 0,0709 K/Wm 4.- Cálculo del aprovechamiento geotérmico Una vez que ya tenemos todos y cada uno de los datos necesarios para realizar el cálculo del aprovechamiento geotérmico de la localidad de Alfaro, tanto la información de la vivienda (necesidad energética del edificio o el estudio de la climatización de la vivienda), como toda la información del terreno (datos geológicos y geofísicos), podemos proceder a calcular las características que deberá tener la instalación de energía geotérmica para satisfacer las necesidades del edificio a climatizar. Las características de la instalación vienen definidas por el número de pozos que se necesitas, su profundidad y la distancia que debe existir entre cada uno de ellos (Blomberg et al., 2008). Para realizar este cálculo, introduciremos todos y cada uno de los datos registrados y calculados en un programa de ordenador conocido como EED (“EarthEnergyDesigner”) (Fig. 4.1). EED es un programa de PC para el diseño de sonda geotérmica vertical. Su facilidad de uso, corta curva de aprendizaje, los tiempos de cálculo y bases de datos rápidas inherentes, hacen del programa EED una herramienta de elevada utilidad en el campo de la geotermia para el diseño de un sistema de bomba de calor geotérmica (GSHP) como el del hotel-restaurante de La Fonda En este programa podemos comprobar cómo determinados factores y datos influirán en mayor medida a la hora de calcular las profundidades de cada pozo; otros, sin embargo, nos ayudarán a determinar la distancia entre cada uno de ellos, y otros datos serán utilizados principalmente para determinar el nº de pozos. Figura 4.1.- Logotipo del programa EED versión 3.0, utilizado en el cálculo de esta instalación. 20 En cuanto a los datos geológicos registrados, la litología del terreno nos ayudará principalmente a conocer la profundidad de cada pozo y el número de pozos necesarios. No obstante, la presencia o ausencia de acuífero en la zona será un gran condicionante de la conductividad térmica del pozo y la difusividad del mismo. Esta información será necesaria para calcular la distancia que debe existir entre cada uno de los pozos. La conductividad térmica también influye en el número de pozos. 4.1 Cálculo del aprovechamiento geotérmico de La Fonda En primer lugar se introducen en el programa las propiedades del suelo o del terreno (Fig. 4.2 A). En este apartado será necesario introducir la conductividad térmica hallada (1,929 W/mK), la capacidad calorífica volumétrica del relleno (para la bentonita es 2,16 M3/(m3 * K), la temperatura de la superficie de la tierra (temperatura ambiente) que es 18,3 ºC y el calor geotérmico del fluido dependiente del lugar de La Tierra en el que nos encontremos. Para toda la cuenca del Ebro se coge el valor de Zaragoza, por lo que este valor es de 0,07 W/m2. En segundo lugar se deben introducir las características del pozo o sondeo inicial y del intercambiador o tanque de calor del TRT (Fig. 4.2 B). En este apartado primeramente se indica el tipo de sonda a utilizar en los sondeos (U-Simple) y su configuración dentro del programa (283), a continuación se expone la profundidad del sondeo inicial (100 m) y el espaciado entre pozos, si es que se tiene alguna preferencia, que no es nuestro caso, por lo que se deja el espaciado que indica el programa como estandar. Después se indica el diámetro del pozo (110 mm). A continuación se indica el estado de contacto entre el tubo y el fluido y se indica, si existe, la presencia de burbujas dentro del tubo. En este cálculo y en la mayoría se supone un contacto perfecto con una resistencia de 0 mK/W. Luego se indica la conductividad térmica del tipo de relleno a utilizar (para la bentonita es 0,6 W/mk). También en este apartado se muestran las características del tubo o sonda (Fig. 4.2 B). De modo que se indican el diámetro exterior (40 mm), el grosor de la pared (3,7 mm), la conductividad térmica del material de la sonda (en nuestro caso se utiliza polietileno, cuya conductividad térmica es de 0,420 W/mK y el espaciado entre los tubos de entrada y de salida (en nuestros sondeos es de 70 mm). Seguidamente se indica la resistencia térmica del pozo (Fig. 4.2 C), cuyo valor lo hemos obtenido, para la zona de estudio, mediante el análisis del TRT: 0,0709 mK/W. Posteriormente se introducen las características y propiedades del fluido a utilizar como fluido caloportador de la instalación. En el programa simplemente se debe indicar el tipo de 21 fluido que se va a utilizar. En nuestro caso utilizaremos agua con glicol cuyas características se muestran en la Fig. 4.2 D. A B C D E 22 F Figura 4.2.- Apartados del programa EED versión 3.0 para la obtención de una correcta y productiva instalación de geotermia en La Fonda. A.- Propiedades del terreno. B.- Pozo e intercambiador o tanque de calor. C.- Resistencia térmica del pozo. D.- Calor fluido portador. E.- Carga base. F.- Carga pico. Seguidamente se introducen los datos de la carga base, es decir, la potencia media máxima que se puede requerir el edificio en cada mes del año a partir del cual se muestra la demanda energética o necesidad de climatización del recinto por meses. Además se expone el calor necesario y el agua caliente doméstica (ACD) anuales. Esta carga base en relación a cada mes del año viene indicada en el gráfico de la Fig. 4.2 E. En el último apartado se muestran la carga pico o potencia máxima que puede requerir el edificio o recinto en un momento dado en cada uno de los meses del año (Fig. 4.2 F). Así comprobamos que teniendo en cuenta la demanda energética del edificio (130 MWh/año) y la potencia de la caldera existente (80 Kw), dicho pico asciende hasta 50 Kw en todos los meses, menos en los de verano. Una vez introducidos todos los datos, se especifican las limitaciones que puede tener el edificio o recinto donde se va a realizar la instalación (Fig. 4.3). En nuestro caso, para La Fonda, ponemos el área que disponemos para la perforación de los pozo (300 x 200 m2). El espaciado entre pozos queremos que sea de 6 metros, ya que para esta zona es la más óptima teniendo en cuenta la litología y propiedades de los materiales y el nivel freático de la zona. La profundidad de los pozos no puede superar los 125 metros, no obstante, deben superar los 100 m para obtener un buen rendimiento de los mismos. Y no se desean realizar más de 11 pozos teniendo en cuenta, principalmente, la accesibilidad de maquinaria a la zona (Fernando Baños, com. pers.). D 23 Figura 4.3.- Apartado final de la Optimización EED. Datos del programa EED versión 3.0 A partir de todos los datos que hemos ido introduciendo, el programa nos ofrece las distintas posibilidades de instalación y, teniendo en cuenta que la vida útil de una instalación de geotermia es de 25 años como mínimo, elegimos aquella más adecuada. Además, cada opción ofrece 2 gráficos; uno donde se expone la temperatura del fluido máxima y mínima en el mes de menor aprovechamiento geotérmico (normalmente el último, el año 25) (Fig. 4.4) y otro donde se expone la temperatura máxima y mínima del fluido a lo largo de los 25 años de vida de la instalación (Fig. 4.5). Figura 4.4.- Gráfica de la Tª del fluido en el año de menor productividad de la opción elegida La elección entre una de las dos opciones que nos proporciona el programa será función de diversos aspectos (Fernando Baños, com. pers.): 30 8.- ANEXOS 86400 25 26,2 23,5 13,8 18,6 2,70 3,49866 24,85 11,37 86460 25,2 26,2 23,5 13,7 18,7 2,70 3,51747 24,85 11,37 86520 24,7 26,2 23,5 13,7 18,6 2,70 3,49866 24,85 11,37 86580 25,2 26,2 23,5 13,7 18,7 2,70 3,51747 24,85 11,37 86640 25 26,2 23,5 13,7 18,7 2,70 3,51747 24,85 11,37 86700 25,2 26,2 23,5 13,6 18,7 2,70 3,51747 24,85 11,37 86760 24,9 26,2 23,5 13,7 18,7 2,70 3,51747 24,85 11,37 86820 24,9 26,2 23,5 13,7 18,7 2,70 3,51747 24,85 11,37 86880 25 26,2 23,5 13,7 18,7 2,70 3,51747 24,85 11,37 86940 24,8 26,1 23,5 13,8 18,8 2,60 3,40530667 24,80 11,37 87000 24,9 26,2 23,6 13,8 18,7 2,60 3,38719333 24,90 11,37 87060 25,4 26,1 23,5 13,8 18,7 2,60 3,38719333 24,80 11,37 87120 25,2 26,2 23,6 13,8 18,8 2,60 3,40530667 24,90 11,38 87180 25 26,1 23,6 13,7 18,7 2,50 3,25691667 24,85 11,38 87240 24,9 26,2 23,6 13,7 18,7 2,60 3,38719333 24,90 11,38 87300 25,4 26,2 23,6 13,8 18,7 2,60 3,38719333 24,90 11,38 87360 25,5 26,1 23,5 13,8 18,7 2,60 3,38719333 24,80 11,38 87420 25,1 26,2 23,6 13,8 18,7 2,60 3,38719333 24,90 11,38 87480 25,5 26,2 23,5 13,8 18,7 2,70 3,51747 24,85 11,38 87540 24,8 26,2 23,6 13,8 18,7 2,60 3,38719333 24,90 11,38 87600 25 26,2 23,5 13,8 18,7 2,70 3,51747 24,85 11,38 87660 25 26,1 23,5 13,8 18,7 2,60 3,38719333 24,80 11,38 87720 24,8 26,2 23,5 13,9 18,7 2,70 3,51747 24,85 11,38 87780 25,1 26,2 23,5 13,8 18,7 2,70 3,51747 24,85 11,38 87840 25,3 26,1 23,6 13,9 18,8 2,50 3,27433333 24,85 11,38 87900 25,9 26,1 23,6 13,9 18,7 2,50 3,25691667 24,85 11,38 87960 24,8 26,2 23,5 14 18,7 2,70 3,51747 24,85 11,38 88020 25,2 26,2 23,6 14 18,8 2,60 3,40530667 24,90 11,39 88080 25,6 26,2 23,6 14 18,8 2,60 3,40530667 24,90 11,39 88140 25,6 26,2 23,6 14 18,8 2,60 3,40530667 24,90 11,39 88200 24,6 26,1 23,6 14 18,8 2,50 3,27433333 24,85 11,39 88260 25,6 26,1 23,6 14 18,8 2,50 3,27433333 24,85 11,39 88320 25,2 26,2 23,6 14 18,8 2,60 3,40530667 24,90 11,39 88380 24,9 26,1 23,6 14 18,8 2,50 3,27433333 24,85 11,39 88440 25,2 26,2 23,6 13,9 18,7 2,60 3,38719333 24,90 11,39 88500 25 26,2 23,6 14 18,8 2,60 3,40530667 24,90 11,39 88560 24,8 26,1 23,6 14 18,8 2,50 3,27433333 24,85 11,39 88620 25,4 26,1 23,6 14 18,8 2,50 3,27433333 24,85 11,39 88680 24,9 26,1 23,6 14 18,7 2,50 3,25691667 24,85 11,39 88740 25,4 26,1 23,6 14 18,8 2,50 3,27433333 24,85 11,39 88800 24,7 26,2 23,6 14 18,8 2,60 3,40530667 24,90 11,39 88860 25,3 26,2 23,6 14 18,7 2,60 3,38719333 24,90 11,39 88920 25,1 26,3 23,6 14 18,7 2,70 3,51747 24,95 11,40 88980 25,3 26 23,6 14 18,7 2,40 3,12664 24,80 11,40 89040 25,4 26,1 23,6 14 18,7 2,50 3,25691667 24,85 11,40 89100 24,6 26,2 23,6 14 18,7 2,60 3,38719333 24,90 11,40 89160 25 26,2 23,6 14 18,8 2,60 3,40530667 24,90 11,40 89220 24,6 26,2 23,6 14 18,8 2,60 3,40530667 24,90 11,40 89280 25 26,2 23,6 14 18,7 2,60 3,38719333 24,90 11,40 89340 24,5 26,2 23,6 14 18,6 2,60 3,36908 24,90 11,40 89400 25,5 26,2 23,6 14,1 18,8 2,60 3,40530667 24,90 11,40 89460 25,3 26,1 23,6 14,1 18,7 2,50 3,25691667 24,85 11,40 89520 25,6 26,2 23,6 14,2 18,7 2,60 3,38719333 24,90 11,40 89580 24,8 26,2 23,6 14,2 18,8 2,60 3,40530667 24,90 11,40 ANEXO I. Resultados del TRT entre las 24 horas y las 26 horas y 40 mins t (s) S1 (ºC) S2 (ºC) S3 (ºC) S4 (ºC) Caudal (l/min) Salto Térmico (S2-S3) Potencia calorífica (W) Temperatura media fluido Tf (ºC) ln (t) 40 minutos 89640 25,1 26,1 23,6 14,2 18,7 2,50 3,25691667 24,85 11,40 89700 25,2 26,1 23,6 14,3 18,7 2,50 3,25691667 24,85 11,40 89760 24,6 26,2 23,6 14,3 18,7 2,60 3,38719333 24,90 11,40 89820 25 26,1 23,6 14,3 18,8 2,50 3,27433333 24,85 11,41 89880 25 26,1 23,6 14,3 18,8 2,50 3,27433333 24,85 11,41 89940 24,9 26,1 23,6 14,4 18,7 2,50 3,25691667 24,85 11,41 90000 25,1 26,2 23,6 14,5 18,8 2,60 3,40530667 24,90 11,41 90060 24,8 26,2 23,6 14,4 18,7 2,60 3,38719333 24,90 11,41 90120 25,8 26,1 23,6 14,4 18,8 2,50 3,27433333 24,85 11,41 90180 25,1 26,2 23,6 14,4 18,8 2,60 3,40530667 24,90 11,41 90240 25,1 26,1 23,6 14,4 18,7 2,50 3,25691667 24,85 11,41 90300 25,2 26,1 23,6 14,4 18,8 2,50 3,27433333 24,85 11,41 90360 25,3 26,2 23,6 14,4 18,7 2,60 3,38719333 24,90 11,41 90420 24,7 26,1 23,6 14,3 18,7 2,50 3,25691667 24,85 11,41 90480 24,9 26,2 23,6 14,3 18,7 2,60 3,38719333 24,90 11,41 90540 24,8 26,1 23,6 14,3 18,8 2,50 3,27433333 24,85 11,41 90600 24,4 26,2 23,6 14,3 18,8 2,60 3,40530667 24,90 11,41 90660 24,9 26,1 23,6 14,2 18,8 2,50 3,27433333 24,85 11,41 90720 24,7 26,1 23,6 14,2 18,7 2,50 3,25691667 24,85 11,42 90780 24,6 26,1 23,6 14,2 18,8 2,50 3,27433333 24,85 11,42 90840 24,6 26,1 23,6 14,2 18,8 2,50 3,27433333 24,85 11,42 90900 24,9 26,1 23,6 14,2 18,7 2,50 3,25691667 24,85 11,42 90960 25,2 26 23,6 14,2 18,8 2,40 3,14336 24,80 11,42 91020 24,7 26,1 23,6 14,2 18,7 2,50 3,25691667 24,85 11,42 91080 24,7 26,2 23,6 14,2 18,8 2,60 3,40530667 24,90 11,42 91140 24,7 26,1 23,6 14,1 18,8 2,50 3,27433333 24,85 11,42 91200 25,2 26,1 23,6 14,1 18,8 2,50 3,27433333 24,85 11,42 91260 24,6 26,2 23,6 14,1 18,7 2,60 3,38719333 24,90 11,42 91320 24,6 26,1 23,6 14,1 18,7 2,50 3,25691667 24,85 11,42 91380 25 26,2 23,6 14 18,7 2,60 3,38719333 24,90 11,42 91440 25,2 26,1 23,6 14 18,6 2,50 3,2395 24,85 11,42 91500 25 26,1 23,5 14 18,6 2,60 3,36908 24,80 11,42 91560 24,9 26,1 23,6 14 18,7 2,50 3,25691667 24,85 11,42 91620 24,7 26,2 23,6 14 18,7 2,60 3,38719333 24,90 11,43 91680 24,5 26,2 23,6 13,9 18,7 2,60 3,38719333 24,90 11,43 91740 25 26,1 23,5 13,9 18,7 2,60 3,38719333 24,80 11,43 91800 24,8 26,1 23,6 13,9 18,6 2,50 3,2395 24,85 11,43 91860 25,2 26,1 23,6 13,9 18,8 2,50 3,27433333 24,85 11,43 91920 25,1 26 23,6 13,9 18,7 2,40 3,12664 24,80 11,43 91980 25,3 26,2 23,6 13,9 18,7 2,60 3,38719333 24,90 11,43 92040 24,7 26,1 23,6 13,8 18,7 2,50 3,25691667 24,85 11,43 92100 25,2 26,1 23,6 13,8 18,7 2,50 3,25691667 24,85 11,43 92160 25,6 26,1 23,6 13,8 18,7 2,50 3,25691667 24,85 11,43 92220 25,7 26,1 23,6 13,8 18,8 2,50 3,27433333 24,85 11,43 92280 24,7 26,2 23,6 13,8 18,8 2,60 3,40530667 24,90 11,43 92340 24,7 26,2 23,6 13,8 18,8 2,60 3,40530667 24,90 11,43 92400 24,5 26,1 23,6 13,8 18,8 2,50 3,27433333 24,85 11,43 92460 24,8 26,2 23,5 13,8 18,8 2,70 3,53628 24,85 11,43 92520 25,4 26,1 23,6 13,8 18,7 2,50 3,25691667 24,85 11,44 92580 25 26,1 23,6 13,7 18,8 2,50 3,27433333 24,85 11,44 92640 24,8 26,2 23,6 13,7 18,7 2,60 3,38719333 24,90 11,44 92700 24,8 26,1 23,6 13,7 18,8 2,50 3,27433333 24,85 11,44 92760 25 26,1 23,6 13,7 18,7 2,50 3,25691667 24,85 11,44 92820 24,8 26,2 23,6 13,6 18,7 2,60 3,38719333 24,90 11,44 92880 25 26,2 23,6 13,6 18,7 2,60 3,38719333 24,90 11,44 92940 24,6 26,2 23,6 13,5 18,7 2,60 3,38719333 24,90 11,44 93000 25 26,2 23,6 13,5 18,8 2,60 3,40530667 24,90 11,44 93060 24,7 26,2 23,6 13,6 18,7 2,60 3,38719333 24,90 11,44 93120 24,7 26,2 23,6 13,6 18,7 2,60 3,38719333 24,90 11,44 93180 24,5 26,1 23,6 13,6 18,7 2,50 3,25691667 24,85 11,44 93240 24,8 26,2 23,6 13,6 18,8 2,60 3,40530667 24,90 11,44 93300 24,8 26,2 23,6 13,6 18,7 2,60 3,38719333 24,90 11,44 93360 24,6 26,2 23,6 13,6 18,7 2,60 3,38719333 24,90 11,44 93420 25,1 26,2 23,6 13,6 18,7 2,60 3,38719333 24,90 11,44 93480 24,5 26,1 23,6 13,5 18,7 2,50 3,25691667 24,85 11,45 93540 24,6 26,1 23,6 13,5 18,8 2,50 3,27433333 24,85 11,45 93600 24,8 26,2 23,6 13,5 19 2,60 3,44153333 24,90 11,45 93660 24,7 26,2 23,6 13,5 18,7 2,60 3,38719333 24,90 11,45 93720 24,8 26,2 23,6 13,5 18,7 2,60 3,38719333 24,90 11,45 93780 24,6 26,2 23,6 13,4 18,7 2,60 3,38719333 24,90 11,45 93840 25,2 26,3 23,6 13,4 18,6 2,70 3,49866 24,95 11,45 93900 24,4 26,2 23,6 13,4 18,7 2,60 3,38719333 24,90 11,45 93960 24,5 26,2 23,6 13,4 18,7 2,60 3,38719333 24,90 11,45 94020 25,3 26,1 23,6 13,4 18,7 2,50 3,25691667 24,85 11,45 94080 24,7 26,2 23,6 13,4 18,7 2,60 3,38719333 24,90 11,45 94140 24,9 26,1 23,6 13,4 18,7 2,50 3,25691667 24,85 11,45 94200 24,5 26,2 23,6 13,3 18,7 2,60 3,38719333 24,90 11,45 94260 25,2 26,2 23,6 13,3 18,8 2,60 3,40530667 24,90 11,45 94320 24,6 26,1 23,6 13,3 18,8 2,50 3,27433333 24,85 11,45 94380 25,5 26,2 23,6 13,3 18,7 2,60 3,38719333 24,90 11,46 94440 24,8 26,2 23,6 13,3 18,7 2,60 3,38719333 24,90 11,46 94500 25,5 26,1 23,6 13,3 18,7 2,50 3,25691667 24,85 11,46 94560 25 26,2 23,6 13,2 18,7 2,60 3,38719333 24,90 11,46 94620 24,7 26,1 23,6 13,2 18,8 2,50 3,27433333 24,85 11,46 94680 25 26 23,6 13,2 18,7 2,40 3,12664 24,80 11,46 94740 24,7 26,1 23,6 13,2 18,7 2,50 3,25691667 24,85 11,46 94800 24,9 26,1 23,6 13,2 18,7 2,50 3,25691667 24,85 11,46 94860 24,7 26,2 23,6 13,1 18,8 2,60 3,40530667 24,90 11,46 94920 24,5 26,1 23,6 13,1 18,8 2,50 3,27433333 24,85 11,46 94980 24,6 26,2 23,6 13,1 18,7 2,60 3,38719333 24,90 11,46 95040 24,8 26 23,6 13,1 18,8 2,40 3,14336 24,80 11,46 95100 25 26,2 23,6 13 18,7 2,60 3,38719333 24,90 11,46 95160 24,5 26,1 23,6 12,9 18,7 2,50 3,25691667 24,85 11,46 95220 24,7 26,1 23,6 12,9 18,8 2,50 3,27433333 24,85 11,46 95280 24,8 26 23,6 13 18,8 2,40 3,14336 24,80 11,46 95340 24,5 26 23,6 13 18,7 2,40 3,12664 24,80 11,47 95400 24,7 25,9 23,6 13 18,8 2,30 3,01238667 24,75 11,47 95460 25 26 23,6 13 18,7 2,40 3,12664 24,80 11,47 95520 24,4 26 23,6 12,9 18,7 2,40 3,12664 24,80 11,47 95580 25,4 26 23,6 12,9 18,7 2,40 3,12664 24,80 11,47 95640 24,4 26,1 23,6 12,9 18,7 2,50 3,25691667 24,85 11,47 95700 24,7 26 23,6 12,9 18,7 2,40 3,12664 24,80 11,47 95760 24,5 26 23,6 12,9 18,8 2,40 3,14336 24,80 11,47 95820 24,7 26,1 23,6 12,9 18,8 2,50 3,27433333 24,85 11,47 95880 24,7 26,1 23,6 12,9 18,7 2,50 3,25691667 24,85 11,47 95940 24,4 26,1 23,6 12,8 18,7 2,50 3,25691667 24,85 11,47 96000 24,4 26 23,6 12,7 18,7 2,40 3,12664 24,80 11,47 ANEXO II. Relación de costes año a año entre la Energía geotérmica y la utilización de gas COSTE ANUAL E. GEOTERMIA COSTE ANUAL GAS OIL COSTE ANUAL GAS COSTE ANUAL REFRIGER. COSTE ACUM. E. GEOTERMIA COSTE ACUM. GAS OIL + ELECTR. COSTE ACUM. GAS + ELECTR. AHORRO ACUMULADO COMPARADO CON GASOIL Año0 130.000,00 50.000,00 50.000,00 0,00 130.000,00 50.000,00 50.000,00 -80.000,00 Año1 8.612,20 23.995,41 16.796,79 0,00 138.612,20 73.995,41 66.796,79 -64.616,79 Año2 8.827,50 25.075,20 17.552,64 0,00 147.439,70 99.070,61 84.349,43 -48.369,09 Año3 9.048,19 26.203,59 18.342,51 0,00 156.487,89 125.274,20 102.691,94 -31.213,69 Año4 9.274,40 27.382,75 19.167,92 0,00 165.762,29 152.656,95 121.859,86 -13.105,34 Año5 9.506,26 28.614,97 20.030,48 0,00 175.268,54 181.271,92 141.890,35 6.003,38 Año6 9.743,91 29.902,65 20.931,85 0,00 185.012,46 211.174,57 162.822,20 26.162,11 Año7 9.987,51 31.248,27 21.873,79 0,00 194.999,97 242.422,83 184.695,98 47.422,87 Año8 10.237,20 32.654,44 22.858,11 0,00 205.237,16 275.077,27 207.554,09 69.840,11 Año9 10.493,13 34.123,89 23.886,72 0,00 215.730,29 309.201,16 231.440,81 93.470,87 Año10 10.755,46 35.659,46 24.961,62 0,00 226.485,75 344.860,62 256.402,43 118.374,87 Año11 11.024,34 37.264,14 26.084,90 0,00 237.510,09 382.124,76 282.487,33 144.614,67 Año12 11.299,95 38.941,02 27.258,72 0,00 248.810,04 421.065,78 309.746,05 172.255,74 Año13 11.582,45 40.693,37 28.485,36 0,00 260.392,49 461.759,15 338.231,41 201.366,66 Año14 11.872,01 42.524,57 29.767,20 0,00 272.264,50 504.283,72 367.998,61 232.019,22 Año15 12.168,81 44.438,18 31.106,72 0,00 284.433,31 548.721,90 399.105,33 264.288,59 Año16 12.473,03 46.437,90 32.506,53 0,00 296.906,34 595.159,80 431.611,86 298.253,45 Año17 12.784,86 48.527,60 33.969,32 0,00 309.691,20 643.687,40 465.581,18 333.996,20 Año18 13.104,48 50.711,34 35.497,94 0,00 322.795,68 694.398,74 501.079,12 371.603,06 Año19 13.432,09 52.993,35 37.095,35 0,00 336.227,77 747.392,09 538.174,47 411.164,32 Año20 13.767,89 55.378,05 38.764,64 0,00 349.995,66 802.770,15 576.939,10 452.774,48