Descripción

Este es un ejemplo donde se compara la metodología tradicional de agendamiento mediante pits anidados y el agendamiento directo de bloques sobre un modelo de bloques específico. 

Agendamiento mediante pits anidados

El agendamiento mediante pits anidados corresponde a utilizar el algoritmo de Lerchs y Grossmann (u otro) sobre una serie de valorizaciones del modelo de bloques, parametrizadas por un multiplicador del precio o "revenue factor". A partir de los pits anidados generados de esta forma, se seleccionan pushbacks o "candidatos a fase", que sirven de guía para el diseño, y a continuación, son agendados para extracción en el tiempo.

Agendamiento directo de bloques

El agendamiento directo de bloques es un enfoque alternativo que pone énfasis en la temporalidad del problema y el costo de oportunidad, buscando resolver desde un principio el problema de en qué momento del tiempo es más conveniente extraer un bloque y cuál es la mejor decisión de procesamiento (incluyendo la opción de lastre). Tradicionalmente esta metodología se apoya en modelos de programación lineal de gran tamaño.

Resumen

En este experimento se compara el agendamiento de la producción de una mina a rajo abierto mediante dos metodologías: (i) utilizando pits anidados y a partir de ahí generando dos estrategias de extracción, una pit por pit de manera incremental, conocida en la literatura como "Best Case" y otra estrategia que extrae banco a banco, conocida como "Worst Case". Luego, (ii) mediante agendamiento directo de bloques, donde se plantea el problema utilizando modelos de programación lineal y se utilizan diversas estrategias para resolverlo, en particular, en este ejemplo mostraremos dos: utilizando una heuristica por ventanas, la cual llamaremos STW, y otra, utilizando una heuristica Toposort, la que llameremos TOPO. En ambas estrategias, se imponen restricciones sobre los tonelajes extraídos. No se hacen consideraciones de diseño.

En los cuatro casos, se analizan sus geometrías de extraccción, tonelajes de mineral y lastre, y valor descontado acumulado a lo largo de la vida del proyecto. En términos de valor, las estrategias de agendamiento directo de bloques muestran el mejor desempeño.

 

 

Estrategia     

Valor acumulado descontado

[MUSD]

 

PITS

ANIDADOS

 


WORST CASE

 

BEST
CASE

1.306



1.372

 

DIRECTO

BLOQUES

 

STW

 

TOPO

1.386

 

1.390

 

A continuación, se muestra el desarrollo del experimento, con todos los detalles de su implementación en DOPPLER.

 

Modelo de Bloques

Nombre

MDB RAJO 

Número de Bloques

151.898 

Dimensiones

5 x 5 x 5 

Rangos de Coordenadas 

X=214.961,125 - 215.251,125

Y=8.080.666,5 - 8.081.201,5

Z=247,5 - 422,5

Columnas en el Modelo

(aparte de coordenadas) 

zone

cu

au

ag

density

Comentarios

Se agregó columna de tonelaje 

VISUALIZACION DEL MODELO

En las siguientes imagenes, se muestran el modelo de bloques con: (i) ley de plata (ag) en cada bloque, (ii) ley de cobre (cu) y (iii) ley de oro (au). Finalmente se muestra (iv) una vista en dos secciones, filtrado por ley de oro.

bm ag

bm cu
(i) (ii)
bm au bm slices au
 (iii)  (iv)

 

CURVAS TONELAJE LEY

 A continuación, se presentan tres curvas tonelaje ley, una para cada tipo de metal: (i) plata (ag), (ii) oro (au) y (iii) cobre (cu) contenidos en el modelo de bloques.

CurvaTonAg ley CurvaTonAu ley CurvaTonCu ley
(i) (ii) (iii)

 

DISTRIBUCION DE LEYES

Para cada tipo de metal en el modelo de bloques: (i) plata (ag), (ii) oro (au) y (iii) cobre (cu), se muestran las distribuciones de ley.

HistogramaAg

HistogramaAu HistogramaCu
(i) (ii) (iii)

 

Pit Final

Valorización y Ángulo de Talud

Datos para valorización:

1) Precio del metal: Ag = 16 [USD/oz]; Au = 1.100 [USD/oz]; Cu = 2,0 [USD/lb] 

