Caracterizacion Geoquimica y Mineralogica de Alteraciones

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 221

UNIVERSIDAD DE CHILE

FACULTAD DE CIENCIAS FSICAS Y MATEMTICAS

DEPARTAMENTO DE GEOLOGA

CARACTERIZACIN GEOQUMICA Y MINERALGICA DE ALTERACIONES


HIDROTERMALES EN PRFIDO CUPRFERO ESCONDIDA

TESIS PARA OPTAR A MAGSTER EN CIENCIAS, MENCIN GEOLOGA

MEMORA PARA OPTAR A TTULO DE GELOGO

EFRAN EDUARDO JUVENAL CRDENAS DEL RO

PROFESOR GUA:

BRIAN TOWNLEY CALLEJAS

MIEMBROS DE LA COMISIN:

JULIN ORTIZ CABRERA

KATJA DECKART

DUBER JOFR PINTO

SANTIAGO DE CHILE

2015
RESUMEN

En la minera del cobre existen minerales de ganga que afectan la recuperacin


de la mena, disminuyendo la eficiencia del proceso extractivo, por lo que predecir la
presencia de estos minerales es fundamental para formular una planificacin minera
que optimice la recuperacin de la mena.

Con el objetivo de crear una herramienta que permita predecir la presencia de


alteraciones hidrotermales a partir de datos geoqumicos, se caracteriza cada una de
las zonas de alteracin cruzando informacin de mapeo geolgico, geoqumica y
mineraloga Qemscan. Debido a que el mtodo de disolucin por agua regia fue
utilizado en la obtencin de la base de datos, se definen las zonas de alteracin
hidrotermal en funcin de la mineraloga que es capaz de disolver este cido.

La metodologa para la clasificacin consta de la definicin de centroides para


cada alteracin hidrotermal en un espacio definido por las concentraciones de
elementos qumicos, para luego calcular un grado de pertenencia a cada uno de las
alteraciones predefinidas. De esta manera, el grupo de elementos que mejor discrimina
entre alteraciones hidrotermales es Al-Mg-K-Ca-Na, logrando definir hasta 3 zonas de
alteracin en donde predominan algunos minerales como biotita, clorita y pirita.

La clasificacin de los datos Qemscan indica que el modelo replica las razones
de minerales que caracterizan las alteraciones, a pesar de que el mapeo no sea
consistente con la mineraloga. As, el modelo de clasificacin tiene como aplicacin
directa ayudar al mapeo de alteraciones sugiriendo las clases e indicando los
principales controles minerales. Adems, es posible utilizar la correlacin entre los
porcentajes de clasificacin y los grados de pertenencia para indicar si una muestra se
encuentra mapeada correctamente segn sus patrones geoqumicos.

El algoritmo utilizado puede ser implementado en trabajos futuros que incorporen


el anlisis de otras variables de inters con el objetivo de integrar la mayor cantidad de
informacin valiosa para la planificacin.

i
A mis padres, Efran y Vernica
A mi hermana, Valentina
A mi sobrino, Felipe
A Juve

ii
Agradecimientos

Quiero agradecer a todas las personas que me apoyaron e hicieron de este trabajo una
experiencia nica.

A mis profesores, Brian Townley y Julin Ortiz, por confiar en m, por el apoyo y por
guiarme en lo largo de este proceso. A los miembros de la comisin, Katja Deckart y
Duber Jofre, por disposicin a entregar su experiencia y conocimiento.

A CONICYT y AMTC, por las becas y viajes. A la gente de BHP y Escondida por su
ayuda.

Al laboratorio ALGES, por la buena onda, la ayuda y la gua. A los compaeros de


oficina Alges VIP: Pa, Roberto, Felipe, Yerko, Rodrigo y Sebastin.

A mis hermanos de Diogenes Magallanes, Gabriel, Daniel, Constanza, Rodrigo, Cristian


y Emilio, por la fuerza, la paciencia y la inspiracin.

A mis amigos por la buena onda, la vibra, las emociones, el apoyo.

A mi familia por estar presente, aguantarme, aconsejarme, y sobre todo, por creer en
m.

iii
Tabla de contenido

1 Introduccin ...............................................................................................................1

1.1 Formulacin del problema ...................................................................................1

1.2 Hiptesis de trabajo.............................................................................................2

1.3 Objetivos .............................................................................................................3

1.4 Metodologa de estudio .......................................................................................3

1.4.1 Caracterizacin geoqumica de zonas de alteraciones hidrotermales ..........3

1.4.2 Caracterizacin mineralgica de las alteraciones hidrotermales ..................4

1.4.3 Correlacin de geoqumica y mineraloga ....................................................4

1.4.4 Clasificacin de muestras a partir de datos geoqumicos .............................4

1.5 Caractersticas relevantes del yacimiento Escondida .........................................4

2 Marco geolgico del yacimiento Escondida ...............................................................7

2.1 Rocas estratificadas ............................................................................................8

2.1.1 Formacin La Tabla (Carbonfero-Prmico) .................................................8

2.1.2 Formacin Agua Dulce (Trisico Superior) ...................................................9

2.1.3 Formacin El Profeta y Santa Ana (Jursico-Cretsico) ...............................9

2.1.4 Formacin Augusta Victoria (Paleoceno-Eoceno inferior) .......................... 10

2.2 Rocas intrusivas ................................................................................................ 10

2.2.1 Complejo Intrusivo Feldesptico Escondida (CIFE).................................... 10

2.2.2 Prfido Rioltico........................................................................................... 11

2.2.3 Stocks diorticos y diques ........................................................................... 11

2.3 Brechas ............................................................................................................. 12

2.4 Gravas ............................................................................................................... 12

2.5 Estructuras ........................................................................................................ 12

2.6 Alteracin hidrotermal y mineralizacin ............................................................. 15

2.6.1 Alteracin potsica ..................................................................................... 16

iv
2.6.2 Alteracin propiltica ................................................................................... 17

2.6.3 Silicificacin ................................................................................................ 17

2.6.4 Vetillas Tempranas ..................................................................................... 18

2.6.5 Clorita-sericitacuarzo ................................................................................ 18

2.6.6 Cuarzo-sericita............................................................................................ 19

2.6.7 Arglica avanzada ....................................................................................... 19

3 Marco Terico .......................................................................................................... 20

3.1 Alteracin hidrotermal ....................................................................................... 20

3.2 Anlisis ICP-AES............................................................................................... 24

3.2.1 Disolucin por Agua Regia ......................................................................... 24

4 Anlisis de Datos ..................................................................................................... 30

4.1 Anlisis geoqumicos y mapeo de sondajes ...................................................... 30

4.2 Definicin del sistema ....................................................................................... 31

4.2.1 Porcentaje de disolucin ............................................................................. 39

4.2.2 Abundancia de minerales ........................................................................... 43

4.3 Anlisis univariable............................................................................................ 48

4.4 Depuracin de base de datos geoqumicos ...................................................... 51

4.5 Anlisis de componentes principales ................................................................ 53

4.5.1 Elementos mayores .................................................................................... 53

4.5.2 Elementos traza .......................................................................................... 57

4.6 Clasificacin Previa ........................................................................................... 60

5 Mineraloga avanzada.............................................................................................. 66

5.1 Base de datos QEMSCAN ................................................................................ 66

5.1.1 Depuracin de base de datos Qemscan ..................................................... 66

5.1.2 Mineraloga Soluble .................................................................................... 68

5.1.3 Mineraloga y geoqumica ........................................................................... 72

v
6 Anlisis preliminar de variables y definicin de grficos de inters ......................... 77

6.1 Seleccin de elementos .................................................................................... 78

6.2 Alteraciones y mineralizacin por sondaje ........................................................ 81

6.3 Grficos de control mineral ............................................................................... 87

7 Modelo de clasificacin de alteraciones hidrotermales ............................................ 92

7.1 Definicin de muestras por alteracin hidrotermal ............................................ 93

7.2 Seleccin de elementos segn desempeo ...................................................... 96

7.3 Anlisis de los clusters por alteracin hidrotermal .......................................... 106

7.4 Nmero ptimo de clusters ............................................................................. 110

7.5 Grados de pertenencia de muestras con una alteracin por tramo ................. 115

7.6 Patrones mineralgicos ................................................................................... 119

7.7 Reclasificacin de alteracin flica indiferenciada ........................................... 124

7.8 Grados de pertenencia de muestras con dos alteraciones por tramo ............. 128

8 Discusin ............................................................................................................... 131

8.1 Limitantes ........................................................................................................ 139

9 Conclusiones y recomendaciones ......................................................................... 141

10 Bibliografa ............................................................................................................. 144

Anexo A ........................................................................................................................ 149

Cantidad de datos por litologa, alteracin y minzone ............................................... 149

Histogramas por alteracin para elementos mayores y traza ................................... 152

Boxplots por alteracin (log base 10) ........................................................................ 155

Histogramas litologas 31 y 50 .................................................................................. 157

Histogramas de minerales segn alteracin ............................................................. 160

Mineraloga vs elementos mayores........................................................................... 161

Mineraloga vs abundancia relativa (tamao indica mayor ndice de abundancia) ... 168

Diagramas de dispersin de elementos mayores ..................................................... 170

vi
Scatter de elementos con menor traslape entre alteraciones ................................... 173

Matrices de confusin entre clasificacin y mapeo de alteraciones considerando


distintos grupos de elementos mayores .................................................................... 183

Matriz de correlacin para elementos mayores [wt%] ............................................... 186

Matriz de correlacin para elementos traza [ppm] .................................................... 187

Matriz de correlacin para elementos traza [ppm] (continuacin) ............................. 188

Componentes principales elementos mayores ......................................................... 189

Valores propios ...................................................................................................... 189

Vectores propios .................................................................................................... 189

Componentes principales elementos traza ............................................................... 190

Valores propios ...................................................................................................... 190

Vectores propios .................................................................................................... 191

Anexo B ........................................................................................................................ 194

Boxplots muestras con una alteracin ...................................................................... 194

Al-Mg-K-Ca-Na ...................................................................................................... 194

Al-Mg-K-Ca ............................................................................................................ 194

Al-Mg-K-Na ............................................................................................................ 195

Boxplots muestras con dos alteraciones ................................................................... 195

Al-Mg-K-Ca-Na ...................................................................................................... 195

Al-Mg-K-Ca ............................................................................................................ 196

Al-Mg-K-Na ............................................................................................................ 196

Plantas de muestras con una alteracin ................................................................... 197

Al-Mg-K-Ca-Na ...................................................................................................... 197

Al-Mg-K-Ca ............................................................................................................ 198

Al-Mg-K-Na ............................................................................................................ 199

Porcentaje de clasificacin versus rangos de porcentaje de alteracin por tramo .... 200

vii
Al-Mg-K-Ca-Na ...................................................................................................... 200

Al-Mg-K-Ca ............................................................................................................ 200

Al-Mg-K-Na ............................................................................................................ 200

Porcentaje de clasificacin versus rangos de grados de pertenencia ....................... 201

Al-Mg-K-Ca-Na ...................................................................................................... 201

Al-Mg-K-Ca ............................................................................................................ 201

Al-Mg-K-Na ............................................................................................................ 201

Anexo C........................................................................................................................ 202

Tabla de Estadgrafos para elementos en ppm......................................................... 202

Tabla de Estadgrafos para mineraloga QemScan .................................................. 204

ndice de Figuras

Figura 1.1. Franja metalognica del Eoceno Superior Oligoceno Inferior. Se


esquematiza el Sistema de Fallas de Domeyko (Maksaev, 2001). ..................................5

Figura 2.1. Mapa geolgico del yacimiento Escondida (modificado de Richards et al.,
2001). ...............................................................................................................................8

Figura 2.2. Esquemas de las estructuras principales de Escondida. A) Geologa y


estructuras de Escondida. B) Geologa y estructuras del distrito Escondida en general
(Padilla et al., 2001). ...................................................................................................... 14

Figura 2.3. Mapa de alteraciones hidrotermales del yacimiento Escondida (Modificado


de Padilla et al., 2001). ................................................................................................... 15

Figura 3.1. Asociaciones de minerales de alteracin comn en sistemas hidrotermales


(Corbett & Leach, 1997). ................................................................................................ 21

Figura 3.2. Modelo de zonacin de alteracin y mineralizacin para prfidos cuprferos


(Sillitoe, 2010). ............................................................................................................... 22

viii
Figura 3.3. Esquema generalizado de la secuencia de mineralizacin y alteracin en un
prfido cuprfero en funcin de la paleoprofundidad y tiempo de formacin (Sillitoe,
2010). ............................................................................................................................. 23

Figura 3.4. Porcentaje de minerales disueltos segn el procedimiento consecutivo de


cido oxlico y agua regia. El grosor de la barra representa el nmero de muestras
utilizada para el estudio (Church et al., 1987). ............................................................... 26

Figura 3.5. Porcentaje de disolucin de varios minerales, obtenidos al comparar


digestin por agua regia relativa la digestin por HF-HNO3-HClO4 (Church et al., 1987).
....................................................................................................................................... 27

Figura 3.6. Porcentaje de recuperacin de elementos por mineral (datos de Snll &
Liljefors, 2000). ............................................................................................................... 29

Figura 4.1. Cantidad de muestras por litologa segn primera y segunda litologa
mapeada por tramo (ver Anexo A). ................................................................................ 32

Figura 4.2. Cantidad de muestras por alteracin para las litologas PF (31), PC (33) y
AN (50) (Ver Anexo A).................................................................................................... 33

Figura 4.3. Cantidad de muestras por zonas de mineralizacin o minzone (ver Tabla 4.1
para codificacin). .......................................................................................................... 34

Figura 4.4. Distribucin espacial de los datos geoqumicos segn la litologa. Arriba:
Perfil EW N = 107350 50 m. Abajo: Planta cota = 2750 50 m. ................................. 37

Figura 4.5. Distribucin espacial de los datos geoqumicos segn alteracin hidrotermal.
Arriba: Perfil EW N = 107350 50 m. Abajo: Planta cota = 2750 50 m....................... 38

Figura 4.6. Distribucin espacial de los datos geoqumicos segn mineralizacin. Perfil
EW N = 107350 50 m. ................................................................................................. 39

Figura 4.7. Boxplot de disolucin total calculada para elementos mayores. .................. 41

Figura 4.8. Porcentaje relativo promedio de elementos disuelto por alteracin. ............ 42

Figura 4.9. Porcentaje normalizado promedio de elementos traza disueltos. ................ 43

Figura 4.10. Esquema de grado de destruccin de los feldespatos. .............................. 44

Figura 4.11. Abundancia relativa de minerales en cada alteracin hidrotermal. ............ 45

ix
Figura 4.12. Abundancia relativa de todos los minerales por alteracin hidrotermal. ..... 47

Figura 4.13. Histogramas de elementos mayores por alteracin hidrotermal................. 49

Figura 4.14. Histogramas de algunos elementos separados por alteracin hidrotermal.


....................................................................................................................................... 50

Figura 4.15. Boxplots con eje logartmico de elementos por alteracin hidrotermal. ...... 51

Figura 4.16. Varianza acumulada segn ACP de elementos mayores. ......................... 54

Figura 4.17. Diagrama de sedimentacin para elementos mayores. ............................. 54

Figura 4.18. Crculos de correlacin para elementos mayores. ..................................... 55

Figura 4.19. Diagrama de dispersin entre PC1 y PC2. Los colores indican alteraciones
hidrotermales.................................................................................................................. 56

Figura 4.20. Varianza acumulada segn ACP de elementos traza. ............................... 58

Figura 4.21. Diagrama de sedimentacin para elementos traza. ................................... 58

Figura 4.22. Diagramas de componentes principales para elementos traza. ................. 59

Figura 4.23. Diagrama de dispersin entre PC1 y PC2. Los colores indican alteraciones
hidrotermales.................................................................................................................. 60

Figura 4.24. Diagrama de clasificacin Al-Fe-Mg (Townley et al., 2012). ...................... 61

Figura 4.25. Diagrama de clasificacin Ca-Fe-S (Townley et al., 2012). ........................ 62

Figura 4.26. Diagrama de discriminacin entre alteraciones usando ndice


(Al+K)/(Na+Ca+Mg) vs Sc*V. Arriba se indica con una lnea roja el lmite para la
clasificacin como alteracin potsica. Abajo se indica la densidad de datos (Townley et
al., 2012). ....................................................................................................................... 63

Figura 4.27. Diagrama Mg-Fe-S mostrando clasificacin de base de datos (Townley et


al., 2012). ....................................................................................................................... 65

Figura 5.1. Concentraciones promedio de minerales Qemscan por alteracin


hidrotermal. .................................................................................................................... 68

Figura 5.2. Proporciones relativas de minerales solubles por alteracin hidrotermal. .... 69

x
Figura 5.3. Histogramas de distribucin de mineraloga cuantitativa Qemscan por
alteracin hidrotermal. .................................................................................................... 70

Figura 5.4. Diagramas de clasificacin segn Townley et al., 2012. El tamao de los
crculos indica abundancia relativa de pirita a la izquierda y clorita a la derecha. .......... 71

Figura 5.5. Diagrama de clasificacin segn Townley et al. (2012) para alteracin K1. El
tamao de los crculos indica abundancia relativa de biotita. ......................................... 72

Figura 5.6. Histogramas de razones de minerales Qemscan. ........................................ 74

Figura 5.7. Grficos molares de dispersin de Al-Mg, Al-K, Mg-K. El tamao de los
crculos indica la razn clorita/sericita. .......................................................................... 75

Figura 5.8. Grficos molares de dispersin de Al-Mg, Al-K, Fe-Mg, Fe-S, Mg-K y Mg-S.
El tamao de los crculos indica la razn clorita/pirita. ................................................... 76

Figura 6.1. Casos posibles de traslapes entre las alteraciones 30 y 31. ........................ 79

Figura 6.2. Boxplot de minzones para Cu total. .............................................................. 81

Figura 6.3. Concentraciones de elementos a lo largo de sondajes resaltando el contacto


entre alteraciones. .......................................................................................................... 82

Figura 6.4. Concentraciones de elementos a lo largo de sondajes resaltando


variaciones no relacionadas al cambio de alteracin. .................................................... 84

Figura 6.5. Concentraciones de elementos a lo largo de sondajes resaltando


variaciones de potasio en una misma alteracin. ........................................................... 85

Figura 6.6. Razn de minerales Qemscan clorita:sericita vs K wt%. Tamao indica


concentracin de biotita Qemscan. ................................................................................ 86

Figura 6.7. Diagrama ternario Cu-Fe-S indicando zonas de pirita-sulfuros de Cu y


enriquecimiento de Fe y S. ............................................................................................. 88

Figura 6.8. Diagrama ternario Mg-Fe-S indicando zonas de control mineral de biotita-
clorita, pirita y transicin entre ambas. ........................................................................... 89

Figura 6.9. Diagrama ternario Mg-K-Ca indicando zonas de Control biotita-clorita y


Control de sericita. ......................................................................................................... 91

Figura 7.1. Abundancia de minerales Qemscan segn alteracin hidrotermal............... 95

xi
Figura 7.2. Diferentes particiones con igual distancia entre centroides pero con distinta
orientacin para igual set de datos (Tomado de Rezaee, 2010). ................................... 98

Figura 7.3. Boxplots grados de pertenencia por alteracin para magnesio. ................. 102

Figura 7.4. Boxplots grados de pertenencia por alteracin para azufre. ...................... 103

Figura 7.5. Promedio de grados de pertenencia por alteracin mapeada. ................... 107

Figura 7.6. Planta Cota 2750 50 m. Tamao de las muestras representa el grado de
pertenencia a un cluster. .............................................................................................. 108

Figura 7.7. Planta Cota 2750 50 m. Tamao de las muestras representa el grado de
pertenencia a un cluster (continuacin). ....................................................................... 109

Figura 7.8. Metodologa de agrupacin de k-means (tomado de Han et al., 2012). ..... 112

Figura 7.9. Media de silhouettes vs nmero de clusters para grupos de elementos


mayores y trazas. ......................................................................................................... 113

Figura 7.10. Boxplots de grados de pertenencia a clusters K, SCC y S calculados para


Al-Mg-K-Ca-Na para muestras con una alteracin por tramo....................................... 116

Figura 7.11. Grados de pertenencia calculados por rango para muestras con solo una
alteracin por tramo para diferentes grupos de elementos. ......................................... 118

Figura 7.12. Diagramas ternarios para las muestras Qemscan. El tamao de los crculos
indica mayor valor de la razn (Biotita+Clorita)/Pirita. .................................................. 121

Figura 7.13. Composicin media de minerales para muestras clasificadas


correctamente y sus respectivos errores. ..................................................................... 122

Figura 7.14. Boxplots composicin media de minerales para muestras clasificadas


correctamente y sus respectivos errores. ..................................................................... 123

Figura 7.15. Boxplots de grados de pertenencia de alteracin flica indiferenciada. .... 125

Figura 7.16. Diagramas ternarios para alteracin flica indiferenciada. ........................ 126

Figura 7.17. Planta Cota= 2750 50 m. Grados de pertenencia a los clusters K, SCC y
S para Al-Mg-K-Ca-Na. ................................................................................................ 127

Figura 7.18. Distribucin espacial de clasificacin de alteracin flica indiferenciada. . 128

xii
Figura 7.19. Boxplots de grados de pertenencia a clusters K, SCC y S calculados para
Al-Mg-K-Ca-Na para muestras con dos alteraciones por tramo. .................................. 129

Figura 7.20. Ejemplo del par SCC S. ........................................................................ 129

Figura 7.21. Porcentaje de clasificacin en funcin de las proporciones de alteracin por


tramo y grado de pertenencia para los pares K SCC y SCC S, considerando el
grupo de elementos Al-Mg-K-Na. ................................................................................. 130

Figura 8.1. Boxplots de concentracin de arcillas por alteracin. ................................. 132

Figura 8.2. Sondaje con mapeos de alteracin incongruente con geoqumica. ........... 133

Figura 8.3. Distribucin espacial de alteraciones hidrotermales. Tamao indica grado de


incertidumbre. Perfil N = 107350 50 m. Planta cota = 2750 50 m. ......................... 136

Figura 8.4. Distribucin espacial de litologas. Tamao indica grado de incertidumbre.


Perfil N = 107350 50 m. Planta cota = 2750 50 m. ................................................. 136

Figura 8.5. Distribucin espacial de zonas de mineralizacin. Tamao indica grado de


incertidumbre. Perfil N = 107350 50 m. Planta cota = 2750 50 m. ......................... 137

Figura 8.6. Distribucin espacial de clasificacin de alteraciones hidrotermales. Tamao


indica grado de incertidumbre. Perfil N = 107350 50 m. Planta cota = 2750 50 m. 138

Figura 8.7. Razn K/Al vs Mg/Al. .................................................................................. 141

ndice de Tablas

Tabla 1.1. Produccin anual de cobre en miles de toneladas en Escondida (BHP Billiton,
2012). ...............................................................................................................................6

Tabla 1.2. Produccin anual de oro y plata en miles de onzas en Escondida (BHP
Billiton, 2012)....................................................................................................................6

Tabla 1.3. Recursos minerales de cobre al 30 de junio de 2012 en Escondida (BHP


Billiton, 2012)....................................................................................................................6

Tabla 1.4. Reservas de mena de cobre al 30 de junio de 2012 en Escondida (BHP


Billiton, 2012)....................................................................................................................6

xiii
Tabla 2.1. Edades de Formaciones de rocas estratificadas. ............................................8

Tabla 2.2. Etapas de alteracin hidrotermal y vetillas asociadas. .................................. 16

Tabla 4.1. Codificacin de zonas mineralgicas. ........................................................... 33

Tabla 4.2. Cantidad de muestras con nica alteracin por tramo. ................................. 35

Tabla 4.3. Cantidad de datos por tramo con mezcla de alteraciones. ............................ 35

Tabla 4.4. Estadgrafos de porcentaje de disolucin calculada. ..................................... 40

Tabla 4.5. Elementos con ndice de correlacin superior o igual a 0.93. ....................... 52

Tabla 4.6. Clasificacin entre clases propuestas por Townley et al. (2012). .................. 64

Tabla 5.1. Minerales por grupo mineral identificados en la base de datos Qemscan ..... 67

Tabla 5.2. Cantidad de muestras de mineraloga QEMSCAN por litologa y alteracin. 67

Tabla 5.3. Frmula estructural de los minerales analizados por QEMSCAN


considerados como solubles en agua regia. .................................................................. 73

Tabla 6.1. Percentiles definidos para el 75% de los datos por alteracin para magnesio.
....................................................................................................................................... 78

Tabla 6.2. Menor porcentaje de traslape obtenido para las combinaciones de


alteraciones indicadas. ................................................................................................... 79

Tabla 6.3. Pares de alteraciones con porcentaje de traslape menor a un 15%. ............. 80

Tabla 6.4. Porcentaje de muestras por zona en sistema Cu-Fe-S. ................................ 88

Tabla 6.5. Porcentaje de muestras por zona en sistema Mg-Fe-S. ................................ 90

Tabla 6.6. Porcentaje de muestras por zona en sistema Mg-K-Ca. ............................... 91

Tabla 7.1. ndices de abundancia relativa de minerales para la definicin de clusters. . 93

Tabla 7.2. ndices de mapeo de abundancia relativa de minerales con mineraloga


Qemscan. ....................................................................................................................... 96

Tabla 7.3. Cantidad muestras por base de datos con una alteracin por tramo. ........... 96

Tabla 7.4. Matriz de error asociada al Mg. Media matriz= 0.172. Media diagonal = 0.60.
..................................................................................................................................... 101

xiv
Tabla 7.5. Errores de clasificacin para elementos mayores. ...................................... 101

Tabla 7.6. Errores de clasificacin para elementos traza. ............................................ 104

Tabla 7.7. Error promedio y promedio de porcentajes mximos de clasificacin para


grupos de elementos mayores. .................................................................................... 105

Tabla 7.8. Error promedio y promedio de porcentajes mximos de clasificacin para


grupos de elementos traza. .......................................................................................... 106

Tabla 7.9. Matrices de confusin de Al-Mg-K-Ca-Na y Al-Mg-K-Na. ............................ 110

Tabla 7.10. Distancia euclidiana entre centros de alteraciones (Al-Mg-K-Ca-Na). ....... 111

Tabla 7.11. Porcentaje de clasificacin de alteraciones en los clusters calculados por K-


means para Al-Mg-K-Ca-Na y Al-Mg-K-Ca................................................................... 113

Tabla 7.12. Distancia entre centros predefinidos y calculados por K-means para Al-Mg-
K-Ca-Na y Al-Mg-K-Ca. ................................................................................................ 114

Tabla 7.13. Porcentaje de clasificacin de alteraciones mapeadas en clusters


redefinidos para distintos grupos de elementos. .......................................................... 115

Tabla 7.14. Matrices de confusin para los datos de mineraloga Qemscan. .............. 120

Tabla 7.15. Resultados de clasificacin de alteracin flica indiferenciada. ................. 124

Tabla 8.1. Porcentaje de clasificacin de alteraciones mapeadas en clusters redefinidos


para distintos grupos de elementos. ............................................................................. 134

xv
1 Introduccin

Este proyecto de tesis de Magster en Ciencias Mencin Geologa, se enmarca


dentro del proyecto MQAlt del laboratorio ALGES (Advanced Laboratory for
Geostatistical Supercomputing) de AMTC (Advanced Mining Technology Center), el cual
tiene como objetivo desarrollar herramientas para el anlisis y modelamiento de las
alteraciones hidrotermales presentes en la mina.

1.1 Formulacin del problema

En la gran minera, uno de los mayores desafos es caracterizar a tiempo el tipo


de roca que ser procesado. No tan solo los minerales de mena son relevantes, sino
que tambin los minerales asociados que, a pesar de no poseer valor econmico,
influyen directamente en el procesamiento de la roca. La presencia de algunos
minerales de ganga en distintas proporciones dentro de la roca pueden controlar su
comportamiento en planta. Si no se conoce ni cuantifica la presencia de estos
minerales, pueden producirse importantes prdidas econmicas durante el
procesamiento metalrgico del material.

Los procesos hidrotermales que producen la mineralizacin econmica en un


sistema tipo prfido cuprfero como Escondida, tambin son responsables de alterar la
roca que la aloja. Esta alteracin corresponde a asociaciones de minerales que dan
informacin acerca de los procesos que han afectado a la roca. A estas asociaciones se
les conoce como alteracin hidrotermal y estn relacionadas directamente con los
procesos mineralizadores.

Trabajos anteriores han demostrado que es posible realizar una clasificacin de


las alteraciones hidrotermales de Escondida en funcin de la geoqumica y los datos de
mapeo geolgico de sondajes (Townley, et al., 2012; Luca, et al., 2012). En el presente
trabajo se propone refinar esta clasificacin al poseer datos de la intensidad de
1
alteracin para cada muestra (leve, moderada e intensa), y datos de mineraloga
avanzada obtenidos por anlisis Qemscan.

El principal problema con el que se debe lidiar para realizar esta clasificacin es
que la base de datos geoqumica no es la ideal para el estudio de la mineraloga de
alteracin, ya que el mtodo de disolucin por agua regia no es total. Este cido no es
capaz de disolver la mayora de los silicatos, sino que solo sulfuros, algunos sulfatos,
xidos de Fe-Mn, carbonatos, entre otros minerales (Chao, 1984; Chao & Sanzalone,
1992; Church, et al., 1987; Hall, 1999; Snll & Liljefors, 2000). Otro problema a
considerar, es que la base de datos se ha generado a lo largo de los aos y por lo tanto
los procedimientos de anlisis geoqumico han podido variar dependiendo de la
cantidad de laboratorios que los han realizado.

Este trabajo busca correlacionar las bases de datos geoqumica, mineralgica y


mapeo de sondajes para crear un modelo de clasificacin de las alteraciones
hidrotermales. Este modelo podr ser usado en estudios posteriores que apunten a
predecir el comportamiento de la roca en la planta a partir de la mineraloga utilizando
criterios geo-minero-metalrgicos.

1.2 Hiptesis de trabajo

La presencia y abundancia relativa de ciertos minerales definen un tipo de


alteracin hidrotermal, la cual imprime composiciones qumicas particulares a la roca.
En el caso de Escondida, es posible definir la presencia de alteraciones hidrotermales a
partir de la base de datos geoqumica, previamente asociada y calibrada a los datos de
logueo de sondajes y mineraloga avanzada Qemscan.

Se postula que correlacionar la informacin de mapeo de sondajes, mineraloga


avanzada Qemscan y base de datos geoqumico, permitir encontrar criterios para
cuantificar la presencia de alteraciones hidrotermales.

2
1.3 Objetivos

El objetivo general de este estudio es definir criterios y parmetros geoqumicos


de caracterizacin mineral cuantitativa para su posterior uso en modelamiento espacial.
De esta manera, se plantea construir un modelo cuantitativo de alteraciones
hidrotermales basado en datos geoqumicos, mineralgicos y mapeo de sondajes para
ser utilizado en la prediccin del comportamiento del material durante el procesamiento
metalrgico.

Los objetivos especficos son los siguientes:

1) Determinar la mineraloga y geoqumica de las zonas de alteracin hidrotermal.


2) Analizar correlacin entre mineraloga y geoqumica de las zonas de alteracin
hidrotermal.
3) Definir zonas de alteracin hidrotermal en funcin de la geoqumica.

1.4 Metodologa de estudio

1.4.1 Caracterizacin geoqumica de zonas de alteraciones hidrotermales

La metodologa para definir la mineraloga de las zonas de alteracin consiste en


una revisin bibliogrfica acerca de la solubilidad de los minerales en agua regia.

El anlisis estadstico de los datos geoqumicos permitir determinar los


elementos que presenten una distribucin preferencial por alteracin hidrotermal. Se
hace la distincin entre elementos mayores y traza, ya que los primeros se
correlacionan con la presencia de minerales de mena y ganga.

3
1.4.2 Caracterizacin mineralgica de las alteraciones hidrotermales

El anlisis estadstico de los datos de mineraloga Qemscan permitir determinar


las asociaciones minerales representativas de cada zona de alteracin considerando
solo los grupos minerales solubles. Adems, se analizarn las principales diferencias
entre alteraciones y se contrasta con las diferencias del protocolo de mapeo.

1.4.3 Correlacin de geoqumica y mineraloga

Considerando la mineraloga y geoqumica definidas para cada zona de


