FuzzyKrig: a comprehensive matlab toolbox for geostatistical estimation of imprecise information

Using kriging has been accepted today as the most common method of estimating spatial data in such different fields as the geosciences. To be able to apply kriging methods, it is necessary that the data and variogram model parameters be precise. To utilize the imprecise (fuzzy) data and parameters, use is made of fuzzy kriging methods. Although it has been 30 years since different fuzzy kriging algorithms were proposed, its use has not become as common as other kriging methods (ordinary, simple, log, universal, etc.); lack of a comprehensive software that can perform, based on different fuzzy kriging algorithms, the related calculations in a 3D space can be the main reason. This paper describes an open-source software toolbox (developed in Matlab) for running different algorithms proposed for fuzzy kriging. It also presents, besides a short presentation of the fuzzy kriging method and introduction of the functions provided by the FuzzyKrig toolbox, 3 cases of the software application under the conditions where: 1) data are hard and variogram model parameters are fuzzy, 2) data are fuzzy and variogram model parameters are hard, and 3) both data and variogram model parameters are fuzzy.     

Fuzzy Modeling for Reserve Estimation Based on Spatial Variability

This article addresses a new reserve estimation method which uses fuzzy modeling algorithms and estimates the reserve parameters based on spatial variability. The proposed fuzzy modeling approach has three stages: (1) Structure identification and preliminary clustering, (2) Variogram analysis, and (3) Clustering based rule system. A new clustering index approach and a new spatial measure function (point semimadogram) are proposed in the paper. The developed methodology uses spatial variability in each step and takes the fuzzy rules from input-output data. The model has been tested using both simulated and real data sets. The performance evaluation indicates that the new methodology can be applied in reserve estimation and similar modeling problems     

Matheronian geostatistics—Quo vadis?

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.     

A parametric study of robustness of kriging variance as a function of range and relative nugget effect for a spherical semivariogram

In geostatistics, an estimation of blocks of a deposit is reported along with the variance of error made in their estimation. This calculation is based on the model chosen for the semivariogram of the deposit so that mistakes in its estimation can manifest themselves in the perception of accuracy with which blocks are known. Changes in kriging variance resulting from various amounts of error in modeling the relative nugget effect and range of the semivariogram are investigated for an extensive set of spherical semivariograms.     

Imprecise (fuzzy) information in geostatistics

A methodology based on fuzzy set theory for the utilization of imprecise data in geostatistics is presented. A common problem preventing a broader use of geostatistics has been the insufficient amount of accurate measurement data. In certain cases, additional but uncertain (soft) information is available and can be encoded as subjective probabilities, and then the soft kriging method can be applied (Journel, 1986). In other cases, a fuzzy encoding of soft information may be more realistic and simplify the numerical calculations. Imprecise (fuzzy) spatial information on the possible variogram is integrated into a single variogram which is used in a fuzzy kriging procedure. The overall uncertainty of prediction is represented by the estimation variance and the calculated membership function for each kriged point. The methodology is applied to the permeability prediction of a soil liner for hazardous waste containment. The available number of hard measurement data (20) was not enough for a classical geostatistical analysis. An additional 20 soft data made it possible to prepare kriged contour maps using the fuzzy geostatistical procedure.This paper was presented at MGUS 87 Conference, Redwood City, California, 14 April 1987.     

New variogram modeling method using MGGP and SVR

A critical step for kriging in geostatistics is estimation of the variogram. Traditional variogram modeling comprise of the experimental variogram calculation, appropriate variogram model selection and model parameter determination. Selecting of the variogram model and fitting of model parameters is the most controversial aspect of geostatistics. Shapes of valid variogram models are finite, and sometimes, the optimal shape of the model can not be fitted, leading to reduced estimation accuracy. In this paper, a new method is presented to automatically construct a model shape and fit model parameters to experimental variograms using Support Vector Regression (SVR) and Multi-Gene Genetic Programming (MGGP). The proposed method does not require the selection of a variogram model and can directly provide the model shape and parameters of the optimal variogram. The validity of the proposed method is demonstrated in a number of cases.     

3D矿床建模技术在数字矿产勘查中的应用

为弥补现有数字矿床建模技术在地质矿产勘查处理和应用中的不足,从原始勘查数据建库及标准化、多指标单工程矿体自动圈定、基于语义识别的剖面矿体的连接与外推、矿体表面和品位建模、基于剖面矿体线框模型构建矿体表面模型及基于TIN+Octree数据结构和地质统计学理论建立矿体的空间属性模型等5个方面总结提出一套面向地质矿产勘查业务处理的矿床建模流程和总体技术解决方案.提高了地质矿产勘查研究精度,为进一步的矿山开采提供可靠的数据模型.      

Conditional Simulation of Random Fields with Bivariate Gamma Isofactorial Distributions

This work focuses on a random function model with gamma marginal and bivariate isofactorial distributions, which has been applied in mining geostatistics for estimating recoverable reserves by disjunctive kriging. The objective is to widen its use to conditional simulation and further its application to the modeling of continuous attributes in geosciences. First, the main properties of the bivariate gamma isofactorial distributions are analyzed, with emphasis in the destructuring of the extreme values, the presence of a proportional effect (higher variability in high-valued areas), and the asymmetry in the spatial correlation of the indicator variables with respect to the median threshold. Then, we provide examples of stationary random functions with such bivariate distributions, for which the shape parameter of the marginal distribution is half an integer. These are defined as the sum of squared independent Gaussian random fields. An iterative algorithm based on the Gibbs sampler is proposed to perform the simulation conditional to a set of existing data. Such ‘multivariate chi-square’ model generalizes the well-known multigaussian model and is more flexible, since it allows defining a shape parameter which controls the asymmetry of the marginal and bivariate distributions.     

