基于遙感與分散流異常綜合找礦預(yù)測(cè)研究

于建華，彭中勤

YUJian-Hua1,PENGZhong-Qin2

(1.天津華北地質(zhì)勘查局,天津300170，2.武漢地質(zhì)調(diào)查中心，武漢430205)

（1.Tianjin North China Geological Exploration Bureau,Tianjin 300170；2.Wuhan Center of China Geological Survey,Wuhan 430205,China）

1 引言

隨著遙感技術(shù)應(yīng)用面不斷擴(kuò)大，利用遙感信息和地球化學(xué)異常找礦綜合分析，進(jìn)而圈定出找礦靶區(qū)，發(fā)展成為一種找礦預(yù)測(cè)新技術(shù)方法[1-5]。該方法基于成礦元素的地球化學(xué)異常場(chǎng)與遙感影像特征之間存在某種確定的專屬聯(lián)系的遙感地球化學(xué)理論。對(duì)于以構(gòu)造控礦為主的成礦區(qū)域，遙感影像資料最能從宏觀上揭示構(gòu)造控制的成礦有利部位，如作為熱液通道的各級(jí)斷裂構(gòu)造及其交匯部位、構(gòu)造容礦機(jī)構(gòu)及伴隨的礦化、蝕變帶地物波譜異常區(qū)等。研究證明，正是這些遙感找礦異常區(qū)通常與地球化學(xué)異常存在著某種相關(guān)性。因此，在找礦預(yù)測(cè)中，應(yīng)用遙感信息和地球化學(xué)異常相互配合，互為補(bǔ)充，就可以提高找礦預(yù)測(cè)的效果與效率。

目前，遙感技術(shù)與應(yīng)用已從地球資源的現(xiàn)狀調(diào)查實(shí)驗(yàn)階段進(jìn)入?yún)⑴c規(guī)劃、決策和評(píng)價(jià)的實(shí)用階段，它正在轉(zhuǎn)向?qū)嵱没蜕虡I(yè)化。空間遙感技術(shù)是21世紀(jì)科學(xué)技術(shù)發(fā)展中最激動(dòng)人心、最有希望的高科技之一。由此可見，遙感技術(shù)與地球化學(xué)信息綜合找礦研究，是將最新技術(shù)與傳統(tǒng)找礦模式的結(jié)合，這將為傳統(tǒng)的地質(zhì)找礦工作帶來(lái)深遠(yuǎn)的影響。

2 地質(zhì)背景

研究區(qū)位于廣西壯族自治區(qū)北部的龍勝縣與湖南省城步縣之間，交通極為不便，利用傳統(tǒng)方法在該區(qū)開展找礦工作，不僅投入人力物力十分巨大，周期相當(dāng)長(zhǎng)，而且效果也不明顯。利用遙感和分散流技術(shù)在該區(qū)開展找礦工作，剛好能彌補(bǔ)傳統(tǒng)方法的不足，更能發(fā)揮其巨大的優(yōu)勢(shì)，取得良好的效果。研究區(qū)南北長(zhǎng)約7.7 km，東西寬約6.6 km。

龍勝地區(qū)大地構(gòu)造位置隸屬于揚(yáng)子準(zhǔn)地臺(tái)南部邊緣，位于華南板塊二級(jí)構(gòu)造單元揚(yáng)子陸塊南緣與華南褶皺帶的交接部位。加里東造山運(yùn)動(dòng)期，該區(qū)處于華南褶皺帶變形的前鋒，構(gòu)造變形強(qiáng)烈，表現(xiàn)為緊密的線性褶皺和推覆構(gòu)造[6]。一條NNE向區(qū)域深斷裂（泗水-五團(tuán)深斷裂)貫穿本區(qū)。本區(qū)區(qū)域深斷裂具有發(fā)生時(shí)間早(晚元古代)、切割深度大(硅鋁-硅鎂層)、活動(dòng)時(shí)間長(zhǎng)(雪峰期-印支燕山期)、普遍發(fā)生脆-韌性剪切變形（如基性-超基性、長(zhǎng)英質(zhì)糜棱巖帶、劈（片）理化帶普遍）等特征。區(qū)域深斷裂及其次級(jí)NNE向、近NS向區(qū)域斷裂及更次級(jí)NE、NW、NEE向剪切斷裂構(gòu)成統(tǒng)一的“面狀”透入性構(gòu)造帶。褶皺構(gòu)造主要為NNE向緊閉背斜，如平等背斜、龍勝背斜、馬海背斜，背斜核部為丹洲群，兩翼為震旦系、寒武系。

