寧墨奐，胡昌松，溫春齊，周 玉，3，費(fèi)光春，何陽陽，李 丹

(1.成都理工大學(xué)地球科學(xué)學(xué)院，四川成都 610059;2.重慶市國土局地質(zhì)調(diào)查院，重慶 401122;3.中國地質(zhì)科學(xué)院礦產(chǎn)綜合利用研究所，四川成都 610041)

1 引言

綜合信息成礦預(yù)測的理論與方法是在找礦新形勢下孕育而生的，國際上實(shí)施了“礦產(chǎn)資源評價(jià)中計(jì)算機(jī)應(yīng)用標(biāo)準(zhǔn)”，推出區(qū)域價(jià)值估計(jì)法、體積估計(jì)法、豐度估計(jì)法、礦床模型法、德菲爾法和綜合方法等六種資源預(yù)測的標(biāo)準(zhǔn)方法(侯翠霞等，2010)，還有Bothol等提出的特征分析法和P.M·康期坦丁諾夫等提出的邏輯信息法，以及Agterberg(1993)提出的證據(jù)權(quán)重法等，特別受到我國地質(zhì)人員的重視?，F(xiàn)在多元統(tǒng)計(jì)方法和計(jì)算機(jī)技術(shù)已廣泛應(yīng)用于成礦預(yù)測中。

多(不雜)-(波)龍礦集區(qū)位于班公湖-怒江縫合帶中西段北側(cè)，羌塘-三江復(fù)合板片南緣，目前已發(fā)現(xiàn)多不雜、波龍、尕爾勤等斑巖型銅礦。但由于交通不便，工作環(huán)境惡劣，地質(zhì)基礎(chǔ)薄弱，以往僅提出過區(qū)域內(nèi)的找礦模型(李玉彬等，2012)，本次研究是在總結(jié)前人工作基礎(chǔ)上，通過對區(qū)內(nèi)相關(guān)典型礦床成礦規(guī)律的充分研究，并在建立綜合信息找礦模型的基礎(chǔ)上，對地、物、化、遙多源信息進(jìn)行優(yōu)化處理，實(shí)現(xiàn)了找礦地質(zhì)特征與找礦數(shù)學(xué)模型的關(guān)聯(lián)與轉(zhuǎn)化，再使用數(shù)學(xué)模型對其關(guān)聯(lián)度求解、排序，得出定量的結(jié)果，使此次預(yù)測具有了功能強(qiáng)、優(yōu)選度高、與礦床經(jīng)濟(jì)意義聯(lián)系密切、便于在工作程度較低的地區(qū)開展應(yīng)用等優(yōu)點(diǎn)。

2 礦區(qū)地質(zhì)特征

2.1 地層

研究區(qū)地層出露比較簡單(見圖1)，主要為中生界中侏羅統(tǒng)曲色組(J2q)、色哇組(J2s)、下白堊統(tǒng)美日切組(K1m)，其次為新生界上第三系康托組(N1k)。主要賦礦地層為中侏羅統(tǒng)曲色組(J2q)和色哇組(J2s)，其中曲色組(J2q)總體上呈近EW向展布，傾向NNE，傾角50°～80°，巖性為長石石英砂巖、粉砂質(zhì)板巖夾硅質(zhì)巖、花崗閃長斑巖、灰綠色玄武巖、基性火山熔巖、英安巖，厚度大于3000m，與下覆地層呈整合接觸。色哇組(J2s)巖性為長石石英砂巖、巖屑砂巖，砂礫巖夾深灰色至深黑色粉砂質(zhì)板巖，有閃長巖瘤、英安巖、流紋巖、火山角礫巖出露，傾角45°～75°，厚度大于3000m，與下覆地層呈整合接觸。

2.2 巖漿巖

研究區(qū)巖漿活動(dòng)十分頻繁、強(qiáng)烈，巖漿活動(dòng)總體上以噴發(fā)、噴溢及超淺成侵入為主，基性、中酸性、酸性巖體均有出露，規(guī)模一般較小，往往是成帶狀、串珠狀展布，成群出現(xiàn)，受斷裂構(gòu)造的控制明顯，并具多期活動(dòng)特征。

花崗閃長斑巖為區(qū)內(nèi)主要含礦斑巖，其呈小巖株分布于在多不雜、波隆、尕爾勤、拿若、拿頓等處，出露面積0.1～0.5km2不等，巖體蝕變強(qiáng)烈，具鉀化、硅化、綠泥綠簾石化、碳酸鹽化，黃鐵礦、黃銅礦化普遍。巖石多呈塊狀構(gòu)造，巖體邊緣有角礫狀構(gòu)造;多為斑狀結(jié)構(gòu)、粒狀結(jié)構(gòu)，也有呈顯微嵌晶結(jié)構(gòu)、蠕蟲狀結(jié)構(gòu)、顯微文象結(jié)構(gòu)以及局部的隱晶質(zhì)結(jié)構(gòu)等。其中斑晶含量約占55%，粒經(jīng)約2～5mm。斑晶主要成分為斜長石，其次為石英、鉀長石、角閃石、黑云母等。

