Tesis Compactación Tranques de Relave

PONTIFICIA UNIVERSIDAD CATÓLICA DE VALPARAÍSO
FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA EN CONSTRUCCIÓN
METODOLOGÍA PARA LA MODELACIÓN DEL CONTROL DE COMPACTACIÓN EN

TRANQUES DE RELAVES.
Tesis para optar al título de Ingeniero

Constructor.
Autor.
Ojeda García, Matías Humberto.
Zamora González, María de los Ángeles.
Profesor Guía.
Gabriel Villavicencio Arancibia.
Octubre 2015
DEDICATORIA.
La finalización de esta grandiosa etapa va dedicado a mis padres Ledy y René por haberme
apoyado en todo momento de mi vida, por sus valores y consejos y por sobre todo hacerme una
persona de bien, a mis hermanas Mary por quererme como si fuese su hijo, desde el primer
inicio de mis días y Yarela, por ser fiel compañera en mi infancia, a mis tíos Tatiana e Iván,
pilares fundamentales e irremplazables en mi desempeño como persona y estudiante en los
primeros años de universidad y estadía en Valparaíso; Camilo, por ser mi hermano cual sean
las circunstancias, Marcos mi fiel compañero de la Universidad y vida, Eunice mi fiel
confidente, mi compañera de tesis María de los Ángeles, fiel compañera de amistad y estudio
y por último, a mi Profesor y amigo Gabriel Villavicencio Arancibia, por ser el impulsor al
inicio, desarrollo y finalización de esta tesis y por haberme guiado al grandioso camino de la
Mecánica de Suelos.
Gracias totales a todos.
Matías Ojeda García.
Tras finalizar esta etapa, me es necesario agradecer a todos los involucrados en el desarrollo
de ésta. Con especial importancia agradezco a mi madre, por todo el apoyo y sacrificio
entregado, por siempre creer en mí y ser un pilar fundamental de amor y motivación, y de
esta forma poder concluir con éxito. A mi padre, que a pesar de no estar físicamente, confío en
que me acompañó en cada uno de los pasos que di.
Agradezco la oportunidad de haber estudiado en esta institución, y específicamente en esta

escuela, en donde conocí y compartí con grandes personas, de las cuales guardaré los mejores
recuerdos.
Por último, agradezco a nuestro profesor Gabriel por brindarnos su tiempo y ayuda en el
transcurso de esta tesis y a la colaboración por parte de la empresa Geotecnia Ambiental,
entregándonos información esencial para el desarrollo de esta tesis.
Mª Ángeles Zamora G.
RESUMEN.
RESUMEN.
El control del proceso de compactación de los tranques de relaves actualmente es puntual,

generando muchas veces una deficiente representación de calidad del proceso de construcción.
Lo anterior provoca que no se evidencien las zonas de baja resistencia dentro de la estructura
del tranque, que puedan generar posibles zonas de fallas.
A partir de lo anterior surge la necesidad de crear una nueva manera de control de compactación,
que se sustenten la hipótesis de que es posible reconstruir la estructura interna del tranque de
relave, mediante la modelación de los parámetros del control de compactación.
En base a lo anterior expuesto, en esta tesis se genera una metodología para la modelación del
control de compactación, la cual se logra mediante la siguiente comprensión de ejes.
– Revisión bibliográfica de los temas acorde al marco teórico y estado del arte de tranques de
arenas de relaves en Chile.
– Análisis de la variabilidad del parámetro resistencia por punta (qd), obtenido del control de
compactación mediante el penetrómetro ligero (PANDA) de una zona de estudio. En donde se
aprecia mediante la utilización de herramientas clásicas de la estadística, el comportamiento
homogéneo del parámetro de resistencia por punta, utilizándolo posteriormente como
herramienta de rectificación de la modelación del control de compactación.
– Estimación de características cualitativas en función de la densidad relativa (DR),

determinando compacidad, comportamiento mecánico y riesgo de licuefacción, evidenciando
zonas de posibles fallas, con compacidad compacta, potencial de licuación de equilibrio y
comportamiento mecánico contractivo, específicamente en sectores altos de la zona de
estudio.
– Modelación de los parámetros del control de compactación, parámetro de resistencia por punta
(qd), apreciando una variación de resistencia por punta de 3 [MPa] a 6 [MPa], siendo valores
aceptables para los controles de compactación que se requieren la norma chilena.
Posteriormente se realiza una modelación de características cualitativas como la compacidad,
entregando en casi la mayoría de las zonas tonalidades azules que indican compacidades
densas a compactas, rectificando ambas modelaciones con la variabilidad de los parámetros
del control de compactación estudiados en el capítulo3.
– Conclusiones que corroboran la hipótesis de trabajo de reconstruir la estructura interna del

tranque de relave, mediante los parámetros del control de compactación, ya sea cuantitativos
resistencia por punta y cualitativos (compacidad, comportamiento mecánico), apreciando
zonas débiles o corroborando calidad del proceso del control de compactación.
ÍNDICE.
RESUMEN.
ÍNDICE.
INTRODUCCIÓN.
1. SÍNTESIS BIBLIOGRÁFICA. 1
1.1. DEPÓSITOS DE RELAVES EN CHILE. 1
1.2. PRINCIPALES MECANISMOS DE FALLAS EN TRANQUES DE RELAVE. 4
1.2.1. Licuefacción o falla de flujo. 4
1.2.2. Inestabilidad del suelo de fundación. 8
1.2.3. Inestabilidad de taludes. 8
1.2.4. Deformaciones. 9
1.2.5. Tubificación o piping. 10
1.2.6. Overtopping. 11
1.2.7. Aluviones. 12
1.3. FALLAS HISTORICAS EN TRANQUES DE RELAVES. 13
1.3.1. Tranque Barahona. 13
1.3.2. Tranque El cobre. 14
1.3.3. Tranque Las Palmas y Tranque Adosado. 14
1.4. PROGRESO DE LA NORMATIVA CHILENA EN TRANQUES DE RELAVE. 15
1.4.1. Decreto Supremo N°86. 15
1.4.2. Decreto supremo N°248. 16
1.5. VARIABILIDAD DE LOS PARÁMETROS DEL CONTROL DE COMPACTACIÓN EN

TRANQUES DE RELAVES. 19
1.5.1. Metodología de análisis de la variabilidad. 20
1.5.2. Descripción de los tranques de relaves experimentales. 21
1.5.3. Análisis de la variabilidad de las características físicas y resistencia in situ. 22
1.5.4. Distribución estadística de los parámetros físicos y mecánicos. 23
1.6. COMPACTACIÓN EN TRANQUES DE RELAVES. 26
1.6.1. Control de compactación en tranques de arenas de relave mediante el

penetrómetro ligero PANDA. 29
1.7. HIPÓTESIS Y METODOS DE MODELACIÓN. 32
1.7.1. Métodos de interpolación espacial software ROCKWORK. 33
1.7.2. Herramientas de visualización del modelo espacial. 34
2. ESTIMACIÓN DE LA VARIABILIDAD DEL PARÁMETRO DE RESISTENCIA A LA

PENETRACIÓN (qd) DEL TRANQUE EN ESTUDIO. 37
2.1. INFORMACIÓN GENERAL DEL TRANQUE DE RELAVE EN ESTUDIO (TRANQUE M-

1). 37
2.2. RECOPILACIÓN DE INFORMACIÓN. 38
2.2.1. Descripción del área de estudio. 39
2.2.2. Emplazamiento de prospecciones. 39
2.3. ANÁLISIS ESTADÍSTICO DE LOS DATOS RECOPILADOS. 40
2.3.1. Resumen estadístico. 40
2.3.2. Distribución estadística del parámetro de resistencia a la penetración (qd). 41
2.3.3. Leyes de probabilidad que se ajustan a la distribución del parámetro de

resistencia a la penetración (qd). 42
2.3.4. Conclusiones de variabilidad Tranque M-1. 43
3. ESTIMACIÓN DE LA COMPACIDAD, COMPORTAMIENTO MECÁNICO Y RIESGO DE

LICUEFACCIÓN EN FUNCIÓN DEL PARÁMETRO RESISTENCIA A LA PENETRACIÓN
(qd). 45
3.1. VARIACIÓN DE LA RESISTENCIA POR PUNTA (qd). 47
3.2. RELACIÓN ENTRE LA RESISTENCIA POR PUNTA NORMALIZADA (qdN1) Y EL

GRADO DE COMPACTACIÓN. 55
3.2.1. Normalización del parámetro resistencia (qd). 56
3.3. CLASIFICACIÓN Y ANÁLISIS DEL GRADO DE COMPACTACION, COMPORTAMIENTO

MECÁNICO Y RIESGO DE LICUEFACCIÓN. 57
3.3.1. Clasificación y análisis estadístico del parámetro densidad relativa (DR) a nivel
global por perfil. 59
3.3.2. Clasificación del parámetro densidad relativa (DR) a nivel local por prospección.
62
3.3.3. Análisis estadístico del parámetro densidad relativa (DR) a nivel local por
prospección. 66
4. MODELACIÓN DE LOS PARAMETROS DEL CONTROL DE COMPACTACIÓN. 73
4.1. METODOLOGÍA DE MODELACIÓN ESPACIAL. 73
4.2. MODELACIÓN Y ANÁLISIS DE LA ZONA DE ESTUDIO. 75
4.2.1. Modelación parámetro resistencia por punta (qd). 76
4.2.2. Modelación estado de compacidad (parámetro cualitativo). 81
4.2.3. Modelación comportamiento mecánico (parámetro cualitativo). 84
CONCLUSIÓN.
BIBLIOGRAFÍA.
INTRODUCCIÓN.
I
INTRODUCCIÓN.
Chile es uno de los países de gran historia minera, con una producción anual aproximada de más
de cinco millones de toneladas, en lo que representa a su mineral más destacado, el cobre lo que
constituye más de 36% de la producción mundial.La extracción de este mineral a lo largo del
tiempo en Chile ha traído consecuencia de disminución en su ley, concepto que hace referencia
a la cantidad de cobre por volumen de material removido para su explotación, que llega al orden
del 1% de cobre por tonelada, porcentaje que genera múltiples problemas asociados a los lugares
de almacenamiento de dichos residuos sobrantes de material mineral residual.
La problemática antes mencionada ha estado presente a lo largo de la historia minera en Chile,

basado en la imposibilidad de reutilizar con un concreto uso estos materiales. Es así como las
empresas mineras en un principio se deshacían de este material vertiéndolos a los cauces de los
ríos, acciones que debieron modificarse a medida que la legislación ambiental fue creciendo,
naciendo así la necesidad de buscar lugares de almacenamiento.
Específicamente, sobre los depósitos de relaves en Chile, existen cinco maneras de

almacenamiento: espesados, filtrados, en pasta, embalses y tranques, siendo estos dos últimos
los con mayor importancia por la cantidad en operación.
Los tranques de relaves están constituidos por material extraído de las fracciones más gruesas
del relave que forman el muro resistente, siendo este almacenamiento el más utilizado en Chile,
con una cantidad en operación de ciento cincuenta y nueve, y con un porcentaje sobre el total de
depósitos del 74%.La gran utilización de este depósito se debe al bajo costo en comparación a
los otros métodos de almacenamiento, sin embargo el que presenta mayor inconvenientes en
relación a la estabilidad que puedan presentar las presas y los riesgos asociados a ello.
Actualmente el control de compactación de los tranques de relave durante la fase de operación,

es realizado en forma puntual y superficial, en relación al gran volumen depositado que forma
parte del muro del tranque de relaves. Esto se traduce en una deficiente representación del
comportamiento interno del tranque, lo que no asegura la calidad del proceso de construcción.
Debido a esta problemática, en esta tesis se genera una metodología para la modelación del
control de compactación, bajo la hipótesis de reconstruir la estructura interna del tranque de
relave, evidenciando así la variabilidad de las características del tranque en forma espacial, y
distinguiendo las zonas con bajas resistencias que provoquen fallas o colapsos.
Para el desarrollo de esta metodología se ejecutan tres objetivos específicos.
‒ Analizar la variabilidad de los parámetros que definen el control de compactación en el tranque

de relaves en estudio, con el objetivo de conocer la variación que poseen las características
del tranque. Lo anterior se realiza a través del análisis estadístico que adoptan ciertas leyes
II
de distribución teóricas, con la finalidad de demostrar la hipótesis de homogeneidad de los

datos recopilados.
‒ Estimar el grado de compacidad, comportamiento mecánico y riesgo de licuefacción en función

del parámetro de resistencia por punta (qd), obtenido durante el control de compactación.
‒ Modelar la base de datos global para reconstruir la estructura interna del tranque de relaves.
El trabajo de tesis se organiza de la siguiente manera.
‒ Revisión de síntesis bibliográfica.
‒ Estimación de la variabilidad del parámetro de resistencia a la penetración (qd) del tranque en

estudio.
‒ Estimación de la compacidad, comportamiento mecánico y riesgo de licuefacción en función

del parámetro resistencia a la penetración (qd).
‒ Modelación de los parámetros que definen el control de compactación en tranques de relave.
‒ Exposición de conclusiones y recomendaciones.

1. SÍNTESIS BIBLIOGRÁFICA.
1. SÍNTESIS BIBLIOGRÁFICA.
1.1. DEPÓSITOS DE RELAVES EN CHILE.
Los depósitos de residuos mineros en Chile, se pueden clasificar en cinco tipos; relaves en pasta,
espesados, filtrados, embalses y tranques de relaves, siendo los dos últimos los con mayor
incidencia o importancia por la cantidad en operación. Los embalses son estructuras donde su
muro es construido a partir de material de empréstito, mientras el tranque de relave tiene como
material de construcción la fracción más gruesa del mismo relave (arenas).
La última información de la situación de los depósitos de relaves en Chile, entregada por Servicio
Nacional de Geología y Minería (SERNAGEOMIN), fue actualizada el año 2015, con la finalidad
de entregar cantidad y distribución por región. En este catastro se evalúa el escenario de cada
uno de los depósitos existentes para cada región, dando a conocer datos como empresa,
ubicación, el mineral producido, el tipo de depósito, coordenadas, entre otros. Esta información
concluye que actualmente existen seiscientos tres depósitos, de los cuales doscientos dieciséis
se presentan en estado activo, doscientos cuarenta y cuatro no activos y ciento cuarenta y tres
sin información.
En el gráfico de la Figura 1.1. se muestra el porcentaje correspondiente a los depósitos de relaves

según su estado operacional.
DISTRIBUCIÓN DEPÓSITOS SEGÚN ESTADO

OPERACIONAL
24%
36% ACTIVOS
NO ACTIVOS
SIN INFORMACION
40%
Figura 1.1. Distribución depósitos según estado operacional.

Fuente. Elaboración propia.
Escuela de Ingeniería en Construcción. 1

Metodología para la modelación del control de compactación en tranques de relaves.
De los seiscientos seis depósitos catastrados, se pueden clasificar según el tipo de depósitos. De
esta forma se puede apreciar la importancia y presencia de cada uno de éstos. En el gráfico de
la Figura 1.2., se muestra el porcentaje correspondiente a los depósitos de relaves totales según
su tipo de almacenamiento.
DISTRIBUCIÓN POR TIPO DE DEPÓSITO

4% 1%
1% 1%
11%
FILTRADO
ESPESADO
TRANQUE
EMBALSE
EN PASTA
S/I
82%
Figura 1.2. Distribución por tipo de depósitos.

Un 82% son tranques de relaves, 11% los embalses, siendo éstos los con mayor importancia en
número en comparación a los otros depósitos, lo que implica un mayor control y estudio de éstos.
Sin embargo, esta tesis pretende enfocarse sólo en los tranques en ejecución o situación activa,
por lo tanto surge la necesidad de visualizar la predominación de los depósitos activos. En el
gráfico de la Figura 1.3., se muestra el porcentaje correspondiente a los depósitos de relaves en
condición activa.
De los doscientos dieciséis depósitos de almacenamientos activos catastrados el 2015, un 74%

son tranques de relaves, 19% embalses, ratificando la incidencia del tranque de relave a lo largo
de la historia en comparación al total de almacenamientos, ya sea en estado activos o no activos.
Por último se expone una comparación entre el catastro del año 2010 con la última actualización
realizada el año 2015, para de esta forma observar el crecimiento por tipo de depósito y así
determinar cuál es el que se encuentra predominando en estos últimos años.

DISTRIBUCIÓN DEPÓSITOS ACTIVOS, AÑO 2015
2% 0% 3%
2%
19% FILTRADO
ESPESADO
TRANQUE
EMBALSE
EN PASTA
S/I
74%
Figura 1.3. Distribución tipo de depósitos activos año 2015.

En el gráfico de la Figura 1.4. se presenta la variación según tipo de depósito entre catastros.
VARIACIÓN SEGÚN TIPO DE DEPÓSITO ENTRE

LOS AÑOS 2010-2015
180
160
140
120
100
AÑO 2010
80
AÑO 2015
60
40
20
0
FILTRADO ESPESADO TRANQUE EMBALSE EN PASTA S/I
Figura 1.4. Variación según tipo de depósitos entre los años 2010-2015.
Fuente. Elaboración Propia.
Una vez más se corrobora la importancia de la construcción de tranques de relaves, los cuales
doblan en cantidad con respecto al catastro anterior; y los embalses tienen un comportamiento

similar. Cabe decir, que en este nuevo catastro todos los depósitos activos presentan información
con respecto a su tipo de depósito, por lo tanto los depósitos activos sin información del año 2010
se encuentran distribuidos entre los depósitos catastrados el año 2015, los cuales suman un valor
cercano a veinte depósitos. Por lo tanto no toda la variación corresponde a un crecimiento en la
construcción de nuevos depósitos, sin embargo, el número que se encontraba sin información no
es muy influyente con respecto a todo el crecimiento.
1.2. PRINCIPALES MECANISMOS DE FALLAS EN TRANQUES DE RELAVE.
Dentro de las principales fallas que se presentan en tranques de relaves se tienen como factores
el método de construcción, mala compactación alto contenido de finos, y saturación. Estos
factores en la fase de diseño, construcción y operaciones repercuten en fallas observadas como:
licuefacción, inestabilidad de taludes, deformaciones, asientos, overtopping, tubificación,
inestabilidad del suelo de fundación, planos de fallas por precipitaciones o erosión eólica y por
último colapso por aluviones.
1.2.1. Licuefacción o falla de flujo.
La licuefacción se puede entender como la excesiva cantidad de presión intersticial ocasionada

