蘇丹紅海州辛卡特金礦遙感找礦方法研究應用

摘要：根據蘇丹紅海州辛卡特礦區遙感影像蝕變異常信息提取，圈定出了鐵染蝕變異常區3處，羥基蝕變異常區3處，鐵染和羥基蝕變異常疊合區1處。然后對異常區進行地質填圖、槽探、民采工程清理工作，圈定出了礦化蝕變帶8條，礦化體23個，取得了較好的找礦成果。本文歸納出了辛卡特礦區礦床成因類型為造山型金礦。總結了地質、遙感找礦標志，提出了找礦前景。并提出了“遙感+地質”的找礦模式，在蘇丹紅海州地區的找礦有較好的示范作用。

關鍵詞：礦床地質特征；遙感解譯；綜合找礦標志；找礦前景；辛卡特；紅海州；蘇丹

中圖分類號：TP79文獻標識碼：Adoi：10.12128/j.issn.1672-6979.2024.09.007

引文格式：張朋，王軍偉，孫功偉，等.蘇丹紅海州辛卡特金礦遙感找礦方法研究應用［J］.山東國土資源，2024，40（9）：50-58. ZHANG Peng， WANG Junwei， SUN Gongwei， et al. Geological Characteristics and Remote Sensing Prospecting Method of Sinkate Gold Deposit in Red Sea State of Sudan［J］.Shandong Land and Resources，2024，40（9）：50-58.

0引言

辛卡特礦區位于蘇丹紅海州重要的Derudeb-Sinkat金礦成礦帶，帶內已在Haya、Derudeb等地發現了金礦床（圖1）。該成礦帶沿“綠巖-蛇綠巖”帶及板塊縫合帶發育，帶內斷裂構造發育，巖漿活動強烈，成礦地質條件優越，地質研究程度和勘探程度很低，但具有較好的找礦前景［1-9］。

近年來，遙感技術在礦業勘查過程中應用越來越廣泛，在金屬找礦領域的應用效果也越來越明顯，尤其是缺少地質資料且覆蓋較少的地區效果明顯［10-14］。通過對獲得的遙感影像解譯提取巖石、構造、礦物蝕變異常信息，優選出找礦有利靶區，能夠提高找礦效率，快速確定普查選區［15］。因此，利用遙感影像解譯提取遙感蝕變異常信息來指導找礦已經成為遙感地質應用的重要領域之一［16-19］。

辛卡特礦區面積大，植被和第四系覆蓋少，地形高差較大，沒有基礎地質資料作為依據，缺乏找礦信息和找礦線索，找礦難度較大，利用遙感蝕變信息解譯作為找礦指導是非常有效的。本文通過對獲取的遙感影像的預處理、去干擾剔除遙感設備、大氣、地表干擾物對遙感影像的干擾，然后通過主成分分析法提取有用的礦化蝕變異常信息［20-21］，圈定鐵染和羥基蝕變異常區，然后進行地質填圖異常查證，通過綜合研究圈定出有利找礦靶區，再進行大比例尺地質填圖、槽探等工作，開展金礦地質找礦工作，并取得了較好的找礦成果。

1礦區地質概況

礦區位于蘇丹東北部，大地構造位置處于阿拉伯—努比亞地盾（ANS）西南部、撒哈拉克拉通以東，東非大裂谷以西，海亞地塊（Haya）東北部［22］，Nakasib縫合帶邊緣。區域上出露的地層主要以前寒武紀結晶基底和第四系為主。前寒武紀結晶基底巖系以中元古代—古元古代粗粒淺色花崗質片麻巖為主，其次為新元古代變質沉積巖和變質火山巖。區域構造以一系列的縫合帶、韌脆性剪切帶和褶皺帶為主。巖漿巖廣布，大面積出露中性火山巖和花崗巖，少量脈巖穿插其中。

