綜合地球物理在荒漠覆蓋區隱伏礦床預測與定位中的應用：以新疆拉伊克勒克銅多金屬礦床為例

嚴加永， 孟貴祥*， 呂慶田， 鄧震， 祁光， 湯賀軍，趙金花， 薛融暉， 王栩

1 中國地質科學院， 北京 100037 2 中國地質調查局中國地質科學院地球深部探測中心， 北京 100037

0 引言

對傳統的礦產勘探空間進行勘查，已滿足不了人類對能源日益增長的需求，必須拓展新的找礦空間.當前，全球大部分的礦產資源勘查工作的目光已向大陸深部和覆蓋層(包括沙漠戈壁、森林、沼澤湖泊等淺覆蓋區)下方聚焦(王欽軍等，2017；崔敏利等，2010).覆蓋區一般指第四紀以前的基巖受沖積、洪積、冰積及風積等沉積作用，最終被冰雪、水體、沼澤、植被、第四紀松散沉積物廣泛掩蓋的地區(孫凱，2018).在大面積的沉積物覆蓋區下，如加拿大、俄羅斯等毗鄰北極的冰川覆蓋區，澳大利亞、非洲的沙漠或紅土覆蓋區等，孕育了北美、格陵蘭、非洲—阿拉伯等全球重要的成礦區帶(梅燕雄等，2009).中國境內廣泛分布大面積的中新生界覆蓋區，僅戈壁、沙漠覆蓋區所占的總面積就高達128萬km2，占全國陸地總面積的12.3%(圖1)，但在戈壁、沙漠覆蓋區內發現的大中型礦床僅為全國大中型礦床總數的2%(申元村等，2013).為進一步拓展我國的找礦空間，在借鑒國外在覆蓋區內進行找礦的眾多成功實踐經驗(Cao et al., 2009；Hosseini et al., 2016；Shan et al., 2013)基礎上，結合我國礦產資源開發的基本現狀，表明在戈壁、沙漠等中新生界覆蓋區開展礦產勘查和綜合預測工作無疑是當前形勢下最佳的選擇之一.

圖1 中國裸地及荒漠覆蓋區示意圖Fig.1 Schematic map showing bare land and desert Gobi covered areas in China

在荒漠戈壁覆蓋區，由于沒有巖石露頭或者露頭出露范圍很小，地表地質的觀察無法獲取更多的含礦信息，極大的增加了地質理論預測的不確定性，傳統的化探方法也難于穿透覆蓋層(黃桂珍，2014；靳職斌等，2014；成秋明，2012).地球物理探測雖然能夠獲取地下的物性結構信息(郝興中等，2013)，但荒漠戈壁區內存在較厚的鹽堿殼、地表干燥等不利因素嚴重影響了地球物理方法特別是電磁法數據的采集，且由于缺少地表地質等已知信息的約束，嚴重制約了地球物理資料的地質解釋，導致在覆蓋區范圍內發現礦床的難度更大(劉光鼎和郝天珧，1995).

急需加強對荒漠戈壁覆蓋區隱伏礦產勘查技術的研究(汪青松等，2021)，探索行之有效的技術方法體系，為保障我國礦產資源的供給提供有力的技術支撐.

瓊河壩地區位于中國新疆哈密地區伊吾縣，東與蒙古人民共和國接壤.構造上位于中國新疆東準噶爾盆地東緣，處于阿爾泰造山帶與天山造山帶交匯部位.該區位于古亞洲巨型成礦帶中，該成礦帶已經成為國際資源勘查的一個熱點，近年已經在蒙古(Perello et al.,2001)、哈薩克斯坦、俄羅斯等國發現8個千萬噸級的銅礦(Seltmann and Porter,2005)，但在中國境內一直沒有大的突破，國內的許多學者對其進行綜合對比研究認為，在瓊河壩地區的找礦潛力十分巨大(王登紅等，2009；楊富全等，2010).近年來，隨著對瓊河壩地區找礦工作的不斷深入，相繼發現了寶山夕卡巖鐵礦、北山構造蝕變巖金礦、綠石溝矽卡巖銅礦、蒙西斑巖銅礦、和爾賽斑巖銅礦、銅華嶺斑巖銅礦、瓊河壩斑巖銅礦及桑南斑巖銅礦等礦床(孟貴祥等，2016).雖然上述礦床的規模并不大，但其表明瓊河壩地區具備尋找大型礦床的前景和成為大型礦集區的潛力(屈訊等，2009)，但該區廣泛分布荒漠戈壁，大大增加了找礦工作的難度.本文采用500 m×100 m網度的地面重力和磁力數據，以多尺度邊緣檢測和三維物性反演技術為主，進行斷裂構造系統劃分和隱伏巖體識別工作，厘定主要斷裂系統的分布，圈定隱伏巖體的分布范圍及三維形態.在此基礎上，結合地球化學資料，進一步的開展隱伏礦找礦靶區預測.通過1∶10000到1∶2000比例尺綜合地球物理探測，成功定位隱伏銅鐵礦體，并對拉伊克勒克靶區進行鉆探查證，獲得銅金屬量118萬噸，評價為大型銅多金屬礦，實現了瓊河壩地區淺覆蓋區的找礦突破.

