內(nèi)蒙古東烏旗遙感構(gòu)造和蝕變信息提取與找礦預(yù)測

錢建平，張 淵，趙小星，趙少杰，李承禮

(桂林理工大學(xué)地球科學(xué)學(xué)院，桂林 541004)

0 引言

近年來，內(nèi)蒙古東部東烏旗一帶Pb-Zn多金屬找礦在一些點(diǎn)上取得了重要成果。如何實(shí)現(xiàn)區(qū)域面上找礦的新突破，很多研究者試圖從區(qū)域遙感地質(zhì)找礦研究入手。遙感是一種經(jīng)濟(jì)快捷、實(shí)用有效的區(qū)域找礦手段。金屬礦床常富集在構(gòu)造帶附近或斷裂交匯部位，后者以線性體或環(huán)形體等結(jié)構(gòu)信息的形式反映在遙感圖像上[1-3]。礦化程度往往與構(gòu)造和蝕變的復(fù)雜程度相關(guān)。礦床具有與圍巖不同的礦物組成和化學(xué)成分，它們能夠以特殊的波譜特征反映在遙感圖像上[4-5]。遙感技術(shù)應(yīng)用于礦產(chǎn)勘查主要是運(yùn)用各種圖像處理技術(shù)增強(qiáng)和提取與成礦有關(guān)的構(gòu)造信息和蝕變信息，以達(dá)到識別礦床的目的[6-7]。另一方面，地球化學(xué)異常和遙感蝕變異常在本質(zhì)上是成礦作用在不同物質(zhì)層次上的表現(xiàn)形式，兩者具有一定的相關(guān)性。因此，遙感、地質(zhì)、地球化學(xué)等多種地學(xué)信息綜合分析已成為一種地質(zhì)找礦的重要手段[8-9]。

前人在本研究區(qū)曾做過部分遙感地質(zhì)工作，但主要局限于遙感線性構(gòu)造解譯方面[10]，尚未開展過區(qū)域遙感蝕變異常信息提取，特別是未開展過多元信息綜合分析和找礦研究工作。本文立足于區(qū)域成礦地質(zhì)背景分析，運(yùn)用Landsat7 ETM+數(shù)據(jù)，對研究區(qū)線性構(gòu)造重新進(jìn)行解譯，并進(jìn)行分形統(tǒng)計(jì);采用主成分分析和閾值分割方法提取遙感蝕變異常信息;總結(jié)線性構(gòu)造和遙感蝕變信息在區(qū)域空間上的分布規(guī)律，探討線性構(gòu)造集中區(qū)和遙感蝕變異常與區(qū)域巖漿巖、地球化學(xué)異常和Pb-Zn-Ag多金屬礦床(點(diǎn))在空間上的相互聯(lián)系，在此基礎(chǔ)上進(jìn)行區(qū)域成礦遠(yuǎn)景分析和找礦預(yù)測。

1 研究區(qū)概況與數(shù)據(jù)源

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

研究區(qū)位于內(nèi)蒙古自治區(qū)中東部、東烏珠穆沁旗額仁高畢蘇木至滿都胡寶格拉蘇木一帶，海拔高度一般在900～1 100 m之間，地形切割較淺，地勢相對平緩，為典型的草原覆蓋區(qū)。區(qū)內(nèi)水系極不發(fā)育，僅在丘陵低洼處散布有少量面積大小不等的濕地。區(qū)域大地構(gòu)造主要位于內(nèi)蒙古—興安嶺地槽褶皺系的東烏旗褶皺束內(nèi)。區(qū)域內(nèi)出露地層主要有古生界和中生界(圖1)，二者為角度不整合接觸。基底地層由奧陶系、志留系、泥盆系和石炭系組成，其上為侏羅系和白堊系不整合覆蓋。其中，中上泥盆統(tǒng)和中石炭統(tǒng)碎屑巖與成礦關(guān)系密切。

圖1 內(nèi)蒙古東烏旗區(qū)域地質(zhì)圖(資料來源:1∶20萬區(qū)域地質(zhì)礦產(chǎn)圖)Fig.1 Regional geological map of Dongwu Banner，Inner Mongolia Autonomous Region