含礦斑巖 SHRIMP鋯石 U-Pb年齡127.8±2.6Ma(曲曉明等，2006)和121.6±1.9Ma(李金祥等，2008)，表明成巖成礦時(shí)代為早白堊世，屬于特提斯洋閉合碰撞階段的產(chǎn)物。據(jù)化學(xué)分析結(jié)果分析，含礦花崗閃長斑巖 SiO2含量為 62.24% ～64.92%，Al2O314.64% ～15.58%，K2O含量為3.18% ～6.60%，Na2O含量為0.44% ～0.79%，屬鈣堿性巖石類型，高鉀高鋁花崗巖。

2.3 構(gòu)造

由于歷次構(gòu)造運(yùn)動(dòng)及巖漿活動(dòng)，區(qū)內(nèi)形成有大量的次級斷層，研究區(qū)南部的班公錯(cuò)-康托-茲格塘錯(cuò)斷裂為一超殼斷裂，是羌塘-三江復(fù)合板片和岡底斯-念青唐古拉板片的分界斷裂，也是班公錯(cuò)-怒江縫合帶的北界斷裂，受其影響，在工作區(qū)內(nèi)形成大量的次級斷層。

工作區(qū)次級斷裂構(gòu)造具長期性、多期次活動(dòng)特征，總體有三組∶①近東西向斷裂構(gòu)造為礦區(qū)內(nèi)規(guī)模大、形成時(shí)間最早、活動(dòng)時(shí)間最長、具有壓扭性特征的斷裂構(gòu)造;②北東向斷裂晚于東西向斷裂，并切割東西向斷裂;③北西向斷裂是形成最晚的一組斷裂，并切割北東向斷裂。這三組斷裂構(gòu)造形成似菱形狀格架。

圖1 多龍礦集區(qū)地質(zhì)簡圖Fig.1 Geological sketch map of the Duolong ore concentration area

3 多元信息特征及變量選擇

通過對多龍礦集區(qū)內(nèi)具有代表性的三個(gè)礦床——多不雜、波龍、尕爾勤銅礦床的地質(zhì)成礦條件進(jìn)行詳細(xì)研究，總結(jié)了如下成礦-找礦規(guī)律(寧墨奐，2012)。

3.1 含礦花崗閃長斑巖特征及其與銅礦化的關(guān)系

多龍礦集區(qū)斑巖銅礦中，與成礦密相關(guān)的巖體均為花崗閃長斑巖?；◢忛W長斑巖出露面積較小，多不雜 0.18km2，波龍 0.1km2，尕爾勤 0.2km2，與圍巖呈侵入接觸。含礦花崗閃長斑巖呈灰白色，塊狀構(gòu)造，斑狀結(jié)構(gòu)，含礦巖石中原生礦物蝕變嚴(yán)重，主要表現(xiàn)為斜長石蝕變?yōu)殁涢L石、綠簾石、高嶺石;鉀長石蝕變?yōu)檎惩恋V物，角閃石、黑云母蝕變?yōu)榫G泥石，氧化為褐鐵礦。而較為新鮮的巖石通常不含礦，或者品位較低。

與礦化有關(guān)的巖體由多個(gè)小巖株(枝)構(gòu)成，它們在深部連接在一起，構(gòu)成一個(gè)形態(tài)不規(guī)則的柱狀礦體，由細(xì)粒浸染狀、細(xì)脈-網(wǎng)脈狀金屬硫化物礦石組成。在同一巖體中，中上部以星點(diǎn)狀、細(xì)脈狀黃鐵礦化為主，偶見石英-黃銅礦礦脈，表現(xiàn)為細(xì)脈狀的黃鐵礦化;巖體中下部以稀疏浸染狀之黃銅礦化、黃鐵礦化為主，黃銅礦多于黃鐵礦，伴生有石英脈、石英-石膏脈產(chǎn)出的黃銅礦化。主成礦元素Cu在走向和垂向上變化不大，總體上礦化較均勻，統(tǒng)計(jì)平均品位為0.422×10-2，品位變化系數(shù)為63.4%。礦(化)體分布在斑巖體及圍巖中，賦礦巖石以花崗閃長斑巖為主，礦體厚度同樣較穩(wěn)定，其變化系數(shù)為42.30%。含礦斑巖SHRIMP鋯石U-Pb年齡均為早白堊世，表現(xiàn)為富鉀貧鈉，堿質(zhì)含量較高，屬于高鋁高鉀鈣堿性系列類型。

3.2 構(gòu)造與銅礦化體關(guān)系

