胡紫豪，唐菊興，張廷斌，吳 華，徐志忠，別小娟

西藏多不雜斑巖銅礦ASTER遙感蝕變異常特征

胡紫豪1，唐菊興2，張廷斌1，吳 華3，徐志忠4，別小娟5

(1.成都理工大學(xué)地球科學(xué)學(xué)院，成都 610059;2.中國地質(zhì)科學(xué)院礦產(chǎn)資源研究所，北京 100037;3.西藏自治區(qū)地質(zhì)調(diào)查院，拉薩 850000;4.西藏地勘局第五地質(zhì)大隊(duì)，格爾木 816000;5.成都理工大學(xué)旅游與城鄉(xiāng)規(guī)劃學(xué)院，成都 610059)

采用ASTER數(shù)據(jù)，利用“掩模+定向主成分分析”方法提取了多不雜斑巖銅礦的遙感蝕變異常。對(duì)蝕變異常的特征、蝕變異常與地面蝕變對(duì)應(yīng)關(guān)系、蝕變異常與已知礦體的關(guān)系等幾個(gè)方面進(jìn)行了初步探討，認(rèn)為利用ASTER數(shù)據(jù)提取的蝕變異常效果較好，為西藏自治區(qū)礦產(chǎn)資源潛力評(píng)價(jià)多龍銅礦預(yù)測(cè)工作區(qū)的最小預(yù)測(cè)區(qū)圈定提供了遙感依據(jù)。

遙感;蝕變異常;ASTER;多不雜斑巖銅礦

0 引言

近礦圍巖蝕變是成礦物質(zhì)在逐步富集成礦過程中留下的印跡。絕大多數(shù)熱液礦床都伴隨有圍巖的交代蝕變現(xiàn)象，而且蝕變帶范圍常常大于礦體分布范圍的數(shù)倍至幾十倍，因此圍巖蝕變是重要的找礦標(biāo)志之一。在全國礦產(chǎn)資源潛力評(píng)價(jià)中，利用遙感技術(shù)提取圍巖蝕變信息(遙感蝕變異常)是一項(xiàng)重要的工作內(nèi)容，采用的基礎(chǔ)數(shù)據(jù)源多為TM/ETM+。

ASTER是由美國和日本聯(lián)合開發(fā)的高性能衛(wèi)星傳感器，從0.52～11.6 μm共設(shè)置有14個(gè)波段，與TM/ETM+相比，光譜范圍覆蓋更寬，輻射分辨率更高。在西藏自治區(qū)礦產(chǎn)資源潛力評(píng)價(jià)典型礦床的研究中，遙感專題組利用ASTER數(shù)據(jù)開展了礦床圍巖蝕變的異常信息提取和研究，取得了較好效果。

1 研究區(qū)地質(zhì)概況

多不雜斑巖型銅金礦是繼西藏驅(qū)龍、玉龍之后發(fā)現(xiàn)的又一大型斑巖型銅金礦床，位于多龍銅礦預(yù)測(cè)工作區(qū)內(nèi)，大地構(gòu)造位置處于班公湖—怒江結(jié)合帶北側(cè)、羌塘地塊南緣的鐵格隆構(gòu)造巖漿帶中。出露的地層比較簡(jiǎn)單，主要是中侏羅統(tǒng)雁石坪群、下白堊統(tǒng)美日切組和新近系康托組。中侏羅統(tǒng)雁石坪群粉砂質(zhì)板巖、長(zhǎng)石石英砂巖和安山巖、英安巖、玄武質(zhì)安山巖是礦床的最重要圍巖［1－2］?；?、中 －酸性淺成侵入巖侵入到中侏羅統(tǒng)雁石坪群和下白堊統(tǒng)美日切組地層中，主要有輝長(zhǎng)巖、輝綠玢巖、閃長(zhǎng)巖、石英閃長(zhǎng)巖、石英閃長(zhǎng)玢巖和花崗閃長(zhǎng)斑巖等，以巖株、巖脈和巖墻等形態(tài)產(chǎn)出。其中，出現(xiàn)于多不雜礦區(qū)的石英閃長(zhǎng)玢巖和花崗斑巖主要呈巖株產(chǎn)出，蝕變強(qiáng)烈，是本礦區(qū)的主要含礦斑巖體［3－4］。銅礦體總體特征是規(guī)模大(含礦斑巖全巖礦化)、Cu品位較高(0.94%)、富 Au(0.21 g/t)［3］。礦石類型為典型的斑巖型礦石，礦石構(gòu)造以細(xì)脈浸染狀構(gòu)造為主，其次為角礫狀、脈狀構(gòu)造，礦石礦物主要為黃銅礦、黃鐵礦、磁鐵礦、斑銅礦和輝鉬礦［1］。