桂北龍勝里市地區(qū)出露地層[7]主要有四堡群、丹洲群、震旦系、寒武系。其中，下-中元古宇四堡群為該區(qū)出露的最老地層，構(gòu)成下部褶皺基底，上元古宇-下古生界地槽型沉積，構(gòu)成上部褶皺基底[8]，其中主要以丹洲群、震旦系為主，約占全區(qū)的80%左右[9]。

3 遙感數(shù)據(jù)的選擇與信息提取

3.1 遙感數(shù)據(jù)的選擇

根據(jù)工作區(qū)遙感圖像資源狀況、地區(qū)自然氣候和植被自然生長(zhǎng)情況，衛(wèi)星遙感圖像月份應(yīng)盡量選擇為晚秋（9~11月上旬）。為了保證圖像解像力和解譯效果，購(gòu)買圖像時(shí)，要盡量選擇無(wú)云或云量較少的圖像。因此，在本次研究中選取了美國(guó)陸地資源衛(wèi)星Landsat-7 ETM+（20001106）數(shù)據(jù)和法國(guó)SPOT-5（20051002）衛(wèi)星數(shù)據(jù)。

本文主要利用Landsat-7 ETM+數(shù)據(jù)的B7、B4、B3的假彩色合成影像來(lái)解譯出研究區(qū)內(nèi)一定數(shù)量的區(qū)域性構(gòu)造信息，宏觀的角度為SPOT-5高分辨率的構(gòu)造信息提取提供有益參考，并利用多光譜波段，采集出一定數(shù)量與我們采集的分散流數(shù)據(jù)點(diǎn)對(duì)應(yīng)的灰度值，制成波譜曲線圖，找出不同波段灰度值與分散流異常之間的聯(lián)系。基于Landsat-7 ETM+遙感數(shù)據(jù)，利用采樣點(diǎn)的最佳示礦波段確定法，開展了基于分散流采樣點(diǎn)的最佳示礦波段選擇研究，利用遙感影像數(shù)據(jù)對(duì)已知分散流采樣點(diǎn)進(jìn)行Z剖面波譜分析，確定出最佳示礦波段。

由于研究區(qū)缺少中、大比例尺的基礎(chǔ)地質(zhì)資料的情況下，利用SPOT-5的高分辨率的特征，解譯出研究區(qū)內(nèi)一定數(shù)量的構(gòu)造信息，為研究區(qū)找礦預(yù)測(cè)提供一定的依據(jù)；并在原有1∶50萬(wàn)地質(zhì)圖的基礎(chǔ)上，運(yùn)用遙感圖像處理方法，圈定出本研究區(qū)各地層之間的界線。最后通過(guò)分散流數(shù)據(jù)分析結(jié)果、Landsat-7 ETM+解譯信息、SPOT-5圖像處理結(jié)果等進(jìn)行綜合疊加、對(duì)比分析，為找礦預(yù)測(cè)研究提供有力的證據(jù)和直觀、精確的定位。

3.2 遙感線性構(gòu)造提取

線性體是指遙感影像上任何天然的線或線性排列，遙感圖像中解譯出的線性構(gòu)造多數(shù)是構(gòu)造要素的反映，并主要與斷裂構(gòu)造（包括節(jié)理、斷層、斷層帶）有關(guān)[10-12]。在一般情況下，線性體的方位多具近于正交的特點(diǎn)，并多具優(yōu)勢(shì)方位，其優(yōu)勢(shì)方位反映了區(qū)域構(gòu)造基本輪廓，而其方位偏差則反映了局部構(gòu)造異常，規(guī)模較大且延續(xù)性強(qiáng)的線性影像異常是深層構(gòu)造在地表的直接表現(xiàn)，因此，通過(guò)遙感影像所提取的線性構(gòu)造大多數(shù)與構(gòu)造要素存在密切聯(lián)系，通過(guò)對(duì)線性影像的統(tǒng)計(jì)分析，有助于分析研究各類構(gòu)造的成生和空間展布規(guī)律，為區(qū)域構(gòu)造及其找礦研究提供重要依據(jù)。

