共查询到20条相似文献，搜索用时 31 毫秒
1.
Kriging without negative weights 总被引：1，自引：0，他引：1
Under a constant drift, the linear kriging estimator is considered as a weighted average ofn available sample values. Kriging weights are determined such that the estimator is unbiased and optimal. To meet these requirements, negative kriging weights are sometimes found. Use of negative weights can produce negative block grades, which makes no practical sense. In some applications, all kriging weights may be required to be nonnegative. In this paper, a derivation of a set of nonlinear equations with the nonnegative constraint is presented. A numerical algorithm also is developed for the solution of the new set of kriging equations. 相似文献
2.
In mineral deposit evaluation, the estimation of block values can involve the solution of as many as 100,000 systems of linear equations. To minimize computing costs it is important that the program to solve the linear systems be as efficient as possible. Standard library routines tend to contain very general stable methods using pivoting and iteration to ensure high accuracy solutions. Here, a particular method has been developed that gives a substantial improvement in run time at the occasional expense of accuracy and stability. The method is not new—just a simple application of Gaussian elimination in the symmetric case. 相似文献
3.
Chris Roth 《Mathematical Geology》1998,30(8):999-1009
Lognormal kriging was developed early in geostatistics to take account of the often seen skewed distribution of the experimental mining data. Intuitively, taking the distribution of the data into account should lead to a better local estimate than that which would have been obtained when it is ignored. In practice however, the results obtained are sometimes disappointing. This paper tries to explain why this is so from the behavior of the lognormal kriging estimator. The estimator is shown to respect certain unbiasedness properties when considering the whole working field using the regression curve and its confidence interval for both simple or ordinary kriging. When examined locally, however, the estimator presents a behavior that is neither expected nor intuitive. These results lead to the question: is the theoretically correct lognormal kriging estimator suited to the practical problem of local estimation? 相似文献
4.
Interval-valued random functions and the kriging of intervals 总被引：1，自引：0，他引：1
Phil Diamond 《Mathematical Geology》1988,20(3):145-165
Estimation procedures using data that include some values known to lie within certain intervals are usually regarded as problems of constrained optimization. A different approach is used here. Intervals are treated as elements of a positive cone, obeying the arithmetic of interval analysis, and positive interval-valued random functions are discussed. A kriging formalism for interval-valued data is developed. It provides estimates that are themselves intervals. In this context, the condition that kriging weights be positive is seen to arise in a natural way. A numerical example is given, and the extension to universal kriging is sketched. 相似文献
5.
Transition probability-based indicator geostatistics 总被引：30，自引：0，他引：30
Traditionally, spatial continuity models for indicator variables are developed by empirical curvefitting to the sample indicator (cross-) variogram. However, geologic data may be too sparse to permit a purely empirical approach, particularly in application to the subsurface. Techniques for model synthesis that integrate hard data and conceptual models therefore are needed. Interpretability is crucial. Compared with the indicator (cross-) variogram or indicator (cross-) covariance, the transition probability is more interpretable. Information on proportion, mean length, and juxtapositioning directly relates to the transition probability: asymmetry can be considered. Furthermore, the transition probability elucidates order relation conditions and readily formulates the indicator (co)kriging equations. 相似文献
6.
Kriging Regionalized Positive Variables Revisited: Sample Space and Scale Considerations 总被引：1，自引：0，他引：1
Frequently, regionalized positive variables are treated by preliminarily applying a logarithm, and kriging estimates are back-transformed
using classical formulae for the expectation of a lognormal random variable. This practice has several problems (lack of robustness,
non-optimal confidence intervals, etc.), particularly when estimating block averages. Therefore, many practitioners take exponentials
of the kriging estimates, although the final estimations are deemed as non-optimal. Another approach arises when the nature
of the sample space and the scale of the data are considered. Since these concepts can be suitably captured by an Euclidean
space structure, we may define an optimal kriging estimator for positive variables, with all properties analogous to those
of linear geostatistical techniques, even for the estimation of block averages. In this particular case, no assumption on
preservation of lognormality is needed. From a practical point of view, the proposed method coincides with the median estimator
and offers theoretical ground to this extended practice. Thus, existing software and routines remain fully applicable. 相似文献
7.
Six different geostatistical estimators (linear kriging, lognormal kriging, and disjunctive kriging, each with and without a nonbias, i.e., universality condition) were compared using data from a polymetallic deposit in Algeria. The differences between estimators with and without the nonbias condition were far more pronounced than between the different kriging methods. This highlights the importance of choosing an appropriate stationarity model for the data. The criterion concerning kriging weight of the mean in simple kriging, proposed by Remacre (1984, 1987) and Rivoirard (1984) was found to be helpful for determining blocks where the choice of the stationarity hypothesis was critical. 相似文献
8.
G. J. McArthur 《Mathematical Geology》1988,20(4):343-366
The Hellyer orebody, a polymetallic massive sulfide deposit, was discovered in western Tasmania by Aberfoyle in 1983. Delineation diamond drilling was carried out in 1984 on a nominal 50-m square grid pattern to outline the resource. Resource estimation methods were influenced by the requirement to develop a regular block model for conceptual mine planning studies. Detailed geological interpretation indicates that the interpolation technique must take into account several important features to retain geological credibility. The deposit has sharp limits defined by visual geological contacts with virtually barren enclosing rocks. Lateral terminations are rapid with no interfingering internal waste. The dip and strike are variable and a major fault with a measurable displacement cuts acutely through the center of the deposit. Ore grades are reasonably correlateable within specific layers from hole to hole indicating a significant across-dip anisotropy. A hanging wall enriched zone is well-defined throughout the deposit. To overcome the variable geometry problems, a stratigraphic coordinate system was defined arbitrarily to replace the normal z coordinate. This allowed variography in stratigraphic layers. Blocks to be estimated were constrained by hand-drawn and subsequently digitized hanging wall and footwall contours. Each block was ascribed a stratigraphic coordinate by calculating its spatial position in relation to nearby stratigraphic unit boundaries within the massive sulfide body. Estimates were generated for each element by ordinary linear kriging. Despite the relatively sparse data in a large massive deposit, the customized technique developed for Hellyer has provided a reliable model of spatial grade distribution by combining conventional geostatistical methods with careful geological observation and interpretation. Some geometry problems remain which are the subject of ongoing studies.This paper was presented at MGUS 87 Conference, Redwood City, California, 14 April 1987. 相似文献
9.
C. Lajaunie 《Mathematical Geology》1990,22(1):123-144
Approximate local confidence intervals can be produced by nonlinear methods designed to estimate indicator variables. The most precise of these methods, the conditional expectation, can only be used in practice in the multi-Gaussian context. Theoretically, less efficient methods have to be used in more general cases. The methods considered here are indicator kriging, probability kriging (indicator-rank co-kriging), and disjunctive kriging (indicator co-kriging). The properties of these estimators are studied in this paper in the multi-Gaussian context, for this allows a more detailed study than under more general models. Conditional distribution approximation is first studied. Exact results are given for mean squared errors and conditional bias. Then conditional quantile estimators are compared empirically. Finally, confidence intervals are compared from the points of view of bias and precision. 相似文献
10.
J. A. Vargas-Guzmán 《Mathematical Geology》2004,36(3):307-322
This paper introduces geostatistical approaches (i.e., kriging estimation and simulation) for a group of non-Gaussian random fields that are power algebraic transformations of Gaussian and lognormal random fields. These are power random fields (PRFs) that allow the construction of stochastic polynomial series. They were derived from the exponential random field, which is expressed as Taylor series expansion with PRF terms. The equations developed from computation of moments for conditional random variables allow the correction of Gaussian kriging estimates for the non-Gaussian space. The introduced PRF geostatistics shall provide tools for integration of data that requires simple algebraic transformations, such as regression polynomials that are commonly encountered in the practical applications of estimation. The approach also allows for simulations drawn from skewed distributions. 相似文献
11.
Geostatistical regionalization of glacial aquitard thickness in northwestern Germany, based on fuzzy kriging 总被引：1，自引：0，他引：1
In many instances hydrogeological parameters obtained by conventional methods for selected localities within an aquifer or an aquitard are not sufficient for adequate regionalization at the scale of the entire layer. Here, we demonstrate an application of the fuzzy kriging method in regionalization of hydrogeological data, in which the set of conventional, crisp values is supplemented by imprecise information subjectively estimated by an expert. It is believed that such an approach eventually may reflect the real-world conditions more closely than a traditional crisp-value approach, because the former does not impose exactness artificially on phenomena which are diffuse by their nature. Spatial interpolation was done for the thickness of one of the major aquitards (till and glaciolacustrine clay) in northwestern Germany. The dataset consists of 329 crisp values from boreholes supplemented by 172 imprecise values defined as fuzzy numbers. It is demonstrated that the reliability of regionalization was higher, compared to regionalization performed with the crisp dataset only. Fuzzy kriging was performed with FUZZEKS (Fuzzy Evaluation and Kriging System) developed at the Ecosystem Research Center at the University of Kiel. 相似文献
12.
胡小荣 《物探化探计算技术》1998,20(3):249-255
泛克立格法是地质统计学的一种重要估值方法,当区域化变量在空间变异几何域内非平稳时常用泛克立格法来估值。然而,用泛克立格方程组计算估计权值时由于未对权的符号作任何限制,从而使得计算出的权值经常出现负权现象,而负权的存在有许多弊端,所以在许多应用中有必要对权值作非负要求。本文基于线性规划方法提出了一种考虑权值非负约束的泛克立格算法,该算法既考虑到了权值的非负约束条件,又利用了线性规划方法求解简便快捷的优点。 相似文献
13.
A geostatistically based approach is developed for the identification of aquifer transmissivities in Yolo Basin, California. The approach combines weighted least-squares with universal kriging and cokriging techniques in an overall scheme that (1)considers a prioriknown information on aquifer transmissivity and specific capacities of wells, (2)considers uncertainties in water level and transmissivity data, and (3)estimates the reliability of the generated transmissivity values. Minimization of a global least-squares function that incorporates calibration and plausibility criteria leads to a transmissivity map that shows a good agreement with pumping-test results. 相似文献
14.
For earthquake ground motion studies, the actual ground motion distribution should be reproduced as accurately as possible. For optimal estimation of ground motion, kriging has been shown to provide accurate estimates. Although kriging is accurate for this application, some estimates it provides are underestimates. This has dire consequences for subsequent design for earthquake resistance. Kriging does not provide enough information to allow an analysis of each estimate for underestimation. For such an application, disjunctive kriging is better applied. This advanced technique quantifies the probability that an estimate equals or exceeds particular levels of ground motion. Furthermore, disjunctive kriging can provide improved estimation accuracy when applied for local estimation of ground motion. 相似文献
15.
Efficient updating of kriging estimates and variances 总被引：2，自引：0，他引：2
This short note presents a method for efficiently updating ordinary kriging estimates and variances when one or more additional samples are incorporated into the kriging system. First, the foundation linear algebra result is presented. Then the update equations are derived. Finally, an illustrative application of updating is briefly discussed. 相似文献
16.
17.
Max D. Morris 《Mathematical Geology》1991,23(7):929-943
In planning spatial sampling studies for the purpose of estimating the semivariogram, the number of data pairs separated by a given distance is sometimes used as a comparative index of the precision which can be expected from a given sampling design. Because spatial data are correlated, this index can be unreliable. An alternative index which partially corrects for this correlation, themaximum equivalent uncorrelated pairs, is proposed for comparing spatial designs. The index is developed under the assumption that the underlying stochastic process is Gaussian and is appropriate when the (population) semivariogram is to be estimated by the sample semivariogram. 相似文献
18.
Adding Local Accuracy to Direct Sequential Simulation 总被引：3，自引：0，他引：3
Jef Caers 《Mathematical Geology》2000,32(7):815-850
Geostatistical simulations are globally accurate in the sense that they reproduce global statistics such as variograms and histograms. Kriging is locally accurate in the minimum local error variance sense. Building on the concept of direct sequential simulation, we propose a fast simulation method that can share these opposing objectives. It is shown that the multiple-point entropy of the resulting simulation is related to the univariate entropy of the local conditional distributions used to draw simulated values. Adding local accuracy to conditional simulations does not detract much from variogram reproduction and can be used to increase multiple-point entropy. The methods developed are illustrated using a case study. 相似文献
19.
纯时间域多元信息地质统计学 总被引：5，自引：0，他引：5
随着人们动态地认识客观世界的要求日益增加,时空域中的多元信息地质统计学技术研究受到重视。本文在时间域内对几种纯时间域克立格技术进行了研究,重点对纯时间域中的普通克立格（ＴＯＫ）、简单克立格（ＴＳＫ）及泛克立格（ＴＵＫ）进行了论述,也对纯时间域中的协同克立格（ＴＣＯＫ）和指示克立格（ＴＩＫ）作了简要的概括。 相似文献
20.
Spatial prediction is a problem common to many disciplines. A simple application is the mapping of an attribute recorded at a set of points. Frequently a nonlinear functional of the observed variable is of interest, and this calls for nonlinear approaches to prediction. Nonlinear kriging methods, developed in recent years, endeavour to do so and additionally provide estimates of the distribution of the target quantity conditional on the observations. There are few empirical studies that validate the various forms of nonlinear kriging. This study compares linear and nonlinear kriging methods with respect to precision and their success in modelling prediction uncertainty. The methods were applied to a data set giving measurements of the topsoil concentrations of cobalt and copper at more than 3000 locations in the Border Region of Scotland. The data stem from a survey undertaken to identify places where these trace elements are deficient for livestock. The comparison was carried out by dividing the data set into calibration and validation sets. No clear differences between the precision of ordinary, lognormal, disjunctive, indicator, and model-based kriging were found, neither for linear nor for nonlinear target quantities. Linear kriging, supplemented with the assumption of normally distributed prediction errors, failed to model the conditional distribution of the marginally skewed data, whereas the nonlinear methods modelled the conditional distributions almost equally well. In our study the plug-in methods did not fare any worse than model-based kriging, which takes parameter uncertainty into account. 相似文献