多不雜斑巖銅礦含礦斑巖和圍巖有強(qiáng)烈的熱液蝕變。蝕變類型有鉀化、硅化、絹英巖化(絹云母化+綠泥石化+硅化)、角巖化、青盤巖化、黃鐵礦化、磁鐵礦化和粘土化等，以絹英巖化、黃鐵礦化、鉀化和硅化最為普遍且蝕變強(qiáng)度大，其次為青磐巖化和高嶺土化［1］。

2 遙感蝕變異常提取

2.1 數(shù)據(jù)源及提取方法

2.1.1 數(shù)據(jù)源

本次研究采用ASTER的AST00226PRDAT0124和AST00160PRDAT024兩景圖像數(shù)據(jù)分別提取遙感羥基異常和遙感鐵染異常信息，數(shù)據(jù)獲取日期分別為2009－01－02和2006－11－25，產(chǎn)品級(jí)別為L(zhǎng)1B級(jí)，已經(jīng)過輻射校正和幾何糾正。ASTER數(shù)據(jù)在光譜覆蓋范圍和空間分辨率等方面均優(yōu)于TM/ETM+，可以提供15 m(可見光—近紅外)、30 m(短波紅外)以及90 m(熱紅外)3種空間分辨率數(shù)據(jù)(表1)。

表1 ASTER和TM/ETM+波段設(shè)置對(duì)比Tab.1 Comparison of the bands between ASTER and TM/ETM+

研究區(qū)的ASTER圖像無云，植被稀疏，能滿足蝕變信息提取的要求(圖1)。

圖1 研究區(qū)ASTER7(R)3(G)1(B)假彩色合成圖像Fig.1 Pseudo － color composite image of ASTER7(R)3(G)1(B)in study area

2.1.2 遙感蝕變異常提取方法

根據(jù)研究區(qū)蝕變礦物組合并結(jié)合ASTER傳感器光譜響應(yīng)特征，在可見光—短波紅外光譜區(qū)間，主要識(shí)別羥基礦物和含F(xiàn)e3+礦物。圖2為研究區(qū)主要蝕變礦物的USGS(United States Geological Survey，美國地質(zhì)調(diào)查局)實(shí)驗(yàn)室光譜曲線。

圖2 研究區(qū)主要蝕變礦物USGS波譜曲線Fig.2 USGS spectral plots of the main alteration minerals in study area

由圖2可以看出，羥基礦物在短波紅外(SWIR)光譜區(qū)間產(chǎn)生吸收特征譜帶，高嶺土、伊利石等含Al羥基標(biāo)志吸收谷在2.2 μm 附近，對(duì)應(yīng)ASTER5，6波段;綠泥石等含Mg羥基在2.3 μm附近產(chǎn)生標(biāo)志吸收谷，對(duì)應(yīng)ASTER8波段;由此可見，羥基礦物特征譜帶的精確位置取決于與OH－相連的金屬離子［5］。鐵在氧化環(huán)境下多以含F(xiàn)e3+礦物形式存在，如褐鐵礦、赤鐵礦、針鐵礦及黃鉀鐵礬等，其分別在0.5 μm 和 0.9 μm 處出現(xiàn)吸收特征，對(duì)應(yīng) ASTER1，3波段。

