共查询到20条相似文献，搜索用时 750 毫秒
1.
Reducing the Impact of Outliers in Ore Reserves Estimation 总被引：1，自引：0，他引：1
Joao Felipe Costa 《Mathematical Geology》2003,35(3):323-345
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. 相似文献
2.
Application of universal kriging for estimation of earthquake ground motion: Statistical significance of results 总被引：1，自引：0，他引：1
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). 相似文献
3.
Xavier Emery 《Mathematical Geology》2007,39(6):607-623
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. 相似文献
4.
Guocheng Pan 《Mathematical Geology》1994,26(1):135-155
Restricted kriging provides a simple and quick remedy for the problem known as the weight independence of data in ordinary kriging. A major consequence of this problem is the effect of over-smearing in estimates, which, in turn, adds one uncertain factor to subsequent mine decisions. A detailed count is reported here on a restricted kriging system that incorporates two restrictions—one for high-grade samples and the other for low-grade samples. The restriction of high grade samples is because of their low priori likelihoods, whereas the main reason to restrict low grade samples is their nature as being waste and low analysis precisions. The two constraints reinforce each other in terms of enhancing the variables of estimates. A detailed case study on an epithermal gold deposit is carried out in terms of both cross validation and block modeling, showing that restricted kriging is superior over OK in mimicking the variables of original data. 相似文献
5.
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 相似文献
6.
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. 相似文献
7.
Estimation or simulation? That is the question 总被引：1，自引：0，他引：1
Jorge Kazuo Yamamoto 《Computational Geosciences》2008,12(4):573-591
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. 相似文献
8.
9.
10.
Compositional data are very common in the earth sciences. Nevertheless, little attention has been paid to the spatial interpolation of these data sets. Most interpolators do not necessarily satisfy the constant sum and nonnegativity constraints of compositional data, nor take spatial structure into account. Therefore, compositional kriging is introduced as a straightforward extension of ordinary kriging that complies with these constraints. In two case studies, the performance of compositional kriging is compared with that of the additive logratio-transform. In the first case study, compositional kriging yielded significantly more accurate predictions than the additive logratio-transform, while in the second case study the performances were comparable. 相似文献
11.
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. 相似文献
12.
Xavier Emery 《Mathematical Geology》2006,38(7):801-819
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. 相似文献
13.
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. 相似文献
14.
Jorge Kazuo Yamamoto 《Mathematical Geology》2000,32(4):489-509
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. 相似文献
15.
Correcting the Smoothing Effect of Ordinary Kriging Estimates 总被引：2，自引：0，他引：2
Jorge Kazuo Yamamoto 《Mathematical Geology》2005,37(1):69-94
The smoothing effect of ordinary kriging is a well-known dangerous effect associated with this estimation technique. Consequently kriging estimates do not reproduce both histogram and semivariogram model of sample data. A four-step procedure for correcting the smoothing effect of ordinary kriging estimates is shown to be efficient for the reproduction of histogram and semivariogram without loss of local accuracy. Furthermore, this procedure provides a unique map sharing both local and global accuracies. Ordinary kriging with a proper correction for smoothing effect can be revitalized as a reliable estimation method that allows a better use of the available information. 相似文献
16.
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. 相似文献
17.
Histograms of observations from spatial phenomena are often found to be more heavy-tailed than Gaussian distributions, which
makes the Gaussian random field model unsuited. A T-distributed random field model with heavy-tailed marginal probability density functions is defined. The model is a generalization
of the familiar Student-T distribution, and it may be given a Bayesian interpretation. The increased variability appears cross-realizations, contrary
to in-realizations, since all realizations are Gaussian-like with varying variance between realizations. The T-distributed random field model is analytically tractable and the conditional model is developed, which provides algorithms
for conditional simulation and prediction, so-called T-kriging. The model compares favourably with most previously defined random field models. The Gaussian random field model
appears as a special, limiting case of the T-distributed random field model. The model is particularly useful whenever multiple, sparsely sampled realizations of the
random field are available, and is clearly favourable to the Gaussian model in this case. The properties of the T-distributed random field model is demonstrated on well log observations from the Gullfaks field in the North Sea. The predictions
correspond to traditional kriging predictions, while the associated prediction variances are more representative, as they
are layer specific and include uncertainty caused by using variance estimates. 相似文献
18.
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 相似文献
19.
Geostatistical estimations of the hydraulic conductivity field (K) in the Carrizo aquifer, Texas, are performed over three regional domains of increasing extent: 1) the domain corresponding
to a three-dimensional groundwater flow model previously built (model domain); 2) the area corresponding to the 10 counties
encompassing the model domain (County domain), and; 3) the full extension of the Carrizo aquifer within Texas (Texas domain).
Two different approaches are used: 1) an indirect approach where transmissivity (T) is estimated first and K is retrieved through division of the T estimate by the screen length of the wells, and; 2) a direct approach where K data are kriged directly. Due to preferential well screen emplacement, and scarcity of sampling in the deeper portions of
the formation (> 1 km), the available data set is biased toward high values of hydraulic conductivities. Kriging combined
with linear regression, simple kriging with varying local means, kriging with an external drift, and cokriging allow the incorporation
of specific capacity as secondary information. Prediction performances (assessed through cross-validation) differ according
to the chosen approach, the considered variable (log-transformed or back-transformed), and the scale of interest. For the
indirect approach, kriging of log T with varying local means yields the best estimates for both log-transformed and back-transformed variables in the model domain.
For larger regional scales (County and Texas domains), cokriging performs generally better than other kriging procedures when
estimating both (log T)∗ and T∗. Among procedures using the direct approach, the best prediction performances are obtained using kriging of log K with an external drift. Overall, geostatistical estimation of the hydraulic conductivity field at regional scales is rendered
difficult by both preferential well location and preferential emplacement of well screens in the most productive portions
of the aquifer. Such bias creates unrealistic hydraulic conductivity values, in particular, in sparsely sampled areas. 相似文献
20.
When do we need a trend model in kriging? 总被引：1，自引：0，他引：1
Under usual estimation practice with local search windows for data and for interpolation situations, universal kriging and ordinary kriging yield the same estimates, using a data set with apparent trend, for both the unknown attribute and its trend component. Modeling the trend matters only in extrapolation situations. Because conditions of the case study presented arise most frequently in practice, the simpler ordinary kriging is the preferred option. 相似文献