Approximate local confidence intervals under change of support

Approximate local confidence intervals are constructed from uncertainty models in the form of the conditional distribution of the random variable Z given values of variables [Z_i, i=1,...,n]. When the support of the variable Z is any support other than that of the data, the conditional distributions require a change of support correction. This paper investigates the effect of change of support on the approximate local confidence intervals constructed by cumulative indicator kriging, class indicator kriging, and probability kriging under a variety of conditions. The conditions are generated by three simulated deposits with grade distributions of successively higher degree of skewness; a point support and two different block supports are considered. The paper also compares the confidence intervals obtained from these methods using the most used measures of confidence interval effectiveness.     

The Bi-Gaussian approach: A simple method for recovery estimation

A new and simple method is proposed to obtain estimates of recovery functions: the Bi-Gaussian approach. Existing methods estimate recovery functions with conditional distributions where the conditioning set is all the data available. Here instead the simple kriging estimate of the Gaussian transform is proposed to be used. Results in the point recovery case are identical to the multi-Gaussian approach of Verly (1983, 1984), whereas in the non-point-support situation, an approximation is derived which saves computer time as compared to employing the strict multi-Gaussian hypothesis. Two examples compare favorably with the well-established disjunctive kriging method (discrete Gaussian model).     

Estimation of linear combinations and co-kriging

Utilizing the matrix formulation of co-kriging developed previously by the author, the relationship between direct kriging of linear combinations and linear combinations of co-kriged variables is developed. Conditions for equality of the estimators and the kriging variances are examined. By presenting the problem in the context of Hilbert spaces the general relationship is clarified.     

Comparison between different kriging estimators

Six different geostatistical estimators (linear kriging, lognormal kriging, and disjunctive kriging, each with and without a nonbias, i.e., universality condition) were compared using data from a polymetallic deposit in Algeria. The differences between estimators with and without the nonbias condition were far more pronounced than between the different kriging methods. This highlights the importance of choosing an appropriate stationarity model for the data. The criterion concerning kriging weight of the mean in simple kriging, proposed by Remacre (1984, 1987) and Rivoirard (1984) was found to be helpful for determining blocks where the choice of the stationarity hypothesis was critical.     

Models for Support and Information Effects: A Comparative Study

The recoverable reserves in an ore deposit depend on several factors, in particular the size of the selective mining units (support effect) and the misclassifications when sending these units to mill or dump according to their estimated grade (information effect). Both effects imply a loss of selectivity and have to be correctly forecasted. In this work, several models are reviewed and applied to a synthetic ore deposit characterized by a highly skewed grade histogram and a spatial connectivity of high grades. The affine correction, mosaic correction, and discrete Gaussian model are compared when assessing the global recoverable reserves, whereas local estimations are performed by indicator kriging with affine correction, bigaussian disjunctive kriging, and multigaussian conditional expectation. Despite their convenience and simplicity, distribution-free methods like affine correction or indicator kriging have a poorer accuracy than the other methods. In the global framework, the discrete Gaussian model is a better alternative and is based on mild assumptions. Local estimations are not accurate and may be improved by resorting to a more suitable parametric model or to conditional simulations.     

An Empirical Comparison of Kriging Methods for Nonlinear Spatial Point Prediction

Spatial prediction is a problem common to many disciplines. A simple application is the mapping of an attribute recorded at a set of points. Frequently a nonlinear functional of the observed variable is of interest, and this calls for nonlinear approaches to prediction. Nonlinear kriging methods, developed in recent years, endeavour to do so and additionally provide estimates of the distribution of the target quantity conditional on the observations. There are few empirical studies that validate the various forms of nonlinear kriging. This study compares linear and nonlinear kriging methods with respect to precision and their success in modelling prediction uncertainty. The methods were applied to a data set giving measurements of the topsoil concentrations of cobalt and copper at more than 3000 locations in the Border Region of Scotland. The data stem from a survey undertaken to identify places where these trace elements are deficient for livestock. The comparison was carried out by dividing the data set into calibration and validation sets. No clear differences between the precision of ordinary, lognormal, disjunctive, indicator, and model-based kriging were found, neither for linear nor for nonlinear target quantities. Linear kriging, supplemented with the assumption of normally distributed prediction errors, failed to model the conditional distribution of the marginally skewed data, whereas the nonlinear methods modelled the conditional distributions almost equally well. In our study the plug-in methods did not fare any worse than model-based kriging, which takes parameter uncertainty into account.     

Compensating for estimation smoothing in kriging

Smoothing is a characteristic inherent to all minimum mean-square-error spatial estimators such as kriging. Cross-validation can be used to detect and model such smoothing. Inversion of the model produces a new estimator—compensated kriging. A numerical comparison based on an exhaustive permeability sampling of a 4-ft² slab of Berea Sandstone shows that the estimation surface generated by compensated kriging has properties intermediate between those generated by ordinary kriging and stochastic realizations resulting from simulated annealing and sequential Gaussian simulation. The frequency distribution is well reproduced by the compensated kriging surface, which also approximates the experimental semivariogram well—better than ordinary kriging, but not as well as stochastic realizations. Compensated kriging produces surfaces that are more accurate than stochastic realizations, but not as accurate as ordinary kriging.     