項(xiàng)目在綜合了對(duì)Landsat-7-ETM+和SPOT-5影像數(shù)據(jù)的解譯基礎(chǔ)上，解譯了141條線性構(gòu)造和6個(gè)弧形構(gòu)造，其中有幾條較大的北東向的線性構(gòu)造：L1、L4、L5穿過(guò)整個(gè)研究區(qū)，伴隨著這幾條大的線性構(gòu)造周圍，發(fā)育一系列的次一級(jí)的北西和北東向線性構(gòu)造。從線性構(gòu)造分布的特征來(lái)看，在全區(qū)范圍內(nèi)均有展布，主要方向?yàn)楸北睎|向，次一級(jí)為北西向。線性構(gòu)造密度分布在不同部位差別明顯，密度最大的區(qū)域主要分布在研究區(qū)SW 方向的東升地區(qū)及其周圍，以及大白山山系帶上，這在很大程度上反映了這兩個(gè)地區(qū)斷裂構(gòu)造的實(shí)際分布狀況。

通過(guò)對(duì)位于大白山山系以東的已知Pb和Sn礦礦點(diǎn)的查證，這個(gè)已知礦點(diǎn)的礦脈展布方向?yàn)镹EE，產(chǎn)于區(qū)域性的大白山山系斷裂構(gòu)造帶的次一級(jí)NEE斷裂構(gòu)造帶中，其形成與大白山山系發(fā)育的構(gòu)造帶有關(guān)。因此，解譯出本研究區(qū)內(nèi)的次一級(jí)、二級(jí)斷裂構(gòu)造有著相當(dāng)重要的意義。在研究區(qū)約50 km2的范圍內(nèi)，解譯出了141線性構(gòu)造和采集了287個(gè)分散流點(diǎn)，經(jīng)實(shí)地查證結(jié)果顯示，遙感解譯的這些線性構(gòu)造大部分都是斷裂構(gòu)造。

4 分散流信息

分散流實(shí)質(zhì)上是次生暈在水系沉積物中的一種特殊表現(xiàn)形式，是在表生作用下形成的一種產(chǎn)物[13]。它可以在次生暈形成的同時(shí)或其后，形成的方式主要是通過(guò)流水作用將遭受風(fēng)化的礦體及其原生暈機(jī)械破碎物和鹽類狀態(tài)的成礦元素及伴生元素進(jìn)行遷移[14]。

地球化學(xué)異常下限的確定是勘查地球化學(xué)的一個(gè)基本問(wèn)題，也是勘查地球化學(xué)應(yīng)用于礦產(chǎn)預(yù)測(cè)時(shí)的一個(gè)關(guān)鍵性環(huán)節(jié)。根據(jù)龍勝里市地區(qū)鎢錫等多金屬勘查區(qū)分散流數(shù)據(jù)分析結(jié)果，通過(guò)地球化學(xué)異常下限的確定方法，即通過(guò)統(tǒng)計(jì)勘查地球化學(xué)數(shù)據(jù)，檢驗(yàn)其是否符合正態(tài)分布或?qū)?shù)正態(tài)分布，如果不符合，則篩除部分異常數(shù)據(jù)，通常是元素的高含量數(shù)據(jù)，直到數(shù)據(jù)符合正態(tài)分布或?qū)?shù)正態(tài)分布為止。本文剔除了六種分散流數(shù)據(jù)的部分高值和低值異常后，將數(shù)據(jù)按對(duì)數(shù)間隔分組，得到數(shù)據(jù)的對(duì)數(shù)直方圖，直到研究區(qū)六種元素分散流異常的對(duì)數(shù)符合對(duì)數(shù)正態(tài)分布，并利用公式[15]：

計(jì)算得到各元素異常信息的平均值（c）和標(biāo)準(zhǔn)離差（σ），并以平均值與1.5倍標(biāo)準(zhǔn)離差之和(表1)作為六種元素異常的下限（c+1.5σ）（我國(guó)目前使用的地球化學(xué)勘查標(biāo)準(zhǔn)中規(guī)定的異常下限是（c+σ）），圈定異常區(qū)域，并將圈得的異常區(qū)域在地理信息系統(tǒng)軟件sufer 8.0中分層保存，經(jīng)投影變換和坐標(biāo)配準(zhǔn)后，進(jìn)行層與層之間的空間疊加，從而得到七種分散流元素的異常區(qū)域空間疊加圖，W和Sn元素之間，Mo、Ag、Pb、Cu和Zn元素異常區(qū)域之間都有很大面積的重疊，顯示了很好的空間相關(guān)性。

表1 地球化學(xué)分散流元素的均值、標(biāo)準(zhǔn)離差和異常下限Table 1 The elements of the mean，standard deviation and anomaly thresholds