alteracin, se estudia el comportamiento conjunto de los elementos mayores y su
correlacin con las principales fases presentes. De esta manera, se obtendr una
aproximacin de las principales fuentes de los elementos qumicos.

1.4.4 Clasificacin de muestras a partir de datos geoqumicos

Considerando la informacin bibliogrfica, la seleccin de elementos qumicos


asociados a la presencia mineral, se crea un modelo de clasificacin de las alteraciones
a partir de los datos geoqumicos ms relevantes. Posteriormente, se compara la
clasificacin con los datos de mapeo de sondaje y estudian las principales fuentes de
error.

1.5 Caractersticas relevantes del yacimiento Escondida

Escondida es un yacimiento tipo prfido Cu-Mo ubicado en el norte de Chile,


perteneciente al cluster Escondida junto a Escondida Norte, Zaldivar y Chimborazo.
Estos yacimientos se encuentran dentro del dominio del Sistema de Fallas de Domeyko
(SFD) en la franja metalognica del Eoceno superior Oligoceno inferior (Ver Figura
4
1.1). Dentro de esta franja se encuentran los depsitos de Cu con mayor importancia
econmica para Chile, adems de Escondida, tales como: Chuquicamata, Collahuasi y
El Salvador (Maksaev, 2001).

Figura 1.1. Franja metalognica del Eoceno Superior Oligoceno Inferior. Se


esquematiza el Sistema de Fallas de Domeyko (Maksaev, 2001).

Minera Escondida Limitada (MEL) pertenece en un 57.5% a BHP Billiton Ltds., 30%
a Rio Tinto Plc, 10% al consorcio corporativo JECO el cual incluye Mitsubishi, Nippon
Mining and Metals Ltd., y 2.5% a International Finance Corp. Minera Escondida es parte
del CSG (Customer Sector Group) de metales base de BHP Billiton, y cuanta con dos
faenas principales a rajo abierto: Escondida y Escondida Norte (Camus, 2003; BHP
Billiton, 2012).

En el ao fiscal del 2012, Escondida produjo 333.8 kt de concentrado de cobre y


172 kt de ctodos de cobre, lo cual es inferior a la produccin de 2010 y 2011. Segn
5
los estndares de produccin actual, Escondida tiene una vida estimada de 54 aos
(BHP Billiton, 2012).

Tabla 1.1. Produccin anual de cobre en miles de toneladas en Escondida (BHP


Billiton, 2012).

2012 2011 2010


Concentrado de cobre [kt] 333.8 390.5 448.1
Ctodos de cobre [kt] 172 179.1 174.2

Tabla 1.2. Produccin anual de oro y plata en miles de onzas en Escondida (BHP
Billiton, 2012).

2012 2011 2010


Oro 50.9 84.7 76.4
Plata 1921 2849 2874

Tabla 1.3. Recursos minerales de cobre al 30 de junio de 2012 en Escondida (BHP


Billiton, 2012).

Medidos Indicados
Tipo de mena Mt %TCu Mt %TCu
xido 104 0.81 58 0.79
Mixto 97 0.82 74 0.65
Sulfuro 3855 0.71 3081 0.55

Tabla 1.4. Reservas de mena de cobre al 30 de junio de 2012 en Escondida (BHP


Billiton, 2012).

Probado Probable
Tipo de mena Mt %TCu Mt %TCu
xido 76 0.9 40 0.88
Sulfuro 2779 0.8 2148 0.58
Lixiviacin de Sulfuro 1150 0.5 827 0.44

6
En sus inicios, la explotacin del yacimiento Escondida se enfoc en la zona de
enriquecimiento suprgeno, caracterizada por altas leyes de cobre. Actualmente,
Escondida se encuentra en una etapa de transicin hacia la zona de sulfuros primarios
en donde la eficiencia del procesamiento mineral mediante flotacin juega un papel
fundamental en la optimizacin de las ganancias en el negocio minero. El proyecto
MQAlt busca crear una herramienta para la caracterizacin mineralgica de la ganga
con el objetivo de predecir su comportamiento en el proceso de flotacin.

2 Marco geolgico del yacimiento Escondida

El yacimiento de Escondida se caracteriza por una secuencia de rocas volcnicas


y sedimentarias intruidas por varios cuerpos magmticos responsables de la
mineralizacin econmica (ver Figura 2.1).

7
Figura 2.1. Mapa geolgico del yacimiento Escondida (modificado de Richards et
al., 2001).

2.1 Rocas estratificadas

Las rocas estratificadas que afloran en el yacimiento Escondida corresponden a


4 formaciones (ver Tabla 2.1). Las caractersticas litolgicas de cada una de ellas se
describen a continuacin.

Tabla 2.1. Edades de Formaciones de rocas estratificadas.

Formacin Edad [Ma] Edad asignada Referencia


Augusta 71.72.9 a Paleoceno - Eoceno
Marinovic et al., 1995
Victoria 41.22.2 inferior
Gallego &
El Profeta y Trisico Superior -
- Covacevich, 1998 y
Santa Ana Jursico Inferior
Cornejo, 2009
1967 a
Agua Dulce Trisico Superior Marinovic et al., 1995
1696
La Tabla 300 a 260 Carbonfero -Prmico Cornejo, 2009

2.1.1 Formacin La Tabla (Carbonfero-Prmico)

Rocas gneas extrusivas e intrusivas intermedias a flsicas correspondientes al


basamento cristalino de Escondida. Solo aflora en el rajo de Escondida Norte (Richards
et al., 2001; Cornejo et al., 2009).

8
2.1.2 Formacin Agua Dulce (Trisico Superior)

Aflora al NE del rajo Escondida Norte y est compuesta por rocas andesticas
con textura porfrica de plagioclasas sausuritizadas inmersas en una masa fundamental
fina con ocasional textura traqutica. Los minerales mficos corresponden
principalmente a hornblenda y clinopiroxeno. Esta roca se encuentra comnmente
afectada por alteracin propiltica.

Las rocas riolticas tienen textura porfrica con fenocristales de cuarzo, feldespatos
y menor biotita inmersos en una masa fundamental vtrea devitrificada, pudiendo
presentar textura eutaxtica y bandeamiento de la masa fundamental (Richards et al.,
2001).

2.1.3 Formacin El Profeta y Santa Ana (Jursico-Cretsico)

Rocas sedimentarias con intercalaciones de rocas volcnicas que afloran al SW


de Escondida. En las zonas donde aflora, El Profeta subyace a la Formacin Agua
Dulce (Richards et al., 2001).

La Formacin El Profeta posee un espesor de alrededor de 1000 m de calizas


fosilferas, calizas brechizada y lentes de arenisca y limo. Puede presentar
intercalaciones con lavas andesticas, la cual presenta fenocristales de plagioclasas,
piroxenos y biotita inmersos en una masa fundamental vtrea devitrificada.

La Formacin Santa Ana se caracteriza por areniscas y limo con horizontes de


carbonato, intercalado con capas de andesitas y riolitas. El nivel inferior marino
corresponde a intercalaciones de areniscas lticas y arcsicas de grano medio a grueso,
presentando bancos de calcarenitas fosilferas finamente estratificadas. El nivel superior
corresponde a intercalaciones de areniscas rojas grano medio, volcarenitas verdes,
conglomerados rojos de grano fino, calizas margosas laminadas, calci-ruditas rojas e
intercalaciones menores de lutitas, conglomerados y brechas con intercalaciones de
limolitas (Richards et al., 2001; Cornejo et al., 2009).

9
2.1.4 Formacin Augusta Victoria (Paleoceno-Eoceno inferior)

Se distribuye principalmente en el sector SW del rajo Escondida y corresponde a


la roca de caja del complejo subvolcnico mineralizador del yacimiento Escondida.

La roca corresponde a lavas andesiticas de textura porfrica con alto contenido en


fenocristales de plagioclasa y clinopiroxeno, inmersos en una masa fundamental
afantica de feldespatos, biotita y menor hornblenda con alteracin a clorita y epidota.
Los fenocristales pueden alcanzar un tamao de hasta 2 mm de largo, y presentar
sericitizacin y argilizacin (Richards et al., 2001; Vliz, 2004).

2.2 Rocas intrusivas

2.2.1 Complejo Intrusivo Feldesptico Escondida (CIFE)

Conformada por 3 fases: Prfido Colorado Grande, Prfido Escondida y Prfido


Rioltico (Camus, 2003).

El CIFE es un stock cuarzo monzontico-granodiortico que intruye a la Formacin


Augusta Victoria y posee una forma elptica con eje mayor de largo mximo de 4.5 km
en direccin N30W a N40W, y eje menor de 2.5 km mximo (Padilla-Garza et al., 2004;
Padilla et al., 2001).

Las 2 primeras fases corresponden a prfidos synmineralizacin cuya mineraloga


se caracteriza por abundantes fenocristales de feldespato con cuarzo subordinado y
fenocristales de biotita en una matriz fina cuarzo-feldesptica. Presentan biotita
secundaria poco desarrollada y, clorita y sericita como alteracin selectiva de biotita y
feldespatos, respectivamente (Richards et al., 2001).

La intrusin Colorado Grande tiene una edad de 37.91.1 Ma (U-Pb en circones,


Richards et al., 1999).

10
La tercera intrusin del CIFE corresponde al prfido-brecha, la cual tiene una
geometra elongada en direccin N10W con eje mayor de hasta 1 km y eje mejor de
250 m. Su composicin es similar a la de la intrusin Escondida pero incluye fragmentos
mineralizados de los intrusivos previos y la andesita. El contenido en fragmentos es de
un 8% en volumen pudiendo encontrarse sectores donde llegara hasta un 60% donde
exhibe una textura de brecha intrusiva (Padilla et al., 2001).

2.2.2 Prfido Rioltico

El Prfido Rioltico es un cuerpo intrusivo hipabisal que se distribuye en la zona


NE y SE del rajo Escondida. Se encuentra en contacto por intrusin con el Prfido
Escondida y la Formacin Augusta Victoria, y tiene una edad de emplazamiento de
34.71.7 Ma. Posee geometra de domo y sera posterior a la mineralizacin del
yacimiento. Petrogrficamente se caracteriza por fenocristales de cuarzo ms
abundante que en los intrusivos synmineralizacin. Tambin presenta fenocristales de
feldespatos y biotita relicta en la matriz de grano fino. Esta intrusin se asocia a los
eventos de silicificacin y alteracin arglica avanzada (Richards et al., 2001; Vliz
2004).

Posteriormente, diques porfricos dacticos y cuarzo-latticos intruyeron los


prfidos anteriores siguiendo la tendencia NW de las fallas postmineralizacin,
produciendo una leve alteracin serictica de edad K-Ar de 31 2.8 Ma (Alpers y
Brimhall, 1988). Estos diques indican un lmite temporal de la actividad magmtica
(Padilla et al., 2001).

2.2.3 Stocks diorticos y diques

Stock diorticos y diques se encuentran intruyendo a la Formacin Augusta


Victoria y el Profeta. Se caracterizan por presentar abundantes fenocristales de

11
hornblenda, y magnetita lo que indica un alto contenido de agua magmtica y alto
estado de oxidacin (Richards et al., 2001).

Finalmente, diques dacticos y cuarzo-latticos con tendencia NW se emplazaron


en fallas postmineralizacin (Richards et al., 1999).

2.3 Brechas

Adems del prfido brecha asociado al CIFE, existen otros 3 tipos de brecha en
Escondida: Brecha de enfriamiento formada en el contacto del domo rioltico y la
andesita, diques de guijarros mineralizados y estriles que postdatan el emplazamiento
de la riolita, y brecha tectnica asociada a los sistemas de fallas (Padilla et al., 2001).

2.4 Gravas

Secuencia sedimentaria continental compuesta por areniscas y brechas


sedimentarias. Textura poco consolidada y mal seleccionada. Los fragmentos son
angulosos a subangulosos, de tamaos variables de milimtricos a decimtricos.
Presenta una matriz compuesta por material detrtico tamao arena, y cemento de yeso,
carbonatos y sales indeterminadas.

Su disposicin es subhorizontal sobreyaciendo en contacto por inconformidad a


todas las unidades anteriores. Presenta zonas con mineralizacin extica de cobre
(Vliz, 2004).

2.5 Estructuras

El yacimiento posee dos sistemas estructurales principales: NNW y NE (ver


Figura 2.2). Los altos ngulos de manteo sugieren que ambos sistemas se asocian al
SFD.
12
El sistema NNW se asocia a la etapa de mineralizacin del yacimiento. En el rajo
Escondida las fallas principales con esta orientacin corresponden a la falla Panadero,
Ortiz y Portezuelo. Las fallas Zaldvar y Ferrocarril se encuentran fuera del pit principal.
Estas estructuras controlan la distribucin de las alteraciones cuarzo-sericita, arglica
avanzada, y las vetillas polimetlicas.

El rumbo de estas estructuras es de N10W con variaciones entre N20E y N20W.


Las fallas de la zona W del rajo mantean hacia el E, mientras que las del E lo hacen en
direccin contraria, lo que le da una geometra tipo graben. La continuidad de las fallas
es de hasta 1 km en andesita, disminuyendo en los intrusivos. La mayora de los diques
de guijarro y diques estriles cuarzo-mozonticos posee una tendencia con igual
orientacin. En particular, la falla Ortiz en el sector E, en el contacto entre el stock y la
andesita se forma una geometra rmbica, la cual en perfil presenta una geometra del
tipo flor positiva (Padilla et al., 2001).

El sistema NE se interpreta como postmineralizacin debido a que corta a los


intrusivos y la mineralizacin. Las fallas son de carcter dextral a lo largo de rumbo NE
presentando un manteo entre 60 y 70 S. Presentan espesor promedio de 0.3 m, de
continuidad lateral de hasta 2 km (Richards et al., 2001; Padilla et al., 2001; Romero et
al., 2010).

13
14
Figura 2.2. Esquemas de las estructuras principales de Escondida. A) Geologa y
estructuras de Escondida. B) Geologa y estructuras del distrito Escondida en general
(Padilla et al., 2001).
2.6 Alteracin hidrotermal y mineralizacin

La alteracin hidrotermal de Escondida se asocia a la intrusin del prfido cuarzo


monzontico a granodiortico del Eoceno Oligoceno emplazado en lavas andesticas
Paleocenas. La distribucin espacial de las alteraciones corresponde a una zonacin
tpica de un sistema tipo prfido cuprfero (ver Figura 2.3).

Figura 2.3. Mapa de alteraciones hidrotermales del yacimiento Escondida


(Modificado de Padilla et al., 2001).

La mineralizacin hipgena de Cu se asocia a la intrusin del prfido dactico


cuarzo-feldespato-biotita, conocido como Prfido Escondida, el cual incluye al Prfido
Colorado Grande. Este intrusivo se habra emplazado a los 37.91.1 Ma siendo
posteriormente intruido por el prfido Rioltico (edad U-Pb de 34.7 1.7 Ma). Durante
este periodo ocurrira el principal evento mineralizador de Escondida (Richards et al.,
1999; Richards et al., 2001).

15
A los prfidos del CIFE asociados a mineralizacin hipgena de Cu-Fe se les
denomina intrusivos synmineralizacin. Estos se caracterizan por presentar una
alteracin potsica con sobreimposicin flica. Por otro lado, el Prfido Rioltico se
denomina intrusivo postmineralizacin, ya que no introduce mayor mineralizacin al
sistema, sino que solo habra removilizado y redistribuido los sulfuros anteriores. A este
se asocia a los eventos de silicificacin y alteracin arglica avanzada, presentando
abundante mineralizacin local de sulfuros en vetillas, pero escasa a escala de roca
total (Richards et al., 1999; Richards et al., 2001).

Se distinguen tres etapas (ver Tabla 2.2) de alteracin hidrotermal sobreimpuestas


caracterizadas por sus asociaciones minerales, distribucin espacial, temporalidad y
mineralizacin asociada (Padilla et al., 2001).

Tabla 2.2. Etapas de alteracin hidrotermal y vetillas asociadas.

Etapa Alteracin hidrotermal Vetillas


Potsica
1 Propiltica Tempranas: 1a, 1b y 1c
Silicificacin
Clorita-sericitacuarzo 2a
2
Cuarzo-sericita 2b y 2c
3 Arglica avanzada 3

2.6.1 Alteracin potsica

Presenta un control litolgico, es decir, el tipo de asociacin mineral depende de


la roca a la cual afecta. En el stock Escondida, la alteracin potsica se caracteriza por
feldespato potsico como reemplazo de plagioclasas en cristales y masa fundamental.
En andesita, se caracteriza por biotita secundaria, menor proporcin de feldespato
potsico y anhidrita diseminadas y en vetillas. La biotita es fina y ocurre como minerales
16
en la masa fundamental y fenocristales de plagioclasa y mficos. Los cristales de
plagioclasa de mayor tamao son alterados parcialmente.

La mineralizacin hipgena asociada corresponde a magnetita, bornita y


calcopirita asociada a biotita y feldespato potsico-cuarzo. En andesita con alteracin
pervasiva de biotita, predomina magnetita, la cual puede ser reemplazada por
calcopirita-bornita. En el stock Escondida, predomina calcopirita por sobre la magnetita.

2.6.2 Alteracin propiltica

Se forma como halo en torno al ncleo potsico presentndose como una


alteracin pervasiva pudiendo destruir total o parcialmente la mineraloga y textura
original de la andesita.

Se define por la presencia de epidota principalmente, y menor proporcin de


clorita, montmorillonita, biotita, grosularia y vetillas de carbonatos y ceolitas. En
andesita, esta asociacin mineral afecta a anfboles y plagioclasas. Las plaglioclasas
pueden ser ocasionalmente reemplazadas por grosularia y presentar un halo albtico.

La mineralizacin de sulfuros se caracteriza por abundante pirita, ocasionalmente


asociada con calcopirita en vetillas.

2.6.3 Silicificacin

Capa de silicificacin en el contacto entre el stock Escondida y la andesita. Es


interpretada como un contacto por metasomatismo que afecta a ambas litologas. Esta
alteracin se define por una capa de cuarzo microcristalino, el cual le otorga una textura
afantica y mayor resistencia a la roca.

17
2.6.4 Vetillas Tempranas

Estas vetillas se componen por cuarzo o cuarzo-ortoclasabiotitaanhidrita. Se


subdividen en tres subgrupos definidos por su temporalidad.

1a: vetillas sinuosas, regulares y discontinuas, compuestas por cuarzo sin halo
de alteracin, principalmente distribuidas en el intrusivo o en andesita cercana al
contacto con este. De largo inferior a 10 cm y espesor entre 1 mm y 2 cm.

1b: vetillas de cuarzo-feldespato potsico con menor proporcin de biotita y


anhidrita. Ms continuas y menos sinuosas que 1a. Presenta halo de feldespato
potsico haciendo los mrgenes irregulares y poco definidos. Espesor variable
entre 1 mm y 3 cm.

1c: vetillas continas de cuarzo sin halo, principalmente distribuidas en el


intrusivo. Mrgenes irregulares. Largo hasta 70 cm y espesor de 1 a 5 cm.

2.6.5 Clorita-sericitacuarzo

Vara segn predominio de clorita o sericita. Puede ocurrir diseminada y en


vetillas afectando a andesita con alteracin biottica. En el stock Escondida, se
caracteriza por la alteracin selectiva de biotita y mficos, y feldespato potsico de la
alteracin potsica alterados a clorita-sericita y menores proporciones de clorita en
vetillas reabiertas o fracturas. Esta alteracin no destruye la textura previa de la roca a
excepcin de las vetillas y en las proximidades de sus respectivos halos.

La relacin sericita/clorita disminuye al aumentar la distancia desde el intrusivo.


Se le asocia el tipo de vetilla 2a.

2a: vetillas que cortan o reabren vetillas tempranas. Largo variable de pocos
centmetros hasta 50 cm con espesor variable entre 1 y 3 cm. Son ms continuas que
las tempranas y presenta bordes irregulares. Son sinuosas y pueden presentar
ramificaciones. En andesita, poseen una zonacin conformada por sulfuros en el centro

18
y bajo contenido de cuarzo, un halo de sericita y clorita subordinada, lo cual cambia
hacia el exterior, llegando a clorita y sericita subordinada.

2.6.6 Cuarzo-sericita

Alteracin presente en el stock Escondida y en andesita cercana al contacto con


este. Ocurre principalmente en vetillas y halos de alteracin. Posee una alta densidad
de vetillas 2b y 2c.

Las vetillas 2b y 2c poseen un largo superior a 20 cm y, en ocasiones, son


mayores a 1 m. Su espesor vara entre 1 y 15 cm con bordes irregulares presentando
alteracin cuarzo-sericita intensa a moderada y dbil. La mineraloga de las vetillas
corresponde a cuarzo y sulfuros con halo serictico pervasivo.

La mineralizacin se asocia a las vetillas intercrecida con la alteracin serictica o


clorita-sericita, y corresponde a calcopirita con menor proporcin de molibdenita y pirita.
En general, el contenido de calcopirita y pirita se encuentran en razn 3:1.

2.6.7 Arglica avanzada

Esta alteracin afecta a la andesita, stock Escondida y riolita. Su distribucin es


controlada por fallas y por el contacto entre riolita y roca caja. La asociacin mineral que
la define corresponde a pirofilita, alunita y cuarzo, adems de vetillas de cuarzo con
abundantes sulfuros. En el contacto entre el domo rioltico y andesita, la alteracin es
pervasiva y se caracteriza por texturas moteada como resultado de la asociacin
pirofilita-sericita.

Se le asocian vetillas de espesor variable entre pocos milmetros hasta 1 o 2


metros con textura brechosa, bandeada y masiva. Las vetillas se caracterizan por la
presencia de cuarzo, abundantes sulfuros y menor alunita. En los mrgenes de estas se

19
presentan bandas de cuarzo azucarado de grano fino y cuarzo fino semitraslcido con
espesor entre variable desde milmetros a 10 cm.

La mineralizacin asociada a esta etapa consiste en una amplia variedad de


sulfuros, entre los cuales se encuentra: calcopirita, bornita, pirita, covelina, calcosina,
enargita, esfalerita, galena y tenantita. Esta mineralizacin ocurre principalmente en
vetillas asociadas a la alteracin o anteriores. La pirita ocurre tanto en diseminacin
como en vetillas, constituyendo el principal relleno de estas. En general, la
mineralizacin de pirita y calcopirita se encuentran en una razn 5:1.

3 Marco Terico

3.1 Alteracin hidrotermal

Una alteracin hidrotermal es el producto de un proceso metasomtico donde la


roca sufre un cambio qumico, mineralgico y morfolgico. Este proceso se debe al
intercambio qumico ocurrido durante la interaccin de un fluido hidrotermal y la roca.

Los cambios producto del proceso de alteracin imprimen asociaciones minerales


que permiten obtener informacin acerca de sus condiciones de formacin como:
composicin del fluido hidrotermal, composicin de la roca, temperatura, pH, Eh y razn
agua/roca.

En la Figura 3.1 se muestra el modelo propuesto por Corbett & Leach (1997)
donde las asociaciones minerales para cada alteracin hidrotermal se ordenan relativas
al pH, temperatura y profundidad de formacin. En este diagrama se visualizan las
posibles alteraciones hidrotermales que pueden encontrarse en un sistema tipo prfido
cuprfero, asumiendo que un fluido inicial, a alta temperatura y de pH neutro alcalino
que produce una alteracin potsica en la roca, evoluciona hacia uno de menor
temperatura y menor pH produciendo una alteracin flica, arglica y arglica avanzada.

20
Figura 3.1. Asociaciones de minerales de alteracin comn en sistemas
hidrotermales (Corbett & Leach, 1997).

En un prfido cuprfero, las alteraciones hidrotermales presentan una zonacin


espacial asociada y centrada en un ncleo. De base a techo, se encuentran las
alteraciones calcosdicas, potsicas, clorita-sericita (SCC), flica y arglica avanzada, y
en las zonas ms distales se encontraran las alteraciones clortica y propiltica (ver
Figura 3.2).

21
Figura 3.2. Modelo de zonacin de alteracin y mineralizacin para prfidos
cuprferos (Sillitoe, 2010).

Las alteraciones hidrotermales son progresivamente ms jvenes de base a


techo, por lo que inevitablemente, las ms someras afectan a las profundas. Es decir, la
alteracin potsica presentar un sobreimposicin clorita-sericita, la que a su vez
presenta una sobreimposicin flica a medida que disminuye la profundidad con
respecto al ncleo. En la Figura 3.3 se muestra la evolucin de la alteracin y
mineralizacin en funcin del aumento de la acidez del fluido hidrotermal producto de la
disminucin de la temperatura. Adems, el aumento en el estado de sulfidizacin del
fluido implica un cambio en las asociaciones de sulfuros, yendo desde calcopirita-
bornita, pasando por calcopirita-pirita y bornita-pirita, hasta pirita-enargita o pirita-
covelina.

22
Figura 3.3. Esquema generalizado de la secuencia de mineralizacin y alteracin
en un prfido cuprfero en funcin de la paleoprofundidad y tiempo de formacin
(Sillitoe, 2010).

En la base de datos geoqumica de Escondida, existen cuatro grupos de


alteraciones hidrotermales: las alteraciones potsicas (30-K1 y 31-K2), clorita-sericita
(40-SCC1 y 41-SCC2), flicas (50-S, 51-S1 y 52-S2) y arglica avanzada suprgena (61-
AA). Las alteraciones potsicas se asocian al ncleo y su mineraloga se caracteriza por
la presencia de biotita y feldespato potsico. A medida que la roca se aleja del ncleo,
se encuentra la alteracin SCC o sericita-clorita-arcillas en donde existe un predominio
de la clorita por sobre la sericita, el cual se invierte al alcanzar la zona de alteracin
flica caracterizada por cuarzo, sericita y pirita. La alteracin arglica avanzada
suprgena, como su nombre lo indica, se asocia a procesos de carcter suprgeno que
afectan a la roca luego de formada la alteracin flica.

Dada la evolucin terica de un sistema tipo prfido cuprfero, es de esperar que


las alteraciones hidrolticas (flica y arglica avanzada) presenten una mayor dispersin
en los datos geoqumicos, ya que se sobreimponen a alteraciones preexistentes.

23
3.2 Anlisis ICP-AES

El ICP-AES (Inductively Coupled Plasma Atomic Emission Spectrometry) es una


de las tcnicas ms utilizadas para el anlisis de elementos qumicos. Su aplicacin
depende directamente de la preparacin de las muestras previa al anlisis. En el caso
de las rocas, la disolucin de la muestra puede ser parcial o total dependiendo del tipo
de reactivo que se utilice.

La base de datos geoqumicos de Escondida fue generada en ICP-AES utilizando


agua regia como reactivo para disolver las muestras. El agua regia disuelve
parcialmente la roca dado que no es capaz de romper el tetraedro de silicio, el cual
constituye la estructura fundamental de los silicatos, a diferencia de otras tcnicas de
disolucin total como la digestin por cido fluorhdrico o digestin por fusin. Es decir,
el agua regia resalta la presencia de algunos grupos minerales como los sulfuros y
carbonatos por sobre los silicatos.

Una de las ventajas es que elementos como el As, Hg, S y Sb se volatilizan


fcilmente al utilizar cido fluorhdrico (HF), mientras que el agua regia los retiene.
Adems, la concentracin de slidos disueltos totales es menor al usar agua regia. Sin
embargo, la ventaja ms importante para los objetivos de este estudio es que la
disolucin por agua regia solo afecta a ciertas fases minerales, lo que disminuye
considerablemente el nmero de variables ya que se acotan las posibles fuentes de los
elementos.

3.2.1 Disolucin por Agua Regia

El agua regia consiste en una mezcla de HCl y HNO3 en razn 3:1 en volumen.
Segn procedimientos estndar, las muestras son tratadas por un perodo entre 1 y 3
horas a una temperatura entre 80 y 95 C (Gill, 1997). Sin embargo, el protocolo de
anlisis varia para cada laboratorio, y dado que se desconocen, no es posible
detallarlos.

24
El mtodo de disolucin por agua regia es parcial ya que no afecta a la mayora de
los silicatos. Sin embargo, la digestin es cercana a la total para algunos minerales
solubles en agua, elementos intercambiables, carbonatos, sulfuros, arsenuros,
selenuros, teluros, oxi-hidrxidos de Fe y Mn, sales, y algunas micas, sulfatos(a
excepcin de la barita) y aluminosilicatos como las arcillas (Chao, 1984; Chao &
Sanzalone, 1992; Church et al., 1987; Hall & Cameron, 2001; Klassen, 2001; Hall, 1999;
Xu et al., 2012).

Los grupos estructurales que contienen sitios con altos nmeros de coordinacin
acomodan a los cationes de mayor radio inico con mayor facilidad. Sin embargo, estos
se encuentran dbilmente enlazados a la estructura cristalina del mineral, por lo que
tienden a entrar en solucin fcilmente, como es el caso de los elementos traza (Church
et al., 1987).

El agua regia tiende a disolver metales presentes en los silicatos que se


encuentren en coordinacin VI o superior (VIII, IX, X y XII), lo que es comn en
inosilicatos y filosilicatos mficos.

Church et al., (1987) realiza un estudio de disolucin por cido oxlico y agua
regia en distintos grupos minerales, donde se comparan los resultados obtenidos con
los de digestin total por HF-HNO3-HClO4. Los resultados para la disolucin por agua
regia muestran que es cercana a la total para el caso de los sulfuros, carbonatos,
algunos xidos de Fe, Mg, Al y Mn, y algunas micas y sulfatos (Ver Figura 3.4).

25
Figura 3.4. Porcentaje de minerales disueltos segn el procedimiento consecutivo
de cido oxlico y agua regia. El grosor de la barra representa el nmero de
muestras utilizada para el estudio (Church et al., 1987).

Los feldespatos, espinel, cromita y muscovita prcticamente no fueron atacados


por el agua regia. Sin embargo, otros minerales como el olivino posee recuperaciones
mayores al 90% para Fe, Mg y elementos traza en posicin octadrica, al contrario del
dipsido y hornblenda que presentan porcentajes bajo el 50%. Mientras que los
resultados para magnetita y biotita, muestran que su disolucin es cercana a la total
para los siguientes grupos de elementos Fe-Mg-Mn-V-Cu-Zn-Ti-Ba-Ce y K-Fe-Mg-Mn-
Al-Ti-Cr-Li-Cu-Zn-Ni-Co-V-Ba-Y-Ce-P-Be, respectivamente (Ver Figura 3.5)

Se recalca la diferencia entre los porcentajes de recuperacin entre muscovita y


biotita, ya que a pesar de que ambos son filosilicatos, difieren en su estructura. La
biotita es trioctadrica y la muscovita es dioctadrica.

26
Figura 3.5. Porcentaje de disolucin de varios minerales, obtenidos al comparar
digestin por agua regia relativa la digestin por HF-HNO3-HClO4 (Church et al.,
1987).

Snll & Liljefors (2000) realizaron un estudio similar para elementos mayores
presentes en los minerales de depsitos de tills. Se analizaron las solubilidades de los
minerales en agua regia y se compar con la disolucin total de la muestra en HF y
HNO3 caliente. Adems, se estudi la dependencia con la granulometra usando
tamaos de <2mm, <63 m y polvo (<10-20 m), concluyendo que a medida que la
fraccin analizada es ms fina, la recuperacin de los elementos aumenta.

27
A tamaos del orden de micrones (m), las partculas tienen importantes
propiedades de superficie. Estas propiedades tienen consecuencias como un aumento
considerable de la solubilidad y un desbalance de la carga superficial lo que las
convierte en potenciales sorbentes para especies disueltas. La importancia de las
propiedades de superficie aumentan proporcionalmente con el rea del material y a su
densidad de carga superficial o nmero de sitios con carga por unidad de superficie.

Las arcillas poseen un alta razn de rea/volumen, debido a su tamao de orden


de micrones, lo que favorece las reacciones con soluciones acuosas. A medida que
disminuye el tamao de las arcillas, la ionizacin de la superficie aumenta debido al
desbalance de cargas ocasionado por la ruptura de enlaces que exponen O 2- y OH-.
Otras razones de la carga neta negativa son las sustituciones isomrfas e
imperfecciones de la red cristalinta. Esta propiedad hace que algunas arcillas
reaccionen rpidamente con las soluciones con las que interactan, hacindolas
solubles en agua regia (Langmuir, 1997).

