礦產資源估算中定量鄰域分析與優化研究

劉金輝，王京彬，李　峰，張漢成，王　翔

(1.北京礦產地質研究院，北京 100012；2.昆明理工大學，云南 昆明650500；3.有色金屬礦產地質調查中心，北京100012；4.中國地質調查局發展研究中心，北京 100037)

礦產資源估算中定量鄰域分析與優化研究

劉金輝1，2，王京彬1，李峰2，張漢成3，王翔4

(1.北京礦產地質研究院，北京 100012；2.昆明理工大學，云南 昆明650500；3.有色金屬礦產地質調查中心，北京100012；4.中國地質調查局發展研究中心，北京 100037)

摘要：在應用線性、非線性地質統計學進行礦產資源/儲量估算或品位條件模擬過程時，鄰域分析是其中非常關鍵和重要的一個環節。通常情況下，在進行資源/儲量估算中鄰域參數通過人為控制和設定，這導致在一定程度上具有主觀因素。研究認為估計鄰域的大小，應由各個方向變異函數變程及搜索橢球體等參數來確定，優化鄰域需選擇足夠多的樣本，且得到最大回歸斜率，同時具有最小負克里金權重?？死锝疣徲騾档暮侠碓O置將直接影響最終塊體估值結果，這為資源量估算和礦床評價提供了理論和實踐基礎。

關鍵詞：地質統計學；回歸斜率；克里金方差；鄰域優化

線性和非線性地質統計學方法，在資源估算和采礦品位控制中的應用已經得到普遍認可，在此過程中如何確定搜索鄰域對于估算結果的質量十分重要。本文以某韌性剪切帶型金礦5#礦體為研究對象，詳細討論了定量鄰域分析與優化對于資源估算結果的重要影響。

1鄰域分析與優化

在進行資源估算或條件模擬之前，要進行搜索鄰域的設置。設置鄰域的參數一般有搜索距離、搜索樣品數、搜索橢球體參數、塊體模型大小等。優化的搜索鄰域是指有足夠多的樣本選擇，且得到最大的回歸斜率并具有最小負克里金權重，其意義是在某個估值點上擁有最大數量的可用的樣品數，且所有的負權重不會超過某個比例。搜索鄰域優化的目的就是指達到樣本域內的局部最大無偏估計，使樣本實際值具有最小的可能誤差，同時條件偏倚達到最小[1-12]。

一般情況下，估計鄰域越大，克里金估計精度越高。考慮樣品的優選級克里金估值特點，其中估計方差、真值與估計值線性表達式的斜率、真值與估計值的相關系數以及克里金估值中負權數目4個結果數值作為估計鄰域優化參數[13-16]。實踐中克里金鄰域優化主要考慮兩個參數：回歸斜率與克里金方差。

1.1回歸斜率

在變異函數有效且回歸為線性條件下，即能計算估計值與真值回歸斜率的主要參數。實際估算中由于不知道真值，但是可以計算真值與估計值之間的協方差進行衡量估計值的精度。如式(1)所示。

(1)

式(1)可由塊體估計方差，真值與估計值的協方差公式所得。

其中塊體品位估計值見式(2)

(2)

塊體的估計方差可以由式(3)推出。

(3)

在這里為樣品位置與和克里金權重相關的表示為和的協方差函數值。由式(4)求出。

(4)

為每一個樣品位置與將要估計的塊V之間的協方差的平均值(式(5)、式(6))

(5)

(6)

式中：μ為拉格朗日乘數，用在解決克里金導入協方差函數。理想狀態下，a的值應該很接近1.0，意味著條件無偏。在這樣的情況下，塊體的真值就接近于估計值。

1.2克里金方差

普通克里金方程組如式(7)所示。由式(7)得出式(8)。將式(8)代入估計方差公式，則引出式(9)。

(7)

(8)

(9)

(10)

克里金方差給出了一個相對的關于數據存儲密度與幾何分布估計質量標準，但是回歸斜率更有用途一些[17]。

2定量鄰域分析

2.1定量鄰域分析步驟

定量克里金鄰域分析主要采用以下步驟：①計算“真”品位與“估計”品位的回歸斜率；②鄰域中的樣品克里金權重分布；③克里金方差。

