山東焦家金成礦帶三維地質建模及成礦預測

，，，，

(1.山東省地質調查院，山東 濟南 250013； 2.中國地質大學(北京)，北京 100083)

0 引 言

21世紀以來，隨著計算機圖形學技術的迅速發展和三維空間數據處理能力的不斷增強，將三維GIS與三維地質建模技術結合起來的三維立體成礦預測技術應運而生。三維地質建模這一概念最早是由加拿大學者Houlding于1994年提出的，即用三維數據模型包括鉆孔數據(Hole data)、圖形數據(Map data)和體元數據(Volume data)以及三維網格數據(3D grid data)為數據源的數據流來展示地質描述。1999年，首屆“國際數字地球”大會提出“數字礦山”的概念。隨后，北歐和美國也建立了“地學信息系統三維建模”的理論。自20世紀80年代至今，我國的礦產資源定量預測評價理論和方法也相繼建立(程裕淇等，1983；翟裕生，1999；趙鵬大，2002；王世稱等，2003；鄧聚龍，2005；王世稱，2010)。

目前，三維地質建模廣泛應用于地下水資源利用與管理、土地資源利用、城市地下空間利用與決策、礦產資源定量預測等領域。在市場需求下,一些礦業軟件，如加拿大的MicroLynx，澳大利亞的Micromine、Surpac Vision和Vulcan，法國的GOCAD以及美國Earth Vision等軟件相繼問世。Surpac作為全球同類軟件中應用最為廣泛的大型三維礦業軟件，具備強大的三維繪圖功能，廣泛應用于資源估算、礦山規劃、生產計劃管理的各個階段，服務于120多個國家的露天和地下開采礦山及勘探項目。

1 礦區地質概況

山東焦家金成礦帶位于著名的膠西北金礦密集區。其中已探明的特大型金礦床3處，分別為新城、焦家、河西金礦床。大、中型金礦床有望兒山、河東、東季、上莊等30處，小型金礦床、金礦點近百處(孔慶友等，2006)。大地構造位置處于華北板塊(Ⅰ)膠北地塊(Ⅱ)膠北隆起及凹陷區(Ⅲ)膠北隆起(Ⅳ)之膠北凸起(Ⅴ)西緣。古老的基底變形變質巖系，多期多成因的巖漿活動和以北東向斷裂為主的構造格架，構成了優越的成礦地質背景。研究區內地層簡單，斷裂構造發育，巖漿巖廣布(圖1)(宋明春等，2012)。

研究區出露的地層有第四系和荊山群祿格莊組。第四系分布于萊州—龍口斷裂以西，荊山群(Pt1J)出露于研究區齊山鎮及招遠城西南，面積僅1.3 km2，呈近東西向殘留體分布于棲霞超單元中，巖性為石榴矽線黑云片巖。

圖1 膠西北地區地質礦產略圖

構造以脆性斷裂構造為特征，主要發育北東、北東東向斷裂。區內的斷裂有龍(口)—萊(州)斷裂、望兒山支斷裂和靈山—北截斷裂。龍(口)—萊(州)斷裂是膠東西北部三大控礦構造之一，在區域上呈“S”型展布，發育有連續穩定的破碎蝕變巖帶，主裂面以厚2～40 cm斷層泥為標志，壓扭性特征明顯。望兒山支斷裂位于焦家主干斷裂下盤，屬分支構造，在朱宋地段與主干斷裂交匯。在研究區的北部，沿玲瓏超單元與郭家嶺超單元接觸帶展布。斷裂走向長10 km，寬80～120 m，總體走向35°，傾向北西，傾角30°～50°。靈山溝—北截斷裂，走向50°左右，傾向南東，傾角65°，發育有破碎蝕變巖帶，斷層泥薄而不連續。3條斷裂均屬壓扭性質，均切割玲瓏超單元，斷裂金礦化比較發育。

區內巖漿巖廣布，以新元古代震旦期玲瓏超單元為主體；新太古代五臺—阜平期馬連莊超單元、棲霞超單元分布于焦家斷裂帶以西；中生代燕山早期郭家嶺超單元侵入玲瓏超單元內，分布于焦家斷裂帶以東。區內派生脈巖較發育，常見閃長玢巖、煌斑巖、花崗斑巖、花崗閃長斑巖等。脈巖走向多為近南北向，少數為北北東向和北東東向，傾向北西或南東，傾角比較陡。

2 地質數據庫的建立

三維資料數據庫是本次研究的數據基礎，涉及大量資料的收集、整理和預處理工作。資料的有效整合是本次研究得以順利進行的第一步和重要前提。收集和應用到的資料包括：勘探線剖面圖(69條)、圖切地質剖面(32條)、可控源音頻大地電磁測量剖面(44條)、井中元素異常剖面圖(11條)、地質勘查規劃平面部署圖、鉆孔編錄數據(680個)、區域地質圖、研究區等高線地形圖和礦田內主要勘查區勘查報告等。這些資料對建立焦家金礦帶三維地質體模型和找礦模型起到了重要的作用(圖2)。

3 三維地質模型的構建

3.1 Surpac軟件建模

Surpac軟件采用實體構模法建立三維地質模型。實體構模法是采用多邊形網格來描述地質和開采過程中形成的形體邊界，并用傳統的塊段模型描述形體內部的品位或質量的分布，實體構模法實質是線框構模法與擴展的塊段構模法的耦合，因此彌補了塊段構模處理邊界的不足。

在系統研究焦家金成礦帶金礦床的基礎上，應用目前主流地質三維建模軟件Surpac，本次研究對研究區巖體、已知礦體、化探、物探異常及鉆孔等進行三維實體建模，從而實現數字模型的建立，并且結合研究區三維單元立方體的劃分，最終建立三維立方體綜合預測模型，指導中比例尺三維成礦預測。研究內容為地質模型分析、三維數字礦山建立、控礦因素分析、地質統計分析及預測結果分析。技術路線如圖3所示。

圖2 焦家金成礦帶實際材料圖

圖3 技術路線圖

圖4 焦家金礦帶巖體實體模型

3.2 三維數字礦山的建立

焦家帶三維數字礦山的建立包括地質體模型及工程模型，其中前者包括巖體、構造、蝕變帶、礦體等，后者包括鉆孔等。

3.2.1 巖體實體模型 區內的巖漿巖以中生代燕山早期玲瓏序列為主體。新太古代五臺—阜平期馬連莊序列、棲霞序列分布于焦家斷裂帶西南側；中生代燕山晚期郭家嶺序列侵入玲瓏序列內，主要分布在焦家斷裂帶以東(圖4)。

3.2.2 構造實體模型 為了清晰準確地反映研究區的構造，依據勘探線剖面圖、中段平面圖和報告中關于斷裂的形態描述，最終得到焦家金礦帶主要斷裂模型(圖5)。研究區內斷裂包括焦家主干斷裂、侯家斷裂、望兒山斷裂、河西斷裂、靈北斷裂、苗家斷裂等。

圖5 焦家金礦帶斷裂實體模型

3.2.3 蝕變帶實體模型 該區金礦床均賦存于斷裂帶所控制的蝕變巖帶內。蝕變巖帶的發育空間、形狀產狀與對應的斷裂帶基本一致。主要礦體多賦存于主裂面下盤的絹英巖化花崗質碎裂巖帶和黃鐵絹英巖化碎裂巖帶內。圖6為焦家金礦帶蝕變帶模型。

圖6 焦家金礦帶蝕變帶模型

3.2.4 已知礦體實體模型 預測研究區內已發現和評價的金礦床(點)基本沿北東向斷裂構造分布。規模較大的焦家主干斷裂帶及其下盤的分支斷裂基本控制了區內主要特大、大型金礦床的分布，次級斷裂控制該區大部分的中、小型金礦床。圖7為焦家金礦帶已知礦體相對位置圖。

3.2.5 地球物理模型 根據收集到的可控源大地音頻電磁測量剖面上的視點祖率等值線，將其在Surpac軟件中連接生成曲面，建立研究區南部朱橋地區的視電阻率等值線曲面模型(圖8)。

圖8 焦家金礦帶視電阻率異常三維模型

3.2.6 鉆孔模型 將孔口坐標表、測斜數據表、樣品數據表按規范的格式要求進行整理后導入Surpac軟件中形成鉆孔數據庫，可以應用Surpac鉆孔三維顯示功能進行瀏覽。圖9顯示，鉆孔集中在礦帶中部地區。通過沿某一剖面(線)兩側一定距離截取切面的方式顯示鉆孔，可以仔細觀察到這一范圍內鉆孔的軌跡和樣品屬性等，也可以通過這一方式進行三維環境下的地質剖面勾繪和礦體的圈定工作。此外，Au元素三維異常分布情況也依鉆孔模型插值分析得到。

圖9 焦家金礦帶鉆孔三維模型

3.2.7 地球化學模型 根據收集的研究區中部地區11條井中元素異常剖面信息，恢復其真三維空間位置，依據剖面圖上圈定的Hg、Ag、Sb、Pb、As、Au、Zn、W、Sn、Mo、Bi、Cu 12種元素的異常界線，將相鄰剖面的異常界線相連接，建立元素異常組合的三維實體模型(圖10)。

圖10 焦家金礦帶元素組合異常三維模型

4 成礦預測及資源量估算

4.1 預測模型建立及預測變量統計

根據對研究區找礦模型的分析及有利成礦信息的提取，結合實際情況，建立了焦家金礦帶預測模型(表1)。

分析方法詳見Chen 等(2007)。

表1 焦家金礦帶定量預測模型

4.2 三維證據權計算

證據權重法是加拿大數學地質學家Agterberg(1990)提出的一種地學統計方法。它采用一種統計分析模式，通過對一些與礦產形成相關的地學信息的疊加復合分析來進行礦產遠景區的預測。其中的每一種地學信息都被視為成礦遠景區預測的1個證據因子，而每一個證據因子對成礦預測的貢獻是由這個因子的權重值來確定的。證據權模型既考慮了地質存在的找礦權重，又考慮了地質因素缺失的找礦權重，實際上，后驗概率就是在先驗概率的基礎上對證據權的正負疊加。

證據權法在本研究中的實現過程如下。

首先，計算先驗概率，計算公式如下：

(1)

(2)

(3)

其次，計算權重值：

(4)

(5)

(6)

(7)

(8)

(9)

最后，后驗概率的計算如下：

(10)

(11)

(12)

以上公式中，Nj為研究區內具有標志A的含礦單元數；N為研究區內的含礦單元數；Sj為研究區內具有標志A的單元數；S為研究區的單元總數。由以上公式計算出來的焦家金礦帶金礦找礦標志權值見表2。

表2 焦家金礦帶金礦找礦標志權值

4.3 三維信息量計算

信息量法應用于區域礦產預測，是由維索科奧斯特羅夫斯卡婭(1968)及恰金(1969)先后提出的(趙鵬大等，1983)。

首先，計算各地質因素、找礦標志所提供的找礦信息量，定量評價各地質因素和標志對指導找礦的作用；其次，計算每個單元中各標志信息量的總和，其大小反映了該單元相對的找礦意義，用以評價找礦遠景區進行預測。找礦標志的找礦信息量用條件概率計算，即:

(13)

式(13)中,IA(B)為A標志有B礦的信息量；P(A/B)為已知有B礦存在的條件下出現A的概率；P(A)為在研究區內出現標志A的概率。

由于概率估計上的困難，以頻率值來估計概率值。此時：

(14)

式(14)中,Nj為研究區內具有標志A的含礦單元數；N為研究區內的含礦單元數；Sj為研究區內具有標志A的單元數；S為研究區的單元總數。

一般二維找礦信息量法用式(15)確定找礦有利標志。

(15)

為了能夠對三維證據權的預測結果進行系統的評價，又利用三維信息量法對數據進行了相關的計算。表3為找礦信息量計算結果。

表3 找礦信息量計算結果

4.4 預測結果

由表2和表3的統計信息量和后驗概率所包含的已知礦塊比例可以直觀地看出本次預測的統計取值。經統計,信息量值分3個級別，分別是A級2.75～3.15、B級3.15～3.95、C級≥3.95；后驗概率值取0.6作為本次預測的最低限制條件。經統計，符合后驗概率和信息量范圍的有利成礦數有106 803個，占研究區總立方塊數(1 193 148)的8.95%，具有很好的預測前景。

將三維證據權法的預測結果與三維信息量法的預測結果進行相交分析(陳建平等，2009)，統計成礦有利地塊的含有信息量高值(包括A、B、C)的立方塊數(圖11)，然后根據找礦信息量的大小及空間位置，劃分出6個預測區，按級別依次為A1-1、B3-1、B3-2、B3-3、C2-1、C2-2(圖12)。

