Minería y Geología

E-ISSN: 1993-8012
revistamg@ismm.edu.cu
Instituto Superior Minero Metalúrgico de
Moa 'Dr Antonio Nuñez Jiménez'
Cuba

Coello-Velázquez, Alfredo Lázaro

Procedimiento para la determinación de la carga circulante en circuitos cerrados de
trituración y molienda
Minería y Geología, vol. 31, núm. 2, abril-junio, 2015, pp. 66-79
Instituto Superior Minero Metalúrgico de Moa 'Dr Antonio Nuñez Jiménez'
Holguín, Cuba

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

Cómo citar el artículo

Número completo
Sistema de Información Científica
Más información del artículo Red de Revistas Científicas de América Latina, el Caribe, España y Portugal
Página de la revista en redalyc.org Proyecto académico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
 Alfredo Lázaro Coello - Velázquez

Con el propósito de disminuir los errores (residuos) que tienen lugar

en la determinación de la carga circulante en los circuitos cerrados de
molienda o trituración, aplicando los métodos de los mínimos
cuadrados y de los multiplicadores de Lagrange a las ecuaciones de
balance, se desarroll ó un procedimiento que minimiza las desviaciones
e incongruencias en el balance de distribución de las clases
granulométricas de los productos del esquema. Tomando como punto
de partida los datos experimentales obtenidos por el autor en 1993 se
validó el procedimiento propuesto. Los resultados muestran que el
procedimiento aplicado permite el mejor ajuste de cálculo de la carga
circulante y a partir de sus valores se recalculan los contenidos de las
clases granulométricas de los productos del esquema que eliminan las
incongruencias y errores de la distribución del contenido de las clases
granulométricas generados en la toma y procesamiento general de las
muestras .

balance material; carga circulante; circuito cerrado;

función de Lagrange; método de los mínimos cuadrados; molienda;
trituración

Recibido: 5 febrero 2015 Aprobado: 19 marzo 2015

 With the objective to eliminate errors (residues) during the estimation
of the circulating load in crusher/grinding closed circuit, a procedure
was developed that minimize the deviatio ns and inconsistencies that
occur in the distribution balance of granulometric classes of the flow
sheet product by applying leas t square and Lagrange multipliers to
balance equations. The proposed procedure is approved taking as start
point the experimental data obtained by the author in 1993. Results
show that applied procedure allows the best fit circulating load
calculation and from her resulting data, the size contents of the flow
sheet products are then re - estimated which eliminates inconsistencies
and errors in the content distribution of the granulometric classes
generated in sampling operations

mass balance ; circu lating load; closed circuit; Lagrange

multiplier; lest square method ; grinding; c rushing

Received : 2 november 2013 Accepted : 17 september 2014

En la tecnología carbonato amoniacal, la molienda de minerales lateríticos
garantiza la superficie específica requerida del mineral para reducir el
níquel y el cobalto distribuidos en la red cristalina de los principales
portadores. Esto último constituye un elemento imprescindible para
conseguir los índices de extractable planificados en el proceso de
lixiviación y la producción final. En este sentido, la mena debe ser
reducida dimension almente hasta un tamaño de 80 - 85 % de la clase
menor a 0,074 μm (Coello - Velázquez 1993).

En este contexto, la molienda del mineral consume aproximadamente el

21 % de toda la energía eléctrica demandada por la tecnología
carbonato - amoniacal (Laborde, Coell o- Velázquez & Angulo 2003),
superada exclusivamente por el proceso de lixiviación.

El esquema tecnológico de la planta de molienda incluye una operación

de molienda en circuito cerrado con un separador de aire y dos etapas de
precipitación de polvos en ba terías de ciclones con un electro - filtro
(Coello - Velázquez 1993). En estos esquemas, el proceso de
fragmentación se ve afectado no solo por las fluctuaciones de las
propiedades físicas y físico - mecánicas del mineral en la alimentación,
propias de su distri bución espacial en el yacimiento, sino también por una
gama de factores tecnológicos que incluyen la carga circulante.

La carga circulante, por concepto, es la relación porcentual de la cantidad

de mineral que retorna al molino para su remolienda (arenas del
separador de aire) y la alimentación fresca al molino.
En los circuitos cerrados de trituración/molienda, la influencia de la carga
circulante está marcada por la incidencia en los niveles de acumulación
de material en el interior del molino sobre el rendimiento del molino y los
indicadores energo - tecnológicos.

La acumulación de mineral en el interior del molino, sin ser una suma

