趙　棟

(西北大學地質學系，陜西 西安 710069)

?

內蒙古石尖山地區(qū)遙感蝕變信息提取

趙棟

(西北大學地質學系，陜西 西安 710069)

[摘要]內蒙古石尖山地區(qū)是內蒙古北部重要的成礦帶，已發(fā)現(xiàn)一批鐵、鋅、鉛、銀、鉬、螢石成礦靶區(qū)，具有優(yōu)越的成礦地質背景和巨大的成礦潛力。石尖山地區(qū)屬于高原山區(qū)，植被覆蓋率高，地面地質工作開展難度較大，因此，在該地區(qū)開展遙感找礦研究是一種經濟便捷且有效的勘察手段。以ETM+數(shù)據(jù)為基礎，在ENVI軟件的支持下，根據(jù)與蝕變礦物組合相關的波譜特征，通過波段比值、主成分分析、彩色空間變換等方法，對石尖山地區(qū)鐵染、羥基、金礦化等蝕變信息進行提取，并結合基礎地質資料，劃定預測靶區(qū)，取得了很好的效果。

[關鍵詞]遙感；波段比值；主成份分析；礦化蝕變信息；石尖山地區(qū)

石尖山地區(qū)位于大興安嶺中段西坡，大興安嶺中生代火山巖區(qū),大地構造背景屬于西伯利亞板塊東南陸緣增生帶，區(qū)內褶皺、斷裂構造和火山機構十分發(fā)育，是內蒙北部重要成礦帶，區(qū)域上已發(fā)現(xiàn)的主要礦產為鐵、鋅、鉛、銀、鉬、螢石，其次有銅、石灰?guī)r、沸石、粘土、煤等，具有優(yōu)越的成礦地質背景和巨大成礦潛力。

調查區(qū)位于內蒙古東部大興安嶺中北段，地貌上屬高原山區(qū)，海拔高度一般800～1 100 m，相對高差較大。屬大陸性亞寒帶季風氣候，溫差較大。工作區(qū)屬森林沼澤景觀區(qū)，水系、沼澤、濕地發(fā)育，基巖露頭差。加之調查區(qū)多為山區(qū)，交通不便，地面地質工作的開展難度大,成本高,因此在該地區(qū)開展遙感找礦工作、劃定研究靶區(qū)，可以大大提高工作效率，節(jié)約成本。本文利用遙感數(shù)據(jù)對內蒙古石尖山地區(qū)進行礦化蝕變信息提取與分析,根據(jù)ETM+數(shù)據(jù)與蝕變礦物組合相關的波譜信息,利用掩膜+波段比值+主成分分析+色度調整的方法進行了蝕變異常信息提取,為該區(qū)成礦預測研究與建立找礦靶區(qū)提供了依據(jù)。

1區(qū)域成礦地質背景

本成礦區(qū)屬西伯利亞板塊南緣大陸邊緣，經歷了陸殼分離、古亞洲洋的生成、發(fā)展和消亡，以及濱西太平洋構造域的疊加，構造運動演化歷史復雜，巖漿作用強烈，因此成礦帶具有多期次疊加特點。

1.1構造控礦作用

本區(qū)中部有北東向特尼河大斷裂通過，南側外圍有北東向海拉爾河大斷裂，西南側有北西向勉渡河大斷裂，這些斷裂相互交匯，次級構造極為發(fā)育。區(qū)內的斷裂構造活動時間長、活動期次多、次級斷裂發(fā)育，常形成不同時代的礦床。晚侏羅世至早白堊世這些斷裂活動達到頂峰，并伴隨有強烈的中—酸性巖漿侵入活動，形成一系列的火山噴發(fā)巖漿帶，并在附近形成多個鐵、鉬、鉛鋅多金屬及螢石礦(床)點。

1.2巖漿巖控礦作用

調查區(qū)內見有泥盆紀、石炭紀和侏羅紀三個旋回巖漿侵入活動，形成近北東向展布的巖漿巖帶。石炭紀的巖漿活動形成于西伯利亞古板塊和華北古板塊的碰撞、擠壓環(huán)境及其后裂解再拉張環(huán)境，有大量的花崗巖類侵入，與本期侵入巖有關的金屬礦產主要分布在調查區(qū)東南部大春日山鉛鋅銀多金屬礦點。

燕山期是本區(qū)巖漿活動的高峰期，形成花崗質的巖漿活動，形成淺成—超淺成相小侵入體，并伴有Cu、Pb、Zn、Ag、Au、Mo等的礦化，區(qū)域上形成許多具有工業(yè)價值的礦床。燕山期侵入巖出露不多，侵入于早石炭世花崗巖體中，形成復合花崗巖帶，與區(qū)域該期巖漿的分布趨勢一致。

