李 希潘振興王 啟高雅寧韓 珂劉樟樂(lè)

(1.長(zhǎng)安大學(xué)地球科學(xué)與資源學(xué)院,陜西西安 710054；2.西北有色地質(zhì)礦業(yè)集團(tuán)有限公司,陜西西安 710054)

高光譜蝕變礦物提取在地質(zhì)找礦中的應(yīng)用
——以吉林某斑巖型鉬礦為例

李 希1，2，潘振興2，王 啟2，高雅寧1，韓 珂1，劉樟樂(lè)2

(1.長(zhǎng)安大學(xué)地球科學(xué)與資源學(xué)院,陜西西安 710054；2.西北有色地質(zhì)礦業(yè)集團(tuán)有限公司,陜西西安 710054)

國(guó)內(nèi)外已廣泛使用高光譜方法提取和識(shí)別巖石和礦物的技術(shù)。本文利用這個(gè)技術(shù)，在地面光譜測(cè)試、數(shù)據(jù)處理和光譜特征分析的基礎(chǔ)上，提取和識(shí)別出斑巖鉬礦區(qū)的蝕變礦物，通過(guò)分析與鉬礦化蝕變的相關(guān)性，從12種主要蝕變礦物中優(yōu)選出5種在礦區(qū)具有代表性的蝕變礦物：絹云母、蒙脫石、伊利石化絹云母、高嶺石(結(jié)晶差)和鐵鎂綠泥石。結(jié)合其空間分布和分帶特征，分析了與鉬礦化的關(guān)系，進(jìn)而劃分出內(nèi)、中和外三個(gè)蝕變帶，預(yù)測(cè)了礦區(qū)的找礦方向，并提供了一種新的找礦思路。

高光譜 光譜特征 斑巖型鉬礦 蝕變礦物

Li Xi, Pan Zhen-xing, Wang Qi, Gao Ya-ning, Han Ke, Liu Zhang-le. Application of alteration minerals extraction based on hyperspectral data to geological prospecting of porphyry molybdenum deposit in Jilin Province [J]. Geology and Exploration,2016,52(3):0489-0496.

0 引言

高光譜蝕變礦物識(shí)別提取和分帶的研究及應(yīng)用在國(guó)內(nèi)已經(jīng)有諸多實(shí)例(甘甫平,2002；雷天賜，2005；周強(qiáng)，2006;鄒林，2006；孫衛(wèi)東，2010)，目前多應(yīng)用于遙感高光譜蝕變礦物識(shí)別提取(田豐，2010；李楠，2010)、高光譜蝕變礦物提取方式方法(雷天賜，2005；吳鳳凰，2013;王亞軍，2015)和對(duì)典型礦物高光譜特性的研究(甘甫平,2003;鄧美容，2010)等方面。已有的高光譜在地質(zhì)礦產(chǎn)方面的應(yīng)用多注重于理論或科研性質(zhì)研究(王晉年,1999；雷天賜，2005；鄧美容，2010；許寧，2013)，少量應(yīng)用于區(qū)域基礎(chǔ)地質(zhì)或礦產(chǎn)調(diào)查(甘甫平,2002；孫衛(wèi)東，2010；朱炳玉，2010；楊志明，2012)。目前，便攜式高光譜礦物分析技術(shù)在加拿大、美國(guó)和澳大利亞等國(guó)已廣泛地應(yīng)用于礦產(chǎn)勘查領(lǐng)域，該方法可以通過(guò)對(duì)礦物種類(lèi)、豐度和成分的識(shí)別，特別是與成礦作用密切相關(guān)的蝕變礦物的識(shí)別有效地圈定熱液礦化蝕變帶，定量或半定量地估計(jì)被測(cè)區(qū)域的蝕變強(qiáng)度和蝕變礦物含量，給出成礦作用規(guī)模和強(qiáng)度等，同時(shí)在輔助勘查決策方面也能提供獨(dú)特的信息。

本文通過(guò)使用便攜式礦物光譜儀對(duì)吉林某斑巖型鉬礦新鮮巖芯進(jìn)行高光譜蝕變礦物提取和空間模型的建立，完成礦化蝕變分帶特征與鉬礦化的關(guān)系，并進(jìn)行找礦預(yù)測(cè)，嘗試高光譜數(shù)據(jù)蝕變礦物提取在地質(zhì)找礦中的應(yīng)用。

1 成礦地質(zhì)環(huán)境

1.1 區(qū)域地質(zhì)背景

礦區(qū)地處華北板塊與佳木斯微板塊之間，張廣才嶺微地塊南緣。地層區(qū)劃屬天山-興安地槽褶皺區(qū)，吉林-延邊分區(qū)的吉林小區(qū)。區(qū)域地層主要有中元古界青龍村群變質(zhì)巖和上元古界一面坡群變質(zhì)巖以及泥盆系和侏羅系火山巖；區(qū)內(nèi)巖漿活動(dòng)時(shí)間漫長(zhǎng)，加里東期、華力西期、印支期、燕山期及喜山期等各個(gè)地質(zhì)時(shí)期巖漿活動(dòng)頻繁，巖漿演化規(guī)律性強(qiáng)，分異明顯；加里東和華力西期構(gòu)造旋回褶皺發(fā)育，軸向北東或北西；斷裂構(gòu)造在區(qū)域發(fā)育，敦化-密山斷裂通過(guò)本區(qū)。

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

