Document Type : Original Research Paper


1 Department of Mining Engineering, Bandung Institute of Technology, Indonesia

2 School of Mines China University of Mining and Technology, P.R.China

3 Baluchistan University of IT, Engineering and Management Sciences, Pakistan

4 Department of Mining Engineering, University of Engineering and Technology, Peshawar, Pakistan


This paper elucidates a new idea and concept for exploration of the gold ore deposits.  The cyanidation method is traditionally used for gold extraction. However, this method is laborious, time-consuming, costly, and depends upon the availability of the processing units. In this work, an attempt is made in order to update the gold exploration method by the Monte Carlo-based simulation. An excellent approach always requires a high quality of the datasets for a good model. A total of 48 incomplete datasets are collected from the Shoghore district, Chitral area of Khyber, Pakhtunkhwa, Pakistan. The cyanidation leaching test is carried out in order to measure the percentage of the gold ore deposits. In this work, the mean, median, mode, and successive iteration substitute methods are employed in such a way that they can compute the datasets with missing attributes. The multiple regression analysis is used to find a correlation between the potential of hydrogen ion concentration (pH), solid content (in %), NaCN concentration (in ppm), leaching time (in Hr), particle size (in µm), and measured percentage of gold recovery (in %). Moreover, the normal Archimedes and exponential distributions are employed in order to forecast the uncertainty in the measured gold ore deposits. The performance of the model reveals that the Monte Carlo approach is more authentic for the probability estimation of gold ore recovery. The sensitivity analysis reveals that pH is the most influential parameter in the estimation of the gold ore deposits. This stochastic approach can be considered as a foundation to foretell the probabilistic exploration of the new gold deposits.


