Teoria Simulacion Molienda

MODELO, CONTROL Y SIMULADOR DE
PLANTA DE MOLIENDA SEMIAUTGENA Y MOLIENDA SECUNDARIA

MEMORIA PARA OPTAR AL TTULO DE INGENIERO CIVIL
ELECTRICISTA

ROBERTO PATRICIO ORELLANA ARAYA

PROFESOR GUA:
GUILLERMO GONZLEZ REES

MIEMBROS DE LA COMISIN:
ALDO CASALI BACELLI
HCTOR AGUSTO ALEGRA

SANTIAGO CHILE
JULIO 2010

UNIVERSIDAD DE CHILE
FACULTAD DE CIENCIAS FSICAS Y MATEMTICAS
DEPARTAMENTO DE INGENIERA ELCTRICA
RESUMEN DE LA MEMORIA
PARA OBTAR AL TTULO DE
INGENIERO CIVIL ELECTRICISTA
POR: ROBERTO ORELLANA A.
FECHA: 08/07/2010
PROF. GUA: GUILLERMO D. GONZLEZ REES

MODELO, CONTROL Y SIMULADOR DE
PLANTA DE MOLIENDA SEMIAUTGENA Y MOLIENDA SECUNDARIA

Esta memoria describe el desarrollo de un modelo y del correspondiente simulador de una
planta de molienda, constituida por un circuito de molienda semiautgena (SAG) y un
circuito inverso de molienda con bolas. El simulador incluye un sistema de control y una
interfaz grfica similar a la de las plantas reales.
En la memoria se explica primero los marcos tericos de cada unidad de la planta de
molienda y se desarrolla la integracin de estos modelos individuales y su representacin
en tiempo discreto. En segundo lugar se describe el desarrollo e implementacin del
simulador en Matlab. Finalmente se describe cmo se ajustaron los parmetros de los
modelos utilizados en el simulador.
El modelo con sus parmetros ajustados, y por lo tanto el simulador, representan una
planta genrica cualitativamente similar al comportamiento de una planta real y no
necesariamente al de una planta particular existente.
El simulador ha sido utilizado en un curso terico para operadores de plantas reales y en
cursos regulares dictados por el Departamento de Ingeniera de Minas de la Universidad
de Chile. El modelo de la planta SAG ha sido tambin utilizado en varios cursos de
posttulo que han usado el sistema Factory-Link de adquisicin y manejo de datos en el
que se ha incluido este modelo.

A mis padres Roberto e Isabel que durante toda mi vida me apoyado
permitindome ser un profesional,
y a Cristal quien me dio el ltimo empujn e hizo posible esta Memoria

i

NDICE
1 INTRODUCCIN................................................................................................ 1
2 MODELOS DE LA PLANTA SAG-MOLIENDA SECUNDARIA........................... 3
2.1 DESCRIPCIN DE LA PLANTA QUE SE QUIERE SIMULAR.....................................3
2.2 MODELOS DE LAS OPERACIONES UNITARIAS .......................................................5
2.2.1 Modelo de molino SAG.............................................................................................6
i) Proceso de los slidos............................................................................................7
ii) Proceso del agua .................................................................................................18
iii) Ecuacin de Estado y Entrada-Salida.................................................................19
2.2.2 Modelo de Harnero..................................................................................................25
i) Proceso de los slidos..........................................................................................25
2.2.3 Modelo de Chancador de pebbles ..........................................................................28
i) Proceso de los slidos..........................................................................................29
2.2.4 Modelo de Pozo ......................................................................................................33
i) Proceso de la pulpa de mineral ............................................................................33
2.2.5 Modelo de Bomba Centrfuga .................................................................................36
2.2.6 Modelo de Hidrocicln.............................................................................................40
i) Descripcin fenomenolgica.................................................................................40
ii) Modelo utilizado en el simulador..........................................................................42
2.2.7 Modelo de Molino de Bolas.....................................................................................46
i) Proceso de los Slidos .........................................................................................47
ii) Proceso del Agua.................................................................................................48
2.3 INTEGRACIN DE LOS MODELOS: MODELO GLOBAL .........................................49
2.3.1 Ecuacin de Estado y Entrada-Salida del circuito de Molienda SAG .....................49
2.3.2 Ecuacin de Estado y Entrada-Salida del circuito de Molienda Secundaria...........56
2.3.3 Modelo Global .........................................................................................................69
2.4 MODELO DE TIEMPO DISCRETO DE LA PLANTA ..................................................71
2.4.1 Molino SAG.............................................................................................................72
2.4.2 Harnero .................................................................................................................74
2.4.3 Chancador de pebbles ............................................................................................75
2.4.4 Pozo .................................................................................................................76
2.4.5 Bomba Centrfuga...................................................................................................78

ii

2.4.6 Hidrocicln ..............................................................................................................78
2.4.7 Molino de Bolas.......................................................................................................80
2.5 CONTROL DE LA PLANTA.........................................................................................82
2.5.1 Perturbaciones (medidas y no medidas).................................................................82
i) Perturbaciones en la Molienda SAG.....................................................................82
ii) Perturbaciones en el molino de bolas..................................................................85
iii) Perturbaciones en el resto de los modelos .........................................................86
2.5.2 Controles automticos en la planta.........................................................................88
i) Lazos de Control ...................................................................................................88
ii) Control PI de hold-up del Molino SAG.................................................................89
iii) Control PID del nivel del Pozo.............................................................................89
iv) Control PID del porcentaje de slidos en la alimentacin a hidrociclones..........90
v) Control PID del +65# en el Producto Final ..........................................................90
3 DESARROLLO DEL SIMULADOR ................................................................... 91
3.1 ESPECIFICACIN DEL SIMULADOR........................................................................91
3.2 PROGRAMACIN DEL SIMULADOR........................................................................91
3.2.1 Lgica constructiva..................................................................................................91
3.2.2 Mdulo DSAG.........................................................................................................92
3.2.3 Mdulo SAG_mod...................................................................................................95
3.2.4 Programacin del Mdulo PHc_mod.......................................................................99
3.2.5 Programacin del Mdulo Mbol_mod....................................................................104
3.2.6 Interfaz Grfica del Usuario e Interaccin con los Modelos..................................105
i) Pantalla de control y supervisin de la molienda SAG.......................................105
ii) Pantalla de control y supervisin de la molienda Secundaria............................108
iii) Pantalla con el Layout de la Planta...................................................................110
3.2.2 Exportacin de Datos............................................................................................112
3.3 AJUSTES DE PARMETROS ..................................................................................115
3.3.1 Antecedentes generales .......................................................................................115
3.3.2 Parmetros del Molino SAG..................................................................................118
i) Datos para la parametrizacin............................................................................118
ii) Determinacin de parmetros
S1
,

S2
.................................................................123
iii) Determinacin de las tasas especficas de molienda k
E
y de rechazo c ......124
iv) Determinacin parmetro de descarga total de slidos o
S
................................132
v) Determinacin de parmetros k
S1
,

k
S2
para descarga de agua .........................132
3.3.3 Parmetros del Harnero........................................................................................132
3.3.4 Parmetros del Chancador de Pebbles ................................................................133
3.3.5 Parmetros de las Bombas Centrfugas ...............................................................134

iii

3.3.6 Parmetros de las bateras de Hidrociclones........................................................136
ii) Determinacin del Rechazo de Agua R
CH
..........................................................138
iii) Determinacin del parmetro k
H2
para determinar el flujo de agua en el
rebose....................................................................................................................139
iv) Determinacin del parmetro k
H1
para determinar presin de pulpa en la
alimentacin...........................................................................................................139
v) Determinacin del parmetro k
H3
para determinar el d
50
...................................140
vi) Determinacin del parmetro de Plitt o
H
para determinar la matriz de
seleccin Y
H
...........................................................................................................140
3.3.7 Parmetros de los Molinos de Bolas.....................................................................143
ii) Determinacin de la tasa de descarga o
M
.........................................................145
iii) Determinacin de la matriz B de tasas de molienda.........................................147
4 CONCLUSIONES Y RECOMENDACIONES.................................................. 156
REFERENCIAS................................................................................................... 158
ANEXOS
ANEXO A Simplificacin del Modelo de chancado de Andersen / Whiten.......................160
ANEXO B Simulador molienda SAG - Manual del Usuario ..............................................164

iv

INDICE DE FIGURAS
Figura 2.1-1 Circuitos del Simulador ............................................................................ 3
Figura 2.2-1 Molino SAG FLSmidth gearless, ring motor ABB..................................... 7
Figura 2.2-2 Diagrama de bloques del Molino SAG..................................................... 8
Figura 2.2-3 Esquema representativo de flujos acumulados retenidos por malla...... 10
Figura 2.2-4 Ejemplo de harnero y diagrama de bloques........................................... 25
Figura 2.2-5 Corte de un Chancador de pebbles y esquema funcional ..................... 28
Figura 2.2-6 Diagrama de bloques del Chancador..................................................... 29
Figura 2.2-7 Diagrama de bloques del Pozo .............................................................. 33
Figura 2.2-8 Circuito de la bomba centrfuga ............................................................. 36
Figura 2.2-9 Diagrama de bloques de la bomba ........................................................ 36
Figura 2.2-10 Batera de hidrociclones y esquema del hidrocicln .............................. 40
Figura 2.2-11 Funcionamiento del Hidrocicln ............................................................. 41
Figura 2.2-12 Circuito y diagrama de bloques del Hidrocicln ..................................... 42
Figura 2.2-13 Diagrama funcional del Molino de Bolas................................................ 46
Figura 2.2-14 Diagrama de bloques del Molino de Bolas............................................. 47
Figura 2.3-1 Diagrama de Bloques del circuito de Molienda SAG ............................. 49
Figura 2.3-2 Diagrama de Bloques del circuito de Molienda Secundaria................... 56
Figura 2.3-3 Diagrama de la Planta como un Bloque................................................. 69
Figura 2.5-1 Diagrama de Bloques de perturbaciones en Molienda SAG.................. 87
Figura 2.5-2 Diagrama de Bloques de perturbaciones en Molienda secundaria........ 87
Figura 2.5-3 Lazos de control ..................................................................................... 88
Figura 3.2-1 Diagrama de flujo general del Simulador ............................................... 94
Figura 3.2-2 Pantalla de Control y Supervisin de la molienda SAG....................... 105
Figura 3.2-3 Pantalla de Control y Supervisin de Molienda Secundaria ................ 108
Figura 3.2-4 Pantalla del Layout de la Planta........................................................... 110
Figura 3.3-1 Diagrama de flujo y balance Concentradora A2 Chuquicamata .......... 120
Figura 3.3-2 Distribucin granulomtrica - Alimentacin SAG 17, Concentradora
A2 Chuquicamata ................................................................................ 121
Figura 3.3-3 Distribucin granulomtrica - Descarga SAG 17, Concentradora A2
Chuquicamata...................................................................................... 122
Figura 3.3-4 Ajuste de parmetros Molino SAG: Porcentaje distribucin
granulomtrica en alimentacin y descarga para los tamaos
representativos .................................................................................... 128
Figura 3.3-5 Curvas de eficiencia de la bomba centrfuga ....................................... 134
Figura 3.3-6 Ajuste de parmetros - Hidrociclones: Rebose acumulado: modelo
v/s datos concentradora A2 Chuquicamata......................................... 142
Figura 3.3-7 Ajuste de parmetros - Hidrociclones: Rechazo acumulado: modelo
v/s datos concentradora A2 Chuquicamata......................................... 142

v

Figura 3.3-8 Ajuste de parmetros Molino de bolas: Datos de distribucin
granulomtrica molino de bolas 17A concentradora A2 ...................... 145
Figura 3.3-9 Ajuste de parmetros Molino de bolas: Comparacin entre holdup
obtenido por matriz de molienda v/s el obtenido con tasa de
descarga .............................................................................................. 154

NDICE DE TABLAS
Tabla 3.2 - 1 Variables guardadas por el simulador ...................................................113
Tabla 3.3 - 1 Definicin de intervalos de tamaos de partculas y el tamao
representativo........................................................................................117
Tabla 3.3 - 2 Ajuste de parmetros Molino SAG: puntos para linealizacin
alimentacin...........................................................................................125
Tabla 3.3 - 3 Ajuste de parmetros Molino SAG: puntos para linealizacin
descarga................................................................................................126
Tabla 3.3 - 4 Ajuste de parmetros molino SAG: Distribucin granulomtrica........127
Tabla 3.3 - 5 Ajuste de parmetros molino SAG: Flujos msicos............................130
Tabla 3.3 - 6 Ajuste de parmetros molino SAG: Determinacin constantes de
molienda................................................................................................131
Tabla 3.3 - 7 Ajuste de parmetros hidrociclones: Distribucin granulomtrica
concentradora A2 de Chuquicamata .....................................................138
Tabla 3.3 - 8 Ajuste de parmetros hidrociclones: Reboses y rechazos
acumulados: Modelo v/s datos concentradora A2.................................141
Tabla 3.3 - 9 Ajuste de parmetros hidrociclones: flujos msicos determinados ....143
Tabla 3.3 - 10 Ajuste de parmetros Molino de bolas: Distribucin
granulomtrica molino de bolas 17A concentradora A2 de
Chuquicamata........................................................................................144
Tabla 3.3 - 11 Ajuste de parmetros Molino de bolas: Determinacin de flujos
de slidos como porcentaje del flujo total ..............................................150
Tabla 3.3 - 12 Ajuste de parmetros Molino de bolas: Determinacin de flujos
msicos de slidos ................................................................................151
Tabla 3.3 - 13 Ajuste de parmetros Molino de bolas: Determinacin valores de
la diagonal de la matriz de molienda .....................................................151

vi

LISTA DE DEFINICIN DE SMBOLOS
A) VARIABLES Y CONSTANTES POR ORDEN ALFABTICO
- A, a
i,j
= Chancador de pebbles - Matriz de clasificacin y sus elementos
- A
P
= Pozo - Superficie del pozo
- B, b
i
= Molino de Bolas - Matriz triangular inferior de molienda
- B, b
i,j
= Chancador de pebbles - Matriz de ruptura y sus elementos
- c, c
i
= SAG - tasas de rechazo de mineral
- d
CE
= Chancador de pebbles - Abertura de la descarga (closed setting) del
chancador
- d
S
= SAG dimetro del molino
d
sp
= Hidrociclones - Dimetro del orificio en la descarga (spigot).
- d
vf
= Hidrociclones Dimetro del buscador de vrtice (vortex finder).
-
A
i
A A
F , f , f = Harnero - flujos de slidos en la alimentacin
-
B
i
B B
F , f , f = Bomba centrfuga - Flujos de slidos en la alimentacin
-
C
i
C C
F , f , f = Chancador de Pebbles Flujos de slidos en la alimentacin
- F
ua
= Harnero - Representa la fraccin de mineral alimentado de tamao menor
que la apertura del harnero.
-
C i C C
F
~
, f
~
, f
~
= Chancador de Pebbles - Flujos de slidos en la alimentacin a
clasificacin
-
P
i
P F
F , f , f = Pozo - Flujos de slidos en la alimentacin
-
H
i
H H
F , f , f = Hidrociclones - Flujos de slidos en la alimentacin
-
M
i
M M
F , f , f = Molino de Bolas - Flujos de slidos en la alimentacin
-
S
i
S S
F , f , f = SAG - Flujos de slidos de alimentacin al Molino
-
H
i
H H
F
~
, f
~
, f
~
= Hidrociclones - Fraccin del flujo de slidos que va a clasificacin
-
S
i
S S
F
~
, f
~
, f
~
= SAG - Flujos de slidos de alimentacin a la cmara de molienda

vii

-
A
S
A
i
S
A
S
F
~
, f
~
, f
~
= SAG - Flujos de slidos acumulados de alimentacin a la cmara de
molienda
-
H
i
H H
F , f , f = Hidrociclones - Flujos volumtricos de la pulpa en la alimentacin
-
P
i
P P
F , f , f = Pozo - Flujo volumtrico de pulpa a la alimentacin del pozo
- G
P
= Pozo - Agua extra que se agrega al pozo
- h
H
= Hidrociclones - Altura de los hidrociclones
- h
P
= Pozo - Altura de la pulpa en el pozo
- J, J
B
= SAG - Porcentaje de llenado del mineral y de las bolas
- k, k
i,j
= SAG - Constantes de Molienda
- k
A
, k
A
i,j
= SAG - Constantes de Molienda acumuladas
- k
E
, k
E
i,j
= SAG - Constantes especficas de Molienda
- ls= SAG - Largo del molino
- n
H
= Hidrociclones - Nmero de hidrociclones
-
A
i
A A
P , p , p = Harnero - Flujos de slidos en la descarga
-
B
i
B B
P , p , p = Bomba centrfuga - Flujos de slidos en la descarga
-
C
i
C C
P , p , p = Chancador de Pebbles - Flujos de slidos en la descarga
-
H
i
H H
P , p , p = Hidrociclones - Flujos de slidos en el rebose
-
M
i
M M
P , p , p = Molino de Bolas - Flujos de slidos en la descarga
-
P
i
P P
P , p , p = Pozo - Flujos de slidos en la descarga
-
S
i
S S
P , p , p = SAG - Descargas de slidos del Molino
-
S
i
S S
P
~
, p
~
, p
~
= SAG - Descargas de slidos de la cmara de molienda
-
A
S
A
i
S
A
S
P
~
, p
~
, p
~
= SAG - Flujos de slidos acumulados en la descarga de la cmara
de molienda
-
P
i
P P
P , p , p = Pozo - Flujo volumtrico de pulpa a la descarga del pozo
- P
1B
= Bomba centrfuga - Presin en la alimentacin de la bomba

viii

- P
2B
= Bomba centrfuga - Presin en la descarga de la bomba
- P
H
= Hidrociclones - Presin en la alimentacin de los Hidrociclones
- q
i
A
= Harnero - Flujos de agua en la alimentacin
- q
i
B
= Bomba centrfuga - Flujos de agua en la alimentacin
- q
i
C
= Chancador de Pebbles - Flujo msico de agua que entra al chancador
- q
i
H
= Hidrociclones - Flujos de agua en la alimentacin
- q
i
M
= Molino de Bolas - Flujos de agua en la alimentacin
- q
i
S
= SAG - Flujo de agua en la alimentacin
- q
i
P
= Pozo - Flujo de agua en la alimentacin
- q
o
A
= Harnero - Flujo de agua en la descarga
- q
o
B
= Bomba centrfuga - Flujo de agua en la descarga
- q
o
C
= Chancador de Pebbles - Flujo de agua en la descarga
- q
o
H
= Hidrociclones - Flujos de agua en el rebose
- q
o
M
= Molino de Bolas - Flujo de agua en la descarga
- q
o
P
= Pozo - Flujo de agua en la descarga
- q
o
S
= SAG - Flujo de agua en la descarga
- q
r
A
= Harnero - Flujo de agua que recircula
- q
r
H
= Hidrociclones - Flujo de agua en el rechazo
- R= SAG - Matriz de rotacin para pasar de Intervalos de tamaos a acumulados
por malla y viceversa.
-
A
i
A A
R , r , r = Harnero - Flujos de slidos que recirculan
-
C
i
C C
R , r , r = Chancador de Pebbles - Flujos de slidos a la descarga de la fase de
ruptura (matriz B) que recircula hacia la fase de clasificacin (matriz A)
-
H
i
H H
R , r , r = Hidrociclones - Flujos de slidos en el rechazo
-
S
i
S S
R , r , r = SAG Flujo de slidos que recircula dentro del molino. Sale de la
fase de clasificacin para recircular a la cmara de molienda.

ix

-
C
i
C C
R
~
, r
~
, r
~
= Chancador de Pebbles - Flujos de slidos en el rechazo de la etapa
de clasificacin que entra a la etapa de ruptura
-
H
i
H H
R
~
, r
~
, r
~
= Hidrociclones - Fraccin del flujo de slidos que va al rechazo
producto de la clasificacin
-
H
i
H H
R
, r
, r
= Hidrociclones - Fraccin del flujo de slidos que es arrastrado

directamente al rechazo por el agua
- R
CH
= Hidrociclones - Rechazo de Agua
- S= SAG - Potencia elctrica del molino
- v
B
= Bomba centrfuga - Velocidad de la bomba
- Vs= SAG - volumen del molino
-
M
i
M M
W , w , w = Molino de Bolas - Hold-up de slidos
- W
MA
= Molino de Bolas - Hold-up de agua
-
S
i
S S
W , w , w = SAG - Hold-up de slidos
-
P P P
i
w , w , W = Pozo - Hold-up de slidos
-
A
S
A
i
S
A
S
W , w , w = SAG - Hold-up de slidos acumulados por malla
- w
SA
= SAG - Hold-up de agua
- w
SB
= SAG Hold-up de bolas
- Y
H
, y
H
i
= Hidrociclones - Matriz diagonal de clasificacin
o
H
= Hidrociclones - Parmetro de Plitt.
- o
M
= Molino de Bolas - tasa descarga de mineral slido
- o
MA
= Molino de Bolas - tasa descarga de agua
o
S
= SAG - tasa descarga de agua
- o
S
= SAG - Constante de proporcionalidad en la descarga
- c= SAG - Porosidad equivalente del mineral
- q
A
= Harnero - Eficiencia del harnero
- k
A1
, k
A2
, k
A3
= Harnero - Parmetros del intrnsecos del modelo
- k
C1
, k
C2
, k
C3
= Chancador de Pebbles - Lmites para el clculo de las matrices A y

x

- k
H1
, k
H2
, k
H3
= Hidrociclones - Parmetros de los hidrociclones
- k
S1
, k
S2
= SAG - Constantes formula descarga de agua
-
s
= SAG - densidad del material equivalente
-
1
,
2
= Bomba centrfuga - Parmetros de la bomba
-
C1

C12
= Chancador de Pebbles - Parmetros del modelo
-
S1
,
S2
= SAG - Constantes frmula de potencia

xi

B) VARIABLES Y CONSTANTES POR TIPO DE PROCESO
B.1) CONCEPTOS GENERALES SOBRE LA NOTACIN
- El primer subndice en variables y constantes identifica al proceso:
- S : Molienda SAG
- A : Harnero
- C : Chancador de Pebbles
- P : Pozo
- B : Bombas Centrfugas
- H : Hidrociclones
- M : Molienda de Bolas
- El segundo subndice representa el intervalo de tamao a que hace referencia la
variable o constante
- f
Y
= Vector de flujos msicos de slidos de entrada en el proceso Y
- f
Y
i
= Flujo msico de slidos en la entrada de tamao X
i
en el proceso Y
- F
Y
= Suma de los flujos msicos de slidos de entrada en el proceso Y
- p
Y
= Vector de flujos msicos de slidos en la descarga en el proceso Y
- p
Y
i
= Flujo msico de slidos en la descarga de tamao X
i
en el proceso Y
- P
Y
= Suma de los flujos msicos de slidos en la descarga en el proceso Y
- r
Y
= Vector de flujos msicos de slidos que recirculan en el proceso Y
- r
Y
i
= Flujo msico de slidos de tamao X
i
que recircula en el proceso Y
- R
Y
= Suma de los flujos msicos de slidos que recirculan en el proceso Y
- q
i
Y
= Flujo msico de agua de entrada en el proceso Y
- q
o
Y
= Flujo msico de agua en la descarga en el proceso Y
- q
r
Y
= Flujo msico de agua que recircula en el proceso Y
B.2) VARIABLES GENERALES A TODOS LOS PROCESOS
- X
i
= Tamao representativo del intervalo i
-
t
= Densidad de la pulpa de mineral
-
m
= Densidad del mineral (slidos)
-
B
= Densidad de las bolas
- c
P
= Porcentaje de slidos en la pulpa de mineral
- d
50
= Dimetro para el cual pasa el 50% de las partculas
- F
80
= Tamao bajo el cual est el 80% del mineral
- f
D
= El suprandice D representa la fraccin del flujo de mineral duro.
- f
B
= El suprandice B representa la fraccin del flujo de mineral blando.
- o = Tasas de descarga

xii

- P = Presiones
- v = Velocidades
- k, = Constantes
- q = Eficiencias
- d = dimetros
- l = largos
- h = alturas
- A = reas, superficies
- V = volmenes
- n = Cantidades, nmero de elementos
- n = Ruidos
-
1
f
F
= f = Granulometra o flujo unitario de mineral
-
i
f,f ,F = Flujos de slidos de la alimentacin externa (alimentacin fresca) que
ingresa a la Planta de Molienda
- q
i

= Flujos de agua de la alimentacin externa (alimentacin fresca) que ingresa a
la Planta de Molienda
- G = Agua extra que se agrega en la alimentacin de la Planta de Molienda
-
i
r,r ,R = Flujos de slidos en la alimentacin de la planta de molienda que
provienen de recirculaciones de mineral en la planta
- q
r

= Flujo de agua a la alimentacin de la Planta de Molienda que provienen de
recirculaciones de mineral en la planta
B.3) MOLIENDA SEMI-AUTGENA (SAG)
-
S
i
S S
F , f , f = Flujos de slidos de alimentacin al Molino
-
S
i
S S
F
~
, f
~
, f
~
= Flujos de slidos de alimentacin a la cmara de molienda
-
S
i
S S
P , p , p = Descargas de slidos del Molino
-
S
i
S S
P
~
, p
~
, p
~
= Descargas de slidos de la cmara de molienda
-
S
i
S S
R , r , r = Rechazos de slidos dentro del SAG
-
A
S
A
i
S
A
S
F
~
, f
~
, f
~
= Flujos de slidos acumulados de alimentacin a la cmara de
molienda
-
A
S
A
i
S
A
S
P
~
, p
~
, p
~
= Flujos de slidos acumulados en la descarga de la cmara de
molienda
-
S
i
S S
W , w , w = Hold-up por intervalo de tamao

xiii

-
A
S
A
i
S
A
S
W , w , w = Hold-up acumulado por malla
- w
SA
= Hold-up de agua
- q
i
S
= Flujo de agua en la alimentacin
- q
o
S
= Flujo de agua en la descarga
- k, k
i,j
= Constantes de Molienda
- k
A
, k
A
i,j
= Constantes de Molienda acumuladas
- k
E
, k
E
i,j
= Constantes especficas de Molienda
- R = Matriz de rotacin para pasar de Intervalos de tamaos a acumulados por
malla y viceversa.
- S = Potencia elctrica del molino
- w
SB
= Carga de bolas
-
s
= Densidad del material equivalente
- d
S
, l
S
, V
S
= Dimetro, largo y volumen del Molino
-
S1
,
S2
= Constantes frmula de potencia
- J, J
B
= Porcentaje de llenado del mineral y de las bolas
- c = Porosidad equivalente del mineral
- c, c
i
= Tasas de rechazo de mineral
- o
S
= Constante de proporcionalidad en la descarga
- o
S
= Tasa descarga de agua
- k
S1
, k
S2
= Constantes formula descarga de agua
B.4) HARNERO
-
A
i
A A
F , f , f = Flujos de slidos en la alimentacin
-
A
i
A A
P , p , p = Flujos de slidos en la descarga
-
A
i
A A
R , r , r = Flujos de slidos que recircula
- q
i
A
= Flujos de agua en la alimentacin
- q
o
A
= Flujos de agua en la descarga
- q
r
A
= Flujos de agua que recircula
- q
A
= Eficiencia del harnero
- k
A1
, k
A2
, k
A3
= Parmetros del sistema.
- F
ua
= Representa la fraccin de mineral alimentado de tamao menor que la
apertura del harnero.

xiv

B.5) CHANCADOR
-
C
i
C C
-
C i C C
F
~
, f
~
, f
~
= Flujos de slidos en la alimentacin a clasificacin
-
C
i
C C
-
C
i
C C
R , r , r = Flujos de slidos de rechazo que sale de la fase ruptura y recircula
hacia la fase de ruptura
-
C
i
C C
R
~
, r
~
, r
~
= Flujos de slidos de rechazo que entra a ruptura
- q
i
C
= Flujo msico de agua que entra al chancador
- q
o
C
= Flujo msico de agua que sale del chancador
- A, a
i,j
= Matriz de clasificacin y sus elementos
- B, b
i,j
= Matriz de ruptura y sus elementos
- k
C1
, k
C2
, k
C3
= Lmites para el clculo de las matrices A y B
-
C1

C12
= Parmetros del modelo
- d
CE
= Abertura de la descarga (closed setting) del chancador
B.6) POZO
-
P
i
P F
-
P
i
P P
- q
i
P
- q
o
P
- h
P
=Altura de la pulpa en el pozo
- G
P
= agua extra que se agrega al pozo
- A
P
= Superficie del pozo
-
P
i
P P
F , f , f = Flujo volumtrico de pulpa a la alimentacin del pozo
-
P
i
P P
P , p , p = Flujo volumtrico de pulpa a la descarga del pozo
-
P P P
i
w , w , W = Hold-up por intervalo de tamao, del tamao i y total

xv

B.7) BOMBA CENTRFUGA
-
B
i
B B
-
B
i
B B
- q
i
B
- q
o
B
- v
B
= Velocidad de la bomba
- P
1B
= Presin en la alimentacin de la bomba
- P
2B
= Presin en la descarga de la bomba
-
1
,
2
=Parmetros de la bomba
B.8) HIDROCICLONES
-
H
i
H H
-
H
i
H H
P , p , p = Flujos de slidos en el rebose
-
H
i
H H
R , r , r = Flujos de slidos en el rechazo
- q
i
H
- q
o
H
= Flujos de agua en el rebose
- q
r
H
= Flujos de agua en el rechazo
-
H
i
H H
F
~
, f
~
, f
~
= Fraccin del flujo de slidos que va a clasificacin
-
H
i
H H
R
, r
, r
= Fraccin del flujo de slidos que es arrastrado directamente por el

agua al rechazo
-
H
i
H H
R
~
, r
~
, r
~
= Fraccin del flujo de slidos que va a rechazo producto de la
clasificacin
-
H
i
H H
F , f , f = Flujos volumtricos de la pulpa en la alimentacin a los Hidrociclones
- P
H
= Presin en la alimentacin del Hidrocicln
- R
CH
= Rechazo de Agua
- Y
H
, y
H
i
= Matriz diagonal de clasificacin
- n
H
= Nmero de hidrociclones
- h
H
= Altura de los hidrociclones
- d
vf
= Dimetro del buscador de vrtice (vortex finder).
d
sp
= Dimetro del orificio en la descarga (spigot).
o
H
= Parmetro de Plitt.
- k
H1
, k
H2
, k
H3
= Parmetros de los hidrociclones

xvi

B.9) MOLIENDA DE BOLAS
-
M
i
M M
-
M
i
M M
- q
i
M
- q
o
M
-
M
i
M M
W , w , w = Hold-up por intervalo de tamao, del tamao i y total
- W
MA
= Hold-up de agua
- B, b
i
= Matriz triangular inferior de molienda
- o
M
= Tasa descarga de Mineral Slido
- o
MA
= Tasa descarga de agua

- 1 -

1 INTRODUCCIN
La presente memoria tiene como objetivo final realizar un simulador de una planta de
molienda constituida por un circuito de molienda semiautgena (SAG) y de molienda
secundaria con bolas. Se incluye un sistema de control y una interfaz grfica similar a la
de las pantallas de operacin de plantas reales. De esta forma el simulador podr ser
utilizado tanto para investigacin como para capacitacin y entrenamiento.
Los objetivos especficos de esta memoria son:
- Describir los distintos modelos matemticos de los equipos involucrados en este
tipo de planta de molienda.
- Generar un modelo general en variables en tiempo discreto para la creacin de un
simulador.
- Definir e implementar sistemas de control
- Programar un simulador que represente este tipo de planta de molienda
La planta de molienda est compuesta por un circuito de molienda SAG (molino semi-
autgeno (SAG), harnero y chancador de pebbles) y un circuito inverso de molienda
secundaria con bolas (pozo de bombeo, bombas centrfugas, hidrociclones y molinos de
bolas). Cada equipo fue modelado y simulado en forma independiente.
El simulador no pretende representar una planta especfica, sino que una planta genrica
cualitativamente adecuada a la realidad. El simulador permitir modificar los parmetros
de los distintos equipos para analizar distintos casos. De esta forma se espera que los
sistemas de control funcionarn adecuadamente en una planta real, cuando
cualitativamente se comporten bien en una simulada segn estos criterios.
Segn lo anterior, los ajustes de parmetros fueron realizados con informacin de
distintas fuentes. La principal fuente de informacin fueron mediciones realizadas en la
concentradora A2 de Chuquicamata.