礦區內出露的地層主要為一套新元古代變質火山巖、變質沉積巖及新生界第四系。新元古代變質火山巖主要為基性灰綠色變質火山熔巖、變質安山質火山巖，夾雜部分層狀變質火山碎屑巖。該套巖系整體遭受區域變質作用的同時也受到造山期剪切片理化作用，但片理化不均勻，是礦區內主要的賦礦圍巖。新元古代變質沉積巖主要為灰色方解石大理巖，亦遭受了剪切片理化作用，片理化較強。區內構造較為發育，以NE向韌—脆性剪切帶為主，寬1 km，長約20 km，是區內主要的控礦構造，其次為NW向的斷裂構造，多為成礦后構造，對蝕變礦化帶有切斷作用。巖漿巖主要出露新元古代灰綠色中粗粒輝長巖、灰綠色粗粒輝石巖、新元古代肉紅色鉀長花崗巖、新元古代花崗閃長巖，大量發育輝綠巖脈、花崗細晶巖脈等脈巖（圖2）。

2遙感數據處理及蝕變信息提取

2.1遙感數據源

本文使用的數據來自Landsat8號衛星OLI遙感數據，面積約300 km2，覆蓋整個礦區。Landsat8衛星發射于2013年，由美國航空航天局發射，裝備有Operation Land Imager（簡稱OLI，陸地成像儀）傳感器及Thermal Infrared Sensor熱紅外遙感器，有9個波段的感應器，覆蓋了從紅外到可見光的不同波長范圍（表1）。

2.2數據預處理

下載的Landsat8數據經專業部門幾何粗校正后，開展數據的預處理工作。為消除傳感器自身等方面導致的輻射信息失真［23］，本文采用ENVI5.3軟件輻射定標功能對Landsant8遙感影像進行輻射定標，將數據的DN值轉換為輻射亮度值，以獲得真實的地物反射率［24］。然后對定標后的影像利用ENVI軟件FLAASH模塊進行大氣校正，消除大氣中的二氧化碳、氧氣、臭氧等成分吸收和散射的影響獲得地物真實的反射率［25］，根據礦區的坐標范圍進行裁剪，得到最終的遙感影像。

2.3去干擾處理

干擾地物通常是指具有較高反射率的地物類型，如云雪、植被、水體、陰影和第四系堆積物等［26］。本次所用的數據為OLI 遙感數據，軌道號為 171047，成像時間為2021年5月21日，該景數據中無云、雪、霧霾等覆蓋，且處于干旱地區地表無水體，主要的干擾因素為少量植被和第四系松散堆積物。本次采用高端切割法去除第四系堆積物的干擾信息，利用第四系坡積物和沖積物在Band6和Band1波段的反射率最高的特征，利用高端切割法去除干擾信息。利用比值法選擇Band5/Band4的植被指數增強植被信息，增強后的植被信息分布在高值區，能被清晰的區分開來，然后通過高端切割法去除植被干擾［27］。

2.4蝕變信息提取

與礦區金礦化有關的蝕變類型主要為硅化、絹云母化、綠泥石化、高嶺土化和碳酸鹽化，礦化類型主要為褐鐵礦化、赤鐵礦化、黃鐵礦化等。上述蝕變礦化類型均含有Fe2+，Fe3+離子或OH-，CO2-3離子，因此鐵染和黏土化是蝕變信息提取的主要目標。

根據USGS波普庫中的相關礦物的波普特征，鐵染礦物在Landsat8數據的波段曲線波普特征，在Band2、Band5波段為吸收特征，在Band4、Band6波段顯示為較高的反射特征，因此選擇Band2、Band4、Band5、Band6波段進行主成分分析提取鐵染異常。Band2、Band5為吸收帶，對應的特征向量符號相同；Band4、Band6對應為反射帶，特征向量符號相同，且Band2與Band4特征向量符號相反，選擇PC4作為含鐵染蝕變異常的主分量進行分析（表2）。

根據USGS波普特征，羥基蝕變礦物在Landsat8數據的Band2處呈現高反射率特征，在Band5處呈現吸收特征，在Band6處的反射率較高，在Band7處羥基黏土礦物出現強烈的吸收特征，因此，選擇Band2、Band5、Band6、Band7波段進行主成分分析，提取羥基蝕變信息。Band2與Band6特征向量符號相同，與Band5和Band7符號相反，選擇PC4作為含羥基蝕變信息的主分量進行分析（表3）。

