INTRODUCCIÓN
En Cuba se genera una elevada cantidad de residuos agroindustriales, los cuales son una fuente de sustrato potencial para la producción de biogás. Entre ellos el residuo de secado del arroz (RS), constituido por paja de arroz, cáscara de arroz y granos de arroz fraccionados, se destaca por su volumen de generación (0,04 t/t de arroz húmedo) y disponibilidad en la industria arrocera, cuyo destino es ser depositado en grandes montones, los cuales se encuentran cercanos al secadero, constituyendo en muchas ocasiones fuente de contaminación ambiental (Contreras, 2013).
En la revisión realizada de la literatura se encontraron pocos artículos que reporten el estudio del RS para la producción de biogás (Contreras et al., 2012; Díaz et al., 2016; López y Heiermann, 2021). Rendimientos de metano similares de 164 y 169 LN CH4/kg Sólidos Volátiles (SV) durante la monodigestión del RS obtuvieron Contreras et al., (2012) y López y Heiermann, (2021) en un sistema discontinuo a régimen mesofílico.
La digestión anaerobia de residuos de la agroindustria arrocera, principalmente la paja de arroz y el residuo del secado de arroz tienen una atención especial debido a su elevado contenido en materia orgánica (Contreras et al., 2012). Sin embargo, su alta relación Carbono/Nitrógeno (C/N) y contenido de lignina, así como su compleja, estable y recalcitrante (difícil de degradar por los microorganismo) estructura lignocelulósica, son consideradas limitantes cuando estos residuos son digeridos como único sustrato (Ngan et al., 2020). De ahí que un mejor balance de nutrientes y alternativas de pretratamiento son estrategias a tener en cuenta para mejorar el proceso de digestión anaerobia para este tipo de residuos.
Existen varios tipos de pretratamientos de materiales lignocelulósicos: físicos (mecánico e irradiación), químicos (álcalis, ácidos, oxidantes y organosolventes), físico-químicos (extrusión, hidrotérmico y explosión de vapor) y biológicos (microorganismos y enzimáticos) (Abraham et al., 2020). Su eficiencia varía según la biomasa y tipo de pretratamiento, por lo que la selección y las condiciones a utilizar deben ser probadas para cada biomasa. Recientemente (López y Heiermann, 2021) obtuvieron un incremento de un 63 % en el rendimiento de metano en comparación a la biomasa sin pretratar.
La codigestión ha sido aplicada con resultados positivos en la mezcla de la paja de arroz y cáscara de arroz con residuos de comida (compuesto por 35 % de raíces vegetables, 20 % arroz, 8 % harina, 6 % carnes, y 31 % de otros) (Zhan-Jiang et al., 2014; Jabeen et al., 2015), excreta porcina (Mei et al., 2016; López-González et al., 2019), excreta vacuna (Wang et al., 2014; Haryanto et al., 2018), hydrilla verticillata (Kainthola et al., 2019) y con gallinaza y excreta porcina (Shen et al., 2018). Recientemente López-González et al., (2019) reportaron el efecto positivo de la codigestión del RS con estiércol porcino (RP), al lograrse un incremento del rendimiento de metano en un 18 %, respecto al valor esperado por la contribución individual en la mezcla.
La codigestión anaerobia proporciona beneficios como estabilización del proceso, pH y capacidad de amortiguación, incorporación de macro y micronutrientes en falta, reduce la concentración de compuestos inhibidores o tóxicos (ácidos grasos volátiles/amoníaco), logrando una mejor sinergia, a través de la mezcla de diversos sustratos, y actividad microbiana, con lo cual se mejora el rendimiento de metano en el proceso (Prajapati et al., 2018).
La mezcla del RS con RP y estiércol vacuno (RV) no ha sido reportada en la literatura, sin embargo esta se presenta como una alternativa atractiva para su implementación en empresas agroindustriales de granos, al tener disponibilidad de ganadería vacuna y porcina cercana a la planta de secado del arroz, lo cual posibilita el uso de la codigestión del RS con estos residuos. Por lo que el objetivo del presente trabajo fue evaluar la cinética de la producción de metano a partir de diferentes mezclas del residuo del secado de arroz con estiércol vacuno y porcino en sistema discontinuo.
MATERIALES Y MÉTODOS
La biomasa utilizada como fuente de inóculo en el ensayo de digestión anaeróbica fue una mezcla de lodo anaerobio procedente de un biodigestor chino con residuo porcino como sustrato y el efluente de un reactor de mezcla con 50 L de capacidad, el cual estaba siendo alimentado con residuo de secadero. El inóculo fue tamizado a 2 mm y mantenido en condiciones anaerobias por una semana para asegurar el agotamiento de la materia orgánica contenida.
El RS fue suministrado por la empresa agroindustrial de granos, perteneciente al municipio La Sierpe, Sancti Spíritus, Cuba. Los RP y RV fueron recolectados en granjas de la cooperativa “Hornos de cal”, Sancti Spíritus, Cuba. Al RS se le aplicó un pretratamiento físico (molinado a 3 mm) para facilitar su homogenización.
Los sólidos totales (ST), sólidos volátiles (SV), cenizas, pH y nitrógeno amoniacal (N-NH4 +), se determinaron de acuerdo a las normas establecidas por los métodos estandarizados (APHA, AWWA, WEF., 1995). La hemicelulosa, celulosa y lignina fueron cuantificadas de acuerdo al método de Van Soest (Van et al., 1991).
La experimentación a escala de laboratorio se realizó en reactores de polietileno de 2 litros de capacidad, a temperatura constante de 37 ± 2 ºC, controlada por un baño de agua con termostato. Los reactores fueron cerrados con tapón de goma y conectados, mediante una manguera, a colectores para la salida del biogás. En la Figura 1 se muestra la configuración del set experimental utilizado.
El líquido desplazado en los colectores consistió en una solución saturada de NaOH al 3%, para facilitar la disolución de CO2 en agua, contabilizándose el metano obtenido mediante el desplazamiento del líquido (VDI-4630, 2006). A cada reactor se le adicionó 900 g de inóculo y una cantidad determinada de sustrato para mantener una relación inóculo/sustrato de 1 (g SV/ g SV), recomendada en el caso de digestión en discontinuo a altas cargas de sólidos (VDI-4630, 2006). Las características físico-químicas del inóculo fueron: 3,41 ST%MF, 1,58 SV%MF y pH: 8,6. Los ensayos se realizaron en triplicado. Todos los reactores fueron agitados manualmente una vez al día. El ensayo anaerobio se realizó por 25 días, cuando la producción de metano fue inferior en un 1 % a la acumulada (VDI-4630, 2006).
El rendimiento de metano PCH4 (LN CH4/kg SV) se calculó como la cantidad de metano acumulada durante el tiempo de digestión (VCH4) por la masa de sustrato alimentada en término de sólidos volátiles (MSV) (Ecuación 1), y se corrigió a condiciones estándar de temperatura (273 K) y presión (101,3 kPa).
Se empleó un diseño de vértice extremo, de acuerdo a la disponibilidad de residuos en la Empresa “Sur del Jíbaro”. Como restricciones para el diseño, se utilizaron valores de 86-92%, 6-12% y 2-4% (en base al porciento de SV de cada componente en la mezcla) para el RS, RV y RP, respectivamente. Estas restricciones constituyen el intervalo de experimentación para cada componente de la mezcla. La cantidad a adicionar de cada sustrato en la mezcla se calculó por la Ecuación 2.
Donde PSV es el porciento de SV del sustrato en la mezcla definidos por la Matriz de diseño experimental (Tabla 1), y SVsustrato es el % de SV contenidos en el sustrato, determinado en los análisis químicos.
Para el procesamiento estadístico de los parámetros evaluados se empleó el paquete estadístico Statgraphics Centurión XV. Se empleó el análisis de varianza de clasificación simple y para comparar las medias se utilizó la prueba de intervalos múltiples de Duncan. La significancia de los modelos se evaluó con un nivel de significación del 5% (p-valor ≤ 0,05).
El comportamiento cinético fue descrito a través de tres modelos de reacciones de primer orden (Ecuaciones 3-5), y el modelo modificado de Hill (Ecuación 6). Estos modelos son los más frecuentemente usados para analizar el comportamiento de la producción acumulada de metano P (LN CH4/kg SV) en un tiempo dado t (d), y permitieron un ajuste adecuado para sustratos similares a los estudiados en la presente investigación (López-González et al., 2019; López y Heiermann, 2021).
Modelo de Roediger (Ecuación 3):
Modelo de Gompertz modificado (Ecuación 4):
Modelo de Función Transferencial (Ecuación 5):
Modelo modificado de Hill (Ecuación 6):
A partir del ajuste de los datos experimentales a los modelos, los parámetros cinéticos obtenidos son: el valor del máximo rendimiento de metano Pm (LN CH4/kg SV), la constante cinética de velocidad aparente k (d-1), la duración del desfase inicial en la producción de metano λ (d), la máxima producción específica de metano Rm (LN CH4/kg SV d), el tiempo al cual la mitad del rendimiento de metano es alcanzado c (d) y b un parámetro para ajuste del modelo.
RESULTADOS Y DISCUSIÓN
3.1 Caracterización de los sustratos
La Tabla 2 presenta los parámetros analizados al RS, RP y RV. El RS tuvo un contenido de materia orgánica, expresado en términos de SV, de 77 % Materia fresca (MF), superior a valores reportados previamente (Contreras et al., 2012; López y Heiermann, 2021), y al obtenido para el RP y RV con 7% MF y 10% MF, respectivamente.
pH | - | 6,07 | 7,57 | 7,85 |
ST | % MF | 83,12 | 12,36 | 9,27 |
SV | % ST | 92,83 | 84,47 | 77,99 |
Cenizas | % ST | 7,17 | 15,53 | 22,31 |
N-NH4 + | g/kg | 0,071 | 0,255 | 1,714 |
Hemicelulosa | % ST | 18,31 | n.a | n.a |
Celulosa | % ST | 19,08 | n.a | n.a |
Lignina | % ST | 6,96 | n.a | n.a |
n.a: No analizado
Sin embargo, el alto contenido de fibra detergente neutra del RS (44,35 %ST), dado por los componentes lignocelulósicos (celulosa, hemicelulosa y lignina) que forman su estructura, le aporta resistencia a este material, y puede afectar su degradación biológica. El contenido de lignina y de fibra detergente ácida (celulosa+lignina) es inferior al reportado por Contreras et al., (2012), lo cual se debe a la variabilidad en la composición de este residuo entre cosecha y durante el procesamiento del arroz.
3.2 Comportamiento cinético de la Digestión anaerobia
En la Tabla 3 se detallan los parámetros cinéticos y sus estadísticos a partir del ajuste de los modelos: de Roediger (MR), de Hill modificado (MH), de Gompertz modificado (MG) y de función Transferencial (MT). Como se puede observar todos los modelos se adecuaron a los datos experimentales del rendimiento de metano acumulativo con un valor de R2aj superior al 98,54 %. Por tanto, los modelos obtenidos explican entre un 98,54 y 99,70 % la variabilidad de los datos.
El error estándar de estimación estuvo entre 4,92 y 10,76, con valores superiores para el modelo de Gompertz modificado.
Para todas las mezclas los valores de Pm para MR, MH y MT son los que presentan un menor error y son los modelos que mejor predicen el comportamiento de estos residuos.
Raj 2: Coeficiente de determinación ajustado %, ES: Error estándar del estimado, P: rendimiento de metano experimental (LN CH4/kg SV), Pm: rendimiento máximo de metano (LN CH4/kg SV), k (d-1), Rm: velocidad máxima de producción de metano (LN CH4/kg SV d), Fase de retardo λ (d).
En la Tabla 3 se puede apreciar que existe una diferencia estadísticamente significativa entre las mezclas (RS92:RV6:RP2-RS86:RV10:RP4), (RS92:RV6:RP2-RS86:RV12:RP2), (RS90:RV6:RP4-RS86:RV10:RP4) y una semejanza entre las mezclas (RS92:RV6:RP2-RS90:RV6:RP4), (RS86:RV12:RP2-RS90:RV6:RP4) y (RS86:RV12:RP2-RS86:RV10:RP4) para el modelo de Hill, mientras que para los tres modelos restantes no hubo diferencias significativas. El método empleado actualmente para discriminar entre las medias es el procedimiento de diferencia mínima significativa (LSD) de Fisher. Con este método hay un riesgo del 5,0% al decir que cada par de medias es significativamente diferente, cuando la diferencia real es igual a cero.
El modelo de Hill tuvo un mejor ajuste a las curvas experimentales obtenidas para las mezclas RS92:RV6:RP2, RS90:RV6:RP4 y RS86:RV10:RP4 por la tendencia de estas al crecimiento para t>20 días, ya que al tener mayor proporción de residuo de secadero la etapa de hidrólisis es más lenta. Sin embargo, para la mezcla RS86:RV12:RP2 se puede observar una tendencia de la curva a alcanzar la meseta más rápido, ya que esta es la de menor proporción de residuo de secado y mayor proporción de residuo vacuno lo que debió favorecer la hidrólisis del RS. El RV contiene bacterias celulolíticas que rompen celulosa, el cual es un material difícil de degradar.
En cuanto al Rm, estimado por dos modelos, el MT y MG, los valores para el MT fueron superiores hasta en un 45 % respecto al MG. El MT obtuvo un valor de R2aj superior por lo que se puede estimar mejor el valor de Rm para el RS y las mezclas. El mejor valor (60,31 LN CH4/kg SV d) se obtuvo para la mezcla RS92:RV6:RP2 lo cual evidencia el efecto sinérgico de la mezcla atribuido a un mejor balance de nutrientes y a una mayor versatilidad de la comunidad microbiana.
La fase de retardo fue despreciable para el RS y las mezclas, con un valor menor a un día. Los valores estimados de λ por MT y MG difieren y resulta difícil llegar a una comparación clara entre su comportamiento, de acuerdo a las mezclas estudiadas.
3.3 Rendimiento de metano para las mezclas estudiadas
En la Figura 2 se representa el rendimiento de metano obtenido experimentalmente y las curvas ajustadas por el modelo modificado de Hill para las cuatro mezclas estudiadas.
Como se puede observar el mayor Pm se obtuvo para la mezcla RS92:RV6:RP2, seguido de la mezcla RS90:RV6:RP4, RS86:RV12:RP2 y RS86:RV10:RP4, con valores de 351,19 LN CH4/kg SV, 334,19 LN CH4/kg SV, 299,88 LN CH4/kg SV y 269,52 LN CH4/kg SV, respectivamente.
Los valores de Pm obtenidos por los modelos fueron superiores al P, con la excepción del modelo modificado de Gomperzt, donde la tendencia de la curva es más aplanada al final.
El Pm fue significativamente superior para las mezclas RS92:RV6:RP2 y RS90:RV6:RP4, con una mayor proporción de RS (Tabla 3). El incremento de Pm para la mezcla RS92:RV6:RP2, de acuerdo a los valores estimados para los cuatro modelos, estuvo entre 10-17% y 17-30% respecto a las mezclas RS86:RV12:RP2 y RS86:RV10:RP4, respectivamente. Este incremento para la mezcla RS90:RV6:RP4 fue entre 11-14% y 16-24%, en comparación con las mezclas RS86:RV12:RP2 y RS86:RV10:RP4, respectivamente. La adición de RV a la digestión anaerobia del RS permitió la mejora de la dinámica microbiana al introducir bacterias hidrolíticas del tipo celulítico lo cual favoreció la fase de hidrólisis. En la literatura, no se encontraron investigaciones en las que se abordara el estudio del RS en codigestión con RV y RP. Solo se encontró un estudio sobre la codigestión del RS y RP (López-González et al., 2019). En este se determinó el rendimiento de metano de los sustratos en monodigestión y para las mezclas RS98:RP2 y RS95:RP2 (en base al % de sólidos volátiles en la mezcla). El máximo rendimiento de metano (401,20 LN CH4/kg SV) y efecto sinérgico de la codigestión (18% respecto al valor esperado), se obtuvo para la mezcla RS95:RP2, superior en 20,61% respecto a la monodigestión del RS.
Reportes de la codigestión de paja de arroz y cáscara de arroz con otros residuos, como residuos porcinos, vacunos, y de cocina, mostraron un efecto positivo en el rendimiento de metano (Kainthola et al., 2019; Kainthola et al., 2020). En la codigestión de la paja de arroz con residuos de alimentos para diferentes relaciones C/N de 25, 30, y 35 se obtuvo similar rendimiento de metano de 294,17 LN CH4/kg SV, con un incremento de un 71% en comparación a la monodigestión de la paja de arroz (Kainthola et al., 2020). Un resultado diferente obtuvo el mismo autor durante la codigestión de la paja de arroz y la hydrilla verticillata (sustrato rico en nitrógeno). En este caso diferentes relaciones C/N revelaron diferencias en la estabilidad del proceso y el rendimiento de metano. El valor máximo se obtuvo para la mezcla de paja de arroz y la hydrilla verticillata con una relación de C/N de 25, con un incremento del 40% respecto a la paja de arroz sola (Kainthola et al., 2019). La DA de la cáscara de arroz en mono y codigestión ha sido reportada por varios autores (Haider et al., 2015; Syafrudin et al., 2016; Sunwanee et al., 2019). En otro trabajo se estudiaron cuatro relaciones de mezcla de residuos de alimentos con cáscara de arroz. Como resultado se obtuvo que la producción de biogás decreció a medida que lo hizo la proporción de residuos de alimentos en la mezcla, demostrando la baja biodegradabilidad de la cáscara de arroz (Haider et al., 2015).
La producción de biogás a partir de aguas residuales del arroz Thai con cáscara de arroz y tres tipos de excretas animales (gallinaza, excreta vacuna y excreta de codorniz), con y sin ceniza fue investigada por Sunwanee et al., (2019). Cada digestor se alimentó con 10 g de cáscara de arroz, 10 g de excreta animal y 200 mL de aguas residuales. El máximo rendimiento de metano (311,2 LN CH4/kg SV eliminados) fue para la mezcla de aguas residuales con gallinaza y ceniza. Este resultado demuestra la necesidad de macronutrientes (N específicamente) y micronutrientes que se requieren para una eficiente DA de la cáscara de arroz, el cual es el componente principal del RS.
Los efectos sinérgicos son atribuidos principalmente a un mejor balance en la relación C/N, macro y micronutrientes, materia orgánica biodegradable y materia seca (Mata-Alvarez et al., 2011). Por tanto, durante el diseño de plantas de biogás para tratar los residuos que se generan en empresas agro-industriales se le debe dar una importancia significativa a la incorporación de co-sustratos en la alimentación, que permitan mejorar el balance de macro y micronutrientes en el proceso de DA. La decisión de la tecnología a utilizar debe basarse en reactores con agitación ya sea del tipo horizontal (digestión seca) o vertical (digestión húmeda), lo cual debe ser previamente evaluado en experimentos a régimen semicontinuo. La reconversión tecnológica de mezcladores existentes en la empresa puede ser opción para la reducción de costos.
CONCLUSIONES
El rendimiento de metano acumulado obtenido experimentalmente ajustó bien a los cuatro modelos utilizados con un valor de R2aj superior al 98,54%.
El mejor rendimiento de metano en sistema discontinuo se obtuvo para la mezcla RS92:RV6:RP2 (351,19 LN CH4/kg SV), pero sin diferencias significativas con la mezcla RS90:RV6:RP4 (LN CH4/kg SV).
Investigaciones futuras deben dirigirse a la evaluación en sistema semicontinuo de la digestión anaerobia del residuo del secado del arroz utilizando reactores con agitación, codigestión con estiércol vacuno y porcino y la incorporación de una fase de pretratamiento para incrementar su biodisponibilidad.