Muy bien entonces en este video, vamos a entrar a la parte geoestadística porque en este video te voy a enseñar a cómo modelar un variograma y eso lo vamos a realizar con el software SGeMS, porque como te dije en principio, RecMin no puede modelar variogramas hasta el momento, se está trabajando en implementar tales algoritmos al RecMin que haga eso pero de momento como no podemos hacerlo con RecMin, lo haremos con el SGeMS, bueno si tú puedes hacerlo con el software GSLIB, en buena hora u otro software libre o si también puede ser un software comercial bueno no hay problema, como te digo la propuesta de utilizar RecMin, que hagas todo tu estudio con o sin tener que comprar un un software comercial, ok¡ Muy bien iniciamos el software SGeMS, para poder modelar el variograma pues necesito los compositos y ya los tengo aquí, si¡ Compositos de óxidos y sulfuros ¿de acuerdo? entonces ¿que es lo que vamos a hacer?, ¿para que hacemos un estudio variográfico y todo esto?, Bueno daré aquí una breve explicación, como digo no pretendo darte una clase de estadística porque supone que tú tienes que saber sus conceptos, si te recomendaría que leas toda la teoría del variograma, voy a ver si pongo información aquí en en los datos que has descargado del módulo sin en caso no están pues me lo solicitas, voy a tratar de ponerlos, también en la página RecMin esta mi monografías de tesis donde habló sobre la teoría del variograma, para que tengas claro si eres estudiante, pues ojalá que hayas cursado la materia de geoestadística pues sería excelente porque estoy seguro que sera más entendible, lo que te enseñare del software SGeMS, en realidad este vídeo es para enseñarte cómo usar el software SGeMS para calcular un variograma como te digo, no pretendo darte clase de geoestadistica, pero aclarare algunas puntuaciones para que más o menos tengas idea, ok! Entonces vamos a suponer que tengo muestras aquí, muestras acá muestras acá y quiero estimar por ejemplo, este bloque, no? Entonces tú vas a utilizar estas muestras y para utilizar estas muestras supongamos que hayan mil muestras. La teoría geoestadística dice que uses las más cercanas pero utilizando un área de influencia y utilizando direcciones de nisontropía, ¿qué significa eso? Que tú utilices un elipsoide de búsqueda, supongamos este elipsoide, y tendrá unos radios respecto del punto del centro de el que vas estimar unos radios ¿y esos radios quien te los dará? Te los dará el estudio variográfico, ¿de acuerdo? El estudio y además este elipsoide tiene una orientación porque podría ser hacia acá, podría ser hacia arriba, entonces ¿quién te dará la dirección? Te lo dará el estudio variográfico, ¿cómo es el estudio variográfico? tienes decenas o cientos de muestras, entonces empieza a calcular variograma en todas las direcciones, para encontrar la mejor área de influencia o la mejor longitud de influencia, ¿cómo se encuentra eso? tú sabes que un variograma, si recuerda la teoría leida, acá tienes la distancia H, acá tienes valores de muestra, entonces empieza a graficar variogramas para ciertas distancias la distancia mínima, por ejemplo 10 metros, 20 metros, 30 metros y así entonces vas descubriendo, vas definiendo unos puntos, que al final eso genera una curva, de repente en algunos lugares la curva saldrá, por ejemplo, constante en otros así nomas, en otros así, entonces tu tienes que buscar la curva más continua que no tenga muchos errores entonces cuando encuentras esa curva en los variogramas, en las diferentes direcciones entonces ese es la dirección principal del yacimiento en que se generó o se llevó a cabo la realización de ese evento magmático, no? Que generó ese yacimiento, ok¡ ¿Que quiere decir cuando tienes esta continuidad? Que la variabilidad por decir es mínima, es decir las muestras en esa dirección tienden a parecerse un poco, tienen familiaridad entonces por eso que sale continúa, entonces en esa dirección, si la muestras tienen una familiaridad en esa dirección, entonces en esa dirección dibujas tu elipsoide, ok¡ y los radios de la elipse te lo va a dar esta opción No se si estudiaste los variogramas, se denomina alcance entonces hay un momento en que está curva se vuelve constante y ahí ya no cambia entonces esta tendencia pronunciada que tiene aquí, te define el área o distancia de influencia entre las muestras en esa zona espacial o en esa área espacial, entonces ese alcance desde el metro 0 hasta ahí es el radio máximo de tu elipsoide, despues tienes que calcular los variogramas en la dirección perpendicular para encontrar también un alcance se supone que eso va a ser menor que éste radio que está aquí o este alcance es para tener el otro radio y también imagínate que este no es un elipse sino un elipsoide ya que está en 3D tendrías que encontrar 3 alcances, ok¡ Con el variograma que vas a realizar el kriging es el variograma, que salió mejor, no? Obviamente lo que vas a calcular son puntos y de tu ahí tienes que ajustarla a un modelo teórico si leiste la teoría pues sabes a qué me refiero, ok¡ Entonces yo te voy a enseñar a calcular el con RecMin, cómo podrías hacer luego como tu eres un geólogo modelador de recurso tienes muchos más conocimientos que yo en este tema y se que vas a sacar mejor provecho, lo que sí, te enseñare a usar las herramientas a lo mejor ya sabes manejar el SGeMS, también puedes manejar el GSLIB, ok¡ Para que veas cómo hacer esto mira, te vas aquí a la opción data análisis y aquí elijes la opción variograma, ok¡ Te sale una ventana paralela perpendicular, entonces, si vas a calcular el variograma por primera vez dejas aquí esta opción que dice computer variograma front scratch, osea desde 0, acá dice si tienes un variograma existente pues continúas en este caso como es primera vez, utilizas la primera opción. Esta ventana si quiere la maximizas para que lo tenga más claro y aquí eliges el grupo de compositos de las muestras que le vas a realizar al variograma, por ejemplo voy a elegir los compositos de óxidos, porque cómo vamos a analizar óxidos y tran, pues voy a elegir compositos óxidos y acá dice que propiedades, perdón no bloques, si no compósitos óxidos y qué propiedades o que elementos vas a utilizar. Tu puedes utilizar hasta 12 elementos a la vez por ejemplo el oro y el cobre pero ¿cuando haces eso? Cuando el estudio exploratorio de datos dice que hay buena correlación por ejemplo el plomo y el zinc tiene un coheficiente mayor a 0 punto 7 con algo así, puedes aplicar una técnica de cokriging si conoces la teoría sabrás a que me refiero cuando hablo de cokriging, pero en este caso solamente vamos a llegar hasta el hasta kriging ordinario tú también puedes aplicar kriging indicativo aquí, pero eso ya depende, pues es que tú eres especialista en geoestadística en estimación de recursos. Como te digo no es mi intención dar una clase de geoestadistica, bueno si de repente eres un geológo quizás soy un alumno para ti. Continuando aquí en esto, elegimos oro ¿por qué solamente vamos a elegir oro?, decimos que en los dos lados solamente voy a tratar el oro, ok¡ hago click en next y aquí tienes una figura como vas a calcular el variograma según la teoría por ejemplo, haremos aquí en una sola línea, si¡ Vamos a poner muestras aquí, otra muestra acá, otras aca, aca, etcétera, son todas sus puntos de muestra, entonces el variograma cuando tu diseñes tu gráfico H, vamos a suponer que esto mide 10 metros, si¡ 10 metros y así tu pones, empiezas tu variograma en 10 metros, 20 metros, 30 metros, 40 metros, etcétera, ok¡ Entonces un primer variograma para la dinstancia H, sería todos los pares de muestras a esa distancia, los pares de muestras ¿Que cosa es un variograma? Es la diferencia de estos dos puntos de muestras elevado al cuadrado, ósea esta muestra que está aquí menos la muestra que está acá, sale el resultado y eso lo elevas al cuadrado mss el otro par de muestra y así en esa dirección nada más, todos los pares de muestras que tengas en esa dirección, ok¡ Pero eso es cuando hay 10 metros, sabes que las campañas de sondaje de exploración, generalmente no son a 10 metros las muestras, puede ser 10, a 12 metros acá puede ser a 13 metros aquí puede ser a 7 metros, entonces es una distribución irregular entonces si tu buscas en 10 metros no vas a encontrar ninguna muestra, por ejemplo en este caso 10 metros llegaría aqui, no encuentra ninguna muestra para hacer un par, para poder hacer la diferencia de cuadrados y elevarlo al cuadrado, entonces lo que sugieren los geoestadísticos es trabajar con tolerancias, tolerancia es que le dices tu, bueno si no encuentras un punto en la distancia definida H, buscas a la mitad de tu distancia más acá o más allá, no? si encuentras puntos, eso lo agregas como si fuera a 10 metros ok¡ Entonces se llama las tolerancias, ok¡ Ambos lados por ejemplo aquí. tengo un par de muestras y 10 metros hasta aquí, y como no encontró, el programa busca 5 metros más allá si hubiera dicho 10 y así no¡ Entonces esa es la tolerancia que se da para poder obtener, pares de muestras en esa dirección y así en las distintas direcciones, luego tenemos la tolerancia angular, puede ser que los sondajes están en una sola dirección sino que de repente este punto de muestra este esta por aca, este por acá, no¡ Está todo desordenado entonces si no encuentra aquí y ademas no encontró en en la tolerancia, pues el programa empieza a buscar hacia ambos lados, no¡ Con una tolerancia angular, tu le dices hasta que ángulo, puedes decir hasta 20 grados, si no encuentras pues no hay, ok¡ Entonces eso es para la distancia 10, el otro variograma sería la dinstancia 20, entonces ya va a ser pares de muestras cada 20 metros, ok¡ Entonces tú tienes que tener claro la separación de muestras que tiene tu proyecto, usaría el promedio para que sea el tamaño de separación de muestras mínimo ya, por ejemplo los sondajes, si te das cuenta están separado cada 40 metros, pero a los costados, están entre 20 y 30 metros, podría escoger 30 metros, ok¡ Como te digo tendrías que leer más la teoría de cómo calcular un variograma, entonces aquí donde dice número de lags, significa cuando tú haces tu variograma, estos son números de lags, hay van 2 lags, alli que sea 30, 40, van cuatro lags, lag, es para cada distancia que calculas un punto de variograma, luego la separación es esta distancia, que en este caso, este ejemplo puse 10 metros, como te digo el tamaño del lag se elige empezando por la mínima separación de un par de muestras. Pero como tu proyecto tiene bastantes pares de muestra y no sabes cual es tu mínima distancia utilizas el promedio, ok¡ más o menos aproximado, ya¡ Y la tolerancia normalmente es la mitad de 10, o del primer lag, porque hay una vez que pase a 20, la tolerancia ¿cuál sería? Sería -5 o mas 5, sería 15 o 25 luego si el siguiente lags es 30 la tolerancia sería 105, sería 35 menos 35, ó 25 y así, número de direcciones acuérdate que si lo miras de planta vamos a poner estas muestras, mirando de planta imagina que estas son muestras de número de planta, tu puedes generar variogramas en la dirección Este variogramas, variogramas en esa dirección, también la dirección norte, o pueden hacer direcciones por ejemplo, norte angulo este o pueden ser direcciones al contrario, también, por ejemplo te mostré 4 direcciones básicas en una vista de planta también puede ser que en una vista de planta pueda ser dirección norte, dirección 0 grados, por ejemplo puede ser cada 15 grados si¡ 15 grados, acá sería 30 grados y entonces en todas las direcciones vas buscando variogramas o calculas hasta encontrar el mejor variograma o la dirección de realización de ese depósito mineral, de acuerdo? En cuanto más direcciones pues muchísimo mejor, imagínate que esto lo hice en planta, imagínate en vista 3D, ok¡ Eso realizaremos con SGeMS, no sé si te acuerdas de la teoría de variogramas o si as leído o en todo caso conoces de esto un variograma normalmente empieza esta línea continua con un efecto pepita que llamamos aqui, no? Que es el error de muestreo que normalmente se refleja en el variograma, entonces para este error de muestreo la mejor forma de encontrarlo es calculando un variograma vertical a los sondajes porque es ahí donde las muestras están más cerca y es hay donde distinguer fácilmente el error, ok¡ Entonces cómo los compositos estan cada 2 metros pues yo voy a realizar un variograma en la vertical, encuentra el valor, ese efecto pepita lo apunto luego calculo el variograma experimental y este cuando pide efecto pepita para no estarlo buscando, entre todos los datos que tengo, pues utiliza el del variograma vertical, es la forma en que yo lo veo, si está de acuerdo bueno, si tienes tu propia teoría adelante, pero lo enseñare de esta forma, es como siempre he visto a los expertos que lo hacen, ok¡ bueno yo no soy un geoestadístico experto, ¿de acuerdo? Muy bien entonces calculamos primeramente el variograma experimental, ok? Cómo sería? número de lags, por ejemplo tu has visto que un sondaje tiene aquí en este proyecto QhuYa como 200 metros cada dos metro tiene una muestra, pues digamos que tienes 100, muestras cada sondeó pero yo creo que con 20 lags o 25 lags es más que suficiente, ok¡ luego como los compositos están separados cada dos metros pues voy a poner cada 2 metros y la tolerancia pues le voy a poner también 2 metros y por si acaso normalmente sería un metro, pero como que no tiene sentido, como todos los compositores estan a 2 metros en esta parte la tolerancia no tiene sentido, tenemos una separación de muestras regular el tema la tolerancia entra cuando calculemos un variograma horizontal en mallas irregulates. ahí si las muestras, digamos que no están a una separación constante aca en los los sondajes, si¡ ok¡ Como es el variograma vertical pues calcularemos en una sola dirección, pondremos la dirección 0, porque no no hay otra el Dip le vamos a poner 90 ¿cómo se si el ángulo es negativo o positivo? Recuerda que para SGeMS en sentido horario el ángulo es positivo, para sentido antihorario es negativo tolerancia se refiere a esta tolerancia angular, no al lag de separación, no a esta toleranci de lag sino a la tolerancia angular, vamos a ponerle 22 punto 5, sí¡ ¿De acuerdo? osea ¿cómo sería eso? puede ser que el sondaje este recto hacia abajo o también hayan sondajes que están así, de acuerdo? Entonces le decimos al programa que venga calculando variogramas con tolerancia de 45 grados máximo 25 para una lado y 25 para el otro lado, bandt with es el ancho de banda que le conocemos acá, tu le asignas una tolerancia vamos a suponer que para diez metros del lag le pones un ángulo de 22 punto 5 más o menos buscaría muestras ahí pero ya cuando el lag sea de 50 metros pues el programa va a seguir buscando una tolerancia, imagínate irá abriendo, ira buscando datos por acá y eso no tiene sentido, entonces tú le pones un tope, ok¡ Para que busquen solamente allí la tolerancia es este ángulo pero en este tope del programa buscara muestras de aquí hasta acá, entonces esa tolerancia pues le vamos a elegir 10 metros, ok¡ eso queda a criterio, si dominas estos conceptos, ok¡ Muy bien nada más hacemos click en next y dice enter a positivo bandt With, ok¡ Quiere decir que el valor, está muy pequeño, para el ángulo está muy pequeños, vamos a ponerle por ejemplo 20 veremos si funciona, next y ya esta, ok¡ El programa ya cálculo, ok¡ Entonces mira aquí cómo ha salido el programa y justamente la imagen que estábamos buscando obviamente las muestras a la distancia más cercana pues se parecen más, entonces de aquí está la muestra por ejemplo estos diez metros es 10 metros, acá es dos metros o sea una muestra separada dos metros, tiene esta variabilidad pero a medida que se va separando la variabilidad va aumentando, ok¡ Aumenta hasta un cierto punto que después ya eso pierde sentido, ok¡ ósea quiere decir que hasta, hasta 20 metros, en este caso estas muestras tienen una tendencia , un comportamiento más o menos a mí similar, ok? Entonces hasta ahí sería tu alcance, ok¡ Esto sería el variograma central central, hay que ajustarlo a un modelo teórico y esto más o menos queda así, esto de acá ya no ¿de acuerdo? Lo que nos interesa aquí son dos cosas, primero que esta línea si la proyecto acá, bueno no tengo una regla, si si lo proyecto tendra un efecto pepita, un lugar donde va a cortar con el eje Y, el eje de los valores, el oro porque yo lo estoy haciendo variograma al elemento oro¡ ok¡ Esto es distancia y me interesa esta distancia del alcance, si esta distancia el alcance que más o menos aquí es el alcance donde hay ese quiebre ahí es el alcance, ok¡ Lo necesito más o menos si esto es 10,esto es 15, 16, 17 más o menos, ok¡ Te enseñare a encontrar ese alcance me va a servir como el radio de ese elipsoide, cuando llega la estimación en la dirección vertical de acuerdo entonces ese elipsoide como radio para buscar más, hacia abajo, hacia arriba en la vertical va a ser 17 y algo, ya lo encontramos y el efecto pepita lo que necesitamos aquí nada más, ok¡ Pero para hacer eso nosotros utilizamos un modelo porque este es el variograma experimental por cierto si haces click derecho sale ahí los números de pares de muestras que ha utilizado cada punto, primer punto ha utilizado 579 pares de muestras, ok¡ El siguiente 555 y todo en esa´ dirección ok¡ Perdón lo cerré vamos a hacerlo de nuevo, ok! lo cerre por error, vamos a hacerlo de nuevo, compositos compositos óxido, elemento oro, next lo haremos de nuevo, habiamos puesto 25 pares de muestras lag 2 y acá pusimos una tolerancia de 1 que era la mitad, aquí hemos puesto 0, acá pusimos 90 y la tolerancia que han opuesto 22 punto 5 y habíamos puesto acá 22 punto 5 ¿por qué? lo hicimos por que queremos que busquen un ángulo 45 esto por esto a 45 grados y tolerancia habíamos puesto 20 ok¡ Click en siguiente, ingresamos bandh with positivo, escribimos 25, creo que está bien, así no vamos a ver creo que acá puse 25 ó 29 no recuerdo, pues que te digo como esto cerró, hasta que tu cambio los datos, y a la cosa cambia, así creo que fue, si esto es o fue¡ ok¡ Muy bien entonces vamos a crear un modelo aquí, si¡ Esto vamos a ver qué modelo puede tener vamos a poner exponencial si no movemos las barritas se mueve la línea del modelo ok¡ Entonces este variograma, cuando tu mueves línea negra es ajustar el modelo del variograma que vamos a utilizar, pero no sabemos si es exponencial esférico, gaussiano. para eso tienes que leer de los modelos de variogramas seegun la teoria geoestadística voy a tratar de poner información aquí en los módulos de datos sino pues lee el manual de lópez jimeno, es un libro bueno en teoría para los principiantes, para entender la geoestadística ok¡ Entonces aquí donde dice meseta vamos a poner por ejemplo, vamos a proyectar esto como meseta sería 2 punto 4 más o menos, ok¡ Entonces vamos a escribir aquí, perdón acá donde dice Sil que significa meseta en español 2 punto 4 y hay más o menos está ok! entonces movemos esto de aquí y te das cuenta que ese modelo exponencial pues no ajusta con el comportamiento que tienen esos datos ahí entonces intentamos con el esférico y eso ajusta mejor mira te das cuenta, encaja mucho mejor, ok¡ Entonces el efecto pepita más o menos vendría a ser aproximadamente de 0 punto 2 Nugget effect vamos a ponerle 0 punto 2 y al aumentar el efecto pepita se suma a la meseta, entonces a esta le bajamos 0 punto 2 ok¡ espero dejarme entender. Repito desde 0 hasta aquí hay 2 punto 4 pero como ya tengo el efecto pepita 0 punto 2 solo agregó 2 punto 2 para completar la meseta, ok¡ Vamos a rechazar de nuevo vamos allí a modelar y ahí estaría el modelo eso se llama ajustar a un modelo teórico ¿de acuerdo? Muy bien entonces ya tenemos el efecto pepita que 0 punto 2 incluso puede ser menos si eso ya queda a criterio del modelador, no se concluye un valor exacto, puedo decir un valor aproximado por ahí, pero cada uno de acuerdo a su criterio decir no para mí el efecto pepita es 0 punto 18, ok¡ Eso queda en manos de los geoestadisticos te invitó que si eres nuevo en este tema pues te informes en la teoría de geoestadística, de cuerdo? Muy bien entonces es un modelo esférico y de esa manera se modela aunque aquí el modelo no interesa, porque no es el utilizaremos en la estimación. solamente sirve para identificar el efecto pepita así que voy a abrir un bloc de notas donde voy apuntando todo, bloc de notas ahí está y voy a escribir aquí efecto pepita igual 0 punto 2 gramos de oro por tonelada ok¡ Es el error de muestreo digamos no y es el programa vertical ok y entonces el radio el alcance si tú estás viendo sale 19, obviamente necesitas tres radios radio mínimo, radio medio, radio máximo, como te digo como es cuestión de criterio de repente para ti puede decir que el radio es 18, otro geólogo puede decir que es 20 otro puede decir 17. Entonces para evitar esos líos el programa usa una especie de intervalo, mínimo medio y máximo algo cualitativo, pero es un valor que se usará en el cálculo seguro que promediará los tres entre el máximo y mínimo. Aquí por casualidad los tres juntos no permanente siempre queda una más adelante otro así no es el programa lo que supongo es que promedia los tres y eso es el alcance utilizando y se obtiene un radio así. No exacto sino así en tolerancia. Hay otro software que brindan el alcance exacto pero éste no más bien te brinda una especie de tolerancia muy bien esto los ángulos no lo manejo muchos y si tú eres geoestadístico y lo conoces entonces lo puedes manejar ok! muy bien entonces el alcance es digamos que es 19 ok y alcance vertical escribimos igual 19 metros entonces ese radio para la elipsoide o tamaño de la vecindad de búsqueda que va a ayudar a la estimación ¿correcto? lo dejas ahí ya está ok si quieres esto lo puedes guardar dice si quiere borrar el modelo guardas el modelo si quiere guardar el variograma experimental también esta opción es abrir modelo esta otra es guardar y también lo puedes imprimir como imagen si tu quieres tener una máquina es fácil, solo eliges plotear imagen hicimos en el curso anterior como plotear imágenes de los histogramas es sencillo hacerlo ok! muy bien entonces cerramos esto ya tenemos los datos que queríamos y ahora sí vamos a calcular el programa experimental eso hagamos en otro video para que no se haga tan extenso este video de acuerdo nos vemos entonces otro vídeo,