為了體現蝕變異常信息強弱不同層次的呈現，采用主分量密度分割法對蝕變信息進行分級處理。一般采用“均值+K倍標準離差（σ）”劃分異常等級［28］。對于鐵染蝕變信息，取K值為1.5、2.0、2.5，對于羥基蝕變信息，取K值2、2.5、3，將鐵染蝕變、羥基蝕變分為一級異常區、二級異常區、三級異常區。對分級后的異常采用5×5中值濾波去除孤立點，消除孤立的噪聲點，得到最終的蝕變異常提取結果。

3野外異常查證與找礦成果

3.1野外異常查證

通過對遙感數據的處理，在礦區范圍內劃出羥基異常區Q1、Q2、Q3（圖3），鐵染異常區T1、T2、T3（圖4）。Q1和T1異常范圍基本重合，為羥基異常與鐵染異常重合區。本次對羥基異常區和鐵染異常區開展了野外異常查證工作。

Q1羥基異常區，分布于礦區的西北部，主要出露巖性為大理巖和變質火山巖，羥基異常分布在大理巖與變質火山巖接觸帶或變質火山巖內部韌性剪切帶中，主要見絹云母化、硅化、高嶺土化等蝕變現象，是引起羥基異常的主要原因。經踏勘見少量民采剝土工程，見黃鐵礦化、褐鐵礦化等礦化現象，取樣分析礦化顯示良好，確定為下一步工作的重點區域。

Q2羥基異常區，位于礦區的東部，主要出露巖性為新元古代中期花崗巖與新元古代早期花崗閃長巖，羥基異常主要分布于2個巖體的接觸帶上，接觸部位見絹云母化、高嶺土化等蝕變現象，是引起羥基異常的主要原因，經踏勘、取樣分析礦化顯示為零。

Q3羥基異常區，位于礦區的東南部，羥基異常零星分布，主要出露的巖性為新元古代中期花崗巖與新元古代早期花崗閃長巖，羥基異常主要分布在局部斷層內，見高嶺土化、白云母化等蝕變現象，是引起羥基異常的主要原因，經踏勘取樣分析未見明顯礦化顯示。

T1鐵染異常區，分布于礦區的西北部，主要出露大面積的大理巖和變質火山巖，鐵染異常分布在變質火山巖內部韌性剪切帶中，主要見褐鐵礦化、黃鐵礦化等礦化現象，是引起鐵染異常的主要原因。該異常區與羥基異常區Q1基本重合，確定為重點找礦工作區。

T2、T3鐵染異常區，分布于礦區的東部，主要出露大面積新元古代中期花崗巖與新元古代早期花崗閃長巖，另外有少量輝石巖、輝長巖捕擄體發育，在捕擄體的邊緣見少量綠泥石化、褐鐵礦化等現象，是引起鐵染異常的主要原因，在踏勘過程中進行取樣分析，未見明顯礦化。

3.2找礦成果

通過異常查證后，選定羥基異常Q1與鐵染異常T1重合區作為重點異常查證區開展地質草測工作。通過地質填圖草測工作，選定出西北異常區為該區的重點靶區，并在該區開展了大比例尺地質填圖、槽探及剝土工程理清編錄工作。

礦區發現礦化蝕變帶8條，賦存于剪切構造內，長500～1 000 m，寬1～5 m，走向30°～40°，傾向NW，傾角50°～70°，產狀與剪切帶產狀基本一致，8條蝕變帶大致平行展布，間距30～100 m不等，帶內共發現金礦體15個，礦化體8個，主要由地表探槽和剝土工程控制，礦體控制長20～80 m，控制斜深0～50 m，金品位1.00～16.61 g/t，平均品位2.83 g/t；厚度0.87～1.20 m，平均厚度0.98 m。礦體賦存于礦化蝕變帶內，產狀與蝕變帶產狀基本一致，嚴格受蝕變帶的控制，呈小透鏡體發育，有尖滅再現、膨脹夾縮的特點。礦體與礦化蝕變帶界限清晰，礦石巖性為黃鐵礦化石英脈、黃鐵礦化硅化變質安山巖，圍巖以硅化變質安山巖為主。礦石含有非金屬礦物為石英、絹云母、斜長石，次為碳酸鹽，含少量綠泥石、綠簾石；金屬礦物含量較少，約占礦物總量的3%，主要為黃鐵礦，次為黃銅礦，含少量自然金、銀金礦。圍巖蝕變主要有硅化、絹云母化、綠泥石化等，礦化主要為黃鐵礦化、褐鐵礦化等，硅化、絹云母化、黃鐵礦化、褐鐵礦化與金礦化關系較為密切。蝕變帶及礦體受韌性剪切帶和變質火山巖分布的控制，礦化體沿走向及傾向均未封閉。

