Tecnología Química

ISSN: 0041-8420
revista.tec.quimica@fiq.uo.edu.cu
Universidad de Oriente
Cuba

Hernández-Pedrera, Carlos; Rivera-Soto, Margarita; Matos-Durán, Rafael; Piñó-Cuenca,

Yanetsis
Estudio de sensibilidad paramétrica en reactores continuos con agitación
Tecnología Química, vol. XXXIV, núm. 1, enero-abril, 2014, pp. 54-61
Universidad de Oriente
Santiago de Cuba, Cuba

Disponible en: http://www.redalyc.org/articulo.oa?id=445543781006

Cómo citar el artículo

Número completo
Sistema de Información Científica
Más información del artículo Red de Revistas Científicas de América Latina, el Caribe, España y Portugal
Página de la revista en redalyc.org Proyecto académico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
 Fecha recepción: sept. 2013; aceptación/publicación: dic. 2013, págs. 54-61

Estudio de sensibilidad paramétrica en reactores

continuos con agitación
Parametric Sensitivity Study on Continuous Reactors with Stirring
Dr. Carlos Hernández-PedreraI, pedrera@fiq.uo.edu.cu, Dra. Margarita Rivera-SotoI,
Dr. Rafael Matos-DuránI , Ing. Yanetsis Piñó-CuencaII
I
Facultad de Ingeniería Química, Universidad de Oriente, IIEmpresa Pedro Sotto Alba, Moa, Holguín, Cuba
En este trabajo se presentan los resultados obtenidos en un estudio de sensibilidad, utilizando un
modelo matemático desarrollado para la simulación del comportamiento de un reactor continuo
con agitación, utilizando como fuente de información datos de operación procedentes de la planta
industrial. El estudio permite valorar el efecto que los cambios en las variables de operación
pueden ocasionar en los resultados del proceso y la posibilidad de que existan o no interacciones
entre las variables analizadas.
Palabras clave: reactores continuos, sensibilidad paramétrica, modelo estadístico, simulación.
_____________________

In this work present the results obtained in a study of sensibility, by using a mathematical model
developed for the simulation of the conduct of an endless reactor with agitation, by using as data
source of information of reasonable operation of the industrial plant. The study permits value the
effect that the changes in the variables of operation can occasion in the results of the process and
the possibility that exists or not interactions between the variables analyzed.
Key words: continuous reactors, parametric sensibility, simulation, statistical model.

Introducción En el presente trabajo el desarrollo matemático

del modelo estadístico, tiene como finalidad
El presente trabajo tiene como bases los presentar el estudio de sensibilidad paramétrica
aspectos teóricos relacionados con el desarrollo realizado únicamente con fines docentes y tomando
de un modelo matemático tipo caja negra, este elemento en consideración se han utilizado
desarrollado con la finalidad de describir la datos de la operación industrial que no
operación de los reactores utilizados en el proceso corresponden al período actual.
de precipitación de sulfuros de níquel y cobalto, a
partir de datos procedentes de la industria. En el
Fundamentación teórica
desarrollo de modelos matemáticos puede ser aplicado
el análisis de sensibilidad paramétrica con diferentes
propósitos, entre los que se pueden citar: la estimación Características del objeto de estudio
de parámetros, discriminación de modelos,
optimización y/o control, entre otros /3/. Las autoclaves de precipitación de sulfuros
La aplicación de métodos estadísticos como el son reactores continuos con agitación; son equipos
análisis de regresión para identificar el modelo cilíndricos situados de forma horizontal, preparados
matemático tipo caja negra que mejor reproduzca el para recibir presiones no muy elevadas y ambiente
comportamiento de los reactores (autoclaves) de moderadamente ácido. Están divididos en tres
precipitación de sulfuros de Ni y Co, resulta de gran secciones: A, B y C, por dos tabiques de tal de
utilidad para realizar el estudio de la sensibilidad del manera que la masa reaccionante pasa de una
sistema, que permita determinar la influencia de las sección a otra por reboso. Los compartimentos (C
variables de operación que tienen un mayor efecto y B) están provistos de agitadores de tipo axial, ya
sobre la calidad del proceso, ya que sus resultados que en ellos se requiere la mayor absorción posible
servirían como base preliminar para el desarrollo de del H2S (g); mientras que el (A) posee uno de tipo
modelos matemáticos más robustos con un mayor radial, debido a que en esta sección la salida de la
nivel de descripción que permitirían realizar estudios pulpa es por el fondo y la suspensión de las
más detallados del proceso. partículas deja de ser un problema.