區(qū)域構(gòu)造線主要呈NE和NW向，少量沿SN和EW向展布，褶皺和斷裂均很發(fā)育，形成格狀構(gòu)造。代表性斷裂有巴潤沙巴爾—朝不楞北大斷裂(F1)、朝不楞西—烏拉蓋斷裂(F2)和白云呼布爾—滿都寶力格大斷裂(F3)。

區(qū)內(nèi)巖漿巖廣泛分布，主要為中酸性、酸性，少量為堿性、基性侵入巖。海西期巖漿巖主要有花崗巖和閃長巖，燕山期巖漿巖主要有花崗巖、花崗斑巖、石英斑巖、正長斑巖、二長斑巖和閃長玢巖。

研究區(qū)位于內(nèi)蒙古—興安嶺晚古生代—中生代Cu，Pb，Zn，Au，Ag，Sn，Cr(Mo)成礦區(qū)，錫林浩特—東烏旗多金屬成礦帶的東段?？傮w來看，區(qū)域構(gòu)造和巖漿活動頻繁，多期次的中酸性巖漿活動為多金屬成礦提供了物質(zhì)來源，中、上泥盆統(tǒng)和石炭系碎屑巖地層是有利的賦礦地層，NE向斷裂及其與NW向斷裂交匯部位是有利的賦礦部位，區(qū)域成礦條件十分優(yōu)越。區(qū)內(nèi)典型礦床有朝不楞Fe多金屬礦、吉林寶力格Ag礦、阿爾哈達(dá)Pb-Zn-Ag礦、達(dá)賽脫P(yáng)b-Zn-Ag礦和額爾登陶勒蓋Cu礦等礦床。主要圍巖蝕變有矽卡巖化、褐鐵礦化、黃鐵礦化、綠泥石化、絹云母化、高嶺土化、螢石化、硅化、次生石英巖化和碳酸鹽化等(表1)。

表1 東烏旗主要金屬礦床(點(diǎn))地質(zhì)特征Tab.1 Geological features of polymetallic ore deposits in Dongwu Banner

續(xù)表

區(qū)內(nèi)多金屬元素地球化學(xué)異常發(fā)育，主要為Cu，Pb，Zn，其次為 Cr。異常區(qū)呈點(diǎn)陣式分布，受 NE和NW向斷裂控制。具體異常主要圍繞燕山期巖體分布，與Cu，Pb，Zn礦床關(guān)系密切。一般出現(xiàn) Cu，Pb，Zn多元素綜合異常，異常強(qiáng)度高、規(guī)模大者，與成礦關(guān)系密切;單元素異常，異常強(qiáng)度小、規(guī)模小者，與成礦關(guān)系較差。

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

研究區(qū)涉及2001年10月24日獲取的1景ETM+圖像及2000年9月24日獲取的1景ETM+圖像。圖像質(zhì)量好，干擾信息較少，適用于地質(zhì)構(gòu)造和蝕變信息提取。

利用ENVI的幾何糾正模塊，選擇圖像對圖像的幾何糾正方法，以區(qū)域地質(zhì)圖為基準(zhǔn)，在采集地面控制點(diǎn)基礎(chǔ)上，選用德洛內(nèi)三角測量(DE los inside triangulation)旋轉(zhuǎn)、縮放和平移(rotation，scaling and translation，RST)糾正和雙線性內(nèi)插重采樣方法，對ETM+圖像進(jìn)行了幾何糾正。為了使彩色合成圖像能最大限度地反映地表信息，通過波段相關(guān)性分析，選擇相關(guān)性較小的ETM+7(R)，4(G)，2(B)波段組合進(jìn)行假彩色合成，得到的假彩色合成圖像色彩飽和，地表信息特征明顯(圖2)。

圖2 內(nèi)蒙古東烏旗ETM+7(R)，4(G)，2(B)假彩色合成圖像Fig.2False color image composed of ETM+7(R)，4(G)，2(B)for Dongwu Banner，Inner Mongolia Autonomous Region

2 線性構(gòu)造解譯與分形研究

2.1 線性構(gòu)造增強(qiáng)處理

