separata unsa 2011
TRANSCRIPT
SEPARATA DEL CURSO DE GEOESTADÍSTICA APLICADAAREQUIPA 2011
Profesor: Reynaldo Canahua Loza
1. INTRODUCCIÓN
La Geoestadística aporta una herramienta denominada Variograma, el cual tiene aplicaciones en las ciencias geológicas y múltiples aplicaciones en otras disciplinas. La Geoestadística fue desarrollada y presentada por George Matheron en la década de los 60 del siglo XX. En efecto la Geoestadística considera que las variables regionalizadas están modeladas en un espacio de
variables aleatorias reales L2 sobre un espacio de probabilidades
Representación de la variable regionalizada en el espacio.En la que se define la función Variograma:
o
1
Y una extensión del Variograma denominado Variograma Cruzado que relaciona dos variables:
En esta pequeña divulgación el objetivo es presentar el Variograma y su extensión Variograma Cruzado mediante algoritmos aritméticos simples, a través de ejemplos, haciendo ver su aporte como herramienta de trabajo en la geología de minas, de petróleo, etc. Aporte que consiste en dar cuenta del aspecto estructural del fenómeno estudiado, aspecto estructural no contemplado por la estadística descriptiva.
2. APLICACIÓN EN LA DIFERENCIACIÓN DE DOS FENÓMENOS
2.1 En una línea de muestreo de la zona A, tenemos los siguientes valores de la
variable regionalizada de plomo en ppm.
Realizamos un análisis estadístico básico.
a) Media aritmética:
b) La varianza:
c) El coeficiente de variación:
2
d) Histograma
2.2 En otra línea de muestreo en la zona B, tenemos los mismos valores de la
variable regionalizada de plomo en ppm, pero dispuesto de la siguiente forma; es
decir, un fenómeno estructuralmente muy diferente, a pesar de tener los mismos
valores de leyes.
Obtenemos la media aritmética, la varianza, el coeficiente de variación y el
histograma, y vemos que da los mismos resultados que los obtenidos en la Zona A.
Es decir que con esta estadística descriptiva no logramos diferenciar dos
fenómenos totalmente diferentes.
2.3 Ahora procedemos a construir los Variogramas de la Zona A y B
3
Graficando el variograma para la zona A.
Graficando el variograma para la zona B.
4
)(h
h
Como se puede observar el semi-Variograma, que más comúnmente se le
denomina variograma, da cuenta de las zonas estructuralmente diferentes.
3. MODELOS DE VARIOGRAMAS MÁS COMUNES
a) Efecto de Pepita Puro
b) Modelo Esférico o de Matherón
5
)(h
h
c) Modelo de Formery o Exponencial
d) Modelo con efecto “HOLE”
e) Modelo Gaussiano
Es de notar que estos modelos y sus combinaciones no son necesariamente los
únicos.
6
4. APLICACIÓN DEL VARIOGRAMA CRUZADO EN EL ESTUDIO DE
CORRELACIONES
4.1 Para visualizar esta función, consideremos un ejemplo de una zona de terreno
explorado del cual hemos obtenidos valores geoquímicos del Oro y Plata y deseamos
estudiar la relación entre los dos valores. Para este efecto aplicaremos el variograma
cruzado.
Siendo la formula:
Aplicando:
Cuya gráfica es la siguiente:
7
h
)(hAuAg
Observamos que cuando hay una correlación positiva alta entre las variables, el
variograma cruzado tiende a tomar valores positivos altos.
4.2 Ahora veamos que pasa en otra zona, donde los valores geoquímicos de la Plata
toman otros valores:
Aplicando la fórmula:
Con su gráfica:
Observamos que cuando hay una correlación negativa alta entre las variables
el variograma cruzado tiende a tomar valores negativos altos.
Al realizar el análisis estadístico descriptivo de las líneas de muestreo de la Zona A y la
Zona B, correspondientes a la variable regionalizada plomo; obtenemos el mismo
resultado estadístico; es decir, que este análisis estadístico descriptivo no nos ayuda a
8
h
)(hAuAg
dar cuenta del aspecto estructural de estos fenómenos totalmente diferentes. En tanto, al
aplicar la función Variograma vemos que sí da cuenta del aspecto estructural de estos
fenómenos totalmente diferentes. Por lo que constituye una herramienta de gran utilidad
en el tratamiento de variables regionalizadas.
Observamos que cuando hay una correlación positiva alta entre las variables
regionalizadas, el variograma cruzado tiende a tomar valores positivos altos; mientras
que, cuando hay una correlación negativa alta entre las variables, el variograma cruzado
tiende a tomar valores negativos altos. Por lo que el variograma cruzado se suma a las
herramientas que permiten estudiar la correlación entre variables.
5. TÉCNICA DEL KRIGING DE MATHERON
Considerando las variables Z(xi) que están cerca de un soporte geométrico a
estimar y dentro de su aureola de influencia. Aureola definida por medio de los
alcances estimados a partir de un estudio de variogramas.
Visualizaremos el procedimiento a partir de la siguiente disposición de las
muestras con respecto a un soporte geométrico V.
9
)( 11 xz
)( 22 xz
)( nn xzV
Disposición de muestras con respecto a un soporte geométrico V.
Y deseamos estimar la variable Z(x) del soporte geométrico V a partir de las
muestras Z(xi). Entonces necesitamos encontrar los pesos para estimar la variable
Zv(x0), a partir de:
En el ejemplo que presentaremos a continuación consideremos que estamos en
condiciones de aplicar un Krigeage ordinario bajo la hipótesis estacionaria de
orden 2, por lo que usaremos el siguiente sistema de ecuaciones.
Este sistema resulta de minimizar la varianza de estimación sujeta a la condición
de universalidad , que hace que nuestro estimador sea insesgado.
El error cometido en este procedimiento de estimación viene dado por la varianza
de Kriging de Matheron siguiente:
Ejemplo:
A partir de los valores de la potencia de un manto de Hematita en los puntos A y
B, se desea estimar la potencia en el punto C. Considerando que la potencia del
manto tiene el siguiente modelo de variograma.
10
n
Vvuvv1
,,
1 a n
n
1
1
h
3
3
21
23
ah
ah
C
C ah
ah ,0
ma 6 y 1C
Es decir un modelo de Matheron o también denominado esférico. A continuación
se muestra la ubicación de las potencias y su orden de magnitud.
Variable regionalizada Potencia en los puntos A y B de un manto de Hematita
Entonces tendremos a partir del sistema de ecuaciones anterior, el siguiente
sistema particular:
Reemplazando:
ó
I.
11
II.
III.
Restando (I) – (II) :
Ahora tenemos el siguiente sistema:
I'.
II'.
Sumando (I’) + (II’) :
Reemplazando en (II’) :
En la primera ecuación inicial:
Despejamos:
12
Parámetro auxiliar que será usado posteriormente en la fórmula de la varianza de Kriging de Matheron.
Siendo la potencia estimada del manto igual a:
Reemplazamos:
Ahora veamos cuál es el error que se comete en esta estimación, para lo cual particularizamos la fórmula de la varianza de Kriging de Matheron dada anteriormente.
Reemplazando:
Este valor es el error cometido en el proceso de estimación realizado.
13
REFERENCIAS
- Matheron G. (1962, 1963) - Traité de Géostatistique Appliquée. Ed. Technip, Paris VOL. 1; VOL. 2.
- Guibal D. (1972) - Simulation de Schémas Intrinsèques. N-291 E.N.S.M.P.
- Journel A. (1977) - Géostatistique Miniere, tomo 1 y 2. E.N.S.M.P
- Maréchal A., Deraisme J., Journel A., Matheron G. (1978) - Cours de Géostatistique non Linéaire. C-74 E.N.S.M.P.
- Marín Suárez A. (1978) - Méthodologie de L'estimation et Simulation Multivariable des Grands Gisements Tridimensionnels. Thèse présentée à I'école Nationale Supérieure des Mines de ParisPara obtener el grado de:Docteur Ingénieur en Sciences et Techniques Minières - Option Géostatistique
- Marín Suárez A. (1986) - Modelo Geoestadístico de Filones de Almadén. Ed. Minas de Almadén S.A., Almadén (España).
- Remy N., Boucher A., Wu J., Journel A. (2009) - Applied Geostatistics with GSLIB and SGEMS. Ed. Cambridge University Press.
14