礦區(qū)主要出露少量地層和大規(guī)模的巖體。礦區(qū)中部見(jiàn)有少量青龍村群變質(zhì)巖，基本呈近南北向斷續(xù)分布，以黑云石英片巖為主。燕山期似斑狀花崗巖體位于礦區(qū)的東北部，嚴(yán)格受北東向構(gòu)造控制呈北東向分布；加里東中期花崗閃長(zhǎng)巖主要分布于工作區(qū)的東南部，呈巖基狀分布；印支期中粗粒花崗巖體分布于礦區(qū)的西部；整個(gè)礦區(qū)內(nèi)脈巖較發(fā)育，主要有輝綠巖，細(xì)粒閃長(zhǎng)巖等。

斷裂為北東向和北西向兩組，早期北東向張性斷裂控制了區(qū)內(nèi)巖體及地層的空間分布，北西向斷裂次之，以壓性近直立斷裂為主。

鉬礦化富集與構(gòu)造和巖漿活動(dòng)有關(guān)。其中北東向張性斷裂的次級(jí)斷裂及青龍村群變質(zhì)巖與燕山期似斑狀花崗巖接觸部位均為鉬元素相對(duì)富集的部位。礦區(qū)圍巖蝕變?yōu)橹?低溫?zé)嵋何g變，主要有硅化、絹云母化、綠泥石化、蒙脫石化、高嶺石化，其次可見(jiàn)黃鐵礦化、磁黃鐵礦化、硅化、螢石化、碳酸鹽化等，各種蝕變相互疊加，分帶不明顯。

土壤地球化學(xué)測(cè)量結(jié)果顯示礦區(qū)內(nèi)具有W、Mo、Sn、Bi異常，礦區(qū)南部礦化已經(jīng)圈閉，但是向北部仍然存在著較高的土壤地球化學(xué)異常，且有繼續(xù)向北延伸的趨勢(shì)。

激電探測(cè)結(jié)果顯示，礦區(qū)內(nèi)見(jiàn)有北北東向激電異常帶，橫穿整個(gè)礦區(qū)，具有南部寬北部窄的特點(diǎn)，向北趨于尖滅。

礦床類(lèi)型屬于典型的斑巖型鉬礦。

2 高光譜地質(zhì)應(yīng)用基本原理

2.1 ASD礦物光譜儀工作參數(shù)

美國(guó)ASD公司開(kāi)發(fā)的TerraSpec系列便攜式礦物光譜儀使用512陣元陣列PDA探測(cè)器和兩個(gè)獨(dú)立的InGaAs探測(cè)器，可探測(cè)波長(zhǎng)范圍為350～2500nm，覆蓋了全部可見(jiàn)光、近紅外、短波紅外波長(zhǎng)區(qū)間；靈敏度高，光譜分辨率強(qiáng)，光譜采樣間隔小，掃描時(shí)間短，每秒可得到10個(gè)光譜，可實(shí)時(shí)測(cè)量并觀察光譜曲線，能實(shí)現(xiàn)蝕變礦物的野外現(xiàn)場(chǎng)快速識(shí)別。

2.2 ASD礦物光譜儀工作原理

利用探測(cè)波長(zhǎng)區(qū)間光線照射被測(cè)樣品，并接收礦物的特征吸收光譜，通過(guò)對(duì)所獲的光譜曲線與標(biāo)準(zhǔn)曲線比對(duì)，反演被測(cè)樣品蝕變礦物信息；ASD礦物光譜儀是利用短波紅外光譜測(cè)量技術(shù)，對(duì)肉眼難以識(shí)別的層狀硅酸鹽、粘土、碳酸鹽以及部分硫酸鹽礦物進(jìn)行快速識(shí)別的實(shí)用儀器。通過(guò)不同地質(zhì)環(huán)境的光譜參考數(shù)據(jù)庫(kù)及計(jì)算機(jī)數(shù)據(jù)處理軟件，輔之野外觀察和必要的巖石學(xué)分析，可以從ASD礦物光譜儀測(cè)量獲得的譜線中得出研究區(qū)蝕變礦物組合的信息。

2.3 高光譜礦物識(shí)別基本原理

礦物光譜特征是礦物識(shí)別提取的前提，在400～2500nm光譜覆蓋范圍內(nèi)，礦物內(nèi)部的離子、分子、基團(tuán)及其它組合的光譜特征基本都具有其唯一性，礦物光譜吸收特征的變化主要是與其物質(zhì)組分以及晶格結(jié)構(gòu)有關(guān)的電子與晶體場(chǎng)的相互作用影響的。相對(duì)于傳統(tǒng)光譜測(cè)量，高光譜測(cè)量精度更好，能夠獲得更為精確的連續(xù)光譜曲線(鄧美容等，2010)。

圍巖蝕變常與礦體伴生且其分布范圍一般比礦體分布范圍廣，因而是一種重要的找礦標(biāo)志。蝕變礦物中的陽(yáng)離子Fe2+與Fe3+和陰離子OH-與CO32-在350～2500nm波段區(qū)間具有明顯的特征吸收譜帶(劉燕君等，1993；孫雨等；2015)，因此其可以被ASD礦物光譜儀所識(shí)別。