構(gòu)造對銅礦化體的控制作用表現(xiàn)在以下幾個(gè)方面，熱液礦床的形成除含礦熱液本身外，首先需要其從來源地點(diǎn)到達(dá)沉淀地點(diǎn)的通道，其次礦液沉淀有賴于一定種類的構(gòu)造裂隙作為有利的空間。常?？梢钥吹揭粋€(gè)礦體規(guī)模的大小，礦石質(zhì)量的優(yōu)劣，取決于容礦構(gòu)造因素的優(yōu)勢(陳國達(dá)，1983)。此外，作為礦液來源的巖漿活動(dòng)，又往往與構(gòu)造運(yùn)動(dòng)具有密切的關(guān)系，因此構(gòu)造對于熱液礦床起到了直接和間接的雙重控制作用。

礦床研究區(qū)最主要的構(gòu)造為南部的班公錯(cuò)-康托-茲格塘錯(cuò)超殼斷裂，該斷裂是羌塘-三江復(fù)合板片和岡底斯-念青唐古拉板片的分界斷裂，也是班公錯(cuò)-怒江縫合帶的北界斷裂，該深大斷裂為深部巖漿提供了通道，可以視為本地區(qū)的導(dǎo)礦構(gòu)造。受該斷裂影響，在工作區(qū)內(nèi)形成大量的次級斷層，眾多的次級斷裂同樣為成礦流體的運(yùn)移提供了通道。

目前研究區(qū)內(nèi)比較具有代表性的三個(gè)礦床都和次級構(gòu)造有一定的關(guān)系，其中多不雜礦區(qū)有近東西向F2、北東向F10兩組斷層;波龍礦區(qū)位于北東向F2和北西西向斷裂F14的復(fù)合交叉點(diǎn)(礦區(qū)北東側(cè));尕爾勤礦區(qū)F1斷層從礦區(qū)北部通過，推測其次級斷裂從礦區(qū)中部穿過。在這些次級構(gòu)造交匯部位，局部空間膨大有利于礦質(zhì)沉淀，斷層部位的巖石破碎有利于熱液交代，在后巖漿成礦階段，富含成礦物質(zhì)的成礦流體在這些成礦有利部位發(fā)生沉淀、交代，形成品位較高的富礦體。另外區(qū)域的構(gòu)造活動(dòng)是含礦熱流體運(yùn)移的重要驅(qū)動(dòng)力因此區(qū)內(nèi)的斷裂構(gòu)造對斑巖銅礦的成礦具有重大的影響。故斷裂構(gòu)造交匯部位及其斷裂帶附近是尋找斑巖的有利部位。

3.3 礦體圍巖蝕變特征及遙感特征

多龍礦集區(qū)各礦化斑巖均有較明顯的礦化蝕變，多不雜礦區(qū)自斑巖體到圍巖可劃分為鉀化帶-石英絹云母化帶-青磐巖化帶。主要銅礦化部位為鉀化帶-石英絹云母帶，該帶中金屬硫化物含量高;波龍礦區(qū)蝕變在地表表現(xiàn)為中心式面型分布，以發(fā)育強(qiáng)烈的鉀化和鉀硅化為特征，以礦體為中心，由內(nèi)向外依次可劃分出鉀硅化帶-絹英巖化帶-青磐巖化帶，泥化僅在局部可見，基本上不存在明顯的泥化帶;尕爾勤銅礦蝕變分帶同多不雜銅礦較為相似，可根據(jù)蝕變類型由蝕變中心往外分為鉀硅化帶-絹英巖化帶-粘土化帶-青磐巖化帶。該地區(qū)斑巖體圍巖蝕變最外圍均存在褐鐵礦化帶，呈環(huán)帶狀分布青磐巖化帶外側(cè)，由遠(yuǎn)離巖體的變長石石英砂巖、泥砂質(zhì)板巖構(gòu)成，褐鐵礦化普遍而強(qiáng)烈，形成大面積的“火燒皮”，具有很強(qiáng)的識別度。

在遙感影像圖上，礦化(蝕變)信息明顯，由于酸性斑巖體發(fā)生高嶺石化等次生變化，色調(diào)比周圍淺，另外還表現(xiàn)為淺黃等色調(diào)信息異常，該種色調(diào)異常是由硫化礦體或硫化物經(jīng)氧化后形成的氧化帶和褐鐵礦化現(xiàn)象引起的，它是尋找硫化礦床和硫化物有關(guān)的銅礦指示標(biāo)志。

3.4 地層接觸帶與銅礦化的規(guī)律

雖然斑巖銅礦成礦物質(zhì)主要來源于斑巖體，但地層的巖性同樣對成礦起到重要的作用。如:地層巖石完整性較好、隔水性好，則可以對成礦熱液起到遮擋作用，使其成礦相對集中;巖石較破碎、隔水性差，則成礦熱液通過裂隙運(yùn)移成礦范圍較大，成礦流體呈現(xiàn)多來源，形成的礦石種類也相對較為復(fù)雜。礦區(qū)地層主要巖石類型有泥巖、粉砂巖、砂巖、玄武巖、安山巖、流紋巖、火山角礫巖等，礦體的圍巖主要為變質(zhì)長石石英砂巖，次為石英粉砂巖。多不雜、波龍、尕爾勤礦區(qū)附近均有美日切組出露，由于不同的地層巖性差距大，其受力變形也不一致容易形成層間裂隙、斷層，不同巖性的層理面作為構(gòu)造薄弱面可以作為含礦熱液流通的通道，因此不同的地層組合熵對成礦的有利度也不一樣，即地層組合越復(fù)雜越易于成礦。

