Traducido del inglés al español - www.onlinedoctranslator.

com

Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Listas de contenido disponibles en ScienceDirect

Revista Internacional de Mecánica de Rocas y Ciencias de la Minería

revista Página de inicio: http://www.elsevier.com/locate/ijrmms

Vibración de campo cercano por voladuras y predicción de daños en rocas con una solución
de campo completo

Santiago Gómez *, Jos´é A. Sanchidrián, Pablo Segarra

Universidad Polit´écnica de Madrid - ETSI Minas y Energía, Ríos Rosas 21, 28003 Madrid, España

INFORMACIÓN DEL ARTÍCULO ABSTRACTO

Palabras clave: Se ha desarrollado y aplicado una solución semi-analítica de campo completo para la radiación de ondas sísmicas para calcular la
Transformada de Fourier
vibración de campo cercano por voladuras. La forma de la función de fuente aplicada a la pared del pozo de voladura se ha obtenido a
Función fuente partir de un modelo numérico de detonación, con un valor máximo determinado como la presión de coincidencia de choque de los
Modelo de radiación sísmica
productos de detonación y la roca. Las señales de velocidad y aceleración de las partículas de cuatro explosiones se han registrado y
factor Q
utilizado con el fin de determinar los factores de calidad (Q) del macizo rocoso, expresado en función de la distancia radial al barreno. La
Voladura de contorno
Holmberg-Persson consistencia del modelo ha sido validada verificando que la velocidad máxima de partícula resultante en la pared del pozo de voladura
coincide con la velocidad de partícula obtenida del cálculo de coincidencia de choque entre explosivo y roca. Las zonas de fractura se
identifican utilizando un criterio de falla basado en las tensiones principales del tensor de tensiones de Cauchy resultante y las
resistencias uniaxiales a la compresión y a la tracción. Se han comparado dos variantes del enfoque estándar de la industria de Holmberg-
Persson con la solución de campo completo.

1. Introducción propagación de ondas en un medio homogéneo e isotrópico resultante de una

carga variable que actúa sobre la pared de un elemento corto de explosivo,
El daño de la roca provocado por la explosión es un problema importante en la documentado por Tubman et al. y Meredith et al.,12-15 y los desarrollos realizados
minería subterránea y de superficie, y en la construcción de túneles. Se debe predecir el por Blair,16-18 se toman como base para lo que Blair llamó solución de campo
daño de la roca para optimizar la perforación y voladura, de modo que la completo (FFS).dieciséis El término FFS significa una solución exacta tanto en el
sobreexcavación, la caída de rocas y las fallas de taludes sean limitadas, mejorando así la campo cercano como en el campo lejano. Heelan dio una solución tradicional
operación desde los ángulos económicos y de seguridad. Existe una variedad de para el campo lejano.19,20 en su obra fundamental, seguido por una larga serie de
métodos disponibles para predecir la extensión de las zonas dañadas debido a las autores, por ejemplo,7,12,16,21-25 pero también criticado por algunos,26,27
voladuras; se puede ver una revisión completa en Ainalis et al.1 Por lo general, utilizan aunque se encontró que produce resultados cercanos a los numéricos y al FFS
modelos de vibración basados en leyes simples de escalado de distancia y peso de 7,12,22 con un costo computacional relativamente bajo. Sin embargo, algunos
carga para calcular las velocidades máximas de las partículas (ppv) que, en combinación supuestos limitan su validez cuando la distancia radial desde el orificio es inferior
con fórmulas simples de esfuerzos planos, permiten calcular el esfuerzo máximo en la a 50 radios, independientemente de la frecuencia.dieciséis
roca a una distancia determinada.2 El conocido enfoque de Holmberg y Persson (HP)3 y Los profesionales rara vez utilizan la solución de campo completo (FFS). Sin
trabajos relacionados del Instituto Nacional de Seguridad y Salud Ocupacional (NIOSH)4,5 embargo, su potencial es único para predecir vibraciones, especialmente en el
proporcionan soluciones de ingeniería sencillas que se aplican con mayor frecuencia a campo cercano, en lo que se refiere al diseño de voladuras de contorno. Incluso si
barrenos para voladuras de contornos y amortiguadores. A pesar de sus defectos físicos es relativamente complejo en comparación con los enfoques HP o de carga
y limitaciones (Blair y Minchinton6,7 criticaron temprano el uso de modelos de predicción puntual, requiere mucho menos esfuerzo computacional que los modelos
de vibraciones basados en las leyes de escalamiento del peso de carga en los años 90 numéricos basados en elementos finitos o métodos de diferencias finitas, y
debido a algunas suposiciones inválidas), han sido los preferidos por la mayoría de los puede ser manejado eficientemente por computadoras personales actuales.
ingenieros de campo en los últimos cuarenta años, e incluso hoy en día se utilizan en El FFS considera el caso de una cavidad seca equivalente presurizada,28
publicaciones de primera línea (ver p. ej.8-11). Estos modelos se describen con cierto con una carga aplicada en la pared para modelar el barreno. Esta carga,
detalle y se comentan en la sección5. originada por la detonación del explosivo, avanza a lo largo del barreno con
En el presente trabajo, la teoría analítica básica para el eje simétrico una velocidad de detonación finita. La energía y la información

* Autor correspondiente.
Dirección de correo electrónico: santiago.gomez@upm.es (S. Gómez).

https://doi.org/10.1016/j.ijrmms.2020.104357
Recibido el 30 de diciembre de 2019; Recibido en forma revisada el 2 de abril de 2020; Aceptado el 25 de abril de 2020
Disponible online el 19 de agosto de 2020
1365-1609 / © 2020 Los Autores. Publicado por Elsevier Ltd. Este es un artículo de acceso abierto bajo la licencia CC BY-NC-ND
(http://creativecommons.org/licenses/by-nc-nd/4.0/).
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

contenido en este intenso pulso de presión se transmite inmediatamente a la roca informó. La metodología seguida para obtener laQ-se describen los valores; Se
y se convierte en ondas de tensión elásticas.29 Generalmente, se aplica una carga comparan las velocidades y aceleraciones medidas y calculadas por FFS. El
radial (para la onda primaria y onda cortante polarizada verticalmente, P-SV, modelo se valida a lo largo de la pared de un barreno genérico, verificando que
problema) o torsional (para la onda cortante polarizada horizontalmente, para diferentes velocidades de detonación elppv es similar a los valores obtenidos
problema SH,), dependiendo del problema a resolver. Sin embargo, en la mayoría por choque de coincidencia. El FFS se aplica a un caso de muestra para evaluar los
de los casos, los problemas de P-SV y SH no se tratan juntos porque existe una patrones de radiación para los problemas P-SV, SH y acoplados (P-SV y SH), a
dificultad fundamental asociada con el componente transversal, ya que la simetría partir del cual se calcula el campo de esfuerzos en la roca y la extensión de las
exige que la radiación de un barreno en un medio viscoelástico uniforme tenga zonas de fractura debido a la incidencia. de ondas corporales se evalúa.
una transversal cero (es decir, cero SH ) componente.13 Esto no es consistente con Finalmente, los resultados se discuten y comparan con prácticas de ingeniería
las observaciones, que siempre muestran un componente transversal distinto de comunes.
cero que es esencialmente aleatorio,18,30
debido a heterogeneidades en el macizo rocoso18; de ahí que se hayan utilizado 2. Formulación del modelo
conjuntamente fuentes sísmicas radiales y de torsión, mediante simple superposición,12
para incluir este componente transversal y analizar su influencia. 2.1. Elemento único de explosivo
El FFS necesita un historial de presión-tiempo aplicado en la roca. Las funciones
analíticas están disponibles para este propósito7,16,17,28,31-42; sin embargo, aunque estas Para evaluar el campo sísmico irradiado por un solo elemento, un sistema de
funciones son una forma sencilla de representar la acción explosiva, a menudo no coordenadas cilíndrico (r, θ, z), con el origen de las coordenadas en el centro de la
modelan adecuadamente el comportamiento del explosivo debido a tiempos de subida carga explosiva y el z Se considera el eje colineal con el eje del barreno. Se supone
prolongados o patrones de caída de presión poco realistas.17 Si se dispone de una que el barreno está rodeado por un espacio infinito de roca homogénea e
caracterización de explosivos adecuada, entonces los resultados de la detonación en el isótropa. En estas condiciones, el campo de desplazamiento se puede representar
pozo de voladura, tanto experimentales como numéricos, pueden utilizarse mediante tres potenciales de desplazamiento φ, ψ y χ, el primero
directamente con ventaja. Este es el enfoque seguido en este trabajo. correspondiente a una onda longitudinal o primaria (onda P) y los otros dos a
Un punto oscuro que a menudo se observa en relación con la presión aplicada en la ondas transversales o transversales (ondas S). Estos potenciales se dan en
interfaz explosivo / roca (es decir, en la pared del pozo de voladura) es que se considera Apéndice A. Una de las ondas transversales está polarizada verticalmente (SV) y la
que alcanza su punto máximo en la presión de von Neumann (PAGVN). De hecho, Blair otra horizontalmente (SH).12 Si bien el P-SV se resuelve en la mayoría de los casos,
16-18,24,36 lo usa en sus funciones fuente. Sin embargo, la presión de von Neumann se por ejemplo,16-18,24 Rara vez se considera el problema de SH, ya que el explosivo
alcanza cuando una onda de choque comprime el explosivo a lo largo de su hugoniot sin no actúa como fuente de torsión. Para el problema P-SV, la axisimetría implica
reaccionar,43 es decir, es un estado del explosivo sin reaccionar y no del espacio de que no hay variación permitida en las propiedades con respecto aθ, por lo tanto el
choque de rocas. En este trabajo, el cálculo de la presión máxima en la pared del pozo de desplazamiento tuθ y los derivados con respecto a θ en las relaciones de potencial
voladura se basa en un esquema de combinación de choque de explosivo a roca que de desplazamiento se desvanecen, es decir, el vector de desplazamiento ūP-SV = (tu
asegura una presión consistente con las propiedades de choque de la roca. r, 0, tuz). El movimiento SH, que es mucho más simple de calcular que el P-SV, se
puede resolver aplicando una fuente de torsión a la pared del barreno que solo
Fuera del barreno, el FFS considera que la roca es puramente elástica.17 La produce radiación de ondas SH, ya que los potenciales de desplazamiento φ y ψ
atenuación se puede introducir por medio de modelos constitutivos simples, que son idénticamente cero.12 es decir, ūSH= (0, tuθ, 0). Apéndice A muestra las
exhiben la capacidad de almacenar y disipar la energía de deformación,44 siendo ecuaciones básicas para calcular los componentes del campo de desplazamiento,
las velocidades de propagación funciones complejas de la frecuencia. Los que dan como resultado:
modelos de atenuación utilizados por el FFS se basan generalmente en la
viscoelasticidad lineal debido a su simplicidad de implementación en el dominio
a2 ∫∫ ∞ S(ω)
de la frecuencia.45 La atenuación intrínseca de la roca se incluye en los cálculos tuζ = PAG(kz )GRAMOζ(kz, ω)miikz zmi- esdkzdω (1)
2π2 - ∞ μ
mediante un factor de calidad (Q) que aumenta con la distancia a la fuente
explosiva y con valores bajos en regiones cercanas al barreno para imitar la fuerte
dónde:
atenuación en esa zona.46 Esto se hace para obtener una variación suave y natural
en la respuesta de la roca a la onda sísmica, cuya amplitud puede aproximarse
ζ denota la dirección correspondiente, es decir, radial (r), azimutal (θ) o
razonablemente por algún factor de desintegración exponencial.17 El constante-Q
axial (z).
(CQ) modelo de atenuación en frecuencia47 se utiliza para simplificar, ya que las
a es el radio del barreno (m);
propiedades de la roca dependientes de la frecuencia son a menudo
S(ω) es la respuesta espectral de la función fuente (Pa⋅s) (ver sección
desconocidas; a pesar de ser una simplificación
2.3);
ción, el modelo CQ ha demostrado ajustarse muy bien al comportamiento de la mayoría de las
ω es la frecuencia angular (rad / s);
rocas.48-50
μ es el módulo de corte (Pa);
Este trabajo proporciona un modelo FFS que se utiliza para estimar el estado
PAG(kz) es el término sinusoidal que es igual a (L/a) ⋅sincL⋅kz), ser
de estrés y a partir de él la extensión de las zonas de daño debido a las ondas
L la mitad de la longitud de un elemento de explosivo (m), y
corporales. Se presentan las formulaciones matemáticas para un solo elemento
kz el número de onda angular axial (rad / m);
de carga y para una columna extendida de explosivo. Se utiliza una formulación
dimensional en lugar de la común independiente de escala (no dimensional) (ver
GRAMOζ(kz, ω) es el término direccional en el ζ-componente, cuya expresión-
p. Ej.17). Esto incluye la función de fuente utilizada y sus características. El modelo
siones son:
se calibra a partir de datos de campo; El diseño de las pruebas, la instrumentación
y las características del explosivo y la roca son

