Universidad Nacional de Piura

1.

Fac. Ingeniera de Minas

EL ANLISIS ESTRUCTURAL

El anlisis estructural o estudio variogrfico segn (Armstrong y Carignan, 1997)

est compuesto por:

El clculo del semivariograma experimental.

El ajuste a este de un modelo terico conocido.

El clculo del semivariograma experimental es la herramienta geoestadstica ms

importante en la determinacin de las caractersticas de variabilidad y correlacin
espacial del fenmeno estudiado (Chica, 1987), es decir, tener conocimiento de
como la variable cambia de una localizacin a otra (Lamorey y Jacobsom, 1995;
Issaks & Co.,1999), representando el til ms importante de que dispone el
geoestadstico para el anlisis del fenmeno mineralizado o de la variable de
distribucin espacial en estudio (Sahin et al.,1998; Genton, 1998a). Este anlisis
tiene como condicionantes: la distribucin estadstica, la existencia de valores
aberrantes o anmalos, la presencia de zonas homogneas o posibles zonaciones
en la distribucin de las leyes.
Puede ser calculado inicialmente el semivariograma medio, global
u
omnidireccional (ver El Semivariograma Experimental), proporcionando una idea
inicial de la variabilidad espacial de los datos, siendo el ms idneo para
representar u obtener una estructura clara y definida. Posteriormente deben ser
calculados los semivariogramas en diferentes direcciones, puede ser calculado en
4 direcciones separadas 45 con tolerancia angular de 22.5, comenzando por 0
(figura 2a) hasta encontrar la direccin de mxima o mnima variabilidad (figura
2b), pueden ser calculados tambin, ms especficamente, en 8 direcciones
separadas por 22.5. Una forma rpida y prctica de visualizar la existencia de
anisotropa es mediante el clculo del Mapa de Variogramas (Frykman y Rogon,
1993; Homand-Etienne et al.,1995; Isaaks & Co.,1999), el cual adems permitir
obtener la direccin inicial aproximada para el clculo de los semivariogramas
direccionales, permitiendo un anlisis adecuado de anisotropa. Posteriormente,
dependiendo de la continuidad espacial, es suficiente slo calcular dos
semivariogramas separados 90.

Geoestadstica
1

Universidad Nacional de Piura

Fac. Ingeniera de Minas

Ahora, el semivariograma experimental obtenido no es utilizado en el proceso de

estimacin, debe ser ajustado a ste uno a varios modelos tericos, obtenindose
un modelo o funcin analtica que caracteriza la continuidad espacial de la variable
estudiada. Los modelos de variograma terico utilizado en el proceso de
estimacin o simulacin deben satisfacer ciertas condiciones, es decir tienen que
ser definido positivo o de tipo positivo (Deutsch, 1994; Myers, 1992; Cressie y
Grondona, 1992) de lo contrario puede existir el riesgo de encontrar varianzas
negativas que no tienen sentido (Armstrong y Carignan, 1997). En general el
ajuste a modelos tericos para la determinacin de los parmetros del
semivariograma se realiza de forma visual. En ocasiones se efectan ajustes
polinomiales por el mtodo de los mnimos cuadrados u otras variantes, que
aunque se encuentra el mejor ajuste, no siempre se verifica la condicin de que el
variograma obtenido sea siempre de tipo positivo, siendo insatisfactorio (Genton,
1998b), por lo que se recomienda el uso de modelos autorizados. Finalmente debe
obtenerse uno o varios modelos de variogramas con los correspondientes valores
de meseta y alcance. El modelo de variograma seleccionado debe representar
fielmente los aspectos que se suponen importantes del variograma experimental
(Wackernagel, 1995), que sern usados posteriormente en el proceso de
estimacin o simulacin.

1.1. EL SEMIVARIOGRAMA EXPERIMENTAL

El variograma se define como la media aritmtica de todos los cuadrados de las
diferencias entre pares de valores experimentales separados una distancia h
(Journel y Huijbregts, 1978), o lo que es lo mismo, la varianza de los incrementos
de la variable regionalizada en las localizaciones separadas una distancia h.

Var{Z(x+h)-Z(x)} = 2(h)
La funcin (h) se denomina semivariograma, la cual puede ser obtenido por la
expresin.

(h)

Np ( h )
1
Z ( x i ) Z ( x i h) 2

2 Np (h) i 1

Geoestadstica
2

Universidad Nacional de Piura

Dnde:

Fac. Ingeniera de Minas

Np(h) es el nmero de pares a la distancia h.

h

es el incremento.

Z(xi)

son los valores experimentales.

xi

localizaciones donde son medidos los valores z(xi).