3.5 銅礦化的物探異常特征

目前區(qū)內(nèi)開展的地球物理找礦手段為磁法和電法找礦，眾所周知，磁異常是地下居里深度以上磁性不均勻的反映，包含著地下目標(biāo)體和非目標(biāo)體的綜合信息(宋新華等，2010)。研究區(qū)內(nèi)的主要巖石具有比較明顯的磁性差異，這些差異正是是磁法找礦的基礎(chǔ)。研究區(qū)內(nèi)主要巖石類型的磁異常為:安山巖、閃長巖、表生礦石具有正磁性，變質(zhì)砂巖、閃長斑巖均為負(fù)磁性。對區(qū)內(nèi)具代表性的三個(gè)礦床(點(diǎn))的磁異常做對比，我們發(fā)現(xiàn)，波龍和多不雜礦區(qū)磁異常非常相似，均為范圍小的串珠狀弱正磁異常，并相鄰有同樣小范圍、低強(qiáng)度的負(fù)磁異常，顯示出該處深部的小型斑巖體和上部被氧化礦體的特征。而尕爾勤礦區(qū)呈現(xiàn)出面積較大的弱強(qiáng)磁異常主要是由于其附近地層中有玄武巖出露，形成較大面積弱正磁異常。

區(qū)內(nèi)多為隱伏礦體，在找礦過程中，如何布置鉆孔一直是地質(zhì)找礦工作中的難點(diǎn)(李帝銓等，2007)。由于隱伏礦體埋深較大、信息弱和干擾大，要求傳統(tǒng)的找礦技術(shù)方法與手段在探測深度、探測靈敏度以及抗干擾方面有較大的改進(jìn)和提高，近年來現(xiàn)代電子技術(shù)、檢測技術(shù)和基礎(chǔ)加工技術(shù)的提高，電法找礦效果有了較大幅度的提升(張作倫等，2008;馬庚杰等，2007;劉國印等，2008;歐陽南，2009)。

區(qū)內(nèi)各巖石電性參數(shù)中電阻率大致可以分為高電阻率、中等電阻率和低電阻率三類。高電阻率的巖石有輝綠巖、安山玢巖、玄武質(zhì)安山巖、花崗斑巖等，低電阻率巖石有變長石石英砂巖、變質(zhì)砂巖、角礫巖、第四紀(jì)殘坡積物等;其它各種巖礦(化)石電阻率屬于中等。各類巖石極化率較低，且比較穩(wěn)定，平均極化率都小于3%;而孔雀石化、黃鐵礦化的巖石極化率略高一些，一般在3% ～6%之間;隨著金屬硫化物含量的增多，極化率在增高，最高可達(dá)22.04%。

結(jié)合本人在西藏邦鋪鉬銅多金屬礦找礦預(yù)測的經(jīng)驗(yàn)(寧墨奐等，2010)，在極化率最高的部位通常無明顯銅礦化，在中高極化率的部位施工鉆孔反而可以取得較好的見礦效果。使用激發(fā)極化法圈定礦體的找礦特征為高電阻率中高激化率。

3.6 銅礦床的地球化學(xué)異常特征

研究區(qū)晝夜溫差大，海拔高，植被稀少，主要風(fēng)化作用為物理風(fēng)化作用，區(qū)內(nèi)相對高差小，地表徑流相對缺乏，風(fēng)化產(chǎn)物搬運(yùn)作用弱，搬運(yùn)距離近，在該地區(qū)化探找礦具有較好的效果，2008年西藏自治區(qū)地質(zhì)調(diào)查院在該區(qū)開展了“西藏班公湖-怒江成礦帶西段銅多金屬資源調(diào)查”工作，在多不雜及周邊地區(qū)開展了1/5萬水系沉積物測量工作，選取Cu、Mo、W、Ag、Au、Bi、Zn、Cd、Pb、Co、Hg 作為異常元素。

在多不雜礦區(qū)，各元素異常套合較好(圖2)，Cu、Mo、W、Bi具濃度分帶，其中銅異常強(qiáng)度高，有六級濃度分帶，異常形態(tài)較規(guī)則，Mo、W、Bi元素異常具外、中、內(nèi)帶，強(qiáng)度大，規(guī)模大，內(nèi)帶就是礦體露頭部位，與地質(zhì)特征對應(yīng)性好。波龍礦區(qū)出現(xiàn)異常的元素為Cu、Mo、W、Bi、Au、Ag、Zn、Cd、Pb、Co、Hg，平面呈不規(guī)則橢圓狀。其中Cu、Bi具有四級濃度分帶，Mo、Zn具有三級濃度分帶，Ag、W具有二級濃度分帶，該異常與已發(fā)現(xiàn)的隱伏含礦斑巖體特征完全吻合，特別是Cu元素異常特征，異常形態(tài)較規(guī)則，套合較好。

