AU2005254579B2

AU2005254579B2 - Method, apparatus and computer program for scheduling the extraction of a resource and for determining the net present value of an extraction schedule

Info

Publication number: AU2005254579B2
Application number: AU2005254579A
Authority: AU
Inventors: Merab Menabde
Original assignee: BHP Billiton Innovation Pty Ltd
Current assignee: BHP Billiton Innovation Pty Ltd
Priority date: 2004-06-21
Filing date: 2005-05-27
Publication date: 2010-03-11
Anticipated expiration: 2025-05-27
Also published as: AU2005254579A1

Description

WO 2005/124104 PCT/AU2005/000761 METHOD, APPARATUS AND COMPUTER PROGRAM FOR SCHEDULING THE EXTRACTION OF A RESOURCE AND FOR DETERMINING THE NET PRESENT VALUE OF AN EXTRACTION SCHEDULE 5 Field of the Invention The present invention relates to a method and apparatus for scheduling the extraction of a resource and for determining the net present value of an extraction schedule. Typically the resource is an ore body which is 10 to be mined in an open cut mine. Background of the Invention Typically, the mining of a resource may take place over a period of 15-30 years before the mine is exhausted to a 15 position where further mining is not economic. The process of mine development and design and long-term scheduling is based on spatial interpretation of drillhole data. Thus, drillholes are drilled over the region to be 20 mined and data relating to the grade of the resource is obtained. The drillholes are generally drilled a reasonable distance apart because this procedure is relatively expensive. A block model of the region to be mined is created and typically, the block model may 25 contain between 50,000 to 1,000,000 blocks which are to be scheduled for mining over the period of 15-30 years. A block is that material enclosed by a rectangular prism in the ground and may contain air to a volumetric percentage strictly less than 100%. A block model is a collection of 30 non-intersecting blocks that are usually, but not necessarily, spatially connected and which contain no less than all material considered to have economic value in a mining enterprise. The objective of the scheduling procedure is to find the block extraction sequence which 35 produces the maximum possible net present value (NPV) and obeys a number of constraints. The constraints include: WO 2005/124104 PCT/AU2005/000761 -2 (a) geotechnical slope constraints which are modelled by a set of precedence rules constraining the order of extraction of individual blocks; (b) mining constraints, i.e. total maximum amount of 5 rock which can be mined in one time period (usually 1 year); (c) processing constraints, i.e. maximum amount or ore which can be processed through a given processing plant in one time period; 10 (d) and the market constraints, i.e. the maximum amount of metal, which can be sold in one time period. (e) Any other constraints salient to the practical mining operation including but not restricted to 15 maximum limits in sinking-rate and available ore. A schedule is a period of extraction for each block and a destination for each block (waste, stockpile or process 20 plant. The ore body model which is built up from the drillhole data is a deterministic model created by spatial interpretation of the data using some kind of so-called 25 Kriging procedure. This enables each of the blocks in the model to be assigned a resource grade (i.e. the amount of the resource present in the block). The resource grade information is then used to determine the scheduling of the mining operation, and also whether a particular block 30 is sent for processing to extract the resource, sent to waste, or stockpiled for later processing. Because the drillholes are generally drilled some distance apart, the drillhole data is usually sparse and therefore 35 this introduces inherent errors in the deterministic block model. To some extent, this can be overcome by providing more data by drilling more holes. However, as is WO 2005/124104 PCT/AU2005/000761 -3 explained above, the drilling of the drillholes is expensive and therefore, this is not desirable. Thus, traditionally open pit mine planning is based on the 5 block model which is built up using some kind of interpolation technique such as the Kriging procedure so that a single model is produced. This single model is assumed to be a fair representation of reality and is used for mine design and optimisation. The design process 10 consists of three main steps: (a) finding the block extraction sequence which produces the best net present value whilst satisfying geotechnical slope constraints; (b) designing the practically mineable mine phases 15 (so-called push backs) which are roughly based on the optimal block sequence; and (c) optimising the mining schedule and cut-off grades. 20 The cut-off grade (COG)is defined as the threshold such that the blocks with a grade above it are sent to the processing plant and with a grade below it are treated as waste. It can be constant for the whole life of mine, or can be variable, i.e. dependent on the period of 25 extraction. In practice, the open mine is divided into a number of the mining phases, which are mined bench by bench, each bench being represented by a horizontal layer of blocks within 30 the given mining phase and having the same elevation. A bench within a mining phase is sometimes referred to as a "panel" (one or more layers of blocks). The mining phases can be mined one by one from top to bottom. However, this kind of schedule is usually sub-optimal. Mining several 35 phases simultaneously and applying a variable cut-off grade can produce much better results. There are several commercially available packages which claim to optimise WO 2005/124104 PCT/AU2005/000761 -4 the schedule and cut-off grade using a single block model representation of the resource. However, it is difficult to estimate their effectiveness as the upper theoretical limit on the net present value remains unknown. 5 The standard optimisation technique widely used in many industrial applications is the linear and integer programming (e.g. Padberg, 2003). However, in order for this program to operate satisfactorily, the problem to be 10 solved needs to be formulated as a linear one. Summary of the Invention A first aspect of the invention is concerned with scheduling the extraction of a resource to improve the net 15 present value. The invention may therefore be said to reside in a method of scheduling the extraction of a resource within a medium comprising the steps of: 20 obtaining drillhole grade data relating to the resource contained within the medium; creating a plurality of different block models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models 25 honouring the drillhole grade data; and scheduling extraction of blocks based on a cut-off grade policy to increase expected net present value having regard to all of the block models. 30 Thus, by using a plurality of different block models, a much larger number of possible representations of the presence of the resource in the medium are obtained. Each of those representations is equally likely to occur because each honours the drillhole grade data which have 35 been measured. Thus, by using all of those possibilities, a more accurate representation of the actual resource distribution within the medium can be obtained, and WO 2005/124104 PCT/AU2005/000761 -5 therefore a better indication of whether a region of the medium represented by a block should be forwarded to processing. For example, if the cut-off grade of the resource in a particular period is 0.6%, the method is 5 able to see that a particular block may in some realities be below this cutoff grade (in which case it is sent to waste), and in some realities is above this cutoff grade (in which case it is processed). By using all of the models, such as by averaging different block valuations 10 (processed or wasted) from different realities we more accurately model the process in reality. This is because at the time of extraction extra blast-hole data are available to enable determination of a more accurate quantification of the block grade and hence the 15 selectivity option encoded in the use of multiple conditional simulations can be gainfully exploited in calculating both the schedule and undiscounted block valuation and hence ultimately in obtaining a more accurate determination of the net present value of each 20 block. Thus, a schedule is obtained which gives an increased and more accurate net present value, compared to existing techniques, of a mine containing the resource and the medium. 25 In the preferred embodiment of the invention, the creation of the plurality of different block models is a conditional simulation technique which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and first and 30 second order statistics of the resource represented, respectively, by probability distribution function and variogram. In one embodiment of the invention, the number of blocks 35 in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction.

WO 2005/124104 PCT/AU2005/000761 -6 Instances of block aggregations are assumed to have homogeneous properties within each grade bin. A collection of one or more blocks that are members of a 5 particular aggregate where the grade attribute of each block falls within defined upper and lower limits is called a grade bin. Grade bins-are used to practically apply cut-off grades, viz: all blocks in grade bins above the cut-off grade are sent to process plant. 10 In one embodiment of the invention, the cut-off grade is a fixed cut-off grade. However, in the preferred embodiment of the invention, the cut-off grade is a variable cut-off grade. 15 In the preferred embodiment of the invention, the plurality of block models are used to provide an average of resource grade and block valuations and the average of the resource grade and block valuations is used to produce 20 the scheduled extraction of blocks and the cut-off grade value to increase expected net present value. In one embodiment of the invention, the scheduled extraction of blocks may be optimised by considering the 25 average valuation produced from the conditional simulations. However, in the preferred embodiment of the invention, all of the conditional simulations are optimized 30 simultaneously to produce a single extraction schedule with a variable cut-off grade policy, and wherein all the conditional simulations are valued using this schedule and the valuations are averaged to produce the expected net present value, which is substantially better than the 35 one obtained by optimizing the average grade block model.

WO 2005/124104 PCT/AU2005/000761 -7 In the preferred embodiment of the invention, the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 5 In the preferred embodiment, the step of scheduling using a variable cut-off grade is converted from a non-linear representation to a linear representation, and a mixed integer programming formulation is used for simultaneous 10 optimisation of the extraction schedule and the cut-off grade. Preferably the scheduling is determined by Maximise -- JE V'x, d' (1) (N n=1 i=1 i=1 ,=1 Y U 15 subject to the following constraints: I R,"x :5 R,*,for allt (2) N n=1 i=1 j-1 1 Qxy, Q, , for all t (3) Nn=1 i=1 j=1 20 y, 1 yi,,foralliandt (4) tG I xy,5y ,,for alli (5) r=1 j=1 G T yi, X Xkjr , for all i, t and k cSi (6) j=1 t=1 G 5 j, =1, for all t (7) j=1 xyt : (j, ,for all i,j, and t (8) 25 wherein; T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; 30 G is the number of all possible cutoff grades; R[" is the total rock in the panel i in simulations n. Qn is the total ore in the panel i, simulation n, when WO 2005/124104 PCT/AU2005/000761 -8 mined with the COG j; Voj is the value of the panel i, simulation n, when mined and processed with the COG j; R,' is the maximum mining capacity in period t; 5 Q," is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; dt is the time discount factor; x,, is the fraction of the panel i is extracted with the 10 COG j in period t; y,, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise; 8, is a binary variable controlling the selection of the 15 COG applied in period t. The invention may also be said to reside in an apparatus of scheduling the extraction of a resource within a medium from which drillhole grade data relating to the resource 20 contained within the medium has been obtained comprising: a processor for creating a plurality of different block models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models honouring the drillhole grade data; and 25 a processor for scheduling extraction of blocks based on a cut-off grade policy to increase expected net present value having regard to all of the block models. In the preferred embodiment of the invention, the 30 processor is for creating conditional simulations which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and first and second order statistics of the resource represented, respectively, by probability distribution 35 function and variogram.

WO 2005/124104 PCT/AU2005/000761 -9 In one embodiment of the invention, the processor is also for reducing the number of blocks in each of the block models by aggregating blocks prior to the step of scheduling extraction. 5 In one embodiment of the invention, the cut-off grade is a fixed cut-off grade. However, in the preferred embodiment of the invention, the cut-off grade is a variable cut-off grade. 10 In the preferred embodiment of the invention, the processor is for providing an average of resource grade and block valuations from the plurality of models, and the average of the resource grade and block valuations is used 15 to produce the scheduled extraction of blocks and the cut off grade value to increase expected net present value. In one embodiment of the invention, the processor is for scheduling extraction of blocks optimised by considering 20 the average valuation produced from the conditional simulations. However, in the preferred embodiment of the invention, all of the conditional simulations are optimized 25 simultaneously to produce a single extraction schedule with a variable cut-off grade policy, and wherein all the conditional simulations are valued using this schedule and the valuations are averaged to produce the expected net present value, which is substantially better than the 30 one obtained by optimizing the average grade block model. In the preferred embodiment of the invention, the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent 35 to waste or stockpiled for later processing.

WO 2005/124104 PCT/AU2005/000761 - 10 In the preferred embodiment, the processor is for converting from a non-linear representation to a linear representation when a variable cut-off grade is used, and a mixed integer programming formulation is used for 5 simultaneous optimisation of the extraction schedule and the cut-off grade. Preferably the scheduling is determined by the processor by I N P G T 10 Maximise -- EEZ V.xn, d' (1) N n=1 i=1 j=1 ,=1 subject to the following constraints: 1 N P G _ZEER"xnR, , for allt (2) n=1 i=1 j-1 15 Q x ZQn Q,* , for all t (3) N n=1 i=1 j=1 yi1:! yi,, for all iand t (4) i G E xu,:5yi, ,for all i (5) r=1 j=1 G T Yt 5 E xk, , for all i, t and k c Si (6) j=1 t=1 G SSj, =1, for all t (7) j=1 20 xY,58j,,foralli,j,and t (8) wherein; T is the number of scheduling periods; N is the number of simulations; 25 P is the total number of panels; G is the number of all possible cutoff grades; Ri' is the total rock in the panel i in simulations n. Q is the total ore in the panel i, simulation n, when mined with the COG j; 30 V," is the value of the panel i, simulation n, when mined and processed with the COG j; WO 2005/124104 PCT/AU2005/000761 - 11 R," is the maximum mining capacity in period t; Q,' is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; 5 dt is the time discount factor; xt is the fraction of the panel i is extracted with the COG j in period t; y,, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 10 otherwise; 6, is a binary variable controlling the selection of the COG applied in period t. The invention may therefore be said to reside in a 15 computer program for scheduling the extraction of a resource within a medium comprising: code for receiving drillhole grade data relating to the resource contained within the medium; code for creating a plurality of different block 20 models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models honouring the drillhole grade data; and code for scheduling extraction of blocks based on a cut-off grade policy to increase expected.net present 25 value having regard to all of the block models. In the preferred embodiment of the invention, the code for creation of the plurality of different block models is code for performing a conditional simulation technique 30 which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and first and second order statistics of the resource represented, respectively, by probability distribution function and variogram. 35 In one embodiment of the invention, the program further comprises code for reducing the number of blocks in each WO 2005/124104 PCT/AU2005/000761 - 12 of the block models by aggregating blocks prior to the step of scheduling extraction. In one embodiment of the invention, the cut-off grade is a 5 fixed cut-off grade. However, in the preferred embodiment of the invention, the cut-off grade is a variable cut-off grade. In the preferred embodiment of the invention, the code for 10 scheduling includes code to provide an average of resource grade and block valuations and the average of the resource grade and block valuations is used to produce the scheduled extraction of blocks and the cut-off grade value to increase expected net present value. 15 In one embodiment of the invention, the code for scheduling further comprises code for optimizing extraction by considering the average valuation produced from the conditional simulations. 20 However, in the preferred embodiment of the invention, the code for scheduling further comprises code for using all of the conditional simulations to optimize simultaneously to produce a single extraction schedule with a variable 25 cut-off grade policy, and wherein all the conditional simulations are valued using this schedule and the valuations are averaged to produce the expected net present value, which is substantially better than the one obtained by optimizing the average grade block model. 30 In the preferred embodiment of the invention, the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 35 In the preferred embodiment, the code for scheduling using a variable cut-off grade converts from a non-linear WO 2005/124104 PCT/AU2005/000761 - 13 representation to a linear representation, and further comprises a mixed integer programming formulation is used for simultaneous optimisation of the extraction schedule and the cut-off grade. 5 Preferably the code for scheduling is to 1 N P G T Maximise -EEE V7xy, d') (1) N n=1 ,=1 i=1 ,=1 subject to the following constraints: 10 -ZERx, R,*,for allt (2) n=1 j=1 j-1 N P G I IQ xy, Q,,forallt (3) n=1 i=1 j=1 Y-1: syi,,.for all iand t (4) G 1] xr sy, ,for alli (5) r=1 j=1 G T 15 yi,5Ixkj,,for all i, tandkcSi (6) j1 =1 G I 5j, =1, for all t (7) j=1 x,, j, ,for all i, j,and t (8) 20 wherein; T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; 25 R," is the total rock in the panel i in. simulations n. Q is the total ore in the panel i, simulation n, when mined with the COG j; " is the value of the panel i, simulation n, when mined and processed with the COG j; 30 R,0 is the maximum mining capacity in period t; Q,0 is the maximum processing rate in period t; S, is the set of panels that must be removed before WO 2005/124104 PCT/AU2005/000761 - 14 starting the panel i; d is the time discount factor; x, is the fraction of the panel i is extracted with the COG j in period t; 5 y, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise;

