SciELO - Scientific Electronic Library Online

 
vol.16 issue32Carbonaceous Ore Characterization and the Identification of the Preg-Robbing EffectEl Volcán, a Magmatic Iron Deposit in NW Mexico: Distribution and General Characteristics author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Epistemus (Sonora)

On-line version ISSN 2007-8196Print version ISSN 2007-4530

Epistemus (Sonora) vol.16 n.32 Hermosillo Jan./Jun. 2022  Epub Sep 23, 2022

https://doi.org/10.36790/epistemus.v16i32.210 

Investigación

Modelamiento y simulación de una planta de recuperación de plata con tiourea

Modeling and Simulation of a Silver Recovery Plant with Thiourea

R. Robles-Ramos1 

M. J. Reyna-Robelo2 

D. Calla-Choque3  * 
http://orcid.org/0000-0003-2848-4388

J.S. Santos-Jallath4 

1Estudiante de Ingeniería de Minas y Metalurgia, Facultad de Ingeniería, Universidad Nacional Autónoma de México, Coyoacán, CDMX, México, ruben.robles@ingenieria.unam.edu

2Estudiante de Ingeniería de Minas y Metalurgia, Facultad de Ingeniería, Universidad Nacional Autónoma de México, Coyoacán, CDMX, México, mauricio67@comunidad.unam.mx

3Doctor en Ciencias en Ingeniería Metalúrgica y Cerámica, Departamento de Ingeniería de Minas y Metalúrgica - Facultad de Ingeniería, Universidad Nacional Autónoma de México, Ciudad de México, México, dcalla@unam.mx

4Doctor en Ingeniería, Departamento de Ingeniería de Minas y Metalúrgica - Facultad de Ingeniería, Universidad Nacional Autónoma de México, Ciudad de México, México, jesantos@unam.mx


RESUMEN

Algunos residuos generados en operaciones metalúrgicas tienen elementos metálicos de interés, por lo que pueden ser reprocesados para darle un valor adicional a la operación. Este trabajo presenta el modelamiento y simulación de tanques de lixiviación para la recuperación de plata con tiourea y oxalato a partir de jarositas que son uno de los principales desechos del procesamiento del zinc. Se usaron modelos de primer orden y se seleccionó el de mejor ajuste para la evaluación de la velocidad de lixiviación en tanques con agitación continua (CSTR). Finalmente se presentan los resultados obtenidos en la simulación del proceso con ayuda del software Matlab-Simulink® por medio de bloques de subsistemas.

Palabras clave: Modelamiento; simulación; plata; tanques de lixiviación

ABSTRACT

Some metallurgical wastes have residual metals of interest generated in metallurgical operations have values ​​of interest that can be reprocessed to give additional value to the operation, this work presents the modeling and simulation of leaching tanks for silver recovery with thiourea and oxalate from jarosite that is principal zinc processing waste for this purpose first-order models were used and the one with the best fit was selected for the evaluation of the leaching rate in tanks with continuous stirring (CSTR). Finally, the results obtained with the process simulation are presented with the help of the Matlab-Simulink® software using subsystem blocks.

Key words: Modeling; simulation; silver; leach tanks

Introducción

Las jarositas son de interés en la industria metalurgia del zinc porque su formación facilita la eliminación de hierro en las soluciones de lixiviación, mientras genera un residuo estable y fácilmente removible del sistema [1]. Además, puede eliminar algunos otros elementos como el plomo (Pb), arsénico (As), cadmio (Cd), plata (Ag) entre otros.

La jarosita con contenido de Ag, usada en este artículo es también llamada argentojarosita o jarosita de Ag, su contenido de plata resulta atractiva por el incremento de los precios de los metales; por lo cual se estudia una alternativa para la recuperación de metales valiosos [2].

Tiourea