對基礎(chǔ)圖像(圖2)進(jìn)行高通卷積濾波增強(qiáng)，以突出線性影像特征;對ETM+1—5和7波段進(jìn)行主成分分析，在得到的主成分中，PC1主要表現(xiàn)地形、地貌信息;PC2突出表現(xiàn)構(gòu)造信息，對PC2進(jìn)行3像元×3像元的低通卷積濾波，得到PC2′;ETM+4反映了與水有關(guān)的地質(zhì)構(gòu)造、地貌和土壤、巖石類型等信息;通過 PC1(R)，ETM+4(G)，PC2′(B)假彩色合成，增強(qiáng)線性構(gòu)造與其兩側(cè)影像間的亮度反差，突出主干構(gòu)造。

2.2 線性構(gòu)造解譯

從圖2和圖3可看出，研究區(qū)內(nèi)線性構(gòu)造十分發(fā)育。線性構(gòu)造的展布方位主要為NE和NW向，次為NNW向;NW向線性構(gòu)造較為密集斷續(xù)，NE向線性構(gòu)造較為稀疏。其中，NE向線性構(gòu)造為成礦有利構(gòu)造;已知礦床(點(diǎn))主要分布于NE向與NW向線性構(gòu)造交匯部位。

圖3 內(nèi)蒙古東烏旗線性構(gòu)造解譯圖Fig.3 Interpretation map of linear structures of Dongwu Banner，Inner Mongolia Autonomous Region

2.3 線性構(gòu)造分形研究

線性構(gòu)造分維值的大小反映了線性構(gòu)造展布結(jié)構(gòu)的復(fù)雜程度、構(gòu)造的發(fā)育程度和構(gòu)造活動的強(qiáng)烈程度。一般認(rèn)為分維值愈大，構(gòu)造的空間愈復(fù)雜[11]，構(gòu)造活動性愈強(qiáng)，愈有利于成礦元素的活化和成礦流體的運(yùn)移和聚集、形成具有較大規(guī)模的礦床。

2.3.1 全區(qū)線性構(gòu)造總體分形特征

目前線性構(gòu)造分維統(tǒng)計(jì)的方法很多，本文采用計(jì)盒維數(shù)法，對研究區(qū)全區(qū)遙感線性構(gòu)造進(jìn)行分形統(tǒng)計(jì)，所采用的標(biāo)度為66～4.125 km。經(jīng)統(tǒng)計(jì)計(jì)算，求得本區(qū)線性構(gòu)造的分維值D=1.737 7，相關(guān)系數(shù)r=0.999 2。當(dāng)顯著性水平取0.05時(shí)，相關(guān)系數(shù)的臨界值為0.878。研究區(qū)的相關(guān)系數(shù)遠(yuǎn)遠(yuǎn)大于臨界值，表明研究區(qū)線性構(gòu)造具有分形特征，其分形結(jié)構(gòu)具有很好的統(tǒng)計(jì)自相似性。

將研究區(qū)的線性構(gòu)造分維值與其他地區(qū)進(jìn)行對比(表2)可知，研究區(qū)線性構(gòu)造較復(fù)雜，活動性較強(qiáng)。

表2 研究區(qū)與其他地區(qū)線性構(gòu)造分形特征比較Tab.2 Comparison of fractal characteristics between the study area and other areas

2.3.2 不同區(qū)段線性體分形特征

以圖3為底圖，將其劃分為64個(gè)邊長為16.5 km的正方形區(qū)域，對每個(gè)正方形區(qū)域采用8.25～0.515 625的5個(gè)標(biāo)度進(jìn)行研究。經(jīng)統(tǒng)計(jì)計(jì)算，求得各小正方形區(qū)域線性構(gòu)造的分維值0.831 6≤D≤1.485 1，相關(guān)系數(shù) 0.951 8≤r≤0.998 4。故每個(gè)小正方形區(qū)域的相關(guān)系數(shù)遠(yuǎn)遠(yuǎn)大于臨界值(0.878)，說明其分形結(jié)構(gòu)具有極好的相似性。利用Surfer軟件繪制出研究區(qū)線性構(gòu)造分維等值線圖(圖4)。