(( ) )
la 2k2 z - k2 β K1(mamá) K1(lr) - 2k2 z K1(la) K1(señor)
GRAMOr(kz, ω) = ( ) (2)
a2 2k2z - k 2β K0(la) K1(mamá) - 2k2 βlaK1(la) K1(mamá) - 4k2
2
z mamáK0(mamá)laK1(la)

2
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

( ( ) )
I akz)-2 a2 2k2 - z k2βK1(mamá) K0(lr) + 2laK1(la)mamáK0(señor)
GRAMOz(kz, ω) = ( (3)
a2 2k2z- k2 Kβ (la)0K ( 1 metro a) - 2k laβ2K (la)1K (mamá
1 ) - 4k2mamá
z K (mamá
0 )laK (1la)

(D.1) a la ecuación. (D.4) en Apéndice D.

K1(señor)
GRAMOθ(kz, ω) = (4) 2.2. Columna extendida de explosivo
mamáK0(mamá) + 2K1(mamá)

El modelo considera un barreno recto por simplicidad y una columna

Ies la unidad imaginaria; explosiva formada por la superposición de J elementos individuales de longitud 2L
t es el tiempo (s); Kι (ξ) son las funciones de Bessel modificadas del cada uno de manera que su contribución se suma en un cierto punto y se retrasa
segundo tipo con argumento adimensional ξ y el orden ι; de acuerdo con la velocidad de detonación (D). La longitud total de la carga está
definida porCL = 2JL. LlamemosFζ (r, z, ω) la transformada de Fourier del
l y metro son los números de onda radial compresional y radial cortante componente ζ de cualquier campo F resultante de un solo elemento que se
2 β) 1 /
(metro-1) respectivamente. Están dados por:l = (k2z-k2α)1 / y metro = (k2z-k2 supone que se inicia instantáneamente.17 El campo total de la detonación de toda
2,dónde kα y kβ son los números de onda de compresión y de corte la columna en una ubicación determinada se puede evaluar utilizando los
(rad / m), kα = ω/α y kβ = ω/β, α y β siendo las velocidades de onda teoremas de cambio y suma.51
compresional y de corte complejas (m / s) respectivamente. por:
∫ {∞ }
Las velocidades complejas de las ondas P y S pueden evaluarse mediante47: ∑J ( ) a2 ∑J ( )
Fζ r, zj, t - Δtj = Fζ r, zj, ω miyo tj mi - Iω
tDω (6)
⃒⃒ ⃒⃒ 2π2
⃒ ⃒ - iπγsgn (ω) ⃒ ω ⃒η iπηsn (ω) j= 1 - ∞ j= 1
α(ω) = V ⃒ ω⃒ ⃒γ mi β(ω) = V ⃒ s⃒ ⃒Ω ⃒ mi - (5)
gramo

2 2
pag ωa ⃒ a El término entre llaves de la Ec. (6) también se puede expresar como:

dónde: ∑J ∫ () ∞ ∑J
Fζ r, zj, ω mies j = Fζ(r, kz, ω) miikz zj miesj dkz (7)
j= 1 - ∞ j= 1
Vpag y Vs son las velocidades de las ondas P y S (m / s) a la frecuencia de
referencia ωa. Haciendo zj = z1+ Δz⋅ (j-1) y tj = Δt⋅ (j-1), z1 siendo la coordenada axial
ωa es la frecuencia promedio de un historial de presión-tiempo particular de P como si se viera desde el centro del primer elemento que detona,
(rad / s). y llamando a Δz = -2L y Δt = 2 L / D, luego Eq. (7) se convierte en:
γy η son una fracción del ángulo de fase entre la tensión y la deformación47 y ( )
⎛ ⎞
se puede calcular mediante γ = (1 / π)⋅broncearse-1(1 /Qpag) y η = (1 / π)⋅ ∑J ( ∫∞ I2LJ - k +ω
) 1 - mi
zD

broncearse-1(1 /Qs), dónde Qpag y Qs son los factores de calidad de la roca para la Fζ r, zj, ω mies j = Fζ(r, kz, ω)miikz z1⎝ ( ) ⎠dkz (8)
j= 1 - ∞ I2L - k +zωD
onda P y la onda S, respectivamente. 1 - mi

Eq. (8) implica una reducción significativa en el tiempo de simulación con

El uso de las Ecs. (1) - (5) permite derivar las relaciones de desplazamiento con
respecto a otras soluciones (ver p. ej.17), ya que ya no depende del número
respecto al tiempo en un punto P de coordenadas (r, z) y obtener los campos de
de elementos de carga. Tenga en cuenta que este formulario solo se aplica a
velocidad y aceleración sin utilizar elementos finitos o códigos de diferencias
D.
finitas. Las velocidades y aceleraciones se pueden calcular multiplicando el
Los vectores unitarios de cada componente ¯mir, ¯miθ y ¯miz se muestran en Figura 1 en el
integrando interno en la ecuación.(1) por -yo y por -ω2, respectivamente. Los
punto de seguimiento P para el caso típico en la explotación de canteras de barrenos inclinados.
campos de deformación y tensión se obtienen de la Ec.(1) utilizando la ecuación.
El pozo de voladura se traza en color rojo y B y C son los valores absolutos

Figura 1. Croquis que ilustra los cálculos de un agujero inclinado. δ es el ángulo formado entre OC y OP.

3
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

coordenadas de la parte inferior y el collar, respectivamente, y O es el obtenido no varía en más del 0,5%. La integral externa se resuelve mediante el
centro de gravedad de la carga. El plano Π contiene el barreno y P cuya algoritmo de transformada rápida de Fourier inversa (IFFT) haciendo uso del
posición en el sistema de barreno (con origen O) es (r, z). hecho de que el integrando tiene simetría hermitiana enω.17
En el caso de que se coloque un sensor uniaxial (geófono o acelerómetro) en
P en la parte inferior de un orificio vertical (línea marrón en Figura 1), es necesario
proyectar los componentes (r, θ, z) en el eje vertical (Z, ver 2.3. Función fuente
Figura 1) para calcular el campo sísmico en esa dirección. Para ello, las
señales obtenidas en elζ-direcciones con Eq. (6) a (8) se multiplican por El explosivo utilizado fue una emulsión pura. Se realizó una caracterización
los directores de coseno (cosγζ) para el eje vertical con respecto a cada detallada en el contexto de un proyecto más amplio; Se realizaron mediciones de
componente ζ. Los angulosγζ se amplían en la parte derecha de Figura velocidad de detonación y pruebas de cilindros a partir de las cuales se obtuvo la
1. Sumando las contribuciones en P de cadaζ-componente e ecuación de estado de los productos de detonación. Estos datos se aplicaron para
inyectando Ecs. (7) y (8) en Eq. (6), lleva a Fv, el componente vertical del modelar numéricamente la detonación a partir de la cual se calculó la presión
campo F: aplicada a la pared del barreno.54,55 La función de presión normalizada se da en
Figura 2. El tiempo de subida es de 4,3 μs y el

( )
⎧ ⎡( ) ⎛ ⎞⎤ ⎫
a2 ∫ ∞⎨∫ ∞ ∑
I2LJ - kz+ω D ⎬
⎣
norte

ikz z1 ⎝1 - mi
Fv = Fζ(r, kz, ω) mi ( ) ⎠⎦dk -
(9)
⎭mi
porque γζ z IωtDω
2π2 - ∞⎩ - ∞ ζ= 1 I2L - k +zωD
1 - mi

la duración aproximada es de 0,6 ms. Se ha agregado una sección de caída

directa a la función original más allá de 0.4985 ms para garantizar que la presión
En el caso del problema P-SV, solo las direcciones r y z son considerados ( se vuelva efectivamente cero.
norte = 2); para el problema acoplado,norte = 3. Eq. (9) puede extenderse por La función de fuente influye significativamente en la forma de la
varios pozos de voladura. onda completa. Para calcular la respuesta espectral S (ω) de la fuente
Aunque el FFS se considera un modelo analítico, se deben aplicar varias función, se aplica la transformada rápida de Fourier (FFT) ya que, a diferencia de
otros trabajos,7161728313233343536373839404142 la función no es analítica.
técnicas numéricas para calcular las integrales en la Ec. (9). La integral interna
(entre llaves), se resuelve aquí mediante un algoritmo de integración numérica La frecuencia angular promedio utilizada para calcular α y β a través de la ecuación. (5)
('integral.m') en MATLAB,52 basado en un método de cuadratura adaptativa (en se evalúa numéricamente como un promedio ponderado (centroide espectral)dieciséis en la
lugar de, por ejemplo, el método de número de onda discreto, DWM,53 utilizado ecuación. (10):
por Blair17). Esto conlleva una reducción significativa del tiempo computacional,
∫∞ ω|S(ω) |dω
ya que no manejamos matrices grandes como con el DWM. Dado que el resultado ωa = ∫0 (10)
∞ |S(ω) |dω
es una señal discreta, se necesita una frecuencia de muestreo lo suficientemente 0

grande para capturar picos con suficiente precisión. Dicha frecuencia depende en La frecuencia angular promedio para la función fuente en Figura 2 es
gran medida de la función de la fuente transitoria; En el caso estudiado, se eligió 77,5 kHz.
una frecuencia de muestreo de 0,7 MHz debido al breve tiempo de subida de la La presión de von Neumann (PAGVN) se utiliza a menudo para escalar la fuente
función de fuente descrita en la sección2.3. Para la integral externa se debe función (ver p. ej.16-18,24); sin embargo, este valor es a menudo especulativo, ya
seleccionar un número suficiente de puntos para incluir hasta la última ola que que generalmente no se dispone de mediciones sobre él; además, el estado de
llega al punto de monitoreo (ver von Neumann pertenece al explosivo chocado sin reaccionar, y no es un estado
Cuadro C2). Los límites de integración se eligen para que el campo sísmico de choque de la roca. Esto nos ha llevado a utilizar un método alternativo que
utiliza información comúnmente disponible para un explosivo con una
caracterización de detonación razonable; la presión que actúa sobre la roca se
deriva de un cálculo de coincidencia de choque entre explosivo y roca. Esto se
resume enapéndice B y su validez se evalúa en la sección 3.3. La función fuente en
Figura 2 se utiliza en este trabajo como fuente radial y torsional para producir
radiación de ondas P-SV y SH, respectivamente, a pesar de que el barreno no
actúa como una fuente de torsión (como el

tabla 1
Propiedades mecánicas del macizo rocoso.

Propiedad Valor

Velocidad de la onda P, Vpag (m / s) 5844

velocidad de la onda S, Vs (m / s) 3300.4
Densidad, ρ (kg / m3) El 2721.2
coeficiente de Poisson, ν 0,24
Módulo joven, mi (GPa) Módulo de cortea, μ (GPa) 74,76
El primer parámetro de LameB, λ (GPa) Resistencia 29,64
a la compresión uniaxial, UCS (MPa) Resistencia a la 33,65
tracción, TS (MPa) 171,21
18.21

a Obtenido de μ = ρV2 s.

Figura 2. Función de fuente normalizada para explosivos en emulsión a partir de B Obtenido de λ = ρV2 pag-μ.

modelos numéricos.

4
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Figura 4. Señales procesadas de A3 en la explosión B3 (izquierda) y G1 en la explosión B9 (derecha).

Tabla 3
D, ρmi y PAGSM valores utilizados como entradas en la sección 3.2.

Explosión D (Sra) ρmi (kg / m3) PAGSM (GPa)

Fig. 3. Disposición de los sensores de monitorización y explosiones.

3 4955,6 1154 9.3941

6 5306 1156 10.5328
7 5254 1152 10.3387
9 5956,4 1300 13.8072
barreno no gira sobre sí mismo). La razón para considerar una fuente
de torsión distinta de cero es investigar la influencia del componente
transversal en el resultado global. Se pueden aplicar dos casos: función
fuente con módulo y dirección aleatorizados, o función fuente con Se utilizaron tres acelerómetros de cizallamiento cerámico uniaxial, A1 (10000

módulo y dirección constantes. El primero produce una contribución g), A2 (5000 g) y A3 (500 g), y un geófono uniaxial, G1 con un rango de 2000 mm /

insignificante, ya que no se produce interferencia constructiva, por lo s. Todos los sensores se orientaron verticalmente y se colocaron en pares en el

que se aplica el segundo. fondo de un pozo vertical de 6 m detrás de las explosiones. A2 / A3 compuso un
par y A1 / G1 el otro. El par A2 / A3 se fijó a un bloque de aluminio atornillado a

3. Calibración y validación de modelos otra pieza de aluminio taponando el extremo de un tubo rígido de PVC de 63 mm;
Tanto el tapón como la parte inferior del tubo se enlecharon a las paredes del

3.1. Descripción del sitio de prueba y recopilación de datos

pozo asegurando una unión firme de los sensores al macizo rocoso, permitiendo
al mismo tiempo recuperar los sensores. Se colocaron primero en un pozo para

Se llevaron a cabo una serie de pruebas de voladuras a gran escala en la monitorear las voladuras B3 y B6 (A2 / A3), y en otro para monitorear B7 y B9

cantera de El Aljibe entre 2017 y 2019. La cantera está ubicada cerca de Almonacid (A2b / A3b en

de Toledo (España). Extrae mylonita, una roca metamórfica de alta resistencia que Fig. 3). A1 / G1 se embebieron en un bloque de lechada y se lecharon

se utiliza como balasto ferroviario convencional y de alta velocidad, y como directamente (sin ningún dispositivo intermedio) al fondo de otro pozo de 6 m de

agregados para la construcción. Las voladuras se ubicaron en la parte suroeste longitud. Este acoplamiento no permitió la recuperación de los sensores. Un

del nivel inferior del tajo, fuera de la zona de producción normal de la cantera. La diseño de las explosiones y los sensores se muestra enFig. 3; en la voladura B9, se

estructura de la roca se puede describir como en bloques. Las propiedades perforó el pozo de monitoreo, en el que se colocó el par A1 / G1, entre los pozos

mecánicas de la roca se determinaron a partir de cinco muestras recolectadas de voladura BH # 4 y BH # 5.

cerca del área de voladura y se muestran entabla 1. Marcadores magenta en Fig. 3 denotar la posición de los orificios de control.

Se instrumentaron cuatro voladuras, cada una compuesta por una fila de Las líneas de colores representan los pozos de cada explosión, siendo los círculos

siete barrenos, con un diámetro de 89 mm, inclinados 15◦–20◦ desde la vertical las posiciones de sus cuellos. Coordenadas locales (r, z) y longitudes de carga,

para hacerlos paralelos a la cara libre del banco. El explosivo utilizado en las CL, calculado utilizando el esquema de Figura 1, se resumen en Tabla 2. Las señales
sísmicas se alimentaron a dos unidades de adquisición MREL DATATRAP II utilizando
pruebas fue la emulsión cuyo historial de presión-tiempo se muestra enFigura 2.
la frecuencia de muestreo máxima disponible, 10 MHz. Las unidades se activaron
El explosivo se inició en el fondo con impulsores de 400 gy detonadores
mediante una tapa electrónica que se disparó simultáneamente con el orificio principal
electrónicos. Las explosiones se realizaron desde el pozo de voladura BH # 7,
BH # 7.
ubicado en la parte noroeste del bloque, hacia el pozo de voladura BH # 1. Se
Las formas de onda medidas, las señales de velocidad y aceleración se pasaron a
utilizó un retardo de 23 ms entre agujeros.
través de un algoritmo de eliminación de tendencias para eliminar los sesgos y minimizar

Tabla 2
Coordenadas locales y longitudes de carga.

BH Explosión 3 Explosión 6 Explosión 7 Explosión 9

A2 / A3 A2 / A3 A2b / A3b A1 / G1 A2b / A3b

r z CL R z CL r z CL r z CL r z CL

1 24,4 7.7 11,9 18 3,5 11,6 28,3 5.8 11,1 14.3 2.4 11,7 28,5 5,6 11,7
2 21,7 8.2 10,6 15,1 4.2 11,4 26,2 5.4 10,3 11 2.6 12 25,8 5.4 12
3 19,3 7.1 11 12,1 3.8 11,1 23,5 8.2 12,4 7.1 2.6 12,4 22,4 5.9 12,4
4 17.3 6,7 11,2 8.8 3.8 10,9 20,9 7.3 12,9 3.6 2.6 12,8 19,5 4.9 12,8
5 15,7 5.7 10,7 7 3,7 10,9 18,2 6 11,9 2,7 1,9 11 16,3 3.9 11
6 14,8 6.5 11,6 5.4 3.1 11,1 15,9 5.7 12,2 4.4 2.1 12,6 13,6 4.4 12,6
7 14,8 7.0 11,9 7.1 3.8 11,4 13,7 9,7 2,7 8.4 2,7 12,8 12,3 5.1 12,8

5
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

las horas extraordinarias del cambio de línea de base. A continuación, se procesaron con una entre los límites inferior y superior de un intervalo de tiempo asociado con la
ventana móvil Tukey (coseno cónico) con un ancho de banda de 23 ms y una fracción de coseno onda P o S. Dado que la onda S no se puede distinguir correctamente en las
del 2,5% para aislar las señales procedentes de cada orificio, y mediante un filtro Butterworth señales registradas, las velocidades de las ondas P y S medidas en el laboratorio
numérico de paso de banda de fase cero. Como ejemplo,Figura 4 (vertabla 1) se utilizan para establecer los tiempos de llegada teóricos como se
muestra dos formas de onda procesadas. muestra en Apéndice C. Cuadro 4 muestra las permutaciones de llegada de las
En las pruebas también se midieron la velocidad de detonación y la densidad del ondas sísmicas a los puntos de monitoreo para las condiciones de prueba.
explosivo. Estos valores se resumen enTabla 3. La presión de coincidencia de choquePAG También se muestran los intervalos de tiempo elegidos. Los subíndices P y S se
SM, obtenido como se muestra en apéndice B, también se da. utilizan para cada tipo de onda, M para la onda sísmica de Mach y U y L para la
fase de descarga o carga.18
Como ejemplo, Figura 5 muestra los intervalos de tiempo para las permutaciones que
3.2. Calibración de valores Q
ocurren en las condiciones de prueba.
Las principales estadísticas del resultado Qw de nuestras señales se muestran en
Las señales sísmicas registradas se utilizan para calibrar el modelo en
Cuadro 5. En términos mediosQpag es el doble de Qs valor. Una prueba de suma de
amplitud ajustando el Q-valores para que las aceleraciones o velocidades
rangos de Wilcoxon de dos caras56 se ha aplicado para evaluar si Qpag y Qs (P-SV) o
pico calculadas se ajusten a los valores medidos. Se consideran dos casos, el
Qpag y Qs (acoplados) provienen de distribuciones con medianas iguales. Lo
problema P-SV solo y el acoplado. Para ello, se ha escrito en MATLAB un
bajopag-valores (6 × 10-6 y 1,36 × 10-6) son una evidencia en contra de la
código basado en un algoritmo de búsqueda de raíces.52 Se ha seguido el
hipótesis de medianas iguales. Si se aplica la misma prueba aQs (P-SV) y
siguiente procedimiento:
Qs (acoplado), el alto pag-valor en este caso (0.545) indica que la hipótesis de que
Qs (P-SV) y Qs (acoplados) provienen de distribuciones con medianas iguales no
1. Coordenadas locales (r, z), longitudes de carga, CLy cosenos de director, γζ, se
pueden rechazarse. La razón de esto es que la componente transversal obtenida
obtienen para cada barreno a partir de las coordenadas absolutas del collar
con una función de fuente de torsión constante tiene una componente vertical
(C) y del fondo (B) y las longitudes de derivación (ver Figura 1 y Tabla 2).
despreciable, incluso con agujeros inclinados. Esto se puede ver enFigura 6.
2. Las horas de llegada se calculan como se describe en Apéndice C utilizando las
propiedades mecánicas del macizo rocoso en tabla 1 y coordenadas locales y
Figura 6 muestra la variación de la resultante Qw con la distancia
longitudes de carga en Tabla 2.
radial. Para describirlos se considera el siguiente modelo de dos
3. Aceleraciones y velocidades máximas de partículas experimentales, ppa y ppv,
coeficientes:
asociadas con cada velocidad de propagación se calculan utilizando los
tiempos de llegada en las señales registradas Qw = a1(1 - mi- a2(r- a)) + 0.5 (12)
4. Se realiza un cálculo de coincidencia de impactos como se describe en
La razón del coeficiente 0.5, que obliga a la Qw valor en la pared del
apéndice B para calcular la presión máxima en la pared del pozo; Se utilizan
pozo, se explica en la sección 3.3.
las propiedades explosivas, la densidad y las constantes de relación lineal de
Para evaluar los coeficientes del modelo, valor atípico Qpag y Qs los valores se
la velocidad de choque a partícula de la roca.
descartan. Las funciones ajustadas (Eq.(12)) están graficados en Figura 6 por Qpag y
5. P y S Q-Los valores se calculan para cada ubicación de barreno de modo que las
Qs. Aumentan monótonamente con la distancia a un valor constante.a1 +
aceleraciones y velocidades calculadas coincidan con los valores medidos.
0,5 en el campo lejano. Las principales estadísticas de los ajustes se muestran enTabla 6.
6. Dos exponenciales Q-Se ajustan las leyes dependientes de la distancia
Los niveles de significancia (pag-valores) de los coeficientes también se dan. A pesar de
radial.
la dispersión, todos están por debajo de 0,05.

El algoritmo de búsqueda de raíces secantes se ha aplicado porque la

derivada analítica de cualquier campo sísmico (ver Ec. (9)) con respecto a Qpag y 3.3. La condición de la pared del pozo de voladura

Qs implica un trabajo laborioso. El método está definido por la relación de

recurrencia: Para evaluar la influencia de usar como presión pico el pico de von Neumann,
estimado como 2.5 veces la presión de Chapman-Jouget (PAGCJ),57 o la presión de
( ) Qnorte- 1, w - Q(norte- 2, w )
(
Qnoroeste = Qnorte- 1, w - gramo Qnorte- 1, w gramo Qnorte- 1, w -) (11) coincidencia del choque (PAGSM), la velocidad de la partícula (ppv) de los puntos
gramo Qnorte- 2, w
de coincidencia para velocidades de detonación entre 4000 y 6500 m / sy
densidad explosiva constante (1178 kg / m3), se han calculado siguiendo la
dónde norte es la iteración actual y w el tipo de onda, es decir, onda P o S, g (
metodología de apéndice B. Están graficados enFigura 7
Qn, w) = |ppxw - ppxw* |, ppxw siendo la velocidad máxima medida (X = v) o
(cruces negras). Se comparan con losppv obtenido del FFS en la pared del
aceleración (X = a) del tipo de onda P o S de cualquier barreno y
pozo de voladura para el problema P-SV, utilizando tanto PAGSM (líneas
ppxw* el valor respectivo obtenido con Eq. (9) en las coordenadas (r, z) del
continuas) y PAGVN (líneas discontinuas), para variar Qw de 0,5 a 50.
punto de seguimiento de la ola procedente de cada hoyo. El número de
En el rango de las condiciones de prueba, que se muestran en Figura 7, velocidades máximas
iteraciones necesarias es relativamente bajo, alrededor de diez, aunque la
del cálculo de coincidencia de choque están en buen acuerdo con los valores obtenidos
tasa de convergencia del algoritmo utilizado es menor que otros métodos.
con el FFS utilizando PAGSM cuando Qw es aproximadamente 0,5. Es por eso que esta
Eq.(11) continúa hasta g (Qw) <ε, dónde ε es de 0,01 mm / so
constante se incluye en la ecuación.(12), forzando Qw ser 0,5 en la pared del barreno. El
0,01 g.
uso dePAGVN sobreestima las velocidades máximas en la pared del pozo de voladura.
Los niveles máximos ppxw y ppxw* se evalúan como el máximo de
Figura 7 también muestra que ppv aumenta con D para todos Q-valores en el
valor absoluto de la señal medida o predicha, respectivamente
caso transónico, es decir, cuando Vs <D <Vpag. Sin embargo, y contrariamente a la
creencia común, comoD enfoques Vpag, los niveles de vibración descienden Q
Cuadro 4
-valores. Esto no ocurre paraQ < 5.
Intervalos de tiempo en función de la permutación a utilizar.

Permutación Orden de llegada Intervalos de tiempo

3.4. Evaluación de las vibraciones máximas previstas
P-Wave S-Wave

2 tPL, tPU, tSU, tSL (tPL, tSL) (tSL, tmax) Los niveles máximos de vibración asociados a cada agujero se han calculado
3 tPL, tSRA, tSL, tPU, tSU (tPL, tSRA) ∪ (tPU, tSU) ( (tSRA, tPU) ∪ (tSU, tmax) (
utilizando la Ec. (9) para voladuras B3, B6, B7 y B9. Datos enTablas 1-3 han sido
4 tPL, tSRA, tPU, tSL, tSU tPL, tSRA) ∪ (tPU, tSL) (t tSRA, tPU) ∪ (tSL, tmax)
considerados y para cada conjunto de datos, los Qw han sido evaluados usando la
6 tPL, tPU, tSRA, tSL, tSU PL, tSRA) (tSRA, tmax)
Ec. (12). Como ejemplo, la forma de onda completa medida por A3 en la explosión
B7 y la obtenida por medio del FFS se muestran enFigura 8. La mitad-

6
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Figura 5. Tiempos de llegada teóricos para señales

individuales y sus intervalos de tiempo utilizados para
calcular ppxw. Las flechas rojas y azules indican los
intervalos de ondas P y S que se consideran para evaluar
ppxw. Tenga en cuenta quetSRA, tPU y tSL
prácticamente se superponen cuando se produce
la cuarta permutación. En caso de que las ondas P
y S se superpongan, el código solo ajustaQpag, ya
que ajustando ambos Qpag y Qs se convierte en un
problema en mal estado. . (Para la interpretación
de las referencias al color en la leyenda de esta
figura, se remite al lector a la versión web de este
artículo).

Cuadro 5 Tabla 6
Media, desviación estándar, 25th percentil, mediana y 75th percentil para Qpag y Análisis de regresión para Q-Modelos de ajuste.
Qs. a1 a2 pag (a1) pag (a2) NRMSEa
Parámetro Significar Mediana
Qpag 49,64 0,12 <0,001 0.0182 0.362
Estándar PAG25 PAG75

Qpag 40,62 20.24 27,46 37,61 55,29 Qs (P-SV) 32,74 0.054 0,0012 0.0479 0,247
Qs (P-SV) 19,79 6.53 14.51 19,88 23.58 Qs (acoplado) 27,73 0,065 <0,001 0.0381 0,259
Qs (Acoplado) 18,69 6.07 14.07 19.31 22.47
a Error cuadrático medio normalizado.

Figura 7. Velocidad máxima de las partículas en la pared del barreno en función de la velocidad
Figura 6. Q-leyes para Qpag (azul), Qs (P-SV) (rojo) y Qs (acoplado) (verde). (Para la
de detonación. Las cruces negras son el resultado del cálculo de coincidencia de impactos. Las
interpretación de las referencias al color en la leyenda de esta figura, se remite al
líneas son los resultados del FFS en la pared del barreno utilizando como pico de la función de
lector a la versión web de este artículo).
fuentePAGSM (líneas continuas) y PAGVN (líneas discontinuas) para diferentes Q-valores.

7
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Figura 8. Comparación entre la señal medida, A3 de la explosión B7 (izquierda) y

calculada a partir de FFS (derecha).

Figura 10. Patrones de radiación para el problema P-SV.

Se nota que los valores calculados con el FFS a distancias escaladas bajas (A1
en B9) son siempre superiores a los medidos (asteriscos verdes). La razon es queQ
-valores obtenidos de la ecuación. (12) para ese rango de distancia, son
sistemáticamente mayores que los obtenidos mediante la Ec. (11) utilizando los
picos medidos (ver Figura 6). En las mediciones con A1 en la voladura B9 esto
también podría explicarse por el hecho de que el macizo rocoso había sido
dañado por las olas provenientes de los pozos anteriores; Además, la densidad
explosiva y la velocidad de detonación fueron particularmente altas en la
explosión B9 (verTabla 3) provocando un valor más alto de PAGSM. Esto implica
que valores más bajos deQ (es decir, peor capacidad de transmisión) para
alcanzar amplitudes similares a las de las señales medidas. La mayoría de los
datos utilizados para calibrar la ecuación.(12) se obtienen a distancias radiales
superiores a 8,5 m donde no se aprecian daños.
Figura 9. Aceleraciones máximas medidas y pronosticadas frente a la distancia escalada.

4. Aplicación de muestra

Se considera un barreno genérico de 10 m compuesto por 11235 elementos con una

Tabla 7
longitud de una centésima parte de su radio (44,5 mm). La distancia radial máxima es de
Propiedades explosivas utilizadas en el estudio de muestra.
6 m, ya que se espera que la roca no se rompa debido a las ondas del cuerpo a una
ρmi (kg/metro3) D (Sra) PAGSM (GPa)
distancia mayor. La malla se extiende 5 m por encima de la parte superior y 1 m por
1193,93 5534.4 11.542 debajo del fondo del barreno, dando un dominio de 16×
6 m. Se utilizó una malla computacional de 10 cm cuadrados, con un total de 9821
puntos de cálculo. La velocidad de detonación y la densidad explosiva elegidas
son la media de todos los valores disponibles obtenidos en las pruebas (verTabla 3
longitud de los elementos de carga L se ha ajustado a una centésima parte ) y se dan en Tabla 7. Las propiedades de las rocas son lastabla 1. Q-los valores se
del radio del barreno (44,5 mm), siguiendo la recomendación de que L < calculan a partir de la ecuación. (12). Los patrones de velocidad de radiación y
0,2a debido a la fase negativa provocada por las ondas de descarga de los aceleración del P-SV se muestran como mapas de contorno enFigura 10. El tiempo
extremos de los elementos.17 de cálculo varía entre 0,5 sy 2 s por nodo y componente de campo (velocidad o
Para obtener el pico de vibración de cada barreno, las señales de aceleración se aceleración) según la longitud de la señal considerada y el campo evaluado, en un
aíslan aplicando una ventana rectangular de ancho de banda de 23 ms. Se trazan frente portátil MSI GF62 8RD de 2,2 GHz.
a la distancia escalada (Dakota del Sur o relación de la distancia al centro de gravedad de La velocidad y la aceleración máximas totales de las partículas se calcularon a partir
un barreno dividido por la raíz cuadrada de la masa explosiva de ese barreno), en escala de los componentes radial y axial como se muestra en la Ec. (13):
logarítmica, en Figura 9. Los valores medidos también se representan como referencia. √̅ ̅
Un ajuste lineal de los valores máximos medidos frente aDakota del Sur en escala ppx =max X2 r + X2z (13)
t
logarítmica del tipo ppa = a1Dakota del Sura2 fue hecho; También se trazan las bandas de
confianza del 95% (líneas discontinuas de color magenta) de la regresión y las nuevas Es interesante notar que el patrón de radiación para la aceleración
predicciones (líneas negras discontinuas). Las señales de velocidad medidas por G1 no se muestra un doblez en la parte central-superior del gráfico. La razón de tal
han utilizado para el ajuste, ya que solo están disponibles cinco puntos.Figura 9 muestra resultado es la interferencia entre la onda S Mach y las ondas PU (ver línea
que los valores predichos están bastante cerca de la línea de ajuste para los valores verde enFigura C2) en el ejemplo del régimen transónico.
medidos y, en la mayoría de los casos, dentro de los límites de confianza de la regresión. El campo de desplazamiento se ha utilizado para obtener el campo de tensión, a
partir del cual derivar un patrón de daño. Dos criterios de falla basados enUCS, TS
y las tensiones principales producidas por la vibración radiada, se han

8
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

utilizado para predecir la extensión de las zonas de fractura. Estos son el criterio
de tensión máxima y el criterio de rendimiento de Drucker-Prager.58,59. Este
último se ha elegido porque el modelo considera un medio isótropo. La
convención de signos adoptada para deformaciones es positiva para
alargamientos y negativa para contracciones. Esto implica un signo positivo para
las tensiones de tracción y uno negativo para la compresión. El criterio de tensión
máxima asume que un material falla cuando la tensión principal máxima (σ1) en
un elemento material excede TS. Alternativamente, el material también falla si la
tensión principal mínima (σ3) es menos que UCS. Entonces se supone que la
región "segura" para el material es:

- UCS <σ3 <σ1 < TS (14)

El criterio de rendimiento de Drucker-Prager conduce a una región "segura" para el
material definida por:
( )
1 1
- (σ1+ σ +2 σ ) 3
TS UCS
1 [ 2 2
2]
+ (σ1- σ2) + (σ2 - σ ) + (3σ3 - σ1) < 1 (15)
2⋅TS⋅UCS
Apéndice D muestra el procedimiento para calcular los componentes del tensor de
tensión a partir del campo de desplazamiento a través de la Ec. (D.1) a (D.6).
Figura 11 muestra los principales campos de tensión a 0,214 ms, 1 ms y 2,143 ms
respectivamente.
Trazar puntos donde UCS y TS se superan, se pueden identificar tres
zonas de rotura (ver Figura 12). El primero (Crush Zone, CZ) es donde el
valor absoluto deσ3 excede UCS. En esa zona la roca se rompe por
compresión. El segundo (Zona de grietas por tracción, TCZ) es donde
σ1 excede TS. En esa zona la roca se rompe bajo tensión. Están basados en la
ecuación.(14). La tercera zona (Zona Drucker-Prager, DPZ) es donde se excede el
criterio de rendimiento Drucker-Prager. Más allá de eso, no hay daño a la roca por
la onda sísmica. Se basa en la ecuación.(15). CZ se extiende 0.36 m (ocho veces el
radio del pozo) a lo largo de todo el pozo y el TCZ se extiende un valor constante
de 1.13 m en la parte central del gráfico a medida que avanza la detonación, y
Figura 11. Componentes principales de la tensión en los tiempos 0,214 ms, 1 ms y 2,143 hasta 2.03 m en la parte superior. zona del agujero debido a la interrupción de la
ms. Los colores rojo y azul indican tensión de tracción y compresión, respectivamente. detonación (ver gráfico superior derecho enFigura 11). El DPZ se extiende hasta
(Para la interpretación de las referencias al color en la leyenda de esta figura, se remite aproximadamente 1,5 m en la parte central del agujero hasta 2,5 m en la parte
al lector a la versión web de este artículo). superior.

5. Discusión: comparación con modelos de ingeniería comunes

A pesar de los defectos obvios en su formulación 7,60, el Holmberg-Persson

El enfoque de la radiación sísmica de una carga lineal es ampliamente
utilizado por ingenieros y consultores en todo el mundo, p. ej.
8,9.
Creemos que la razón fundamental para utilizar el modelo HP es su
enorme sencillez. Sin embargo, teníamos curiosidad por saber cómo se
comparan sus predicciones con un modelo más exacto y físicamente sólido
como el que se presenta aquí que, a su vez, es más complejo y requiere
habilidades computacionales no despreciables. Además del modelo de HP,
consideraremos el modificado por NIOSH.5
El enfoque HP utiliza la ecuación básica de propagación de vibraciones
desde un barreno:

ppv =K⋅Wα/ Rβ (dieciséis)

Dónde K, α y β son constantes específicas del sitio (α y β se utilizan para preservar

la notación de uso común, que no deben confundirse con las velocidades de onda
complejas en la sección 1.1), W es la carga explosiva total (kg) y R es el

Tabla 8
Constantes para funciones de atenuación HP y NIOSH, Ecs. (16) y (20). R (metro), Q
Figura 12. Zonas de daño. CZ (rojo), TCZ (amarillo), DPZ (azul claro) y sin daños (kg), ppv (mm / s).
(azul oscuro). (Para la interpretación de las referencias al color en la leyenda de
Acercarse α K ARCE
esta figura, se remite al lector a la versión web de este artículo).
HP 0,752 218.613 23,68
Modificado por NIOSH 1.203 187.521 23,77

9
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Figura 13. Mapas de contorno de velocidades máximas de FFS (izquierda), el modelo HP (centro) y el enfoque modificado por NIOSH (derecha). Isolíneas deppv están etiquetados en m / s.

Figura 14. Mapas de contorno llenos de diferencias relativas para el modelo HP (gráfico de la izquierda) y para el modelo modificado por NIOSH (gráfico de la derecha) con respecto a FFS.

distancia (m). Dividiendo la columna explosiva en segmentos y sumando la dónde q es la densidad de carga lineal (kg / m); yzI y zF son las coordenadas
contribución de cada elemento de explosivo a la vibración en un punto P (r0, z0), axiales de la parte inferior y superior de la carga explosiva,
Eq. (dieciséis) se convierte en: respectivamente. La solución de la integral en la ecuación.(17) es:
⎡ ⎤α
∫ zF
dz ⎦ ∫ zF ( )⃒
ppv =K⎣q (17) dz 1 β 3 (z - z )2 ⃒zF 0
(r02 + (z - z0)2)
β /2α
= r-0β / α(z - z0)2F1 2 2α, 2 ; ; - ⃒ (18)
zI β /2α r2 ⃒
zI (r2 0 + (z - z0)2) 0 zI

10
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

distancias reales. Existe una zona justa con patrones similares para ambos
dónde 2F1(a,B;C;s) es la función hipergeométrica de Gauss definida en.61
modelos, coaxial con la carga hasta 1–2 m de distancia radial, ensanchándose
Configuración β = 2α por simplicidad, Eq. (18) se reduce a Eq. (19):
hacia arriba y hacia abajo. Más allá de la zona justa, ambos modelos sobreestiman
∫ ( )⃒
zF
dz 1 z - z0 ⃒z⃒ F las vibraciones en una especie de 'islotes de error' donde los errores aumentan a
norte- 1
(19)
r0 r0 ⃒ zI
= ejército de reserva
β /2α
zI (r20 + (z - z0)2) casi 3 veces para el modelo HP y 4 veces para el modificado por NIOSH.

El enfoque modificado por NIOSH reescribe la ecuación. (dieciséis) como sigue:

6. Conclusiones

ppv =K⋅Wα/ Rβ (20)

El FFS de la radiación de ondas sísmicas de un barreno, a pesar de su solidez
física, no ha sido utilizado a fondo por los ingenieros debido a su relativa
dóndeRes la distancia promedio al punto de observación para todos los
complejidad, requiriendo mayores habilidades matemáticas y de programación y
elementos explosivos; se puede calcular mediante la ecuación.(21):
tiempos computacionales más largos que los modelos más simples, aunque
∫ √̅z̅ F
1 físicamente cuestionables. El FFS ha sido reevaluado en este trabajo, reformulado
R= r02 + (z - z0)2dz (21)
CL zI y resuelto por un nuevo método de cálculo que, a través de una transformación
de la ecuación de campo total (Eq.(8)), reduce significativamente el esfuerzo
Con el fin de proporcionar un escenario en el mejor de los casos para la
computacional, permitiendo así su implementación en un portátil y su uso como
comparación de los enfoques de HP y los modificados por NIOSH con el FFS, las
herramienta de ingeniería.
constantes específicas del sitio en la Ec. (dieciséis) y Eq. (20), con β = 2α se estiman
Las mediciones de vibración de las pruebas de voladura se han utilizado para
mediante un ajuste lineal en escala logarítmica utilizando las velocidades máximas
desarrollar una metodología para calibrar el FFS, notablemente la atenuación
calculadas por el propio FFS, para el cálculo de la muestra en la sección 4. Las principales
Q-factores para las ondas P y S, para las cuales se han establecido
estadísticas de los ajustes se muestran enTabla 8. En todos los casos, los coeficientes
funciones dependientes de la distancia al barreno.
obtenidos son significativos a un nivel de 0,05. La métrica de error utilizada es el Error
La salida directa del historial de presiones en la pared del pozo a partir del
Logarítmico de Porcentaje Absoluto Medio (MAPLE) y se calcula mediante la Ec.
modelado numérico de la detonación (calibrada con datos experimentales del
(22):
comportamiento explosivo) se ha utilizado para obtener la forma de la función
⃒ ⃒
100 ∑⃒ ⃒l⃒og ppvFFS - Iniciar sesión ppvHP⃒
⃒ fuente. La presión de coincidencia de choque de la detonación en la pared de roca
ARCE = ⃒ ⃒ (22)
norte FFS ppv
Iniciar sesión se usa como presión máxima en lugar del pico poco realista de von Neumann que
se usa a menudo, un valor no respaldado por la física ya que pertenece al
Tenga en cuenta que el exponente de atenuación para el modelo HP es
hugoniot explosivo sin reaccionar, y no al hugoniot de choque de la roca.
similar al valor predeterminado clásico α = 0,7 mientras K es inferior a los 700
mm / sy similar a los valores obtenidos por Onederra y Esen62 para explosivos a
El campo de desplazamiento del FFS se utiliza para obtener el campo de
granel en el campo cercano, que varió de 150 a 470 para K y 0,87 a 1,12 para α.
tensión, del cual se pueden obtener fácilmente las zonas dañadas. Criterios de
Figura 13 muestra ppv mapas de contorno para cada modelo.
falla basados enUCS, TS y se han aplicado los principales esfuerzos producidos
Para facilitar la comparación, la razón de la ppv de enfoques
por la vibración radiada.
modificados por HP y NIOSH a la ppv a partir de la FFS se calcula.
El modelo de vibración de campo cercano Holmberg-Persson clásico y
Figura 14 muestra los mapas de contorno rellenos de ppvHP/ ppvFFS (relación entre ppv
ampliamente utilizado y su modificación NIOSH se han comparado con la
de los modelos HP y FFS). El código de color se muestra a la derecha; por
solución de campo completo. Las principales diferencias con HP y FFS son la
ejemplo, en el área magenta el modelo HP o modificado por NIOSH
ausencia de direccionalidad delppv contornos en el primero. Los contornos
subestima los niveles de vibración, arrojando resultados entre 1/4 y 1/3 de
son elipsoidales para el modelo HP y casi esféricos para la modificación
los FFS; lo contrario ocurre en el área roja donde los resultados de HP son
NIOSH, ya que no hay detonación hacia arriba ni interferencias de ondas.
de 3 a 4 veces la FFS.
Los valores tres veces o un tercio de los valores de FFS son comunes
A primera vista, ambos enfoques de HP subestiman (severamente en el caso
alrededor del pozo, incluso más allá de esos límites en algunas áreas.
modificado por NIOSH) los niveles de vibración con respecto al FFS en la zona
superior desplazada del pozo. Esto se debe a que los cálculos de HP no pueden
predecir el componente dominante hacia arriba, ya que la detonación se Declaración de intereses en competencia
considera instantánea. Por el contrario, ambos modelos sobreestiman las
velocidades máximas en las zonas superior e inferior de la roca, en el eje del pozo. Los autores declaran que no tienen intereses económicos en
El enfoque modificado por NIOSH subestima gravemente las vibraciones en las competencia o relaciones personales conocidas que puedan haber influido
proximidades del barreno; Este modelo utiliza una distancia de viaje promedio a en el trabajo informado en este documento.
la posición de observación que, cerca del barreno, es mucho más larga que la
distancia real, de ahí lappv la predicción es muy baja; el criterio de distancia Agradecimientos
promedio da como resultado un patrón de radiación pseudoesférico (verFigura 13
). El modelo HP no muestra este mal funcionamiento (excepto en una sección muy Este trabajo se ha realizado en el marco del proyecto SLIM financiado por el
delgada adyacente al orificio), ya que los componentes de vibración de los programa de investigación e innovación Horizonte 2020 de la Unión Europea en virtud
elementos de carga cercanos al punto de observación se calculan en su del acuerdo de subvención núm. 730294.

Apéndice A. Ecuaciones

Los potenciales de desplazamiento son12 (vea la notación en el texto principal):

∫∫) ∞ (
1- S(ω) - ρω2 + 2μk2 z K1(mamá) ⋅2LsincLkz) K0(lr)
φ = (2π)2 ( )2 miikz zmi- es dkzdω (A.1)
- ∞ - ρω2 + 2μk2z K0(la) K1(mamá) - 2μlρω2Ka1(la) K1(mamá) - 4μ2k2 z lmK0(mamá) K1(la)

