克里格法在離子吸附型稀土礦勘查儲量估算中的應用

趙 汀，王登紅，王瑞江，鄧茂春，陳為光

(1.中國地質科學院礦產資源研究所，國土資源部成礦作用與資源評價重點實驗室，北京 100037；2.江西省地質礦產勘查開發局贛南地質調查大隊，江西 贛州 341000)

離子吸附型稀土礦是我國獨特的優勢礦產資源，近年來取得了很大的找礦進展，但資源消耗也是很快的[1-2]。如何客觀評價我國離子吸附型稀土資源的儲量，乃當務之急。在我國地質礦產勘查工作中，傳統的儲量估算方法有：算術平均法、地質塊段法、多角形法、斷面法(包括垂直剖面法和水平斷面法)及距離平方反比法等[3]。這些傳統方法在我國多年的地質找礦工作中已經成為國家規范，絕大部分礦產勘查報告儲量是用傳統方法估算的，但存在的不足之處在于：①把計算儲量的主要參數(如：品位)視成固定變量，忽略了這類變量的隨機性；②在變量的空間相關關系上也僅僅考慮了距離，忽略了變量的地質變異規律性；③無法給出估計精度。克里格插值法的區域化變量理論及其變差函數工具，正好有效地彌補了傳統的儲量估算方法的缺點。

1 克里格插值法的基本原理

克里格插值法是法國馬特隆教授以南非礦山地質工程師D. G .克里格的名字命名的一種方法[4]。它是以區域化變量理論為基礎，以變差函數為主要工具，在保證估計值滿足無偏性條件和最小方差條件的前提下求得估計值。設區域化變量f(x)滿足二階平穩假設或本征假設，則待插點P的估計值為：