Search Documents

This paper describes a detailed geostatistical study of producing copper orebodies at Mount Isa Mines Limited in Queensland, Australia. As a result of this work, kriging has been adopted for grade estimation for final stope design at Isa Mine. Extensive variogram and back analysis studies indicate that one relative variogram model can be used to give optimal kriged grade estimates from production sampling throughout the 100 million tonne copper orebody complex. The kriging method replaces a computer technique used at [sa Mine based on inverse distance weighting using parameters chosen from experience.Comparative stope designs based on 10 million tonnes of ore indicate a 16% increase in minable ore tonnage at a grade well above the cut-off, as a result of the geostatistical work.Although this paper is written as a practical description with simple explanations of geostatistical concepts, the theory and use of relative variograms has been expanded. A possible link is suggested between relative variograms and first order kinetics of natural mineral concentration processes.

Jan 1, 1982

By P Chapman, B Moosapoor

The objective of this paper is to present the potential applications of geostatistical analyses in planning CPTu investigations, as well as interpolation of CPTu data. The paper presents a case study of a geostatistical analysis of CPTu data undertaken between 2004 and 2017 for an upstreamraised TSF located in the Goldfields region of Western Australia. An initial model is developed using the 2004, 2009 and 2017 CPTu data. The model is used to predict in situ cone tip resistance (qc) measurements, as the interested primary variable, at the selected locations where 2019 CPTu were undertaken, and the results of the prediction and the actual data are compared.

Jul 1, 2021

By Leach B. G

A geostatistical study was made of diamond drill hole intersections within 18CC stope block at the C.S.A. Mine, Cobar, N.S.W. The stope block is located in 18CC orebody which is a massive sulphide - multiple vein, Copper-Zinc type, lacking clear definition on gangue mineralogy, base-metal assays and total sulphide considerations. Of the variograms constructed for inter- section width and copper, lead and zinc grades and accumulations, only zinc displayed behaviour amenable to model fitting and kriging procedures. Zinc grades were calcu- lated for sub-blocks comprising the total stope block and comparison with equivalent grades determined by traditional methods showed them to be generally less extreme in the highest and lowest grade sub-blocks. The total-block grades derived by the two methods were almost the same. It is suggested that the failure of copper to show variogram continuity might reflect the difficulty in defining the actual 18CC orebody, a proposition which is supported by the poor continuity in the variogram of horizontal width. In the case of lead the poor continuity may be due to the frequent manner of occurrence of galena as coarse grained patches of limited extent. The results of the study, did not provide a convincing reason to depart from the traditi- onal methods of grade estimation at the C.S.A. Mine.

Jan 1, 1979

By Royle A. G

The paper describes a geostatistical structural analysis of an 80 square mile lease of the Athabasca tar sands. Each step of the analysis is presented with the ultimate objective of condensing the principal structural characteristics of the chosen regionalised variables into the opera- tional form of the variogram model. The variogram model is then used to form kriged local and global estimates, with the corresponding estimation errors.

Jan 1, 1977

Knowledge of true block grades is obtained from a mathematical model of an orebody, simulated with sample grades normally distributed and spatially correlated according io a spherical semi-variogram model.The superiority of kriging to the polygonal method of estimation is demonstrated on this modelled deposit by comparing the true values with the kriged and polygonal estimates obtained from sampling the deposit. Finally, selection of blocks of the deposit on the basis of the estimated values, shows that if kriging is used, what is predicted is largely mined, whilst if the polygonal method is used large biases can occur depending on the position of the cut-off grade.

Jan 1, 1978

By Armstrong M

The highest grades of the Porgera gold deposit, Papua New Guinea, are contained in Zone 7 due to the presence of the Roamane fault, which acted as channelway of the gold-bearing fluids. Prior to the development and underground mining of this fault-controlled high-grade zone it was critical to provide grade estimates of relatively small mining units, and of the proportion of mining units above cutoff. To be applicable to Zone 7 of Porgera a grade estimation method should account for the following:1. the highly skewed and apparently lognormal distribution of the grades,2. the rapid drop off of the local mean grade and the variance as the distance to the Roamane fault increases,3. the presence of high grade pockets in the footwall, away from that fault, and4. the existence of variograms for the logs with sills parallel to the fault and without sills perpendicular to it.

Jan 1, 1990

In the estimation of reserves the known values of the variables of interest from surrounding drill holes are used to compute or estimate the unknown value of the orebody over a given region or at a given point of interest. The classical means of doing this are: 1. The Polygonal Method û where the region is divided into polygons centred on bore holes, each variable is then assumed to have constant value over each polygon. 2. Inverse distance methods, where the variable is assumed to vary linearly with the distance between the known values of the drill holes. 3. Inverse distance squared method û a quadratic variation is assumed for the variables between the bore holes.

Jan 1, 1983

This study begins with the semivario- gram modelling of variables rutile grade and clay content from a beach sand rutile deposit. This is followed by two dimensional block kriging. The deposit was divided into blocks of 30 x 30 x 1.5 m dimensions and the kriged estimate together with the associated kriging variances for each variable within each block were determined. Finally, using the results of the krig- ing program and considering the acceptable limit of clay content together with the cutoff grade for rutile, a method for optimization of mineable thickness and hence the tonnage of ore was applied and the results were dis- Cussed.

Jan 1, 1986

By W Patton, I Minniakhmetov, H Talebi, U Mueller

The sustainability of mining projects is linked to informed investment decisions based on public reporting of exploration and mineral resource estimation results. In Australia, public reporting guidelines are established by the Joint Ore Reserves Committee reporting code through the JORC Code (2012). Although the assessment of uncertainty in the results reported is a requirement, this is often communicated qualitatively and evaluated subjectively. This can become a liability particularly in the early stages of mining projects when spatial domains of geological interpretation and mineralisation envelopes inform resource estimations’ reliability. A recent review of JORC reports found no reports used quantitative assessments of the geological interpretation or mineral resource estimation envelopes, and 27 per cent of reports did not address the question of quality of geological interpretation. Reports that did address the quality used 19 different terms to communicate quality. This work presents methodologies for quantitative uncertainty assessment and communication and explores how they could be applied in public reporting practice. Quantitative uncertainty reporting would benefit mine planning and financial modelling as they normally assume the geology is 100 per cent right (since no error is reported or only a single deterministic model is provided). The complexity, cost, and additional work of doing a quantitative assessment could hinder a straightforward implementation. This could be overcome if mining companies budget for quantitative uncertainty assessment and associated professional development. A compulsory requirement for the inclusion of uncertainty assessments in public reporting or adopting standardised subjective language would improve industry practice.

Mar 22, 2022

By C V. Deutsch

"Geometallurgical variables often cause problems to conventional geostatistical workflows. There are many variables; some are compositional and some are non-additive. They often show: complex multivariate relationships undersampling relative to conventional grade variables unequal sampling; that is, all variables are not available at all data locations.A variety of approaches have been developed over the years, but geostatisticians remain mostly focused on modelling rock types and a few grade variables with an abundance of data. The problems are reviewed and a chained workflow is presented where appropriate techniques are employed in sequence to permit reliable spatial prediction of geometallurgical variables. The spatial models may be used to assess the value of additional data, for prediction of process performance, for improved mine planning and other optimisation tasks. The chained workflow will include data transformations such as ratios and normal scores that make the variables amenable to geostatistical modelling. Data imputation techniques to account for missing values are considered with special consideration to spatial correlation with the available data. Decorrelation techniques such as the principal components, min/max autocorrelation factors and projection pursuit multivariate transformation are employed to manage massively multivariate modelling. Full distribution modelling must be considered to avoid biases in the prediction of non-additive variables. The post-processing of geometallurgical models must be streamlined for robustness and efficiency. The precise workflow varies according to the specific purpose of the study. Guidelines and thoughts on best practice are presented.CITATION:Deutsch, C V, 2013. Geostatistical modelling of geometallurgical variables - problems and solutions, in Proceedings The Second AusIMM International Geometallurgy Conference (GeoMet) 2013, pp 7-16 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Sep 29, 2013

By P Landeros, D Benado, J Cornejo, A Pinochet, C Caviedes

In 2005, El Teniente mine began preconditioning the primary rock mass by hydraulic fracturing (HF). The major perceived benefits of this process are a decrease in the magnitude of the expected maximum seismic event and improved cave propagation. In addition, HF has resulted in new information on the propagation pressures of each fracture. Previous studies have attempted to relate propagation pressure to minor principal stress or confining stress.Considering the spatial continuity of propagation pressure and the significant quantity of collected data since the incorporation of the HF technique, minimum requirements for geostatistical modelling have been reached, permitting robust estimates inside the representative volume with acceptable errors (Cornejo, Muñoz and Rojas, 2012; Cornejo and Landeros, 2013). These results are then used to identify zones with higher magnitudes, which are mostly related to more complex geotechnical zones (intrusive and major faults) and caving singularities, making it possible to estimate the spatial extension of these zones.Finally, these results are integrated with other geomechanics aspects to support risk analysis of the mining plan.CITATION:Cornejo, J, Landeros, P, Benado, D, Caviedes, C and Pinochet, A, 2016. Geostatistical modelling of hydraulic fracturing pressures at El Teniente Mine, in Proceedings Seventh International Conference and Exhibition on Mass Mining (MassMin 2016), pp 183–190 (The Australasian Institute of Mining and Metallurgy: Melbourne).

May 9, 2016

By Raymond G. F

Geostatistics - a tool used by geologists and mine planning engineers to estimate grades and rapidly evaluate alternative mining strategies has reached an advanced stage of development for the silver-lead-zinc orebodies at Isa and Hilton Mines. Geostatistics has been used for over a decade in the massive copper orebodies but the stratiform lead-zinc orebodies are much more complex and provided a completely different set of challenges. Recent increases in computer processing speed and power provided by the Integrated Mine Planning System (IMPS) has aided the development of the geostatistics package. Also the formation of a project team and the subsequent in-house development which has drawn together the expertise of geologists, engineers and system analysts has contributed significantly to the success of the project.The most important feature of any approach to modelling the lead-zinc orebodies is that the kriging search must consider estimate and sample points according to their positions relative to bedding surfaces. By transforming co-ordinates as perpendicular offsets from a plane representing the local strike and dip, it is possible to elongate the ellipsoidal influence to parallel bedding. By interactively modelling surfaces using CAD facilities and then specifying these surfaces as 'hard' or 'soft' boundaries, the geologist has complete control over the kriging process.

Jan 1, 1993

By Leahey T. A

The variability of gold mineralization in the Warrego Gold Pod is described by classical statistical techniques using sample distrib- utions and analysis of variance; and by the use of geostatistical variogram analysis to define the spatial variability and anisotropy of the mineralization. Comparison is made between declared mine production from the Nineteen Stope and a series of ore reserve estimates derived by manual, numerical and disjunctive kriging methods.

Jan 1, 1979

By Sarkar B. C, Sen A. K

India is well endowed with vast reserves of banded hematite ore, about 12 billion tonnes, estimated by some conventional methods. Various authors have assessed this huge reserve using some advanced geostatistical methods. A Banded Iron Formation (BIF) belt named the `Singhbhum-Keonjhar-Bonai' belt in India has been chosen as a case study. In this paper, the experiences of Chiria iron ore, in particular, have been described. This study will finally create a dynamic mineral inventory after standardising the geostatistical parameters which will help to cater for the needs of the stringent quality requirements of the Indian steel industries.

Jan 1, 1995

By Sirayanone S, Rashad MZ

An important problem occurring in the iron and steel industry of Egypt is the presence of manganese oxide and chlorine with different percentages in the Bahariya Iron Ore Deposit, Western Desert, Egypt. In this work we used a geostatistical method known as the Volume- Variance Relationship to predict grade fluctua- tions of Fe% , Mn% and Cl% in the ore after blasting a certain volume from the ore deposit, on the other hand, if we have a certain limit for each of these chemical components, we can select the bench blasting dimension to give these limits.

Jan 1, 1986

Types of structure found in 'Jest Australian beach sands are described. Their recognition by geostatistical methods follows nth a discussion of the effectiveness of sampling patterns. Examples of structures are given, and an account of their effect on valuation procedures.

Jan 1, 1977

By M Field

Kimberlite pipes, as opposed to dykes, sill and secondary deposits, are the primary target for diamond exploration companies because they have simple geometries and can contain large volumes of potentially diamond-bearing ore. Once discovered, important decisions have to be made regarding the sampling of these geological entities to establish the total volume of ore present, and to define the internal geology, which can be complex in nature and controls the distribution of diamond grade. It is essential that as part of this process the uncertainty in constructed geological models are assessed. In this paper the Orapa AK1 kimberlite pipe in Botswana is used as a case study to explore a potential methodology for defining this uncertainty. The uncertainty of the rock volumes of different types can be characterised by a geostatistical approach. Geostatistical simulations, based upon a plausible representation of the geological bodies, provide different possible images of the orebody that can be considered as different realities. By testing different borehole layouts and drilling the orebody virtually on the computer, many different, but equi-probable realisations of the pipes can be produced. It is then possible to estimate the rocksÆ volumes and compare them with the actual ones. A statistical analysis of the relationship between estimation error and the number of boreholes can then be used to optimise the next sampling campaign. The rock volumes are obtained from a two-step simulation process. Firstly, the pipe geometry is simulated by means of introducing variations around the geological block model, by an original method called centre point simulation (CPS), which has been developed for this purpose. Secondly, the internal geology is then simulated within the pipe. The methodology for simulating the internal geological zones has to be adapted to the level of information and to the geological structure itself. The upper portion at Orapa AK1 is filled with sedimentary crater facies. Here plurigaussian simulation techniques can be applied that are adapted to the simulation of sedimentary facies. The central zones of the pipe are typically subvertically-oriented massive volcaniclastic deposits and breccias and the lower regions comprise a complex root zone. CPS is appropriate for these zones.

Jan 1, 2006

By Shrivastava P

The iron orebodies of the Pilbara region are formed by supergene enrichment of banded iron formation and rarely do such processes produce orebodies with uniform grade distrib- utions. The orebodies being mined or developed have to be worked to grade/quality constraints imposed by contractual agreement rather than natural characteristics. Hence it is imperative to use a non-biased ore reserve calculation method. The geostatistical study of a typical iron orebody is presented. Emphasised are the geological characteristics of the orebody manifested as a variable grade distribution of iron and impurities; the determination of variogram model; subsequent uses of the model for computation of ore reserves and kriged lock estimates; computation of dispersion variance and conditional simulation of orebody.

Jan 1, 1977

By Dimitrakopoulos R

Ore textures in base metal deposits are an important factor in the liberation of the economic components of the ore, metal recovery and reagent consumption during beneficiation. Optimising the mine profit requires knowing the relationship between ore texture and mineral processing behaviour, and the ability to predict the metal grade and the ore texture at a suitable scale. A new method that generates distributions of mesotextures for any mining block size or shape is proposed and includes: I . establishing a protocol for acquiring mesotexture data: 2. measuring the mesotexture spatial correlations using a multi-texture variogram; 3. generating multiple simulations of mesotextures; and 4. upscaling the simulations to generate the expected mining block/stope mesotexture distributions. Models of the ore textures built using the proposed methodology can be used to assess resources and reserves, enhance ore/waste discrimination and predict the metallurgical behaviour of parcels of ore.

Jan 1, 1997

By Ramani R. V

In this paper, a comparative performance evaluation of geostatistical, spectral, and fractal methods is made for short-scale variability prediction. The algorithms compared are: sequential Gaussian simulation method; a random fractal algorithm based on the generalised stochastic subdivision method; and a spectral simulation approach adopted for direct multi-dimensional simulation. The data for this study have been obtained from a mined out coal property. As this property was being mined, a 7315 m by 7315 m (24000 ft by 24000 ft) coal block was extensively sampled at an approximate spacing of 107 m (350 ft) and the 60-ton samples were analysed for sulphur in the 1.50 float gravity fraction. This data set is used here with the three methods to generate the simulated realisations, which are analysed for semivariogram reproduction and for simulation errors. Although all three methods converge to the true semivariogram in an average sense, the coefficient of variation of the average semivariogram for the Fourier integral method is minimum. The realisations generated by the sequential Gaussian method have minimum mean simulation error. The choice of a simulation method based on variability reproduction and simulation errors remains difficult.

Jan 1, 1995