圖11 焦家金礦帶靶區各信息量區間立方塊數統計圖

圖12 焦家金礦帶遠景預測區立體圖

4.5 資源量估算

本次資源量估算采用了體積估計法和豐度估計法2種方法。

體積估計法就是把已知地區有代表性的單位體積礦產平均含量估計值外推到研究地區體積內的資源， 這種估計方法稱體積估計法。具體做法是，用控制區地殼單位體積內的礦產平均含量乘以所研究地區的總體積，得到礦產資源估計值。公式如下：

(16)

式(16)中，T1為已知地區礦床儲量，t；T2為預測區礦床儲量，t；V1為已知地區礦床體積，m3；V2為預測區礦床體積，m3。

豐度估計法就是通過求出已知地區成礦元素的富集系數公式(17)，并外推到預測區去的辦法來求預測資源，其計算公式見式(18)。公式中，A值取自研究區內玲瓏序列Au的豐度值0.001 79 g/t (李士先等，2006)。

(17)

(18)

式(18)中，S為已知地區面積，km2；h為已知地區深度，km；G為含該元素的巖石密度，t/m3；A為已知地區該元素的豐度，g/t；T為已知地區礦床儲量，t。

運用體積估計法計算得出的找礦有利區總量706.107 2 t，6個靶區資源總量為382.482 8 t。運用豐度值估計法計算得出各成礦有利區資源總量為699.064 5 t，6個靶區資源總量為382.466 3 t。2種方法計算的結果很接近，具有一致性。

5 結 論

(1) 本次研究基于Surpac軟件成功地建立了山東焦家金礦帶的三維數字礦山模型，并分析了該成礦帶的成礦要素特征值，建立了預測模型。使用地質學統計方法計算得出成礦有利區域的信息值，按照信息值的等級依次圈定了6個成礦靶區。本次研究首次把焦家金礦帶的二維地質資料有效整合，建立了山東焦家金成礦帶的三維數字礦山模型，填補了這一重要金成礦帶三維數字化的空白，對該地區隱伏礦體的找礦突破具有指導意義。

(2) 通過體積估計法和豐度值估計法2種完全不同的資源量計算方法，并且通過查閱收集完全不同的數據來源計算出來的找礦有利區資源量都約為700 t，找礦靶區資源量約為380 t，證明本次資源量的估算結果較為準確。

程裕淇，陳毓川，趙一鳴.1983.再論礦床的成礦系列問題：兼論中生代某些礦床的成礦系列[J].地質評論，29(2)：127-139.

陳建平，尚北川，呂鵬，等.2009.云南個舊礦區某隱伏礦床大比例尺三維預測[J].地質科學，44(1):324-337.

鄧聚龍.2005.灰色系統基本方法[M].武漢：華中科技大學出版社.

孔慶友，張天禎，于雪峰,等.2006.山東礦床[M].濟南：山東科學技術出版社.

李士先，王建收.2006.膠東“玲瓏-焦家式”金礦資源潛力與找礦[J].山東國土資源,22(12):36-41.

宋明春,劉明渭,張淑芳,等.2012.山東省地質礦產圖集[M].濟南:山東省地圖出版社.

王世稱，劉玉強，伊丕厚.2003.山東省金礦床及金礦床密集區綜合信息成礦預測[M].北京：地質出版社.

王世稱.2010.綜合信息礦產預測理論與方法體系新進展[J].地質通報，29(10)：1399-1403.

趙鵬大，胡旺亮，李紫金.1983.礦床統計預測[M].北京：地質出版社.

翟裕生.1999.論成礦系統[J].地學前緣，6(1)：13-27.

趙鵬大.2002.三聯式資源定量預測與評價：數字找礦理論與實踐探討[J].地球科學：中國地質大學學報，27(5)：482-489.

AGTERBERG F P.1990. Combining indicator patterns for mineral Resources evaluation[C]// Proceedings of International Workshop on Statistical Prediction of Mineral resources. Wuhan, China:China University of Geoscienices(Wuhan), 1:1-15.

CHEN JIANPING, LU PENG, WU WEN,et al. 2007.A 3D prediction method for blind orebody based on 3D visualization model and its application[J]. Earth Science Frontiers, 14(5): 54-61.

HOULDING S W.1994.3D Geosciences Modeling-computer Techniques for Geological Characterization[M].Berlin:Springer-Verlag.