Geoestadística III

GLG-213

Universidad Mayor de San Andrés

 Variogramas
• Ya se vio que la geoestadística consiste en la
dependencia espacial entre las muestras
• El variograma expresa básicamente este
postulado
• En la práctica, describe el promedio de las
diferencias al cuadrado de los valores que son
calculados en función de la distancia que los
separa

Universidad Mayor de San Andrés

Variogramas

Universidad Mayor de San Andrés

Variogramas

Universidad Mayor de San Andrés

Variogramas

Universidad Mayor de San Andrés

 Variogramas
n

 x  x 
1
 (h)  i ih
2

2n( h ) i 1

xi son los valores (leyes)

xi+h es el valor del dato a una distancia h de xi
n(h) es el número total de pares de valores

Universidad Mayor de San Andrés

 Variogramas
n

 x  x 
1
 (h)  i ih
2

2n( h ) i 1

•Variogramas son ideales para identificar características anisotrópicas en depósitos

minerales asi como cualquier dependencia espacial
•En la primera fase variogramas son calculados en una dirección, en una segunda
etapa son calculados en dirección perpendicular

Universidad Mayor de San Andrés

Variogramas
200 m

400 m
Universidad Mayor de San Andrés
Variogramas

Universidad Mayor de San Andrés

 Variogramas
No. de diferencia xi xi+h xi-xi+h (xi-xi+h)2
∆1 60 85 -25 625
∆2 85 105 -20 400
∆3 60 75 -15 225
∆4 85 92 -7 49
∆5 75 92 -17 289
1 n
 ( 200)  
 i ihx  x 2 1588
2  5 i 1
Sin embargo se necesitan al menos 5 valores
1 de diferencias para obtener un variograma
 ( 200)  1588  158 .8
25 significativo
 ( 400) 
1
60  105 2  105  120 2  562 .5
2 2

Universidad Mayor de San Andrés

Variograma transitivo

Universidad Mayor de San Andrés

  h  C0  Ch
Variograma

Variograma

Covarianza

Distancia
Para funciones aleatorias estacionarias se tiene que
Universidad Mayor de San Andrés
 Variograma transitivo
La varianza de una suma o diferencia de variables puede
ser generada de la misma manera que una ecuación cuadrática
Para la varianza de la diferencia entre xi y xi+h,  xi  xih  deriva de:
2

 2 x  xi ih    2 xi   2 xih  2 xi   xi  h

Debido a que xi y xi+h derivan de la mima población tienen la misma varianza

 2 x  x
i ih   2 2 xi  2Covxi , xi  h 

 xi  xih    2 xi  Covxi , xi  h 
1 2
2
Esto es lo mismo que  (h ) por lo tanto:

 2 x  x
i ih   2 2 xi  2Covxi , xi  h 
La covarianza es 0 si xi y xi+h son
 xi  xih    ( h )   2 xi  Covxi , xi  h 
1 2 variables independientes y la distancia
2 h es mayor que el rango a.
 ( h )   2 xi
Universidad Mayor de San Andrés
 Variograma transitivo

Rango a
Efecto pepita Co
Meseta C1
Co+C1 el valor total de la meseta
La tangente de la curva intercepta la extensión de la meseta a 2/3 del rango a

Universidad Mayor de San Andrés

 Variograma transitivo
• Cada deposito mineral tendrá un variograma
típico
• El efecto pepita corresponde a distancias
prácticamente igual a cero
• Depósitos muy irregular como aurífero o
pegmatita tendrá un efecto pepita y/o un
pequeño rango
• Depósitos uniformes (de carbón, Pb-Zn
sedimentarios) tienen 0 efecto pepita y un gran
rango “a”

Universidad Mayor de San Andrés

 Variogramas
valor acumulado GT ley G (en variograma gradado
Depósito espesor
(ley G *espesor T) espesor constante)

a Co/C a Co/C a Co/C

15 m (buzamiento abajo)
Niquel magmático Australia 30 m 0.14 0.3
30 (a lo largo del rumbo)
Oro en vetas de cuarzo 40-60 m (buzamiento abajo)
0.3
Australia 20-30 (a lo largo del rumbo)
Pb, Zn en carbonatos 100 m (a lo largo del rumbo) 0.4
10
Canada 20-30 (perpendicular al rumbo) 0.33
Pb, Zn en vetas
10 m 0.4
Alemania
Fosfatos, Israel 700 m 0 600 m 0
Cu volcanogenico
400 m 0.35 300 m 0.35
Portugal
Cu/Zn volcanogenico
110 m 0.18 70 m 0.2
Turquia
50-90 m 0.01
Bauxita, Hungria
30-50 m 0.15