durante una carga dinámica, producida durante un movimiento sísmico. Este aumento de presión
de poros momentáneo, está en condiciones de velocidad de drenaje por debajo de la velocidad
de aplicación de la carga, lo que ocasiona que la resistencia al corte de las partículas tienda a
cero, produciendo en pocos minutos grandes deformaciones de la masa de suelo [1], [2].
Este tipo de falla es la principal en tranques de relaves en Chile, ya que cumple con todas las
características para que ocurra, uniformidad en la granulometría y redondez de las partículas
material, añadiendo factores de operación como mala compactación, estado permanente de
saturación y por sobre todo la carga cíclica al considerar a Chile un país con comportamiento
activamente sísmico.
Para evaluar este fenómeno es posible a través de ensayos triaxiales monótonos y cíclicos en
condición no drenada sobre muestras representativas del relave. Con los resultados de los
ensayos en forma cíclica se obtiene la resistencia cíclica del relave, con la cual es posible estimar
el potencial de licuefacción del depósito.
En la Figura 1.5. se visualiza las características de la licuefacción en tranques de relave.

Figura 1.5. Visualización de licuefacción en tranques de relaves.
Fuente. Villavicencio, Espinace, Palma, Fourie, & Valenzuela, 2014, pp.449-469.
Existirá una subdivisión para este tipo de falla, la cual constará de la falla de flujo o licuación
verdadera y la movilidad cíclica. Sin el aumento en la presión de poros no existirá licuefacción,
será el responsable de que ocurra ya sea como falla de flujo o movilidad cíclica. Para una mejor
visualización de estos conceptos, se presenta en la Figura 1.6. un gráfico que muestra cuando
ocurre cada uno de estos casos.
Figura 1.6. Falla de flujo y movilidad cíclica.

Fuente. Asenjo,(2010).

Ocurre falla de flujo, cuando el corte estático inicial esté sobre la resistencia última no drenada,
como el caso de A en la Figura 1.6., ya que el suelo no resiste y fluye. El segundo caso ocurrirá
cuando el corte estático inicial sea inferior a la resistencia última no drenada, es decir, como el
caso B, la movilidad cíclica se producirá por el efecto de una carga cíclica que va generando
progresivamente un aumento en la presión de poros.
Existen factores que afectan la resistencia cíclica de las arenas, los más importantes se nombran
a continuación.
– Densidad. Se demostró que una arena en estado suelto, requerirá una menor cantidad de
ciclos que una arena densa, para alcanzar el estado de licuefacción, a un mismo nivel de corte
cíclico solicitante [3]. También el suelo es capaz de obtener una presión efectiva igual a cero
sin deformaciones excesivas y continuar aumentando progresivamente la deformación. Por el
contrario si se encuentra en estado suelto, cuando alcanza la presión efectiva igual a cero
también ocurren deformaciones excesivas.
– Presión de Confinamiento. Otro factor importante es la presión de confinamiento, la cual es

explicada a través de un gráfico que muestra la relación entre la razón de vacíos y la presión
efectiva de consolidación.
Figura 1.7. Relación Vacía – Presión de consolidación.

Fuente. Asenjo,(2010).
Se pueden apreciar dos casos, A y B, los cuales representan el estado inicial del suelo para la
ejecución de los ensayos. El gráfico indica que luego de consolidar, ambos casos obtienen la
misma densidad y resistencias no drenadas casi iguales (A= 20,4 y B= 20,8). El ensayo A para

una presión de consolidación de 2 [kg/cm²] necesita cincuenta y siete ciclos para sufrir
licuación, en cambio el ensayo B requiere de diez mil ciclos para licuar con una presión de
consolidación de 12 [kg/cm²]. Por lo tanto, se concluye que para un corte cíclico solicitante de
la misma magnitud, mientras mayor sea la presión efectiva de confinamiento se necesitan
cantidades superiores de ciclos para que el material licue, es decir, la resistencia a tensiones
cíclicas es mayor al aumentar su presión efectiva [3].
– Contenido de finos. Éste factor tendrá importancia en la resistencia cíclica de las arenas, ya
que por el contrario los suelos finos gracias a su cohesión provocan que los granos sean más
difíciles de separar, por lo tanto aumenta la resistencia a la licuefacción. Esto se demuestra en
la Figura 1.8.
Figura 1.8. Efecto del índice de plasticidad en la resistencia cíclica de arenas con presencia de
finos.
Fuente. Ishihara, (1996).
Es importante mencionar que la presencia de finos no plásticos con acompañamiento de

minerales no inducen este mejoramiento en la resistencia, por el contrario hacen el material
más susceptible a que sufra licuefacción. Un ejemplo de esto es el relave, ya que la presencia
de estos minerales logra que sea más fácil la separación de las partículas, por ser una
superficie seca [4].
– Frecuencia de la carga cíclica. Este elemento será importante considerarlo, ya que existirán
movimientos sísmicos de distintas magnitudes. Luego de estudios se concluye que la
resistencia a la licuefacción no es influenciada por la frecuencia de la carga cíclica para rangos
entre 0,1 [Hz] y 20 [Hz] [5].

1.2.2. Inestabilidad del suelo de fundación.
La capacidad portante del suelo de fundación es la capacidad que tiene el suelo de fundación
para soportar las cargas de trabajo que solicite una estructura sobre él. Cuando ve sobrepasada
dicha capacidad portante el suelo sufre fallas por rotura o asentamientos.
Es por ello que al generar grandes cargas de trabajo por el peso del muro de relave, pueden
generar asentamientos diferenciales en el suelo de fundación o planos de rotura considerables
generando inestabilidad global en todo el depósito del tranque de relave.
Figura 1.9. Inestabilidad Suelo de Fundación.

Fuente. Espinace, (2008).
1.2.3. Inestabilidad de taludes.
Un movimiento ocurre cuando la resistencia al esfuerzo cortante del suelo es excedida por los
esfuerzos cortantes que se producen en una superficie relativamente continua. Por lo tanto, las
fallas localizadas en un solo punto de la masa de tierra no indican necesariamente que la masa
sea inestable. La inestabilidad se produce como resultado de la falla de esfuerzo de corte en una
serie de puntos que definen una superficie, a lo largo de la cual se produce el movimiento. Esta
falla puede ser del tipo superficie plana (deslizamiento plano), rotura circular o rotura compuesta
[6].
Entre los factores que podemos distinguir de la inestabilidad de los taludes se pueden clasificar
en internos y externos.
– Internos.
‒ Capa de formación del tranque de relave que presenten algún grado de heterogeneidad en

el grado de compactación la cual ocasione comportamiento desigual o que presente alguna
variación en la calidad de sus propiedades mecánicas.
‒ Aumento de la carga de presiones hidrostáticas, producto del mal diseño del sistema de
drenaje producto de aumentos de nivel freático ocasionado por lluvias excesivas.
– Externos.
‒ Sobrecargas sísmicas.
‒ Sobrecargas.
‒ Geometría inapropiada.
Cabe señalar que aunque la saturación en la masa no sea la necesaria para originar una patología
de licuación sísmica, el movimiento dinámico producido por la carga sísmica puede ocasionar
zonas de debilitamiento, reduciendo las resistencias al corte y el reacomodo de las partículas
ocasionando zonas fallas y colapso de la estructura interna del tranque de relave.
Figura 1.10. Falla circular de taludes.

Fuente.Espinace, (2008).
1.2.4. Deformaciones.
Las deformaciones que se producen en los tranques de relaves pueden provocar inestabilidad
posterior en taludes del mismo, o muchas veces ocasionar el colapso total de la estructura
soportante.
Entre los tipos de asientos se encuentran: asientos por contracción, deformaciones por tensión
de tracción y baja resistencia al corte por condiciones sísmicas.

Los asientos por contracción corresponden a la disminución del contenido en agua, que hace
contraer el material hasta su límite de plasticidad, y luego agrietarse, esto es producto que el
contenido de humedad al momento de ser depositadas las arenas va a ir disminución producto
del drenaje natural o la evaporación por efectos del cambio de temperatura.
Las deformaciones bajo tensiones de tracciones son causadas por la propagación de las ondas
sísmicas, particularmente al nivel de la corona donde la cantidad de materiales es mínima y pues
donde la resistencia por confinamiento es débil. Las fisuras se propagan sobre 1 a 2 [m] de
profundidad y no constituyen riesgo potencial cuando son longitudinales. Cuando ondas de
compresión sísmicas aparecen en el eje longitudinal de la obra, tensiones de tracción pueden
crear las fisuras perpendiculares a la corona, en el eje de la pendiente. Estas fisuras son
peligrosas, porque pueden crear zonas de rotura por inestabilidad de pendiente [6].
Las deformaciones por tensión de corte sísmico corresponden al desplazamiento de las

superficies agrietadas por las acciones de una aceleración sísmica. Fisuras de 1 [cm] a 1 [m] se
producen al principio de las superficies de deslizamiento pero el bloque en movimiento no se
desprende completamente. La amplitud del desplazamiento del bloque depende de la duración
de la aceleración sísmica [7].
Figura 1.11. Grietas producto de deformaciones.

Fuente. Villavicencio, (2014).
1.2.5. Tubificación o piping.
La tubificación o piping puede erosionar tanto los terraplenes de las presas, ya sea de tierra como
de relaves, como los sitios de paso de las estructuras de hormigón, socavando los terraplenes o

sus fundaciones.
El flujo de agua surge por una diferencia de energía (dada por la carga o diferencia en la cota de
agua) entre dos puntos, ejemplificando esto en la diferencia que existe entre el nivel de la cubeta
en tranque de depósitos de relaves y el muro aguas abajo.
El flujo tiende a ejercer presión sobre las partículas del suelo, y la fuerza de filtración actuante,
definida por unidad de volumen, puede arrastrarlas en el sentido de la corriente.
Las fuerzas que se oponen al arrastre del flujo es la trabazón entre las partículas (estructura) y
su cohesión, pero en el caso particular de muros hechos con arenas de relaves no existe cohesión
y la trabazón de las partículas se ve limitada por la naturaleza del suelo al ser arenas de fracción
fina y limos, además de deficientes procedimientos constructivos lo que lleva a una mala
compactación o deficiente sistemas de drenaje, generan como producto que el agua se infiltre
por las grietas internas ocasionando arrastre de finos, produciendo arrastre de finos y posterior
colapso de la estructura soportante.
Figura 1.12. Esquematización de Tubificación en presas.

1.2.6. Overtopping.
El Overtopping o desbordamiento puede producirse por dos causas: uno es por desbordamiento
del fluido ya sea agua o relave por un aumento en el caudal por precipitaciones intensas o un
error de depósitos de residuos, a causa de una deficiente construcción del muro resistente o de
las cantidades de volumen incorporado de relave a la cubeta.
Al ser sobrepasado el nivel más alto del muro resistente (cresta o playa), puede ocasionar zonas
de fallas en talud, además de generar una condición de saturación del talud de aguas abajo
generando aumentos en la probabilidad de riesgos de licuefacción.

Figura 1.13. Overtopping.
1.2.7. Aluviones.
Los aluviones corresponden a un tipo de movimiento brusco de suelo o roca mezclado con agua.
Se caracterizan por sus flujos rápidos y violentos capaces de arrastrar rocas y otros materiales
que descienden por una quebrada o lecho de río. Éstos ocurren cuando el agua se acumula
rápidamente en el suelo a raíz de una lluvia intensa o deshielos repentinos, convirtiendo el terreno
en un caudaloso río de lodo o barro (USGS).
Estas corrientes fluyen rápidamente por una quebrada, destruyendo todo a su paso con poca o
nula advertencia. Pueden extenderse varios kilómetros desde su punto de origen, aumentando
considerablemente de tamaño a medida que arrastran árboles, rocas, y otros materiales que
encuentren en su recorrido.
Esto se coloca en el contexto de los últimos acontecimientos de lluvias de gran intensidad que
han afectado al norte de Chile, ocasionando arrastres de masas de suelo y roca que destruyen
todo a su paso. Es así como pueden ocasionar fallas y colapsos totales en tranques de relaves,
al ser sometidos a solicitaciones externas no contempladas en diseño, produciendo desastres
catastróficos para la naturaleza y las personas.

Figura 1.14. Escurrimiento flujo con suelo, rocas y varios.
Fuente. http://www.sochige.cl/, (2015).
1.3. FALLAS HISTÓRICAS EN TRANQUES DE RELAVES.
A lo largo de la historia de la minería en Chile, han ocurrido numerosas fallas en los tranques de
relaves, producto de inestabilidad mecánica. Los eventos más importantes registrados fueron los
años 1928 en el Tranque Barahona de la mina El Teniente en Rancagua; año 1965 Tranque el
Cobre de la mina El soldado en Nogales; y el año 2010 los principales tranques con fallas fueron
Tranque Las Palmas en Pencahue y Tranque adosado en Alhué, ambos casos en la región del
Maule. Todos estos hechos fueron producto de sismos de gran magnitud que ha sufrido Chile.
1.3.1. Tranque Barahona.
El tranque Barahona se emplaza en la división El Teniente de CODELCO, en la ciudad de

Rancagua. El hecho ocurrió luego del sismo del 1 de octubre de 1928, con epicentro en la ciudad
de Talca y una magnitud de ocho grados en la escala de Richter. El movimiento produjo
aproximadamente el avance de cuatro millones de toneladas, las cuales arrasaron con el valle
aguas abajo, matando a cincuenta y cuatro personas.
Las fallas detectadas fueron licuefacción con falla de flujo, producida por las solicitaciones de
esfuerzo dinámico del sismo.

Figura 1.15. Colapso del muro del tranque Barahona, producto de terremoto de 1928.
Fuente. Http://www.relaves.org/, (2015).
1.3.2. Tranque El cobre.
Ubicado a unos 10 [km] de la calera cerca de la hacienda El Melón. El tranque cede al primer
movimiento sísmico, durante el terremoto del año 1965 en la ciudad de la Ligua, con una magnitud
de 7,6 grados Richter, ocasionando que la mezcla de minerales, de alrededor de diez millones de
metros cúbicos, avanzaran aguas abajo, cubriendo con una capa de hasta 5 [m] al poblado
colindante, lo que deja casi en un 100% víctimas fatales. Con este terremoto otros tranques se
vieron afectados, pero sin ocasionar víctimas fatales, como: los Maquis, Cerro Negro, Bellavista,
La africana, El cerrado y La patagua. Las causas del colapso fueron fallas de licuefacción con
flujo e inestabilidad de taludes con deformaciones sísmicas inducida.
1.3.3. Tranque Las Palmas y Tranque Adosado.
En febrero del año 2010 ocurre el sexto terremoto más grande registrado en el último siglo, el
cual tuvo su epicentro en la comuna de Cobquecura, con una magnitud de 8,8 grados Richter.
Dentro de los casos más importantes, el que tuvo mayor daño fue el Tranque Las Palmas, ubicado
en la región del Maule. Éste colapsó, desplazándose aproximadamente 500 [m], donde la causa
fue falla de licuefacción y flujo. La primera de éstas, ocurre en la base del muro, producto de la
saturación del subsuelo. A causa de este colapso fallecen cuatro personas aledañas al sector del
tranque.
Con respecto al caso del tranque Adosado, éste también tuvo como causante de falla licuefacción

de las arenas del muro, el cual había sido construido bajo el método de aguas abajo. Sin embargo,
la ejecución de su compactación fue deficiente. La geometría de éste era de altura entre 8 y 13
[m], con una pendiente de relación 1:4 (V: H). Según información, las principales causas de
licuefacción fueron las siguientes.
‒ Un alto grado de saturación, debido a la filtración proveniente de la laguna de aguas claras y

debido al sistema de rociadores utilizado como una medida de mitigación para evitar la erosión
del viento.
‒ A través de ensayos de penetración con el equipo PANDA, se obtuvo que el estado de

compactación de las arenas fue clasificado como muy baja a baja (densidad relativa menor a
45%), con una alta probabilidad de sufrir licuefacción.
1.4. PROGRESO DE LA NORMATIVA CHILENA EN TRANQUES DE RELAVE.
A partir de los eventos de sismos ocurridos y sus respectivas fallas en la estructura de los tranques
de relaves, surge la necesidad de mejorar la normativa vigente en ese entonces. En primera
instancia, luego del terremoto del año 1965 nombrado en el apartado anterior, es cuando se
requiere una evaluación de los métodos de construcción para tranques, el cual es anunciado el
año 1970. Posteriormente el año 2007 es modificado el decreto supremo n° 86 del año 1970, el
cual es reemplazado por el decreto supremo n°248 denominado Reglamento de Construcción y
Operación de Tranques e Relaves.
1.4.1. Decreto Supremo N°86.
Se origina a partir de la necesidad de precisar ideas y exigencias tecnológicas para aplicar los
conceptos más avanzados en la construcción de los tranques de relave. Hasta esa instancia, el
método más utilizado para la construcción de los muros de tranques de relaves era utilizando las
arenas provenientes del mismo relave, siendo además una opción más económica. Sin embargo,
la mala ejecución y operación de éstos provocó las diversas fallas luego del sismo. Por esta razón
es que se profundiza en el comportamiento interno de los tranques; recomendaciones
constructivas y control, para incluirlo en el decreto.
Los principales requerimientos establecidos son los siguientes.
‒ En relación a la mantención y operación, se especifica que los usuarios deberán entregar en

