王明明，羅清威，張克川，李新超

(1.河南省有色金屬礦產(chǎn)探測工程技術研究中心，河南鄭州 450016；2.河南省有色金屬地質礦產(chǎn)局第七地質大隊，河南鄭州 450018)

0 引言

近些年來，隨著遙感技術的迅速發(fā)展，其具有快捷、經(jīng)濟、易于數(shù)據(jù)綜合分析等優(yōu)點，已被廣泛應用于傳統(tǒng)地質工作中(成功等，2016；趙玉梅等，2016；梁軍等，2017)。結合計算機圖像處理技術，可從獲得的遙感影像中識別并提取不同巖石、礦物蝕變異常信息，為地質找礦的前期工作圈定有利靶區(qū)(常睿春等，2012；陳富強等，2017；涂其軍等，2018；陳琪等，2021)。因此，利用遙感影像提取遙感蝕變異常信息來指導找礦已成為遙感地質應用的重要領域之一。

坦桑尼亞環(huán)維多利亞湖太古宙綠巖型金礦成礦帶是著名的金礦成礦帶之一，而位于其東部的穆索馬-馬拉綠巖帶是環(huán)維多利亞湖“八大”太古宙綠巖帶之一，該綠巖帶構造發(fā)育，巖漿活動強烈，具備良好的成礦條件。目前該綠巖帶已發(fā)現(xiàn)有馬拉、布亨巴、伊昆古等大中型金礦床，這些金礦床的成礦類型和成礦地質條件大致相似，均為剪切帶或構造破碎帶控制的構造蝕變巖型金礦，其礦化蝕變特征主要為黃鐵礦化、褐鐵礦化及硅化等，次為綠泥石化、絹云母化及碳酸鹽化等(崔小軍等，2014；王建光等，2017；楊友等，2017；李俊生等，2018)。

本文以PL7233金礦權區(qū)為例，通過對獲取的遙感影像采用掩膜技術剔除水體、植被等干擾信息提取有用的礦化蝕變異常信息(陳三明等，2010；錢建平等2014；穆媛芮等，2017；梁棪先等，2020)，圈定有利找礦靶區(qū)，并結合地質資料、物探等手段進行綜合分析研究，再對圈定的有利靶區(qū)利用山地工程(探槽、鉆探)進行異常成果查證，取得了較好的找礦成果。

1 研究區(qū)地質概況

1.1 區(qū)域地質特征

穆索馬-馬拉綠巖型金礦成礦帶位于坦桑尼亞環(huán)維多利亞湖金礦成礦區(qū)東部，該綠巖帶主要由太古界尼安薩群、元古界布科班群、新生界新近系和第四系組成，其地層為花崗巖-綠巖體，下部為太古界尼安薩群，上部為元古界布科班群的碎屑沉積巖所覆蓋，呈不整合接觸。巖性主要以長英質、鐵鎂質火山巖及BIF等為主，夾有少量凝灰?guī)r和變火山巖。其中太古界尼安薩群為金礦(化)體主要賦存層位。區(qū)域構造主要為近EW、NE和NW向的3組近平行剪切構造帶組成，其延伸長度超過200 km，貫穿整個穆索馬-馬拉綠巖帶。此外，區(qū)內次生褶皺、斷裂帶和剪切帶等不同類型、不同規(guī)模的構造較為發(fā)育，局部地層遭受強烈變質變形作用，這些強烈的變質變形構造活動為區(qū)內金礦床的形成的有利的成礦地質環(huán)境(羅清威等，2019)。

研究區(qū)位于穆索馬-馬拉綠巖帶東部，區(qū)內出露的地層有太古界尼安薩群、新生界新近系和第四系(見圖1)。其中太古界尼安薩群主要分布在北部和中西部，巖性為斜長片麻巖、斜長角閃片麻巖、黑云斜長片麻巖和花崗片麻巖等；新近系主要分布在東部，巖性為粗面巖；第四系分布在西南部、南部和北東部，多為片麻巖、粗面巖等殘積物及沖-洪積物等。地層總體走向為NE-NEE向，南傾，傾角為50°～70°(羅清威等，2019)。

圖1 研究區(qū)區(qū)域位置及地質圖

區(qū)內斷裂構造發(fā)育，主要為NE、NW和近EW向，其中NW、NWW向剪切構造為主要的控礦構造，剪切構造帶發(fā)育在片麻巖系中，巖石多見硅化、褐鐵礦化、鐵染等蝕變，在構造帶兩側的巖石中可見碎裂結構、角礫狀構造及網(wǎng)脈狀構造等。

區(qū)內分布有少量正長斑巖脈、石英脈及輝綠巖脈等。

1.2 礦床地質特征

