首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
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.  相似文献   

2.
Conditioning realizations of stationary Gaussian random fields to a set of data is traditionally based on simple kriging. In practice, this approach may be demanding as it does not account for the uncertainty in the spatial average of the random field. In this paper, an alternative model is presented, in which the Gaussian field is decomposed into a random mean, constant over space but variable over the realizations, and an independent residual. It is shown that, when the prior variance of the random mean is infinitely large (reflecting prior ignorance on the actual spatial average), the realizations of the Gaussian random field are made conditional by substituting ordinary kriging for simple kriging. The proposed approach can be extended to models with random drifts that are polynomials in the spatial coordinates, by using universal or intrinsic kriging for conditioning the realizations, and also to multivariate situations by using cokriging instead of kriging.  相似文献   

3.
Kriging Prediction Intervals Based on Semiparametric Bootstrap   总被引:1,自引:0,他引:1  
Kriging is a widely used method for prediction, which, given observations of a (spatial) process, yields the best linear unbiased predictor of the process at a new location. The construction of corresponding prediction intervals typically relies on Gaussian assumptions. Here we show that the distribution of kriging predictors for non-Gaussian processes may be far from Gaussian, even asymptotically. This emphasizes the need for other ways to construct prediction intervals. We propose a semiparametric bootstrap method with focus on the ordinary kriging predictor. No distributional assumptions about the data generating process are needed. A simulation study for Gaussian as well as lognormal processes shows that the semiparametric bootstrap method works well. For the lognormal process we see significant improvement in coverage probability compared to traditional methods relying on Gaussian assumptions.  相似文献   

4.
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.  相似文献   

5.
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.  相似文献   

6.
The numerical stability of linear systems arising in kriging, estimation, and simulation of random fields, is studied analytically and numerically. In the state-space formulation of kriging, as developed here, the stability of the kriging system depends on the condition number of the prior, stationary covariance matrix. The same is true for conditional random field generation by the superposition method, which is based on kriging, and the multivariate Gaussian method, which requires factoring a covariance matrix. A large condition number corresponds to an ill-conditioned, numerically unstable system. In the case of stationary covariance matrices and uniform grids, as occurs in kriging of uniformly sampled data, the degree of ill-conditioning generally increases indefinitely with sampling density and, to a limit, with domain size. The precise behavior is, however, highly sensitive to the underlying covariance model. Detailed analytical and numerical results are given for five one-dimensional covariance models: (1) hole-exponential, (2) exponential, (3) linear-exponential, (4) hole-Gaussian, and (5) Gaussian. This list reflects an approximate ranking of the models, from best to worst conditioned. The methods developed in this work can be used to analyze other covariance models. Examples of such representative analyses, conducted in this work, include the spherical and periodic hole-effect (hole-sinusoidal) covariance models. The effect of small-scale variability (nugget) is addressed and extensions to irregular sampling schemes and higher dimensional spaces are discussed.  相似文献   

7.
 A thorough understanding of the characteristics of transmissivity makes groundwater deterministic models more accurate. These transmissivity data characteristics occasionally possess a complicated spatial variation over an investigated site. This study presents both geostatistical estimation and conditional simulation methods to generate spatial transmissivity maps. The measured transmissivity data from the Dulliu area in Yun-Lin county, Taiwan, is used as the case study. The spatial transmissivity maps are simulated by using sequential Gaussian simulation (SGS), and estimated by using natural log ordinary kriging and ordinary kriging. Estimation and simulation results indicate that SGS can reproduce the spatial structure of the investigated data. Furthermore, displaying a low spatial variability does not allow the ordinary kriging and natural log kriging estimates to fit the spatial structure and small-scale variation for the investigated data. The maps of kriging estimates are smoother than those of other simulations. A SGS with multiple realizations has significant advantages over ordinary kriging and even natural log kriging techniques at a site with a high variation in investigated data. These results are displayed in geographic information systems (GIS) as basic information for further groundwater study. Received: 27 August 1999 · Accepted: 22 February 2000  相似文献   

8.
This paper studies vector (multivariate, multiple, or multidimensional) random fields in space and/or time with second-order increments, for which the variogram matrix is an important tool to measure the dependence within each component and between each pair of distinct components. We introduce an efficient approach to construct Gaussian or non-Gaussian vector random fields from the univariate random field with higher dimensional index domain, and particularly to generate a class of variogram matrices.  相似文献   

9.
Looking at kriging problems with huge numbers of estimation points and measurements, computational power and storage capacities often pose heavy limitations to the maximum manageable problem size. In the past, a list of FFT-based algorithms for matrix operations have been developed. They allow extremely fast convolution, superposition and inversion of covariance matrices under certain conditions. If adequately used in kriging problems, these algorithms lead to drastic speedup and reductions in storage requirements without changing the kriging estimator. However, they require second-order stationary covariance functions, estimation on regular grids, and the measurements must also form a regular grid. In this study, we show how to alleviate these rather heavy and many times unrealistic restrictions. Stationarity can be generalized to intrinsicity and beyond, if decomposing kriging problems into the sum of a stationary problem and a formally decoupled regression task. We use universal kriging, because it covers arbitrary forms of unknown drift and all cases of generalized covariance functions. Even more general, we use an extension to uncertain rather than unknown drift coefficients. The sampling locations may now be irregular, but must form a subset of the estimation grid. Finally, we present asymptotically exact but fast approximations to the estimation variance and point out application to conditional simulation, cokriging and sequential kriging. The drastic gain in computational and storage efficiency is demonstrated in test cases. Especially high-resolution and data-rich fields such as rainfall interpolation from radar measurements or seismic or other geophysical inversion can benefit from these improvements.  相似文献   