Esta expresin de (h) representa el til ms importante en todo estudio

geoestadstico (Armstrong y Carignan, 1997; Weerts, y Bierkens, 1993; Chica,
1987). Su clculo no consiste en una simple evaluacin de su expresin, segn se
plantea en (Krajewski y Gibbs, 1993; Journel y Huijbregts, 1978; David, 1977; Xie
y Myers, 1995a; Pannatier, 1993) esta operacin est relacionada con los
elementos siguientes:

La direccin en la que ser calculado el semivariograma, uno o dos ngulos

que definen una direccin en el espacio y/o con tolerancias angulares d
y/o d. El semivariograma calculado usando tolerancia angular de 90 se
denomina semivariograma medio, global u omnidireccional como ya se
indic.

El incremento o paso en el clculo del semivariograma h y su tolerancia lineal

dh, se recomienda que el valor de dh sea la mitad del incremento inicial.

Una distancia, que representa la distancia mxima a que pueden estar alejados
los segundos puntos del par con respecto a la lnea que define la direccin de
clculo, conocido como ancho de banda.

La distancia Lmax hasta la cual ser calculado del semivariograma. Se

recomienda que sta sea la mitad de la distancia entre las muestras ms
alejadas (Armstrong y Carignan, 1997; Krajewski y Gibbs, 1993), aunque
dependiendo de la geometra del fenmeno regionalizado en algunos casos
puede ser calculado hasta una distancia superior.

Definido los elementos anteriores, se evala la expresin del semivariograma para

todos los pares de localizaciones separadas a la distancia h que cumplan las
siguientes condiciones:

1.- La distancia entre las localizaciones xi y xi+h sea mayor que h-dh y menor que
h+dh, o lo que es lo mismo, el segundo punto del par est incluido en el espacio
definido por h-dh y h+dh encontrndose el primer punto del par en el origen o
(figura 3), este origen se mueve entre las muestras a analizar.

Geoestadstica
3

Universidad Nacional de Piura

Fac. Ingeniera de Minas

2.- El ngulo formado entre la lnea que une los dos puntos del par y la direccin
0o debe estar incluido entre -d y +d (figura 4).

3.- La distancia entre el segundo punto del par y la lnea que define la direccin de
clculo del semivariograma no debe superar el ancho de banda (Deutsch y
Journel, 1998) (figura 5).

Geoestadstica
4

Universidad Nacional de Piura

Fac. Ingeniera de Minas

Finalmente se representan grficamente los valores de (h) en funcin de h.

El grfico de (h) tiene las siguientes caractersticas segn (Armstrong y Carignan,

1997; Krajewski y Gibbs, 1993; Curran y Atkinson, 1998) (figura 6).

Pasa por el origen (para h=0, (h)=0)

Es en general una funcin creciente de h.

En la mayor parte de los casos (h) crece hasta cierto lmite llamado meseta, en
otros casos puede crecer indefinidamente. El comportamiento en el origen puede
tener diferentes formas, las cuales son segn Journel y Huijbregts (1978),
Armstrong y Carignan (1997), Chica (1987) (figura 7):
Parablico: Caracteriza a una variable muy regular, siendo continua y
diferenciable.
Lineal: Caracteriza a una variable continua, pero no diferenciable, es decir menos
regular.
Discontinuidad en el origen: Efecto de pepita, es el caso en que (h) no tiende
a cero cuando h tiene a cero. Representa a una variable muy irregular.
Discontinuo puro: Llamado tambin ruido blanco, representa el caso de mayor
discontinuidad, siendo el caso limite de ausencia de estructura, donde los valores
de dos puntos cualesquiera no tienen correlacin alguna.

Geoestadstica
5

Universidad Nacional de Piura

1.1.1.
2D

Fac. Ingeniera de Minas

Construccin del semivariograma experimental en

A continuacin se presentan ocho pasos para la construccin del semivariograma

experimental para datos distribuidos en dos dimensiones, resultado del anlisis
realizado en la bibliografa consultada.
Sea Z(x) una funcin aleatoria con N variables aleatorias regionalizadas Z(x i)
donde x = x, y es la localizacin y Z(xi) es el valor medido correspondiente.
Dados una direccin a travs de un ngulo en la cual se desea calcular el
semivariograma, d una tolerancia angular, dh una tolerancia lineal y el ancho de
banda.

Se proponen los siguientes pasos:

1.- Calcular la cantidad de pares de datos posibles por: Np = N(N-1)/2

2.- Para cada par, calcular la distancia entre las localizaciones correspondientes
por:
di

X 1 X 2 2 Y1 Y2 2

i = 1, . . . , Np

