Krigin Puntual y Bloques
Krigin Puntual y Bloques
Krigin Puntual y Bloques
FACULTAD DE INGENIERA
ESCUELA ACADMICO PROFESIONAL DE INGENIERA DE MINAS
CURSO : Geoestadstica.
INTEGRANTES :
CASTILLO MORENO, Keylor.
MOSQUEIRA VILLAR, Paola
MUJICA CABRERA, Juan
PERALTA RODAS, Pool.
TERAN TOLEDO, Fernando.
20 DE JULIO DE 2016
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
NDICE
1 INTRODUCCIN......................................................................................................................3
2 OBJETIVOS ............................................................................................................................ 4
8 CONCLUSIONES ................................................................................................................... 37
9 REFERENCIAS BIBLIOGRFICAS............................................................................................ 38
2
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
1 Introduccin
La geoestadstica ofrece un mtodo de estimacin de reservas, usando el variograma,
llamado krigeado, se utiliza en la evaluacin de yacimientos para estimar el valor de una
variable regionalizada, en un punto o en un bloque, a partir del uso de factores de
ponderacin.
Se considera al mtodo de kriging del tipo MELI (Mejor Estimador Lineal Insesgado)
o ELIO (Estimador Lineal Insesgado ptimo): es lineal porque sus estimaciones son
combinaciones lineales ponderadas de los datos existentes; y es insesgado porque procura
que la media de los errores sea nula; es el mejor porque los errores de estimacin tienen
una variancia mnima.
3
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
2 Objetivos
4
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
3 Marco Terico.
Lineales:
- Simple.
- Ordinario.
- Universal.
- Residual
5
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
No lineales:
- Disyuntivo
- Indicador
- Probabilstico
Puntual.
En bloques.
Paramtrico:
- Multigaussiano.
- Disyuntivo.
- Lognormal.
No paramtrico:
- Simple
- Ordinario
- Universal
- Residual
- Indicador
- Probabilstico.
6
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
El objetivo del Kriging es interpolar la ley a partir de muestras que lo rodean, asignando
pesos a cada muestra que reduzca al mnimo la varianza de estimacin, se interpreta la ley
con la mxima precisin posible.
Imagen N 1. Kriging
4 Kriging puntual.
El krigeage puntual o krigeage ordinario; a diferencia del kriageage simple que asume
que las medias locales son relativamente constantes y de valor muy semejante a la media
de la poblacin que es conocida, la media de la poblacin es utilizada para cada
7
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
estimacin local, en conjunto con los puntos vecinos establecidos como necesarios para
la estimacin. En el Kriageage puntual las medias locales no son necesariamente
prximas de la media de la poblacin, usndose apenas los puntos vecinos para la
estimacin. Es el mtodo ms ampliamente utilizado en los problemas ambientales.
Es un interpolador exacto en el sentido de que las ecuaciones calculadas las utiliza para
interpretar los valores y habr una coincidencia exacta entre los valores interpolados y
los puntos de datos originales.
El mayor inters que presenta el krigeage puntual es que permite obtener una serie de
valores estimados a los cuales, posteriormente, se les puede ajustar una serie de lneas
de isocontenidos, por ejemplo para las leyes o potencias del nivel mineralizado.
En algunas ocasiones, en vez de estimar la ley media de un bloque V, interesa estimar la
ley en un punto x0 (problema de interpolacin).
8
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Imagen N 2. Krigeado del punto x0. Se puede generar una grilla de valores interpolados al hacer variar x0. Esta
tcnica tiene aplicacin en la cartografa automtica y en la simulacin de leyes.
5 Ejercicio resuelto.
El cuadrillado representa sondajes en los puntos X1, X2, X3, X4, con leyes de 3,2%; 2,5%;
4,0% y 1,5% respectivamente. Se desea conocer la ley en el punto X0, sabiendo que el
variograma asociado corresponde a un modelo esfrico con las siguientes caractersticas:
Co = 2, C = 20 y a = 200 m.
9
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Solucin
El problema se resuelve de la siguiente manera:
[ ] = []1 [2]
10
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Se tiene que:
11
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
[ ] = []1 [2]
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Al resolver se obtiene:
1 = 0.518
2 = 0.022
3 = 0.089
4 = 0.371
= 0.915
Se verifica que:
d) Calcular la ley.
6 Kriging en bloques.
El mtodo ms usado en la modelacin de recursos consiste en la discretizacin del espacio
3D en bloques o celdas tridimensionales. Cada celda contiene los atributos (litologa, tipo
de mineralizacin etc.) y las mediciones (leyes, propiedades fsico mecnicas) del dominio
geolgico en que se encuentra. Los atributos de los bloques se determinan sobre la base
de la interseccin con el modelo geolgico o su posicin respecto a una superficie
triangulada y las leyes a travs de la estimacin con tcnicas de interpolacin espacial.
13
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
2.8 = (90 ; 0)
El primer modelo de bloque fue utilizado a comienzos de los aos 60 por la Kennecott
Koper Corporation en un depsito de prfido cuprfero. Se emple para describir la
distribucin espacial de las leyes y no la geometra de los dominios geolgicos.
Cada bloque debe contener toda la informacin disponible en las fases de desarrollo de
un proyecto: litologa-mineraloga, contenidos de metales, calidades en el caso del carbn
y rocas industriales, contenidos de contaminantes, parmetros geomecnicos, datos
hidrogeolgicos, etc.
14
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Posicin del modelo: se especifica a partir de las coordenadas del centroide del
bloque llave (key block).
Con el objetivo de alcanzar una mayor resolucin del modelo de bloque en los lmites de
los cuerpos minerales se utilizan bloques (sub bloques) con dimensiones menores que los
originales. El modelo de bloques puede ser rotado y orientado de manera que se ajuste a la
estructura geolgica y respete los elementos de yacencia del yacimiento estudiado.
15
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Imagen N 5. Parmetros que definen el modelo de bloque: Posicin del bloque llave (Xmorg, Ymorg, Zmorg),
dimensiones del bloque (dx, dy, dz), extensin del modelo (X (nx), Y (ny), Z (nz).
Este enfoque terico en muchos casos no es prctico desde el punto de vista tcnico
(demasiados sub-bloques para poder respetar los lmites del modelo geolgico y lograr una
16
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
17
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
() = ()
=1
18
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
La varianza del error de prediccin del kriging en bloques est dada por:
Isaaks y Srivastava (1989) muestran a travs de ejemplos que el kriging en bloques coincide
con el promedio de predicciones hechas por kriging ordinario sobre cada uno de los puntos
del enmallado dentro del bloque. As mismo indican que en la prctica es suficiente con un
enmallado cuadrado (6x6) para obtener estimaciones estables en los bloques.
19
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
2K = 0,1888 2K = 0,2162
0,25 0,265
0,25
0,129
50 50
2K = 0,1668 2K = 0,2248
0,247 0,074
50 0,233 0,233 50
50 0,426 0,426
0,174
0,080 0,074
0,033
50 50
20
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
2K = 0,1668 2K = 0,2107
0,215 0,3437
0,111 0,016
50 0,242 0,106 50 0,2674 0,013
0,111 0,016
0,215
0,3437
50 50 100 150
2 2
K = 0,0827 K = 0,1206
0,208 0,1044
50 50
0,042 0,1456
50 50
(h) = 0,2 + 0,8 Sph(100) (h) = 0,7 + 0,3 Sph(100)
Interpolacin exacta: la estimacin en un sitio con dato es igual al valor del dato y la
varianza de kriging en este sitio vale 0
21
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Al tener sesgo condicional, se incurre en una mala apreciacin del negocio. La ley media
del material mandado a planta (material cuya estimacin supera una ley de corte) es
inferior a la ley media estimada de este material, mientras que la ley media del material
mandado a botadero es superior a la ley media estimada de este material.
22
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Vecindad mvil: se usa slo los datos cercanos al sitio (bloque) a estimar
Para validar los parmetros del kriging (modelo de variograma, vecindad elegida), se puede
usar los siguientes mtodos:
23
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Jack-knife: se divide la muestra inicial en dos partes (por ejemplo, cuando hay dos
campaas de sondajes), y se estima una parte a partir de la otra
Luego, se hace un estudio estadstico de los errores cometidos para saber si el kriging fue
satisfactorio (buena precisin, poco sesgo condicional)
Criterios de validacin:
Medias de los errores y de los errores estandarizados: deben ser cercanas a cero
estimador sin sesgo
7 Ejercicio resuelto.
1. Mediante el mtodo del Kriging Ordinario, se tienen las ecuaciones que debemos
resolver:
24
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Resolvemos para V:
25
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
26
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
27
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
2. Considere una regionalizacin de dos dimensiones que se caracteriza por una funcin
aleatoria de punto () que, como una primera aproximacin, se supone que es intrnseca
con un semivariograma isotrpico estacionario: () = (||).Por ejemplo, Z (u, v) puede
ser el grosor vertical de una veta sedimentaria.
28
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
= ( )
=1
Con
1
(2 ) = [(4 )+(5 )]
2
29
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Donde H (L, l) denota una funcin auxiliar, que da la media del semi-variograma () si un
extremo del vector se fija en una esquina de un rectangular de tamao Lxl y el otro
extremo se mueve uniformemente por el rea de la rectangular. Esta funcin se puede dar
analticamente para variogramas simples, los valores ms complicados pueden ser ledos
desde diagramas grficos.
Donde (, ) es una funcin auxiliar similar H; que representa la media de () si los dos
extremos de se mueven uniformemente a travs del rea del rectangular de tamao
() .
30
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
31
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
32
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Lo que da
Podemos ver que en este caso de efecto de pepita puro, los pesos son
proporcionales al tamao del soporte. La ubicacin no es importante.
33
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
34
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
1 + 2 + 3 = 1
Con la solucin:
24.
35
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
Para este propsito, es necesario caracterizar en alguna forma la variabilidad espacial del
fenmeno en estudio. El mtodo geoestadstico es utilizar la funcin estructural () con
la que se puede calcular la varianza de estimacin de cualquier estimador lineal insesgado.
.344 .182
.172 .091
.323 .553
.344 .182
.172 .091
.296 .339
.344 .182
.172 .091
.27 .225
36
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
8 Conclusiones
El krigeage puntual a diferencia del kriageage simple que asume que las medias
locales son relativamente constantes y de valor muy semejante a la media de la
poblacin que es conocida.
37
UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERA
Escuela Acadmico Profesional de Ingeniera de Minas
9 Referencias bibliogrficas
Proyecto de Evaluacin de Yacimientos - Clase N 6
DE MXICO.
http://geoestadistica.com/objetivo_kriging.htm.
38