1 瓊河壩覆蓋區隱伏礦預測

1.1 成礦信息提取技術方法

目前在瓊河壩地區發現的礦床主要有三類：斑巖型銅金礦(蒙西、和爾賽、桑南、瓊河壩和銅華嶺銅礦)、矽卡巖型鐵銅礦(寶山磁鐵礦和綠石溝鐵銅礦)、構造蝕變巖型金礦(北山金礦).這三類礦床的形成均嚴格受巖體和斷裂構造控制，查明斷裂構造和巖體的分布特征及三維形態是尋找這三類礦床的關鍵(王曉地等，2006；董連慧等，2009；孟貴祥等，2015).另一方面，由于該區地質研究程度相對較低，且荒漠戈壁覆蓋區分布范圍較廣，地表巖石露頭較少，對隱伏的巖體和斷裂系統缺乏深入研究，嚴重制約了瓊河壩地區找礦預測研究工作的開展.因此，查明區域斷裂構造和隱伏巖體的分布對找礦預測有著重要的現實意義，蝕變及礦化也是找礦預測的重要線索.本文將瓊河壩地區成礦要素與信息歸納為：構造、巖體和礦化蝕變(周文月等，2021).針對具體成礦信息的特征，設計了如圖2所示的技術路線.

圖2 瓊河壩地區成礦信息提取與隱伏礦預測技術路線Fig.2 Technology route of metallogenic information extraction and concealed ore deposit forecasting in Qiongheba area

要提取礦集區尺度特別是覆蓋區構造、隱伏巖體控礦地質體形態特征等成礦信息，必須要有覆蓋面廣、精度高的數據，同時，還需要有針對性的數據處理、反演和可視化方法(俞貴平等，2020).雖然反射地震具有較大的探測深度和較高的垂向分辨率，但受其施工成本的影響，現階段在金屬礦礦集區探測中仍以二維剖面探測為主(呂慶田等, 2010)，很少開展大面積的三維面積性探測，因此，獲得的結果主要是建立礦集區的“骨架”結構(呂慶田等, 2014).而重力和磁力方法由于施工成本相對較低，完全可以實現高密度的數據全覆蓋.目前，除西藏和青海的局部地區以外，大比例尺的航磁和地面重力數據已經基本覆蓋了我國大部分地區(熊盛青等，2016).對重磁數據進行有效的處理，是現階段實現礦集區尺度隱伏巖體和斷裂厘定最有可能成功的途徑.

瓊河壩地區成礦信息可以從結構和屬性兩方面開展.結構主要指斷裂構造和巖體結構特征，以1∶50000高精度地面重力和磁力資料為基礎，采用多尺度邊緣檢測技術，結合重磁場特征分析及導數變換等技術，識別和劃分斷裂構造.以重磁三維反演和巖性識別技術為主，輔助多尺度邊緣檢測技術，圈定和識別巖體及隱伏巖體的平面投影位置和三維形態，厘定不同基性程度巖體的分布范圍和三維形態.屬性主要指礦化和識別信息，采用1∶50000大功率激電中梯和地球化學測量，從極化率特征和地球化學異常組合綜合分析，提取礦化和蝕變信息.在此基礎上，結合地質資料，開展了成礦預測工作，圈定找礦靶區.

1.2 斷裂構造厘定和巖體識別

由于地質體邊界兩側一般存在密度或磁化率差異，重磁異常在地質構造邊界附近表現為梯度變化帶，這些梯度帶的實質就是場源的邊界，因此，重磁場構造信息提取的主要方法是從重磁異常圖中提取場源邊界(張壹等，2015).通常的做法是對重磁異常中的場源邊界信息進行增強，然后利用某種邊緣檢測的手段確定邊界位置(孟小紅等，2012;Guo et al.,2015).重磁異常邊緣增強檢測的方法有很多種，如：斜導數法、斜導數水平梯度法、Theta圖法等(鄭強等，2019).每種方法的原理和應用的前提條件不盡一致，相同的一組重磁數據，用不同方法處理后的結果也有所不同.因此，有必要對不同方法的檢測效果進行對比分析，優選出效果較好的方法.通常，在不同地區或不同地質背景的重磁異常解釋過程中，也需要對各種檢測增強方法獲取的結果進行對比分析，篩選出某一種與實際地質情況吻合較好的方法，結合已有資料綜合分析，再給出合理的地質解釋(張興東等，2018). 湯井田等(2019)、胡雙貴等(2019)將重力梯度張量曲率應用到重力數據的邊界識別中，通過理論模型和實際數據詳細分析和比較了各種曲率在重力邊界識別中的應用效果，認為局部坐標系下所計算的高斯曲率進行邊界識別能夠較好的圈定地下地質體的邊界.嚴加永等(2011，2015)通過模型試驗對比了應用較為多的總水平梯度法、斜導數法、Theta圖法和多尺度邊緣檢測法(簡稱Worms法)，發現多尺度邊緣檢測法比較適合大尺度重磁數據的邊緣提取，并在長江中下游成礦帶、銅陵礦集區構造信息提取中發揮了作用.