圖4 內(nèi)蒙古東烏旗線性構(gòu)造分維等值線圖Fig.4 Counter map of fractal dimension value of linear structures of Dongwu Banner，Inner Mongolia Autonomous Region

由圖4可知，研究區(qū)線性構(gòu)造分維值與其發(fā)育數(shù)量緊密相關(guān)，各正方形塊子區(qū)內(nèi)的線性構(gòu)造數(shù)量越多、密度越大，分維值也就越大，反之越小。

將已知礦床(點(diǎn))疊置在圖4上，發(fā)現(xiàn)已知礦床(點(diǎn))均分布在高分維值區(qū)或附近。由此可見，線性構(gòu)造分維值對預(yù)測礦點(diǎn)有一定的參考價(jià)值。因此，從分形理論的角度來說，研究區(qū)中分維值D＞1.25的強(qiáng)構(gòu)造變形域是成礦的優(yōu)勢場所。

3 遙感蝕變信息提取

3.1 蝕變礦物反射光譜特征

圖5 幾種常見蝕變巖的反射波譜曲線Fig.5 Several common reflective curves of alteration rocks

研究區(qū)內(nèi)礦化蝕變發(fā)育，幾種常見蝕變巖的反射波譜曲線如圖5所示。與金屬礦化有關(guān)的中低溫?zé)嵋何g變類型主要有黃鐵礦化、綠簾石化、螢石化、綠泥石化、硅化和粘土礦化等，這些近礦圍巖蝕變礦物分別富含 Fe3+和 Fe2+離子、水(H2O)、羥基(OH-)或碳酸根(CO2-3)等基團(tuán)。上述結(jié)構(gòu)離子的電子振動過程，會使富含這些離子或基團(tuán)的礦物產(chǎn)生特征光譜[15]。

3.2 蝕變信息提取方法

提取礦化蝕變信息的過程是對圖像中所有像元信息統(tǒng)計(jì)歸類分析的過程。研究區(qū)受植被、水系和陰影等諸多因素的影響，所以在蝕變信息提取前先要去除這些干擾因素，以使受干擾的像元不參與統(tǒng)計(jì)分析。利用ENVI軟件中的掩模運(yùn)算(Masking)功能對干擾信息進(jìn)行掩模處理。

對于上述去除干擾后的圖像，采用主成分分析方法提取蝕變異常信息。提取鐵染信息用ETM+1，3，4，5波段做主成分分析，異常主成分PC4具有的特點(diǎn)是ETM+1和ETM+4的貢獻(xiàn)系數(shù)與ETM+3的貢獻(xiàn)系數(shù)符號相反(表3)。

表 3 ETM+1，3，4，5 波段的 PCA 統(tǒng)計(jì)Tab.3 PCA statistics of ETM+1，3，4，5

從表3可以看出，在PC4中ETM+3與ETM+1，ETM+4具有相反的貢獻(xiàn)值，且ETM+4具有強(qiáng)負(fù)載值(表現(xiàn)為正值)，故PC4反映了鐵染信息，即

提取羥基信息用 ETM+1，4，5，7波段做主成分分析，異常主成分PC4具有的特點(diǎn)是ETM+4，ETM+7的貢獻(xiàn)系數(shù)與ETM+5的貢獻(xiàn)系數(shù)符號相反(表4)。

表 4 ETM+1，4，5，7 波段的 PCA 統(tǒng)計(jì)Tab.4 PCA statistics of ETM+1，4，5，7

從表4可以看出，在PC4中ETM+5與ETM+7具有相反的貢獻(xiàn)值，且ETM+5具有強(qiáng)負(fù)載值(表現(xiàn)為正值)，故PC4反映了羥基信息，即

根據(jù)上述判斷準(zhǔn)則，鐵染異常應(yīng)該包含在ETM+1，3，4，5 主成分分析的 PC4 中，羥基蝕變應(yīng)該包含在ETM+1，4，5，7主成分分析的PC4中。根據(jù)概率密度分布曲線的數(shù)學(xué)含意，可以把統(tǒng)計(jì)均值理解為主成分分析結(jié)果的區(qū)域背景值，利用Y=X+Kσ來劃分異常等級，其中，Y為異常等級;X為均值;σ為標(biāo)準(zhǔn)離差;K為常數(shù)，對于羥基異常K一般取2～3，對于鐵染異常 K 一般取1.5～2.5。據(jù)此，對所包含主要信息的主成分進(jìn)行閾值分割，得到分級異常信息圖像，并圈出7個(gè)遙感蝕變異常區(qū)(圖6)。

