張悅秋，居維偉，劉小胡，肖 娥，牛 琳，梅 芳

（江蘇省有色金屬華東地質(zhì)勘查局，江蘇 南京 210007）

鋯鈦?zhàn)鳛閼?zhàn)略性新興礦產(chǎn)，被廣泛應(yīng)用于航空、航天、化工、電力等領(lǐng)域，需求量不斷增大，而國內(nèi)鋯鈦資源有限，每年鋯、鈦礦進(jìn)口需求量分別達(dá)到90%和70%。為擺脫對ILUKA、TRONOX和Rio Tinto三大供應(yīng)商的依賴，保障我國未來鋯鈦供應(yīng)，國內(nèi)相關(guān)企業(yè)一方面積極拓展進(jìn)口渠道；另一方面也在通過綠地勘查、購買礦權(quán)等多種方式積極參與海外鋯鈦砂礦資源的開發(fā)。

研究區(qū)位于西非某國西部沿海地區(qū)的綠地，總面積約2100平方公里，地質(zhì)工作程度較低，河網(wǎng)密布，沼澤叢生，交通極為不便，常規(guī)找礦方法效率低下。筆者所在團(tuán)隊(duì)利用Landsat-8的OLI/TIRS多光譜數(shù)據(jù)對研究區(qū)進(jìn)行了遙感地質(zhì)特征的提取，并通過遙感地質(zhì)解譯、遙感異常篩選和綜合分析，結(jié)合區(qū)域地質(zhì)背景及成礦條件，進(jìn)而確定找礦遠(yuǎn)景區(qū)，取得了較好的階段性成果。

1 區(qū)域地質(zhì)背景及成礦條件

研究區(qū)位于西非克拉通博韋（Bove）向斜盆地的西南翼、大西洋近岸沉降區(qū)中。區(qū)內(nèi)出露地層簡單，主要為新生界第四系及古生代地層，其中第四系地層又可分為第四系濱海和海相沉積物層以及陸相沉積物層。濱海和海相沉積物層在區(qū)內(nèi)出露范圍最廣，巖性為砂層、砂土及淤泥；陸相沉積物層巖性為砂礫、粘土及砂層；古生代地層巖性為中粗粒石英砂巖。區(qū)域上巖漿巖較為發(fā)育，從上太古代到中生代均有出露，巖石類型包括角閃巖、花崗巖以及橄欖巖等，并以花崗巖類為主。新構(gòu)造運(yùn)動(dòng)發(fā)育不明顯，對重砂礦的形成和分布沒有明顯的控制作用。

根據(jù)礦石元素分析儀對黑云母花崗巖及石英砂巖露頭速測的結(jié)果，顯示露頭中的鋯鈦含量有一定富集，見表1。因研究區(qū)屬于熱帶季風(fēng)氣候，終年高溫多雨，基巖的物理-化學(xué)風(fēng)化作用強(qiáng)烈，有利于重礦物從巖石中分離，經(jīng)風(fēng)化剝蝕后，可為重砂礦的形成提供豐富的物質(zhì)來源。

表1 工區(qū)內(nèi)露頭的礦石元素分析速測結(jié)果一覽表

區(qū)內(nèi)的地貌類型包括海灘、階地、水下沙壩、沙洲及濕地紅樹林，其中砂質(zhì)海灘對成礦較為有利。

綜上所述，區(qū)內(nèi)物源豐富、氣候條件較好、有利于形成重砂礦的地形地貌，故存在形成重砂礦的可能性。

2 遙感影像的選取和預(yù)處理

選擇Landsat-8的OLI/TIRS影像數(shù)據(jù)做為數(shù)據(jù)源，軌道號分別為202～053、202～54、203～52、203～53，多光譜波段的空間分辨率為30m，全光譜波段的空間分辨率為15m。

使用ENVI軟件自帶的FLAASH大氣校正模塊，有效去除水蒸氣/氣溶膠散射效應(yīng)，同時(shí)基于像素級的校正，矯正目標(biāo)像元和鄰近像元交叉輻射的“鄰近效應(yīng)”，并對由于人為抑止而導(dǎo)致波譜噪聲進(jìn)行光譜平滑處理。

由于該地區(qū)開發(fā)程度低，多為無人區(qū)，很難尋找到居民點(diǎn)和標(biāo)志性建筑物作為校正采樣點(diǎn)，并且道路也多沿著河流分布，基本失去幾何校正的采點(diǎn)作用。若按照通常采點(diǎn)，控制點(diǎn)不足，很難保證均勻采點(diǎn)，故通過對比影像圖、已知地形圖和地質(zhì)圖及規(guī)模較大的河流和巖石、植被的界線和交叉部位，以“地貌點(diǎn)”和“地形點(diǎn)”做為幾何校正控制點(diǎn)。在多項(xiàng)式和克呂金插值數(shù)學(xué)模型下幾何校正總誤差控制在1個(gè)像元之內(nèi)。

