基于GIS的證據(jù)權(quán)重法在桂西南地區(qū)優(yōu)質(zhì)錳礦成礦預(yù)測中的應(yīng)用

張寶一，毛先成，周尚國，胡 超，閆 芳

（1．中南大學(xué)有色金屬成礦預(yù)測教育部重點(diǎn)實(shí)驗(yàn)室，地球科學(xué)與信息物理學(xué)院，長沙410083；2．中國冶金地質(zhì)總局，北京100025）

0 引言

廣義上說，“優(yōu)質(zhì)錳礦”就是低鐵（Mn／Fe≥4～6）、低磷（P／Mn≤0．003～0．005）的工業(yè)錳礦石，其內(nèi)核是按工業(yè)利用要求限定了礦石中雜質(zhì)磷和鐵的含量。優(yōu)質(zhì)錳礦已被列為國家戰(zhàn)略資源和緊缺礦種的序列，是近期地質(zhì)勘查的目標(biāo)礦種之一。駱華寶［1］從優(yōu)質(zhì)錳礦成礦地質(zhì)背景、成礦條件、地質(zhì)工作程度和找礦潛力分析，將桂西南地區(qū)作為近期優(yōu)質(zhì)錳礦勘查的5個(gè)重要遠(yuǎn)景區(qū)之一，區(qū)內(nèi)產(chǎn)有亞洲最大的錳礦——廣西下雷超大型錳礦。

桂西南地區(qū)地質(zhì)構(gòu)造上屬右江盆地，含錳建造為上泥盆統(tǒng)榴江組硅質(zhì)巖和五指山組硅質(zhì)、泥質(zhì)碳酸鹽巖建造，下石炭統(tǒng)大塘組硅質(zhì)巖、泥質(zhì)碳酸鹽巖建造和下三疊統(tǒng)北泗組碳酸鹽巖建造［1］。這3個(gè)建造自南向北分布，構(gòu)成NE－SW向延伸的3個(gè)礦帶，礦源層經(jīng)過氧化、淋積、遷聚往往形成大礦。優(yōu)質(zhì)錳礦床的主要類型有：產(chǎn)于上泥盆統(tǒng)五指山組含錳硅質(zhì)巖的“下雷式”熱水沉積－次生氧化型優(yōu)質(zhì)錳礦床；產(chǎn)于下三疊統(tǒng)北泗組含錳泥質(zhì)巖的“東平式”淋積型優(yōu)質(zhì)錳礦床；產(chǎn)于下石炭統(tǒng)大塘組含錳碳酸鹽巖的“寧干式”正常沉積－次生氧化型優(yōu)質(zhì)錳礦床。

1 優(yōu)質(zhì)錳礦形成條件

優(yōu)質(zhì)錳礦是錳在沉積環(huán)境中在走向上和層（時(shí)）序上發(fā)生錳、磷、鐵的分離而形成的，分離的關(guān)鍵因素是pH值和Eh值；錳質(zhì)主要在成錳盆地中蘊(yùn)集成礦，錳礦富集在層序地層格架中的特定部位［1］。

桂西南地區(qū)錳礦的形成經(jīng)歷了原生沉積、構(gòu)造重建及熱液疊加改造、表生風(fēng)化3個(gè)階段，各階段物理化學(xué)環(huán)境的改變都可發(fā)生錳、鐵、磷分離及錳富集，并形成優(yōu)質(zhì)（富）錳礦。何知禮［2－3］從地洼學(xué)說 的角度分析，地臺期形成的礦源層總厚度較大，加上地洼期強(qiáng)烈的熱液疊加改造并帶來部分錳質(zhì)，對于形成優(yōu)質(zhì)復(fù)成菱錳礦－硫錳礦礦床具有決定意義；對于次生氧化優(yōu)質(zhì)錳礦床，則礦源層的總厚度、錳含量以及風(fēng)化富集的程度特別重要。因此，不論原生沉積優(yōu)質(zhì)錳礦還是次生氧化優(yōu)質(zhì)錳礦，含錳巖系或錳礦層均受到同沉積走滑拉斷盆地、地層層位、沉積相位及其微相組合類型、后期剝離斷層及分劃性順層韌性走滑斷層、構(gòu)造熱液和巖漿熱液及表生風(fēng)化等多種地質(zhì)因素的控制。

