遙感技術在格咱島弧地區斑巖銅礦勘查中的應用

劉峻杉，胡 濱，何政偉，b

（成都理工大學a.地球科學學院，b.地質災害防治與地質環境保護國家重點實驗室，成都 610059）

遙感技術在格咱島弧地區斑巖銅礦勘查中的應用

劉峻杉a，胡 濱a，何政偉a，b

（成都理工大學a.地球科學學院，b.地質災害防治與地質環境保護國家重點實驗室，成都 610059）

遙感數據能夠提取圍巖蝕變信息，為斑巖銅礦的勘查提供依據。研究區地處格咱島弧地區，地質研究程度較低，成礦潛力巨大。采用主成分分析法處理OLI數據得到研究區的鐵染和羥基蝕變信息，采用該方法從ASTER數據中提取了絹云母信息。利用ASTER數據的SWIR波段來提取絹云母信息，能夠抑制植被信息，突出礦物信息。經過野外踏勘發現了斑巖銅礦點，證實了遙感技術的有效性，該方法能夠為云南地區斑巖銅礦的勘查提供借鑒。

遙感技術；OLI；ASTER；蝕變信息；格咱島弧地區

0 引言

多光譜遙感數據在蝕變礦物填圖方面具有不可替代的優勢，采用遙感技術進行鐵染與羥基類蝕變礦物、碳酸鹽礦物和硅酸鹽礦物填圖已經被廣泛采用。前人采用TM／ETM＋數據進行礦化蝕變信息提取的技術趨于成熟，所提取的蝕變圍巖是重要的找礦標志［1－3］。ASTER數據有3個波段在VNIR區間，空間分辨率是15m；有6個波段在SWIR區間，空間分辨率是30m。前人采用ASTER數據來識別明礬石、高嶺石、方解石、綠泥石、白云母等熱液蝕變礦物［4－8］，以上蝕變礦物對于斑巖銅礦的勘查具有重要的指示作用。

格咱島弧位于西南三江地區義敦島弧的南段，是甘孜——理塘洋殼向西俯沖而形成的。格咱島弧內發生過多期的構造運動和巖漿活動，豐富的成礦作用發育于島弧環境中［9］。格咱島弧熱液活動劇烈，成礦作用復雜，該火山巖漿弧以安山巖——英安巖及超淺成侵位的中酸性巖體為主要特征。該地區氣候復雜，交通不便，野外條件艱苦，地貌上具有山高谷深地形起伏大的特點，遙感技術在格咱地區的礦產勘查上具有獨特的技術優勢。

1 區域地質概況

格咱島弧是晚三疊世甘孜——理塘洋殼向德格——中甸陸塊俯沖的產物，是我國最近發現的重點銅多金屬成礦帶。格咱島弧東邊是甘孜－理塘結合帶，西邊是中咱——香格里拉地塊。西邊的邊界是格咱斷裂，東邊的邊界是黑水潭斷裂。

格咱島弧地區內出露地層主要為晚三疊系地層，為巖漿弧發育時期，該地區發育巨厚的碎屑巖——碳酸鹽巖——火山巖建造，巖性主要為砂板巖夾灰巖、安山玄武巖——安山巖、英安巖，可以劃分為喇嘛埡組（T3lm）、拉納山組（T3l）、圖姆溝組（T3t）、曲嘎寺組（T3q）。該地區北西向斷層控制了印支晚期擠壓型鈣堿性系列鉀質中酸性火山巖，并有大量同源中酸性淺成斑巖及次火山巖分布。格咱地區內出露的巖體是與格咱島弧火山巖同源、同時異相巖漿活動的產物。印支期島弧型中酸性淺成、超淺成巖漿巖是格咱島弧區域內最重要的巖漿巖類，該時期的巖漿巖以閃長斑巖、石英二長斑巖、花崗閃長斑巖為主，部分巖體蝕變強烈，產斑巖銅礦。研究區的地質圖如圖1所示。

圖1 研究區地質圖Fig.1 Geological map of the research area

2 OLI數據蝕變信息提取

