UNIVERSIDAD NACIONAL DE CAJAMARCA

FACULTAD DE INGENIERÍA
ESCUELA ACADÉMICO PROFESIONAL DE INGENIERÍA DE MINAS

“Informe final – Programación lineal en las etapas de minería”

ASIGNATURA: INVESTIGACION DE OPERACIONES II

DOCENTE: MENDOZA SANTOS CARLOS EDUARDO

PRESENTADO POR:
✓ ALCANTAR CUZCO, DANNY
✓ CORNADO VARELA, ALEX JACINTO
✓ MARIN SÁNCHEZ, RONY
✓ PAISIG ARRASCUE, JAER DERMALI
✓ QUITO NEIRA, JAIR JESÚS
✓ SOLANO VARGAS, RICARDO ALFREDO

CAJAMARCA, PERÚ

2023

1
 INDICE

INDICE............................................................................................................................. 2
Indice de figuras ............................................................................................................... 5
Indice de tablas ................................................................................................................. 5
I. INTRODUCCION .................................................................................................... 7
II. OBJETIVOS ............................................................................................................. 7
2.1. Objetivo general ................................................................................................. 7
2.2. Objetivo especifico ............................................................................................ 8
III. DESARROLLO ..................................................................................................... 8
IV. OPTIMIZACIÓN DE LAS VIBRACIONES EN LA ETAPA DE VOLADURA 9
4.1. Los Parámetros De Voladura En Minería Superficial. ....................................... 9
4.2. Variables Que Pueden Influir En Las Características De La Vibración Del Suelo
Durante La Voladura ................................................................................................... 10
4.2.1. Variables controlables ............................................................................... 10
4.2.2. Variables no controlables:..........................................................................11
4.3. Parámetros Relacionados Con El Explosivo. ...................................................11
4.4. Parámetros Relacionados Con Las Propiedades Geomecánicas De La Roca. 12
4.5. Parámetros Básicos en la Geometría De Voladura........................................... 13
4.6. Parámetros Para Optimización De La Voladura .............................................. 14
4.7. Vibraciones Inducidas Por La Voladura De Rocas. ......................................... 14
4.7.1. Ondas sísmicas inducidas por la voladura. ............................................... 14
4.7.2. Ondas internas. ......................................................................................... 15
4.7.3. Ondas superficiales. .................................................................................. 15
4.8. Modelos Matemáticos Usados ......................................................................... 15
4.9. Caso De Estudio............................................................................................... 17
4.9.1. Monitoreo de las vibraciones para la ciudad cercana al tajo abierto de la
mina. 17
V. APLICACIÓN DE LA PROGRAMACIÓN LINEAL EN LAS OPERACIONES DE
CARGUIO Y ACARREO. ............................................................................................. 29
5.1. Introducción ..................................................................................................... 29
5.2. Realidad Problemática ..................................................................................... 29
5.3. Formulación Del Problema .............................................................................. 30
5.4. Justificación ..................................................................................................... 30
5.4.1. Justificación teórica .................................................................................. 30
5.4.2. Justificación práctica ................................................................................ 30
5.4.3. Justificación valorativa ............................................................................. 30
5.4.4. Justificación académica ............................................................................ 30

2
 5.5. Limitaciones..................................................................................................... 30
5.6. Objetivos .......................................................................................................... 30
5.6.1. Objetivo General....................................................................................... 30
5.6.2. Objetivos Específicos ............................................................................... 31
5.7. Bases Teóricas .................................................................................................. 31
5.7.1. Programación Lineal................................................................................. 31
5.7.2. Formulación de un modelo de Programación Lineal................................ 31
5.8. Descripción De Solver ..................................................................................... 32
5.8.1. Cargar el Complemento SOLVER ............................................................ 32
5.9. Número De Ciclos ........................................................................................... 33
5.10. Material Y Métodos ...................................................................................... 34
5.10.1. Unidad de estudio. ................................................................................ 34
5.10.2. Población. ............................................................................................. 34
5.10.3. Muestra. ................................................................................................ 34
5.10.4. Técnicas de recolección y análisis de datos .......................................... 35
5.10.5. Instrumentos.......................................................................................... 35
5.11. Procedimientos ............................................................................................. 35
5.12. Ubicación ..................................................................................................... 35
5.13. Geología Local Y Mineralización ................................................................ 36
5.14. Clasificación Del Mineral ............................................................................ 38
5.15. Optimización De Las Operaciones De Carguío Y Acarreo .......................... 38
5.15.1. Equipos y Productividad ....................................................................... 38
5.15.2. Descripción del Proceso de Carguío y Acarreo .................................... 38
5.16. Base De Datos Tomados De Campo ............................................................ 40
VI. APLICACIÓN DE CADENAS DE MARKOV Y SISTEMAS DE COLAS EN LA
OPERACIÓN UNITARIA DE ACARREO. .................................................................. 44
6.1. Introducción ..................................................................................................... 44
6.2. Objetivos .......................................................................................................... 44
6.2.1. Objetivo General....................................................................................... 44
6.2.2. Objetivos Específicos ............................................................................... 44
6.3. Caso de Aplicación .......................................................................................... 44
6.3.1. Análisis De Resultados ............................................................................. 50
6.3.2. Sistema De Colas ...................................................................................... 50
6.3.3. Análisis De Resultados En La Línea De Espera....................................... 52
6.3.4. Conclusiones............................................................................................. 52
VII. OPTIMIZACIÓN DEL PROCESO DE CHANCADO MEDIANTE
PROGRAMACIÓN LINEAL EN UNA PLANTA CONCENTRADORA DE COBRE 53

3
4
 INDICE DE FIGURAS

Figura 1 Diagrama conceptual que sirvió para postular el modelo matemático de

Holmberg & Persson, para distancias >= 50 metros ................................................. 17
Figura 2 Muestra la relación Velocidad pico de partícula (mm/s) Vs. La Distancia
Escalada (m/Kg^1/2). ................................................................................................... 18
Figura 3 Muestra la relación de la velocidad pico de partícula (mm/s) Vs. La
distancia (m). ................................................................................................................. 19
Figura 4 Muestra la relación de carga de mezcla explosiva comercial (kg) Vs. La
distancia (m). ................................................................................................................. 19
Figura 5: Muestra la relación vibración máxima (mm/s) Vs. Frecuencia (Hz). ..... 20
Figura 6 Muestra la relación de K Vs. Frecuencia (Hz)............................................ 21
Figura 7 Muestra la relación de Amplitud Vs. El tiempo de retardo entre taladros
(ms)................................................................................................................................. 22
Figura 8 Muestra la ubicación del sismógrafo con respecto a una voladura primaria.
........................................................................................................................................ 23
Figura 9 El software muestra el PPV (mm/s) de la onda elemental en dirección de la
ciudad Vs. El tiempo (ms). ........................................................................................... 23
Figura 10 Muestra el diseño de carga de un taladro. ................................................ 24
Figura 11 El software muestra el resultado obtenido de PPV (mm/s). .................... 24
Figura 12 Muestra la ubicación del sismógrafo con respecto a una voladura
primaria. ........................................................................................................................ 25
Figura 13 El software muestra el PPV (mm/s) de la onda elemental en dirección de
la ciudad Vs. El tiempo (ms). ....................................................................................... 25
Figura 14 Muestra el diseño de carga de un taladro. ................................................ 26
Figura 15 El software muestra el resultado obtenido de PPV (mm/s). .................... 26
Figura 16 Muestrala relación vibración máxima (mm/s) Vs. Frecuencia (Hz).
27
Figura 17 Muestra la relación de K Vs. Frecuencia (Hz). ........................................ 28
Figura 18 Configuración de Solver ............................................................................. 33
Figura 19 Metodología de la Investigación de Operaciones ..................................... 35
Figura 20 Mapa de Ubicación de Lagunas Norte ...................................................... 36
Figura 21 Geología Local Lagunas Norte .................................................................. 37
Figura 22 Clasificación del mineral Barrick Gold Corporation – Lagunas Norte. 38
Figura 23 Sistema Tradicional de Carguío - Acarreo................................................ 39

INDICE DE TABLAS

Tabla 1 Muestra el resumen de los resultados obtenidos. ......................................... 18

