個舊高松礦區三維地質建模及儲量估算

睢 瑜，陳建國，趙江南

（1.江蘇省有色金屬華東地質勘查局，南京 210007；2.中國地質大學（武漢）資源學院，武漢 430074）

0 引言

三維地質建模和可視化概念最早是由加拿大Simon.W.Houlding（1993）提出的。所謂三維地質建模（3D Geosciences Modeling）①是指采用適當的數據結構在計算機中建立能反映地質構造的形態和各要素之間關系以及地質體物理、化學屬性空間分布等地質特征的數學模型。3DGM是由勘探地質學、數學地質、地球物理、礦山測量、礦井地質、GIS、圖形圖像學、科學可視化等學科交叉而形成的一門新型學科。隨著計算機軟硬件技術的飛速發展，三維地質建模技術備受關注，并得到了廣泛的研究和應用［1］。本文以云南個舊高松礦區蘆塘壩10號礦群為主要研究對象（圖1），通過地質體三維數字化建模，提取地質控礦作用定量化指標，建立三維礦床模型，并對已知礦體進行儲量計算，為礦山的進一步開采提供參考依據。

圖1 高松礦區地質圖Fig.1 Geological map of Gaosong mining district

1 研究區地質概況

個舊礦區被近SN向的斷裂分割為西區和東區兩部分。西區目前工作程度尚低，僅發現有?。行偷腻a、鉛礦床。但錫、鉛、鋅、鎢、錫、鋁、銻等多種有色金屬礦點星羅棋布，是尋找大型或超大型錫多金屬礦床的遠景區［2］。東區是超大型錫多金屬礦床分布的主要區域，其地表分布有三疊系厚達千米的碳酸鹽巖建造，深部隱伏有燕山中晚期的花崗巖。NE向的一級褶皺構造五子山復式背斜控制了所有礦床的分布，橫跨五子山復式背斜之上的近EW向斷裂控制了五大礦田的分布，即馬拉格、松樹腳、高松、老廠和卡房礦田。礦田內次級褶皺和花崗巖株控制著礦段或礦床。以錫為主的各類礦床的產出主要受有利地層層位、構造、花崗巖三者有機組合聯合控制。

圖2 基于Micromine的三維地質體建模流程圖Fig.2 Flow sheet for the Micromine software－based 3D modeling of geological body

2 Micromine及三維地質建模流程

Mircomine是澳大利亞Micromine礦業軟件有限公司開發的深受國際好評的大型礦業專用軟件。它主要用于地質勘探、資源評估、儲量計算及露天礦和地下礦礦山設計和開采。Micromine提供了與其他數據庫和相關軟件接口的功能，使該系統的數據可被其他數據庫管理系統和相關軟件查詢和編輯，能夠實現各種工程和礦體的三維立體顯示和成圖，并根據地質統計學的方法和原理提供進行礦體品位和儲量估值的各種方法，同時還可以進行礦山的開采設計以及數字地形模型的建立②。

Micromine軟件提供的多種插值方法進行空待估塊段模型的品位插值（如距離反比加權法、泛/普通克里格法及多重指示克里格法），由此建立的三維品位模型是自動生成礦體儲量報告的基礎工作，相比傳統的手動估算儲量有著巨大的優越性。

基于Micromine的三維地質建模流程見圖2。

3 鉆孔及地形模型的三維可視化

鉆孔三維可視化能夠直觀地顯示鉆孔的空間位置及空間狀態，有利于對鉆孔的三維空間變化進行分析，進而研究樣品的變化規律，為建立礦體三維模型做準備。建立鉆孔地質數據庫后，Micromine軟件提供了強大的諸多數據庫應用功能，除了可以實現基本的數據庫管理更新功能外，還可以生成平面鉆孔軌跡，根據需要提取采樣區間、組合區間，生成井中、區間坐標，鉆孔計算與編錄等。此外還可實現鉆孔的三維可視化，同時根據樣品分析數據范圍劃分樣品區間并賦以不同的顏色以表示不同的品位等級。本次研究共收集研究區鉆孔198個，共包含分析化驗樣品43 975個，并根據高松礦田的地形文件構建DTM模型，建立研究區地形三維可視化圖形（圖3）。

圖3 蘆塘壩礦段DTM加載鉆孔三維可視化圖Fig.3 3D map of DTM loading drill hole in fractures

4 三維地質體模型建立

自然界中地質體（礦體、巖體、地質構造等）的形態復雜多變，用規則的幾何體很難進行描述，因此需要用不規則的幾何模型來近似模擬地質體空間分布形態。目前大多數三維礦業軟件都提供了比較方便、靈活、快速的建立不規則幾何模型的功能，從而能形象直觀地描繪地質體的形態和邊界特性［3］。

4.1 地層模型的建立