3 光譜數(shù)據(jù)采集

采集樣品的光譜數(shù)據(jù)時(shí)，使用了ASD礦物光譜儀自帶的FieldSpec Pro RS3程序。由于樣品的光譜測(cè)試結(jié)果受多種因素影響，例如：光照條件(太陽(yáng)高度角、太陽(yáng)方位角)、大氣特性(風(fēng)、云)、樣品的物理性質(zhì)及幾何形狀、掃描速度等。因此為了得到測(cè)試樣品的反射率，需要采用基于統(tǒng)計(jì)學(xué)模型的平場(chǎng)域法進(jìn)行光譜重建工作獲得反射率數(shù)據(jù)，首先應(yīng)對(duì)標(biāo)準(zhǔn)白板進(jìn)行反射光譜測(cè)量，令其為參考光譜即平場(chǎng)域，樣品的相對(duì)光譜反射率數(shù)據(jù)是樣品反射譜測(cè)量所得的高光譜數(shù)據(jù)與平場(chǎng)域平均光譜數(shù)據(jù)的比值。標(biāo)準(zhǔn)白板與樣品的反射光譜測(cè)試環(huán)境假定完全一致的情況下，可以假設(shè)標(biāo)準(zhǔn)白板對(duì)各波段的反射率均為一個(gè)恒定值1，那么樣品的相對(duì)光譜反射率數(shù)據(jù)可以認(rèn)為其是絕對(duì)反射率數(shù)據(jù)。

為了消除因?yàn)椴僮魃袭a(chǎn)生的偶然誤差，本次工作采用了對(duì)同一采樣點(diǎn)進(jìn)行多次數(shù)據(jù)采集的方法，確保最終得到的數(shù)據(jù)的準(zhǔn)確性。

本次ASD礦物光譜儀數(shù)據(jù)采集的對(duì)象為礦區(qū)新近施工的金剛石鉆孔巖芯，數(shù)據(jù)提取位置不包括每個(gè)鉆孔開(kāi)孔位置提取的第四系浮土、轉(zhuǎn)石等，數(shù)據(jù)采集點(diǎn)間隔為1米，本次共采集巖芯長(zhǎng)度18661米，樣品采集點(diǎn)18664個(gè)。

4 蝕變礦物識(shí)別與結(jié)果驗(yàn)證

對(duì)采集到的礦物光譜曲線，采用TSG(The Spectral Geologist)軟件進(jìn)行識(shí)別的，該軟件由澳大利亞聯(lián)邦科學(xué)與工業(yè)研究組織(CSIRO)研發(fā)，其中儲(chǔ)存了大量的特征蝕變礦物標(biāo)準(zhǔn)譜線。解譯時(shí)軟件自動(dòng)將采集的光譜曲線與數(shù)據(jù)庫(kù)中的標(biāo)準(zhǔn)譜線相對(duì)比，生成光譜質(zhì)量評(píng)價(jià)圖，并對(duì)礦物進(jìn)行自動(dòng)識(shí)別。

表1 蝕變礦物含量總表

4.1 主要蝕變礦物

經(jīng)過(guò)解譯，本次共識(shí)別出蝕變礦物38種(表1)，其中12種蝕變礦物單個(gè)礦物在樣品中的平均含量大于1.5%，為主要的蝕變礦物，按照含量的多少依次為：絹云母、蒙脫石、黑云母、鐵鎂綠泥石、金云母、伊利石化絹云母、鐵綠泥石、角閃石、高嶺石(結(jié)晶好)、綠簾石、坡縷石、高嶺石(結(jié)晶差)。

4.2 蝕變礦物光譜特征

上述12種蝕變礦物，應(yīng)挑選出具有代表性，且與鉬礦化蝕變具有一定相關(guān)性的蝕變礦物，才可以更有效的進(jìn)行找礦預(yù)測(cè)。角閃石和黑云母與巖漿作用關(guān)系密切，常見(jiàn)于各種中、酸性侵入巖中，且礦區(qū)及周邊地區(qū)均廣泛分布；坡縷石主要形成于火山沉積、內(nèi)陸湖相沉積和表生淋濾等，其與熱液礦床無(wú)關(guān)；綠簾石廣泛發(fā)育于巖漿巖、變質(zhì)巖和沉積巖熱液交代后的圍巖蝕變中，其在礦區(qū)周邊多與基性巖脈有關(guān)，與礦化關(guān)系不大；高嶺石(結(jié)晶好)、金云母和鐵綠泥石在樣品中分布不均勻，在部分鉆孔中含量極高，但在多數(shù)見(jiàn)礦鉆孔中均未見(jiàn)到，因此上述七種礦物并不具有代表性，可以排除。所以主要考慮其余5種蝕變礦物與鉬礦化的關(guān)系。對(duì)此5種蝕變礦物通過(guò)TSG軟件中自帶的標(biāo)準(zhǔn)蝕變礦物光譜特征進(jìn)行對(duì)比和解譯。其主要光譜特征如下：