2) Costo de venta: Ag = 1,6 [USD/oz]; Au = 100 [USD/oz]; Cu = 0,3 [USD/lb] 

3) Costo de proceso: Ag = 8,0 [USD/Ton]; Au = 8,2 [USD/Ton]; Cu = 8,0 [USD/Ton]

4) Recuperación metalúrgica: Ag = 0,77; Au = 0,76; Cu = 0,85

5) Costo de Mina: 1,5 [USD/Ton]

Datos para ángulo de talud:

1) 48° (global)

Visualización de pit final

A continuación, se muestra el pit final obtenido a partir de los modelos de bloque y económico proporcionados, (i) vista isométrica y (ii) vista en sección.

 nofinal pit iso  nofinal pit slice
(i)  (ii) 

 

Pits Anidados

Valorización y Ángulo de Talud

Se utilizaron los mismos datos para valorización del caso anterior de pit final, pero utilizando una serie de ponderadores sobre el precio (revenue factors), donde para cada uno de ellos, se obtiene un pit anidado. 

La familia de revenue factors (RF) seleccionada va desde 0,00 hasta 1,00, utilizando un paso de 0,05 entre factores sucesivos. 

 

Gráfico Pit By Pit

centella edited

 

Visualización de pits anidados

En las siguientes figuras, se muestra la familia de pits anidados obtenida a partir de la elección de los revenue factors: (i) muestra una vista isométrica de los pits anidados dentro de la envolvente económica, mientras que (ii) muestra una vista en sección. Cabe destacar que el primer pit que "flota" es aquel que corresponde a un RF= 0,15.

 np bm  

np bm slice

(i)   (ii)

 

Best & Worst Cases

PLAN DE PRODUCCIÓN BEST CASE

Este caso corresponde cuando la extraccción se realiza siguiendo el orden de los pits anidados (pit por pit). Se impuso un límite en el tonelaje extraído de 400.000 [Ton], con el cual se obtuvo un horizonte de 18 periodos (anuales). No se dispone de la opción de restringir el mineral a planta en este módulo. Note que la producción alcanzó el máximo en varios periodos, dejando para el final del horizonte las mayores porciones de lastre: (i) muestra el plan de producción, y (ii) muestra una vista en sección del agendamiento, ambos utilizando una estrategia "best case", pit por pit.

NPV Best Case: 1.372 [MUSD]

PP BC bench pit
(i) (ii)

 

PLAN DE PRODUCCIÓN WORST CASE

Este caso corresponde cuando la extraccción se realiza siguiendo el orden banco a banco. Igual que antes, se impuso un límite en el tonelaje extraído de 400.000 [Ton] y se obtuvieron 18 periodos. No se dispone de la opción de restringir el mineral a planta en este módulo. Nótese que tanto la producción como el lastre en este caso se distribuyen de manera más uniforme a lo largo del horizonte del proyecto. Sin embargo, no se aprecian diferencias más allá de las mencionadas entre este caso y el anterior: (i) muestra el plan de producción, y (ii) muestra una vista en sección del agendamiento, ambos utilizando una estrategia "worst case", banco a banco.

NPV Worst Case: 1.306 [MUSD] 

PP WC bybench slice
(i) (ii)

 

Agendamiento Directo de Bloques

En esta parte, se usarán modelos de agendamiento directo de bloques. Para su resolución, mostraremos dos estrategias utilizadas: (i) heurística por ventanas (STW), la cual resuelve de manera incremental el problema, utilizando una ventana de periodos (menor al horizonte) y fijando parte, o completamente, esa ventana como solución, y (ii) heurística Toposort (TOPO), la cual explota la relajación continua del problema y define los tiempos esperados para realizar un agendamiento utilizando el orden dado por los tiempos esperados y las capacidades de producción.

Parámetros de Entrada

Número de periodos: 18

Tasa de descuento: 10%

Ángulo de talud: 48° (global)

Valorizaciones por destino:

  • Valor a planta: cada bloque se valoriza como si su destino final es planta.
  • Valor a botadero: cada bloque se valoriza como si su destino final es botadero.