研究區(qū)內已發(fā)現(xiàn)一條金礦(化)體K1，該金礦(化)體由剪切構造帶控制，鉆孔中可見蝕變輝綠巖脈穿插其中。礦(化)體控制長約550 m，賦存標高1127～1326 m，控制斜深約130 m。

礦體產(chǎn)狀與Fb1構造蝕變帶基本一致，金品位1.00×10-6～38.85×10-6，平均品位8.69×10-6，品位變化系數(shù)為2.67；厚度為3.68～5.76 m，平均厚度4.75 m。此外，區(qū)內還揭露有Au2、Au3金礦(化)體，金品位0.53×10-6～1.04×10-6，平均厚度1.78 m。金礦(化)體形態(tài)較簡單，整體呈脈狀、透鏡狀產(chǎn)于構造帶中，礦(化)體圍巖蝕變主要有硅化、綠泥石化、絹云母化等，金屬礦化主要有黃鐵礦化、磁黃鐵礦化、褐鐵礦化等，其中硅化、綠泥石化、黃鐵礦化、褐鐵礦化與金礦化的關系較為密切(李水平等，2016)。

1.3 礦床成因類型分析

根據(jù)環(huán)維多利亞湖東部綠巖帶型金礦床成礦規(guī)律和成礦模式，結合研究區(qū)的礦床地質特征等，初步分析認為研究區(qū)內金礦床為受剪切構造帶控制的構造蝕變巖型金礦床。

現(xiàn)有資料表明，太古界尼安薩群地層中的超鎂鐵質、鎂鐵質火山巖的原始含金量偏高，為礦床的形成提供了天然礦源層。此外，當來自后期的中酸性侵入體攜帶的成礦流體在上升過程中穿過這套綠巖帶地層時，成礦流體與圍巖相互作用形成了Eh值較低且富硫的礦源體。區(qū)域上的剪切構造帶及次級斷裂、褶皺為成礦流體提供了通道和富集空間，在特定的溫度和壓力環(huán)境下，成礦物質發(fā)生沉淀后形成礦(化)體，而后期熱液活動可能使成礦物質進一步富集，使得區(qū)內礦(化)體在時間尺度上具有多期次、多階段的成礦特點(崔小軍等，2015；孫宏偉等，2015；張耀等，2017；李俊生等，2018)。

2 Crosta方法基于TM數(shù)據(jù)蝕變信息提取

目前，在遙感地質中常用的蝕變異常信息提取方法主要有波段比值法、LS-Fit法、Crosta法及匹配濾波法等。Crosta法具有無需實測光譜，效果穩(wěn)定等優(yōu)點，且經(jīng)過Crosta技術處理后的TM/ETM和ASTER等遙感影像具備較好的找礦指示能力(彭光雄等，2013)。因此，Crosta法是目前遙感礦化蝕變異常信息提取應用最主要的方法(薛曉剛，2009；鄧會娟等，2014；汪子義等，2016；趙玉梅等，2016)。

2.1 遙感數(shù)據(jù)源選取

根據(jù)研究區(qū)植被、地表水體、降雨及地形等特征(汪樹棟等，2014)，本次選取坦桑尼亞旱季的Landsat-5遙感影像(TM數(shù)據(jù))，日期為2016年10月28號，軌道號為170/061。

2.2 遙感影像預處理

研究區(qū)位于坦桑尼亞東北部，區(qū)內植被較為發(fā)育，地表水體多為季節(jié)性沖溝。為了盡量減少蝕變異常信息提取時干擾因素的影響，本次對選取的遙感影像通過ENVI5.3平臺進行一系列數(shù)據(jù)預處理操作過程，具體見圖2。

圖2 遙感影像數(shù)據(jù)預處理流程圖

2.3 遙感礦化蝕變信息提取

利用Crosta法對預處理后的TM數(shù)據(jù)進行處理。該方法是通過選取第1、3、4、5波段組合利用主成分分析法(PCA)提取鐵染蝕變異常信息，其特征向量矩陣見表1；選取第1、4、5、7波段組合利用主成分分析法提取羥基蝕變異常信息(于彩虹等，2015；梁軍等，2017；趙慧童和王萍，2017)，其特征向量矩陣見表2。

從表1可以看出，在PC4分量中TM3系數(shù)與TM1及TM4系數(shù)符號相反，且與TM5系數(shù)符號相同。根據(jù)代表鐵染物主分量的判斷準則(張玉君等，2003)，選擇PC4分量作為鐵染蝕變異常信息提取的主分量。

表1 TM 1、3、4、5波段PCA特征向量矩陣

表2 TM 1、4、5、7波段PCA特征向量矩陣

從表2可以看出，在PC4分量中TM5系數(shù)與TM4及TM7系數(shù)符號相反，且與TM1系數(shù)符號相同。根據(jù)代表羥基化物主分量的判斷準則(張玉君等，2003)，選擇PC4分量作為羥基蝕變異常信息提取的主分量。

