trabajo final kriging ordinario
TRANSCRIPT
1UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
2UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
3UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
4UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
5UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
6UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
7UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
8UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
9UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
10UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
11UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
12UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
13UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
14UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
15UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
16UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
17UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
18UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
4. PARTE PRACTICA
APLICACIÓN DEL KRIGING ORDINARIO
1) iniciamos el programa y se aprecia el siguiente entorno
19UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
2) Luego damos click en Import File y escogemos el archivo oro.txt
3) Se aperturara la siguiente ventana y le damos click en ok, para modificar los atributos de cada columna
20UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
4) Se renombraran las columnas: nro de muestra, coordenada x, coordenada y z
21UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
5) Los datos ya se mostraran ordenados
22UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
6) ahora procedemos a realizar el análisis estadístico unidimensional mediante el menú Data/Summary stadistics/ Primate variables
23UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
7) Aparecerá la siguiente ventana y hacemos click en x,y, coordinates
8) Y veremos el plano de ubicación de las muestras con las respectivas leyes
24UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
9) Dando click en “Z” y se apreciara la distribución de frecuencias
0 2 5 7 9
0.04 0.78 1.51 2.25
Fre
quen
cy
Au oz/TM
NonTransformed
0.0
25.0
50.0
75.0
100.0
0.04 0.78 1.51 2.25
Cum
ulat
ive
Pct
Au oz/TM
NonTransformed
Y las demás medidas estadísticas
25UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
10) Una vez realizado el análisis, procedemos a calcular el respectivo variograma, mediante el menú Autocorrelation/Variogram/Primate variable
26UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
11) Le damos click en calcúlate, y se calcularan los siguientes variogramas:
0.000
0.026
0.053
0.079
0.106
0.00 222.10 444.21 666.31
Sem
ivar
ianc
e
Separation Distance (h)
Au oz/TM: Isotropic Variogram
Spherical model (Co = 0.00600; Co + C = 0.08700; Ao = 42.00; r2 = 0.152; RSS = 7.204E-04)
0.000
0.030
0.061
0.091
0.121
0.00 222.10 444.21 666.31
Sem
ivar
ianc
e
Separation Distance (h)
Au oz/TM: Anisotropic Variogram (0º)
Linear model (Co = 0.08700; Co + C = 0.20933; AMajor = 1156100.00; AMinor = 1156000.00; r2 = 0.053; RSS = 0.0122)
27UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
0.000
0.030
0.061
0.091
0.121
0.00 222.10 444.21 666.31
Sem
ivar
ianc
e
Separation Distance (h)
Au oz/TM: Anisotropic Variogram (45º)
Linear model (Co = 0.08700; Co + C = 0.20933; AMajor = 1156100.00; AMinor = 1156000.00; r2 = 0.053; RSS = 0.0122)
0.000
0.030
0.061
0.091
0.121
0.00 222.10 444.21 666.31
Sem
ivar
ianc
e
Separation Distance (h)
Au oz/TM: Anisotropic Variogram (90º)
Linear model (Co = 0.08700; Co + C = 0.20933; AMajor = 1156100.00; AMinor = 1156000.00; r2 = 0.053; RSS = 0.0122)
28UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
0.000
0.030
0.061
0.091
0.121
0.00 222.10 444.21 666.31
Sem
ivar
ianc
e
Separation Distance (h)
Au oz/TM: Anisotropic Variogram (135º)
Linear model (Co = 0.08700; Co + C = 0.20933; AMajor = 1156100.00; AMinor = 1156000.00; r2 = 0.053; RSS = 0.0122)
Los 4 ultimos variogramas son modelados con el modelo linear.
12) Procedemos a la estimación mediante el kriging
29UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
Escogemos Kriging ordinate y marcamos point kriging. Le damos click en calcúlate
El proceso de estimación tarda unos minutos de acuerdo a la cantidad de muestras.
30UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
Los resultados en 2D y 3d respectivamente
31UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
ORO; Gilkey field
Au oz/TM
1.000.870.730.600.470.330.200.07-0.07-0.20-0.33-0.47-0.60-0.73-0.87-1.00
32UNIVERSIDAD MACIONAL JORGE BASADRE GROHMANN Facultad de Ingeniería de MinasGeoestadistica
5. CONCLUSIONES
El kriging ordinario establece una condición adicional al sistema de ecuaciones del kriging simple para filtrar el valor desconocido de la media.
Generalmente el valor de la media m es desconocido y por lo tanto no se puede utilizar el kriging simple.
El paquete GS+ es una herramienta fuerte para el cálculo de vario gramas, kriging, cokriging y simulaciones.
Con los menús interactivos, es posible elegir el modelo de variograma y ver el de mayor correlación
El paquete tiene la opción aparte de calcular el kriging puntual, calcula el kriging por bloques, dándoles sus respectivas medidas y graficar los mapas de leyes de acuerdo a la ley cutt off ingresada.
La interpolación se puede realizar por otros métodos: como el cokriging y el inverso a la distancia, todos ellos se pueden realizar con el GS+.
Para la validación del kriging se uso el método de la validación cruzada, el cual estima sucesivamente cada dato considerando solamente los datos restantes.