①標(biāo)準(zhǔn)的絹云母礦物的Al-OH的吸收谷值波長(zhǎng)從2190nm到2225nm，樣品中的Al-OH的吸收谷值波長(zhǎng)是2201～2217nm，與標(biāo)準(zhǔn)數(shù)據(jù)庫(kù)中的光譜線吻合很好。短波吸收谷值接近2190nm(也就是特征吸收谷值越短)，表明蝕變礦物中鈉含量越高，而短波吸收峰谷越長(zhǎng)，接近2206nm，蝕變礦物中鉀含量越高，樣品的的吸收谷波長(zhǎng)介于2201～2217nm，主要位于2205nm，靠近2206nm(圖1a)，礦區(qū)絹云母是相對(duì)富鉀的。

②綠泥石光譜Fe-OH的吸收谷波長(zhǎng)為2250～2260nm。在綠泥石成分中富鐵到富鎂含量的變化，會(huì)影響光譜中Fe-OH吸收谷的位置，一般隨著Fe和Mg的二價(jià)離子在綠泥石晶體中八面體含量的變化而變化，富鐵綠泥石的吸收谷一般在2260nm，而富鎂綠泥石的吸收谷一般在2250nm,隨著鐵、鎂含量的不同，2250～2260nm處的峰值變化呈現(xiàn)出一種線狀的變化趨勢(shì)(Scottetal.,1998)。在斑巖型或熱液礦床中，富鎂的綠泥石往往更靠近礦化帶。樣品的吸收谷為介于2248～2256nm之間，主要為2253nm附近(圖1b)，解譯為鐵鎂綠泥石。

③對(duì)于短波紅外光蒙脫石在2208nm波長(zhǎng)顯示強(qiáng)烈的吸收，同時(shí)在1410nm和1910nm也有明顯的吸收，但不對(duì)稱，這主要是因?yàn)閹r石中的水分吸收了部分的紅外光，與絹云母相比，蒙脫石在2345nm和2440nm沒(méi)有明顯的吸收，本次樣品吸收谷介于1407～1418nm、1900～1918nm和2200～2219nm之間，主要為1414nm、1909nm和2208nm附近(圖1c)，解譯為蒙脫石。

④ 絹云母中Na和Al含量較多時(shí)，光譜曲線吸收谷在2190nm附近，當(dāng)Al和K含量較高時(shí)，吸收谷在2200nm附近，當(dāng)貧Al且富K時(shí)，吸收谷在2210～2225nm之間；伊利石屬于含水層狀硅酸鹽礦物，其含有Al-OH和H2O兩種含氫基團(tuán)，Al-OH基團(tuán)在2200nm附近有明顯的特征吸收谷，H20在1900 nm附近有明顯的吸收谷。樣品在1891～1911nm和2200～2211nm附近有明顯的吸收谷，多數(shù)位于1902nm和2205nm附近(圖1d)，其與伊利石吸收特征相似，解譯為伊利石化絹云母。

圖1標(biāo)準(zhǔn)礦物與樣品反射光譜對(duì)比圖Fig.1 Comparison of reflection spectra between standard minerals and samples

⑤高嶺石主要以1400nm和2200nm附近的吸收雙峰為特征。第一個(gè)吸收峰值在1400～1450nm，第二個(gè)吸收峰值在2166～2206nm。本次光譜采集與現(xiàn)場(chǎng)礦物識(shí)別對(duì)比顯示，該礦區(qū)光譜曲線雙峰平緩程度與高嶺石的結(jié)晶程度有關(guān)，結(jié)晶差的高嶺石，吸收波峰的短波部分相對(duì)平緩，而結(jié)晶好的高嶺石，吸收波峰的短波部分相對(duì)陡峭，本次解譯過(guò)程中高嶺石少量樣品的雙峰中短波部分相對(duì)陡峭(圖1e)，大多雙峰中短波部分相對(duì)平緩(圖1f)。

4.3 識(shí)別結(jié)果驗(yàn)證

基于數(shù)據(jù)可靠性，本次挑選巖芯完整且蝕變礦物含量較多的巖芯段進(jìn)行巖石標(biāo)本采集，共采集了5件巖石標(biāo)本進(jìn)行巖礦鑒定，來(lái)對(duì)比檢驗(yàn)ASD礦物光譜儀測(cè)試數(shù)據(jù)是否準(zhǔn)確(表2)。從表中可見(jiàn)采用ASD礦物光譜儀采集樣品的光譜數(shù)據(jù)，在經(jīng)過(guò)TSG軟件對(duì)光譜進(jìn)行解譯后，得出的蝕變礦物信息是準(zhǔn)確的。

表2 巖礦鑒定與ASD礦物光譜儀蝕變礦物對(duì)比

5 蝕變礦物分布特征及與鉬礦化的關(guān)系

5.1 蝕變礦物與鉬礦化的相關(guān)性

首先查明蝕變礦物和鉬礦化的關(guān)系，隨機(jī)統(tǒng)計(jì)鉆孔中蝕變礦物含量與鉬礦化品位的相關(guān)性系數(shù)。從統(tǒng)計(jì)結(jié)果來(lái)看(表3)，絹云母、蒙脫石、伊利石化絹云母、高嶺石(結(jié)晶差)和鉬礦化表現(xiàn)出一定程度的正相關(guān)，鐵鎂綠泥石與鉬礦化表現(xiàn)出一定程度的負(fù)相關(guān)。

5.2 蝕變礦物含量下限的確定

