IMM 3310.

PLANIFICACION MINERA

CLASE 3. Teoría de los recursos

Naturales. Pit Final. Optimización.

Profesor: Felipe Bernal R.

Centro de Minería PUC

Marzo, 2012
Modelo de Beneficio Económico
 Suponga una función de beneficio del tipo:
 t  pt q t  C ( q t )
qt Ritmo de producción de metal en periodo t

pt Precio del producto en periodo t

C (qt ) Costo de producción como función de la cantidad producida

Modelo Económico Simple
 Se simplifican las restricciones para realizar la
extracción total
 La utilidad marginal descontada debe ser
constante a través de la vida del proyecto
T  T
 T

Max t t  L   t  t     R0   qt 
 t 1  t 1  t 1 
sa T
L
q t  R0  t t'    0
t 1 qt
L T
 0   qt  R0
1'   2 2'  ...   n n'    t 1
Teoría de Hotelling
 Si existe mercado para solamente comercializar qt
 Mercado de competencia perfecta
 No existe asimetría ni de información ni tecnológica
 La extracción mineral para que sea óptima debe
generar un perfil de precios que aumente en valor
absoluto en la misma tasa que la tasa de descuento
utilizada en la industria En otras palabras cualquier programa
de producción fija el valor de la mina
 t'  pt  p1   2 p2  ...   n pn en forma tal que el ingreso en cualquier
tiempo, después de tener presente el
agotamiento, es exactamente igual al
interés sobre el valor de la inversión en
ese tiempo.
Modelo a Optimizar
 Si se desea optimizar el valor presente neto
del proyecto o plan de producción
T 
Max t  pt qt  C (qt )
 t 1 
sa
Rt 1  Rt  qt t  1..T

t Factor de actualización de los flujos obtenidos en el periodo t

Rt Recursos Remanentes en el periodo t
Solución
 Supuestos
R1 Conocida la envolvente de reserva base inicial
T Es un parámetro en la optimización y no una variable

 Langrangiano
T T
L   t  pt qt  C (qt )   t  qt  Rt  Rt 1 
t 1 t 1
Solución del Langrangiano
T T
L   t  pt qt  C (qt )   t  qt  Rt  Rt 1 
t 1 t 1

L

 t pt  Cq   t  0
qt
pt  Cq   t
L
 0  Rt 1  Rt  qt
t
El valor óptimo de producción se
produce cuando la utilidad marginal
Es igual al costo marginal más un
costo de depreciación del activo
mineral
Definición de Mineral
 Mineral es todo aquel porción de un yacimiento minero que
paga sus costos de producción y el costo de oportunidad
 Definición económica
 Si embargo se debe diseñar con una envolvente económica
que pudiese contener material estéril en su interior
Roca de Caja

Cuerpo Mineralizado
o Mena
Definición
 Criterio optimizante que pretende incorporar
el costo de oportunidad tanto en la definición
de la envolvente económica como en la
definición de mineral en el tiempo
Interpretación de Lane
dV *   * dV *  
 Maxw c( w, t )   V  
dR   dT 

 * dV * 
c( w, t )   V  0
 dT 

 * dV * 
 V    t
 dT 
Enlazando la Teoría de Lane
 El óptimo de producción se produce en el
punto donde la estrategia de producción paga
los costos marginales y de oportunidad

 * dV  *
 V    t
 dT 
Parámetros formulación algoritmo de
Lane
f costos fijos $/año
Etapas del proceso Capacidad máxima Costos unitarios
p precios $/metal
Mina M m gc Ley de corte %metal
gm Ley media %metal
Proceso H h
y Recuperación
Mercado K k d Costo de capital

x, proporción de mineral con respecto al

material removido

Flujo de caja C  ( p  k)  g  y  x  m  h  x  f t
 Formulación algoritmo de Lane
x, proporción de mineral con respecto al
x, t C, flujo de caja material removido
Flujo de caja

C  ( p  k )  g  y  x  m  h  x  f  t (1)
Valor V Valor
Presente Remanente
W Sea V el valor máximo de los flujos de
caja actualizado. W el valor máximo
de los flujos de caja actualizados
después de extraer Q
C W
V (2)
(1  d )t
Si t es pequeño, entonces
(1  d )t  1  d  t (3)
Formulación algoritmo de Lane
Incremento marginal del VAN 3 en 2
v  V  W  C  d V  t (4)

Incremento marginal del VAN 4 en1

v  ( p  k )  g  y  x  m  h  x  ( f  d V )  t
•Lo anterior representa el incremento marginal del VAN que conduce a
maximizarlo siguiendo una estrategia determinada