Tabla 2 Muestra el resumen de resultados de monitoreo de vibraciones en la ciudad.
........................................................................................................................................ 20
Tabla 3 Muestra los niveles de percepción ................................................................. 21
Tabla 4 Muestra el resumen de nuevos resultados de monitoreo de vibraciones en la
ciudad............................................................................................................................. 27
Tabla 5 Registro de tiempos NIVEL 4160. PALA hacia la chancadora ................... 40
Tabla 6 Registro de tiempos NIVEL 4070.CARGADOR hacia la chancadora....... 41
Tabla 7 Registro de tiempos NIVEL 4040. (CARGADOR) hacia chancadora ....... 42
Tabla 8 Cuadro de producción de plan semanal ........................................................ 42

5
Tabla 9 Cuadro de tiempos Sistema Tradicional ....................................................... 42
Tabla 10 Costos de Acarreo.......................................................................................... 43

6
 I. INTRODUCCION

La investigación de operaciones es una disciplina que utiliza métodos cuantitativos para

resolver problemas complejos relacionados con la toma de decisiones en empresas y
organizaciones. En la industria minera, la investigación de operaciones desempeña un
papel importante en la mejora de la eficiencia y la rentabilidad de las operaciones mineras.

La minería es una actividad que implica la extracción y procesamiento de minerales para

su posterior uso en una variedad de industrias. La minería es un proceso complejo que
involucra la coordinación de múltiples factores, como la geología, la ingeniería, la
planificación, la logística y la gestión de riesgos.

La investigación de operaciones puede ayudar a la industria minera a mejorar la toma de

decisiones en áreas clave como la planificación minera, la programación de la producción,
el control de calidad, la gestión de riesgos, la logística y la optimización del proceso.

La aplicación de técnicas de investigación de operaciones puede proporcionar soluciones

más precisas y eficientes a los problemas que enfrenta la industria minera, lo que puede
mejorar la rentabilidad, reducir los costos y mejorar la seguridad en el lugar de trabajo.

Por lo tanto, la investigación de operaciones es una herramienta valiosa para la industria

minera y puede ayudar a las empresas mineras a enfrentar los desafíos que surgen en un
entorno cada vez más competitivo y exigente.

II. OBJETIVOS

2.1. Objetivo general

• Identificar los principales desafíos y oportunidades en la industria minera que
pueden ser abordados mediante la aplicación de técnicas de investigación de
operaciones.

7
2.2. Objetivo especifico
• Describir las técnicas de investigación de operaciones más relevantes y útiles para
la industria minera, y cómo pueden aplicarse en diferentes áreas del proceso
minero.

• Analizar casos de estudio o ejemplos de aplicaciones exitosas de investigación de

operaciones en la industria minera, y evaluar los beneficios y limitaciones de estas
aplicaciones.

• Evaluar el impacto de la investigación de operaciones en la rentabilidad, la

eficiencia y la seguridad en la industria minera, y proporcionar recomendaciones
para su implementación efectiva.

III. DESARROLLO

En el presente trabajo buscamos hacer un análisis de programación en las diferentes

etapas de la minería, para esto hemos dividido el análisis en los siguientes temas:
• Cadena de markov y simulación de Montecarlo en la operación unitaria de
perforación.
• Optimización de las vibraciones en la etapa de voladura.
• Aplicación de programación lineal en las operaciones de carguío y acarreo.
• Aplicación de cadenas de markov y sistemas de colas en la operación unitaria de
acarreo.
• Optimización del proceso de chancado mediante programación lineal.

8
 IV. OPTIMIZACIÓN DE LAS VIBRACIONES EN LA

ETAPA DE VOLADURA

La voladura en minería superficial es una técnica utilizada para extraer materiales de la

superficie de la tierra, como minerales, rocas y tierra. Es una de las técnicas más comunes
y efectivas utilizadas en la minería a cielo abierto. La voladura implica la detonación
controlada de explosivos colocados en agujeros perforados en la roca o el suelo. Los
explosivos liberan energía que fragmenta la roca y permite la extracción de los materiales.
Para realizar una voladura en minería superficial, se debe seguir un proceso cuidadoso y
controlado que incluye la planificación de la voladura, la perforación de los agujeros de
carga, la carga de explosivos, el retardo de las detonaciones y la limpieza del área después
de la voladura. Es importante tener en cuenta que la voladura en minería superficial puede
tener efectos negativos en el medio ambiente y en la salud y seguridad de los trabajadores
y comunidades cercanas. Por lo tanto, se deben seguir prácticas adecuadas y cumplir con
las regulaciones y normativas locales y nacionales para minimizar los impactos negativos
y maximizar la seguridad y la sostenibilidad en la operación minera.

4.1. Los Parámetros De Voladura En Minería Superficial.

Son aquellos que se utilizan para planificar y ejecutar de manera segura y eficiente
la voladura en una operación minera a cielo abierto. Estos parámetros pueden variar
dependiendo de la naturaleza del material a extraer, las condiciones geológicas y las
características del terreno, pero algunos de los parámetros comunes incluyen:
• Diámetro de perforación: Es el diámetro de los agujeros que se perforan en la
roca para colocar los explosivos. Este diámetro varía según la carga explosiva a
utilizar, el tipo de roca, la profundidad de la perforación y otros factores.

• Longitud de perforación: Es la longitud total del agujero que se perfora en la

roca. Esta longitud también varía según los factores mencionados anteriormente,
y es importante para asegurar que la carga explosiva se distribuya adecuadamente
en la roca.

• Espaciamiento de perforación: Es la distancia entre cada uno de los agujeros de

perforación. Este parámetro también varía según los factores mencionados
anteriormente, y es importante para garantizar una fragmentación adecuada de la
roca y evitar daños en la pared del pozo o la superficie de trabajo.

• Carga explosiva: Es la cantidad de explosivos que se coloca en cada agujero de

perforación. Esta carga varía según la naturaleza del material a extraer y las
características del terreno. Es importante asegurarse de que la carga explosiva no
sea excesiva ni insuficiente para garantizar una voladura eficiente y segura.

9
 • Tiempo de retardo: Es el tiempo que transcurre entre la detonación de cada carga
explosiva. Este tiempo se ajusta de forma cuidadosa para que la secuencia de
explosiones proporcione una fragmentación adecuada de la roca y minimice las
vibraciones y el impacto en la superficie.

• Tipo de explosivo: Existen diferentes tipos de explosivos, cada uno con sus
propias características y aplicaciones. La elección del tipo de explosivo depende
del material a extraer, las características del terreno y los requisitos de seguridad.

• Control de vibraciones: Es importante controlar las vibraciones generadas por la

voladura para evitar daños a las estructuras cercanas y reducir la emisión de ruido.
Los parámetros de voladura se ajustan cuidadosamente para minimizar las
vibraciones y el ruido generados por la voladura.
En resumen, los parámetros de voladura en minería superficial son importantes para
garantizar una voladura segura y eficiente en una operación minera a cielo abierto.
La elección y ajuste adecuado de estos parámetros depende de varios factores,
incluyendo el material a extraer, las características del terreno y los requisitos de
seguridad y ambientales.

4.2. Variables Que Pueden Influir En Las Características De La Vibración

Del Suelo Durante La Voladura
4.2.1. Variables controlables
• Cantidad y tipo explosivo: La elección de la cantidad y tipo de explosivo
utilizados en la voladura puede ser controlada y ajustada para minimizar los
efectos de las vibraciones del suelo. La reducción de la cantidad de explosivo o el
uso de explosivos de baja vibración puede ayudar a disminuir la intensidad de las
vibraciones.

• Diseño de la voladura: El diseño de la voladura, incluyendo la configuración de

los taladros, la profundidad y el diámetro de los taladros, la carga y la secuencia
de detonación, puede ser controlado para minimizar los efectos de las vibraciones
del suelo. El uso de técnicas de voladura controlada, como la reducción de la carga
por taladro, puede ayudar a reducir la intensidad de las vibraciones.

• Control de la secuencia de detonación: El control de la secuencia de detonación

de los explosivos puede ser una variable importante para controlar las
características de la vibración del suelo. La secuencia de detonación puede ser
controlada para evitar que varias cargas detonen simultáneamente y produzcan
vibraciones excesivas.

• Distancia de la fuente de la voladura: La distancia entre la fuente de la voladura

y el punto de medición es una variable que puede ser controlada. Cuanto mayor

10
 sea la distancia entre la fuente de la voladura y el punto de medición, menor será
la intensidad de las vibraciones.

• Monitoreo de las vibraciones: La monitorización de las vibraciones del suelo en

