Cálculo de variogramas promedio Rubén Medinaceli O.

Tema 4
CÁLCULO DE VARIOGRAMAS PROMEDIO (𝜸
̅)

En todas las aplicaciones prácticas de la geoestadística, es necesario calcular el

variograma promedio entre puntos de una muestra o block con soporte 𝑊 con puntos de
una muestra o block con soporte 𝑊′, denotado por 𝛾̅ (𝑊; 𝑊′).
Esta presentación incluye algunos métodos para el cálculo de variogramas promedio.
En la siguiente figura:

𝑊
𝑊′
●●
𝑧′

●
𝑧

𝑊 es una muestra o block con soporte 𝑊

𝑊′ es una muestra o block con soporte 𝑊′
𝑧 es un punto del soporte 𝑊
𝑧′ es un punto del soporte 𝑊′

Por definición,

1
𝛾̅ (𝑊; 𝑊 ′ ) = ∫ ∫ 𝛾(𝑧; 𝑧 ′ )𝑑𝑧𝑑𝑧′
𝑊 𝑊′
𝑧 𝑒𝑛 𝑊 𝑧 ′ 𝑒𝑛 𝑊 ′

Esta integral puede ser resuelta por diferentes métodos, algunos de ellos se describen a
continuación.

1. Aproximación numérica
Es posible utilizar métodos convencionales de integración numérica para aproximar la
integral anotada en el punto anterior. Para ello es necesario discretizar los soportes
𝑊 (𝑡𝑜𝑚𝑎𝑛𝑑𝑜 𝑛 𝑝𝑢𝑛𝑡𝑜𝑠 𝑟𝑒𝑝𝑟𝑒𝑠𝑒𝑛𝑡𝑎𝑡𝑖𝑣𝑜𝑠) y 𝑊 ′ (𝑡𝑜𝑚𝑎𝑛𝑑𝑜 𝑛′ 𝑝𝑢𝑛𝑡𝑜𝑠 𝑟𝑒𝑝𝑟𝑒𝑠𝑒𝑛𝑡𝑎𝑡𝑖𝑣𝑜𝑠).
La integral puede, en general, aproximarse por la siguiente expresión:

1
Cálculo de variogramas promedio Rubén Medinaceli O.

𝑛 𝑛′
1
𝛾̅ (𝑊; 𝑊′) = ∑ ∑ 𝛾(𝑧𝑖 ; 𝑧𝑗′ )
𝑛 𝑛′
𝑖=1 𝑗=1

Considérese el siguiente caso particular,

𝑾′

𝑾 ● 𝑧1′ ● 𝑧2′

● 𝑧1 ● 𝑧2 ● 𝑧3 ● 𝑧3′

● 𝑧4 ● 𝑧5 ● 𝑧6 ● 𝑧7 ● 𝑧4′

● 𝑧8 ● 𝑧9 ● 𝑧10 ● 𝑧11 ● 𝑧12 ● 𝑧5′

● 𝑧6′

El soporte 𝑊 ha sido discretizado tomando 12 puntos representativos y el soporte 𝑊 ′ ha

sido discretizado con seis (6) puntos representativos.

En este caso, la aproximación de la integral viene dada por:

12 6
1
𝛾̅ (𝑊; 𝑊 ′ ) = ∑ ∑ 𝛾(𝑧𝑖 ; 𝑧𝑗′ )
12(6)
𝑖=1 𝑗=1

Ejercicio

Problema:

El modelo de variograma aceptado para un determinado depósito mineralizado es el que

se muestra en la siguiente figura.

El modelo es esférico con una meseta 𝐶 = 1,0 y un rango o alcance 𝑎 = 3 𝑢 (léase tres
unidades)

2
Cálculo de variogramas promedio Rubén Medinaceli O.

𝛾(ℎ⃗)

1,0

0,5

1 2 3 4
⃗⃗⃗ (𝑢)
ℎ

Calcular 𝛾̅ (𝑊; 𝑊 ′ ) en los siguientes casos:

Caso a)

𝑊 𝑊′
●𝑧 ● 𝑧′

2𝑢

Caso b)
𝑊

● 𝑧1 ● 𝑧2 ● 𝑧3

● 𝑧4 ● 𝑧5 ● 𝑧6 3𝑢

● 𝑧7 ● 𝑧8 ● 𝑧9

𝑧 ′■
𝑊′
3𝑢

3
Cálculo de variogramas promedio Rubén Medinaceli O.

Caso c)

● 𝑧1 ● 𝑧2 ● 𝑧3

● 𝑧4 𝑊 ′ ■ 𝑧5 ● 𝑧6 3𝑢
𝑧′

● 𝑧7 ● 𝑧8 ● 𝑧9

3𝑢

Caso d)

● 𝑧1 ● 𝑧2 ● 𝑧3

● 𝑧4 ● 𝑧5 ● 𝑧6 3𝑢

● 𝑧7 ● 𝑧8 ● 𝑧9

3𝑢

Solución:
Caso a)

𝑊 𝑊′
●𝑧 ● 𝑧′

2𝑢

4
Cálculo de variogramas promedio Rubén Medinaceli O.

Tanto 𝑊 como 𝑊 ′ son muestras puntuales separadas una distancia ℎ = 2 𝑢. 𝑊 es un punto

y 𝑊 ′ también es un punto. Por tanto, el variograma promedio de puntos en 𝑊 con puntos
en 𝑊 ′ , es igual al valor del variograma para una distancia ℎ = 2. Como el variograma es
esférico,

3 ℎ 1 ℎ 3 3 2 1 2 3
𝛾̅ (𝑊; 𝑊 ′ ) = 𝛾(𝑧; 𝑧 ′ ) = 𝛾(2) = 𝐶 { ( ) − ( ) } = 1,0 { ( ) − ( ) } = 0,852
2 𝑎 2 𝑎 2 3 2 3

Caso b)

● 𝑧1 ● 𝑧2 ● 𝑧3

● 𝑧4 ● 𝑧5 ● 𝑧6 3𝑢

● 𝑧7 ● 𝑧8 ● 𝑧9

𝑧 ′■
𝑊′
3𝑢

𝑊 es un block cuadrado de 3𝑢 × 3𝑢. Este block ha sido discretizado y queda representado

por nueve (9) puntos; cada punto es el centro de un block pequeño de 1𝑢 × 1𝑢 tal como
muestra el gráfico correspondiente.

𝑊 ′ es una muestra puntual ubicada en un vértice del block 𝑊 y está representada por un
solo punto (𝑧 ′ ).

El variograma promedio de puntos en 𝑊 y puntos en 𝑊 ′ viene dado por:

9
1 1
𝛾̅ (𝑊; 𝑊 ′ ) = ∑ 𝛾(𝑧𝑖 ; 𝑧 ′ ) = {𝛾(𝑧1 ; 𝑧 ′ ) + 𝛾(𝑧2 ; 𝑧 ′ ) + 𝛾(𝑧3 ; 𝑧 ′ ) + ⋯ + 𝛾(𝑧9 ; 𝑧 ′ )}
9×1 9
𝑖=1

El término 𝛾(𝑧𝑖 ; 𝑧 ′ ) es igual al valor del variograma para la distancia entre los puntos 𝑧𝑖 y 𝑧 ′ .

Por ejemplo,

𝛾(𝑧1 ; 𝑧 ′ ) = 𝛾 (√0,52 + 2,52 ) = 𝛾(2,55)

Si el modelo de variograma es esférico,

5
Cálculo de variogramas promedio Rubén Medinaceli O.

3 ℎ 1 ℎ 3 3 2,55 1 2,55 3
𝛾(2,55) = 𝐶 { ( ) − ( ) } = 1,0 { ( )− ( ) } = 0,968
2 𝑎 2 𝑎 2 3 2 3

La siguiente tabla muestra todos los valores de 𝛾(𝑧𝑖 ; 𝑧 ′ )

(𝒛𝒊 ; 𝒛′ ) 𝒉 𝜸(𝒉)
(𝑧1 ; 𝑧 ′ ) √0,52 + 2,52 = 2,55 0,968
(𝑧2 ; 𝑧 ′ ) √1,52 + 2,52 = 2,91 0,999
(𝑧3 ; 𝑧 ′ ) √2,52 + 2,52 = 3,53 1,000
(𝑧4 ; 𝑧 ′ ) √0,52 + 1,52 = 1,58 0,717
(𝑧5 ; 𝑧 ′ ) √1,52 + 1,52 = 2,12 0,883
(𝑧6 ; 𝑧 ′ ) √2,52 + 1,52 = 2,91 0,999
(𝑧7 ; 𝑧 ′ ) √0,52 + 0,52 = 0,71 0,348
(𝑧8 ; 𝑧 ′ ) √1,52 + 0,52 = 1,58 0,717
(𝑧9 ; 𝑧 ′ ) √2,52 + 0,52 = 2,55 0,968
Ʃ 7,599

