雷天賜，羅士新，王　磊

（中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心，武漢430205）

基于多元信息的證據(jù)權(quán)法在南嶺鎢錫多金屬成礦預測中的應用

雷天賜，羅士新，王磊

（中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心，武漢430205）

以南嶺成礦帶1∶50萬地質(zhì)數(shù)據(jù)庫、物化探數(shù)據(jù)為信息源，通過數(shù)據(jù)處理、分析與推斷，編制相關(guān)成果圖件。以GIS為平臺、地質(zhì)異常理論為指導，基于多元信息成礦耦合分析與異常提取，優(yōu)選找礦要素有利區(qū)間并構(gòu)建了13個有效證據(jù)層，基于Morpas、Arc-wofe證據(jù)權(quán)模型開展了成礦定量預測與評價。以0.024、0.078和0.098為后驗概率臨界值劃分了三級成礦遠景區(qū)，共圈定一級遠景區(qū)41個、二級遠景區(qū)70個和三級遠景區(qū)85個，為該區(qū)找礦工作部署及資源規(guī)劃提供了參考。

南嶺成礦帶；優(yōu)益度；單位礦產(chǎn)當量；證據(jù)權(quán)法；成礦預測

1　前言

南嶺成礦帶屬于我國20個重點成礦帶之一，跨湘、粵、桂、贛四?。▍^(qū)），地理坐標范圍：107°～116°E、24°～27°N，面積約30萬km2。構(gòu)造上位于揚子板塊、華夏板塊交匯處，區(qū)內(nèi)巖漿活動強烈、成礦條件十分優(yōu)越，是我國有色、稀有及貴金屬的重要生產(chǎn)基地，也是世界鎢和原生錫礦床分布最密集、最有潛力的地區(qū)之一，擁有世界上主要鎢、錫礦床類型。

南嶺地區(qū)鎢錫大規(guī)模成礦作用發(fā)生于中晚侏羅紀（165～150 Ma），古太平洋板塊向大陸俯沖時形成一系列NE向伸展帶與EW向古斷裂的交匯部位是巖漿活動和成礦作用中心區(qū)［1］，鎢錫多金屬礦床（點）具有成帶、成群分布特點，主要沿古板塊結(jié)合帶、大型隆起區(qū)與坳陷區(qū)結(jié)合部、深大斷裂（帶）等3個部位產(chǎn)出原生錫礦，包括云英巖型、變花崗巖型、矽卡巖型、石英脈型、破碎帶蝕變巖型和斑巖型等6個主要類型，與H型和鋁質(zhì)A型花崗巖關(guān)系密切［2-3］。

證據(jù)權(quán)法是一種離散的多元統(tǒng)計方法，運用的是相似類比理論。其最早應用是在上世紀60年代的醫(yī)學領域，80年代末由Agterberg等［4］通過改進后首次引入地學領域。隨后，作為一種新的數(shù)據(jù)模型方法，廣泛應用于國內(nèi)外礦產(chǎn)資源預測中并取得巨大成功［5-6］。目前，國內(nèi)應用主要是基于Morpas、MRAS軟件及Arc-Wofe模塊。

本文以南嶺成礦帶為研究區(qū)、Morpas為平臺、鎢錫多金屬礦為預測對象，通過地、物、化、遙數(shù)據(jù)處理與分析、異常信息提取與成礦耦合關(guān)系研究，參照單位礦產(chǎn)當量指標優(yōu)選找礦要素有利區(qū)間構(gòu)建有效證據(jù)層，通過證據(jù)權(quán)重模型自動圈定并劃分遠景區(qū)。

2　證據(jù)權(quán)法及其理論基礎

證據(jù)權(quán)法是一種基于Bayes的統(tǒng)計方法，它通過對一些與礦產(chǎn)形成相關(guān)的地學信息的疊加復合分析來進行礦產(chǎn)遠景區(qū)的預測，其中的每一種地學信息都被視為成礦遠景區(qū)預測的一個證據(jù)因子，而每一個證據(jù)因子對成礦預測的貢獻都是由這個因子的權(quán)重值來確定的。

