Metodos de Geoestadistica PDF

98
195
METODO
GEOESTADISTICO
Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade
196
15
10 X = 0,93
5
S = 1,20
0
5 10 15 20 25
15
X = 0,93
10
S = 1,20
5
0
5 10 15 20 25
INTERCADE www.intercade.org
CONSULTANCY & TRAINING
99
197
¿Cómo poder relacionar los valores con sus

posiciones en el espacio?
¿Cómo relacionar dichos valores entre sí?
GEOESTADISTICA
Teoría de la variable regionalizada
La variable toma un valor en cada punto del espacio.
198
DISTRIBUCION ESTADISTICA
DE LAS LEYES KRIGEADO DE BLOQUES:
2
Zi + G E
YACIMIENTO
KRIGEADO
NIVEL H
NIVEL K
h x+h
x
z+ VALOR MEDIO ESTIMADO Z

2
Y VARIANZA ESTIMADA 6 E
CORRELACION ESPACIAL
VARIOGRAMA Z = ESTIMADOR LINEAL

LEYES DE LOS TESTIGOS OPTIMO E INSESGADO
DE LOS SONDEOS DE UN BLOQUE O UN
PUNTO, BASADO EN Z
LOS VALORES
CIRCUNDANTES Y EN
1 EL VARIOGRAMA
(h) = (f(x) - f(x+h)) 2
KRIGEADO PUNTUAL:
h MAPA DE ISOPLETAS
100
199
El variograma liga el valor de una variable con su posición

en el espacio, asumiendo que las muestras están
correlacionadas.
γ (h) = Σ [Z(xi) - Z(xi + h)]2 / 2n

(h)
SiII
Rango a
h(distancia)
• h: distancia entre los pares
• n: número de pares
• Z(xi): localización y valor de la muestra
200
Pero ¿cómo se construye un variograma?
100'
100 39 44 40 42 40 39 37 36
44 40 42 40 37 36
42 42 41 40 36
42 43 42 39 39 41 40 38 43 39 39
*(100)= [(40 - 42) 2 + (42 - 40) 2 + (40 - 39) 2 + (39 - 37) 2
+(37 - 36)2 + (43 - 42)2 + (42 - 39)2 + (39 - 39)2
37 35 38 37 37 34
+(39 - 37) 2 + (41 - 40)2 + (40 - 38)2 + (37 - 37)2
37 37 37 35 38 37 37 33 34 37 37 33
+(37 - 37) 2 + (37 - 35)2 + (35 - 38)2 + (38 - 37)2
100'
+(37 - 37)2 + (37 - 33)2 + (33 - 34)2 + (35 - 38)2

35 38 35 37 38 36 35 35 38 36 37 38 36 35 +(35 - 37)2 + (37 - 36)2 + (36 - 36)2 + (36 - 35)2
+(36 - 35)2 + (35 - 36) 2 + (36 - 35)2 + (35 - 34)2
+(34 - 33) 2 + (33 - 32)2 + (32 - 29)2 + (29 - 28)2
36 35 36 35 34 33 32 29 28 36 35 36 35 34 33 32 29 +(38 - 37)2 + (37 - 35) 2 + (29 - 30) 2 +
+(30 - 32)2 ] / ( 2x 36)
*(100)= 1.46(%) 2
38 37 35 30 29 30 32 38 37 35 30 29 30 32
44 40 42 40 39 37 36 (h)
*(200)= [(44 - 40) 2 + (40 -40 ) 2 + (42 - 39) 2 + (40 - 37) 2
2
Variograma experimental (%)
42 43 42 39 39 41 40 38 +(39 - 36) 2 + (42 - 43) 2 + (43 - 39) 2 + (42 - 39) 2

+(39 - 41) 2 + (39 - 40) 2 + (41 - 38) 2 + (37 - 37) 2
+(37 - 35) 2 + (39 - 40) 2 + (35 - 37)2 + (38 - 37)2 3.0
37 37 37 35 38 37 37 33 34
+(37 - 33) 2 + (37 - 34) 2 + (38 - 35) 2 + (35 - 36) 2
35 38 35 37 38 36 35
+(37 - 36) 2 + (36 - 35) 2 + (36 - 36) 2 + (35 - 35) 2
+(36 - 34) 2 + (35 - 33) 2 + (34 - 32) 2 + (33 - 29) 2
1.5
+(35 - 28) 2 + (38 - 35) 2 + (35 - 30) 2 + (30 - 29) 2
36 35 36 35 34 33 32 29 28
+(29 - 32) 2 ] / (2 x 33)
*(200)= 3.30(%) 2
38 37 35 30 29 30 32
100 200 h
γ (h) = Σ [Z(xi) - Z(xi + h)]2 / 2n

