Recoverable reserves and support effects when optimizing open pit mine designs

Orebody modelling, support effects and the estimation of recoverable reserves are key parts of open pit optimization studies. A case study is presented on the estimation of recoverable reserves using an implementation of indicator kriging where metal quantity is used to select cutoffs, and support corrections founded on a conditional simulation approach. Mining selectivity is explored in the subsequent optimization study to compare results from indicator kriging of grade estimates on a regular size blocks and indicator kriging estimates on small size blocks. The use of indicator kriging models adjusted for a given selectivity and the use of grade proportions in each block for the optimization study, provide a presentation of the expected ore recovery for a predefined level of selectivity. The case study shows that indicator kriging estimation with full accounting of block grade distributions generates substantially better results in the pit optimization study. In addition, the adverse effects of small blocks and over-smoothing on optimization results are illustrated.

Conditional simulation of random fields by successive residuals

This paper presents a new approach to the LU decomposition method for the simulation of stationary and ergodic random fields. The approach overcomes the size limitations of LU and is suitable for any size simulation. The proposed approach can facilitate fast updating of generated realizations with new data, when appropriate, without repeating the full simulation process. Based on a novel column partitioning of the L matrix, expressed in terms of successive conditional covariance matrices, the approach presented here demonstrates that LU simulation is equivalent to the successive solution of kriging residual estimates plus random terms. Consequently, it can be used for the LU decomposition of matrices of any size. The simulation approach is termed conditional simulation by successive residuals as at each step, a small set (group) of random variables is simulated with a LU decomposition of a matrix of updated conditional covariance of residuals. The simulated group is then used to estimate residuals without the need to solve large systems of equations.

Successive nonparametric estimation of conditional distributions

Spatial characterization of non-Gaussian attributes in earth sciences and engineering commonly requires the estimation of their conditional distribution. The indicator and probability kriging approaches of current nonparametric geostatistics provide approximations for estimating conditional distributions. They do not, however, provide results similar to those in the cumbersome implementation of simultaneous cokriging of indicators. This paper presents a new formulation termed successive cokriging of indicators that avoids the classic simultaneous solution and related computational problems, while obtaining equivalent results to the impractical simultaneous solution of cokriging of indicators. A successive minimization of the estimation variance of probability estimates is performed, as additional data are successively included into the estimation process. In addition, the approach leads to an efficient nonparametric simulation algorithm for non-Gaussian random functions based on residual probabilities.