本次研究采用的OLI數據的軌道號是132－40，獲取日期是2014年11月22號，數據級別是L1T，該區域的OLI數據只有少量的積雪覆蓋，影像質量合格。對該景數據進行輻射定標、FLAASH大氣校正、投影變換、裁剪、掩膜等處理，通過對OLI2、OLI4、OLI5、OLI6四個波段進行主成分變換來提取鐵染礦化蝕變信息，通過分析主成分變換后的特征向量矩陣（表1），可知PC4為鐵染圖像。對PC4進行3×3中值濾波，然后進行統計分析，在統計分析時要添加掩膜波段，可以避免背景值對圖像的特征值產生影響。通過均值＋N×標準差進行主分量異常等級閾值切割，采用N＝3、2.5、2進行一級、二級、三級鐵染異常的分割。從OLI數據中得到研究區的鐵染蝕變信息如圖2所示。

表1 OLI2，4，5，6波段主成分變換特征向量矩陣Tab.1 Eigenvector matrix of principal component transform for OLI band 2，4，5and 6

圖2 研究區鐵染蝕變信息Fig.2 Iron－stain alteration information of the research area

通過對OLI2、OLI5、OLI6、OLI7四個波段進行主成分變換來提取羥基礦化蝕變信息，通過分析主成分變換后的特征向量矩陣（表2），可知PC4為羥基圖像。對PC4進行3×3中值濾波，然后進行統計分析，在統計分析時要添加掩膜波段，通過均值＋N×標準差進行主分量異常等級閾值切割，采用N＝3、2.5、2進行一、二、三級羥基異常的分割。從OLI數據中得到研究區的羥基蝕變信息（圖3）。

表2 OLI2，5，6，7波段主成分變換特征向量矩陣Tab.2 Eigenvector matrix of principal component transform for OLI band 2，5，6and 7

圖3 研究區羥基蝕變信息Fig.3 Hydroxyl alteration information of the research area

由圖2與圖3可知，在研究區下部的查拉牛場一帶中鐵染異常信息和羥基異常信息呈面狀分布和北西向條帶狀分布，遙感蝕變信息通常分布于巖體或者隱伏巖體的周圍，與巖體以及巖漿熱液活動關系密切。遙感異常強度與巖體的埋藏深度有關系，巖體埋藏深度越淺，蝕變強度越大。在查拉牛場區域中一級鐵染異常信息與一級羥基異常信息呈大面積面狀分布，推測該區域存在埋藏深度較淺的巖體。查拉牛場區域附近的蝕變呈北西向條帶狀分布，與該地區的構造方向一致，遙感蝕變信息受到地形、地貌的影響，容易出現蝕變信息不連續的現象。

3 ASTER數據蝕變信息提取

ASTER數據在SWIR區間有6個波段，有利于識別不同類型的蝕變礦物［10］。白云母在2.20 μm處具有強吸收特征，在2.35μm處有次一級的吸收特征；高嶺石和明礬石的吸收谷位于2.17μm 和2.20μm處［11］。植被在短波紅外區間受到植物含水量的影響，反射率下降。采用ASTER的SWIR區間的6個波段來提取蝕變信息，能夠壓抑植被信息，突出礦物信息，適合中高植被覆蓋區蝕變信息提取。短波紅外數據的信噪比較低，采用最小噪聲分離（Minimum Nosie Fraction，MNF）處理能夠有效的去除噪聲，進而提取蝕變礦物信息［12－14］。這里首先對最小噪聲分離處理后的數據進行彩色合成來增強蝕變信息，然后分別采用主成分分析法和最小能量約束法（Constrained Energy Minimization，CEM）提取絹云母信息［15－18］，最后根據蝕變信息進行野外驗證。

對ASTER數據進行輻射定標、FLAASH大氣校正、幾何校正、裁剪、掩膜等預處理。采用最小噪聲分離處理ASTER數據的SWIR區間的6個波段，根據圖4可知，前三個波段的特征值遠大于后三個波段的特征值，前三個波段中集中了大部分的地物光譜信息。采用變換后的3、2、1波段進行彩色合成結果如圖5所示。

圖4 MNF特征值曲線Fig.4 MNF eigenvalue curve

對研究區的ASTER1、ASTER3、ASTER5、ASTER6四個波段進行主成分分析，其特征向量矩陣如表3所示。絹云母信息在ASTER 6波段有強吸收特征，由表3可知，PC4波段中高亮度信息表示絹云母信息，對PC4進行取反操作，然后進行3×3中值濾波、統計分析，通過均值＋N×標準差進行主分量異常等級閾值切割，采用N＝3、2.5、2進行一、二、三級絹云母信息的分割。從ASTER數據中得到研究區的絹云母蝕變信息如圖6所示。

