Obtención de Pit Final Óptimo
Obtención de Pit Final Óptimo
Obtención de Pit Final Óptimo
FACULTAD DE INGENIERA
DEPARTAMENTO DE INGENIERA EN MINAS
OBTENCIN DE
PIT FINAL PTIMO
MEDIANTE WHITTLE Y ALGORITMO
MILAWA
PREPARADO POR:
CAMILO CONTRERAS
VALESKA DAZ
REVISADO POR:
CRISTIAN VERA
19-12-2012
NDICE DE CONTENIDOS
OBJETIVOS ................................................................................................................................................... 2
ALGORITMO WHITTLE .................................................................................................................................. 3
ALGORITMO MILAWA .................................................................................................................................. 7
MODO NPV.............................................................................................................................................. 7
MODO BALANCEANDO ................................................................................................................................ 7
APLICACIN ................................................................................................................................................. 8
CONCLUSIONES .......................................................................................................................................... 16
|1
OBJETIVOS
Explicar claramente cmo se determina la optimizacin del pit final mediante Whittle y
algoritmo de Milawa
|2
ALGORITMO WHITTLE
Es una herramienta computacional especializada en la planificacin a nivel estratgico,
perteneciente a la familia de optimizadores de rajos que trabajan con el algoritmo de Lerchs &
Grossmann (3D y 4D).
La determinacin del pit final mediante este criterio se realiza por medio de la sensibilizacin
del precio, el cual es una de las variables ms influyente en la valorizacin que utiliza Whittle.
Esta sensibilizacin se logra con el Revenue Factor RF, el cual entrega un rango de valores que
permiten ir aumentando el precio, este aumento se puede determinar mediante delta RF y
definiendo el nmero de pit.
Adems, Whittle calcula los efectos de la variacin de los diferentes factores como costos,
precios, leyes o recuperaciones en la forma y tamao del pit.
La variacin que realiza Whittle est definida de la siguiente forma:
Donde:
- metal= Unidades de producto contenidas en la mena.
- recup= Fraccin del producto a recuperar.
- precio = Precio obtenible por unidad de producto.
- mena = Toneladas de mena en el bloque (para cualquier roca que contenga producto, independiente
de la ley).
- Cp = Costo de procesar una tonelada de mena + costo minero extra C, sobre o bajo el costo de extraer
el mineral como estril (C = costo ext. mena - costo ext. estril)
- roca = Toneladas de roca (estril y mineral) en el bloque.
- Cm = Costo de extraccin mina por tonelada de estril.
- Cv = Costo de refinacin por unidad de producto (debe incluir gastos de comercializacin del
producto)
Por lo tanto, por cada RF, se tendr un precio, y por cada precio un pit. Cada uno de estos pits
son ptimos por s solos para las condiciones dadas. De esta forma se crea una serie de pits
anidados.
|3
b) Best Case: Se considera la secuencia de extraccin pit a pit. Bajo el punto de vista
econmico es el mejor caso, pero operacionalmente no es viable debido a anchos
operacionales de equipos (pits muy cerca), ngulos de talud, factor de seguridad, etc.
|4
c) Real Case: Consiste en considerar tanto los criterios econmicos como operacionales.
Se toma el Van mayor del Best Case y el van mayor del Worst Case para acotar un
rango donde se encontrar el pit final.
|5
El nmero de fases depende de la vida del pit final escogido y de los mximos
horizontes planificados por fases.
Ya habiendo definido el Pit Final ptimo, slo resta definir parmetros para aplicar el
algoritmo de Milawa en Whittle y obtener la secuencia de extraccin de fases.
|6
ALGORITMO MILAWA
Este algoritmo es parte del mdulo de anlisis de Whittle. En
forma genrica busca una envolvente final con una adecuada
estrategia de consumo de reservas que entregue el mejor valor
presente neto y que adems asegure la vialidad del proyecto.
Se busca la mejor determinacin del pit final y una secuencia de explotacin que sern la base
de la planificacin a largo plazo de un rajo, siendo ac donde actan las herramientas
computacionales especializadas en la planificacin minera a nivel estratgico.
Si bien es cierto los pasos del algoritmo no son conocidos, se sabe que trabaja con dos
mdulos:
Modo NPV
Este modo busca un plan minero de extraccin de los bloques que optimice el valor
presente neto (VAN) del proyecto, sin considerar las unidades de la faena (mina-plantarefinera), por lo que no se preocupa de balancear la secuencia de extraccin de acuerdo a
las capacidades de cada una de dichas unidades. El orden en que los bancos de las fases son
explotados, es determinado por el algoritmo.
Modo Balanceando
Busca un plan que maximice el valor presente neto del proyecto, privilegiando el equilibrio
entre etapas (mina-planta-refinera). En este modo el algoritmo busca aumentar
tempranamente al mximo el uso de sus instalaciones de produccin en la vida de la mina,
en vez de aumentar al mximo el NPV.
Adems, el algoritmo de Milawa utiliza tres rutinas:
1) La primera toma un conjunto de variables y genera un posible programa entre ellas. El
nmero de variables depende del nmero de:
- Bancos en el Pit Final
- Fases
- Nmero de perodos en la vida de la Mina
2) La segunda rutina es una evaluacin la cual calcula el VPN de cada plan o programa
establecido.
3) La tercera rutina busca dentro de los diferentes programas el que posee el ms alto
VPN o el mejor balance entre etapas de produccin.
|7
APLICACIN
El primer paso para realizar la corrida Whittle es tener los datos de entrada que se utilizarn
en el anlisis. En este caso, se utilizaron los siguientes:
ngulo ()
48
Costo Mina
(US$/T min)
Costo proceso
(US$/T min)
Precio Cu (US$/lb)
Precio Mo (US$/lb)
Precio Venta Mo
(US$/lb)
Precio Venta Cu
(US$/lb)
Tasa descuento (%)
Capacidad Mina
(Mt/ao)
Capacidad Planta
(Mt/ao)
Recuperacin Mo
(%)
Recuperacin Cu
(%)
3.5
11.25
2.5
15
2.35
1.64
12
30
8
60
30
Modelo de Bloques:
Dimensin Bloques
20x20x18
Bloques X
75
Bloques Y
66
Bloques Z
64
Nmero Total
316800
Previamente, se haba realizado el estudio con un ngulo de 58. Sin embargo, se cambi a 48
debido a que el cono original presentaba muchas dificultades operacionales para lograr su
extraccin completa. Estas dificultades tenan relacin principalmente con la incorporacin
gigantesca de estril para llegar a las cotas inferiores, por lo que los flujos asociados al cono ya
no tenan relacin con el pit operativizado. Esto debido a la gran inclinacin que presenta con
|8
un ngulo empinado de 58. Segn esto, y luego de ingresar los parmetros al software, es
necesario realizar un anlisis de pits ptimos utilizando el Revenue Factor. Para esto se
necesita determinar el paso del factor, y hasta que valor se har la iteracin, es decir, el
mximo valor de RF hasta el que se iterar, y cada cunto se cambiar a otro pit. Cuando se
ingresan estos valores, adems de indicarle si se quieren pits anidados o fases de bancos, el
software realiza la iteracin, dando como resultado una tabla como sigue:
Pit
Rev
Factor
Rock
Ore
Strip
Max
Min
Bench Bench
Mo Units
Mo
Grade
Cu Units
Cu
Grade
0.96
115735000
28933000
54
29
3737826
0.13
27927990
0.97
0.98
122298000
31396000
2.9
54
29
3991194
0.13
29899692
0.95
144262000
38102000
2.79
54
28
4717440
0.12
35688600
0.94
1.02
166658000
45085000
2.7
54
26
5440932
0.12
41670288
0.92
1.04
687761000
227950000
2.02
55
15
23143320
0.1
180628056
0.79
1.06
731988000
246224000
1.97
55
15
24682644
0.1
193841262
0.79
1.08
761715000
258795000
1.94
55
15
25701300
0.1
202747374
0.78
1.1
795167000
270751000
1.94
55
15
26702298
0.1
211930506
0.78
1.12
831327000
285379000
1.91
55
15
27829656
0.1
222348402
0.78
10
1.14
858115000
296606000
1.89
55
14
28647594
0.1
230159232
0.78
11
1.16
885065000
312239000
1.83
55
14
29635794
0.09
239479902
0.77
12
1.18
908227000
324227000
1.8
55
14
30412422
0.09
246628800
0.76
13
1.2
942600000
337106000
1.8
55
14
31341330
0.09
255224844
0.76
14
1.22
951418000
341966000
1.78
55
14
31634712
0.09
257869332
0.75
15
1.24
997676000
356189000
1.8
55
14
32701320
0.09
267978294
0.75
16
1.26
1029948000 365666000
1.82
55
13
33402780
0.09
274723974
0.75
17
1.28
1039624000 368485000
1.82
55
13
33609978
0.09
276675426
0.75
18
1.3
1051520000 372989000
1.82
55
13
33887484
0.09
279527760
0.75
19
1.32
1075989000 382126000
1.82
55
13
34442172
0.09
285241824
0.75
20
1.34
1087842000 388670000
1.8
55
13
34770870
0.09
288713970
0.74
En esta tabla, podemos observar que el pit final es el pit N3, debido a que en este valor el
Revenue Factor es 1. Los pits que se encuentren bajo el RF (en este caso, pits 1 y 2), podran
ser considerados como fases de explotacin del pit 3.
|9
Una vez obtenido el Pit Final con el RF, se procede a determinar una secuencia de extraccin,
con sus flujos de caja asociados. Para realizar esto, es necesario especificar el caso que ser
generado: Best, Wost o Specified. que se ve representado por la siguiente imagen: (curvas
azul, roja y verde respectivamente)
Luego, se elige la opcin para la construccin del plan de produccin, la que puede ser Fixed
Lead o Milawa (Balanced y NPV). Adems, es necesario definir los pushbacks para cada pit, en
el sentido de cuantos bancos de un pit como mnimo y mximo puede bajar en un perodo
(slo para Milawa NPV o Balanceado). Cuando se definen estos parmetros, el resultado es el
siguiente:
Final
pit
Open pit
cashflow
best
$ disc
Open pit
cashflow
specified
$ disc
Open pit
cashflow
worst
$ disc
tonne
input
best
Waste
best
tonne
Mine
life
years
best
Mine
life
years
specified
Mine
life
years
worst
Internal
rate of
return
best %
1
2
3
-36505129
-35786788
-35553840
-36505129
-46324780
-71229176
-36505129
-46324780
-71229176
29371000
31574000
38102000
86364000
90724000
106160000
5.93563
6.21103
7.0271
5.93563
6.50111
7.70522
5.93563
6.50111
7.70522
0
0
0
Internal
rate of
return
specified
%
0
0
0
| 10
Internal
rate of
return
worst
%
0
0
0
Como se observa en la tabla, los 3 pits a los que se les realiz el anlisis presentan flujos de
caja descontados menores que 0. Por lo tanto, si el proyecto se evala en este momento, no
flotara, ya que de ningn cono se extraen beneficios, para los parmetros dados. Entonces,
para que el anlisis que sigue tenga algn sentido, se realiza una sensibilizacin del precio de
venta del Mo. Si el precio se baja, obviamente que los flujos que entregar el software seguirn
siendo negativos, ya que a ese precio ser mucha menor la ganancia obtenida. Por lo tanto, es
necesario realizar un aumento del precio del Mo. Se escoge el precio 17 US$/lb, ya que en este
precio se obtiene 1 pit con flujo descontado positivo con un valor suficientemente grande
como para trabajarlo. Tambin se analiz el caso de 18 US$/lb, para determinar el
comportamiento para un valor de venta relativamente pequeo.
Rev
Factor
Rock
Ore
Strip
Max
Bench
Min
Bench
Mo
Units
Mo
Grade
Cu
Units
Cu
Grade
0.9
122298000
31331000
2.9
54
29
3988440
0.13
29842506
0.95
0.92
133673000
34798000
2.84
54
28
4369302
0.13
32817474
0.94
0.94
161492000
43303000
2.73
54
26
5268240
0.12
40227030
0.93
0.96
653241000
215492000
2.03
55
15
22005918
0.1
170749944
0.79
0.98
711612000
238448000
1.98
55
15
24019740
0.1
187785216
0.79
760917000
257450000
1.96
55
15
25634394
0.1
201911130
0.78
1.02
789669000
267559000
1.95
55
15
26492670
0.1
209725848
0.78
1.04
830549000
282852000
1.94
55
15
27719172
0.1
221011254
0.78
1.06
855334000
294500000
1.9
55
14
28527228
0.1
228900006
0.78
10
1.08
887116000
312514000
1.84
55
14
29676942
0.1
239717880
0.77
11
1.1
919540000
327985000
1.8
55
14
30706128
0.09
249162642
0.76
12
1.12
942681000
338548000
1.78
55
14
31402080
0.09
255649932
0.76
13
1.14
973490000
348284000
1.8
55
14
32139990
0.09
262277028
0.75
14
1.16
1005056000
357842000
1.81
55
14
32851008
0.09
269219214
0.75
15
1.18
1037391000
367060000
1.83
55
13
33541938
0.09
275807268
0.75
16
1.2
1048335000
370105000
1.83
55
13
33764526
0.09
278004312
0.75
17
1.22
1064089000
375824000
1.83
55
13
34118334
0.09
281685114
0.75
18
1.24
1084041000
384199000
1.82
55
13
34597044
0.09
286648470
0.75
19
1.26
1102779000
394924000
1.79
55
13
35115606
0.09
292328676
0.74
20
1.28
1121406000
404174000
1.77
55
13
35583138
0.09
297186570
0.74
| 11
Para este caso, el RF =1 se encuentra en el pit 6, con un tonelaje de mineral y total mucho
mayor que el realizado para un precio de 15 US$.
En el anlisis de los casos, se observa que a partir del pit 3, el valor comienza a disminuir, por
lo que es en este punto que el proyecto ya no gana ms valor.
Realizando la secuencia de extraccin, el resultado fue:
Final
pit
1
2
3
4
5
6
7
8
9
10
Open pit
cashflow
best
$ disc
11666581
14879737
23356460
-162857698
-162954842
-163111041
-163150130
-163285701
-16356498
-163451080
Open pit
cashflow
specified
$ disc
11666581
-1400606
-21549299
-543873823
-578772257
-605327153
-619240173
-638108775
-650040087
-662626002
Open pit
cashflow
worst
$ disc
11666581
-1400606
-21549299
-543873823
-578772257
-605327153
-619240173
-638108775
-650040087
-662626002
tonne
input
best
Waste
best
tonne
32562000
36304000
44825000
220774000
240376000
257450000
266555000
279272000
286076000
295180000
89736000
97368000
116667000
432468000
471236000
503467000
523114000
551277000
569258000
591936000
Mine
life
years
best
6.6054
7-07318
8-13833
37.68533
40.45025
43.10147
44.23952
46.29077
47.14127
48.34281
Mine
life
years
specified
6.6054
7.36125
8.60318
37.9613
41.18705
43.96661
45.45416
47.5971
48.79025
50.33174
Mine
life
years
worst
6.6054
7.36125
8.60318
37.9613
41.18705
43.96661
45.45416
47.5971
48.79025
50.33174
Internal
rate of
return
best %
0
0
0
0
0
0
0
0
0
0
Internal
rate of
return
specified %
0
0
0
0
0
0
0
0
0
0
Como se ve en la tabla, se obtuvo un pit con flujos de caja descontados positivos. Por lo tanto,
como no es suficiente, se realiz el anlisis con 18 US$/lb
| 12
Internal
rate of
return
worst %
0
0
0
0
0
0
0
0
0
0
En la imagen, se observa que en el pit 3 es donde existe la inflexin de las curvas de best case,
specified case y worst case. Esto, coincide con la lista de flujos de caja descontados, ya que son
estos pits los que presentan flujos de caja positivos.
Final
pit
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
Open pit
cashflow
best
$ disc
50542136
64167095
75809348
-86496461
-86561561
-86529848
-86530169
-86576582
-86600907
-86671935
-86726655
-86824371
-86878381
-86972469
-87062351
-87089164
-87160940
-87185174
-87233014
-87283242
Open pit
cashflow
specified
$ disc
50542136
29089920
12145887
-497022829
-556220600
-580430242
-597885774
-617854572
-628923636
-645301829
-654132541
-668152166
-676081256
-689687394
-700877748
-704504823
-712482994
-714990437
-720750836
-726789767
Open pit
cashflow
worst
$ disc
50542136
29089920
12145887
-497022829
-556220600
-580430242
-597885774
-617854572
-628923636
-645301829
-654132541
-668152167
-676081256
-689687394
-700877748
-704504823
-712482993
-714990437
-720750836
-726789767
tonne
input
best
Waste
best
tonne
30407000
36531000
45085000
209741000
241396000
256122000
267057000
279839000
286141000
297027000
303734000
312903000
319043000
329103000
338013000
340346000
346907000
349304000
353646000
358004000
85328000
97142000
116408000
408102000
468110000
495805000
517833000
546065000
562172000
588704000
604477000
629746000
646020000
673971000
699378000
707990000
727597000
734656000
749133000
765930000
Mine
life
years
best
6.02
6.78545
8.0692
35.80173
40.66196
42.75531
44.12219
46.23747
47.0252
48.83427
49.67262
5108083
51.8483
53.25753
54.65971
54.95131
55.82013
56.11982
56.67144
57.54472
Mine
life
years
specified
6.02
7.3896
8.63558
36.11687
41.27352
43.64145
45.44814
47.5784
48.68945
50.51583
51.66391
53.24394
54.31189
56.06223
57.63096
58.08142
59.25412
59.65956
60.47905
61.32961
Mine
life
years
worst
6.02
7.3896
8.63558
36.11688
41.27352
43.64145
45.44814
47.5784
48.68945
50.51583
51.66391
53.24394
54.31189
56.06223
57.63096
58.08142
59.25412
59.65956
60.47905
61.32961
Internal
rate of
return
best %
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Internal
rate of
return
specified %
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
| 13
Internal
rate of
return
worst %
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
En este caso se observa que los flujos de caja descontados, asociados a los pits son positivos
hasta el pit 3. Sin embargo, es el pit 1 el que presenta el mayor flujo de caja, por lo que ste se
elige como pit final, y es ste el que se trabaja.
Ahora bien, como se tiene un cono con flujos de caja descontados positivos, se le pueden
realizar los anlisis de secuencia de extraccin como se mencionaron anteriormente.
Entonces, en este caso especfico se entrega el anlisis realizado para el precio de 18 US$/lb
Mo. De acuerdo a eso, los valores obtenidos, utilizando Milawa Balanceado, son los siguientes:
Period
tonne
input
Waste
tonne
Strip
ratio
Grade
input
Mo
Grade
input
Cu
30000000
999.99
30000000
999.99
3
4
5
6
7
6247000
8000000
8000000
8000000
160000
23753000
1199000
367000
8000
0
3.8
0.15
0.05
0
0
0.1051
0.1142
0.1258
0.1502
0.1949
0.5424
0.8895
1.0706
1.1687
1.1829
Open pit
Open pit
cashflow
cashflow
$
$ disc
-93750000
105000000
-83705357
105000000
-20101322 -14307724
107458355 68291727
137814647 78199732
183958670 93199187
5172411
2.614571
| 14
Como se observa en el grfico y en la tabla, se tiene un plan de produccin dado para 7 aos, el
cual considera 2 perodos para la extraccin exclusiva de estril, para en el tercer perodo
comenzar la extraccin de mineral. Adems, considera un vector de leyes de Mo creciente,
como lo indica la curva roja.
Por lo tanto, se establece un plan de produccin para el cono 1, asociado a un precio de Mo de
18US$/lb, utilizando todos los factores antes mencionados.
| 15
CONCLUSIONES
El proyecto dado, utilizando un ngulo de 48 no flota, debido a que ningn pit obtuvo un flujo
de caja descontado positivo en la corrida whittle. Sin embargo, si se vara el precio del
molibdeno (elemento a considerar), el proyecto puede flotar, y se pueden realizar todos los
dems anlisis. Esta variacin fue slo demostrativa, no est relacionada a ninguna
proyeccin de precios realizada por alguna institucin competente. Sin embargo, los precios
histricos del molibdeno dan la posibilidad de que estos precios se presenten, pero sin tener
una certeza de ello.
Por otra parte, luego de esta sensibilizacin, el proyecto es capaz de flotar, entregando un plan
de produccin, asociado a un pit que tiene un flujo de caja descontado de US$ 50,542,136 y
utilizando Milawa Balanceado, se obtiene un plan de produccin, el cual muestra la extraccin
de estril en los primeros perodos, para luego comenzar con la extraccin de mineral,
sumado con un vector de leyes creciente, tanto para el Mo, como para el Cu.
| 16