為了更直觀的顯示各蝕變礦物分布和鉬礦化(體)空間關(guān)系，在使用Micromine軟件建立三維模型過(guò)程中，需要確定蝕變礦物含量的下限，通過(guò)眾值法統(tǒng)計(jì)蝕變礦物的含量符合正態(tài)分布特征(圖2)，眾值和中位值都在20%左右，因此蝕變帶圈定的原則是：以蝕變礦物含量20%作為邊界值圈連蝕變體，但是如果蝕變帶在空間連續(xù)性較好時(shí)，可根據(jù)實(shí)際數(shù)據(jù)進(jìn)行適當(dāng)調(diào)整。

圖2 蝕變礦物含量直方圖Fig.2 Histogram of content of alteration minerals

5.3 蝕變礦物與鉬礦化空間分布特征

礦體與蝕變礦物空間分布形態(tài)及位置(圖3)對(duì)比結(jié)果主要為：絹云母化蝕變體規(guī)模較大(圖3b)，基本遍布整個(gè)礦區(qū)，其包括了礦區(qū)內(nèi)的所有鉬礦化(體)，絹云母化發(fā)育較強(qiáng)的部位鉬礦化也較強(qiáng)；礦區(qū)蒙脫石化較強(qiáng)，蝕變體分布范圍較廣，空間上基本上覆蓋所有鉬礦化(體)，其蝕變體內(nèi)部和兩側(cè)分布有較強(qiáng)的鉬礦化，蒙脫石化較強(qiáng)的部位也是鉬礦化較強(qiáng)的部位(圖4c)；伊利石化絹云母的蝕變相對(duì)弱，僅見(jiàn)細(xì)小條帶狀蝕變體零星散布于空間各位置中，在該蝕變體周?chē)３霈F(xiàn)強(qiáng)鉬礦化，而蝕變體內(nèi)部鉬礦化反而相對(duì)較弱(圖4d)；高嶺石(結(jié)晶差)蝕變中等，在空間上主要為細(xì)條帶狀零散斷續(xù)分布，蝕變體空間分布特征與鉬礦體的空間分布特征相似，高嶺石(結(jié)晶差)化相對(duì)發(fā)育的地段，鉬礦化也相對(duì)較發(fā)育(圖4e)；鐵鎂綠泥石化蝕變體分布在鉬礦化帶的上下兩側(cè)，與鉬礦體之間有一定的距離，在鐵鎂綠泥石化出現(xiàn)的地方，一般不出現(xiàn)鉬礦化(圖4f)；總體而言，上述五種蝕變?cè)诳臻g上均與鉬礦化關(guān)系密切。

5.4 蝕變礦物分帶特征分析

按照礦化蝕變與礦體的遠(yuǎn)近可以把礦化蝕變劃分為內(nèi)、中、外三個(gè)蝕變帶。內(nèi)帶的蝕變礦物組合是高嶺石(結(jié)晶差)、蒙脫石和絹云母，絹云母化蝕變最強(qiáng)，蒙脫石化蝕變次之，高嶺石(結(jié)晶差)化蝕變相比較弱，鉬礦化常常分布在三種蝕變體之中，蝕變相對(duì)較強(qiáng)的區(qū)域，鉬礦化也相對(duì)較強(qiáng)；中帶的蝕變礦物組合是伊利石化絹云母、蒙脫石和絹云母，伊利石化絹云母在此三種蝕變中最弱，鉬礦化的周?chē)秶鷥?nèi)均能發(fā)現(xiàn)這三種蝕變，且這三種蝕變相對(duì)越好，其周?chē)f礦化也相對(duì)較好；外帶的蝕變礦物組合是鐵鎂綠泥石和絹云母，鐵鎂綠泥石化出現(xiàn)在鉬礦化帶的兩側(cè)，與鉬礦化帶有一定的間距，絹云母在整個(gè)礦區(qū)均有分布，發(fā)現(xiàn)鐵鎂綠泥石和絹云母蝕變，那么在一定范圍內(nèi)應(yīng)會(huì)有一定程度的鉬礦化。

表3 樣品的蝕變礦物與鉬礦化的相關(guān)系數(shù)

圖3 蝕變礦物與鉬礦化空間分布特征Fig.3 Spatial distribution characteristics of molybdenum mineralization and alteration minerals a-礦體；b-絹云母；c-蒙脫石；d-伊利石化絹云母；e-高嶺土(結(jié)晶差)；f-鐵鎂綠泥石a-orebody;b-sericite;c-montmorillonite;d-illitized sericite;e-kaolinite(poorly crystallized);f-brunsvigite

通過(guò)上述分析，內(nèi)、中、外三個(gè)蝕變帶可以分別表述為絹云巖化-弱粘土化帶、粘土化-絹云巖化、絹云巖化-綠泥石化(青磐巖化)帶。其與大多數(shù)貧氟的斑巖鉬礦蝕變分帶由內(nèi)向外常常具有鉀化核、石英-絹云母-黃鐵礦、粘土化帶和青磐巖化帶等特征(Westraetal., 1981)是類(lèi)似。

三維模型數(shù)據(jù)顯示，上述五種蝕變體空間上在礦區(qū)北部均未封閉，而是繼續(xù)開(kāi)放，表明礦區(qū)北部應(yīng)還存在一定程度的礦化異常，這與現(xiàn)有化探數(shù)據(jù)中礦區(qū)北部存在較高的化探異常是一致的。