研究區(qū)分散流測(cè)量顯示，W 分散流值為：1.3~69.2μg/g。W異常高值出現(xiàn)在研究區(qū)的西北角上，地理位置上主要位于東升地區(qū)的西北部，在東升的南部和東部也有一小片面積的W異常出現(xiàn)。

Pb的分散流值為：0.17~9.98μg/g，Mo異常也出現(xiàn)在大白山山系南部；

Zn的分散流值為：20~245μg/g，Zn異常出露的位置大致與Mo相同。

綜合分析W、Sn、Zn、Ag、Pb、Cu和Mo七個(gè)元素，顯示出這些元素之間的高度相關(guān)性。W元素的異常高值區(qū)流值為8.8~56μg/g，分散流測(cè)量Pb異常在研究區(qū)出露面積廣，在東升的西部和東部，大白山山系的南部都有大面積的出露，是所有分散流數(shù)據(jù)分析樣品中出露面積最廣的；

Cu的分散流值為：12.1~566.2μg/g，Cu的異常出露面積也很廣，但高值區(qū)面積較小，主要的高值點(diǎn)位于大白山山系南部；

Ag的分散流值為：0.034~0.67μg/g，Ag異常主要出現(xiàn)在大白山山系南部；

Sn的分散流值為：1.2~27.1μg/g，Sn異常與W異常出露的位置大致相同。二者表現(xiàn)為很好的相關(guān)性；

Mo的分散，Sn元素也同樣出現(xiàn)幾乎重合的異常高值現(xiàn)象。顯示出了W和Sn高度相關(guān)性，而在這兩個(gè)元素的高值區(qū)，其他幾種元素除Pb以外幾乎無(wú)高于異常下限值出現(xiàn)。同樣的情況也出現(xiàn)在Mo、Ag、Zn、Cu元素的異常區(qū)。只有Pb在同時(shí)出現(xiàn)在兩種異常區(qū)中，說(shuō)明在研究區(qū)內(nèi)，元素伴生和組合存在規(guī)律性。

5 找礦預(yù)測(cè)研究

5.1 波譜信息提取與最佳示礦波段的確定

本次研究的波譜信息提取是在Landsat-7 ETM+遙感數(shù)據(jù)上進(jìn)行的。由于Landsat-7 ETM+遙感數(shù)據(jù)的分辨率為30 m，而由MapGIS中采集的數(shù)據(jù)點(diǎn)坐標(biāo)數(shù)據(jù)是精確到小數(shù)點(diǎn)后兩位。所以，在我們利用分散流數(shù)據(jù)點(diǎn)的坐標(biāo)值在遙感圖像上運(yùn)用象索定位找取對(duì)應(yīng)點(diǎn)時(shí)，只能找取與坐標(biāo)值最臨近的點(diǎn)作為波譜信息采集點(diǎn)，但是，從測(cè)量的結(jié)果表明，采取的最臨近點(diǎn)的波譜值基本能夠反應(yīng)出分散流數(shù)據(jù)點(diǎn)的波譜特性。在采集的287個(gè)分散流數(shù)據(jù)點(diǎn)中，隨機(jī)抽取了10個(gè)點(diǎn)作為本次波譜信息測(cè)定的點(diǎn)，波譜測(cè)定結(jié)果如表2所示。根據(jù)ETM+的不同波段的特征和用途，本文選擇Band-1、Band-2、Band-3、Band-4、Band-5和Band-7的波譜反射率進(jìn)行對(duì)比。

表2 Landsat-7 ETM+波譜信息測(cè)定結(jié)果Table 2 The results of the spectrum information of Landsat-7 ETM+

由表2波譜測(cè)定結(jié)果、表3的分散流值和圖1的分散流數(shù)據(jù)點(diǎn)的波譜曲線特征圖可以看出，分散流點(diǎn)3、點(diǎn)8和點(diǎn)9的Mo、Ag、Cu、Pb、Zn的光譜分析值較低，對(duì)于光譜曲線圖中，這幾個(gè)點(diǎn)的光譜曲線都表現(xiàn)為相對(duì)起伏不大，點(diǎn)2的分散流光譜分析結(jié)果中Cu的值相對(duì)較高，點(diǎn)4、點(diǎn)6、點(diǎn)7和點(diǎn)10的分散流光譜分析結(jié)果中表現(xiàn)出Pb的值較高，而恰好在這幾個(gè)點(diǎn)的位置的Band-4的灰度值也較高，相對(duì)于Band-3與Band-5出現(xiàn)了異常的高值現(xiàn)象。在點(diǎn)5的分散流光譜分析結(jié)果中表現(xiàn)出了Cu和Pb的值均較高的現(xiàn)象，該點(diǎn)處的Band-4的灰度值較Band-3和Band-5有一定的提高，但相對(duì)與點(diǎn)4、點(diǎn)6、點(diǎn)7和點(diǎn)10的提高效果反而并沒(méi)

