共查询到20条相似文献，搜索用时 31 毫秒
1.
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. 相似文献
2.
Hybrid Estimation of Semivariogram Parameters 总被引：1，自引：0，他引：1
Two widely used methods of semivariogram estimation are weighted least squares estimation and maximum likelihood estimation.
The former have certain computational advantages, whereas the latter are more statistically efficient. We introduce and study
a “hybrid” semivariogram estimation procedure that combines weighted least squares estimation of the range parameter with
maximum likelihood estimation of the sill (and nugget) assuming known range, in such a way that the sill-to-range ratio in
an exponential semivariogram is estimated consistently under an infill asymptotic regime. We show empirically that such a
procedure is nearly as efficient computationally, and more efficient statistically for some parameters, than weighted least
squares estimation of all of the semivariogram’s parameters. Furthermore, we demonstrate that standard plug-in (or empirical)
spatial predictors and prediction error variances, obtained by replacing the unknown semivariogram parameters with estimates
in expressions for the ordinary kriging predictor and kriging variance, respectively, perform better when hybrid estimates
are plugged in than when weighted least squares estimates are plugged in. In view of these results and the simplicity of computing
the hybrid estimates from weighted least squares estimates, we suggest that software that currently estimates the semivariogram
by weighted least squares methods be amended to include hybrid estimation as an option. 相似文献
3.
On unbiased backtransform of lognormal kriging estimates 总被引：4，自引：0，他引：4
Jorge Kazuo Yamamoto 《Computational Geosciences》2007,11(3):219-234
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. 相似文献
4.
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. 相似文献
5.
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. 相似文献
6.
James R. Carr 《Mathematical Geology》1990,22(8):1063-1079
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. 相似文献
7.
P. S. N. Murthy 《Mathematical Geology》1989,21(4):443-461
This paper presents the results of disjunctive kriging applied to a supergene iron ore deposit of Bailadila Range of India. Disjunctive kriging is applied firstly to compare estimates of the blocks by ordinary kriging and secondly to estimate benchwise local recoverable reserves of the orebody. Good agreement exists between block estimates by ordinary kriging and disjunctive kriging except for peripheral blocks with less borehole information. Estimation of benchwise reserves shows that the behavior of the distribution of grades is different in various benches. The study shows that disjunctive kriging can be applied successfully for estimation of local recoverable reserves in the case of a good grade hematite iron ore deposit. 相似文献
8.
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. 相似文献
9.
Pierre Goovaerts 《Mathematical Geosciences》2008,40(1):101-128
This paper presents a methodology to conduct geostatistical variography and interpolation on areal data measured over geographical units (or blocks) with different sizes and shapes, while accounting for heterogeneous weight or kernel functions within those units. The deconvolution method is iterative and seeks the point-support model that minimizes the difference between the theoretically regularized semivariogram model and the model fitted to areal data. This model is then used in area-to-point (ATP) kriging to map the spatial distribution of the attribute of interest within each geographical unit. The coherence constraint ensures that the weighted average of kriged estimates equals the areal datum.This approach is illustrated using health data (cancer rates aggregated at the county level) and population density surface as a kernel function. Simulations are conducted over two regions with contrasting county geographies: the state of Indiana and four states in the Western United States. In both regions, the deconvolution approach yields a point support semivariogram model that is reasonably close to the semivariogram of simulated point values. The use of this model in ATP kriging yields a more accurate prediction than a naïve point kriging of areal data that simply collapses each county into its geographic centroid. ATP kriging reduces the smoothing effect and is robust with respect to small differences in the point support semivariogram model. Important features of the point-support semivariogram, such as the nugget effect, can never be fully validated from areal data. The user may want to narrow down the set of solutions based on his knowledge of the phenomenon (e.g., set the nugget effect to zero). The approach presented avoids the visual bias associated with the interpretation of choropleth maps and should facilitate the analysis of relationships between variables measured over different spatial supports. 相似文献
10.
Goovaerts P 《Mathematical Geology》2008,40(1):101-128
This paper presents a methodology to conduct geostatistical variography and interpolation on areal data measured over geographical units (or blocks) with different sizes and shapes, while accounting for heterogeneous weight or kernel functions within those units. The deconvolution method is iterative and seeks the pointsupport model that minimizes the difference between the theoretically regularized semivariogram model and the model fitted to areal data. This model is then used in area-to-point (ATP) kriging to map the spatial distribution of the attribute of interest within each geographical unit. The coherence constraint ensures that the weighted average of kriged estimates equals the areal datum.This approach is illustrated using health data (cancer rates aggregated at the county level) and population density surface as a kernel function. Simulations are conducted over two regions with contrasting county geographies: the state of Indiana and four states in the Western United States. In both regions, the deconvolution approach yields a point support semivariogram model that is reasonably close to the semivariogram of simulated point values. The use of this model in ATP kriging yields a more accurate prediction than a na?ve point kriging of areal data that simply collapses each county into its geographic centroid. ATP kriging reduces the smoothing effect and is robust with respect to small differences in the point support semivariogram model. Important features of the point-support semivariogram, such as the nugget effect, can never be fully validated from areal data. The user may want to narrow down the set of solutions based on his knowledge of the phenomenon (e.g., set the nugget effect to zero). The approach presented avoids the visual bias associated with the interpretation of choropleth maps and should facilitate the analysis of relationships between variables measured over different spatial supports. 相似文献
11.
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. 相似文献
12.
Spatial analyses of groundwater levels using universal kriging 总被引：6，自引：0，他引：6
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. 相似文献
13.
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). 相似文献
14.
P. Goovaerts 《Mathematical Geology》1998,30(5):511-534
This paper presents both estimation and simulation as optimization problems that differ in the optimization criteria, minimization of a local expected loss for estimation and reproduction of global statistics (semivariogram, histogram) for simulation. An intermediate approach is proposed whereby an initial random image is gradually modified using simulated annealing so as to better match both local and global constraints. The relative weights of the different constraints in the objective function allow the user to strike a balance between smoothness of the estimated map and reproduction of spatial variability by simulated maps. The procedure is illustrated using a synthetic dataset. The proposed approach is shown to enhance the influence of observations on neighboring simulated values, hence the final realizations appear to be better conditioned to the sample information. It also produces maps that are more accurate (smaller prediction error) than stochastic simulation ignoring local constraints, but not as accurate as E-type estimation. Flow simulation results show that accounting for local constraints yields, on average, smaller errors in production forecast than a smooth estimated map or a simulated map that reproduces only the histogram and semivariogram. The approach thus reduces the risk associated with the use of a single realization for forecasting and planning. 相似文献
15.
Xavier Emery 《Mathematical Geology》2005,37(3):295-319
Multigaussian kriging is used in geostatistical applications to assess the recoverable reserves in ore deposits, or the probability for a contaminant to exceed a critical threshold. However, in general, the estimates have to be calculated by a numerical integration (Monte Carlo approach). In this paper, we propose analytical expressions to compute the multigaussian kriging estimator and its estimation variance, thanks to polynomial expansions. Three extensions are then considered, which are essential for mining and environmental applications: accounting for an unknown and locally varying mean (local stationarity), accounting for a block-support correction, and estimating spatial averages. All these extensions can be combined; they generalize several known techniques like ordinary lognormal kriging and uniform conditioning by a Gaussian value. An application of the concepts to a porphyry copper deposit shows that the proposed “ordinary multigaussian kriging” approach leads to more realistic estimates of the recoverable reserves than the conventional methods (disjunctive and simple multigaussian krigings), in particular in the nonmineralized undersampled areas. 相似文献
16.
A parametric study of robustness of kriging variance as a function of range and relative nugget effect for a spherical semivariogram 总被引：1，自引：0，他引：1
P. I. Brooker 《Mathematical Geology》1986,18(5):477-488
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. 相似文献
17.
Restricted kriging for mixture of grade models 总被引：2，自引：0，他引：2
A modified type of kriging, referred to as restricted kriging (RK), is proposed in this study. The method incorporates constraints on different grade classes to restrict the influence of the samples having different likelihoods in estimation. RK is motivated by the estimation of mineral reserves when grades have highly skewed distributions. Ordinary kriging tends to produce an overly smoothed interpolated surface by underestimating high grades and overestimating low grades. The fact that ordinary kriging gives a uniform prior treatment to all samples independent of their values is a major factor associated with this smoothing effect. The new approach differentiates each grade portion by preselected cutoffs. RK is developed for a single cutoff and then extended into a general form for any finite number of cutoffs. Restricted cokriging (RCK) is also formulated to simultaneously estimate a set of random functions with restriction conditions. Methods are suggested for determination of the probabilities of occurrence of different grade portions. Finally, the new approach is demonstrated on a case study of an epithermal gold deposit. 相似文献
18.
《Applied Geochemistry》2005,20(1):157-168
In monitoring a minor geochemical element in groundwater or soils, a background population of values below the instrumental detection limit is frequently present. When those values are found in the monitoring process, they are assigned to the detection limit which, in some cases, generates a probability mass in the probability density function of the variable at that value (the minimum value that can be detected). Such background values could distort both the estimation of the variable at nonsampled locations and the inference of the spatial structure of variability of the variable. Two important problems are the delineation of areas where the variable is above the detection limit and the estimation of the magnitude of the variables inside those areas. The importance of these issues in geochemical prospecting or in environmental sciences, in general related with contamination and environmental monitoring, is obvious. In this paper the authors describe the two-step procedure of indicator kriging and ordinary kriging and compare it with empirical maximum likelihood kriging. The first approach consists of using a binary indicator variable for estimating the probability of a location being above the detection limit, plus ordinary kriging conditional to the location being above the detection limit. An estimation variance, however, is not available for that estimator. Empirical maximum likelihood kriging, which was designed to deal with skew distributions, can also deal with an atom at the origin of the distribution. The method uses a Bayesian approach to kriging and gives intermittency in the form of a probability map, its estimates providing a realistic assessment of their estimation variance. The pros and cons of each method are discussed and illustrated using a large dataset of As concentration in groundwater. The results of the two methods are compared by cross-validation. 相似文献
19.
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 相似文献
20.
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. 相似文献