Técnico Científica

Predicción y estudio de vibraciones en

las preparaciones para la optimización
del minado Bench and Fill

Por: ñan las mallas de perforación durante la retardos, secuencia de iniciación, tipo de
Ing. Edgar Roque Ortiz, fase de minado, teniendo como variables detonador) que están relacionados con el
Ing. Edwin Cereceda Espinoza controlables: geométricas (diámetro, longi- efecto sísmico de la roca.
Compañía de Minas Buenaventura (2016) tud de carga, burden, espaciamiento); Fí-
sico-Químicas o explosivos (tipo de explo- El desarrollo de este trabajo está en-
sivo, energía, sistemas de cebado, VOD, focado a predecir y optimizar los efectos
carga operante); de Tiempo (tiempos de de la voladura en los costos de operación

Resumen
El método de explotación por Bench
and Fill tiene importantes ventajas como
menores costos y mejores rentabilidades. Abstract
Asimismo, entre los factores que caracte-
The Bench and Fill mining method has important advantages such as lower
rizan a este método son la dilución y la
costs and better profitability. Also, among the factors characterizing this me-
fragmentación con respecto a otros méto-
thod are the dilution and fragmentation, compared to other conventional me-
dos convencionales como corte y relleno,
thods such as cutting and filling, shrinkage, chambers and pillars, etc. This is
shrinkage, cámaras y pilares, etc. Es por
why special importance should be given to dilution and fragmentation control,
ello, que se le debe dar una importancia
also considering the variables generated during blasting. These controls are
especial al control de la dilución y frag-
normally carried out during Bench and Fill blasting, generalizing results for a
mentación, tomando en consideración va-
complete level, which are subject to trial and error, without considering certain
riables generadas en la voladura. Normal-
rock mass characteristics in different areas of that level. Most mining operations
mente estos controles se efectúan durante
perform the blasting study and redesign the drill pattern during the mining pha-
la voladura en Bench and Fill, generalizan-
se, having as controllable variables: geometric variables (diameter, loading
do los resultados para un nivel completo,
length, burden, spacing); physical-chemical or explosive variables (type of ex-
los cuales están sujetos a prueba y error,
plosive, power, priming systems, VOD, operating load); time variables (delay
sin considerar ciertas características del
times, initiation sequence, type of detonator), which are related to the seismic
macizo en diferentes zonas de ese nivel.
effect of rock.
La mayoría de las operaciones mineras
realizan el estudio de voladura y redise-

Figura 1. Muestra el centro volcánico Tambomayo, se trata de un estratovolcán que tuvo por lo menos dos conductos de emisión, y la cima del
cono debe haber superado los 6,700 m de altitud (Fuente: BVN, 2012).

MINERÍA la mejor puerta de acceso al sector minero MINERÍA / ABRIL 2019 / EDICIÓN 499 41
 n Minimización del impacto al macizo ro- cánicas peculiares para cada veta; Paola
coso con un rango de 0.40 m. con un RMR 33 y las rocas encajonantes
n Reducción de la voladura secundaria de un RMR 37, Mirtha con un RMR 50 y las
5% planteada originalmente en un 1%. rocas encajonantes un RMR de 42, se
optó por explotar estas vetas mediante el
Este trabajo también permitirá a los método de Bench and Fill, para minimizar
planificadores de minería, antes de dise- los costos y obtener una mejor rentabilidad
ñar, perforar y cargar un bloque de pro- del negocio.
ducción, evaluar los patrones de voladura
previstos en el futuro para minimizar el Dado el método de minado, el control
potencial de dilución de mineral con des- de la perforación y voladura es un pará-
monte y la probabilidad de dañar las es- metro clave que influye en la dilución e
Figura 2. Veta Mirtha Nv. 4740 (fuente: BVN). tructuras (roca encajonante). impacto hacia la roca encajonante.

(acarreo, transporte, relleno y chancado), Introducción La dilución y el impacto a la roca en-

evaluando el comportamiento del macizo La unidad de Tambomayo de Compa- cajonante van ligados de forma directa
rocoso con los parámetros obtenidos en ñía de Minas Buenaventura es un depósito ya que un mayor impacto en el macizo
los estudios de respuesta sísmica en la fase de oro y plata, ubicado a 4,800 msnm en rocoso repercute con el descaje de estas
de preparación de los tajeos (ventanas, el distrito de Tapay provincia de Caylloma y por consiguiente un incremento en la
cruceros, subniveles), zonificando los nive- a 271 Km al suroeste de la ciudad de Are- dilución, esta relación es una preocu-
les a minar de acuerdo a los parámetros quipa, considerado uno de los principales pación en todas las minas del mundo.
de respuesta sísmica: factor de velocidad proyectos de oro de 2016, con recursos Los ingenieros buscan reducir sus efec-
(k) y constante de atenuación (α), paráme- y reservas de 2.6M MT, 700k Au Oz y tos desde la parte inicial del proceso la
tros que son constantes por tipo de roca, 22.9M Ag Oz, las estructuras principales perforación y voladura con el sentido
explosivo y nivel, dando una característica del depósito corresponden a los sistemas de aumentar sus ganancias y reducir sus
especial a la U.E.A. Tambomayo de Com- de Mirtha y Paola, las cuales se encuen- costos.
pañía de Minas Buenaventura. tran cortando secuencias volcánicas com-
puestas principalmente por lavas, tobas y Objetivo
Las metodologías que se aplicaron tobas brechas de composición andesítica. n Generar y establecer un modelo predic-
para este estudio y evaluación son de dos tivo con el fin de determinar un nuevo
tipos: la primera referente al análisis del Las potencias de las vetas varían entre diseño de voladura para la aplicación
registro de vibraciones producto de las vo- 0.5 m y 15 m, con características geome- del método de minado Bench and Fill.
laduras en la fase de preparación, donde
se cuantificó su comportamiento a través
del macizo rocoso y su efecto sobre el
perímetro (entorno) y la fragmentación, to-
mando en consideración las propiedades
geomecánicas, estructurales y litológicas
del macizo rocoso, y la segunda relacio-
nada con la simulación e implementación
de voladuras en Bench and Fill, analizan-
do el efecto del impacto tanto al perímetro
como a labores e infraestructuras continuas

Gracias a los trabajos realizados, se

logró obtener los siguientes resultados:
n Una fragmentación con un P80 de 8
pulgadas.
Figura 3. Esquema de levantamiento para la zonificación de acuerdo a los parámetros de vibra-
n Se obtuvo una dilución promedio del
ción (fuente propia).
5%.
n Se logró zonificar cada nivel a minar
tomando como característica los pará- Tabla 1: Resumen de los Casos de Estudio. Validación del Método Propuesto
metros de vibración K y α propios del
macizo teniendo como característica el
tipo de roca de la unidad (lavas, tobas
propilíticas, tobas argilíticas, tobas bre-
chadas, arcillas).
n Predecir y controlar el impacto ocasio-
nado a las labores continuas permanen-
tes como Raise Climber 7, rampa 100
Tabla de resumen realizado por Onederra y Esen (2003), donde se obtiene k y α para distintos
positiva.
tipos de explosivo para un determinado tipo de roca.

42 MINERÍA / ABRIL 2019 / EDICIÓN 499 www.mineriaonline.com.pe

n Reducir el porcentaje de dilución a través troladas que repercuten en mayor costo de con el objeto de motivar a todos los in-
de diseños de la malla de perforación, sostenimiento, mayor plazo de ejecución de volucrados hacia mejores resultados, en
distribución y optimización del consumo la obra; así como la generación de daño por congruencia con la política de la mina.
de explosivos, basado en el estudio de vibraciones en estructuras adyacentes; más
vibraciones de la respuesta sísmica pro- aún en explotación por Bench and Fill donde Condiciones geológicas
ducto de la interacción explosivo – roca, es crucial controlar la sobre rotura para evitar Tambomayo es un depósito epitermal
con el fin de obtener resultados óptimos; altas diluciones y zonas inseguras de trabajo. de sulfuración intermedia vetiforme con mi-
menor daño, control de la estabilidad neralización de Au, Ag, Pb y Zn, siendo
de las cajas, zonas más estables y se- Alcance Au y Ag los principales elementos de inte-
guras cumpliendo con el D.S 055-2010- A todo el personal de la compañía rés económico.
EM, a un menor costo de minado. inmerso en la explotación de Bench and
Fill. Identificar los aspectos y parámetros En la base de la secuencia estratigrá-
Planteamiento del problema significativos, enfatizando aquellos que fica están los afloramientos del Mesozoi-
En muchas minas subterráneas se rea- pueden crear problemas para conjugar co (Grupo Yura). Estos afloran en la parte
lizan excavaciones con el empleo de ex- soluciones en el menor tiempo. sur-oeste de la región conformando los
plosivos que van del más al menos denso; flancos de los valles Andagua-Orcopam-
con mayores velocidades de detonación a Implementar y gestionar en forma in- pa y Colca. Estas unidades litológicas
menores; voladuras convencionales y no con- terna y externa los resultados obtenidos datan desde el Jurásico hasta el Cretáceo
y están constituidas por intercalaciones de
secuencias de areniscas, areniscas cuarzo-
Tabla 2. Propiedades Geomecánicas que se utilizarán para la Estimación
sas, calizas, calizas arenosas, delgados
de la Rotura para un Explosivo Determinado
horizontes de lutitas y limonitas.
Modulo
Tipo de Rcompresion Rtension Ratio Densidad
Rocas
Labor RMR
(Mpa) (Mpa)
de Young
Poisson (m.t./m3) Geología local
(Gpa) El centro volcánico Tambomayo se en-
Mineral SN 205 37 49.0 9.34 17.03 0.28 2.54 cuentra rodeado por otros centros volcá-
nicos más antiguos y más jóvenes, como
Mineral SN 206 37 49.5 9.28 17.03 0.28 2.54
son: hacia el lado oeste con la caldera
Mineral SN 986 37 50.0 8.16 17.03 0.27 2.54 Chinchón y Huayta (20 y 11 Ma). Al no-
Mineral SN 206 41 59.0 10.2 19.60 0.26 2.54 roeste con la caldera Caylloma (2.5 a 4.4
Ma). Al sureste el centro volcánico Mismi
T-And VNT 261 32 46.2 4.48 17.03 0.24 2.65
(1.5 a 1.0 Ma) y al sur con el centro volcá-
T-And CX 206 31 69.8 5.21 24.41 0.2 2.65 nico Huaca Hualca (1.0 a 0.5 Ma).
T-And VNT 263 32 69.8 7.98 24.41 0.2 2.65
T-And CX 205 35 49.0 5.21 17.03 0.19 2.65
El centro volcánico Tambomayo registra una
larga y compleja historia eruptiva; las secuen-
T-And SN 318 35 69.0 8.15 21.00 0.21 2.65 cias piroclásticas, efusivas o intrusivas muestran
Levantamiento de parámetros geomecánicos de la zona de estudio (Tajo Mirtha Nv. 4740 - diferentes eventos que gradan su composición
4772) (fuente: área Geomecánica – Vicor). desde traquiandesitas hasta dacitas.

Tabla 3. Parámetros del Tipo de la Zona de Estudio

Parametros SN 205 Nivel 4772
Clasificacion general Propiedades Locales
Zona Mirtha Resistencia a la Compresion Estatica (Mpa) 49
Tipo Volcánica Resistencia a la Compresion Dinamica (Mpa) 77
Localizacion Tapay Resistencia a la Tracción Estatica (Mpa) 9.34
Unidad Tambomayo Resistencia a la Tracción Dinamica (Mpa) 12.3
Propiedades Globales Estructura general y geologia insitu
Modulo de Young (Gpa) 17.03 Conjunto de juntas three joint sets 50
Ratio de Poisson 0.28 Geología T-And Mineral TBx-And
Densidad de la roca (m.t./m3) 2.54 Q(10-100) 43 JRC(10-100) 43
Indice de Fractura (%) 38 RMR(10-100) LOW 37 HICH 45
Velocidad de la Onda P(m/s) 3162.76
Evaluacion realizada de la primera zona junio-julio 2016, para el
Velocidad de la Onda S(m/s) 2635.63
minado del tajo piloto.
Velocidad pico particula critica (mm/s) 910.01

Tabla resumen de los parámetros geomecánicos y sísmicos del subnivel 4772 (fuente área de Geomecánica – Vicor).

MINERÍA la mejor puerta de acceso al sector minero MINERÍA / ABRIL 2019 / EDICIÓN 499 43
 Considerando la ubicación de las uni- vector máximo desplazamiento al tiempo (Equipo de sismografía Instantel), con el ob-
dades volcánicas respecto al centro volcáni- de la mineralización en un rango que va jetivo de determinar el impacto al entorno,
co Tambomayo y el grado de conservación desde NW/SW 30° a NW/SE 70° que análisis de frecuencias (software de simula-
de sus flancos del extinto cono volcánico indica una evolución tectónica del mismo ción Matlab – Simulink) para estudiar y de-
se puede asumir que se desarrolló entre el al tiempo de la mineralización. terminar el comportamiento de las frecuen-
Mioceno medio y Plioceno inferior. cias, todos estos estudios van determinados
Fuera de las estructuras geológicas a obtener las constates empíricas k y α.
Mineralogía pueden considerarse rocas intermedias a
Las estructuras del sistema de Mirtha y malas, como lo indica la alta frecuencia Los estudios se realizaron en la veta
Paola son formaciones epitermales recono- de discontinuidades por metro y la baja Mirtha, en el proceso de preparación de
cidas como un sistema de cuarzo-adularia, resistencia de la roca intacta. El UCS pro- subniveles con sección 3.5 m x 3.5 m, con
enriquecidas en oro, plata y metales base. medio del mineral masivo es de aproxima- puente de 12m de altura. Para una mayor
Los espacios abiertos fueron rellenados damente 99 MPa. La resistencia de las precisión se realizó el levantamiento en tor-
por texturas de bandeamiento simétrico, lavas y andesita es de aproximadamente no a los bloques a minar (Bench and Fill),
estructuras crustiforme, brechas con mine- 70 MPa, con módulos de Young bajos. ubicando los geófonos en los subniveles
rales de mena y ganga. de preparación y labores continuas.
Metodología
Distribución de vetas La metodología del estudio considera el Estimación de la velocidad pico
En el comportamiento estructural de la uso de modelos y simuladores de voladu- partícula
veta Mirtha y sus estructuras asociadas se ra (software JKSimBlast), levantamiento en La extracción del mineral requiere
puede ver que los esfuerzos presentes indi- campo de la respuesta sísmica de la roca realizar voladuras con explosivos y estas
can una extensión regional o dirección del frente a determinadas cargas explosivas producen vibraciones que se propagan

Figura 4. Zonificación de los puntos de estudio en la unidad Tambo- Figura 5. Esquema de levantamiento de voladuras tajo Mirtha, Paola
mayo de acuerdo a las características litológicas y geomecánicas del Nv. 4740 para la zonificación de acuerdo a los parámetros de vi-
macizo rocoso. (Fuente propia). bración (Fuente propia).

Tabla 4. Cálculo de la Ley de Atenuación SN 205 Zona 1, Nivel 4772

Datos de Monitoreo
Q (Kg) (Peso
Nº Turno Fecha Lecturas Labor D(m) De (m) VPP (mm/s)
max/Retardo)
1 Día 01/07/2016 2 SN 205 16.50 10.64 6.4 17.02
2 Día 02/07/2016 2 SN 205 16.50 10.64 6.4 18.90
3 Noche 01/07/2016 2 SN 205 19.00 11.70 6.8 15.30
4 Noche 02/07/2016 2 SN 205 19.00 11.70 6.8 14.90
5 Día 22/07/2016 2 SN 205 13.60 12.77 5.0 22.60
6 Día 22/07/2016 2 SN 205 13.60 12.77 5.0 23.12
7 Noche 29/07/2016 2 SN 205 17.20 12.77 6.0 21.04
8 Noche 30/07/2016 2 SN 205 17.20 12.77 6.0 19.54
9 Día 04/08/2016 2 SN 205 17.80 12.10 6.4 17.23
10 Día 06/08/2016 2 SN 205 17.80 12.10 6.4 16.90
11 Noche 14/08/2016 2 SN 205 19.65 4.20 11.7 6.80
12 Noche 16/08/2016 2 SN 205 19.65 4.15 11.8 7.34

44 MINERÍA / ABRIL 2019 / EDICIÓN 499 www.mineriaonline.com.pe

por el macizo rocoso. La velocidad pico veles de vibración generados a diferentes lejano), para estimar las vibraciones pro-
partícula (VPP), a distintas frecuencias es distancias y para diversas cantidades del ducidas por la voladura.
el parámetro utilizado internacionalmente mismo explosivo. Las constantes K, α, β,
para predecir el nivel de daño y los efec- van a depender de la geología del terreno Modelo que permite conocer los com-
tos causados por la voladura. y del tipo de explosivo a utilizar. ponentes vertical, radial y transversal de la
velocidad de partícula, que pueden ser de-
Existen modelos que estiman el com- Una de las ventajas de determinar la terminadas a través de la siguiente fórmula:
portamiento de la velocidad de partícu- velocidad pico partícula producto de las
la bajo ciertas condiciones de voladura. voladuras es la posibilidad de desarrollar VPP = K(D⁄√Q) α (2)
Estos relacionan tres variables: velocidad un criterio de daño y grado de fragmenta-
máxima de partícula (VP,mm/s), carga ción que se ajuste a cada tipo de trabajo K: factor de amplitud, α: factor de ate-
máxima por retardo (Q, Kg) y distancia del en particular, desde el punto de vista de nuación.
punto de interés a la voladura (D,m). las características del macizo rocoso, de
los explosivos utilizados, los diseños apli-
VP = K(D⁄Q α)-β (1) cados y el nivel de daño generado.

