A parametric study of robustness of kriging variance as a function of range and relative nugget effect for a spherical semivariogram

In geostatistics, an estimation of blocks of a deposit is reported along with the variance of error made in their estimation. This calculation is based on the model chosen for the semivariogram of the deposit so that mistakes in its estimation can manifest themselves in the perception of accuracy with which blocks are known. Changes in kriging variance resulting from various amounts of error in modeling the relative nugget effect and range of the semivariogram are investigated for an extensive set of spherical semivariograms.     

重磁场的变差函数特征及应用

用变差函数研究重磁场的区域变化特征.变差函数的变程反映重磁场的相干范围, 块金效应反映随机干扰, 基台值反映变异程度.重磁场的理论模拟说明: 重力场的相干范围大于磁场, 重磁场变程主要取决于场源深度, 浅源重磁场变差函数近似为球状模型或指数模型, 深源重磁场近似为连续性更好的高斯模型.磁场场源深度近似等于变程的一半, 重力场场源深度近似等于变程的四分之一.湖北大冶铁矿垂直分量磁异常具有几何各向异性, 北西-南东走向, 变差函数推测磁铁矿平均深度为250m.磁异常小波多尺度分解细节和逼近部分磁场具有协调几何各向异性, 变差函数的各阶场源深度估计结果与功率谱估计结果吻合.      

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.     

Notes on the robustness of the kriging system

The robustness of the kriging system with respect to uncertainty of the theoretical variogram is investigated. Inequalities for possible changes of the kriging estimator and the estimation variance are derived. Results of a numerical study show that changes of kriging weights can be predicted partly with the help of the maximal kriging weight.     

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.     

含水层渗透性空间分布的指示克立格估值

详细介绍了指示克立格估值计算的理论和方法。以指示变异函数为基本工具分析了华北某地区第四系含水层渗透性空间分布的结构特征,结果表明该地区含水层渗透性存在明显的各向异性特征。水平方向上,X轴方向的相关性较Y轴方向的好,Z轴的相关性最差。用指示克立格法对未采样点处进行估值,估值结果显示含水层渗透性由山前向滨海逐渐变低,在垂直方向上,渗透性变化不明显,浅部比深部略好;同时给出了估计精度,并认为对估计精度不高的区域可通过增加适当的工程加以控制。最后用交叉验证法对估值结果进行了检验,证明建立的指示变异函数模型合理且估值效果较好。这一实际应用表明指示克立格法可以很好地描述第四系含水层渗透性的空间分布规律。     

Another look at anisotropy in geostatistics

A thorough geostatistical data analysis includes a careful study of how the data's second-order variation, as characterized by the semivariogram, depends on the relative orientation of data locations. If the semivariogram depends on only the (Euclidean) distance between locations, then the semivariogram is isotropic; otherwise, it is anisotropic. In this article, I take another look at the modeling of anisotropy in geostatistics. A new, more specific classification of types of anisotropy is proposed. More importantly, some heretofore inadequately understood implications of the dependence of various semivariogram attributes on direction are discussed, and the wisdom of some current practices for modeling the direction-dependence of these attributes is questioned.     

Spatial analyses of groundwater levels using universal kriging

For water levels, generally a non-stationary variable, the technique of universal kriging is applied in preference to ordinary kriging as the interpolation method. Each set of data in every sector can fit different empirical semivariogram models since they have different spatial structures. These models can be classified as circular, spherical, tetraspherical, pentaspherical, exponential, gaussian, rational quadratic, hole effect, K-bessel, J-bessel and stable. This study aims to determine which of these empirical semivariogram models will be best matched with the experimental models obtained from groundwater-table values collected from Mustafakemalpasa left bank irrigation scheme in 2002. The model having the least error was selected by comparing the observed water-table values with the values predicted by empirical semivariogram models. It was determined that the rational quadratic empirical semivariogram model is the best fitted model for the studied irrigation area.     

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.     

Distribution of kriging error and stationarity of the variogram in a coal property

If a particular distribution for kriging error may be assumed, confidence intervals can be estimated and contract risk can be assessed. Contract risk is defined as the probability that a block grade will exceed some specified limit. In coal mining, this specified limit will be set in a coal sales agreement. A key assumption necessary to implement the geostatistical model is that of local stationarity in the variogram. In a typical project, data limitations prevent a detailed examination of the stationarity assumption. In this paper, the distribution of kriging error and scale of variogram stationarity are examined for a coal property in northern West Virginia.     

A comparison between cokriging and ordinary kriging: Case study with a polymetallic deposit

Cokriging is applied to the estimation of mineral resources in a polymetallic deposit. Several major steps, which should be taken in using cokriging, are highlighted as necessary practical considerations. The case study is related to an ultramafic copper-nickel deposit. Six elements, Cu, Ni, Au, Ag, Pt, and Pd, occurring in the deposit, are partitioned into three subgroups and the elements within each group are simultaneously estimated on the basis of over 4000 drill assays. A comparison was made between ordinary kriging and cokriging methods through cross-validation. The results show that cokriging has significantly improved the estimates of resources by reducing the overall estimation error by over 15% and the variance of error by over 20%.     

Application of spatial filter theory to kriging

Data-processing requirements for remotely sensed, digital images include spatial filtering to suppress image noise, enhance edges/contacts, and improve image clarity. Spatial filter theory demonstrates that the addition of a high-pass filtered image to a low-pass filtered image yields the original digital image. Application of this principle in kriging can be accomplished by using the same covariance matrix to solve for two weighting vectors to yield a result analogous to low- and high-pass filtering. The addition of kriged estimates calculated using both weighting vectors is analogous to summing high-, and low-pass filtered digital images. This modified method of kriging yields estimates associated with less smoothing compared to ordinary kriging. Statistical moments of original sample data are better preserved through estimation by this method.     

Adding bounds to kriging

The ordinary kriging interpolation algorithm is extended by the inclusion of explicit lower and upper bounds on the estimate. The associated estimation variance is written as the ordinary kriging variance plus a non-negative correction term.     

Efficient updating of kriging estimates and variances

This short note presents a method for efficiently updating ordinary kriging estimates and variances when one or more additional samples are incorporated into the kriging system. First, the foundation linear algebra result is presented. Then the update equations are derived. Finally, an illustrative application of updating is briefly discussed.     

The variogram sill and the sample variance

The relationship between the sill of the variogram and the sample variance is explored. The common practice of using the sample variance as an estimate of the variogram sill is questioned, and a conceptual framework for determining the appropriateness of this heuristic is constructed.     

On unbiased backtransform of lognormal kriging estimates

Lognormal kriging is an estimation technique that was devised for handling highly skewed data distributions. This technique takes advantage of a logarithmic transformation that reduces the data variance. However, backtransformed lognormal kriging estimates are biased because the nonbias term is totally dependent on a semivariogram model. This paper proposes a new approach for backtransforming lognormal kriging estimates that not only presents none of the problems reported in the literature but also reproduces the sample histogram and, consequently, the sample mean.     

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.     

Determination of support in soil sampling

Short-range variation of measurements of a soil characteristic is a function of the support (i.e., size, shape, and orientation of the physical sample taken at a sample point) of the soil samples. This short-range variance often shows up as a substantial component of the nugget effect in experimental semivariograms or estimated generalized covariance functions. Proper choice of support may substantially reduce short-range variance (or nugget) of measurements made at sample points and thereby reduce estimation variance in point or block kriging based on these measurements. Appropriate variance formulas and an example are given for determination of support when support is an array of vertical cores (or bucket auger samples) of equal depth and diameter.