基于OLI的特征，采用不同的波段和波譜計(jì)算、比較和相關(guān)性分析，對區(qū)內(nèi)的鐵礦化信息、羥基、碳酸根離子信息和矽卡巖等蝕變信息進(jìn)行定性和半定量分析。將OLI564和OLI754波段組合作為地層和蝕變巖解譯的基礎(chǔ)影像圖，地層、巖土、水系、植被等都有較好的影像反映。

3 遙感影像特征

區(qū)內(nèi)遙感影像特征較為簡單（見圖1），按成因形態(tài)分類原則，可分10種主要地貌單元。在衛(wèi)片上，各地貌單元有以下影像特征：

（1）低臺地影像特征：灰到深灰色，臺面較平整，由陸向海緩傾斜，近海地段常被風(fēng)成砂或海成砂堤覆蓋，多被河脊海灣，瀉湖分割而呈不連續(xù)狀，常以低矮山包似饅頭狀鑲嵌于濱海平原上。低臺地巖性為寒武系變質(zhì)巖和花崗巖，海拔高度35m～80m。

（2）洪積階地、洪積扇、洪積裙影像特征：色調(diào)多變，與鄰近地貌單元的色調(diào)常有區(qū)別。形態(tài)較規(guī)整。或扇形為洪積扇，或幾個(gè)扇形相連呈裙?fàn)顬楹榉e裙。

（3）沖積階地影像特征：影紋細(xì)膩，沿河流兩側(cè)樹枝狀展布。

（4）三角洲平原影像特征：平原面開闊平坦，略向海作2°～3°微傾斜，溝渠縱橫，牛扼湖等河流改道遺跡顯示清楚。

（5）海積平原影像特征：較大規(guī)模的海積平原面平坦，微向海傾，北東向不連續(xù)的幾列砂堤位于平原上，砂堤內(nèi)灰白色調(diào)，樹木密集時(shí)呈黑色。海積平原高程小于7m。

（6）瀉湖平原影像特征：影像細(xì)膩柔和，積水色調(diào)更深。

（7）海成砂堤影像特征：海成砂堤是區(qū)內(nèi)成礦地貌單元之一，形成時(shí)代早為Q3者稱老砂堤，以后形成者稱新砂堤。老砂堤色調(diào)較深，灰色為主，部分灰白色(受風(fēng)成砂影響)，粗糙。老砂堤離海遠(yuǎn)，位于海岸線內(nèi)側(cè)，延長方向與海岸線一致。規(guī)模大小不一，長約500m～800m，寬約100m～1000m不等，標(biāo)高一般在3m～20m之間，常背靠臺地，或位于瀉湖平原、河?；旌隙逊e平原之上，砂島上也見產(chǎn)出。形態(tài)多樣，有規(guī)則或不規(guī)則長條狀，有短條狀或圓狀。新砂堤色調(diào)淺于老砂堤，呈灰白色，形成時(shí)代越新，調(diào)越淺，目前尚在繼續(xù)擴(kuò)展的砂堤色調(diào)最淺。新砂堤在海邊廣泛發(fā)育，規(guī)模大小不一，空間位置及形態(tài)特征在航片上表現(xiàn)最為清楚，毗岸砂堤背靠平原或臺，呈彎月狀或紡睡狀、細(xì)條狀、開叉狀，堤間常夾暗色調(diào)洼地；攔灣砂堤籍丘體或臺地作依托，向一端發(fā)展。攔封海灣，逼灣口位移；砂咀頂部鉤形特征明顯。

（8）砂礫質(zhì)海灘影像特征：在潮間帶地貌中色調(diào)較淺，結(jié)構(gòu)較粗糙。有基巖出露海灘見黑色凸起的斑點(diǎn)或斑塊。

（9）泥質(zhì)海灘影像特征：灰結(jié)構(gòu)細(xì)膩，灰朦朦似霧狀，有密如蛛網(wǎng)的潮水細(xì)溝。

（10）紅樹林海灘影像特征：密布淡黑色小斑點(diǎn)，紅樹林越茂盛，色調(diào)越深。灘面潮水溝也較發(fā)育。

4 鐵染信息提取

