Laboratorio de investigacin ICQ-392 Primer semestre 2013

Deteccin de fallas operacionales y de medicin en circuito de tres

celdas de flotacin

Profesor tutor: Luis Bergh O
a
, Alumno: Claudio Ramrez A
b
, Ayudante: Isabel Panire N
b

a
Ph.D, Master Engineer, Ing. civil qumico, Universidad tcnica Federico Santa mara, profesor ramo ICQ-392, tutor del alumno autor.
b
Estudiantes de ingeniera civil qumica, Universidad tcnica Federico Santa Mara, Casa Central





INFORMACIN DEL ARTCULO

Historial de artculo:
Fecha entrega, lunes 8 de junio
Fecha presentacin, martes 9 de junio

Palabras clave:
Control PI
Celdas mecnicas
Flotacin de minerales
Simulador metalrgico
Mtodo PCA






R E S U M E N

Un control de calidad aplicado a un proceso complejo y sensible como la
flotacin de minerales depende de la precisin de las mediciones de las
variables de proceso y una correcta toma de decisiones tanto en la
modelacin, parmetros de control automtico y operacin. Lo anterior puede
en algn momento causar posibles fallas operacionales y de instrumentacin,
modificando la recuperacin de mineral y ley esperada. Se realiz un estudio
en una planta piloto conformada por tres celdas mecnicas y mediante el uso
de tecnologa computacional y fenomenologa sobre flotacin se logr
implementar un simulador metalrgico con registro de datos on-line. Lo
anterior en conjunto con la aplicacin de mtodos estadsticos multivariados,
en particular PCA result ser una til herramienta para reducir la
dimensionalidad y dificultad del problema lograndose identificar las causas
de variabilidad en el proceso.



1. Introduccin

La flotacin de minerales, como proceso fisicoqumico de
tres fases (slido, lquido y gas) tiene por objetivo la
separacin selectiva mediante la adhesin de partculas de
mineral valioso a burbujas de aire, siendo esto posible
gracias a las caractersticas hidrofbicas de las partculas de
inters siendo este el principio bsico en el que se
fundamenta el proceso de flotacin. Actualmente en la
industria existen variados tipos de procesos en el campo de
concentracin de minerales, destacando las columnas de
flotacin, celdas mecnicas, celdas neumticas, flotacin
flash, flotacin centrifuga y celdas J ameson.
(Jameson,1985).

En la flotacin de minerales de cobre se utilizan
mayoritariamente las etapas rougher, cleaner, scavenger,
cleaner-scavenger y recleaner cuyos objetivos involucran
una alta recuperacin de las especies tiles con la mayor
selectividad posible. Para cumplir con estos objetivos los


circuitos estn divididos en etapas, ordenndose en bancos
de celdas y columnas de flotacin.

La etapa primaria de flotacin (etapa rougher) se alimenta
con el rebalse de los hidrociclones de un circuito cerrado
molienda/clasificacin. La etapa Cleaner limpia la corriente
de concentrado Rougher cuyo objetivo es aumentar la ley,
desde esta etapa se obtiene el concentrado final y la
corriente de cola se recircula a la alimentacin de la etapa
Rougher. La etapa Scavenger se encarga de la corriente de
cola Rougher cuyo objetivo es aumentar la recuperacin,
desde ac se genera la corriente de cola final y la de
concentrado se recircula a la alimentacin de la etapa
Rougher En relacin a las celdas de flotacin utilizadas en
los circuitos, las celdas mecnicas son utilizadas en las
etapas rougher, scavenger y cleaner-scavenger, mientras
que, columnas de flotacin se aplican a las etapas cleaner y
recleaner.

Figura 1: Ejemplo de circuitos de flotacin
Respecto al desarrollo actual en el rea de control de
flotacin las mejoras corresponden a acrecentamientos en
diferentes reas: mediciones, modelos funcionales y
dinmicos, estrategias y algoritmos de control. Desde la
dcada de 70, uno de los principales avances en control de
flotacin estuvo relacionado con el uso de mediciones en
lnea de leyes de mineral, de ah la oportunidad para
automatizar el control de medicin directa de un objetivo
metalrgico.