1．1 原生含錳巖系沉積古地理?xiàng)l件

桂西南地區(qū)的原生沉積優(yōu)質(zhì)錳礦或次生氧化優(yōu)質(zhì)錳礦的“礦胚層”均屬于海相沉積型碳酸錳礦床。沉積礦床的形成是一個(gè)復(fù)雜的成礦地質(zhì)過程，必須具備一系列待定的成礦地質(zhì)條件，諸如構(gòu)造狀況、地層特點(diǎn)、沉積相和古地理環(huán)境、物質(zhì)來源等等［4］。海相錳礦的成礦作用過程是在特定的海洋－環(huán)大陸地質(zhì)構(gòu)造、水圈和氣候帶演化的動態(tài)體制中，受各種地質(zhì)地球化學(xué)因素相互制約（配合）的條件下進(jìn)行的［5］。因此，分析研究錳與磷、鐵分異并形成優(yōu)質(zhì)錳礦的巖相古地理環(huán)境是優(yōu)質(zhì)碳酸錳礦正確評價(jià)的前提。

（1）古地理?xiàng)l件。被動大陸邊緣、板塊離散和會聚環(huán)境下的成錳盆地控制錳礦帶，成錳盆地的封閉條件是錳質(zhì)富集成工業(yè)礦床的一個(gè)重要因素。一般來說，成錳盆地的范圍不大，常為臺地間狹長的溝槽和其間的拉斷盆地或緊靠臺地邊緣的拗陷，水體較深，具有特殊的地球化學(xué)和物理化學(xué)環(huán)境，有利于錳質(zhì)的沉積和富集。

（2）同沉積構(gòu)造。區(qū)域性基底斷裂和同沉積斷裂控制的臺緣和臺間盆地控制著錳礦床的形成和分布［1］。桂西南地區(qū)晚泥盆世存在2條平行的NE向同沉積走滑斷裂（下雷—東平斷裂和龍邦—地州斷裂），是控制晚泥盆世錳礦床（點(diǎn)）的同生走滑斷裂構(gòu)造。桂西南地區(qū)早三疊世晚期的弧后盆地前期斷裂拗陷也是錳礦形成的有利空間，控制著早三疊世晚期錳礦的形成和分布。

（3）沉積相條件。錳礦層一般賦存在半封閉、欠補(bǔ)償?shù)某练e相帶中。桂西南地區(qū)晚泥盆世為碳酸鹽臺地相區(qū)，可分為半局限臺地和臺溝相帶，分別形成半局限臺地相碳酸鹽巖沉積和臺溝相泥質(zhì)、硅質(zhì)、灰質(zhì)巖沉積，曾友寅［6－7］在對下雷層狀錳礦床噴氣熱液沉積的特征和成因研究中發(fā)現(xiàn)，熱液沉積錳礦層只產(chǎn)于淺海深水臺溝沉積相中，而且主要集中分布于下雷—東平臺溝相帶西南部的拉斷盆地（溝底）亞相中。早石炭世大塘期時(shí)桂西南的印茶—愛屯沉積盆地沉積了灰?guī)r、含燧石灰?guī)r、泥質(zhì)灰?guī)r、硅質(zhì)巖、含錳碳酸鹽巖。早三疊世晚期，桂西南東平地區(qū)處于田東陸棚相帶邊緣，沉積環(huán)境為水體較深、寧靜的淺海較深水域，封閉條件良好，錳質(zhì)易于聚集，又離古陸剝蝕區(qū)不遠(yuǎn)，海相火山活動強(qiáng)烈，形成了大型海相沉積型貧碳酸錳沉積礦床［8］。

