白 卉

吉林省地質(zhì)調(diào)查院，吉林 長(zhǎng)春 130061

吉林臨江銅礦地區(qū)遙感礦化蝕變信息提取

白 卉

吉林省地質(zhì)調(diào)查院，吉林 長(zhǎng)春 130061

以吉林省臨江銅礦附近地區(qū)為研究區(qū)，選取OLI多光譜數(shù)據(jù)為基本數(shù)據(jù)源，充分利用遙感數(shù)據(jù)多時(shí)相、多波段、多空間分辨率等特點(diǎn)，結(jié)合已知地質(zhì)資料和已知礦點(diǎn)，采取主成分分析和比值法相結(jié)合，提取遙感礦化蝕變信息，得到具有成礦潛力的區(qū)域。

臨江銅礦；遙感；礦化蝕變信息；OLI數(shù)據(jù)

0　引言

森林覆蓋區(qū)的大面積遙感地質(zhì)礦產(chǎn)調(diào)查與評(píng)價(jià)一直以來(lái)都是地質(zhì)研究的熱點(diǎn)，特別是利用遙感方法在不適用于常規(guī)地質(zhì)調(diào)查方法的植被覆蓋區(qū)確定找礦遠(yuǎn)景區(qū)更是其中的難點(diǎn)。遙感礦化蝕變信息提取是利用光譜異常為依據(jù)進(jìn)行提取。鐵的氧化物、氫氧化物和硫酸鹽三類(lèi)蝕變礦物理論上能被OLI圖像識(shí)別，而蝕變礦物中大多含有羥基（OH-），鐵染（Fe3+）則主要與礦化有關(guān)。

1　研究區(qū)地質(zhì)背景

研究區(qū)共涉及四幅1/5萬(wàn)圖幅（圖1），分別為螞蟻河幅（K52E014005）、東北岔幅（K52E014006）、六道溝公社幅（K52E015005）和臨江銅礦幅（K52E015006），地處鴨綠江斷裂帶的西段，六道溝—臨江銅礦斷裂帶呈北東向通過(guò)研究區(qū)中部。 區(qū)內(nèi)地層出露較豐富，自下而上主要為古元古界老嶺（巖）群的一套海相碎屑-碳酸鹽巖和新元古界青白口系長(zhǎng)石石英砂巖及頁(yè)巖組成；古生界地層主要為灰?guī)r和白云質(zhì)灰?guī)r組成；中生界晚三疊世長(zhǎng)白組出露范圍較廣，下段為安山巖和安山質(zhì)凝灰?guī)r，上段為流紋巖和流紋質(zhì)凝灰?guī)r；區(qū)內(nèi)西部為大面積第四系玄武巖覆蓋。

2　遙感礦化蝕變信息提取

首先選取OLI數(shù)據(jù)11631景，鑒于東北地區(qū)的特點(diǎn)，選取接收日期為2013年10月18日數(shù)據(jù)，盡可能減少植被和冰雪的干擾。為了保持?jǐn)?shù)據(jù)的原始性，不進(jìn)行幾何校正及掩膜處理。

3.1主成分分析

首先利用主成分分析法得到表1、表2：

經(jīng)過(guò)主成分分析，由上表可知，第一主成分可代表7個(gè)波段全部信息的91.571%，第二、三主成分貢獻(xiàn)率分別為4.272%和3.599%，均較高，三者之和占全部信息的99.442%，其代表的信息可作為背景信息；而第六和第七主成分貢獻(xiàn)率低，二者之和僅占全部信息的0.032%，不具實(shí)際意義，因此，該地區(qū)的礦化信息主要集中在第四和第五主成分中。

圖1 臨江銅礦附近地區(qū)地質(zhì)圖Fig.1　The geological map of Linjiang copper mine

表1　臨江銅礦附近地區(qū)OLI數(shù)據(jù)PCA特征值、貢獻(xiàn)率Table 1 Eigenvalue and contribution rate of OLI image principal components analysis of in Linjiang copper mine

表2　臨江銅礦附近地區(qū)OLI數(shù)據(jù)PCA特征向量矩陣Table 2 Eigenvector matrix of OLI image principal components analysis in Linjiang copper mine