•El parámetro dV representa el costo de oportunidad de postergar los flujos

de caja futuros en el periodo t. Integrando el valor del dinero en el tiempo

•La formula anterior sin embargo considera un escenario económico

constante en el tiempo
Resolución de la formulación anterior
 El valor óptimo del VAN (V) es función de la
ley de corte
 La ley de corte a través de v es función del
valor óptimo del VAN (V)
 La manera de resolver este problema es a
través de múltiples iteraciones asumiendo un
valor V=0, en la primera iteración
Formulación algoritmo de Lane
Incorporando las condiciones variables del mercado en
la formula anterior
 dV 
v  ( p  k )  g  y  x  m  h  x  ( f  d  V   )t
 dT 
dV
F  d V 
dT

• F: corresponde al costo de oportunidad.

• Este costo se compone de dos partes :
• El costo de asignar capital a un bien en vez de invertir
en un bien alternativo (dV)
• La variación de valor de este bien debido a
condiciones variables en el mercado (-dV/dT)
Datos
 Costo mina 1.32$/t
 Costo planta 3.41 $/t
 Capacidad de mina 12 millones toneladas al año
 Capacidad de tratamiento 3.9 millones toneladas al
año
 Capacidad de venta 900 toneladas al año
 Costos fijos 11.9 millones al año
 Precio 60$/kg
 Recuperación 87%
 Costo de oportunidad 15.2 millones
Ley de corte (capacidad máxima mina)
0
 dV 
v  ( p  k )  g  y  x  m  h  x  ( f  d  V   )t
 dT 
1
t
( f  F) M
v  ( p  k)  g  y  x  m  h  x 
M h
gm 
x  f ( g ( gc )) ( p  k)  y
Max[ NPV ]  v  0

3.41
gm   0.07kg / t
60 0.87
Ley de corte (capacidad máxima
planta)
( f  F )x
v  ( p  k)  g  y  x  m  h  x 
H
x
t
H
( f  F)
Max[ NPV ]  v  0 ( h) 
gh  H
( p  k)  y

3.41  (11.9  15.2)

gh  3.9  0.2kg / t
60 0.87
Ley de corte una capacidad máxima
mercado
( f  F )g  y  x
v  ( p  k)  g  y  x  m  h  x 
H
gk  y  x
t
K
h
gk 
( f  F)
(pk  ) y
K

3.41
gk   0.13kg / t
(60 - (11.9  15.2) )0.87
0.9
COG de equilibrio (1)
COG de equilibrio (2)
COG de equilibrio (3)
COG de equilibrio (4)
Rajo Subterránea
Tomar en consideración que cuando se trata
de una mina subterránea todo el costo de
mina es considerado proceso ya que cada
tonelada extraída se procesa por lo tanto h
incluye el costo de mina
Ejemplo de aplicación
 Una mina mediana Plomo-Zinc con una capacidad de planta
de 1.3 millones de toneladas/año
 El cuerpo mineralizado es potente, semi vertical y consistente
en la corrida por varios kilómetros
 Método de explotación Sub Level Caving
 El valor de plomo es secundario y se convierte la ley de
plomo a Zinc dividiendo por 2
 Ley media Zinc equivalente 14%
 Las recursos son estimadas de esta manera y el control se
realiza basada solamente utilizando Zinc equivalente
 La dilución es del orden del 15% y la recuperación minera del
orden del 88%
Ejemplo de aplicación
 Existen dos tipos de leyes de corte
 Aplicada a las paredes del caserón. Define el
material a perforar
 Aplicada en el punto de extracción definiendo el
punto de extracción óptimo
 El potencial valor del proyecto es 150M$,
utilizando una tasa de descuento 15% antes de
impuestos
Datos necesarios

Perforación y tronadura $/t 3.6

Carguío $/t 1.5
Transporte $/t 3.2
Acarreo (pique) $/t 6.8
Molienda $/t 11.4

Total $/t 26.5

Costos fijos $/año 15.3

Capacidad de planta t/año 1.3
Recuperación metallúrgica % 81
Precio Zn $/conc. 650 Primer año
Precio Zn $/conc. 750 luego

Seguros y transporte $/t ZN 35

Resolución
 Tonelaje modelado debe ser corregido por
1.15 y la ley por 0.88/1.15=0.77 de modo de
simular la situación en la mina
 Sobre estos valores se deben calcular las leyes
de corte
Resolución
Cálculo de F para precios constantes