不同采样点数量下土壤有机质含量空间预测方法对比

运用普通克里格、泛克里格、协同克里格和回归克里格4种方法,结合由DEM获取的高程因子以及土壤全氮和阳离子交换量（CEC）,预测了黑龙江省海伦市耕地有机质含量的空间分布。不同样点数量下海伦市土壤有机质含量的空间变异结构分析表明,样点数量多并不一定能够识别土壤有机质含量的结构性连续组分,最优化的布置采样点位置可能比单纯增加...     

Reducing the Impact of Outliers in Ore Reserves Estimation

Mining applications commonly faces surprising high values designated as outliers. These values impact dramatically statistical analysis and interpretation. A comprehensive analysis on the causes for the presence of unexpected high values was recommended. However, if an erroneous value was accepted as a part of the solution, some form of correction is recommended. A methodology based on the robust kriging (RoK) algorithm is proposed to be used in exploratory data analysis and also to deal with problems associated with the presence of outliers in the sample data set. The efficiency of RoK method as an interpolator is tested in different types of mineralizations. Importantly, the parent population from which the data was sampled is available, thus allowing direct quantitative assessment of the effectiveness of the estimation technique. The performance of the method is tested in the context of ore reserves estimation. RoK model is compared to models generated by ordinary kriging, median indicator kriging, and lognormal kriging. RoK proved to be more accurate and more precise than those methods reducing substantially the number of misclassified blocks.     

Normal and lognormal estimation

A comprehensive theoretical study of the problem of estimation of regionalized variables with normal or lognormal distribution is presented. Unbiased linear estimators are derived, under both assumptions that the population mean is known and unknown, and their error variance is calculated. The minimum variance kriging estimators are studied in more detail and are compared with the conditional expectations. The emphasis is on the study of lognormally distributed variates. The derived mathematical formulas are applicable to the optimal contouring of sample values with the appropriate distribution, as well as the optimal estimation of blocks of ore in mineral deposits.     

Kriging Regionalized Positive Variables Revisited: Sample Space and Scale Considerations

Frequently, regionalized positive variables are treated by preliminarily applying a logarithm, and kriging estimates are back-transformed using classical formulae for the expectation of a lognormal random variable. This practice has several problems (lack of robustness, non-optimal confidence intervals, etc.), particularly when estimating block averages. Therefore, many practitioners take exponentials of the kriging estimates, although the final estimations are deemed as non-optimal. Another approach arises when the nature of the sample space and the scale of the data are considered. Since these concepts can be suitably captured by an Euclidean space structure, we may define an optimal kriging estimator for positive variables, with all properties analogous to those of linear geostatistical techniques, even for the estimation of block averages. In this particular case, no assumption on preservation of lognormality is needed. From a practical point of view, the proposed method coincides with the median estimator and offers theoretical ground to this extended practice. Thus, existing software and routines remain fully applicable.     

Determining the best search neighbourhood in reserve estimation, using geostatistical method: A case study anomaly No 12A iron deposit in central Iran

Ordinary kriging and non-linear geostatistical estimators are now well accepted methods in mining grade control and mine reserve estimation. In kriging, the search volume or ‘kriging neighbourhood’ is defined by the user. The definition of the search space can have a significant impact on the outcome of the kriging estimate. In particular, too restrictive neighbourhood, can result in serious conditional bias. Kriging is commonly described as a ‘minimum variance estimator’ but this is only true when the neighbourhood is properly selected. Arbitrary decisions about search space are highly risky. The criteria to consider when evaluating a particular kriging neighbourhood are the slope of the regression of the ‘true’ and ‘estimated’ block grades, the number of kriging negative weights and the kriging variance. Search radius is one of the most important parameters of search volume which often is determined on the basis of influence of the variogram. In this paper the above-mentioned parameters are used to determine optimal search radius.     

Comparative performance of indicator algorithms for modeling conditional probability distribution functions

This paper compares the performance of four algorithms (full indicator cokriging. adjacent cutoffs indicator cokriging, multiple indicator kriging, median indicator kriging) for modeling conditional cumulative distribution functions (ccdf).The latter three algorithms are approximations to the theoretically better full indicator cokriging in the sense that they disregard cross-covariances between some indicator variables or they consider that all covariances are proportional to the same function. Comparative performance is assessed using a reference soil data set that includes 2649 locations at which both topsoil copper and cobalt were measured. For all practical purposes, indicator cokriging does not perform better than the other simpler algorithms which involve less variogram modeling effort and smaller computational cost. Furthermore, the number of order relation deviations is found to be higher for cokriging algorithms, especially when constraints on the kriging weights are applied.     

Two Ordinary Kriging Approaches to Predicting Block Grade Distributions

