SUBSIDENCIA EN LA MINA EL TENIENTE

 HARO JARA, LUIS-FERNANDO JOEL

 HUAYTALLA TENORIO, LUIS MIGUEL
 TORRES ROJAS, DARWIN
 TORRES ORTIZ, PAVEL
 VALDEZ CHILLITUPA, YURY
 INDICE

 Introducción al Caving
 Fenómeno de Subsidencia
 Subsidencia en Minería
 Causas y efectos
 Monitoreo y Control/Instrumentación Geomecánica
 Modelamiento
 Simple/Metodologías de trabajo
 Análisis Predictivo
 Conclusiones y Recomendaciones
 SUBSIDENCIA ASOCIADA A MINERÍA POR HUNDIMIENTO

La minería por métodos de Efecto de la posición relativa entre un nuevo

hundimiento consiste en inducir el sector productivo y los sectores superiores
quiebre y hundimiento del macizo ya explotados
rocoso, lo que permite la extracción del
mineral quebrado y al mismo tiempo Los nuevos sectores en explotación se
genera una cavidad por encima del piso ubican siempre dentro de los límites de
del hundimiento, la cual termina por un sector ya explotado.
conectarse a superficie, generando un
cráter de subsidencia.

Hundimiento en subniveles
Los nuevos sectores en explotación se
• Conexión a superficie en corto plazo
extienden más allá de los límites del sector
ya explotado.
Hundimiento de bloques
Se pierde:
• Conexión a superficie en mediado y
largo plazo
• Efecto sombra
• Efecto escalón
 EFECTOS EN LA FORMACION DE CRÁTER EN SUBSIDENCIA

EFECTO DE LA TOPOGRAFÍA EFECTO DE LA GEOLOGÍA

Paredes empinadas:
Sectores de Sectores de • Roca primaria > Roca secundaria > Material
mayor altura menor altura
tipo talus o de cobertura
Pared del cráter Tendida Empinada • Cuerpos litológicos masivos > Macizos
Inestabilidad en las Mayor volumen Menor volumen
rocosos fracturados
paredes de material de material • Sistemas estructurales y/o estructuras mayores
de manteo empinado a subvertical: Rumbo
Escorrentías de Hacia el cráter Hacia otro lado
superficie similar a la estructura
 MORFOLOGÍA DE UN CRÁTER DE SUBSIDENCIA

Tipo circular Tipo herradura Tipo rosquilla

• Angulo de desplome (α)

• Angulo de influencia (β)
• Zona de influencia
• Ancho de la zona de influencia
• Ancho basal del cráter
• Ancho perimetral del crater
• Altura de la pared del cráter (H)
• Columna de material quebrado (h)
• Distancia a la zona agrietada
• Ancho de la zona agrietada
 MECÁNICA DE LA FORMACIÓN DEL CRÁTER DE
SUBSIDENCIA

1. Se desarrolla la socavación del 2. El material quebrado comienza (a) “Volumen activo” en condición sin (b) “Volumen activo” en condición con
UCL y se genera un volumen de a rellenar la cavidad. La avance de la socavación. avance de la socavación.
macizo “activo” que tiende a extracción provoca la propagación
desplazarse verticalmente hacia hacia arriba del volumen de
la cavidad. macizo “activo”. 3. El proceso continúa. Desplazamientos horizontales en la parte
inferior de la pared de la cavidad.
 MECÁNICA DE LA FORMACIÓN DEL CRÁTER DE
SUBSIDENCIA

4. La cavidad tiende a tomar una

forma de “catedral” y puede : (1)
estabilizarse y no continuar
creciendo; o (2) continuar
creciendo más rápidamente y
definir un crown-pillar en
superficie.
 MECÁNICA DE LA FORMACIÓN DEL CRÁTER DE
SUBSIDENCIA

5. Se produce la ruptura del crown- 6. La parte superior de las 7. Continua el crecimiento del 8. Finalmente se alcanza una
pillar; formándose un cráter tipo paredes se derrumba y se perímetro del cráter y la condición de equilibrio donde el
chimenea, de paredes subverticales. ensancha el perímetro del cráter disminución de la inclinación de cráter no sigue creciendo. Esto
El macizo tiende a desplazarse en superficie, disminuye la sus paredes, pero la condición es puede ocurrir varios años
horizontalmente hacia la cavidad. inclinación de sus paredes. cada vez más estable. después del término de la
extracción de mineral.
 EVALUACIÓN DE LA SUBSIDENCIA ASOCIADA A MINERIA POR
MÉTODOS DE HUNDIMIENTO

•Definir en forma empírica valores: para los ángulos de

