http www iro umontreal caferland http www researchandpractise

  • Slides: 56
Download presentation
http: //www. iro. umontreal. ca/~ferland/ http: //www. researchandpractise. com/vrp/

http: //www. iro. umontreal. ca/~ferland/ http: //www. researchandpractise. com/vrp/

Capacitated Open Pit Mining Problem Semya Elaoud, Sfax University (Tunesia) Jacques A. Ferland, University

Capacitated Open Pit Mining Problem Semya Elaoud, Sfax University (Tunesia) Jacques A. Ferland, University of Montreal Jonathan Bellemare, University of Montreal Jorge Amaya, University of Chile

RIOT Mining Problem web site: http: //riot. ieor. berkeley. edu/riot/Applications/OPMInteractive. html the net value

RIOT Mining Problem web site: http: //riot. ieor. berkeley. edu/riot/Applications/OPMInteractive. html the net value of extracting block i objective function Maximal Open Pit problem: to determine the maximal gain expected from the extraction

Maximal pit slope constraints to identify the set Bi of predecessor blocks that have

Maximal pit slope constraints to identify the set Bi of predecessor blocks that have to be removed before block i

Maximal pit slope constraints to identify the set Bi of predecessor blocks that have

Maximal pit slope constraints to identify the set Bi of predecessor blocks that have to be removed before block i

Maximal pit slope constraints to identify the set Bi of predecessor blocks that have

Maximal pit slope constraints to identify the set Bi of predecessor blocks that have to be removed before block i

Use the open pit graph G = (V, A) to specify the maximal pit

Use the open pit graph G = (V, A) to specify the maximal pit slope constraints The maximal pit slope constraints:

 • (MOP) equivalent to determine the maximal closure of G = (V, A)

• (MOP) equivalent to determine the maximal closure of G = (V, A) • Equivalent to determine the minimum cut of the associated Picard’s graph where The maximal open pit is equal to N* = (S – {s})

Scheduling block extraction Account for operational constraints: Ct the maximal weight that can be

Scheduling block extraction Account for operational constraints: Ct the maximal weight that can be extracted during period t and for the discount factor during the extracting horizon: discount rate period

pi weight of block i the net value of extracting block i N can

pi weight of block i the net value of extracting block i N can be replaced by the maximal open pit N* = (S – {s})

Scheduling block extraction ↔ RCPSP • Open pit extraction ↔ project • Each block

Scheduling block extraction ↔ RCPSP • Open pit extraction ↔ project • Each block extraction ↔ activity • Precedence relationship derived from the open pit graph

Scheduling block extraction ↔ RCPSP

Scheduling block extraction ↔ RCPSP

Scheduling block extraction ↔ RCPSP

Scheduling block extraction ↔ RCPSP

Scheduling block extraction ↔ RCPSP • Open pit extraction ↔ project • Each block

Scheduling block extraction ↔ RCPSP • Open pit extraction ↔ project • Each block extraction ↔ activity • Precedence relationship derived from the open pit graph • Reward associated with activity (block) i depends of the extraction period t

Solution encoding and decoding

Solution encoding and decoding

Decoding a block list into a schedule Serial decoding • Initiate the first extraction

Decoding a block list into a schedule Serial decoding • Initiate the first extraction period t = 1 • During any period t: - The next block to be extracted is the first block in the rest of the block list (including the blocks not extracted yet) having all their predecessors already extracted such that the capacity Ct is not exceeded by its extraction. Include this block in the newsol block list. - If no such block exists, then a new extraction period (t + 1) is initiated.

Metaheuristic solution approach

Metaheuristic solution approach

Outline of the solution approach

Outline of the solution approach

Initial solution

Initial solution

First neighborhood moving one ore block

First neighborhood moving one ore block

Unique ore block process on cluster i

Unique ore block process on cluster i

Second neighborhood moving multiple ore blocks

Second neighborhood moving multiple ore blocks

Multiple ore block process on cluster i

Multiple ore block process on cluster i

Implementation of a metaheuristic procedure

Implementation of a metaheuristic procedure

Numerical experimentation

Numerical experimentation

Second encoding of the solution and Particle Swarm Solution Approach

Second encoding of the solution and Particle Swarm Solution Approach

Genotype representation of solution Similar to Hartman’s priority value encoding for RCPSP priority of

Genotype representation of solution Similar to Hartman’s priority value encoding for RCPSP priority of scheduling block i extraction

Decoding of a representation PR into a solution x • Serial decoding to schedule

Decoding of a representation PR into a solution x • Serial decoding to schedule blocks sequentially one by one to be extracted • To initiate the first extraction period t = 1: remove the block among those having no predecessor (i. e. , in the top layer) having the highest priority. • During any period t, at any stage of the decoding scheme: the next block to be removed is one of those with the highest priority among those having all their predecessors already extracted such that the capacity Ct is not exceeded by its extraction. If no such block exists, then a new extraction period (t + 1) is initiated.

Priority of a block • Consider its net value bi and impact on the

Priority of a block • Consider its net value bi and impact on the extraction of other blocks in future periods • Block lookahead value (Tolwinski and Underwood) determined by referring to the spanning cone SCi of block i

Genotype priority vector generation • Several different genotype priority vectors can be randomly generated

Genotype priority vector generation • Several different genotype priority vectors can be randomly generated with a GRASP procedure biased to give higher priorities to blocks i having larger lookahead values • Several feasible solutions can be obtained by decoding different genotype vectors generated with the GRASP procedure.

Particle Swarm Procedure • Evolutionary process evolving in the set of genotype vectors to

Particle Swarm Procedure • Evolutionary process evolving in the set of genotype vectors to converge to an improved feasible solution • Initial population P of M genotype vectors (individuals) generated using GRASP • Denote the best achievement of the individual k up to the current iteration the best overall genotype vector achieved up to the current iteration

Particle Swarm Procedure • Denote the best achievement of the individual k up to

Particle Swarm Procedure • Denote the best achievement of the individual k up to the current iteration the best overall genotype vector achieved up to the current iteration Modification of the individual vector k at each iteration

Particle Swarm Procedure • Denote the best achievement of the individual k up to

Particle Swarm Procedure • Denote the best achievement of the individual k up to the current iteration the best overall genotype vector achieved up to the current iteration Modification of the individual vector k at each iteration

Numerical experimentation

Numerical experimentation