Por tanto,

1
𝛾̅ (𝑊; 𝑊 ′ ) = (7,599) = 0,844
9

Caso c)

● 𝑧1 ● 𝑧2 ● 𝑧3

● 𝑧4 𝑊 ′ ■ 𝑧5 ● 𝑧6 3𝑢
𝑧′

● 𝑧7 ● 𝑧8 ● 𝑧9

3𝑢

𝑊 es un block cuadrado de 3𝑢 × 3𝑢. Este block ha sido discretizado y queda representado

por nueve (9) puntos; cada punto es el centro de un block pequeño de 1𝑢 × 1𝑢 tal como
muestra el gráfico correspondiente.

6
Cálculo de variogramas promedio Rubén Medinaceli O.

𝑊 ′ es una muestra puntual ubicada en el centro del block 𝑊 y está representada por un
solo punto (𝑧 ′ ).

El variograma promedio de puntos en 𝑊 y puntos en 𝑊 ′ viene dado por:

9
′)
1 1
𝛾̅ (𝑊; 𝑊 = ∑ 𝛾(𝑧𝑖 ; 𝑧 ′ ) = {𝛾(𝑧1 ; 𝑧 ′ ) + 𝛾(𝑧2 ; 𝑧 ′ ) + 𝛾(𝑧3 ; 𝑧 ′ ) + ⋯ + 𝛾(𝑧9 ; 𝑧 ′ )}
9×1 9
𝑖=1

La siguiente tabla muestra todos los valores de 𝛾(𝑧𝑖 ; 𝑧 ′ )

(𝒛𝒊 ; 𝒛′ ) 𝒉 𝜸(𝒉)
(𝑧1 ; 𝑧 ′ ) √1,02 + 1,02 = 1,41 0,653
(𝑧2 ; 𝑧 ′ ) = 1,00 0,481
(𝑧3 ; 𝑧 ′ ) √1,02 + 1,02 = 1,41 0,653
(𝑧4 ; 𝑧 ′ ) = 1,00 0,481
(𝑧5 ; 𝑧 ′ ) = 0,00 0,000
(𝑧6 ; 𝑧 ′ ) = 1,00 0,481
(𝑧7 ; 𝑧 ′ ) √1,02 + 1,02 = 1,41 0,653
(𝑧8 ; 𝑧 ′ ) = 1,00 0,481
(𝑧9 ; 𝑧 ′ ) 2 2
√1,0 + 1,0 = 1,41 0,653
Ʃ 4,536

Por tanto,

1
𝛾(𝑊; 𝑊 ′ ) = (4,536) = 0,504
9

Caso d)

● 𝑧1 ● 𝑧2 ● 𝑧3

● 𝑧4 ● 𝑧5 ● 𝑧6 3𝑢

● 𝑧7 ● 𝑧8 ● 𝑧9

3𝑢

7
Cálculo de variogramas promedio Rubén Medinaceli O.

𝑊 es un block cuadrado de 3𝑢 × 3𝑢. Este block ha sido discretizado y queda representado

por nueve (9) puntos; cada punto es el centro de un block pequeño de 1𝑢 × 1𝑢 tal como
muestra el gráfico correspondiente.

El variograma promedio de puntos en el block 𝑊 con puntos en el mismo block 𝑊, viene

dado por:

9 9
1
𝛾̅ (𝑊; 𝑊) = ∑ ∑ 𝛾(𝑧𝑖 ; 𝑧𝑗 )
9 ×9
𝑖=1 𝑗=1

1
𝛾̅ (𝑊; 𝑊) = {𝛾(𝑧1 ; 𝑧1 ) + 𝛾(𝑧1 ; 𝑧2 ) + ⋯ + 𝛾(𝑧1 ; 𝑧9 ) + 𝛾(𝑧2 ; 𝑧1 ) + ⋯ + 𝛾(𝑧2 ; 𝑧9 ) + ⋯
81
+ 𝛾(𝑧9 ; 𝑧1 ) + ⋯ + 𝛾(𝑧9 ; 𝑧9 )}

Nótese que:

𝛾(𝑧1 ; 𝑧1 ) = 𝛾(𝑧2 ; 𝑧2 ) = 𝛾(𝑧3 ; 𝑧3 ) = 𝛾(𝑧4 ; 𝑧4 ) = 𝛾(𝑧5 ; 𝑧5 ) = 𝛾(𝑧6 ; 𝑧6 ) = 𝛾(𝑧7 ; 𝑧7 ) = 𝛾(𝑧8 ; 𝑧8 )

𝛾(𝑧9 ; 𝑧9 ) = 𝛾(0)

𝛾(𝑧1 ; 𝑧2 ) = 𝛾(𝑧2 ; 𝑧1 ) = 𝛾(𝑧2 ; 𝑧3 ) = 𝛾(𝑧3 ; 𝑧2 ) = 𝛾(𝑧4 ; 𝑧5 ) = 𝛾(𝑧5 ; 𝑧4 ) = 𝛾(𝑧5 ; 𝑧6 ) = 𝛾(𝑧6 ; 𝑧5 ) =

𝛾(𝑧7 ; 𝑧8 ) = 𝛾(𝑧8 ; 𝑧7 ) = 𝛾(𝑧8 ; 𝑧9 ) = 𝛾(𝑧9 ; 𝑧8 ) = 𝛾(𝑧1 ; 𝑧4 ) = 𝛾(𝑧4 ; 𝑧1 ) = 𝛾(𝑧2 ; 𝑧5 ) = 𝛾(𝑧5 ; 𝑧2 ) =
𝛾(𝑧3 ; 𝑧6 ) = 𝛾(𝑧6 ; 𝑧3 ) = 𝛾(𝑧4 ; 𝑧7 ) = 𝛾(𝑧7 ; 𝑧4 ) = 𝛾(𝑧5 ; 𝑧8 ) = 𝛾(𝑧8 ; 𝑧5 ) = 𝛾(𝑧6 ; 𝑧9 ) = 𝛾(𝑧9 ; 𝑧6 ) =
𝛾(1)

𝛾(𝑧1 ; 𝑧3 ) = 𝛾(𝑧3 ; 𝑧1 ) = 𝛾(𝑧4 ; 𝑧6 ) = 𝛾(𝑧6 ; 𝑧4 ) = 𝛾(𝑧7 ; 𝑧9 ) = 𝛾(𝑧9 ; 𝑧7 ) = 𝛾(𝑧1 ; 𝑧7 ) = 𝛾(𝑧7 ; 𝑧1 ) =

𝛾(𝑧2 ; 𝑧8 ) = 𝛾(𝑧8 ; 𝑧2 ) = 𝛾(𝑧3 ; 𝑧9 ) = 𝛾(𝑧9 ; 𝑧3 ) = 𝛾(2)

𝛾(𝑧1 ; 𝑧5 ) = 𝛾(𝑧5 ; 𝑧1 ) = 𝛾(𝑧2 ; 𝑧4 ) = 𝛾(𝑧4 ; 𝑧2 ) = 𝛾(𝑧2 ; 𝑧6 ) = 𝛾(𝑧6 ; 𝑧2 ) = 𝛾(𝑧3 ; 𝑧5 ) = 𝛾(𝑧5 ; 𝑧3 ) =

𝛾(𝑧4 ; 𝑧8 ) = 𝛾(𝑧8 ; 𝑧4 ) = 𝛾(𝑧5 ; 𝑧7 ) = 𝛾(𝑧7 ; 𝑧5 ) = 𝛾(𝑧5 ; 𝑧9 ) = 𝛾(𝑧9 ; 𝑧5 ) = 𝛾(𝑧6 ; 𝑧8 ) = 𝛾(𝑧8 ; 𝑧6 ) =
𝛾 (√12 + 12 ) = 𝛾(√2) = 𝛾(1,41)

𝛾(𝑧1 ; 𝑧6 ) = 𝛾(𝑧6 ; 𝑧1 ) = 𝛾(𝑧3 ; 𝑧4 ) = 𝛾(𝑧4 ; 𝑧3 ) = 𝛾(𝑧4 ; 𝑧9 ) = 𝛾(𝑧9 ; 𝑧4 ) = 𝛾(𝑧6 ; 𝑧7 ) = 𝛾(𝑧7 ; 𝑧6 ) =

𝛾(𝑧1 ; 𝑧8 ) = 𝛾(𝑧8 ; 𝑧1 ) = 𝛾(𝑧2 ; 𝑧7 ) = 𝛾(𝑧7 ; 𝑧2 ) = 𝛾(𝑧2 ; 𝑧9 ) = 𝛾(𝑧9 ; 𝑧2 ) = 𝛾(𝑧3 ; 𝑧8 ) = 𝛾(𝑧8 ; 𝑧3 ) =
𝛾 (√12 + 22 ) = 𝛾(√5) = 𝛾(2,24)

