對Surpac軟件與地質統計學的礦床資源量估算結合應用的探究

■陳媛娟光有明

（1西寧特殊鋼股份有限公司有限公司青海西寧 810005；2青海請齊鑫地質礦產勘查股份有限公司青海西寧 810016）

對Surpac軟件與地質統計學的礦床資源量估算結合應用的探究

■陳媛娟1光有明2

（1西寧特殊鋼股份有限公司有限公司青海西寧 810005；2青海請齊鑫地質礦產勘查股份有限公司青海西寧 810016）

利用Surpac軟件將某銅礦床的鉆探工程數據進行處理，建立其地質數據庫，銅礦體的三維模型及品位變異函數模型，結合塊體估值技術選擇合理的地質統計學方法估算其資源量，并與距離冪次反比法估算結果進行對比研究，分析銅礦體不同方向上的品位變化趨勢。結果表明，基于Surpac軟件與地質統計學方法對該礦床進行資源量估算，結果可靠，操作簡便，工作效率高，有利于推動地質統計學在國內礦產資源評價中的應用。

Surpac軟件 地質統計學 礦床資源量 估算 結合

1　地質數據庫

地質數據庫是管理地質數據信息的數據庫，其記錄勘探工程的空間位置、樣品化驗數據、巖性等信息。強大的數據庫管理功能是Surpac軟件的一大特色，用戶可根據需要對數據進行一系列處理，如查看、編輯、基礎統計分析等。本文依據Surpac軟件格式要求，建立了所研究礦床的地質數據庫，結構見表1。

表1　地質數據庫結構

2　樣品組合與統計分析

區域化變量具有相同的承載是得到參數的無偏估計量的前提。樣品組合即將空間不等長的樣長和品位，量化到一些離散點。在地質數據庫中提取全部樣品長度，基礎統計分析窗口。絕大多數樣品長度為2.00m左右，故選擇2.00m為組合樣長度，采用長度加權法組合樣品。

確定組合樣長度后，在Surpac中以2.00m對Cu品位進行組合，并對組合后的Cu品位進行基礎統計分析，得到其直方分布圖。Cu品位近似服從正太分布，故選擇普通克立格法對礦體品位進行插值。該礦床屬于斑巖型銅礦床，Cu品位分布均勻，變化系數為0.58，特高品位的樣品較少，無需處理。

3　塊體模型與品位插值

3.1塊體模型

Surpac軟件將塊體模型與地質統計學相結合，根據構建的變異函數模型對其他塊體進行品位插值。對地質數據庫中的Cu品位進行解譯，以銅礦床工業指標一般要求圈定礦體（見圖1）。以礦體模型為約束條件，在Surpac中創建塊體模型，塊尺寸為10m×10m× 5m，次級分塊尺寸為5m×5m×2.5m。塊體模型的功能與數據庫相似，記載礦體的品位、巖礦類型、比重等諸多屬性，可根據需要進行相關操作。

3.2品位插值

根據地質統計學原理，距離待估塊體越近的已知點對待估塊體的影響越大，即插值時所分配的權系數越高。在插值過程中，Surpac軟件根據變異函數模型與各向異性橢球體參數分配已知點對待估塊體的權系數。通過分析，本文選擇普通克立格法對每個塊體進行品位插值，并將估值結果與距離冪次反比法估值結果進行對比分析。通過分析可得，兩種估算方法的Cu品位幾乎相等，金屬量的差異也非常微小，相對誤差在1%內。這表明了普通克立格法估值的準確性，也驗證了變異函數模型的合理性。

圖1　礦體模型

4　校驗及評估

任何一種估值都不可能與實際情況完全相符，運用普通克立格法對礦體品位進行估值亦然。因此，對估值結果進行校驗及評估異常重要。其目的是再次檢驗所構建的變異函數模型是否合理，同時驗證選擇的估值方法是否得當。本次工作選擇品位-噸位曲線和品位趨勢分析對估值結果進行評估。品位-噸位曲線從Surpac軟件中，可直接導出礦體的資源量報告，根據報告作出礦體的品位-噸位曲線（見圖2）。從圖2可以看出，普通克立格法與距離冪次反比法估值結果異常相近，僅在高品位區間（Cu品位大于0.60%）存在微弱差異。兩種估值方法Cu品位在0～0.20%和0.60%～1.40%的金屬量很少，大部分金屬量品位落在Cu品位0.20%～0.60%之間。本例是一個斑巖型銅礦床，品位分布較均勻。因而，這也與礦體的地質意義相吻合。

圖2　品味—噸位曲線

5　結語

本文分析并處理了某銅礦床的鉆探工程數據，利用Surpac軟件建立其三維地質模型（地質數據庫、礦體模型等），對Cu品位進行了統計學分析，構建了其變異函數模型，選擇普通克立格法對該礦床進行品位插值，估算其資源量，將估算結果與距離冪次反比法估算結果進行對比研究，從側面證實了資源量估算工作的可靠性。同時，本文通過品位-噸位曲線揭示了兩種資源量估算方法得出的結果相吻合，且符合礦床本身的地質意義；通過品位變化趨勢分析揭示了兩種方法估值后，礦體Cu品位在北、東兩個方向的變化趨勢保持高度一致且與組合樣品位變化趨勢保持較好的協調性，進一步校驗了估算結果的可靠性。綜上所述，基于Surpac軟件與地質統計學方法對該礦床進行資源量估算，結果可靠，操作簡便，工作效率高，有利于推動地質統計學在國內礦產資源評價中的應用。

P208[文獻碼]B

1000-405X（2016）-7-269-1