Multigaussian kriging aims at estimating the local distributions of regionalized variables and functions of these variables (transfer or recovery functions) at unsampled locations. In this paper, we focus on the evaluation of the recoverable reserves in an ore deposit accounting for a change of support and information effect caused by ore/waste misclassifications. Two approaches are proposed: the multigaussian model with Monte Carlo integration and the discrete Gaussian model. The latter is simpler to use but requires stronger hypotheses than the former. In each model, ordinary multigaussian kriging gives unbiased estimates of the recoverable reserves that do not utilize the mean value of the normal score data. The concepts are illustrated through a case study on a copper deposit which shows that local estimates of the metal content based on ordinary multigaussian kriging are close to the optimal conditional expectation when the data are abundant and are not dominated by the global mean when the data are scarce. The two proposed approaches (Monte Carlo integration and discrete Gaussian model) lead to similar results when compared to two other geostatistical methods: service variables and ordinary indicator kriging, which show strong deviations from conditional expectation.     

中天山乌拉斯台地区矿化异常提取及评价

以多重分形理论为基础,对中天山乌拉斯台地区铜多金属元素的岩屑测量数据,采用C-A法获得铜多金属的异常下限值,将其作为阈值进行指示克里格插值,绘制研究区的铜多金属地球化学异常图。研究显示,基于该方法获得的Cu矿化异常高值区主要集中在华力西早期第三侵入次的石英闪长岩和花岗闪长岩岩体中,受北西向和次级北东向断裂构造控制明显,该异常区可以作为寻找热液型铜多金属矿产的重要远景区。该方法对于地球化学数据空间变异性强烈的地区,较之普通克里格插值法具有更好的地球化学异常识别能力和高值信息重建能力,所得结果的最高累计频率值范围与已知矿化点的空间位置吻合度更高,在地球化学异常信息提取工作中具有推广意义。     

Comparing the robustness of ordinary kriging and lognormal kriging: Outlier resistance

Ordinary kriging is well-known to be optimal when the data have a multivariate normal distribution (and if the variogram is known), whereas lognormal kriging presupposes the multivariate lognormality of the data. But in practice, real data never entirely satisfy these assumptions. In this article, the sensitivity of these two kriging estimators to departures from these assumptions and in particular, their resistance to outliers is considered. An outlier effect index designed to assess the effect of a single outlier on both estimators is proposed, which can be extended to other types of estimators. Although lognormal kriging is sensitive to slight variations in the sill of the variogram of the logs (i.e., their variance), it is not influenced by the estimate of the mean of the logs.This paper was presented at MGUS 87 Conference, Redwood City, California, 14 April 1987.     

On the equivalence of kriging and maximum entropy estimators

This study compares kriging and maximum entropy estimators for spatial estimation and monitoring network design. For second-order stationary random fields (a subset of Gaussian fields) the estimators and their associated interpolation error variances are identical. Simple lognormal kriging differs from the lognormal maximum entropy estimator, however, in both mathematical formulation and estimation error variances. Two numerical examples are described that compare the two estimators. Simple lognormal kriging yields systematically higher estimates and smoother interpolation surfaces compared to those produced by the lognormal maximum entropy estimator. The second empirical comparison applies kriging and entropy-based models to the problem of optimizing groundwater monitoring network design, using six alternative objective functions. The maximum entropy-based sampling design approach is shown to be the more computationally efficient of the two.     

Indicator principal component kriging

An alternative to multiple indicator kriging is proposed which approximates the full coindicator kriging system by kriging the principal components of the original indicator variables. This transformation is studied in detail for the biGaussian model. It is shown that the cross-correlations between principal components are either insignificant or exactly zero. This result allows derivation of the conditional cumulative density function (cdf) by kriging principal components and then applying a linear back transform. A performance comparison based on a real data set (Walker Lake) is presented which suggests that the proposed method achieves approximation of the conditional cdf equivalent to indicator cokriging but with substantially less variogram modeling effort and at smaller computational cost.     

Robustness of variograms and conditioning of kriging matrices

Current ideas of robustness in geostatistics concentrate upon estimation of the experimental variogram. However, predictive algorithms can be very sensitive to small perturbations in data or in the variogram model as well. To quantify this notion of robustness, nearness of variogram models is defined. Closeness of two variogram models is reflected in the sensitivity of their corresponding kriging estimators. The condition number of kriging matrices is shown to play a central role. Various examples are given. The ideas are used to analyze more complex universal kriging systems.Research performed while on leave at Centre de Geóstatistique et de Morphologie Mathématique, Fontainebleau.     

Mapping by simple indicator kriging

The problem of predicting the type of an unsampled point in a two-type map is addressed using a procedure called simple indicator kriging. This procedure estimates the conditional probability that the point is of one type or the other given the types of sample points. Simple indicator kriging is applied to examples of certain map models. It is shown to perform well, and in some cases to be exact in a specified sense. Use of simple indicator kriging to estimate multivariate normal quandrant probabilities is explored.