Antiguamente era comn ver en los circuitos de flotacin
con numerosas cmaras filmadoras de video que permitan
a un operador de forma remota vigilar cada celda y
observar el rebalse de la espuma, tomando decisiones
respecto a lo que observaba pero con el pasar de los aos y
el desarrollo de la tecnologa se ha logrado una
caracterizacin en lnea de la superficie de espumas
industriales, donde la medicin de recuperacin en lnea
resulta imperativa para optimizar la operacin.
Actualmente se ha logrado un control efectivo en celdas
mecnicas de flotacin usando un sistema de control a lazo
cerrado que coordina la manipulacin del nivel de pulpa, el
flujo de aire, el flujo de la salida de cola y el flujo inicial de
pulpa que entra al banco de celdas, con el fin de mantener
la recuperacin en un rango determinado.

2. Descripcin de la planta piloto.

2.1 Sobre el circuito de celdas.
El circuito de flotacin piloto est ubicado en el tercer piso
del laboratorio de operaciones unitarias de la universidad
tcnica Federico Santa Mara. Est conformado por tres
celdas de acrlico de forma paraleleppeda, cada una con 31
[cm] de ancho, 31 [cm] de largo y 40 [cm] de alto, con cada
placa de acrlico con 1 [cm] de grosor, cada celda con 36
[L] de capacidad. Estas fueron construida de manera de
canalizar el rebalse de fluido al sobrepasar los 40 [cm],
existiendo una canaleta que descarga de forma directa en el
estaque de alimentacin. Las celdas se encuentran ubicadas
en desnivel de forma descendente de la primera a la tercera.
Estas son del tipo mecnicas, debido a la presencia de
agitador y el ingreso de aire corresponde al tipo forzado. El
fluido utilizado en la planta piloto corresponde a agua
potable mezclado con espumante MIBC (Metil-isobutil
carbinol en concentracin de 15 [ppm]. El agua es
bombeada mediante una bomba peristltica desde un
estanque ingresando a la primera celda, en donde el flujo de
cola de la celda uno corresponde a la alimentacin de la
celda dos y el flujo de cola de la celda dos corresponde a la
alimentacin de la celda tres. Finalmente el flujo de cola de
la celda tres retorna al contenedor inicial siendo descargada
a presin atmosfrica al estanque de alimentacin
formando un circuito cerrado. Existen tres vlvulas de bola
manuales de forma paralela a las vlvulas de control de
colas, su uso se sita cuando ocurre saturacin de los
recursos a grandes flujos entrada.

2.2 Instrumentacin en la planta piloto.
La instrumentacin en el circuito de celdas est
conformada por 6 sensores y 7 actuadores, los que se
detallan a continuacin:

Tabla 1: Sensores instalados en la planta piloto
Instrumento Sensor Medicin
Sensor ultrasnico LI 1 Nivel celda 1
Sensor ultrasnico LI 2 Nivel celda 2
Sensor ultrasnico LI 3 Nivel celda 2
Flujometro msico FI 1 Flujo aire celda 1
Flujometro msico FI 2 Flujo aire celda 2
Flujometro msico FI 3 Flujo aire celda 3

Respecto a la tabla 1 se puede mencionar que cada celda de
flotacin utiliza un sensor de nivel ultrasnico que permite
registrar el nivel de espuma. El flujometro msico de aire
registra el flujo que ingresa a cada una de las celdas.



















Respecto a la tabla 2, los actuadores, las vlvulas de
control de las colas corresponden a vlvulas pinch y las de
entrada de aire a vlvulas tipo globo.Cada una de estas
vlvulas tiene asociado un transductor I/P que es el
encargado de traducir una seal elctrica proveniente del
Programador Lgico Controlable (PLC) de 4-20 [mA] a
una seal neumtica de 3-15 [Psi], lo que permite modificar
el porcentaje de apertura de las vlvulas actuadoras.
Respecto al resto de la planta se puede mencionar un
regulador de presin manual que modifica el ingreso de
aire, adems de los agitadores mecnicos asociados a cada
celda, en el cual solo se posee la capacidad de prenderlos y
apagarlos a de forma remota debido a que las [RPM] de
estos se modifica de manera manual.

Tabla 2: Actuadores asociados a la planta piloto
Instrumento Controlador Recurso
Bomba peristltica B-1 F. alimentacin
V. de control aire FIC-1 F. aire celda 1
V. de control aire FIC-2 F. aire celda 2
V. de control aire FIC-3 F. aire celda 3
V. de control colas LIC-1 F. cola agua C1.
V. de control colas LIC-2 F. cola agua C2.
V. de control colas LIC-3 F. cola agua C3.

3. Realizacin de la interfaz HMI.

3.1 PLC.
Un PLC es un dispositivo electrnico muy usado en
automatizacin industrial cuya funcin principal es
controlar la lgica de funcionamiento y realizar estrategias
de control, tales como controladores PID (Proporcional
Integral y Derivativo) de un proceso.

Los PLC son diseados para mltiples seales de entrada y
de salida, rangos de temperatura ampliados, inmunidad al
ruido elctrico y resistencia a la vibracin y al impacto. En
la planta piloto el PLC asignado corresponde al n2, siendo
este un modelo de las series 90-30 de Ge Fanuc.

3.2 Identificacin de conexiones elctricas.

El primer paso para la realizacin de una interfaz humano
mquina (HMI) corresponde a corroborar las conexiones
fsicas (cableado) desde cada uno de los sensores y
actuadores hasta el PLC. Se realiz un seguimiento fsico
de los cables utilizando un tester en modo de continuidad,
verificando a que canal y modulo del PLC corresponda
cada salida de la bornera, se identificaron y se registraron
por numeracin y color correspondiente para
posteriormente conocer el TAG que corresponde a una
mtodo de combinacin alfanumrica para asignar y
referenciar ubicacin de memoria en el PLC.

3.3 Comunicacin Software PLC

La comunicacin entre los sensores, PLC, ordenador y
usuario es escencial en todo control de procesos, es por este
motivo que se utilizan cuatro principales programas que
permiten el intercambio de informacin, estos son
VersaPro, Matlab ,Intouch y Microsoft excel. El primero de
estos programas permite, con una programacin previa de
los racks y canales del PLC, mostrar los valores que estn
midiendo los sensores en bits. En el programa Intouch
Wonderware, se pueden desarrollar HMIs que permiten
interactuar de manera fcil al usuario con el proceso en
unidades ingenieriles realizando un escalamiento de estas.
Matlab, por otra parte, permite la programacin tanto del
Figura 2: P&ID planta piloto, circuito celdas de flotacin


0
5
10
15
20
0 10000 20000 30000 40000
P
r
e
s
i

n

[
P
s
i
]

Bits
Transductores I/P vlvulas de aire
V. Aire celda 1
V. Aire celda 2
V. Aire celda 3
simulador metalurgico como del algoritmo por varianzas
utilizado para deteccin de estado estacionario (Rhinehart,
1995). Microsofts Excel permite un registro de datos en
tiempo real mediante una conexin directa al Software
intouch previa programacin en visual basic.
4. Calibracin de la instrumentacin.

La evaluacin y correcto funcionamiento de las
condiciones generales de los equipos e instrumentacin
depende de una correcta calibracin e inspeccin de los
instrumentos involucrados. Es comn que la
instrumentacin industriales pierda efectividad en el
tiempo, por lo tanto, no aseguran la calidad y certeza de las
mediciones y actuaciones que son claves en el proceso. Por
lo tanto, previo a la implementacin del sistema de control
se evalu las condiciones de los equipos y se caracteriz y
calibr la instrumentacin de la siguiente manera:

i. Para la calibracin de los sensores de nivel
ultrasnicos se midi la distancia utilizando
huincha mtrica ajustando el valor mnimo en 4
[mA] y el mximo de la celda en 20 [mA].

ii. Para los medidores de flujo se caracteriz con
ayuda de un rotmetro.

iii. Para la caracterizacin de la bomba se envi
seales en dentro del rango de bits (0 a 32000) y
se midi el caudal mediante aforo, posteriormente
se ajust mediante el anlisis numrico de
mnimos cuadrados ponderados.

iv. Para verificar el correcto funcionamiento de las
bombas mediante un multitester se enviaron
seales de 0 y 32000 [bits] donde se verific si la
apertura o bloqueo total de estas coincida con 3 y
15 [psi] respectivamente.


Figura 3: Calibracin Vlvulas de aire.



Figura 4: Calibracin Bomba peristltica de alimetancin
5. Parametros de control.

Mediante el programa VersaPro se realiz la configuracin de
la lgica de funcionamiento para ser cargada en el PLC. La
parte ms importante corresponde a 6 bloques PID-IND, donde
se ingresaron parmetros proporcional e integral descartndose
el parmetro derivativo debido a que el uso de este es ms
adecuada cuando hay retraso entre el movimiento de la vlvula
de control y su repercusin a la variable controlada, no siendo
este el caso. En la siguiente tabla hace referencia a los
parmetros, obtenidos en un comienzo, mediante registro de
datos en lazo abierto con matlab y posterior uso de PID tool,
obtenindose parmetros de partida, luego segn ensayo y error
se definieron los siguientes parmetros para las 6 vlvulas de
control.

Tabla 3: Parmetros ingresados a bloques de control.
Flujo de aire Nivel de agua
Celda 1 P: 0 I: 0.33 P: 1.5 I: 0.5
Celda 2 P: 0 I: 0.33 P: 5.0 I: 0.010
Celda 3 P: 0 I: 0.33 P: 1.75 I: 0.010

6. Modelacin y simulacin

6.1 Caracterizacin individual de cada celda.

Para efectos de modelado, anlisis y correcta
implementacin del simulador es necesario identificar dos
zonas importantes en cada una de las celdas, la zona de
coleccin y la zona de limpieza (espuma). En la primera
mencionada es donde se forma el agregado partcula-
burbuja y llevado a la interface pulpa-espuma. En la zona
de limpieza localizada entre la interfaz de pulpa-espuma y
el desbordamiento de concentrado, donde existe la
0 0.5 1 1.5 2 2.5 3 3.5
x 10
4
0
1
2
3
4
5
6
7
8
9
Bits
C
a
u
d
a
l

[
]


Datos
Reg ponderada
Reg no ponderada


posibilidad de que partculas arrastradas vuelvan a la zona
de coleccin, por lo tanto esta zona se describe segn la
recuperacin de espuma.

(1

)
(1)

La ecuacin n1 describe la recuperacin global por celda
(Rg), relacionando esta segn la recuperacin en la zona de
coleccin (Rc) y la zona de limpieza (Rf), esta relacin
proviene del balance general de masa en torno a la celda.
(Finch y Dobby, 1990;. Finch et al, 2008).

6.2 Zona de coleccin.
La recuperacin de mineral en esta zona (Rc) puede ser
descrita en general por la ec. n2. (Polat y Chander,2000;.
Yianatos, 2007)

(1

)()() (2)

Para la recuperacin en la zona de espuma R
MAX
representa
la mxima recuperacin por flotacin en un tiempo infinito,
E(t) corresponde a la funcin de distribucin de tiempo de
residencia para un proceso continuo, por otro lado F(k) es
la funcin de distribucin de constantes cinticas para cada
especie. El modelo considera la recuperacin de mineral
como proceso con cintica de primer orden.

Para la distribucin de tiempos de residencias, se considera
la expresin mostrada en la ecuacin n3, que considera
una distribucin de tiempos de residencia de primer orden
(Yianatos, 2005).

() =
1

(3)

Para la funcin de distribucin de constante cinticas para
cada especie se consider el modelo rectangular de Klimpel
(1980)

() =
1

, 0 < <
,

0 , >
,

(4)

Si se asocian las ecuaciones presentadas e integrando la
ecuacin n2 se obtiene una expresin final simplificada
que simula a la recuperacin en la zona de coleccin.
(Yianatos, 2006)

1
ln(1

(5)

6.3 Zona de limpieza.
La zona de limpieza fue caracterizada segn las principales
variables que afectan la recuperacin de mineral, siendo
estas la altura de espuma, el flujo de aire (o velocidad
superficial del gas, J
g
[cm/s]), Rmax y tiempo de
residencia, segn la siguiente ecuacin usualmente
empleada (Gorain et al.,1998; Mathe et al., 1998; Zheng et
al., 2004):

,
=

(6)

De la ec. n6 se comenta que corresponde a la
recuperacin mxima en la zona de coleccin,
corresponde al factor de estabilidad de la espuma y
f

corresponde al tiempo de residencia del gas en la espuma
de la celda ensima. Debido a que las tres celdas poseen las
mismas caractersticas los parmetros son equivalentes para
cada una de ellas. Se observa que a menor tiempo de
residencia del gasse logra maximizar la recuperacin de
mineral. El trmino
f
viene dado segn la ecuacin n9
(Yianatos et al., 2008):

,
=

(7)

Donde H
f,n
es la altura de espuma y
g
corresponde al hold-
up.

6.4 Recuperacin de agua.
Para fines de simulacin, el modelo utilizado para la
recuperacin de agua, definiendo esta como la cantidad de
agua que se obtiene en el concentrado en la celda ensima
sobre la cantidad de agua que ingresa a esta, fue
determinado de forma emprica (Yianatos y Bergh ,2012)
en funcin de la altura de espuma, el factor de estabilidad
de la espuma y velocidad superficial del gas para cada
celda.

,
=

(8)

En la ec. N8, corresponde a la mxima recuperacin de
agua, adems y son parmetros de ajuste..

6.4 Recuperacin de Ganga.
Durante el proceso de flotacin existe la posibilidad de que
ganga (mineral sin valor) sea arrastrada por el agua desde
la zona de coleccin a la zona de limpieza por lo tanto la
recuperacin de ganga est directamente relacionada con la
recuperacin de agua en la celda ensima, mediante la
siguiente expresin (Yianatos y Contreras, 2010):

,
(9)



Donde EF
i
corresponde al factor de arrastre por clase de
tamao, obtenido de manera emprica.

6.5 Obtencin de la ley de cobre en el concentrado.
Por definicin la ley se define como una medida que
describe el grado de concentracin de recursos naturales
valiosos (metales o minerales), disponibles en la mena, es
decir como la fraccin de material til contenido en el
concentrado. La expresin para la celda ensima se
presenta segn (Yianatos, 2005):

,
=


(10)

7. Deteccin de estado estacionario

Debido a que el simulador metalrgico funciona slo en
estado estacionario, surge la necesidad de un test en lnea
de deteccin del estado del sistema. (Rhinehart, 1995)
propuso un mtodo que consiste en calcular la varianza de
una variable de proceso mediante dos mtodos diferentes.
Idealmente si el proceso se encuentra en estado
estacionario ambas varianzas calculadas por mtodos
diferentes debieran arrojar el mismo valor y por la tanto la
razn entre ellas debiera ser igual a 1. En la razn entre las
varianzas, el numerador es estimado a partir de
desviaciones entre informacin individual y valores
filtrados, mientras que el denominador es estimado de
desviaciones entre informacin sucesiva

2

1 (11)

Segn Rhinehart, la relacin de varianzas puede ser
calculada de la siguiente manera:

= (2

2
(12)
Segn literatura, un valor aceptable para

es de 1.6, por lo
tanto mediante mediciones del sistema se ajust los valores
de

necesarios para que esta condicin se cumpla. Se

realiz una medicin de datos cada 2 segundos, ajustando
los valores de

para las variables ms difciles de

controlar (niveles de agua). Se realiz una suavizacin de

2
y adems disminuir el valor de
3
recomendado a la
mitad, siguiendo una tendencia parecida a
2
lo que se
traduce en un menor valor de

, los resultados se pueden

apreciar en la figura 5. Se obtuvieron finalmente los
siguientes valores para el mtodo de identificacin de E.E.

Tabla 4: Valores de para el mtodo de identificacin de E.E

1

2

3

0.01 0.02 0.05


Figura 5: Estadsticos auxiliares para deteccin de estado
estacionario

Figura 6: comportamiento r
3
con = 0.05
8. Mtodo de anlisis de componentes
principales:

El mtodo PCA es una tcnica utilizada para reducir la
dimensionalidad de un conjunto de datos siendo til la
tcnica para hallar las causas de la variabilidad de un
conjunto de datos. La aplicacin de mtodos estadsticos
multivariados, y en particular PCA, es una herramienta
para construir modelos lineales que contienen los
elementos esenciales del proceso de fenmenos con el
nmero mnimo de variables latentes. Si se posee un
problema de N variables, se traduce en un problema de N
dimensiones, aumentando la complejidad de este, adems
al estar las variables correlacionadas agregan un problema
de estimacin de parmetros para el modelo. Lo anterior se
puede aprovechar como una oportunidad para reducir la
dimensin del problema, es decir proyectarlo en una
dimensin menor (y disminuir la complejidad). Esto se
logra creando un conjunto menor de variables auxiliares
lineales e independientes entre si creando una condicin
optima de estimacin de parmetros para un modelo
0 500 1000 1500 2000 2500 3000
0
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
0.18
Dato


sigma
2
rho
2
0 500 1000 1500 2000 2500 3000
0
0.5
1
1.5
2
2.5
3
3.5
4
Dato
r


8.1 PCA en deteccin de fallas de operacin y medicin

En una aplicacin tpica de anlisis de datos industriales
con k cantidad de variables de proceso son medidos
simultneamente y son combinados con un vector asociado
con cada observacin de las k variables. Si se realizan N
observaciones, se obtendr una matriz de tamao X=(N*k)
donde cada fila corresponde a una observacin para las
variables k y cada columna corresponde a la variable
especifica en la observacin N. La descomposicin de la
matriz X viene dada segn:

=

+ (13)

Para esta descomposicin las variables latentes son las
columnas de la matriz T, es decir en cada una de las N
observaciones de la X se calcula una fila de la matriz T,
siendo llamadas estas comnmente como Scores y
calculadas por medio de mapeo matemtico d las variables
originales medidas a travs de la matriz de carga P. El
factor presentado corresponde al modelo residual como
resultado de esta asignacin.

El modelo PCA puede ser visto como una descomposicin
valor propio-vector propio de la matriz de informacin XX
de los datos medidos. La dimensin del modelo (m) se
encuentra mediante la aplicacin de criterios de seleccin
significativos que conservan slo los grandes valores
propios de estas matrices y rechazan el resto de los
modelos. La dimensin de las variables latentes resultantes
se compone por un subespacio de los datos originales
medidos. En el mtodo PCA la matriz de carga P est
conformada por los vectores propios que son retenidos en
el modelo resultante, y la matriz T almacena la proyeccin
de los datos originales X definido por estos vectores
propios. (Bergh et al., 2008)

8.1 herramientas de inferencia estadstica:
La ventaja principal de este mtodo (PCA) es que cualquier
comportamiento que pueda ser considerado como anormal
se puede detectar mediante combinaciones de los mtodos
SPE (Square Predicted error) y distribucin de T
2
de
Hotteling.

2
=

=1
(14)

La ec. n14 presenta el mtodo de clculo de T de
Hotteling para cada observacin registrada.

2
es la
varianza de los datos
2

de partida. El lmite de control

superior de este test se obtiene utilizando la ec. n15.

=
(
2
1)
( )

,(),
(15)
De la ecuacin anterior se define a m como el nmero de
componentes, F corresponde al a distribucin F y al nivel
de significancia para la hiptesis multivariada.

Respecto a la prueba SPE, este es calculado como la suma
cuadrada a travs de los resultados de cada una de las de las
observaciones (i) indicada para cada variable (j)

=
2

=1
(14)

El lmite de control superior para este test con una
significancia est dado por.

2

2

,
(15)

Con
2
, distribucin chi cuadrado, es la varianza de la
muestra y m corresponde al promedio de cada SPE.

8.2 Criterios de identificacin de fallas de medicin y
operacin:
Los resultados de las pruebas basadas en los registros on-
line otorgados meditando HMI y simulador metalrgico se
basaron en los siguientes criterios (Bergh et al., 2008):

(i) Se considera operacin normal si el nuevo set de datos
satisface los test T
2
y SPE.
(i) (ii) Se considera operacin anormal si falla T
2
y se cumple
SPE.
(ii) (iii) Problemas de medicin o de representacin de modelo
PCA si falla SPE y se cumple T
2

(iii) (iv) Si ambos test fallan el modelo PCA no es apropiado o
se ha producido problemas de medicin.
(iv)
9. Implementacin del anlisis PCA:

El procedimiento general para la generacin y obtencin de
resultados incluye los siguientes pasos:

i) Se llev el circuito de celdas mediante la interfaz
HMI y los lazos de control activados a un estado
estacionario, siendo los niveles medidos por el
sensor, verificando el estado estacionario tanto por
el usuario como por la alarma de deteccin por
varianzas.
ii) Se gener una matriz de datos considerando todas
las variables tanto reales (medibles de la planta
piloto) como virtuales (generadas mediante el
simulador metalrgico). El nmero de variables
consiste en 73 para las cuales se obtuvieron
21.827 observaciones siendo el metro de registro



el software Microsofts Excel mediante la
programacin de una Macro en lenguaje visual
basic.

iii) Posteriormente los datos fueron importados desde
Excel al software ProSensus Multivariate 13.02
versin 2013. Mediante este se realiz el anlisis
PCA para la obtencin de resultados.

