En la industria minera, las simulaciones geoestadísticas se utilizan para cuantificar la incertidumbre en las leyes de mineral y predecir los recursos recuperables sobre una determinada ley de corte. Para estos efectos, uno de los modelos más utilizados es el modelo multigaussiano, el cual permite simular las leyes de un material de interés a soporte puntual, para luego rebloquear los valores puntuales y construir un modelo de leyes a tamaño de bloques. En este trabajo, se propone estudiar el modelo gaussiano discreto, que permite realizar directamente las simulaciones a tamaño de bloques sin pasar por simulaciones puntuales, con los beneficios en tiempos de cálculo que esto significa, pero también con algunas aproximaciones en las que incurre con respecto al modelo multigaussiano. La primera parte del trabajo apunta a establecer las condiciones bajo las cuales el modelo gaussiano discreto entrega resultados similares al modelo multigaussiano. Así, se realiza un estudio de sensibilidad al variograma, que modela la variabilidad espacial de la ley, así como al tamaño de los bloques a utilizar en el modelo, a la cantidad de datos condicionantes y a la asimetría que presente el histograma de estos. En la segunda parte del trabajo, se aplica el modelo multigaussiano y el gaussiano discreto a dos bases de datos de leyes reales con distintas características, con el objetivo de evaluar en la práctica los resultados de estos modelos. En particular, se compara las estimaciones de recursos mediante curvas tonelaje-ley y las medidas de incertidumbre asociadas a cada estimación. Los resultados obtenidos indican que, con una base de datos con un histograma de asimetría leve o moderada, el modelo gaussiano discreto aproxima de buena forma al modelo multigaussiano, observándose distribuciones de leyes simuladas y curvas tonelaje-ley prácticamente idénticas, además de medidas de incertidumbre similares. Sin embargo, al aplicar el modelo gaussiano discreto en una base de datos con un histograma de asimetría considerable, la similitud de resultados con el modelo multigaussiano se pierde: se observa diferencias importantes en la estimación de recursos recuperables, especialmente en las leyes medias obtenidas para distintas leyes de corte, y en las medidas de incertidumbre obtenidas con cada modelo. Las conclusiones de este trabajo indican que aplicar el modelo gaussiano discreto a datos con un histograma muy asimétrico no es recomendable. En estos casos, el modelo gaussiano discreto tiende a sobrestimar la cola de altas leyes en el histograma de frecuencias de leyes simuladas, y además presenta diferencias importantes en las frecuencias de leyes menores, todo esto en comparación al modelo multigaussiano. Sin embargo, con una base de datos de histograma con asimetría leve o moderada, el modelo gaussiano discreto es perfectamente aplicable, además de requerir tiempos de cálculo considerablemente menores.

La evaluación de recursos minerales es esencial para el diseño y la planificación minera, dado que cuantifica la distribución de elementos de interés, subproductos y contaminantes dentro de un depósito minero. Tradicionalmente, los modelos de recursos se elaboran mediante ponderación del inverso de la distancia o kriging, considerando una variable a la vez e ignorando las correlaciones espaciales entre las especies minerales. Métodos de estimación multivariable, como el cokriging, siguen siendo poco utilizados en la práctica, debido principalmente a la dificultad de ajuste de un modelo variográfico. Este trabajo aborda el problema de la predicción multivariable de recursos recuperables en depósitos mineros. Con este objetivo, un modelo geoestadístico (el modelo Gaussiano discreto) es utilizado para cosimular las leyes de bloques y así determinar sus distribuciones locales conjuntas, considerando además el caso en que las leyes medias son desconocidas. De las distribuciones locales obtenidas, es posible calcular, para cada bloque (unidad de selectividad minera), el valor esperado de cualquier función de las leyes y evaluar los recursos que pueden ser recuperados por sobre leyes de corte dadas. La metodología es aplicada a un caso de estudio que consiste en datos de producción (pozos de tronadura) de un depósito de lateritas niquelíferas. Las variables de interés corresponden a leyes de níquel, fierro, cromo, alúmina y sílice. Se demuestra que el modelo Gaussiano discreto permite estimar el beneficio esperado de unidad de selectividad minera y su mejor destino (planta de procesamiento o botadero). Estas estimaciones dependen de la probabilidad de superar leyes de corte dadas y de la razón entre las leyes de sílice y magnesio, la cual juega un papel importante en el procesamiento metalúrgico para obtener ferroníquel. Los resultados de esta metodología son comparados con una estimación tradicional mediante cokriging ordinario, arrojando un 24.6% de bloques clasificados de forma diferente (entre estéril y mineral) y una discrepancia de un 22.5% en la relación estéril/mineral estimada y de 8.2 [MUS$] en el beneficio estimado para el caso de estudio. Estas discrepancias se explican por la propiedad de suavizamiento del cokriging, que no reproduce la variabilidad espacial de las leyes y produce sesgos en la estimación de variables no aditivas. En cambio, la metodología propuesta entrega resultados teóricamente insesgados y permite el análisis de escenarios y el cálculo de la respuesta más probable, en particular cuando hay involucradas variables no aditivas.