有那么大，正是由于該點(diǎn)Cu含量的增加反而降低了該點(diǎn)的灰度值。點(diǎn)1的光譜分析值比較特殊，雖然在該點(diǎn)的光譜分析值較低，卻在Band-4中的灰度值同其他數(shù)據(jù)點(diǎn)一樣出現(xiàn)了異常的高值現(xiàn)象。

表3 龍勝里市地區(qū)部分分散流值Table 3 Some results of the spectral analysis in the point of dispersion flow data

圖1 部分分散流數(shù)據(jù)點(diǎn)的波譜曲線特征圖Fig.1 Thespectral curve in thepoint of dispersion flow data

總的來(lái)說(shuō)，分散流異常與遙感灰度值之間有著一定的聯(lián)系，根據(jù)不同波段波譜反射率對(duì)比結(jié)果來(lái)看，主要表現(xiàn)為在Band-4波段處，分散流異常高值點(diǎn)的遙感影像灰度值也相應(yīng)出現(xiàn)異常高值，Band-5和Band-2波段次之，在其他波段表現(xiàn)卻并不明顯。所以，在本特定的研究區(qū)，可認(rèn)為L(zhǎng)andsat-7 ETM+數(shù)據(jù)的Band-4波段為最佳示礦波段，Band-5和Band-2次之。

5.2 遙感波譜信息與分散流數(shù)據(jù)回歸分析

地球化學(xué)分散流異常和遙感信息雖然在其信息專業(yè)屬性上不同，但都是根據(jù)地表層的地質(zhì)體或地質(zhì)現(xiàn)象所獲得的，都是以不同形式的化學(xué)量和物理量表現(xiàn)出來(lái)的與找礦有關(guān)的信息，兩者之間必然存在某種相關(guān)性，很多研究結(jié)果也證實(shí)了這一點(diǎn)[16-21]。

為了找出研究中遙感采集的波譜信息與分散流數(shù)據(jù)之間的相關(guān)性，本研究引入了回歸分析和相關(guān)分析[22]。對(duì)任何兩個(gè)變量X和Y的一組實(shí)驗(yàn)數(shù)據(jù)（Xi，Yi），i=1,2,…，n，都可以按上述步驟配出一條直線。但是，在實(shí)際中，只有當(dāng)Y和X之間存在某種線性關(guān)系時(shí)配出的直線才有意義。

檢驗(yàn)回歸線有無(wú)意義，在數(shù)學(xué)上給出了一種輔助辦法，引進(jìn)了相關(guān)系數(shù)的量，用r表示。

相關(guān)系數(shù)由下式定義：

式中，r為相關(guān)系數(shù)，其絕對(duì)值越接近于1，X和Y的線性關(guān)系挺好。如果它近于0，就可以認(rèn)為X和Y之間沒(méi)有線性關(guān)系，這時(shí)有兩種情況，或者兩者沒(méi)有關(guān)系，或者兩者有非線性關(guān)系。算出的相關(guān)系數(shù)，可以查閱相應(yīng)相關(guān)系數(shù)表來(lái)檢驗(yàn)兩種數(shù)據(jù)相關(guān)與否（表4）。結(jié)合本次研究中的分散流波譜測(cè)量值（Bi）與遙感采集的波譜值（Wn），得出：

式中，Bi為不同波段遙感采集的波譜值；Wn為不同元素分散流波譜測(cè)量值；rBi-W為兩者的相關(guān)系數(shù)。

表4 相關(guān)系數(shù)檢驗(yàn)表Table 4 Checklists of Correlation coefficient

表4為相關(guān)系數(shù)的起碼值，求出的相關(guān)系數(shù)要大于表上的數(shù)，才能考慮用直線來(lái)描述X和Y之間的關(guān)系，本研究中，采集了10個(gè)相對(duì)應(yīng)的點(diǎn)，查對(duì)應(yīng)n-2=8的一行，相應(yīng)的數(shù)為0.632（5%）和0.765（1%）。并求出了W、Sn、Mo、Ag、Cu、Pb、Zn對(duì)應(yīng)于Landsat-7 ETM的B1、B2、B3、B4、B5、B7波段的相關(guān)系數(shù)（表5）。