9.1 Caractersticas estado estacionario:
El E.E utilizado para la medicin consisti en un set point
de 36, 34 y 32 [cm] de nivel para la celda 1,2 y 3
respectivamente. Esta disminucin para emular los niveles
de espuma reales en la industria. Se utiliz flujo de aire de
45 [l/min] lo que se traduce en un J
g
de 0.99 [cm/s] en cada
celda y flujo de bomba de 1.5 [l/min]. En un comienzo se
pens en utilizar flujos mayores, entre 6 y 8 [l/min], dado la
existencia de vlvulas de bola utilizadas como by-pass para
evitar la saturacin de los recursos sin embargo solo se
alcanzaron estados meta-estacionarios con una variacin
muy grande en las aperturas de las vlvulas de nivel sub
amortiguando el comportamiento de cada una de estas.

10. Resultados:

Los siguiente resultados fueron realizados en base a
recopilacin de datos en un estado estacionario .


Figura 7: Modelo PCA con 6 componentes principales.
Respecto a la figura 7, se observa que el modelo es capaz
de explicar cerca del 97% de la variacin de los datos, no
existiendo una diferencia significativa en la eficacia del
modelo en relacin a la presencia de 5 o 6 componentes
principales. A continuacin en las figuras 8 y 9 se
presentan ambos test por separado.