101
201
Variograma experimental
202
El crecimiento de la curva indica el grado de

continuidad de la mineralización.
700
c
525
125
(h) 350
100
c VARIOGRAMA
175
Co
75
0
0 18 36 54 72 90 106 126 m
(h) 50
(x 10') Figura 21 - Variograma de un pórfido cuprífero
VARIOGRAMA
25
0 100
0 9 18 27 36 45 54 63 m c
(x 10')
Figura 20 - Variograma de un yacimiento estratiforme (h) 75
(x 10 3 )
VARIOGRAMA
50
22.5 45 67.5 90 112.5 135 157.5 m

Figura 22 - Variograma de un yacimiento de oro
102
203
Los variogramas en varias direcciones ponen de

manifiesto la anisotropía en la distribución de leyes.
POCILLOS
GALERIAS
SECTOR ANISOTROPICO MODELOS SECTOR ISOTROPICO
160 852
Variabilidad (Oz Ag/ ton) 2
Variabilidad (Oz Ag/ ton) 2

GALERIAS
135 710
CHIMENEAS
110 568
GALERIAS
85 426
60 284
CHIMENEAS
35 142
(h) (h)
8 16 24 32 40 (m.) Distancia 8 16 24 32 40 (m.) Distancia
204
Variograma teórico
Variograma experimental
103
205
Variograma teórico
Meseta: valor máximo de variabilidad
Alcance: área de influencia de la correlación
Efecto pepita: discontinuidad en el origen
(h)
SiII
Efecto Pepita
Rango a
h(distancia)
206
VARIOGRAMAS TEORICOS
Esférico Gaussiano Exponencial Potencial (a<1)
1.2
0.8
0.6
0.4
0.2
0
0 50 100 150 200 250 300
h(m)
104
207
Modelo esférico
γ(h) = C [(3/2)(h/a) – 1/2(h/a)3] + Co h ≤ a
GAMMA (H)
MESETA
C + Co h>a
ALCANCE
DISTANCIA
Modelo exponencial
γ(h) = C [1 - Exp(-|h|/a)] + Co h ≤ a
C + Co h>a
GAMMA (H)
DISTANCIA
208
Modelo gaussiano
GAMMA (H)
γ(h) = C [1 - Exp(-|h|2/a2)] + Co h ≤ a
C + Co h>a
DISTANCIA
Modelo logarítmico
γ(h) = log h
GAMMA (H)
DISTANCIA
105
209
Efecto pepita
GAMMA (H)
γ(h) = C0
DISTANCIA
Modelo lineal
γ(h) = Ah
GAMMA (H)
DISTANCIA
210
El variograma expresa la
CORRELACION ESPACIAL
correlación espacial entre
(h)
los valores muestreados.
VARIOGRAMA Z
(h) =
1
(f(x) - f(x+h)) 2
Es la función de partida
h
para el método
geoestadístico.
106
211
KRIGEADO
Estimar valores desconocidos a

partir no solo de los conocidos,
z+ sino también de su estructura y
continuidad espacial.
Z + = ESTIMADOR LINEAL
OPTIMO E INSESGADO
DE UN BLOQUE O UN
PUNTO, BASADO EN Krigeage o kriging
LOS VALORES
CIRCUNDANTES Y EN
EL VARIOGRAMA
212
El valor de la variable se calcula como una combinación

lineal de los valores que presenta dicha variable en los
puntos vecinos.
Z* = ∑ ai * zi KRIGEADO
• Z*: valor de la variable a estimar

z+
• zi: valor de la variable en cada punto
• ai: pesos asignados en las muestras
Z+ = ESTIMADOR LINEAL
OPTIMO E INSESGADO
DE UN BLOQUE O UN
PUNTO, BASADO EN
LOS VALORES
CIRCUNDANTES Y EN
EL VARIOGRAMA
107
213
El krigeado es el mejor estimador de ai.
Insesgado: E [Z* - Z] = 0
∑ ai = 1
Varianza mínima: Var [Z* - Z] = σk2
VAR(Z* -Z) = σk2 = σv2 - 2 Σ ai σv Xi + Σ Σ ai aj σ Xi Xj

i=1 i=1 j=1
σ=γ
• σv2: varianza del bloque del volumen V

• σv Xi: covarianza entre el bloque y cada una de las muestras
• σ Xi Xj: covarianza entre las muestras
214
Sistema de ecuaciones
11 12 ... 1n a1 vx1
21 22 ... 2n a2 vx2
= A = D =
n1 r2 ... nn an vxn
... 0
Las varianzas y covarianzas son datos obtenidos del

variograma.
σ=γ
Las incógnitas son los coeficientes de ponderación.

a1, a1 … an
Se introducen en la siguiente ecuación:

Z* = ∑ ai * zi
108
215
VENTAJAS E INCONVENIENTES
La aplicabilidad es muy amplia.

Evita la ponderación arbitraria.
Tiene gran cantidad de información necesaria.
Los datos deben estar distribuidos uniformemente.
No se comprueban los resultados.
Estos modelos no tienen la participación de expertos.
Con datos incorrectos o insuficientes se tienen
variogramas no fiables.
216
EJEMPLO 1
5 A 2
1
109
217
¿Cómo determinamos la ley en el punto A?
A
5 2
a = 100 m 1
El variograma determina qué muestras puede tener una

influencia real en la estimación.
218
Los cálculos son realizados por programas

especializados.
Las opciones relativamente económicas son

© ©
programas como Surfer8 o EcoSSe , que no permiten
el diseño ni el estudio de bloques.
Puede realizar una buena modelización de

variogramas experimentales y desarrollar kriging
puntual.
Se pueden obtener mapas, donde la interpolación de

valores en el espacio XY está controlada por la función
γ(h).
110
0 0
219
0 0 0
0 0 1 0
0 1
0 2 3
15
0 0 1
0 2
0 4
0 2 0 0
1 0
10 3 2
0
3 0 1
1 3 3
0
2 2 0 3
5 3 2
1 2 0 0
0 2 1
0 1 0 0
0 0 0
0
0
0 5 10 15 20 25
Column C
Direction 0.0 Tolerance 90.0
2.2
1.8
1.6
1.4
Variogram
1.2
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9 10
Lag Distance
220
Column C
Direction 0.0 Tolerance 90.0
2.2
1.8
1.6
1.4
Variogram
1.2
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9 10
Lag Distance
111
221
222
112
223
EJEMPLO 2
Se tiene un conjunto de cuatro muestras de un
yacimiento de plomo, cuyas leyes se muestran en la
tabla. Se quiere calcular la ley en X0. Los datos del
variograma se muestran en la siguiente tabla.
Muestra Ley (%) 8.2 S1

S1 8,2
S2 9,6
S3 13,1
S4 6,4 X0 9.6 S 2
Variograma 13.1 S 3
Tipo Esférico
Alcance a 250 m 6.4 S4
Efecto pepita C 0 17 100 m

Meseta C 66
224
γ(h) = C [(3/2)(h/a) – 1/2(h/a)3] + Co h ≤ a

C + Co h>a
γ(h) = 66 [(3/2)(h/250) - 1/2(h/250)3] + 17 h ≤ 250

66 + 17 = 83 h > 250
11 a1 + 12 a2 + 13 a3 + 14 a4 + = 01
21 a1 + 22 a2 + 23 a3 + 24 a4 + = 02
31 a1 + 32 a2 + 33 a3 + 34 a4 + = 03
41 a1 + 42 a2 + 43 a3 + 44 a4 + = 04
a1 + a2 + a3 + a4 + 0 = 1 (22)
113
225
La varianza σij entre las muestras Si y Sj viene dada por

el variograma γ(h), siendo h la distancia entre las
mismas.
8.2 S1
12 = 21 = 14 = 41 = ( 2002 + 100 2) =
2 2 3
200 + 100 ( 2002 + 100 2 (
X0 9.6 S 2 = 66 3/2 1/2 + 17 = 81,94
250 250 3
13.1 S 3
6.4 S4
100 m
13 = 31 = ( 1002 + 150 2 ) = 76,02 34 = 43 =

( 2002 + 50 2 ) = 80,13 ( h > 250)
= 01 = (100) = 54,48
23 = 32 = 30 = 03 = ( 1002 + 50 2) = 58,32 10
20 = 02 = (200) = 79,30
24 = 42 = ( 3002 + 1002) = 83,00 (h > 250)
40 = 04 =
( 1002 + 1002) = 67,02
226
Calculados los valores σij, se puede resolver el sistema

de ecuaciones, obteniendo los siguientes valores de ai:
a1 = 0,393
a2 = 0,022
a3 = 0,329
a4 = 0,256
Se introducen los valores de a y z en la siguiente

ecuación:
Z* = ∑ ai * zi
Gx = a1 S1 + a2 S2 + a3 S3 + a4 S4 =
= 0,393 (8,2) + 0,022 (9,6) + 0,329 (13,1)
+0,256 (6,4) = 9,38% de Pb
114
227
7. Cálculo de reservas
1) Datos de cada celda: F, P, L, d…
2) Se calcula el volumen.
Vb = Fb * Pb V = ∑ Vb
3) Se calculan las reservas de mineral.

Rb = db * Vb Rb = ∑ Rb
4) Se calculan las reservas de metal.