6 結(jié)論

(1)礦區(qū)礦化蝕變可分為內(nèi)、中、外三個(gè)蝕變帶，分別為絹云巖化-弱粘土化(高嶺石(結(jié)晶差)和蒙脫石)帶、粘土化(伊利石化絹云母和蒙脫石)-絹云巖化帶、絹云巖化-綠泥石化帶。

(2)蝕變礦物空間分布和分帶特征顯示鉬礦化向北部具有繼續(xù)延伸的趨勢(shì)，其與礦區(qū)土壤地球化學(xué)異常指示方向是一致的。

(3)通過(guò)對(duì)礦區(qū)已經(jīng)施工鉆孔巖芯中蝕變礦物的提取和識(shí)別，分析蝕變礦物在空間上的分布和分帶特征，更為清晰的展現(xiàn)蝕變礦物在空間上與礦化之間的關(guān)系，可以快速有效的進(jìn)行找礦預(yù)測(cè)，尤其在大型礦山的預(yù)測(cè)中可以有效的利用已有工程進(jìn)行深部和外圍找礦預(yù)測(cè)。

[注釋]

① Scott K M， Yang K， Huntington J F. 1998. The application of spectral re-flectauce studies to chlorites in exploration [R].Australia： CSIRO，Exploration & Mining Report 545R

Deng Mei-rong，Zhao Yun-sheng,Liu Lan-fang.2010.The Research of Hyper-spectral Polarized Characteristics of Several Typical Minerals[J].Journal of Hengyang Normal University,31(6)：97-101(in Chinese with English abstract)

Du Xiao-di,Tang Yue, Liu De-chang,Li Zhao,Wang Hai-da.2015.The application of aerial hyperspectral technology to the geological survey for oil and gas in the East Junggar Basin[J].Geology in China，42(1)：275-287(in Chinese with English abstract)

Gan Fu-ping,Wang Run-sheng，Yang Su-ming.2002.Studying on the alteration minerals identification using hyperion data[J].Remote Sensing for Land &Resources,(4):44-51(in Chinese with English abstract)

Gan Fu-ping,Wang Run-sheng，Ma Ai-nai.2003.Spectral identification tree(sit) for mineral extraction based on spectral characteristics of minerals[J].Earth Science Frontiers,10(2):445-454(in Chinese with English abstract)

Guo Na.2012.The altered mineral assemblage research in Jiama porphyryskarn copper deposit-based on Hyperspectral high frequency wave infrared technology[D].Cheng Du:Chengdu University of Technology:1-122(in Chinese with English abstract)

Hou Yi,Wang Mao-zhi,Zhang Xiao-dong.2014.Application of SASI Hyperspectral Data in Oil and Gas Exploration[J].Journal of Chengdu University of Technology(Science & Technology Edition),33(4)：347-350(in Chinese with English abstract)

Lei Tian-ci，Zhu Ming-qiang，Zhou Wan-peng.2005. Application on data mining of hyperspectrum to Identification and Extraction of alteration minerals[J].Resour,ces Environment & Engineering,19(3)：213-218(in Chinese with English abstract)

Li Nang, Xiao Ke-yang, Chen Xi-qiu, Lou De-bo. 2010.Mineralizing alteration extraction based on Hyperion hyper-spectral data-taking Tanyaokou deposit, Inner Mongolia, China as as example. Geological Bulletin of China, 29(10):1558-1563(in Chinese with English abstract)

Li Zhi-zhong Z, Yang Ri-hong, Dang Fu-xing, Zhang Xian-feng, Tan Bbing-xiang, Zhao Hui-jie. 2009.The hyperspectral remote sensing technology and its application[J]. Geological Bulletin of China, 28(2/3) ;270-277(in Chinese with English abstract)

Liu Yan-jun，Jin Li-fang.1993. Remote sensing ground model for mineral information[M].Beijing:Geological Publishing House:1-157(in Chinese)

Sun Wei-dong,Chen Jian-ming,Wang Run-shen,Yan Ba-Kun,Yu Hao.2010.The application and research of hyperion for hyperspectral mineral mapping on east of TianShan mountain[J].Xin Jiang Geology,28(2)：214-217(in Chinese with English abstract)

Sun Yu,Nie Jiang-tao,Tian Feng, Qin Kai,Yang Guo-fang,Wang Jian.2015.Alteration mineral mapping of the Xiangshan uranium core using HySpex imaging hyperspectral data and its geological significance[J].Geology and Exporation, 51(1) :0165-0174(in Chinese with English abstract)

Tian Feng,Dong Li-na, Yang Su-ming, Wang Run-sheng .Application of combined spectra of mixed minerals to mapping altered minerals A case study in theYunnan region based on Hyperion data[J]. Geology and Exploration 2010,46 (2): 0331-0337(in Chinese with English abstract)

Wang Cheng-hui，Song Quan-heng，Wang Deng-hong，Li Li-xing，YU Cheng，Wang Zhi-gang，Qu Wen-jun，Du An-dao，Ying Li-juan.2009.Re-Os Isotopic Dating of Molybdenite from the Daheishan Molybdenum Deposit of Jilin Province and Its Geological Significance[J].Rock and Mineral analysis,28(3)269-273(in Chinese with English abstract)