El cálculo de la ley de corte para una capacidad de planta fija

( f  F)
( h  m) 
gh  H
F  d V
( p  k)  y
F  0.15 150  22.5M $

El cálculo de la ley de corte para HW y FW

(15.3  22.5)
(26.5) 
1.3 55.6
gh    9.6%
(750  35)  0.81 579
Resolución
Ley de corte en los puntos de extracción. Perforación y tronadura es un
costo incurrido
Perforación y tronadura $/t 3.6
Carguío $/t 1.5 (15.3  22.5)
Transporte $/t 3.2 (22.9) 
1.3 52
Acarreo (pique) $/t 6.8 gh    9.0%
Molienda $/t 11.4 (750  35)  0.81 579
Total $/t 22.9

Los valores de leyes de corte son altos producto que el valor residual (150M$)
es alto, es decir, la oportunidad es alta. Si estuviésemos al fin de la vida de la
mina

(15.3)
(22.9) 
gh  1.3  34.7  6.0%
(750  35)  0.81 579
Resolución
Mas aún si la planta está ociosa la ley de corte a utilizar es aquella limitada
por la capacidad de la mina
(22.9)
gm   4.0%
(750  35)  0.81

Efecto del cambio de condiciones económicas o del mercado. Supongamos

que el precio del concentrado disminuye a 650 $. El cambio de utilidad anual
producto de la caída del precio es

(750  35) * 0.14 * 0.81*1.3  (650  35) * 0.14 * 0.81*1.3  14.7M $

Resolución
Valor presente del depósito a 650$

Valor presente al esperar un año para partir la explotación es 150 M$. Entonces

dV
 (150  137.2)  12.8M $
dT

dV 14.7
F  d V  150   137.2M $
dT (1  0.15)
F  0.15 137.2  12.8  7.8M $

La ley de corte de los puntos de extracción es ahora

(15.3  7.8) Producto de la baja de precios,

(22.9) 
gh  1.3 
40.7
 8.2%
la ley de corte debe bajar de
(650  35)  0.81 498 modo de ser consecuente con
el modelo económico
 Método iterativo para la optimización
del VAN
Closed dpts Active  Abrir nuevos recursos
t, $, gi  Calcular reservas utilizando leyes de corte
T, $ New  Cerrar reservas que no pagan ley de corte
 Simular la extracción del periodo
 Calcular el valor remanente del depósito
 Pasar al próximo periodo del plan
Remaining Value  Al final calcular vector de costos de
Vi oportunidad
 Adherir costos de oportunidad a la
estructura de costos
 Repetir el programa de producción
Iteración 1
Step Current tons $Value Cu Used-new Revenue Profit Wt Oct

1 183,050 15.9 1.996 3 2,919,365 2,919,365 10,730,724 5.86

2 183,050 16.0 1.997 0 2,921,512 2,921,512 8,882,285 4.85
3 183,050 7.6 1.300 0 1,391,415 1,391,415 8,379,099 4.58
4 183,050 2.5 0.877 1 462,398 462,398 8,754,610 4.78
5 183,050 5.8 1.153 1 1,067,322 1,067,322 8,562,749 4.68
6 183,050 7.9 1.327 1 1,451,375 1,451,375 7,967,649 4.35
7 183,050 7.5 1.288 0 1,364,235 1,364,235 7,400,179 4.04
8 183,050 4.4 1.030 0 799,055 799,055 7,341,142 4.01
9 183,050 3.0 0.913 0 540,757 540,757 7,534,498 4.12
10 183,050 5.1 1.092 2 934,640 934,640 7,353,308 4.02
11 183,050 7.8 1.318 0 1,430,197 1,430,197 6,658,442 3.64
12 270,000 14.48159 1.873 2 3,910,029 3,910,029 3,414,257 1.87
13 270,000 11.1 1.6 0 2,991,085 2,991,085 764,598 0.42
14 230,494 2.844505 0.904 0 655,643 655,643 185,415 0.10
15 81,625 2.49869 0.875 0 203,957 203,957 0 0.00

a 0.1
DV 0
Cap 183,050
NPV 12,409,171
Net DV costs $6.00
Iteración 2
Step Current tons $Value Cu Used-new Revenue Profit Wt Oct

1 183,050 15.9 1.996 3 2,919,365 2,919,365 11,816,720 6.46