Ru = Rb * L / 100 Ru = ∑ Ru
228
CURVAS LEYES-
RESERVAS-RENTABILIDAD
115
229
Geometría del yacimiento ► ley de corte ► factores

económicos.
Cada ley de corte ► yacimiento de geometrías

distintas ► cantidad de reservas ► ley media.
La rentabilidad se calcula basándose en estos tres

parámetros.
Se debe calcular la rentabilidad a partir de distintas

leyes de corte.
230
Curvas leyes-reservas-rentabilidad
Establecen la relación entre cada ley de corte, con sus
reservas, calidad asociada y rentabilidad.
Reservas
+ +
+ + R
++
+ +
+ +
Rentabilidad Rmáx. Lm Ley media
+
+ +
Lc +
+
Ley de corte
116
231
En las fases avanzadas de investigación, se pueden

usar para afinar la selección de la parte explotable.
Cuando la ley de corte está muy definida, lo que se

tantea son variantes de explotación.
M
Rmáx.
+ C
B+
Rentabilidad +D
+ A + E
R Reservas
232
CUBICACION EN MINERIA
117
233
La definición de yacimiento no puede separarse de su

valoración y consideración como futura mina.
Por ello, las cubicaciones se deben realizar atendiendo

los tipos de minería.
La aplicación del método de explotación supone