54 TECNOLOGÍA QUÍMICA Vol. XXXIV, No. 1, enero-abril, 2014

 Por el compartimiento (C) entra el licor con un sulfuros de níquel y cobalto se encuentran: la
flujo que es regulado automáticamente, mientras temperatura, el pH, el grado de agitación, la
que el H2S(g) es introducido al autoclave por el presión parcial de H2S y las semillas.
compartimiento (A), y al ponerse en contacto con
Para determinar el efecto que estas variables
la solución se producen las reacciones en las que
pueden tener sobre el proceso, es necesario hacer
se obtienen los sulfuros de níquel y cobalto, esto
un estudio donde se manifieste su papel en el
provoca que la solución se sobresature y comience
el proceso de precipitación en forma de núcleos de desarrollo de las reacciones /1/. A continuación
pequeño tamaño, incrustando las partes fijas y móviles se analiza Ia influencia de diferentes variables de
del equipo y haciendo crecer las partículas que trabajo sobre la operación de los reactores utilizados
fueron alimentadas como semilla, produciendo una en el proceso:
pulpa en solución ácida de alrededor de 2,5 % de
sólido a la salida de la secuencia de reactores. En Efecto de la temperatura
esta etapa se produce del 90 al 97 % de la reacción.
Las reacciones que se realizan en estos En este sistema el efecto de la temperatura es
reactores son representadas a continuación: el resultante de sus efectos sobre dos fenómenos
distintos: la absorción del gas en la fase líquida y
H 2 S ( gas) → H 2 S (aq) (1) la reacción química. Un aumento de temperatura
produce un aumento en la velocidad de la reacción,
H 2 S (aq) → H + (aq) + HS − (aq) (2) pero al mismo tiempo, disminuye la solubilidad de
HS − (aq). → H + (aq) + S −2 (aq) (3)
los gases; este efecto negativo del aumento de
temperatura puede ser contrarrestado aumentando
Ni +2 (aq) + S −2 (aq) → NiS(crist) (4) la presión del sistema /1/. Por esta razón, el efecto
resultante en el sistema pudiera no resultar
Zn +2 (aq) + S −2 (aq) → ZnS(crist) (5) significativo al introducir pequeñas variaciones en
esta variable.
Co+2 (aq) + S −2 (aq) → CoS(crist) (6)
Por otro lado, un efecto secundario al operar a
SO4−2 (aq) + H + (aq) → HSO4− (aq) (7) una elevada temperatura es que provoca que los
gases de purga salgan más caliente y con mayor
Fenómenos que están presentes en la operación
humedad, provocando problemas en el
de los reactores de precipitación.
deshumidificador y los compresores.
Los fenómenos que ocurren en estos reactores son:
Pedrera /1/ señala que en mucha literatura se
- Transferencia de cantidad de movimiento
estudian los procesos de cristalización en sistemas
- Transferencia de calor
con enfriamiento para lograr la sobresaturación
- Transferencia de masa
de la solución, o por evaporación del solvente con
- Reacción química
la misma finalidad; pero el proceso que se estudia
no se ajusta a ninguno de estos casos.
Influencia de diferentes variables de trabajo
sobre la operación de los reactores de
precipitación de sulfuros de Ni y Co Efecto del pH

Para lograr una eficiencia de precipitación de

Influencia de las condiciones de trabajo 90 % o más del Ni y Co se debe neutralizar el licor
sobre los procesos de precipitación ácido lixiviado, llevándolo a un pH de 2 a 3 antes
de la precipitación /1/. El aumento de la acidez
Entre los principales factores que tienen mayor inicial disminuye la solubilidad del H2S en el licor,
influencia en el proceso de precipitación de los aumentando la concentración residual del Ni.