表5 遙感波譜采集值與分散流數(shù)據(jù)的相關(guān)系數(shù)表Table 5 The correlation coefficients for the spectrum value and the spectral analysis data

由表5可以看出，分散流Zn的光譜測(cè)量值與Band-2、Band-4、Band-5、Band-7呈高水平正相關(guān)，與Band-1、Band-3呈低水平正相關(guān)；Pb的光譜測(cè)量值與Band-4、Band-5呈高水平正相關(guān)，與Band-2呈低水平正相關(guān)；Ag的光譜測(cè)量值與Band-4呈低水平正相關(guān)；Sn的光譜測(cè)量值與Band-1呈高水平負(fù)相關(guān)，其余各個(gè)元素與波譜值之間均未相關(guān)。

研究中，綜合相關(guān)性分析系數(shù)統(tǒng)計(jì)和回歸方程，求出了具有相關(guān)性的分散流波譜測(cè)量值（Bi）與遙感采集的波譜值（Wn）的回歸方程。

其中，呈高水平相關(guān)元素的回歸方程：

Zn=-90.77+2.82(ETM+2)

Zn=-356.17+4.14(ETM+4)

Zn=-186.09+3.52(ETM+5)

Zn=-70.68+2.81(ETM+7)

Pb=-1162.79+10.97(ETM+4)

Pb=-646.60+8.51(ETM+5)

低水平相關(guān)的回歸方程：

Zn=-43.48+2.61(ETM+1)

Zn=-36.56+2.45(ETM+3)

Pb=-364.76+6.02(ETM+2)

Sn=47.67+2.08(ETM+1)

由所求出了分散流光譜測(cè)量值（Bi）與遙感采集的波譜值（Wn）的回歸方程，得到兩者之間的相關(guān)關(guān)系，可以看出兩者的線性關(guān)系。因此，在不考慮其他因素的影響下，可以通過(guò)采集研究區(qū)遙感影像上的每個(gè)象元點(diǎn)上的的波譜值（Bi），求出全區(qū)的分散流光譜測(cè)量值（Wn）。但是，由于分散流是通過(guò)流水作用將遭受風(fēng)化的礦體及其原生暈機(jī)械破碎物和鹽類狀態(tài)的成礦元素及伴生元素進(jìn)行遷移的產(chǎn)物，在應(yīng)用回歸方程求出分散流的光譜值時(shí)，必須考慮分散流的成因、來(lái)源范圍以及地形條件等方面的影響。但是，回歸方程為遙感找礦提供了一種定量化的手段，具有十分重要的意義。

5.3 找礦預(yù)測(cè)區(qū)優(yōu)選

根據(jù)遙感所解譯的線性構(gòu)造分布和分散流地球化學(xué)測(cè)量的異常等值線圖，進(jìn)行數(shù)據(jù)綜合分析，在研究區(qū)內(nèi)圈定了兩個(gè)一級(jí)找礦預(yù)測(cè)區(qū)和兩個(gè)二級(jí)找礦預(yù)測(cè)區(qū)，獲得研究區(qū)線性構(gòu)造、分散流異常綜合找礦圖（圖版Ⅱ）。找礦預(yù)測(cè)區(qū)圈定的條件和原則主要有以下幾個(gè)點(diǎn)：

第一、同時(shí)擁有遙感解譯的線性構(gòu)造高密度區(qū)和分散流異常的高值區(qū)的位置，定義為一級(jí)找礦預(yù)測(cè)區(qū)；

第二、根據(jù)本研究區(qū)內(nèi)構(gòu)造控礦的特點(diǎn)，找礦預(yù)測(cè)區(qū)位于遙感等密度分析結(jié)果的線性構(gòu)造高密度區(qū)，而沒(méi)有出現(xiàn)分散流異常的高值，定義為二級(jí)找礦預(yù)測(cè)區(qū)；

第三、找礦預(yù)測(cè)區(qū)內(nèi)具有分散流光譜測(cè)定出的高值點(diǎn)，而沒(méi)有出現(xiàn)線性構(gòu)造的高密度區(qū)，也定義為二級(jí)找礦預(yù)測(cè)區(qū)。