表3 ASTER 1，3，5，6波段主成分變換特征向量矩陣Tab.3 Eigenvector matrix of principal component transform for Aster band 1，3，5and 6

圖5 ASTER短波紅外數據MNF變換后的3、2、1波段合成圖像Fig.5 RGB color composites of MNF Eigen images 1，2and 3extracted from the SWIR subsystem of ASTER data

圖6 采用主成分分析法從ASTER數據中提取的絹云母信息Fig.6 Sericite information extracted from ASTER data by using principle component analysis

將MNF1、MNF2波段作為X、Y軸構成二維散點圖，選擇散點圖周圍凸出部分區域，再獲取這個區域相應原圖上的平均波譜作為端元波譜。從ASTER數據中提取的絹云母波譜曲線如圖7所示。

根據ASTER數據中提取的絹云母端元波譜，采用最小能量約束法得到單波段的信息，最小能量約束法所得結果是每個端元波譜對比每個像元得到的灰度圖像，像元值越大表示越接近目標［19－20］。對最小能量約束法得到的單波段進行高端切割，將圖像亮度值大于0.52的區域提取出來得到研究區的絹云母分布信息如圖8所示。圖8中最小能量約束法所提取的絹云母信息，與圖6中主成分分析法所提取的絹云母信息分布范圍基本一致，這也表明不同方法在提取絹云母信息時的有效性。

4 野外驗證與結論

在1∶200 000貢嶺幅地質圖中，查拉牛場區域沒有巖體出露。受研究區地質工作程度的影響，小巖體容易被遺漏，嚴重制約了礦化斑巖體的發現［21］。查拉牛場成礦遠景區處于早堆——磨莫亞——亞雜——卓瑪銅成礦帶的下段，從遙感影像上所提取的圍巖蝕變信息呈大面積面狀分布，范圍較大，成礦潛力巨大。通過野外實地踏勘，在查拉牛場地區發現了強烈了圍巖蝕變現象，根據野外踏勘結果可知該區域存在一個填圖錯過的巖體。

圖7 從ASTER數據中獲取的絹云母波譜曲線Fig.7 Image average spectral curve of sericite from ASTER

根據最新的大比例尺地質填圖結果（圖10），查拉牛場地區存在石英閃長玢巖和石英二長斑巖組成的復式巖體，巖體呈NW－SE向延伸，石英閃長玢巖出露面積較大，長大于10km，寬1km～2km。石英二長斑巖侵位于石英閃長玢巖，出露面積約0.8km2。羥基蝕變信息主要分布于絹英巖化帶中，少部分分布于青磐巖化帶。該地區巖石經過了區域低溫動力變質作用和熱接觸變質作用，巖體外接觸帶形成角巖或角巖化巖石。

圖8 CEM提取的絹云母信息Fig.8 Sericite information extracted by CEM

圖9 實地圍巖蝕變信息Fig.9 The image of mineralized altered rocks

通過以上研究，作者得到以下結論：

1）云南地區地處我國西部，地貌上具有山高谷深地形起伏大的特點，地質研究程度較低，在區域地質調查過程中容易遺漏礦化斑巖。斑巖銅礦具有圍巖蝕變強烈的特點，采用遙感技術來提取圍巖蝕變信息，可以為斑巖銅礦的勘查提供依據。

2）采用傳統的主成分分析法從OLI數據中提取了鐵染蝕變信息和羥基蝕變信息，并分析了蝕變信息的特點。研究區具有植被發育的特點，采用ASTER數據SWIR波段來提取絹云母化信息，能夠有效地壓抑植被信息，突出蝕變礦物信息。

圖10 查拉牛場地質圖Fig.10 Geological map of Chalaniuchang

3）通過野外踏勘和大比例尺的地質填圖結果，在查拉牛場地區發現了礦化巖體信息，證明了遙感技術在云南地區斑巖銅礦勘查中的實用性與有效性。

［1］ 姚偉，劉亮明.遙感技術在北非努比亞地盾Wadi Halfa地區金礦勘查中的應用［J］.地質與勘探，2014，50（1）：167－172

YAO W，LIU L M.Application of the remote sensing technique to the gold exploration in Wadi Halfa district，Nubia shield of north Africa［J］.Geology and Exploration，2014，50（1）：167－172.（In Chinese）

［2］ 楊建民，張玉君，姚佛軍，等.遙感找礦信息在新疆羅東鎳礦發現中的主導作用［J］.巖石學報，2007，23（10）：2647－2652.

