共查询到20条相似文献，搜索用时 62 毫秒
1.
Indicator principal component kriging 总被引：1，自引：0，他引：1
An alternative to multiple indicator kriging is proposed which approximates the full coindicator kriging system by kriging the principal components of the original indicator variables. This transformation is studied in detail for the biGaussian model. It is shown that the cross-correlations between principal components are either insignificant or exactly zero. This result allows derivation of the conditional cumulative density function (cdf) by kriging principal components and then applying a linear back transform. A performance comparison based on a real data set (Walker Lake) is presented which suggests that the proposed method achieves approximation of the conditional cdf equivalent to indicator cokriging but with substantially less variogram modeling effort and at smaller computational cost. 相似文献
2.
James R. Carr 《Mathematical Geology》1994,26(5):605-621
Consideration of order relations is key to indicator kriging, indicator cokriging, and probability kriging, especially for the latter two methods wherein the additional modeling of cross-covariance contributes to an increased chance of violating order relations. Herein, Gaussian-type curves are fit to estimates of the cumulative distribution function (cdf) at data quantiles to: (1) yield smoothed estimates of the cdf; and (2) to correct for violations of order relations (i.e., to correct for situations wherein the estimate of the cdf for a larger quantile is less than that for a smaller quantile). Smoothed estimates of the cdf are sought as a means to improve the approximation to the integral equation for the expected value of the regionalized variable in probability kriging. Experiments show that this smoothing yields slightly improved estimation of the expected value (in probability kriging). Another experiment, one that uses the same variogram for all indicator functions, does not yield improved estimates.Presented at the 25th Anniversary Meeting of the IAMG, Prague, Czech Republic, October 10–15, 1993. 相似文献
3.
Probability kriging is implemented in a general cokriging procedure (c.f. Myers, 1982) for estimatingboth the indicator and uniform transforms. Paired-sum semi-variograms are used to facilitate the modeling of the cross-covariance between the uniform transform and each indicator transform. Estimates of the uniform transform are averaged over all cutoffs, the average used to derive an estimate of the original data. This estimate can be biased with respect to the mean data value, but is unbiased with respect to the data median. 相似文献
4.
On the Equivalence of the Cokriging and Kriging Systems 总被引：2，自引：0，他引：2
Simple cokriging of components of a p-dimensional second-order stationary random process is considered. Necessary and sufficient conditions under which simple cokriging is equivalent to simple kriging are given. Essentially this condition requires that it should be possible to express the cross-covariance at any lag series h using the cross-covariance at |h|=0 and the auto-covariance at lag series h. The mosaic model, multicolocated kriging and the linear model of coregionalization are examined in this context. A data analytic method to examine whether simple kriging of components of a multivariate random process is equivalent to its cokriging is given 相似文献
5.
Indicator cokriging is an alternative to disjunctive kriging for estimation of spatial distributions. One way to determine which of these techniques is more accurate for estimation of spatial distributions is to apply each to a particular type of data. A procedure is developed for evaluation of disjunctive kriging and indicator cokriging for such an application. Application of this procedure to earthquake ground motion data found disjunctive kriging to be at least as accurate as indicator cokriging for estimation of spatial distributions of peak horizontal acceleration. Indicator cokriging was superior for all other types of earthquake ground motion data. 相似文献
6.
The problem of estimating a regionalized variable in the presence of other secondary variables is encountered in spatial investigations. Given a context in which the secondary variable is known everywhere (or can be estimated with great precision), different estimation methods are compared: regression, regression with residual simple kriging, kriging, simple kriging with a mean obtained by regression, kriging with an external drift, and cokriging. The study focuses on 19 pairs of regionalized variables from five different datasets representing different domains (geochemical, environmental, geotechnical). The methods are compared by cross-validation using the mean absolute error as criterion. For correlations between the principal and secondary variable under 0.4, similar results are obtained using kriging and cokriging, and these methods are superior slightly to the other approaches in terms of minimizing estimation error. For correlations greater than 0.4, cokriging generally performs better than other methods, with a reduction in mean absolute errors that can reach 46% when there is a high degree of correlation between the variables. Kriging with an external drift or kriging the residuals of a regression (SKR) are almost as precise as cokriging. 相似文献
7.
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. 相似文献
8.
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. 相似文献
9.
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. 相似文献
10.
Environmental studies require multivariate data such as chemical concentrations with space-time coordinates. There are two
general conditions related to such data: the existence of correlations among the coregionalized variables and the differences
in numbers of data which occur because of insufficient data caused by measurement error or bad weather conditions. This study
proposes geostatistical techniques for space-time multivariate modeling that take into consideration these correlations and
data absences. These techniques consist of suitable modeling of semivariograms and cross-semivariograms for quantifying correlation
structures among multivariables and of extending standardized ordinary cokriging. The tensor product cubic smoothing surface
method is used for space-time semivariogram modeling. These methods are applied to the chemical component data of the Ariake
Sea, a typical closed sea in southwest Japan. In order to clarify environmental changes in the Ariake Sea, the concentration
data of four nutritive salts (NO2–N, NO3–N, NH4–N, and PO4–P) at 38 stations over 25 years are used as environmental indicators. For each of the kinds of data, there are spaces and
times for which there is no data available. The effectiveness of the modeling of space-time semivariograms and the high estimation
capability of the extended cokriging are demonstrated by cross-validation. Compared with ordinary kriging for a single variable,
multivariate space-time standardized ordinary cokriging can provide a more detailed concentration map of nutritive salts and
while elucidating their temporal changes over sparsely spaced data areas. In the space-time models by ordinary kriging, on
the other hand, smooth trends are obvious. 相似文献
11.
Under the intrinsic coregionalization model if both primary and secondary measurements are available at all sample locations, the conventional geostatistical wisdom is that cokriging provides exactly the same solution as univariate kriging on the primary process alone. However, recent eamples have been given where nonzero secondary cokriging weights have accurred under this spatial dependence structure. This note identifies the conditions under which secondary information is useful under the assumption of intrinsic coregionalization. An illustration is given using a dataset of plutonium and americium concentrations collected from a region of the Nevada Test Site. 相似文献
12.
Three approaches for estimating the hydraulic conductivity (K) of the Trifa aquifer, Morocco were investigated: (1) kriging of the K values obtained from pumping tests, (2) cokriging of the pumping test data with electrical resistivity data as a secondary variable, and (3) cokriging of the pumping test data with the slope of the water table. Gauss-transformed values of the variables are used because they provide more robust variograms and transformed values of the primary and secondary variables show correlations higher than the raw values, which is beneficial in cokriging. In cokriging with electrical resistivity, two zones are considered since the geological deposits are different from the north to the south of the aquifer, which is reflected in different correlations between the variables. Comparison of the three approaches is based mainly on the estimation errors, and to a lesser degree on the cross-validations of the corresponding variogram models and general considerations, like the measurements’ reliability and aquifer make-up. The best-estimated K is given by cokriging with the slope of the water table and is therefore preferred for further use in groundwater flow modeling. Thus, electrical resistivity or the slope of the water table can both be used as secondary variables to estimate K, especially in heterogeneous aquifers with lateral variations in lithology, as is the case of the Trifa aquifer. 相似文献
13.
地质统计学是数学地质领域最为活跃而实用的分支,它是以区域化变量理论为基础,以变异函数为基本工具,研究那些在空间分布上既具有随机性又具有结构性的自然现象的科学。在第四纪研究中的很多特征(变量)均可看成区域化变量进行地质统计学分析。作者在讨论了经典概率论及数理统计方法简单地应用于第四纪研究可能出现的问题后,着重介绍了用于第四纪研究中的若干地质统计学方法及基本理论,同时,对地质统计学方法应用于第四纪研究中的前景进行了分析。 相似文献
14.
Covariance models provide the basic measure of spatial continuity in geostatistics. Traditionally, a closed-form analytical model is fitted to allow for interpolation of sample Covariance values while ensuring the positive definiteness condition. For cokriging, the modeling task is made even more difficult because of the restriction imposed by the linear coregionalization model. Bochner's theorem maps the positive definite constraints into much simpler constraints on the Fourier transform of the covariance, that is the density spectrum. Accordingly, we propose to transform the experimental (cross) covariance tables into quasidensity spectrum tables using Fast Fourier Transform (FFT). These quasidensity spectrum tables are then smoothed under constraints of positivity and unit sum. A backtransform (FFT) yields permissible (jointly) positive definite (cross) covariance tables. At no point is any analytical modeling called for and the algorithm is not restricted by the linear coregionalization model. A case study shows the proposed covariance modeling to be easier and much faster than the traditional analytical covariance modeling, yet yields comparable kriging or simulation results. 相似文献
15.
Ordinary Cokriging Revisited 总被引：12，自引：0，他引：12
P. Goovaerts 《Mathematical Geology》1998,30(1):21-42
This paper sets up the relations between simple cokriging and ordinary cokriging with one or several unbiasedness constraints. Differences between cokriging variants are related to differences between models adopted for the means of primary and secondary variables. Because it is not necessary for the secondary data weights to sum to zero, ordinary cokriging with a single unbiasedness constraint gives a larger weight to the secondary information while reducing the occurrence of negative weights. Also the weights provided by such cokriging systems written in terms of covariances or correlograms are not related linearly, hence the estimates are different. The prediction performances of cokriging estimators are assessed using an environmental dataset that includes concentrations of five heavy metals at 359 locations. Analysis of reestimation scores at 100 test locations shows that kriging and cokriging perform equally when the primary and secondary variables are sampled at the same locations. When the secondary information is available at the estimated location, one gains little by retaining other distant secondary data in the estimation. 相似文献
16.
比较岩性模型建立方法。首先,在高分辨率层序地层学的指导下,最大限度地应用地质、露头、三维地震、测井等静态资料,发挥井点资料垂向分辨率高,地震资料横向信息丰富的优势,在地质规律约束下建立不同时间的高精度等时地层格架模型。然后,在精细格架模型的基础上,以测井解释得到的岩相数据作为条件数据,分别采用指示克里格、截断高斯模拟、Object-modeling算法、贯指示模拟建立砂体展布模型。最后,通过抽稀检验评价不同算法对模拟结果的影响,实现算法及其参数的优选,从而指导整个区块不同开发阶段,不同井网密度时全区三维精细地质模型的建立,也可为具有相似地质环境的油田建立三维地质模型提供参考。通过比较,优选出指示克里格、序贯指示模拟两种算法都能较好表征本研究区地质情况。 相似文献
17.
Modeling Spatial Variability with One and Multidimensional Continuous-Lag Markov Chains 总被引：10，自引：0，他引：10
The continuous-lag Markov chain provides a conceptually simple, mathematically compact, and theoretically powerful model of spatial variability for categorical variables. Markov chains have a long-standing record of applicability to one-dimensional (1-D) geologic data, but 2- and 3-D applications are rare. Theoretically, a multidimensional Markov chain may assume that 1-D Markov chains characterize spatial variability in all directions. Given that a 1-D continuous Markov chain can be described concisely by a transition rate matrix, this paper develops 3-D continuous-lag Markov chain models by interpolating transition rate matrices established for three principal directions, say strike, dip, and vertical. The transition rate matrix for each principal direction can be developed directly from data or indirectly by conceptual approaches. Application of Sylvester's theorem facilitates establishment of the transition rate matrix, as well as calculation of transition probabilities. The resulting 3-D continuous-lag Markov chain models then can be applied to geo-statistical estimation and simulation techniques, such as indicator cokriging, disjunctive kriging, sequential indicator simulation, and simulated annealing. 相似文献
18.
Giuseppe Raspa Massimiliano Moscatelli Francesco Stigliano Antonio Patera Fabrizio Marconi Daiane Folle Roberto Vallone Marco Mancini Gian Paolo Cavinato Salvatore Milli Joo Felipe Coimbra Leite Costa 《Engineering Geology》2008,101(3-4):251-268
We are presenting an attempt to evaluate the spatial variability of geotechnical parameters in the upper Pleistocene–Holocene alluvial deposits of Roma (Italy) by means of multivariate geostatistics.The upper Pleistocene–Holocene alluvial deposits of Roma are sensitive to high levels of geohazard. They occupy a sizable and significant part of the city, being the foundation for many monuments, historical neighborhoods, and archaeological areas, and the main host of the present and future subway lines. We have stored information from more than 2000 geotechnical boreholes crossing the alluvial deposits into a relational database. For the present study, only the boreholes with lithologic/textural interpretation and geotechnical information were selected. The set includes 283 boreholes and 719 samples, which have a set of geotechnical information comprising physical properties and mechanical parameters.Techniques of multivariate statistics and geostatistics were combined and compared to evaluate the estimation methods of the mechanical parameters, with special reference to the drained friction angle from direct shear test (φ′). Principal Component Analysis was applied to the dataset to highlight the relationships between the geotechnical parameters. Through cross-validation analysis, multiple linear regression, kriging, and cokriging were tested as estimators of φ′. Cross-validation demonstrates that the cokriging with granulometries as auxiliary variables is the most suitable method to estimate φ′. In addition to proving that cokriging is a good estimator of φ′, cross-validation demonstrates that input data are coherent and this allows us to use them for estimation of geotechnical parameters, although they come from different laboratories and different vintages.Nevertheless, to get the same good results of cross-validation in estimation, it is necessary for granulometries to be available at grid points. Since this information being not available at all grid points, it is expected that, in the future, textural information can be derived in an indirect way, i.e., from lithologic/textural spatial reconstructions. 相似文献
19.
Large cokriging systems arise in many situations and are difficult to handle in practice. Simplifications such as simple kriging,
strictly collocated and multicollocated cokriging are often used and models under which such simplifications are, in fact,
equivalent to cokriging have recently received attention. In this paper, a two-dimensional second-order stationary random
process with known mean is considered and the redundancy of certain components of the data at certain locations vis-à-vis
the solution to the simple cokriging system is examined. Conditions for the simple cokriging weights of these components at
these locations are set to zero. The conditions generalise the notion of the autokrigeability coefficient and can, in principle,
be applied to any data configuration. In specific sampling situations such as the isotopic and certain heterotropic configurations,
models under which simple kriging, strictly collocated, multicollocated and dislocated cokriging are equivalent to simple
cokriging are readily identified and results already available in the literature are obtained. These are readily identified
and the results are already available in the literature. The advantage of the approach presented here is that it can be applied
to any data configuration for analysis of permissible simplifications in simple cokriging. 相似文献