Universidad Mayor de San Andrés

Co/C1 vs. rango a

Universidad Mayor de San Andrés

Ejemplos de correcta determinación
de h para el variograma

Universidad Mayor de San Andrés

 Calculando variogramas experimentales direcionales

Puntos en 2-D o 3-D regular o irregularmente espaciados

Especificación de dirección (irregular): Especificación de dirección (regular):

Universidad Mayor de San Andrés

Calculando variogramas experimentales direccionales

Ejemplo: Empezamos con una distancia lag (por ej. #4)

Empezar en un punto y compara

los valores con todos los demás
puntos que caen en la distancia o
(lag) y el ángulo de tolerancia

Universidad Mayor de San Andrés

Calculando variogramas experimentales direccionales

Nos movemos al siguiente punto

Universidad Mayor de San Andrés

 Calculando variogramas experimentales direccionales

Ahora repetimos para todos los puntos

y repetimos para todas las lags

Sin correlación
Variograma, γ (h)

Se incrementa
la variabilidad

Distancias o lags (h)

Universidad Mayor de San Andrés

 Calculando variogramas experimentales direccionales

Selección de puntos P1 y P2 Vector (P1, P2, 45°) que

en un dominio espacial cae entre 30° y 60°

Universidad Mayor de San Andrés

 Calculando variogramas experimentales direccionales

Puntos P1 y P2 Puntos P3 y P4
caen en el la primera clase o lag Caen fuera del limite de ancho de banda

Universidad Mayor de San Andrés

Universidad Mayor de San Andrés
 Calculando variogramas experimentales direccionales

• Clases de dirección :

Para cada dirección  se define una tolerancia 

y se utilizan únicamente los puntos que se
encuentran entre las direcciones    y   
Calculando variogramas experimentales direccionales

  

puntos aceptados

puntos descartados

  
 Calculando variogramas experimentales direccionales

   


  

b puntos aceptados

puntos descartados

b = ancho de banda
Calculando variogramas experimentales direccionales y dirección



clase de distancia h

clase de distancia 2h

clase de distancia 3h
 Modelos de variogramas
• Modelo esférico o modelo Matheron
Tiene un crecimiento rápido en el origen, pero los incrementos marginales
decrecen con distancias grandes, hasta que para distancias superiores al
rango los incrementos son nulos
• Su expresión matemática es:

• Este modelo ocurre naturalmente

en todos los depósitos donde las menas
se tornan independientes entre si con la distancia (ej.: pórfidos de Cu,
depósitos sedimentarios, depósitos de Fe, Ag, Zn, Au , bauxita, Ni en lateritas,
U y fosfatos).
 Modelos de variogramas
• Modelo exponencial
Este modelo se aplica cuando la dependencia espacial tiene un creciemiento
exponencial con respecto a la distancia. El valor del rango es igual a la
distancia para la cual el semivariograma toma
un valor igual a 95% de la meseta
• Su expresión matemática es:

Universidad Mayor de San Andrés

 Modelos de variogramas
• Modelo gaussiano
Al igual que el modelo exponencial, la dependencia espacial se desvanece
solo en una distancia que tiende al infinito. Este modelo se distingue por
la forma parabólica en el orígen

• Su expresión matemática es:

Universidad Mayor de San Andrés

 Modelos de variogramas

Carencia de correlación No alcanzan la meseta

espacial No tienen estacioneriedad
en alguna dirección
También llamados De Wijsian
típico en depósitos de vetas

Universidad Mayor de San Andrés

 Anisotropías
• Cuando los variogramas en diferentes direcciones varían se dice que
existen anisotropías.

[Link]
 Anisotropía geométrica o elíptica
• Cuando variogramas en diferentes direcciones tienen la misma meseta y se
diferencian en la longitud del rango se dice que presentan una anisotropía
geométrica o elíptica.

Anisotropía elíptica o geométrica (a1 y a2 son los rangos 1 y 2, respectivamente).

Universidad Mayor de San Andrés

 Los ejes de anisotropía coinciden con los ejes de coordenadas

Elipses mostrando los ejes mayores y menores en caso de anisotropías

Geométricas. a1 y a2 son los rangos 1 y 2, respectivamente.

Universidad Mayor de San Andrés

 Modelo Anisotrópicos

Y
Los ejes de anisotropía NO coinciden con los
Y’ ejes de coordenadas
Rx
𝑢 𝑐𝑜𝑠𝛼 −𝑠𝑒𝑛𝛼 𝑥
Ry
 ℎ` = =
𝑣 𝑠𝑒𝑛𝛼 𝑐𝑜𝑠𝛼 𝑦
X 𝑢 = 𝑐𝑜𝑠𝛼 𝑥 − 𝑠𝑒𝑛𝛼 𝑦

𝑣 = 𝑠𝑒𝑛𝛼 𝑥 + 𝑐𝑜𝑠𝛼 𝑦

X’ 𝑟𝑥 𝑟 0 𝑢
𝑇ℎ` = 𝑦 =
0 1 𝑣
1) Transformar los puntos del sistema Luego de la corrección geométrica, será usado el variograma
de menor contuinidad como el variograma isotrópico
de coordenadas XY al sistema de
(Kazuo y BARBOSA, 2015).
coordenadas X’Y’
3) Evaluar el variograma isotrópico en el
h'  Rh R = matriz de rotación resultado.