3　多源信息提取、評價及證據(jù)權(quán)層優(yōu)取

3.1構(gòu)建成礦信息評價因子

對于成礦要素而言，評價其是否有利于成礦不僅看其賦礦數(shù)量，而且還要充分考慮其在預測中的有效性。因此，構(gòu)建單位礦產(chǎn)當量（KN）評價因子作為統(tǒng)一衡量標準具有一定的現(xiàn)實意義，其數(shù)學表達式為［7］：

（l）式中，S為賦礦因子出露面積，單位km2；N為由（2）式計算得到的礦產(chǎn)當量，單位個。礦產(chǎn)當量計算表達式如下：

（2）式中，N0、N1、N2、N3、N4、N5分別為超大、大、中、小型礦床、礦點和礦化點的個數(shù)；K0、K1、K2、K3、K5是超大、大、中、小型礦床和礦化點相應的權(quán)系數(shù)。在南嶺成礦帶多元信息因子找礦有利度分析中，W_Sn_Mo_Bi礦產(chǎn)當量（N）計算取K0=625、K1=125、K2=25、K3=5、K5=0.2，總礦產(chǎn)當量為8636.4個。

3.2地質(zhì)信息提取與優(yōu)選

3.2.1控礦構(gòu)造

南嶺地區(qū)斷裂多形成于前加里東期、加里東期、印支期和燕山期，構(gòu)造特征主要為壓扭性和張扭性；深大斷裂控制了礦床的形成與分布，具有多期活動的特點，以NE向為主，伴有近SN向、近EW向和NW向［8］；控礦模式空間上表現(xiàn)為橫成排、縱成列，交叉部位最優(yōu)越［9］。NE向斷裂主要為燕山期以前各時期及不同期次的疊加復合。前燕山期形成的為基礎斷裂，對后期的形成起著明顯的控制作用；燕山期的是重要的導巖導礦斷裂，其活動對花崗巖的侵位、鎢錫及多金屬礦的形成有著明顯的控制作用。近SN向深大斷裂大多形成于印支期和燕山期，印支期的主要表現(xiàn)為壓性或壓扭性，與蓋層中的波狀弧形褶皺走向一致；燕山期的主要表現(xiàn)為波狀弧形，以壓扭-張扭活動為主。

本次分析斷裂主要來源于重磁推斷、遙感解譯及1∶50萬地質(zhì)數(shù)據(jù)庫。為全面解析其控礦特征，依據(jù)線性構(gòu)造地質(zhì)意義或巖石力學意義引入4個評價指標:線性構(gòu)造buffer區(qū)、優(yōu)益度、等密度和方位異常度［10］；再參照它們的單位礦產(chǎn)當量值［11］，優(yōu)選了各自的有利區(qū)間段分別構(gòu)建了有效證據(jù)層。

線性構(gòu)造buffer分析

斷裂構(gòu)造對礦床（點）形成起著重要的作用，尤其是深大斷裂，而且其對控礦的影響范圍不一。結(jié)合本區(qū)構(gòu)造控礦特征，選擇深大斷裂和次級斷裂做半徑（r）1～5 km的緩沖區(qū)分析，計算各緩沖區(qū)的單位礦產(chǎn)當量值。根據(jù)各緩沖區(qū)r-KN統(tǒng)計圖（圖1a）看，KN與r呈反比、KN值隨r的增大而降低；當3 km＜r＜4 km時曲線變化趨勢急劇放緩并出現(xiàn)明顯拐點，說明該區(qū)3 km范圍是斷裂對礦產(chǎn)分布影響的拐點和最佳極值。

線性構(gòu)造優(yōu)益度分析

一種有利儲礦的構(gòu)造特征函數(shù)，理論依據(jù)是多組斷裂交匯處或斷裂密集部位為內(nèi)生礦產(chǎn)有利容礦構(gòu)造部位。它以線性構(gòu)造兩兩間的夾角和線性構(gòu)造的方位的控礦程度加權(quán)的構(gòu)造密度的度量，通過斷層定量分析來對有利成礦區(qū)進行評價。其表達式如下：