主分量門限化技術是以PCA法得到的主分量數(shù)據(jù)為基礎，利用標準離差對蝕變異常強弱進行量化分級。該技術是把圖像灰度均值(X)作為背景值，以標準離差(σ)作為量化尺度，用標準離差倍數(shù)作為閾值，即利用公式(X+kσ)對各級異常進行等級劃分，其中k值一般取1.0～3.0(錢建平等，2013；吳繼煒，2015)。通過ENVI5.3平臺的Compute Statistics功能得到各異常蝕變主分量的基本統(tǒng)計數(shù)據(jù)，見表3。

表3 各異常蝕變主分量的基本數(shù)據(jù)統(tǒng)計表

其中，無異常為最小值x+2σ，一級異常為(x+2σ)～(x+2.5σ)，二級異常為(x+2.5σ)～(x+3σ)，三級異常為x+3σ最大值，其蝕變異常等級劃分結果見圖3、圖4。

圖3 研究區(qū)鐵染蝕變異常信息提取結果

圖4 研究區(qū)羥基蝕變異常信息提取結果

通過與鄰近的馬拉金礦對比可以看出，研究區(qū)共圈定4處鐵染蝕變異常區(qū)，3處羥基蝕變異常區(qū)，其中位于研究區(qū)西北部的Ⅰ號鐵染蝕變異常與Ⅰ號羥基蝕變異?；局睾?，該區(qū)域為重點找礦靶區(qū)。

3 物探方法及裝置

為了進一步縮小找礦靶區(qū)，查證遙感影像蝕變異常信息，對異常區(qū)域的找礦潛力作出評價，對研究區(qū)內遙感影像蝕變異常圈定的找礦靶區(qū)開展物探測量工作。

3.1 方法選擇及裝置設置

現(xiàn)有資料表明，在合適的地電條件下，時間域激發(fā)極化法可以作為異常查證方法進行金礦(體)勘查(彭俊等，2017；李水平等，2019)。本次物探工作選擇激電中梯剖面與激電測深工作對異常區(qū)進行剖析，了解其水平展布方向及垂向延深情況。

激電中梯剖面工作采用時間域激發(fā)極化法，工作裝置為中間梯度裝置，供電設備為重慶奔騰WDZF-5型激電發(fā)射系統(tǒng)，接收機為WDJS-2數(shù)字直流激電接收機，選擇供電周期16 s，延時150 ms，基本積分間隔40 ms，采用雙向短脈沖供電模式，占空比1:1，觀測參數(shù)I、△U、ηs。

激電測深工作采用時間域激發(fā)極化法，工作裝置為對稱四極裝置，供電設備為重慶奔騰WDFZ-5型激電發(fā)射系統(tǒng)，接收機為WDJS-2數(shù)字直流激電接收機，選擇供電周期16 s，延時150 ms，基本積分間隔40 ms，采用雙向短脈沖供電模式，占空比1:1，觀測參數(shù)I、△U、ηs。其中，激電測深工作電極距離排列表見表4。

表4 激電測深工作電極距離排列表(m)

3.2 物探數(shù)據(jù)處理及資料解譯驗證

3.2.1 研究區(qū)巖(礦)石標本物性特征

根據(jù)不同巖石及礦化蝕變會導致巖石電性參數(shù)等物理性質存在差異的原理(張同中等，2010)，首先對區(qū)內常見的巖(礦)石電性參數(shù)進行了測量統(tǒng)計，詳見表5。

表5 研究區(qū)巖(礦)石電性參數(shù)特征統(tǒng)計表

從表5的統(tǒng)計數(shù)據(jù)中可以看出，構造蝕變極化率最高，電阻率最低；粗面巖電阻率最高，極化率次之；石英脈、輝綠巖等電阻率較高，極化率最低。表明構造蝕變巖具有“低阻、高極化”特征，即研究區(qū)內的構造蝕變帶與圍巖之間存在明顯的電性參數(shù)差異，為在研究區(qū)開展激發(fā)極化法提供了工作前提，利用該方法可以探尋有利成礦位置，圈定目標靶區(qū)。

3.2.2 成果解譯驗證

結合地質資料及取得的成果，在認為成礦有利部位布設物探剖面，開展了物探工作，通過后期鉆孔驗證情況來看，選擇的物探方法對深部礦體有較好的反映，具有很好的指示意義。

