1 - Definición de Las Unidades de Estimación
1 - Definición de Las Unidades de Estimación
1 - Definición de Las Unidades de Estimación
1. Propsitos
Determinar regiones bien definidas dentro del depsito denominas unidades de estimacin que
constituirn la base para la estimacin de los recursos mineros mediante la tcnica de kriging.
2. Responsabilidades
Es responsabilidad del Ingeniero Jefe de la Unidad de Evaluacin y Diseo de Mina, velar por la
difusin y cumplimiento de este procedimiento a todo el personal relacionado con la ejecucin
del control de calidad y de la validacin de los resultados obtenidos, tanto trabajadores propios
como contratistas de empresas de servicios, el presente instructivo y/o procedimiento.
3. Definiciones
3.1 Unidad de Estimacin
Se denomina unidad de estimacin (UE) a aquella porcin del espacio contenida dentro de la
regin en estudio y se caracteriza por la estacionariedad de la variable a estimar, es decir, la
esperanza y la varianza son independientes de la posicin espacial dentro de la unidad de
estimacin.
Dependiendo de la escala de trabajo, las distribuciones de la mineralizacin reflejan el
comportamiento general de una funcin aleatoria (por ejemplo la ley de hierro). Estas
distribuciones varan drsticamente de un lugar a otro, especficamente a pequea escala. A
pesar de esto, estas distribuciones exhiben cierto grado de permanencia al interior de fronteras
bien definidas en el depsito que se denominan unidades de estimacin.
3.2 Variograma Cruzado
El variograma cruzado consiste en el promedio de diferencia de leyes entre dos unidades
geolgicas distintas a la distancia . Esta herramienta permite cuantificar cuanto se
correlacionan espacialmente muestras pertenecientes a ambas unidades.
Al igual que para el caso de una sola variable, se puede calcular el variograma cruzado
experimental por medio de un estimador de dicha funcin, que se define como:
()
1
=
[ ( + ) ( )][ ( + ) ( )]
2
=1
5. Descripcin de la actividad
5.1 Anlisis de las poblaciones de los datos
El primer criterio para definir las unidades de estimacin se basa en la distribucin estadstica
de los valores de Fe de los compsitos, donde se pueden observar distintas poblaciones de datos
en la grfica de probabilidad.
5.1.1
Una vez copiada la columna con leyes de hierro de los compsitos al software Minitab se debe
acceder a la opcin Grfica de probabilidad en la pestaa Grfica de la barra de herramientas,
ver figura 1.
Luego se debe seleccionar la opcin Grfica Individual y aceptar. Aparecer un cuadro de
dialogo como se muestra en la figura 2 donde debe seleccionar la variable FE (u otro nombre
dependiendo de cmo haya sido nombrada).
Grfica de probabilidad de FE
Grfica de probabilidad de FE
Normal - 95% de IC
Normal
99.99
Media
Desv .Est.
N
AD
Valor P
99
80
99
95
80
Porcentaje
Porcentaje
95
99.99
22.49
11.71
2605
23.567
<0.005
50
20
50
20
0.01
-20
-10
10
20
FE
30
40
50
60
0.01
70
10
20
30
FE
40
50
60
Figura 3. Grfica de probabilidad. A la izquierda grfica por defecto y a la derecha la grfica modificada para
visualizar mejor las poblaciones.
5.1.2
Grfica de probabilidad de FE
Normal
99.99
99
Porcentaje
95
80
50
20
5
1
0.01
10
20
30
FE
40
50
60
De acuerdo a lo anterior, debemos ser flexibles y acomodar las poblaciones con las unidades
geolgicas definidas, por lo que, se puede considerar que las UG 25 y UG 45 (UE 2545) tienen
las misma pendiente y se podran considerar como una unidad de estimacin. Pero se observa
que existe un cambio de pendiente significativo en el grfico en los valores con leyes de hierro
prximos a 15% por lo que se podran considerar las UG 10 y UG 15 como unidades de estimacin
separadas (UE 10 y UE 15).
Unidades de Estimacin Unidades Geolgicas
UE 10
UG 10
UE 15
UG 15
UE 2545
UG 25 + UG 45
Tabla 2. Definicin de las unidades de estimacin segn las
poblaciones observadas en el grfico de probabilidad.
Teniendo ya creado el archivo de los compsitos (figura 5) debemos dividir este archivo en
funcin de las unidades geolgicas (UG). Las UG fueron definidas espacialmente en base a las
leyes de hierro del cuerpo, pero esto no quiere decir que por ejemplo que la UG 45 contiene
exclusivamente compsitos con leyes mayores a 45%, sino ms bien que esta UG contiene
predominantemente compsitos con leyes superiores a 45% pero tambin contiene leyes
menores. Es por tanto que la separacin del archivo compsitos no es funcin de la variable Fe
sino de las coordenadas x,y,z del compsito.
Necesitamos entonces cargar en el software los archivos DXF que contienen el modelo
tridimensional de las unidades geolgicas. Para ello debemos hacer click en la opcin Import
External Data Into Project como se muestra en figura 6 y luego seleccionar Advanced DXF/DWG
(AutoCAD). Una vez hecho esto aparecer un cuadro de dialogo donde se deber buscar el
archivo DXF deseado. Repetir el proceso hasta cargar todas las unidades geolgicas en el
software. Se crearan dos archivos por cada UG, uno contiene los puntos o vertices del slido y
el segundo contiene las triangulaciones.
Para dividir los compsitos en funcin de las unidades geolgicas debemos ir a la opcin
Wireframes en la barra de herramientas y luego seleccionar la opcin Wireframing Processes y
finalmente a la opcin Select Data Using Wireframe.
En la pestaa Fields se deben asignar el nombre que poseen las coordenadas este (X), norte (Y)
y la cota (Z) de los compsitos dentro del archivo de entrada en la rutina.
En la pestaa Parameters en el campo SELECT se debe ingresar el nmero 3, esta opcin
seleccionar los compsitos que se encuentran dentro del slido. Finalmente presione aceptar,
deber repetir el proceso y generar un archivo de compsitos para cada UG.
Para que al realizar el variograma cruzado el software reconozca que cada compsito pertenece
a una u otra UG debemos ingresar un cdigo a la tabla de los compsitos por cada unidad
geolgica para su identificacin. Para ello debemos abrir la tabla de los compsitos recin creada
haciendo doble click sobre su nombre en el panel desplegado a la izquierda de la pantalla, los
archivos se encuentran en la categora Assays. Se desplegar la tabla mostrada en la figura 10.
Para crear una nueva columna debemos ir al final de la tabla y seleccionar la ltima columna
llamada RADIUS (N), luego seleccionar la opcin Insert y despus Column en la barra de
herramientas. Debemos asignar el nombre CODE a la nueva columna y en Default Value
debemos asignarle un nombre para diferenciar esta unidad geolgica de las otras, en nuestro
caso 45. Presione OK.
Se agregar la columna CODE (N) antes de la columna RADIUS (N), esto es muy importante
porque la columna CODE (N) debe estar en la misma posicin en todos los archivos de
compsitos. Haga click en guardar y le aparecer un cuadro de dialogo donde deber presionar
la opcin Studio. Repetir el proceso para cada archivo de compsitos y asignar un cdigo distinto
para identificarlos.
Ahora debemos unir los archivos de compsitos pertenecientes a dos unidades geolgicas en
contacto, como ejemplo, uniremos las unidades geolgicas UG25 y la UG45. Slo es posible unir
dos unidades geolgicas a la vez, ya que el variograma cruzado trabaja con pares de datos. Al
haber asignado previamente un cdigo a cada UG el nuevo archivo identificar claramente a que
unidad geolgica pertenece cada compsito. Para realizar esta tarea debemos ir a la opcin
Applications, File Manipulation Processes, Merge y finalmente Append.
Aparecer la interfaz mostrada en la figura 13. En la pestaa Files deber ingresar los archivos
de compsitos a unir en los campos IN1 e IN2. En el campo OUT deber ingresar un nombre
para el archivo de salida y presionar Aceptar.
5.2.2
Para crear un variograma cruzado, herramienta que nos indicar la correlacin espacial entre
los compsitos perteneciente a cada una de las UG debemos ir a la opcin Applications,
Statistical Processes, Variograms y finalmente Calculate.
Dentro de la rutina variograma, en la pestaa Files se debe ingresar el archivo que contiene los
datos de los compsitos en el campo IN y se debe asignar un nombre al archivo de salida en el
campo OUT.
En la pestaa Fields debe ingresar el nombre de las coordenadas X,Y,Z del compsito dentro del
archivo, por defecto es el mismo nombre. En el campo F1 se ingresa la variable a analizar su
correlacin espacial mediante un variograma, en nuestro caso hierro (FE).
En la pestaa Parameters (Figura 18) se deben ingresar los parmetros para la construccin de
o los variogramas experimentales. En el campo LAG se debe ingresar la distancia () del paso del
variograma, es decir, cada cuantos metros va a buscar datos, generalmente se usa 10 m. El
campo LAGTOL se ingresa la tolerancia del lag () es decir entre que intervalos de distancia
del lag buscar datos. Por defecto el programa utiliza la mitad del lag. En el campo NLAGS se
ingresa la cantidad de pasos que se realizaran, se recomienda utilizar 25 pasos, lo que significa
que en total se buscaran puntos hasta una distancia mxima de 250 metros (considerando que
el lag mide 10 m). En el campo AZI se ingresa la primera direccin de forma azimutal donde se
calcular el variograma, se recomienda empezar desde 0. En el campo HORANG se debe
ingresar la tolerancia angular en grados alrededor de la direccin para buscar los datos, se
recomienda usar 15. En el campo HORINC se debe ingresar el valor del siguiente azimut para
realizar el variograma, se recomienda usar 30, que corresponde a dos veces la tolerancia
angular. En el campo NUMHOR se ingresa la cantidad de direcciones o azimut donde se calcular
el variograma, para cubrir los 360 se requieren 14 variogramas. En el campo DIP se ingresa el
ngulo de manteo del primer variograma, se recomienda el valor de 0. En el campo VERANG se
ingresa la tolerancia angular en la vertical, se recomienda 15. En el campo VERINC se ingresa la
siguiente direccin vertical donde se calcular el variograma, se recomienda 30. En el campo
NUMBER se ingresa la cantidad de variogramas a realizar en la direccin vertical, para cubrir
toda la vertical se requieren 4 variogramas. Para especificar en la rutina que estamos realizando
un variograma cruzado debemos en el campo KEYMETH ingresar la opcin 2 que corresponde a
la construccin de variogramas solo usando pares de datos con distinto cdigo.
Ahora debemos volver a la pestaa Fields e ingresar en el campo KEY la variable que distingue
cada compsito a la que llamamos CODE. Finalmente presionar Aceptar.
El software calcular cada variograma y luego de eso nos preguntar si continuar, debe hacer
click sobre el ticket verde, deber hacerlo aproximadamente 55 veces hasta que aparezca el
mensaje VGRAM COMPLETE.
Una vez terminado el proceso se crear un archivo con el nombre asignado. Las filas asignadas
con un guin ( - ) en las columnas azimut y dip corresponden al variograma omnidireccional. Las
columnas AVE.DIST (N) y VGRAM (N) corresponden a la distancia () y la varianza (())
respectivamente.
Cuando comience a ejecutarse la rutina le preguntar si desea continuar, presione el ticket verde
una vez, tal como lo hizo para calcular el variograma. Luego se generar una lista con todos los
variogramas numerados, el software le pedir que ingrese el nmero del variograma que desea
visualizar. Para nuestro anlisis utilizaremos el variograma omnidireccional que es el primero de
la lista, por lo tanto ingresaremos el valor 1 en el cuadro de comandos y presionamos Enter.
La versin 3.0.1742.6 del software presenta una interfaz de visualizacin limitada. Versiones
superiores han mejorado esta opcin.
En el variograma observamos que existe una correlacin entre los datos hasta una distancia
mayor a 100 metros. Esta versin de Datamine es limitada y no permite ajustar fcilmente un
modelo de variograma para determinar exactamente el alcance. Dado que el alcance del
variograma es alto, es decir, existe una correlacin espacial entre los datos de las dos unidades
geolgicas, con esto validamos lo planteado con el mtodo de las poblaciones estadsticas de
que las UG 25 y UG 45 deberan considerarse como una sola unidad de estimacin (UE 2545).
Ahora debemos analizar la influencia entre el resto de las unidades geolgica, para ello debes
repetir el proceso anterior.
5.2.3
Las imgenes siguientes muestran los variogramas obtenidos para las distintas unidades
geolgicas en contacto.
Para caracterizar la frontera entre ambas unidades geolgicas se deben copiar las coordenadas
de los puntos o vrtices del slido de cualquiera de las dos unidades de estimacin, figura 30,
este archivo fue generado en Datamine cuando se importaron los cuerpos geolgicos desde un
archivo DXF.
Si existen datos anteriores presione el botn Limpiar para borrarlos automticamente. Luego
debe ingresar la regin del espacio desea estudiar limitada por sus coordenadas mximas y
mnimas, ver figura 31. El tiempo de clculo depende de la cantidad de datos, por lo tanto se
recomienda primero seleccionar un rea pequea no superior a los 200200200 metros
cbicos y luego ir aumentado progresivamente el tamao de la regin de acuerdo a los
resultados obtenidos.
Finalmente presione el botn Buscar en la Zona y espere hasta que el proceso termine, puede
tardar varios minutos dependiendo de la cantidad de datos. Los resultados se desplegarn en
una tabla a la derecha de la hoja de clculo. Para visualizar la grfica obtenida vaya a la pestaa
Grfico Anlisis de Frontera y ver algo similar a lo siguiente.
Anlisis de Frontera
70
60
Ley Fe (%)
50
40
30
20
10
0
-60
-40
-20
20
40
60
En la grfica se observa a la izquierda del eje (lnea azul punteada) la ley de hierro de los
compsitos que pertenecen a la UG 45 mientras que a la derecha los pertenecientes a la UG 25.
Se aprecia claramente que no existe un quiebre marcado en cuanto a las leyes de los compsitos
a cada lado de la frontera, es cierto que en la UG 45 hay menor cantidad de compsitos y estos
son de mayor ley, pero se debe principalmente a dos razones, la primera es que la UG 45 es ms
pequea y posee un menor nmero compsitos que la UG 25 y tambin puede influenciar en la
cantidad de datos la eleccin de la regin a analizar donde puede que se concentren
mayoritariamente compsitos de una de las unidades geolgicas. No obstante podemos concluir
que existe una transicin suave entre ambas unidades geolgicas y por lo tanto podemos
considerarla como solo una unidad de estimacin.
6. Referencias
M. Alfaro, 2013. Tendencias en la estimacin de recursos mineros.
J. Mateu, I. Morell, 2003. Geoestadstica y modelos matemticos en hidrogeologa.
E. Tulcanaza, 1999. Evaluacin de recursos y negocios mineros.