forma periódica.
‒ Información de los trabajos de mantención mensualmente.
‒ Información de los trabajos de operación mensualmente.

‒ Densidad relativa a ⅓ y ⅔ de altura del prisma resistente cada 100[m] lineales mínimo en
forma mensual.
‒ Granulometría sobre 48, 60,100 y 150 mallas Tyler y bajo 150 mallas Tyler en porcentaje, una
vez al mes.
‒ Subtotales de arenas y lamas agregadas en el mes, año y total.
‒ Para la construcción del prisma resistente se aceptarán sólo dos alternativas.
– Prisma resistente que se construye colocando la fracción arenosa hacia afuera de la

acumulación lamosa. Al ir incrementando la cota el equipo repartidor se irá desplazando
hacia el exterior de la acumulación lamosa. Este sistema se denomina aguas afuera.
– El prisma resistente se construye manteniendo el coronamiento dado, mientras va
aumentando la cota se mantiene la vertical sobre ese mismo punto. Este sistema se
denomina mixto.
‒ El vaciamiento de las arenas será depositándolas al exterior sobre el talud externo,

extendiéndolas de la manera más uniforme posible, compactando con capas no mayores de
30 [cm], en condiciones de humedad cercanas a la óptima que indique el ensayo Proctor, hasta
lograr una densidad relativa o un porcentaje de Proctor máximo igual o superior al límite fijado.
‒ Los tranques de relaves deberán ser diseñados de acuerdo a los principios de la Mecánica de
Suelos e Ingeniería de Fundaciones.
‒ La necesidad de conseguir un aumento de la permeabilidad de la base de apoyo del prisma

resistente, mediante un sistema de drenaje en la base, piezómetros para el control.
‒ Se utiliza un factor de seguridad igual o superior a 1,2.
‒ Para la seguridad de los poblados colindantes aguas abajo de la zona lamosa, se establece
una distancia peligrosa, la cual está determinada por la pendiente y la distancia que recorrerían
las lamas y arenas producto de licuefacción sísmica.
En relación a la geometría, monitoreo y sistema de revestimiento no existen especificaciones.
A pesar de los perfeccionamientos de la normativa y los controles realizados con ésta,

nuevamente surgieron daños en los terremotos ocurridos posterior a la puesta en práctica de este
decreto. Una de las principales causas fue que no se dejó claro la prohibición del método de
construcción aguas arriba, en el cual la secuencia de crecimiento del prisma resistente se ejecuta
hacia el interior de la zona lamosa, quedando sobre una base de mala calidad. Fallas importantes
ocurrieron luego del terremoto de 1981, 1985 y 1997, lo que se traduciría en un nuevo decreto.
1.4.2. Decreto supremo N°248.

El año 2007 entra en vigencia este nuevo reglamento, el cual es actualmente por el que se rigen
los tranques de relaves en Chile. Surge por la necesidad de proteger la salud y seguridad tanto
de las personas como del medio ambiente y recursos naturales. Además los tipos de depósitos y
métodos de diseño y construcción han tenido un avance en el tiempo, por lo que es necesario
precisar las exigencias técnicas para éstos. Éste se denomina Reglamento para la Aprobación de
Proyectos de Diseño, Construcción, Operación y Cierre de los depósitos de relaves.
En primera instancia para presentar el proyecto para la aprobación de las organizaciones

correspondientes, se especifican los siguientes antecedentes con respecto a la construcción del
prisma resistente.
‒ Capacidad final del depósito en metros cúbicos y toneladas, el ritmo en que se deposita, vida
útil y crecimiento anual. Para tranques de relaves además incluir el porcentaje del relave que
será utilizado en el muro.
‒ Antecedentes geológicos, geotécnicos, hidrológicos, hidrogeológicos, sísmicos,

meteorológicos y topográficos del sector de emplazamiento.
‒ Descripción del método constructivo, utilizando sólo los métodos de construcción de eje central
o aguas abajo, dejando en forma clara la prohibición del uso del método aguas arriba. Además
las características del muro de arena, granulometría, humedad y forma de deposición de las
lamas. También definir el valor mínimo de compactación para el ensayo Proctor y el porcentaje
de arenas bajo la malla 200 que se puede utilizar en el muro.
‒ Descripción de las dimensiones del depósito tanto en altura y largo de muro, como de área y
volumen del depósito, como también su Plan de Crecimiento.
‒ Determinación de los parámetros de diseño geotécnicos, incluyendo lo que sigue.
– Resistencia al corte.
– Compresibilidad.
– Permeabilidad.
– Granulometrías.
– Pesos unitarios.
– Pesos específicos.
– Plasticidad.
‒ Para el control y monitoreo del comportamiento hidráulico y estructural del interior del prisma,
se utiliza lo siguiente.
– Presiones de poros.
– Niveles freáticos.
– Desplazamientos.

– Asentamientos.
– Filtraciones.
– Aceleraciones sísmicas.
‒ Análisis de estabilidad de taludes para las etapas de operación y cierre. Para esto se divide
en cuatro diferentes fases, según los posibles riesgos que pudiera presentar el depósito. Las
fases son las siguientes.
– Fase I. Análisis Pseudo-estático, considerando licuefacción total de los relaves de la cubeta.

– Fase II. Análisis Pseudo-estático, con determinación simplificada de las presiones de poros.
– Fase III. Análisis dinámico, basado en ensayos dinámicos de suelos y cálculo de
desplazamientos.
– Fase IV. Análisis para el cierre, considerando eventos solicitantes máximos y efectos del
tiempo en las propiedades del depósito.
Para las fases I y II el factor de seguridad obtenido debe ser menor a 1,2. Cuando el muro
tenga una altura inferior a 15 [m], si cumple con el factor de seguridad no será requisito la fase
III.
‒ Para cálculo dinámico el sismo de diseño debe obtenerse a partir de estadísticas de la zona
correspondiente y estimar su aceleración máxima respectiva.
‒ Determinación de la distancia peligrosa y el análisis de ésta en terreno y además incluir

plano de la hoya hidrográfica que pudiera afectar la trayectoria de ésta.
Esta nueva normativa fija ciertos criterios respecto al control de los depósitos de relaves. Los
puntos más importantes se mencionan a continuación.
‒ En los tranques de relaves, la laguna de aguas claras deberá estar lo más alejada del muro
para evitar aumentar la humedad del muro y su saturación, ya que esto se traduce en un
aumento de la presión intersticial de la estructura interna y posible colapso por licuefacción.
Será necesario un constante drenaje de las aguas claras.
‒ La revancha deberá ser, como mínimo, un metro. Sin perjuicio de considerar los fenómenos
climáticos que exigieren una mayor revancha.
‒ El coronamiento debe tener un ancho que asegure la estabilidad del muro, el cual debe ser,
como mínimo, dos metros.
‒ El proyecto deberá incluir la descripción de las obras, planos, memorias de cálculo,

especificaciones de construcción, planes de control y monitoreo, normas de operación, etc.
Las etapas del proyecto deben ejecutarse en forma secuencial, continua e integrada, de
manera de permitir la evaluación, la realimentación y la optimización de las obras
‒ El muro deberá contar con un sistema drenante en su base. Además en su inicio deberá tener

una altura mínima de 1/10 de la altura final proyectada y con un mínimo de dos metros.
También deberá poseer un sistema de impermeabilización en su coronamiento y talud interior.
‒ El porcentaje de arena más gruesa de los relaves, usada para la construcción del muro en un
tranque, no deberá ser más de un 20% de partículas menores a la malla 200.
‒ El sistema de impermeabilización del fondo de la cubeta deberá contemplar un tratamiento

previo del terreno, con el fin de impedir filtraciones de aguas contaminadas al suelo o aguas
subterráneas. En tranques de relaves se prohíbe el uso de geomembranas.
1.5. VARIABILIDAD DE LOS PARÁMETROS DEL CONTROL DE COMPACTACIÓN EN

TRANQUES DE RELAVES.
El estudio de las arenas de relaves que provienen de los residuos mineros, se pueden clasificar
en una única clase. Sin embargo, en la práctica existe una variación de las características físicas
y geotécnicas. Entonces para adaptar una condición de homogeneidad es posible ejecutar un
análisis estadístico de aquellos parámetros. Mediante este análisis será posible acotar la
variación de las variables nombradas y detectar en forma previa el comportamiento global del
tranque.
Existen dos tipos de variabilidad en las arenas de relaves; la variabilidad material y la variabilidad
estructural. La primera hace referencia a las propiedades índices del material y la estructural al
proceso de construcción empleado.
‒ Variabilidad material. Se refiere al origen mineralógico y el proceso de fabricación de las

arenas de relave (molienda, flotación, ciclonaje y depositación), es decir, una variación en la
densidad de las partículas, granulometría y plasticidad de las partículas finas.
‒ Variabilidad estructural. Tiene relación con el proceso de construcción empleado y su

compactación, teniendo como factores el grado de compactación y el espesor de capa.
En ambos casos se subclasifican en dos análisis; análisis a escala espacial, el cual se asocia a
la dispersión de las características geotécnicas de las arenas de relaves y la estructura en el
espacio. El segundo análisis será a escala temporal, la cual es fundamentalmente el
envejecimiento que sufren estos materiales, producto de la cementación de las partículas,
cambios en el estado hídrico, entre otros [8]. La Figura 1.16., muestra una tabla con el resumen
de esta clasificación.
Efectos sobre los residuos

Origen Tipos de variabilidad
mineros

Uno o varios Homogeneidad o heterogeneidad de
yacimientos mineros los residuos
Tipo de mineral y calidad Variabilidad de los minerales

mineralógica constituyentes
Material
Distribución de las tallas de las
Especie mineralógica Espacial
partículas
Contenido de corte del
Y añadido de polímeros y de cal,
yacimiento
características mineralógicas
Metodología de ciclón y Variación de las características

Estructural
de construcción físicas y mecánicas
Efecto del
Cohesión por cementación Material Temporal
envejecimiento
Figura 1.16. Tipos de variabilidad existente en los tranques de relaves.

Fuente. Durant, (2012).
1.5.1. Metodología de análisis de la variabilidad.
Para el estudio de la variabilidad de las propiedades físicas y geotécnicas en tranques de relaves,

será necesario contar con una batería de datos y gestionarlas de forma óptima. Para esto se
requerirá una bibliografía con relación a las características de los materiales usados en tranques
de relaves y a los análisis estadísticos empleados en tranques de relaves. Para el desarrollo de
esta metodología es necesario realizar los siguientes pasos.
‒ Recopilación de información de los tranques de relaves a analizar, como ubicación de

explotación, origen mineralógico, método constructivo, compactación, entre otros.
‒ Gestión de datos del control de compactación, los cuales serán parámetros de naturaleza, de
estado y resistencia a la penetración. Los de naturaleza son obtenidos mediante ensayos
realizados en laboratorio; los de estado se refieren a los datos obtenidos in situ y la resistencia
a la penetración mediante ensayos de penetrómetros ligeros in situ.
‒ Estimación de los parámetros estadísticos descritos en el punto anterior, entregando valores

máximos, mínimos, promedio, desviación estándar, coeficiente de variación, siendo éste último
el principal dato que mostrará la variación del tranque.
‒ Construcción de histogramas de frecuencia para cada uno de los parámetros analizados.
‒ Finalmente se adopta una ley de distribución, obtenida como base de los histogramas y

pruebas estadísticas (Kolmogorov-Smirnov).
1.5.2. Descripción de los tranques de relaves experimentales.
Los tranques de arenas de relaves experimentales en los que se basa el estudio de la variabilidad,
presentan características típicas de los depósitos de la minería chilena y sus métodos de
construcción empleados. Para esto, se definen tres tranques de estudio, los cuales se presentan
a continuación. En la tabla correspondiente a la Figura 1.17., se muestran las características
generales con respecto al sector que pertenecen (estatal o privado), descripción de las arenas
del tranque y características del muro resistente.
Características del muro

Industria minera Arenas de relaves
resistente
Tranque Origen Talud Talud Alt

Producción
Escala Sector mineralógi Tipo A/A A/AA act.
[ton/mes]
co (V:H) (V:H) (m)
Diversos T/AV
1 Mediana Estatal 100.000 1 :3.5 1:2 27.0
yacimientos et E/C
Diversos
2 Mediana Estatal 60.000 E/C 1 :2.7 1:2 20.0
yacimientos
Yacimiento
3 Mediana Privado 110.000 E/C 1 :3.0 1:2 25.0
único
Figura 1.17. Resumen datos tranque experimentales para el análisis de la variabilidad.

Luego se presenta un resumen con los principales elementos de la metodología de construcción

empleada en cada uno de los tranques seleccionados.
Los controles de compactación realizados durante la fase operacional de cada tranque y las
campañas de ensayos de penetración para el control geotécnico, se resumen a continuación en
la tabla de la Figura 1.18.

Control
Control de compactación
geotécnico
Tranque Nivel de
Densidad Ensayo Ensayo de
Granulometría compactación
penetración
especificado
in –situ Proctor
DCPT et
1 x x x 95% OPN
PANDA
DCPT et
2 x x x 95% OPN
PANDA
3 x x x 95% OPN PANDA
Figura 1.18. Controles de compactación realizados a los tranques experimentales.
1.5.3. Análisis de la variabilidad de las características físicas y resistencia in situ.
En primer lugar se definen cuáles serán las características del tranque para evaluar su
variabilidad. Con respecto a las características físicas, éstas son obtenidas desde el control de
calidad realizado en la etapa de operación del tranque, donde los parámetros serán granulometría
(porcentaje finos y diámetro medio partículas), densidad seca y humedad in situ, densidad
máxima y humedad óptima Proctor. Las características de resistencia in situ, corresponderán a
los datos obtenidos por ensayos de penetración dinámica.
Para cada una de estas variables se realizará un análisis estadístico, a partir del cálculo del
coeficiente de variación (CV%), el cual muestra en forma excelente la dispersión de las
propiedades del suelo [9]. Luego de esto, se escogerá una función de distribución, la que mejor
se adapte a la variación espacial de cada parámetro.
Para el desarrollo del análisis estadístico, se deben clasificar los parámetros en tres niveles:
parámetros de naturaleza, parámetros de estado y parámetros de penetración.
– Parámetros de naturaleza. Éstos como se menciona previamente son los datos obtenidos a
partir de ensayos de laboratorio. En este caso corresponden al porcentaje de finos (menos a
80 [µm] ), diámetro medio partículas (D50), diámetro máximo (Dmáx), densidad máxima seca
(γdmax) y humedad óptima Proctor (Wop).
Luego del cálculo de la variación de los parámetros, se pudo observar que los tranques n°1 y
n°2 son los que presentan mayor coeficiente de variación, por lo que esto demuestra que han
sido construidos a partir de arenas de diferentes orígenes mineralógicos y sometidos a
procesos distintos de molienda y flotación. Sin embargo, el tranque n°3 presenta índices de

variación más pequeños lo que se traduce en que las arenas son más homogéneas, teniendo
un único origen mineralógico y los mismo procesos.
– Parámetros de estado. Los parámetros medidos in situ fueron la densidad natural medida
mediante cono de arena, humedad natural y densidad seca asociada.
Los resultados obtenidos mostraron una clara relación entre d y dmax, lo cual demuestra que
el proceso de compactación ha sido muy uniforme y los dos ensayos no presentan una
variación suplementaria. Sin embargo, estos controles son hechos en forma superficial, por lo
que no se puede concluir que existe ausencia de dispersión en mayor profundidad asociado a
la compactación.
– Resistencia a la penetración. Los datos fueron obtenidos a partir de ensayos de penetración

DCPT y PANDA, los valores arrojados son NDCPT y qd, respectivamente.
Los coeficientes de variación bordean el 50%, lo que puede deberse a distintas causas del
origen de esta dispersión. Una opción es que los ensayos se realizan en forma puntual, por lo
que puede ser una dispersión natural por el nivel de escala. Además puede existir una
variación en profundidad de este parámetro por la estructura inter-estratificada que presentan
los tranques, debido al proceso de construcción. Como tercera opción, la resistencia puede
presentar variación en profundidad y en función del tiempo, producto del efecto de
envejecimiento.
Estos ensayos permiten tener una imagen de la variabilidad espacial de las propiedades
mecánicas en un determinado momento.
1.5.4. Distribución estadística de los parámetros físicos y mecánicos.
La distribución de los parámetros en un macizo de suelo no es conocido perfectamente, ya que

es imposible determinarlo en cada punto. El uso de leyes de distribución experimentales no es
práctica, por esto se modela a partir de leyes de distribución teóricas [9].
Para modelar los parámetros geotécnicos mencionados previamente se utilizaron variables

aleatorias continuas (X), de promedio (u) o de esperanza matemática E(X) y de varianza [σ²].
Se realizaron histogramas representando los parámetros del estudio, en los cuales diferentes
autores han señalado que el logaritmo de los valores, siguen una distribución normal [10]. Para
el estudio de estos tranques experimentales se analizaron los valores log (qd) y log (NCPT60).
Para cada tranque se realizaron histogramas en el análisis de sus parámetros. En la Figura 1.19.,
y Figura 1.20., se muestran a modo de ejemplo algunos histogramas de la zona de estudio.

400 400
Frecuencia
Frecuencia
300 300
200 200
100 100
0 0
1,4 1,6 1,8 2 2,2 2,4 1,5 1,75 2 2,25 2,5
Densidad seca in-situ (grs/cc) Densidad natural (grs/cc)
Figura 1.19. Histogramas de frecuencia parámetro densidad en tranque n°1.

Figura 1.20. Histograma de frecuencia resistencia por punta en tranque n°3.

Los histogramas arrojados a partir de las densidades se presentan de dos formas; densidad seca
in situ y densidad natural, mostrando el primero de éstos un comportamiento heterogéneo que se
observa a través de dos grandes familias de materiales. El gran tamaño de la muestra permite
representar en forma global la dispersión de las arenas del tranque n°1. La heterogeneidad se
puede explicar por posibles factores.
‒ El tranque es del sector estatal de la mediana y pequeña minería, por lo tanto los minerales
provienen de distintos yacimientos. Sin embargo, esto no implica que provengan de sólo dos
familias.
‒ El proceso de ciclonaje es el responsable del porcentaje de finos. El resultado puede ser