𝛾(𝑧1 ; 𝑧9 ) = 𝛾(𝑧9 ; 𝑧1 ) = 𝛾(𝑧3 ; 𝑧7 ) = 𝛾(𝑧7 ; 𝑧3 ) = 𝛾(√22 + 22 ) = 𝛾(√8) = 𝛾(2,83)

Consecuentemente,

1
𝛾̅ (𝑊; 𝑊) = {9𝛾(0) + 24𝛾(1) + 12𝛾(2) + 16𝛾(1,41) + 16𝛾(2,24) + 4𝛾(2,83)}
81

Siendo el modelo de variograma esférico con meseta (𝐶) = 1,0 y rango o alcance (𝑎) = 3 𝑢

1
𝛾̅ (𝑊; 𝑊) = {9(0) + 24(0,481) + 12(0,852) + 16(0,653) + 16(0,912) + 4𝛾(0,995)}
81

𝛾̅ (𝑊; 𝑊) = 0,627

8
Cálculo de variogramas promedio Rubén Medinaceli O.

Ejercicio

Problema:

Dada la siguiente información,

𝑊′ 30 m.

40 m.

30 m.

𝑊 = Block mineralizado en un depósito tipo veta

𝑊 ′ = Chimenea en el block mineralizado superior

El variograma experimental de depósito mineralizado es:

𝒉(𝒎) 10 20 30 40 50 60 70
𝜸(𝒉) 0,75 1,55 2,00 3,25 3,50 4,75 5,00
𝑵 20 25 30 35 30 25 20

Calcular 𝛾̅ (𝑊 ′ ; 𝑊)

Solución:

Inicialmente, se debe obtener un modelo de variograma.

El siguiente gráfico muestra el variograma experimental.

9
Cálculo de variogramas promedio Rubén Medinaceli O.

Variograma experimental
ϒ(h) 6

3
5 5,3
2

0
10 20 30 40 50 60 70
h (m)

70

La forma del variograma experimental sugiere un modelo lineal.

Vale decir,

𝛾 (ℎ ) = 𝛼 ℎ

Donde,

5,3 − 0
𝛼= = 0,076
70 − 0

Por tanto,

𝛾(ℎ) = 0,076 ℎ

La siguiente tabla muestra valores del variograma experimental y los correspondientes del
modelo lineal:

10
Cálculo de variogramas promedio Rubén Medinaceli O.

𝒉 𝑽𝒂𝒓𝒊𝒐𝒈𝒓𝒂𝒎𝒂 𝑴𝒐𝒅𝒆𝒍𝒐 𝒅𝒆 𝑵
𝒆𝒙𝒑𝒆𝒓𝒊𝒎𝒆𝒏𝒕𝒂𝒍 𝒗𝒂𝒓𝒊𝒐𝒈𝒓𝒂𝒎𝒂
10 0,75 0,76
20 1,55 1,52
30 2,00 2,28
40 3,25 3,04
50 3,50 3,80
60 4,75 4,56
70 5,00 5,32

El siguiente gráfico muestra la bondad del ajuste efectuado,

Variograma experimental vs. modelo

6
ϒ(h)

0
10 20 30 40 50 60 70
h (m)

El ajuste efectuado es bueno.

Ahora si es posible obtener el variograma promedio solicitado.

Para el efecto, ambos soportes deben ser discretizados tal como muestra la siguiente figura,

11
Cálculo de variogramas promedio Rubén Medinaceli O.

𝑊′

■𝑧3′ 10 𝑚

■ 𝑧2′ 10 𝑚

■ 𝑧1′ 10 𝑚
𝑊

■ ■ ■ ■ 10 𝑚
𝑧1 𝑧2 𝑧3 𝑧4

■ ■ ■ ■ 10 𝑚
𝑧5 𝑧6 𝑧7 𝑧8

■ ■ ■ ■ 10 𝑚
𝑧9 𝑧10 𝑧11 𝑧12

10 𝑚 10 𝑚 10 𝑚 10 𝑚

12 3
1
𝛾̅ (𝑊; 𝑊′) = ∑ ∑ 𝛾(𝑧𝑖 ; 𝑧𝑗′ )
12 × 3
𝑖=1 𝑗=1

1
𝛾̅ (𝑊; 𝑊 ′ ) = {𝛾(𝑧1 ; 𝑧1′ ) + 𝛾(𝑧2 ; 𝑧1′ ) + 𝛾(𝑧3 ; 𝑧1′ ) + 𝛾(𝑧4 ; 𝑧1′ ) + 𝛾(𝑧5 ; 𝑧1′ ) + 𝛾(𝑧6 ; 𝑧1′ ) + 𝛾(𝑧7 ; 𝑧1′ )
36
+ 𝛾(𝑧8 ; 𝑧1′ ) + 𝛾(𝑧9 ; 𝑧1′ ) + 𝛾(𝑧10 ; 𝑧1′ ) + 𝛾(𝑧11 ; 𝑧1′ ) + 𝛾(𝑧12 ; 𝑧1′ ) + 𝛾(𝑧1 ; 𝑧2′ )
+ 𝛾(𝑧2 ; 𝑧2′ ) + 𝛾(𝑧3 ; 𝑧2′ ) + 𝛾(𝑧4 ; 𝑧2′ ) + 𝛾(𝑧5 ; 𝑧2′ ) + 𝛾(𝑧6 ; 𝑧2′ ) + 𝛾(𝑧7 ; 𝑧2′ )
+ 𝛾(𝑧8 ; 𝑧2′ ) + 𝛾(𝑧9 ; 𝑧2′ ) + 𝛾(𝑧10 ; 𝑧2′ ) + 𝛾(𝑧11 ; 𝑧2′ ) + 𝛾(𝑧12 ; 𝑧2′ ) + 𝛾(𝑧1 ; 𝑧3′ )
+ 𝛾(𝑧2 ; 𝑧3′ ) + 𝛾(𝑧3 ; 𝑧3′ ) + 𝛾(𝑧4 ; 𝑧3′ ) + 𝛾(𝑧5 ; 𝑧3′ ) + 𝛾(𝑧6 ; 𝑧3′ ) + 𝛾(𝑧7 ; 𝑧3′ )
+ 𝛾(𝑧8 ; 𝑧3′ ) + 𝛾(𝑧9 ; 𝑧3′ ) + 𝛾(𝑧10 ; 𝑧3′ ) + 𝛾(𝑧11 ; 𝑧3′ ) + 𝛾(𝑧12 ; 𝑧3′ )}

12
Cálculo de variogramas promedio Rubén Medinaceli O.

1
𝛾̅ (𝑊; 𝑊 ′ ) = {𝛾(√52 + 102 ) + 𝛾(√152 + 102 ) + 𝛾(√252 + 102 ) + 𝛾(√352 + 102 )
36
+ 𝛾(√52 + 202 ) + 𝛾(√152 + 202 ) + 𝛾(√252 + 202 ) + 𝛾(√352 + 202 )
+ 𝛾(√52 + 302 ) + 𝛾(√152 + 302 ) + 𝛾(√252 + 302 ) + 𝛾(√352 + 302 )
+ 𝛾(√52 + 202 ) + 𝛾(√152 + 202 ) + 𝛾(√252 + 202 ) + 𝛾(√352 + 202 )
+ 𝛾(√52 + 302 ) + 𝛾(√152 + 302 ) + 𝛾(√252 + 302 ) + 𝛾(√352 + 302 )
+ √52 + 402 + 𝛾(√152 + 402 ) + 𝛾(√252 + 402 ) + 𝛾(√352 + 402 )
+ √52 + 302 + 𝛾(√152 + 302 ) + 𝛾(√252 + 302 ) + 𝛾(√352 + 302 )
+ 𝛾(√52 + 402 ) + 𝛾(√152 + 402 ) + 𝛾(√252 + 402 ) + 𝛾(√352 + 402 )
+ 𝛾(√52 + 502 ) + 𝛾(√152 + 502 ) + 𝛾(√252 + 502 ) + 𝛾(√352 + 502 ) }

1
𝛾̅ (𝑊; 𝑊 ′ ) = {𝛾(11,18) + 𝛾(18,03) + 𝛾(26,93) + 𝛾(36,40) + 𝛾(20,61) + 𝛾(25,00)
36
+ 𝛾(32,02) + 𝛾(40,31) + 𝛾(30,41) + 𝛾(33,54) + 𝛾(39,05) + 𝛾(46,10)
+ 𝛾(20,61) + 𝛾(25,00) + 𝛾(32,02) + 𝛾(40,31) + 𝛾(30,41) + 𝛾(33,54)
+ 𝛾(39,05) + 𝛾(46,10) + 𝛾(40,31) + 𝛾(42,72) + 𝛾(47,17) + 𝛾(53,15)
+ 𝛾(30,41) + 𝛾(33,54) + 𝛾(39,05) + 𝛾(46,10) + 𝛾(40,31) + 𝛾(42,72)
+ 𝛾(47,17) + 𝛾(53,15) + 𝛾(50,25) + 𝛾(52,20) + 𝛾(55,90) + 𝛾(61,03)}