瓊河壩地區的重磁多尺度邊緣檢測結果如圖3所示，重力多尺度邊緣檢測以線性構造為主，反映了構造斷裂、地層和巖體邊界.化極磁異常多尺度邊緣檢測結果顯示低延拓高度大多呈現環形異常，高延拓高度大多呈現線狀異常，此是由于淺部磁異常的分布廣泛形成了較多的小型環形異常.因此，小型環形異常主要反映了淺部磁性體的邊界，較大的環形異常可能反映了中基性巖體的邊界，而規模較大的線狀構造則反映了斷裂構造.本文從重磁異常劃分斷裂構造的主要依據是重磁多尺度邊緣檢測，同時結合重磁場平面特征和地質信息，開展斷裂構造的劃分，具體依據如下：重磁多尺度邊緣檢測：在同一圖幅內，將不同延拓高度檢測結果疊加到一起(圖3)，采用不同顏色、大小的線型表示不同深度尺度的信號邊界形跡，不同深度邊界信號形跡在分布位置上所體現的相似性即可獲知該構造形跡的發育深度及傾向特征：線束越密集表示邊界構造切割深度越大，線束稀疏則表示其切割發育的深度較淺；線束組合越寬表示該構造傾向越緩，反之則表示該邊界構造發育產狀較陡、傾角較小.

圖3 瓊河壩地區重力多尺度邊緣檢測(a)與化極磁異常多尺度邊緣檢測(b)結果線束從藍色到紅色反映構造深度從淺到深的變化.Fig.3 Multiple-scale edge detection result of Bouguer gravity (a) and reduction to the pole magnetic anomalies (b) overlapped on geological map of Qiongheba areaThe color lines from blue to red reflect the change of structural depth from shallow to deep.

按照上述的原則，以重磁多尺度邊界檢測為主要手段，分別根據重力和磁力多尺度邊緣檢測結果劃分了斷裂構造，然后根據地質填圖和其他相關信息，對斷裂系統進行了修正，得到了瓊河壩地區斷裂構造簡圖(圖4).根據構造特征分析，三條骨干斷裂將本區劃分為四個大的構造分區(圖4).F1斷裂由南北兩段追蹤而成，走向從北部的北西向逐漸轉換為北西西向，北部產狀較陡，到中部逐漸變為北東傾向，這可能反映了研究區中部的大巖基是從北東方向侵入過來的.F1斷裂將研究區北部分為I和II兩個差別明顯的構造分區，I區構造方向以北東和近東西為主，II區則以北西方向為主，I區巖體侵入時代以志留系為主，II區巖體以石炭紀為主，I區出露地層為奧陶、石炭、二疊和侏羅，奧陶紀地層對應為高重力異常，而II區主要出泥盆系地層，表現出瓊河壩中央高重力異常帶.F2斷裂和F3斷裂將南北分為II、III和IV區，F2斷裂是橫亙研究區中部一條大斷裂，無論從重力、磁力還是地形、地質分析都可清晰看到這條斷裂，斷裂走向為北西西，總體朝北東反向傾，西端產狀較緩，東段產狀較陡.F3走向北東東，傾向朝東，這是一條推測斷裂，在磁力多尺度邊緣檢測上顯示明顯.F3及F2斷裂將研究區南部劃分為III和IV區，III區斷裂以北西和北東兩組方向為主，IV區斷裂則以北西為主，這些斷裂與F2小角度交匯，說明F2斷裂具有右行特征.III中巖體以基性為主，侵入時代多為石炭紀，IV區巖體以中酸性為主，侵入時代多為志留系，且多為隱伏巖體，推測存在隱伏的大巖基.

圖4 瓊河壩地區斷裂識別與構造分區Fig.4 Fault identification and structural division in Qiongheba area

1.3 隱伏巖體識別

嚴加永等(2009，2014a，b)在長江中下游成礦帶銅陵礦集區、廬樅礦集區和沙溪礦床用重磁三維反演研究巖漿巖的三維形態，根據這些實例，探索了利用重磁三維反演方法識別巖漿巖體的基本原理和方法.首先對重磁數據進行處理：對布格重力異常采用位場分離技術分離區域場與局部場，提取研究目標深度內的剩余重力異常；對磁異常數據先進行化極處理，再采用位場分離技術分離區域場與局部場，提取目標深度內的磁異常信息；然后將剩余重力和磁異常作為反演數據，采用相同的網格剖分，進行三維物性反演，得到地下半空間磁化率和密度差三維數據體；最后，根據研究區巖體與圍巖的磁性和密度參數差異，結合其他地質信息，確定圈定巖體的磁化率和密度差閾值，滿足磁化率和密度差的閾值組合部分即可視為是巖漿巖體的反映，通過三維可視化的交互分析，可以確定巖漿巖體三維空間形態.本次工作采用加拿大大不列顛哥倫比亞大學開發的三維物性反演軟件UBC Mag3D和Grav3D，軟件的核心算法基于Li和Oldenburg(1996，1998)提出的重磁反演算法，該軟件提供了開展先驗信息約束的接口，可以開展帶先驗信息約束反演(Williams，2008).