Wang Jin-nian，Zhang Bing，Liu Jian-gui，Tong Qing-xi and Zheng Lan-fen.1999.Hyperspectral Data Mining-Toward Target Recognition and Classification[J].Journal of Image and Graphics,4(11):957-964(in Chinese with English abstract)

Wang Ya-jun,Wang Qin-jun,Chen Yu,Hu Fang,Xu Ru,Lin-Qi-Zhong.2015.A Method to Extract Content of Minerals Based on Measured Hyperspectral Data[J].Spectroscopy and Spectral Analysis，35(6)：1700-1704(in Chinese with English abstract)

Westra G, Keith S B. 1981. Classification and genesis of stockwork molybdenum deposits[J].Economic Geology,76(4):844-873

Wu Feng-huang，Jiang Xiping，Jiang Yu.2013.Hyper-spectral Mineral Identification Based on Improved ICA Algorithm[J]. Journal of Henan Normal University(Natural Science)，2013，41(5)：48-52(in Chinese with English abstract)

Xu Ning,Hu Yu-xin,Lei Bin,Zhang Cong,Wang Da-ming,Shi Tao.2013.Automated mineral information extraction based on PPI algorithm for hyperspectral imagery[J].Science of Surveying and Mapping,38(4)：138-141(in Chinese with English abstract)

Yan Shou-xun，Zhang Bing，Zhao Yong-chao，Zheng Lan-fen,Tong Qing-xi，Yang Kai.2003.Summarizing the VIS-NIR Spectra of Minerals and Rocks[J].Remote Sensing Technology and Application,18(4)：191-201(in Chinese with English abstract)

Yang Zhi-ming,Hou Zeng-qian,Yang Zhu-sen,Qu Huan-chun,Li Zhen-qing,Liu Yun-fei.2012.Application of short wavelength infrared(SWIR) technique in exploration of poorly eroded porphyry Cu district: A case study of Niancun ore district, Tibet[J].Mineral Deposits,31(4):699-717(in Chinese with English abstract)

Zhang Ting，Ding Jian-li，WANG Fei.2010.Simulation of Image Multi-Spectrum Using Field Measured Endmember Spectrum[J].Spectroscopy and Spectral Analysis，30(11)：2889-2892(in Chinese with English abstract)

Zhou Qiang, Gan Fu-ping, Chen JianPing，Yan Bo-kun，Wang Run-sheng.Application of hyperspectral imaging data to mineral identification in Tudun Cu-Ni deposit, Hami，Xinjiang[J].Mineral Deposits，2006，25(3)：329-336(in Chinese with English abstract)

Zou Lin,Yang Zi-an,Zhu Gu-chang,Wu De-wen,Xiang Ai-qin.2006.New method of Alteration information extraction from Multi-Spectrum remote sensing data[J].Geology and Exploration,2006，42(6)：71-76(in Chinese with English abstract)

[附中文參考文獻(xiàn)]

鄧美容，趙云升，劉蘭芳.2010.幾種典型礦物的偏振高光譜特性研究[J].衡陽(yáng)師范學(xué)院學(xué)報(bào)，31(6)：97-101

杜小弟，唐 躍，劉德長(zhǎng)，李 昭，王海達(dá).2015.航空高光譜探測(cè)技術(shù)在準(zhǔn)噶爾盆地東部地區(qū)油氣調(diào)查中的應(yīng)用[J].中國(guó)地質(zhì)，42(1)：275-287

甘甫平，王潤(rùn)生，楊蘇明.2002.西藏Hyperion數(shù)據(jù)蝕變礦物識(shí)別初步研究[J].國(guó)土資源遙感，(4):44-51

甘甫平，王潤(rùn)生，馬藹乃.2003.基于特征譜帶的高光譜遙感礦物譜系識(shí)別[J].地學(xué)前緣，10(2):445-454

郭 娜.2012.甲瑪斑巖一矽卡巖型銅礦床蝕變礦物組合研究-基于高光譜短波紅外技術(shù)[D].成都：成都理工大學(xué)：1-122

侯 毅，王茂芝，張小東.2014.SASI高光譜數(shù)據(jù)在油氣勘探中的應(yīng)用[J].成都大學(xué)學(xué)報(bào)(自然科學(xué)版)，33(4)：347-350

雷天賜，祝明強(qiáng)，周萬(wàn)篷.2005.高光譜數(shù)據(jù)挖掘在蝕變礦物識(shí)別與提取中的應(yīng)用[J].資源環(huán)境與工程，19(3)：213-218

李 楠，肖克炎，陳析璆，婁德波.2010.基于Hyperion高光譜數(shù)據(jù)的礦物蝕變提取—以內(nèi)蒙古西部狼山地區(qū)炭窯口礦床為例[J].地質(zhì)通報(bào)，2010，29(10)：1558-1563

李志忠，楊日紅，黨福星，張顯峰，譚炳香，趙慧潔.2009.高光譜遙感衛(wèi)星技術(shù)及其地質(zhì)應(yīng)用[J].地質(zhì)通報(bào)，28(2-3)：270-277