11
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

∫∫ ∞
1 - S(ω) 2μikzlK1(la) ⋅2LsincLkz) K1(señor)
ψ= ( ) miikz zmi- esdkzdω (A.2)
(2π)2 - ∞- ρω2 + 2μk 2 z2 K0(la) K1(mamá) - 2μlρω2Ka1(la) K1(mamá) - 4μ2k2 z lmK0(mamá) K1(la)

1 ∫∫ ∞ 2LS(ω) sinc (Lkz) K0(señor)

χ= ( ) miikz zmi- esdkzdω (A.3)
(2π)2 -
∞ μ metro2K0(mamá) + 2metroK1(mamá
a
)

Y el desplazamiento de un solo elemento, acoplando los problemas P-SV y SH en una ubicación P (r, z), viene dado por el vector ū = (tur, tuθ, tuz), cuyos
componentes se evalúan como:

∂φ ∂ψ ∂φ 1 ∂ (rψ) ∂χ
tur = - tuz = + tuθ = - (A.4)
∂r ∂z ∂z r ∂r ∂r

Apéndice B. Coincidencia de descargas

Presión y velocidad de las partículas en la interfaz explosivo / roca (PAGSM y vSM) están determinadas por la intersección del hugoniot directo de la roca
con el hugoniot reflejado de los productos de detonación en el punto Chapman-Jouguet (CJ). Esto se basa en las relaciones que gobiernan la ruptura de una
discontinuidad arbitraria resultante del reflejo de una onda de choque en la superficie de contacto entre dos medios (ver p. Ej.63,64). Figura B.1 muestra
gráficamente el procedimiento para calcular PAGSM y vSM.

