Grade and Tonnage Uncertainty Analysis of an African Copper Deposit Using Multiple-Point Geostatistics and Sequential Gaussian Simulation

Spatial uncertainty analysis is a complex and difficult task for orebody estimation in the mining industry. Conventional models (kriging and its variants) with variogram-based statistics fail to capture the spatial complexity of an orebody. Due to this, the grade and tonnage are incorrectly estimated resulting in inaccurate mine plans, which lead to costly financial decision. Multiple-point geostatistical simulation model can overcome the limitations of the conventional two-point spatial models. In this study, a multiple-point geostatistical method, namely SNESIM, was applied to generate multiple equiprobable orebody models for a copper deposit in Africa, and it helped to analyze the uncertainty of ore tonnage of the deposit. The grade uncertainty was evaluated by sequential Gaussian simulation within each equiprobable orebody models. The results were validated by reproducing the marginal distribution and two- and three-point statistics. The results show that deviations of volume of the simulated orebody models vary from ? 3 to 5% compared to the training image. The grade simulation results demonstrated that the average grades from the different simulation are varied from 3.77 to 4.92% and average grade 4.33%. The results also show that the volume and grade uncertainty model overestimates the orebody volume as compared to the conventional orebody. This study demonstrates that incorporating grade and volume uncertainty leads to significant changes in resource estimates.     

An Analysis of the Published Mineral Resource Estimates of the Haji-Gak Iron Deposit,Afghanistan

The Haji-Gak iron deposit of eastern Bamyan Province, eastern Afghanistan, was studied extensively and resource calculations were made in the 1960s by Afghan and Russian geologists. Recalculation of the resource estimates verifies the original estimates for categories A (in-place resources known in detail), B (in-place resources known in moderate detail), and C₁ (in-place resources estimated on sparse data), totaling 110.8 Mt, or about 6% of the resources as being supportable for the methods used in the 1960s. C₂ (based on a loose exploration grid with little data) resources are based on one ore grade from one drill hole, and P₂ (prognosis) resources are based on field observations, field measurements, and an ore grade derived from averaging grades from three better sampled ore bodies. C₂ and P₂ resources are 1,659.1 Mt or about 94% of the total resources in the deposit. The vast P₂ resources have not been drilled or sampled to confirm their extent or quality. The purpose of this article is to independently evaluate the resources of the Haji-Gak iron deposit by using the available geologic and mineral resource information including geologic maps and cross sections, sampling data, and the analog-estimating techniques of the 1960s to determine the size and tenor of the deposit.     

Geological Modeling of Gold Deposit Based on Grade Domaining Using Plurigaussian Simulation Technique

Mineral resource evaluation requires defining grade domains of an ore deposit. Common practice in mineral resource estimation consists of partitioning the ore body into several grade domains before the geostatistical modeling and estimation at unsampled locations. Many ore deposits are made up of different mineralogical ensembles such as oxide and sulfide zone: being able to model the spatial layout of the different grades is vital to good mine planning and management. This study addresses the application of the plurigaussian simulation to Sivas (Turkey) gold deposits for constructing grade domain models that reproduce the contacts between different grade domains in accordance with geologist’s interpretation. The method is based on the relationship between indicator variables from grade distributions on the Gaussian random functions chosen to represent them. Geological knowledge is incorporated into the model by the definition of the indicator variables, their truncation strategy, and the grade domain proportions. The advantages of the plurigaussian simulation are exhibited through the case study. The results indicated that the processes are seen to respect reproducing complex geometrical grades of an ore deposit by means of simulating several grade domains with different spatial structure and taking into account their global proportions. The proposed proportion model proves as simple to use in resource estimation, to account for spatial variations of the grade characteristics and their distribution across the studied area, and for the uncertainty in the grade domain proportions. The simulated models can also be incorporated into mine planning and scheduling.     

Monte Carlo Simulations of Product Distributions and Contained Metal Estimates

Estimation of product distributions of two factors was simulated by conventional Monte Carlo techniques using factor distributions that were independent (uncorrelated). Several simulations using uniform distributions of factors show that the product distribution has a central peak approximately centered at the product of the medians of the factor distributions. Factor distributions that are peaked, such as Gaussian (normal) produce an even more peaked product distribution. Piecewise analytic solutions can be obtained for independent factor distributions and yield insight into the properties of the product distribution. As an example, porphyry copper grades and tonnages are now available in at least one public database and their distributions were analyzed. Although both grade and tonnage can be approximated with lognormal distributions, they are not exactly fit by them. The grade shows some nonlinear correlation with tonnage for the published database. Sampling by deposit from available databases of grade, tonnage, and geological details of each deposit specifies both grade and tonnage for that deposit. Any correlation between grade and tonnage is then preserved and the observed distribution of grades and tonnages can be used with no assumption of distribution form.     

