next up previous contents
Next: Estimation of Resources Up: Regularization and Estimation Variance Previous: Regularization and Estimation Variance   Contents

Regularization by Cores along a Bore-hole

This type of regularization, very frequent in practice, corresponds to the construction of the variogram of the mean grades of core samples along the length of a bore-hole, cf. Figure 5.8. It is assumed that all the core samples have the same length $l$ and the same cross-sectional area $s$.

Figure 5.8: Core Samples Aligned along a Bore-hole.
/begin{figure}/begin{center}
/mbox
{/beginpicture
/setcoordinatesystem units <1c...
...0.6 0.7 1.25 /
/put {/small$s$} at 0.7 1.4
/endpicture}
/end{center}/end{figure}

. Let $Z({/boma x})$ be a point random function which is either second-order stationary or at least intrinsic with a semi-variogram $/gamma({/boma h})$.

The random function regularized on the support $v = s /times l$ of the core sample is written as

/begin{displaymath}Z_{v}(/boma{x})=/frac{1}{v} /int_{v(/boma{x})} Z(/boma{y})dy/ / / ,/end{displaymath}

where the sign $/int_{v({/boma x})}$ represents, in fact, a triple integral over the volume $v$. When the diameter of the core is small compared to length $l$, the regularization effect of the cross-sectional area $s$ of the core sample can be neglected. The mean value over the length $/ell$ of the core sample can then be written in terms of a single integral in the direction of the length of the core as

/begin{displaymath}Z_{v}(/boma{x}) /cong Z_/ell (/boma{x})= /frac{1}{/ell} /int_{/ell(/boma{x})}
Z(y)dy/ / / ./end{displaymath}

When the cross-sectional area $s$ of the core sample is negligible, two core samples can be considered as two aligned segments $/ell$ and $/ell_h$ of the same length $/ell$ and separated by a distance $h$, cf. Figure 5.8.

The regularized semi-variogram can then be written

/begin{displaymath}/gamma_{/ell}(/boma{h})=/frac{1}{2}E[Z_/ell (/boma{x}+/boma{h...
...{/gamma}(/ell,/ell_{/boma{h}}) - /bar{/gamma}(/ell,/ell)/ / / ./end{displaymath}

This regularization by cores assimilated to a segment of length $/ell$ will be illustrated on the simplest theoretical point variogram model, the linear model

/begin{displaymath}/gamma(/boma{h})= /; /vert /boma{h} /vert /; = r/ / / ./end{displaymath}

If $u$ denotes the coordinate along the length of the bore-hole, the first term $/bar{/gamma}(/ell,/ell_{/boma{h}})$ for the regularized semi-variogram is written

/begin{displaymath}/bar{/gamma}(/ell,/ell_{/boma{h}})=/frac{1}{/ell^2} /int^/ell_o
/int^{r+/ell}_{r} /vert
u-u' /vert du'du/ / / ./end{displaymath}

(Recall here that vector $/boma h$ is parallel to the direction of the bore-hole.)

According to whether or not $r = /vert h/vert$ is less than $/ell$, i.e., whether or not the two segments $/ell$ and $/ell_h$ overlap, we have the following two expressions (see Figure 5.9):

Figure 5.9: Geometrical Allocations of the Core Samples.
/begin{figure}/begin{center}
/mbox
{/beginpicture
/setcoordinatesystem units <1c...
... at 1.8 -2.25
/put {$u$} [l] at 5.1 0.0
}
/endpicture}
/end{center}/end{figure}

(i)
$/vert /boma{h} /vert /; /leq /; /ell:$

/begin{displaymath}/ell^{2}/bar{/gamma}(/ell,/ell_{/boma{h}})= /int_{o}^{r}
/int...
.../int_r^u (u-u')du'+ /int^{r+/ell}_{u} (u'-u)du'/right]du/ / / ./end{displaymath}

These integrals can be evaluated quite simply to give

/begin{displaymath}/bar{/gamma}(/ell,/ell_{/boma{h}})=/frac{1}{3} /frac{r^2}{/ell^2}
(3/ell - r) +
/frac{/ell}{3}/end{displaymath}

and, when $/vert /boma{h}/vert/ = 0 :
/bar{/gamma}(/ell,/ell)=/frac{/ell}{3}.$
(ii)
$/vert /boma{h} /vert /; /geq /; /ell:$

/begin{displaymath}/ell^2 /bar{/gamma}(/ell,/ell_{/boma{h}})=/int_{o}^/ell /int^{r+/ell}_{r} (u'-u)du'du =
r/ell^{2}/ / / ,/end{displaymath}

and, thus, $/bar{/gamma}(/ell,/ell_{/boma{h}})=r /; /mbox{for} /; /vert /boma{h} /vert /geq
/ell.$

Collecting the formulas, the regularized semi-variogram is

/begin{displaymath}/gamma_/ell (/boma{h}) = /left/{ /begin{array}{ll}
/frac{r^2}...
...rall r /geq /ell/ / / / / / / / / / / / / .
/end{array} /right./end{displaymath}

see Figure 5.10).

Figure 5.10: Regularized Variogram.
/begin{figure}/begin{center}
/mbox{/beginpicture
/setcoordinatesystem units <1.0...
...at 1.8 3.5
/put {$/gamma_l(h)$} at 3.3 2.8
/endpicture}
/end{center}/end{figure}

For distances $/vert{/boma h}/vert /leq /ell$, the regularized semi-variogram is parabolic. The cores of length $/ell$ overlap each other and the regularizing effect is strong enough to change the behavior of the model at the origin from linear to parabolic.

For distances $/vert{/boma h}/vert /geq /ell$, which, in practice, are the only distances observable, the linear behavior of the point model is preserved and the regularized model differs from the point model by a constant value

/begin{displaymath}/gamma_/ell({/boma h}) = /gamma({/boma h}) - /bar{/gamma}(/ell,/ell) =
r - /ell/3,
/ / r=/vert h/vert /geq /ell/ / ./end{displaymath}

If the regularized curve for $/vert{/boma h}/vert /geq /ell$ is projected towards ${/boma h} = {/boma 0}$, the ordinate axis is intercepted at $-/ell/3$ and this negative value is called a ``pseudo-negative nugget effect'' due to regularization.


next up previous contents
Next: Estimation of Resources Up: Regularization and Estimation Variance Previous: Regularization and Estimation Variance   Contents
Rudolf Dutter 2003-03-13