2) Proceder como antes para ajustar la

longitud de los ejes de anisotropía  h  s 1  TRh 
Th' T = matriz para transformar Es un variograma anisotrópico en la dirección 
las distancias
con eje mayor igual a Rx y eje menor igual a R y
 Modelo Anisotrópicos

Y
Los ejes de anisotropía NO coinciden con los
Y’ ejes de coordenadas
Rx
Ry

X

X’

1) Transformar los puntos del sistema 3) Evaluar el variograma isotrópico en el

de coordenadas XY al sistema de resultado.
coordenadas X’Y’

h'  Rh R = matriz de rotación  h  s 1  TRh 

2) Proceder como antes para ajustar la
longitud de los ejes de anisotropía
Es un variograma anisotrópico en la dirección 
con eje mayor igual a Rx y eje menor igual a R y
Th' T = matriz para transformar
las distancias
Anisotropía zonal o estratificada
• Cuando la meseta no es la misma en diferentes direcciones,
(estratos son más homogéneos lateralmente que verticalmente).
Se divide en 2:

Universidad Mayor de San Andrés

Ejemplo de aplicación a variogramas

Universidad Mayor de San Andrés

Ejemplo de aplicación a variogramas

Universidad Mayor de San Andrés

Universidad Mayor de San Andrés
Ejemplos de variogramas (esférico)

Universidad Mayor de San Andrés

Ejemplos de variogramas (esférico)

Universidad Mayor de San Andrés

Ejemplos de variogramas (esférico)

Universidad Mayor de San Andrés

Universidad Mayor de San Andrés
Universidad Mayor de San Andrés
• [Link] = variogram(log(lead)~1, meuse)
• [Link]
• plot([Link])
• [Link] = [Link]([Link], model = vgm(1, "Sph", 900, 1))
• [Link]
• plot([Link], [Link])

• ## VARIOGRAMS SEVERAL DIRECTION (ANISOTROPY)

• [Link] = variogram(log(zinc)~1, meuse, alpha = c(0, 45, 90, 135))

• plot([Link])
• [Link] = vgm(.59, "Sph", 900, .05, anis = c(45, .4))
• plot([Link], [Link], [Link] = TRUE)

Universidad Mayor de San Andrés

Universidad Mayor de San Andrés
Ejemplos de variogramas (esférico)

Universidad Mayor de San Andrés

Ejemplos de variogramas (esférico)

Universidad Mayor de San Andrés

Ejemplos de variogramas (esférico)

Universidad Mayor de San Andrés

Ejemplos de variogramas (De Wijsian)

Universidad Mayor de San Andrés

Ejemplos de variogramas (De Wijsian)

Universidad Mayor de San Andrés

 Derechos de autor
El material presentado ha sido extraído de varias fuentes:
• Introducción a la geoestadistica: Teoría y Aplicación. Universidad Nacional de Colombia. Sede
Bogotá Facultad de Ciencias Departamento de Estadística
• Geoestatistical ore reserve estimation. Michel David (1977).
• Basic linear Geoestatistics. Armstrong M. (1998).

Universidad Mayor de San Andrés