通過優(yōu)益度等值線圖ε-KN統(tǒng)計分析（圖1b）看，按單位礦產(chǎn)當量值選擇0～11的4個區(qū)間為優(yōu)益度證據(jù)層。

線性構(gòu)造方位異常度分析

方位異常度是一種排除區(qū)域優(yōu)勢方位構(gòu)造影響的程度指數(shù)，評價局部性構(gòu)造的指標。一般方位異常度指數(shù)越高表示其受優(yōu)勢方位構(gòu)造的影響越小，反映局部構(gòu)造特征越強。通過單元網(wǎng)格構(gòu)建、以22.5°為間距在-90°～90°進行方位劃定，編制了方位異常圖。

基于線性構(gòu)造方位異常圖的方位異常度-單位礦產(chǎn)當量統(tǒng)計圖（圖1c）看，優(yōu)選了NS、NNE和NE向3個方位區(qū)為方位異常度有利成礦區(qū)。

線性構(gòu)造等密度分析

指統(tǒng)計網(wǎng)格單元格中各線狀構(gòu)造的長度之和，一定程度上揭示了線性構(gòu)造空間分布結(jié)構(gòu)特征及深部構(gòu)造信息和找礦線索［12］。其表達式為：

基于線性構(gòu)造等密度圖的L-KN統(tǒng)計分析結(jié)果（如圖1d），優(yōu)選了0～19間的4個區(qū)段為構(gòu)造等密度有利成礦區(qū)。

3.2.2賦礦地層

南嶺地區(qū)鎢錫多金屬礦的主要類型有石英脈型、花崗巖型（蝕變花崗巖型、云英巖型、花崗巖巖脈型、斑巖型）和矽卡巖型，表現(xiàn)為“多位一體”的復合礦床模式。其中，石英脈型分布最廣、花崗巖型和矽卡巖型規(guī)模最大。

區(qū)內(nèi)矽卡巖性礦床成礦方式為燕山期花崗巖類侵入于上寒武統(tǒng)、奧陶系、泥盆系、石炭系及二疊系等地層中，與地層中的灰?guī)r、鈣質(zhì)砂巖等發(fā)生接觸交代形成矽卡巖型鎢錫多金屬礦床。礦體呈層狀、似層狀、透鏡狀、囊包狀、不規(guī)則狀；成礦作用明顯具多階段性，并呈順向分帶。主要礦物組合：錫石、白鎢礦、方鉛礦、閃鋅礦、黃銅礦、黃鐵礦等。當多種成因類型復合時，則形成復控型礦床，往往形成大-超大型礦床。

利用GIS空間分析與統(tǒng)計功能，南嶺地區(qū)優(yōu)選出了3個W-Sn-Mo-Bi成礦元素的賦礦地層組：Nh-O、D-P和J-Qb，其中Nh-O、D-P地層組包含的礦床（點）數(shù)約占總礦床（點）數(shù)的90%。根據(jù)賦礦地層巖性單位礦產(chǎn)當量值，最終提取了Nh-O、D-P地層組構(gòu)建有效證據(jù)層。

3.2.3巖漿巖

研究表明，南嶺地區(qū)成礦巖體主要為燕山期侵入花崗巖體［13］。成礦花崗巖體多呈巖株、巖枝、巖筒、巖脈、巖瘤、巖盤或巖床產(chǎn)出，與圍巖間接觸界線清楚，其礦化富集部位常常在巖體的頂部或內(nèi)外接觸帶；部分成礦花崗巖體具有隱爆特征，因而常有角礫巖伴生，為超淺成或次火山相產(chǎn)物。

圖1　南嶺地區(qū)線性構(gòu)造各參數(shù)單位礦產(chǎn)當量統(tǒng)計圖Fig.1 Histogram of unit ore equivalence of all parameters on linear structure in Nanling region

