Aliaga Aw
Aliaga Aw
Aliaga Aw
TESIS
“REDUCCIÓN DEL DAÑO AL MACIZO ROCOSO
CIRCUNDANTE OCASIONADO POR LA VOLADURA
DE ROCAS, USANDO EL MONITOREO,
MODELAMIENTO Y ANÁLISIS DE VIBRACIONES,
CASO MINA UCHUCCHACUA”
ELABORADO POR
WILLAM LIONEL ALIAGA ALIAGA
ASESOR
M.Sc. ING. CARMEN ROSALÍA MATOS AVALOS
LIMA-PERÚ
2016
DEDICATORIA
A MI PADRE, ROGER.
ALMUDENA
AGRADECIMIENTO
Agradezco a todas las personas e instituciones, en especial la empresa EXSA S.A., que
hicieron posible la recopilación de los datos, con la instrumentación adecuada, para
la elaboración del presente trabajo.
Un especial agradecimiento a mis asesores, M.Sc. Ing. Carmen Rosalía Matos Avalos
y M.Sc. Ing. Fidel Julio Hidalgo Mendieta, por contribuir con sus enseñanzas y
observaciones para con el presente trabajo.
iii
ÍNDICE GENERAL
DEDICATORIA 1
AGRADECIMIENTO 2
RESUMEN 12
ABSTRACT 13
ANTECEDENTES REFERENCIALES 14
OBJETIVOS 16
HIPÓTESIS 16
METODOLOGÍA 17
INTRODUCCIÓN 19
2.3.2.3 Sostenimiento 91
2.3.2.4 Relleno 91
CONCLUSIONES 121
RECOMENDACIONES 123
BIBLIOGRAFÍA 124
ANEXOS 127
viii
RELACIÓN DE FIGURAS
Figura 1.1: Ciclo de esfuerzo sobre la roca, compresión seguido de tensión ........... 23
Figura 1.9: Onda generada por la detonación de una taladro cargado con anfo ....... 40
Figura 1.13: Instrumentación para medir las vibraciones de una voladura .............. 50
Figura 1.14: Acelerómetro piezoeléctrico ............................................................. 52
Figura 1.19: Forma geométrica del modelo de Holmberg y Persson, near-field ...... 69
Figura 1.21: Registro de las tres componentes de una forma de onda producto de la
voladura de una labor de desarrollo en mina Uchucchacua .................................... 76
Figura 1.24: Posiciones relativas de los puntos de registro (Wiss y Linehan) .......... 79
Figura 3.4: Ley de atenuación en el Crucero 6685 con anfo (Examón Q) ............... 99
Figura 3.5: Ley de atenuación en el Crucero 6685 con anfo (Examón P) ...............100
x
Figura 3.10: Malla de perforación y voladura para el carguío con emulsión ..........111
Figura 3.12: Perforación de los taladros de la corona (09 taladros en total) ...........113
RELACIÓN DE TABLAS
Tabla 1.1: Valores de las constantes 𝐾𝐾 y 𝛼𝛼 para distintos tipo de roca .................... 68
Tabla 1.2: Velocidad pico de partícula crítica para distintos tipo de roca ................ 74
Tabla 3.7: Velocidad pico de partícula crítica de las labores monitoreadas ............101
Tabla 3.9: Criterio de daño del modelo inicial versus el modelo replanteado .........115
Tabla 3.11: Criterio de daño Examón Q versus Examón P - Crucero 6685 ............116
Según el MEM entre los años 2000 y 2013 se han producido 837 accidentes mortales
de minas apliquen técnicas científicas de vanguardia para crear zonas de trabajo más
de vibraciones para determinar las leyes de atenuación que gobiernan cada zona de
estudio. Con los valores de las contantes de la roca para cada tipo de explosivo se
formularán modelos matemáticos que nos simularán el daño en zonas próximas a las
perforación y voladura para cada tipo de roca, tomando especial énfasis en el diseño
According to MEM between 2000 and 2013 there have been 837 fatal accidents en the
reduce the damage caused mainly by rock blasting and in this way create work areas
determine the attenuation laws that governing each study area. With the values of the
constants of the rock mass for each type of used explosive will be made mathematical
models that will simulate us the damage in close areas to the excavation and they will
Our vibrations study will allow us to redesign the mesh drilling and blasting
parameters for each rock type, taking particular emphasis on the design contour
blasting. Our results will be compared with the vibrations study of the year 2012,
no se han realizado más estudios hasta mediados del año 2014 donde la compañía
circundante en la zona baja de la mina (entre los niveles 4060 y 3780), llamada mina
Socorro.
PLANTEAMIENTO DEL PROBLEMA
En una operación minera uno de los grandes problemas con el que se tiene que lidiar
a diario son las eficiencias de los avances, ya que de éste principalmente depende el
Pero es importante detenernos y analizar el daño causado por las vibraciones después
de cada voladura, para determinar un límite permisible de vibración que nos permita
Por ende nos planteamos la pregunta: ¿Es posible determinar un modelo matemático
de vibraciones que nos permita predecir y reducir el daño al macizo rocoso circundante
deben de considerar para llevarlo a cabo y sobre qué dominio se puede aplicar?
Devine (1962) y Holmberg y Persson (1978). Para nuestro caso, por las limitaciones
El siguiente estudio se realizó entre los meses de mayo y noviembre del año 2014 en
1. Objetivo General:
análisis de vibraciones.
2. Objetivos Específicos:
vibraciones,
HIPÓTESIS
1. Hipótesis General:
circundante.
2. Hipótesis Específicas:
cercanas a la excavación,
• Con la reducción del daño al macizo rocoso, se obtendrán labores más seguras y
METODOLOGÍA
1. Diseño metodológico:
Se realizará una tesis del tipo experimental, siguiendo los siguientes pasos:
2. Instrumentos:
3. Población:
secciones de 4 m x 4 m y que son clasificadas geo mecánicamente como roca tipo III,
ubicadas entre los niveles 4060 y 3780 de la zona baja, denominada mina Socorro.
4. Muestra:
Para nuestro estudio se definió una muestra de tres labores: Galería 661, Crucero 6685
y Rampa7046.
INTRODUCCIÓN
los operadores. En la actualidad y cuando las exigencias ambientales son cada vez más
aplicación de excesiva carga explosiva, que puede ir más allá de lo requerido para
fragmentar la roca.
voladura superficial o subterránea. En este plano deben ser definidas una adecuada
de rocas, son efectivamente aprovechados para fragmentar la roca. Por lo que la mayor
De esos efectos, son cuatro los que destacan: Vibraciones transmitidas a los macizos
causa más común de las preocupaciones y de protestas de las personas afectadas en las
vecindades de los trabajos de voladura. En algunos casos estos reclamos pueden ser
21
disparo.
niveles de las vibraciones, las que pueden ser modificadas alterando la carga explosiva,
Las vibraciones son un movimiento cíclico que ocurre dentro de un medio, debido al
paso de fases alternativas de compresión y tensión. Con respecto a las vibraciones por
Consideremos un caso simple, de un solo taladro de voladura con una pequeña carga
mayor o menor medida como un resorte, una vez que la fuerza de compresión es
removida, la roca se relaja y vuelve a su estado inicial, pasando más allá del mismo
Como la roca se mueve más allá de su posición original, esto crea una fase de tensión
Ya que la roca se comporta como un resorte al paso de las ondas de vibraciones (el
módulo de Young representa la rigidez del resorte), el ritmo al cual la roca se relaja
puede ser diferente al ritmo al cual es sometida por un pulso de corta duración. En este
caso, no es inusual ver en rocas débiles con módulos de Young relativamente bajos,
La ecuación que explica las características del movimiento de una partícula por el paso
Donde:
La frecuencia de las vibraciones indica el número de veces por segundo que la onda
de propagación pasa por un ciclo completo de compresión y tensión. El factor que tiene
una gran influencia sobre las frecuencias de las vibraciones son las dimensiones de las
cargas, columnas grande de carga tienden a producir bajas frecuencias. Sin embargo,
otros factores incluyen los módulos de la roca y la razón de carga producida por la
cortas distancias. Es importante saber que una onda con una frecuencia única, y que se
El hecho que las ondas registradas a grandes distancias tengan menores frecuencias a
contienen un amplio rango de frecuencias, y que las ondas de altas frecuencias son
continua, en la cual cada ciclo de propagación tiene la misma forma, un valor único es
Las unidades de amplitud dependen del tipo de sensor utilizado para detectar el paso
cualquiera de éstas, teóricamente podría permitir el cálculo de las otras dos. Los
dispositivos más baratos y fáciles de usar para medir las vibraciones son los geófonos,
considerada como el mejor indicador del esfuerzo inducido en el macizo rocoso, y por
voladura y la distancia del punto de monitoreo a la voladura. Para asegurar que el valor
máximo de la velocidad de vibraciones generado por una voladura sea registrado y que
un completo detalle de las ondas de vibraciones, y un tiempo después del paso de las
cuerpo. A 500 m la onda de vibración puede ser de 500 a 1000 ms más larga que la
duración de la voladura.
La longitud de onda, 𝜆𝜆, se puede calcular a partir de una onda de vibración con na
frecuencia única, 𝑓𝑓, es decir una onda armónica simple, por la fórmula 𝜆𝜆 = 𝑉𝑉𝑃𝑃 /𝑓𝑓,
través de la roca. Esta velocidad puede ser medida utilizando dos geófonos ubicados a
permitir un cálculo más preciso. En la Figura 1.5, los geófonos están separados 300 m, y
usando la ecuación simple, 𝑉𝑉𝑃𝑃 = 𝑠𝑠/∆𝑡𝑡, donde 𝑠𝑠 es la distancia entre los geófonos y ∆𝑡𝑡
29
La mayoría de las rocas tienen una velocidad de propagación entre 3,000 m/s y 5,000
m/s. Mediciones de la velocidad de propagación en roca menores que 1,500 m/s son
medición antes de aceptar un valor tan bajo, ya que ellas implican un grado de
los 6,000 m/s son consideradas también poco confiables, y nuevamente los sistemas
tan altos.
La onda de choque generada por la detonación de cargas explosivas crea tensiones que
todos los sentidos) y transfiere una energía vibracional al macizo rocoso que es
a través de ondas.
Los tres tipos de onda que se pueden observar cuando se monitorean las vibraciones
como ondas de cuerpo. Las ondas de superficie, sin embargo, se transmiten muy
que 1.5 a 2 veces su longitud de onda (aunque esto pueda representar algunos cientos
de metros).
de propagación. El término primario se origina en que esta onda tiene una gran
produce un segundo tipo de onda que corta o tiende a cambiar la forma del
ondas 𝑆𝑆 viajan sólo en los sólidos, puesto que su existencia depende de módulos
Ambas ondas, 𝑃𝑃 y 𝑆𝑆, viajan en trayectorias esféricas desde el cráter, a través del
cuerpo de los materiales, por dicha razón a esta clase de ondas se les denomina
ondas de cuerpo.
de las ondas 𝑃𝑃 y 𝑆𝑆. Cuando las ondas de cuerpo alcanzan la superficie de la tierra,
La primera onda que llega al monitor será siempre la onda 𝑃𝑃, ya que ésta, de todos los
32
tipos de onda, es la que viaja más rápido (entre 30 y 50 % más rápido que la onda 𝑆𝑆).
Sin embargo ésta puede no ser la onda de mayor amplitud, de tal manera que no
siempre es fácil de identificar. La Figura 1.6 presenta dos ondas de vibraciones, donde
onda, donde se hace difícil determinar con claridad el tiempo de arribo y su amplitud.
La primera onda en llegar será siempre la onda 𝑃𝑃 la cual viaja una distancia más corta
su amplitud puede ser muy pequeña, y otra onda que realice un recorrido más largo a
33
través de una roca más competente llegará con un pequeño atraso y con una amplitud
cuando la amplitud del primer arribo es muy baja y difícil de identificar. Cuando
utilizados para medir vibraciones son los geófonos, los cuales miden la velocidad del
utilizan geófonos, y por lo tanto producen una señal de voltaje que varía con el tiempo,
Si la onda de vibración es continua y armónica simple, (es decir una única frecuencia
explosivo por taladro, ya que esto controla la fuente de energía que generan
posteriormente las vibraciones. El otro factor que tiene una gran influencia en la
ya que esto controla cuan rápidamente se disipa la energía. La atenuación del macizo
(las grietas generan una gran pérdida de energía) y el módulo de la roca. Rocas con un
bajo módulo tienden a presentar una deformación plástica, resultando en una pérdida
de energía, mientras que rocas con un alto módulo se comportan de una manera elástica
vibración es continua y armónica simple (es decir una única frecuencia sinusoidal, 𝑓𝑓),
muestra a continuación:
𝑎𝑎0
𝑣𝑣(𝑡𝑡 ) = � 𝑎𝑎(𝑡𝑡 ). 𝑑𝑑𝑑𝑑 = 𝑠𝑠𝑠𝑠𝑠𝑠(2𝜋𝜋𝜋𝜋𝜋𝜋) (Ecuación 1.6)
2𝜋𝜋𝜋𝜋
36
𝑠𝑠𝑠𝑠𝑠𝑠(2𝜋𝜋𝜋𝜋𝜋𝜋 ) = 1). Hay que notar que el valor pico de la aceleración estará desfasado
con respecto a la velocidad pico, es decir, el valor pico de la aceleración ocurre cuando
aceleración, pero fuera de fase con la velocidad de la partícula, es decir el valor pico
aceleración.
𝑎𝑎0
𝑠𝑠 (𝑡𝑡 ) = � 𝑎𝑎 (𝑡𝑡 ). 𝑑𝑑𝑑𝑑 = − cos(2𝜋𝜋𝜋𝜋𝜋𝜋 ) (Ecuación 1.7)
4𝜋𝜋 2𝑓𝑓 2
Cercano a la voladura, los niveles de aceleración son muy altos, y usualmente mayores
son bastante realistas cuando las mediciones son hechas muy cercanas al frente de
El fundamento para usar el valor pico de velocidad de partículas, VPP, es debido a que
la mayoría de los análisis consideran que el daño, o potencial daño, está relacionado
velocidad de vibración tiene una correlación directa con el esfuerzo, cuando la onda
está viajando como una onda plana, y la deformación es elástica, tal como muestra la
37
siguiente ecuación:
𝑉𝑉𝑉𝑉𝑉𝑉 𝜎𝜎
𝜀𝜀 = = (Ecuación 1.8)
𝑉𝑉𝑃𝑃 𝐸𝐸
𝜎𝜎𝑚𝑚á𝑥𝑥 𝑉𝑉𝑃𝑃
𝑃𝑃𝑃𝑃𝑃𝑃𝑚𝑚á𝑥𝑥 = 𝜀𝜀𝑚𝑚á𝑥𝑥 𝑉𝑉𝑃𝑃 = (Ecuación 1.9)
𝐸𝐸
Puesto que las vibraciones viajan con una componente sinusoidal de compresión y
esfuerzo a la tensión. Este valor es difícil de medir, por lo tanto, es normal estimar la
indirecta tal como el Índice de Resistencia Brasileño. Como resultado de esto, uno
puede estimar la velocidad de partícula que probablemente causará una ruptura por
conservadores (es decir bajos valores de VPP) y pueden ser más apropiados en zonas
donde la estabilidad es crítica y donde existe incertidumbre acerca del valor real de la
resistencia a la tensión.
LAS VOLADURAS
pero no pueden ser bien descritas por una ecuación armónica única que involucre una
de 120 taladros.
Figura 1.8: Registro de onda de vibración de una voladura de producción, se muestra la variación
de la velocidad de partícula respecto al tiempo
Fuente: Blastware Help File
mejor empezar con la onda producida por un solo taladro (onda elemental), y proceder
a partir de ese punto a un análisis de la onda de vibración de una voladura que contiene
más simple que puede ser analizada en una operación normal de voladura; sin
embargo, puede ser interesante comparar la onda de un único taladro de voladura con
la onda generada por una carga esférica (una carga de longitud igual al diámetro). En
voladuras que tienen muchos taladros, se considera generalmente que cada taladro
genera una onda elemental con las mismas características, y que se encuentran
separadas en el tiempo por los retardos y en el espacio por las dimensiones geométricas
de la malla de perforación.
El pulso tiene una forma única, pero a pesar de que ésta es una oscilación natural no
se asemeja a una onda sinusoidal, siendo más representativa a una forma de onda no
compresión (+) y la tensión (-) de la onda pueda indicar la diferencia entre el grado de
Figura 1.9: Onda generada por la detonación de una taladro cargado con anfo
Fuente: Blastware Help FileFrecuencia de las Vibraciones
Es claro que el pulso de vibración de un único taladro no puede ser descrito usando
efectivo para el análisis de frecuencia, aunque cuando se utiliza este método es muy
difícil representar la frecuencia dominante por un único valor, puesto que se observa
de compresión (360 mm/s) es mucho más alto que el esfuerzo de tensión (115 mm/s).
considera que el gran esfuerzo de compresión cercano a los taladros es la causante del
intenso quebrantamiento alrededor del taladro, mientras en zonas alejadas del taladro
Es importante indicar que el área bajo la fase de compresión debe ser igual al área bajo
dirección más que en la otra. Para revisar si las áreas son iguales, se debe integrar la
aproximadamente 4 ms. El tiempo adicional es aquel para que el macizo rocoso vuelva
roca, pero también por el tiempo de retención de los gases de detonación en el taladro.
Se puede ver que la duración de los pulsos de vibración serán más y más largos a
ondas de superficie.
Los retardos son usados en las voladuras de manera de asegurar que las vibraciones
solo taladro. A pesar de los retardos usados, es inevitable algún grado de interferencia
43
sitio específico, normalmente el valor pico es el único dato requerido. Sin embargo, en
términos del daño potencial que la onda de vibración pueda causar a las estructuras en
una valor importante, como por ejemplo una onda con una amplitud constante de 100
mm/s de 2 s de duración, se esperaría que tenga una daño mayor que una onda con un
limitado número de mediciones, y algunas veces sin medición alguna que nos muestre
una idea básica de la forma de la onda. Esto requiere tener una comprensión razonable
vez más grande. Sin embargo, la energía total no puede aumentar y en un medio
frente de onda esférico. Por lo tanto, la energía por unidad de área debe decrecer,
debido a que la energía total es irradiada en una superficie total en aumento. Para un
distancia de propagación, por lo tanto la energía por unidad de área debe de crecer
𝐸𝐸 1 1
∝ − = − (Ecuación 1.11)
𝐴𝐴 𝐴𝐴 4𝜋𝜋𝑟𝑟 2
expansión es esférico. Para una forma general del frente de onda, se espera que la
45
velocidad de vibración a cierta distancia 𝑟𝑟, 𝑉𝑉(𝑟𝑟) , esté relacionada con la velocidad 𝑉𝑉0 a
energía por fricción significa que la energía de vibración por unidad de área no se
La pérdida friccional ocurre debido a que las partículas están sometidas a esfuerzos de
compresión y tensión, y una suposición común es que la energía perdida por ciclo de
∆𝐸𝐸
− = 𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 (Ecuación 1.13)
𝐸𝐸
se refieren a esta constante como 𝑄𝑄−1 , donde 𝑄𝑄 está referido al factor de pérdida de
46
energía. En rocas competentes, 𝑄𝑄 puede tener un valor de 70, es decir, 1/70 ó 1.4 % se
blandas, 𝑄𝑄 puede ser 10, es decir, se pierde un 10% de la energía por ciclo.
como:
𝜋𝜋𝜋𝜋 𝜋𝜋𝜋𝜋𝜋𝜋
− −
𝑉𝑉 (𝑟𝑟) = 𝑉𝑉0𝑒𝑒 𝑄𝑄𝑄𝑄 = 𝑉𝑉0 𝑒𝑒 𝑄𝑄𝑉𝑉 𝑃𝑃 (Ecuación 1.15)
Donde:
100 metros, una onda de 10 Hz pasará 0.25 ciclos en una roca con una velocidad de
onda 𝑃𝑃 de 4,000 m/s, pero a 100 Hz la onda pasará a través de ella 2.5 ciclos. Por lo
tanto, las ondas de frecuencias altas perderán más energía por unidad de distancia que
GEOMÉTRICA Y FRICCIÓNAL
Para tener un sentido general del significado de este efecto, es mejor considerar el
de 100 Hz, en dos tipos de roca extremas, una dura casi elástica con 𝑄𝑄 = 100 y 𝑉𝑉𝑃𝑃 =
5,000 m/s, y una fracturada y plástica con 𝑄𝑄 = 10 y 𝑉𝑉𝑃𝑃 = 2,000 m/s. La amplitud en la
roca dura decrecerá a 0.94 veces su valor inicial, es decir, perderá aproximadamente
Ahora consideremos en la roca fracturada, la onda pasará los 100 m con 5 ciclos, y la
amplitud de la onda decrecerá a 0.21 veces su valor inicial, es decir, perderá 79% de
48
término de la distancia por la raíz cuadrada del peso del explosivo. Sobre distancias
observar que los parámetros de la ecuación son diferentes para mediciones de campo
En la Figura 1.12, se muestran tres ondas de 10 Hz, 50 Hz y 250 Hz, con la misma
10 Hz.
amplitud de la onda.
La instrumentación que se utiliza para medir las vibraciones (Figura 1.13) inducidas
análisis de la información.
tiene directa relación con el tipo de sensor utilizado y la técnica empleada para su
instalación y orientación.
Los dos tipos básicos de sensores usados para medir las vibraciones del terreno son:
1.4.1.1 Acelerómetros
• Su unidad de medida es 𝑔𝑔 (𝑔𝑔 = 9.81 𝑚𝑚/𝑠𝑠 2), con rangos de 0 hasta 250 𝑔𝑔,
52
1.4.1.2 Geófonos
la velocidad de partícula de la onda. Esta tensión es producida por una bobina móvil
ubicada dentro del campo de un imán fijo (Figura 1.15). Son los sensores más
datos antes mencionados. Como regla general se puede establecer que resulta
inadecuado usar un transductor de velocidad del tipo bobina móvil, cuando las
En primera instancia el equipo de vibraciones debe ser tal que los geófonos asociados
sensibilidad es la adecuada. Para tal efecto se debe conocer los niveles máximos
presentan a continuación:
los de campo cercano y los de campo lejano, los cuales poseen un rango de
infinito, por lo que los geófonos poseen una característica que se refiere a su
generadora de cada vibración debe ser conocida. Se requiere para esto una
en estudios de modelamiento.
cada uno de los cuales puede ser amplificado o atenuado. Después de esta
principales:
Minimate PlusTM (Figura 1.13); pero en el mercado existen una variedad de equipos
adecuado para el análisis de la forma de onda, y donde las principales tareas que se
• Cursor móvil sobre la señal para un análisis acucioso de los tiempos y las
amplitudes,
• Filtro de frecuencia,
En nuestro caso el análisis se realizó con el software Blastware III, compatible con el
de ocurrencia,
función directa con la carga del explosivo (𝑄𝑄) y en función inversa de la distancia (𝐷𝐷)
60
Donde:
𝐾𝐾, 𝛼𝛼 y 𝛽𝛽: Son constantes empíricas que dependen de la geología del terreno, la
y la secuencia de la voladura.
Por lo que es importante tener en cuenta que para cada zona de estudio los valores de
Los modelos mostrados a continuación son sólo referencias utilizadas para poder
Una de las primeras ecuaciones de propagación fue sugerida por Morris (1950) y
obedece a la expresión:
�𝑄𝑄
𝐴𝐴 = 𝐾𝐾. (Ecuación 1.18)
𝐷𝐷
61
Donde:
𝐾𝐾: Es la constante característica del lugar que varía desde 0.57 a 3.40.
Leconte (1967), en una revisión de las técnicas de control de las vibraciones, sugiere
siguiente:
�𝑄𝑄
𝑉𝑉 = 𝐾𝐾𝑣𝑣 . (Ecuación 1.19)
𝐷𝐷
Entre los trabajos posteriores más rigurosos destacan los de Blair y Duvall (1954),
Suponiendo una simetría esférica de la carga explosiva, la conclusión fue que cualquier
dimensión lineal debe de ser corregida por la raíz cúbica de la carga del explosivo.
(1971).
la siguiente ley:
62
𝐷𝐷 −𝛼𝛼
𝑉𝑉 = 𝐾𝐾. � 1/3� (Ecuación 1.20)
𝑄𝑄
Donde:
las distancias deben de ser corregidas por la raíz cuadrada de la carga, Devine (1962),
𝐷𝐷 −𝛼𝛼
𝑉𝑉 = 𝐾𝐾. � 1/2� (Ecuación 1.21)
𝑄𝑄
Esta expresión ha sido una de las más usadas hasta la actualidad por numerosos
Otros autores como Atewel et al (1965), Holmberg y Persson (1978), Shoop y Daemen
general:
Por lo que para estimar la vibración resultante en puntos muy próximos a la carga, se
63
𝑄𝑄 1/2
𝑉𝑉 = 𝐾𝐾. � � (Ecuación 1.23)
𝐷𝐷𝐷𝐷 3/2
Lundborg (1971), basándose en datos del U.S. Bureau of Mines (Nicholls et al, 1971),
log 𝑉𝑉 = 4.08 + 0.14 log 𝑄𝑄 − 2.06 log 𝐷𝐷 + 0.22 log 𝑄𝑄 . log 𝐷𝐷 (Ecuación 1.24)
Las investigaciones efectuadas en los últimos años han permitido realizar una mejor
Aunque la caída exponencial fue desde hace tiempo reconocida, e.g. Duvall y Petkof
𝑄𝑄 −1 −𝛼𝛼(𝐷𝐷/𝑄𝑄 1/3 )
𝑉𝑉 = 𝐾𝐾. � � . 𝑒𝑒 (Ecuación 1.26)
𝐷𝐷1/3
Suponen que las ondas internas son las que predominan y que existe divergencia
esférica.
En nuestro caso los modelos matemáticos que utilizaremos para modelar nuestra ley
Cuando nos referimos a Campo Lejano nos estamos refiriendo a puntos de monitoreo
carga del explosivo (𝐷𝐷 > 3𝐿𝐿𝐿𝐿), donde 𝐿𝐿𝐿𝐿 representa la longitud de la carga del
a aquellos puntos alejados a las voladuras, donde la carga máxima por retardo se
mundial por las empresas fabricadoras de explosivos, dicho modelo está expresado por
la Ecuación 1.27.
−𝛼𝛼
𝐷𝐷
𝑉𝑉𝑉𝑉𝑉𝑉 = 𝐾𝐾. � � (Ecuación 1.27)
�𝑄𝑄
65
Donde:
𝐾𝐾 y 𝛼𝛼: Son las contantes empíricas que dependen de la características del macizo.
El término (𝐷𝐷/�𝑄𝑄) es lo que se conoce como la distancia escalar (𝐷𝐷𝐷𝐷), de éste modo
Esta ecuación es la que utilizaremos para predecir el nivel de vibración y que nos
A pesar de que éste modelo tiene buenos resultados en el campo lejano, presenta las
siguientes restricciones:
estructuras predominantes,
• Los valores de las constantes 𝐾𝐾 y 𝛼𝛼 son muy variables para distintos tipos de
estos datos desde una mina a otra o desde un sector a otro en la misma mina
partícula en zonas muy cercanas a la voladura, debido que considera la carga explosiva
Por esta razón, Holmberg y Persson (1979) introdujeron un modelo modificado, que
al contrario del modelo de campo lejano, considera que la velocidad pico de partícula
densidad lineal constante, sobre toda la longitud de carga; remplazando de esta manera
el término de la carga operante del modelo de Devine por el de densidad lineal de carga
explosiva.
campo cercano como un modelo aceptable para predecir la vibración en distancias muy
Donde:
En la Figura 1.19 se muestra con más claridad todos los parámetros geométricos que
−𝛼𝛼
𝑙𝑙 𝐻𝐻 + 𝑥𝑥𝑠𝑠 − 𝑥𝑥0 𝑥𝑥0 − 𝑥𝑥𝑠𝑠
𝑉𝑉𝑉𝑉𝑉𝑉 = 𝐾𝐾 � �tan−1 � � + tan−1 � ��� (Ecuación 1.30)
𝑟𝑟0 𝑟𝑟0 𝑟𝑟0
experimentalmente con el modelo de Devine, por las limitaciones que se tiene con el
determinación de las constantes específicas del lugar en el que se vaya a disparar. Estas
vibraciones generadas por el disparo de taladros cargados con cierta mezcla explosiva.
1995); (Villaescusa et al., 1997); (García Bermúdez y López Jimeno, 1999); se refleja
en la Tabla 1.1.
distancia que debería ser aplicada entre la fuente de vibración y el punto de monitoreo.
del centro de la voladura puede ser designado como el punto de la fuente, ya que la
distancia entre cargas es mucho menor que la distancia desde la fuente hasta el punto
elegido para determinar el modelo de regresión. Estas distancias incluyen las distancias
69
frecuencia en el collar de carga. La intuición nos sugerirá que la distancia que más
adquiridos lo más cercano posible a una carga explosiva típica, con geófonos o
Por otra parte, en el proceso de análisis se ajustan los datos de vibraciones al modelo
de campo cercano propuesto por Holmberg y Persson, que considera la carga explosiva
distancia a la carga, longitud del taco, densidad lineal del explosivo, etc., a diferencia
del modelo tradicional propuesto por Devine que utiliza una sola distancia para el total
Figura 1.20: Comparación entre los modelos de Devine (far-field) y Holmberg y Persson (near-
field)
Fuente: “Curso de Monitoreo, Interpretación y Modelamiento de Vibraciones”,
SCHERPENISSE CARLOS, Codelco Chile, 1999
zona más cercana a la fuente donde se aprecia como el modelo de Devine sobre estima
de una definición. Los primeros que enunciaron definiciones más precisas al respecto
siguiente:
de revoque,
ni la estabilidad de la estructura.
golpe de una puerta puede causar deformaciones de hasta 0.00014 ", equivalente a la
deformación causada por una vibración de 12.5 mm/s. El daño asociado a la voladura
que degrade la calidad geotécnica del macizo rocoso y que afecta el rendimiento de las
operaciones. Puede ser generado por un mal diseño de la voladura o una mala
72
en la roca una vez que se supera la velocidad crítica, extensión y apertura de fracturas
vecindad inmediata a las cargas explosivas, donde el impacto de las voladura es más
intenso y los niveles de esfuerzos inducidos son similares a los esfuerzos necesarios
𝑉𝑉𝑉𝑉𝑉𝑉
𝜀𝜀 = (Ecuación 1.31)
𝑉𝑉𝑃𝑃
deformación inducida 𝜀𝜀, para una roca con una velocidad de onda de compresión 𝑉𝑉𝑃𝑃 .
73
Esta ecuación supone una estabilidad lineal del material a través del cual la vibración
se está propagando y permite una estimación razonable para la relación entre la roca
soportada por la roca antes de que ocurra la falla por tracción, puede ser estimada en
𝜎𝜎𝜏𝜏 . 𝑉𝑉𝑃𝑃
𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 = (Ecuación 1.32)
𝐸𝐸𝑑𝑑
Donde:
expresada en 𝑀𝑀𝑀𝑀𝑀𝑀,
macizo rocoso en cuestión y que las ondas son del tipo armónicas planas.
• Hay extensión de fracturas preexistentes si 𝑽𝑽𝑽𝑽𝑽𝑽𝑽𝑽 > 𝑽𝑽𝑽𝑽𝑽𝑽 ≥ 𝟎𝟎. 𝟐𝟐𝟐𝟐 ∗ 𝑽𝑽𝑽𝑽𝑽𝑽𝑽𝑽.
ASP Blastronic en función a los datos proporcionados por diferentes minas, empleó la
74
partícula crítica para algunos tipos de roca más frecuentes (Tabla 1.2), donde los
estimación y deben constantemente ajustarse con mediciones in-situ, que den cuenta
cuantitativamente del daño y que permitan establecer con mayor precisión la capacidad
Tabla 1.2: Velocidad pico de partícula crítica para distintos tipo de roca
Velocidad Pico de
TIPO DE Resistencia a la Velociad de la Módulo de
Partícula Crítica
ROCA Tracción (Mpa) Onda P (m/s) Young (Gpa)
(mm/s)
Cuarzo-1 14.00 5102.00 55.50 1286.00
Brecha-1 7.30 4298.00 30.20 1037.00
Milonita 1.90 2940.00 14.60 380.00
Ox-1 6.70 4373.00 35.20 836.00
Ox-2 7.20 4804.00 44.40 774.00
Cuarzo-2 8.40 4207.00 42.70 823.00
Brecha-2 3.90 4041.00 39.30 401.00
Andesita-1 14.90 4975.00 67.30 1100.00
Diorita 13.20 4650.00 48.60 1260.00
Brecha-3 11.30 4650.00 58.30 900.00
Pórfido-1 6.70 3829.00 31.20 823.00
Pórfido-2 5.10 3661.00 46.60 401.00
Fuente: Manual de Vibraciones ASP Blastronics
desplazamiento máximo y duración del evento. Existe una gran variedad de criterios
MONITOREO DE VIBRACIONES
Estas medidas deben describir de la mejor forma el evento de vibraciones, para lo cual
El registro obtenido entrega una onda producto de la vibración, medida en sus tres
Onda Longitudinal
Onda Vertical
Onda Transversal
Tiempo (ms)
Figura 1.21: Registro de las tres componentes de una forma de onda producto de la voladura de
una labor de desarrollo en mina Uchucchacua
VIBRACIONES
A partir de estos registros se puede obtener información del nivel de vibraciones que
Sin embargo a este nivel de vibraciones que se identifique, se debe verificar que el
utilizando la función Integrate sobre la onda de velocidad, la cual entrega los valores
de desplazamiento de la onda completa, tal y como muestran las Figuras 1.22 y 1.23.
77
La magnitud de las vibraciones terrestres y aéreas en un punto varía según la carga del
La primera etapa del análisis de la señal es determinar qué carga representa cada
Una detonación con baja eficiencia puede producirse por diferentes causas, algunas de
las cuales son: Tiros con burden y/o espaciamiento muy pequeño, cargas con algún
grado de insensibilización, proyección del material del taco, etc. En estos casos se
Si una lectura del sismógrafo nos muestra picos de velocidad de partícula bastante
que origina que la energía del explosivo no sea utilizada eficientemente en fragmentar
la roca, sino al contrario, se disipe en forma de vibración a través del macizo rocoso
circundante.
diferencia de los tiempos de llegada de los pulsos generados por la detonación de los
La Figura 1.24 muestra el caso de una fila de taladros con diferentes posiciones
𝑆𝑆. cos ∅
𝑡𝑡𝑒𝑒 = 𝑡𝑡𝑛𝑛 − (Ecuación 1.33)
𝑉𝑉𝑉𝑉
Donde:
en grados sexagesimales.
80
Figura 1.25: Direcciones preferentes de cooperación de ondas en una voladura múltiple (Wiss y
Linehan)
Fuente: Manual de Perforación y Voladura de Rocas, Carlos López Jimeno, Emilio López Jimeno
y Pilar García Bermúdez, 2003
El ángulo crítico de la posición relativa donde las ondas llegan al mismo tiempo y, por
lo tanto, puede producirse una cooperación entre dichas ondas sísmicas, será aquel
𝑆𝑆. 𝑡𝑡𝑛𝑛
∅ = cos −1 � � (Ecuación 1.34)
𝑉𝑉𝑉𝑉
En la Figura 1.25 se representa una voladura múltiple y las direcciones para las que
existe una mayor probabilidad de refuerzo o cooperación de las ondas de acuerdo con
constructivas o con efectos sumatorios, en los primeros estudios realizados por Duvall
Langerfors (1963) señala que con intervalos mayores a 3 veces el periodo de vibración
real de disparo con respecto al tiempo nominal. Muchos, en el pasado, han asumido
las normas indican que siempre y cuando haya una diferencia de por lo menos 8 ms
82
que los iniciadores dispararán de manera retardada uno del otro. Uno debe de ser
En general, uno puede asumir que dispararán con una distribución normal, donde la
Los iniciadores deben tener, bajo buenas condiciones, una desviación o dispersión
iniciadores nuevos al salir de la fábrica. ¿Qué sucede con iniciadores viejos que fueron
fabricados uno, dos o hasta tres años de que fueran utilizados?, se sabe que el tiempo
tiempo de retardo que son tanto inesperados cómo desapercibidos. Los traslapes,
todas direcciones. El mismo fenómeno puede ocurrir también con las ondas en el aire
si los taladros se retardan de tal manera que la onda de sonido de uno alcance a la otra
UBICACIÓN
de Lima (ver Anexo 1); a una altitud de 4,300 a 5,000 m.s.n.m. Ubicada alrededor de
las siguientes coordenadas: 10° 36´34´´ Latitud Sur y 76°59´56´´ Longitud Oeste.
La Tabla 2.1 muestra que el acceso a la mina puede ser por una ruta principal que
GEOLOGÍA
Las rocas que hospedan la mineralización son mayormente calizas y margas del
capas rojas del Santoniano. Los sedimentos cretácicos fueron fuertemente plegados y
fallados antes de la deposición de las rocas volcánicas terciarias del Grupo Calipuy
(Cobbing y Garay, 1971; Cobbing, 1973; Romaní, 1982 y Bussell et al. 1990; Petersen
et al. 2000).
Los volcánicos Calipuy son posteriores a los plegamientos y a las fallas de sobre
2.2.2.1 Pliegues
Las fases compresivas han plegado los sedimentos cretácicos formando los anticlinales
El área ha sido afectada por numerosas fallas en diversas etapas, a nivel regional se
ese sentido, aunque también los movimientos verticales son importantes. Las
Esta falla y sus estructuras asociadas están íntimamente ligadas a los procesos
buzamiento.
86
de rumbo E-W; e indistintamente en las tres zonas existen fracturas NW-SE. Todas
intenso de calcita.
mineral.
especialmente de zinc.
se extrae zinc, se observa además una amplia gama de minerales de ganga muchos de
rara naturaleza.
que tienen una reserva geológica de 258,755 TCS con una ley promedio de 19.3
2.2.3.1 Mineralogía
12’055,570 oz de plata, 7,705 TMS finas de plomo y 6,349 TMS finas de zinc, en
4’171,605 TMS con 15.87 oz/t de plata, 1.11 % de plomo, 1.50 % de zinc y 9.32 % de
88
Los recursos minerales del distrito minero, al 31 de diciembre de 2014, que incluyen
a los de Uchucchacua y Yumpag, totalizaron 3’848,555 TMS con 16.88 oz/t de plata,
oz finas de plata, 51,429 TMS finas de plomo y 68,787 TMS finas de zinc.
METODO DE EXPLOTACIÓN
Figura 2.1: Representación gráfica del método de explotación de corte y relleno ascendente
Fuente: http://intrawww.ing.puc.cl
neumáticas,
• Sostenimiento con pernos split sets, malla electro soldada, shotcrete, cuadros de
• Sistema de bombeo en los niveles 3990, 3920 y 3850 para drenar el agua a través
del túnel Patón, ubicado en el nivel 4120, con una longitud de 4,560 m,
utilizan jumbos “upper drill” en los tajeos de mayor dimensión; así como perforadoras
45 mm (Anexo 4, 5 y 6).
La voladura en los tajeos se lleva a cabo con anfo (Examón P) explosivo, dinamita
En las labores de avance se utiliza anfo, emulsiones (Emulex 80) en zonas con
coronas.
Los factores de carga, en los primeros 3 trimestres del año 2014, en los tajeos fue en
promedio 0.34 kg/t para jumbos y 0.36 kg/t para stoper. En las labores de avance, se
m/disparo y un factor de carga lineal de 37.70 kg/m para perforaciones con barrenos
de 14 pies.
de capacidades de 4.1 yd3, 3.5 yd3 y 2.2 yd3, que trabajan principalmente en la mina
Carmen y Socorro.
dumper hasta el echadero del nivel 3990, luego a través de locomotoras hasta el Pique
Luz y Pique Master, para que finalmente sean izados hasta el nivel 4450 y
2.3.2.3 Sostenimiento
Generalmente casi todas las labores utilizan pernos split set combinados con malla
electro soldada, estos van dispuestos en forma de rombo a una distancia de separación
También es frecuente en labores de roca tipo IV, principalmente en los tajeos, shotcrete
2.3.2.4 Relleno
relleno detrítico. Se recicla el 100 % del desmonte generado por las labores de
desarrollo y exploraciones.
CAPÍTULO 3: MONITOREO, MODELAMIENTO Y ANÁLISIS DE
VIBRACIONES
• Los geófonos se instalan sobre roca sólida, casi siempre pegado a los hastiales,
donde es más fácil de encontrar roca no removida. Si por las condiciones del
cemento para crear una base plana que permita que el geófono quede fijo al
macizo rocoso a través de un perno y una tuerca. Para obtener mejores resultados
Como lo describimos al inicio, las labores que se eligieron para ser monitoreas son
Labor con condiciones de agua subterránea casi nula, por lo que se pudo hacer una
para el primer caso el explosivo usado fue emulsión encartuchada tipo Emulex 80 de
1 7
1 " x 12 " en los taladros de producción, dinamita tipo Semexa 45 de "x7"y
2 8
1
Semexa 65 de 1 "x 7" en los taladros de la corona.
8
El monitoreo se realizó en el mes de mayo del año 2014, según la Tabla 3.1.
1
Se usó anfo tipo Examón P para los taladros de producción, Emulex 80 de 1 4 " x 12 "
El monitoreo se realizó en el mes de junio del 2014, tal como muestra la Tabla 3.2.
Se usó Examón P para los taladros de producción, Emulex 80 como cebo y dinamita
en la corona.
El monitoreo se realizó en el mes de julio del 2014, tal como muestra la Tabla 3.3.
El monitoreo se realizó en el mes de julio del 2014, tal como muestra la Tabla 3.4.
Tabla 3.5.
Una vez tomado la data en campo con los equipos de medición, se recolecta la
información para encontrar las leyes de atenuación que gobierna cada una de las
Con los datos que muestra la Tabla 3.1, se obtuvo la siguiente ley de atenuación:
−𝟏𝟏.𝟕𝟕𝟕𝟕
𝑫𝑫
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟔𝟔𝟔𝟔𝟔𝟔. 𝟑𝟑𝟑𝟑 � � (Ecuación 3.1)
�𝑸𝑸
97
−𝟏𝟏.𝟕𝟕𝟕𝟕
𝒍𝒍 𝑯𝑯 + 𝒙𝒙𝒔𝒔 − 𝒙𝒙𝟎𝟎 𝒙𝒙𝟎𝟎 − 𝒙𝒙𝒔𝒔
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟔𝟔𝟔𝟔𝟔𝟔. 𝟑𝟑𝟑𝟑 � �𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � � + 𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � ��� (Ecuación 3.2)
𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎
Modelo de Devine:
−𝟏𝟏.𝟔𝟔𝟔𝟔
𝑫𝑫
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟏𝟏𝟏𝟏𝟏𝟏. 𝟎𝟎𝟎𝟎 � � (Ecuación 3.3)
�𝑸𝑸
−𝟏𝟏.𝟔𝟔𝟔𝟔
𝒍𝒍 𝑯𝑯 + 𝒙𝒙𝒔𝒔 − 𝒙𝒙𝟎𝟎 𝒙𝒙𝟎𝟎 − 𝒙𝒙𝒔𝒔
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟏𝟏𝟏𝟏𝟏𝟏. 𝟎𝟎𝟎𝟎 � �𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � � + 𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � ��� (Ecuación 3.4)
𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎
98
Modelo de Devine:
−𝟏𝟏.𝟒𝟒𝟒𝟒
𝑫𝑫
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟏𝟏𝟏𝟏𝟏𝟏. 𝟑𝟑𝟑𝟑 � � (Ecuación 3.5)
�𝑸𝑸
−𝟏𝟏.𝟒𝟒𝟒𝟒
𝒍𝒍 𝑯𝑯 + 𝒙𝒙𝒔𝒔 − 𝒙𝒙𝟎𝟎 𝒙𝒙𝟎𝟎 − 𝒙𝒙𝒔𝒔
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟏𝟏𝟏𝟏𝟏𝟏. 𝟑𝟑𝟑𝟑 � �𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � � + 𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � ��� (Ecuación 3.6)
𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎
Modelo de Devine:
−𝟏𝟏.𝟖𝟖𝟖𝟖
𝑫𝑫
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟐𝟐𝟐𝟐𝟐𝟐. 𝟑𝟑𝟑𝟑 � � (Ecuación 3.7)
�𝑸𝑸
99
−𝟏𝟏.𝟖𝟖𝟖𝟖
𝒍𝒍 𝑯𝑯 + 𝒙𝒙𝒔𝒔 − 𝒙𝒙𝟎𝟎 𝒙𝒙𝟎𝟎 − 𝒙𝒙𝒔𝒔
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟐𝟐𝟐𝟐𝟐𝟐. 𝟑𝟑𝟑𝟑 � �𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � � + 𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � ��� (Ecuación 3.8)
𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎
Modelo de Devine:
−𝟏𝟏.𝟐𝟐𝟐𝟐
𝑫𝑫
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟐𝟐𝟐𝟐𝟐𝟐. 𝟓𝟓𝟓𝟓 � � (Ecuación 3.9)
�𝑸𝑸
−𝟏𝟏.𝟐𝟐𝟐𝟐
𝒍𝒍 𝑯𝑯 + 𝒙𝒙𝒔𝒔 − 𝒙𝒙𝟎𝟎 𝒙𝒙𝟎𝟎 − 𝒙𝒙𝒔𝒔
𝑽𝑽𝑽𝑽𝑽𝑽 = 𝟐𝟐𝟐𝟐𝟐𝟐. 𝟓𝟓𝟓𝟓 � �𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � � + 𝐭𝐭𝐭𝐭𝐭𝐭−𝟏𝟏 � ��� (Ecuación 3.10)
𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎 𝒓𝒓𝟎𝟎
100
NIVEL 3780
Las características geomecánicas del macizo rocoso es un factor muy importante para
paso de un tren de ondas por una roca competente será muy diferente al paso del mismo
a través de una roca fracturada, e incluso diferente cuando pase por la misma roca
las labores monitoreadas. Para la Galería 661 se subdivide en 02 zonas (tanto para la
101
Con los valores de la Tabla 3.6, y usando la Ecuación 1.32, podemos calcular la
velocidad pico de partícula crítica para cada labor. Dicho valor será la base para
Tabla 3.7 nos muestra los valores de la velocidad pico de partícula crítica para cada
de daño podemos saber el daño que generan las vibraciones a distancias cercanas a la
excavación; y que nos permitirá obtener nuevos parámetros de diseño. Antes de entrar
una tarea delicada, que exige el reconocimiento riguroso de los mecanismos que
Los criterios de prevención de las vibraciones generadas por voladuras han sido objeto
de numerosos estudios desde comienzos del siglo XX. Entre ellos destacan los de
Rockwell en 1927, Thoenen y Winders en 1942 que utilizaban como parámetro más
cimentaban las estructuras y proponían unos criterios de ámbito más general. Durante
esos niveles de prevención de los distintos tipos de construcciones, tal como hizo
además del tipo de roca donde se asienta la edificación y el tipo de estructura que se
Australia, la DIN (1983), etc. En la Figura 3.6 se hace una comparación de las
Maik (1979), Walker, Young y Davey (1981), Sisking, Stagg, Koop y Dowding
(1981), etc., dirigieron sus esfuerzos hacia la correlación de las respuestas de las
estructuras con los daños producidos por distintas intensidades de vibración a través
del análisis de los espectros sísmicos. Un hecho que se puso de manifiesto en estos
trabajos es la importancia cada día mayor que poseen las bajas frecuencias. Pero, a
organismos, haciéndose máxima cuando los estudios tienen un carácter local. Además,
hay que añadir que en muy contadas ocasiones se llegan a dar recomendaciones o bases
utilizando como base la velocidad pico de partícula crítica (𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉), que está
relacionado con las propiedades elásticas y dinámicas del macizo rocoso en estudio.
Se tomó como línea base el estudio de vibraciones realizado por la empresa Exsa en el
año 2012 (ver Anexo 3), cuya ley de atenuación se muestra en la Figura 3.8.
La Figura anterior nos muestra una ley de atenuación con valores de las contantes
𝐾𝐾 = 234.75 y 𝛼𝛼 = 1.42, para tipos de roca III. El criterio de daño en la corona y los
muestra en el Anexo 4 y 5.
El Anexo 6 nos muestra el diseño de la carga explosiva para cada uno de los taladros
del nuevo modelo de perforación y voladura, donde se observa que: Para longitudes
perforadas de 3.30 m la densidad lineal de carga (𝑙𝑙) en la corona y los hastiales, tanto
para el modelo con anfo y emulsión, es de 0.17 kg/m y 1.47 kg/m respectivamente.
𝐻𝐻
Entonces, el criterio de daño para el nuevo modelo, en puntos (𝑥𝑥𝑠𝑠 + 2
, 𝑟𝑟 ), es el
siguiente:
a.1. Corona:
a.2. Hastiales:
b.1. Corona:
b.2. Hastiales:
c.1. Corona:
c.2. Hastiales:
d.1. Corona:
d.2. Hastiales:
e.1. Corona:
e.2. Hastiales:
controlada que se está aplicando; ya que los taladros de contorno, que generalmente
detonan últimos, son los que están más próximos a la excavación. Además la
energía de vibración de los demás taladros, desde el arranque hasta las ayudas, han
debilitado y cambiado las propiedades iniciales del macizo rocoso, por lo que
cuando detonan los últimos taladros se encuentran con un macizo mucho más
(kilogramos de explosivo por longitud cargada); esta es una razón del porque en
desde el arranque hasta la corona (ver Anexos 4 y 5 y 6). Los cambios realizados son
los siguientes:
• Los 4 taladros vacíos en la corona servirán como cara libre para amortiguar el
daño de los demás taladros y para definir mejor la forma del contorno de la
excavación,
• Se tuvo que diseñar una nueva malla de perforación y voladura para el carguío
Figura 3.9: Malla de perforación y voladura inicial (derecha) y replanteada (izquierda), para anfo
como carga principal
y 4 de 102 mm que sirven como cara libre. La particularidad de esta malla es la forma
PARÁMETROS
Φ Taladro cargado 45 mm
Φ Taladro de alivio 102 mm
Longitud del barreno 3.6 m
EFICIENCIAS
DATOS UNIDADES RESULTADO
TOTAL DE EXPLOSIVO Kg 113.75
LONGITUD DE TALADRO m 3.30
AVANCE TEORICO (90%) m/disp 3.00
SECCION m2 16.00
VOLUMEN m3 48.05
TONELADAS t 129.73
Kg/m 37.88
FACTOR DE CARGA Kg/t 0.88
Kg/m3 2.37
Para una longitud de perforación efectiva de 3.30 m, se necesita una carga total
van cargados con dinamita encartuchada tipo Semexa 65 como cebo y Semexa 45 en
la columna de carga.
cada taladro, como muestra la Figura 3.11. Pero para asegurar la detonación total de la
carga, se instala un cable de cordón detonante sobre toda la longitud del taladro.
generará una presión en las paredes del taladro menor con respecto a los demás tipos
excavación final.
de Holmberg y Persson, para cada labor, simulamos el daño que se origina alrededor
La simulación y medición en campo nos determinó que se obtiene menor daño cuando:
114
1
1. La carga de fondo es de: 1 cartucho de Semexa 65 de 1 8 " x 7 ",
7
2. La columna de carga tiene: 5 cartuchos de Semexa 45 8 " x 7 ",
La simulación y medición en campo nos determinó que se obtiene menor daño cuando:
1
1. La carga de fondo es de: 2 cartucho de Semexa 65 de 1 8 " x 7 ",
7
2. La columna de carga tiene: 5 cartuchos de Semexa 45 8 " x 7 ",
roca de tipo III con anfo como explosivo principal y dinamita para la voladura
controlada, podemos hacer una comparación de los resultados de dicho estudio con
Tabla 3.9: Criterio de daño del modelo inicial versus el modelo replanteado
Línea Base Modelo Replanteado Variación (m)
Criterio de Daño
Corona Hastiales Corona Hastiales Corona Hastiales
Intenso fracturamiento 0.125 0.270 0.107 0.225 0.02 0.05
Creación de nuevas fracturas 0.215 0.455 0.082 0.182 0.13 0.27
Extensión de fracturas preexistentes 0.570 1.205 0.150 0.353 0.42 0.85
Total 0.910 1.930 0.338 0.760 0.572 1.170
La Tabla 3.9 nos muestra como hemos reducido el daño al macizo rocoso en las tres
logrado reducir el daño en 0.02 m en la corona y 0.05 m en los hastiales, para la zona
116
daño generado en la Galería 661 cuando el carguío de los taladros se realiza con
La Tabla 3.10 nos muestra dicha comparación, donde observamos que el daño
posvoladura fue 0.26 m menos en la corona y 0.49 m menos en los hastiales cuando
Analizando de la misma manera para el Crucero 6685, se puede hacer una comparación
del daño al macizo rocoso cuando el carguío es con anfo tipo Examón P y Examón Q.
La Tabla 3.11 muestra que en la corona el daño fue 0.08 m menos cuando se usó
Examón P, sin embargo, en los hastiales el daño fue 0.07 menos solamente en la zona
de intenso fracturamiento, más no en las otras dos zonas donde se obtuvo un menor
Si deberíamos elegir entre uno de estos dos tipos de explosivo, deberíamos tener en
cuenta que al momento del carguío con Examón Q, éste genera una excesiva polución
del aire debido a la gran cantidad de partículas finas que posee, y que además aumenta
el remanente de dicho explosivo después del carguío (las mediciones nos arrojaron un
SISMOGRAFÍA
denominado Proyecto Masi y que fue llevado a cabo por la Empresa Exsa durante los
El proyecto Masi tenía como objetivo principal diseñar una malla de perforación y
voladura estándar de acuerdo al tipo de roca, que permita optimizar en forma global el
avance lineal (metros por disparo), pero a la vez obtener un mínimo daño al macizo
Otro factor importante que se consideró en dicho proyecto es el tema referido al daño
principalmente por: La excesiva altura, medida desde el piso de cada labor, del
Los resultados de dicho proyecto arrojaron un ahorro, entre los meses de setiembre y
noviembre del año 2014, de 37,572 dólares americanos, de los cuales 25,110 fue por
y malla electro soldada en el 100 % de las labores que presentan roca tipo III, en rocas
voladura replanteada para el caso del carguío con anfo y con la nueva malla diseñada
para el caso del carguío con emulsión, es cuantificando el ahorro que se obtuvo al
reducir el daño al tipo de sostenimiento para roca tipo III (el daño pos voladura de la
malla electro soldada). Se consideró este método por las limitaciones de acceder a la
ascendió a 7,862 dólares americanos entre los meses de setiembre y noviembre del año
a la línea base.
en la corona debe ser de 0.17 kg/m para perforaciones con barrenos de 12 pies y
• El Anfo tipo Examón P generó menos daño en la corona que el Anfo tipo
reducir la sobre excavación o rotura excesiva en los contornos, que nos permitirá
de vibraciones, ya que de eso depende los resultados y análisis para nuestra toma
de decisiones.
• En todo experimento mientras más ensayos se tenga mejor serán los resultados,
Devine para hallar los valores de las constantes 𝐾𝐾 y 𝛼𝛼. Con el valor de dichas
• Caldas J., Sauri J., Farfán C., Enero 1998, “Levantamiento geológico de la
Región Uchucchacua”, Informe Técnico BISA.
• Contreras W., 2009, “Selección del explosivo adecuado y carga máxima por
retardo usando el monitoreo, modelamiento y análisis de vibraciones”, Tesis,
Lima, Perú, cap. 1, p. 12 - 34.
• Exsa, “Manual Práctico de Voladura”, Lima, Perú, cap. 16, p. 353 - 357.
• Fribla H., Pérez J., 2005, “Propuesta de estudio sobre el uso de las micro
fracturas en el negocio minero”, IV Jornadas de Tronadura.
• López J. C., López J. E., García B. P., 2003, “Manual de perforación y voladura
de rocas”, Madrid, cap. 25, p. 502 - 516 y cap. 33, p. 663 - 699.
• Navarro V., Marques P., 2004, “El Blastware III y Milinreg III como
herramientas para la prevención y control ambiental de vibraciones en
voladuras”, V Simposium Internacional de Tecnología de la Información
Aplicada a la Minería, Lima, Perú, p. 10 - 11.
• www.bvl.com.pe
• https://es.scribd.com/doc/285901240/Voladura-de-Rocas-EXSA
• http://intrawww.ing.puc.cl
• www.buenaventura.com.pe
ANEXOS
Geometría
Total
Longitud EFICIENCIAS
UBICACIÓN Burden Espaciamiento Taladro
de Carga DATOS
(m) (m) s UNIDADES RESULTAD
TOTAL DE EXPLOSIVO kg 114.36
Arranque 3.1 0.25 0.25 4
LONGITUD DE TALADRO m 3.30
Ayud. Arranq 1 3.0 0.50 1.00 4
Ayud. Arranq 2 3.0 0.57 1.49 4 AVANCE TEORICO (90%) m/disp 3.00
Ayud. Arranq 3 3.0 0.74 2.00 4 SECCION m2 16
Ayud. Corona 3.0 0.60 1.00 3 VOLUMEN m3 48.048
Corona ( 7 Cañas) 3.0 0.50 0.63 7 TONELADAS 129.7296
Hastial 2.5 0.80 0.80 6 kg/m 38.08
Arrastre 3.0 0.70 1.10 5
FACTOR DE CARGA kg/t 0.88
Alivios 4
kg/m3 2.38
Total Taladros Perforados 41
Total Taladros Cargados 37
Figura 4: Diseño geométrico y cantidad de explosivo para una sección de 4.0 x 4.0 m
ANEXO 4: Modelo de Perforación y Voladura Replanteado para Anfo
Geometría
Longitud
Total
UBICACIÓN de Carga Burden Espaciamiento EFICIENCIAS
Taladros
(m) (m) (m) DATOS UNIDADES RESULTADO
Arranque 3.0 0.25 0.25 3 TOTAL DE EXPLOSIVO kg 109.93
Ayud. Arranq 1 2.9 0.23 0.67 4 LONGITUD DE TALADRO m 3.30
Ayud. Arranq 2 2.7 0.33 1.00 4 AVANCE TEORICO (90%) m/disp 3.00
Ayud. Arranq 3 2.7 0.65 1.52 4 SECCION m2 16
Ayud. Arranq 4 2.6 0.79 2.20 4
VOLUMEN m3 48.048
Ayud. Corona 2.3 0.65 1.13 3
Corona ( 5 Cañas) 3.1 0.40 0.59 9 TONELADAS 129.7296
Hastial 2.5 0.80 0.90 4 kg/m 36.61
Arrastre 2.5 0.65 1.10 5 FACTOR DE CARGA kg/tn 0.85
Alivios 4 kg/m3 2.29
Total Taladros Perforados 44
Total Taladros Cargados 36
Figura 5: Diseño geométrico y cantidad de explosivo para una sección de 4.0 x 4.0 m
ANEXO 5: Modelo de Perforación y Voladura Diseñado para Emulsión
Geometría
Longitud
Total
UBICACIÓN de Carga Burden Espaciamiento
Taladros EFICIENCIAS
(m) (m) (m)
DATOS UNIDADES RESULTADO
Arranque 3.0 0.25 0.25 3 TOTAL DE EXPLOSIVO kg 113.75
Ayud. Arranq 1 2.8 0.33 0.74 4
Ayud. Arranq 2 2.2 0.37 1.00 4
LONGITUD DE TALADRO m 3.30
Ayud. Arranq 3 2.2 0.55 1.49 5 AVANCE TEORICO (90%) m/disp 3.00
Ayud. Arranq 4 2.2 0.74 2.20 4 SECCION m2 16
Ayud. Corona 2.2 0.80 1.14 3 VOLUMEN m3 48.048
Corona ( 5 Cañas) 3.1 0.50 0.59 9
Hastial 2.2 0.80 0.90 4
TONELADAS 129.7296
Arrastre 2.5 0.50 1.10 5 kg/m 37.88
Alivios 4 FACTOR DE CARGA kg/tn 0.88
Total Taladros Perforados 45 kg/m3 2.37
Total Taladros Cargados 37
Figura 9: Onda resultante de vibraciones en el Crucero 6685 por tipo de taladro y con Examón P
ANEXO 10: Simulación del Daño al Macizo Rocoso Circundante
Figura 10: Simulación del daño con el software JKsimblast 2D Face en la Galería 661
usando emulsión
Figura 12: Simulación del daño de 01taladro de la corona para el modelo con
emulsión
Figura 13: Simulación del daño de 01taladro de los hastiales para el modelo con
emulsión
ANEXO 10: Simulación del Daño al Macizo Rocoso Circundante (continuación)
Figura 14: Simulación del daño de 01taladro de la corona para el modelo con anfo
Figura 15: Simulación del daño de 01taladro de los hastiales para el modelo con
anfo
ANEXO 11: Fotografías del Monitoreo de Vibraciones
Foto 3: Voladura controlada antes del Foto 4: Voladura controlada después del
estudio, Rampa 7046 estudio, Rampa 7046
Foto 13: Tipo de emulsión utilizada Foto 14: Tipo de anfo utilizado