（4）地層層位。地層與沉積錳礦的關(guān)系反映在地層時(shí)代與錳礦成礦作用的關(guān)系、地層產(chǎn)狀和空間分布對錳礦層的控制、巖性及巖石組合與錳礦層的關(guān)系等方面。桂西南地區(qū)主要成錳時(shí)代為晚泥盆世、早石炭世大塘期和早三疊世晚期。晚泥盆世沉積錳礦主要與硅質(zhì)巖－碳酸鹽巖組合有關(guān)；下石炭統(tǒng)大塘階錳礦一般產(chǎn)于含硅質(zhì)的碳酸鹽巖或硅質(zhì)頁巖中；早三疊世錳礦夾于含硅質(zhì)的泥灰?guī)r中。雖然各成錳時(shí)代地層的巖石和巖石組合不盡相同，但就其巖石組分和巖石組合成分來說，均可納入通常所說的硅、灰、泥組合系列［4］。

1．2 構(gòu)造重建及熱液疊加條件

一般來說，構(gòu)造重建及熱液改造對于優(yōu)質(zhì)錳礦的形成沒有什么特別明顯的意義，但是剝離斷層、分劃性壓剪切構(gòu)造帶可以使原生沉積錳礦中的元素活化、遷移，構(gòu)造－巖漿活動帶中的侵入體和熱液流體也對原生沉積錳礦有疊加改造作用。

（1）剝離斷層。剝離斷層是伸展構(gòu)造在脆－韌性過渡域中的主要表現(xiàn)形式，在桂西南地區(qū)的錳質(zhì)活化－遷移、疊加富集和蝕變過程中起了重要的動力驅(qū)動作用。由于主剝離斷層的影響，在含錳巖系和普通錳礦層內(nèi)產(chǎn)生一系列次級剝離斷層，斷層碎裂帶為含礦液的運(yùn)移、充填交代提供了良好的空間。

（2）分劃性壓剪切構(gòu)造帶（S面理帶）。斷裂帶中元素的遷移及其產(chǎn)生的化學(xué)變化（演化）是在構(gòu)造動力作用下同步產(chǎn)生的，離子半徑大的元素和一些離子半徑小而不穩(wěn)定的元素向壓剪切面以外運(yùn)移，而中等離子半徑的元素（如Mn，Si，F(xiàn)e）則向壓剪切面內(nèi)運(yùn)移，并在分劃性壓剪切帶附近富集。

（3）熱液疊加條件。在某些中生代構(gòu)造－巖漿活動帶中，侵入體和熱液流體攜帶的化學(xué)組分參與原生沉積錳礦床的疊加改造作用，并形成一些與多金屬礦床共生的“層控型”或沿?cái)嗔褞Ы淮?、充填氧化錳－碳酸錳－硅酸錳礦床［5］。

1．3 次生富集氧化條件

在次生氧化錳礦成礦作用過程中，元素的化學(xué)分異和Fe，Mn的遷移富集，乃至礦體的空間定位形 式，都與發(fā)生在表生帶中的2種化學(xué)風(fēng)化機(jī)制（減積機(jī)制和加積機(jī)制）有關(guān)，并受礦源層性質(zhì)、氣候，以及與區(qū)域新構(gòu)造運(yùn)動穩(wěn)定程度有關(guān)的地形地貌和排水條件等因素的控制［5］。

（1）礦源層條件。次生氧化錳礦床必須以原生沉積碳酸錳礦床作為礦源層，其礦床規(guī)模及礦石質(zhì)量受原生碳酸錳礦床的規(guī)模和質(zhì)量控制。

（2）地質(zhì)構(gòu)造條件。桂西南地區(qū)次級褶皺構(gòu)造發(fā)育，使含錳層在地表反復(fù)褶曲，造成埋深淺、分布廣的有利構(gòu)造條件，易形成大型氧化錳礦床，含錳礦床多分布在復(fù)式背向斜構(gòu)造的兩翼［8］。

（3）地形地貌條件。當(dāng)層狀原生礦的產(chǎn)狀與地形坡度一致，上覆厚度不大時(shí)，往往有利于錳帽型氧化錳的形成，礦體在平面上呈大面積分布，形成所謂“鋪山礦”。毛先成等［9］也發(fā)現(xiàn)坡度小、起伏度小、切割微弱的地帶有利于錳礦物的次生氧化富集。

（4）近代地殼運(yùn)動及水文條件。當(dāng)?shù)貧ど仙俣容^慢且與潛水面的下降速度大致相等時(shí)，潛水面較低的礦床氧化帶的深度及礦石的氧化程度較大，若原生礦石品位較高有利于形成較富的氧化錳。