asociado a dos metodologías usadas en el tranque, pero no existe registro de esto.
‒ La metodología de compactación no ha sido realizada en forma homogénea, pero esto no

explica la variación observada en estos parámetros. No obstante, la metodología de
compactación puede ser una variable importante en el estado de los materiales depositados.
Esto se aprecia en la gran dispersión de la resistencia a la penetración.
A diferencia de los tranques n°1 y n°2, el tranque n°3 presenta una alta regularidad, simétricos
con una única familia y reducida dispersión, es decir, se considera un material homogéneo en
relación a la naturaleza y depositación de las arenas. Esto ocurre por los siguientes factores.
‒ El tranque es del sector privado, con arenas provenientes de un solo yacimiento. Es por esto
que la mineralogía es más constante.
‒ El proceso de construcción es adecuadamente controlado, con depósitos de arenas en capas

uniformes y compactadas mecánicamente a través de rodillos liso vibratorios.
Como último paso de esta metodología, luego de realizar histogramas, éstos deben ser
adaptados a una ley de distribución teórica. Para esto se analizan los parámetros según las
familias de datos formadas, a través de coeficientes que muestran la asimetría y aplastamiento.
Los parámetros se adoptan a una ley de distribución normal, donde los datos de resistencia a la
penetración son ingresados con el logaritmo de los valores. En la Figura 1.21., se muestra la tabla
resumen.
Parámetro Tranque Nº 1 Tranque Nº 2 Tranque Nº 3

% des finos Normal Normal Normal
D50 Normal Normal Normal
d No Normal Normal
h No Normal Normal
log (qd) Normal Normal Normal
log (N60) Normal Normal No
qd Log-Normal Log-Normal Log-Normal
N60 Log-Normal Log-Normal Normal
Figura 1.21. Resumen datos tranque experimentales para el análisis de la variabilidad.
Luego de todo el análisis de variabilidad de los tranques experimentales realizado por

Villavicencio, se realiza un resumen en la tabla de la Figura 1.22., que se muestra a continuación.

Tranque 1 2 3
Escala Mediana Mediana Mediana
Minería
Sector Estatal Estatal Privado
Diversos Diversos
Origen Yacimiento Único
yacimientos yacimientos
CV (%) 28,7 26,3 10 ,00
FC%
Ley Dist. Nor. Nor. Nor.
Propiedades
CV (%) 19,03 15,2 8,67
físicas D50
Ley Dist. Nor. Nor. Nor.
CV (%) 6,56 8,18 2,92

Ley Dist. - Nor. Nor.
Arenas Parámetros CV (%) 6,74 8,21 3,51
de estado 
de Ley Dist. - Nor. Nor.
relaves Propiedades CV (%) 50,63 45,89 52,75
de qd
Ley Dist. Log-Nor. Log-Nor. Log-Nor.
resistencia
CV (%) 62,49 58,89 -
a la N60
penetración Ley dist. Log-Nor. Log-Nor. -
Figura 1.22. Resumen del análisis de variabilidad de los tranques experimentales.
1.6. COMPACTACIÓN EN TRANQUES DE RELAVES.
La compactación es el procedimiento de aplicar energía al suelo eliminando los huecos,

aumentando así su densidad y en consecuencia, su capacidad de soporte y estabilidad, entre
otras propiedades.
Cabe señalar, que un suelo logre una eficaz compactación dependerá de la elección de la
maquinaria a emplear, esto se puede lograr sabiendo las características del mismo mediante
ensayos y procedencia de su composición mineralógica.
En Chile existen principalmente dos ensayos para determinar las densidades máximas o mínimas
que un suelo puede alcanzar dependiendo de su clasificación, ya sea granular o cohesivo, para
luego poder medir el grado de compactación requerida. Por un lado, existe el ensayo que
relaciona densidades v/s humedades, con una cierta energía de aplicación llamado usualmente
Proctor modificado, el cual permite obtener la densidad máxima seca a través de la humedad
óptima y la energía de compactación que puede alcanzar un suelo cohesivo (12% de partículas
mayor a 0,08 [mm]). Por otro lado para la determinación en suelos granulares, no cementados,
de flujo libre, se utiliza un procedimiento el cual permite conocer mediante compactación por

vibrado su densidad máxima y la densidad mínima mediante vaciado, que posteriormente
relacionándolas con la variable de compacidad en estado natural del suelo, entrega un valor
porcentual llamado densidad relativa, concepto que clasifica el suelo en cinco categorías de
manera cualitativa (muy suelto, suelto, medio, denso, muy denso.). De esta manera se puede
tener un parámetro de control para la compactación dependiendo del tipo de suelo empleado.
Se deben conocer los siguientes datos para la elección de la maquinaria adecuada en la

compactación de tranques de relave.
‒ Granulometría.
‒ Homogeneidad del suelo.
‒ Especificaciones del proyecto (densidad máxima, humedad óptima, número de pasadas, etc.).
‒ Tiempo necesario para realizar trabajo de compactación.
Para la compactación de las arenas de tranques de relaves se utiliza principalmente dos

maquinarias tipo rodillo compactador dinámico en tándem con bulldozer y rodillos
autopropulsados. Las características de peso estático, amplitud, frecuencia, varían de acuerdo
al tipo de depósito o tipo específico de modelo de maquinaria utilizada.
Figura 1.23. Rodillo liso vibratorio propulsado por Bulldozer para compactación muro tranque
de relave.
La metodología para el depósito y posterior compactación de arenas de tranques de relave en la

actualidad viene dada por los siguientes pasos.
‒ La depositación de arenas de relave se efectúa por medio de unas tuberías que transportan
dicho material para su colocación por debajo de la línea central, que es la ubicación del prisma

resistente o muro de partida del tranque de relave. Desde ahí se sitúan aguas abajo en paños
aproximadamente de 350 [m], y espesor de óptimo 30 [cm], pero que en muchas ocasiones
llegan a capas de 90 [cm], este espesor es controlado mediante levantamiento referencial
topográfico.
‒ La fracción más gruesa, es decir, las arenas a depositar, según decreto deben poseer un
contenido mayor al 20% sobre la malla N° 200 (0,08 [mm]).
‒ Antes de iniciar la compactación de las arenas, es necesario que se encuentren en su

humedad óptima, es por ello que una vez sean depositadas se espera una cantidad de tiempo
suficiente para que estas puedan drenar naturalmente y de esta forma alcancen humedades
cercanas al punto óptimo. Esta humedad óptima está definida por el ensayo de relación
humedad/ densidad (Proctor normal o modificado).
‒ La compactación se realiza primero pasando un rodillo compactador dinámico propulsado por

un Bulldozer, donde las orugas son de gran importancia. La dirección de compactación la
mayoría de las veces se efectúa perpendicular al eje del coronamiento del tranque, variando
en sentido paralelo sólo en el coronamiento de éste. El traslape que se efectúa entre pasada
es el mínimo, teniendo que garantizar la homogeneidad de la compactación. Muchas veces se
finaliza la compactación generando una pasada con un rodillo liso autopropulsado. El número
de pasadas se determina en canchas de prueba, donde se analiza la cantidad necesaria para
alcanzar la densidad requerida en proyecto.
‒ Si después de haber aplicado las pasadas correspondientes, al tomar el control no se llega a

densidades mayores o iguales al 95% (DMCS), se debe generar una pasada adicional para
así llegar a la densidad requerida. Si se verifica que la humedad in situ está por debajo de
humedades óptimas de compactación será necesario regar la cancha y volver a aplicar todo
el procedimiento completo anteriormente descrito.
Los controles utilizados al finalizar la compactación de arenas de tranques de relaves, serán los
dispuestos en el decreto supremo N° 248, de 2007 Ministerio de Minería (Reglamento para la
Aprobación de Proyectos de Diseño, Construcción, Operación y Cierre de los Depósitos de
Relaves), y los que estipule pertinente el ente a cargo de la faena de compactación.
Con relación a los controles dispuestos por el decreto supremo N° 248, de 2007 Ministerio de
Minería (Reglamento (Reglamento para la Aprobación de Proyectos de Diseño, Construcción,
Operación y Cierre de los Depósitos de Relaves), se piden los siguientes controles.
‒ Arenas de relave con un porcentaje de finos (0,08 [mm]) menor o igual al 20%.
‒ Tasa de compactación señalada por el ensayo de Proctor normal o modificado.
‒ Criterios para definir el N ° y la ubicación de los puntos de control de compactación no se

definen.
‒ La entrega trimestral de los informes de control de compactación.
Sin embargo, en la práctica lo que se utiliza son los siguientes criterios ingenieriles.
‒ Proctor normal o modificado a alcanzar en la compactación mayor a un 95%.
‒ Densidad relativa entre 60% y 70%.
‒ Controles topográficos para ver espesor de capa y geometría del talud.
‒ Ensayos de penetración dinámica (PANDA).
1.6.1. Control de compactación en tranques de arenas de relave mediante el

penetrómetro ligero PANDA.
Una de las herramientas que se ha adaptado bastante bien a la caracterización del suelo en el
último tiempo tiene relación a la familia de los penetrómetros dinámicos, más específicamente
penetrómetros ligeros o de energía variable, debido a su alcance en las profundidades 6 [m] y
su maniobrabilidad en terreno. Si bien en un principio se utilizó para identificación de terrenos se
adaptó posteriormente para su control de compactación. Es por ello que representa una gran
alternativa para el control de compactación en tranques de relaves debido a las características
que posee.
Principio teórico del control de compactación mediante penetrómetros ligeros.El

principio teórico para el control de compactación con penetrómetros dinámicos se fundamenta
en la relación existente entre la resistencia de punta (qd) y la densidad seca (d) de suelo, la
que depende principalmente de la naturaleza del material y de su estado hídrico. Lo anterior
se fundamenta en hipótesis derivadas de la micromecánica de medios granulares. Ellas
postulan que la caracterización del comportamiento mecánico de un material se realiza a partir
de la descripción de su microestructura [11], para lo cual es posible identificar tres familias de
datos: naturaleza del grano, interacción entre dos granos y fluidos, conjunto de muchos granos.
Considerando lo anterior, fue propuesta la siguiente clasificación de los parámetros empleados
habitualmente en geotecnia, con el objetivo de describir los suelos según una lógica de la
micromecánica [12], [13].
– Parámetros de naturaleza. Mineralogía, granulometría y forma de grano.

– Parámetros de conjunto de partículas. Densidad y anisotropía.
– Parámetros de interacción. Porcentaje de unidad y plasticidad de las partículas finas.
Si el parámetro a medir in situ es la resistencia a la penetración, bajo la hipótesis de

reversibilidad (granos irrompibles, resistentes al desgaste y reversibilidad de la succión) y
conociendo los parámetros de naturaleza e interacción, es posible determinar in situ la

densidad seca (gd) en profundidad del material,Chaigneau[10] ha establecido la relación entre
qd y d, empleando el penetrómetro PANDA. En la fórmula 1.1. se presenta la relación antes
mencionada.
d   A  lnqd    B (1.1.)
Las constantes A y B son valores que dependen del tipo de suelo y de su estado hídrico
[14][15]. Demostraron que el penetrograma resultante de un ensayo de penetración realizado
en un medio granular homogéneo, de densidad y contenido de agua constante, puede ser
esquematizado en el espacio log(qd)-z a partir de dos rectas que definen tres parámetros
característicos; la resistencia de punta a nivel superficial (qdo), la profundidad critica (zc) más
allá de la cual la resistencia de punta es constante y la resistencia de punta constante (qd). En
la Figura 1.24., se muestran una señal penometrica bruta obtenida en el espacio log (qd)-z.
Figura 1.24. A) Señal penetrométrica bruta obtenida en el espacio log (qd)-z para un medio
granular homogéneo. B) Parámetros característicos de la curva de referencia.
Fuente. Chaigneau, (2001).
– Penetrómetro ligero Panda. Para realizar el control de compactación en tranques de relave,

determinando el grado de compactación y el espesor final de capa compactado y considerando
la variabilidad que presenta in situ este tipo de materiales, una alternativa posible de emplear
es el penetrómetro dinámico ligero PANDA (Penetrometre Autonome Numerique Dynamique
Assiste). Comparado con las herramientas habitualmente empleadas para el control de
compactación (método cono de arena y densímetro nuclear), el penetrómetro PANDA presenta
una serie de ventajas, tales como: rapidez en la ejecución del ensayo, fácil repetitividad,

versatilidad del equipo, fiabilidad de los resultados obtenidos, sustentabilidad ambiental, bajos
costos y facilidad de transporte lo que permite analizar lugares de difícil acceso como los
taludes de tranques de relaves [16].
El equipo PANDA es un penetrómetro dinámico ligero a energía variable cuyo principio básico
de funcionamiento consiste en hincar en el suelo, mediante el impacto de un martillo de masa
estándar (2,0 [kg]), un tren de barras de 14,0 [mm] de diámetro provisto en su extremo de
punta cónica metálica de 2,0 [cm2 ] de diámetro.
El impacto de martillo genera una onda de esfuerzo en la cabeza de golpeo, la cual se propaga
hacia la punta del penetrómetro generando la penetración. Para cada golpe de martillo es
registrada de manera continua la penetración alcanzada (e) y la resistencia de punta del suelo
(qd), se obtiene con la adaptación de la fórmula 1.2. la Holandesa.
1 E M
qd   
Ac e M  P (1.2.)
Donde M es la masa del martillo empleado, E la energía aplicada, e la penetración de la punta

en el suelo por cada golpe de martillo, Ac es el área de la punta y P la masa del tren de barras
y de la cabeza de golpe. Una vez finalizado el ensayo, es posible observar de manera
inmediata la señal penetrometrica resultante y las capas de suelo detectadas, en términos de
la resistencia qd. La profundidad máxima de penetración en modalidades control de
compactación es del orden de 1,5 [m]. Los valores máximos para la resistencia de punta qd
varían entre 20,0 a 30,0 [MPa] [17].
Principalmente el control de compactación mediante la utilización del penetrómetro ligero

Panda, se basa en la comparación, de señales de resistencia por punta obtenidas en terreno
y con curvas obtenidas en laboratorio teniendo en consideración calibración del instrumento,
grado de compactación y estado hídrico. A continuación se muestran dos figuras
correspondientes a las zonas establecidas para el control de compactación y señales
penetrométricas.

Figura 1.25. A) Principio para el control de compactación empleando el penetrómetro PANDA.
B) Ejemplo de control de compactación a partir de un penetrograma in-situ.
Fuente. Sol Solution, (2009).
1.7. HIPÓTESIS Y METODOS DE MODELACIÓN.
Para poder generar una representación de la estructura interna de un tranque de relave mediante
una gestión de parámetros de la zona de estudio, es necesario realizar una de modelación de
dichas variables, donde se conocerá la homogeneidad en la masa de suelo o posibles zonas
débiles de esta. Se trabaja con la entrada de parámetros de control de compactación modelando
así espesores de capas, ensayos de toma de muestras y parámetros de penetración [qd], estos
últimos serán la base para clasificar la compacidad, comportamiento y posibles riesgo de fallas
en el tranque. La modelación tiene por hipótesis el comportamiento homogéneo de las
características del tranque, la cual es demostrada por un análisis estadístico de adaptación de
una ley de distribución teórica [18], además que se puede reconstruir la estructura interna de un
tranque de relave con la interpolación de características debidamente referenciadas en el espacio
del tranque.
Se va a entender modelación de un sólido o masa de suelo como verdadero proceso de

cuadriculación de dimensiones, que se utiliza para crear tres dimensiones mediante nodos los
cuales forman una matriz regularmente espaciada. Como serán conocidas las coordenadas del
emplazamiento espacial de los ensayo realizados (X, Y, y Z), como la profundidad de los mismos.

La modelación será bajo la base de interpolación espacial, que mediante el ingreso de datos de
ensayos de campo podrá generar una representación de la masa de suelo en estudio.
El programa utilizado para generar realizar la modelación será el Rockwork que mediante datos
de entrada, levantados por ensayos de campo y utilizando diferentes métodos de interpolación
espacial, podrá representar de mejor manera la estructura interna del tranque de relave.
1.7.1. Métodos de interpolación espacial software ROCKWORK.
Hay varios métodos ofrecidos para llevar a cabo la interpolación espacial, cada uno opera de
manera diferente, y cada uno tiene fortalezas y debilidades. Se describirán a continuación los
métodos de interpolación que podrán ser adoptados para este trabajo.
– Punto más cercano. El método de modelación de sólidos más básico, en la que el valor de
un punto del nodo se establece para que sea igual al valor del punto de datos más cercano.
Es fiel a un control puntos, pero pueden generar modelos con cambios bruscos. El Figura
1.26., se muestra una esquematización de la modelación al punto más cercano.
Figura 1.26. Diagrama de modelación del sólido punto más cercano.

Fuente. Rockworks, (2014).
– Distancia al punto. Este método asigna a cada nodo de modelo sólido de un valor igual a la
Distancia al punto de control más cercano, generando así un modelo basado en la confianza
distancia. La distancia se registra en su X, Y, Z unidades. En la Figura 1.27., se muestra una
esquematización de la modelación distancia al punto.
– Distancia inversa. Este método se asigna un valor de nodo basado en el promedio ponderado
de puntos de datos cercanos o vecinos, o bien todos los puntos o los direccionalmente ubicado,
usando exponentes de ponderación fija o variable.
– Isotrópico. El programa utilizará todos los puntos de datos disponibles al calcular un valor del

nodo dado, útil cuando se modelan los datos distribuidos de manera uniforme En los no-
entornos estratiformes. Inconvenientes: alturas se concentran alrededor de los pozos de
sondeos. Ponderación exponente = 2.
Figura 1.27. Diagrama de modelación del sólido distancia al punto.