對(duì)遙感蝕變異常的提取方法，很多學(xué)者都進(jìn)行過相關(guān)研究［6－10］。本次工作借鑒全國礦產(chǎn)資源潛力評(píng)價(jià)遙感異常提取技術(shù)流程，基于ENVI 4.4遙感影像處理系統(tǒng)，采用“掩模+定向主成分分析”的方法，以“景”為單位開展了研究區(qū)遙感羥基與鐵染異常信息的提取(其中主要對(duì)河流、第四系等干擾因素進(jìn)行了“掩?！碧幚?，以排除非礦化信息的干擾)。

2.2 羥基異常提取

ASTER數(shù)據(jù)在短波紅外波段擁有更高的光譜分辨率，這使得識(shí)別不同類型的含羥基蝕變礦物成為可能。在利用ASTER不同波段組合的數(shù)據(jù)進(jìn)行主成分分析、提取不同的含羥基蝕變礦物方面已有過不少的研究，例如毛曉長(zhǎng)等［11］和楊佳佳等［12］等利用ASTER1，4，6，7波段提取高嶺石化、絹云母化等蝕變信息;叢麗娟等［5］和耿新霞等［13］等利用ASTER1，3，4，8波段提取黑云母化、綠泥石化等蝕變信息;呂鳳軍等［14］和黃照強(qiáng)等［15］等利用ASTER 1，3，4，5(6)波段提取方解石、白云石等碳酸鹽化的蝕變信息。本文則分別采用ASTER1，3，4，5(6)波段，ASTER1，3，4，8 波段和 ASTER1，4，6，7 波段3 種波段組合的主成分分析方案提取了研究區(qū)含羥基蝕變信息。通過與已知地面蝕變信息的對(duì)比認(rèn)為，采用ASTER1，3，4，8波段主成分分析提取的蝕變信息效果較好。

采用ASTER1，3，4，8波段提取羥基異常信息(主要是綠泥石化、方解石化、黑云母化的蝕變信息)。對(duì)進(jìn)行掩模運(yùn)算后的圖像選取ASTER1，3，4，8波段進(jìn)行定向主成分分析，其特征向量見表2。

表2 ASTER1，3，4，8波段主成分分析特征向量Tab.2 PCA eigenvector statistics for ASTER1，3，4 and 8

在第三主成分(PC3)中，ASTER4波段對(duì)應(yīng)的特征向量值最大(－0.977)，反映了此類羥基蝕變礦物在該波段的高反射特性。但由于該值為負(fù)值，應(yīng)先對(duì)PC3進(jìn)行取反處理，取反后得到的圖像可作為羥基異常圖像;再對(duì)羥基異常圖像進(jìn)行低通濾波和分級(jí)處理(即首先根據(jù)“平均值+n倍標(biāo)準(zhǔn)離差”的值確定出異常下限;最后對(duì)異常值進(jìn)行平均分配，確定出異常的一、二、三級(jí))。分級(jí)的結(jié)果為:164～194為羥基三級(jí)異常，195～225為羥基二級(jí)異常，226～255為羥基一級(jí)異常。

2.3 鐵染異常提取

對(duì)進(jìn)行掩模處理后的圖像選取ASTER1，2，3，4波段進(jìn)行定向主成分分析，提取鐵染異常信息，其特征向量見表3。

表3 ASTER1，2，3，4波段主成分分析特征向量Tab.3 PCA eigenvector statistics for ASTER 1，2，3 and 4