Métodos desplome, α, y de influencia, β;
empíricos •Correlacionar el volumen de la subsidencia con el
volumen de mineral extraído.

Método de
análisis •Predecir crecimiento de un cráter de subsidencia en la
equilibrio límite medida que aumenta la profundidad de la minería
(Hoek, 1974)

•Incluir el efecto de una

Brown & superficie de terreno inclinada
Ferguson (1979) y, también, la posible presencia
de aguas subterráneas.
 ANÁLISIS DE SUBSIDENCIA

Fenómeno de inestabilidad donde los efectos son locales y tienen lugar en niveles superiores. Para una
desarrollar un amplio conocimiento es necesario conocer:

• Ángulos que acontecen en un cráter. Son dos: Ángulo de desplome (Alpha) y

Ángulo de subsidencia ángulo de fracturamiento (Beta)

• Área delimitada de un sector productivo, donde las reservas deben ser

Pilar de protección excluidas para proteger la zona de los límites de subsidencia.

Monitoreos • Control del movimiento y dirección de la subsidencia

 MONITOREOS EN SUBSIDENCIA

Para determinar el movimiento de la subsidencia, es necesario recurrir a técnicas para conocer el comportamiento de
movimiento de la subsidencia a través de diferentes metodologías.

Inspecciones Visuales
• Usada durante años, determina los efectos de subsidencia.

Instrumentación
• Uso de herramientas técnicas, que permiten recopilar una información más precisa.
 TIPOS DE MONITOREOS

Monitoreos Visuales Monitoreos Instrumentales

Cable TDR

Medidor de Desplazamiento
Mapeo de daños en labores

Inspección de fisiografía del cráter Indicador de falla de corte

Levantamiento topográfico Extensómetro

Uso de perforaciones con cámara Inclinómetro

Sísmica
 Metodología de trabajo

Diferentes métodos se han desarrollado durante la vida de Mina de El Teniente, entre ellos tenemos tres grandes
períodos:

Período A Período B Período C

Explotación por BC, se Explotación por PC, ángulos de Explotación por PC, los ángulos
obtuvieron ángulos de subsidencia obtenidos por análisis de subsidencia se obtienen de
subsidencia empíricamente. de inspecciones visuales e herramienta computacional.
instrumentales (cables TDR)
 METODOLOGÍA DE TRABAJO

ETAPA DEL
PARÁMETROS METODOLOGÍA B METODOLOGÍA C
PROCEDIMIENTO
Áreas hundidas, fisiografía del cráter y
Geometría del sector Topografía
ubicación de los daños.
Geología estructural
Geología Hidrogeología Litología
Calidad Geotécnica
Monitoreos
Cables TDR
Instrumentación Sísmico
Medidores de desplazamiento
Propiedades del MR
Geotecnia/Geomecánica No se toma en cuenta
Modos de falla
Otros No se toma en cuenta Experiencia en otras minas

Ángulos de subsidencia Ángulo de desplome& fracturamiento Para todo el sector sin considerar efectos Para todo el sector y se considera efectos

Efectos en galerías Infraestructura en los cráteres

Eefectos de subsidencia
en niveles superiores Inestabilidad de cuña
Efectos en superficie
Hoya hidrográfica
Buscar alternativas de reemplazo de infraestructura
Acciones a tomar Determinación de pilares de protección
Definir un plan de monitoreo
 ANÁLISIS DE DESMONTE EN UN CRÁTER
ACTIVO

El análisis indirecto de estabilidad es manejado donde los ratios de subsidencia en la zona donde se encuentran y se ubican los
botaderos de desmonte.

Características

 Desintegración progresiva
 Representa impactos en la
dilatación y volumen
 Se ha identificado una tasa de 40
m/año

Modelamiento numérico

Algoritmo de hundimiento: Hoek

Representa: suavizado, ajuste de densidad, dilatación,
propiedades
TIPOS DE FRACTURAMIENTO
INDUCIDO
Fracturamiento leve o moderado
Ocurre en el limite entre la roca no perturbada y
aquella con los primeros indicios de daño por
subsidencia.
Caracterizado por pequeñas fracturas en la
matriz de la roca asociado a niveles de presión
bajos.
Fracturamiento severo
Localizada entre la roca totalmente fractura y
la roca perturbada por el minado.
Esta caracterizado por fracturas de tamaño
considerable y la apertura de estructuras
relevantes, formando bloques inestable que
caerán en el cráter de subsidencia.
 MODELAMIENTO SIMPLE DEL
FRACTURAMIENTO INDUCIDO

Características – Modelo Simple Geometría del sector Este de la Pipa Braden