2蝕變信息提取

通過對目前各種常用方法的試驗以及優(yōu)劣性分析，我們認為針對不同蝕變類型先用特征波段比值、選擇主成分分析、彩色空間變換增強，然后經過閥值分割、異常疊加、假彩色合成、非監(jiān)督分類等方法的處理，在此基礎上結合主分量門限化模型對蝕變信息提取并分級，最終生成分級預測圖效果較好。基于該方法，利用ETM+數(shù)據(jù)對研究區(qū)礦化蝕變信息進行了提取，獲取了多種礦化蝕變信息成果。

2.1數(shù)據(jù)來源

本次研究采用Landsat-7 ETM+數(shù)據(jù)，景號為122/25和122/26。原始影像已經經過了幾何校正與輻射校正。為保持蝕變信息提取的完整性，本次研究對圖像數(shù)據(jù)進行了圖像鑲嵌與裁剪處理。在提取的蝕變圖像中，由于有第四系覆蓋、植被等干擾信息的影響，可以采用掩膜處理的方法對這些干擾信息進行剔除。

2.2綠泥石化蝕變信息提取

綠泥石在中低溫熱液礦床中廣泛分布。其為Mg-OH類蝕變礦物的典型代表，在2.32 μm附近具有一個明顯的吸收谷，在2.17 μm附近有一個反射峰。

選取ETM+1、4、5、7 四個波段，對其進行主成份分析。泥化蝕變分別在 ETM+5 和 ETM+7 有較強的反射和吸收，符號相反且絕對值較大，由此獲得“羥基圖像”，它反映富泥化蝕變的遙感異常。對代表泥化物主分量的判斷準則是：構成該主分量的本征向量，其中ETM+5系數(shù)應與ETM+7及ETM+4 的系數(shù)符號相反，ETM+1一般與ETM+5系數(shù)符號相同(張玉君等，2003)。據(jù)此，對 ETM+1、4、5、7 四個波段主成分分析結果進行統(tǒng)計，得出 PC4 代表了泥化信息。

對于主成分變換的第四分量 PC4 進行統(tǒng)計。據(jù)此，對 PC4 進行閾值分割，得到泥化信息圖像。最后，對結果進行中值濾波處理，得到最終的綠泥石化信息圖像。將生成的綠泥石化信息圖像矢量化；然后與單波段圖像疊加，生成綠泥石化信息提取遙感影像圖(見圖1)。

2.3羥基蝕變信息提取

在ETM+影像上，羥基信息在ETM+5 (1.55～1.75 μm)處表現(xiàn)出較強的反射，在ETM+7(2.08～2.35 μm)處表現(xiàn)出較強的吸收；此外，在ETM+1(0.45～0.52 μm)和 ETM+4(0.76～0.90 μm)處還有兩個較弱反射和吸收帶。鑒于羥基信息這一光譜特征，我們采用主成分分析的方法對工作區(qū)羥基蝕變信息進行遙感提取。

選擇ETM+1、4、5、7波段主成分分析，只有PC4這一主分量中ETM+1、ETM+5符號相同且與ETM+4和ETM+7符號相反，符合羥基蝕變信息的ETM+光譜特征，所以PC4是羥基信息的表征分量，我們利用PC4提取羥基信息。在第4主成分中ETM+5波段的權重因子最大，主要反映了羥基蝕變信息在第5波段的高反射性。由于第4主成分中還含有少量植被、河流及第四紀等干擾因素，需要將其去除。根據(jù)公式：PC4(1457)- ETM+4/3 - PC2(123457)(張廷斌，2009)對干擾因素進行去除。

對處理結果進行低通濾波處理后，統(tǒng)計分析得到最終的羥基蝕變信息圖像。將生成的羥基蝕變信息圖像矢量化；然后與單波段圖像疊加，生成羥基蝕變信息提取遙感影像圖(見圖1)。

2.4鐵染蝕變信息提取

在ETM+影像上，鐵染信息在ETM+1(0.4～0.5 μm)和ETM+4(0.8～1.0 μm)波段處有兩個明顯的吸收帶，在ETM+3(0.63～0.69 μm)波段處附近反射相對較高，可將此獨特的波譜特征作為鐵染信息蝕變提取的依據(jù)。在不考慮符號的情況下，ETM+3波段的特征向量載荷因子最大的主分量影像反映了含鐵離子蝕變礦物信息，如果ETM+3波段的特征向量載荷因子為負，則鐵離子分布區(qū)域在主分量影像中為低值異常。