一級(jí)預(yù)測(cè)區(qū)I面積為2760 m2，位于里市的西側(cè)，主要受北北東向區(qū)域斷裂影響較大，該斷裂位于研究區(qū)北西角上，具體位置在Landsat-7 ETM+影像中能很清楚地看到，該區(qū)的W和Sn異常主要產(chǎn)于區(qū)域斷裂構(gòu)造旁側(cè)的次一級(jí)斷裂構(gòu)造帶中。區(qū)內(nèi)主要發(fā)育北北東向斷裂構(gòu)造，是遙感解譯的線性構(gòu)造的高密度區(qū)。在該預(yù)測(cè)區(qū)內(nèi)，也是分散流異常的高值區(qū)，主要為W、Sn和Pb的異常，而其他金屬元素在此區(qū)中的異常表現(xiàn)不明顯。在該區(qū)所采集的土壤樣品分析結(jié)果與預(yù)測(cè)效果相吻合，也驗(yàn)證了本次研究中找礦預(yù)測(cè)的準(zhǔn)確性。

一級(jí)預(yù)測(cè)區(qū)II的面積為6320 m2，該區(qū)位于田邊東側(cè)。地形條件復(fù)雜，是本次研究區(qū)內(nèi)地勢(shì)最高的地帶。主要受北北東向構(gòu)造控礦，區(qū)內(nèi)發(fā)育北北東向、北西向線性構(gòu)造，也是遙感解譯的線性構(gòu)造高密度區(qū)。該區(qū)同時(shí)也是分散流異常的高值區(qū)，主要為Pb、Mo、Ag、Zn、Cu的異常，而W和Sn在此區(qū)的異常并不明顯。

此外，還圈出了兩個(gè)二級(jí)預(yù)測(cè)區(qū)，分別位于東升和大白山以東，這兩個(gè)預(yù)測(cè)區(qū)的面積較小，二級(jí)預(yù)測(cè)區(qū)I內(nèi)，主要發(fā)育北北東向、北西向線性構(gòu)造，分散流測(cè)量結(jié)果高于異常下限；二級(jí)預(yù)測(cè)區(qū)II內(nèi)主要發(fā)育北北東向線性構(gòu)造，少量南北向線性構(gòu)造。在該區(qū)內(nèi)我們尚未布置分散流采樣點(diǎn)，只是根據(jù)其南部的已知礦點(diǎn)的特征進(jìn)行推測(cè)的結(jié)果，有待下一步工作進(jìn)行查證。

6 結(jié)論

從遙感技術(shù)與地球化學(xué)分散流的結(jié)合來(lái)看，地球化學(xué)分散流為礦產(chǎn)勘查提供了直接信息，但其對(duì)具體目標(biāo)的指示性不明確。而記錄地物光譜特征的遙感圖像具有視域廣、直觀性強(qiáng)、綜合信息豐富的特點(diǎn)，其對(duì)地面地質(zhì)條件(地層、巖性、構(gòu)造等)、地形地貌條件等信息均可直接提取，是查明地球化學(xué)分散流異常形成因素的有利手段。因此，結(jié)合成礦地質(zhì)背景，對(duì)這兩種異常信息加以綜合分析并利用，將會(huì)使地表找礦信息更加專業(yè)化，從而提高地質(zhì)找礦水平和找礦效率。

通過(guò)對(duì)分散流異常和遙感信息的綜合找礦預(yù)測(cè)研究，可以得出以下幾點(diǎn)認(rèn)識(shí)：

（1）遙感構(gòu)造信息和地球化學(xué)分散流異常帶的分布同構(gòu)造體系的展布特征有很強(qiáng)的一致性，分散流異常高值較多的聚集在遙感所解譯的線性/環(huán)形構(gòu)造交切部位，證明了本研究區(qū)構(gòu)造控礦的特點(diǎn)；

（2）利用遙感影像數(shù)據(jù)對(duì)已知分散流采樣點(diǎn)進(jìn)行波譜分析，確定出最佳示礦波段為Band-4，從而揭示出Band-4最佳示礦功能，因而其可以作為本區(qū)遙感找礦基本波段資源；

（3）通過(guò)對(duì)區(qū)內(nèi)遙感高光譜與分散流數(shù)據(jù)之間的回歸分析，確定出高水平相關(guān)波段-元素對(duì)是：Zn-ETM+5、Zn-ETM+7、Zn-ETM+4、Zn-ETM+2，Pb-ETM+4、Pb-ETM+5、Sn-ETM+1。

圖版Ⅰ

圖版Ⅱ

[1]吳昀昭,田慶久,季峻峰，陳駿.遙感地球化學(xué)研究[J].地球科學(xué)進(jìn)展,2003,18（2）：228-235.

