Science">
Tesis Huancavelica Minas
Tesis Huancavelica Minas
Tesis Huancavelica Minas
FACULTAD DE INGENIERÍA
PRESENTADO POR:
1
DEDICATORIA:
Roberto.
AGRADECIMIENTOS
la Facultad de Ingeniería
iii
3
ÍNDICE
Pág.
Dedicatoria ii
Agradecimientos iii
Índice iv
Introducción viii
Resumen x
Capítulo I: PROBLEMA
MARCO TEÓRICO
2. Antecedentes 14
1.
2. Consideraciones geomecánicas para diseñar la voladura 15
2. 2.2.1. Clasificación geomecánicas de Protodyakonov 15
2.2.2. Clasificación geomecánicas de Bieniawski 15
2.2.3. Uso de las clasificaciones geomecánicas del macizo
rocoso para la estimación del GSI
22
2. Caracterización del macizo rocos para el diseño de voladura 26
3. 2.3.1. Factor e volabilidad de Borquez (1981) 26
2.3.2. Propuesta de la compañía ”Steffen Robertson and 27
Kirsten ltd”
2.3.3. Volumetric joint count “Jv” y la ecuación de Palsmtrom 27
2.3.4. Propuesta de Ashby (1977) 28
2.3.5. Índice de volabilidad de Lilly (1986-1992) 29
2.3.6. Impedancia y sísmica de refracción 31
4
2.4. Diseño y cálculo de voladura en galerías y túneles (Metodología sueca) ……32
2.4.1. Sistemas de avance
…………………………………………………………………………………………..32
2.4.2. Esquemas de voladuras
32
2.4.3. de arranque con taladros paralelos
34
2.5. Calculo de voladuras 38
2.5.1. Avance por disparo
38
2.5.2. Arranque y corte de cuatro secciones
39
2.5.3. Diseño de arrastres
43
2.5.4. Diseño de núcleo
44
2.5.5. Diseño de contorno
45
2.6. Factor de seguridad “FS” 46
2.6.1. Determinación del factor de seguridad en voladura subterránea 46
2.6.2. Análisis de fragmentación
47
2.6.3. El modelo Kuz-Ram
47
2.7. Las condiciones geológicas y sus efectos en la voladuras 48
2.7.1. Estratificación o bandeamiento 48
2.7.2. Esquistosidad
48
2.7.3. Fractura
49
2.7.4. 49
2.7.5. Contactos
49
2.7.6. Pocas estructuras o estructuras ampliamente separadas
49
2.7.7. Estructuras apretadas
50
2.7.8. Estratificación plana u horizontal
50
2.7.9. Estructuras en trabajos subterráneos
51
2.8. Principales criterios para seleccionar explosivos
53
2.8.1. Propiedades geomecánicas
53
2.8.2. cas masivas resistentes
53
2.8.3. cas muy fisuradas
53
2.8.4. Rocas formadas por bloques
53
2.8.5. cas porosas
54
2.8.6. Volumen de roca a volar
534
v
2.8.7. Condiciones atmosféricas 55
2.8.8. Presencia de agua 55
2.8.9. 55
2.8.10 Transmisión o simpatía 56
2.8.11. Vibraciones 56
2.9. Hipótesis 56
2.9.1. Hipótesis general 56
2.9.2. ipótesis específicos 56
2.10. Definición de términos 57
2.11. Identificación de variables 58
2.11.1.Variable independiente 58
2.11.2. Variable dependiente 58
2.12. Definición operativa de variables e indicadores 58
vi
6
3.7.2. Los instrumentos 62
3.8. Procedimientos de recolección de datos 62
vi
3.9. Técnicas de procesamiento y análisis de datos 62
Conclusiones 87
Recomendaciones 88
Referencias bibliográficas (Modelo Vancouver). 89
Anexos
7
vii
INTRODUCCIÓN
como subterráneas.
remover.
8
viii
respectivos objetivos de investigación, seguido de la justificación e
importancia.
viii
bibliografía actualizada y finalmente los definición de términos básicos
utilizados en la investigación, seguido de la hipótesis y el sistema de
variables.
análisis y procesamiento de la información tanto del pre test como del post
Los autores.
9
RESUMEN
de Kaushik Dey & Phalguni Sen (India escuela de minas). Este principio nos
La mina piloto donde se recolecta datos de campo para las pruebas de esta
roca para diferentes tipos de roca y encontró una relación lineal con el
para diseño de frentes, por ser este método más aceptado en el campo de
Para concluir se hace una programación en Excel para una rápida toma de
10
CAPÍTULO I
PROBLEMA
aun hoy en día las fallas en la voladura es un desafío que todos los mineros
RMR. Pero esto solo de manera empírica, sin ninguna ecuación matemática
11
c. ¿Qué tan necesario es involucrar las variables de las
de la voladura.
importante.
2.1 ANTECEDENTES
(2).
Retamas” (3).
2.2 CONSIDERACIONES GEOMECÁNICAS PARA DISEÑAR LA
VOLADURA
Cuadro N° 01
Siendo:
(Ecu. 2.1)
macizo rocoso:
definiéndose unos valores para dichos parámetros, cuya suma, en cada caso
nos da el índice de Calidad del RMR que varía entre 0 – 100. (5)
macizo rocoso.
problemas de ingeniería.
Cuadro Nº 02
Primer procedimiento
Estimación de la Resistencia Compresiva mediante el martillo Schmidt
de Dureza.
Segundo procedimiento
Determinación de la Resistencia Compresiva mediante el Ensayo de
Tercer procedimiento
Determinación de la Resistencia Compresiva mediante el Ensayo de
Primer procedimiento
Se calcula midiendo y sumando el largo de todos los trozos de testigo
(Ecu. 2.2)
Segundo procedimiento
Comprende el cálculo del RQD en función del número de fisuras, por
metro lineal, determinadas al realiza r el levantamiento litológico-
minera.
Fórmula matemática:
(Ecu. 2.3)
Siendo:
(Ecu. 2.4)
Tercer procedimiento
Comprende el cálculo del RQD en función del número de fisuras, por metro
Siendo:
Así resulta que un material rocoso de alta resistencia de 100 a 200 MPa,
DESCRIPCIÓN
ESPACIO DE TIPO MACIZO
ESPACIAMIENTO JUNTAS
>3m ROCOSO
Muy ancho Ancho Sólido Masivo
1-3m
Moderadamente En bloques
0.3 - 1 m 50 -
Cerrado Fracturado
300 mm
Muy cerrado machacado
Fuente: “Engineering Rock Mass Classifications”, Bieniawski, Z. T.
a. Apertura
Descripción Separación
Abierta > 5 mm
Moderadamente 1 - 5 mm
abierta
Cerrada 0.1 – 1 mm
ejecutar.
Cuadro N° 05
(Ecu. 2.7)
2.2.3Uso de las clasificaciones geomecánicas del macizo rocoso para la
Según Hoek and Brown (1980), propusieron utilizar para la estimación de las
consideran para determinar la calidad del macizo rocoso. Parte del cuadro
Lunde. (10)
utilizada para estimar el valor GSI de una manera similar a lo descrito para
orientación de Juntas. Nótese que el valor mínimo que se puede obtener con
final, llamada RMR89, puede ser utilizada para estimar el valor de GSI.
Según Barton, N., Lien, R., and Lunde, J. (1988), Para estimar el valor de
(Ecu. 2.10)
Este Valor de Q´ puede ser utilizado para estimar el valor GSI a partir de:
Figura Nº 02
(Ecu. 2.11)
Figura Nº 03
Según Barton, N., Lien, R., and Lunde, J. (1988), determina el factor de
Cuadro N° 06
Resistencia de
Figura Nº 05
Factor de volabilidad (Kv)en funcion del indice de calidad RQDE
Fuente: Manual de perforacion y voladura- Lopez Jimeno
2.3.2 esta de la compañía “Steffen Robertson and kirsten ltd.”
(1985)
compañía Steffen Robertson and kirsten ltd, lanza una manera de calcular
Figura Nº 06
Consumo específico en función a diversos parámetros geomecánicos
cúbico, obtenido al sumar las juntas presentes por metro para cada una de
de la siguiente forma:
disparo.
Figura Nº 07
Fuente: Block size and shape, Arild Palsmtrom
Figura Nº08
Correlación entre la frecuencia de fracturación y el consumo específico de
explosivo
Este índice se aplicó por primera vez en las minas de hierro de Pilbara,
densidad es de 4 t/m3.
Figura Nº 09
que los rangos de factor de carga se mueve entre: 0.25 y 0.3 kg/ton de
de propagación. (18)
Figura Nº 11
Correlación entre velocidad sísmica y “CE”
Fuente: Manual de perforación y voladura- López Jimeno
Como puede observarse, conforme aumenta la velocidad sísmica se
15%.
(METODOLOGÍA SUECA)
Tipo de roca
Tipo de sostenimiento
Sistema de ventilación
donde la sección resulta demasiado grande para ser cubierta por el equipo
parte de la sección.
Figura Nº 11
Zonas de una voladura en un frente
Los taladros de contorno son los que establecen la forma final de una
elige normalmente esta última, ya que se evita la caída libre del material,
fragmentado.
2.4.3 de arranques con taladros paralelos
El mismo Carlos López Jimeno/Emilio López Jimeno, (2003), dice que las
el confinamiento de las cargas alto. Por otro lado, las dimensiones del
burden en el arranque son pequeñas, por lo que los explosivos deben ser lo
elevada, superior a los 3000 m/s, para evitar el efecto canal en los
cartuchos por delante del frente de la onda de choque destruyendo así los
Debido a que cada día las minas se mecanizan más, se emplea equipos
ángulo sería algo complejo, por tanto la perforación con taladros paralelos
Todos los taladros dentro del arranque se sitúan muy próximos alineados y
paralelos, por lo que es muy habitual usar jumbos dotados con paralelismo
automático.
tipos de arranques cilíndricos que se han utilizado con éxito y están bien
experimentados.
crítica o de muerte.
propicias a los fallos de los arranques que las de grano grueso, debido al
mayor volumen del hueco de alivio que se precisa para la salida del
material.
espacio disponible para la expansión de los taladros que aún no han salido,
cargas en espiral.
Figura Nº 13
Arranque cilíndrico en espiral
espiral. Los taladros 1-2, 3-4 y 5-6 se corresponden en cada una de sus
espirales respectivas.
Figura Nº 14
Arranque y corte de doble espiral
(57 mm), que constituyen el hueco libre en forma de “8” para las primeras
Figura Nº 15
Arranque Coromant
mientras que otros se dejan vacíos, al ser tan elevadas las concentraciones
como ocurre con los arranques cilíndricos. Los avances son reducidos y no
Figura Nº 17
Ejemplo de arranques quemados
taladros cargados. Siempre que esta última se mantenga por debajo del 2%
(Ecu. 2.18)
Figura Nº 18
Arranque de cuatro secciones
del arranque, por eso se recomienda que el burden se calcule sobre la base
de:
(Ecu. 2.19)
Figura Nº 19
Resultados para diferentes distancias de los taladros cargados a los vacíos
y diámetros de estos.
Dónde:
burden igual a vez y media del diámetro del taladro vacío la concentración
Dónde:
= Concentración lineal de carga (kg/m)
(Ecu. 2.21)
= Diámetro de perforación (m) = Diámetro del taladro vacío (m) B
= Burden
C = Constante de roca
RWSANFO = Potencia relativa en peso del explosivo referida al ANFO.
carga están bastante limitados, por cuanto no existe una variedad amplia
complejo.
Figura Nº 20
Concentración lineal de carga en función al burden máximo para diferentes
diámetros de broca
por lo que:
Figura Nº 21
(Ecu. 2.23)
(Ecu. 2.24)
Este valor tiene que reducirse con la desviación de los taladros para
B2 = B – EP (Ecu.2.25)
B 2 ≤ 2 Ah
Para que no se produzca solo la deformación plástica. Si esto no se cumple,
(Ecu. 2.26)
la fragmentación.
El ángulo de apertura debe ser también menor de 1.6 radianes (90°) ,pues
significa que:
B2 > 0.5 Ah
Gustafsson (1973), sugiere que el burden para cada sección se calcule con.
longitud del lado de la última sección “B” no sea menor que la raíz cuadrada
del avance. El método de cálculo del resto de las secciones es el mismo que
(Ecu. 2.29)
Dónde:
f = Factor de fijación, generalmente se toma 1.45 para tener en cuenta el
Figura Nº 22
Geometría de los taladros de arrastres
FACTO
DIRECCIÓN DE SALIDA RELACI
R DE ÓN
DE LOS TALADROS
FIJACI S/B
Hacia arriba y 1.45 1.25
horizontalmente
Hacia abajo 1.20 1.25
Fuente: Manual de perforación y voladura- López Jimeno
se calculan de acuerdo con lo indicado para los taladros de arrastres con los
siguientes valores:
= 1.25
4.5 x 4.5 m
Si:
(Ecu. 2.31)
Despejando “FS”
(Ecu. 2.32)
arranque.
Figura Nº 24
determinada malla.
(Ecu. 2.33)
(Ecu. 2.34)
Dónde:
fracturadas.
Si se sabe que:
(Ecu. 2.35)
(Ecu. 2.36)
También:
(Ecu. 2.37)
Cunningham (1983) indica que en su experiencia el límite más bajo para A
Según Antonio Karzulovic, (2006), los planos que dividen a las capas o
2.7.2 Esquistocidad
Se rompen fácilmente.
2.7.3 Fractura (joints, fisuras o juntas)
(disyunción) y otras.
Las rocas son propicias a sobrerotura (over break, back break) junto a los
planos de falla.
2.7.5 Contactos
Pueden ser una desventaja para la fragmentación por los siguientes motivos:
estratos (incompetencia).
con mejor fragmentación y control del disparo. Las rocas con baja
controlarse la vibración.
Estructuras predecibles
La perforación perpendicular a estratos horizontales.
reduce la probabilidad de que se traben o agarren los barrenos.
Los taladros son verticales y rectos ya que estos planos no afectan por
desviación.
En estas condiciones son factibles de aplicar opciones técnicas en mallas,
taladros.
2.7.9 Estructuras en trabajos subterráneos
subsuelo. Caso especial son los túneles, galerías, rampas y piques donde los
Los sistemas dominantes clasificados con relación al eje del túnel son tres:
Figura Nº 26
Fuente: Manual práctico de voladura -EXSA
2.7.9.3 Sistema de fracturas o juntas en ángulos variables con
relación al eje del túnel
En estos casos usualmente los taladros de un flanco trabajan mejor que los
del otro. Puede decirse que los del lado favorable trabajan “a favor del
impulsores.
Figura Nº 27
Los explosivos con una alta “ET” tienen en esos macizos muy poca
que posean una elevada energía de los gases “EG”, como es el caso del
ANFO.
2.8.4 Rocas conformadas por bloques
siguientes medidas para retener los gases dentro de los taladros el mayor
tiempo posible.
Figura Nº 28
-20°C.
líquido.
labores con pendientes negativos como rampas, para lo cual habrá que
que son muy recomendables para asegurarse que el taladro este sin agua.
(Emulex, Emulnor u otros), pero esto siempre con criterio es decir según
2.8.9 Humos
que con los hidrogeles a granel hay que tomar ciertas precauciones, lo
mismo que con el ANFO que produce una elevada concentración de gases
2.8.11 Vibraciones
2.9 HIPÓTESIS
independientes y cartografiables.
hidráulico.
macizo rocoso, con una resistencia a la tracción nula o muy baja. (Genera
anisótropo).
hace exploraciones,
luego perforaciones diamantinas y, si se encuentra mineral, se procede a
hacer una mina. Se construye una galería principal con una entrada que se
llama bocamina.
2.11IDENTIFICACIÓN DE VARIABLES
Cuadro N° 09
VARIABLE TIPO DE
DIMENSI INDICADORES
INDEPENDIE VARIABLE
ÓN
- Conocimiento de la
aplicación de las
clasificaciones
geomecánica en los
Evaluación de X1
Cognitivo distintos tipos de
la voladura
macizos rocoso en la
basada en la Variable
X2 minería subterránea.
clasificación Independie
de voladura, las cuales
incidencia en la
productividad
Fuente: Elaboración propia y
Cuadro N° 10
as de Índices sobre
campañas La tasa actitudes ambientales
de crecimiento de
poblacional
Necesi
de El totalen
influye deelladesarrollo
población carece de de
sostenible
dades acceso real a las necesidades básicas.
Social necesida En la actualidad, el mineral tiene un
básicas Valor valor económico
Crecimi Los costos de la significativo.
voladura están
basados en el tipo de roca.
Econó ento económic
METODOLOGÍA DE LA INVESTIGACIÓN
Región: La Libertad.
científica”.
sucederá?”. (32)
Asimismo se hará uso del método estadístico. Según Ary, Donald y otros
3.6.POBLACIÓN Y MUESTRA
3.6.1 La Población
3.6.2 Muestra
El mismo Oseda, Dulio (2008:122) menciona que “la muestra es una parte
población”. (35)
3.7.1 Técnicas
investigador”. (36)
siguientes cálculos:
Softwar Aplicativo:
- Procesadores de Texto.
- Hojas de Cálculo.
- Bases de Datos.
RESULTADOS
cada tramo de avance por voladura (mapeo de líneas por detalle) en una
Dónde:
(0.00088IR 1.01)
c (Ecu. 4.1)
10
IR: Índice de rebote del martillo Schmidt
: Peso Específico de la roca
caso)
(0.0008824.6326.461.01)
10 c
Valoración: 4
4.1.1.2 RQD (Rock Quality Designation)
Dónde:
0.1
RQD 100 e 0.11
(Ecu. 4. 2)
: Nº de discontinuidades por metro lineal
0.115.05
RQD 100 e 0.115.05 1
RQD = 55.6 %
Valoración: 13
Figura Nº 29
HISTOGRAMA DE ESPACIADO
mayor a 2000
100
600 a 2000
200 a 600
60 a 200
80
60
40 1
ESPACIADO
20
0
FRECUENCIA
Fuente: Elaboración propia
Cuadro Nº 09
Espaciado Espacia Nº de
(mm) do datos
mayor a 2000 1 34
600 a 2000 1 55
200 a 600 1 96
60 a 200 1 49
TOTAL 200
Fuente: Elaboración propia
4.1.1.5Aguas Subterráneas
Valoración: 20
La zona presenta goteo incesante en la mayoría de los tramos analizados
Valoración: 7
RMR: 54
(Ecu. 4.3)
Dónde: RQD: 55.6%
Jn: 9
Jr: 1.5
Ja: 2
Reemplazando en la ecuación Nº 4.3
Q = 4.63
Este Valor de Q´ puede ser utilizado para estimar el valor GSI a partir de:
Entonces:
GSI = RMR89 - 5 (Ecu. 4.5)
GSI = 49
4.2 ADAPTACIÓN DE LOS ÍNDICES GSI, RMR Y RQD EN LA
DE EXPLOSIVO
(Ecu. 4.6)
Dónde:
Paso 1.
(Ecu. 4.8)
(Ecu. 4.11)
(Ecu. 4.12)
GSI, será
de la siguiente manera.
Dónde:
FACTORES DE VOLADURA
Calculamos las siguientes constantes a partir del C.E, los cuales son
C = 0.31 kg/m3
A = 9.3
Cunningham (1983) indica que en su experiencia el límite más bajo para “A”
incluso en tipos de roca muy débiles es A=8 y el límite superior es A = 12.
geomecánicamente incompetente.
FE = 10.5 MJ/ton
Variables de diseño
RMR = 55
GSI = 49.99
RQD = 55.6 %
Resistencia a la compresión = 38.33 MPa
Diámetro de broca ( ) = 45 mm
% de acoplamiento = 80 %
Eficiencia de perforación = 95%
Longitud del barreno (LB) = 12 pies
Ancho de labor = 3.5 m
Alto de labor =3.5 m
Distancia a una zona critica = 100 m
Diámetro de broca rimadora( )= 101.6 mm
Error de enboquille “ ” = 20 mm
Especificaciones técnicas de los explosivos a usar:
Cuadro Nº 10
Peso/ unid
Explosivo Dimensiones (diam x Densidad
Emulex long)
1 1/4” x 8” = 31.80 x 1.03 3)
(gr/cm (kg)
0.195
45% 200.3 mm
Exadit 45% 7/8” x 7” = 22.23 x 1 0.076
177.8 mm
Exablock 7/8” x 7” = 22.23 x 0.78 0.072
45 % 177.8 mm
Fuente: Manual práctico de voladura-EXSA
Cuadro Nº 11
Gelatina
Descripci Emulex Exadit Semexa Exsabl
ón especial ock
80 65 45 65 45 90 75 80 65 45 45%
Unid 1.1
Densidad g/c % %
1.1 %
1.0 %
1.0 %
1.0 %
1.4 %
1.3 %
1.1 %
1.1 %
1.0 0.78
VDO m3
m/s 4
50 2
50 3
45 5
36 0
34 2
60 8
55 8
45 2
42 8
38 2000
Presión 00 00 00 00 00 00 00 00 00 00
Kba 87 85 63 53 44 136 125 86 70 60 30
de r
Kcal/
Energía 12 110 70 85 80 110 105 10 95 90 600
kg 00 0 0 0 0 0 0 00 0 0
Volume
lit/k 83 91 93 94 94 86 87 91 93 93
n normal
RWS ANF g % 0 0 0 0 5 7 8 6 2 9
13 121 77 93 88 121 115 110 10 99
O
RBSANF 2
18 16 95 121 10 71 196 16 4
14 13
O
5 7 8 0 4 2
Resistencia
hora 72 72 72 2 2 11 11 12 6 4
al agua
Fuente: Manual práctico de voladura-EXSA
4.4.2 Diseño del arranque
la voladura:
(Ecu. 4.18)
Dónde:
L = (Long. Barra) x (eficiencia perforación) = (12 pies x 0.95) = 11.4pies = 3.
juntos
(Ecu. 4.19)
Reemplazando la ecuación con dos taladros juntos perforados tenemos:
0.1016 = 0.1437 m
Ahora tenemos:
L = 0.15 + 34.1 (0.1437) – 39.4 (0.1437)2 = 4.236m
Con dos taladros de expansión es más que suficiente para alcanzar hasta
Figura N° 30
Ilustración del taladro vacío equivalente
(Ecu. 4.21)
Dónde:
Ashby:
C.E=0.35 Kg/m3
Se sabe que:
= 9.32
Cunningham (1983) indica que en su experiencia el límite más bajo para “A”
geomecánicamente incompetente.
(Ecu. 4.24)
Dónde:
(m) = 0.226
Entonces:
Longitud de taco = 10 (Ecu. 4.27) Longitud de taco = 0.45 m
N cart = 12 cartuchos/tal
(Ecu. 4.28)
Paso 5. Predicción de la fragmentación Aquí entra a tallar el modelo Kuz-
Ram.
Dónde:
= tamaño medio de los fragmentos, cm. = 8” = 20 cm A = factor de roca =
9.32
resultado de manera técnica se podría decir que el 50% del material roto
(Ecu. 4.30)
Cuando existe un error de perforación la superficie libre “Ah” difiere de la
= (Ecu. 4.31)
Ah = 0.24m
B2= 0.28 m
El nuevo: = (
Ah = 0.49m
Figura N° 32
Tercer cuadrante:
= 0.4 m
Ah = 0.83 m
Figura N° 33
= 0.52 m
Ah = 1.25 m
Figura N° 34
EXPERIMENTAL)
Experimental se tiene:
Cuadro N° 01
FALLAS DE VOLADURA
2010 2011
ENERO 15
FREBERO 18
MARZO 22
ABRIL 19
MAYO 21
JUNIO 17
JULIO 15
AGOSTO 18
SETIEMBRE 17 11
OCTUBRE 19 10
NOVIEMBRE 22 11
DICIEMBRE 20 8
TOTAL 223 40
Gráfico N° 01
del grupo control (2010) los cuales se puede apreciar gráficamente. Cabe
peculiaridades, motivo por el cual era necesario decidir por uno de ellos
iguales.
Hipótesis Alterna:
diferentes.
α=0,05.
gl = ge+gc-2.
gl = 12+4-2=14.
posibilidades ( H1 : 1 2
ó 1 2 ) se aplicó la prueba bilateral, o sea a dos colas.
X 1 X 2
tc (Ecu. 4.33)
S2 S2
1 2
n1n2
tc
8.82
Figura N° 36
la empresa.
óptimo en la voladura.
A.BIBLIOGRAFÍA:
(1) Kaushik Dey & Phalguni Sen (India escuela de minas) los principios de
(12) Barton, N., Lien, R., and Lunde, J., “Engineering Classification
pp. 98
(13) Barton, N., Lien, R., and Lunde, J., “Engineering Classification
pp. 103
(15) Barton, N., Lien, R., and Lunde, J., “Engineering Classification
pp. 195
(16) Carlos López Jimeno/ Emilio López Jimeno. “Manual de
2006. pp. 67. (26)EXSA (2009) “Manual práctico de voladura” 3ra Edición.
pp. 236
B. PÁGINAS WEB:
del 2011.
94
TABLA DE CUANTILES DE LA DISTRIBUCIÓN DE LA T DE STUDENT