董麗娜,張 微,2,王 雪,陳 玲,楊金中,莫子奮
(1.中國國土資源航空物探遙感中心,北京 100083;2.北斗航天衛星應用科技集團,北京 100070;3.中國地質大學(北京)地球科學與資源學院,北京 100083;4.江西省核工業地質局265大隊,鷹潭 335000)
20世紀70年代以來,隨著航天技術和計算機技術的飛速發展,遙感技術在地質、礦產、國土、環境等調查領域的應用越來越受到重視[1-2];其宏觀性、現勢性和多時效性等特點,更在近年來的區域地質調查、礦產資源勘查、工程地質調查和地質災害動態監測等工作中得到了廣泛應用[3]。經過十幾年的應用與發展,隨著空間分辨率和光譜分辨率的不斷提高[3],利用多光譜衛星遙感數據開展礦產資源勘查與資源評價已經成為一種快速經濟的找礦手段[4-5]。本文以 ALOS,SPOT5 和 ETM+等多源遙感數據為信息源,在完成圖像幾何糾正與配準和圖像增強的基礎上,對江西盛源火山盆地進行1∶5萬比例尺遙感地質解譯和遙感蝕變異常信息提取與篩選;并結合已有地質資料分析尋找鈾礦的前景,為該區的礦產資源勘查和資源評價提供基礎地質資料和參考依據。
研究區盛源火山盆地位于江西省東北部,隸屬于鷹潭市。浙贛、皖贛及鷹廈等鐵路和320國道等多條高等級公路穿越研究區,交通十分便利(圖1)。

