Krigagem Ordinária Aplicada à Densidade: Fundamentos, Variograma e Validação
A krigagem ordinária (Ordinary Kriging – OK) é o método geoestatístico mais robusto e amplamente aceito para estimar densidade aparente (bulk density) em modelos de recursos minerais. Ela pertence à família dos estimadores lineares não-viesados com variância mínima (BLUE) e foi desenvolvida por Georges Matheron a partir dos trabalhos pioneiros de D.G. Krige. Quando aplicada à densidade, a OK não estima apenas um número — ela estima a variável que converte volume em tonelagem, tornando-a crítica para o cálculo de metal contido, planejamento de frota, reconciliação e viabilidade econômica.
Diferente da estimativa de teores (que são composicionais e frequentemente lognormal), a densidade é uma variável aditiva em volume, mas com características próprias: apresenta continuidade espacial forte dentro de domínios geológicos homogêneos, mas alta variabilidade entre domínios (ex.: saprolito vs. rocha fresca). Por isso, a krigagem de densidade exige domínios bem delimitados — exatamente como se faz com teores. Fundamentos aplicados à densidade A equação da krigagem ordinária para estimar a densidade 𝑍∗(𝑥0) em um bloco ou ponto não amostrado é:
𝑍∗(𝑥0) = ∑𝜆𝑖 𝑖=1 𝑛 𝑍(𝑥𝑖)com∑𝜆𝑖 𝑖=1 =1
Os pesos 𝜆𝑖 são obtidos resolvendo o sistema de krigagem, que minimiza a variância do erro de estimativa e garante não-viesamento. Na prática, isso significa que a densidade estimada em um bloco de 10×10×10 m será ponderada pelas amostras vizinhas, respeitando a continuidade espacial real do depósito. O resultado direto é a tonelagem do bloco:
Tonelagem = Volume do bloco × 𝑍∗(𝑥0)
Essa multiplicação simples explica por que um erro de apenas 0,05–0,10 g/cm³ na densidade pode gerar diferenças de milhões de toneladas no inventário de recursos.
Variograma aplicado à densidade O variograma é o coração da krigagem de densidade. Ele quantifica como a diferença de densidade aumenta com a distância:
y(ℎ) = 1 2𝑁(ℎ) 𝑁(ℎ) ∑[𝑍( 𝑖=1 𝑥𝑖) −𝑍(𝑥𝑖 +ℎ)]2
Em dados reais de densidade, observa-se frequentemente:
– Efeito pepita (nugget) moderado a alto em zonas supergênicas, saprolitos ou materiais fraturados (devido à micro-variabilidade e erro de amostragem).
– Modelo esférico ou exponencial como os mais adequados.
– Anisotropia forte: maior continuidade ao longo do mergulho ou da foliação, menor na direção perpendicular.
O alcance típico varia de 40–150 m em depósitos de ouro ou cobre e pode chegar a 300 m em minério de ferro bandeado. O variograma deve ser modelado separadamente por domínio (litologia ou geometalúrgico). Krigar densidade sem separar domínios é um erro grave — gera suavização excessiva e viés na tonelagem.
Validação específica para densidade
Antes de usar o modelo, a validação é obrigatória:
1. Validação cruzada (leave-one-out): remove-se cada amostra, estima-se sua densidade e compara-se com o valor real. Métricas principais: erro médio próximo de zero, variância normalizada próxima de 1 e correlação > 0,75–0,85.
2. Swath plots: comparação real × estimado ao longo de faixas (N-S, L-E e vertical). Ideal para densidade porque revela viés sistemático em zonas de transição litológica.
3. Kriging Efficiency (KE) e Slope of Regression (SLOR): valores próximos de 1 indicam boa performance e baixa suavização.
Discussão científica: é válido krigar densidade?
Existe um debate recorrente na literatura geoestatística sobre a possibilidade de krigar densidade diretamente. Alguns autores argumentam que:
– Densidade não é estritamente aditiva como o teor (a média ponderada por volume não é sempre linear).
– Em depósitos com forte variação de porosidade ou umidade, a estacionariedade pode ser violada.
– Há risco de “smearing” (suavização) que distorce o cálculo de tonelagem em blocos pequenos.
No entanto, a prática consolidada (JORC Code, NI 43-101, SAMREC) e dezenas de estudos publicados (ex.: Lipton, 2001; Abzalov, 2016; Nwaila et al., 2025; Rossi & Deutsch, 2014) demonstram que a krigagem ordinária de densidade é não apenas válida, mas recomendada, desde que:
– Seja feita dentro de domínios homogêneos bem delimitados.
– O variograma seja robusto e anisotrópico.
– A validação cruzada confirme baixa viés e boa reprodutibilidade.
– Seja usada a densidade aparente (bulk) correta para cada domínio (seca, úmida ou in situ).
Quando essas condições são atendidas, a OK supera métodos simples (média por domínio ou IDW) e reduz drasticamente o erro na tonelagem final. A alternativa de não krigar densidade (usar média constante por domínio) só é aceitável em depósitos muito homogêneos ou em fases iniciais de exploração.
Krigar densidade não é apenas tecnicamente correto — é uma responsabilidade técnica e econômica. Um modelo de densidade bem krigado transforma um volume geológico em tonelagem confiável, impactando diretamente frota, pilhas, reconciliação, carga de caminhões/navios e, consequentemente, o valor do projeto. Quando bem executada, a krigagem ordinária de densidade garante que o número que sai do software reflita a realidade do maciço, e não apenas uma média arbitrária.
É por isso que, nos relatórios de recursos de classe Measured e Indicated, a krigagem de densidade com validação rigorosa é hoje padrão de mercado.