Fuente.Rockworks, (2014).
– Anisotrópico. En lugar de utilizar todos los puntos de control disponibles para la distancia
inversamodelado, el programa buscará el punto más cercano en cada sector de noventa
grados alrededor del nodo, útil para el modelado de datos en función del sondeo realizado en
depósitos estratiformes. Tiende a conectar hasta altas y bajas a la misma altura. Ponderación
ex = 2.
– Polarización horizontal. Método que funciona similar al método inverso de la distancia,

excepto que el usuario puede definir una distancia vertical desde cada nodo de la cuadricula
más allá del cual los puntos ya no se usa en el cálculo del valor de nodo.
Luego de saber cuál es el funcionamiento teórico de cada algoritmo para la modelación y

representación espacial, se adopta el que represente de mejor manera la estructura interna del
tranque, en consideración con la variabilidad de las características del tipo del material
depositado.
1.7.2. Herramientas de visualización del modelo espacial.
El Rockwork como programa de modelación de variables reconstruyendo una estructura espacial,

contiene varias formas de visualización del modelo, para poder entregar al usuario una correcta
y variada gama de opciones a la hora de tomar un juicio del modelo entregado.
– Modelo en tres dimensiones (3D). Se representa la interpolación en forma de paralelepípedo,

indicando las dimensiones predispuestas en la base de datos del proyecto de modelación.

Además se indica las tonalidades definidas a las clasificaciones cuantitativas de datos de
ingreso. En la Figura 1.28., se muestra la modelación en tres dimensiones.
Figura 1.28. Modelación tres dimensiones (3D).

– Visualización en 2D de cortes al modelo espacial. Para poder identificar de forma correcta

las características internas del sólido a modelar, se pueden efectuar cortes en cualquier
dirección del modelo visualizando de buena manera el espacio modelado. En la Figura 1.29.,
se muestra la visualización del modelo en dos dimensiones.
Figura 1.29. Perfil en dos dimensiones (2D) sobre el modelo creado.

– Visualización en modo de Fence o rejillas al modelo espacial. La visualización de modo

Fence o rejillas, realiza en el modelo múltiples cortes identificando y comparando de mejor
manera diferentes puntos de la estructura modelada ya sea en su interior como en el exterior.

Figura 1.30. Visualización en modo Fence al modelo creado.

PENETRACIÓN (qd) DEL TRANQUE EN ESTUDIO.
PENETRACIÓN (qd) DEL TRANQUE EN ESTUDIO.
Cabe señalar que se puede definir un tranque de relave con un comportamiento homogéneo por
sus parámetros naturaleza y de estado, debido a que la procedencia de sus arenas de relaves
que corresponden a una familia de origen mineralógico y el método de construcción utilizado. No
obstante siempre existirá una variabilidad de sus parámetros los cuales deben ser estudiados
para cualificar la manera en como varían sus características.
Para eso, será necesario realizar análisis estadístico de los datos recopilados mediante varias
campañas de ensayos realizadas del tranque de estudio, en donde se busca adquirir parámetros
estadísticos los definirán de mejor la manera las variables estudiadas.
El objetivo es poder adaptar leyes estadísticas de distribución que describen lo mejor posible el
comportamiento aleatorio de las variables. Así se puede utilizar estas leyes para formular una
hipótesis de trabajo en donde indica la homogeneidad del tranque de relave en estudio. Estos
análisis serán usados posteriormente como base para el modelo de representación de
parámetros de control de compactación en donde se sustentara las características del tranque de
relave como homogéneo.
Se enfoca el estudio hacia la variación de las características físicas, específicamente las de

resistencia por punta [qd], las que fueron entregadas mediante ensayos de penetración dinámica.
Se recopila información de campañas de ensayos del tranque de relave en estudio, acotado en

un periodo comprendido del año 2009 al 2014, en donde se analizan señales penetrómetricas del
resultado de los ensayos de penetración ligera. El modo a analizar se hará mediante herramientas
clásicas estadísticas las cuales seguirán el siguiente orden de trabajo.
– Obtención de los parámetros estadísticos descriptivos de los datos, generando una batería
amplia de información de entrada al modelo estadístico.
– Construcción de los histogramas de frecuencia.
– Obtención de los parámetros de forma de los datos, momentos estadísticos de orden tres
(Asimetría) y de orden cuatro (Aplanamiento).
– Adaptación de una ley de distribución a partir de pruebas estadísticas (Kolmogorov-Smirnov).
2.1. INFORMACIÓN GENERAL DEL TRANQUE DE RELAVE EN ESTUDIO (TRANQUE M-1).
Para poder distinguir y dimensionar las características globales del tranque en estudio, se realiza
la siguiente tabla.

Tabla 2.1. Información general Tranque de Relave M-1.
INFORMACIÓN GENERAL TRANQUE DE RELAVES M-1

Tranque de relaves Minera SAP - Tranque M-1
Tipo Agua abajo
Almacenamiento del material Almacenamiento hidráulico
Construcción Por capas de espesor igual a 50 - 60 [cm]
Compactación Mecánica mediante de rodillo liso vibratorio de 10 [t]
Densidad máxima seca de las arenas

18,1 [kN/m3]
de relaves
Grado de compactación apuntado 95% OPN
Coronamiento Anchura mínimo: 10,0 [m]
Desquite mínimo 5 [m]
Altura media 100 [m]
Talud agua abajo 1:4,3 (V:H)
Talud agua arriba 1:2,0 (V:H)
Sistema de drenaje Alfombra que drena
Impermeabilización del talud agua
Si
arriba
2.2. RECOPILACIÓN DE INFORMACIÓN.
El Análisis de datos entregados se realiza a partir de controles de calidad a los depósitos arenas
de relaves, efectuado mediante ensayos de penetración dinámica (Panda), toma de muestras in
situ y ensayos de laboratorio, en el tranque de relave M-1, minera SAP. Este control se realiza
para verificar el cumplimiento de los estándares de calidad definidos por Minera SAP para la
operación del depósito en donde se acotan los siguientes parámetros de control.
– Arenas con porcentaje de finos (80[m]) menor o igual a 20%.

– Control de compactación mediante ensayo Proctor normal o modificado.
– Ensayos de penetración dinámica.
Para el análisis de esta batería de datos se realiza la siguiente metodología.
– Definición de sectores a analizar, modos de situación de ensayos dentro del campo espacial
y georeferenciación de ellos.
– Procesamiento de los penetrogramas obtenidos en cada una de las prospecciones definidas,

mediante el análisis de la señal penetrómetricas obtenida en cada punto de prospección
identificando resistencia por punta (qd).
– Realización de alisadura de datos penetrometricos cada 0,025 [m], asociando resistencia por
punta.
2.2.1. Descripción del área de estudio.
La campaña de trabajos de terreno fue desarrollada en los periodos de años 2009 – 2014, en los
sectores denominados S-1, S-2, S-3del tranque de relave M-1 de Minera SAP, en donde se
realizaron ensayos de penetración dinámica. En la Figura 2.1., se caracteriza la sectorización con
su respectiva nomenclatura del muro resistente.
Figura 2.1. Sectorización del área de estudio.

2.2.2. Emplazamiento de prospecciones.
Para el emplazamiento de los puntos de prospección definidos, se consideró como referencia la

ubicación de un sistema de estacado, posicionado en forma de cuadrícula permanente en el muro
resistente del tranque M-1.
El total de prospecciones mediante ensayos de penetración ligera (PANDA) llega a un total de

ciento treinta y cinco, entregando características de resistencia (qd) del tranque de relave,
generando una base de datos bastante grande para poder analizar la variación mediante
evaluación probabilística de características de los parámetros de control de compactación.
En la figura siguiente se presenta de manera esquemática la ubicación de los ensayos PANDA

realizados, empleando la modalidad de reconocimiento del terreno.
Figura 2.2. Ubicación ensayos PANDA.

2.3. ANÁLISIS ESTADÍSTICO DE LOS DATOS RECOPILADOS.
2.3.1. Resumen estadístico.
El conjunto de las señales penetrómetricas obtenidas en el momento de campañas

correspondientes a los periodos 2009 – 2014, arrojaron los siguientes datos estadísticos del
tranque M-1, representados en la Tabla 2.2.
Tabla 2.2. Resumen estadístico valores logarítmicos del parámetro qd Tranque M-1.
qd [MPa]
Tranque Desviació Ley de

Periodo Coef. de Mínimu Máximu
Datos Media n distribució
variación m m
estándar n
Tranque
M-1– 2009 a Log-
Minera 36910 0,68 0,188 27,45 % 0 2,93
2014 normal
SAP
El valor del coeficiente de variación obtenido no es un valor alto, teniendo en cuenta que es
posible contemplar dispersión en este parámetro (qd), por lo que las arenas de relaves del tranque

presentan una condición de homogeneidad con respecto al parámetro de resistencia a la
penetración.
Adicionalmente se realizan estos análisis estadísticos a los valores de la resistencia a la

penetración (qd) en forma global para cada sector. Los resultados obtenidos se muestran en la
Tabla 2.3.
Tabla 2.3. Resumen estadístico por sectores Tranque M-1.
DATOS S-1 S-2 S-3

Recuento 5078 11245 4654
Promedio 5,13179 5,0761 4,73357
Desv. Estándar 2,24286 1,82012 1,72554
Coef. Variación 43,7052% 35,8567% 36,4532%
Mínimo 0,74931 0,283415 0,51472
Máximo 15,4461 16,1675 11,1889
Al igual que los resultados para la totalidad del tranque, se pueden apreciar coeficientes de
variación inferiores al 50%, lo que evidencia la homogeneidad de las arenas de la zona de estudio.
Existen valores mínimos de resistencia cercanos a cero, sin embargo estos valores son en menor
frecuencia y pueden ser despreciables por encontrarse en su mayoría próximos a la superficie.
2.3.2. Distribución estadística del parámetro de resistencia a la penetración (qd).
La distribución de los valores de resistencia a la penetración en el tranque en estudio no pueden

ser conocidos en toda su extensión sin embargo es posible utilizar leyes de distribución.
Se mostrara la distribución de los datos obtenidos mediante histogramas de frecuencia. Como se

mencionó en el marco teórico se ha demostrado que el logaritmo de los valores puede seguir una
ley de distribución normal.
Para esto el parámetro en estudio será analizado variables aleatorias continuas, esperanza
matemática y varianza. En la Figura 2.3., se presenta el histograma de frecuencia de los datos
recopilados mediante ensayos de penetración.
El histograma obtenido muestra una forma muy regular, simétrica a un solo comportamiento y
una reducida dispersión, comprobando el coeficiente de variación obtenido anteriormente. Esto
representa un material homogéneo producto de la naturaleza del material (de un solo yacimiento)
y de la depositación, ya que se utiliza un mismo proceso de compactación.

Figura 2.3. Histograma de frecuencia resistencia a la penetración, tranque M-1.
Como se observa en el gráfico, existe una mayor concentración de datos en el centro de éste
debido a la solicitación de calidad con respecto a la resistencia que se debe obtener en el proceso
de compactación.
2.3.3. Leyes de probabilidad que se ajustan a la distribución del parámetro de resistencia

a la penetración (qd).
En la práctica no se utilizan leyes experimentales de distribución, debido a esto se modela la

frecuencia de los datos arrojados, adoptando una ley de distribución teórica que mejor represente
el comportamiento del tranque.
Figura 2.4.Histograma ajustamiento de la distribución del parámetro resistencia a la

penetración, tranque M-1.

Para esto se ejecuta la prueba de bondad de ajuste a través del test de Kolmogorov- Smirnov. La
distribución de ajuste adoptada, que representa de manera global la forma general del histograma
se presenta a continuación en la Figura 2.4.
Además se presenta la función de densidad de la distribución entregada anteriormente, el cual

es esencialmente un histograma suavizado de la muestra, mostrando la concentración de datos
para cada valor de resistencia. En la Figura 2.5., se presenta la función densidad del tranque
estudiado.
Figura 2.5. Distribución de la resistencia a la penetración tranque M-1.

En figura anterior deja en evidencia la condición de homogeneidad a partir de la reducida

dispersión que presentan los datos graficados, demostrando una vez el origen único de las arenas
del tranque y del proceso de depositación.
2.3.4. Conclusiones de variabilidad Tranque M-1.
Luego del análisis estadístico de los datos del parámetro de resistencia a la penetración, mediante
la metodología presentada, es posible señalar las siguientes conclusiones.
– Se detecta una única familia de datos, la cual fue representada en el histograma de frecuencia.
Esto producto de la procedencia del material y la forma en que ha sido depositado. El
coeficiente de variación de los datos analizados fue de un 27% aproximadamente, el cual es
un adecuado valor para este tipo de parámetro.
– El logaritmo de los valores del parámetro de resistencia a la penetración, muestra un ajuste

global a la ley de distribución normal, describiendo la variabilidad del parámetro.

– Debido a este análisis, es posible demostrar la homogeneidad que presentan las arenas de
relaves de las prospecciones realizadas.

LICUEFACCIÓN EN FUNCIÓN DEL PARÁMETRO DE RESISTENCIA A LA
PENETRACIÓN [qd].
LICUEFACCIÓN EN FUNCIÓN DEL PARÁMETRO RESISTENCIA A LA PENETRACIÓN
[qd].
La estimación de la compacidad, comportamiento mecánico y riesgo de licuefacción en la zona

de estudio, se realiza con la finalidad de verificar el control de calidad en profundidad del tranque
de arenas de relaves, depositadas y compactadas. Esto se ejecuta para el cumplimiento de los
estándares de calidad exigidos según normativa y los de la propia empresa minera.
Para el desarrollo de este análisis se efectuará a partir de la siguiente metodología.
– Definir las áreas de emplazamiento en que se divide el tranque M-1 para el posterior análisis.
– Ubicación de los ensayos de penetración ligera (PANDA) en cada uno de los sectores definidos
en el tranque, con sus respectivas fechas de ejecución.
– Determinación de la variabilidad a través de análisis estadísticos para los parámetros

resistencia por punta (qd) y espesores de capa, en forma global usando todos los datos
obtenidos y en forma local por cada prospección.
– Determinación de la relación entre la resistencia por punta normalizada y el estado de

compactación, a través de la densidad relativa (DR).
– Clasificación y análisis del grado de compactación, comportamiento mecánico y riesgo de

licuefacción.
El tranque M-1se encuentra dividido en tres sectores (S-1, S-2, S-3), para efectos de logística y
análisis de información de depositación y control.
Para el desarrollo de la estimación de la compacidad, comportamiento mecánico y riesgo de

licuefacción se considerarán los datos de las campañas en terreno del período de años entre
2010 y 2013. La cantidad de ensayos realizados en cada período para cada sector se muestra
en la Tabla 3.1.
Tabla 3.1. Cantidad de ensayos por sectores, Tranque M-1.
Sectores
Año
S-1 S-2 S-3
2010 6
2011( Marzo) 3 6 3
2011( Junio) 3 6 3
2011 ( Noviembre) 3 6 3
2012 ( Febrero) 3 6 3

Sectores
Año
S-1 S-2 S-3
2012 ( Junio) 3 6 3
2012 ( Agosto) 3 6 2
2012 ( Noviembre) 3 6 3
2013 ( Noviembre) 6 6 6
Total por sector 27 54 26
Total tranque 107
Además se determinaron tres perfiles longitudinales que definirán los ensayos en cada sector y
de esta forma poder representar las características del tranque. El perfil N°1 corresponde a todas
las prospecciones realizadas en el sector 1; el perfil N°2 corresponde a las prospecciones del
sector S-2 y el perfil N°3 queda determinado por las prospecciones del sector 3. La distribución
de los perfiles en el tranque se presenta a continuación.
Figura 3.1. Emplazamiento prospecciones.

Tabla 3.2. Ubicación de las prospecciones según sector.

PROSPECCIÓN SECTOR
P1 S-1
P2 S-1
P3 S-1
P4 S-1
P5 S-1
P6 S-1
P7 S-2
P8 S-2
P9 S-2
P10 S-2
P11 S-2
P12 S-2
P13 S-3
P14 S-3
P15 S-3
P16 S-3
P17 S-3
P18 S-3
3.1. VARIACIÓN DE LA RESISTENCIA POR PUNTA (qd).
En el capítulo N°2 se presentó el análisis de la variabilidad de la resistencia por punta a nivel

global, con la totalidad de datos obtenidos a partir de los ensayos de penetración ligera mediante
el equipo PANDA. Ahora se realizará este análisis en dos niveles, de forma global por perfiles y
de manera local a través de cada prospección.
Tabla 3.3. Análisis estadístico global del parámetro resistencia por punta, por perfil.
PARÁMETROS ESTADÍSTICOS [MPa]

PERFIL
MEDIA DESV. ESTÁNDAR COEF. VARIACIÓN (%)
Perfil 1 5,13 2,24 43,71
Perfil 2 5,08 1,82 35,86
Perfil 3 4,73 1,73 36,45
Tabla 3.4. Análisis estadístico local del parámetro resistencia por punta por prospección, perfil 2

año 2010.

PERFIL 2 (2010)
P1 3,45 1,09 32
P2 4,21 1,16 28
P3 5,67 1,87 33
P4 4,54 1,02 23
P5 4,5 1,28 28
P6 5,96 1,6 27
año 2011 (marzo).

PERFIL 1 (2011-MARZO)
P1 3,50 1,09 31
P2 4,6 1,5 33
P3 3,66 1,36 37
año 2011 (marzo).

P1 4,42 1,74 39
P2 3,4 1,15 34
P3 5,05 1,79 35
P4 4,15 1,47 35
P5 4,59 1,65 36
P6 2,92 1,2 41

año 2011 (marzo).

P1 2,12 0,97 46
P2 4,5 1,59 35
P3 4,4 1,53 35
año 2011 (junio).

PERFIL 1 (2011-JUNIO)
P1 3,86 1,76 46
P2 5,95 1,36 23
P3 5,39 1,53 28
año 2011 (junio).

