Introducción
El desarrollo de la sociedad humana se ha relacionado con la capacidad de producir y aprovechar la energía. Los primeros humanos, usaron la energía para sostener las necesidades más básicas, y tan pronto se desarrolló la tecnología para producción de energía, el nivel de vida aumentó.
El consumo de energía para satisfacer las demanda de origen social, el transporte, la agricultura y la industria se requiere la generación de cantidades crecientes de energía, lo que depende del uso eficiente de todas las fuentes disponibles, incluida la gasificación de la biomasa [1].
La gasificación es la conversión de un caudal sólido en un combustible gaseoso y/o líquido que puede ser quemado para liberar energía o ser utilizado para la producción de sustancias químicas valiosas y constituye uno de los recursos que permite sustituir una parte importante de los combustibles fósiles; requiere de un agente gasificante (aire, oxígeno o vapor) para reordenar la estructura molecular de la biomasa y convertirla en un combustible útil con una relación mayor hidrógeno--carbón (H/C) [2]. El gas de salida depende de la cantidad del carbono que se produce durante la pirolisis y del agente gasificante utilizado en el proceso. La calidad del gas producido en la instalación depende no solo de las características del flujo de gas, del agente gasificante utilizado, del tipo de biomasa experimental y del diseño del reactor, sino también de los parámetros de operación [3].
La relación de equivalencia tiene un fuerte efecto en el proceso de gasificación.numerosos autores [4, 5]. refieren esta relación como factor preponderante en el aumento de temperatura en el proceso de gasificación y su relación inversamente proporcional al contenido de CO y H2,.Por lo tanto, los estudios de esta relación se reducen a la concentración de CO y H2[6].
En el estudio bibliográfico realizado sobre el modelado en instalaciones tipo downdraft, relacionados de alguna manera con la operación del proceso se destacan:
Modelos de equilibrio termodinámico, utilizado para predecir la composición del gas de síntesis suponiendo una mezcla de reactivos en un tiempo infinito, los que clasifican en modelos estequiométricos y no estequiométricos, los primeros se basan en constantes de equilibrio [7], para una identificación final de los componentes del gas a través de las diferentes reacciones existentes; los segundos minimizan la energía libre de Gibbs [8] y no toman en cuenta la presencia de las diferentes reacciones.
Los modeloscinéticos. se enfocan en realizar un análisis del gas de síntesis producto del proceso y predecir su rendimiento en función de su composición en un volumen de control en un tiempo dado [9], con este modelo se predicen los perfiles tanto de composición del gas de síntesis como los de temperatura dentro del reactor en condiciones de operación dadas.
Modelos basados en la dinámica de fluidos computacional (CFD), los que permiten realizar una simulación rápida y eficaz del flujo de fluidos y los procesos de transferencia de calor y masa dentro del reactor, con lo que se puede calcular las fuerzas que actúan en el fluido y ayuda a comprender el impacto del gas en el rendimiento de la instalación de gasificación y se basan en la solución de ecuaciones simultaneas obtenidas de los balances de energía, masa y momento en cualquier zona del reactor downdraft, para predecir y formular perfiles de temperatura. Se han utilizado para caracterizar el comportamiento de varios tipos de gasificadores.
Modelos basados en redes neuronales artificiales a partir de datos experimentales: Con el uso de esta herramienta se han desarrollado algunos trabajos en el campo de la gasificación de la biomasa y el carbón. Es así que en el año 2010 María Puig et al. [10] desarrolla un trabajo referente a la revisión y análisis de los modelos de gasificación de biomasa, los autores analizan varios modelos del proceso de gasificación, los que se basan en el modelo de equilibrio termodinámico, cinético y de redes neuronales para predecir el comportamiento del reactor enfocado a la producción del gas y sus componentes, hacen una descripción del proceso de creación de la red neuronal multicapa (MFNN) que sirve como un estimador de proceso no medido en parámetros difíciles de modelar, los autores se refieren a este diseño como función de aproximación universal, con la capacidad de aproximar cualquier función, sin conocer la estructura de la función que ha ser aproximada [11-12].
El estudio del estado del arte muestra muy pocas publicaciones dedicadas al modelado por medio de técnicas de regresión a partir de un diseño experimental adecuado a la tarea que se quiere resolver, entre ellos se encuentra Vladimir Kirsanovs et al (2017) [13-14]. En este trabajo el autor realiza una investigación experimental del proceso de gasificación en una planta real de 400 kW con el fin de caracterizar indicadores de eficiencia de la operación. Se hallan cuatro modelos para definir el calórico del gas de síntesis, la capacidad del gasificador y la eficiencia en frio y en caliente del gas. No obstante, carece de un análisis sistémico previo que permita definir el conjunto de relaciones que se requiere modelar con el fin de operar racionalmente el proceso.
El objetivo del presente trabajo consiste en la elaboración de modelos de los indicadores de desempeño requeridos para la operación de instalaciones de gasificación tipo downdraft a partir de los datos de un plan experimental correspondientes al análisis sistémico para la operación de este tipo de instalaciones utilizando el análisis de regresión no lineal.
Materiales y métodos
Instalación experimental:
Está constituida por un gasificador tipo downdraft con una capacidad de 10 kW conectado a un motor de combustión interna de 0.8 kW de potencia nominal.
Como resultado del análisis sistémico realizado en correspondencia con la metodología de análisis y síntesis de sistema de ingeniería expuesta en [15], se obtuvo la siguiente composición de variables para la tarea de operación de la instalación de gasificación downdraft.
Variables de coordinación:
Generación deseable de flujo horario de gas que se requiere de la instalación, con el fin de satisfacer las necesidades del consumidor determinada por el nivel superior.
Disponibilidad máxima de los diferentes tipos de biomasa para ser consumidos en él.
Indicadores de eficiencia:
Eficiencia de la instalación Ef, Temperatura de la zona de pirolisisTzp, poder calórico del gas PCG, caudal másico de gas generado MG.
Variables de decisión:
Humedad Hu
Cantidad de biomasa Cbio
Caudal de aire de combustión Cau.
Modelo matemático conceptual para la operación del proceso: Se determina por la el modelo de optimización dados por las ecuaciones (1,2,3,4,5)
Minimizar:
Deben, por tanto, ser identificadas las ecuaciones que permitan cuantificar las relaciones (6,7,8,9)
Una vez identificadas las relaciones, el modelo (1) - (5) pasa a ser un modelo de programación estocástica, dado el origen experimental de las relaciones componentes, relativamente sencillo.El plan experimental con vistas a la identificación de los modelos, se concibió con una estructura 3N (alto, medio y bajo, por variables) con tres réplicas, con el fin de tomar en consideración el carácter necesariamente no lineal de las relaciones requeridas y que queden reflejados los factores casuales, incluyendo los errores de medición, para modelar los indicadores de desempeño como función de las variables de operación, con ayuda de modelos de regresión no lineales. Para la recopilación automatizada de la información experimental, se diseñó e implementó en la instalación una tarjeta de adquisición de datos myRIO y la interface hombre-máquina con un programa desarrollado en LabView.
Para estas variables se seleccionaron los niveles 12, 16, 20 %; 1, 3, 5 Kg; 0.1666, 0.2098, 0.2531 g/s; respectivamente. El intervalo de Hu se determina por el valor mínimo alcanzable mediante secado al sol, único admisible para las comunidades agrícolas y el mayor se selecciona por el mayor admisible para la ignición adecuada de la instalación. El valor mayor del intervalo de Cbio se selecciona por la capacidad de la instalación experimental y el mínimo por el requerido para que se produzca la generación de gas, obtenido como resultado de experimentos auxiliares realizados. El intervalo de valores de Cau se determina por los valores mínimos requeridos para que se produzca una combustión incompleta con la generación de gas con el máximo poder calórico posible, para el intervalo de valores de Cau adoptado. La materia prima seleccionada obedece a la necesidad de experimentar con materiales de bajo, medio y alto poder calórico.
La instrumentación utilizada consta de una alarma digital de detección de CO marca Kidde con una intervalo de medición de 30-999 PPM con error ±30 % de la magnitud medida en condiciones normales; un medidor de humedad MD914, con un intervalo de medición de 2 a 60%, una resolución de 0.5%, un sensor de nivel de combustible controlado mediante una paleta para censar el nivel de biomasa, 8 termopares tipo chromel - alumel con cable de compensación tipo K, con intervalo de medición de 0 a 1250°C de temperatura, con un error máximo de ± 2.2% °C, ubicados de la siguiente forma: cuatro para las zonas de reducción, pirolisis, combustión y secado, una en la descarga de cenizas, y tres para censar las temperaturas del gas a la entrada del ciclón, a la entrada del intercambiador de calor y a la salida del intercambiador de calor.
Se instalaron, además, dos manómetros de agua tipo U con intervalo de medición de 200 mm de H2O, con error de medición de ± 1 mm de H2O uno de ellos toma la diferencia de presión existente en el reactor y el otro en la placa orificio para medición de caudal másico de aire de ingreso en las diferentes aperturas de válvula. La disposición de toda la instrumentación se ilustra en la figura 1.
Para elaborar las ecuaciones de regresión que se ajusten adecuadamente a los resultados experimentales se utilizan las técnicas de linealización de modelos en el esquema que se propone en el libro Sistemas de Ingeniería[16].
En particular, se linealiza el modelo monomial dado por la ecuación (10).
Con ayuda de la búsqueda de las potencias de mejor ajuste del modelo dado por la ecuación (11):
Se generan términos posibles de interacción x αi x αk …, con 1, 2, 3 o 4 potencias, permaneciendo en el modelo solo aquellos que resultan significativas, lográndose finalmente modelos signomialescon relación significativa de las variables de la salidacon los términos utilizados.
Resultados
Siguiendo la metodología de selección de modelos expuesta anteriormente se hallaron los siguientes modelos,ayudados por la metodología de elaboración de modelos expuestas en materiales y métodos. Como resultado fueron encontrados los siguientes modelos para la estimación de los indicadores de desempeño para las biomasas estudiadas en el artículo.
Biomasa eichhorniacrassipes (Lechugin)
Para la estimación de la eficiencia (Ef) se halló el modelo dado por la ecuación (12):
Como se puede observar en la figura 2, el R2indica que el modelo matemático de la eficiencia (Ef) es de 61.38% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 59.36%.El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 1.98.
Para la estimación del caudal másico de gas generado (MG) se halló el modelo dado por la ecuación (13):
Como se puede observar en la figura 3, el R2 indica que el modelo matemático del caudal másico del gas generado (MG) es de 92.01% de variabilidad.
La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 91.91%. .El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 0.18.
Para la estimación de la temperatura de la zona de pirolisis (Tzp) se halló el modelo dado por la ecuación (14):
Como se puede observar en la figura 4, el R2 indica que el modelo matemático de la temperatura de la zona de pirolisis (Tzp) es de 83.40% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 82.53%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 4.34.
Para la estimación del poder calórico del gas (PCG):se halló el modelo dado por la ecuación (15):
Como se puede observar en la figura 5, el R2 indica que el modelo matemático del poder calórico del gas (PCG) es de 86.20% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 85.48%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 0.71.
Biomasa Cascarilla
Para la estimación de la eficiencia (Ef) se halló el modelo dado por la ecuación (16):
Como se puede observar en la figura 6, el R2 indica que el modelo matemático del poder calórico del gas (Ef) es de 60.05% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 57.95%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 2.03.
Para la estimación del caudal másico de gas generado (MG) se halló el modelo dado por la ecuación (17):
Como se puede observar en la figura 7, el R2 indica que el modelo matemático del caudal másico de gas generado (MG) es de 93.02% de variabilidad. La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 92.75%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 0.18.
Para la estimación de la temperatura de la zona de pirolisis (Tzp) se halló el modelo dado por la ecuación (18):
Como se puede observar en la figura 8, el R2 indica que el modelo matemático de la temperatura de la zona de pirolisis (Tzp) es de 84.43% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 83.39%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 4.24.
Para la estimación del poder calórico del gas (PCG) se halló el modelo dado por la ecuación (19):
Como se puede observar en la figura 9, el R2 indica que el modelo matemático del poder calórico del gas (PCG) es de 84.69% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 83.66%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 0.77.
Biomasa Leña
Para la estimación de la eficiencia (Ef) se halló el modelo dado por la ecuación (20):
Como se puede observar en la figura 10, el R2 indica que el modelo matemático del poder calórico del gas (Ef) es de 60.87% de variabilidad. La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 58.81%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 1.99.
Para la estimación del caudal másico de gas generado (MG) se halló el modelo dado por la ecuación (21):
Como se puede observar en la figura 11, el R2 indica que el modelo matemático del caudal másico de gas generado (MG) es de 89.94% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 89.68%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 0.21.
Para la estimación de la temperatura de la zona de pirolisis (Tzp) se halló el modelo dado por la ecuación (22):
Como se puede observar en la figura 12, el R2 indica que el modelo matemático de la temperatura de la zona de pirolisis (Tzp) es de 80.29% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 79.24%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 4.74.
Para la estimación del poder calórico del gas (PCG) se halló el modelo dado por la ecuación (23):
Como se puede observar en la figura 13, el R2 indica que el modelo matemático del poder calórico del gas (PCG) es de 84.80% de variabilidad, La R2 ajustada que es más adecuada para comparar modelos con diferentes números de variables independientes como son el caudal de aire de combustión (Cau), cantidad de biomasa (Cbio) y Humedad (Hu), es de 83.78%. El error estándar de la estimación muestra que la desviación estándar de los residuos es del valor de 0.76.
Conclusiones
El trabajo realizado permite extraer las siguientes conclusiones principales:
En el caso concreto de la operación de las instalaciones de gasificación tipo downdraft se requieren identificar la eficiencia de la instalación (Ef), temperatura de la zona de pirolisis (Tzp), poder calórico del syngas (PCG), caudal másico del gas (MG) como función de la humedad (Hu), cantidad de biomasa (Cbio) y caudal de aire de combustión (Cau).Las variables de salida se corresponden con los indicadores de eficiencia y las funciones restringidas y las independientes por las variables de decisión y el dato de entrada Hu. derivados todos del análisis externo realizado.
En ausencia de modelos que satisfagan las necesidades del análisis sistémico se requiere realizar planes experimentales con suficiente número de ensayos para lograr la identificación suficientemente precisa mediante los métodos de procesamiento de efectividad comprobados, entre los que se encuentran los modelos de regresión no lineal de los parámetros energéticos de desempeño de un gasificador tipo downdraft.
La aplicación de los modelos de regresión no lineal de los parámetros energéticos al caso concreto del desempeño de gasificadores tipo downdraft para las biomasas lechugin, cascarilla de arroz y leña se obtienen modelos individuales por cada indicador de desempeño de la operación como sonla eficiencia de la instalación (Ef), temperatura de la zona de pirolisis (Tzp), poder calórico del syngas (PCG), caudal másico (MG).
Dada la poca cantidad de información primaria que se puede obtener de los experimentos planificados en una instalación experimental, los modelos obtenidos se deben considerar como iniciales para implementar los procedimientos de operación, los definitivos se deben de obtener como resultado de la operación permanente del proceso, mediante su renovación paulatina en un proceso de adaptación permanente.