Figura B.1. Parcela de coincidencia de choque explosivo-a-roca. Rojo: rock hugoniot. Azul: productos hugoniot. Amarillo: línea de Rayleigh. Verde: productos reflejados hugoniot.

El hugoniot de la velocidad de la partícula de presión de la roca se puede calcular mediante:

( )
PAG = ρr C + svpag vpag (B.1)

dónde ρr es la densidad inicial de la roca, 2721 kg / m3, vpag es la velocidad de la partícula (m / s) y C y s son las constantes de la relación lineal de la velocidad del choque a la partícula.

No hay muchas referencias confiables con datos duros sobre golpes de rocas competentes. Algunas referencias interesantes para granitos con propiedades
mecánicas no muy diferentes a la roca involucrada en este estudio son.65-68 Cuadro B.1 y Figura B.2 mostrar los datos obtenidos de estas referencias. En todos los casos, la
pendientes tiene un valor diferente dependiendo del rango de velocidades de las partículas con un cambio alrededor de vpag = 2000 m / s para que los datos se puedan
ajustar con dos funciones lineales por partes, en las que la velocidad de partícula más baja tiene la pendiente más pequeña.

12
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Cuadro B.1
Constantes de relación lineal de velocidad de choque a partícula para granito.

Referencias ρr (g / cm3) Vpag (km / s) Vs (km / s) C (km / s) s vpag rango (km / s)