圖6 內(nèi)蒙古東烏旗遙感蝕變異常信息圖Fig.6 Abnomaly map of remote sensing alteration of Dongwu Banner，Inner Mongolia Autonomous Region

4 綜合信息找礦預(yù)測

4.1 多種地學(xué)信息綜合分析

將遙感蝕變異常、化探異常、礦床(點(diǎn))與巖漿巖分布等多種地學(xué)信息疊加(圖7)，可以看出如下規(guī)律:

1)與區(qū)域的格狀構(gòu)造一致，巖漿巖體、蝕變異常區(qū)呈明顯的點(diǎn)陣式排列。

圖7 內(nèi)蒙古東烏旗找礦預(yù)測圖Fig.7 Prediction map of mineral prospection of Dongwu Banner，Inner Mongolia Autonomous Region

2)遙感蝕變異常、化探異常位于巖漿巖體之上或其附近，已知礦床(點(diǎn))大多位于蝕變異常區(qū)內(nèi)。這表明，作為成礦物質(zhì)運(yùn)移通道和驅(qū)動力，區(qū)域構(gòu)造和巖漿活動是本區(qū)成礦的基本因素。巖體周邊熱能量和物質(zhì)交換最為活躍，導(dǎo)致遙感蝕變異常、化探異常最為顯著，它們是成礦物質(zhì)活化—遷移—再分布作用在不同物質(zhì)層次上的表現(xiàn)形式，因此二者均可作為找礦預(yù)測依據(jù)。

3)由NE向SW，鐵染異常強(qiáng)度漸減弱，羥基異常漸增強(qiáng)，相應(yīng)地礦床(點(diǎn))增加;且蝕變類型與礦床類型具有一定關(guān)聯(lián)性，即北段分布鐵或鐵多金屬礦床較多，南段分布多金屬礦床較多(例外的是，阿爾哈達(dá)多金屬礦床位于北段Ⅱ號預(yù)測區(qū)，但仔細(xì)分析可以看出，事實(shí)上它亦位于Ⅱ號預(yù)測區(qū)局部羥基異常中)。

4.2 找礦預(yù)測

根據(jù)區(qū)內(nèi)線性構(gòu)造分維值、遙感蝕變信息、巖體分布特征、區(qū)域地球化學(xué)異常和已知礦床(點(diǎn))分布等多種地學(xué)信息的綜合分析，探討區(qū)域找礦前景，劃分了三級找礦遠(yuǎn)景區(qū)(圖7)，具體如下:

一級找礦遠(yuǎn)景區(qū)蝕變強(qiáng)度強(qiáng)或較強(qiáng)，D＞1.25，具有多元素化探異常，區(qū)內(nèi)有已知礦床(點(diǎn))分布。

二級找礦遠(yuǎn)景區(qū)蝕變強(qiáng)度較強(qiáng)，線性構(gòu)造分維值中等，具有多元素化探異常，區(qū)內(nèi)有已知礦床(點(diǎn))分布。

三級找礦遠(yuǎn)景區(qū)蝕變強(qiáng)度較強(qiáng)或較弱，線性構(gòu)造分維值中等或較低，具有多或單元素化探異常，目前未發(fā)現(xiàn)礦點(diǎn)。

按上述原則，將研究區(qū)劃分為4個(gè)一級找礦遠(yuǎn)景區(qū)(即Ⅰ，Ⅱ，Ⅴ，Ⅵ號找礦遠(yuǎn)景區(qū))、1個(gè)二級找礦遠(yuǎn)景區(qū)(即Ⅳ號找礦遠(yuǎn)景區(qū))和2個(gè)三級找礦遠(yuǎn)景區(qū)(即Ⅲ，Ⅶ號找礦遠(yuǎn)景區(qū))(表5)。

