puntes sobre estimación de recursos y reservas (página 3)

Enviado por MSc. Elmidio Est�vez Cruz

Partes: 1, 2, 3

Coeficiente de variación Interpretación

Distribución simple y simétrica de la ley. La estimación de recursos es fácil

0-25 %
y cualquier método brinda buenos resultados.

Distribución asimétrica. Dificultad moderada para la estimación de recursos.

25-100%
La distribución es típicamente lognormal

Distribución marcadamente asimétrica con un amplio rango de valores.

100 - 200 %
Dificultad para la estimación local de los recursos.

Distribución asimétrica y muy errática o presencia de poblaciones

( 200 %
complejas. La estimación local de la ley es difícil o imposible

Kurtosis: Es una medida de cuan ¨afilado¨ es el pico de la distribución de los datos.

Tradicionalmente el valor de este coeficiente se compara con cero, el cual es el valor para el
caso de una distribución normal. Un valor mayor que cero indica una distribución más
picuda que la normal y un valor menor que cero refleja una distribución más achatada que
la normal. Sin un número de muestras grandes, el uso de este coeficiente es de dudoso
valor.
Si el número de muestras es mayor que 25 entonces es posible hacer comparaciones de la
media y el coeficiente de variación de las variables de interés entre los distintos dominios
geológicos presentes en el yacimiento.
Noble, 2000 propone las siguientes reglas generales para evaluar las diferencias entre
las leyes medias:

Diferencia de la ley media entre dominios

Interpretación
geológicos (%)
 La diferencia es mínima y no es necesario diferenciar las
0-25 %
poblaciones en el modelo de recursos

Las poblaciones requieren ser diferenciadas en el modelo de

25-100% recursos si están separadas por discontinuidades(fallas etc.) o el
variograma en ambos dominios es diferente

Las poblaciones tienen que ser necesariamente separadas.

Mayor de 100 % Diferencias mayores de 100 % pueden indicar presencia de
poblaciones estériles o de alta ley

Métodos de estimación de reservas asistido

por computadoras
Los métodos de estimación de recursos asistido por computadoras se han desarrollado
ampliamente en los últimos años debido al desarrollo vertiginoso que han tenido los
ordenadores y los softwares de aplicación. Los primeros intentos estuvieron dirigidos a
informatizar los métodos clásicos o geométricos (método de los polígonos y de los perfiles)
posteriormente con el desarrollo de la geoestadística aparecieron métodos más potentes y
con una filosofía totalmente diferente de trabajo.
Los métodos asistido por computadoras permiten realizar estimaciones en bloques más
pequeños (estimación bloque a bloque, que idealmente deben ser del tamaño de la unidad
de selección minera) definidos en función del objetivo de la estimación y la densidad de
la red de exploración. En esto radica precisamente la diferencia con los métodos clásicos o
geométricos los cuales, como se conoce, definen el tamaño del bloque sobre la base de
conceptos de área o volumen de influencia que comúnmente son determinados
empíricamente, o también de acuerdo con la disposición de las intersecciones de
exploración. Estos bloques así definidos son de grandes dimensiones y se evalúan
generalizando la variable estudiada a partir de la media aritmética o la media ponderada
según sea el caso.
Los métodos computarizados para la estimación de recursos se basan
en procedimientos matemáticos de interpolación local y solamente emplean los datos de los
pozos vecinos al bloque para realizar la estimación de la variable estudiada. Un esquema
simplificado de la estimación de recursos asistida por computadora se muestra en la figura
4.1. El yacimiento o la zona mineral definido por la geología se discretiza en una matriz de
bloques bidimensional o tridimensional según el caso y cada bloque se estima utilizando los
datos localizados dentro del área o volumen de búsqueda.
Los métodos de ventanas móviles (inverso de la distancia y kriging) son las técnicas más
empleadas en la estimación de recursos asistida por computadora.
El procedimiento general de los métodos computarizados es el siguiente.
 Confección de la base de datos con toda la información relevante de la exploración
del yacimiento
 Análisis exploratorio de datos y variografía
 Creación del modelo geológico
 Modelo de recurso - División del yacimiento mineral en una matriz de bloques
regulares –modelo de bloque o capa.
 Estimación en cada bloque de las variables de interés (contenido, masa volumétrica
etc.) empleando una técnica de interpolación espacial (funciones de extensión). El valor
estimado de la ley en cada celda se calcula por la siguiente fórmula:
Z*(x) = (i=1,n Wi Z (xi) i= 1,2,3,.......n
Donde: Wi es el peso o coeficiente de ponderación asociada a cada muestra y su valor
depende del método de estimación espacial empleado, n es la cantidad de muestras
seleccionadas para hacer la estimación.

Figura 4.1 Representación bidimensional del caso general de la estimación de recursos

asistida por computadoras (Según Sinclair y Blackwell, 2002)
A continuación se describe detalladamente cada una de estas etapas. Mucho de los aspectos
aquí discutidos para los métodos asistidos por computadoras son igualmente válidos para
los tradicionales o geométricos
4.1 Confección de la base de datos
La estimación de recursos/reservas a partir de un número limitado de pozos es una de las
tareas más complejas que enfrenta un geólogo explorador. En la medida que avanza la
exploración y se van recopilando datos relevantes del yacimiento estos se van organizando
en una base de datos. El geólogo debe garantizar la integridad de la misma y tomar las
medidas necesarias para evitar los errores de trascripción o digitalización.
La precisión y exactitud de la estimación de reservas depende de la confiabilidad de los
datos sobre los cuales se basa y ningún método de estimación, por muy sofisticado que sea,
puede compensar la mala calidad de los datos (Burn, 1981). Queda claro que contar con una
buena base de datos es la garantía de una correcta estimación de reservas. Por esta razón,
desde el mismo comienzo de los trabajos es necesario implementar procedimientos que
garanticen un muestreo representativo y un estricto control sobre la calidad de los
resultados del laboratorio.
La base de datos debe integrar todo tipo de datos incluyendo litología, documentación,
resultados de los análisis, registros de carotaje, datos hidrogeológicos y geotécnicos.
Los datos de exploración (pozos, galerías, trincheras etc.) generalmente se guarden en 3 o 4
tablas
Tabla 1 recoge la posición espacial del pozo y está conformada por los siguientes campos:
Nombre del pozo (DDH;BH), coordenadas X;Y;Z de la boca del pozo (collar), profundidad y
otros datos adicionales como tipo de perforación, año de la campaña etc.
Tabla 2 recoge los datos de inclinometría que permiten determinar la posición de las
muestras en el espacio. Campos: Nombre del pozo, profundidad de la medición, azimut y
ángulo de inclinación o ángulo cenital
Tabla 3 recoge los datos referentes a las muestras tomadas (testigos) su descripción
(litología codificada) y los resultados de los análisis químicos. Campos: Nombre del pozo,
desde, hasta, # de la muestra, litología, análisis químicos. Esta tabla en ocasiones se
desdobla en dos, una para los análisis químicos y otra para la litología, pues la longitud de
los intervalos generalmente no coincide.
Una vez que la base de datos ha sido compilada el geólogo debe validarla e identificar y
eliminar todos los posibles errores. La mejor manera de comprobar la integridad de la base
de datos es mostrando la información en perfiles y comparar con el modelo que se tiene del
yacimiento. También se acostumbra a comparar un porcentaje (10%) de los registros con la
información original.
La etapa siguiente tiene como objetivo caracterizar estadísticamente cada una de las
variables estudiadas y cuantificar su variabilidad espacial. Estos aspectos son tratados en
capítulos independientes (capítulos 3 y 5).
4.2 Creación del modelo geológico
La experiencia ha demostrado que el problema principal en la estimación de recursos no
esta relacionado directamente con el método de estimación empleado sino con la correcta
aplicación de los principios geológicos. La cuestión medular a resolver antes de la
estimación de recursos propiamente dicha es establecer la continuidad de la mineralización
y la ley dentro del yacimiento. Un muestreo representativo, análisis confiables y una
coherente interpretación geológica son los componentes principales de la estimación de
recursos (Arseneau y Roscoe, 1997)
La interpretación geológica siempre se basa en los datos y en el conocimiento que posee el
geólogo sobre el yacimiento o del modelo de yacimiento que se estudia. Esto conduce
inevitablemente a que los mismos datos sean interpretados de forma diferente por distintos
especialistas, y por consiguiente resultados diferentes de la estimación de recursos.
La interpretación geológica tradicionalmente ha descansado en la construcción de planos y
secciones en los cuales se representa la morfología, dimensiones y propiedades del
yacimiento (Popoff,1966). La interpretación de todos los datos recopilados durante la
exploración se basa en 3 enfoques principales:
 Interpretación basada en perfiles y secciones
 Interpretación basada en planos de isolíneas
 Principio de analogía o inferencia geológica.
En la actualidad aunque el proceso se ha informatizado el procedimiento continúa siendo el
mismo y se basa en los enfoques anteriormente descritos.
En los métodos de estimación asistidos por computadora el conocimiento geológico que
posee el geólogo sobre la continuidad de la estructura y la ley se traduce en la definición de
zonas o dominios geológicos que poseen una forma geométrica única.
Los dominios geológicos no son más que zonas geológicamente y estadísticamente
homogéneas. Lo que realmente se hace cuando se crea el modelo geológico del yacimiento
es subdividir el mismo en subpoblaciones que cumplan o se aproximen a la hipótesis de
estacionaridad (ver capitulo 5). La división del yacimiento en dominios siempre debe
basarse en el conocimiento geológico y el sentido común. Guibal (1997) señala que la
selección de los dominios geológicos debe estar respaldada y validada por la estadística y la
variografía.
La definición de los límites o contornos de los dominios geológicos, los cuales determinan la
geometría de los cuerpos y zonas, es el método básico para aplicar control geológico a la
estimación durante la modelación de recursos.
Los límites o contornos de los dominios geológicos se clasifican en difusos o gradacionales y
físicos. En los yacimientos que están definidos por limites físicos, (carbón, yacimientos
sedimentarios, filones de oro etc.), el trazado de los contornos es relativamente fácil, pues
estos coinciden con los planos o contactos geológicos. En el caso de los yacimientos
gradacionales (Cobre porfirico) los límites se definen sobre la base de una ley económica o
cutoff.
4.2.1 Métodos para la definición de los contornos
Existen distintos métodos para correlacionar los límites de los dominios y conformar las
superficies o sólido 3D que define la geometría de la unidad geológica en el espacio
tridimensional.
4.2.1.1 Sólidos a partir de los modelos de alambre o wireframe.
El wireframe es la manera más versátil de definir los límites de los dominios geológicos.
Cualquier límite o contorno interpretado en perfiles, planos o directamente en 3D puede ser
combinado para formar el modelo de alambre y finalmente generar el sólido de la unidad
geológica.
Para generar el sólido a través de este método se digitaliza interactivamente en la pantalla
de la computadora los límites o contornos de las unidades geológicas que constituyen el
depósito (fig. 4.2). La interpretación de cada unidad en los perfiles es representada por un
polígono cerrado que posee una orientación y posición espacial, un código que indica
el dominio geológico que encierra y una potencia o zona de influencia lo que hace que cada
perímetro englobe cierto volumen.
Existen herramientas que garantizan que los contornos de unidades contiguas en un perfil
coincidan, de forma tal que no queden huecos o espacios vacíos en el modelo geológico.

Figura 4.2 Interpretación de los límites de los dominios geológicos a partir de los datos
mostrados en perfiles.
Una vez que el mismo dominio geológico ha sido interpretado en los distintos perfiles
entonces se procede a correlacionar la unidad geológica en perfiles contiguos para lograr
una representación 3D del cuerpo geológico (fig. 4.3). Los perímetros (contornos)
digitalizados se van uniendo por los puntos de inflexión correspondientes en los perfiles a
través de líneas de enlace (tie lines) hasta completar el wireframe que encierra el volumen
geológico 3D. Esta correlación sección a sección se requiere cuando los contornos son
irregulares (bifurcaciones, fallas, acuñamientos, etc.). Si las unidades geológicas muestran
una clara continuidad en 3D entonces no es necesaria la correlación interactiva y se puede
prescindir de las líneas de enlace, en su lugar se usa un algoritmo de interpolación lineal de
los límites en las secciones
Los modelos sólidos creados pueden ser cortados por planos de cualquier orientación y
espaciamiento, esto nos permite corroborar nuestras interpretaciones en otra dirección
cualquiera.
Figura 4.3 Creación del modelo sólido 3D a partir del perfiles previamente interpretados.
Como toda la información en los softwares de modelación y estimación de recursos es
manejada en el espacio no es necesario realizar la interpretación de los limites a partir de
perfiles y planos, esta puede hacerse directamente en el espacio 3D.
4.2.1.2 Sólidos generados partir de superficies (Método de superficies)
Muchas veces es posible confeccionar los sólidos modelando las superficies estructurales
que limitan por encima y por abajo los cuerpos geológicos. Este enfoque puede ser
empleado en yacimientos estratiformes de poca complejidad estructural
Este método se realiza en 3 pasos fundamentales:
 Extracción de los puntos que yacen en las superficies (techo y piso) a partir de la
intercepción de la traza del pozo con los planos estructurales)
 Modelación de las superficies (triangulación o gridding)
 Combinación de las superficies y generación del sólido
Existen 2 variantes de modelación de superficies:
Modelación de superficie por triangulación
La superficie de contacto se modela usando una red optimizada de triángulos irregulares
(TIN) que conectan los puntos con coordenadas conocidas. La triangulación se realiza en un
plano de referencia que debe ser paralelo a la superficie considerada, el mosaico de
triángulos se genera sobre la base de la coordenadas X;Y mientras que la Z define la
topología. Para generar el TIN se utiliza el método de Delaunay (Voronoi, Thiessen).
Cualquier superficie puede ser modelada (superficie topografica, planos de fallas etc.),
interceptada con el plano del perfil (viewplane) y ser utilizada para controlar la
interpretación geológica. La modelación de estas superficies controladoras es el primer
paso en la modelación geológica en 3D.
Para la triangulación se debe tener el máximo de información posible sobre la superficie a
modelar. Esta información debe estar libre de errores o inconsistencias como pueden ser
puntos con las mismas coordenadas, curvas de nivel que se corten (strings) o puntos que no
residan en la superficie.
El método de triangulación no suaviza y respeta la información original, aspecto muy
importante cuando la Z del punto se conoce con bastante certidumbre (superficie
topográfica, planos de falla), puede ser aplicado a superficies de cualquier orientación
(horizontal, inclinado o vertical), da la posibilidad de representar las discontinuidades
(breaklines) y no genera superficie más allá de los puntos extremos.
Modelación de superficie por interpolación (malla –gridding)
Representa la superficie como un malla o matriz bidimensional (formato RASTER) y
requiere la interpolación de la cota en los nodos de la red. El método tiene el inconveniente
que la representación es solo en el plano horizontal, lo que limita su uso para modelar
superficies inclinadas, además suaviza los datos originales al tener que interpolar. Para
generar la malla se utilizan diferentes técnicas de interpolación espacial (inverso de la
distancia, kriging, vecino más cercano, etc.).
La combinación de las superficies controladoras permiten el sólido que representa el
dominio geológico estudiado.
4.3 Modelo de recursos
Una vez obtenido el modelo geológico entonces resta establecer la variación espacial de las
distintas variables de interés (leyes), además de determinar los recursos, ley media y
cantidad de metal en cada dominio geológico y en el yacimiento en general.
4.3.1 Método de bloques

El método más usado en la modelación de recursos consiste en la discretización del espacio

3D en bloques o celdas tridimensionales (voxels) (fig.4.5). Cada celda contiene los atributos
(litología, tipo de mineralización etc.) y las mediciones (leyes, propiedades físico mecánicas)
del dominio geológico en que se encuentra. Los atributos de los bloques se determinan
sobre la base de la intersección con el modelo geológico o su posición respecto a una
superficie triangulada y las leyes a través de la estimación con técnicas de interpolación
espacial.
Figura 4.5 Modelo de bloque
El primer modelo de bloque fue utilizado a comienzos de los años 60 por la Kennecott
Koper Corporation en un depósito de pórfido cuprífero. Se empleó para describir la
distribución espacial de las leyes y no la geometría de los dominios geológicos.
Cada bloque debe contener toda la información disponible en las fases de desarrollo de
un proyecto: litología-mineralogía, contenidos de metales, calidades en el caso del carbón
y rocas industriales, contenidos de contaminantes, parámetros geomecánicos, datos
hidrogeológicos, etc.
Para definir el modelo de bloque es necesario establecer los siguientes parámetros (Fig.4.6)
 Posición del modelo: se especifica a partir de las coordenadas del centroide del
bloque llave (key block).
 Extensión del modelo en las distintas direcciones X, Y, Z (debe ser lo
suficientemente grande para enmarcar la región de interés)
 Dimensiones de las celdas o bloques por la X, Y y Z.
 Orientación del modelo definido (ángulo de inclinación y el azimut)