almacenando para cada i:

- P1: Nmero del primer punto del par,
- P2: Nmero del segundo punto del par,
Geoestadstica
6

Universidad Nacional de Piura

Fac. Ingeniera de Minas

- d: Valor de la distancia entre los dos puntos del par.

- Angulo que fija la direccin de la recta que pasa por los dos puntos del
par.

3.- Ordenar ascendentemente el grupo de datos anteriores por la distancia.

4.- Calcular la amplitud mxima del semivariograma L max como Lmax = Dmax/2, donde
Dmax es la distancia a que estn separadas las localizaciones ms lejanas. Esto es
la mxima distancia calculada en el paso (2), o lo que es lo mismo, el ltimo valor
despus del ordenamiento del paso anterior.
5.- Fijar una distancia h inicial conocida como paso o incremento del
semivariograma, se recomienda la distancia promedio entre las muestras
contiguas. Para los mltiplos de esta distancia ser calculada (h), por la
expresin del semivariograma. Esto indica la cantidad de puntos a procesar en el
semivariograma, el cual se puede obtener como L max / h
6.- Calcular la expresin del semivariograma para todos los pares almacenados en
el paso (2) que cumplan las condiciones siguientes:
a) La distancia d sea mayor que h-dh y menor que h+dh, es decir, h-dh d
h+dh. Si esta condicin se cumple examinar la condicin b, de lo contrario
continuar con la distancia siguiente.
b) El ngulo formado entre las lneas que parten del primer punto del par en la
direccin 0o y la que pasa por los dos puntos del par en la direccin positiva, es
decir, en contra de las manecillas del reloj, sea mayor que -d y menor que
+d, es decir, -d +d. Si esta condicin se cumple examinar la
condicin c, de lo contrario continuar con la distancia siguiente
c) La distancia entre el segundo punto del par y la lnea que pasa por el primer
punto en la direccin no supere el ancho de banda.
Observaciones:
*0
Note que como los datos almacenados en el paso (2) estn ordenados
ascendentemente por la distancia, este paso se interrumpe cuando la distancia
siguiente sea mayor que h+dh, y aqu precisamente, comienza la prxima
iteracin.
*1
Al interrumpir este paso calcular el semivariograma con los pares que
cumplieron las condiciones a, b y c, as obtenemos un valor de (h)
correspondiente al incremento h actual.
7.- Incrementar la distancia h en su propio valor, es decir, h ser el prximo
mltiplo del h inicial. Si el nuevo valor de h no supera el valor de L. Regresar al
paso (6) de lo contrario continuar el siguiente paso.
Geoestadstica
7

Universidad Nacional de Piura

Fac. Ingeniera de Minas

8.- Al finalizar el paso (7) debemos tener para cada valor transitado por h un valor
calculado de (h), los cuales sern representados en un grfico X-Y donde en la
abscisa representan los valores de h y en la ordenada los de (h). Obteniendo as
el semivariograma experimental o emprico para una direccin, incremento y
tolerancias definidas.

1.1.2.
3D

Construccin del semivariograma en tres dimensiones

Para la construccin del semivariograma 3D es necesario incorporar a la direccin

del clculo un nuevo ngulo que permita fijar unido al ngulo una direccin en
el espacio tridimensional. El ngulo debe variar entre -90o y 90o, teniendo en
cuenta que los valores extremos coinciden con la direccin vertical y son
independientes de la direccin del ngulo

La construccin del semivariograma 3D es similar al 2D con cambios en dos de

sus pasos presentados anteriormente:

1.- En el paso 2: el clculo de la distancia se sustituye por:

di

X 1 X 2 2 Y1 Y2 2 Z 1 Z 2 2

Almacenar para cada i adems: Otro ngulo que fija junto al ngulo la
direccin de la recta que pasa por los dos puntos del par en tres dimensiones.

2.- En el paso 6 punto b, la direccin que contiene a los dos puntos del par debe
estar incluida en el ngulo slido formado por la direccin del clculo del
semivariograma y la tolerancia d, con centro en el primer punto del par.
En el caso del clculo del semivariograma en tres dimensiones, an cuando
tericamente pueden ser calculados, en la prctica nos encontramos una direccin
que juega un rol diferente a la del resto (Armstrong y Carignan, 1997). En el caso
minero las variaciones a travs de los estratos es diferente a su comportamiento a
lo largo de un estrato, esto unido a la forma en que se realiza la exploracin,
varios pozos distanciados decenas de metros y cada uno contiene un conjunto de
muestras mineralizadas con una longitud del orden de 1 m. Es recomendable
entonces analizar est direccin por separado y desarrollar un anlisis de
variabilidad espacial en la direccin vertical, es decir, perpendicular a la estructura
geolgica y otro anlisis en la direccin horizontal, a lo largo de la estructura
Geoestadstica
8