（5）氣候條件。我國北緯23°帶內(nèi)是氧化錳礦分布的主要地帶，長期濕熱的氣候及充沛的雨水是加速碳酸錳分解形成氧化富錳礦的重要條件。

（6）微生物條件。一方面微生物可以直接對錳進(jìn)行氧化，使原生碳酸錳礦轉(zhuǎn)變?yōu)檠趸i礦；另一方面微生物又可以改變環(huán)境的pH值，促進(jìn)錳的化學(xué)氧化［11］。

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

證據(jù)權(quán)重法是加拿大數(shù)學(xué)地質(zhì)學(xué)家F．P．Agterberg提出的一種地學(xué)統(tǒng)計(jì)方法。它采用一種統(tǒng)計(jì)分析模式，通過對一些與成礦相關(guān)的地學(xué)信息的疊加復(fù)合分析進(jìn)行礦產(chǎn)遠(yuǎn)景區(qū)的預(yù)測。每種地學(xué)信息都是成礦預(yù)測的一個(gè)證據(jù)因子，而每個(gè)證據(jù)因子對成礦預(yù)測的貢獻(xiàn)是由這個(gè)因子的權(quán)重來確定的，進(jìn)而計(jì)算空間任意位置的礦產(chǎn)發(fā)育的概率值，以圈定不同級別的預(yù)測靶區(qū)［11］。其數(shù)學(xué)公式包括各地質(zhì)證據(jù)因素的找礦信息量前概估算和預(yù)測單元的后概計(jì)算兩部分。假設(shè)研究區(qū)T被劃分成面積相等的單元，其中D為含礦單元，對于任意一個(gè)證據(jù)因子，其權(quán)重定義為：

式中，W＋，W－分別為證據(jù)因子存在區(qū)和不存在區(qū)的權(quán)重值，對于原始數(shù)據(jù)缺失區(qū)域權(quán)重值為0；B為因子存在區(qū)的單元，為因子不存在的單元。如果W＋和W－不等于0，那么將永遠(yuǎn)取相異的符號。正權(quán)值與負(fù)權(quán)值之差C值的大小表示該證據(jù)因子的找礦指示性優(yōu)劣。

預(yù)測錳礦礦床產(chǎn)出的后驗(yàn)概率比（odds）的對數(shù)值等于先驗(yàn)概率比的對數(shù)值與各種地質(zhì)標(biāo)志的權(quán)系數(shù)之和。該法源于概率論中的貝葉斯關(guān)系式，在計(jì)算成礦后驗(yàn)概率時(shí)，需要用到有限個(gè)獨(dú)立隨機(jī)事件的概率乘法公式，因此，每一種控礦地質(zhì)因素相對于礦床產(chǎn)出這一概率事件來說，都必須是條件獨(dú)立的。對于n個(gè)證據(jù)因子，若它們關(guān)于礦點(diǎn)條件獨(dú)立，研究區(qū)任一單元k內(nèi)含有礦點(diǎn)的可能性，即后驗(yàn)幾率用O表示：

式中，b為基值：

證據(jù)權(quán)重法既考察地質(zhì)因素存在時(shí)的找礦權(quán)重，也考慮因素不存在的找礦權(quán)重，而且當(dāng)因素未知時(shí)也可以處理。

研究區(qū)任一單元k為含礦單元的后驗(yàn)概率為：

證據(jù)權(quán)重法已發(fā)展為一種比較成熟的成礦預(yù)測方法，國內(nèi)基于GIS的礦產(chǎn)預(yù)測軟件大多提供了該方法，如 MOPAS［12］，MRAS［13］和 GeoDAS［14］等。近年來，證據(jù)權(quán)重法不僅應(yīng)用于礦產(chǎn)資源評價(jià)［15－18］，還在地下水資源評價(jià)［19］、滑坡危險(xiǎn)性評價(jià)［20－23］、邊坡穩(wěn)定性評價(jià)［24］、污染危險(xiǎn)性評價(jià)［25］等領(lǐng)域得到了應(yīng)用。