programarlo en el tiempo y en el espacio.
Hay que considerar todos los factores.
Las reservas geológicas son distintas de las explotables.
234
CUBICACIONES EN MINERIA
A CIELO ABIERTO
El 70% de los minerales procede de minas a cielo
abierto.
• Alta productividad
• Mayor concentración de operaciones
• Mayor producción por explotación
• Menor inversión por tonelada
• Menos costes de extracción
• Posibilidad de yacimientos de baja ley
• Mejor conocimiento geológico
• Menor limitación en el tamaño y peso de los equipos
• Operaciones sencillas
• Mayor recuperación y menor dilución
118
235
1
Estéril 1
2 2
3
Estéril Talud
3
Mineral
Mineral
1
2
3
236
LOS RATIOS DE DESMONTE:

CRITERIOS ECONOMICOS
Ratio. Es el cociente entre la cantidad de estéril

arrancado y la de mineral extraído (t/t, m3/m3 o m3/t).
Connotaciones económicas: profundidad y extensión.
Tipos de ratio
• Ratio medio
• Ratio límite
119
237
Ratio medio. Es el cociente entre todo el estéril y todo

el mineral extraídos en un momento dado.
Rm = E / M
2
Estéril
3
Mineral
238
Ratio límite. Es la relación puntual máxima de estéril a

mineral a la que puede explotarse una unidad de este
para que proporcione un beneficio prefijado.
RI = E/M
120
239
Para introducir el componente económico se realiza la

siguiente operación:
Pv = ž * L * Pu = ∑CT+ B
∑CT
• Coste del arranque carga y transporte CE = E * Rm
E: coste unitario de extracción de m3 de estéril
Rm: ratio medio
• Resto de costes ∑C
Pv = ž * L * Pu = E * Rm + ∑C+ B
Rm tiene connotaciones económicas y geométricas.

Hay dos incógnitas Rm y L.
L = a * Rm + b
240
CALCULO DEL RATIO MEDIO
1. Se modeliza el yacimiento.
1
1.5
2
2.5
2. Se trazan perfiles geológicos.
121
241
CALCULO DEL RATIO MEDIO
3. Se hacen tanteos.
E
M
Rm = aceptable
Rm = E/M
242
Cuando se busca un ratio límite determinado…

1. Se modeliza geométricamente el yacimiento.
2. Se trazan los perfiles geológicos.
3. Sobre cada perfil, se tantean puntualmente ratios
límite.
• Se construyen mapas de isorratios.
• Se deducen las zonas más favorables.
122
243
EJEMPLO
Se ha cubicado un yacimiento aurífero a cielo abierto,