Universidad Nacional de Piura

Fac. Ingeniera de Minas

geolgica, utilizando en este caso compsitos de la zona de inters, realizando

adems, un anlisis de anisotropa. Elementos que permitirn describir la
variabilidad en tres dimensiones.

1.1.3.
Problemas ms comunes encontrados en el clculo
de semivariograma
De lo expresado hasta aqu, adems de lo planteado en muchos textos de
geoestadstica, se puede obtener la impresin de que es fcil el clculo del
semivariograma experimental (Armstrong y Carignan, 1997). La fuente de
problemas que se pueden presentar en la realizacin del un anlisis estructural es
muy variada, lo que est en correspondencia con la variedad de casos que se
presentan en la naturaleza. Algunos de los problemas ms comunes discutidos en
Armstrong y Carignan (1997) son:

El valor idneo del incremento h: Una inadecuada seleccin de h puede

proporcionar un semivariograma errtico, aunque no se puede dar un criterio
exacto o aproximado sobre cual el mejor valor de h, es recomendable recalcular
(h) para distintos valores de h, hasta encontrar una forma suavizada del mismo.

Distribuciones con valores extremos: La existencia de valores extremos, altos o

bajos, en una distribucin, puede conducir a la obtencin de un variograma
fuertemente errtico. En este caso la solucin puede ser simple, eliminar los datos
extremos, porque pueden ser ocasionados por errores, en otros casos pueden
encontrarse en zonas geogrficamente distintas y pueden ser tratados de manera
separada.
Una herramienta til para la deteccin de valores extremos y encontrar el
incremento adecuado puede ser, calculado la Nube de Variogramas (Armstrong
y Carignan, 1997), el cual consiste en representar los valores de [Z(xi+h)-Z(xi)] 2/2
contra h, para cada par posible de la informacin inicial.

La existencia de poblaciones mixtas: Existen datos que pueden mostrar

diferentes poblaciones, los cuales pueden estar estadsticamente diferenciados.
En muchos casos las poblaciones estn geogrficamente diferenciadas, donde se
recomienda tratar las zonas por separado. En otros casos las poblaciones se
presenten mezcladas geogrficamente, en este caso una solucin puede ser un
cambio de escala, con lo que se logra reducir la diferencia de los valores
extremos.
Geoestadstica
9

Universidad Nacional de Piura

Fac. Ingeniera de Minas

En Krajewski y Gibbs (1993) se presentan otras razones por los que los
semivariogramas son errticos, las cuales son: 1.- No hay suficientes muestras, 2.Las muestras no son representativas del fenmeno, 3.- Las clasificaciones de las
muestras no son vlidas, 4.- El rea estudiada es no homognea, 5.- Pequeos o
largos conjuntos de datos son necesarios, 6.- Pequeas o largas distancia deben
ser calculadas, 7.- Ms o menos distancias deben ser calculadas, 8.- Pequeas
tolerancias son necesarias, 9.- Las muestras pueden tener localizaciones
incorrectas, 10.- Los valores muestreados pueden ser errneos.
El problema fundamental en la obtencin de un semivariograma correcto es, la
eleccin adecuada de los intervalos de distancias para los cuales ser calculado el
semivariograma, de modo que en stos la cantidad de pares encontrados sea
suficiente desde el punto de vista estadstico.

1.2. ANLISIS DE ANISOTROPA

Conviene aqu realizar un anlisis sobre el comportamiento de la variabilidad del
atributo en estudio. Se conoce que el semivariograma describe las caractersticas
de continuidad espacial de la variable regionalizada en una direccin, pero este
comportamiento pueden variar segn la direccin que se analice, como se discute
en Journel y Huijbregts (1978), David (1977), Zimmerman (1993), Krajewski y
Gibbs (1993). Se exige por este motivo un anlisis del comportamiento de la
continuidad en distintas direcciones, el Anlisis de Anisotropa.
Cuando el semivariograma calculado en diferentes direcciones (norte-sur, esteoeste, y en direcciones intermedias de 45 o de 22.5, con tolerancia de 22.5 o),
muestra similar comportamiento, se dice que el fenmeno es Isotrpico, cuando
muestran diferentes comportamientos es Anisotrpico (Krajewski y Gibbs, 1993).
Los tipos de anisotropas ms comunes son la Geomtrica y la Zonal. (Krajewski y
Gibbs, 1993; Journel y Huijbregts, 1978; Armstrong y Carignan, 1997)
Anisotropa Geomtrica: Est presente cuando los semivariogramas en
diferentes direcciones tiene la misma meseta pero distintos alcance (figura 15).
Anisotropa Zonal: Est presente cuando los semivariogramas en diferentes
direcciones tiene diferentes mesetas y alcances (figura 16).

