共查询到20条相似文献，搜索用时 78 毫秒
1.
The classical variogram estimator proposed by Matheron can be written as a quadratic form of the observations. When data have an elliptically contoured distribution with constant mean, the correlation between the classical variogram estimator at two different lags is a function of the spatial design matrix, the covariance matrix, and the kurtosis. Several specific cases are studied closely. A subclass of elliptically contoured distributions with a particular family of covariance matrices is shown to possess exactly the same correlation structure for the classical variogram estimator as the multivariate independent Gaussian distribution. The consequences on variogram fitting by generalized least squares are discussed. 相似文献
2.
Marc G. Genton 《Mathematical Geology》2000,32(1):127-137
The classical variogram estimator proposed by Matheron can be written as a quadratic form of the observations. When data have an elliptically contoured distribution with constant mean, the correlation between the classical variogram estimator at two different lags is a function of the spatial design matrix, the covariance matrix, and the kurtosis. Several specific cases are studied closely. A subclass of elliptically contoured distributions with a particular family of covariance matrices is shown to possess exactly the same correlation structure for the classical variogram estimator as the multivariate independent Gaussian distribution. The consequences on variogram fitting by generalized least squares are discussed. 相似文献
3.
Four variogram models for regional groundwater geochemical data are presented. These models were developed from an empirical study of the sample variograms for more than 10 elements in groundwaters from two geologic regions in the Plainview quandrangle, Texas. A procedure is given for the estimation of the variogram in the isotropic and anisotropic case. The variograms were found useful for quantifying the differences in spatial variability for elements within a geologic unit and for elements in different geologic units. Additionally, the variogram analysis enables assessment of the assumption of statistical independence of regional samples which is commonly used in many statistical procedures. The estimated variograms are used in computation of kriged estimates for the Plainview quadrangle data. The results indicate that an inverse distance weighting model was superior for prediction than simple kriging with the particular variograms used. 相似文献
4.
Positive definiteness is not enough 总被引：2，自引：0，他引：2
M. Armstrong 《Mathematical Geology》1992,24(1):135-143
Geostatisticians know that the mathematical functions chosen to represent spatial covariances and variograms must have the appropriate type of positive definiteness, but they may not realize that there are restrictions on the types of covariances and variograms that are compatible with particular distributions. This paper gives some examples showing that (1) the spherical model is not compatible with the multivariate lognormal distribution if the coefficient of variation is 2.0 or more (even in 1-D), and (2) the Gaussian covariance and several other models are not compatible with indicator random functions. As these examples concern quite different types of random functions, it is clear that there is a general problem of compatibility between spatial covariance models (or variograms) and a specified multivariate distribution. The problem arises with all distributions except the multivariate normal, and not just the two cited here. The need for a general theorem giving the necessary and sufficient conditions for a covariance or a variogram to be compatible with a particular distribution is stressed. 相似文献
5.
Johannes Bruining Diederik van Batenburg Larry W. Lake An Ping Yang 《Mathematical Geology》1997,29(6):823-848
Random field generators serve as a tool to model heterogeneous media for applications in hydrocarbon recovery and groundwater
flow. Random fields with a power-law variogram structure, also termed fractional Brownian motion (fBm) fields, are of interest
to study scale dependent heterogeneity effects on one-phase and two-phase flow. We show that such fields generated by the
spectral method and the Inverse Fast Fourier Transform (IFFT) have an incorrect variogram structure and variance. To illustrate
this we derive the prefactor of the fBm spectral density function, which is required to generate the fBm fields. We propose
a new method to generate fBm fields that introduces weighting functions into the spectral method. It leads to a flexible and
efficient algorithm. The flexibility permits an optimal choice of summation points (that is points in frequency space at which
the weighting function is calculated) specific for the autocovariance structure of the field. As an illustration of the method,
comparisons between estimated and expected statistics of fields with an exponential variogram and of fBm fields are presented.
For power-law semivariograms, the proposed spectral method with a cylindrical distribution of the summation points gives optimal
results. 相似文献
6.
Yu-Pin Lin Hone-Jay Chu Yu-Long Huang Bai-You Cheng Tsun-Kuo Chang 《Environmental Earth Sciences》2011,62(2):299-311
This study proposes the method of simulating spatial patterns and quantifying the uncertainty in multivariate distribution
of heavy metals (Cr, Cu, Ni, and Zn) by sequential indicator simulation (SIS) combined with conditional Latin hypercube sampling
(cLHS) in Changhua County, Taiwan. The cLHS is used for a sampling then for SIS mapping and assessing uncertainties of heavy
metal concentrations. The indicator variogram results indicate that the 700 cLHS samples replicate statistical multivariate
distribution and spatial structure of the 1,082 samples. Moreover, the SIS realizations based on 700 cLHS samples are more
conservative and reliable than those based on 1,082 samples for delineating soil contamination by all heavy metals with the
exception of Zn. Given adequate sampling, soil contamination simulation provides sufficient information for delineating contaminated
areas and planning environmental management. 相似文献
7.
Earth scientists and land managers often wish to group sampling sites that are both similar with respect to their properties and near to one another on the ground. This paper outlines the geostatistical rationale for such spatial grouping and describes a multivariate procedure to implement it. Sample variograms are calculated from the original data or their leading principal components and then the parameters of the underlying functions are estimated. A dissimilarity matrix is computed for all sampling sites, preferably using Gower's general similarity coefficient. Dissimilarities are then modified using the variogram to incorporate the form and extent of spatial variation. A nonhierarchical classification of sampling sites is performed on the leading latent vectors of the modified dissimilarity matrix by dynamic clustering to an optimum. The technique is illustrated with results of its application to soil survey data from two small areas in Britain and from a transect. In the case of the latter results of spatially weighted classifications are compared with those of strict segmentation. An appendix lists a Genstat program for a spatially constrained classification using a spherical variogram as an example. 相似文献
8.
Principal component analysis (PCA) is commonly applied without looking at the spatial support (size and shape, of the samples and the field), and the cross-covariance structure of the explored attributes. This paper shows that PCA can depend on such spatial features. If the spatial random functions for attributes correspond to largely dissimilar variograms and cross-variograms, the scale effect will increase as well. On the other hand, under conditions of proportional shape of the variograms and cross-variograms (i.e., intrinsic coregionalization), no scale effect may occur. The theoretical analysis leads to eigenvalue and eigenvector functions of the size of the domain and sample supports. We termed this analysis growing scale PCA, where spatial (or time) scale refers to the size and shape of the domain and samples. An example of silt, sand, and clay attributes for a second-order stationary vector random function shows the correlation matrix asymptotically approaches constants at two or three times the largest range of the spherical variogram used in the nested model. This is contrary to the common belief that the correlation structure between attributes become constant at the range value. Results of growing scale PCA illustrate the rotation of the orthogonal space of the eigenvectors as the size of the domain grows. PCA results are strongly controlled by the multivariate matrix variogram model. This approach is useful for exploratory data analysis of spatially autocorrelated vector random functions. 相似文献
9.
Non-stationary Geostatistical Modeling Based on Distance Weighted Statistics and Distributions 总被引：3，自引：0，他引：3
A common assumption in geostatistics is that the underlying joint distribution of possible values of a geological attribute at different locations is stationary within a homogeneous domain. This joint distribution is commonly modeled as multi-Gaussian, with correlations defined by a stationary covariance function. This results in attribute maps that fail to reproduce local changes in the mean, in the variance and, particularly, in the spatial continuity. The proposed alternative is to build local distributions, variograms, and correlograms. These are inferred by weighting the samples depending on their distance to selected locations. The local distributions are locally transformed into Gaussian distributions embedding information on the local histogram. The distance weighted experimental variograms and correlograms are able to adapt to local changes in the direction and range of spatial continuity. The automatically fitted local variogram models and the local Gaussian transformation parameters are used in spatial estimation algorithms assuming local stationarity. The resulting maps are rich in nonstationary spatial features. The proposed process implies a higher computational effort than traditional stationary techniques, but if data availability allows for a reliable inference of the local distributions and statistics, a higher accuracy of estimates can be achieved. 相似文献
10.
The variogram is a critical input to geostatistical studies: (1) it is a tool to investigate and quantify the spatial variability of the phenomenon under study, and (2) most geostatistical estimation or simulation algorithms require an analytical variogram model, which they will reproduce with statistical fluctuations. In the construction of numerical models, the variogram reflects some of our understanding of the geometry and continuity of the variable, and can have a very important impact on predictions from such numerical models. The principles of variogram modeling are developed and illustrated with a number of practical examples. A three-dimensional interpretation of the variogram is necessary to fully describe geologic continuity. Directional continuity must be described simultaneously to be consistent with principles of geological deposition and for a legitimate measure of spatial variability for geostatistical modeling algorithms. Interpretation principles are discussed in detail. Variograms are modeled with particular functions for reasons of mathematical consistency. Used correctly, such variogram models account for the experimental data, geological interpretation, and analogue information. The steps in this essential data integration exercise are described in detail through the introduction of a rigorous methodology. 相似文献
11.
Teacher''s Aide Variogram Interpretation and Modeling 总被引：13，自引：0，他引：13
The variogram is a critical input to geostatistical studies: (1) it is a tool to investigate and quantify the spatial variability of the phenomenon under study, and (2) most geostatistical estimation or simulation algorithms require an analytical variogram model, which they will reproduce with statistical fluctuations. In the construction of numerical models, the variogram reflects some of our understanding of the geometry and continuity of the variable, and can have a very important impact on predictions from such numerical models. The principles of variogram modeling are developed and illustrated with a number of practical examples. A three-dimensional interpretation of the variogram is necessary to fully describe geologic continuity. Directional continuity must be described simultaneously to be consistent with principles of geological deposition and for a legitimate measure of spatial variability for geostatistical modeling algorithms. Interpretation principles are discussed in detail. Variograms are modeled with particular functions for reasons of mathematical consistency. Used correctly, such variogram models account for the experimental data, geological interpretation, and analogue information. The steps in this essential data integration exercise are described in detail through the introduction of a rigorous methodology. 相似文献
12.
Xiao-Lin Sun Yun-Jin Wu Hui-Li Wang Yu-Guo Zhao Gan-Lin Zhang 《Mathematical Geosciences》2014,46(4):429-443
Information on the spatial distribution of soil particle-size fractions (psf) is required for a wide range of applications. Geostatistics is often used to map spatial distribution from point observations; however, for compositional data such as soil psf, conventional multivariate geostatistics are not optimal. Several solutions have been proposed, including compositional kriging and transformation to a composition followed by cokriging. These have been shown to perform differently in different situations, so that there is no procedure to choose an optimal method. To address this, two case studies of soil psf mapping were carried out using compositional kriging, log-ratio cokriging, cokriging, and additive log-ratio cokriging; and the performance of Mahalanobis distance as a criterion for choosing an optimal mapping method was tested. All methods generated very similar results. However, the compositional kriging and cokriging results were slightly more similar to each other than to the other pair, as were log-ratio cokriging and additive log-ratio cokriging. The similar results of the two methods within each pair were due to similarities of the methods themselves, for example, the same variogram models and prediction techniques, and the similar results between the two pairs were due to the mathematical relationship between original and log-ratio transformed data. Mahalanobis distance did not prove to be a good indicator for selecting an optimal method to map soil psf. 相似文献
13.
14.
Variogram Fitting by Generalized Least Squares Using an Explicit Formula for the Covariance Structure 总被引：4，自引：0，他引：4
Marc G. Genton 《Mathematical Geology》1998,30(4):323-345
In the context of spatial statistics, the classical variogram estimator proposed by Matheron can be written as a quadratic form of the observations. If data are Gaussian with constant mean, then the correlation between the classical variogram estimator at two different lags is a function of the spatial design matrix and the variance matrix. When data are independent with unidimensional and regular support, an explicit formula for this correlation is available. The same is true for a multidimensional and regular support as can be shown by using Kronecker products of matrices. As variogram fitting is a crucial stage for correct spatial prediction, it is proposed to use a generalized least squares method with an explicit formula for the covariance structure (GLSE). A good approximation of the covariance structure is achieved by taking account of the explicit formula for the correlation in the independent situation. Simulations are carried out with several types of underlying variograms, as well as with outliers in the data. Results show that this technique (GLSE), combined with a robust estimator of the variogram, improves the fit significantly. 相似文献
15.
On the Use of Non-Euclidean Distance Measures in Geostatistics 总被引：4，自引：0，他引：4
Frank C. Curriero 《Mathematical Geology》2006,38(8):907-926
In many scientific disciplines, straight line, Euclidean distances may not accurately describe proximity relationships among
spatial data. However, non-Euclidean distance measures must be used with caution in geostatistical applications. A simple
example is provided to demonstrate there are no guarantees that existing covariance and variogram functions remain valid (i.e.
positive definite or conditionally negative definite) when used with a non-Euclidean distance measure. There are certain distance
measures that when used with existing covariance and variogram functions remain valid, an issue that is explored. The concept
of isometric embedding is introduced and linked to the concepts of positive and conditionally negative definiteness to demonstrate
classes of valid norm dependent isotropic covariance and variogram functions, results many of which have yet to appear in
the mainstream geostatistical literature or application. These classes of functions extend the well known classes by adding
a parameter to define the distance norm. In practice, this distance parameter can be set a priori to represent, for example,
the Euclidean distance, or kept as a parameter to allow the data to choose the metric. A simulated application of the latter
is provided for demonstration. Simulation results are also presented comparing kriged predictions based on Euclidean distance
to those based on using a water metric. 相似文献
16.
Guillaume Larocque Pierre Dutilleul Bernard Pelletier James W. Fyles 《Mathematical Geology》2007,39(3):263-288
Coregionalization analysis has been presented as a method of multi-scale analysis for multivariate spatial data. Despite an
increasing use of this method in environmental and earth sciences, the uncertainty associated with the estimation of parameters
in coregionalization analysis (e.g., sills and functions of sills) is potentially high and has not yet been characterized.
This article aims to discuss the theory underlying coregionalization analysis and assess the robustness and limits of the
method. A theoretical framework is developed to calculate the ergodic and fluctuation variance-covariance matrices of least-squares
estimators of sills in the linear model of coregionalization. To adjust for the positive semidefiniteness constraint on estimated
coregionalization matrices, a confidence interval estimation procedure for sills and functions of sills is presented. Thereafter,
the relative importance of uncertainty measures (bias and variance) for sills and structural coefficients of correlation and
determination is assessed under different scenarios to identify factors controlling their uncertainty. Our results show that
the sampling grid density, the choice of the least-squares estimator of sills, the positive semidefiniteness constraint, the
presence of scale dependence in the correlations, and the number and range of variogram models, all affect the level of uncertainty,
sometimes through multiple interactions. The asymptotic properties of variogram model parameter estimators in a bounded sampling
domain impose a theoretical limit to their accuracy and precision. Because of this limit, the uncertainty was found to be
high for several scenarios, especially with three variogram models, and was often more dependent on the ratio of variogram
range to domain extent than on the sampling grid density. In practice, in the coregionalization analysis of a real dataset,
the circular requirement for sill estimates in the calculation of uncertainty measures makes the quantification of uncertainty
very problematic, if not impossible. The use of coregionalization analysis must be made with due knowledge of the uncertainty
levels and limits of the method. 相似文献
17.
Factorial Kriging (FK) is a data- dependent spatial filtering method that can be used to remove both independent and correlated
noise on geological images as well as to enhance lineaments for subsequent geological interpretation. The spatial variability
of signal, noise, and lineaments, characterized by a variogram model, have been used explicitly in calculating FK filter coefficients
that are equivalent to the kriging weighting coefficients. This is in contrast to the conventional spatial filtering method
by predefined, data-independent filters, such as Gaussian and Sobel filters. The geostatistically optimal FK filter coefficients,
however, do not guarantee an optimal filtering effect, if filter geometry (size and shape) are not properly selected. The
selection of filter geometry has been investigated by examining the sensitivity of the FK filter coefficients to changes in
filter size as well as variogram characteristics, such as nugget effect, type, range of influence, and anisotropy. The efficiency
of data-dependent FK filtering relative to data-independent spatial filters has been evaluated through simulated stochastic
images by two examples. In the first example, both FK and data-independent filters are used to remove white noise in simulated
images. FK filtering results in a less blurring effect than the data-independent fillers, even for a filter size as large
as 9 × 9. In the second example, FK and data-independent filters are compared relative to the extraction of lineaments and
components showing anisotropic variability. It was determined that square windows of the filter mask are effective only for
removing Isotropie components or white noise. A nonsquare windows must be used if anisotropic components are to be filtered
out. FK filtering for lineament enhancement is shown to be resistant to image noise, whereas data-independent filters are
sensitive to the presence of noise. We also have applied the FK filtering to the GLORIA side-scan sonar image from the Gulf
of Mexico, illustrating that FK is superior to the data-independent filters in removing noise and enhancing lineaments. The
case study also demonstrate that variogram analysis and FK filtering can be used for large images if a spectral analysis and
optimal filter design in the frequency domain is prohibitive because of a large memory requirement. 相似文献
18.
Nonstationarity of the Mean and Unbiased Variogram Estimation: Extension of the Weighted Least-Squares Method 总被引：1，自引：0，他引：1
When concerned with spatial data, it is not unusual to observe a nonstationarity of the mean. This nonstationarity may be modeled through linear models and the fitting of variograms or covariance functions performed on residuals. Although it usually is accepted by authors that a bias is present if residuals are used, its importance is rarely assessed. In this paper, an expression of the variogram and the covariance function is developed to determine the expected bias. It is shown that the magnitude of the bias depends on the sampling configuration, the importance of the dependence between observations, the number of parameters used to model the mean, and the number of data. The applications of the expression are twofold. The first one is to evaluate a priori the importance of the bias which is expected when a residuals-based variogram model is used for a given configuration and a hypothetical data dependence. The second one is to extend the weighted least-squares method to fit the variogram and to obtain an unbiased estimate of the variogram. Two case studies show that the bias can be negligible or larger than 20%. The residual-based sample variogram underestimates the total variance of the process but the nugget variance may be overestimated. 相似文献
19.
Robust Resampling Confidence Intervals for Empirical Variograms 总被引：1，自引：0，他引：1
The variogram function is an important measure of the spatial dependencies of a geostatistical or other spatial dataset. It
plays a central role in kriging, designing spatial studies, and in understanding the spatial properties of geological and
environmental phenomena. It is therefore important to understand the variability attached to estimates of the variogram. Existing
methods for constructing confidence intervals around the empirical variogram either rely on strong assumptions, such as normality
or known variogram function, or are based on resampling blocks and subject to edge effect biases. This paper proposes two
new procedures for addressing these concerns: a quasi-block-bootstrap and a quasi-block-jackknife. The new methods are based
on transforming the data to decorrelate it based on a fitted variogram model, resampling blocks from the decorrelated data,
and then recorrelating. The coverage properties of the new confidence intervals are compared by simulation to a number of
existing resampling-based intervals. The proposed quasi-block-jackknife confidence interval is found to have the best properties
of all of the methods considered across a range of scenarios, including normally and lognormally distributed data and misspecification
of the variogram function used to decorrelate the data. 相似文献
20.
Geochemical characterization of heavy metal contaminated area using multivariate factorial kriging 总被引：1，自引：0，他引：1
Joaquim C. B. Queiroz José R. Sturaro Augusto C. F. Saraiva Paulo M. Barbosa Landim 《Environmental Geology》2008,55(1):95-105
This paper describes a geostatistical method, known as factorial kriging analysis, which is well suited for analyzing multivariate
spatial information. The method involves multivariate variogram modeling, principal component analysis, and cokriging. It
uses several separate correlation structures, each corresponding to a specific spatial scale, and yields a set of regionalized
factors summarizing the main features of the data for each spatial scale. This method is applied to an area of high manganese-ore
mining activity in Amapá State, North Brazil. Two scales of spatial variation (0.33 and 2.0 km) are identified and interpreted.
The results indicate that, for the short-range structure, manganese, arsenic, iron, and cadmium are associated with human
activities due to the mining work, while for the long-range structure, the high aluminum, selenium, copper, and lead concentrations,
seem to be related to the natural environment. At each scale, the correlation structure is analyzed, and regionalized factors
are estimated by cokriging and then mapped. 相似文献