- 2 -

Los modelos utilizados representan el comportamiento de la planta durante una operacin
normal y no son vlidos para estados transitorios como partidas de equipos o paradas de
planta. Los ajustes fueron realizados para un punto de operacin, y la validez de la
simulacin es en torno a este.
El simulador fue implementado en Matlab versin 4.1, dado la capacidad y versatilidad
matemtica as como sus posibilidades grficas.
La memoria est estructurada de la siguiente manera:
- Captulo 2: Se explica todo el marco terico de los modelos de cada equipo, se
determina el modelo integrado y se lleva a variables en tiempo discreto. Luego se
describen los controles considerados para la planta.
- Captulo 3: Se detalla cmo se realiz la programacin del simulador
- Captulo 4: Se describe y explica cmo fueron ajustados los parmetros de los
modelos de los equipos involucrados en el proceso de molienda.

- 3 -

2 MODELOS DE LA PLANTA SAG-MOLIENDA SECUNDARIA
2.1 DESCRIPCIN DE LA PLANTA QUE SE QUIERE SIMULAR
La planta a simular consiste en dos circuitos de molienda hmeda: circuito de molienda
semiautgena (circuito 1) o SAG (Semi-Autogenous Grinding) y un circuito de molienda
secundaria con bolas (circuito 2). La Figura 2.1-1 representa la planta de molienda a
simular.

Figura 2.1-1 Circuitos del Simulador
El propsito del circuito de molienda SAG es el de moler el mineral hasta dejarlo de un
tamao mximo definido por el harnero vibratorio (4). El propsito del circuito de molienda

- 4 -

secundaria con bolas es el de entregar una pulpa de mineral de tamao de partculas
adecuado para su procesamiento en una planta de flotacin. Estos dos circuitos estn
conectados en serie, a travs de un pozo (7).
El circuito SAG est conformado por una correa transportadora de alimentacin de
mineral (1), una lnea de agua (2), un molino SAG (3), un harnero vibratorio (4), un
sistema de correas transportadoras de retorno de mineral grueso a molienda (5) y un
chancador de pebbles (6).
El mineral que ingresa al circuito de molienda SAG proviene del proceso de chancado
primario. Al mineral se le agrega agua al ingreso de la molienda SAG. El mineral de
alimentacin al circuito de molienda SAG se denominar alimentacin fresca. Una vez
procesada la alimentacin fresca por el molino, entra en la etapa de clasificacin por el
harnero vibratorio, en donde aquellas partculas que no han alcanzado el tamao
especificado, vuelven al proceso de molienda. Este se puede hacer de dos maneras:
pasando o no por el chancador de pebbles. El material que ha alcanzado el tamao
adecuado, sigue su camino hacia el pozo, en donde entra en el segundo circuito.
El circuito de molienda secundaria est conformado por un pozo (7), una lnea de agua
(8), dos bombas centrfugas (11), dos molinos de bolas (9) y dos bateras de seis
hidrociclones cada una (10). La configuracin mostrada corresponde a un circuito inverso
de molienda secundaria dado que el mineral entra primero a una fase de clasificacin en
los hidrociclones y luego a molienda de bolas.
El pozo recibe el producto de la molienda SAG y las descargas de los molinos de bolas.
En el pozo se agrega agua para formar una pulpa de mineral con una densidad adecuada
para alimentar, mediante las bombas centrfugas (11), los hidrociclones (10). En los
hidrociclones la pulpa de mineral que contiene las partculas de menor tamao (12) es
enviada a los procesos de flotacin y la pulpa con partculas de tamaos mayores es
enviada a molinos de bolas (9) para una mayor conminucin. Finalmente, el producto de
los molinos de bolas vuelve al pozo de bombeo de mineral (7).

- 5 -

2.2 MODELOS DE LAS OPERACIONES UNITARIAS
El modelo de la planta a simular est compuesto de modelos de operaciones unitarias
interrelacionadas mediante sus flujos de entrada y salida. Las operaciones unitarias son:
el molino SAG, el harnero, el chancador de pebbles, el pozo, las bombas centrfugas, los
hidrociclones y los molinos de bolas. En este captulo se analizar cada proceso en forma
independiente y se formularan sus modelos matemticos.
En general se llamar flujos de alimentacin a las entradas de mineral y agua a los
distintos procesos; descargas, a las salidas de mineral y agua de los procesos; y
rechazo, a aquel material, que por alguna razn, recircula en el proceso.
Notaciones adicionales empleadas son las siguientes:
- El primer subndice en variables y constantes identifica al proceso:
- S : Molienda SAG
- A : Harnero
- C : Chancador de Pebbles
- P : Pozo
- B : Bombas Centrfugas
- H : Hidrociclones
- M : Molienda de Bolas
- Con la letra f se representan flujos msicos de alimentacin de mineral slido
- Con la letra q
i
se representan flujos msicos de alimentacin de agua
- Con la letra p se representan descargas de mineral slido
- Con la letra q
o
se representan descargas de agua
- Con las letras r se representan los rechazos
- Las sumas de las componentes de un vector se representan con letras
maysculas
El listado de las variables y constantes utilizadas se encuentra en el Glosario de
Variables y Constantes al final de esta memoria.

- 6 -

Los modelos matemticos empleados son adaptaciones de modelos existentes. En ellos,
tanto las entradas como las salidas corresponden a flujos msicos [t/h] de mineral slido y
de agua. Los flujos msicos de slidos son modelados como vectores de flujos msicos
de partculas que estn dentro de un intervalo de tamaos definido. Cada intervalo de
tamaos est representado por un tamao representativo que corresponde a la media
geomtrica del intervalo. Es decir:
Sea a un vector de flujo msico:
| |
1 2 n
a a a a =
en donde a
1
es el flujo para partculas de tamao representativo X1
en donde
1 2
X1 x x = es el intervalo de las partculas _ tal que:
1 2
x x s _ <
En los molinos (tanto el SAG como el de bolas) los flujos msicos por intervalo de tamao
en la descarga, dependen tanto de los flujos de alimentacin como de la variacin del
tamao del mineral dentro del molino debido a la conminucin. Los tiempos de residencia
del mineral para cada intervalo de tamao, as como las tasas de descarga para cada
tamao dependern de las caractersticas propias del molino. Para representar
matemticamente estos procesos se utilizaron modelos dinmicos como se ver ms
adelante.
En el harnero, chancador de pebbles, bomba centrfuga e hidrociclones, se consideran
despreciables los tiempos de residencia del mineral (para cada intervalo de tamao) por lo
que se utilizan modelos estticos. El pozo es considerado un mezclador perfecto, de
modo que la distribucin granulometra en la alimentacin (mezcla de las descargas del
molino SAG y los molinos de bolas) es igual a la distribucin granulomtrica en la
descarga (alimentacin de las bombas) y por lo tanto tambin se considera un modelo
esttico para los flujos msicos.
2.2.1 Modelo de molino SAG
El Molino SAG est constituido por un cuerpo cilndrico-cnico mayor dimetro que largo,
puesto en forma horizontal y con un sistema de rotacin en torno a su eje horizontal.

- 7 -

Figura 2.2-1 Molino SAG FLSmidth gearless, ring motor ABB
La Molienda Semiautgena se basa en el principio de que las colpas de mineral
contenidas en la alimentacin fresca producen su conminucin as como la de partculas
menores por el impacto entre ellas mismas. En los molinos SAG, a diferencia de los
molinos AG (autgenos), se agregan algunas bolas de acero para ayudar la molienda. En
los molinos de bolas o de barras son estos elementos los que producen la conminucin.
A la alimentacin fresca se le agrega agua. Como se ver ms adelante la descarga de
agua depende del contenido de slidos dentro del molino. En este captulo se analizar
primero en comportamiento de los slidos y luego el del agua dentro del molino.
i) Proceso de los slidos
El proceso de los slidos en el molino SAG se puede dividir en dos etapas:
- Molienda: Esta ocurre cuando dentro del molino, el choque de las colpas entre s
y con partculas menores, ayudadas por el impacto de bolas de acero, reducen
tamao del mineral contenido en el molino (colpas y partculas menores).
Diremos que este proceso ocurre en la cmara de molienda

- 8 -

Clasificacin: El molino tiene en el extremo de descarga una parrilla interna y
alzadores de pulpa. De esta manera el mineral que descarga el molino es una
fraccin del mineral contenido en el. El mineral que es rechazado por la parrilla
permanece dentro del molino continuando su reduccin de tamao.
A continuacin se desarrolla un modelo matemtico que representa el proceso de
molienda que experimentan los slidos dentro del molino SAG. Se utiliza una forma
simplificada de los modelos desarrollados por Weymont [1], Austin [2] y
Barahona [3].
El simulador trabaja con vectores de flujos msicos de mineral por intervalos de tamao.
Se define X
i
como el tamao representativo de partculas que perteneces a un cierto
intervalo de tamaos. Como convencin X
1
es el tamao mayor y X
n
es el menor
(
1 i n
i 1 i n X X X )
Se define a continuacin los flujos msicos de slidos considerados:

Figura 2.2-2 Diagrama de bloques del Molino SAG

S
f : Vector de flujos de slidos por intervalo de tamaos en la alimentacin del
Molino

S
p : Vector de flujos de slidos por intervalo de tamaos en la descarga del
Molino
Cmara de
Molienda
Clasificacin
S
f
S
f
~
S
p
~
S
p
S
r

- 9 -

-
S
f
~
: Vector de flujos de slidos por intervalo de tamaos en la alimentacin a la
cmara de molienda.
-
S
p
~
: Vector de flujos de slidos por intervalo de tamaos en la salida de la
cmara de molienda.
- r
S
: Vector de flujos de slidos por intervalo de tamaos que son rechazados por
la parrilla dentro del Molino
Molienda
Se analiza el proceso de molienda del mineral dentro del molino. Con fines de
modelamiento se dice que este proceso ocurre en la Cmara de Molienda.
Los modelos de Weymont [1], Austin [2] y Barahona [3] utilizan vectores de flujos de
mineral acumulados retenidos por malla. Esto significa que, por ejemplo para la
alimentacin de slidos a la cmara de molienda, el flujo
A
1
S
f
~
representa la suma de los
flujos de mineral de tamaos representativos mayores o iguales a X
1
. Como X
1
es el
tamao mayor,
A
1
S
f
~
contiene slo los flujos de tamao X
1
, es decir
A
S S
1 1
f f =

. De forma
anloga,
A
2
S
f
~
representa la suma de los flujos de mineral de tamaos representativos
mayores o iguales a X
2
(
A
S S S
2 1 2
f f f = +

),
A
S
3
f
representa la suma de los flujos de mineral

de tamaos representativos mayores o iguales a X
3
, etc.
La Figura 2.2-3 muestra en forma esquemtica los flujos retenidos por malla para la
alimentacin a la cmara de molienda.

- 10 -

Figura 2.2-3 Esquema representativo de flujos acumulados retenidos por malla
Se define a continuacin los flujos acumulados por malla en funcin de los flujos por
intervalos de tamao:
- Flujos en la entrada de la cmara de molienda:
1
i
A
S S
k
i k
f f
=
=

] [
A
n
S
A
3
S
A
2
S
A
1
S
A
S
f
~
... f
~
f
~
f
~
f
~
= (2.2-1)
- Flujos de descarga de la cmara de molienda:

=
=
i
1 k
k
S
A
i
S
p
~
p
~
] [
A
n
S
A
3
S
A
2
S
A
1
S
A
S
p
~
... p
~
p
~
p
~
p
~
= (2.2-2)
El mineral acumulado dentro del molino se denomina hold-up y se representa por w
S
=
[w
S
1
w
S
2
w
S
n
] en que w
S
i
es la cantidad de mineral de tamao representativo X
i
que se
encuentra dentro del molino.
X
1
X
2
X
3
......... X
n

S
1
f

S
2
f

S
3
f
S
n
f
A
S S
1 1
f f =

A
S S S
2 1 2
f f f = +

A
S S S S
3 1 2 3
f f f f = + +

n
A
S S
k 1
n k
f f
=
=

- 11 -

Se define el hold-up de mineral acumulado por malla para el tamao i como el hold-up
de mineral con partculas de tamao mayor o igual al tamao i, es decir:
=
=
i
1 k
k
S
A
i
S
w w
1 2 3
A A A A A
S S S S S
n
w w w w ... w = [ ] ( 2.2-3)
En la cmara de molienda se tiene que la variacin en el tiempo del mineral de tamao
superior o igual a X
i
contenido dentro del molino, ser igual al mineral de tamao superior
o igual a X
i
que entre, menos el que salga y menos el mineral de tamao superior o igual
a X
i
que producto de la molienda pase a tamaos inferiores a X
i
. Si k
i
A
es la tasa con que
partculas de tamao i o superior pasan a tamaos inferiores a i, podemos escribir:

A
i
S
A
i
A
i
S
A
i
S
A
i
S
w k p
~
f
~
w
dt
d
= [ 2.2-4]
y en forma matricial se tiene:
A A A A A
S S S S
d
w f p k w
dt
=

[ 2.2-5]
donde k
A
es una matriz diagonal que contiene las tasas de molienda de mineral
acumulado por malla.
Por conservacin de masa se tiene que k
n
A
= 0 ya que, considerando la totalidad del
material, este no puede desaparecer y luego la variacin del hold-up total ser la
diferencia entre los flujos de entrada y salida. Es decir:

A
n
S
A
n
S
A
n
p
~
f
~
w
dt
d
= ( 2.2-6)
Vale decir: La variacin en el tiempo de todo el mineral contenido en el molino es igual al
mineral total que entra al molino menos el mineral total que sale.
Dado que el simulador ocupa flujos de mineral por intervalo de tamaos, se transforma las
ecuaciones obtenidas para flujos acumulados por malla en ecuaciones por intervalos de
tamao:

- 12 -

Para el tamao X
1
:
A
1
S
A
1
A
1
S
A
1
S
A
1
S
w k p
~
f
~
w
dt
d
=
como:
1
S
A
1
S
w w =
1
S
A
1
S
f
~
f
~
=
1
S
A
1
S
p
~
p
~
=
se tiene:

1
S
A
1
1
S
1
S
1
S
w k p
~
f
~
w
dt
d
= ( 2.2-7)
Para el tamao X
2
:
A
2
S
A
2
A
2
S
A
2
S
A
2
S
w k p
~
f
~
w
dt
d
=
como:
2
S
1
S
A
2
S
w w w + =
2
S
1
S
A
2
S
f
~
f
~
f
~
+ =
2
S
1
S
A
2
S
p
~
p
~
p
~
+ =
se tiene:
( ) ( )
2
S
1
S
A
2
2
S
1
S
2
S
1
S
2
S
1
S
w w k p
~
p
~
f
~
f
~
w
dt
d
w
dt
d
+ + + = +
utilizando [2.2-7] se obtiene:
( )
2
S
A
2
1
S
A
1
A
2
2
S
2
S
2
S
w k w k k p
~
f
~
w
dt
d
= ( 2.2-8)
Para el tamao X
3
:
A
3
S
A
3
A
3
S
A
3
S
A
3
S
w k p
~
f
~
w
dt
d
=
como:
3
S
2
S
1
S
A
3
S
w w w w + + =
3
S
2
S
1
S
A
3
S
f
~
f
~
f
~
f
~
+ + =
3
S
2
S
1
S
A
3
S
p
~
p
~
p
~
p
~
+ + =

- 13 -

se tiene:
( ) ( )
3
S
2
S
1
S
A
3
3
S
2
S
1
S
3
S
2
S
1
S
3
S
2
S
1
S
w w w k p
~
p
~
p
~
f
~
f
~
f
~
w
dt
d
w
dt
d
w
dt
d
+ + + + + + = + +
utilizando [2.2-7] y [2.2-8] se obtiene:
( ) ( )
3
S
A
3
2
S
1
S
A
2
A
3
3
S
3
S
3
S
w k w w k k p
~
f
~
w
dt
d
+ = ( 2.2-9)
Para el tamao X
i
:
A
i
S
A
i
A
i
S
A
i
S
A
i
S
w k p
~
f
~
w
dt
d
=
como:
=
=
i
1 k
k
S
A
i
S
w w
=
=
i
1 k
k
S
A
i
S
f
~
f
~

=
=
i
1 k
k
S
A
i
S
p
~
p
~

procediendo de al misma forma que antes se obtiene:

i
S
A
i
1 i
1 k
k
S
A
1 i
A
i
i
S
i
S
i
S
w k w ) k k ( p
~
f
~
w
dt
d
=

(2.2-10)

Para el tamao X
n
:
dado que k
n
A
= 0 se tiene:

=
1 n
1 k
k
S
A
1 n
A
n
n
S
n
S
n
S
w ) k k ( p
~
f
~
w
dt
d
( 2.2-11)
Escribiendo las ecuaciones en forma matricial se obtiene:
=

S S S S
d
w f p k w
dt
(2.2-12)
en donde los elementos de la matriz k tienen los siguientes valores:

- 14 -

>
=
<
=
j i k k
j i k
j i 0
k
A
1 i
A
i
A
i j , i
(2.2-13)
Realizando un manejo algebraico se obtiene la siguiente relacin:

=
1 A
k R k R (2.2-14)
en donde R es una matriz cuyos elementos tienen los siguientes valores:

>
<
=
j i 1
j i 0
R
j , i
(2.2-15)
Hasta el momento se ha desarrollado ecuaciones que representan la variacin del hold-up
de mineral por intervalo de tamaos en funcin de los flujos de alimentacin y descarga, y
tasas de molienda. Falta determinar los valores de las tasas de molienda acumuladas por
malla k
i
A
para poder determinar las tasas k
i
.
Se ha observado experimentalmente que las tasas k
i
A
dependen fuertemente de:
- El flujo de mineral
- El hold-up de mineral
- La potencia elctrica consumida por el motor del Molino
- Caractersticas intrnsecas del mineral, como su dureza
El flujo de mineral y el hold-up de mineral estn relacionados por las ecuaciones antes
desarrolladas. En consecuencia se puede escribir:

S S S
P
~
F
~
W
dt
d
= (2.2-16)
con:
=
S S
w W
=
S S
f
~
F
~

=
S S
p
~
P
~

- 15 -

L. G. Austin [4] determin que el producto entre el material retenido dentro del molino y las
tasas de molienda son directamente proporcionales a la potencia consumida por el
molino.
Si se denomina k
E
al vector con las constantes de proporcionalidad, entonces:
=
A E
S
W k S k (2.2-17)
En que S es la potencia consumida por el motor del molino.
Luego:
=
A E
i i
S
S
k k
W
(2.2-18)
Las constantes k
E
dependen slo de caractersticas del mineral y sern identificadas para
cada caso en forma experimental. Las constantes k
E
se denominan: tasas especficas de
molienda.
Cuando el molino gira con velocidad constante, la potencia elctrica que este consuma
depender de la carga interna del molino. Austin [4] desarroll una frmula emprica para
determinar la potencia consumida por el motor en funcin de la carga de mineral, agua y
bolas; y de caractersticas intrnsecas del molino y del material que se est procesando.
Las ecuaciones se basan en balances de torque en torno al eje de rotacin del molino. Se
considera que el material dentro del molino es de un solo tipo. Este material equivalente
se determina por el hold-up de mineral, de agua y de bolas. La densidad de este material
equivalente ser:

S SA SB
S
eq
W w w
V
+ +
= (2.2-19)

- 16 -

En donde:
- w
SA
: hold-up de agua
- w
SB
: carga de bolas
- V
eq
: Volumen que ocupa el material equivalente
Si V
S
es el volumen interno del molino se puede escribir esta ecuacin como sigue:

+ +
=
S SA SB
S
S
W w w
J V
(2.2-20)
En que J es la fraccin del volumen del molino ocupado por el material equivalente.
Asumiendo que el molino es cilndrico la ecuacin emprica para el clculo de la potencia
consumida por el motor elctrico del molino, desarrollada por Austin, es la siguiente:
=
2.5
S1 S S S2 S
S d l J (1 J) (2.2-21)
donde:
- d
S
: dimetro del molino
- l
S
: Largo del molino
-
S1
,
S2
: Parmetros identificables para cada molino
Para determinar J, Austin desarroll una frmula emprica en que relaciona W
S
con el
material equivalente y el espacio que ocupa este en el molino. La ecuacin es la
siguiente:

| |
1 1 0 4 = c
S S m B m
W V ( ) J ( . ) J (2.2-22)

- 17 -

En que:
- J
B
: Porcentaje de carga de bolas, el que es conocido porque se sabe cuntas bolas
fueron puestas en el molino. No se considera en este modelo el consumo de bolas.
- c : Porosidad equivalente del mineral
-
m
: Densidad del mineral
Se considera una porosidad estndar de 0.4 para las bolas.
Luego se tiene:

1
0 6
1
| |
= +
|
c
\ .
S
B
m S
W
J . J
V
(2.2-23)
Clasificacin - descarga de slidos
Del mineral que se enfrenta a la parrilla interna del molino, una parte ser descargada del
molino y otra ser rechazada y permanecer en este.
Se define las tasas de rechazo por intervalo de tamao c
i
de manera que:
=

S i S
i i
r c p (2.2-24)
=

S i S
i i
p (1 c ) p (2.2-25)
En donde: r
S
i

- r
S
i
: Flujo de rechazo de mineral para el tamao representativo Xi
- p
S
i
: Flujo de descarga del molino para el tamao representativo X
i

- 18 -

Las ecuaciones matriciales son:
=

S S
r c p (2.2-26)
( ) I =

S S
p c p (2.2-27)
Para calcular las tasas de descarga se utiliza la siguiente relacin: la descarga unitaria de
la cmara de molienda es directamente proporcional al hold-up unitario del molino.
=
S S
S S
1 1
p w
W P
(2.2-28)
con:
=
i
S S
p
~
P
~

Empricamente se determina la siguiente relacin:
= o
S S S
P W (2.2-29)
En que o
S
es una constante identificable para cada molino.
ii) Proceso del agua
De manera similar a como fueron establecidas las ecuaciones para los slidos, se puede
establecer la ecuacin que rige al estado del hold-up de agua:

S
o
S
i SA
q q w
dt
d
= (2.2-30)
donde q
i
S
es el flujo de entrada de agua, y q
o
S
el flujo en la descarga.
Como el agua no es fraccionable, en la forma aplicada a los slidos, no se tiene una
ecuacin matricial. Por ello mismo tambin carece de sentido definir cambios, por lo que
no aparece ninguna matrz k de transicin de estados. En todo lo dems, se aprecia la

- 19 -

similitud de las formas de esta ecuacin y la ecuacin de estado para los slidos. w
SA

representa, entonces, el hold-up de agua en el molino.
La descarga de agua, q
o
S
, es proporcional al hold-up de agua en el molino w
SA
. Sea o
S

esta tasa de descarga de agua.
= o
o S SA
S
q w (2.2-31)
La tasa de descarga de agua, se obtiene empricamente de la ecuacin:

2
1
4
S
S S
S
W
k
o = k + (2.2-32)
en que k
S1
, k
S2
son constantes identificables para cada molino.
iii) Ecuacin de Estado y Entrada-Salida
Un sistema dinmico descrito mediante ecuaciones de Estado y Entrada-Salida permite
conocer inequvocamente la evolucin en el tiempo de este, conociendo el estado inicial
del sistema y las variables de entrada.
Sea:
X: Vector de variables de estado del sistema
U: Vector de variables de entrada del sistema
Y: Vector de variables de salida del sistema
u: Vector de constantes del sistema
El sistema se dice que se puede representar mediante ecuaciones de estado y Entrada-
Salida si este se puede escribir de la siguiente forma:
( )
d
dt
=
X
F X, U, u
( ) = Y H X, U, u

- 20 -

El modelo del molino SAG es un sistema de ecuaciones de estado y Entrada-Salida, en
este punto ordenaremos las ecuaciones vistas en los puntos anteriores para mostrar en
forma explcita este hecho.
La ecuacin (2.2-12) representa el comportamiento del hold-up de mineral dentro de la
cmara de molienda del molino SAG, esta ecuacin es:
=

S S S S
d
w f p k w
dt

Del diagrama de la figura 2.2-2 se deduce que el mineral que entra a la cmara de
molienda es igual a la suma de la alimentacin fresca al molino SAG ms el rechazo
mineral de la parrilla interna del molino:

S S S
f f r = +
(2.2-33)
La ecuacin (2.2-26) relaciona el rechazo del harnero interno con el mineral en la
descarga de la cmara de molienda del molino, mediante tasas de clasificacin.
=

S S
r c p
La ecuacin (2.2-28) muestra la relacin directa entre la descarga de la cmara de
molienda y el hold-up de mineral.
=
S S
S S
1 1
p w
W P

donde:
=
S S
w W
La relacin emprica (2.2-29) entrega el valor de la descarga total de la cmara de
molienda.
= o
S S S
P W

- 21 -

Reuniendo las ecuaciones anteriores tenemos:
o o
= +

S S S
S S S S
S S
dw
f c w w k w
dt
w w
(2.2-34)
Utilizando las ecuaciones 2.2-14 y 2.2-17 se obtiene:
1
o o
= +

E S S S
S S S S
S
S S
dw S
f c w w R k R w
dt w
w w

( )
1
I

