INTRODUCCIÓN
Las principales fuentes de contaminación ambiental con metales pesados son los residuos procedentes de los vertidos industriales. Una vez liberados al ambiente, los metales pesados circulan entre los ciclos bióticos y abióticos, se acumulan en diferentes compartimentos de la cadena trófica o en el ambiente, y pueden alcanzar concentraciones tóxicas para animales, plantas, microorganismos e inclusive el hombre (Akhtar et al. 2004, Özer et al. 2008).
Los métodos utilizados para el tratamiento de efluentes con metales pesados, que incluyen precipitación química, tratamiento electroquímico, intercambio iónico, adsorción y tecnologías de membrana, entre otros, la mayoría de las veces son inviables técnica y/o económicamente, en especial cuando los metales disueltos están contenidos en grandes volúmenes de agua y en concentraciones relativamente bajas (Gupta et al. 2012, Imran 2012). De esta forma, la mayoría de las investigaciones se han orientado a la búsqueda de nuevos procesos y materiales, entre ellos la bioadsorción, una de las tecnologías más prometedoras no sólo por su bajo costo sino porque posibilita el tratamiento de grandes volúmenes de agua con bajas concentraciones de metal en forma eficaz (Ali 2014). Además, la posibilidad de emplear biomasa muerta o productos derivados de su metabolismo supera problemas de toxicidad e incluso permite la regeneración y reutilización del biomaterial por varios ciclos de adsorción/desorción (Gupta et al. 2012, Imran 2012).
La mayor parte de los estudios sobre bioadsorción que aparecen en la literatura se han realizado en sistemas discontinuos, aunque desde un punto de vista práctico los procesos de bioadsorción a gran escala se llevan a cabo de un modo continuo, generalmente en columnas de lecho fijo. Para diseñar unidades efectivas de lecho empacado con adsorbente sólido, se requieren modelos matemáticos que puedan simular satisfactoriamente las curvas de ruptura obtenidas experimentalmente (Ali et al. 2012). Estos modelos deben predecir la dinámica del proceso de adsorción o de intercambio iónico para facilitar el desarrollo de aplicaciones novedosas, tales como el diseño de columnas empacadas con material biosorbente para remover metales pesados de aguas residuales industriales (Alomá et al. 2013, Fallah et al. 2015).
Se han reportado diversos estudios en sistemas discontinuos que utilizan cáscara de cacao para la remoción de Pb+2 y Cd+2 en solución acuosa, los cuales han determinado la influencia del pH de la disolución, la dosis de adsorbente, el tiempo de equilibro y los parámetros termodinámicos y cinéticos del equilibrio de adsorción (Olu-owolabi et al. 2012, Odoemelam et al. 2011). Sin embargo, no se han reportado estudios en columnas utilizando la cáscara de cacao como bioadsorbente, por lo que el objetivo del presente trabajo fue la remoción de Pb+2 y Cd+2 de soluciones sintéticas en columna de lecho fijo y ajustar los datos experimentales a varios modelos que describen la curva de ruptura para sistemas mono y bicomponentes en función del volumen del efluente.
MARCO TEÓRICO
Bioadsorción de metales pesados en sistemas continuos
El sistema continuo más estudiado es el de columna de lecho fijo, ya que su implementación a escala industrial es más factible (Yan y Viraraghavan 2001, Khoramzadeh et al. 2013). El funcionamiento y la eficiencia de las columnas de lecho fijo se describen mediante el concepto de curva de ruptura, que resulta al graficar la concentración normalizada (C/Co) en función del tiempo o del volumen tratado, siendo éste el menos utilizado en investigaciones recientes. El tiempo de operación y la curva de ruptura son características muy importantes para determinar la respuesta de una columna de bioadsorción y extrapolar los datos a escala industrial (Bhatnagar y Minocha 2010, Rodríguez et al. 2014 ).
Parámetros de funcionamiento de la columna
Volumen del efluente
donde Q es el caudal que circula por la columna (mL/min) y ttotal el iempo total (min).
Capacidad máxima de la columna
El área bajo la curva de ruptura, entre los límites apropiados, representa la cantidad total de metal retenido, q total (mm), para una determinada concentración de la alimentación, y puede determinarse por integración:
donde C R es la concentración de metal retenido (mg/L).
Cantidad total de metal que pasa por la columna
Porcentaje total de metal retenido durante la operación
Longitud de lecho no utilizada
El cálculo de la longitud de lecho no utilizada constituye un método para evaluar la capacidad de adsorción de bioadsorbentes en columnas de flujo continuo:
donde H es la altura del lecho (cm), t s el tiempo en que se alcanza la saturación del bioadsorbente (min), y t r el tiempo en que se alcanza el punto de ruptura (min).
Modelado matemático del proceso de bioadsorción en columnas
Para conocer la capacidad de adsorción de una biomasa es importante describir la cinética del proceso. Con este fin existen varios modelos matemáticos, los cuales son útiles para el diseño experimental y suministran las predicciones para la ejecución del proceso de bioadsorción en cuestión, en diferentes condiciones operativas.
Entre los modelos matemáticos más utilizados para describir el comportamiento dinámico de la eliminación de contaminantes en columna de lecho fijo en función del volumen del efluente están el de Thomas y el dosis-respuesta. Éstos se aplicarán a los datos experimentales obtenidos para bioadsorción de Pb+2 y Cd+2 con cáscara de cacao.
Modelo de Thomas
El modelo de Thomas es uno de los más generales y utilizados para describir el comportamiento del proceso de bioadsorción en columnas de lecho fijo. Este modelo se usa para conocer la máxima capacidad de adsorción de un adsorbente y obtener el rendimiento teórico de la columna.
El modelo se describe mediante:
donde K th es la constante de velocidad (mL/min mg), q o la capacidad de adsorción en el modelo de Thomas (mg/g), m B la masa del bioadsorbente (g), V ef el volumen del efluente (mL) y Q el caudal de alimentación (mL/min)
Modelo de dosis-respuesta
Este modelo ha sido comúnmente utilizado en la farmacología para describir diferentes tipos de procesos. En la actualidad se emplea para describir los procesos de bioadsorción en columna, y su importancia radica en que describe la curva completa de ruptura con gran exactitud, además de que puede minimizar los errores que arroja el modelo de Thomas, especialmente con tiempos de remoción bajos o altos.
En función del volumen del efluente,
Considerando
MATERIALES Y MÉTODOS
Materiales bioadsorbentes
La cáscara de cacao se obtuvo en el cantón de Portovelo, ubicado en las estribaciones de la cordillera occidental de los Andes, al sureste de la provincia de El Oro. La altitud del cantón varía de 600 a 3000 msnm y sus coordenadas geográficas son 9594200-9578900 lat. N, 680900-650500 long. E. La cáscara fue lavada, secada, triturada y tamizada, para determinar propiedades físicas química como densidad, porosidad, sitios ácidos y básicos, punto de carga cero. Se trabajó con partículas >1 mm.
Preparación de las soluciones
Las soluciones se prepararon a partir de nitrato de plomo y cadmio, adquiridos como estándares certificados de 1000 ppm trazables del National Institute of Standards and Technology (NIST). Se utilizó agua destilada para las diluciones correspondientes a pH 5, que garantiza la especiación iónica de los metales.
Análisis de infrarrojo
Se realizó la técnica de espectroscopia de infrarrojo con transformada de Fourier (FTIR) en un espectrofotómetro marca Perkin Elmer, modelo 20 SXB, para identificar los grupos funcionales.
Cada muestra fue sometida a un proceso de secado a 50-60 ºC por 24 h.
Para realizar el análisis de infrarrojo se prepararon pastillas de KBr, con fragmentos de cada una de las muestras trituradas en un mortero.
Se obtuvo el espectro por transmitancia realizando 25 barridos en la frecuencias de 4000 a 450 cm-1 con una resolución de 4 cm-1.
Pruebas hidráulicas y selección de los parámetros de operación
Los experimentos de bioadsorción se realizarons en dos columnas de vidrio 1.5 y 1.6 cm de diámetro interno, ambas con una altura de 50 cm. La preparación de las columnas se llevó a cabo como se describe a continuación:
En la parte inferior de las columnas se colocó gravilla con un diámetro de 2 mm con el objetivo de que las partículas del bioadsorbente no fueran arrastradas en el proceso de bioadsorción.
A continuación, se introdujo el bioadsorbente en el interior de las columnas, hasta la altura de lecho correspondiente a cada experimento.
Finalmente, la parte superior de la columna se rellenó con gravilla con un ϕ de 5 mm para evitar la flotación del lecho y la pérdida de material.
La operación del sistema se realizó con flujo descendente.
Los experimentos se ejecutaron a temperatura ambiente (17 ºC).
Cuantificación de las muestras
Una vez estudiada la bioadsorción de Pb2+ y Cd2+ en la columna de lecho fijo, se aplicaron varios modelos a los datos experimentales de las curvas de ruptura y la determinación de los parámetros. Las concentraciones de plomo y cadmio se determinaron mediante espectrometría de absorción atómica utilizando un espectrofotómetro Perkin Elmer AAnalyst 400.
RESULTADOS Y DISCUSIÓN
A continuación se presentan los resultados obtenidos en la bioadsorción de Pb+2 y Cd+2 en columnas utilizando la cáscara de cacao como bioadsorbente. El proceso inicia con un análisis infrarrojo para determinar los grupos funcionales del bioadsorbente.
Análisis de infrarrojo
La figura 1 muestra el espectro infrarrojo obtenido para la cáscara de cacao antes y después de la bioadsorción de los metales. Se puede observar la banda ancha entre las longitudes de 3000 a 3600 cm-1. En esta región aparecen vibraciones del grupo -OH, el más abundantes tanto en la celulosa como en el ácido péctico. El pico a 2915 cm-1 se atribuye a enlaces C-H alifáticos.
Pueden observarse una serie de señales superpuestas en la banda centrada en 1600 cm-1, en la cual se consideran las vibraciones de los grupos carbonilos (C=O) provenientes de ésteres, amidas y, en menor proporción, de ácidos carboxílicos presentes en las moléculas de ácido péctico. La banda centrada en 1477 cm-1 es la vibración estrecha de un ácido carboxílico dímero (el cual forma puentes de hidrógeno intermoleculares con otro ácido carboxílico) que puede estar presente en la celulosa. Finalmente, la banda centrada en 1052 cm-1 corresponde a la vibración del enlace -C-O.
Después de la bioadsorción de Pb+2 y Cd+2, el espectro ha cambiado poco. Se muestra una disminución en la intensidad de las bandas característica de los grupos éster alifáticos, lo que sugiere que hay una fuerte interacción de estos grupos funcionales con los metales.
Estudio hidrodinámico
Se realizó el estudio hidrodinámico en columnas para conocer el flujo a través de las mismas, variando la altura del lecho fijo, el diámetro de las columnas y el diámetro de las partículas (Cuadro I).
Cáscara de cacao | ||||
Diámetro de la columna (cm) | h (cm) | Masa de bioadsorbente (gramos) | Q (mL/min) | |
ɸ partícula > 1 mm | ɸ partícula > 0.594 mm | |||
1.5 | 10.5 | 0 | 5 | 0.34 |
1.5 | 10.5 | 5 | 0 | 0.72 |
1.5 | 10.5 | 5 | 0 | 0.72 |
1.6 | 10.5 | 5 | 0 | 1.25 |
1.6 | 15 | 7 | 0 | 2.12 |
Q: caudal de alimentación
El estudio hidrodinámico a escala de laboratorio permitió conocer que las mejores condiciones de operación, en las que no hay inundación ni fragmentación de la columna fueron las siguientes: diámetro de la columna = 1.6 cm; tamaño de partícula >1 mm (fracción gruesa); altura del lecho de 15 cm, y 7 g de biomasa para el relleno, con un flujo a través de la columna de 2.12 mL/min, recomendado por varios autores (Aksu y Gonen 2004, Calero et al. 2009, Ali 2014).
Se han llevado a cabo estudios para remoción de Pb+2 de soluciones acuosas sintéticas utilizando aserrín de Pinus sylvestris en una columna empacada, con las siguientes condiciones hidrodinámicas: 1.5 cm de diámetro, 15 cm de largo y 2.51 mL/min (Taty et al. 2005).
Ensayos de bioadsorción en columnas. Análisis de la curva de ruptura
Para analizar la remoción de Pb2+ y Cd2+ en columna de lecho fijo se graficó la relación de concentración final e inicial (C/C0) en mg/L vs. el tiempo (t), siendo las concentraciones iniciales de 10 mg/L y pH de 5. Los resultados se ilustran en las figuras 2 y 3 para sistemas mono y bicomponentes, respectivamente.
Para los sistemas mono y bicomponentes, el área bajo la curva es mayor para el Cd+2 que para el Pb+2, por lo que existe mayor afinidad de los sitios activos por el Pb+2.
En la mayoría de los estudios de bioadsorción en sistemas con mezclas de metales se ha encontrado que los grupos funcionales presentan mayor afinidad por unos iones que por otros, lo que parece estar relacionado, entre otros factores, con determinadas características iónicas del metal como el radio iónico, la electronegatividad o el índice covalente que relaciona ambos parámetros (Liu et al. 2012,Senthilkumar et al 2010). En general, los bioadsorbentes suelen presentar mayor afinidad por los iones con mayor radio iónico o índice covalente y se produce una mayor competitividad entre iones con características semejantes (Wang y Chen 2009, Vijayaraghavan y Balasubramanian 2013). En este caso, el Pb+2 posee mayor radio iónico y electronegatividad que el Cd+2.
En el cuadro II se muestran los resultados de los parámetros obtenidos con la curva de ruptura en ambos sistemas.
Parámetros | Sistema monocomponente | Sistema bicomponente | |||
Pb+2 | Cd+2 | Pb+2 | Cd+2 | ||
Curvas de ruptura | tr (min) | 328 | 306 | 254 | 232 |
ts (min) | 643 | 522 | 641 | 463 | |
V ef (mL) | 629.76 | 587.52 | 624.84 | 570.72 | |
q o (mg/g) | 1.181 | 0.787 | 0.837 | 0.765 | |
LNU (cm) | 4.34 | 12.69 | 9.06 | 7.48 | |
M total (mg) | 15.06 | 15.06 | 19.29 | 19.29 | |
q t (mg) | 13.71 | 13.58 | 17.01 | 17.4 | |
% adsorción | 91 | 90 | 88 | 90 |
t r : tiempo de ruptura, t s : tiempo de saturación, V ef : volumen del efluente, q o : capacidad de adsorción, LNU: lecho no utilizado, M total : cantidad total de metal que pasa por la columna, q t : cantidad total de metal retenido
En sistemas mono y bicomponentes se logran resultados superiores con el Pb+2 en cuanto a los parámetros, tiempos de ruptura y saturación de la columna (t r y t s ), volumen tratado al alcanzar el punto de ruptura (V r ), y la cantidad de metal retenida por gramos de bioadsorbente (q t ).
El lecho no utilizado (LNU) para el Pb+2 es menor que para el Cd+2 en sistemas monocomponentes, a diferencia de los sistemas bicomponentes. El porcentaje de remoción en ambos sistemas es muy similar para el Pb+2 y el Cd+2.
Modelos aplicados en función del volumen del efluente en sistema monocomponentes. Modelos de Thomas y dosis-respuesta
De todos los modelos que se mencionan en la metodología, sólo el de Thomas y el de dosis-respuesta se encuentran en función del volumen del efluente. Es recomendable destacar la importancia de predecir las curvas de ruptura en función del volumen del efluente, ya que para fines prácticos permite conocer la cantidad de volumen del efluente que se puede tratar antes de iniciar la operación. En las figuras 4 y 5 se presenta la modelación de las mismas.
El modelo de dosis-respuesta se ajusta adecuadamente a los datos experimentales de adsorción de Pb+2 y Cd+2 en toda la curva, a diferencia del modelo de Thomas, donde se ve una variación en la curva que no se ajusta en todo el experimento realizado, con valores de coeficientes de determinación (R2) para dosis-respuesta de 0.98 para el Pb+2 y 0.99 para el Cd+2.
En el cuadro III se muestran los valores obtenidos de la constante de velocidad del modelo de Thomas (K th ) y la capacidad de adsorción (q o ), que son mayores para el Pb+2 que para el Cd+2, lo cual ratifica la afinidad del Pb+2 por los sitios activos del cacao. En el modelo de dosis-respuesta la capacidad de adsorción es mayor para el Pb+2 que para el Cd+2; sin embargo, los coeficientes de determinación (R2) para el Cd+2 muestran mayor relación debido a su ajuste a la forma de la curva de ruptura.
Modelos aplicados en función del volumen del efluente en sistemas bicomponentes. Modelos de Thomas y de dosis-respuesta
Para sistemas bicomponentes los modelos de Thomas y de dosis-respuesta de aplicaron por separado para apreciar mejor los ajustes de las curvas de ruptura modeladas con los resultados experimentales.
La figura 6 ilustra las curvas de ruptura experimentales de Pb+2-Cd+2 con las respectivas curvas del modelo de Thomas, donde se evidencia que no hay un buen ajuste.
Como se observa en el cuadro IV, la constante de velocidad del modelo de Thomas (K th ) es mayor para el Cd+2 que para el Pb+2, en tanto que la capacidad de adsorción (q o ) es mayor para el Pb+2. En la figura 7 se comparan las curvas de ruptura experimentales del Pb+2-Cd+2 con las respectivas curvas de ajustes obtenidas del modelo de dosis-respuesta (en función del volumen del efluente), mientras que en la cuadro V se encuentran los valores de los parámetros encontrados del modelo.
Parámetros | Pb+2 | ||
Pb+2 | Cd+2 | ||
Thomas | K th (L/min mg) | 1.1228 | 1.8612 |
q o (mg/g) | 1.458 | 1.254 | |
R2 | 0.96 | 0.97 |
K th : constante de velocidad, q o : capacidad de adsorción, R2: coeficiente de determinación
Parámetros | Metal | ||
Pb+2 | Cd+2 | ||
Dosis-respuesta | a | 3.964 | 6.1284 |
q o (mg/g) | 1.21 | 1.03 | |
R2 | 0.95 | 0.97 |
a: constante del modelo, q o : capacidad de adsorción, R2: coeficiente de determinación
En la figura 7 se observa que el modelo de dosis-respuesta reproduce la curva de ruptura completa con valores de R2 de 0.95 y 0.97 para el Pb+2 el Cd+2, respectivamente. Varios investigadores han obtenido este ajuste a los modelos de dosis-respuesta en la adsorción del Cu+2 con corteza de pino y Cr 6+ con bagazo de caña de azúcar (Calero et al. 2009, Rodríguez et al. 2014).
CONCLUSIONES
El residuo de cáscara de cacao es un material de bajo costo que se propone como alternativa para la tecnología de adsorción de Pb+2 y Cd+2 en procesos de tratamiento de aguas residuales utilizando columnas de lecho fijo.
En el estudio de los sistemas monocomponentes para Pb+2 y Cd+2 existe un porcentaje de remoción de 91 % para Pb+2 y de 90 % para Cd+2, mientras que en sistema bicomponentes (Pb+2-Cd+2) la remoción fue de 88 % para Pb+2 y de 90 % para Cd+2.
El modelo de dosis-respuesta en sistemas mono y bicomponentes fue el que mejor reprodujo las curvas de ruptura experimentales en todo el rango medido. A través de los parámetros encontrados con este modelo es posible obtener una expresión que reproduzca el comportamiento de las columnas en otras condiciones experimentales, sin necesidad de realizar experimentos adicionales, manteniendo las mismas condiciones de operación del sistema de lecho fijo.