首先對重力和磁力進行匹配濾波，提取10 km以淺的局部重磁場，對重磁數據分別反演，獲取三維密度數據體和磁化率數據體，然后根據各類巖體對應的不同密度和磁化率組合，對密度體和磁化率體交集運算，滿足交集的部分即視為是對應類型的巖體(Yan et al.，2019).圖5a為中基性巖體三維分布特征，圖5b為酸性巖體三維分布特征，由于中基性巖體具有一定磁性和較高的密度，與圍巖的物性差異明顯，圈定的結果與地表實際情況吻合度高.而酸性巖體無磁性或磁性較弱，密度較低，與圍巖的物性差異相對較小，所以圈定效果與實際情況有一定偏差，如和爾賽—銅華嶺—拉伊克勒克一帶的隱伏巖體未能有效的刻畫出來，但總體來看，圈定的巖體與地表出露巖體實際情況基本吻合.在三維可視化軟件中，通過旋轉、切片、調整透明的等操作，進行人機交互分析巖體的三維形態，以及巖體與礦床、地層之間的關系，極大的方便了地質解釋工作的開展.

圖5 基于重磁三維反演識別的巖體分布立體圖(a) 中基性巖體； (b) 酸性巖體，底圖為地質圖.Fig.5 3D perspective diagram (view from SE to NW) of intrusions derived from gravity and magnetic 3D inversion(a) Intermediate-basic intrusion; (b) Acid intrusion. The base map is geological map.

1.4 蝕變信息提取

地球化學異常直接反映了蝕變礦化的可能，而極化率間接的反映了地質體蝕變特征.我們利用新疆物化探大隊采集的1∶50000激電中梯面積性數據進行處理和分析.由于1∶50000激電測量由于供電極距大，測量范圍大，參數設置不一定能完全合乎工區的各類地表環境，所以極化率受近地表地質條件影響較大，在低阻區容易出現測不準的情況，造成出現畸變點甚至假異常，而在高阻區往往會出現無法改善接地條件、供電不足導致一次電位太低，也可能出現畸變點，給區域解釋帶來不便，為突出區域極化率的分布特征，對原始極化率進行滑動窗口濾波處理，獲得相對清晰的區域極化率分布圖(圖6).

圖6 瓊河壩地區極化率異常及編號Fig.6 Induced polarization chargeability anomalies and number in Qiongheba area

總體來看上強度較大或背景較高的極化率場分布在研究區北部，一般呈面狀分布，規模較大，這種特征主要與普查區北部分布的奧陶系荒草坡群有密切的關系，由于該地層受到不同時期的侵入巖體的熱液活動使其硫化物(主要黃鐵礦)富集所致.其他多數地區極化率背景不高，形成的異常規模較小，這種特征正是反映研究區巖性主要以一套海相火山巖、火山噴發巖以及分布較為廣泛的中酸性侵入巖為主的地質環境.

根據異常特征，可劃分出編號IP-1到IP-17共17個主要極化率異常，可分為以下幾種類型：

(1)東北部與荒草坡群地層密切相關的IP-1-IP-4異常.該異常區極化率異常走向以北東東向為主，異常強度大，多呈面狀分布.異常帶出露的地層以奧陶系荒草坡群為主，巖性與以小面積的下泥盆系托讓格庫都克下段為輔，巖性為灰綠色-淺紫紅色玄武巖和中酸性的凝灰巖.后期侵入巖體的期次較多，主要為石炭紀的黑云母二長花崗巖、黑云母花崗閃長巖，志留世第一侵入期次的石英閃長巖、第二、第三、第四侵入期次的英云閃長巖等.強極化率背景上疊加了礦(化 )致異常，蒙西銅礦、201金礦、桑南銅礦和瓊河壩銅礦位于IP-1中，鑫源鐵礦位于IP-4異常中.

(2)與泥盆紀地層出露區的異常(IP-7、IP-10).這幾個極化率異常位于泥盆系火山巖出露區，對應中部的高重力異常區，異常呈面狀分布，這可能與泥盆紀火山巖中廣泛發育的黃鐵礦化、褐鐵礦化有關，礦化富集的地段形成礦床，如IP-10中的北山金礦、IP-7中的寶山鐵礦、綠石溝銅礦，IP-9中礦化點分布較多，如瓊東銅礦、瓊河壩鐵礦等.