選擇ETM+1、3、4、5波段主成分分析，只有PC4這一主分量中ETM+3系數(shù)應與ETM+1及 ETM+4的系數(shù)符號相反，且ETM+3與ETM+5系數(shù)符號相同。在第4主成分中ETM+3波段的權重因子最大，主要反映了羥基蝕變信息在第3波段的高反射性。ETM+3/1、ETM+5/1和ETM+5/4波段比值圖像不同程度的增強了黃鐵礦化、褐鐵礦化以及鉀化等具體與鐵染蝕變相關的蝕變信息；同時，波段比值ETM+4/3可以一直植被信息，因此選擇ETM+3/1、ETM+5/1、ETM+5/4、PC4(1345)和ETM+4/3圖像進行二次主成分分析，將結果中的第一主分量作為鐵染蝕變基礎圖，進行蒂通濾波處理后，統(tǒng)計分析結果得到最終的鐵染蝕變信息圖像(張廷斌，2009)。將生成的鐵染蝕變信息圖像矢量化，生成鐵染蝕變信息提取遙感影像圖(見圖1)。

2.5金礦化蝕變信息提取

與金礦礦化蝕變有關的礦物主要為一些含羥基(OH-)的粘土礦物(絹云母化、綠泥石化等)和富含三價鐵離子(Fe3+)的礦物。因此，在ETM+影像上，只有羥基礦物和三價鐵離子礦物同時顯示異常的區(qū)域才有可能是與金礦化有關的部位。

選擇上述羥基蝕變信息提取ETM+1、4、5、7波段主成分分析結果中的PC4分量、上述鐵染蝕變信息提取ETM+1、3、4、5波段主成分分析結果中的PC4分量、ETM+3/ETM1比值圖像(可以增強鐵染蝕變信息)以及ETM5/ETM7比值圖像(可以增強羥基蝕變信息)進行主成分分析，主要信息集中在 PC1 和 PC2 兩個組分中。其中 PC1 組分的信息主要來自“鐵組分”和“羥基組分”，且符號相反，是4個PC 組分中能最好地突出礦化蝕變信息的特征組分。因而選擇信息量最豐富的 PC1 組分進行金礦化蝕變信息提取。

對 PC1 特征組分進行線性反差增強、低通濾波、密度分割和分類處理，壓抑乃至去除與蝕變無關的信息，從而形成了只有金礦化蝕變信息圖像。將以上生成的圖像矢量化，生成金礦化信息矢量提取圖像，生成金礦化信息提取遙感影像圖(見圖1)。

3找礦預測

礦產的形成是多種成礦因素綜合作用的結果，研究區(qū)地質構造、地層巖性等都是地質歷史演化的結果，也是礦化信息提取的重要參考因素。因此，作者基于研究區(qū)遙感礦化蝕變提取圖并結合研究區(qū)構造、巖性、地球物理、地球化學和已知礦點等資料，圈定了幾個成礦預測區(qū)。同時，考慮到地質成礦因素與提取區(qū)的重合情況，將其分為了三個等級(見圖1)。

圖1　內蒙古石尖山地區(qū)ETM+遙感綜合信息提取預測圖

一級預測區(qū)：包括兩個預測區(qū)：Ⅰ1-1位于調查區(qū)西北部，呈北東向展布，預測區(qū)面積約140.07 km2，在調查區(qū)范圍內面積約112.06 km2；Ⅰ1-2位于調查區(qū)中北部，預測區(qū)面積約89.48 km2，兩異常區(qū)都與土壤異常、磁異常相重疊，且均位于線環(huán)構造的交切部位及兩側，區(qū)內已知礦點多，礦化信息提取結果表現(xiàn)為多種蝕變礦物的高異常。

二級預測區(qū)：包括兩個預測區(qū)：Ⅱ2-1位于調查區(qū)東南部，預測區(qū)面積約86.82 km2；Ⅱ2-2位于調查區(qū)東南角，預測區(qū)面積約41.98 km2，在調查區(qū)范圍內面積約27.99 km2，兩異常區(qū)都與土壤異常、磁異較為重疊，且均位于線環(huán)構造的交切部位及兩側，有已知礦化點分布，礦化信息提取結果表現(xiàn)為多種蝕變礦物的中異常。

