Modelamiento de volúmenes geológicos basado en la estimación y simulación estocástica de distancias al contacto

Abstract

El modelamiento geológico es uno de los primeros pasos en la estimación de recursos en un yacimiento. Este modelamiento se realiza a partir de una interpolación basada en un número limitado de muestras obtenidas en sondajes y tiene por objetivo reproducir de la manera más exacta lo que ocurre en el subsuelo. Esta representación es de suma importancia debido a que influye directamente sobre la toma de decisiones y, por lo tanto, en la metodología para explotar el yacimiento, la rentabilidad y otros asuntos operacionales. El objetivo de este trabajo es implementar una nueva metodología para el modelamiento de cuerpos geológicos utilizando geoestadística para estimar y simular la distancia de las muestras al contacto entre el cuerpo mineralizado y la roca caja que lo contiene. Esta nueva metodología permite incorporar información interpretativa o blanda, en los algoritmos de estimación mediante kriging ordinario y simulación secuencial gaussiana. El trabajo consistió en modificaciones a GSLIB, ésta es una librería que contiene los algoritmos más usados en geoestadística. Las rutinas modificadas fueron dos: \textit{kt3d}, que corresponde al algoritmo de kriging, y \textit{sgsim}, que corresponde a la simulación secuencial gaussiana. Estos algoritmos se modificaron para incluir la información blanda en dos modalidades. La primera consiste en sondajes artificiales, o pseudosondajes, creados por el usuario y que explicitan lo que el experto infiere con respecto a la continuidad y conectividad del cuerpo y/o de la roca caja. La segunda es usar direcciones variables, las cuales definen las direcciones de continuidad en el algoritmo. Estas direcciones se utilizan de dos maneras: orientando la búsqueda y el variograma según el campo de direcciones o bien, utilizándolas en una nueva definición de distancia no euclidiana que las incorpora. Los resultados obtenidos, para una prueba en un yacimiento real, muestran una mejora de entre 1% y 5% para el acierto total. Si bien, para dicho indicador es una leve mejora, sigue siendo una mejora, mas aún considerando que el yacimiento no posee plegamientos importantes que permitan explotar todo el potencial de estos algoritmos. De los algoritmos de estudiados, cabe destacar el comportamiento de la simulación utilizando anisotropía variable, pues mostró una mejora importante en la precisión y en el calce con el volumen real del cuerpo, aunque acompañada de una baja considerable en la cobertura. Esta característica hace de este método una buena alternativa para evaluar yacimientos de una manera conservadora, cuando la inversión se ve restringida y el riesgo debe ser el menor posible, es decir, no tanto alcance, pero con gran presición. El uso de estos nuevos algoritmos constituyen una nueva alternativa para el modelamiento, donde el usuario puede asistir con su conocimiento experto, la creación del modelo, sobre todo en yacimientos donde exista estructuras complejas que dificulten el proceso normal. La ventaja de éstos sobre otras técnicas que persiguen el mismo objetivo es lo intuitivo de las ideas detrás de las modificaciones y de la información adicional que se genera. Estos algoritmos contituyen parte del backend de una aplicación que permitirá al usuario interactuar con los datos para guiar el modelamiento de cuerpos geológicos.