根據(jù)含F(xiàn)e3+礦物的波譜特征，對(duì)代表鐵染蝕變信息的主成分圖像的判斷準(zhǔn)則為:構(gòu)成該主成分的特征向量，其ASTER4的波段系數(shù)應(yīng)該與ASTER1，3的波段系數(shù)反向，與ASTER2的波段系數(shù)同向。經(jīng)過比較，ASTER4符合這一標(biāo)準(zhǔn)，但由于ASTER4在PC4上表現(xiàn)為負(fù)值，應(yīng)先對(duì)PC4進(jìn)行取反處理，取反后得到的圖像就可作為鐵染異常圖;再對(duì)鐵染異常圖像進(jìn)行低通濾波和分級(jí)處理(首先根據(jù)“平均值+n倍標(biāo)準(zhǔn)離差”的值確定出異常下限，然后對(duì)異常值進(jìn)行平均分配、確定出異常的一、二、三級(jí))。分級(jí)的結(jié)果為:191～212為鐵染三級(jí)異常，213～234為鐵染二級(jí)異常，235～255為鐵染一級(jí)異常。

3 遙感蝕變異常特征

從總體上看，遙感羥基異常與鐵染異常絕大部分重合，鐵染異常的分布范圍較羥基異常的大;羥基異常呈斷續(xù)的團(tuán)塊狀分布，鐵染異常則呈連續(xù)的帶狀分布。遙感鐵染異常與羥基異常都位于野外已知蝕變帶范圍之內(nèi)，與已知蝕變帶分布范圍基本吻合。

3.1 羥基異常特征

羥基異常以三級(jí)異常為主，二級(jí)異常次之，一級(jí)異常缺失。從整體上看，羥基異常分布的平面形態(tài)呈向西開口的半圓形，半圓北部的異常最多，兩個(gè)大型的團(tuán)塊狀異常呈北西向展布;半圓南部的異常少而分散，略呈北北東向展布(圖3)。

圖3 羥基蝕變異常特征Fig.3 The characteristics of hydroxyl alteration anomalies

遙感羥基異常全部落在野外調(diào)查見到的蝕變帶上，其中面積最大的團(tuán)狀異常的主體落在硅化蝕變帶內(nèi)，其他異常則落在了絹云母化+硅化(絹英巖化帶)蝕變帶內(nèi)和鉀長(zhǎng)石化+高嶺土化蝕變帶內(nèi);位于半圓南部呈帶狀分布的羥基異常都落在青磐巖化蝕變帶內(nèi)，且異常走向與蝕變帶界線的走向基本一致;此外在區(qū)內(nèi)的東北角有一面積很小的羥基異常落在絹云母化+硅化(絹英巖化帶)蝕變帶內(nèi)。

3.2 鐵染異常特征

鐵染異常主要分布在礦床的中部和中南部(圖4)。鐵染異常以三級(jí)為主，二級(jí)次之，一級(jí)鐵染異常缺失。鐵染異常分布的平面形態(tài)可以分為兩部分，北部呈似“碗”狀將已知礦體圍限;南部具有與羥基異常類似的分布特點(diǎn)，但異常分布范圍較羥基異常大。

圖4 鐵染蝕變異常特征Fig.4 The characteristics of iron alteration anomalies

3.3 遙感異常的應(yīng)用

通過對(duì)多不雜銅礦床遙感蝕變異常與野外調(diào)查見到的蝕變帶、已探明礦體的對(duì)比研究發(fā)現(xiàn)，已探明礦體基本上被遙感蝕變異常所圍限，類似的情況也出現(xiàn)在玉龍斑巖型銅礦床［16］。

通過對(duì)遙感蝕變異常的展布特征、遙感影像特征、蝕變帶實(shí)地分布及已探明礦體特征的綜合分析，本次工作推測(cè)出3處新的找礦靶區(qū)(見圖5中的A，B，C 區(qū))。

圖5 研究區(qū)遙感地質(zhì)特征解譯及推測(cè)找礦靶區(qū)分布圖Fig.5 The map of remote sensing interpretation and prospecting targets in study area

上述3處找礦靶區(qū)的推測(cè)依據(jù)見表4。

表4 找礦靶區(qū)推測(cè)依據(jù)Tab.4 The basis for prospecting targets in study area

