劉文凱，張雄天

(蘭州有色冶金設計研究院有限公司, 甘肅 蘭州 730000)

0 引言

Micromine礦業(yè)軟件是一款大型三維礦業(yè)軟件，該軟件廣泛運用于礦山地質勘探、資源量估算、礦山開采設計等。Micromine礦業(yè)軟件通過對礦體進行三維地質建模，實現了礦山數字化、可視化，對于礦山設計及促進礦山數字化建設具有十分重要的意義。

本文利用Micromine軟件對某銀多金屬礦進行三維地質建模及其資源量估算。

1 地質概況

某銀多金屬礦大地構造位置屬于內蒙古地槽褶皺系（Ⅱ）溫都爾廟-翁牛特旗加里東地槽褶皺帶。礦區(qū)內出露地層有侏羅系上統(tǒng)滿克頭鄂博組（J3mk）、第四系全新統(tǒng)（Qh）。礦區(qū)內主要構造以北東-南西向為主，構造類型以斷裂構造為主，褶皺微弱。礦體賦存于構造裂隙中，與北西西向（310°～330°）節(jié)理裂隙密切相關。礦床屬于與中生代火山熱液有關的淺成低溫熱液礦床。

礦床內共圈定11條礦體，礦區(qū)規(guī)模最大的為I號礦體，其呈陡傾斜的大脈狀產于地表136 m以下，是礦區(qū)最主要的礦體。礦體東西長634 m，傾向延深673 m，賦存標高1342～720 m，沿走向波折延伸，傾向上亦呈舒緩、波狀產出。礦體走向109°～165°，傾向南199°～255°，傾角為65°～89°，平均為77°。礦體厚度為0.1～11.26 m，平均厚度為2.52 m。I號礦體平均品位Ag 330.06 g/t、Cu 0.31%、Pb 7.07%、Zn 6.80%。

I號礦體按礦石類型可分為兩個礦層：中心為 大脈狀銀鉛鋅硫化物塊狀礦層，走向、傾向上連續(xù)性較好，產狀穩(wěn)定，有用組分含量較高，編號為I-1號礦體，礦石量占I號礦體的22.47%，金屬量占73.83%。緊鄰該富礦層頂、底板產出有厚度不等的浸染狀銀鉛鋅硫化物貧礦層，編號為I-2號礦體，是近礦圍巖具不同程度的礦化而形成。

礦石中主要有用元素為銀、鉛、鋅、銅。礦石結構簡單，塊狀礦石主要為半自形—它形晶粒狀結構、交代結構，浸染狀礦石主要為它形晶粒狀結構、交代侵蝕結構。

礦石構造主要有塊狀構造、稀疏浸染構造、脈狀、網脈狀構造等。

礦石自然類型均為原生礦石，礦石工業(yè)類型為含銅的銀鉛鋅硫化礦石。

2 鉆孔數據庫建立

鉆孔數據庫是三維地質建模的基礎，主要由4個Excel電子表格組成，分別為孔口表、測斜表、化驗表、巖性表（前3個為必要文件）。此4個表格中的基礎數據主要是依據地質資料而創(chuàng)建，其中孔口表為定位文件、測斜表為軌跡文件、化驗表為基本化驗分析文件、巖性文件為不同區(qū)段的巖性信息。該銀多金屬礦共錄入孔口數據178個、測斜數據1173個、化驗分析數據4721個。

鉆孔數據庫建立后，將勘探工程的空間軌跡關聯起來，形成了鉆孔三維模型（見圖1），實現了圖形與屬性數據的可視化，為下一步礦體三維建模提供準確的數據基礎。

圖1 鉆孔三維模型

3 實體模型的建立

3.1 礦體剖面解譯

將鉆孔數據庫進行三維顯示后，需在每條勘探線上進行礦體輪廓線解譯，即礦體的圈定。本礦礦體圈定原則與地質報告確定的原則保持一致，所采用的工業(yè)指標如下：

邊界品位：Ag≥40×10-6、Pb≥0.30%、Zn≥0.50%、Cu≥0.20%；

工業(yè)品位：Ag≥80×10-6、Pb≥0.70%、Zn≥1.00%、Cu≥0.40%；

最小可采厚度≥1 m；

夾石剔除厚度≥2 m；