YANG J M，ZHANG Y J，YAO F J，et al.Dominant role of remote sensing mineral exploration information in the discovery of the Luodong Ni deposit，Xinjiang ［J］.Acta Petrologica Sinica，2007，23（10）：2647－2652.（In Chinese）

［3］ 張玉君，楊建民.基巖裸露區蝕變巖遙感信息的提取方法［J］.國土資源遙感，1998，（2）：46－53

ZHANG Y J，YANG J M.The method of abstracting remote sensing information of alterated rocks in theuncovered bedrocks area［J］.Remote Sensing for Land＆Resources，1998（2）：46—53.（In Chinese）

［4］ ROWAN L C，HOOK S J，ABRAMS M J，et al.Mapping hydrothermally altered rocks at Cuprite，Nevada，using the advanced spaceborne thermal emission and reflection radiometer（ASTER）［J］.A new satellite－imaging system：Economic Geology and the Bulletin of the Society of Economic Geologists，2003，98 （5）：1019－1027.

［5］ SAFWAT GABR，ABDUWASIT GHULAM，TIMOTHY KUSKY，Detecting areas of high—potential gold mineralization using ASTER data［J］.Ore Geology Reviews 2010，38：59－69.

［6］ SALEM S M，SOLIMAN N M，RAMADAN T M.Exploration of new gold occurrences in the alteration zones at the Barramiya district，central eastern desert of egypt using ASTER data and geological studies［J］.Arab J Geosci，2014（7）：1717－1731

［7］ 代晶晶，曲曉明，辛洪波.基于ASTER遙感數據的西藏多龍礦集區示礦信息的提取［J］.地質通報，2010，29 （5）：752－759

DAI J J，QU X M，XIN H B.Extraction of alteration mineral information using ASTER remote sensing data in Duolong area，Tibet，China［J］.Geological Bulletin of China，2010，29（5）：752—759.（In Chinese）

［8］ 毛曉長，劉文燦，杜建國，等.ETM＋和ASTER數據在遙感礦化蝕變信息提取應用中的比較—以安徽銅陵鳳凰山礦田為例［J］.現代地質，2005，19（2）：309—314

MAO X CH，LIU W C，DU J G，et al.Comparison between ETM＋and ASTER data for extraction of alteration information：a case study of fenghuangshan orefield，Tongling Anhui province［J］.Geoscience，2005，19（2）：309－314.（In Chinese）

［9］ 劉學龍.中甸島弧巖漿巖地球化學特征及斑巖型銅礦成礦研究［D］.昆明：昆明理工大學，2009

LIU X L.Zhongdian island arc magmatic rock geochemical characteristics and porphyry copper metallogenic research［D］.Kunming：Kunming University of Science and Technology，2009.（In Chinese）

［10］ROWAN，L.C.，MARS，J.C.，Lithologic mapping in the Mountain Pass，California，area using advanced spaceborne thermal emissionand reflection radiometer（ASTER）data［J］.remote sensing of environment，2003，84（3）：350－366.

［11］JOHN C，MARS，LAWRENCE C.ROWAN.Regional mapping of phyllic－and argillic－altered rocks in the Zagros magmatic arc，Iran，using Advanced Spaceborne Thermal Emission and Reflection Radiometer（ASTER）data and logical operator algorithms［J］.Geosphere，2006，2（3）：161－186.

［12］王愛云，王安建，李麗輝.基于最大噪聲分量變換（MNF）和礦物標識的植被區蝕變信息提取［J］.地質與勘探，2011，47（4）：710－718

WANG AI Y，WANG AN J，LI L H.mapping mineralogical alteration using MNF transform and mineral in vegetated areas［J］.Geology and Exploration，2011，47 （4）：710－718.（In Chinese）

［13］胡濱，何政偉，王銳，等.遙感技術在西藏朱諾地區巖性識別中的應用［J］.桂林理工大學學報，2014，34（4）：123－131

HU B，HE ZH W，WANG R，et al.Application of remote sensing technique for lithologic identification of the Zhunuo region，Tibet［J］.Journal of Guilin University of Technology，2014，34（4）：123－131.（In Chinese）

［14］GREEN A A，BERMANM，SWITZER P.A transformation for ordering multispectral data in terms of image quality with implications for noise removal［J］.IEEE Transactions on Geoscience and Remote Sensing，1988，26（1）：65－74.