obteniendo los resultados de la tabla. Evaluar si la
configuración de mina cubicada (para beneficio nulo) se
ajusta a lo que resulta de aplicar los parámetros
económicos de la tabla.
Datos
Reservas de mineral 1,90 Mt
Ratio medio 6,52 t/t
Ley media 1,97 g/t
Estimaciones
Coste de estéril 2,42 $/t
Resto de costes 17,13 $/t
Rendimiento del concentrador 0,92
Precio de venta 3,71 $/onza
244
EJEMPLO
a) Cálculo del ratio medio de beneficio nulo
ž * L * Pu = E * Rm + ∑C + B
Datos
E 2,42
∑C 17,13 Rm = 1,94 t/t
ƞ 0,92
Pu 3,71/31,104=11,92
B 0
L 1,97
123
245
EJEMPLO
b) Cálculo de la ley crítica
ž * L * Pu = E * Rm + ∑C+ B
Datos
E 2,42
L = 2,96 g/t
∑C 17,13
ƞ 0,92
Pu 3,71/31,104=11,92
B 0
Rm 6,52 t/t
246
EJEMPLO
c) Rentabilidad de la explotación
Ingresos Costes
R 1,9 Esteril Otros
L 1,97 R 1,9 1,9
ƞ 0,92 Rm 6,52
Pv 371
C 2,42 17,13
(R*L* ƞ*Pv)/31,104
R*Rm*C R*C
41.520.296 $
$ 20.978.960 $ 32.547.000
Total
Diferencia
$ 62.525.960
21 M$
11 $/t
124
247
CUBICACIONES EN MINERIA
SUBTERRANEA
Este tipo de minería, no está afectada por el ratio.

En etapas iniciales, se cubica el yacimiento en conjunto.
El método de explotación es fundamental.
No se aprovecha todo el mineral.
Se deben aplicar coeficientes reductores.
Al aumentar la investigación, se definen zonas
explotables.
Existe el fenómeno de dilución.
Hay reducción de las reservas debido a las fracturas.
248
CASOS PRACTICOS
CON METODOS
TRADICIONALES
125
249
DISEÑO OPTIMO DE UNA MINA

DE CARBON A CIELO ABIERTO
250
Yacimiento: Sinclinal con dos capas explotables de

hulla.
Objetivo: Definir el diseño óptimo para programar su

