共查询到20条相似文献，搜索用时 46 毫秒
1.
Sample schemes used in geostatistical surveys must be suitable for both variogram estimation and kriging. Previously schemes
have been optimized for one of these steps in isolation. Ordinary kriging generally requires the sampling locations to be
evenly dispersed over the region. Variogram estimation requires a more irregular pattern of sampling locations since comparisons
must be made between measurements separated by all lags up to and beyond the range of spatial correlation. Previous studies
have not considered how to combine these optimized schemes into a single survey and how to decide what proportion of sampling
effort should be devoted to variogram estimation and what proportion devoted to kriging
An expression for the total error in a geostatistical survey accounting for uncertainty due to both ordinary kriging and variogram
uncertainty is derived. In the same manner as the kriging variance, this expression is a function of the variogram but not
of the sampled response data. If a particular variogram is assumed the total error in a geostatistical survey may be estimated
prior to sampling. We can therefore design an optimal sample scheme for the combined processes of variogram estimation and
ordinary kriging by minimizing this expression. The minimization is achieved by spatial simulated annealing. The resulting
sample schemes ensure that the region is fairly evenly covered but include some close pairs to analyse the spatial correlation
over short distances. The form of these optimal sample schemes is sensitive to the assumed variogram. Therefore a Bayesian
approach is adopted where, rather than assuming a single variogram, we minimize the expected total error over a distribution
of plausible variograms. This is computationally expensive so a strategy is suggested to reduce the number of computations
required 相似文献
2.
Gilles Bourgault 《Mathematical Geosciences》2014,46(7):841-868
A multivariate probability transformation between random variables, known as the Nataf transformation, is shown to be the appropriate transformation for multi-Gaussian kriging. It assumes a diagonal Jacobian matrix for the transformation of the random variables between the original space and the Gaussian space. This allows writing the probability transformation between the local conditional probability density function in the original space and the local conditional Gaussian probability density function in the Gaussian space as a ratio equal to the ratio of their respective marginal distributions. Under stationarity, the marginal distribution in the original space is modeled from the data histogram. The stationary marginal standard Gaussian distribution is obtained from the normal scores of the data and the local conditional Gaussian distribution is modeled from the kriging mean and kriging variance of the normal scores of the data. The equality of ratios of distributions has the same form as the Bayes’ rule and the assumption of stationarity of the data histogram can be re-interpreted as the gathering of the prior distribution. Multi-Gaussian kriging can be re-interpreted as an updating of the data histogram by a Gaussian likelihood. The Bayes’ rule allows for an even more general interpretation of spatial estimation in terms of equality for the ratio of the conditional distribution over the marginal distribution in the original data uncertainty space with the same ratio for a model of uncertainty with a distribution that can be modeled using the mean and variance from direct kriging of the original data values. It is based on the principle of conservation of probability ratio and no transformation is required. The local conditional distribution has a variance that is data dependent. When used in sequential simulation mode, it reproduces histogram and variogram of the data, thus providing a new approach for direct simulation in the original value space. 相似文献
3.
以浅剖数据为源数据,钻孔实测数据为验证数据,利用普通克里金法对海底地层厚度进行空间插值得到地层分布特征,采用3种半变异函数模型和不同取样间距对某井场3组地层厚度进行普通克里金插值并验证其插值效果。结果表明:普通克里金是一种有效的海底地层厚度预测方法;结构分析最佳的模型不一定是误差最小的模型,应对不同模型下的插值结果进行综合分析来选择最合适的模型,并提出球状模型在该井场厚度估计中最优,高斯模型次之;对于球状模型,增大取样间距对地层厚度变化剧烈的地层回归效果影响较小,对地层厚度变化不大的地层回归效果影响较大;同时,SE预测值变化率分析表明对于地层厚度变化剧烈的地层,减小取样间距可以大幅度地减少插值误差,而对于地层厚度变化不大的地层,减小取样间距对插值精度提高的意义不大。 相似文献
4.
On Visualization for Assessing Kriging Outcomes 总被引：7，自引：0，他引：7
James R. Carr 《Mathematical Geology》2002,34(4):421-433
Extant opinion about kriging is that all weights should be positive. Visualizations rendered by converting kriged grids to digital images are presented to show that negative weights may be beneficial to some spatial problems. In particular, variogram models with zero-valued nuggets, already well known to minimize smoothing through kriging, result in a visual resolution substantially superior to that from kriging with a variogram model having a nonzero nugget value in application to satellite acquired data. Negative weights are more likely when using variogram models with zero-valued nuggets, but resultant visualizations often show a smoother transition between extreme data values. This is true even when a variogram model having a nugget value of zero is not optimum with respect to mean square error, as is demonstrated using a nitrate data set. An analogy to digital image processing is used to suggest that the influence of negative weights in kriging is similar to a high-boost kernel. 相似文献
5.
Subhash Lele 《Mathematical Geology》1995,27(5):673-692
Two important problems in the practical implementation of kriging are: (1) estimation of the variogram, and (2) estimation of the prediction error. In this paper, a nonparametric estimator of the variogram to circumvent the problem of the precise choice of a variogram model is proposed. Using orthogonal decomposition of the kriging predictor and the prediction error, a method for selecting, what may be considered, a statistical neighborhood is suggested. The prediction error estimates based on this scheme, in fact, reflects the true prediction error, thus leading to proper coverage for the corresponding prediction interval. By simulations and a reanalysis of published data, it is shown that the proposals made in this paper are useful in practice. 相似文献
6.
Comparison of kriging techniques in a space-time context 总被引：1，自引：0，他引：1
Patrick Bogaert 《Mathematical Geology》1996,28(1):73-86
Space-time processes constitute a particular class, requiring suitable tools in order to predict values in time and space, such as a space-time variogram or covariance function. The space-time co-variance function is defined and linked to the Linear Model of Coregionalization under second-order space-time stationarity. Simple and ordinary space-time kriging systems are compared to simple and ordinary cokriging and their differences for unbiasedness conditions are underlined. The ordinary space-time kriging estimation then is applied to simulated data. Prediction variances and prediction errors are compared with those for ordinary kriging and cokriging under different unbiasedness conditions using a cross-validation. The results show that space-time kriging tend to produce lower prediction variances and prediction errors that kriging and cokriging. 相似文献
7.
The Second-Order Stationary Universal Kriging Model Revisited 总被引：3，自引：0，他引：3
Universal kriging originally was developed for problems of spatial interpolation if a drift seemed to be justified to model the experimental data. But its use has been questioned in relation to the bias of the estimated underlying variogram (variogram of the residuals), and furthermore universal kriging came to be considered an old-fashioned method after the theory of intrinsic random functions was developed. In this paper the model is reexamined together with methods for handling problems in the inference of parameters. The efficiency of the inference of covariance parameters is shown in terms of bias, variance, and mean square error of the sampling distribution obtained by Monte Carlo simulation for three different estimators (maximum likelihood, bias corrected maximum likelihood, and restricted maximum likelihood). It is shown that unbiased estimates for the covariance parameters may be obtained but if the number of samples is small there can be no guarantee of good estimates (estimates close to the true value) because the sampling variance usually is large. This problem is not specific to the universal kriging model but rather arises in any model where parameters are inferred from experimental data. The validity of the estimates may be evaluated statistically as a risk function as is shown in this paper. 相似文献
8.
含水层渗透性空间分布的指示克立格估值 总被引：3，自引：0，他引：3
详细介绍了指示克立格估值计算的理论和方法。以指示变异函数为基本工具分析了华北某地区第四系含水层渗透性空间分布的结构特征,结果表明该地区含水层渗透性存在明显的各向异性特征。水平方向上,X轴方向的相关性较Y轴方向的好,Z轴的相关性最差。用指示克立格法对未采样点处进行估值,估值结果显示含水层渗透性由山前向滨海逐渐变低,在垂直方向上,渗透性变化不明显,浅部比深部略好;同时给出了估计精度,并认为对估计精度不高的区域可通过增加适当的工程加以控制。最后用交叉验证法对估值结果进行了检验,证明建立的指示变异函数模型合理且估值效果较好。这一实际应用表明指示克立格法可以很好地描述第四系含水层渗透性的空间分布规律。 相似文献
9.
Assessment of the sampling variance of the experimental variogram is an important topic in geostatistics as it gives the uncertainty of the variogram estimates. This assessment, however, is repeatedly overlooked in most applications mainly, perhaps, because a general approach has not been implemented in the most commonly used software packages for variogram analysis. In this paper the authors propose a solution that can be implemented easily in a computer program, and which, subject to certain assumptions, is exact. These assumptions are not very restrictive: second-order stationarity (the process has a finite variance and the variogram has a sill) and, solely for the purpose of evaluating fourth-order moments, a Gaussian distribution for the random function. The approach described here gives the variance–covariance matrix of the experimental variogram, which takes into account not only the correlation among the experiemental values but also the multiple use of data in the variogram computation. Among other applications, standard errors may be attached to the variogram estimates and the variance–covariance matrix may be used for fitting a theoretical model by weighted, or by generalized, least squares. Confidence regions that hold a given confidence level for all the variogram lag estimates simultaneously have been calculated using the Bonferroni method for rectangular intervals, and using the multivariate Gaussian assumption for K-dimensional elliptical intervals (where K is the number of experimental variogram estimates). A general approach for incorporating the uncertainty of the experimental variogram into the uncertainty of the variogram model parameters is also shown. A case study with rainfall data is used to illustrate the proposed approach. 相似文献
10.
Denis Marcotte 《Mathematical Geology》1995,27(6):749-762
Conditional simulation with data subject to measurement error has received little attention in the geostatistical literature. The treatment of measurement error in simulation must be different from its treatment in estimation. Two approaches are examined: pre- and post-simulation filtering of data measurement error. The pre-simulation filtering is shown to be inefficient. The post-simulation filtering performs best. It is done by factorial kriging and a modified version of factorial kriging which ensures predetermined theoretical variance for the filtered data. It also is shown that the theoretical variogram of the filtered data reproduces the underlying variogram (i.e., without noise) almost perfectly. A simulation with a high level of correlated noise is used for validation and comparison. The post-simulation filtered values show an experimental variogram in agreement with the previously identified underlying variogram. Moreover, the filtered image compares well with the true image. The theoretical variogram corresponding to the post-simulation filter can be computed beforehand. Thus, the size of the simulation grid and of the filter neighborhood can be adjusted to ensure good reproduction of the underlying variogram. 相似文献
11.
Comparing the robustness of ordinary kriging and lognormal kriging: Outlier resistance 总被引：1，自引：0，他引：1
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. 相似文献
12.
Saeed Soltani-Mohammadi 《Earth Science Informatics》2016,9(2):235-245
Using kriging has been accepted today as the most common method of estimating spatial data in such different fields as the geosciences. To be able to apply kriging methods, it is necessary that the data and variogram model parameters be precise. To utilize the imprecise (fuzzy) data and parameters, use is made of fuzzy kriging methods. Although it has been 30 years since different fuzzy kriging algorithms were proposed, its use has not become as common as other kriging methods (ordinary, simple, log, universal, etc.); lack of a comprehensive software that can perform, based on different fuzzy kriging algorithms, the related calculations in a 3D space can be the main reason. This paper describes an open-source software toolbox (developed in Matlab) for running different algorithms proposed for fuzzy kriging. It also presents, besides a short presentation of the fuzzy kriging method and introduction of the functions provided by the FuzzyKrig toolbox, 3 cases of the software application under the conditions where: 1) data are hard and variogram model parameters are fuzzy, 2) data are fuzzy and variogram model parameters are hard, and 3) both data and variogram model parameters are fuzzy. 相似文献
13.
Robustness of variograms and conditioning of kriging matrices 总被引：1，自引：0，他引：1
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. 相似文献
14.
A. G. Journel 《Mathematical Geology》1988,20(4):459-475
The projection or minimum error norm algorithm does not require that the distance measure be a variogram. In non-Gaussian cases, the traditional variogram distance measure leading to minimization of an error variance offers no definite advantage. Other distance measures, more outlierresistant than the variogram, are proposed which fulfill the condition of the projection theorem. The resulting minimum error norms provide the same data configurations ranking as traditionally obtained from kriging variances. A case study based on actual digital terrain data is presented.This paper was presented by title at MGUS 87 Conference, Redwood City, California, 14 April 1987. 相似文献
15.
16.
Geological data frequently have a heavy-tailed normal-in-the-middle distribution, which gives rise to grade distributions that appear to be normal except for the occurrence of a few outliers. This same situation also applies to log-transformed data to which lognormal kriging is to be applied. For such data, linear kriging is nonrobust in that (1)kriged estimates tend to infinity as the outliers do, and (2)it is also not minimum mean squared error. The more general nonlinear method of disjunctive kriging is even more nonrobust, computationally more laborious, and in the end need not produce better practical answers. We propose a robust kriging method for such nearly normal data based on linear kriging of an editing of the data. It is little more laborious than conventional linear kriging and, used in conjunction with a robust estimator of the variogram, provides good protection against the effects of data outliers. The method is also applicable to time series analysis. 相似文献
17.
Notes on the robustness of the kriging system 总被引：3，自引：0，他引：3
Andras Bardossy 《Mathematical Geology》1988,20(3):189-203
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. 相似文献
18.
Geological data frequently have a heavy-tailed normal-in-the-middle distribution, which gives rise to grade distributions that appear to be normal except for the occurrence of a few outliers. This same situation also applies to log-transformed data to which lognormal kriging is to be applied. For such data, linear kriging is nonrobust in that (1)kriged estimates tend to infinity as the outliers do, and (2)it is also not minimum mean squared error. The more general nonlinear method of disjunctive kriging is even more nonrobust, computationally more laborious, and in the end need not produce better practical answers. We propose a robust kriging method for such nearly normal data based on linear kriging of an editing of the data. It is little more laborious than conventional linear kriging and, used in conjunction with a robust estimator of the variogram, provides good protection against the effects of data outliers. The method is also applicable to time series analysis. 相似文献
19.
Changik Han Jiyang Wang Mingguo Zheng Ende Wang Jianming Xia GwangSu Li Sunchol Choe 《Earth Science Informatics》2016,9(2):197-213
A critical step for kriging in geostatistics is estimation of the variogram. Traditional variogram modeling comprise of the experimental variogram calculation, appropriate variogram model selection and model parameter determination. Selecting of the variogram model and fitting of model parameters is the most controversial aspect of geostatistics. Shapes of valid variogram models are finite, and sometimes, the optimal shape of the model can not be fitted, leading to reduced estimation accuracy. In this paper, a new method is presented to automatically construct a model shape and fit model parameters to experimental variograms using Support Vector Regression (SVR) and Multi-Gene Genetic Programming (MGGP). The proposed method does not require the selection of a variogram model and can directly provide the model shape and parameters of the optimal variogram. The validity of the proposed method is demonstrated in a number of cases. 相似文献
20.
应用Kriging方法研究格尔木河流域地下水位动态观测网的优化配置 总被引：4，自引：0，他引：4
地下水动态是水文地质研究的一个重要方面，如何获满足一定精度要求全面而合理的地下水位动态资料具有现实经济效益。本文研究格尔木河流域地下水位动态测网，通过选取理想的变差函，运有Ｋｒｉｇｉｎｇ方法进行地下水位的线性无偏最优估计和计算估计误差的标准差，结合给定的允许误差限评判地下水位动态观测网的配置是否合适，并提出调整现有观测网的方案 相似文献