3 優(yōu)質(zhì)錳礦成礦預(yù)測

為了滿足桂西南地區(qū)優(yōu)質(zhì)錳礦定量預(yù)測的需求，借助GIS空間分析技術(shù)，以錳礦多元地學(xué)空間數(shù)據(jù)庫［26］為基礎(chǔ)，引入成礦影響場和成礦距離場的 概念［27］來表達(dá)線狀地質(zhì)體對成礦的作用，通過錳礦床（點(diǎn)）地質(zhì)變量的資源量線性回歸來輔助進(jìn)行變量二值化，計(jì)算了證據(jù)因子的權(quán)重值，得到了桂西南地質(zhì)優(yōu)質(zhì)錳礦成礦有利度后驗(yàn)概率（圖1）。

3．1 網(wǎng)格單元的劃分

按照1／20萬區(qū)域礦產(chǎn)預(yù)測精度要求，將研究區(qū)劃分為2km×2km的規(guī)則單元（圖上面積1cm2），每個(gè)單元采用其行列號（ix，iy）來標(biāo)識。后續(xù)步驟的點(diǎn)、線、面地質(zhì)體的柵格化，以及含礦性因子、證據(jù)因子的組織，均采用統(tǒng)一的網(wǎng)格單元?jiǎng)澐?。然后，我們選取局部勘探程度高、含有較多找礦信息的區(qū)域內(nèi)的單元作為證據(jù)權(quán)重計(jì)算的訓(xùn)練區(qū)。

3．2 證據(jù)因子的建立

參與證據(jù)權(quán)重法計(jì)算的因子主要有2類：含礦性因子和證據(jù)因子。含礦性因子可直接通過將優(yōu)質(zhì)錳礦床（點(diǎn)）投射到網(wǎng)格單元即可得到，證據(jù)因子要由地質(zhì)變量二值化來得到。

圖1 基于證據(jù)權(quán)法的優(yōu)質(zhì)錳礦成礦預(yù)測流程圖Fig．1 Flow sheet of weight evedonce－based prediction of quality Mn ore

在建立證據(jù)因子之前，首先要構(gòu)建控制優(yōu)質(zhì)錳礦分布的地質(zhì)變量。對于面狀地質(zhì)體可以直接進(jìn)行柵格化，構(gòu)建的地質(zhì)變量表達(dá)了地質(zhì)類型的分布。對于線狀地質(zhì)體，本文摒棄了先進(jìn)行緩沖區(qū)分析、再將緩沖區(qū)柵格化的傳統(tǒng)方法，而是引入場論概念，用場模型（包括成礦影響場和成礦距離場）［27－28］來表達(dá)線狀地質(zhì)體對成礦的影響作用。圖2給出了研究區(qū)晚泥盆世同沉積斷裂的成礦影響場示例。采用場模型來表達(dá)地質(zhì)變量（還可對其進(jìn)行非線性變換）有效地解決了緩沖區(qū)方法將證據(jù)因子限定在線狀地質(zhì)體的一定距離范圍內(nèi)的弊病，使證據(jù)因子的分布更符合地質(zhì)實(shí)際情況。

地質(zhì)變量構(gòu)建完畢后，要對其進(jìn)行二值化得到證據(jù)因子。在證據(jù)因子建立過程中，我們在訓(xùn)練區(qū)內(nèi)對優(yōu)質(zhì)錳礦資源量與地質(zhì)變量進(jìn)行了線性回歸分析，如果線性關(guān)系顯著，則用該地質(zhì)變量構(gòu)建證據(jù)因子；如果錳礦資源量與地質(zhì)變量呈非線性關(guān)系，則對地質(zhì)變量進(jìn)行非線性變換，圖3給出了研究區(qū)晚泥盆世同沉積斷裂成礦影響場變換后的結(jié)果（變換公式為X′＝－｜X－25．75｜0．43），驗(yàn)證變換后的地質(zhì)變量與錳礦資源量線性關(guān)系顯著時(shí)，使用變換后的地質(zhì)變量構(gòu)建證據(jù)因子；對于那些與錳礦資源量沒有表現(xiàn)出相關(guān)關(guān)系的地質(zhì)變量，則不能用于構(gòu)建證據(jù)因子。通過線性回歸分析達(dá)到篩選證據(jù)因子的目的，避免了證據(jù)因子設(shè)置的盲目性和主觀性。