異常規(guī)律為:礦體上方出現(xiàn)了斑巖銅礦床特征的近礦指示元素Mo、W、Bi異常和銅的前緣指示元素Pb、Zn異常，異常規(guī)模大，含量高，連續(xù)性好，具有明顯的濃集中心。銅異常下限＞2000×10-6，并具有多級濃度分帶的部位，與鉆孔見礦對應(yīng)性好。

3.7 礦體氧化帶特征

地質(zhì)工作遵循由表及里，由已知到未知的工作方法，地表填圖通常是地質(zhì)找礦最先開展的工作，地表的礦化部位也是實(shí)施深部工程的重點(diǎn)區(qū)域，故礦體氧化帶特征具有重要的找礦意義。在尕爾勤和多不雜礦區(qū)，地表可見銅礦體的表生氧化帶，其代表性金屬礦物為孔雀石、藍(lán)銅礦，孔雀石呈翠綠色，葡萄狀、皮殼狀產(chǎn)出，藍(lán)銅礦呈深藍(lán)色，多呈包膜狀產(chǎn)出?？兹甘退{(lán)銅礦可以作為尋找斑巖銅礦體的一個(gè)找礦標(biāo)志。

圖2 多不雜礦區(qū)Cu、Mo、W、Bi化探異常圖②Fig.2 Map showing geochemical anomalies of the Duobuza deposit②

4 找礦模型的建立

多源地學(xué)信息包括地質(zhì)礦產(chǎn)資料、物化探資料、遙感資料等。對這些資料進(jìn)行信息提取，獲得遙感異常、化探異常、物探異常、控礦地層、控礦構(gòu)造、與礦化有關(guān)的巖漿巖、礦床(點(diǎn))、礦化蝕變帶等異常信息，進(jìn)一步研究異常(包括單一異常、組合異常)與成礦的關(guān)系，獲得區(qū)域找特征，并總結(jié)找礦模型，進(jìn)行找礦預(yù)測。

綜合多龍礦集區(qū)各礦體的含礦巖體類型、構(gòu)造關(guān)系、圍巖蝕變、地層巖石類型、物化探異常等各種找礦信息，構(gòu)成了一個(gè)統(tǒng)一的受深部地質(zhì)構(gòu)造制約的多龍礦集區(qū)斑巖銅礦找礦模型(表1)。

表1 多龍礦集區(qū)斑巖銅礦多元信息找礦模型Table 1 A multi-information prospecting model of porphyry copper ore deposits in the Duolong ore concentration area

5 加權(quán)特征分析與找礦靶區(qū)定量預(yù)測

5.1 加權(quán)特征分析法

找礦預(yù)測方法種類繁多，目前還沒有找到普遍行之有效的統(tǒng)一方法，故應(yīng)根據(jù)具體地質(zhì)情況和預(yù)測對象、預(yù)測過程的不同階段及研究程度等因素合理選擇預(yù)測方法。與基于GIS的成礦預(yù)測方法相比，加權(quán)特征分析法對硬件及對計(jì)算機(jī)操作能力要求較低，方便實(shí)用，在地質(zhì)基礎(chǔ)薄弱的地區(qū)，工作在找礦第一線的地質(zhì)單位特別適用，并在找礦領(lǐng)域取得巨大成功。一般特征分析使用的變量是離散型的邏輯變量，把找礦的綜合信息轉(zhuǎn)化為二態(tài)(1，0)或三態(tài)(1，0，-1)賦值的變量，以便地學(xué)數(shù)據(jù)信息提取，并利用數(shù)學(xué)方法進(jìn)行區(qū)域成礦預(yù)測(Carranz et al.，2002)。

成礦預(yù)測加權(quán)特征分析與一般特征分析的不同之處就在于“加權(quán)”，權(quán)值在各種礦床類型或勘查程度不相同的礦區(qū)可以具有不同的地質(zhì)意義?？辈槌潭容^高的地區(qū)可以用單礦體儲(chǔ)量作為權(quán)值，是顯示礦床、礦體成礦規(guī)模大小的權(quán)(李英杰，2007)。但是在地質(zhì)勘查程度較低的地區(qū)，如本研究區(qū)，無法滿足該這種應(yīng)用條件。

本次研究中，用具代表性的單工程的平均品位×礦體厚度來定義權(quán)，稱為品位權(quán)。由于數(shù)學(xué)模型在變量匹配矩陣中增加了品位權(quán)，使加權(quán)特征分析模型在這種地質(zhì)基礎(chǔ)薄弱的地區(qū)，最大限度的實(shí)現(xiàn)了功能強(qiáng)、優(yōu)選度高、與礦床經(jīng)濟(jì)意義聯(lián)系密切的特點(diǎn)。

5.2 數(shù)學(xué)模型及計(jì)算方法

一個(gè)模型單元或預(yù)測單元，有多個(gè)成礦-找礦的綜合信息地質(zhì)變量，但對控制礦床的形成或指示礦床的有無而論，各成礦-找礦的綜合信息地質(zhì)變量的重要性或權(quán)是不等的。在已知有礦的模型統(tǒng)計(jì)單元中選取有代表性的典型單元n個(gè)，在每個(gè)單元內(nèi)取m個(gè)狀態(tài)變量(l，0)，用乘積矩陣求變量權(quán)。

(1)原始數(shù)據(jù)及加權(quán)矩陣

設(shè)取n個(gè)模型統(tǒng)計(jì)單元，m個(gè)變量，n組單工程品位，用離散型的變量數(shù)據(jù)，建立了二個(gè)數(shù)據(jù)陣。

原始數(shù)據(jù)陣:

其中i為模型單元，j為選取的變量，而 X=［Xji］(n × m)

品位對數(shù)加權(quán)矩陣:

ndi表示第i各單元的品位對數(shù)值。

(2)變量權(quán)及關(guān)聯(lián)度計(jì)算

相應(yīng)品位加權(quán)矩陣:Z=X’［D］2X

利用品位加權(quán)乘積矩陣，分別計(jì)算變量權(quán):aj=為矩陣Z中橫列數(shù)值)。