McQueen y col. 1967sesenta y cinco 2,67 5.33 3,28 2100 1,63 2,1 - 4
Anderson y Larson 197766 2,65 6.5 ± 0,34 - 5196 0,54 0,1 - 0,6
4904 0,34 0,9 - 1,9
Marsh 198067 2.627 5.33 3,28 2428 1,54 2-4
4400 0,6 0 - 2,1
Sekine y col. 199568 2,63 5.36 ± 0,1 3.3 ± 0,2
2660 1,49 1.8 - 4.8

Figura B.2. Choque a la velocidad de la partícula hugoniot para rocas. Datos de referencias: azul,67 rojo,68 verdesesenta y cinco y negro.66 Magenta: c = 4562 m / sys = 0.6.2

Cuadro B.2
Productos de detonación coeficientes hugoniot.

A (GPa) B (Pa⋅s⋅metro-1) C (kg / m3)

3.419 890144.25 2231,75

Para las funciones de baja velocidad, la constante independiente C es la velocidad del sonido a granel.43 En nuestro caso, utilizando datos de tabla 1, C = (V2 pag-4/3Vs2 )1 2/ = 4562
m / s, un valor bastante cercano a C-valores a bajas velocidades en Cuadro B.1. La pendientes ha sido elegido igual que el de Sekine68 pendiente de rango bajos = 0,6). Esta velocidad de
choque - velocidad de partícula hugoniot se muestra (línea magenta) enFigura B.2.
Los productos de detonación hugoniot obtenidos a partir de cálculos termodinámicos con W-DETCOM69 para la emulsión utilizada, cuya densidad y velocidad de
detonación son 1177,5 kg / m3 y 5448 m / s, respectivamente, se pueden representar mediante el polinomio:

PAG = A + Bvpag + CV2 pag (B.2)

donde los coeficientes se dan en Cuadro B.2.

Los productos hugoniot reflejados a través del punto CJ son:
( ) ()
PAG = A + B 2vp, CJ - vpag + C 2vp, CJ - v 2pag (B.3)

dónde vp, CJ es la velocidad de la partícula en el punto CJ. Intersección de la ecuación.(B.1) con Eq. (B.3) da el estado de la interfaz, es decir, la presión de la pared del barreno y
velocidad de partícula. La presión en el punto CJ se puede calcular a partir de la teoría simple70 por:
/
PAGCJ = ρ0D2 (Γ + 1) (B.4)

dónde:

• ρ0 es la densidad inicial del explosivo.

• D es la velocidad de detonación.
• Γ es el coeficiente adiabático, igual a 3.4 del estado CJ conocido (PAGCJ = 7,9404 GPa, D = 5448 m / s).

Dado que diferentes valores de D y ρ0 se midieron en las diferentes explosiones (aunque se utilizó el mismo explosivo), el hugoniot (Eq. (B.2)) debiera ser
permitido variar con D y ρ0. Dejar:

13
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

PAG = A'+ B'vpag + C'v2 pag (B.5)

ser los productos de detonación hugoniot en D′ y ρ0'. A′ es la presión hugoniot cuando vpag = 0 y tú = tú0. La ecuación de energía es.70:

mi(P, υ) - mi(PAG0, υ0) = 1 / 2 (PAG + PAG0) (tú0 - tú) + Qv (B.6)

dónde mi es la energía interna de los productos de detonación. Qv es el calor de la reacción explosiva, esencialmente constante para una composición explosiva constante. PAG
y PAG0 son las presiones en los estados final e inicial, respectivamente. tú y tú0 son los volúmenes específicos (inversos de la densidad) en los estados final e inicial,
respectivamente.

La energía interna se puede calcular asumiendo mi(PAG0, tú0) = 0 como.

mi ≈ Pυ / (Γ + 1) (B.7)

Por lo tanto, para tú = tú0 = 1 /ρ0.

PAG ≈ (Γ + 1)ρ0Qv (B.8)

Entonces, asumiendo Γ y el calor de la reacción Qv ser constante A es proporcional a la densidad del explosivo (ver Eq. (B.9)):
()
ρ0'
A' ≈ A (B.9)
ρ0
Para obtener B′ y C', se pueden utilizar las siguientes condiciones:
a) Conservación del momento en el punto CJ:

PAG
CJ = 'A' + B'v' p, CJ + C'v' p, CJ = ρ' 0D'v'p, CJ (B.10)

b) Tangencia del hugoniot y la línea de Rayleigh en el punto CJ:

B' + 2C'v' p, CJ = ρ' 0D' (B.11)

Resolviendo el sistema de Ecs. (B.10) y (B.11),B′ y C' están:

2A' A'
B' = ρ' 0D' - C' = v«2 (B.12)
vp,' CJ p, CJ

Apéndice C. Horas de llegada y permutaciones de llegada