2 183,050 16.0 1.997 0 2,921,512 2,921,512 10,076,881 5.50
3 183,050 12.8 1.734 2 2,343,786 2,343,786 8,740,783 4.78
4 183,050 7.9 1.327 1 1,451,375 1,451,375 8,163,486 4.46
5 183,050 7.5 1.292 0 1,374,377 1,374,377 7,605,457 4.15
6 183,050 7.5 1.293 1 1,375,122 1,375,122 6,990,881 3.82
7 183,050 10.2 1.520 1 1,874,756 1,874,756 5,815,213 3.18
8 183,050 13.9 1.825 1 2,543,706 2,543,706 3,853,028 2.10
9 183,050 15.9 1.993 1 2,913,460 2,913,460 1,324,871 0.72
10 74,355 19.6 2.300 0 1,457,358 1,457,358 0 0.00
1,721,805 1.692 0 0 0 0 0.00

a 0.1
DV 0
Cap 183,050
NPV 13,396,441
Net DV costs $6.88
Iteración 4
Step Current tons $Value Cu Used-new Revenue Profit Wt Oct

1 183,050 15.9 1.996 3 2,919,365 2,919,365 12,119,382 6.62

2 183,050 16.0 1.997 0 2,921,512 2,921,512 10,409,809 5.69
3 183,050 12.8 1.734 2 2,343,786 2,343,786 9,107,003 4.98
4 183,050 7.9 1.327 1 1,451,375 1,451,375 8,566,328 4.68
5 183,050 7.5 1.292 0 1,374,377 1,374,377 8,048,584 4.40
6 183,050 7.5 1.293 1 1,375,122 1,375,122 7,478,320 4.09
7 183,050 10.2 1.520 1 1,874,756 1,874,756 6,351,396 3.47
8 183,050 13.9 1.825 1 2,543,706 2,543,706 4,442,830 2.43
9 183,050 10.3 1.528 0 1,892,513 1,892,513 2,994,600 1.64
10 170,790 11.1 1.594 1 1,901,043 1,901,043 1,393,016 0.76
11 115,680 11.2 1.600 0 1,294,969 1,294,969 237,349 0.13
12 74,612 2.799999 0.900 0 208,914 208,914 52,170 0.03
13 20,495 2.8 0.9 0 57,387 57,387 0 0.00
2,029,028 1.577 0 0 0 0 0.00
0 0 0 0 0.00

a 0.1
DV 0
Cap 183,050
NPV 13,671,588
Net DV costs $6.84
Costo de oportunidad para distintas
iteraciones
7.0

6.0
Opportunity Cost ($/t)

5.0

4.0

3.0
Optimum Shut Off
2.0
Value Policy
1.0

0.0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
OC_1 OC_2 OC_3 OC_4 Period #
Programas de producción resultantes
2.50

2.00
Grade-Scheduled %Cu

1.50

1.00

0.50

0.00
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

Grades Cu_1 Grades Cu_2 Grades Cu_3 Grades Cu_4 Period #

Múltiples iteraciones y su
influencia en el programa de
producción
 Incorporación del componente dV/dT
Evalúa recursos
remanentes Vt Wt dV dV/dT Oct
utilizando precios
12,721,023 12,474,513 1,272,102 -246,510 8.30
del periodo t+1 9,979,814 10,051,010 997,981 71,196 5.06
8,103,285 8,181,601 810,328 78,316 4.00
7,024,005 7,565,354 702,401 541,348 0.88
6,310,133 6,905,616 631,013 595,483 0.19
5,808,757 6,863,295 580,876 1,054,538 -2.59
6,048,845 7,208,837 604,885 1,159,991 -3.03
6,661,790 8,070,175 666,179 1,408,385 -4.05
7,723,333 9,272,557 772,333 1,549,223 -4.24
7,877,226 8,920,520 787,723 1,043,294 -1.40
7,389,207 8,536,830 738,921 1,147,623 -2.23
4,460,392 4,807,741 446,039 347,349 0.54
1,335,172 1,335,172 133,517 0 0.73
364,283 364,283 36,428 0 0.20
0 0 0 0 0.00
Efecto de variación de precios en una
política de leyes de corte

4.0
16
Revenue Factor 3.5
14
Revenue Factor $/%Cu

3.0
12
2.5
10

Grade-%Cu
8 OC(without -dV/dt)_Grade 2.0

6 1.5

4 1.0
2 OC(with -dV/dT)-Grade 0.5
0 0.0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Period #
Ejercicios
 Suponga que la función de costo está dad por
C (q)  a  bq 2
 Y el precio está dado por
Pt  p
 Obtenga una relación para el costo de
oportunidad en el tiempo. Explique
Ejemplos de la industria
 Planes de LP
 Planes de CP
 Requerimiento de equipos