10.
This paper provides a comparison between linear (universal) and nonlinear (disjunctive) kriging estimators when they are computed from small samples chosen randomly on simulated stationary and nonstationary fields. Point estimation results are reported. In all cases considered, kriging estimators were found better than a local mean estimator, with universal kriging either better than or as good as disjunctive kriging. The latter, which is suited to handle stationary fields, did not provide more accurate estimates because the use of small samples led to inconsistencies in the assumed bivariate model. Universal kriging was particularly better with nonstationary fields.  相似文献   

11.
在冲积含水层中,由于岩相的非均质分布,渗透系数一般呈现出明显的非高斯特性(例如砂和黏土两种岩相),非高斯特性给地下水模型参数的推估带来了困难与挑战。目前广泛使用的集合平滑数据同化方法(ESMDA)虽然有效且计算成本较低,但仅适用于高斯场。多点地质统计方法虽已广泛用于模拟非高斯场,但其无法融入动态观测数据推估参数。基于多点地质统计方法中的直接采样法(DS)与集合平滑数据同化方法,构建一种新的数据同化框架(ESMDA-DS),既可保持参数场的非高斯特性,又可融合多源数据精确推估非高斯场。构建三个理想算例验证ESMDA-DS对非高斯参数场的推估效果,并探讨了不同类型观测数据对推估效果、水位与浓度预测精度的影响。三个理想算例包括仅融合水位数据(Case 1),同时融合水位与浓度数据(Case 2),同时融合水位、浓度与对数渗透系数数据(Case 3)。结果表明:ESMDA-DS方法结合了ESMDA与DS的各自优势,能有效融合观测数据推估渗透系数场,并保持参数场的非高斯特性。通过对比三个算例推估结果,Case 3的参数场推估效果最好,水位与浓度预测精度最高,Case 2次之,Case 1最差,表明融合多源数据可改善推估效果,提高预测精度。  相似文献   

12.
In this study, machine learning methods such as neural networks, random forests, and Gaussian processes are applied to the estimation of copper grade in a mineral deposit. The performance of these methods is compared to geostatistical techniques, such as ordinary kriging and indicator kriging. To ensure that these comparisons are realistic and relevant, the predictive accuracy is estimated on test instances located in drill holes that are different from the training data. The results of an extensive empirical study in the Sarcheshmeh porphyry copper deposit in Southeastern Iran illustrate that specially designed Gaussian processes with a symmetric standardization of the spatial location inputs and an anisotropic kernel yield the most accurate predictions. Furthermore, significant improvements are obtained when, besides location, information on the rock type is included in the set of predictor variables. This observation highlights the importance of carrying out detailed studies of the geological composition of the deposit to obtain more accurate ore grade predictions.  相似文献   

13.
Estimation or simulation? That is the question   总被引:1,自引:0,他引:1  
The issue of smoothing in kriging has been addressed either by estimation or simulation. The solution via estimation calls for postprocessing kriging estimates in order to correct the smoothing effect. Stochastic simulation provides equiprobable images presenting no smoothing and reproducing the covariance model. Consequently, these images reproduce both the sample histogram and the sample semivariogram. However, there is still a problem, which is the lack of local accuracy of simulated images. In this paper, a postprocessing algorithm for correcting the smoothing effect of ordinary kriging estimates is compared with sequential Gaussian simulation realizations. Based on samples drawn from exhaustive data sets, the postprocessing algorithm is shown to be superior to any individual simulation realization yet, at the expense of providing one deterministic estimate of the random function.  相似文献   

14.
Universal kriging is compared with ordinary kriging for estimation of earthquake ground motion. Ordinary kriging is based on a stationary random function model; universal kriging is based on a nonstationary random function model representing first-order drift. Accuracy of universal kriging is compared with that for ordinary kriging; cross-validation is used as the basis for comparison. Hypothesis testing on these results shows that accuracy obtained using universal kriging is not significantly different from accuracy obtained using ordinary kriging. Tests based on normal distribution assumptions are applied to errors measured in the cross-validation procedure;t andF tests reveal no evidence to suggest universal and ordinary kriging are different for estimation of earthquake ground motion. Nonparametric hypothesis tests applied to these errors and jackknife statistics yield the same conclusion: universal and ordinary kriging are not significantly different for this application as determined by a cross-validation procedure. These results are based on application to four independent data sets (four different seismic events).  相似文献   