圖1 研究區交通位置圖Fig.1 Traffic position of study area
從區域上看,鄱陽湖及其周邊地區的地貌類型復雜,山地、丘陵、崗地、平原和湖泊洼地兼而有之。研究區內則以崗地、平原、湖泊和洼地為主,主要由上更新統亞粘土直接覆蓋于前第四系之上,并組成二級階地與一級階地;信江自東向西橫貫研究區中北部,沿江地區由全新世晚期堆積物組成河漫灘。
江西盛源火山盆地位于揚子準地臺與華南褶皺系接壤處、信江凹陷中段南緣,處在近EW向萍鄉—廣豐深斷裂、NE向鷹潭—安遠深斷裂與NW向鷹潭—瑞昌大斷裂的復合部位。基底主要由元古宙—早寒武世淺變質巖系和加里東期花崗巖組成,蓋層為上侏羅統火山巖地層和白堊系砂巖紅層。在新近紀特別是燕山期,太平洋板塊開始向歐亞板塊俯沖,被動大陸邊緣充分發育的硅鋁殼發生從熔改造,過渡巖漿房的巖漿分異結晶,形成強分異的富鈾巖漿;隨著俯沖速度的減慢和俯沖帶逐漸向大洋方向遷移,大陸邊緣因受到強烈擠壓而逐漸向拉張轉化,沿結合帶發育的斷裂帶開始復活;過渡巖漿房的富鈾巖漿及其聚集的揮發份和流體迅速上移和噴出,形成大面積的裂隙式火山巖;在多組構造的結合帶附近發生廣泛的鈾化[1],聯合控制了帶內的陸相火山活動和鈾成礦作用。其中萍鄉—廣豐深斷裂不但控制了2個Ⅰ級大地構造單元的地殼演化,而且控制了帶內的構造-巖漿活動及礦化活動,使得盛源盆地及其鄰區成為贛—杭構造火山帶內的1個重要的產鈾火山盆地,具有良好的找鈾礦前景。研究區及其周邊地區礦產資源豐富,已探明儲量的有金、銀、鉛、鋅、鈾、瓷土、石膏、硅等30余種,其中銀、石膏、瓷土和鈾礦均十分豐富[6]。
本文采用 ALOS,SPOT5和ETM+衛星數據,首先以控制點-多項式擬合校正方法對原始數據進行了幾何糾正和配準;然后利用基于地理坐標的無縫鑲嵌方法對圖像進行鑲嵌,并運用亮度匹配等方法確保鑲嵌圖像的質量[7]。
圖像增強主要是突出不同地物之間波譜特征的差別,以及空間形態、邊緣、線條及紋理結構特征(如蝕變帶、接觸帶、破碎構造帶、線性體及地貌特征等)[8]。本文選用圖像融合[9-11]、主成分變換[5]、比值分析[6]和濾波增強等不同的圖像處理方法,通過對圖像處理結果的對比分析及相互驗證,達到提高區內巖性、構造信息提取效果的目的(圖2)。
而對局部重點地區,則針對不同的地物特征信息,運用圖像糾正、彩色合成與空間變換、圖像增強處理、圖像鑲嵌等技術手段,以達到最佳效果[12]。
盛源火山盆地位于揚子準地臺與華南褶皺系接壤處、信江凹陷中段南緣,處在近EW向萍鄉—廣豐深斷裂、NE向鷹潭—安遠深斷裂與NW向鷹潭—瑞昌大斷裂復合部位。其中萍鄉—廣豐深斷裂是1條自中元古宙以來長期發展的活動帶,它不但控制了2個Ⅰ級大地構造單元的地殼演化,而且控制了帶內的構造-巖漿以及礦化活動;特別是在燕山期,該活動帶與太平洋板塊向歐亞板塊俯沖,聯合控制了帶內的陸相火山活動和鈾成礦作用。故遙感地質解譯中先利用ETM+圖像建立研究區遙感構造格架(圖3),進行區域地質背景綜合分析。
再利用ALOS,SPOT5等空間分辨率優于2.5 m的數據,對重點區進行詳細的巖性和構造解譯(圖4)。
依據各地質體的影像特征,結合已有地質資料,在研究區內共解譯出3組不同方向的線性構造;并分別采用紋理分析和卷積對NE,NW和近EW向構造進行增強處理,解譯結果對推測隱伏斷裂構造具有很大意義。為了更好地對斷裂構造的規模和空間展布方向進行解譯,還利用ALOS圖像與空間分辨率為30m的DEM進行疊加顯示,對生成的三維遙感影像圖進行詳細解譯。最后,對解譯出的典型構造和巖性進行了野外實地驗證[13],并根據驗證結果完善了解譯標志、提高了解譯精度。
此外,通過濾波增強處理,還發現若干疑似火山機構,在三維遙感影像圖上有較好的反映(圖5)。

圖5 疑似火山機構的環形構造影像圖Fig.5 Circular structure images of suspected volcanic institutions
研究區地表覆被嚴重,尤其是植被、陰影區、云、水體和沼澤地等對遙感蝕變異常信息提取的干擾很大。結合以往經驗,通過改進“去干擾-主成分分析-異常篩選”技術[14-15],成功地去除了上述干擾信息,實現了不同水平遙感異常的精確定位;在此基礎上,利用ETM+數據進行了遙感蝕變異常信息的提取、分類和分級[16](圖 6)。

圖6 研究區遙感異常圖Fig.6 Remote sensing anomaly map of study area
在以上處理的基礎上,以遙感影像為基礎,綜合圖像處理、地質解譯和蝕變異常信息提取結果,在楊梅塘、壩口和金源地區圈定出3處環形構造,初步推測為火山機構(在下文的綜合找礦信息及鈾礦前景分析中作具體論述);并對需要補充開展的驗證工作提出了建議。
以從遙感圖像獲取的解譯結果和異常信息為主要預測依據,結合已有相關資料,分別從地層、構造和遙感異常等方面[17-18]進行綜合找礦預測和鈾礦前景分析,提出曲溪、余壩和盛源3處遙感找礦有利區,分述如下:

圖7 曲溪遙感找礦預測圖Fig.7 Rem ote sensing prospecting p rediction m ap for Quxi zone
1)曲溪遙感找礦有利區(圖7)。該區段位于研究區的西部、盆地西緣外側的推覆體內。地層以下白堊統贛州群(K1g)為主,西部發育有少量的上侏羅統打鼓頂組(J3d)和鵝湖嶺組(J3e)。區段內地質構造十分發育,其中EW向葉家—盛源斷裂與NE向壩口—貴溪斷裂、NW向余干—鷹潭斷裂聯合控制著火山坳陷西邊緣。賦礦地層主要為上侏羅統打鼓頂組(J3d)熔結凝灰巖,受上述3組斷裂構造影響,形成了規模較大的順層破碎帶,現已在該破碎帶中發現大型鈾礦床。該區段遙感圖像的影紋較平滑,色調較均勻單一,地形起伏不大,形成一系列低矮山丘,局部地區環形構造較為明顯。遙感蝕變異常較為發育,主要以羥基異常為主。
經研究發現,該區段內已發現的礦床因受后期盆地西部抬升作用的影響,在重力和伸展作用的影響下,向西滑覆到現在位置。通過不同尺度遙感圖像、遙感異常對比分析和疊加顯示發現,礦床東部發育有較為明顯的環形構造,且環形構造與遙感異常套合較好,因此建議在該區東部進一步開展找礦工作。

圖8 余壩遙感找礦預測圖Fig.8 Remote sensing prospecting prediction map for Yuba zone
2)余壩遙感找礦有利區(圖8)。該區段位于盛源火山盆地的南部,為火山巖分布區的三角地帶。該區段火山盆地基底為加里東期花崗巖(γ3)和震旦—寒武系海相火山碎屑巖富鈾層,地層以上侏羅統打鼓頂組(J3d)和鵝湖嶺組(J3e)為主,北部含少量的下白堊統贛州群(K1g)。該區段地處近EW向壩口—莊家斷裂帶與NE向壩口—貴溪斷裂帶、金源—流口斷裂帶的復合部位,地層展布明顯受斷裂控制,遙感影像特征明顯。經圖像增強處理,發現若干環形構造和疑似火山隆起的影像。在這些環形構造周邊已有次花崗斑巖出露,其深部可能隱伏有次花崗斑巖體。該區段內遙感異常顯著,以羥基異常為主;局部地區羥基異常和鐵染異常套合較好,鐵染異常多出現在羥基異常邊部,分布較為稀疏,強度較弱;異常總體分布與已知破碎蝕變帶吻合較好。
綜上所述,推測深部可能存在與次火山巖有關的礦化(體)。該區段遙感找礦預測圖上的蝕變異常的空間分布特征與已知礦床形態相似,是尋找深部火山巖中低溫熱液型鈾礦床的有利地段,顯現出良好的找礦遠景。

