UNIVERSIDAD DE CHILE FACULTAD DE CIENCIAS FÍSICAS Y MATEMÁTICAS DEPARTAMENTO DE INGENIERÍA DE MINAS
MODELO DE PROGRAMACIÓN MATEMÁTICA PARA SUSTENTAR LA TRANSICIÓN RAJO SUBTERRÁNEA.
TESIS PARA OPTAR AL GRADO DE MAGÍSTER EN MINERÍA MEMORIA PARA OPTAR AL TÍTULO DE INGENIERO CIVIL DE MINAS
ANDRÉS ALEJANDRO SOLAR DROGUETT
PROFESOR GUÍA: ENRIQUE RUBIO ESQUIVEL MIEMBROS DE LA COMISIÓN: XAVIER EMERY DANIEL ESPINOZA GONZÁLEZ MAURICIO LARRAÍN MEDINA
SANTIAGO DE CHILE AGOSTO 2010
Caballeros, ¿Debo recordarles que mis probabilidades de éxito aumentan en cada nuevo intento? John Nash
II
Caballeros, ¿Debo recordarles que mis probabilidades de éxito aumentan en cada nuevo intento? John Nash
II
Resumen. Varias de las grandes operaciones a rajo abierto se están enfrentando al desafío de diseñar y planificar una transición del método de explotación para sostener el futuro a largo plazo de sus operaciones. El método subterráneo deseado es el Block Caving, dados los altos ritmos de producción y bajos costos de operación, que pueden competir contra los costos asociados a la extracción de las últimas fases de un rajo. De este modo, existe una pérdida de valor al explotar las últimas fases del rajo por el método de superficie en vez de extraerlas por un método subterráneo, con la consiguiente reducción de la cantidad de estéril que debe ser removida. Lo propuesto al momento de planificar una transición, corresponde a definir el fondo del rajo sujeto a una opción subterránea, es decir, una optimización simultánea entre el rajo y la subterránea. Tradicionalmente, la forma de abordar la transición rajo subterránea ha sido a través de un método que consiste en decidir secuencialmente el pit final usando la optimización de Lerchs y Grossmann para encontrar el máximo pit , y luego correr simulaciones para encontrar el mejor piso para la opción subterránea. Esta investigación está relacionada con el desarrollo de un modelo de programación lineal entero, que apunta a calcular simultáneamente el máximo VAN de la opción combinada entre el pit final del rajo y el piso de la subterránea. El modelo utiliza pits previamente diseñados a diferentes leyes de corte de manera de discretizar el inventario de recurso de mineral. Este conjunto inicial se preocupa de la geometría y los principales efectos sobre la estabilidad de las paredes y además del área mínima del piso. El modelo ha sido probado usando dos yacimientos, Palabora y Andina Sur-Sur, representando dos de las actuales operaciones en el proceso de transición de una operación rajo a una subterránea y también dos tipos de mineralización distintos. Tanto los resultados obtenidos en Palabora como en Sur-Sur indican que el modelo de optimización permite definir un diseño que genera una mejora en el VAN del proyecto, sujeto principalmente a un cambio en la profundidad final del rajo. Esto es, las últimas fases o incrementos del rajo son extraídos mediante el Block Caving , de manera que generan mejores beneficios. Esto, además, permite que el valor de la mina subterránea se vea afectado de menor manera por la acción de la tasa de descuento. Con la metodología propuesta, en el caso de Palabora se logró una mejora en el VAN del proyecto de hasta un 7%, mientras que en el caso de Andina Sur-Sur se alcanzó una mejora de un 16%. En consecuencia, la aplicabilidad del modelo de optimización es probada en dos tipos de yacimientos distintos, distintos, mejorando notablemente el valor obtenido por una aproximación tradicional. Asimismo, permite generar un análisis de sensibilidad asociado a parámetros críticos vinculados con el proyecto.
III
Abstract. Several large open pit mines are facing the challenge to design and plan an underground transition to sustain the long term future of their operations. The desired underground method is Block Cave because its high productivity and low operating cost that can compete against the open pit cost to mine the final open pit pushbacks. It has been found at several operations where this transition is taking place that the pit limit at which transition is performed is rather deep. Thus, there is a substantial loss in value by mining the final open pit pushbacks by surface methods instead of mining them by underground methods, reducing the amount of waste to be extracted. This suggests that several transition projects may have started late in the life of an operation, producing an over stress in the transition project and reducing the flexibility in the design of the overall infrastructure. Traditionally the approach to plan open pit underground transition has been taken over by a method consisting of sequentially deciding the final pit using a Lerchs and Grossmann optimization to find the overall maximum pit shell and then running block cave simulations to find the best footprint for the underground option. The research summarized in this thesis is related to the development of an integer linear programming model that aims to compute simultaneously the maximum VAN of the combined option final limit open pit and block cave underground footprint. The model uses predesigned pit shells at different cut-off grades in order to discretize the mineral resource inventory. This initial aggregation takes care of geometries and main effects on wall stabilities and minimum underground footprint. The results obtained in Palabora and Sur-Sur, show that the integer lineal program allows to achieve a better NPV, due to define another transition design which proposes a less final pit depth. The final pushbacks or open pit marginal increments are mined by block caving, this change in the mining method improves the profit obtained. Furthermore, the discount rate has not a high impact over the value of underground mine. The proposed methodology makes an improvement on the NPV of 7% in Palabora, meanwhile; in Andina Sur-Sur the improvement achieve an 16%. Finally, the applicability of the integer linear program has been proved in two different ore bodies, and makes an improvement in the NPV of the project compare with the NPV of a traditional method. In addition, the model allows to makes an easy and quick study of different parameters concerned to the open pit underground transition.
IV
Tabla de Abreviaturas. col US$ tpa tpd lb t VAN m m3 LC
Columnas. Dólares americanos. Toneladas por año. Toneladas por día. Libras. Toneladas. Valor actual neto. Metros. Metros cúbicos. Ley de corte. Existencia. S.R. Sin Restricción. HOD Altura de columna extraíble (Height of draw) dpt Punto de extracción (Drawpoints)
V
Agradecimientos. Me gustaría comenzar agradeciendo especialmente a mis padres, Orlando y Aída, quienes siempre se han preocupado y esmerado en darme lo que he necesitado y me han apoyando en todo momento, también a mi hermano Christian, Carolina, Francisca y Cristóbal, que siempre están pendientes y me brindan su apoyo. No puedo dejar de mencionar también a mi profesor guía, Enrique Rubio, quien me permitió realizar éste trabajo y me dio la oportunidad de crecer tanto profesionalmente como personalmente, con cada uno de sus consejos, anécdotas o datos útiles. También aprovecho la oportunidad de dar las gracias a Xavier Emery, Daniel Espinoza y Mauricio Larraín, quienes aceptaron participar en el desarrollo de éste trabajo y que me permitieron mejorarlo e incorporar nuevas perspectivas e ideas. Además quisiera agradecer a mis compañeros del Laboratorio de Planificación “DELPHOS”, Sebastián y Camilo, por su buena onda, apoyo y con los cuales el trabajo durante los meses de verano se hizo menos difícil de sobrellevar. Así como también a mis amigos del Departamento de Ingeniería de Minas y de otras carreras, que hicieron que estos años de Universidad fueran menos duros. Quiero agradecer también, a Claudia, por tu constante ayuda, cariño, comprensión y por acompañarme durante ésta etapa de mi vida.
VI
Índice de contenidos. 1.
Introducción. ______________________________________________________________ 1
1.1.
Objetivos. __________________________________________________________________ 4
1.2.
Alcances.___________________________________________________________________ 5
1.3.
Metodología de trabajo. ______________________________________________________ 6
1.3.1. 1.3.2. 1.3.3.
2.
3.
4.
Método secuencial. _________________________________________________________________ 6 Método incremental por fases. _______________________________________________________ 7 Método de transición ILP. ____________________________________________________________ 9
Revisión bibliográfica. ______________________________________________________ 9
2.1.
Transición de minería a rajo abierto a subterránea. ________________________________ 9
2.2.
Planificación en minas a rajo abierto. ___________________________________________ 13
2.3.
Planificación de minas subterráneas de block caving. ______________________________ 14
2.4.
Programación matemática asociada a la transición rajo subterránea. _________________ 15
Modelo de optimización para la transición rajo – subterránea. ____________________ 17
3.1.
Supuestos. ________________________________________________________________ 17
3.2.
Generación de las unidades de reserva para el modelo de optimización. ______________ 17
3.3.
Descripción del modelo. _____________________________________________________ 19
Información disponible para el caso de estudio de Palabora. _____________________ 255
4.1.
Inventario de recursos. ______________________________________________________ 25
4.2.
Parámetros técnico-económicos para Palabora. _________________________________ 255
4.3.
Datos de entrada para el modelo de optimización. ________________________________ 27
5. Análisis de la transición para el caso de estudio de Palabora. ______________________ 28
6.
5.1.
Resultados para la metodología secuencial. ______________________________________ 28
5.2.
Resultados para la metodología incremental por fases. ____________________________ 31
5.3.
Resultados para la metodología de transición ILP en Palabora. ______________________ 35
5.4.
Comparación entre las alternativas de transición para Palabora. _____________________ 52
Información disponible para el caso de aplicación en Andina Sur-Sur. __________ _____ 54
6.1.
Inventario de recursos. ______________________________________________________ 54
6.2.
Parámetros técnico-económicos para Andina Sur-Sur. _____________________________ 54
6.3.
Datos de entrada para el modelo de optimización. ________________________________ 55
7. Análisis de la transición para el caso de aplicación en Andina Sur-Sur. _______________ 57
7.1.
Resultados para la metodología secuencial. ______________________________________ 57
7.2.
Resultados para la metodología incremental por fases. ____________________________ 61
7.3.
Resultados para la metodología de transición ILP en Andina Sur-Sur. _________________ 65
VII
7.4.
Comparación entre las alternativas de transición para Andina Sur-Sur. ________________ 71
8.
Conclusiones y Recomendaciones. ____________________________________________ 73
9.
Bibliografía.______________________________________________________________ 77
10. Anexos. _________________________________________________________________ 79
10.1.
Flujos de caja descontados para los estratos extraídos por rajo y subterránea. Caso Palabora. _________________________________________________________________ 79
10.2.
Flujos de caja descontados para los estratos extraídos por rajo y subterránea. Caso Sur Sur. _________________________________________________________________________ 80
10.3.
Valorización de estratos ILP, extraídos vía rajo. Caso Palabora. ______________________ 81
10.3.1. 10.3.2. 10.3.3.
10.4.
Max Pit. _______________________________________________________________________ 81 Med Pit._______________________________________________________________________ 86 Min Pit. _______________________________________________________________________ 91
Valorización de estratos ILP, extraídos vía subterránea. Caso Palabora. _______________ 96
10.4.1. 10.4.2.
Discretización A. ________________________________________________________________ 96 Discretización B. ________________________________________________________________ 99
10.5.
Valorización de estratos ILP, extraídos vía rajo. Caso Andina Sur-Sur. ________________ 103
10.6.
Valorización de estratos ILP, extraídos vía subterránea. Caso Andina Sur-Sur. _________ 108
10.7.
Paper presentado en APCOM 2009, Vancouver, Canadá. __________________________ 113
VIII
Índice de tablas. Tabla 1. Valores de las variables para cumplir el criterio de marginalidad de la transición. __________ 3 Tabla 2. Parámetros económicos para el rajo. Caso Palabora. ________________________________ 26 Tabla 3. Parámetros económicos para el block caving. Caso Palabora. _________________________ 26 Tabla 4. Valorización de pits Whittle al ritmo de producción óptimo. Caso Palabora. _____________ 29 Tabla 5. Características del pit final elegido para la metodología secuencial en Palabora. _________ 29 Tabla 6. Corrida de footprint finder para Palabora. ________________________________________ 29 Tabla 7. Resultados para la metodología secuencial en Palabora. _____________________________ 30 Tabla 8. Estandarización de cotas. ______________________________________________________ 31 Tabla 9. Cota del fondo de los pits Whittle para Palabora. ___________________________________ 31 Tabla 10. Valorización de los estratos extraidos vía rajo. Caso Palabora. _______________________ 32 Tabla 11. Ley y mineral para cada estrato de Palabora extraíbles vía subterránea. _______________ 32 Tabla 12. Valorización de los estratos extraídos vía subterránea. Caso Palabora. ________________ 33 Tabla 13. Flujos de caja descontados para los estratos extraídos por rajo y subterránea. Caso Palabora. _______________________________________________________________ 33 Tabla 14. Resultados para la metodología incremental por fases en Palabora. __________________ 34 Tabla 15. Parámetros de evaluación de los estratos para el ILP extraíbles vía rajo. Caso Palabora. __ 35 Tabla 16. Valorización de los estratos para el ILP extraíbles vía rajo. Caso Palabora.______________ 36 Tabla 17. Parámetros de evaluación de los estratos para el ILP extraíbles vía subterránea. Caso Palabora. _______________________________________________________________ 36 Tabla 18. Costo de inversiones para la subterránea. Caso Palabora. ___________________________ 36 Tabla 19. Valorización de los estratos para el ILP extraíbles vía subterránea. Caso Palabora. _______ 37 Tabla 20. Vector de parámetros de inicio del minado del estrato vía rajo. Caso Palabora.__________ 37 Tabla 21. Resultados del VAN para ley de corte fija predefinida y dos discretizaciones de cilindros distintos (A y B). Caso Palabora. _________________________________________________ 38 Tabla 22. Resultados de la geometría de transición para ley de corte fija predefinida y dos discretizaciones de cilindros distintos (A y B). Caso Palabora. __________________________ 39 Tabla 23. Resultados del VAN para ley de corte fija y dos discretizaciones de cilindros distintos (A y B). Caso Palabora. _______________________________________________________________ 40 Tabla 24. Resultados de la geometría de transición para ley de corte fija y dos discretizaciones de cilindros distintos (A y B). Caso Palabora. __________________________________________ 41 Tabla 25. Resultados del VAN para ley de corte variable y dos discretizaciones de cilindros distintos (A y B). Caso Palabora. _________________________________________________________ 42 Tabla 26. Política de leyes de corte para la extracción de los estratos. Discretización MAX PIT. _____ 43 Tabla 27. Política de leyes de corte para la extracción de los estratos. Discretización MED PIT ______ 44 Tabla 28. Política de leyes de corte para la extracción de los estratos. Discretización MIN PIT ______ 45 Tabla 29. Resultados de la geometría de transición para ley de corte variable y dos discretizaciones de cilindros distintos (A y B). Caso Palabora. __________________________________________ 46 Tabla 30. Resultados con respecto al VAN para la corrida del modelo de optimización. Caso Palabora. _______________________________________________________________ 47 Tabla 31. Variación máxima porcentual del VAN debido al cambio en la discretización del rajo. Caso Palabora. _______________________________________________________________ 47 Tabla 32. Variación máxima porcentual del VAN debido al cambio en la altura de columna máxima extraíble. Caso Palabora _______________________________________________________ 47 Tabla 33. Resultados respecto a la geometría de la transición obtenidos a partir del modelo de optimización. Caso Palabora. ____________________________________________________ 48 Tabla 34. Máxima diferencia de la cota del fondo de pit sujeto a un cambio en la discretización del rajo. Caso Palabora. _______________________________________________________________ 48
IX
Tabla 35. Máxima diferencia de la cota del footprint sujeto a un cambio en la discretización del rajo. Caso Palabora. _______________________________________________________________ 48 Tabla 36. Máxima diferencia de la cota del fondo de pit sujeto a un cambio en la altura de columna máxima extraíble. Caso Palabora. ________________________________________________ 49 Tabla 37. Máxima diferencia de la cota del footprint sujeto a un cambio en la altura de columna máxima extraíble. Caso Palabora ________________________________________________ 49 Tabla 38. Análisis de sensibilidad para el VAN con respecto a la variación del precio. Caso Palabora. 49 Tabla 39. Análisis de sensibilidad para la geometría de transición con respecto a la variación del precio. Caso Palabora. _______________________________________________________________ 50 Tabla 40. Análisis de sensibilidad para el VAN con respecto a la variación del ritmo de producción del rajo. Caso Palabora. ___________________________________________________________ 50 Tabla 41. Análisis de sensibilidad para la geometría de transición con respecto a la variación del ritmo del producción del rajo. Caso Palabora. ___________________________________________ 50 Tabla 42. Análisis de sensibilidad para el VAN con respecto a la variación del ritmo de producción del block caving. Caso Palabora. ____________________________________________________ 50 Tabla 43. Análisis de sensibilidad para la geometría de transición con respecto a la variación del ritmo del producción del block caving. Caso Palabora. _____________________________________ 51 Tabla 44. Comparación entre metodologías según VAN. Caso Palabora.________________________ 52 Tabla 45. Comparación entre metodologías según el diseño de la transición. Caso Palabora. _______ 52 Tabla 46. Valor económico en competencia entre el rajo y la subterránea. Caso Palabora. _________ 53 Tabla 47. Parámetros económicos para el rajo. Caso Andina Sur-Sur. __________________________ 55 Tabla 48. Parámetros económicos para el block caving. Caso Andina Sur-Sur. ___________________ 55 Tabla 49. Valorización de los estratos extraídos vía rajo. Caso Andina Sur-Sur. __________________ 58 Tabla 50. Características del pit final elegido para la metodología secuencial en Andina Sur-Sur. ___ 58 Tabla 51. Corrida de footprint finder para Andina Sur-Sur. ___________________________________ 59 Tabla 52. Resultados para la metodología secuencial en Andina Sur-Sur. _______________________ 60 Tabla 53. Fondo de pits Whittle. Caso Andina Sur-Sur. ______________________________________ 61 Tabla 54. Cota estandarizada para los estratos extraíbles vía rajo. Caso Andina Sur-Sur. __________ 61 Tabla 55. Valorización de los estratos extraídos vía rajo. Caso Andina Sur-Sur. __________________ 61 Tabla 56. Corrida de footprint finder para Andina Sur-Sur. ___________________________________ 62 Tabla 57. Valorización de los estratos extraídos vía subterránea. Caso Sur Sur. __________________ 62 Tabla 58. Flujos de caja descontados para los estratos extraídos por rajo y subterránea. Caso Andina Sur-Sur. __________________________________________________________ 63 Tabla 59. Resultados para la metodología incremental por fases en Andina Sur-Sur. ______________ 63 Tabla 60. Parámetros de evaluación de los estratos para el I LP extraíbles vía rajo. Caso Andina Sur-Sur. __________________________________________________________ 66 Tabla 61. Parámetros de evaluación de los estratos para el ILP extraíbles vía subterránea. Caso Andina Sur-Sur. __________________________________________________________ 66 Tabla 62. Costo de inversiones para la subterránea. Caso Andina Sur-Sur. ______________________ 66 Tabla 63. Valorización de los estratos para el ILP extraíbles vía rajo. Caso Andina Sur-Sur. _________ 67 Tabla 64. Valorización de los estratos para el ILP extraíbles vía subterránea. Caso Andina Sur-Sur. __ 68 Tabla 65. Resultados del VAN para ley de corte variable y la discretización de la subterránea tipo A. Caso Andina Sur-Sur. __________________________________________________________ 69 Tabla 66. Resultados de la geometría de transición para ley de corte variable y la discretización de la subterránea tipo A. Caso Andina Sur-Sur. __________________________________________ 69 Tabla 67. Política de leyes de corte para la extracción de los estratos. Discretización MIN PIT. ______ 70 Tabla 68. Comparación entre metodologías según VAN. Caso Andina Sur-Sur. ___________________ 71 Tabla 69. Comparación entre metodologías según el diseño de la transición. Caso Andina Sur-Sur. __ 71 Tabla 70. Valor económico en competencia entre el rajo y la subt erránea. Caso Andina Sur-Sur. ____ 71 Tabla 71. Tabla comparación de resultados para la transición. Caso Palabora. __________________ 73
X
Tabla 72. Valor económico en competencia para el caso de Palabora. _________________________ 75 Tabla 73. Tabla comparación de resultados para la transición. Caso Sur Sur. ____________________ 75
Índice de figuras. Figura 1. Unidades de reserva que discretizan la explotación a rajo abierto. ____________________ 17 Figura 2. Vista en planta de cilindros que discretizan la explotación subterránea. ________________ 18 Figura 3. Esquema del diseño de la transición según la metodología secuencial para Palabora. _____ 30 Figura 4. Esquema del diseño de la transición según la metodología incremental por fases para Palabora. __________________________________________________________________________ 34 Figura 5. Esquema del diseño de la transición según la metodología secuencial para Andina Sur-Sur. 60 Figura 6. Esquema del diseño de la transición según la metodología incremental por fases para el Sur Sur _______________________________________________________________________________ 64
Índice de gráficos. Gráfico 1. Curvas de VAN incrementales. __________________________________________________ 8 Gráfico 2. Curva tonelaje - ley para Palabora. _____________________________________________ 25 Gráfico 3. Valorización de pits Whittle. Caso Palabora. _____________________________________ 28 Gráfico 4. Curva tonelaje – ley para Andina Sur-Sur. ________________________________________ 54 Gráfico 5. Valorización de los pits Whittle. Caso Andina Sur-Sur. ______________________________ 57
XI
1. Introducción. Uno de los grandes problemas al que se ven enfrentados los grandes yacimientos mineros explotados en la actualidad, es el no conocer en detalle toda la información referida al modelo de bloques. Esta situación hace que se utilice un determinado diseño de explotación que no permite una recuperación óptima de los beneficios. De esta manera, si al comienzo de la vida de un proyecto se estudia y se toma la decisión de llevar a cabo una transición del método de explotación a rajo abierto hacia una mina subterránea, entonces se enfrentará un nuevo escenario, en el cual el VAN y la vida útil del proyecto serán distintas y probablemente generen mejores expectativas. El método subterráneo indicado para soportar un proyecto de transición corresponde a los de caving, ya que éstos son los que devengan los menores costos de operación y alcanzan los mayores ritmos de producción. Esta visión permite, además, dar una mayor flexibilidad al diseño subterráneo. Sin embargo, ¿cómo se define el pit final en el cual se detiene la minería a rajo abierto y se da paso a un método subterráneo? ¿Está condicionado el diseño del pit final al proponer una opción subterránea en profundidad? En la actualidad no se tiene plena claridad acerca de cómo abordar este problema. Las técnicas de diseño y planificación utilizadas obtienen el máximo de cada una de las operaciones involucradas, pero es sabido que no se alcanza necesariamente el máximo global de un sistema si se maximizan cada una de sus componentes. Una metodología que permite abordar el problema de la transición corresponde a aquella que está incluida dentro del software de optimización de rajo Whittle de Gemcom. Ésta consta de cargar los beneficios de extraer un bloque vía subterránea como un costo para él. De esta manera, al momento de generar la envolvente óptima, se extraerán sólo aquellos bloques que generen un valor económico positivo, es decir, aquellos en que la diferencia entre el beneficio de extraer el bloque vía rajo y vía subterránea sea positiva (C.L. Wharton, 1996). Sin embargo, ésta forma de abordar el problema no considera la temporalidad de los beneficios ni tampoco realiza los descuentos asociados a la explotación subterránea. Además, esta metodología se basa en la elección de un piso para la subterránea, sin considerarlo variable. Es por éste motivo, y dado el número de variables que se encuentran involucradas en la decisión, que un método basado en un modelo de optimización es muy útil para encontrar una solución al problema planteado, permitiendo de ésta manera tomar la decisión simultánea de cómo debe ser extraído un determinado estrato del yacimiento, si por vía rajo o por vía subterránea. Ahora bien, para probar el supuesto que un yacimiento puede sostener una transición de rajo a subterránea, debe cumplirse que al menos algún estrato cumpla con el supuesto de que su valor económico descontado sea igual para la opción rajo y la opción subterránea.
1
Matemáticamente la condición de marginalidad se puede expresar de la siguiente forma: Condición marginal de transición
v
o t k
o k
v
u t k
u k
k 1 z
1
u t z
Donde; v
o k
RF l
o k
o 1 e C m
v
u k
RF l
u k
u C m
v
o, u k
Cp
Cp
Valor del estrato k explotado vía o/u.
o, u t k
1 (1
o, u r ) t k
Factor actualización por tiempo de minado o/u para el estrato k.
r
Tasa de descuento (10%).
R
Recuperación Metalúrgica. Factor conversión lb/ton. (Precio – Costo(fundición&refinacion))*R*α = Revenue Factor.
RF l
o, u k
e
Razón estéril/mineral.
o, u C m Cp
Ley de corte del estrato k para o/u.
Costo mina o/u. Costo planta.
2
Dada la condición de marginalidad anteriormente presentada, es posible formular un “problema de programación lineal”, el cual es posible resolver usando la herramienta Solver de
Excel. De esta manera, la condición de marginalidad de la transición se cumple si se dan los siguientes valores para las variables: Precio [US$/lb] 1.95 e 3.01
Recuperación [%] 80 Lk Rajo [%] 0.31
RF [US$/t/%] 28.97 Lk Subte [%] 0.61
Cm Rajo [US$/t]
Cm Subte [US$/t]
CP [US$/t]
1.02
5
5
Ritmo de producción Subte [ktpd]
Ritmo de producción Rajo [ktpd]
150
35
C(fundición y refinación) [US$/lb] 0.31
Tabla 1. Valores de las variables para cumplir el criterio de marginalidad de la transición.
Como los valores encontrados para que la condición de marginalidad de la transición se cumpla son comunes en las operaciones mineras, entonces es alta la posibilidad de que exista un estrato que satisfaga la igualdad en la realidad. Este ejercicio demuestra la existencia de un estrato marginal que define la transición desde el método a rajo abierto, por lo tanto, el estudio de encontrarlo en un caso real es totalmente válido y con una solución existente.
3
1.1. Objetivos. Los objetivos que persigue la presente tesis son los siguientes: 1. Probar a una escala de proyecto minero y, además, construir un modelo de programación matemática, para soportar la decisión de definir la transición en el método de explotación de una mina a rajo abierto a una de block caving subterránea. 2. Utilizar dos metodologías tradicionales para abordar el problema del diseño de la transición del método de explotación, y además el modelo de programación entera. 3. Cuantificar el efecto que se produce en la planificación de una transición en el método de explotación al utilizar una optimización secuencial y una optimización simultanea de los métodos de explotación.
4
1.2. Alcances. Se trabajará con dos modelos de bloques, el de la mina Palabora y el de Andina Sur-Sur. La mina Palabora está localizada en Sudáfrica, a 360 km al noreste de Pretoria, cerca del parque nacional Kruger. Palabora es el productor líder de cobre en Sudáfrica, y también la mayor fuente de vermiculita y oxido de zirconio. Los principales accionistas de Palabora Mining Co. son Rio Tinto plc (57.7%) y Anglo-American. Palabora comenzó sus operaciones a rajo abierto en 1964 y terminó en el 2002 cuando el pit alcanzó su máxima profundidad económica. El desarrollo de una mina subterránea para extraer el mineral remanente bajo el fondo del pit comenzó durante los últimos 5 años de la producción del rajo, a un costo de US$465m, con la esperanza de extender la operación de la mina en unos 20 años 1. La mina Andina, correspondiente a la división Andina de CODELCO Chile, opera el yacimiento Río Blanco, cuya riqueza era conocida desde 1920. Pero los intentos por iniciar su explotación no se concretaron hasta medio siglo después, en 1970. Está ubicada a ochenta kilómetros al noreste de Santiago, Chile, entre 3.700 y 4.200 metros sobre el nivel del mar. En la actualidad esta división realiza la explotación de minerales en la mina subterránea de Río Blanco y en la mina a rajo abierto Sur Sur. Andina produce unas 219.554 toneladas métricas anuales de concentrados de cobre que son materia prima fundamental para obtener el metal refinado. Además coloca en los mercados 2.133 toneladas métricas de molibdeno al año 2. En el desarrollo de las metodologías tradicionales se utilizará el software Whittle para definir la envolvente óptima para la mina a cielo abierto, y por otro lado, para definir la mina subterránea de block caving, se utilizará el módulo PCBC del programa GEMS. Para el diseño de las minas involucradas en el proceso de transición, se abordará como si se tratara de un proyecto de largo plazo. Para la estimación del VAN, se cubicarán las reservas extraíbles y se determinará su ley media, con lo cual se obtendrán los ingresos asociados a una determinada geometría del complejo minero. Una vez obtenidos los resultados asociados a las metodologías tradicionales presentadas, se compararán tomando en cuenta la geometría asociada al complejo minero y el VAN entregado. Para que los resultados sean comparables entre sí, se utilizarán parámetros de diseño y planificación comunes para todos los métodos probados, como por ejemplo un ángulo de talud de 45° y una altura de columna máxima de aproximadamente 300 metros. El modelo de optimización funcionará sobre pits y cilindros pre definidos que discretizarán el rajo y la subterránea respectivamente. Se asume una capacidad de producción constante para el rajo y la subterránea e igual al de las metodologías tradicionales. Se construirá
1 2
http://www.mining-technology.com/projects/palabora/ http://www.codelco.com/la_corporacion/fr_division_andina.html
5
una base de datos conveniente para ambos proyectos, de manera de realizar el estudio con el modelo matemático y con el diseño secuencial independiente. Se utilizará una estructura de costos e inversión conocidos además de una tasa de descuento fija e igual para el rajo y la subterránea, para cualquiera de las metodologías aplicadas. Finalmente, cabe hacer presente que el plan de producción del rajo y la mina subterránea se encuentra fuera del alcance del objeto de estudio del presente trabajo.
1.3. Metodología de trabajo. 1.3.1. Método secuencial.
6
El método secuencial consta principalmente de ejecutar un diseño y una planificación de largo plazo para el rajo y, posteriormente, con las reservas remanentes, definir el piso óptimo para la mina subterránea. Los pasos a seguir para el desarrollo del diseño y la planificación a largo plazo de la metodología son los siguientes: 1. Cargar el modelo de bloques en el software Whittle. 2. Definir el ángulo de talud para el rajo. 3. Ingresar los parámetros económicos necesarios para ejecutar el algoritmo de Lerchs & Grossman (L&G). 4. Valorizar los pits shells obtenidos del L&G con los parámetros económicos definidos para el rajo. 5. Generar un gráfico pit vs VAN y determinar el pit final. 6. Definir el valor del proyecto a rajo abierto. 7. Modificar el modelo de bloques cargado en GEMS, de manera que sólo queden los bloques remanentes. Es decir, modificar las leyes y la densidad de los bloques ya extraídos por el rajo. 8. Definir el valor económico para cada uno de los bloques, considerando los parámetros económicos definidos para la subterránea. 9. Sobre el modelo de bloques resultante, ejecutar corridas de Footprint Finder de manera de determinar la cota a la cual se debe colocar éste. 10. Calcular el VAN asociado a la mina subterránea. 11. Considerando el desfase en los ingresos de la mina subterránea, definidos por la vida útil del rajo, se calcula el VAN del proyecto.
1.3.2. Método incremental por fases.
El método incremental por fases consta en utilizar los volúmenes incrementales de roca para una serie de pits anidados definidos al usar el algoritmo de L&G, los cuales posteriormente serán evaluados económicamente considerando su extracción vía rajo o subterránea. De ésta manera, se construyen curvas de beneficios marginales de los incrementos para el rajo y para la 7
subterránea. En el caso de la mina subterránea, se tendrán tantas curvas paralelas como niveles del piso sean probados. Luego, el pit final será aquel en el cual se intercepten las curvas de beneficios marginales de los métodos de explotación. A continuación se muestra un gráfico en el que se representan las curvas antes señaladas.
Elevación del piso
VAN Incremental Subte
Gráfico 1. Curvas de VAN incrementales.
Los pasos a seguir para el desarrollo del diseño y la planificación a largo plazo de la metodología son los siguientes: 1. Cargar el modelo de bloques en el software Whittle. 2. Definir el ángulo de talud para el rajo. 3. Ingresar los parámetros económicos necesarios para ejecutar el algoritmo de L&G. 4. Calcular el tonelaje incremental entre los pits para estéril y para mineral. 5. Valorar los incrementos según los parámetros económicos definidos para el rajo y para la subterránea. 6. Graficar los beneficios marginales tanto de la mina rajo como la subterránea en un mismo gráfico. 7. Determinar el pit final asociado a una elevación del footprint y encontrar el VAN para la opción rajo con el pit final elegido y de la opción subterránea con la elevación del footprint determinada. 8. Determinar el VAN del complejo minero considerando el desfase que poseen los flujos de la mina subterránea. 9. Elegir la configuración que maximice el VAN del complejo. 8
1.3.3. Método de transición ILP.
Este método contempla utilizar un modelo de optimización entero. El modelo cuenta con una serie de variables binarias, las cuales definen la forma de minar cada estrato entregado. Esta decisión la realiza tomando en cuenta el valor económico del estrato y además es simultánea para un mismo estrato. Finalmente, la función objetivo de éste modelo es maximizar el VAN asociado, considerando todas las inversiones y costos involucrados. Los pasos a seguir para utilizar esta metodología son los siguientes: 1.
En base a los pits anidados por el software Whittle, se encuentran conos que representan el pit más grande, el más chico y un caso promedio.
2.
Para cada uno de los conos encontrados en el paso anterior, se replica en todo el modelo de bloques y se crean pits con la misma geometría.
3.
Los pits antes construidos son discretizados en estratos que siguen la forma del pit y luego son cubicados incrementalmente en GEMS.
4.
Se determinan los inputs para el modelo: a) Determinar el tiempo que se requiere para minar un estrato por la opción rajo. b) Determinar el tiempo que se requiere para minar un estrato por la opción subterránea. c) Valor económico del estrato al extraerlo por la mina rajo (P1). d) Valor económico del estrato al extraerlo por la mina subterránea (P2). e) Costo de construir la mina subterránea hasta el nivel donde se ubicará el footprint (cost_su).
5.
Cargar el modelo y los datos en CPLEX 10 y ejecutar el modelo.
2. Revisión bibliográfica. 2.1. Transición de minería a rajo abierto a subterránea. En la actualidad, existe una gran cantidad de minas con explotación a rajo abierto en las que los cuerpos mineralizados no se encuentran acotados en profundidad. Por esta razón se ha comenzado a tomar conciencia acerca de que al término de la vida útil de una mina a rajo 9
abierto es factible que continúe en producción, a través de la implementación de una minería subterránea bajo el rajo. Esta idea de plantear una opción de minería subterránea bajo un rajo, se presenta como una alternativa válida cuando no se vislumbra un pit final en un mediano plazo, es decir 5 a 15 años. Pero no se queda tan sólo con abrir la posibilidad de una minería subterránea, sino también se sugiere un diseño de block caving como “sucesor natural” dado sus “altos ritmos de producción, altos niveles de mecanización y sus bajos costos” (Fuentes, 2004). Sin embargo, el
hecho de no concebir al mismo tiempo la mina a cielo abierto y la opción de término subterránea, hace que se posterguen los estudios relacionados con ésta. Ello posteriormente llevará a realizar el estudio en un período más estrecho de tiempo, y “usualmente ‘algunos años’ no son suficientes para desarrollar seguro y correctamente un análisis en esta materia, principalmente cuando es necesario recolectar y definir la información base para los propósitos de una mina subterránea”. Por último , existen otros problemas no técnicos que afectan el proceso de transición; el cambio en el método de explotación requerirá de una capacitación de la mano de obra, un cambio de los proveedores de los equipos mineros y en el tipo de restricciones a las cuales se encuentra sujeta la mina (Fuentes, 2004). Además, existe un impacto en la definición del pit final de una mina a rajo abierto al barajar dentro de las opciones su futura transición a subterránea. El principal tema son los riesgos asociados a ésta decisión, que consisten en cambiar o adaptar la infraestructura, habilidades y conocimiento, y la habilidad de la organización para aceptar el cambio (Fuentes y Cáceres, 2004). Una forma en la cual se puede determinar el límite de la explotación a rajo abierto para luego pasar a una subterránea corresponde al análisis del costo diferencial. En dicho análisis se toma en cuenta el costo de oportunidad de la multiplicidad de los métodos para extraer un bloque, un estudio fase a fase, considerando que el volumen mineralizado encerrado dentro de las últimas fases del rajo fuera explotado por un método subterráneo, y, finalmente, un análisis iterativo de planta a planta, con el cual, se pretende encontrar la planta que representa la mejor transición, la cual reporta el mejor mineral para extraerlo por rajo y el mejor mineral remanente para ser extraído por el método subterráneo. Dada la simplicidad y los resultados entregados, la metodología de análisis fase a fase es aquella que mejor responde para realizar un análisis de este tipo (Avendaño, 2005).
10
Un caso real de diseño de una mina subterránea bajo un rajo corresponde al proyecto Chuqui Subterráneo. El méto do que se utilizará en este proyecto es el panel caving ya que, “es altamente mecanizable, con altos ritmos de producción y bajos costos”. Dentro de los retos técnicos que se plantean para el diseño de la mina subterránea se destacan los siguientes: “la profundidad de las excavaciones”, “la localización de los accesos”, “el número de frentes
mineras y la secuencia de los sectores de producción para lograr una producción de 125.000 tpd” (Arancibia y Flores, 2004). Dentro de las razones que justifican el realizar una transición de rajo a subterránea se encuentran la disponibilidad de una gran cantidad de recursos geológicos bajo el pit final, la infraestructura disponible en superficie, la trayectoria que tiene CODELCO manejando grandes yacimientos a través de métodos de block/panel caving, la disponibilidad de recursos financieros y el mejoramiento del valor presente neto de la división Codelco Norte y de Codelco Chile (Arancibia y Flores, 2004). Otros temas importantes al momento de analizar la transición de las operaciones a rajo abierto a una subterránea, corresponden a la geología, hidrogeología, la estimación de reservas, los estudios geotécnicos, el acceso a la mina y el diseño de ésta. Con respecto al cálculo de reservas, es aceptado en la mayoría de las minas de block/panel caving el uso del modelo PCBC del software GEMS. Este modelo simula el proceso subterráneo y entrega como resultado una ley y un tonelaje asociado a la opción, incluyendo en el análisis estimaciones de toppling de mineral desde las paredes del rajo a la zona de caving (Brannon, Casten, Johnson, 2004). Un tema que guarda relación con la transición rajo subterránea corresponde a los estudios geotécnicos relacionados con la estabilidad de las paredes del pit final. Por ejemplo, en los estudios realizados sobre las paredes del pit final de la mina Palabora, hacia finales del 2003 y comienzos del 2004, una vez que ya habían comenzado las operaciones subterráneas, y por ende el caving de la roca bajo el fondo del pit. Se encontraron movimientos anómalos en las estructuras presentes en los taludes, principalmente de la pared norte del rajo, por lo que era esperable que ésta fallara. Dado lo anterior, se inició el estudio y el modelamiento numérico de la falla con el software 3DEC. El objetivo del estudio acabado de los movimientos anómalos en las paredes del rajo era determinar si existía una relación directa entre la inestabilidad en las paredes del rajo y el desarrollo de la mina subterránea. Finalmente, como conclusión del modelo numérico cargado en el software 3DEC, se determinó que la inestabilidad de la pared norte es causada por las fracturas presentes en ella, las cuales forman cuñas y provocan la inestabilidad, por sobre el hecho de la aparición del caving en superficie. Además se concluyó a partir del modelo, que tomando en consideración las características del macizo rocoso en profundidad (presencia de agua, fallas, fracturas, etc), es factible que la pared norte del rajo se desplace si el mineral entre el fondo del rajo y el nivel de producción de la subterránea es removido (Brummer, Li, Moss, 2006).
11
Dentro de los principales puntos que se deben tener en cuenta al desarrollar un proyecto de transición, se destacan tres desafíos: incluir los aspectos geotécnicos relevantes para el diseño y la planificación de la mina, definir la infraestructura apropiada en la superficie como en profundidad, y definir los plazos del proyecto y las metas de producción. Además, otros temas importantes a considerar dentro de la transición son: la optimización del rajo, el plan de cierre del rajo, la posible interacción entre el rajo y la subterránea, y el diseño de la mina subterránea. (Olavarría, Adriasola y Karzulovic, 2006). Respecto a los aspectos geotécnicos relevantes en un estudio de esta materia, se considera que la presencia de un rajo en superficie produce zonas de concentración de esfuerzo y de bajo confinamiento. Este esfuerzo inducido afectará la propagación del caving y debe ser considerado en la forma en que se espera que éste se propague a través de la columna de mineral que será extraída. Además, en el caso de Chuquicamata, el estéril es una roca blanca, débil y muy fracturada, la cual podría ser fuente de una temprana dilución si la cavidad alcanza la falla oeste antes de conectarse con el fondo del rajo. Un tema no menor corresponde a aquellos esfuerzos inducidos, los cuales se espera que sean relativamente altos y generen sismicidad durante el desarrollo de la mina subterránea, la cual eventualmente podría generar estallidos de roca (Olavarría, Adriasola y Karzulovic, 2006). Es importante considerar al comienzo de una mina subterránea de caving, el lugar donde se emplazará el nivel de hundimiento, esto a raíz de las características constitutivas de la roca, la cual puede afectar positiva o negativamente el comienzo del quiebre de la misma. Además, se debe prestar atención a la interacción que se genere al momento que el caving alcance el fondo de rajo, en cuanto a la estabilidad de las paredes como en el comportamiento del caving. La idea de mantener una operación simultánea del rajo y la subterránea está sujeta a un diseño robusto y estable del crownpillar que las separa (Olavarría, Adriasola y Karzulovic, 2006). Como se puede observar, la mayoría de los trabajos relacionados con la transición de los métodos de explotación tienden a concentrarse en la estabilidad de las infraestructuras ya construidas así como también del pit final alcanzado. Sin embargo, no se realiza un estudio temprano del uso de una opción subterránea como una continuación de la vida útil del yacimiento. Es más, al momento de definir o probar nuevos recursos en profundidad se presenta el problema de la existencia de una infraestructura asociada a la operación del rajo dentro de los límites de una posible expansión. De hecho, un estudio acabado de transición podría llegar a los 20 años en algunos casos de minas muy grandes, para poder considerar todos los aspectos importantes tales como la definición del yacimiento, las características del macizo rocoso, el estudio de los métodos subterráneos aplicables, entre otros (Stacey, Terbrugge, 2000). De esta manera se observa que la decisión de la transición no debe responder a una alternativa de último recurso, ya que existen una serie de consideraciones que requieren de tiempo de estudio al momento de realizar la transición.
12
2.2. Planificación en minas a rajo abierto. La planificación minera se puede dividir en tres clasificaciones: la planificación estratégica, la planificación conceptual y la planificación operativa. La planificación estratégica tiene como objetivo el conectar el estado del mercado de minerales con la misión y objetivos del dueño de la compañía minera (Rubio, 2008). Por su parte, la planificación conceptual define los recursos existentes para alcanzar las promesas productivas definidas como parte de la planificación estratégica del proyecto. Finalmente, la planificación conceptual involucra el desarrollo de una envolvente económica, el diseño del proyecto, la construcción del programa de producción, la clasificación entre recursos y reservas del proyecto, la valorización y la estructura de costos asociada, y la determinación de los indicadores de desarrollo sustentable. La metodología a través de la que habitualmente se diseñan los pits finales de las minas de operación a rajo abierto corresponde a la aplicación algoritmos que permitan encontrar un conjunto de bloques que, sujetos a algunas restricciones como por ejemplo el ángulo de talud, obtengan el máximo beneficio. Entre estos algoritmos se encuentra el de Lerchs y Grossmann, 1965. Este método permitirá encontrar la envolvente óptima para un pit. El algoritmo se encuentra programado en la mayoría de los actuales softwares y se afirma que este procedimiento es favorable y, a la vez, eficiente computacionalmente hablando. Dentro del problema de la obtención del pit final o envolvente económica, tenemos que el pit final está definido por el ángulo de talud, la ley crítica de diseño y un modelo valorado de bloques (Magñin, 2002). Esta serie de datos son interpretados como un input que, al ser tratados bajo diferentes algoritmos computacionales de optimización, entregan un cono en el cual se encuentran contenidos los bloques potencialmente extraíbles para el proyecto. Posteriormente, se valorizan cada una de las geometrías obtenidas anteriormente y se procede a la definición del pit final, decisión que está conectada directamente con el VAN del proyecto. 2.2.1.
Metodología tradicional: definición del pit final
La metodología para determinar los límites económicos de un rajo se basa principalmente en definir las reservas mineras a partir de los recursos geológicos. Para ello se llevan a cabo dos etapas de planificación: definir el diseño del pit final y luego construir el plan de producción. Para realizar el diseño del pit final se necesita tener un modelo de bloques, donde cada bloque posea la información relevante que represente al macizo rocoso, por ejemplo, la densidad y cantidad de elementos presentes en el, tales como la especie de interés, contaminantes, entre otros. Además, se deben especificar las restricciones que tendrá el diseño, dentro de los que destacan el ángulo de talud, el fondo de rajo operativo, etc. Dada la información antes descrita, se ejecuta algún algoritmo para generar múltiples geometrías, las que posteriormente serán valorizadas y permitirán definir el pit final.
13
El algoritmo más conocido y utilizado por las herramientas actuales es el algoritmo de grafos de Lerchs - Grossman (1965). Su función objetivo busca maximizar el beneficio total de un pit, basado en el beneficio neto de cada bloque y su ubicación física en la mina. El beneficio neto de un bloque es la diferencia entre el valor total de extraer el bloque y el costo de extracción del material de la mina y el procesamiento del bloque en la planta. La restricción de precedencia permite asegurar que un bloque no se explote si es que los bloques sobre él no han sido explotados anteriormente. El límite del pit final representa un límite estático de bloques que maximizan el valor no descontado del mineral en el cuerpo mineralizado. Hochbaum y Chen (2000) aseguran que la heurística también es usada para obtener los límites del pit final. Sin embargo, muchas soluciones son erróneas o no son capaces de limitar el valor de la función objetivo que derivan. Pana (1965) creó el cono flotante heurístico, también conocido como el cono móvil o el cono heurístico dinámico, donde para cada bloque con beneficio positivo se genera un cono invertido. El algoritmo busca los conos tal que el beneficio total sea positivo. Estos conos son agrupados para generar el pit final. El supuesto de que cada cono en el pit final es rentable es incorrecto, en efecto, un pit óptimo podría consistir en la colección de conos donde ninguno por sí solo tiene un valor positivo, pero en conjunto sí lo tienen. Además se mencionan otros algoritmos heurísticos como el de Koborov (1974) que corresponde a una modificación del cono móvil. Hochbaum y Chen realizan un estudio comparativo del algoritmo de Lerchs-Grossman versus el algoritmo de Push-Relabel. Según se indica en el trabajo de estos autores, el algoritmo de Push-Relabel fue propuesto por primera vez de manera genérica por A.Goldberg (1987). Luego, Goldberg y Tarjan (1988) describen una implementación del algoritmo usando árboles dinámicos, obteniendo una muy buena eficiencia con respecto al tiempo de ejecución. 2.3. Planificación de minas subterráneas de block caving. Para la mayoría de los depósitos que cumplen ciertas condiciones de tamaño y ley, principalmente los llamados pórfidos, es casi indudable que un método de caving puede ser aplicado. Es esencial además que el macizo rocoso involucrado sea de una baja competencia, de manera que se posibilite la propagación de caving y de ésta manera la eficiencia del método (Tobie, Julin, 1982). Dentro de los métodos que interesan para el desarrollo de una transición corresponde a los de hundimiento, específicamente el block c aving, “método de explotación masivo en el cual un bloque de mineral en algunos casos representando el área basal del cuerpo mineralizado se corta en su base y luego a partir de la extracción se produce la propagación del hundimiento” (Rubio, 2007), y el panel caving, el cual “es una forma del método de hundimiento en que
bloques consecutivos se hunden en forma continua de modo de evitar la dilución lateral y los esfuerzos de relajación producidos en el método convencional de block caving” (Rubio, 2007).
14
En la actualidad, la forma en la cual se diseña y planifica el desarrollo de una mina subterránea es a través de la utilización del módulo PCBC del software GEMS. En este módulo se encuentra incluida la herramienta Footprint Finder. Esta herramienta realiza una optimización sobre el modelo de bloques y, además, considerá una mezcla vertical del tipo volumétrica de Laubsher (Diering, 2000). Para determinar el valor económico de cada una de las columnas, utiliza una tasa de descuento vertical, la cual es aplicada dado un parámetro de extracción vertical denominado Vertical Mining Rate (VMR). Finalmente, se compara el valor de cada una de las columnas con el costo de desarrollar un nivel y alcanzar esta columna, de manera de determinar si la extracción de dicha columna genera ingresos para, posteriormente, buscar la altura en la cual se produce el máximo valor de la columna. (Pais, 2009). 2.4. Programación matemática asociada a la transición rajo subterránea. En la actualidad, dada la alta contingencia del problema de la transición rajo subterránea se han planteado varias investigaciones orientadas a de definir el límite final óptimo del pit. La mayoría de estas investigaciones se basan en una optimización sobre el valor económico global de la explotación. Los modelos que se han tratado son algoritmos de programación lineal y un modelo heurístico. En el caso del modelo heurístico, se utiliza el valor económico de cada uno de los bloques al ser extraído por el método rajo y subterráneo, y luego se busca la combinación de estratos al ser extraídos mediante rajo y subterránea, de manera de maximizar el valor del proyecto (Bakhtavar, Shahriar y Oraee, 2008). El procedimiento anterior se realiza dos veces para cada corrida. La primera, con un tamaño de bloques mayor (50 m x 50 m), busca encontrar de manera superficial la ubicación del estrato de transición. Posteriormente, el estrato encontrado es refinado a un tamaño de bloque menor (12.5 m x 12.5 m) de manera de refinar el fondo del rajo. Cabe destacar que éste método se encuentra aplicado en forma bidimensional. Por otro lado, se encuentran los modelos de programación entera. En estos, se tiene una función objetivo (maximización del VAN) sujeta a ciertas restricciones, las que buscan de alguna manera representar la operación minera (precedencias, capacidad planta, etc.). En relación al modelo de programación lineal, se tiene el cuerpo mineralizado discretizado en estratos, los cuales pueden ser extraídos vía rajo o subterránea. Se espera que como resultado del proceso de optimización se determinen los estratos que deben ser minados por cada método. Para realizar esta elección, el modelo compara el ingreso económico actualizado para cada estrato. Dado lo anterior, se busca maximizar el valor del proyecto considerando cualquiera de los estratos menos el costo de desarrollo de la infraestructura de la mina subterránea, sujeto a restricciones mineras y de secuencia. (Newman, Rubio, Yano, 2008). Dentro de las restricciones más representativas que se pueden dar en el modelo subterráneo se tienen la conexidad de los bloques y la subsidencia. Por otro lado, el rajo está restringido por la definición del ángulo de talud global, además de las precedencias para extraer los bloques. Finalmente, una restricción que compete a ambos diseños corresponde a la capacidad de tratamiento, el cual viene dado por cada uno de los métodos utilizados (Azcarategui, 2005). 15
La forma en la cual se pueden cargar los datos de los modelos de bloques a los métodos de optimización son múltiples, pasando desde considerar bloques individuales u agrupaciones de ellos (macro bloques) que da cuenta del estéril y el mineral (Cáceres, 2006), o bien, considerando volúmenes con las formas ad hoc a los métodos de explotación determinados, como por ejemplo conos para discretizar el rajo y cilindros para la subterránea (Newman, Rubio, Yano, 2008). Aquellos trabajos que utilizan el modelo de bloques íntegramente hace que el resolver el problema sea inmanejable dada la cantidad de variables. Es por esto que las pruebas que se han hecho en éstos trabajos sean sólo sobre un perfil (Azcarategui, 2005). Por otro lado, aquellos que generan un rebloqueo generando bloques de un tamaño mayor, pierden la sensibilidad de la ubicación exacta del nivel de transición. Considerando el trabajo que compara el resultado obtenido por un modelo de optimización que usa una agregación en macrobloques, se podría llegar a obtener hasta un 50% de mejora en el VAN de un proyecto desarrollado de la manera tradicional, dado principalmente porque al optimizar conjuntamente la opción rajo y la subterránea, permite obtener mayores beneficios asociados a las mejores leyes extraídas por la mina subterránea (Cáceres, 2006). La mayoría de los trabajos presentados anteriormente, salvo el de Cáceres (Cáceres, 2006) y Newman, Rubio y Yano (Newman, Rubio, Yano, 2008), han sido desarrollados en base a modelos de bloques ficticios y aplicándolo sólo a secciones 2D. Si bien han alcanzado un resultado, que involucra definir un nivel de transición, su aplicabilidad en un proyecto real es poco viable y además, al ser un trabajo de investigación, no posee un desarrollo convencional para abordar la transición con el cual ser comparado.
16
3. Modelo de optimización para la transición rajo – subterránea. 3.1. Supuestos. En la formulación del modelo de programación lineal entera (ILP) existen algunos supuestos. El primero de ellos corresponde a la forma en la cual se introducen los datos. El modelo de optimización utiliza una pre agregación del modelo de bloque en unidades de reserva, que corresponden a una representación aproximada de la forma en la cual será explotada la mina por los métodos mineros utilizados (rajo – block caving). Para el caso del rajo, éstas unidades adquieren una geometría de cono, y en el caso de la mina subterránea, son cilindros. De esta manera, el modelo de optimización toma decisiones sobre extraer o no cada una de éstas unidades de reserva que discretizan el modelo de bloques. Cada una de éstas unidades de reserva son cubicadas y valorizadas externamente al modelo. El modelo de optimización se encuentra definido sólo para realizar la transición desde el método a rajo abierto a una mina subterránea de block caving. Tanto como el ritmo de producción de los métodos de explotación, los costos de desarrollo y los demás parámetros económicos se asumen conocidos. 3.2. Generación de las unidades de reserva para el modelo de optimización. Utilizando los pits generados por la corrida de Whittle sobre un modelo de bloques, se eligen arbitrariamente las geometrías que correspondan al máximo pit, al mínimo pit y a un caso promedio. Esto se obtiene definiendo el fondo correspondiente y proyectando el talud a un ángulo de 45° hasta interceptar la topografía. De esta manera, se obtienen tres sets de pits con forma de conos que discretizarán el yacimiento en estratos extraíbles a través del rajo.
Figura 1. Unidades de reserva que discretizan la explotación a rajo abierto.
17
Para discretizar la extracción del mineral a través del block caving se utiliza una corrida de footprint finder, a partir de la cual se definen cilindros no regulares, que encierran las columnas favorables económicamente (A), o bien las columnas con ley media sobre una ley de corte especificada (0.2, 0.3, 0.4, 0.5, 0.6, 0.8 o 1.0) (B).
Figura 2. Vista en planta de cilindros que discretizan la explotación subterránea.
Existen dos alternativas para realizar la valorización de los estratos. La primera de ellas consiste en predefinir una ley de corte para el rajo (0.4 %Cu) y una para el block caving (0.6 %Cu). La segunda alternativa corresponde a realizar la valorización de cada uno de los estratos para un set preestablecido de leyes de corte (0.2, 0.3, 0.4, 0.5, 0.6, 0.8 o 1.0), es decir, generando una curva beneficio ley para cada estrato.
Para cada una de las valorizaciones se consideran parámetros económicos y ritmos de producción conocidos e iguales a los utilizados en las otras metodologías de manera de que los resultados sean comparables.
18
3.3. Descripción del modelo. 3.3.1. Función objetivo.
3.3.2. Sets.
Los sets en los cuales se guardan los índices relevantes para el modelo de programación lineal son los siguientes: 1. S: Corresponde al conjunto donde se guardan todos los estratos. 2. T: Corresponde al conjunto donde se tienen todos los períodos disponibles para realizar minería. 3. G: Corresponde al conjunto de leyes de corte para las cuales se desea evaluar el estrato. 4. TS{s, g}: Corresponde al conjunto donde se guardan los períodos de inicio del minado por rajo de cada estrato asociado a una ley de corte g.
3.3.3. Parámetros.
Los parámetros que son utilizados en el modelo de programación lineal son los siguientes: 1. Discount_rate: Tasa de descuento con que se evaluará el proyecto. 2. Rho1{s, g}: Número de períodos que se requieren para minar el estrato s, a una ley de
corte g, con la opción rajo. 3. Rho2{s, g}: Número de períodos que se requieren para minar el estrato s, a una ley de
corte g, con la opción subterránea.
19
4. P1{s, g}: Beneficio correspondiente a minar el estrato s, con una ley de corte g, por la
mina rajo. 5. P2{s, g}: Beneficio correspondiente a minar el estrato s, con una ley de corte g, por la
mina subterránea. 6. F{g}: Costo adicional por extraer una ley de corte en particular. 7. Shaft_cost{s, g}: Costo de construir la infraestructura necesaria para producir con la
mina subterránea. 8. H: Número de estratos que conforman la altura máxima de extracción para las
columnas.
3.3.4. Variables.
Las variables utilizadas en el modelo de programación lineal son las siguientes:
20
3.3.5. Restricciones.
Explotar sólo una vez cada estrato.
Secuencia explotación por rajo.
21
Secuencia explotación subterránea.
Secuencia relativa entre métodos de explotación. El rajo antes que el block caving.
No exista extracción simultánea.
No existan estratos sin explotar entre aquellos que son explotados.
Por sobre la base del crownpillar no debe existir extracción por el método subterráneo.
Bajo el crownpillar sólo es admisible la extracción por método subterráneo.
Sobre el comienzo del crownpillar sólo es admisible la extracción por método rajo.
El estrato por sobre la base del footprint debe ser minado por método subterráneo.
22
Los estratos bajo la base del footprint no pueden ser minados por el método subterráneo.
Sólo se puede definir el comienzo del crownpillar una vez.
Sólo se puede definir la base del crownpillar una vez.
Sólo se puede definir el comienzo del footprint una vez.
Sólo se puede definir la base del footprint una vez.
Si el crownpillar existe, entonces se debe definir el estrato de comienzo y de final.
Si el footprint existe, entonces se debe definir el estrato de comienzo y de final.
El crownpillar debe estar por sobre el footprint.
23
Sólo se puede usar una ley de corte para cada método de explotación.
Sólo es posible minar estratos que tengan la misma ley elegida para un determinado método de explotación.
En todos los períodos debe ser extraído un estrato, a menos que pertenezca a un gap.
Cuando se definen los gaps, en adelante, todos los periodos serán gaps.
El pique de acceso debe estar construído antes que se inicie la explotación subterránea.
El pique debe llegar hasta el estrato definido como la base del footprint.
La mina subterránea posee una restricción de estratos que pueden ser explotados (altura de columna de extracción).
24
4. Información disponible para el caso de estudio de Palabora. 4.1. Inventario de recursos. En el caso del modelo de bloques de Palabora se tiene lo siguiente; Tamaño del Bloque: X: 40 m, Y: 20 m, Z: 60 m. Numero de Bloques: X: 71, Y: 149, Z: 37. Densidad promedio: 2.65 [t/m 3] Además, se debe considerar que la topografía es plana y se encuentra ubicada en la cota 296 m. La curva tonelaje ley para el modelo de bloques de Palabora se encuentra expresada en la siguiente tabla.
Gráfico 2. Curva tonelaje - ley para Palabora.
4.2. Parámetros técnico-económicos para Palabora. Se considera que el ángulo de talud global para el pit final que se diseña corresponde a un ángulo constante para todas las paredes del rajo, con un valor de 45°. Se asume también una dilución minera del 5% y una recuperación minera del 95%. Además, no se aplicará una estructura de ajuste sobre el costo mina (CAF), así como tampoco el costo de proceso (PCAF).
25
Se utilizará una serie de revenue factors en el software Whittle entre 0.3 y 1.5 con un paso de 0.015, reproduciendo 81 valores distintos. La tasa de descuento utilizada para el cálculo del VAN será fija con un valor del 10%. Para realizar la evaluación de los estratos al ser minados por el método rajo, se consideran los siguientes parámetros económicos. Parámetros económicos para el rajo. Precio 2 US$/lb Costo de 5 US$/t Proceso Costo Mina 1 US$/t Recuperación 90 % Costo de Venta 0.26 US$/lb Tabla 2. Parámetros económicos para el rajo. Caso Palabora.
Además, el ritmo de producción de la mina subterránea alcanza las 23,500 tpd. Este valor se encuentra al definir que aproximadamente se tendrán 290 columnas en el footprint y dado un vertical mining rate de 14, donde este parámetro corresponde a los metros de roca que se extraerán aproximadamente en la vertical sobre cada una de las columnas, entonces se tiene que el ritmo de producción anual por columna alcanza a 29,680 tpa/col, que es equivalente a producir 23,500 tpd. Con respecto a los parámetros económicos utilizados en las corridas de footprint finder, se considerará para el desarrollo de la mina subterránea un costo de US$ 2000 por metro cuadrado y además para realizar la evaluación de los estratos se consideran los siguientes parámetros económicos. Parámetros económicos para el block caving Precio 2 US$/lb Costo de 5 US$/t Proceso Costo Mina 5 US$/t Recuperación 90 % Costo de Venta 0.26 US$/lb Tabla 3. Parámetros económicos para el block caving. Caso Palabora.
26
4.3. Datos de entrada para el modelo de optimización. Para poder correr el modelo de optimización deben definirse de antemano cada uno de los sets y parámetros involucrados. A continuación, se presentan los valores constantes para cada una de las formas en que se corrió el modelo.
Discount rate = 10% Por otro lado, para el caso de la ley de corte fija, se consideró como el set de leyes el siguiente.
Y para el caso de la ley de corte variable, se utilizó el siguiente vector de leyes de corte.
Un parámetro importante para el cálculo de la transición corresponde a la altura de columna máxima extraíble. Para este parámetro, se definió el siguiente vector, además de considerar que no existía restricción.
Finalmente, el modelo fue desarrollado en el lenguaje AMPL y se resolvió usando CPLEX Parallel 10.2 con 4 procesos. Los tiempos de solución fueron calculados en base a un computador con 4 procesadores (IntelXeon Cuadcore
[email protected] Ghz, 4 Gb RAM).
27
5. Análisis de la transición para el caso de estudio de Palabora. Lo que se espera obtener de cada una de las metodologías corresponde a una geometría asociada al complejo minero, esto es, una profundidad máxima para el rajo y una cota en la cual se debe ubicar el footprint de la mina subterránea. Además, se obtendrán los VAN de cada una de éstas, parámetro que funcionará como criterio de comparación. 5.1. Resultados para la metodología secuencial. De acuerdo a la metodología, se obtuvo un conjunto de pits anidados, cada uno de ellos valorizados a distintos ritmos de producción y a los mismos parámetros económicos ya definidos en la sección de parámetros técnico-económicos. De esta manera, se tiene el siguiente gráfico que resume los valores de los pits para cada ritmo de producción. Además se adjunta la tabla para el caso elegido de ritmo de producción de 180,000 tpd. Valorización de pits a diferentes ritmos de producción 4,000.00
3,000.00
2,000.00
$ 1,000.00 S U M N A V
-1,000.00
-2,000.00
-3,000.00 1 2 3 4 5 6 7 8 9 101112131415161718192021222324252627282930313233343536373839404142434445464748
Pits 60000 tpd
80000 tpd
100000 tpd
120000 tpd
140000 tpd
160000 tpd
Gráfico 3. Valorización de pits Whittle. Caso Palabora.
28
180000 tpd
190000 tpd
Ritmo de producción 180000 tpd Pit Final VAN [MUS$] Pit Final VAN [MUS$] Pit Final VAN [MUS$] Pit Final VAN [MUS$] Pit Final VAN [MUS$] 1 -2,140 16 2,379 31 3,483 46 3,593 61 3,558 2 -2,101 17 2,545 32 3,508 47 3,593 62 3,557 3 -2,064 18 2,666 33 3,524 48 3,593 63 3,556 4 -2,030 19 2,862 34 3,543 49 3,588 64 3,555 5 -1,957 20 2,908 35 3,546 50 3,587 65 3,553 6 -1,888 21 3,043 36 3,561 51 3,586 66 3,552 7 -1,815 22 3,074 37 3,570 52 3,584 67 3,550 8 -482 23 3,141 38 3,578 53 3,584 68 3,550 9 75 24 3,250 39 3,580 54 3,579 69 3,548 10 389 25 3,289 40 3,581 55 3,576 70 3,548 11 978 26 3,332 41 3,588 56 3,574 71 3,548 12 1,414 27 3,347 42 3,588 57 3,567 72 3,547 13 1,697 28 3,413 43 3,591 58 3,566 73 3,545 14 1,901 29 3,424 44 3,593 59 3,564 74 3,545 15 2,227 30 3,466 45 3,593 60 3,561 75 3,543
Tabla 4. Valorización de pits Whittle al ritmo de producción óptimo. Caso Palabora.
De esta manera, se elige el pit 45 como pit final. La información referida a este pit se encuentra en la siguiente tabla. Pit Final
Fondo
VAN [MUS$]
Mineral [Mt]
Estéril [Mt]
Vida útil [año]
REM
Ley Media [%Cu]
45
-1,147
3,593
2,146
8,604
33
4
0.483
Tabla 5. Características del pit final elegido para la metodología secuencial en Palabora.
Una vez definido el pit final, siguiendo con la metodología, se obtuvieron los resultados del footprint finder para determinar el nivel donde se debe ubicar la base del block caving. Nivel Tonelaje [kt] VAN [kUS$] Área [m2] Max H [m] Ley media %Cu 176 116 56 -4 -64 -124 -184 -244 -304 -364 -424 -484 -544 -604 -664 -724 -784 -844 127 148 800 180 0.91 -904 254 1,737 800 300 0.98 -964 3,180 12,441 8,800 300 0.82 -1,024 10,303 34,169 27,200 300 0.79 -1,084 14,831 57,245 28,800 300 0.80 -1,144 21,230 84,038 42,400 300 0.81 -1,204 28,327 98,179 43,200 300 0.79 -1,264 117,799 352,812 213,600 300 0.76 -1,324 153,173 468,147 238,400 300 0.75 -1,384 177,685 545,551 245,600 300 0.75 -1,444 192,223 587,124 248,800 300 0.75 -1,504 196,675 602,703 248,000 300 0.76
Tabla 6. Corrida de footprint finder para Palabora.
29
Dado que en el nivel -1,504 se obtiene el máximo VAN asociado al método subterráneo, entonces se determina que el footprint debe estar ubicado en esta cota. Sin embargo, el VAN reportado en éste piso es de MUS$ 602, el cual, desplazado en 33 años que corresponde a la extracción del rajo, reporta un VAN de MUS$ 26.8. Adicionalmente, cabe destacar que el costo de la infraestructura mayor de las labores subterráneas no se encuentran consideradas. A modo de resumen se presenta a continuación una tabla donde se sintetiza el valor del proyecto completo. Pit Final N° 45 Fondo del pit VAN [MUS$] Mineral [Mt] Estéril [Mt] Vida útil [años] -1,147 3,593 2,146 8,604 33 Footprint Nivel del FP VAN [MUS$] Mineral [Mt] Ley de Cu % Vida útil [años] -1,504 26.8 197 0.76 23 Proyecto de transición VAN [MUS$] 3,620
Tabla 7. Resultados para la metodología secuencial en Palabora.
A continuación, se muestra un esquema del diseño del complejo minero obtenido a través de la metodología secuencial.
Figura 3. Esquema del diseño de la transición según la metodología secuencial para Palabora.
Finalmente, la profundidad del rajo es de 851 [m] y el footprint del block caving se encuentra en el nivel -1504 con una altura de columna de 300 [m]
30
5.2. Resultados para la metodología incremental por fases. Para esta metodología, se estandariza la cota de los estratos según se presenta en la siguiente tabla. Estrato Cota Estrato Cota 1 236 16 -664 2 176 17 -724 3 116 18 -784 4 56 19 -844 5 -4 20 -904 6 -64 21 -964 7 -124 22 -1,024 8 -184 23 -1,084 9 -244 24 -1,144 10 -304 25 -1,204 11 -364 26 -1,264 12 -424 27 -1,324 13 -484 28 -1,384 14 -544 29 -1,444 15 -604 30 -1,504 Tabla 8. Estandarización de cotas.
De esta manera, el conjunto de pits valorizados a un ritmo de producción de 180,000 tpd y a los parámetros económicos establecidos, conformarán la base para construir los estratos que son extraídos por el método de superficie y se ajustarán a los niveles antes presentados. El parámetro para consolidar los estratos corresponde al nivel del fondo del pit, los cuales son presentados a continuación. Pit Cota del fondo 1 2 3 4 5 6 171 7 54.13 8 -12.21 9 -67.5 10 -126.5 11 -181.7 12 -218.6 13 -248.1 14 -277.6 15 -310.7
Pit Cota del fondo Pit Cota del fondo Pit Cota del fondo 16 -369.7 31 -787.2 46 -1,147 17 -369.7 32 -820.4 47 -1,178.3 18 -432.4 33 -848 48 -1,178.3 19 -432.4 34 -848 49 -1,178.3 20 -487.7 35 -848 50 -1,238.1 21 -528.2 36 -911 51 -1,238.1 22 -546.7 37 -911 52 -1,238.1 23 -546.7 38 -911 53 -1,269.1 24 -606.9 39 -965 54 -1,301 25 -609 40 -1,025.2 55 -1,301 26 -636.3 41 -1,025.2 56 -1,327.1 27 -667.7 42 -1,086.7 57 -1,356.7 28 -667.7 43 -1,086.7 58 -1,387.4 29 -726.4 44 -1,086.7 59 -1,388.5 30 -726.4 45 -1,147 60 -1,412.1 Tabla 9. Cota del fondo de los pits Whittle para Palabora.
31
Pit Cota del fondo 61 -1,412.1 62 -1,412.1 63 -1,412.1 64 -1,446.5 65 -1,504 66 -1,504 67 -1,504 68 -1,504 69 -1,504 70 -1,504 71 -1,504 72 -1,504 73 -1,504 74 -1,504 75 -1,504
Como resultado de lo anterior, se tiene la siguiente tabla que resume el valor de los estratos extraídos por el método a rajo abierto. Cota 236 176 116 56 -4 -64 -124 -184 -244 -304 -364 -424 -484 -544 -604
VAN del Tiempo VAN del Tiempo estrato rajo alimentación a Cota estrato rajo alimentación a [MUS$] planta [años] [MUS$] planta [años] 0 0 -664 163 5 0 0 -724 53 3 0 0 -784 17 1 -1,815 7 -844 63 6 1,333 2 -904 32 3 557 1 -964 3 1 314 1 -1,024 7 4 589 1 -1,084 5 4 436 2 -1,144 0 3 812 4 -1,204 -5 5 318 2 -1,264 -4 3 318 3 -1,324 -10 5 46 1 -1,384 -10 4 233 4 -1,444 -12 7 109 2 -1,504 -12 13 Tabla 10. Valorización de los estratos extraidos vía rajo. Caso Palabora.
Por otro lado, el valor de cada uno de los estratos extraídos por el método subterráneo son los siguientes. Cota Mineral [kt] Ley (% Cu) Cota Mineral [kt] Ley (% Cu) 236 127 1.06 -664 609,917 0.67 176 5,851 0.74 -724 646,296 0.67 116 13,356 0.72 -784 657,872 0.68 56 32,945 0.79 -844 717,909 0.68 -4 108,755 0.71 -904 744,494 0.68 -64 175,535 0.69 -964 763,955 0.69 -124 223,108 0.68 -1,024 849,561 0.68 -184 276,785 0.68 -1,084 913,415 0.68 -244 327,283 0.67 -1,144 946,106 0.68 -304 393,172 0.66 -1,204 974,598 0.68 -364 444,180 0.66 -1,264 1,030,945 0.68 -424 463,132 0.66 -1,324 1,066,814 0.69 -484 511,850 0.66 -1,384 1,121,382 0.69 -544 521,643 0.67 -1,444 1,148,986 0.69 -604 583,333 0.67 -1,504 1,189,560 0.69 Tabla 11. Ley y mineral para cada estrato de P alabora extraíbles vía subterránea.
32
Dada la información antes obtenida, es posible valorizar cada uno de los estratos considerando el costo de desarrollo y además el costo de la construcción del pique de acceso.
Cota
VAN estrato BC [MUS$]
Tiempo de explotación [año]
Cota
VAN estrato BC [MUS$]
Tiempo de explotación [año]
236 176 116 56 -4 -64 -124 -184 -244 -304 -364 -424 -484 -544 -604
1 29 48 237 818 592 399 456 415 542 401 118 369 58 527
1 1 1 2 9 8 6 6 6 8 6 2 6 1 7
-664 -724 -784 -844 -904 -964 -1,024 -1,084 -1,144 -1,204 -1,264 -1,324 -1,384 -1,444 -1,504
189 289 87 579 232 168 922 615 265 228 524 305 503 230 359
3 4 1 7 3 2 10 7 4 3 7 4 6 3 5
Tabla 12. Valorización de los estratos extraídos vía subterránea. Caso Palabora.
Sin embargo, para comparar el valor de los estratos entre rajo y subterránea, éstos últimos deben encontrarse descontados al inicio del proyecto, considerando la previa extracción del rajo sobre ellos. De esta manera, se tendrá la siguiente tabla de comparación de los valores descontados para cada estrato. Cota 116 56 -4 -64 -124 -184 -244 -304 -364 -424 -484 -544 -604 -664 -724 -784 -844 -904 -964 -1,024 -1,084 -1,144 -1,204 -1,264 -1,324 -1,384 -1,444 -1,504
Flujo de caja Flujo de caja descontado para el descontado para el estrato rajo [kUS$] estrato BC [kUS$] 0 0 -1,815,100 0 1,332,868 1 557,412 2 313,632 3 589,314 7 436,299 12 812,489 36 317,724 49 317,761 17 45,595 99 232,516 16 109,478 256 163,076 139 52,962 286 16,784 80 63,318 1,161 31,606 679 2,656 555 7,150 7,396 5,233 10,144 312 8,141 -4,821 8,046 -4,329 38,598 -9,865 -9,795 -11,644 -11,840
Tabla 13. Flujos de caja descontados para los estratos extraídos por rajo y subterránea. Caso Palabora.
33
Los valores económicos para las dos minas antes descritas se encuentran en la tabla dispuesta a continuación, considerando que el valor de la subterránea se calculó descontando la extracción previa del rajo. Pit Final Fondo del pit [m] VAN [MUS$] Mineral [Mt] Estéril [Mt] Vida útil [años] -964 3,580 1,757 5,587 27 Footprint Nivel del FP [m] VAN [MUS$] Mineral [Mt] Ley de Cu % Vida útil [años] -1,264 72 764 0.67 32 Proyecto de Transición VAN [MUS$] 3,653
Tabla 14. Resultados para la metodología incremental por fases en Palabora.
A continuación se muestra un esquema del diseño del complejo minero obtenido a través de la metodología incremental por fases.
Figura 4. Esquema del diseño de la transición según la metodología incremental por fases para Palabora.
Por lo tanto, se concluye que el rajo debe ser explotado hasta la cota -964, alcanzando una profundidad máxima de 668 [m], para inmediatamente pasar a una explotación de block caving con una altura de columna de 300 m y su footprint ubicado en la cota -1,264.
34
5.3. Resultados para la metodología de transición ILP en Palabora. 5.3.1. Aspectos previos. i) Escenarios para definir la transición con el modelo ILP.
Para la prueba que se realizará con el modelo, se considerarán las siguientes combinaciones de valorizaciones, decisiones del modelo y configuración de parámetros.
ii) Valorización de los estratos.
Una vez construidos las unidades de reserva, éstas son cubicadas en GEMS y posteriormente evaluadas tomando en cuenta los siguientes parámetros. Parámetros de evaluación de los estratos tipo rajo. Ritmo de 180 ktpd Producción Precio 2 US$/lb Costo de Proceso
5
US$/t
Costo Mina
1
US$/t
Recuperación
90
%
Costo de Venta
0.26
US$/lb
Tabla 15. Parámetros de evaluación de los estratos para el ILP extraíbles vía rajo. Caso Palabora.
35
De esta manera, se calcula el valor de cada estrato considerando una tasa de descuento del 10%, además del tiempo que toma extraer el estrato mediante un método a rajo abierto con un ritmo de producción de 180 ktpd. A continuación, a modo de ejemplo, se presenta la tabla de valorización para una ley de corte de 0.6 %Cu. Ley de corte (%Cu) 0.60
Estrato
Mineral (kt)
Estéril (kt)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
29.7 5,468 16,169 24,092 27,560 29,441 29,340 29,488 28,875 29,211 30,139 29,852 31,514 30,926 30,592 32,072 30,864 30,331 31,602 30,957 30,152 31,798 32,535 32,531 32,726 33,267 33,523 33,489 33,122 33,330
21,352 25,067 32,755 46,787 68,563 96,122 129,407 166,175 206,773 250,788 298,065 349,370 403,276 463,121 526,491 590,775 662,313 737,266 812,897 895,049 981,263 1,067,636 1,159,706 1,256,253 1,356,220 1,458,608 1,565,954 1,677,357 1,791,465 1,910,020
Ley de Cu (%Cu) 0.64 0.78 0.79 0.80 0.79 0.80 0.79 0.80 0.81 0.82 0.81 0.80 0.79 0.80 0.80 0.80 0.81 0.80 0.80 0.80 0.81 0.81 0.82 0.82 0.81 0.81 0.81 0.81 0.81 0.81
Potencia del estrato (m) 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60 60
Rho1 (años) 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
P1 (MUS$) -21 89 311 471 521 536 499 470 426 397 361 294 266 202 137 101 14 -80 -129 -223 -316 -368 -438 -537 -636 -726 -833 -946 -1,066 -1,183
Tabla 16. Valorización de los estratos para el ILP extraíbles vía rajo. Caso Palabora.
Por otro lado, se tienen los estratos que discretizan la mina subterránea, los cuales son valorizados siguiendo los parámetros presentados a continuación. Adicionalmente, se debe considerar una inversión que corresponde al desarrollo del acceso y la infraestructura necesaria para la producción. Estos costos se encuentran en la siguiente tabla. Parámetros de evaluación de los estratos tipo subterráneo. Ritmo de 23.5 ktpd producción Precio
2
US$/lb
Costo de Proceso
5
US$/t
Costo Mina
5
US$/t
Recuperación
90
%
Costo de Venta
0.26
US$/lb
Costo de inversiones para la subterránea. Costo de desarrollo
520
k$/dpt
Tramming
5
kUS$/m
Costo del pique 8 kUS$/m Tabla 18. Costo de inversiones para la subterránea. Caso Palabora.
Tabla 17. Parámetros de evaluación de los estratos para el ILP extraíbles vía subterránea. Caso Palabora.
36
A continuación, y a modo de ejemplo, se presenta una tabla que contiene la cubicación y la valorización de los estratos subterráneos a una ley de corte de 0.6 % Cu. Cil_A Lc=0.6 Estrato
Mineral (Mt)
Ley de Cu (%)
Footprint (m2)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
33 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35
0.21 0.44 0.62 0.64 0.65 0.65 0.68 0.70 0.71 0.72 0.72 0.72 0.70 0.70 0.69 0.70 0.71 0.71 0.76 0.75 0.75 0.76 0.76 0.75 0.75 0.76 0.76 0.76 0.76 0.76
204,197 209,571 211,362 212,258 212,795 213,153 213,409 213,601 213,750 213,870 213,967 214,049 214,118 214,177 214,228 214,273 214,312 214,347 214,379 214,407 214,433 214,456 214,477 214,497 214,515 214,531 214,546 214,561 214,574 214,586
Costo de desarrollo (MUS$) 409 420 423 425 426 426 427 427 428 428 428 428 428 428 428 429 429 429 429 429 429 429 429 429 430 430 430 430 430 430
Costo de sostenimiento (MUS$) 7 8 10 11 12 13 15 16 17 18 19 21 22 23 24 25 27 28 29 30 31 33 34 35 36 37 39 40 41 42
Costo de tramming (MUS$) 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
Costo de transición (MUS$) 418 430 435 438 440 442 444 445 447 448 450 451 453 454 455 457 458 459 460 461 463 464 465 466 468 469 470 472 473 474
Rho 2 (año)
P2 (MUS$)
4 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
-89 174 400 424 429 433 466 490 509 517 516 523 499 494 481 492 507 510 560 555 557 561 566 556 548 569 568 569 570 568
Tabla 19. Valorización de los estratos para el ILP extraíbles vía subterránea. Caso Palabora.
Una vez valorizadas las agregaciones necesarias para correr el modelo, solo resta definir el set TS, el cual, define el período de inicio de la extracción del estrato, cubicado a una cierta ley de corte, por el método de rajo abierto. Este vector se obtiene al acumular el tiempo de explotación de cada uno de los estratos. A continuación se presenta un vector de ejemplo. set TS[1,0.6]:=1; set TS[11,0.6]:=11; set TS[21,0.6]:=21; set TS[2,0.6]:=2; set TS[12,0.6]:=12; set TS[22,0.6]:=22; set TS[3,0.6]:=3; set TS[13,0.6]:=13; set TS[23,0.6]:=23; set TS[4,0.6]:=4; set TS[14,0.6]:=14; set TS[24,0.6]:=24; set TS[5,0.6]:=5; set TS[15,0.6]:=15; set TS[25,0.6]:=25; set TS[6,0.6]:=6; set TS[16,0.6]:=16; set TS[26,0.6]:=26; set TS[7,0.6]:=7; set TS[17,0.6]:=17; set TS[27,0.6]:=27; set TS[8,0.6]:=8; set TS[18,0.6]:=18; set TS[28,0.6]:=28; set TS[9,0.6]:=9; set TS[19,0.6]:=19; set TS[29,0.6]:=29; set TS[10,0.6]:=10; set TS[20,0.6]:=20; set TS[30,0.6]:=30; Tabla 20. Vector de parámetros de inicio del minado del estrato vía rajo. Caso Palabora.
Finalmente todos los datos calculados anteriormente se consolidan en un archivo de texto (.dat) con el formato correspondiente para ser leído por Cplex. 37
iii) Resultados del modelo de optimización para Palabora.
En las siguientes tablas se encuentran los resultados obtenidos para las corridas realizadas con el modelo de optimización, utilizando las diferentes agregaciones definidas anteriormente. Para cada una de las corridas se obtuvo un VAN y una geometría asociada, las cuales se compararán más adelante con los resultados obtenidos en las metodologías tradicionales. (a) Resultados para ley de corte fija predefinida y dos discretizaciones de cilindros distintos (A y B).
Para la discretización máxima, promedio y mínima del pit final se obtuvo el siguiente resultado con respecto al VAN, considerando la ley de corte para el rajo de 0.4 %Cu y para la subterránea de 0.6 %Cu. MAX PIT MED PIT MIN PIT HOD = 120 Convergencia VAN HOD = 120 Convergencia VAN HOD = 120 Convergencia VAN [m] [s] MUS$ [m] [s] MUS$ [m] [s] MUS$ A 11 1,997 A 8 2,613 A 8 2,519 B 12 1,952 B 9 2,587 B 18 2,519 HOD = 180 Convergencia VAN HOD = 180 Convergencia VAN HOD = 180 Convergencia VAN [m] [s] MUS$ [m] [s] MUS$ [m] [s] MUS$ A 8 2,067 A 12 2,652 A 12 2,529 B 16 1,993 B 12 2,610 B 18 2,515 HOD = 240 Convergencia VAN HOD = 240 Convergencia VAN HOD = 240 Convergencia VAN [m] [s] MUS$ [m] [s] MUS$ [m] [s] MUS$ A 14 2,108 A 8 2,674 A 9 2,543 B 23 2,022 B 17 2,628 B 21 2,524 HOD = 300 Convergencia VAN HOD = 300 Convergencia VAN HOD = 300 Convergencia VAN [m] [s] MUS$ [m] [s] MUS$ [m] [s] MUS$ A 12 2,134 A 15 2,692 A 13 2,535 B 12 2,039 B 13 2,637 B 24 2,532 HOD = 360 Convergencia VAN HOD = 360 Convergencia VAN HOD = 360 Convergencia VAN [m] [s] MUS$ [m] [s] MUS$ [m] [s] MUS$ A 8 2,140 A 9 2,701 A 17 2,546 B 16 2,047 B 14 2,637 B 31 2,532 Convergencia VAN Convergencia VAN Convergencia VAN H=S.R. H=S.R. H=S.R. [s] MUS$ [s] MUS$ [s] MUS$ A 30 2,169 A 23 2,704 A 9 2,559 B 22 2,056 B 0 2,646 B 30 2,532 Tabla 21. Resultados del VAN para ley de corte fija predefinida y dos discretizaciones de cilindros distintos (A y B). Caso Palabora.
38
Para la discretización máxima, promedio y mínima del pit final se obtuvo el siguiente resultado con respecto a la geometría. MAX PIT
MED PIT
MIN PIT
HOD = 120 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 120 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 120 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
360
1,080
2
A
720
1,380
2
A
1,140
-
-
B
360
-
-
B
720
-
-
B
1,140
-
-
HOD = 180 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 180 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 180 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
360
1,440
3
A
720
1,380
3
A
1,200
1,440
3
B
360
1,440
3
B
720
1,260
3
B
1,200
-
-
HOD = 240 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 240 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 240 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
360
1,080
4
A
720
1,140
4
A
1,140
1,440
4
B
360
1,260
4
B
720
1,440
4
B
1,200
1,440
4
HOD = 300 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 300 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 300 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
360
1,080
5
A
720
1,440
5
A
1,200
1,440
4
B
360
1,440
5
B
720
1,440
5
B
1,140
1,440
5
HOD = 360 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 360 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 360 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
360
960
6
A
720
1,440
6
A
1,020
1,260
4
B
360
1,440
6
B
720
1,440
5
B
1,140
1,440
5
H=S.R.
Fondo del Pit [m]
Cota del footprint [m]
HOD real
H=S.R.
Fondo del Pit [m]
Cota del footprint [m]
HOD real
H=S.R.
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
360
1,440
8
A
720
1,140
7
A
960
1,440
8
B
360
1,440
18
B
720
1,440
12
B
1,140
1,440
5
Tabla 22. Resultados de la geometría de transición para ley de corte fija predefinida y dos discretizaciones de cilindros distintos (A y B). Caso Palabora.
39
(b) Resultados para ley de corte fija y dos discretizaciones de cilindros distintos (A y B).
Al correr el modelo sin la restricción de la ley de corte predefinida, fijando la ley de corte elegida para cada método y considerando dos tipos de discretizaciones para la subterránea. Un cilindro que encierra las columnas con ley de corte sobre la elegida (A) y un cilindro que encierra las columnas económicas del modelo de bloques, definidas por el footprint finder (B). Se tienen los siguientes resultados con respecto al VAN. MAX PIT
MED PIT
MIN PIT
HOD = 120 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 120 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 120 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
A
151
2,093
0.40
0.30
A
64
2,667
0.40
0.30
A
179
2,540
0.40
0.30
B
165
1,990
0.40
0.40
B
111
2,609
0.40
0.40
B
190
2,526
0.40
0.40
HOD = 180 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 180 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 180 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
A
157
2,129
0.40
0.50
A
84
2,667
0.40
0.30
A
209
2,551
0.40
0.50
B
232
2,035
0.40
0.60
B
164
2,618
0.40
0.80
B
257
2,520
0.40
0.40
HOD = 240 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 240 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 240 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
A
217
2,157
0.40
0.50
A
80
2,709
0.40
0.50
A
139
2,559
0.40
0.50
B
283
2,059
0.40
0.60
B
284
2,634
0.40
0.60
B
445
2,536
0.40
0.80
HOD = 300 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 300 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 300 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
A
199
2,192
0.40
0.50
A
91
2,723
0.40
0.50
A
257
2,562
0.40
0.50
B
393
2,095
0.40
0.80
B
216
2,670
0.40
0.80
B
398
2,545
0.40
0.80
HOD = 360 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 360 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
HOD = 360 [m]
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
A
192
2,192
0.40
0.50
A
260
2,722
0.40
0.50
A
303
2,564
0.40
0.50
B
280
2,120
0.40
0.80
B
178
2,684
0.40
0.80
B
733
2,550
0.40
0.80
H=S.R.
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
H=S.R.
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
H=S.R.
T.S. [s]
VAN MUS$
LC RAJO %Cu
LC SUBTE %Cu
2,219
0.40
0.50
A
248
2,724
0.40
0.90
A
435
2,564
0.40
0.50
2,161
0.40
0.80
B
227
2,712
0.40
0.90
B
1,9 79
2,554
0.40
0.80
A B
5,1 79 1,3 58
Tabla 23. Resultados del VAN para ley de corte fija y dos discretizaciones de cilindros distintos (A y B). Caso Palabora.
40
Los resultados obtenidos relacionados con la geometría se presentan en la tabla siguiente. MAX PIT MED PIT MIN PIT Fondo Cota del Fondo Cota del Fondo Cota del HOD = HOD HOD = HOD HOD = HOD del Pit footprint del Pit footprint del Pit footprint 120 [m] real 120 [m] real 120 [m] real [m] [m] [m] [m] [m] [m] A 360 1,080 2 A 720 1,080 2 A 1,140 1,260 2 B 360 1,320 2 B 720 1,440 2 B 1,140 1,440 2 Fondo Cota del Fondo Cota del Fondo Cota del HOD = HOD HOD = HOD HOD = HOD del Pit footprint del Pit footprint del Pit footprint 180 [m] real 180 [m] real 180 [m] real [m] [m] [m] [m] [m] [m] A 360 1,260 3 A 720 1,080 2 A 1,020 1,260 3 B 360 1,080 3 B 720 1,440 3 B 1,140 1,320 2 Fondo Cota del Fondo Cota del Fondo Cota del HOD = HOD HOD = HOD HOD = HOD del Pit footprint del Pit footprint del Pit footprint 240 [m] real 240 [m] real 240 [m] real [m] [m] [m] [m] [m] [m] A 360 1,080 4 A 720 1,260 4 A 960 1,320 4 B 360 1,260 4 B 720 1,260 3 B 1,020 1,320 4 Fondo Cota del Fondo Cota del Fondo Cota del HOD = HOD HOD = HOD HOD = HOD del Pit footprint del Pit footprint del Pit footprint 300 [m] real 300 [m] real 300 [m] real [m] [m] [m] [m] [m] [m] A 360 1,440 5 A 720 1,440 5 A 1,080 1,440 5 B 360 1,320 5 B 720 1,440 5 B 1,080 1,440 5 Fondo Cota del Fondo Cota del Fondo Cota del HOD = HOD HOD = HOD HOD = HOD del Pit footprint del Pit footprint del Pit footprint 360 [m] real 360 [m] real 360 [m] real [m] [m] [m] [m] [m] [m] A 360 1,440 5 A 720 1,080 5 A 1,080 1,440 6 B 360 1,140 6 B 720 1,440 6 B 1,020 1,440 6 Fondo Cota del Fondo Cota del Fondo Cota del HOD HOD HOD H=S.R. del Pit footprint H=S.R. del Pit footprint H=S.R. del Pit footprint real real real [m] [m] [m] [m] [m] [m] A 360 1,440 18 A 720 1,440 12 A 1,080 1,440 6 B 360 1,440 9 B 720 1,440 12 B 960 1,440 8 Tabla 24. Resultados de la geometría de transición para ley de corte fija y dos discretizaciones de cilindros distintos (A y B). Caso Palabora.
41
(c) Resultados para ley de corte variable y dos discretizaciones de cilindros distintos (A y B).
Finalmente, utilizando una ley de corte variable para los estratos del rajo y considerando dos tipos de discretizaciones para la subterránea. Un cilindro que encierra las columnas con ley de corte sobre la elegida (A) y un cilindro que encierra las columnas económicas del modelo de bloques, definidas por el footprint finder (B). Se tienen los siguientes resultados con respecto al VAN. MAX PIT
MED PIT
MIN PIT
HOD = 120 [m] T.S. (s) VAN MUS$ HOD = 120 [m] T.S. (s) VAN MUS$ HOD = 120 [m] T.S. (s) VAN MUS$ A
382
2,822
A
381
2,819
A
124
2,617
B
582
2,765
B
500
2,775
B
270
2,599
HOD = 180 [m] T.S. (s) VAN MUS$ HOD = 180 [m] T.S. (s) VAN MUS$ HOD = 180 [m] T.S. (s) VAN MUS$ A
1,603
2,766
A
318
2,821
A
277
2,622
B
635
2,722
B
870
2,793
B
774
2,607
HOD = 240 [m] T.S. (s) VAN MUS$ HOD = 240 [m] T.S. (s) VAN MUS$ HOD = 240 [m] T.S. (s) VAN MUS$ A
647
2,784
A
333
2,836
A
252
2,631
B
3,206
2,736
B
651
2,807
B
606
2,610
HOD = 300 [m] T.S. (s) VAN MUS$ HOD = 300 [m] T.S. (s) VAN MUS$ HOD = 300 [m] T.S. (s) VAN MUS$ A
1,113
2,794
A
637
2,865
A
222
2,635
B
3,637
2,750
B
513
2,806
B
579
2,614
HOD = 360 [m] T.S. (s) VAN MUS$ HOD = 360 [m] T.S. (s) VAN MUS$ HOD = 360 [m] T.S. (s) VAN MUS$ A
2,566
2,800
A
401
2,871
A
143
2,637
B
1,201
2,761
B
1,728
2,833
B
689
2,624
H=S.R.
T.S. (s) VAN MUS$
H=S.R.
T.S. (s) VAN MUS$
H=S.R.
T.S. (s) VAN MUS$
A
11,781
2,806
A
1,671
2,877
A
334
2,644
B
2,481
2,782
B
1,017
2,843
B
564
2,628
Tabla 25. Resultados del VAN para ley de corte variable y dos discretizaciones de cilindros distintos (A y B). Caso Palabora.
42
La secuencia de leyes de corte por estrato obtenida para esta corrida se encuentra en la tabla inferior. MAX PIT HOD = 120 [m] 1
LC A %Cu
LC B %Cu
0.3
0.2
HOD = 180 [m] 1
LC A %Cu
LC B %Cu
0.3
0.2
HOD = 240 [m] 1
LC A %Cu
LC B %Cu
0.3
0.2
HOD = 300 [m] 1
LC A %Cu
LC B %Cu
0.3
0.2
HOD = 360 [m] 1
LC A %Cu
LC B %Cu
H= FREE
0.3
0.2
1
LC A %C u 0.3
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
0.2
0.2
5
0.4
0.4
5
0.4
0.4
5
0.4
0.4
5
0.4
0.4
5
0.4
0.4
5
0.4
0.4
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
12
0.5
0.5
12
0.5
0.5
12
0.5
0.5
12
0.5
0.5
12
0.5
0.5
12
0.5
0.5
13
0.5
0.5
13
0.5
0.5
13
0.5
0.5
13
0.5
0.5
13
0.5
0.5
13
0.5
0.5
14
0.5
0.5
14
0.5
0.5
14
0.5
0.5
14
0.5
0.5
14
0.5
0.5
14
0.5
0.5
15
0.5
0.5
15
0.5
0.5
15
0.5
0.5
15
0.5
0.5
15
0.5
0.5
15
0.5
0.5
16
0.5
0.5
16
0.5
0.5
16
0.5
0.5
16
0.5
0.5
16
0.5
0.5
16
0.5
0.5
17
0.5
0.5
17
0.5
0.5
17
0.5
0.5
17
0.5
0.5
17
0.5
0.5
17
0.5
0.5
18
0.5
0.5
18
0.5
0.5
18
0.5
0.5
18
0.5
0.5
18
0.5
0.5
18
0.5
0.5
19
0.5
0.5
19
0.5
0.5
19
0.5
0.5
19
0.5
0.5
19
0.5
0.5
19
0.5
0.5
LC B %Cu 0.2
20
20
20
20
20
20
0.5
0.9
21
21
21
21
21
21
0.5
0.9
22
0.3
0.4
22
22
22
22
0.9
22
0.5
0.9
23
0.3
0.4
23
23
23
23
0.9
23
0.5
0.9
24
24
24
24
24
0.5
0.9
24
0.5
0.9
25
25
25
25
0.5
0.8
25
0.5
0.9
25
0.5
0.9
26
26
0.5
0.6
26
0.5
0.7
26
0.5
0.8
26
0.5
0.9
26
0.5
0.9
27
27
0.5
0.6
27
0.5
0.7
27
0.5
0.8
27
0.5
0.9
27
0.5
0.9
28
28
0.5
0.6
28
0.5
0.7
28
0.5
0.8
28
0.5
28
0.5
0.9
29
29
29
0.5
0.7
29
0.5
0.8
29
0.5
29
0.5
0.9
30
30 30 30 30 30 Tabla 26. Política de leyes de corte para la extracción de los estratos. Discretización MAX PIT.
43
MED PIT HOD = 120 [m]
LC B %Cu
HOD = 180 [m]
LC B %Cu
HOD = 240 [m]
LC B %Cu
HOD = 300 [m]
LC B %Cu
HOD = 360 [m]
LC A %Cu
LC A %Cu
LC A %Cu
LC A %Cu
LC A %Cu
LC B %Cu
H= FREE
LC A %Cu
LC B %Cu
1
0.4
0.2
1
0.4
0.2
1
0.4
0.2
1
0.4
0.2
1
0.4
0.2
1
0.4
0.2
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
0.2
0.2
5
0.2
0.2
5
0.2
0.2
5
0.2
0.2
5
0.2
0.2
5
0.2
0.2
5
0.2
0.2
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
12
0.4
0.4
12
0.4
0.4
12
0.4
0.4
12
0.4
0.4
12
0.4
0.4
12
0.4
0.4
13
0.4
0.4
13
0.4
0.4
13
0.4
0.4
13
0.4
0.4
13
0.4
0.4
13
0.4
0.4
14
0.4
0.4
14
0.4
0.4
14
0.4
0.4
14
0.4
0.4
14
0.4
0.4
14
0.4
0.4
15
0.4
0.4
15
0.4
0.4
15
0.4
0.4
15
0.4
0.4
15
0.4
0.4
15
0.4
0.4
16
0.4
0.4
16
0.4
0.4
16
0.4
0.4
16
0.4
0.4
16
0.4
0.4
16
0.4
0.4
17
0.4
0.4
17
0.4
0.4
17
0.4
0.4
17
0.4
0.4
17
0.4
0.4
17
0.4
0.4
18
0.5
0.5
18
0.5
0.5
18
0.5
0.5
18
0.5
0.5
18
0.5
0.5
18
0.5
0.5
19
0.5
0.5
19
0.5
0.5
19
0.5
0.5
19
0.5
0.5
19
0.5
0.5
19
0.5
0.8
0.5
20
0.5
0.5
20
0.5
20
20
20
0.5
0.8
21
0.5
0.7
21
0.5
21
21
21
0.5
0.8
20 21 22
0.3
0.4
22
0.7
22
22
22
22
0.5
0.8
23
0.3
0.4
23
0.7
23
23
23
23
0.5
0.8
24
24
24
0.5
0.8
24
0.5
0.8
25
0.5
0.9
25
0.5
0.8
25
0.5
0.8
24
24
25
25
0.4
25
26
26
0.4
26
0.5
0.7
26
0.5
0.9
26
0.5
0.8
26
0.5
0.8
27
27
0.4
27
0.5
0.7
27
0.5
0.9
27
0.5
0.8
27
0.5
0.8
28
28
28
0.5
0.7
28
0.5
0.9
28
0.5
0.8
28
0.5
0.8
29
29
29
0.5
0.7
29
0.5
0.9
29
0.5
0.8
29
0.5
0.8
30
30 30 30 30 30 Tabla 27. Política de leyes de corte para la extracción de los estratos. Discretización MED PIT
44
HOD = 120 [m] 1
LC A %Cu
LC B %Cu
0.2
0.2
HOD = 180 [m] 1
2
0.2
0.2
3
0.2
4
LC A %Cu
LC B %Cu
0.2
MIN PIT HOD LC B = %Cu 300 [m] 0.2 1
0.2
0.2
HOD = 360 [m] 1
LC A %Cu
LC B %Cu
0.2
0.2
HOD = 240 [m] 1
2
0.2
0.2
2
0.2
0.2
2
0.2
0.2
0.2
3
0.2
0.2
3
0.2
0.2
3
0.2
0.2
0.2
4
0.2
0.2
4
0.2
0.2
4
5
0.2
0.2
5
0.2
0.2
5
0.2
0.2
6
0.4
0.4
6
0.4
0.4
6
0.4
7
0.4
0.4
7
0.4
0.4
7
8
0.4
0.4
8
0.4
0.4
9
0.4
0.4
9
0.4
10
0.4
0.4
10
11
0.4
0.4
12
0.4
13
LC A %Cu
LC B %Cu
H= FREE
LC A %Cu
LC B %Cu
0.2
0.2
1
0.2
0.2
2
0.2
0.2
2
0.2
0.2
0.2
3
0.2
0.2
3
0.2
0.2
0.2
0.2
4
0.2
0.2
4
0.2
0.2
5
0.2
0.2
5
0.2
0.2
5
0.2
0.2
0.4
6
0.4
0.4
6
0.4
0.4
6
0.4
0.4
0.4
0.4
7
0.4
0.4
7
0.4
0.4
7
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
8
0.4
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
9
0.4
0.4
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
10
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
11
0.4
0.4
0.4
12
0.4
0.4
12
0.4
0.4
12
0.4
0.4
12
0.4
0.4
12
0.4
0.4
0.4
0.4
13
0.4
0.4
13
0.4
0.4
13
0.4
0.4
13
0.4
0.4
13
0.4
0.4
14
0.4
0.4
14
0.4
0.4
14
0.4
0.4
14
0.4
0.4
14
0.4
0.4
14
0.4
0.4
15
0.4
0.4
15
0.4
0.4
15
0.4
0.4
15
0.4
0.4
15
0.4
0.4
15
0.4
0.4
16
0.4
0.4
16
0.4
0.4
16
0.4
0.4
16
0.4
0.4
16
0.4
0.4
16
0.4
0.4
17
0.4
0.4
17
0.4
0.4
17
0.4
0.4
17
0.4
0.4
17
0.4
0.4
17
0.4
0.4
18
0.4
0.4
18
0.4
0.4
18
0.4
0.4
18
0.4
0.4
18
0.4
0.4
18
0.4
0.4
19
0.4
0.4
19
0.4
0.4
19
0.4
0.4
19
0.4
0.4
19
0.4
0.4
19
0.4
0.4
20
0.4
0.4
20
0.4
0.4
20
0.4
0.4
20
0.4
0.4
20
0.4
0.4
20
0.4
0.4
21
0.4
0.4
21
0.4
0.4
21
0.4
0.4
21
0.4
0.4
21
0.4
0.4
21
0.4
0.4
22
0.4
0.4
22
0.4
22
0.4
0.4
22
0.4
0.4
22
0.4
0.4
22
0.5
0.9
0.4
23
0.4
23
0.4
23
0.4
0.4
23
0.4
23
0.5
0.9
0.4
24
0.4
24
0.5
24
0.5
0.9
25
0.8
25
0.5
0.8
25
0.5
0.9
23 24
24
0.5
24 25
25
0.3
0.4
25
0.5
26
0.3
0.4
26
0.5
LC A %Cu
0.6
26
0.5
0.9
26
0.5
0.8
26
0.5
0.8
26
0.5
0.9
27
27
0.6
27
0.5
0.9
27
0.5
0.8
27
0.5
0.8
27
0.5
0.9
28
28
0.6
28
0.5
0.9
28
0.5
0.8
28
0.5
0.8
28
0.5
0.9
29
29
29
0.5
0.9
29
0.5
0.8
29
0.5
0.8
29
0.5
0.9
30
30 30 30 30 30 Tabla 28. Política de leyes de corte para la extracción de los est ratos. Discretización MIN PIT
45
El resultado de la corrida anterior con respecto a la geometría se encuentra contenido en la siguiente tabla. MAX PIT
MED PIT
MIN PIT
HOD = 120 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 120 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 120 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
840
1,080
2
A
840
1,080
2
A
1,020
1,260
2
B
840
1,080
2
B
900
1,080
2
B
1,080
1,260
2
HOD = 180 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 180 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 180 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
840
1,380
3
A
960
1,320
3
A
960
1,260
2
B
840
1,380
3
B
900
1,080
3
B
1,080
1,380
3
HOD = 240 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 240 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 240 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
840
1,440
4
A
960
1,440
4
A
1,140
1,440
4
B
840
1,440
4
B
840
1,440
4
B
1,020
1,440
3
HOD = 300 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 300 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 300 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
840
1,440
5
A
840
1,440
5
A
1,080
1,440
4
B
840
1,440
5
B
840
1,440
5
B
1,140
1,440
4
HOD = 360 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 360 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
HOD = 360 [m]
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
840
1,440
6
A
840
1,440
6
A
1,080
1,440
6
B
840
1,320
6
B
840
1,440
6
B
1,020
1,440
5
H=S.R.
Fondo del Pit [m]
Cota del footprint [m]
HOD real
H=S.R.
Fondo del Pit [m]
Cota del footprint [m]
HOD real
H=S.R.
Fondo del Pit [m]
Cota del footprint [m]
HOD real
A
780
1,440
10
A
840
1,440
6
A
960
1,440
6
B
840
1,440
10
B
780
1,440
9
B
960
1,440
7
Tabla 29. Resultados de la geometría de transición para ley de corte variable y dos discretizaciones de cilindros distintos (A y B). Caso Palabora.
46
iv) Consolidación de resultados para el método ILP.
Para cada uno de los escenarios definidos por los tipos de agregación y criterios de aplicación de leyes de corte para los estratos, se calcula un pit final, la elevación del footprint y la altura de columna para el block caving. De esta manera, las tres configuraciones óptimas factibles de estas metodologías son las siguientes: Metodología VAN [MUS$] Discretización Rajo Discretización Subte HOD [m] (a) 2,701 MED PIT A 360 (b) 2,723 MED PIT A 300 (c) 2,871 MED PIT A 360 Tabla 30. Resultados con respecto al VAN para la corrida del modelo de optimización. Caso Palabora.
Claramente se observa que la elección del tipo de agregación del rajo como de la subterránea no es una decisión trivial. La máxima diferencia porcentual que se alcanzan por el concepto de tener una agregación distinta del rajo para cada discretización subterránea y altura de columna se detalla a continuación. Metodología Discretización Subte HOD = 120 [m] HOD = 180 [m] HOD = 240 [m] HOD = 300 [m] HOD = 360 [m]
(a) A 24 % 22 % 21 % 21 % 21 %
B 25 % 24 % 23 % 23 % 22 %
(b) A 22 % 20 % 20 % 20 % 19 %
B 24 % 22 % 22 % 22 % 21 %
(c) A 7% 7% 7% 8% 8%
B 6% 7% 7% 7% 7%
Tabla 31. Variación máxima porcentual del VAN debido al cambio en la discretización del rajo. Caso Palabora.
Considerando ahora la variación de la altura de columna para una discretización del rajo elegida, se tienen las siguientes diferencias máximas porcentuales. Metodología Discretización Subte MAX PIT MED PIT MIN PIT A 7% 3% 1% (a) B 5% 2% 1% A 4% 2% 1% (b) B 6% 3% 1% A 2% 2% 1% (c) B 2% 2% 1% Tabla 32. Variación máxima porcentual del VAN debido al cambio en la altura de columna máxima extraíble. Caso Palabora
47
Por otro lado, además de tener en cuenta el valor económico del proyecto de transición, también es importante analizar las cotas que determinan la transición, es decir, el fondo del pit y el footprint. En la tabla a continuación se presentan estos valores. Metodología Fondo del Pit [m] Cota del footprint [m] HOD (a) -720 -1440 360 (b) -720 -1440 300 (c) -840 -1440 360 Tabla 33. Resultados respecto a la geometría de la transición obtenidos a parti r del modelo de optimización. Caso Palabora.
Al igual que el valor económico del proyecto, las cotas anteriormente presentadas se ven modificadas por efecto de la agregación del pit. A continuación se presenta la máxima diferencia en metros de la cota del fondo del pit y la elevación del footprint para cada combinación de metodología, agregación subterránea y altura de columna. Máxima diferencia de la cota del fondo de Pit Metodología Discretización Subte
(a) A
(b) B
A
(c) B
A
B
H=2
780 780 780 780 180 240
H=3
840 840 660 780 240 240
H=4
780 840 600 660 180 180
H=5
840 780 720 720 300 300
H=6 660 780 720 660 180 180 Tabla 34. Máxima diferencia de la cota del fondo de pit sujeto a un cambio en la discretización del rajo. Caso Palabora.
Máxima diferencia de la cota del footprint Metodología Discretización Subte
(a) A
(b) B
H=2
A
(c) B
A
B
180 120 180 180
H=3
60
180 360 120 300
H=4
360 180 240
60
0
0
H=5
360
120
0
0
0
0
H=6 300 0 360 300 0 120 Tabla 35. Máxima diferencia de la cota del footprint sujeto a un cambio en la discretización del rajo. Caso Palabora.
48
Por otro lado, al variar la restricción de la altura de columna también se generan diferencias en las cotas del fondo del pit y la elevación del footprint. Máxima diferencia de la cota del fondo de Pit Discretización Metodología MAX PIT MED PIT MIN PIT Subte A 0 0 180 (a) B 0 0 60 A 0 0 180 (b) B 0 0 180 A 60 120 180 (c) B 0 120 180 Tabla 36. Máxima diferencia de la cota del fondo de pit sujeto a un cambio en la altura de columna máxima extraíble. Caso Palabora. Máxima diferencia de la cota del footprint Discretización Metodología MAX PIT MED PIT MIN PIT Subte A (a) B A 360 360 180 (b) B 240 180 120 A 360 360 180 (c) B 360 360 180 Tabla 37. Máxima diferencia de la cota del footprint sujeto a un cambio en la altura de columna máxima extraíble. Caso Palabora
Además del tipo de agregación y de la altura de columna máxima admisible es posible generar una sensibilidad respecto a parámetros económicos, tales como el precio del cobre o el ritmo de producción de los métodos utilizados. En las tablas puestas a continuación se presenta cada uno de los análisis antes mencionados para cada uno de los casos óptimos factibles definidos anteriormente. Variación respecto al precio de evaluación del Cu. Variación del Precio de evaluación Metodología -20% -10% 10% 20% (a) -69% -34% 33% 69% (b) -68% -34% 34% 69% (c) -68% -34% 35% 71% Tabla 38. Análisis de sensibilidad para el VAN con respecto a la variación del precio. Caso Palabora.
49
Variación del Precio
Fondo del Pit
Metodología (a)
Cota del Footprint -20%
-720
-1,260
HOD [m] (c)
-720
360 -720
Fondo del Pit
-1,440
Cota del footprint
-720
-1,260
-720
-1,440
-1,140
Cota del footprint
-720
-1,440
-1,440
-720
-1,440
240 -720
-1,080
300 -840
Cota del footprint 20%
-1,260
360 -720
Fondo del Pit
10%
-1,440
300 -780
Fondo del Pit
0%
360
300 -720
Cota del Footprint -10%
-1,440
HOD [m] (b)
Fondo del Pit
360 -720
-1,440
300 -900
-1,440
300 -1,080
-1,440
HOD [m] 300 360 360 360 360 Tabla 39. Análisis de sensibilidad para la geometría de transición con respecto a la variación del precio. Caso Palabora.
Variación con respecto al ritmo de producción del rajo. Variación del ritmo de producción del rajo Metodología -50% -25% 25% (a) -82% -82% -12% (b) -44% -20% -12% (c) -40% -17% -9%
50% -30% -2% -5%
Tabla 40. Análisis de sensibilidad para el VAN con respecto a la variación del ritmo de producción del rajo. Caso Palabora.
Variación del ritmo de producción del rajo
Fondo del Pit
Metodología (a) HOD [m]
HOD [m]
HOD [m]
-1,440
-180
Fondo del Pit
-1,440
-840
-1,440 360
-1,440
-720
Fondo del Pit
-1,260
-720
-1,140
-1,440
-1,080
300
Fondo del Pit
-1,440
-1,080
-1,440
-1,440
-1,080
-1,440
360
-1,440 300
-1,440 360
-1,140
300 -1,140
Cota del footprint 50%
360
300 -840
Cota del footprint 25%
360
300 -840
Cota del footprint 0%
360
300 -720
Cota del Footprint -25%
360 -720
(c)
Fondo del Pit
-50% -180
(b)
Cota del Footprint
-1,440 300
-1,140
-1,440 300
Tabla 41. Análisis de sensibilidad para la geometría de transición con respecto a la variación del ritmo del producción del rajo. Caso Palabora.
Variación con respecto al ritmo de producción del block caving. Variación del ritmo de producción del block caving Metodología -50% -25% 25% (a) -3% -1% 1% (b) -3% -2% 1% (c) -3% -1% 1%
50% 3% 2% 2%
Tabla 42. Análisis de sensibilidad para el VAN con respecto a la variación del ritmo de producción del block caving. Caso Palabora.
50
Variación del ritmo de producción del block caving
Fondo del Pit
Metodología (a)
-720
-720
HOD [m]
-1,440
-720
Fondo del Pit
-1,440
-720
-1,260 360
-1,440
-720
Fondo del Pit
-1,260
-720
-1,440
-1,440
-720
360
Fondo del Pit
-1,440
-720
-1,440 360
-1,080
-720
-1,440
-1,200 300
-1,440 360
-720
300 -840
Cota del footprint 50%
360
300 -840
Cota del footprint 25%
360
300 -840
Cota del footprint 0%
360
180 -900
Cota del Footprint -25%
300
HOD [m] (c)
Fondo del Pit
-50%
HOD [m] (b)
Cota del Footprint
-1,440 300
-900
-1,260 360
Tabla 43. Análisis de sensibilidad para la geometría de transición con respecto a la variación del ritmo del producción del block caving. Caso Palabora.
De esta manera, se observa que el proyecto de transición eleva su valor si es que el precio de evaluación se aumenta, ya que permite extraer una mayor cantidad de mineral que antes no era económicamente extraíble, variando así el diseño de la transición. Con respecto al análisis de sensibilidad relacionado a los ritmos de producción de los métodos implicados, resulta que el valor del proyecto es más sensible a un cambio en el ritmo de producción del rajo, dado principalmente porque a un tratamiento menor, el pit final económico será de menor tamaño y por ende de un valor menor. Por otra parte, un mayor tratamiento implica una mayor inversión inicial, la cual también perjudica el VAN del proyecto.
51
5.4. Comparación entre las alternativas de transición para Palabora. El criterio que será utilizado para comparar cada una de las metodologías corresponde al VAN asociado al proyecto completo. Sin embargo, los VAN reportados por el modelo de optimización permiten realizar comparaciones y análisis sólo entre los resultados de esta metodología, pero no son comparables con los obtenidos por las otras dos metodologías debido a que los supuestos y las bases de cálculo de las valorizaciones entregadas al modelo son distintas a las utilizadas en los métodos tradicionales. De esta manera, para hacer comparables todas las alternativas de transición considerando el VAN, se deben evaluar los resultados entregados por el método ILP con la valorización de estratos utilizada en la metodología incremental por fases, respetando las restricciones de las leyes de corte establecidas por el modelo. Considerando lo anterior, se presentan los VAN para el rajo y para el block caving, actualizados al inicio del proyecto, es decir, el valor de la subterránea se encuentra descontado considerando la extracción previa del pit. Valor del proyecto
VAN rajo [MUS$] VAN Block Caving[MUS$] VAN del proyecto [MUS$]
Secuencial
3,593
27
3,620
Incremental por fases
3,580
72
3,653
ILP
3,754
124
3,878
Tabla 44. Comparación entre metodologías según VAN. Caso Palabora.
Por otro lado, se tienen los resultados asociados al diseño que debe ser implementado en la transición propuesta por cada metodología. Diseño de la transición
Cota fondo del rajo Cota del Footprint
Secuencial
-1,147
-1,504
Incremental por fases
-964
-1,264
ILP
-840
-1,440
Tabla 45. Comparación entre metodologías según el diseño de la transición. Caso Palabora.
Se observa que existe una coincidencia entre las metodologías para un cierto set de estratos que deben ser minados por el rajo. Por lo tanto, es posible determinar el valor económico que está en competencia, entre el rajo y el block caving, para los estratos por debajo de la última cota en común para el rajo en las tres metodologías. Este valor económico en competencia es determinado al evaluar las últimas fases del rajo y sumarle el valor del proyecto subterráneo. Dicho valor se encontrará descontado sólo por el tiempo que requiera extraer las últimas fases del rajo más el de la subterránea.
52
Valor económico en Últimas fases del rajo Extracción rajo Valor del BC Extracción competencia [MUS$] [años] [MUS$] BC [años] Secuencial 85 11 215 23 Incremental por fases 34 4 648 32 ILP 0 0 839 14 Tabla 46. Valor económico en competencia entre el rajo y la subterránea. Caso Palabora.
Total [MUS$] 300 682 839
De esta manera, se observa que el máximo valor en competencia es obtenido en la opción ILP, donde no existe un desarrollo de fases marginales para el rajo, por lo que los mejores valores económicos entregados por la subterránea se encuentran afectados por un menor descuento. Por lo tanto, la metodología que entrega un mejor resultado con respecto al VAN del proyecto corresponde a la ILP, logrando una mejora promedio de 7% con respecto a las otras dos metodologías.
53
6. Información disponible para el caso de aplicación en Andina Sur-Sur. 6.1. Inventario de recursos. En el caso del modelo de bloques de Andina Sur-Sur se tiene lo siguiente: Tamaño del Bloque: X: 15 m, Y: 15 m, Z: 16 m. Numero de Bloques: X: 312, Y: 331, Z: 121. Densidad promedio: 2.65 [t/m 3] La curva tonelaje ley para el modelo de bloques de Andina Sur-Sur está expresada en la siguiente tabla.
Gráfico 4. Curva tonelaje – ley para Andina Sur-Sur.
6.2. Parámetros técnico-económicos para Andina Sur-Sur. Se considera que el ángulo de talud global para el pit final que se diseña corresponde a un ángulo constante para todas las paredes del rajo, con un valor de 45°. Se asume también una dilución minera del 5% y una recuperación minera del 95%. Además no se aplicará una estructura de ajuste sobre el costo mina (CAF), así como tampoco el costo de proceso (PCAF). Se utilizará una serie de revenue factors en el software Whittle entre 0.4 y 1.0, con un paso de 0.02, reproduciendo 31 valores distintos. La tasa de descuento utilizada para el cálculo del VAN será fija con un valor del 10%.
54
Para realizar la evaluación de los estratos a ser minados por el método rajo se consideran los siguientes parámetros económicos. Parámetros económicos para el rajo. Precio 1.5 US$/lb Costo de 5.5 US$/t Proceso Costo Mina 1.8 US$/t Recuperación 85 % Costo de Venta 0.32 US$/lb Tabla 47. Parámetros económicos para el rajo. Caso Andina Sur-Sur.
Además, el ritmo de producción de la mina subterránea alcanza las 23,000 tpd. Este valor se encuentra al definir que aproximadamente se tendrán 1,008 columnas en el footprint. De esta manera, dado un vertical mining rate de 14, donde este parámetro corresponde a los metros de roca que se extraerán aproximadamente en la vertical sobre cada una de las columnas, se tiene que el ritmo de producción anual por columna alcanza a 8,348 tpa/col, lo cual es equivalente a una producción de 23,000 tpd. Con respecto a los parámetros económicos utilizados en las corridas de footprint finder, se considerará para el desarrollo de la mina subterránea un costo de US$ 2,000 por metro cuadrado. Para realizar la evaluación de los estratos a ser minados por el método subterráneo, se consideran los siguientes parámetros económicos. Parámetros económicos para el block caving. Precio 1.5 US$/lb Costo de 5.5 US$/t Proceso Costo Mina 5 US$/t Recuperación 85 % Costo de Venta 0.32 US$/lb Tabla 48. Parámetros económicos para el block caving. Caso Andina Sur-Sur.
6.3. Datos de entrada para el modelo de optimización. Para poder correr el modelo de optimización, deben ser cuidadosamente definidos cada uno de los sets y parámetros definidos en él. A continuación, se presentan los valores constantes para cada una de las formas en que se corrió el modelo.
Discount rate = 10% 55
Para el caso de la ley de corte variable, se utilizó el siguiente vector de leyes de corte.
Un parámetro importante para el cálculo de la transición corresponde a la altura de columna máxima extraíble. Para dicho parámetro, se definió el siguiente vector, además de considerar que no existía restricción.
Finalmente, el modelo fue desarrollado en el lenguaje AMPL y se resolvió usando CPLEX Parallel 10.2 con 4 procesos. Los tiempos de solución fueron calculados en base a un computador con 4 procesadores (IntelXeon Cuadcore
[email protected] Ghz, 4 Gb RAM).
56
7. Análisis de la transición para el caso de aplicación en Andina Sur-Sur. Lo que se espera obtener de cada una de las metodologías corresponde a una geometría asociada al complejo minero, esto es, una profundidad máxima para el rajo y una cota en la cual se debe ubicar el footprint de la mina subterránea. Adicionalmente, se obtendrán los VAN de cada una de éstas, parámetro que funcionará como criterio de comparación. 7.1. Resultados para la metodología secuencial. Cargado el modelo de bloques en Whittle y habiendo ingresado los parámetros técnicoeconómicos definidos para esta corrida y la serie de revenue factors, se obtienen las siguientes curvas de VAN para el rajo variando el ritmo de producción de la mina. Valorización de pits a diferentes ritmos de produ cción. 2,500
2,000
1,500
$ S U M N A V 1,000
500
0 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
Pits 80000 tpd
100000 tpd
120000 tpd
140000 tpd
160000 tpd
180000 tpd
200000 tpd
220000 tpd
Gráfico 5. Valorización de los pits Whittle. Caso Andina Sur-Sur.
57
230000 tpd
250000 tpd
30
31
A partir del análisis de los anteriores resultados se encuentra que el ritmo de producción que maximiza el VAN del proyecto es de 230,000 tpd. A continuación se detalla el valor de cada uno de los pits para el ritmo antes mencionado. Pit Final VAN [MUS$] Pit Final VAN [MUS$] 1 -2,339 17 1,229 2 -2,103 18 1,549 3 -2,029 19 1,680 4 -1,767 20 1,771 5 -1,686 21 1,862 6 -1,555 22 1,945 7 -1,345 23 1,972 8 -1,055 24 1,981 9 -812 25 1,985 10 -621 26 1,986 11 -507 27 1,970 12 -291 28 1,970 13 44 29 1,970 14 178 30 1,969 15 546 31 1,969 16 872 Tabla 49. Valorización de los estratos extraídos vía rajo. Caso Andina Sur-Sur.
De esta manera, se define que el pit final para el rajo corresponde al número 26, que posee las siguientes características: Fondo del pit VAN Mineral Vida útil Ley Media Estéril [Mt] REM [m] [MUS$] [Mt] [años] [%Cu] 26 2,971 1,986 4,776 3,979 57 4 0.483 Tabla 50. Características del pit final elegido para la metodología secuencial en Andina Sur-Sur.
Final Pit
Posteriormente, la siguiente etapa del procedimiento corresponde a realizar una corrida de footprint finder, considerando la extracción del rajo, para determinar el nivel en el cual debe ser ubicada la base de la mina subterránea.
58
A continuación se detallan estos resultados para cada una de las cotas de la pata de los bloques del modelo. Cota 3,340 3,324 3,308 3,292 3,276 3,260 3,244 3,228 3,212 3,196 3,180 3,164 3,148 3,132 3,116 3,100 3,084 3,068 3,052 3,036 3,020 3,004 2,988 2,972 2,956 2,940
Tonelaje [kt] 49 364 2,284 3,544 4,803 5,974 7,212 8,696 10,546 12,160 13,566 16,322 18,572 21,029 23,487 31,888 40,515 44,326 46,300 49,672 978 1,874 4,222 8,597 10,624 62,949
VAN [kUS$] 7 242 1,789 3,881 6,300 8,989 11,955 14,859 17,825 20,821 24,277 28,239 32,380 35,940 38,254 42,993 47,236 51,190 55,429 60,247 3,741 4,189 6,029 10,092 14,186 100,096
Área Max H Ley media Tonelaje VAN Área Cota [m2] [m] [%Cu] [kt] [kUS$] [m2] 225 160 1.04 2,924 68,555 109,102 101,925 1,350 176 1.00 2,908 81,198 123,952 125,325 7,425 192 0.97 2,892 106,169 166,711 168,750 10,800 208 0.98 2,876 120,073 193,530 187,875 13,725 224 0.98 2,860 127,094 213,961 194,400 15,975 240 0.98 2,844 136,774 234,603 204,300 17,775 256 0.98 2,828 146,749 251,833 214,650 20,250 272 0.98 2,812 157,907 271,597 225,000 23,400 288 0.98 2,796 138,544 255,471 196,875 25,650 288 0.97 2,780 221,357 375,152 313,425 27,225 288 0.97 2,764 233,930 404,400 325,800 31,725 288 0.97 2,748 249,710 434,656 342,675 34,875 288 0.97 2,732 24,495 35,728 31,500 38,250 288 0.96 2,716 374,904 588,734 505,575 40,725 288 0.95 2,700 18,437 28,489 23,625 50,625 288 0.91 2,684 430,692 683,530 570,825 61,650 288 0.88 2,668 448,590 715,337 593,775 66,825 288 0.88 2,652 540,769 827,054 717,075 68,175 288 0.89 2,636 562,689 884,036 744,525 73,125 288 0.89 2,620 607,327 956,218 802,125 5,175 288 1.36 2,604 631,507 1,001,623 832,725 7,650 288 1.17 2,588 672,512 1,058,612 894,600 12,825 288 1.02 2,572 701,150 1,115,146 932,175 23,625 288 0.97 2,556 730,623 1,176,763 968,400 26,775 288 0.97 2,540 748,447 1,236,899 987,975 92,925 288 0.92 2,524 762,623 1,297,800 1,003,275 Tabla 51. Corrida de footprint finder para Andina Sur-Sur.
Max H Ley media [m] [%Cu] 288 0.92 288 0.91 288 0.92 288 0.92 288 0.92 288 0.93 288 0.92 288 0.92 288 0.93 288 0.91 288 0.91 288 0.91 288 0.91 288 0.90 288 0.91 288 0.90 288 0.90 288 0.89 288 0.89 288 0.90 288 0.90 288 0.90 288 0.90 288 0.90 288 0.91 288 0.91
De esta manera se encuentra que la cota del footprint a la cual se maximiza el valor de la mina subterránea extraída por block caving es la 2,524. En dicha cota se reporta un VAN de MUS$ 1,298, sin embargo, éste valor debe ser castigado considerando el descuento del tiempo de extracción de la mina rajo previamente, que corresponde a 57 años. De esta manera, se obtiene que el VAN de la opción subterránea descontado asciende a MUS$ 5.7. Es importante señalar también que el costo de las infraestructuras mayores de la mina subterránea no se encuentran consideradas en éste análisis.
59
En resumen, los resultados de la metodología secuencial para el modelo de Andina SurSur se encuentran en la tabla siguiente, considerando que el valor de la subterránea se encuentra descontada a la extracción del rajo. Pit Final N° 26 Fondo del Pit [m] VAN [MUS$] Mineral [Mt] Estéril [Mt] Vida útil [años] 2,971 1,986 4,776 3,979 57 Footprint Nivel del FP [m] VAN [MUS$] Mineral [Mt] Ley de Cu % Vida útil [años] 2,524 5.7 763 0.91 89 Proyecto de transición VAN [MUS$] 1,992
Tabla 52. Resultados para la metodología secuencial en Andina Sur-Sur.
Se presenta además una vista isométrica esquemática de la geometría final del proyecto.
Figura 5. Esquema del diseño de la transición según la metodología secuencial para Andina Sur-Sur.
Por lo tanto, la profundidad del rajo alcanza los 1041 [m], mientras que el piso de la mina subterránea se encuentra ubicada en la cota 2524 con una altura de columna máxima de extracción de 300 [m].
60
7.2. Resultados para la metodología incremental por fases. Para esta metodología, se estandariza la cota de los estratos según el fondo de los pits obtenidos desde las corridas de Whittle, los cuales se presentan en la siguiente tabla. Pit 1 2 3 4 5 6 7 8 9 10 11 12
Cota [m] 4,204 3,932 3,900 3,868 3,820 3,772 3,740 3,676 3,612 3,548 3,484 3,420
Pit 13 14 15 16 17 18 19 20 21 22 23
Cota [m] 3,356 3,196 3,132 3,068 3,004 2,972 2,684 2,620 2,588 2,556 2,524
Tabla 53. Fondo de pits Whittle. Caso Andina Sur-Sur.
Sin embargo, dado que entre algunos de los estratos quedaban diferencias de cota sobre los 150 m., se cambió la serie de revenues factors, de manera de refinar los fondos de los pits. Con esto se fijaron los siguientes estratos: Estrato 1 2 3 4 5 6 7 8 9 10 11 12
Cota [m] 4,012 3,980 3,948 3,884 3,852 3,820 3,788 3,756 3,660 3,596 3,532 3,452
Estrato 13 14 15 16 17 18 19 20 21 22 23
Cota [m] 3,404 3,372 3,340 3,308 3,180 3,036 3,020 2,924 2,860 2,604 2,524
Tabla 54. Cota estandarizada para los estratos extraíbles vía rajo. Caso Andina Sur-Sur.
Considerando los pits anteriormente analizados, se define la valorización de cada uno de los estratos al ser extraídos por el método rajo. Cota [m]
VAN del estrato rajo [MUS$]
4,000 3,968 3,936 3,872 3,840 3,808 3,776 3,744 3,648 3,584 3,520 3,392
-1,686,000,000 131,500,000 209,500,000 289,900,000 243,036,881 304,573,382 216,431,654 469,483,070 368,019,268 325,422,146 357,460,590 319,176,794
Cota [m]
VAN del estrato rajo [MUS$]
3,360 3,328 3,296 3,168 3,024 3,008 2,912 2,880 2,848 2,592 2,524
131,734,843 90,375,481 91,782,447 91,782,447 27,730,725 8,908,304 3,205,171 1,792,531 875,348 -8,894,646 -97,021
Tabla 55. Valorización de los estratos extraídos vía rajo. Caso Andina Sur-Sur.
61
Para realizar la valorización de los estratos por el método subterráneo se procedió a correr un footprint finder sin restricción de altura de columna y utilizando el modelo de bloques completo. Sin embargo, se presentaron diferencias notables entre pares de estratos contiguos. Cota [m] 4,012 3,980 3,948 3,884 3,852 3,820 3,788 3,756 3,660 3,596 3,532 3,404 3,372 3,340 3,308 3,180
Mineral [kt] 20,744 31,433 44,273 67,736 85,649 97,010 111,251 123,939 195,213 216,447 276,717 386,410 405,422 419,443 441,410 603,947
Ley [%Cu] 1.65 1.59 1.53 1.42 1.34 1.30 1.27 1.25 1.13 1.13 1.09 1.05 1.05 1.05 1.05 1.00
VAN [kUS$] 216,541,776 287,426,784 349,372,640 381,210,720 376,714,816 371,109,504 366,700,448 350,596,576 366,751,264 377,816,544 435,617,824 520,492,864 535,306,816 531,678,752 548,037,248 616,987,072
Área [m2] 75,150 87,750 102,150 112,275 127,800 129,150 133,875 135,675 170,550 166,725 189,450 220,725 221,400 219,375 221,175 269,775
3,036
1,047,563
0.89
711,057,728
476,550
3,020
1,061,067
0.89
697,270,784
478,350
2,924
605,230
0.89
457,968,640
263,925
2,860
635,498
0.87
449,069,504
271,350
2,604
2,346,299
0.78
1,281,557,760
918,675
2,524
2,882,696
0.79
1,603,775,744
1,057,275
Tabla 56. Corrida de footprint finder para Andina Sur-Sur.
Esto evidencia la presencia de diferentes niveles. De esta manera, se aplicó la metodología a cada nivel por separado, para así obtener el valor de cada incremento. De este modo, se obtiene la siguiente tabla para los valores de los estratos extraídos vía subterránea, realizando todas las consideraciones económicas asociadas. VAN del estrato BC VAN del estrato BC Cota [m] [MUS$] [MUS$] 4,012 288.22 3,404 1,198 3,980 99.63 3,372 127 3,948 96.88 3,340 108 3,884 158.02 3,308 131 3,852 74.61 3,180 2,882 3,820 31.51 3,036 408 3,788 72.60 3,020 50 3,756 18.96 2,924 22 3,660 435.57 2,860 302 3,596 90.41 2,604 418,329 3,532 471.08 2,524 46,313 3,452 508.44 Tabla 57. Valorización de los estratos extraídos vía subterránea. Caso Sur Sur. Cota [m]
62
Ahora bien, para que los valores de los estratos minados por vía rajo y subterránea sean comparables, deben encontrarse descontados al inicio del proyecto. Esto implica que el valor de los estratos extraídos por la mina subterránea deben encontrarse descontados considerando la previa extracción del rajo. A continuación se presenta el escenario en el cual se cumple la transición rajo subterránea asociado a las restricciones de la subterránea y la profundidad del rajo. Flujo de caja Flujo de caja Flujo de caja descontado Cota Flujo de caja descontado descontado para el rajo descontado para el rajo para la subterránea [kUS$] [m] para la subterránea [kUS$] [kUS$] [kUS$] 4,012 -1,686,000 3 3,404 131,735 783 3,980 131,500 1 3,372 90,375 154 3,948 209,500 1 3,340 91,782 158 3,884 289,900 3 3,308 91,782 233 3,852 243,037 2 3,180 27,731 14,890 3,820 304,573 1 3,036 8,908 6,271 3,788 216,432 2 3,020 3,205 0 3,756 469,483 1 2,924 1,793 0 3,660 368,019 23 2,860 875 0 3,596 325,422 8 2,604 -8,895 0 3,532 357,461 120 2,524 -97 0 3,452 319,177 130 Tabla 58. Flujos de caja descontados para los estratos extraídos por rajo y subterránea. Caso Andina Sur-Sur. Cota [m]
Por lo tanto, se define que el rajo debe ser explotado hasta la cota 3,308, para luego pasar a una explotación de block caving, con una altura de columna de 272 m y su footprint ubicado en la cota 3,036. De esta manera se alcanza una profundidad del rajo de 704 [m]. Los valores económicos para las dos minas antes descritas se encuentran en la tabla dispuesta a continuación. Pit Final Fondo del Pit [m] VAN [MUS$] Mineral [Mt] Estéril [Mt] Vida útil [años] 3,308 1,954 3,354 3,230 40 Footprint Nivel del FP [m] VAN [MUS$] Mineral [Mt] Ley de Cu % Vida útil [años] 3,036 21 227 0.99 26 Proyecto de transición VAN [MUS$] 1,975
Tabla 59. Resultados para la metodología incremental por fases en Andina S ur-Sur.
63
A continuación se muestra un esquema del diseño del complejo minero obtenido a través de la metodología incremental por fases.
Figura 6. Esquema del diseño de la transición según la metodología incremental por fases para el Sur Sur
64
7.3. Resultados para la metodología de transición ILP en Andina Sur-Sur. 7.3.1. Aspectos previos.
Para el caso de la mina Sur-Sur, se construyeron discretizaciones de las mismas características para el rajo que las construidas en el caso de estudio de Palabora, es decir, agregaciones para el máximo, mínimo y el pit promedio, conservando un ángulo de talud constante de 45°. Con respecto a la discretización de la mina subterránea, se construyó un cilindro que encierra las columnas económicas obtenidas a partir de la corrida de footprint finder utilizada para el desarrollo de las metodologías tradicionales. i) Escenarios para definir la transición con el modelo ILP.
Para la prueba que se realizará con el modelo, se considerará la siguiente valorizacion, decisiones del modelo y configuración de parámetros.
65
ii) Valorización de los estratos.
Cada una de las discretizaciones construidas son cubicadas en GEMS y posteriormente son valorizados siguiendo los mismos parámetros técnico-económicos utilizados en las metodologías tradicionales antes desarrolladas. A continuación, se presenta una tabla con los valores utilizados para la valorización de rajo y los estratos subterráneos. Parámetros de evaluación de los estratos tipo rajo. Ritmo de 230 ktpd producción Precio 1.5 US$/lb Costo de Proceso
5.5
US$/t
Costo Mina
1.8
US$/t
Recuperación
85
%
Costo de Venta
0.32
US$/lb
Tabla 60. Parámetros de evaluación de los estratos para el ILP ex traíbles vía rajo. Caso Andina Sur-Sur.
Parámetros de evaluación de los estratos tipo subterráneo. Ritmo de 23 ktpd producción Precio 1.5 US$/lb Costo de Proceso
5.5
US$/t
Costo Mina
5
US$/t
Recuperación
85
%
Costo de Venta
0.32
US$/lb
Tabla 61. Parámetros de evaluación de los estratos para el ILP ex traíbles vía subterránea. Caso Andina Sur-Sur.
Costo de inversiones para la subterránea. Costo de desarrollo
520
k$/dpt
Tramming
5
kUS$/m
Costo del pique 8 kUS$/m Tabla 62. Costo de inversiones para la subterránea. Caso Andina Sur -Sur.
66
A continuación, se presenta una tabla que detalla la valorización de los estratos que son minados por el rajo y por la mina subterránea, respectivamente. Lc (%Cu) Estrato Mineral (t) Estéril (t)
Cu Potencia del (%Cu) estrato (m)
0.6
Rho1 (año)
P1 (MUS$)
1 0 345 0.00 60 1 -621,261.3417 2 1 4,076 0.63 60 1 -7,329,040.782 3 205 23,060 0.79 60 1 -39,433,603.78 4 11,038 60,738 1.13 60 1 86,940,091.18 5 29,954 88,088 1.26 60 1 457,827,630 6 34,425 124,313 1.19 60 1 433,541,278.6 7 42,059 162,518 1.07 60 1 392,075,021.9 8 55,161 197,268 0.94 60 1 391,826,696.6 9 68,088 227,795 0.90 60 1 446,210,662.6 10 78,382 256,331 0.87 60 1 474,996,202.3 11 84,441 293,522 0.88 60 2 493,802,141 12 84,675 340,617 0.88 60 2 418,751,229.6 13 89,666 388,804 0.88 60 2 384,983,588.8 14 95,520 450,462 0.88 60 2 357,058,373.8 15 101,237 514,293 0.87 60 2 275,162,702.8 16 99,741 590,600 0.84 60 2 65,713,199.68 17 103,830 661,282 0.84 60 2 -27,107,856.7 18 121,364 722,219 0.84 60 2 65,208,105.52 19 133,639 790,884 0.84 60 2 75,951,971.48 20 139,853 868,023 0.84 60 2 23,751,813.08 21 145,780 950,103 0.83 60 2 -114,209,738.2 22 149,242 1,031,762 0.80 60 2 -305,406,265.2 23 149,240 1,132,702 0.78 60 2 -554,432,978.4 24 148,075 1,230,798 0.79 60 2 -720,425,223 25 145,310 1,334,875 0.82 60 2 -839,185,038.2 26 151,718 1,431,920 0.84 60 2 -855,460,172.2 27 161,655 1,535,359 0.87 60 2 -843,920,297.4 28 166,812 1,644,541 0.90 60 2 -850,014,430.3 29 179,992 1,752,751 0.92 60 3 -803,554,895 30 184,882 1,869,219 0.94 60 3 -873,700,858 Tabla 63. Valorización de los estratos para el ILP extraíbles vía rajo. Caso Andina Sur-Sur.
67
Cil_Lc A=0.6 Estrato 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Costo de Costo de Costo de Costo de Rho 2 desarrollo sostenimiento tramming transición (año) (MUS$) (MUS$) (MUS$) (MUS$) 89 0.99 547,874 1,096 7 4 1,107 11 84 0.98 533,866 1,068 8 4 1,080 11 74 0.96 507,453 1,015 10 4 1,028 9 65 0.94 480,671 961 11 3 976 8 58 0.91 456,438 913 12 3 929 7 53 0.88 435,021 870 13 3 887 7 56 0.84 422,137 844 15 3 862 7 62 0.82 417,446 835 16 3 854 8 69 0.85 418,696 838 17 3 858 9 74 0.88 422,329 845 18 3 866 9 75 0.89 425,938 852 19 3 875 9 73 0.89 428,080 856 21 3 880 9 70 0.89 428,424 857 22 3 882 9 65 0.87 426,425 853 23 3 880 8 64 0.86 424,178 849 24 3 876 8 65 0.87 422,839 846 25 3 875 8 58 0.87 418,965 838 27 3 868 7 47 0.87 411,921 824 28 3 855 6 39 0.90 402,863 806 29 3 838 5 37 0.92 394,157 788 30 3 822 5 39 0.88 386,756 774 31 3 808 5 35 0.91 378,863 758 33 3 794 5 31 0.95 370,719 742 34 3 778 4 27 1.13 362,281 725 35 3 763 4 26 1.23 354,110 708 36 3 747 4 24 1.22 346,240 693 37 3 733 3 10 1.13 335,628 671 39 3 713 2 0 0.64 323,649 647 40 3 690 1 0 0.77 312,496 625 41 3 669 1 0 0.00 302,080 604 42 3 649 0 Tabla 64. Valorización de los estratos para el ILP extraíbles vía subterránea. Caso Andina Sur- Sur. Mineral (Mt)
Ley de Cu (%)
Footprint (m2)
P2 (MUS$) 1,019 940 782 665 563 479 446 474 571 653 694 675 643 566 542 567 506 418 368 363 348 333 329 396 425 399 141 0 0 0
Finalmente, se correrá el modelo de optimización para el Sur Sur siguiendo aquella metodología que obtuvo el mayor valor de VAN en la prueba en el caso de estudio que corresponde a la iii y es aquella que utiliza una curva beneficio – ley y deja la ley de corte variable para cada estrato.
68
iii) Resultados del modelo de optimización para Andina Sur Sur.
Utilizando la evaluación de los estratos que serán extraídos por el rajo y por la subterránea y con una restricción de altura de columna máxima de 360 m., se obtuvieron los siguientes resultados con respecto al VAN para la corrida del modelo de optimización en el caso de Andina Sur-Sur. MAX PIT
MED PIT
MIN PIT
H=2
T.S. (s)
VAN MUS$
H=2
T.S. (s)
VAN MUS$
H=2
T.S. (s)
VAN MUS$
A
0
0
A
0
0
A
634
389
H=3
T.S. (s)
VAN MUS$
H=3
T.S. (s)
VAN MUS$
H=3
T.S. (s)
VAN MUS$
A
0
0
A
0
0
A
1,392
412
H=4
T.S. (s)
VAN MUS$
H=4
T.S. (s)
VAN MUS$
H=4
T.S. (s)
VAN MUS$
A
0
0
A
0
0
A
2,783
437
H=5
T.S. (s)
VAN MUS$
H=5
T.S. (s)
VAN MUS$
H=5
T.S. (s)
VAN MUS$
A
0
0
A
0
0
A
1,031
457
H=6
T.S. (s)
VAN MUS$
H=6
T.S. (s)
VAN MUS$
H=6
T.S. (s)
VAN MUS$
A
0
0
A
0
0
A
2,780
472
H=S.R.
T.S. (s)
VAN MUS$
H=S.R.
T.S. (s)
VAN MUS$
H=S.R.
T.S. (s)
VAN MUS$
A
0
0
A
0
0
A
3,411
517
Tabla 65. Resultados del VAN para ley de corte variable y la discretización de la subterránea tipo A. Caso Andina Sur-Sur.
Por otro lado, con respecto a la geometría que se obtiene para el pit final y para el nivel donde se ubica el nivel de hundimiento, se obtienen éstos resultados.
H=2 A H=3 A H=4 A H=5 A H=6 A H=S.R. A
MAX PIT Fondo del Cota del Pit [m] footprint [m]
HOD real
Fondo del Pit [m] Fondo del Pit [m] Fondo del Pit [m] Fondo del Pit [m] -
Cota del footprint [m] Cota del footprint [m] Cota del footprint [m] Cota del footprint [m] -
HOD real HOD real HOD real HOD real -
Fondo del Pit [m] -
Cota del footprint [m] -
HOD real -
H=2 A H=3 A H=4 A H=5 A H=6 A H=S.R. A
MED PIT Fondo del Cota del Pit [m] footprint [m]
HOD real
Fondo del Pit [m] Fondo del Pit [m] Fondo del Pit [m] Fondo del Pit [m] -
Cota del footprint [m] Cota del footprint [m] Cota del footprint [m] Cota del footprint [m] -
HOD real HOD real HOD real HOD real -
Fondo del Pit [m] -
Cota del footprint [m] -
HOD real -
H=2 A H=3 A H=4 A H=5 A H=6 A H=S.R. A
MIN PIT Fondo del Cota del Pit [m] footprint [m]
HOD real
3,304 Fondo del Pit [m] 3,304 Fondo del Pit [m] 3,304 Fondo del Pit [m] 3,304 Fondo del Pit [m] 3,304
Cota del footprint [m] 2,764 Cota del footprint [m] 2,764 Cota del footprint [m] 2,764 Cota del footprint [m] 2,764
HOD real 3 HOD real 4 HOD real 5 HOD real 6
Fondo del Pit [m] 3,304
Cota del footprint [m] 2,704
HOD real 10
Tabla 66. Resultados de la geometría de transición para ley de corte variable y la discretización de la subterránea tipo A. Caso Andina Sur-Sur.
69
Finalmente, la política de leyes de corte que permite alcanzar cada uno de los VAN óptimos para las combinaciones de discretización de pit y altura de columna, son las siguientes. MIN PIT H=2 LC A %Cu H=3 LC A %Cu H=4 LC A %Cu H=5 LC A %Cu H=6 LC A %Cu H=FREE LC A %Cu 1
0.2
1
0.2
1
0.2
1
0.2
1
0.2
1
0.2
2
0.3
2
0.3
2
0.3
2
0.3
2
0.3
2
0.3
3
0.2
3
0.2
3
0.2
3
0.2
3
0.2
3
0.2
4
0.2
4
0.2
4
0.2
4
0.2
4
0.2
4
0.2
5
0.2
5
0.2
5
0.2
5
0.2
5
0.2
5
0.2
6
0.2
6
0.2
6
0.2
6
0.2
6
0.2
6
0.2
7
0.2
7
0.2
7
0.2
7
0.2
7
0.2
7
0.2
8
0.4
8
0.4
8
0.4
8
0.4
8
0.4
8
0.4
9
0.4
9
0.4
9
0.4
9
0.4
9
0.4
9
0.4
10
0.5
10
0.5
10
0.5
10
0.5
10
0.5
10
0.5
11
0.6
11
0.6
11
0.6
11
0.6
11
0.6
11
0.6
12
0.6
12
0.6
12
0.6
12
0.6
12
0.6
12
0.6
13
0.7
13
0.7
13
0.7
13
0.7
13
0.7
13
0.7
14
0.7
14
0.7
14
0.7
14
0.7
14
0.7
14
0.7
15
0.7
15
0.7
15
0.7
15
0.7
15
0.7
15
0.7
16
16
16
16
16
16
0.9
17
17
17
17
17
17
0.9
18
18
18
18
18
18
0.9
19
19
19
19
19
0.9
19
0.9
20
20
20
20
0.8
20
0.9
20
0.9
21
21
21
0.6
21
0.8
21
0.9
21
0.9
22
22
0.6
22
0.6
22
0.8
22
0.9
22
0.9
23
23
0.6
23
0.6
23
0.8
23
0.9
23
0.9
24
24
0.6
24
0.6
24
0.8
24
0.9
24
0.9
25
25
25
25
25
25
0.9
26
26
26
26
26
26
27
27
27
27
27
27
28
28
28
28
28
28
29
29
29
29
29
29
30
30 30 30 30 30 Tabla 67. Política de leyes de corte para la extracción de los estratos. Discretización MIN PIT.
70
7.4. Comparación entre las alternativas de transición para Andina Sur-Sur. En las siguientes tablas se encuentran los resultados obtenidos para las corridas realizadas con el modelo de optimización. Para cada una de las corridas se obtuvo un VAN y una geometría asociada, las cuales se compararán más adelante con los resultados obtenidos en la metodología secuencial y con la incremental por fases. De la misma manera con lo sucedido en Palabora, el procedimiento de comparación de los resultados del modelo se realizará valorizando el resultado final entregado por el modelo con la información obtenida en la metodología incremental por fases, considerando además la ley de corte asociada a cada unidad de reserva. Considerando que se utilizará el criterio del mayor VAN para elegir la metodología que genera el proyecto del complejo minero, se presentan a continuación los valores alcanzados por cada una. El valor asociado al rajo y al block caving se encuentran actualizados al inicio del proyecto. Valor del proyecto
VAN rajo [MUS$] VAN Block Caving[MUS$] VAN del proyecto [MUS$]
Secuencial 1,986 5,6 1,992 Incremental por fases 1,954 21 1,975 ILP 2,181 132 2,313 Tabla 68. Comparación entre metodologías según VAN. Caso Andina S ur-Sur.
Además, la geometría obtenida para las diferentes metodologías de transición se encuentran resumidas en la siguiente tabla. Diseño Cota fondo fondo del rajo Cota del del Footprint de la transición Secuencial 2,971 2,524 Incremental por fases 3,308 3,036 ILP 3,304 2,764 Tabla 69. Comparación entre metodologías según el diseño de la transición. Caso Andina Sur-Sur.
Si se examinan las cotas que definen el diseño de la transición por las tres metodologías, se observa que existe una cota sobre la cual los tres métodos determinan que el único método posible de extracción es el rajo, en éste caso corresponde a la cota 3,304. De ésta manera, existe un valor conformado por la valorización de la subterránea y de la fase marginal por debajo de la cota 3,304, que determina la potencial ganancia o perdida de mantener el método rajo por debajo de esta cota. Este valor puede ser denominado como un valor en competencia entre la última fase del rajo y la subterránea. Finalmente, el valor en competencia se encontrará descontado sólo por el tiempo que requiera extraer la última fase del rajo más el de la subterránea. Valor económico Últimas fases del rajo Extracción rajo Valor del BC Extracción Total en competencia [MUS$] [años] [MUS$] BC [años] [MUS$] Secuencial 124 36 41 89 165 Incremental por 0 0 958 26 958 fases ILP 0 0 936 13 936 Tabla 70. Valor económico en competencia entre el rajo y la subterránea. Caso Andina Sur-Sur.
71
Con el análisis del valor económico en competencia, se aprecia claramente que la mejor opción no corresponde a alcanzar el máximo individual de cada uno de los métodos de explotación, ya que éste valor en competencia es el menor para la metodología secuencial. Por otro lado, las alternativas que contemplan una simultaneidad de la evaluación de los métodos, obtienen mejores resultados. Este análisis deja en evidencia que la opción subterránea para el resultado de la metodología incremental por fases e ILP son comparables. Sin embargo, la ganancia en valor del proyecto, que permite definir como opción óptima al resultado del modelo de optimización, es la posibilidad de generar una política de leyes de corte para el rajo que aumenta el valor de la opción rajo y disminuye la vida útil del rajo. Lo cual hace que la valorización de la subterránea sea afectada de menor manera por la tasa de descuento.
72
8. Conclusiones y Recomendaciones. Para el caso de Palabora, se encontraron los siguientes resultados para cada una de las metodologías.
Valor del proyecto VAN rajo [MUS$] VAN Block Caving[MUS$] VAN del proyecto [MUS$] Secuencial 3,593 27 3,620 Incremental por fases 3,580 72 3,653 ILP 3,754 124 3,878 Tabla 71. Tabla comparación de resultados para la transición. Caso Palabora.
Cada una de las metodologías estudiadas a lo largo de este trabajo, han permitido enfrentar de manera eficaz y detallada el problema de realizar una transición de una minería a rajo abierto a una mina subterránea de block caving. Cada una ha entregado una geometría y un VAN asociado a su desarrollo. Con respecto a las metodologías tradicionales, si bien dan una buena definición del proyecto, el tiempo que requiere realizar el estudio y, además, generar sensibilidad de aquellos factores importantes en la definición de la geometría de explotación, toma más tiempo en comparación con el utilizado para cambiar los parámetros de entrada que son requeridos para correr el modelo de optimización. Si se considera el resultado final obtenido a través de los dos tipos de metodologías aplicados en el presente trabajo, en el procedimiento tradicional se obtiene un pit final y una determinada elevación donde se debe ubicar el nivel de hundimiento del block caving. Sin embargo, la definición de cada uno de éstos está sujeta a la información disponible al instante de la evaluación de cada uno de los métodos de explotación, es decir, la optimización del rajo se hará teniendo en cuenta todos los recursos minerales, para posteriormente, evaluar la subterránea considerando sólo los recursos que no fueron tomados en cuenta dentro del pit final. Por otro lado, el utilizar el modelo de optimización permite evaluar simultáneamente si el siguiente estrato de evaluación debe ser extraído por el método subterráneo o por el rajo, pero, además de definir el estrato óptimo en el cual se realiza el cambio del método de explotación, es posible generar una política de leyes de corte que mejora el VAN del proyecto, como se ve representado en los resultados de la metodología correspondiente.
73
Si se comparan los resultados correspondientes a la geometría de la transición obtenidos con las metodologías tradicionales y con el modelo de optimización, se obtiene un resultado interesante en el caso de Palabora, esto es, que existe una reducción del tamaño del rajo, es decir, una cota del fondo del rajo no tan profunda. Ello implica que el rajo definido por la metodología secuencial extraía estratos que generan menos valor que los que serían obtenidos al realizar la transición del método de explotación. Si se valoriza el resultado entregado por el modelo de optimización con los valores de los estratos calculados por la metodología incremental por fases se obtiene una considerable mejora mejora en el VAN del proyecto con respecto a cualquiera de las metodologías tradicionales utilizadas a priori. De ésta manera se demuestra que el resultado del modelo de optimización mejora notablemente el diseño y, por ende, el valor económico de la transición rajo subterránea. En el caso de estudio de Palabora, se realizó un análisis de sensibilidad para diferentes factores que podrían de una u otra manera incidir en la valorización o en el tiempo de extracción de los estratos, y por ende, afectar a los resultados entregados por el modelo de optimización. Con respecto a estos, se tiene que una variación en la discretización del rajo puede llegar a alcanzar una diferencia del 25% del valor económico si se utiliza un método de elección fijo y predefinido de la ley de corte. Por otro lado, se obtiene que al variar la altura de columna se genera una diferencia máxima del 7%. Ahora bien, con respecto a la geometría de la transición, también es afectada por el efecto de la agregación del pit y por la altura de columna. Con respecto a la primera de éstas, el fondo de pit puede llegar a cambiar hasta en 840 metros y además definir la existencia o no de una mina subterránea bajo el rajo, según la agregación subterránea y el método de definición de la ley de corte. Con respecto a los cambios en la altura de columna de extracción, se alcanza una diferencia de hasta 180 m en el fondo del pit y nuevamente la existencia o no de la mina subterránea como continuación del pit. Considerando el precio de evaluación del cobre, se obtiene que el valor del proyecto se comporta de manera directamente proporcional con respecto al precio, es decir, que al aumentar el precio, el valor del proyecto aumenta. A modo de ejemplo, si se aumenta el precio en un 20%, el valor del proyecto aumentará en un 69%. Por su parte, en el análisis de sensibilidad con respecto al ritmo de producción del rajo y de la mina subterránea, se tiene que los resultados son más sensibles a un cambio en el ritmo de producción del rajo y más insensible con respecto a la producción del block caving.
74
Finalmente, el análisis en relación con el valor económico en competencia, permite definir si es realmente rentable extraer las últimas fases del rajo. En el caso de Palabora, se observa claramente que éstas últimas fases no permitían alcanzar el mejor valor asociado a la explotación subterránea, y por ende, el del proyecto global. Adicionalmente es posible determinar si la transición es realmente rentable, dado que se determina el valor de la subterránea al momento de generar la decisión de realizar el cambio en el método de explotación. De esta manera, el monto de la ganancia que está en juego en la decisión de la transición llega a alcanzar los MUS$ 500 en el peor caso para Palabora. Valor económico en Últimas fases del rajo Extracción rajo Valor del BC Extracción competencia [MUS$] [años] [MUS$] BC [años] Secuencial 85 11 215 23 Incremental por fases 34 4 648 32 ILP 0 0 839 14 Tabla 72. Valor económico en competencia para el caso de Palabora.
Total [MUS$] 300 682 839
Por otro lado, la comparación de las tres metodologías aplicadas para el caso del Sur Sur se encuentra en la siguiente tabla. Valor del proyecto
VAN rajo [MUS$] VAN Block Caving[MUS$] VAN del proyecto [MUS$]
Secuencial 1,986 5,6 1,992 Incremental por fases 1,954 21 1,975 ILP 2,181 132 2,313 Tabla 73. Tabla comparación de resultados para la transición. Caso Sur Sur.
Al aplicar la valorización de los estratos obtenidos en la metodología de incremental por fases en los resultados encontrados por el modelo de optimización y considerando la política de leyes de corte, se encuentra que mejora el valor del proyecto de transición. Principalmente por el hecho de optimizar el rajo a través de la política de leyes de corte y de generar una vida útil del rajo menor, lo cual hace que el valor de la subterránea no se vea tan afectada con el efecto de la tasa de descuento, como sucede en las otras dos metodologías estudiadas. En relación a la geometría propuesta para la transición, la metodología incremental por fases coincide en la definición del fondo del pit final del modelo, sin embargo, difiere en la cota de la ubicación del footprint. La principal razón de ésta diferencia se encuentra en el análisis del valor económico en competencia entre las últimas fases del rajo Sur-Sur más la subterránea entre las alternativas. Si bien, el valor de las subterráneas puede ser comparable en órdenes de magnitud, la decisión de la metodología incremental por fases no considera la inversión en infraestructura de la mina subterránea, el que si es considerado por el modelo de optimización. Además, se desprende del análisis del valor económico en competencia para el Sur-Sur, que el monto en ganancia que conlleva la decisión de realizar el proyecto de la transición podría llegar hasta los MUS$ 800.
75
Por lo tanto, dada la diferencia en geometría entre los tipos de mineralización de los modelos de bloques que se utilizaron para probar las diferentes metodologías, se puede aseverar que el modelo de optimización es ampliamente aplicable y que permite generar un diseño alternativo a las que pueden ser encontradas aplicando las herramientas de planificación convencionales y que a la luz de una valorización minera genera un proyecto rentable. Uno de los principales problemas a los que se ven enfrentadas las operaciones a rajo abierto en grandes yacimientos corresponde a la planificación de fases cada vez más grandes a medida que el mineral se va profundizando, es decir, presentan un aumento drástico de la razón estéril mineral. Por lo tanto, se presentan problemas al momento de la estimación de flota de transporte, estabilidad geotécnica de las paredes del rajo, disposición de botaderos, entre otros y consecuentemente una disminución de los beneficios asociados a la explotación de estas últimas fases incluidas dentro del pit final definido como óptimo por Whittle. Sin embargo, dados los resultados obtenidos a partir del modelo de optimización, se concluye que los beneficios generados por las labores subterráneas logran competir con los generados por los incrementos marginales del rajo, es más, son mayores y permiten mejorar el VAN del proyecto. En consecuencia, el modelo de optimización entrega al planificador una herramienta de análisis que permite sensibilizar el resultado de la transición tomando en cuenta la variabilidad de los parámetros económicos y de diseño, definiendo un proyecto minero en que los beneficios producidos por la extracción del mineral son óptimos y con el método más adecuado.
76
9. Bibliografía. 1. ARANCIBIA E., FLORES G. Design for Underground Mining at Chuquicamata orebody. Scoping Engineering Stage; Proceedings MassMin 2004; 603-609; 2004; Santiago Chile. 2. AVENDAÑO A. Límite explotación rajo abierto v/s minería subterránea; Memoria de Ingeniería; Universidad de Chile; 2005; Santiago Chile. 3. AZCARATEGUI B. Modelo de optimización para la transición de minería a rajo abierto a minería subterránea; Tesis de Magister; Universidad de Chile; 2005; Santiago Chile. 4. BAKHTAVAR E., SHAHRIAR K. and ORAEE K. A model for determining optimal transition depth over from open-pit to underground mining; 393-400; Proceedings Massmin; 2008; Luleå, Sweden. 5. BRANNON C., CASTEN T., JOHNSON M. Design of the Grasberg block cave mine; Proceedings MassMin 2004; 623-628; 2004; Santiago Chile. 6. BRUMMER et Al. The Transition from open pit to underground mining: an unusual slope failure mechanism at Palabora; Proceedings of The SAIMM International Symposium on Stability of Rock Slopes in Open Pit Mining and Civil Engineering; 2006. 7. CÁCERES J. “Open Pit to Underground Transition, How much value has been left behind? ”, Proceedings Minin 2006, Santiago, Chile, 2006. 8. C. L. WHARTON. “What they don’t teach you in mining school – Tips and Tricks with Pit Optimizers”; Surface Mining 1996. Johannesburg, South African Institute of Mining and
Metallurgy. 1996, pp. 17-22. 9. DIERING T. “PC-BC: A block cave design and draw control system ”. MassMin, (pp. 469484). Brisbane. Australia, 2000. 10. FUENTES S. “Going to an underground (UG) mining method”; Proceedings MassMin 2004; 633-635; 2004; Santiago Chile. 11. FUENTES S., CÁCERES J. “Block/panel caving pressing final open pit limit”; Technical paper, CIM Bulletin, Vol 97, N°1082, 2004. 12. GOLDBERG A. V. “Efficient graph algorithms for sequential and parallel computers”, PhD thesis, MIT, 1987. 13. GOLDBERG A. V. and TARJAN R. E.; “A new approach to the maximum flow problem”, J. Assoc. Comput. Mach.; 921-940; 1988. 77
14. HOCHBAUM. D, CHEN. A. “Performance Analysis and Best Implementations of Old and New Algorithms for the Open-Pit Mining Problem, INFORMS, Vol 48, No 6, pp. 894-914, Noviembre - Diciembre 2000. 15. KOVOROV S. “Method for determining optimal open pit limits”; Rapport Technique ED 74R-4; Dept. of Mineral Engineering; Ecole Polytechnique de Montreal; Canada; 1974. 16. LERCHS. H, I. F GROSSMAN. “Optimum design of Open-pit Mines”, Transactions of the Canadian Institute of Mining and Metallurgy, Volume LXVII, pp. 17-24, 1965. 17. MAGÑIN J.M. “Reformulación de una planificación Minera de Largo Plazo en una Mina a Rajo Abierto”, Memoria de Ingeniería, Universidad de Chile, 2002. 18. NEWMAN A., RUBIO E. Y YANO C. Optimizing the Transition from Surface to Underground Mining; Presentation in Workshop on Operations research in Mining, 2008, Viña del Mar. 19. OLAVARRÍA S., ADRIASOLA P., KARZULOVIC A. Transition from open pit to underground mining at Chuquicamata, Antofagasta, Chile; Proceedings of The SAIMM International Symposium on Stability of Rock Slopes in Open Pit Mining and Civil Engineering; 2006. 20. PAIS G. Apuntes del Laboratorio de Planificación Minera del Diplomado de Planificación Minera; 2008; Universidad de Chile; Santiago, Chile. 21. PANA M. T. The simulation approach to open pit design. Preceedings 5 th APCOM Symposium, Tucson Arizona, 1965. 22. RUBIO E. Apuntes del curso de diseño de minas subterráneas; 2007; Universidad de Chile; Santiago Chile. 23. RUBIO E. Apuntes del Curso “Tópicos Avanzados en Planificación Minera”, Universidad de Chile, 2008. 24. TOBIE R., JULIN D., “Block caving. General description”, Underground mining method handbook, 967-972, New York, 1982.
78
10. Anexos. 10.1. Flujos de caja descontados para los estratos extraídos por rajo y subterránea. Caso Palabora.
79
10.2. Flujos de caja descontados para los estratos extraídos por rajo y subterránea. Caso Sur Sur.
80
10.3.
Valorización de estratos ILP, extraídos vía rajo. Caso Palabora.
10.3.1. Max Pit.
81
82
83
84
85
10.3.2. Med Pit.
86
87
88
89
90
10.3.3. Min Pit.
91
92
93
94
95
10.4.
Valorización de estratos ILP, extraídos vía subterránea. Caso Palabora.
10.4.1. Discretización A.
96
97
98
10.4.2. Discretización B.
99
100
101
102
10.5.
Valorización de estratos ILP, extraídos vía rajo. Caso Andina Sur-Sur.
103
104
105
106
107
10.6. Valorización de Caso Andina Sur-Sur.
estratos
108
ILP,
extraídos
vía
subterránea.
109
110
111
112
10.7.
Paper presentado en APCOM 2009, Vancouver, Canadá.
Open Pit Underground Transition Supported by Integer Linear Programming Modeling.
Andrés Solar D. Master of Mining Engineering Candidate. Universidad de Chile.
Enrique Rubio E. Assistant Professor Mining Engineering Department Universidad de Chile.
Alexandra Newman. Associate Professor Business Department Colorado School of Mines.
Abstract. Several large open pit mines are facing the challenge to design and plan an underground transition to sustain the long term future of their operations. The desired underground method is Block Cave since its natural high productivity and low operating cost that could compete against the open pit cost to mine the final open pit pushbacks. It has been found at several operations where this transition is taking place that the pit limit at which transition is performed is rather deep. This aims to think that perhaps several transition projects have started late in the life of an operation producing an over stress in the transition project and reducing the flexibility in the design of the overall infrastructure and a substantial lost in the complete project NPV. Traditionally the approach to plan open pit underground transition consists of sequentially deciding the final pit using a Lerchs and Grossman type of optimization to find the overall maximum pit shell and then running block cave simulations to find the best footprint for the underground option. The research summarized in this paper is related to the development of an integer linear programming model that aims to compute simultaneously the maximum NPV of the combined option final limit open pit and block cave underground footprint. The model has been tested using two large block models representing two of the current operations in the process of transiting from open pit to underground operations. Finally, the integrated open pit underground transition model is going to be compared to the traditional sequential approach in order to define which method achieves the best NPV. The paper illustrates both case studies in depth.
113
Introduction. One of the important problems that large open pits are facing in nowadays is doesn’t have block model information enough. This makes that the mine planning support one kind of mining method, that not necessary let to achieve the optimum NPV of the project. So, if at the beginning of mine’s life it makes a study that consider a transition of the mining method, from open pit to underground, produce a new scenario where the NPV and the life of mine raises. The underground method that fit appropriately in transition are the caving methods, because these ones have lowest operation costs and achieve the highest production rates. Nowadays haven’t plenty knowledge of how aboard the transition problem. The actual planning and design techniques achieve the maximum of each operations, but it knows that the maximum of a system isn’t necessary achieve optimizing each part of the system.
But, to test that an ore body could support a mine transition, it must satisfy that one slice of rock achieve the same value if it’s evaluate with the open pit or underground method. It could formulate the following mathematical equations to probe the mine transition in a generic ore body. k 1
v ok
ρ
o t k
u t k
v uk
ρ
v k
o
RF l k
o
1
v uk
RF l uk
C um
o, u
v k ρ
z 1
e
Transition condition o
Cm
Cp
Cp
Value of stratum k mined by openpit/underground.
1
u t o, k
u t o, k
(1 r)
r
u
tz
Discount factor considering openpit/underground mine time for slice k.
Discount rate (10%).
RF
(Price – Selling Cost)*Recovery*22,04 = Revenue Factor.
u l o, k
e
Cut off grade for stratum k mined by openpit/underground. waste/ore rate.
o ,u
C m
Cp
Mining cost. (openpit/underground) Processing Cost.
114
The marginal condition what defines the transition its true, if the variable takes these values; Table 1. Variables value what defines the condition true. Mine Cost OP Production Rate OP Price [US/lb] Gradek OP [%] [US$/t] [tpd] 1.95 1.02 150.000 0.31 Mine Cost UG Recovery [%] Waste/Ore Rate Gradek UG [%] [US$/t] 80 5 3.01 0.61 Selling Cost Processing Cost Production Rate UG [US$/lb] [US$/t] [tpd] 0.31 5 35.000 RF [US$/t/%] 28.97
Due to the values that define true the condition it could be found in a real ore body, then exist a highest possibility that a transition slice be defined in a real mine. In this way, is totally valid to study this problem and try to find an intelligent way to find the slice where occurred the transition. Block model Information. To take place this investigation, it used a real mine block model (pipe form) with a plane topography at level 296 and whose have this characteristics; Block Size Number of blocks Cu Grade (%Cu) Average density
X Y 40 20 71 149 Max Min 1.6 0 2.65[t/m ]
115
Z 60 37
Technical parameters and assumptions. It assumed a 5% of mining dilution fraction and 95% of recovery mining fraction. By the other hand isn’t applied a cost adjustment factor (CAF) or processing cost adjustment factor (PCAF). It considered a 45° constant pit slope, a 10% of discount rate. To evaluate the open pit slices it uses the following economical assumptions. Table 2. OP Economical Assumptions. OP Economical Assumptions 2 [US$/lb] Price 5 [US$/t] Processing Cost 1 [US$/t] Mining Cost 90 % Recovery 0.26 [US$/lb] Selling Cost
On the other hand, underground production rate is set in 23,500 tpd and a development cost of 2000 dollars per square meter. To evaluate the underground slices it uses the following economical assumptions. Table 3. UG Economical Assumptions. UG Economical Assumptions 2 US$/lb Price 5 US$/t Processing Cost 5 US$/t Mining Cost 90 % Recovery 0.26 US$/lb Selling Cost
Methodology. Using the mine block model, it presents two traditional approaches what use the actual planning and design mine tools. Open pit planning used Whittle software and underground planning and design used GEMS PCBC Footprint finder software.
The first traditional approach is a sequential mine plan and design, in other words, a long term planning and design for underground and open pit mine, considering that the open pit take place first. The second one present use incremental rock volumes, what in the case of the open pit look like a inverted concentric cones what are defined with the L&G algorithm what be evaluate considering the OP economical assumptions, and the underground increments looks like a concentric cylinder of different radius what are evaluates considering the UG economical assumptions. To make easier the comparison between the slices the incremental volumes are cut by levels, leaving a 60 meter width single slice. These increments make possible to build incremental discounted cash flow curves for open pit and underground mine (illustration 1). Then, the transition slice will be the one where the incremental discounted cash flow curves cut what define that the revenues goes down for apply an incompetent method.
116
Change final pit
OP NPV Change footprint elevation
Incremental OP discounted cash flow
Incremental UG discounted cash flow.
Ilustration 1: Graphical method to find the transition.
The last methodology presents to avoid the open pit underground transition, its use an integer programming model. This model use binary variables which defines the mine method for each rock slice. This decision considers the slice OP and UG economic value simultaneously. The data what feed the model take count the geometry of the increments for the OP and UG, an invert concentrically cone and a cylinder. These volumes could represent the biggest and smallest pit or footprint area. Finally, the goal function is maximizing the project NPV.
Results. Traditional approach #1 results.
According to the methodology, an 180,000 tpd open pit production rate achieves the maximum open pit NPV and an -1504 m elevation footprint maximizes the block caving NPV. Table 4. Methodology 1 results. Final Pit N° 45 Pit Bottom VAN Ore Waste LOM [m] [MUS$] [Mt] [Mt] [yr] -1147 3,593 2,146 8,604 33 Best FP FP Elevation VAN Ore Cu LOM [m] [MUS$] [Mt] Grade % [yr] -1504 603 197 0.76 23 Entire Project VAN [MUS$] 3,620
117
Traditional approach #2 results.
Evaluating the rock slices defined using the toe block elevations for OP/UG method, it found the follow results. Table
5.
Methodology
2
results.
Final Pit Pit Bottom [m] -964
VAN Waste Ore [Mt] [MUS$] [Mt] 3,580 1,757 5,587 Best FP FP Elevation VAN Cu Ore [Mt] [m] [MUS$] Grade % -1,264 72 764 0.67 Entire Project VAN [MUS$]
LOM [yr] 27 LOM [yr] 32
3,653
ILP methodology results.
Assuming a height of draw for the block caving of 300 m, it obtained the follow results. The model was implemented in AMPL and solved using CPLEX Parallel 10.2 (8 threads).
OP Shape
Big pit
Medium Pit
Small Pit
Table 6. Methodology 3 results. Footprint Area [m2] Pit Bottom Footprint Elevation -844 -1,144 A (Biggest) -964 -1,264 B -964 -1,204 C -964 -1,264 D -964 -1,264 E -964 -1,444 F -904 -1,444 G H (Smallest) -904 -1,204 -964 -1,264 A (Biggest) -964 -1,264 B -964 -1,264 C -964 -1,264 D -1,024 -1,324 E -1,024 -1,264 F -904 -1,204 G H (Smallest) -1,024 -1,324 -1,024 -1,384 A (Biggest) -1,024 -1,324 B -904 -1,264 C -964 -1,264 D -904 -1,204 E -1,024 -1,324 F -1,024 -1,264 G H (Smallest) -1,024 -1,324
118
NPV [MUS$] 5,979 6,011 5,976 6,051 6,040 6,040 6,071 6,089 5,804 5,803 5,801 5,840 5,835 5,800 5,857 5,883 5,609 5,610 5,574 5,642 5,607 5,640 5,640 5,686
Results Analysis. The comparison criterion is the NPV of the entire project. Below its present the open pit and underground PV discounted to the Project start time. Methodology Traditional Approach Graphical Approach ILP Approach
Table 6. Project NPV summary. OP NPV [MUS$] UG NPV [MUS$] 3,593 3,580 -
Project NPV [MUS$] 27 3,620 72 3,653 6,089
It easy to find out that exist a highest NPV when the project is studied using the ILP approach. This only confirms that the maximum of a system is not necessary achieve maximizing the parts separately. Furthermore, it must consider the transition project design obtained with the ILP and the traditional approaches. Exist big differences, approximately 240 m in the pit depth between the ILP and the worst case of the traditional approach. This depth decrease make the life of mine shorter and doesn’t present geotechnical problems in pit walls. Table 7. Project design summary. Design Comparison
Pit Bottom [m]
Traditional Approach
FP elevation [m]
-1,147
-1,504
Graphical Approach
-964
-1,264
ILP Approach
-904
-1,204
Conclusion and finals remarks. Each one of the methodologies presents have made possible to face up in an intelligent way the open pit underground transition. It obtained a design and a NPV what could be expected. In this investigation it found that the best NPV it achieve using the ILP methodology (MUS$ 6,089), this shows and probe that the maximum of a system is not necessary achieve for maximize each part. Another improvement obtained for use the third methodology is haven’t a too much deep pit, a shorter life of mine and an underground mine not too deep, the last improve makes to avoid the too deep problems like rock blasting. The ILP find the first transition slice and don’t care about the ore in deep, this support the idea to place another underground mine below the first one defined with the ILP. This problem could be called as multi-lifting problem, and the solution its part of a simple underground design with the ore left.
An important conclusion what find out from the ILP analysis is that the deepest benches of the optimum traditional pit are mined by the block caving. That make us think that the traditional methodology isn’t find the optimum transition final pit (final pit of the complex) and this make to lose value of that ore. Although, that benches report a greater than zero discounted cash flow, by underground method will be even greater.
119