劉燕君，金麗芳.1993.礦產(chǎn)信息的遙感地面模式[M].北京：地質(zhì)出版社：1-157

孫衛(wèi)東，陳建明，王潤(rùn)生，閻柏昆，于 浩.2010.阿爾金地區(qū)高光譜遙感礦物填圖方法及應(yīng)用研究[J].新疆地質(zhì)，28(2)：214-217

孫 雨，聶江濤，田 豐，秦 凱，楊國(guó)防，王 健.2015.相山鈾礦巖芯HySpex成像高光譜數(shù)據(jù)蝕變礦物提取及其地質(zhì)意義[J].地質(zhì)與勘探，51(1)：165-174

田 豐，董麗娜，楊蘇明，王潤(rùn)生.2010.混合礦物組合光譜在蝕變礦物填圖中的應(yīng)用一以云南香格里拉地區(qū)Hyperion數(shù)據(jù)蝕變礦物填圖為例[J].地質(zhì)與勘探，46(2)：331-337

王成輝，松權(quán)衡，王登紅，李立興，于 城，汪志剛，屈文俊，杜安道，應(yīng)立娟.2009.吉林大黑山超大型鉑礦輝鉑礦錸-鋨同位素定年及其地質(zhì)意義[J].巖礦測(cè)試，28(3)269-273

王晉年，張 兵，劉建貴，童慶禧，鄭蘭芬.1999.以地物識(shí)別和分類(lèi)為目標(biāo)的高光譜數(shù)據(jù)挖掘[J].中國(guó)圖像圖形學(xué)報(bào)，4(11):957-964

王亞軍，王欽軍，陳 玉，胡 芳，徐 茹，藺啟忠.2015.基于實(shí)測(cè)高光譜數(shù)據(jù)的礦物含量提取方法研究[J].光譜學(xué)與光譜分析，35(6)：1700-1704

吳鳳凰，蔣夕平，蔣 昱.2013.基于優(yōu)化ICA算法的高光譜礦物識(shí)別[J].河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，41(5)：48-52

許 寧，胡玉新，雷 斌，張 聰，汪大明，石 濤.2013.一種基于PPI的高光譜數(shù)據(jù)礦物信息白動(dòng)提取方法[J].測(cè)繪科學(xué)，38(4)：138-141

燕守勛，張 兵，趙永超，鄭蘭芬、童慶禧，楊 凱.2003.礦物與巖石的可見(jiàn)-近紅外光譜特性綜述.遙感技術(shù)與應(yīng)用，18(4)：191-201

楊志明，侯增謙，楊竹森，曲煥春，李振清，劉云飛.2012.短波紅外光譜技術(shù)在淺剝蝕斑巖銅礦區(qū)勘查中的應(yīng)用—以西藏念村礦區(qū)為例[J].礦床地質(zhì)，31(4):699-717

張 婷，丁建麗，王 飛.2010.基于實(shí)測(cè)端元光譜的多光譜圖像光譜模擬研究[J].光譜與光譜分析，30(11)：2889-2892

周 強(qiáng)，甘甫平，陳建平，閆柏琨，王潤(rùn)生.2006.哈密土墩礦區(qū)高光譜影像蝕變礦物識(shí)別初步研究[J].礦床地質(zhì)，25(3)：329-336

鄒 林，楊自安，朱谷昌，吳德文，相愛(ài)芹.2006.多光譜遙感蝕變信息提取新方法研究[J].地質(zhì)與勘探，42(6)：71-76

Application of Alteration Mineral Extraction base on Hyperspectral Data to Geological Prospecting of a Molybdenum Deposit in Jilin Province

LI Xi1,2，PAN Zhen-xing2，WANG Qi2，GAO Ya-ning1，HAN Ke1，LIU Zhang-le2

(1.SchoolofEarthScienceandResources,Chang’anUniversity,Xi’an,Shaanxi710054; 2.NorthwestNonferrousGeologicalandMiningGroupCo.,Ltd.,Xi’an,Shaanxi710054)

The hyperspectra method for extracting and recognizing rocks and minerals has been commonly used in the world.By means of this technology, this study builds on ground spectrum test, data processing and analysis of spectrum characteristics to extract and recognize the altered minerals in a porphyry molybdenum ore district of Jilin Prorince.Then, by analyzing the relevance of Mo mineralization with alteration,the following five representative altered minerals of this area have been selected:sericite,montmorillonite,illite-sericite,kaolinite(poorly crystallized) and brunsvigite. In addition,in combination with the characteristics of spatial distribution and zonation of those five altered minerals above,this work analyzed their relation to Mo-mineralization,and suggested that the altered belt can be further divided into inner, middle and outer zones.The results permit to predict the direction for ore-search and to provide a new line of thought for mine exploration.

hyperspectral,spectrum characteristics,porphyry molybdenum deposit,alteration minerals

2015-12-30；

2016-03-20；[責(zé)任編輯]陳英富。

陜西省自然科學(xué)基金(2016JM4001) 、中央高?；究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金(310827151056，310827153408)資助。

李 希(1986年-)，男，在讀博士生，工程師，從事金屬礦產(chǎn)勘查和礦山數(shù)字化工作，E-mail:lixikey@163.com

P579 [文獻(xiàn)標(biāo)示碼]A

0495-5331(2016)03-0489-08