15.
边少锋  Menz.J 《地球科学》2000,25(2):195-200
首先引入利用旋转面作为基函数的函数逼近概念, 在此基础上经过复杂的矩阵推导证明泛克立格法可表示为传统的带权最小二乘多项式拟合与以旋转面作为基函数的函数逼近, 并在一定条件下(随机场高度连续无块金效应) 论证了协方差(即旋转面) 的参数可通过数学分析的方法确定, 给出了以高斯函数为例确定协方差函数的两个准则.   相似文献   

16.
Compensating for estimation smoothing in kriging   总被引:2,自引:0,他引:2  
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-ft2 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.  相似文献   

17.
For earthquake ground motion studies, the actual ground motion distribution should be reproduced as accurately as possible. For optimal estimation of ground motion, kriging has been shown to provide accurate estimates. Although kriging is accurate for this application, some estimates it provides are underestimates. This has dire consequences for subsequent design for earthquake resistance. Kriging does not provide enough information to allow an analysis of each estimate for underestimation. For such an application, disjunctive kriging is better applied. This advanced technique quantifies the probability that an estimate equals or exceeds particular levels of ground motion. Furthermore, disjunctive kriging can provide improved estimation accuracy when applied for local estimation of ground motion.  相似文献   

18.
An important aspect in mineral resource evaluation is the reduction of variance when post-processing the grade distributions defined on the support (volume) of the available data into distributions defined on the support of the proposed selective mining units. Although the volume-variance relationship is well understood for the estimation of global grade distributions, it is still an unsolved issue for local estimation studies based on non-parametric geostatistical methods, such as indicator kriging, for which the support correction is not inherent to the method. To clarify this relationship, the local change of support problem is examined in the scope of two parametric models (multi-Gaussian and discrete Gaussian models). It is shown that the variance reduction factor between point and block-support local distributions depends on the block being considered and is less than the global variance reduction factor. As a consequence, post-processing the local point-support grade distributions on the basis of the latter systematically understates the importance of the change of support at the local scale and makes selective mining appear more economically attractive than it really is. In the light of these results, a methodology is proposed to post-process the local point-support distributions obtained via non-parametric (indicator) methods into block-support distributions. An application to simulated data indicates that this methodology provides an accurate estimation at the block support when dealing with diffusion-type random fields.  相似文献   

19.
Most significant iron ore deposits in Iran are located in Central Iran Zone. These deposits belong to the Bafq mining district. The Bafq mining district is located in the Early Cambrian Kashmar-Kerman volcanic arc of Central Iran. Linear estimation of regionalized variables (for example by inverse distance weighting or ordinary Kriging) results in relatively high estimation variances, i.e. the estimates have very low precision. Assessment of project economics (or other critical decision making) based on linear estimation is therefore risky. Non-linear estimation methods like disjunctive kriging perform better and the lower estimation variance allows less risky economic decision-making. Another advantage of disjunctive kriging is that it allows estimation of functions of the primary variable, which here is the grade (Fe %) of the ore. In particular it permits estimation of indicator functions defined using thresholds on the primary variable. This paper is devoted to application of disjunctive kriging method in Choghart North Anomaly iron ore deposit in Central Iran, Yazd province, Iran. In this study, the Fe concentration of Choghart North Anomaly iron ore deposit was modelled and estimated. The exploration data consists of borehole samples measuring the Fe concentration. A Gaussian isofactorial model is fitted to these data and disjunctive kriging was used to estimate the regionalized variable (Fe %) at unsampled locations and to assess the probabilities that the actual concentrations exceed a threshold value at a given location. Consequently a three dimensional model of probability of exceeding a threshold value and the estimated value are provided by disjunctive kriging to divide the ore into an economic and uneconomic part on the basis of estimation of indicator functions using thresholds grades defined on point support. The tools and concepts are complemented by a set of computer programs that are applied to the case study. The study showed that disjunctive kriging can be applied successfully for modeling the grade of an ore deposit. Results showed that the correlation between the estimated value and real value at locations close to each other is 81.9%.  相似文献   

20.
Environmental, engineering and industrial modelling of natural features (e.g. trees) and man-made features (e.g. pipelines) requires some form of fitting of geometrical objects such as cylinders, which is commonly undertaken using a least-squares method that—in order to get optimal estimation—assumes normal Gaussian distribution. In the presence of outliers, however, this assumption is violated leading to a Gaussian mixture distribution. This study proposes a robust parameter estimation method, which is an improved and extended form of vector algebraic modelling. The proposed method employs expectation maximisation and maximum likelihood estimation (MLE) to find cylindrical parameters in case of Gaussian mixture distribution. MLE computes the model parameters assuming that the distribution of model errors is a Gaussian mixture corresponding to inlier and outlier points. The parameters of the Gaussian mixture distribution and the membership functions of the inliers and outliers are computed using an expectation maximisation algorithm from the histogram of the model error distribution, and the initial guess values for the model parameters are obtained using total least squares. The method, illustrated by a practical example from a terrestrial laser scanning point cloud, is novel in that it is algebraic (i.e. provides a non-iterative solution to the global maximisation problem of the likelihood function), is practically useful for any type of error distribution model and is capable of separating points of interest and outliers.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号