Document Type : Original Research Paper


1 Faculty of Mining, Petroleum & Geophysics Engineering, Shahrood University of Technology, Shahrood, Iran

2 School of Mining and Petroleum Engineering, University of Alberta, Alberta, Canada


In mining projects, all uncertainties associated with a project must be considered to determine the feasibility study. Grade uncertainty is one of the major components of technical uncertainty that affects the variability of the project. Geostatistical simulation, as a reliable approach, is the most widely used method to quantify risk analysis to overcome the drawbacks of the estimation methods used for an entire ore body. In this work, all the algorithms developed by numerous researchers for optimization of the underground stope layout are reviewed. After that, a computer program called stope layout optimizer 3D is developed based on a previously proposed heuristic algorithm in order to incorporate the influence of grade variability in the final stope layout. Utilizing the sequential gaussian conditional simulation, 50 simulations and a kriging model are constructed for an underground copper vein deposit situated in the southwest of Iran, and the final stope layout is carried out separately. It can be observed that geostatistical simulation can effectively cope with the weakness of the kriging model. The final results obtained show that the frequency of economic value for all realizations varies between 6.7 M$ and 30.7 M$. This range of variation helps designers to make a better and lower risk decision under different conditions.


Main Subjects

[1]. Ross, J.G. (2004). Risk and uncertainty in portfolio characterization. Journal of Petroleum Science and Engineering. 44(1): 41-53.
[2]. Dominy, S.C. Annels, A.E. and Noppe, M. (2002). Errors and uncertainty in Ore Reserve estimates-Operator beware. In Proceedings, Underground Operators Conference. (pp. 121-126).
[3]. Morley, D. (1999). Financial impact of resource/reserve uncertainty. Journal of the Southern African Institute of Mining and Metallurgy. 99(6): 293-301.
[4]. Dimitrakopoulos, R. (1990). Conditional simulation of intrinsic random functions of order k. Mathematical geology. 22(3): 361-380.
[5]. Armstrong, M. and Dowd, P.A.(Eds.). (2013). Geostatistical Simulations: Proceedings of the Geostatistical Simulation Workshop, Fontainebleau, France. 27–28 May 1993 (Vol. 7). Springer Science & Business Media.
[6]. Smith, M. and Dimitrakopoulos, R. (1999). The influence of deposit uncertainty on mine production scheduling. International Journal of Surface Mining, Reclamation and Environment. 13(4): 173-178.
[7]. Kumral, M. and Dowd, P.A. (2001). Short-term scheduling for industrial minerals using multi-objective simulated annealing, APCOM 2001, Phoenix, Arizona.
[8]. Ramazan, S. and Dimitrakopoulos, R. (2007). Stochastic optimization of long-term production scheduling for open pit mines with a new integer programming formulation. Ore-body Modelling and Strategic Mine Planning, Spectrum Series. 14: 359-365.
[9]. Koushavand, B. Askari-Nasab, H. and Deutsch, C.V. (2014). A linear programming model for long-term mine planning in the presence of grade uncertainty and a stockpile. International Journal of Mining Science and Technology. 24(4): 451-459.
[10]. Dimitrakopoulos, R. Farrelly, C.T. and Godoy, M. (2002). Moving forward from traditional optimization: grade uncertainty and risk effects in open-pit design. Mining Technology. 111(1): 82-88.
[11]. Dimitrakopoulos, R. Martinez, L. and Ramazan, S. (2007). A maximum upside/minimum downside approach to the traditional optimization of open pit mine design. Journal of Mining Science. 43(1): 73-82.
[12]. Leite, A. and Dimitrakopoulos, R. (2007). Stochastic optimization model for open pit mine planning: application and risk analysis at copper deposit. Mining Technology. 116(3): 109-118.
[13]. Montiel, L. Dimitrakopoulos, R. and Kawahata, K. (2016). Globally optimizing open-pit and underground mining operations under geological uncertainty. Mining Technology. 125(1): 2-14.
[14]. Grieco, N. and Dimitrakopoulos, R. (2007). Managing grade risk in stope design optimization: probabilistic mathematical programming model and application in sublevel stoping. Mining technology. 11
[15]. Dimitrakopoulos, R. and Grieco, N. (2009). Stope design and geological uncertainty: quantification of risk in conventional designs and a probabilistic alternative. Journal of Mining Science. 45(2): 152-163.
[16]. Riddle, J.M. (1977). A dynamic programming solution of a block-caving mine layout. 14th International Symposium on Application of Computers and Operations Research in the Minerals Industries. Society for Mining, Metallurgy and Exploration (pp. 767-780).
[17]. Deraisme, J. De Fouquet, C. and Fraisse, H. (1984). Geostatistical ore-body model for computer optimization of profits from different underground mining methods. 18th APCOM IMM.
[18]. Cheimanoff, N.M. Deliac, E.P. and Mallet, J.L. (1989). GEOCAD: an alternative CAD and artificial intelligence tool that helps moving from geological resources to mineable reserves. In 21st International Symposium on the Application of Computers and Operations Research in the Mineral Industry. (pp. 471-478).
[19]. GEOCAD, (2017).
[20]. Ataee-Pour, M. (2005). A critical survey of the existing stope layout optimization techniques. Journal of Mining Science. 41(5): 447-466.
[21]. Alford, C, (1995), Optimization in underground mine design. Proceedings of the 25th International APCOM Symposium. The Australasian Institute of Mining and Metallurgy, Melbourne.
[22]. CAE Studio 3. (2014).
[23]. Ovanic, J. and Young, D.S. (1999). Economic optimization of open stope geometry. In 28th international APCOM symposium. Colorado school of Mines, Golden, Colorado. USA (pp. 855-862).
[24]. LINGO and optimization modeling, (2017).
[25]. GAMS Development Corporation, (2017). Washington DC, USA.
[26]. ILOG CPLEX's mathematical optimization technology, (2017). IBM Corporation
[27]. Mirzaeian, Y. and Ataee-pour, M. (2011). Optimization of Stope Geometry Using Piecewise Linear Function and MIP Approach. Journal of Civil & Environmental Engineering and Science Technology. 43(1): pp.79–87. (in Persian).
[28]. MPS, Ketron Optimization mathematical programming system, (2011).
[29]. Ataee-pour, M. (2000). A Heuristic Algorithm to Optimize Stope Boundaries, Ph.D. Thesis, University of Wollongong, Australia.
[30]. Ataee-Pour, M. (2004). Optimization of stope limits using a heuristic approach. Mining Technology. 113(2): 123-128.
[31]. Carwse, I. (2001). multiple pass floating stope process, In Proceedings of the Fourth Biennial Conference: Strategic Mine Planning. Melbourne: Australasian Institute of Mining and Metallurgy.
[32]. Jalali, S.E. and Ataee-pour, M. (2004). A 2D Dynamic Programming Algorithm to Optimize Stope Boundaries, Proceedings of the 13th Symposium on Mine Planning and Equipment Selection, Rotterdam, Balkema. pp.45–52.
[33]. Jalali, S.E. Ataee-pour, M. and Shahriar, K. (2007). Rigorous Algorithms to Optimize Stope Boundaries; Capabilities, Restrictions and Applications, Modern Management of Mine Producing, Geology and Environmental Protection. Albena, Bulgaria.
[34]. Jalali, S.E. Ataee-pour, M. Shahriyar, K. and Elahi-Zeyni, E. (2007). A Computer Program to Optimize Stope Boundaries Using Probable Stope Algorithm. Iranian Journal of Mining Engineering (IRJME). 2(3): pp.7–14, (in Persian).
[35]. Jalali, S.E. Ataee-pour, M. Shahriar, K. Elahi-Zeyni, E. and Nikbin, V. (2016). Computer Based Optimization of Underground Mining Area. Journal of Civil & Environmental Engineering and Science Technology. 48(4): pp.475–489. (in Persian).
[36]. Grieco, N. and Dimitrakopoulos, R. (2007). Managing grade risk in stope design optimization: probabilistic mathematical programming model and application in sublevel stoping. Mining technology. 116(2): 49-57.
[37]. Dimitrakopoulos, R. and Grieco, N. (2009). Stope design and geological uncertainty: quantification of risk in conventional designs and a probabilistic alternative. Journal of Mining Science. 45(2): 152-163.
[38]. Jalali, S.E. and Hosseini, H. (2009). Optimization of Extraction Range in Underground Mining Using a Greedy Algorithm. Journal of Science and Research in Mining Engineering. 4(9): pp.1–11, (in Persian).
[39]. Cormen, T.H. Leiserson, C.E. Rivest, R.L. and Stein, C. (2001). Greedy algorithms. Introduction to algorithms. 1: 329-355.
[40]. Topal, E. and Sens, J. (2010). A new algorithm for stope boundary optimization. Journal of Coal Science and Engineering (China). 16(2): 113-119.
[41]. Bai, X. Marcotte, D. and Simon, R. (2013). Underground stope optimization with network flow method. Computers & Geosciences. 52: 361-371.
[42]. Sandanayake, D.S.S. Topal, E. and Asad, M.W.A. (2015). A heuristic approach to optimal design of an underground mine stope layout. Applied Soft Computing. 30: 595-603.
[43]. Sandanayake, D.S.S., Topal, E. and Asad, M.W.A. (2015). Designing an optimal stope layout for underground mining based on a heuristic algorithm. International Journal of Mining Science and Technology. 25(5): 767-772.
[44]. Nikbin, V., Ataee-pour, M., Shahriar, K. and Pourrahimian, Y. (2018). A 3D approximate hybrid algorithm for stope boundary optimization. Computers & Operations Research.
[45]. Visual Studio, (2013).
[46]. Sotoudeh, F. Kakaie, R. and Ataei, M. (2017). 'Development of a computer program for underground mine stope optimisation using a heuristic algorithm', in M Hudyma & Y Potvin (eds). Proceedings of the First International Conference on Underground Mining Technology. Australian Centre for Geomechanics. Perth. pp. 689-700.
[47]. Manchuk, J. and Deutsch, C. (2008). Optimizing stope designs and sequences in underground mines. SME Transactions. 324. 67-75.
[48]. Ashgari, O. and Esfahani, N.M. (2013). A new approach for the geological risk evaluation of coal resources through a geostatistical simulation. Arabian Journal of Geosciences. 6(3): 929-943.
[49]. Yunsel, T.Y. (2012). Risk quantification in grade variability of gold deposits using sequential Gaussian simulation. Journal of Central South University. 19(11): 3244-3255.
[50]. Dimitrakopoulos, R. (1998). Conditional simulation algorithms for modelling ore body uncertainty in open pit optimization. International Journal of Surface Mining, Reclamation and Environment. 12(4): 173-179.
[51] Asghari, O. Soltni, F. and Amnieh, H.B. (2009). The comparison between sequential gaussian simulation (SGS) of Choghart ore deposit and geostatistical estimation through ordinary kriging. Australian Journal of Basic and Applied Sciences. 3(1): 330-341.
[52]. Soltani, F. Afzal, P. and Asghari, O. (2013). Sequential Gaussian simulation in the Sungun Cu porphyry deposit and comparing the stationary reproduction with ordinary kriging. Universal Journal of Geoscience. 1(2): 106-113.
[53]. Sotoudeh, F. Ataei, M. Kakaie, R. (2016). Risk Assessment of Chehlkoreh Copper Deposit Based on Grade Variability Using Geostatistical Simulation; 4th International Mine and Mining Industries Congress, Tehran, Iran.
[54]. Rezaee, H. Asghari, O. Koneshloo, M. and Ortiz, J. M. (2014). Multiple-point geostatistical simulation of dykes: application at Sungun porphyry copper system, Iran. Stochastic environmental research and risk assessment. 28(7): 1913-1927.
[55]. Safikhani, M. Asghari, O. and Emery, X. (2016). Assessing the accuracy of Sequential Gaussian Simulation through statistical testing. Stochastic Environmental Research and Risk Assessment, 1-11.
[56]. Emery, X. (2004). Testing the correctness of the sequential algorithm for simulating Gaussian random fields. Stochastic Environmental Research and Risk Assessment. 18(6): 401-413.
[57]. Webster, R. and Oliver, M.A. (2007). Geostatistics for environmental scientists. John Wiley and Sons.
[58]. SGEMS: The Stanford Geostatistical Earth Modeling Software, (2008).
[59]. Emery, X. and Peláez, M. (2011). Assessing the accuracy of sequential Gaussian simulation and cosimulation. Computational Geosciences. 15(4): 673-689.
[60] Tatiya, R.R. (2005). Classification Stoping Methods. In Surface and Underground Excavations’, London: Taylor & Francis Group. 415–431.
[61] Villaescusa, E. (2000). A review of sublevel stoping. MassMin. 577–590