礦體有限外推時尖推兩工程間距的1/2作為礦體邊界，無限外推時按相應資源類別基本工程間距的1/2尖推確定礦體邊界，采用米·克噸值（或米·百分值）圈定且工程位于邊部時不再外推。

鑒于I號主礦體分為兩個礦層：中心為大脈狀銀鉛鋅硫化物塊狀富礦層，外圍為浸染狀銀鉛鋅硫化物貧礦層，在礦體解譯時兩礦層分別進行了解譯，避免了后續(xù)資源量估值時貧富礦層混合估值而產生偏差。

解譯后的礦體輪廓線見圖2。

3.2 模型建立

礦體模型是一個三維空間三角網連接而成的實體，礦體模型的建立實現了礦體的三維可視化。在建立礦體模型時，不同的礦體分別建立不同的線框模型，利用各剖面解譯完成的礦體輪廓線進行連接，連接方法選擇等邊三角形法。其中I號主礦體分別建立了I-1號富礦層和I-2號貧礦層兩個線框。在I-1、I-2線框創(chuàng)建過程中，局部地段礦體邊界三角網產生交叉，此種情況在當兩個線框相隔較近或成包含關系時常有發(fā)生，不過在Micromine中，兩個線框產生輕微交叉是允許的，并不影響后續(xù)的品位插值，可在礦塊模型線框賦值時，通過賦值的先后順序予以解決。但當交叉發(fā)生在同一個線框中時是不允許的，必須將線框修正至零交叉、零無效邊后方可進行后續(xù)操作。

圖2 礦體輪廓線

建立完成的礦體三維模型見圖3。

4 資源量估算

利用Micromine軟件進行資源量估算的流程：生成井中坐標→線框賦值→樣長組合（特高品位處理）→建立空塊模型→建立搜索橢球體→品位插值（資源類別劃分）→模型報告。

4.1 生成井中坐標

生成井中坐標的目的是將樣品文件中的品位 數據轉換至估值時所需要的三維數據，具體的實現過程是在樣品文件中添加了東坐標、北坐標、高程三個字段，三個字段的數值是由軟件根據孔口數據、測斜數據以及樣品文件中的從、至字段自動計算生成。

圖3 礦體三維模型

4.2 線框賦值

在估值的過程中，某一個線框的品位估值只能利用本線框內的樣品數據來完成，因此，需要對已經生成完井中坐標的樣品文件進行線框賦值，賦值的過程即對號入座的過程，每一個圈入礦體的樣品數據都會添加一個礦體名稱的屬性。該礦I號主礦體分為富礦層和貧礦層，當線框賦值完畢后，富礦層內的樣品只能對富礦體進行估值，貧礦層內的樣品只能對貧礦體進行估值。

4.3 樣品組合

為了確保用于估值的品位按照相同的樣長進行加權平均，保證估值過程不出現偏差，需要對已知采樣樣品進行樣長組合。樣長組合是把已知采樣樣品的長度全部轉化成同一長度，目的在于避免由于樣品長度大使得其品位也高，在應用樣品進行空間插值的時候，使得該樣品分配的計算權值過大。

利用Micromine統(tǒng)計/直方圖功能分別對I-1號富礦體及其他礦體內樣品的樣長進行統(tǒng)計分析，其直方圖分別見圖4、圖5。

從圖4、圖5可以看出，I-1號富礦體內樣長中值為0.6 m，均值為0.65 m，樣長組合時按照0.5 m進行組合；其他礦體內樣長中值為1.0 m，均值為0.95 m，樣長組合時按照1.0 m進行組合。

4.4 建立空塊模型

為了估算整個礦體的平均品位，需將礦體劃分成若干大小相同的塊段模型。建立空塊模型的目的在于，通過模型中已知的采樣點品位，利用特定的數學方法對其進行空間插值，直到每一個塊段都有品位。

圖4 I-1號礦體樣長統(tǒng)計直方圖

圖5 其他礦體樣長統(tǒng)計直方圖

創(chuàng)建空塊模型時，需要注意以下幾點：

（1）空塊模型的范圍要包括全部的礦體；

（2）空塊模型要選擇合適的尺寸；

（3）根據礦體產狀對空塊模型進行適當的旋轉，以便塊更好地落在礦體內；

（4）使用地表DTM模型、礦體線框模型等已知數據，對塊體模型進行賦值，以區(qū)分塊體是否位于礦體內、是否位于地表下；如有必要，也需要區(qū)分該塊體是氧化礦還是原生礦；

