An Alternative Measure of the Reliability of Ordinary Kriging Estimates

This paper presents an interpolation variance as an alternative to the measure of the reliability of ordinary kriging estimates. Contrary to the traditional kriging variance, the interpolation variance is data-values dependent, variogram dependent, and a measure of local accuracy. Natural phenomena are not homogeneous; therefore, local variability as expressed through data values must be recognized for a correct assessment of uncertainty. The interpolation variance is simply the weighted average of the squared differences between data values and the retained estimate. Ordinary kriging or simple kriging variances are the expected values of interpolation variances; therefore, these traditional homoscedastic estimation variances cannot properly measure local data dispersion. More precisely, the interpolation variance is an estimate of the local conditional variance, when the ordinary kriging weights are interpreted as conditional probabilities associated to the n neighboring data. This interpretation is valid if, and only if, all ordinary kriging weights are positive or constrained to be such. Extensive tests illustrate that the interpolation variance is a useful alternative to the traditional kriging variance.     

土壤水分空间插值的克里金平滑效应修正方法

地统计学的普通克里金法是研究土壤水分空间变异特性和描绘其空间分布的有效方法。但与其它建立在最小二乘标准上的插值方法一样,普通克里金法也存在着平滑效应问题,即估计值的变异程度比实际要小,从而导致估计值往往不能反映出土壤水分真实的空间变化特征。结合实际的土壤水分监测数据,采用Yamamoto提出的一套针对普通克里金估计值进行后处理的方法,较好地解决了普通克里金法平滑效应的问题,在保证局部估计值精度的同时,重现了土壤含水率在空间的分布与变化特征。     

Correcting the Smoothing Effect of Estimators: A Spectral Postprocessor

The postprocessing algorithm introduced by Yao for imposing the spectral amplitudes of a target covariance model is shown to be efficient in correcting the smoothing effect of estimation maps, whether obtained by kriging or any other interpolation technique. As opposed to stochastic simulation, Yao's algorithm yields a unique map starting from an original, typically smooth, estimation map. Most importantly it is shown that reproduction of a covariance/semivariogram model (global accuracy) is necessarily obtained at the cost of local accuracy reduction and increase in conditional bias. When working on one location at a time, kriging remains the most accurate (in the least squared error sense) estimator. However, kriging estimates should only be listed, not mapped, since they do not reflect the correct (target) spatial autocorrelation. This mismatch in spatial autocorrelation can be corrected via stochastic simulation, or can be imposed a posteriori via Yao's algorithm.     

Robustness of Kriging weights to non-bias conditions

In this paper, the effect on Kriging weights of non-bias conditions, when the same residual covariance model is used, has been studied by the l₂ norm of the weights difference between Ordinary Kriging and Kriging with a trend model. Four covariance models, in 1-D and 2-D, and in interpolation and extrapolation conditions are examined. Situations in which both algorithms yield the same results are pointed out.     

拓扑克里格法与普通克里格法在区域洪水频率分析中的比较研究

依据北江(珠江流域支流)流域6个水文测站年最大洪峰流量资料,分别用Top-kriging(拓扑克里格法)和普通克里格法进行区域洪水频率估计。采用均方根误差作为频率分布线型拟合优度指标。运用线性矩法进行单站洪水频率分析,确定10、50、100、1000年一遇设计洪水值。在此基础上,从Topkriging和普通克里格法设计洪水估计不确定性和相对线性矩法单站洪水频率的估计误差两个方面比较Top-kriging和普通克里格法。结果表明:(1)Top-kriging法是更好的线性无偏估计,相比普通克里格法更适合于区域洪水频率估计;(2)Top-kriging法设计洪水估计不确定性明显小于普通克里格法;(3)Top-kriging法设计洪水估计结果更接近线性矩法单站洪水频率分析结果。     

Smoothing and interpolation by kriging and with splines