La expresión (D⁄Q α) es conocida como Para este trabajo utilizaremos el mo-

distancia escalar y busca relacionar los ni- delo de Devine 1966 (modelo a campo

Tabla 5. Resultados de las Constante k y α Obtenidos de 68 Voladuras

#
Labor Nivel Voladuras Zona Veta K α RMR Litología
Analizadas
Figura 6. Modelo ley de atenuación ecu. (1),
SN 205 4772 12 Zona 1 Mirtha 399.00 1.10 47 Min VPPc=240.08(D⁄√Q)1.21, con un nivel de corre-
SN 206 4760 8 Zona 1 Mirtha 384.00 1.31 47 Min lación de 0.9769 en la zona 1, nivel 4772.
Este modelo corresponde a una muestra del
SN 986 4785 4 Zona 1 Mirtha 404.00 1.12 47 Min total que se realizó (fuente propia).
SN 206 4760 7 Zona 1 Mirtha 399.00 1.14 51 Min
VNT 261 4760 4 Zona 2 Paola 218.00 1.05 32 T-And
CX 206 4760 9 Zona 2 Paola 218.00 1.05 31 T-And
VNT 263 4772 4 Zona 2 Paola 196.00 1.14 32 T-And
CX 205 4772 14 Zona 3 Paola 244.00 1.14 35 T-And
SN 318 4760 6 Zona 3 Mirtha 249.00 1.21 35 T-And

Valores de k y α para una terminada zona, Mirtha zona 1 rango para k: 350-400, α: 1.00-1.5.
Las voladuras se realizaron durante el periodo junio-octubre de 2016, obteniéndose 136 lecturas
con el equipo Instantel de doble geófono.

De acuerdo a la Tabla 4, se pudo dividir en tres zonas el tajo a minar, tomando en considera- Figura 7. Esquema de zonificación de acuerdo
ción como punto de partida los valores k y α, el tipo de roca y la litología del macizo rocoso a los parámetros k y α. Nivel 4740 – 4772
(fuente propia). (fuente propia).

Tabla 6. Límites de Vibración (velocidad pico partícula crítica) por Zona

Módulo de Velocidad de
Zona Veta Labor Nivel Litología RMR RC (Mpa) VPPc (mm/s)
Young (Gpa) Onda P (m/s)
Zona 1 Mirtha SN 205 4772 Min 47 17.03 49.0 3162.76 1049.00
Zona 1 Mirtha SN 206 4760 Min 47 17.03 49.5 3162.76 1059.80
Zona 1 Mirtha SN 986 4785 Min 47 17.03 50.0 3162.76 1070.50
Zona 1 Mirtha SN 206 4760 Min 51 19.60 59.0 3357.33 1155.45
Zona 2 Paola VNT 261 4760 T-And 32 17.03 46.2 2914.97 790.79
Zona 2 Paola CX 206 4760 T-And 31 24.41 69.8 2880.21 823.59
Zona 2 Paola VNT 263 4772 T-And 32 24.41 69.8 2914.97 833.53
Zona 3 Paola CX 205 4772 T-And 35 17.03 49.0 3068.20 877.35
Zona 3 Mirtha SN 318 4760 T-And 35 21.00 69.0 3068.20 1008.12
Valores de velocidad pico partícula crítica (VPPc), propio de la roca, para un determinado tipo de explosivo, levantamiento realizado en los niveles
4740-4772 (fuente propia)

MINERÍA la mejor puerta de acceso al sector minero MINERÍA / ABRIL 2019 / EDICIÓN 499 45
 Para una predicción de la zona a mi- Tabla 7. Cuadro de Densidades de Velocidad Ef: energía requerida para fragmentar
nar es necesario obtener las constantes Pico Partícula (VPP) (MJ), V: Volumen (m3), σt: resistencia a la
K(i,j) y α(i,j), donde i representa a una de- tracción (MPa), E: módulo de Young (MPa).
terminada zona y j un nivel, con otros dos
parámetros, el coeficiente de correlación y Energía consumida como efecto
la desviación estándar que proporcionan sísmico
la precisión del ajuste. Un coeficiente de La energía necesaria requerida para
correlación cercano a uno indica que el fe- vencer la resistencia de la roca para ser
nómeno corresponde al modelo aplicado. fragmentada es considerablemente menor
a la energía entregada por los explosivos
Velocidad pico partícula (VPP) industriales, este porcentaje de energía
Una partícula dentro del macizo roco- adicional se consume como vibración.
so está limitada a un determinado movi-
miento y velocidad, de acuerdo a sus ca- Para determinar la energía es necesa-
racterísticas geomecánicas y litológicas, rio registrar la amplitud, frecuencia domi-
cuando la velocidad pico partícula supera nante, duración de la vibración, distancia
este límite se genera la creación de nue- de la fuente y las características de la
vas fracturas, la velocidad pico partícula roca:
crítica es un indicador para medir el nivel
de daño ocasionado por las voladuras y Es = 124.d.ρr.Vs.a.f.t v (6)
determinar una distancia límite de acuerdo
a la carga operante. Con el cuadro se pudo evaluar el impacto Donde:
en cada zona a determinadas distancias, y d : distancia al punto de monitoreo (m).
σt×Vp cuánto influye una determinada cantidad de
ρr: densidad roca (Kg/m3).
VPPc = explosivo, en la fase de preparación de subni-
E (3) veles (fuente propia). Vs: velocidad del sonido en la roca (m/s).
a: amplitud de la vibración (m).
Donde σt es la resistencia a la trac- f: frecuencia de la vibración (Hz).
ción, Vp velocidad de onda p, E módulo confinamiento del explosivo, y la energía tv: duración de la vibración (s).
de Young. de explosión del dispositivo a utilizar.
Onederra y Esen (2003) realizaron
Selección del material explosivo Etr = FI.Fc.En (4) un estudio donde hallan las constantes de
Para la obtención de una adecuada Holmberg & Persson (K y α) para simular y
fragmentación debemos considerar las Etr: energía transferida a la roca (MJ). modelar la velocidad pico partícula (PPV)
dos variables controlables como es la se- FI: Factor de impedancia Fc: Factor de a campo cercano para un determinado
lección del explosivo y el secuenciamien- confinamiento En: Energía de explosión tipo de explosivo. El método analítico pro-
to. Para la selección de la carga explosiva (MJ). puesto reduce la instrumentación y el moni-
se tomó como base la energía necesaria toreo necesario para estimar los halos de
para fragmentar un determinado volumen Energía requerida para impacto producto de la voladura. La vali-
de roca, y a la vez generar el menor daño fragmentar la roca dación se dio en cuatro casos de estudio.
posible. Cada roca tiene propiedades mecáni- En la Tabla 1, se muestra un resumen del
cas particulares, como la resistencia a la trabajo descrito.
Energía transferida a la roca tracción y módulo de elasticidad que de-
La transferencia de energía al macizo terminan la cantidad de energía que es ne- Varios autores han realizado una inte-
rocoso está en función de las característi- cesaria para generar la rotura del macizo. gración de las ondas y frecuencias pro-
cas del explosivo y la roca, tomando en σt×V venientes desde cada zona de la carga,
consideración la impedancia, el grado de Et = (5) estimando la diferencia de tiempos en que
2E
Tabla 9. Distancias Óptimas vs Esquemas
Tabla 8: Modelo de Atenuación para Cada Zona a Minar de Carguío

Zona Distancias Óptimas

Zona 1 D= 1.17
VPP/396.5 × Q

Zona 2 D= 1.08
VPP/210.67 × Q

Zona 3 D= 1.18
VPP⁄246.50 × Q

Modelo para calcular las distancias óptimas

De acuerdo a los estudios se pudo elaborar un modelo de atenuación, con el objetivo de ser uti-
de impacto para una determinada cantidad
lizado para el cálculo y predicción de diseños de voladura en Bench and Fill. Para la obtención
de explosivo y una velocidad pico partícula
de estos valores se consideró el promedio de los k y α por zona (fuente propia).
crítica.

46 MINERÍA / ABRIL 2019 / EDICIÓN 499 www.mineriaonline.com.pe

las ondas llegan a un punto P. Sin embar- mografía Instantel, Vibracord. Las variables Establecimiento de los modelos
go, la velocidad pico partícula el punto de carga operante, consumo de explosivo de vibraciones para el minado
máximo de estrés en la roca no ocurre y retardos se obtienen del carguío del fren- con Bench and Fill
cuando las primeras ondas producto de la te. Estudios realizados en las labores de La determinación de los valores k y
voladura llegan al punto P. Esto acontece preparación (subniveles, cruceros), de los α, realizada a partir de nuestros ensayos
cuando las ondas se dispersan a través tajos para el minado con taladros largos, de junio a octubre de 2016, nos permitió
del macizo rocoso y se encuentran con desde el nivel 4740 hasta el nivel 4840. establecer correlaciones mucho mejores
una cara libre, lo cual genera una onda entre los niveles de vibraciones medidos
de tensión reflejada. Obtención de las constantes K y α a la proximidad de la labor y los criterios
del modelo de atenuación de impacto hacia el macizo. Estos pará-
Entonces, las variaciones de tiempo Con el levantamiento de los equipos metros nos permiten realizar la simulación
entre las llegadas de ondas son tan cortas, de sismografía se obtuvo los parámetros y aplicación de voladuras por Bench and
que simplemente no se consideran como que se muestran en la Tabla 4. Fill, para la obtención del burden y espa-
un problema para los cálculos. ciamiento óptimo de acuerdo al tipo de
Obtención de los límites de explosivo a utilizar en la búsqueda de la
Propiedades geomecánicas de la vibración fragmentación óptima.
roca (unidad Tambomayo) De acuerdo al levantamiento realizado
Se realizó el levantamiento en los sub- en campo y al cálculo de los parámetros Diseño de mallas de perforación
niveles de los parámetros geomecánicos para la simulación de la voladura, factor para Bench and Fill
necesarios para el estudio, con apoyo del de velocidad (K) y factor de atenuación Los parámetros de vibración obteni-
área de Geomecánica y la consultora Vi- (α), para la comparación y el cálculo de dos por cada zona, nos permiten selec-
cor. la distancia a la cual se genera daño, es cionar un adecuado tipo de explosivo,
necesario calcular el límite de vibración secuenciamiento, un esquema de carguío
Desarrollo del trabajo definido por la velocidad pico partícula adecuado, con el objetivo de obtener
El levantamiento en campo de las labo- crítica del macizo. una fragmentación planeada, controlar la
res de preparación brinda la información
necesaria para estimar las ecuaciones de
los modelos antes descritos, estos valores
han sido obtenidos con los equipos de sis-

Tabla 10. Simulación de Densidades de

Velocidad Pico Partícula para Diseños de
Voladura en Bench and Fill

Figura 8. Curvas de respuesta sísmica, para velocidades pico partícula (VPP) y distancias varia-
bles, manteniendo una determinada carga de explosivo constante (fuente propia).

Con los modelos obtenidos se calcula el daño

que puede producir una voladura de Bench
and Fill a una determinada distancia, de
acuerdo al cálculo se tiene un daño a 0.60
m para una carga operante de 12 kg, y una
VPPc de 1,500 mm/s (fuente propia).

Tabla 11. Parámetros para Diseño de Malla

Q D B S
K α
Kg m m m

350 1.31 12 0.6 1.28 1.29 Figura 9. Simulación, comportamiento de la densidad de frecuencias de voladuras simultáneas
para 3 y 6 filas. Software Matlab (fuente propia).

MINERÍA la mejor puerta de acceso al sector minero MINERÍA / ABRIL 2019 / EDICIÓN 499 47
dilución y, por ende, reducir el daño al Asimismo, para una mayor precisión se Con estos parámetros-rangos de vibra-
macizo rocoso. realizó el análisis de frecuencias, simulan- ciones, y utilizando las ecuaciones de con-
do distintas situaciones. Para este caso se minución, logramos obtener un diseño de
Se realizó la simulación y análisis de la tomó el ejemplo de ver cuantas filas son malla adecuado para un determinado tra-
zona 1, para su posterior implementación, las adecuadas a disparar. mo a minar por el método de Bench and
manteniendo constante los parámetros de Fill, obteniendo para la primera zona a mi-
vibración k= 350 y α= 1.31, variando El análisis del rango de densidades de nar, con la emulsión Emulex 80 1 1/2x12,
las distancias y la carga operante, con el frecuencias, donde se relaciona la con- un burden de 1.28 m y un espaciamiento
objetivo de optimizar la malla de perfora- centración de frecuencias versus un deter- de 1.29 m.
ción, para esa determinada zona. minado tiempo, para distintas filas a volar
en Bench and Fill. Resultados
De acuerdo con la Tabla 10, el rango Con los estudios y análisis desarrolla-
óptimo, está en el intervalo de 11.1 kg, a Se realizó el análisis comparativo del dos en la primera parte determinamos el
13.8 kg de carga operante a una distan- comportamiento de las frecuencias con el diseño de carga, utilizándose el modelo
cia promedio de 0.60 m. uso del Emulex 80 1 1/2x12, con el mis- de conminución para la obtención del bur-
mo tipo de secuenciamiento, y la variación den y espaciamiento adecuado, para pos-
Con estos rangos podemos deducir de de la carga operante, obteniendo para la teriormente ponerlo a prueba en campo
cuánto debe ser el burden y espaciamien- voladura de 3 filas un promedio de 150 con el fin de contrastar los datos obtenidos
to óptimo, para un tipo de explosivo. Hz y para 6 filas un promedio de 180 Hz. de los parámetros empíricos k y α.

Tabla 12. Diseño de Carga. Indicadores de Perforación

y Voladura Bench and Fill. Fila 15-16

Figura 10. Izquierda, diseño y esquema de carguío para Bench and Fill,
con aplicación de Air Deck. Derecha, simulación del comportamiento de
los halos de energía haciendo uso del software Jk Simblast (fuente propia). (Fuente propia).

Figura 11. Zona perturbada por los disparos realizados en los 20 m

de abertura del tajo Mirtha 11 nivel 4740. Distancia de daño 0.8 m Figura 12. Análisis de fragmentación tajo Mirtha 11 – nivel 4740 (fuen-
(fuente propia). te propia).

48 MINERÍA / ABRIL 2019 / EDICIÓN 499 www.mineriaonline.com.pe

Figura 13. Levantamiento con equipo escáner, para la obtención de la Figura 14. Comportamiento de los indicadores (factor de potencia Kg/
dilución y recuperación (fuente propia). Ton versus fragmentación P80). Por fila disparada (fuente propia).

De acuerdo al estudio con un k Kg/Ton, el caso inverso es el incremento braciones, que ha permitido mejorar la
= 350 y α = 1.31 se ha obtenido un del P80 hasta lograr 5.5 pulgadas. comprensión y el control del daño oca-
burden de 1.28 m y espaciamiento de sionado por las voladuras, la eficacia
1.29 m para la veta Mirtha con una Conclusiones de este estudio depende de objetivos
condicionante de carga operante de 1. Se ha obtenido un modelo predictivo claros y un sólido apoyo de gestión en
12 kg, considerando como material ex- para cada zona a minar, para los nive- la obtención de equipos de instrumenta-
plosivo a usar la emulsión, el objetivo les 4740-4772. ción, conocimientos en la instalación de
trazado fue: instrumentos, recopilación de datos, con
n Un factor de potencia de 0.45 Kg/Ton. técnicas detalladas para análisis y así
n Una fragmentación (P80) de 6 pulga- poder obtener las variables buscadas.
das.
n Daño al macizo de 0.60 m. 8. Estos estudios son difíciles de implemen-
tar en algunas minas por la falta de
Para llegar a estos objetivos se im- equipos y de un área de perforación y
plementó el uso de cámaras de aire (Air voladura, pero se demuestra en el es-
Deck), para reducir la carga explosiva en tudio que se puede controlar mejor el
la columna y, por ende, el impacto a las efecto de la voladura sobre la fragmen-
rocas encajonantes. tación, dilución y daño al macizo.

Los primeros disparos se desarrollaron 2. Se ha logrado una dilución de 0% y una Recomendaciones

en la veta Mirtha Tajo 11 nivel 4760, con recuperación del 90%, se está trabajan- n Se recomienda realizar más estudios a
una abertura de 20 metros, con un tiempo do en nuevos diseños para incrementar detalle considerando otros tipos de ex-
de auto sostenimiento de 45 días. la recuperación en base al radio de in- plosivos y más parámetros geomecáni-
fluencia del explosivo. cos – estructurales.
El resultado después de las pruebas de
Referencias
16 filas es: 3. El daño a la roca encajonante está en
n Factor de potencia de 0.42 Kg/Ton. promedio de 0.8 m. Centro de Innovación Tecnológica de Explosivos
Enaex (Cintex). 2006). Manual de Monitoreo de
n Una fragmentación de 5.5 pulgadas.
Vibraciones Generadas por Tronaduras, Análisis
n Un daño al macizo de 0.80 m. 4. El factor de potencia obtenido es 0.42
y Modelamiento.
Kg/Ton con una fragmentación (P80) de
Valores cercanos a los objetivos trazados. 5.5 pulgadas. Jack de la Vergne. McIntosh Engineering, North
Bay Ontario, Hard Rock Miner's Handbook.
Se ha obtenido en los últimos dispa- 5. Los valores obtenidos k y α son propios Ministerio de Industria y Energía – Secretaría de
ros una fragmentación con un P80 de 5.5 de un determinado tramo de bloque a la Energía y Recursos Minerales. 1985. Estudio
pulgadas. minar sujeto a un tipo de explosivo. de la influencia del seccionado de cargas en
un mismo barreno para disminuir el nivel de
Los valores obtenidos en la dilución y 6. El uso de un modelo predictivo permite vibraciones. España.
recuperación mediante el levantamiento reducir el riesgo de una voladura defi- Pablo García, Martín. 2014. Análisis del
con escáner son: ciente en el método de minado Bench contenido en frecuencias de vibraciones por
n Dilución: 0%. and Fill. voladura, Departamento de Ingeniería Química
n Recuperación: 90%. y Combustibles.
7. En la unidad de Tambomayo, se logró Vicor Ingenieros 15-002-02. 2015. Estudio
Se observa como el factor de potencia obtener un modelo predictivo con un Geomecánico del Proyecto Minero Tambomayo,
ha disminuido, hasta haber obtenido 0.42 programa exitoso de monitoreo de vi- Compañía de Minas Buenaventura.

MINERÍA la mejor puerta de acceso al sector minero MINERÍA / ABRIL 2019 / EDICIÓN 499 49