以三價(jià)鐵和二價(jià)鐵離子的光譜反射波譜特征為出發(fā)點(diǎn)，尋找出特征譜帶所對應(yīng)的OLI波段，采用適當(dāng)?shù)倪b感增強(qiáng)方法提取到鐵礦化信息。含三價(jià)鐵（Fe3+）礦物的反射光譜曲線圖見圖2。

圖2 含F(xiàn)e3+礦物的反射光譜曲線圖

OLI的2、4、5和6波段是反射谷和反射峰交替出現(xiàn)，光譜亮度值反差較大，將這四個(gè)波段進(jìn)行組合，有利于鐵化信息的提取。為方便后期與其他地質(zhì)信息進(jìn)行綜合分析，使獲得的類鐵礦化信息半量化處理，即根據(jù)影像圖對應(yīng)的灰度直方圖分布特征和信息分布規(guī)律，設(shè)定閾值，將各類鐵礦化信息分級量化處理，再將量化后的礦化信息統(tǒng)一轉(zhuǎn)化到ARCGIS軟件平臺上，不斷地和已知的、已獲取的構(gòu)造、巖性等信息進(jìn)行疊加綜合分析，來正確的判斷蝕變信息的分布特征、剔除一些假的或無用的礦化信息。獲得的OLI的類鐵化蝕變分布及強(qiáng)度見圖3。

圖3 鐵染異常分布圖（藍(lán)色區(qū)域?yàn)殍F染異常區(qū)）

5 遙感異常查證

圖4 砂灘上的鐵染異常區(qū)照片

根據(jù)遙感影像特征及鐵染異常分布形態(tài)的不同，從圖3可見，藍(lán)色的鐵染異?？煞譃楹０毒€附近條帶狀異常區(qū)、河流入?？谔幜骶€狀或條帶狀異常區(qū)以及內(nèi)陸團(tuán)塊狀、星點(diǎn)狀異常區(qū)。通過現(xiàn)場調(diào)查發(fā)現(xiàn)前兩種鐵染異常區(qū)地表通常為灰白色、土黃色夾薄層狀、星點(diǎn)狀黑色砂層或混合砂層（見圖4），而內(nèi)陸的異常區(qū)地表通常為鐵紅色、紅褐色或土黃色砂土所覆蓋（見圖5）。利用礦石元素分析儀對部分內(nèi)陸鐵染異常區(qū)內(nèi)的樣品進(jìn)行測試，測試結(jié)果表明樣品中鐵元素含量一般2%～5%，最高可達(dá)24.19%，故推測為鐵質(zhì)氧化所致，非鋯鈦砂礦致異常。

圖5 農(nóng)田中的鐵染異常區(qū)照片

6 遙感找礦預(yù)測與找礦靶區(qū)圈定

鑒于重礦物鋯英石與鈦鐵礦的含量存在正相關(guān)性，且砂質(zhì)海灘對鋯鈦砂礦成礦較為有利，在高分遙感解譯及野外調(diào)查驗(yàn)證的基礎(chǔ)上，初步判斷海岸線附近遙感影像為亮白色、灰白色調(diào)，且重疊或附近有流線狀或條帶狀鐵染異常分布的區(qū)域成礦可能性較大?；谝陨显瓌t，圈定2處鋯鈦砂礦找礦靶區(qū)（圖3中紅色圖框范圍），并利用人工沖擊鉆在靶區(qū)1內(nèi)施工了6個(gè)鉆孔進(jìn)行驗(yàn)證，經(jīng)樣品采集及礦石元素分析儀速測，發(fā)現(xiàn)一中等規(guī)模的鋯鈦砂礦體，取得了較好的找礦效果。

7 結(jié)論

（1）采用中高分辨率的Landsat-8的OLI/TIRS多光譜數(shù)據(jù)，利用ENVI遙感圖像處理平臺，在地質(zhì)工作程度低、交通困難的區(qū)域開展遙感地質(zhì)工作，具有快速、經(jīng)濟(jì)和高效的特點(diǎn)。

（2）根據(jù)濱海型砂礦的控礦要素及礦化特征，以O(shè)LI/TIRS多光譜數(shù)據(jù)為基礎(chǔ)，確定地形地貌及鐵染異常等解譯重點(diǎn)，對找礦靶區(qū)的圈定具有明確的指導(dǎo)意義。

（3）在總結(jié)區(qū)域地質(zhì)背景和成礦條件的基礎(chǔ)上，結(jié)合遙感解譯成果，開展成礦預(yù)測，圈定2個(gè)找礦靶區(qū)，并施工少量鉆孔對其中1個(gè)靶區(qū)進(jìn)行驗(yàn)證，發(fā)現(xiàn)一中等規(guī)模的鋯鈦砂礦體。