Implemetacion Del Modelo Adm1 en Matlab y Simulink

6.
Implementaci on del modelo ADM1 en Matlab/Simulink

Uno de los objetivos principales de este estudio es realizar la comparaci on de los resultados experimentales obtenidos mediante el desarrollo de los ensayos de biodegradabilidad anaer obica sobre la mezcla de residuos c tricos y esti ercol bovino, respecto a los resultados obtenidos mediante la implementaci on del modelo ADM-1 para la simulaci on del proceso de digesti on anaer obica realizado bajo condiciones similares. La importancia de la implementaci on de un modelo especico del proceso de digesti on anaer obica radica en que los modelos matem aticos pueden simular las condiciones reales de operaci on de las instalaciones a escala de laboratorio o escala de planta piloto con una reducci on notable del tiempo y la energ a requerida para la obtenci on de estos resultados [44]. De esta manera mediante la aplicaci on del modelo ADM-1 es posible realizar la simulaci on del desempe no del proceso de digesti on anaer obica de residuos org anicos biodegradables, a trav es de la generaci on de las curvas que describen la evoluci on temporal de la concentraci on de los metabolitos presentes en el sustrato durante el desarrollo del proceso. De igual forma mediante la implementaci on del modelo es posible establecer la producci on acumulada de biog as generado, y la composici on de esta mezcla de gases en t erminos de las fracciones de metano y di oxido de carbono; lo cual es importante desde el punto de vista de la valorizaci on energ etica del biog as generado en el proceso. Por otra parte la rentabilidad de una aplicaci on a gran escala se basa en la selecci on apropiada del tama no del equipo y de los par ametros de operaci on, estas caracter sticas pueden ser establecidas mediante el desarrollo de un modelo matem atico generalizado sobre el proceso global de digesti on anaer obica. No obstante algunos estudios est aticos desarrollados en reactores a escala de laboratorio no presentan las mismas caracter sticas que pueden presentarse en un sistema de ujo continuo como en las instalaciones a escala piloto o escala industrial [69]. El modelamiento del proceso de digesti on anaer obica permite evaluar el potencial de biometanizaci on o biodegradabilidad anaer obica de los sustratos org anicos biodegradables, bajo las consideraciones de inhibici on del proceso debido a las interacciones biol ogicas y sicoqu micas de los constituyentes primarios del sustrato [98].
6.1.
Condiciones generales
El modelo ADM-1 emplea una gran cantidad de par ametros, constantes y coecientes en las ecuaciones que describen los procesos bioqu micos y sicoqu micos. Adem as requiere que se establezcan las concentraciones iniciales de las especies solubles e insolubles, las cuales representan en conjunto tanto las concentraciones de los compuestos presentes en el sustrato as como las concentraciones de las poblaciones microbianas globalizadas. Por lo tanto se
72
6 Implementaci on del modelo ADM1 en Matlab/Simulink
procedi o de acuerdo a la metodolog a descrita en algunas investigaciones [17, 20, 39], en las cuales se desarroll o la implementaci on del modelo mediante la utilizaci on del conjunto de valores establecidos para estas variables y par ametros, denominado caso base de ADM-1, conforme a como se presentan en el reporte t ecnico original emitido por por la IWA-ADM1 [14]. Algunas investigaciones [20, 131, 44, 113, 97], se nalan la factibilidad de la utilizaci on de los valores reportados como el caso base del modelo ADM-1 para este conjunto de constantes y par ametros, a pesar de que estos valores son el resultado de la mediciones exhaustivas durante el tratamiento anaer obico sobre fangos cloacales mediante un sistema UASB (Manto de lodo anaer obico con ujo ascendente, por sus siglas en ingl es). Adicionalmente existen estudios en los cuales se argumenta la limitada variabilidad de los valores de estos par ametros respecto a la composici on del sustrato especico sobre el cual se realiza el proceso de digesti on anaer obica [115, 28, 111, 13].
La implementaci on desarrollada en este trabajo fue realizada mediante la plataforma computacional Matlab/Simulink la cual fue escogida para realizar la simulaci on con base en estudios previos en los cuales se ha evaluado el modelo ADM-1 bajo diferentes plataformas entre las cuales se incluyen WEST y Aquasim [17, 147]; de las cuales se concluy o que Matlab/Simulink es la plataforma m as exible para la simulaci on de condiciones din amicas durante el proceso de digesti on anaer obica tales como aumento de la raz on de carga org anica e inclusive el cambio de la composici on del sustrato que esta siendo digerido. Otras investigaciones como la conducida por Parker [105], consider o la implementaci on del modelo ADM-1 sobre una variedad de sistemas para la digesti on anaer obica en las cuales los resultados mostraron que en la mayor a de casos el modelo es capaz de reproducir las tendencias de los resultados experimentales. As mismo en la investigaci on desarrollada por Feng et al. [39], se encontr o que el modelo ADM-1 no es sensible a la distribuci on de las razones de carbohidratos, prote nas y l pidos, sin embargo se apreci o que la fracci on de acidos grasos de cadena larga presentes en el inuente genera alta sensibilidad en los resultados.
En algunas aplicaciones previas de ADM-1 los par ametros del modelo fueron estimados a partir del ajuste de los datos experimentales [148, 144, 39], sin embargo tal estimaci on no es factible en el proceso de co-digesti on anaer obica, debido a la dicultad para encontrar un conjunto u nico de par ametros aplicable a todas las posibles combinaciones de las fracciones de carbohidratos, prote nas y l pidos presentes en el sustrato inuente, as como de las diferentes razones de residuo org anico biodegradable respecto a la biomasa anaer obica presente en el reactor [105]. En la revisi on bibliogr aca se encontraron aplicaciones en las cuales se realizaron variaciones al modelo ADM-1 original, o no se consideraron algunos factores como los procesos sicoqu micos o las funciones de inhibici on, no obstante en estas investigaciones obtienen resultados relativamente cercanos a las aplicaciones del modelo en su extensi on original [97, 131, 113].
6.2 Desarrollo de la simulaci on del modelo ADM-1
73
6.2.
Desarrollo de la simulaci on del modelo ADM-1
La gr aca 6-1 muestra la metodolog a establecida para el desarrollo de la implementaci on del modelo bajo los lineamientos de la estructura general del modelo original [14]. En la gura se observan las entradas del modelo como son las concentraciones iniciales de las variables que representan las especies solubles e insolubles (qu micas o biol ogicas), adem as del conjunto de coecientes de consumo, par ametros sicoqu micos, constantes cin eticas, entre otros, los cuales est an presentes en las ecuaciones que conforman el modelo. Como se muestra en la gr aca, estas entradas del modelo permiten el c alculo de las mismas variables en un tiempo posterior a trav es de las ecuaciones din amicas que representan los procesos de conversi on bioqu mica, adem as generan los valores de las concentraciones de las especies qu micas que afectan el equilibrio i onico ( acidos grasos vol atiles) y los gases disueltos que est an presentes en la fase l quida. A partir de estas concentraciones, el modelo realiza la iteraci on correspondiente para el c alculo de las concentraciones de los gases en la fase gaseosa (ecuaciones de transferencia l quido-gas) y de la concentraci on de los aniones y cationes presentes en la fase l quida (ecuaci on de balance de carga para el c alculo del pH). Finalmente el modelo realiza el control del proceso mediante las funciones de inhibici on las cuales calculan el factor de inhibici on que afecta algunas ecuaciones especicas, mediante las concentraciones de los metabolitos que generan el efecto inhibitorio. Para el desarrollo de la implementaci on del modelo ADM-1 en la plataforma mencionada, se utiliz o la metodolog a presentada en la gr aca 6-1 para representar las ecuaciones que describen los procesos de conversi on bioqu mica as como para los procesos sicoqu micos de equilibrio acido-base y transferencia de masa entre las fases l quida y gaseosa mediante ecuaciones diferenciales ordinarias y ecuaciones diferenciales algebraicas, desarrolladas en forma de diagramas de bloques, utilizando la versatilidad de la plataforma Matlab/Simulink. Los valores de los par ametros de entrada del modelo fueron tomados de los valores presentados en el reporte t ecnico del modelo ADM-1 (anexo E).
6.3.
Validaci on del modelo respecto al caso base ADM-1
Para realizar la validaci on del modelo se utiliz o otra implementaci on generada en la misma plataforma mediante la cual pueda establecerse la precisi on num erica de la implementaci on desarrollada mediante diagramas de bloques. Esta segunda implementaci on desarrollada por Rosen et al. [116], utiliza funciones MEX y c odigo C para la resoluci on del modelo ADM-1 incorporando un algoritmo especico para calcular la concentraci on del hidr ogeno en fase l quida (Sh2Solver) y el valor del pH (pHSolver). Los datos de la caracterizaci on del sustrato, la estequiometr a y los par ametros cin eticos son escritos en Excel de Microsoft y son le dos desde Matlab. Este ultimo utiliza vectores para almacenar los valores le dos para el conjunto total de par ametros. Simulink act ua como un diagrama de ujo en el cual las
74
Funciones de inhibici on debido al H2 y Ninorganico
Funciones de inhibici on debido al NH3 y pH
Concentraciones iniciales de las especies solubles e insolubles
Ecuaciones Din amicas 12 estados referentes a los componentes solubles (Si ) 5 estados referentes a los componentes insolubles (Xi )
Estados relacionados con la asociaci on y disociaci on de iones Concentraci on de CH4 , CO2 y H2 en la fase liquida
Calculo de la concentraci on de iones en la ecuaci on de equilibrio de cargas
Par ametros bioqu micos, sicoqu micos y de operaci on del reactor
7 estados referentes a los grupos de microorganismos (Xi )
Ecuaciones din amicas de transferencia de masa en la interfaz liquido-gas para CH4 , CO2 y H2
Figura 6-1.: Esquema general para el desarrollo de la implementaci on y simulaci on de los procesos bioqu micos y sicoqu micos que ocurren durante la digesti on anaer obica con base en el modelo ADM-1.
diferentes partes del modelo son conectadas mediante las razones de ujo. No obstante, el modelo implementado mediante la generaci on de las funciones MEX C-S en Matlab/Simulink presenta diferencias respecto al modelo original ADM-1 sugerido por Batstone et al. [14]; el cual fue implementado en esta investigaci on mediante diagramas de bloques, estas diferencias son descritas a continuaci on [116]. Se realizo la extensi on de las ecuaciones estequiom etricos para garantizar la consistencia de los balances de masa para nitr ogeno y carbono durante el proceso de digesti on anaer obica.
1
Se hizo la modicaci on de los valores dados por defecto para los par ametros (fxi,xc ,fli,xc , NI ,Nbac y Nxc ) con el n de corregir el desbalance inherente de nitr ogeno presente en la conguraci on original del modelo ADM-1. Se modico el valor por defecto para el par ametro(Cxc - contenido de carbono presente en el material compuesto) para corregir el desbalance inherente de carbono en el modelo ADM-1 original. Se modicaron las ecuaciones acido-base para mejorar las propiedades num ericas durante la implementaci on generada mediante los sistemas de ecuaciones diferenciales ordinarias que conforman el modelo ADM-1. Se utilizaron diferentes funciones de inhibici on (funciones Hill, [57] ) para evitar discontinuidades y los problemas num ericos asociados a estas.
6.4 Vericaci on del modelo respecto al caso experimental
75
Mediante esta implementaci on se realiz o la simulaci on para la obtenci on de las curvas de comportamiento de la concentraci on de algunas de las especies qu micas. Estas fueron comparadas respecto a la implementaci on desarrollada mediante diagramas de bloques para las ecuaciones diferenciales ordinarias y algebraicas que conforman el modelo. De esta forma se realizo la evaluaci on del comportamiento del modelo bajo las dos implementaciones mencionadas, utilizando los valores reportados para el caso base de ADM-1 [14]. De acuerdo los resultados de esta comparaci on, (presentados en la gura 6-2), se puede observar la precisi on y la convergencia relativa de los valores obtenidos mediante las dos implementaciones; esto permite asegurar la conveniencia de utilizar la implementaci on del modelo ADM-1 basada en diagramas de bloques, con el n de simular el proceso de biodegradabilidad anaer obica desarrollado en la etapa experimental de esta investigaci on.
CH (CB-A)
4
CH4 (CB-B) CO2 (CB-A) Producci on de metano [g DQO-CH4 /L] 1 .2 1 .1 1 .0 0 .9 0 .8 0 .7 0 .6 0 .5 0 .4 0 .3 0 .2 0 .1 0 0 5 10 15 20 25 30 Tiempo [d as] (CB-A) (CB-B) 35 40 Fracci on Molar de metano y Di oxido de carbono CO2 (CB-B) 0 .9 0 .8 0 .7 0 .6 0 .5 0 .4 0 .3 0 .2 0 .1 0 0 5 10 15 20 25 30 35 40 Tiempo [d as]
(a)
(b)
Figura 6-2.: Curvas generadas para la DQO-CH4 (a) y la composici on aproximada del biog as generado (b), mediante la simulaci on num erica en Matlab/Simulink del caso base del modelo ADM-1 utilizando la implementaci on con diagramas de bloques [A] y la implementaci on a trav es de funciones C-S MEX [B].
6.4.
Vericaci on del modelo respecto al caso experimental
Con el n de realizar la 1 implementaci on del modelo ADM-1 para simular las condiciones 1 en las cuales se desarroll o el proceso de digesti on anaer obica durante la fase experimental de este estudio, se hizo necesario establecer la composici on de la mezcla de residuos
76
org anicos biodegradables (esti ercol bovino y residuos c tricos). Debido a la alta complejidad de los residuos org anicos utilizados, result o fuera del alcance de este estudio la realizaci on de los an alisis de composici on de estos residuos. Por lo tanto como primera aproximaci on se consult o la literatura para establecer los par ametros de entrada del modelo (contenidos de carbohidratos, prote nas, l pidos, azucares, amino acidos, acidos grasos de cadena larga, entre otros). Algunos estudios presentan resultados obtenidos para las composiciones estimadas de algunos residuos org anicos biodegradables con base en el desarrollo de mediciones especicas realizadas sobre el sustrato. Dentro de estos estudios se puede se nalar el realizado por Zaher et al. [148], quienes desarrollaron un modelo integrado para la simulaci on del proceso de digesti on anaer obica de diversos residuos org anicos mediante la implementaci on del modelo ADM-1. En ese estudio previamente se caracterizaron los residuos utilizando mediciones estandarizadas dentro de las que se incluyen demanda qu mica de ox geno, pH, alcalinidad total y parcial, contenidos totales de carbono y nitr ogeno total Kjeldahl. Con base en estas mediciones Zaher et al. evaluaron los resultados obtenidos mediante la implementaci on del algoritmo que transforma los valores de estas mediciones realizadas sobre el sustrato y los valores obtenidos mediante un estudio de caracterizaci on de los residuos en el cual se establecieron los contenidos de carbohidratos, prote nas y l pidos entre otros. Para la obtenci on de la composici on de los residuos org anicos como el esti ercol bovino estos autores realizaron el montaje de varias t ecnicas de medici on dentro de las cuales se incluyen la extracci on secuencial usando detergente acidicado y neutral, seguido de una extracci on utilizando un acido fuerte para la estimaci on del contenido de carbohidratos; el contenido prot eico fue estimado mediante el montaje de la t ecnica colorim etrica de Lowry, calibrada a trav es una curva patr on de alb umina extra da del suero bovino. De igual forma el contenido de l pidos presentes en el esti ercol bovino fue medido por medio de la t ecnica de extracci on Soxhlet con eter como solvente. La tabla (6-1) resume los resultados de la composici on del esti ercol bovino obtenidos por Zaher et al. [148]. Adicionalmente se hizo necesario establecer la composici on aproximada de los desechos c tricos utilizados como co-sustrato durante las pruebas de biodegradabilidad anaer obica. En la revisi on bibliogr aca se encontr o un estudio en el cual se utiliz o espec camente el residuo c trico procedente de la extracci on del jugo de la naranja variedad valencia, para la elaboraci on de un suplemento dietario para la alimentaci on de especies menores en una explotaci on agroindustrial [119]. De acuerdo a la metodolog a empleada en dicho estudio, se realiz o un pretratamiento sobre los residuos c tricos; los cuales fueron escaldados por un tiempo de 5 o minutos a 105 C, posteriormente se enfriaron empleando agua a 20o C durante 3 minutos y nalmente se realiz o la operaci on de secado en un t unel de aire forzado con una temperatura o promedio de 60 C y una velocidad del aire de 2 m/s lo cual permiti o que el residuo alcanzara un contenido de humedad m aximo de 25 %. Finalmente el residuo c trico fue molido hasta obtener un tama no de grano entre 500m y 600m. La tabla 6-2 resume la composici on qu mica del residuo c trico de acuerdo a los resultados reportados en el estudio en menci on.
77
Tabla 6-1.: Caracter sticas sicoqu micas esti ercol bovino [148].
Par ametro Demanda Qu mica de Ox geno Total (DQOT ) Demanda Qu mica de Ox geno Particulada (DQOP ) Demanda Qu mica de Ox geno Soluble (DQOS ) DQOS sin DQO-AGV (DQOS -AGV) Acidos Grasos Vol atiles (AGV) Carbono Total (CT) Carbono Org anico Total (COT) Carbono Inorg anico Total (CIT) Nitr ogeno Total Kjeldahl (NTK) Nitr ogeno Org anico Total (NOT) Nitr ogeno Amoniacal Total(NAT) F osforo Total (FT) F osforo Org anico Total (FT-OrtoP) Orto-Fosfato (OrtoP) Alcalinidad Total S olidos Totales (ST) S olidos Fijos (SF) S olidos Totales Vol atiles (STV) Carbohidratos Prote nas L pidos Unidades mg/L mg/L mg/L mg/L mg/L mg C/L mg C/L mmol HCO 3 /L mg N/L mg N/L mg N/L mg P/L mg P/L mg P/L milieq./L mg/L mg/L mg/L mg/L mg/L mg/L Valor medido 27217 23550 3667 2521 1146 10064 9340 60 882 598 284 219 187 32 60 20697 5397 15300 10924428 4069367 30661.2
Con base en estas composiciones obtenidas de la literatura se intent o establecer una relaci on directa entre el contenido de s olidos y las fracciones de cada especie qu mica contenidas en el sustrato que se utiliz o para el desarrollo de las pruebas de biodegradabilidad anaer obica. Para este n se utiliz o la aproximaci on presentada en un estudio que desarrolla la implementaci on del modelo ADM-1 para la simulaci on del proceso de digesti on anaer obica de los residuos org anicos provenientes de una destiler a de vinos, tomando como componente principal el etanol [113].
Tabla 6-2.: Caracter sticas sicoqu micas residuos c tricos [119].