5.3 多龍礦集區(qū)找礦靶區(qū)定量預(yù)測優(yōu)選

取具代表性的見礦工程10個(gè)，選取7個(gè)變量，運(yùn)用離散型二態(tài)0、1變量數(shù)據(jù)，將10個(gè)見礦工程單元7個(gè)變量以及每個(gè)見礦工程的品位和自然對數(shù)分別列入表2。

表2 多龍礦集區(qū)斑巖銅礦原始數(shù)據(jù)表Table 2 Original data of porphyry copper ore deposit in the Duolong metallogenic district

計(jì)算儲(chǔ)量加權(quán)矩陣:

計(jì)算變量權(quán)，為便于計(jì)算，所得結(jié)果除以100得:

各找礦靶區(qū)關(guān)聯(lián)度計(jì)算:

現(xiàn)將各找礦靶區(qū)地質(zhì)特征列于表3，以便通過公式計(jì)算關(guān)聯(lián)度。

關(guān)聯(lián)度計(jì)算:

通過上述計(jì)算，求出關(guān)聯(lián)度排序，分別是f4＞f3＞f2＞f7＞f5＞f1＞f6，關(guān)聯(lián)度排名最高的靶區(qū)為 NR，關(guān)聯(lián)度21.09，其次為SJ，關(guān)聯(lián)度19.43，再次為SN，關(guān)聯(lián)度17.76，上述三個(gè)靶區(qū)具有很好的找礦潛力，在上述靶區(qū)內(nèi)采用高精度的物探方法查明礦體大概的方位、產(chǎn)狀，不僅可以提高找礦效率，還能夠節(jié)約經(jīng)費(fèi)，指導(dǎo)鉆孔布設(shè)(沈遠(yuǎn)超等，2008;王慶乙等，2009;曹新志，2009)。

表3 找礦靶區(qū)地質(zhì)特征數(shù)據(jù)表Table 3 Geological feature data of prospecting target zones

6 結(jié)論

通過大量野外實(shí)地觀察和室內(nèi)綜合分析研究，建立了研究區(qū)內(nèi)的斑巖型銅礦的成礦-找礦模型，并且針對區(qū)內(nèi)地質(zhì)工作基礎(chǔ)較差的現(xiàn)實(shí)情況，首次將單工程品位和礦體厚度之積作為權(quán)值應(yīng)用到多龍礦集區(qū)找礦靶區(qū)預(yù)測中，優(yōu)選出10個(gè)具代表性的鉆孔作為模型單元，對7個(gè)未知預(yù)測單元進(jìn)行了定量預(yù)測，最終確定了最優(yōu)找礦靶區(qū)，指出了進(jìn)一步找礦方向。

將加權(quán)特征分析方法應(yīng)用到多龍礦集區(qū)找礦靶區(qū)預(yù)測中，既是在一定程度上豐富和發(fā)展了找礦靶區(qū)定量預(yù)測方法，同時(shí)也為多龍礦集區(qū)，甚至班公湖-怒江成礦帶這類地質(zhì)工作薄弱的地區(qū)的找礦預(yù)測提供了新的認(rèn)識和經(jīng)驗(yàn)，具有較為重要的理論與實(shí)際意義。

［注釋］

① 西藏自治區(qū)地質(zhì)礦產(chǎn)勘查開發(fā)局第五地質(zhì)大隊(duì).2010.西藏自治區(qū)改則縣多不雜銅礦2009年度普查報(bào)告［R］.