El tiempo que tarda la detonación en viajar a lo largo de la columna de explosivo es:

C
tD = L (C.1)
D

Y los tiempos de llegada de las ondas P procedentes de la parte inferior (tPL) y parte superior (tPU) de la columna explosiva son:
√̅ ̅ √̅ ̅
(r - a)2 + (- CL/2 - z)2
(r - a)2 + (CL/2 - z)2
tPL = tPU = tD + (C.2)
Vpag Vpag

mientras que el tiempo de llegada de las ondas S desde el fondo (tSL) y parte superior (tSU) de la columna explosiva son:
√̅ ̅ √̅ ̅
(r - a)2 + (- CL/2 - z)2
(r - a)2 + (CL/2 - z)2
tSL = tSU = tD + (C.3)
Vs Vs
L y U denotan las fases de carga y descarga de la ola correspondiente 18,24.
Las ondas de Mach se generan cuando hay interferencia constructiva, lo que ocurre cuando el número de Mach de detonación METROw = D/Vw es mayor que uno. De
Figura C.1, se puede determinar el tiempo de viaje de las ondas Mach. Si un punto de control P (r, z) está dentro de la zona definida por z ∈ [zinf, zsorber], entonces el
tiempo de viaje de la onda Mach consiste en el tiempo de detonación desde la parte inferior de la carga hasta un punto M y luego como una onda Mach con un vector
normal en el φcw dirección al punto P (r, z). La ecuación del tiempo de viaje está dada por Meredith12 (ver Figura C.1):

14
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Figura C.1. Tiempo de viaje para las ondas Mach (modificado de Meredith12).

z + CL/2 + (r - a)broncearse φcw

tMw = t1w - t2w + t3w = (C.4)
D

dónde:
√̅ ̅
broncearse φ 2w- 1
cw = METRO (C.5)