Let scalar measurements at distinct points x₁, , x_n be y₁, , y_n.We may look for a smooth function f(x)that goes through or near the points (x_i, y_i).Kriging assumes f(x)is a random function with known (possibly estimable) covariance function (in the simplest case). Splines assume a definition of the smoothness of a nonrandom function f(x).An elementary explanation is given of the fact that spline approximations are special cases of the solution of a kriging problem.     

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.     

Some limitations in the geostatistical evaluation of ore deposits

Conclusions The foregoing discussion indicates that geostatistical estimation of ore deposits is not local; it is not objective; it is not sensitive to local data trends; and it is not unrestrained by the range of data values.Kriging, as an interpolation method, is a variant of IDW least squares linear fit. As such, it suffers from the limitations of all IDW linear interpolation methods that employ only data values.The estimation variance, currently used to calculate the confidence limits of values for individual mining blocks, is hypothetical and globally derived. It is more closely related to sampling density than to local variation in the data set.Geostatistical methods, of course, have a real place in ore deposit assessment, e.g. global, comparative evaluation to assist decisions on development and investment. What is questioned here is the validity of employing a global method to assess detail (mining blocks) within an ore deposit.     

An Experimental Comparison of Ordinary and Universal Kriging and Inverse Distance Weighting

A factorial, computational experiment was conducted to compare the spatial interpolation accuracy of ordinary and universal kriging and two types of inverse squared-distance weighting. The experiment considered, in addition to these four interpolation methods, the effects of four data and sampling characteristics: surface type, sampling pattern, noise level, and strength of small-scale spatial correlation. Interpolation accuracy was measured by the natural logarithm of the mean squared interpolation error. Main effects of all five factors, all two-factor interactions, and several three-factor interactions were highly statistically significant. Among numerous findings, the most striking was that the two kriging methods were substantially superior to the inverse distance weighting methods over all levels of surface type, sampling pattern, noise, and correlation.     

Evaluation of Interpolation Accuracy of Neural Kriging with Application to Temperature-Distribution Analysis

An interpolation method based on a multilayer neural network (MNN), has been examined and tested for the data of irregular sample locations. The main advantage of MNN is in that it can deal with geoscience data with nonlinear behavior and extract characteristics from complex and noisy images. The training of MNN is used to modify connection weights between nodes located in different layers by a simulated annealing algorithm (one of the optimization algorithms of the network). In this process, three types of errors are considered: differences in values, semivariograms, and gradients between sample data and outputs from the trained network. The training is continued until the summation of these errors converges to an acceptably small value. Because the MNN trained by this learning criterion can estimate a value at an arbitrary location, this method is a form of kriging and termed Neural Kriging (NK). In order to evaluate the effectiveness of NK, a problem on restoration ability of a defined reference surface from randomly chosen discrete data was prepared. Two types of surfaces, whose semivariograms are expressed by isotropic spherical and geometric anisotropic gaussian models, were examined in this problem. Though the interpolation accuracy depended on the arrangement pattern of the sample locations for the same number of data, the interpolation errors of NK were shown to be smaller than both those of ordinary MNN and ordinal kriging. NK can also produce a contour map in consideration of gradient constraints. Furthermore, NK was applied to distribution analysis of subsurface temperatures using geothermal investigation loggings of the Hohi area in southwest Japan. In spite of the restricted quantity of sample data, the interpolation results revealed high temperature zones and convection patterns of hydrothermal fluids. NK is regarded as an interpolation method with high accuracy that can be used for regionalized variables with any structure of spatial correlation.     

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

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

On the Equivalence of the Cokriging and Kriging Systems

Simple cokriging of components of a p-dimensional second-order stationary random process is considered. Necessary and sufficient conditions under which simple cokriging is equivalent to simple kriging are given. Essentially this condition requires that it should be possible to express the cross-covariance at any lag series h using the cross-covariance at |h|=0 and the auto-covariance at lag series h. The mosaic model, multicolocated kriging and the linear model of coregionalization are examined in this context. A data analytic method to examine whether simple kriging of components of a multivariate random process is equivalent to its cokriging is given     