Figura 8: Resultado Test T de Hotteling.

Figura 9: Resultado Test SPE.
Los lmites superiores corresponden a 13 para el test de T
de Hotteling y de 3.25 para el test SPE, ambos con un 95%
de confianza, los lmites inferiores en ambos casos es 0.
10.1 Comparacin de test para anlisis de operacin
normal:
El siguiente grfico representa un anlisis grfico de SPE
v/s T de hotteling. Este resulta muy til para utilizar los
criterios de inferencia estadstica propuestos en el punto
8.2. Se observa que casi la totalidad de los datos se
concentra bajo el 95% de ambos test, lo que se considera
como operacin normal, sin embargo se debe destacar la
presencia de datos dispersos que segn los criterios
propuestos se pueden clasificar como fallas de operacin o
fallas de medicin.



Figura 10: SPE v/s HT2.
10.1.1 Identificacin de falla de operacin en anlisis de
datos.
Se considera falla de operacin si en alguna medicin falla
T
2
y se cumple SPE. Respecto a la figura n9 notamos que
el caso ms claro que cumple este criterio corresponde a la
medicin nmero 102. Esto se pude corroborar analizando
especficamente tal medicin. En donde obtiene lo
siguiente:

Figura 11: Contribucin P en observacin nmero 102.
Lamentablemente la gran cantidad de variables asociadas al
proceso impide una correcta lectura de la figura, sin
embargo se opt por dejarla completa para hacer una
comparacin de la contribucin P en la medicin nmero
102. Las tres barras azules corresponden a a los flujos
msicos de aire para las celdas 1, 2 y 3, por lo tanto se
ratifica una falla operacional ocurrida en el instante de
medicin 102, corroborando el criterio. (La imagen
ampliada se encuentra en el anexo A)