tiempo real puede ser utilizada para controlar las características de la vibración
del suelo y ajustar la voladura en consecuencia.
Es importante tener en cuenta estas variables controlables al diseñar y planificar la
voladura para minimizar los efectos de las vibraciones en el suelo y reducir los
impactos ambientales.

4.2.2. Variables no controlables:

• Geología del suelo: Las características geológicas del suelo, como la densidad, la
rigidez, la profundidad de la capa superficial y la presencia de fallas y grietas, son
variables no controlables que influyen en la propagación y la amplitud de las
ondas de vibración generadas por la voladura.

• Topografía y configuración del terreno: La topografía y la configuración del

terreno en el sitio de la voladura son variables no controlables que pueden influir
en la propagación de las ondas de vibración y en la magnitud de las vibraciones
registradas.

• Condiciones ambientales: Las condiciones ambientales, como la humedad, la

temperatura y la presión atmosférica, son variables no controlables que pueden
afectar la propagación y la amplitud de las ondas de vibración.

• Interferencias externas: Las interferencias externas, como la presencia de otras

fuentes de vibración, pueden afectar las mediciones de las vibraciones del suelo y
complicar la evaluación de los impactos de la voladura.

4.3. Parámetros Relacionados Con El Explosivo.

• Tipo de explosivo: El tipo de explosivo utilizado es un factor importante a
considerar en la voladura de rocas. Los tipos de explosivos comunes incluyen
explosivos de alta velocidad, emulsiones, geles, ANFO (nitrato de amonio y fuel
oil) y emulsiones sensibles al agua. Cada tipo de explosivo tiene características
diferentes en términos de velocidad de detonación, densidad, energía, resistencia
a la humedad, estabilidad y costo.

11
 • Carga de explosivo: La cantidad de explosivo utilizada en cada taladro es un
parámetro crítico. Una carga insuficiente puede resultar en una voladura
incompleta, mientras que una carga excesiva puede generar vibraciones excesivas
y dañar el medio ambiente. La carga de explosivo se determina en función de la
densidad del explosivo y la longitud y diámetro del taladro.

• Velocidad de detonación: La velocidad de detonación del explosivo afecta la

eficiencia de la voladura y la producción de energía. Un explosivo de alta
velocidad proporciona una onda de choque más rápida y puede generar una mayor
cantidad de energía en un período más corto. Por otro lado, un explosivo de baja
velocidad puede proporcionar una voladura más controlada con menos
vibraciones.

• Sensibilidad a la iniciación: La sensibilidad a la iniciación del explosivo es un

parámetro importante que influye en la seguridad de la voladura. Algunos
explosivos pueden ser muy sensibles a la iniciación y pueden detonar
accidentalmente, lo que puede provocar daños a las personas y al equipo.

• Estabilidad del explosivo: La estabilidad del explosivo es un parámetro crítico a

considerar en el almacenamiento y transporte de los explosivos. Los explosivos
pueden ser sensibles a la humedad, la temperatura y otros factores ambientales, lo
que puede afectar su estabilidad y su capacidad para detonar.

4.4. Parámetros Relacionados Con Las Propiedades Geomecánicas De La

Roca.
• Resistencia a la compresión uniaxial (UCS): La resistencia a la compresión
uniaxial de la roca es un parámetro importante que influye en la capacidad de la
roca para soportar cargas. La UCS se mide mediante pruebas de laboratorio y se
utiliza para seleccionar la cantidad de explosivos y la carga adecuada en la
voladura.

• Resistencia a la tracción: La resistencia a la tracción de la roca es otro parámetro

importante que influye en la voladura. La resistencia a la tracción se mide
mediante pruebas de laboratorio y se utiliza para determinar la cantidad de
explosivos y la carga necesarias para fracturar la roca.

• Densidad de la roca: La densidad de la roca es un parámetro importante que

afecta la cantidad de energía necesaria para la voladura. Una roca más densa
requiere más energía para la voladura.

12
 • Tipo de roca: El tipo de roca es un factor importante a considerar en la voladura.
Las diferentes rocas tienen diferentes propiedades geomecánicas y pueden
requerir diferentes cantidades de explosivos y carga para la voladura.

• Presencia de fracturas: La presencia de fracturas en la roca puede afectar la

voladura. Las fracturas pueden actuar como puntos débiles en la roca y pueden
provocar que la voladura se desvíe de su trayectoria deseada.

• Características de la discontinuidad: Las discontinuidades son otro factor

importante a considerar en la voladura. Las discontinuidades pueden afectar la
voladura al actuar como puntos débiles en la roca o al provocar la reflexión y la
dispersión de la onda de choque.

• Orientación de las discontinuidades: La orientación de las discontinuidades es

importante para determinar la trayectoria de la voladura. Las discontinuidades
pueden afectar la voladura al actuar como puntos débiles en la roca o al provocar
la reflexión y la dispersión de la onda de choque.

4.5. Parámetros Básicos en la Geometría De Voladura

Son aquellos que determinan la forma y tamaño de la cavidad creada en la roca
después de la detonación. Algunos de los parámetros básicos más importantes son:
• Longitud de carga (L): La longitud de carga se refiere a la distancia entre el
extremo superior y el extremo inferior de la carga explosiva en la perforación.
Una longitud de carga más larga puede crear una cavidad más larga en la roca.
• Diámetro de la carga (D): El diámetro de la carga se refiere al diámetro de la
carga explosiva en la perforación. Un diámetro de carga más grande puede crear
una cavidad más grande en la roca.
• Espaciamiento de perforación (S): El espaciamiento de perforación se refiere a
la distancia entre dos perforaciones adyacentes en un patrón de perforación. Un
espaciamiento de perforación más estrecho puede resultar en una cavidad más
pequeña y más precisa, mientras que un espaciamiento de perforación más amplio
puede resultar en una cavidad más grande y menos precisa.
• Diámetro de perforación (d): El diámetro de perforación se refiere al diámetro
del agujero perforado en la roca. Un diámetro de perforación más grande puede
permitir una carga explosiva más grande y, por lo tanto, una cavidad más grande.
• Ángulo de perforación (α): El ángulo de perforación se refiere al ángulo formado
por la línea central del agujero y la superficie de la roca. Un ángulo de perforación
más estrecho puede resultar en una cavidad más precisa y controlada, mientras
que un ángulo de perforación más amplio puede resultar en una cavidad más
grande y menos precisa.

13
 • Patrón de carga: El patrón de carga se refiere a la disposición de los explosivos
en la perforación. Pueden utilizarse diferentes patrones de carga, como patrones
en línea recta, en abanico o en espiral, para crear diferentes geometrías de
voladura.

4.6. Parámetros Para Optimización De La Voladura

De todos los parámetros mencionados anteriormente no se podría tomar todos al
mismo tiempo para optimizar la voladura, pero si al menos más de uno de estos
parámetros nos ayudará en una voladura más optima en esta atapa importante de la
minería, en este caso de estudio los parámetros que se manipularon para la
optimización son principalmente el Tiempo de retardo y Control de vibraciones.

4.7. Vibraciones Inducidas Por La Voladura De Rocas.

Se conoce por vibraciones un fenómeno de transmisión de energía mediante la
propagación de un movimiento ondulatorio a través de un medio. El fenómeno
vibracional queda caracterizado por una fuente o emisor o sea un generador de
vibraciones; y por un receptor de las mismas. El fenómeno vibracional es
manifestado mediante un movimiento ondulatorio. Las vibraciones causadas por la
voladura de rocas, tienen el potencial de afectar estructuras y provocar molestias en
las personas. Éstas presentan parámetros que la definen, como amplitud de
aceleración, frecuencia y desplazamiento, pero es la amplitud de velocidad pico de
partícula a una frecuencia predominante la que mejor se relaciona con un potencial
de daño estructural. En general la energía vibracional es directamente proporcional a
la cantidad del explosivo detonado y matemáticamente se puede expresar de la
siguiente manera:

4.7.1. Ondas sísmicas inducidas por la voladura.

Cuando una mezcla explosiva comercial cargada dentro de un taladro es detonada,
se produce una onda de choque, la cual crea una serie de tensiones las que
producen la fractura miento del macizo | Reiterando que la voladura de rocas
genera vibraciones en el macizo rocoso, las cuales pueden tener diferentes efectos

14
 negativos sobre las infraestructuras; y, además, ocasiona daños a la misma roca
remanente, como fallas y diaclasas, las cuales podrían dar origen a deslizamientos
en los bancos.

4.7.2. Ondas internas.