②魏啟榮，任利民，樊俊青，魏俊浩，楊奇荻，郝偉，梁云漢，孫驥，張小強(qiáng)，王明志，付樂兵，李歡，秦雅東，賀新星，謝瑜.2009.西藏1:5萬班公湖-怒江成礦帶西段多龍地區(qū)四幅(I-44E019022、I-44E019023/I-44E020022、I-44E020023)900km2土壤與水系沉積物地球化學(xué)調(diào)查報(bào)告［R］.

Agterberg F P，Bonham-Carter G F，Cheng Q M.1993.Weights of evidence modeling and weighted logistic regression for mineral potential mapping［C］.Computers in Geology:25 Years of Progress.New York:Oxford University Press:13-32

Agterberg F P，Cheng Q M，Wright D F.1993.Fractal modeling of mineral deposits［C］.Proceedings，APOM XXIV international symposium on the application of computers and operations research in the mineral industries.Montreal，Quebec，Canada:Canadian Institute of Mining and Metallurgy:43-53

Cao Xin-zhi，Zhang Wang-sheng，Sun Hua-shan.2009.Progress in the Study of Deep Exploration in China［J］.Geological Science and Technology Information，28(5):104-109(in Chinese with English Abstract).

Carranza E J M，Hale M.2002.Where are porphyry copper deposits spatially localized?A case study in Berguet province，Philippines［J］.Natural Resource Research，l1(1):45-60

Chen Guo-da.1983.The research method of Constraints on metallogenic tectonic settings［M］.Beijing:Geological Publishing House:1-412(in Chinese)

Hou Cui-xia，Liu Xiang-chong，Zhang Wen-bing，Shen Wei.2010.New method and theory of metallogenic prediction［J］.Geological Bulletin of China，29(6):953-960(in Chinese with English Abstract)

Li Di-quan，Wang Guang-jie，Di Qing-yun，F(xiàn)u Chang-ming，Shi Kun-fa，Li Ying-xian.2007.The application of high-power induced polarization in the middle section of eerguna metallogenic belt［J］.Progress in Geophysics，(22)5:1621-1626(in Chinese with English Abstract)

Li Jin-xiang，Li Guang-ming，Qin Ke-zhang，Xiao Bo.2008.Geochemistry of porphyries and volcanic rocks and ore-forming geochronology of Duobuza gold rich porphyry copper deposit in Bangonghu belt，Tibet:Constraints on metallogenic tectonic settings［J］.Acta Petrologica Sinica，24(3):531-543(in Chinese with English Abstract)

Li Ying-jie.2007.Application of weighted characteristic analysis in target prospecting in the Dashui gold deposit［D］.Shijiazhuang:Shijiazhuang University of Economics:1-56(in Chinese with English Abstract)

Li Yu-bin，Duo Ji，Zhong Wan-ting，Li Yu-chang，Qiangba Wang-dui，Chen Hong-qi，Liu Hong-fei，Zhang Jin-shu ，Zhang Tian-ping ，Xu Zhi-zhong，F(xiàn)an An-hui，Suolang Wang-mu.2012.An exploration model of the Duobuza porphyry Cu-Au deposit in Gaize Country，northern Tibet［J］.Geology and Exploration，48(2):274-87(in Chinese with English Abstract)

Liu Guo-yin，Yan Chang-hai，Zhao Jian-min，Wang Ji-zhong，Li Zhongmng.2008.Applications of combined microgravitysurvey and CSAMT sounding in search for concealed bauxite deposits in western Henan，China［J］.Geological Bulletin of China，27(5):641-648(in Chinese with English Abstract)

Ma Geng-jie，Zhang Rui，Zheng Xiao-li，Liu Hong-tao.2007.Application of VLF-EM method in the mineral exploration of Longtoushan Ag-Pb-Zn Deposit，Alukerqinqi，Inner Mongolia［J］.Progress in Geophysics，(22)5:1627-1636(in Chinese with English Abstract)

Ning Mo-huan，Wen Chun-qi.2010.Geological and geophysical characteristics and analysis on enlarging the prospect for prospecting work of Tibetan Bangpu porphyry-type molybdenum-copper mine area［J］.Mineral Resources and Geology，(24)6:542-503(in Chinese with English Abstract)

Ning Mo-huan.2012.Prospecting prediction and Characteristics of por-phyry copper deposit in the Duolong metallogenic district，Tibet［D］.Chengdu:Chengdu University of Technology:1-78(in Chinese with English Abstract)

Ouyang Nan.2009.Analysis of Anomalies of No 2，3，4 Low Resistance on EH-4 High Frequency Magnetotelluric Depth in Deep Prospecting in Tongshankou Ore District［J］.Resources Environment＆ Engineering，34(4):500-503(in Chinese with English Abstract)

Qu Xiao-ming，Xin Hong-bo.2006.Ages and tectonic environment of the Bangong Co porphyry copper belt in western Tibet，China［J］.Geological Bulletin of China，25(7):792-799(in Chinese with English Abstract)

