Modelamiento-geologico-y-estimacion-de-recursos-yacimiento-La-Niña-Tongoy - Cementos Melón
Modelamiento-geologico-y-estimacion-de-recursos-yacimiento-La-Niña-Tongoy - Cementos Melón
Modelamiento-geologico-y-estimacion-de-recursos-yacimiento-La-Niña-Tongoy - Cementos Melón
PROFESOR GUÍA:
BRIAN TOWNLEY CALLEJAS
MIEMBROS DE LA COMISIÓN:
VICTOR MAKSAEV JURCHUC
KATJA DECKART
SANTIAGO DE CHILE
2017
RESUMEN DE LA MEMORIA PARA
OPTAR AL TITULO DE: Geólogo
POR: Oscar Naranjo Cabello
FECHA: 20/10/2017
PROFESOR GUÍA: Brian Townley
i
Agradecimientos
ii
Tabla de Contenido
1. Introducción. 1
1.1 Fundamento de estudio 1
1.2 Objetivos. 3
1.2.1 Objetivos principales. 3
1.2.3 Objetivos específicos. 3
1.3 Área de estudio 4
1.4 Metodología 5
1.5 Alcance del trabajo 6
2. Marco teórico 7
2.1 Química y fabricación del cemento 7
2.1.1 Materias primas 8
2.1.2 Fabricación del clinker 8
2.1.3 Relación entre componentes del clinker y propiedades del cemento 10
2.1.4 Componentes minoritarios del clinker 12
2.1.5 Fraguado y endurecimiento 14
2.2 Problemática del Cl 15
2.3 Modelamiento Geológico 16
3. Marco geológico 17
3.1 Marco regional 17
3.1.1 Formaciones geológicas 21
3.2 Geología local 23
4. Marco Metodológico 29
4.1 Campaña sondajes Tongoy 2017 29
4.1.1 Objetivos campaña 29
4.1.2 Armado malla de sondajes 30
4.1.3 Recuperación de testigos 31
4.1.4 Descripción, armado de muestra y análisis químico. 33
4.2 Modelamiento geológico 34
iii
4.2.1 Base de datos 34
4.2.2 Modelamiento geológico 36
4.3 Estimación de recursos 42
4.3.1 Compositación muestras 43
4.3.2 Estadística descriptiva 43
4.3.3 Análisis variográfico 54
4.3.4 Criterios de estimación 61
4.3.5 Modelo de bloques 62
5. Resultados 63
5.1 Estratigrafía Local 63
5.2 Estimación del modelo de bloques. 65
5.3 Validaciones estimación 71
5.4 Calculo de recursos 74
6. Discusiones 76
6.1 Comparación Modelo Geológico SIGA (2005) 76
6.2 Comparación estimación SIGA (2005) 79
7. Conclusiones 81
8. Bibliografía 83
iv
Tabla de Cuadro
v
Tabla de Figura
FIGURA 1.1 UBICACIÓN MINA LA NIÑA ............................................................... 4
FIGURA 2.1 CURVAS DE RESISTENCIA A LA COMPRESIÓN FRENTE AL
TIEMPO DE CURADO. ................................................................................... 11
FIGURA 2.2 ELEVACIÓN DE TEMPERATURA EN MUESTRAS DE HORMIGÓN.
........................................................................................................................ 11
FIGURA 3.1 MAPA GEOLÓGICO ÁREA TONGOY .............................................. 19
FIGURA 3.2 PERFIL ESQUEMÁTICO TRANSVERSAL A LA CUENCA
CENOZOICA ................................................................................................... 20
FIGURA 3.3 GEOLOGÍA LOCAL CAMPAÑA 2001 (S1-S21), SOTO (2005), VISTA
EN PLANTA. ................................................................................................... 25
FIGURA 3.4 SECCIÓN E-O ACUÑAMIENTO LATERAL MANTOS DE
CONCHUELA, SOTO (2005) .......................................................................... 25
FIGURA 3.5 SECCIÓN N-S, ACUÑAMIENTO LATERAL DEL MANTO AL SUR,
SOTO 2005) .................................................................................................... 26
FIGURA 3.6 PALEO AMBIENTE BAHÍA DEL MIOCENO, SOTO (2005) ............. 26
FIGURA 3.7 DISTRIBUCIÓN DE LAS LEYES DE CACO3, SOTO (2005), VISTA EN
PLANTA. ......................................................................................................... 27
FIGURA 3.8 CURVA TONELAJE-LEY PARA CONCHUELA ALTA LEY(SOTO205)
........................................................................................................................ 28
FIGURA 4.1 MALLA SONDAJES CAMPAÑA 2017, VISTA EN PLANTA. ............ 30
FIGURA 4.2 MONTAJE SONDAJE SC-02 ............................................................ 31
FIGURA 4.3 DESCRIPCIÓN Y COMPOSITACIÓN DE MUESTRAS CON CORTE
DIAMETRAL .................................................................................................... 33
FIGURA 4.4 SONDAJES DESPLEGADOS POR LITOLOGÍA, 3D. ...................... 36
FIGURA 4.5 SECCIONES NORTE (AZUL) Y SECCIONES E (ROJO), VISTA EN
PLANTA. ......................................................................................................... 37
FIGURA 4.6 PERFIL LITOLÓGICO AA', PARA CONA Y CONB. ......................... 38
FIGURA 4.7 PERFILES NORTE, CONA Y CONB, 3D. ........................................ 38
FIGURA 4.8 PERFIL LITOLÓGICO BB', PARA CONA Y CONB. ......................... 39
FIGURA 4.9 PERFILES ESTE, CONA Y CONB, 3D. ............................................ 39
FIGURA 4.10 LIMITES MODELO GEOLÓGICO, TOPOGRAFÍA ACTUAL,
BATOLITO, PROPIEDAD MINERA Y AGRÍCOLA.3D. ................................... 40
FIGURA 4.11 MODELO GEOLÓGICO CONA, 3D................................................ 41
FIGURA 4.12 MODELO GEOLÓGICO CONCHUELA BAJA. 3D.......................... 41
FIGURA 4.13 MODELO GEOLÓGICO YACIMIENTO LA NIÑA, 3D. .................... 42
FIGURA 4.14 BOXPLOT CAO .............................................................................. 43
FIGURA 4.15 LOGNORMAL CAO, CONA ............................................................ 44
FIGURA 4.16 LOGNORMAL CAO, CONB ............................................................ 44
FIGURA 4.17 BOXPLOT CL ................................................................................. 45
FIGURA 4.18 LOGNORMAL CL, CONA ............................................................... 45
vi
FIGURA 4.19 LOGNORMAL CL, CONB ............................................................... 46
FIGURA 4.20 CORRELACIÓN CAO-SIO2 EN CONA Y CONB............................ 47
FIGURA 4.21 CORRELACIÓN CAO-FE2O3 EN CONA Y CONB ........................ 48
FIGURA 4.22 CORRELACIÓN CAO-AL2O3 PARA CONA Y CONB .................... 49
FIGURA 4.23 CORRELACIÓN CAO-MGO PARA CONA Y CONB ...................... 50
FIGURA 4.24 CORRELACIÓN CAO-CL PARA CONA Y CONB .......................... 51
FIGURA 4.25 CORRELACIÓN CL-SO3 PARA CONA Y CONB ........................... 52
FIGURA 4.26 CARACTERIZACIÓN VARIOGRAMAS .......................................... 55
FIGURA 4.27 VARIOGRAMAS CAO PARA CONA .............................................. 55
FIGURA 4.28 VARIOGRAMA CAO PARA CONB ................................................. 56
FIGURA 4.29 VARIOGRAMAS SIO2 PARA CONA .............................................. 56
FIGURA 4.30 VARIOGRAMAS SIO2 PARA CONB .............................................. 57
FIGURA 4.31 VARIOGRAMAS CL PARA CONA.................................................. 57
FIGURA 4.32 VARIOGRAMAS CL PARA CONB.................................................. 58
FIGURA 4.33 MODELO DE BLOQUES, CONA (A) Y CONB (B ........................... 62
FIGURA 5.1 ESTRATIGRAFÍA YACIMIENTO LA NIÑA ....................................... 64
FIGURA 5.2 DISTRIBUCIÓN ESPACIAL DEL CACO3 EN CONA, VISTA EN
PLANTA. ......................................................................................................... 66
FIGURA 5.3 ESTIMACIÓN MODELO DE BLOQUES PARA CACO3 EN CONB,
VISTA EN PLANTA ......................................................................................... 66
FIGURA 5.4 ESTIMACIÓN MODELO DE BLOQUES PARA CACO3 Y CL, VISTA
EN SECCIÓN E-W .......................................................................................... 67
FIGURA 5.5 ESTIMACIÓN MODELO DE BOQUES PARA CL EN CONA, VISTA EN
PLANTA .......................................................................................................... 68
FIGURA 5.6 ESTIMACIÓN MODELO DE BLOQUES PARA CL EN CONB, VISTA
EN PLANTA .................................................................................................... 68
FIGURA 5.7 ESTIMACIÓN MODELO DE BLOQUES PARA CL, VISTA EN PERFIL
E-W CON SC-04 DESPLEGADO.................................................................... 69
FIGURA 5.8 CATEGORÍA DE LA ESTIMACIÓN DEL MODELO DE BLOQUES
PARA CONA. .................................................................................................. 70
FIGURA 5.9 CATEGORÍA DE LA ESTIMACIÓN DEL MODELO DE BLOQUES
PARA CONB ................................................................................................... 70
FIGURA 5.10 DATA DRIFT CAO EN DIRECCIÓN N-S ........................................ 71
FIGURA 5.11 DATA DRIFT PARA CAO EN DIRECCIÓN E-W ............................ 72
FIGURA 5.12 DATA DRIFT PARA CAO EN COTA VERTICAL ............................ 72
FIGURA 5.13 VALIDACIÓN CRUZADA PARA CAO ............................................ 73
FIGURA 5.14 CURVA TON-LEY CACO3, PARA CONA ...................................... 74
FIGURA 6.1 MODELOS GEOLÓGICOS PARA CONA MELÓN 2017 Y SIGA 2005,
3D. ................................................................................................................... 76
FIGURA 6.2 MODELOS GEOLÓGICOS PARA CONB MELÓN 2017 Y SIGA 2005
........................................................................................................................ 77
vii
FIGURA 6.3 MODELOS GEOLÓGICOS SUPERPUESTOS, MELÓN 2017 Y SIGA
2005, 3D. ......................................................................................................... 77
FIGURA 6.4 DIFERENCIA SUPERFICIES ENTRE MODELOS MELÓN 2017 Y
SIGA 2005, PERFIL E-W ................................................................................ 78
FIGURA 6.5 CURVA TONELAJE-LEY CACO3 EN CONA, SIGA 2005 ................ 79
viii
1. Introducción.
1
coquinas del Yacimiento La Niña, sin embargo, no hay información de la variable Cl
en este yacimiento, debido a que, en las campañas de sondajes anteriores se utilizó
agua salobre para recuperar los testigos, contaminando con Cl marino las muestras,
invalidando los resultados obtenidos para este elemento.
Ante el poco conocimiento geológico del yacimiento La Niña, y la falta de
información sobre la variable Cl, es necesario realizar una nueva campaña de
sondajes utilizando agua potable, que tendrá como objetivos principales aumentar
la información que se tiene sobre el yacimiento, para así poder rehacer el modelo
geológico de manera mucho más precisa y poder limitar el cuerpo calcáreo con
mayor confiabilidad.
Con los resultados de este trabajo, la empresa Cementos Melón SA, tendrá una
base mucho más sólida para realizar una producción eficiente de clinker, además
de todos los procesos relacionados con la extracción de minerales del yacimiento
La Niña, ya sea el diseño de la mina, planificaciones mineras, evaluaciones de
proyecto, etc., lo que evidencia la importancia de realizar este trabajo.
2
1.2 Objetivos.
.
Ante la problemática existente, se plantean los siguientes objetivos principales y
secundarios para este trabajo, los cuales se realizarán en el yacimiento La Niña,
ubicado en la zona de Tongoy, Región de Coquimbo.
3
1.3 Área de estudio
Este yacimiento tiene como componente principal una roca sedimentaria que
se comporta como un manto subhorizontal de potencia 20 m aproximadamente,
compuesta por bioclastos de tamaño arcilla y arena principalmente, la cual se
encuentra parcialmente cementada por calcita a lo largo de toda su extensión.
Este manto de calizas o coquinas, al cual desde aquí en adelante se le
llamara “Conchuela”, sobreyace a arenas silíceas y subyace a un manto de arcillas
eólicas. Este paquete pertenece a la serie marina Formación Coquimbo de edad
Miocena- Pleistocena (Le Roux et al, 2006).
4
1.4 Metodología
Para lograr los objetivos propuestos, es necesario dejar en claro los limites en
cuanto a información que se utilizarán, como también su delimitación espacial y
profundización sobre el tema en estudio.
En cuanto a la información que se desea abarcar, este trabajo se limitara a los
datos recopilados en las dos campañas de sondajes realizadas en el año 2001 y
2005, además de la información obtenida en la campaña de sondajes que se realizó
en este trabajo y al levantamiento topográfico del techo del manto calcáreo
observado en las paredes interiores del rajo en explotación
El yacimiento se separará en 5 unidades geológicas: Arcillas, Coquina de Alta
Ley (Conchuela Alta > 80% en CaCO3), Coquina de Baja Ley (Conchuela Baja < 80
% en CaCO3), Arenas y batolito basal.
La metodología de modelamiento geológico se limitará a las propiedades minera
y agrícola perteneciente a la empresa cementos Melón SA en donde se ubica el
yacimiento La Niña, sin extenderse a zonas aledañas. El modelo geológico se
realizará en base a la construcción de perfiles tanto para la dirección N-S como E-
W, utilizando el software Vulcan 10.
En cuanto a la estimación de recursos que será en función de CaO, solo
comprenderá el sólido correspondiente al manto calizo ubicado dentro de las
partencias mineras y agrícolas de la empresa. En relación con la caracterización del
𝐶𝑙, al igual que para el caso de 𝐶𝑎𝐶𝑂3, solo abarcara las litologías de interés
económico.
La estimación de recursos se pretendía realizar mediante un kriking ordinario o
simple, se realizará finalmente con el método geoestadísticos de inverso a la
distancia al cuadrado, para ambas unidades de calizas. Esto debido al poco número
de muestras y baja densidad de sondajes que se tiene para desarrollar estos
procesos. Para encontrar los parámetros y radios de búsqueda de estimación de
recursos, se realizará un análisis de estadística descriptiva y variográfico
respectivamente, para las dos unidades de conchuela.
Para comprender de mejor manera la variable 𝐶𝑙 se realizará un análisis
univariante para determinar poblaciones (génesis). Para determinar la distribución
espacial de la variable 𝐶𝑙 se estimará con inverso a la distancia al cuadrado de la
misma forma que para 𝐶𝑎𝐶𝑂3.
6
2. Marco teórico
7
2.1.1 Materias primas
Los elementos principales que debe contener la mezcla o crudo son el Óxido
de Calcio (𝐶𝑎𝑂), el Dióxido de Silicio (𝑆𝑖𝑂2), el Óxido de Aluminio (𝐴𝑙2 𝑂3), y el Óxido
de Hierro (𝐹𝑒2 𝑂3), los cuales tienen que estar relacionados entre sí en proporciones
pre establecidas, con el objeto de dar determinadas características al clinker. Estos
pueden provenir de diferentes minerales. Las calizas que aportan el 𝐶𝑎𝑂 están
fundamentalmente compuestas por carbonato de calcio y óxido de silicio, teniendo
niveles muy bajos de carbonato de magnesio, alúmina y óxido de hierro. La
denominada caliza alta es rica en carbonato de calcio, presentando niveles de sílice
relativamente bajos. A medida que disminuye el contenido de calcio, y aumenta el
contenido de sílice, la caliza pasa a ser tipificada como intermedia y finalmente como
baja. Las arcillas y cenizas aportan 𝑆𝑖𝑂2 y el 𝐴𝑙2 𝑂3 que también puede provenir de
bauxitas. La pirita o hematita aporta el 𝐹𝑒2 𝑂3, etc.; teniendo que proceder a
mezclarlos previamente.
Sílice 𝑆𝑖𝑂2 20 %
Óxido de 𝑁𝑎 y 𝐾 0.3 %
8
En forma simplificada podemos decir que en los distintos tipos de cemento
su composición varía entre:
• 55 – 65% de 𝐶𝑎𝑂 + 𝑀𝑔𝑂
• 20 - 25% de 𝑆𝑖𝑂2 + 𝑆𝑂3
• 8 - 15% de 𝐴𝑙2 𝑂3 + 𝐹𝑒2 𝑂3
• Resistencia mecánica
La figura 2.1 muestra la variación de la resistencia a la compresión en pastas de
cada uno de los cuatro componentes principales del clinker. La resistencia en los
primeros tiempos se debe principalmente al C3S, pero, cuando prosigue la
hidratación, el efecto del C2S se vuelve cada vez más importante. Una alta
resistencia inicial puede obtenerse así, aumentando el contenido de C3S.
10
10000
6000 C2S
C3A
4000 C4AF
2000
0
0 90 180 270 360
Tiempo (días)
• Desprendimiento de calor
La hidratación del cemento Portland es un proceso exotermico y la velocidad de
desprendimiento de calor y la cantidad de calor liberado pueden ser importantes,
especialmente en la construccion de presas y otras estructuras grandes de
hormigón. Las fases C3A Y C3S tienen las mayores contribuciones al calor de
hidratación; en los cementos Portland de bajo calor, las proporciones de estos
componentes son reducidos, sin embargo, al tener menores concentraciones de
C3A y C2S disminuye su resistencia mecanica inicial, pero no la final. La figura 2.2
muestra resultados típicos correspondientes a las elevaciones de temperatura en
hormigones curados adiabaticamente para cemento Portland de endurecimiento
rápido (A), cemento Portland normal (B) y cemento Portland de bajo calor (C).
50°
Elevación de temperatura (°C)
40°
30° A
B
20°
C
10°
0°
0 7 14 21 28
Tiempo (días)
11
• Resistencia a los sulfatos
El hormigón fabricado en base al cemento Portland, es susceptible a ser
atacado por soluciones de sulfatos, que pueden causar su expansión y
desintegración. Algunas corrientes de aguas superficiales, del mismo modo que las
aguas del mar contienen suficientes sulfatos y plantean serios problemas. La
facilidad de ataque depende primariamente del contenido de C3A. Los cementos
resistentes a los sulfatos contienen proporciones reducidas de dicho componente,
o no lo tienen en absoluto; un rápido enfriamiento del clinker cuando deja el horno,
favorece la resistencia a los sulfatos, ya sea por que origina menor proporción de
C3A o porque afecta el tamaño o la perfección de los cristales.
• Finura
La disolución de un sólido en un líquido es tanto más rápida y tanto más
completa cuanto más fino es el sólido. Esto es lo que sucede con la hidratación del
cemento. Cuanto más fino sea un cemento, mayor es el número de partículas
aglomerantes presentes, y por lo tanto el número de uniones, dando lugar a un
conjunto más resistente e impermeable. La velocidad de hidratación incrementa al
aumentar la finura, y, por ende, el desarrollo de resistencias es más rápido, con una
mayor liberación de calor.
Aproximadamente un 95% del clinker está compuesto por una mezcla de los
cuatro minerales, C3S, C2S, C3A y C4AF, mientras que el resto está constituido por
cal libre (𝐶𝑎𝑂 libre), óxido de magnesio (𝑀𝑔𝑂) y sulfatos alcalinos.
• Cal libre
El contenido de cal libre en el clinker está relacionado con el grado de cocción
del mismo. Un elevado contenido de cal libre implica un clinker mal quemado, es
decir, con un insuficiente grado de cocción. Esto resulta obviamente, en un nivel de
C3S menor, con el consiguiente deterioro de la resistencia. Además, puede llevar a
problemas de expansión, durante la hidratación de morteros u hormigones. Muchas
veces, el contenido de cal libre del clinker se utiliza como parámetro de control del
proceso, ya que el mismo refleja el grado de cocción obtenido
• Oxido de magnesio
El magnesio (periclasa), generalmente acompaña al calcio en la naturaleza.
Las calizas, están formadas por una mezcla de carbonato de calcio (calcita) y
12
carbonato de magnesio (dolomita). Contrariamente a lo que ocurre con el calcio, el
magnesio no se combina durante el proceso de clinkerización, sino que permanece
“inerte”. Sin embargo, desde el punto de vista de la calidad del cemento es un
componente peligroso, ya que al reaccionar con el agua forma el hidróxido de
magnesio, de mayor tamaño, existiendo por lo tanto riesgo de expansión. Por este
motivo, las normas técnicas establecen un límite al contenido posible de óxido de
magnesio en el cemento.
• Cloruros
Los cloruros, al igual que los álcalis y el azufre, forman compuestos volátiles
que son arrastrados por los gases calientes del horno, que pueden condensar en
las partes frías del sistema, recircular y depositarse gradualmente hasta formar
obstrucciones que lleven a parar el horno. Además, los cloruros facilitan la corrosión
de las armaduras de hierro, lo cual es muy peligroso y perjudicial para la elaboración
de hormigón. Por estos motivos, las normas técnicas limitan el contenido posible de
cloruros en el cemento.
13
2.1.5 Fraguado y endurecimiento
El yeso retarda el fraguado inicial del cemento (de molerse el clinker sin yeso
fragua instantáneamente), lo cual se debe a que el aluminatotricálcico, en presencia
del yeso, no se hidrata de inmediato, sino que reacciona con el yeso produciendo
con ello una doble sal cristalina o sulfoaluminato (3𝐶𝑎𝑂. 𝐴𝑙2 𝑂3 . 3𝐶𝑎𝑆𝑂4 . 3𝐻2 𝑂),
manteniéndose así la concentración de aluminatos tricálcico por debajo de la
necesaria, para que así empiecen a separarse sus hidratos cristalinos. La adición
de yeso no debe superar el 3%, porque de lo contrario, se dilata dentro de él
produciendo fisuras.
14
2.2 Problemática del Cl
Esta problemática obliga a Cementos Melón SA. realizar una nueva campaña
de sondajes, que se realizará con agua potable, donde se recalculará el modelo
geológico y su estimación de recursos con mayor precisión, incluyendo esta vez el
elemento Cl, para así poder modelar las concentraciones de este elemento dentro
del Yacimiento La Niña.
15
2.3 Modelamiento Geológico
19
Qc: Depósitos coluviales (Cuaternario); Qa: Depósitos aluvionales (Cuaternario); Qm(a): depósitos de playas
marinas (Cuaternario); Qf2: Depósitos fluviales recientes (Cuaternario); Qf1: Depósitos fluviales antiguos
(Cuaternario); Qe: Depósitos eólicos (Cuaternario); MPlj: Depósitos eólicos de Quebrada el Jardín (Mioceno-
Pleistoceno); MPlc(a): Formación Confluencia, conglomerados (Mioceno-Pleistoceno); MPlc(b): Formación
Confluencia, brechas (Mioceno-Pleistoceno); MPlcq(a): Formación Coquimbo, areniscas finas limosas
(Mioceno-Plesitoceno); MPlcq(b): Formación Coquimbo, calizas bioclásticas (Mioceno-Pleistoceno); Kh:
Intrusivos Hipabisales andesíticos basálticos-dioríticos (Cretácico Superior); Kigd: Intrusivos graníticos a
dioríticos (Cretácico Inferior); Kig(a): Intrusivo de granodioritas (Cretácico Inferior); Kig(b): Intrusivo de
sienogranitos (Cretácico Inferior); Jmg: Intrusivos monzodioríticos-granodioríticos (Jurásico Medio-Superior);
JKd: Intrusivos dioríticos (Jurásico-Cretácico); Jas: Complejo Volcánico Agua Salada (Jurásico); Trat:
Complejo Plutónico Altos de Talinay (Triásico Superior); Pzch: Complejo Metamórfico del Choapa (Devónico-
Carbonífero). Fuente: Emparan y Pineda (2006).
20
3.1.1 Formaciones geológicas
22
ii. Constituye una costra cementada de superficies aterrazadas con
desarrollo de palolíneas de playa (25 m de espesor frente a Puerto
Velero), y sobreyace a la facies (a) con contacto por erosión; se interpreta
como parte de un evento regresivo. Para estos depósitos, se ha medido
un espesor ´máximo´ de 50 m en la ladera oeste de la quebrada Pachingo
(López, 1965).
Los depósitos de esta formación representan un ambiente litoral marino en
sus diversas facies. En general, la asociación de facies (a) representaría
depósitos de planicie y de estuarios en régimen macromareal, donde los
sedimentos gruesos corresponderían a depósitos de crecidas de canales de los
ríos confluentes. Las areniscas laminadas se interpretan como depósitos de
playa con los horizontes de “fierrillo” (concentración de minerales pesados)
desarrollados por el oleaje intramareal. La asociación de facies (b) se interpreta
como depósitos de playas dominadas por el oleaje y de canales de erosión. Los
depósitos de la asociación de facies (b), que forman líneas de paleoplayas sobre
diversos niveles estratigráficos de esta formación, se interpretan como
regresivos. Lo que representa el retiro paulatino del mar, debido a un alzamiento
(tectónico) continental, iniciado tentativamente en el Pleistoceno medio (Leonard
y Wehmiller, 1992; Frassinetti, 2001b).
24
FIGURA 3.3 GEOLOGÍA LOCAL CAMPAÑA 2001 (S1-S21), S OTO
(2005), VISTA EN PLANTA.
25
FIGURA 3.5 SECCIÓN N-S, ACUÑAMIENTO LATERAL DEL MANTO AL SUR , S OTO 2005)
27
Los resultados del cálculo de recursos son presentados en la tabla 3.2. Estos
valores se restringen a los modelos geológicos de las conchuelas de alta y baja ley
que se encuentran al interior de la propiedad minera que controla Cemento Melón
S.A.
TABLA 3.1 LEYES Y TONELAJES , SEGÚN UNIDADES GEOLÓGICAS (S OTO 2005)
Categoría Ley CaCO3 (%) Ton (Millones)
T_conchuela_alta_este_prop.00t Probadas 87,9 15,2
T_conchuela_alta_este_prop.00t Probables 86,0 7,2
Total 87,3 22,4
T_conchuela_baja_este_1_prop.00t Probadas 73,1 1,8
T_conchuela_baja_este_1_prop.00t Probables 73,9 0,6
Total 73,3 2,4
T_conchuela_baja_este_2_prop.00t Probadas 73,4 1,0
T_conchuela_baja_este_2_prop.00t Probables 75,2 1,5
Total 74,5 2,5
Total Unidades Controladas por Melón 84,9 27,3
30
25
20
Toneladas (millones)
15
10
0
84 86 88 90 92 94 96 98
Ley CaCO3 (%)
FIGURA 3.8 CURVA TONELAJE -LEY PARA CONCHUELA ALTA LEY (S OTO205)
29
4.1.2 Armado malla de sondajes
SC13
30
4.1.3 Recuperación de testigos
Durante el mes de febrero y marzo, la empresa Basalto Drilling SA, que fue
adjudicada para la ejecución de la campaña planeada por Geología de Minera
Melón SA, realizó los 12 sondajes presupuestados inicialmente, llamados
“SONDAJES COMPLEMENTARIOS” o SC-01, SC-02, …. al SC-12, habiéndose
ejecutado un sondaje extra SC-13 para completar el conocimiento general del
yacimiento que actualmente se explota para la provisión de Conchuela de Alta Ley
con destino exclusivo la Planta de cemento de La Calera.
La perforación en esta campaña fue realizada utilizando una Sonda canadiense
Versadrill Fordia Modelo Eider 1100 con motor Cummins. Capacidad hasta 600
metros en diámetro HQ (61 mm de testigo obtenido), pesando 8 Ton. Montada sobre
chasis a medida, con bomba de lodos FMC Bean, como se observa en la Figura
4.2. La metodología de perforación se realizó de dos formas, con tricono destructivo
sin recuperación de testigos para sobrecarga, y con sistema HQ3 para la
recuperación de testigos del manto calcáreo.
El Equipo Profesional al cual se le encomendó esta campaña, permaneció el 100 %
del tiempo controlando la calidad de las muestras, del agua con la que se produjeron
los lodos de perforación, el encajonado de los testigos e identificación, por otro lado,
el personal de Basalto Drilling SA, demostró grandes habilidades y conocimiento
para realizar la tarea encomendada, alcanzando un 100 % en las maniobras de
rescate de la barra porta-testigos. Con esto, se garantiza la calidad, validez y
representatividad del trabajo ejecutado.
.
FIGURA 4.2 MONTAJE SONDAJE SC-02
31
Resumen de las operaciones
• 310,18 m totales perforados en arcillas de escarpe.
o 268,8 m totales perforados con tricono destructivo para pasar las
arcillas eólicas superiores (escarpe).
o 42,18 m totales perforados con sistema HQ3 para garantizar el
contacto geológico exacto entre arcillas y conchuela.
• 261, 97 m totales perforados de conchuela alta y baja ley obtenidos como
testigos indisturbados HQ3.
• 31,7 m totales perforados de arenas de playa inferiores (sobre las que se
depositó la conchuela), estériles en 𝐶𝑎𝐶𝑂3, obtenidas como testigos
indisturbados HQ3 para garantizar contacto exacto entre conchuela y arenas
inferiores.
• Ningún sondaje de esta nueva campaña toco basamento, o algún bloque de
granito.
• Total, metros perforados con ambos métodos (tricono destructivo y testigos
con triple barril HQ3) =604,65 m.
• NOTA 1: Todos los sondajes atravesaron el 100% del espesor real de la capa
de conchuelas, tanto de alta ley como de baja ley, desde su escarpe hasta
las arenas basales de playa que ya no contienen conchuelas ni bioclastos.
• NOTA 2: La recuperación de los testigos alcanzo el 100% del espesor de las
capas perforadas con barra triple HQ3.
La Tabla 4.2, resume lo encontrado en cada uno de los 13 sondajes.
32
4.1.4 Descripción, armado de muestra y análisis químico.
Una vez recuperados los testigos, son compositadas las muestras según
cambios litológicos notables. Las premisas usadas para este procedimiento son, si
la roca era similar en un tramo muy largo, se definía una muestra cada 1.5 m de
largo o espesor. Si, por otro lado, el cambio litológico en forma de intercalación fina
no superaba los 30 cm, se la integraba a la muestra mayor. De lo contrario (>0.3 m)
se la definió como una muestra independiente. La Figura 4.3 muestra los testigos
cortados longitudinalmente para ser compositados y descritos según los criterios
propuestos.
33
4.2 Modelamiento geológico
34
TABLA 4.5 BASE DE DATOS COLLAR
35
4.2.2 Modelamiento geológico
36
Para definir las superficies de los techos y bases de ConA y ConB, se
procede a la confección perfiles geológicos que delimitaran los contactos geológicos
de interés, estos perfiles se realizaran en dos juegos de secciones (Figura4.5), uno
en la dirección Norte, con secciones cada 30 m, y el otro en la dirección W, con
perfiles cada 60 m, esto debido a que la malla de sondajes tiene una distancia
aproximada de 160 m entre cada sondaje. Estas distancias son propuestas para
alcanzar una buena densidad de perfiles y así un modelo más preciso.
Es posible observar en las paredes del rajo del yacimiento La Niña, un manto
de conchuela alta y baja continuo, además del claro contacto continuo entre
Conchuela Alta y Arcilla a lo largo de toda la mina, esto sumado a la similar
estratigrafía encontrada en todos los sondajes, son argumentos suficientes para
asegurar una continuidad litológica de ConA y ConB a la hora de realizar cada uno
de los perfiles geológicos.
El primer paso es elaborar los perfiles Norte para las secciones propuestas,
donde solo nos concentramos en delimitar los contactos entre las dos unidades ya
especificadas. Esta delimitación de litologías para cada sección se presenta en la
Figura 4.6 (perfil AA’ pertenece a la Figura 4.5), donde ConA aparece delimitado por
el color celeste, y ConB por el color rosa.
37
FIGURA 4.6 PERFIL LITOLÓGICO AA', PARA CON A Y CONB.
Quedando finalmente las secciones Norte, tanto para ConA y ConB como se
presenta en la Figura 4.7.
38
Con los interceptos definidos, se pueden crear los perfiles para cada sección
Este. La Figura 4.9 presenta la construcción del perfil BB’ (Figura 4.5) para ambas
litologías. De la misma manera que para las secciones Norte, se construyen todas
las secciones Este, obteniendo el juego de perfiles Este.
39
Modeladas las superficies de cada unidad geológica, se procede a juntar
techos y bases respectivos a cada litología para generar los sólidos o modelos
geológicos de ambas unidades geológicas, sin embargo, para poder obtener los
modelos geológicos finales, faltaría limitar los modelos de ConA y ConB (Figura
4.10) a la topografía actual del yacimiento (superficie café), al batolito (superficie
roja), y tanto a la propiedad minera como agrícola (línea roja y rosa
respectivamente).
Debido a la actividad y explotación del yacimiento La Niña, se limita a la
topografía presente el mes de mayo del 2017. En cuanto al batolito, fue modelado
como una superficie irregular, en base a secciones Norte cada 10 m, las cuales
fueron elaboradas a partir de los datos de sondajes y mapeo del batolito en
superficie.
40
N
41
N
43
FIGURA 4.15 LOGNORMAL CAO, CON A
44
FIGURA 4.17 BOXPLOT CL
45
FIGURA 4.19 LOGNORMAL CL, CONB
47
FIGURA 4.21 CORRELACIÓN CAO-F E2O3 EN CON A Y CONB
48
FIGURA 4.22 CORRELACIÓN CAO-AL2O3 PARA CON A Y CONB
49
FIGURA 4.23 CORRELACIÓN CAO-MGO PARA CON A Y CONB
50
FIGURA 4.24 CORRELACIÓN CAO-CL PARA C ON A Y CONB
51
FIGURA 4.25 CORRELACIÓN CL-SO3 PARA CON A Y CONB
52
Un resumen de los coeficientes de correlaciones se presenta a continuación
en la Tabla 4.7.
TABLA 4.7 RESUMEN CORRELACIONES
53
4.3.3 Análisis variográfico
54
FIGURA 4.26 CARACTERIZACIÓN VARIOGRAMAS
55
FIGURA 4.28 V ARIOGRAMA C AO PARA CONB
56
FIGURA 4.30 V ARIOGRAMAS SI O2 PARA CONB
57
FIGURA 4.32 V ARIOGRAMAS CL PARA CONB
58
Donde cada parámetro del variograma mencionado, hace referencia a una
característica especifica en el cono de búsqueda del variograma, como se muestra
en la Figura A.2 del capítulo de anexos.
Construidos los variogramas, y sus modelos esféricos ajustados, se presenta
la Tabla 4.9, que resume los alcances de los radios de búsqueda para cada una de
las tres direcciones, de los compuestos químicos en estudio, en ambas unidades
geológicas.
TABLA 4.9 ALCANCES VARIOGRAMAS
59
un alcance de 60 m en esta dirección. Así, al igual que en ConA, se descarta la
posibilidad de realizar un Kriking ordinario, debido a los bajos alcances de búsqueda
y la falta de un variograma aceptable en el eje Semi.
Ahora, en los variogramas encontrados para 𝑆𝑖𝑂2, en ConA, se observa un
variograma pobre en el Major (Figura 4.29.A), aunque es posible modelar el
variograma y estimar un alcance de 100 m. El eje Semi (Figura 4.29.B), presenta un
variograma muy malo que no llega a la meseta, estimando un alcance de 50 m. Para
el eje Minor (Figura 4.29.C), se observa un muy buen variograma que nos otorga un
alcance de 16 m. Es posible percibir que los variogramas nos otorgan resultados
muy similares para el eje Major y Minor con el caso de CaO, lo cual tiene sentido
según lo visto en las correlaciones entre estos dos elementos presentados en la
estadística descriptiva. Para los variogramas de 𝑆𝑖𝑂2 en ConB, tanto para el eje
Major y Semi (Figura 4.30.A y 4.30.B), se tienen variogramas bastante pobres,
obteniendo alcances de 100 m y 50 m respectivamente, sin embargo, para el eje
Minor (Figura 4.30.C), se obtiene un variograma bastante decente, resultando un
alcance de 6 m. Al igual que para ConA, los alcances son muy similares al caso de
CaO en ConB, lo cual tiene sentido con las correlaciones entre estos elementos
para esta unidad geológica. De esta información recopilada, es posible desprender
que ante los bajos alcances obtenidos y los pobres variogramas resultantes, se
descarta la posibilidad de realizar un Ko, también se percibe que la correlación
presente entre 𝑆𝑖𝑂2 y CaO, respeta de forma aproximada, los alcances obtenidos
de los variogramas para estos dos elementos en cada unidad geológica.
60
4.3.4 Criterios de estimación
Con los criterios definidos para realizar una estimación optima del yacimiento,
se procede a la construcción del modelo de bloques, que será el soporte para la
estimación de nuestras variables en estudio. Las dimensiones del modelo de bloque
se definen en función de la planificación minera del yacimiento La Niña, el cual
considera bancos para el rajo de 5 m de altura, ante esto, el modelo de bloques que
se construye para la estimación de recursos, se define con un tamaño de bloques
de 20x20x5 m, con sub bloques de 5x5x1.25 m, para delimitar de mejor manera las
irregularidades de las superficies de los modelos geológicos. La Figura 4.33 ilustra
los modelos de bloques para ambas litologías.
63
• Arenas Basales de Playa: Subyacen a la conchuela de baja ley, estériles
en 𝐶𝑎𝐶𝑂3. Fueron perforados hasta unos 5 metros para garantizar que no
hubiese un resto de conchuela útil por debajo.
La Figura 5.1, muestra la columna estratigráfica del sector Sur del yacimiento La
Niña, esta columna característica al yacimiento en su totalidad, solo con la diferencia
que el espesor de la sobrecarga disminuye hacia el Norte.
64
5.2 Estimación del modelo de bloques.
65
N
66
La Figura 5.4, presenta una vista en perfil para entender el comportamiento
de la estimación en un modelo de bloques en profundidad, donde es posible
observar una clara cercanía entre el valor estimado de CaCO3 para cada bloque
(N° celeste) y la litología asociada al sondaje SC-02 desplegado con el valor del
composito correspondiente (N° rojo). Esto sirve como validación visual del modelo
de bloques estimado. También, es posible observar como los sub bloques nos
ayudan a afinar los límites de la superficie de las unidades geológicas modeladas.
Para el caso de la estimación del Cl, tanto para ConA como ConB (Figura 5.5
y 5.6), se observa falta de información en el sector Norte del yacimiento, esto es
debido a la ausencia de análisis de Cl en los sondajes Norte realizados en las
campañas de sondajes anteriores, donde los análisis de Cl fueron errados. A esto
se complementa, que solo 3 sondajes de los 13 de la campaña del 2017 arrojaron
muestras con altos contenidos de Cl, los cuales se evidencian en las pequeñas
zonas con altos contenidos de Cl.
67
N
68
La Figura 5.7, muestra el sondaje que tiene más compositos con altos
contenidos de Cl (> 25 ppm), donde se observa que no existe una continuidad
vertical de estas altas concentraciones, corroborando la suposición de que el Cl se
presenta de forma aleatoria en el yacimiento, afirmación que es apoyada por la
estadística descriptiva y variogramas realizados sobre la variable Cl. Sin embargo,
es posible observar una muy buena similitud entre los valores obtenido entre la
estimación del modelo de bloque de la variable Cl, y los valores de las
concentraciones de Cl de cada composito (N° rojo), lo cual nos da una validación
visual de la estimación realizada.
FIGURA 5.7 ESTIMACIÓN MODELO DE BLOQUES PARA CL, VISTA EN PERFIL E-W CON
SC-04 DESPLEGADO
69
N
70
5.3 Validaciones estimación
71
FIGURA 5.11 DATA DRIFT PARA CAO EN DIRECCIÓN E-W
72
En cuanto a la validación cruzada presentada en la Figura 5.13, se observa
una correlación del 73% entre la estimación y el valor real medido para cada
composito. Así, mediante este método se corroboraría que la estimación representa
con una confiabilidad aceptable en relación con los datos obtenidos del análisis de
Fluorescencia de Rayos X de CaO de todos los compositos medidos.
73
5.4 Calculo de recursos
74
La tabla 5.1, nos muestra las toneladas correspondientes a ConA y ConB,
donde resulta que las toneladas presentes en la unidad geológica Conchuela de alta
ley, es de 15.02 millones de toneladas, con una ley media del 90.6%. Para
Conchuela de baja ley, se obtienen 13.9 millones de toneladas con una ley media
del 75.3%, resultando un total de 28.94 toneladas de conchuela. Sin embargo, ante
el requerimiento por parte de Cementos Melón SA. de una ley media del 92%
aproximadamente en La Niña para una eficiente producción de clinker, la Fig.23
como se mencionó en el párrafo anterior, indica que para obtener una ley media del
92% en ConA, se necesita un cutoff del 86,5%, donde obtendríamos un recurso de
12.3 millones de toneladas de 𝐶𝑎𝐶𝑂3.
En la práctica, para alcanzar estas leyes, sectores con leyes > al 92%, son
mezclados en el día a día de la faena con conchuela de baja ley, para así, enviar
camiones cargados con producto lo más homogéneo (poco variable) posible, de
esta forma se mejora el uso de esta materia prima para la fabricación de un clinker
uniforme.
TABLA 5.1 RESUMEN TONELADAS FINALES CON A Y CONB
75
6. Discusiones
FIGURA 6.1 MODELOS GEOLÓGICOS PARA CON A MELÓN 2017 Y S IGA 2005, 3D.
76
N
FIGURA 6.2 MODELOS GEOLÓGICOS PARA CONB MELÓN 2017 Y SIGA 2005
FIGURA 6.3 MODELOS G EOLÓGICOS SUPERPUESTOS, MELÓN 2017 Y SIGA 2005, 3D.
77
En la Figura 6.3.A, es posible observar las discrepancias entre las superficies
de ConA, esto es debido a las diferencias en la elaboración de las secciones cuando
se modela el sólido, como se muestra en la Figura 6.4.
En cuanto a las diferencias de modelos geológicos para ConB (Figura 6.3.B),
el modelo geológico SIGA 2005, se corta a la mitad en ConB norte y ConB sur, esto
se debe a la falta de interpretación entremedio producto de la falta de sondajes. Con
la nueva campaña de sondajes 2017 se tiene más información del techo y piso de
ConB, pudiendo definir la continuidad norte – sur de la misma.
FIGURA 6.4 DIFERENCIA SUPERFICIES ENTRE MODELOS MELÓN 2017 Y SIGA 2005,
P ERFIL E-W
78
6.2 Comparación estimación SIGA (2005)
TABLA 6.1 RESUMEN TONELADAS FINALES DE CON A Y CONB PARA SIGA 2015
79
Según la estimación de recursos recalculada para el modelo de bloques de
SIGA 2005, resulta que ConA tiene 15,7 millones de toneladas con ley media de
90.5% de 𝐶𝑎𝐶𝑂3. Ahora, según esta estimación de recursos, la curva Tonelaje- Ley
(Figura 6.5), arroja para una ley optima media del 92% con ley de corte del 88,5%
de 𝐶𝑎𝐶𝑂3, una cantidad de recursos de 11.6 millones de toneladas.
La tabla 6.2, presenta la comparación final entre el cálculo de recursos de
SIGA 2005 y Melón 2017, de las toneladas existentes en ConA para ambos modelos
geológicos, tanto para ConA total, o sea, ley de corte 80% debido a que con esta
ley mínima fue definida ConA (Tabla 6.2.A) como para las toneladas con ley optima
media del 92% (Tabla 6.2.B)
80
7. Conclusiones
81
Al comparar el cálculo de recursos con el realizado con la consultora SIGA 2005,
se observa una pérdida de 700 toneladas en ConA, producto del aumento de
volumen del batolito modelado en este trabajo, sin embargo, los recursos para una
ley media optima del 92% de CaCO3 en ConA, aumentan en 700 toneladas, esto se
debe a las altas leyes encontradas en la campaña de sondajes 2017, lo que provocó
un aumento en la concentración de 𝐶𝑎𝐶𝑂3 en los bloques al ser re evaluados.
En cuanto a la problemática del Cl, las concentraciones de este elemento se
localizan en áreas proporcionalmente reducidas y son susceptibles de manejo
durante la explotación minera, mezclando estas zonas de elevadas concentraciones
con zonas reconocidas con niveles bajos en Cl.
Desde el punto de vista práctico, a raíz de este trabajo, se podrá entregar a los
jefes de faena planos (en Auto CAD) de los diferentes niveles de explotación, con
las leyes esperables sector por sector, que funcionaría como herramienta de control
químico de la conchuela para una correcta explotación, garantizando precisión en
el control de esta.
Como recomendación, sin dudas, otra pequeña campaña de sondajes en unos
7 años más será necesaria para cubrir el sector Sur - Sur del Rajo, puesto que, para
la campaña que terminó, no hubo recursos económicos suficientes como para
abarcar esa área distal. Por otro lado, se propone utilizar métodos de “micro-
minería” (ataque a pequeños sectores con maquinaria de pequeño porte), para
recuperar el 100 % de las reservas de conchuela en el fondo irregular de la faena.
Por otro la
82
8. Bibliografía
83
• López, M.C. 1965. Estudio de los depósitos marinos de la Bahía de Tongoy.
Memoria para optar al título de geólogo (Inédito), Universidad de Chile, Escuela
de Geología, 174 p. Santiago.
• Marco Matías Pfeiffer Jackob (2011). Evolución y Genesis de calcretas
pedogénicas en la Paleobahía de Tongoy.
• Martínez, R.; Caro, R. 1980. Microfósiles silíceos de las diatomitas de Tongoy,
Provincia de Coquimbo, Chile: su significado biocronoestratigráfico,
biocronogeológico, paleoecológico y paleogeográfico. Revista Geológicade
Chile, No. 10, p. 33- 53.
• Paskoff, R. 1970. Le Chili semi-aride: recherches geómorphologiques. Biscaye
Frères, Imprimeurs, 420 p. Bordeaux, France.
• Pfeiffer, M. Evolución y génesis de calcretas pedogénicas en la paleobahía de
Tongoy. Tesis para optar al grado de magister en ciencias, mención Geología.
Universidad de Chile. Santiago, Chile.
• Pineda, G.; Emparan, C. En edición. Geología del área Vicuña-Pichasca. Región
de Coquimbo. Servicio Nacional de Geología y Minería, Carta Geológica de
Chile, Serie Geología Básica, 1 mapa escala 1:100.000.
• Saillard, M., Hall, S.R., Audin, L., Farber, D.L., Hérail, G., Martinod, J., Regard,
V., Finkel, R.C., Bondoux, F., 2009. Non-steady long-term uplift rates and
Pleistocene marine terrace development along the Andean margin of Chile
(31°S) inferred from 10Be dating. Earth and Planetary Science Letters, 277: 50–
63.
• SERNAGEOMIN (2006). Carta geológica Andacollo-Puerto Aldea, Región de
Coquimbo.
• SIGA 2001-2004, Informe Proyecto La Niña, Modelamiento y estimación de
recursos.
84
A. Anexo
85
• Toma de muestras, tiene importancia en aspectos como cantidad,
distribución y representatividad de las muestras.
• Visualizar muestras, determinar ubicación de las muestras.
• Análisis exploratorio de datos, Análisis mediante herramientas de
estadísticas tradicional de las muestras
• Desagrupamiento, definir las zonas de estudio. Una división del campo en
varias sub-zonas puede ser relevante si uno observa cambios abruptos en la
distribución espacial de valores, o si la geología del fenómeno lo indica.
• Predicciones (Simulaciones con Inverso a la distancia o Kriking)
86
A.1.1.3 Valor esperado
• 𝐸(𝑋+𝑌) = 𝐸(𝑋)+𝐸(𝑌)
• 𝐸(𝑐𝑋) = 𝑐+𝐸(𝑥), 𝑝𝑎𝑟𝑎 𝑡𝑜𝑑𝑜 𝑛𝑢𝑚𝑒𝑟𝑜 𝑟𝑒𝑎𝑙 𝑘
• 𝐸(𝑐) = 𝑐, 𝑝𝑎𝑟𝑎 𝑡𝑜𝑑𝑜 𝑛𝑢𝑚𝑒𝑟𝑜 𝑟𝑒𝑎𝑙 𝑐
• 𝐸(𝑎∗𝑋+𝑏) = 𝑎∗𝐸(𝑋)+𝑏, 𝑝𝑎𝑟𝑎 𝑡𝑜𝑑𝑜 𝑝𝑎𝑟 𝑑𝑒 𝑛𝑢𝑚𝑒𝑟𝑜 𝑟𝑒𝑎𝑙𝑒𝑠 𝑎 𝑦 𝑏
Esperanza del producto:
• 𝐸(𝑋∗𝑌) = 𝐸(𝑋)∗𝐸𝑛 𝑒𝑙 𝑐𝑎𝑠𝑜 𝑞𝑢𝑒 𝑋 𝑒 𝑌 𝑠𝑒𝑎𝑛 𝑣𝑎𝑟𝑎𝑖𝑏𝑙𝑒𝑠 𝑎𝑙𝑒𝑎𝑡𝑜𝑟𝑖𝑎𝑠 𝑖𝑛𝑑𝑒𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒
A.1.2 Momentos
87
A.1.2.2 Momento de segundo Orden:
A.1.3 Regionalización
88
• La distribución bivariable solo depende de la separación entre las muestras
(h)
Estos parámetros muestran la interacción que existe entre dos valores, luego dan
una descripción sintética de la continuidad espacial de la variable regionalizada.
La covarianza y el correlograma indican qué tan semejantes son los valores entre
dos sitios, mientras que el variograma indica qué tan semejantes son. Siendo
constante la esperanza, se tiene también.
• La varianza es igual a la función de covarianza evaluada para el vector H=0
A.1.4 Variograma
89
La fórmula de variograma:
Donde:
Z(x)= valor de la variable en el sitio x
Z(x+h) = valor de la variable en el sitio x+h
N= número de pares formados
90
Además de las direcciones tendrá inclinaciones tales como azimuth, dip, y
plunge entre otros parámetros necesarios para su utilización en Vulcan, todos estos
parámetros están explicados en la Figura A.2. Cada una de estas direcciones tiene
todos sus parámetros bien definidos.
• El modelo esférico:
Este modelo define que el efecto pepita es de alcance infinitamente corto o
pequeño. Es el modelo más comúnmente usado para depósitos minerales. Se
ajusta a variogramas que exhiben un comportamiento casi lineal desde el origen,
elevándose rápidamente mientras se va curvando gradualmente.
92
ejemplo; p es la potencia a la que se eleva el inverso de la distancia y que en general
está entre 1 y 2; z (𝑢𝛼 ) es el valor de la variable en el punto muestreado 𝑢𝛼 .
93
A.2 Análisis químicos
TABLA A.1 ANÁLISIS QUÍMICOS CAMPAÑAS ANTERIORES (HOLE ID S) Y
CAMPAÑA 2017 (HOLE ID SC)
94
95
96
97
98
99
100
101
102
103
104
105
A.3 Descripción mapeo sondajes.
106
TABLA A.3 MAPEO SONDAJE SC-03
107
TABLA A.4 MAPEO SONDAJE SC-05
108
TABLA A.5 MAPEO SONDAJE SC-10
109
TABLA A.7 MAPEO SONDAJE SC-04
110
TABLA A.8 MAPEO SONDAJE SC-06
111
TABLA A.9 MAPEO SONDAJE SC-01
112
TABLA A.10 MAPEO SONDAJES SC-02
113
TABLA A.12 MAPEO SONDAJE SC-11
114
TABLA A.13 MAPEO SONDAJE SC-08
115
TABLA A.14 MAPEO SONDAJE SC-07
116