圖2 桂西南地區(qū)晚泥盆世同沉積斷裂成礦影響場Fig．2 Influence field of Late Devonian syn－sedimetary faults in Sowthwest Guang xi

圖3 桂西南地區(qū)晚泥盆世同沉積斷裂成礦影響場的變換結(jié)果Fig．3 Transform effect of influence field of Late Devonian syn－sedimentan fault in southwest Guangxi

對于與錳礦資源量呈線性關(guān)系的地質(zhì)變量（或變換后的地質(zhì)變量），綜合考察C值及證據(jù)因子中資源量所占的比例，來選取合適的分界值將地質(zhì)變量二分。如圖4所示，晚泥盆世同沉積斷裂成礦影響場（變換后）沿分界值－1．1分割進(jìn)行二值化，所取得的C值存在一個(gè)峰值，并且證據(jù)因子內(nèi)所占的錳礦資源量比例也比較高。在選取分界值時(shí)不僅僅考慮C值最大，還綜合考慮證據(jù)因子內(nèi)錳礦資源量所占的比例，可有效避免二值化時(shí)信息的損失。

設(shè)地質(zhì)變量X在訓(xùn)練區(qū)內(nèi)的取值范圍為［Xmin，Xmax］，二值化分界值為X0。若地質(zhì)變量與錳礦資源量為正相關(guān)，證據(jù)因子在單元K取值為：

反之，地質(zhì)變量與錳礦資源量呈負(fù)相關(guān)，則證據(jù)因子在單元K取值為：

圖4 桂西南晚泥盆世同沉積斷裂影響場二值化分界值的選取Fig．4 Boundary value selection of binaryzation of influence of Late Devonian syn－sedimentary fault in southwest Guangxi

由式（7），式（8）所得到的桂西南晚泥盆世同沉積斷裂影響場證據(jù)因子分布如圖5所示。

圖5 桂西南地區(qū)晚泥盆世同沉積斷裂影響場證據(jù)因子分布圖Fig．5 Map showing distribution of evedence factor of influence field of late Devonian syn－sedimentary fault in southwest Guangxi

表1 桂西南地區(qū)優(yōu)質(zhì)錳礦證據(jù)因子一覽表Table 1 Evedence factor schedule for quality Mn ore in southwest Guangxi

3．3 證據(jù)權(quán)重的計(jì)算

在上述優(yōu)質(zhì)錳礦形成條件分析的基礎(chǔ)上，選取優(yōu)質(zhì)錳礦的地層與巖性、同沉積斷裂、沉積盆地、沉積相、巖性組合、地層等厚線、切割度、航磁異常等證據(jù)因子共15個(gè)，并在訓(xùn)練區(qū)范圍內(nèi)計(jì)算了證據(jù)因子的正負(fù)權(quán)重值，如表1所示。

3．4 成礦有利度的計(jì)算

圖6 桂西南地區(qū)優(yōu)質(zhì)錳礦成礦預(yù)測圖Fig．6 Map showing prediction of quality Mn ore in southwest Guangxi

利用貝葉斯疊加模型計(jì)算了研究區(qū)每個(gè)單元的成礦有利度后驗(yàn)概率，從預(yù)測結(jié)果來看，優(yōu)質(zhì)錳礦主要位于龍邦—下雷—那利—足榮—東平一帶，一直延伸到田東縣境內(nèi)，是晚泥盆世、早石炭世大塘期、早三疊世晚期的主要成錳沉積盆地分布區(qū)域（圖6）。

在成礦有利度分布圖上，劃分了1個(gè)Ⅰ類遠(yuǎn)景區(qū)：Ⅰ－1龍邦—下雷—土湖成礦遠(yuǎn)景區(qū)（成錳時(shí)代為晚泥盆世）；2個(gè)II類遠(yuǎn)景區(qū)：II－1把荷—那利成礦遠(yuǎn)景區(qū)（成錳時(shí)代為晚泥盆世和早石炭世大塘期），II－2足榮—龍懷—東平成礦遠(yuǎn)景區(qū)（成錳時(shí)代為早三疊世）。