區(qū)內(nèi)的鎢錫多金屬礦床在成因上與燕山期的復式花崗巖，尤其是燕山中晚期酸性、中酸性巖體密切相關(guān)。這些含礦巖體或產(chǎn)在大型巖基的邊緣部位、或是在隱伏巖基的頂部，一般都呈中、小型巖株。巖漿巖控礦特征：巖性上，含礦巖體的主要巖石類型為黑云母花崗巖、白云母花崗巖和二長花崗巖等；在成礦時間上，成礦年齡與成巖年齡比較接近；在空間上，礦床圍繞巖體分布，且自巖體中心向外礦化具有明顯的分帶性，依次出現(xiàn)Mo、W、Bi、Sn、Cu、Pb、Zn等成礦特征。

根據(jù)重磁推斷酸性巖體、巖漿巖帶成果及1∶50萬地質(zhì)圖，編制了區(qū)內(nèi)酸性巖體、巖漿巖帶分布圖，選擇成鎢、成錫花崗巖體組合［14］及邊緣外擴3km區(qū)構(gòu)建有效證據(jù)層。

3.3重磁信息提取與優(yōu)選

3.3.1重磁數(shù)據(jù)處理與信息提取

以南嶺地區(qū)1∶20萬全國重磁掃面數(shù)據(jù)為對象，RGIS、Mapgis為平臺，通過濾波、單元網(wǎng)格化和克里格插值等處理編制了布格重力異常、航磁ΔT等值線平面圖等地球物理系列圖件［15］。本區(qū)地質(zhì)背景與地球物理場聯(lián)系密切，礦化集中區(qū)主要分布在重力梯度帶偏負值一側(cè)，尤其是重力梯度帶的彎曲處、舌狀重力低或短軸狀重力異常邊緣，如南丹礦化集中區(qū)、都龐嶺礦化集中區(qū)；航磁方面，礦化集中區(qū)主要位于負磁異常內(nèi)及彎曲處，正負磁場的交變帶或短軸狀正磁力異常的周緣。

本次重磁找礦有利度分析主要圍繞布格剩余重力異常和航磁ΔT異常開展。首先，通過GIS空間分析與統(tǒng)計功能，求出重磁各區(qū)間面積及對應的礦產(chǎn)當量值；然后，計算各區(qū)間的單位礦產(chǎn)當量，如圖2；再參考單位礦產(chǎn)當量指標優(yōu)選重磁有利區(qū)間段構(gòu)建有效證據(jù)層。

剩余重力異常單位礦產(chǎn)當量分布在起始值為（-14～4）×10-5/s2間，明顯呈3個區(qū)間段；高值分布在梯度帶附近偏負重力場一側(cè)，該區(qū)間可進一步優(yōu)選為2個區(qū)間：（-14～-12）×10-5/s2和（-10～-7）×10-5/s2。航磁ΔT異常單位礦產(chǎn)當量分布在-200～320nT/Km范圍間，呈不規(guī)則齒狀，提取的3個有效證據(jù)層區(qū)間分別為：-160nT/km、20～140nT/km和220～320nT/km。

3.3.2重磁異常地質(zhì)推斷

重磁主要推斷了區(qū)內(nèi)斷層、巖體和巖漿巖帶（圖3）。依據(jù)重磁異常特征（如線性梯度帶、異常特征的分界線、異常的錯動、等值線的規(guī)則性扭曲、異常寬度突變和串珠狀異常等），全區(qū)推斷深大斷裂12條、中層斷裂40條和淺層斷裂若干條。這些斷裂不僅控制了本區(qū)地質(zhì)構(gòu)造基本格架，與成巖成礦有著密不可分的關(guān)系，而且大多數(shù)深大斷裂（帶）還構(gòu)成了大地構(gòu)造單元的分界線。

圖2　重磁異常區(qū)單位礦產(chǎn)當量統(tǒng)計圖（左：剩余重力異常，右：航磁△T異常）Fig.2 Histogram of unit ore equivalence of gravity and magnetic anomaly zone（Left:residual gravityanomaly；Right:aeromagnetic anomaly）

利用重磁推斷巖體118個，包括已知巖體49個、半隱伏巖體28個和隱伏巖體41個，其中部分巖體表現(xiàn)出良好的重磁同源異常和成礦相關(guān)性。在剩余布格重力異常圖上，出露花崗巖體呈低重異常，密集封閉的環(huán)狀等值線反映了產(chǎn)狀較陡的接觸帶；布格重力異常垂向二階導中，負異常較好的反映了花崗巖體，零值線為巖體邊界。此外，依據(jù)上延高度不同的重力異常特征，在南嶺成礦帶初步劃分13條大型巖漿巖帶。