1
𝛾̅ (𝑊; 𝑊 ′ ) = {0,850 + 1,370 + 2,047 + 2,766 + 1,566 + 1,900 + 2,433 + 3,064 + 2,311
36
+ 2,549 + 2,698 + 3,504 + 1,566 + 1,900 + 2,433 + 3,064 + 2,311 + 2,549
+ 2,968 + 3,504 + 3,064 + 3,247 + 3,585 + 4,039 + 2,311 + 2,549 + 2,968
+ 3,504 + 3,064 + 3,247 + 3,585 + 4,039 + 3,819 + 3,967 + 4,248 + 4,638}

𝛾̅ (𝑊; 𝑊 ′ ) = 2,867

2. Funciones auxiliares

Afortunadamente, se dispone de funciones (expresiones algebraicas o ábacos) que

permiten calcular variogramas promedio, 𝛾̅ (𝑊; 𝑊 ′ ), para formas comunes de los soportes
𝑊 𝑦 𝑊 ′ y modelos de variograma más frecuentes.

Estas funciones reciben el nombre de funciones auxiliares. Se cuenta con funciones

auxiliares en una dimensión, dos dimensiones y tres dimensiones.

2.1. Funciones auxiliares en una dimensión

Si ̅̅̅̅
𝐴𝐵 es un segmento de recta de longitud 𝐿; es posible definir las siguientes funciones
auxiliares:

2.1.1. Función auxiliar X

𝐿
1
̅̅̅̅)
𝑋(𝐿) = ∫ 𝛾(𝐴; 𝑧 ′ )𝑑𝑧 ′ = 𝛾̅ (𝐴; 𝐴𝐵
𝐿
0

13
Cálculo de variogramas promedio Rubén Medinaceli O.

𝐴 𝑧′ 𝐵

La función auxiliar 𝑿(𝑳) permite obtener el variograma promedio de un extremo del

̅̅̅̅ (𝐴 ó 𝐵) con puntos en el segmento 𝐴𝐵
segmento 𝐴𝐵 ̅̅̅̅.

2.1.2. Función auxiliar F

𝐿 𝐿
1
𝐹(𝐿) = 2 ∫ ∫ 𝛾(𝑧; 𝑧 ′ )𝑑𝑧 𝑑𝑧 ′ = 𝛾̅ (̅̅̅̅
𝐴𝐵; ̅̅̅̅
𝐴𝐵 )
𝐿
0 0

𝐴 𝑧 𝑧′ 𝐵

La función auxiliar 𝑭(𝑳) permite obtener el variograma promedio de puntos en el segmento

̅̅̅̅
𝐴𝐵 con puntos en el mismo segmento ̅̅̅̅𝐴𝐵.

2.1.3. Cálculo de las funciones auxiliares en una dimensión

Las expresiones algebraicas que permiten calcular las funciones auxiliares en una
dimensión son:

Modelo de variograma lineal

𝑋(𝐿) = 𝐿⁄2

𝐹(𝐿) = 𝐿/3

Modelo de variograma esférico

3 𝐿 1 𝐿 3
𝑋(𝐿) = 𝐶 { − ( ) } ; ∀𝐿 ≤𝑎
4𝑎 8 𝑎

3𝑎
= 𝐶 {1 − } ; ∀𝐿 >𝑎
8 𝐿

14
Cálculo de variogramas promedio Rubén Medinaceli O.

1 𝐿 1 𝐿 3
𝐹(𝐿) = 𝐶 { − ( ) } ; ∀𝐿≤𝑎
2 𝑎 20 𝑎

3𝑎 1 𝑎 2
= 𝐶 {1 − + ( ) } ;∀ 𝐿 > 𝑎
4 𝐿 5 𝐿

Donde:

𝐶 = Meseta del variograma

𝑎 = Rango o alcance del variograma

Modelo de variograma exponencial

𝑎 𝐿
𝑋(𝐿) = 𝐶 {1 − (1 − 𝑒 − 𝑎 )}
𝐿
2𝑎 𝑎 𝐿
𝐹(𝐿) = 𝐶 {1 − [1 − (1 − 𝑒 − 𝑎 )]}
𝐿 𝐿

Donde:

𝐶 = Meseta del variograma

𝑎 = Rango o alcance del variograma

Ejercicio

Problema:

En un depósito mineralizado tipo veta, el modelo de variograma es esférico con una meseta
𝐶 = 0,5 y un rango o alcance 𝑎 = 10:

𝛾(ℎ⃗)

0,50

0,25

5,0 10,0 15,0

⃗⃗⃗ (𝑚)
ℎ

15
Cálculo de variogramas promedio Rubén Medinaceli O.

Obtener el variograma promedio, 𝛾(𝑊; 𝑊 ′ ), para los casos que se muestran a continuación:

Caso a)

𝑊 𝑊′
●

𝐿 = 15 𝑚

Caso b)

𝑊 𝑊′
●

𝐿⁄2 = 7,5

𝐿 = 15 𝑚

Caso c)

𝑊 𝑊′
●

𝐿⁄3 = 5 𝑚 𝐿 = 15 𝑚

Solución:

Caso a)

𝑊 𝑊′
●

𝐿 = 15 𝑚

La muestra 𝑊 se encuentra en un extremo de una galería (corrida) 𝑊 ′ de 15 𝑚 de

longitud; consecuentemente,

16
Cálculo de variogramas promedio Rubén Medinaceli O.

3𝑎 3 10
𝛾̅ (𝑊; 𝑊 ′ ) = 𝑋(𝐿) = 𝑋(15) = 𝐶 {1 − } = 0,5 {1 − } = 0,375
8𝐿 8 15

Caso b)

𝑊 𝑊′
●

𝐿⁄2 = 7,5 𝑚

𝐿 = 15 𝑚

La muestra 𝑊 se encuentra al centro de una galería (corrida) de 15 𝑚 de longitud.

No existe una función auxiliar directa. Sin embargo, es posible recurrir a una
descomposición del caso para poder aplicar la función auxiliar 𝑋(𝐿) y luego efectuar
una ponderación tomando como factor las longitudes correspondientes.

𝑊′ 𝑊
● 𝑾′
●

𝐿⁄2 = 7,5 𝑚 𝐿⁄2 = 7,5 𝑚

𝐿 𝐿 𝐿 𝐿 𝐿 𝐿
𝑋 ( ) + 𝑋 ( ) 2 { 𝑋 ( )} 𝐿 3 𝐿 1 𝐿 3
2 2 2 2 2 2
𝛾̅ (𝑊; 𝑊 ′ ) = = = 𝑋 ( ) = 𝑋(7,5) = 𝐶 { − ( ) }
𝐿 𝐿 𝐿 2 4 𝑎 8 𝑎
2+2

′)
3 7,5 1 7,5 3
𝛾̅ (𝑊; 𝑊 = 0,5 { − ( ) } = 0,255
4 10 8 10

Caso c)

𝑊 𝑊′
●

𝐿⁄3 = 5 𝑚 𝐿 = 15 𝑚

17
Cálculo de variogramas promedio Rubén Medinaceli O.

La muestra se encuentra alejada 5𝑚 del extremo de una galería (corrida) de 15𝑚 de

longitud. Nuevamente, se recurre a una composición y posterior ponderación para
posibilitar la aplicación de la función auxiliar 𝑋(𝐿); esto es,