P1 4,67 1,56 33
P2 4,85 1,44 30
P3 5,25 1,51 29
P4 5,3 1,48 28
P5 3,89 1,88 48
P6 4,93 1,67 34
Tabla 3.10. Análisis estadístico local del parámetro resistencia por punta por prospección, perfil
3 año 2011 (junio).

P1 4,12 1,2 29
P2 4,6 1,71 37
P3 4,66 1,92 41

1 año 2011 (noviembre).

PERFIL 1 (2011-NOV)
P1 5,60 1,98 35
P2 6,87 2,47 36
P3 7,39 1,91 26

PERFIL 2 (2011-NOV)
P1 3,58 1,53 43
P2 5,31 2,35 44
P3 3,39 1,56 46
P4 5,6 1,98 35
P5 6,87 2,47 36
P6 7,39 1,91 26

PERFIL 3 (2011-NOV)
P1 3,58 1,53 43
P2 5,31 2,35 44
P3 3,39 1,56 46
1 año 2012 (febrero).

PERFIL 1 (2012-FEBR)
P1 6,49 1,83 28
P2 5,05 2,23 44
P3 3,56 1,38 39


P1 4,28 1,23 29
P2 5,67 1,64 29
P3 5,15 2,12 41
P4 4,5 1,62 36
P5 4,4 1,66 38
P6 4,13 1,56 38

P1 4,10 1,59 39
P2 5,18 2,14 41
P3 4,25 1,66 39

P1 3,74 1,74 46
P2 2,62 1,11 42
P3 5,14 2,45 48


P1 5,33 1,65 31
P2 4,32 1,74 40
P3 4,3 1,19 28
P4 4,84 1,55 32
P5 4,49 1,19 27
P6 6,37 2,3 36

P1 3,43 1,17 34
P2 5,47 1,9 35
P3 3,43 1,17 34
Tabla3.20. Análisis estadístico local del parámetro resistencia por punta por prospección, perfil 1
año 2012 (agosto).

PERFIL 1 (2012-AGOSTO)
P1 6,01 2,49 41
P2 3,96 1,84 46
P3 4,2 1,8 43
Tabla3.21. Análisis estadístico local del parámetro resistencia por punta por prospección, perfil 2
año 2012 (agosto).

P1 5,12 2,63 51
P2 3,7 1,24 34
P3 3,21 1,36 42
P4 4,86 2,37 49

P5 5,14 1,88 37
P6 4,93 2,4 49
3 año 2012 (agosto).
PERFIL 3 (2012- PARÁMETROS ESTADÍSTICOS [MPa]

AGOSTO) MEDIA DESV. ESTÁNDAR COEF. VARIACIÓN (%)
P1 4,85 1,85 38
P2 3,63 1,33 37

PERFIL 1 (2012-NOV)
P1 8,07 3,57 44
P2 6,02 1,68 28
P3 3,1 0,91 29

PERFIL 2 (2012-NOV)
P1 4,14 1,12 27
P2 4,55 1,6 35
P3 5,47 1,89 35
P4 5 2,22 45
P5 5,02 1,65 33
P6 5,25 1,59 30


PERFIL 3 (2012-NOV)
P1 4,45 1,78 40
P2 4,21 1,12 27
P3 5,31 2,18 41

PERFIL 1 (2013-NOV)
P1 5,20 1,1 21
P2 5,7 1,2 21
P3 6,8 1,6 23
P4 6,5 2,7 41
P5 6,4 2,5 40
P6 4,2 1,6 37

PERFIL 2 (2013-NOV)
P1 4,00 1,8 44
P2 3,7 1,2 34
P3 4,8 1,5 32
P4 5,5 2,6 49
P5 4,7 2,6 54
P6 4,8 2,8 58


PERFIL 3 (2013-NOV)
MEDIA DESV. ESTÁNDAR COEF. VARIACIÓN [%]
P1 5,10 1,8 35
P2 4,9 1,2 24
P3 5 1,5 30
P4 6 2,4 41
P5 5,8 2,2 38
P6 6 3,7 61
Los resultados obtenidos muestran coeficientes de variación entre 35% y 43% en forma global,
lo cual también se ve reflejado en las prospecciones locales. Estos valores generalmente son
inferiores al 50% por lo tanto representan un estado de homogeneidad en el interior del tranque,
y respaldan el análisis de variabilidad realizado en el capítulo 2.
Con respeto a los resultados de la media de la resistencia por punta, se tienen valores cercanos
a 5 [MPa] en el análisis global y similares en el análisis local, con resultados bajo 3 [MPa] en
casos muy puntuales, los cuales se podrían asociar a las capas superiores al momento de hacer
la prospección, las cuales por factores externos como el transporte de la maquinaria y la erosión,
no son representativas para el análisis. Por el contrario, los valores sobre 3[MPa] representan un
nivel de compacidad de compacto a denso, el cual es producto del envejecimiento que sufren las
arenas de relaves por efecto de pequeños sismo que provocan densificación y por la depositación
consecutiva de las capas de arenas de relaves.
3.2. RELACIÓN ENTRE LA RESISTENCIA POR PUNTA NORMALIZADA (qdN1) Y EL GRADO

DE COMPACTACIÓN.
El nivel de compactación que se expresa en función de los intervalos de los valores de densidad
relativa (DR), ha sido correlacionado a partir de la resistencia por punta normalizada (qdN1),
empleando la relación propuesta por Villavicencio [18], calibrada para arenas de relaves del
Tranque M-1, a partir de la fórmula 3.1., que se muestra a continuación.
𝐷𝑅(%) = 𝐴 ∙ ln(𝑞𝑑𝑁1 ) − 𝐵 (3.1.)
Las constantes A y B dependen del origen mineralógico, de las características físicas en relación
a la distribución granulométrica, y del estado hídrico que presentan las arenas de relaves in situ.
Cabe señalar que posteriormente esta relación se ha ajustado en forma más exacta, a través de

datos de las propiedades del tranque M-1, las cuales fueron las siguientes.
‒ Arena Limosa SM.

‒ D50 = 0,19[mm].
‒ % menor a la malla #200= 14-19.
‒ Gs: 2,6-2,7.
‒ DMCS= 1,68[t/m3].
‒ W óptima= 17%.
‒ Estado Hídrico= seco.
A partir de estas propiedades mostradas, se actualizan las constantes A y B para esta zona de
estudio, obteniéndose la fórmula 3.2., que se presenta a continuación.
DR(%) = 19,25 ∙ ln(qdN1 ) − 18,48 (3.2.)
El parámetro de qdN1, se obtiene como se muestra en la fórmula 3.3., que se presenta a

continuación.
qdN1 = CN ∗ qd (3.3.)
3.2.1. Normalización del parámetro resistencia (qd).
El parámetro de resistencia a la penetración normalizado qdn1depende del esfuerzo vertical

efectivo (σ’) y la densidad relativa (DR), entre otros factores. La normalización se realiza con el
factor CN que en un principio dependía de la relación entre el esfuerzo efectivo (σ’) y una
atmosfera [atm], considerando una DR constante (DR= 54 %). Investigaciones han demostrado
que es necesario realizar iteraciones en función a la relación de presión y esfuerzo efectivo y
además a la densidad relativa, que al aumentar su valor en profundidad, se tiene que considerar
como parámetro variable según lo descrito por Boulanger, R. W. Idriss, I. M. [19]. Es por ello que
fue propuesta la fórmula 3.4., para obtener el factor CN.
Pa 0,784−0,521DR
CN = ( σ’
) (3.4.)
De esta manera es posible obtener la densidad relativa, en función del parámetro de resistencia
normaliza (qdN1), como se muestra en la fórmula 3.5.
DR = 0,478( qdN1)0,264 − 1,063 (3.5.)
Para efectos más prácticos se llegó a la fórmula 3.6., que se presenta a continuación.
Pa 1,338−0,249(qdN1) 0,264
CN = ( σ’
) (3.6.)
Para la resolución del parámetro CN es necesario realizar una iteración, que se puede efectuar
con cualquier software. Para efectos de esta tesis se realizó con el software Microsoft Excel 2013,

aplicando la herramienta referencia circular, dejando como variables la densidad relativa y el
parámetro qdN1.
3.3. CLASIFICACIÓN Y ANÁLISIS DEL GRADO DE COMPACTACION, COMPORTAMIENTO

MECÁNICO Y RIESGO DE LICUEFACCIÓN.
Para estimar el grado de compactación en arenas de relaves, se basa en una clasificación por
intervalos en función de la densidad relativa, la cual es calculada a partir del parámetro q d1,
mediante la normalización presentada en el apartado previo. Esta clasificación fue propuesta por
Espinace [20] para tranques de relaves chilenos, la cual se presenta en la tabla de la Figura 3.2.
DENSIDAD
ESTADO DE COMPORTAMIENTOMEC POTENCIAL DE
qdN1 RELATIVA
COMPACIDAD ÁNICO LICUACIÓN
[%]
<7,5 < 20 Muy Suelto Contractivo Muy Alto
7,5 - 27 20 - 45 Suelto Contractivo Alto
27 - 35 45 - 50 Compacto Contractivo Equilibrio
Compacto a
35 - 75 50 - 65 Límite Bajo
Denso
75 - 215 65 - 85 Denso Dilatante Muy Bajo
>215 85 - 100 Muy Denso Dilatante Nulo
Figura 3.2. Clasificación del nivel de compacidad y estimación del potencial de licuación, en
función de qdN1 y (N1)60.
Además, R. Verdugo [21] demostró que arenas de relaves con un porcentaje de finos (menor a
80 [μm]) entre un 10% y 25%, como el caso de las arenas del muro del Tranque M-1, un grado
de compactación de 95% OPN es equivalente a una densidad relativa (DR) de entre un 60% y
65%. Esto corresponde al límite de un comportamiento mecánico dilatante o contractivo [22].
Por lo tanto, es necesario definir el significado de cada uno de estos estados, en donde el
comportamiento podrá ser contractivo o dilatante, dependiendo de su presión de confinamiento y
de la densidad del material. Esto podrá desencadenar en el fenómeno de licuefacción, el cual se
subdividirá en Licuación verdadera o falla de flujo y Licuación o movilidad cíclica. El primer caso
será cuando la masa de suelo pierde su resistencia última no drenada de manera repentina y sus
tensiones solicitantes son mayores. El fenómeno de movilidad cíclica o licuación, ocurrirá cuando
la masa de suelo sea solicitada por cargas cíclicas como por ejemplo sismos [23], provocando un
aumento progresivo de deformaciones y aumento de la presión de poros hasta un 100%

(momento en el cual ocurre la licuación). Se ha demostrado que tanto materiales finos en estado
de compacidad sueltos como densos son susceptibles de sufrir este fenómeno, dependiendo de
la duración y magnitud del evento que lo generó.
En la Figura 3.3., se grafican los eventuales casos de fallas de flujo, en que se detallan cuatro
casos posibles, ya sea en forma drenada o no drenada.
En la figura 3.3., se puede apreciar las gráficas que muestran el esfuerzo al corte solicitante con
respecto a las resistencias máximas en condiciones drenadas o no drenadas. Al considerar
ambos valores podemos determinar un factor de seguridad para cada caso, el cual será mayor o
igual a la unidad si es que el corte solicitante sea inferior a la resistencia máxima. Se considera
el caso no drenado cuando la solicitación ocurre en forma rápida, como un sismo, movilizando la
resistencia no drenada.
Figura 3.3. Casos posibles para comportamiento Dilatante y Contractivo.

Fuente. Verdugo, (1995).
Se proceden a detallar los cuatro casos mostrados en la Figura 3.3.
– Caso A. Se puede observar que la resistencia no drenada es mayor a la drenada y ésta mayor
al corte solicitante. Por lo tanto, se está en presencia de un factor de seguridad mayor a 1. En
caso de perderse la resistencia no drenada producto de un sismo, el factor de seguridad
seguirá siendo superior a 1 y mayor que el caso drenado, por ser la resistencia también mayor.
La falla ocurre cuando se ubica bajo la gráfica, siendo del tipo dilatante.
– Caso B. La diferencia al caso anterior, es que la resistencia máxima drenada es superior a la

no drenada y ésta al corte solicitante. Ambas resistencias siguen siendo superior al corte
solicitante, por lo tanto se tiene un factor de seguridad superior a 1. Sin embargo, si se moviliza
la resistencia máxima no drenada, a través de un efecto rápido, el factor seguirá siendo mayor
a 1, pero inferior que en el caso drenado. Esto se traduce en una menor dilatancia que el caso
anterior, es decir, sigue estando ubicado bajo la curva pero más cerca de ésta.
– Caso C. Al igual que el caso B, la resistencia máxima drenada es mayor a la no drenada y

esta última superior al corte solicitante. Como se nombra previamente, el comportamiento
dependerá de la densidad y del confinamiento del material, en este caso este conjunto se ubica
sobre la gráfica, por lo tanto pertenece a un comportamiento contractivo. El factor de seguridad
sigue siendo superior a la unidad, por lo tanto la falla no ocurre.
– Caso D. Éste caso está definido por una resistencia máxima drenada superior al corte
solicitante, pero este es mayor a la resistencia máxima no drenada. Es por esto, que si el corte
solicitado es alcanzado en condiciones no drenadas es mayor a la resistencia movilizada
drenada y el factor de seguridad es superior a 1. Sin embargo, si se aplica una carga rápida
(sismo), se moviliza la resistencia no drenada, la cual es menor al corte solicitante,
obteniéndose un factor inferior a 1, teniendo como resultado final la ocurrencia de falla de flujo
con grandes deslizamientos de material.
Además, al igual que para el parámetro de resistencia por punta (qd), se realizó un análisis
estadístico de los resultados obtenidos de densidad relativa (DR) en forma global para todos los
datos y por perfil y además un análisis en formal local a partir de cada una de las prospecciones
PANDA, con su respectiva clasificación del comportamiento mecánico asociado.
A continuación se muestra la tabla con los resultados a partir de todos los datos.
Tabla 3.29. Clasificación y análisis del parámetro DR en forma global para Tranque M-1.
COMPORTAMIENTO POTENCIAL DE
COMPACIDAD DR (%) % VALORES
MECÁNICO LICUACIÓN
Muy Suelto <20 0% Contractivo Muy Alto
Suelto 20-45 7% Contractivo Alto
Compactado 45-50 6% Contractivo Equilibrio
Compactado a Denso 50-65 51% Límite Bajo
Denso 65-85 36% Dilatante Muy Bajo
Muy Denso >85 0% Dilatante Nulo
3.3.1. Clasificación y análisis estadístico del parámetro densidad relativa (DR) a nivel
global por perfil.

Se realiza la clasificación de las arenas de relaves a partir del parámetro densidad relativa (DR),
utilizando la tabla propuesta por R. Espinace A., que se mostró previamente. Además se
analizaron los datos, mostrando los porcentajes de valores contenidos en cada intervalo, para
esto relacionarlo con su respectivo comportamiento mecánico y riesgo de licuefacción.
Posteriormente se mostraran los histogramas que grafican estos análisis por perfil.
Tabla 3.30. Clasificación y análisis del parámetro DR perfil N°1 a escala global.
COMPACIDAD DR (%) % VALORES COMPORTAMIENTO POTENCIAL DE


Figura 3.4. Distribución valores densidad relativa para perfil 1 en escala global.



Se puede observar tanto en el análisis global para la totalidad de datos, como para el análisis
global por perfil que la mayor parte de los resultados de densidad relativa se concentran en el
intervalo de compacto a denso, seguido por el intervalo denso, sumando casi un 90% de la
totalidad. Esto representa un comportamiento mecánico asociado, como límite o dilatante, con un
respectivo riesgo de licuefacción bajo o muy bajo. Por el contrario, se obtienen resultados bajo el
10% para arenas de relaves en estado de compacidad suelto, los cuales se pueden asociar a los
espesores superiores que no son representativos como se menciona en el análisis del parámetro
resistencia por punta o por zonas puntuales que presentan una deficiente compactación.
3.3.2. Clasificación del parámetro densidad relativa (DR) a nivel local por prospección.
Al igual que el análisis global, se realizó la misma clasificación planteada por R. Espinace, en
donde se mostrarán ejemplos de prospecciones para cada campaña realizada, entre los años
2010 y 2013.
Tabla 3.33. Clasificación DR prospección n°1, año 2010.
LÍMITE LÍMITE ESPESOR ESTADO DE COMPORTAM. POTEN.

SUPERIOR INFERIOR CAPA COMPACIDAD MECÁNICO. LICUACIÓN.
0,0 0,9 0,9 Denso Dilatante Muy Bajo
2,2 3 0,8 Compacto a Denso Límite Bajo

3 3,5 0,5 Compacto a Denso Límite Bajo
5,1 5,5 0,4 Compacto Contractivo Equilibrio
Tabla 3.34. Clasificación DR prospección n°8 sector 3, Junio 2011.

0,5 0,9 0,4 Compacto Límite Bajo
0,9 1,6 0,6 Compacto a Denso Límite Bajo
5,0 6,0 1 Compacto a Denso Límite Bajo
Tabla 3.35. Clasificación DR prospección n°1 sector 2, marzo 2011.

3,5 4,5 1 Denso Dilatante Muy Bajo
Tabla 3.36. Clasificación DR prospección n°10 sector 2, noviembre 2011.

1,6 2 0,4 Denso Dilatante Muy Bajo
2 2,7 0,7 Denso Dilatante Muy Bajo
Tabla 3.37. Clasificación DR prospección n°3 sector 2, febrero 2012.

Tabla 3.38. Clasificación DR prospección n°2 sector 1, junio 2012.

0,0 1 1 Compacto Contractivo Muy Bajo
1 1,8 0,8 Compacto a Denso Límite Muy Bajo
1,8 2,4 0,6 Compacto a Denso Límite Muy Bajo
5,5 6 0,5 Suelto Contractivo Alto
Tabla 3.39. Clasificación DR prospección n°2 sector 3, agosto 2012.