4 結(jié)論

在對優(yōu)質(zhì)錳礦形成條件分析的基礎(chǔ)上，利用GIS空間分析方法構(gòu)造了優(yōu)質(zhì)錳礦的地層與巖性、同沉積斷裂、沉積盆地、沉積相、巖性組合、地層等厚線、切割度、航磁異常等證據(jù)因子共15個(gè)；采用證據(jù)權(quán)重法對桂西南地區(qū)優(yōu)質(zhì)錳礦成礦有利度進(jìn)行了預(yù)測，并劃分了1個(gè)Ⅰ類和2個(gè)Ⅱ類優(yōu)質(zhì)錳礦成礦遠(yuǎn)景區(qū)。取得的主要認(rèn)識如下：

（1）引入場模型來表達(dá)線狀地質(zhì)體（同沉積斷裂、盆地邊界）對成錳的影響作用，解決了傳統(tǒng)緩沖區(qū)方法將證據(jù)因子限定在線狀地質(zhì)體的一定距離范圍內(nèi)的弊病，使證據(jù)因子的分布更符合地 質(zhì)實(shí)際情況。

（2）在對地質(zhì)變量二值化前，對優(yōu)質(zhì)錳礦資源量與地質(zhì)變量進(jìn)行線性回歸分析，根據(jù)線性關(guān)系是否顯著來篩選證據(jù)因子，避免了證據(jù)因子設(shè)置的盲目性和主觀性。

（3）在對地質(zhì)變量二值化選取分界值時(shí)不僅僅考慮C值最大，還綜合考慮證據(jù)因子內(nèi)錳礦資源量所占的比例，避免了二值化時(shí)信息的損失。

致謝：教育部GIS軟件及其應(yīng)用工程中心與中南大學(xué)共建“MAPGIS實(shí)驗(yàn)室”為本項(xiàng)研究工作提供了MAPGIS軟件，在此表示感謝。

［1］駱華寶．我國優(yōu)質(zhì)錳礦的勘查方向［J］．地質(zhì)與勘探，2002，38（4）：8－11．

［2］何知禮．從地洼成礦作用看中國某些優(yōu)質(zhì)錳礦的成因與遠(yuǎn)景［J］．大地構(gòu)造與成礦學(xué)，1990，14（3）：203－210．

［3］He Zhili．Study on genesis and perspective of some high quality manganese ore deposits in China［J］．Geotectonica et Metallogenia，1995，19（2－3）：28－38．

［4］茹廷鏘，韋靈敦，樹皋．廣西錳礦地質(zhì)［M］．北京：地質(zhì)出版社，1992．

［5］侯宗林，薛友智，黃金水，等．揚(yáng)子地臺周邊錳礦［M］．北京：冶金工業(yè)出版社，1997．

［6］曾友寅．廣西下雷晚泥盆世錳礦床沉積學(xué)研究［J］．沉積學(xué)報(bào)，1992，9（1）：73－79．

［7］曾友寅．廣西下雷錳礦床沉積學(xué)研究［C］∥范德廉．錳礦床地質(zhì)地球化學(xué)研究．北京：氣象出版社，1994：69－75．

［8］劉騰飛．廣東東平表生富集型錳礦床地質(zhì)特征及成礦條件初步研究［J］．地質(zhì)找礦論叢，1996，11（4）：42－55．

［9］毛先成，曾文波，周尚國，等．桂西—滇東南地區(qū)數(shù)字地形特征與次生富集錳礦空間分布關(guān)系研究［J］．地質(zhì)與勘探，2009，45（3）：292－298．

［10］陳祥軍，周眉成，曲力群．微生物在錳的氧化富集過程中的作用——以廣西湖潤錳礦為例［J］．地質(zhì)與勘探，2003，39（1）：23－26．

［11］肖克炎，張曉華，陳鄭輝，等．成礦預(yù)測中證據(jù)權(quán)重法與信息量法及其比較［J］．物探化探計(jì)算技術(shù)，1999，21（3）：223－226．