𝐿 𝐿 𝐿 𝐿
(3 + 𝐿) 𝑋 (3 + 𝐿) − 3 𝑋 (3) 20 𝑋(20) − 5 𝑋 (5)
′)
𝛾̅ (𝑊; 𝑊 = =
𝐿 𝐿 15
3+𝐿−3

3 10
𝑋(20) = 0,5 {1 − } = 0,406
8 20

3 5 1 5 3
𝑋(5) = 0,5 { − ( ) } = 0,179
4 10 8 10

20(0,406) − 5(0,179)
𝛾̅ (𝑊; 𝑊 ′ ) = = 0,482
15

Ejercicio

Problema:

Obtener el 𝛾̅ (𝑊; 𝑊 ′ ) para el caso que se muestra en el siguiente gráfico.

El modelo de variograma es esférico con una meseta 𝐶 = 1,5 y un rango o alcance 𝑎 =

20 𝑚.

𝑊 𝑊′

20 𝑚 10 𝑚 20 𝑚

Solución:

No existe una función auxiliar directa para calcular el variograma promedio solicitado. En
este caso, es posible combinar el método de la discretización con el uso de funciones
auxiliares; en esta línea, se ha discretizado el soporte 𝑊 el mismo que queda representado
por los puntos 𝑃1 , 𝑃2 , 𝑃3 𝑦 𝑃4 .

18
Cálculo de variogramas promedio Rubén Medinaceli O.

𝑊 𝑊′
● ● ● ●
𝑃1 𝑃2 𝑃3 𝑃4

5𝑚 5𝑚 5𝑚 5𝑚 10 𝑚 20 𝑚

El variograma promedio de puntos en el soporte 𝑊 con puntos en el soporte 𝑊 ′ es

aproximado por el promedio ponderado del variograma de cada punto 𝑃𝑖 con el soporte 𝑊 ′ ;
el factor de ponderación será la distancia de influencia de cada punto 𝑃𝑖 ; vale decir,

5 𝛾̅ (𝑃1 ; 𝑊 ′ ) + 5 𝛾̅ (𝑃2 ; 𝑊 ′ ) + 5 𝛾̅ (𝑃3 ; 𝑊 ′ ) + 5 𝛾̅ (𝑃4 ; 𝑊 ′ )

𝛾̅ (𝑊; 𝑊 ′ ) =
5+5+5+5

47,5 𝑋(47,5) − 27,5 𝑋(27,5)

𝛾̅ (𝑃1 ; 𝑊 ′ ) =
47,5 − 27,5

3 20
𝑋(47,5) = 1,5 {1 − } = 1,263
8 47,5

3 20
𝑋(27,5) = 1,5 {1 − } = 1,090
8 27,5

47,5 (1,263) − 27,5 (1,090)

𝛾̅ (𝑃1 ; 𝑊 ′ ) = = 1,500
47,5 − 27,5

42,5 𝑋(42,5) − 22,5 𝑋(22,5)

𝛾̅ (𝑃2 ; 𝑊 ′ ) =
42,5 − 22,5

3 20
𝑋(42,5) = 1,5 {1 − } = 1,235
8 42,5

3 20
𝑋(22,5) = 1,5 {1 − } = 1,000
8 22.5

42,5(1,235) − 22,5(1,000)
𝛾̅ (𝑃2 ; 𝑊 ′ ) = = 1,499
42,5 − 22,5

37,5 𝑋(37,5) − 17,5 𝑋(17,5)

𝛾̅ (𝑃3 ; 𝑊 ′ ) =
37,5 − 17,5

3 20
𝑋(37,5) = 1,5 {1 − } = 1,200
8 37,5

19
Cálculo de variogramas promedio Rubén Medinaceli O.

3 17,5 1 17,5 3
𝑋(17,5) = 1,5 { − ( ) } = 0,858
4 20 8 20

37,5 (1,200) − 17,5 (0,858)

𝛾̅ (𝑃3 ; 𝑊 ′ ) = = 1,499
20

32,5 𝑋(32,5) − 12,5 𝑋(12,5)

𝛾̅ (𝑃4 ; 𝑊 ′ ) =
32,5 − 12,5

3 20
𝑋(32,5) = 1,5 {1 − } = 1,153
8 32,5

3 12,5 1 12,5 3
𝑋(12,5) = 1,5 { − ( ) } = 0,658
4 20 8 20

32,5 (1,153) − 12,5 (0,658)

𝛾̅ (𝑃4 ; 𝑊 ′ ) = = 1,462
20

Finalmente,

5(1,500) + 5(1,499) + 5(1,499) + 5(1,462)

𝛾̅ (𝑊; 𝑊 ′ ) = = 1,490
5+5+5+5

2.2. Funciones auxiliares en dos dimensiones

Si 𝐴𝐵𝐶𝐷 es un block rectangular con vértices 𝐴, 𝐵, 𝐶 𝑦 𝐷 tal que,

̅̅̅̅
𝐴𝐵 = ̅̅̅̅
𝐶𝐷 = 𝐿

̅̅̅̅ = 𝐵𝐷
𝐴𝐶 ̅̅̅̅ = 𝑑

Es posible definir las siguientes funciones auxiliares:

20
Cálculo de variogramas promedio Rubén Medinaceli O.

2.2.1. Función auxiliar 𝜶

̅̅̅̅ ; 𝐵𝐷
𝛼(𝐿; 𝑑) = 𝛾̅ (𝐴𝐶 ̅̅̅̅ )

𝐴 𝐿 𝐵

𝑧′

𝐶 𝐷

𝛼(𝑑; 𝐿) = 𝛾̅ (̅̅̅̅
𝐴𝐵; ̅̅̅̅
𝐶𝐷 )

𝐴 𝑧 𝐿 𝐵

𝐶 𝑧′ 𝐷

En general, la función auxiliar 𝜶 permite obtener el variograma promedio de puntos en una

arista con puntos en la arista opuesta; las aristas son identificadas por el orden de los
parámetros de la función.

2.2.2. Función auxiliar 𝑿

𝑋(𝐿; 𝑑) = 𝛾̅ (̅̅̅̅
𝐴𝐶 ; 𝐴𝐵𝐶𝐷) = 𝛾̅ (̅̅̅̅
𝐵𝐷 ; 𝐴𝐵𝐶𝐷)

𝐴 𝐿 𝐵

𝑑 𝑧′

𝐶 𝐷
̅̅̅̅; 𝐴𝐵𝐶𝐷) = 𝛾̅ (𝐶𝐷
𝑋(𝑑; 𝐿) = 𝛾̅ (𝐴𝐵 ̅̅̅̅ : 𝐴𝐵𝐶𝐷)

21
Cálculo de variogramas promedio Rubén Medinaceli O.

𝐴 𝑧 𝐿 𝐵

𝑑 𝑧′

𝐶 𝑧′ 𝐷

En general, la función auxiliar 𝑿 permite obtener el variograma promedio de puntos en una

arista con puntos en todo el block; la arista es identificada por el orden de los parámetros
de la función.

2.2.3. Función auxiliar 𝑯

𝐻(𝐿; 𝑑) = 𝐻(𝑑; 𝐿) = 𝛾̅ (̅̅̅̅

𝐴𝐶 ; ̅̅̅̅
𝐴𝐵 ) = 𝛾̅ (̅̅̅̅
𝐵𝐷 ; ̅̅̅̅
𝐷𝐶 )

𝐴 𝐿 𝑧′ 𝐵

𝐶 𝐷

𝐻(𝐿; 𝑑) = 𝐻(𝑑; 𝐿) = 𝛾̅ (𝐴; 𝐴𝐵𝐶𝐷) = 𝛾̅ (𝐵; 𝐴𝐵𝐶𝐷) = 𝛾̅ (𝐶; 𝐴𝐵𝐶𝐷) = 𝛾̅ (𝐷; 𝐴𝐵𝐶𝐷)

𝑧
𝐴 𝐿 𝐵
●

𝑑 𝑧′

𝐶 𝐷

22
Cálculo de variogramas promedio Rubén Medinaceli O.

La función auxiliar 𝑯 es múltiple. Permite obtener el variograma promedio de puntos en una

arista con puntos en una arista adyacente; igualmente, permite obtener el variograma
promedio de un vértice con puntos en todo el block. Nótese que el orden de los parámetros
de esta función no es importante.

2.2.4. Función auxiliar 𝑭

𝐹(𝐿; 𝑑) = 𝐹(𝑑: 𝐿) = 𝛾̅ (𝐴𝐵𝐶𝐷; 𝐴𝐵𝐶𝐷)

𝐴 𝐿 𝐵

𝑧′
𝑑

𝐶 𝐷

La función auxiliar 𝑭 permite obtener el variograma promedio de puntos en el block con

puntos en el mismo block. El orden de los parámetros de esta función auxiliar no es
importante.

2.2.5. Estimación de las funciones auxiliares en dos dimensiones

Los valores de las funciones auxiliares en dos dimensiones se obtienen de los ábacos que
se muestran en las siguientes páginas. Se incluyen ábacos para el modelo de variograma
esférico y para el modelo de variograma exponencial. La fuente de dichos ábacos es el libro
Mining Geostatistics de A.G. Journel y Ch.J. Huijbregts, publicado por Academic Press en
1978.

23
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 1

𝑳 𝜶(𝟎; 𝒍)
𝟏 𝟏
𝒍 𝜸 (𝒍) = 𝑪 {𝜶(𝑳; 𝒍) − 𝜶(𝟎; 𝒍)}
𝑪 𝒍
𝑀′ 𝑙

Modelo de variograma esférico

24
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 2

𝑳
𝟏 𝑳 𝒅
𝒅 𝑿 ( ; )
𝑪 𝒂 𝒂
𝑧′

Modelo de variograma esférico

25
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 3

𝑳 𝑧′ 𝑳
𝟏 𝑳 𝒅
𝒅 𝒅 𝑯 ( ; )
𝑪 𝒂 𝒂
𝑧′

𝑧 𝑧

Modelo de variograma esférico

26
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 4

𝑳
𝟏 𝑳 𝒅
𝒅 𝑧′ 𝑭 ( ; )
𝑪 𝒂 𝒂

Modelo de variograma esférico

27
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 5

𝑳 𝜶(𝟎; 𝒍)
𝟏 𝟏
𝒍 𝜸 (𝒍) = 𝑪 {𝜶(𝑳; 𝒍) − 𝜶(𝟎; 𝒍)}
𝑪 𝒍
𝑀′ 𝑙

Modelo de variograma exponencial

28
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 6

𝑳
𝟏 𝑳 𝒅
𝒅 𝑿 ( ; )
𝑪 𝒂 𝒂
𝑧′

Modelo de variograma exponencial

29
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 7

𝑳 𝑧′ 𝑳
𝟏 𝑳 𝒅
𝒅 𝒅 𝑯 ( ; )
𝑪 𝒂 𝒂
𝑧′

𝑧 𝑧

Modelo de variograma exponencial

30
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 8

𝑳
𝟏 𝑳 𝒅
𝒅 𝑧′ 𝑭 ( ; )
𝑪 𝒂 𝒂

Modelo de variograma exponencial

31
Cálculo de variogramas promedio Rubén Medinaceli O.

Ejercicio

Problema:

Estimar 𝛾̅ (𝑊; 𝑊 ′ ) para el caso que se muestra en la siguiente figura. Nótese que el soporte
𝑊 ′ es un block de 30 𝑚 × 40 𝑚 sobre una veta mineralizada y el soporte 𝑊 es una
chimenea lateral de 30 𝑚 de longitud.

El variograma es esférico con alcance 𝑎 = 10 𝑚 y meseta 𝐶 = 0,5.

𝐿 = 40 𝑚

𝑑 = 30 𝑚 𝑊′

Solución:

El variograma solicitado puede obtenerse directamente con una función auxiliar, esto es:

1 𝐿 𝑑 1 40 30 1
𝛾̅ (𝑊; 𝑊 ′ ) = 𝑋(𝐿; 𝑑) = 𝐶 { 𝑋 ( ; )} = 𝐶 { 𝑋 ( ; )} = 𝐶 { 𝑋(4,0; 3,0)}
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

1
El valor de { 𝑋(4,0; 3,0)} se obtiene del ábaco correspondiente (Ábaco No. 2); tal como
𝐶
muestra la figura de la siguiente página:

Se puede ver en el ábaco mencionado que:

1
𝐶{ 𝑋(4,0; 3,0)} = 0,9755
𝐶

Por tanto,

𝛾̅ (𝑊; 𝑊 ′ ) = 0,5(0,9755) = 0,488

32
Cálculo de variogramas promedio Rubén Medinaceli O.

𝑳
𝟏 𝑳 𝒅
𝒅 𝑿( ; )
𝑪 𝒂 𝒂
𝑧′
𝑧

33
Cálculo de variogramas promedio Rubén Medinaceli O.

Ejercicio

Problema

Estimar el 𝛾̅ (𝑊; 𝑊 ′ ) para el caso que se muestra en la siguiente figura.

El variograma es esférico con una meseta 𝐶 = 0,5 y un rango o alcance 𝑎 = 10 𝑚

𝐿 = 40 𝑚 0m

𝑊 𝑊′

𝑑 = 30 𝑚 = 𝑙

Solución:

Tal como se puede observar, el soporte 𝑊 es una chimenea en uno de los extremos de un
block de 30 𝑚 × 40 𝑚 sobre una veta mineralizada. El soporte 𝑊 ′ representa otra chimenea
en el extremo opuesto del block.

El ábaco a ser utilizado en este caso es el que se muestra en la siguiente página (Ábaco
No. 1).

𝛾̅ (𝑊; 𝑊 ′ ) = 𝛼(𝐿; 𝑑) = 𝛼(40,30)

El valor de esta función auxiliar se obtiene de la siguiente manera. Para,

𝐿 40
= = 4,0
𝑎 10
𝑙 30
= = 3,0
𝑎 10

El ábaco de la siguiente página proporciona el siguiente valor,

1
𝛾 (𝑙) = 0,22
𝐶 𝑑

Como 𝐶 = 0,5

𝛾𝑑 (𝑙) = 0,5(0,22) = 0,11

Por otro lado, se tiene que 𝛼(0; 𝑙) = 𝐹(𝑙) = 𝐹(30)

3𝑎 1 𝑎 2 3 10 1 10 2
𝐹(30) = 𝐶 {1 − + ( ) } = 0,5 {1 − + ( ) } = 0,386
4 𝑙 5 𝑙 4 30 5 30

34
Cálculo de variogramas promedio Rubén Medinaceli O.

𝑳 𝜶(𝟎; 𝒍)
𝟏 𝟏
𝒍 𝜸 (𝒍) = 𝑪 {𝜶(𝑳; 𝒍) − 𝜶(𝟎; 𝒍)}
𝑪 𝒍
𝑀′ 𝑙

35
Cálculo de variogramas promedio Rubén Medinaceli O.

Reemplazando valores en la ecuación,

𝛾𝑑 (𝐿) = 𝛼(𝐿; 𝑑) − 𝛼(0; 𝑑)) = 𝛾𝑙 (𝐿) = 𝛼(𝐿; 𝑙) − 𝛼(0; 𝑙)