航磁ΔT局部磁異常大體上反映了巖漿巖活動強弱狀況。以茶陵-彬州-恭城為界，西側(cè)巖體較少，局部磁異常也較少；東側(cè)出露巖體眾多，局部磁異常也顯著增多。

3.4化探信息提取與優(yōu)選

3.4.1地球化學數(shù)據(jù)處理與元素選擇

化探分析數(shù)據(jù)來源于全國1∶20萬水系沉積物數(shù)據(jù)庫，共包括39種元素及化合物。本次成礦預測礦種元素為W、Sn，根據(jù)華南地區(qū)燕山期重熔型花崗巖類有關(guān)的成礦系列的特征［16］，選擇了W、Sn、Cu、Pb、Zn、Mo、Bi、Ag、As、F、B、Hg、Sb等13種與成礦相關(guān)的元素做R型分層聚類分析，以歐氏距離作為衡量變量間相似程度的統(tǒng)計量，聚類方法采用離差平方和法，聚類譜系圖見圖4。

圖3　南嶺地區(qū)重磁推斷主要斷裂、酸性巖體及巖漿帶Fig.3 Map of interpretation of main faults，magmatic rocks and magmatic belt bygravity and magneticin Nanling region

相似系數(shù)取11.25時，W、Sn、Bi、Mo為一組，為高溫礦物。其中W、Sn為區(qū)內(nèi)主成礦元素，Bi、Mo為副成礦元素。結(jié)合成礦元素相關(guān)性分析，選擇了W、Sn、Bi、Mo 4種元素開展地球化學分析與異常提取。

圖4　南嶺地區(qū)13種元素聚類分析龍骨圖Fig.4 Diagram of cluster analysis of 13 elements in Nanling region

3.4.2地球化學異常信息提取

首先，對數(shù)據(jù)進行中值低通濾波處理，填補原始數(shù)據(jù)缺失和消除孤立的噪聲點，保證數(shù)據(jù)完整性；其次，通過數(shù)據(jù)處理與統(tǒng)計分析，采用迭代法、累頻率法分別確定各元素異常下限，并求取元素相應的均值、標準方差、中位數(shù)及變異系數(shù)等（表1）；再與礦產(chǎn)疊加分析，優(yōu)選異常下限求法：迭代法以均值+2倍方差確定異常下限、累頻法以累積頻率的85%確定異常下限值；最后，迭代法按異常下限、2倍異常下限和4倍異常下限劃分異常區(qū)內(nèi)帶、中帶和外帶，累頻法以累積頻率的85%、95%和98.5%為閥值，劃分異常區(qū)內(nèi)帶、中帶和外帶，編制地球化學異常圖（圖5）并作為有效證據(jù)層。

表1　南嶺地區(qū)成礦元素異常參數(shù)Table 1 Anomaly parameter of metallogenic elements in Nanling region

4　基于Morpas的證據(jù)權(quán)法模型的建立

根據(jù)預測比例尺及精度要求，工作區(qū)選擇4 km×4 km大小進行網(wǎng)格單元創(chuàng)建，盡量保證一個網(wǎng)格單元內(nèi)只有一個礦床（點），全區(qū)共劃分網(wǎng)格單元19465個。基于Morpas證據(jù)權(quán)重法預測模型與流程［17］（圖6），分別導入南嶺地區(qū)地質(zhì)、物化探證據(jù)層專題圖件，構(gòu)建了證據(jù)層并分別計算出各證據(jù)層證據(jù)權(quán)重值和成礦的相關(guān)系數(shù)（表2），對研究區(qū)各網(wǎng)格單元進行成礦后驗概率計算。

圖5　南嶺地區(qū)W_Sn_Mo_Bi地球化學組合異常圖Fig.5 Anomaly map of combined elements of W_Sn_Mo_Bi in Nanling region