TECNOLOGÍA QUÍMICA Vol. XXXIV, No. 1, enero-abril, 2014

55
 El uso de soluciones con valores de pH diferentes Identificación de un modelo físico adecuado
a los reportados para que las reacciones ocurran en y un modelo matemático conveniente
las mejores condiciones, podría provocar la ocurrencia
de otras reacciones no deseadas. Los modelos estadísticos como modelos
empíricos se basan en relaciones estadísticamente
Efecto del grado de agitación significativas entre las variables que involucra, a
diferencia de los fenomenológicos cuyas bases
La agitación es de gran importancia para descansan en las leyes que rigen los procesos. Sin
mantener los sólidos en suspensión, evitando la embargo el modelo estadístico se caracteriza por
tener un alto poder predictivo, ya que reproducen
acumulación de los mismos en los reactores, lo
el funcionamiento del sistema razonablemente
cual conduciría a un sobreconsumo de energía en
bien, a pesar de que no permiten saber por qué
los agitadores /1/. Con el aumento de la intensidad
funciona así.
de la agitación se mejora la solubilidad del H2S de
la fase gaseosa en la solución, incrementando el Las ecuaciones que describen un modelo
proceso de disociación y la formación de S2- lo estadístico no son físicamente o dimensionalmente
cual incrementa la velocidad de precipitación de consistentes ni universales, ya que en rigor sólo
los sulfuros /2/. son válidas para el contexto espacio-temporal en
el que se calibraron /2/. En estos modelos la falta
de conocimiento acerca del sistema se compensa
Efecto de la presión parcial de H 2S con datos de calidad y en cantidad suficiente de las
variables de entrada y de salida; a partir de estos
La concentración de un gas en un líquido es datos se puede construir un modelo empírico que
proporcional a la presión parcial (ley de Henry) permita, posteriormente, obtener las variables de
C H2S = K*P H2S, por lo que manteniendo una presión salida partiendo de un nuevo conjunto de valores
parcial de H2S elevada se garantiza la concentración para las variables de entrada. Se conocen también
necesaria para la reacción. Algunos autores plantean como modelos caja negra y su uso depende de la
que la variación de presión entre 344,7 y 689,3 kPa finalidad que se persiga en el análisis de procesos.
no ocasiona cambios significativos en la eficiencia
Tomando en consideración los elementos
de precipitación /1/.
anteriores, se realiza un estudio de simulación del
comportamiento de los reactores continuos con
Efecto de las semillas agitación empleados en el proceso de precipitación
de sulfuros de Ni y Co, utilizando un modelo tipo
Varios autores, según cita Pedrera /1/, plantean caja negra, desarrollado y validado con este fin;
el uso de semilla como un mecanismo para sus resultados sólo serán válidos para el período
favorecer la obtención de cristales de mayor en el cual fue realizado y puede ser utilizado como
tamaño lo que favorece su posterior separación. estudio preliminar para futuros trabajos en los
Sin la adición de semilla el producto que se forma cuales se analice el comportamiento de los
es extremadamente fino, siendo muy difícil la reactores de precipitación utilizando un modelo de
sedimentación y lavado, mientras que con el uso mayor complejidad o nivel de descripción, como
de éstas se logra disminuir considerablemente el es el caso de un modelo fenomenológico.
por ciento de finos reduciendo con ello las pérdidas
por suspensión en los sedimentadores.
Análisis de sensibilidad paramétrica
La semilla es de gran importancia en el
desarrollo del proceso, ya que ejerce un efecto El análisis de sensibilidad paramétrica es una
catalítico a bajas temperaturas y presión técnica que resulta de gran utilidad en el análisis
atmosférica, sin embargo éste se ve disminuido del comportamiento de cualquier proceso industrial,
para altas temperatura y presión. puede ser desarrollado a partir de un modelo

56 TECNOLOGÍA QUÍMICA Vol. XXXIV, No. 1, enero-abril, 2014

 matemático, que describa y permita simular el tiene solución única y continua en x y ö