Se tiene que,

0,11 = 𝛼(𝐿; 𝑑) − 0,386

De donde,

𝛾̅ (𝑊; 𝑊 ′ ) = 𝛼(𝐿; 𝑑) = 0,11 + 0,386 = 0,496

Ejercicio

Problema:

a) Calcular 𝛾̅ (𝑊; 𝑊 ′ ) para el caso que se muestra en la siguiente figura, si el

modelo de variograma es esférico con meseta 𝐶 = 2,0 y rango o alcance 𝑎 = 10 𝑚.

b) Calcular 𝛾̅ (𝑊; 𝑊 ′ ) para el caso que se muestra en la siguiente figura, si el

modelo de variograma es exponencial con meseta 𝐶 = 2,0 y rango o alcance 𝑎 = 10 𝑚.

𝑊′

𝑎 = 2𝑚

■ 𝑊
5𝑚

𝑏 = 3𝑚

𝑐 = 3𝑚 𝑑 = 7𝑚

10𝑚

Solución:

a)

Tal como se observa, el soporte 𝑊 es una muestra puntual, un taladro de diamantina que
atraviesa el block 𝑊 ′ ubicado sobre una veta mineralizada. El block es 5 𝑚 × 10 𝑚.

36
Cálculo de variogramas promedio Rubén Medinaceli O.

En este caso, no es posible aplicar directamente una función auxiliar para calcular el
variograma promedio solicitado; sin embargo, una división conveniente del block 𝑊 ′ , que
se muestra en la figura superior, posibilita la aplicación de funciones auxiliares.

(𝑎 × 𝑐) 𝐻(𝑎; 𝑐) + (𝑎 × 𝑑) 𝐻(𝑎; 𝑑) + (𝑏 × 𝑐) 𝐻(𝑏; 𝑐) + 𝑏 × 𝑑 𝐻(𝑏; 𝑑)

𝛾̅ (𝑊; 𝑊 ′ ) =
(𝑎 × 𝑐) + (𝑎 × 𝑑) + (𝑏 × 𝑐) + (𝑏 × 𝑑)

1 𝑎 𝑐 1 2 3 1
𝐻 (𝑎; 𝑐 ) = 𝐻 (2; 3) = 𝐶 { 𝐻 ( ; )} = 𝐶 { 𝐻 ( ; )} = 2,0 { 𝐻 (0,2; 0,3)} = 2,0(0,285) = 0,570
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

1 𝑎 𝑑 1 2 7 1
𝐻 (𝑎; 𝑑 ) = 𝐻 (2; 7) = 𝐶 { 𝐻 ( ; )} = 𝐶 { 𝐻 ( ; )} = 2,0 { 𝐻 (0,2; 0,7)} = 2,0(0,520) = 1,040
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

1 𝑏 𝑐 1 3 3 1
𝐻 (𝑎; 𝑐 ) = 𝐻 (3; 3) = 𝐶 { 𝐻 ( ; )} = 𝐶 { 𝐻 ( ; )} = 2,0 { 𝐻 (0,3; 0,3)} = 2,0(0,340) = 0,680
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

1 𝑏 𝑑 1 3 7 1
𝐻 (𝑏; 𝑑 ) = 𝐻 (3; 7) = 𝐶 { 𝐻 ( ; )} = 𝐶 { 𝐻 ( ; )} = 2,0 { 𝐻 (0,3; 0,7)} = 2,0(0,550) = 1,100
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

Los valores como,

1
{ 𝐻(0,2; 0,3)} = 0,285
𝐶

han sido obtenidos del ábaco que se muestra en la siguiente figura.

37
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 3

𝑳 𝑧′ 𝑳
𝟏 𝑳 𝒅
𝒅 𝒅 𝑯 ( ; )
𝑪 𝒂 𝒂
𝑧′