從找礦要素有利區(qū)間構(gòu)建的有效證據(jù)層權(quán)重值（W+和W-）及相關(guān)系數(shù)（C=W+-W-）判斷，地球化學手段系列證據(jù)層影響最大，是本次找礦預測最重要變量，對找礦具有重要的指示意義；W、Sn、Mo、Bi組合異常高權(quán)重值反映了該區(qū)典型的成礦系列特征。燕山期酸性巖體證據(jù)層相關(guān)系數(shù)大于2，反映了其在系列控礦要素中成礦優(yōu)越性。賦礦地層和控礦構(gòu)造系列參數(shù)證據(jù)層影響相對較小，分析原因有二：一、由于成礦后期地質(zhì)構(gòu)造過于復雜或構(gòu)造活動太強，導致賦礦層位和礦體遭受不同程度破壞和剝蝕；二、地層以組合（或建造）形式構(gòu)建的證據(jù)層，線性構(gòu)造分解參數(shù)過多，降低了主要賦礦層和控礦構(gòu)造的影響。而重磁主要用來推斷和驗證構(gòu)造和巖體等地質(zhì)體，是一種間接的找礦信息，直接成礦指示作用不明顯。

條件獨立性檢驗對證據(jù)權(quán)法來說至關(guān)重要，如果地質(zhì)找礦要素不滿足條件獨立性就會引起后驗概率估計上的偏差［18］。

因此，計算后驗概率前，基于Arc-wofe模塊開展了條件獨立性檢驗，在顯著性水平為0.02時13個證據(jù)層彼此滿足條件獨立。

圖6　證據(jù)權(quán)重法模型流程圖Fig.6 Flow chart of weights of evidence model

5　預測結(jié)果及評價

成礦遠景區(qū)圈定與分級是基于后驗概率臨界值定的，后驗概率臨界值一般采用做含礦單元成礦概率值-累積面積對數(shù)曲線圖求取。根據(jù)經(jīng)驗，一般視含礦單元成礦概率值-累積面積對數(shù)曲線變形點為成礦后驗概率拐點值，即為成礦概率臨界值［19］。借助Arcview軟件，通過Arc-Wofe模塊自動繪出后驗概率值-累積面積關(guān)系變化曲線圖（圖7）；參考Boleneus D E等［20］觀點，總結(jié)建立了如下4點確定成礦概率臨界值標準：①曲線變化圖中曲線自然變形點；②盡可能包含更多的已發(fā)現(xiàn)的成礦元素礦床；③成礦預測區(qū)應當含有較高后驗概率值；④成礦預測區(qū)應當存在具有最高成礦相關(guān)度（C值）的專題證據(jù)層。

表2　南嶺地區(qū)鎢錫多金屬礦預測證據(jù)權(quán)重值Table 2 Weight of evidence factor for W_Sn polymetallic deposits in Nanlingregion

參照上述標準，結(jié)合成礦帶地質(zhì)背景及成礦規(guī)律，本區(qū)確定了3個臨界值：P=0.024、P=0.078和P=0.098，按臨界值大小劃分4個區(qū)域，劃分標準為：P＜0.024為非成礦區(qū)域、0.024＜P＜0.078為可能成礦區(qū)三級遠景區(qū)、0.078＜P＜0.098為可能成礦區(qū)二級遠景區(qū)、P＞0.098為有利成礦區(qū)一級遠景區(qū)，見表3。按求取的后驗概率臨界值，圈定一級遠景區(qū)41個、二級遠景區(qū)70個和三級遠景區(qū)85個。最終，按成礦帶構(gòu)造特征、巖漿活動、地層巖性和成礦特征，將196個遠景區(qū)整合為10個潛在成礦區(qū)域（Ⅰ～X），并編制了南嶺成礦帶鎢錫多金屬礦成礦遠景區(qū)預測圖（圖8）。

圖7　后驗概率值與累積面積關(guān)系變化曲線圖Fig.7 Curve diagram of relation between Post-probability values and cumulative area

表3　后驗概率臨界值及遠景區(qū)分級Table 3 Critical value of post-probability and classification of prospect