exacta de la productividad inicial y la carga circulante, define el desarrollo
del proceso no solo en el orden cualitativo, sino también cuantitativo
(Coello - Velázquez, Marrero - Ramírez & Hernández - Flores 1994). Así que
la calidad y cantidad del producto final es función de dicha variable. Para
garantizar los niveles máximos de la productividad del producto molido es
necesario m antener los niveles de acumulación en valores cercanos al
valor crítico.
Wills y Napier - Munn (2006) plantearon que, probablemente, son las
fluctuaciones de la densidad y la dureza en el mineral alimentado los que
con mayor fuerza influyen en el balance de clases de la molienda. Para
las lateritas de Moa, Coello - Velázquez, Menéndez - Aguado & Laborde
(2008) establecieron que la relación serpentina/limonita influye
decisivamente en el rendimiento y la distribución de las clases en el
balance del molino.
Comúnm ente, en la práctica industrial, la carga circulante se calcula a
partir de los resultados del muestreo de la descarga del molino y los
productos del clasificador, según la ecuación 1.
(1)

Donde:

C - es la carga circulante, unidades;

f´, c y p –porcentaje en peso de la clase i en la descarga del molino,

retorno del separador y su producto fino, respectivamente .
No siempre los cálculos por la ecuación 1 pueden ser considerados como
válidos, en la ma yoría de los casos están sujetos a los errores graves,
casuales y sistemáticos que tienen lugar durante el muestro. En el valor
de las variables f´, c y p se resumen los errores que se cometen en las
operaciones de tratamiento de las muestras, de su anális is
(granulométrico en este caso), los errores en la toma de muestras y las
fluctuaciones en la alimentación. La presencia de valores negativos de la
carga circulante delata la existencia de fluctuaciones en el balance de las
clases de los productos de la i nstalación.
Wills y Napier - Munn (2006), para resolver estos problemas en el balance
material de las operaciones de beneficio de minerales, aplicaron el
multiplicador de Lagrange para optimizar los residuos. El procedimiento
expuesto por los autores antes c itados, en su sentido expreso, es
aplicable a los circuitos abiertos de trituración o molienda con cribado o
clasificación preliminar, para la determinación de los flujos de los
productos, no así para los circuitos cerrados, donde el balance material
de cl ases es un tanto diferente.
Para el circuito abierto la ecuación principal del balance material de
cualquier instalación se rige por la ecuación 2.

(2)

La ecuación 2, que aparece en Wills y Napier - Munn (2006), es para el

cálculo del balance material en operaciones de separación mecánica en
condiciones de “estado de equilibrio del proceso”, considerando que el
error del cierre del balance es igual a cero. Al adecuar esta ecuación a un
esquema de circuito cerrado de tritura ción/molienda es obvio que los
parámetros cambian su concepto: se sustituye el contenido de los
componentes útiles (f, p y c) por el contenido de las clases
granulométricas.
Por otra parte, como F’ (alimentación a la operación de clasificación) es la
suma de la alimentación fresca en el molino (F) y la recirculación al
molino (carga circulante), por tanto, la alimentación fresca no es
precisamente la suma de los productos del esquema, sino que es igual a
P (flujo del producto final). Esto hace que no sea p osible emplear
directamente la ecuación 2 en el cálculo del esquema de la Figura 1.

El objetivo del presente trabajo es exponer un procedimiento para el

cálculo de la carga circulante optimizando los errores residuales.

La propuesta considera el cálculo de la carga circulante en los circuitos

cerrados de molienda/trituración a partir de los datos experimentales,
resultantes del muestreo en los puntos de descarga del molino (f´), y en
el reboso del clasificador (p) y en su s arenas (c), para luego ajustar la
data a los valores estimados de la carga circulante.

De modo general, el balance material de los procesos de reducción de

tamaño está formado por dos ecuaciones básicas: la ecuación de salida
de los productos y la ecuac ión de contenido de las clases. Para la
instalación de la Figura 1 el balance material toma la siguiente forma:

(3)

En consecuencia:

(4)

(5)

Donde:
F, F´, P y C - son los rendimientos (salidas) de los flujos de alimentación
fresca, descarga del molino, producto fino y retorno al molino (carga
circulante) relativos a la alimentación fresca, unidades.
 - es el residuo o error.
De la ecuación (5) podemos obtener la ecuación de la distribución de las
clases de tamaño de la instalación:

(6)