(3)與巖體、隱伏巖體對應的極化率異常(IP-6、IP-8、IP-12、IP-13、IP-14).這些異常與巖體(隱伏巖體)空間關系密切，可能與對應巖體中硫化物分布廣泛有關，此類異常區是尋找斑巖銅礦的有利地段.IP-5和IP-6異常經過新疆物化探大隊的查證，發現了多處孔雀石化、黃鐵礦化，實施的鉆孔未發現工業規模礦體，但黃鐵礦廣泛分布，說明這兩處巖體中熱液活動強烈，在合適部位仍有可能富集成礦.IP-8地表為石炭系黑石頭組，但異常中心部位為二長巖，推測極化率與后期巖漿活動有關，且成為多條斷裂的交匯部位，也有組合的化探異常疊加，是找礦的有利地段.IP-11對應為賽北巖體，雖然目前還沒見到工業品位礦體，但仍有良好的找礦前景.IP-12為新發現的拉伊克勒克銅多金屬礦區，其東段為第四系覆蓋，西端出露泥盆系地層，目前的驗證工作表明，此異常為硫化物引起，局部地段銅礦化富集成礦，由于該異常規模較大且排除了碳質等其他因素的干擾，該異常具有良好的找礦前景，初步推測能形成中型銅(鉬、鋅)礦.IP-13為銅華嶺銅礦，IP-14為和爾賽銅礦，這兩個礦已經完成了初步勘探，求出了儲量.

(4)與含碳火山巖或煤系地層有關的異常(IP-11、IP-16、IP-17).此類異常由石炭系和侏羅系中的煤系地層和碳質引起，IP-16為橫干河北的煤礦化引起，IP-17為侏羅溝附近的煤礦化所引起.此類異常主要對應石炭系那林卡拉組中酸性火山熔巖夾中酸性沉凝灰巖，以及侏羅紀砂巖、泥質砂巖等.石炭系那林卡拉組地層中的極化率異常相對較少，多數異常分布在侏羅紀地層之中，從目前發現的礦點看普查區侏羅系地層除煤礦以外，僅有菱鐵礦點分布其中，因此這些激電異常對尋找金屬礦床的意義不大.

(5)與明礬石相關的異常(IP-5、IP-9、IP-15).此類異常與明礬石及硫鐵礦中的大量黃鐵礦化關系密切，如IP-15淖毛湖北山明礬石礦地表褐鐵礦化蝕變強烈，巖石裂隙中黃鐵礦化廣泛發育.本區激電極化率異常特征十分明顯，強度較高的極化率異常有四類主要源體引起：①黃鐵礦化普遍的奧陶系火山凝灰巖；②石炭紀黃鐵礦化發育的中酸性火山巖(明礬石礦點)；③中泥盆統北塔山組局部的含碳泥質火山片巖；④石炭—侏羅系含碳地層及含煤系地層.這其中，前兩類地質因素值得重視，在奧陶系火山巖中，與侵入斑巖體密切相關的蒙西大型銅礦；在與火山—侵入巖密切相關的北山明礬石礦附近綜合化探異常顯示有進一步尋找隱伏銅金礦前景.中等強度的激電異常影響因素較多，如北山金礦、綠石溝銅礦及和爾賽、銅華嶺及拉伊克勒克附近的激電異常，引起異常的原因多是與成礦有直接聯系的黃鐵礦化密切相關.

金屬因子(極化率與電阻率的比值)與地質情況也有一定的對應關系(圖7).受碳質影響，金屬因子最強(大于2.82)的為IP-16和IP-17對應的石炭系和侏羅系中的極化率異常，金屬因子第二強(1.4～2.28)的多為覆蓋區，這是由于覆蓋區電阻率低造成的，此級別金屬因子還有一部分是隱伏斑巖銅礦，如拉伊克勒克礦區、蒙西銅礦等；金屬因子在1～1.25的區域中一部分與礦化巖體、礦床關系密切，如銅華嶺、和爾賽、瓊河壩銅礦等地；金屬因子在0.85以下的區域反映了蝕變礦化均不強烈的巖體.從金屬因子分析，在瓊河壩地區尋找斑巖及熱液型硫化物金屬礦的有利地段主要在1～1.25及1.4～3之間范圍之間.

圖7 瓊河壩地區激電中梯金屬因子圖Fig.7 Metallic factors of central gradient array induced polarization in Qiongheba area

1.5 靶區預測

找礦靶區預測主要是根據推斷劃分的斷裂構造、圈定的巖體和化探異常，并結合已知的地質資料綜合分析研究，進行找礦預測(王巧云等，2020).在瓊河壩鐵銅金礦集區，矽卡巖-斑巖型銅礦、矽卡巖型鐵銅金礦及淺成低溫熱液構造蝕變巖型金礦是本區最主要的礦床類型.孟貴祥等(2021)對寶山磁鐵礦、蒙西斑巖型銅礦和北山構造蝕變巖金礦進行地球物理、地球化學和礦床地質解剖，總結了瓊河壩地區這三類礦床的地球物理和構造標志.

斑巖銅金礦找礦標志：(1)與隱伏巖體關系密切，物探異常組合具有明顯的高極化、高重力、中等電阻率、中等磁異常的“兩高、兩過渡帶”特征.(2)構造標志：北西向、東西向主干斷裂的次級斷裂破碎帶是含礦斑巖侵入的通道.(3)化探異常標志：礦床指示元素組合：以Au、Cu、Mo元素異常為主的Au、Cu，Au、Cu、Zn，Cu、Mo、Au和Cr、Co、Ni綜合異常組合.

