INTRODUCCIÓN
La validación de los modelos de estimación de leyes constituye una etapa imprescindible del proceso de estimación de recursos mineros que permite evaluar la calidad y confiabilidad de los resultados como base para la definición de los recursos y reservas minerales y el desarrollo del proceso de planificación y explotación mineras.
Arisaca (2012) considera que deben conocerse y aplicarse metodologías de validación que aseguren la calidad de los resultados de una estimación de recursos. Rossi y Deutsch (2014) plantean que existen razones importantes que obligan a chequear y validar los modelos de recursos, que aseguren la consistencia interna del modelo estimado y proporcionen, en lo posible, precisión con respecto a las variables estimadas.
De la misma manera, Abzalov (2016) cree que el modelo de recursos debe ser validado antes de que este sea publicado y utilizado para la evaluación de proyectos mineros y las tareas de planificación de la producción. De acuerdo con Freites (2016), los resultados de la técnica de validación cruzada no son adaptables al modelo interpolado, debiéndose aplicar una auditoría, comparando el modelo de bloques interpolado con respecto a los datos duros o compósitos.
Entre las técnicas más utilizadas para la validación de los modelos de estimación de leyes a escala global se encuentra el análisis comparativo que se realiza con las curvas tonelaje-ley a partir del modelo de bloques, con curvas patrones predictivas donde se investiga el grado de “suavizamiento” o degradación de la ley estimada, usando la corrección por cambio de soporte a partir del modelo gaussiano discreto (Pincock 2010; Rossi y Deutsch 2014; Abzalov 2016; Sama 2016), mediante el método de polinomios de Hermite (Journel y Huijbregts 1978), la denominada validación “HERCO”, o a través de las simulaciones condicionales, sirviendo en este último caso para la cuantificación del riesgo y la valoración de la incertidumbre de las inversiones mineras (Lopes, Fernandes y Vanzela 2011).
El objetivo fundamental de la investigación que aquí se reporta fue aplicar un procedimiento para la validación de los modelos de leyes y discutir los resultados en el sector Three Hills, del yacimiento aurífero Golden Hill, el cual es utilizado como caso de estudio.
MATERIALES Y MÉTODOS
Ubicación y características del Sector Three Hills del yacimiento Golden Hill
El yacimiento aurífero Golden Hill se encuentra ubicado en el municipio Jobabo, unos 8 km al NE del poblado del mismo nombre en la provincia de Las Tunas (Figura 1).
Desde el punto de vista geológico, la mineralización aurífera se encuentra asociada a la zona de oxidación (gossan) de un depósito epitermal de Au-Cu de alta sulfuración (Vila-Placeres et al. 2018) emplazado en rocas volcánicas alteradas, tobas andesíticas fundamentalmente, conformada por cuerpos oxidados que sobreyacen las menas primarias de sulfuros masivos y semimasivos enriquecidos en oro y cobre.
El gossan es de consistencia blanda a medianamente dura, manchado de hematita color rojo oscuro, localmente atravesado por venas centimétricas de óxidos de hierro que cubren superficies de fracturas y como matriz en brechas angulares. Localmente se encuentran óxidos de hierro masivo (goethita y hematita) de textura porosa, con potencias mayores de un metro, formando corazas en afloramientos y en grandes bloques.
El cuerpo mineral se destaca por el predominio de rocas cuarzo-caoliníticas fuertemente oxidadas que forman cuerpos subverticales de 10 m a 20 m de espesor. Hacia los niveles inferiores aumenta el contenido de arcilla, presentándose, entre otras, montmorillonita y nontronita. El perfil de oxidación alcanza entre 35 m y 40 m y la ley de oro aumenta con la profundidad (López-Kramer et al. 2021).
Información de partida
Como información de partida para realizar el procesamiento estadístico, geoestadístico y la construcción de gráficos y tablas se utilizaron los modelos de estimación de la ley de Au de 788 unidades de selección minera (USM) de 5m x 5m x 5m, estimadas a partir de los datos provenientes de la exploración de explotación de los bancos +80,0 m; +75,0 m; +70,0 m; +65,0 m y +60,0 m, del depósito aurífero oxidado Three Hills.
Dos modelos de estimación de la ley de Au fueron utilizados: el modelo de estimación de la etapa 1, denominado como EE1 y el modelo de estimación de la etapa 2, denominado como EE2 (Vila-Placeres et al. 2018). El modelo de estimación EE1 se hizo utilizando todos los compósitos, mientras que en EE2 se consideraron solo los compósitos dentro de la envoltura mineral con ley de borde 0,3 g/t de oro.
Para validar ambos modelos se emplearon en total 173 compósitos de 5 m de longitud ubicados dentro de la envoltura mineral. Como soporte informático para el procesamiento estadístico y geoestadístico se utilizaron los softwares GSLIB (Deutsch y Journel 1998) y SGeMS (Remy, Boucher y Wu 2009).
Procedimiento para validar los modelos de stimación EE1 Y EE2 del sector Three Hills
El procedimiento desarrollado para validar los modelos de estimación constó de 3 fases (Figura 2):
Fase 1: Simulación geoestadística de la ley de Au para 10 escenarios equiprobables y cálculo de la ley media en bloques de 5m x 5m x 5m mediante cambio de soporte.
Fase 2: Construcción de las curvas: Tonelaje (𝑇), Ley media de Au (𝑀) y cantidad de metal (Q) según leyes de corte.
Fase 3: Cálculo de la pérdida o el exceso de metal a partir de las curvas construidas para tonelaje-ley y cantidad de metal.
Fase 1. Simulación de la ley de Au para 10 escenarios equiprobables
Simulación geoestadística
La simulación geoestadística es el proceso de construcción de modelos equi-probables alternativos de la variable en cuestión, donde la reproducción de las características globales (la textura), estadística (el histograma) y variabilidad espacial (el variograma) toman precedencia sobre la exactitud local (Deutsch y Journel 1998; Arias del Toro 2008).
Simulación Secuencial Gaussiana
La simulación secuencial gaussiana (SGS) se basa en el modelo de las funciones aleatorias multigaussianas, las cuales poseen media igual a cero y varianza unitaria (σ = 1). Este método de simulación se realiza según el siguiente algoritmo (Deutsch y Journel 1998):
Transformación de los datos originales de la variable a simular al espacio normal (1).
Establecimiento de una red de puntos (nodos) que cubran el espacio del dominio a simular.
Definición de un camino aleatorio en el cual los puntos de la red previamente establecida sean visitados solamente una vez.
Estimación en cada punto u de la red, de la media 𝜇 𝑢 y desviación estándar 𝜎 𝑢 de la función de distribución condicional acumulada (FDCA), mediante Krigeaje a partir de los datos reales y los previamente simulados.
Obtención del valor simulado 𝑦 𝑢 a partir de la función de distribución condicional acumulada local (FDCA) determinada en el paso 4.
Adición del valor simulado a la serie de datos.
Los pasos 4 y 5 son repetitivos, hasta que todos los puntos (nodos) hayan sido simulados. Al final del proceso de simulación se transforman los valores gaussianos simulados al espacio real. La simulación secuencial gaussiana garantiza que los datos reales sean honorados en sus localizaciones. Cualquier número de tales realizaciones puede ser obtenido por repetición del proceso de simulación n veces con caminos aleatorios diferentes. La fase 1 se realizó según los siguientes pasos:
Análisis exploratorio de datos a los compósitos de la ley de Au (histogramas, diagramas de probabilidad, acotamiento de la ley de Au)
Anamorfosis gaussiana (normalización de los datos originales).
Variografía (cálculo y ajuste de los semi-variogramas experimentales de las gaussianas).
Simulación de la ley de Au en nodos de 2,5m x 2,5m x 2,5m (10 escenarios equiprobables) empleando el método secuencial gaussiano (Remy, Boucher y Wu 2009).
Validación (reproducción del histograma y los semi-variogramas gaussianos).
Cambio de soporte. Cálculo de la ley media del bloque de 5m x 5m x 5m m por promedio aritmético de los nodos de 2,5m x 2,5m x 2,5m ubicados dentro del bloque de 5m x 5m x 5m.
Fase 2. Construcción de las curvas
Inicialmente se calcularon las curvas de la ley (M) y el tonelaje (T) para 21 leyes de corte (desde 0 g/t hasta 3,0 g/t), de cada uno de los 10 escenarios equiprobables simulados de la ley de Au del sector Three Hills. Se hizo uso del programa “GTcurve.exe” de la librería geoestadística GSLIB (Deutsch y Journel 1998) calculada de la siguiente manera:
El Tonelaje (T) se calculó como la fracción del total de bloques del modelo, cuya ley de Au es igual o superior a la ley de corte (Zc), donde IZC es la variable indicadora, que es 1 si el valor Au ≥ Zc; si no es cero y se calcula como (2):
La ley de Au (M) se determinó como la sumatoria del producto de la variable indicadora IZC por la ley del bloque (mZC) y se expresa como (3):
Donde n es la cantidad de bloques mayores o iguales a Zc, y N la cantidad total de bloques del modelo.
Finalmente, se obtuvo la curva promedio del tonelaje y la ley de Au para los 10 escenarios equiprobables, promediando los valores obtenidos para cada ley de corte Zc de cada curva simulada tonelaje-ley. La cantidad de metal se obtuvo como el producto del tonelaje medio y la ley media.
Fase 3. Cálculo de la pérdida o el exceso de metal
En esta fase se evalúa la correspondencia entre el modelo simulado y el estimado, para lo cual se propone calcular la pérdida o el exceso del metal, además del método gráfico de comparación de las curvas. Para calcular el porcentaje de la pérdida o el exceso de metal a partir de las curvas obtenidas de la cantidad de metal, fue empleada la siguiente expresión (4):
Donde: Qest- Cantidad de Metal Estimado, Qsim- Cantidad de Metal Simulado (curva promedio para 10 escenarios equiprobables simulados)
RESULTADOS Y DISCUSIÓN
El histograma de frecuencias de los valores de la ley de Au de los compósitos (Figura 3) evidencia que la distribución estadística presenta asimetría positiva, revelando una cola de valores altos hacia la derecha del gráfico, típica de los metales preciosos como el oro y la plata, con alto coeficiente de variación.
Como valor umbral para el acotamiento de la distribución estadística del oro se determinó el valor 5,06 g/t, que corresponde con el 98 % de probabilidad acumulada (Figura 4). Las leyes de Au que se encontraban por encima del umbral fueron remplazadas por ese valor.
El coeficiente de variación de la ley de Au cortada a 5,06 g/t mostró una disminución considerable con respecto a ese mismo parámetro determinado para la ley de Au de las muestras compósito originales (Figura 5).
El análisis variográfico realizado a los valores de la ley de Au transformados al espacio gaussiano a través de los correlogramas, no permitió definir direcciones de anisotropía evidentes, obteniéndose un alcance máximo observado de 30 m y un efecto pepita asumido relativamente bajo de alrededor del 20 % (0,2) con respecto a la meseta total (1,0).
El modelo teórico general tridimensional ajustado se compone de un efecto pepita y de dos estructuras esféricas con radios máximos de 10 m y 30 m (Figura 6).
La validación de los modelos simulados de la ley de Au para los 10 escenarios equiprobables mediante el método secuencial gaussiano, según la red de nodos de 2,5m x 2,5m x 2,5m, indica que fueron reproducidos satisfactoriamente la estadística (Tabla 1), las curvas de frecuencia acumulada de la ley de Au (Figura 7), y el patrón de la variabilidad espacial del sector Three Hill, al coincidir favorablemente el modelo teórico ajustado al correlograma experimental con el correlograma gaussiano medio simulado de la ley de Au (Figura 8).
3.1. Análisis de las curvas ley- tonelaje y cantidad de metal
La comparación de las curvas tonelaje-ley de los modelos de estimación EE1 y EE2 con la curva media obtenida por simulación denota que el modelo de estimación de EE1, tanto en tonelaje como en ley, está subvalorado debido a que no fue usada una ley de borde para definir la envoltura mineral en el sector Three Hills en la etapa 1, por lo cual fueron utilizadas todas las muestras compósitos para la estimación de la ley de Au del modelo EE1, incluso las de muy bajo contenido de Au (menor de 0,1 g/t), diluyendo significativamente la ley de Au estimada.
Sin embargo, para el modelo EE2 la estimación se realizó dentro de la envoltura mineral 0,3 g/t de Au definida previamente, usándose exclusivamente los compósitos ubicados dentro de esa envoltura. De esa manera el tonelaje calculado para EE2 es más consistente, si se compara con la curva del tonelaje de la simulación media. No obstante, la ley media de Au diverge por encima de la curva media patrón obtenida por simulación indicando una sobrevaloración para el modelo EE2, debido a que se utilizó un valor umbral de 15,1 g/t de oro para crear el modelo de estimación en la etapa 2, valor de corte significativamente superior al empleado para limitar la ley de Au de los compósitos utilizados en la simulación que fue de 5,06 g/t de Au (Figura 9).
Consecuentemente, si se comparan las curvas de la cantidad de metal de los modelos de estimación en dos etapas (EE1 Y EE2) con la curva patrón de la cantidad de metal obtenida a partir del promedio de los 10 escenarios equiprobables simulados, existe para el modelo EE1 una subvaloración de la cantidad de metal, mientras que para el modelo EE2 existe una sobrevaloración de esta (Figura 10).
De acuerdo con las pérdidas o el exceso de metal calculados usando la expresión 4, el modelo EE1 exhibe pérdidas en todos los valores de cortes que van desde el 34,56 % hasta 97,85 % con respecto a la curva promedio de la cantidad de metal simulada (Tabla 2).
Contrariamente, el modelo de estimación EE2 muestra para todos los valores de corte un exceso de metal estimado a escala del yacimiento que va desde 4,20 % hasta 29,38 % (Tabla 3).
CONCLUSIONES
El tonelaje, la ley de Au y la cantidad de metal obtenidos según el modelo de estimación EE1, están subestimados si se comparan con las curvas patrones obtenidas al promediar 10 escenarios equiprobables simulados, existiendo una degradación de la ley de Au con la consecuente pérdida de la cantidad de metal que alcanza cifras superiores al 30 %. El modelo no debe ser considerado para la definición de los recursos de oro del sector Three Hills.
Las curvas de la ley de oro y la cantidad de metal obtenidos del modelo EE2, aunque se aproximan a las curvas patrones simuladas, muestran un exceso de metal que va desde 4,20 % al valor porcentual 29,38 % como consecuencia de la sobrevaloración de la ley de Au; por tanto, el modelo EE2 de la ley de Au debe ser reelaborado, limitando los valores de la ley de los compósitos a utilizar en el proceso de interpolacion. El valor de 15,1 g/t de Au resulta excesivo para propiciar un estimado aceptable de la ley de ese metal.
La comparación de las curvas tonelaje-ley y cantidad de metal de los modelos EE1 y EE2 estimados para la ley de Au con las curvas obtenidas por simulación secuencial gaussiana sobre la base de la metodología desarrollada, demuestra ser una vía eficaz para validar los modelos de recursos del sector Three Hills.