comportamiento de los equipos utilizados en el representada por:
proceso /3, 5/. yo = y ( x , ö)
El estudio de sensibilidad paramétrica realizado
en un equipo permite al ingeniero determinar si los Índice de sensibilidad local de primer orden
límites de temperatura, presión u otras
restricciones, se sobrepasan o no cuando se Se identifica simplemente como "sensibilidad
producen perturbaciones de las condiciones límites local" de la variable. Con relación al parámetro de
o de las entradas, aun cuando el proceso pueda entrada ϕ j, se define por /3/:
ser estable /2/. En el mismo se puede analizar
tanto la sensibilidad de una misma respuesta del ∂ y ( x ,ϕ )
j
sistema con relación a cada parámetro del modelo, S ( y, ϕ ) =
como la sensibilidad de cada respuesta frente a un j ∂ϕ
j
mismo parámetro /3/.
Si el sistema tiene n parámetros, se tendrá
también n índices de sensibilidad. Esta técnica,
Sensibilidad paramétrica
aplicada en el análisis de un proceso industrial,
permite identificar los parámetros que ejercen el
Está dada por la sensibilidad que muestra el
mayor efecto sobre cada uno de los resultados de
comportamiento de un sistema ante los cambios
interés y da la posibilidad de seleccionar aquellos
en los valores de sus parámetros, siendo dos
que requieren una estimación rigurosa. El análisis
elementos a destacar los siguientes /3-5/:
de sensibilidad paramétrica puede ser empleado
- La sensibilidad del sistema no es la misma como una efectiva herramienta para identificar las
frente a los cambios en los diferentes variables de operación, cuyos cambios ejercen mayor
parámetros. efecto en el comportamiento del equipo o proceso
- La sensibilidad del sistema frente a cambios de analizado y para determinar la posibilidad de que
igual magnitud en un mismo parámetro depende ocurran interacciones entre las variables, en
además, del nivel que tenía el parámetro antes dependencia del nivel de operación en que se trabaje
de sufrir el cambio. en los reactores. Estas valiosas informaciones
permiten establecer acciones correctivas para el
Índices de sensibilidad control y la dirección eficaz del proceso.

En términos matemáticos la sensibilidad se

Métodos utilizados y condiciones
expresa mediante índices.
experimentales
Índice de sensibilidad local. Definición general
Para una variable se define por la expresión Como herramienta de trabajo principal se
general /3/ utilizaron los programas STATISTIC y EXCEL,
para el procesamiento de datos procedentes de la
dy industria (correspondientes a un período histórico
= f ( y, x , ϕ ) de la producción), la interpretación y análisis de
dx
los resultados.
y: variable dependiente
X: variable independiente Resultados y discusión
ö: vector que contiene los parámetros de entrada
al sistema Se identificaron dos modelos matemáticos que
f: es una función continua y diferenciable en correlacionan la eficiencia de precipitación de
todos sus argumentos, la cual para yo = y (o), sulfuros de Ni y Co como función de las variables

TECNOLOGÍA QUÍMICA Vol. XXXIV, No. 1, enero-abril, 2014

57
de operación de la planta; ambos con elevados Se seleccionaron dos corridas base: una para
coeficientes de regresión: 99,97 % para el caso el nivel mínimo y otra para el nivel máximo de
de la eficiencia de precipitación del níquel ( ENI ) operación; estos niveles fueron identificados dentro
y 99,93 % para la eficiencia de lixiviación del del intervalo de variación de las condiciones de
trabajo en una muestra formada por 183 juegos de
cobalto ( ECO ). Las pruebas de validación
datos y para todo el período estudiado, con la
arrojaron resultados muy satisfactorios, ya que
finalidad de analizar posteriormente la sensibilidad
ambos modelos los valores reportados por el
que tendría el sistema ante pequeños cambios
laboratorio y los predichos por el modelo,
introducidos en cada una de las variables de operación
discrepaban con un error menor a la unidad.
señaladas y la posibilidad de que existieran o no
Los resultados obtenidos permitieron interacciones entre algunas de ellas. A continuación
comprobar que los modelos son adecuados para se analiza el comportamiento observado.
realizar diferentes estudios de simulación, entre
ellos un análisis de sensibilidad paramétrica que Análisis de la influencia de pH, concentración
permita valorar el efecto que producirían los
cambios de algunas variables de trabajo sobre la
de manganeso y temperatura sobre la
operación de los reactores de precipitación. eficiencia de precipitación del níquel
El análisis se realizó como un estudio preliminar En las figuras 4, 5 y 6 se puede apreciar el
y como tal, teniendo en cuenta el tipo de modelo comportamiento del sistema ante los cambios en
desarrollado, se encaminó a analizar el efecto estas variables y el efecto de estas sobre la
que, sobre la eficiencia de precipitación del níquel eficiencia de precipitación del níquel. En la figura
y del cobalto, ocasionarían los cambios en las 1 se puede apreciar que iguales cambios
variables: pH, concentración de manganeso en el introducidos en el pH, ocasionan efectos de
licor de entrada (CMn(Le)) y temperatura. Se consideró diferente magnitud, en dependencia del nivel de
que el efecto de la presión parcial de H2S y de las presión con que se trabaje en la planta, lo que
semillas debe ser estudiado con un modelo que tenga demuestra que existen interacciones entre las
mayor nivel de descripción y robustez. variables analizadas.