Song Xin-hua，Yin Bing-xi，Yan Hong，Yan Shu-qiang.2010.The Relationship between the Aeromagnetic weak anomalies and the polymetallic ore deposits(ore spots)in the Beishan mountain，Weining［J］.Geophysical and Geochemical Exploration，34(3):289-293(in Chinese with English Abstract)

Shen Yuan-chao，Shen Ping，Liu Tie-bing，Li Guang-ming，Zeng Qingdong.2008.Prediction of hidden gold orebodies in depleted mines by the Stratagem EH4 system［J］.Progress in Geophysics，23(1):55-56(in Chinese with English Abstract)

Wang Qing-yi，Li Xue-sheng，Xu Li-zhong.2009.The High-precision borehole three-component magnetic system:An effective tool in search for deep deposits in mines［J］.Geophysical and Geochemical Exploration，33(3):235-244(in Chinese with English Abstract)

Zhang Zuo-lun，Zeng Qing-dong，Ye Jie，Jia Chang-shun，Li Wentao.2008.Application of VLF-EM in forecasting of ore exploration［J］.Geology and Prospecting，44(1):67-69(in Chinese with English Abstract)

［附中文參考文獻(xiàn)］

曹新志，張旺生，孫華山.2009.我國深部找礦研究進(jìn)展綜述［J］.地質(zhì)科技情報(bào)，28(5):104-109

陳國達(dá).1983.成礦構(gòu)造研究法［M］.北京:地質(zhì)出版社:1-412

侯翠霞，劉向沖，張文斌，申維.2010.成礦預(yù)測理論與方法新進(jìn)展［J］. 地質(zhì)通報(bào)，29(6):953-960

李帝銓，王光杰，底青云，付長民，石昆法，李英賢.2007.大功率激發(fā)極化法在額爾古納成礦帶中段找礦中的應(yīng)用［J］.地球物理學(xué)進(jìn)展，(22)5:1621-1626

李金祥，李光明，秦克章，肖 波.2008.班公湖帶多不雜富金斑巖銅礦床斑巖-火山巖的地球化學(xué)特征與時(shí)代:對成礦構(gòu)造背景的制約［J］.巖石學(xué)報(bào)，24(3):531-543

李英杰.2007.加權(quán)特征分析在大水金礦靶區(qū)預(yù)測中得應(yīng)用研究［D］.石家莊:石家莊經(jīng)濟(jì)學(xué)院:1-56

李玉彬，多吉，鐘婉婷，李玉昌，強(qiáng)巴旺堆，陳紅旗，劉鴻飛，張金樹，張?zhí)炱剑熘局?，范安輝，索朗旺欽.2012.西藏改則縣多不雜斑巖型銅金礦床勘查模型［J］.地質(zhì)與勘探，48(2):274-287

劉國印，燕長海，趙建敏，王紀(jì)中，李中明.2008.微重力法與可控源音頻大地電磁法組合在豫西尋找隱伏鋁土礦中的應(yīng)用［J］.地質(zhì)通報(bào)，27(5):641-648

馬庚杰，張銳，鄭曉禮，劉紅濤.2007.甚低頻電磁法龍頭山銀鉛鋅礦床勘查中的應(yīng)用［J］.地球物理學(xué)進(jìn)展，22(5):1627-1636

寧墨奐，溫春齊.2010.西藏邦鋪斑巖型鉬銅礦區(qū)地質(zhì)及地球物理特征與擴(kuò)大找礦前景分析［J］.礦產(chǎn)與地質(zhì)，(24)6:542-546

寧墨奐.2012.西藏多龍礦集區(qū)斑巖銅礦床地質(zhì)特征及找礦預(yù)測［D］.成都:成都理工大學(xué):1-78

歐陽南.2009.銅山口礦區(qū)深部找礦EH-4高頻大地電磁測深低阻2、3、4 異常解析［J］.資源環(huán)境與工程，23(4):500-503

曲曉明，辛洪波.2006.藏西班公湖斑巖銅礦帶的形成時(shí)代與成礦構(gòu)造環(huán)境［J］.地質(zhì)通報(bào)，25:792-799

宋新華，尹秉喜，閆紅，顏樹強(qiáng).2010.航磁資料在衛(wèi)寧北山尋找多金屬礦中的應(yīng)用［J］.物探與化探，34(3):289-293

沈遠(yuǎn)超，申 萍，劉鐵兵，李光明，曾慶棟.2008.EH4在危機(jī)礦山隱伏金礦體定位預(yù)測中的應(yīng)用研究［J］.地球物理學(xué)進(jìn)展，23(1):55-56

王慶乙，李學(xué)圣，徐立忠.2009.高精度井中三分量磁測是礦山深部找礦的有效手段［J］.物探與化探，33(3):235-244

張作倫，曾慶棟，葉杰，賈長順，李文濤.2008.甚低頻電磁法在礦體勘查中的應(yīng)用［J］.地質(zhì)與勘探，44(1):67-69