𝑧 𝑧

Modelo de variograma esférico

38
Cálculo de variogramas promedio Rubén Medinaceli O.

Reemplazando valores, se tiene:

(2 × 3)(0,570) + (2 × 7)(1,040) + (3 × 3)(0,680) + (3 × 7)(1,100)

𝛾̅ (𝑊; 𝑊 ′ ) = = 0,944
(2 × 3) + (2 × 7) + (3 × 3) + (3 × 7)

b)

Al igual que en el inciso a), se tiene que,

(𝑎 × 𝑐) 𝐻(𝑎; 𝑐) + (𝑎 × 𝑑) 𝐻(𝑎; 𝑑) + (𝑏 × 𝑐) 𝐻(𝑏; 𝑐) + 𝑏 × 𝑑 𝐻(𝑏; 𝑑)

𝛾̅ (𝑊; 𝑊 ′ ) =
(𝑎 × 𝑐) + (𝑎 × 𝑑) + (𝑏 × 𝑐) + (𝑏 × 𝑑)

1 𝑎 𝑐 1 2 3 1
𝐻 (𝑎; 𝑐 ) = 𝐻 (2; 3) = 𝐶 { 𝐻 ( ; )} = 𝐶 { 𝐻 ( ; )} = 2,0 { 𝐻 (0,2; 0,3)} = 2,0(0,175) = 0,350
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

1 𝑎 𝑑 1 2 7 1
𝐻 (𝑎; 𝑑 ) = 𝐻 (2; 7) = 𝐶 { 𝐻 ( ; )} = 𝐶 { 𝐻 ( ; )} = 2,0 { 𝐻 (0,2; 0,7)} = 2,0(0,300) = 0,600
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

1 𝑏 𝑐 1 3 3 1
𝐻 (𝑎; 𝑐 ) = 𝐻 (3; 3) = 𝐶 { 𝐻 ( ; )} = 𝐶 { 𝐻 ( ; )} = 2,0 { 𝐻 (0,3; 0,3)} = 2,0(0,205) = 0,410
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

1 𝑏 𝑑 1 3 7 1
𝐻 (𝑏; 𝑑 ) = 𝐻 (3; 7) = 𝐶 { 𝐻 ( ; )} = 𝐶 { 𝐻 ( ; )} = 2,0 { 𝐻 (0,3; 0,7)} = 2,0(0,325) = 0,650
𝐶 𝑎 𝑎 𝐶 10 10 𝐶

Los valores como,

1
{ 𝐻(0,2; 0,3)} = 0,175
𝐶

han sido obtenidos del ábaco que se muestra en la siguiente figura.

39
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 7

𝑳 𝑧′ 𝑳
𝟏 𝑳 𝒅
𝒅 𝒅 𝑯 ( ; )
𝑪 𝒂 𝒂
𝑧′

𝑧 𝑧

Modelo de variograma exponencial

40
Cálculo de variogramas promedio Rubén Medinaceli O.

Reemplazando valores, se tiene:

(2 × 3)(0,350) + (2 × 7)(0,600) + (3 × 3)(0,410) + (3 × 7)(0,650)

𝛾̅ (𝑊; 𝑊 ′ ) = = 0,557
(2 × 3) + (2 × 7) + (3 × 3) + (3 × 7)

Ejercicio

Problema:

Calcular 𝛾̅ (𝑊; 𝑊 ′ ) para el caso que se muestra en la siguiente figura.

𝑊′

𝑎 = 2𝑚

5𝑚

𝑏 = 3𝑚
𝑊
■

𝑐 = 3𝑚 𝑑 = 7𝑚

10𝑚

El modelo de variograma es esférico con meseta 𝐶 = 2,0 y rango o alcance 𝑎 = 10 𝑚.

Solución:

Tal como se observa, el soporte 𝑊 es una muestra puntual, un taladro de diamantina

alejado del block 𝑊 ′ ubicado sobre una veta mineralizada. El block es 5 𝑚 × 10 𝑚.

En este caso, tampoco es posible aplicar directamente una función auxiliar para calcular el
variograma promedio solicitado; sin embargo, una composición de áreas conveniente, para
relacionar la muestra 𝑊 y el block 𝑊 ′ , como se muestra en la figura superior, posibilita la
aplicación de funciones auxiliares a la estimación del variograma promedio solicitado.

(𝑎 + 𝑏)(𝑐 + 𝑑)𝐻(𝑎 + 𝑏; 𝑐 + 𝑑) − 𝑏(𝑐 + 𝑑)𝐻(𝑏; 𝑐 + 𝑑) − 𝑐(𝑎 + 𝑏)𝐻(𝑐; 𝑎 + 𝑏) + (𝑏 × 𝑐)𝐻(𝑏, 𝑐)

𝛾̅ (𝑊; 𝑊 ′ ) =
(𝑎 + 𝑏)(𝑐 + 𝑑) − 𝑏(𝑐 + 𝑑) − 𝑐(𝑎 + 𝑏) + (𝑏 × 𝑐)

(𝑎 + 𝑏)(𝑐 + 𝑑)𝐻(𝑎 + 𝑏; 𝑐 + 𝑑) − 𝑏(𝑐 + 𝑑)𝐻(𝑏; 𝑐 + 𝑑) − 𝑐(𝑎 + 𝑏)𝐻(𝑐; 𝑎 + 𝑏) + (𝑏 × 𝑐)𝐻(𝑏, 𝑐)

𝛾̅ (𝑊; 𝑊 ′ ) =
𝑎×𝑑

41
Cálculo de variogramas promedio Rubén Medinaceli O.

1 5 10 1
𝐻 (𝑎 + 𝑏; 𝑐 + 𝑑 ) = 𝐻 (5; 10) = 𝐶 { 𝐻 ( ; )} = 2,0 { (0,5; 1,0)} = 2,0(0,730) = 1,460
𝐶 10 10 𝐶

1 3 10 1
𝐻 (𝑏; 𝑐 + 𝑑 ) = 𝐻 (3; 10) = 𝐶 { 𝐻 ( ; )} = 2,0 { (0,3; 1,0)} = 2,0(0,675) = 1,350
𝐶 10 10 𝐶

1 3 5 1
𝐻 (𝑐; 𝑎 + 𝑏) = 𝐻 (3; 5) = 𝐶 { 𝐻 ( ; )} = 2,0 { (0,3; 0,5)} = 2,0(0,445) = 0,890
𝐶 10 10 𝐶

1 3 3 1
𝐻 (𝑏; 𝑐 ) = 𝐻 (3; 3) = 𝐶 { 𝐻 ( ; )} = 2,0 { (0,3; 0,3)} = 2,0(0,340) = 0,680
𝐶 10 10 𝐶

Finalmente,

5(10)(1,460) − 3(10)(1,350) − 3(5)(0,890) + 3(3)(0,680)

𝛾̅ (𝑊; 𝑊 ′ ) = = 1,805
2(7)

Ejercicio

Problema:

Calcular 𝛾̅ (𝑊; 𝑊 ′ ) para el caso que se muestra en la siguiente figura:

𝑊 𝑊′ 30 𝑚

20 𝑚 20 𝑚 40 𝑚

El modelo de variograma es esférico con una meseta 𝐶 = 2,0 y un rango o alcance 𝑎 =

40 𝑚.

Solución:

El soporte 𝑊 es un block de 30 𝑚 × 20 𝑚 y el soporte 𝑊 ′ es otro block de 30 𝑚 × 40 𝑚. Los

blocks están separados entre si 20 𝑚.

No existe una función auxiliar directa para estimar el variograma promedio solicitado. Para
posibilitar la estimación del variograma promedio inicialmente se discretiza el block 𝑊 y
luego se aplican las funciones auxiliares correspondientes, tal como muestra la siguiente
figura:

42
Cálculo de variogramas promedio Rubén Medinaceli O.

ℎ2 ℎ1 𝑊′ 30 𝑚

10𝑚 10𝑚 20 𝑚 40 𝑚

El block 𝑊 ha sido discretizado utilizando las rectas ℎ1 y ℎ2 ; en otras palabras, estas dos
rectas representan al block 𝑊. Por tanto,

(10 × 30) 𝛾̅ (ℎ1 ; 𝑊 ′ ) + (10 × 30) 𝛾̅ (ℎ2 ; 𝑊 ′ ) 𝛾̅ (ℎ1 ; 𝑊 ′ ) + 𝛾̅ (ℎ2 ; 𝑊 ′ )

𝛾̅ (𝑊; 𝑊 ′ ) = =
(10 × 30) + (10 × 30) 2

(65 × 30)𝑋(65; 30) − (25 × 30)𝑋(25; 30)

𝛾̅ (ℎ1 ; 𝑊 ′ ) =
(65 × 30) − (25 × 30)