根據(jù)上表中第四和第五主成分的特征向量矩陣可以看出，一波段（B1）和四波段（B4）在第四主成分中存在序偶關(guān)系；六波段（B6）和七波段（B7）在第五主成分中存在序偶關(guān)系。

3.2二維散點(diǎn)圖結(jié)構(gòu)層次分離

利用B4和B1，B6和B7分別形成二維散點(diǎn)圖（圖2、圖3）如下：

圖3　臨江銅礦附近地區(qū)OLI B6/B7二維散點(diǎn)圖Fig.3　Scatter plot of OLI B6/B7 of Linjiang copper mine

根據(jù)二維散點(diǎn)圖層次結(jié)構(gòu)分析，礦化異常必須滿(mǎn)足B4/B1（或B6/B7）大于1，礦化異常區(qū)應(yīng)位于圖中Y=X線(xiàn)上方的點(diǎn)陣，同時(shí)又高于背景值即主橢圓，經(jīng)過(guò)試驗(yàn)分析，最終確定礦化異常區(qū)大致位置。

3.3蝕變異常圖

分別對(duì)B1和B4波段，B6和B7波段做主成分分析，并分別對(duì)二者分析后的第二主成分采用“密度分割”，對(duì)蝕變信息作進(jìn)一步的分離，最終得到蝕變異常圖（圖4、圖5）。

圖4　臨江銅礦附近地區(qū)鐵染蝕變信息異常圖Fig.4　Fe3+alteration information of Linjiang copper mine

圖5　臨江銅礦附近地區(qū)羥基蝕變信息異常圖Fig.5　OH-alteration information of Linjiang copper mine

4　結(jié)論

由圖4、圖5可見(jiàn)臨江銅礦附近鐵染和羥基異常都有一定的反映，且二者吻合性較好；羥基異常在六道溝—臨江銅礦構(gòu)造帶附近呈線(xiàn)性分布，同時(shí)八道溝兩側(cè)羥基異常較好，局部有鐵染異常，與臨江銅礦附近提取結(jié)果具有相似性，具有較好的找礦潛力，其他地區(qū)吻合程度較低。

該方法未對(duì)研究區(qū)進(jìn)行幾何校正和掩膜處理，較好的保持了數(shù)據(jù)的原始性和完整性，且礦化蝕變信息的提取結(jié)果并未受到植被、水體、陰影和玄武巖等較大影響，對(duì)研究區(qū)高植被覆蓋的礦化信息提取具有重要的意義；但該方法也忽略了巖性本身的背景信息，如何能更好的結(jié)合地層信息進(jìn)行礦化信息提取還有待進(jìn)一步研究。

［1］ 楊長(zhǎng)保，姜琦剛.遼東地區(qū)礦化蝕變遙感信息提取的研究和應(yīng)用［J］ .遙感應(yīng)用，2007（4）.

［2］ 楊金中.遙感異常篩選的方法和意義［C］ .全國(guó)國(guó)土資源與環(huán)境遙感技術(shù)應(yīng)用交流會(huì)論文文集，2004.

［3］ 張遠(yuǎn)飛，等.基于物理意義的二維散點(diǎn)圖類(lèi)型劃分與遙感蝕變信息提取［J］ .國(guó)土資源遙感，2013（2）.

［4］ 張玉君，等.ETM+（TM）蝕變遙感異常提取方法技術(shù)［J］.航天遙感圖像數(shù)據(jù)索引及應(yīng)用實(shí)例，2004：37-59.

Remote sensing mineralization and alteration information extraction in the Linjiang copper mine area of Jilin Province

BAI Hui
Institute of Geological Survey of Jilin Province， Changchun 130061， Jilin， China

This paper studies in the Linjiang copper mine， Jilin Province， useing OLI multi-spectral data， which are characterized by multi temporal， multiband and multi-spatial resolution.According to the geological data and deposits，took the method of principal components analysis and ratio to extract mineralized alteration information with remote sensing and carried out the metallogenic potential area.

Linjiang copper mine； remote sensing； mineral alteration information； OLI image

P627

B

1001—2427（2015）04 - 103 -4

2015-07-05；

2015-11-17

白 卉（1987—），女，吉林長(zhǎng)春人，吉林省地質(zhì)調(diào)查院助理工程師.