其中，已知礦體平面投影范圍落在A區(qū)，根據(jù)遙感異常特征推測(cè)的A區(qū)范圍比已知礦體范圍大，因地面鉆探尚未完全控制礦體，A區(qū)可為地面勘探工程的布設(shè)提供重要的指導(dǎo)方向;B區(qū)和C區(qū)則為多不雜銅金礦床外圍(多龍銅礦預(yù)測(cè)工作區(qū))找礦工作縮小了找礦范圍。

4 結(jié)論

1)班公湖—怒江縫合帶是西藏的第三個(gè)重要成礦帶，由于受到自然條件的限制，目前的礦產(chǎn)勘查工作都是在點(diǎn)上進(jìn)行的。通過研究發(fā)現(xiàn)，利用ASTER數(shù)據(jù)、采取“掩模+定向主成分分析”的方法提取的遙感蝕變異常與野外調(diào)查見到的礦化蝕變范圍基本吻合。

2)在西藏礦產(chǎn)資源潛力評(píng)價(jià)的工作中，本次研究提取的遙感蝕變異常為多龍銅礦預(yù)測(cè)工作區(qū)最小預(yù)測(cè)區(qū)??最小預(yù)測(cè)區(qū)指根據(jù)遙感提取的構(gòu)造和礦化蝕變異常信息圈定的找礦遠(yuǎn)景區(qū)。的圈定提供了重要的遙感依據(jù)。為進(jìn)一步縮小(鎖定)找礦靶區(qū)、加快礦產(chǎn)資源勘查進(jìn)度，建議在班—怒成礦帶全面開展基于ASTER數(shù)據(jù)的遙感蝕變異常提取掃面工作。

［1］ 佘宏全，李進(jìn)文，豐成友，等.西藏多不雜斑巖銅礦床高溫高鹽度流體包裹體及其成因意義［J］.地質(zhì)學(xué)報(bào)，2006，80(9):1434－1447.

［2］ 佘宏全，李進(jìn)文，馬東方，等.西藏多不雜斑巖銅礦床輝鉬礦Re－Os和鋯石U－Pb SHRIMP測(cè)年及地質(zhì)意義［J］.礦床地質(zhì)，2009，28(6):737 －746.

［3］ 李光明，李金祥，秦克章，等.西藏班公湖帶多不雜超大型富金斑巖銅礦的高溫高鹽高氧化成礦流體:流體包裹體證據(jù)［J］.巖石學(xué)報(bào)，2007，23(5):935 －952.

［4］ 李金祥，秦克章，李光明，等.西藏班公湖帶多不雜富金斑巖銅礦床中金紅石的特征及其意義［J］.礦床地質(zhì)，2008，27(2):209－219.

［5］ 叢麗娟，岑 況，朱 所，等.利用ASTER數(shù)據(jù)提取蝕變異常方法研究—以內(nèi)蒙古朱拉扎嘎金礦為例［J］.河南理工大學(xué)學(xué)報(bào):自然科學(xué)版，2007，12(6):652 －658.

［6］ 楊金中，方洪賓，張玉君，等.中國西部重要成礦帶遙感找礦異常提取的方法研究［J］.國土資源遙感，2003(3):50－53.

［7］ 楊金中.多光譜遙感異常提取技術(shù)方法體系研究［J］.國土資源遙感，2007(4):43－46.

［8］ 張廷斌，唐菊興，黃丁發(fā).礦化蝕變信息提取的TM/ETM+遙感影像模式［J］.遙感信息，2009(2):47－51.

［9］ 陳建明，孫衛(wèi)東，閆柏琨，等.ASTER多光譜遙感異常提取在新疆天湖鐵礦中的應(yīng)用［J］.新疆地質(zhì)，2009，12(4):368 －372.

［10］張廷斌，鐘康惠，易桂花，等.東昆侖五龍溝金礦集中區(qū)遙感地質(zhì)信息提取與找礦預(yù)測(cè)［J］.地質(zhì)與勘探，2009，7(4):444 －449.