explotación temporal.
Parámetro: Beneficio bruto global.
126
251
1. Cubicación de variantes
2. Determinación del precio de venta
3. Estimación de costes
4. Cálculo de la rentabilidad
252
CUBICACION DE VARIANTES
Sucesivas campañas de investigación
Modelización espacialmente de las capas de carbón
Realización de tanteos de hueco de explotación
Criterios de explotación
Variantes 1, 2 y 3; variantes 4 y 5; y variante 6
Taludes en función informe geotécnico
127
253
CUBICACION DE VARIANTES
Variante Reservas de carbón (t) Ratio medio (m3/t)
1 572.878 23,7
2 530.096 21,8
3 544.355 22,1
4 632.708 17,6
5 687.830 18,9
6 434.551 15,7
254
PRECIO DE VENTA
Fórmula oficial de precios para hulla y antracita
P = (8974/1000) * [1000 + 7(V – 20 + 20(25 – C)] * (88 – H/78)
• P: precio tonelada vendible (PTA/t)

• V: contenido en volátiles sobre muestra seca (%)
• C: contenido en cenizas sobre muestra seca (%)
• H: humedad (%)
128
255
PRECIO DE VENTA
Volátiles: 18,2%
Cenizas: 26,5%
Humedad: 5,4%
Azufre: 0,71%
Poder calorífico: 5722 kcal/kg
V = 18,2 * 100 / (100 – 5,4) = 19,2%

C = 26,5 * 100 / (100 – 5,4) = 28%
P = 8.800 PTA/t
256
ESTIMACION DE COSTES
COSTES
Arranque, carga y transporte estéril 1 m3 250 PTA/m3
Arranque, carga y transporte de carbón 1t 450 PTA/t
Transporte a central térmica 1t 650 PTA/t
Restauración 300 PTA/t
Amortización Coste de la operación
Variante 1, 2 y 3 221,55 MPTA
Resto
Variante 4, 5 y 6 257,55 MPTA
Propiedad 200 PTA/t
129
257
ESTIMACION DE COSTES
VARIANTE RESERVAS RATIO COSTE POR TONELADA PRECIO BENEFICIO

DE
CARBON MEDIO CONTRATA (PTA/t) PROPIEDAD TOTAL VENTA en GLOBAL
(t) (m 3 /t) (PTA/t)
(PTA / t) (PTA/t) PTA/t BRUTO
ESTERIL CARBON TRANSPOR. RESTAUR. AMORTIZACION TOTAL (MPTA)
1 572,878 23,7 5,925 450 650 300 387 7.712 200 7.912 8.880 968 554
2 530,096 21,8 5,450 450 650 300 418 7.268 200 7.468 8.880 1,412 748
3 544,355 22,1 5,525 450 650 300 407 7.332 200 7.532 8.880 1,348 743
4 632,708 17,6 4,400 450 650 300 407 6.207 200 6.407 8.880 2,473 1,565
5 687,830 18,9 4,725 450 650 300 374 6.499 200 6.699 8.880 2,181 1,500
6 434,551 15,7 3,925 450 650 300 593 5.918 200 6.118 8.880 2,762 1,200
258
CALCULO DE LA RENTABILIDAD
Bbg (MPTA)
1600
1500
1400
1300
1200
1100
R(kt)
420 460 500 540 580 620 660 700
Am (m3/t)
20
5
19 E
4
18 D
C
17 B
A
6
16
15
14
R(k1)
420 460 500 540 580 620 660 700
130
259
Bbg = Pv * R - ∑C * R
Bbg = Pv * R – (250 + 1600 + 257,55 * 106) * R
Beneficio bruto global

Variante Reservas de carbón (t) Ratio medio (m3/t)
(MPTA)
A 500.000 16,05 1.376
B 540.000 16,40 1.460
C 580.000 16,85 1.522
D 600.000 17,10 1.545
E 650.000 18,00 1.549
260
Bbg(MPTA)
1600 D 4 E
C
5
1500 B
A
1400
1300
6
1200
1100
R(k1)
420 460 500 540 580 620 660 700
131
261
EVALUACION DE RESERVAS
EXPLOTABLES A CIELO ABIERTO
EN UN YACIMIENTO
DE MAGNESITA
262
Yacimiento: Estratiforme, en un nivel de gran potencia y

cambios laterales muy acusados.
Buzamiento: 70º.
Composición: magnesita (CO3Mg), sílice (SiO2) y cal

(CaO).
Información
• Planos geológicos, perfiles
• 44 sondeos (41 horizontales)
• 10 800 contenido en sílice y cal
• Estudio geotécnico
132
263
1. Definición de la ley de corte
2. Definición de los criterios de selectividad
3. Definición del método de cubicación
4. Evaluación de reservas
264
DEFINICION DE LA LEY DE CORTE
Material para fabricación de ladrillos
Contenidos en cal y sílice penalizados

Cal: < 8%
Sílice < 6%
Mena: magnesita con % inferiores

Estéril: magnesita con % superiores rocas presentes
133
265
CRITERIOS DE SELECTIVIDAD
Potencia mínima explotable de mineral: 3 m
Potencia mínima separable de estéril: 3 m
Dilución
• Densidad magnesita = densidad de estéril
• Mineral bruto +7% sílice y +10% cal que el mineral
in situ.
266
METODO DE CUBICACION
Método de perfiles: 10
Perfil dividido en rebanadas horizontales de 15 m
Cada rebanada fue denominada bancada, y enumerada
por su cota inferior
Metodología
• Medida de superficies de estéril y mineral
• Cubicación de bloques limitados por perfiles
• Cálculo de leyes por bloque
• Cubicación total
134
267
1. Superficies de estéril y mineral
PERFIL 7 BANCADA 675
PERFIL 7 BANCADA 725 PERFIL 7 BANCADA 710 SUPERFICIE TRAMOS SUPERF. SUPERF. %CONTENIDO
SUPERFICIE
TOTAL
TRAMOS
DE
SUPERF.
MINERAL
SUPERF.
ESTERIL
%CONTENIDO SUPERFICIE
TOTAL
TRAMOS
DE
SUPERF.
MINERAL
SUPERF.
ESTERIL
%CONTENIDO
TOTAL DE MINERAL ESTERIL
(EST-MIN) MINERAL (m) (m) CaO SiO 2 (EST-MIN) MINERAL (m) (m) CaO SiO2
3770 1 30 60 6.50 5.50 3610 1 60 3550 6.50 5.50

(EST-MIN) MINERAL (m) (m) CaO SiO 2
2 120 130 6.50 3.00 2 130 3420 6.50 3.00
3 280 250 3 250
4 220 160 4 160 3010
5 880 1020 2.00 2.10 5 1020 1990 2.00 2.10
6 190 110 3.00 3.50 6 100 1900 3.00 3.50
7
8
9
5
0
0
0
0
0
0.00 0.00
7
8
9
0
0
0
1800 0.00 0.00 3400 1 60 3350 6.50 5.50
10
11
0
0
0
0 0.00 0.00
10
11
0
0
1800
0.00 0.00 2 110 3240 6.50 3.00
TOTAL 1725 2045 2,95 2,54 TOTAL 1810 1800 3,06 2,55
3 250 2990
PERFIL 7 BANCADA 675 PERFIL 7 BANCADA 680 4 200 2790
SUPERFICIE
TOTAL
(EST-MIN)
TRAMOS
DE
MINERAL
SUPERF.
MINERAL
SUPERF.
ESTERIL
(m)
%CONTENIDO SUPERFICIE
TOTAL
(EST-MIN)
TRAMOS
DE
MINERAL
SUPERF.
MINERAL
SUPERF.
ESTERIL
(m)
%CONTENIDO
5 2710 2.90 2.10
(m) CaO SiO 2 (m) CaO SiO 2
3400 1 30 3360 6.50 5.50 2910 1 30 2880 6.50 5.50

6 910 1800 2.90 2.10
2 120 3240 6.50 3.00 2 120 2750 6.50 3.00
3
4
280
220
2990
2790
3
4
200
210
2550
2340 7 180 1620
5 880 2710 2.90 2.10 5 2240 2.90 2.10
6
7
190
5
1800
1620
2.90 2.10 6
7
920
110
1320
1210
2.90 2.10
8 0 0.00 0.00
8 0 1620 0.00 0.00 8 0 0.00 0.00
9
10
11
0
0
0
1620
1620 0.00 0.00
9
10
11
0
0
0
1210
1210 0.00 0.00
9 0
TOTAL 1780 1620 3,01 2,52 TOTAL 1700 1210 3,06 2,42 10 0 1620
11 0 1620 0.00 0.00
TOTAL 1760 1620 3,01 2,52
268
2. Cubicación por bloques
BLOQUE - F
PERFIL - 6
LOTE
PERFIL - 6
BANCADA SUPERFICIE SUPERFICIE % CONTENIDO BANCADA SUPERFICIE SUPERFICIE % CONTENIDO CaO SO2
DE DE DE DE
ESTERIL MINERAL 2
ESTERIL MINERAL
(m2 ) (m2 ) CaO SO2 (m2 ) (m2 ) CaO SO2
SUP -35 0 0 0.00 0.00 - 0 0 0.00 0.00 -
935-920 0 0 0.00 0.00 - 0 0 0.00 0.00 -
SUP -35 0 0 0.00 0.03 SUP -35 0 0 0.00 0.00
920-906 0 0 0.00 0.00 - 0 0 0.00 0.00 -
035-920 0 0 0.00 0.03 935-920 0 0 0.00 0.00
920-905 0 0 0.00 0.03 920-906 905-890 62843 0 0.00 0.00 - 52943 0 0.00 0.00 -
0 0 0.00 0.00
906-830 580 0 0.00 0.00 905-890 810-875 71179 0 0.00 0.00 - 134421 0 0.00 0.00 -
400 0 0.00 0.00
830-875 480 0 0.00 0.00 810-875 875-880 73744 0 0.00 0.00 - 207765 0 0.00 0.00 -
650 0 0.00 0.00
875-880 470 0 0.00 0.00 875-880 680 0 0.00 0.00 880-845 85285 0 0.00 0.00 - 293951 0 0.00 0.00 -
880-845 720 0 0.00 0.00 880-845 610 0 0.00 0.00 845-830 83775 10901 5,21 4,76 8,24 382826 10901 5,21 4,76 35,12
845-830 740 170 5.21 4.76 845-830 660 0 0.00 0.00 830-815 83383 22045 5,28 4,18 3,61 468189 33985 5,26 4,37 13,72
830-815 640 330 5.31 4.24 830-815 660 30 5.00 3.50 815-800 79838 54827 4,90 3,50 1,46 546024 88812 5,03 3,45 6,15
815-830 620 640 4,94 3,65 815-800 625 245 4,78 3,20 800-785 97505 105135 4,80 2,91 0,90 643929 193977 4,85 3,34 3,30
830-785 740 100 5,04 3,01 800-785 740 740 4,27 2,73
785-770 74054 123751 4,04 3,04 0,60 714994 317738 4,54 3,22 2,25
785-770 615 110 4,42 2,94 785-770 540 1320 3,71 3,12
770-755 97791 134983 3,83 3,05 0,72 812784 452721 4,32 3,17 1,80
770-755 935 115 4,04 2,82 770-755 590 1190 3,68 3,22
755-740 196543 160313 3,42 2,52 1,23 1009328 613033 4,09 3,00 1,65
755-740 2075 180 4,16 2,74 755-740 190 1520 2,94 2,33
740-725 294854 186810 3,29 2,55 1,55 1303982 832843 3,90 2,90 1,62
740-725 2550 1235 3,74 2,53 740-725 2345 1725 2,94 2,54
725-710 1460 1500 3,76 2,43 725-710 1810
725-710 209058 212254 3,38 2,50 0,96 1513029 1015096 3,79 2,41 1,49
1800 3,06 2,55
710-695 1230 1560 3,84 2,33 710-695 1620 1780 3,01 2,52
710-695 182756 216101 3,45 2,41 0,85 1695784 1231197 3,72 2,75 1,38
695-680 970 1260 3,77 2,43 695-680 1210 1700 3,05 2,42 695-680 139793 189810 3,36 2,44 0,74 1835578 1421007 3,67 2,71 1,29
680-665 742 1068 4,17 2,97 680-665 1500 1850 3,01 2,37 680-665 111706 174292 3,46 3,61 0,64 1947284 1595256 3,65 2,70 1,22
680-650 30 260 4,35 3,58 665-650 750 1480 3,00 2,25 665-650 50018 110295 3,29 2,51 0,45 1997301 1706693 3,62 2,68 1,17
650630 0 0 0,00 0,00 660-630 790 1820 2,74 2,19 660-630 50859 110706 2,74 2,19 0,43 2047960 1822300 3,56 2,65 1,12
TOTAL 15577 11728 4,15 2,81 TOTAL 16360 16890 3,15 2,54
TOTAL 2047960 1822304 3,58 2,85 1,12 2047960 1822300 3,56 2,65 1,12
135
269
2. Cubicación por bloques
P-11
920 905 890 875

d1
920
d2 P-11
905
d3
d3
d2
890
d1
d4
878
860
PERFIL PLANTA
270
EVALUACION DE RESERVAS
Los tonelajes de mineral y estéril se evaluaron
sumando los tonelajes de los bloques.
Las calidades resultaron de la ponderación de los
valores de los bloques con sus tonelajes individuales

Metodos de Geoestadistica PDF

Cargado por

Copyright:

Formatos disponibles

Metodos de Geoestadistica PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Metodos de Geoestadistica PDF

Cargado por

Copyright:

Formatos disponibles

98

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

¿Cómo poder relacionar los valores con sus

¿Cómo relacionar dichos valores entre sí?

Teoría de la variable regionalizada

La variable toma un valor en cada punto del espacio.

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

z+ VALOR MEDIO ESTIMADO Z

VARIOGRAMA Z = ESTIMADOR LINEAL

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

El variograma liga el valor de una variable con su posición

γ (h) = Σ [Z(xi) - Z(xi + h)]2 / 2n

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Pero ¿cómo se construye un variograma?

+(37 - 37)2 + (37 - 33)2 + (33 - 34)2 + (35 - 38)2

42 43 42 39 39 41 40 38 +(39 - 36) 2 + (42 - 43) 2 + (43 - 39) 2 + (42 - 39) 2

γ (h) = Σ [Z(xi) - Z(xi + h)]2 / 2n

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

El crecimiento de la curva indica el grado de

22.5 45 67.5 90 112.5 135 157.5 m

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Los variogramas en varias direcciones ponen de

Variabilidad (Oz Ag/ ton) 2

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Esférico Gaussiano Exponencial Potencial (a<1)

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

γ(h) = C [(3/2)(h/a) – 1/2(h/a)3] + Co h ≤ a

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Estimar valores desconocidos a

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

El valor de la variable se calcula como una combinación

• Z*: valor de la variable a estimar

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

El krigeado es el mejor estimador de ai.

Varianza mínima: Var [Z* - Z] = σk2

VAR(Z* -Z) = σk2 = σv2 - 2 Σ ai σv Xi + Σ Σ ai aj σ Xi Xj

• σv2: varianza del bloque del volumen V

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Las varianzas y covarianzas son datos obtenidos del

Las incógnitas son los coeficientes de ponderación.

Se introducen en la siguiente ecuación:

La aplicabilidad es muy amplia.

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

¿Cómo determinamos la ley en el punto A?

El variograma determina qué muestras puede tener una

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Los cálculos son realizados por programas

Las opciones relativamente económicas son

Puede realizar una buena modelización de

Se pueden obtener mapas, donde la interpolación de

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade

Dra. Natalia Caparrini Marín - [email protected] - Consultora Intercade