Fig. 1 Efecto de los cambios de pH sobre la eficiencia precipitación del Ní.

En la figura 1 se puede notar que la eficiencia nivel e igualmente aumenta cuando se mantiene
de precipitación del níquel aumenta cuando se el pH constante en el máximo nivel; sin embargo
incrementa la presión desde (512,30 a 977,10) la magnitud del cambio que experimenta el por
kPa, manteniendo el pH constante en el mínimo ciento de eficiencia al aumentar el pH es mucho

58 TECNOLOGÍA QUÍMICA Vol. XXXIV, No. 1, enero-abril, 2014

 mayor cuando se opera a presiones elevadas eficiencia de precipitación del Níquel, siendo el nivel
(P = 977, 10 kPa), lo que se refleja claramente de operación máximo el más recomendable.
por la pendiente de la línea superior. Esto
indica que en el período analizado pueden existir La figura 2 refleja que un aumento de la
interacciones entre estas variables. concentración de Manganeso tiene un efecto
favorable sobre la eficiencia de precipitación del
De acuerdo con los resultados observados, como Níquel. Puede apreciarse que en ambos niveles
muestra la figura 1, es evidente que el aumento de de presión al aumentar la (CMn(Le)) aumenta el por
presión y de pH son condiciones favorables para la ciento de eficiencia.

Fig. 2 Efecto de los cambios de la CMNLE sobre la eficiencia de precipitación del Ni.

Por otra parte se observa claramente que los Lo anterior permite afirmar que existe interacción
resultados dependen del nivel de presión en que entre estas variables y que no es recomendable
estén trabajando los reactores, ya que el efecto operar los reactores en el nivel de operación con
presiones máximas cuando la concentración de
que produce el aumento de (C Mn(Le)) sobre la
Manganeso es elevada, porque el resultado será la
eficiencia es menor en el nivel de presión mínima, disminución en la eficiencia de precipitación del
comparativamente con el de presión máxima, níquel, lo que afectaría el rendimiento y eficiencia de
nivel en el cual el efecto de la variable es mucho la operación con las consecuentes pérdidas
más significativo. económicas que se pueden derivar de ello.

Fig. 3 Efecto de los cambios de la temperatura sobre la eficiencia de precipitación del Ni.

TECNOLOGÍA QUÍMICA Vol. XXXIV, No. 1, enero-abril, 2014

59
 La figura 3 muestra que los cambios de del cobalto (ECO), como respuesta del sistema,
temperatura no ejercen efecto alguno sobre la ante los cambios en las variables: pH, concentración
eficiencia de precipitación del níquel, de Manganeso y temperatura, respectivamente.
independientemente del nivel de presión en el que
operen los reactores de precipitación, lo que En la figura 4 se puede observar que la
demuestra que entre estas variables no existen eficiencia de precipitación del cobalto aumenta
interacciones. El resultado confirma que los cuando se incrementa la presión desde 512,3 a
cambios en esta variable puede tener un efecto 977,10 kPa, manteniendo el pH constante en el
resultante no significativo para el sistema. mínimo nivel e igualmente aumenta cuando se
mantiene el pH constante en el máximo nivel.
Sin embargo, a diferencia de lo observado en el
Análisis de la influencia de pH, concentración caso del níquel, el cambio que experimenta la
de manganeso y temperatura sobre la respuesta ECO al aumentar el pH es
eficiencia de precipitación del cobalto aproximadamente de igual magnitud,
independientemente del nivel de presión con
En las figuras 4, 5 y 6 se puede apreciar el que se trabaje, lo que indica que no existen
comportamiento de la eficiencia de precipitación interacciones entre estas variables.