1 65 30 1
𝑋(65; 30) = 𝐶 { 𝑋 ( ; )} = 𝐶 { 𝑋(1,625; 0,75)} = 2,0(0,820) = 1,640
𝐶 40 40 𝐶

1 25 30 1
𝑋(25; 30) = 𝐶 { 𝑋 ( ; )} = 𝐶 { 𝑋(0,625; 0,75)} = 2,0(0,600) = 1,200
𝐶 40 40 𝐶

1950(1,640) − 750(1,200)
𝛾̅ (ℎ1 ; 𝑊 ′ ) = = 1,915
1200

(75 × 30)𝑋(75; 30) − (35 × 30)𝑋(35; 30)

𝛾̅ (ℎ2 ; 𝑊 ′ ) =
(75 × 30) − (35 × 30)

1 75 30 1
𝑋(75; 30) = 𝐶 { 𝑋 ( ; )} = 𝐶 { 𝑋(1,875; 0,75)} = 2,0(0,855) = 1,710
𝐶 40 40 𝐶

1 35 30 1
𝑋(35; 30) = 𝐶 { 𝑋 ( ; )} = 𝐶 { 𝑋(0,875; 0,75)} = 2,0(0,690) = 1,380
𝐶 40 40 𝐶

2250(1,710) − 1050(1,380)
𝛾̅ (ℎ2 ; 𝑊 ′ ) = = 2,000
1200

Finalmente,

𝛾̅ (ℎ1 ; 𝑊 ′ ) + 𝛾̅ (ℎ2 ; 𝑊 ′ ) 1,915 + 2,000

𝛾̅ (𝑊; 𝑊 ′ ) = = = 1,957
2 2

43
Cálculo de variogramas promedio Rubén Medinaceli O.

2.3. Funciones auxiliares en tres dimensiones

𝑷 es un paralelepípedo con vértices 𝐴, 𝐵, 𝐶, 𝐷, 𝐸, 𝐹, 𝐺, 𝐻, tal como se muestra en la siguiente

figura:

𝐵 𝒅 𝐷

𝑳 𝒅
𝑷
𝐹 𝐻

𝐴 𝐶

𝐸 𝐺
Nótese que,

̅̅̅̅
𝐴𝐶 = ̅̅̅̅
𝐴𝐸 = ̅̅̅̅
𝐶𝐺 = ̅̅̅̅
𝐸𝐺 = 𝑑

̅̅̅̅ = 𝐵𝐹
𝐵𝐷 ̅̅̅̅ = 𝐷𝐻
̅̅̅̅ = 𝐹𝐻
̅̅̅̅ = 𝑑

̅̅̅̅
𝐴𝐵 = ̅̅̅̅
𝐶𝐷 = ̅̅̅̅
𝐸𝐹 = ̅̅̅̅
𝐺𝐻 = 𝐿

Es posible definir cuatro funciones auxiliares.

2.3.1. Función auxiliar 𝜶

𝛼(𝐿; 𝑑 2 ) = 𝛾̅ (𝐴𝐶𝐸𝐺; 𝐵𝐷𝐹𝐻)

2.3.2. Función auxiliar 𝑿

𝑋(𝐿; 𝑑 2 ) = 𝛾̅ (𝐴𝐶𝐸𝐺; 𝑃)

𝑋(𝐿; 𝑑 2 ) = 𝛾̅ (𝐵𝐷𝐹𝐻; 𝑃)

2.3.3. Función auxiliar 𝑯

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐴𝐶𝐸𝐺; 𝐴𝐶𝐵𝐷)

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐴𝐶𝐸𝐺; 𝐴𝐸𝐵𝐹)

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐴𝐶𝐸𝐺; 𝐸𝐺𝐹𝐻)

44
Cálculo de variogramas promedio Rubén Medinaceli O.

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐴𝐶𝐸𝐺; 𝐶𝐺𝐷𝐻)

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐵𝐷𝐹𝐻; 𝐵𝐷𝐴𝐶)

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐵𝐷𝐹𝐻; 𝐵𝐹𝐴𝐸)

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐵𝐷𝐹𝐻; 𝐹𝐻𝐸𝐺)

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐵𝐷𝐹𝐻; 𝐷𝐻𝐶𝐺)

̅̅̅̅ ; 𝑃)
𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐴𝐶

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (̅̅̅̅
𝐴𝐸 ; 𝑃)

̅̅̅̅ ; 𝑃)
𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐶𝐺

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (̅̅̅̅
𝐸𝐺 ; 𝑃)

̅̅̅̅ ; 𝑃)
𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐵𝐷

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (̅̅̅̅
𝐵𝐹 ; 𝑃)

̅̅̅̅ ; 𝑃)
𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (𝐷𝐻

𝐻(𝐿; 𝑑 2 ) = 𝛾̅ (̅̅̅̅
𝐹𝐻 ; 𝑃)

2.3.4. Función auxiliar F

𝐹(𝐿; 𝑑 2 ) = 𝛾̅ (𝑃; 𝑃)

2.3.5. Estimación de las funciones auxiliares en tres dimensiones

Los valores de las funciones auxiliares en tres dimensiones se obtienen de los ábacos que
se muestran en las siguientes páginas. Se incluyen ábacos para el modelo de variograma
esférico y para el modelo de variograma exponencial. La fuente de dichos ábacos es el libro
Mining Geostatistics de A.G. Journel y Ch.J. Huijbregts, publicado por Academic Press en
1978.

45
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 9

𝑳
𝟏 𝑳 𝒅
𝑭(𝑳; 𝒅𝟐 ) 𝑭 (𝒂 ; 𝒂)
𝑪
𝒅

Modelo de variograma esférico

46
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 10

𝑳
𝟏
{𝜶(𝑳; 𝒍𝟐 ) − 𝜶(𝟎, 𝒍𝟐 )}
𝑪
𝒅

Modelo de variograma esférico

47
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 11

𝑳
𝟏 𝑳 𝒅
𝑭(𝑳; 𝒅𝟐 ) 𝑭 (𝒂 ; 𝒂)
𝑪
𝒅

Modelo de variograma exponencial

48
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 12

𝑳
𝟏
{𝜶(𝑳; 𝒍𝟐 ) − 𝜶(𝟎, 𝒍𝟐 )}
𝑪
𝒅

Modelo de variograma exponencial

49
Cálculo de variogramas promedio Rubén Medinaceli O.

Ejercicio

Problema:

Se tiene un block 𝑊 de 20 𝑚 × 25 𝑚 × 20 𝑚 𝑑𝑒 𝑎𝑙𝑡𝑢𝑟𝑎 en un depósito masivo de 𝐶𝑢. El

variograma es esférico con una meseta 𝐶 = 2,0 (%)2 y un rango o alcance 𝑎 = 20𝑚.

Obtener 𝛾̅ (𝑊; 𝑊)

𝑑 = 20 𝑚

𝐿 = 25 𝑚

𝑑 = 20 𝑚

Solución:

Una primera solución puede obtenerse discretizando el block 𝑊, tal como muestra
la siguiente figura; en este caso,
8 8
1
𝛾̅ (𝑊; 𝑊 ) = ∑ ∑ 𝛾(𝑧𝑖 ; 𝑧𝑗 )
8×8
𝑖=1 𝑗=1

Inicialmente, se calculan las distancias entre los puntos discretos; y posteriormente,

los valores del variograma para las diferentes distancias obtenidas. Sin embargo,
nótese que la doble sumatoria incluye 64 términos que, son reiterativos y sencillos,
calcularlos manualmente tomaría mucho tiempo.

Es importante reiterar que cualquier software que calcula variogramas promedio,

utiliza la discretización.

50
Cálculo de variogramas promedio Rubén Medinaceli O.

𝑧3 𝑧4

𝑧1 𝑧2

𝑧7 𝑧8

𝑧5 𝑧6

Una alternativa, es recurrir a la función auxiliar tridimensional 𝐹(𝐿; 𝑑 2 ); esto es,

1 𝐿 𝑑 1 25 20 1
𝛾̅ (𝑊; 𝑊) = 𝐹(𝐿; 𝑑 2 ) = 𝐶 { 𝐹 ( ; )} = 𝐶 { 𝐹 ( ; )} = 𝐶 { 𝐹(1,25; 1,0)}
𝐶 𝑎 𝑎 𝐶 20 20 𝐶
1
Para obtener 𝐹(1,25; 1,0) se recurre al ábaco que muestra la siguiente página. Donde,
𝐶

1
{𝐹(1,25; 1,0)} = 0,820
𝐶

Consecuentemente,

𝛾̅ (𝑊; 𝑊) = 2,0(0,820) = 1,640

51
Cálculo de variogramas promedio Rubén Medinaceli O.

ÁBACO No. 9

𝑳
𝟏 𝑳 𝒅
𝑭(𝑳; 𝒅𝟐 ) 𝑭 ( ; )
𝑪 𝒂 𝒂
𝒅

Modelo de variograma esférico

52