[2]吳德文，袁繼明，張遠(yuǎn)飛，朱谷昌.遙感與化探數(shù)據(jù)融合處理技術(shù)方法及應(yīng)用研究[J].國(guó)土資源遙感，2005,3（65）：44-47.

[3]蔣立軍，邢立新，梁一鴻，潘軍，梁立恒，黃競(jìng)鋮.融合化探信息的遙感異常提取[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版)，2011,41(3):932-935.

[4]楊自安，彭省臨，劉悟輝，石菲菲，和志軍，鄒林.高寒山區(qū)遙感與化探綜合找礦信息的提取[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,40(4):1121-1126.

[5]吳虹，彭中勤，賈志強(qiáng).ETM%2b多光譜遙感-化探分散流數(shù)據(jù)回歸反演模型——以桂北龍勝里市地區(qū)為例[C]//第七屆成像光譜技術(shù)與應(yīng)用研討會(huì)論文集，2010：92-99.

[6]郭令智，施央申，馬瑞士.華南大地構(gòu)造格架與地殼演化[A]//第26屆國(guó)際地質(zhì)大會(huì)國(guó)際交流學(xué)術(shù)論文集（一）[C].北京：地質(zhì)出版社,1980：109-116.

[7]廣西壯族自治區(qū)地質(zhì)礦產(chǎn)局.廣西區(qū)域地質(zhì)志[M].1985：6.

[8]郭進(jìn)義，楊忠芳.淺議地球化學(xué)的外延[J].地球科學(xué)進(jìn)展，1998,13（1）：78-80.

[9]陳毓串，毛景文.桂北地區(qū)礦床成礦系列和成礦歷史演化軌跡[M].南寧：廣西科學(xué)技術(shù)出版社.1995.

[10]孫家炳.遙感原理與應(yīng)用[M].武漢：武漢大學(xué)出版社，2003.

[11]成建勛，葉松青，丁楓，唐菊興.藏東馬牧普地區(qū)斑巖型銅金銀多金屬礦化帶地球化學(xué)異常評(píng)價(jià)[J].世界地質(zhì)，2005，24(4):334-337.

[12]Turner A K.Three-Deimension Modeling with Geoscientific Information Systems [M].Dordrecht:Kluwer Academic Publishers,NATOASINo.354,1992.

[13]幸世軍.化探分散流異常查證方法探討[J].礦產(chǎn)與地質(zhì),2001,15(Z1):480-482.

[14]Raper JF.Three dimensional GISfor the 1990's[J].ITCJournal,1994,(1):64-65.

[15]S1att R M,Thomasson M R,Romig Jr PR,Pastemack E S,Boulanger A,Anderson R N,Nelson Jr H R.Visualization technology for the oil and gas industry:Today and tomorrow[J].AAPGBulletin.1996,80(4):453-458.

[16]何延波,楊琨,候英雨.淺談遙感地球化學(xué)[J].地質(zhì)地球化學(xué),1997,4:98-103.

[17]吳昀昭,田慶久,季峻峰,陳駿.遙感地球化學(xué)研究[J].地球科學(xué)進(jìn)展，2003,18（2）：228-235.

[18]陳三明,吳虹,譚泛，張振.基于ASTER的遙感地球化學(xué)統(tǒng)計(jì)預(yù)測(cè)模型及應(yīng)用-以金川銅鎳礦床外圍找礦為例[J].桂林理工大學(xué)學(xué)報(bào)，2010,30（4）：480-488.

[19]陳勇敢，劉桂閣,路彥明,王美娟,韓先菊，張慧玉，常春郊.遙感地球化學(xué)異常信息提取研究-以青海省白日其利溝地區(qū)為例[J].礦床地質(zhì)，2010，29：35-40.

[20]張華平,劉麗,王增輝,趙西強(qiáng).遙感在多目標(biāo)地球化學(xué)調(diào)查中的應(yīng)用研究 [J].國(guó)土資源遙感,2012,95（4）：132-136.

[21]朱莉莉，洪金益，張金良，張昌俊,張思穎,段宗恩.遙感信息與化探在湘南地區(qū)找礦預(yù)測(cè)中的應(yīng)用[J].物探化探計(jì)算技術(shù)，2014,36（2）：241-246.

[22]中國(guó)科學(xué)院數(shù)學(xué)研究所統(tǒng)計(jì)組.常用數(shù)理統(tǒng)計(jì)方法[M].北京：科學(xué)出版社，1974:11，82-88.