高松礦田內分布地層為中三疊統個舊組卡房段和馬拉格段（T2g32，T2g51）碳酸鹽巖，加上隱伏花崗巖層。高松礦區地表未見有花崗巖出露，但深部有隱伏的花崗巖（γ35）分布，根據鉆探工程及勘探線地質剖面圖上可以看出，該區花崗巖體產出標高為800～1 200 m，距地表約1 200～1 400 m，由地層實體模型（圖4）可知，其表面形態為南、東、西三面較高，北部較低，即為NE向的凹槽，該巖體巖性與老卡巖體相似，為中細粒黑云母花崗巖，該花崗巖體是區內成礦母巖［2］。

圖4 研究區地層實體模型Fig.4 Real stratigrphic 3D model of the study area

4.2 斷裂模型的建立

本文研究范圍是大馬蘆礦段北東之蘆塘壩10號礦群，主要分布在麒麟山斷裂、大菁東斷裂周圍。

（1）麒麟山斷裂。產于高松礦田中部，走向長＞6 km，傾斜延伸＞800 m，破碎帶寬6～30 m，傾向NE，傾角70°～83°。該斷裂明顯將NE向斷裂錯斷。在其深部礦化漸強，局部有錫礦體賦存，起導礦及容礦作用。并與蘆塘壩斷裂和層間剝離滑動構造聯合控制著蘆塘壩10號礦群的分布。

（2）大菁東斷裂。分布在礦田中部，走向長2～3 km，破碎帶寬5～10 m，局部達30 m，走向310°～320°，傾向NE，傾角66°～88°。角礫巖帶明顯，同時可見碎斑巖及碎粒巖，地表具弱赤鐵礦化。該斷裂往南東走向漸變為EW向，向下含礦性逐漸加強，1 720中段以下更為明顯，區內最大的脈狀錫礦體即產于其中。

根據勘探線地質剖面圖上的斷裂分布，建立了斷裂三維模型（圖5）。

圖5 蘆塘壩礦段斷裂三維立體模型Fig.5 3D model of fractures in Lutangba ore domain

4.3 礦體模型的建立

研究區內礦體的產出受斷裂構造及地層層間構造控制特征明顯，礦床形成與花崗巖密切相關，一類為花崗巖接觸帶硫化礦，而另一類則屬錫石－硫化物礦床經氧化后而成的層間氧化礦床，按產狀及形態可分為2種類型（圖5）：

（1）陡傾斜脈狀礦體：該類型礦體是成礦熱液沿大菁東斷裂及其派生的斷裂破碎帶充填交代而形成的，礦體產狀較陡，傾角一般60°～85°，走向長數十米至數百米，傾斜延深數十米至600余m，厚0.5～20 m不等；礦體在斷裂帶中呈大小不等的脈狀或透鏡狀產出；礦石以土狀、土塊狀赤褐鐵礦為主，一般含錫鉛品位較高，同時伴生有鋅、銦等有用元素。

（2）緩傾斜似層狀礦體：該類型礦體主要分布在NE向蘆塘壩斷裂旁側的T2g61層位中，少數產于T2g51層位，受主斷裂旁側的次級層間剝離或層間破碎帶控制，礦體呈平行多層狀疊瓦式產出，產狀與地層基本一致，走向近EW向，傾角較緩（1°～30°），呈似層狀、層狀、透鏡狀、長條狀等形態。礦層數較多，一般為3～5層，最多可達9層；礦體規模較大，長數十米至數百米，最長可達千米，寬為100～200 m；礦石以土狀赤褐鐵礦為主，少數膠狀褐鐵礦，局部還夾有未完全氧化的硫化礦殘塊。

從蘆塘壩10號礦群的礦體三維模型（圖6）可以比較直觀地看出，蘆塘壩主要層間氧化礦體分布在標高1 450～1 720 m標高之間，礦體多成層狀分布；由上至下，礦體的范圍呈逐漸減小的趨勢。而在標高1 200 m處，可見夕卡巖硫化礦體，位于隱伏花崗巖與碳酸鹽巖接觸帶成礦，常在花崗巖凹兜部位形成礦體。

圖6 蘆塘壩礦段礦體三維立體模型Fig.6 3D model of ore bodies in Lutangba ore domain

5 儲量估算

5.1 空塊模型的建立

建立蘆塘壩10號礦群三維實體模型后，通過Micromine軟件建立礦體的空待估塊模型。待估塊段模型亦是三維實體，因此就需要以三維物體的定義方式確定待估塊段模型的各種參數［4］。按照1∶2000礦區勘探線，結合礦體的形態、走向、傾向和空間分布特征，建立的三維實體的東坐標、北坐標及高程信息結合蘆塘壩10號礦群空間分布特點與工業生產要求，確定了建模的范圍和基本參數。模型單元塊的行、列、層為6 m×6 m×2 m，次分塊的規格為3 m×3 m×l m。在建立三維網格模型后，可以將找礦數字模型所確定的預測參數作為屬性賦給每一個單元塊。運行后就可得到相應實體模型的空待估塊段模型。由于鉆孔的終孔深度沒有到達賦存于花崗巖實體的硫化礦體，所以本文中僅對10號礦群的其他礦體做儲量計算。創建的蘆塘壩10號礦群的待估空塊模型見圖7。