3.3礦床成因類型分析

造山型金礦指產于區域上各個地體中、在時間和空間上與增生構造有關的脈型金礦床系列，礦床形成于增生或碰撞造山帶的匯聚板塊邊界上的擠壓和扭壓作用過程中。前寒武紀形成的造山型金礦多發育于綠巖帶和被動大陸邊緣沉積巖系內。曾國平等［29］對努比亞地盾造山帶型金礦床與膠東造山帶型金礦床進行比對，總結了努比亞地盾造山型金礦床的特征，認為東北非造山型金礦床形成時代為新元古代（約600 Ma），并闡述了蘇丹東北部造山型金礦床主要受剪切帶的控制，并且具有多級構造分級控礦的特點，蝕變巖型礦化產于韌性剪切帶的韌性區域內，石英脈型礦化則產于脆性區域內，次級斷裂構造、片理化構造對金礦體控制作用明顯；綠巖帶及其中的中—基性變質火山巖與金礦化有密切關系。張龍升等［30］對蘇丹東北部金礦床地質特征總結，認為Haweit、Galat Sufar、Qbgih、Gebeit等金礦床是造山型金礦床，造山時限界定為750～550 Ma。

礦區位于東非造山帶邊緣，Haya地塊北部Nakasib縫合帶南部，受多期次造山作用影響，經歷了一系列的弧弧拼貼、弧陸碰撞的造山運動，形成了擠壓、走滑剪切為主的構造環境。礦體的發育嚴格受NE向脆-韌性剪切帶的控制，圍巖主要為新元古代變質安山質火山巖，呈透鏡狀產出，賦礦巖性為黃鐵礦化硅化變質安山巖和石英脈；礦化有黃鐵礦化、褐鐵礦化、黃銅礦化；蝕變以硅化、絹云母化、綠泥石化、碳酸鹽化為主；原生礦石中金屬硫化物以黃鐵礦為主，含量lt;3%。與蘇丹東北部地區已知典型造山型金礦床特征相比較具有相似性，認為礦區礦床成因類型屬于造山型金礦床（表4）。

4礦區找礦標志

4.1地質找礦標志

（1）地層標志：蘇丹東北部地區大部分金礦床產于前寒武紀變質火山巖中，且對金礦床的分布有明顯的層控作用。礦區內廣泛出露的變質基性火山巖是金礦找礦重要的地層標志。

（2）構造標志：造山型金礦多受剪切帶控制，礦區內發育大型的韌—脆性剪切帶，帶內已發現多個礦化體，并對礦體的產狀有明顯的控制作用，剪切帶構造是金礦找礦重要的構造標志。

（3）蝕變礦物標志：產于韌性剪切帶內的金礦，熱液蝕變十分發育，熱液蝕變主要有硅化、絹云母化和綠泥石化等，礦化有黃鐵礦化、褐鐵礦化、黃銅礦化等，蝕變和礦化礦物是最直接的找礦標志（圖5）。

（4）礦區及周邊民采工程及砂金采礦點等采礦活動。

4.2遙感解譯蝕變標志

通過對遙感數據的處理，用主成分分析法提取蝕變異常信息，礦區內圈定出羥基異常和鐵染異常區各3處，根據鐵染與羥基的異常疊加選定出重點異常查證區，把其他異常區作為次重點區開展地質踏勘工作。根據填圖工作證實，鐵染異常和羥基異常重疊區作為重點找礦靶區可作為找礦預測的標志。