A partir de los resultados obtenidos por Snll & Liljefors (2000), se concluye que la
solubilidad de algunos elementos mayores es cercana a la total en las micas
trioctadricas: biotita y clorita. Estos elementos corresponderan a K, Mg, Fe, Ti y Al
para biotita y Mg, Fe y Al para clorita. Cabe destacar la solubilidad de Ca y Ti en titanita,
y de Fe y Ti en ilmenita, las cuales son cercanas al 60% en ambos casos. Para todos
los dems minerales, feldespatos, epidota, hornblenda y mica blanca, el porcentaje de
recuperacin de los elementos es inferior al 35% por lo que se consideran
prcticamente insolubles en agua regia. Los errores varan entre un 10% y un 15%
debido a variaciones composicionales de las muestras analizadas (ver Figura 3.6).

Estudios de disolucin en agua regia ms recientes, realizados en lutitas, indican


que disuelve casi la totalidad de Mg y Fe, mientras que solo una fraccin muy baja de Al
y K, debido a que se encuentran en fases con poca solubilidad como arcillas, micas y
feldespatos. Adems, concluye que las fases residuales no disueltas por agua regia
corresponden principalmente a silicatos, entre otros minerales (Xu et al., 2012).

28
a Aluminio Hierro

Muscovite Ilmenite
Biotite
Biotite
Fe-Chlorite
Mg-Chlorite Fe-Chlorite

Epidote Mg-Chlorite
Hornblende
Epidote
Plagioclase
Microcline Hornblende

0 20 40 60 80 100 0 20 40 60 80 100

Powder 30500 m <63 m <2 mm 0.52 mm Powder 30500 m <63 m <2 mm 0.52 mm

Titanio Magnesio

Anatase Hornblende

Rutile
Biotite

Ilmenite

Fe-Chlorite
Sphene

Biotite Mg-Chlorite

0 20 40 60 80 100 0 20 40 60 80 100

Powder <63 m <2 mm 0.52 mm Powder 30500 m <63 m <2 mm 0.52 mm

Potasio Calcio

Hornblende Sphene

Microcline Epidote

Muscovite Hornblende

Biotite Plagioclase

0 20 40 60 80 100 0 20 40 60 80 100

Powder <63 m <2 mm 0.52 mm Powder <63 m <2 mm

Figura 3.6. Porcentaje de recuperacin de elementos por mineral (datos de


Snll & Liljefors, 2000).

29
4 Anlisis de Datos

4.1 Anlisis geoqumicos y mapeo de sondajes

La base de datos consiste en anlisis qumicos de muestras de sondajes del


yacimiento Escondida. Los datos fueron obtenidos por anlisis mediante el mtodo
analtico ICP-AES por disolucin en agua regia de la muestra. Cada muestra
corresponde a un compsito de 14 y 16 metros consecutivos, de tal manera de tener
muestras puntuales cada 15 metros a lo largo del sondaje. Se cuenta con un total de
15408 datos georreferenciados para 53 elementos mayores y traza. Esta cantidad
corresponde a solo una fraccin de la base de datos total de ms de 35.000 muestras,
la cual fue procesada y depurada previamente en el proyecto denominado
Modelamiento Geoqumico Distrito Escondida: Discriminacin de Ambientes
Hidrotermales, Fase 2: Estudio de Distribucin Espacial y Determinacin de Criterios de
Discriminacin de Ambientes Hidrotermales en el Yacimiento Escondida (Townley et
al., 2012; Luca et al., 2012).

La base de datos fue compilada a partir de datos preexistentes y no fue generada


para este estudio en particular. Tampoco se cuenta con informacin acerca de los
distintos laboratorios que han producido los datos, sus respectivos protocolos y/o
control de ciertos parmetros tales como: tiempo de disolucin, pH, temperatura y
granulometra.

Los parmetros mencionados anteriormente controlan la estabilidad qumica de


las fases minerales en contacto con el agua regia. Si el tiempo de disolucin es muy
prolongado, las variaciones del pH pueden llevar a la formacin de coloides, los cuales
a su vez pueden actuar como sustrato para la sorcin o reprecipitacin de metales. La
temperatura influye en la cintica de la reaccin, ya que un aumento de esta implica un
aumento en la energa y por lo tanto, una mayor interaccin entre las partculas (Kelley
et al., 2003).

La base de datos geoqumicos se cruza con la informacin de logueo de sondajes


de tal manera de contar con caractersticas geolgicas como: litologa, alteracin, zona
mineral y abundancia relativa de minerales de alteracin. Como el procedimiento de
30
logueo de sondajes no considera longitudes constantes para la caracterizacin, se debe
realizar un compsito de la informacin y ajustarla a los tramos definidos para los
anlisis qumicos.

Este compsito de informacin categrica se realiza considerando las dos


longitudes que representen el mayor porcentaje del tramo definido por la base de datos
geoqumica. La suma de estas longitudes no necesariamente debe ser 100% del tramo,
pero se considerarn aquellas muestras cuya suma sea mayor o igual a un 93% con el
objetivo de disminuir la mezcla de informacin geoqumica. Por ejemplo, si el tramo
correspondiente a un anlisis qumico es de 14 m y el mapeo indica que posee 3 tipos
de alteracin abarcando 7, 6 y 1 m, respectivamente, solo se considerarn los 2
primeros dado que suman 93% de la longitud total (13 m /14 m).

4.2 Definicin del sistema

Con la finalidad de acotar las variables del sistema, se definen las caractersticas
litolgicas de las muestras con las que se realizar el anlisis (ver Anexo A). Como
criterio general, se utilizaron aquellas categoras de litologa, alteracin hidrotermal y
zona de mineralizacin que presentan mayor cantidad de datos como primer y segundo
mapeo por muestra (ver Figura 4.1).

31
Cantidad de datos por litologa
8000

7000

6000

5000

4000

3000

2000

1000

PFT

TR

RM
D
BP
PG
TD
AV
AN

BI

PNO

AB
PF
PC

BH
PT

BT

DT

PR
DR

RA
RS

PDA
GR

UVS
1 litologa 2 litologa

Figura 4.1. Cantidad de muestras por litologa segn primera y segunda litologa
mapeada por tramo (ver Anexo A).

A pesar de que la litologa PC (Prfido Cuarcfero - 33) presenta una gran cantidad
de muestras, solo posee alteraciones sericita-cuarzo y algunas para clorita-sericita, por
lo cual se deja fuera del anlisis (ver Figura 4.2).

32
Cantidad de datos de alteracin por litologa
2500

2000

1500

1000

500

PF PC AN

Figura 4.2. Cantidad de muestras por alteracin para las litologas PF (31), PC (33)
y AN (50) (Ver Anexo A).

Otro parmetro importante al momento de definir el sistema en el que se trabajar


es la mineralizacin. La mineralizacin predominante corresponde a la variable
denominada minzone y se agrupa en 3 categoras: alto enriquecido (HE), bajo
enriquecido (LE) y primario (PR) (ver Figura 4.3). La zona de HE se caracteriza por la
presencia de calcosina, covelina y pirita. La zona de LE se caracteriza por calcosina,
covelina, calcopirita, posible bornita y pirita. La zona PR se caracteriza por pirita,
calcopirita y/o bornita (ver Tabla 4.1).

Tabla 4.1. Codificacin de zonas mineralgicas.

Zona mineralgica Abrev. Cdigo


Cc + Py HE1 50
Cc + Cv + Py HE2 51
Cv + Py HE3 52
Cc + Cpy + Py LE1 60
Cc + Cv + Cpy + Py LE2 61
Cv + Cpy + Py LE3 62
Cv + Cpy + Bo (+/- Cc, +/- Py) LE4 63
Py PR1 72
Cpy - Py PR2 70
Bo - Cpy (+/- Py) PR3 71
33
Cantidad de datos por minzone
7000
6000
5000
4000
3000
2000
1000
0

CPCCPY
PR2

PR1
LX

SM

OXN
PR3

Plx

OXV

OXA
PlxE

LxH
CPOX

sk
HE1

HE2

HE3

CPMX

CPCu

POXV
LE2

LE3
LE1

MX

LE4

OXVN
1 Minzone 2 Minzone

Figura 4.3. Cantidad de muestras por zonas de mineralizacin o minzone (ver


Tabla 4.1 para codificacin).

Dado que las muestras pueden tener un mximo de dos alteraciones


predominantes, se considerarn solo aquellas cuyos tramos sumen ms de un 93%, es
decir, las 2 alteraciones representen casi la totalidad de la muestra. Por ejemplo, si el
tramo correspondiente a una muestra es de 15 metros y el mapeo indica que posee 3
alteraciones que abarcan 5 metros cada una, esta muestra es descartada porque el
porcentaje mximo que alcanzan a cubrir las alteraciones es del 66,6% solamente.

Considerando los criterios anteriormente descritos, se obtiene la base de datos para


43 elementos con un total de 9738 muestras ubicadas espacialmente en el yacimiento
Escondida. De estas, 8379 muestras corresponden a tramos de sondajes que
presentan solo una alteracin mapeada. Las 1359 muestras restantes presentan
mezcla de alteraciones, es decir, presentan dos alteraciones mapeadas por tramo (ver
Tabla 4.2).

34
Tabla 4.2. Cantidad de muestras con nica alteracin por tramo.

Alteracin hidrotermal Abrev. Cdigo Alteracin


Potsico (Slo Bt secundaria) K1 30 264
Feld K > Bt K2 31 122
Clorita-Sericita-Arcillas SCC1 40 1111
Clorita-Sericita-Cuarzo SCC2 41 1417
Flica S 50 1151
Sericita-Cuarzo S1 51 2551
Sericita-Cuarzo-Arcilla S2 52 1706
Argilizacin Suprgena AA 61 57
Total 8379

En la Tabla 4.3 se muestran los datos cuyos tramos presentan dos alteraciones
distintas. Las alteraciones etiquetadas como 1 alteracin corresponden a aquellas que
representan un mayor porcentaje del tramo en cuestin y no tienen ninguna
connotacin gentica. Se destaca que las alteraciones que presentan mayor cantidad
de muestras con 2 alteraciones son las tipo S, S1, S2, SCC1 y SCC2.

Tabla 4.3. Cantidad de datos por tramo con mezcla de alteraciones.

2 alteracin
K1 K2 SCC1 SCC2 S S1 S2 AA Total
K1 0 2 49 38 17 6 2 2 116
1 alteracin

K2 3 0 4 11 12 3 1 0 34
SCC1 43 1 0 17 122 5 30 19 237
SCC2 31 12 24 0 19 109 35 0 230
S 12 7 94 23 0 0 0 22 158
S1 7 2 10 120 0 0 169 4 312
S2 1 3 33 41 0 158 0 5 241
AA 1 0 5 1 16 4 4 0 31
1359

35
La alteracin S (cdigo 50) corresponde a un cdigo de mapeo antiguo el cual
caracteriza a la alteracin flica sin hacer distincin de la presencia de arcillas como lo
hacen la S1 y S2, es decir, posee informacin agrupada de la alteracin flica, por lo
que ser reclasificada al finalizar el anlisis.

De esta manera, la base de datos geoqumicos se reduce a 8240 muestras, 7225


para tramos con alteracin nico y 1015 para tramos con mezcla de 2 alteraciones. Esta
base de datos corresponde aproximadamente a un 23% de la base de datos original
compilada por Townley et al. (2012) y Luca et al.(2012), dado que solo se consideran
las muestras de la zona sur del rajo Escondida que cumplen con los criterios descritos
anteriormente.

La distribucin espacial de los datos muestra tendencias segn litologa, alteracin


hidrotermal y mineralizacin. Segn litologa (ver Figura 4.4) el prfido feldesptico se
distribuye preferentemente en la zona noreste del yacimiento pudiendo alcanzar mayor
profundidad que la andesita, la cual tiende a encontrarse en la zona oeste.

36
W E

Figura 4.4. Distribucin espacial de los datos geoqumicos segn la litologa.


Arriba: Perfil EW N = 107350 50 m. Abajo: Planta cota = 2750 50 m.

Segn la alteracin hidrotermal (ver Figura 4.5) es posible identificar tres zonas:
Oeste, Central y Este. En la zona Oeste, se encuentran preferentemente las
alteraciones SCC1 y SCC2, asociadas a la andesita. La alteracin K1 se distribuye en
la zona noroeste del yacimiento. En la zona Este, se encuentran las alteraciones S1 y
S2 asociadas al prfido feldesptico. La zona Central correspondera a una transicin
entre la Este y Oeste. En esta se encuentran las alteraciones SCC2 y K2 en contacto
en profundidad con las alteraciones S1 y S2. La alteracin AA es posible encontrarla en
superficie en cualquiera de las tres zonas descritas.

37
W E

Figura 4.5. Distribucin espacial de los datos geoqumicos segn alteracin


hidrotermal. Arriba: Perfil EW N = 107350 50 m. Abajo: Planta cota = 2750 50 m.

Al igual que la alteracin hidrotermal, es posible definir 3 zonas en funcin de la


distribucin espacial de la mineralizacin (ver Figura 4.6): alto enriquecido (HE), bajo
enriquecido (LE) y primario (PR). Estas zonas se caracterizan por presentarse
gradualmente a mayor profundidad. La zona de HE (50, 51 y 52) es de baja profundidad
y tiende a asociarse al prfido feldesptico en la zona Este. La zona de LE (60, 61, 62 y
63) se encuentran a profundidad intermedia y no presentan una distribucin preferencial
38
segn litologa o alteracin. Por ltimo, la zona PR (70, 71 y 72) se presenta a
profundidad alta y corresponde a la ms abundante ya que se caracteriza por la
presencia predominante de pirita. Al igual que LE, no presenta una tendencia a
asociarse con alguna litologa o alteracin en particular.

W E

Figura 4.6. Distribucin espacial de los datos geoqumicos segn mineralizacin.


Perfil EW N = 107350 50 m.

4.2.1 Porcentaje de disolucin

Los valores calculados de la disolucin muestran que en promedio se disuelve tan


solo un 7.5% de la roca, pudiendo llegar hasta 30% en algunos casos (ver Tabla 4.4).
Es decir, la suma de las concentraciones de elementos mayores no corresponde a un
100% como en el caso de los anlisis de roca total utilizando cidos ms fuertes como
el cido fluorhdrico.
39
Cabe sealar, que se realiza la siguiente distincin entre elementos mayores y
trazas. Los elementos considerados como mayores son aquellos que constituyen
minerales que conforman una roca dentro de un sistema en particular. Por otro lado, los
elementos traza son aquellos que pueden reemplazar a los elementos mayores en la
estructura cristalina sin afectar significativamente las fases. Adems, se destaca que
algunos elementos pueden comportarse como mayores o traza dependiendo el sistema
que se defina (Hofmann, 2003). Para este caso en particular, dada la mineraloga
soluble en agua regia, el grupo de los elementos mayores corresponde a: Al, Mg, K, Ca,
Na, Fe, S, Cu y Ti. Los elementos de P y Mn, pese a constituir parte de la mineraloga
accesoria, sus concentraciones son del orden de ppm por lo que solo se consideran al
analizar las variaciones conjuntas y no su valor absoluto (ver Anlisis de componentes
principales).

Tabla 4.4. Estadgrafos de porcentaje de disolucin calculada.

Estadgrafo % Total disuelto elementos mayores % Total disuelto


Media 7,54 7,73
Desviacin Estndar 3,24 3,31
Mnimo 1,88 1,98
Mximo 30,21 30,27

En la Figura 4.7 se muestran los boxplots de la disolucin total calculada para


elementos mayores por alteracin, evidenciando que el porcentaje de disolucin parece
no estar controlado por el tipo de alteracin hidrotermal asociada, ya que los rangos de
variacin se sobreponen casi totalmente.

40
Figura 4.7. Boxplot de disolucin total calculada para elementos mayores.

Los elementos representan un mayor porcentaje relativo de disolucin en orden


decreciente son Fe, S, Al, CuT, Mg, K, Ca, Na y Ti. Para el caso del Al, Mg y Ca, el
porcentaje de disolucin es mayor para las muestras caracterizadas por alteraciones
neutras (K1 30 a SCC2 41) y menor para las cidas (S1 51, S2 52 y AA 61),
mientras que el azufre se comporta de manera opuesta (ver Figura 4.8).

41
Figura 4.8. Porcentaje relativo promedio de elementos disuelto por alteracin.

La Figura 4.9 muestra que los elementos traza que presentan un mayor porcentaje
de disolucin corresponde a P, Zn y Mn, lo que puede deberse a la presencia de
minerales altamente solubles en agua regia como fosfatos, blenda y xidos de
manganeso.

42
% normalizado de disolucin de elementos traza
8.00%

7.00%

6.00%

5.00%

4.00%

3.00%

2.00%

1.00%

0.00%
Sr

Pb
Ba

Bi

Sn
As
As

Se

Ag

Sc
Sb

Be
Ga
Zn

Th

Te

Tl
Cr

Ni
P

V
Ce
Co
Rb

Cd

Cs

Re
U

In

Hg
Mn

Mo

La
Figura 4.9. Porcentaje normalizado promedio de elementos traza disueltos.

La diferencia entre los porcentajes de disolucin calculados por alteraciones se


debera al tipo digestin utilizada. El agua regia disuelve ciertos minerales que pueden
darse en distintas proporciones dependiendo de la alteracin hidrotermal. As, la
alteracin potsica biottica puede presentar mayor porcentaje de disolucin que la
potsica feldesptica, puesto que la biotita es disuelta fcilmente por el agua regia a
diferencia del feldespato y la sericita.

4.2.2 Abundancia de minerales

La cartilla de logueo de sondajes de Escondida cuenta con una seccin en la cual


se debe entregar un ndice de intensidad de presencia de algunos minerales. La escala
consiste en los valores nominales 0, 1, 2 y 3 que indican intensidad nula, dbil,
moderada e intensa, respectivamente. Esta seccin del mapeo de sondajes representa
la abundancia relativa entre los minerales de alteracin y en ningn caso la intensidad
de un tipo de alteracin en particular. Por ejemplo, la intensidad de slice puede

43
representar la generada por silicificacin, alteracin potsica, cuarzo-sericita, entre
otras.

Otros parmetros que se encuentran en esta base de datos son: Alteracin de los
mficos (AltMaf) y Destruccin de feldespatos (Desfel). Al igual que la intensidad
mineral, la alteracin de los feldespatos consta de 4 categoras: ausente (0), biotitizados
(1), cloritizados (2) y blanqueados (3). La categora de biotitizados y cloritizados son
asociados comnmente a las alteraciones potsica biottica y clortica, respectivamente.
La categora de blanqueados se asocia a las alteraciones flicas y arglicas suprgenas.

La destruccin de feldespatos consta de 5 categoras: ausente (0), nulo o fresco (1),


pseudomorfo (2), destruccin parcial (3) y destruccin total (4). La Figura 4.10
esquematiza el grado de destruccin de los feldespatos.

Figura 4.10. Esquema de grado de destruccin de los feldespatos.

44
Slice Clorita
100% 100%

80% 80%

60% 60%

40% 40%

20% 20%

0% 0%
30 31 40 41 51 52 61 30 31 40 41 51 52 61

Biotita Feldespato potsico


100% 100%

80% 80%

60% 60%

40% 40%

20% 20%

0% 0%
30 31 40 41 51 52 61 30 31 40 41 51 52 61

Sericita Arcillas
100% 100%

80% 80%

60% 60%

40% 40%

20% 20%

0% 0%
30 31 40 41 51 52 61 30 31 40 41 51 52 61

Figura 4.11. Abundancia relativa de minerales en cada alteracin hidrotermal.

45
En la Figura 4.11 se muestran las abundancias minerales mapeadas para un total
de 3906 muestras. Esta cantidad se debe a que el mapeo no es por tramos de longitud
constante y, por lo tanto, para hacerlos coincidir con la base de datos geoqumicos se
han considerado solo los mapeos de abundancia que representen la totalidad del
soporte geoqumico. En los grficos, la abundancia relativa de minerales por alteracin
se normaliza por la cantidad de datos.

La biotita y clorita presentan un mapeo de abundancia medio a alto


preferentemente en las alteraciones potsicas y clorita-sericita. Se destaca que la
presencia de clorita es esencial para distinguir entre las alteraciones SCC (40 y 41) y S
(51 y 52). La intensidad de la biotita es alta en la alteracin biottica. En la alteracin
arglica suprgena la biotita tiene un mapeo moderado a intenso lo que pone en duda la
clasificacin relativa de minerales. La biotita es un mineral inestable bajo las
condiciones fisicoqumicas de la alteracin arglica suprgena, lo que hace improbable
encontrar esta asociacin en la mayora de los datos de alteracin AA.

La sericita es mapeada con abundancia media a alta en las alteraciones hidrolticas


(flicas y arglica avanzada), mientras que solo con intensidad media a nula en las
alteraciones potsicas (30 y 31). Sin embargo, la intensidad media y baja de sericita en
las alteraciones potsica indica que existe una superposicin de las alteraciones
potsica y serictica a pesar de haber filtrado la base de datos por tramos con solo una
alteracin.

La arcilla es mapeada con abundancia baja a media en las alteraciones 40 y 52, las
cuales por definicin poseen arcillas. Mientras que en la alteracin arglica suprgena
se presenta con intensidad media a alta dado que es un mineral clave para definir esta
alteracin.

El slice es mapeado con intensidades relativamente homogneas en todas las


alteraciones, mientras que el feldespato solo se mapea en la alteracin potsica
feldesptica.

46
K1 K2
100% 100%
80% 80%
60% 60%
40% 40%
20% 20%
0% 0%

SCC1 SCC2
100% 100%
80% 80%
60% 60%
40% 40%
20% 20%
0% 0%

S1 S2
100% 100%
80% 80%
60% 60%
40% 40%
20% 20%
0% 0%

AA
100%
80%
60%
40%
20%
0%

Figura 4.12. Abundancia relativa de todos los minerales por


alteracin hidrotermal.

47
La alteracin potsica biottica (K1) se caracteriza por la destruccin de feldespatos
parcial a total; moderada a intensa biotita y clorita; moderada a baja de sericita y
arcillas; y mficos cloritizados y biotitizados, principalmente (ver Figura 4.12).

La alteracin potsica feldesptica (K2) se caracteriza por destruccin de


feldespatos pseudomorfa a parcial, moderada a baja sericita, biotita y clorita, moderado
a intenso feldespato potsico, nula arcillas y mficos cloritizados a blanqueados.

La alteracin clorita-sericita-arcillas (SCC1) se caracteriza por destruccin de los


feldespatos parcial a total; clorita, sericita y arcilla moderada a intensa; biotita moderada
a nula; nulo feldespato potsico; y mficos biotitizados y cloritizados.

La alteracin clorita-sericita-cuarzo (SCC2) se caracteriza por destruccin de


feldespatos parcial a total; clorita, sericita y slice moderada a intensa; nula a baja
presencia de arcillas, biotita y feldespato potsico; y mficos cloritizados.

La alteracin cuarzo-sericita (S1) se caracteriza por destruccin parcial a total de


feldespatos; nula clorita, biotita, feldespato potsico y arcillas; moderada a abundante
sericita y slice; y mficos blanqueados.

La alteracin cuarzo-sericita-arcillas (S2) se caracteriza por destruccin parcial a


total de los feldespatos; nula clorita, biotita y feldespato potsico; baja a moderada
slice, sericita y arcillas; y mficos blanqueados.

La alteracin arglica suprgena (AA) se caracteriza por una destruccin total a


parcial de los feldespatos; nula clorita y feldespato potsico; baja a moderada slice y
sericita; moderada a intensa arcilla y biotita; y mficos blanqueados.

4.3 Anlisis univariable

Los histogramas (ver Anexo A: Histogramas por alteracin para elementos mayores
y traza) muestran que la distribucin de los elementos se asemeja a una lognormal.
Adems, existen algunos elementos mayores y traza que presentan dos poblaciones
bien definidas y que coinciden con dos grupos de alteraciones hidrotermales:

48
alteraciones neutras y cidas. Las primeras la conforman las alteraciones potsicas (30
y 31) y las sericita-clorita (40 y 41), y en las segundas se encuentran las alteraciones
cuarzo-sericita (50, 51 y 52) y arglica suprgena (61).

Entre los elementos mayores destacan Mg, Al, Ca, P y Mn (ver Figura 4.13), y en
los elementos traza lo hacen Zn, Be, Ga, Sc, V e Y (ver Figura 4.14).

Figura 4.13. Histogramas de elementos mayores por alteracin hidrotermal.

49
Figura 4.14. Histogramas de algunos elementos separados por alteracin
hidrotermal.

Los grficos de caja (boxplot) se caracterizan por mostrar una caja, una lnea, un
punto y dos bigotes, uno inferior y otro superior. Dentro de la caja se encuentran los
valores entre el primer cuartil y el tercero (Q1 a Q3), la lnea corresponde a la mediana,
el punto es la media, el bigote inferior corresponde al 5% de los datos, y el bigote
superior al 95% de estos. Los datos graficados como crculos son aquellos que se
encuentran fuera de rango entre 5 y 95%.

Las alteraciones a pH neutro presentaran mayores concentraciones de Al, Mg, Ca,


Na, Ti, Mn, P, Zn, Ce, Rb, V, Y, Th y Be, mientras que las alteraciones cidas presentan
mayor concentracin de S (ver Figura 4.15).

50
Figura 4.15. Boxplots con eje logartmico de elementos por alteracin hidrotermal.

4.4 Depuracin de base de datos geoqumicos

Segn estudios previos realizados con la misma base de datos (Townley et al.,
2012; Luca et al., 2012), existen variables que fueron eliminadas debido a no cumplir
ciertas condiciones.

51
En este estudio, se procede a realizar un anlisis exploratorio de datos usando las
siguientes herramientas: Histogramas de frecuencia y acumulados, clculo de
estadgrafos bsicos y medidas de dispersin, y matriz de correlacin.

Los histogramas de frecuencia (ver Anexo A) evidencian que existen ciertos


elementos cuyas concentraciones son casi constantes, lo que se puede observar
tambin en los estadgrafos de rango y rango intercuartil cercanos a cero. De esta
manera, se descartan los siguientes elementos: Au, B, Hf, Hg, Nb, Ta y Zr. Adems,
algunos elementos, Li y Au, poseen ms del 20% de sus datos en blanco, lo que quiere
decir que no se posee informacin y son eliminados (ver Anexo ).

De este modo, la base de datos geoqumicos final consta de un total de 45


variables. Cabe destacar que existen elementos que poseen ms de una columna como
Cu_ppm, CuS y CuT y As_ppm y As_Tot. Con lo cual se contara con informacin de 42
elementos por cada muestra. En el caso del Cu, se decide trabajar con el CuT porque
representa el contenido total de cobre en la muestra y es uno de los principales
elementos presente en los minerales solubles en agua regia, mientras que para el As se
considera el soluble (As ppm) ya que el total no est disponible para todas las
muestras.

Comnmente, otro factor importante a considerar para disminuir el nmero de


variables es la correlacin entre elementos. Una correlacin muy alta indicara que hay
elementos que presentan una relacin lineal entre ellos y estaran entregando
informacin redundante. Los pares de elementos con mayor correlacin se muestran en
la Tabla 4.5. Sin embargo, se prefiere no eliminar estos elementos antes del anlisis ya
que podran presentar un comportamiento favorable para el objetivo de ste estudio.

Tabla 4.5. Elementos con ndice de correlacin superior o igual a 0.93.

Elementos ndice de correlacin


Ce-La 0.98
Ga-Mg 0.94
Li-Mg 0.93
V-Sc 0.94

52
Finalmente, la base de datos geoqumicos constara de un total de 8240 muestras
en el espacio con informacin de concentracin para 42 elementos qumicos, 10
mayores: Al, Ca, Na, K, Mg, Fe, S, Cu, P y Mn; 32 elementos traza y su informacin de
mapeo: litologa, alteracin hidrotermal, minzone y abundancia relativa de minerales.

La mineraloga que se considerar para el estudio de las alteraciones hidrotermales


ser detallada ms adelante en la seccin de anlisis mineralgico.

4.5 Anlisis de componentes principales

El anlisis de componentes principales (ACP) busca reducir el nmero de variables


determinando vectores propios cuyos valores propios describan la mxima varianza en
todas las variables (ver Anexo A: Componentes principales). El ACP permite describir y
analizar las relaciones entre varias variables, jerarquizando la informacin en factores
ortogonales. Este anlisis se realizar para los elementos mayores y trazas por
separado.

4.5.1 Elementos mayores

En la Figura 4.16, se muestra que las 4 primeras componentes principales explican


un 79.7% de la varianza, mientras que la Figura 4.17 muestra que solo las primeras 3
poseen valores propios asociados mayores a 1.

53
Porcentaje de varianza acumulada
100.0 99.0 99.6100.0
94.8 97.3
90.0 91.6
87.7
80.0 79.7
70.0 71.1
60.0 61.5
50.0
44.8
40.0
30.0
20.0
10.0
0.0
PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 PC9 PC10 PC11

Figura 4.16. Varianza acumulada segn ACP de elementos mayores.

Diagrama de sedimentacin
6.00

5.00 4.93

4.00

3.00

2.00
1.83

1.00 1.07 0.94


0.89
0.42 0.36 0.27
0.00 0.19 0.07 0.04
PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 PC9 PC10 PC11

Figura 4.17. Diagrama de sedimentacin para elementos mayores.

Los coeficientes de correlacin entre los elementos mayores indican que existe una
fuerte correlacin entre dos grupos de elementos (ver Anexo A). Por un lado se
encuentran Al, Mg, K, Ti, Na, P y Mn, y por otro, Fe, S y Ca. El primer y segundo grupo
se correlaciona fuertemente con las componentes principales 1 y 2, respectivamente.
Los crculos de correlacin muestran una proyeccin de la esfera de las variables en
dos componentes, representadas por un punto dentro de un crculo unitario. Si la

54
distancia entre la proyeccin de las variables es baja entonces su correlacin es alta,
sin embargo, no siempre la proximidad en las proyecciones implica una proximidad en
la esfera, a menos que las proyecciones se encuentren cerca de la circunferencia del
crculo.

En la Figura 4.18, se resalta la relacin entre los elementos mayores. Se destaca el


grupo de elementos conformado por Al, Mg, K, Ti, Na, Mn y P, los cuales muestran una
alta correlacin entre ellos asocindose con el primer factor. Adems, K y Ti muestran
un comportamiento conjunto pudiendo alejarse de las dems proyecciones. El grupo de
S, Fe y Ca tienen una correlacin baja con el grupo anterior, y se asocia con el segundo
factor. Por su parte, el Cu parece tener un comportamiento inverso con el resto de los
elementos, y no se explica completamente por las primeras 3 componentes principales.

Figura 4.18. Crculos de correlacin para elementos


mayores.

55
El diagrama de dispersin (ver Figura 4.19) entre el primer y segundo factor resalta
que la primera componente es capaz de clasificar los datos en 3 grupos. El primer
grupo presenta los valores ms altos dentro del PC1 y corresponde principalmente a la
alteracin K1, mientras que el segundo grupo corresponde a los valores de PC1
menores a cero y se caracteriza por las alteraciones flicas S1 y S2. El tercer grupo se
encuentran las alteraciones SCC1, SCC2 y K2, y corresponderan a una zona
transicional entre los dos grupos mencionados anteriormente.

Figura 4.19. Diagrama de dispersin entre PC1 y PC2. Los colores indican
alteraciones hidrotermales.

El comportamiento conjunto de los elementos mayores puede reflejar directamente


la mineraloga de la cual provienen. El grupo de elementos asociados con la PC1,
muestran un comportamiento capaz de separar las alteraciones neutras de las
hidrolticas. El Al, Mg, K, Ti, Mn y P se encuentran en mayor proporcin en las
alteraciones potsicas y sericita-clorita, que en las alteraciones flicas. Esto puede
reflejar el control mineralgico que ejerceran los filosilicatos de biotita, clorita y xidos
de hierro como magnetita. Por otro lado, el grupo de Fe, S y Ca no muestran una
clasificacin asociada a la PC2. Esto puede deberse principalmente a que los
56
elementos no muestran una correlacin total con esta componente, lo que se explicara
por las mltiples proveniencias de estos elementos, como sulfuros, sulfatos, carbonatos,
xidos e hidrxidos de hierro, entre otros. De esta manera, el comportamiento de los
elementos mayores se relacionara con la composicin de aquellos minerales propios
de algunas alteraciones hidrotermales debido a los efectos de la digestin por agua
regia.