表5 內(nèi)蒙古東烏旗各找礦預(yù)測區(qū)特征Tab.5 Characteristics of mineral prediction areas of Dongwu Banner，Inner Mongolia Autonomous Region

5 結(jié)論

1)區(qū)域遙感線性構(gòu)造解譯表明，研究區(qū)內(nèi)線性構(gòu)造十分發(fā)育，按線性構(gòu)造的方位主要為NE和NW向，其次為SN和EW向，構(gòu)成格狀構(gòu)造系統(tǒng)。

2)研究區(qū)內(nèi)線性構(gòu)造具有分形特征，其分形結(jié)構(gòu)具有良好的統(tǒng)計(jì)自相似性，其分維值高于中國大陸許多地區(qū)的斷裂構(gòu)造分維值。線性構(gòu)造分維值D＞1.25的強(qiáng)構(gòu)造變形域有利于成礦作用的發(fā)生，是成礦優(yōu)勢區(qū)域。

3)研究區(qū)內(nèi)具有顯著的鐵染和羥基遙感蝕變異常，已知礦床(點(diǎn))大多位于遙感蝕變異常區(qū)內(nèi)。與區(qū)域的格狀構(gòu)造一致，巖漿巖體、遙感蝕變異常區(qū)呈明顯的點(diǎn)陣式排列。遙感蝕變異常和化探異常位于巖漿巖體之上或其附近。作為成礦物質(zhì)運(yùn)移通道和驅(qū)動力，區(qū)域構(gòu)造和巖漿活動是本區(qū)成礦的基本因素。

4)研究區(qū)內(nèi)由NE向SW，鐵染異常強(qiáng)度漸減弱，羥基異常漸增強(qiáng)，相應(yīng)的礦床(點(diǎn))亦增多;蝕變類型與礦床類型具有一定關(guān)聯(lián)性，即北段分布鐵或鐵多金屬礦床較多，南段分布多金屬礦床較多。

5)通過遙感信息和已知礦點(diǎn)對比研究，證實(shí)了遙感解譯線性構(gòu)造和遙感蝕變異常信息可以作為找礦預(yù)測的標(biāo)志。綜合分析線性構(gòu)造、遙感蝕變異常信息、地球化學(xué)異常、巖漿巖和礦床(點(diǎn))分布等多種地學(xué)信息，確定了三級找礦遠(yuǎn)景區(qū)(其中Ⅰ，Ⅱ，Ⅴ，Ⅵ號為一級遠(yuǎn)景區(qū);Ⅳ號為二級遠(yuǎn)景區(qū);Ⅲ，Ⅶ號為三級遠(yuǎn)景區(qū))，可為研究區(qū)的地質(zhì)找礦工作提供線索。

[1]況 忠，龍勝清，曾禹人，等.黔西南地區(qū)遙感構(gòu)造與金礦的關(guān)系及找礦預(yù)測[J].國土資源遙感，2012，24(1):160-165.Kuang Z，Long S Q，Zeng Y R，et al.The relationship between remote sensing structures and gold deposits and ore-prospecting prognosis in Southwest Guizhou[J].Remote Sensing for Land and Resources，2012，24(1):160-165.

[2]趙少杰，錢建平，陳宏毅.遙感線性構(gòu)造分形統(tǒng)計(jì)和蝕變信息提取在桂東地區(qū)金鉛鋅錫多金屬成礦預(yù)測中的應(yīng)用[J].大地構(gòu)造與成礦學(xué)，2011，35(3):364-371.Zhao S J，Qian J P，Chen H Y.Application of fractal statistics of linear structure and alteration information extraction of remote sensing on the Au，Pb，Zn，Sn polymetallic minerogenetic prognosis in eastern Guangxi[J].Geotectonica et Metallogenia，2011，35(3):364-371.

[3]余 勇，錢建平，袁愛平.高龍金礦區(qū)高分辨率遙感線性構(gòu)造分形特征及綜合成礦預(yù)測[J].礦產(chǎn)與地質(zhì)，2005，19(2):209-211.Yu Y，Qian J P，Yuan A P.Fractal characteristics of high resolution remote sensing lineations in Gaolong gold deposit and integrated ore deposit prediction[J].Mineral Resources and Geology，2005，19(2):209-211.