2.2測試塊克里金鄰域分析

搜索鄰域設置為一個橢球體，該橢球體三個軸的方向與變異函數的方向一致。最大主軸方向為礦體的礦化最連續方向，最短主軸方向為垂直礦化體的厚度方向。在5#礦體中選擇一個測試塊段，估計鄰域結果。5#礦體選擇塊段中心坐標為(X=436455.00，Y=497875.00，Z=309.00)，最大搜索鄰域設置為60m，表測試結果如表1所示。表1中顯示并不是搜索的樣品數量越多，估計精度就越高，當OK32時，回歸斜率最接近于1，并且出現最小的負權重，當大于這個值時候，克里金方差基本保持不變，回歸斜率保持不變，但是負權重比例增加很快，這種測試結果表明在OK32附近斜率出現開始趨于平穩的規律。

圖1、圖2分別為5#礦體中OK32鄰域中樣品距離與所屬權重分布立體圖與權重與距離之間的關系圖。從圖1和圖2分析可知，權重與樣品距離有一定關系，但并不是距離越近，權重就一定越高，并沒有形成隨著距離增加，近距離樣品貢獻越大，而遠距離樣品貢獻很小，因此此時的估計鄰域最為理想。

表1　5#礦體搜索鄰域測試結果

3鄰域大小與鄰域優化關系

鄰域大小與克里金估計方差，均值、真值與斜率，真值與估計值的相關系數，以及克里金估值中負權重樣品數等之間的關系，表面上顯示估計鄰域越大，估計方差、均值權重及克里金估值中的負權數目越小，真值與估計值的回歸斜率以及真值與估計值的相關系數就越大，當估計鄰域大到一定程度，這些參數都穩定在一定數值上，此時的估計鄰域最為理想。通過對鄰域分析結果研究表明，鄰域大小并不是決定鄰域是否優化的關鍵因素，對于鄰域相關參數的影響主要是由于不同鄰域大小會選擇不同量的樣品數。一般在進行鄰域搜索之前已經確定了最大樣品數，當估計鄰域小的時候，則達不到所設定的最大樣品數，隨著估計鄰域增大，樣品數增多，則估計鄰域的各個參數如估計方差和回歸斜率等趨于穩定，當達到最大樣品數時，各個參數值均位于穩定值范圍內，即使搜索鄰域再擴大，鄰域各個參數亦穩定在某范圍內。

進行品位估值之前，估計鄰域的大小可以由各個方向變異函數的變程及搜索橢球體參數等確定。由此可以得出結論，即優化的搜索鄰域是指要在某個估值點上擁有最大數量的可利用樣品，最少的負權和條件偏倚達到最小。

4克里金鄰域分析對估值結果的影響

測試點結果顯示，樣品品位估計值與鄰域有著重要的關系，不同的鄰域樣品數所產生的估值結果不盡相同，變化也很大。如圖3所示，估計鄰域對于測試點的估計值有著重要影響，樣品數量在80個之前，估計值呈現非線性遞增現象，當樣品數達到一定程度后，再增加樣品數則呈現估計值稍微遞減的現象。

圖4和圖5分別為5#礦體不同鄰域分析下的邊界品位-平均品位圖和品位-噸位曲線圖。綜合研究圖3～5可知，塊體中品位估計值直接影響著整個礦化域的資源量估算，如果不選擇優化鄰域，則可能對估算結果造成較大影響，尤其對于貴金屬礦產，品位和噸位直接關系礦床的評估價值。

圖5　5#礦化域品位-噸位模型

5結論

鄰域優化為礦床資源量估算，以及相應礦床評價提供了相對精確的理論依據。通過研究5#礦體相關領域參數，顯示鄰域優化對于克里金估值具有重要影響，主要結論如下所示。

1) “真”品位與“估計”品位的回歸斜率、克里金權重和克里金方差是評價塊體估值精確度的三個基本要素，估計鄰域的大小由各個方向變異函數的變程及搜索橢球體參數等確定。

2) 通過定量研究克里金鄰域分析，確定了優化鄰域是指選擇足夠多的樣本，且得到最大回歸斜率，同時具有最小負克里金權重。