Este tipo de ondas se propagan por el interior del macizo rocoso dentro de estas
se tienen:
• Ondas longitudinales, de compresión o principales (P).
• Ondas transversales, de cizallamiento o secundaria (S).
• Ondas verticales (V)

4.7.3. Ondas superficiales.

Entre las principales, se tienen las siguientes:
• Ondas rayleigh (R)
• Ondas love (L)

4.8. Modelos Matemáticos Usados

• Definición: Varios son los investigadores que han postulado modelos
matemáticos, que sirven para cuantificar el daño producido por las vibraciones
inducidas por la voladura de rocas.
Breve descripción de dichos modelos matemáticos:
• Modelo general.

El factor “D” distancia escalada toma en cuenta de la influencia de la distancia en

(m) y la cantidad de mezcla explosiva comercial (Kg.) por etardo.
• Devine y Duvall

15
• Langerfors

• Hendron y Dowding

• Holmberg y Persson.

Para campo lejano D >= 50 metros se deben aplicar los otros modelos matemáticos
mencionados anteriormente.

16
 Figura 1 Diagrama conceptual que sirvió para postular el modelo matemático de Holmberg & Persson, para distancias >= 50
metros

Fuente: Persson & Holmberg

4.9. Caso De Estudio

4.9.1. Monitoreo de las vibraciones para la ciudad cercana al tajo abierto de la
mina.

Para este caso se consideró siete voladuras primarias. El monitoreo de las

vibraciones se realizó de la siguiente manera:
• Tres monitoreos fueron realizados en la ciudad cercana a la mina.
• Cuatro monitoreos fueron realizados en el tajo en dirección de la ciudad.
La voladura primaria tuvo en promedio 150 taladros (distribuidas en 10 filas de
15 taladros cada una). Para llevar a cabo esta operación minera unitaria de
voladura de rocas se usó la mezcla explosiva comercial AN/FO, cuya cantidad por
taladro varía desde 330 Kg. hasta 650 Kg. (esta variación está en función directa
de la dureza del macizo rocoso). Usándose los siguientes tiempos de retardo:
• Los tiempos de retardo entre taladros utilizados fueron de 17 ms y 25 ms.
• Los tiempos de retardo entre filas utilizados fueron de 107 ms y 130 ms.
• Las distancias de las voladuras primarias hacia los puntos de monitoreo
varían desde 80.0 metros hasta 1550.0 metros.
• El equipo para el monitoreo de vibraciones se utilizó un sismógrafo
(Minimate Plus) con un geófono triaxial de superficie. dicho geófono es el
encargado de captar las vibraciones, a su vez las registra en la memoria,
para finalmente transmitir toda la data a una computadora. Los datos
resultantes son proporcionados en términos de velocidad y frecuencia. La
velocidad resultante se calcula en cada intervalo de tiempo de acuerdo a
las normas técnicas, de tal manera que se puede obtener una velocidad
resultante máxima.

17
 Tabla 1 Muestra el resumen de los resultados obtenidos.

Con los datos obtenidos y mostrados en la tabla, se empieza a crear el modelo

predictivo de campo lejano para la ciudad, para el caso se seleccionó y uso el
modelo matemático de Devine; el cual relaciona las variables de PPV suma
(mm/segundo), distancia al punto de monitoreo (metros) y la carga de la mezcla
explosiva por retardo (Kilogramos), cuyos resultados se muestra en el gráfico
siguiente.

Figura 2 Muestra la relación Velocidad pico de partícula (mm/s) Vs. La Distancia Escalada (m/Kg^1/2).

Con el modelo matemático de Devine, ecuación de modelo predictivo de campo

lejano obtenida al 85% de confiabilidad, se puede construir la gráfica de PPV vs
Distancia, con la cual se puede calcular la distancia segura a la que se puede
realizar una voladura primaria (según los Kilogramos de carga por retardo a
utilizar), de tal manera de no sobrepasar el valor límite de 50.0 mm/segundo.
Dicho límite ha sido tomado del criterio de daño USBM, que es la norma más
conservadora postulada a la fecha.

18
 Figura 3 Muestra la relación de la velocidad pico de partícula (mm/s) Vs. La distancia (m).

Por otro lado, con la ecuación del modelo predictivo de campo lejano obtenida
(85% de confiabilidad), se puede construir el gráfico siguiente, donde se muestra
la carga por retardo vs Distancia. Enfatizando se debe tener presente siempre como
referencia el valor límite de 50.0 mm/segundo; es decir, con este valor como
referencia se puede calcular la distancia segura a la que se puede realizar una
voladura primaria, según los Kilogramos de carga por retardo que se van a utilizar.

Figura 4 Muestra la relación de carga de mezcla explosiva comercial (kg) Vs. La distancia (m).

Por otro lado, los resultados de monitoreo de vibraciones en la ciudad han sido
relacionados con la norma internacional de percepción humana “USBM–
Goldman”.
A continuación, se presenta el cuadro resumen de resultados de monitoreo y
gráfico de los resultados relacionados a la norma internacional.

19
 Tabla 2 Muestra el resumen de resultados de monitoreo de vibraciones en la ciudad.

Figura 5: Muestra la relación vibración máxima (mm/s) Vs. Frecuencia (Hz).

De la figura se observa que el 55.5 % de los datos se encuentran en la zona

“Perceptible” por el ser humano y el 44.5 % de los datos se encuentran en la zona
“Imperceptible” por el ser humano.
A continuación, se presenta otro gráfico N° 10.6, donde muestra los resultados de
monitoreo relacionados a la norma internacional de percepción humana según
“Steffens”:
Cuya relación matemática, es la siguiente:

Donde:
V = Vibración máxima (mm/s
F = Frecuencia (Hz)

20
 Tabla 3 Muestra los niveles de percepción

VALOR DE K NIVEL DE PERCEPCIÓN

< 0,1 No se siente

0,1 Comienza a percibirse

0,25 Escasamente perceptible

0,63 Perceptible

1,6 Fácilmente perceptible

4,0 Fuertemente detectable

10,0 Muy fuertemente detectable

Figura 6 Muestra la relación de K Vs. Frecuencia (Hz).

21
 De la figura se observa que el 66.7 % de los datos se encuentran en la zona
“Escasamente Perceptible” por el ser humano, y el 33.3 % de los datos se
encuentran en la zona “Comienza a percibirse” por el ser humano.
Para mejorar y reducir la percepción humana por las vibraciones inducidas por la
voladura de rocas se utilizó el método de onda elemental.
Se capturó una onda elemental en dirección de la ciudad y se calculó el tiempo de
retardo entre taladros adecuado para reducir las vibraciones producidas por efectos
de la voladura, dicho tiempo es de 27 ms., para su mejor compresión, se muestra
en el siguiente gráfico.

Figura 7 Muestra la relación de Amplitud Vs. El tiempo de retardo entre taladros (ms).

Una vez obtenido el tiempo de retardo entre taladros, seguidamente se realiza un

comparativo de resultados de PPV suma.

• Voladura Primaria N° 01:

Para esta voladura, se cuenta con la siguiente información:
• Cantidad de taladros = 150 (10 filas, de 15 taladros).
• Tiempo entre taladros = 17 ms
• Tiempo entre filas = 107 ms.
• Distancia de la voladura hacia el sismógrafo “A” = 60 m. (sismógrafo
ubicado en dirección de la ciudad).

22
 Figura 8 Muestra la ubicación del sismógrafo con respecto a una voladura primaria.

En el software de simulación de vibraciones denominado “Shotplus” se introduce

la onda elemental (vector suma), cantidad de carga por taladro (Kg), constantes
halladas en el modelo de Devine (K y α) y la velocidad de la onda “P”, con estos
datos se obtiene un resultado de PPV suma = 155.3 mm/segundo.

Figura 9 El software muestra el PPV (mm/s) de la onda elemental en dirección de la ciudad Vs. El tiempo (ms).

23
 Figura 10 Muestra el diseño de carga de un taladro.

Figura 11 El software muestra el resultado obtenido de PPV (mm/s).

• Voladura primaria N° 02:

Para esta voladura, se cuenta con la siguiente información:
• Cantidad de taladros = 150 (10 filas de 15 taladros).
• Tiempo entre taladros = 27 ms. (Cálculo de onda elemental)
• Tiempo entre filas = 170 ms.
• Distancia de la voladura hacia el sismógrafo “A” = 60 m. (sismógrafo
ubicado en dirección de la ciudad).

24
 Figura 12 Muestra la ubicación del sismógrafo con respecto a una voladura primaria.

En el software de simulación de vibraciones denominado “Shotplus”, se introduce

la onda elemental (vector suma), cantidad de carga por taladro (Kg), constantes
obtenidas en el modelo matemático de Devine (K y α) y la velocidad de la onda
“P”, con estos datos se obtiene un resultado de PPV suma = 134.9 mm/segundo.

Figura 13 El software muestra el PPV (mm/s) de la onda elemental en dirección de la ciudad Vs. El tiempo (ms).

25
 Figura 14 Muestra el diseño de carga de un taladro.

Figura 15 El software muestra el resultado obtenido de PPV (mm/s).

Por lo tanto, al utilizar el tiempo de retardo entre taladros de 27 ms (cálculo de

onda elemental) se reduce las vibraciones en un 13 %.
Utilizando también el tiempo de retardo entre taladros de 27 ms, se realizó las
simulaciones para obtener nuevos resultados de PPV y frecuencia en puntos de
monitoreos ubicados en la ciudad y se obtuvo la siguiente información que se
muestra en la tabla.

26
 Tabla 4 Muestra el resumen de nuevos resultados de monitoreo de vibraciones en la ciudad.

Los nuevos resultados de monitoreo de vibraciones en la ciudad, han sido

relacionados a la norma internacional de percepción humana “USBM –
Goldman”, a continuación, se presenta el gráfico de los resultados relacionados a
la norma internacional.

Figura 16 Muestra la relación vibración máxima (mm/s) Vs. Frecuencia (Hz).

De la figura se observa que el 11.1 % de los datos se encuentran en la zona

“Perceptible” por el ser humano, y el 88.9 % de los datos se encuentran en la zona
“Imperceptible” por el ser humano. A continuación, se presenta otro gráfico de
resultados de monitoreo de vibraciones relacionados a la norma internacional de
percepción humana según “Steffens”.

27
 Figura 17 Muestra la relación de K Vs. Frecuencia (Hz).

De la figura se observa que el 33.3 % de los datos se encuentran en la zona

“Escasamente Perceptible” por el ser humano y el 66.7 % de los datos se
encuentran en la zona “Comienza a percibirse” por el ser humano.

28
V. APLICACIÓN DE LA PROGRAMACIÓN LINEAL EN

LAS OPERACIONES DE CARGUIO Y ACARREO.

5.1. Introducción
Campos y Ricra (2018) El objetivo de la optimización de procesos es ayudar a las
empresas a rediseñar sus procesos de negocio para reducir costos y aumentar la
eficiencia, logrando así los máximos beneficios mediante el uso de herramientas de
mejora adecuadas, cuyo uso es fundamental para lograr resultados efectivos. Un
proceso defectuoso conduce a sobrecostos y, si se trata de un "proceso central",
reduce la competitividad y, en última instancia, afecta la rentabilidad. La hoja de
cálculo EXCEL ® (Microsoft Corporation) con el complemento SOLVER puede
resolver problemas de optimización de procesos con base en un modelo matemático
de programación lineal llamado SIMPLEX y se utilizó como herramienta para
resolver el problema de optimización en el ciclo de extracción de carga de la mina
Lagunas Norte. El objetivo de este trabajo consiste en la aplicación de técnicas de
programación lineal SOLVER ® para garantizar un rendimiento óptimo del ciclo de
transporte de carga de minerales, lo que significa un flujo más dinámico y eficiente
de minerales al sistema de trituración. Nuestra hipótesis de optimizar el proceso de
carga-transporte de mineral en mina Laguna Norte mediante la aplicación del
software SOLVER® ha sido confirmada por los resultados obtenidos, las principales
conclusiones son: Mejora significativa de mineral y ahorros en costos operativos de
USD 18,533.00 durante la carga y descarga en mina Laguna Norte, 12,35% de los
costes del proceso. El uso de SOLVER ha demostrado ser una herramienta técnica y
útil para modelar y optimizar procesos mineros. SOLVER es parte de una familia de
comandos, a veces llamados herramientas de análisis, que encuentran el mejor valor
(mínimo o máximo) para una fórmula en una celda (llamada celda de destino) dentro
de las restricciones de otras celdas en la fórmula. Un límite o limitación de valor. en
la hoja de cálculo.
PALABRAS CLAVE: carga, carguío-acarreo, modelo matemático, programación
lineal, SOLVER, optimización.

5.2. Realidad Problemática

En línea con las tendencias globales, las empresas mineras de todo el mundo han
implementado estrategias significativas de reducción de costos, con disciplinas de
gestión de capital más estrictas que reemplazan las inversiones en proyectos,
mientras que las empresas mineras han simplificado carteras, renegociado la deuda
y cerrado operaciones marginales. En mi país, las empresas mineras continúan
enfrentando múltiples presiones, la calidad del mineral continúa decayendo, los
costos aumentan y los conflictos sociales se agravan aún más; dada esta realidad, las
mejoras incrementales ya no son suficientes. (Campos, 2018)
El cargue y transporte es el costo más alto (45%-65%) de esta operación minera, por
lo que es responsable de garantizar el mejor desempeño para todas las partes:
Material, equipos, instalaciones y la mano de obra. Usando métodos PL, podemos

29
 formular estrategias de optimización para estos procesos. Para ello utilizaremos una
hoja de cálculo EXCEL y el complemento SOLVER. (Ricra, 2018)
5.3. Formulación Del Problema
¿Cuál es el impacto de la Programación Lineal con el uso de SOLVER en la
optimización de las operaciones de carguío-acarreo de mineral en la Mina Lagunas
Norte, la Libertad?
5.4. Justificación
5.4.1. Justificación teórica
La programación lineal (PL) utilizada en procesos industriales ha demostrado ser
una herramienta muy poderosa y útil para la optimización y reducción de costos,
especialmente en la industria minera, y este estudio pretende facilitar su
incorporación y análisis en el proceso productivo minero empresarial.

5.4.2. Justificación práctica

El objetivo de este estudio es optimizar los recursos de las empresas mineras a
través de una gestión técnica más eficiente y tecnológico de las operaciones de
carguío y transporte de mineral.

5.4.3. Justificación valorativa

Las empresas mineras que han luchado durante años podrían lograr resultados
diferentes con nuevos enfoques para la gestión de operaciones, lo que podría
incluir invertir en herramientas como la programación lineal que crea valor para
la empresa e impulsa una mejor toma de decisiones.

5.4.4. Justificación académica

La finalidad del desarrollo del presente trabajo es ampliar nuestros conocimientos
teóricos en la aplicación de la programación lineal en las operaciones de carguío-
acarreo y como se aplica en nuestro campo laboral de ingeniería de minas.

5.5. Limitaciones
El desarrollo de este artículo sobre la optimización del carguío y acarreo en la mina
Laguna Norte corresponde a los procesos específicos de Laguna Norte utilizados en
sistemas de minería a cielo abierto y procesos de carguío y acarreo operados por
camiones grandes, cargadores frontales y palas.
5.6. Objetivos
5.6.1. Objetivo General
Determinar el Impacto de la programación lineal con el uso de SOLVER en la
optimización de las operaciones de carguío- acarreo de mineral en la mina lagunas
norte, la libertad, 2017.

30
 5.6.2. Objetivos Específicos

• Determinar el impacto de la decisión que se toma a través del programa

SOLVER en el desempeño de las palas hidráulicas y cargadores frontales.
• Conocer cuál es la influencia de la decisión que se toma a través del
programa SOLVER en el tiempo de carguío.
• Determinar cuál es el impacto de las restricciones identificadas por el
SOLVER en el desempeño de los camiones KOMATSU.
• Conocer el impacto de las restricciones identificadas por el SOLVER en el
tiempo de acarreo.

5.7. Bases Teóricas

5.7.1. Programación Lineal
La programación lineal utiliza un modelo matemático para describir el problema.
El adjetivo lineal significa que todas las funciones matemáticas del modelo deben
ser funciones lineales. En este caso, la palabra programación no se refiere aquí a
términos computacionales; en esencia es sinónimo de planeación. Por lo tanto, la
programación lineal involucra la planeación de actividades para obtener un
resultado óptimo; esto es, el resultado que mejor alcance la meta especificada de
acuerdo con el modelo matemático entre todas las alternativas factibles.
(Taha,2012).
5.7.2. Formulación de un modelo de Programación Lineal
En el proceso de formulación de un modelo de programación lineal hay que dar
los siguientes pasos:
Determinación de las variables de decisión.
Representan los elementos del sistema a modelar que son controlables por el
decisor. En los modelos lineales continuos estas variables toman como valores
números reales y se representan por letras con subíndices como se acostumbra a
hacer con las variables matemáticas, o literales alusivos a su significado: peso,
valor, etc. En el primer caso también se utiliza la representación como vector de
un conjunto indexado de variable:

𝑥 = (𝑥1 , 𝑥2 , … )
Determinación de las restricciones.
Representan las limitaciones prácticas de determinados recursos o imposiciones
físicas de la realidad. Se expresan como ecuaciones e inecuaciones lineales de las
variables de decisión. Matemáticamente adoptan una de las siguientes formas:

Formulación de la función objetivo.

31
 Se trata de la función que mide la calidad de la solución y que hay que optimizar
(maximizar un beneficio o minimizar un coste). También es una función lineal de
todas o parte de las variables de decisión.

Problema General de Programación Lineal

Para un caso de maximización según Verma (2010) es:

(Función objetivo)

(Restricciones explícitas)

(Restricciones de no-negatividad)

Donde:
Z=El valor de la medida global del rendimiento.
Xn =Niveles de actividad (1, 2…n); X1, X2, … Xn son variables de decisión.
ai,j=Cantidad de recursos i consumidos por cada unidad de actividad j.
bi=Cantidad de recursos que son disponibles para la asignación de actividades
(j=1, 2…m).

5.8. Descripción De Solver

5.8.1. Cargar el Complemento SOLVER
El complemento SOLVER es un programa de complemento de Microsoft Office
Excel que está disponible cuando se instala Microsoft Office o Excel, para usar el
complemento SOLVER, sin embargo, primero debe cargarlo en Excel.
1. En Excel Professional Plus 2016 Ir archivo > opciones
2. Haga clic en Complementos y, en el cuadro Administrar, seleccione
Complementos de Excel.
3. Haga clic en Ir.
4. En el cuadro Complementos disponibles, active la casilla Solver y, a
continuación, haga clic en Aceptar.

32
 a. Sugerencia Si el complemento Solver no aparece en el cuadro
Complementos disponibles, haga clic en Examinar para buscar el
complemento.
b. Si se le indica que el complemento Solver no está actualmente
instalado en su equipo, haga clic en Sí para instalarlo.
5. Una vez cargado el complemento Solver, el comando Solver está
disponible en el grupo análisis de la ficha datos.

Figura 18 Configuración de Solver

Fuente: Software Excel Solver v 16.0

5.9. Número De Ciclos

Ciclos efectivos de operación realizados por el camión minero. Operaciones de
Carguío-Acarreo
Operaciones que se llevan a cabo con el fin de retirar el material volado del frente y
transportarlo adecuadamente a su lugar de destino; la cual tiene la siguiente
secuencia:
1. Preparación de la zona de trabajo.
2. Posicionamiento de equipos.
3. Retirar el material volado desde el frente de trabajo.

33
 4. Traspaso del material al equipo de transporte dispuesto para el traslado.
5. Transporte del material a su lugar de destino (planta, acopio, botaderos, etc.).
6. Descarga del material.
7. Retorno del equipo de transporte al punto de carguío.
Ésta secuencia es cumplida hasta que haya sido retirado el material requerido del
frente. (Estrella, César 2016)
✓ Optimización
Optimización es el arte y la ciencia de determinar la óptima colocación de
recursos escasos y que conlleven al mejor efecto posible. Las técnicas de
optimización se ponen en juego todos los días en las cuestiones de
planificación, asignación de recursos, programación, toma de decisiones
industriales, etc. (Chinneck, J. 2000).

✓ Parámetro
Variable cuyo valor no está determinado por la persona a cargo de tomar
las decisiones, sino que se establece en forma exógena. (Eppen, G. 2000).

✓ Programación Lineal
Modelo determinístico formado por ecuaciones lineales, y que contiene
una sola medida de desempeño (función objetivo) por optimizar, sujeta a
satisfacer un conjunto dado de restricciones. (Eppen, G. 2000).

✓ Solver ®
Solver es un programa complementario de Microsoft Excel que se utiliza
para resolver problemas de Optimización. Se usa para encontrar un valor
óptimo (máximo o mínimo) para una fórmula en una celda - llamada célula
objetivo - sujeta a restricciones o límites, sobre los valores de otras celdas
de fórmula de una hoja de cálculo. (Microsoft Corporation© 2017).

5.10. Material Y Métodos

5.10.1. Unidad de estudio.

Reporte mensual de equipos de carguío - acarreo (cargadores frontales, palas y
camiones) utilizados en la Mina Lagunas Norte, la Libertad.

5.10.2. Población.
Reportes de la cantidad de equipos (09 equipos) de carguío-acarreo utilizados en
la mina Lagunas Norte la Libertad.

5.10.3. Muestra.
Se va a trabajar con 03 equipos (Cargador frontal, Pala y Camión) de carguío-
acarreo en el nivel 4160 de la mina Lagunas Norte la Libertad.

34
 5.10.4. Técnicas de recolección y análisis de datos
Las técnicas usadas en la presente investigación son la recolección de información
en el Internet, libros de Investigación de Operaciones, Análisis de Reportes
Técnicos etc.
5.10.5. Instrumentos
✓ Reportes Técnicos Independientes realizados al proyecto Lagunas Norte
✓ Memorias Anuales de Barrick Gold Corporation y otros.
✓ Revistas especializadas de la industria.
✓ Páginas web de Optimización y Programación Lineal.
✓ Página web de Barrick Gold Corporation.
✓ Uso del software EXCEL SOLVER® para la optimización del proceso
✓ La data recolectada se ordena, analiza y procesa a través de las
herramientas informáticas de procesamiento de textos y hojas de cálculo
Word® y Excel® respectivamente.

5.11. Procedimientos
Hacemos uso de la metodología de la Investigación de Operaciones:

Figura 19 Metodología de la Investigación de Operaciones

Fuente: www.gestiondeoperaciones.net

5.12. Ubicación
Minera Barrick Misquichilca S.A. (MBM) UEA Lagunas Norte se ubica en un área
aproximada de 185 Km2 del Alto Chicama, Distrito de Quiruvilca, Provincia de
Santiago de Chuco, Departamento de La Libertad en la parte Central Norte del Perú
aproximadamente a 90 km al este de la ciudad de Trujillo y a 11 km. al noreste del
pueblo de Quiruvilca. La mina está localizada a 7º57’ Latitud Sur y 78º15’ W
Longitud Oeste y se sitúa en la Cordillera Occidental de los Andes Peruanos a una
altitud aproximada de 4000 m.s.n.m. a 4260 m.s.n.m.

35
 Figura 20 Mapa de Ubicación de Lagunas Norte

Fuente: Google Images

5.13. Geología Local Y Mineralización

La mineralización de Lagunas Norte es del tipo epitermal de alta sulfuración, se
encuentra diseminada y hospedada en tufos y brechas volcánicas y sedimentarias. La
mineralización aflora en la superficie y ha sido definida sobre un área de 1600 m de
largo por 750 metros y hasta 300 metros de profundidad.
El 85% corresponde a óxidos y el 15% restante son sulfuros. Cerca del 80% de las
reservas están contenidas en los sedimentos Chimú y el 20% restante corresponde a
las rocas volcánicas del Grupo Calipuy.
El depósito consiste en un yacimiento de mineral de plata y oro, muy diseminado y
de baja ley, que se encuentra intercalado con lechos de carbón y lutitas carbonáceas,
así como areniscas, esquistos de barro y litologías volcánicas que contienen material
sulfuroso con contenido de oro. Aunque parte del carbón y de los esquistos de barro
carbonáceos contienen oro, su extracción no es práctica debido a factores técnicos y
económicos; por lo tanto, este material se almacena en forma separada para su uso
potencial en el futuro.

36
 Figura 21 Geología Local Lagunas Norte

Fuente: Barrick Gold Corporation – Lagunas Norte Technical Report NI-101

37
5.14. Clasificación Del Mineral

Figura 22 Clasificación del mineral Barrick Gold Corporation – Lagunas Norte

5.15. Optimización De Las Operaciones De Carguío Y Acarreo

En las tres zonas de mineralización de la mina (Alexa, Josefa y Dafne) obtenemos
mineralizaciones con contenidos de oro complejos y no todos los tipos son enviados
a chancadora antes de la lixiviación, algunos son enviados directamente y otras a
stock.

5.15.1. Equipos y Productividad

La productividad depende de factores como tipo de frente de minado, antigüedad
del equipo y pericia del operador principalmente. La mina cuenta con volquetes
de acarreo KOMATSU, modelo 730E, motor Diésel y capacidad de carga de 190
Tn.

5.15.2. Descripción del Proceso de Carguío y Acarreo

El proceso de envío de material a la chancadora se realiza generalmente con dos
equipos de carguío para cumplir con los planes de producción; estos son pala –
camión y cargador – camión.

38
 Figura 23 Sistema Tradicional de Carguío - Acarreo

Fuente: (Toribio Atencio, Poma Romero, Rodriguez Damas, & Villanera Perez, 2017)

Sistemas de carguío en Mina Lagunas Norte

Carguío con pala /cargador frontal

39
5.16. Base De Datos Tomados De Campo

Tabla 5 Registro de tiempos NIVEL 4160. PALA hacia la chancadora

40
Tabla 6 Registro de tiempos NIVEL 4070.CARGADOR hacia la chancadora

41
Tabla 7 Registro de tiempos NIVEL 4040. (CARGADOR) hacia chancadora

Tabla 8 Cuadro de producción de plan semanal

Tabla 9 Cuadro de tiempos Sistema Tradicional

42
Y considerando el costo operativo de cada camión de $250/h, desarrollamos la
siguiente. tabla tomando en cuenta el cuadro de tiempo del sistema tradicional.

Tabla 10 Costos de Acarreo

43
 VI. APLICACIÓN DE CADENAS DE MARKOV Y

SISTEMAS DE COLAS EN LA OPERACIÓN

UNITARIA DE ACARREO.

6.1. Introducción
Un ejemplo simple de un proceso estocástico es una sucesión de ensayos de
Bernoulli, por ejemplo, una sucesión de lanzamientos de un dado. Sin embargo, en
la mayoría de los procesos estocásticos, cada resultado depende de lo que sucedió en
etapas anteriores del proceso. El caso más simple de un proceso estocástico en que
los resultados dependen de otros, ocurre cuando el resultado en cada etapa sólo
depende del resultado de la etapa anterior y no de cualquiera de los resultados
previos.
Tal proceso se denomina proceso de Markov o Cadena de Markov. Finalmente, uno
de los objetivos de este informe es construir un modelo aplicado donde se seleccionen
los equipos de movimiento de tierra para una mina subterránea en particular, rampas
y camiones, estos, deben corresponder en función de sus características. La rampa
debe ser de un tamaño adecuado en relación con la altura y la anchura de los camiones
en los cuales se ira a cargar combustible. El número de camiones necesarios para
cumplir con los requisitos de producción y maximizar la eficiencia es difícil de
determinar, y el número de camiones necesarios cambiará con el tiempo a medida
que avanza la minería y las rutas de larga distancia se hacen más largas.
6.2. Objetivos
6.2.1. Objetivo General
Aplicar la teoría fundamental de cadenas de Markov y teoría de colas para
determinar el comportamiento de la flota de camiones en la carga de combustible
en una mina subterránea.

6.2.2. Objetivos Específicos

• Mostrar que el proceso es conjuntamente una cadena de Markov y líneas
de espera.
• Construir la matriz de transición.
• Mostrar que los estados son accesibles, se comunican, son recurrentes o
aperiódicos.
• Presentar las probabilidades de estado estable.
• Presentar los tiempos entre llegadas exponenciales y tiempos de servicios
exponenciales de la flota de camiones.

6.3. Caso de Aplicación

En la mina subterránea La Consentida se ha realizado un estudio demedición
de tiempos de transporte a su flota de camiones. La mina La Consentida es una mina
subterránea que explota vetas de cuarzo – oro y plata epitermales albergados en

44
 volcánicos y brechas. Típicamente los depósitos de mena exhiben una gran
orientación y varían de más de 24 m en espesor con buzamiento casi vertical. Este
mineral es explotado por el sistema de Corte y Relleno (Cut and Fill), por lo que a
medida que el mineral es extraído, el material estéril es introducido para soportar y
continuar su explotación ascendente por lo que requiere de una gran cantidad de
equipos de transporte para la operación.
La mina cuenta con 3 Bombas de abastecimiento de combustible que pueden ser
utilizadas a su vez por otros equipos como camiones aljibes, camionetas, Bobcat, etc.
Para un mejor orden se ha dispuesto que los diferentes equipos mineros carguen
combustible en determinados horarios. Disponiendo así que los equipos de transporte
carguen combustible al finalizar su turno. Pero dicha condición no siempre se cumple
alterando la disposición de las bombas para atender la demanda.
El camión puede cargar combustible en cualquiera de las Bombas. Para esto debe
situarse justo en frente de una y esperar el verde para ingresar a cargar combustible.
El situarse en una determinada bomba no implica el ser atendido por la misma.
El departamento de operaciones minas solicitó a su ingeniero realizar un estudio
sobre la probabilidad de que un camión se abastezca en una determinada Bomba. La
decisión fue entonces enviar por turnos camiones a la Bomba A, otro turno a la
Bomba B y otro turno a la Bomba C ya que el ciclo del problema es creciente a través
de los turnos. Se ordenó también que en el caso que el tiempo de espera se extienda
más de 15 minutos el camión se traslade a las Naves de mantención y quede pendiente
su abastecimiento de combustible para el siguiente Turno.

• Número de equipos = 10 camiones

• Tiempo medio de abastecimiento = 15 minutos
• Tasa de llegada = 1 camión / 6 minutos

Figura 22 Distribución de bombas de abastecimiento

45
Fuente: (Araya Alfaro, Córdova Becerra, Cavieres Liempí, & Jaime Jaime,
2014)
Las mediciones arrojaron los siguientes resultados:
• Al ir a la bomba A
o Probabilidad de abastecerse en A: 50%
o Probabilidad de abastecerse en B: 30%
o Probabilidad de abastecerse en C: 15%
o Probabilidad de retornar a Naves: 5%

• Al ir a la bomba B
o Probabilidad de abastecerse en A: 15%
o Probabilidad de abastecerse en B: 55%
o Probabilidad de abastecerse en C: 20%
o Probabilidad de retornar a Naves: 10%

• Al ir a la bomba C
o Probabilidad de abastecerse en A: 5%
o Probabilidad de abastecerse en B: 10%
o Probabilidad de abastecerse en C: 85%
o Probabilidad de retornar a Naves: 0%

46
• ESTADOS
E1: abastecimiento en bomba A
E2: abastecimiento en bomba B
E3: abastecimiento en bomba C
E4: retorno a nave

47
 • MATRIZ DE TRANSICIÓN

Determinamos la matriz de estados al cabo de 5 y 10 turnos y analizamos los

resultados obtenidos:

Tabla 29 Matriz de estados luego de 5 turnos

E1 E2 E3 E4

E1 0.1423 0.2321 0.3973 0.2283

E2 0.1233 0.2098 0.3927 0.2742

E3 0.1180 0.2072 0.5934 0.0814

E4 0 0 0 1

Fuente: (Araya Alfaro, Córdova Becerra, Cavieres Liempí, &

Jaime Jaime, 2014)

Tabla 30 Matriz de estados luego de 10 turnos

E1 E2 E3 E4

E1 0.0958 0.1641 0.3835 0.3567

E2 0.0898 0.1540 0.3644 0.3918

E3 0.1124 0.1938 0.4804 0.2135

E4 0 0 0 1

48
 Cuál es la probabilidad de que al cabo de 4 turnos un camión esperando en la bomba
A sea atendido en la bomba C.

Tabla 31 Matriz de estados luego de 4 turnos

E1 E2 E3 E4

E1 0.1684 0.2618 0.3761 0.1936

E2 0.1372 0.2372 0.3820 0.2436

E3 0.1128 0.2004 0.6310 0.0557

E4 0 0 0 1

Fuente: (Araya Alfaro, Córdova Becerra, Cavieres Liempí, & Jaime

Jaime, 2014)

Entonces la respuesta es 37.61% de probabilidad.

Porcentaje de utilización de las bombas al cabo de 4 turnos: 𝑓4 = (𝑃𝑇)4 ∗ 𝑓0

Tabla 32 Porcentaje de utilización desde el estado inicial y hasta el cuarto turno

𝑓0 𝑓1 𝑓2 𝑓3 𝑓4

E1 0.15 0.1350 0.1284 0.1241 0.1204

E2 0.20 0.2150 0.2160 0.2123 0.2070

E3 0.60 0.5725 0.5499 0.5299 0.5114

E4 0.05 0.0775 0.1058 0.1338 0.1612

49
 Tabla 33 Comportamiento del porcentaje de utilización para 10 turnos

Porcentaje de utilización para 10 turnos

0,7

0,
6

0,
5

0, f0 f1 f2 f3 f4 f5 f6 f7 f8 f9 f10
E1 E2 E3 E4

6.3.1. Análisis De Resultados

Al iniciar los turnos se debe tener en cuenta que hay un 16% de probabilidades de
tener que ir a abastecer equipos del turno anterior. Información que debe tener
presente el jefe de turno de operaciones de camiones. De esta forma sacar dicho
estado de “imprevistos” y agregarlos a “programados”
El tiempo de abandono de Bombas (retorno a naves) tiende a aumentar en el
tiempo por lo que es necesario exponer a los distintos departamentos de operación
el problema y plantear un mejor sistema de coordinación y control de tiempos y
horarios de abastecimiento.
Se recomienda a la administración de Bombas de Abastecimiento de Combustible
que: La tendencia de camiones es utilizar la Bomba C por lo que se recomienda
optimizar el tiempo de abastecimiento y maniobras, junto con mantener la Bomba
con suministro constante de combustible y así disminuir la cantidad de abandonos.

6.3.2. Sistema De Colas

Tabla 34 Datos de servidores y usuarios

50
 𝜇 4 Camión/hora (Tasa media del servicio)

𝜆 10 Camiones/hora (Tasa media de arribo al servicio)

𝑐 3 Cantidad de bombas (servidores)

𝑛 10 Número de camiones total (usuarios)

𝜆
2.5 Factor de utilización del sistema
𝜌=
𝜇

Fuente: (Araya Alfaro, Córdova Becerra, Cavieres Liempí, & Jaime Jaime, 2014)
Para el cálculo de los indicadores va a ser necesario hacer uso de las Siguientes
ecuaciones:

Tabla 35 Resultados obtenidos con la teoría de colas

𝑃0 4.49% Probabilidad que el sistema se encuentre en ocio

𝑃𝑛 3.26% Probabilidad de encontrar n clientes en el sistema

𝐿𝑠 6.0112 Cantidad esperada de camiones en el sistema

51
 𝐿𝑞 3.5112 Cantidad esperada de camiones en una cola

𝑊𝑠 0.6011 Tiempo de espera en el sistema (horas)

𝑊𝑞 0.3511 Tiempo de espera anticipado en cola (horas)

Fuente: (Araya Alfaro, Córdova Becerra, Cavieres Liempí, &

Jaime Jaime, 2014)

6.3.3. Análisis De Resultados En La Línea De Espera

En base a los resultados obtenidos se recomienda dividir en 2 el horario de

abastecimiento de camiones, al inicio y al término del turno debido ya que 3.5
camiones de la flota se encuentra en fila al momento de ir a cargar combustible.
Se recomienda también, en caso de persistir el problema, destinar más de 0.6011
horas (36.07 minutos) al abastecimiento, y de esta forma dejar sin pendientes al
turno inmediatamente siguiente. Si no se destina el tiempo necesario para el
abastecimiento se aumentará el número de retornos a naves.

6.3.4. Conclusiones
El desarrollo de este informe brinda soporte para asegurar estas dos técnicas que
nos brinda la Investigación de Operaciones, nos sirven para modelar el
comportamiento de una flota de camiones cuando se dirige a cargar combustible
a la rampa en minas subterráneas. Estos modelos son capaces de analizar sistemas
de transporte como existen en la actualidad y que se pueden utilizar para evaluar
la eficiencia de las operaciones sobre la base de sus tamaños actuales de la flota.
Esto también se puede combinar con datos de costos para los equipos en uso para
averiguar cuánto dinero se gasta para operar equipos que no están contribuyendo
directamente a la producción.
Los datos obtenidos fueron consistentes a una operación a una mina subterránea.
Tiempos entre llegadas exponenciales y tiempos de servicio exponenciales de la
mina que fue elegida, por lo que los supuestos del modelo son válidos para algunas
operaciones.
Para determinar un tamaño de flota óptima para un diseño de la mina dado y
configuración de carga sin correr una simulación completa, que puede ser más útil
el uso de modelos que implican simulación estocástica, para incorporar rutas de
larga distancia, tiempos de viajes y tamaño de las flotas. Esto permitiría a varios
tamaños y configuraciones de la flota a ser comparados sin tener que hacer
cambios del mundo real para adquirir entradas adicionales para el modelo, como
sería necesario para el modelo de cola. Este tipo de modelo estocástico también
podría ser utilizado en situaciones en las que la teoría de colas no es aplicable, por
ejemplo, en minas en las que las operaciones de carga no se ajustan a las
distribuciones necesarias para la teoría de colas para ser aplicado.

52
VII. OPTIMIZACIÓN DEL PROCESO DE CHANCADO

MEDIANTE PROGRAMACIÓN LINEAL EN UNA

PLANTA CONCENTRADORA DE COBRE

El proceso de chancado es una etapa crítica en la producción de concentrado de cobre en

una planta concentradora. La eficiencia de esta etapa afecta directamente la producción y
los costos operativos de la planta. En este trabajo se propone un modelo de programación
lineal para optimizar el proceso de chancado en una planta concentradora de cobre en
Perú.

En la metodología, se utilizaron datos de la planta concentradora de cobre, tales como la

capacidad de producción de la planta, la ley de cabeza del mineral, el tamaño de
alimentación, la distribución de tamaños del mineral, la tasa de reducción deseada, entre
otros. Estos datos se usaron para construir el modelo matemático de programación lineal
y para definir las restricciones del proceso de chancado, tales como la capacidad máxima
de la planta, la tasa de alimentación de mineral, el tamaño máximo de alimentación, la
tasa de reducción deseada, entre otras.

Las ecuaciones utilizadas en un modelo de programación lineal para la optimización de

chancado fueron:

53
 • Balance de masa en el alimentador: Qf = Qc + Qr, donde Qf es la tasa de
alimentación del mineral al alimentador, Qc es la tasa de descarga del alimentador
al chancador y Qr es la tasa de retorno del chancador al alimentador.
• Balance de masa en el chancador: Qc = Qp + Qd, donde Qc es la tasa de
alimentación del chancador, Qp es la tasa de producción de producto y Qd es la
tasa de descarga del chancador a la siguiente etapa del proceso.
• Restricción de tamaño máximo de alimentación: F80 <= Fmax, donde F80 es el
tamaño de partícula del mineral que pasa el 80% del total de la alimentación y
Fmax es el tamaño máximo de partícula que se puede alimentar al chancador.
• Restricción de tamaño de producto final: P80 <= Pmax, donde P80 es el tamaño
de partícula del producto que pasa el 80% del total del producto y Pmax es el
tamaño máximo de partícula del producto final.
• Restricción de tasa de reducción: R <= Rmax, donde R es la tasa de reducción
del chancador y Rmax es la tasa máxima de reducción permitida.

El modelo matemático se compone de varias ecuaciones de balance de masa, que

describen el flujo de material en la planta de chancado. Estas ecuaciones están sujetas a
las restricciones de la planta y a las especificaciones del producto final deseado, que se
definen en términos de la ley de cabeza del mineral y del tamaño de la partícula del
producto final.
Para resolver el modelo matemático, se utilizó el software de programación lineal
LINGO, que es un software de optimización muy utilizado en la industria. El software
permitió resolver el modelo en un tiempo razonable y obtener los resultados óptimos para
las variables de decisión del modelo, como la tasa de alimentación de mineral a cada etapa
de chancado, la configuración de los equipos de chancado, la ley de cabeza del mineral
en cada etapa, entre otras.

Los resultados del modelo mostraron que, al aplicar las configuraciones óptimas de
chancado, se podía aumentar la producción de concentrado de cobre en un 6% y reducir

54
los costos operativos en un 10%. Además, el modelo permitió identificar las variables
más importantes del proceso de chancado y establecer los puntos de operación óptimos
para cada una de ellas.
En las conclusiones del trabajo, se señala que la aplicación de herramientas de
investigación de operaciones, como la programación lineal, puede ser de gran utilidad en
la optimización del proceso de chancado en la minería. Además, se sugiere seguir
investigando y aplicando estas herramientas en otras etapas del proceso productivo
minero.

55