矽卡巖型鐵銅金礦找礦標志：(1)巖體邊部的高磁異常，在三維反演結果表現為高磁性高密度地質體.(2)地面物探組合特征：區內3000～4000 nT以上的磁異常主要是磁鐵礦所引起，1000～3000 nT的磁異常多為玄武巖及玄武質凝灰巖引起，部分2000～3000 nT磁異常套合重力低緩異常部位是深部可能存在隱伏盲礦體所引起，需要根據地表出露的巖石加以區別.

構造蝕變巖型金礦找礦標志：(1)區域上北西向主干深大斷裂附近的次級北東、北西向、近南北和近東西向構造蝕變帶是直接找礦標志.(2)巖漿巖標志：隱伏的花崗閃長巖、閃長巖是間接找礦標志，重磁三維反演識別的中酸性巖體邊部是重要找礦地段.地面物探組合異常是高極化、低阻異常、弱-負磁異常及位于重力梯級帶中，即“一高兩低一過渡”異常組合是尋找該類礦床的物探標志.(3)化探標志：分散流Au、Ag、As、Co、Cu、Zn異常分布范圍指示礦化分布范圍.(4)地層標志：中泥盆統北塔山組，特別是隱伏巖體之上的第二亞組(D2b2)地層是間接找礦標志.

根據上述標志，結合本文給出的構造、巖體分布，圈定了7處找礦靶區(圖8)，順通靶區：斑巖型銅礦；壩西靶區：斑巖型銅礦；綠石溝—寶山靶區：矽卡巖型鐵銅礦；橫干河靶區：矽卡巖型鐵銅礦；北山外圍靶區：構造蝕變巖型金礦；拉伊克勒克靶區：斑巖型-矽卡巖型鐵銅礦；獅子溝靶區：構造蝕變巖型金礦.

圖8 瓊河壩地區巖體、構造與找礦靶區預測圖Fig.8 Intrusions，faults and prospecting targets in Qiongheba area

2 拉伊克勒克隱伏礦床定位

2.1 隱伏礦定位與勘查技術

拉伊克勒克靶區(圖8中6號靶區)，地表絕大部分為第四系覆蓋，根據巖體識別結果和斷裂劃分結果，該區下面存在隱伏的偏酸性巖體，還有2條北西南東向斷裂發育，并且有地球化學組合異常和高極化率異常，“構造+巖體+礦化蝕變”信息俱全，是尋找矽卡巖、斑巖型礦床的有利靶區.由于拉伊克勒克靶區大部分地段為第四系覆蓋，常規的方法難于進一步甄別礦在何處，為定位隱伏礦體和擴大礦床規模，我們采用圖9所示的技術路線開展進一步的工作.

圖9 拉伊拉克覆蓋區隱伏礦定位技術路線Fig.9 Location technology route for concealed ore deposits in Layikeleke cover area

2.2 隱伏礦發現—預查

預查階段前期主要任務是圈定找礦重點地段，為獲取足夠的信息，我們在拉伊克勒克靶區布設了1∶5000的大比例尺重力、磁力和1∶10000激電面積測量(中間梯度裝置)，圈定出異常后，在異常區進一步加密測線，局部地段加密到1∶2000比例尺.通過重磁三維反演、結合地球化學異常，圈定出找礦重點區段.然后在再部署電磁測深、激電測深等工作量，結合重磁交互正演等手段，進行定量(半定量)分析，確定礦體可能賦存的位置，實施鉆孔驗證.

(1)矽卡巖鐵銅礦的發現

拉伊克勒克靶區地表為第四系覆蓋，根據巖體識別和斷裂劃分結果，該區下面存在著隱伏的偏酸性巖體和發育有兩條北西—南東向斷裂.開展1∶5000的重力、1∶2000的磁力測量和1∶10000大功率激電中梯面積性測量，并在異常組合有利地段部署了音頻大地電磁測深和對稱四極激電測深，對靶區開展查證工作(圖10).

圖10 拉伊克勒克礦床極化率異常、地球化學異常與礦區地質圖1—第四系；2—新近系上新統；3—石炭系下統姜巴斯套組；4—泥盆系中統北塔山組；5—泥盆系下統托讓格庫都克組；6—英云閃長斑巖；7—閃長巖;；8—輝綠巖脈；9—英云閃長玢巖脈；10—礦床點；11—地球化學異常；12—鉆孔.Fig.10 Map showing geology, induce polarizability anomalies and geochemical anomalies in the Layikeleke deposit1—Quaternary; 2—Neogene poliocene; 3—Jiangbasitao Formation of Lower Carboniferous series; 4—Beitashan Formation of Middle Devonian series; 5—Tuorangekuduke Formation of Lower Devonian series; 6—Quartz mica dioritic porphyry; 7—Diorite; 8—Diabase dikes; 9—Quartz mica dioritic porphyrite vein; 10—Deopsit point; 11—Geochemical anomalies; 12—Borehole.