Basic concepts in three-part quantitative assessments of undiscovered mineral resources

Since 1975, mineral resource assessments have been made for over 27 areas covering 5×10⁶ km² at various scales using what is now called the three-part form of quantitative assessment. In these assessments, (1) areas are delineated according to the types of deposits permitted by the geology,(2) the amount of metal and some ore characteristics are estimated using grade and tonnage models, and (3) the number of undiscovered deposits of each type is estimated.Permissive boundaries are drawn for one or more deposit types such that the probability of a deposit lying outside the boundary is negligible, that is, less than 1 in 100,000 to 1,000,000.     

Non-stationary Geostatistical Modeling: A Case Study Comparing LVA Estimation Frameworks

Incorporating locally varying anisotropy (LVA) in geostatistical modeling improves estimates for structurally complex domains where a single set of anisotropic parameters modeled globally do not account for all geological features. In this work, the properties of two LVA-geostatistical modeling frameworks are explored through application to a complexly folded gold deposit in Ghana. The inference of necessary parameters is a significant requirement of geostatistical modeling with LVA; this work focuses on the case where LVA orientations, derived from expert geological interpretation, are used to improve the grade estimates. The different methodologies for inferring the required parameters in this context are explored. The results of considering different estimation frameworks and alternate methods of parameterization are evaluated with a cross-validation study, as well as visual inspection of grade continuity along select cross sections. Results show that stationary methodologies are outperformed by all LVA techniques, even when the LVA framework has minimal guidance on parameterization. Findings also show that additional improvements are gained by considering parameter inference where the LVA orientations and point data are used to infer the local range of anisotropy. Considering LVA for geostatistical modeling of the deposit considered in this work results in better reproduction of curvilinear geological features.

     

Typing Mineral Deposits Using Their Grades and Tonnages in an Artificial Neural Network

A test of the ability of a probabilistic neural network to classify deposits into types on the basis of deposit tonnage and average Cu, Mo, Ag, Au, Zn, and Pb grades is conducted. The purpose is to examine whether this type of system might serve as a basis for integrating geoscience information available in large mineral databases to classify sites by deposit type. Benefits of proper classification of many sites in large regions are relatively rapid identification of terranes permissive for deposit types and recognition of specific sites perhaps worthy of exploring further.Total tonnages and average grades of 1,137 well-explored deposits identified in published grade and tonnage models representing 13 deposit types were used to train and test the network. Tonnages were transformed by logarithms and grades by square roots to reduce effects of skewness. All values were scaled by subtracting the variable's mean and dividing by its standard deviation. Half of the deposits were selected randomly to be used in training the probabilistic neural network and the other half were used for independent testing. Tests were performed with a probabilistic neural network employing a Gaussian kernel and separate sigma weights for each class (type) and each variable (grade or tonnage).Deposit types were selected to challenge the neural network. For many types, tonnages or average grades are significantly different from other types, but individual deposits may plot in the grade and tonnage space of more than one type. Porphyry Cu, porphyry Cu-Au, and porphyry Cu-Mo types have similar tonnages and relatively small differences in grades. Redbed Cu deposits typically have tonnages that could be confused with porphyry Cu deposits, also contain Cu and, in some situations, Ag. Cyprus and kuroko massive sulfide types have about the same tonnages. Cu, Zn, Ag, and Au grades. Polymetallic vein, sedimentary exhalative Zn-Pb, and Zn-Pb skarn types contain many of the same metals. Sediment-hosted Au, Comstock Au-Ag, and low-sulfide Au-quartz vein types are principally Au deposits with differing amounts of Ag.Given the intent to test the neural network under the most difficult conditions, an overall 75% agreement between the experts and the neural network is considered excellent. Among the largestclassification errors are skarn Zn-Pb and Cyprus massive sulfide deposits classed by the neuralnetwork as kuroko massive sulfides—24 and 63% error respectively. Other large errors are the classification of 92% of porphyry Cu-Mo as porphyry Cu deposits. Most of the larger classification errors involve 25 or fewer training deposits, suggesting that some errors might be the result of small sample size. About 91% of the gold deposit types were classed properly and 98% of porphyry Cu deposits were classes as some type of porphyry Cu deposit. An experienced economic geologist would not make many of the classification errors that were made by the neural network because the geologic settings of deposits would be used to reduce errors. In a separate test, the probabilistic neural network correctly classed 93% of 336 deposits in eight deposit types when trained with presence or absence of 58 minerals and six generalized rock types. The overall success rate of the probabilistic neural network when trained on tonnage and average grades would probably be more than 90% with additional information on the presence of a few rock types.     