4.5.2 Elementos traza

La mayora de los elementos traza muestran un comportamiento que no se


correlaciona con las caractersticas de logueo de sondajes como la litologa, alteracin
o minzone, y parecen distribuirse de forma aleatoria en el yacimiento. Sin embargo,
existen algunos elementos (ver Anlisis univariable) como Zn, Be, Ga, Sc, V e Y que
presentan una particin preferencial segn la alteracin presente en la roca. Se
eliminaron del anlisis aquellos elementos que presentan un ndice de correlacin
mayor a 0.94% (ver Anexo ), de esta manera solo se consideran 29 elementos.

El comportamiento errtico de la mayora de los elementos traza se ve reflejado en


la cantidad de factores necesarios para explicar la varianza. En la Figura 4.20 se
muestra que el porcentaje de varianza acumulada por las primeras 11 componentes
solo alcanza el 72.2%. Adems, en la Figura 4.21 muestra que los valores propios son
inferiores a la unidad a partir del dcimo factor.

57
Porcentaje de varianza acumulada
80.0
70.0 72.2
69.3
66.1
60.0 62.5
58.8
54.8
50.0 50.7
45.5
40.0 39.6
30.0 32.8

20.0 19.9
10.0
0.0
PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 PC9 PC10 PC11

Figura 4.20. Varianza acumulada segn ACP de elementos traza.

Diagrama de sedimentacin
7.0

6.0
5.8
5.0

4.0
3.8
3.0

2.0 2.0
1.7 1.5
1.0 1.2 1.2 1.1 1.1 0.9 0.9
0.0
PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 PC9 PC10 PC11

Figura 4.21. Diagrama de sedimentacin para elementos traza.

A pesar de la gran cantidad de variables, los crculos de correlacin muestran que


existen elementos que se relacionan fuertemente entre s y con los primeros factores.
En la Figura 4.22, se observa que las proyecciones de algunos elementos dentro de los
crculos de correlacin se relacionan fuertemente entre ellos y con la PC1, como: Ce,
Be, Ce, Th, Y, Cs, Ga, Rb y V. El Zn, que es el elemento que posee una mayor
concentracin relativa a los dems, se correlaciona levemente con la PC1, sin embargo,
tambin lo hace con otras componentes lo que habla de su gran dispersin dentro del
yacimiento.
58
Figura 4.22. Diagramas de componentes principales para
elementos traza.

59
Al igual que el primer factor del ACP para elementos mayores, el primer factor de
los elementos clasifica las alteraciones en tres grupos: potsica biottica (K1), flicas (S1
y S2) y transicionales (SCC1, SCC2 y K2) (ver Figura 4.23).

Figura 4.23. Diagrama de dispersin entre PC1 y PC2. Los colores indican
alteraciones hidrotermales.

4.6 Clasificacin Previa

En estudios anteriores (Townley et al., 2012) esta base de datos fue clasificada
utilizando parmetros construidos a partir de la representacin grfica del
comportamiento de los elementos qumicos. Esta clasificacin reagrupa las alteraciones
hidrotermales segn una mineraloga terica. A continuacin se describen los 3 pasos
de clasificacin.

60
1. Diagrama de clasificacin Al-Fe-Mg: La lnea caracterizada por la razn
Al/(Mg+Al)=80% divide el tringulo en dos zonas, denominadas zona 1 (>biotita,
>clorita, <sericita) y zona 2 (>sericita, arcillas), ambos con pirita. Esta etapa es
capaz de clasificar la base de datos geoqumicos con una eficiencia del 92,4%
(ver Figura 4.24).

Figura 4.24. Diagrama de clasificacin Al-Fe-Mg (Townley et al., 2012).

2. Diagrama de clasificacin Ca-Fe-S: Se divide el tringulo en dos zonas segn los


siguientes criterios:
Ca
10%
(Ca Fe S)
Fe
50%
Ca Fe S

61
De esta manera, es posible distinguir aquellas muestras que tericamente
contendran anhidrita de las que no (ver Figura 4.25).

Figura 4.25. Diagrama de clasificacin Ca-Fe-S (Townley et al., 2012).

3. Diagrama de clasificacin (Al+K)/(Na+Ca+Mg) vs Sc*V: Este diagrama destaca 4


zonas bien definidas (ver Figura 4.26): zona de alteracin sericita-cuarzo-arcillas,
clorita-sericita-potsica, potsica y con anhidrita/yeso. El criterio propuesto es
que aquellas muestras cuyo producto Sc*V sea mayor a 250 seran clasificados
como alteraciones potsicas.

62
Figura 4.26. Diagrama de discriminacin entre alteraciones usando ndice
(Al+K)/(Na+Ca+Mg) vs Sc*V. Arriba se indica con una lnea roja el lmite para la
clasificacin como alteracin potsica. Abajo se indica la densidad de datos
(Townley et al., 2012).
63
De acuerdo a la clasificacin propuesta por Townley et al., 2012, las muestras de la
base de datos geoqumicos de Escondida fueron catalogadas en 5 clases distintas.
Como se observa en la Tabla 4.6, la alteracin potsica biottica fue clasificada en su
mayora como K1 y SCC, mientras que la alteracin potsica feldesptica fue clasficada
como SCC y SCC+Anh, pudiendo incluso caer en el campo de las alteraciones flicas y
arglica avanzada. Las alteraciones SCC fueron clasificadas preferentemente como
SCC, K1 y SCC+Anh. Las alteraciones S1, S2 y AA fueron clasificadas preferentemente
como S+AA pudiendo caer en el campo de las SCC.

Tabla 4.6. Clasificacin entre clases propuestas por Townley et al. (2012).

% de clasificacin por geoqumica


K1 SCC SCC+Anh S+AA S+Anh Cantidad total
K1 83% 15% 2% 0% 0% 264
Alteracin Mapeada

K2 2% 71% 20% 4% 2% 122


SCC1 21% 77% 1% 1% 0% 1108
SCC2 15% 68% 15% 2% 0% 1417
S1 1% 11% 1% 86% 1% 2551
S2 0% 9% 1% 90% 0% 1706
AA 0% 7% 0% 91% 2% 57

A medida que la alteracin a clasificar presenta ms sobreimposicin de


alteraciones su dispersin en relacin a los datos geoqumicos aumenta. Las
alteraciones potsicas muestran una mejor clasificacin y un menor rango de variacin
a diferencia de las alteraciones flicas y arglica avanzada que pueden ser clasificadas
frecuentemente como alteraciones SCC. Esto habla de que existe una transicin del
grado de alteracin desde un ncleo potsico a un halo flico pasando por alteracin
clorita-sericita-arcillas. En la Figura 4.27, se evidencia que las alteraciones potsicas

64
caen cercanas a la composicin de la clorita-biotita, mientras que las alteraciones
clasificadas como S+AA y S+Anh se encuentran cerca de la composicin de la pirita,
mineral que predomina en esta alteracin. Por su parte, las alteraciones SCC tienden a
graficar entre estos dos extremos lo que coincide con una etapa de transicin entre
ambas.

Figura 4.27. Diagrama Mg-Fe-S mostrando clasificacin de base de datos


(Townley et al., 2012).

65
5 Mineraloga avanzada

5.1 Base de datos QEMSCAN

QEMSCAN (Quantitative Evaluation Minerals by Scanning electron microscopy) es


una tcnica analtica configurada para medir variaciones mineralgicas basado en la
geoqumica a escala microscpica. Esta tcnica realiza mediciones de dos seales:
backscattered electron (BSE) y Energy Dispersive Spectra (EDS). La primera es usada
para localizar partculas y la segunda para identificar minerales. Posteriormente, la
seal EDS de cada anlisis es clasificada mediante la comparacin con una base de
datos de minerales conocida (Pirrie & Rollinson, 2011).

Los datos de mineraloga cuantitativa Qemscan utilizados en este estudio


corresponden a proporciones de 47 minerales de inters definidos para el sistema tipo
prfido cuprfero de Escondida. Cada muestra cuenta con informacin para distintos
rangos de tamao de partcula entre: 1000-150 m, 150-53 m, 53-5 m, 5-0 m y la
combinacin de todos los anteriores.

5.1.1 Depuracin de base de datos Qemscan

Al igual que en la base de datos geoqumicos, la base de datos de mineraloga


avanzada de Qemscan se eliminaron un total de 30 minerales con rango intercuartil
cercano a cero (<0.1) y con ms del 5% de sus muestras nulas. De esta manera, la
base de datos mineralgica consta de 75 muestras con informacin de concentracin
para 16 minerales pertenecientes a 6 grupos (ver Tabla 5.1), cruzados con la base de
datos geoqumicos e informacin de mapeo: litologa, alteracin, zona mineral e
intensidad de presencia de minerales.

66
Tabla 5.1. Minerales por grupo mineral identificados en la base de datos Qemscan

Grupo Minerales
Tectosilicatos cuarzo, ortoclasa, albita y plagioclasa
Filosilicatos biotita, clorita, mica blanca/sericita, caolinita y pirofilita/esmectita
Sulfuros pirita, calcopirita, calcosina
Hidroxisulfato alunita
xidos xidos de Fe
Fosfatos apatito y svanbergita

La distribucin de la cantidad de muestras por alteracin y litologa no es


homognea (ver Tabla 5.2). Las alteraciones K2 y AA solo presentan una muestra cada
una y se encuentran en tramos con superposicin de alteraciones. Como es de esperar,
la alteracin potsica biottica (K1) se encuentra solo en la andesita. Las alteraciones
SCC y S se encuentran preferentemente en la andesita y prfido feldesptico,
respectivamente. A excepcin de las SCC, estas alteraciones presentan un nmero de
muestras mayor por litologa.

Tabla 5.2. Cantidad de muestras de mineraloga QEMSCAN por litologa y


alteracin.

Alteracin Prfido Feldesptico Andesita


K1 0 4
K2 1 0
SCC1 1 11
SCC2 5 11
S1 13 6
S2 13 10
AA 1 0
Total 35 42

67
5.1.2 Mineraloga Soluble

Dado que los tectosilicatos no son solubles en agua regia, su presencia no se ver
reflejada en la geoqumica, y por lo tanto, no se considerarn para el anlisis
mineralgico. De esta manera, la abundancia de minerales por alteracin hidrotermal
muestra un predominio de los filosilicatos por sobre los dems minerales (ver Figura
5.1).

Figura 5.1. Concentraciones promedio de minerales Qemscan por alteracin


hidrotermal.

En la Figura 5.2, se observan las proporciones minerales promedio para 6


alteraciones hidrotermales. Los minerales ms abundantes son la sericita, biotita,
clorita, pirofilita/esmectita y caolinita, seguido por xidos de hierro, sulfuros, fosfatos y
alunita.

68
100%
90%
Svanbergite
80% Alunite
Porcentaje relativo

70% Chalcocite

60% Apatite
Chalcopyrite
50%
Pyrite
40% Fe Oxides
30% Kaolinite
20% Pyrophyllite/Smectite
Biotite
10%
Chlorite
0%
K1 SCC1 SCC2 S1 S2 AA Muscovite/Sericite
Alteraciones hidrotermales

Figura 5.2. Proporciones relativas de minerales solubles por alteracin


hidrotermal.

Segn la clasificacin de las alteraciones hidrotermales en funcin de su pH,


temperatura y profundidad (Corbett & Leach, 1997), es posible distinguir dos familias de
alteraciones: neutras y cidas. En la primera se encuentran las alteraciones potsica
(K1 y K2) y clorita-sericita (SCC1 y SCC2), mientras que en la segunda familia se
encuentra las alteraciones flicas (S1 y S2) y arglica suprgena (AA).

La sericita y pirita estn presentes en todas las alteraciones, pero en mayor


proporcin en las alteraciones cidas. Por otro lado, la biotita, clorita y xidos de hierro
se encuentran preferentemente en las alteraciones neutras. La alunita es un mineral
que se solo se encuentra en concentraciones importantes en la alteracin arglica
avanzada, y por lo tanto, es un mineral ndice. Los fosfatos presentan un
comportamiento opuesto, el apatito y la svanbergita se particionan preferentemente en
las alteraciones neutras y cidas, respectivamente. El grupo de las arcillas, conformado
por pirofilita/esmectita y caolinita, y los sulfuros de cobre, calcopirita y calcosina, no
presentan una distribucin preferencial entre las alteraciones (ver Figura 5.3).

69
Figura 5.3. Histogramas de distribucin de mineraloga cuantitativa Qemscan por
alteracin hidrotermal.

Los minerales que mejor se particionan en funcin de las alteraciones son biotita,
clorita, sericita, xidos de Fe, pirita, alunita, apatito y svanbergita. Tomando en cuenta
su distribucin es posible dividir las alteraciones en neutras y cidas. Las alteraciones
neutras se caracterizan por presentar mayor contenido de biotita, clorita, xidos de
hierro y apatito, mientras que las cidas poseen mayor proporcin de pirita, sericita,
alunita y svanbergita. Es decir, estos minerales presentan un comportamiento opuesto.
Los dems minerales no muestran este comportamiento.

Los diagramas de clasificacin de alteraciones hidrotermales propuestos por


Townley et al. (2012) para las muestras de mineraloga avanzada evidencian que la
separacin de los datos se relaciona directamente con la mineraloga presente en la
roca. Las zonas dentro de los diagramas Al-Mg-Fe separan los datos en dos grupos,
dominados por minerales de pirita y biotita-clorita. El diagrama de enriquecimiento en
anhidrita/yeso muestra que no existen datos presencia significativa de estos minerales,
lo que concuerda con la mineraloga Qemscan (ver Figura 5.4). Finalmente, el diagrama
(Al+K)/(Na+Ca+Mg) vs Sc*V muestra que efectivamente los datos con mayor
proporcin de biotita-clorita se encuentran en la zona definida para esta mineraloga de
alteracin (ver Figura 5.5).
70
Figura 5.4. Diagramas de clasificacin segn Townley et al., 2012. El tamao de
los crculos indica abundancia relativa de pirita a la izquierda y clorita a la
derecha.

71
Figura 5.5. Diagrama de clasificacin segn Townley et al. (2012) para alteracin
K1. El tamao de los crculos indica abundancia relativa de biotita.

5.1.3 Mineraloga y geoqumica

Para analizar la mineraloga y la geoqumica se estudiar el comportamiento


conjunto de los elementos mayores que componen los minerales solubles en agua
regia. De esta manera, se determinar qu minerales son aquellos cuya presencia se
refleja en la base de geoqumica. Adems, se pretenden estudiar el comportamiento de
los elementos para determinar cules son las fuentes minerales ms probables de cada
uno.

En la Tabla 5.3, se muestra la estequiometra de los 14 minerales considerados


solubles. Como se observa, los elementos pueden encontrarse en ms de una especie
mineral, como el caso del aluminio y hierro.

72
Tabla 5.3. Frmula estructural de los minerales analizados por QEMSCAN
considerados como solubles en agua regia.

Mineral Frmula
calcosina Cu2S - Cu1.4S
calcopirita CuFeS2
pirita FeS2
xidos de Fe FeO(OH) / Fe2O3 / Fe3O4
muscovita/sericita KAl3Si3O10(OH,F)2
caolinita Al2Si2O5(OH)4
pirofilita/esmectita Al2Si4O10(OH)2 / A0.3D2-3[T4O10]Z2.nH2O
biotita K(Mg,Fe2+)3AlSi3O10(OH,F)2
clorita (Mg,Fe2+)5Al2Si3O10(OH)8
alunita (Na,K)Al3(SO4)2(OH)6
svanbergita (Sr,Pb)Al3(PO4)(SO4)(OH)6
apatito Ca5(PO4)3(OH,F,Cl)
yeso/anhidrita CaSO4.2H2O / CaSO4
carbonatos CaCO3 / (Mg,Ca)CO3

Las razones de concentraciones entre minerales de alteracin Qemscan muestran


que existe una particin de estos en funcin de la alteracin hidrotermal predominante
en la roca. Las razones clorita/sericita y biotita/sericita son mayores en las alteraciones
potsicas biotticas (K1) y clorita-sericita (SCC1 y SCC2), a pesar de que la sericita sea
un mineral predominante en la base de datos. Lo mismo ocurre con las razones de
estas micas con pirita (ver Figura 5.6).

73
Figura 5.6. Histogramas de razones de minerales Qemscan.

Es posible identificar una correlacin entre la mineraloga y geoqumica mediante la


visualizacin de grficos de concentraciones molares de elementos. En ellos se utilizan
lneas de control mineral, calculadas a partir de sus frmulas tericas, y se vara el
tamao de los crculos que identifican cada muestra para indicar la razn relativa de
minerales Qemscan: (clorita+biotita)/sericita y (clorita+biotita)/pirita. De esta manera, las
muestras de alteracin potsica biottica y clorita-sericita, con mayor concentracin
clorita-biotita relativa a sericita, tienden a graficar entorno a las lneas de tendencia de
estos minerales. Lo contrario ocurre para las alteraciones flicas y arglica avanzada
donde sta relacin se invierte (ver Error! No se encuentra el origen de la
referencia. y Figura 5.8).

74
Figura 5.7. Grficos molares de dispersin de Al-Mg, Al-K, Mg-K. El tamao de los
crculos indica la razn clorita/sericita.

75
Figura 5.8. Grficos molares de dispersin de Al-Mg, Al-K, Fe-Mg, Fe-S, Mg-K y
Mg-S. El tamao de los crculos indica la razn clorita/pirita.

A partir de esta correlacin y del hecho de que algunos grupos minerales son
altamente solubles en agua regia, es posible identificar la proveniencia ms probable de
los elementos mayores en funcin de la estequiometra terica de cada mineral. A
continuacin se proponen los minerales fuente de los elementos mayores usando la
estequiometra de la Tabla 5.3.

El aluminio proviene preferentemente de las micas que definen las alteraciones:


biotita, clorita, arcillas y, probablemente, sericita.

76
El magnesio proviene de biotita, clorita y en menor proporcin de carbonatos
(dolomita). Estos minerales son caractersticos de las alteraciones potsicas y SCC, por
lo que se correlacionan con altas concentraciones de magnesio.

El potasio tiene como principal fuente terica la biotita, sericita y alunita. Sin
embargo, la bibliografa y las lneas de control mineral indican que este elemento
proviene principalmente de la biotita y no de la sericita o alunita. Es solo debido a la
ausencia de biotita que las muestras donde predomina la sericita grafiquen en su
respectiva lnea.

El azufre proviene preferentemente de pirita (sulfuro predominante), calcopirita y


calcosina-covelina, en menor proporcin. Adems, puede provenir de algunos sulfatos
como la anhidrita/yeso y la alunita.

El hierro posee varias fuentes, ya que puede provenir de pirita, calcopirita, xidos
de hierro (magnetita-hematita-oxihidrxidos), biotita y clorita. A excepcin de la pirita,
todos los minerales mencionados se encuentran preferentemente en las alteraciones
potsicas y sericita-clorita.

El calcio es un elemento que se encuentra en bajas concentraciones y se


correlaciona directamente con la presencia de anhidrita/yeso y calcita. Ambos minerales
son altamente solubles en agua regia.

6 Anlisis preliminar de variables y definicin de grficos de inters

El anlisis previo indica que existen minerales cuya qumica se ve reflejada en el


comportamiento conjunto de los elementos. Algunos de estos minerales presentan una
particin preferencial segn la alteracin hidrotermal a la que se asocien, como el caso
de la biotita y clorita en las alteraciones K y SCC. Del mismo modo, existen elementos
que presentan una distribucin relativamente homognea entre alteraciones y no son
de inters para la clasificacin.

77
En esta seccin, se seleccionan los elementos que presenten una mayor
diferenciacin entre alteraciones hidrotermales, para luego utilizarlos en la construccin
de la clasificacin.

6.1 Seleccin de elementos

Una forma de seleccionar los elementos es mediante el clculo de traslape entre


alteraciones para cada uno de ellos, como se explica a continuacin.

Por cada elemento, se consideran el 75% de los datos con solo una alteracin
hidrotermal por tramo, lo cual define un rango de valores comprendidos entre los
percentiles 12,5 y 87,5 (ver Tabla 6.1).

Tabla 6.1. Percentiles definidos para el 75% de los datos por alteracin para
magnesio.

Mg wt%
Alteracin Cdigo percentil 12,5 percentl 87,5
30 K1 1,20 1,89
31 K2 0,19 0,50
40 SCC1 0,80 1,69
41 SCC2 0,30 1,62
51 S1 0,02 0,16
52 S2 0,02 0,09
61 AA 0,02 0,12

Para cada rango de valores por alteracin, se calcula el porcentaje de datos


contenidos en las dems categoras. En la Figura 6.1, se esquematizan los 3 casos
posibles de traslapes. En A, el traslape es total independiente de que exista parte del
rango de la alteracin 30 que no est contenido en 31. En B el traslape es nulo, lo que
indica que este elemento separa efectivamente las alteraciones. En C, el traslape es
parcial (zona roja) y se calcula el porcentaje de traslape de 30 en 31 y viceversa, para
luego obtener un promedio entre ambas poblaciones.

78
Figura 6.1. Casos posibles de traslapes entre las alteraciones 30 y 31.

De esta manera, se cuantifica cuan separables son las 21 combinaciones de pares


de alteraciones hidrotermales por cada elemento qumico. En el Tabla 6.2, se muestran
las combinaciones de alteraciones que poseen un porcentaje de traslape mayor al 15%.
En particular, las alteraciones con mayor traslape mnimo corresponderan a SCC1-
SCC2 (40-41) y S1-S2 (51-52).

Tabla 6.2. Menor porcentaje de traslape obtenido para las combinaciones de


alteraciones indicadas.

Alteraciones Cdigos Elemento Mnimo traslape


30-40 K1-SCC1 V (Sc, Rb, Ca) 34%
30-41 K1-SCC2 V (Sc, Na, Tl, K) 19%
40-41 SCC1-SCC2 Co (Al, Na) 67%
51-52 S1-S2 Mn 68%
51-61 S1-AA Mn 55%
52-61 S2-AA As 50%

Los elementos con un menor porcentaje mnimo de traslape son los que presentan
una mayor diferenciacin entre alteraciones. Los elementos mayores que permiten
diferenciar entre alteraciones son Mg, Al, Fe, Ca, K, Na, P y Mn. Los elementos traza
son V, Sc, Rb, Ga, Y, Be y Co. Esta cantidad puede reducirse an ms ya que existen
elementos con una alta correlacin por lo que simplemente introducen informacin
redundante como es el caso de Mg-Ga y V-Sc (ver Tabla 6.3).

79
Tabla 6.3. Pares de alteraciones con porcentaje de traslape menor a un 15%.

Mayores Traza mnimo traslape


30-31 Mg, Al, Fe, P V, Sc, Rb, Ga, Y, Co 0%
30-51 Mg, Al, Ca, K, Na, P, Mn V, Sc, Rb, Ga, Y, Co 0%
30-52 Mg, Al, Ca, K, Na, P, Mn V, Sc, Rb, Ga, Y, Be 0%
30-61 Mg, Al, Ca, Na, P, Mn V, Sc, Rb, Ga 0%
31-40 Mg, Al, Fe Sc, Co 0%
31-41 Al Co 14%
31-51 Mg, Mn Ga 0%
31-52 Mg, Mn Ga 0%
31-61 Mg, Mn - 0%
40-51 Mg, Al, P, Mn V, Sc, Ga, Y, Be 0%
40-52 Mg, Al, P, Mn V, Sc, Ga, Y, Be 0%
40-61 Mg, Mn V, Ga 0%
41-51 Mg, Mn Ga 0%
41-52 Mg, Mn Ga 0%
41-61 Mg, Mn - 0%

El clculo de traslapes entre poblaciones permite concluir que existen alteraciones


hidrotermales que no pueden ser diferenciadas entre s: SCC1-SCC2 y S1-S2-AA. Del
mismo modo, es posible diferenciar las alteraciones K-SCC de las S-AA.

Al igual que para las alteraciones hidrotermales, se buscan los elementos que mejor
clasifiquen las zonas de mineralizacin utilizando el clculo del traslape entre las
poblaciones. Siguiendo la misma metodologa, los elementos mayores que presentan
un menor traslape entre zonas de mineralizacin son Cu, Fe, S y Al. Los elementos
traza son Pb, Co, Sc, V y Cd.

A pesar de que estos elementos presentan un bajo traslape, la mayora de ellos no


son capaces de separar grupos de clases. Solo en el caso del cobre es posible definir
una diferencia significativa entre las muestras que presentan mineralizacin y las que
no (ver Figura 6.2).

80
Figura 6.2. Boxplot de minzones para Cu total.

6.2 Alteraciones y mineralizacin por sondaje

Para analizar el comportamiento de las alteraciones y minzones en funcin de la


qumica dentro de un sondaje, se eligen sondajes que presenten dos o ms
alteraciones mapeadas y se grafican las concentraciones de elementos mayores a lo
largo de estos.

Las alteraciones hidrotermales potsicas biotita y sericita-clorita se relacionan a


altas concentraciones de Al, Mg y K, con respecto a las alteraciones flicas. Es
importante destacar que el contacto entre zonas de alteracin puede presentarse en la
qumica como variaciones progresivas de las concentraciones o como variaciones
puntuales (ver Figura 6.3).

81
Figura 6.3. Concentraciones de elementos a lo largo de sondajes resaltando el
contacto entre alteraciones.

82
De igual forma, existen sondajes con variaciones qumicas no asociadas a
cambios en la alteracin, litologa o mineralizacin, o que presentan un comportamiento
contrario al descrito anteriormente (ver Figura 6.4), lo que puede deberse a variaciones
de mineraloga no identificadas durante el mapeo de los sondajes. Las muestras que
presentan esta incongruencia seran de especial inters y deberan ser revisadas con
mayor detalle en futuras campaas de mapeo, con el objetivo de verificar si la
geoqumica se correlaciona efectivamente con la presencia de minerales ndices de las
alteraciones o son variaciones propias de cada una de ellas.

83
Figura 6.4. Concentraciones de elementos a lo largo de sondajes resaltando
variaciones no relacionadas al cambio de alteracin.

84
Figura 6.5. Concentraciones de elementos a lo largo de sondajes resaltando
variaciones de potasio en una misma alteracin.

85
Las concentraciones de potasio, en general, son ms altas en las muestras con
alteracin K1, y muy bajas en las muestras con alteraciones tipo S. Las alteraciones
SCC pueden presentar cualquiera de los dos patrones (ver Figura 6.5). La Figura 6.6
indica que las muestras con mayor razn de minerales Qemscan clorita:sericita poseen
una mayor concentracin de K, asociada al contenido de biotita representado por el
tamao de los puntos. Este comportamiento indicara que la principal fuente de potasio
es la biotita.

Figura 6.6. Razn de minerales Qemscan clorita:sericita vs K wt%. Tamao indica


concentracin de biotita Qemscan.

Las variaciones de las alteraciones SCC pueden ser explicadas por la razn
entre clorita y sericita. Una roca con alteracin SCC y una razn chl/ser baja tendera a
presentar bajas concentraciones de K debido al bajo contenido en biotita remanente de
la alteracin potsica. Por otro lado, una roca con alteracin SCC y con una razn
chl/ser alta tendera a presentar altas concentraciones de K debido a que las biotitas del
sistema solo se encontrara parcialmente alteradas a clorita. Esto destaca la
importancia que tienen las alteraciones cercanas a una muestra en particular. Las
86
transiciones entre alteraciones y las relaciones entre los minerales de alteracin pueden
imprimir patrones geoqumicos distintos a los que presentan las alteraciones bien
definidas.

6.3 Grficos de control mineral

Al considerar las variables que permiten distinguir entre alteraciones y zonas de


mineralizacin, es posible identificar rangos de valores que entregan distinta
informacin. Estas zonas sern definidas a partir de los elementos mayores debido a
que se correlacionan directamente con la mineraloga. Esta informacin es
complementaria a la clasificacin previa (Townley et al., 2012).

Los elementos mayores permiten distinguir muestras con un enriquecimiento


relativo asociado a la presencia de minerales ndices. Las variables Al, Mg y K permiten
diferenciar entre zonas de alta y baja presencia de biotita y clorita. Los elementos Fe y
S indican si las muestras se encuentran bajo el control mineral de la pirita, o si
presentan un enriquecimiento/empobrecimiento con respecto a esta. Si el
enriquecimiento es en Fe, la muestra se asocia a alteraciones con mayor contenido en
biotita, clorita u xidos de hierro, y si el enriquecimiento es de S, entonces la muestra
presentara un mayor contenido en sulfatos como anhidrita/yeso. Adems, el CuT
entrega informacin de si la mineralizacin corresponde a solo pirita (minzone PR1) o
presenta mineralizacin econmica asociada a sulfuros de cobre como calcopirita,
bornita, covelina, calcosina u otros.

La Figura 6.7 muestra el sistema Cu-Fe-S en donde se resaltan tres zonas distintas.
La zona central se caracteriza por el control de pirita y sulfuros de cobre. Estos
minerales son disueltos totalmente por el agua regia y por lo tanto se ven reflejados en
la qumica de estos elementos. Las otras dos zonas de enriquecimiento relativo de
hierro y azufre destacan la presencia de minerales como biotita, clorita y xidos de Fe
para el caso del hierro, y sulfatos para el caso del azufre. El porcentaje de datos
asignados en cada zona se muestra en la

Tabla 6.4.

87
Figura 6.7. Diagrama ternario Cu-Fe-S indicando zonas de pirita-sulfuros de Cu y
enriquecimiento de Fe y S.

Tabla 6.4. Porcentaje de muestras por zona en sistema Cu-Fe-S.

Alto Fe Pirita - sulfuros de Cu Alto S


K1 93% 7% 0%
K2 26% 70% 4%
SCC1 88% 12% 0%
SCC2 64% 30% 6%
S1 4% 95% 1%
S2 4% 96% 0%
AA 0% 100% 0%

88
El sistema Mg-Fe-S (ver Figura 6.8) muestra una particin de los datos entre dos
nodos definidos por biotita-clorita y pirita. Las alteraciones tipo S y AA son ms
cercanas al nodo de pirita, y las K1 y SCC al de biotita-clorita. Tambin es posible
definir una zona de transicin entre ambos nodos en la cual se encuentra un porcin
importante de las alteraciones SCC y K2 (ver Tabla 6.5).

Figura 6.8. Diagrama ternario Mg-Fe-S indicando zonas de control mineral de


biotita-clorita, pirita y transicin entre ambas.

89
Tabla 6.5. Porcentaje de muestras por zona en sistema Mg-Fe-S.

Control Pirita Transicin Control Biotita - Clorita


K1 1% 8% 91%
K2 28% 54% 15%
SCC1 4% 28% 68%
SCC2 17% 40% 43%
S1 92% 7% 1%
S2 94% 5% 1%
AA 100% 0% 0%

El sistema Mg-K-Ca (ver Figura 6.9) muestra la particin entre los nodos de biotita-
clorita y sericita, pudindose definir una zona enriquecida en calcio asociada a la
presencia de anhidrita/yeso o calcita. De esta manera, se definen la zona de control
Biotita-Clorita y la zona de control de Sericita. A pesar de que la biotita es la principal
fuente de potasio, el diagrama muestra que el enriquecimiento relativo de este elemento
sera producto de la presencia de sericita. Sin embargo, este comportamiento se explica
por un empobrecimiento de Mg en las alteraciones flicas y no a que la sericita sea la
fuente de K. La particin de los datos en estos grupos se muestra en la Tabla 6.6.

90
Figura 6.9. Diagrama ternario Mg-K-Ca indicando zonas de Control biotita-clorita y
Control de sericita.

Tabla 6.6. Porcentaje de muestras por zona en sistema Mg-K-Ca.

Control Sericita Control Biotita - Clorita


K1 0% 100%
K2 20% 80%
SCC1 1% 99%
SCC2 7% 93%
S1 90% 10%
S2 92% 8%
AA 95% 5%

91
Estos diagramas entregan informacin preliminar til para la clasificacin acerca de
los controles ejercidos por la mineraloga soluble en agua regia. Como se aprecia, estos
diagramas permiten distinguir dentro de un mximo de solo tres agrupamientos
principales de los datos.

7 Modelo de clasificacin de alteraciones hidrotermales

Una vez identificados los elementos y parmetros que mejor diferencian entre
alteraciones hidrotermales, se propone construir un modelo de clasificacin
supervisada, que permita cuantificar las alteraciones predominantes en una muestra. El
principal supuesto es que las muestras con una misma alteracin presentan un patrn
geoqumico similar.