對磁異常進行化極處理后獲得化極磁異常(石磊等，2012)，異常總體呈北西走向，局部高磁異常為北東向，磁異常最高為12000 nT.為提取淺部密度體信息，對重力數據進行了匹配濾波，獲得了剩余重力異常，重力異常形態與磁異常基本一致，說明二者為同源異常，高磁異常與高重力異常位置基本吻合.大功率中梯激電獲取的視極化率異常現了一處主體異常和小異常，異常為北西向，由于背景極化率很低(0.3%)，0.7%以上的極化率即可視為高異常，推測極化率異常為金屬硫化物所致.中梯視電阻率走向也為北西向，由于地表覆蓋較厚，整體電阻率較低，最高僅400 Ωm，高電阻率反映了地層或巖體的隆起.高磁高重地段的異常組合為“三高一中”，即高重力異常、高磁異常、高極化率異常和中等電阻率異常，極高的磁異常反映地下肯定有強磁性物質的存在，加之高重力異常，高極化率異常，通過綜合分析排除了基性巖體的可能，推測該異常組合可能由磁鐵礦和金屬硫化物礦體所引起.重磁三維反演揭示了該組合異常下方存在高磁性高密度的地質體，各類異常位置基本一致(圖11).根據上述物探異常特征和重磁三維反演結果，部署了zk1號鉆孔，開孔揭穿12 m厚的第四系覆蓋層之后即發現了矽卡巖型富銅鐵礦體，礦體連續視厚度達46 m，為該隱伏礦床的發現打開了局面.

圖11 拉伊克勒克礦床綜合地球物理異常三維顯示圖(修改自嚴加永等，2017)Fig.11 Three-dimensional display of comprehensive geophysical anomalies in the Layikeleke deposit(modified from Yan et al.,2017)

(2)斑巖銅礦的發現

在zk1孔發現夕卡巖銅鐵礦后，可初步確定極化率異常與成礦關系密切.通過對極化率的進一步分析，經激電中梯等測量發現了一個形態較規整的極化率異常，它以1%等值線劃分的異常大致呈北西向帶狀展布，長度大于5 km，異常寬度500～1000 m，異常形態不規則，在異常中部有南北分叉現象(圖10).該極化率異常西北段對應HS-7號化探異常，該異常為Cu、Au、Zn、As、Mn、Co等多元素的綜合異常，尤其以Cu、Au、Zn較為富集.該激電異常對應為低磁異常和低重力異常，局部地段為高極化率、低重力、低磁力和低電阻率的異常組合，故推測為由礦化巖體所引起.2011年決定對該極化率異常進行鉆探驗證，首先在垂直異常走向的零線上，布置了zk2鉆孔(圖10)，在40 m厚的第四系和泥盆系蓋層之下，查明存在一厚達近千米強蝕變似斑狀英云閃長巖體，并在其中發現了較好的銅(鉬)礦化，揭開了拉伊克勒克大型銅礦的面紗.

2.3 礦床規模擴大-普查詳查

進入普查、詳查階段后，為擴大礦床規模，進行合理部署鉆孔.在該階段，我們通過預查鉆孔巖性物性分析，將勘探線剖面與地球物理探測結果進行綜合對比分析，優選出有效的方法技術組合，為普查、詳查的鉆孔布設提供依據.

首先我們對比了激電中梯、測深與礦體的對應關系.在激電中梯面積測量的基礎上，在276線開展了對稱四極激電測深試驗(測線位置如圖10所示)，發現激電測深高極化率中心與中梯激電高極化率中心偏差300 m(圖12)，在中梯極化率異常中心和北側實施了276-1和276-2鉆孔.276-1鉆孔見到了不同程度的銅、鉬礦化體，但礦體連續性較差，未能圈出較厚礦體(表1).276-2鉆孔見到泥巖等圍巖地層，云英閃長巖范圍較小，只有零星蝕變和礦化，為發現具有工業價值的礦體.激電測深發現ZK276-1孔雖然在中梯激電異常中心，但在測深極化率斷面上明顯極化體的外側，ZK276-2則離極化體更遠，故未鉆遇礦體，分析認為這是ZK276-1和ZK276-2孔未見連續礦體的主要原因.依據激電極化率測深成果，在極化體中心部署的ZK276-3孔揭露了9層銅礦體(表1)，其中最厚一層厚達近157 m，銅平均品位大于0.3%，表明激電測深可以揭示厚覆蓋下的隱伏礦體，而中梯激電反映的為覆蓋較淺地段的極化體，同時，也說明了激電測深結果可作為擴大礦體規模階段鉆孔部署的重要依據.

圖12 276線綜合地球物理探測剖面(a) 重力異常剖面; (b) 磁異常剖面; (c) 激電中梯剖面; (d) 高頻大地電磁電阻率斷面; (e) 對稱四極測深極化率斷面.Fig.12 Integrated geophysical profiles of Line 276(a) Gravity anomaly; (b) Magenetic anomaly; (c) Gradient induced polarization; (d) High frequency magnetotelluric resistivity; (e) Symmetric four-pole IP sounding.