Dual Kriging with Local Neighborhoods: Application to the Representation of Surfaces

Ordinary kriging, in its common formulation, is a discrete estimator in that it requires the solution of a kriging system for each point in space in which an estimate is sought. The dual formulation of ordinary kriging provides a continuous estimator since, for a given set of data, only a kriging system has to be estimated and the resulting estimate is a function continuously defined in space. The main problem with dual kriging up to now has been that its benefits can only be capitalized if a global neighborhood is used. A formulation is proposed to solve the problem of patching together dual kriging estimates obtained with data from different neighborhoods by means of a blending belt around each neighborhood. This formulation ensures continuity of the variable and, if needed, of its first derivative along neighbor borders. The final result is an analytical formulation of the interpolating surface that can be used to compute gradients, cross-sections, or volumes; or for the quick evaluation of the interpolating surface in numerous locations.     

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 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.     

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.     

指数平滑技术在斜坡位移预测中的应用

应用趋势型数据指数平滑模型预测了链子崖危岩体监测点GA的位移量.根据其观测数据呈线性趋势的特点,选取趋势型二次指数平滑的线性预测公式和合适的平滑常数值进行了计算预测.预测结果表明位移预测值与实际观测值之间的误差很小,说明该模型可很好地应用于斜坡变形位移的预测.     

Comparison of Kriging and Neural Networks With Application to the Exploitation of a Slate Mine

To carry out an efficient and effective exploitation of a slate mine, it is necessary to have detailed information about the production potential of the site. To assist us in estimating the quality of slate from a small set of drilling data within an unexploited portion of the mine, the following estimation techniques were applied: kriging, regularization networks (RN), multilayer perceptron (MLP) networks, and radial basis function (RBF) networks. Our numerical results for the test holes show that the best results were obtained using an RN (kriging) which takes into account the known anisotropy. Differing deposit configurations were obtained, depending on the method applied. Variations in the form of pockets were obtained when using a radial pattern with RBF, RN, and kriging models while a stratified pattern was obtained with the MLP model. Pockets are more suitable for a slate mine, which indicates that the selection of a technique should take account of the specific configuration of the deposit according to mineral type.     

Kriging法在区域滑坡危险性评价中的应用

区域滑坡危险性评价是进行区域滑坡风险性研究的基础.由于滑坡演变机制的复杂性,使得目前基于独立分析各因素对滑坡影响的“白箱”型评价模式具有一定风险性,同时这类评价方法要求对滑坡演变和研究区地质地理背景进行非常细致的监测和调查.为了克服这些问题,文章提出了一种基于Kriging插值理论的“黑箱”型评价方法.在利用该方法对历史滑坡点的规模进行评价的基础上,利用Kriging插值法获取研究区的滑坡危险性区划,并以四川省苍溪县为例,验证了运用该方法进行区域滑坡危险性评价的可行性.     

基于灰色系统的黄龙钙华景观演化研究

黄龙的景观是在数万年的岩溶地质作用下的产物。近年来黄龙钙华出现了干涸、变黑、沙化等现象严重地影响了景观的观赏性。水资源在黄龙钙华发展变化中起着重要的作用，笔者通过收集并分析黄龙地区的监测数据，针对监测系统不健全．时序数据缺乏的特点，选用对时序数据要求不高、预测效果较好的灰色系统模型，以岩溶水体的pH值为指标预测了钙华未来的发展情况。采用地质统计学空间分析的克里金插值法，对预测结果进行插值获得了整个景区的钙华预测结果。笔者按pH值将钙华演化情况划为强侵蚀、弱侵蚀、堆积3种类型，指出黄龙钙华景观目前正处于动态平衡与消亡重组阶段。