La ecuación es válida solo si el error o residuo es igual a cero ( , lo

cual es inamisible en términos experimentales, por lo tanto:
(7)

El objetivo del método es encontrar el valor de la carga circulante que

minimiza suma de los cuadrados de los residuos de cada parámetro en la
ecuación 6. Entonces se puede escribir que:
(8)
∑

Sustituyendo en la ecuación 8 se obtiene que:

(9)
∑

(10)
∑ ∑ ∑

De la ecuación (10) se deduce que para cualquier valor de C, la función S

es igual a cero solo para la condición en que los valores de sus
parámetros f´, p, c sean coherentes. No obstante, la función S tiene su
mínimo si su primera derivada es igual a cero, por lo tanto, el valor de C
para esta condición, minimiza los residuos.
(11)

(12)
∑ ∑

(13)
∑ ∑

Entonces el valor de Ċ que minimiza →min se puede calcular por la

fórmula 14:
∑ (14)
∑
 que satisfacen el mejor ajuste de Ċ.
La ecuación 6 puede ser rescrita como:
(15)

Escribiendo la ecuación 15 en función de la carga circulante reajustada Ċ

por la fórmula 14.

(16)

Entonces:

(17)

(18)

(19)

Donde:
- son los valores ajustados de los parámetros de la ecuación
(14), respectivamente;

- son los valores de los residuos de los parámetros de la

ecuación (14).
Siendo así, la ecuación 15 puede escribirse de la siguiente manera:

(20)

Si restamos la ecuación 15 y la 16 entonces se obtiene que:

(21)

Sustituyendo las fórmulas 17, 18 y 19 en la ecuación 21 se obtiene:

(22)

Para la ecuación 21, la suma de los cuadrados que hace mínimo cada
parámetro toma la siguiente forma:

(23)
∑

De acuerdo con Boizán (1984), la función de Lagrange o lagraniana

(ecuación 24) puede ser utilizada para la determinación de las posibles
localizaciones del óptimo (máximo o mínimo). Wills y Napier - Munn
(2006) aseguraron que los problemas de minimización de los errores
pueden ser resueltos con el método de Lagrange.
 ∑ (24)

Donde:

F – es la función de Lagrange;

λ – es el multiplicador de Lagrange;
R(x) - es la restricción.

Para este caso particular, la ecuación 25 es la sustracción de las dos

restricciones formadas por las ecuaciones 15 y 16:

(25)
Entonces para el caso en cuestión la función de Lagrange toma la
siguiente forma:

(26)
∑ ∑

Para determinar el valor de los residuos basta con derivar e

igualar a cero la ecuación 25 (función de Lagrange) por cada uno de
ellos:
(27)

(28)

(29)
Siguiendo el mismo proce dimiento obtenemos:

(30)

(31)

Como el valor de landa (multiplicador de Lagrange) se desconoce, si

derivamos la ecuación 25 por el mismo, y sustituimos los valores de los
residuos de cada parámetro obtenemos:

(32)

(33)

El valor de puede calcularse por la ecuación 15.

La instalación del CEDINIQ (Centro de Investigaciones del Níquel) cuenta
con un circuito cerrado de molienda formado por un molino cilíndrico -
cónico 1,85 X 0,56 m y un separador de aire con ventilación interna de
3,5 t/h (Figura 1), donde f es el contenido de la clase de cálculo en la
alimentación fresca; f´, p, c mantienen significado y la saeta señala la
dirección del flujo.

Figura 1. Esquema de la instalación de molienda del CEDINIQ (Coello -

Velázquez 1993).
Para el desarrollo del procedimiento los dato s iniciales fueron tomados de
Coello - Velázquez (1993). Los experimentos se llevaron a cabo en las
instalaciones del Centro de Investigaciones del Níquel, en período marzo -
junio de 1990, bajo la dirección del autor.

En la Tabla 1 apar ecen los valores de la carga de circulante calculados

por la fórmula 1 para diferentes niveles de la alimentación. A juzgar por
los valores de la tabla, existen valores de la carga circulante que poseen
signo negativo.
Estos valores de carga circulante pi erden todo sentido físico y
básicamente aparecen por fluctuaciones aleatorias de las propiedades
físicas y físico - mecánicas del material en la alimentación fresca o por
errores, casi siempre graves, en las operaciones de muestreo y las
determinaciones anal ísticas que provocan desviaciones en el balance de
las clases de tamaño.
El desbalance de la distribución de las clases de tamaño de manera
natural se reflejará en los valores ponderados de la carga circulante; en
este caso particular se debe a la distribu ción del contenido de las clases
granulométricas en los productos de descarga del molino, la recirculación
y el producto final. En este sentido, se ha desplazado hacia un
incremento de la clase intermedia 0,074 - 0,044 mm en el producto de
descarga del molin o con respecto al producto de recirculación, condición
que hace negativa el término de la derecha de la fórmula 1.

En la Tabla 2 se muestran los valores de la carga circulante determinados

por las fórmulas 1 y 14. En esta dirección, existen desviaciones
importantes entre los de valores Ċi y Ci, para todos los regímenes
tecnológicos. Vale recordar que los parámetros que conforman ambas
fórmulas son variables aleatorias, supeditadas no solo a los cambios de
valor (magnitud) en cada medición, dados por estar sujetos a la fuerte
influencia de errores casuales, sistemáticos e incluso graves, que tiene su
expresión de manera importante en su determinación y, por ende, de
forma muy particular, en el cálculo de los valores de la carga circulante
por la fórmula 1.
En los valores de Ci están representadas todas las fluctuaciones de los
parámetros de la fórmula 1, mientras que en los valores de Ċi, por
concepto, se minimiza el efecto de estas fluctuaciones e incongruencias
en los valores del contenido de los diferentes productos.
En el caso de la variable carga circulante, calculada por la fórmula 14
(Ċi), quedó demostrado matemáticamente en el epígrafe 2.1 que sus
valores hacen mínimo la suma cuadrática de los residuos, y por ende,
resultan más cercanos a los valores reales de esta variable aleatoria.
Dicho en otras palabras, los valores resultantes de “mejor ajuste” de la
carga circulante.
En la Tabla 3 aparecen los valores recalculados de los contenidos de las
clases granulométricas de los productos del esquema de la Figura 1.
Estos valores coinciden con el “mejor ajuste” de la carga circulante. El
cálculo del residuo ( ) con e stos “nuevos contenidos” por la fórmula 15
confirma la demostración matemática de la minimización de los residuos
en el cálculo de la carga circulante por la fórmula 14 por el procedimiento
propuesto.
 Tabla 1. Valores de la carga de circulante calculad os por la fórmula 1

>0,074 26,5 33,4 11,8 6,2 2,13

0,074 - 0,044 8,7 11,7 9,3 - 4,3 - 0,20
1 1 276,3 0,86
0,044 - 0,0 64,8 54,9 73,9 3,1 0,92
- 0,074 - 0 73,5 66,6 83,2 - 1,2 1,41
>0,074 51,6 56,3 22,1 11,3 6,28
0,074 - 0,044 10,6 13,3 13,6 - 13,5 - 1,11
2 2 234,2 0,86
0,044 - 0,0 37,8 30,4 64,3 2,2 3,58
0,074 - 0,0 48,4 43,7 77,9 - 11,3 6,28
>0,074 59,1 62,5 27,9 3,7 9,18
0,074 - 0,044 10,0 7,7 11,7 16,9 0,74
3 2 839,1 0,86
0,044 - 0,0 30,9 29,4 60,4 - 17,4 19,67
0,074 - 0,0 40,9 37,1 72,1 - 0,5 8,21
>0,074 33,3 58,3 22,1 - 0,6 0,45
0,074 - 0,044 20,8 9,3 14,3 11,9 - 0,57
4 1 276,3 3,71
0,044 - 0,0 45,9 32,1 63,6 - 11,2 1,28
0,074 - 0,0 66,7 41,4 77,9 0,7 0,44
>0,074 37,6 54,8 23,7 2,0 0,81
0,074 - 0,044 20,8 9,8 14,3 14,1 - 0,59
5 2 234,2 3,71
0,044 - 0,0 45,8 37,3 65,6 - 13,9 2,33
0,074 - 0,0 66,6 47,1 79,9 0,2 0,68
>0,074 43,9 54,7 28,6 - 2,7 1,42
0,074 - 0,044 10,4 8,0 5,5 8,9 - 2,04
6 2 839,1 3,71
0,044 - 0,0 45,7 37,3 65,9 - 6,2 2,40
0,074 - 0,0 56,1 45,3 71,4 2,7 1,42

Tabla 2. Valores de la carga circulante calculados por las fórmulas 1 (Ci) y 14

(Ċi)

Ċi
1 1 276,3 0,86 1,41 1,23 0,18
2 2 234,2 0,86 6,28 3,88 2,40
3 2 839,1 0,86 8,21 8,08 0,13
4 1 276,3 3,71 0,44 0,47 - 0,03
5 2 234,2 3,71 0,68 0,69 - 0,01
6 2 839,1 3,71 1,42 1,67 - 0,75
(*) Carga circulante por la clase de cálculo -0,074 mm.
 Tabla 3. Valores ajustados de los contenidos de las clases granulométricas de
los productos de la Figura 1

>0,074 24,6 34,4 12,6 0

0,074 -
10,0 11,0 8,7 0
0,044
1 1276,3 0,86
0,044 - 0,0 63,9 55,4 74,3 0

- 0,074 - 0 73,9 66,4 83,0 0

>0,074 50,2 57,4 22,4 0
0,074 -
12,2 12,0 13,3 0
2 2234,2 0,86 0,044
0,044 - 0,0 37,5 30,6 64,4 0
0,074 - 0,0 49,8 42,6 77,6 0
>0,074 58,9 62,7 27,9 0
0,074 -
9,0 8,6 11,8 0
3 2839,1 0,86 0,044
0,044 - 0,0 32,0 28,5 60,3 0
0,074 - 0,0 40,9 37,1 72,1 0
>0,074 33,5 58,2 21,9 0
0,074 -
15,6 11,0 17,8 0
4 1276,3 3,71 0,044
0,044 - 0,0 50,8 30,5 60,3 0
0,074 - 0,0 66,4 41,5 78,1 0
>0,074 36,8 55,1 24,2 0
0,074 -
15,3 12,0 17,6 0
5 2234,2 3,71 0,044
0,044 - 0,0 51,2 35,1 62,4 0
0,074 - 0,0 66,5 47,1 79,9 0
>0,074 44,6 54,3 28,3 0
0,074 -
8,2 9,4 6,3 0
6 2839,1 3,71 0,044
0,044 - 0,0 47,2 36,4 65,3 0
0,074 - 0,0 55,4 45,7 71,7 0

1. El procedimiento presentado permite el cálculo de la carga

circulante y el recálculo de los contenidos de las clases de los
diferentes productos en un circuito cerrado de molienda o
trituración , a partir de los datos experimentales generados por el
muestro y determinación del contenido de las clases de los
productos del esquema.
 2. La aplicación del procedimiento propuesto ha sido validado en la
práctica , sobre la base de datos experimentales , resultado del
muestreo y la determinación de los contenidos de las clases
granulométricas mediante el análisis de tamiz en las instalaciones
semi - industriales del C entro de Investigaciones de Níquel .
3. Los valores de la carga circulante determinados por el
procedimiento y los contenidos de las clases granulométricas
resul tantes de su ajuste al parámetro citado representan una
posibilidad concreta de minimizar los errores cometidos durante la
toma y el tratamiento de las muestras, y la determinación de los
contenidos de las clases granulométricas.

4. El procedimiento propuest o es válido para el cálculo de la carga

circulante y la determinación de los contenidos de las clases de los
productos de los circuitos cerrados de trituración o molienda a
cualquier escala.

Quiero agradecer al ingeniero Gilberto Font, director del CIL en el

momento del desarrollo de las pruebas, por su inmenso apoyo,
comprensión, confianza e inteligencia demostrada. A los ingenieros
Arnulfo Carcasés, Idalmis Peña y Mirelis Molina (estudiantes en aquel
momento) por su participación directa en la toma de las muestras.

BOIZÁN , M. 1984: Optimización (sistemas estáticos). Editorial Oriente ,

Santiago de Cuba, 174 p.
COELLO- VELÁZQUEZ , A. L. 1993: Sobersenstobovanie tecnologii izmelchenia
lateritovij rud. Dis. na izk. uch. step. Kan. tecn. nauk. Inst. Mejanobr.
Saint Petersburg.
COELLO- VELÁZQUEZ , A. L.; MARRERO- RAMÍREZ , S. & H ERNÁNDEZ - FLORES, A.
1994: Control de la carga interior del molino a través de terceros
parámetros. Minería y Geología 11 (1): 35 - 39.
COELLO- VELÁZQUEZ , A. L.; MENÉNDEZ - AGUADO , J. M. & LABORDE, R. 2008:
Grindability of lateritic nickel ore in Cuba. Powder Technology 182 (1) :
113 –115 .
LABORDE, R.; COELLO- VELÁZQUEZ , A. L. & ANGULO , O. 2003 : Racionalización
del consumo de ener gía en la molienda de los minerales lateríticos.
En: Seminario de Ingeniería Eléctrica SIE. UCLV. Las Villas, Cuba. 1 -
4.
W ILLS , B. & NAPIER- MUNN , T. 2006: Mineral Processing Technology: An
introduction to the practical aspect of ore treatment and mineral
recovery. Seventh Edition. Elsevier Science & Technology Books, 450
p.