三級預測區(qū)：包括三個預測區(qū)：Ⅲ3-1位于調查區(qū)西部，預測區(qū)面積約：46.75 km2，在調查區(qū)范圍內面積約36.73 km2；Ⅲ3-2位于調查區(qū)中西部，預測區(qū)面積約29.49 km2；Ⅲ3-3位于調查區(qū)東北部，預測區(qū)面積約95.78 km2，在調查區(qū)范圍內面積約38.31 km2，此預測區(qū)屬于未發(fā)現(xiàn)礦化點或礦化點較少的區(qū)域，但這些區(qū)域多種蝕變礦物集中分布，是今后找礦的潛在研究區(qū)。

4結語

從蝕變信息提取圖像來看，整個石尖山地區(qū)的礦化蝕變明顯是受到斷裂、環(huán)狀等構造現(xiàn)象所控制的，一級預測區(qū)往往出現(xiàn)在線環(huán)構造的交切部位，這類地區(qū)往往是熱液礦床的產出部位，而劃出的預測區(qū)也與已知的礦點、礦化點有很好的吻合，說明本次遙感蝕變信息提取結果有較好的可信性。

在常規(guī)找礦方法有較大難度的石尖山地區(qū)，本次研究主要采用基于波段比值、選擇主成分分析、彩色空間變換增強的方法，然后經過閥值分割、異常疊加、假彩色合成、非監(jiān)督分類等方法的處理，對研究區(qū)進行了蝕變信息提取，提取結果與區(qū)域斷裂展布方向基本一致，說明本次研究有一定的參考價值。但由于一些植被水體等干擾因素的存在，加之同物異譜與異物同譜現(xiàn)象的干擾，肯定存在一些錯漏信息，因此在進行礦化信息選擇時，須結合研究區(qū)的其它地質資料綜合分析，去除干擾信息后方可進行下一步工作。

參考文獻

[1]趙元洪,張福祥,陳南峰.波段比值的主成分復合在熱液蝕變信息提取中的應用[J].國土資源遙感.1991,(3):12-17.

[2]馬建文.利用TM數(shù)據(jù)快速提取含礦蝕變帶方法研究[J].遙感學報.1997,1(3):208-213.

[3]張遠飛,吳健生.基于遙感圖像提取礦化蝕變信息[J].有色金屬礦產與勘查.1999,8(6):604-606.

[4]劉慶生,燕守勛,馬超飛，等.內蒙哈達門溝金礦區(qū)山前鉀化帶遙感信息提取[J].遙感技術與應用.1999,14(3):7-11.

[5]劉素紅,馬建文,藺啟忠.通過Gram-Schmidt投影方法在高山區(qū)提取TM數(shù)據(jù)中含礦蝕變帶信息[J].地質與勘探.2000,36(5):62-65.

[6]劉成,王丹麗,李笑梅.用混合像元線性模型提取中等植被覆蓋區(qū)的粘土蝕變信息[J].遙感技術與應用.2003,18(2):95-98.

[7]鄒林,楊自安,朱谷昌，等.多光譜遙感蝕變信息提取新方法研究[J].地質與勘探. 2006, 42(6): 71-76.

[8]薛云,戴塔根,鄒艷紅，等. 基于主成分分析的SVM礦化信息提取研究——以青海黃南州阿哇地區(qū)為例[J]. 遙感信息.2007,06:32-35+44+104.

[9]張玉君,楊建民,姚佛軍. 多光譜遙感技術預測礦產資源的潛能——以蒙古國歐玉陶勒蓋銅金礦床為例[J]. 地學前緣.2007,05:63-70.

[10]陳建平,王倩,董慶吉，等. 青海沱沱河地區(qū)遙感蝕變信息提取[J]. 地球科學(中國地質大學學報).2009,02:314-318.

[11]張廷斌，等，東昆侖五龍溝金礦集中區(qū)遙感地質信息提取與找礦預測[J].地質與勘探.2009,04:444-449.

[12]張玉君,曾朝銘. 西藏岡底斯地區(qū)斑巖銅礦識別的最佳多光譜遙感異常判據(jù)研究[J]. 礦床地質.2012,04:671-698.

[13]張玉君,曾朝銘,陳薇. ETM+(TM)蝕變遙感異常提取方法研究與應用——方法選擇和技術流程[J]. 國土資源遙感.2003,02:44-49+78.

[14]劉曉玲,孫巖,陳建平，等. 基于GIS的沱沱河地區(qū)銅礦綜合信息成礦預測[J]. 物探化探計算技術.2010,04:408-416+338.

[中圖分類號]P617

[文獻標識碼]B

[文章編號]1004-1184(2016)02-0218-03

[作者簡介]趙棟(1991-)，男，山西晉中人，在讀碩士研究生，主攻方向：構造地質、遙感地質。

[收稿日期]2015-12-02