10.1.2 Identificacin y anlisis en falla de medicin.
Se considera como falla de medicin (instrumentacin) si
falla SPE y se cumple T
2
. Volviendo a la figura n9 se
observa que el punto 297 est bajo este criterio. Analizando
ese punto en particular mediante la contribucin de P en tal
observacin:


Figura 12: Contribucin P en observacin nmero 297.
En la figura 10, se busc entregar una visin general de la
contribucin P de todas las variables del sistema, debido a
la gran cantidad de variables se opt por realizar una figura
explicativa y adjuntar la imagen completa en el anexo A.
El primer ovalo de la izquierda (1), muestra la
contribucin de las variables nivel_2 y nivel_3, que
corresponden a la altura de nivel de las celdas 2 y 3. El
ovalo (2) muestra la contribucin de las variables V_cola2
y Vcola_3, siendo estas las vlvulas de control asociadas a
las mismas celdas. El ovalo nmero 3 muestra la
repercusin que tuvo lo ocurrido en (1) y (2) en las
variables virtuales del simulador que estn en funcin de
los valores reales medidos en la planta piloto. Nuevamente
un uso de criterio para anlisis PCA logra justificar en este
caso un punto con falla de medicin.
11. Conclusiones:
La calidad del control de un banco de celdas de flotacin
depende en gran medida de la precisin de las mediciones y
estimaciones realizadas. Caracterizando la flotacin como
un proceso complejo es difcil decidir cunto y cuales
variables de deben modificar con el fin de impulsar el
proceso a una operacin normal y beneficiosa. La
aplicacin de mtodos estadsticos multivariales (PCA) es
una herramienta poderosa parar construir modelos lineales
independientes entre s con el mnimo nmero de variables
latentes. Respecto a la planta piloto se logr cumplir con
los objetivos planteados, desde la implementacin de
control automtico, comunicacin PC-PLC y utilizacin de
software para la realizacin de HMI. La implementacin
del modelo PCA a un circuito de celdas de flotacin fue
realizada con xito, utilizando 6 componentes principales,


corroborando los criterios que permiten detectar fallas de
medicin y operacin, obtenindose coherencia en los
resultados obtenidos. La seleccin de datos fue bajo una
operacin normal de planta obtenindose resultados
congruentes y de buena calidad, esto debido al buen control
que se logr de la planta adems de buena dinmica y
calidad de la instrumentacin.

11. Referencias:

R. Russel Rhinehart (1995) a novel method for automated
identification of steady state

Luis G. Bergh, J . Yianatos, S. Acosta (2008) Online
monitoring of a pilot column flotation columns by using
PCA models

Luis G. Bergh Monitoreo y diagnstico de la operacin
de procesos de flotacin usando tcnicas de proyeccin

L,Bergh, J , Yianatos, L.Vinnett (2012) Modelling and
simulation of Rougher flotation circuits

L. Bergh, J . Yianatos (2013) Control of rougher flotation
circuits aided by industrial simulator

J .Yianatos, (2005) Flotacin de minerales





































































ANEXO A: Ampliacin de imgenes.



Figura 11 AMPLIADA: Contribucin P en observacin nmero 102.


Figura 12 AMPLIADA : Contribucin P en observacin nmero 297.