No abstract

The definition of the search neighbourhood in kriging can have a significant impact on the resulting estimates. Stationary domains are usually estimated using a unique search strategy for the entire domain. However, the use of a global search neighbourhood ignores the local variations within each domain, i.e. all blocks are interpolated using a unique search strategy. In this paper, localised kriging parameter optimisation (LKPO) is proposed as an alternative methodology that considers the best 'local estimation parameter settings' block by block. The optimisation process is based on absolute error minimisation obtained in crossvalidation. Two datasets are presented, the first is a synthetic mineral deposit (2D) and the second is a gold deposit (3D). A wide variety of validation checks show that the use of local kriging parameters significantly improves the grade estimation, obtaining more precise and accurate results than the methodologies currently available in the geostatistical literature.

No abstract

Estimation of some mineral deposits involves chemical species or a granulometric mass balance that constitute a closed constant sum (e.g., 100%). Data that add up to a constant are known as compositional data (CODA). Classical geostatistical estimation methods (e.g., kriging) are not satisfactory when CODA are used, since bias is expected when estimated mean block values are back-transformed to the original space. CODA methods use nonlinear transformations, and when the transformed data are interpolated, they cannot be returned directly to the space of the original data. If these averages are back-transformed using the inverse function, bias is generated. To avoid this bias, this article proposes geostatistical simulation of the isometric logratio ratio (ilr) transformations back-transforming point simulated values (instead of block estimations), with the averaging being postponed to the end of the process. The results show that, in addition to maintaining the mass balance and the correlations among the variables, the means (E-types) of the simulations satisfactorily reproduce the statistical characteristics of the grades without any sort of bias. A complete case study of a major bauxite deposit illustrates the methodology.

Geostatistical modelling of grades in mineral deposits often requires the simulation of multiple related variables that have sum and fraction constraints. Sum constraints occur when the sum of some variables may not exceed or must equal a given constant. Fraction constraints occur when a variable may not exceed another variable. In this case, the geostatistical simulations should reproduce the histograms, variograms, multivariate relationships and honour the constraints. We present a methodology for the geostatistical simulation of multiple related variables that considers sum and fraction constraints. The methodology is illustrated in a case study. The data were obtained from a bauxite deposit. First, the original variables were re-expressed as ratios. Second, the variables were transformed using the Projection Pursuit Multivariate Transform (PPMT). Then the PPMT transformed variables were simulated independently using sequential Gaussian simulation and back-transformed. The simulations reproduced the histograms, variograms and bivariate relationships and honoured the sum and fraction constraints.