Print version ISSN 0038-223X
J. S. Afr. Inst. Min. Metall. vol.112 no.2 Johannesburg 2012
A new mathematical programming model for long-term production scheduling considering geological uncertainty
J. GholamnejadI; E. MoosaviII
IDepartment of Mining and Metallurgy Engineering, Yazd University, Yazd, Iran
IIDepartment of Mining Engineering, Islamic Azad University, Tehran, Iran
Determination of the optimum production schedules over the life of a mine is a critical mechanism in open pit mine planning procedures. Long-term production scheduling is used to maximize the net present value of the project under technical, financial, and environmental constraints. Mathematical programming models are well suited for optimizing long-term production schedules of open pit mines. There are two approaches to solving long-term production problems: deterministic- and uncertainty- based approaches. Deterministic-based models are unable to deal with grade and geological uncertainties, which are two important sources of risk in mining industries. This may lead to discrepancies between actual production obtained by these algorithms and planning expectations. In this paper, a new binary integer programming model was developed for long-term production scheduling that incorporates geological uncertainty within the orebody. Then, traditional and uncertainty-based models are applied to an iron ore deposit. Results showed that the uncertainty-based approach yields more practical schedules than traditional approaches in terms of production targets.
Keywords: open pit mine, long-term production scheduling, integer programming, geological uncertainty.
Long-term production scheduling involves the determination of sequences of ore and waste blocks to give the maximum net present value (NPV). This problem is subject to many aspects being satisfied, such as grade constraints, tonnage requirement for the plants, wall slope restriction, equipment capacities, etc. Long-term production scheduling determines the distribution of cash flow over life of mines. The NPV of each project depends on the grade and tonnage of the deposit, economic issues (operating and capital costs, and commodity price), and technical mining specifications (slope constraints, excavation capacities, etc.) In reality, some of the abovementioned issues can vary within certain limits, and the planner should make his decision on the production plan before knowing the exact values of the data. Among these uncertainties, those that are related to the orebody model lead to considerable deviation from production targets during the extraction process. Orebody model uncertainty has two important components1:
Geological uncertainty (tonnage uncertainty), which reflects the uncertainty related to the ore/waste contacts
Grade uncertainty, which reflects the error associated with the ore block grade estimation.
The geological and grade uncertainty cannot be eliminated; therefore, the best solution is to quantify uncertainties, reduce these uncertainties as far as investment permits, and finally manage the associated risk during the scheduling procedure. The latter can be achieved by explicit incorporation of uncertainty into the pit design and production scheduling process.
Many researchers have incorporated grade uncertainty into the optimization process during production scheduling2-8. However, geological uncertainty has not been considered yet during the scheduling process. In this paper, a binary integer programming formulation is developed for long-term production scheduling based on geological uncertainty to present an optimal solution under geological uncertainty. This formulation considers a probability index (PI) that reflects the probability that a block is an ore block. This index can be determined by the use of Indicator Kriging method. This model is then applied to an iron ore deposit.
Deterministic model for long-term production scheduling
Here the traditional integer programming formulation for long-term production scheduling is presented. Let:
|T||= maximum number of scheduling periods|
|t||= period in a scheduling span T|
|ijk||= indices correspond to the row, column, and level of blocks in the model|
|Γ||= the set of blocks|
|yijkt||= a binary variable that is equal to 1 if block ijk is mined in period t, otherwise 0|
|cijk||= the net value of block ijk|
|d1||= economic rate of interest|
|d2||= the risk discount rate|
|Zijk||= average grade of a block ijk|
|OTijk||= ore tonnage in the block ijk|
|WTijk||= tonnage of waste material within the block ijk|
|PU t||= maximum processing capacity in any period, t|
= minimum amount of ore that should be processed in period t
|MUt||= maximum available equipment capacity for each period t|
|MLt||= minimum amount of material (waste and ore) that should be removed in any period t|
|Ztmax||= maximum average grade of ore sent to the mill in any period t|
|Ztmin||= minimum average grade of ore sent to the mill in any period t|
|x||= index for 9 overlaying blocks.|
Then, the deterministic IP formulation's for long-term production scheduling can be written as follows:
Objective function [1) will maximize the NPV of a mining project. Constraints [2) and  that are grade blending constraints and control the feed quality sent to the mill. Constraints  and  are processing capacity constraints and ensure that milling capacities hold. Constraint  ensures that the equipment capacities hold. Constraint  can control the stripping ratio of ore. Constraints  are the wall slope restriction and constraint  ensures that the entire block in the block model should be removed in one period only.
In order to incorporate the geological uncertainty into this model, this kind of uncertainty should first, be quantified. In the next section, the Indicator Kriging method is introduced as a powerful tool to quantify the geological uncertainty.
Quantification of geological uncertainty by the use of Indicator Kriging
Indicator Kriging (IK) as a non-parametric technique in resource estimation is over fifteen years old. IK was introduced by Journel in 19839, and since then, despite the relative difficulty in its application, it has grown to become one of the most widely-applied grade estimation techniques in the minerals industry. Its appeal lies in the fact that it makes no assumptions about the distribution underlying the sample data, and, indeed, it can handle moderate mixing of diverse sample populations. However, despite the elegant and simple theoretical basis for IK, there are many practical implementation issues that affect its application and require serious consideration. These include aspects of order relations and their correction, the change of support, issues associated with highly skewed data, and the treatment of the extremes of the sample distribution when deriving estimates10. It is the main non-linear geostatistical technique used today in the minerals industry.
The essence of the indicator approach is the binomial coding of data into either 1 or 0 depending upon its relationship with a cut-off value, Zc. For a given value Z(x), ik (x) = 1 if Z(x) > Zc and ik(x) = 0 otherwise11. This is a non-linear transformation of the data value into either a 1 or a 0. Values that are much greater than a given cut-off, Zc, will receive the same indicator value as those values which are only slightly greater than that cut-off. Thus indicator transformation of data is an effective way of limiting the effect of very high values. It should be noted that this transformation can be also performed according to the presence or absence of a rock type (direct binary data requiring no transform), or a series of lithological or facies codes, or mineral sands hardness values. Simple or Ordinary Kriging of a set of indicator-transformed values will provide a resultant value between 0 and 1 for each point estimate, which can be interpreted either as12:
Proportions-the proportion of the block above the specified cut-off on data support
Probabilities-the probability that the grade is above the specified indicator. This probability is considered as the Probability Index (PIijk) for block ijk that is used in the optimization process in this paper. The blocks with those higher of probability have less risk than with the lower probability.
In the next section, this probability index is incorporated into the traditional deterministic model in order to obtain an uncertainty-based mathematical model for long-term production scheduling.
Production scheduling under geological uncertainty
In this section, an integer programming-based model is developed in order to taking account the geological uncertainty. In this approach, a probability is assigned to each block (PIik) which represents the probability that is produced from ik for each block in the block model. Now, we should set our objective function in a way that the blocks with higher certainty are mined in earlier production periods, leaving uncertain blocks for later periods, when additional information usually becomes available. Therefore, another objective function is added to the objective function of traditional model in the form of:
where d2 is the risk discount rate, which is used to discriminate the block extraction preference between time periods. If a higher rate is used, the differences in the probabilities between different periods are expected to be higher.
It can be seen that we have two objective functions  and  that should be optimized simultaneously. The process of optimizing systematically and simultaneously a collection of objective function is called multi-objective optimization. There are several approaches to solve multi objective optimization problems13, but the most common one is the weighted sum method. If we have M objective function Fl (x), l = 1,2,...,M, then the utility function (utility function is an amalgamation of the individual functions and is a mathematical expression that attempts to model the decisionmaker's preferences) can be expressed as:
Here, Wl are weights typically set by the decision-maker. These weights reflect the relative importance of each objective. There are several ways to select Wl which can be seen in Marler and Arora13. If all of the weights are positive, the maximum of  is Pareto-optimal14; i.e. maximizing  is sufficient for Pareto optimality. In addition, as it is clear, objective functions should be transformed such that they are dimensionless. One of the robust approaches to transform objective functions, regardless of their original range, is given as follows13:
: transformed objective function with a lower limit of unity
: utopia point (point that optimizes the objective function within the feasible design space). In this model we set W1 = W2 = 1, therefore the final objective function of this model can be written as:
Therefore Equation  can be rewritten as:
This objective function is subject to the constraints  to .
Application of deterministic and uncertainty based model in an iron ore deposit
A case study has been conducted on a central Iranian iron orebody to compare the results of the suggested algorithm with the deterministic one. The annual ore production rate of this deposit is 7 Mt of iron ore. The deposit contains about 317 Mt of ore with an average grade of 53% Fe and 1% P.
The block model contains 17 921 blocks with these dimensions 25 m × 25 m × 12.5 m. In order to determine the PI of each block, Indicator Kriging analysis was performed using Surpac 6.1.2 software. Cut-off grade is set at 25% Fe15.
Then the optimal pit limit has been determined using Lerches and Grossman method15 by the aid of Surpac16. This pit contains 232 Mt of ore with the stripping ratio of 1.33 and Fe and P average grades of 56.2% and .086%, respectively. Because there are a large number of blocks within the ultimate pit limit, the pit can be divided into a series of subpits commonly called push-backs, cut-backs, or phases. These push-backs are designed with haul road access and act as a guide during the yearly scheduling process. Five pushbacks are designed and the blocks within the first and second push-backs were submitted to the deterministic (optimization) formulation in Equation  and constraints in Equations  through ) and uncertainty-based algorithm (optimization) formulation in Equation  and constraints in Equations  through . This push-back, which contains 6770 blocks (4100 waste blocks and 2670 ore blocks) is depicted in Figure 1. Production schedule was fulfilled for a 5-year mine life. This study considers that ore material sent to the processing plant during each production period should have an average Fe grade of more than 52% and average P grade of less than 1%. d1 and d2 are assumed to be equal to 10%. The deterministic and suggested models were then solved by using the Risk Solver Platform V9.6 software18. Figure 2 shows the scheduling pattern generated from the proposed uncertainty-based model in the form of a bench plan. In order to compare the performance of these two algorithms, the number of blocks with the probability of more than 0.5 (PI>50%) is calculated in each period. This number is called the confidence index (CI). A low CI reflects high risk in achieving planned ore production, and a high CI a low risk.
Table I shows the summary results of production scheduling using traditional and uncertainty-based approaches.
As shown in Table I, the produced schedule using the uncertainty-based model has the highest CI in the first year (467), a lower CI in the second year (402) and the lowest in the last year (301). This means that the uncertainty-based model seeks the more certain areas of deposit in the first periods of operation and leaves less certain areas for the later periods, when additional information usually becomes available. Therefore, the probability of deviation from ore production target becomes lower by the use of the suggested model. The NPV (not mentioned in Table I) of the schedule generated by the traditional model is 9% higher than that of the suggested model, because the suggested model tends to maximize NPV and minimize the geological risk. This leads to some high-grade blocks are left for the latter periods. Therefore, it generates a realistic NPV, which is the best under conditions of uncertain geology. However, increasing the NPV will generally increase the risk of not meeting production targets.
In this paper, a geological uncertainty-based model is presented for open pit long-term production scheduling. The suggested integer programming model generates schedules that explicitly maximize the NPV of a project and reduce geological risk during the early production periods. Geological uncertainty was quantified by the use of the Indicator Kriging method. Using this method, a probability index can be obtained that reflects the probability that a block is an ore block. This index is then incorporated into the mathematical formulation for long-term production scheduling. The proposed model was applied on an iron ore deposit. Comparing the results from the proposed model with those from the deterministic one shows that in the proposed model, the total number of blocks with higher confidence in the first periods of exploitation is greater than that of the deterministic model while all the production constraints are satisfied. This can result in decreasing the risk of not meeting production targets during the first period of the mining process.
1. JACKSON, S., FREDERICKSEN, D., STEWART, M., VANN, J., BURKE, A., DUGDALE, J., and BERTOLI, O. Geological and grade risk at the Golden Gift and Magdala gold deposits, Stawell, Victoria, Australia. 5th International Mining Geology Conference, Bendigo, Victoria. Australasian Institute of Mining and Metallurgy, Melbourne, 17-19 November 2003. pp. 207-214. [ Links ]
2. DOWD, P.A. Risk in minerals projects: analysis, perception and management. Transaction of the Institution of Mining and Metallurgy A, vol. 106, January-April 1997. pp. A9-A18. [ Links ]
3. DENBY, B. and SCHOFIELD, D. Inclusion of risk assessment in open pit design and scheduling. Transaction of the Institution of Mining and Metallurgy A, vol. 104, January-April 1995. pp. A67-A71. [ Links ]
4. DIMITRAKOPOULOS, R. and RAMAZAN, S. Managing risk and waste mining in long-term production scheduling of open pit mine. 2003 SME Annual Meeting & Exhibit, Feb. 24-26, Cincinnati Ohio. Preprint: 03-151. [ Links ]
5. GODOY, M. and DIMITRAKOPOULOS, R. Managing risk and waste mining in long-term production scheduling of open pit mine. 2003 SME Annual Meeting & Exhibit, Feb. 24-26, Cincinnati Ohio. [ Links ]
6. RAMAZAN, S. and DIMITRAKOPOULOS, R. Traditional and new MIP models for production scheduling with in situ grade variabulity. International Journal of Surface Mining, Reclamation and Environment, vol. 18, no. 2, 2003. pp. 85-98. [ Links ]
7. GHOLAMNEJAD, J., OSANLOO, M., and KARIMI, B. A chance-constrained programming approach for open pit long-term production scheduling in stochastic environments. Journal of the Southern African Institute o f Mining and Metallurgy, vol. 106, 2006. pp. 105-114. [ Links ]
8. GHOLAMNEJAD, J. and OSANLOO, M. Using chance constrained binary integer programming in optimizing long term production scheduling for open pit mine design. Mining Technology: IMM Transactions section A, vol. 116, no. 2, June 2007, pp. 58-66. [ Links ]
9. JOURNEL, A.G. Non-parametric estimation of spatial distributions, Mathematical Geology, vol. 15, no. 3, 1983. pp. 445-468. [ Links ]
10. GLACKEN, I.M. and BLACKNEY, P.C.J. A practitioners implementation of indicator Kriging. Beyond Ordinary Kriging Seminar 30 October, 1998, Perth, Western Australia. The Geostatistical Association of Australia. [ Links ]
11. DIMITRAKOPOULOS, R. Indicator Kriging course notes, W.H. Bryan Centre, University of Queensland, 1997. [ Links ]
12. ISAAKS, E.H. and SRIVASTAVA, R.M. An introduction to applied geostatistics. Oxford University Press, New York. 1989. [ Links ]
13. MARLER, R.T. and ARORA, J.S. Survey of multi-objective optimization methods for engineering, Structural and Multidisciplinary Optimization, vol. 26, no. 6, April 2004, pp. 369-395. [ Links ]
14. PARETO, V. Manual of Political Economy. Schwier, A.S. and Page, A.N. (eds.). A.M. Kelley, New York, 1971. [ Links ]
15. LERCHS, H. and GROSSMAN, F. Optimum design of open-pit mines. CIM Transactions, vol. 58, no. 633, 1965. pp. 47-54. [ Links ]
16. Surpac Users Manual, Version V5.2-D. Surpac Minex Group. 2006. [ Links ]
17. GHOLAMNEJAD, J. and OSANLOO, M. Incorporation of ore grade uncertainty into the push back design process, Journal of the Southern African Institute of Mining and Metallurgy, vol. 107, 2007. pp. 177-185. [ Links ]
Paper received Aug. 2010; revised paper received Sep. 2011.
© The Southern African Institute of Mining and Metallurgy, 2012. SA ISSN 0038-223X/3.00 + 0.00.