Support and Information Effect Modeling for Recoverable Reserve Estimation of a Beach Sand Deposit in India

Application of geostatistics in estimating recoverable reserves of beach sand deposit is rare. This paper made an attempt to estimate local recoverable reserves using disjunctive kriging and discrete Gaussian model considering support and information effects for a beach sand deposit located in the eastern part of India. The dependence of different selective mining unit (SMU) sizes and different production sampling strategies on the estimated tonnage, metal quantity, and the ore tonnage versus metal quantity relationships has been examined. The results of the study show that nonlinear geostatistics should be used for more precise assessment of the grade, ore tonnage, and metal quantity and their relationships, which are necessary for recoverable reserve estimation. In selective mining operation, both support and information effects have significant influence on recoverable reserve. Recoverable reserve estimation based on SMU involves estimating grade distributions of mining unit with much bigger support than the available drill core sample data. Information effect comes into picture from the real scenario where the actual grades of the blocks remain unknown even during mining. At the mining stage, discrimination of ore and waste blocks is carried out based on estimated grades of the production samples and it is likely that the blocks might be misclassified as either ore or waste and thus sent to wrong destination. Information effect modeling makes the estimation more reliable by taking care of misclassification.     

Probabilistic Estimates of Number of Undiscovered Deposits and Their Total Tonnages in Permissive Tracts Using Deposit Densities

Empirical evidence indicates that processes affecting number and quantity of resources in geologic settings are very general across deposit types. Sizes of permissive tracts that geologically could contain the deposits are excellent predictors of numbers of deposits. In addition, total ore tonnage of mineral deposits of a particular type in a tract is proportional to the type’s median tonnage in a tract. Regressions using size of permissive tracts and median tonnage allow estimation of number of deposits and of total tonnage of mineralization. These powerful estimators, based on 10 different deposit types from 109 permissive worldwide control tracts, generalize across deposit types. Estimates of number of deposits and of total tonnage of mineral deposits are made by regressing permissive area, and mean (in logs) tons in deposits of the type, against number of deposits and total tonnage of deposits in the tract for the 50th percentile estimates. The regression equations (R ² = 0.91 and 0.95) can be used for all deposit types just by inserting logarithmic values of permissive area in square kilometers, and mean tons in deposits in millions of metric tons. The regression equations provide estimates at the 50th percentile, and other equations are provided for 90% confidence limits for lower estimates and 10% confidence limits for upper estimates of number of deposits and total tonnage. Equations for these percentile estimates along with expected value estimates are presented here along with comparisons with independent expert estimates. Also provided are the equations for correcting for the known well-explored deposits in a tract. These deposit-density models require internally consistent grade and tonnage models and delineations for arriving at unbiased estimates.     

Sequential Simulation Approach to Modeling of Multi-seam Coal Deposits with an Application to the Assessment of a Louisiana Lignite

There are multiple ways to characterize uncertainty in the assessment of coal resources, but not all of them are equally satisfactory. Increasingly, the tendency is toward borrowing from the statistical tools developed in the last 50 years for the quantitative assessment of other mineral commodities. Here, we briefly review the most recent of such methods and formulate a procedure for the systematic assessment of multi-seam coal deposits taking into account several geological factors, such as fluctuations in thickness, erosion, oxidation, and bed boundaries. A lignite deposit explored in three stages is used for validating models based on comparing a first set of drill holes against data from infill and development drilling. Results were fully consistent with reality, providing a variety of maps, histograms, and scatterplots characterizing the deposit and associated uncertainty in the assessments. The geostatistical approach was particularly informative in providing a probability distribution modeling deposit wide uncertainty about total resources and a cumulative distribution of coal tonnage as a function of local uncertainty.     

A Chart for Judging Optimal Sample Spacing for Ore Grade Estimation

