Algoritmo para la Interpolacion Espacial de Krigeado Ordinario

Abstract

El calculo de las reservas y recursos minerales es de capital importancia para las empresas geológico-mineras, tanto para la evaluación como para la planificación de la explotación de las reservas. Se puede realizar utilizando diversas técnicas, predominando los métodos geoestadísticos, los cuales tienen en cuenta las características de variabilidad y correlación espacial de los datos originales. La interpolación de Krigeado Ordinario es uno de los métodos geoestadísticos más frecuentemente usados para el calculo de las reservas y recursos minerales. Su objetivo consiste en encontrar el Mejor Estimador Lineal Insesgado a partir de los datos disponibles, los cuales generalmente son insuficientes debido al costo de su obtención. Se caracteriza por costosas operaciones de álgebra lineal que repercuten en altos tiempos de ejecución, fundamentalmente la resolución de grandes sistemas de ecuaciones lineales. Con el objetivo de disminuir los tiempos asociados a la interpolación espacial de Krigeado Ordinario, se propuso un algoritmo basado en el uso de técnicas de programación paralela, mas como métodos optimizados de búsqueda espacial; que permita resolver los problemas desestimación en tiempos razonables, fundamentalmente en el campo de las Geociencias. Este algoritmo fue implementado usando C++11 como lenguaje de programación, OpenMP 4.8.2 como biblioteca de programación paralela en memoria compartida, y Atlas CLapack como biblioteca de álgebra lineal optimizada para los cálculos matriciales. El algoritmo propuesto ofrece menores tiempos en los procesos de estimación de reservas de minerales útiles, a partir del aprovechamiento eficiente de los recursos computacionales.