［12］陳永清，夏慶霖，黃靜寧，等．“證據(jù)權(quán)”法在西南“三江”南段礦產(chǎn)資源評價(jià)中的應(yīng)用［J］．中國地質(zhì)，2007，34（1）：132－141．

［13］肖克炎，王勇毅，薛群威，等．中國銅礦數(shù)字礦床模型評價(jià)系統(tǒng)的開發(fā)［J］．礦床地質(zhì)，2003，22（4）：425－429．

［14］成秋明，陳志軍，Khaled A．模糊證據(jù)權(quán)方法在鎮(zhèn)沅（老王寨）地區(qū)金礦資源評價(jià)中的應(yīng)用［J］．地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào)，2007，32（2）：175－184．

［15］Wang Z ，Cheng Q．GIS－Based（W＋－W－）Weight of Evidence Model and Its Application to Gold Resources Assess－ment in Abitibi，Canada［J］．Journal of China University of Geosciences，2006，17（1）：71－78．

［16］徐善法，陳建平，葉繼華．證據(jù)權(quán)法在三江北段銅金礦床成礦預(yù)測中的應(yīng)用研究［J］．地質(zhì)與勘探，2006，42（2）：54－59．

［17］鄧勇，邱瑞山，羅鑫．基于證據(jù)權(quán)重法的成礦預(yù)測——以廣東省鎢錫礦的成礦預(yù)測為例［J］．地質(zhì)通報(bào)，2007，26（9）：1228－1234．

［18］Ford A，Blenkinsop T G．Combining fractal analysis of mineral deposit clustering with weights of evidence to evaluate patterns of mineralization：Application to copper deposits of the Mount Isa Inlier，NW Queensland，Australia［J］．Ore Geology Reviews，2008，33（3－4）：435－450．

［19］Corsini A，Cervi F，Ronchetti F．Weight of evidence and artificial neural networks for potential groundwater spring mapping：an application to the Mt．Modino area（Northern Apennines，Italy）［J］．Geomorphology，2009，111（1－2）：79－87．

［20］Neuh?user B，Terhorst B．Landslide susceptibility assessment using"weights－of－evidence"applied to a study area at the Jurassic escarpment （SW－Germany）［J］．Geomorphology，2007，86（1－2）：12－24．

［21］Dahal R K，Hasegawa S，Nonomura A，et al．Predictive modelling of rainfall－induced landslide hazard in the Lesser Himalaya of Nepal based on weights－of－evidence［J］．Geomorphology，2008，102（3－4）：496－510．

［22］Song R－H，Hiromu D，Kazutoki A，et al．Modeling the potential distribution of shallow－seated landslides using the weights of evidence method and a logistic regression model：a case study of the Sabae Area，Japan［J］．International Journal of Sediment Research，2008，23（2）：106－118．

［23］Regmi N R，Giardino J R，Vitek J D．Modeling susceptibility to landslides using the weight of evidence approach：Western Colorado，USA［J］．Geomorphology，2010，115（1－2）：172－187．

［24］Zahiri H，Palamara D R，F(xiàn)lentje P，et al．A GIS－based Weights－of－Evidence model for mapping cliff instabilities associated with mine subsidence［J］．Environmental Geology，2006，51（3）：377－386．

［25］Morales－Caselles C，Riba I，DelValls T á．A weight of evidence approach for quality assessment of sediments impacted by an oil spill：The role of a set of biomarkers as a line of evidence［J］．Marine Environmental Research，2009，67（1）：31－37．

［26］張寶一，毛先成，周尚國，等．礦產(chǎn)資源預(yù)測評價(jià)數(shù)據(jù)庫的設(shè)計(jì)與實(shí)現(xiàn)——以桂西—滇東南錳礦為例［J］．地質(zhì)與勘探，2009，45（6）：697－703．

［27］毛先成，周尚國，張寶一，等．基于場模型的成礦信息提取方法研究——以桂西—滇東南錳礦為例［J］．地質(zhì)與勘探，2009，45（6）：704－715．

［28］毛先成，汪凡云．基于場論的動態(tài)緩沖區(qū)生成方法與應(yīng)用［J］．測繪科學(xué)，2009，34（6）：11－14．