[1]. Ellis, L., Hlongwa, N., Kemp, H., Kershoff, G., Lemboe, C., Smit, B. and Van der Wath, N. (2011) First Quarter.
[2]. Martens, E., Zhang, H., Prommer, H., Greskowiak, J., Jeffrey, M. and Roberts, P. (2012). In situ recovery of gold: Column leaching experiments and reactive transport modeling. Hydrometallurgy, 125: 16-23.
[3]. Ramesh, A., Hasegawa, H., Sugimoto, W., Maki, T. and Ueda, K. (2008). Adsorption of gold (III), platinum (IV), and palladium (II) onto glycine modified crosslinked chitosan resin. Bioresource technology, 99 (9): 3801-3809.
[4]. Wang, H. and Ren, Z.J. (2014). Bio-electrochemical metal recovery from wastewater: a review. Water research, 66, 219-232.
[5]. Schodde, R.C. (2004). Discovery performance of the western world gold industry over the period 1985–2003. In Proceedings of PACRIM 2004 Congress.
[6]. Vearncombe, J.R. and Phillips, G.N. (2020). The importance of brownfields gold exploration. Mineralium Deposita, 55 (2): 189-196.
[7]. Chapman, R.J., Leake, R.C., Bond, D.P.G., Stedra, V. and Fairgrieve, B. (2009). Chemical and mineralogical signatures of gold formed in oxidizing chloride hydrothermal systems and their significance within populations of placer gold grains collected during reconnaissance. Economic Geology, 104 (4): 563-585.
[8]. Townley, B.K., Hérail, G., Maksaev, V., Palacios, C., de Parseval, P., Sepulveda, F. and Ulloa, C. (2003). Gold grain morphology and composition as an exploration tool: application to gold exploration in covered areas. Geochemistry: Exploration, Environment, Analysis, 3 (1): 29-38.
[9]. Chapman, R.J., Mortensen, J.K. and LeBarge, W.P. (2011). Styles of lode gold mineralization contributing to the placers of the Indian River and Black Hills Creek, Yukon Territory, Canada as deduced from micro-chemical characterization of placer gold grains. Mineralium Deposita, 46 (8): 881-903.
[10]. Lalomov, A.V., Chefranov, R.M., Naumov, V.A., Naumova, O.B., LeBarge, W. and Dilly, R.A. (2017). Typomorphic features of placer gold of Vagran cluster (the Northern Urals) and search indicators for primary bedrock gold deposits. Ore Geology Reviews, 85: 321-335.
[11]. Moles, N.R., Chapman, R.J. and Warner, R.B. (2013). The significance of copper concentrations in natural gold alloy for reconnaissance exploration and understanding gold-depositing hydrothermal systems. Geochemistry: Exploration, Environment, Analysis, 13 (2): 115-130.
[12]. Omang, B.O., Suh, C.E., Lehmann, B., Vishiti, A., Chombong, N.N., Fon, A.N. and Shemang, E.M. (2015). Micro-chemical signature of alluvial gold from two contrasting terrains in Cameroon. Journal of African Earth Sciences, 112: 1-14.
[13]. Youngson, J.H. and Craw, D. (1993). Gold nugget growth during tectonically induced sedimentary recycling, Otago, New Zealand. Sedimentary geology, 84 (1-4): 71-88.
[14]. Jankovic, S. (1984). Metallogeny and mineral potential of Northern Pakistan (north of Indus Suture zone). A preliminary assessment.
[15]. Halfpenny, R. and Mazzucchelli, R.H. (1999). Regional multi-element drainage geochemistry in the Himalayan Mountains, northern Pakistan. Journal of Geochemical Exploration, 67 (1-3): 223-233.
[16]. Ali, L., Williamson, B.J., Moon, C.J., Shah, M.T. and Khattak, S.A. (2015). Distribution of gold in different size fractions of stream sediments as a guide to bedrock gold mineralization along the Shyok Suture Zone and adjacent areas of the Kohistan Island Arc, Pakistan. Arabian Journal of Geosciences, 8 (4): 2227-2235.
[17]. Miandad, S., Shah, M.T., Khan, S.D. and Ahmad, L. (2014). Investigation for gold and base metals mineralization and petrochemical characteristics of the rocks of upper parts of Bagrot valley, Gilgit-Baltistan, Pakistan. Journal of Himalayan Earth Sciences, 47 (2): 29.
[18]. Shah, M.T., Khan, S.D. and Ahmad, L. (2015). Investigation of gold and base metal mineralization and the petrochemical studies of the associated host rocks, Shagari Bala area, Skardu Gilgit-Baltistan, Pakistan. Journal of Himalayan Earth Science, 48 (1).
[19]. Alam, M., Li, S.R., Santosh, M. and Yuan, M.W. (2019). Morphology and chemistry of placer gold in the Bagrote and Dainter streams, northern Pakistan: Implications for provenance and exploration. Geological Journal, 54 (3): 1672-1687.
[20]. Alam, M., Li, S.R., Santosh, M., Shah, A., Yuan, M.W., Khan, H. and Zeng, Y.J. (2019). Morphological, thermo-electrical, geochemical and isotopic anatomy of auriferous pyrite from the Bagrote valley placer deposits, North Pakistan: Implications for ore genesis and gold exploration. Ore Geology Reviews, 112: 103008.
[21]. Nasri, S., Kalate, A.N., Kahoo, A.R. and Monfared, M.S. (2020). New insights into the structural model of the Makran subduction zone by fusion of 3D inverted geophysical models. Journal of Asian Earth Sciences, 188: 104075.
[22]. Radfar, A., Chakdel, A.R., Nejati, A. and Soleimani, M. (2019). New insights into the structure of the South Caspian Basin from seismic reflection data, Gorgan Plain, Iran. International Journal of Earth Sciences, 108 (2): 379-402.
[23]. Idris, M.A. (2014). Probabilistic stability analysis of underground mine excavations (Doctoral dissertation, Luleå Tekniska Universitet).
[24]. Wattimena, R.K., Kramadibrata, S., Sidi, I.D., Arif, I. and Azizi, M.A. (2012, January). Probabilistic analysis of single bench using new slope stability curves. In ISRM Regional Symposium-7th Asian Rock Mechanics Symposium. International Society for Rock Mechanics and Rock Engineering.
[25]. Hajsadeghi, S., Asghari, O., Mirmohammadi, M., Afzal, P., and Meshkani, S.A. (2020). Uncertainty-volume fractal model for delineating copper mineralization controllers using geostatistical simulation in Nohkouhi volcanogenic massive sulfide deposit, Central Iran. Maden Tetkik ve Arama Dergisi, 161 (161): 1-11.
[26]. Maleki, M. and Madani, N. (2016). Multivariate geostatistical analysis: An application to ore body evaluation. Iranian Journal of Earth Sciences, 8 (2): 173-184.
[27]. Aslam, B., Maqsoom, A., Alaloul, W.S., Musarat, M.A., Jabbar, T., and Zafar, A. (2020). Soil erosion susceptibility mapping using a GIS-based multi-criteria decision approach: Case of district Chitral, Pakistan. Ain Shams Engineering Journal.
[28]. Azizi, A., Ghaedrahmati, R., Ghahramani, N. and Rooki, R. (2016). Modelling and simulation of the cyanidation process of Aghdareh gold ore using artificial neural network and multiple linear regression. International Journal of Mining and Mineral Engineering, 7 (2): 139-154.
[29]. Yilmaz, T., Seckin, G. and Yuceer, A. (2010). Modeling of effluent COD in UAF reactor treating cyanide containing wastewater using artificial neural network approaches. Advances in Engineering Software, 41 (7-8): 1005-1010.
[30]. Zhang, J., Mao, Z.Z. and He, D.K. (2015). Real-time optimization based on a serial hybrid model for gold cyanidation leaching process. Minerals Engineering, 70: 250-263.
[31]. Jun, Z., Hua, Y., Hongxia, Y., Zhongda, T. and Runda, J. (2019). Gold recovery modeling based on interval prediction for a gold cyanidation leaching plant. IEEE Access, 7, 71511-71528. [32Paulin, F. and Santhakumaran, A. (2010). An Algorithm to Reconstruct the Missing Values for Diagnosing the Breast Cancer. Global Journal of Computer Science and Technology.
[33]. Twala, B., Cartwright, M. and Shepperd, M. (2005, November). Comparison of various methods for handling incomplete data in software engineering databases. In 2005 International Symposium on Empirical Software Engineering, 2005. (pp. 10-pp). IEEE.
[34]. Khotimah, B.K., Miswanto, M. and Suprajitno, H. (2020). Optimization of feature selection using genetic algorithm in naïve Bayes classification for incomplete data. Int. J. Intell. Eng. Syst, 13 (1): 334-343.
[35]. Valdiviezo, H.C. and Van Aelst, S. (2015). Tree-based prediction on incomplete data using imputation or surrogate decisions. Information Sciences, 311: 163-181.
[36]. Zhang, S., Jin, Z. and Zhu, X. (2011). Missing data imputation by utilizing information within incomplete instances. Journal of Systems and Software, 84 (3): 452-459.
[37]. Yuen, K.V., Beck, J.L. and Katafygiotis, L.S. (2006). Efficient model updating and health monitoring methodology using incomplete modal data without mode matching. Structural Control and Health Monitoring: The Official Journal of the International Association for Structural Control and Monitoring and of the European Association for the Control of Structures, 13 (1): 91-107.
[38]. Liu, F. (2011). Direct mode-shape expansion of a spatially incomplete measured mode by a hybrid-vector modification. Journal of Sound and Vibration, 330 (18-19): 4633-4645.
[39]. Gromski, P.S., Xu, Y., Kotze, H.L., Correa, E., Ellis, D.I., Armitage, E.G. and Goodacre, R. (2014). Influence of missing values substitutes on multivariate analysis of metabolomics data. Metabolites, 4 (2): 433-452.
[40]. Soleimani, M. and Shokri, B.J. (2016). Intrinsic geological model generation for chromite pods in the Sabzevar ophiolite complex, NE Iran. Ore geology reviews, 78: 138-150.
[41]. Rajesh, G. and Sunoj, S.M. (2019). Some properties of cumulative Tsallis entropy of order $$\alpha $$ α. Statistical Papers, 60 (3): 933-943.
[42]. Karkevandi-Talkhooncheh, A., Sharifi, M. and Ahmadi, M. (2018). Application of hybrid adaptive neuro-fuzzy inference system in well placement optimization. Journal of Petroleum Science and Engineering, 166: 924-947.
[43]. Koopialipoor, M., Nikouei, S.S., Marto, A., Fahimifar, A., Armaghani, D.J. and Mohamad, E.T. (2019). Predicting tunnel boring machine performance through a new model based on the group method of data handling. Bulletin of Engineering Geology and the Environment, 78 (5): 3799-3813.
[44]. Aghajani-Bazzazi, A., Osanloo, M. and Azimi, Y. (2009, September). Flyrock prediction by multiple regression analysis in Esfordi phosphate mine of Iran. In Proceedings of the 9th international symposium on rock fragmentation by blasting. Granada, Spain (pp. 649-57).
[45]. Wang, Y., Mazin, H.E., Xu, W. and Huang, B. (2016). Estimating harmonic impact of individual loads using multiple linear regression analysis. International Transactions on Electrical Energy Systems, 26 (4): 809-824.
[46]. Liang, W., Sari, A., Zhao, G., McKinnon, S.D. and Wu, H. (2020). Short-term rock-burst risk prediction using ensemble learning methods. Natural Hazards, 104 (2): 1923-1946.
[47]. Wattimena, R.K., Kramadibrata, S., Sidi, I.D., Arif, I. and Azizi, M.A. (2012, January). Probabilistic analysis of single bench using new slope stability curves. In ISRM Regional Symposium-7th Asian Rock Mechanics Symposium. International Society for Rock Mechanics and Rock Engineering.
[48]. Azizi, M.A., Kramadibrata, S., Wattimena, R.K. and Sidi, I.D. (2013). Probabilistic analysis of physical models slope failure. Procedia Earth and Planetary Science, 6: 411-418.
[49]. Solver, F. (2010). Premium solver platform. User Guide, Frontline Systems.
[50]. EPA, U. (1997). Environmental Protection Agency. Guiding principles for Monte Carlo analysis. EPA/630/R-97/001.
[51]. Ghasemi, E., Sari, M. and Ataei, M. (2012). Development of an empirical model for predicting the effects of controllable blasting parameters on fly-rock distance in surface mines. International Journal of Rock Mechanics and Mining Sciences, 52: 163-170.
[52]. Mahdiyar, A., Hasanipanah, M., Armaghani, D.J., Gordan, B., Abdullah, A., Arab, H. and Abd Majid, M.Z. (2017). A Monte Carlo technique in safety assessment of slope under seismic condition. Engineering with Computers, 33 (4): 807-817.
[53]. Bianchini, F. and Hewage, K. (2012). Probabilistic social cost-benefit analysis for green roofs: A lifecycle approach. Building and Environment, 58, 152-162.
[54]. Lima, H.R.B.R., Nascimento, D.S., Bispo, G.F.C., Teixeira, V.C., Valério, M.E.G. and Souza, S.O. (2014). Production and characterization of spodumene dosimetric pellets prepared by a sol–gel route. Radiation Physics and Chemistry, 104: 93-99.
[55]. Morin, M.A. and Ficarazzo, F. (2006). Monte Carlo simulation as a tool to predict blasting fragmentation based on the Kuz–Ram model. Computers & geosciences, 32 (3): 352-359.
[56]. Kamran, M. (2021). A Probabilistic Approach for Prediction of Drilling Rate Index using Ensemble Learning Technique. Journal of Mining and Environment.
[57]. Wu, S., Wu, Z. and Zhang, C. (2019). Rock-burst prediction probability model based on case analysis. Tunneling and Underground Space Technology, 93: 103069.