Stochastic Distance Based Geological Boundary Modeling with Curvilinear Features

Obtaining accurate geological boundaries and assessing the uncertainty in these limits are critical for effective ore resource and reserve estimation. The uncertainty in the extent of an ore body can be the largest source of uncertainty in ore resource estimation when drilling is sparse. These limits are traditionally interpreted deterministically and it can be difficult to quantify uncertainty in the boundary and its impact on ore tonnage. The proposed methodology is to consider stochastic modeling of the ore boundary with a distance function recoding of the available data. This technique is modified to incorporate non-stationarities in the form of a locally varying anisotropy field used in kriging and sequential Gaussian simulation. Implementing locally varying anisotropy kriging retains the geologically realistic features of a deterministic model while allowing for a stochastic assessment of uncertainty. A case study of a gold deposit in Northern Canada is used to demonstrate the methodology. The proposed technique generates realistic, curvilinear geological boundary models and allows for an assessment of the uncertainty in the model.     

Two Ordinary Kriging Approaches to Predicting Block Grade Distributions

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.     

Generalized sequential kriging (GSK) for heterogeneous cell size: property modeling without matrices

Sequential kriging avoids the use of matrices and resolves the issue of unstable solutions. It allows for stepwise ways to get joint estimations and cosimulations that are equivalent to the simultaneous solution. The approach is proposed as the solution for geocellular modeling with variable cell size from heterogeneous structural properties (HSPs) as required for modeling with structural constraints. Rock properties are controlled by structural domains, regions, and structural geology parameters. In some cases, rock properties are cross-correlated to formation thickness, curvature of structures, and other structural attributes. Cell thickness may be proportional to formation thickness and may enter as a conditioning property in the estimation of rock property parameters for simulation. In addition, cell volume controls the upscaling of covariance structures (i.e., regularized variograms). Structural properties are priorly modeled. Perturbation response functions (PRFs) are computed for each cell vs all possible sample point locations to facilitate sequential kriging. Upscaled PRFs are modified following conditional updating after each new data value is included in the estimation of parameters. Generalized sequential kriging is expected to become the main tool for real-time spatial modeling of 3D cellular models with HSP. In addition, some new developments related to the sequential kriging algorithm are included. Sequential kriging can be used for the estimation of parameters for simulation in the so-called unstructured grids.     

地质统计学方法在井间储层参数预测中的应用

本文探讨了利用地质统计学中的三种克立格方法进行井间储层参数预测的原理及方法,并在某油田中应用,通过最优性检验,建立起了适合研究区地质特点的各种储层参数的最佳估计模型,揭示了储层参数在三维空间的变化规律,为建立三维储层预测模型开辟了新的方法。     

某金矿技术经济评价及品位指标优化研究

根据某金矿床的地质特征,利用地质统计学方法对矿床的地质储量进行了评估,并给出了矿石量、金属量及品位的数学模型。在矿山现行生产能力的基础上,对矿山未来开发经济效益进行了２７５个方案的评价,并利用模糊数学理论对各方案进行了多目标优化决策 。     

基于3DMine软件的地质统计学在河北羊崖山铁矿床资源量估算中的应用

现有矿业软件的储量估算模块核心方法主要基于地质统计学。文章以河北迁安羊崖山铁矿床为例,通过3DMine软件在数据资料准备、变异函数确定、块体模型构建及储量报告等多方面的应用,规范了矿业软件在资源量估算方面的工作流程。利用普通克里金法和距离幂次反比法分别进行品位插值,完成资源量估算;并与传统储量估算结果进行比较,分析产生误差原因,建立了开采境界模型,计算了矿山开采储量和保有储量;利用赋值参数分析矿体控制程度,指导未来探矿工程重点部位。     

矿床数学—经济模型及其在元江金矿的应用

利用地质统计学的原理和方法，编制出矿床数学－经济模型软件包，它作为一种储量计算和矿床经济评价方法。介绍了该软件包各模块的功能，并以元江金矿应用为例，展示了该软件包的功能及效果。     

Simple and Ordinary Multigaussian Kriging for Estimating Recoverable Reserves

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.     

地质统计学中变差函数参数估计的新方法

遗传算法是一种模拟生物进化规律的全局优化算法。对传统的遗传算法进行改进,并应用于地质统计学变差函数参数估计中。实例分析表明,该方法简便、通用,具有较高拟合精度,是非线性、不连续可微模型参数估计的方法     

Successive Nonparametric Estimation of Conditional Distributions

Spatial characterization of non-Gaussian attributes in earth sciences and engineering commonly requires the estimation of their conditional distribution. The indicator and probability kriging approaches of current nonparametric geostatistics provide approximations for estimating conditional distributions. They do not, however, provide results similar to those in the cumbersome implementation of simultaneous cokriging of indicators. This paper presents a new formulation termed successive cokriging of indicators that avoids the classic simultaneous solution and related computational problems, while obtaining equivalent results to the impractical simultaneous solution of cokriging of indicators. A successive minimization of the estimation variance of probability estimates is performed, as additional data are successively included into the estimation process. In addition, the approach leads to an efficient nonparametric simulation algorithm for non-Gaussian random functions based on residual probabilities.     

Application of disjunctive kriging for estimating economic grade distribution in an iron ore deposit: A case study of the Choghart North Anomaly,Iran

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%.