表1 276線鉆孔編錄表Table 1 Logging data of boreholes in Line 276

為進一步指導鉆孔部署以擴大礦床規模，垂直圖10中的中梯高極化異常帶部署了對稱四極激電測深，完成了421個激電測深點.圖13為拉伊克勒克鐵銅礦區激電測深獲取的極化率和電阻率經三維插值后獲取的礦區極化體和電阻率體的三維空間分布示意圖，為詳查、勘探階段鉆孔部署提供了直接依據.鉆探驗證結果表明，見礦較好的地段為高極化體南西側、高電阻率體北東側，也即二者重合部位是目前發現礦體較好的地段(圖13).

圖13 拉伊克勒克礦區高極化體(a)和高電阻率體(b)與鉆孔三維空間特征Fig.13 Three-dimensional characteristics of Layikeleke deposit with (a) high-polarization body and (b) high-resistivity body

根據綜合地球物理方法資料的綜合分析下，在拉伊克勒克礦區共部署施工鉆孔53個，有效進尺共計22446.4 m，發現并評價為大型銅多金屬礦.斑巖銅鉬礦主體產于英云閃長巖體中，少部分延伸于巖體外的火山巖地層及黑云母二長花崗巖體內.據工程控制目前連接礦體的特點，斑巖型銅鉬礦帶平面分布呈北西向延長的不規則紡綞狀，垂向上呈中間延深大，東西兩端薄而揚起的“船”形(圖14).礦體在走向上變化較大，相鄰孔見礦層數和厚度差別較大.礦體在垂向上大體呈串狀相疊，稀疏不均.橫向上主體延伸比較連續和平穩，部分有斷續和波狀起伏，整體具有成層組集中的特點.銅鉬礦體中主要為銅礦體，共生少量鉬礦體，分為同體共生或異體共生，以同體共生為主.斑巖型銅鉬礦多為同體共生，鉬礦體分布范圍較廣但礦體分散且規模均較小.矽卡巖型銅鐵礦產于中基性火山巖地層之中，為盲礦體.矽卡巖型銅鐵礦基本為同體共生礦體，但不完全重疊.其中，鐵礦整體賦存于火山巖地層之中，礦體呈沿傾向略長的馬鞍形透鏡體狀，平面投影范圍呈南北向略長的不規則橢圓狀.銅礦體為主要工業銅礦體，與磁鐵礦礦體賦存位置基本相同，但規模略小.礦體為透鏡體狀，局部有波狀彎曲和起伏，中部向上拱起，北段平緩近于水平，南段向南下部傾伏.通過塊段法，全區求得銅礦金屬量(333+334)118萬噸，全區銅礦伴生銀礦1404.44噸，鉬礦(333+334)6376噸，高品位磁鐵礦礦石量251萬噸.

圖14 拉伊克勒克礦床礦體三維形態MFe-磁鐵礦體；Mo-鉬礦體；Cu-銅礦體.Fig.14 Three-dimensional shape of the ore body in Layikeleke depositMFe-magnetite ore body；Mo-molybdenum ore body；Cu-copper ore body.

3 結論

(1)在總結新疆瓊河壩地區三類主要礦床均與巖體和構造有著密切關系的基礎上，通過1∶50000重磁資料的處理和挖掘，厘定了斷裂系統，識別了巖體分布，結合1∶5激電中梯的極化率和化探信息，預測了7處找礦靶區.對其中的拉伊克勒克靶區進行了查證，從預查到詳查，通過地球物理的創新組合，在地表沒有礦化線索的荒漠戈壁，新發現和評價了拉伊克勒克大型銅多金屬礦，包括大型斑巖銅鉬礦和小型矽卡巖型富鐵銅礦2個礦段，在瓊河壩地區荒漠戈壁覆蓋區實現找礦突破.

(2)合理的技術組合加正確的找礦思路是覆蓋區隱伏礦發現的關鍵，研究結果表明本文提出的荒漠覆蓋區預覆蓋預測與定位技術，在類似荒漠覆蓋景觀區具有一定的示范意義.

(3)拉伊克勒克大型銅多金屬礦的發現，說明在新疆、內蒙等戈壁荒漠覆蓋層下找礦潛力巨大，是國家戰略性礦產勘查的重要方向.

致謝新疆地礦局物化探大隊提供了瓊河壩地區1∶5重磁和激電數據，參加拉伊克勒克地球物理數據采集的有張振林、夏福成、趙建豪、劉有峰、郭建華、胡浩、丁峰、邵陸森等，參加鉆孔編錄、儲量計算等地質工作的有邵星月、楊超、張樹德、許亞寧、童英、吳健、李高峰、呂博、黃偉等，項目開展期間得到了新疆地勘基金得大力支持，中國地質調查局及西安地質調查中心和中國地質科學院礦產資源研究所領導專家多次到野外指導，拉伊克勒克銅多金屬礦的發現離不開上述個人和單位支持和幫助，三位匿名審稿人提出了建設性意見，在此一并表示衷心感謝！