Available in Isatis: a new efficient way of classifying mineral resources
During Geostats2016, Mines ParisTech Geostatistics Group presented an innovative way, easy and quick to use, of classifying mineral resources. A beta version of this innovation is available in Isatis 2016 on demand and will be integrated as a new tool in Isatis 2017.
In the mining industry, it is common to classify estimated resources as inferred, indicated or measured according to their confidence level which depends on the sample spacing and the expected production volumes.
This confidence level is classically derived either from the linear kriging variance or from a post-processing of conditional simulations. Both methods have their own limitations, the first one because it does not take into account data values (and the error distribution is rarely gaussian by nature), the second one because it is rather time-consuming and requires strong hypotheses. Isatis offers a third classification technique allowing the direct calculation of a Confidence Interval attached to the estimation of blocks in the framework of the Discrete Gaussian Model. Whilst producing CI’s faster than when resorting to conditional simulation, the mapping of confidence levels is less powerful than with simulations (where one can produce complex risk indices).
Whatever the method used, the confidence is computed for each individual block from the block model. The drawback is that resource categories depend on the size of the blocks and can vary quickly at small scale because of possible rapid variations
from a block to the next one (which is often referred to as “the spotted dog syndrome”).
A post-process is then required to smooth categories. Another drawback is that the block support is not appropriate for long-term management implying larger production volumes.
The Geostatistics Group from the French school Mines ParisTech proposes a new geostatistical classification of mineral resources which consists of two steps:
- Firstly, the spatial sampling density is characterized throughout the deposit using a specific volume quantity similar to the inverse of the classical density of sample points in space but taking into account the variographic structure of the grades of interest. The method requires the kriging at smu scale within “super” blocks with a size much larger than the actual block size and usually corresponding to specific production periods;
- Secondly, this volume is converted into a coefficient of variation corresponding to large production volumes, typically the expected annual volume to be mined. Then, a risk index is obtained by thresholding this coefficient of variation.
WHY USING THIS NEW INDEX?
- Robust classification independant from sample number, block size and domain extension
- Appropriate for mid and long-term resource management
- Coherent variation of the index in space
- Comparable index from one deposit to another
- Integrates uncertainties related to the sampling mesh and the spatial variability
- Objective resource classification criteria
- Accessible index, quick and easy to get, does not require post-processing
- Simplified communication between teams
- Facilitated comparison of different deposits
The examples below show and compare this risk index to classical classification on 2 cases:
Copper grade with irregular sampling pattern – Courtesy Mines ParisTech/Codelco:
Iron ore proportion with irregular sampling pattern – Courtesy Mines ParisTech/Vale:
Isatis 2016 offers a beta version of the tool Sampling Density Variance allowing the calculation of the spatial sampling density. It is part of Isatis Plugins which can be downloaded from Isatis download page.
The complete tool enabling the direct calculation of the Sampling Density Risk Index will be available in Isatis 2017.
- A Specific Volume to Measure the Spatial Sampling of Deposits – Rivoirard J. & Renard D. (Centre de Géosciences, Mines ParisTech) – in Mathematical Geosciences – October 2016, Volume 48, Issue 7, pp 791–809
- From the Spatial Sampling of a Deposit to Mineral Resources Classification – Rivoirard J. and Al – Presented at Geostats2016, Valencia