ZA200502772B - System and method(s) of mine planning, desing and processing - Google Patents
System and method(s) of mine planning, desing and processing Download PDFInfo
- Publication number
- ZA200502772B ZA200502772B ZA200502772A ZA200502772A ZA200502772B ZA 200502772 B ZA200502772 B ZA 200502772B ZA 200502772 A ZA200502772 A ZA 200502772A ZA 200502772 A ZA200502772 A ZA 200502772A ZA 200502772 B ZA200502772 B ZA 200502772B
- Authority
- ZA
- South Africa
- Prior art keywords
- blocks
- constraints
- determining
- accordance
- block
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims description 160
- 238000012545 processing Methods 0.000 title claims description 21
- 239000000463 material Substances 0.000 claims description 89
- 238000013461 design Methods 0.000 claims description 63
- 239000002699 waste material Substances 0.000 claims description 39
- 238000000605 extraction Methods 0.000 claims description 24
- 230000000644 propagated effect Effects 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 12
- 230000004931 aggregating effect Effects 0.000 claims description 10
- 230000001902 propagating effect Effects 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 4
- 239000012535 impurity Substances 0.000 claims 2
- 238000013102 re-test Methods 0.000 claims 1
- 239000000203 mixture Substances 0.000 description 79
- 238000009472 formulation Methods 0.000 description 78
- 239000000243 solution Substances 0.000 description 59
- 238000005065 mining Methods 0.000 description 35
- 238000013459 approach Methods 0.000 description 30
- 230000002776 aggregation Effects 0.000 description 24
- 238000004220 aggregation Methods 0.000 description 24
- 238000005520 cutting process Methods 0.000 description 23
- 230000008569 process Effects 0.000 description 19
- 238000004422 calculation algorithm Methods 0.000 description 13
- 230000008901 benefit Effects 0.000 description 10
- 238000010420 art technique Methods 0.000 description 8
- 230000003247 decreasing effect Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 230000001965 increasing effect Effects 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 6
- 239000011435 rock Substances 0.000 description 5
- 238000012163 sequencing technique Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 231100000817 safety factor Toxicity 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- RTZKZFJDLAIYFH-UHFFFAOYSA-N Diethyl ether Chemical compound CCOCC RTZKZFJDLAIYFH-UHFFFAOYSA-N 0.000 description 1
- 101100328463 Mus musculus Cmya5 gene Proteins 0.000 description 1
- 208000009989 Posterior Leukoencephalopathy Syndrome Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000005267 amalgamation Methods 0.000 description 1
- 239000010953 base metal Substances 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 230000003467 diminishing effect Effects 0.000 description 1
- 238000011143 downstream manufacturing Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
SYSTEM AND METHOD(S) OF MINE PLANNING, DESIGN AND PROCESSING
) The present invention relates to the field of extracting resource(s) from a particular location. In particular, the present invention relates to the planning, . 5 design and processing related to a mine location in a manner based on enhancing the extraction of material considered of value, relative to the effort “os and/or time in extracting that material.
In the mining industry, once material of value, such as ore situated below the surface of the ground, has been discovered, there exists a need to extract that material from the ground.
In the past, one more traditional method has been to use a relatively large open cut mining technique, whereby a great volume of waste material is removed from the mine site in order for the miners to reach the material considered of 16 value. For example, referring to Figure 1, the mine 101 is shown with its valuable material 102 situated at a distance below the ground surface 103. In the past, most of the (waste) material 104 had to be removed so that the valuable material 102 could be exposed and extracted from the mine 101. In the past, this waste material was removed in a series of progressive layers 105, which are ever diminishing in area, until the valuable material 102 was exposed for extraction.
This is not considered to be an efficient mining process, as a great deal of waste material must be removed, stored and returned at a later time to the mine site 101, in order to extract the valuable material 102. It is desirable to reduce the volume of waste material that must be removed prior to extracting the valuabie material.
The open cut method exemplified in Figure 1 is viewed as particularly inefficient where the valuable resource Is located to one side of the pit 105 of a desirable mine site 101. For example, Figure 2 illustrates such a situation. The valuable material 102 is located to one side of the pit 105. In such a situation, itis _ not considered efficient to remove the waste material 104 from region 206, that is where the waste material is not located relatively close to the valuable material 102, but it is considered desirable to remove the waste material 104 from region : 207, that Is where it Is located nearer to the valuable material 102. This then rings other considerations to the fore. For example, it would be desirable to determine the boundary between regions 206 and 207, so that not too much } undesirable waste material is removed (region 2086), yet enough is removed to ensure safety factors are considered, such as cave-ins, etc. This then leads to a further consideration of the need to design a ‘pit’ 105 with a relatively optimal design having consideration for the location of the valuable material, relative to the waste material and other issues, such as safety factors.
This further consideration has led to an analysis of pit design, and a technique of removing waste material and valuable material called ‘pushbacks’.
This technique is illustrated in Figure 3. Basically, the pit 105 is designed to an extent that the waste material 104 to be removed is minimised, but still enabling extraction of the valuable material 102. The technique uses ‘blocks’ 308 which represent smaller volumes of material. The area proximate the valuable material is divided into a number of blocks 308. It is then a matter of determining which blocks need to be removed in order to enable access to the valuable material 102. This determination of ‘blocks 308’, then gives rise to the design or extent of the pit 105.
Figure 3 represents the mine as a two dimensional area, however, it should be appreciated that the mine is a three dimensional area. Thus the blocks 308 to be removed are determined in phases, and cones, which represent more accurately a three dimensional ‘volume’ which volume will ultimately form the pit 105.
Further consideration can be given to the prior art situation lllustrated In
Figure 3. Consideration should be given to the scheduling of the removal of blocks. In effect, what is the best order of block removal, when ather business aspects such as timefvalue and discounted cash flows are taken into account?
There is a need to find a relatively optimal order of black removal which gives a relatively maximum value for a relatively minimum effort/time.
Attempts have been made in the past to find this "optimum’ block order by determining which block(s) 308 should be removed relative to a "violation free’ N order. Turning to the illustration in Figure 4, a pit 105 is shown with valuable : material 102. For the purposes of discussion, if it was desirable to remove block - 414, then there Is considered to be a ‘violation’ if we determined a schedule of block removal which started by removing block 414 or blocks 414, 412 & 413 before blocks 409, 410 and 411 were removed. In other words, a violation free schedule would seek to remove other blacks 409, 410, 411, 412 and 413 before block 414. (it is important to note that the block number does not necessarily indicate a preferential order of block removal).
It can also be seen that this block scheduling can be extended to the entire pit 105 in order to remove the waste material 104 and the valuable material 102.
With this violation free order schedule in mind, prior art attempts have been made. Figure 5 illustrates one such attempt. Taking the blocks of Figure 4, the blocks are numbered and sorted according to a ‘mineable block order’ having regard to practical mining techniques and other mine factors, such as safety etc and is illustrated by table 515. The blocks In table 515 are then sorted 516 with regard to Net Present Value (NPV) and is based on push back design via Life-of- mine NPV sequencing, taking Into account obtaining the most value block from the ground at the earliest time. To illustrate the NPV sorting, and turning again to
Figure 4, there is a question as which of blocks 409, 410 or 411 should be removed first. All three blocks can be remaved from the point of view of the ability to mine them, but it may, for example, be more economic to remove block 410, before block 409. Removing blocks 409, 410 or 411 does not lead to 'violations' thus consideration can be given to the order of block removal which is more economic.
The NPV sorting Is conducted in a manner which does not lead to violations of the ‘violation free order’, and provides a table 517 listing an ‘executable block order’. In other words, this prior art technique leads to a listing of blocks; in an order which determines their removal having regard to the ability to mine them, and the economic return for doing so.
Furthermore, a number of prior art techniques are considered to take a relatively simple view of the probfems confronted by the mine designer in a ‘real world" mine situation. For example, the size, complexity, nature of blocks, grade, _ slope and other engineering constraints and time taken fo undertake a mining operation is often not fully taken into account in prior art techniques, leading to computational problems or errors in the mine design. Such errors can have . significant financial and safety implications for the mine operator.
With regard to size, for example, prior art techniques fail to adequately take account of the size of a "block’. Depending on the size of the overall project, a ‘block’ may be quite large, taking some weeks, months or even years to mine. If . this is the case, many assumptions made in prior art techniques fall to give sufficient accuracy for the modern day business environment, : Given that many of the mine designs are mathematically and computational complex, according to prior art techniques, If the size of the blocks were raduced for greater accuracy, the result will be that either the optimisation techniques used will be time in feasible ( that is they will take an inordinately long time to complete), or other assumptions will have to be made conceming aspects of the mine design such as mining rates, processing rates, etc which will result in a decrease the accuracy of the mine design solution.
Some examples of commercial software do use mixed integer programming engines, however, the method of aggregating blocks requires 16 further improvement. For example, it is considered that product ‘ECSI Maximiser by ECS International Pty Ltd uses a form of integer optimisation in thelr pushback design, but the optimisation is local in time, and it's problem formulation is considered too large to optimise globally over the life of a mine. Also the product ‘MineMax’ by MineMAX Ptd Ltd may be used to find a rudimentary optimal block sequencing with a mixed integer programming engine, however it is considered that it's method of aggregation does not respect slopes as is required in many situations. ‘MineMax’ also optimises locally in time, and not globally. Thus, where there are a large number of variables, the user must resort to subdividing the pit into separate sections, and perform separate optimisations on each section, and thus the optimisation is not global over the entire pit. itis considered desirable to have an optimisation that is global in both space and time.
Dynamic Programming Approach
The Lerchs-Grossman graph-theoretic algorithm (H. Lerchs & I. Grossman, "Optimum Design of Open-Pit Mines”, Transactions CIM, 1965) has been proved _ to give a relatively exact solution to the ultimate pit problem for an open-cut mine in three dimensions. Lerchs and Grossman also presents a dynamic programming approach to the problem in two dimensions, which has since been extended to three dimensions. However, solution of the three-dimensional graph theoretic algorithm is computationally inefficient in practical cases.
Linear Programming Approach ) There Is a linear program (LP), as presented by Underwood and Tolwinski 9 (R. Underwood & B. Tolwinski, “A mathematical programming viewpoint for : solving the ultimate pit problem”, EJOR, 1998). The availability of CPLEX (by llog, www.ilog.com) as a powerful LP solver motivates investigation of the LP approach to the ultimate pit problem.
The ultimate pit problem can be modelled as an integer program (IP), where a value of 1 is assigned to blocks included in the ultimate pit, and a value of 0 is assigned otherwise. The IP formulation for the problem is then as follows.
Let %= 1, ifblockiis included in the ultimate pit 0, otherwise 16 Then max 3 vx,
SL. <x, VieP@) x {0,1} Vi ......equation 1 where vi is the value assigned to block i x is the decision variable that designates whether block i is included in the ultimate pit or not
P(i) is the set of predecassor blocks of block i.
One objective is to maximise the net value of the material removed from the pit. Consider that the only constraints are precedence constraints, which enforce the requirement of safe wall slopes in the mine. In fact, this IP formulation has the property of total unimodularity. That is, the solution of the LP relaxation of this formulation will be integral (i.e. a set of -0's and 1's). This is -an extremely - desirable property for an integer program. It allows the IP to be solved as an LP using the Simplex method. This leads to greatly increased solution efficiency In terms of both CPU time and memory requirements. The exact mathematical formulation of the linear programming approach to the ultimate pit problem Is therefore max Pp vx,
Il 8.2. x, $x, VjelP(@) 0<x, <1 Wi ......equation 2
This is the ideal approach to solve the problem, and is considered to give the optimal solution in every case. Unforlunately, implementation of this exact formulation in CPLEX fails to solve for mining projects of realistic size. Since the optimisation is carried out at the block level, and there is a constraint for every precedence arc for each block, a very large number of constraints are applied.
For example, if a mine has 198,917 blocks, and after CPLEX performs pre- processing on the formulation, the resulting reduced LP still has 1,676,003 constraints. CPLEX attempts to solve this formulation using the dual simplex method, generally recognized as the most efficient method for solving linear programs of this size. However, in the case of the example mine, CPLEX was found to crash during the solution process due to the very large number of constraints. Inversion of a constraint matrix of this magnitude (as required for : converting solutions obtained from the dual simplex method back into primal space) is considered to place too great a memory requirement on the system.
There still exists a need, however, to improve prior art techniques. Given that mining projects, on the whole, are relatively large scale operations, even small improvements in pror art techniques can represent millions of dollars in savings, and / or greater productivity and / or safety. it is desirable to provide an improved mine design.
An object of the present invention is to provide an improved method of pit design, which takes into account slope constraints. _
Another object of the present invention is to provide an improved method of determining a cluster.
A further object of the present invention is to determine which blocks of a mine pit provide a relative maximum net value of material, also having regard to practical limitations, such as slope constraints.
Yet another object of the present invention is to alleviate at least one disadvantage of the prior art.
Any discussion of documents, devices, acts or knowledge in this specification is included to explain the context of the invention. It should not be taken as an admission that any of the material forms a part of the prior art base or the common general knowledge in the relevant art in Australia or elsewhere on or before the priority date of the disclosure and claims herein.
The present invention provides, In a first Inventive aspect, a method of and apparatus for determining slope constraints related to a design configuration for extracting material from a particular location, the method including the steps of determining a selected volume of material to be extracted, dividing at least a portion of the selected volume into blacks, forming a plurality of cones, at least one cone from each block, and determining from the cones, a clump having a corresponding slope constraint.
Preferably, the cone is propagated upwards using precedence arcs.
The present aspect also provides a method of determining slope constraints related to a design configuration for extracting material from a particular location, in which precedent arcs emanating from a selected block(s) are used to establish, at least in part, slope constraints.
The present aspect also provides a mine designed in accordance with the method as disclosed herein.
The present aspect further provides a computer program product including a computer usable medium having computer readable program code and computer readable system code embodisd on said medium for determining slope constraints related to a design configuration for extracting material from a _ particular location within a data processing system, the computer program product including computer readable code within said computer usable medium for performing the method as disclosed herein.
In essence, the present invention, referred to as Propagation of clusters and formation of clumps, forms relatively minimal Inverted cones with clusters at their apex and intersects these cones to form clumps, or aggregations of blocks that respect slope constraints. Advantageously, it has been found that $5 aggregating the small blocks in an intelligent way serves to reduce the number of atoms” variables to be fed into the mixed integer programming engine. The clumps allow relatively maximum flexibility In potential mining schedules, while keeping variable numbers to a minimum. The collection of clumps has three important properties. Firstly, the clumps allow access to all the targets as quickly as possible (minimality), and secondly the clumps allow many possible orders of access to the identified ore targets (flexibility). Thirdly, because cones are used, and due to the nature of the cone(s), an extraction ordering of the clumps that is feasible according to the precedence arcs will automatically respect and accommodate minimum slope constraints. Thus, the slope constraints are 16 automatically built into this aspect of invention.
In other words, the present invention provides that clumps are determined from the overlap of cones. The cones are preferably ‘minimal’.
The present invention provides, in a second inventive aspect, a method of and apparatus for determining a cluster of material, the method including allocating at least a portion of the material between a plurality of blocks, detemmining a first attribute related to co-ordinates corresponding to each block, assigning the first attribute to each corresponding block, determining a second attribute related to the plurality of blocks, and aggregating at least two of the plurality of blocks in accordance with the first attribute and the second attribute.
In essence, the second related aspect of invention, referred to as Initial
Identification of Clusters, aggregates a number of blocks Into collections or clusters. The clusters preferably more sharply identify regions of high-grade and a low-grade materials, while maintaining a spatial compactness of a cluster. The clusters ara formed by blocks having certain x, vy, z spatial coordinates, combined with another coordinate, representing a number of selected values, such as grade . or value. The advantage of this is to produce inverted cones that are relatively tightly focused around regions of high grade so as not to necessitate extra stripping.
In other words, where there is an ore bady having a number of blocks, the present invention deals with building cones and clumps etc from the information known about the ore body and it's blocks.
The present invention provides, in a third inventive aspect, a method and . apparatus of determining characteristics of a selected portion of material, the method includi ng determining the contents of the selected portion of material, and identifying region(s) of material within the selected portion according to at least one of a plurality of characteristic(s).
In essence, a third related aspect of invention, referred to as splitting of waste and ore in clumps, is based on the realisation that clumps contain both ore blocks and waste blocks. Many Integer programs assume that the value is distributed uniformly within a clump. This is, however, not true. Typically, clumps will have higher value near their base. This is because most of the value is lower underground while closer to the surface: one tends to have more waste blocks.
By splitting the clump into relatively pure waste and desirable material, the assumption of uniformity of value for each portion of the clump is more accurate.
In other words, the present invention reflects the consideration to determine, where necessary, block ‘grade’. If the ore is above a certain value, then the cone may be divided into smaller canes, and re-iterated for more precise determination and extraction.
The present invention provides, in a fourth inventive aspect, a method of and apparatus for analysing a selected volume of material, the material being at least partially comprised of a plurality of blocks, the methad including the steps of clumping a number of blocks together, and analysing the selected volume of material based on the clumped blocks.
In essence, a fourth related aspect of invention, referred to as Aggregation of blocks into clumps; high-level ideas, reduces the number of variables io a relatively manageable amount for use in current technology of integer i ' programming engines. Advantageously, this aspect enables the use of an integer programming engine and the ability to Incorporate further constraints such es : mining, processing, and marketing capacities, and grade constraints.
The present invention provides, in a fifth inventive aspect, a method of determining a selected group of blocks of a mine pit which are capable of being mined, the method including the steps of selecting a plurality of blocks, and determining a relative value and canstraints applicable to the selected blocks in accordance with any one of the equations 3, 4 or 9 as disclosed herein.
The present invention also provides the method as described above and including the further step of testing for violations.
The present invention also seeks to reiterate the selection and determination of value and constraints of blocks in order to obtain a group of blocks which have a relative optimal mining value.
In essence, the present aspect, in one form, utilises aggregating algorithm(s) to determine a selected group of blocks which are to be mined, where the selection of blocks to be included into the group of blocks is made relative to value and constraints applicable to the blocks. The present invention, in another aspect further tests for violations, and iteratively recaiculates until substantially all violations are removed. Given a block model of an ore body containing value-in-ground and designated slope constraints, the ultimate pit problem concerns the determination of the shape of the final pit of the mine. it is assumed that all the material can be removed at once. That is, the effect of time on the value of the ore body is not considered, In terms of mine scheduling, the ultimate pit can be used as the initial collection of blocks on which a scheduling algorithm is run. In this respect, the ultimate pit is the largest possible final pit that can be realised following scheduling of removal of the ore body. The case considered throughout this disclosure is that of base metals but aiso has application to blended products or stochastic elements of open-pit mining.
In other words, the present invention is used to determine how to split a relatively large ore body into clump(s). The present invention can be used to : ensure that the clump or ore body is not too large, computationally. for example for practical consideration with the use of existing algorithms.
Other related aspects of invention, include: B
In essence, one related aspect of invention, referred to as Generic
Klumpking, is a method of mine design that firstly, is considered a clever choice of aggregation to reduce the number of varisbles via a spatialfvalue clustering
“and propagation to form clumps. Secondly, the Inclusion of mining and . processing constraints in an integer program based around the clump variables to ultimately produce an optimal block sequence. Thirdly, the rapid loop of clustering ) blocks in this optimal sequence according to spacefime of extraction and propagating these clusters to form pushbacks, interrogating them for value and ’ mineability, and adjusting clustering parameters as needed.
In essence, another relaied aspect of invention, referred to as
Determination of a block ordering from a clump ordering, turns a clump ordering into an ordering of blocks. This is, in effect, a de aggregation. Using techniques disclosed herein, the integer program engine was used on the relatively small number of clumps, and thus the result can now be transtated back into the large number of small blocks.
In essence, still another related aspect of invention, referred to as ‘fuzzy clustering; second identification of clusters for pushback design, clusters blocks according to their spatial position and their time of extraction. This is considered necessary because if pushbacks were formed from the block sequence in its raw form, the pushbacks would be generally highly fragmented and considered non- mineable. The clustering gives control over the connectivity and mineability of the resulting pushbacks.
In essence, still another related aspect of invention, referred to as fuzzy clustering; alternative 1, clusters blacks according to their spatial position and their time of extraction. The clusters may be controlled to be a certain size, or have a certain rock tonnage or ore tonnage. The shapes of the clusters may be controlled through parameters that balance the space and the time coordinate.
The advantage of shape control is to. produce pushbacks that are mineable and not fragmented. The advantage of size contro) is the ability to control stripping ratios in years where the mill may be operating under capacity.
In essence, a further related aspect of Invention, referred to as fuzzy clustering; alternative 2, propagates inverted cones from the clusters identified in the secondary clustering. The clusters in the secondary clustering are time - ordered, and the propagation occurs In this time order, with no intersections of inverted cones allowed. Advantageously, this provides the ability to extract pushbacks from the block ordering that are well connected and mineable, while retaining the bulk of the NPV optimality of the block sequence.
In essence, still a further related aspect of invention, referred fo as fuzzy : clustering; alternative 3, provides the creation of a feedback loop of clustering,
S propagating to find pushbacks, valuing relatively quickly, and then feeding this : information back into the choice of clustering parameters. The advantage of this is that the effect of different dustering parameters may be very quickly checked for NPV and mineability. It Is heretofore been virtually impossible to evaluate a pushback design for NPV and mineability before it has been constructed, and the fast process loop of this aspect allows many high-quality pushbacks designs fo be constructed and evaluated (by the human eye in the case of mineability).
Other aspects and preferred aspects are disclosed in the specification and } or defined in the appended claims.
The method(s), systems and techniques disclosed in this application may be used in conjunction with prior art integer programming engines. Many aspects of the present disclosure serve to improve the performance of the use of such engines and the use of other known mine design techniques.
The present invention may be used, for example, by mine planners to design relatively optimal pushbacks for open cut mines. Advantageously, the present invention is considered is different to prior art pushback design software in that: ’ . The praesent Invention does not use either of the most common pit design algorithms (Lerchs-Grossmann or Floating Cone) but instead uses a unique concept of optimal “clump” sequencing to develop an optimal block sequence that Is then used as a basis for pushback design. . The design is relatively optimal with respect to properly discounted block values. No other pushback design software is considered to correctly allow for the effect of time (viz: block value discounting) in the pushback design step. Traditional phase designs ignore medium grade ore pods close to the surface with good NPV whilst focussing on higher value pods B that may be deeply buried. . The present invention can properly address the so-called ‘“Whittle-gap” problem where consecutive Lerchs-Grossmann shells can be very far apart, offering little temporal information. The present invention obtains relatively complete and accurate temporal information on the block ordering. ® Process and mining constraints can be explicitly incorporated into the pushback design step. . The planner can rapidly design and value pushbacks that have different topologies, the trade-off being between pits with high NPV, but with difficult-to-mine (eg: ring) pushback shapes, and those with mare mineable pushback shapes but lower NPV. The advantage of the more mineable pushback shapes is that much less NPV will be wasted in enforcing minimum mining width and In accommodating pit access (roads and berms). . The ability to quickly generate and evaluate a number of different sets of candidate pushback designs is a feature not allowed in traditional pushback design software where dasign options are usually fairly limited (eg: the amalgamation of adjacent Whittle shells into a single pushback) ) Various aspects of the present invention also serve to improve the use of existing integer programming engines, such as “cplex” by ILOG.
Throughout the specification: 1. a ‘collection’ is a term for a group of objects, 2. a ‘cluster is a collection of ore blocks or blocks of otherwise desirable material that are relatively close to one another in terms of space and / or other attributes, 3. a ‘clump’ Is formed from a cluster by first producing a substantially minimal 25° inverted cone extending from the cluster to the surface of the pit by propagating all blocks in the cluster upwards using the arcs that describe the minimal slope constraints. Each cluster will have its own minimal inverted cone. These minimal Inverted cones are then intersect with one another and the intersections form clumps, and 4. an ‘aggregation’ is a term, although mostly applied to collections of blocks N that are spatially connected (no “holes” in them). For example, a clump may. be an aggregation, or may be "Super blocks" that are larger cubes made by joining together smaller cubes or blocks.
5. reference to block constraints equally implies reference to arc constraints. 6. a block may also refer to a number of blocks.
Further disclosure, objects, advantages and aspects of the present application may be bettar understoad by those skilled in the relevant art with reference to the following description of preferred embodiments taken in ~. conjunction with the accompanying drawings, in which:
Figures 1 to 5 illustrate prior art mining techniques,
Figure 6 illustrates, schematically, a flow chart outlining the overall process according to one aspect of invention,
Figure 7 illustrates schematically the identification of clusters,
Figure 8 illustrates schematically cone propagation in pit design,
Figure 9 illustrates schematically the splitting or ore from waste material,
Figure 10 illustrates an example of ‘fuzzy clustering’ in a mine site,
Figures 11a, 11b and 11c¢ illustrate a secondary ¢lustering, propagation, and NPV valuation process,
Figure 12 illustrates a comparison between outcomes of equations 2 and 4,
Figure 13 illustrates a vertical cross-section of a pit design using equation 2,
Figure 14 illustrates a vertical cross-section of a pit design using equation 4,
Figure 15 illustrates an example portion of a pit,
Figures 16 and 18 illustrate a plane view through a pit using the cutting plane formulation {equation 9), and : Figures 17 and 19 illustrate the same view as that of Figures 16 and 18 but for the use of the LP relaxation of the aggregated formulation (equation 4).
In order to more fully describe the present invention, a number of related aspects will also be described. In this way, the reader can gain a better a. understanding of the context and scope of the present invention.
1. Generic KlumpKing
Figure 6 illustrates, schematically an overall representation of one aspect of Invention.
Although specific aspects of various elements of the overali flow chart are discussed below in more detail, it may be helpful to provide an outline of the flow chart lilustrated in Figure 6.
Block model 601, mining and processing parameters 602 and slope constraints 603 are provided as input parameters. When combined, precedence arcs 604 are provided. For a given block, arcs will point to other blocks that myst be removed before the given block can be removed.
As typically, the number of blocks can be very large, at 605, blocks are aggregated into larger collections, and clustered. Cones are propagated from respective clusters and clumps are then created 606 at intersections of cones.
The number of clumps is now much smaller than the number of blocks, and clumps Include slope constraints. At 607, the clumps may then be scheduled in a manner according to specified criteria, for example, mining and processing constraints and NPV. It is of great advantage that the scheduling occurs with clumps (which number much less than blocks). It is, in part, the reduced number of clumps that provides a relative degree of arithmetic simplicity and / or reduced requirements of the programming engine or algorithms used to determine the schedule. Following this, a schedule of individual block order can be determined from the clump schedule, by de-aggregating. The step of polish at 608 is optional, but does improve the value of the block sequence.
From the block ordering, pushbacks can be designed 608. Secondary clustering can be undertaken 610, with an additional fourth co-ordinate. The fourth co-ordinate may be time, for example, but may also be any other desirable value or parameter. From here, cones are again propagated from the clusters, but in a sequence commensurate with the fourth co-ordinate. Any blocks already assigned to previously propagated cones are not included In the next cone propagation. Pushbacks are formed 611 from these propagated cones. -
Pushbacks may be viewed for mineability 612. An assessment as to a balance between mineability and NPV can bs made at 613, whether in accordance with a
Predetermined parameter or not. The pushback design can be repeated if necessary via path 614.
Other consideration can also be taken into account, such as minimum mining width 615, and validation 616. Balances can be taken into account for mining constraints, downstream processing constraints and / or stockpiling options, such as blending and supply chain determination and / or evaluation.
The following description focuses on a number of aspects of invention which reside within the overall flow chart disclosed above. For the purposes of
Figure 6, sections 2 and 5 are associated with 605, sections 3, 4 and 5 are associated with 606, sections 4, 6 are assoclatad with 607, sections 7 and 7.3 are associated with 610, sections 7.2 and 7.3 are associated with 611, section 7.3 is associated with 612, 613 and 614, and sections 7, 7.1, 7.2 and 7.3 are associated with 609. 1.1 Inputs and preliminaries
Input parameters include the block model 601, mining and processing parameters 602, and slope constraints 603. Slope regions (eg. physical areas or zones) are contalned in 601; slope parameters (eg. slopes and bearings for each zane) are contained in 602.
The block madel 601 contains information, for example, such as the value of a block in dollars, the grade of tha block in grams per tonne, the tonnage of rock in the block, and the tonnage of ore in the block.
The mining and processing parameters 602 are expressed in terms of tonnes per year that may be mined or processed subject to capacity constraints.
The slope constraints 603 contain information about the maximal slope around in given directions about a particular block.
The slope constraints 603 and the block mode! 601 when combined give rise to precedence arcs 604. For a given block, arcs will point from the given block to all other blocks that must be removed before the given block. The - number of arcs Is reduced by storing them in an inductive, where, for example, in Co two dimensions, an inverted. cone of blocks may be described by every block pointing to the three blocks centred immediately above it. This principle can also be applied to three dimensions. If the inverted cone is large, for example having a depth of 10 : the number of arcs required would be 100; one for each block.
However, using the inductive rule of "point to the three blocks centred directly above you", the entire inverted cone may be described by only three arcs instead of the 100. In this way the number of arcs required to be stored is greatly reduced. As block models typically contain hundreds of thousands of blocks, with each block containing hundreds of arcs, this data compression is considered a ) significant advantage. 1.2 Producing an optimal block ordering
The number of blocks In the block model 601 is typically far too large to schedule individually, therefore it is desirable to aggregate the blocks into larger collections, and then to schedule these larger collections. To proceed with this aggregation, the ore blocks are clustered 605 (these are typically located towards the bottom of the pit. In one preferred form, those blocks with negative value, which are taken to be waste, are not clustered}. The ore blocks are clustered spatially (using their x, y, z coordinates) and in terms of their grade or value. A balance is struck between having spatially compact clusters, and clusters with similar grade or value within them. These clusters will form the kernels of the atoms of aggregation.
From each cluster, an (imaginary) inverted cone is formed, by propagating upwards using the precedence arcs. This Inverted cone represents the minimal amount of material that must be excavated before the entire cluster can be extracted. Ideally, for every cluster, there is an inverted cone. Typically, these cones will intersect. Each of these intersections (including the trivial intersections of a cone intersecting only itself) will form an atom of aggregation, which is call a clump. Clumps are created, represented by 608.
The number of clumps produced is now far smaller than the original number of blocks. Precedence arcs between clumps are induced by the precedence arcs between the Individual blocks. An extraction ordering of the clumps that is feasible according to these precedence arcs will automatically respect minimum slope constraints. It is feasible to schedule these clumps to find a substantially NPV maximal, clump schedule 607 that satisfies all of the mining oC and processing constraints.
Now that there is a schedule of clumps 607, this can be turned into a schedule of individual blocks. One method is 10 consider all of those clumps that are begun in a calendar year one, and to excavate these block by block starting from the uppermost level, proceeding level by level to the lowermost level. Other methods are disclosed in Section 6 of this specification. Having produced this block ordering, the next step may be to optionally Polish 608 the block ordering to further improve the NPV.
In a more complex case, the step of polish 608, can be bypassed. If itis desirable, however, palishing can be performed to improve the value of the block sequence. 1.3 Balanced NPV optimal / mineable pushback design from block ordering
From this block ordering, we can produce pushbacks, via pushback design 609. Advantageously, the present invention enables the creation of pushbacks that allow for NPV optimal mining schedules. A pushback is a large section of a pit in which trucks and shovels will be concentrated to dig, sometimes for a period of time, such as for one or more years. The block ordering gives us a guide as to where one should begin and end mining. In essence, the block ordering is an optimal way to dig up the pit. However, often this block ordering is not feasible because the ordering suggested is too spatially fragmented. in an aspect of invention, the block ordering is aggregated so that large, connected portions of the pits are obtained (pushbacks). Then a secondary clustering of the ore blocks can be undertaken 610. This time, the clustering is spatial (x, y, z) and has an additional 4th coordinate, which represents the block extraction time ordering.
The emphasis of the 4th coordinate of time may be increased and decreased.
Decreasing the emphasis produces clusters that are spatially compact, but ignore the optimal extraction sequence. Increasing the emphasis of the 4" coordinate produces clusters that are more spatially fragmented but follow the optimal extraction sequence more closely.
Once the clusters have been selected (and ordered in time), inverted cones are propagated upwards in time order. That is, the earliest cluster (in time) is propagated upwards to form an inverted cone. Next, the second earliest N cluster is propagated upwards. Any blocks that are already assigned to the first cone are not included in the second cone and any subsequent cones. Likewise, any blocks assigned to the second cone are not included in any subsequent cones. These propagated cones or parts of cones form the pushbacks 611. This secondary clustering, propagation, and NPV valuation Is relatively rapid, and the intention is that the user would select an emphasis for the 4th coordinate of tims, perform the propagation and valuation, and view the pushbacks for mineability 612. A balance between mineability and NPV can be accessed 613, and if necessary the pushback design steps can be repeated, path 614. For example, if mineability is too fragmented, the emphasis of the 4th coordinate would be reduced. If the NPV from the valuation is too low, the emphasis of the 4th coordinate would be increased.
Once a pushback design has been selected, a minimum mining width routine 615 is run on the pushback design to ensure that a minimum mining width
Is maintained between the pushbacks and themselves, and the pushbacks and the boundary of the pit. An example in the open literature is “The effect of minimum mining width on NPV” by Christopher Wharton & Jeff Whittle, “Optimizing with Whittle” Conference, Perth, 1997. 1.4 Further valuation
A more sophisticated valuation method 616 is possible at this final stage that balances mining and processing constraints, and additionally could take into account stockplling optlons, such as blending and supply chain determination and [or evaluation. 2 Initial identification of clusters
It has been found that the number of blocks in a block model is typically far too large to schedule individually, therefore in accordance with one related aspect of invention, the blocks are aggregated into larger collections. These larger collections are then preferably scheduled. Scheduling means assigning a clump io be excavated in a particular period or periods.
To proceed with the aggregation, a number of ore blocks are clustered.
Ore blocks are identified as different from waste material. The waste material is to be removed to reach the ore blocks. The ore blocks may contain substantially only ore of a desirably quality or quantity and / or be combined with other material - or even waste material. The ore blocks are typically located towards the bottom of the pit, but may be located any where in the pit. In accordance with a preferred aspect of the present invention, the ore blocks which are considered to be waste are given a negative value, and the ore blocks are not clustered with a negative value. It is considerad that thoss blocks with a positive value, present themselves as possible targets for the staging of the open pit mine. This approach Is built around targeting those blacks of value, namely those blocks with positive value.
Waste blocks with a negative value are not considered targets and are thersfore this aspect of invention does not cluster those targets. The ore blocks are clustered spatially (using their x, y, z coordinates) and in terms of their grade or value. Preferably, limits or predetermined criteria are used in deciding the clusters. For example, what is the spatial limit to be applied to a given cluster of blocks? Are blocks spaced 10 meters or 100 meters apart considered one cluster? These criteria may be varied depending on the particular mine, design and environment. For example, Figure 7 illustrates schematically an ore body 701. Within the ore body are a number of blocks 702, 703, 704 and 705. (The ore body has many blocks, but the description will only refer to a limited number for simplicity) Each block 702, 703, 704 and 705 has its own individual Xx, y, z coordinates. If an aggregation is to be formed, the coordinates of blocks 702, 703, 704 and 705 can be analysed according to a predetermined criteria. If the criteria is only distance, for example, then blocks 702, 703 and 704 are situated closer than block 705. The aggregation may be thus formed by blocks 702, 703 and 704. However, if, in accordance with this aspect of invention, another criteria is also used, such as grade or value, blocks 702, 703 and 705 may be considered an aggregation as defined by line 706, even though block 704 is situated closer to blocks 702 and 703. A balance is struck bstween having spatially compact clusters, and clusters with similar grade or value within them. These clusters will form the kernels of the atoms of aggregation. It is important that there is control over spatial compactness versus the grade/value similarity. If the clusters are too spatially separated, the inverted cone that we will ultimately propagate up from the cluster (as will be described below) will be too wide and contain superfluous stripping. If the clusters internally contain too much grade or value variation, there will be dilution of value. It is preferable for the clusters to B substantially sharply identify regions of high grade and low-grade separately, while maintaining a spatial compactness of the clusters. Such clusters have been found to produce high-quality aggregations.
Furthermore, where a relatively large body of ore Is encountered, the ore body may be divided into a relatively large number of blocks. Each block may have substantially the same or a different ore grade or value. A relatively large number of blocks will have spatial difference, which may be used io define aggregates and clumps in accordance with the disclosure above. The ore body, in this manner may be broken up into separate regions, from which individual : cones can be defined and propagated. 3 Propagation of clusters and formation of clumps
From each cluster, an inverted cone (imaginary) is formed. A cone is referred to as a manner of explaining visually to -the reader what occurs.
Although the collection of blocks forming the cone does look like a discretised cone to the human eye. In a practical embodiment, this step would be simulated mathematically by computer. Each cone is preferably a minimal cone, that is, not over sized. This cone Is represented schematically or mathematically, but for the ~ 15 purposes of explanation it is helpful to think of an inverted cone propagating upward of the aggregation. The inverted cone can be propagated upwards of the atom of aggregation using the precedence arcs. Most mine optimisation software packages use the idea of precedence arcs. The cone is preferably three dimensional. The inverted cone represents the minimal amount of material that must be excavated before the entire cluster can be extracted. In accordance with a preferred form of this aspect of invention, every cluster has a corresponding inverted cons.
Typically, these cones will Intersect another cone propagating upwardly from an adjacent aggregation. Each intersection (including the trivial intersactions of a cone Intersecting only itself) will form an atom of aggregation, which is call a ‘clump’, in accordance with this aspect. Precedence arcs batween clumps are induced by the precedence arcs between the Individual blocks.
These precedence arcs are Important for identifying which extraction ordering of clumps are physically feasible and which are not. Extraction orderings must be consistent with the precedence arcs. This means that if block/clump A points to - block/clump B, then block/clump B must be excavated earlier than block/ciump A.
With reference to Figure 8, illustrating a pit 801, In which there are ore bodies 802, 803, and 804. Having identified the important "ore targets” in the stage of initial identification of clusters, as described above, the procedure of propagation and formation of clumps goes on to produce mini pits (clumps) that are the most efficient ways access these "ore targets”. The clumps are the regions formed by an intersection of the cones, as well as the remainder of cones once the intersected areas are removed. In accordance with the embodiment ) aspect, intersected areas must be removed before any others, eg. 814 must be dug up before either 805 or 808, In Figure 8. In accordance with the description above, cones 805, 806 and 807 are propagated (for the purposes of illustration) from ore bodies to be extracted. The cones are formed by precedence arcs 808, 809, 810, 811, 812 and 813. In Figure 8, for example, clumps are designated regions 814 and 815. Other clumps are also designated by what is left of the inverted cones 805, 806 and 807 when 814 and 815 have been removed. The clump area is the area within the cone. The overlaps, which are the intersections of the cones, are used to allow the excavation of the inverted cones in any particular order. The collection of clumps has three important properties. Firstly, the clumps allow access to the all targets as quickly as possible (minimality), and secondly the clumps allow many possible orders of access to the identified ore targets (flexibility). Thirdly, because cones are used, an extraction ordering of the clumps that is feasible according to the precedence arcs will automatically respect and accommodate minimum slope constraints. Thus, the slope constraints are automatically built into this aspect of invention. 4 Splitting of waste and ore In clumps
Once the initial clumps have been formed, a search is performed from the lowest level of the clump upwards. The highest level at which ore is contained in the clump is identified; everything above this level is considered to be waste. The option is given to split the clump into two pieces; the upper piece contains waste, and the lower piece contains a mixture of waste and ore. Figure 9 illustrates a pit 801, in which there is an ore body 962. Fram the ore body, precedence arcs 903 and 904 define a cone propagating upward. in accordance with this aspect of
Invention, line 905 Is identified as the highest level of the clump 902. Then 906 can designate ore, and 907 can designate waste. This splitting of waste from ore designations is considered to allow for a more accurate valuation of the clump.
Many techniques assume that the value within a clump is uniformly distributed,
however, in practice this is often not the case. By splitting the clump into two pleces, one with pure waste and the other with mostly ore, the assumption of homogeneity is more likely to be accurate. More sophisticated splitting based on finer divisions of value or grade are also possible in accordance with predetermined criteria, which can be set from time to time or in accordance with a particular pit design or location. 5 Aggregation of blocks into clumps: high-level ideas
The feature of ‘clumping blocks together’ may be viewed for the purpose of arithmetic simplicity where the number of blocks are too large. The number of clumps produced is far smaller than the original number of blocks. This allows a . mixed integer optimisation engine to be used, otherwise the use of mixed integer engines would be considered not feasible. For example, Cplex by 1LOG may be used. This aspect has beneficial application to the invention disclosed in pending provisional patent application no. 2002951892, titled “Mining Process and Design” filed 10 October 2002 by the present applicant, and which is herein incorporated by reference. This aspect can be used to reduce problem and calculation size for other methods (such as disclosed in the co-pending application above).
The number of clumps produced is far smaller than the original number of blocks. This allows a mixed integer optimisation engine to be used. The advantage of such an engine is that a truly optimal (in terms of maximising NPV) schedule of clumps may be found in a (considered) feasible time. Moreover this optimal schedule satisfies mining and processing constraints. Allowing for mining and processing constraints, the ability to find truly optimal solutions represents a significant advance aver currently available commercial software.
The quality of the solution will depend on the quality of the clumps that are input to the optimisation engine. The selection procedures to identify high quality clumps have been outlined in the sections above.
Some commercial software, as noted in the background section of this specification, do use mixed integer programming engines, however, the method of aggregating blocks is different either in method, or In application, and we — believe of lower-quality. For example, it is considered that ‘ECSI Maximiser uses a form of integer optimisation in their pushback design, and restricts the time window for each block, but the optimisation is local in time, and it's problem formulation is considered too large to optimise globally over the life of a mine. in contrast, in accordance with the present invention, a global optimisation over the entire life of mine is performed by allowing clumps to be taken at any time from : start of mine life to end of mine life. ‘MineMax’ may be used to find rudimentary
S optimal block sequencing with a mixed integer programming engine, however it is : considered that it's method of aggregation does not respect slopes as is required in many situations. ‘MineMax’ also optimises locally in time, and not globally. In use, there is a large huge number of variables, and the user must therefore resort to subdividing the pit to perform separate optimisations, and thus the aptimisation
Is not global aver the entire pit. The present invention is global in both space and time. 6 Determination of a block ordering from a clump ordering
Now that there is a schedule of clumps, it is desirable to turn this into a schedule of individual blocks. One method is to consider all of those clumps that are begun in year one, and to excavate these block by block starting from the uppermost level, proceeding level by level to the lowermost level. One then moves on to year two, and considers all of those clumps that are begun in year two, excavating all of the ‘blocks contained in those clumps level by level from the top level through to the bottom level. And 50 on, until the end of the mine life.
Typically, some clumps may be extracted over a period of several years.
This method just described Is not as accurate as may be required for some situations, because the block ordering assumes that the entire clump is removed without stopping, once it is begun. Another method is to consider the fraction of the clump that is taken in each year. This methad begins with year one, and extracts the blocks in such a way that the correct fractions of each clump for year one are taken in approximately year one. The Integer programming engine assigns a fraction of each clump to be excavated in each period/year. This fraction may also be zero. This assignment of clumps to years or periods must be turned into a sequence of blocks. This may be done as follows. If half of the _ clump A is taken in year one, and one third of ciump B is taken in year one, and all other fractions of clumps in year one are zero, the blocks representing the upper half of clump A and the blocks representing the upper one-third of clump B are joined together. This union of blocks is then ordered from the uppermost bench to the lowermost bench and forms the beginning of the blocks sequence (because we are dealing with year one). One then moves on to year two and repeats the procedure, concatenating the blocks with those already In the sequence.
Having produced this block ordering, block ordering may be in a position to be optionally Polished to further improve the NPV. The step of Polishing is similar to the method disclosed in co-pending application 2002951892 (described above, and incorporated hereln by reference) but the starting condition is different.
Rather than best value to lowest value, as is disclosed in the co-pending application, in the present aspect, the start is with the block sequence obtained from the clump schedule. 7 Second Identification of clusters for pushback design 74 Fuzzy clustering; alternative 1 (spaceitime clustering of block sequence)
From this block ordering, we must produce pushbacks. This is the ultimate goal of KiumpKing - to produce pushbacks that allow for NPV optimal mining schedules. A pushback is a large section of a pit in which trucks and shovels will be concentrated for one or more years to dig. The block ordering gives us a guide as to where one should begin and end mining. In principle, the block ordering is the optimal way to dig up the pit. However, it is not feasible, because the ordering is too spatially fragmented. It is desirable to aggregate the block . ordering so that large, connected portions of the pits are obiained (pushbacks). A secondary clustering of the ore blocks is undertaken. This time, clustering is spatially (x, y, Z) and as a 4th coordinate, which is used for the block extraction time or ordering. The emphasis of the 4th coordinate of time may be increased or decreased. Decreasing the emphasis produces clusters that. are spatially compact, but tend to ignore the optimal extraction sequence. Increasing the emphasis produces clusters that are more spatially fragmented but follow the optimal extraction sequence more closely.
Once the clusters have been selected, they may be ordered In time. The B clusters are selected based on a known algorithm of fuzzy clustering, such as JC
Bezdek, RH Hathaway, MJ Sabin, WT Tucker. “Convergence Theory for Fuzzy c- . means: Counterexamples and Repairs”. IEEE Trans. Systems, Man, and
Cybematics 17 (1987) pp 873-877. Fuzzy clustering is a clustering routine that triss to minimise distances of data points from a cluster centre. In this Inventive aspect, the cluster uses a four-dimensional space; (x, vy, 2, v), where x, y and z give spatial coordinates or references, and ‘v’ is a variable for any one or a
S combination of time, value, grade, ore type, time or a period of time, or any other desirable factor or attribute. Other factors to control are cluster size (in terms of ore mass, rock mass, rock volume, $value, average grade, homogeneity of grade/value), and cluster shape {In terms of irregularity of boundary, spherical- ness, and connectivity). In one specific embodiment, 'v' represents ore type. In another embodiment, clusters may be ordered in time by accounting for 'v' as representing clusters according to their time centres.
There is also the alternative embodiment of controling the sizes of the clusters and therefore the sizes of the pushbacks. "Size®™ may mean rock tonnage, are tonnage, total value, among other things. In this aspect, there is provided a fuzzy clustering algorithm or method, which in operation serves to, where if a pushback Is to begin, its corresponding cluster may be reduced in size by reassigning blocks according to thelr probability of belonging to other clusters.
There Is also another embodiment, where there is an algorithm or method that is a farm of ‘crisp’, as opposed to fuzzy, clustering, specially tailored for the particular type of size control and time ordering that are found in mining applications. This ‘crisp’ clustering is based on a method of slowly growing clusters while continually shuffling the blocks between clusters to improve cluster quality. 7.2 Fuzzy clustering; alteinative 2 (Propagation of clusters)
Having disclosed clustering, above, another related aspect of invention is to then propagate these clusters in a time ordered way without using intersections, to produce the pushbacks.
Referring to Figure 10, a mine site 1001 is schematically represented, in which there is an ore body of 3 sections, 1002, 1003, and 1004.
Inverted cones are then propagated upwards in a time order, as - represented in Figure 10, by lines 1005 and 1006 for cone 1. That is, the earliest cluster (in time) is propagated upwards to form an inverted cone. Next, the second earliest cluster is propagated upwards, as represented in Figure 10 by lines 1007 and 1008 (dotted) for cone 2, and lines 1009 and 1010 (dotted) for cone 3. Any blocks that are already assigned to the first cone are not included in the second cone. .This is represented in Figure 10 by the area between lines 1008 and 1005. This area remains a part of cone 1 according to this inventive aspect. Again, in Figure 10, the area between lines 1010 and 1007 remains a part of cone 2, and not any subsequent cone. This method is applied to any subsequent cones. Likewise, any blocks assigned to the second cone are not included in any subsequent cones. These propagated cones or parts of cones form the pushbacks. 7.3 Fuzzy Clustering; alternative 3 (Feedback loop of pushback design)
In this related aspect, there is a process loop of clustering, propagating to find pushbacks, valuing relatively quickly, and then feeding this information back ‘Into the choice of clusiering parameters.
This secondary clustering, propagation, and NPV valuation is relatively rapid, and the intention is that there would be an iterative evaluation of the result, either by computer or user , and accordingly the emphasis for the 4th coordinate can be selected, the propagation and valuation can be considered and performed, and the pushbacks for mineability can also be considered and reviewed. I the result is considered too fragmented, the emphasis of the 4th coordinate may be reduced. If the NPV from the valuation is too low, the emphasis of the 4th coordinate may be increased.
Referring to Figure 11a, there is illustrated in plan view a two dimensional slice of a mine site. In the example there are 15 blocks, but the number of blocks may be any number. [n this example, blocks have been numbered to correspand with extraction time, where 1 is earliest extraction, and 15 is latest extraction time.
In the example illustrated, the numbers indicate relatively optimal extraction ordering.
In accordance with the aspect disclosed above, Figure 110 illustrates an example of the result of clustering where there is a relatively high fudge factor and relatively high emphasis on time. Cluster number 1 is seen to be fragmented, has a relatively high NPV but Is not considered mineable.
In accordance with the aspect disclosed above, Figure 11c illustrates an example of the result of clustering where there is a lower emphasis on time, as compared to Figure 11b. The result fliustrated is that both clusters number one and two are connected, and ‘rounded’, and although they have a slightly lower ] NPV, the clusters are considered mineable. 8. Aggregation of Precedence Constraints
An approach in accordance with a first aspect of invention is to aggregate the precedence constraints as follows: max. », vx, t 5.t. nx, s > x; 16) x, e€{0,1} Vi where n, =| P(i)] ceveeen..@quation 3
In this first aspect approach, the number of constraints is reduced to one for every block below the surface (there are no precedence constraints for the blocks on the top bench of the pit). In this case each constraint enforces the rule that a block can only be extracted if all of its predecessor blocks are extracted.
However, the total unimodularity property of the exact (disaggregated) formulation is not preserved In this first approach formulation. Hence, the integrality constraints on the decision variables must be enforced. Equation 3 manifests therefore as an integer program, and must be solved using the method of branch- and-bound, rather than the Simplex method. This solution method takes a relatively long time in terms of computation time and can also require a relatively large amount of memory for storage of the decision tree. In particular, obtaining the truly optimal solution (as opposed to a solution within a specified percentage of the optimal solution) may take a relatively long time.
When the aggregated formulation (equation 3) is LP-relaxed and solved In
CPLEX, the decision variables may take fractional values, and the outcome is _ expressed in equation 4 following:
max PRES st. . nx < > x,
FP) 0<x isglIVi . where n, =| P(i)] equation 4 : Consider the case of a relatively small first example of a mine (16,049 blocks) that is provided as an example with the Whitlle software package (by
Whittle Pty Ltd, www,whittle.com.au ). Figure 12 shows the view from above of a comparison of the optimal solutions found by the exact formulation (equation 2) and the LP relaxation of the aggregated formulation (equation 4). The blocks 10 are those that are set to 1 by both the exact formulation {equation 2) and the aggregated formulation (equation 3). The blocks 11 around the outside of this pit are those blocks which are included (set to 1) in the uitimate pit found by the exact formulation {equation 2), but are not included (sat to 0) in the solution found by the LP relaxation of the aggregated formulation (equation 4). It is evident that there are a number of blocks that are included in the true ultimate pit that are not included by the LP relaxation of the aggregated formulation (equation 4). The blocks 12 are waste.
A comparison of a vertical cross-section of the pit design using the exact formulation (equation 2) and the LP relaxation of the aggregated formulation (equation 4) for this first mine example is illustrated in Figure 13 when compared with Figure 14.
Figure 13 shows a plane through the example pit from the view of the solution using the exact formulation (equation 2). The area 20 is the ultimate pit and the area 21 Is waste. Referring to Table 1, below, the total value of this pit is found to be $1.43B85E+09, and CPLEX requires 29.042 seconds to obtain this solution. —
Figure 14 shows the equivalent view when the LP relaxation of the aggregated formulation (equation 4) for the ultimate pit Is used. The area 20 is blocks set to 1, area 21 is waste (blocks set to 0) and area 22 is material which may be further interrogated in ardar to decide whether It Is Included (or not) in the ultimate pit (set to a value between 0 and 1). The total value of this pit is found to } be $1.54268E+09, and found in a CPU time.of 0.992 seconds. Note that the solution of the aggregated formulation (equation 3) (where integrality constraints i 5 are imposed on the decision variables) glves a total value of the ultimate pit to be $1.43591E+09 (using a branch-and-bound stopping criteria of 1% from optimal), which Is similar to the value as that given by equation 2, and a CPU time of 1675.18 seconds was required to obtain this solution.
First example mine Total Blocks [16049
A FS EE
Formulation 1 xact LG (equation 2 I EE
I A EE otal Number of Precedence Constraints 264859 —
TotalValue | tasessewo
CPU Time (Seconds) | 20402
No. Blocks in Ultimate Pit HE % of Total Blocks __ sss8
EE A E——
Aggregated LG (equation 3 IE
Total Number of Precedence Constraints — sa077____
Total Value | taseowesoe
CPU Time (Seconds) | _ tersas
No. Blocks in Ultimate Pit oro %of Total Blocks | eos
Final Gap (from optimal oaes% | 0000
Aggregated LG (equation 4 I © lereweton | 1 0 tal Number of Precedence Constraints_|____ 14077 1542806408
PU Time (Seconds ose
No. Blocks In Ultmate Pt | 7ee0 oof Total Blocks | aes peers equation 9, below — constraints otal Number of Precedence Constraints 34819 1]
TotalVaue | 1assesewo9 or6565 a oofTotalBlocks | sess | 0]
Numberofterations | eo ~~ 1 ~~~
Table 1: Summary of results for first mine example. . It is evident that CPLEX, when using this relaxed aggregated formulation for the problem, provides a relatively higher valued ultimate pit to be found, but does so in a relatively shorter time. This relatively higher value results, in part, from a relaxation of the predecessor constraints, thus allowing a fraction of a block to be taken sven when all of its predecessor blocks have not been taken.
By way of illustration of the reason for finding a relatively higher pit value using equation 4, consider the situation shown in Figure 15. The number within each block represents the value assigned to the decision variable (x) for that block by the LP relaxation of the aggregated formulation (equation 4).
In the case illustrated in Figure 15, Blocks 2 and 3 are predecessors of
Block 1. Block 1 is represented by x4, block 2 by x2 and block 3 by x; in the equations below. In the exact formulation (equation 2), the constraints for this 16 situation illustrated are x, £x,
XN SX ......6quation 5
The solution given (x1 = 0.5, x2 = 0, x3 = 1) is infeasible for the exact formulation (equation 2), since * =0.5>x, 50 equation 6
However, in the LP relaxation of the aggregated formulation (equation 4), the relevant constraint is £7; +X ......equation 7
In this case the solution from Figure 15 is considered feasible (since 2 x 25 05=1<=0+1=1). 2x1 <0+1 - 2 ......equation 8
Hence if Blocks 1 and 3 were ore blocks and had positive value, while
Block 2 was a waste block with negative value, the LP relaxation of the aggregated formulation (equation 4) can take all of Block 3 and 0.5 of Block 1 without incurring the penalty of taking the negative valued Block 2. Hence the aggregated formulation (equation 4) can take fractions of positive blocks that ) otherwise would not have been teken in the exact formulation (equation 2). This leads to a solution of greater value than in the disaggregated case. ) 9. Cutting Plane Method
The LP relaxation of the aggregated formulation (equation 4) can be modified to overcome this solution of arlificially greater value. The result is equation 9 below, namely: max > vx, sd. : 10 nx < > x,
JeE(D 0x islivi where », =| P(i) loop over all arcs {ifi— j, and x > x in sofution, then add the constraint x, = } ceewe-...€quation 9
This approach as expressed by equation 9 is considered a second aspect of Invention termed a ‘cutting plane method’. In this second aspect, an initial (reduced) problem is solved to give an upper bound on the optimal value, and then any consiraints from the overall (Master) problem that are violated by this solution are added, and the problem is re-sclved. This is repeated until substantially no constraints from the Master problem are found to be violated. In this second aspect, the linear program for the aggregated formulation (equation 4) is run and a solution, call it £ is obtained. Each element of the vector £ represents the value (possibly fractional) assigned to each block. Within £ there will be instances of pairs of individual blocks where the constraint that the successor block cannot be taken until the entire predecessor block has been ~ taken {from the exact formulation) is violated. For example, In Figure 15, the constraint in the exact formulation that block 1 is assigned an | value of 0.5 and i is assigned a value of 0 xX SX ......equation 10 is violated, since x1 = 0.5 and x2 = 0. . Thus, in the case of Figure 15, i has a value greater than j and the constraint is added and the solution re-run. The result will be the violation posed . § by Figure 15 as far as blacks 1 and 2, will be removed. Some individual block constraints can be added to the LP relaxation of the aggregated formulation - (equation 4) to make it feasible for the ultimate pit problem. It is possible to perform the following iteration.
For each element of ¥, compare its value with that of each of its predecessor blocks in turn. Whenever there Is a situation where the successor block has a greater value than the predecessor block, add the relative single block constraint to the formulation. For example, in the situation from Figure 15, the constraint xX, SX, : will be added to the LP relaxation of the aggregated formulation (equation 4).
After checking the relationship for all pairs of predecessors, re-solve the problem, subject to the aggregated constraints as well as the added single block precedence constraints. Again, the solution may be infeasible, so the process may have to be repeated. This process should be repeated until the step of - checking single block dependencies reveals that substantially no single block precedence relationships are violated. The solution at this point has been found to be the same as the optimal solution, found by solving the exact formulation (equation 2). : It is considered that the number of constraints needed to obtain the solution using this second aspect approach is significantly less than the number used in the disaggregated formulation. Since the initial aggregated solution gives a reasonable approximation to the ultimate pit, it has been found that only a small percentage of the total number of single block precedence constraints for the problem should need to be added to the formulation. In this way, the _ computational requirement in terms of memory (storage and manipulation of the constraint matrix) to find the optimal solution should be significantly reduced.
However, the cost of this approach is that the process of checking and identification of viclated constraints will require more time than the prior art method of equation 2. When equation 9 is applied to the first mine example ) referred to above, this second approach found the total value of the pit to be : $1.43885E+09, the same as the solution to the problem using the disaggregated formulation (equation 2). The computation time required to achieve this second approach was 976.565 seconds.
A brief comparison of these two methods for the ultimate pit problem at the first example mine is given in Table 1, above. . 10. Aggregation — Cutting Plane and added Blocks and Arc Constraints
It is evident that the trade off between the prior art approach and the approaches of the first and second aspects is time against memory, as illustrated in Table 1, above). The exact formulation (equation 2) finds the optimal solution in 29.402 seconds, while the cutting plane formulation (equation 9) takes 976.565 seconds to find the optimal solution. This is due, in part, to the fact that the cutting plane formulation re-solves a large LP a number of times in the process of solving the problem. [n addition, the process of searching through and checking the entire arcs file (which Is completed as a part of each iteration) takes a significant amount of time. However, the exact formulation {equation 2) solves a model with 264,859 precedence constraints (requiring a significant amount of memory), compared with 34,819 precedence constraints in the cutting plane formulation (equation 5). This is a decrease of 87%. It is expected that the number of constraints in the model is proportional to the memory required to store and solve the problem, in particular, to perform the inversion on the final constraint matrix once the optimal solution has been found. Thus, advantageously, a solution of the cutting plane formulation (equation 8) may be possible in cases where CPLEX runs out of memory when trying to solve the exact formulation (equation 2).
In a second example mine, which has 38,612 blocks, the same approach was taken to that above, with similar results, as shown in Table 2.
Example Mine 2 Total Blocks 38612
SR i RE
Exact LG (equation 2)
Total Number of Precedence
Constraints 1045428
Total Value 1.870640+009 ]
CPU Time (Seconds) 223.762 ]
No. Blocks in Ultimate Pit 33339 ]
Aggregated LG (Cutting
Plane) (equation 9) (LP relaxation + add arc o single block constraints)
Total Number of Precedence
Constraints 159832
Total Value 1.87064E+09 ]
CPU Time (Seconds) 12354.3 ]
Table 2: Summary of results for second mine example. —
In particular, referring to Table 2 above, the exact formulation (equation 2) contains 1,045,428 constraints, while the final model following implementation of . the cutting plane algorithm (equafion 9) requires only 159,832 constraints.
However, the cutting plane method (equation 9) takes 12,354.3 seconds to find the solution, while the exact formulation {equation 2) requires 223.762 seconds of
CPU time. :
Further testing of the alternative mixed integer program approaches to the pit design was carried out on a third mine example, as detailed in Table 3 below.
The block model! for the third mine example contains 198,917 blocks.
Initially, the exact formulation {equation 2) was trailed. This resulted in
CPLEX attempting lo solve a linear program with 3,526,057 single block : constraints. The size of thls constraint matrix caused CPLEX to run out of memory when trying to apply the dual simplex algorithm to solve the problem.
Thus, the exact solution to the pit design In the case of this third mine example is unable to be determined by this approach.
The aggregate formulation (equation 3) was next trailed. This resulted in 188,082 constraints, a value of $3.34125E+09, and a CPU time of 33298.5 seconds.
The next trail was to run the LP relaxation of the aggregated formulation (equation 4). It is expected that the solution to this problem will give an upper bound on the optimal value of the ultimate pit, as was described above. This is due to the fact that CPLEX includes fractions of blocks without necessarily taking their entire precedence set. In this trail, the model had 188,082 constraints, The optimal solution was found to have a value of $3.40296E+09, and this was found in 12.989 seconds of CPU time.
otal Blocks
Exact LG (equation 2 1] otal Number of Precedence Constraints | 352607 . TotalValve | 1]
CPU Time (Seconds omtofmemory o. Blocks in Ultimate Pit I of Total Blacks I
Aggregated LG (equation 3 I
EE BE EE otal Number of Precedence Constraints 188082 ] 334128E+09
CPU Time (Seconds | 33085 | 0]
No. Blocks in Ultimate Pit 97221 ] % of Total Blocks 4.888 {
Final Gap (from optimal 0.99% 1
Aggregated LG (equation 4 oo] (LPretaxationy | 01 otal Number of Precedence Constraints | 1eso82 | ~~ otal Value 340206E+09
CPU Time (Seconds 12.989 ]
No. Blocks in Ultimate Pit ews22 5% of Total Blocks _4e01 | 0]
EE A
I EE
Rr |] (equation 9 constraints otal Number of Precedence Constraints 285598 I [Total Value 337223+09
PU Time (Seconds 19703.8 1]
No. Blocks in Uttimate Pit _ — oesas | 0] % of Total Blocks 49.69 ]
Number of iterations EE —
Table 3: Summary of results for third mine example.
The cutting plane formulation (equation 9) was alsa trailed on this example third mine. This is the method where the solution to the LP relaxation of the } aggregated formulation is used as a starting solution, and then violated single block constraints are added to the model and then again resoived. This process is repeated until no more single block constraints are violated, and thus the solution is similar to that for the exact formulation. The solution to this equation 9 is considered to be the correct solution to the problem. When equation 9 was run, it was found that CPLEX was able to handle the size of the problem, and the exact ultimate pit was found. The solution contained 285,538 constraints, & : reduction of 92% on the exact formulation. The optimal value of the pit design was found to be $3.37223E+09, and the CPU time required to find this solution was 19703.8 seconds.
Thus the cutting plane algorithm (equation 9) has been found to provide an improved solution within the memory limits of a practical implementation of the present invention, using computers and / or computer madelling, where the exact formulation (equation 2) could not. Again, the saving in memory is offset by a longer computation time.
As in the case of the first mine example, a comparison of a vertical cross- section of the solution to the ultimate pit problem using the cutting plane formulation and the LP relaxation of the aggregated formulation for the third mine example is illustrated in the Figures. Figures 16 and 18 show a plane view through the pit using the cutting plane formulation (equation 9). The area 20 is the ultimate pit and the area 21 is waste. Figures 17 and 19, on the other hand, show the same view, but for the LP relaxation of the aggregated (equation 4). Again, areas 20 are the pit and areas 21 are waste. Again, it is evident that the LP relaxation of the aggregated (equation 4) takes fractions of blacks that are infeasible for the exact formulation.
This result is considered to confirm that solution of the cutting plane - formulation (equation 9) may be possible in cases where CPLEX runs out of memory when trying to solve the exact formulation (equation 2). 7
A summary of the results for the third mine example is found in Table 3.
11. Variations On The Cutting Plane Method 11.1 First variation
Since it was found that adding all violated constraints at once causes additional loading on the cutting plane approach (equation 9), due to the very large number of constraints added by the first iteration, one variation of the cutting plana method is to add the constraints incrementally. Initially, the effect of : adding the most violated constraints first, and then re-solving the formulation was investigated. This method was thoroughly tested on the first mine example. The approach taken was as follows. At each lteratlon of the method, a lower bound on the size of the vialation of the single block constraint was specified (e.g. 0.5, 08, ...). For example, Figure 15 illustrates violations for each block. In this example Figure 15, the violation = x; — xj, and so the ‘size’ of the violation is 0.5 ~ 0 = 0.5. Constraints that were violated by an amount greater than this tolerance were added to the formulation, and the problem was re-solved. However, using this approach the optimisation process compisted before the optimal solution was found. This occurs because this method of adding constraints does not identify and add all single block constraints that are violated, only those that are violated by more than a certain amount. In this way, not all of the necessary single block constraints are added fo the formulation, and the truly optimal solution is not reached. To alleviate this problem, violation(s) greater than a selected lower bound is added to at least the first iteration. This approach enables an optimal solution is still obtained. 11.2 Second variation
Another approach is to add the most violated constraints, but to decrease the amount of violation required at each iteration until a certain number of constraints have been added. For example, it may be designated that a minimum of 5000 constraints should be added at each iteration. Say the initial violation parameter is set to 0.6 (that is, only single block constraints that are violated by 0.6 or more are added to the formulation). it may be the case that 1200 B constraints are added. Then, before re-solving the formulation, the violation parameter couid be decreased to 0.5. This may result in a further 3000 constraints being added to the model. Since there are still less than 5000 . constraints added, the violation parameter is further decreased to 0.4, and more single block constraints are added. This may result in 2000 constraints being added to the farmulation, and the problem is now re-solved since the minimum of 5000 constraints has been reached. The process is then repeated until the optimal solution is obtained. 11.3 Third variation
Altemnatively, the tolerance could be reduced on a smaller incremental level (say 0.01 at a time instead of 0.1) in an attempt to reduce the size of the overshoot on the number of constraints added compared with the prescribed minimum number of constraints. 11.4 Fourth variation
A further alternative is simply to add a specified number of constraints to the model before the formulation is re-solved. In any approach where a minimum number of constraints are added, the determination of the appropriate number of constraints fo add at each iteration is a non-frivial matter. This element of the 16 problem may itself require optimisation. It is expected that the maximum size of the problem that is able to be stored in memory and handled by CPLEX will affect this value. Consideration of this fact may allow a test to be built in to the program for solving the ultimate pit problem. The form of the test procedure could proceed as follows. If the size of the constraint matrix following the first iteration is less than the maximum size able to be solved by CPLEX, (with a margin to allow more constraints fo be added in subsequent iterations based on the general proportion of constraints added after the initial loop — it appears that approximately 90% of the constraints that are required are added in the first loop), take the path of adding all violated constraints. If the size of the constraint matrix following the first iteration is greater than the maximum able to be solved, restart the iteration process using one of ths alternative constraint-adding processes described above.
The approaches described above were tested on the first mine example above. In this case, the approach that performed the best was to add single block constraints that were violated by more than 0.6 in the first 5 loops, and in BN subsequent loops, add all violated constraints. This approach found the optimal solution in 2152.24 seconds. This was significantly longer than the standard a1 cutting plane procedure, which required 976.565 seconds (compare with statement below). 11.5 Fifth variation
Another approach for adding constraints incrementally takes advantage of ; 5 the specific geometry of the mine. In this case, a vector containing the z coordinate {or “height®) for each block is stored. Using this information, violated single block constraints are added from the largest z coordinate (corresponding to the top of the pit) down, decreasing by block height, in each loop. The constraint adding process stops either once a specified number of constraints have been added, or after a specified number of z coordinates have been descended. By adding violated single block constraints from the largest z coordinate down, itis hoped that the subsequent optimisation steps will force more single block constraints from lower in the pit to be satisfled before they need to be explicitly added to the formulation in a cutting plane iteration. That Is, once decisions regarding the uppermost benches of the pit have been made, the precedence constraints within the formulation could force these decisions to propagate down the pit. Subsequently, less single block constraints may need to be added through the cutting plane iterations pefore the problem is solved to optimality.
This approach was particularly effective in the case of the third mine example. The optimal solution to the problem was found in 2664.11 seconds : when constraints were added from the top z coordinate down in each iteration, with ten z coordinates descended in each iteration. This compares very favourably with the standard cutting plane formulation, which requires 19,703.8 seconds to find the optimal solution. "While this Invention has been described in connection with specific embodiments thereof, it will be understood that it Is capable of further modification(s). This application is intended to cover any variations uses or adaptations of the invention following in general, the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains and as may be — applied to the essential features hereinbefore set forth.
The present Invention may be embodied In several forms without departing from the spirit of the essential characteristics of the invention, it should be understood that the above described embodiments are not to limit the present invention unless otherwise specified, but rather should be construed broadly within the spirit and scope of the invention as defined in the appended claims.
Various modifications and equivalent arrangements are iniended to be included within the spirit and scope of the invention and appended claims. Therefore, the specific embodiments are to be understood to be illustrative of the many ways in which the principles of the present invention may be practiced. In the following claims, means-plus-function clauses are intended to cover structures as performing the defined function and not only structural equivalents, but also equivalent structures. For example, although a nail and a screw may not be structural equivalents in that a nail employs a cylindrical surface to secure ' wooden parts together, whereas a screw employs a helical surface to secure waoden parts together, in the environment of fastening wooden parts, a nail and a screw are equivalent structures.
Claims (1)
- THE CLAIMS DEFINING THE INVENTION ARE AS FOLLOWS:1. A method of aggregating a plurality of blocks related to a design configuration for extracting material from a particular location, the method including the steps of: detenmining a selected volume of material to be extracted, dividing at least a portion of the selected volume into blocks, forming at [east one cluster, and : propagating from at least one cluster, a cone. 2, A method as claimed in claim 1, further including the step of determining from Intersections of the cones, a clump.3. A method as claimed in claim 1 or 2, wherein the cone is propagated upwards using precedence arcs. 4, A method as claimed in claim 1, 2 or 3, wherein the cone is three dimensional. 8, A method as claimed in any one of claims 1 to 4, wherein the cone is minimal.6. A method as claimed in any one of claims 1 to 5, wherein the cone includes a number of blocks. :7. A method as claimed in any one of claims 2 to6, wherein the precedence arcs relate ta an extraction order of corresponding clump(s) and / or block(s).8. A method of determining slope constraints related to a design configuration for extracting material from a particular location, the method including the steps = of: using precedent arcs of blocks emanating from a selected clump lo. establish, at least in part, clump precedents.9. A computer program product including: a computer usable medium having computer readable program code and computer readable system code embodied on said medium for determining slope constraints related to a design configuration for extracting material from a particular location within a data processing system, said computer program product Including: computer readable code within said computer usable medium for performing the method as claimed in any one of claims 1 10 8.10. Apparatus adapted to aggregating a plurality of blocks related to a design configuration for extracting material from a particular location, the apparatus including: first means for determining a selected volume of material to be extracted, second means for dividing at least a portion of the selected volume into blocks, third means for forming at least one cluster, and fourth means for propagating from at least one cluster, a cone.11. Apparatus adapted to determining slope constraints related to a design configuration for extracting material from a particular location, the apparatus including: - first means for using precedent arcs of blocks emanating from a selected clump to establish, at least in part, clump precedents.12. Apparatus adapted to aggregating a plurality of blocks related to a design configuration for extracting material from a particular location, the apparatus including: processor means adapted to operate in accordance with a predetermined instruction set, said apparatus, in conjunction with said instruction set, being adapted to _ perform the method as claimed in any one of claims 1 to 8.13. Apparatus adapted to deterniining slope constraints related to a design configuration for extracting material from a particular location, the apparatus including: processor means adapted to operate in accordance with a predetermined instruction set, said apparatus, in conjunction with said instruction set, being adapted to perform the method as claimed in any one of claims 1 to 8.14. A method of determining a cluster of material, the method including the steps of: allocating at least a portion of the material between a plurality of blocks, determining a first attribute related to co-ordinates corresponding to each block, assigning the first attribute to each comesponding block, determining a second and { or at least one further attribute(s) related to the plurality of blocks, and aggregating at least two of the plurafity of blocks In accordance with the first attribute and the second attribute.15. A method as claimed in claim 14, wherein the second and / or further atiribute(s) provide for distinction between material considered of value, and material considered to be of lesser or little value.16. A method as claimed in claim 14 or 15, wherein the second and / or further attribute(s) corresponds to a positive value.17. A method as claimed in claim 14 or 16, wherein the second and / or further attribute(s) corresponds to a spatial value.18. A method as claimed in claim 14 or 15, wherein the second and / or further attribute(s) corresponds to grade of the material.19. A method as claimed In any one of claims 14 fo 18, wherein the at least two blocks are aggregated into a larger collection of blocks.20. A method as claimed in any one of claims 14 to 19, further including repeating the steps defined to determination a further cluster from material already aggregated.21. A method of removing material from a predetermined location, the method including the steps of: determining a cluster in accordance with any ane of claims 14 to 20 , and scheduling the collections into ona or more periods.22. A computer program product including: a computer usable medium having computer readable program code and computer readable system cade embodied on said medium for determining slope constraints related to a design configuration for extracting material from a particular location within a data processing system, said computer program product including: computer readable codes within said computer usable medium for performing the method as claimed in any one of claims 14 to 21.23. Apparatus adapted to determining a cluster of material, the apparatus including: first means for allocating at least a portion of the material between a plurality of blocks, second means for determining a first attribute related to co-ordinates corresponding to each block, third means for assigning the first attribute to each corresponding block, fourth means for determining a second and / or at least one further attribute(s) related to the plurality of blocks, and a fifth means for aggregating at least two of the plurality of blocks in accordance with the first atiribute and the second attribute.24. Apparatus adapted to deiermine the removal of material from a predetermined location, the apparatus including: first means for determining a cluster in accordance with any one of claims 14 to 20, and second means for scheduling the collections into one or more periods.25. Apparatus adapted to determining a cluster of material, sald apparatus including: processor means adapted to operate in accordance with a predetermined instructions set, said apparatus, in conjunction with said Instruction set, being adapted to perform the method as claimed in any one of claims 14 to 20.26. Apparatus adapted to determine the removal of material from a predetermined location, said apparatus including: processor means adapted to operate in accordance with a predetermined instructions set, : sald apparatus, in conjunction with said instruction set, being adapted to perform the method as clalmed in any one of claims 14 to 20.27. A method of determining characteristics of a selected portion of material, the method including the steps of: determining the contents of the selected portion of material, and identifying region(s) of material within the selected portion according to at least ona of a plurality of characteristic(s).28. A method as claimed in claim 27, wherein the determination of contents is performed from substantially the lowest level of the selected portion of material, and then upwards.20. A method as claimed in claim 27 or 28, wherein the selected portion of : material is a clump.30. A method as claimed in claim 27, 28 or 29, wherein the characteristics are predetermined.31. A method as claimed in any one of claims 27 to 30, wherein at least two of the characteristics reflect one or more grade(s) or value(s) of material.32. A method as claimed In claim 31, wherein the region(s) distinguish between ore considered of value, waste material and / or impurity.33. A method as claimed in claim 32, wherein the region(s) of ore considered of value includes a portion of waste material and / or impurity.34. A method as claimed in any one of claims 27 to 33, further including the step of re-iterating an evaluation and / or extraction analysis of the selected material based on the identification of regions.35. A mine analysed in accordance with the method as claimed In any one of claims 27 to 34.36. Material extracted from a mine as claimed in claim 35.37. Apparatus adapted to determine characteristics of a selected portion of } material, the apparatus including: first means to determine the contents of the selected portion of material, and means fo Identify region(s) of material within the selected portion according to at least one of a plurality of characteristic(s).38. Apparatus including processor means adapted to operate in accordance with a predetermined instruction set, said apparatus, in conjunction with the instruction set, being adapted to perform the method as claimed in any one of claims 27 to 34.39. A computer program product including: computer usable medium having computer readable program code and computer readable system code embodied on said medium for determining slope constraints related to a design configuration for extracting material from a particular location within a data processing system, said computer program product including: computer readable code within said computer usable medium for performing the method as claimed in any one of claims 27 to 34.40. A method of analysing a selected volume of material, the material being at least partially comprised of a plurality of blocks, the method including the steps of: clumping a number of blocks together, analysing the selected volume of material based on the clumped blocks.41. A method as claimed in claim 40, wherein a mixed integer optimisation engine is used to analyse the selected volume of material.42. A method as claimed in claim 41, wherein further constraints are incorporated into the engine.43. A method as claimed in claim 42, wherein the further constraints are mini ng, processing, marketing capabilities and/or grade constraints.44. A method as claimed in any one of claims 40 to 42, wherein analysis of the selected volume is performed relatively globally both in space and time. =45. A mine analysed In accordance with the method as claimed in any one of i claims 40 to 44.46. Material extracted from a mine as claimed in claim 45.47. A computer program product including: computer usable medium having computer readable program code and computer readable system code embodied on said medium for determining slape constraints related to a design configuration for extracting material from a particular location within a data processing system, said computer program praduct including: computer readable code within said computer usable’ medium for performing the method as claimed in any one of claims 40 to 43.48. Apparatus adapted to analyse a selected volume of material, the material being at least partially comprised of a plurality of blocks, the apparatus including: first means for clumping a number of blocks together, : second means for analysing the selected volume of material based on the clumped blocks. 49, Apparatus including processor means adapted to opsrate in accordance with a predetermined instruction set, said apparatus, in conjunction with the instruction set, being adapted to perform the method as claimed in any one of claims 40 to 44.50. A method of determining a selected group of blocks of a mine pit which are capable of being mined, the method including the step of selecting a first plurality of blocks, and determining a relative value and constraints applicable to the selected first : plurality of blocks in accordance with:max vx,!s.t. nx, > x, . JEP) x, {0,1} Vi where n, =|P(i}| ceeererenn.equation 351. A method of determining a selecied group of blocks of a mine pit which are capable of being mined, the method including the step of selecting a first plurality of blocks, and determining a relative value and constraints applicable to the selected first plurality of blocks in accordance with: max PRE? {sd. nx =< y, x, JaP(1) 0<x isl vi where n, =| P()] eens n@quUation 452. A method of determining a selected group of blocks of a mine pit which are capable of being mined, the method including the step of selecting a first plurality of blocks, and determining a relative value and constraints applicable to the selected first plurality of blocks In accordance with:max Sux,st. nx, < 3 x; JeP() 0<x i<1Vi : where n, =| P(§) loop over all arcs : {ifi— j, and xi > xj in solution, then add the constraint xi =x] } veer.....€Quation 953. A method as claimed in claim 52, further including the step of re-testing for violations.54. A method as claimed in any one of claims 50 to 53, further including the further step of: selecting a second plurality of blocks, determining a relative value and constraints applicable to the selected second plurality of blocks in accordance with any one of equations 3, 4 or 9, and determining whether the first or second plurality of block have a greater value.55. A method as claimed in any one of claims 50 to 53, further including the step of repeating the further steps as claimed in claim 5 untii substantially no constraints from the problem are found to be violated.56. A method as claimed In claim 53, wherein a violation greater than or equal to a selected lower bound is added in the first iteration.57. A method as claimed in claim 56, wherein, on subsequent iteration(s), a constraint with-a-reduced-lower bound is.added.58. A method as claimed in claim 53, wherein a tolerance value Is used to determine the number of constraints.50. A method as claimed in claim 53, wherein the number of constraints added is determined In accordance with:a. If the number of constraints does not exceed a memory limit, then add substantially all the constraints, or b. If the number of constraints exceeds a memory limit, then add some (only) of the constraints.60. A method as claimed in claim 53, wherein the constraint added is in accordance with a z co-ordinate.81. A method as claimed in claim 60, wherein the constraint added Is in accordance with the largest z coordinate first.62. A method of determining a selected group of blocks of a mine pit which are capable of being mined, the method including the step of selecting a first plurality of blocks, and detemmining a relative value and constraints applicable to the selected first plurality of blocks in accordance with: a CPLEX method as defined by equation 2, disclosed herein, and if the CPLEX method times out, redetermining the relative dimension and / or size of a mine pit using the method as claimed in any one of claims 50 to 61.63. Apparatus adapted to determining a selected group of blocks of a mine pit which are capable of being mined, the apparatus including: processor means adapted to operate in accordance with a predetermined instruction set, said apparatus, in conjunction with sald instruction set, being adapted to perform the method as claimed in any one of claims 50 to 62.64. A computer program product including: a computer usable medium having computer readable program code and aN computer readable system code embodied on sald medium for determining the relative dimension and / or size of a mine pit using a data processing system, said : computer program product including:computer readable code within said computer usable medium for determining a selection of blocks of a mine pit which are capable of being mined in accordance with the method as claimed in any one of claims 50 to 62.65. A method as claimed in any one of claims 1 to 8, 14 to 21, 27 to 34, 40 to 44 and 50 to 62 substantially as herein described with reference to the accompanying drawings.66. An apparatus as claimed in any one of claims 12 to 15, 23 to 26, 37, 38, 48, 49 and 63, substantially as herein described with reference to the accompanying drawings.67. A mine design In accordance with the method as claimed in any one of claims 1 to 8, 14 to 21, 27 to 34, 40 to 44 and 50 to 62.68. Material extracted from a mine in accordance with the design as claimed in any one of claims 12 fo 15, 23 to 26, 37, 38, 48, 49 and 63.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| AU2002951891A AU2002951891A0 (en) | 2002-10-09 | 2002-10-09 | System and Method(s) of Mine Planning, Design and Processing |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| ZA200502772B true ZA200502772B (en) | 2006-06-28 |
Family
ID=28679524
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ZA200502772A ZA200502772B (en) | 2002-10-09 | 2006-04-06 | System and method(s) of mine planning, desing and processing |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN1723334B (en) |
| AU (1) | AU2002951891A0 (en) |
| ZA (1) | ZA200502772B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106640081A (en) * | 2016-08-25 | 2017-05-10 | 中国黄金集团内蒙古矿业有限公司 | Strip mine drilling and blasting design, construction and management method |
| CN109458179B (en) * | 2019-01-04 | 2020-05-08 | 天地科技股份有限公司上海分公司 | Underground mechanized continuous mining method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN85107072A (en) * | 1985-09-25 | 1987-04-08 | 卡托维兹煤炭协会 | Be used in the quarry bed process of safeguarding permanent control surface deformable, especially in influenced scope, dig up mine |
| CN1114712A (en) * | 1994-06-04 | 1996-01-10 | 崔正洙 | Mining and ore dressing process new modle for mine |
| CN1165237A (en) * | 1997-04-24 | 1997-11-19 | 张中科 | Boundary go up (down) hill along ampty tunnelling type mining method |
-
2002
- 2002-10-09 AU AU2002951891A patent/AU2002951891A0/en not_active Abandoned
-
2003
- 2003-10-02 CN CN200380105361.2A patent/CN1723334B/en not_active Expired - Fee Related
-
2006
- 2006-04-06 ZA ZA200502772A patent/ZA200502772B/en unknown
Also Published As
| Publication number | Publication date |
|---|---|
| AU2002951891A0 (en) | 2002-10-24 |
| CN1723334A (en) | 2006-01-18 |
| CN1723334B (en) | 2011-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7957941B2 (en) | System and method(s) of mine planning, design and processing | |
| AU2009202456B2 (en) | System and method(s) of blended mine planning, design and processing | |
| Topal et al. | A new algorithm for stope boundary optimization | |
| US20210208305A1 (en) | Mine planning method and system | |
| Yashar et al. | A multi-step approach for block-cave production scheduling optimization | |
| Little et al. | Integrated optimization of underground mine design and scheduling | |
| Manríquez et al. | A simulation–optimization framework for short-term underground mine production scheduling | |
| Khodayari et al. | Mathematical programming applications in block-caving scheduling: a review of models and algorithms | |
| Pourrahimian | Mathematical programming for sequence optimization in block cave mining | |
| US20070185750A1 (en) | Method, apparatus and computer program for scheduling the extraction of a resource and for determining the net present value of an extraction schedule | |
| Smith et al. | Medium-term production scheduling of the Lumwana mining complex | |
| e Faria et al. | Integrated stochastic optimization of stope design and long-term underground mine production scheduling | |
| Grieco et al. | Grade uncertainty in stope design—improving the optimisation process | |
| ZA200502772B (en) | System and method(s) of mine planning, desing and processing | |
| Malaki et al. | An application of mathematical programming and sequential Gaussian simulation for block cave production scheduling | |
| AU2003266820B2 (en) | System and method(s) of mine planning, design and processing | |
| Hall et al. | The role of mine planning in high performance | |
| Afum | Open pit-underground mining options and transitions planning: a mathematical programming framework for optimal resource extraction evaluation | |
| Villa et al. | Multi-lifts selection using scenario simulation and financial metrics | |
| O'Sullivan | An optimization-based decomposition heuristic for solving complex underground mine scheduling problems | |
| Sikora | Choosing the location of the opening cut to expose brown coal deposits–problem solving and decision making with the use of multiple-criteria decision analysis (MCDA) | |
| Quigley | The stochastic optimization of long and short-term mine production schedules incorporating uncertainty in geology and equipment performance | |
| Sotoudeh et al. | A New Mathematical Programming Formulation for Production Scheduling Optimisation of Sublevel Stoping Operations in the Presence of Pre-concentration Systems | |
| Khodayari et al. | A Review of Mathematical Models and Algorithms in Block-Caving Scheduling | |
| Dimitrakopoulos | An Extended Stochastic Optimization Method for Multi-Process Mining Complexes L. Montiel |