Mineral deposit grades are usually estimated using data from samples of rock cores extracted from drill holes. Commonly, mineral deposit grade estimates are required for each block to be mined. Every estimated grade has always a corresponding error when compared against real grades of blocks. The error depends on various factors, among which the most important is the number of correlated samples used for estimation. Samples may be collected on a regular sampling grid and, as the spacing between samples decreases, the error of grade estimated from the data generally decreases. Sampling can be expensive. The maximum distance between samples that provides an acceptable error of grade estimate is useful for deciding how many samples are adequate. The error also depends on the geometry of a block, as lower errors would be expected when estimating the grade of large-volume blocks, and on the variability of the data within the region of the blocks. Local variability is measured in this study using the coefficient of variation (CV). We show charts analyzing error in block grade estimates as a function of sampling grid (obtained by geostatistical simulation), for various block dimensions (volumes) and for a given CV interval. These charts show results for two different attributes (Au and Ni) of two different deposits. The results show that similar errors were found for the two deposits, although they share similar features: sampling grid, block volume, CV, and continuity model. Consequently, the error for other attributes with similar features could be obtained from a single chart.     

New Tabu Algorithm for Positioning Mining Drillholes with Blocks Uncertainty

This paper proposes a new approach to the mining exploration drillholes positioning problem (DPP) that incorporates both geostatistical and optimization techniques. A metaheuristic was developed to solve the DPP taking into account an uncertainty index that quantifies the reliability of the current interpretation of the mineral deposit. The uncertainty index was calculated from multiple deposit realizations obtained by truncated Gaussian simulations conditional to the available drillholes samplings. A linear programming model was defined to select the subset of future drillholes that maximizes coverage of the uncertainty. A Tabu Search algorithm was developed to solve large instances of this set partitioning problem. The proposed Tabu Search algorithm is shown to provide good quality solutions approaching 95% of the optimal solution in a reasonable computing time, allowing close to optimal coverage of uncertainty for a fixed investment in drilling.

     

Public attitudes and policies toward mineral resources on the brink of the 21st century

In this issue, we feature an article by W. David Menzie, a research geologist with the U.S. Geological Survey, Reston, Virginia. Dr. Menzie is a leading expert on quantitative mineral-resource assessment. He has made significant contributions to quantitative assessment methodologies through the development of spatial mineral deposit density models, grade and tonnage models, and the design of metrics for describing mineral deposit occurrences. He has also studied the geology and mineral resources of the Circle quadrangle, Alaska. Dr. Menzie earned a B.S. degree in geology from Dickinson College, an M.S. in geology, an M.A. in statistics, and a Ph.D. in Geology from the Pennsylvania State University.     

An Objective Approach to the Definition of Some Mining Terms: Application to the Hasançelebi Iron Ore Deposit

In the mining industry, definitions of terms such as geologic resource, geologic reserve, mineral resource, mineral, mineral reserve, ore, ore reserve, mineable reserve, and industrial minerals always have been debated, and have caused much confusion. The process of arriving at these definitions requires complicated exploration, calculation, and evaluation. Based on such work, the definitions about the mineral property will be distinctly different. The aims of this paper is to discuss and compare these definitions, and then contrast the differences among them in the example of the Hasançelebi iron are deposit, which is an important source of iron ore in Turkey.     

Geostatistical Algorithm Selection for Mineral Resources Assessment and its Impact on Open-pit Production Planning Considering Metal Grade Boundary Effect

Mine planning is influenced by many sources of uncertainty. Significant sources of geological uncertainty in mine planning include uncertainty in layout of geological domains and uncertainty in metal grades. These two sources of uncertainty cannot be modeled separately because the distribution of the grade is controlled usually by geological domains. Two approaches exist for combining these two sources of uncertainty: the joint simulation approach and the cascade approach. In this paper, these two approaches were compared using a real case study. To this end, uncertainty in iron grade (quantitative variable) and ore zones (qualitative variable) was modeled using both approaches. There were some considerable differences in the results obtained by each approach, which confirm the importance of choosing the most appropriate approach with consideration of the dominate features of a deposit.

     

Sparse Data Division Using Data Segmentation and Kohonen Network for Neural Network and Geostatistical Ore Grade Modeling in Nome Offshore Placer Deposit