［15］耿新霞，楊建民，張玉君，等.ASTER數據在淺覆蓋區蝕變遙感異常信息提取中的應用—以新疆西準葛爾包古圖斑巖銅礦巖體為例［J］.地質論評，2008，54（2）：184－191.

GENG X X，YANG J M，ZHANG Y J，et al.The application of ASTER remote sensing data for extraction of alteration anomalies information in shallow Overburden area—a case study of the Baoguto porphyry copper deposit intrusion in western Junggar，Xinjiang［J］.GEOLOGICAL REVIEW.2008，54（2）：184－191.（In Chinese）

［16］黃照強，李祥強.基于ASTER和ETM＋數據的蝕變信息提取比較研究—以西藏澤當礦田為例［J］.地質與勘探，2009，45（5）：606－611

HUANG ZH Q，LI X Q.Extraction of alteration anomalies based on ASTER and ETM data：a case study of Tsedang orefield，Tibet［J］.Geology and Exploration，2009，45（5）：606—611.（In Chinese）

［17］郭娜，郭科，胡敬仁，等.多源遙感分析技術在西藏唐噶幅找礦中的應用［J］.國土資源科技管理，2012，29（6）：70－74.

GUO N，GUO K，HU J R，et al.Application of integrated remote sensing techniques to mineral exploration in Tangga，Tibet［J］.Scientific and Technological Management of Land and Resources，2012，29（6）：70－74.（In Chinese）

［18］張玉君，曾朝銘.西藏岡底斯地區斑巖銅礦識別的最佳多光譜遙感異常判據研究［J］.礦床地質，2012，31（4）：671－698

ZHANG Y J，ZENG ZH M.Study of optimum discriminator of multi—channel remote sensing anomaly for recognition of porphyry Cu deposits in Gangdese belt ［J］.MINERAL DEPOSITS，2012，31（4）：671－698.（In Chinese）

［19］SAFWAT GABR，ABDUWASIT GHULAM，TIMOTHY KUSKY.Detecting areas of high－potential gold mineralization using ASTER data［J］.Ore Geology Reviews，2010，38：59－69.

［20］FARRAND W H，HAESANYI J C.Mapping distributed geological and botanical targets through constrained energy minimization［C］.Proceedings，10th Thematic Conference on Geological Remote Sensing，San Antonio，1994：1419－1429.

［21］尹芳，李前志，周軍，等.新疆加甫沙爾蘇地區中酸性小巖體型礦化蝕變遙感異常提取與找礦［J］.遙感技術與應用，2014，29（1）：54－60

YIN F，LI Q ZH，ZHOU J，et al.Remote sensing extraction of intermediate－acid small intrusions related alteration anomalies and mineral exploration in Jiafushaersu area，Xinjiang［J］.Remote Sensing Technology and Application，2014，29（1）：54－60.（In Chinese）

The application of remote sensing technique for porphyry copper exploration in the Gezan arc region

LIU Jun－shana，HU Bina，HE Zheng－weia，b
（Chengdu University of Technology a.College of Earth Science，b.State Key Laboratory of Geohazard Prevention and Geoenvironment Protection，Chengdu 610059，China）

Remote sensing data can be used to extract alteration information，which provides the basis for porphyry copper ore exploration.The study area is located in the Gezan arc region.The level of geological research is low and the metallogenic potential is tremendous.Principal component analyses（PCA）were applied to get the iron－stain alteration and hydroxyl alteration information by making use of OLI data.Principal component analyses（PCA）were applied to get sericite information from ASTER data.SWIR band of ASTER data is used to extract sericite information，the method can inhibit the vegetation information，highlight the mineral information.A porphyry mineralization point was discovered during ground inspection which has confirmed the validity of remote sensing technology.The method can provide a reference for the mineral exploration in Yunnan region.

remote sensing technology；OLI；ASTER；alteration information；Gezan arc region

TP 79

：ADOI：10.3969／j.issn.1001－1749.2015.06.20

1001－1749（2015）06－0790－07

2015－08－23改回日期：2015－10－10

國家自然科技基金（40972225）；中國地質調查局項目（12120113095400）；教育部博士點基金（20095122110003）；成都理工大學礦物學巖石學礦床學國家重點（培育）學科建設項目（SZD0407）

劉峻杉（1980－），男，博士，主要研究方向為地球探測與信息技術，E－mail：77729166＠qq.com。