Componente Porcentaje Carbohidratos 26.2 %0.1 Prote nas 6.7 %0.2 L pidos 0.9 %0.1 F DI Relaci on( F DS ) Componente Cenizas Fibra Dietaria Insoluble (FDI) Fibra Dietaria Soluble (FDS) 5.34:1 Porcentaje 2.8 %0.1 53.4 %1.1 10 %0.4
78
En ese estudio el autor transforma las ecuaciones del modelo ADM-1 con el n de incorporar la conversi on bioqu mica del etanol contenido en este residuo especico. De esta forma este modelo presenta el desarrollo de las ecuaciones en base molar a diferencia del modelo original ADM-1 el cual esta desarrollado con base en unidades de DQO [g DQO/L]. La tabla 6-3 fue obtenida de la investigaci on mencionada y muestra la composici on elemental de cada una de las especies qu micas del modelo ADM-1 con base molar. A partir del valor relacionado en la tabla para la demanda qu mica de ox geno de cada especie biol ogica o qu mica, se estableci o la concentraci on de cada variable de estado de acuerdo a los valores medidos para los contenidos de s olidos presentes en la mezcla de residuos. Con base en los resultados obtenidos para los valores iniciales de cada una de las especies (variables de estado) y utilizando los valores reportados para cada uno de los par ametros y constantes que est an presentes en el modelo de acuerdo al caso base de ADM-1, se desarrollo la simulaci on del proceso de digesti on anaer obica. Los resultados obtenidos mediante esta simulaci on dieren bastante de las mediciones desarrolladas en el laboratorio para todos los par ametros de control del proceso. De acuerdo a esta comparaci on de resultados se inere que los valores de las concentraciones iniciales de las especies qu micas, deben determinarse mediante otra metodolog a. No se obtienen resultados congruentes al relacionar los contenidos de s olidos (ST, STV, SST, SSV) con el contenido de carbohidratos, prote nas, l pidos, azucares, amino acidos o acidos grasos vol atiles de cadena larga presentes en el residuo. Esto se debe a la presencia de materiales complejos tales como la celulosa, hemicelulosa y lignina en el esti ercol bovino y como las bras dietarias solubles e insolubles en el desecho c trico, las cuales afectan el valor de la medici on del contenido de s olidos y no permiten la determinaci on directa de los contenidos de los compuestos considerados en la implementaci on del modelo ADM-1. Por lo tanto se hizo necesario establecer otra metodolog a para la determinaci on de las concentraciones iniciales de los metabolitos mencionados, con el n de realizar la simulaci on del proceso de digesti on anaer obica desarrollado en la etapa de experimentaci on. La metodolog a que permite establecer estos contenidos se hall o en un estudio desarrollado sobre el comportamiento del proceso de digesti on anaer obica cuando se realiza la codigesti on de aguas de desecho de origen dom estico y desechos org anicos de cocina [144]. En ese estudio se presenta una consideraci on importante sobre la cual se desarrolla el modelamiento; de acuerdo a este autor, es posible fraccionar la medici on de la demanda qu mica de ox geno total (DQOT otal ) en componentes menores como la demanda qu mica de ox geno particulada (DQOP ) y la demanda qu mica de ox geno soluble (DQOS ). A su vez cada una de estas fracciones presenta algunos componentes principales. La demanda qu mica de ox geno particulada es diferenciada en fracciones biodegradable y no biodegradable. Similarmente la demanda qu mica de ox geno soluble es diferenciada en sus componentes acidicada y no acidicada. En ese estudio se establecen los porcentajes que representan cada uno de estos componentes respecto a
79
la fracci on principal. Aunque estos porcentajes han sido establecidos con base en mediciones sobre un sustrato en particular, fueron utilizados como una aproximaci on para determinar las concentraciones iniciales de los compuestos que est an presentes en la mezcla de residuos bajo estudio. Tabla 6-3.: Composici on elemental de las especies modelo ADM-1 base molar [113].
Especie Glucosa Amino acidos (Globalizados) AGCL Acido Val erico Acido But rico Acido Propi onico Acido Ac etico Hidr ogeno Metano Di oxido de Carbono Amon aco Inertes Solubles Sustrato Particulado Carbohidratos S olidos Prote nas Solidas L pidos S olidos Consumidoras de Glucosa Consumidoras de Amino acidos Consumidoras de AGCL Consumidoras de Butirato Consumidoras de Propionato Consumidoras de Acetato Consumidoras de Hidr ogeno S olidos Inertes S mbolo Ssu Saa Sfa Sva Sbu Spro Sac Sh2 Sch4 Sic Sin Si Xc Xch Xpr Xli Xsu Xaa Xfa Xbu Xpro Xac Xh2 Xi C 6 1 1 5 4 3 2 0 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 H 12 1.9 1.9 10 8 6 4 2 4 0 3 1.946 1.5 2 1.9 1.9 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 O 6 0.6 0.1 2 2 2 2 0 0 2 0 0.6754 0.5324 1 0.6 0.1 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4024 N 0 0.225 0 0 0 0 0 0 0 0 1 0.1429 0.0895 0 0.2253 0 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.1429 (gDQO/mol) 192 32.19 45.6 208 160 112 64 16 64 0 0 33.3 33.3 32 32.2 45.6 32 32 32 32 32 32 32 33.3
La gr aca 6-3 presenta la relaci on existente entre la demanda qu mica de ox geno total y sus componentes particulada y soluble. Adem as muestra la distribuci on porcentual de los metabolitos que conforman cada uno de estos componentes. A partir de las mediciones desarrolladas en la fase experimental de este estudio como la demanda qu mica de ox geno (total y soluble) y contenido de acidos grasos vol atiles, se estimaron las concentraciones de cada uno de los compuestos considerados en el modelo ADM-1 (especies qu micas) utilizando los porcentajes establecidos en la gr aca 6-3, para cada raz on de carga respectivamente. Los valores de cada variable de entrada del modelo, que representan un metabolito contenido en el sustrato fueron calculados para cada una de las razones de carga ensayadas, mientras que las variables de estado que representan las especies biol ogicas fueron tomadas del reporte original de ADM-1. Con base en estos nuevos valores calculados se realizaron las simulaciones
80
para cada una de las razones de carga evaluadas mediante el desarrollo experimental. De esta forma fue posible obtener un conjunto de gr acas que muestran las variaciones de la concentraci on de estos productos intermedios durante el desarrollo del proceso de digesti on anaer obica. Los resultados obtenidos con esta nueva simulaci on se encuentran m as cercanos a los resultados medidos experimentalmente para los par ametros de control del proceso como pH, alcalinidad, acidos grasos vol atiles, demanda qu mica de ox geno total y soluble as como para la producci on total de biog as y producci on de metano. Estas gr acas obtenidas mediante la simulaci on se compararon con las gr acas obtenidas a partir del conjunto de datos medidos durante el desarrollo de los ensayos de biodegradabilidad anaer obica y se muestran en el capitulo concerniente a los resultados y an alisis.
DEMANDA QU IMICA DE OXIGENO TOTAL (DQOT )
Demanda Qu mica de Ox geno Particulada (DQOP )
Demanda Qu mica de Ox geno Soluble (DQOS )
Degradables (63 % (DQOP ) )
Inertes (37 % (DQOP ))
Acidos Grasos Vol atiles (AGV)
Material Disuelto Excepto AGV
Carbohidratos (32 %) de Degradables
Proteinas (32 %) de Degradables
L pidos (36 %) de Degradables
Monosac aridos (33 %) de Disueltos no AGV
Amino acidos (33 %) de Disueltos no AGV AGCL (33 %) de Disueltos no AGV
Inertes Disueltos (18 %) de Inertes
Inertes Particulados (82 %) de Inertes
Acido ac etico (63 %) de AGV
Acido Propi onico (24 %) de AGV
Acido But rico (8 %) de AGV
Acido Val erico (5 %) de AGV
AGCL: Acidos grasos de cadena larga
Figura 6-3.: Fraccionamiento de la demanda qu mica de ox geno total -DQOT , a partir de la cual se obtuvieron los valores de las condiciones iniciales para la implementaci on del modelo ADM-1 [144].
1
7. Resultados y An alisis
7.1. Resultados Experimentales
pH, alcalinidad y acidos grasos vol atiles (AGV)

La tabla 7-1, presenta los valores obtenidos mediante las mediciones realizadas de los par ametros de control del proceso de acuerdo a la descripci on dada en el capitulo correspondiente al desarrollo experimental. Dentro de este conjunto de par ametros medidos se encuentran el pH, concentraci on de acidos grasos vol atiles y alcalinidad. Como se estableci o anteriormente, se tomaron cinco (5) mediciones en conjunto para cada uno de los par ametros mencionados para las cargas org anicas correspondientes a los valores (1:1) y (3:1); para la carga evaluada (2:1) se realizaron tres (3) mediciones de cada par ametro debido a que mediante la planeaci on experimental se deni o la realizaci on exclusiva de las mediciones inicial, intermedia y nal. La gura 7-1-(a) presenta los perles obtenidos del conjunto de datos experimentales para el pH del sustrato contenido en los reactores correspondientes a cada una de las razones de carga org anica utilizadas en este estudio. A partir de las curvas establecidas por este conjunto de mediciones se observa que los valores del pH se mantuvieron en el rango entre pH 6.25 y pH 7.5, el cual se encuentra dentro del rango optimo reportado para el desarrollo eciente del proceso de digesti on anaer obica (pH entre 6.5 y 7.5) [69, 82]. A pesar de la limitada variabilidad que presenta el pH, en la gr aca 7-1-(a) se observa un ligero decremento en el valor del pH para todas las razones de carga durante los primeros d as de desarrollo de los ensayos de biodegradabilidad anaer obica. Como puede observarse en la gura 7-1-(b) este descenso inicial del valor de pH es consecuente con el aumento gradual de los acidos grasos vol atiles durante este mismo periodo. Estas tendencias observadas para el pH y los acidos grasos vol atiles son ndices de la presencia de constituyentes f acilmente biodegradables en el sustrato, los cuales son r apidamente utilizados por los microorganismos fermentativos presentes en el medio anaer obico [30]. Hacia el periodo nal del ensayo se observa un aumento del pH para todas las razones de carga ensayadas. El comportamiento observado para el pH es coherente con los resultados reportados por Geogarcakis et al. [47], quienes determinaron tendencias similares para el pH al desarrollar la digesti on anaer obica de esti ercol bovino. De acuerdo a ese estudio el valor del pH depende de la resistencia que 2 + presentan los sistemas buer (HCO como el contenido de acidos 3 /CO3 y NH4 /NH3 ) as grasos vol atiles. Cuando el contenido de estos acidos disminuye producto de la formaci on del metano, se observa un leve incremento del pH como consecuencia de la acci on de los sistemas amortiguadores.
82
7 Resultados y An alisis
Tabla 7-1.: Valores medidos para el pH, alcalinidad y acidos grasos vol atiles durante el ensayo de biodegradabilidad anaer obica. Raz on de carga 1:1 Tiempo (d as) 0 10 20 30 40 0 20 40 0 10 20 30 40 pH [Unidades de pH] 7.5 7.2 6.56 6.7 6.9 7.4 6.46 6.6 7.1 6.54 6.52 6.53 6.6 AGV [mg CH3 COOH/L] 2110.1 2340.7 1796.1 1246.2 1034.2 2508.1 1834.3 867.1 Raz on de carga 3:1 3624.2 4321.3 1303.3 873.8 628.1 3918 4441 3657 3135 3396 0.92 0.97 0.35 0.28 0.18 Alcalinidad [mg CaCO3 /L] 3135 2873 2351 2612 2873 2873 2351 2612 AGV/Alcalinidad [mg CH3 COOH /mg CaCO3 ] 0.67 0.8 0.76 0.47 0.36 0.87 0.78 0.33
Raz on de carga 2:1
La gr aca 7-1-(b) muestra la evoluci on en funci on del tiempo de la concentraci on de acidos grasos vol atiles totalizados como la concentraci on de acido ac etico. De acuerdo con las concentraciones iniciales de los acidos grasos vol atiles, estos acidos est an presentes normalmente en el sustrato al inicio de los ensayos de biodegradabilidad anaer obica, y son gradualmente descompuestos durante el desarrollo del proceso. En la gr aca se pueden identicar los patrones de aumento de la concentraci on de estos acidos org anicos en las razones de carga (1:1) y (3:1) durante los primeros diez d as despu es de iniciado el ensayo, lo cual es evidencia del desarrollo eciente de las etapas iniciales de hidr olisis y fermentaci on. La razones de carga (3:1 y 1:1) presentan la m axima y la m nima concentraci on de acido grasos respectivamente; mientras que para la raz on de carga (2:1) no es posible determinar una tendencia de aumento en ninguna de las fases del proceso debido a que u nicamente se tomaron las mediciones de AGV en tres oportunidades como se explic o en la descripci on de la etapa experimental. La raz on de carga (3:1) muestra un r apido decremento de la concentraci on de los acidos una vez alcanza el m aximo, esto podr a ser consecuencia de la alta disponibilidad de estos acidos para las bacterias acetog enicas y metanog enicas que los utilizan directamente como sustrato.
7.1 Resultados Experimentales
83
La raz on de carga (1:1) presenta una menor raz on de disminuci on de la concentraci on de los acidos org anicos despu es de alcanzado el m aximo, lo cual es atribuible a la falta de alg un nutriente requerido para el desarrollo optimo de la biomasa microbiana responsable del consumo de estos acidos. La disminuci on del contenido de acidos grasos vol atiles hacia el periodo nal de los ensayos esta en concordancia con la disminuci on progresiva de la generaci on de metano y biog as como se muestra posteriormente. Esta tendencia en el comportamiento de las curvas de concentraci on de los acidos grasos vol atiles es muy similar a la mostrada en la gura 5-2-(d) [112], la cual muestra la variaci on de este par ametro durante ensayos de biodegradabilidad anaer obica sobre residuos agroindustriales similares al utilizado en esta investigaci on. Con base en la tendencia general de la curva se observa que el proceso de degradaci on de la materia org anica se desarrollo de forma adecuada hasta el proceso de la metanog enesis, alcanzando valores similares para la concentraci on nal de estos acidos en las tres razones de carga org anica ensayadas. El comportamiento de la alcalinidad en cada una de las cargas utilizadas en el ensayo de biodegradabilidad anaer obica es mostrado en la gura 7-1-(c). En esta gr aca se observa la diferencia en los valores iniciales de la alcalinidad para cada una de las razones de carga evaluadas. Este par ametro presenta uctuaciones importantes debido a la generaci on de metabolitos como los acidos grasos vol atiles durante el desarrollo de la digesti on anaer obica; los aniones procedentes de los acidos d ebiles generados interact uan con los iones hidr onio + [H ], lo cual da lugar a reacciones de neutralizaci on en el medio anaer obico. De acuerdo con la gura 7-1-(c) la raz on de carga (3:1) presenta una capacidad amortiguadora considerable la cual se reeja en el valor inicial de la alcalinidad y el aumento progresivo de esta durante los primeros d as de desarrollo del proceso de digesti on anaer obica. Sin embargo esta raz on de carga presenta la mayor acumulaci on de acidos grasos vol atiles. Esta tendencia puede ser explicada debido al complejo equilibrio existente entre algunos componentes de la fase l quida como los aniones procedentes de los acidos grasos vol atiles, adem as de los aniones 2+ HCO3 , CO3 y la presi on parcial de CO2 en la fase gaseosa, los cuales afectan directamente la alcalinidad [63]. La raz on de carga (1:1) presenta un ligero descenso de la alcalinidad hacia el inicio del ensayo hasta aproximadamente la mitad del tiempo del ensayo y un aumento gradual hacia el periodo nal de la prueba, lo cual es consistente con el aumento del pH y la disminuci on del contenido de acidos grasos vol atiles en este mismo periodo.
Demanda qu mica de ox geno soluble

La tabla 7-2 muestra los valores obtenidos en la medici on de la demanda qu mica de ox geno soluble (DQOS ) sobre cada uno de los reactores que se descartaron con una periodicidad de 10 d as para cada una de las razones de carga org anicas utilizadas en este estudio. Con base en los valores medidos para este par ametro se estableci o el porcentaje m aximo de remoci on de la demanda qu mica de ox geno soluble, como un indicador de la efectividad del tratamien-
84
8.00 7.75 pH [Unidades de pH] 7.50
(1:1-A) (2:1-B) (3:1-C)
5000 Acidos Grasos Vol atiles [mg CH3 COOH/L] 4500 4000 3500 3000 2500 2000 1500 1000 500 0 0 5 10 15 20 25 Tiempo [d as] 30
(1:1) (2:1) (3:1)
7.25
7.00
6.75
6.50
6.25 6 0 5 10 15 20 25 30 35 40 Tiempo [d as]
35
40
(a)
5000
(b)
Demanda qu mica de ox geno removida [mg O2 /L] 3000
4500 4000 Alcalinidad [mg CaCO3 /L] 3500 3000 2500 2000 1500 1000 500 0 0 5 10 15 20 25 Tiempo [d as] 30 35 40 (1:1) (2:1) (3:1)
(1:1) (2:1) (3:1)
2500
2000
1500
1000
y = 41.2x
1
500 y = 32.5x 0 0 5 10 15 20 25 Tiempo [d as] 30 35 40
(c)
(d)
Figura 7-1.: Variaci on de la concentraci on de los par ametros de control del proceso durante el desarrollo de las pruebas de biodegradabilidad anaer obica. (a). pH, (b). Acidos grasos vol atiles, (c.) Alcalinidad, (d). Variaci on de la DQOS -Removida en los reactores de tratamiento.
to anaer obico sobre la estabilizaci on del residuo org anico biodegradable. De acuerdo con los 1 valores calculados para este ndice, se obtuvieron porcentajes m aximos de remoci on de la DQOS de 56.3 %, 51.3 % y 48.0 % para las cargas org anicas analizadas (1:1), (2:1) y (3:1) 1 respectivamente. A partir de las diferencias en los valores medidos para la demanda qu mica de ox geno soluble se pueden establecer los valores para la demanda qu mica de ox geno removida DQOS-REM ; la gr aca 7-1-(d) presenta estos valores as como las ecuaciones de las rectas que pasan por los puntos en los cuales se encuentran los valores de las pendientes m aximas de estas curvas.
7.2 Producci on de metano y biog as
85
Tabla 7-2.: Variaci on de la demanda qu mica de ox geno total y soluble (DQOT y DQOT ) durante los ensayos de biodegradabilidad anaer obica. N.D.: No determinado, % Rem.= Porcentaje de remoci on. Raz on de carga 1:1 Tiempo (d as) 0 10 20 30 40 DQOT 3075 2749 2539 2209 1227 DQOS 2786 2446 2373 2012 1215 % Rem. [ %] 0 12.2 14.8 27.7 56.3 [mg O2 /L] Raz on de carga 2:1 DQOT 4422 N.D. 2217 N.D. 1933 DQOS 3751 N.D. 2155 N.D. 1826 % Rem. [ %] 0 N.D. 42.5 N.D. 51.3 [mg O2 /L] Raz on de carga 3:1 DQOT 5427 4929 4245 3235 2892 DQOS 5244 4832 4115 3032 2724 % Rem. [ %] 0 7.85 21.5 42.1 48.0 [mg O2 /L]
La tasa m axima de remoci on de la DQOS es un par ametro que permite establecer la velocidad con la cual los microorganismos descomponen el sustrato org anico. Este par ametro se estableci o como la pendiente m axima de la parte inicial de la curva de remoci on de la DQOS , dividido en la concentraci on inicial de s olidos suspendidos vol atiles (SSV) de cada una de las cargas org anicas utilizadas en el ensayo. La tabla 7-3 presenta los valores calculados para la tasa m axima de remoci on de la DQOS para las razones de carga org anica (1:1 y 3:1), obtenidos a partir de la pendientes dadas en la gr aca 7-1-(d) y los valores del contenido de SSV dados en la misma tabla. De acuerdo a estos resultados se observa que la carga ensayada (3:1) presenta el menor porcentaje de remoci on de la DQOS , sin embargo revela la mayor tasa de remoci on de la DQOS . Por su parte la raz on carga (1:1) presenta la menor velocidad de remoci on y la mayor tasa de remoci on. Es decir el tratamiento realizado con la raz on de carga (1:1) es m as efectivo para la reducci on de la demanda qu mica de ox geno soluble, no obstante el tratamiento realizado con la raz on de carga (3:1) realiza la reducci on de este par ametro a una mayor velocidad. Pese a que los resultados obtenidos para esta raz on de remoci on se encuentran dentro del rango recomendado en la literatura (0.1 - 1.0 [g DQOREM /(g SSVd a)]); estos se encuentran lejos del promedio ideal de (0.5 [g DQOREM /(g SSVd a)]) [110].
7.2.
Producci on de metano y biog as
La gura 7-2-(a) muestra el comportamiento en funci on del tiempo de la producci on acumulada neta de metano y biog as durante el ensayo de biodegradabilidad anaer obica para cada una de las cargas evaluadas, los valores presentados corresponden al promedio de las mediciones realizadas en los reactores disponibles durante el tiempo de desarrollo del ensayo. De
86
Tabla 7-3.: Tasa m axima de remoci on de la DQOS en los ensayos de biodegradabilidad anaer obica. Raz on de carga (g DQO/g SSV) 1:1 3:1 DQOS-REMOVIDA [mg O2 /Ld a] 32.5 41.2 SSV [mg/L] 1943.2 2230 Tasa m axima de remoci on [g DQOREMOVIDA /(g SSV d a)] 0.016 0.018
acuerdo a la gura se puede observar el retraso que se presenta en los reactores de tratamiento para todas las cargas, debido a que la biomasa microbiana presente en el lodo anaer obico requiere de un periodo de adaptaci on para iniciar los procesos subsecuentes de la degradaci on anaer obica del residuo org anico. Despu es de este periodo inicial de adaptaci on, se observa un r apido incremento en la producci on de biog as y gas metano para todas las cargas, previo a la fase nal en la cual la raz on de producci on tiende a disminuir considerablemente. De acuerdo con las tendencias presentadas por cada una de las curvas, se observa que la raz on de carga evaluada (1:1) presenta la tasa m as baja de generaci on de metano y biog as, no obstante mantiene esta misma raz on de generaci on por casi m as de 20 d as; al nal de este periodo se presenta una disminuci on de la producci on de biog as y metano respectivamente, lo que indica el nal del proceso de digesti on anaer obica para esta raz on de carga. En la gr aca tambi en puede observarse que esta raz on de carga es la que presenta el tiempo de adaptaci on m as largo. Esta caracter stica puede atribuirse a que esta raz on carga presento la menor concentraci on de acidos grasos vol atiles, lo cual implica que la producci on de metano se retrasa debido al requerimiento de la producci on de estos acidos org anicos en una etapa anterior a la etapa de metanog enesis. En los reactores de tratamiento correspondientes a la raz on de carga (2:1) se observan las mayores tendencias de generaci on de metano y biog as, las cuales inician aproximadamente a partir del s eptimo d a del ensayo, y contin uan hasta aproximadamente el d a quince; posteriormente estas razones de producci on decrecen hasta el periodo de culminaci on de la prueba. De igual forma raz on de carga (3:1) presenta una velocidad intermedia para la generaci on de metano y biog as la cual inicia a partir del d a quince, y se mantiene hasta un tiempo cercano al d a 22 de ensayo, tiempo en el cual comienza a decrecer la producci on tanto de biog as como de metano. Es notable la similitud que presentan en conjunto estas gr acas con las gr acas mostradas anteriormente en la gura 5-1-(a) [35], las cuales representan el comportamiento de la producci on de biog as a partir de la digesti on anaer obica de un residuo similar.
Tasa m axima y Porcentaje de biodegradabilidad anaer obica

