Calculo reservas de yacimiento de estaño

[pic]
UNIVERSIDAD SIMÓN BOLÍVAR

DEPARTAMENTO DE CIENCIAS DE LA TIERRA

COORDINACIÓN DE INGENIERA GEOFÍSICA

YACIMIENTOS MINERALES GC-4114

Cálculo de reservas de yacimiento de Estaño(MatLab)

Jesús A. Márquez P.

Carnet: 06-39862

Codigo en matlab:

clear all
clc

%En primer lugar se procede a cargar los datos de las coordenadas de los puntos con el comando loaddatos=load('YacimientosSn.txt');
densidad=2.5;
tenor=0.012;
Fc=0.85;
dolar=18740;

% Para asignar un vector a cada columna se hace lo siguiente:

x=datos(:,1);
y=datos(:,2);
z=datos(:,3);

% Setoman los valores máximos y mínimos de las columnas (x,y) y se hace un espaciamiento en intervalos de 10

xi=min(x):10:max(x);
yi=min(y):10:max(y);
yi=yi';

% Haciendo uso del comando griddatase realiza la interpolación usando el comando “nearest” y generamos la grafica 3d

[XI,YI,ZI]=griddata(x,y,z,xi,yi,'nearest');
surf(XI,YI,ZI)
shading interp
colorbar

% Se multiplican losvalores del area de cada espaciado por los valores de espesor ZI:

N=(100.*ZI);

% Se utiliza el comando sum para realizar la suma de las filas de la matriz N para luego sumar todas las columnas de lamatriz resultante y obtener el volumen de la reserva:

sumfilas=sum(N,1);
sumcolumnas=sum(sumfilas,2);
volumen=sumcolumnas;

% Aplicando la formula de calculo de reserva:reserva=volumen*densidad*tenor*Fc;
valordolar=reserva*dolar;
fprintf('el valor de la reserva es %e toneladas, \nla reserva en dólares es de %e$',reserva,valordolar)

% Finalmente obtenemos el valor de la reserva entoneladas y en dolares:

% Luego el grafico en 3d generado es el siguiente:

[pic]

[pic]

|Valor de la reserva en toneladas: |Valor de la reserva en dólares:|
|1.034255e+006 toneladas |1.938193e+010$ |

Haciendo la comparación con el método...