Usando criterios estadsticos y geolgicos, se eligen muestras que mejor definan


las alteraciones hidrotermales. Se definen clusters para cada alteracin caracterizada
por su media y varianza experimental.

Luego, se calcula la distancia a cada uno de los clusters con el objetivo de


cuantificar el grado de pertenencia a cada uno de ellos. Por ejemplo, una muestra
puede tener un porcentaje calculado de 50% K1, 40% SCC1 y 10% S1.

Finalmente, este resultado ser comparado con los porcentajes de alteraciones por
tramo calculadas al cruzar la base de datos geoqumicos con el mapeo geolgico.

Este algoritmo de clasificacin est basado en el algoritmo fuzzy c-means, el cual


calcula grados de pertenencia a una cantidad c de clusters para cada una de las
muestras dentro del espacio definido por las variables consideradas. A diferencia de
fuzzy c-means, los clusters son predefinidos y no aleatorios, lo que introduce las
caractersticas geoqumicas del sistema a la clasificacin (Rojas et al., 2008; Liu et al.,
2009; Pimentel & M.C.R de Souza, 2013; Nazari et al., 2013; Xing & Ha, 2014).

92
7.1 Definicin de muestras por alteracin hidrotermal

Los criterios que deben cumplir las muestras para definir los clusters por alteracin
hidrotermal son los siguientes:

1. Solo una alteracin mapeada por tramo.

2. Mapeo de abundancia relativa de minerales por tramo, que ser utilizado como
criterio para validar el mapeo geolgico.

3. Espacialmente, deben encontrarse cercanas unas de otras, y a su vez, lo ms


lejanas posibles de los lmites con otras alteraciones.

En la seccin 4.2.2, se muestran las abundancias relativas de minerales,


destruccin de feldespatos y alteracin de mficos. Con esta informacin, es posible
definir muestras cuyo mapeo de minerales sea ms representativo de un tipo de
alteracin hidrotermal. Por ejemplo, la alteracin potsica biottica (K1) debe poseer un
alto ndice de biotita (2 o 3) y uno bajo de sericita (0 o 1), ya que de lo contrario, es muy
probable que exista sobreimposicin de alteraciones.

En la Tabla 7.1 se muestran los ndices de abundancia relativa que se considerarn


para definir una alteracin hidrotermal tipo. Estos ndices se proponen a partir de los
minerales que definen cada alteracin. Aquellos minerales que pueden tomar cualquier
ndice de abundancia relativa se indican con un guin.

Tabla 7.1. ndices de abundancia relativa de minerales para la definicin de


clusters.

Alteracin
desfel silice clorita biotita feld K sericita altmaf arcilla
hidrotermal
K1 - - 0-1 2-3 0 0-1 1-2 0
K2 - - - - 1-3 0-1 - 0
SCC1 3-4 0-1 1-3 - 0 1-3 - 1-3
SCC2 3-4 1-2 1-3 - 0 1-3 1-2 0
S1 3-4 1-3 0 0 0 2-3 - 0
S2 3-4 1-3 0 0 0 2-3 - 1-3
AA 4 - - - - 2-3 - 2-3

93
Al cruzar la base de datos Qemscan (ver Figura 7.1) con el mapeo de abundancia
relativa mineral (Tabla 7.2), se observa que existe una correlacin entre ambos para el
caso de los minerales de slice (cuarzo), clorita, biotita y muscovita/sericita. Sin
embargo, para el caso de las arcillas, representado por los minerales de caolinita y
pirofilita/esmectita, no existe dicha correspondencia, ya que las alteraciones SCC1 y S2
son definidas por la presencia de estas y no presentan una concentracin
evidentemente mayor a la de las dems alteraciones.

94
Figura 7.1. Abundancia de minerales Qemscan segn alteracin hidrotermal.

95
Tabla 7.2. ndices de mapeo de abundancia relativa de minerales con mineraloga
Qemscan.

Alteracin cantidad slice clorita biotita sericita arcillas feldK


K1 2 1 1 2 1 0 0
SCC1 10 0 1-3 0-1 1-2 1-2 0
SCC2 7 1-2 1-2 0 1-3 0 0
S1 17 1-3 0 0 2-3 0 0
S2 20 1-2 0 0 2-3 1-3 0

Se divide la base de datos geoqumicos en dos grupos. Cada uno de ellos con
aproximadamente el 50% de los datos por alteracin. La primera mitad se usa para
definir los clusters para cada alteracin y la segunda es utilizada para la clasificacin.
La Tabla 7.3 muestra la cantidad de datos asignados a cada grupo de datos.

Tabla 7.3. Cantidad muestras por base de datos con una alteracin por tramo.

Alteracin Definicin clusters Clasificacin Total


K1 32 53 85
K2 25 21 46
SCC1 133 148 281
SCC2 303 491 794
S1 809 810 1619
S2 565 507 1014
AA 10 10 20

7.2 Seleccin de elementos segn desempeo

Con los datos restantes de la base de datos, se evalua el desempeo de


clasificacin de los elementos qumicos filtrados anteriormente mediante la comparacin
de los traslapes por alteracin hidrotermal.

La metodologa para la evaluacin de los elementos por s solos y agrupados se


presenta a continuacin:

96
1. Definicin de cada alteracin hidrotermal considerando las muestras
seleccionadas segn criterios asociados al mapeo de sondajes.

2. Clculo de distancias de cada muestra a clasificar a cada uno de los centroides


definidos por las medias de las alteraciones en el paso 1.

3. Clculo de grado de pertenencia a cada uno de los centroides en funcin de las


distancias. Estos valores de grado de pertenencia pertenecen al intervalo
comprendido entre 0 y 1, y la suma de todos ellos es igual a la unidad para una
muestra.

4. Clculo de error entre la alteracin calculada y mapeada.

De esta manera, es posible acotar el nmero de variables que se analizarn para


muestras con ms de una alteracin por tramo de anlisis qumico.

Cada cluster se define por la media y varianza experimental de las alteraciones. La


media es el punto al cual se calculan las distancias de las muestras por clasificar. Esta
distancia debe considerar la covarianza existente entre las variables, dado que la
dispersin caracteriza la geometra de cada cluster en el espacio definido por las
variables afectando directamente la clasificacin. La Figura 7.2 muestra dos particiones
para el mismo set de datos, en donde un dato puede ser clasificado dentro de uno u
otro grupo si se consideran distintas covarianza entre las variables (Rezaee, 2010).

97
Figura 7.2. Diferentes particiones con igual distancia entre centroides pero con
distinta orientacin para igual set de datos (Tomado de Rezaee, 2010).

Una manera de considerar este efecto es utilizando la distancia de Mahalanobis


para calcular la distancia entre el centroide, correspondiente a la media, y la muestra
(Falasconi et al., 2010; Liu et al., 2009). Esta distancia se define como:

d2 (xj , ai ) = (xj ai )t 1 (xj ai )

Donde xj corresponde a una muestra, ai al centroide de la alteracin i, y es la


matriz de covarianza de las muestras que definen la alteracin i. Cuando la matriz de
covarianza es igual a la matriz identidad, la distancia de Mahalanobis es igual a la
Euclidiana.

Para el clculo del grado de pertenencia, se utiliza un paso del algoritmo conocido
como Fuzzy C-Means, el cual permite calcular grados de pertenencias a un nmero C
de clusters. El grado de pertenencia para un nmero n de muestras y c de cluster
corresponde a una matriz U de n filas y c columnas, donde la sumatoria de las
columnas para cada muestra es igual a la unidad (Rojas et al., 2008; Pimentel & M.C.R
de Souza, 2013; Nazari et al., 2013; Xing & Ha, 2014).

98
U = [uij ] , uij [0,1] , i = 1,2, , c, j = 1,2, , n
nxc

Donde,

uij = 1
i=1

El grado de pertenencia queda definido como:

1
d2 (xj, ai )
uij = [ ]
d2 (xj, ak )
=1

Para cada centro definido por un grupo de muestras representativas se espera que
las muestras mapeadas como este tipo de alteracin sean clasificadas como tal. Por
ejemplo, se espera que las muestras mapeadas como cuarzo-sericita, cuarzo-sericita-
arcillas y arglicas suprgenas, estn muy distantes al centro definido por las muestras
de alteracin potsica biottica y que por lo tanto no sean clasificadas como tales.

En la etapa de evaluacin de las variables segn la clasificacin de las alteraciones,


se calculan los errores con respecto a un comportamiento esperado. Segn las
distancias obtenidas anteriormente, la menor de ellas definir a qu alteracin se
asignar, es decir, se utilizar el algoritmo del vecino ms cercano, asignndoles el
grado de pertenencia igual a 1 al mayor de ellos y 0 a los dems. En etapas posteriores
a este anlisis, se utilizarn los grados de pertenencias en lugar de los lmites duros.

De esta manera, se calcula el error cuadrtico de la clasificacin comparando las


alteraciones mapeadas y las alteraciones asignadas. Para evitar que las alteraciones
con mayor cantidad de muestras adquieran un mayor peso en el clculo, el error
asociado a cada una de las alteraciones mapeadas es dividido por el nmero de
muestras de la siguiente manera:

99
nk
1 2
MSE(k, i) = (j j )
nk
j

Donde k es la alteracin mapeada, nk la cantidad de muestras pertenecientes a


esta alteracin, i es el cluster definido para la alteracin i-sima, j es una modificacin
del grado de pertenencia j , tal que es 1 para el grado de pertenencia mayor y 0 para
los dems, y j se define como:

1 , alteracin mapeada xj = alteracin definida para cluster i


j = {
0 , ~

Finalmente, se obtiene una matriz de error para cada grupo de elementos


considerados para el anlisis. En cada columna de esta matriz de 7x7 se encuentra el
error asociado al clster definido por la alteracin i-sima, y en cada fila, el error con
respecto al comportamiento esperado. Por ejemplo, la primera columna representa el
clster definido por la alteracin potsica K1, y cada fila corresponde al error asociado a
todas las muestras mapeadas como alteracin K1, K2, SCC1, SCC2, S1, S2 y AA. En
este caso, j es igual a 1 para las muestras mapeadas como K1, y 0 para los dems.
En la Tabla 7.4, se muestra la matriz de error para la clasificacin obtenida usando solo
el elemento magnesio. En la primera columna, el error es nulo para las alteraciones K2,
S2 y AA, es decir, no existen muestras mapeadas como estas alteraciones que sean
asignadas al cluster K1. Las alteraciones mapeadas como K1 asignadas al cluster K1,
es decir, cuyo valor de j es igual a 1, tienen un error igual a 0.21.

100
Tabla 7.4. Matriz de error asociada al Mg. Media matriz= 0.172. Media diagonal =
0.60.

Alteraciones por cluster


K1 K2 SCC1 SCC2 S1 S2 AA
K1 0.21 0.00 0.15 0.06 0.00 0.00 0.00
Alteraciones K2 0.00 0.24 0.00 0.05 0.00 0.00 0.19
mapeadas SCC1 0.41 0.14 0.85 0.29 0.01 0.00 0.01
SCC2 0.20 0.41 0.12 0.81 0.07 0.00 0.01
S1 0.01 0.09 0.00 0.01 0.93 0.04 0.77
S2 0.00 0.08 0.00 0.02 0.07 0.96 0.80
AA 0.00 0.10 0.00 0.00 0.00 0.10 0.20

En la Tabla 7.5, se muestra el error total de clasificacin asociados a los elementos


mayores. La segunda columna de la tabla corresponde al error asociado solo a las
muestras con alteracin mapeada igual a la definida para el cluster, es decir, para el
valor de j =1. Cada columna corresponde a la media de la matriz de error y la media de
la matriz de error, respectivamente. El comportamiento del magnesio (ver Tabla 7.4) se
reduce a estos dos valores de la Tabla 7.5.

Tabla 7.5. Errores de clasificacin para elementos mayores.

Media matriz de Media diagonal


Elemento [wt%]
error matriz de error
Mg 0.172 0.600
Al 0.199 0.695
K 0.213 0.745
Ca 0.213 0.747
Na 0.213 0.747
CuT 0.215 0.751
Fe 0.215 0.752
S 0.217 0.760

101
En las Figura 7.3 y Figura 7.4, se muestran boxplots de los clusters definidos por
las muestras representativas de cada alteracin (K1, K2, SCC1, SCC2, S1, S2 y AA).
En cada uno de los boxplots, se grafican los grados de pertenencia al cluster
correspondiente de las muestras mapeadas por alteracin (30, 31, 40, 41, 51, 52 y 61).

Figura 7.3. Boxplots grados de pertenencia por alteracin para magnesio.

102
Figura 7.4. Boxplots grados de pertenencia por alteracin para azufre.

La principal diferencia entre los resultados obtenidos para Mg y S, es que el Mg


muestra una mejor separacin de las alteraciones si observamos las medias, mientras
que el S tiende a presentar valores similares en la mayora de los clusters. Por ejemplo,
en el caso de la alteracin S1, los datos tienen aproximadamente el mismo grado de
pertenencia en este cluster al considerar solo el azufre, mientras que poseen una clara
distincin al considerar el magnesio (ver Figura 7.3 y Figura 7.4).

Mediante el anlisis del error y la visualizacin del comportamiento de los grados de


pertenencia calculados para cada alteracin, se concluye que los elementos mayores
que mejor reflejan la separabilidad de las alteraciones hidrotermales definidas son Mg,
K, Al, Na y Ca. Los dems elementos mayores no sern considerados en el anlisis
posterior.

Para el caso de los elementos traza, el tratamiento fue idntico. Los errores se
muestran en la Tabla 7.6.

103
Tabla 7.6. Errores de clasificacin para elementos traza.

Media matriz de Media diagonal


Elemento [ppm]
error matriz de error
V 0.178 0.622
Mn 0.185 0.648
P 0.189 0.662
Rb 0.194 0.678
Y 0.198 0.692
Be 0.209 0.733
Co 0.219 0.766

Mediante el anlisis del error y la visualizacin del comportamiento de los grados de


pertenencia calculados para cada alteracin, se concluye que los elementos traza que
mejor reflejan la separabilidad de las alteraciones hidrotermales definidas son: V, Mn, P,
Rb e Y.

Una vez obtenidos los elementos mayores y traza que individualmente reflejan
mejor el comportamiento esperado para obtener una clasificacin ptima en funcin de
las alteraciones definidas previamente, se analiza el comportamiento de todos los
grupos posibles de elementos.

En la Tabla 7.7, se muestran los errores de los grupos de elementos mayores en


orden creciente, y el promedio de los porcentajes mximos de clasificacin por
alteracin. Este ltimo es una medida de qu grupos de elementos tienden a concentrar
ms los datos por cluster. De esta manera, se eligen los grupos de elementos que
posean un error bajo y un alto promedio de clasificacin. Los que cumplen mejor estos
criterios son los grupos de: Al-Mg-K-Ca-Na, Al-Mg-K-Ca y Al-Mg-K-Na.

104
Tabla 7.7. Error promedio y promedio de porcentajes mximos de clasificacin
para grupos de elementos mayores.

Promedio de porcentajes
Grupos Error promedio
mximos de clasificacin
Al - Mg - Ca - Na 0.164 64%
Al - Mg 0.164 50%
Al - Mg - Na 0.165 60%
Al - Mg - Ca 0.167 57%
Al - Mg - K - Na 0.168 67%
Al - Mg - K - Ca - Na 0.168 67%
Mg - K - Na 0.169 62%
Al - Mg - K - Ca 0.170 68%
Mg - K 0.171 58%
Al - K - Ca - Na 0.172 59%
Al - Mg - K 0.172 66%
Mg - K - Ca - Na 0.172 66%
Mg - Na 0.174 56%
Mg - Ca - Na 0.175 60%
Mg - K - Ca 0.175 67%

En general, para el caso de los elementos traza, se alcanza un menor error


promedio pero tambin un menor promedio de porcentajes mximos de clasificacin
que los elementos mayores. Usando los mismos criterios anteriores, los grupos de
elementos trazas elegidos son: P-Rb, P-Y-Rb y P-V-Y-Rb (ver Tabla 7.8).

105
Tabla 7.8. Error promedio y promedio de porcentajes mximos de clasificacin
para grupos de elementos traza.

Promedio de
Grupos Error promedio porcentajes mximos
de clasificacin
P - Rb 0.154 57%
P - Y - Rb 0.156 57%
P - V - Rb 0.161 54%
Mn - V - Rb 0.161 57%
P - V - Y - Rb 0.161 58%
P-V-Y 0.162 51%
Mn - Y - Rb 0.163 63%
Mn - V - Y - Rb 0.164 64%
V - Y - Rb 0.165 62%
Mn - V - Y 0.166 65%
Mn - P - V - Y 0.168 55%
Y - Rb 0.168 59%
Mn - P - Y - Rb 0.169 57%
V - Rb 0.169 49%
P-V 0.170 50%

7.3 Anlisis de los clusters por alteracin hidrotermal

Cada muestra clasificada utilizando los elementos mayores Al-Mg-K-Ca-Na tiene


como resultado un grado de pertenencia a cada uno de los clusters definidos por las
muestras representativas.

En la Figura 7.5 se muestran los promedios de los grados de pertenencia a cada


cluster. Existen alteraciones, como el caso de K1 y AA, en las que los grados de
pertenencia asociados a la misma clase son los mayores. Sin embargo, en el caso de
las alteraciones K2 (Potsica Feldesptica), SCC (Sericita-Clorita-Cuarzo) y S (Cuarzo-
Sericita), los grados de pertenencia mayores pueden ser compartidos y traslaparse
entre dos o tres alteraciones, como el cluster definido para SCC2, donde las muestras
mapeadas como SCC2 presentan el mayor grado de pertenencia, pero dentro del
mismo rango caen los valores de SCC1 y K2. Adems, las alteraciones flicas (S1 y S2)
presentan valores similares para los clusters definidos para ambas alteraciones.

106
Figura 7.5. Promedio de grados de pertenencia por alteracin mapeada.

La Figura 7.6 y Figura 7.7 muestran plantas de la zona de estudio de la mina


Escondida. El tamao del crculo asociado a las muestras se relaciona con el grado de
pertenencia al cluster indicado. Nuevamente, el cluster definido para K1 muestra una
alta correlacin con las muestras de dicha alteracin. Por otro lado, las alteraciones S1
y S2 muestran un comportamiento similar entre ellas, al igual que las SCC1 y SCC2. K2
es clasificada preferentemente como SCC2. La AA presenta un comportamiento
levemente distinto a las alteraciones flicas.

107
Figura 7.6. Planta Cota 2750 50 m. Tamao de las muestras representa el grado
de pertenencia a un cluster.

108
Figura 7.7. Planta Cota 2750 50 m. Tamao de las muestras representa el grado
de pertenencia a un cluster (continuacin).

109
7.4 Nmero ptimo de clusters

A pesar de que el error es un indicador de la calidad del clasificador, las matrices


de confusin muestran que existen alteraciones que pueden estar clasificadas dentro
de ms de 2 clases a la vez como se muestra en la Tabla 7.9.

Tabla 7.9. Matrices de confusin de Al-Mg-K-Ca-Na y Al-Mg-K-Na.

Al-Mg-K-Ca-Na
Alteracin Clasificada
K1 K2 SCC1 SCC2 S1 S2 AA
K1 89% 0% 8% 4% 0% 0% 0%
Alteracin Mapeada

K2 19% 14% 0% 48% 19% 0% 0%


SCC1 19% 1% 43% 34% 2% 0% 0%
SCC2 11% 6% 18% 58% 6% 0% 0%
S1 1% 1% 2% 15% 78% 2% 1%
S2 1% 0% 2% 9% 82% 5% 1%
AA 0% 0% 10% 20% 70% 0% 0%
Al-Mg-K-Na
Alteracin Clasificada
K1 K2 SCC1 SCC2 S1 S2 AA
K1 87% 0% 8% 6% 0% 0% 0%
Alteracin Mapeada

K2 14% 14% 0% 67% 5% 0% 0%


SCC1 15% 3% 49% 32% 1% 0% 0%
SCC2 9% 9% 21% 54% 7% 0% 0%
S1 0% 2% 1% 15% 76% 3% 2%
S2 0% 1% 2% 11% 77% 8% 2%
AA 0% 0% 10% 30% 60% 0% 0%

Las matrices de confusin muestran que las alteraciones flicas (S1 y S2) son
clasificadas principalmente como S1, pudiendo ser asignadas tambin al cluster SCC2.
Este comportamiento se repite en otras alteraciones, K2 que es clasificada
principalmente como SCC2, S1 y K1. Esto se explica principalmente al calcular las
distancias entre los centroides, las cuales muestran que las alteraciones con mayor
confusin se encuentran ms cercanas (ver Tabla 7.10).

110
Tabla 7.10. Distancia euclidiana entre centros de alteraciones (Al-Mg-K-Ca-Na).

K1 K2 SCC1 SCC2 S1 S2 AA
K1 - 3.52 0.27 0.50 4.46 4.64 4.27
K2 - - 3.42 1.78 0.78 0.86 1.00
SCC1 - - - 0.82 3.60 3.74 3.30
SCC2 - - - - 3.22 3.39 3.18
S1 - - - - - 0.00 0.04
S2 - - - - - - 0.04
AA - - - - - - -

De lo anterior, se plantea la pregunta de cuntos son efectivamente el nmero de


clusters propios de la base de datos.

Para encontrar el nmero ptimo de clusters, sin imponerlo previamente, se


plantea utilizar el coeficiente silhouette para el mtodo de clasificacin por particin: k-
means. Este algoritmo organiza los datos en una cantidad k definida de grupos o
clusters. Estos clusters se forman de tal manera que optimicen un criterio de particin
como una funcin de disimilitud basada en la distancia, tal que los datos dentro de un
cluster sean similares entre s y dismiles a los de los dems, en funcin de las variables
utilizadas. Cada cluster es caracterizado por su centroide, el cual corresponde al
promedio de los valores asignados a cada grupo.

En k-means, la funcin objetivo a optimizar es la suma de los errores cuadrticos


de los clusters. Al minimizar esta funcin, se intentan definir los clusters lo ms
compactos y separables entre s, es decir, que la distancia entre datos de un mismo
cluster sea mnima, y la distancia entre clusters sea la mxima posible. Esta funcin se
define de la siguiente forma:

E = dist(x, ci )2
i=1 x Ci

Donde x es un dato en el espacio, ci es el centroide del cluster Ci y dist(x, ci ) es


una funcin de distancia. En este caso, se usar la distancia euclidiana (Han et al.,
2012).

111
El algoritmo de particin k-means necesita como input un set de datos definidos
por una serie de variables, y un nmero entero definido de cluster k, y entrega como
output los clusters a los que pertenece cada muestra. La metodologa utilizada se
describe a continuacin (ver Figura 7.8):

i) Elegir k datos del set de datos como los centroides iniciales.


ii) Asignar cada dato al cluster al que se asemeje ms. En este caso se asigna
al cluster ms cercano utlizando la distancia definida.
iii) Actualizar los centroides, calculando el valor promedio de los valores
pertenecientes a cada cluster.
iv) Repetir pasos ii y iii hasta que no haya variacin en los centroides.

Figura 7.8. Metodologa de agrupacin de k-means (tomado de Han et al., 2012).

Una vez asignados los datos a cada cluster, es posible usar distintos parmetros
para evaluar las caractersticas de estos. Uno de los parmetros ms utilizados para
este mtodo es el valor de silhouttes (Pedregosa et al., 2011), el cual consiste en
calcular dos valores por cada dato, para finalmente obtener un valor promedio que
caracterice la clasificacin como se muestra a continuacin.

a: distancia promedio entre una muestra y todos los puntos en un cluster

b: distancia promedio entre una muestra y todos los puntos en el cluster ms


cercano.

s: silhouette de una muestra

ba
s=
max(a, b)

112
As se definido por las variables, obtiene un valor que evale la geometra de los
clusters en el espacio mediante el clculo del promedio de los valores de silhouette. En
el caso de los grupos Al-Mg-K-Ca-Na y Al-Mg-K-Ca, este valor alcanza su mximo para
un nmero de clusters igual a 3 (ver Figura 7.9).

Figura 7.9. Media de silhouettes vs nmero de clusters para grupos de elementos


mayores y trazas.

El agrupamiento de las alteraciones en los clusters calculados por el mtodo K-


means para los distintos grupos de elementos indica nuevamente que las alteraciones
S1 S2 y SCC1 SCC2 son afines (ver Tabla 7.11).

Tabla 7.11. Porcentaje de clasificacin de alteraciones en los clusters calculados


por K-means para Al-Mg-K-Ca-Na y Al-Mg-K-Ca.

Al-Mg-K-Ca-Na Al-Mg-K-Ca
C1 C2 C3 C1 C2 C3
K1 6% 4% 91% 5% 73% 22%
K2 91% 7% 2% 96% 0% 4%
SCC1 12% 0% 88% 9% 48% 43%
SCC2 36% 16% 47% 38% 26% 36%
S1 97% 1% 2% 97% 1% 2%
S2 98% 0% 2% 97% 0% 3%
AA 100% 0% 0% 95% 0% 5%

113
La distancia entre los centros de los clusters definidos para los grupos de
elementos mayores y los centros de los clusters predefinidos con las muestras
representativas de cada alteracin, muestran que para un mismo grupo de elementos
los centros ms cercanos entre s son aquellos que se clasifican dentro de un mismo
cluster (ver Tabla 7.12).

Tabla 7.12. Distancia entre centros predefinidos y calculados por K-means para
Al-Mg-K-Ca-Na y Al-Mg-K-Ca.

Al-Mg-K-Ca-Na Al-Mg-K-Ca
C1 C2 C3 C1 C2 C3
K1 2.12 2.61 0.42 2.04 0.48 0.68
K2 0.91 2.49 1.88 0.25 2.26 1.21
SCC1 1.9 2.96 0.11 1.89 0.54 0.51
SCC2 1.81 2.1 0.87 1.5 0.93 0.11
S1 0.03 3.33 1.96 0.02 2.43 1.39
S2 0.04 3.37 2 0.03 2.46 1.42
AA 0.2 3.39 1.89 0.2 2.34 1.3

De igual forma, la clasificacin supervisada muestra que las alteraciones pueden


ser clasificadas indistintamente como una u otra clase, como sucede con S1, S2 y AA.

Los estudios previos con esta base de datos (Townley et al., 2012), el anlisis de
las zonas de inters y el clculo de las silhouttes para K-means, indican que existen 3
grupos de datos principalmente. Segn la clasificacin propuesta, los grupos de
alteraciones hidrotermales se pueden definir de la siguiente manera.

El grupo de las alteraciones flicas, o grupo S, se define por S1, S2 y AA. El


grupo de las alteraciones clorita-sericita, o grupo SCC, se define por K2, SCC1 y SCC2.
El grupo K se define por la alteracin K1.

114
La alteracin K2, si bien muestra una menor distancia entre centros con S1 y S2,
posee un mayor porcentaje de clasificacin dentro de la alteracin SCC2 para los
grupos de elementos estudiados.

Segn esta nueva definicin de clases, las alteraciones hidrotermales quedan


distribuidas segn se muestra en la Tabla 7.13. Los 3 grupos de elementos
considerados, si bien tienden a clasificar las alteraciones de forma similar, difieren en
los porcentajes de datos. Por ejemplo, la alteracin potsica biottica (K1) muestra un
porcentaje muy alto para el primer grupo de elementos y uno bajo para la alteracin
potsica feldesptica (K2). Mientras que para el segundo grupo de elementos, la
tendencia es inversa.

Tabla 7.13. Porcentaje de clasificacin de alteraciones mapeadas en clusters


redefinidos para distintos grupos de elementos.

Al-Mg-K-Ca-Na Al-Mg-K-Ca Al-Mg-K-Na


Alteracin Mapeada K SCC S K SCC S K SCC S
K1 94% 6% 0% 77% 23% 0% 92% 8% 0%
K2 14% 76% 10% 10% 86% 5% 10% 86% 5%
SCC1 18% 81% 1% 16% 83% 1% 17% 82% 1%
SCC2 9% 86% 5% 9% 80% 11% 10% 83% 7%
S1 0% 20% 80% 0% 14% 85% 0% 20% 79%
S2 0% 13% 86% 0% 11% 89% 0% 15% 85%
AA 0% 40% 60% 0% 40% 60% 0% 40% 60%

7.5 Grados de pertenencia de muestras con una alteracin por tramo

Los grados de pertenencia calculados a partir de los datos con solo una alteracin
mapeada muestran los mayores valores para las alteraciones que fueron definidas. Es
decir, la alteracin K1 muestra sus mayores valores en el cluster K; SCC1, SCC2 y K2
al cluster SCC; y S1, S2 y AA al cluster S (ver Figura 7.10).

115
Figura 7.10. Boxplots de grados de pertenencia a clusters K, SCC y S calculados
para Al-Mg-K-Ca-Na para muestras con una alteracin por tramo.

116
La Figura 7.11 se construye calculando el porcentaje de clasificacin en 10 rangos
del grado de pertenencia, es decir, se consideran todos los datos pertenecientes a un
intervalo del grado de pertenencia a un cluster y se calcula el porcentaje de los datos
clasificados correctamente. Este grfico muestra el comportamiento creciente de la
distribucin de la clasificacin conforme aumenta el grado de pertenencia. La
distribucin de la clasificacin en funcin de los grados de pertenencia puede ser
utilizada para entregar informacin acerca de la incertidumbre con la que se asigna una
muestra a un cluster.

117
Figura 7.11. Grados de pertenencia calculados por rango para muestras con solo
una alteracin por tramo para diferentes grupos de elementos.

118
Un ejemplo del clculo del porcentaje de clasificacin en funcin de los grados
de pertenencia se da a continuacin, en donde se distingue entre 3 rangos de valores
separados por 2 inflexiones notorias para el grupo de elementos Al-Mg-K-Ca-Na. Cada
uno de estos rangos de valores presenta un porcentaje de clasificacin asignado de la
siguiente manera:

0, u < 0.4
% Clasificacin = {1.8 u 0.68, 0.4 u < 0.8
> 0.8, u 0.8

Donde u corresponde al grado de pertenencia asignado por la clasificacin. Esta


regla simple seala una de las utilidades del modelo propuesto utilizando las distancias
a los cluster de cada alteracin, pudiendo ser utilizado en el mapeo de sondajes
tomando mayor atencin a aquellos que posean un porcentaje de clasificacin menor.

7.6 Patrones mineralgicos

Con los modelos geoqumicos descritos se clasifican las muestras de mineraloga


avanzada Qemscan, y se analiza la razn de los errores con respecto al mapeo. Segn
las matrices de confusin, las alteraciones que tienden a ser asignadas errneamente a
otra clase son principalmente las alteraciones K1, SCC1 y SCC2. Dependiendo de los
elementos considerados en la clasificacin, entre 9 y 12 muestras de alteracin SCC1 y
SCC2 son asignadas al cluster K, 2 muestras de alteraciones K1 son asignadas al
cluster SCC, y 2 o 3 de las S1 y S2 al cluster SCC (ver Tabla 7.14).

119
Tabla 7.14. Matrices de confusin para los datos de mineraloga Qemscan.

Al-Mg-K-Ca-Na Al-Mg-K-Ca Al-Mg-K-Na


K SCC S K SCC S K SCC S Total
K1 2 2 0 2 2 0 2 2 0 4
Alt. Mapeada

SCC1 5 7 0 5 7 0 5 7 0 12
SCC2 6 10 0 7 9 0 4 12 0 16
S1 0 1 18 0 1 18 0 1 18 19
S2 0 2 21 0 1 22 0 2 21 23
AA 0 0 1 0 0 1 0 0 1 1

Los diagramas ternarios de elementos permiten agrupar las alteraciones en


funcin de los minerales predominantes. El tamao de los crculos asociadas a las
muestras indica la razn (Biotita+Clorita)/Pirita, y se asocian inversamente a la distancia
a los nodos definidos por la composicin de estos minerales (ver Figura 7.12).

120
Figura 7.12. Diagramas ternarios para las muestras Qemscan. El tamao de los
crculos indica mayor valor de la razn (Biotita+Clorita)/Pirita.

Utilizando los datos de mineraloga cuantitativa se comparan los patrones de


concentraciones de minerales de las muestras clasificadas dentro de los clusters
esperados y fuera de estos. El grupo de elementos utilizado para la clasificacin fueron
Al-Mg-K-Ca-Na. Debido a la cantidad de muestras, solo se analizar el caso de las
alteraciones K1 clasificadas como SCC, y las alteraciones SCC1 y SCC2 clasificadas
como K (ver Tabla 7.14).

121
La Figura 7.13 muestra las concentraciones de minerales Qemscan promedio de
los datos clasificados correctamente dentro de los cluster y de los clasificados
errneamente, segn el mapeo geolgico. La alteracin K1 posee concentraciones de
biotita mayores a las de clorita, mientras que ocurre lo contrario en SCC1 y SCC2. En
estas ltimas, la sericita es altamente predominante, mientras que en K1 presenta una
concentracin similar a la de biotita. Las muestras clasificadas dentro de otro cluster
presentan un patrn claramente distinto para K1 y SCC2. Para SCC1, el patrn de las
muestras con error es levemente menos pronunciado.

Figura 7.13. Composicin media de minerales para muestras clasificadas


correctamente y sus respectivos errores.

La Figura 7.14 compara las distribuciones de tres razones entre biotita, clorita y
sericita. En el caso de K1, los datos clasificados correctamente muestran un marcado
predominio de biotita+clorita sobre sericita, y de biotita sobre clorita, mientras que los
122
clasificados como SCC muestran el comportamiento opuesto. En SCC1, los datos mal
clasificados muestran un predominio de biotita+clorita sobre sericita. En SCC2,
muestran un predominio de biotita+clorita sobre sericita, y de biotita sobre clorita.

Figura 7.14. Boxplots composicin media de minerales para muestras clasificadas


correctamente y sus respectivos errores.

Los patrones mineralgicos de los datos pertenecientes a una alteracin


correctamente clasificada son distintos a los patrones de los datos clasificados
incorrectamente. Es decir, las relaciones mineralgicas de las muestras se ven
reflejadas directamente en su clasificacin segn el modelo propuesto para elementos
mayores. De esta manera, el error asociado a la clasificacin de datos Qemscan se
debera a que presentan un patrn de minerales de alteracin diferente a los de la clase
mapeada.

123
A partir de este anlisis se concluye que el modelo de clasificacin propuesto en
funcin de las concentraciones de elementos mayores es capaz de distinguir entre
patrones de proporciones modales de minerales. Finalmente, se evidencia que el
mapeo geolgico no siempre refleja la mineraloga de la roca.

7.7 Reclasificacin de alteracin flica indiferenciada

La alteracin flica indiferenciada (cdigo 50) corresponde a un cdigo de mapeo


antiguo, el cul no considera las posibles variaciones dentro de sta categora como,
por ejemplo, el contenido en arcilla o la razn clorita:sericita. Se reclasifican las
muestras con el modelo propuesto considerando el grupo de elementos Al-Mg-K-Ca-Na
(ver Tabla 7.15).

Tabla 7.15. Resultados de clasificacin de alteracin flica indiferenciada.

Cantidad de
datos
K 37
SCC 351
S 643
Total 1031

Si se compara la distribucin de los grados de pertenencia de esta clasificacin


(ver Figura 7.15) a los obtenidos anteriormente (ver Figura 7.11), se pueden asociar a
porcentajes de clasificacin de las alteraciones si se asume un comportamiento similar
de los datos. Es decir que es posible modelar el error de la clasificacin usando los
grados de pertenencia, ya que estos tienen un comportamiento inverso.

124
Figura 7.15. Boxplots de grados de pertenencia de alteracin flica indiferenciada.

Por otro lado, los diagramas ternarios ponen en evidencia que algunas muestras
poseen un mayor control asociado a biotita-clorita, en lugar de pirita y sericita que es lo
esperado para la alteracin flica (ver Figura 7.16).

125
Figura 7.16. Diagramas ternarios para alteracin flica indiferenciada.

Al comparar la distribucin espacial de la alteracin flica indiferenciada con la de


las muestras clasificadas anteriormente, se observa que el modelo tambin refleja una
correlacin en estas variables. Las muestras asignadas al tipo K se encuentran en la
zona NW, las SCC en la zona central y SW, y las tipo S en el Este (ver Figura 7.17 y
Figura 7.18).

126
Figura 7.17. Planta Cota= 2750 50 m. Grados de pertenencia a los clusters K,
SCC y S para Al-Mg-K-Ca-Na.

127
Figura 7.18. Distribucin espacial de clasificacin de alteracin flica
indiferenciada.

7.8 Grados de pertenencia de muestras con dos alteraciones por tramo

A diferencia de los datos con una alteracin, las muestras con dos alteraciones no
presentan una diferencia significativa entre los grados de pertenencia por cluster (ver
Figura 7.19). Este comportamiento es de esperar debido a que cada muestra contiene
las dos alteraciones en distintas proporciones, es decir, existe una mezcla tanto fsica
como qumica. Cabe destacar que esta separacin de la base de datos se realiza a
partir de la informacin de mapeo y no considera el efecto de la superposicin natural
de las alteraciones hidrotermales, en donde la roca presenta una serie de alteraciones
superpuestas una sobre la otra lo que se refleja en la presencia de minerales
caractersticos de cada una de ellas.

128
Figura 7.19. Boxplots de grados de pertenencia a clusters K, SCC y S calculados
para Al-Mg-K-Ca-Na para muestras con dos alteraciones por tramo.

Para analizar la clasificacin se separan los datos entre pares de alteraciones con
sus respectivas proporciones, por ejemplo para el par SCC S, una proporcin igual a
0.3 equivale a una muestra con 30% de alteracin S y 70% SCC (ver Figura 7.20).

Figura 7.20. Ejemplo del par SCC S.

Al igual que la seccin anterior, se calcula el porcentaje de clasificacin para los


pares de alteraciones por rango de proporcin entre ellas y grado de pertenencia para
una de las alteraciones. En este caso, el grupo de elementos utilizado es Al-Mg-K-Na
ya que es el que presenta un comportamiento ms definido (ver Anexo: Porcentaje de
clasificacin versus rangos de porcentaje de alteracin por tramo). Adems, el par K-S
no se incluye en el anlisis debido a que la clasificacin es tal que no contiene la
129
cantidad de datos necesarios que los vinculen, y por la confusin que existe entre las
alteraciones SCC y K, las cuales tienden a ser asignadas con gran error.

El porcentaje de clasificacin entre pares de alteraciones alcanza un mnimo


cuando la diferencia en la proporcin de alteraciones por sondaje disminuye, y aumenta
cuando una de las dos se hace predominante. Debido a que el modelo tiende a
confundir entre las alteraciones K y SCC, el porcentaje de clasificacin disminuye
cuando una de ellas est en mayor proporcin. Por otro lado, la clasificacin presenta
un comportamiento similar en funcin de los grados de pertenencia, pudiendo tener
mayores variaciones a medida que la diferencia entre los grados de pertenencia es
menor (ver Figura 7.21).

Figura 7.21. Porcentaje de clasificacin en funcin de las proporciones de


alteracin por tramo y grado de pertenencia para los pares K SCC y SCC S,
considerando el grupo de elementos Al-Mg-K-Na.

130
8 Discusin

El modelo de clasificacin propuesto permite distinguir entre 3 zonas de alteracin:


K, SCC y S. La clasificacin considera los elementos mayores que se distribuyen en
funcin de las alteraciones. A su vez, estos elementos se correlacionan directamente
con la presencia de minerales solubles en agua regia: biotita, clorita, xidos de hierro y
pirita, principalmente.

En el cluster K, el principal control mineral lo ejercen biotita y clorita sobre sericita, y


biotita por sobre clorita. En S, la sericita y pirita predominan ampliamente por sobre los
dems minerales. En SCC, la clorita y sericita ejercen el control principal de la
mineraloga de alteracin, pudiendo presentar distintas proporciones entre ellas.

La seleccin de variables mediante el clculo de traslapes entre las alteraciones


hidrotermales, permite concluir que existen grupos de alteraciones que no pueden ser
diferenciados a partir de la geoqumica de elementos solubles en agua regia. Estos
grupos corresponden a SCC1 SCC2 y S1 S2 AA. Adems, si bien la diferencia de
mapeo entre estas alteraciones es la presencia de arcillas, los datos de mineraloga
avanzada indican que ni la caolinita o la pirofilita presentan una particin preferencial en
ninguna de estas alteraciones estudiadas (ver Figura 8.1).

131
Figura 8.1. Boxplots de concentracin de arcillas por alteracin.

Los datos Qemscan muestran que los minerales predominantes en estas


alteraciones son cuarzo y sericita, indicando que la distincin entre ellas en el protocolo
de mapeo no se corresponde totalmente con la mineraloga cuantitativa. De esta
manera, se pone en duda el supuesto de que el mapeo geolgico determina la
alteracin hidrotermal predominante de la roca, lo que introduce un significativo grado
de error a la clasificacin.

Estos errores son esperables debido a que comnmente en el mapeo de sondajes


no se considera la sobreimposicin de alteraciones. Por ejemplo, un tramo de sondaje
puede ser mapeado como alteracin potsica feldesptica, debido a que presenta
feldespatos potsicos de gran tamao, pero no se considera que puedan estar
altamente sericitizados, lo que la asociara ms a una alteracin flica. Es decir, el
anlisis visual que pasa por alto este tipo de detalles, puede ser incongruente con el
comportamiento geoqumico. Un ejemplo de esto es el caso del sondaje D-5233
mostrado anteriormente. Este sondaje presenta muestras mapeadas como alteracin
flica S2 con concentraciones altas de Al y Mg, por lo que el modelo de clasificacin las
asigna al cluster SCC (ver Figura 8.2).

132
Figura 8.2. Sondaje con mapeos de alteracin incongruente con geoqumica.

Adems, el modelo clasifica las muestras obedeciendo a los patrones de razones


entre minerales de alteracin como se demuestra en la seccin 7.6. En esta seccin, las
muestras de mineraloga Qemscan clasificadas errneamente segn el mapeo
presentan razones de biotita, clorita y sericita distintas a las definidas para cada
alteracin. De esta manera, el modelo de clasificacin tiene una aplicacin en el mapeo
de sondajes, al sugerir alteraciones en funcin de la qumica, en zonas donde exista
mayor incertidumbre para el gelogo.

Las zonas de contacto entre alteraciones pueden presentar distintos patrones de


transicin segn la geoqumica: progresiva o a saltos. Esto depender de si es que
existe o no una gradacin de las razones de minerales. En el caso de los contactos
entre los clusters K SCC, la transicin se dara cuando la clorita comience a
predominar sobre la biotita, y los xidos de hierro, como magnetita, disminuyan su
proporcin. Sin embargo, para el caso de la transicin SCC-S, el lmite se vuelve ms
difuso ya que, si bien las primeras poseen mayor proporcin de clorita que la segunda,
en ambos casos el cuarzo y la sericita tienden a ser los minerales predominantes. De
esta manera, los mayores errores de clasificacin ocurren principalmente en las
alteraciones K2, SCC1 y SCC2, asignadas al cluster SCC, y las alteraciones S1 y S2
asignadas al cluster S (ver Tabla 8.1). Este error de mapeo se vincula directamente a
las concentraciones de los elementos qumicos provenientes de los minerales solubles
133
en agua regia. Es por ello, que las alteraciones sericita-clorita pueden ser clasificadas
como potsicas biotticas, si las concentraciones de Al-Mg-K son relativamente altas o
como flica, si son bajas.

Tabla 8.1. Porcentaje de clasificacin de alteraciones mapeadas en clusters


redefinidos para distintos grupos de elementos.

Al-Mg-K-Ca-Na Al-Mg-K-Ca Al-Mg-K-Na


Alteracin Mapeada K SCC S K SCC S K SCC S
K1 94% 6% 0% 77% 23% 0% 92% 8% 0%
K2 14% 76% 10% 10% 86% 5% 10% 86% 5%
SCC1 18% 81% 1% 16% 83% 1% 17% 82% 1%
SCC2 9% 86% 5% 9% 80% 11% 10% 83% 7%
S1 0% 20% 80% 0% 14% 85% 0% 20% 79%
S2 0% 13% 86% 0% 11% 89% 0% 15% 85%
AA 0% 40% 60% 0% 40% 60% 0% 40% 60%

Una aproximacin para la definicin del contacto entre las alteraciones SCC S, se
muestra en la seccin 7.8, donde se calculan los porcentajes de clasificacin por rangos
de proporcin y por rangos de grado de pertenencia entre ambas alteraciones. Los
resultados de este anlisis indican que cuando el grado de pertenencia se encuentra
entre 30-40% del cluster S, la clasificacin es mnima, y aumenta a medida que una de
las dos se hace predominante. Para el caso de las proporciones fsicas entre
alteraciones, no es claro dnde se encuentra el lmite entre ambas, pero a medida que
la proporcin de S aumenta considerablemente, el porcentaje de clasificacin aumenta
por sobre el 70%.

En el caso de los lmites entre K SCC, el comportamiento del porcentaje de


clasificacin es similar pero no tan pronunciado. Por un lado, cuando el grado de
pertenencia a SCC aumenta sobre un 70%, la clasificacin lo hace junto con ella por
sobre un 80%. Mientras que la clasificacin se mantiene entre un 40% y 70% para
cualquier rango de proporcin fsica entre alteraciones. Esto puede atribuirse a un error
en el mapeo, ya que como se evidencia en la clasificacin de los datos Qemscan, la

134
proporcin biotita clorita no siempre est correctamente determinada por el ojo
geolgico. Es por ello, que estas dos alteraciones tienden a ser mapeadas
indistintamente, ya que el reemplazo de biotita por clorita no siempre es evidente a
escala de muestra de mano.

Los grados de pertenencia para un grupo de elementos determinado entregan


informacin de cun cercanos son los datos a los centroides definidos para una
alteracin en particular. As, es posible asociar los grados de pertenencia a un grado de
incertidumbre. Se tendr una incertidumbre alta cuando los grados de pertenencia sean
muy similares entre s, y una baja cuando uno de estos valores sea considerablemente
mayor a los dems. En los siguientes perfiles y secciones, se destacan aquellas
muestras que poseen un mayor grado de incertidumbre asociado al tamao del crculo
que las representa. Un mayor tamao indica mayor incertidumbre debido a que la
diferencia entre los grados de pertenencia mximo e intermedio es baja. Como
evidencian las figuras, las zonas con mayor incertidumbre corresponderan a la central y
Este, en donde se ubican principalmente las alteraciones sericita-clorita y flica,
respectivamente. Esta incertidumbre puede deberse a las razones de minerales de
biotita-clorita con respecto a sericita, los cuales pueden estar en mayor o menor
proporcin debido a la propia litologa o a una alteracin clortica propiamente tal. Cabe
mencionar que la alteracin flica puede presentar clorita en proporciones considerable
(ver Figura 8.3, Figura 8.4 y Figura 8.5).

Se destaca que la clasificacin emula la distribucin espacial de las alteraciones, la


cual se caracteriza por una zona de alteracin potsica biottica al NW, una zona de
alteracin flica al E, y una zona de alteracin sericita-clorita en la zona central (ver
Figura 8.6).

135
Figura 8.3. Distribucin espacial de alteraciones hidrotermales. Tamao indica
grado de incertidumbre. Perfil N = 107350 50 m. Planta cota = 2750 50 m.

Figura 8.4. Distribucin espacial de litologas. Tamao indica grado de


incertidumbre. Perfil N = 107350 50 m. Planta cota = 2750 50 m.
136
Figura 8.5. Distribucin espacial de zonas de mineralizacin. Tamao indica grado
de incertidumbre. Perfil N = 107350 50 m. Planta cota = 2750 50 m.

137
Figura 8.6. Distribucin espacial de clasificacin de alteraciones hidrotermales.
Tamao indica grado de incertidumbre. Perfil N = 107350 50 m. Planta cota =
2750 50 m.

La metodologa puede ser utilizada de varias maneras tiles considerando que las
diferencias entre las clases predefinidas se relacionan a los minerales solubles en agua
regia. Cada clase que se defina debera ser estudiada considerando las variables
crticas. Desde el punto de vista del procesamiento minero, se deben considerar todas
las etapas para estudiar las caractersticas que definen cada clase y las hace distintas
para su tratamiento. Un enfoque integral considerara las variables mineralgicas,
geoqumicas, texturales, geotcnicas y metalrgicas para definir zonas con distintos
protocolos de procesamiento para la recuperacin de mena. En particular, los
filosilicatos son el grupo mineral que causa mayor impacto en el proceso de flotacin,
especialmente la muscovita y arcillas como la vermiculita y montmorillonita. Su
presencia puede causar: disminucin de la recuperacin por la adhesin de ganga fina
sobre partculas mineralizadas o en las burbujas; aumento del consumo de reactivos
debido a alta rea superficie de los minerales; y floculacin en celdas de flotacin
debido a que la presencia de arcillas disminuye la capacidad espumante (Bulatovic
2007; Jorjani et al., 2011; Ndlovu et al., 2013).

138
8.1 Limitantes

El modelo no considera una serie de variables que pueden afectar a la clasificacin


como estructuras y texturas, principalmente. Por ejemplo, una roca puede tener una
alteracin biottica de fondo con una densidad alta de vetillas de pirita con halo
serictico, lo que pone en duda cul alteracin es predominante y cul subordinada
desde el punto de vista qumico. Para las texturas el efecto sera similar. Una textura de
reemplazo de minerales o inclusin de otras litologas, como sera el caso de las
brechas (no incluidas en el modelo), puede llevar a que una muestra sea clasificada
dentro de una clase totalmente distinta. Por ejemplo, en una brecha hidrotermal con
matriz de turmalina, los minerales que se disuelven en agua regia corresponderan
nicamente a los que se encuentren en los fragmentos y no a la matriz, lo que
conducira a clasificar la muestra segn los fragmentos predominantes. Es por ello, que
el modelo necesita ser calibrado segn las clases esperadas y ser utilizado como
complemento del mapeo.

A pesar de que parezca una desventaja que la base de datos haya sido obtenida
en agua regia y no disuelva la mayora de los minerales, este hecho acota las posibles
fuentes de los elementos. Una disolucin mediante otros mtodos como 4-cidos o
cido fluorhdirico habra significado que los elementos varen considerablemente sus
concentraciones y tratar determinar su proveniencia sera un problema litogeoqumico
ms complejo.

Los grupos minerales solubles en agua regia son sulfuros, micas trioctadricas,
sulfatos, carbonatos y oxi-hidrxidos de Fe-Mn. Los dems grupos minerales, como las
micas blancas, feldespatos, hornblenda y caolinita, mantienen sus proporciones
constantes posteriores al agua regia. Otros minerales, como la albita, pueden disolverse
solo en fracciones muy bajas por lo que la proveniencia de los elementos disueltos no
puede ser determinada directamente de la presencia de estos (Ryan et al., 2002; Ryan
et al., 2008; Tarvainen et al., 2009; Macfarlane et al., 2005).

139
El grupo mineral soluble en agua regia a partir del cual es posible definir
alteraciones hidrotermales son los filosilicatos. Las micas y arcillas pertenecientes al
grupo de los filosilicatos presentan una estructura cristalina consistente en capas de
cationes coordinados con O2- y/u OH-. Estas capas son tetradricas y octadricas. Las
tetradricas consisten en Si4+ o Al3+ en coordinacin tetradrica con O2- u OH-. Las
octadricas corresponden a cationes, tpicamente a Al, Fe o Mg, en coordinacin
octadrica (nmero de coordinacin 6) con O2- u OH-. Entre las capas, pueden
encontrarse por cationes adsorbidos o molculas de agua. Las micas pueden ser
clasificadas en funcin de la proporcin de los sitios octadricos ocupados. Si 2 de cada
3 sitios octadricos estn ocupados se dice que la mica es dioctadrica, y si se
encuentran totalmente ocupados, entonces se denomina trioctadrica (Langmuir, 1997).

Dentro del grupo de las micas trioctadricas, se destaca el comportamiento de las


biotitas, cloritas y esmectitas, dado que son disueltas totalmente y se asocian a la
liberacin de Mg, Fe y Al, lo que indicara una disolucin preferencial de las capas
octadricas por sobre las tetradricas. Adems, las cloritas y esmcticas ricas en Mg
presentan una tasa de disolucin mayor que las ricas en Fe (II) y Al, respectivamente.
Esto se debera a que los enlaces Mg-O son ms dbiles que los Fe-O o Al-O en este
tipo de minerales (Ryan et al., 2008).

Los minerales que poseen la mayor proporcin de K son los feldespatos potsicos
y filosilicatos. Sin embargo, la fraccin extrable por agua regia corresponde a los
filosilicatos excluyendo a las micas blancas, es decir, el K provendra de micas
trioctadricas, illita y filosilicatos de capa mixta (Andrist-Rangel, et al., 2010). En el
anlisis de la mineraloga y la geoqumica, las altas concentraciones de K se
correlacionan a la presencia de biotita y no a la de sericita, por lo que se corrobora la
proveniencia principal de este elemento.

El conocimiento de la proveniencia de los elementos qumicos es de vital


importancia para el anlisis de alteraciones porque permite asociar los datos
geoqumicos con zonas de enriquecimiento relativo de minerales. De lo contrario, las
alteraciones pueden ser definidas errneamente al asociar elementos a minerales
equivocados. Un ejemplo de esto se tiene al revisar parmetros artificiales como K/Al y
Mg/Al (ver Figura 8.7), donde las muestras con alteracin S1, S2 y AA grafican
140
cercanas al nodo definido para sericita, y las alteraciones K1, SCC1 y SCC2 lo hacen
cercanas al nodo de clorita. A pesar de que los histogramas y la mineraloga
cuantitativa muestran que K proviene de las alteraciones K1, SCC1 y SCC2
principalmente, el grfico da a entender que la sericita es una de las fuentes de este
elemento. Sin embargo, esto se debe a que el rango de variacin del aluminio es
mucho mayor que el del potasio y por lo tanto la razn K/Al muestra valores altos para
las alteraciones flicas. En este caso la conclusin es correcta pero la hiptesis de la
fuente del elemento es errnea, mostrando que los parmetros artificiales pueden tener
un comportamiento opuesto al esperado.

Figura 8.7. Razn K/Al vs Mg/Al.

9 Conclusiones y recomendaciones

La correlacin entre la mineraloga, geoqumica y mapeo geolgico permite


identificar zonas de alteracin hidrotermal. Una vez caracterizada cada zona, es posible
definir zonas de alteracin en funcin de la geoqumica. Los datos geoqumicos se
correlacionan con la presencia de los minerales solubles en agua regia como: biotita,
clorita, xidos de hierro, sulfuros, sulfatos y carbonatos, principalmente.
141
Las variables geoqumicas que permiten diferenciar entre alteraciones son los
elementos mayores ya que tienen directa relacin con la mineraloga. Estos elementos
son Al, Mg, K, Ca, Na, S, Fe y CuT. Algunos elementos traza presentan un
comportamiento similar en funcin de las alteraciones, pero su presencia no se
correlaciona directamente con la mineraloga.

El anlisis de traslape entre variables y la clasificacin de muestras a partir de la


definicin de clusters por alteracin, indican que existen alteraciones que no pueden ser
diferenciadas de otras a partir de la geoqumica como: S1 S2 AA y SCC1 SCC2.
Esto puede tener dos explicaciones posibles no excluyentes. La primera se debera a
que el agua regia no disuelve las arcillas que corresponden a los minerales que las
diferencian por definicin en el protocolo de mapeo. La segunda, es que no existe una
diferencia mineralgica entre ellas, como muestran los datos de mineraloga cuantitativa
Qemscan, ya que ninguna de las alteraciones definidas por la presencia de arcillas
presenta concentraciones significativamente mayores que las que no. Se recomienda
realizar un anlisis de las variables mineralgicas previo a la definicin de clases de
alteraciones hidrotermales con la finalidad de disminuir la incertidumbre en el mapeo,
sobre todo en minerales que son difciles de diagnosticar como las arcillas.

Los datos geoqumicos y mineralgicos indican que las alteraciones hidrotermales


tendran una evolucin propia de un sistema de tipo prfido cuprfero caracterizadas por
un ncleo potsico con sobreimposicin flica a medida que la distancia aumenta a
partir de ste. La alteracin SCC es una zona de alteracin caracterizada por
variaciones de la razn clorita:sericita con la distancia al ncleo potsico, aumentando
la predominancia de la sericita (Sillitoe, 2010). Este comportamiento se aprecia en la
unidad de lavas andesticas, en la cual se predomina la alteracin flica a medida que
aumenta la distancia de la unidad de prfido feldesptico.

El modelo propuesto permite distinguir entre tres zonas de alteracin hidrotermal


en funcin de la mineraloga soluble dominante. La alteracin potsica se distingue con
facilidad de las alteraciones flicas. Sin embargo, las alteraciones clorita-sericita
presentan el mayor problema de clasificacin ya que dependiendo del predominio de
uno de estos dos minerales, pueden ser clasificadas como flicas o potsicas. Adems,
el modelo permite caracterizar cada muestra en funcin de diversos controles minerales
142
utilizando diagramas ternarios, pudiendo diferenciar entre muestras con control de
biotita-clorita, sericita, sulfatos, pirita y otros sulfuros de cobre.

Una manera de mejorar la clasificacin es considerar las muestras espacialmente


ms cercanas clasificadas con menor error. Este error puede ser calculado utilizando
los rangos de grados de pertenencia asociados a un porcentaje de clasificacin. En
trabajos futuros se recomienda la utilizacin de puntos de control, y la utilizacin de la
distribucin espacial de las muestras para mejorar la clasificacin.

Esta metodologa puede ser utilizada para obtener un ndice de alteracin, es decir,
para estimar cun alterada se encuentra una roca. Para ello, se debe considerar una
lnea base o background, es decir, es necesario contar con la caracterizacin
geoqumica y mineralgica de las litologas que no presenten alteracin para poder
calcular un ndice de cun alterada est la roca con respecto a su estado no alterado.
Una forma de hacer esta estimacin sera definir un cluster con las litologas frescas y
utilizar los grados de pertenencia a los dems cluster de alteraciones.

La aplicacin principal de este modelo es sugerir una alteracin hidrotermal en lugar


de determinarla, ya que existen una serie de parmetros que no se consideran dentro
de ste anlisis, y puede conducir a errores si se utiliza en zonas para las que no se ha
calibrado. Por ello, el modelo entrega informacin acerca del enriquecimiento o
empobrecimiento relativo de los minerales en la roca utilizando los diagramas ternarios
de elementos, adems de una alteracin probable dadas las clases previamente
definidas.

Se recomienda revisin detallada de las muestras cuya clasificacin geoqumica sea


inconsistente con el mapeo geolgico para determinar la fuente de las variaciones de
los elementos mayores.

Como en petrografa descriptiva, la identificacin de minerales no se realiza


mediante solo una de sus caractersticas, sino que es la suma de ellas la que permite
determinar su presencia en la roca.

143
10 Bibliografa
Amilibia, A. y otros, 2008. The role of inherited tectono-sedimentary architecture in the
development of the central Andean mountain belt: Insight from Cordillera de Domeyko.
Journal of Structural Geology, pp. 1520-1539.

Andrist-Rangel, Y. y otros, 2010. Assesing potassium reserves in northern temperate


grassland soils: A perspective based on quantitative mineralogical analysis and aqua-
regia extractable potassium. Geoderma, pp. 303-314.

BHP Billiton, 2012. Annual Report, s.l.: s.n.

Bulatovic, S., 2007. Handbook of Flotation Reagents Chemistry, Theory and Practice:
Flotation of Sulphide Ores. Elsevier.

Camus Infanta, F., 2003. La Faja del Eoceno Superior - Oligoceno. En: Geologa de los
Sistemas Porfricos en los Andes de Chile. Santiago: Servicio Nacional de Geologa y
Minera, pp. 196-201.

Camus, F. & Dilles, J. H., 2001. A Special Issue Devoted to Porphyry Copper Deposits
of Northern Chile. Economic Geology, Volumen 96, pp. 233-237.

Chao, T. T., 1984. Use of partial dissolution techniques in geochemical exploration.


Journal of Geochemical, pp. 101-135.

Chao, T. T. & Sanzalone, R. F., 1992. Decomposition techniques. Journal of


Geochemical Exploration, pp. 65-106.

Charrier, R., Faras, M. & Maksaev, V., 2009. Evolucin tectnica, paleogeogrfica y
metalognica durante el cenozoico en los Andes de Chile norte y central e
implicaciones para las regiones adyacentes de Bolivia y Argentina. Revista de la
Asociacin Geolgica Argentina, Volumen 65, pp. 5-35.

Church, S. E., Mosier, E. L. & Motooka, J. M., 1987. Mineralogical Basis for the
Interpretation of Multi-Element (ICP-AES), Oxalic Acid, and Aqua Regia Partial
Digestions of Stream Sediments for Reconnaissance Exploration Geochemistry. Journal
of Geochemical Exploration, pp. 207-233.

144
Corbett, G. & Leach, T., 1997. Short course manual: Southwest Pacific rim gold-copper
systems: Structure, alteration and mineralization, s.l.: s.n.

Cornejo, P., Mpodozis, C., Rivera, O. & Matthews, S., 2009. Carta Exploradora. Regin
de Antofagasta y Atacama. Carta Geolgica de Chile. Serie Geolgica Bsica, Issue
119, pp. 10-11.

Falasconi, M. y otros, 2010. A stability based validity method for fuzzy clustering. Pattern
Recognition, Volumen 43, p. 12921305.

Gallego, O. F. & Covacevich, V., 1998. Conchostracos trisicos de las Regiones de


Antofagasta, Atacama y Coquimbo, Chile. Revista geolgica de Chile, 25(2).

Gill, R., 1997. Modern Analytical Geochemistry. An introduction to quantitative chemical


analysis for earth, environmental and materials scientists. London: Longman.

Hall, G., 1999. Near-total acid digestions. Explore, Volumen 104, pp. 15-19.

Hall, G. E. & Cameron, E. M., 2001. Comparison of Selective Extractions of Soils Over
The Spence Deposit, Chile. Part A: Sequential Extractions of a Soil Profile, s.l.: s.n.

Han, J., Kamber, M. & Pei, J., 2012. Data Mining: Chapter10 Cluster Analysis Basic
Concepts and Methods. s.l.:Elsevier.

Hofmann, A., 2003. Sampling Mantle Heterogeneity through Oceanic Basalts: Isotopes
and Trace Elements. En: A. M. Davis, ed. Treatise On Geochemistry. s.l.:Elsevier
Pergamon.

Jorjani, E., Barkhordari, H., Tayebi Khorami, M. & Fazeli, A., 2011. Effects of
aluminosilicate minerals on coppermolybdenum flotation from Sarcheshmeh porphyry
copper. Minerals Engineering, Volumen 24, p. 754759.

Kelley, D. L. y otros, 2003. The use of partial extraction geochemistry for copper
exploration in northern Chile. Geochemistry: Exploration, Environment, Analysis,
Volumen 3, pp. 85-104.

145
Klassen, R. A., 2001. The interpretation of background variation in regional geochemical
surveys - an example from Nunavut, Canada. Geochemistry: Exploration, Environment,
Analysis, pp. 163-173.

Langmuir, D., 1997. Aqueous environmental geochemistry. New Jersey: Prentice-Hall.

Liu, H.-C., Jeng, B.-C., Yih, J.-M. & Yu, Y.-K., 2009. Fuzzy C-Means Algorithm Based on
Standard Mahalanobis Distances. Proceedings of the 2009 International Symposium on
Information Processing, pp. 422-427.

Luca, R., Lpez, L. F. & Townley, B., 2012. Geoqumica de Asociaciones de Alteracin
Hidrotermal en el Sistema tipo Prfido Cuprfero Escondida, Santiago: s.n.

Macfarlane, W. R. y otros, 2005. Continuous leach inductively coupled plasma mass


spectrometry: applications for exploration and environmental geochemistry.
Geochemistry: Exploration, Environment, Analysis, Volumen 5, pp. 123-134.

Maksaev, V., 2001. Resea Metalognica de Chile y de los Procesos que Determinan la
Metalognesis Andina, s.l.: s.n.

Marinovic, N. y otros, 1995. Hoja Aguas Blancas. Regin de Antofagasta, Santiago:


Servicio Nacional de Geologa y Minera.

Nazari, M., Shanbehzadeh, J. & Sarrafzadeh, A., 2013. Fuzzy C-means based on
Automated Variable Feature Weighting. Proceedings of the International
MultiConference of Engineers and Computer Scientists, Volumen I.

Ndlovu, B., Farrokhpay, S. & Bradshaw, D., 2013. The effect of phyllosilicate minerals
on mineral processing industry. International Journal of Mineral Processing, Volumen
125, p. 149156.

Padilla-Garza, R. A., Titley, S. R. & Eastoe, C. J., 2004. Hypogene Evolution of the
Escondida Copper Deposit, Chile. Society of Economic Geologists, Volumen Special
Publication, pp. 141-165.

Padilla, R. A., Titley, S. & Pimentel, F., 2001. Geology of the Escondida Porphyry
Copper Deposit, Antofagasta Region, Chile. Economic Geology, Volumen 96, pp. 307-
324.
146
Pedregosa, F. y otros, 2011. Scikit-learn: Machine Learning in Python. Journal of
Machine Learning Researc, pp. 2825-2830.

Pimentel, B. & M.C.R de Souza, R., 2013. A multivariate fuzzy c-means method. Applied
Soft Computing, Volumen 13, p. 15921607.

Pirrie, D. & Rollinson, G. K., 2011. Unlocking the applications of automated. The
Geologists Association & The Geological Society of London, Geology Today,
Noviembre-Diciembre.27(6).

Rezaee, B., 2010. A cluster validity index for fuzzy clustering. Fuzzy sets and systems,
Volumen 161, p. 30143025.

Richards, J., Boyce, A. & Pringle, M., 2001. Geologic Evolution of the Escondida Area,
Northern Chile: A Model for Spatial and Temporal Localization of Porphyry Cu
Mineralization. Economic Geology, Volumen 96, pp. 271-305.

Richards, J. P., Noble, S. R. & Pringle, M. S., 1999. A Revised Late Eocene Age for
Porphyry Cu Magmatism in the Escondida Area, Northern Chile. Economic Geology,
Volumen 94, pp. 1231-1248.

Rojas, J., Chavarro, J. C. & Moreno, R., 2008. Tcnicas de lgica difusa aplicadas a la
minera de datos. Scientia et Technica Ao XIII, Universidad Tecnolgica de Pereira.

Romero, B. y otros, 2010. Molybdenite Mineralization and Re-Os Geochornology of the


Escondida and Escondida Norte Porphyry Deposits, Northern Chile. Resource Geology,
61(1), pp. 91-100.

Ryan, P., Hillier, S. & Wall, A., 2008. Stepwise effects of the BCR sequential chemical
extraction procedure on dissolution and metal release from common ferromagnesian
clay minerals: A combined solution chemistry and X-ray powder diffraction study.
Science of the total environment, pp. 603-614.

Ryan, P., Wall, A., Hillier, S. & Clark, L., 2002. Insights into sequential chemical
extraction procedures from quantitative XRD: a study of trace metal partitioning in
sediments related to frog malformities. Chemical Geology, Volumen 184, pp. 337-357.

147
Sillitoe, R., 2010. Porphyry Copper Systems. Economic Geology, Volumen 105, pp. 3-
41.

Sillitoe, R. H. & McKee, E. H., 1996. Age of Supergene Oxidation and Enrichment in the
Chilean Porphyry Copper Province. Economic Geology, Volumen 91, pp. 164-179.

Snll, S. & Liljefors, T., 2000. Leachability of major elements from mineral stron acids.
Journal of Geochemical Exploration, pp. 1-12.

Tarvainen, T., Jarva, J. & Kahelin, H., 2009. Geochemical baselines in relation to
analytical methods in the It-Uusimaa and Pirkanmaa regions, Finlad. Geochemistry:
Exploration, Environment, Analysis, Volumen 9, pp. 81-92.

Townley, B., Muoz, M. & Luca, R., 2012. Modelamiento Geoqumico Distrito
Escondida: Discriminacin de Ambientes Hidrotermales. Fase 2: Estudio de Distribucin
Espacial y Determinacin de Criterios de Discriminacin de Ambientes Hidrotermales en
el Yacimiento Escondida, Santiago: UNTEC.

Vliz, W., 2004. Relacin espacio-temporal del sistema prfido cuprfero y epitermal en
el yacimiento Escondida, Provincia de Antofagasta, Segunda Regin, Chile.Tesis para
optar al grado de magster en geologa econmica mencin exploracin., Antofagasta:
Universidad Catlica del Norte.

Xing, H.-J. & Ha, M.-H., 2014. Further improvements in Feature-Weighted Fuzzy C-
Means. Information Sciences.

Xu, G. y otros, 2012. Digestion methods for trace element measurements in shales:
Paleoredox proxies examined. Chemical Geology, Volumen 324-325, pp. 132-147.

148
Anexo A
Base de datos

Cantidad de datos por litologa, alteracin y minzone

N 1 N 2
Litologa Abrev Cdigo
litologa litologa
Relleno Artificial RA 1 1 4
Relleno Acopio RS 2 0 0
Grava GR 6 4 2
Roca Moteada RM 10 0 0
Dique (en el sentido amplio) D 11 0 0
Brecha Hidrotermal BH 21 287 197
Brecha Tectnica BT 22 42 68
Brecha Magmtica-Ignea-
BI 23 337 254
Contacto
Brecha Pebbles Dikes BP 24 0 0
Prfido Feldesptico PF 31 5888 2016
Prfido Rioltico (ENorte) PR 32 7 3
Prfido Cuarcfero (Mina) PC 33 1824 484
Prfido Grueso (ENorte) PG 34 0 0
Prfido Dactico (tardia) DT 35 19 22
Diorita DR 40 6 10
Andesita AN 50 6690 2617
Autobrecha (Volcnica) AB 51 24 16
Tobas Dacticas TD 53 0 0
Tobas Riolticas TR 54 3 1
Arenisca Volcnica AV 64 0 0
Porfido Feldesptico Tardio PFT 70 18 2
Prfido Negro PNO 71 46 10
Porfido Dactico de Anfbol PDA 72 0 0
Porfido Temprano PT 73 48 9
Unidad Volcano-
UVS 81 0 0
Sedimentaria

149
Alteracin Abrev. Cdigo Lito 31 Lito 33 Lito 50
Fresco: sin Alteracin F 10 1 0 1
Propiltico P 20 0 0 5
Potsico (Slo Bt secundaria) K1 30 9 0 504
Feld K > Bt K2 31 184 0 16
Feld K < Bt K3 32 0 0 2
(Bt sec +/- Feld K) > Cl K4 33 3 0 0
Clorita-Sericita-Arcillas SCC1 40 106 6 1794
Clorita-Sericita-Cuarzo SCC2 41 790 66 1107
Clorita-Biotita +/- Feld K SCC3 42 3 0 4
Clorita SCC4 43 0 0 44
Flica S 50 1100 92 1192
Sericita-Cuarzo S1 51 2314 1160 883
Sericita-Cuarzo-Arcilla S2 52 1270 500 992
Sericita-Cuarzo S3 53 25 0 0
Sericita Gris Verde S4 54 7 0 0
Argilizacin Suprgena AA 61 62 0 123
Argilizacin Avanzada AAV 62 4 0 7
Silicificacin Q 70 10 0 16

150
1 2
Zonas mineralgicas Abrev. Cdigo
Minzone Minzone
Lixiviado LX 10 313 224
Oxidos Verdes (Sulfatos) OXV 21 16 10
Oxidos Azules (Silicatos) OXA 22 0 0
Oxidos Negros OXN 23 0 0
Cuprita + Ox Cu CPOX 24 12 14
Oxidos Verdes+Negro+Arcillas
OXVN 25 0 2
Cuprferas
Cuprita / Cobre Nativo CPCu 26 1 1
Fosfatos + Oxidos Verdes POXV 27 0 0
Parcial Lixiviado (Oxidacin
Plx 31 46 28
Parcial)
Parcial Lixiviado (Limonitas
PlxE 32 78 54
Exticas)
Mezcla OxCu + Sulf MX 40 81 65
Mixto + Cuprita y/o Nativo CPMX 41 22 24
Slfuro + Cuprita y/o Nativo CPCCPY 42 5 11
Cc + Py HE1 50 2061 702
Cc + Cv + Py HE2 51 891 344
Cv + Py HE3 52 81 55
Cc + Cpy + Py LE1 60 928 451
Cc + Cv + Cpy + Py LE2 61 2446 1021
Cv + Cpy + Py LE3 62 975 427
Cv + Cpy + Bo (+/- Cc, +/- Py) LE4 63 22 20
Py PR1 72 446 185
Cpy - Py PR2 70 6348 1738
Bn - Cpy (+/- Py) PR3 71 269 162
Lixiviacin Hipgena (Mg-Esp) LxH 73 0 0
Sin Mineralizacin SM 80 3 6
skarnificacin sk 90 0 0

151
Histogramas por alteracin para elementos mayores y traza

152
153
154
Boxplots por alteracin (log base 10)

155
156
Histogramas litologas 31 y 50

157
158
159
Histogramas de minerales segn alteracin

160
Mineraloga vs elementos mayores

161
162
163
164
165
166
167
Mineraloga vs abundancia relativa (tamao indica mayor ndice de abundancia)

168
169
Diagramas de dispersin de elementos mayores

170
171
172
Scatter de elementos con menor traslape entre alteraciones

173
174
175
176
177
178
179
180
181
182
Matrices de confusin entre clasificacin y mapeo de alteraciones considerando
distintos grupos de elementos mayores
Alteracin Clasificada Alteracin Clasificada
Al -Mg Al -Mg -K
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 58% 0% 32% 9% 0% 0% 0% 74% 0% 17% 9% 0% 0% 0%


K2 0% 33% 0% 43% 10% 0% 14% 10% 14% 0% 67% 5% 0% 5%
SCC1 27% 4% 44% 24% 1% 1% 0% 14% 3% 49% 33% 1% 0% 0%
SCC2 21% 14% 16% 36% 12% 0% 0% 9% 11% 19% 49% 12% 0% 0%
S1 1% 3% 0% 4% 63% 18% 10% 1% 3% 1% 7% 76% 3% 9%
S2 0% 2% 1% 4% 63% 23% 6% 0% 1% 2% 6% 77% 8% 6%
AA 0% 0% 10% 0% 30% 20% 40% 0% 0% 10% 10% 70% 0% 10%
Al -K Al -Mg -Ca
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 77% 0% 13% 9% 0% 0% 0% 66% 0% 4% 30% 0% 0% 0%


K2 24% 14% 0% 14% 33% 0% 14% 0% 33% 0% 43% 10% 0% 14%
SCC1 13% 4% 55% 24% 3% 0% 1% 28% 4% 40% 27% 1% 0% 0%
SCC2 9% 9% 25% 29% 12% 1% 16% 22% 15% 15% 40% 8% 0% 0%
S1 3% 7% 1% 12% 58% 5% 14% 1% 4% 0% 5% 78% 4% 8%
S2 2% 6% 1% 7% 65% 10% 10% 0% 2% 1% 6% 84% 3% 5%
AA 0% 0% 10% 40% 50% 0% 0% 0% 0% 10% 10% 50% 0% 30%
Al -Ca Al -Mg -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 74% 0% 15% 11% 0% 0% 0% 85% 0% 9% 6% 0% 0% 0%


K2 0% 0% 0% 29% 71% 0% 0% 10% 33% 0% 48% 5% 0% 5%
SCC1 32% 0% 25% 34% 9% 0% 1% 22% 3% 43% 31% 1% 1% 0%
SCC2 17% 7% 15% 33% 24% 2% 3% 19% 12% 18% 41% 8% 1% 0%
S1 1% 1% 1% 6% 72% 6% 14% 1% 4% 1% 10% 72% 10% 2%
S2 0% 0% 0% 6% 78% 4% 12% 0% 1% 2% 8% 74% 12% 3%
AA 0% 0% 0% 50% 50% 0% 0% 0% 0% 10% 10% 60% 10% 10%
Al -Na Al -K -Ca
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 66% 0% 19% 15% 0% 0% 0% 83% 0% 2% 15% 0% 0% 0%


K2 10% 33% 0% 52% 0% 0% 5% 29% 14% 0% 5% 52% 0% 0%
SCC1 16% 5% 42% 34% 1% 1% 0% 16% 4% 41% 34% 5% 0% 0%
SCC2 14% 11% 23% 29% 17% 2% 4% 10% 12% 21% 36% 18% 0% 2%
S1 0% 11% 1% 13% 59% 10% 6% 5% 6% 1% 8% 70% 1% 9%
S2 1% 10% 1% 8% 64% 11% 5% 2% 4% 0% 7% 79% 2% 7%
AA 0% 0% 10% 40% 40% 10% 0% 10% 0% 0% 40% 50% 0% 0%
Mg -K Al -K -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 79% 0% 13% 8% 0% 0% 0% 87% 0% 8% 6% 0% 0% 0%


K2 10% 33% 0% 33% 24% 0% 0% 29% 33% 0% 38% 0% 0% 0%
SCC1 14% 3% 39% 42% 1% 0% 0% 16% 5% 45% 32% 3% 0% 0%
SCC2 7% 15% 26% 38% 13% 0% 0% 13% 11% 25% 33% 14% 1% 3%
S1 0% 4% 1% 5% 67% 20% 3% 4% 10% 1% 18% 60% 5% 2%
S2 0% 3% 0% 5% 65% 23% 3% 2% 9% 1% 10% 65% 9% 4%
AA 0% 0% 0% 0% 80% 20% 0% 10% 0% 10% 30% 50% 0% 0%

183
Alteracin Clasificada Alteracin Clasificada
Mg -Ca Al -Ca -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 77% 0% 9% 13% 0% 0% 0% 70% 0% 15% 15% 0% 0% 0%


K2 0% 52% 0% 24% 19% 5% 0% 29% 33% 0% 33% 5% 0% 0%
SCC1 42% 5% 14% 38% 1% 0% 0% 25% 3% 40% 28% 3% 1% 0%
SCC2 21% 26% 12% 32% 9% 0% 0% 20% 11% 17% 37% 12% 3% 0%
S1 1% 5% 1% 3% 73% 17% 0% 3% 6% 1% 11% 66% 9% 3%
S2 0% 4% 0% 4% 80% 11% 0% 4% 3% 1% 6% 71% 11% 4%
AA 0% 0% 0% 0% 50% 50% 0% 10% 0% 10% 30% 50% 0% 0%
Mg -Na Mg -K -Ca
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 91% 0% 8% 2% 0% 0% 0% 81% 0% 8% 11% 0% 0% 0%


K2 19% 43% 0% 19% 10% 10% 0% 10% 29% 0% 38% 24% 0% 0%
SCC1 42% 3% 16% 37% 1% 1% 0% 17% 3% 26% 53% 1% 0% 0%
SCC2 24% 18% 12% 34% 11% 1% 0% 7% 17% 19% 47% 10% 0% 0%
S1 1% 6% 0% 5% 59% 29% 0% 0% 4% 1% 5% 85% 4% 1%
S2 0% 4% 0% 6% 60% 30% 0% 0% 3% 0% 5% 88% 3% 1%
AA 0% 0% 0% 0% 40% 60% 0% 0% 0% 0% 0% 80% 20% 0%
K -Ca Mg -K -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 83% 0% 8% 9% 0% 0% 0% 87% 0% 9% 4% 0% 0% 0%
K2 29% 38% 10% 5% 19% 0% 0% 24% 29% 0% 33% 14% 0% 0%
SCC1 20% 12% 28% 12% 24% 2% 1% 20% 3% 39% 37% 1% 0% 0%
SCC2 9% 29% 18% 18% 23% 1% 2% 9% 13% 26% 42% 10% 0% 0%
S1 7% 30% 5% 11% 38% 2% 7% 0% 4% 1% 9% 78% 7% 1%
S2 2% 31% 1% 6% 48% 2% 8% 0% 3% 1% 7% 76% 12% 1%
AA 10% 10% 10% 40% 30% 0% 0% 0% 0% 0% 0% 80% 20% 0%
K -Na Mg -Ca -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 89% 0% 11% 0% 0% 0% 0% 91% 0% 8% 2% 0% 0% 0%


K2 33% 29% 19% 19% 0% 0% 0% 19% 33% 0% 29% 10% 10% 0%
SCC1 24% 13% 26% 31% 5% 1% 0% 45% 2% 22% 29% 2% 0% 0%
SCC2 16% 26% 19% 22% 15% 1% 2% 27% 13% 11% 38% 10% 0% 0%
S1 17% 18% 2% 12% 41% 3% 6% 1% 2% 0% 7% 75% 14% 0%
S2 7% 18% 2% 8% 54% 3% 9% 1% 1% 0% 6% 78% 13% 0%
AA 20% 0% 0% 30% 50% 0% 0% 0% 0% 0% 0% 60% 40% 0%
Ca -Na K -Ca -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 81% 0% 15% 4% 0% 0% 0% 91% 0% 8% 2% 0% 0% 0%


K2 33% 10% 19% 24% 10% 5% 0% 67% 19% 5% 10% 0% 0% 0%
SCC1 29% 11% 29% 16% 14% 1% 0% 34% 7% 26% 24% 9% 1% 0%
SCC2 18% 13% 21% 25% 22% 2% 0% 19% 15% 18% 32% 18% 0% 0%
S1 3% 6% 3% 11% 75% 3% 0% 19% 6% 4% 10% 59% 0% 1%
S2 2% 3% 2% 8% 83% 2% 0% 9% 5% 3% 6% 76% 0% 1%
AA 0% 0% 10% 0% 80% 10% 0% 20% 0% 10% 0% 60% 0% 10%

184
Alteracin Clasificada Alteracin Clasificada
Al -Mg -K -Ca Al -Mg -K -Ca -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 81% 0% 4% 15% 0% 0% 0% 89% 0% 8% 4% 0% 0% 0%


K2 10% 14% 0% 62% 14% 0% 0% 19% 14% 0% 48% 19% 0% 0%
SCC1 17% 3% 42% 36% 2% 0% 0% 19% 1% 43% 34% 2% 0% 0%
SCC2 10% 11% 17% 53% 9% 0% 0% 11% 6% 18% 58% 6% 0% 0%
S1 1% 2% 1% 9% 81% 0% 5% 1% 1% 2% 15% 78% 2% 1%
S2 0% 1% 1% 8% 85% 2% 4% 1% 0% 2% 9% 82% 5% 1%
AA 0% 0% 10% 10% 70% 0% 10% 0% 0% 10% 20% 70% 0% 0%
Al -Mg -K -Na Al -K -Ca -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 87% 0% 8% 6% 0% 0% 0% 91% 0% 6% 4% 0% 0% 0%
K2 14% 14% 0% 67% 5% 0% 0% 43% 29% 0% 19% 10% 0% 0%
SCC1 15% 3% 49% 32% 1% 0% 0% 22% 4% 44% 26% 3% 0% 0%
SCC2 9% 9% 21% 54% 7% 0% 0% 17% 9% 19% 42% 12% 0% 0%
S1 0% 2% 1% 15% 76% 3% 2% 9% 5% 1% 13% 68% 3% 1%
S2 0% 1% 2% 11% 77% 8% 2% 5% 3% 1% 6% 75% 7% 3%
AA 0% 0% 10% 30% 60% 0% 0% 30% 0% 10% 10% 50% 0% 0%
Al -Mg -Ca -Na Mg -K -Ca -Na
K1 K2 SCC1 SCC2 S1 S2 AA K1 K2 SCC1 SCC2 S1 S2 AA
Alteracin mapeada

K1 85% 0% 8% 8% 0% 0% 0% 87% 0% 8% 6% 0% 0% 0%
K2 14% 29% 0% 48% 5% 0% 5% 24% 24% 0% 38% 14% 0% 0%
SCC1 25% 3% 40% 31% 1% 0% 0% 27% 3% 31% 37% 2% 0% 0%
SCC2 19% 8% 16% 50% 6% 1% 0% 15% 9% 19% 49% 9% 0% 0%
S1 1% 2% 1% 11% 75% 8% 1% 1% 2% 0% 11% 83% 4% 0%
S2 1% 1% 2% 8% 78% 9% 1% 1% 1% 1% 7% 86% 4% 0%
AA 0% 0% 10% 10% 70% 0% 10% 0% 0% 0% 0% 80% 20% 0%

185
Matriz de correlacin para elementos mayores [wt%]

Al Mg K Fe S Ca Na Cu Ti
Al 1.00 0.94 0.59 0.38 -0.23 0.32 0.70 -0.22 0.48
Mg 0.94 1.00 0.55 0.38 -0.23 0.36 0.64 -0.23 0.51
K 0.59 0.55 1.00 0.13 -0.18 0.28 0.49 -0.15 0.71
Fe 0.38 0.38 0.13 1.00 0.64 0.15 0.21 -0.15 0.15
S -0.23 -0.23 -0.18 0.64 1.00 0.25 -0.32 0.04 -0.22
Ca 0.32 0.36 0.28 0.15 0.25 1.00 0.15 -0.16 0.24
Na 0.70 0.64 0.49 0.21 -0.32 0.15 1.00 -0.16 0.57
Cu -0.22 -0.23 -0.15 -0.15 0.04 -0.16 -0.16 1.00 -0.12
Ti 0.48 0.51 0.71 0.15 -0.22 0.24 0.57 -0.12 1.00

186
Matriz de correlacin para elementos traza [ppm]

Ag Mo Re Cd Pb Zn Te Bi Sb Co Ni Se As Ba Be Ce Cr
Ag 1.00 0.02 0.12 0.18 0.37 0.37 0.32 0.68 0.46 -0.05 0.07 0.41 0.46 -0.11 -0.05 0.02 0.15
Mo 0.02 1.00 0.30 0.02 -0.01 -0.03 -0.03 0.04 0.01 0.04 0.02 0.18 0.03 -0.08 -0.01 0.07 -0.07
Re 0.12 0.30 1.00 0.15 -0.04 -0.07 0.03 0.16 0.04 0.13 -0.05 0.21 0.10 -0.18 -0.25 -0.14 0.00
Cd 0.18 0.02 0.15 1.00 0.24 0.32 0.03 0.08 0.08 0.08 0.08 0.07 0.04 -0.01 0.07 0.09 -0.01
Pb 0.37 -0.01 -0.04 0.24 1.00 0.58 0.07 0.10 0.27 0.02 0.20 0.04 0.07 0.06 0.15 0.15 0.01
Zn 0.37 -0.03 -0.07 0.32 0.58 1.00 0.04 0.16 0.19 0.04 0.19 0.05 0.07 0.09 0.37 0.36 0.00
Te 0.32 -0.03 0.03 0.03 0.07 0.04 1.00 0.42 0.24 0.02 0.01 0.28 0.32 -0.15 -0.15 -0.21 0.07
Bi 0.68 0.04 0.16 0.08 0.10 0.16 0.42 1.00 0.29 0.00 0.02 0.39 0.44 -0.16 -0.15 -0.10 0.14
Sb 0.46 0.01 0.04 0.08 0.27 0.19 0.24 0.29 1.00 -0.01 0.06 0.15 0.47 -0.07 -0.07 -0.10 0.05
Co -0.05 0.04 0.13 0.08 0.02 0.04 0.02 0.00 -0.01 1.00 0.13 0.05 0.00 -0.11 0.02 -0.12 -0.10
Ni 0.07 0.02 -0.05 0.08 0.20 0.19 0.01 0.02 0.06 0.13 1.00 -0.02 -0.02 0.05 0.10 0.04 -0.15
Se 0.41 0.18 0.21 0.07 0.04 0.05 0.28 0.39 0.15 0.05 -0.02 1.00 0.19 -0.23 -0.18 -0.12 0.18
As 0.46 0.03 0.10 0.04 0.07 0.07 0.32 0.44 0.47 0.00 -0.02 0.19 1.00 -0.11 -0.14 -0.15 0.05
Ba -0.11 -0.08 -0.18 -0.01 0.06 0.09 -0.15 -0.16 -0.07 -0.11 0.05 -0.23 -0.11 1.00 0.26 0.29 -0.03
Be -0.05 -0.01 -0.25 0.07 0.15 0.37 -0.15 -0.15 -0.07 0.02 0.10 -0.18 -0.14 0.26 1.00 0.58 -0.11
Ce 0.02 0.07 -0.14 0.09 0.15 0.36 -0.21 -0.10 -0.10 -0.12 0.04 -0.12 -0.15 0.29 0.58 1.00 -0.01
Cr 0.15 -0.07 0.00 -0.01 0.01 0.00 0.07 0.14 0.05 -0.10 -0.15 0.18 0.05 -0.03 -0.11 -0.01 1.00
Cs 0.04 -0.07 -0.16 0.09 0.19 0.31 -0.03 -0.02 -0.01 0.07 0.05 -0.23 -0.03 0.24 0.52 0.36 -0.14
Cu 0.37 0.22 0.22 0.07 -0.05 -0.04 0.12 0.28 0.09 -0.14 -0.07 0.56 0.15 -0.19 -0.21 -0.05 0.10
Ga -0.05 -0.06 -0.24 0.07 0.18 0.40 -0.14 -0.14 -0.07 0.10 0.19 -0.23 -0.14 0.30 0.77 0.52 -0.13
In 0.45 0.05 0.28 0.27 0.20 0.46 0.18 0.38 0.19 0.07 0.04 0.23 0.20 -0.10 -0.06 0.00 0.08
La 0.03 0.07 -0.11 0.08 0.14 0.33 -0.21 -0.09 -0.10 -0.14 0.02 -0.10 -0.15 0.28 0.52 0.98 0.01
Mn 0.05 -0.03 -0.19 0.07 0.33 0.51 -0.08 -0.05 -0.02 0.03 0.15 -0.23 -0.10 0.25 0.58 0.47 -0.08
P -0.06 0.01 -0.23 0.08 0.19 0.39 -0.13 -0.15 -0.06 0.08 0.17 -0.23 -0.13 0.24 0.77 0.60 -0.14
Rb -0.03 0.02 -0.08 0.02 0.13 0.23 -0.12 -0.09 -0.04 0.07 0.10 -0.15 -0.09 0.41 0.43 0.35 -0.14
Sc -0.07 -0.01 -0.16 0.05 0.11 0.27 -0.12 -0.13 -0.06 0.15 0.19 -0.18 -0.12 0.28 0.57 0.35 -0.14
Sn 0.23 0.05 0.07 -0.01 -0.02 -0.04 0.41 0.38 0.16 0.03 0.03 0.20 0.24 -0.10 -0.11 -0.19 0.02
Sr -0.01 0.01 -0.07 0.01 -0.06 -0.07 0.01 0.00 -0.02 -0.01 -0.06 0.04 0.02 0.00 0.27 -0.01 -0.04
Th -0.04 -0.03 -0.26 0.03 0.11 0.29 -0.22 -0.14 -0.10 -0.16 0.05 -0.15 -0.17 0.28 0.57 0.73 0.04
Tl 0.12 0.14 0.54 0.09 -0.01 -0.04 0.06 0.14 0.06 0.27 0.02 0.14 0.12 -0.05 -0.16 -0.17 -0.06
U 0.28 0.07 0.18 0.21 0.05 0.14 0.01 0.25 0.07 0.03 -0.02 0.11 0.09 0.04 0.07 0.23 0.06
V -0.09 -0.01 -0.18 0.05 0.12 0.30 -0.12 -0.14 -0.07 0.15 0.20 -0.22 -0.12 0.29 0.62 0.38 -0.15
W 0.07 0.01 0.00 0.00 0.01 0.00 0.03 0.05 0.05 0.01 0.02 0.03 0.06 0.00 -0.01 -0.03 0.00
Y -0.04 0.03 -0.20 0.14 0.13 0.34 -0.14 -0.13 -0.07 0.01 0.12 -0.12 -0.14 0.22 0.75 0.60 -0.11

187
Matriz de correlacin para elementos traza [ppm] (continuacin)

Cs Cu Ga In La Mn P Rb Sc Sn Sr Th Tl U V W Y
Ag 0.04 0.37 -0.05 0.45 0.03 0.05 -0.06 -0.03 -0.07 0.23 -0.01 -0.04 0.12 0.28 -0.09 0.07 -0.04
Mo -0.07 0.22 -0.06 0.05 0.07 -0.03 0.01 0.02 -0.01 0.05 0.01 -0.03 0.14 0.07 -0.01 0.01 0.03
Re -0.16 0.22 -0.24 0.28 -0.11 -0.19 -0.23 -0.08 -0.16 0.07 -0.07 -0.26 0.54 0.18 -0.18 0.00 -0.20
Cd 0.09 0.07 0.07 0.27 0.08 0.07 0.08 0.02 0.05 -0.01 0.01 0.03 0.09 0.21 0.05 0.00 0.14
Pb 0.19 -0.05 0.18 0.20 0.14 0.33 0.19 0.13 0.11 -0.02 -0.06 0.11 -0.01 0.05 0.12 0.01 0.13
Zn 0.31 -0.04 0.40 0.46 0.33 0.51 0.39 0.23 0.27 -0.04 -0.07 0.29 -0.04 0.14 0.30 0.00 0.34
Te -0.03 0.12 -0.14 0.18 -0.21 -0.08 -0.13 -0.12 -0.12 0.41 0.01 -0.22 0.06 0.01 -0.12 0.03 -0.14
Bi -0.02 0.28 -0.14 0.38 -0.09 -0.05 -0.15 -0.09 -0.13 0.38 0.00 -0.14 0.14 0.25 -0.14 0.05 -0.13
Sb -0.01 0.09 -0.07 0.19 -0.10 -0.02 -0.06 -0.04 -0.06 0.16 -0.02 -0.10 0.06 0.07 -0.07 0.05 -0.07
Co 0.07 -0.14 0.10 0.07 -0.14 0.03 0.08 0.07 0.15 0.03 -0.01 -0.16 0.27 0.03 0.15 0.01 0.01
Ni 0.05 -0.07 0.19 0.04 0.02 0.15 0.17 0.10 0.19 0.03 -0.06 0.05 0.02 -0.02 0.20 0.02 0.12
Se -0.23 0.56 -0.23 0.23 -0.10 -0.23 -0.23 -0.15 -0.18 0.20 0.04 -0.15 0.14 0.11 -0.22 0.03 -0.12
As -0.03 0.15 -0.14 0.20 -0.15 -0.10 -0.13 -0.09 -0.12 0.24 0.02 -0.17 0.12 0.09 -0.12 0.06 -0.14
Ba 0.24 -0.19 0.30 -0.10 0.28 0.25 0.24 0.41 0.28 -0.10 0.00 0.28 -0.05 0.04 0.29 0.00 0.22
Be 0.52 -0.21 0.77 -0.06 0.52 0.58 0.77 0.43 0.57 -0.11 0.27 0.57 -0.16 0.07 0.62 -0.01 0.75
Ce 0.36 -0.05 0.52 0.00 0.98 0.47 0.60 0.35 0.35 -0.19 -0.01 0.73 -0.17 0.23 0.38 -0.03 0.60
Cr -0.14 0.10 -0.13 0.08 0.01 -0.08 -0.14 -0.14 -0.14 0.02 -0.04 0.04 -0.06 0.06 -0.15 0.00 -0.11
Cs 1.00 -0.23 0.54 0.07 0.33 0.56 0.50 0.51 0.47 -0.04 0.14 0.27 0.01 0.11 0.51 0.00 0.44
Cu -0.23 1.00 -0.27 0.19 -0.03 -0.27 -0.26 -0.12 -0.18 0.08 -0.01 -0.13 0.22 0.12 -0.22 0.00 -0.10
Ga 0.54 -0.27 1.00 -0.06 0.44 0.67 0.81 0.67 0.87 -0.08 0.03 0.51 -0.05 -0.01 0.91 -0.01 0.72
In 0.07 0.19 -0.06 1.00 0.02 0.03 -0.08 -0.04 -0.09 0.05 0.07 -0.05 0.17 0.35 -0.11 0.03 -0.08
La 0.33 -0.03 0.44 0.02 1.00 0.42 0.52 0.30 0.28 -0.19 -0.02 0.71 -0.16 0.25 0.30 -0.03 0.53
Mn 0.56 -0.27 0.67 0.03 0.42 1.00 0.66 0.42 0.50 -0.07 -0.06 0.43 -0.10 0.02 0.57 -0.01 0.50
P 0.50 -0.26 0.81 -0.08 0.52 0.66 1.00 0.48 0.67 -0.08 0.12 0.51 -0.13 -0.01 0.73 0.00 0.73
Rb 0.51 -0.12 0.67 -0.04 0.30 0.42 0.48 1.00 0.75 -0.03 -0.02 0.27 0.24 0.01 0.72 0.00 0.45
Sc 0.47 -0.18 0.87 -0.09 0.28 0.50 0.67 0.75 1.00 -0.04 0.00 0.35 0.07 -0.03 0.95 0.00 0.62
Sn -0.04 0.08 -0.08 0.05 -0.19 -0.07 -0.08 -0.03 -0.04 1.00 -0.01 -0.17 0.09 0.01 -0.04 0.25 -0.09
Sr 0.14 -0.01 0.03 0.07 -0.02 -0.06 0.12 -0.02 0.00 -0.01 1.00 0.06 -0.07 0.12 -0.01 0.06 0.08
Th 0.27 -0.13 0.51 -0.05 0.71 0.43 0.51 0.27 0.35 -0.17 0.06 1.00 -0.30 0.19 0.36 -0.02 0.55
Tl 0.01 0.22 -0.05 0.17 -0.16 -0.10 -0.13 0.24 0.07 0.09 -0.07 -0.30 1.00 0.10 0.03 0.01 -0.12
U 0.11 0.12 -0.01 0.35 0.25 0.02 -0.01 0.01 -0.03 0.01 0.12 0.19 0.10 1.00 -0.06 0.02 0.10
V 0.51 -0.22 0.91 -0.11 0.30 0.57 0.73 0.72 0.95 -0.04 -0.01 0.36 0.03 -0.06 1.00 0.00 0.65
W 0.00 0.00 -0.01 0.03 -0.03 -0.01 0.00 0.00 0.00 0.25 0.06 -0.02 0.01 0.02 0.00 1.00 -0.02
Y 0.44 -0.10 0.72 -0.08 0.53 0.50 0.73 0.45 0.62 -0.09 0.08 0.55 -0.12 0.10 0.65 -0.02 1.00

188
Componentes principales elementos mayores

Valores propios

Valores propios Porcentaje % Acumulado


PC1 3.87 43.00 43.00
PC2 1.77 19.71 62.71
PC3 0.95 10.55 73.26
PC4 0.90 10.04 83.30
PC5 0.71 7.94 91.24
PC6 0.41 4.51 95.75
PC7 0.25 2.81 98.56
PC8 0.08 0.87 99.43
PC9 0.05 0.57 100.00

Vectores propios

PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 PC9


Al 0.46 0.03 0.02 -0.19 0.35 0.24 -0.04 0.16 0.74
Mg 0.46 0.04 -0.01 -0.14 0.35 0.27 0.27 0.30 -0.63
K 0.40 -0.08 0.14 0.28 -0.42 0.46 -0.56 -0.15 -0.13
Fe 0.18 0.61 0.17 -0.40 -0.13 0.00 0.16 -0.60 -0.04
S -0.12 0.69 0.14 0.08 -0.17 -0.06 -0.22 0.63 0.02
Ca 0.20 0.31 -0.20 0.74 0.39 -0.21 0.03 -0.26 0.01
Na 0.41 -0.14 0.08 -0.24 0.05 -0.75 -0.42 0.04 -0.12
Cu -0.15 -0.10 0.92 0.15 0.30 0.00 0.02 -0.05 -0.01
Ti 0.38 -0.11 0.18 0.24 -0.54 -0.22 0.60 0.16 0.15

189
Componentes principales elementos traza

Valores propios

Valores propios Porcentaje % Acumulado


PC1 8.10 23.13 23.13
PC2 4.20 11.99 35.12
PC3 2.36 6.74 41.86
PC4 2.00 5.70 47.57
PC5 1.64 4.68 52.24
PC6 1.30 3.72 55.97
PC7 1.26 3.61 59.58
PC8 1.20 3.43 63.01
PC9 1.09 3.12 66.12
PC10 0.98 2.79 68.91
PC11 0.95 2.71 71.62
PC12 0.90 2.59 74.21
PC13 0.83 2.37 76.58
PC14 0.79 2.26 78.84
PC15 0.72 2.04 80.88
PC16 0.67 1.91 82.79
PC17 0.61 1.75 84.53
PC18 0.57 1.62 86.15
PC19 0.55 1.56 87.71
PC20 0.48 1.38 89.09
PC21 0.46 1.31 90.40
PC22 0.44 1.27 91.67
PC23 0.43 1.22 92.89
PC24 0.37 1.07 93.96
PC25 0.35 1.00 94.96
PC26 0.33 0.93 95.90
PC27 0.31 0.87 96.77
PC28 0.24 0.67 97.44
PC29 0.22 0.63 98.07
PC30 0.20 0.58 98.65
PC31 0.19 0.53 99.18
PC32 0.17 0.48 99.67
PC33 0.07 0.19 99.86
PC34 0.04 0.11 99.97
PC35 0.01 0.03 100.00

190
Vectores propios

PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 PC9 PC10 PC11 PC12
Ag -0.05 -0.40 0.08 0.07 0.03 -0.08 0.05 0.10 0.06 0.03 -0.09 -0.02
Mo -0.02 -0.06 0.01 -0.33 0.12 -0.12 -0.10 -0.44 -0.06 -0.09 -0.02 0.34
Re -0.11 -0.12 -0.09 -0.46 -0.08 0.03 -0.19 -0.02 -0.09 -0.12 0.12 0.15
Cd 0.03 -0.17 0.04 -0.14 -0.31 0.15 0.10 -0.03 -0.03 0.22 0.02 -0.04
Pb 0.08 -0.22 0.02 0.15 -0.39 -0.16 0.07 -0.08 -0.02 0.19 -0.02 0.22
Zn 0.16 -0.27 0.08 0.07 -0.33 -0.07 0.16 0.01 -0.03 0.07 0.04 0.17
Te -0.09 -0.21 -0.11 0.19 0.21 0.02 0.23 0.12 -0.07 -0.27 -0.12 0.12
Bi -0.09 -0.34 0.01 0.05 0.18 0.02 0.11 0.18 -0.04 -0.20 -0.11 -0.05
Sb -0.05 -0.27 -0.06 0.25 -0.04 -0.08 -0.33 -0.18 0.18 0.01 0.08 -0.10
Hg -0.05 -0.27 -0.07 0.18 -0.02 0.01 -0.40 -0.25 0.07 0.06 0.18 -0.08
Co 0.01 -0.03 -0.30 -0.17 -0.16 0.18 0.17 -0.10 -0.06 -0.25 0.33 -0.37
Ni 0.06 -0.07 -0.15 0.04 -0.21 -0.16 0.17 -0.30 -0.21 0.02 -0.40 -0.51
Se -0.12 -0.22 0.08 -0.17 0.24 -0.18 0.31 -0.09 0.12 0.18 0.03 -0.08
As -0.09 -0.27 -0.08 0.19 0.13 0.04 -0.31 -0.06 0.16 -0.17 0.03 -0.09
Ba 0.15 0.05 0.02 0.05 0.01 -0.07 -0.33 0.32 -0.10 0.25 -0.34 -0.09
Be 0.29 -0.04 0.04 0.02 0.11 0.18 0.08 -0.14 0.11 0.00 0.05 0.03
Ce 0.25 -0.05 0.34 -0.11 0.07 -0.07 -0.12 -0.04 -0.13 -0.17 0.02 -0.04
Cr -0.05 -0.05 0.19 0.06 0.09 -0.15 0.10 0.35 0.05 0.17 0.61 -0.17
Cs 0.22 -0.08 -0.10 0.04 -0.02 0.24 -0.05 0.19 0.06 -0.07 -0.09 0.25
Cu -0.11 -0.17 0.14 -0.27 0.25 -0.25 0.12 -0.08 0.19 0.28 -0.15 0.00
Ga 0.32 -0.05 -0.15 0.01 0.07 -0.04 0.06 0.01 0.06 0.02 0.07 -0.04
In -0.03 -0.30 0.08 -0.13 -0.23 0.22 0.10 0.18 0.00 0.04 -0.04 0.07
La 0.23 -0.05 0.37 -0.12 0.06 -0.07 -0.14 -0.03 -0.15 -0.17 0.00 -0.04
Mn 0.26 -0.08 -0.03 0.09 -0.11 -0.05 0.07 0.06 -0.08 -0.14 0.06 0.26
P 0.30 -0.04 -0.04 0.02 0.06 0.04 0.08 -0.15 0.03 -0.08 0.10 0.04
Rb 0.23 -0.05 -0.24 -0.13 0.10 -0.12 -0.16 0.24 0.05 0.17 -0.13 0.02
Sc 0.28 -0.03 -0.27 -0.07 0.12 -0.10 0.03 0.06 0.09 0.11 0.04 -0.11
Sn -0.07 -0.16 -0.17 0.13 0.33 0.05 0.11 0.03 -0.45 -0.11 -0.05 0.12
Sr 0.02 -0.01 0.05 0.03 0.19 0.62 0.08 -0.18 0.31 0.28 -0.12 0.04
Th 0.24 0.00 0.32 0.02 0.09 -0.02 -0.06 -0.04 -0.08 -0.06 0.03 -0.15
Tl -0.05 -0.12 -0.29 -0.42 -0.03 -0.02 -0.20 0.17 0.00 0.01 0.03 0.01
U 0.02 -0.18 0.22 -0.22 -0.02 0.36 -0.11 0.18 -0.11 -0.10 -0.12 -0.32
V 0.29 -0.03 -0.27 -0.04 0.11 -0.10 0.05 0.04 0.07 0.07 0.06 -0.07
W -0.01 -0.06 -0.06 0.08 0.16 0.18 -0.09 -0.13 -0.63 0.49 0.21 0.05
Y 0.28 -0.05 0.04 -0.05 0.14 0.01 0.08 -0.14 0.06 -0.01 0.05 -0.06

191
Valores propios elementos traza (Continuacin)
PC13 PC14 PC15 PC16 PC17 PC18 PC19 PC20 PC21 PC22 PC23 PC24
Ag 0.10 0.16 -0.01 0.05 -0.19 0.09 0.24 -0.14 0.02 0.24 0.09 0.00
Mo 0.24 -0.41 -0.18 0.43 -0.02 0.14 -0.07 -0.05 0.12 0.10 -0.03 0.09
Re -0.01 -0.08 -0.17 -0.36 0.07 -0.12 0.15 -0.13 -0.25 -0.06 0.21 -0.29
Cd -0.72 -0.31 0.04 0.07 -0.04 0.29 0.05 0.04 0.06 0.09 0.11 -0.07
Pb 0.18 -0.13 0.26 -0.09 -0.38 -0.18 0.11 -0.22 0.19 -0.23 -0.04 0.03
Zn 0.06 0.07 0.00 0.04 0.24 -0.10 0.05 0.06 0.05 -0.06 -0.21 0.06
Te -0.13 -0.28 0.20 -0.16 0.19 -0.04 -0.51 -0.43 -0.02 -0.01 0.03 0.13
Bi 0.05 0.09 -0.14 0.08 -0.05 0.15 0.34 -0.04 0.07 0.10 0.47 0.23
Sb -0.03 -0.05 0.07 0.00 -0.13 -0.27 -0.20 0.10 -0.17 0.62 -0.03 -0.21
Hg -0.05 -0.02 0.00 -0.03 0.15 -0.08 -0.12 0.26 -0.05 -0.46 0.34 0.36
Co 0.20 0.00 0.50 0.25 0.08 0.09 0.08 0.08 0.06 0.09 -0.01 0.08
Ni 0.22 -0.13 -0.29 -0.30 -0.02 0.22 -0.15 0.10 -0.01 0.00 -0.01 -0.05
Se 0.09 0.00 0.29 0.06 0.10 -0.05 -0.04 0.08 -0.24 -0.25 0.17 -0.50
As -0.08 0.08 -0.05 0.05 0.14 0.40 0.17 -0.22 0.15 -0.26 -0.48 -0.30
Ba 0.16 -0.39 0.27 0.14 0.32 0.08 0.23 -0.04 -0.31 0.05 0.00 0.08
Be 0.02 -0.03 0.02 -0.13 0.03 -0.02 0.12 -0.01 -0.23 0.04 -0.09 0.10
Ce -0.02 0.05 0.17 -0.14 0.05 0.07 -0.06 0.03 0.20 0.07 0.06 -0.04
Cr 0.21 -0.40 -0.24 -0.14 -0.08 0.18 -0.09 0.07 0.04 0.03 -0.09 0.06
Cs 0.10 0.09 0.07 0.01 -0.26 0.39 -0.35 0.38 -0.09 -0.05 0.10 -0.13
Cu -0.09 0.21 0.10 -0.01 -0.04 0.11 -0.11 0.16 -0.10 0.01 -0.26 0.33
Ga -0.04 0.03 -0.12 0.02 0.06 -0.09 0.04 -0.08 -0.02 -0.02 0.04 -0.02
In 0.12 0.13 -0.26 0.11 0.51 -0.18 -0.16 0.15 0.11 0.13 -0.11 0.00
La 0.00 0.05 0.19 -0.14 0.04 0.07 -0.08 0.04 0.23 0.07 0.08 -0.06
Mn 0.15 0.04 -0.07 -0.01 -0.05 0.18 0.03 0.06 -0.39 -0.07 -0.05 -0.07
P -0.02 -0.01 -0.03 -0.09 0.04 0.03 0.12 -0.09 -0.12 0.08 -0.04 0.07
Rb 0.04 -0.01 0.05 0.07 -0.05 -0.06 -0.15 0.07 0.35 -0.02 0.08 -0.13
Sc -0.08 0.02 -0.15 0.13 -0.01 -0.14 -0.03 -0.10 0.15 -0.04 0.06 -0.07
Sn -0.14 -0.22 0.02 -0.08 -0.10 -0.27 0.26 0.49 0.11 -0.04 -0.26 -0.07
Sr 0.23 -0.19 0.04 -0.25 0.01 -0.02 0.10 -0.02 0.22 0.04 0.05 -0.04
Th -0.01 -0.01 0.00 -0.06 0.07 -0.16 0.03 0.00 0.12 -0.08 0.02 -0.08
Tl 0.05 0.07 0.15 -0.41 -0.09 0.03 0.02 -0.05 0.05 0.07 -0.21 0.25
U 0.02 -0.04 -0.13 0.30 -0.40 -0.26 -0.17 -0.15 -0.21 -0.23 -0.17 0.03
V -0.07 0.01 -0.15 0.11 0.00 -0.08 -0.01 -0.09 0.06 -0.02 0.05 -0.04
W 0.05 0.30 0.01 0.04 0.02 0.15 -0.13 -0.25 -0.05 0.10 0.04 -0.02
Y -0.19 -0.01 -0.02 -0.02 -0.04 -0.05 0.04 -0.05 -0.23 0.01 -0.13 0.24

192
Valores propios elementos traza (Continuacin)
PC25 PC26 PC27 PC28 PC29 PC30 PC31 PC32 PC33 PC34 PC35
Ag 0.03 -0.11 0.05 -0.23 0.25 -0.46 -0.47 -0.03 -0.02 -0.01 0.00
Mo -0.02 0.13 0.05 -0.04 0.07 0.01 -0.04 0.01 0.03 0.01 0.00
Re 0.15 -0.35 0.17 0.24 -0.02 -0.11 0.04 -0.05 0.02 0.01 0.00
Cd -0.15 0.05 0.00 -0.07 0.01 -0.03 0.04 0.02 -0.01 0.00 0.00
Pb 0.17 -0.15 0.12 -0.19 -0.12 0.18 0.21 0.06 0.00 -0.01 0.00
Zn -0.01 0.19 -0.13 0.66 0.18 0.02 -0.20 -0.10 0.03 0.01 0.00
Te -0.06 -0.07 0.06 0.04 0.04 -0.08 -0.02 -0.01 -0.01 0.00 0.00
Bi 0.00 0.16 0.03 0.20 -0.17 0.29 0.30 0.01 0.00 0.00 0.00
Sb -0.06 0.04 0.01 0.10 -0.08 0.14 0.10 0.02 0.01 0.00 0.00
Hg -0.05 -0.02 -0.08 -0.09 0.06 -0.12 -0.09 -0.01 -0.02 -0.01 0.00
Co -0.03 -0.19 0.13 0.06 -0.06 -0.09 0.03 0.00 0.01 0.01 0.00
Ni 0.04 0.01 0.00 0.02 -0.01 -0.03 0.02 0.01 0.00 0.00 0.00
Se 0.04 0.30 -0.13 -0.12 0.00 0.11 -0.03 -0.01 0.00 -0.01 0.00
As 0.06 -0.01 0.04 0.02 -0.07 0.05 0.10 0.03 0.01 0.00 0.00
Ba 0.07 -0.10 0.03 -0.01 0.03 0.10 -0.05 -0.01 0.01 0.01 0.00
Be 0.15 0.18 0.08 -0.02 0.44 -0.27 0.48 0.31 -0.24 -0.05 -0.01
Ce 0.05 -0.10 -0.23 -0.02 -0.01 0.06 -0.03 0.13 0.03 0.01 -0.73
Cr 0.06 -0.03 -0.02 0.03 -0.01 -0.02 0.02 -0.01 0.00 0.01 0.00
Cs 0.26 -0.04 0.24 0.04 0.16 0.18 -0.10 -0.10 0.03 0.03 0.00
Cu -0.16 -0.42 0.12 0.13 -0.02 0.04 0.19 -0.05 0.03 0.02 0.01
Ga -0.09 -0.05 0.01 -0.11 0.17 0.06 0.06 0.13 0.81 0.28 0.05
In 0.13 -0.06 0.06 -0.42 -0.14 0.06 0.12 0.06 -0.05 -0.02 0.00
La 0.04 -0.11 -0.24 0.00 0.00 0.06 -0.03 0.18 -0.03 -0.03 0.68
Mn -0.59 -0.04 -0.10 -0.08 -0.36 -0.13 0.00 0.20 -0.11 0.01 0.00
P 0.02 -0.09 -0.31 -0.22 0.05 0.12 0.10 -0.77 -0.10 0.05 0.04
Rb -0.01 0.14 -0.07 0.13 -0.25 -0.53 0.30 -0.20 0.03 -0.04 0.00
Sc -0.03 -0.15 0.01 0.04 0.06 0.21 -0.17 0.18 -0.49 0.54 0.00
Sn 0.01 -0.09 -0.02 -0.04 0.03 -0.01 -0.04 0.03 0.01 0.01 0.00
Sr -0.19 -0.10 -0.08 0.12 -0.16 0.03 -0.15 0.05 0.04 -0.01 0.00
Th -0.33 0.17 0.72 -0.03 0.02 0.05 -0.05 -0.25 -0.05 -0.04 0.00
Tl -0.17 0.45 -0.03 -0.18 0.08 0.20 -0.14 0.04 -0.01 -0.01 0.00
U -0.07 0.01 -0.18 0.06 0.06 0.01 0.04 -0.05 0.02 -0.01 -0.01
V -0.05 -0.18 -0.02 0.01 0.10 0.22 -0.13 0.12 -0.04 -0.78 -0.01
W 0.02 0.03 0.03 0.04 -0.02 0.02 0.04 0.00 0.00 0.00 0.00
Y 0.48 0.18 0.13 0.04 -0.56 -0.09 -0.27 0.08 0.05 0.00 0.03

193
Anexo B

Clasificacin usando grupos de elementos con menor traslape

Boxplots muestras con una alteracin


Al-Mg-K-Ca-Na

Al-Mg-K-Ca

194
Al-Mg-K-Na

Boxplots muestras con dos alteraciones

Al-Mg-K-Ca-Na

195
Al-Mg-K-Ca

Al-Mg-K-Na

196
Plantas de muestras con una alteracin

Al-Mg-K-Ca-Na

197
Al-Mg-K-Ca

198
Al-Mg-K-Na

199
Porcentaje de clasificacin versus rangos de porcentaje de alteracin por tramo

Al-Mg-K-Ca-Na

Al-Mg-K-Ca

Al-Mg-K-Na

200
Porcentaje de clasificacin versus rangos de grados de pertenencia

Al-Mg-K-Ca-Na

Al-Mg-K-Ca

Al-Mg-K-Na

201
Anexo C
Tabla de Estadgrafos para elementos en ppm
cantidad 1 2 4 Desviacin Rango
Media 3cuartil Mnimo Mximo Moda Rango
blancos cuartil cuartil cuartil Estndar intercuartil
Au ppm 4983 0.21 0.2 0.2 0.2 4.5 0.2 4.5 0.2 0.08 4.3 0
Ag ppm 0 3.09 1.27 2.24 3.83 83.1 0.01 83.1 1 3.45 83.09 2.56
Mo ppm 0 91.68 35.3 64.6 114.26 2120 0.1 2120 107.5 109.96 2119.9 78.96
Re ppm 0 0.59 0.14 0.3 0.61 13.95 0 13.95 0.1 0.96 13.95 0.46
Cd ppm 0 4.22 0.34 1.13 3.36 601 0.01 601 0.01 13.9 600.99 3.02
Pb ppm 0 70.33 4.83 20 64.4 10000 0.2 10000 1.2 226.88 9999.8 59.58
Zn ppm 0 555.96 30 242 763.75 10000 2 10000 10 912.75 9998 733.75
Te ppm 0 0.59 0.18 0.32 0.58 52.6 0.01 52.6 0.13 1.29 52.59 0.4
Bi ppm 0 1.22 0.31 0.7 1.46 59.1 0.01 59.1 0.24 1.9 59.09 1.15
Sb ppm 0 1.79 0.13 0.25 0.77 274 0.05 274 0.05 8.44 273.95 0.64
Hg ppm 0 0.03 0.01 0.01 0.03 6.87 0.01 6.87 0.01 0.09 6.86 0.02
Co ppm 0 13.4 7.4 12.1 16.9 349 0.1 349 7.4 10.24 348.9 9.5
Ni ppm 0 21.03 4.9 8.3 18.4 502 0.2 502 4.4 40.04 501.8 13.5
Se ppm 0 3.44 2.4 3.2 4.1 19.4 0.2 19.4 2.6 1.55 19.2 1.7
Al wt% 0 0.97 0.45 0.68 1.42 4.7 0.06 4.7 0.39 0.7 4.64 0.97
As ppm 0 37.68 2.7 7.8 25.88 6640 0.1 6640 1.3 142.97 6639.9 23.18
B ppm 0 10.08 10 10 10 30 10 30 10 0.91 20 0
Ba ppm 0 40.04 30 38 50 270 0.2 270 30 21.94 269.8 20
Be ppm 0 0.27 0.09 0.21 0.4 1.57 0.05 1.57 0.05 0.21 1.52 0.31
Ca wt% 0 0.27 0.03 0.09 0.28 7.41 0.01 7.41 0.02 0.58 7.4 0.25
Ce ppm 0 17.41 10.85 16.3 23.3 46.3 0.02 46.3 21.8 8.5 46.28 12.45
Cr ppm 0 103.65 38 74 138 1560 1 1560 20 97.18 1559 100
Cs ppm 0 1.13 0.67 0.95 1.41 8.81 0.05 8.81 0.65 0.66 8.76 0.74
Cu ppm 0 5999.61 4010 5770 8140 10000 66.9 10000 10000 2624.29 9933.1 4130
Ga ppm 0 2.63 0.94 1.69 3.99 13.35 0.05 13.35 0.7 2.15 13.3 3.05
Ge ppm 0 0.06 0.05 0.05 0.07 1.16 0.05 1.16 0.05 0.03 1.11 0.02

202
cantidad 1 2 4 Desviacin Rango
Media 3cuartil Mnimo Mximo Moda Rango
blancos cuartil cuartil cuartil Estndar intercuartil
Hf ppm 0 0.02 0.02 0.02 0.02 0.19 0.02 0.19 0.02 0.01 0.17 0
In ppm 0 0.28 0.08 0.15 0.32 9.25 0.01 9.25 0.05 0.42 9.25 0.25
K wt% 0 0.23 0.17 0.21 0.26 1.18 0.01 1.18 0.19 0.1 1.17 0.09
La ppm 0 8.52 5.4 8 11.3 26.7 0.2 26.7 6.8 4.1 26.5 5.9
Li ppm 2058 7.27 1 4 13 39 1 39 1 7.51 38 12
Mg wt% 0 0.47 0.03 0.08 0.89 3.13 0.01 3.13 0.02 0.6 3.12 0.86
Mn ppm 0 294.7 21 55 425 6540 5 6540 15 452.22 6535 404
Na wt% 0 0.06 0.04 0.06 0.08 0.5 0.01 0.5 0.05 0.04 0.49 0.04
Nb ppm 0 0.06 0.05 0.05 0.05 0.53 0.05 0.53 0.05 0.03 0.48 0
P ppm 0 473.36 90 240 850 3060 10 3060 50 475.73 3050 760
Rb ppm 0 12.6 8.9 10.9 14.1 76.3 0.1 76.3 10.2 6.55 76.2 5.2
S wt% 0 2.22 1.23 1.85 2.7 10 0.14 10 10 1.51 9.86 1.47
Sc ppm 0 1.73 0.5 0.9 2.3 16.5 0.1 16.5 0.4 1.96 16.4 1.8
Sn ppm 0 0.72 0.3 0.5 0.8 73.8 0.2 73.8 0.4 1.29 73.6 0.5
Sr ppm 0 93.2 28 55.45 114 1240 0.2 1240 16.6 107.13 1239.8 86
Ta ppm 0 0.01 0.01 0.01 0.01 0.26 0.01 0.26 0.01 0 0.25 0
Th ppm 0 1.85 1 1.8 2.5 9.2 0.2 9.2 0.8 0.96 9 1.5
Ti wt% 0 0.01 0.01 0.01 0.01 0.32 0.01 0.32 0.01 0.02 0.31 0.01
Tl ppm 0 0.18 0.11 0.15 0.21 2.06 0.02 2.06 0.12 0.11 2.04 0.1
U ppm 0 0.43 0.27 0.36 0.5 5.56 0.05 5.56 0.3 0.29 5.51 0.23
V ppm 0 18.23 3 6 26 196 1 196 2 24.21 195 23
W ppm 0 0.96 0.25 0.44 0.74 540 0.05 540 0.25 8.57 539.95 0.49
Y ppm 0 4.94 0.6 2.43 8.8 74.2 0.05 74.2 0.5 5.54 74.15 8.2
Zr ppm 0 1 1 1 1 4 1 4 1 0.08 3 0
CuT wt% 0 0.71 0.4 0.58 0.82 6.48 0.01 6.48 0.4 0.55 6.47 0.42
CuS wt% 123 0.04 0.01 0.02 0.04 0.71 0.01 0.71 0.01 0.04 0.7 0.03
Fe wt% 0 2.53 1.63 2.27 3.13 16.55 0.1 16.55 1.71 1.34 16.45 1.5
As ppm 1 43.46 20 20.75 33.5 2938.71 -99 2938.71 20 103.81 3037.71 13.5

203
Tabla de Estadgrafos para mineraloga QemScan
cantidad 1 2 4 Desviacin Rango
Mineral Media 3cuartil Mnimo Mximo Rango
blancos cuartil cuartil cuartil Estndar intercuartil
Chalcocite 0 0.29 0.07 0.15 0.32 1.59 0.01 1.59 0.35 1.58 0.25
Covellite 0 0.12 0.07 0.10 0.14 0.42 0.00 0.42 0.09 0.42 0.07
Chalcopyrite 0 1.12 0.57 1.07 1.62 4.08 0.01 4.08 0.77 4.07 1.05
Bornite 0 0.06 0.00 0.01 0.06 0.77 0.00 0.77 0.13 0.77 0.05
Enargite/Tennantite 16 0.03 0.00 0.00 0.01 0.46 0.00 0.46 0.08 0.46 0.01
Native Copper 76 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Cuprite 70 0.00 0.00 0.00 0.00 0.01 0.00 0.01 0.00 0.01 0.00
Malachite/Azurite 31 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Brochantite 1 0.01 0.00 0.00 0.01 0.03 0.00 0.03 0.01 0.03 0.00
Atacamite 11 0.00 0.00 0.00 0.00 0.01 0.00 0.01 0.00 0.01 0.00
Turquoise 72 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Psuedomalachite 15 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Cu-Mn Wad 54 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Chrysocolla 0 0.00 0.00 0.00 0.01 0.02 0.00 0.02 0.00 0.01 0.00
Cu-bearing Clays 3 0.00 0.00 0.00 0.00 0.01 0.00 0.01 0.00 0.01 0.00
Cu-bearing Fe Ox/OH 53 0.00 0.00 0.00 0.00 0.01 0.00 0.01 0.00 0.01 0.00
Other Cu Mins 22 0.00 0.00 0.00 0.00 0.01 0.00 0.01 0.00 0.01 0.00
Pyrite 0 2.39 0.53 1.66 3.34 15.58 0.01 15.58 2.62 15.57 2.82
Galena 18 0.00 0.00 0.00 0.00 0.11 0.00 0.11 0.01 0.11 0.00
Sphalerite 7 0.06 0.00 0.03 0.08 0.52 0.00 0.52 0.08 0.52 0.07
Molybdenite 0 0.02 0.01 0.01 0.03 0.09 0.00 0.09 0.02 0.09 0.02
Fe Oxides 0 0.32 0.05 0.08 0.27 1.75 0.02 1.75 0.49 1.73 0.22
Ilmenite 1 0.04 0.00 0.01 0.04 0.64 0.00 0.64 0.08 0.64 0.04
Rutile 0 0.41 0.35 0.40 0.45 0.73 0.25 0.73 0.10 0.48 0.10
Corundum 0 0.04 0.00 0.01 0.02 1.24 0.00 1.24 0.14 1.24 0.02
Quartz 0 35.60 26.81 37.22 43.23 50.15 17.32 50.15 9.40 32.83 16.42
Orthoclase 0 5.33 3.17 4.54 6.35 13.52 1.86 13.52 2.81 11.66 3.18

204
cantidad 1 2 4 Desviacin Rango
Mineral Media 3cuartil Mnimo Mximo Rango
blancos cuartil cuartil cuartil Estndar intercuartil
Albite 0 10.75 0.14 5.75 21.22 36.59 0.04 36.59 11.72 36.55 21.09
Plagioclase 0 2.52 0.04 0.09 1.81 20.15 0.01 20.15 4.80 20.14 1.77
Muscovite/Sericite 0 27.70 17.71 30.80 36.42 45.68 3.25 45.68 11.08 42.44 18.71
Kaolinite 0 2.62 1.17 1.85 3.01 9.83 0.57 9.83 2.22 9.26 1.84
Pyrophyllite/Smectite 0 3.79 1.68 3.38 4.50 18.35 0.62 18.35 3.07 17.74 2.82
Biotite 0 2.22 0.08 0.24 3.87 12.11 0.01 12.11 3.01 12.10 3.78
Chlorite 0 3.16 0.03 0.18 7.38 11.05 0.01 11.05 3.84 11.04 7.34
Hornblende 0 0.03 0.00 0.01 0.02 0.22 0.00 0.22 0.05 0.22 0.02
Epidote 76 0.05 0.00 0.00 0.00 3.65 0.00 3.65 0.42 3.65 0.00
Calcite/CO3 0 0.01 0.01 0.01 0.01 0.03 0.00 0.03 0.01 0.03 0.01
Gypsum 0 0.01 0.01 0.01 0.02 0.05 0.00 0.05 0.01 0.05 0.01
Alunite 0 0.51 0.05 0.08 0.49 5.77 0.02 5.77 1.15 5.75 0.43
Jarosite 0 0.01 0.00 0.00 0.01 0.06 0.00 0.06 0.01 0.06 0.00
Mg-SO4 38 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Barite 6 0.00 0.00 0.00 0.00 0.03 0.00 0.03 0.00 0.03 0.00
Svanbergite 0 0.38 0.03 0.21 0.54 1.74 0.00 1.74 0.45 1.74 0.52
Fe-Al PO4 1 0.01 0.00 0.00 0.01 0.25 0.00 0.25 0.03 0.25 0.00
Apatite 0 0.21 0.00 0.01 0.52 0.84 0.00 0.84 0.29 0.84 0.52
Monazite 0 0.05 0.04 0.05 0.05 0.08 0.02 0.08 0.01 0.06 0.02
Zircon 1 0.01 0.00 0.00 0.02 0.16 0.00 0.16 0.02 0.16 0.02
Sphene 34 0.02 0.00 0.00 0.01 0.45 0.00 0.45 0.07 0.45 0.01
Others 0 0.08 0.05 0.06 0.08 0.24 0.02 0.24 0.05 0.22 0.03

205

También podría gustarte