王江霞 ,陳建平,張 瑩,鄭永琴,陳東越

(1. 中國地質(zhì)大學(北京)地球科學與資源學院，北京 100083；2. 中國地質(zhì)大學國土資源與高新技術(shù)研究中心，北京 100083；3. 中國地質(zhì)大學北京市國土資源信息開發(fā)研究重點實驗室，北京 100083；)

基于GIS的證據(jù)權(quán)重法在冀東地區(qū)多元信息成礦預測中的應用

王江霞1,2,3,陳建平1,2,3,張 瑩1,2,3,鄭永琴1,2,3,陳東越1,2,3

(1. 中國地質(zhì)大學(北京)地球科學與資源學院，北京 100083；2. 中國地質(zhì)大學國土資源與高新技術(shù)研究中心，北京 100083；3. 中國地質(zhì)大學北京市國土資源信息開發(fā)研究重點實驗室，北京 100083；)

本文通過詳細分析和系統(tǒng)總結(jié)前人的研究成果，充分利用目前新的成礦預測理論和礦產(chǎn)資源勘查與評價理論和方法技術(shù)，結(jié)合GIS技術(shù)分別對冀東地區(qū)主要的沉積變質(zhì)型鐵礦床成礦規(guī)律進行分析總結(jié)，并對與該區(qū)成礦相關(guān)的地質(zhì)信息、物探信息、化探信息等多元化信息歸納，總結(jié)出該區(qū)找礦概念模型。繼而在證據(jù)權(quán)重法的基礎上采用網(wǎng)格單元法，圈定了找礦遠景區(qū)，實現(xiàn)了對冀東地區(qū)沉積變質(zhì)型鐵礦資源的綜合預測研究。

冀東 沉積變質(zhì)型鐵礦 證據(jù)權(quán)重 多元信息

Wang Jiang-xia, Chen Jian-ping, Zhang Ying, Zheng Yong-qin, Chen Dong-yue. Application of evidence weight method basing on GIS to metallogenic prediction in East Hebei Province[J]. Geology and Exploration, 2014, 50(3):0464-0474.

0 引言

冀東地區(qū)是我國僅次于鞍(鞍山)本(本溪)的第二大鐵礦資源基地，河北省沉積變質(zhì)型鐵礦的96%資源量在冀東變質(zhì)巖區(qū)，而80%的資源量又集中遷安水廠-灤縣司家營一帶(河北省地質(zhì)調(diào)查院，2010)。迄今為止,前人主要在冀東地區(qū)金礦成礦有深入的研究并取得了重要進展，王自力(2008)、孔德鑫(2013)等對峪耳崖金礦成礦物質(zhì)進行了研究；宋揚(2011)、牛樹銀(2012)等對金廠峪金礦成礦模式和預測做了詳細研究；楊彥辰(1995)、李俊健(2002)等在冀東金礦的控礦條件和找礦方向做了深入研究。而對冀東鐵礦的研究相對較少，本文綜合地質(zhì)信息、地球物理信息以及地球化學信息的多元化信息，分析總結(jié)冀東沉積變質(zhì)型鐵礦床的成因及規(guī)律，建立找礦模型，并在證據(jù)權(quán)方法和網(wǎng)格單元法的基礎上，完整而全面地進行礦產(chǎn)資源預測。

1 區(qū)域成礦要素研究

冀東地區(qū)位于北京以東，青龍-山海關(guān)以西；北界在承德至平泉一帶，南界在寶坻至灤南一線，面積約為40000km2。1∶25萬圖幅有承德市幅、青龍縣幅、秦皇島幅、天津市幅；地理坐標E117°～119°30′；N39°20′～41°20′。河北省變質(zhì)鐵礦礦產(chǎn)地117處，累計探明資源量69.23億噸。冀東沉積變質(zhì)礦產(chǎn)地103處，累計探明資源量66.94噸，占全省沉積變質(zhì)型鐵礦資源總量的96.7%。因此選取了該區(qū)最主要的沉積變質(zhì)型鐵礦床為研究對象。

冀東地區(qū)比較典型的沉積變質(zhì)型鐵礦有：遵化石人溝、遷安水廠、司家營、青龍柞欄杖子等，其成礦過程是：沉積含鐵物質(zhì)，盆地下沉壓實成巖成礦，經(jīng)麻粒巖、角閃巖相及綠片巖相多次疊加變質(zhì)，多期褶皺、斷裂變形；區(qū)域性隆起與凹陷，風化剝蝕?？偨Y(jié)分析這些礦床的成礦要素研究得出區(qū)域成礦要素：沉積作用及大地構(gòu)造環(huán)境、構(gòu)造作用、變質(zhì)作用。

沉積作用及大地構(gòu)造環(huán)境：接受含鐵物質(zhì)沉積盆地的大地構(gòu)造位置，是主要成礦要素。沉積盆地周圍一定距離內(nèi)有豐富的含鐵物質(zhì)來源，盆地有利于接收這些沉積物。

水廠式鐵礦分布于遷安水廠一帶出露的地層為遷西巖群，是冀東地區(qū)主要含鐵巖系。主要巖性為一套閃巖相-麻粒巖相變質(zhì)表殼巖系，由二長斜長麻粒巖、紫云蘇麻粒巖夾石榴淺粒巖、磁鐵石英巖、二輝斜長角閃巖和磁鐵石英巖等組成，其原巖為一套火山巖-沉積巖系，巖石類型為拉斑玄武巖、安英巖及雜砂巖、高鋁質(zhì)巖。大地構(gòu)造環(huán)境為初始克拉通火山盆地。

石人溝式鐵礦出露的遵化巖群，由角閃巖相變質(zhì)相系的斜長角閃巖、斜長透輝閃巖、黑云斜長變粒巖、條帶狀磁鐵石英巖等組成，其原巖為拉斑-鈣堿性玄武巖加雜砂巖。大地構(gòu)造環(huán)境為陸塊邊緣及弧后盆地。

司家營鐵礦含礦地層為晚太古代灤縣群,含礦系為黑云變粒巖，局部有少量角閃巖、黑云長石石英巖和大理巖，低角閃巖相變質(zhì)巖系。原巖為拉斑玄武巖夾少量酸性火山巖、泥質(zhì)沉積巖。大地構(gòu)造環(huán)境為陸塊邊緣島弧及弧后盆地。

柞欄杖子式鐵礦的含礦地層為晚太古代朱杖子群變質(zhì)巖，為一套綠片巖相變質(zhì)巖系。主要巖石為復成分變質(zhì)礫巖、黑云變粒巖、斜長角閃巖、磁鐵石英巖。原巖為泥質(zhì)、粉砂質(zhì)沉積巖及拉斑玄武巖，為新太古代末期膠遼、晉冀和豫皖古陸塊拼接帶。

構(gòu)造作用:長期多次的褶皺、斷裂、隆起、凹陷等構(gòu)造作用，使得原始簡單的層狀巖、礦層，形成復雜多樣的礦帶格局和礦體形態(tài)，同時使得現(xiàn)今地表能看到從古代到新太古代長達10億年的各個沉積建造形成的礦床。冀東地區(qū)變質(zhì)基底最主要的構(gòu)造作用就是褶皺作用，礦體在平面上和剖面上都呈現(xiàn)著不同尺度的褶皺，是礦帶分布和礦體形態(tài)的控制因素，全區(qū)普遍是向斜部分保存礦體，背斜部分礦體剝蝕，褶皺樞紐總體向南傾伏，向北抬起，褶皺的轉(zhuǎn)折端和同斜緊密褶皺使礦體加厚。

變質(zhì)作用：變質(zhì)作用可使含鐵物質(zhì)變質(zhì)為磁鐵礦，增加礦石可選性，提高礦床的經(jīng)濟價值。冀東地區(qū)橫跨兩個一級變質(zhì)地質(zhì)單元。東部為膠遼變質(zhì)域的秦皇島新太古代巖漿弧，變質(zhì)地層沿北東向青龍河斷裂帶局限分布，為低角閃巖相-綠片巖相變質(zhì)巖石，西部主體是華北變質(zhì)域的晉冀古陸塊的遷西陸核，是我國最古老的變質(zhì)結(jié)晶基底，由于中新太古代的多期次中酸性巖漿侵入、麻粒巖相到高角閃巖相變質(zhì)和多期次變形作用，使得變質(zhì)地層呈殘留體或包體形式“漂浮”在變質(zhì)深成雜巖中。

2 GIS多元成礦信息分析

GIS(地理信息系統(tǒng))，具有對地球空間數(shù)據(jù)和屬性數(shù)據(jù)進行采集、存儲、檢索、建模、分析、輸出的功能，利GIS技術(shù)可以很好地將各種地質(zhì)、礦產(chǎn)、物探、化探、遙感等數(shù)據(jù)實現(xiàn)分層存儲和管理， 進行各種必要的空間分析。因此，此項研究是在建立了較為齊全的空間數(shù)據(jù)庫的基礎上進行的，其中包括礦點、地質(zhì)圖、物探圖件(重力、航磁)、化探圖件等。要進行這些分析工作，首先要有充分的數(shù)據(jù)基礎和較高質(zhì)量的各種圖件，并能從中推斷、進而挖掘隱含的成礦信息， 這也是應用證據(jù)權(quán)重法的一個基本條件。

2.1 成礦地質(zhì)信息分析

2.1.1 賦礦地層

沉積變質(zhì)型鐵礦賦存于太古界變質(zhì)巖系中，嚴格受地層層位控制。主要含礦層位為遷西群三屯營組、馬蘭峪組、單塔子群白廟子組和朱杖子群上白城子組、悖羅臺組(任樹祥等，2010)。其中遷西群為最古老的結(jié)晶基底。經(jīng)礦點與地層疊加分析發(fā)現(xiàn)，礦點分布在跑馬廠組(Ar3p)中占17.65 %，拉馬溝組(Ar3l)中占14.07 %，遵化雜巖(Ar3zh)中有13.56 %，單塔子雜巖(Ar3d)中有12.53 %，遷西雜巖(Ar2qx)有4.9 %。地層組合熵反映了地層構(gòu)造的復雜度，也反映事物發(fā)生的不確定度。一般說來，地質(zhì)構(gòu)造特征越復雜，不確定程度越高，熵值越高(孫巖，2010)。將得到的熵值與礦點進行疊加統(tǒng)計可以看出，74.17%礦點落在熵值的次高區(qū)(59,80]。

2.1.2 控礦構(gòu)造

成礦前斷裂一般具有控巖、控礦的作用，經(jīng)?？刂频V體的總體展布格局及成礦期斷裂活動的范圍及特性， 成礦期斷裂是控制礦化富集的主導因素， 成礦后斷裂活動常常對礦體起破壞作用(趙鵬大，2006)。沉積變質(zhì)型鐵礦床的成礦時代主要是前寒武紀，所以挑選出限制在太古代和元古代地層內(nèi)的構(gòu)造，作為與成礦有關(guān)的構(gòu)造。

區(qū)域內(nèi)主干斷裂就是斷裂等密度與構(gòu)造頻數(shù)的比值大的部位，斷裂等密度是單位面積中斷裂長度加和，而斷裂頻數(shù)單位面積中條數(shù)的加和，其比值大的部分就是密度大而頻數(shù)小的部位，即單位面積內(nèi)斷裂長而條數(shù)少，表現(xiàn)了區(qū)域主干斷裂的特征。經(jīng)礦點與主干斷裂疊加統(tǒng)計選取[1.25，5.625]為區(qū)域主干斷裂的有利區(qū)間。

圖1 有利地層疊加礦點Fig.1 Beneficial stratigraphic and ore occurrences 1-Ar3p;2-Ar3l;3-Ar3zh;4-Ar3d;5-Ar2 qx

圖2 地層組合熵有利區(qū)間疊加礦點Fig.2 Superimposed map of strata combination entropy and ore occurrences 1-地名; 2-礦點; 3-斷裂; 4-地層組合熵有利區(qū)間1-geographic name; 2-mineral occurrence; 3-fracture; 4-strata combination entropy

圖3 主干斷裂有利區(qū)間分布圖Fig.3 The beneficial range diagram of the main fracture 1-地名; 2-礦點; 3-斷裂; 4-主干斷裂有利區(qū)間1-geographic name; 2-mineral occurrence;3-fracture;4-beneficial range diagram of the main fracture

圖4 斷裂優(yōu)益度有利區(qū)間分布圖Fig.4 The beneficial range diagram of the priority degree for the fracture 1-地名; 2-礦點; 3-斷裂; 4-斷裂優(yōu)益度1-geographic name; 2-mineral occurrence; 3-fracture; 4-priority degree for the fracture

圖5 構(gòu)造2000 m緩沖區(qū)與礦點疊加圖Fig.5 Structure 2km buffer overlying mines 1-地名;2-礦點;3-斷裂;4-斷裂緩沖區(qū)1-geographic name; 2-mineral occurrence; 3-fracture; 4-structure buffer

圖6 研究區(qū)變質(zhì)巖與礦點疊加圖Fig.6 Metamorphic rock over erlying mines 1-地名;2-礦點;3-變質(zhì)巖1-geographic name;2-mineral occurrence;3-metamorphic rock

斷裂優(yōu)益度是以線性構(gòu)造(斷裂)兩兩之間的夾角與線性構(gòu)造方位的控礦度加權(quán)的構(gòu)造密度的量度，其代表研究區(qū)主干構(gòu)造方向成礦的優(yōu)越性，高值區(qū)多為成礦有利地段(董慶吉等，2010)。經(jīng)礦點與優(yōu)益度疊加統(tǒng)計，發(fā)現(xiàn)85.77%的礦點落在 [1.75，10.5]區(qū)間，并與主干斷裂有利區(qū)間圖進行對比發(fā)現(xiàn)，得出的結(jié)果基本相同，都沿著主干斷裂走向，且都分布在斷裂旁側(cè)，這充分說明了該方法的可行性。

區(qū)域性的斷裂構(gòu)造對成礦起著至關(guān)重更要的作用，探討區(qū)域性斷裂構(gòu)造展布特征能夠更好地指明區(qū)域找礦方向。通過該類型礦床成礦規(guī)律發(fā)現(xiàn)，礦體均分布在斷裂的旁側(cè)。根據(jù)實際情況對斷裂作一定范圍內(nèi)的緩沖區(qū)處理，從礦點與斷裂距離統(tǒng)計發(fā)現(xiàn)，在距離2000m內(nèi)含有174個礦點，占總礦點的54.4%，而在距離2500m中僅有180個，可見在緩沖區(qū)2000～2500m之間礦點分布較少，因此，認為本區(qū)構(gòu)造最佳影響區(qū)域為斷裂緩沖2000m范圍。

2.1.3 變質(zhì)作用

沉積變質(zhì)型鐵礦床主要產(chǎn)出在前寒武紀變質(zhì)巖系中(代堰锫等，2012)，對該研究區(qū)的地質(zhì)圖進行巖性分類為沉積巖、侵入巖和變質(zhì)巖，并分別依次和礦點疊加分析發(fā)現(xiàn)，變質(zhì)巖中沉積變質(zhì)型鐵礦點有275個，占總礦點的71.6%，而在沉積巖中的僅有52個，因此，選擇變質(zhì)巖為成礦的有利巖性。

2.2 地球物理成礦信息分析

冀東重力場由南部平原區(qū)向北西山區(qū)呈級次降低趨勢，其幅值為32～-80×10-5m/s2,重力場總體呈東西向。有兩條東西向的梯級帶，一是遵化-山海關(guān)，二是承德大廟-娘娘廟。重力勘探雖然沒有直接用于找礦工作，但是所獲取的重力場信息，通過結(jié)合已有的地質(zhì)資料、鉆孔資料、礦產(chǎn)資料以及航磁、遙感、化探資料進行綜合分析、深入研究，可以發(fā)現(xiàn)和找到與鐵礦有關(guān)的侵入巖和具有一定控礦作用的斷裂構(gòu)造和盆地構(gòu)造等，進而為礦產(chǎn)預測提供深部地質(zhì)構(gòu)造等有價值的信息。

布格重力異常資料是重力勘探的基礎資料，可以通過區(qū)域布格重力異常推斷巖體構(gòu)造情況。將重力布格異常圖與已知礦化點疊加分析，發(fā)現(xiàn)礦點落在布格異常值為(-280,140]區(qū)間內(nèi)的占 62.97 %，該異常區(qū)東北部斷裂構(gòu)造發(fā)育，伴隨著基底褶皺形成的緯向主干斷裂及其派生的北西、北東向次級斷裂，控制了本區(qū)基性-超基性雜巖體及有關(guān)的釩鈦鐵磷礦等礦產(chǎn)的形成。該異常區(qū)主要是太古代變質(zhì)巖系地層，斷裂構(gòu)造發(fā)育，巖漿活動強烈，主要是中酸性巖體，新生界沉積盆地廣泛。重力場宏觀地反映了這些地質(zhì)構(gòu)造現(xiàn)象。

航磁異?？傮w是一個正、負異常相間排列的高值異常區(qū)。從磁場的總體特征看,正的背景場反映了基底的隆起，因為該區(qū)基底巖層大部分為中基性變質(zhì)建造具有較高的磁性，故產(chǎn)生高背景磁場。而磁場的不穩(wěn)定則反映了結(jié)晶基底巖石磁性的不均勻性。將航磁化極異常、化極垂向一階導與礦點疊加分析，航磁化極的含礦區(qū)間為(-280,100]，一階導數(shù)區(qū)間為(-260,280]，將航磁化極異常、化極垂向一階導兩個含礦區(qū)間進行疊加，得到航磁異常成礦有利區(qū)間，其中礦點有176個,占總礦點數(shù)的55.6%。

2.3 地球化學成礦信息分析

化探異常反映地表或近地表區(qū)域內(nèi)成礦元素的富集程度(董英君等，2005)，地球化學異常是一種比較直觀的找礦標志，異常本身的特征在很大程度上反映了礦體的特征。鐵族元素的高低背景分布對應于地質(zhì)單元的展布形態(tài)，太古宇、古元古界及中基性巖漿巖分布區(qū)為高背景，長城系下部、侏羅系、白堊系酸性巖漿巖分布為低背景區(qū)，薊縣系、青白口系、寒武-奧陶系分布區(qū)為過渡區(qū)，地質(zhì)化學特征是地質(zhì)建造類型的良好反映，選取代表鐵元素異常的氧化鐵分布等值線為沉積變質(zhì)型鐵礦的控礦因子，將氧化鐵元素異常圖與礦點疊加統(tǒng)計發(fā)現(xiàn)在(4.56,8.76]礦點有66.53 %。

3 基于模型的證據(jù)權(quán)重方法預測

3.1 找礦模型的建立

礦床的找礦模型是指特定類型中某一典型礦床或同一類型礦床的地質(zhì)-地球物理-地球化學特征、找礦標志與找礦方法組合的基本概述與表述(戎景會等，2012)。找礦模型可突出重要的控礦因素，抓住找礦的關(guān)鍵信息，提出獲得關(guān)鍵信息的有效方法組合，總結(jié)主要找礦標志組合，因而簡化了找礦的實際過程，是提高預測可信程度的主要依據(jù)(曹鐘清，2004；劉洪滔等，2008)。礦產(chǎn)勘查過程中，需要將地、物、化等信息進行有機結(jié)合，并總結(jié)找礦標志，建立找礦預測模型，優(yōu)選有效技術(shù)手段，以此對礦產(chǎn)資源進行有效勘查(葉天竺等，2007)。同樣，預測模型也需要不斷地驗證，并根據(jù)其實際運用效果來及時修正(王登紅等，2010)，通過對冀東地區(qū)沉積變質(zhì)型鐵礦特征的分析研究，同時對比典型礦床成礦地質(zhì)背景和成礦要素，綜合分析歸納并建立了該區(qū)沉積變質(zhì)型鐵礦的找礦預測模型(表1)。

表1 研究區(qū)找礦預測模型Table 1 Prospecting model of the study area

3.2 證據(jù)權(quán)重法成礦預測

本次資源預測統(tǒng)計方法采用證據(jù)權(quán)重法，它是由加拿大數(shù)學地質(zhì)學家Agterberg提出的一種基于二值圖像的地學統(tǒng)計方法(Agterberg，1989)。它采用一種統(tǒng)計分析模式，通過對一些與礦產(chǎn)形成相關(guān)的地學信息的疊加分析來進行礦產(chǎn)遠景區(qū)的預測，其中每一種地學信息都被視為成礦遠景區(qū)預測的一個證據(jù)因子，而每一個證據(jù)因子對成礦預測的貢獻是由這個因子的權(quán)重值來確定的。對于證據(jù)權(quán)，為了便于解釋預測(證據(jù))圖通常采用二態(tài)賦值形式。應用地質(zhì)判斷或統(tǒng)計方法能夠?qū)⑦@種形式主觀地轉(zhuǎn)換成其它形式以確定臨界值，其臨界值能夠最大限度的揭示二態(tài)賦制圖成果模式與數(shù)據(jù)模型的空間組合關(guān)系。證據(jù)權(quán)權(quán)值計算后各證據(jù)因子之間相對于礦點分布滿足條件獨立(徐善法等，2006)，最終結(jié)果是以權(quán)的形式或以后驗概率圖的形式表達的組合圖；證據(jù)權(quán)法的優(yōu)點在于權(quán)的解釋是相對直觀的，并能夠獨立地確定，易于產(chǎn)生重現(xiàn)性。

圖7 布格重力異常成礦有利區(qū)間圖Fig.7 The beneficial range diagram of bouguer gravity anomaly 1-地名;2-礦點;3-布格重力異常1-Geographic name;2-Mineral occurrence;3-Bouguer gravity anomaly

圖8 航磁化極異常成礦有利區(qū)間分布圖Fig.8 The beneficial range diagram of the pole of aeromagnetic prospecting 1-地名;2-礦點;3-航磁化極異常1-geographic name;2-mineral occurrence;3-pole of aeromagnetic prospecting

圖9 氧化鐵元素異常有利區(qū)間圖Fig.9 The beneficial range diagram of iron oxide element abnormal 1-地名;2-礦點;3-氧化鐵異常1-geographic name;2-mineral occurrence;3-iron oxide element abnormal

在運用證據(jù)權(quán)方法進行權(quán)值計算之前，要先對研究區(qū)進行統(tǒng)計單元劃分，單元劃分對最終的預測結(jié)果擁有重大的影響。網(wǎng)格單元法是把研究區(qū)按著一定的間隔劃分成面積相等、形狀相同的若干個單元。此法是由阿萊斯(1957)首先提出的，他用網(wǎng)格法建立了撒哈拉沙漠地區(qū)礦床的泊松分布模型,后又用網(wǎng)格單元建立了判別分析模型，成功地進行了礦產(chǎn)預測工作。網(wǎng)格法劃分出的單元是兩兩相鄰的，構(gòu)成了單元網(wǎng)格，可以覆蓋整個研究區(qū)。本次預測按2000m×2000m對整個冀東地區(qū)進行標準網(wǎng)格單元劃分，然后在MRAS軟件中分別計算各證據(jù)層的證據(jù)權(quán)值(表2)并以此來對研究區(qū)內(nèi)的各個網(wǎng)格單元進行成礦后驗概率計算(劉曉玲等，2010)。運用證據(jù)權(quán)重法預測模型計算得到的后驗概率等值線圖(圖10)。根據(jù)后驗概率的大小與分布，依據(jù)找礦靶區(qū)分類原則，以圖為基礎，劃分出A類找礦遠景區(qū)區(qū)4處；B類靶區(qū)4處；并確定出各靶區(qū)的地理位置。

4 結(jié)論

(1) 在全面收集整理冀東地區(qū)以往地質(zhì)礦產(chǎn)資料的基礎上，系統(tǒng)分析了研究區(qū)地質(zhì)構(gòu)造演化歷史與成礦關(guān)系，對該區(qū)典型礦床的成因類型進行了深入的剖析，最后選取沉積變質(zhì)型鐵礦床作為該區(qū)典型礦床。

(2) 運用GIS技術(shù)，對收集到的物探、化探信息進行綜合分析，分別推斷出構(gòu)造。將地質(zhì)信息、化探信息、物探信息等進行成礦有利分析。依據(jù)典型礦床的成礦要素，歸納總結(jié)出冀東沉積變質(zhì)型礦床的成礦規(guī)律和找礦概念模型。

(3) 運用證據(jù)權(quán)重法結(jié)合網(wǎng)格單元法對該區(qū)沉積變質(zhì)型鐵礦床進行預測；最終根據(jù)勘查學中分別對沉積變質(zhì)型鐵床圈定A類遠景區(qū)4處；B類靶區(qū)4處；得到的預測結(jié)果與已有礦點有很好的疊合度，證明該研究方法的可行性和準確性。預測的結(jié)果表明該區(qū)仍有一定的找礦空間。

表2 研究區(qū)沉積變質(zhì)型鐵礦主要證據(jù)層權(quán)值參數(shù)表Table 2 The main evidence layer weights parameterlist of the Sedimentary-metamorphic iron

圖10 冀東沉積變質(zhì)型鐵礦預測后驗概率圖Fig.10 Posterior probability of the Sedimento-metamorphic iron ore deposits in Eastern Hebei 1-地名;2-礦點;3-斷裂;4-后驗概率1-geographic name;2-mineral occurrence;3-fracture;4-posterior probability

致謝 工作過程中得到“華北克拉通前寒武紀重大地質(zhì)事件與成礦”項目組大力支持，在此表示衷心的感謝!

[注釋]

① 河北省地質(zhì)調(diào)查院，河北省冀東鞍山式沉積變質(zhì)型鐵礦典型示范成果報告[R].

Bonham-Carter GF, Agterberg FP, Wright DF. 1989. Weights of evidence modeling：a new approach to mapping mineral potential：In Statistical Applications in the Earth Sciences, Eds[J]. Agterberg FP and Bonham-Carter GF, Geological Survey of Canada, 89(9)：171-183

Cao Zhong-qing. 2004. Geological characteristics and prospecting model of Dajishan niobium-tantalum deposit[J]. Geology and Exploration, 40(6):34-37 (in Chinese with English abstract)

Chen Yong-qing, Wang Shi-cheng. 1995. The Basic principles and methods ore-forming series prognosis of comprehensive information [J]. Shandong Geology, 11(1):55-62(in Chinese with English abstract)

Dai Yan-pei, Zhang Liang-chang, Wang Chang-le, Liu Li, Cui Min-li, Zhu Ming-tian, Xiang Peng. 2012. Genetic type, formation age and tectonic setting of the Waitoushan banded iron formation, Benxi, Liaoning Province[J]. Acta Petrologica Sinica, 28(11)：3574-3594 (in Chinese with English abstract)

Dong Qing-ji, Xiao Ke-yan, Chen Jian-ping, Cong Yuan. 2010. The quantitative analysis of regional metallogenic fault in the northern segment of the Sanjiang metallogenic belt, southwestern China[J]. Geological Bulletin of China, 29(10)：1479-1485 (in Chinese with English abstract)

Dong Ying-Jun, Zhang De-quan, Xu Wen-yi, She Hong-quan. 2005. The mineral resources assessment system and its application in East Kunlun[J]. Acta Geoscientica Sinica, 01:83-88 (in Chinese with English abstract)

Emmanuel John M. Carranza. 2004. Weights of evidence modeling of mineral potential：A case study using small number of prospects, Abra, Philippines[J]. Natural Resources Research, 13(3)：173-187

Liu Hong-tao, Li Jia-sheng， Yan Jian-guo. 2008. Integrated information ore-finding model and direction in the south-central segment of the “Sanjiang melallogenic belt”, Southwestern China[J]. Geology in China, 35(1)：101-110 (in Chinese with English abstract)

Liu Xiao-ling, Sun Yan, Chen Jian-ping,Wang Qian. 2010. Copper minerogenetic prediction based on GIS in Tuotuo River region of Qinghai[J]. Computing Techniques for Geophysical and Geochemical Exploration, 32(4)：408-416 (in Chinese with English abstract)

Qian Xiang-lin, Cui Wen-yuan, Wang Shi-qi, Wang Guan-yu. 1985. Geology of precambrian iron ores in Eastern Hebei province, China[M]. Shijiazhuang：Hebei Science and Technology Press: 1-273 (in Chinese)

Rong Jing-hui, Chen Jian-ping, Shang Bei-chuan. 2012. The prospecting model based on the Yunnan old a deep concealed ore body 3D forecast [J]. Geology and Exploration, 01:191-198 (in Chinese with English abstract)

Sun Yan, Wang Xun-lian, Chen Jian-ping, Wang Qian. 2010. Lead-zinc-silver metallogenic prediction based on GIS in Tuotuo River region, Qinghai, China[J]. Geological Bulletin of China, 29(4)：556-564 (in Chinese with English abstract)

Xu Shan-fa, Chen Jian-ping, Ye Ji-hua. 2006. Application of evidence weight method in the copper-gold mineral resources prediction in the north section of the Sanjiang region[J]. Geology and Prospecting, 42(2)：54-59 (in Chinese with English abstract)

Ye Tian-zhu, Xiao Ke-yan, Yan Guang-sheng. 2007. Deposit model of comprehensive geological information of prediction [J]. Earth Science Frontiers, 05:11-19(in Chinese with English abstract)

Wang Deng-hong，Tang Ju-xing，Ying Li-juan. 2010. Application of “Five levels+Basement” Model for Prospecting Deposits into Depth[J]. Joural of Jilin University: Earth Science Edition，40(4)：733-738

Zhao Peng-da. 2004. Quantitative Geoscience methods and its applications[M]. Beijing: Higher Education Press: 1-30 (in Chinese)

Zhao Peng-da, Wang Jing-gui, Rao Ming-hui. 1995. Geologic anomaly of China[J]. Earth Science, 20 (2) : 117-127 (in Chinese)

[附中文參考文獻]

陳建平，陳 勇，王全明.2008.基于GIS的多元信息成礦預測研究——以赤峰地區(qū)為例[J]. 地學前緣，15(4)：18-26

曹鐘清. 2004. 大吉山鉭鈮鎢礦床地質(zhì)特征及找礦模型[J]. 地質(zhì)與勘探, 40(6)：34-37

陳永清，王世稱. 1995. 綜合信息成礦系列預測的基本原理和方法[J].山東地質(zhì),11(1):55-62

代堰锫,張連昌,王長樂,劉 利,崔敏利,朱明田,相 鵬. 2012. 遼寧本溪歪頭山條帶狀鐵礦的成因類型、形成時代及構(gòu)造背景[J]. 巖石學報, 28(11)：3574-3594

董慶吉,肖克炎,陳建平,叢 源. 2010. 西南“三江”成礦帶北段區(qū)域成礦斷裂信息定量化分析[J].地質(zhì)通報, 29(10)：1479-1485

董英君,張德全,徐文藝,佘宏全. 2005. 礦產(chǎn)資源評價系統(tǒng)及其在東昆侖的應用[J]. 地球?qū)W報,01:83-88

劉洪滔,李家盛,晏建國. 2008. “三江”中南段綜合信息找礦模型及找礦方向[J]. 中國地質(zhì), 35(1)：101-110

劉曉玲,孫巖,陳建平,王 倩. 2010. 基于GIS的沱沱河地區(qū)銅礦綜合信息成礦預測[J]. 物探化探計算技術(shù), 32(4)：408-416

錢祥麟,崔文元,王時麒,王關(guān)玉. 1985. 冀東前寒武紀鐵礦地質(zhì)[M]. 石家莊：河北科學技術(shù)出版社： 1-273

任樹祥,賈正海,張貴寶,宋立軍. 2010. 河北省鐵礦成礦規(guī)律研究及資源潛力預測[J]. 礦床地質(zhì),S1:103-104

戎景會,陳建平,尚北川. 2012. 基于找礦模型的云南個舊某深部隱伏礦體三維預測[J]. 地質(zhì)與勘探,01:191-198

孫巖,王訓練,陳建平,王倩. 2010. 基于證據(jù)權(quán)重法的沱沱河地區(qū)鉛鋅銀礦成礦預測[J].地質(zhì)通報, 29(4)：556-564

王登紅，唐菊興，應立娟. 2010. “五層樓-地下室”找礦模型的適用性及其對深部找礦的意義[J].吉林大學學報：地球科學版,40(4)：733-738

徐善法,陳建平,葉繼華. 2006. 證據(jù)權(quán)法在三江北段銅金礦床成礦預測中的應用研究[J]. 地質(zhì)與勘探, 42(2)：54-59

葉天竺,肖克炎,嚴光生. 2007. 礦床模型綜合地質(zhì)信息預測技術(shù)研究[J]. 地學前緣,05:11-19

趙鵬大. 2004. 定量地學方法及應用[M]. 北京：高等教育出版社： 1-30

趙鵬大, 王京貴, 饒明輝. 1995. 中國地質(zhì)異常[J]. 地球科學—中國地質(zhì)大學學報, 20 (2) : 117-127

Application of Evidence Weight Method Basing on GIS to Metallogenic Prediction in East Hebei province

WANG Jiang-xia1,2,3,CHEN Jian-ping1,2,3*,ZHANG Ying1,2,3,ZHENG Yong-qin1,2,3,CHEN Dong-yue1,2,3

(1. Land resources and high-tech research center of China University of Geosciences,Beijing 1000832. Key Laboratory of Land and resources information development research in Beinjing,Beijing 1000833. School of Earth Science and Resource, China University of Geosciences,Beijing 100083)

Based on the detailed analysis and system summary of the previous research results ,and make full use of the new forecasting theory , mineral resource exploration and evaluation theory and method , we make a systematically arrangement and induction to the information and data etc in the study area. At the same time, combined with GIS technology ,the metallogenic regularity of the main sedimentary metamorphic iron ore was summarized ,also geological information、Geophysical prospecting information and geochemical information related to the ore-forming of the study area was analysed detailed, thus generalizing the prospecting model of diversified information in the region. Then we use the grid cell method on the basis of the method of Evidence weights, circling the metallogenic prospects area and realizing the comprehensive prediction research of the sedimentary metamorphic iron ore resources in Jidong.

East Hebei province,sedimentary metamorphic ore type, multivariate information,metallogenic prediction

2012-05-03；

2012-11-15；[責任編輯]郝情情。

國家重點基礎研究發(fā)展計劃973項目(2012CB416605)資助。

王江霞(1987年-)，女，碩士在讀，研究方向為礦產(chǎn)資源定量預測與評價。E-mail:jiangxia.0709@163.com。

P628+.1

A

0495-5331(2014)03-0464-11