Este trabajo presenta el estudio del acoplamiento de una planta solar con ciclo combinado. Se estudian diferentes configuraciones considerando la planta solar como el economizador o el sobrecalentador de la caldera de recuperación de calor. El objetivo es obtener el diseño óptimo, desde el punto de vista termoeconómico de la caldera, incluyendo la planta solar, determinando los parámetros de diseño optimizados para ambos sistemas. Se aplica una metodología empleada en trabajos previos para la optimización de ciclos combinados, pero incluyendo ahora la planta solar. Asimismo, se realiza un análisis de sensibilidad con respecto a la variación de radiación solar, tomando un día promedio del mes de mayo. Como resultado se obtiene el rendimiento y costo óptimos de las configuraciones analizadas.
This paper presents the study of the coupling of a solar plant with a combined cycle. Different configurations are analyzed considering the solar plant as economizer or superheated boiler with heat recovery. The goal is to obtain the thermoeconomic optimal design of the boiler, including the solar plant, so determining the optimized design parameters for both systems. The methodology used in previous papers to optimize the combined cycle was used, but now including the solar plant. Also, a sensitivity analysis with respect to the variation of solar radiation is conducted, taking the average day of May. As a result, the best performance and cost of the obtained configurations is analyzed.
Área de apertura del colector solar m2
Área de transferencia de calorm2
Approach point K
Alta presión
Baja presión
Razón de concentración
Costo M€
Ciclo combinado
Calor específico a presión constante kJ/kg K
Costo Total M€
Costo de generación por kWh €/kWh
Condensador
Caldera de Recuperación de Calor
Desgasificador
Diferencia terminal de temperaturas en el sobrecalentador K
Economizador
Evaporador
Factor de remoción de calor
Factor de eficiencia del colector
Irradiancia directa instantánea W/m2
Generación directa de Vapor
Horas de operación de la planta h
Heat Recovery Steam Generator
Integrated solar combined cycle
Constate de costo correspondiente a una sección de la CRC Flujo másico kg/s
Pinch Point K
Calor útil W
Sobrecalentador
Temperatura K
Turbina de vapor
Turbina de gas
Coeficiente global de transferencia de calor del colector W/m2K
Potencia MW Potencia media anual de la planta MW
Absortancia
Eficiencia térmica
Eficiencia óptica
Transmitancia
Interior
Ambiente
Alta presión
Baja presión
Sección de la caldera de recuperación de calor
La electricidad generada a partir de energía solar-térmica (Solar Thermal Electricity, STE) es una tecnología probada con más de 30 años de experiencia esperando que sea competitiva con las centrales térmicas de combustible fósil para el 2020 (Llorente et al., 2011; IEA, 2011). Esta tecnología está disponible principalmente a partir de plantas tipo CSP (Concentrating Solar Power) instaladas en regiones áridas o semi-áridas.
Actualmente existen cuatro tecnologías comerciales disponibles de sistemas de concentración solar (CSP); dos de ellas, el canal parabólico y el concentrador Fresnel linear reflejan la radiación solar hacia un concentrador lineal, mientras que las otras dos la concentran en un punto.
A finales de 2010, había 1318 MW de capacidad instalada acumulada de CSP en todo el mundo. Hasta ese año, España fue el líder mundial en instalaciones termosolares, con 450 MW de capacidad adicional, acumulando 55.4% de la capacidad instalada en todo el mundo. Mientras tanto, Estados Unidos añadió 78 MW de potencia CSP, para un total de 38.5% de la capacidad acumulada mundial, mientras que otros países como: Irán (5.0%), Israel (0.5%), Australia (0.2%) y Alemania (0.1%) juntos sumaron un porcentaje menor (U.S. DE, 2011).
La tecnología CSP del tipo canal parabólico es la única que ha demostrado su funcionamiento comercial a largo plazo, esto se debe a su corto tiempo de puesta en marcha y gran tiempo de operación (más de 30 años). Actualmente existen alrededor de 29 plantas en operación y más de 1220 MWe instalados, lo que corres-ponde a 96% del total de sistemas CSP instalados (Lorente et al., 2011).
Así pues, el potencial de los CSP, principalmente de canal parabólico y torre central, es muy grande y se tie-ne la expectativa de que su capacidad instalada se incremente considerablemente en los próximos años, debido a la disminución de precios. De acuerdo con IEA-ETSAP e IRENA (2013), se espera que para 2015, derivado de los desarrollos tecnológicos y la economía de escala, el costo de generación con este tipo de tecnologías disminuya 15%, y para el 2020 podría disminuir entre 30 y 50%. Este mismo reporte sugiere que el éxito de los sistemas CSP se verá sumamente apoyado con la hibridación de estos con plantas térmicas que emplean combustibles fósiles, particularmente ciclo combinado con turbina de gas y de vapor (conocidos como ISCC), como el que se presenta en este trabajo. Por su parte, el GEF (2006), en su programa operacional incluye la demostración de la viabilidad de estos sistemas con la instalación de una planta con 31MW solares y 530MW de ciclo combinado en la región de Agua Prieta, Sonora, México.
Con respecto a la hibridación de sistemas con plantas térmicas de combustibles fósiles, se realizaron diversas investigaciones (Montes et al., 2011; Nezammahalleh et al., 2010) sobre la combinación de plantas solares con ciclo combinado (ISCC), las cuales emplean la energía solar como energía suplementaria a la que proviene de la turbina de gas; esto permite un incremento de la potencia generada por la turbina de vapor, y una compensación de la disminución del rendimiento de la turbina de gas bajo ciertas condiciones ambientales (Valdés et al., 2006a). En la mayoría de las configuraciones propuestas el campo de canal parabólico se puede utilizar como el economizador de la caldera de recuperación de calor (Montes, 2008).
Tomando los trabajos mencionados como precedente, el objetivo de este artículo es presentar un análisis de diferentes configuraciones de plantas de canal parabólico con generación directa de vapor (GDV) integradas con ciclo combinado.
Ahora bien, se propone la generación directa de vapor porque se considerará al campo solar como el economizador o el sobrecalentador de la planta, por lo que no habrá flujo bifásico en el absorbedor del concentrador. Además, existen varios trabajos que han estudiado la generación directa de vapor (Lentz y Almanza, 2006; Zarza, 2006), que muestran que esta es factible tanto a baja como a alta presión. También está el de Montes (2008), en el cual se compara el desempeño de una planta solar GDV con aquella que utiliza aceite Therminol VP-1 como fluido de trabajo, concluyendo que la GDV presenta un mayor rendimiento energético y exergético porque no requiere un intercambiador de calor.
En cuanto al análisis térmico de la planta solar, distintos trabajos lo realizan, como el de Bakos et al. (2001), que muestran la variación del rendimiento del colector parabólico como función del fluido de transferencia de calor. Asimismo, el de Tyagi et al. (2007) que está relacionado con el análisis por segunda ley de este tipo de sistemas. Sin embargo, existen pocos trabajos relacionados con el análisis termoeconómico y la optimización de los mismos (Baghernejad y Yaghoubi, 2011). Por esta razón se eligió realizar este análisis, considerando además que en trabajos previos (Valdés et al., 2003; Durán et al., 2013) se ha realizado este tipo de análisis a otros sistemas y se han obtenido resultados más realistas.
En trabajos previos realizados por los autores (Durán, 2004), se desarrolló un modelo de optimización termoeconómica con algoritmos genéticos que se aplicó a plantas de ciclo combinado (CC), por lo que se propone aplicar este mismo modelo a la planta de ciclo combinado, integrada con la planta solar.
DesarrolloDescripción y parámetros de la planta solar integrada con ciclo combinadoDescripción: en este trabajo se realiza el análisis de diferentes configuraciones de plantas de canal parabólico integradas con un ciclo combinado (ISCC, Integrated Systems with Combined Cycle). Estas plantas constan de un ciclo combinado común pero incluyen una planta solar que sustituye a una sección de la caldera, consta de los siguientes elementos:
Turbina de gas: transforma la energía térmica de los gases que ingresan a alta temperatura y alta presión provenientes de la cámara de combustión en trabajo. Los gases de escape de este elemento aún tienen suficiente entalpía que puede ser aprovechada en un ciclo de baja temperatura (ciclo de vapor).
Caldera de recuperación de calor: es el elemento de unión entre el ciclo de gas y de vapor, es el encargado de recuperar la energía térmica de los gases provenientes de la turbina de gas para producir vapor. Consta de tres secciones por cada nivel de presión:
- •
Economizador (EC). En esta sección el agua se calienta hasta unos grados (Approach Point) menos que el punto de saturación a la presión de trabajo de la caldera. Para CRC de dos niveles de presión se suele tener uno o dos economizadores por cada nivel de presión.
- •
Evaporador (EV). En esta sección de la caldera presenta el cambio de fase y se lleva al agua desde líquido saturado hasta vapor saturado a la presión nominal. Se suele tener un evaporador por cada nivel de presión de la CRC.
- •
Sobrecalentador (EV). En esta última sección se lleva el vapor saturado a sobrecalentado. Se suele tener un sobrecalentador por cada nivel de presión de la CRC.
En el presente trabajo el campo solar sustituirá, según el caso, al economizador o al sobrecalentador de un nivel de presión de la caldera.
Planta solar. Se trata de una planta de canal parabólico, con generación directa de vapor (GDV) en el presente trabajo sustituirá alguna sección de la caldera de recuperación de calor. La planta solar está constituida por un número lazos de concentradores solares conectados en paralelo, cada lazo es un conjunto de concentradores solares en serie.
Turbina de vapor. Es el elemento de potencia del ciclo de baja temperatura, recibe al vapor de la CRC y lo transforma en trabajo adicional.
En el presente trabajo se estudian 3 configuraciones con 2 niveles de presión sin recalentamiento y sin almacenamiento térmico en el campo solar:
Parámetros de diseño y consideraciones
Para el análisis que aquí se presenta los parámetros de diseño que se consideran durante la optimización corresponden al ciclo de vapor, particularmente a la CRC y a la planta solar. Los parámetros de diseño más importantes de la caldera son los siguientes:
- •
Presión nominal. Es la presión de operación de cada nivel de presión de la caldera.
- •
Pinch Point (PP). Corresponde a la diferencia de temperaturas entre el agua líquida saturada a la entrada del evaporador y los productos de la combustión a la salida del mismo. Este valor en gran medida determinará el tamaño de la caldera y por lo tanto del campo solar. Si se diseñara una caldera únicamente considerando la optimización térmica, este valor tendería a ser muy pequeño, sin embargo, esto ocasionaría un incremento sumamente importante en el costo de la caldera de recuperación de calor o del campo solar en su caso.
- •
Approach Point (AP). Corresponde a la diferencia de temperaturas entre el agua que sale del economizador y la temperatura de saturación de la misma a la presión nominal. Este valor es un parámetro de seguridad que permitirá evitar evaporación en el economizador y por lo tanto en el campo solar cuando sustituya a esta sección de la caldera.
- •
Diferencia de temperaturas en el sobrecalentador. Corresponde a la diferencia de temperaturas entre el vapor a la salida del sobrecalentador y los gases a la entrada del mismo. Esta diferencia determina el tamaño del sobrecalentador y del campo solar si este está acoplado en esta sección.
Ahora bien, en cuanto a los parámetros de diseño de la turbina de gas y de vapor, estos no se considerarán como variables independientes en el proceso de optimización. Se elige una turbina de gas comercial cuyas características de diseño se presentan en la tabla 1. Por su parte, los parámetros relacionados con la simulación del ciclo de vapor tampoco se incluirán dentro de la optimización, se empleará un valor fijo, el cual se presenta en la tabla 2. En el presente trabajo se describe la optimización del ciclo centrada en los parámetros de diseño de la CRC y el campo solar; esta optimización se realizará aplicando la herramienta de algoritmos genéticos, para lo cual debe considerarse un intervalo de variación de los parámetros de diseño (tabla 3). Se puede observar en esta tabla que se proponen valores límite del pinch point y approach point entre 3 y 15 K. Se propone esto porque es lo sugerido para centrales de ciclo combinado y calderas de recuperación de calor (Horlock, 1991).
Parámetros de diseño del ciclo de gas
| Parámetro | Valor |
|---|---|
| Relación de compresión | 30 |
| Pérdida de carga en la admisión | 20mbar |
| Rendimiento isentrópico del compresor | 0.85 |
| Rendimiento de la cámara de combustión | 0.95 |
| Rendimiento isentrópico de la turbina | 0.91 |
| Temperatura de entrada a la turbina | 1430 K |
| Temperatura de salida | 710 K |
| Flujo másico de aire | 120.2 kg/s |
| Pérdida de carga en el escape | 40 mbar |
| Resultados de la simulación | |
| Potencia nominal de la turbina de gas | 38.8 MWe |
| Rendimiento del ciclo de gas | 40.12% |
| Temperatura de entrada a la CRC | 710.7 K |
| Temperatura de la mezcla | 1222.19K |
Parámetros para la optimización termoeconómica con su intervalo de variación en el algoritmo
| Parámetro de diseño | Intervalo de variación |
|---|---|
| Presión en colector de baja | 3-15 bar |
| Pinch Point de la sección de baja presión | 3-15 K |
| Approach Point de la sección de baja | |
| presión | 3-15 K |
| Diferencia de temperatura en el | |
| sobrecalentador de baja presión | 20-70 K |
| Presión en el colector de alta | 50-100 bar |
| Pinch Point de la sección de alta presión | 3-15 K |
| Approach Point de la sección de alta | |
| presión | 3-15 K |
| Diferencia de temperatura en el sobrecalentador de alta presión | 20-80 K |
La optimización del sistema se realiza considerando una irradiancia directa promedio mensual correspondiente a la zona de Cerro Prieto, en Baja California, México, de acuerdo con datos de GoogleTM Earth, la ubicación geográfica del lugar tiene una latitud de 32.40¿ N y una longitud de 115.24¿ O. Las condiciones de Cerro Prieto se presentan en la tabla 4.
Datos geográficos y ambientales de Cerro Prieto, Baja California
| Parámetro | Valor |
|---|---|
| Latitud | 32,40¿ N |
| Longitud | 115,24¿ O |
| Altitud | 12 msnm |
| Temperatura ambiente máxima normal | |
| promedio anual | 32,2 ¿C |
| Irradiación directa normal diaria promedio | |
| mensual (Almanza y Estrada, 2005). | 23,79 MJ/m2 |
| Irradiancia directa instantánea | 797W/m2 |
El rendimiento térmico (η) se define como la razón de la ganancia térmica útil (Qu) en un periodo de tiempo (dt), entre la energía solar incidente (integral de la irradiancia global, GT) multiplicada por el área de apertura del concentrador (AA) en el mismo intervalo de tiempo.
(1)
El periodo de tiempo típico es de 1 h, pero si se tienen datos de irradiación global (IT) de un día, entonces el periodo de tiempo serán las horas pico del día en promedio para un mes o para un año.
Asimismo se define otra variable que relaciona la ganancia útil de energía instantánea de un concentrador entre la ganancia útil si toda la superficie del absorbedor estuviera a la temperatura del fluido, a esto se le llama factor de remoción de calor del concentrador (FR).
El modelo termodinámico aplicado determina el rendimiento térmico instantáneo (ηi) del colector parabólico en función del factor de remoción de calor (Duffie y Beckman, 2006), en el que se toma en cuenta la irradiancia directa instantánea (Gb)
(2)
donde
UL = coeficiente global de transferencia de calor
C = razón de concentración solar
τ = transmitancia
α = absortancia del material del absorbedor
Ti = temperatura de entrada del fluido al colector
Ta = temperatura ambiente
En esta ecuación el factor de extracción de calor (FR) se determina mediante la ecuación
(3)
Donde el factor de eficiencia del colector (F’) está dado como la razón de la resistencia de transferencia de calor del absorbedor y el aire ambiente (Uo) entre la resistencia de transferencia de calor del fluido de trabajo y el ambiente (UL), esto significa que la mejora en el rendimiento de captación de la energía en el concentrador es función del coeficiente global de transferencia de calor y del factor transferencia de calor por convección del fluido. Por su parte es el flujo másico (Cp) es el calor específico a presión constante y AA es el área de apertura total del concentrador (Tiwari, 2002).
El rendimiento del colector determinado mediante la ecuación (2) para un periodo diario promedio anual, se empleará para obtener el tamaño de la planta solar en el proceso de optimización.
Análisis termodinámico del ciclo combinadoSe aplicó un programa desarrollado en macros de Visual Basic para Excel para simular y optimizar el ciclo combinado integrado con el campo de canal parabólico. Este programa emplea el modelo de costo-flujo de caja, propuesto por Rovira (2004).
En el método aplicado se realiza la simulación del ciclo de gas operando en condiciones de diseño, aplicando el modelo descrito en Muñoz (2004) y Facchinni y Stecco (1999). La caldera de recuperación de calor, el campo solar y el ciclo de vapor se simulan aplicando las correlaciones de la IAPSW (Wagner y Kretzchmar, 1988).Wagner y Kretzschmar, 1988
Modelo económicoLa optimización del sistema se realiza minimizando el costo de generación, el cual corresponde al costo total anual dividido por la energía generada anualmente, es decir
(4)
En la ecuación (3), el costo total anual (CTot) incluye el costo de amortización, el costo del combustible y el de operación y mantenimiento de todo el sistema, considerando al campo solar. Por su parte,
Para llevar a cabo el análisis económico del ciclo combinado se consideraron las funciones de costo de la turbina de gas, turbina de vapor y caldera de recuperación de calor que se describen en Durán (2004).
- •
Costo de la turbina de gas: este costo se determinó del análisis de tendencias del costo de turbinas de gas (Valdés et al., 2006b), en el cual se obtiene una expresión del costo de la turbina de gas en función de la potencia de la misma, expresada en MW.
- •
(5)
- •
Costo de la turbina de vapor: este costo se determinó del análisis de tendencias del costo de turbinas de vapor (Valdés et al., 2006b), en el cual se obtiene una expresión del costo de la turbina de vapor en función de la potencia de la misma, expresada en MW.
(6)
- •
Costo de la CRC: incluye varios componentes como el costo del diseño, de la estructura e instalación y el costo de los materiales de la superficie de intercambio de calor. Sería deseable conocer todos estos costos, sin embargo, se cuenta con muy poca información al respecto. Es necesario desarrollar un modelo que sea capaz de predecir el costo de dicha caldera en función de sus parámetros termodinámicos, sin conocer el diseño geométrico de la misma. En la bibliografía se encuentran varios modelos que desarrollan correlaciones para predecir este costo. Algunos autores (Attala et al., 2001) y (Valero, et al. 1994) definen el costo de cada sección de la caldera en función del producto UA de la forma:
(7)
donde el exponente (x) al que se encuentra elevado el producto UA se incluye debido a la economía de escala, es decir, que cuanto mayor sea el área de la caldera, más barata será la unidad de UA. El valor empleado en este trabajo es el propuesto por Valero et al. (1994) y Tsatsaronis y Pisa (1994) quienes consideran este exponente como 0.8 y proponen el costo de la caldera de la siguiente forma
(8)
Por otra parte, para realizar el análisis económico de la planta solar se tomaron en cuenta los costos establecidos por Montes et al. (2011). Vea la tabla 5.
En cuanto a los parámetros financieros, se tomó una tasa de interés de 10% y un tiempo de operación de 25 años, se consideró este valor porque es el tiempo promedio de operación de una planta de ciclo combinado.
El método de optimización que se aplicó es el de algoritmos genéticos (AG), ya que este se aplicó en trabajos previos relacionados con la optimización de ciclos combinados, donde presentaron resultados bastante aceptables. El algoritmo de optimización se explica mediante el diagrama de la figura 3. Como se indica, para cada generación el AG genera aleatoriamente un número establecido de individuos que incluyen las variables de optimización, con ellas se simula el ciclo y la planta solar y se calcula el costo de generación para las condiciones dadas. Finalmente se selecciona la combinación de variables que proporcione el menor costo de generación. Mayor información de esta metodología y su aplicación en la optimización de sistemas térmicos puede encontrarse en Toffolo y Lazzareto (2002) y Valdés et al. (2003).
Discusión y análisis de resultadosParámetros considerados para la simulación de la planta solar
Con el fin de realizar la evaluación numérica de los resultados, se consideró un captador parabólico comercial modelo Eurotrough ET-100. Las especificaciones técnicas de este tipo de captador se presentan en la tabla 6. Especificaciones más detalladas de este tipo de colectores, así como un análisis de su desempeño se puede encontrar en Geyer et al. (2002).
Parámetros de la planta de canal parabólico
| Parámetros del modulo solar | |
|---|---|
| Modelo del colector solar | ET100 |
| Radio del absorbedor | 0.035 m |
| Área de la superficie reflejante | 545 m2 |
| Apertura | 5.77m |
| Longitud del colector | 99.5 m |
| Reflectividad del espejo | 94% |
| Rendimiento óptico pico | 80% |
| Razón de concentración | 23.23 |
| Longitud de cada modulo | 12 m |
| Flujo másico en el absorbedor | 1.4 kg/s |
Es importante señalar que se consideró un flujo másico dentro del absorbedor de 1.4 kg/s porque es el valor sugerido para evitar grandes pérdidas de carga en el mismo (Montes, 2009).
Resultados de la optimizaciónTomando en cuenta las consideraciones mencionadas en los apartados anteriores, se aplicó el programa de optimización a las configuraciones analizadas. Los resultados de la optimización se obtuvieron considerando el mínimo costo de generación (tabla 7).
. Resultados de la optimización para las configuraciones analizadas
| Parámetros termodinámicos | Config. A CCPP s/canal parabólico optimizada | Config. B ISCC | Config. C ISCC |
|---|---|---|---|
| Presión de baja | 324kPa | 389kPa | 371kPa |
| PP de baja | 3.01K | 14.5K | 9.4k |
| AP de baja | 4.07K | 5.91K | 3.09K |
| DT de baja | 84.8K | 69.24K | 51.55K |
| Presión de alta | 6621 kPa | 9639 kPa | 10168 kPa |
| PP de alta | 3.98K | 12.455K | 14.83K |
| AP de alta | 7.6K | 4.25 K | 7.59K |
| DT de alta | 43.145K | 69.95K | 50.5K |
| Costo/kWh | 0.05298€/kWh | 0.04602€/kWh | 0.052€/kWh |
| Rendimiento | 54.68% | 56.0% | 57.9% |
| Potencia del CC | 52986.82kW | 54286.33kW | 53083.0kW |
| Flujo másico de alta presión | 10.89kg/s | 10.95kg/s | 9.26kg/s |
| Flujo másico de baja presión | 5.32kg/s | 9.20 kg/s | 5.84kg/s |
| Número de bucles totales | 0 | 26 | 9 |
| Número de bucles en paralelo | 0 | 7 | 7 |
Es importante aclarar que se trata de una optimización termoeconómica, si se tratara solo de una optimización térmica podría esperarse que los valores del pinch point y el approach point tendieran al mínimo del intervalo, sin embargo, al involucrar el costo el programa de optimización encuentra el valor óptimo con el que se obtiene el menor costo de generación. Por esta razón los valores de la tabla 7 no corresponden al límite inferior del intervalo propuesto.
Primero se optimizó la planta de ciclo combinado de dos niveles de presión, sin recalentamiento y sin el acoplamiento del campo solar, estos resultados se muestran en la segunda columna de la tabla 7. Se observa que esta configuración presenta un rendimiento de 54.68% y un costo de generación de 0.529€/kWh. Estos valores se encuentran en las tendencias de diseño de estos sistemas (Valdés et al., 2006b).
En la columna 3 se presenta la optimización de la ISCC, considerando al campo solar como parte del economizador de alta presión. Se observa que esta configuración es la que presenta un menor costo de generación, ya que el acoplamiento del campo solar en la sección de alta presión deja más energía disponible de los gases calientes en la CRC en la zona de baja presión de la misma, permitiendo incrementar el flujo másico en esta última sección y aumentando así la potencia en la turbina de baja.
Por su parte, en la configuración donde el campo solar es el sobrecalentador de la sección de baja presión de la CRC (Configuración C, ISCC), se observa que, si bien, el rendimiento es mayor que en los otros dos casos, el costo de generación es mayor que el obtenido cuando el campo solar está acoplado en el economizador de alta presión, aunque menor que el obtenido con la configuración A. Esto se debe a que el acoplamiento con el campo solar sí tiene un efecto en el incremento del flujo másico de vapor del ciclo de baja, pero es mucho menor porque la energía que se transfiere en este sobrecalentador es mucho menor.
En cuanto al número de bucles totales obtenidos en cada configuración, es importante aclarar que la diferencia es grande cuando se acopla el campo solar al economizador de alta, contrario a cuando se acopla al sobrecalentador de baja. Esta diferencia tan marcada radica en que el área de captación obtenida para el sobrecalentador de baja es muy pequeña porque este elemento demanda mucha menos energía que el economizador de alta presión. Mayor información acerca del tamaño de las superficies de intercambio de una CRC se puede encontrar en Durán et al. (2013).
Por otro lado, se realizó un análisis de sensibilidad en función de irradiancia directa, únicamente para la configuración B, que fue la que dio mejores resultados. Se observó que la irradiancia tiene una influencia muy importante sobre el flujo másico de agua en el canal parabólico, que al pasar por las secciones de la caldera se convierte en vapor. Esto se muestra en la gráfica de la figura 4. Esta tendencia se debe a que durante el análisis las temperaturas de entrada y salida de los economizadores permanecen constantes, entonces quien varía en función de la irradiancia es el flujo másico.
El análisis se realizó para las condiciones del mes de mayo, ya que es el que mejores condiciones de irradiancia presenta, se observó más de 6 horas que la irradiancia promedio es superior a la promedio considerada; entonces durante este tiempo se tiene un incremento considerable en el flujo de vapor de alta y por lo tanto en el rendimiento del ciclo. Durante las horas de menor irradiancia la falta de energía solar se compensa con la de la turbina de gas, debido a que el vapor que no genera el campo solar, lo genera la CRC en su sección del economizador de alta presión; esto implica una disminución en la temperatura de los gases que llegan a la sección de baja presión, sin embargo, no tiene un efecto notable en el rendimiento como puede observarse en la figura 5. Se observa también que la variación del flujo másico es muy similar a la de la irradiancia, esto se debe a que, como se explicó anteriormente, las temperaturas de entrada y salida del economizador permanecen constantes y quien varía con respecto a la radiación es solo el flujo másico.
Por otro lado, también se realizó el análisis de la variación del costo de generación (gráfica con línea continua en la figura 5) y el rendimiento (gráfica con línea punteada en la figura 5) en función de la irradiancia directa. Se observó que se presenta una variación importante en estos dos parámetros, principalmente en el rendimiento para valores de irradiancia entre 0.76 y 0.86 kW/m2, intervalo en el cual el rendimiento cambia entre 55.6% y 57.2%. Posteriormente, para valores de irradiancia menores a 0.76kW/m2, la turbina de gas compensa la falta de energía del campo solar y por esta razón la variación del rendimiento disminuye considerablemente. Por otro lado, se observa una pendiente constante en la variación del costo de generación con la irradiación, de hecho se obtuvo una correlación lineal con pendiente negativa que se presenta en la misma figura. De acuerdo con la gráfica y la ecuación obtenida se observa que aún con la disminución en el costo de generación, este permanece menor que el obtenido para la configuración de ciclo combinado sin la integración con el campo solar (tabla 7) aún para valores de irradiancia bajos.
ConclusionesEn este trabajo se aplicó un modelo de optimización termoeconómica con algoritmos genéticos para la optimización de diferentes configuraciones de ISCC. Se observó que la metodología es aplicable y brinda resultados interesantes.
Se aplicó el modelo a dos configuraciones de ISCC y se compararon los resultados con un CC convencional; se observa que la planta optimizada de ciclo combinado integrada con la planta solar (ISCC) presenta mejores resultados que la misma configuración de ciclo combinado optimizada, esto empleando el valor de irradiación directa promedio de la región de Cerro Prieto, B.C.
Asimismo, al considerarse la planta de canal parabólico como el economizador o el sobrecalentador de la CRC y tomando en cuenta un Approach Point mayor que 3¿, no se espera tener evaporación en el absorbedor. Lo cual favorece el empleo de GDV y evita el uso de un aceite como fluido de trabajo.
Se observa que es preferible acoplar la planta de canal parabólico en el economizador de la sección de alta presión de la CRC; la razón es que de esta forma se tiene más energía disponible de los gases en la sección de baja presión de la misma y, por lo tanto, se puede incrementar el flujo másico en esta sección y la potencia del ciclo.
Por otro lado, se observó que la variación de la radiación directa tiene un efecto importante sobre el rendimiento y costo de generación; sin embargo, aún con esta variación este costo resulta menor que el obtenido con la configuración que no está integrada con el canal parabólico. Asimismo, resulta interesante resaltar que la gráfica de variación del costo de generación en función de la radiación, sigue una correlación lineal con pendiente negativa.
Se considera importante analizar otras configuraciones con y sin acoplamiento con la planta solar. Esto permitirá establecer tendencias en el diseño de este tipo de sistemas, además de que los diseñadores y tomadores de decisiones de empresas generadoras de energía tendrán una idea más clara de hacia dónde deben dirigirse los valores de los parámetros de diseño y qué desempeño y costo se espera tener, previo a seleccionar un diseño o adquirir una central llave en mano.
La GDV es una nueva opción que bajará los costos, ya que se eliminan el intercambiador de calor y los aceites sintéticos (fluido de transferencia de calor).
Este artículo se cita:
Citación estilo Chicago
Durán-García, María Dolores, Rafael Almanza-Salgado, Ivan Galileo Martínez-Cienfuegos. Estudio del acoplamiento de ciclo combinado con planta solar de canal parabólico. Ingeniería Investigación y Tecnología, XVI, 02 (2015): 253-263.
Citación estilo ISO 690
Durán-García M.D., Almanza-Salgado R., Martínez-Cienfuegos I.G. Estudio del acoplamiento de ciclo combinado con planta solar de canal parabólico. Ingeniería Investigación y Tecnología, volumen XVI (número 2), abril-junio 2015: 253-263.
Semblanzas de los autores
María Dolores Durán-García. Profesor-investigador de tiempo completo en la Facultad de Ingeniería de la Universidad Autónoma del Estado de México. Es ingeniera mecánica de profesión con maestría en técnicas de energías renovables y doctorado en ingeniería térmica. Ha realizado proyectos relacionados con la optimización termoeconómica de sistemas térmicos y plantas de potencia. Asimismo ha realizado proyectos y escrito artículos referentes a concentración solar, principalmente, plantas de canal parabólico.
Rafael Almanza-Salgado. Investigador titular C en el Instituto de Ingeniaría, UNAM, ha trabajado en proyectos de energía solar por más de 38 años. Tiene nivel 3 en el Sistema Nacional de Investigadores y obtuvo el Premio Universidad Nacional 2006, Innovación Tecnológica y Desarrollo Industrial.
Iván Galileo Martínez-Cienfuegos. Profesor-investigador de tiempo completo en la Facultad de Ingeniería de la Universidad Autónoma del Estado de México. Ingeniero químico de carrera con maestría en ingeniería en energía orientada a procesos y uso eficiente de la energía. Doctor en ingeniería en energía, con el proyecto de investigación de flujo bifásico de agua en sistemas termosolares de canal parabólico para generación de electricidad. Ha realizado otros proyectos en la línea de integración de sistemas termosolares y fotovoltaicos para desalación de agua, así como geoquímica del agua de mar en descargas de plantas desaladoras. Sus líneas de investigación son: generación directa de vapor para sistemas termosolares, materiales para concentración solar y modelos de radiación solar.