Figura C.2 muestra los lugares geométricos de los puntos que satisfacen tPU = tSL (cian), tPU = tPL (amarillo), tSU = tSL (naranja), tPL = tSU (magenta), tPL = tSRA (negro y tPU = tSRA
(verde) para el transónico y supersónico29 casos con un pozo de 10 my velocidades de propagación Vpag = 5884 m / sy Vs = 3300,4 m / sy velocidades de detonación de
4500 m / s (izquierda) y 7000 m / s (derecha). Estas líneas son los límites entre las permutaciones de llegada. La curva cian coincide con el bucle exterior de un óvalo de
Descartes, que se convierte en un limaçon de Pascal cuandoD = Vs. Las curvas naranja, amarilla, verde y negra son hipérbolas. Las líneas rojas y azul oscuro son los límites
donde se desarrollan las ondas Mach. Gráficos deFigura C.2 también son útiles para elegir la ubicación de los puntos de seguimiento. Cuadro C.1 muestra el orden de
llegada de las ondas provenientes de la carga explosiva en números romanos dependiendo de la posición del punto de monitoreo, que corresponde a una de las
permutaciones de llegada indicadas en Figura C.2. Las permutaciones de llegada en las condiciones de prueba de este trabajo están marcadas con asteriscos enCuadro
C.1 Las permutaciones de llegada para cada hoyo y voladura se muestran en Cuadro C.2. Figura C.2 puede orientar el posicionamiento de los puntos de monitoreo de
manera que se ubiquen lejos de los límites, a fin de que los tiempos de llegada de diferentes olas se puedan detectar claramente en los registros de vibraciones.

15
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Figura C.2. Límites entre las permutaciones de llegada para los casos transónico (izquierda) y supersónico (derecha) considerando un barreno de 10 m.

Cuadro C.1
Ranking de tiempos de llegada de diferentes olas al punto de monitoreo para diferentes permutaciones.

Permutación tMP tSRA tPL tSL tPU tSU

1 - - I II III IV
2(*) - - I III II IV
3(*) - II I III IV V
4(*) - II I IV III V
5 - II I V III IV
6(*) - III I IV II V
7 - III I V II IV
8 - - I IV II III
9 I III II IV V VI
10 I II III IV V VI
11 I II III V IV VI
12 I III II V IV VI
13 I IV II V III VI
14 I II III VI IV V
15 I III II VI IV V
dieciséis I IV II VI III V
17 I II IV VI III V
18 I II V VI III IV
19 I III IV VI II V
20 I III V VI II IV
21 I IV III VI II V
22 I - II V III IV
23 I - III V II IV
24 I - IV V II III
25 - - III IV I II
26 - - II IV I III
(*) Permutaciones de llegada en las condiciones de prueba de este trabajo.

dieciséis
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

Cuadro C.2
Permutación de llegada para cada hoyo y voladura.

BH Explosión 3 Explosión 6 Explosión 7 Explosión 9

A2 / A3 A2 / A3 A2 / A3 A1 / G1 A2 / A3

1 2 2 2 2 2
2 2 2 2 6 2
3 2 2 2 3 2
4 2 6 2 3 2
5 2 6 2 3 2
6 2 3 2 3 6
7 2 6 6 4 6

Apéndice D. Relaciones tensión-tensión

Las relaciones deformación-desplazamiento en el sistema de coordenadas cilíndricas para el problema P-SV son:

∂tur tu ∂tuz
εr = εθ= r εz = (D.1)
∂r r ∂z
( )
1 ∂tur + ∂tuz
εrθ = 0 εzθ = 0 εrz = (D.2)
2 ∂z ∂r
Aplicando la ley de Hooke, el estado de tensión para el problema P-SV es:

pagr = λΔ + 2μεrpagθ = λΔ + 2μεθpagz = λΔ + 2μεz (D.3)

pagrθ = 0pagzθ = 0pagrz = 2μεrz (D.4)

donde Δ es la dilatación, es decir, el cambio volumétrico, dado por el primer invariante principal del tensor de deformación de Cauchy:

Δ = εr + εθ + εz (D.5)

El tensor de tensión del problema P-SV desacoplado es:

⎛⎞
pagr 0 pagrz
σPAG- SV =⎝ 0 pagθ 0⎠ (D.6)
pag 0rz
pagz

Referencias 14 Tubman KM, Cheng CH, Toksöz MN. Registros acústicos sintéticos de forma de onda completa en
pozos entubados.Geofísica. 1984; 49 (7): 1051–1059.https://doi.org/10.1190/1.1441720. Meredith JA,
15 Toksöz MN, Cheng CH. Ondas de corte secundarias de los pozos fuente.
1 Ainalis D, Kaufmann O, Tshibangu JP, Verlinden O, Kouroussis G. Modelado de la fuente de
Perspectiva de Geophys. 1993; 41 (3): 287–312.https://doi.org/10.1111/j.1365-2478.1993.
voladuras para la simulación numérica de vibraciones del suelo inducidas por voladuras:
tb00571.x.
una revisión. Rock Mech Rock Eng. 2017; 50 (1): 171–193.https://doi.org/10.1007/s00603-
dieciséis Blair DP. Una comparación de Heelan y soluciones exactas para la radiación sísmica de una
016-1101-2.
carga cilíndrica corta.Geofísica. 2007; 72 (2): E33 – E41.https://doi.org/10.1190/
2 Persson PA, Holmberg R, Lee J. Ingeniería de explosivos y voladura de rocas. Prensa CRC;
1.2424543.
1993.
17 Blair DP. Radiación sísmica de una columna explosiva.Geofísica. 2010; 75 (1):
3 Holmberg R, Persson PA. El enfoque sueco para la voladura de contornos.Actas de la IV
E55 – E65.https://doi.org/10.1190/1.3294860.
Conferencia sobre explosivos y técnica de voladura. 1978: 113–127. Hustrulid W, Lu W. El
18 Blair DP. Dependencia de la vibración de la explosión en la longitud de la carga, la velocidad de detonación y los
4 enfoque Lu-Hustrulid para calcular la velocidad máxima de partículas causada por
medios estratificados.Ciencia ficción Int J Rock Mech Min. 2014; 65: 29–39.https://doi.org/10.1016/j.
voladuras. En:Técnica de explosivos y voladuras. Prensa CRC; 2003: 291–300.https://doi.org/
ijrmms.2013.11.007.
10.1201/9781439833476.ch36.
19 Heelan PA. Radiación de una fuente cilíndrica de longitud finita.Geofísica. 1953; 18
5 Iverson S, Kerkering C, Hustrulid W. Aplicación del enfoque de Holmbergpersson modificado por
(3): 685–696.
NIOSH al diseño de voladuras perimetrales. Actas de la 34a Conferencia sobre explosivos y técnica de
20 Heelan PA. Sobre la teoría de las ondas de la cabeza.Geofísica. 1953; 18 (4): 871–893.
voladura. 2008; 2: 1–33. Blair DP.Algunos problemas asociados con las leyes de escalamiento de
https: // doi.org/10.1190/1.1437941.
6 vibraciones del peso de carga estándar. En: Int. Symp. sobre fragmentación de rocas por voladura.
21 JE blanco, Sengbush RL. Cortar ondas de fuentes explosivas.Geofísica. 1963; 28 (6):
Brisbane. 1990: 149-158.https: // doi. org / 10.1016 / 0148-9062 (92) 91889-d.
1001–1019.https://doi.org/10.1190/1.1439296.
22 Blanco J. Sonido subterráneo: aplicaciones de ondas sísmicas. Londres: Elsevier Science;
7 Blair DP, Minchinton A. Sobre la zona de daños que rodea un solo barreno.
1983.
Fragblast. 1997; 1 (1): 59–72.https://doi.org/10.1080/13855149709408390. Yilmaz O, Unlu T.
23 Triviño LF, Mohanty B, Milkereit B. Formas de onda sísmicas de fuentes explosivas ubicadas
8 Una aplicación del enfoque modificado de Holmberg-Persson para el diseño de voladuras
en pozos e iniciadas en diferentes direcciones. J Appl Geophys. 2012; 87: 81–93.https://
de túneles.Tunn Undergr Space Technol. 2014; 43: 113–122.https: // doi. org / 10.1016 /
doi.org/10.1016/j.jappgeo.2012.09.004. Blair DP. La influencia de la superficie libre sobre la
j.tust.2014.04.009.
24 vibración de la explosión.Ciencia ficción Int J Rock Mech Min. 2015; 77: 182-191.https://
9 Silva J, Worsey T, Lusk B. Evaluación práctica del daño de la roca debido a voladuras. Int J Min Sci
doi.org/10.1016/j.ijrmms.2015.04.006.
Technol. 2019; 29 (3): 379–385.https://doi.org/10.1016/j.ijmst.2018.11.003. Gou Y, Shi X, Zhou J, Qiu X,
25 Couceiro P, Pascual J, Lopez Cano M. Distribución de energía de campo cercano en voladuras de
10 Chen X, Huo X. Evaluación de la atenuación de las vibraciones inducidas por explosiones derivadas de
control de muros: un enfoque de ingeniería. En: 12º Simposio Internacional sobre Fragmentación de
una mina subterránea.Ciencia ficción Int J Rock Mech Min. 2020; 127:
Roca por Voladura.Luleå. 2018: 113–124.
104220. https://doi.org/10.1016/j.ijrmms.2020.104220.
26 Jordan D. La onda de tensión de una fuente explosiva cilíndrica finita. J Math Mech.
11 Salum AH, Murthy VMSR. Optimización de la extracción de explosiones y control de la zona de daño de
1962; 11 (4): 503–551.http://www.jstor.org/stable/24900912. Abo-Zena AM.
excavación inducida por explosiones en la excavación de túneles a través de diversas clases de rocas.Tunn
27 Radiación de una fuente explosiva cilíndrica finita.Geofísica. 1977:
Undergr Space Technol. 2019; 85: 307–318.https://doi.org/10.1016/j.tust.2018.11.029.
42. https://doi.org/10.1190/1.1440799.
12 Meredith JA. Modelado numérico y analítico de fuentes sísmicas de fondo de pozo: el campo cercano
28 Sharpe JA. La producción de ondas elásticas por presiones de explosión.I. Teoría y
y lejano. Tesis doctoral. 1990. Instituto de tecnología de Massachussetts.
observaciones empíricas de campo. Geofísica. 1942; 7 (2): 144-154.https://doi.org/10.1190/
13 Tubman KM. Registros acústicos de forma de onda completa en perforaciones en capas radiales. Tesis doctoral.
1.1445002.
Instituto de Tecnología de Massachussetts; 1984.

17
S. Gómez y col. Revista Internacional de Mecánica de Rocas y Ciencias de la Minería 134 (2020) 104357

29 Rossmanith HP, Uenishi K, Kouzniak N. Propagación de ondas expansivas en macizos 50 Blair DP. Una comparación directa entre la resonancia vibratoria y los datos de transmisión
rocosos — Parte I: medio monolítico. Fragblast. 1997; 1 (3): 317–359.https://doi.org/ de pulsos para evaluar la atenuación sísmica en la roca.Geofísica. 1990; 55 (1): 51–60.
10.1080/ 13855149709408402. https://doi.org/10.1190/1.1442771.
30 Blair DP. Modelos estadísticos de vibración del suelo y chorro de aire.Fragblast. 1999; 3 (4): 51 Bracewell R. La transformada de Fourier y sus aplicaciones. McGraw-Hill;
335–364.https://doi.org/10.1080/13855149909408056. Duvall WI. Formas de ondas de 52 1965. MathWorks Inc. MATLAB R2018b.
31 tensión en rocas cercanas a explosiones.Geofísica. 1953; 18 (2): 310–323.https://doi.org/ 53 Bouchon M. Una revisión del método del número de onda discreto. Geofis de aplicación pura. 2003;
10.1190/1.1437875. 160 (3-4): 445–465.https://doi.org/10.1007/PL00012545.
32 Ricker N. La forma y leyes de propagación de ondas sísmicas. Geofísica. 1953; 18 54 Castedo R, Natale M, López LM, et al. Estimación de los parámetros de Jones-Wilkins-Lee de
(1): 10–40. https://doi.org/10.1190/1.1437843. explosivos en emulsión mediante pruebas de cilindros y su validación numérica.Ciencia ficción Int J
33 Uenishi K, Rossmanith HP. Propagación de ondas expansivas en macizos rocosos. Parte II: Rock Mech Min. 2018; 112: 290–301.https://doi.org/10.1016/j.ijrmms.2018.10.027. Marcha.
medios estratificados.Fragblast. 1998; 2 (1): 39–77.https://doi.org/
34 10.1080/13855149809408879. Vanbrabant F, Chacón EP, Quiñones LA P, Mach S. Ondas 55 Gómez S, Bermejo M, Navarro J, et al. Modelo sísmico de campo cercano para daños en
generadas por la detonación de una carga explosiva cilíndrica - experimentos y rocas. Madrid: SLIM Deliverable Report D3.1, Programa de Investigación e Innovación
simulaciones.Fragblast. 2002; 6 (1): 21–35.https://doi.org/10.1076/frag.6.1.21.8849. Horizonte 2020 de la Unión Europea; 2019. Gran Acuerdo No. 730294.
35 Cho SH, Miyake H, Kimura T, Kaneko K. Efecto de la forma de onda de la presión aplicada sobre el 56 Gibbons JD, Chakraborti S. Inferencia estadística no paramétrica. En: Lovric M, ed.
proceso de fractura de roca en un modelo de cara libre. Sci Technol Energetic Mater. 2003; 64 Enciclopedia internacional de ciencia estadística. Berlín, Heidelberg: Springer Berlin
(3): 116–125. Heidelberg; 2011: 977–979.https://doi.org/10.1007/978-3-642-04898-2_420. Pastor JE,
36 Blair DP. Una solución rápida y eficiente para la radiación de las olas de un barreno presurizado. 57 Teodorczyk A, Knystautas R, Lee JH. Ondas de choque producidas por detonaciones
Fragblast. 2003; 7 (4): 205–230.https://doi.org/10.1076/frag.7.4.205.23533. Rossmanith HP, reflejadas.Prog Astronaut Aeronaut. 1991; 134: 244–264.https://doi.org/
37 Kouzniak N. Detonación supersónica en macizo rocoso - Parte 2: desplazamientos de 10.2514 / 5.9781600866074.0244.0264.
partículas y campos de velocidad para pozos de detonación únicos y múltiples no 58 Drucker DC, Prager W. Mecánica de suelos y análisis de plásticos o diseño de límites. Q Matemáticas de
retardados y retardados.Fragblast. 2004; 8 (2): 95-117.https://doi.org/10.1080/ aplicación. 1952; 10 (2): 157–165.
13855140412331336179. 59 Christensen RM. Una teoría completa de la fluencia y el fracaso de los materiales
38 Cho SH, Kaneko K. Control de fragmentación de rocas en voladuras. Mater Trans. 2004; 45 isotrópicos.J Eng Mater Technol Trans ASME. 2007; 129 (2): 173–181.https://doi.org/
(5): 1722-1730.https://doi.org/10.2320/matertrans.45.1722. 10.1115 / 1.2712847.
39 Triviño L, Mohanty B, Munjiza A. Patrones de radiación sísmica de cargas explosivas 60 Blair DP, Minchinton A. Modelos de vibración por voladura de campo cercano. Proc 8th Int Symp
cilíndricas por métodos combinados analíticos y combinados de elementos finitos- Rock Fragm por Blasting-Fragblast. 2006; 8: 152–159.
discretos. Proc 9th Int Symp Rock Fragm por voladura. 2009: 415–426. 61 Abramowitz M, Stegun IA. Manual de funciones matemáticas con fórmulas, gráficos y tablas
40 Blair DP. Voladura de control de pared.11th Int Symp Rock Fragm de Blasting. 2015: 13–26. Ahn JK, matemáticas.Departamento de Comercio de Estados Unidos. Serie 55 de la Oficina Nacional de
41 Park D. Predicción de la atenuación de la onda de campo cercano debido a una fuente de explosión Estándares de Matemáticas Aplicadas. 1965: 555–566.
esférica.Rock Mech Rock Eng. 2017; 50 (11): 3085–3099.https://doi.org/10.1007/ s00603-017-1274-3. 62 Onederra I, Esen S. Un enfoque alternativo para determinar las constantes de Holmberg-Persson
para modelar la atenuación de la velocidad máxima de las partículas de campo cercano. Fragblast.
42 Yi C, Johansson D, Greberg J. Efectos de las tensiones in situ sobre la fracturación de la roca por 2004; 8 (2): 61–84.https://doi.org/10.1080/13855140412331336151.
voladura. Computación Geotech. Mayo de 2017; 2018 (104): 321–330.https://doi.org/10.1016/j. 63 Meyers MA. Comportamiento dinámico de materiales.John Wiley & Sons, Inc.. 1994: 179–201.
compgeo.2017.12.004. 64 Zel'dovich YB, Raizer YP. Física de ondas de choque y fenómenos hidrodinámicos de alta
43 Minchinton A. Sobre la influencia de los detónicos fundamentales en la práctica de voladuras. 11th temperatura. Publicaciones de Dover; 2002.
Int Symp Rock Fragm de Blasting. 2015: 41–54. McQueen RG, Marsh SP, Fritz JN. Ecuación de Hugoniot de estado de doce rocas.
sesenta y cinco

44 Meza-Fajardo KC, Lai CG. Relaciones causales explícitas entre la relación de amortiguamiento del J Geophys Res. 1967; 72 (20): 4999–5036.https://doi.org/10.1029/
material y la velocidad de fase a partir de soluciones exactas de las ecuaciones de dispersión de la JZ072i020p04999.
viscoelasticidad lineal.Geophys J Int. 2007; 171 (3): 1247–1257.https://doi.org/10.1111/ 66 Anderson GD, Larson DB. Estudios de ondas de choque planas de medios geológicos.https: // inis.
j.1365-246X.2007.03590.x. iaea.org/search/search.aspx?orig_q=RN:9394856; 1977. Consultado el 19 de diciembre de
45 Zhu T, Carcione JM, Harris JM. Aproximación de la propagación sísmica de Q constante en 2019.
el dominio del tiempo.Perspectiva de Geophys. 2013; 61 (5): 931–940.https://doi.org/ 67 Marsh SP. Datos de LASL Shock Hugoniot. Berkeley, CA. 1980.
10.1111/ 1365-2478.12044. 68 Sekine T, Duffy TS, Rubin AM, Anderson WW, Ahrens TJ. Choque de compresión y liberación
46 Minster JB, Day SM. Desintegración de los campos de ondas cerca de una fuente explosiva debido a isentrópica de granito.Geophys J Int. 1995; 120 (2): 247–261.https://doi.org/
una atenuación no lineal de alta deformación.J Geophys Res. 1986; 91: 2113–2122.https://doi.org/ 10.1111 / j.1365-246X.1995.tb01817.x.
10.1029/ JB091iB02p02113. 69 Sanchidrián JA, López LM. Cálculo de la energía de explosivos con un modelo de reacción
47 Kjartansson E. Propagación y atenuación constantes de la onda Q. Tierra sólida de J parcial. Comparación con datos de prueba de cilindros.Propelentes, Explos Pyrotech.
Geophys Res. 1979; 84 (B9): 4737–4748.https://doi.org/10.1029/JB084iB09p04737. 2006; 31 (1): 25–32.https://doi.org/10.1002/prep.200600003.
48 Kjartansson E. Atenuación de ondas sísmicas en rocas y aplicaciones en exploración 70 Fickett W, WC Davis. Detonación. En: Berkeley, ed.Los Ángeles - Londres. Prensa de la Universidad de
energética. Tesis doctoral. Universidad Stanford; 1980. California; 1979.
49 Blair DP, Spathis AT. Atenuación del pulso generado por explosión en macizos rocosos.
Tierra sólida de J Geophys Res. 1982; 87 (B5): 3885–3892.https://doi.org/10.1029/
JB087iB05p03885.

18