3) 噸位-品位曲線模型和邊界品位-平均品位模型分析說明，克里金鄰域參數的合理設置直接影響最終塊體估值結果，這對資源量估算和礦床評價十分重要。

參考文獻

[1]Goovaerts P.Geostatistics for Natural Resources Evaluation [M].New York：Oxford University Press，1997.

[2]王正帥，顧和和.數字化礦山資源儲量計算方法[J].煤田地質

與勘探，2010，38(3)：6-9.

[3]Diamond，P.and M.Armstrong，Robustness of Variograms and conditioning of krigingrices[J].Mathematical Geology，1984，16(8)：809-822.

[4]陳伯茂.地質統計學在礦產含量分級中的應用[J].地質與勘探，1986.22(2)：24-30.

[5]Srivastava，R.M.Minimum variance or maximum profitability[J].CIM Bulletin，1987，80(901)：63-68.

[6]Srivastava，R.M.The visualization of spatial uncertainty.Volume 3 of AAPG Computer Applications in Geology，1994，339-345.

[7]Thayer，William C.，Daniel A.Griffith，Philip E.Goodrum，Gary L.Diamond，and James M.Hassett.Application of Geostatistics to Risk assessment [J].Risk Analysis，2003，23(5):945-960.

[8]馬洪濱，熊俊楠.基于地質統計學的儲量估算系統[J].煤田學報，2007，32(3)：267-271.

[9]Goovaerts，P.Geostatistics for Natural Resource Evaluation [M].New York：Oxford University Press，1997:167.

[10]向永生，孔愛云．克里金估計鄰域大小的確定方法[J]．黃金地質，1998，4(3) ：76-79．

[11]孫進，張佳榮，侯斌.礦產儲量計算經典統計與地質統計學方法的對比分析[J].采礦技術，2005，5(2)：80-81.

[12]A.G.Geostatistics for conditional simulation of orebodies[J].Economic Geology，1974，69，673-687.

[13]郭奇峰，劉金輝．析取克里金法在可回采資源量的總體估算的應用[J]．中國礦業，2012，21(8)：20-23．

[14]侯景儒，黃競先，等譯．礦業地質統計學[M]．北京：冶金工業出版社，1982．

[15]王仁鐸，胡光道．線性地質統計學[M]．北京：地質出版社，1988．

[16]侯景儒，尹鎮南，李維明．實用地質統計學[M]．北京：地質出版社，1998．

[17]孫洪泉．地質統計學及其應用[M]．徐州：中國礦業大學出版社，1990．

Quantified kriging neighborhood analysis and optimization in mineral resource estimation

LIU Jin-hui1，2，WANG Jing-bin1，LI Feng2，ZHANG Han-cheng3，WANG Xiang4

(1.Beijing Institute of Geology for Mineral Resources，Beijing 100012，China;2.Kunming University of Science and Technology，Kunming 650500，China;3.China Non-Ferrous Metals Resource Geological Survey，Beijing 100012，China;4.Development and Research Center，China Geological Survey，Beijing 100037，China)

Abstract:During the mineral Resources/Reserves estimation by linear，non-linear geostatistics as well as mineral grade conditional simulation，neighborhood analysis is of the very critical and important process.However，the neighborhood parameters for the mineral Resources/Reserves estimation are commonlyto be defined manually，which has subjective factorsin part and might lead to the estimated value mismatch with actual value，once the neighborhood parameters are to be set arbitrarily.The searching distance of neighborhood should be defined by the directional variogram and search ellipsoid，and the optimized neighborhood is to search sufficient samples，and get maximum regression of slop with the minimum negative kriging weight.The parameters of kriging neighborhood directly affect the final block estimated value，therefore，selecting reasonable and optimized neighborhood parameters will provide the theoretical and practical basis for the mineral resource estimation and mineral deposit valuation.

Key words：geostatistics；regression of slope； Kriging variance； neighborhood optimization

收稿日期：2015-04-25

作者簡介：劉金輝(1978-)，男，在站博士后，綜合信息成礦預測專業，主要從事礦產資源評價研究。E-mail:jhliu0922@qq.com。

中圖分類號：TD15

文獻標識碼：A

文章編號：1004-4051(2016)03-0154-04