Estimación de Recursos
Una de las etapas mas importantes para evaluar un depósito o parte de un depósito de minerales

Teoría Geoestadística: Variografía - Interpretación - Modelamiento


DEFINICIÓN

Como se explicó en la formulación teórica, el variograma de una función intrínseca aleatoria esta expresado por:



En donde Z(x+h) y Z(x) son los valores de leyes en el punto x+h y x respectivamente, considerando que las variables estacionarias e intrínsicas, el promedio de Z(x+h) - Z(x) es cero, el variograma resulta ser el promedio de las diferencias.



Expresado en términos numéricos para aplicación práctica, tenemos:



Los puntos x y x+h se debe entender que pueden estar en un espacio de n dimensiones como n=1, 2 ó 3. En todos los casos es necesario determinar gama(h) en todas las direcciones posibles para identificar la orientación del comportamiento de la mineralización.

Si bien la expresión matematica indica que el variograma es función de h, en donde h es la distancia entre pares de muestras, esto significa que la función se creará en base al valor promedio de la diferencia de pares de muestras que se encuentren distanciados a h metros.


VARIOGRAMA A 1 DIMENSIÓN

Por ejemplo para el caso de cálculo del variograma en una dirección se tiene en el siguiente ejemplo, Fig. 1, leyes distanciadas cada dos metros. En la Fig. 2 se describe el cálculo de cada punto del variograma para distancias h=2m, h=4m, h=6m, etc.


Fig. 1

Fig. 2


VARIOGRAMA A DOS DIMENSIONES

Para el cálculo del variograma a dos dimensiones es necesario precisar la dirección de cálculo, esta dirección se indica en base a un rango de ángulo que deberán formar los pares de muestras. Cualquier otro par de muestras que no cumplen con el rango de orientación no formará parte del cálculo del variograma. Por ejemplo en el gráfico siguiente se presentan muestras en X e Y.


Fig. 3

En la Fig. 3 se observa que el cálculo del variograma se realizará para la dirección N80ºE considerando un ángulo de tolerancia de +/- 15º, que indica que los pares de muestras que se encuentren que intervienen en el cálculo deben una orientación entre N65ºE y N95ºE.

Como el variograma es una función de "h", se debe indicar su primer valor, en este caso se fijará en 6 m, el siguiente valor 12 m, el tercer valor 18 m, y así sucesivamente.

La forma de cálculo del variograma en esta dirección se presenta como sigue (Fig. 4)


Fig. 4

Estos nueve resultados generan el gráfico del variograma, que servirá para interpretar su comportamiento en la dirección calculada, para otras direcciones se deberán volver a realizar lo cálculos tomando pares de muestras que tengan la orientación deseada.

Para un cálculo a tres dimensiones el rango del ángulo de orientación de los pares se definirá en un cono teniendo como orientación principal el eje del cono.

Para un depósito a tres dimensiones se recomienda calcular los variogramas en todas las direcciones a fin de encontrar la dirección preferencial de la mineralización y las anisotropías


ELECCIÓN DEL MEJOR VARIOGRAMA

Una vez calculado el variograma en todas las direcciones posibles, se seleccionan los mejores que tengan características que reflejen el mejor comportamiento de correlación espacial entre las muestras. Entre las formas que adoptan los variogramas, se pueden definir las que se indican en la Fig. Nº 5:


Fig. 5

Los tipos de Variograma 1, 3 y 5 tienen una misma característica, no tienen correlación espacial entre las muestras, por lo tanto la estimación de leyes puede ser realizada con cualquiera de los métodos de estimación tradicionales (icd, triangulación, poligonación, etc.). Esto indica que no será posible medir la precisión de la estimación. Sin embargo es importante destacar que la información del variograma sobre el nivel de aleatoriedad (alto, medio o baja) permite determinar el nivel de precisión que se obtendría al aplicar los métodos tradicionales.