5.2 儲量計算

本次建模分別以0.1％，0.2％，0.5％，1.0％和10％為邊界品位，分別計算礦體中錫的儲量；對于＜0.1％和＞10％的品位值做為異常值剔除。對蘆塘壩10號礦群進行儲量計算時，礦石的體重取2.355 t/m3。本文研究工作進行儲量計算采用距離反比加權方法、普通克里格法對前面建立的待估塊模型進行品位插值。

圖7 10號礦群待估塊段模型Fig.7 Model of blocks that are not yet evaluarted in No.10 ore body cluster

圖8 錫品位累積頻率圖Fig.8 Accumulative frequency plot of Sn grade

由于研究區收集到的數據未進行過特異值處理，因此在通過Micromine中的過濾功能對組合樣品進行特異值處理，然后再進行組合樣品統計分析，圖8為錫品位數據的累積頻率圖，除了在高值區較為平緩，其余基本呈一直線，反映樣品品位基本服從正態分布。研究中選擇距離反比加權法和普通克里格法對待估塊段模型插值［5］并進行儲量估算［6］，旨在驗證研究區采用品位插值方法計算錫鉛儲量的可行性。

（1）采用距離加權反比法。在數據搜索菜單中的半徑文本框中填入搜索半徑長度，通常該長度為2條勘探線間距的1.25倍，它是一個經驗值。扇區選擇8，扇區可以減少叢聚效應，選擇橢球體，在每一扇面最多點數文本框中輸入數值，該數值決定了每個扇區內參與估算的點的個數，如果每個扇區點數大于該數值，那么只取離估算點最近的幾個點。最少點數（總和）的含義是每個扇區內不少于該數值才進行估算。本文中設搜索半徑分別為40，80，120，其中搜索半徑為120時所有的塊模型全部塊段都被插值完畢。

（2）采用普通克里格法：首先要擬合試驗變異函數，變異函數的作用是用來對區域變量進行結構分析并對樣品品位進行估值，因此必須將得到的試驗變異函數與理論變異函數進行擬合［7］，其中球狀模型為常用的理論變異函數模型，故選擇球狀模型對試驗變異函數進行擬合，得出變異函數的參數值（塊金值、基臺值和變程）。

圖9 錫水平方向的變異函數曲線Fig.9 Variation function curve of Snin horizontal direction

通過2種不同的插值方法計算出來的結果見圖9和圖10。

通過距離乘方反比和普通克里格2種方法的計算，得到蘆塘壩10號礦群錫鉛的儲量，為了比較2種計算結果的差異，分別對2種方法計算所得到的儲量結果進行比較（表1）。從表1可以看出，2種儲量計算方法所得到的累計金屬量的差值比值較小，說明所采用的儲量計算方法可靠，所建立的數學模型比較準確，能滿足儲量計算的要求。

圖10 鉛垂直方向的變異函數曲線Fig.10 Variation function curve of Sn in vertical direction

表1 2種計算方法的結果對比Table 1 Result comparison for the two calculation methods

6 結論

通過對高松礦田蘆塘壩礦段的三維數字化建模，建立了礦體、斷裂、地層的三維地質體模型，將成礦地質信息進行了空間直觀表達，把礦體以數字形式存放的抽象數據用圖形、圖像、動畫等形式、形象表示出來，更直觀地表達礦體分布幾何形態、便于分析礦體分布規律。同時通過不同的插值方法，計算了蘆塘壩10號礦群的錫鉛金屬量，對今后進一步的礦山生產具有一定的指導意義。

注釋：

① 吳立新.3DGM理論與方法.博士生學位課程講義，1999－ 2000版.

② Micromine 11.0培訓教程，2008.

［1］ 肖克炎，趙鵬大.試論大比例尺成礦預測的基本問題及研究途徑［J］.有色金屬礦產與勘查，1994，3（1）：49-56.

［2］ 莊永秋，王任重，楊樹培，等.云南個舊錫銅多金屬礦床［M］.北京：地震出版社，1996.

［3］ 修群業，王軍，高蘭，等.云南金項礦床礦體三維模型的建立及其研究意義［J］.礦床地質，2005，24（5）：501-507.

［4］ 惠勒A J，斯托克斯P C.侯運炳（譯）.塊段模型和線框模型在地下采礦中的應用［J］.Mining Magazine.1988，3：209-213.

［5］ 任利明.用子塊距離加權反比法計算礦石儲量［J］.露天采煤技術，1998（3）：47-48.

［6］ 秦來勇.地質統計學在地質及礦業中的應用與發展［J］.有色冶煉，2002，12（6）：104-106.

［7］ 王仁鐸.地質統計學與礦產資源統計預測［C］.武漢：礦產資源統計預測國際學術討論會（IWSPMR＇90），1990.