圖9 盛源遙感找礦預測圖Fig.9 Remote sensing prospecting prediction map for Shengyuan zone
3)盛源重點遙感找礦有利區(圖9)。該區段位于盛源火山盆地東部邊緣,NNE向漆工—圳上斷裂帶和弋陽—資溪斷裂的復合部位。該區段中部以上白堊統南雄組(K2n)為主,遙感圖像中紋理結構十分發育;而兩側則以下白堊統贛州群(K1g)為主,邊部發育有少量的上侏羅統打鼓頂組(J3d)和鵝湖嶺組(J3e)。該區段內地質構造十分發育,有NE,NW和近EW向3組,其中近EW向基底斷裂對礦床控制作用明顯。受EW向和NE向次一級斷裂構造的影響,順層破碎帶十分發育。侏羅系上統火山巖出露齊全,火山碎屑巖與火山碎屑沉積巖相間出現,底部為中元古界變質巖,外圍有燕山期花崗巖體;遙感異常與斷裂構造關系較為緊密,具有較明顯的近EW向斑狀分布特征,基本上沿近EW向線性構造(帶)和礦帶斷續展布。
遙感異常以鐵染異常為主,有規律地分布在不同巖石地層單元的邊界處。該區段較多,成礦地質條件優越。建議在已知礦床周邊,針對遙感地質解譯結果與遙感異常套合較好且分布集中的有利地段,開展進一步礦產詳查工作。
1)本文利用多源遙感數據,采用圖像融合、主成分變換、比值分析和濾波增強等方法,通過對地物信息的增強處理,突出了不同影像單元的對比度,有效地圈定了巖石和地層界線,取得了較好的地質效果;同時結合數字高程模型(DEM)制作三維立體模型,提高了不同地物信息的對比度和可識別程度,更加有效地提取出線性構造、環形構造和火山機構,為該地區的找礦預測和資源評價提供了基礎數據和參考依據。
2)盡管研究區地表破壞嚴重,且大部分為白堊系紅層所覆蓋,但通過遙感異常信息提取,結合已有礦產資料,發現和圈定的遙感找礦有利區的鐵染異常均有規律地呈斑狀分布在火山碎屑巖與火山碎屑沉積巖出露帶上;且基本沿重要含礦地層單元界線呈線性展布,受構造控制特征明顯,與鈾礦成礦關系密切。
3)受客觀條件限制,本文的找礦預測方法和成果還有待進一步補充和完善,如遙感異常篩選、圖像增強等技術手段的針對性有待加強。在后續工作中將繼續進行更深入的遙感找礦研究。
[1] 張廷斌,唐菊興,李志軍,等.西藏尕爾窮銅金礦多光譜遙感地質特征與外圍找礦預測[J].國土資源遙感,2014,26(4):170-178.doi:10.6046/gtzyyg.2014.04.27.Zhang T B,Tang JX,Li Z J,etal.Remote sensing geological characteristics and ore prediction in the Ga'erqiong Cu-Au deposit,Tibet[J].Remote Sensing for Land and Resources,2014,26(4):170-178.doi:10.6046/gtzyyg.2014.04.27.
[2] 張云峰,焦超衛,李劍斌,等.陜西柞水冷水溝地區遙感地質特征及找礦預測[J].國土資源遙感,2014,26(4):187-194.doi:10.6046/gtzyyg.2014.04.29.Zhang Y F,Jiao CW,Li JB,et al.Remote sensing geological features and ore prediction in Lengshuigou,Zhashui County,Shaanxi Province[J].Remote Sensing for Land and Resources,2014,26(4):187-194.doi:10.6046/gtzyyg.2014.04.29.
[3] 關 震,吳 虹,曹 翠,等.基于ETM+6-γ 示礦信息反演的花山花崗巖鈾礦預測[J].國土資源遙感,2014,26(3):92-98.doi:10.6046/gtzyyg.2014.03.15.Guan Z,Wu H,Cao C,et al.Uranium ore prediction based on inversion of ETM+6-γmineral information in Huashan granite area[J].Remote Sensing for Land and Resources,2014,26(3):92-98.doi:10.6046/gtzyyg.2014.03.15.
[4] 鐘江文,彭 翼.小秦嶺地區遙感線性構造密集帶與金礦關系分析及找礦預測[J].國土資源遙感,2014,26(2):148-153.doi:10.6046/gtzyyg.2014.02.24.Zhong JW,Peng Y.Analysis of relationship between belts of concentrated remote sensing linear structures and gold deposits aswell as prospecting prognosis in Xiaoqinling region[J].Remote Sensing for Land and Resources,2014,26(2):148-153.doi:10.6046/gtzyyg.2014.02.24.
[5] 鄧會娟,姚聿濤,彭光雄,等.青海納日貢瑪斑巖型Cu-Mo礦床遙感蝕變異常提取與找礦預測[J].國土資源遙感,2014,26(2):154-161.doi:10.6046/gtzyyg.2014.02.25.Deng H J,Yao Y T,Peng G X,et al.Extraction of remote sensing alteration anomalies and prospecting prediction of porphyry Cu-Mo deposits in Narigongma,Qinghai province[J].Remote Sensing for Land and Resources,2014,26(2):154-161.doi:10.6046/gtzyyg.2014.02.25.
[6] 張 凱,高振兵.EH4成像系統在鈾礦地區的研究應用[J].科學之友,2013(8):32-33.Zhang K,Gao Z B.Application research of EH4 electromagnetic imaging system to a uranium deposition region[J].Friend of Science Amateurs,2013(8):32-33.
[7] 陳 功,舒 寧,王 歡,等.控制點影像庫在遙感影像幾何糾正中的應用[J].地理空間信息,2006,4(1):58-61.Chen G,Shu N,Wang H,et al.Application of image database for ground control point to geography correction of images[J].Geospatial Information,2006,4(1):58-61.
[8] 湯國安,張友順,劉詠梅,等.遙感數字圖像處理[M].北京:科學出版社,2004.Tang G A,Zhang Y S,Liu Y M,et al.Remote Sensing Digital Image Processing[M].Beijing:Science Press,2004.
[9] Ghosh A.Atmospheric correction and fusion[J].Geosciences and Remote Sensing Symposium,2010(2):112-119.
[10] Chen X,Wu J,Zhang Y.Comparison of fusion algorithms for ALOS panchromatic and multispectral images[C]//2008 International Workshop on Education Technology and Training and 2008 InternationalWorkshop on Geoscience and Remote Sensing.Shanghai:IEEE,2008:167-170.
[11] Das A,Revathy K.Image fusion techniques for remote sensed SPOT images[C]//International Conference on Computational Intelligence and Multimedia Applications.Sivakasi,Tamil Nadu:IEEE,2007:135-139.
[12] 劉玉英,郝福江.遙感地質學[M].北京:地質出版社,2010.Liu Y Y,Hao F J.Remote Sensing Geology[M].Beijing:Geological Publishing House,2010.
[13] Pirasteh S,Pradhan B,Safari H O,et al.Coupling of DEM and remote-sensing-based approaches for semi-automated detection of regional geostructural features in Zagrosmountain,Iran[J].Arabian Journal of Geosciences,2013,6(1):91-99.
[14] 況 忠,黃欣欣,況順達,等.貴州弱礦化蝕變遙感信息的分布特征[J].國土資源遙感,2014,26(2):140-147.doi:10.6046/gtzyyg.2014.02.23.Kuang Z,Huang X X,Kuang SD,et al.Distribution characteristics of remote sensing information on weakmineralization and alteration in Guizhou[J].Remote Sensing for Land and Resources,2014,26(2):140-147.doi:10.6046/gtzyyg.2014.02.23.
[15] 劉文蘭,張 微.遙感構造蝕變異常信息提取及找礦預測——以老撾為例[J].國土資源遙感,2012,24(2):68-74.doi:10.6046/gtzyyg.2012.02.13.Liu W L,ZhangW.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.doi:10.6046/gtzyyg.2012.02.13.
[16] 張 微,楊金中,方洪賓,等.東昆侖—阿爾金地區遙感地質解譯與成礦預測[J].西北地質,2010,43(4):288-294.ZhangW,Yang JZ,Fang H B,et al.Remote sensing interpretation and metallogenic prediction in themetallogenic belts of East Kunlun and Altun Mts[J].Northwestern Geology,2010,43(4):288-294.
[17] 張萬良.盛源盆地火山巖地層的時代歸屬[J].地質評論,1999,1(5):469-469.Zhang W L.Stratigraphic attribution of volcanic series in Shengyuan basin[J].Geological Review,1999,1(5):469-469.
[18] 張 微,張 偉,劉世英,等.基于核PCA方法的高分辨率遙感圖像自動解譯[J].國土資源遙感,2011,23(3):82-87.doi:10.6046/gtzyyg.2011.03.15.Zhang W,ZhangW,Liu SY,etal.Automatic interpretation of high resolution remotely sensed images by using Kernelmethod[J].Remote Sensing for Land and Resources,2011,23(3):82-87.doi:10.6046/gtzyyg.2011.03.15.