Tabla 3.40. Clasificación DR prospección n°1 sector 1, noviembre 2012.

0,0 1 1 Denso Dilatante Muy Bajo
4,9 6 1,1 Compacto Contractivo Equilibrio
Tabla 3.41. Clasificación DR prospección n°1 sector Canelo, noviembre 2013.


5,4 6,1 0,7 Compacto Contractivo Bajo
Los resultados obtenidos, en general presentan en mayor magnitud un estado de compacidad

compacto a denso, es decir, un comportamiento mecánico límite. Esto se puede verificar con el
análisis estadístico realizado previamente, en donde se pudo observar en los histogramas que
más del 50% de los datos pertenecen al intervalo de comportamiento límite, es decir, con un
porcentaje DR entre el 50% y 65%.
3.3.3. Análisis estadístico del parámetro densidad relativa (DR) a nivel local por
prospección.
Para complementar el análisis global realizado para cada perfil, se presenta el análisis estadístico
para cada una de las prospecciones realizadas, clasificadas según perfil y fecha de ejecución. A
continuación se presentan tablas resúmenes con valores de media, desviación estándar y el
respectivo coeficiente de variación (CV).
Tabla 3.42. Análisis estadístico nivel local parámetro densidad relativa perfil 2, año 2010.
PARÁMETROS ESTADÍSTICOS DR
PERFIL 2 (2010)
P1 58 10 18
P2 61 8 13
P3 67 7 10
P4 63 7 11
P5 64 5 8
P6 68 7 11
Tabla 3.43. Análisis estadístico nivel local parámetro densidad relativa perfil 1, año 2011 (marzo).
P1 59 7 12
P2 63 7 11
P3 59 11 19

P1 62 9 15
P2 57 9 16
P3 65 10 15
P4 61 10 16
P5 63 11 18
P6 53 11 20
P1 48 10 21
P2 62 7 12
P3 62 10 16
Tabla3.46. Análisis estadístico nivel local parámetro densidad relativa perfil 1, año 2011 (junio).
P1 59 8 14
P2 69 8 12
P3 67 8 12
Tabla 3.47. Análisis estadístico nivel local parámetro densidad relativa perfil 2, año 2011 (junio).
P1 63 10 16
P2 64 9 14
P3 66 8 12
P4 66 9 13
P5 60 12 19
P6 65 8 13

P1 62 8 14
P2 63 8 13
P3 63 12 20
Tabla 3.49. Análisis estadístico nivel local parámetro densidad relativa perfil 1, año 2011
(noviembre).
PERFIL 1 (2011-NOV)
P1 67 11 17
P2 70 8 12
P3 72 7 10
(noviembre).
PERFIL 2 (2011-NOV)
P1 57 10 18
P2 63 11 18
P3 56 9 15
P4 67 11 17
P5 70 8 12
P6 72 7 10
(noviembre).
PERFIL 3 (2011-NOV)

P1 57 10 18
P2 63 11 18
P3 56 9 15
Tabla 3.52. Análisis estadístico nivel local parámetro densidad relativa perfil 1, año 2012 (febrero).
P1 69 8 11
P2 64 8 13
P3 58 7 12
P1 62 7 11
P2 67 7 10
P3 65 8 13
P4 61 9 15
P5 61 9 15
P6 60 9 15
P1 60 7 12
P2 63 10 15
P3 60 8 13
P1 57 11 17
P2 51 10 20
P3 63 9 15

P1 65 8 12
P2 61 6 10
P3 62 8 13
P4 65 7 11
P5 63 8 13
P6 68 8 13
P1 58 9 16
P2 66 9 13
P3 58 9 16
(noviembre).
PERFIL 1 (2012-NOV)
P1 71 13 18
P2 68 9 13
P3 56 9 15
(noviembre).
PERFIL 2 (2012-NOV)
P1 61 5 9
P2 62 9 14

PERFIL 2 (2012-NOV)
P3 66 8 13
P4 63 9 15
P5 63 8 13
P6 65 8 13
(noviembre).
PERFIL 3 (2012-NOV)
P1 60 10 17
P2 61 7 12
P3 64 11 18
(noviembre).
PERFIL 1 (2013-NOV)
P1 65 5 8
P2 66 6 9
P3 70 5 8
P4 67 11 17
P5 67 11 16
P6 70 8 13
(noviembre).
PERFIL 2 (2013-NOV) PARÁMETROS ESTADÍSTICOS DR

P1 58 11 19
P2 57 9 15
P3 62 9 15
P4 64 11 17
P5 59 14 24
P6 60 12 20
(noviembre).
PERFIL 3 (2013-NOV)
P1 64 8 13
P2 63 6 9
P3 63 8 12
P4 66 8 13
P5 66 9 13
P6 63 16 25
Al igual que en el análisis global, se observan resultados de densidad relativa que bordean el 60
o superior a este valor, por lo tanto se encuentran en un estado de compacidad compactado a
denso y denso, superando generalmente el límite para este tipo de material (DR=60), lo que
implica un bajo riesgo de licuefacción. Es decir, el análisis local compatibiliza con los resultados
obtenidos en el análisis global.
Además se obtuvieron coeficientes de variación que se encuentran entre un 8% y un 20%

aproximadamente, esto se puede relacionar con la gran concentración de datos en valores de
densidad relativa DR entre 50 y 85 que representan casi un 90%, obteniendo un nivel muy bajo
de datos en estado suelto, que provocan este bajo coeficiente de variación.

4. MODELACIÓN DE LOS PARÁMETROS DEL CONTROL DE COMPACTACIÓN.
4. MODELACIÓN DE LOS PARAMETROS DEL CONTROL DE COMPACTACIÓN.
La modelación de los parámetros de la compactación se efectúa mediante la hipótesis de trabajo

que se puede reconstruir la estructura interna del tranque de relave conociendo sus parámetros
y dicha variabilidad producto de su origen mineralógico y proceso constructivo. La gestión de los
parámetros del control de compactación se efectuó mediante los controles de calidad periódicos
realizados del tranque en estudio, debido a esto se conoce la compacidad del tranque,
comportamiento mecánico y riesgo de licuefacción asociadas a una densidad relativa (DR), por
medio de las relaciones estipuladas anteriormente por parámetro de resistencia por punta (qd),
realizando así una batería de datos para la modelación espacial de dichas variables.
Se reconstruirá la estructura interna del tranque de relave para así poder identificar zonas débiles
producto del proceso de compactación que puedan ocasionar zonas de posibles fallas o colapsos
en la estructura interna de éste.
Para la modelación de las variables se utiliza la herramienta computacional Rockwork 15-16, la

cual mediante metodologías de iteración espacial busca representar de manera global el
comportamiento del sólido a modelar. Es por ello que se realiza un análisis de cuál es el método
representa de manera más correcta la estructura interna del tranque de relave.
Como ya bien se ha mencionado, nuestro sector experimental tiene por nombre Tranque
experimental con nombre ficticio M-1, el cual pertenece a la Minera SAP, Es formado por una
pared de 1,4 [km] de longitud que alcanzará, al fin del ciclo de vida del tranque, 237 [m] de altura,
es construido según la técnica talud aguas abajo.
La modelación de la zona de estudio se realizará en tres dimensiones, mediante el uso de un

software basado en un algoritmo anisótropo inverso. Para la ejecución de esto, fueron recogidos
datos de las campañas de ensayos realizados entre los años 2010 y 2013, generando una batería
de datos esencial para la representación de la estructura interna del tranque. Esta modelación
tuvo un supuesto, el cual considera una altura de 31[m] del tranque, producto de la diferencia de
cota entre el año 2010 y el año 2013.
4.1. METODOLOGÍA DE MODELACIÓN ESPACIAL.
La metodología de modelación será una herramienta que no sólo servirá para representar
tranques de relaves, sino que podrá ser implementada en cualquier presa de elevada magnitud.
Para esto será necesario contar con una abundante cantidad de datos, y así establecer una base
gestionada de datos útiles para la futura modelación de éstos. El programa ocupado para el
desarrollo de este modelo es el software Rockworks 15-16, el cual permite la representación y la

interpolación de señales en profundidad. Los diferentes métodos de interpolación han sido
presentados anteriormente y retenemos el método inverso de interpolaciones de la distancia, bajo
condición Anisotrópico.
Se plantea la siguiente metodología para la modelación, basada en los parámetros del control
de compactación.
‒ Definición del área de trabajo, para poner límites a la zona de estudio que se pretende modelar.
Posteriormente determinar la ubicación de las prospecciones realizadas en la campaña de
ensayos y de la profundidad de estas. Cabe señalar que solo es posible generar un sólido
paralelepípedo en el programa lo que es una limitación a la hora de poder representar de
buena forma las características reales del tranque. La distribución de las prospecciones las
define un mapeo que genera el propio programa. A continuación se muestra la ubicación de
los puntos en el área de trabajo.
Figura 4.1. Mapeo de prospecciones en Tranque M-1.


‒ Recopilar los datos de cada una de las prospecciones realizadas, vinculadas con la ubicación
mostrada en la etapa anterior. Se deberá realizar un tratamiento de los datos para ordenar el
ingreso de ellos al software de modelación.
‒ Generar el modelo que más represente a las características de variabilidad estudiadas

anteriormente, pudiendo alisar los datos y que sea claro la representación de dicha
modelación. A continuación se presenta la modelación del sólido correspondiente a el
parámetro de resistencia por punta (qd), en los intervalos de datos 2010-2013, realizando el
supuesto que cada prospección inicia y finaliza en otra campaña, pudiendo interponer los
planos generando alturas aproximadamente de 70 [m] y en algunos periodos se traslapan las
prospecciones por la cercanía en las fechas de la ejecución de los ensayos.
‒ Luego de obtener la modelación completa, se deberá escoger las visualizaciones que mejor
entreguen la información, para que sea útil para detectar las características de la zona. Las
opciones serán realizar perfiles, fence o el sólido en forma completa.
‒ Como última etapa se identificarán las zonas con posibles riesgos de falla, las cuales podrán
verificarse mediantes los análisis realizados ya sea en forma global o local. Otra opción será
que la compañía minera se encargue de verificar mediante ensayos estos sectores de bajas
resistencia o riesgos de fallas arrojados por la modelación, sin embargo, no es alcance de esta
tesis.
4.2. MODELACIÓN Y ANÁLISIS DE LA ZONA DE ESTUDIO.
Para la ejecución del modelo se realizó una grilla de trabajo en la cual se localizaron las ciento
siete prospecciones, utilizando la metodología descrita previamente. Las dimensiones de esta
superficie quedó determinada por un ancho de 460 [m] por un largo de 530 [m] y una altura de
100[m], la cual corresponde a la diferencia de cota entre los años de las campañas de ensayos
consideradas para la modelación. La grilla de la zona de estudio se presenta en la siguiente figura.

Figura 4.2.Representación de la zona de trabajo, con las respectivas prospecciones vistas
desde planta y el Datum de inicio.
4.2.1. Modelación parámetro resistencia por punta (qd).
Se mostrará la representación del tranque en tres dimensiones, en la cual se apreciaran a través

de las distintas tonalidades de colores el nivel de resistencia a la penetración que presenta. Los
tonos anaranjados a rojizos representan bajos valores de qd, y por otro lado los tonos azules
representan altos valores de resistencia del parámetro.
A continuación se presenta las distintas formas de representación del modelamiento espacial del
parámetro de resistencia (qd).

– Modelación sólido.
Figura 4.3. Modelación sólido tres dimensiones (3D).

– Modelación apreciando Vortex(Vértices).
Figura 4.4.Modelación sólido tres dimensiones(3D) apreciando vértices.


– Modelación sólido apreciando malla de puntos.
Figura 4.5. Modelación sólido tres dimensiones (3D) malla de puntos.

Entrando en forma más detallada, se realizan perfiles longitudinales y transversales para cada
uno de los sectores en que se subdivide el tranque. En ellos se puede apreciar las zonas
internas de cada uno de los estratos, pudiendo identificar generalmente resistencias que
varían entre los 3 [MPa] y los 6 [MPa]. Esto concuerda con los datos obtenidos mediante los
análisis estadísticos realizados en el capítulo anterior y con la clasificación cualitativa asociada,
validando la representación del tranque. Existen pequeñas zonas de color rojo, las cuales se
identifican en la parte superior, lo que puede traducirse como espesores sueltos producto de
la maquinaria, erosión, etc., por lo tanto no son considerados datos representativos.

– Perfiles longitudinales.
Figura 4.6. Perfil longitudinal Sector 1 (S-1).




– Perfiles Transversales.
Figura 4.9. Perfil transversal zona superior.

Figura 4.10. Perfil trasversal zona intermedia.

– Visualización del modelo a través de herramienta Fence (Rejillas). Este tipo de modelo
representa la zona en estudio en forma de rejillas, para de esta forma poder observar la
estructura interna del tranque en forma fraccionada. Al igual que todas las representaciones,
esta concuerda con las resistencias en los distintos perfiles, sólo detectando bajas resistencias
en sectores localizados y superficiales.

Figura 4.11. Modelo del parámetro (qd) visualización Fence del tranque M-1.
4.2.2. Modelación estado de compacidad (parámetro cualitativo).
Al igual que para el parámetro resistencia por punta, se realizó la modelación de la clasificación
cualitativa del estado de compacidad, dependiendo si se encuentra suelto o denso. Para esto
existe una simbología con distintas tonalidades para identificar cada uno de estos estados, siendo
en color rojo una capa en estado suelto y en tonos azules cuando se encuentre frente a un
material denso, como se presenta en la Figura 4.12.
Figura 4.12. Simbología estado de compacidad.


A partir de los valores de densidad relativa y de la variación de estos en profundidad para cada
prospección, se estimaron los espesores de capas con su respectiva clasificación del estado de
compacidad. En primera instancia se muestra en forma general la modelación del tranque en tres
dimensiones en la Figura 4.13.
Figura 4.13. Modelación sólido tres dimensiones (3D), compacidad.

Fuente. Elaboración Propia.
Según la simbología empleada, se puede observar que básicamente se presentan capas en

estado compacto a denso y denso, sólo en pequeños sectores se exhiben tonos rosados que
representan zonas en estado compacto, sin embargo, no se aprecian zonas que se encuentren
totalmente sueltas, las que se identificarían en tono rojo. La modelación nuevamente coincide con
los demás análisis realizados a la zona de estudio, ya que en el capítulo 3, cuando se realizó el
análisis estadístico del parámetro densidad relativa, se obtuvieron casi el 90% de los datos en
estado compacto a denso y denso y siempre bajo el 10% de datos en estado suelto y ninguno en
estado muy suelto.
Al igual que el parámetro de resistencia, se realiza en forma detallada para cada perfil la
modelación de esta característica, con la misma simbología utilizada para el modelo en tres
dimensiones. Las figuras 4.14., 4.15. y 4.16., muestran los modelos para cada uno de los
sectores.

Figura 4.14. Perfil longitudinal S-1 compacidad.



Comparando los resultados obtenidos en forma más localizada con respecto al modelo en 3D, se
puede apreciar que nuevamente existe una mayor concentración de zonas en tonos azulados,
los cuales representan zonas compactas a densas y densas, excepto algunos tonos rosados a
rojizos que se pueden expresar básicamente a capas superficiales. Específicamente se marca
una zona roja en el perfil S-3 en la parte inferior de aguas arribas, esta tonalidad se da porque
ese punto es el vértice del paralelepípedo, el cual genera una modelación no efectiva por la
ausencia y distancia en esa parte de datos ingresados.
En resumen se aprecia que el tranque responde a una condición de compacidad de acuerdo a

los estándares exigidos (compacto a denso- denso), esto a través de verificación mediante la
visualización de la modelación de la estructura interna.
4.2.3. Modelación comportamiento mecánico (parámetro cualitativo).
Para complementar la modelación de la zona de estudio, también se incluyó representar el

comportamiento mecánico para cada una de las capas encontradas en cada prospección, el cual
está asociado al parámetro modelado previamente. La simbología empleada se muestra en la
Figura 4.17.
Figura 4.17. Simbología modelación comportamiento mecánico.

El primer modelo se realizó en 3D, el cual se visualiza en la Figura 4.18., que presenta al igual
que el modelo del parámetro estado de compacidad, colores azulados en su mayoría y pequeños
sectores en tonos suavemente rojizos que podrían evidenciar zonas débiles y susceptibles de
sufrir fallas. Sin embargo, como se menciona en los análisis previos, estos sectores pueden ser
despreciables por motivos de su ubicación, ya que se aprecia que se encuentran cercanos a la
superficie, por lo tanto son capas que se encuentran en constante contacto con maquinarias o
sufren de erosión y de esta forma no son representativas del tranque en estudio.

Figura 4.18. Modelación sólido tres dimensiones (3D), comportamiento mecánico.
Completando la modelación, de la misma forma que los demás parámetros se presentan en las
figuras 4.19. a la 4.21., los perfiles longitudinales para cada uno de los sectores en donde se
ubicaron las prospecciones realizadas.
Figura 4.19. Perfil longitudinal S-1 comportamiento mecánico.


Figura4.20. Perfil longitudinal S-2 comportamiento mecánico.
Figura 4.21. Perfil longitudinal S-3 comportamiento mecánico.

Los resultados obtenidos de estos modelos, evidencian una similitud con los del resto de los
parámetros. Se puede apreciar que los perfiles S-1 y S-2 exhiben sectores casi en su totalidad
en tonos azulados, lo que según la simbología empleada representa sectores de un
comportamiento dilatante, es decir, con un bajo o nulo nivel de riesgo de licuefacción. A diferencia
de esto, el perfil S-3 se aprecia zonas con una tonalidad violeta, indicio de características de
zonas de transición del comportamiento mecánico de un estado límite a contractivo, sacando
como conclusión zonas con un potencial de licuación de equilibrio.