La tiourea ha sido propuesta como una alternativa para el tratamiento de minerales complicados en la recuperación de metales como, los sulfurosos, minerales que consumen cianuro (cianicidas) y para lugares donde las preocupaciones ambientales hacen que el uso del cianuro sea restringido. La tiourea es un reactivo relativamente no tóxico, que se comporta como fertilizante para algunas plantas, y podría ser una alternativa atractiva al cianuro en ambientes sensibles.

Por otra parte, se sospecha que es cancerígeno y es capaz de disolver metales pesados, además del Au y Ag, presenta algunos problemas ambientales similares al cianuro en el manejo y disposición de efluentes, sin embargo, su descomposición y tratamiento de las soluciones con tiourea es mucho más sencilla por métodos químicos [2].

Modelos cinéticos

Los modelos matemáticos son un instrumento muy común en el estudio de sistemas de toda índole. Para el caso que se estudia la recuperación de elementos o especies valiosas, un parámetro de importancia es la velocidad de la recuperación, por lo que se busca un modelo matemático para la evaluación y el análisis del proceso de concentración o lixiviación, así como para la optimización y simulación [3].

El primer modelo cinético presentado es el de García Zúñiga (1), se basa en el supuesto que la velocidad de colisión de partículas es de primer orden con respecto al número de partículas presentes y la concentración de la solución en el sistema permanece constante [4].

La ecuación cinética se expresa de la siguiente forma [5]:

Rt= Rmax(1-e-kt) (1)

El segundo modelo cinético presentado es el Rápido-Lento (2), matemáticamente este modelo se describe como:

Rt=Rmax[(θr(1-e-krt)+θl1-e-klt]

θr+θl=1 (2)

El tercer modelo cinético presentado es el de Klimpel (3), su ecuación es la siguiente:

Rt=Rmax1-1kt1-e-kt (3)

Donde R es la recuperación del elemento útil, k es la constante de velocidad de primer orden, t es el tiempo de lixiviación, Rmax = Recuperación máxima [6].

El cuarto modelo cinético presentado es el de recuperación de orden n, linealizado, donde CA representa la concentración de plata en solución (4).

Log CA-Log CA0=-kt (4)

Finalmente, el modelo cinético presentado es el de modelo de núcleo decreciente (5) que involucra un proceso cinético controlado por reacción química, además describe la velocidad de disolución a diferente tamaño de partícula en la constante de velocidad, k [7].

Su ecuación es la siguiente:

1-[1-X]13=kt (5)

Para el presente trabajo se desarrollaron los cinco modelos descritos para hacer el análisis y simulación con los datos experimentales.

Metodología

Los siguientes datos experimentales de recuperación de plata a partir de jarositas industriales fueron obtenidos con las pruebas de lixiviación con tiourea y oxalato [1] bajo las siguientes condiciones, pH=1 a una temperatura de 50°C:

Tabla 1: Datos experimentales de la recuperación de plata en jarositas con tiourea. 

Tiempo (h) 0.00 0.17 0.33 0.50 1.00 1.50 2.00
%R (Ag) 4.18% 43.92% 53.80% 60.63% 75.95% 88.80% 94.62%

Posteriormente se usaron los diferentes modelos cinéticos para después comparar cuál presenta un mejor ajuste a los datos experimentales. Los modelos considerados para este fin son los siguientes modelos de primer orden, García Zúñiga, Rápido-Lento, Klimpel, ecuación de orden n y núcleo decreciente.

Con el modelo cinético de mejor ajuste se determinó el número de tanques y el volumen de los tanques (6) de agitación continua (CSTR), con una concentración inicial CA0=158gt Ag y un F0=30m3h .

Vn=CA0*Flujo pulpa*(Rn-Rn-1)-rAn (6)

Además, se procedió a determinar el volumen promedio que debe tener cada tanque (7).

V=V1+V2+V3+V4+V55=VT5 (7)

Con todos los datos obtenidos de las expresiones anteriores, se procedió a la simulación en el software Matlab-Simulink®, para ello es necesario determinar concentración CA , se partió de la ecuación (8) y de ahí la ecuación diferencial (9).

A=F0CA=Vd(CA)dt (8)

dCA1dt=F0CA0-F1CA1V1-kCA1n (9)

Llegado a este punto, en el reactor CSTR se obtiene una concentración en el primer tanque con la ecuación (9) y para el cálculo de la recuperación, la ecuación (10).

%R=CA0-CA1CA1 (10)

Colocando un segundo reactor en serie, el modelo para el segundo reactor se presenta en la ecuación (11) y la recuperación la ecuación (12).

dCA2dt=F1CA1-F2CA2V (11)

%R=CA0-CA2CA2 (12)

Y se procede a hacer el mismo calculo para el tercer, cuarto y quinto tanque.

Por último, se determina el promedio de los volúmenes de cada tanque (12).

V=V1+V2+V3+V4+V55 (13)

Resultados

Se presentan los resultados de la elección del modelo cinético y ajuste de datos de recuperación de plata a partir de jarositas con tiourea y oxalato (tabla 2).

Tabla 2: Modelos cinéticos en la recuperación de plata a partir de jarositas con tiourea. 

  Rápido-Lento Garcia-Zuñiga Klimpel Orden n Núcleo decreciente
Rmax 0.888 0.888 0.981 - -
K - 0.047 0.059 0.013 0.050
θr,θl 0.500 - - - -
kr,kl 0.047 - - - -
θr + θl 1.000 - - - -
Error total 0.027 0.027 0.016 0.016 0.012

Del ajuste de los resultados experimentales a los cinco modelos presentados, rápido-lento, Garcia-Zuñiga, orden n y núcleo decreciente se obtuvieron los errores al ajuste de 0.027 rápido-lento y Garcia-Zuñiga, 0.016 Klimpel y Orden, y 0.012 para el modelo de núcleo decreciente por lo que este modelo muestra una mejor descripción de los resultados experimentales obtenidos bajo las condiciones experimentales para la lixiviación de la plata (Figura 1) con un orden n=1   y con una constante de velocidad, k=0.05.[/p]

Figura 1: Modelo del núcleo decreciente.  

Sustituyendo en la ecuación (6), se obtuvieron los volúmenes de cada tanque y con la ecuación (9) se determina la concentración en cada tanque (ver tabla 3).

Tabla 3: Obtención de volumen y concentración en cada tanque. 

Tanque Volumen[m3]
dCAdt
Ag, CA
gt
TK-1 809
30 m3h158 gt-30 m3h(CA1)287.3 m3-0.05(CA1)
101.3
TK-2 287
30 m3h101.3 gt-30 m3h(CA2)287.3 m3-0.05(CA2)
65.04
TK-3 139
30 m3h65.04  gt-30 m3h(CA3)287.3 m3-0.05(CA3)
42.11
TK-4 149
30 m3h42.11  gt-30 m3h(CA3)287.3 m3-0.05(CA4)
27.9
TK-5 52 
30 m3h27.9 gt-30 m3h(CA5)287.3 m3-0.05(CA5)
19.4

El volumen promedio de cada reactor que se obtuvo de la ecuación (12) para procesar 30 m3h pulpa, con una ley de oro de 158 gt Ag.

V=1437 [m3]5=287.30  [m3]

Las recuperaciones en cada tanque se presentan en la tabla 4 a un tiempo de 47.88 horas:

Tabla 4: Recuperación de plata en jarositas con tiourea y oxalato en g/t por tanque. 

Recuperación de Plata (%)

TK-1

TK-2

TK-3

TK-4

TK-5

35.92

58.94

73.35

82.34

87.72

Los resultados que se obtuvieron de la simulación del proceso con ayuda del software Matlab-Simulink® en bloques de subsistemas se presentan en la figura 2 y la gráfica de recuperación en cada tanque con respecto al tiempo se muestra en la figura 3.

Figura 2: Modelamiento y simulación en el software Matlab-Simulink® 

Figura 3: Gráfica de las recuperaciones en cada tanque 

En la figura 2, se muestra parte de la simulación obtenida con ayuda del software Matlab-Simulink con un flujo constante en cada tanque de 30 m3h .y una concentración inicial de 158 gt de Ag; en el primer tanque la concentración de plata que se obtiene es de 101.3 gt , el segundo tanque tiene una concentración 65.04 gt , el tercer tanque 42.11 gt , el cuarto 27.9 gt y el último tanque permite obtener una concentración de 19.4 gt .

Para las recuperaciones obtenidas en cada tanque, la figura 3 muestra el porcentaje de lixiviación de plata con tiourea y oxalato a partir de jarositas en un tiempo de 47.88 horas. En el primer tanque se obtiene una recuperación del 35.92%, el segundo 58.94%, el tercero 73.35%, el cuarto de 82.34% y el último tanque es de 87.92%.

Conclusión

De los cinco modelos presentados, el modelo con un mejor ajuste para la lixiviación con tiourea y oxalato de las jarositas para obtener plata, bajo una temperatura 50 [°C] y un pH=1 fue el de núcleo decreciente, debido a que representa una mejor descripción de los resultados experimentales con un error total de 0.012 para una constante de velocidad de 0.05.

Se determinó que con flujo constante igual a 30m3h en los 5 reactores CSTR se obtiene una recuperación del 87.72% de plata a partir de las jarositas con contenidos de plata, equivalente a 175.54 ozdía .

Se presenta la simulación de una planta de agitación en tanques de lixiviación mediante el software Matlab-Simulink®, así como la representación gráfica del modelado de la recuperación en cada tanque con respecto al tiempo.

Reconocimientos

Investigación realizada gracias al Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica (PAPIIT) de la Universidad Nacional Autónoma de México IA 106221.

Referencias

D. Calla-Choque, G.T. Lapidus, “Jarosite dissolution kinetics in the presence of acidic thiourea and oxalate media”, Hidrometallurgy, vol. 206, pp. 1-6, January 2021. [ Links ]

D. Calla-Choque , “ Tratamiento de los residuos del proceso jarosita de la industria metalurgica del zinc, con la finalidad de mitigar este pasivo ambiental”, Tesis de maestria, Universidad Nacional de Ingenieria, UNI, Lima, Perú, 2012. [ Links ]

M. Alvarez-Silva, L. Vinnett , R. Langlois, K.E. Waters, “A comparison of the predictability of batch flotation kinetic models”, Mineral Engineering , vol 174, August 2016. [ Links ]

O. Peña i Lillo, “Principios de Flotación Selectiva aplicados a minerales complejos”, Boletin minero de la Sociedad Nacional de Minería, vol. 47, pp.113, Febrero de 1935. [ Links ]

R. Klimpel, “Selección de reactivos químicos para flotación”, SME, vol.33 , pp16, Enero 1980. [ Links ]

T. Coleman, M.A. Branch, A. Grace.Optimization Toolbox. 2nd ed. USA:The Math Works, Inc.1996. [ Links ]

K.C. Lillell “Shrinking core models in hydrometallurgy: What students are not being told about the pseudo-steady approximation”, Hidrometallurgy, vol. 79, pp. 62-68, September 2005. [ Links ]

Cómo citar este artículo:

Robles Ramos, R., Reyna Robledo, M. J., Calla Choque, D., & Santos Jallath, J. E. (2022). MODELAMIENTO Y SIMULACIÓN DE UNA PLANTA DE RECUPERACIÓN DE PLATA CON TIOUREA. EPISTEMUS, 16(32). https://doi.org/10.36790/epistemus.v16i32.210

Recibido: 02 de Noviembre de 2021; Aprobado: 25 de Marzo de 2022; Publicado: 29 de Abril de 2022

*Autor de Correspondencia: Dandy Calla Choque, dcalla@unam.mx

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons