利用證據(jù)權(quán)重法進(jìn)行綜合信息礦產(chǎn)資源預(yù)測(cè)—以大興安嶺北部地區(qū)為例

滕菲，張燕，張素榮，邢怡

（中國(guó)地質(zhì)調(diào)查局天津地質(zhì)調(diào)查中心，天津 300170）

利用證據(jù)權(quán)重法進(jìn)行綜合信息礦產(chǎn)資源預(yù)測(cè)—以大興安嶺北部地區(qū)為例

滕菲，張燕，張素榮，邢怡

（中國(guó)地質(zhì)調(diào)查局天津地質(zhì)調(diào)查中心，天津 300170）

證據(jù)權(quán)重法是一種重要的礦產(chǎn)資源預(yù)測(cè)方法，它是基于二值圖像的地學(xué)統(tǒng)計(jì)方法，通過(guò)對(duì)與礦產(chǎn)資源相關(guān)的地學(xué)信息的疊加分析進(jìn)行預(yù)測(cè)，其優(yōu)點(diǎn)是數(shù)據(jù)驅(qū)動(dòng)、易于編程實(shí)現(xiàn)。筆者通過(guò)對(duì)大興安嶺北部地區(qū)金礦預(yù)測(cè)實(shí)驗(yàn)，證明證據(jù)權(quán)重法數(shù)學(xué)預(yù)測(cè)模型對(duì)區(qū)域成礦遠(yuǎn)景區(qū)預(yù)測(cè)有較高的準(zhǔn)確性和科學(xué)性,同時(shí)在與ARCGIS系統(tǒng)結(jié)合后,可對(duì)多源信息進(jìn)行空間位置綜合處理并形成新的綜合變量,能更全面、客觀的反映地質(zhì)體空間特征,使成礦預(yù)測(cè)更加準(zhǔn)確。

地理信息系統(tǒng);證據(jù)權(quán)重法;成礦預(yù)測(cè);證據(jù)層

隨著大量的多源地學(xué)信息數(shù)據(jù)的積累，以及成礦預(yù)測(cè)學(xué)科自身的不斷發(fā)展，如何適應(yīng)形勢(shì)的變化，充分利用這些豐富的資料，快速、經(jīng)濟(jì)有效地進(jìn)行成礦預(yù)測(cè)，就成了地質(zhì)工作者面臨的一個(gè)重大問(wèn)題。當(dāng)前地理信息系統(tǒng)（GIS）技術(shù)具備了在計(jì)算機(jī)輔助下對(duì)地質(zhì)、地理、地球物理、地球化學(xué)和遙感等多源地學(xué)信息進(jìn)行集成管理、有效綜合與分析的能力，而空間信息的綜合分析方法是決定GIS應(yīng)用水平的關(guān)鍵[1-4]。目前應(yīng)用的方法有很多種，包括布爾邏輯、代數(shù)方法、模糊邏輯神經(jīng)網(wǎng)絡(luò)、秩相關(guān)分析、信息量、判別分析、聚類(lèi)分析和證據(jù)權(quán)重法等[5-7]。證據(jù)權(quán)重法是基于二值圖像的地學(xué)統(tǒng)計(jì)方法，通過(guò)對(duì)與礦產(chǎn)資源相關(guān)的地學(xué)信息的疊加分析進(jìn)行預(yù)測(cè)，其中每一種地學(xué)信息都作為預(yù)測(cè)的一個(gè)證據(jù)因子，而每一個(gè)證據(jù)因子對(duì)礦產(chǎn)資源預(yù)測(cè)的貢獻(xiàn)由其權(quán)重值確定。相對(duì)于其他方法而言，其運(yùn)行由數(shù)據(jù)驅(qū)動(dòng)、亦可專(zhuān)家干預(yù)[8]，對(duì)證據(jù)因子權(quán)重的解釋通俗易懂，其計(jì)算結(jié)果后驗(yàn)概率即為成礦概率，直觀明了，顯示出實(shí)用、易用等優(yōu)勢(shì)，是國(guó)際地學(xué)領(lǐng)域使用GIS技術(shù)進(jìn)行礦產(chǎn)資源評(píng)價(jià)預(yù)測(cè)一種非常受重視的計(jì)算方法。筆者以大興安嶺北部地區(qū)為例，對(duì)證據(jù)權(quán)重法的特點(diǎn)方法、實(shí)現(xiàn)過(guò)程等進(jìn)行了詳細(xì)探討，希望能為區(qū)域礦產(chǎn)預(yù)測(cè)工作提供可借鑒的方法經(jīng)驗(yàn)。

1 證據(jù)權(quán)重法

1.1 方法原理

證據(jù)權(quán)重法是加拿大數(shù)學(xué)地質(zhì)學(xué)家Agterberg提出的一種地學(xué)統(tǒng)計(jì)方法，通過(guò)對(duì)一些與礦產(chǎn)形成相關(guān)的地學(xué)信息圖層的加權(quán)疊加復(fù)合來(lái)進(jìn)行礦產(chǎn)遠(yuǎn)景區(qū)的預(yù)測(cè)。其基礎(chǔ)是將地質(zhì)評(píng)價(jià)模型轉(zhuǎn)換為“網(wǎng)格模型”，針對(duì)每一網(wǎng)格上信息數(shù)據(jù)進(jìn)行權(quán)值計(jì)算，結(jié)合GIS技術(shù)對(duì)各種有利成礦因素（證據(jù)層）進(jìn)行有效綜合的數(shù)學(xué)預(yù)測(cè)評(píng)價(jià)。只要所輸入的圖層合理，選擇條件適當(dāng)，計(jì)算結(jié)果就客觀可信[9]。

數(shù)學(xué)原理如下：假設(shè)研究區(qū)被劃分成為T(mén)個(gè)等大小的網(wǎng)格（區(qū)域劃分網(wǎng)格單元要以保證一個(gè)網(wǎng)格中至多出現(xiàn)一個(gè)已知礦點(diǎn)為基礎(chǔ)），提取相關(guān)信息。其中有D個(gè)網(wǎng)格單元含有已知礦（床）點(diǎn)，提取的每個(gè)變量看作獨(dú)立的證據(jù)因子（以下稱(chēng)證據(jù)因子）。對(duì)于任意一個(gè)證據(jù)因子，其對(duì)成礦貢獻(xiàn)的權(quán)重定義為：分別為證據(jù)因子存在區(qū)和不存在區(qū)的權(quán)重值，如果原始數(shù)據(jù)缺失，那么權(quán)重值為0；其中B為這個(gè)因子存在的網(wǎng)格單元總數(shù)，-B為這種因子不存在的網(wǎng)格單元總數(shù)。-B=T-B;-D為不含礦單元數(shù)；條件概率：

利用證據(jù)權(quán)重法計(jì)算出各個(gè)變量的W＋和W－值；控礦地質(zhì)因素與礦床產(chǎn)出狀態(tài)之間的關(guān)聯(lián)性強(qiáng)弱，可以通過(guò)正負(fù)權(quán)的差值大小來(lái)度量，即：Cj=W+j+W-j, Cj值越大表示該地質(zhì)標(biāo)志的找礦指示性越好,若Cj= 0，表示該找礦標(biāo)志對(duì)有礦與無(wú)礦無(wú)指示意義；Cj＞0表示有利于成礦，Cj＜0表示不利于成礦。

對(duì)于n個(gè)證據(jù)因子，若它們都關(guān)于礦點(diǎn)條件獨(dú)立，則研究區(qū)任一網(wǎng)格單元k為礦點(diǎn)的可能性，用后驗(yàn)概率對(duì)數(shù)表示為：

根據(jù)該公式，后驗(yàn)概率為:P=O/(1+O)。后驗(yàn)概率P代表了單個(gè)單元格內(nèi)找礦的有利度，根據(jù)后驗(yàn)概率提出找礦遠(yuǎn)景區(qū)。

1.2 方法優(yōu)越性

證據(jù)權(quán)重法礦產(chǎn)資源評(píng)價(jià)相對(duì)于其它基于GIS的礦產(chǎn)資源評(píng)價(jià)方法相比，還具有如下優(yōu)點(diǎn)：易于編程，易于組合多圖件模式；預(yù)測(cè)過(guò)程中以礦產(chǎn)產(chǎn)出的空間位置為基礎(chǔ)，并在空間位置最有效、最直接的表達(dá)形式，即可在圖形上進(jìn)行操作；對(duì)證據(jù)因子權(quán)重的解釋通俗易懂，正負(fù)權(quán)重的差值直接指示了證據(jù)因子（找礦標(biāo)志）與礦的相關(guān)程度；證據(jù)因子的處理簡(jiǎn)單，只需離散到二值圖像即可，調(diào)整閥值使正負(fù)權(quán)重的差值C達(dá)最大，則用該閥值離散成二值圖像為最佳；能處理缺失信息的區(qū)域，只需設(shè)其權(quán)重為0即可，不需對(duì)數(shù)據(jù)空白區(qū)進(jìn)行人為的插值；計(jì)算結(jié)果是后驗(yàn)概率，代表了單個(gè)單元格內(nèi)找礦的成礦概率，直觀明了；可采用單純數(shù)據(jù)驅(qū)動(dòng)方法運(yùn)行，也可采用人工干涉，即專(zhuān)家證據(jù)權(quán)重法，在運(yùn)行過(guò)程中加入經(jīng)驗(yàn)、認(rèn)識(shí)，綜合預(yù)測(cè)，驅(qū)動(dòng)方式靈活[10-12]；可實(shí)現(xiàn)基于GIS的證據(jù)權(quán)重法礦產(chǎn)資源評(píng)價(jià)的軟件豐富。

2 區(qū)域地質(zhì)概況

研究區(qū)處于古西伯利亞板塊南緣，額爾古納斷褶帶、內(nèi)蒙—大興安嶺褶皺系和上黑龍江斷陷的接合部位（圖1）?；鹕交顒?dòng)強(qiáng)烈，被燕山期火山巖大片覆蓋。區(qū)內(nèi)地層分布廣泛但比較零星，除中上寒武統(tǒng)，中石炭統(tǒng)外，自元古界至第四系全新統(tǒng)均有出現(xiàn)，以中、新生界最為發(fā)育。區(qū)內(nèi)巖漿活動(dòng)十分頻繁，具有多期次的特點(diǎn)。侵入巖以燕山期中酸性、酸性花崗巖類(lèi)最發(fā)育，分布南北存在差異性，南部分布廣泛，北部分布相對(duì)很少。區(qū)內(nèi)礦產(chǎn)資源豐富，以金、銅、鉛、鋅等礦床為主。

圖1 研究區(qū)大地構(gòu)造位置Fig.1 Geotectonic position o f the study area

3 區(qū)域找礦模型

筆者在分析了研究區(qū)地層、構(gòu)造、巖漿巖等地質(zhì)資料的基礎(chǔ)上,對(duì)該區(qū)金礦床的地質(zhì)特征和控礦因素進(jìn)行了分析和總結(jié),建立了金礦區(qū)域找礦模型[13]-[17]。

（1）賦礦地層：二十二站組、甘河組、倭勒根群大網(wǎng)子組、興華渡口期片麻巖、修峰組、塔木蘭溝組、光華組。

（2）控礦巖漿巖（侵入巖）：主要以燕山期中酸性侵入巖為主，其它時(shí)期中酸性侵入巖次要。

（3）主要控礦構(gòu)造：NE和NW及WE向。

（4）指示有礦的重力異常以密度異常區(qū)為主，即主要為高密度體內(nèi)的低密度體或低密度體內(nèi)的高密度體。

（5）指示有礦的航磁異常為正異常。

（6）與礦化有關(guān)的地球化學(xué)異常主要有Au、Cu、Pb、Zn。

4 預(yù)測(cè)流程

4.1 預(yù)測(cè)要素的提取

（1）建立研究區(qū)礦床(點(diǎn))和成礦作用有關(guān)的證據(jù)因子專(zhuān)題數(shù)據(jù)庫(kù)。

分別建立地質(zhì)數(shù)據(jù)庫(kù)（包括地層、構(gòu)造、巖漿巖、礦床(點(diǎn))數(shù)據(jù)）、地球化學(xué)數(shù)據(jù)庫(kù)（被確認(rèn)與金礦有密

切共生及伴生關(guān)系的Au、Cu、Pb、Zn等元素）、金礦床(點(diǎn))數(shù)據(jù)庫(kù)、遙感數(shù)據(jù)庫(kù)（利用TM遙感圖像對(duì)區(qū)內(nèi)的線環(huán)構(gòu)造進(jìn)行校正、解譯）、地球物理數(shù)據(jù)庫(kù)（重力、磁力異常），作為基本證據(jù)層統(tǒng)一投影存儲(chǔ)在MAPGIS、ARCGIS地圖庫(kù)中。

（2）對(duì)于每一個(gè)證據(jù)因子專(zhuān)題數(shù)據(jù)庫(kù)，提取與試驗(yàn)區(qū)成礦關(guān)系密切的控礦因素，生成新的圖層，即證據(jù)層。以地質(zhì)數(shù)據(jù)庫(kù)為例：根據(jù)所建立金礦區(qū)域地質(zhì)找礦模型，針對(duì)可能與礦產(chǎn)相關(guān)地質(zhì)要素，在MAPGIS平臺(tái)下，初步劃分出57個(gè)地質(zhì)圖層(表1)，其中控礦巖漿巖主要以燕山期中酸性侵入巖為主，控礦構(gòu)造以北東向、北西、東西向?yàn)橹鳎c成礦關(guān)系密切地層主要有中生界的侏羅系和白堊系甘河組等6組及興華渡口期花崗片麻巖，其它成礦相關(guān)圖層還包括火山巖、潛火山巖、糜棱巖、脈巖、礦產(chǎn)地等。

（3）提取礦點(diǎn)圖層，劃分網(wǎng)格單元

單元格面積的設(shè)置要保證一個(gè)單元格中最多只出現(xiàn)一個(gè)礦床（點(diǎn)）（訓(xùn)練點(diǎn)）。故單元格面積的確定是根據(jù)已知礦床（點(diǎn)）的多少?zèng)Q定的。實(shí)驗(yàn)區(qū)所有的金礦床（點(diǎn)）11個(gè)，將礦床（點(diǎn)）的文件轉(zhuǎn)換為Arcview下的shape格式，并將建立模型的單元格面積設(shè)置為1 km2，搜索含礦網(wǎng)格單元，生成含礦網(wǎng)格圖層。

4.2 預(yù)測(cè)模型的建立

（1）證據(jù)因子的選擇、定量化

地質(zhì)證據(jù)因子：將所提取賦礦地層、控礦巖漿巖、控礦斑巖、脈巖轉(zhuǎn)換成grid格式，考慮到控礦巖漿巖、控礦斑巖、脈巖等熱液活動(dòng)有一定影響范圍，利用A rc-Wofe的緩沖分析功能，分別進(jìn)行緩沖分析，緩沖半徑以1 km為間隔，并計(jì)算各緩沖半徑的C（正負(fù)權(quán)重的差值），燕山期中酸性侵入巖最大的C（1.0014）緩沖半徑為7 km，因此選緩沖半徑7 km作為其影響區(qū)域。同理，選緩沖半徑5 km作為控礦斑巖、脈巖的影響區(qū)域。下面控礦構(gòu)造以此為例，介紹相關(guān)地質(zhì)要素圖層緩沖半徑計(jì)算過(guò)程?？氐V構(gòu)造是線文件，利用MAPGIS的空間檢索功能檢索出構(gòu)造線方向?yàn)镹E和NW及EW向的所有斷裂構(gòu)造，從遙感影像解譯的線、環(huán)構(gòu)造與原有構(gòu)造合為一個(gè)證據(jù)層。再將其轉(zhuǎn)換為shape格式，利用Arc-Wofe的功能，對(duì)其進(jìn)行緩沖分析，緩沖半徑以1 km為間隔，并計(jì)算各緩沖半徑的C（表2），最大的C（0.591）緩沖半徑為8 km。

重磁異常證據(jù)因子：從區(qū)域找礦模型來(lái)看，局部重磁異常與礦關(guān)系密切，故對(duì)其進(jìn)行人工提取，作為一個(gè)證據(jù)因子，并轉(zhuǎn)換為grid格式。根據(jù)重力異常等值線圖，按10毫伽為間隔計(jì)算各重力等值線為閥值的異常區(qū)與礦相關(guān)程度C值，結(jié)果顯示當(dāng)?shù)戎稻€為-30毫伽時(shí)，C值獲得最大值3.9373（表3），故用此等值線為閥值，生成重力異常二值圖像并轉(zhuǎn)換成grid格式。根據(jù)航磁異常等值線圖，按100 nT為間隔計(jì)算各重力等值線為閥值的異常區(qū)與礦相關(guān)程度C值，結(jié)果顯示當(dāng)?shù)戎稻€為700 nT時(shí)，C值獲得最大值2.0628，故用此等值線為閥值，生成重力異常二值圖

像并轉(zhuǎn)換成grid格式。

表1 地質(zhì)證據(jù)層分類(lèi)表Table 1 Geo logica levidence layers system

表2 斷裂構(gòu)造與金礦床(點(diǎn))相關(guān)程度Table 2 Pertinency between fau lts and go ld deposits

化探異常證據(jù)因子：依據(jù)區(qū)域找礦模型，將研究區(qū)Cu、Au、Pb、Zn元素化探異常分別作為一個(gè)證據(jù)因子，并轉(zhuǎn)換為grid格式?；疆惓Ｗ鳛榈乇碓馗患谋硐?對(duì)成礦遠(yuǎn)景區(qū)的圈定具有直接指示意義。對(duì)化探異常下限的設(shè)定，同樣采用C值最大來(lái)進(jìn)行優(yōu)選[18]。

表3 布格重力異常與金礦床（點(diǎn)）相關(guān)程度Table 3 Pertinency between bouguer gravity anom a lies and go ld deposits

（2）金礦證據(jù)因子的條件獨(dú)立性檢驗(yàn)及其合并或篩選

利用Arc-Wofe條件獨(dú)立性檢驗(yàn)功能，對(duì)上述所有證據(jù)因子進(jìn)行關(guān)于礦床（點(diǎn)）是否為條件獨(dú)立進(jìn)行檢驗(yàn)，結(jié)果有多個(gè)證據(jù)因子之間不存在條件獨(dú)立，經(jīng)過(guò)反復(fù)對(duì)各證據(jù)因子的條件獨(dú)立性檢驗(yàn)，剔除或合并條件獨(dú)立的證據(jù)因子，最后選取如下12個(gè)證據(jù)因子建立金礦證據(jù)權(quán)重模型，它們之間關(guān)于已知11個(gè)礦床（點(diǎn)）都是條件獨(dú)立的，也即通過(guò)了條件獨(dú)立性檢驗(yàn)（表4）。概率值＜0.05表示對(duì)應(yīng)的證據(jù)因子之間不存在條件獨(dú)立。

表4 證據(jù)因子條件獨(dú)立性檢驗(yàn)結(jié)果表Table 4 Testing resu lt o f cond itiona l dependency o f evidence-gene

4.3 預(yù)測(cè)結(jié)果及評(píng)價(jià)

（1）證據(jù)因子正負(fù)權(quán)重計(jì)算

在選定的證據(jù)因子通過(guò)條件獨(dú)立性檢驗(yàn)后，計(jì)算它們的正負(fù)權(quán)重，其計(jì)算結(jié)果見(jiàn)表5。根據(jù)其結(jié)果進(jìn)行分析可得出以下認(rèn)識(shí):本實(shí)驗(yàn)區(qū)侏羅紀(jì)和早白堊世地層與成礦關(guān)系密切;燕山期的中酸性巖漿巖在本區(qū)所占面積不大,但與金礦化的關(guān)系比較密切;構(gòu)造破碎帶與遙感解譯的線性構(gòu)造與金礦化的關(guān)系都比較密切;航磁、重力異常與金礦化也有一定的關(guān)系。

（2）后驗(yàn)概率計(jì)算

筆者根據(jù)每種條件的權(quán)重值計(jì)算得到實(shí)驗(yàn)區(qū)的后驗(yàn)概率等值線圖。并根據(jù)后驗(yàn)概率值劃分3級(jí)有利預(yù)測(cè)成礦區(qū)（圖2）:把值域范圍（0.0046～0.0295）的預(yù)測(cè)單元作為Ⅰ級(jí)遠(yuǎn)景區(qū)，值域范圍（0.0015～0.0046)的預(yù)測(cè)單元作為Ⅱ級(jí)遠(yuǎn)景區(qū)，值域范圍(0.0005～0.0015)的預(yù)測(cè)單元作為III級(jí)遠(yuǎn)景區(qū)。分析表明已知礦點(diǎn)基本落入后驗(yàn)概率出現(xiàn)的Ⅰ、Ⅱ級(jí)分布范圍內(nèi)，說(shuō)明預(yù)測(cè)結(jié)果對(duì)已知條件沒(méi)有出現(xiàn)漏判；預(yù)測(cè)成礦遠(yuǎn)景區(qū)分布走向基本與區(qū)域內(nèi)的控礦構(gòu)造、巖漿巖分布一致，成礦作用明顯受NE和NW向的斷裂控制，預(yù)測(cè)好的區(qū)域基本都在NE和NW向斷裂的交匯處，出露地層均為實(shí)驗(yàn)區(qū)內(nèi)有利含礦地層，表現(xiàn)出成礦受巖漿活動(dòng)和層位以及斷裂的控制，這與區(qū)域內(nèi)已知金礦床特征相吻合；除了覆蓋已知

礦床（點(diǎn)）的區(qū)域外還有幾個(gè)預(yù)測(cè)區(qū)尚未發(fā)現(xiàn)金礦床（點(diǎn)）。研究分析表明這一帶都位于較大斷裂帶的交匯帶上和大型巖體出露附近，并與已知礦床(點(diǎn))具有相似的成礦條件，所以這一帶尚未發(fā)現(xiàn)礦床（點(diǎn)）的預(yù)測(cè)區(qū)應(yīng)該在下一步找礦勘探工作中引起足夠的重視。

表5 證據(jù)因子正負(fù)權(quán)重參數(shù)表Table 5 W eigh t pa ram e te rs o f evidence-gene

圖2 大興安嶺北部地區(qū)金成礦預(yù)測(cè)圖Fig.2 Go ld p rognosis m ap o f No rthDaxingan ling

5 結(jié)論

證據(jù)權(quán)重法預(yù)測(cè)模型是根據(jù)已知礦床（點(diǎn)）與各種控礦成礦條件之間的條件概率來(lái)確定每種條件的權(quán)重值，然后推廣到全區(qū)，其原理簡(jiǎn)單而且易于實(shí)現(xiàn)，既可以對(duì)多金屬礦床進(jìn)行預(yù)測(cè)，也可以對(duì)特定成因礦種進(jìn)行綜合預(yù)測(cè)。由于后驗(yàn)概率是在大量地質(zhì)、地球物理和地球化學(xué)等圖層疊加操作的基礎(chǔ)上計(jì)算出來(lái)的，因此，其結(jié)果綜合反映了各種控礦因素和礦化信息對(duì)礦床的控制和指示意義。同時(shí)由于證據(jù)權(quán)重法著重對(duì)各控礦要素和礦床（點(diǎn)）間相互空間位置進(jìn)行分析，其和相關(guān)GIS軟件的結(jié)合滿(mǎn)足大量的空間疊加、復(fù)合分析的需求。

筆者在大興安嶺北部地區(qū)建立的基于GIS的證據(jù)權(quán)重法礦產(chǎn)資源評(píng)價(jià)模型，反映該方法能夠在研究程度相對(duì)較低，但數(shù)據(jù)資料豐富的地區(qū)快速而有效應(yīng)用的特點(diǎn)，可供其它地區(qū)礦產(chǎn)資源評(píng)價(jià)工作借鑒。筆者認(rèn)為證據(jù)權(quán)重法礦產(chǎn)資源的評(píng)價(jià)工作，既可以在專(zhuān)家的指導(dǎo)下，用其實(shí)現(xiàn)專(zhuān)家的思想，也可以發(fā)揮證據(jù)權(quán)重法數(shù)據(jù)驅(qū)動(dòng)對(duì)豐富資料充分利用的功能，作為專(zhuān)家思想的補(bǔ)充或旁證，為地質(zhì)找礦提供靈活、有效的預(yù)測(cè)工具，實(shí)現(xiàn)隱伏礦床的找礦突破。

[1]朱思才，吳家齊，劉和發(fā).GIS技術(shù)在區(qū)域礦產(chǎn)資源勘查評(píng)價(jià)中的應(yīng)用[J].有色金屬礦產(chǎn)與勘查，1999，8(6):615-618.

[2]向運(yùn)川，任天祥，楊竹溪.開(kāi)發(fā)利用地理信息系統(tǒng)(GIS)綜合分析地學(xué)信息進(jìn)行礦產(chǎn)預(yù)測(cè)[J].物探與化探，1996，20(1):1-13.

[3]黃海峰.GIS在成礦預(yù)測(cè)中的應(yīng)用[J].甘肅地質(zhì)學(xué)報(bào)2002，11(1):89-96.

[4]張建龍，游麗君，秦舉禮.利用GIS空間分析功能進(jìn)行礦產(chǎn)資源評(píng)價(jià)的方法研究[J].物探化探計(jì)算技術(shù)，1997，19(3): 283-286.

[5]葛詠.地學(xué)數(shù)據(jù)集成及空間決策支持的方法與應(yīng)用[J].地球信息科學(xué)，2006，8(1):16-20.

[6]肖克炎，葛艷.成礦預(yù)測(cè)中證據(jù)權(quán)重法與信息量法及其比較[J].物探化探計(jì)算技術(shù),1999,21(3):223-226.

[7]張健挺，邱友良.人工智能和專(zhuān)家系統(tǒng)在地學(xué)中的應(yīng)用綜述[J].地理科學(xué)進(jìn)展，1998，17(1):44-51.

[8]丁清峰，孫豐月.專(zhuān)家證據(jù)權(quán)重法及其在東昆侖地區(qū)的應(yīng)用[J].地質(zhì)與勘探2005，41(4):88-94.

[9]袁峰，周濤發(fā)，岳書(shū)倉(cāng).GIS礦產(chǎn)資源預(yù)測(cè)中的證據(jù)權(quán)重法[J].黃金地質(zhì)，2003，9(3):75-77.

[10]黃海峰，姚書(shū)振，丁振舉.基于GIS的證據(jù)權(quán)重法在成礦預(yù)測(cè)中的應(yīng)用—以甘肅省岷縣—禮縣地區(qū)的金礦預(yù)測(cè)為例[J].地質(zhì)科技情報(bào)，2003，22(3):77-82.

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

[12]楊茂森，黎清華，楊海巍.GIS支持下的證據(jù)權(quán)法在膠東金礦集中區(qū)預(yù)測(cè)中的應(yīng)用[J].地球?qū)W報(bào)，2005，26(5): 487-491.

[13]鄒光華，歐陽(yáng)宗圻，李慧，等.中國(guó)主要類(lèi)型金礦床找礦模型[M].北京:地質(zhì)出版社，1996.

[14]王鐘，邵孟林，肖樹(shù)建.隱伏有色金屬礦床綜合找礦模型[M].北京:地質(zhì)出版社，1996.

[15]趙鵬大，孟憲國(guó).地質(zhì)異常與礦產(chǎn)預(yù)測(cè)[J].地球科學(xué)-中國(guó)地質(zhì)大學(xué)學(xué)報(bào)，1993，18(1):39-47.

[16]趙鵬大，胡旺亮，李紫金.礦床統(tǒng)計(jì)預(yù)測(cè)[M].北京:地質(zhì)出版社，1983.

[17]王世稱(chēng)，成秋明，范繼璋.金礦資源綜合信息評(píng)價(jià)方法[M].長(zhǎng)春:吉林科學(xué)技術(shù)出版社，1990.

[18]丁汝福.國(guó)內(nèi)外尋找隱伏礦化探新方法研究進(jìn)展[J].地質(zhì)與勘探，1999，3，35(3).

P628+.2

：A

：1672－4135（2014）04－0269-05

2014-07-31

中國(guó)地質(zhì)調(diào)查項(xiàng)目：區(qū)域地球物理成果集成與方法技術(shù)研究（1212011120915）

滕菲（1982-），女，工程師，2007年7月畢業(yè)于吉林大學(xué)數(shù)字地質(zhì)科學(xué)專(zhuān)業(yè)，現(xiàn)主要從事應(yīng)用地球物理研究，E-mail:tengbeibei2006@163.com。