5礦區找礦前景分析

實踐證明，大多數情況下受剪切構造帶控制的金礦床特征表現為大斷裂成大礦，小斷裂成小礦的特點，特別是特大型金礦床往往位于區域構造系的復合交叉及轉折部位［25］。礦區金礦床類型為造山型金礦床，目前已發現并初步控制的礦體均賦存于剪切帶構造中，剪切帶規模較大，剪切帶內礦化蝕變較好，根據槽探及地表民采工程控制，地表礦化相對較為連續，且局部民采工程已采至50 m以下，認為該剪切帶在深部有較好的成礦潛力。

6結論

（1）通過與周邊金礦床地質特征對比，將辛卡特金礦成因類型歸于造山型金礦床。其地表發現的蝕變帶沿走向金礦化較為連續，并嚴格受NE向剪切帶和變質火山巖控制，局部民采工程已達50 m以下，礦體沿傾向向深部延續，建議向深部施工鉆探工程進一步驗證礦體的連續性，以期有更大的找礦突破。

（2）利用獲得的遙感影像對礦區進行了地質構造及巖性信息解譯，對礦化蝕變信息提取，圈定出鐵染異常和羥基異常各3處。通過遙感異常疊加并進行異常查證，發現強羥基蝕變異常與弱鐵染蝕變異常疊合區容易發現金礦化體。因此，“強羥基蝕變異常+弱鐵染蝕變異常”特征疊合區，可作找礦重點關注的遙感異常區。另外，遙感地質解譯的線性、環形構造對尋找到剪切帶、火山口是較好的金礦找礦標志。

（3）通過對礦區的遙感數據的解譯，獲取羥基和鐵染異常，然后進行地質異常查證，取得了較好的找礦效果，這種找礦模式在蘇丹東北部地區缺少基礎地質資料的勘查區找礦有較好的示范作用。

參考文獻：

［1］李杰，倪師軍，李世勇，等.蘇丹北部拜尤達地區金礦地質特征及找礦遠景［J］.地質找礦論叢，2012，27（2）：214-221.

［2］趙鵬，曲為貴，劉建權，等.蘇丹共和國哈馬迪金礦成礦規律及成礦作用［J］.新疆有色金屬，2018，41（5）：1-4.

［3］雷義均，李勇，姚華舟，等.蘇丹東部地區BIF型鐵礦和CID型鐵礦的發現及其找礦意義［J］.中國地質調查，2014，1（1）：38-45.

［4］徐九華，王建雄，向鵬，等.極富CO2流體的造山型金礦：蘇丹哈馬迪金礦［J］.巖石學報，2015，31（4）：1040-1048.

［5］薩拉（SALIH HUSSAIN SALAH MOHAMED ）.蘇丹東北部紅海山Ariab塊狀硫化物礦床地球化學特征及其找礦意義［D］.北京：中國地質大學（北京），2020.

［6］ABDALLA M S A.蘇丹東北部哈米薩納（Hamisana）剪切帶新元古代巖漿作用對努比亞地盾（Nubian Shield）的構造演化的啟示［D］.北京：中國地質大學，2022.

［7］胡建勇，楊倫，李輝，等.蘇丹東北部金礦成礦特征及成礦遠景淺析［J］.地質與勘探，2011，47（3）：505-511.

［8］楊學明，秦正永，劉文奎，等.地面高精度磁法測量在北蘇丹國哈馬迪金礦勘查中的應用［J］.地質調查與研究，2013，36（3）：233-238.

［9］鮑中義，楊立新，譚建剛，等.蘇丹哈佳吉金礦區地質特征及控礦因素淺析［J］.山東國土資源，2011，27（10）：16-19.

［10］欒卓然，呂鳳軍，趙志遠，等.干旱半干旱地區降雨對遙感地質信息提取影響研究［J］.地質找礦論叢，2021，36 （3）：342-354.

［11］韓學林，王秀芬，郭寶奎，等.山東省砂石資源遙感解譯成果及應用［J］.山東國土資源，2023，39（11）：59-68.