Fig. 4 Efecto del pH sobre la eficiencia de precipitación del Co.

Fig. 5 Efecto de los cambios de la concentración de manganeso (CMn(Le)) sobre la eficiencia de

precipitación de Co.

60 TECNOLOGÍA QUÍMICA Vol. XXXIV, No. 1, enero-abril, 2014

 En la figura 5 se muestra que un aumento de descrito anteriormente es mucho más
la concentración de manganeso ejerce un efecto significativo.
positivo sobre la eficiencia de precipitación del La figura 6 muestra que los cambios de
cobalto, se puede apreciar que en ambos niveles temperatura no ejercen efecto alguno sobre la
de presión al aumentar la (C Mn(Le)) aumenta el eficiencia de precipitación del cobalto,
por ciento de eficiencia, además se aprecia independientemente del nivel de presión en el que
claramente que los resultados dependen del operen los reactores de precipitación, lo que
nivel de presión en que se esté trabajando, ya demuestra que entre estas variables no existen
que el efecto que produce el aumento de (C Mn(Le)) interacciones, de manera similar a lo que ocurre en
sobre la respuesta ECO es menor en el nivel de el caso del níquel, como se explicó anteriormente, lo
presión mínima, comparativamente con el de que permite afirmar que la temperatura no es una
presión máxima, nivel en el cual el efecto variable significativa en ninguno de los dos casos.

Fig. 6 Efecto de los cambios de la temperatura sobre la eficiencia de precipitación del Co.

De manera general, el comportamiento ENI: Eficiencia de precipitación del níquel (%)

observado será válido siempre que las variaciones ECO: Eficiencia de precipitación del cobalto (%)
de cada variable se introduzcan dentro del intervalo
establecido por las normas para dicho parámetro Bibliografía
en el período analizado, para el cual fue validado
el modelo matemático. 1 HERNÁNDEZ P. C.: "Modelo de simulación para el
comportamiento de los autoclaves de precipitación
de sulfuros de níquel y cobalto". Tesis en opción
Conclusiones al Grado Científico de Doctor en Ciencias Técnicas.,
Universidad de Oriente, 2002.
1 El estudio de la interacción entre las variables 2 PIÑO C. Y.: "Análisis de la influencia de las variables
del proceso seleccionadas para este análisis, de operación sobre el comportamiento de los
reactores del proceso de precipitación de sulfuros
demuestra que las variables que ejercen mayor de Ni y Co". Trabajo de diploma, Facultad de
influencia sobre la eficiencia de precipitación Ingeniería Química., Universidad de Oriente, 2012.
del níquel y cobalto son el pH y la concentración 3 RIVERA S. M. y col. "Análisis de sensibilidad paramétrica
del manganeso en el licor de entrada. en reactores de reformación catalítica de LPG". vol. XXVI
N0: 2. Revista: Tecnología Química, 2005.
4 RIVERA S. M.; HERNÁNDEZ P. C.; MATOS D. R.;
Nomenclatura VIERA B. A. R: "Sensibilidad paramétrica relacionada
con la operación de reactores de reformación catalítica
A, B, C: Compartimientos del reactor continuo de LPG". Revista Tecnología Química, Vol.29, Nº 1, 2009.
con agitación utilizado en el proceso 5 Varma, A.; Morbidelli, M. and Wu, H.: "Parametric
C Mn(Le): Concentración de Manganeso en el Sensitivity in Chemical Systems". Cambridge
licor de entrada University Press, New York, 1999.

TECNOLOGÍA QUÍMICA Vol. XXXIV, No. 1, enero-abril, 2014

61