La gura 7-2-(b) muestra las curvas obtenidas para la producci on neta de metano en los reactores de tratamiento anaer obico, para cada una de las razones de carga en t erminos
87
de demanda qu mica de ox geno biometanizada (DQOCH4 ) [g DQOCH4 /L]. Esta gr aca se obtuvo a partir de los datos consignados en la gr aca 7-2-(a), mediante la aplicaci on del factor de conversi on establecido para el c alculo de la demanda qu mica de ox geno de un volumen de metano bajo las condiciones de presi on y temperatura a las cuales se realizo el ensayo. Adicionalmente la gr aca muestra las ecuaciones de las rectas establecidas sobre los puntos de la curva en los cuales se presenta la mayor raz on de generaci on de metano. A partir de las pendientes de estas rectas se establecieron los valores del par ametro de tasa m axima de biodegradabilidad anaer obica o tasa m axima de producci on de metano (MMPR, por sus siglas en ingl es), la cual corresponde a la pendiente m axima de la parte lineal al inicio de la curva de producci on de metano expresada en t erminos de DQO [mg DQO-CH4 /Ld a] divido en el contenido de los s olidos suspendidos vol atiles (SSV) de cada raz on de carga [98]. Este par ametro es un ndice de la velocidad m axima de generaci on de metano respecto al contenido de biomasa celular en el sustrato. La tabla 7-4 muestra los valores calculados de este par ametro para cada una de las razones de carga; como se observa en esta tabla, las razones de carga org anica (3:1 y 1:1) presentan la mayor y la menor tasa de biodegradabilidad anaer obica respectivamente. Este resultado es consecuente con los valores obtenidos para el par ametro tasa m axima de remoci on de la DQOS , seg un el cual la raz on de carga (3:1) y (1:1) presentaron la mayor y la menor tasa de remoci on de la DQOS respectivamente. El porcentaje de biodegradabilidad anaer obica mostrado en la tabla 7-6, se c alculo mediante la utilizaci on del valor reportado para la m axima producci on de metano a partir de la digesti on anaer obica del esti ercol bovino [148 mL CH4 /g STV] [91], y los valores m aximos establecidos en la gura 7-2-(c) para la producci on de metano de cada una de las tres razones de carga evaluadas. De acuerdo a este par ametro se obtuvo el mayor porcentaje de biodegradabilidad anaer obica a partir de la raz on de carga (2:1), sin embargo las dem as razones de carga presentan valores similares para este ndice. Tabla 7-4.: Tasa m axima de biodegradabilidad anaer obica para cada una de las cargas org anicas ensayadas. Raz on de carga (g DQO/g SSV) 1:1 2:1 3:1 DQO-CH4 [g DQO/Ld a] 23.1 42.2 67.5 SSV [mg SSV/L] 1943.2 2120.1 2230.4 MMPR [g DQO-CH4 /(g SSV d a)] 0.011 0.019 0.03
88
900 850 800 750 Producci on neta de metano (CH4 ) y Biog as [mL] 700 650 600 550 500 450 400 350 300 250 200 150 100 50 0 0 5 10 15 20 25 30 35 40 Tiempo [d as]
(1:1-A) (1:1-B) (2:1-A) (2:1-B) (3:1-A) (3:1-B)
Producci on de metano [mg DQO-CH4 /L]
1000 900 800 700 600 500 400 300 200 100 0 0
(1:1) (2:1) (3:1)
y=67,5(x)-892,24
y=42,2(x)-255,31
y=23,1(x)-222,1
10
15
20
25
30
35
40
Tiempo [d as]
(a)
Producci on neta de metano (mL CH4 /g STVadicionados ) 125
(b)
(1:1) (2:1) (3:1)
100
75
50
25
0 0 5 10 15 20 25 30 35 40
Tiempo [d as]
(c)
Figura 7-2.: Generaci on de productos nales en los reactores de tratamiento durante los ensayos de biodegradabilidad anaer obica. (a). Producci on volum etrica de metano y biog as [A]= biog as, [B]= metano (b). Producci on de metano en t erminos de [DQO-CH4 ]. (c.) Producci on neta de metano en t erminos de [mL CH4 /g STV].
89
Factor de conversi on de sustrato a metano

Este par ametro permite establecer la producci on de metano generada por cada gramo de demanda qu mica soluble de ox geno removida. La tabla 7-5 presenta los valores calculados para este par ametro en funci on de los valores obtenidos previamente para la producci on de metano en t erminos de DQOCH4 y la DQOS-REM . De acuerdo a los valores consignados en la tabla, el mayor valor y el menor valor para este par ametro lo presentan la raz on de carga org anica (2:1) y la raz on de carga org anica (1:1) respectivamente. Por consiguiente de acuerdo a estos datos la raz on de carga (2:1) presentar a un mejor desempe no en el proceso de digesti on anaer obica al obtener un ndice m as alto de producci on de metano por cada gramo de demanda qu mica de ox geno soluble removida. Sin embargo para la raz on de carga (3:1) el valor para este par ametro es muy cercano al valor de la raz on de carga (2:1), por lo tanto esta carga presentar a un desempe no similar en cuanto a la producci on de metano. Tabla 7-5.: Factor de conversi on del sustrato a metano en los ensayos de biodegradabilidad anaer obica. Raz on de carga (g DQO/g SSV) 1:1 2:1 3:1 DQO-CH4 (M axima) [mg O2 /L] 519.36 746.88 926.4 DQOS (Inicial) [mg O2 /L] 2786 3751 5244 DQOS (Final) [mg O2 /L] 1215 1826 2724 DQOS (Removida) [mg O2 /L] 1571 1925 2520 Factor [g DQO-CH4 / g DQO-(Removida)] 0.33 0.38 0.36
Balance de la demanda qu mica de ox geno

La tabla 7-6 presenta la evaluaci on de los par ametros de eciencia del proceso de digesti on anaer obica; los cuales son desarrollados a partir de los resultados obtenidos de DQOS , DQOAGV y DQOCH4 discutidos anteriormente. De forma similar las gr acas mostradas en la gura 7-3-(a, b, c) presentan estos mismos resultados en forma de lineas de tendencia, mediante las cuales se observa la relaci on existente entre la demanda qu mica de ox geno soluble y sus componentes: la demanda qu mica de metano generado y la demanda qu mica de acidos grasos vol atiles presentes en el medio anaer obico. En estas gr acas se observa como la demanda qu mica de ox geno soluble (DQOS ) disminuye, mientras que la demanda qu mica de ox geno correspondiente al metano generado aumenta a medida que avanza el proceso de digesti on anaer obica. Este comportamiento permite establecer la capacidad que presenta la digesti on anaer obica para disminuir la carga contaminante del sustrato al mismo tiempo genera un recurso energ etico renovable.
90
Tabla 7-6.: Balance de la DQOS y porcentaje de biodegradabilidad anaer obica. Raz on de carga org anica 1:1 %M. 0 1.37 8.63 15.54 18.64 Raz on de carga org anica 2:1 %M. %AGV %A. %Remoci on DQOS 0 71.11 71.11 0 13.71 52.16 65.88 42.54 19.91 24.66 44.57 51.31 Raz on de carga org anica 3:1 %M. %AGV %A. %Remoci on DQOS 0 73.73 73.73 0 1.56 87.91 89.48 7.85 8.76 26.51 35.28 21.52 16.63 17.76 34.39 42.18 17.66 12.77 30.44 48.05 [mL CH4 / gDQOadicionada ] 0 71.44 103.7 [mL CH4 / gDQOadicionada ] 0 8.16 45.7 86.7 92.08 80.8 79.02 68.78 47.71 39.59 80.8 80.4 77.41 63.26 58.23 %AGV %A. DQO del residuo %Remoci on [mL CH4 / DQOS gDQOadicionada ] 0 0 12.2 7.1 14.82 45 27.78 81.1 56.38 97.3 %Biodegradabilidad anaer obica 0 4.79 30.4 54.8 65.7 %Biodegradabilidad anaer obica 0 48.2 70.0 %Biodegradabilidad anaer obica 0 5.5 30.8 58.5 62.2
Tiempo (d as) 0 10 20 30 40 DQOS [mg O2 /L] 3751 2155 1826 DQOS [mg O2 /L] 5244 4832 4115 3032 2724
Balance nal de DQO neta DQO-CH4 DQO-AGV DQOS [mg O2 /L] [mg O2 /L] [mg O2 /L] 0 2251.3 2786 38.4 2201.6 2446 240.5 1916.3 2373 433.15 1329.4 2012 519.36 1103.2 1215
Tiempo (d as) 0 20 40
DQO-CH4 [mg O2 /L] 0 514.36 746.88
DQO-AGV [mg O2 /L] 2667.5 1956.8 925.1
Tiempo (d as) 0 10 20 30 40
DQO-CH4 [mg O2 /L] 0 82.17 459.84 872.25 926.4
DQO-AGV [mg O2 /L] 3866.8 4610.5 1390.3 931.5 670
7.3 Comparaci on de la simulaci on num erica y los resultados experimentales
91
Composici on aproximada del biog as generado

Las gr acas mostradas en la gura 7-3-(d, e, f), presentan la composici on aproximada del biog as generado durante el desarrollo de las pruebas de biodegradabilidad anaer obica de la mezcla de residuos org anicos bajo estudio. Las gr acas muestran la existencia de una etapa inicial de aproximadamente tres d as para las cargas (1:1) y (2:1) en la cual la producci on de biog as es muy baja, de forma que la fracci on molar de los componentes del biog as no puede estimarse. Las gr acas fueron desarrolladas mediante la suposici on de que el biog as generado durante los ensayos de biodegradabilidad anaer obica corresponde a una mezcla gaseosa con solo dos componentes: metano y di oxido de carbono. Aunque las gr acas para la composici on del biog as generado por cada raz on de carga org anica exhiben tendencias similares, existen diferencias apenas apreciables en los porcentajes m aximos medidos para el metano generado, de forma que la razones de carga (1:1 y 3:1) presentan el mayor contenido de metano (65 % aproximadamente) en un tiempo pr oximo a la mitad del ensayo, mientras que el contenido m aximo que alcanz o la raz on de carga (2:1) fue de (60 %) para un tiempo de 25 d as despu es del inicio del ensayo.
7.3.
Comparaci on de la simulaci on num erica y los resultados experimentales
Con base en las tres simulaciones desarrolladas, una para cada raz on de carga org anica, se generaron las gr acas presentadas en la gura 7-4 para los par ametros de pH y contenido de acidos grasos vol atiles, en las cuales se realiza la comparaci on de las curvas establecidas mediante el conjunto de mediciones experimentales y las curvas obtenidas a trav es de la simulaci on num erica del modelo ADM-1. De acuerdo a la tendencia de las curvas simuladas en la gr aca de la parte (a), el pH se mantiene en un rango entre (pH 7.1 y 7.9). Sin embargo esta tendencia es diferente a la que presentan las curvas obtenidas mediante la medici on del pH durante el desarrollo experimental. La raz on de estas diferencias podr a atribuirse al efecto que ejercen sobre el pH las concentraciones de solutos contenidos en la soluci on del medio Balch y la soluci on reductora que se adiciona durante la instalaci on de los ensayos, las cuales no son consideradas durante el desarrollo de la simulaci on. Por su parte el modelo realiza una estimaci on del valor del pH a partir de una ecuaci on algebraica para el balance de cargas la cual es una aproximaci on al complejo proceso de equilibrio qu mico que se suscita durante el desarrollo del proceso de digesti on anaer obica. No obstante la diferencia m axima entre los valores de pH para las curvas experimentales y simuladas es de aproximadamente una unidad de pH. En la gura 7-4-(b,c,d) se presentan las gr acas que comparan las curvas experimentales y simuladas para el contenido de acidos grasos vol atiles totalizados como contenidos de acido ac etico para cada una de las razones de carga.
92
3000 Demanda Qu mica de Ox geno[mg O2 /L]
2500
(DQO-AGV) (DQO-CH4)
Demanda Qu mica de Ox geno[mg O2 /L]
(DQOS )
4000 3500 3000 2500 2000 1500 1000 500 0 0 5 10 15 20 25 Tiempo [d as]
(DQOS ) (DQO-AGV) (DQO-CH4 )
2000
1500
1000
500
0 0 5 10 15 20 25 Tiempo [d as] 30 35 40
30
35
40
(a)
6000 Demanda Qu mica de Ox geno[mg O2 /L] 5500 5000 4500 4000 3500 3000 2500 2000 1500 1000 500 0 0 5 10 15 1 20 25 30 35 40 Tiempo [d as] (DQOS ) (DQO-AGV) (DQO-CH4 )
(b)
Fracci on molar de metano y di oxido de carbono
1 .0 0 .9 0 .8 0 .7 0 .6 0 .5 0 .4 0 .3 0 .2 0 .1 0 0 5 10 15 20 25 30
CH4 (1:1) CO2 (1:1)
35
40
Tiempo [d as]
(c)
Fracci on molar de metano y di oxido de carbono 1 .0 0 .9 0 .8 0 .7 0 .6 0 .5 0 .4 0 .3 0 .2 0 .1 0 0 5 10 15 20 25 1 Tiempo [d as] 30 35 40 Fracci on molar de metano y di oxido de carbono CH4 (2:1) CO2 (2:1) 1 .0 0 .9 0 .8 0 .7 0 .6 0 .5 0 .4 0 .3 0 .2 0 .1 0 0 5 10
(d)
CH4 (3:1) CO2 (3:1)
15
20 25 1 Tiempo [d as]
30
35
40
(e)
(f)
Figura 7-3.: (a, b, c). Balance de la DQOS para cada una de las razones de carga org anica evaluadas. (d, e, f). Variaci on de la concentraci on de metano y di oxido de carbono presente en el biog as obtenido en los reactores de tratamiento durante los ensayos de biodegradabilidad anaer obica.
93
En cada una de las gr acas mencionadas se puede observar la escasa similitud entre las tendencias de las curvas experimentales y simuladas. Este resultado puede atribuirse al conjunto de valores de los par ametros cin eticos utilizados en el desarrollo de la simulaci on. No obstante en algunos intervalos de la simulaci on los valores de estas concentraciones se encuentran en un rango relativamente cercano. En la gura 7-5 se presentan las gr acas de comparaci on entre las curvas experimentales y simuladas para la alcalinidad (a,b,c) y la demanda qu mica de ox geno soluble (d). En las gr acas para la alcalinidad puede observarse alguna similitud entre estos dos tipos de curvas para cada una de las razones de carga. Aunque en cada gr aca las curvas se encuentran en un rango cercano, existe una ligera discrepancia adjudicada a que en el desarrollo del modelamiento de la alcalinidad se requiere relacionar los contenidos de cada uno de los acidos y bases en t erminos de DQO/L con una concentraci on total en t erminos de mg CaCO3 /L, de forma que deben utilizarse las relaciones expresadas en la tabla 6-3, las cuales est an basadas en estimaciones te oricas y no tienen en cuenta los equilibrios complejos que se desarrollan en la interfaz l quido-gas del medio anaer obico. En esta misma gura en la gr aca (d) puede observarse como la curva obtenida de cada una de las implementaciones del modelo tiende a representar el comportamiento observado mediante las curvas experimentales para el par ametro demanda qu mica de ox geno soluble. No obstante que se presenta cierta similitud en alguna de las curvas como la que representa la raz on de carga (1:1), existen diferencias apreciables entre los dos tipos de curvas durante el periodo considerado para la simulaci on. Esta divergencia en los resultados es consecuencia del efecto que ejercen algunos compuestos org anicos presentes en el residuo los cuales afectan la medici on de la DQO, y no son tomados en consideraci on durante el desarrollo del modelamiento. En la gura 7-6-(a) se muestra la gr aca en la cual se comparan la curvas simuladas y las curvas establecidas mediante la medici on de gas metano generado durante el desarrollo experimental. En la gr aca se observan diferencias entre ambos tipos de curvas; especialmente en la velocidad de generaci on del metano para cada raz on de carga y la duraci on de la fase inicial de adaptaci on. Sin embargo en esta gr aca tambi en se distingue la relativa convergencia que presentan ambos tipos de curvas espec camente en el valor del volumen m aximo generado de metano. En la misma gura tambi en se observan las gr acas (b,c,d) en las cuales se presentan en contraste las curvas obtenidas para la composici on aproximada del biog as generado, tal como fueron medidas durante el desarrollo de los ensayos de biodegradabilidad anaer obica y las curvas generadas mediante la simulaci on. Con base en estas gr acas se puede establecer la relativamente alta precisi on que presenta el modelo al predecir muy cercanamente la composici on porcentual del biog as generado, espec camente para la raz on de carga (1:1). No obstante para las otras razones de carga se observa que la diferencia m axima para las fracciones molares de cada componente entre los dos tipos de curvas es de alrededor del 10 %.
94
9.00 8.75 8.50 8.25 pH [Unidades de pH] 8.00
5000
(1:1-E) (2:1-E) (2:1-S) (3:1-E) (3:1-S)

Acidos Grasos Vol atiles [mg CH3 COOH/L]
(1:1-E) (1:1-S)
(1:1-S)
4500 4000 3500 3000 2500 2000 1500 1000 500 0 0 5 10 15 20 25 Tiempo [d as] 30
7.75
7.50
7.25
7.00
6.75 6.50 6.25 6 0 5 10 15 20 25 Tiempo [d as] 30 35 40
35
40
(a)
5000 Acidos Grasos Vol atiles [mg CH3 COOH/L] 4500 4000 3500 3000 2500 2000 1500 1000 500 0 0 5 10 15 20 25 Tiempo [d as] (2:1-E) (2:1-S) 5000 Acidos Grasos Vol atiles [mg CH3 COOH/L] 4500 4000 3500 3000 2500 2000 1500 1000 500 0 30 35 40 0 5 10
(b)
(3:1-E) (3:1-S)
30
35
40
(c)
(d)
Figura 7-4.: Comparaci on de las curvas generadas mediante la simulaci on num erica y las curvas obtenidas mediante la medici on de la concentraci on del pH (a) y del contenido de acidos grasos vol atiles (b,c,d) para cada una de las razones de carga org anica durante los ensayos de biodegradabilidad anaer obica. [E=Experimental; S=Simulaci on].
95
4000 3500 Alcalinidad [mg CaCO3 /L] 3000 2500 2000 1500 1000 500 0 0 5 10 15 20 25 Tiempo [d as] 30 35 40 (1:1-E) (1:1-S) Alcalinidad [mg CaCO3 /L]
4000 3500 3000 2500 2000 1500 1000 500 0 0 5 10 15 20 25 Tiempo [d as] 30 35 40 (2:1-E) (2:1-S)
(a)
6000 5500 Demanda Qu mica de Ox geno [mg O2 /L] 5000 4500 4000 3500 3000 2500 2000 1500 1000 (3:1-E) (3:1-S) 0 5 10 15 20 25 Tiempo [d as] 30 35 40 500 0 0 5 10
(b)
(1:1-E) (1:1-S) (2:1-E) (2:1-S) (3:1-E) (3:1-S)
6000 5500 5000 Alcalinidad [mg CaCO3 /L] 4500 4000 3500 3000 2500 2000 1500 1000 500 0 1
30
35
40
(c)
(d)
Figura 7-5.: Comparaci on de las curvas generadas mediante la simulaci on num erica y las curvas obtenidas mediante la medici on de la alcalinidad (a,b,c) y la DQOS para cada una de las razones de carga org anica durante los ensayos de biodegradabilidad anaer obica. [E=Experimental; S=Simulaci on].
96
1000 900 CH4 (1:1-E) Producci on de metano [mg DQO-CH4 /L] 800 700 600 500 400 300 200 100 0 0 5 10 15 20 25 Tiempo [d as] 30 (1:1-E) (1:1-S) (2:1-E) (2:1-S) (3:1-E) (3:1-S) 35 40 CH4 (1:1-S) CO2 (1:1-E) Fracci on molar de metano y di oxido de carbono 1 .0 0 .9 0 .8 0 .7 0 .6 0 .5 0 .4 0 .3 0 .2 0 .1 0 0 5 10 15 20 25 30 35 40 Tiempo [d as] CO2 (1:1-S)
(a)
CH4 (2:1-E) CH4 (2:1-S) CO2 (2:1-E) Fracci on molar de metano y di oxido de carbono 1 .0 0 .9 0 .8 0 .7 0 .6 0 .5 0 .4 0 .3 0 .2 0 .1 0 0 5 10 15 20 25 30 35 40 Tiempo [d as] 1 Fracci on molar de metano y di oxido de carbono CO2 (2:1-S)
(b)
CH4 (3:1-E) CH4 (3:1-S) CO2 (3:1-E) CO2 (3:1-S) 0 .9 0 .8 0 .7 0 .6 0 .5 0 .4 0 .3 0 .2 0 .1 0 0 5 10 15 20 25 30 35 40 Tiempo [d as] 1
(c)
(d)
Figura 7-6.: Comparaci on de las curvas generadas mediante la simulaci on num erica y las curvas obtenidas mediante la medici on de la producci on de metano (a) y la composici on aproximada del biog as generado (b,c,d) para cada una de las razones de carga org anica durante los ensayos de biodegradabilidad anaer obica. [E=Experimental; S=Simulaci on].
8. Estimaci on del potencial de generaci on energ etico

En la gran mayor a de las plantas de tratamiento de residuos org anicos provenientes del desarrollo de las actividades agropecuarias se utiliza el esti ercol l quido como material base para la generaci on de biog as, este residuo es mezclado con otros tipos diferentes de residuos org anicos biodegradables hasta formar una pasta homog enea la cual es introducida al interior de los reactores de tratamiento mediante la utilizaci on de bombas hidr aulicas y equipo auxiliar. Para lograr un desarrollo eciente del proceso de digesti on anaer obica se requiere realizar operaciones tales como calentamiento y mezclado del sustrato contenido en el reactor. De esta forma una parte de la energ a obtenida mediante la valorizaci on energ etica del biog as generado debe utilizarse durante el desarrollo continuo del proceso de digesti on anaer obica. Con el n de realizar el planteamiento de un modelo que permita establecer el balance energ etico de una instalaci on a escala piloto para la generaci on de biog as a partir de la codigesti on anaer obica de residuos c tricos y esti ercol bovino, se utilizan los valores reportados en la literatura para los consumos energ eticos de los equipos auxiliares utilizados en instalaciones similares para el desarrollo del proceso [81]. Mediante el modelo de balance energ etico se establecen las formas en las que se desarrolla la producci on y consumo neto de la energ a del sistema. La gura 8-1 muestra el diagrama de procesos de un sistema a escala piloto del cual se han extra do los par ametros de funcionamiento para cada equipo, mientras que la gura 8-2 muestra los ujos energ eticos de entrada y salida de la planta piloto.
8.1.
Planteamiento del modelo energ etico
La ecuaci on (8-1) representa el cambio en funci on del tiempo de la energ a neta que se genera en el sistema para la producci on de biog as. En esta ecuaci on este cambio de la energ a neta es igual a la suma algebraica de t erminos que representan la producci on y consumo de energ a durante el desarrollo del proceso de digesti on anaer obica [81]. dEneta prod E perd E perd ) + (E prod perd perd prod = (E term Econd Ecalen + Ecalenmicr )[kWh/dia] agit elect bomb dt
(8-1)
prod es la Donde Eneta es la producci on neta de energ a del reactor en t erminos de [kWh], E elect producci on de energ a el ectrica mediante un sistema CHP (Generador combinado de potencia y calor, por sus siglas en ingl es), a partir de la combusti on del gas metano presente en el
98
8 Estimaci on del potencial de generaci on energ etico
Figura 8-1.: Esquema general de funcionamiento de la planta piloto que se utilizo como referencia para establecer el modelo energ etico [40].
Ecalentamiento microbial Ecalentamiento Ebombeo Eagitacin ETransferencia de calor Eelctrica Etermica
Figura 8-2.: Esquema general para los ujos energ eticos de entrada y salida de la planta piloto que se utilizo como referencia para establecer el modelo energ etico [40].
perd representa el requerimiento energ biog as generado, E etico del sistema de bombeo del bomb perd sustrato hacia el interior del reactor, Eagit es la energ a que consume el sistema de agiprod taci on del sustrato dentro del reactor, Eterm es la energ a t ermica producida mediante el perd aprovechamiento del biog as generado. Econd es la perdida energ etica debido a la transferencia perd de calor por conducci on y convecci on desde el reactor hacia el ambiente circundante. E calen corresponde a la energ a requerida para el calentamiento del sustrato durante el desarrollo prod del proceso de digesti on anaer obica y E a interna que se genera calen-micr representa la energ
8.1 Planteamiento del modelo energ etico
99
mediante el desarrollo de los procesos bioqu micos. A continuaci on se muestran las ecuaciones para el c alculo de los t erminos correspondientes a los consumos energ eticos en la instalaci on a escala piloto.
prod = QG PCH4 HC elect [kWh/dia] E elect
(8-2)
prod E term = QG PCH4 HC term [kWh/dia]
(8-3)
Donde QG es la producci on de biog as [Nm3 /d a], PCH4 es la fracci on molar de metano contenido en el biog as [ %], HC es el poder calor co inferior del gas metano [kWh/Nm3 ], elect y term , son la eciencia de generacion el ectrica y la eciencia t ermica del sistema utilizado para la generaci on de potencia y calor. El consumo energ etico de la operaci on de bombeo se calcula mediante la siguiente expresi on: 1 bombatornillo
perd = Qin (H)(g)(tp ) E bomb
[kWh/dia]
(8-4)
Donde Qin es el caudal de residuo que ingresa al tanque reactor [m3 /s], H es la altura a la cual debe elevarse el residuo para el ingreso al reactor, es la densidad del residuo org anico que esta siendo bombeado [kg/m3 ], g es la aceleraci on de la gravedad [m/s2 ], tp es el tiempo requerido para el bombeo diario del residuo hacia el interior del reactor [h/d a] y bomba-tornillo es la eciencia mec anica del equipo de bombeo utilizado. El termino de perdida de energ a debido al consumo de potencia por parte del sistema de agitaci on puede ser calculado mediante la siguiente relaci on:
perd = Vliq (S)(ts )[kWh/dia] E agit
(8-5)
Donde Vliq es el volumen l quido del tanque reactor [m3 ], S es la potencia especica de consumo del sistema mec anico de agitaci on [kW/m3 ] y ts es el tiempo requerido para realizar la agitaci on completa del sustrato contenido en el reactor [h/d a]. La energ a requerida para realizar la agitaci on del sustrato es funci on de la geometr a y el volumen del reactor, as como de la viscosidad del sustrato. Las perdidas energ eticas causadas por la transferencia de calor por conducci on y convecci on desde el reactor hacia el ambiente puede ser estimada mediante las expresiones mostradas en la tabla 8-1. Donde RTotal es la resistencia t ermica total a la transferencia de calor, L1 y L2 corresponden al espesor de la pared de acero inoxidable de tanque reactor y el espesor del aislante t ermico (lana mineral) del reactor respectivamente; k1 y k2 corresponden a las
100
Tabla 8-1.: Ecuaciones transferencia de calor por conducci on y convecci on. perd = (Treactor Tambiente ) E cond RT otal RT otal = 1 ln (R2 /R1 ) ln (R3 /R2 ) 1 + + + h1 A1 2K1 L1 2K2 L2 h2 A2
conductividades t ermicas del acero inoxidable y el aislante del reactor. R1 , R2 y R3 corresponden al radio interno del reactor, radio de contacto entre la pared del reactor y el aislante, y el radio externo de la capa cil ndrica de aislante t ermico respectivamente. Los factores h1 y h2 corresponden a los coecientes de transferencia de calor desde el reactor para las partes interna y externa respectivamente. La gura 8-3-(b) muestra la geometr a del reactor. La energ a requerida para desarrollar el calentamiento del sustrato contenido en el reactor puede ser estimada mediante la siguiente relaci on:
perd = Qin Cp [(Treact Tsustrato )( 1 )[kWh/dia] E calen 3,6
(8-6)
Donde Qin es el caudal de residuo que ingresa al tanque reactor [m3 /s], Cp es el calor especico del sustrato [kJ/(kg K)], Treact es la temperatura del sustrato contenido en el reactor [K], Tsustrato es la temperatura del sustrato fresco que va ingresar al reactor [K]. El calor generado por las reacciones bioqu micas que se desarrollan durante el proceso de digesti on anaer obica puede ser estimado mediante la siguiente relaci on:
1 prod a] E (Ej fj j Vliq ( )) [kWh/d calenmicr = 3,6 j =512 (8-7) Donde Vliq es el volumen l quido del tanque reactor [L], Ej es la energ a interna liberada hacia el ambiente debido al proceso (j) de actividad microbiana [kJ/mol], fj es la masa molar por gDQO de metabolito producido mediante los procesos bioqu micos que se suscitan durante la digesti on anaer obica [mol/g DQO], j es la raz on cin etica del proceso bioqu mico (j) [g DQO/(L d a)].
8.2 Desarrollo del modelo
101
8.2.
Desarrollo del modelo
El modelo energ etico fue implementado como un subsistema de ecuaciones en la plataforma utilizada para la simulaci on del modelo ADM-1; debido a que algunos de los t erminos que conforman la ecuaci on de balance energ etico requieren la interacci on con las variables que se eval uan dentro del modelo del proceso de digesti on anaer obica. Para la implementaci on de la ecuaci on de balance de energ a para el sistema compuesto por el tanque reactor con agitaci on semicontinua y los sistemas auxiliares de calentamiento, alimentaci on y generaci on combinada de potencia y calor, se requiere establecer puntualmente los valores de cada una de las variables que est an presentes en los t erminos de la ecuaci on de balance, los cuales son mostrados en la tabla 8-2. Los t erminos que determinan la producci on de energ a el ectrica y energ a t ermica a partir de la combusti on del metano generado se calculan tomando como referencia el valor del caudal gas ) y la fracci de biog as (V on de metano (XCH4 ) presente en el biog as, las cuales se estiman mediante la implementaci on desarrollada del modelo ADM-1. Para los valores poder calor 3 co inferior del metano se utilizo el valor (10.26 kWh/m ). De acuerdo con Lubken et al. [81], las unidades CHP generalmente presentan una eciencia de conversi on de energ a mecanica/el ectrica de alrededor de 30 % (elect =0.3) y la eciencia de generaci on de energ a t ermica es de aproximadamente 50 % (term =0.5); por lo tanto se estima que aproximadamente el 20 % de la energ a contenida en el gas metano se disipa en forma de irreversibilidades del sistema electromec anico CHP. Estas eciencias de conversi on se muestran en el diagrama de Sankey gura 8-3-(b) desarrollado para la distribuci on de energ a a trav es del sistema CHP. En el mismo diagrama se observan los porcentajes relativos de los consumos energ eticos requeridos para el funcionamiento de la instalaci on, estos fueron calculados con base en la fracci on que representa cada termino respecto a la energ a el ectrica total y la energ a t ermica total generadas por el sistema CHP. Para el termino potencia consumida durante la alimentaci on del sustrato al reactor mediante una bomba hidr aulica, se utilizaron los valores referidos en la tabla (8-2) para el caudal de ingreso al reactor, altura del punto de ingreso del sustrato en el reactor, la densidad del sustrato medida durante la etapa experimental. El tiempo promedio de bombeo y la eciencia de la bomba tornillo se estimaron en 15 min/d a y (Bombatornillo =0,5), respectivamente [81]. El termino energ a consumida durante la agitaci on del sustrato fue establecido con base en el volumen l quido del reactor establecido en la tabla 8-2 y los par ametros reportados para la energ a especica de agitaci on y el tiempo promedio diario que dura la agitaci on 3 (0.005 kW/m y 8h/d a respectivamente). El termino energ a requerida para el calentamiento del sustrato se estableci o mediante los par ametros de caudal de ingreso al reactor y el calor especico del residuo obtenida de una investigaci on desarrollada sobre un sustrato
102
similar [37]. Estos par ametros se muestran en la tabla 8-2. Se estimo adicionalmente una temperatura ambiental promedio de 18o C para las condiciones locales y una temperatura del reactor de 30o C de acuerdo a las condiciones en las cuales se desarrollo el ensayo de biodegradabilidad anaer obica. El termino correspondiente a la energ a perdida debido al calor transferido por conducci on y convecci on desde el reactor se c alculo utilizando los valores 17 W/mK y 0.55 W/mK como los coecientes de conductividad t ermica para los materiales acero inoxidable y lana mineral respectivamente. Estos materiales se asumieron como capas cil ndricas con un espesores de 0.0127 m (1/2) para la pared de acero inoxidable y 0.1 m para la lana mineral. De igual forma se utilizaron los valores 1 W/m2 K y 4 W/m2 K para los coecientes de transferencia de calor desde el sustrato hacia la pared interna del reactor y entre la pared externa del reactor al aire exterior [18]. Como resultado de la implementaci on del modelo energ etico se obtuvo la gr aca 8-3-(a) en la cual se muestra la evoluci on de la variable que representa la energ a neta generada y raz on de cambio de la energ a generada por el sistema de tratamiento a escala piloto. A partir de estas gr acas se puede inferir la notable capacidad de generaci on de energ a que presenta este sistema de tratamiento de residuos org anicos biodegradables. De acuerdo a estos resultados, en una instalaci on a escala piloto con el funcionamiento descrito anteriormente, el sistema genera aproximadamente 200 KWh/mes, la cual corresponde a la energ a promedio de consumo mensual de una vivienda en el area urbana, y aproximadamente 4 viviendas en el sector rural en Colombia [31]. Tabla 8-2.: Caracter sticas del sistema a escala piloto.
Variable Volumen total del reactor Volumen l quido del reactor Tiempo de retenci on hidr aulica Caudal de ingreso al reactor Altura del reactor Area efectiva del reactor Densidad del sustrato Calor especico del sustrato S mbolo VT VL HRT in [VL /HRT] V H A Cp
3
Valor 1.89 m (500 Galones) 1.51 m3 (400 Galones) 40 dias 37.85 L/dia 2.7 m 1.97 m2 1145 Kg/m3 4.13 kJ/kg
8.2 Desarrollo del modelo
103
315 Energ a neta y variaci on de la energ a neta generada [kWh, KWh/d a] 295 275 255 235 215 195 175 155 135 115 95 75 55 35 15 5 0 5 10 15 20 25 Tiempo [d as] 30 35 40 o Eneta neta E
(a)
(b)
(c)
Figura 8-3.: (a). Curva obtenida mediante la simulaci on num erica para la energ a neta y raz on de cambio de la energ a neta generada. (b). Geometr a del tanque reactor simulado. (c). Diagrama de Sankey para la energ a del sistema a escala piloto.
9. Conclusiones
En la actualidad existen muchos residuos org anicos biodegradables como los desechos c tricos y el esti ercol bovino, sobre los cuales no se realiza un tratamiento eciente y terminan dispuestos en rellenos sanitarios o en las vertientes h dricas. La digesti on anaer obica es un proceso biotecnol ogico que presenta potenciales alternativas para la soluci on de aspectos relevantes de la conservaci on medio ambiental como el tratamiento efectivo de los residuos biodegradables y la generaci on de energ a renovable. El primer aspecto se evidencia en la recuperaci on de los nutrientes contenidos en el residuos, los cuales pueden nuevamente aprovecharse mediante la aplicaci on controlada del fango euente del sistema de tratamiento anaer obico, sobre los terrenos cultivables. Por su parte el segundo aspecto subyace de la valorizaci on energ etica de la fracci on de metano contenida en el biog as generado, a trav es de sistemas de transformaci on energ etica como los sistemas de generaci on combinada de potencia y calor (CHP). Los ensayos de biodegradabilidad anaer obica est an basados en metodolog as estandarizadas internacionalmente. Estos procedimientos permiten establecer en un tiempo relativamente corto los par ametros de funcionamiento del tratamiento anaer obico sobre un sustrato en particular. A trav es de estos ensayos se calculan ndices de desempe no como el volumen m aximo de generaci on de metano, la tasa m axima de generaci on de metano, el porcentaje m aximo de remoci on de la DQO, y la tasa m axima de remoci on de la DQO, entre otros. Estos ndices son necesarios en las etapas de dise no, puesta a punto, operaci on y control de las instalaciones a gran escala para el tratamiento de residuos org anicos biodegradables. Con base en los resultados obtenidos en este estudio mediante los ensayos de biodegradabilidad anaer obica se puede inferir que es factible la producci on de biog as con un alto contenido de metano a partir de la codigesti on anaer obica de los residuos c tricos y el esti ercol bovino. El modelamiento del proceso de digesti on anaer obica permite establecer el comportamiento del sistema de tratamiento de forma que puedan predecirse los resultados que se obtendr an de un sistema a escala piloto para el tratamiento de un sustrato en particular. Dentro del conjunto de modelos que han sido publicados, el modelo ADM-1, propuesto por Batstone et al. [14], ha demostrado ser una herramienta ecaz para la simulaci on del proceso de digesti on anaer obica, estableciendo un balance positivo entre la complejidad y la aplicabilidad. No obstante la implementaci on de este modelo requiere la determinaci on de varias decenas de par ametros espec cos y la composici on qu mica del sustrato bajo an alisis.
105 En este estudio se encontr o que la implementaci on del modelo ADM-1 utilizando los valores de los par ametros reportados en el caso base y las concentraciones de los metabolitos estimadas mediante las metodolog as que relacionan las mediciones experimentales con el contenido de estos compuestos en el sustrato, present o un ajuste cercano entre los datos experimentales y los datos obtenidos a trav es de la simulaci on. Sin embargo se evidencian algunas restricciones que se generan sobre el modelo cuando se ha realizado una limitada caracterizaci on del sustrato. Debido a la complejidad de las t ecnicas de medici on necesarias para establecer la composici on qu mica del residuo, es importante el desarrollo e implementaci on de algunas metodolog as [144, 147, 148] que intentan relacionar mediciones sicoqu micas que se realizan rutinariamente en el an alisis de aguas residuales y desechos org anicos con los resultados obtenidos de la caracterizaci on exhaustiva del sustrato. Con base en los resultados obtenidos a partir de las simulaciones establecidas, se observa que el modelo intenta reproducir la tendencia general de cada una de las curvas establecidas mediante las mediciones experimentales. Sin embargo existen algunas discrepancias entre los valores representados por cada tipo de curva debido a que el modelo presenta importantes simplicaciones para la descripci on de cada uno de los subprocesos m as relevantes que conforman la digesti on anaer obica; de igual manera suprime la inuencia que ejercen diversos factores como el potencial de reducci on-oxidaci on de cada reacci on bioqu mica. Por consiguiente es importante establecer la necesidad de realizar una mayor caracterizaci on del sustrato y una serie de ensayos experimentales m as detallados con el n de calibrar el modelo y obtener una mayor precisi on en los resultados de la simulaci on para un sustrato especico. De acuerdo al modelo energ etico desarrollado en conjunto con la implementaci on del modelo ADM-1, se puede establecer el potencial que presenta esta tecnolog a de tratamiento de los residuos org anicos biodegradables, como fuente de energ a renovable. Como se observo en los resultados obtenidos mediante este modelo, durante los primeros estadios del desarrollo del proceso el sistema presenta un balance energ eticamente negativo, sin embargo a medida que se desarrolla el proceso de digesti on anaer obica del sustrato y se inicia la producci on del gas metano, se obtiene una generaci on neta de energ a que supera ampliamente el requerimiento energ etico del sistema a escala piloto que fue modelado. El aporte de esta investigaci on se basa en la determinaci on de la biodegradabilidad anaer obica de la mezcla de residuos, de acuerdo con un protocolo estandarizado, con lo cual se demuestra la factibilidad que presenta la tecnolog a de digesti on anaer obica para la generaci on de energ a renovable a partir de estos residuos org anicos. Este estudio es una primera aproximaci on al desarrollo de una estrategia optima tanto para el tratamiento de residuos org anicos biodegradables de origen agroindustrial, como para la obtenci on biocombustibles como bioetanol y biog as a partir de procesos biotecnol ogicos como la fermentaci on alcoh olica y la digesti on anaer obica.
A. Anexo: Preparaci on de las soluciones qu micas utilizadas

Soluci on stock de minerales
Contenido de soluto por cada litro de la soluci on correspondiente. Soluci on 1: 27.2 g KH2 PO4 Soluci on 2: 36.5 g Na2 HPO4 2H2 O Soluci on 3A: 24 g NH4 Cl; 24 g NaCl; 8 g MgCl2 6H2 O Soluci on 3B: 11 g CaCl2 2H2 O*(Agrafar y guardar bajo refrigeraci on (4o C)) Soluci on 4: 80 g Na2 HCO3 Soluci on 5: 240.2 g Na2 S9H2 O*(Preparar con agua destilada realizando intercambio de gases con N2 , proteger de la luz directa y guardar en refrigeraci on (4o C)) Soluci on 6: 0.5 g Resarzurina*(Proteger de la luz directa y guardar en refrigeraci on o (4 C))
Soluci on stock acida (I)

Contenido de soluto por cada litro de la soluci on. Tabla A-1.: Concentraci on de las sustancias que componen la soluci on acida I.
Concentraci on 50 mM HCL 1 mM H3 BO3 0.5 mM MnCl2 7.5 mM FeCl2 0.5 mM CoCl2 0.1 mM NiCl2 0.5 mM ZnCl2 masa de soluto 1.80 g HCL 61.80 mg H3 BO3 61.25 mg MnCl2 943.50 mg FeCl2 64.50 mg CoCl2 12.86 mg NiCl2 67.70 mg FeCl2
107
Soluci on stock alcalina (II)

Contenido de soluto por cada litro de la soluci on. Tabla A-2.: Concentraci on de las sustancias que componen la soluci on alcalina II.
Concentraci on 10 mM NaOH 0.1 mM Na2 SeO3 0.1 mM Na2 WO4 0.1 mM Na2 MoO4 masa de soluto 400 mg NaOH 17.3 mg Na2 SeO3 29.4 mg Na2 WO4 20.5 mg Na2 MoO4
Medio mineral
Para la preparaci on de un litro de soluci on de medio mineral se agregan: 15 mL Soluci on 1 15 mL Soluci on 2 12.5 mL Soluci on 3A 1 mL Soluci on 5 1 mL Soluci on Stock Acida (I) 1 mL Soluci on Stock Alcalina (II) Aforar hasta un litro con Agua destilada Calentar hasta ebullici on para liberar el ox geno disuelto (No hervir por mas de 20 segundos para evitar la evaporaci on excesiva) Enfriar bajo ujo o burbujeo de nitr ogeno para remover el ox geno desde el medio.
Soluci on est eril-reductora

En una botella de 160 mL se adicionan contenidos de las soluciones 4 y 5 de forma que la relaci on de vol umenes sea 50:1, teniendo en cuneta el volumen total requerido para el montaje de todos los reactores que se incluyen en el montaje de la prueba de biodegradabilidad anaer obica: 50 mL de la soluci on (4).
108 1 mL de la soluci on (5).
A Anexo: Preparaci on de las soluciones qu micas utilizadas
Esta botella se sella mediante agrafes y se realiza intercambio de gases con nitr ogeno o y se guarda en refrigeraci on (4 C) para evitar la oxidaci on, la cual se manifestara por oscurecimiento de la soluci on o por la formaci on de precipitados.
Soluci on de hidr oxido de sodio NaOH al 5 %-[m/v]

Para cada 100 mL de soluci on se adicionan 5 g de NaOH, por consiguiente para cada litro de soluci on se requieren 50 g de NaOH. Se agrega un volumen conocido de agua destilada en un vaso de precipitado, se agrega la totalidad del soluto correspondiente, esta adici on se realiza bajo agitaci on mediante un agitador magn etico y dentro de la c amara de extracci on debido a que la reacci on es altamente exot ermica. Posteriormente cuando se encuentre diluido se agrega el volumen restante para completar la soluci on, y se agregan algunas gotas de indicador fenolftale na, con lo cual la soluci on adquiere un color fucsia (la fenolftale na indica un pH alcalino); nalmente se dispensan los vol umenes de la soluci on en cada botella utilizada para medici on y se sellan mediante tap on de caucho y agrafe de aluminio.
C alculo de la densidad de la soda

Se requiere conocer el valor de la densidad de la soluci on de NaOH utilizada para la medici on del volumen de metano generado en los reactores anaer obicos, la estimaci on de la densidad de la soluci on se realiza de la siguiente forma; Pesar un bal on aforado de 100 mL (V) Agregar al bal on el volumen de soda hasta que se complete el aforo y registrar este peso (w) Calcular la densidad como =w/V
B. Anexo: An alisis de s olidos

Determinaci on de la concentraci on de s olidos totales, suspendidos y vol atiles
Materiales y equipos
Bomba de vac o y Filtros Milipore Probeta, pinzas, c apsulas de porcelana Desecador Mua/Estufa para desecado Balanza anal tica Ba no Mar a
Procedimiento s olidos totales

1. Desecar las c apsulas durante una hora en la mua a una temperatura de 105o C, dejar enfriar y pesar las c apsulas. (Registrar este valor de peso como P1 ) 2. Colocar en las c apsulas un volumen conocido (V en mL) de la muestra homogeneizada y evaporar en el ba no mar a la humedad. 3. Secar en la mua a una temperatura entre 103o C y 105o C durante una hora. 4. Transferir las c apsulas al desecador hasta que estas se enfr en y posteriormente pesarlas. (Registrar este valor de peso como P2 ). 5. Incinerar en la mua a 550o C +/- 50o C durante una hora. 6. Colocar las c apsulas en el desecador y dejar que se enfr en. 7. Pesar las c apsulas. (Registrar este valor como P3 )
C alculos
ST T (mg/L) = ST V (mg/L) = (P2 P1 )(1X 106 ) V
(P2 P3 )(1X 106 ) V ST F (mg/L) = ST T ST V
110
B Anexo: An alisis de s olidos
Donde P1 =Peso de las c apsulas desecadas (g) P2 =Peso de las c apsulas luego de secado a una temperatura entre 103o C y 105o C (g) P3 =Peso de las c apsulas luego de incineraci on a una temperatura de 550o C +/- 50o C (g) (1X106 )=Factor de Conversi on de gramos a miligramos y de mililitros a Litros
Procedimiento s olidos suspendidos

1. Desecar las ltros durante una hora en la mua a una temperatura de 105o C, Dejar enfriar y pesar los ltros. (Registrar este valor de peso como P1 ) 2. Colocar el ltro en la bomba de vac o, humedecerlo con agua destilada y cerrar el sistema. 3. Filtrar un volumen conocido (V en mL) de la muestra homogeneizada. 4. Secar los ltros en la mua a una temperatura entre 103o C y 105o C durante una hora. 5. Transferir los ltros al desecador hasta que estos se enfr en y posteriormente pesarlos. (Registrar este valor de peso como P2 ). 5. Incinerar en la mua a 550o C +/- 50o C durante una hora. 6. Colocar los ltros en el desecador y dejar que se enfr en. 7. Pesar los ltros. (Registrar este valor como P3 )
C alculos
SST (mg/L) = SSV (mg/L) = (P2 P1 )(1X 106 ) V
(P2 P3 )(1X 106 ) V SSF (mg/L) = SST SSV
Donde P1 =Peso del ltro desecado (g) P2 =Peso del ltro luego de secado a una temperatura entre 103o C y 105o C (g) P3 =Peso del ltro luego de incineraci on a una temperatura de 550o C +/- 50o C (g) (1X106 )=Factor de Conversi on de gramos a miligramos y de mililitros a Litros
C. Anexo: Determinaci on de AGV mediante m etodo colorim etrico

Determinaci on de AGV mediante m etodo colorim etrico
Este procedimiento se basa en una t ecnica descrita por Montgomery et al. [94], seg un la cual los acidos carbox licos (ac etico, but rico, val erico y propi onico) son estericados con etanodiol y los esteres resultantes reaccionan con hidroxilamina para formar acidos. Adicionalmente, la reacci on con cloruro f errico da como resultado la formaci on de hidroximatos ferrosos de color p urpura. La intensidad del color es determinada a trav es del espectrofot ometro y es proporcional al contenido de acidos grasos vol atiles en la muestra. La intensidad de color observada en este an alisis se compara con la intensidad que presenta una soluci on patr on de acido ac etico, por lo tanto la concentraci on calculada a trav es de esta t ecnica se expresa como concentraci on de acido ac etico.
Reactivos
Acido sulf urico al 50 % (V/V):Adicionar con cuidado y con agitaci on constante 50 mL +/- 0.5 mL de acido sulf urico (d2 0 1.84) a 50 mL +/- 0.5 mL de agua y dejar enfriar. Hidr oxido de sodio al 18 % (m/V): Disolver 18 g de NaOH en 80 mL de agua destilada, dejar enfriar y diluir con agua destilada hasta completar 100 mL. Reactivo Acido Etanodiol: Mezclar 15 mL +/- 0.5 mL de etanodiol (etilenglicol) con 2 mL +/- 0.05 mL de acido sulf urico al 50 % (V/V). Sulfato Hidroxiamonio al 10 % (m/V): Disolver 10 g +/- 0.1 g de sulfato hidroxiamonio en 80 mL de agua destilada, diluir posteriormente hasta completar 100 mL. Este reactivo se debe preparar mensualmente y debe guardarse refrigerado. Reactivo de Hidroxilamina: Mezclar 20 mL +/- 0.5 mL de NaOH al 18 % (m/V) con 5 mL +/- 0.1 mL de sulfato hidroxiamonio al 10 % (m/V). Esta soluci on debe prepararse inmediatamente antes de ser utilizada. Reactivo de Cloruro f errico acido: Disolver 10 g +/- 0.05 g de cloruro f errico hexahidratado en 250 mL de agua destilada, agregar 10 mL +/- 0.05 mL de acido sulf urico (d2 0 1.84) y diluir con agua destilada hasta 500 mL (ltrar si es necesario). Guardar refrigerado, y protegido de la luz directa. Esta soluci on tiene un pH de alrededor de 1.6 para evitar la precipitaci on de las sales f erricas.
112
C Anexo: Determinaci on de AGV mediante m etodo colorim etrico
Materiales y Equipos
Espectrofot ometro (Rango visible), celda para espectrofot ometro. Tubos de ensayo, gradilla, pipetas, ltros milipore (o centr fuga), balones aforados de 25 mL y 100 mL.
Soluciones patr on de acido ac etico

Soluci on A: 1 mL de soluci on debe contener 10 mg de acido ac etico Pesar 1 g de acido ac etico glaciar (d2 0 1.05) en un recipiente previamente tarado y transferir cuantitativamente a un bal on aforado de 100 mL. Diluir con agua destilada hasta completar 100 mL. Preparar soluciones frescas cada mes. Agregar en una serie de balones de 100 mL vol umenes de 0, 1, 2.5, 5, 10, 25 y 50 mL de la soluci on A y completar con agua destilada hasta 100 mL. Dispensar estas soluciones en recipientes limpios y secos, estos recipientes deben guardarse en refrigerador. Cada una de estas soluciones contiene 0, 100, 250, 1000, 2500 y 5000 mg/L de acido ac etico respectivamente. Curva de calibraci on de acido ac etico La gura (C-1) muestra una de las curvas de calibraci on que se utilizaron para la determinaci on de la concentraci on de los acidos grasos vol atiles del sustrato contenido en los reactores de tratamiento durante los ensayos de biodegradabilidad anaer obica.
5000 Concentraci on de acido ac etico [mg CH3 COOH/L] 4500 4000 3500 3000 2500 2000 1500 1000 500 0 0 0 .1 0 .2 0 .3 0 .4 0 .5 0 .6 Absorbancia 0 .7 0 .8 0 .9 y = 7164.2x; R2 =0,9992
Figura C-1.: Curva de calibraci on para la concentraci on de acido ac etico mediante la t ecnica emp rica espectrofotom etrica.
113
Procedimiento
Procedimiento Experimental Construcci on de la curva de calibraci on. Se siguen los pasos 2 a 5, utilizando como muestra cada una de las soluciones patr on preparadas con anterioridad. Gracar la Absorbancia contra la concentraci on de acido ac etico en (mg/L) Pipetear 0.5 mL +/- 0.01 mL de la muestra en un tubo de ensayo seco. Agregar 1.7 mL +/- 0.05 mL de acido etanodiol y mezclar con fuerza. Calentar en ba no mar a por 3 min +/- 10 segundos (B) Enfriar inmediatamente el tubo en agua fr a. Adicionar 2.5 mL +/- 0.1 mL de hidroxilamina y mezclar fuertemente, dejar reposar por 10 s (C) Adicionar 10 mL +/- 0.1 mL de cloruro f errico acido en un bal on de 25 mL calibrado. Transferir cuantitativamente la soluci on del tubo de ensayo al bal on. (D) Completar el volumen con agua destilada y agitar vigorosamente. No debe formarse ning un precipitado. Dejar reposar durante 5 minutos. (E) Medir la absorbancia de la soluci on a una longitud de onda de 500 nm usando celdas de 40 mm, esta medici on se debe realizar dentro de la primera hora de realizado el paso 5. (F); (G) Se debe correr un blanco y un patr on con cada lote de muestras usando los mismos reactivos, para realizar estas muestras se deben repetir los pasos 2 a 5, utilizando 0.5 mL +/- 0.01 mL de agua destilada (Blanco) y uno de los patrones en lugar de muestra. El blanco se utiliza para calibrar el 100 % de la transmitancia del equipo. (H) Calcular la concentraci on de AGV de acuerdo con la ecuaci on de la curva de calibraci on obtenida (I) Notas Absorbancia=2-log( %Transmitancia)
B.Evitar el contacto directo del tubo de ensayo con las paredes y resistencias del ba no mar a.
C.El an alisis en tandas (Batch) se facilita con el uso de la gradilla D. Utilizar agua destilada para la remoci on de las ultimas trazas en el tubo de ensayo. E. Esto facilita el escape de los gases disueltos.
F. Tener cuidado de evitar la formaci on de burbujas dentro de la celda. G. Se pueden utilizar celdas de diferente tama no H. El patr on se corre con el objetivo de chequear la curva de calibraci on. La curva de calibraci on es normalmente lineal hasta 10000 mg/L de acido ac etico. La linealidad de la curva es funci on del instrumento utilizado. I. Utilizar la misma celda que se utiliz o para la construcci on de la curva de calibraci on.
Ssu
Saa
Sf a
Sva
va
dSsu Km,su Ssu =Khydch ch +(1-ff a,li )Khydli li su I1 dt Ks,Ssu + Ssu dSaa Km,aa Saa =Khydpr pr aa I1 dt Ks,Saa + Saa dSf a Km,f a Sf a =(ff a,li )Khydli li f a I2 dt Ks,Sf a + Sf a dSva Km,aa Saa Sva =(1-Yaa )(fva,aa ) I1 - KKm,c4+ c4 1 Sbu I2 dt Ks,Saa + Saa aa Sva 1+ S s,Sva Sbu aa I1 - KKm,c4+ c4 S
s,Sbu bu
Sbu
Km,su Ssu Saa dSbu =(1-Ysu )(fbu,su ) I1 +(1-Yaa )(fbu,aa ) KKm,aa dt Ks,Ssu + Ssu su + Saa s,Saa 1 I2 va 1+ S Sbu
Spro
Km,su Ssu Km,aa Saa dSpro =(1-Ysu )(fpro,su ) su I1 +(1-Yaa )(fpro,aa ) dt Ks,Ssu + Ssu Ks,Saa + Saa Km,c4 Sva 1 Km,pro Spro (1-Yc4 )(0.54) c4 I2 pro I2 Sbu Ks,Sva + Sva K 1+ S s,Spro + Spro
aa I1 -
va
Sac
Km,su Ssu Saa Sf a dSac =(1-Ysu )(fac,su ) su I1 +(1-Yaa )(fac,aa ) KKm,aa aa I1 +(1-Yf a )(0.7) KKm,f a dt Ks,Ssu + Ssu s,Saa + Saa s,Sf a + Sf a Km,c4 Sva Km,c4 Sbu 1 (1-Yc4 )(0.31) I2 +(1-Yc4 )(0.8) c4 1Sva I2 + Sbu Ks,Sva + Sva c4 1 + S Ks,Sbu + Sbu 1+ S
va bu
f a I2 +
Sh2
Km,pro Spro Sac (1-Ypro )(0.57) pro I2 - KKm,ac ac I3 Ks,Spro + Spro s,Sac + Sac Km,su Ssu Saa Sf a dSh2 =(1-Ysu )(fh2,su ) su I1 +(1-Yaa )(fh2,aa ) KKm,aa aa I1 +(1-Yf a )(0.3) KKm,f a dt Ks,Ssu + Ssu + S + Sf a s,Saa aa s,Sf a Km,c4 Sva 1 Km,c4 Sbu 1 (1-Yc4 )(0.15) c4 I2 +(1-Yc4 )(0.2) I2 + Sbu va Ks,Sva + Sva Ks,Sbu + Sbu c4 1 + S 1+ S S
va bu
f a I2 +
Sch4
Sic
Km,pro Spro Sh2 (1-Ypro )(0.43) pro I2 - KKm,h2 h2 I1 Ks,Spro + Spro s,Sh2 + Sh2 Km,ac Sac S h2 dSch4 =(1-Yac ) ac I3 +(1-Yh2 ) KKm,h2 h2 I1 dt Ks,Sac + Sac s,Sh2 + Sh2 dSIC =Ci vi,5 Ci vi,6 Ci vi,10 Ci vi,11 Ci vi,12 dt i=19,1124 i=19,1124 i=19,1124 i=19,1124 i=19,1124
Sin
D. Anexo: Ecuaciones modelo ADM-1
Tabla D-1.: Ecuaciones modelo digesti on anaer obica (ADM-1)-especies solubles.
SI
dSIN =-(Ysu )Nbac +Naa -(Yaa )Nbac -(Yf a )Nbac -(2)(Yc4 )Nbac -(Ypro )Nbac -(Yac )Nbac -(Yh2 )Nbac dt dSI =(fSi,Xc )Kdis,Xc dt
ch
pr
li
su
aa
aa
f a
c4
ac
h 2
Tabla D-2.: Ecuaciones modelo digesti on anaer obica (ADM-1)-especies insolubles
dc =-Kdis c +Kdec,Xsu su +Kdec,Xaa aa +Kdec,Xf a f a +Kdec,Xc4 c4 +Kdec,Xpro pro + dt Kdec,Xac ac +Kdec,Xh2 h2 dch =(fch,Xc ) Kdis c -Khyd,ch ch dt dpr =(fpr,Xc ) Kdis c -Khyd,pr pr dt dli =(fli,Xc ) Kdis c -Khyd,li li dt Km,su Ssu dsu =(Ysu ) su I1 -Kdec,Xsu su dt Ks,Ssu + Ssu Km,aa Saa daa =(Yaa ) aa I1 -Kdec,Xaa aa dt Ks,Saa + Saa daa Km,aa Saa =(Yaa ) aa I1 -Kdec,Xaa aa dt Ks,Saa + Saa Km,f a Sf a df a =(Yf a ) f a I1 -Kdec,Xf a f a dt Ks,Sf a + Sf a Km,c4 Sva 1 Km,c4 Sbu 1 dc4 =(Yc4 ) c4 I2 +(Yc4 ) c4 I2 -Kdec,Xc4 c4 Sbu Sva dt Ks,Sva + Sva Ks,Sbu + Sbu 1 + 1+ S Sbu va Km,ac Sac dac =(Yac ) ac I3 -Kdec,Xac ac dt Ks,Sac + Sac Km,h2 Sh2 dh2 =(Yh2 ) h2 I1 -Kdec,Xh2 h2 dt Ks,Sh2 + Sh2 dI =(fXi,Xc ) Kdis,Xc dt
115
116
D Anexo: Ecuaciones modelo ADM-1
Tabla D-3.: Ecuaciones modelo digesti on anaer obica (ADM-1)-Procesos acido-base (DAE)
Numero 1 2 3 4 5 6 7 8
Ecuaci on
+ +SH + -S -2S 2 Scat+ +SN H4 HCO3 CO3
Variable SH + SOH Sva Sbu Spro Sac

SHCO3
SAc Spro Sbu Sva -SOH -San =0 64 112 160 208 KW SH +
(SOH )(Sva )(Sbu )(Spro )-
Ka,va Sva,total =0 Ka,va + SH + Ka,bu Sbu,total =0 Ka,bu + SH +
Ka,pro Spro,total =0 Ka,pro + SH + Ka,ac Sac,total =0 Ka,ac + SH + Ka,CO2 SIC =0 Ka,CO2 + SH + =0
(Sac )-
)(SHCO3
2 )(SCO3
S Ka,HCO3 HCO3 + SH + Ka,HCO3
2 SCO3
9 10 11
)-(S 2 )=0 (SIC )-(SCO2 )-(SHCO3 CO3
SCO2
+ SN H4
+ )(SN H4
SH + SIN =0 Ka,N H4 + SH +
+ )=0 (SIN )-(SN H3 )-(SN H4
SN H 3
Tabla D-4.: Ecuaciones modelo digesti on anaer obica (ADM-1), reacciones acido-base (DE). 5b Sbu Spro
SHCO3
4b Sva -1 -1 -1 -1 -1 -1
6b
j A1 A2 A3 A4 A5 A6
Componente i Procesos Acido-Base Valerato Acido-Base Butirato Acido-Base Propionato Acido-Base Acetato Acido-Base Carbono Inor. Acido-Base Nitr ogeno Inor.
7b Sac
10b
11b SN H3
Razon Cin etica KB,va (Sva (SH + + Ka,va ) Ka,va Sva ) KB,bu (Sbu (SH + + Ka,bu ) Ka,bu Sbu ) KB,pro (Spro (SH + + Ka,pro ) Ka,pro Spro ) KB,ac (Sac (SH + + Ka,ac ) Ka,ac Sac ) (SH + + Ka,CO ) Ka,CO SIC ) KB,CO2 (SHCO3 2 2 KB,N H3 (SN H3 (SH + + Ka,N H3 ) Ka,N H3 SIN )
Tabla D-5.: Ecuaciones complementarias modelo digesti on anaer obica (ADM-1), procesos acido-base. Variable SH +
+ SN H4
Ecuaci on 2 + 4KW SH + =- + 2 2 SAc Spro Sbu Sva + +SH + -S -2S 2 =Scat+ +SN H4 -SOH -San HCO3 CO3 64 112 160 208 + =SIN -SN H SN H4 3
SCO2 =SIC -SCO3
SCO2
117
118
D Anexo: Ecuaciones modelo ADM-1
Tabla D-6.: Ecuaciones para la transferencia de masa en la interfaz l quido-gas. dSgas,h2 Sgas,h2 Qgas Vliq Sgas,h2 =+ (KL a)(Sliq,h2 -16KH,h2 Pgas,h2 ) dt Vgas Vgas dSgas,ch4 Sgas,ch4 Qgas Vliq Sgas,ch4 =+ (KL a)(Sliq,ch4 -64KH,ch4 Pgas,ch4 ) dt Vgas Vgas dSgas,co2 Sgas,co2 Qgas Vliq Sgas,co2 + (KL a)(Sliq,co2 -KH,co2 Pgas,co2 ) =dt Vgas Vgas
Tabla D-7.: Razones cin eticas para la transferencia de masa en la interfaz gaseosa. T,8 T,9 T,10 T,8 =(KL a)(Sliq,h2 -16KH,h2 Pgas,h2 ) T,9 =(KL a)(Sliq,ch4 -64KH,ch4 Pgas,ch4 ) T,10 =(KL a)(Sliq,co2 -KH,co2 Pgas,co2 )
Tabla D-8.: Ecuaciones algebraicas para la transferencia de masa en la interfaz l quido-gas. RT(OP ) Pgas,h2 Pgas,h2 = Sgas,h2 16 RT(OP ) Pgas,ch4 Pgas,ch4 = Sgas,ch4 64 Pgas,co2 Qgas Pgas,co2 =RT(OP ) Sgas,co2 RT(OP ) T,8 T,9 Qgas = Vliq ( + +T,10 ) Patm Pgas,H 2O 16 64
ANAEROBICA MODELO DIGESTION ADM-1 2 Saa fsi,xc Sfa Sva Sbu Spro Sac Sh2 Sch4 Sic Sin Si 3 4 5 6 7 8 9 10 11 12
Estados i
Procesos
Ssu
Desintegraci on
2 1 ffa,li (1-Ysu )fbu,su (1-Ysu )fpro,su (1-Ysu )fac,su (1-Ysu )fh2,su -1 -1 -1 -1 -1 -1 -1 (1-Yh2 ) (1-Yac ) (1-Ypro )0.57 (1-Ypro )0.43 (1-Yc4 )0.8 (1-Yc4 )0.2 -Yc4 Nbac Ci vi,10 -Ypro Nbac Ci vi,11 -Yac Nbac Ci vi,12 -Yac Nbac (1-Yc4 )0.54 (1-Yc4 )0.31 (1-Yc4 )0.15 -Yc4 Nbac (1-Yfa )0.7 (1-Yfa )0.3 -Yfa Nbac (1-Yaa )fva,aa (1-Yaa )fbu,a (1-Yaa )fpro,aa (1-Yaa )fac,aa (1-Yaa )fh2,aa Ci vi,6 -Yaa Nbac Ci vi,5 -Ysu Nbac
Hidr olisis de Carbohidratos
Hidr olisis de prote nas
Hidr olisis de l pidos
1-ffa,li
Consumo de azucares
-1
Consumo de amino acidos
Consumo de LCFA
Consumo de Valerato
Consumo de Butirato
10
Consumo de Propionato
11
Consumo de Acetato
12
Consumo de Hidr ogeno
13
Decaimiento de Xsu
14
Decaimiento de Xaa
15
Decaimiento de Xfa
16
Decaimiento de Xc4
17
Decaimiento de Xpro
18
Decaimiento de Xac
19
Decaimiento de Xh2
Butirato total
Propionato total
LCFA
Acetato total
Gas metano
Carbono Inorg anico [mol C/L]
Nitr ogeno Inorg anico [mol N/L]
Monosac aridos
Amino acidos
Valerato total
E. Anexo: Par ametros caso base ADM-1
Tabla E-1.: Matriz de Petersen: Coecientes bioqu micos i y tasas cin eticas i para los componentes solubles [14].
Inertes Solubles Gas hidr ogeno
Descripci on de los Estados [gDQO/L]
ANAEROBICA MODELO DIGESTION ADM-1
120
Estados i Procesos Desintegraci on -1 -1 -1 Ysu Yaa Yf a Y c4 Y c4 Ypro Yac Y h2 1 1 1 1 1 1 1 -1 -1 -1 -1 -1 -1 -1 Khyd,pr Xpr Khyd,ch Xch -1 fch,xc fpr,xc fli,xc fxi,xc Kdis Xc Xc Xch Xpr Xli Xsu Xaa Xf a X c4 Xpro Xac X h2 XI
13
14
15
16
17
18
19
20
21
22
23
24
Tasa Cin etica[gDQO/L]
2 Hidr olisis de prote nas Hidr olisis de l pidos Consumo de azucares
Hidr olisis de Carbohidratos
6 Consumo de LCFA Consumo de Valerato Consumo de Butirato
Consumo de amino acidos
10 Consumo de Acetato Consumo de Hidr ogeno Decaimiento de Xsu Decaimiento de Xaa Decaimiento de Xfa Decaimiento de Xc4 Decaimiento de Xpro Decaimiento de Xac Decaimiento de Xh2
Consumo de Propionato
11
12
13
Khyd,li Xli Ssu Km,su Xsu I1 Ks + Ssu Saa Km,aa Xaa I1 Ks + Saa Sf a Km,f a Xf a I2 Ks + S f a Sva 1 Km,c4 Xc4 I2 Ks + Sva 1 + Sbu /Sva Sbu 1 Km,c4 Xc4 I2 Ks + Sbu 1 + Sva /Sbu Spro Xpro I2 Km,pro Ks + Spro Sac Km,ac Xac I3 Ks + Sac S h2 Km,h2 Xh2 I1 Ks + Sh2 Kdec,Xsu Xsu Kdec,Xaa Xaa Kdec,Xf a Xf a Kdec,Xc4 Xc4 Kdec,Xpro Xpro Kdec,Xac Xac Kdec,Xh2 Xh2
14
15
16
17
18
19
l pidos
Consumidoras de propionato
Compuestos complejos
Carbohidratos
Prote nas
Consumidoras de amino acidos
Consumidoras de azucares
Consumidoras de LCFA
Consumidoras de Valerato y Butirato
Consumidoras de acetato
Descripci on de los Estados [gDQO/L]
Consumidoras de Hidr ogeno
Particulado Inerte
Factores de Inhibici on I1 =IpH IIN,lim I2 =IpH IIN,lim IH2 I3 =IpH IIN,lim INH3 ,Xac
Tabla E-2.: Matriz de Petersen: Coecientes bioqu micos i y tasas cin eticas i para los componentes insolubles en el modelo ADM-1. E Anexo: Par ametros caso base ADM-1
121
Tabla E-3.: Par ametros de la cin etica extracelular de desintegraci on e hidr olisis [d1 ].
Par ametro Kdis Khydch Khydli Khydpr Valor 0.5 10 10 10 Descripci on Tasa de desintegraci on de los compuestos complejos Tasa de hidr olisis de carbohidratos para la formaci on de monosac aridos Tasa de hidr olisis de l pidos para la formaci on de AGCL Tasa de hidr olisis de prote nas para la formaci on de amino acidos
Tabla E-4.: Par ametros de la cin etica biol ogica: Tasas m aximas de crecimiento bacteriano.
Par ametro Kmsu Kmaa Kmf a Kmc4 Kmpro Kmac Kmh2 Valor 30 50 6 20 13 8 35 Descripci on Biomasa consumidora de monosac aridos Biomasa consumidora de amino acidos Biomasa consumidora de acidos grasos Biomasa consumidora de valerato y butirato Biomasa consumidora de propionato Biomasa consumidora de acetato Biomasa consumidora de hidr ogeno
Tabla E-5.: Par ametros de la cin etica biol ogica: Tasas de producci on de la biomasa.
Par ametro Ysu Yaa Yf a Yc4 Ypro Yac Yh2 Valor 0.1 0.08 0.06 0.06 0.04 0.05 0.06 Descripci on Biomasa consumidora de monosac aridos Biomasa consumidora de amino acidos Biomasa consumidora de acidos grasos Biomasa consumidora de valerato y butirato Biomasa consumidora de propionato Biomasa consumidora de acetato Biomasa consumidora de hidr ogeno
Tabla E-6.: Par ametros de la cin etica biol ogica: Constantes de saturaci on media.
Par ametro Kssu Ksaa Ksf a Ksc4 Kspro Ksac Ksh2 Valor 0.5 0.3 0.4 0.2 0.1 0.15 0.000007 Descripci on Biomasa consumidora de monosac aridos Biomasa consumidora de amino acidos Biomasa consumidora de acidos grasos Biomasa consumidora de valerato y butirato Biomasa consumidora de propionato Biomasa consumidora de acetato Biomasa consumidora de hidr ogeno
122
E Anexo: Par ametros caso base ADM-1
Tabla E-7.: Par ametros de la cin etica biol ogica: Tasas de decaimiento de la biomasa microbiana.
Par ametro Kdecxsu Kdecxaa Kdecxf a Kdecxc4 Kdecxpro Kdecxac Kdecxh2 Valor 0.02 0.02 0.02 0.02 0.02 0.02 0.02 Descripci on Biomasa consumidora de monosac aridos Biomasa consumidora de amino acidos Biomasa consumidora de acidos grasos Biomasa consumidora de valerato y butirato Biomasa consumidora de propionato Biomasa consumidora de acetato Biomasa consumidora de hidr ogeno
Tabla E-8.: Par ametros de la cin etica biol ogica: Par ametros para la inhibici on debida al pH.
Par ametro pHul pHll pHul ac pHll ac pHul h2 pHll h2 Valor 5.5 4 7 6 6 5 Descripci on Limite superior de inhibici on por pH en los procesos de consumo de monosac aridos, amino acidos, acidos grasos, valerato y butirato. Limite inferior de inhibici on por pH en los procesos de consumo de monosac aridos, amino acidos, acidos grasos, valerato y butirato. Limite superior de inhibici on por pH en el proceso de consumo de acetato. Limite inferior de inhibici on por pH en el proceso de consumo de acetato. Limite superior de inhibici on por pH en el proceso de consumo de hidr ogeno. Limite inferior de inhibici on por pH en el proceso de consumo de hidr ogeno.
Tabla E-9.: Par ametros de la cin etica biol ogica: Par ametros para la inhibici on debida a la concentraci on del H2 .
Par ametro Ki h2 fa Ki h2 c4 Ki h2 pro Valor 0.000005 0.00001 0.0000035 Descripci on Bacterias que consumen acidos grasos Bacterias que consumen valerato y butirato Bacterias que consumen propionato
Tabla E-10.: Par ametros de la cin etica biol ogica: Par ametros para la inhibici on debida a la concentraci on del NH3 .
Par ametro Ki nh3 Valor 0.000005 Descripci on Bacterias que consumen acetato
123
Tabla E-11.: Par ametros de la cin etica biol ogica: Par ametros para la inhibici on debida a la concentraci on de Nitr ogeno.
Par ametro Ki in3 Valor 0.0001 Descripci on Inhibici on debida a la limitaci on del Nitr ogeno en los nutrientes
Tabla E-12.: Contenido de carbono en las variables de estado del modelo.

Par ametro Csu Cbu Cpro Cac Cbiom Caa Cva Cch4 Valor 0.03125 0.0250 0.03125 0.03125 0.0313 0.03 0.024 0.015625 Descripci on Contenido de Contenido de Contenido de Contenido de Contenido de Contenido de Contenido de Contenido de carbono carbono carbono carbono carbono carbono carbono carbono en en en en en en en en los monosac aridos el butirato el propionato el acetato los 7 grupos de bacterias los amino acidos el valerato el metano
Tabla E-13.: Constantes de Henry para la transferencia de masa en la interfaz l quido-gas.

Par ametro KH ch4 KH co2 KH h2 Valor 0.001161 0.002713 0.0007384 Descripci on Constante de Henry para el gas metano Constante de Henry para el gas carb onico Constante de Henry para el gas hidr ogeno
Tabla E-14.: Coecientes de equilibrio acido-base.

Par ametro Ka ac Ka bu Ka co2 Kw Ka nh4 Ka pro Ka va Valor 1.7378e5 1.4454e5 4.9375e7 2.0803e14 1.1110e9 1.3182e5 1.5848e5 Descripci on Equilibrio HAc/AcEquilibrio HBu/BuEquilibrio CO2 /HCO 3 Equilibrio H2 O/[OH ][H+ ] Equilibrio NH+ 4 /NH3 Equilibrio HPro/Pro Equilibrio HVa/Va
124
E Anexo: Par ametros caso base ADM-1
Tabla E-15.: Par ametros estequiom etricos de las reacciones bioqu micas.
Par ametro fsi xc fxi xc fch xc fpr xc i xc a li fh2 su fbu su fpro su fh2 aa fva aa fbu aa fpro aa fac aa fac su Valor 0.1 0.2 0.2 0.2 0.3 0.95 0.19055 0.1328 0.269 0.06 0.23 0.26 0.05 0.4 0.40755 Descripci on Tasa de formaci on de los solubles inertes obtenidos del consumo de los compuestos complejos (Xc) Tasa de formaci on de los particulados inertes obtenidos del consumo de los compuestos complejos (Xc) Tasa de formaci on de los carbohidratos obtenidos del consumo de los compuestos complejos (Xc) Tasa de formaci on de los prote nas obtenidos del consumo de los compuestos complejos (Xc) Tasa de formaci on de los l pidos obtenidos del consumo de los compuestos complejos (Xc) Tasa de formaci on de los acidos grasos obtenidos del consumo de los l pidos Tasa de formaci on de hidr ogeno obtenidos del consumo de los monosac aridos Tasa de formaci on del butirato obtenidos del consumo de los monosac aridos Tasa de formaci on del propionato obtenidos del consumo de los monosac aridos Tasa de formaci on de hidr ogeno obtenidos del consumo de los amino acidos Tasa de formaci on de valerato obtenidos del consumo de los amino acidos Tasa de formaci on de butirato obtenidos del consumo de los amino acidos Tasa de formaci on de propionato obtenidos del consumo de los amino acidos Tasa de formaci on de acetato obtenidos del consumo de los amino acidos Tasa de formaci on de acetato obtenidos del consumo de los monosac aridos
Tabla E-16.: Contenido de nitr ogeno en las variables de estado del modelo.
Par ametro Nbac Naa Valor 0.00625 0.007 Descripci on Contenido de carbono en los 7 grupos de bacterias Contenido de carbono en los amino acidos
125
Tabla E-17.: Concentraciones iniciales de las especies para la operaci on en r egimen estacionario, caso base modelo (ADM-1).
Par ametro Ssu i Saa i Sf a i Sva i Sbu i Spro i Sac i Sh2 i Sch4 i SIC i SIN i SI i Xc i Xch i Xpr i Xli i Xsu i Xaa i Xf a i Xc4 i Xpro i Xac i Xh2 i XI i SCAT i SAN i Valor 0.01261 0.00563 0.1098 0.01172 0.01485 0.01796 0.04923 2.55e7 0.01261 0.07506 0.03244 1.6707 1.040 0.0103 0.0103 0.01551 0.3541 0.3568 0.3906 0.1448 0.0609 0.4703 0.2237 19.9414 0.04 0.00357 Descripci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on Concentraci on de de de de de de de de de de de de de de de de de de de de de de de de de de monosac aridos amino acidos acidos grasos valerato butirato propionato acetato gas hidr ogeno gas metano carbono inorg anico nitr ogeno inorg anico solubles inertes compuestos complejos carbohidratos prote nas l pidos bacterias consumidoras de monosac aridos bacterias consumidoras de amino acidos bacterias consumidoras de acidos grasos bacterias consumidoras de valerato y butirato bacterias consumidoras de propionato bacterias consumidoras de acetato bacterias consumidoras de hidr ogeno particulado inerte cationes de bases fuertes[equivalentes] aniones de acidos fuertes[equivalentes]
Tabla E-18.: Par ametros de operaci on del reactor anaer obico ADM-1 utilizado en el caso base del modelo.
Par ametro Q V liq V gas Kla T Pgas Pgas h2o Valor 170 3400 300 200 303 1.013 0.0557 Descripci on Caudal[m3 /d a] Volumen de la fase l quida en el reactor [m3 ] Volumen de la fase gaseosa en el reactor [m3 ] Coeciente de transferencia de masa multiplicado por el area especica de transferencia [d1 ] Temperatura promedio en el reactor [K] Presi on total en el volumen ocupado por los gases [bar] Presi on del vapor de agua [bar]
Bibliograf a
[1] Ahring, B. ; Angelidaki, I. ; Johansen, K.: Anaerobic treatment of manure together with industrial waste. En: Water Science and Technology 25 (1992), p. 311318 n, G.: Low temperature anaerobic digestion of mixtures of llama, Alvarez, R. ; Lide cow and sheep manure for improved methane production. En: Biomass and Bioenergy 33 (2009), p. 527533 Andrews, J. ; Graef, S.: Dynamic modeling and simulation of the anaerobic digestion process. En: Anaerobic Biological Treatment Processes. Advances in Chemistry Series-Amer. Chemical Soc., Washington 105 (1971), p. 126162 Angelidaki, I. ; Ahring, B.: Eect of the clay mineral bentonite on ammonia inhibition of anaerobic thermophilic reactors degrading animal waste. En: Biotech. Bioeng. 3 (1993), p. 409414 Angelidaki, I. ; Ellegaard, L. ; Ahring, B.: A mathematical model for dynamic simulation of anaerobic digestion of complex substrates: Focusing on ammonia inhibition. En: Biodegradation 42 (1993), p. 159166 APHA: Standard Methods for the Examination of Water and Wastewater. 21th ed. Washington, DC. : American Public Health Association, 2005 ASTM: Standard Test Method for Determining the Anaerobic Biodegradation Potential of Organic Chemicals. En: American Society for Testing and Materials, West Conshohocken, PA. (1992) Atkinson, B.: Reactores Bioquimicos. Madrid : Reverte, 2008 Bailey, J.: Mathematical modeling and analysis in biochemical engineering: Past accomplishments and future opportunities. En: Biotechnology Program 14 (1998), p. 820
[2]
[3]
[4]
[5]
[6]
[7]
[8] [9]
[10] Baily, J. ; Ollis, D.: Biochemical engineering fundamentals, 2nd ed. New York : McGraw-Hill, 1986 [11] Balch, W. ; Fox, G. ; Magrum, L. ; Woese, C. ; Wolfe, R.: Methanogens: Reevaluation of a unique biological group. En: Microbiol. Rev. 43 (1979), p. 260269 [12] Banks, C. ; Stentiford, E.: Biodegradable municipal solid waste: biotreatment options. En: Water resource magnagement. 160 (2007), p. 1118
Bibliograf a
127
[13] Batstone, D.: High rate anaerobic treatment of complex wastewater, Universidad de Queensland, Brisbane, Tesis de Doctorado, 2000 [14] Batstone, D. ; Keller, J. ; Angelidaki, I. ; Kalyuzhnyi, S. ; Pavlostathis, S. ; Rozzi, A. ; Sanders, W. ; Siegrist, H. ; Vavilin, V.: Anaerobic Digestion Model1,(ADM1). En: IWA Task Group for Mathematical Modelling of Anaerobic Digestion Processes. IWA Publishing London. (2002) [15] Battersby, N. ; Wilson, V.: Evaluation of a serum bottle technique for assessing the anaerobic biodegradability of organic chemicals under methanogenic conditions. En: Chemosphere 17 (1988), p. 24412460 [16] Baumann, U. ; Muller, M.: Determination of anaerobic biodegradability with a simple continuous xed-bed reactor. En: Wat. Res. 31 (1997), p. 15131517 [17] Blumensaat, F. ; Keller, J.: Modelling of two-stage anaerobic digestion using the IWA Anaerobic Digestion Model No. 1 (ADM1). En: Water Research 39 (2005), p. 171183 [18] Bohn, I. ; Bjornsson, L. ; Mattiasson, B.: The energy balance in farm scale anaerobic digestion of crop residues at 11-37o C. En: Process Biochemistry 42 (2007), p. 5764 [19] Bolivar, Efr en: Determinaci on de la biodegradabilidad anaerobica del residuo l quido de una industria de sabaj on, Universidad Nacional de Colombia, Tesis de Grado, 2004 [20] Boubaker, F. ; Cheikh, B.: Modelling of the mesophilic anaerobic co-digestion of olive mill wastewater with olive mill solid waste using anaerobic digestion model No. 1 (ADM1). En: Bioresource Technology 99 (2008), p. 65656577 [21] Brinkman, J.: Anaerobic digestion of mixed waste slurries from kitchens, slaughterhouses and meat processing industries. En: Proceedings of the II International symposium on Anaerobic Digestion of Soild Waste., Barcelona (1999), p. 190191 [22] Callaghan, F. ; Wase, J. ; Thayanithy, K. ; Forster, C.: Co-digestion of waste organic solids: batch studies. En: Bioresource Technology 67 (1999), p. 117122 [23] Castells, X.: Tratamiento y valorizacion energetica de residuos. Madrid : Ediciones Diaz de Santos, 2005 [24] Chanakya, H. ; Srikumar, K. ; Anand, V. ; Modak, J. ; K., Jagadish: Fermentation properties of agro-residues, leaf biomass and urban market garbage in a solid phase biogas fermenter. En: Biomass and Bioenergy 16 (1999), p. 417429
128
Bibliograf a
[25] Chen, Y. ; Cheng, J. ; Creamer, K.: Inhibition of anaerobic digestion process: A review. En: Bioresource Technology. 99 (2008), p. 40444064 [26] Chen, Y. ; Cheng, J. ; Creamer, K.: Inhibition of anaerobic digestion process: A review. En: Bioresource Technology 99 (2008), p. 40444064 [27] Collazos, C.: Biodegradabilidad anaerobia de euentes cerveceros, Universidad Nacional de Colombia, Tesis de Grado (Magister en Ingenier a Ambiental), 2000 [28] Costello, D. ; Greenfield, P. ; Lee, P.: Dynamic modelling of a single-stage high-rate anaerobic reactor - I. Model derivation. En: Wat. Res. 25 (1991), p. 847858 [29] Cuellar, S. ; Cuesta, M. ; Medina, J. ; Rodriguez, M. ; Ruiz, M. ; Villareal, A.: Perl de producto - Naranja. En: Ministerio de Agricultura y Desarrollo Rural Corporacion Colombia Internacional 9 (2000), p. 112 [30] Cun-fang, L. ; Xing-zhong, Y. ; Guang-ming, Z. ; Wen-wei, J.: Prediction of methane yield at optimum pH for anaerobic digestion of organic fraction of municipal solid waste. En: Bioresource Technology 99 (2008), p. 882888 [31] Deloitte, Asesores y Consultores L. Informe ejecutivo de gesti on Codensa S.A. E.S.P. Junio 2009 [32] Diaz Baez, M. ; Espitia, S. ; Molina, F.: Digesti on anaerobia - Una aproximaci on a la tecnolog a. Bogot a : Universidad Nacional de Colombia, Instituto de biotecnolog a. Colciencias. Unibiblos, 2002 [33] Eastman, J. ; Ferguson, J.: Solubilization of particulate organic carbon during the acid phase anaerobic digestion. En: J Wat Pollut Control Fed 53 (1981), p. 352366 [34] El-Mashad, H. ; Zeeman, G. ; Wilko, K. ; Bot, G. ; Lettinga, G.: Eect of temperature and temperature uctuation on thermophilic anaerobic digestion of cattle manure. En: Bioresource Technology 95 (2004), p. 191201 [35] El-Mashad, H. ; Zhang, R.: Biogas production from co-digestion of dairy manure and food waste. En: Bioresource Technology 101 (2010), p. 40214028 [36] Eldem, N. ; Ozturk, I. ; Soyer, E. ; Calli, B. ; Akgiray, O.: Ammonia and pH inhibition in anaerobic treatment of wastewaters, Part I: Experimental. En: Journal of Environment science and health 39 (2004), p. 24052420 [37] Fan, L. ; Chen, L. ; Mehta, C. ; Chen, Y.: Energy and available energy contents of cattle manure and digester sludge. En: Agricultural Wastes 13 (1985), p. 239249 [38] Fannin, K.: Start-up, operation, stability, and control in Anaerobic digestion of biomass. En: Elsevier applied science 29 (1987), p. 101106
Bibliograf a
129
[39] Feng, Y. ; Behrendt, J. ; Wendland, C. ; Otterpohl, R.: Parameters analysis and discussion of the anaerobic digestion model no. 1 (ADM-1). En: Water Sci. Technol 54 (2006), p. 139147 [40] Fernandez, J. ; Omil, F. ; Mendez, R. ; Lema, J.: Anaerobic treatment of breboard manufacturing wastewaters in a pilot scale hybrid USBF reactor. En: Wat. Res. 35 (2001), p. 41504158 [41] Field, J.: Par ametros operativos del manto anaerobios de lodos de ujo ascendente Arranque y operaci on de sistemas de ujo ascendente con mantos de lodo UASB. En: Univalle-Universidad Agr cola de Wageningen. Cali, Colombia. 3 (1987), p. 409414 [42] Field, J. ; Sierra, R.: Biodegradability and toxicity lecture series. En: International course on Anaerobic Wastewater treatment. Agricultural University Wageningen. The Netherlands. 7 (1989) iga, Juan: Evaluacion de la actividad metanogenica de lodos, [43] Forero, Wiliam ; Zun Universidad Nacional de Colombia, Tesis de Grado, 2003 [44] Gal , A. ; Benabdallah, T. ; Astals, S. ; Mata-Alvarez, J.: Modied version of ADM-1 model for agro-waste application. En: Bioresource Technology 100 (2009), p. 27832790 [45] GATE, German Agency for Technical C.: Biogas Digesters. Eschborn : GTZ, 1990 [46] Gavala, H. ; Skiades, I. ; Bozinis, N. ; Lyberatos, G.: Anaerobic co-digestion of agricultural industries wastewaters. En: Water Science and Technology 34 (1996), p. 6775 [47] Georgacakis, D. ; Sievers, D. ; Ianotti, E.: Buer stability in manure digesters. En: Agric. Wastes 4 (1982), p. 427441 [48] Gottschalk, G.: Bacterial Metabolism, second ed. New York : Springer-Verlag, 1986 evara, M.: Proyecto para la producci [49] Gu on y comercializaci on de naranja. En: Centro de Comercio Internacional, UNCTAD/GATT., Ginebra (1994), Nr. 3, p. 112 [50] Gujer, W. ; Zehnder, A.: Conversion processes in anaerobic digestion. En: Wat Sci Technol 15 (1983), p. 127167 [51] Gunaseelan, N.: Anaerobic digestion of biomass for methane production: a review. En: Biomass Bioenergy 13 (1997), p. 83114 [52] Gunaseelan, N.: Biochemical methane potential of fruits and vegetable solid waste feedstocks. En: Biomass Bioenergy 26 (2004), p. 389399
130
Bibliograf a
[53] Gunaseelan, N.: Regression models of ultimate methane yields of fruits and vegetable solid wastes, sorghum and napiergrass on chemical composition. En: Bioresource Technology 98 (2007), p. 12701277 [54] Hamzawi, N. ; Kennedy, K. ; Mclean, D.: Review of applications of high-solids anaerobic digestion to solid waste types. En: Journal of Solid Waste Techno. and Manag. 26 (1999), p. 119132 [55] Hansen, T. ; Schmidt, J. ; Angelidaki, I. ; Jensen, J. ; Mosbaek, H.: Method for determination of methane potentials of solid organic waste. En: Waste Manage 24 (2004), p. 393400 [56] Hashimoto, A. ; Chen, Y.: Methane from cattle wastes: eects of temperature, hidraulic retention time and inuent substrate concentration on kinetic parameter (K). En: Biotech Bioeng 24 (1982), p. 20392052 [57] Hill, D.: A comprehensive dynamic model for animal waste methanogenesis. En: Trans. ASAE 25 (1982), p. 13741380 [58] Hongwei, Y. ; Zhanpeng, J. ; Shaoqi, S.: Dehydrogenase activity test for assessing anaerobic biodegradability of organic compounds. En: Ecotoxicol. Environ. Saf. 53 (2002), p. 416421 [59] Hongwei, Y. ; Zhanpeng, J. ; Shaoqi, S.: Biodegradability of nitrogenous compounds under anaerobic conditions and its estimation. En: Ecotoxicology and Environmental Safety 63 (2006), p. 299305 [60] ISO: Water quality-evaluation of the ultimate anaerobic biodegradability of organic compounds in digested sludge. Method by measurement of the biogas production. En: International Organization for Standardization, Seattle. (1995) [61] Jae, K. ; Soon, P.: Biochemical methane potential and solid state anaerobic digestion of korean food wastes. En: Bioresource Technology 52 (1995), p. 245253 [62] Jain, S. ; Lala, A. ; Bhatia, S. ; Kudchadker, A.: Modelling of hydrolysis controlled anaerobic digestion. En: J. Chem. Technol. Biotechnol. 53 (1992), p. 337344 [63] Jianguo, J. ; Xuejuan, D. ; Siio, N. ; Chang, Z.: Comparison of atmospheric pressure eects on the anaerobic digestion of municipal solid waste. En: Bioresource Technology 101 (2010), p. 63616367 [64] Jih-Gaw, L. ; Ying-Shih, M. ; Allen, C. ; Cheng-Lung, H.: BMP test on chemically pretreated sludge. En: Bioresource Technology 68 (1999), p. 187192
Bibliograf a
131
[65] Kameya, T. ; Murayama, T. ; Kitano, M. ; Urano, K.: Testing and classication methods for the biodegradabilities of organic compounds under anaerobic conditions. En: The Science of the Total Environment 170 (1995), p. 3141 [66] Kaparaju, P. ; Rintala, J.: Thermophilic anaerobic digestion of industrial orange waste. En: Environmental Technology 27 (2006), p. 623633 [67] Karim, K. ; Hoffmann, R. ; Klasson, T. ; Al-Dahhan, M.: Anaerobic digestion of animal waste: Eect of mode of mixing. En: Water Research 39 (2005), p. 35973606 [68] Kayhanian, M. ; Rich, D.: Pilot-scale high solids thermophilic anaerobic digestion of municipal solid waste with an emphasis on nutrient requirements. En: Biomass and Bioenergy 26 (1995), p. 433444 [69] Khanal, S.: Anaerobic biotechnology for bioenergy production - principles and applications. Iowa, USA : Blackwell Publishing, 2008 [70] Klesov, A. ; Rabinovich, M.: Enzymatic hydrolysis of cellulose. En: Itogi nauki i technik 12 (1978), p. 499 [71] Klinger, I. ; Alexander, R. ; Marchaim, U: A suggested solution for slaughterhouse wastes: Uses of the residual materials after anaerobic digestion. En: Bioresource Technology 37 (1991), p. 127134 [72] Koster, I.: Liquefaction and acidogenesis of tomatoes in an anaerobic two-phase solid-waste treatment system. En: Agric Wastes 11 (1984), p. 241252 [73] Kugelman, I. ; Chin, K.: Toxicity, sinergism and antagonism in anaerobic waste treatment processes. En: Anaerobic treatment Biological processes, ACS-Washington D.C. (1971), p. 100105 [74] Kuo, W. ; Cheng, K.: Use of respirometer in evaluation of process and toxicity of thermophilic anaerobic digestion for treating kitchen waste. En: Bioresource Technology 98 (2007), p. 18051811 [75] Lapuerta, M. ; Hernandez, J. ; Pazo, A. ; Lopez, J.: Gasication and cogasication of biomass wastes: eect of the biomass origin and the gasier operating conditions. En: Fuel Process. Technol. 89 (2008), p. 828837 [76] Lay, J. ; Lee, J. ; Noike, T.: Feasibility of biological hydrogen production from organic fraction of municipal solid waste. En: Water Res. 33 (1999), p. 25792586 [77] Lehtomaki, A. ; Vavilin, V. ; Rintala, J.: Kinetic analysis of methane production from energy crops. En: Proceedings of the Fourth International Symposium on Anaerobic Digestion of Solid Waste Copenhagen, Denmark 64 (2005), p. 6772
132
Bibliograf a
[78] Lesteur, M. ; Bellon-Maurel, V. ; Gonzalez, C. ; Latrille, E.: Alternative methods for determining anaerobic biodegradability: A review. En: Process Biochemistry 45 (2010), p. 431440 [79] Lettinga, G. ; McMahon, K.D. ; Mackie, R.I. ; Raskin, L.: Anaerobic digestion and wastewater treatment systems. En: Rev. Antonie Van Leeuwenhoek 67 (1995), p. 115123 [80] Lin, A.: Eect of heavy metals on volatile fatty acid degradation in anaerobic digestion. En: Water Research 26 (1992), p. 177183 [81] Lubken, M. ; Wichern, M. ; Schlattmann, M. ; Gronauer, A. ; Horn, H.: Modelling the energy balance of an anaerobic digester fed with cattle manure and renewable energy crops. En: Water Research 41 (2007), p. 40854096 [82] Madigan, M. ; Martinko, J. ; Parker, J.: Biologia de los microorganismos 10a edicion. Illinois : Pearson - Prentice hall, 1999 [83] Mamma, D. ; Kourtoglou, E. ; Christakopoulos, P.: Fungal multienzyme production on industrial by-products of the citrus-processing industry. En: Bioresour. Technol. 99 (2008), p. 23732383 [84] Marchaim, U.: Biogas processes for sustainable development. En: Food and Agriculture Organization of the United Nations (FAO), FAO Agricultural Service Bulletin 95 (1992) s, P. ; Pavan, P.: Performance of digesters [85] Mata-Alvarez, J. ; Cecchi, F. ; Llabre treating the organic fraction of municipal solid waste dierently sorted. En: Biol Wastes 33 (1990), p. 181199 [86] McCarty, P.: Anaerobic waste treatment fundamentals. En: Public Works 95 (1964), p. 9194 [87] McCarty, P. ; McKinney, R.: Salt toxicity in anaerobic digestion. En: Journal of w ater pollution control federation. 33 (2004), p. 399415 [88] McCarty, P. ; Mosey, F.: Modelling of anaerobic digestion process (A discussion of concepts). En: Water Research 24 (1991), p. 1733 [89] Metcalf, L. ; Eddy, H.: Wastewater Enginnering. 4th. New York : Mc Graw-Hill, Inc., 2003 [90] Mladenovska, Z. ; Dabrowski, S. ; Ahring, B.: Anaerobic digestion of manure and mixture of manure with lipids: biogas reactor performance and microbial community analysis. En: Water Science and Technology 48 (2003), p. 271279
Bibliograf a
133
[91] Moller, H. ; Sommer, S. ; Ahring, B.: Methane productivity of manure, strawand solid fractions of manure. En: Biomass and Bioenergy 26 (2004), p. 485495 [92] Molnar, L. ; Bartha, I.: Factors inuencing solid-state anaerobic digestion. En: Biol. wastes 28 (1989), p. 1524 [93] Monod, J.: La technique de culture continue: theorie et applications. En: Ann. Inst. Pasteur, Paris 79 (1950), p. 390410 [94] Montgomery, H. ; Dymock, J. ; Thom, N.: The rapid colorimetric determination of organic acids and their salts in sewage-sludge liquor. En: The Analyst. 87 (1962), p. 949952 [95] Montusiewicz, A. ; Lebiocka, M. ; Pawlowska, M.: Characterization of the biomethanization process in selected waste mixtures. En: Archives of Environmental Protection 34 (2008), p. 4961 [96] Motteta, A. ; Franc ois, E. ; Latrille, E. ; Steyer, J. ; Deleris, S. ; Vedrenne, F. ; Carrere, H.: Estimating anaerobic biodegradability indicators for waste activated sludge. En: Chemical Engineering Journal 160 (2010), p. 488496 [97] Myung-Yeol, L. ; Chang-Won, S. ; Yong-Tae, A. ; Hang-Sik, S.: Variation of ADM-1 by using temperature-phased anaerobic digestion (TPAD) operation. En: Bioresource Technology 100 (2009), p. 28162822 [98] Neves, L. ; Oliveira, R. ; Alves, M.: Inuence of inoculum activity on the biomethanization of a kitchen waste under dierent waste/inoculum ratios. En: Process Biochemistry 39 (2004), p. 20192024 [99] Nichols, C.: Overview of anaerobic digestion technologies in Europe. En: Biocycle. 45 (2004), p. 4763 [100] Nordberg, A. ; Edstrom, M.: Codigestion of energy crops and the source sorted organic fraction of municipal solid waste. En: Water science and technology. 52 (2005), p. 217222 [101] Owen, W. ; Stuckey, D. ; Healy, J. ; Young, L. ; McCarty, P.: Bioassay for monitoring biochemical methane potential and anaerobic toxity. En: Wat.Res. 13 (1978), p. 485492 [102] Owens, J. ; Chynoweth, D.: Biochemical methane potential of MSW components. En: Water Sci. Technol. 27 (1993), p. 114 [103] Palmowski, L. ; Muller, J.: Inuence of the size reduction of organic waste on their anaerobic digestion. En: Water Sci. Technol. 41 (2000), p. 155162
134
Bibliograf a
[104] Palmowski, L. ; Muller, J.: Anaerobic degradation of organic materials - signicance of the substrate surface area. En: Water Sci. Technol. 47 (2003), p. 231238 [105] Parker, W.: Application of the ADM-1 model to advanced anaerobic digestion. En: Bioresour. Technol. 96 (2005), p. 18321842 mez, E.: Kinetics of anaerobic treatment: a critical [106] Pavlostathis, S. ; Giraldo-Go review. En: Critical reviews in environmental Control 21 (1991), p. 411490 [107] Petersen, E.: Chemical Reaction Analysis. Englewood Clis, NJ, USA : PrenticeHall, 1965 [108] Prescott, L. ; Harley, J. ; Klein, D.: Microbiology 5th edition. New York : McGrawHill, 2002 al, A. ; Franco, A. ; Carrasco, F. ; Rodriguez, J. ; Roca, E. ; Lema, [109] Pun J.: Advanced Monitoring and control of anaerobic wastewater treatment plants. En: System analysis and computing in water quality management 13 (2000), p. 109124 [110] Rajakumar, R. ; Meenambal, T. ; Rajesh Banu, J. ; Yeom, T.: Treatment of poultry slaughterhouse wastewater in upow anaerobic lter under low upow velocity. En: International Journal of Environment Science and Technology 8 (2010), p. 149158 [111] Ramsay, I.: Modelling and control of high-rate anaerobic wastewater anaerobic treatment systems. Brisbane, University of Queensland, Tesis de Doctorado, 1997 [112] Raposo, F. ; Banks, C. ; Siegert, I. ; Heaven, S. ; Borja, R.: Inuence of inoculum to substrate ratio on the biochemical methane potential of maize in batch tests. En: Process Biochemistry 41 (2006), p. 14441450 [113] Rodriguez, J. ; Ruiz, G. ; Roca, E. ; Lema, J.: Modication of the IWA-ADM1 for the application to anaerobic treatment of ethanolic wastewater from wine factories [Telemac contribution-4]. En: 10th IWA World Congress Anaerobic Digestion 28th August - 2nd September 2004, Montreal (Canada) rez, M. ; Sales, D.: Digesti [114] Romero, L. ; de la Rubia, M. ; Pe on anaer obica termof lica de lodos de EDAR. En: Revista T ecnica de Medio Ambiente-RETEMA 15 (1991), p. 9097 [115] Romli, M. ; Keller, J. ; Lee, P. ; Greenfield, P.: Model prediction and verication of a two-stage high-rate anaerobic wastewater treatment system subject to shock loads. En: Process Saf. Environ. Protect. 73 (1995), p. 151154
Bibliograf a
135
[116] Rosen, C. ; Vrecko, D. ; Gernaey, K. ; Pons, M. ; Jeppsson, U.: Implementing ADM-1 for plant-wide benchmark simulations in Matlab/Simulink. En: Wat. Sci. Tech. 54 (2006), p. 1120 [117] Sanchez, J. ; Valle, L. ; Rodriguez, F. ; Moriago, M. ; Borrego, J.: Inhibition of methanogenesis by several heavy metals using pure cultures. En: Letters in applied microbiology 23 (1996), p. 439 444 [118] Schober, G. ; Schafer, J. ; Schmid-Staiger, U. ; Trosch, W.: One and two stage digestion of solid organic waste. En: Water Research 33 (1999), p. 854 860 enz, C. ; Este vez, M. ; Sanhueza, S.: Utilizaci [119] Sa on de residuos de la industria de jugos de naranja como fuente de bra diet etica en la elaboraci on de alimentos. En: Organo Ocial de la Sociedad Latinoamericana de Nutrici on, Depto. de Agroindustria y Enolog a. Facultad de Ciencias Agron omicas. Universidad de Chile. Santiago, Chile 57 (2007), p. 186191 [120] Shelton, D. ; Tiedje, J.: General method for determining anaerobia biodegradation potential. En: Applied and enviromental microbiology 47 (1984), p. 850857 [121] Siegrist, H.: Mathematical model for meso- and thermophilic anaerobic sewage sludge digestion. En: Environmental Science and Technology. 36 (2002), p. 11131123 [122] Siles, J. ; Li, Q. ; Thompson, I.: Biorenery of waste orange peel. En: Crit. Rev. Biotechnol. 30 (2010), p. 6369 [123] Sleat, R. ; Mah, R.: Hydrolytic bacteria in Anaerobic digestion of biomass. En: Elsevier applied science 12 (1987), p. 226237 [124] Speece, R.: Anaerobic Biotechnology for Industrial Wastewaters. En: Nashville, TN: Archae Press 3 (1996) [125] Spreen, T.: The citrus industries of the United States and Mexico after Nafta. En: Revista Chapingo, Serie Horticultura 6 (2000), p. 145152 [126] Stams, A.: Metabolic interactions between anaerobic bacteria in methanogenic environments. En: Antonie van Leeuwenhoek. 66 (1994), p. 271294 [127] Stewart, J. ; Bhattacharya, S. ; Madura, R. ; Mason, S. ; Schonberg, J.: Anaerobic treatability of selected organic toxicants in petrochemical wastes. En: Water Research 29 (1995), p. 27302738 [128] Stronach, S. ; Rudd, T. ; Lester, J.: Anerobic digesti on processes in wastewater treatment. En: Berlin: Springer. (1986)
136
Bibliograf a
[129] Stroot, P. ; McMahon, K. ; Mackie, R. ; Raskin, L.: Anaerobic codigestion of municipal solid waste and biosolids under various mixing condition: I. Digester performance. En: Water Research 24 (2001), p. 18041816 [130] Stryer, L.: Biochemistry. New York : Freeman and Company, 1988 [131] Thamsiriroj, T. ; Murphy, J.: Modelling mono-digestion of grass silage in a 2-stage CSTR anaerobic digester using ADM1. En: Bioresource Technology 102 (2011), p. 948959 a, M.: Anaerobic biodegradability and [132] Vallecillo, A. ; Garcia-Encina, P. ; Pen toxicity of chlorophenols. En: Water Sci. Technol. 40 (1999), p. 161168 [133] Van Lier, J. ; Martin, J. ; Lettinga, G.: Eect of temperature on the anaerobic thermophilic conversion of volatile fatty acids by dispersed and granular sludge. En: Water Res. 30 (1996), p. 199207 [134] Vandevivere, P. ; De Baere, L. ; Verstraete, W.: Biometanization of the organic fraction of municipal solid wastes. En: IWA publishing company. 10 (2003) [135] Vavilin, V. ; Rytov, S. ; Lokshina, L. ; Himmel, M.: A description of hydrolysis kinetics in anaerobic degradation of particulate organic matter. En: Bioresource Technol. 56 (1996), p. 229237 [136] Vavilin, V. ; Vasiliev, V. ; Ponomarev, A.: Simulation model methane as a tool for eective Biogas production during anaerobic conversion of complex organic matter. En: Bioresource Technology 52 (1994), p. 171177 [137] Veeken, B. ; Hamelers, B. ; Johansen, K.: Eect of temperature on hydrolysis rate of selected biowise components. En: Bioresour Technol 69 (1999), p. 249254 [138] Vidal, G. ; Carvalho, A. ; Mendez, R. ; Lema, J.: Inuence of the content in fats and proteins on the anaerobic biodegradability of dairy wastewaters. En: Bioresource Technology 74 (2000), p. 231239 [139] Viswanath, P. ; Devi, S. ; Krishna, N.: Anaerobic digestion of fruit and vegetable processing wastes for biogas production. En: Bioresource Technology 40 (1992), p. 4348 [140] Vogt, G. ; Liu, H. ; Kennedy, K. ; Holbein, B.: Super blue box recycling (SUBBOR) enhanced two-stage. Anaerobic digestion process for recycling municipal solid waste: laboratory pilot studies. En: Bioresource Technology 85 (2002), p. 291299
Bibliograf a
137
[141] Von Munch, E. ; Keller, J. ; Lant, P. ; Newell, R.: Mathematical modelling of prefermenters-I; model development and verication. En: Water Resource 33 (1994), p. 27572768 [142] Ward, A. ; Hobbs, P. ; Holliman, P. ; Jones, D.: Optimisation of the anaerobic digestion of agricultural resources. En: Bioresource Technology 99 (2008), p. 79287940 [143] Watanabe, T. ; Kimura, K. ; Asakawa, S.: Distinct members of a stable methanogenic archaeal community transcribe mcrA genes under ooded and drained conditions in Japanese paddy eld soil. En: Soil Biology and Biochemistry 41 (2009), p. 276285 [144] Wendland, C.: Anaerobic Digestion of Blackwater and Kitchen Refuse, Technischen Universit at Hamburg-Harburg, Tesis de Doctorado, 2008 [145] Wilkins, M. ; Widmer, W. ; Grohmann, K.: Simultaneous saccharication and fermentation of citrus peel waste by Saccharomyces cerevisiae to produce ethanol. En: Process Biochem. 42 (2007), p. 16141619 [146] Wittmann, C. ; Zeng, A. ; Deckwer, W.: Growth inhibition by ammonia and use of pH controlled feeding strategy for the eective cultivation of microbacterium chlorophenolicum. En: Applied microbiology and biotechnology. 44 (1995), p. 519525 [147] Zaher, U. ; Grau, P. ; Benedetti, L. ; Ayesa, E. ; Vanrolleghem, P.: Transformers for interfacing anaerobic digestion models to pre- and post-treatment processes. En: Water Environment Research 81 (2005), p. 407415 [148] Zaher, U. ; Rongping, L. ; Jeppsson, U. ; Steyer, J. ; Chen, S.: GISCOD: General Integrated Solid Waste Co-Digestion model. En: Water Research 43 (2009), p. 27172727
Declaraci on
Me permito armar que he realizado la presente tesis de manera aut onoma y con la u nica ayuda de los medios permitidos y no diferentes a los mencionados en la propia tesis. Todos los pasajes que se han tomado de manera textual o gurativa de textos publicados y no publicados, los he reconocido en el presente trabajo. Ninguna parte del presente trabajo se ha empleado en ning un otro tipo de tesis. Bogot a, D.C., dd.mm.aaaa
(Nombre del autor)

Implemetacion Del Modelo Adm1 en Matlab y Simulink

Cargado por

Copyright:

Formatos disponibles

Implemetacion Del Modelo Adm1 en Matlab y Simulink

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Implemetacion Del Modelo Adm1 en Matlab y Simulink

Cargado por

Copyright:

Formatos disponibles

6.

Implementaci on del modelo ADM1 en Matlab/Simulink

6 Implementaci on del modelo ADM1 en Matlab/Simulink

6.2 Desarrollo de la simulaci on del modelo ADM-1

Desarrollo de la simulaci on del modelo ADM-1

Validaci on del modelo respecto al caso base ADM-1

6 Implementaci on del modelo ADM1 en Matlab/Simulink

Funciones de inhibici on debido al H2 y Ninorganico

Funciones de inhibici on debido al NH3 y pH

Concentraciones iniciales de las especies solubles e insolubles

Calculo de la concentraci on de iones en la ecuaci on de equilibrio de cargas

Par ametros bioqu micos, sicoqu micos y de operaci on del reactor

7 estados referentes a los grupos de microorganismos (Xi )

6.4 Vericaci on del modelo respecto al caso experimental

Vericaci on del modelo respecto al caso experimental

6 Implementaci on del modelo ADM1 en Matlab/Simulink

6.4 Vericaci on del modelo respecto al caso experimental

Tabla 6-2.: Caracter sticas sicoqu micas residuos c tricos [119].

6 Implementaci on del modelo ADM1 en Matlab/Simulink

6.4 Vericaci on del modelo respecto al caso experimental

6 Implementaci on del modelo ADM1 en Matlab/Simulink

DEMANDA QU IMICA DE OXIGENO TOTAL (DQOT )

Demanda Qu mica de Ox geno Particulada (DQOP )

Demanda Qu mica de Ox geno Soluble (DQOS )

Degradables (63 % (DQOP ) )

Inertes (37 % (DQOP ))

Acidos Grasos Vol atiles (AGV)

Material Disuelto Excepto AGV

Carbohidratos (32 %) de Degradables

Proteinas (32 %) de Degradables

L pidos (36 %) de Degradables

Monosac aridos (33 %) de Disueltos no AGV

Amino acidos (33 %) de Disueltos no AGV AGCL (33 %) de Disueltos no AGV

Inertes Disueltos (18 %) de Inertes

Inertes Particulados (82 %) de Inertes

Acido ac etico (63 %) de AGV

Acido Propi onico (24 %) de AGV

Acido But rico (8 %) de AGV

Acido Val erico (5 %) de AGV

AGCL: Acidos grasos de cadena larga

pH, alcalinidad y acidos grasos vol atiles (AGV)

Raz on de carga 2:1

7.1 Resultados Experimentales

Demanda qu mica de ox geno soluble

8.00 7.75 pH [Unidades de pH] 7.50

(1:1-A) (2:1-B) (3:1-C)

(1:1) (2:1) (3:1)

6.25 6 0 5 10 15 20 25 30 35 40 Tiempo [d as]

(1:1) (2:1) (3:1)

500 y = 32.5x 0 0 5 10 15 20 25 Tiempo [d as] 30 35 40

7.2 Producci on de metano y biog as

Producci on de metano y biog as

Tasa m axima y Porcentaje de biodegradabilidad anaer obica

7.2 Producci on de metano y biog as

(1:1-A) (1:1-B) (2:1-A) (2:1-B) (3:1-A) (3:1-B)

Producci on de metano [mg DQO-CH4 /L]

(1:1) (2:1) (3:1)

(1:1) (2:1) (3:1)

7.2 Producci on de metano y biog as

Factor de conversi on de sustrato a metano