（5）次分塊方法：次分塊需要分的足夠小，以便能更好地模擬邊界（如礦體邊界），但是同時也需要考慮計算機的運算性能。

鑒于本礦床礦體厚度較薄，空塊模型大小確定為5 m×1 m×5 m，次分塊尺寸為2.5 m×0.25 m× 2.5 m，空塊模型經次分后共包含459 865個單元塊。

4.5 建立搜索橢球體

搜索橢球包含與礦塊品位估值相關的樣品品位，不包括無關的（不需要的）品位。根據礦體的產狀設置搜索橢球體的參數，以便更快更準地進行插值。

在建立搜索橢球體時，先將搜索橢球體的顯示中心落在礦體中間，以便能更好地比對礦體形態(tài)；搜索橢球體第一回次的搜索半徑一般為勘探間距的1.2～1.5倍，本礦選擇1.2倍即60 m為初次搜索半徑；扇區(qū)選擇為8個，每個扇面的最多點數為6。

4.6 品位插值

根據創(chuàng)建的空塊模型，設置的搜索橢球體參數，即可對每個空塊進行品位插值。品位插值的方法有多種，主要有距離冪次反比法、克里格法、多重指示克里格等，其中以距離冪次反比法應用最為廣泛。由于各樣品距離待估點的距離不同，其品位對待估點的影響程度也不同。距離待估點越近的樣品，其品位對待估點的影響也就越大。因而在計算中，離待估點近的樣品的權值應比離待估點遠的樣品的權值大。距離冪次反比法就是基于這一思想產生的。

在估值過程中，采用1次估值往往不能將待估塊段全部完成估值，往往需要2次或3次才能完成。在估值完成的塊模型中不應有缺失的品位值，在礦體內部的塊應該檢查是否全部都已插值。

在Micromine軟件中，資源類別的劃分與傳統(tǒng)塊段劃分方法有所不同。在軟件中，當建立搜索橢球體估值參數時，可通過設置搜索半徑和參與估值的工程數對資源類別進行劃分，但這往往與地質報告中傳統(tǒng)地質塊段法統(tǒng)計出來的各類別資源量相差較大，且現階段來講，通過軟件劃分的資源類別還未能普遍運用于報告編制或礦山設計中。因此，為減少不同類別資源量的誤差，有必要圈定不同類別資源量的線框模型，資源類別線框模型可參考地質報告進行圈定。資源類別線框圈定完成后，通過線框賦值即可實現對每個單元塊資源類別賦值。

該銀多金屬礦品位插值完成后的塊體模型見圖6、圖7。從Ag品位分區(qū)間顯示可以看出，I-1號富礦體品位明顯比其他礦體要高。

進行資源類別賦值后的塊體模型見圖8。從圖8可以看出，控制程度較高的地段集中于上部由坑道控制的1200 m、1100 m中段。

圖6 I-1號富礦體塊體模型

圖7 其他礦體塊體模型

圖8 賦值后的塊體模型

4.7 模型報告

依據插值模型即可進行資源量報告。資源量報告時，可設置各類參數分類進行報告，如分礦體、資源類別、礦石類型、不同中段標高、動用/保有等，其報告過程非常便捷。

5 資源量報告及對比

本礦通過建立Micromine三維模型采用距離冪次反比法估算的資源量與地質報告采用傳統(tǒng)的地質塊段法估算的資源量結果進行對比分析，見表1。

由表1可知，核實報告與模型計算的資源量在礦石量、品位、金屬量方面的相對誤差均在5%以內，從不同資源類別的礦石量誤差來看，其差別也是非常小的，由此可見，模型估算的資源量結果是可靠的，可以作為下階段設計的依據。

表1 資源量對比

6 結論

（1）利用Micromine礦業(yè)軟件建立該礦三維地質模型，采用距離冪次反比法進行品位估值，通過與地質報告資源量對比，其誤差在規(guī)范允許范圍內，說明模型估算的資源量是可靠的。

（2）礦床三維模型可全方位立體地展現礦體的三維空間形態(tài)及礦體品位分布特征，實現了礦山的可視化、數字化；通過估值完成的塊體模型，可快速地計算各種約束條件下礦體的資源量；通過礦體三維空間的顯示，可指導下階段設計中采礦系統(tǒng)及開拓工程的布置，可大大提高礦山設計工作效率。