［11］毛曉長(zhǎng)，劉文燦，杜建國，等.ETM+和ASTER數(shù)據(jù)在遙感礦化蝕變信息提取應(yīng)用中的比較——以安徽銅陵鳳凰山礦田為例［J］.現(xiàn)代地質(zhì)，2005，6(2):309 －314.

［12］楊佳佳，姜琦剛，趙 靜，等.基于ASTER和ETM+數(shù)據(jù)的遙感蝕變信息提取——以內(nèi)蒙古塔日根敖包地區(qū)為例［J］.吉林大學(xué)學(xué)報(bào):地球科學(xué)版，2008，38(增刊):153 －154.

［13］耿新霞，楊建民，張玉君，等.ASTER數(shù)據(jù)在淺覆蓋區(qū)蝕變遙感異常信息提取中的應(yīng)用——以新疆西準(zhǔn)噶爾包古圖斑巖銅礦巖體為例［J］.地質(zhì)論評(píng)，2008，54(2):184 －190.

［14］呂鳳軍，郝躍生，石 靜，等.ASTER遙感數(shù)據(jù)蝕變遙感異常提取研究［J］.地球?qū)W報(bào)，2009，30(2):271 －276.

［15］黃照強(qiáng)，李祥強(qiáng).基于ASTER和ETM+數(shù)據(jù)的蝕變信息提取比較研究——以西藏澤當(dāng)?shù)V田為例［J］.地質(zhì)與勘探，2009，25(5):606－611.

［16］西藏自治區(qū)地質(zhì)調(diào)查院.西藏自治區(qū)礦產(chǎn)資源潛力評(píng)價(jià)玉龍銅礦典型礦床遙感鐵染異常說明書［R］.拉薩:西藏自治區(qū)地質(zhì)調(diào)查院，2010.

Characteristics of Remote Sensing Alternation Anomalies from ASTER in the Duobuza Porphyry Copper Deposit

HU Zi－h(huán)ao1，TANG Ju －xing2，ZHANG Ting－bin1，WU Hua3，XU Zhi－zhong4，BIE Xiao－juan5
(1.College of Earth Sciences，Chengdu University of Technology，Chengdu 610059，China;2.Institute of Mineral Resources，Chinese Academy of Geological Sciences，Beijing 100037，China;3.Institute of Geological Survey of Tibet，Lhasa 850000，China;4.No.5 Geological Party of Tibet Bureau of Geology and Mineral Resources，Geermu 816000，China;5.College of Tourism and Urban and Rural Planning，Chengdu University of Technology，Chengdu 610059，China)

Using the“mask and directed principal component analysis”approach，the authors extracted the remote sensing alternation anomalies of the Duobuza porphyry copper deposit on the basis of ASTER data.The characteristics of alteration anomalies，the correlation between alteration anomalies and surface alteration and the relationship between alteration anomalies and known ore bodies were preliminarily discussed in this paper.The results show that the alteration anomalies extracted by ASTER data are relatively satisfactory and can provide the remote sensing basis for delineation of the smallest prognostic target within the forecasting work area of the Duolong copper deposit in the mineral resource potential assessment of Tibet.

remote sensing;alternation anomaly;ASTER;Duobuza porphyry copper deposit

TP 79

A

1001－070X(2012)01－0150－05

10.6046/gtzyyg.2012.01.26

2011－02－24;

2011－05－23

中國地質(zhì)調(diào)查局全國礦產(chǎn)資源潛力評(píng)價(jià)西藏自治區(qū)礦產(chǎn)資源潛力評(píng)價(jià)項(xiàng)目(編號(hào):1212010813025)資助。

胡紫豪(1987－)，男，成都理工大學(xué)在讀碩士生，主要從事遙感地質(zhì)方面的研究。E－mail:mr.zihao@gmail.com。

張廷斌(1978－)，男，碩士，副教授，主要從事遙感教學(xué)與研究工作。E－mail:zhangtb@cdut.edu.cn。

(責(zé)任編輯:劉心季)