[4]Loughlin W P.Principal component analysis for alteration mapping[J].Photogrammetric Engineering and Remote Sensing，1991，57(9):163-1169.

[5]鄭桂香，池天河，藺啟忠.分形在巖性分類及蝕變信息提取中的應(yīng)用[J].國土資源遙感，2012，24(2):110-115.Zheng G X，Chi T H，Lin Q Z.Fractal application in lithological classification and alteration extraction[J].Remote Sensing for Land and Resources，2012，24(2):110-115.

[6]Rajesh H M.Mapping proterozoic unconformity-related uranium deposits in the Rockhole area，northern Territory，Australia using Landsat ETM+[J].Ore Geology Reviews，2008，33(3/4):382-396.

[7]鄧吉秋，謝 楊，張寶一，等.ETM+圖像錳礦化蝕變信息提取與找礦預(yù)測[J].國土資源遙感，2011，23(1):102-105.Deng J Q，Xie Y，Zhang B Y，et al.The extraction of the manganese mineralization alteration information from the ETM+image and ore prognosis[J].Remote Sensing for Land and Resources，2011，23(1):102-105.

[8]劉文蘭，張 微.遙感構(gòu)造蝕變異常信息提取及找礦預(yù)測——以老撾為例[J].國土資源遙感，2012，24(2):68-74.Liu W L，Zhang W.Remote sensing structural alteration information extraction and ore prognosis:A case study of Laos[J].Remote Sensing for Land and Resources，2012，24(2):68-74.

[9]錢建平，伍貴華，陳宏毅.現(xiàn)代遙感技術(shù)在地質(zhì)找礦中的應(yīng)用[J].地質(zhì)找礦論叢，2012，27(3):355-360.Qian J P，Wu G H，Chen H Y.The application of modern remote sensing technology to geology and ore exploration[J].Contributions to Geology and Mineral Resources Research，2012，27(3):335-360.

[10]李繼宏，張萬昌.內(nèi)蒙古東烏旗草原覆蓋區(qū)遙感構(gòu)造信息提取[J].地質(zhì)找礦論叢，2005，20(2):137-141.Li J H，Zhang W C.Tectonic information extract in the Dong Ujimqin county grassland district of Inner Mongolia[J].Contributions to Geology and Mineral Resources Research，2005，20(2):137-141.

[11]陳建平，胡明銘，李巨初.康滇地軸中南段區(qū)域構(gòu)造格架的遙感地質(zhì)統(tǒng)計(jì)分析[J].成都理工學(xué)院學(xué)報(bào):自然科學(xué)版，1999，26(1):78-81.Chen J P，Hu M M，Li J C.Statistical analysis for regional structure by remote sensing geology in Kang-Dian Axis of the Earth[J].Journal of Chengdu University of Technology:Science and Technology Edition，1999，26(1):78-81.

[12]朱曉華.中國主要地貌與地質(zhì)災(zāi)害的空間分維及其關(guān)系研究[D].南京:南京師范大學(xué)，2002.Zhu X H.The study on the relationship between landform or geological hazard and spatial fractal dimension of China[D].Nanjing:Nanjing Normal University，2002.

[13]Hirata T.Fractal dimension of fault systems in Japan，F(xiàn)ractal structure in rock fracture geometry at various scales[J].Pure and Applied Geophysics，1989，131(12):157-170.

[14]孔凡臣，丁國瑜.線性構(gòu)造分維值的含義[J].地震，1991(5):33-37.Kong F C，Ding G Y.The implications of the fractal dimension values of lineaments[J].Earthquake，1991(5):33-37.

[15]張玉君，楊建民，陳 薇.ETM+(TM)蝕變遙感異常提取方法研究與應(yīng)用——地質(zhì)依據(jù)和波譜前提[J].國土資源遙感，2002，14(4):30-36.Zhang Y J，Yang J M，Chen W.A study of the method for extraction of alteration anomalies from the ETM+(TM)data and its application:Geologic basis and spectral precondition[J].Remote Sensing for Land and Resources，2002，14(4):30-36.