共查询到20条相似文献，搜索用时 30 毫秒
1.
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. 相似文献
2.
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. 相似文献
3.
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. 相似文献
4.
Ali Akbar Daya 《Journal of the Geological Society of India》2014,83(5):567-576
Most significant iron ore deposits in Iran are located in Central Iran Zone. These deposits belong to the Bafq mining district. The Bafq mining district is located in the Early Cambrian Kashmar-Kerman volcanic arc of Central Iran. Linear estimation of regionalized variables (for example by inverse distance weighting or ordinary Kriging) results in relatively high estimation variances, i.e. the estimates have very low precision. Assessment of project economics (or other critical decision making) based on linear estimation is therefore risky. Non-linear estimation methods like disjunctive kriging perform better and the lower estimation variance allows less risky economic decision-making. Another advantage of disjunctive kriging is that it allows estimation of functions of the primary variable, which here is the grade (Fe %) of the ore. In particular it permits estimation of indicator functions defined using thresholds on the primary variable. This paper is devoted to application of disjunctive kriging method in Choghart North Anomaly iron ore deposit in Central Iran, Yazd province, Iran. In this study, the Fe concentration of Choghart North Anomaly iron ore deposit was modelled and estimated. The exploration data consists of borehole samples measuring the Fe concentration. A Gaussian isofactorial model is fitted to these data and disjunctive kriging was used to estimate the regionalized variable (Fe %) at unsampled locations and to assess the probabilities that the actual concentrations exceed a threshold value at a given location. Consequently a three dimensional model of probability of exceeding a threshold value and the estimated value are provided by disjunctive kriging to divide the ore into an economic and uneconomic part on the basis of estimation of indicator functions using thresholds grades defined on point support. The tools and concepts are complemented by a set of computer programs that are applied to the case study. The study showed that disjunctive kriging can be applied successfully for modeling the grade of an ore deposit. Results showed that the correlation between the estimated value and real value at locations close to each other is 81.9%. 相似文献
5.
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. 相似文献
6.
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. 相似文献
7.
Evaluation of spatial interpolation methods for groundwater level in an arid inland oasis,northwest China 总被引：3，自引：1，他引：2
Liqiang Yao Zailin Huo Shaoyuan Feng Xiaomin Mao Shaozhong Kang Jin Chen Jijun Xu Tammo S. Steenhuis 《Environmental Earth Sciences》2014,71(4):1911-1924
The Wuwei oasis, situated in the upper reaches of the Shiyang River basin in the arid inland of northwest China, is intensively cultivated using both groundwater and irrigation water originating from the Qilian Mountains. Groundwater levels are declining due to overuse of irrigation water. To estimate the decline over the entire Wuwei oasis, eight different interpolation methods were used for interpolating groundwater levels over 3 years, i.e. starting in 1983, followed by 1988 and ending with 1992. Cross-validation and orthogonal-validation were applied to evaluate the accuracy of the different methods. Root mean squared error and the correlation coefficient (R 2) were calculated for each of the interpolation methods and years. Three kriging methods (simply, ordinary, and universal) gave the best fit. Modified ordinary kriging was found better than simple and universal kriging methods with a smaller number of points having large differences (>50 m) between estimated and predicted values. Based on the groundwater surfaces determined by the ordinary kriging as modified by Yamamoto, the groundwater decline was found from 1983 to 1992 to be a modest 2.1 m in average. 相似文献
8.
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. 相似文献
9.
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 相似文献
10.
《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. 相似文献
11.
V. Kameshwara Rao C. R. Rao A. C. Narayana 《Journal of the Geological Society of India》2014,83(5):549-554
Conventional methods of ore deposit estimates are time consuming, whereas geostatistical methods provide quick and reliable estimates with minimum variance. Geostatistical tools, semi-variograms and kriging, have been used for estimation of grades of an iron ore deposit in the present study. In order to model the deposit and estimate grade, 4537 samples collected from 93 boreholes were used in the study. 3-D data have been converted to 2-D for analyzing the variation of Fe within the boreholes. For each borehole, the weighted mean of Fe grade and its coefficient of variation (CV) are calculated and further analysis is carried out for these two variables. Semi-variogram model suggests that the deposit extends over a zone of influence up to 700 m. Grade maps of kriged estimates reveal that the iron ore deposit is distributed in three distinct zones. 相似文献
12.
An Experimental Comparison of Ordinary and Universal Kriging and Inverse Distance Weighting 总被引：17，自引：0，他引：17
A factorial, computational experiment was conducted to compare the spatial interpolation accuracy of ordinary and universal kriging and two types of inverse squared-distance weighting. The experiment considered, in addition to these four interpolation methods, the effects of four data and sampling characteristics: surface type, sampling pattern, noise level, and strength of small-scale spatial correlation. Interpolation accuracy was measured by the natural logarithm of the mean squared interpolation error. Main effects of all five factors, all two-factor interactions, and several three-factor interactions were highly statistically significant. Among numerous findings, the most striking was that the two kriging methods were substantially superior to the inverse distance weighting methods over all levels of surface type, sampling pattern, noise, and correlation. 相似文献
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.
Geostatistics provides a suite of methods, summarized as kriging, to analyze a finite data set to describe a continuous property of the Earth. Kriging methods consist of moving window optimum estimation techniques, which are based on a least-squares principle and use a spatial structure function, usually the variogram. Applications of kriging techniques have become increasingly wide-spread, with ordinary kriging and universal kriging being the most popular ones. The dependence of the final map or model on the input, however, is not generally understood. Herein we demonstrate how changes in the kriging parameters and the neighborhood search affect the cartographic result. Principles are illustrated through a glaciological study. The objective is to map ice thickness and subglacial topography of Storglaciären, Kebnekaise Massif, northern Sweden, from several sets of radio-echo soundings and hot water drillings. New maps are presented. 相似文献
15.
Geostatistics provides a suite of methods, summarized as kriging, to analyze a finite data set to describe a continuous property of the Earth. Kriging methods consist of moving window optimum estimation techniques, which are based on a least-squares principle and use a spatial structure function, usually the variogram. Applications of kriging techniques have become increasingly wide-spread, with ordinary kriging and universal kriging being the most popular ones. The dependence of the final map or model on the input, however, is not generally understood. Herein we demonstrate how changes in the kriging parameters and the neighborhood search affect the cartographic result. Principles are illustrated through a glaciological study. The objective is to map ice thickness and subglacial topography of Storglaciären, Kebnekaise Massif, northern Sweden, from several sets of radio-echo soundings and hot water drillings. New maps are presented. 相似文献
16.
砂岩型铀矿是绿色经济可采的重要能源矿种,目前是世界上重要的铀矿勘探类型之一。二连盆地芒来铀矿床、鄂尔多斯盆地纳岭沟铀矿床铀矿体多呈板状产出,板状铀矿体成因备受关注。为研究砂岩型铀矿内部U、Ra和铀镭平衡系数（Kp）分布的垂向分布规律,本文以二连盆地芒来铀矿床为例,采用具有网格精度高优点的克里格插值法研究U、Ra和Kp垂向分布规律,分别利用放射性样品分析数据和定量伽马测井五点式反褶积法反演U含量数据,通过克里格法精细刻画砂岩型矿体中U、Ra和Kp分布的垂向分布规律。研究发现,该矿床板状矿体内部具有卷状的特征,Kp分布形态可以用来判断含氧含铀水的运移方向以及氧化强弱。该方法对研究铀矿体形态、铀矿成矿规律和后续地浸开采具有重要意义。 相似文献
17.
以浅剖数据为源数据,钻孔实测数据为验证数据,利用普通克里金法对海底地层厚度进行空间插值得到地层分布特征,采用3种半变异函数模型和不同取样间距对某井场3组地层厚度进行普通克里金插值并验证其插值效果。结果表明:普通克里金是一种有效的海底地层厚度预测方法;结构分析最佳的模型不一定是误差最小的模型,应对不同模型下的插值结果进行综合分析来选择最合适的模型,并提出球状模型在该井场厚度估计中最优,高斯模型次之;对于球状模型,增大取样间距对地层厚度变化剧烈的地层回归效果影响较小,对地层厚度变化不大的地层回归效果影响较大;同时,SE预测值变化率分析表明对于地层厚度变化剧烈的地层,减小取样间距可以大幅度地减少插值误差,而对于地层厚度变化不大的地层,减小取样间距对插值精度提高的意义不大。 相似文献
18.
This paper presents a modified ordinary kriging technique referred to as the “Area Influence Kriging” (AIK). The method is a simple and practical tool to use for more accurate prediction of global recoverable ore resources in any type of deposit. AIK performs well even in deposits with skewed grade distributions when the ordinary kriging (OK) results are unreasonably smooth. It is robust and globally unbiased like OK. The AIK method is not intended to replace OK, which is a better estimator of the average grade of the blocks. Rather it aims to complement OK with its excellent performance in predicting recoverable resources that have been the major pitfalls of OK in many resource estimation cases. The paper details the methodology of AIK with a couple of examples. It also reports the results from its application to a gold deposit. 相似文献
19.
Components of geostatistical estimation, developed as a method for ore deposit assessment, are discussed in detail. The assumption that spatial observations can be treated as a stochastic process is judged to be an inappropriate model for natural data. Problems of semivariogram formulation are reviewed, and this method is considered to be inadequate for estimating the function being sought. Characteristics of bivariate interpolation are summarized, highlighting kriging limitations as an interpolation method. Limitations are similar to those of inverse distance weighted observations interpolation. Attention is drawn to the local bias of kriging and misplaced claims that it is an “optimal” interpolation method. The so-called “estimation variance,” interpreted as providing confidence limits for estimation of mining blocks, is shown to be meaningless as an index of local variation. The claim that geostatistics constitutes a “new science” is examined in detail. Such novelties as exist in the method are shown to transgress accepted principles of scientific inference. Stochastic modeling in general is discussed, and purposes of the approach emphasized. For the purpose of detailed quantitative assessment it can provide only prediction qualified by hypothesis at best. Such an approach should play no part in ore deposit assessment where the need is for local detailed inventories; these can only be achieved properly through local deterministic methods, where prediction is purely deductive. 相似文献
20.
Marzeihe Shademan Khakestar Hassan Madani Hossein Hassani Parvizz Moarefvand 《Journal of the Geological Society of India》2013,81(4):581-585
Ordinary kriging and non-linear geostatistical estimators are now well accepted methods in mining grade control and mine reserve estimation. In kriging, the search volume or ‘kriging neighbourhood’ is defined by the user. The definition of the search space can have a significant impact on the outcome of the kriging estimate. In particular, too restrictive neighbourhood, can result in serious conditional bias. Kriging is commonly described as a ‘minimum variance estimator’ but this is only true when the neighbourhood is properly selected. Arbitrary decisions about search space are highly risky. The criteria to consider when evaluating a particular kriging neighbourhood are the slope of the regression of the ‘true’ and ‘estimated’ block grades, the number of kriging negative weights and the kriging variance. Search radius is one of the most important parameters of search volume which often is determined on the basis of influence of the variogram. In this paper the above-mentioned parameters are used to determine optimal search radius. 相似文献