從成礦帶面積分配來看，有利成礦區(qū)占全區(qū)面積的3.75%、可能成礦區(qū)（0.024＜P＜0.098）占全區(qū)面積的0.06%，而非成礦區(qū)域占全區(qū)面積的96.19%。

從礦床數(shù)量統(tǒng)計數(shù)據(jù)來看，共有571個W、Sn、Mo、Bi礦床（點）落入可能成礦預測區(qū)，占總礦床（點）數(shù)的86%。其中：一級遠景區(qū)內(nèi)分布168個礦床（點）、二級遠景區(qū)內(nèi)分布277個礦床（點）、三級遠景區(qū)內(nèi)分布126個礦床（點）。

因此，根據(jù)礦床（點）空間分布特征、區(qū)域面積分配與成礦概率關(guān)系判斷表明，基于證據(jù)權(quán)法的鎢錫多金屬成礦預測結(jié)果是可行的。

6　結(jié)論與問題

圖8　南嶺地區(qū)鎢錫成礦預測圖Fig.8 Prognostic map of W_Sn polymetallic deposits in Nanling region

（1）利用地質(zhì)、物化探資料，通過數(shù)據(jù)處理、分析及推斷，提取了南嶺地區(qū)鎢錫多金屬礦的重要控礦要素。基于成礦耦合分析與異常提取，依據(jù)單位礦產(chǎn)當量指標優(yōu)選并構(gòu)建了找礦要素有效證據(jù)層。采用證據(jù)權(quán)法開展了南嶺地區(qū)鎢錫多金屬礦定量預測與評價，圈定并分級找礦遠景區(qū)196個，經(jīng)初步檢驗預測結(jié)果基本可行。

（2）基于GIS平臺的多元信息集成與證據(jù)層權(quán)重預測，不僅數(shù)據(jù)集成度高、綜合性強，而且主次突出、成果表達上可視性強；與傳統(tǒng)人工預測相比具有更客觀、更準確、更高效優(yōu)勢，是今后開展礦產(chǎn)資源快速預測和評價的重要方向。

（3）后驗概率臨界值還是一個基于數(shù)理統(tǒng)計法得出的經(jīng)驗值，其與成礦內(nèi)在關(guān)聯(lián)性不明，還有待進一步的機理研究。

［1］毛景文，謝桂青，郭春麗，陳毓川.南嶺地區(qū)大規(guī)模鎢錫多金屬成礦作用:成礦時限及地球動力學背景［J］.巖石學報，2007，23（10）:2329-2338.

［2］付建明，程順波，盧友月，馬麗艷.南嶺地區(qū)鎢錫多金屬礦成礦規(guī)律及找礦方向 ［J］.地球科學進展，2012，27（增刊）: 162-165.

［3］魏道芳，潘仲芳，金光富.南嶺錫礦調(diào)查評價主要進展及找礦前景分析［J］.華南地質(zhì)與礦產(chǎn)，2005，（2）:2-11.

［4］Bonham-Carter G F，Agterberg F P and Wright D F.Weights of evidence modeling:a new approach to mapping mineral potential［A］.Statistical Applications in the Earth Sciences［C］. Geological Survey of Canada，1989:89-9，171-183.

［5］陳建平，王功文，侯昌波，唐菊興.基于GIS技術(shù)的西南三江北段礦產(chǎn)資源定量預測與評價［J］.礦床地質(zhì)，2005，24（1）: 15-24.

［6］丁清峰，孫豐月.基于專家證據(jù)權(quán)重法的成礦遠景區(qū)劃與評價以東昆侖地區(qū)金礦為例［J］.地質(zhì)科技情報，2006，25（1）:41-46.

［7］池順都，周順平，吳新林.GIS支持下的地質(zhì)異常分析及金屬礦產(chǎn)經(jīng)驗預測 ［J］.地球科學-中國地質(zhì)大學學報，1997，22（1）:99-103.

［8］付建明，徐德明，楊曉君，馬麗艷，等.南嶺錫礦［M］.武漢:中國地質(zhì)大學出版社，2011:11-14.