Los tipos de Variograma 2, 4, y 6, definitivamente indican ventaja superior en la estimación de leyes con kriging, la diferencia entre ellas se encuentra en la diferencia de precisión de la estimación. Por ejemplo para baja aleatoriedad indicado con el variogama 6, la varianza de estimación o el error relativo de estimación será menor que el que se obtenga con los variogramas 2 y 4.

Otros aspectos importantes para elegir los variogramas mas representativos del comportamiento de la mineralización son el número de pares utilizados para determinar cada punto del variograma y la altura o valores de la meseta del variograma. La meseta es la línea horizontal máxima que indica el límite en promedio del variograma.

Por lo general los primeros puntos del variograma deben tener alto número de pares, a mayor número de pares en cada punto del variograma el valor de éste será mas representativo.

Una meseta mas baja en un variograma con respecto a otro, indica que la varianza o error relativo en la estimación de leyes será menor, por lo tanto la estimación será mas precisa.


MODELOS AUTORIZADOS PARA EL MODELAMIENTO

Para el modelamiento de un variograma experimental se deben aplicar fórmulas o modelos matemáticos autorizados, estos modelos matemáticos tienen la característica de ser una función siempre positiva para cualquier valor de |h|. A continuación presentamos los modelos autorizados siguientes:

Efecto de Pepita (Co): Corresponde a un fenómeno netamente aleatorio, sin correlación entre valores, y sin importar que tan próximos se encuentren ellos. Se aplica por lo general en el origen para h=0. En la práctica se obtiene al inferir y determinar en que punto cruza el variograma experimental con el eje vertical.

La interpretación del término "efecto de pepita" se aproxima a imaginar leyes de alto y bajo valor distantes pocos milímetros o centímetros, al aplicar la fórmula del variogram las diferencias entre ellas generan un alto valor del variograma a una distancia h de casi cero metros.

Modelo Esférico: Corresponde a un comportamiento del variograma de crecimiento gradual similar a la figura Nº xx. La expresión matemática es:



para valores de h < a



para valores de h >= a

en donde "a" es el alcance, que es la distancia "h" en donde el variograma alcanza la meseta.

Modelo Exponencial: Corresponde a un comportamiento del variograma de crecimiento muy gradual similar a la figura Nº yy. La expresión matemática es:



en donde "a" es el alcance que equivale en este modelo a un tercio de la distancia que se alcanza a la meseta.

Modelo Potencial: Corresponde a la expresión matemática



en donde el exponente "n" puede adoptar valores entre cero y 2, reproduciendo lo indicado en la figura adjunta. (el signo "^n" se debe interpretar como potencia "n")

Tambien se tienen definidos los modelos Gausiano, Cúbico, función Seno. Sin embargo los tres primeros modelos descritos son los mas utilizados para el modelamiento, ya sean estos en forma independiente o combinados entre si.


MODELAMIENTO DEL VARIOGRAMA

El modelamiento del variograma se realiza una vez que se eligieron los variogramas representativos del depósito para las direcciones principales.

El modelamiento consiste en elegir el variograma experimental mas representativo de acuerdo a las características mencionadas en el punto anterior. Una vez elegido, el variograma experimental suele presentar algunas formas similares a los tipos Nº 2, 4, y 6.

En el modelamiento corresponde aplicar modelos matemáticos "autorizados" para encontrar la fórmula que representen los variogramas experimentales. Así por ejemplo en la Fig. Nº 8 se observa un variograma experimental superpuesto con el modelo autorizado que mas se le aproxima. El modelo autorizado aplicado es el "exponencial", el modelo final encontrado es el que indica para Gama(h). En donde el Efecto de Pepita Co = 0.02


Modelo de variograma de dos componentes

"Dist.Prom" es el promedio de las distancias de los pares de muestras utilizados en el cálculo del variograma



Modelo de variograma de tres componentes

En este variograma experimental fue necesario aplicar tres componentes en el modelo


Esta figura presenta un variograma ajustado con tres componentes, el primer componente "efecto de pepita" Co, el segundo componente esférico, y también el tercer componente esférico.