CONCLUSIÓN.
I
CONCLUSIÓN.
La metodología de modelación para el control de compactación, fue capaz de reconstruir la

estructura interna de la zona de estudio. A partir de la recopilación y gestión de datos, obtenidos
a través de campañas de ensayos de penetración ligera (PANDA). La modelación a la base de
los resultados de este estudio, será una herramienta global para visualizar zonas de posibles
fallas o colapsos,permitiendo tomar precauciones frente a sectores que puedan originar riesgos
a corto, mediano o largo plazo.
Para la ejecución de esta metodología, fue necesario realizar análisis de variabilidad de los datos
recopilados, tanto en forma global como en forma local, utilizando pruebas estadísticas teóricas
para verificar la hipótesis de homogeneidad del tranque. Estos análisis, permiten acotar los
rangos y la forma en que se distribuyeron los datos, a través de la realización de histogramas
para posteriormente comprobar si los datos obtenidos en la modelación concuerdan con estas
distribuciones.
Los análisis se fundamentan en una base de datos correspondientes a campañas realizadas entre
los años 2010 y 2013, las cuales fueron un total de ciento siete prospecciones, distribuidas en
tres sectores: S-1, S-2 y S-3. Por lo tanto, la modelación presentaba una gran cantidad de batería
de datos para generar un análisis de los parámetros del control de compactación.
Luego de estos análisis, se modelaron tres parámetros: la resistencia a la penetración (qd) en

forma cuantitativa y en forma cualitativa el estado de compacidad junto al comportamiento
mecánico. Los resultados obtenidos para el parámetro resistencia por punta bordearon los valores
entre 3 [MPa] y 6 [MPa], lo que evidencia un buen cumplimiento de la compactación mínima
requerida para que sea ejecutada de forma óptima. En menor porcentaje se observan valores
inferiores a los 3 [MPa], los cuales deberían traducirse en zonas de baja resistencia, en cuanto a
su compactación, y que traen consigo un riesgo de falla frente a solicitaciones externas que
pudiera sufrir el tranque. Sin embargo, la mayoría de estos resultados pueden ser desechados o
despreciables, ya que pertenecen a capas superficiales, que por efecto de factores externos a su
construcción se encuentran en un estado menos denso, tales como el recorrido constante de
maquinarias sobre estos sectores o la erosión a la que se encuentran sujetas estas capas más
externas. Con relación a los parámetros cualitativos representados, se puede mencionar que se
observaron casi en su totalidad tonalidades azuladas en la modelación, que según la simbología
propuesta representan capas en estado compacto a denso. Cabe destacar que sólo se detectan
pequeñas zonas en tonos más cercanos al rojo (violetas), las cuales coinciden con los sectores
de menor resistencia identificados en el modelo de la resistencia por punta, lo que podría deberse
a los mismos factores recién nombrados o a la falta de ensayos para el control de compactación,
que posteriormente puedan ser ingresados al modelo, generando información con una base de
II
datos más amplia y que abarcara toda la zona de estudio.
Es importante señalar que el modelo que muestra el comportamiento mecánico asociado a cada
una de estas capas, trae consigo un potencial de licuefacción, como fue previsto en el capítulo 3.
El modelo de los parámetros cualitativos, en consecuencia beneficia la visualización del estado
en que se encuentra cada uno de los sectores, debido a la mayor utilidad que este modelo entrega
al encargado de control de calidad para detectar zonas con cierto estado de compacidad,
comportamiento mecánico, riesgo de licuefacción u otro parámetro, en comparación a datos
cuantitativos que requieren mayor comprensión y conceptos teóricos.
Mediante la modelación del control de compactación, se pudieron identificar zonas de posibles

fallas, las cuales pueden ser comprobadas a través de los distintos análisis realizados con la base
de datos. Además, la modelación entregará información útil para el control de calidad de la
empresa minera, ya que de esta forma se podrán considerar las siguientes recomendaciones.
‒ Identificar las zonas débiles localizadas en la modelación, mediante estacado en terreno.
‒ Comprobar mediante ensayos de penetración dinámica (SPT) en profundidad las zonas de

fallas detectadas.
‒ Utilizar ensayos de geofísica para rectificar características de variación de la masa de suelo,

especialmente en zonas en donde se identificaron zonas de fallas o cambios de características
ya sea cuantitativas o cualitativas.
‒ Utilizar métodos de control de deformaciones como inclinómetros y control de saturación,

mediante el uso de piezómetros, para disminuir las probabilidades de deformaciones y posibles
riesgos de licuefacción en esas zonas.
Para el uso de la metodología descrita en el último capítulo, se tomaron ciertas consideraciones

que permitan su correcta utilización. Las más importantes fueron.
‒ Tener precisión en la ubicación de las prospecciones en la zona de estudio, a través de la

georeferenciación de los puntos.
‒ Trabajar bajo la hipótesis que se pueden reconstruir la estructura interna de la masa del suelo,
bajo análisis de homogeneidad de los datos analizados
‒ Se deberá tener un correcto tratamiento de los datos, y el usuario deberá contar con los
conocimientos necesarios para el ingreso ordenado de los datos en el software.
BIBLIOGRAFÍA.
BIBLIOGRAFÍA.
– Referencias.
– [1] Blight, G., Troncoso, J., Fourie, A.B. and Wolski, W. (2000, November). Issues in the
geotechnics of mining wastes and tailings. In Proceedings of GeoEng 2000, An International
Conference on Geotechnical and Geological Engineering, 1253-1285, Melbourne.
– [2] Verdugo, R. (2011, January). Seismic performance of slopes and earth and tailings dams
(2010 Maule Earthquake). Fifth International Conference on Geotechnical Earthquake
Engineering (5-ICEGE), 10-13, Santiago de Chile.
– [3] Lee, K. L., & Seed, H. B. (1966). Liquefaction of Satured Sands During Cyclic Loading.
Journal of the Soil Mechanics and Foundations Division, 92 (SM6), 105-134.
– [4] Ishihara, K. (1996).Soil Behavior in Earthquake Geotechnics (1ª Ed).Oxford, Inglaterra:
Oxford Universit Press. Recuperado de: http://www.amazon.com/Behaviour-Earthquake-
Geotechnics-Engineering-Science/dp/0198562241.
– [5] Riemer, et. (1994). Effects of loading frequency and control on the liquefaction behavior
of clean sands. (University of California, Berkeley, California).
– [6] Villavicencio, G., Espinace, R., Palma, J., Fourie, A., & Valenzuela, P. (2014).Failures of
sand tailings dams in a highly seismic country. Canadian Geotechnical Journal, 51, 449–
464.
– [7] Villavicencio, G., Bacconnet, C., Breul, P., Boissier, D., & Espinace, R. (2010, September
16-17). Probabilistic evaluation of the parameters governing the stability of the tailing dams.
1rst International Conference on Information Technology in Geo-Engineering (ICITG),
Shanghai.
– [8] Villavicencio, G. Breul, Espinace, R. & Valenzuela, P. (2012).Control de compactación
con penetrómetro ligero en tranques de relaves, considerando su variabilidad material y
estructural. Revista de la Construcción, 11(1), 119–133.
– [9] Magnan, J (n.d). Les méthodes statisques et probabilistes en mécanique des sols.
Presses des Ponts et Chaussées.
– [10] Chaigneau, L. (2001).Caractérisation des Milieux Granulaires de Surface à l’aide d’un
Pénétromètre (Tesis doctoral, Blaise Pascal-Clermont II Univ, Francia).
– [11] Cambou, B. (1987). Mécanique des milieux granulaires: l’approche microstructurale,
Rhéologie des géomatériaux, éd. F. Darve, presse E.N.P.C, 261–278.
– [12] Biarez, J. & Favre, J.L. (1977). Statistical estimation and extrapolation from
observations. Reports of Organisers IX ICSMFE, (3), 505-509.
– [13] Favre, J. (1980). Milieu continu et milieu discontinu: Mesure statistique indirecte des
paramètres rhéologiques et approche probabiliste de la sécurité’, (Thèse de Doctorat d’Etat,
Université Pierre et Maire Curie, Paris).
– [14] Gourvès, R. & Richard, B. (1995). The Panda Ultra-light Dynamics Penetrometer.
Proceedings of ECSMF, 83–88, Copenhagen, Denmark.
– [15] Gourvès, R. & Zhou, S. (1997). The In Situ Characterization of the Mechanical
Properties of Granular Media with the Help of Penetrometer. 3rd International conférence
on micromécanique of granular média, Duram, 57–60.
– [16] Espinace, R. Villavicencio, G. Palma, J. Peña, A. Bacconnet, C. Gourvés, R. & Bossier,
D. (2007). Nueva Alternativa para el Control de Compactación de Tranques de Relave. El
Penetrómetro PANDA. Revista de la Construcción, 6(2), 3–41.
– [17] Benz, M. (2009).Mesures dynamiques lors du battage du pénétromètre Panda 2 (Tesis
doctoral, Blaise Pascal-Clermont II Univesité, France).
– [18] Villavicencio, G. (2009).Méthodologie pour évaluer la stabilité des barrages de résidus
miniers (Tesis doctoral, Blaise Pascal-Clermont II Université, Francia).
– [19] Boulanger, R; Idriss, I (2003).State Normalization of Penetration Resistance and the
Effect of Overburden Stress on Liquefaction Resistance.
– [20] Espinace, R., Villavicencio, G., & Lemus, L. (2013, August 28-30). The PANDA
technology applied to design and operation of tailings dams. Seminario expuesto en la First
International Seminar on Tailings Management, Santiago, Chile. Recuperado
de:http://www.geotecniaambiental.cl.
– [21] Verdugo, R. (1997, October). Compactación de Relaves IV Chilean Congress of
Geotechnical Engineering. Chilean Society of Geotechnics, Santiago. Universidad Federico
Santa María. Valparaíso. Chile.
– [22] Troncoso, J. (1986).Envejecimiento y estabilidad sísmica de un depósito de residuos
minerales en condición de abandono. ISSN-0716-0348, (22), 147-158.
– [23] Osorio, B. (2009). Resistencia estática y cíclica de relaves integrales (Memoria para
optar al título de ingeniero civil, Universidad de Chile, Santiago, Chile).
– Tesis.
– Asenjo, R. (2010). Efecto del esfuerzo de corte estático inicial en la resistencia cíclica en
una arena de relaves (Memoria para optar al título de ingeniero civil, Universidad de Chile,
Santiago, Chile).
– Durant, B. (2012). Estudio de presas de residuos mineros en el contexto probabilístico (Tesis
para optar al título de ingeniero constructor, Pontificia Universidad Católica de Valparaíso,
Valparaíso, Chile).
– Libros.
– Ishihara, K. (1996).Soil Behavior in Earthquake Geotechnics (1ª Ed).Oxford, Inglaterra:

Oxford Universit Press. Recuperado de: http://www.amazon.com/Behaviour-Earthquake-
Geotechnics-Engineering-Science/dp/0198562241.
– Artículos.
– Seed, H. (1967). Analysis of Soil Liquefaction: Niigata Earthquake. Journal of the Soil
Mechanics and Foundations Division, 93, 83-108. Recuperado de: http://cedb.asce.org/.
– Troncoso, J. (1986).Envejecimiento y estabilidad sísmica de un depósito de residuos
minerales en condición de abandono. ISSN-0716-0348, (22), 147-158.
– Verdugo R. (1995). “Zonificación Sísmica y Geotecnia”, X Congreso Panamericano de
Mecánica de Suelos e Ingeniería de Cimentaciones. Guadalajara, México, 4,647-663.
– Documentos Gubernamentales.
– Ministerio de la Minería. (2007).Decreto supremo nº 248, reglamento para la aprobación de

proyectos, de diseño, construcción, operación y cierre, de los depósitos de relaves.
Santiago, Chile: Servicio Nacional de Geología y Minería del Gobierno de Chile.
– SERNAGEOMIN, (2015).Catastro de Depósitos de Relaves. Depósitos de Relaves Activos
y No Activos 2015. Santiago, Chile: Servicio Nacional de Geología y Minería del Gobierno
de Chile.
– Presentaciones en Power Point.
– Espinace, R. (2008).Normativa y tecnologías innovadoras aplicadas a depósitos de relaves.

[Diapositivas en Power Point]. Recuperado de: www.geotecnicaambiental.cl
– Manuales.
– Rockworks (2014).
– Sol Solutions S.A (n.d)
– Recursos Web.
– Página de información sobre el acontecer nacional de la contaminación de tranques de

relaces en Chile. Catástrofes producto de fallas y colapsos de tranques de relaves en Chile.
Tranque Barahona.
www.relaves.org/caso-tranque-barahona.
[email protected]
Twitter : @relaves_org
– Sociedad Chilena de geotecnia. Como principal propósito, impulsar y divulgar el

conocimiento, desarrollo e investigación de la geotecnia y del amplio espectro de sus
especialidades, donde se incluyen la mecánica de suelos, dinámica de suelos, fundaciones,
mecánica de rocas, geología aplicada a la ingeniería, ingeniería de presas, ingeniería de
túneles, geotecnia ambiental, geotecnia minera, geotecnia vial y geotecnia de puertos, entre
otras.
www.sochige.cl
[email protected]

Tesis Compactación Tranques de Relave

Cargado por

Copyright:

Formatos disponibles

Tesis Compactación Tranques de Relave

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tesis Compactación Tranques de Relave

Cargado por

Copyright:

Formatos disponibles

PONTIFICIA UNIVERSIDAD CATÓLICA DE VALPARAÍSO

METODOLOGÍA PARA LA MODELACIÓN DEL CONTROL DE COMPACTACIÓN EN

Tesis para optar al título de Ingeniero

Gracias totales a todos.

Matías Ojeda García.

Agradezco la oportunidad de haber estudiado en esta institución, y específicamente en esta

El control del proceso de compactación de los tranques de relaves actualmente es puntual,

– Estimación de características cualitativas en función de la densidad relativa (DR),

– Conclusiones que corroboran la hipótesis de trabajo de reconstruir la estructura interna del

1.1. DEPÓSITOS DE RELAVES EN CHILE. 1

1.2. PRINCIPALES MECANISMOS DE FALLAS EN TRANQUES DE RELAVE. 4

1.2.1. Licuefacción o falla de flujo. 4

1.2.2. Inestabilidad del suelo de fundación. 8

1.2.3. Inestabilidad de taludes. 8

1.2.5. Tubificación o piping. 10

1.3. FALLAS HISTORICAS EN TRANQUES DE RELAVES. 13

1.3.1. Tranque Barahona. 13

1.3.2. Tranque El cobre. 14

1.3.3. Tranque Las Palmas y Tranque Adosado. 14

1.4. PROGRESO DE LA NORMATIVA CHILENA EN TRANQUES DE RELAVE. 15

1.4.1. Decreto Supremo N°86. 15

1.4.2. Decreto supremo N°248. 16

1.5. VARIABILIDAD DE LOS PARÁMETROS DEL CONTROL DE COMPACTACIÓN EN

1.5.1. Metodología de análisis de la variabilidad. 20

1.5.2. Descripción de los tranques de relaves experimentales. 21

1.5.3. Análisis de la variabilidad de las características físicas y resistencia in situ. 22

1.5.4. Distribución estadística de los parámetros físicos y mecánicos. 23

1.6. COMPACTACIÓN EN TRANQUES DE RELAVES. 26

1.6.1. Control de compactación en tranques de arenas de relave mediante el

1.7. HIPÓTESIS Y METODOS DE MODELACIÓN. 32

1.7.1. Métodos de interpolación espacial software ROCKWORK. 33

1.7.2. Herramientas de visualización del modelo espacial. 34

2. ESTIMACIÓN DE LA VARIABILIDAD DEL PARÁMETRO DE RESISTENCIA A LA

2.1. INFORMACIÓN GENERAL DEL TRANQUE DE RELAVE EN ESTUDIO (TRANQUE M-

2.2. RECOPILACIÓN DE INFORMACIÓN. 38

2.2.1. Descripción del área de estudio. 39

2.2.2. Emplazamiento de prospecciones. 39

2.3. ANÁLISIS ESTADÍSTICO DE LOS DATOS RECOPILADOS. 40

2.3.1. Resumen estadístico. 40

2.3.2. Distribución estadística del parámetro de resistencia a la penetración (qd). 41

2.3.3. Leyes de probabilidad que se ajustan a la distribución del parámetro de

2.3.4. Conclusiones de variabilidad Tranque M-1. 43

3. ESTIMACIÓN DE LA COMPACIDAD, COMPORTAMIENTO MECÁNICO Y RIESGO DE

3.1. VARIACIÓN DE LA RESISTENCIA POR PUNTA (qd). 47

3.2. RELACIÓN ENTRE LA RESISTENCIA POR PUNTA NORMALIZADA (qdN1) Y EL

3.2.1. Normalización del parámetro resistencia (qd). 56

3.3. CLASIFICACIÓN Y ANÁLISIS DEL GRADO DE COMPACTACION, COMPORTAMIENTO

4. MODELACIÓN DE LOS PARAMETROS DEL CONTROL DE COMPACTACIÓN. 73

4.1. METODOLOGÍA DE MODELACIÓN ESPACIAL. 73

4.2. MODELACIÓN Y ANÁLISIS DE LA ZONA DE ESTUDIO. 75

4.2.1. Modelación parámetro resistencia por punta (qd). 76

4.2.2. Modelación estado de compacidad (parámetro cualitativo). 81

4.2.3. Modelación comportamiento mecánico (parámetro cualitativo). 84

La problemática antes mencionada ha estado presente a lo largo de la historia minera en Chile,

Específicamente, sobre los depósitos de relaves en Chile, existen cinco maneras de

Actualmente el control de compactación de los tranques de relave durante la fase de operación,

Para el desarrollo de esta metodología se ejecutan tres objetivos específicos.

‒ Analizar la variabilidad de los parámetros que definen el control de compactación en el tranque

de distribución teóricas, con la finalidad de demostrar la hipótesis de homogeneidad de los