In this paper, sparse data problem in neural network and geostatistical modeling for ore-grade estimation was addressed in the Nome offshore placer gold deposit. The problem of sparse data arises because of the random data division into training, validation, and test subsets during ore-grade modeling. In this regard, the possibility of generating statistically dissimilar data subsets by random data division was also explored through a simulation exercise. A combined approach of data segmentation and application of a Kohonen network then was used to solve the data division problem. Two neural networks and five kriging models were applied for grade modeling. The neural network was trained using an early stopping method. Performance evaluation of the models was carried out on the test data set. The study results indicated that all the models that were investigated in this study performed almost equally. It was also revealed that by using the secondary variable watertable depth the neural network and the kriging models slightly improved their prediction precision. Further, the overall R ² of the models was poor as a result of high nugget (noisy) component in ore-grade variation.     

Selection of Optimal Thresholds for Estimation and Simulation Based on Indicator Values of Highly Skewed Distributions of Ore Data

This study strives to outline a geostatistics model for estimation and simulation of the Qolqoleh gold ore deposit located in Saqqez, NW of Iran. Considering that this gold deposit contains high-grade values, accurate evaluation of such values is of high importance, and therefore different methods based on indicator values, such as full indicator kriging (FIK) and sequential indicator simulation (SIS), have been employed to improve the accuracy of estimation and simulation of high-grade values. FIK and SIS cover the full range of grades based on several thresholds on the indicator data. The cumulative distribution function (CDF) is typically used for selection of threshold values. Given the highly skewed distribution of gold grade and its intense fluctuations, the number of thresholds is increased using CDF, which in turn results in a whole lot of calculations. To reduce the volume of calculations, the number–size (N–S) fractal model has been used to select thresholds. From such a model, all optimal thresholds are chosen with respect to geology and the unnecessary thresholds are excluded from selection. Thus, a study of the selection of optimal thresholds for estimation and simulation of a gold ore resource by means of FIK and SIS, respectively, based on thresholds selected using the N–S fractal model is presented. Finally, it is proved that results of these geostatistical methods based on thresholds selection from the N–S model appear to be better-positioned to explain ore grade variability compared to thresholds selected from the CDF and threshold selection from the N–S model is more effective for reducing the volume of required calculations.     

Statistical sampling of the distribution of uranium deposits using geologic/geographic clusters

The concept of geologic/geographic clusters was developed particularly to study grade and tonnage models for sandstone-type uranium deposits. A cluster is a grouping of mined as well as unmined uranium occurrences within an arbitrary area about 8 km across. A cluster is a statistical sample that will reflect accurately the distribution of uranium in large regions relative to various geologic and geographic features. The example of the Colorado Plateau Uranium Province reveals that only 3 percent of the total number of clusters is in the largest tonnage-size category, greater than 10,000 short tons U₃O₈, and that 80 percent of the clusters are hosted by Triassic and Jurassic rocks. The distributions of grade and tonnage for clusters in the Powder River Basin show a wide variation; the grade distribution is highly variable, reflecting a difference between roll-front deposits and concretionary deposits, and the Basin contains about half the number in the greater-than-10,000 tonnage-size class as does the Colorado Plateau, even though it is much smaller. The grade and tonnage models should prove useful in finding the richest and largest uranium deposits.     

Geostatistical Modeling of Ore Grade Distribution from Geomorphic Characterization in a Laterite Nickel Deposit

Due to growing consumption of nickel (Ni) in a range of industries, the demand for Ni has increased rapidly around the world. This trend requires a more precise estimation of available Ni grade deposits and an identification of factors controlling the grade distribution. To achieve these requirements, this study applies geostatistical techniques to spatial modeling of the Ni grade in a laterite Ni deposit, with reference to geomorphic features such as slope gradient and the thickness of limonite and saprolite zones. The Sorowako area in Sulawesi Island, Indonesia, was chosen as a case study area because it has a representative laterite Ni deposit with large reserves. Chemical content data from drillhole cores at 294 points were used for the analysis. The slope gradient was found to have a remarkable correlation with the thickness of the limonite zone, but there was no correlation between the thickness of the limonite and the saprolite zones above the bedrock. One important feature was a general correlation between the thickness of the saprolite zone and the maximum Ni grade in this zone: the grade increases with the thickness of the zone. Co-kriging was adopted to incorporate this correlation into estimating the maximum Ni grade in the saprolite zone. As a result, the maximum Ni grade in the saprolite zone tends to be high mainly in areas of slight slope. The Ni accumulation at this topographic feature probably originates from deep weathering by groundwater infiltrating through well-developed rock fractures.