Capacidades por período:

  • Mina: 400.000 [Ton]
  • Planta: 350.000 [Ton]

 

Agendamiento eligiendo el proceso óptimo 

RESULTADOS HEURISTICA POR VENTANAS (STW) (ventana:1, paso: 1)

En las figuras a continuación se muestran distintas vistas del agendamiento en el tiempo obtenido mediante STW: (i) vista isométrica, (iii) vista en planta y (v) vista en sección. Note que esta estrategia requirió 17 periodos para agendar la producción. Por otra parte, se muestra el agendamiento por destino: (ii) vista isométrica, (iv) vista en planta, y (vi) vista en sección. En amarillo (1) se muestra lo que es enviado a planta, en rojo (2) lo que es enviado a botadero, y en azul (0) lo que nunca es extraído.

 stw18 iso sched  stw18 iso dest
(i)  (ii) 
 stw18 z sched stw18 z dest 
 (iii)  (iv)
stw18 slice sched stw18 slice dest
(v) (vi)

 

PLAN DE PRODUCCIÓN USANDO HEURISTICA POR VENTANAS (STW)

Se muestran los planes de producción (tonelajes de mineral y lastre), además de las leyes medias de cada metal por periodo.

stw18 pps

 

RESULTADOS HEURISTICA TOPOSORT (TOPO)

En las figuras a continuación se muestran distintas vistas del agendamiento en el tiempo obtenido mediante TOPO: (i) vista isométrica, (iii) vista en planta y (v) vista en sección. Note que esta estrategia requirió los 18 periodos para agendar la producción, a diferencia de STW, que lo hizo en un periodo menos. Por otra parte, se muestra el agendamiento por destino: (ii) vista isométrica, (iv) vista en planta, y (vi) vista en sección. En amarillo (1) se muestra lo que es enviado a planta, en rojo (2) lo que es enviado a botadero, y en azul (0) lo que nunca es extraído. Note las diferencias de asignación entre ambas estrategias, STW y TOPO: muchos de los bloques que STW envía a botadero, TOPO los clasifica como mineral a planta.

topo18 iso sched topo18 iso dest
(i) (ii)
topo18 z sched topo18 z dest
(iii) (iv)
topo18 slice sched topo18 slice dest
(v) (vi)

 

PLAN DE PRODUCCIÓN USANDO HEURISTICA TOPOSORT (TOPO)

Se muestran los planes de producción (tonelajes de mineral y lastre), además de las leyes medias de cada metal por periodo.

topo18 pps

 

Comparación y conclusión

Considerando los resultados obtenidos por ambas heurísticas, se ve claramente que la heurística toposort agendó la producción en 18 periodos, tal como se pidió, mientras que la heurística por ventanas acortó la vida de la mina a 17 años. Mirando los planes de producción, TOPO pudo retrasar la entrada de lastre en los primeros periodos, clasificando como mineral todo el material extraído entre los periodos 1 y 4, además del 6, y desde ahí extrayendo de manera gradual el lastre hasta completar la capacidad de mina en los periodos finales. Por su parte, STW hace una extracción homogenea, saturando tanto las capacidades de mina y planta en la mayoría de los periodos (salvo 16 y 17). Por una parte, la característica "míope" de STW no le permite tomar decisiones mirando lo que podría hacer en periodos posteriores, tomando en cada periodo la mejor decisión posible. Por otra parte, TOPO es capaz de tomar una decisión "global" a partir de la información obtenida de la solución continua del problema original, por lo que en cierto sentido toma una decisión más robusta, pues para decidir sobre un periodo particular, lo hace mirando qué pasará en periodos posteriores.  

En términos de valor descontado acumulado (VAN), la estrategia TOPO muestra un VAN ínfimamente mayor al VAN de STW, de hecho, un 0,23% mayor. 

VAN STW  : 1.386 [MUSD]

VAN TOPO: 1.390 [MUSD]

En general, este tipo de estrategias no pueden asegurar encontrar una solución óptima al problema de agendamiento, sin embargo, son una buena alternativa considerando que en problemas reales de gran tamaño (en cuanto a número de bloques, periodos y destinos) se hace impracticable obtener la solución óptima por tratamiento directo del problema.