［9］陸慶嵐.南嶺及鄰區(qū)深部構(gòu)造與鎢礦分布［J］.河北地質(zhì)學院學報，1993，16（4）:333-344.

［10］趙不億，秦小光.遙感構(gòu)造的定量分析方法［J］.地質(zhì)科技情報，1988，7（1）:127-136.

［11］朱建東，池順都.新疆阿舍勒銅礦GIS預測時的找礦有利度和空間相關(guān)性分析［J］.地質(zhì)與勘探，2006，42（2）:60-64.

［12］王永立.天津北部山區(qū)遙感找礦信息提取及找礦遠景分析［J］.地質(zhì)調(diào)查與研究，2012，35（3）：229-235.

［13］席斌斌，張德會，周利敏，張文淮，王成.江西省全南縣大吉山鎢礦成礦流體演化特征 ［J］.地質(zhì)學報，2008，82（7）: 956-966.

［14］程順波，付建明，馬麗艷，盧友月.南嶺地區(qū)成鎢、成錫花崗巖組合的幾個判別標志 ［J］.華南地質(zhì)與礦產(chǎn)，2014，30 （4）:352-360.

［15］趙希剛，吳漢寧，柏冠軍，王靖華.重磁異常解釋斷裂構(gòu)造的處理方法及圖示技術(shù) ［J］.地球物理學進展，2008，23 （2）:414-421.

［16］翟裕生，姚書振，崔彬，程小久，秦長興，金偉.成礦系列研究［M］.武漢:中國地質(zhì)大學出版社，1995:59-71.

［17］蔣小芳.川滇黔地區(qū)玄武巖銅礦資源評價（碩士論文）［D］.武漢:中國地質(zhì)大學（武漢），2006.

［18］劉世翔，薛林福，郄瑞卿，張旭東，孟慶龍.基于GIS的證據(jù)權(quán)重法在黑龍江省西北部金礦成礦預測中的應用［J］.吉林大學學報（地球科學版），2008，37（5）:889-894.

［19］嚴冰，陽正熙，王小春.基于Arc-Wofe系統(tǒng)的證據(jù)權(quán)法在層控型鉛鋅成礦預測中的應用［J］.新疆地質(zhì)，2005，23（3）: 300-303.

［20］Holeneus D E，Raines G L，Causey J D，Bookstrom A A，F(xiàn)rost P T and Hyndman P C.Assessment method for epithermal gold deposits in northeast Washingtong State:Using weights of evidence GIS modeling［M］.U.S.Geological Survey，2001:30-39.

Lei T C，Luo S X and Wang L.Application of weights of evidence method for prediction of W-Sn polymetallic deposits based on multi-source information in Nanling metallogenic belt.Geology abd Mineral Resources of South China ，2016，32（1）：34-42.

LEI Tian-Ci，LUO Shi-Xin，WANG Lei

（Wuhan Centre of China Geological Survey，Wuhan 430205，China）

According to1:500000 geological database，geophysical and geochemical exploration data in Nanling metallogenic belt，all kinds of related maps were compiled by data-processing，analyzing and interpreting.On basis of the GIS platform，guided by the geological anomaly theory，all favorable ranges of prospecting factors were selected through multi-information metallogenic coupling analysis and anomalous information extraction.Finally，13 effective evidence layers were built.Ore-forming quantitative prediction and assessment was done based on weights of evidence of Morpas and Arc-wofe model.According to the critical value of 0.024、0.078 and 0.098，41 for prospecting，70 secondary prospect areas and 85 triple prospect areas were submited，it was significanct to plan work and plan mineral resources in the region.

Nanling metallogenic belt；benefits；unit of mineral equivalent；weights of evidence；metallogenic prediction

中圖分類法：P622+.6；P618.44；P618.67A

1007-3701（2016）01-034-09

10.3969/j.issn.1007-3701.2016.01.005

2015-12-30；

2016-02-09.

中國地質(zhì)調(diào)查局項目“武當-桐柏-大別成礦帶多元信息提取及深部找礦評價”（12120113068900）.

雷天賜（1977-），男，高級工程師，從事地、物、化、遙等方面的生產(chǎn)與科研工作，E-mail:leitianci@126.com.