RE: [ai-geostats] Adding kriging variances.
- Hi Digby,
The variance of a sum of random variables is not equal to the sum of their variances, except if they are independent.. Consequently, the kriging variances cannot be simply combined to compute the variance of the global estimator, see Mining Geostatistics, Page 323 (Journel and Huijbregts, 1978).
The most flexible approach is to use stochastic simulation to generate a set of realizations of the block values, aggregate them, and use the empirical distribution of aggregated block values as a model of uncertainty.
You can find an example in my publication http://home.comcast.net/~goovaerts/karen.pdf and there is a whole
book devoted to the use of stochastic simulation in mining.
Journel and Kyriakidis. 2004. Evaluation of Mineral reserves. A simulation approach. Oxford University Press.
From: Digby Millikan [mailto:digbym@...]
Sent: Mon 11/7/2005 1:09 AM
To: AI Geostats mailing list
Subject: [ai-geostats] Adding kriging variances.
If you had a set of blocks is it possible to add their kriging variances
to get a standard error for the mine head grade, is this something that is
in practice or has been done in the past to estimate possible deviations
kriged head grade?
- Dear Digby!If your samples are Normal, your variance is chi-squared!!!You can combine (not add together) kriging variances for separate blocks but, unless they are widely separated or big enough to be considered independent this is a non-trivial task. Basically you need to store how much weight is given to each sample for each block, geometry of all blocks and compute a humungeous estimation variance.Generally for head grade purposes, stopes or blocks are widely separated enough to combine the variances. Weight by size and shape of block (or tonnage) as a proportion of the total -- squared.Isobel
Digby Millikan <digbym@...> wrote:* By using the ai-geostats mailing list you agree to follow its rules
If you had a set of blocks is it possible to add their kriging variances together,
to get a standard error for the mine head grade, is this something that is done
in practice or has been done in the past to estimate possible deviations from a
kriged head grade?
( see http://www.ai-geostats.org/help_ai-geostats.htm )
* To unsubscribe to ai-geostats, send the following in the subject or in the body (plain text format) of an email message to sympa@...