(
o
(
= +
(

E S S
S S
S
S
dw S
f c R k R w
dt w
w
(2.2-35)
La potencia S consumida por el motor del molino se puede determinar por la ecuacin
2.2-21:
=
2.5
S1 S S S2 S
S d l J (1 J)
La densidad del material equivalente del molino
S
y la fraccin volumtrica en el molino
J, estn dadas por (2.2-20) y (2.2-23).
+ +
=
S SA SB
S
S
W w w
J V

1
0 6
1
| |
= +
|
c
\ .
S
B
m S
W
J . J
V

Luego se ve que la potencia del molino es funcin del hold-up total de slidos W
S
, del
hold-up de agua w
SA
y de las constantes:
S1
,
S2
, d
S
, l
S
, V
S
, w
SB
, c,
m
y J
B

- 22 -

Luego se tiene:
( )
( )
1
I

(
o
(
= +
(

S SA
E S S
S S
S
S
S w , w
dw
f c R k R w
dt w
w
(2.2-36)
Por otro lado se tiene que la variacin del hold-up de agua en el molino es igual a la
diferencia entre el flujo de agua que entra al molino y la que sale de este, este hecho est
representado por la ecuacin 2.2-30.

S
o
S
i SA
q q w
dt
d
=
Utilizando las ecuaciones 2.2-31 y 2.2-32 podemos escribir la ecuacin de estado del
hold-up de agua de la siguiente forma:

( )
2
1
4
S
SA i S SA
S
S
d
w q w
dt
w
| |
k
|
= k +
|
\ .
(2.2-37)
Sea:
X
S
: Vector de variables de estado del Molino SAG
U
S
: Vector de variables de entrada del Molino SAG
Y
S
: Vector de variables de salida del Molino SAG
u
S
: Vector de constantes del Molino SAG

S
S
SA
w
w
(
(

X =
S
S
i
S
f
q
(
(
(

U =
S
S
o
S
p
q
(
(
(

Y =

1 2 1 2
T
E
S S S S S S S S S SB m B
c R k d l V w J ( = o k k c

u

- 23 -

Luego:
( )
S
S S S S
d
dt
=
X
F X , U , u
( )
( )
( )
( )
1
2
1 4
I

( (
o
( (
+
( (
(
(
= =
(
(
| |

k
(
|
k +
(
|
(
\ .
X
= F X , U ,
S SA
E S
S S
S
S
S
S
S S S S
SA
S
i S SA
S
S
S w , w
f c R k R w
w
w
w
d d
w dt dt
q w
w
u (2.2-38)
La ecuacin (2.2-27) representa la descarga de slidos del molino:
( ) I =

S S
p c p
Y dado que segn (2.2-28) y (2.2-29), se tiene:
=
S S
S S
1 1
p w
W P

= o
S S S
P W
Se concluye que:
( ) I
S
S S
S
p c w
w
o
=
(2.2-39)
La descarga de agua est dada por las ecuaciones (2.2-31) y (2.2-32).
o S SA
S
q w = o
2
1
4
S
S S
S
W
k
o = k +

- 24 -

Luego:

( )
2
1 4
S
o S SA
S
S
q w
w
| |
k
|
= k +
|
\ .
(2.2-40)
Finalmente obtenemos:
( )
( )
( )
2
1 4
I
S
S
S
S
S S S S
o
S
S SA
S
S
c w
w
p
q
w
w
o (
(
(
(
( = = =
(
| |
( ( k

|
k +
(
|
(
\ .
Y H X , u (2.2-41)

- 25 -

2.2.2 Modelo de Harnero
Para este modelo, se emplea la adaptacin de A. Casali al modelo de Karra. Primero se
utiliza en el mbito del procesamiento de los slidos que pasan por el harnero, y luego se
extienden sus resultados al caso del agua.
En esencia el harnero es una malla con un reticulado tal que permite el paso de las
partculas que tengan un tamao menor o igual al de la ranura.
Los flujos msicos utilizados en la modelacin son los siguientes:

Figura 2.2-4 Ejemplo de harnero y diagrama de bloques
- f
A
: Vector de flujos msicos por intervalo de tamaos en la alimentacin al
harnero
- p
A
: Vector de flujos msicos por intervalo de tamaos en la descarga del harnero
- r
A
: Vector de flujos msicos por intervalo de tamaos en el rechazo del harnero
- q
i
A
: Flujo msico de agua en la alimentacin al harnero
- q
o
A
: Flujo msico de agua en la descarga del harnero
- q
r
A
: Flujo msico de agua en el rechazo del harnero
Debido a la estructura plana de la malla del harnero, la seleccin del mineral que
atravesar solo puede hacerse en dos dimensiones. Por esta razn pueden existir
partculas en que dos de sus dimensiones sean menores que las ranuras del harnero y
Clasificacin
q
A
A
f
A
p
A
r
q
i
A
q
o
A

q
r
A

- 26 -

aunque su tercera dimensin sea mayor, estas pasaran el harnero. As tambin, puede
ocurrir que partculas que tienen todas sus dimensiones menores que las dimensiones de
las ranuras, no pasen el harnero por chocar con los bordes de la malla o con otras
partculas. Luego el modelamiento deber ser probabilstico.
Una manera de evitar tener que lidiar con los aspectos formales de la teora de la
probabilidad, es observar los resultados de experiencias realizadas en terreno. En este
campo, entramos en los dominios de la estadstica, en donde se pueden encontrar
representaciones que simplifican la modelacin.
Una de las simplificaciones que se usa tiene que ver con la eficiencia, definida como la
proporcin de rocas que cumpliendo con un determinado tamao no pasa por el harnero.
Esto puede deberse a que la roca no tiene el tamao definido por el reticulado o est
siendo obstaculizada en su intento de traspaso de la malla.
Karra propone la siguiente expresin para la ecuacin de eficiencia para el tamao i:

A2
i
A A1
50
i
X
1 exp
d
k
(
| |
( q = k
|
(
\ .

(2.2-42)
En que el d
50
es dimetro para el cual pasa el 50% de las partculas. Para la
determinacin de esta variable se considera la simplificacin de Casali:

0.148
A ua
50
A3
F F
d
1.6
| |
=
|
k
\ .
(2.2-43)
donde:
- k
A1
, k
A2
, k
A3
son parmetros del sistema.
- F
ua
representa la fraccin de mineral alimentado de tamao menor que la
apertura del harnero.
De este modo, la cantidad de rocas que no pasa por el harnero, o rechazo es:

A A A
i i i
r f = q (2.2-44)
Luego, la salida del harnero ser:

- 27 -

p
A
= f
A
r
A
(2.2-45)

o escrito para cada tamao:

i
A
i
A
i
A
r f p = (2.2-46)
Dado que los slidos y lquidos entran y salen mezclados en el proceso, es que las
cantidades de slidos y agua deben ser proporcionales. Dado que la mayor parte del
agua alimentada se va a la descarga, en el rechazo del harnero slo queda el agua
correspondiente a la humedad del mineral, que se estima en un 10% en promedio y por lo
tanto:

n
r A
k 1
A k
q 0.1 r
=
=

(2.2-47)
y por lo tanto, la salida de agua es:

A
r
A
i
A
o
q q q = (2.2-48)

- 28 -

2.2.3 Modelo de Chancador de pebbles
Un chancador es un equipo que tiene una cavidad con un eje cnico excntrico que tritura
el mineral mediante el choque de este con las paredes de la cavidad. En la siguiente
Figura 2.2-5 se esquematiza este proceso:

Figura 2.2-5 Corte de un Chancador de pebbles y esquema funcional
El modelo empleado para el Chancador es el de Andersen & Whiten que ha sido
adaptado por Casali [5].
En el chancador se ha modelado considerando dos procesos independientes: clasificacin
de mineral, y ruptura (Figura 2.2-6).

- 29 -

Figura 2.2-6 Diagrama de bloques del Chancador
Los flujos msicos utilizados en la modelacin son los siguientes:
- f
C
: Vector de flujos de slidos en la entrada al chancador
-
C
f
~
: Vector de flujos de slidos en la entrada a clasificacin
- p
C
: Vector de flujos de slidos en la descarga del chancador
- r
C
: Vector de flujos de slidos de rechazo que sale de la fase ruptura y recircula
hacia la fase de ruptura
-
C
r
~
: Vector de flujos de slidos de rechazo que entra a ruptura
- q
i
C
: Flujo msico de agua que entra al chancador
- q
o
C
: Flujo msico de agua que sale del chancador
Del mineral que entra al chancador una porcin de este pasar directamente a la
descarga y otro tanto ser fraccionado. Como se ve, al contrario de lo que se tiene en el
Clasificacin
A
Ruptura
B
C
f
C
f
~
C
p
C
r
Chancador
C
r
~

- 30 -

molino SAG, en el chancador comprende primero por un proceso de clasificacin y luego
por uno de ruptura.
El material que entra a la fase de clasificacin
C
f
~
es igual a la suma del flujo de mineral
fresco que entra al chancador producto del rechazo del harnero f
C
(f
C
= p
A
) ms el flujo
de mineral que no fue descargado por el chancador y que pas por ruptura r
C
.

C C C
r f f
~
+ = (2.2-49)
El flujo msico de mineral que no pasa a la descarga se calcula en funcin de tasas de
descarga. Sea A la matriz diagonal con las tasas de descarga por intervalo de tamaos.
Luego el flujo msico de mineral que no pas a descarga es:

C C
r A f =

(2.2-50)
El mineral rechazado, pasar a una fase de ruptura. Si B es la matriz que contiene las
eficiencias de ruptura, se tiene que el vector de rechazo de mineral r
C
vale:

C C C
r B r B A f = =

(2.2-51)
reemplazando en (2.2-4)] se tiene:

C C C
f f B A f = +

( )
C C
B A f f I =
(2.2-52)
( )
1
C C
f B A f
= I

luego la ecuacin para el flujo msico de descarga del chancador es:

C C C
p f A f =

( ) ( )
1
C C
p A B A f
= I I (2.2-53)

- 31 -

La clasificacin del material depende del flujo de mineral de entrada y de caractersticas
fsicas del Chancador.
Para el clculo de las tasas de clasificacin para cada intervalo de tamao (las que estn
contenidas en la matriz A) experimentalmente se ha determinado la siguiente ecuacin:

k <
k < < k
k k
k
k >
=
1 C i
2 C i 1 C
1 C 2 C
i 2 C
2 C i
i
X 0
X
X
1
X 1
a
(2.2-54)
Las variables auxiliares k
C1
y k
C2
se calculan con las siguientes ecuaciones:
k
C1
=
C0
+
C1
d
CE
+
C2
F
C
+
C3
F
80
(2.2-55)
k
C2
=
C4
+
C5
d
CE
+
C6
F
C
+
C7
F
80
(2.2-56)
Donde:
- d
CE
: abertura de descarga (closed setting) del Chancador en [mm]
-
=
i
C C
f F : flujo total de entrada al Chancador
- F
80
: tamao representativo X
i
tal que la suma de los flujos de alimentacin de
tamaos inferiores o iguales a X
i
constituyen el 80% de la alimentacin total
del Chancador
-
C0

C7
: constantes caractersticas del Chancador
En la fase de ruptura se supone que en el flujo msico de partculas con tamao
representativo X
i
en la salida de la fase de ruptura r
C
i
est compuesta slo por el aporte
de partculas del flujo de entrada
C
r
~
de tamaos mayores a X
i
. Luego se tiene que no
hay partculas del tamao mayor X
1
en el rechazo del chancador. Adems, por ley de
conservacin de la materia, se tendr que el flujo msico de las partculas de los ltimos
dos tamaos X
n-1
y X
n
en
C
r
~
pasarn ntegramente a r
C
. Es decir, la matriz de
ruptura est determinada por:

- 32 -

> =
< < =
< < < <
= < <
s
=
1 n j n i 1
1 n j 1 n i b 1
i j 1 n i 2 b
1 j n i 1
j i 0
b
1 n
1 k
j , k
) 1 j ( ), 1 i (
1 i i
j , i
(2.2-57)
1 2
2 3 1 2
3 4 2 3 1 2
4 5 3 4 2 3 1 2
2 1 3 2 4 3 5 4 2 3 1 2
1 1 1 2 1 3 1 4 1 3 1 2
0 0 0 0 0 0
0 0 0 0 0
0 0 0 0
0 0 0
0 0
0 0
1 1 1 1 1 1 1 1
n n n n n n n n
n n n n
B
( ) ( ) ( ) ( ) ( ) ( )

(
(
(
(
(
(
(
(
(
(
(
(
(

En que:

3 C
8 C
1
i
i
X
X
1 1
k
(
(
|
|
.
|
\
|
= (2.2-58)
k
C3
=
C9
+
C10
d
CE
+
C11
F
C
+
C12
F
80
(2.2-59)

La mayor parte del agua contenida en la pulpa de mineral que sale del molino SAG pasa
el harnero junto con los slidos menores hacia el pozo, luego el agua que entra al
chancador es mnima y est poco adherida a los pebbles y por lo tanto la mayor parte
pasa directamente a la descarga.
Por lo anterior el modelo supone que el flujo de agua que entra al chancador es igual al
flujo de agua que sale y que no existe ningn tipo de recirculacin. Es decir:

C
i
C
o
q q = (2.2-60)

- 33 -

2.2.4 Modelo de Pozo
El circuito de molienda que se est simulando contiene un pozo de acumulacin de
mineral (Figura 2.2-7), segn se describi en 2.1. Los flujos msicos considerados en este
modelo son los siguientes:

Figura 2.2-7 Diagrama de bloques del Pozo
- f
P
: Vector de flujos de slidos en la entrada al pozo
- p
P
: Vector de flujos de slidos en la salida del pozo
- q
i
P
: Flujo msico de agua que entra al pozo
- q
o
P
: Flujo msico de agua que sale del pozo
- G
P
: Flujo msico de agua externo que se agrega al pozo
i) Proceso de la pulpa de mineral
A diferencia de los modelos vistos anteriormente, en el modelo del pozo no tiene sentido
separar los procesos de slidos y del agua. Todos los fenmenos que ocurren en el pozo
consideran una pulpa de mineral uniforme.
Si
m
es la densidad del mineral de entrada al pozo, entonces el flujo volumtrico de
entrada al pozo es el siguiente:
Pozo
q
i
P

p
P

G
P

- 34 -

P
P i P
m
P
F
f q G = + +
(2.2-61)
donde:

P P
F f =

Si A
P
es el rea transversal equivalente del pozo, en el punto de operacin, y h
p
es la
altura o nivel de la pulpa en el pozo, la variacin en el tiempo del volumen de pulpa en el
pozo, ser igual al flujo volumtrico de pulpa en la alimentacin al pozo menos el flujo
volumtrico de pulpa en su descarga.
Llamando
P
p al flujo volumtrico de salida del pozo, se tiene:

p
p P P
dh
A f p
dt
= (2.2-62)
Si
t
p
es la densidad de la pulpa en el pozo, la ecuacin con flujos msicos es la siguiente:

( ) p p t P i P P t
P P P
d
A h F q G p
dt
= + + (2.2-63)

Desarrollando:
t
p
p p p t P i P P t
P
P P P
d
dh
A h A F q G p
dt dt
+ = + +
Reemplazando el valor de
dt
dh
p
, en (2.2-62), se tiene:
t
P
p p p t i P P P i P P t
p m
P
P P P P
d
F 1
A h A q G p F q G p
dt A
| |
+ + + = + +
|
\ .

- 35 -

Simplificando y despejando
dt
d
P
t
se tiene:
( )( )
t t
P t i P
p p m
P P
P P
d
1
1 F 1 q G
dt A h
( | |
= + + ( |
|
(
\ .
(2.2-64)
El porcentaje de slidos en la pulpa es:
1
1
t
m
p
m t
P
P
c

=

(2.2-65)
Luego los flujos msicos de slidos y agua en la descarga del pozo son:
P P t p
P
P p c = (2.2-66)
o P t p
P P
q p (1 c ) = (2.2-67)
La distribucin granulomtrica del flujo msico de slidos en la descarga del pozo p
P
es la
misma distribucin que tiene el material contenido en el pozo w
P
y est dada por la
siguiente funcin de estado:
P P P
d
w f p
dt
= (2.2-68)

- 36 -

2.2.5 Modelo de Bomba Centrfuga
En el simulador se considera una bomba centrfuga de velocidad variable que extrae la
pulpa del pozo y la eleva hasta las bateras de hidrociclones, como se ilustra a
continuacin (Figura 2.2-8):
h
p P
1
P
2
P
H
h
H

Figura 2.2-8 Circuito de la bomba centrfuga
El modelo de la bomba considera flujos volumtricos y presiones. En la bomba el flujo
volumtrico de alimentacin es igual al de descarga
B
f y es funcin de las presiones en
la alimentacin P
B1
y en la descarga P
B2
de la bomba, y de la velocidad v
B
de la
bomba.

Figura 2.2-9 Diagrama de bloques de la bomba
Bomba
k
B1
, k
B2
, v
B
P
B1
P
B2

B
f
B
f

- 37 -

El flujo volumtrico de pulpa
B
f que impulsa la bomba es proporcional a la velocidad de la
bomba e inversamente proporcional a la diferencia entre las presiones de alimentacin y
de descarga:

| |
B B1 B B2 B2 B1
f = k v k P P (2.2-69)
En donde k
B1
y k
B2
son constantes de proporcionalidad propias de cada bomba.
La presin en la alimentacin P
B1
es causada por la columna de pulpa en el pozo. La

densidad de la pulpa es la densidad de esta en el pozo de mineral
t
p
, luego se tiene:

1 B t P
P
g h = P (2.2-70)
donde g es la aceleracin de gravedad.
La presin en la descarga de la bomba ser igual a la suma de la presin ejercida en la
alimentacin a los hidrociclones, P
H
, ms el peso de la columna de pulpa en las caeras
que alimentan a los hidrociclones, ms las prdidas por roce de la pulpa en las caeras
P
B3
y ms las prdidas singulares P
B4
causadas por el enfrentamiento del flujo de pulpa
con quiebres (como codos) en las caeras. Es decir:

B2 H t H B3 B4
P
g h = + + + P P P P (2.2-71)
Las prdidas por friccin dependen de la velocidad de la pulpa v
P
, del dimetro d
T
y
longitud L
T
de la tubera, y de un factor de prdidas de friccin k
TF
que es caracterstico
de la tubera. La ecuacin que representa ests prdidas es:

2
T P
B3 TF
T
L
d 2 g
v
= k P (2.2-72)
Las prdidas por singularidades dependen de la velocidad de la pulpa v
P
y de las
caractersticas topolgicas del sistema de tuberas, las que se representan en un factor de
prdidas singulares k
TS
. La ecuacin que representa ests prdidas es:

- 38 -

2
P
B4 TS
2 g
v
= k P (2.2-73)
Reemplazando se obtiene:

2 2
T P P
B2 H t H TF TS
T
P
L
g h
d 2 g 2 g
v v
= + + k + k P P (2.2-74)
Si A
T
es el rea transversal de la caera, entonces la velocidad de la pulpa est dada por:

B
P
T
f
A
v = (2.2-75)
Luego:

2 2
T B B
B2 H t H TF TS
T T T
P
L f f 1 1
g h
d 2 g A 2 g A
| | | |
= + + k + k
| |
\ . \ .
P P

2 TS T TF
B2 H t H B 2 2
T T T
P
L 1
g h f
A d 2 g A 2 g
| | k k
= + + +
|
\ .
P P (2.2-76)
Sea:

TS T TF
B3 2 2
T T T
L 1
A d 2 g A 2 g
k k
k = + (2.2-77)
en que k
B3
slo depende de parmetros geomtricos y las constantes del sistema.
Obtenemos:

2
B2 H t H B3 B
P
g h f = + + k P P (2.2-78)

- 39 -

Reemplazando en (2.2-60) el valor de las presiones se obtiene el valor de la presin en la
entrada a los hidrociclones:
( ) ( )
2
B B1 B B2 H t H B3 B t P
P P
f g h f g h
(
= k v k + + k

P
( )
2
B B1 B B2 H B2 t H P B3 B
P
f g h h f
(
= k v k k + k

P
( )
2 B1 B
H B t H P B3 B
B2 B2
P
f
g h h f
k
= v + k
k k
P (2.2-79)

- 40 -

2.2.6 Modelo de Hidrocicln
i) Descripcin fenomenolgica
El Hidrocicln es un clasificador de mineral muy eficiente en partculas pequeas. Su
principio de funcionamiento se basa en aplicar fuerza centrfuga a las partculas entrantes.
En la Figura 2.2-10 se muestra una batera de hidrociclones y un esquema, como se
puede apreciar el hidrocicln est compuesto de un cuerpo cnico abierto unido en la
parte superior con una seccin cilndrica que tiene una entrada en involuta de
alimentacin de mineral. En la parte inferior del cono est el orificio de salida del mineral
grueso, el que se denomina Spigot. La parte superior de la seccin cilndrica est
cerrada excepto en el centro donde entra axialmente una caera de descarga del mineral
fino. La caera se extiende hacia dentro del hidrocicln por medio de una seccin corta y
removible llamada vortex finder, la que previene que el mineral de alimentacin pase
directamente a la descarga de finos.

Figura 2.2-10 Batera de hidrociclones y esquema del hidrocicln
La pulpa de mineral ingresa tangencial a las paredes del Hidrocicln y con alta presin
(sobre 7 psi) producindose un vrtice dentro de este. Las partculas de mineral se ven
sometidas a fuerzas centrfugas crendose una zona de baja presin alrededor del eje
vertical (ver Figura 2.2-11).
Spigot
Alimentacin
de Mineral
Vortex Finder
Descarga
de Mineral
Descarga de
Mineral Grueso

- 41 -

La fuerza centrfuga ser mayor sobre las partculas de mayor tamao, ya que aunque la
densidad de todas las partculas no es la misma, como el tamao incide mucho ms que
la densidad en la masa de las partculas (exponente 3 v/s exponente 1), luego las
partculas mayores tendrn mayor masa. De esta forma, en la zona cercana a las paredes
cnicas del Hidrocicln habr preferentemente partculas de tamaos mayores ya que
pesan ms, adems esta zona est bajo alta presin y la velocidad tangencial es menor,
por lo cual estas partculas escurren hacia la parte inferior del Hidrocicln.
Alrededor del eje vertical del Hidrocicln se genera una zona de baja presin con
partculas de menor tamao, ms livianas y con mayor velocidad. Cuando el Hidrocicln
est funcionando en forma adecuada, la baja presin en la zona axial, induce que aire
ingrese por el spigot para igualar presiones, generando un vrtice, el que arrastra las
partculas menores hacia la salida superior. El mineral grueso saldr por los bordes
externos del spigot, quedando la columna de aire en el eje axial del Hidrocicln. La
columna de aire est compuesta por aire que ingresa desde el exterior ms el aire disuelto
en la pulpa de mineral.
El flujo de mineral grueso que sale por la parte inferior del Hidrocicln se denomina
rechazo. El mineral fino que sale por la parte superior del Hidrocicln se denomina
rebose.

Figura 2.2-11 Funcionamiento del Hidrocicln
Rechazo
C
o
l
u
m
n
a

d
e

A
i
r
e

Alimentacin
Rebose

- 42 -

ii) Modelo utilizado en el simulador
En el simulador se utiliza el modelo de Hidrocicln desarrollado por Lynch y Rao [6] con la
ecuacin de eficiencia de Plitt [7]. En este se consideran por separado los flujos de slidos
y de agua. El diagrama en bloques del modelo de Hidrocicln es el siguiente:

Figura 2.2-12 Circuito y diagrama de bloques del Hidrocicln
De acuerdo a los procesos cinticos descritos en el captulo anterior, la densidad de agua
en la pulpa es distinta en las diferentes zonas interiores del hidrocicln. No todo el mineral
fino es arrastrado hacia el rebose, mineral fino disuelto en una pulpa con un cierto
contenido de agua mnimo es conducido al rechazo. El modelo representa este hecho
separando el mineral en el rechazo en dos, uno que slo es funcin de la densidad de
agua en la pulpa (el que se dice que es arrastrado por el agua) y otro que es producto
de clasificacin.
De esta forma el mineral que es arrastrado por el agua al rechazo es:

H CH H
r f = R (2.2-80)
R
CH

Y
H

H
r
H
r
r
r
H
q
H
f

i
H
q
H
f

H
p
o
H
q
H
f

i
H
q
H
r
r
H
q
H
p
o
H
q

- 43 -

En donde R
CH
es el Rechazo de Agua o Cortocircuito de Agua y est definido como el
cuociente entre el flujo de agua en la alimentacin q
i
H
y en el rechazo q
r
H
.
=
r
CH
i
H
H
q
q
R (2.2-81)
Se tiene que:
q
r
H
=
q
i
H
-
q
o
H
(2.2-82)
En forma emprica [6] se encontr una relacin para el flujo de agua en el rebose, la cual
depende del porcentaje de slidos en la alimentacin c
p
FH
y el nmero de Hidrociclones
(n
H
).

2
1 1 3 96 = + k
FH
o i P H H
H H
q . q ( . c )n (2.2-83)
Donde k
H2
es un parmetro identificable para cada Hidrocicln.
Luego el mineral que va a clasificacin es:
1
H H H CH H
f f r ( )f = =
R (2.2-84)
El flujo msico de mineral que va al rechazo producto de clasificacin es:

H H H
r Y f =

(2.2-85)
En donde Y es la matriz diagonal de clasificacin, y est dada por:

(
(
(
(
(
=
n
H
2
H
1
H
H
y 0 0
0
y 0
0 0 y
Y
(2.2-86)

- 44 -

con [7]:

i
H
H
50
i
x
0.693
d
y 1 e
| |
|
|
\ .
o
= (2.2-87)
Donde:
- x
i
: Tamao representativo del intervalo de tamao i.
- o
H
: Parmetro de Plitt.
- d
50
: Dimetro para el cual el 50% de las partculas que son clasificadas pasan al
rechazo.
Para determinar el d
50
se desarroll [6] la siguiente ecuacin emprica, la que es funcin
de caractersticas fsicas, de la presin de alimentacin al Hidrocicln y del agua en el
rebose:

50 3
6 6 8 9 10 7 47 2
o
sp
vf H
H
H
H
q
d
d
d exp
. . . . n
| |
= + + k
|
|
\ .
P
(2.2-88)
Donde:
- d
vf
: Dimetro del buscador de vrtice (vortex finder).
- d
sp
: Dimetro del orificio de salida (spigot).
- P
H
: Presin de la pulpa de mineral en la alimentacin de los hidrociclones.
- k
H3
: Parmetro del hidrocicln.
La presin de la pulpa de mineral en la alimentacin de los hidrociclones, se determina
por la siguiente ecuacin emprica [6]:

( )
2
0 125
1
13 725
1 48 1
H
H .
vf H p H
f
.
. d n c
| |
|
=
|
k
\ .
P (2.2-89)
En donde k
H1
es un parmetro del sistema y
H
f el flujo volumtrico total de pulpa en la
alimentacin al Hidrocicln. Se puede observar que:

- 45 -

H
H H
H i
m t P
F F
f q
c
= + =

(2.2-90)
El rechazo de mineral, el que ser enviado a los molinos de bolas, es:
H H H
r r r = +
(2.2-91)
En que
m
es la densidad del mineral,
t
es la densidad de la pulpa de mineral en la
entrada a los Hidrociclones (la que es igual al de la salida del pozo) y F
H
es el flujo msico
total de mineral en la alimentacin a los Hidrociclones.
Por ltimo tenemos que el vector de flujo msico de mineral por intervalos de tamaos en
el rebose es:

H H H
p f r =

(2.2-92)

- 46 -

2.2.7 Modelo de Molino de Bolas
El Molino de Bolas est constituido por un cuerpo cilndrico puesto en forma horizontal,
con un sistema que le permite girar en torno a su eje horizontal, en cuyo interior ingresa el
mineral a ser molido.
Dentro del Molino se colocan bolas de acero. Al girar, el molino levanta por sus paredes la
pulpa de mineral y las bolas. Por gravedad estos caen y las bolas al golpear el mineral lo
reduce de tamao.

Figura 2.2-13 Diagrama funcional del Molino de Bolas
El modelo utilizado para el molino de bolas es una adaptacin del modelo utilizado por
Gonzlez [8] [9].

Alimentacin
Descarga
Mineral
Bola de Acero

- 47 -

El diagrama en bloques del molino es el siguiente:

Figura 2.2-14 Diagrama de bloques del Molino de Bolas
i) Proceso de los Slidos
La molienda se modela como un sistema de ecuaciones de estado para el contenido
msico de mineral slido dentro del Molino o hold-up w
M
por intervalo de tamaos X
i
.
La variacin del hold-up para el tamao X
i
en el molino, w
M
i
, ser igual al mineral de
tamao X
i
que ingresa al molino, f
M
i
, menos el mineral de tamao X
i
que sale del Molino,
p
M
i
, menos el mineral de tamao X
i
, contenido dentro del Molino, que a causa de la
molienda pasa a tamaos menores que X
i
y ms el mineral de tamaos mayores que X
i

que a causa de la molienda pasaron al tamao X
i
. Es decir:

i 1
M i,i M i,k M
k 1
M M
i i i i k
d
w f p b w b w
dt

=
= +
(2.2-93)
En forma matricial se expresa como sigue:

M M M M
d
w f p B w
dt
= + (2.2-94)
B es la matriz triangular inferior de molienda que contiene las tasas de molienda por
intervalo de tamaos. En la diagonal se encuentran las tasas de desaparicin del tamao
X
i
, luego estos valores son negativos. Los valores bajo la diagonal representan los aportes
al tamao X
i
de los tamaos superiores y luego estos valores son positivos. b
n,n
es igual a
Molino de Bolas
B, o
M
f
M
p
M
q
i
M
q
o
M

- 48 -

cero (0) ya que el mineral del tamao menor no puede molerse. La matriz B tiene la
siguiente forma:

(
(
(
( =
(
(
(

1,1
2,1 2,2
3,1 3,2 3,3
n,1 n,2 n,3
b 0 0 0
b b 0 0
b b b B
0
b b b 0
(2.2-95)
El Molino de Bolas no posee un harnero interno como en el caso del molino SAG, por lo
que la tasa de descarga de mineral es independiente del tamao del mineral. Se
considera la misma tasa de descarga de slidos o
M
para todos los intervalos de tamao,
luego el vector de flujos msicos de descarga del Molino de Bolas p
M
es:
M M M
p w = o
M M M
i i
p w i = o (2.2-96)
ii) Proceso del Agua
El comportamiento del agua dentro del Molino se trata en forma independiente a los
slidos. Esta forma de modelar el Molino, cuando los parmetros estn bien identificados,
es una buena representacin de la realidad en un amplio rango dentro del punto de
operacin.
La variacin del agua dentro del molino w
MA
ser igual al flujo de agua entrante, q
i
M
,
menos el flujo saliente, q
o
M
. Es decir:

MA i o
M M
d
w q q
dt
= (2.2-97)
El flujo de agua saliente ser proporcional al contenido de agua dentro del Molino. La tasa
de descarga de agua la llamamos o
MA
. La ecuacin de descarga de agua es la siguiente:

MA MA o
M
q w = o (2.2-98)

- 49 -

2.3 INTEGRACIN DE LOS MODELOS: MODELO GLOBAL
Se denomina Modelo Global al modelo que interconecta los modelos individuales,
descritos en los captulos anteriores, de modo que representan la planta de Molienda que
se simular.
Como se vio en 2.1 la planta se puede dividir en dos circuitos de molienda: Molienda
Semiautgena y Molienda Secundaria. Se analizarn en primer lugar la interconexin de
modelos que permite generar un modelo para cada circuito de molienda y luego la
interconexin de estos dos circuitos.
Los modelos de Molienda Semiautgena, Molienda Secundaria y Modelo Global, se
pueden expresar como ecuaciones de estado y ecuaciones de Entrada-Salida, y desde
este punto de vista ser desarrollada la interconexin de submodelos.
2.3.1 Ecuacin de Estado y Entrada-Salida del circuito de Molienda SAG
El circuito de molienda SAG est constituido por un molino SAG que descarga en un
harnero, este clasifica el mineral de forma que el mineral que sale del ciclo de molienda
SAG tenga tamaos inferiores a la apertura del harnero. El mineral rechazado por el
harnero retorna al molino SAG con la opcin de pasar por un Chancador de pebbles. En la
alimentacin del molino SAG se agrega agua a la mezcla segn necesidad. El diagrama
en bloques del circuito de molienda SAG se presenta en la Figura 2.3-1

Figura 2.3-1 Diagrama de Bloques del circuito de Molienda SAG
Molino
SAG
Harnero
Chancador
G
q
o
S

p
S
q
i
A

f
A
q
i
S

f
S
f

q
o
A

p
A
q
i
C

f
C
q
o
C

p
C
q
r
A
r
A
q
r
r

- 50 -

Como se observa, los flujos de slidos y agua en la entrada de cada modelo estn
vinculados de la siguiente forma:
Molino SAG
- f
S
= f + r (2.3-1)
- q
i
S

=

G + q
r

Harnero
- f
A
= p
S
(2.3-2)
- q
i
A

= q
o
S

Chancador (si est conectado)
- f
C
= r
A
(2.3-3)
- q
i
C

= q
r
A

Flujos de realimentacin
-
C
A
Si el Chancador est conectado
Si el Chancador no est conectado
p
r
r
(2.3-4)
-
r
o
C
r
A
Si el Chancador est conectado
Si el Chancador no est conectado
q
q
q

El flujo msico f y el flujo de agua G son variables de entrada que son manejadas por el
operador de la planta.
Se observa de la Figura 2.3-1 que el vector de variables de entrada a la molienda
Semiautgena es:

SAG
f
G
(
(
(
(

U =
Ch

- 51 -

en donde:

1
0
Chancador conectado
Chancador noconectado
si
si
Ch
Todas las variables de salida de los submodelos del modelo de molienda Semiautgena
son de inters, luego el vector de variables de salida queda definido por:

S
S
o
A
A
o
o
A
A
r
SAG
r
C
C
o
o
r
r
S
A
A
A
A
C
C
p
q
p
p
q
q
r
r
q
q
p
p
q
q
r
r
q
q
(
( (
( (
( (
( (
( (
( (
( (
=
( (
( (
( (
( (
( (
( (
( (

(

Y
Y =
El vector de variables de estado es:

S
SAG S
SA
w
w
(
(

X = X =
Determinaremos el valor de las variables de salida.
Harnero
Segn la ecuacin (2.2-44) el rechazo de mineral del harnero est dado por:

i
A
i
A
i
A
f r q =
I
A A A
r f = q

- 52 -

Aqu
T
A
q es el vector con las tasa de eficiencia del harnero para cada intervalo de tamao.
Como se vio en el captulo 2.2.2, estas tasas dependen del d
50
(dimetro para el cual
pasa el harnero el 50% de las partculas) y de los parmetros k
A1
, k
A2
, k
A3
.
Definimos el vector de constantes del harnero como sigue:

1
2
3
A
A A
A
k (
(
= k
(
( k

u
Segn (2.3-2) se tiene que:
I
A A S
r p = q (2.3-5)
El rechazo de agua del harnero est dado por la ecuacin (2.2-47)

1
0 1
n
r A
k
A k
q . r
=
=

Luego segn (2.3-2) se tiene que:
( ) 0 1 I
r A S
A
q . p = q
(2.3-6)
Segn [2.2-45] se tiene que:
p
A
= f
A
r
A

luego:
I
A S A S
p p p = q
( ) I I
A A S
p p = q (2.3-7)

- 53 -

Segn [2.2-48] se tiene que:

A
r
A
i
A
o
q q q =
luego:
( ) 0 1 I
o o A S
A S
q q . p = q
(2.3-8)
Chancador de pebbles
El modelo del Chancador fue desarrollado en el punto 2.2.4.
Segn la ecuacin (2.2-53) la descarga del Chancador es:
( ) ( )
1
C C
p A B A f
= I I
Segn se aprecia en el punto 2.2.4, las matrices de clasificacin A y ruptura B
dependen del flujo de alimentacin al Chancador f
C
, del F
80
y de los parmetros del
Chancador.
Definimos el vector de constantes del Chancador como sigue:

| |
1 2 3 0 12
T
C CE C C C C C
d = k k k u
Segn [2.3-3] obtenemos que:
( )( )
1
C A
p A B A Ch r
= I I
luego:
( )( )
1
I
C A S
p A B A p
= I I q Ch (2.3-9)

- 54 -

Segn [2.2-60] se tiene:

C
i
C
o
q q =
( ) 0 1 I
o r A S
A C
q q . p = = q
Ch Ch (2.3-10)

Realimentacin del ciclo de molienda SAG
Segn (2.3-4) se tiene que:
( ) 1
C A
r p r = + Ch
( ) ( ) ( )
1
I 1 I
A S A S
r A B A p p
= I I q + q Ch Ch
( ) ( ) ( )
1
I 1 I
A A S
r A B A p
(
= I I q + q

Ch Ch (2.3-11)
Por otro lado se tiene que:
( ) ( ) 1 1
r o r r r r
A A A A C
q q - q q - q q = + = + = Ch Ch Ch
( ) 0 1 I
r A S
q . p = q
( 2.3-12)

- 55 -

Ecuacin de Estado
Segn lo visto anteriormente tenemos que el vector de constantes es el vector de 35
constantes siguiente:
| |
| |
1 2 1 2
1 2 3
1 2 3 0 12
T
E
S S S S S S S S SB m B
S
T
SAG A A A A
T
C
CE C C C C C
c R k d l V w J
d
(
( o k k c

( (
( (
= = k k k
( (
( (
k k k
(

u
u u
u
(2.3-13)
Dado que:
S
SAG S
SA
w
w
(
(

X = X =
se tiene que:
( )
1 0 0
0 1 0
SAG
SAG SAG SAG SAG S SAG SAG S
r
r
d
q dt
| | | | ( (
= = +
| |
( (
|
\ . \ .
X
F X , U , F X , U , u u ( 2.3-14)
Por ltimo el vector de variables de salida queda resuelto mediante la siguiente ecuacin:
( )
( )
( )
( )
( )
( )
( ) ( )
( )
( )( )
2
1 4
1
1
I
I I
0 1 I
I
0 1 I
I
0 1 I
I 1
S
S
S
S
S o
S SA
S
A
o
A S
o A S
A
SAG SAG SAG SAG SAG
r
A S
A S
C
o
A S
r A S
r
A
S
A
S
A
C
A
c w
w
p
q
w
w
p
q
p
q . p
r
q
p
. p
p
q
A B A p
q . p
r
q
A B A
(
| | (
k
| (
k +
| (
\ .
(
(
q
(
q
(
= = =
(
q
(
(
q
(
I I q (
(
= q
(
(

I I q +
Y H X , U , u
Ch
Ch Ch
Ch C ( )
( )
I
0 1 I
A S
A S
p
. p
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
q
(
(
q

h

( 2.3-15)

- 56 -

2.3.2 Ecuacin de Estado y Entrada-Salida del circuito de Molienda Secundaria
En el circuito de molienda secundaria, la descarga de mineral del circuito de molienda
SAG entra a un pozo de bombeo. La pulpa de mineral que est en el pozo es impulsada
por medio de bombas centrfugas a hidrociclones, estos clasifican el mineral de forma que
los de tamao menor salen del proceso (bajos tamao que van al proceso de flotacin) y
los mayores son enviados a molinos de bolas. El producto de los molinos de bolas ingresa
al pozo de bombeo. Adicionalmente, a la pulpa de mineral contenida en el pozo de
bombeo, se le puede agregar agua para manejar la densidad de la pulpa. El diagrama en
bloques del circuito de molienda secundaria es mostrado en la Figura 2.3-2.

Figura 2.3-2 Diagrama de Bloques del circuito de Molienda Secundaria
Los vnculos entre los distintos flujos mostrados son los siguientes:
Pozo de Bombeo
- f
P
= f
Sec
+ p
M
( 2.3-16 )
- q
i
P

=

q
i
Sec
+ q
o
M

Bomba Centrfuga
-
P
B
t
o
P
p
f q = +
(flujo volumtrico por la bomba)

( 2.3-17 )
Pozo Hidrocicln
Molino de
Bolas
G
P
q
o
P

p
P
q
i
H

f
H
q
i
P

f
P
q
i
Sec

f
Sec
q
o
H

p
H
q
i
M

f
M
q
o
M

p
M
q
r
H
r
H
Bomba
Centrifuga
B
f
B
f
v
B

- 57 -

El flujo volumtrico por la bomba se regula mediante la velocidad de estas v
B
. Adems el
flujo volumtrico resultante depender de las presiones en los dos extremos de la bomba
y de las prdidas de energa por roce en las caeras y singularidades (codos, T, etc.). La
ecuacin [2.2-69] muestra la dependencia entre el flujo volumtrico de pulpa por la bomba
y las presiones en la alimentacin y descarga de la bomba, las que dependen de las
condiciones del pozo y los hidrociclones, y de las prdidas en caeras.
Hidrocicln
- f
H
= p
P
(2.3-18)
- q
i
H

= q
o
P

Molino de Bolas
- f
M
= r
H
(2.3-19)
- q
i
M

= q
r
H

El flujo G
P
es una variable fijada por el operador de la planta o mediante algn algoritmo
de control automtico para optimizar el proceso.
Luego se observa que el vector de variables de entrada de entrada al la molienda
secundaria es:
Sec
Sec
P
B
i
Sec
f
q
G
(
(
(
(
(
( v

U = (2.3-20)
Todas las variables de salida de los submodelos del modelo de molienda secundaria que
se ven en la Figura 2.3-2 son de inters. Tambin sern de inters, como se ver ms
adelante, el porcentaje de slidos en la pulpa de mineral c
p
y las presiones en los
extremos de la bomba centrfuga P
B1
y P
B2
, y en la entrada al Hidrocicln P
H
. Luego el
vector de variables de salida queda definido por:
1 2
(
=

Y
T
Sec P i P o B H o H r M o P B B H
P P H H M
f q p q f p q r q p q c P P P (2.3-21)

- 58 -

Determinaremos el valor de las variables de salida.
Pozo
Segn se vio en el captulo 2.2.4 el sistema dinmico del pozo est dado por el estado de
las variables nivel de la pulpa en el pozo h
p
y densidad de la pulpa en el pozo
t
p
. Las
ecuaciones (2.2-62) y (2.2-64).
p
p P P
dh
A f p
dt
=
( )( )
t t
P t i P
p p m
P P
P P
d
1
1 F 1 q G
dt A h
( | |
= + + ( |
|
(
\ .

en donde:

P
P i P
m
P
f
f q G = + +

segn (2.2-61)
P B
p f = est determinado por la bomba centrfuga. Este valor ser determinado
ms adelante
A
P
,
m
corresponden al rea transversal equivalente del pozo y densidad del
mineral. Estos valores son parmetros del modelo
Luego tenemos las siguientes ecuaciones de estado:

( )( )
1
1
1 1
P
i P B
P m
P
t
t
P t i P
p p m
P
P
P
P P
f
q G f
A
h
d
dt
f q G
A h
( ( | |
+ + ( ( |
|
( ( (
\ .
=
( (
( | | (
(
+ + ( |
(
|
(
( \ .
( 2.3-22 )

- 59 -

La ecuacin [2.2-65] entrega el valor del porcentaje de slidos en la pulpa:

1
1
t
m
p
m t
P
P
c

=

Las ecuaciones (2.2-66) y (2.2-67) entregan los valores de los flujos msicos de slidos y
agua en la salida del pozo.
P P t p
P
P p c =
1
o B t p
P P
q f ( c ) =
Luego remplazando c
p
tenemos:
( )
1
1
m
P B t
m
P
P f

=

( 2.3-23)
( )
1 1
1
m
o B t
m
P P
q f
(
=
(

( 2.3-24)
Por ltimo la ecuacin (2.2-68) entrega la ecuacin de estado para el mineral contenido
en el pozo:

P P P
d
w f p
d t
=

Tenemos que la ecuacin de estado para el pozo es:

( )( )
1
1
1 1
P P
P
P
i P B
P m
P
P
t
t
P t i P
p p m
P
P
P
P P
f p
f
w
q G f
d d
A
h
dt dt
f q G
A h
(
(
( | |
( (
+ +
( |
( ( |
( = = \ .
( (
( (
( | |

(
+ + ( |
|
(
(
\ .
X (2.3-25)

- 60 -

Bomba
La ecuacin (2.2-69) entrega el valor del flujo volumtrico de la pulpa a travs de la
bomba centrfuga
B
f , en funcin de la velocidad de giro de la bomba v
B
, de las presiones
en la entrada y salida P
B1
y P
B2
, y de constantes de la bomba k
B1
y k
B2
.

| |
B B1 B B2 B2 B1
f = k v k P P
La ecuacin (2.2-70) entrega el valor de la presin en la entrada de la bomba.

1 B t P
P
g h = P
La ecuacin (2.2-76) entrega el valor de la presin en la salida de la bomba.

2 TS T TF
B2 H t H B 2 2
T T T
P
L 1
g h f
A d 2 g A 2 g
| | k k
= + + +
|
\ .
P P
Donde P
H
es la presin en la entrada de los hidrociclones. Este valor ser determinado
ms adelante.
Son parmetros identificables del modelo y constantes los siguientes trminos:
- L
T
: Largo de la tubera entre la bomba e Hidrociclones
- d
T
: Dimetro de la tubera
- A
T
: rea transversal de la tubera
- k
TF
, k
TS
: Factor de prdidas por friccin y por singularidades
- k
B1
, k
B2
: Constantes caratersticas de la bomba

- 61 -

- g: Aceleracin de gravedad
Definimos la siguiente constante:

TS T TF
B3 2 2
T T T
L 1
A d 2 g A 2 g
k k
k = +
Luego (2.2-76) se puede escribir de la siguiente forma:
2
B2 H t H B3 B
P
g h f = + + k P P ( 2.3-26)
Hidrociclones
La ecuacin (2.2-83) entrega el valor del flujo agua en la salida del Hidrocicln:

2
1 1 3 96
o i P H H
H H
q . q ( . c )n = + k
La ecuacin (2.2-82) entrega el valor del rechazo de agua:
q
r
H
=
q
i
H
-
q
o
H

Luego se tiene:
2
0 1 3 96
r i P H H
H H
q . q ( . c ) n = + k ( 2.3-27)
La ecuacin (2.2-91) entrega el valor del flujo de mineral por intervalo de tamaos del
rechazo del Hidrocicln:

H H H
r r r = +

Utilizando (2.2-80) y (2.2-84) se obtiene:

H CH H H H
r f Y f = +

R

- 62 -

Y
H
es la matriz de clasificacin, esta matriz se determin en el captulo 2.2.6.
Utilizando (2.2-81) y (2.2-85) se obtiene:
1
i i
H H H H
r r
H H
H H
q q
r f Y f
q q
| |
| = +
|
\ .
(2.3-28)
La salida del hidrocicln est dada por la ecuacin (2.2-91):

H H H
r
~
f
~
p =
reemplazando se obtiene:
( ) ( ) 1 1
H CH H H CH H
p f Y f = R R
( ) I 1
i
H H H
r
H
H
q
p Y f
q
| |
| =
|
\ .
(2.3-29)
Por otro lado la presin en la entrada a los hidrociclones est dada por la ecuacin
(2.2-89)

( )
2
0 125
1
13 725
1 48 1
H
H .
vf H p H
f
.
. d n c
| |
|
=
|
k
\ .
P
Se define la siguiente constante
( )
4 2
1
13 725
1 48
H
vf H H
.
. d n
k =
k
(2.3-30)
luego se tiene.
( )
2
4
0 25
1
H H
H .
p
f
c
k
=
P (2.3-31)

- 63 -

De la Figura 2.3-2 se observa que:

H B
f f =
luego:
( )
2
4
0 25
1
H B
H .
p
f
c
k
=
P ( 2.3-32)
Reemplazando (2.3-32) en (2.3-26) se tiene:
( )
2
H4 B 2
B2 t H B3 B 0.25
p
P
f
g h f
1 c
k
= + + k
P

( )
2 H4
B2 B3 B t H 0.25
p
P
f g h
1 c
| |
k
|
= + k +
|
\ .
P
Reemplazando en (2.2-69) se tiene que:

( )
2 H4
B B1 B B2 B3 B t H t p 0.25
p
P P
f f g h g h
1 c
( | |
k
( |
= k v k + k +
( |
\ .

Definimos la siguiente funcin g, la que es una funcin de c
p
.
( )
( )
H4
p B2 B3 0.25
p
c
1 c
| |
k
|
= k + k
|
\ .
g (2.3-33)
luego:

( ) ( )
2
B p B B1 B B2 t H P
P
f c f g h h = + k v k g

- 64 -

resolviendo la ecuacin de segundo grado se obtiene:
( ) ( )
( )
p B1 B B2 t H P
B
p
P
1 1 4 c g h h
f
2 c
(
+ + k v k

=
g
g
(2.3-34)
Molino de Bolas
En el captulo 2.2.7 se desarrollo el modelo dinmico del molino de bolas. Segn este
modelo el hold-up de mineral est determinado por la ecuacin (2.2-94):

M M M M
d
w f p Bw
dt
= +
La descarga de mineral est dada por la ecuacin (2.2-96)

M M M
p w = o
De la Figura 2.3-2 se deduce que:

M H
f r =
Luego tenemos que la ecuacin de estado del molino se puede escribir como sigue:

M H M M M
d
w r w B w
dt
= o +
( ) I
M H M M
d
w r B w
dt
= o ( 2.3-35)
Por otro lado las ecuaciones (2.2-97) y (2.2-98) describen el proceso del agua en el
molino. Estas ecuaciones son las siguientes:

MA i o
M M
d
w q q
dt
=
MA MA o
M
q w = o

- 65 -

Y dado que:

i r
M H
q q =
Se tiene que:

MA MA MA r
H
d
w q w
dt
= o ( 2.3-36 )
Luego los vectores de variables de estado, variables de entrada, variables de salidas y
constantes del Molino de Bolas son:

M
M
MA
w
w
(
(

X = ( 2.3-37 )

H
M
r
H
r
q
(
(
(

U = ( 2.3-38 )

M
M
o
M
p
q
(
(
(

Y = ( 2.3-39 )

M M
MA
B (
(
o =
(
(
o

u ( 2.3-40 )
Las ecuaciones de estado y entrada-salida son:

( )
( ) I
H M M M
M
M M M M
MA MA MA r
H
r B w w
d d
q w w dt dt
( o
(
= =
(
(
o
(

X
= F X , U , u
( 2.3-41 )
( )
M
M M
M M M M
o
MA MA
M
p
w
q
w
( o (
= = =
(
(
o
(

Y H X , u ( 2.3-42 )

- 66 -

Molienda Secundaria
De lo visto anteriormente se obtiene el vector de variables de estado:

P
P
P
t
Sec
P
M
M
MA
w
h
w
w
(
(
(
(
(
=
(
(

(
(

X
X =
X
( 2.3-43 )
El vector de variables de entrada es (2.3-20)

Sec
Sec
P
B
i
Sec
f
q
G
(
(
(
(
(
( v

U =
El vector de variables de salida es (2.3-21)

1 2
(
=

Y
T
Sec P i P o B H o H r M o P B B H
P P H H M
f q p q f p q r q p q c P P P
De lo visto anteriormente se deduce que el vector de constantes es:
|
|
1 2
3 1 2 3 4
Sec M P T T T TF TS B B
T
B H H H H H H vf M MA
A g L d A
n Y d B
= k k k k
k k k k k o o
u
(2.3-44 )
De la Figura 2.3-2 se deduce que:

P Sec M
f f p = +

i i o
P M Sec
q q q = +

- 67 -

Luego las ecuaciones de estado y entrada-salida de la Molienda Secundaria son las
siguientes:
( )
( )( )
( )
1
1
1 1
I
P P
P
P
i P B
P m
P
Sec
t
t Sec Sec Sec Sec
P
P t i P
p p m M
MA
H M M
MA MA
P
P
P P
r
H
f p
f
w
q G f
A
h
d d
dt dt f q G
A h w
w
r B w
q w
(
(
( | |
(
(
+ +
( |
( |
(

(
\ .
(
(
(
( ( = = | |
(
(
+ + ( |
|
(
(
(
\ .
(
(
( o
(
o
(

X
= F X , U , u
( 2.3-45 )
La ecuacin matricial de entrada-salida de la Molienda Secundaria es:

- 68 -

( )
( )
( )
( ) ( )
( )
1 2
1
2
1
1
1 1
1
1 1 4
2
I
Sec M M
i MA MA
m
B t
m
m
B t
m
p B B B t H P
p
P
i
P
o
B
H
o
Sec Sec Sec Sec
H
r
M
o
P
B
B
H
Sec
P
P
P
P
P
H
H
M
f w
q w
f
f
c g h h
c
f
q
p
Y
q
f
p
q
r
q
p
q
c
+ o
+ o

(

(

(
+ + k v k

(
(
(
(
(
(
(
(
(
(
(
(
= = =
(
(
(
(
(
(
(
(
(
(
(
(

Y H X ,
g
g
P
P
P
u
( )
( )
( )
( ) ( )
( )
2
2
1 1
1
1
1 1 1 1 3 96
1
1 1 1 1
1 1
1
0 1 1 1 3 96
1
m
B t
m
H P
r
m
B t P H H
m
m m
B t B t
m m
P H P
r r
m
B t P H H
m
M M
P
H
P
P P
H H
P
f
p
q
. f ( . c )n
f f
p Y p
q q
. f ( . c )n
w
| | (

|
(

|
|
|
|
\ .
(
+ k
(

| | ( (

|
( (

|
+
|
|
|
\ .
(
+ k
(

o
o
( )
( )
2
4 2
3 0 25
2
4
0 25
1
1
1
1
MA MA
t
m
m t
t P
H B
t H B B .
p
H B
.
p
P
P
P
P
w
gh
f
g h f
c
f
c
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(

(
(

(
(
(
( k
+ + k
(
(
(
k (
(
(

( 2.3-46 )

- 69 -

2.3.3 Modelo Global
El modelo global sale de interconectar los dos circuitos de molienda, la interconexin de
flujos de slidos y agua es:
- f
Sec
= p
A
(2.3-47)
- q
i
Sec

= q
o
A

La planta, vista como un slo bloque, recibe mineral grueso del chancado primario y
produce una pulpa de mineral que ir al proceso de flotacin.

Figura 2.3-3 Diagrama de la Planta como un Bloque
Las ecuaciones de estado y entrada-salida, as como los vectores de estados, entrada,
salida y constantes estn dados por la suma de las ecuaciones y variables de la molienda
SAG y la molienda Secundaria. Luego:

SAG
Planta
Sec
(
(

X
X =
X
(2.3-48)

SAG
Planta
Sec
(
(

U
U =
U
(2.3-49)

SAG
Planta
Sec
(
(

Y
Y =
Y
(2.3-50)
Planta de Molienda SAG
y Molienda Secundaria
q
o
H

q
i

p
H
f

- 70 -

SAG
Planta
Sec
(
=
(

u
u
u
(2.3-51)
( )
( )
( )
( (
= =
( (

F X , U , X X
= F X , U ,
F X , U , X
SAG SAG SAG SAG SAG Planta
Planta Planta Planta Planta
Sec Sec Sec Sec Sec
d d
dt dt
u
u
u
(2.3-52)
( )
( )
( )
( (
= = =
( (

Y H X ,
Y H X ,
Y H X ,
SAG SAG SAG SAG
Sec Planta Planta Planta
Sec Sec Sec Sec
u
u
u
(2.3-53)
Para implementar el modelo de la planta en el simulador se utiliz el modelo equivalente
en tiempo discreto (ver 2.4). El detalle de la implementacin se presenta en 3.2.

- 71 -

2.4 MODELO DE TIEMPO DISCRETO DE LA PLANTA
Todos los modelos descritos en los captulos anteriores son modelos matemticos
generales y luego aplicables en tiempo continuo. Sin embargo para fines de programacin
del simulador de la Planta de Molienda, debemos utilizar modelos en tiempo discreto.
Para la determinacin de las ecuaciones en tiempo discreto de los modelos se utiliz una
discretizacin lineal tipo Euler, es decir si se tiene el sistema lineal siguiente:
( )
d
dt
=
X
F X, U, u
( ) = Y H X, U, u
El sistema discretizado es:
( )
1 n n n n n
t
+
= A X X + F X , U , u
( )
n n n n
= Y H X , U , u
( )
0
0 t = = X X
Describiremos a continuacin los modelos en tiempo discreto de cada uno de los
submodelos. El modelo en tiempo discreto de la planta se obtiene interconectando los
modelos de los procesos individuales de la misma forma que se interconectan los
modelos en tiempo continuo (ver 2.3)

- 72 -

2.4.1 Molino SAG
En 2.2.1 se describi el modelo de un molino SAG, segn este el proceso se puede dividir
en dos: Molienda y Clasificacin (ver Figura 2.2-2).
La etapa de molienda se describe mediante la ecuacin de estado del hold-up de mineral
dentro del molino. La ecuacin (2.2-12) es:

S S S S
d
w (t) f (t) p (t) k(t)w (t)
dt
=

En la representacin en tiempo discreto se tiene que la variacin del hold-up, para
cualquier tamao de mineral, en el instante de tiempo n, es igual al mineral que ingresa a
la cmara de molienda menos el mineral que sali de la cmara de molienda y menos el
mineral que cambi de tamao, es decir:

S S
S
S S S
n 1
n
n 1
n
n 1 n 1
w w
w
f p k w
t t

A
= =
A A

(2.4-1)
en que:

S S n n
w w (t t ) = =

n 1
n 1
k k(t t )
= =

S S n n
f f (t t ) = =

S S n 1
n 1
p p (t t )
= =

El mineral que ingresa a la cmara de molienda es la suma del mineral fresco que ingresa
al molino en el instante n ms el mineral que no se descargo en el instante n-1, es decir:

S S S
n n
n 1
f f r

= +
(2.4-2)

- 73 -

Despejando
Sn
w en [2.4-1] obtenemos el hold-up en n:

( ) S S S S S
n 1 n n
n 1 n 1 n 1
w w f p k w t

= + A

(2.4-3)
En el instante n, el mineral que ingresa a la fase de clasificacin, segn (2.2-28) y (2.2-
29), est dado por:

S
S S
S
n
n
n
w
p
W
= o
(2.4-4)
en que:

S Si
i
n n
W w =

Luego la descarga del molino SAG y el mineral rechazado por la parrilla interna del molino
son:
( ) I
S S n n n
p c p =

(2.4-5)

S S n n n
r c p =

(2.4-6)
De forma anloga tenemos que las ecuaciones para el agua (2.2-30) y (2.2-31) son:

S
o
S
i SA
q q w
dt
d
=

o S SA
S
q w = o
En tiempo discreto se expresan como:

SA SA
SA
i o
n
n 1
Sn Sn 1
w w
w
q q
t t

A
= =
A A
(2.4-7)

( ) SA SA i o
n
n 1 Sn Sn 1
w w q q t

= + A (2.4-8)

o S SAn
Sn
q w = o (2.4-9)

- 74 -

2.4.2 Harnero
En el captulo 2.2.2 se describi el modelo del Harnero. Este modelo considera una
ecuacin de eficiencia para la determinacin de los flujos de mineral que se descargan y
los que son rechazados, la ecuacin de eficiencia para el mineral de tamao X
i
est dada
por la ecuacin (2.2-33):
A2
i
A A1
50
i
X
(t) 1 exp
d (t)
k
(
| |
( q = k
|
(
\ .

En tiempo discreto esta ecuacin se expresa como:

A2
i
A n A1
50n
i
X
1 exp
d
k
(
| |
(
q = k
|
|
(
\ .

(2.4-10)
El dimetro para el cual pasa el 50% de las partculas se calcula con la siguiente
ecuacin:

0.148
An ua
50n
A3
F F
d
1.6
| |
=
|
k
\ .
(2.4-11)
En que:
=
i
n n A
f F
i
A

Luego la ecuacin de rechazo de mineral es:

A n A n A n
i i i
r f i = q ( 2.4-12 )
Luego, la descarga del harnero ser:
p
A n
= f
A n
r
A n
( 2.4-13 )

- 75 -

Segn se vio en 2.2.2 el rechazo de agua en el harnero corresponde a la humedad del
mineral, la que se estima en un 10%. La ecuacin en tiempo discreto para el rechazo de
agua es:

r A n
i
An i
q 0.1 r =

( 2.4-14 )
y por lo tanto, la salida de agua es:

n A n A n A
r i o
q q q = ( 2.4-15 )
2.4.3 Chancador de pebbles
En el captulo 2.2.3 se desarroll el modelo del Chancador de pebbles utilizado en el
simulador.
El modelo considera que el mineral que ingresa al Chancador primero pasa por una fase
de clasificacin representada por la matriz A, en que las partculas de menor tamao
pasan directamente a la descarga, y luego el flujo de partculas rechazadas pasa a una
fase de molienda representada por la matriz B.
La descarga del chancador est dada por la ecuacin (2.2-44)
( ) ( )
1
C C
p (t) A(t) B(t)A(t) f (t)
= I I
Escribiendo esta ecuacin en tiempo discreto tenemos:
( ) ( )
1
Cn n n n Cn
p A B A f
= I I ( 2.4-16 )
Las matrices A y B dependen de parmetros caractersticos del Chancador a identificar, y
del flujo total de entrada al Chancador F
C
(sumatoria de los flujos de entrada al Chancador
por intervalo de tamaos).

- 76 -

El flujo de agua en la alimentacin al Chancador es igual a la descarga en todo momento,
luego:

Cn Cn
=
o i
q q ( 2.4-17 )
2.4.4 Pozo
En el captulo 2.2.4 se desarrollo el modelo matemtico del pozo utilizado en el simulador,
segn este el flujo de salida de este proceso, tanto de slidos como de agua, dependen
de la ecuacin de estado para la densidad de mineral en el pozo. Esta ecuacin es (2.2-
55):

( )( )
t t
P t i P
p p m
P P
P P
d (t) (t)
1
1 F (t) 1 (t) q (t) G (t)
dt A h (t)
( | |
= + + ( |
|
(
\ .

En tiempo discreto diremos que la densidad de mineral en el intervalo n depender del
nivel del pozo en el intervalo n, los flujos de slidos y agua que ingresen al pozo en n y la
densidad de mineral en el pozo que se tena en n-1, es decir:
( )( )
t t t
Pn t i P
n
n
p p m
n
Pn Pn-1 Pn-1
Pn-1 P
1
1 F 1 q G
t A h
( | |
( | = + +
|
A
(
\ .
( 2.4-18 )
Luego la densidad del mineral en el pozo en el intervalo n es:
( )( )
t
t t Pn t i P
n
n
p p m
n
Pn-1
Pn Pn-1 Pn-1 P
1
1 F 1 q G t
A h
( | |

( | = + + + A
`
|
(

\ . )
( 2.4-19 )
El nivel del pozo en t est dado por la frmula [2.2-53], vista en 2.2.4, esta es:

p
p P P
dh (t)
A f (t) p (t)
dt
=

- 77 -

Como se ve el nivel del pozo depende de los flujos volumtricos de entrada y salida del
pozo. En tiempo discreto decimos que la diferencia del nivel del pozo entre el intervalo n y
n-1 depende del flujo volumtrico que entra en n y el flujo volumtrico que sali en n-1, es
decir:

( )
p p
n n 1
P P
n n 1
p
h h
1
f p
t A
=
A

( ) p p P P
n n 1 n n 1
p
1
h h f p t
A

(
= + A
(
(

( 2.4-20 )
El flujo volumtrico de salida lo determina el estado de operacin de la bomba.
Luego los flujos de slidos y agua en la salida del pozo, en el intervalo n, son:

P P t p
n n n Pn
p p c = ( 2.4-21 )

o P t p
n n Pn Pn
q p (1 c ) = ( 2.4-22 )
En donde el porcentaje de slidos en la pulpa en todo instante esta dado por la ecuacin
[2.2-56], es decir:

1
1
t
m
p
m t
P
P
c

=

El hold-up de mineral dentro del pozo est dado por la discretizacin de la ecuacin (2.2-
59), esto es:

( )
1 1
b b b b
n n n n
w w f p t

= + A

- 78 -

2.4.5 Bomba Centrfuga
En el captulo 2.2.5 se describi el modelo de bomba centrfuga utilizada en el simulador.
La ecuacin que rige el comportamiento es (2.2-59), aqu se determina el flujo volumtrico
de pulpa que mueve la bomba en funcin de su velocidad y de las presiones en los dos
extremos de la bomba.

| |
B B1 B B2 B2 B1
f (t) (t) (t) (t) = k v k P P
Esta ecuacin escrita en tiempo discreto queda como sigue:

B B1 B B2 B2 B1
n n n n
f
(
= k v k
(

P P ( 2.4-23 )
Las presiones en los extremos de la bomba quedan determinadas por (2.2-60) y (2.2-68).
Las ecuaciones equivalentes en tiempo discreto son las siguientes:

B1 t P
n n
p
n
g h = P ( 2.4-24 )

2
B2 H t H B3 B
n n 1 n 1
p
n
g h f

= + + k P P ( 2.4-25 )
2.4.6 Hidrocicln
En el Captulo 2.2.6 se desarroll el modelo de Hidrocicln utilizado en el simulador,
resumiendo las ecuaciones se obtiene que:

Hidrocicln
Y
H q
i
H

f
H
q
o
H

p
H
q
r
H

r
H

- 79 -

| |
( )
( ) I ( ) 1 ( )
( )
i
H H H
o
H
H
q t
p t Y t f t
q t
(
= (
(

( )
( ) ( ) ( ) ( )
( )
i
H H H H
o
H
H
q t
r t f t Y t f t
q t
= +
( ) ( ) ( ) =
r i o
H H H
q t q t q t
( )
2
( ) 1.1 ( ) 3.96 ( )
o i P H H
H H
q t q t c t k n = +
En donde:
- Y
H
(t) es la matriz de clasificacin la que slo depende q
i
H
(t) y de
( ) ( ) =
H H
F t f t

segn se vio en 2.2.6.
- k
H
es una constante identificable para cada hidrocicln y n
H
es el nmero de
hidrociclones cuando se tiene una batera de hidrociclones.
- c
P
(t) es la densidad de los slidos en la alimentacin.
Como se aprecia todas ecuaciones mencionadas slo dependen de las variables de
entrada. Las ecuaciones equivalentes en tiempo discretos son:
I 1
i
H H H
o
H
n
n n n
H
n
q
p Y f
q
(
(
( =

(

( 2.4-26 )
i
H H H H
o
H
n
n n n n
H
n
q
r f Y f
q
= + ( 2.4-27 )
r i o
H H H
n n n
q q q = ( 2.4-28 )
( ) 2
1.1 3.96
o i P H H
n
H H
n n
q q c k n = + ( 2.4-29 )

- 80 -

2.4.7 Molino de Bolas
En el Captulo 2.2.7 se desarroll el modelo del Molino de Bolas que se utiliza en el
simulador. El modelo considera, tanto para los slidos como para el agua, una ecuacin
de estado para el hold-up y una descarga proporcional al hold-up. Las ecuaciones
desarrolladas en 2.2.7 se resumen a continuacin:
M M M M
d
w (t) f (t) p (t) B w (t)
dt
= +
M M M
p (t) w (t) = o
MA i o
M M
d
w (t) q (t) q (t)
dt
=
MA MA o
M
q (t) w (t) = o
B es una matriz triangular inferior con las tazas de molienda para cada intervalo de
tamao de los slidos. o
M
es la tasa de descarga de slidos y es la misma para todos los
intervalos de tamao. o
MA
es la tasa de descarga de agua.
Las ecuaciones equivalentes en tiempo discreto para los slidos son las siguientes:
1
1 1
M M
M
M M M
n n
n n n
w w
w
f p B w
t t

A
= = +
A A
(2.4-30)
( )
1 1 1
M M M M M
n n n n n
w w f p B w t

= + + A (2.4-31)
M M M
n n
p w = o (2.4-32)

- 81 -

Las ecuaciones equivalentes en tiempo discreto para el agua se indican a continuacin:
1
1
MA MA
MA
i o
n n
M M
n n
w w
w
q q
t t

A
= =
A A
(2.4-33)
1
1
MA MA i o
n n M M
n n
w w q q t

| |
= + A
|
\ .
(2.4-34)
o MA MA
n
M
n
q w = o (2.4-35)

- 82 -

2.5 CONTROL DE LA PLANTA
2.5.1 Perturbaciones (medidas y no medidas)
El mineral que llega a la planta de molienda no es siempre el mismo, vara por diversas
razones tanto en dureza (perturbacin no medida) como en la granulometra que tiene
(perturbacin medida). Si se asume que el proceso minero es el mismo, es decir las
tronaduras se realizarn de la misma forma (por ejemplo misma cantidad de explosivo), el
cmo carga la pala los camiones, cmo estos transportan y descargan el mineral en el/los
chancador(es) primario(s), y el cmo el mineral chancado cae a las correas
transportadoras que llegan al molino SAG es realizado siempre de la misma forma,
dureza y granulometra son variables no independientes. Mineral ms duro implicar una
granulometra con ms gruesos.
Mayor dureza afecta fuertemente al molino SAG ya que la eficiencia de las rocas para
molerse en los choques entre ellas disminuye y ms pebbles son generados aumentando
la recirculacin por el chancador de pebbles.
En los molinos de bolas la dureza tambin influye, pero debido a que la molienda es va
choque contra bolas de acero el efecto es menor.
i) Perturbaciones en la Molienda SAG
Tanto que el mineral sea ms duro como que su granulometra sea mayor se puede
interpretar, desde el punto de vista de la modelacin, como una prdida en la eficiencia
para moler del Molino SAG. De acuerdo a (2.2-11) tenemos que:
=

S S S S
d
w f p k w
dt

Con:
S S S
r f f
~
+ =
S S
r c p =

( ) I
S S
p c p =

- 83 -

La eficiencia del Molino se observa (se puede medir) en la descarga p
S
, esta depende de
cunto mineral r
S
recircule dentro del molino a causa de la parrilla interna (modelada
mediante la constante de descarga c) y de la eficiencia de la molienda dada por la matriz
k. La matriz k, como se vio en 2.2.1, se determin de la siguiente manera:
A 1
k R k R
=
A E
i i
S
S
k k
W
=
As tenemos que la eficiencia de la molienda depende de las constantes especficas de
molienda k
E
.
Se determinan tasas especficas de molienda para un mineral de dureza media
E
k (ver
captulo 3.3.3). Cuando el mineral sea el ms duro posible, se asumir que las tasas
especficas de molienda son un 30% menos eficientes que las tasas para mineral de
dureza media y que cuando el mineral es lo ms blando posible, las tasas especficas de
molienda sern un 30% ms efectivas.
Sea la proporcin de la dureza del mineral, con e [0,1]. = 0 es el mineral ms blando
posible, = 1 es el mineral lo ms duro posible y = 0.5 representa al mineral de dureza
media, luego k
E
es en cada iteracin:
( )
( )
( )
= +
E E E
k 0.7k 1 1.3k ( 2.5-1 )
Para simular los cambios en la granulometra, se definen tres granulometras posibles
para la alimentacin fresca de mineral al Molino SAG, una granulometra para mineral
duro f
D
, uno para el mineral de dureza media f
N
y uno para el mineral blando f
B
. De esta
forma la granulometra de alimentacin fresca f ser:

- 84 -

< s
= < s
< s
D
N
B
si 0.8 1
si 0.3 0.8
si 0 0.3
f
f f
f
( 2.5-2 )
Y el flujo de alimentacin:
S S
f F = f ( 2.5-3 )
La proporcin de dureza se define de la siguiente forma:
( ) n = +
n1
1
1
2
( 2.5-4 )
En que n es un valor aleatorio y
n1
es una constante que determina cuanto afecta n a la
dureza del mineral fresco.
La variable n da una variacin aleatoria (ruido). Sin embargo n no se puede definir como
ruido blanco ya que la variacin en el tiempo de la dureza no es tan rpida como la que
obtendramos al ocupar directamente este tipo de ruido. Para crear un ruido de menor
frecuencia definimos al ruido n como la suma de un ruido blanco n
B
ms una fraccin de
la variacin en el tiempo que tiene el ruido n. Con esta definicin obtenemos una inercia
del ruido con sus valores pasados. Se define n como sigue:

B n2
d
dt
= + n n n ( 2.5-5 )
La constante
n2
permite regular la dependencia que tiene n con sus valores pasados.
n
B
se define entre [-1,1], y
n2
entre [0,1), luego n estar tambin en el rango [-1,1].
En la molienda secundaria (como se ver ms adelante) se utilizan dos flujos de slidos,
uno de mineral duro y otro de mineral blando. La relacin entre la dureza del mineral en el

- 85 -

SAG dada por y los dos flujos en la molienda secundaria es realizada en la descarga del
harnero (ver 2.5.2 iii). Si bien el rechazo del harnero, ya sea en forma directa o va el
chancador de pebbles, se suma a la alimentacin fresca en la entrada al SAG, se
considera que no afecta a . Esto es vlido ya que el rechazo del harnero es mucho
menor que la alimentacin fresca (aproximadamente un 10%), de no serlo se estara muy
alejado del punto de operacin en equilibrio.
ii) Perturbaciones en el molino de bolas
Para representar la dureza del mineral en los molinos de bolas, se defini dos matrices de
molienda. Una de estas matrices contiene tasas de molienda que representan la eficiencia
del molino para moler mineral duro (B
D
) y la otra matriz est conformada por tasas de
molienda para la conminucin de mineral blando (B
B
).
Como se mencion anteriormente, en la molienda secundaria se considera dos flujos de
mineral, uno de mineral duro y otro de mineral blando en cada modelo. En el molino de
bolas el flujo de mineral duro se representa por
D
M
f y el de mineral blando por
B
M
f .
La porcin del flujo de alimentacin conformada por mineral duro es tratada por la matriz
de molienda para mineral duro y la conformada por mineral blando por la matriz para
mineral blando. Si
D
M
w es el hold-up de mineral duro y
B
M
w el hold-up de mineral blando,
entonces:
D D D D D
M M M M
d
w f p B w
dt
= + (2.5-6)
B B B B B
M M M M
d
w f p B w
dt
= + (2.5-7)
Como la granulometra de mineral duro y blando est representada en los flujos de
entrada y hold-ups, as como la conminucin es diferenciada con diferentes matrices de
molienda; se consider una sola tasa de descarga (o
M
) para ambos flujos.
D D
M M M
p w = o (2.5-8)
B B
M M M
p w = o (2.5-9)

- 86 -

iii) Perturbaciones en el resto de los modelos
Los modelos del harnero y del hidrocicln son basados en clasificacin y no consideran
ninguna conminucin, por lo tanto lo que importa es la distribucin granulomtrica y no la
dureza.
Sin embargo, como la descarga del harnero es la alimentacin fresca del circuito de
molienda secundaria, es necesario definir los flujos de slidos duros y blandos en funcin
del parmetro de dureza. Estos flujos se definen como sigue:
=
D
A A
p p ( 2.5-10 )
( ) =
B
A A
p 1 p ( 2.5-11 )
Para el chancador de pebbles se asume que la eficiencia de conminucin no vara con los
cambios en la dureza del mineral.
En el pozo se mantiene en forma independiente el hold-up de duros y de blandos en
funcin de los flujos de duros y blandos que recibe. El control de nivel del pozo se hace
con la suma de los hold-ups.
El modelo de la bomba centrfuga se basa en flujos volumtricos y presiones. Los flujos
volumtricos y presiones son determinados por la suma de los flujos volumtricos de
mineral duro y blando. La distribucin granulomtrica se considera que no vara entre la
entrada y la salida.
En las Figura 2.5-1 y Figura 2.5-2, se presentan los diagramas de bloques que indican
todos los flujos considerando la dureza del mineral cuando corresponde

- 87 -

Figura 2.5-1 Diagrama de Bloques de perturbaciones en Molienda SAG

Figura 2.5-2 Diagrama de Bloques de perturbaciones en Molienda secundaria
SAG
k
E
(), c
Harnero
q
A
Chancador
A, B
G
q
o
S

p
S
q
i
A

f
A
q
i
S
+
q
r

f
S
=f+r f=F f()

q
o
A

p
A
D
= p
A

p
A
B
= (1-)p
A
q
i
C

f
C
q
o
C

p
C
q
r
A
r
A q
r
r

Pozo
w
p
D
, w
p
B

Hidrocicln Molino de
Bolas
G
P
q
i
H

=
+
=
+

H B t p
D B
P P
H H D B
P P
F f c
w w
f F
w w

q
i
P

f
P
q
i
Sec
=
q
o
A

f
D
Sec
= p
A
D

f
B
Sec
= p
A
B

q
o
H

p
H
q
i
M

f
M
D
= r
M
D

f
M
B
= r
M
B

q
o
M

p
M
D

p
M
B

q
r
H

Bomba
k
B1
, k
B2

B
f
v
B
v k k
B B B1 B1
f ( , , , ) P
=
+
=
+

D
D P
H H D B
P P
B
B P
H H D B
P P
w
r r
w w
w
r r
w w

- 88 -

2.5.2 Controles automticos en la planta
i) Lazos de Control
El simulador cuenta con cuatro controles automticos:
- Control PI de hold-up del Molino SAG el que acta sobre la alimentacin fresca al
molino.
- Control PID de nivel del pozo que acta sobre la velocidad de la bomba
centrfuga del pozo.
- Control PID de porcentaje de slidos en la alimentacin a hidrociclones que acta
sobre el agua que se agrega al pozo.
- Control PID de +65# en el producto final que acta en cascada con el control
anterior variando el set-point del control de Cph.
Los lazos de control se resumen en el esquema de la Figura 2.5-3.
LIC
AIC
LT
AT
Lp*
Cph*
+65#
*
Velocidad Bomba
DT DIC
Set-point de Cph
Agua al pozo
WT WIC
Hold-up*
Alimentacin Fresca
LIC LIC
AIC AIC
LT LT
AT AT
Lp*
Cph*
+65#
*
Velocidad Bomba
DT DT DIC DIC
Set-point de Cph
Agua al pozo
WT WT WIC WIC
Hold-up* Hold-up*
Alimentacin Fresca

Figura 2.5-3 Lazos de control

- 89 -

ii) Control PI de hold-up del Molino SAG
En el caso del SAG el control implementado es la versin discreta incremental del
controlador PI. El hold-up del molino es controlado variando la alimentacin fresca al
molino.
1
1
1 1 1
S S C
i
e (t) e (t 1) e (t)
F (t) F (t 1) K
t T
| |

= + + |
|
A
\ .
( 2.5-12 )
Donde F
S
(t) es el flujo de slidos en la alimentacin en [t/h],
1
C
K es la constante
proporcional,
1
i
T el tiempo integral y e
1
(t) el error entre el hold-up medido y el set-point

1 S S
e (t) w (t) w (t)
-
= ( 2.5-13 )
con
S
w (t)
-
es el set-point
iii) Control PID del nivel del Pozo
El control PID del nivel del pozo (h
p
) se hace controlando la velocidad de las bombas
centrfugas.
2 2
2
2 2 2 2 2 2
B B C d 2
i
e (t) e (t 1) e (t) e (t) 2e (t 1) e (t 2)
(t) (t 1) K T
t T t
| |
+
v = v + + + |
|
A A
\ .
( 2.5-14 )
en que:
2 P P
e (t) h (t) h (t)
-
= ( 2.5-15 )
con
P
h (t)
-
es el set-point
v
B
: velocidad de la bomba centrfuga en [RMS]
2 2 2
C i d
K ,T ,T : Constante proporcional, tiempo derivativo y tiempo integral

- 90 -

iv) Control PID del porcentaje de slidos en la alimentacin a hidrociclones
El control PID del porcentaje de slidos en alimentacin de los hidrociclones (c
p
) se realiza
controlando la alimentacin de agua fresca al pozo:
3 3
3
3 3 3 3 3 3
P P C d 2
i
e (t) e (t 1) e (t) e (t) 2e (t 1) e (t 2)
G (t) G (t 1) K T
t T t
| |
+
= + + + |
|
A A
\ .
( 2.5-16 )
en que:
3 P P
e (t) c (t) c (t)
-
= ( 2.5-17 )
con
P
c (t)
-
es el set-point
G
P
: Flujo msico de agua al pozo en [T/h]
3 3 3
C i d
v) Control PID del +65# en el Producto Final
El control PID del +65# en el producto final de los hidrociclones se realiza controlando el
set-point del control PID del porcentaje de slidos en la alimentacin de los hidrociclones:
4 4
4
4 4 4 4 4 4
P P C d 2
i
e (t) e (t 1) e (t) e (t) 2e (t 1) e (t 2)
c (t) c (t 1) K T
t T t
- -
| |
+
= + + + |
|
A A
\ .
( 2.5-18 )
en que:
| |
4
e (t) 65# (t) 65#(t)
-
( = + +

( 2.5-19 )
Con 65# (t)
-
+ es el set-point
P
c
-
: Set-point del control PID del porcentaje de slidos (c
P
) en la alimentacin a los
hidrociclones.
4 4 4
C i d

- 91 -

3 DESARROLLO DEL SIMULADOR
3.1 ESPECIFICACIN DEL SIMULADOR
El simulador desarrollado representa la operacin de una planta de molienda
semiautgena (SAG) con un circuito inverso de molienda secundaria.
El simulador cuenta con un entorno grfico que es similar a las pantallas de control que se
utilizan en las plantas reales. En el punto 3.2.6 se detalla el ambiente grfico del
simulador.
Para programar el simulador se utiliz el software MatLab versin 4.1 debido a ser este un
lenguaje de programacin matemticamente muy poderoso.
3.2 PROGRAMACIN DEL SIMULADOR
3.2.1 Lgica constructiva
El simulador est constituido por 6 sub-programas o mdulos que interactan entre s.
Estos son:
- DSAG: Programa maestro que controla el resto de los subprogramas
- DSAG_ini: Programa que inicializa el entorno grfico y crea las constantes y
variables de cada modelo con sus valores iniciales
- SAG_mod: Programa que simula la molienda SAG, es decir el funcionamiento del
Molino SAG, del harnero y del chancador de pebbles

- 92 -

- PHc_mod: Programa que simula el funcionamiento del pozo, las bombas
centrfugas y las bateras de hidrociclones
- MBol_mod: Programa que simula el funcionamiento del Molino de Bolas.
- DSAG_set: Programa para cambiar los valores de las constantes de cada
modelo y de los controles automticos
Los controles automticos que tiene implementados el simulador estn programados
dentro de los subprogramas que contienen los modelos con las variables a controlar. Es
decir, el control del hold-up del molino SAG est en el mdulo SAG_mod y los controles
del nivel del pozo, porcentaje de slidos en la alimentacin de los hidrociclones y el +65#
en el rebose de los hidrociclones estn en el mdulo PHc_mod.
3.2.2 Mdulo DSAG
La Figura 3.2-1 muestra el diagrama de flujo del mdulo DSAG, el que adems
corresponde al diagrama general del simulador. A continuacin se explica el diagrama de
flujo.
El simulador comienza llamando al subprograma DSAG_ini, el que inicializa el entorno
grfico (ver captulo 3.2.6) y crea las variables y constantes que ocupa el simulador.
Inicializadas todas las variables y constantes, el programa queda en Pausa hasta que se
da la orden de Comenzar con un botn que est en el entorno grfico. Dada la orden de
Comenzar se comienza el ciclo de simulacin de la planta.
Para que el usuario del programa pueda interactuar con el simulador durante la
simulacin, cada vez que se termina de correr un sub-programa, DSAG ve si existe una
orden de interrumpir el programa. Las rdenes de interrupcin son causadas cuando el
usuario presiona algn botn, cambia el valor de una casilla de datos o mueve una barra
de desplazamiento.
El funcionamiento del simulador es como sigue:

- 93 -

- Se llama al subprograma SAG_mod en el que se simula el funcionamiento del
molino SAG, el harnero y el chancador de pebbles.
- Se retorna a DSAG, se ve si existe alguna interrupcin y si no hay se llama al
subprograma PHc_mod donde se simula el funcionamiento del pozo, las bombas
centrfugas y los hidrociclones.
- Se retorna a DSAG, se ve si existe alguna interrupcin y si no hay se llama al
subprograma MBol_mod donde se simula el funcionamiento de los molinos de
bolas.
- Si existi alguna interrupcin sucede lo siguiente: si la interrupcin no fue
Pausa, sino que alguna de las otras que se pueden ver en el diagrama de flujo
(procesos 1, 2, 3 y 4), se desarrolla la accin solicitada y se vuelve al proceso de
simulacin. Si la interrupcin fue Pausa, se detiene la simulacin y se espera
una orden del usuario. Las rdenes pueden ser alguno de los procesos 1, 2, 3 y 4
en cuyo caso se desarrolla el proceso y se vuelve al estado Pausa hasta que se
d la orden Comenzar.
Cuando se ordena la accin 4 se corre el programa DSAG_set, donde se puede cambiar
el valor de constantes de los modelos y de los controles automticos.

- 94 -

Figura 3.2-1 Diagrama de flujo general del Simulador

D
S
A
G
_
i
n
i

I
n
i
c
i
a
l
i
z
a

v
a
r
i
a
b
l
e
s

S
e

p
u
l
s

c
o
m
e
n
z
a
r

S
A
G
_
m
o
d

S
A
G
,

H
a
r
n
e
r
o
,

C
h
a
n
c
a
d
o
r

P
H
c
_
m
o
d

P
o
z
o
,

B
o
m
b
a
,

H
i
d
r
o
c
i
c
l
n

M
b
o
l
_
m
o
d

M
o
l
i
n
o

d
e

b
o
l
a
s

S
e

r
e
a
l
i
z
a

a
c
c
i
n

s
o
l
i
c
i
t
a
d
a

I
n
t
e
r
r
u
p
c
i
n

S
e

r
e
a
l
i
z
a

a
c
c
i
n

s
o
l
i
c
i
t
a
d
a

I
n
t
e
r
r
u
p
c
i
n

S
e

r
e
a
l
i
z
a

a
c
c
i
n

s
o
l
i
c
i
t
a
d
a

I
n
t
e
r
r
u
p
c
i
n

S
e

a
c
t
u
a
l
i
z
a
n

l
o
s

g
r
f
i
c
o
s

S
i

N
o

S
i

N
o

S
i

N
o

S
i

S
e

c
a
m
b
i
a

a

p
a
n
t
a
l
l
a

d
e

m
o
l
i
e
n
d
a

S
A
G
,

s
e
c
u
n
d
a
r
i
a

o

l
a
y
o
u
t

S
e

p
u
l
s

P
a
u
s
a

N
o

S
i

S
e

p
u
l
s

C
o
n
t
i
n
u
a
r

S
i

S
o
l
i
c
i
t
u
d

c
a
m
b
i
o

p
a
n
t
a
l
l
a

N
o

S
i

S
e

r
e
a
l
i
z
a

e
l

c
a
m
b
i
o

e
n

e
l

c
o
n
t
r
o
l

s
o
l
i
c
i
t
a
d
o
:

c
a
m
b
i
o

e
n

f
l
u
j
o
s

d
e

s
l
i
d
o
s

o

a
g
u
a
,

s
e
t
t
i
n
g
s

d
e

c
o
n
t
r
o
l
e
s

P
I

o

P
&
I
D
s

i
m
p
l
e
m
e
n
t
a
d
o
s
,

a
c
t
i
v
a
c
i
n

d
e
l

c
h
a
n
c
a
d
o
r

S
o
l
i
c
i
t
u
d

a
c
c
i
n

e
n

c
o
n
t
r
o
l
e
s

N
o

S
i

R
e
a
l
i
z
a

e
l

c
a
m
b
i
o

e
n

l
a

v
e
l
o
c
i
d
a
d

d
e

a
c
u
e
r
d
o

a
l

s
e
t
t
i
n
g

d
a
d
o

C
a
m
b
i
o

v
e
l
o
c
i
d
a
d

d
e

s
i
m
u
l
a
c
i
n

N
o

S
i

D
S
A
G
_
s
e
t

C
a
m
b
i
a

e
l

v
a
l
o
r

d
e

l
o
s

p
a
r
m
e
t
r
o
s

d
e
l

s
i
m
u
l
a
d
o
r

S
o
l
i
c
i
t
u
d

C
a
m
b
i
o

c
o
n
f
i
g
u
r
a
c
i
n

S
i

A
c
c
i
o
n
e
s

p
r
o
d
u
c
t
o

d
e

u
n
a

i
n
t
e
r
r
u
p
c
i
n

- 95 -

3.2.3 Mdulo SAG_mod
Para facilitar la lectura de variables se presenta la figura 2.5-1.

SAG_mod comienza asignando las perturbaciones a la moliendabilidad del mineral. Como
se vio en el captulo 2.5.1 esto se realiza determinando el factor de dureza de mineral
en funcin de un ruido coloreado n. Con se determina las tasas especficas de molienda
k
E
y a la granulometra de la alimentacin fresca f.
Si el control PI del hold-up del SAG (w
S
) est activado, se lee el set-point definido por el
usuario y el control PI determina el valor total del flujo de alimentacin fresca (F), si no
est activado el usuario define este valor. El flujos de slidos por intervalos de tamao se
determina como f = F f.
Si esta activada la razn fija de agua, G toma el valor especificado por el usuario, si no
G ser un 30% de F.
El flujo de slido al SAG f
S
ser el valor de f en la iteracin n ms el flujo de rechazo r
calculado en la iteracin n-1, del mismo modo
S
i
q ser la suma de G en iteracin n ms q
r

en n-1.
SAG
k
E
(), c
Harnero
q
A
Chancador
A, B
G
q
o
S

p
S
q
i
A

f
A
q
i
S
+
q
r

f
S
=f+r f=F f()

q
o
A

p
A
D
= p
A

p
A
B
= (1-)p
A
q
i
C

f
C
q
o
C

p
C
q
r
A
r
A q
r
r

- 96 -

Con w
S
en (n-1) se calcula el porcentaje de ocupacin de slidos en el molino J para n-1,
y con este la potencia S en (n-1). Con w
S
y J en n-1 y k
E
en n, se determina las constantes
de molienda k (en n)
Se determina la descarga de la cmara de molienda para n-1, la que depende de w
S
en n-
1. Adems de obtiene la alimentacin a la cmara de molienda en n, la que es la
alimentacin fresca en n ms el rechazo de la iteracin anterior (n-1).
Con los valores anteriores se calcula el hold-up de slidos y agua para n. Finalmente se
determinan las descargas y rechazos para el intervalo n.
El flujo de slidos al harnero ser igual a la descarga del SAG. Se determina el d
50
. Con el
d
50
se determina el valor de la eficiencia del harnero. Con la eficiencia se determina el
rechazo de mineral y con este la descarga de slidos. En forma paralela se determina el
rechazo de agua y la descarga de agua. La descarga del harnero es la alimentacin
fresca del circuito de molienda secundaria, para esto se definen los flujos de slidos duros
y blandos de acuerdo a [2.5-6] y [2.5-7]
=
D
A A
p p
( ) =
B
A A
p 1 p
Si el est activado el uso del Chancador de Pebbles el rechazo del harnero se convierte
en la alimentacin del Chancador, con esta se determinan las matrices de descarga A y
de ruptura B y luego la descarga del Chancador. La descarga del chancador de pebbles
ser el flujo de recirculacin de la molienda SAG. Si no est activado el Chancador el flujo
de recirculacin ser igual al rechazo del harnero.
A continuacin se muestra el diagrama de flujo del ciclo de molienda SAG.

- 97 -

( ) n = +
n1
1
1
2

E
n
n
k
f

Control PI
de W
S

F
n
por usuario
f
n
= F
n
f
n

Set-point W
S
*
por usuario
F
n
por control PI
f
n
= F
n
f
n
Razn fija de
Agua
Si
No

G
n
= 0.3 F
n
G
n
por usuario
No

f
s
n
= f
n
+ r
n -1
q
is
n

= G
n
q
r
n -1

J
n - 1
= J(w
n - 1
)

S
n - 1
= S(J
n - 1
)
( )
1 1
E
n n n n
k k w , J , k

=
SAG
Si

- 98 -

Chancador
Harnero
SAG
( )
1 1
s s s
n n
p p w

=

1
s s s
n n n
f f r
= +

( )
1 1 1 1
s s
n n n n n n
s s s
w w f p k w t

A = +

( )
1
1
SA SA
n n
i o
s s
n n
w w t q q

A = +

s s
n n
r c p =

( )
I
s s
n n
p c p =

o s SA
s n
n
q w = o

( )
s s s
n n
p p w =

= =
A s iA oS
n n n n
f p q q

( )
50 50 A
n n
d d f =

A A A
n n n
p f r =

o i r
A A A
n n n
q q q =

A A A
n n n
r f = q

0 1 =
r A
A
n
n
q . r

( )
50 A A
n n
d q q =

Chancador
conectado
C A
n n
f r =

i r
C A
n n
q q =

( )
n
n
C
A A f =

( )
n
n
C
B B f =

( ) ( )
1
=
C C
n n
n n n
I A I B A p f
o i
C C
n n
q q =

=
C
n
n
r p

r o
n C
n
q q =

=
A
n
n
r r

r r
n A
n
q q =

Si
No
Flujos de recirculacin
Flujos de recirculacin

= =
= =
=
Sec A
n n n
n n n
n n
D D
Sec A A
B D
Sec A A
i o
f p p
f p (1 )p
q q

Flujos a molienda secundaria

- 99 -

3.2.4 Programacin del Mdulo PHc_mod
Para facilitar la lectura de variables se presenta la figura 2.5-2.

El mdulo PHc_mod cubre el funcionamiento del pozo, bombas centrfugas y la batera de
hidrociclones.
Los flujos de entrada a PHc_mod son las descargas de slidos duros y blandos del
harnero, y la descarga de agua del harnero. Al pozo llegan los flujos antes descritos ms
las descargas del molino de bolas. Con los flujos de mineral duro se calcula el hold-up de
duros en el pozo (w
P
D
) y con los flujos de blandos el hold-up de blandos (w
P
B
).
Si est conectado el control PID de +65# en producto final (descarga de hidrociclones), el
controlador lee el set-point y determina el set-point del control PID del porcentaje de
slidos en la alimentacin a los hidrociclones (c
p
). Si el control PID de +65# est
conectado, entonces el control de PID de c
p
se activa y se desactiva el set-point de c
p
por
usuario.
Si est conectado el control PID del porcentaje de slidos en la alimentacin a los
hidrociclones (c
p
) y el control de +65# no est conectado, entonces el set-point de c
p
es
Pozo
w
p
D
, w
p
B

Hidrocicln Molino de
Bolas
G
P
q
i
H

=
+
=
+

H B t p
D B
P P
H H D B
P P
F f c
w w
f F
w w

q
i
P

f
P
q
i
Sec
=
q
o
A

f
D
Sec
= p
A
D

f
B
Sec
= p
A
B

q
o
H

p
H
q
i
M

f
M
D
= r
M
D

f
M
B
= r
M
B

q
o
M

p
M
D

p
M
B

q
r
H

Bomba
k
B1
, k
B2

B
f
v
B
v k k
B B B1 B1
f ( , , , ) P
=
+
=
+

D
D P
H H D B
P P
B
B P
H H D B
P P
w
r r
w w
w
r r
w w

- 100 -

ingresado por el usuario. Ya sea que el set-point est ingresado por el usuario o por el
control PID de +65#, el control determina el agua fresca que ingresa al pozo (G
P
).
Si el control PID de c
p
no est activado (lo que implica que el PID de +65# no est
activado) el agua fresca al pozo (G
P
) es ingresada por el usuario.
Se calcula el flujo volumtrico de alimentacin al pozo (
P
f ) para el intervalo n, para esto
se suma los flujos msicos de duros y blandos. Se determina el flujo volumtrico del a
descarga del pozo (
P
p ) de la iteracin anterior (n-1). Con estos datos se determina la
altura de la pulpa de mineral en el pozo, luego la densidad de la pulpa de mineral en el
pozo (
P
t
) y luego el porcentaje de slidos en la pulpa de mineral (
P
p
c ).
Si el control PID de nivel de la pulpa en el pozo est activado, el usuario ingresa el set-
point y el control PID determina la velocidad de la bomba centrfuga (v
B
). Si no est
activado el control, la velocidad es ingresada por el usuario. Luego se calcula el flujo
volumtrico de la pulpa hacia los hidrociclones.
Con los flujos volumtricos, densidad de la pulpa y porcentaje de slidos en la iteracin n,
se determinan los flujos, presiones y variables del modelo para la iteracin n en el
hidrocicln. Finalmente se determinan los flujos de mineral duro y blando en el rechazo y
rebose, considerando que para cada intervalo de tamao estos son proporcionales a los
hold-up de mineral duro y blando para cada intervalo de tamao en el pozo.

- 101 -

Si
No
No
Control PID
+65#
Set-point +65#
*
por usuario
Control PID
C
p
Si
1
B B B
p Sec M
n n n
f f p
= +
1
D D D
p Sec M
n n n
f f p
= +
( )
1 1
B B B B
p p p p
n n n n
t w w f p

A = + +
( )
1 1
D D D D
p p p p
n n n n
t w w f p

A = + +
1
i i o
p Sec M
n n n
q q q
= +

Set-point c
P
*

por usuario
Set-point c
P
*
por
control PID de +65#

G
P
n
por
control PID

n
n p n
n
p
p i p
m
F
f q G = + +

1
1
1 1
n
n
n n
H
p
tp p
F
p
c

=

= =
D D B B
Sec A Sec A
n n n n
f p f p

=
Sec o
n A
n
i
q q

G
P
n
por usuario
Pozo
Alimentacin desde molienda SAG
Control PID
+65#
Si
No

- 102 -

( ) ( )
1 1
1
n n n n p
p p p p A
h h f p t

= + A
( )
( )
1
1 1
1
1 1
tp
n
n n n n n n p m
n
p i p p
p
tp tp tp A h
F q G t

= + + A +
(

1
1

=
tp
n
m
Pn m tp
n
p
c
No
Control PID
Nivel Pozo
Si
Set-pint de h
P
*
por usuario

v
B
por control PID
( )
1 2 2 1
n
B B B
n n n
B B B
f = k v k P P -
v
B
por usuario

n
H
n n n
B tp p
F f c =
( )
1
i
n n n H
n
B tp p
q f c =
( ) 2
1 1 3 96
o i H H
n H H
n n
p
q , q , c K n = +
r i o
H H H
n n n
q q q =
( ) B p
n n n
H H
f c , = P P
( ) o
n n H
n
H H H
Y Y q , = P
Pozo
Bomba
Hidrocicln

- 103 -

( )
n
n n n
n n
H
B D
H B D p p
p p
F
f w w
w w +
+ =

( )
I 1
n n
H H
H
n
n
H
n
i
H
o
q
p Y F
q
+
| |
| =
|
\ .
n n n
H H H
H
n
n
H
n
i
H
o
q
r f Y f
q
+ =

=
+
B
P
B
H B D
P P
i
n
i i
n n
i i
n n
H
w
w w
r r
=
+
D
P
D
H B D
P P
i
n
i i
n n
i i
n n
H
w
w w
r r
H
n
r
q
Hidrocicln

=
+
B
P
B
H B D
P
i
n
i i
n n
i i
n n
H
P
w
w w
p p
=
+
D
P
D
H B D
P P
i
n
i i
n n
i i
n n
H
w
w w
p p
H
n
o
q
Rechazo alimentacin al molino de Bolas Rebose producto final

- 104 -

3.2.5 Programacin del Mdulo Mbol_mod
Mbol_mod determina los flujos de entrada los que son iguales a los flujos en el rechazo de
los hidrociclones.
Se determinan los hold-up de mineral duro (w
M
D
) y blando (w
M
B
) considerando las matrices
de molienda para duros (B
D
) y blandos (B
B
). Adems se determina el hold-up de agua
(W
MA
).
Finalmente con la tasa de descarga o
M
se determinan los flujos de descarga de slidos
duros, slidos blandos y agua. Estos flujos son descargados al pozo.

B B
M
n n
H
f p =
D D
M
n n
H
f p =
i o
H M
n n
q q =

( )
1 1 1
B B B B B B
M M M M M
n n n n n
t w w f p B w

A = + +
( )
1 1 1
D D D D D D
M M M M M
n n n n n
t w w f p B w

A = + +
( )
1
1
M M
MA MA i o
n n
n n
q q t w w

A = +

B B
M
n n
M M
p w = o
D D
M
n n
M M
p w = o
o MA MA
M n
n
q w = o

Molino de Bolas

- 105 -

3.2.6 Interfaz Grfica del Usuario e Interaccin con los Modelos
El simulador presenta tres pantallas: Una de control y supervisin del circuito de molienda
SAG, otra de control y supervisin del circuito de molienda secundaria y una tercera
donde se muestra un layout de la planta y los principales flujos de mineral y de otras
variables de la planta.
i) Pantalla de control y supervisin de la molienda SAG
Esta pantalla est constituida por dos monitores con grficos de tendencia, una zona para
manejar las variables manipulables, una subpantalla con informacin de la Molienda
Secundaria y una zona con los sistemas de control del simulador. La pantalla de control y
supervisin de la molienda SAG se muestra en la Figura 3.2-2

Figura 3.2-2 Pantalla de Control y Supervisin de la molienda SAG

- 106 -

Grficos
El grfico superior muestra la potencia consumida por el molino SAG en [kW]. El valor
numrico instantneo de la potencia se muestra justo sobre el grfico. Tanto el grfico de
tendencia como el valor de la potencia aparecen en rojo.
El grfico inferior muestra en azul el hold-up de mineral dentro del molino en [t], en violeta
la alimentacin fresca de slidos en [t/h] y en verde el set-point del control automtico de
hold-up del molino SAG. El hold-up y el set-point utilizan la escala de la izquierda que est
en [t], la alimentacin de slidos utiliza la escala en violeta que est al lado derecho y que
indica [t/h].
Controles
Al costado derecho en color caqui tenemos los controles del SAG.
En la parte superior se encuentra una barra de desplazamiento con la cual se controla la
alimentacin fresca de slidos al SAG.
Despus se encuentra una barra de desplazamiento con la cual se controla el agua en la
alimentacin al SAG y una casilla para fijar o no una razn de 30% de agua en la
alimentacin.
Ms abajo hay una casilla para activar el control PI de hold-up, una barra de
desplazamiento para fijar el set-point de hold-up y, en azul, el valor actual del hold-up en
[T].
Por ltimo hay una casilla para activar el Chancador de pebbles.
Subpantalla
Sobre la zona de controles, se encuentra una subpantalla en color verde agua con datos
importantes de la molienda secundaria.

- 107 -

En la subpantalla se muestra el valor numrico instantneo de:
- El +65# del producto final o rebose de los hidrociclones.
- El porcentaje de slidos (Cph [%]) en la alimentacin de los hidrociclones.
- La carga circulante (C.C. [%]).
Comandos del simulador
Debajo de los grficos estn los comandos del simulador.
A la izquierda estn los botones para comenzar y para poner en pausa la simulacin, ms
a la derecha hay una barra de desplazamiento para variar la velocidad aparente de la
simulacin.
A continuacin hay un botn para realizar un refresco de los grficos, en caso de que sea
necesario.
Despus se tiene un men emergente para elegir el intervalo de tiempo mostrado en el
grfico. Los grficos muestran 4 horas de operacin.
En el extremo derecho se muestra el tiempo de simulacin en minutos.
Alarma de sobrecarga
La potencia del molino SAG aumenta a medida que aumenta el holdup hasta un punto
donde esto se revierte. Cuando al aumentar el holdup aumenta la potencia se est en una
condicin de operacin estable, por el contrario, si al aumentar el holdup la potencia
disminuye se est en una condicin inestable. Se dice que el molino est en sobrecarga
cuando se opera en la zona inestable.
En el extremo superior izquierdo de la pantalla de control y supervisin se encuentra una
alarma que se activa cuando el molino entra en sobrecarga.
Produccin
En la parte superior se indica, en negro, la produccin total en toneladas del circuito de
molienda desde que se inici la simulacin hasta el tiempo indicado abajo.

- 108 -

Botones SAG, SEC y Layout
Sobre la subpantalla Molienda Secundaria se encuentran tres botones los que llevan a
las ventanas de molienda SAG (botn SAG), a la molienda secundaria (botn SEC) y al
layout de la planta (botn Layout). Si uno se encuentra en la pantalla de molienda SAG el
botn SAG est desactivado.
ii) Pantalla de control y supervisin de la molienda Secundaria
Esta pantalla est constituida por dos monitores con grficos de tendencia, una zona para
manejar las variables manipulables, una subpantalla con informacin de la Molienda SAG
y una zona con los sistemas de control del simulador. La pantalla de control y supervisin
de la molienda Secundaria se muestra en la Figura 3.2-3

Figura 3.2-3 Pantalla de Control y Supervisin de Molienda Secundaria

- 109 -

Grficos
En el grfico superior se presenta en rojo el porcentaje de slidos en la alimentacin de
los hidrociclones, en verde el set-point del control PID de este porcentaje de slidos, y en
violeta el agua que se agrega al pozo, que es una variable manipulada del control.
El grfico inferior muestra en azul el porcentaje de gruesos +65# en el producto final, es
decir en el rebose de los hidrociclones, y en verde el set-point del control PID del
porcentaje de +65#. Este control PID es un control en cascada con el control PID de
porcentaje de slidos, en que la variable manipulada es el set-point del control PID de
porcentaje de slidos.
Sobre los grficos se puede observar la alarma de sobre carga en el molino SAG, la
produccin total y la carga circulante en el circuito de molienda secundaria.
Controles
En el panel de controles se tiene, en primer lugar, una barra de desplazamiento con la
que se fija la cantidad de agua que se agrega al pozo.
Despus se encuentra la casilla para activar el control de porcentaje de slidos en la
alimentacin a hidrociclones (Cph), la barra de desplazamiento para fijar el set-point y el
valor actual del porcentaje de slidos. Como el control acta variando la cantidad de agua
que se agrega el pozo, si est activo la barra de desplazamiento de cantidad de agua al
pozo no aparecer y el valor no estar accesible.
Luego se ubica la casilla para activar el control PID del porcentaje de +65# en la salida de
los hidrociclones, la barra de desplazamiento para fijar el set-point y el valor actual del
porcentaje. Como este control est en cascada con el de porcentaje de slidos, si este
control est activo no est accesible el set-point del control PID de porcentaje de slidos
Por ltimo aparece el control PID de nivel del pozo, este es un control que est activado
en forma predeterminada y que controla el nivel del pozo manipulando la velocidad de la
bomba centrfuga. Aqu adems de estar la barra de desplazamiento para el set-point del
nivel del pozo estn los valores actuales del nivel del pozo y de la velocidad de la bomba.

- 110 -

Subpantalla
En la subpantalla que est sobre los controles, aparecen valores importantes de la
molienda SAG: el valor de la potencia, del hold-up del molino y el estado del chancador
(conectado/desconectado).
Botones SAG, SEC y Layout, y controles
Los botones SAG, SEC y Layout as como el rea con los controles de simulador, son
los mismos descritos en el punto anterior.
iii) Pantalla con el Layout de la Planta
Esta pantalla se muestra un layout de la planta y los valores actuales de los principales
flujos y variables de la Planta. La pantalla del layout de la Planta se muestra en la figura
Figura 3.2-4

Figura 3.2-4 Pantalla del Layout de la Planta

- 111 -

Siguiendo el layout, los flujos y variables mostradas en la pantalla son las siguientes:
- F: flujo msico de slidos en [T/h] en la alimentacin del SAG.
- Fa: flujo msico de agua en [T/h] en la alimentacin del SAG.
- Hold-up: el hold-up del molino SAG en [T].
- Pot: potencia consumida por el motor del Molino SAG en [KW].
- P: flujo msico de slidos en [T/h] en la descarga del harnero.
- Pa: flujo msico de agua en [T/h] en la descarga del harnero.
- R: flujo msico de slidos en [T/h] en el rechazo del harnero.
- Ra: flujo msico de agua en [T/h] en el rechazo del harnero.
- hp: nivel del pozo en [m].
- G: flujo msico de agua en [T/h] que es adicionada al pozo.
- V: velocidad de la bomba centrfuga en [RPM].
- Fc: flujo msico de slidos en [T/h] en la alimentacin a hidrociclones.
- Fac: flujo msico de agua en [T/h] en la alimentacin a hidrociclones.
- Cph: porcentaje de slidos en [%]en la alimentacin a hidrociclones.
- Pc: flujo msico de slidos en [T/h] en rebose de hidrociclones.

- 112 -

- Pac: flujo msico de agua en [T/h] en rebose de hidrociclones.
- +65#: +65# en [%] en rebose de hidrociclones.
- Fb: flujo msico de slidos en [T/h] en alimentacin del molino de bolas.
- Fab: flujo msico de agua en [T/h] en alimentacin del molino de bolas.
- C.C.: Carga circulante en [%].
- Pb: flujo msico de slidos en [T/h] en descarga del molino de bolas.
- Pab: flujo msico de agua en [T/h] en descarga del molino de bolas.
En el extremo inferior izquierdo adems de los botones para continuar y poner en pausa
la simulacin, hay una casilla para utilizar o no el Chancador de pebbles.
El circuito de molienda secundaria tiene un pozo y dos circuitos paralelos formados cada
uno por una bomba centrfuga, una batera de hidrociclones y molino de bolas. Por
razones visuales en el layout de la planta se ha representado un circuito equivalente a los
dos circuitos reales, tal como se indica en el extremo inferior derecho de la pantalla. Por lo
tanto, los flujos mostrados son la suma de los flujos de los dos circuitos.
3.2.2 Exportacin de Datos
El simulador permite grabar en un archivo electrnico los valores de los flujos y otras
variables de inters obtenidas durante la simulacin, de manera que puedan hacerse
anlisis posteriores con estos datos.
Con la opcin Guardar del men Simulador del simulador, se genera un archivo ASCII
con los valores en cada instante de las variables. El archivo es tipo tabla, en la primera fila

- 113 -

estn los nombres de las variables y en las restantes los valores. Las columnas son de un
ancho fijo de 9 caracteres.
Planillas electrnicas como Excel permiten fcilmente transformar estos datos ASCII en
una tabla de valores para analizarla. En el ANEXO B: Simulador molienda SAG - Manual
del Usuario, se explica cmo exportar los datos a una hoja Excel.
Los datos que se guardan en forma electrnica va la opcin Guardar son los que se
listan a continuacin. Se muestra a la izquierda la descripcin de la variable que aparece
en el archivo ASCII y a la derecha su explicacin.
Variable Descripcin
F[T/h] Flujo de slidos en alimentacin al SAG en [T/h]
Fa[T/h] Flujo de agua en alimentacin al SAG en [T/h]
P[T/h] Flujo de slidos en descarga del harnero en [T/h]
Pa[T/h] Flujo de agua en descarga del harnero en [T/h]
R[T/h] Flujo de slidos en rechazo del harnero en [T/h]
Ra[T/h] Flujo de agua en rechazo del harnero en [T/h]
Pot[KW] Potencia consumida por el motor del SAG en [KW]
W[T] Hold-up del SAG en [T]
PID_W Estado del control PI del SAG (0 desconectado 1 conectado)
spW[T] Set-point del control PI de hold-up en [T]
+65#_P +65# del bajo tamao del harnero en [/1]
Chancado Estado del chancador de pebbles (0 desconectado 1 conectado)
Tabla 3.2 - 1 Variables guardadas por el simulador (continua en pgina siguiente)

- 114 -

Variable Descripcin
Raz.Agua Estado del control de razn fija de agua en un 30% en alimentacin
fresca al SAG (0 desconectado 1 conectado)
hp[m] Nivel del pozo en [m]
PID_hp Estado del control PID de nivel del pozo (0 desconectado 1 conectado)
sphp[m] Set-point del control PID de nivel del pozo en [m]
V[RPM] Velocidad de la bomba centrfuga del pozo en [RMP]
G[T/h] Flujo de agua agregada al pozo en [T/h]
Cph Porcentaje de slidos en alimentacin a hidrociclones en [/1]
PID_Cph Estado del control PID de porcentaje de slidos en alimentacin a
hidrociclones (0 desconectado 1 conectado)
spCph Set-point del control PID de Cph en [/1]
Fc[T/h] Flujo de slidos en alimentacin a hidrociclones en [T/h]
Pc[T/h] Flujo de slidos en rebose de hidrociclones en [T/h]
+65#_Pc +65# en rebose de hidrociclones en [/1]
PID_+65# Estado del control PID de +65# en rebose de hidrociclones (0
desconectado 1 conectado)
sp+65#Pc Set-point del control PID de +65# en rebose de hidrociclones en [/1]
CpPc Porcentaje de slidos en producto final (rebose de hidrociclones) en
[/1]
+65#_Pb +65# en descarga del molino de bolas en [/1]
C.C. Carga circulante del circuito de molienda secundaria en [%] (razn entre
la descarga del harnero del molino SAG y las descargas de los molinos
de bolas)
mu Dureza del mineral (0 = el ms blando, 1 = el ms duro)
Tabla 3.2-1 Variables guardadas por el simulador (continuacin)

- 115 -

3.3 AJUSTES DE PARMETROS
3.3.1 Antecedentes generales
Como se mencion anteriormente, el objetivo final buscado en esta memoria es crear un
simulador de una planta de molienda que sea cualitativamente adecuado a la realidad, y
no al ajuste de una planta especfica. Un modelo que cualitativamente se comporta como
se esperara de una planta real permite implementar y estudiar en forma adecuada
distintos sistemas de control. Adems un modelo abierto permite estudiar la variacin al
cambiar algn subsistema, por ejemplo un chacador de pebbles, si se cuenta con alguno
especfico parametrizado.
En funcin de lo anterior, la mayor parte de los ajustes realizados fue considerando
caractersticas y puntos de operacin de la planta de molienda A2 de Chuquicamata, pero
no todos los subsistemas fueron ajustados en base a esta.
El ajuste de parmetros fue realizado en forma independiente para cada subsistema:
molino SAG, harnero, chancador de pebbles, bomba centrfuga, hidrociclones y molino de
bolas. El modelo del pozo es un balance volumtrico considerndose como mezclador
perfecto, y no hay parmetros que identificar.
Las constantes de cada modelo fueron determinadas en punto de operacin de equilibrio.
Posteriormente se realizaron sesiones de trabajo para hacer otros ajustes menores. En
estas sesiones se utiliz los conocimientos expertos de A. Casali para lo concerniente a
los modelos de la planta de molienda y de G. Gonzlez en lo referente al sistema de
control.
Para la identificacin de los parmetros de los molinos y de los hidrociclones, la principal
fuente utilizada fueron los muestreos realizados en la Concentradora A2 de Chuquicamata
el 30 de junio de 1996. La informacin obtenida de estos muestreos corresponde a flujos
totales y distribuciones granulomtricas en el molino SAG 17, molinos de bolas 17A y
17B, y las bateras de hidrociclones 17A y 17B.

- 116 -

El modelo desarrollado en esta memoria est compuesto, al igual que la concentradora
A2, de dos lneas de molienda de bolas las que se asumen son iguales, es decir el valor
de los parmetros son los mismos. Las medidas realizadas en la concentradora A2,
difieren levemente en cada una de las lneas de molienda de bolas. Para la determinacin
de los parmetros de los molinos de bolas e hidrociclones se utiliz la informacin del
molino de bolas 17A y de la batera de hidrociclones 17A.
Para la identificacin de los parmetros del harnero del molino SAG y de las bombas
centrfugas se utiliz informacin de catlogos ms ajustes posteriores realizados tras
anlisis de los resultados obtenidos del simulador.
La concentradora A2 no tena un chancador de pebbles cuando se realizaron los
muestreos utilizados para la identificacin de los parmetros de los molinos e
hidrociclones. Para la identificacin de parmetros del chancador de pebbles se utiliz en
forma inicial los valores recomendados por A. Casali en el modelo de chancador
desarrollado por l basado en una simplificacin del modelo de Andersen y Whiten.
Posteriormente se realizaron, al igual que con el harnero, ajustes a estos parmetros va
anlisis de resultados obtenidos del simulador.
El modelo desarrollado considera flujos msicos de mineral para doce intervalos de
tamao de slidos. Para cada intervalo se define el tamao representativo (X
1
a X
12
) que
corresponde a la medida geomtrica del intervalo, exceptuando el tamao X
12
(el tamao
menor) que es calculado como la media aritmtica del intervalo y X
1
(el tamao mayor)
para el que se asume 8.
La Tabla 3.3 - 1 muestra los tamaos representativos utilizados en el simulador y el
intervalo que representan.

- 117 -

Tamao Intervalo Intervalo en [mm] Tamao Rep. X
i
[mm]
X
1
+6 >152.4 203.2
X
2
-6 , +4 (101.6, 152.4] 124.4341
X
3
-4 , +2.1 (53.3, 101.6] 73.5886
X
4
-2.1 , +3/4 (19.1, 53.3] 31.9066
X
5
-3/4 , +1/2 (12.7, 19.1] 15.5747
X
6
-1/2 , +6# (3.35, 12.7] 6.5227
X
7
-6# , +10# (1.7, 3.35] 2.3864
X
8
-10# , +20# (0.85, 1.7] 1.2021
X
9
-20# , +35# (0.425, 0.85] 0.6010
X
10
-35# , +65# (0.212, 0.425] 0.3002
X
11
-65# , +150# (0.106, 0.212] 0.1499
X
12
-150# [0, 0.106] 0.053
Tabla 3.3 - 1 Definicin de intervalos de tamaos de partculas y el tamao representativo

- 118 -

3.3.2 Parmetros del Molino SAG
i) Datos para la parametrizacin
Los parmetros a identificar en el molino SAG son los siguientes
- o
S
: parmetro para determinar la descarga de slidos, ver frmula (2.2-29)
-
S1
,

S2
: parmetros para determinar la potencia consumida por el molino, ver
frmula (2.2-21)
- c: matriz con las tasas de rechazo de slidos, ver frmula (2.2-24)
- k
E
: matriz con las tasas especificas de molienda ver ecuaciones (2.2-12),
(2.2-14) y (2.2-17). Estas tasas son utilizadas para determinar las tasas de
molienda k del molino
- k
S1
,

k
S2
: parmetros para determinar la tasa de descarga de agua del molino, ver
ecuacin (2.2-32)
Las dimensiones fsicas del molino SAG utilizado son:
- Dimetro: d
S
= 32 (9.75 [m])
- Largo: l
S
= 17 (5.18 [m])
- Volumen interno: V
S
= 359 [m
3
]
El punto de operacin utilizado est determinado por los siguientes valores:
- Alimentacin total: F
S
= 1332 [t/h]. Fuente: diagrama de flujo y balance de la
concentradora A2 (ver Figura 3.3-1)

- 119 -

- Distribucin granulomtrica de la alimentacin al molino: de acuerdo a muestreo
de alimentacin al SAG 17 de Chuquicamata realizada el 19/06/96 (ver
Figura 3.3-2)
- Distribucin granulomtrica de la descarga del molino: de acuerdo a muestreo de
la concentradora A2 de Chuquicamata realizada el 30/06/96 (ver Figura 3.3-3)
Los siguientes valores utilizados corresponden a valores tpicos esperados para el punto
de operacin indicado arriba. Estos valores fueron obtenidos de la experiencia de
A. Casali.
- Potencia consumida: S = 6500 [kW]
- Porcentaje del volumen interno del molino ocupado por los slidos: J = 40%
- Porcentaje del volumen interno del molino ocupado por bolas: J
B
= 8%
- Carga de agua: W
SA
= 60 [t]
Adems se consideraron los siguientes valores tpicos:
- Densidad del mineral:
m
= 2.7 [t/m
3
]
- Densidad de las bolas:
B
= 7.1 [t/m
3
]
- Porosidad equivalente del mineral: c = 0.4

- 120 -

Figura 3.3-1 Diagrama de flujo y balance Concentradora A2 Chuquicamata

- 121 -

Figura 3.3-2 Distribucin granulomtrica - Alimentacin SAG 17, Concentradora A2 Chuquicamata

- 122 -

Figura 3.3-3 Distribucin granulomtrica - Descarga SAG 17, Concentradora A2 Chuquicamata

- 123 -

ii) Determinacin de parmetros
S1
,

S2

De la ecuacin de hold-up de slidos de Austin (2.2-22) se tiene:
| |
1 1 0 4
186
S S m B m
S
W V ( )J ( . )J
W T
= c
=

El hold-up de bolas lo determinamos de la siguiente forma:
( )
SB B S B
SB
W J V
W 204 T
=
=

De la ecuacin [2.2-20] se tiene que la densidad del material equivalente es:
3
3
186 60 204
0 4 359
3 134
S SA SB
S
S
S
W w w [T] [T] [T]
J V . [m ]
. [T / m ]
+ + + +
= =
=

Por otro lado, analizando la ecuacin de Austin para la potencia [2.2-21]:
2.5
S1 S S S2 S
S d l J(1 J) =
se tiene que como S0 entonces:
S2 S2 S2
1
1 J 0 2.50
J
> s s

- 124 -

Inicialmente se consider para este parmetro el valor que se tiene del molino piloto de
Chuquicamata:
S2
1.31 =
Reemplazando en la ecuacin de potencia:
2.5
S1 S S S2 S
2.5
S1
S1
S d l J(1 J)
6500 9.75 5.18 0.4(1 1.31 0.4)3.134
7.08
=
=
=

Tras anlisis posteriores con el simulador se ajust
S2
a 1.10, de esta forma los
parmetros en el simulador son:
S1 S2
6.01 1.10 = =
iii) Determinacin de las tasas especficas de molienda k
E
y de rechazo c
Los grficos de distribucin granulomtrica para la carga y descarga del molino SAG de
Chuquicamata (ver Figura 3.3-2 y Figura 3.3-3) fueron analizados grficamente de forma
de determinar una cantidad suficiente de puntos que permitiese una buena aproximacin
al considerar una funcin lineal entre ellos.
El anlisis grfico es el siguiente:

- 125 -

Alimentacin Molino SAG
Puntos considerados en la linealizacin por tramos:
Tamao %
13.600 10%
20.000 20%
27.000 30%
30.000 35%
36.000 40%
53.000 50%
80.000 60%
130.000 70%
210.000 80%
Tabla 3.3 - 2 Ajuste de parmetros Molino SAG: puntos para
linealizacin alimentacin

- 126 -

Descarga Molino SAG

Puntos considerados en la linealizacin por tramos:
Tamao %
45 19%
109 27%
235 37%
480 51%
1.170 71%
2.440 83%
5.000 88%
10.000 91%
13.000 91%
30.000 98%
Tabla 3.3 - 3 Ajuste de parmetros Molino SAG: puntos
para linealizacin descarga

- 127 -

Los grficos de granulometras analizados, corresponden al acumulado pasante, es
decir para un tamao X, se tiene el valor de X ms todos los tamaos menores. Lo que se
necesita son los valores para cada tamao representativo y para los acumulados
retenidos por malla (los de tamao X y mayores).
Del anlisis grfico se determina el porcentaje acumulado pasante para cada tamao
representativo, con lo que se determina el porcentaje para cada tamao representativo y
el porcentaje retenido por malla.

Tamao X
i
Acumulado pasante Valor para tamao Xi Acumulado por malla
Descripcin [m] Alimentacin Descarga Alimentacin Descarga Alimentacin Descarga
X
1
203.200 100,0% 100,0% 31,1% 0,0% 31,1% 0,0%
X
2
124.434 68,9% 100,0% 11,3% 0,0% 42,4% 0,0%
X
3
73.616 57,6% 100,0% 21,1% 0,0% 63,4% 0,0%
X
4
31.877 36,6% 100,0% 23,5% 8,0% 86,9% 8,0%
X
5
15.575 13,1% 92,0% 8,3% 2,9% 95,2% 10,9%
X
6
6.523 4,8% 89,1% 3,0% 6,6% 98,2% 17,5%
X
7
2.386 1,8% 82,5% 0,9% 11,2% 99,1% 28,7%
X
8
1.202 0,9% 71,3% 0,4% 16,9% 99,6% 45,6%
X
9
601 0,4% 54,4% 0,2% 13,9% 99,8% 59,4%
X
10
300 0,2% 40,6% 0,1% 10,4% 99,9% 69,8%
X
11
150 0,1% 30,2% 0,1% 10,2% 100,0% 80,0%
X
12
53 0,0% 20,0% 0,0% 20,0% 100,0% 100,0%
Tabla 3.3 - 4 Ajuste de parmetros molino SAG: Distribucin granulomtrica

- 128 -

Distribucin granulomtrica Alimentacin y Descarga SAG
Porcentaje acumulado retenido por malla
0%
20%
40%
60%
80%
100%
10 100 1,000 10,000 100,000 1,000,000
Al i mentaci n
Descarga

Figura 3.3-4 Ajuste de parmetros Molino SAG: Porcentaje distribucin granulomtrica en
alimentacin y descarga para los tamaos representativos
Las distribuciones granulomtricas analizadas corresponden a los flujos de entrada y
salida del SAG y no a los flujos de entrada a la cmara de molienda. Recordando la figura
2.2-1:

Se tiene que los flujos analizados corresponden a
S
f y
S
p , y no
S
f
y
S
p .
Las tasas de rechazo son debidas a la parrilla interna del molino SAG. De los grfico de
distribucin granulomtrica se ve que no hay descarga de slidos para los tamaos X
1,
X
2

Cmara de
Molienda
Clasificacin
S
f
S
f
~
S
p
~
S
p
S
r

- 129 -

y X
3
, lo que quiere decir que son rechazados en un 100% y luego su tasas de rechazo son
igual a 1. Del mismo modo, el 100% del tamao X
12
es descargado del molino, por lo que
su tasa de rechazo es 0. Las tasas de rechazo para los tamaos restantes fueron
definidas de acuerdo a valores tpicos esperados. Las tasas de rechazo consideradas
son:
c = [1.00, 1.00, 1.00, 0.90, 0.40, 0.35, 0.30, 0.25, 0.20, 0.15, 0.01, 0.00]
El flujo total en el punto de equilibrio analizado es de 1332 T/h, luego los flujos para cada
tamao representativo para la carga y descarga del molino corresponden a la distribucin
granulomtrica determinada arriba por 1332 T/h, luego:
f = [414.4, 149.9, 280.7, 312.7, 110.4, 40.5, 11.6, 5.9, 2.9, 1.5, 0.9, 0.5] T/h
p = [0.0, 0.0, 0.0, 106.8, 38.9, 87.5, 149.1, 224.8, 184.8, 138.1, 135.6, 266.4] T/h
Tenemos que:
( )
S
i
S i S S
i i i
i
p
p 1 c p p
1 c
= =

Luego tenemos los rechazos y los flujos de entrada a la cmara de molienda:
S i S
i i
S S S
i i i
r c p
f f r
=
= +

Sin embargo para los tamaos X
1,
X
2
y X
3
no se puede usar la misma metodologa ya que
sus tasas de rechazo son igual a 1, y luego no se pueden determinar las descargas de la
cmara de molienda para ellos. Para estos flujos se consider que cada uno tena una
carga circulante de 300%. Es decir:

- 130 -

Para j = 1, 2, 3
j j
j j j
j j
f 3 f
p f f
r p
=
=
=

Como determinamos todos los flujos de descarga de la cmara de molienda, utilizando la
ecuacin (2.2-28) se pueden determinar los hold-up dentro de la cmara de molienda para
cada tamao representativo.
S
S S S S
i i i i
S S S
W 1 1
p w w p
P W P
= =

Los resultados son:
X
i
S
f [T/h]
S
f
[T/h]
S
p [T/h]
S
p [T/h]
S
r [T/h]

S
w [T]
X
1
414.4 1,243.3 828.9 0.0 828.9 36.2
X
2
149.9 449.6 299.7 0.0 299.7 13.1
X
3
280.7 842.0 561.3 0.0 561.3 24.5
X
4
312.7 1,273.6 1,067.7 106.8 960.9 46.6
X
5
110.4 136.4 64.9 38.9 26.0 2.8
X
6
40.5 87.6 134.6 87.5 47.1 5.9
X
7
11.6 75.5 212.9 149.1 63.9 9.3
X
8
5.9 80.8 299.7 224.8 74.9 13.1
X
9
2.9 49.2 231.0 184.8 46.2 10.1
X
10
1.5 25.8 162.5 138.1 24.4 7.1
X
11
0.9 2.3 137.0 135.6 1.4 6.0
X
12
0.5 0.5 266.4 266.4 0.0 11.6
Total = 1,332 4,267 4,267 1,332 2,935 186
Tabla 3.3 - 5 Ajuste de parmetros molino SAG: Flujos msicos
Dado que determinamos los flujos de entrada y de salida a la cmara de molienda as
como el hold-up para cada intervalo, podemos determinar los flujos y hold-ups retenidos
por malla. Con estos valores determinamos las tasas de molienda para los flujos retenidos
por malla (ecuacin (2.2-4)) y por lo tanto las constantes especficas de molienda
(ecuacin (2.2-18))

- 131 -

A A A A A
S S S i S
i i i i
d
w f p k w
dt
=

Como estamos en un punto de equilibrio
A
S
i
d
w 0
dt
= , luego:
A A
S S
A i i
i A
S
i
f p
k
w

A E E A S
i i i i
S
W S
k k k k
W S
= =
Los resultados son:
X
i
A
S
f
[T/h]
A
S
p [T/h]
A A
S S
f p
[T/h]
A
S
w [T]
A
k [1/h]
E
k [T/kW h]
X
1
1,243.3 828.9 414.4 36.2 11.46 0.33
X
2
1,692.9 1,128.6 564.3 49.2 11.46 0.33
X
3
2,534.9 1,689.9 845.0 73.7 11.46 0.33
X
4
3,808.5 2,757.6 1,050.9 120.3 8.73 0.25
X
5
3,944.9 2,822.5 1,122.4 123.1 9.11 0.26
X
6
4,032.5 2,957.1 1,075.4 129.0 8.34 0.24
X
7
4,108.0 3,170.0 938.0 138.3 6.78 0.19
X
8
4,188.8 3,469.8 719.1 151.4 4.75 0.14
X
9
4,238.0 3,700.8 537.2 161.5 3.33 0.10
X
10
4,263.8 3,863.3 400.6 168.6 2.38 0.07
X
11
4,266.2 4,000.3 265.9 174.5 1.52 0.04
X
12
4,266.7 4,266.7 0.0 186.2 0.00 0.00
Tabla 3.3 - 6 Ajuste de parmetros molino SAG: Determinacin constantes de molienda

- 132 -

iv) Determinacin parmetro de descarga total de slidos o
S

De acuerdo a (2.2-29)
S S S
S
S
P W
4267 186
312.72
= o
= o
o =

v) Determinacin de parmetros k
S1
,

k
S2
para descarga de agua
La tasa de descarga de agua o
S
queda determinada por la ecuacin (2.2-32):
2
1
4
S
S S
S
W
k
o = k +
Dado que no se contaba con informacin de la concentradora A2, se consider los valores
del molino piloto, ms ajustes posteriores en base a pruebas del simulador. Los
parmetros considerados son:
1
2
3 77
0 71
k =
k =
S
S
.
.

3.3.3 Parmetros del Harnero
Los parmetros a identificar en el harnero son los siguientes:
- k
A1
, k
A2
: parmetros para determinar la eficiencia q
A
, ver ecuacin (2.2-42)
- k
A3
: parmetro para determinar el d
50
La apertura del harnero considerado es (12.7 mm), que corresponde al harnero del
molino SAG de la concentradora A2 de Chuquicamata.

- 133 -

Se consideraron los valores utilizados por Casali en su adaptacin del modelo de Karra.
Los valores son:
- k
A1
= -0.638
- k
A2
= 8.98
- k
A3
= 0.00013
3.3.4 Parmetros del Chancador de Pebbles
Los parmetros a identificar en el Chancador de Pebbles son los siguientes:
-
C0
,
C1
,
C2
,
C3
,
C4
,
C5
,
C6
,
C7
: parmetros para determinar las tasas de
clasificacin. Ver ecuaciones (2.2-54), (2.2-55), (2.2-56).
-
C8
,
C9
,
C10
,
C11
,
C12
: parmetros para determinar la matriz de ruptura. Ver
ecuaciones (2.2-57), (2.2-58), (2.2-59).
Se consideraron los valores utilizados por Casali en su adaptacin del modelo de
Andersen & Whiten [5]. Los valores son:
-
C0
= 0.1
-
C1
= 1
-
C2
= 0
-
C3
= 0
-
C4
= 0.1
-
C5
= 2
-
C6
= 0
-
C7
= 0
-
C8
= 0.9
-
C9
= 0.03
-
C10
= 1.2
-
C11
= 0.05
-
C12
= 0.01

El ajuste corresponde a un Chancador de pebbles con una abertura de descarga (closed
setting) de:
- d
CE
= (12.7 mm)

- 134 -

3.3.5 Parmetros de las Bombas Centrfugas
Los parmetros a identificar en las bombas centrfugas son los siguientes:
- k
B1
, k
B2
: parmetros de la ecuacin de flujo volumtrico en funcin de velocidad y
presin, ver ecuacin (2.2-69)
Se consider las curvas de eficiencia de bombas centrfugas mostradas en SME Mineral
Processing Handbook [11] captulo 4, pag. 10-157. La siguiente grfica con las curvas de
desempeo de la bomba:

Figura 3.3-5 Curvas de eficiencia de la bomba centrfuga

- 135 -

De (2.2-69) tenemos:
| |
B B1 B B2 B2 B1
f = k v k P P
La presin en la bomba es:
B B2 B1
= P P P
Luego:
( )
B B B B1 B B2 B
f , v = k v k P P
De la ecuacin de la recta para 500 rpm es:
( )
1 1
B B B B B
f 500, m m f = + P P P
Se considera puntos en torno al flujo msico de equilibrio considerado en la
concentradora A2 de Chuquicamata, esto es un flujo volumtrico en las bombas de
2400 m
3
/hr. Del grfico se consideran los siguientes puntos:
| | | |
| | | |
1 1
2 2
B B
B B
80 ft f 12000 gpm
88 ft f 10000 gpm
= =
= =
P
P

Se tiene que:
2 1
3
2 1
B B
ft mhr
gpm
m
B B
m 250 2236
f f

= = = ( (

P P

Se observa que:
B2
m k = -
3
mhr
B2 m
2236

k (

=

- 136 -

Se observa que:
| |
1 1
B1 B B
500 rpm m f k + = - P
3
gpm
m s
rpm hr B1
16 218 k = ( (

=
3.3.6 Parmetros de las bateras de Hidrociclones
Los parmetros a identificar en el modelo de hidrociclones son los siguientes:
- R
CH
: Rechazo de Agua ver ecuacin (2.2-81)
- k
H2
: parmetro para determinar el flujo de agua en el rebose q
o
H
, ver ecuacin (2.2-
83)
- o
H
: parmetro de Plitt para determinar matriz de clasificacin Y
H
, ver ecuacin (2.2-
87)
- k
H3
: parmetro para determinar el d
50
- k
H1
: parmetro para determinar la presin de la pulpa en la alimentacin de la
batera de hidrociclones P
H
Las caractersticas de la batera de hidrociclones utilizada para el ajuste de parmetros
son las de la batera de hidrociclones 17A Concentradora A2 Chuquicamata:
- Nmero de hidrociclones, n
H
= 6
- Dimetro del buscador de vrtice (vortex finder), d
vf
= 26.416 cm
- Dimetro del orificio de salida (spigot), d
sp
= 15.24 cm
Los siguientes valores tpicos de operacin fueron considerados:

- 137 -

- Presin de pulpa de mineral en alimentacin hidrociclones, P
H
= 3 bar
- Dimetro para el cual el 50% de las partculas pasan al rechazo, d
50
= 30 m
El punto de operacin utilizado est determinado por los siguientes valores:
Del diagrama de flujo y balance de la concentradora A2 (ver Figura 3.3-1)
- Alimentacin total de slidos: F
H
= 2400 [T/h]
- Rebose total de slidos: P
H
= 600 [T/h]
- Rechazo total de slidos: R
H
= 1800 [T/h]
Del muestreo de la concentradora A2 de Chuquicamata 30/06/96
- Densidad pulpa en la alimentacin:
t FH
= 1.78 [T/m
3
]
- Densidad pulpa en el rebose:
t PH
= 1.53 [T/m
3
]
- Densidad pulpa en el rechazo:
t RH
= 1.97 [T/m
3
]
- Porcentaje de slidos en la alimentacin: c
p FH
= 68%
- Porcentaje de slidos en el rebose: c
p PH
= 55%
- Porcentaje de slidos en el rechazo: c
p RH
= 78%

- 138 -

- Distribucin granulomtrica:
Alimentacin f
H
Rechazo r
H
Rebose p
H

Tamao
[mm]
Malla
Tyler
%
Retenido
%
Acumulado
retenido
%
Acumulado
pasante
%
Retenido
%
Acumulado
retenido
%
Acumulado
pasante
%
Retenido
%
Acumulado
retenido
%
Acumulado
pasante
6700 3 0,2 0,2 99,8 1,9 1,9 98,1 0,0 0,0 100,0
4750 4 0,4 0,6 99,4 2,3 4,2 95,8 0,0 0,0 100,0
3350 6 0,9 1,5 98,5 3,2 7,4 92,6 0,0 0,0 100,0
2360 8 1,5 3,0 97,0 4,3 11,7 88,3 0,0 0,0 100,0
1700 10 1,9 4,9 95,1 4,3 16,0 84,0 0,0 0,0 100,0
1180 14 4,0 8,9 91,1 8,6 24,6 75,4 0,0 0,0 100,0
850 20 4,7 13,6 86,4 9,8 34,4 65,6 1,1 1,1 98,9
600 28 6,3 19,9 80,1 12,9 47,3 52,7 3,5 4,6 95,4
425 35 7,4 27,3 72,7 12,5 59,8 40,2 8,9 13,5 86,5
300 48 6,5 33,8 66,2 8,9 68,7 31,3 11,0 24,5 75,5
212 65 5,6 39,4 60,6 6,7 75,4 24,6 10,0 34,5 65,5
150 100 4,1 43,5 56,5 4,5 79,9 20,1 9,0 43,5 56,5
106 150 3,2 46,7 53,3 3,2 83,1 16,9 7,7 51,2 48,8
75 200 2,9 49,6 50,4 2,5 85,6 14,4 7,0 58,2 41,8
53 270 2,1 51,7 48,3 1,8 87,4 12,6 5,3 63,5 36,5
45 325 1,1 52,8 47,2 0,9 88,3 11,7 2,9 66,4 33,6
23 -325 47,2 100,0 0,0 11,7 100,0 0,0 33,6 100,0 0,0
Tabla 3.3 - 7 Ajuste de parmetros hidrociclones: Distribucin granulomtrica
concentradora A2 de Chuquicamata
ii) Determinacin del Rechazo de Agua R
CH

Con los datos de flujo de slidos y porcentaje de slidos en la pulpa, se determina los
flujos de agua:
H
p FH
i H
p FH
1 c
q F 1129[T / h]
c
= =
H
p PH
o H
p PH
1 c
q P 491[T / h]
c
= =

- 139 -

H
p RH
r H
p RH
1 c
q R 508[T / h]
c
= =
Luego, utilizando (2.2-81), se tiene:
0 45 = =
r
CH
i
H
H
q
.
q
R
iii) Determinacin del parmetro k
H2
para determinar el flujo de agua en el rebose
Utilizando (2.2-83)
( )
2
2
1 1 3 96
1
1 1 3 96 122 55
= + k
k = =
o i p FH H H
H i o p FH H
H
H H
H H
q . q ( . c )n
. q q . c n .
n

iv) Determinacin del parmetro k
H1
para determinar presin de pulpa en la
alimentacin
La ecuacin (2.2-89) determina la presin en la alimentacin a la batera de hidrociclones
en funcin del flujo volumtrico en la alimentacin
H
f su porcentaje de slidos y
constantes. La ecuacin es:
( )
2
0 125
1
13 725
1 48 1
H
H .
vf H p H
f
.
. d n c
| |
|
=
|
k
\ .
P
El flujo volumtrico en la alimentacin es:
= =
3 H
H
t FH p FH
F
f 1983 [m h]
c

- 140 -

Se considera una presin tpica a la entrada de la batera de hidrociclones de
P
H
= 3 bar. Luego:
( )
3
0 5
1 0 125
13 725
20 64
1 48 1
(
k = =

.
H
m
H . h kg
H
vf H p FH
f
.
.
. d n c
P

v) Determinacin del parmetro k
H3
para determinar el d
50

Para determinar k
H3
se considera un valor tpico de d
50
= 30 m. Luego utilizando la
ecuacin [2.2-88] se obtiene:
| |
50 3
3
30
6 6 8 9 10 7 47 2
2 56
| |
= + + k =
|
|
\ .
k =
o
sp
vf H
H
H
H
H
q
d
d
d exp m
. . . . n
.
P

vi) Determinacin del parmetro de Plitt o
H
para determinar la matriz de seleccin Y
H

Si se asume un valor cualquiera para el parmetro de Plitt, con los otros parmetros y
datos que tenemos es posible determinar todos los flujos (alimentacin, rechazo y
rebose).

( )
=
i i
H H H
f f % F datos
H CH H
r f = R
=
i i i
H H H
f f r
=

i i i
H H H
r y f
=

i i i
H H H
p f r

R
CH

Y
H

H
r
H
r
r

H
f

H
f

H
p

- 141 -

Los valores de las tasas de clasificacin salen de la ecuacin (2.2-87)
i
H
H
50
i
x
0.693
d
y 1 e
| |
|
|
\ .
o
=
Se calculan los flujos en el rebose y el rechazo para los tamaos utilizados en la medida
de distribucin granulomtrica de la batera de hidrociclones 17A. Luego se calculan los
flujos acumulados por malla y se vara el parmetro de Plitt o
H
de modo que el error
cuadrtico medio entre el valor calculdo y el medido, para cada tamao, sea mnimo.
El resultado obtenido al minimizar los errores cuadrticos medios se muestra ms abajo, y
es obtenido para un parmetro de Plitt igual a:
o
H
= 0.378
Tamao Reboses y rechazos acumulados: Modelo v/s datos concentradora A2
[m]

A
H
p
Modelo
A
H
p
Dato
A

A
H
r
Modelo
A
H
r
Dato
A

6700 0,01 0 4,79 34,2 86%
4750 0,06 0 14,34 75,6 81%
3350 0,25 0 35,75 133,2 73%
2360 0,79 0 71,21 210,6 66%
1700 1,83 0 115,77 288 60%
1180 5,12 0 208,48 442,8 53%
850 10,47 6,6 -59% 315,93 619,2 49%
600 20,17 27,6 27% 457,43 851,4 46%
425 34,98 81 57% 620,22 1076,4 42%
300 51,41 147 65% 759,79 1236,6 39%
212 68,75 207 67% 876,85 1357,2 35%
150 83,91 261 68% 960,09 1438,2 33%
106 97,75 307,2 68% 1023,05 1495,8 32%
75 112,14 349,2 68% 1078,26 1540,8 30%
53 123,89 381 67% 1116,91 1573,2 29%
45 130,37 398,4 67% 1136,83 1589,4 28%
23 463,55 600 23% 1936,45 1800 -8%
Error cuadrtico medio 59,95% 50,80%
Tabla 3.3 - 8 Ajuste de parmetros hidrociclones: Reboses y rechazos acumulados: Modelo
v/s datos concentradora A2

- 142 -

Figura 3.3-6 Ajuste de parmetros - Hidrociclones: Rebose acumulado: modelo v/s datos
concentradora A2 Chuquicamata

Figura 3.3-7 Ajuste de parmetros - Hidrociclones: Rechazo acumulado: modelo v/s datos
concentradora A2 Chuquicamata

- 143 -

Con los resultados anteriores, los flujos para cada tamao representativo utilizado en el
simulador, son:
Tamao Flujos en [T/h]
Descripcin [m] H
f
H
H
f
H
y
H
r
H
p
H
r
X
1
203.200 0,00 0,00 0,00 1,00 0,00 0,00 0,00
X
2
124.434 0,00 0,00 0,00 1,00 0,00 0,00 0,00
X
3
73.616 0,00 0,00 0,00 1,00 0,00 0,00 0,00
X
4
31.877 0,00 0,00 0,00 1,00 0,00 0,00 0,00
X
5
15.575 0,00 0,00 0,00 1,00 0,00 0,00 0,00
X
6
6.523 14,40 6,47 7,93 0,99 7,89 0,04 14,36
X
7
2.386 57,60 25,89 31,71 0,97 30,86 0,85 56,75
X
8
1.202 141,60 63,65 77,95 0,94 73,18 4,76 136,84
X
9
601 264,00 118,67 145,33 0,88 128,42 16,91 247,09
X
10
300 333,60 149,96 183,64 0,81 148,53 35,11 298,49
X
11
150 232,80 104,65 128,15 0,72 92,26 35,89 196,91
X
12
53 1.356,00 609,55 746,45 0,58 430,35 316,10 1.039,90
TOTAL = 24.00,00 1.078,85 1.321,15 911,49 409,66 1.990,34
Tabla 3.3 - 9 Ajuste de parmetros hidrociclones: flujos msicos determinados
3.3.7 Parmetros de los Molinos de Bolas
Los parmetros a identificar en el molino de Bolas son los siguientes
- o
M
: Tasa de descarga, ver frmula (2.2-96)
- B: Matriz de tasas de molienda, ver frmula (2.2-94) y (2.2-93)
Las dimensiones fsicas del molino SAG utilizado son:
- Dimetro: d
M
= 18 (5.49 [m])
- Largo: l
M
= 26 (7.92 [m])

- 144 -

El punto de operacin utilizado est determinado por los siguientes valores obtenidos del
muestreo de alimentacin al molino de Bolas 17A de Chuquicamata realizada el 19/06/96:
- Alimentacin total: F
M
= 1249.6 [t/h]
- Densidad de la pulpa:
t M
= 1.92 [T/m
3
]
- Porcentaje de slidos: c
p M
= 76.1%
- Distribucin granulomtrica:
Alimentacin f
M
Descarga p
M

X[mm] Tyler % Retenido
% Acumulado
retenido
% Acumulado
pasante
% Retenido
% Acumulado
retenido
% Acumulado
pasante
6700 3 1.90% 1.90% 98.10% 0.30% 0.30% 99.70%
4750 4 2.30% 4.20% 95.80% 0.70% 1.00% 99.00%
3350 6 3.20% 7.40% 92.60% 1.40% 2.40% 97.60%
2360 8 4.30% 11.70% 88.30% 2.20% 4.60% 95.40%
1700 10 4.30% 16.00% 84.00% 2.80% 7.40% 92.60%
1180 14 8.60% 24.60% 75.40% 6.10% 13.50% 86.50%
850 20 9.80% 34.40% 65.60% 7.10% 20.60% 79.40%
600 28 12.90% 47.30% 52.70% 9.50% 30.10% 69.90%
425 35 12.50% 59.80% 40.20% 11.10% 41.20% 58.80%
300 48 8.90% 68.70% 31.30% 9.70% 50.90% 49.10%
212 65 6.70% 75.40% 24.60% 8.50% 59.40% 40.60%
150 100 4.50% 79.90% 20.10% 6.10% 65.50% 34.50%
106 150 3.20% 83.10% 16.90% 4.80% 70.30% 29.70%
75 200 2.50% 85.60% 14.40% 4.30% 74.60% 25.40%
53 270 1.80% 87.40% 12.60% 3.20% 77.80% 22.20%
45 325 0.90% 88.30% 11.70% 1.70% 79.50% 20.50%
23 -325 11.70% 100.00% 0.00% 20.50% 100.00% 0.00%
Tabla 3.3 - 10 Ajuste de parmetros Molino de bolas: Distribucin granulomtrica
molino de bolas 17A concentradora A2 de Chuquicamata

- 145 -

Figura 3.3-8 Ajuste de parmetros Molino de bolas: Datos de distribucin granulomtrica
molino de bolas 17A concentradora A2
ii) Determinacin de la tasa de descarga o
M

De acuerdo a (2.2-96) la tasa de descarga es la misma para todos los tamaos de
partculas, y est determinado por:
= o
M M M
i i
p w i
Determinaremos el holdup total de slidos (W
M
) del molino de Bolas de la siguiente
manera:

- 146 -

Se asume que el volumen interno del molino es igual al 95% del volumen geomtrico de
este y que el volumen que ocupa la pulpa de mineral en el punto de equilibrio analizado
es de un 40%.
| |
| |
|
= t =
|
|
\ .
\ .
2
M 3
M M
d
V 0.95 l 178 [m ]
2

= =
3
Pulpa M
V 0.95 V 71.2 [m ]
Luego con el valor de densidad de la pulpa y el porcentaje de slidos se determina el
holdup total de slidos.
= =
Pulpa t M M
W V 136.7 [T]
= =
M p M pulpa
W c W 104.0 [T]
Luego:
= o
M M M
i i
p w i
= o
M M M
i i
p w
o = = o =
M M
M M
M M
P F
12.01
W W

- 147 -

iii) Determinacin de la matriz B de tasas de molienda
La variacin del holdup de slidos en el molino est dada por (2.2-93)
=
= +
i 1
M i,i M i,k M
k 1
M M
i i i i k
d
w f p b w b w
dt

En equilibrio se tiene que:
=
= +
i 1
i,i M i,k M
k 1
M M
i i i k
0 f p b w b w (3.3-1)
= =
=

i 1 n
i,i M i,k M
k 1 i 1
M M
i i i k
f p b w b w
= = =
=

n n i 1
M M i,i M i,k M
i 1 i 2 k 1
i k
F P b w b w
Y como se est en equilibrio el flujo total de entrada es igual al total de salida, luego
= = =
=

n n i 1
i,i M i,k M
i 1 i 2 k 1
i k
0 b w b w
= = =
=

n n i 1
i,i M i,k M
i 1 i 2 k 1
i k
b w b w

= = = = = =
= + + + + + +

n 1 2 3 i 1 n 1
i,i M 2,k M 3,k M 4,k M i,k M n,k M
i 1 k 1 k 1 k 1 k 1 k 1
i k k k k k
b w b w b w b w b w b w

- 148 -

( ) ( ) ( )
( )
1 1 2 1 2
1 2 n 1
n
i,i M 2,1 M 3,1 M 3,2 M 4,1 M 4,2 M 4,3 M
i 1
n,1 M n,2 M n,n 1 M
i 3
b w b w b w b w b w b w b w
b w b w b w

=
= + + + + + +
+ + + +

( ) ( )
1 2
k n 1
n
i,i M M 2,1 3,1 4,1 n,1 M 3,2 4,2 5,2 n,2
i 1
n
M i,k M n,n 1
i k 1
i
b w w b b b b w b b b b
w b w b
= +
= + + + + + + + + + +
| |
+ + +
|
\ .

= = = = = + =
= + + + + + +

1 n 1
n n n n n n
i,i M M i,1 M i,2 M i,3 M i,k M i,n 1
i 1 i 2 i 3 i 4 i k 1 i n
i 2 3 k
b w w b w b w b w b w b
Por lo tanto:
= +
=

n
i,i k,i
k i 1
b b (3.3-2)
Es decir, cada trmino de la diagonal de la matriz de molienda es igual a la suma de los
trminos bajo l.
La matriz de molienda (ecuacin (2.2-95)) es:
(
(
(
( =
(
(
(

1,1
2,1 2,2
3,1 3,2 3,3
n,1 n,2 n,3
b 0 0 0
b b 0 0
b b b B
0
b b b 0

- 149 -

Basados en (3.3-2), ocuparemos primero el siguiente caso especial:

(
(
(
(
=
(
(
(
(
(

1,1
1,1 2,2
2,2 3,3
3,3
n 1,n 1
n 1,n 1
b 0 0 0 0
b b 0 0 0
0 b b 0 0
B
0 0 b
b 0
0 0 0 0 b 0

Luego:
= +
M M M
0 f p B w

( (
(
(

(
(
(
(
(
(
=
(
(
(
(
(
(
(
(

(

1
1 1
2
2 2
i i
n n n
1,1 M
2,2 M 1,1 M
i,i M i 1,i 1 M
M n 1,n 1 M
M M
M M
1
M M
i i 1
M M
n 1
b w
f p
b w b w
f p
b w b w f p
f p 0 w b w

Aqu se ve que:
=
1 1
1
M M
1,1
M
f p
b
w

( ) ( )
+
+
= =
2 2 1 1
2 2 1
2 2
M M M M M M 1,1 M
2,2
M M
f p f p
f p b w
b
w w

( ) ( ) ( )
+ +
+
= =
3 3 2 2 1 1
3 3 2
3 3
M M M M M M
M M 2,2 M
3,3
M M
f p f p f p
f p b w
b
w w

- 150 -

En general:
=

=

k k
i
i
M M
k 1
i,i
M
f p
b
w
(3.3-3)
Usando la informacin del molino de Bolas 17A de Chuquicamata y la constante de
descarga o
M
determinada se obtendrn los valores de la diagonal de la matriz de
molienda.
Se determinarn primero los flujos de slidos de entrada y salida del molino de bolas para
los tamaos representativos utilizados en el simulador, en base a los datos de distribucin
y granulomtrica y flujo total del molino de bolas 17A.
El valor del flujo de slidos para un tamao representativo X
i
se extrapola linealmente de
los valores medidos del tamao inmediatamente superior e inferior. Considerando los
valores acumulados pasante, se tiene:
Datos Molino Bolas 17A Interpolacin para tamaos utilizados en el simulador
% acumulado pasante % acumulado pasante Distr. granulomtrica
Tamao
X
i
[m] Alimentacin Descarga
Intervalo

[m]
Tamao
X
i
[m] Alimentacin Descarga f
M
[%] p
M
[%]
6.700 98,1% 99,7% [6700, [ 15.575 100,0% 100,0% 2,1% 0,4%
4.750 95,8% 99,0% [4750, 6700[ 6.523 97,9% 99,6% 9,5% 4,2%
3.350 92,6% 97,6% [3350, 4750[
2.360 88,3% 95,4% [2360, 3350[ 2.386 88,4% 95,5% 12,6% 8,7%
1.700 84,0% 92,6% [1700, 2360[
1.180 75,4% 86,5% [1180, 1700[ 1.202 75,8% 86,8% 23,0% 16,8%
850 65,6% 79,4% [850, 1180[
600 52,7% 69,9% [600, 850[ 601 52,8% 69,9% 21,4% 20,8%
425 40,2% 58,8% [425, 600[
300 31,3% 49,1% [300, 425[ 300 31,3% 49,1% 11,2% 14,6%
212 24,6% 40,6% [212, 300[
150 20,1% 34,5% [150, 212[
106 16,9% 29,7% [106, 150[ 150 20,1% 34,5% 7,5% 12,3%
75 14,4% 25,4% [75, 106[
53 12,6% 22,2% [53, 75[ 53 12,6% 22,2% 12,6% 22,2%
45 11,7% 20,5% [45, 53[
23 0,0% 0,0% [23, 45[
Tabla 3.3 - 11 Ajuste de parmetros Molino de bolas: Determinacin de flujos de slidos como
porcentaje del flujo total

- 151 -

Luego, los flujos de slidos para los tamaos representativos utilizados en el simulador
son:
Tamao Distrib. granulomtrica Flujo slidos
Descripcin [m] fM [%] p
M
[%] fM [T/h] p
M
[T/h]
X
1
203.200 0,0% 0,0% 0,0 0,0
X
2
124.434 0,0% 0,0% 0,0 0,0
X
3
73.616 0,0% 0,0% 0,0 0,0
X
4
31.877 0,0% 0,0% 0,0 0,0
X
5
15.575 2,1% 0,4% 26,4 4,5
X
6
6.523 9,5% 4,2% 118,4 52,2
X
7
2.386 12,6% 8,7% 158,1 108,7
X
8
1.202 23,0% 16,8% 287,6 210,2
X
9
601 21,4% 20,8% 267,9 260,2
X
10
300 11,2% 14,6% 140,2 182,7
X
11
150 7,5% 12,3% 93,6 153,6
X
12
53 12,6% 22,2% 157,4 277,4
Total = 1.249,6 1.249,6
Tabla 3.3 - 12 Ajuste de parmetros Molino de bolas: Determinacin de flujos msicos de
slidos
Finalmente determinamos los valores de la diagonal de la matriz de molienda:
Tamao
Descripcin [m]
M
i
f
M
i
p
M M
i i
f p
=
o
M
i
M
i
M
p
w i,i
b
X
1
203.200 0,0 0,0 0,0 0,0 0,0
X
2
124.434 0,0 0,0 0,0 0,0 0,0
X
3
73.616 0,0 0,0 0,0 0,0 0,0
X
4
31.877 0,0 0,0 0,0 0,0 0,0
X
5
15.575 26,4 4,5 -21,8 0,4 57,7
X
6
6.523 118,4 52,2 -66,2 4,3 20,3
X
7
2.386 158,1 108,7 -49,4 9,0 15,2
X
8
1.202 287,6 210,2 -77,4 17,5 12,3
X
9
601 267,9 260,2 -7,7 21,7 10,3
X
10
300 140,2 182,7 42,5 15,2 11,8
X
11
150 93,6 153,6 59,9 12,8 9,4
X
12
53 157,4 277,4 120,0 23,1 0,0
Tabla 3.3 - 13 Ajuste de parmetros Molino de bolas: Determinacin valores de la diagonal de
la matriz de molienda

- 152 -

Hasta el momento se han obtenido los valores de la diagonal de la matriz de molienda.
Como se vio, de acuerdo a (3.3-2) se puede hacer una matriz de tasas de molienda
particular en la que cada partcula slo se muele a su tamao inmediatamente inferior.
Si podemos determinar la distribucin en los distintos tamaos al moler una partcula de
tamao Xi, como tenemos los valores de la diagonal de la matriz B y dado (3.3-2) se
obtendr la matriz de molienda. Para determinar esta distribucin se utiliz la teora
estadstica de quiebre primario desarrollada por R. R. Klimpel and L. G. Austin [10]. De
acuerdo a esta teora, el porcentaje de partculas de tamao menor a y producto del
quiebre de una partcula de tamao x es:
| | | |
| | | | | |
= | |
| | |
| |
\ . \ . \ .
\ . \ .
S V
L
r r
r 2 3
y y y
C(x, y) 1 1 1 1
x x x

Donde r
L
, r
S
, y r
V
son constantes que representan las fallas por bordes, reas y volmenes
respectivamente. Los valores de r
L
, r
S
, y r
V
pueden ser fraccionarios mayores a 1 o cero
pero la suma de los tres debe ser a lo menos 1.
Con los valores de C(x,y) se puede construir una matriz de distribucin acumulada de
ruptura, que tiene la siguiente forma:
(
(
(
(
=
(
(
(

2,1
3,1 3,2
n,1 n,2 n,3
1 0 0 0
c 1 0 0
c c 1 C
0
c c c 1

En que:
| | | |
| | | | | |
= | |
| | |
| |
\ . \ . \ .
\ . \ .
S V
L
r r
r 2 3
i,j
Xj Xj Xj
c 1 1 1 1
Xi Xi Xi

- 153 -

Luego se puede obtener una matriz de distribucin de ruptura por tamao:
(
(

(
( | =
(
(
(

2,1
2,1 3,1 3,2
3,1 4,1 3,2 4,2 4,3
n,1 n,2 n,3
1 c 0 0 0
c c 1 c 0 0
c c c c 1 c
0
c c c 1

Y dado los valores b
i,j
y la ecuacin (3.3-2) se obtiene la matriz de molienda:

( (
( (

( (
( ( =
( (
( (
( (

2,1 1,1
2,1 3,1 3,2 2,2
3,1 4,1 3,2 4,2 4,3
n 1,n 1
n,1 n,2 n,3
1 c 0 0 0 b 0 0 0
c c 1 c 0 0 0 b 0 0
c c c c 1 c B 0 0
0 b 0
c c c 1 0 0 0 0

Obtenida la matriz de molienda, se determina el holdup utilizando (3.3-1) y los flujos de
alimentacin y descarga para cada intervalo de tamao determinados con los datos del
molino 17A. Se sencibiliza con los valores de r
L
, r
S
, y r
V
de modo de minimizar el error
cuadrtico medio entre el holdup determinado con la matriz de molienda y el obtenido
directamente ocupando la tasa de descarga o
M
. Los resultados son:
r
L
= 1.02277
r
S
= 0.03462
r
V
= 1.60714

- 154 -

0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 -28,57 0 0 0 0 0 0 0
B = 0 0 0 0 19,23 -8,59 0 0 0 0 0 0
0 0 0 0 4,74 4,59 -9,29 0 0 0 0 0
0 0 0 0 2,31 2,06 5,08 -7,45 0 0 0 0
0 0 0 0 1,15 0,98 2,21 4,08 -6,22 0 0 0
0 0 0 0 0,57 0,48 1,02 1,77 3,40 -7,17 0 0
0 0 0 0 0,37 0,31 0,64 1,05 1,88 4,94 -3,82 0
0 0 0 0 0,20 0,17 0,35 0,56 0,94 2,23 3,82 0

El siguiente grfico (Figura 3.3-9) muestra la comparacin entre los valores del holdup
obtenidos con la matriz B versus los obtenidos con la tasa de descarga.

Figura 3.3-9 Ajuste de parmetros Molino de bolas: Comparacin entre holdup obtenido por
matriz de molienda v/s el obtenido con tasa de descarga

- 155 -

Posteriormente, mediante pruebas con el simulador, se generaron las matrices para tratar
mineral duro y blando (ver captulo 2.5.1) modificando la matriz anterior de modo que, al
haber slo mineral duro la carga circulante sea 350% (matriz B
D
) y de haber slo mineral
blando la carga circulante de sea 250% (matriz B
B
). El simulador trabaja con dos flujos
uno de mineral duro tratado con la matriz B
D
y otro de mineral blando tratado con la matriz
B
B
.

- 156 -

4 CONCLUSIONES Y RECOMENDACIONES
Se realiz en forma exitosa un simulador de una planta de molienda semi-autgena con
circuito inverso de molienda secundaria de bolas. En reuniones de trabajo se revis el
buen funcionamiento tanto de la planta simulada como de los sistemas de control
implementados.
El simulador fue utilizado en un curso terico orientado a operadores de plantas de
molienda reales y en cursos de la Escuela de Ingeniera en Minas. En estos seminarios y
cursos el simulador tuvo una buena evaluacin tanto por su versatilidad como por el
ambiente grfico que permita un intuitivo manejo de este. Adjunto se presenta el Manual
del Usuario que fue entregado en estos cursos.
Adems este simulador fue utilizado como fuente para disear un simulador en
FactoryLink.
Es importante mencionar que la informacin contenida en esta memoria da un compendio
de los principales modelos contenidos en una planta de molienda semi-autogena, adems
de la integracin de estos e implementacin en tiempo discreto.
Por otro lado se explica metodologas para determinar los parmetros de los modelos
utilizando la informacin tpica que se obtiene de plantas reales. Esto permite el
determinar los parmetros de equipos de caractersticas diferentes (tamao, modelos,
etc), ingresarlos en el simulador va un ambiente grfico y estudiar el comportamiento de
la planta con esta modificacin as como el funcionamiento de los sistemas de control.
Las siguientes son algunas recomendaciones para uso y complementacin del simulador:
- Agregar otros sistemas de control como control predictivo, controles con lgica
difusa, etc.

- 157 -

- Utilizar este simulador en cursos tanto de Ingeniera en Minas como de Sistemas
de Control.
- Integrar este simulador con simuladores del proceso aguas arriba: Chancado
Primario y transporte de mineral; y aguas abajo: plantas de concentracin (Cobre:
flotacin; Oro/Plata: lixiviacin, CCD; etc.)

- 158 -

REFERENCIAS
[1] Weymont, N.P., 1979. The Analysis and Simulation of Autogenous Grinding Systems.
Ph.D. thesis, The Pennsylvania State University, the Graduate School, Department of
Mineral Engineering.
[2] Austin, L.G., Menacho, J. y Pearcy, F., 1987. A General Model for Semiautogenous
Milling. ACOM 87 Proceeding of Twentieth International Symposium on the Application of
Computers and Mathematics in the Mineral Industries, Vol. 2, Metalurgy, pp 107 126.
AIME, Johannesburg.
[3] Barahona, C.A. 1984. Modeling and Simulation of Semiautogenous Grinding Systems.
Master of Science Thesis, Pennsylvania State University.
[4] Austin, L.G., 1990. A mill power equation for SAG mills, Mineral and Metallurgical
Processing, 7:(1),57-62
[5] A Casali. Simplificacin del Modelo de Chancado de Andersen / Whiten (ver ANEXO A)
[6] A. Lynch and T. C. Rao, Modeling and scale-up of hydrocyclone classifiers, Proc. 11
th

International Mineral Processing Congress, Cagliari, Italy. 4-25.
[7] L. R. Plitt, A Mathematical Model of Hydrocyclone Classifier CIM Bulletin December
1976 p. 114.
[8] Gonzlez, G. Informe Final del Proyecto CONICYT N124: Control ptimo de una planta
concentradora de minerales de cobre. Departamento de Electricidad, U. de Chile, Sept.
1971
[9] Gonzlez, G.D. "Adaptive Optimal Equilibrium Control for an Ore Concentrator Model".
Tesis de Ph.D., University of Michigan, U.S.A., 1981.

- 159 -

[10] R. R. Klimpel and L. G. Austin. The statistical theory of primary breakage distributions
for brittle materials. Society of Mining Engineers, Transactions, March 1965.
[11] J. M. Link, E. J. Wasp, C. A. Horne. Hydraulics SME Mineral Processing Handbook,
N. L. Weiss editor, 1985.
[12] R. Amstica, G. Gonzlez, J. Barra , L. Mange J. Menacho, and O. Castro (1993). "A
SAG mill circuit dynamic simulator based on a simplified mechanistic model". in Proc.
XVIII International Mineral Processing Congress, Sydney, Australia, 1993, vol. 1, pp.
117-130

ANEXO A
SIMPLIFICACIN DEL MODELO DE CHANCADO DE ANDERSEN / WHITEN

- 161 -

MODELO DE CHANCADO

-Modelo de Andersen/whiten

-Explicado, simplificado y procesado por A Casali

1.- Descripcin del modelo

Ecuaciones de balance:
X = f + B C x (1)

X = p + C x (2)

Donde:
f es el vector columna con los flujos de alimentacin por
tamaos.

P es el vector columna con los flujos de producto por
tamaos

C es la matriz diagonal de clasificacin, con elementos
en la diagonal Ci correspondiente a la clasificacin
para facturara de cada tamao i.

B es la matriz triangular inferior de ruptura, con
elementos bij correspondientes a la fraccin de mineral
de tamao j que al chancarse termina de tamao i.

CLASIFICACION
RUPTURA
f x
Cx
BCx
P

- 162 -

Combinando las ecuaciones (1) y (2) se obtiene el modelo de Whiten:

2.- Matriz de clasificacin

Elementos de la matriz C:

A) Ci = 1 si di > K2

b) Ci = 0 si di < K1

c) Para K1 < di < K2

2 ,3
K 2 - d i
C i = 1 -( )
K 2 - k 1

Donde:

di [mm], es la media Geom. del intervalo de tamao i.

Ecuaciones para K1 y K2

K1 = A0 + A1*Sc + A2*Gs + A3*F80
K2 = A4 + A5*Sc A6*Gs + A7*F80

Donde:

Sc es el setting cerrado del chancador [mm]
(12,7 mm, para el caso del simulador)

Gs es el flujo mineral alimentado [t/h]

F80 es el tamao bajo el cual est el 80% de la
alimentacin al chancador [mm]

Valores iniciales:

p = (I C) * (I B C)
-1
* f (3)
Parmetros : A0 A1 A2 A3 A4 A5 A6 A7
Valores iniciales : 0,1 1 0 0 0,1 2 0 0

- 163 -

3.- Matriz de ruptura

Elementos de la matriz B: bij

a) 1a. columna:

bi1 = 0 si i = 1

bi1 = R(i-1)1 Ri1 si i N (N = ultima fraccin)

bN1 = 1 - E(bkj)(j-1) si i > j con i N (*)

b) resto de la matriz (otra columna)

bij = 0 si i j

bij = b(i-1)(j-1) si i > j, con i N

bNj = 1 -E(bkj) , con k N

(*) NOTA: cada elemento de la ltima fila de la matriz B
corresponde siempre a 1 la suma del resto de su
columna respectiva.

Modelo para el ri1

A 8 K 3
d i
R i 1 = 1 -( 1 -( ) )
d i

Ecuacin para K3

K3 = A9 A10*SC A11*Gs A12*F80

Valores iniciales:

Parmetros : A8 A9 A10 A11 A12
Valores iniciales : 0,9 0,03 1,2 0,05 0,01

ANEXO B
SIMULADOR MOLIENDA SAG - MANUAL DEL USUARIO

M
M
M
a
a
a
n
n
n
u
u
u
a
a
a
l
l
l

d
d
d
e
e
e
l
l
l

U
U
U
s
s
s
u
u
u
a
a
a
r
r
r
i
i
i
o
o
o

D
D
D
S
S
S
A
A
A
G
G
G

Ctedra Phelps Dodge de Procesamiento de Minerales
Proyecto DOCENMIN

MAYO DE 1999

Universidad de Chile
Departamento de Ingeniera de Minas
Departamento de Ingeniera Elctrica
Simulador DSAG Manual del Usuario

Pgina i
Resumen
El presente manual explica como utilizar el simulador DSAG, que es uno de los
simuladores desarrollados en el proyecto DOCENMIN, destinado a mejorar la
docencia en control y simulacin de plantas de procesamiento de minerales, en
el marco de la ctedra Phelps-Dodge de Procesamiento de Minerales, del
Departamento de Ingeniera en Minas de la Facultad de Ciencias Fsicas y
Matemticas de la Universidad de Chile.
El simulador DSAG implementa una planta de molienda semiautgena que
consta de un molino SAG con harnero y chancador de pebbles, y dos circuitos
inversos de molienda secundaria, cada uno con una bomba centrfuga de
velocidad variable, una batera de hidrociclones y un molino de bolas.
El simulador DSAG fue desarrollado en Matlab 4.1 y fue programado de forma
de crear un entorno grfico como el que existe en una planta real y simular las
distintas condiciones de operacin que pueden encontrarse en estas plantas.

Pgina iii
Indice
Resumen............................................................................................................... i
Indice..................................................................................................................iii
1 Introduccin................................................................................................... 1
2 Para comenzar................................................................................................ 2
3 Descripcin general ....................................................................................... 4
3.1 Pantalla SAG.................................................................................................... 5
3.1.1 Grficos ................................................................................................................. 5
3.1.2 Controles ............................................................................................................... 6
3.1.3 Subpantalla............................................................................................................ 6
3.1.4 Comandos del simulador ....................................................................................... 6
3.1.5 Alarma de sobrecarga y produccin ...................................................................... 7
3.1.6 Botones SAG, SEC y Layout ................................................................................ 7
3.2 Pantalla SEC .................................................................................................... 8
3.2.1 Grficos ................................................................................................................. 8
3.2.2 Controles ............................................................................................................... 9
3.2.3 Subpantalla............................................................................................................ 9
3.3 Pantalla Layout .............................................................................................. 10
4 Sesin de simulacin.................................................................................... 12
5 Controles Automticos................................................................................. 14
5.1 Cambio de parmetros .................................................................................. 14
5.2 Lazos de control ............................................................................................. 16
5.2.1 Control PI de hold-up.......................................................................................... 16
5.2.2 Controles PID en molienda secundaria ............................................................... 17
6 Guardando la sesin .................................................................................... 20
6.1 Recuperando la sesin desde EXCEL.......................................................... 21
7 Configuracin .............................................................................................. 24
8 ndice de Ilustraciones................................................................................. 27
9 Referencias................................................................................................... 28


Pgina 1
1 Introduccin
El simulador DSAG simula el comportamiento de una planta de molienda que
incluye molienda semiautgena y convencional. Para esto se realiz la unin de
diversos modelos matemticos. El modelo del molino SAG utilizado es el
desarrollado por Amstica, Gonzlez et al. [1][2]. El modelo del harnero es el
de Karra [4] adaptado por A. Casali, 1997. El modelo de chancador es el de
Andersen y Whiten [5][6] adaptado por A. Casali, 1997. El modelo de
hidrocicln es el de Lynch y Rao [7] y modificado por Plitt [8]. El modelo de
molino de bolas es una adaptacin del modelo utilizado por Gonzlez [9].
El presente manual indica cmo trabajar con el simulador, ya sea desde el punto
de vista de supervisin o desde el punto de vista de control automtico.
Para un estudio ms detallado del simulador ver Manual Tcnico de este
simulador [10] o ver Orellana [3].
Manual del Usuario Simulador DSAG

Pgina 2
2 Para comenzar
El simulador DSAG fue desarrollado en Matlab 4.1 que es un lenguaje de
programacin poderoso matemticamente. El simulador requiere de Matlab
versin 4.1 o superior para ser ejecutado.
Para iniciar una sesin con el simulador cargue Matlab. Luego, desde la ventana
de comando vaya al directorio donde se encuentran los archivos del simulador.
Por ejemplo, si instal Matlab en forma estndar, el editor Matlab comienza en
el directorio c:\matlab\bin y si el simulador se encuentra en c:\dsag escriba
cd \dsag para cambiarse a este directorio (el comando cd cambia al directorio
especificado; cd por si slo muestra el directorio actual; cd .. retrocede hacia el
directorio inmediatamente superior).
Una vez que est en el directorio del simulador escriba: dsag

Figura 1. Partiendo el simulador desde Matlab
Esto iniciar el simulador. Al partir se presenta la pantalla de supervisin del
molino SAG (Figura 2).

Pgina 3

Figura 2. Pantalla inicial al partir el simulador
El simulador consta de tres ventanas:
- Supervisin y control del molino sag.
- Supervisin y control de la molienda secundaria
- Layout de la planta y valor de los flujos de mineral y de las principales
variables.


Pgina 4
3 Descripcin general
El circuito de molienda simulado es una planta que consta de una etapa de
molienda semiautgena y de una etapa con molienda secundaria de bolas en
circuito inverso. Tanto para definir la forma del simulador como para el ajuste
de sus parmetros, se utiliz informacin de varias plantas concentradoras de
minerales de cobre.
En la figura N3 se observa el layout del circuito de molienda que cuenta con
un molino SAG, un chancador de pebbles, dos molinos de bolas y dos bateras
de hidrociclones.

Figura 3. Layout del circuito de molienda

Pgina 5
3.1 Pantalla SAG
La primera pantalla que muestra el simulador es la de control y supervisin del
molino SAG.
Tanto esta pantalla como la de la molienda secundaria estn constituidas por
dos monitores con grficos de tendencia, una zona para manejar las variables
manipulables, una subpantalla con informacin de la otra etapa de molienda y
una zona con los sistemas de control del simulador.
La pantalla de control y supervisin del SAG se muestra en la figura 4.

Figura 4. Pantalla de molienda SAG
3.1.1 Grficos
El grfico superior muestra la potencia consumida por el molino SAG en [KW].
La potencia se grafica en rojo. El valor numrico instantneo de la potencia se
da, tambin en rojo, justo sobre el grfico.
El grfico inferior muestra en azul el hold-up de mineral dentro del molino en
[T], en violeta la alimentacin fresca de slidos en [T/h] y en verde el
set-point de hold-up del control automtico de hold-up del molino SAG. El
hold-up y el set-point de hold-up utilizan la escala de la izquierda que est en

Pgina 6
[T], la alimentacin de slidos utiliza la escala en violeta que est al lado
derecho y que indica [T/h].
3.1.2 Controles
Al costado derecho en color caqui tenemos los controles del SAG.
Lo primero que se observa es un slider para controlar la alimentacin fresca de
slidos al SAG.
Despus tenemos un slider para controlar el agua en la alimentacin al SAG y
una casilla para fijar o no una razn fija de 30% de agua en la alimentacin.
Ms abajo hay una casilla para activar el control PI de hold-up, un slider para
fijar el set-point de hold-up y en azul el valor actual del hold-up en [T].
Por ltimo hay una casilla para activar el chancador de pebbles.
3.1.3 Subpantalla
Sobre la zona de controles, se encuentra una subpantalla en color verde agua
con datos importantes de la molienda secundaria.
En la subpantalla se muestra el valor numrico instantneo de:
- El +65# del producto final o rebose de los hidrociclones.
- El porcentaje de slidos (Cph [%]) en la alimentacin de los hidrociclones.
- La carga circulante (C.C. [%]).
3.1.4 Comandos del simulador
Debajo de los grficos estn los comandos del simulador.
A la izquierda tenemos los botones para comenzar y para poner en pausa la
simulacin, mas a la derecha hay un slider para variar la velocidad de la
simulacin.
A continuacin hay un botn para realizar un refresco de los grficos, en caso
de que sea necesario.

Pgina 7
Despus tenemos un men emergente para elegir el intervalo de tiempo
mostrado en el grfico, lo que es til cuando se est en pausa y se desea ver la
tendencia 4 horas hacia atrs o ms.
Por ltimo, se muestra el tiempo de simulacin en minutos de planta.
3.1.5 Alarma de sobrecarga y produccin
En el extremo superior izquierdo encontramos una alarma que se activa cuando
el molino entra en sobrecarga.
En la parte superior se indica, en negro, la produccin total en toneladas del
circuito de molienda desde que se inici la simulacin hasta el tiempo indicado
abajo.
3.1.6 Botones SAG, SEC y Layout
Sobre la subpantalla Molienda Secundaria se encuentran tres botones que nos
llevan a las ventanas de molienda SAG (desactivada en este momento pues ya
estamos en esa ventana), a la molienda secundaria (botn SEC) y al layout de la
planta (botn Layout).

Pgina 8
3.2 Pantalla SEC
Apretando el botn SEC llegamos a esta pantalla la que podemos observar en la
figura 5.

Figura 5. Pantalla de molienda secundaria (SEC)
3.2.1 Grficos
En el grfico superior se presenta en rojo el porcentaje de slidos en la
alimentacin de los hidrociclones, en verde el set-point del control PID de este
porcentaje de slidos, y en violeta el agua que se agrega al pozo, que es una
variable manipulada del control.
El grfico inferior muestra en azul la tendencia del +65# en el producto final
(rebose de los hidrociclones) y en verde el set-point del control PID del +65#.
Este control PID es un control en cascada en que la variable manipulada es el
set-point del control PID de densidad de la alimentacin a los hidrociclones.
Sobre los grficos se puede observar la alarma de sobre carga en el molino
SAG, la produccin total y la carga circulante en el circuito de molienda
secundaria.

Pgina 9
3.2.2 Controles
En el panel de controles tenemos un slider para fijar el agua que se agrega al
pozo. Cuando no hay controles conectados, el valor del agua al pozo se muestra
en color violeta al igual que en el grfico.
Despus puede observarse la casilla para activar el control de porcentaje de
slidos en la alimentacin a hidrociclones (Cph), el slider para fijar el set-point
y el valor actual del Cph.
Sigue el control PID en cascada de +65# a la salida de los hidrociclones, al
igual que el anterior tenemos el slider para fijar el set-point y el valor puntual
del +65#.
Por ltimo aparece el control PID del pozo, este es un control que est activado
en forma predeterminada y que controla el nivel del pozo manipulando la
velocidad de la bomba centrfuga. Aqu adems de estar el slider para el set-
point de nivel del pozo estn los valores puntuales del nivel y de la velocidad de
la bomba.
3.2.3 Subpantalla
En la subpantalla que est sobre los controles, se presenta ahora valores
importantes de la molienda SAG: el valor de la potencia, del hold-up del molino
y el estado del chancador (conectado/desconectado).


Pgina 10
3.3 Pantalla Layout
Apretando el botn layout llegamos a la pantalla de layout (figura 6).

Figura 6. Pantalla de Layout
En esta pantalla se presenta el valor actual de las principales variables del
simulador.
Siguiendo el layout tenemos:
- F: flujo msico de slidos en [T/h] en la alimentacin del SAG.
- Fa: flujo msico de agua en [T/h] en la alimentacin del SAG.
- Hold-up: el hold-up del molino SAG en [T].
- Pot: potencia consumida por el motor del Molino SAG en [KW].
- P: flujo msico de slidos en [T/h] en la descarga del harnero.
- Pa: flujo msico de agua en [T/h] en la descarga del harnero.
- R: flujo msico de slidos en [T/h] en el rechazo del harnero.
- Ra: flujo msico de agua en [T/h] en el rechazo del harnero.
- hp: nivel del pozo en [m].
- G: flujo msico de agua en [T/h] que es adicionada al pozo.
- V: velocidad de la bomba centrfuga en [RPM].
- Fc: flujo msico de slidos en [T/h] en la alimentacin a hidrociclones.

Pgina 11
- Fac: flujo msico de agua en [T/h] en la alimentacin a hidrociclones.
- Cph: porcentaje de slidos en [%]en la alimentacin a hidrociclones.
- Pc: flujo msico de slidos en [T/h] en rebose de hidrociclones.
- Pac: flujo msico de agua en [T/h] en rebose de hidrociclones.
- +65#: +65# en [%] en rebose de hidrociclones.
- Fb: flujo msico de slidos en [T/h] en alimentacin del molino de bolas.
- Fab: flujo msico de agua en [T/h] en alimentacin del molino de bolas.
- C.C.: Carga circulante en [%].
- Pb: flujo msico de slidos en [T/h] en descarga del molino de bolas.
- Pab: flujo msico de agua en [T/h] en descarga del molino de bolas.
En el extremo inferior izquierdo adems de los botones para continuar y poner
en pausa la simulacin, hay una casilla para utilizar o no el chancador de
pebbles.
El circuito de molienda secundaria tiene un pozo y dos circuitos paralelos
formados cada uno por una bomba centrfuga, una batera de hidrociclones y un
molino de bolas. Por razones visuales en el layout de la planta se ha
representado un circuito equivalente a los dos circuitos reales, tal como se
indica en el extremo inferior derecho de la pantalla. Por lo tanto, los flujos
mostrados son la suma de los flujos de los dos circuitos.

Pgina 12
4 Sesin de simulacin
Para iniciar una nueva sesin presione el botn Comenzar. La simulacin parte
desde un estado de equilibrio. Este estado inicial slo puede ser cambiado en
forma externa (ver Configuracin en la pag. 25).
Los controles del simulador estn debajo de los grficos. Con el slider
velocidad simulacin se har ms rpida o ms lenta la simulacin. Al apretar
Pausa, el simulador se detiene y se podr observar el estado de las diferentes
variables al pasar de una ventana a otra con los botones SAG, SEC y Layout.
Para ver datos que quedaron fuera del intervalo de tiempo del grfico actual,
utilice el men emergente Interv. tiempo. Los intervalos de tiempo mostrados
en los grficos son de 4 horas, es decir de 0 a 4 hr, de 4 a 8, etc. hasta llegar a
20 a 24 hr.
Puede utilizar zoom para explorar ms en detalle los grficos. Ponga el mouse
sobre el punto donde desea hacer zoom, presione el botn izquierdo del mouse
y el eje del tiempo se expandir al doble centrado en el punto donde puso el
mouse. Para realizar un zoom out siga el procedimiento anterior pero
apretando el botn derecho del mouse.
Figura 7. Realizando Zoom
Estando en Pausa si uno se cambia a otros programas, puede ocurrir que al
volver al simulador no aparezca parte o todo el grfico, si esto ocurriese,
presione Redibuja, esto realizar una actualizacin del grfico. Si presiona
Redibuja cuando se encontraba el grfico en zoom, se volver al estado
original sin zoom.

Pgina 13
Terminada las primeras 24 hr de simulacin el simulador quedar en pausa y le
preguntar si desea guardar la sesin. Despus de responder, si presiona
Continuar seguir la simulacin, pero no podr recorrer los valores anteriores.
Si guard la sesin podr revisarla posteriormente en una planilla electrnica
(Ver Guardando la sesin en la pag. 21).
Para terminar la simulacin elija la opcin salir en el men simulador.

Pgina 14
5 Controles Automticos
El simulador cuenta con cuatro controles automticos:
- Control PI de hold-up del Molino SAG el que acta sobre la alimentacin
fresca al molino.
- Control PID de nivel del pozo que acta sobre la velocidad de la bomba
centrfuga del pozo.
- Control PID de porcentaje de slidos en la alimentacin a hidrociclones que
acta sobre el agua que se agrega al pozo.
- Control PID de +65# en el producto final que acta en cascada con el
control anterior variando el set-point del control de Cph.
Para activar los controles se debe seleccionar la casilla correspondiente en el
panel de control del simulador.
5.1 Cambio de parmetros
Para cambiar los parmetros de los controles, en el men Configurar elija la
opcin PID (figura 8).

Figura 8. Llamando a la ventana de parmetros de los PID

Pgina 15
Esto desplegar una ventana (figura 9) donde Ud. puede editar y cambiar el
valor de los parmetros de los controladores.
Para cambiar el valor pngase con el mouse en la casilla de edicin
correspondiente y apriete el botn izquierdo, as la editar. Escriba el valor
deseado. Para que el simulador asuma el nuevo valor presione el botn Aceptar.

Figura 9. Parmetros de los PID
El botn Cerrar cierra la ventana, si realiz cambios y no apret Aceptar, al
cerrar la ventana no se realiza ningn cambio en los parmetros.
El botn Guardar, cambia el valor predeterminado del control por el que usted
est especificando, as la prxima vez que utilice el simulador tendr el ltimo
valor guardado.
El botn Restaurar preguntar si desea volver a los ltimos valores guardados
(es decir de la ltima vez que apret Guardar) o los valores predeterminados del
software.

Pgina 16
5.2 Lazos de control
Los lazos de control se resumen en el esquema de la figura 10.
LIC
AIC
LT
AT
Lp*
Cph*
+65#
*
Velocidad Bomba
DT DIC
Set-point de Cph
Agua al pozo
WT WIC
Hold-up*
Alimentacin Fresca

Figura 10. Lazos de control
A continuacin se da una explicacin de cada uno de los controles automticos
y como manejarlos.

5.2.1 Control PI de hold-up
Este es un control PI que controla el hold-up del molino SAG manipulando la
alimentacin de slidos. Por lo tanto, al conectar el control uno pierde la opcin
de variar manualmente el valor de la alimentacin y el slider correspondiente
desaparecer del panel de control.

Pgina 17

Figura 11. Trabajando con el control PI de hold-up
El control implementado es la versin discreta incremental del controlador PI:
|
.
|
\
|
+
A

+ =
Ti
) t ( e
t
) 1 t ( e ) t ( e
K ) 1 t ( F ) t ( F
C

Donde F es el flujo de slidos en la alimentacin en [T/h], Kc es la constante
proporcional, Ti el tiempo integral y e(t) el error, dado por
e(t) = Hold-up
*
(t)-Hold-up(t). [Hold-up
*
es el set-point].
5.2.2 Controles PID en molienda secundaria
Estos controles tambin son discretos y se activan en las casillas
correspondientes. Al igual que en el caso del control PI de hold-up, las
variables manipuladas de los PID dejarn de estar disponibles al usuario
mientras est conectado el control (los sliders desaparecen y las casillas de
valor no se pueden acceder).

Pgina 18

Figura 12. Trabajando con los controles PID de la molienda secundaria
Los controles son los siguientes:
Control PID del nivel del pozo (hp)
|
.
|
\
|
A
+
+ +
A

+ =
2
d C
t
) 2 t ( e ) 1 t ( e 2 ) t ( e
T
Ti
) t ( e
t
) 1 t ( e ) t ( e
K ) 1 t ( V ) t ( V
en que:
e(t) = hp
*
(t) hp(t)
V : Velocidad de la bomba en [RMS]
T
d
: tiempo derivativo.
Control PID del porcentaje de slidos en alimentacin de Hidrociclones (Cph)
|
.
|
\
|
A
+
+ +
A

+ =
2
d C
t
) 2 t ( e ) 1 t ( e 2 ) t ( e
T
Ti
) t ( e
t
) 1 t ( e ) t ( e
K ) 1 t ( G ) t ( G
en que:
e(t) = Cph
*
(t) Cph(t)
G : Flujo de agua al pozo en [T/h]

Pgina 19

Control PID del +65# en producto final
|
.
|
\
|
A
+
+ +
A

+ =
- -
2
d C
t
) 2 t ( e ) 1 t ( e 2 ) t ( e
T
Ti
) t ( e
t
) 1 t ( e ) t ( e
K ) 1 t ( Cph ) t ( Cph
en que:
e(t) = +65#
*
(t) +65#(t)
Cph
*
: es el set-point del control PID de Cph
Este control de granulometra en el producto final es un control en cascada con
el control PID de porcentaje de slidos en la alimentacin a hidrociclones, es
decir, el set-point del control de porcentaje de slidos es la variable manipulada
para el controlador de +65#.


Pgina 20
6 Guardando la sesin
Para guardar la sesin con el objetivo de algn anlisis posterior, en el men
Simulador elija la opcin Guardar (figura 13).

Figura 13. Men Simulador Guardar
El simulador desplegar una ventana proponiendo un nombre, que ser de la
forma Da_Hora_Minuto (Ej: propondr el nombre 12_16_23 si en el
momento de guardar son las 16:23 del da 12). Sin embargo usted puede poner
el nombre que dese al archivo, siempre y cuando utilice 8 dgitos y smbolos
permitidos. Si coloca un nombre de ms de 8 caracteres, se guardar la sesin
con los primeros 8 caracteres del nombre que dio.
El archivo con la sesin quedar en el directorio donde se corri el simulador y
tendr extensin sag (Ej: 12_16_23.sag). Este archivo tiene formato texto.
Despus de guardar Ud. puede seguir simulando. Puede guardar todas las veces
que dese, la sesin de simulacin.
Para recuperar la informacin de la sesin de simulacin, abra el archivo desde
una planilla electrnica.

Pgina 21
6.1 Recuperando la sesin desde EXCEL
Si corri el simulador desde el directorio DSAG, despus de guardar la sesin
al observar este directorio ver el archivo guardado, en este caso 3_17_32.sag.

Figura 14. Archivo de sesin de simulacin
Inicie EXCEL y en el men Archivo elija abrir, seleccione el archivo de la
sesin de simulacin, EXCEL comenzar una ayuda para importacin de datos.

Figura 15. Rescatando con EXCEL archivo de sesin de simulacin

Pgina 22
Si sigue las opciones propuestas por EXCEL, abrir en forma exitosa el
archivo.
Obtendr una hoja EXCEL como la siguiente:

Figura 16. Archivo de sesin de simulacin abierto desde EXCEL
En la primera fila est el encabezado de cada variable. La nomenclatura de
variables es la siguiente:
- t[h]: Tiempo en horas
- F[T/h]: Flujo de slidos en alimentacin al SAG en [T/h]
- Fa[T/h]: Flujo de agua en alimentacin al SAG en [T/h]
- P[T/h]: Flujo de slidos en descarga del harnero en [T/h]
- Pa[T/h]: Flujo de agua en descarga del harnero en [T/h]
- R[T/h]: Flujo de slidos en rechazo del harnero en [T/h]
- Ra[T/h]: Flujo de agua en rechazo del harnero en [T/h]
- Pot[KW]: Potencia consumida por el motor del SAG en [KW]
- W[T]: Hold-up del SAG en [T]
- PID_W: Estado del control PI del SAG (0 desconectado 1 conectado)
- spW[T]: Set-point del control PI de hold-up en [T]
- +65#_P: +65# del bajo tamao del harnero en [/1]
- Chancado: Estado del chancador de pebbles (0 desconectado 1 conectado)

Pgina 23
- Raz.Agua: Estado del control de razn fija de agua en un 30% en
alimentacin fresca al SAG (0 desconectado 1 conectado)
- hp[m]: Nivel del pozo en [m]
- PID_hp: Estado del control PID de nivel del pozo (0 desconectado 1
conectado)
- sphp[m]: Set-point del control PID de nivel del pozo en [m]
- V[RPM]: Velocidad de la bomba centrfuga del pozo en [RMP]
- G[T/h]: Flujo de agua agregada al pozo en [T/h]
- Cph: Porcentaje de slidos en alimentacin a hidrociclones en [/1]
- PID_Cph: Estado del control PID de porcentaje de slidos en alimentacin a
hidrociclones (0 desconectado 1 conectado)
- spCph: Set-point del control PID de Cph en [/1]
- Fc[T/h]: Flujo de slidos en alimentacin a hidrociclones en [T/h]
- Pc[T/h]: Flujo de slidos en rebose de hidrociclones en [T/h]
- +65#_Pc: +65# en rebose de hidrociclones en [/1]
- PID_+65#: Estado del control PID de +65# en rebose de hidrociclones (0
desconectado 1 conectado)
- sp+65#Pc: Set-point del control PID de +65# en rebose de hidrociclones en
[/1]
- CpPc: Porcentaje de slidos en producto final (rebose de hidrociclones) en
[/1]
- +65#_Pb: +65# en descarga del molino de bolas en [/1]
- C.C.: Carga circulante del circuito de molienda en [%]
- mu: Dureza del mineral (0 = el ms blando, 1 = el ms duro)
Hay que tener cuidado en la configuracin regional que tenga el computador, el
punto debe ser separador de decimales, si no tendr problemas al abrir el
archivo de sesin de simulacin. La configuracin regional se establece en el
Panel de Control de Windows. Para abrir el panel de control seleccione
Inicio, Configuracin, Panel de control.

Pgina 24
7 Configuracin
Para cambiar el valor de las variables del simulador, en el men Configurar
elija la opcin Modelo. Aparecer una ventana preguntando por la contrasea
(password) de acceso. Ingrese la password y presione aceptar.

Figura 17. Ventana de password
La primera ventana que se despliega es la llamada Mineral que trae el set-up
de parmetros y variables relacionadas con el tipo de mineral. Adems aqu es
donde puede Ud. cambiar la password de acceso (para un detalle de cada una de
estas variables vea el Manual Tcnico).

Figura 18. Ventana de Configuracin de Mineral

Pgina 25
Para cambiar un valor, cmbielo en el cuadrado de texto y presione Aceptar.
El botn Cerrar cierra la ventana, el botn Guardar cambia el valor
predeterminado de la variable, el botn Restaurar preguntar si desea
redefinir las variables con los valores de la ltima vez que se guardaron o con
los valores de fabricacin.
El botn ayuda despliega la ventana de ayuda correspondiente.
Si presiona el botn SAG, se desplegar las variables relacionadas con el set-up
del molino SAG, el harnero y el chancador de pebbles.

Figura 19. Ventana de Configuracin del SAG
Los botones inferiores tienen la misma funcin que en la ventana Mineral.
Todos los botones, excepto Cerrar, actan sobre las variables desplegadas.
As, si Ud. cambi un valor en la ventana Mineral y no apret Guardar,
despus se cambi a la ventana SAG y cambi otro valor y luego apret
Guardar, slo guardar el cambio en la variable de la ventana SAG.
Por ltimo, apretando el botn SEC (secundaria) se despliega la ventana con
variables y parmetros relacionados con el pozo, la bomba centrfuga, las
bateras de hidrociclones y los molinos de bolas (figura 20).

Pgina 26

Figura 20. Ventana de Configuracin del SEC
Una vez cerrada la ventana de configuracin el simulador asume que el usuario
actual tiene acceso a la configuracin, luego la prxima vez que uno elija del
men Configurar la opcin Mineral en la misma sesin, no se preguntar
por la password. Para quitar este privilegio al usuario, en el men Configurar
elija Logout.
Existen algunos valores que utiliza el simulador y que no estn disponibles en
ests ventanas, principalmente porque su cambio no debera ser necesario para
los usuarios corrientes. Para realizarlo se requiere de un profundo conocimiento
tanto de los modelos como de la programacin. Tal es el caso de las matrices de
molienda secundaria. Para mayor informacin al respecto ver el Manual
Tcnico de este simulador [10].


Pgina 27
8 Indice de Ilustraciones
Figura 1. Partiendo el simulador desde Matlab................................................ 2
Figura 2. Pantalla inicial al partir el simulador ................................................ 3
Figura 3. Layout del circuito de molienda ....................................................... 4
Figura 4. Pantalla de molienda SAG................................................................ 5
Figura 5. Pantalla de molienda secundaria (SEC)............................................ 8
Figura 6. Pantalla de Layout .......................................................................... 10
Figura 7. Realizando Zoom............................................................................ 12
Figura 8. Llamando a la ventana de parmetros de los PID........................... 14
Figura 9. Parmetros de los PID .................................................................... 15
Figura 10. Lazos de control.............................................................................. 16
Figura 11. Trabajando con el control PI de hold-up ........................................ 17
Figura 12. Trabajando con los controles PID de la molienda secundaria........ 18
Figura 13. Men Simulador Guardar ........................................................ 20
Figura 14. Archivo de sesin de simulacin.................................................... 21
Figura 15. Rescatando con EXCEL archivo de sesin de simulacin ............. 21
Figura 16. Archivo de sesin de simulacin abierto desde EXCEL................ 22
Figura 17. Ventana de password...................................................................... 24
Figura 18. Ventana de Configuracin de Mineral............................................ 24
Figura 19. Ventana de Configuracin del SAG............................................... 25
Figura 20. Ventana de Configuracin del SEC................................................ 26


Pgina 29
9 Referencias
[1] Amstica Valenzuela, Rodrigo. Memoria para optar al ttulo de Ingeniero
Civil Electricista: Nuevos Modelos de Molienda Semiautgena para
Simulacin y Control 1995.
[2] R. Amstica, G. Gonzlez, J. Menacho and J. Barra (1996). A
mechanistic state equation model for a semiatogenous grinding mill.
International Journal of Mineral Processing, v. 44-45, 349-360.
[3] Orellana Araya, Roberto. Memoria para optar al ttulo de Ingeniero Civil
Electricista: Simulador DSAG 1999.
[4] V. K. Karra, Development of a Model for Predicting the Screening
Performance of a Vibrating Screen CIM Bulletin, April 1979, pp. 167
171.
[5] JKSIMMET, Steady State Mineral Processing Simulator, User Manual,
JKMRC, 1991. pp.: A34-A49
[6] Mular A. L. & Bhauppu, R. B. 1978 Mineral Processing Plant Desing.
pp.: 335-336.
[7] A. Lynch and T. C. Rao, Modeling and scale-up of hydrocyclone
classifiers, Proc. 11
th
International Mineral Processing Congress, Cagliari,
Italy. 4-25.
[8] L. R. Plitt, A Mathematical Model of Hydrocyclone Classifier CIM
Bulletin December 1976 p. 114.
[9] Gonzlez, G. Informe Final del Proyecto CONICYT N124: Control
ptimo de una planta concentradora de minerales de cobre. Departamento
de Electricidad, U. de Chile, Sept. 1971
[10] Orellana, Roberto Manual Tcnico DSAG Proyecto DOCENMIN
Departamento de Ingeniera en Minas y Departamento de Ingeniera
Elctrica, Universidad de Chile, 1999.

Teoria Simulacion Molienda

Cargado por

Copyright:

Formatos disponibles

Teoria Simulacion Molienda

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

Teoria Simulacion Molienda

Cargado por

Copyright:

Formatos disponibles

MODELO, CONTROL Y SIMULADOR DE

PLANTA DE MOLIENDA SEMIAUTGENA Y MOLIENDA SECUNDARIA

= Hidrociclones - Fraccin del flujo de slidos que es arrastrado

= Fraccin del flujo de slidos que es arrastrado directamente por el

representa la suma de los flujos de mineral

es causada por la columna de pulpa en el pozo. La

(flujo volumtrico por la bomba)

También podría gustarte