［12］李明波，武斌，顏丙鵬，等.新疆甜水海地區鉛鋅礦成礦遠景分析［J］.山東國土資源，2022，38（1）：13-18.

［13］王昌宇.遙感技術在地質找礦中的應用及發展前景［J］.世界有色金屬，2024（2）：96-99.

［14］葉發旺，張杰林，張川，等.航空高光譜遙感技術在鈾礦找礦中的典型應用：以新疆雪米斯坦地區為例［J］.世界核地質科學，2024，41（2）：233-249.

［15］趙玉梅，劉磊，周軍，等.新疆烏什塔拉地區遙感、化探、地質綜合找礦研究［J］.地質找礦論叢，2016，31（1）：135-141.

［16］吳志春，葉發旺，郭福生，等.主成分分析技術在遙感蝕變信息提取中的應用研究綜述［J］.地球信息科學學報，2018，20 （11）：1644-1656.

［17］魏英娟，劉歡.北衙金礦床遙感礦化蝕變信息提取及找礦預測［J］.自然資源遙感，2021，33 （3）：156-163.

［18］趙超.現代遙感技術在地質找礦中的應用［J］.中國金屬通報，2023（12）：37-39.

［19］覃旭，王振琦.現代遙感技術在青海木里礦地質找礦中的應用研究［J］.中國金屬通報，2021（3）：27-28.

［20］梁棪先，馬玉輝，思積勇.遙感技術在青海省哈得爾甘地區找礦應用研究［J］.中國錳業，2020，38（1）：70-74.

［21］張磊，孟令華，周龍濤.新疆溫泉縣牙馬特一帶1∶5萬區域地質調查的遙感技術應用［J］.山東國土資源，2020，36（3）：66-73.

［22］周佐民，李勇，劉曉陽，等.蘇丹紅海州新元古代A型花崗巖的地球化學特征及構造意義［J］.華北地質，2023，46（1）：71-79.

［23］張勝龍，劉林，耿海軍，等.基于Landsat 8 OLI數據甘肅隴南趙家莊地區遙感蝕變信息提取［J］.鈾礦地質，2020，36 （6）：535-540.

［24］賀金鑫，姜天，董永勝，等.基于Landsat 8的遼寧弓長嶺區遙感蝕變信息提取［J］.吉林大學學報（地球科學版），2019，49（3）：894-902.

［25］唐超，周可法，張楠楠，等.基于Landsat 8 OLI和ASTER數據集成和融合的礦化蝕變信息提取：以包古圖斑巖型銅礦為例［J］.地質科技情報，2018，37（6）：211-217.

［26］茹菲娜·阿力木江，陳川，高玲玲，等.基于ASTER數據的遙感蝕變信息提取：以西天山穆龍套地區為例［J］.地球學報，2022，43（2）：235-245.

［27］別小娟，張廷斌，孫傳敏，等.藏東羅布莎蛇綠巖遙感巖礦信息提取方法研究［J］.國土資源遙感，2013，25（3）：72-78.

［28］杜曉川，婁德波，張長青，等.基于GF-5、Landsat8與GF-2遙感數據的蝕變信息提取研究：以四川寧南鉛鋅礦集區為例［J］.礦床地質，2022，41（4）：839-858.

［29］曾國平，向文帥，吳發富，等.東北非與中國膠東造山型金礦對比及對非洲找礦勘查的啟示［J］.地質通報，2022，41（1）：48-59.

［30］張龍升，張志強，林芳梅，等.蘇丹東北部哈維特金礦地質特征及控礦構造研究［J］.地質與勘探，2022，58（1）：207-216.

［31］王明明，羅清威，張克川，等.坦桑尼亞穆索馬-馬拉綠巖帶型金礦遙感、物探、地質綜合找礦研究［J］.地質與勘探，2022，58（5）：1117-1127.

收稿日期：20240530；修訂日期：20240725；編輯：陶衛衛

作者簡介：張朋（1984—），男，山東淄博人，高級工程師，主要從事地質礦產勘查工作；E-mail：270770617@qq.com

*通訊作者：王軍偉（1978—），男，山東萊蕪人，工程師，主要從事地質礦產勘查工作；E-mail：lywhwjw@163.com