Geoestadstica
10

Universidad Nacional de Piura

Fac. Ingeniera de Minas

Al respecto en (Zimmerman, 1993), se hace un estudio profundo de los tipos de

anisotropa, proponiendo una nueva terminologa. En estos casos conviene
realizar transformaciones de coordenadas con el objetivo de obtener modelos
Isotrpicos (Journel y Huijbregts, 1978; Chica, 1987; Armstrong y Carignan, 1997).

1.2.1.

Efecto proporcional

Cuando en el clculo del semivariograma se detecta que existe una relacin linear
entre el valor medio de las muestras usadas en el clculo de cada (h) y la
desviacin estndar correspondiente, se dice que existe un efecto proporcional
(heterosedasticidad). Este efecto se puede detectar ploteando los valores de X m
contra , es decir, que el coeficiente de variacin (/Xm) sea aproximadamente
constante, ocurre cuando los datos presentan una distribucin lognormal (Journel
y Huijbregts, 1978) (figura 17).

La solucin a este problema propuesta por David (1977) consiste en dividir cada
valor del semivariograma local por el cuadrado de la media local, y obtener lo que
se conoce como semivariograma relativo (David, 1977).

Geoestadstica
11

Universidad Nacional de Piura

Fac. Ingeniera de Minas

F(h) = (h)/Xm2(h)
1 N ( h ) Z (x i ) Z (x i h )
F (h )

2 N (h ) i 1 Z (x ) Z (x h) 2
i
i
2

Puede ser calculado usando los pasos anteriormente presentados para el clculo
de los semivariogramas tradicionales.
Existen otras medidas de la continuidad espacial descritas en Journel y Huijbregts
(1978) y Pannatier (1993), las cuales permiten un anlisis estructural detallado con
diferentes objetivos.

1.2.2.

Problemas en el modelaje de semivariogramas

Los problemas ms comunes al modelar semivariogramas que complican este

proceso se presentan en Krajewski y Gibbs (1993). Se analizan los siguientes
casos.

1.- La anisotropa geomtrica est presente: Indica que los semivariogramas

direccionales tienen la misma meseta pero diferentes alcances, sta puede ser
corregida a travs de una transformacin linear de coordenadas que permita
reducir una elipse a un circulo.

2.- La anisotropa zonal est presente: indica que tanto las mesetas como los
alcances son diferentes para los semivariogramas direccionales, puede ser
corregido separando el semivariograma en sus componentes isotrpicos horizontal
y anisotrpico vertical.

3.- La tendencia de los datos est presente: indica que los valores medidos
aumentan o disminuyen dramticamente en la zona estudiada con el aumento de
la distancia. Esto puede ser resuelto aplicando polinomios a la ecuacin del
semivariograma, es decir un anlisis de tendencia.

4.- El efecto proporcional est presente: Indica que la desviacin estndar local
es proporcional al cuadrado de la media local y que los datos presentan una
Geoestadstica
12

Universidad Nacional de Piura

Fac. Ingeniera de Minas

distribucin lognormal, puede ser resuelto dividiendo cada valor del

semivariograma local por el cuadrado de la media local, es decir usando
semivariogramas relativos.

5.- Existencia de estructuras anidadas: indica que diferentes procesos operan a

diferentes escalas, como por ejemplo alguno o todos los siguientes: A muy
pequeas distancias la variabilidad puede estar presente debido a cambios de una
composicin mineral a otra. A pequeas distancias la variabilidad puede estar
presente debido a errores. A grandes distancia la variabilidad puede estar presente
debido a casos transitorios de desgaste mineral. El cual puede ser resuelto
aplicando varios modelos simultneamente.

6.- Existencia de efecto hueco: indica que muy pocos pares estn disponible
para la comparacin a una distancia especfica. Y puede ser resuelto recuperando
ms casos para la distancia definida.

7.- La periodicidad est presente: indica que el comportamiento del

semivariograma repite por s mismo periodicidades, por ejemplo: El valor de la
meseta puede aumentar o disminuir sistemticamente, o un caso en que los
valores son tomados alternativamente a travs de diferentes estratos, como
piedras areniscas, esquistos, etc. Esto puede ser resuelto si es un problema real y
no un antifaz del anlisis, la periodicidad puede ser tambin un fenmeno real
mostrado por zonal ricas y pobres repetidas a espacios similares.

Geoestadstica
13