Kriging indicador con surfer 8

Cuando la variable no sigue una tendencia normal, inclusive después de realizar una transformación, se pude utilizar kriging indicador para analizar los datos ya que es una técnica geoestadística que no requiere ninguna hipótesis sobre la distribución estadística de la variable.

Los datos originales de una variable “X” deben ser transformados en un conjunto de datos codificados, en base a límites específicos definidos previamente, como resultado se obtiene una nueva variable indicadora. Para utilizar Kriging indicador se debe:

Definir los límites. Los límites pueden estar asociados a las características de los datos como límites de calidad de agua, capacidad de la roca para almacenar agua, profundidad del nivel del agua, concentración de un parámetro, niveles de ruido o para ciertos percentiles de la estadística de los datos originales.

Realizar la transformación binaria de los datos, dando a cero(0) a todos los valores inferiores o iguales al valor límite y uno (1) a los valores que están por encima. La nueva variable se construye de acuerdo a la expresión:

• 0 si X < valor límite
• 1 si X > valor límite

Con el kriging indicador obtendremos las probabilidades de que los valores estén por encima o por debajo de cierto valor límite, que se debe escoger cuidadosamente, por ejemplo podríamos obtener un mapa de la probabilidad de que el agua esté contaminado por nitrato y no sea apta para consumo humano si la concentración de este parámetro es superior a 50 mg/l.

Con esta técnica podemos definir varios límites y construir varios mapas de probabilidad para un mismo parámetro.

De igual forma nos permite trabajar con datos no mesurable, como por ejemplo podríamos tener varias unidades litológicas l1, l2, l3, l4 y l5 y, clasificarlas por la capacidad de almacenar agua, aunque no conozcamos la continuidad de cada unidad, obtendríamos un mapa que me da la probabilidad de encontrar agua en un acuífero.

Para ilustrar, el uso de esta técnica de la geoestadística, utilizaremos surfer 8 y datos de la conductividad eléctrica del acuífero. Definiremos nuestro límite de transformación 1900 uS/cm, es decir aquellas muestras que tengan la conductividad eléctrica por debajo de 1900 uS/cm les daremos una calificación de “0” y las restantes “1”. Para realizar la transformación, en Excel realizamos el siguiente condicional.

El resultado se muestra a continuación.

Una vez realizada la transformación, ingresamos a surfer 8 y construimos el variograma tal como se indica en Construcción del Variograma en Surfer 8.

El siguiente paso es construir nuestro mapa de probabilidad y para ellos debemos seguir los pasos indicados en Creación del mapa de superficie en Surfer

Finalmente obtenemos el siguiente resultado .

Recortamos el mapa, tal como se explica en Recortando contornos en Surfer. Obtenemos el siguiente resultado.

Como se observa, acabamos de obtener un mapa que nos muestra la probabilidad de que la conductividad eléctrica sea mayor a 1900 uS/cm. Es un ejemplo sencillo pero que nos ilustra la utilidad de esta metodología que podría ser aplica para analizar la litología o en minería.

Bueno finalmente quiero decirles que para aplicar kriging indicador con surfer 8, requiere el repaso y la aplicación de los siguientes temas vistos anteriormente:

• Construimos el variograma
• Creamos un Archivo GRD en Surfer
• Creación del mapa de superficie en Surfer
• Construimos el archivo bln
• Recortamos el mapa

No olvides dejar tu comentario y compartir agua y sig.