• Aborda una única unidad geotécnica.
• Geometría simplificada para la pipa y el área de
minado.
• Empleado, principalmente, para realizar análisis de
sensibilidad de los parámetros e interpretar estos
resultados.
Altura del sector
Altura total del hundido = 600 m
modelo = 1000 m
Resultados – Modelo Simple
• Sirve como criterio para la posterior construcción
de los modelos detallados.
• Base para el análisis de los resultados e
identificación de las variables clave que influyen
en el fracturamiento.
 MODELAMIENTO SIMPLE DEL
FRACTURAMIENTO INDUCIDO

Contornos de deformación “total”

• El análisis de cavidades totalmente
hundidas tuvieron mejor resultado que
aquellas con llenas de material roto.
• El único parámetro con mucha
importancia para el fracturamiento de la
= brecha fue la fuerza de tensión de la
masa rocosa.
Las concentraciones • Se decidió que el contorno de 5e-3 de
de estos contornos son esfuerzo total es similar al limite de
buenos indicadores fracturamiento identificado en el
del ángulo de falla de
una masa rocosa.
campo.
 MODELAMIENTO DETALLADO DEL
FRACTURAMIENTO INDUCIDO

Características – Modelo Detallado Forma de las cavidades por el minado 2005 - 2025
• Aborda una cantidad mayor de unidades
geotécnicas. 2005 2025
• Posee una geometría considerablemente detallado
para la Pipa y el área de minado.

Resultados – Modelo Detallado

• Identificación de la subsidencia actual y asociación
de la misma a los niveles de esfuerzo.
• Simulación del futuro minado y predicción del
comportamiento de la subsidencia hasta el año
2025.

Cavidades huecas por

Pipa Braden el minado
 MODELAMIENTO DETALLADO DEL
FRACTURAMIENTO INDUCIDO
Comparación de la subsidencia real vs
modelamiento numérico en 3DEC La profundidad del daño en la Pipa Braden
depende muchos factores, entre ellos los
siguientes:
Limite de • Fuerza de la masa rocosa.
Fracturamiento • Profundidad del minado.
3DEC • Efecto de estructuras de control.
• Efecto de la inclinación de la pared de la
Este limite coincide
para la roca que se Pipa.
encuentra en
condiciones promedio La estabilidad de la Pipa y el impacto que podría
o “normales”. tener en las instalaciones que contiene, son temas
de largo plazo y de vital importancia para la mina.
Por esta razón se busca predecir su comportamiento
a lo largo de la vida de la mina; pero, debido a los
factores anterior mencionados, no se puede
establecer un modelo que cumpla fielmente lo
ocurrido en el campo.
 MODELAMIENTO BIDIMENSIONAL DETET-1
DESCRIPCIÓN DEL MODELO
HIPÓTESIS
DETET-1
Es un modelo bidimensional de equilibrio límite, el cual
considera que la estabilidad de la pared del cráter (ángulo
Condición de de inclinación α), queda definida por un plano de
equilibrio límite deslizamiento (ángulo de inclinación β), afectado
por una grieta de tracción.
Inestabilidad por
falla planar
OBJETIVOS
El problema es
estático  Determinar gráficamente los ángulos de
desplome α
No existen fuerzas  Determinar gráficamente el ángulo de
externas fracturamiento β.
 Los ángulos α y β en función del RMR y la
profundidad donde se encuentra el nivel de
Se considera las Material Hundido hundimiento.
características del
material Macizo Rocoso
MODELAMIENTO BIDIMENSIONAL DETET-1
PARÁMETROS DE ENTRADA Estimación de resistencia al corte según el RMR_
Bieniawski
VALOR DEL CARACTERIZACIÓN
COHESIÓN ÁNGULO DE
INDICE CUALITATIVA DEL

Parámetros RMR
0-20
MACISO ROCOSO
MUY MALO
(kPa)

<100
FRICCIÓN

<15

Geométricos 21-40
41-60
MALO
REGULAR
100-200
200-300
15-25
25-35
61-80 BUENO 300-400 35-45

Propiedades del 81-100 MUY BUENO >400 >45

Macizo Rocoso
Propiedades del
material hundido
Parámetros
Geomecánicos
 MODELAMIENTO BIDIMENSIONAL DETET-1
PARÁMETROS DE
SALIDA DEL MODELO

Observación
Cada curva se diseña según un ángulo
de terreno , evaluado en condición seca
y con agua.

Por ejemplo: Una grieta de tracción

llena con agua y drenaje al pie de talud
Relación entre el ángulo de desplome y
Curvas de diseño para estimar el ángulo de desplome fracturamiento
 MODELAMIENTO BIDIMENSIONAL DETET-1
CALIBRACIÓN DE LOS RESULTADOS DEL
MODELO