J

5 , is a binary variable controlling the selection of the COG applied in period t. 10 The invention may therefore be said to reside in a method of mining a mine comprising a resource within a medium, the method comprising: removing medium and resource from the mine using an 15 extraction schedule created by the following steps; obtaining drillhole grade data relating to the resource contained within the medium; creating a plurality of different block models, each formed of a plurality of blocks, each block in each model 20 having a resource grade, each of the block models honouring the drillhole grade data; and scheduling extraction of blocks based on a cut-off grade policy to increase expected net present value having regard to all of the block models. 25 In the preferred embodiment of the invention, the creation of the plurality of different block models is a conditional simulation technique which provides the generation of a plurality of equally probable block model 30 realisations, all of which honour the data and first and second order statistics of the resource represented, respectively, by probability distribution function and variogram. 35 In one embodiment of the invention, the number of blocks in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction.

WO 2005/124104 PCT/AU2005/000761 - 15 Instances of block aggregations are assumed to have homogeneous properties within each grade bin. In one embodiment of the invention, the cut-off grade is a 5 fixed cut-off grade. However, in the preferred embodiment of the invention, the cut-off grade is a variable cut-off grade. In the preferred embodiment of the invention, the 10 plurality of block models are used to provide an average of resource grade and block valuations and the average of the resource grade and block valuations is used to produce the scheduled extraction of blocks and the cut-off grade value to increase expected net present value. 15 In one embodiment of the invention, the scheduled extraction of blocks may be optimised by considering the average valuation produced from the conditional simulations. 20 However, in the preferred embodiment of the invention, all of the conditional simulations are optimized simultaneously to produce a single extraction schedule with a variable cut-off grade policy, and wherein all the 25 conditional simulations are valued using this schedule and the valuations are averaged to produce the expected net present value, which is substantially better than the one obtained by optimizing the average grade block model. 30 In the preferred embodiment of the invention, the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 35 In the preferred embodiment, the step of scheduling using a variable cut-off grade is converted from a non-linear representation to a linear representation, and a mixed WO 2005/124104 PCT/AU2005/000761 - 16 integer programming formulation is used for simultaneous optimisation of the extraction schedule and the cut-off grade. 5 Preferably the scheduling is determined by (1 N P G T Maximise -, EEEI V" x , d' (1) N n=1 i =I y,11 subject to the following constraints: 1N PG 10 1_ E R,"xt R,,for allt (2) N n=I i=1 j-1 N P G :EZQn x, M: Qt, for all t (3 ) N n=1 i=1 j=I yi _ s: y, for all i and t ( 4) G Z xij, 5 y, for all i ( 5) r=1 j=I G T yit Z I x, , for all i, t and k c Si (6) j=I t=1 G 15 (5j, =1, for all t (7) j=I Xy 5j, ,,foralli,j,andt (8) wherein; 20 T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; R," is the total rock in the panel i in simulations n. 25 Q is the total ore in the panel i, simulation n, when mined with the COG j; V" is the value of the panel i, simulation n, when mined and processed with the COG j; R,4 is the maximum mining capacity in period t; 30 Q,4 is the maximum processing rate in period t; S, 'is the set of panels that must be removed before starting the panel i; WO 2005/124104 PCT/AU2005/000761 - 17 dt is the time discount factor; x0 is the fraction of the panel i is extracted with the COG j in period t; y 1 is a binary variable equal to 1 if the extraction of 5 the panel i has started in periods 1 to t, and equal to 0 otherwise; og is a binary variable controlling the selection of the COG applied in period t. 10 The invention may therefore be said to reside in a method of determining the net present value of a mine comprising a resource within a medium, the method comprising the steps of: obtaining drillhole grade data relating to the 15 resource contained within the medium; creating a plurality of different block models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models honouring the drillhole grade data; and 20 scheduling extraction of blocks based on a cut-off grade policy to provide the net present value having regard to all of the block models. In the preferred embodiment of the invention, the creation 25 of the plurality of different block models is a conditional simulation technique which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and first and second order statistics of the resource represented, 30 respectively, by probability distribution function and variogram. In one embodiment of the invention, the number of blocks in each of the block models is reduced by aggregating 35 blocks prior to the step of scheduling extraction.

WO 2005/124104 PCT/AU2005/000761 - 18 In one embodiment of the invention, the cut-off grade is a fixed cut-off grade. However, in the preferred embodiment of the invention, the cut-off grade is a variable cut-off grade. 5 In the preferred embodiment of the invention, the plurality of block models are used to provide an average of resource grade and block valuations and the average of the resource grade and block valuations is used to produce 10 the scheduled extraction of blocks and the cut-off grade value to increase expected net present value. In one embodiment of the invention, the scheduled extraction of blocks may be optimised by considering the 15 average valuation produced from the conditional simulations. However, in the preferred embodiment of the invention, all of the conditional simulations are optimized 20 simultaneously to produce a single extraction schedule with a variable cut-off grade policy, and wherein all the conditional simulations are valued using this schedule and the valuations are averaged to produce the expected net present value, which is substantially better than the 25 one obtained by optimizing the average grade block model. In the preferred embodiment of the invention, the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent 30 to waste or stockpiled for later processing. In the preferred embodiment, the step of scheduling using a variable cut-off grade is converted from a non-linear representation to a linear representation, and a mixed 35 integer programming formulation is used for simultaneous optimisation of the extraction schedule and the cut-off grade.

WO 2005/124104 PCT/AU2005/000761 - 19 Preferably the scheduling is determined by 1 N P G T Maximise -EIZ ii X Uxt d' (1) Nn=1 i=1 j=1 t=1 5 subject to the following constraints: 1 N P 11 ZRx1 R, ,for allt (2) n=1 i=1 j-I 1 :1Qn x,, Qt, ,for all t (3 ) Nn=1 i=1 j=1 yi11! y, ,for all iand t (4) t G 10 1:1:x,: sy,, ,for all i (5) r=1 j=1 G T yit ! E xk,, , for all i, t and k c Si (6) j=I i=1 G Z 8, =1, for allt (7) j=1 x,, Sj, , for all i, j,and t (8) 15 wherein; T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; 20 G is the number of all possible cutoff grades; R," is the total rock in the panel i in simulations n. Q is the total ore in the panel i, simulation n, when mined with the COG j; V is the value of the panel i, simulation n, when mined 25 and processed with the COG j; Ro is the maximum mining capacity in period t; Q,' is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; 30 dt is the time discount factor; xy, is the fraction of the panel i is extracted with the COG j in period t; WO 2005/124104 PCT/AU2005/000761 - 20 y, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise; 6, is a binary variable controlling the selection of the 5 COG applied in period t. A further aspect of the invention concerns reducing the difficulty of processing information in order to provide extraction schedules and net present value estimations. 10 According to a further aspect, the invention may be said to reside in a method for scheduling the extraction of a resource contained within a medium, comprising the steps of: 15 obtaining drillhole grade data relating to the resource contained within the medium; creating a block model with each block in the model having a resource grade; scheduling extraction of blocks based on a variable 20 cut-off grade value to increase expected net present value, the variable cut-off grade being represented by a non-linear function; and converting the non-linear representation to a linear representation and using a mixed integer programming 25 formulation for a simultaneous optimisation of an extraction sequence provided by the scheduled extraction and cut-off grade policy. Thus, the method according to this aspect of the invention 30 simplifies processing because of the linearisation of the problem and also therefore enables a variable cut-off grade to be used in the optimisation instead of a fixed cut-off grade value. 35 In the preferred embodiment of the invention, the non linear function is converted to a plurality of linear values by producing a plurality of step values relating to WO 2005/124104 PCT/AU2005/000761 - 21 the non-linear function, and selecting one of the values to provide the cut-off grade decision by the mixed integer programming function. 5 In one embodiment of the invention, the block model is a deterministic model built by spatial interpretation of the data using a Kriging procedure. Thus, in this embodiment, only a single model is provided. 10 However, in the preferred embodiment of the invention, the block model is a conditional simulation producing a plurality of different block models providing a number of equally probably realisations of block models honouring the drillhole data and first and second order statistics 15 of the resource distributed in the medium represented, respectively, by probability distribution and variogram. In one embodiment of the invention, the number of blocks in each of the block models is reduced by aggregating 20 blocks prior to the step of scheduling extraction. In the preferred embodiment of the invention, the plurality of block models are used to provide an average of resource grade valuations and the average of the 25 resource grade valuations is used to produce the scheduled extraction of blocks and the cut-off grade value to increase expected net present value. In one embodiment of the invention, the scheduled 30 extraction of blocks may be optimised by considering the average valuation produced from the conditional simulations. However, in the preferred embodiment of the invention, all 35 of the conditional simulations are optimised to produce a plurality of extraction schedules, each having valuations and the valuations are averaged to produce the net present WO 2005/124104 PCT/AU2005/000761 - 22 value, and the extraction schedule comprises the extraction schedule which is closest to that giving the net present value. 5 In the preferred embodiment of the invention, the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 10 Preferably the scheduling is determined by SlN P G T Maximise -((((V" x, d' t Maximise N n=1 i=1 j=1 ,=1 JXtd subject to the following constraints: 15 NZEE Rx Rt, ,for allt (2) N n=1 i=1 j-1 N P G -: E 1 Q" xy, 5 Q,* , for all t (3) N n=1 i=1 j=1 y ,_, &! yi, , for all i and t ( 4) t G E E xi, y , for all i (5) r=1 j=1 G T yi, ! Z x , for all i, t and k cSi (6) j=1 t=1 .G 20 S=j,=1,forallt (7) j=1 xU1 , ,,for all i,j, and t (8) wherein; 25 T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; R," is the total rock in the panel i in simulations n. 30 Qg is the total ore in the panel i, simulation n, when mined with the COG j; V," is the value of the panel i, simulation n, when mined WO 2005/124104 PCT/AU2005/000761 - 23 and processed with the COG j; Ro is the maximum mining capacity in period t; Q,0 is the maximum processing rate in period t; S, is the set of panels that must be removed before 5 starting the panel i; dt is the time discount factor; x., is the fraction of the panel i is extracted with the COG j in period t; y,, is a binary variable equal to 1 if the extraction of 10 the panel i has started in periods 1 to t, and equal to 0 otherwise; og is a binary variable controlling the selection of the COG applied in period t. 15 According to a further aspect, the invention may be said to reside in an apparatus for scheduling the extraction of a resource contained within a medium, comprising: a processor for receiving drillhold grade data relating to the resource contained within the medium; 20 a processor for creating a block model with each block in the model having a resource grade; a processor for scheduling extraction of blocks based on a variable cut-off grade value to increase expected net present value, the variable cut-off grade being 25 represented by a non-linear function; and a processor for converting the non-linear representation to a linear representation and using a mixed integer programming formulation for a simultaneous optimisation of an extraction sequence provided by the 30 scheduled extraction and cut-off grade policy. In the preferred embodiment of the invention, the processor is for converting the non-linear function to a plurality of linear values by producing a plurality of 35 step values relating to the non-linear function, and selecting one of the values to provide the cut-off grade decision by the mixed integer programming function.

WO 2005/124104 PCT/AU2005/000761 - 24 In one embodiment of the invention, the block model is a deterministic model built by spatial interpretation of the data using a Kriging procedure. Thus, in this embodiment, 5 only a single model is provided. However, in the preferred embodiment of the invention, the block model is a conditional simulation producing a plurality of different block models providing a number of 10 equally probably realisations of block models honouring the drillhole data and first and second order statistics of the resource distributed in the medium represented, respectively, by probability distribution and variogram. 15 In one embodiment of the invention, the number of blocks in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction. In the preferred embodiment of the invention, the 20 processor is for averaging of resource grade valuations from the block models and the average of the resource grade valuations is used to produce the scheduled extraction of blocks and the cut-off grade value to increase expected net present value. 25 In one embodiment of the invention, the scheduled extraction of blocks may be optimised by considering the average valuation produced from the conditional simulations. 30 However, in the preferred embodiment of the invention, the processor is for optimizing all of the conditional simulations to produce a plurality of extraction schedules, each having valuations and the valuations are 35 averaged to produce the net present value, and the extraction schedule comprises the extraction schedule which is closest to that giving the net present value.

WO 2005/124104 PCT/AU2005/000761 - 25 In the preferred embodiment of the invention, the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent 5 to waste or stockpiled for later processing. Preferably the scheduling is determined by the processor by the following: i N P G T 10 Maximise - l Vxl, d' (1) N n= =i=t= subject to the following constraints: 1N P G !z R,"x, R 0 .for allt (2) n=1 i=1 j-1 1N P G 15 - E Z Q x,,:5Q, ,for all t (3) N n=1 i=1 j=1 yi t1: y,, , for all i and t ( 4) G xE XY: y,, , for all i ( 5) r=1 j=1 G T Y,, Z Z xkj, , for all i, t and k c Si (6) j=1 t=1 G Z 3, =1,for allt (7) j=1 20 xj 5j ,,foralli,j,andt (8) wherein; T is the number of scheduling periods; N is the number of simulations; 25 P is the total number of panels; G is the number of all possible cutoff grades; R," is the total rock in the panel i in simulations n. Q is the total ore in the panel i, simulation n, when mined with the COG j; 30 V7 is the value of the panel i, simulation n, when mined and processed with the COG j; R,4 is the maximum mining capacity in period t; WO 2005/124104 PCT/AU2005/000761 - 26 Q, is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; dt is the time discount factor; 5 xj, is the fraction of the panel i is extracted with the COG j in period t; y, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise; 10 8, is a binary variable controlling the selection of the COG applied in period t. According to a further aspect, the invention may be said to reside in a computer program for scheduling the 15 extraction of a resource contained within a medium, comprising the steps of: code for receiving drillhole grade data relating to the resource contained within the medium; code for creating a block model with each block in 20 the model having a resource grade; code for scheduling extraction of blocks based on a variable cut-off grade value to increase expected net present value, the variable cut-off grade being represented by a non-linear function; and 25 code for converting the non-linear representation to a linear representation and using a mixed integer programming formulation for a simultaneous optimisation of an extraction sequence provided by the scheduled extraction and cut-off grade policy. 30 In the preferred embodiment of the invention, the non linear function is converted to a plurality of linear values by producing a plurality of step values relating to the non-linear function, and selecting one of the values 35 to provide the cut-off grade decision by the mixed integer programming function.

WO 2005/124104 PCT/AU2005/000761 - 27 In one embodiment of the invention, the block model is a deterministic model built by spatial interpretation of the data using a Kriging procedure. Thus, in this embodiment, only a single model is provided. 5 However, in the preferred embodiment of the invention, the block model is a conditional simulation producing a plurality of different block models providing a number of equally probably realisations of block models honouring 10 the drillhole data and first and second order statistics of the resource distributed in the medium represented, respectively, by probability distribution and variogram. In one embodiment of the invention, the number of blocks 15 in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction. In the preferred embodiment of the invention, the plurality of block models are used to provide an average 20 of resource grade valuations and the average of the resource grade valuations is used to produce the scheduled extraction of blocks and the cut-off grade value to increase expected net present value. 25 In one embodiment of the invention, the scheduled extraction of blocks may be optimised by considering the average valuation produced from the conditional simulations. 30 However, in the preferred embodiment of the invention, all of the conditional simulations are optimised to produce a plurality of extraction schedules, each having valuations and the valuations are averaged to produce the net present value, and the extraction schedule comprises the 35 extraction schedule which is closest to that giving the net present value.

WO 2005/124104 PCT/AU2005/000761 - 28 In the preferred embodiment of the invention, the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 5 Preferably the scheduling is determined by ( 1 N P G T Maximise -±Z ( V"xY, 1 d' (1) N ,= = =1= subject to the following constraints: 10 -JEJ: R,"x, R ,forallt (2) Nn=1 i=1 j-1 N P G 1 1Qx ,Q ,for allt (3) N n=1 i=1 j=1 y,,.5 yi,,foralliand t (4) G xy,:5y,, ,for all i (5) r=1 j=1 G T 15 Yi, EXkir , for all i, t and k cSi (6) j=1 1=1 G Z 5 j, =1, for all t (7) j=1 xj, <oj, ,for all i,j, and t (8) 20 wherein; T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; 25 R," is the total rock in the panel i in simulations n. Qn is the total ore in the panel i, simulation n, when mined with the COG j; Vj is the value of the panel i, simulation n, when mined and processed with the COG j; 30 R,0 is the maximum mining capacity in period t; Q, is the maximum processing rate in period t; S, is the set of panels that must be removed before WO 2005/124104 PCT/AU2005/000761 -29 starting the panel i; dt is the time discount factor; x, is the fraction of the panel i is extracted with the COG j in period t; 5 yi, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise; 9, is a binary variable controlling the selection of the COG applied in period t. 10 Brief Description of the Drawings Embodiments of the invention will be described, by way of example, with reference to the accompanying drawings in which: 15 Figure 1 is a diagram showing an open cut mine in which a resource is distributed in a medium (such as an ore body in earth and rock) in which a block model has been overlaid; Figure 2 is a graph showing COG against time; 20 Figure 3 is a block diagram according to one embodiment of the invention; Figure 4 is a graph showing the non-linear representation of variable cut-off grade; Figure 5 is a block diagram showing net present value 25 estimations according to embodiments of the invention; Figures 6A and 6B are diagrams showing various valuations of net present value; Figure 7 is a graph showing mining schedule optimised with a marginal cut-off grade; 30 Figure 8 is a graph showing net present value of the schedule optimisation with the marginal cut-off grade; Figure 9 is a graph showing a mining schedule optimised with a mean grade model; Figure 10 is a graph showing net present value of the 35 scheduled optimisation with the mean grade model; Figure 11 is a graph showing mining schedule optimised with a set of conditional simulations; WO 2005/124104 PCT/AU2005/000761 - 30 Figure 12 is a graph showing net present value of the schedule optimised with the set of conditional simulations. 5 Detailed Description of the Embodiments With reference to Figure 1 an open cut mine 10 is shown which is comprised of a resource such as an ore body which is distributed in a medium such as earth and rock. 10 The open cut mine 10 is to be mined over a significant period, such as 10-30 years, and the preferred embodiment of the invention is concerned with optimising return from the mine by scheduling mining phases to maximum net present value of the mine 10. 15 The mine 10 is initially surveyed by drilling a plurality of drillholes represented by the lines 11 to provide data relating to the resource contained within the mine. The data comprises the grade of the resource at particular 20 locations of the drill holes 11. Typically the data will provide a value which represents the grade of the resource as a percentage such as 0.1%, 0.6%, 1%, etc. A block model of the mine 10 is created which is formed 25 from a plurality of blocks 12. Typically, the block model may contain up to half a million blocks. The blocks 12 provide an ore body model.of blocks built up by the spatial interpolation of the data obtained from the 30 drillholes 11. Typically, it is desirable to mine the most valuable parts of the mine as quickly or as early as possible. Also typically, the most valuable parts of the mine are likely 35 to be well below ground level, which means that a large amount of rock and earth needs to be removed before the most valuable deposits are uncovered. The preferred WO 2005/124104 PCT/AU2005/000761 - 31 embodiments of the invention therefore seek to optimise the mining schedule so that regions of the mine represented by various ones of the blocks are removed to enable the most valuable deposits to be uncovered as 5 quickly as possible. The blocks which are mined prior to those blocks need to be considered to determine whether they are to be forwarded to waste, stockpiled for later processing, or forwarded direct to processing. By interpolating the grade of the ore in each of the blocks, 10 a decision can be made at a particular time as to whether it is desirable to forward the mined material for processing or whether that is uneconomical and the material should merely be sent to waste or stockpiled for later processing. 15 In order to determine whether a block is forwarded for processing, stockpiled or sent to waste, a cut-off grade value is established. For example, the cut-off grade at a particular time may be 0.6%. Thus, any blocks which have 20 a grade above 0.6% are sent for processing. However, if the grade is below 0.6%, a decision is made as to whether the block is simply sent to waste or whether it is stockpiled for later processing. Figure 2 is a graph which illustrates cut-off grade with respect to time and 25 shows three blocks which may need to be mined at particular times. For example, block 12a which is mined at time TO may have a grade which does not justify processing at time TO but which will justify processing at a later time Ti. Thus, the block 12a can be stockpiled 30 for processing at the later time T1. Block 12b which is mined at T2 has a grade which justifies processing at the time T2, and therefore may be sent immediately to processing. Block 12c has a grade which is below the cut off grade and at no future time has a grade which equals 35 the cut-off grade and therefore, is sent to waste. The preferred embodiment of the invention targets all phases of mine optimisation, including the net present value WO 2005/124104 PCT/AU2005/000761 - 32 optimal block extraction sequence, pushback design, and simultaneous cut-off grade and mining schedule optimisation. 5 One embodiment of the invention is shown in Figure 3. With reference to this figure, the data acquired from the bore holes 11 is used to create a deterministic block model 201 to thereby provide the blocks 12 having assigned grade values. The deterministic block model 201 is 10 created by spatial interpolation of drillhole data, preferably using a Kriging procedure. The number of blocks in the block model 301 may be reduced by aggregation of the blocks at step 302. The aggregation 15 of the blocks is disclosed in our International Patent Application Nos. PCT/AU2003/001298 and PCT/AU2003/001299. The contents of these two International applications are incorporated into this specification by this reference. 20 However the invention is applicable to any sort of aggregation that is spatially connected and respects geometrical block extraction precedence constraints. At step 303, a block removal schedule is created based on 25 a variable cut-off grade to maximise net present value of the mine 10. The use of a variable cut-off grade introduces a significant disadvantage in that optimisation of the 30 schedule to obtain maximised net present value using variable cut-off grade in its direct formulation leads to a non-linear problem. This is much harder to solve than a linear problem or linear value which occurs if a fixed variable cut-off grade is used in the formulation. To 35 overcome this disadvantage, the preferred embodiment of the invention provides a linearisation of the problem, making it possible to use a mixed integer programming WO 2005/124104 PCT/AU2005/000761 - 33 formulation for a simultaneous optimisation of the extraction sequence and cut-off grade. Figure 4.is a graph representing the non-linear cut-off 5 grade 15 which is a plot of total tonnes in a block with a grade greater than grade plotted against grade. In order to linearise the problem, the function 15 is divided into a number of linear steps represented at 17, each providing a discrete grade value such as G1, as shown in Figure 4. 10 The grade value G1 is selected using a mixed integer programming formulation for the simultaneous optimisation of the extraction sequence and the cut-off grade. Thus, this enables a schedule of block removal to be 15 determined using a variable cut-off grade which maximises the net present value of the mine. Thus, return from the mining operation is optimised at any particular time. Figure 5 is a block diagram showing a still further 20 embodiment of the invention, and the most preferred embodiment of the invention. In Figure 5 the data previously described is used to create multiple block models because the drillhole data is typically too sparse to support a unique and deterministic block model in 25 accordance with the embodiment of Figure 3. Thus, the likelihood of the block model produced with reference to Figure 3 matching actual reality is uncertain and questionable. Thus, a better approach is to use multiple block models which are created by conditional simulation 30 techniques. The conditional simulation techniques allow the generation of a number of equally probably block models, all of which honour the drillhole data obtained from the drillholes 11 and first and second order statistics of the ore body represented, respectively, by 35 probably distribution and variogram (e.g. Isaaks and Srivastava, 1989).

WO 2005/124104 PCT/AU2005/000761 - 34 Once again, the number of blocks in each of the models 501 may be reduced by aggregation 502 in the manner referred to in the aforementioned International applications. The block aggregation formed according to the algorithm set 5 forth in the aforesaid International application, preserves the required slope constraints and is very flexible allowing the user to fully control the size and shape of the aggregation. The optimisation according to the embodiments of the invention can be applied to any 10 aggregation of blocks with a set of precedence arcs, prescribing which blocks should be extracted before a particular given block. At step 303, a schedule 503 of block removal is created 15 dependent upon a cut-off grade value to maximise the net present value of the mine. The cut-off grade value which is used may be a fixed cut-off grade value 503a shown in Figure 5 or, in the most preferred embodiment of the invention, a variable cut-off grade value 503b which is 20 the same as the variable cut-off grade value described with reference to Figure 3. The simplest and most straight forward use of the multiple block models 501 is to estimate the variability in the 25 projected net present value associated with the ore body uncertainty by valuing the optimal schedule obtained from the Kriged model through each of the conditionally simulated realisations. This is shown in Figure 6A in which, for example, ten conditional simulations are 30 generated and those ten conditional simulations are immediately averaged at step 601. Figures 3 and 5 schematically show an apparatus in the form of a processor which may be a personal computer or 35 the like for performing the method of the preferred embodiment.

WO 2005/124104 PCT/AU2005/000761 - 35 Figures 6A and 6B are an actual example of valuations of a known pit. Figure 6A shows the preferred embodiment of the invention 5 and is a schedule determination which gives the best estimate of the net present value of a mine. In Figure 6A it is noted that the optimization uses all ten simulations of the conditional simulation represented in Figure 6A. 10 Figure 6B shows optimization using only an average grade rather than all ten simulations, and shows that an increased net present value is still obtained if the average of the ten valuations of the conditional simulations are used later rather than earlier in the 15 process. In Figure 6B the term E-type represents the average of the conditional simulations, equivalent to the Kriged model. Here, we explain the difference between the $671M and 20 $739M valuations that use identical extraction and cut-off grade (COG) schedules. In following this explanation, it is important to realise that each conditional simulation is an equally likely representation of the true block grade. 25 The extra value arises because the conditional simulations can be valued individually using the prescribed COG, then their value averaged, as compared to the case with the "average" E-type model where the averaging is done 30 earlier, at the grade level. A very simple single block model with no discounting will be used to illustrate this. Single Block Model Example The following data is used-the average value corresponds 35 to the E-type block model. This is a single block model with only three conditional simulations.

WO 2005/124104 PCT/AU2005/000761 - 36 Block mass (tonnes) 10,000 Number of blocks Conditional simulations 3 Grade in each conditional simulation (%Cu) [0.04, 0.06, 0.11] Average block grade (%Cu) [E-type value] 0.07 Tonnes of Cu in each conditional simulation [4, 6, 11] Average block tonnes of Cu [E-type value] 7 Revenue per tonne of Cu ($) 2,000 Block value as waste ($) (Zero mining cost!) 0 In valuing this block, we will use a marginal COG (Processing Cost = Revenue), and look at four situations, 5 represented by four different processing costs. The marginal COG is calculated thus: (Block Tonnes) * COG * (Revenue/tonne) = (Block Tonnes) * (Processing Cost/tonne) 10 COG = (Processing Cost/tonne) / (Revenue/tonne) 1. COG is less than all conditional- simulation grades. COG = 0.03% Cu, equivalent to a processing cost of 0.60 $/tonne. 15 2. COG is greater than only one of the conditional simulation grades. COG = 0.05% Cu, equivalent to a processing cost of 1.00 $/tonne. 20 3. COG is greater than two of the conditional simulation grades. COG = 0.08% Cu, equivalent to a processing cost of 1.60 $/tonne. 4. COG is greater than all conditional simulation 25 grades. COG = 0.12%Cu, equivalent to a processing cost of 2.40 $/tonne. The block valuation will be performed for each situation.

WO 2005/124104 PCT/AU2005/000761 - 37 Situation COG Block Coditional Simulaion Value($) Conditional E4ypeValue Conditional (%Cu) Pocessing 1 2 3 Simulation ($) Simulation Cost($) Avwg Value ($) Value lmium ($) 1 0.03 6000 2000 6000 16000 8000 8000 0 2 0.05 10000 0 2000 12000 4667 4000 667 3 0.08 16000 0 0 6000 2000 0 2000 4 0.12 24000 0 0 0 0 0 0 This example demonstrates that whenever the conditional simulated grades for a block are found to be distributed 5 on both sides of the COG, the block valuation for the average of the conditional simulations will be greater than the valuation of the average grade (Kriged) model. The following observations may be helpful: 10 e In situation 2, the E-type valuation effectively attributes negative value to the first conditional simulation, whereas when the conditional simulations are valued separately this instance is attributed exactly zero value - its value as waste. 15 * In situation 3, the E-type valuation ignores the value of one of the conditional simulations. The valuation using the conditional simulations is more realistic because when a block is mined, blast hole assays 20 can be used to determine a close estimate of the actual block grade. With this near-perfect information about block grade, we can make an accurate assessment of whether to process or waste a block. This is equivalent to valuing each conditional simulation individually as if 25 each were an actual reality. Thus, by using the multiple simulations to produce multiple valuations and then averaging the multiple valuations in order to produce the net present value, a 30 higher expected net present value is obtained than if the multiple simulations are first averaged and simply that WO 2005/124104 PCT/AU2005/000761 - 38 average is then used to produce a valuation. As previously mentioned, the preferred embodiment of the invention uses mixed integer programming formulation to 5 include the multiple conditional simulations and variable cut-off grades. This approach allows one to estimate the gap between the obtained solution and the upper theoretical limit. To illustrate this, an example of how optimisation is generated will be given by reference to a 10 simple case in which one rock-type contains one metal type, which can be processed through one processing plant, generalisation to the case of multiple rock types, metals and processing streams is cumbersome but straight forward. For simplicity, we consider below only the case of a 15 discrete set of cut-off grades that is possible to generalise the results to the continuous cut-off grade case. The following notations are used: T is the number of scheduling periods; 20 N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; R7 is the total rock in the panel i in simulations n. Q" is the total ore in the panel i, simulation n, when 25 mined with the COG j; V," is the value of the panel i, simulation n, when mined and processed with the COG j; R,' is the maximum mining capacity in period t; Q, is the maximum processing rate in period t; 30 S, is the set of panels that must be removed before starting the panel i; dt is the time discount factor; xYr is the fraction of the panel i is extracted with the COG j in period t; 35 y, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 WO 2005/124104 PCT/AU2005/000761 - 39 otherwise; S, is a binary variable controlling the selection of the COG applied in period t; 5 The MIP formulation is: Maximize -E 0X(1" xy, d' (1) N n=1 , =1= subject to the following constraints: 10 - Rx RO, for allt (2) N n=1 i=1 j-1 SN P G E Z Q x, Q, for all t (3) N n=1 i=1 j=1 yi ,_1 & y,, , for all i and t ( 4) G 1: xu, 5 yi, , for all i (5) r=1 j=1 G T 15 y,, E E xkjr , for all i, t and k cSi (6) j=1 t=1 G ZS*j, =1, for allt (7) j=1 Xyt j, ,for alli,j,andt (8) The objective function (1) represents the discounted cash 20 flow. Constraints (2) and (3) enforce the mining and processing limits on average. Constraints (4) - (6) enforce the panel extraction precedence constraints, and constraints (7) and (8) ensure that the same COG is applied to all panels extracted in any given time period. 25 This MIP formulation is solved by the commercially available software package CPLEX version 9.0, by ILOG Inc. To test the algorithm we have chosen 10 conditional simulations of a block model containing one type of metal 30 and using one processing plant. Because of confidentiality requirements all the economic parameters were rescaled and do not represent reality. However, all WO 2005/124104 PCT/AU2005/000761 - 40 the relative characteristics which demonstrate the potential of the new method are not affected by the rescaling. The ultimate pit for the design was chosen by using the Lersch-Grossmann algorithm (Lersch and 5 Grossmann, 1965) and the procedure similar to that used in Whittle Four-X software. The ultimate pit contains 191 million tonnes of rock and 62.9 ± 2.7 million tonnes of ore (above the marginal COG = 0.6 %). The undiscounted value in the ultimate pit (if processed with the marginal 10 COG) is $ (1,316 ± 99) million. It was divided into 6 mining phases and scheduled over 12 years. The mining rate was set to 30 MT/year and the processing rate to 5 MT/year. The initial capital investment was assumed to be $300 million, and the discount rate 10%. The base case 15 optimisation was done using the marginal COG and produced the discounted cash flow $(704 ± 31) million, and the NPV was $(404 ± 31) million. The mining schedule and the NPV are shown, respectively, in Fig. 7 and 8. The second optimisation was done using the variable COG, but was 20 based on the mean grade block model, i.e. it was similar to the one which can be generated by using one deterministic model. The schedule was then evaluated against all 10 realisations of ore body model and produced the NPV = $(485 ± 40) million, an increase of 20 % over 25 the base case. The results are shown in Fig. 9 and 10. The third optimisation was done using the algorithm described above, and produced the NPV = $(505 ± 43) million, a further increase of 4.1% over the case of mean grade based optimisation. The results are shown in Fig. 11 and 12. 30 The relative variability of NPV in all cases was roughly the same, about 8%. Another important result of the variable COG policy is that the pay-back period (defined here as the time when the cummulative NPV becomes equal to zero) is decreased from 5 to 3 years. 35 The increase of 4.1% in NPV may be not seen as a very substantial, but it should be mentioned that the block WO 2005/124104 PCT/AU2005/000761 - 41 model considered does not have a high variability. The relative variance in the undiscounted value of the ultimate pit is only 7.6%. There are many deposits which can have variability of the order of 20 - 30%. For these 5 kind of deposits the potential improvement in the expected NPV may be substantially high. Thus, according to the preferred embodiment of the invention, an extraction sequence and cut-off grade policy 10 can be determined which, when evaluated through the whole set of conditionally simulated ore bodies, will produce the best possible expected net present value. The degree of accuracy of this optimised schedule can be estimated precisely, in contrast to prior art techniques. 15 In the claims which follow and in the preceding description of the invention, except where the context requires otherwise due to express language or necessary implication, the word "comprise", or variations such as 20 "comprises" or "comprising", is used in an inclusive sense, ie. to specify the presence of the stated features but not to preclude the presence or addition of further features in various embodiments of the invention. 25 It is to be understood that the prior art publications referred to herein, do not constitute an admission that that the publication forms a part of the common general knowledge in the art, in Australia or in any other country. 30 Since modifications within the spirit and scope of the invention may readily be effected by persons skilled within the art, it is to be understood that this invention is not limited to the particular embodiment described by 35 way of example hereinabove.

Claims

1. A method of scheduling the extraction of a resource within a medium comprising the steps of: 5 obtaining drillhole grade data relating to the resource contained within the medium; creating a plurality of different block models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models 10 honouring the drillhole grade data; and scheduling extraction of blocks based on a cut-off grade policy to increase expected net present value having regard to all of the block models. 15 2. The method of claim 1 wherein the creation of the plurality of different block models is a conditional simulation technique which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and first and second order 20 statistics of the resource represented, respectively, by probability distribution function and variogram.

3. The method of claim 1 wherein the number of blocks in each of the block models is reduced by aggregating blocks 25 prior to the step of scheduling extraction.

4. The method of claim 1 wherein the cut-off grade is a variable cut-off grade. 30 5. The method of claim 1 wherein the plurality of block models are used to provide an average of resource grade and block valuations and the average of the resource grade and block valuations is used to produce the scheduled extraction of blocks and the cut-off grade value to 35 increase expected net present value. WO 2005/124104 PCT/AU2005/000761 - 43 6. The method of claim 1 wherein the scheduled extraction of blocks is optimised by considering the average valuation produced from the conditional simulations. 5

7. The method of claim 1 wherein all of the conditional simulations are optimized simultaneously to produce a single extraction schedule with a variable cut-off grade policy, and wherein all the conditional simulations are 10 valued using this schedule and the valuations are averaged to produce the expected net present value, which is substantially better than the one obtained by optimizing the average grade block model. 15 8. The method of claim 1 wherein the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 20 9. The method of claim 1 wherein the step of scheduling using a variable cut-off grade is converted from a non linear representation to a linear representation, and a mixed integer programming formulation is used for simultaneous optimisation of the extraction schedule and 25 the cut-off grade.

10. The method of claim 1 wherein the scheduling is I N P G T determined by Maximise - E V. xn, d' (1) (N n=1i1jIt1 Y Y 30 subject to the following constraints: 1 N P G0 -QEE R,"x,! R, ,for allt (2) Nn=1 i=1 j-1 i N P G -EZJ Q,)x., < Q," for all t ( 3) n=1 i=I j= Y yo-, 5 y,, , for all i and t (4) WO 2005/124104 PCT/AU2005/000761 - 44 G I xy, &g yi,, for all i ( 5) r=1 j=1 G T y,:, 1 Z xkj , for all i, t and k cSi (6) j=1 1=1 G Z 6 t =1, for all t (7) j=1 xYt 1 ,, for all i, j,and t (8) 5 wherein; T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; 10 G is the number of all possible cutoff grades; R," is the total rock in the panel i in simulations n. Qn is the total ore in the panel i, simulation n, when mined with the COG j; V, is the value of the panel i, simulation n, when mined 15 and processed with the COG j; R, is the maximum mining capacity in period t; Q, is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; 20 dt is the time discount factor; xy, is the fraction of the panel i is extracted with the COG j in period t; y,, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 25 otherwise; Sj, is a binary variable controlling the selection of the COG applied in period t.

11. An apparatus of scheduling the extraction of a 30 resource within a medium from which drillhole grade data relating to the resource contained within the medium has been obtained comprising: WO 2005/124104 PCT/AU2005/000761 - 45 a processor for creating a plurality of different block models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models honouring the drillhole grade data; and 5 a processor for scheduling extraction of blocks based on a cut-off grade policy to increase expected net present value having regard to all of the block models.

12. The apparatus of claim 11 wherein the processor is 10 for creating conditional simulations which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and first and second order statistics of the resource represented, respectively, by probability distribution function and 15 variogram.

13. The apparatus of claim 11 wherein the processor is also for reducing the number of blocks in each of the block models by aggregating blocks prior to the step of 20 scheduling extraction.

14. The apparatus of claim 11 wherein the cut-off grade is a variable cut-off grade. 25 15. The apparatus of claim 11 wherein the processor is for providing an average of resource grade and block valuations from the plurality of models, and the average of the resource grade and block valuations is used to produce the scheduled extraction of blocks and the cut-off 30 grade value to increase expected net present value.

16. The apparatus of claim 11 wherein the processor is for scheduling extraction of blocks optimised by considering the average valuation produced from the 35 conditional simulations. WO 2005/124104 PCT/AU2005/000761 - 46 17. The apparatus of claim 11 wherein all of the conditional simulations are optimized simultaneously to produce a single extraction schedule with a variable cut off grade policy, and wherein all the conditional 5 simulations are valued using this schedule and the valuations are averaged to produce the expected net present value, which is substantially better than the one obtained by optimizing the average grade block model. 10 18. The apparatus of claim 11 wherein the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 15 19. The apparatus of claim 11 wherein the processor is for converting from a non-linear representation to a linear representation when a variable cut-off grade is used, and a mixed integer programming formulation is used for simultaneous optimisation of the extraction schedule 20 and the cut-off grade.

20. The apparatus of claim 11 wherein the scheduling is determined by the processor by i N P G T Maximise -1nx V"x, d' (1) (N n=1 i=1 j=I ,=I 25 subject to the following constraints: iN P G nX. 0 (EE Rxy R,* ,for allt (2) n=1 i j-1 iN P G 1-ZE Q x., <Q,',forallt (3) n=1 i=1 j=1 30 y _,s:5yi,, for all iand t (4) i G I xi, yi, ,for alli (5) r=I j=I G T Yi u X xk, , for all i, t and k cSi (6) j1l 1=1 WO 2005/124104 PCT/AU2005/000761 - 47 G Y,,= 1, for all t (7) j=1 xj 5 j ,foralli,j,andt (8) wherein; 5 T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; R2' is the total rock in the panel i in simulations n. 10 Qg is the total ore in the panel i, simulation n, when mined with the COG j; V,' is the value of the panel i, simulation n, when mined and processed with the COG j; Ro is the maximum mining capacity in period t; 15 Q,' is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; dt is the time discount factor; x., is the fraction of the panel i is extracted with the 20 COG j in period t; y,, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise; 3, is a binary variable controlling the selection of the 25 COG applied in period t.

21. A computer program for scheduling the extraction of a resource within a medium comprising: code for receiving drillhole grade data relating to 30 the resource contained within the medium; code for creating a plurality of different block models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models honouring the drillhole grade data; and WO 2005/124104 PCT/AU2005/000761 - 48 code for scheduling extraction of blocks based on a cut-off grade policy to increase expected net present value having regard to all of the block models. 5 22. The program of claim 21 wherein the code for creation of the plurality of different block models is code for performing a conditional simulation technique which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and 10 first and second order statistics of the resource represented, respectively, by probability distribution function and variogram.

23. The program of claim 21 wherein the program further 15 comprises code for reducing the number of blocks in each of the block models by aggregating blocks prior to the step of scheduling extraction.

24. The program of claim 21 wherein the code for 20 scheduling further comprises code for optimizing extraction by considering the average valuation produced from the conditional simulations.

25. The program of claim 21 wherein the code for. 25 scheduling further comprises code for using all of the conditional simulations to optimize simultaneously to produce a single extraction schedule with a variable cut off grade policy, and wherein all the conditional simulations are valued using this schedule and the 30 valuations are averaged to produce the expected net present value, which is substantially better than the one obtained by optimizing the average grade block model.

26. The program of claim 21 wherein the extraction 35 schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. WO 2005/124104 PCT/AU2005/000761 - 49 27. The program of claim 21 wherein the code for scheduling using a variable cut-off grade converts from a non-linear representation to a linear representation, and 5 further comprises a mixed integer programming formulation is used for simultaneous optimisation of the extraction schedule and the cut-off grade.

28. The program of claim 21 wherein the code for 10 scheduling is to Maximise -VE V x, dt (1) NI n=1 ,=1 j=1 ,=i subject to the following constraints: 1 N PG -E ERi x: sR,,forallt (2) N n=1 i=1 j-I i N P G 1 5 - Z Q xJ, & Q,, for all t (3) N n=1 i=1 j=1 yi ,_1:5 y,,, for all i and t ( 4) G 1:( x,,5 yi, , for all i ( 5) ,=1 j=I G T Yi a xkj, , for all i, t and k cSi (6) j=1 1=1 G E c5 =1, for all t (7) j=1 20 x,55,,foralli,j,andt (8) wherein; T is the number of scheduling periods; 25 N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; R" is the total rock in the panel i in simulations n. Q is the total ore in the panel i, simulation n, when 30 mined with the COG j; V,, is the value of the panel i, simulation n, when mined and processed with the COG j; WO 2005/124104 PCT/AU2005/000761 - 50 R,0 is the maximum mining capacity in period t; Q, is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; 5 dt is the time discount factor; xy, is the fraction of the panel i is extracted with the COG j in period t; y,, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 10 otherwise; 8, is a binary variable controlling the selection of the COG applied in period t.

29. A method of mining a mine comprising a resource 15 within a medium, the method comprising: removing medium and resource from the mine using an extraction schedule created by the following steps; obtaining drillhole grade data relating to the resource contained within the medium; 20 creating a plurality of different block models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models honouring the drillhole grade data; and scheduling extraction of blocks based on a cut-off 25 grade policy to increase expected net present value having regard to all of the block models.

30. The method of claim 29 wherein the creation of the plurality of different block models is a conditional 30 simulation technique which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and first and second order statistics of the resource represented, respectively, by probability distribution function and variogram. 35 WO 2005/124104 PCT/AU2005/000761 - 51 31. The method of claim 29 wherein the number of blocks in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction. 5 32. The method of claim 29 wherein the plurality of block models are used to provide an average of resource grade and block valuations and the average of the resource grade and block valuations is used to produce the scheduled extraction of blocks and the cut-off grade value to 10 increase expected net present value.

33. The method of claim 29 wherein the scheduled extraction of blocks is optimised by considering the average valuation produced from the conditional 15 simulations.

34. The method of claim 29 wherein all of the conditional simulations are optimized simultaneously to produce a single extraction schedule with a variable cut-off grade 20 policy, and wherein all the conditional simulations are valued using this schedule and the valuations are averaged to produce the expected net present value, which is substantially better than the one obtained by optimizing the average grade block model. 25

35. The method of claim 29 wherein the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 30

36. The method of claim 29 wherein the step of scheduling using a variable cut-off grade is converted from a non linear representation to a linear representation, and a mixed integer programming formulation is used for 35 simultaneous optimisation of the extraction schedule and the cut-off grade. WO 2005/124104 PCT/AU2005/000761 - 52 37. The method of claim 29 wherein the scheduling is i N P G T' determined by Maximise -J E V'"xud' (1) (N n=1 i=1 j=1 ,I subject to the following constraints: 5 1N PG RxyR, ,foralit (2) n=1 j=1 j-, iN P G - E E Q" x, <Q,* , for all t (3) Nn=1 i=1 j=I y _ 1 ! yi, ,for all i and t (4) G 1: xy, ! y, , for all i (5) r= j=I G T 10 Ya:,ZX x,,for all i, tandk cS (6) j=I t=1 G (j, =1, for allt (7) j=1 xij, 5, , for all i,j, and t (8) 15 wherein; T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; 20 R," is the total rock in the panel i in simulations n. Qu is the total ore in the panel i, simulation n, when mined with the COG j; V is the value of the panel i, simulation n, when mined and processed with the COG j; 25 R,0 is the maximum mining capacity in period t; Q,' is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; dt is the time discount factor; 30 x,, is the fraction of the panel i is extracted with the COG j in period t; y,, is a binary variable equal to 1 if the extraction of WO 2005/124104 PCT/AU2005/000761 - 53 the panel i has started in periods 1 to t, and equal to 0 otherwise; J, is a binary variable controlling the selection of the COG applied in period t. 5

38. A method of determining the net present value of a mine comprising a resource within a medium, the method comprising the steps of: obtaining drillhole grade data relating to the 10 resource contained within the medium; creating a plurality of different block models, each formed of a plurality of blocks, each block in each model having a resource grade, each of the block models honouring the drillhole grade data; and 15 scheduling extraction of blocks based on a cut-off grade policy to provide the net present value having regard to all of the block models.

39. The method of claim 38 wherein the creation of the 20 plurality of different block models is a conditional simulation technique which provides the generation of a plurality of equally probable block model realisations, all of which honour the data and first and second order statistics of the resource represented, respectively, by 25 probability distribution function and variogram.

40. The method of claim 38 wherein the number of blocks in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction. 30

41. The method of claim 38 wherein the cut-off grade is a variable cut-off grade.

42. The method of claim 38 wherein the plurality of block 35 models are used to provide an average of resource grade and block valuations and the average of the resource grade and block valuations is used to produce the scheduled WO 2005/124104 PCT/AU2005/000761 - 54 extraction of blocks and the cut-off grade value to increase expected net present value.

43. The method of claim 38 wherein the scheduled 5 extraction of blocks is optimised by considering the average valuation produced from the conditional simulations.

44. The method of claim 38 wherein all of the conditional 10 simulations are optimized simultaneously to produce a single extraction schedule with a variable cut-off grade policy, and wherein all the conditional simulations are valued using this schedule and the valuations are averaged to produce the expected net present value, which is 15 substantially better than the one obtained by optimizing the average grade block model.

45. The method of claim 38 wherein the extraction schedule determines whether resource and medium 20 corresponding to each block is sent for processing, sent to waste or stockpiled for later processing.

46. The method of claim 38 wherein the step of scheduling using a variable cut-off grade is converted from a non 25 linear representation to a linear representation, and a mixed integer programming formulation is used for simultaneous optimisation of the extraction schedule and the cut-off grade. 30 47. The method of claim 38 wherein the scheduling is i N P G T determined by Maximise -ZZ E V" x, d'J (1) (N n=1 i=1 j=I1 =1 Y U subject to the following constraints: 1 N P G 35 -JEE RxO,: R, ,for allt (2) n=1 i=l j-1 WO 2005/124104 PCT/AU2005/000761 - 55 :ZZQn x,,5 Q, , for all t (3 ) n=1 i=1 j=1 y, 1 yi, , for all i and t (4) t G ZEx r 5 yi, , for all i ( 5) r=I j=1 G T Yi 5 E Z xkj, , for all i, t and k cS (6) j=1 t=1 G 5 E5,=1,for allt (7) j=1 xi,,, t, for all i,j, and t (8) wherein; 10 T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; R, is the total rock in the panel i in simulations n. 15 Qn is the total ore in the panel i, simulation n, when mined with the COG j; V," is the value of the panel i, simulation n, when mined and processed with the COG j; R, is the maximum mining capacity in period t; 20 Q,' is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; dt is the time discount factor; x., is the fraction of the panel i is extracted with the 25 COG j in period t; y,, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise; 9, is a binary variable controlling the selection of the 30 COG applied in period t.

48. A method for scheduling the extraction of a resource contained within a medium, comprising the steps of: WO 2005/124104 PCT/AU2005/000761 - 56 obtaining drillhole grade data relating to the resource contained within the medium; creating a block model with each block in the model having a resource grade; 5 scheduling extraction of blocks based on a variable cut-off grade value to increase expected net present value, the variable cut-off grade being represented by a non-linear function; and converting the non-linear representation to a linear 10 representation and using a mixed integer programming formulation for a simultaneous optimisation of an extraction sequence provided by the scheduled extraction and cut-off grade policy. 15 49. The method of claim 48 wherein the non-linear function is converted to a plurality of linear values by producing a plurality of step values relating to the non linear function, and selecting one of the values to provide the cut-off grade decision by the mixed integer 20 programming function.

50. The method of claim 48 wherein the block model is a deterministic model built by spatial interpretation of the data using a Kriging procedure. 25

51. The method of claim 48 wherein the block model is a conditional simulation producing a plurality of different block models providing a number of equally probably realisations of block models honouring the drillhole data 30 and first and second order statistics of the resource distributed in the medium represented, respectively, by probability distribution and variogram.

52. The method of claim 48 wherein the number of blocks 35 in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction. WO 2005/124104 PCT/AU2005/000761 - 57 53. The method of claim 48 wherein the plurality of block models are used to provide an average of resource grade valuations and the average of the resource grade valuations is used to produce the scheduled extraction of 5 blocks and the cut-off grade value to increase expected net present value.

54. The method of claim 48 wherein the scheduled extraction of blocks is optimised by considering the 10 average valuation produced from the conditional simulations.

55. The method of claim 48 wherein all of the conditional simulations are optimised to produce a plurality of 15 extraction schedules, each having valuations and the valuations are averaged to produce the net present value, and the extraction schedule comprises the extraction schedule which is closest to that giving the net present value. 20

56. The method of claim 48 wherein the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent to waste or stockpiled for later processing. 25

57. The method of claim 48 wherein the scheduling is i N P G T determined by Maximise -- V nx, d' (1) (N n=1 ,i .= ,.1 Y U subject to the following constraints: 30 1 N P G -ZEE Rx R, , forallt (2) N n=I ,=1 j IN P G NZIE Q;x ,Q, forallt (3) N n=1 j=1 j=1 yi ,_, & y,, ,for all i and t ( 4) G E xy, ! y, , for all i (5) r=1 j=I WO 2005/124104 PCT/AU2005/000761 - 58 G T Yi, E < xkj, , for all i, t and k cS (6) j= t=1 G Si5jf =, 1, for allt (7) j=1 xy, 3S ,for all i, j,and t (8) 5 wherein; T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; 10 G is the number of all possible cutoff grades; R n is the total rock in the panel i in simulations n. Q is the total ore in the panel i, simulation n, when mined with the COG j; V is the value of the panel i, simulation n, when mined 15 and processed with the COG j; Ro is the maximum mining capacity in period t; Q,O is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; 20 dt is the time discount factor; x, is the fraction of the panel i is extracted with the COG j in period t; Y,, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 25 otherwise; 9, is a binary variable controlling the selection of the COG applied in period t.

58. An apparatus for scheduling the extraction of a 30 resource contained within a medium, comprising: a processor for receiving drillhold grade data relating to the resource contained within the medium; a processor for creating a block model with each block in the model having a resource grade; WO 2005/124104 PCT/AU2005/000761 - 59 a processor for scheduling extraction of blocks based on a variable cut-off grade value to increase expected net present value, the variable cut-off grade being represented by a non-linear function; and 5 a processor for converting the non-linear representation to a linear representation and using a mixed integer programming formulation for a simultaneous optimisation of an extraction sequence provided by the scheduled extraction and cut-off grade policy. 10

59. The apparatus of claim 58 wherein the processor is for converting the non-linear function to a plurality of linear values by producing a plurality of step values relating to the non-linear function, and selecting one of 15 the values to provide the cut-off grade decision by the mixed integer programming function.

60. The apparatus of claim 58 wherein the block model is a deterministic model built by spatial interpretation of 20 the data using a Kriging procedure.

61. The apparatus of claim 58 wherein the block model is a conditional simulation producing a plurality of different block models providing a number of equally 25 probably realisations of block models honouring the drillhole data and first and second order statistics of the resource distributed in the medium represented, respectively, by probability distribution and variogram. 30 62. The apparatus of claim 58 wherein the number of blocks in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction. 35 63. The apparatus of claim 58 wherein the processor is for averaging of resource grade valuations from the block models and the average of the resource grade valuations is WO 2005/124104 PCT/AU2005/000761 - 60 used to produce the scheduled extraction of blocks and the cut-off grade value to increase expected net present value. 5 64. The apparatus of claim 58 wherein the scheduled extraction of blocks is optimised by considering the average valuation produced from the conditional simulations. 10 65. The apparatus of claim 58 wherein the processor is for optimizing all of the conditional simulations to produce a plurality of extraction schedules, each having valuations and the valuations are averaged to produce the net present value, and the extraction schedule comprises 15 the extraction schedule which is closest to that giving the net present value.

66. The apparatus of claim 58 wherein extraction schedule determines whether resource and medium corresponding to 20 each block is sent for processing, sent to waste or stockpiled for later processing.

67. The apparatus of claim 58 wherein the scheduling is determined by the processor by the following: 25 N P G T Maximise - E V, xn, d' (1) (N n=1 i=1 j=1 ,=1 subject to the following constraints: 1 N P G 30 -ZEE R,"x,& R, ,for allt (2) Sn=1 i=1 j-1 1 N P G QE x, <Q, ,for allt (3) n=1 i=1 j= ylt-I yi , , foralliandt (4) t G ExY5 yi, ,for alli (5) r=1 j=1 WO 2005/124104 PCT/AU2005/000761 - 61 G T Y, ZZ Xk, , for all i, t and k c Si (6) j=1 t=1 G i c5, =1, for all t (7) j=1 xyt:8,, ,,for all i,j, and t (8) 5 wherein; T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; 10 R, is the total rock in the panel i in simulations n. Q, is the total ore in the panel i, simulation n, when mined with the COG j; Vo" is the value of the panel i, simulation n, when mined and processed with the COG j; 15 R, 0 is the maximum mining capacity in period t; Q,' is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; de is the time discount factor; 20 x, is the fraction of the panel i is extracted with the COG j in period t; yi, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise; 25 8, is a binary variable controlling the selection of the COG applied in period t.

68. A computer program for scheduling the extraction of a resource contained within a medium, comprising the steps 30 of: code for receiving drillhole grade data relating to the resource contained within the medium; code for creating a block model with each block in the model having a resource grade; WO 2005/124104 PCT/AU2005/000761 - 62 code for scheduling extraction of blocks based on a variable cut-off grade value to increase expected net present value, the variable cut-off grade being represented by a non-linear function; and 5 code for converting the non-linear representation to a linear representation and using a mixed integer programming formulation for a simultaneous optimisation of an extraction sequence provided by the scheduled extraction and cut-off grade policy. 10

69. The program of claim 68 wherein the non-linear function is converted to a plurality of linear values by producing a plurality of step values relating to the non linear function, and selecting one of the values to 15 provide the cut-off grade decision by the mixed integer programming function.

70. The program of claim 68 wherein the block model is a deterministic model built by spatial interpretation of the 20 data using a Kriging procedure.

71. The program of claim 68 wherein the block model is a conditional simulation producing a plurality of different block models providing a number of equally probably 25 realisations of block models honouring the drillhole data and first and second order statistics of the resource distributed in the medium represented, respectively, by probability distribution and variogram. 30 72. The program of claim 68 wherein the number of blocks in each of the block models is reduced by aggregating blocks prior to the step of scheduling extraction.

73. The program of claim 68 wherein the plurality of 35 block models are used to provide an average of resource grade valuations and the average of the resource grade valuations is used to produce the scheduled extraction of WO 2005/124104 PCT/AU2005/000761 - 63 blocks and the cut-off grade value to increase expected net present value.

74. The program of claim 68 wherein the scheduled 5 extraction of blocks is optimised by considering the average valuation produced from the conditional simulations.

75. The program of claim 68 wherein all of the 10 conditional simulations are optimised to produce a plurality of extraction schedules, each having valuations and the valuations are averaged to produce the net present value, and the extraction schedule comprises the extraction schedule which is closest to that giving the 15 net present value.

76. The program of claim 68 wherein the extraction schedule determines whether resource and medium corresponding to each block is sent for processing, sent 20 to waste or stockpiled for later processing.

77. The program of claim 68 wherein the scheduling is 1 N P G T determined by Maximise 112: Vonxy, d' (1) (Nn=I i=1 j=I1 =1 25 subject to the following constraints: N P G0 -1:]R,"x,, & R, for all t (2 ) Nn=1 i=1 j-I N P G N 2Q, xY :Q, ,for all t (3) Nn=l i=1 j=I yH-.I y, ,,foralliandt (4) , G 30 1 Z xy, sy, ,for alli (5) ,=I j=J G 7' yi, Z t >Xkjr , for all i, t and k c Si (6) j=1 ,=1 G Z8, =1, for allt (7) j= WO 2005/124104 PCT/AU2005/000761 - 64 xyt: , ,,for all i,j, and t (8) wherein; 5 T is the number of scheduling periods; N is the number of simulations; P is the total number of panels; G is the number of all possible cutoff grades; R," is the total rock in the panel i in simulations n. 10 Q, is the total ore in the panel i, simulation n, when mined with the COG j; Vo" is the value of the panel i, simulation n, when mined and processed with the COG j; R,' is the maximum mining capacity in period t; 15 Q,4 is the maximum processing rate in period t; S, is the set of panels that must be removed before starting the panel i; dt is the time discount factor; x, is the fraction of the panel i is extracted with the 20 COG j in period t; y, is a binary variable equal to 1 if the extraction of the panel i has started in periods 1 to t, and equal to 0 otherwise; 3, is a binary variable controlling the selection of the 25 COG applied in period t.