共查询到20条相似文献,搜索用时 15 毫秒
1.
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. 相似文献
2.
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. 相似文献
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.
《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. 相似文献
5.
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). 相似文献
6.
A number of criteria based on kriging variance calculations may be used for infill sampling design in geologic site characterization. Searching for the best new sample locations from a set of candidate locations can result in excessive computation time if these criteria and the naive rekriging are used. The relative updated kriging estimate and variance for universal kriging estimation are demonstrated as a simple kriging estimate and variance, respectively. The updated kriging variance is demonstrated as the multiplication of two kriging variances. Using these updated kriging variance equations can increase the computational speed for selecting the best new sample locations. The application results for oil rock thickness in an oilfield indicate that minimizing the average relative updated kriging variance is a useful alternative to the other criteria based on kriging variance in optimal infill sampling design for geologic site characterization. 相似文献
7.
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 相似文献
8.
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. 相似文献
9.
A. G. Journel 《Mathematical Geology》1977,9(6):563-586
In the last few years, an increasing number of practical studies using so-called kriging estimation procedures have been published. Various terms, such as universal kriging, lognormal kriging, ordinary kriging, etc., are used to define different estimation procedures, leaving a certain confusion about what kriging really is. The object of this paper is to show what is the common backbone of all these estimation procedures, thus justifying the common name of kriging procedures. The word kriging (in French krigeage) is a concise and convenient term to designate the classical procedure of selecting, within agiven class of possible estimators, the estimator with a minimum estimation variance (i.e., the estimator which leads to a minimum variance of the resulting estimation error). This estimation variance can be seen as a squared distance between the unknown value and its estimator; the process of minimization of this distance can then be seen as the projection of the unknown value onto the space within which the search for an estimator is carried out. 相似文献
10.
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. 相似文献
11.
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. 相似文献
12.
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. 相似文献
13.
Chunfa Wu Jiaping Wu Yongming Luo Haibo Zhang Ying Teng Stephen D. DeGloria 《Environmental Earth Sciences》2011,63(5):1093-1103
It was not unusual in soil and environmental studies that the distribution of data is severely skewed with several high peak
values, which causes the difficulty for Kriging with data transformation to make a satisfied prediction. This paper tested
an approach that integrates kriging and triangular irregular network interpolation to make predictions. A data set consisting
of total Copper (Cu) concentrations of 147 soil samples, with a skewness of 4.64 and several high peak values, from a copper
smelting contaminated site in Zhejiang Province, China. The original data were partitioned into two parts. One represented
the holistic spatial variability, followed by lognormal distribution, and then was interpolated by lognormal ordinary kriging.
The other assumed to show the local variability of the area that near to high peak values, and triangular irregular network
interpolation was applied. These two predictions were integrated into one map. This map was assessed by comparing with rank-order
ordinary kriging and normal score ordinary kriging using another data set consisting of 54 soil samples of Cu in the same
region. According to the mean error and root mean square error, the approach integrating lognormal ordinary kriging and triangular
irregular network interpolation could make improved predictions over rank-order ordinary kriging and normal score ordinary
kriging for the severely skewed data with several high peak values. 相似文献
14.
基于不同地质统计方法的渗透系数场对污染物运移的影响 总被引:1,自引:0,他引:1
渗透系数场的空间变异性是影响污染物运移结果的决定因素,而地质统计方法是解决渗透系数空间变异性的主要技术手段。本文利用野外场地实测数据,采用普通克里格法和指示克里格法、顺序高斯模拟法和顺序指示模拟法四种地质统计方法,插值估测和模拟再现随机渗透系数场,进而对比研究四种渗透系数场对大尺度污染物运移的影响。研究结果表明,污染羽的质心位置(一阶矩)主要由渗透系数的平均值来决定;污染羽在空间上的展布范围(二阶矩)主要受渗透系数空间变异方差的影响;条件模拟克服了估计法的平滑效果,较好地再现真实曲线的波动性,渗透系数( lnK)估计方差与污染羽空间二阶矩随着条件模拟次数的增加而减小,并且顺序指示模拟程度更加明显。 相似文献
15.
Pierre Petitgas Mathieu Woillez Mathieu Doray Jacques Rivoirard 《Mathematical Geosciences》2018,50(2):187-208
Marine research survey data on fish stocks often show a small proportion of very high-density values, as for many environmental data. This makes the estimation of second-order statistics, such as the variance and the variogram, non-robust. The high fish density values are generated by fish aggregative behaviour, which may vary greatly at small scale in time and space. The high values are thus imprecisely known, both in their spatial occurrence and order of magnitude. To map such data, three indicator-based geostatistical methods were considered, the top-cut model, min–max autocorrelation factors (MAF) of indicators, and multiple indicator kriging. In the top-cut and MAF approaches, the variable is decomposed into components and the most continuous ones (those corresponding to the low and medium values) are used to guide the mapping. The methods are proposed as alternatives to ordinary kriging when the variogram is difficult to estimate. The methods are detailed and applied on a spatial data set of anchovy densities derived from a typical fish stock acoustic survey performed in the Bay of Biscay, which show a few high-density values distributed in small spatial patches and also as solitary events. The model performances are analyzed by cross-validating the data and comparing the kriged maps. Results are compared to ordinary kriging as a base case. The top-cut model had the best cross-validation performance. The indicator-based models allowed mapping high-value areas with small spatial extent, in contrast to ordinary kriging. Practical guidelines for implementing the indicator-based methods are provided. 相似文献
16.
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 相似文献
17.
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. 相似文献
18.
《Applied Geochemistry》1999,14(1):133-145
Three univariate geostatistical methods of estimation are applied to a geochemical data set. The studied methods are: ordinary kriging (cross-validation), factorial kriging, and indicator kriging. These techniques use the probabilistic and spatial behaviour of geochemical variables, giving a tool for identifying potential anomalous areas to locate mineralization. Ordinary kriging is easy to apply and to interpret the results. It has the advantage of using the same experimental grid points for its estimates, and no additional grid points are needed. Factorial kriging decomposes the raw variable into as many components as there are identified structures in the variogram. This, however, is a complex method and its application is more difficult than that of ordinary or indicator kriging. The main advantages of indicator kriging are that data are used by their rank order, being more robust about outlier values, and that the presentation of results is simple. Nevertheless, indicator kriging is incapable of separating anomalous values and the high values from the background, which have a behaviour different to the anomaly. In this work, the results of the application of these 3 kriging methods to a set of mineral exploration data obtained from a geochemical survey carried out in NW Spain are presented. This area is characterised by the presence of Au mineral occurrences. The kriging methods were applied to As, considered as a pathfinder of Au in this area. Numerical treatment of Au is not applicable, because it presents most values equal to the detection limit, and a series of extreme values. The results of the application of ordinary kriging, factorial kriging and indicator kriging to As make possible the location of a series of rich values, sited along a N–S shear zone, considered a structure related to the presence of Au. 相似文献
19.
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. 相似文献
20.
Goovaerts P 《Mathematical Geosciences》2010,42(5):535-554
A common issue in spatial interpolation is the combination of data measured over different spatial supports. For example,
information available for mapping disease risk typically includes point data (e.g. patients’ and controls’ residence) and
aggregated data (e.g. socio-demographic and economic attributes recorded at the census track level). Similarly, soil measurements
at discrete locations in the field are often supplemented with choropleth maps (e.g. soil or geological maps) that model the
spatial distribution of soil attributes as the juxtaposition of polygons (areas) with constant values. This paper presents
a general formulation of kriging that allows the combination of both point and areal data through the use of area-to-area,
area-to-point, and point-to-point covariances in the kriging system. The procedure is illustrated using two data sets: (1)
geological map and heavy metal concentrations recorded in the topsoil of the Swiss Jura, and (2) incidence rates of late-stage
breast cancer diagnosis per census tract and location of patient residences for three counties in Michigan. In the second
case, the kriging system includes an error variance term derived according to the binomial distribution to account for varying
degree of reliability of incidence rates depending on the total number of cases recorded in those tracts. Except under the
binomial kriging framework, area-and-point (AAP) kriging ensures the coherence of the prediction so that the average of interpolated
values within each mapping unit is equal to the original areal datum. The relationships between binomial kriging, Poisson
kriging, and indicator kriging are discussed under different scenarios for the population size and spatial support. Sensitivity
analysis demonstrates the smaller smoothing and greater prediction accuracy of the new procedure over ordinary and traditional
residual kriging based on the assumption that the local mean is constant within each mapping unit. 相似文献