Es la verificación y/o comparación de los parámetros de salida del

modelo con la información recolectada en la etapa de Recopilación
de antecedentes.

Con el fin de calibrar o ajustar las curvas

teóricas y obtener las curvas de diseño
operacional más ajustadas al diseño real.

CURVAS DE DISEÑO OPERACIONAL PARA

ESTIMAR EL ÁNGULO DE DESPLOME

75
 SUBSIDENCIA EN LA MINA EL TENIENTE

Objetivo
Predecir el ángulo de desplome y la
extensión de la zona de influencia DEFINICIÓN DE LAS SECCIONES DE INTERÉS​

Procedimiento

• Se revisaron antecedentes 
• Se definieron secciones de interés
• Caracterización geotécnica y análisis
numérico
• Contrastar resultados
• Definición de criterios para definir los
objetivos.
• Definición de curvas de predicción 
• Desarrollo de una metodología para el
uso de estas curvas
 SUBSIDENCIA EN LA MINA EL TENIENTE

CARACTERIZACION GEOMECÁNICA Y ANÁLISIS NUMÉRICO

ESTADO TENSIONAL IN SITU

Procedimiento Se supuso que el estado tensional in situ queda definido por los siguientes esfuerzos
principales: 
• Un esfuerzo vertical de tipo gravitacional 
• Un esfuerzo horizontal en el plano de la sección analizada igual a 1.4 veces el
esfuerzo vertical.  
• Un esfuerzo horizontal en el plano normal a la sección analizada igual a 1.2 veces
el esfuerzo vertical.
PROPIEDADES GEOMECÁNICAS DE LAS ESTRUCTURAS
MAYORES
• La rigidez normal, k n sea igual a su rigidez tangencial, k t , y ambas tienen una
magnitud de 1011 Pa/m.
• Ángulo de fricción, φ j, era de 30º.
• cohesión, c j ,era de 150 kPa.
• resistencia en tracción sea suficientemente alta como para que el activamiento de
las estructuras se deba exclusivamente a una falla por corte (con el fin de evitar
problemas numéricos por intersección de varias estructuras)
 SUBSIDENCIA EN LA MINA EL TENIENTE

CARACTERIZACION GEOMECÁNICA Y ANÁLISIS NUMÉRICO

Procedimiento

“El mejor parámetro para evaluar el efecto del cráter era el desplazamiento en el
punto de interés”
 SUBSIDENCIA EN LA MINA EL TENIENTE

EVALUACIÓN DEL ANGULO DE DESPLOME

Procedimiento ALCANCES:

UCL: Piso del nivel de socavación

FW: Costado opuesto o alejado de la Pipa Braden,

usualmente correspondiente a la zona de montaña.

HW: Costado adyacente o vecino de la Pipa

Braden 
 SUBSIDENCIA EN LA MINA EL TENIENTE

EVALUACIÓN DE LA ZONA DE INFLUNCIA

Procedimiento

t i (z) : ancho de la zona de influencia, medido desde la pared del cráter, a una altura z sobre el piso del UCL.
d(z) : valor relativo del ancho de la zona de influencia a
la altura z, el cual se obtiene de la curva correspondiente al sector en evaluación;
C: coeficiente de proporcionalidad que depende del sector que se analice
 SUBSIDENCIA EN LA MINA EL TENIENTE

METODOLOGÍA PROPUESTA

Procedimiento
• Definir el sector donde se interesa evaluar la subsidencia
• Desarrollar una sección vertical que muestre la geología geotécnica del sector que interesa, indicando la
presencia de los siguientes tipos de macizo rocoso : roca primaria, roca secundaria, y Pipa Braden; así como la
profundidad del UCL.
• Seleccionar la sección analizada en este estudio que esté más próxima a la sección que interesa evaluar.
• Seleccionar entre la curva correspondiente a dicha sección para evaluar el ángulo de desplome 
• Seleccionar entre la curva correspondiente a dicha sección para evaluar el ancho de la zona de influencia
• Definir una faja de seguridad de 40 m respecto al límite de la zona de influencia así definida.
 CONCLUSIONES

• Para el caso de la mina el Teniente , el parámetro más influyente para la evaluar el efecto de la subsidencia , es
midiendo el desplazamiento a partir de un punto de interés.

• En los resultados para diferentes alternativas de “profundidad de minado”, en este caso entre el nivel Teniente 7 y
el nivel Teniente 8, con una diferencia aproximada de 70 metros, no muestran una diferencia significativa en el
efecto de subsidencia.

• Emplear monitoreos por instrumentación y visuales, generan una mejor recepción de información, por lo que es
clave trabajar ambos monitoreos.