激電中梯剖面顯示，視極化率曲線南端較平緩，異常較弱，對應為粗面巖區(qū)，在構造蝕變帶南側視極化率有相對較強異常顯示，極值為1.15%。而電阻率在構造蝕變帶兩端為高阻，中間為低阻異常。根據(jù)激電中梯剖面異常特征，布設了12個激電測深點，其視電阻率特征表現(xiàn)為淺層呈現(xiàn)低阻異常，在04～08號點深部存在高阻異常，在16～28號點處存在低阻異常帶，且一直向下延深，異常形態(tài)向南傾斜，推測與構造蝕變帶有關；視極化率在20～36號點間存在一相對高極化異常體，隨極距增大異常相對封閉，極化體形態(tài)向南傾斜，以20號點為界，南側視極化率值明顯高于北側，形成相對較高極化率異常，這與視電阻率低阻異常推斷的構造蝕變帶基本吻合，其后期鉆孔驗證成果見圖5。

圖5 研究區(qū)典型地物綜合剖面圖

4 找礦標志及找礦前景

4.1 研究區(qū)找礦標志

(1)地質找礦標志

研究區(qū)金礦(化)體賦存在太古界尼安薩群的片麻巖、片巖中，且含礦性與圍巖成分有一定的關系，出現(xiàn)蝕變輝綠巖對成礦最有利。金礦(化)體均嚴格受剪切構造帶控制，剪切構造帶即為礦體賦存空間，是較為直接的找礦標志，特別是在次級斷裂構造密集發(fā)育、交匯或復合部位。構造蝕變帶具有明顯的礦化蝕變特征，綜合分析認為硅化、黃鐵礦化、褐鐵礦化、綠泥石化與金礦化關系最為密切。此外，地表褐鐵礦化及硅化蝕變的存在也是較好的找礦標志。

(2)遙感解譯蝕變標志

對獲取的遙感數(shù)據(jù)進行預處理后，利用Crosta法可有效提取蝕變異常信息，研究區(qū)內共圈定4處鐵染異常蝕變、3處羥基異常蝕變，均呈線狀、面狀展布，其中Ⅰ號鐵染蝕變異常與Ⅰ號羥基蝕變異常基本重合，呈線狀分布，表現(xiàn)出“弱鐵染蝕變異常+強羥基蝕變異?！蹦Ｊ?。通過已有地質資料和已知礦(化)點對比研究，證實遙感蝕變異常信息可以作為找礦預測的標志。

(3)地球物理異常標志

研究區(qū)內金礦(化)體與圍巖相比，表現(xiàn)為“低阻、相對高極化”的異常特征。因此，針對性地利用時間域激發(fā)極化法可以幫助查明金礦(化)體的展布及延深情況，以指導后期鉆孔施工。

(4)其他標志

研究內或附近的民采坑、豎井等采礦活動。

4.2 研究區(qū)找礦前景分析

實踐證明，大多數(shù)情況下受剪切構造帶控制的金礦床特征表現(xiàn)為大斷裂成大礦，小斷裂只能形成中、小型礦床(點)，特別是特大型金礦床往往位于區(qū)域構造系的復合交叉及轉折部位(劉軍和朱谷昌，2012)。穆索馬-馬拉綠巖帶中產(chǎn)出的眾多金礦床多成群、成帶出現(xiàn)，且其控礦和容礦構造均為剪切構造帶，較好地反映了地殼深部的構造活動(司建濤等，2017；王建光等，2017；楊友等，2017)。研究區(qū)金礦床類型為受剪切構造帶控制的構造蝕變巖型金礦，目前已發(fā)現(xiàn)并初步控制的K1金礦(化)體賦存于NWW向的斷裂構造蝕變帶中。且現(xiàn)有資料表明，K1金礦(化)體南側存在一條全隱伏的推測構造蝕變帶，其具有相似的成礦條件，結合已有資料和構造蝕變巖型金礦床特征，認為該推測構造蝕變帶具有較大的成礦潛力。因此，研究區(qū)找礦勘查應在點上有所突破的情況下，根據(jù)有效勘查手段圈定的異常區(qū)，進一步查明對成礦有利的構造體系，以擴大找礦成果。

5 結論

(1)利用獲得的遙感影像，對研究區(qū)進行了礦化蝕變異常信息提取，圈定了4處鐵染蝕變異常、3處羥基蝕變異常，其中“弱鐵染蝕變異常+強羥基蝕變異?！蹦Ｊ骄哂泻芎玫恼业V指示意義。

(2)通過后期鉆孔施工驗證，分析認為研究區(qū)金礦(化)體具有“低阻、相對高極化”異常特征，其中低電阻率推測與構造破碎帶有關，高極化率推測是由黃鐵礦等金屬硫化物所引起，且礦(化)體的含金性與視電阻率、視極化率的高低具有一定的相關性。因此，“低阻、相對高極化”激電異常體為成礦有利部位。

(3)研究區(qū)金礦床類型是受剪切帶控制的構造蝕變巖型金礦。根據(jù)構造蝕變巖型金礦類型特征，結合區(qū)內已有資料表明，研究區(qū)具有良好的找礦前景，特別是剪切帶的轉折部位及次級斷裂分支交叉復合部位，是將來擴大找礦成果的一個重要方向。