閩西南大排多金屬礦區(qū)地球物理響應(yīng)特征與找礦啟示

孟小紅， 王俊*， 劉國峰， 陳召曦， 倪建輝， 張達(dá)， 譚捍東， 李淑玲， 方圓， 鄭師警

1 中國地質(zhì)大學(xué)(北京)地球物理與信息技術(shù)學(xué)院， 北京 100083 2 福建省第八地質(zhì)大隊， 福建龍巖 364000

0 引言

大排多金屬礦區(qū)位于武夷山成礦帶，其所處的閩西南地區(qū)是特提斯東西向構(gòu)造與環(huán)太平洋北東向構(gòu)造時空演化交替的典型地區(qū)，也是中國東南部重要的成礦聚集區(qū).在長期的地質(zhì)構(gòu)造演化過程中，研究區(qū)的中上地殼受到多期構(gòu)造事件改造，形成了不同時期的區(qū)域性隆起及坳陷，復(fù)雜多變的構(gòu)造形變格局以及不同時期頻繁的火山-侵入活動，使區(qū)內(nèi)蘊(yùn)藏著豐富的Fe、Cu、Pb、Zn、Au、Sn等多金屬礦床(點(diǎn))，具有良好的找礦前景.

大排多金屬礦區(qū)的研究歷史相對較長，礦產(chǎn)調(diào)查和科研工作經(jīng)歷了由淺到深的長期積累過程.20世紀(jì)70年代，大排多金屬礦區(qū)及其外圍就已經(jīng)進(jìn)行過了多次礦產(chǎn)勘查工作，當(dāng)時的研究工作認(rèn)為大排礦床類型為矽卡巖型，礦體規(guī)模較小、儲量少，埋藏有一定深度不易開采.自1999年實(shí)施國土資源大調(diào)查和2010年馬坑-湯泉鐵多金屬成礦區(qū)被列為國家首批整裝礦產(chǎn)勘查區(qū)以來，地質(zhì)礦產(chǎn)部門在大排礦區(qū)投入了多個礦產(chǎn)資源評價項目，開展了成礦地質(zhì)背景、成礦地質(zhì)條件及成礦規(guī)律等研究.比如，在成礦時代方面，趙希林等(2016)采集大排礦區(qū)的花崗閃長巖和花崗閃長斑巖樣品，分別進(jìn)行鋯石U-Pb測年和輝鉬礦Re-Os同位素測試，得到Re-Os等值線年齡為135±4.1 Ma，鋯石U-Pb年齡是127±1.8 Ma.袁遠(yuǎn)等(2013)測得大排似斑狀花崗巖鋯石U-Pb年齡為132.35±0.83 Ma及輝鉬礦Re-Os年齡為133±2 Ma，與趙希林等測得數(shù)據(jù)相近，由此推斷燕山期是本研究區(qū)的一個重要的巖體侵位和成礦時期.在控礦構(gòu)造方面，許乃政等(2008)認(rèn)為逆沖推覆構(gòu)造是大排礦區(qū)的主要控礦構(gòu)造，礦體主要賦存在由逆沖推覆構(gòu)造形成的破碎帶內(nèi)部.倪建輝(2012)通過對大排礦區(qū)詳盡的地質(zhì)調(diào)查，認(rèn)為經(jīng)畬組-棲霞組地層和礦區(qū)發(fā)育的淺層次級逆沖推覆構(gòu)造是主要的賦礦層位和控礦構(gòu)造；在礦床成因方面，張振杰等(2014)提出閩西南成礦與早白堊中酸性巖漿相關(guān)的礦床主要有鐵、鉛、鋅、銅、鉬等多金屬，大多屬于層控矽卡巖型礦床.倪建輝(2012)認(rèn)為大排鐵鉛鋅多金屬礦在成因上與區(qū)域上的馬坑鐵礦存在相似性，他們同處在閩西南鐵銅多金屬成礦帶上，屬于海相火山沉積-改造礦床.毛建仁等(2014)以鉆探揭露為主要手段，對礦區(qū)進(jìn)行了補(bǔ)充詳查地質(zhì)工作，投入鉆探工作量16267.00 m/38個，探獲鉛鋅資源量從原來的中型擴(kuò)大為大型，伴生的鐵礦、銅礦、鉬礦、銀礦、錳礦等資源量達(dá)到中型，認(rèn)為礦床成因為層控-矽卡巖型→層控-熱液疊加改造型.

上述研究工作一方面使得對于大排多金屬礦區(qū)的成礦地質(zhì)背景、成礦地質(zhì)條件和礦床成因有了深入的了解，同時也建立了該區(qū)域新的找礦模型(倪建輝，2012；杜一波，2018；鄭明泉，2019)和“賦礦層位—控礦構(gòu)造—燕山期巖漿活動”三位一體的找礦思路：即研究區(qū)內(nèi)上石炭船山組至下二疊棲霞組的一套碳酸鹽建造是鉛鋅礦床重要的賦礦層位和礦源層(具多層性)，經(jīng)畬組碎屑巖與碳酸鹽建造是鐵礦化的主要層位和礦源層；鉛鋅多金屬礦體受一系列規(guī)模不等、傾向南東東的疊瓦狀淺層逆沖推覆斷裂和層間裂隙(破碎帶)控制，礦體產(chǎn)狀與地層產(chǎn)狀或斷裂產(chǎn)狀基本一致；燕山晚期二長花崗斑巖從深部帶來了部分成礦物質(zhì)，同時對已有的礦體起到熱液交代和改造作用.這一新的找礦模型和三位一體的找礦思路，將大排礦區(qū)外圍和深部礦產(chǎn)預(yù)測工作與推覆構(gòu)造界面、隱伏巖體分布和賦礦巖層的揭示緊密地聯(lián)系起來.前人研究表明，地球物理方法在識別斷裂構(gòu)造、發(fā)現(xiàn)隱伏巖體和確定標(biāo)志地層等方面具有獨(dú)特優(yōu)勢(Wang and Meng, 2019).但由于大排礦區(qū)地表地形復(fù)雜、植被覆蓋嚴(yán)重，以往很少進(jìn)行高精度大比例尺的綜合地球物理調(diào)查，特別是在本次研究之前，尚未開展過反射地震工作，因此，對隱伏的控礦構(gòu)造、賦礦層位和火成巖體分布的認(rèn)識遠(yuǎn)不能滿足深部找礦的需求.本次研究通過實(shí)施剖面反射地震、可控源音頻大地電磁和大比例尺航磁測量，在成礦理論和找礦模型的指導(dǎo)下采用反射地震真地表偏移技術(shù)、可控源音頻大地電磁二維帶地形反演技術(shù)和航磁數(shù)據(jù)精細(xì)處理與三維反演技術(shù)，提取識別了淺層逆沖推覆控礦構(gòu)造和賦礦層位的反射地震特征，與成礦關(guān)系密切的燕山期火成巖體的地震、電磁特征，以及控礦構(gòu)造和隱伏巖體平面分布的航磁異常特征.在地表地質(zhì)和已知鉆孔巖心及物性的約束下，構(gòu)建了剖面綜合地球物理地質(zhì)解釋模型，包括淺層推覆構(gòu)造和控礦地層分布，推斷解釋大斷裂6條、次級斷裂16條和3處較大的隱伏巖體分布，為深入認(rèn)識研究區(qū)的礦床成因和分布規(guī)律，指導(dǎo)在外圍和深部找礦提供了地球物理有效響應(yīng)特征.

1 地質(zhì)背景

大排多金屬礦區(qū)位于福建省西南部，在大地構(gòu)造位置上屬于華南地塊內(nèi)的華夏地塊(圖1).政和—大埔斷裂將福建省分為兩部分，其中西部為永安—梅州凹陷帶，東部為閩東火山斷坳帶，本次研究區(qū)則位于永安—梅州凹陷帶的西南部分，是東西向特提斯構(gòu)造與東北向太平洋大陸邊緣構(gòu)造重疊的典型區(qū)域(Yang et al., 2017).

圖1 研究區(qū)構(gòu)造地質(zhì)圖(修改自李培等，2019)Fig.1 Simplified structural geological map of the study area (modified from Li et al., 2019)

研究區(qū)目前的構(gòu)造格局是不同時期、不同構(gòu)造作用疊加改造后的結(jié)果.根據(jù)前人的總結(jié)基本可概括如下：晉寧運(yùn)動期，晚太古代的古老陸塊與早元古代麻源群發(fā)生強(qiáng)烈變形、變質(zhì)作用，形成了較穩(wěn)定的結(jié)晶基底；加里東運(yùn)動期，該區(qū)域的造山作用使其結(jié)晶基底發(fā)生了變質(zhì)及褶皺變化，伴隨有大量殼源型花崗巖的侵入，形成了該區(qū)域的褶皺基底；海西—印支運(yùn)動期是區(qū)內(nèi)沉積蓋層形成的主要階段，在該時期，地層以升降作用為主，形成了張性的較為封閉的裂谷盆地.期間，從晚泥盆紀(jì)至早三疊紀(jì)接受到了來自其北部和西部的陸源碎屑，形成了巨厚的海相沉積建造，伴隨有基性海相火山活動，形成了研究區(qū)晚古生代拗陷帶.目前，在石炭紀(jì)地層中發(fā)現(xiàn)了大量Pb、Zn、Mn礦床及礦點(diǎn)，Cu、Pb、Zn、Ag、Sb礦床大多產(chǎn)于石炭紀(jì)和早侏羅紀(jì)地層中，在二疊紀(jì)地層中局部也有礦床分布.早三疊紀(jì)晚期，印支造山運(yùn)動導(dǎo)致洋盆開始閉合，區(qū)內(nèi)由拉張的裂谷環(huán)境轉(zhuǎn)化為以擠壓為主的造山環(huán)境，以特提斯構(gòu)造域為主體的構(gòu)造格局漸變?yōu)橐原h(huán)太平洋活動大陸邊緣為主的構(gòu)造格局；燕山運(yùn)動期，研究區(qū)全面轉(zhuǎn)入中國東部大陸邊緣活動階段，以斷陷盆地發(fā)育和中酸性火山-侵入巖漿活動為主要特征，奠定了現(xiàn)今的構(gòu)造格局，該期構(gòu)造活動對先期形成的礦體起到熱液疊加和改造作用，使礦床明顯受斷裂和裂隙控制(許乃政等，2008)，圖4為研究區(qū)地質(zhì)圖.

根據(jù)地層巖性、巖相、變質(zhì)特征及接觸關(guān)系，研究區(qū)的地層可劃分為前泥盆紀(jì)基底巖系、晚古生代-中三疊紀(jì)以海相沉積為主的蓋層巖系以及中新生代陸相碎屑-火山巖系(張達(dá)，1999；林東燕，2011)，即從泥盆系到第四系地層均有發(fā)育.出露的地層主要有石炭系船山組，二疊系棲霞組、文筆山組、童子巖組和翠屏山組地層，較老的地層則在鉆孔中揭露.其中，船山組-棲霞組地層是經(jīng)歷了海相-淺海相沉積形成的碳酸鹽地層，巖性為灰?guī)r和大理巖等，頂部和上覆的文筆山組、童子巖組地層以斷層或整合接觸；文筆山組地層是在瀉湖相環(huán)境下沉積形成的，巖性為泥質(zhì)粉砂巖、細(xì)砂巖等，與下伏船山組-棲霞組碳酸鹽巖地層的斷層接觸面常形成較寬的破碎帶，并發(fā)生熱液蝕變和鐵、錳、鉛鋅等礦化.童子巖組地層巖性主要是泥巖、粉砂巖和砂巖等，巖性較軟，與下伏文筆山組地層斷層或整合接觸，在斷層面和層間滑動面上可看到硅化、黃鐵礦化、絹云母化等蝕變現(xiàn)象；翠屏山組地層巖性為粉砂巖、細(xì)砂巖和泥巖夾煤層的海陸交互相碎屑巖，與下伏的童子巖組地層之間常發(fā)育斷層，斷層面可見硅質(zhì)和鐵錳礦化.圖2是研究區(qū)的主要出露地層的巖性柱狀圖.

圖2 大排礦區(qū)主要地層巖性柱狀圖(修改自杜一波，2018)Fig.2 Lithology histogram of the main strata in Dapai deposit (modified from Du, 2018)

研究區(qū)位于中國東南大陸中生代火山—侵入巖帶西部，巖漿活動具多旋回性，十分頻繁，每個構(gòu)造發(fā)展階段都有相應(yīng)的巖漿巖組合出現(xiàn)，在大排礦區(qū)侵入巖的巖性主要以酸性-中酸性為主，在地表出露較少，主要是隱伏巖體，經(jīng)過鉆孔巖心取樣和觀察(圖3)，巖性有二長花崗巖、花崗閃長巖斑巖、花崗斑巖、輝綠巖等以及少量的基性巖，其中，似斑狀二長花崗巖和花崗閃長巖在侵位空間和侵位時代上與大排礦區(qū)成礦關(guān)系密切.同時，研究區(qū)位于閩西南逆沖推覆構(gòu)造體系的南平—龍巖逆沖推覆構(gòu)造帶上，屬于王莊—樟坑逆沖推覆構(gòu)造帶的南部.燕山早期形成的由一系列北北東向或近南北向逆沖斷層組成的淺層次級推覆構(gòu)造，為巖漿熱液提供了通道，其形成的構(gòu)造斷裂面是礦質(zhì)重新運(yùn)移、沉淀、富集的重要部位，鉛鋅多金屬礦多賦存于該系列斷層及其平行的裂隙中.本次研究重點(diǎn)針對控礦構(gòu)造、賦礦層位和與礦產(chǎn)關(guān)系密切的隱伏巖體開展了地球物理特征提取和推斷解釋.

圖3 大排礦區(qū)主要侵入巖巖心(修改自杜一波，2018)Fig.3 Main invade rock cores of Dapai deposit (modified from Du, 2018)

2 地球物理數(shù)據(jù)采集與處理

為了能夠更有針對性地開展地球物理工作，我們對研究區(qū)內(nèi)的主要巖石類型的物性資料進(jìn)行了系統(tǒng)的收集和整理，同時對采集的巖心標(biāo)本進(jìn)行了密度、P波速度測量(表1)，總體來說研究區(qū)內(nèi)與物探關(guān)系密切的大類巖性中各類砂巖、粉砂巖、灰?guī)r的磁化率呈弱磁性或無磁性，平均密度為2.6 g·cm-3，平均P波速度為5894 m·s-1，電阻率為102～104Ωm；泥巖的磁化率也呈弱磁性或無磁性，平均密度較灰?guī)r略低，為2.57 g·cm-3，平均P波速度為5923 m·s-1，泥巖隨著泥質(zhì)的增加，電阻率也不斷降低，平均為10～102Ωm；花崗巖類具有相對較高的磁化率，平均密度與泥巖接近，為2.57 g·cm-3，P波速度較低，為5490 m·s-1，電阻率較高，為6×(102～105) Ωm；研究區(qū)內(nèi)，與其他巖類相比，二長花崗巖、花崗閃長斑巖、輝綠巖以及矽卡巖具有較高的磁化率(80～2200(4π×10-6SI))和較高的電阻率(6×(102～105) Ωm)，與灰?guī)r伴生的大理巖或大理巖化的巖石具有較低的P波速度(3947 m·s-1).這些大類巖性之間的物性差異，為開展物探工作奠定了地質(zhì)地球物理物性基礎(chǔ).

表1 巖心物性參數(shù)表Table 1 Physical properties of drilling core

本次研究實(shí)施了約100 km2的1∶10000直升機(jī)航磁測量，約13.5 km長的剖面反射地震測量和可控源音頻大地電磁測量，對實(shí)測數(shù)據(jù)進(jìn)行了精細(xì)處理和真地表偏移成像和反演.圖4是研究區(qū)地質(zhì)圖和航磁測區(qū)范圍以及剖面地震、可控源音頻大地電磁測線位置分布.

圖4 研究區(qū)地質(zhì)圖和測區(qū)測線位置圖(地震和CSAMT剖面)Fig.4 Geological map with survey line location of Dapai deposit (seismic and CSAMT profiles)

2.1 航磁數(shù)據(jù)采集與處理

航磁數(shù)據(jù)采集使用的是AS-350B3型直升機(jī)搭載航空磁測系統(tǒng)，配備有CS-3型高精度銫光泵磁力儀，PEI MAG Comp航磁補(bǔ)償系統(tǒng)、Hemisphere R320導(dǎo)航定位系統(tǒng)和TRA-3000雷達(dá)高度計等，地面日變觀測站使用GB-4A氦光泵磁力儀，空中數(shù)據(jù)收錄系統(tǒng)采用航空操作平臺IMPAC，實(shí)時顯示測量曲線、監(jiān)測各參量的工作質(zhì)量.采用沿地形起伏飛行方式，飛行方向為EW向，平均飛高為112 m，主測線間距為100 m，測量區(qū)域基本覆蓋了礦床的整個范圍.同時測量了總長40.2 km的重復(fù)線進(jìn)行精度控制.經(jīng)各項改正及調(diào)平后全區(qū)航磁異?？偩葹椤?.02 nT.利用Geosoft公司的Oasis montaj軟件對航測數(shù)據(jù)進(jìn)行了網(wǎng)格化處理，網(wǎng)度為100 m×100 m，圖5a為實(shí)測總場磁異常圖，由于研究區(qū)地表地形復(fù)雜、植被茂密覆蓋率高，地面面積性磁測工作幾乎難以進(jìn)行，本次采集的航磁數(shù)據(jù)消除了地表地形和隨機(jī)干擾的影響，是研究區(qū)首次獲得的高精度高質(zhì)量大比例尺的磁測數(shù)據(jù).

由于磁性地質(zhì)體受地磁場斜磁化的影響，所產(chǎn)生的磁異常往往不對稱且極值位置與異常體的中心位置相偏移，給解釋工作帶來困難，因此，磁異常數(shù)據(jù)處理首先需要把斜磁化的磁異常轉(zhuǎn)化為垂直磁化的磁異常，以消除斜磁化對磁異常形態(tài)的影響.本次研究采用非對稱因子化極計算方法(Guo et al., 2013)，選擇測區(qū)中心位置的磁傾角37°25′和磁偏角-3°35′，對實(shí)測航磁數(shù)據(jù)進(jìn)行了化極處理(圖5b).

實(shí)測磁異常是下半空間所有磁性地質(zhì)體產(chǎn)生的磁異常的疊加，為了將區(qū)域磁異常和局部剩余磁異常有效分離，我們采用自動趨勢面分析方法對化極后的磁異常進(jìn)行了異常分離處理(Obasi et al., 2016).利用多項式擬合區(qū)域背景，通過與研究區(qū)內(nèi)大型地質(zhì)結(jié)構(gòu)的形態(tài)對比，確定多項式近似的階次為二次，擬合區(qū)域異常(圖5c)，然后將實(shí)測異常減去區(qū)域異常獲得了局部剩余異常(圖5d).定性分析認(rèn)為，區(qū)域磁異常主要由變質(zhì)巖基底形態(tài)引起，總體呈“南高北低”趨勢，形態(tài)上類似于“傾斜臺階”異常，反應(yīng)了研究區(qū)基底構(gòu)造形態(tài)表現(xiàn)為特提斯東西向構(gòu)造特征，在區(qū)域異常西部疊加了一個形狀類似于“柱體”的異常，推測深部存在有磁性火山機(jī)構(gòu)；剩余磁異常整體上表現(xiàn)為環(huán)太平洋北東向構(gòu)造特征，測區(qū)西部較大的正異常，推測為深部磁性火山機(jī)構(gòu)在淺部的延續(xù)，測區(qū)中部和東部可明顯看出有北西向展布的正異常區(qū)和北西、北東向交叉展布的正異常區(qū)，推測為北西向淺層推覆構(gòu)造和北東向火山機(jī)構(gòu)造成的巖漿通道所致.剩余磁異常主要由基底以上磁性地質(zhì)體和不同深度的線性構(gòu)造(包括斷裂)所引起，因此，下面的異常特征提取和三維反演主要針對剩余磁異常進(jìn)行.

圖5 (a) 航磁總場磁異常； (b) 航磁化極異常； (c) 區(qū)域磁異常； (d) 剩余磁異常Fig.5 (a) Aeromagnetic total field magnetic anomaly; (b) Reduced-to-pole magnetic anomaly; (c) Regional magnetic anomaly; (d) Residual magnetic anomaly

2.2 反射地震數(shù)據(jù)采集與處理

本次研究采用的是Sercel 428XLTM數(shù)字地震儀，所采集的地震數(shù)據(jù)是在閩西南大排多金屬礦區(qū)首次開展的反射地震探測.因工區(qū)常年植被覆蓋，地形陡峭，施工難度極大，且無經(jīng)驗可循.在現(xiàn)場反復(fù)試驗后，確定的采集系統(tǒng)和采集參數(shù)為：2D直線觀測，中間激發(fā)兩邊接收，64次覆蓋，單井炸藥震源，平均藥量8 kg，井深14 m，采樣間隔0.5 ms，采樣長度5 s，炮點(diǎn)距60 m，檢波點(diǎn)距20 m.

由于地表的非均勻性和地表高差變化較大等自然因素，所采集的地震數(shù)據(jù)信噪比較低，幾乎在野外單炮記錄中看不到有效的地震信號，采用常規(guī)處理流程處理后，疊前時間偏移亦未獲得很好的成像效果，剖面上很難識別有效的同相軸，難以進(jìn)行構(gòu)造和巖性等地質(zhì)解釋，因此我們投入了大量精力，針對復(fù)雜地表情況，開展了改進(jìn)成像效果的方法技術(shù)研究.

通過分析前人研究和前期工作積累，我們認(rèn)識到能否真實(shí)地保留近地表信息是影響復(fù)雜地表條件下金屬礦地震勘探成像效果的重要因素之一.地震常規(guī)處理針對的是目的層較深的情況，經(jīng)常將近地表信息視為干擾，通過建立浮動基準(zhǔn)面靜校正后，一般更注意深部的成像結(jié)果，不太關(guān)注淺層500 m甚至更深的近地表信息，而這些淺層信息在隱伏固體礦產(chǎn)勘探中經(jīng)常是要向下追蹤的控礦構(gòu)造，有時甚至是直接要尋找的礦體的有效信息.因此，金屬礦地震資料處理除了要獲取深層地質(zhì)信息外，還需要真實(shí)地保留近地表信息，才能夠保證目標(biāo)深度成像結(jié)果的清晰性和與地表地質(zhì)情況的連續(xù)性.

另外，地震資料處理分為基于時間偏移的處理流程和基于深度偏移的處理流程，前者通過靜校正和時間偏移實(shí)現(xiàn)聚焦，后者通過深度域速度模型和深度偏移實(shí)現(xiàn)地下位置的正確歸位，前者的優(yōu)點(diǎn)是便于聚焦，但構(gòu)造成像位置不準(zhǔn)，后者的優(yōu)點(diǎn)是構(gòu)造成像位置準(zhǔn)確，但在沒有地震同相軸的原始資料上不能直接使用.因此，深度域速度模型是實(shí)現(xiàn)深度偏移，進(jìn)而實(shí)現(xiàn)精確構(gòu)造成像的關(guān)鍵.深度域速度模型的建立，首先要保證在時間域成像已經(jīng)聚焦的同相軸在深度域成像中仍然可以聚焦，即要求深度域速度模型和時間域速度模型等效，同時又要求深度域速度接近真實(shí)，只有在精確地構(gòu)建了全深度域速度模型的基礎(chǔ)上，進(jìn)行從真地表開始的深度偏移成像，才能保證在復(fù)雜地表條件下深度偏移成像結(jié)果聚焦和位置準(zhǔn)確.基于以上認(rèn)識，我們開展了保留近地表信息、構(gòu)建全深度域速度模型、進(jìn)行從真地表開始的深度偏移成像方法研究與應(yīng)用.

首先，我們開展了從真地表開始的疊前深度偏移成像研究(Liu et al., 2018)，這需要解決兩個問題：一是數(shù)據(jù)消除靜校正量后回到地表；二是結(jié)合近地表速度建模建立全深度域的速度模型.我們將進(jìn)行完反褶積后的地震數(shù)據(jù)重新恢復(fù)到地表面上，得到了包含了整個處理流程在地表上的地震數(shù)據(jù).然后是速度建模方面，主要也包含兩部分：一是近地表速度建模，采用初至波波形反演獲得；二是深部速度，采用常規(guī)處理手段獲得.然后將二者拼合成統(tǒng)一的速度模型，即全深度域速度模型，進(jìn)一步應(yīng)用全波形反演修改完善全深度域速度模型，再進(jìn)行真地表疊前深度偏移成像.

圖6是基于全深度域速度模型和全波形反演的真地表深度偏移成像結(jié)果，可以看出，剖面上界面反射同相軸、層位產(chǎn)狀變化、斷面波形態(tài)等非常清楚，清晰地反應(yīng)了3000 m深度上的地層結(jié)構(gòu)和斷裂構(gòu)造特征.

圖6 反射地震真地表深度偏移成像Fig.6 Depth migration imaging of reflection seismic data from real surface

2.3 可控源大地電磁數(shù)據(jù)采集與處理

為了克服大地電磁法中天然場源的隨機(jī)性及信號微弱的缺點(diǎn)，加拿大學(xué)者提出了利用CSAMT方法，該方法使用接地導(dǎo)線或不接地回線為場源，在波區(qū)測量相互正交的電場、磁場切向分量，并根據(jù)測量結(jié)果計算出卡尼亞視電阻率及阻抗相位，目前已應(yīng)用于普查、金屬礦產(chǎn)、石油勘探等方面.

圖7 (a) CSAMT數(shù)據(jù)采集觀測系統(tǒng)； (b) CSAMT源與測試位置圖Fig.7 (a) Data acquisition and observation system of CSAMT; (b) The source of CSAMT and its location

本次研究CSAMT數(shù)據(jù)采集設(shè)備采用GDP-32，采集系統(tǒng)如圖7所示，接收機(jī)位于排列中央，兩側(cè)各布置三個不極化電極，與接收機(jī)處不極化電極形成六道電道，不極化電極間點(diǎn)距80 m，排列總長480 m.采用兩個有限長導(dǎo)線發(fā)射源，赤道裝置，最大電流分別為15 A和14 A，發(fā)射源和采集線間距離6 km，每個采集點(diǎn)包含6個接收電極和一個與之垂直的接收磁極，采集頻率為8192～0.125 Hz，17個頻點(diǎn)采集數(shù)據(jù)，手動換頻，發(fā)射和接收采用石英鐘XMT對時，保證時間上嚴(yán)格一致.

在CSAMT數(shù)據(jù)處理中，進(jìn)行了干擾數(shù)據(jù)剔除、靜校正和過渡區(qū)校正等常規(guī)處理.反演計算中用有限差分法模擬考慮地形的CSAMT響應(yīng)，用偽δ函數(shù)來代替場源直接計算總場，采用自主研發(fā)的帶地形的反演方法(Lin et al., 2018; Zhang et al., 2016)進(jìn)行了二維反演.帶地形的CSAMT反演與地震數(shù)據(jù)的真地表偏移成像具有相同的成像空間，為后續(xù)的地質(zhì)解釋奠定了基礎(chǔ)，圖8是剖面的二維帶地形反演結(jié)果，可以看出在接近地表(深度500 m以內(nèi))處以小范圍的高阻體和低阻體交錯分布為主，隨著深度的增加，視電阻率逐漸增加，呈現(xiàn)出以高阻為主的態(tài)勢.剖面上從西向東視電阻率由相對低向相對高變化，分布有三個規(guī)模較大、深度較深的高阻異常體，推測這三個高阻體可能與巖漿活動有關(guān).剖面東部有一個規(guī)模較大形狀倒置的高阻異常體，推測可能與推覆體構(gòu)造有關(guān).

圖8 CSAMT數(shù)據(jù)反演的電阻率剖面Fig.8 Inverted resistivity profile from CSAMT data

3 地球物理有效特征提取與綜合解釋

在獲得了高質(zhì)量的一手實(shí)測數(shù)據(jù)和處理結(jié)果后，我們首先利用航磁異常數(shù)據(jù)在平面上的分辨能力，對局部剩余磁異常進(jìn)行了線性特征提取和三維反演，推斷解釋了22條平面斷裂分布和3處隱伏磁性巖體分布.

在對剖面成像結(jié)果進(jìn)行地質(zhì)解釋時，我們首先對地震剖面附近的兩口深井進(jìn)行了巖心采樣和物性測量，建立了井旁水平層狀介質(zhì)模型，通過與模型正演對比，對地震真地表深度偏移成像剖面的有效反射層位和地層進(jìn)行了標(biāo)定和識別，同時利用地震剖面上的斷面波形態(tài)以及層位產(chǎn)狀變化和視電阻率剖面上的高阻與低阻交界位置進(jìn)行了斷裂構(gòu)造交互識別.進(jìn)一步綜合地震剖面上的無反射特征、視電阻率剖面上的高阻異常、航磁三維反演結(jié)果給出的高磁化率異常和鉆井資料，進(jìn)行了剖面隱伏巖體識別.最終構(gòu)建了包括巖性、斷層和巖體的剖面地質(zhì)地球物理綜合解釋模型.

3.1 基于航磁線性構(gòu)造特征和三維物性反演的平面構(gòu)造體系劃分和隱伏巖體分布推斷

磁異常線性特征提取與增強(qiáng)是突出與斷裂構(gòu)造或地質(zhì)體邊界相對應(yīng)的構(gòu)造信息的一種重要方法.Wang和Meng(2019)利用傾斜角方法對大排礦區(qū)的磁異常進(jìn)行了斷裂推斷，本次研究在其基礎(chǔ)上將多種線性構(gòu)造增強(qiáng)算法(總水平導(dǎo)數(shù)、θ圖、歸一化標(biāo)準(zhǔn)差)應(yīng)用于剩余磁異常，綜合分析開展了基于航磁數(shù)據(jù)的隱伏斷裂推斷.共推斷大斷裂6條、次級斷裂16條，其中大斷裂與已有研究結(jié)果基本一致，次級斷裂則存在著一定的差別.

圖9是總水平導(dǎo)數(shù)法處理后的剩余磁異常線性特征與推斷斷裂的疊合圖，其中，標(biāo)號F的斷裂(F1—F16)為研究區(qū)已知斷裂，標(biāo)號D的斷裂(D1—D6)和標(biāo)號S的斷裂(S1—S16)為本次研究推斷的大斷裂和次級斷裂.

圖9 研究區(qū)推斷斷裂和磁異常線性特征疊合Fig.9 The overlap of fault inference with magnetic linear features of the study area

將已知斷裂和推斷斷裂分別疊合在區(qū)域磁異常和剩余磁異常圖上(圖10)，從區(qū)域磁異常疊合圖10a中可看出，已知斷裂中的F1、F2、F4、F10，推斷次級斷裂中的S1、S2和S3 等北東向斷裂與區(qū)域磁異常北東向分布趨勢具有較好的對應(yīng)性，反映出北東向斷裂與研究區(qū)磁性基底構(gòu)造特征具有很好的繼承性.從剩余磁異常疊合圖10b中可看出，北西向主要斷裂D1、D2、D5，南北向主斷裂D3、D4等基本位于磁性巖體的邊界，北西向和南北向次級斷裂中的S6、S7、S8、S9、S16等也都位于磁性巖體的邊界位置或和磁性巖體的走向一致(后期的反演結(jié)果也可看出)，說明由磁性特征推斷的北西向和南北向斷裂與研究區(qū)蓋層中的磁性巖漿活動密切相關(guān)，很好地反映了磁性巖體的巖漿侵入通道和賦存狀態(tài).

圖10 研究區(qū)斷裂推斷圖(a)在區(qū)域磁異常上的疊加圖；(b)在剩余磁異常上的疊加圖Fig.10 The deduced faults in the study area overlay image on (a) regional magnetic map; (b) residual magnetic map

圖11 研究區(qū)航磁3D反演結(jié)果(a) 3D反演結(jié)果透視圖，其中紅線、藍(lán)線是推斷斷裂，黑線是剖面AB位置; (b) 3D反演結(jié)果在剖面AB處的切片.Fig.11 3D inversion result of magnetic data in the study area(a) Perspective view of the inversion result. The red and blue lines are the inferred faults, and the black line is profile AB; (b) Slices along profile AB of the 3D inverse result.

綜合來看研究區(qū)內(nèi)已知和由磁異常線性特征推斷的隱伏斷裂主要有北東向、近南北向和北西向三組，其中中部為近南北走向斷裂，西部和東部為北西、北東向斷裂，斷裂的性質(zhì)和成因為燕山早期形成的逆沖斷層所組成的淺層次級推覆構(gòu)造，這些斷裂構(gòu)造反映了與成礦關(guān)系密切的巖漿侵入通道和巖體賦存狀態(tài)，根據(jù)礦床成因和找礦模式研究，其形成的構(gòu)造斷裂面是礦質(zhì)運(yùn)移、沉淀和富集的重要部位，對在外圍和深部找礦具有參考意義.

磁異常三維反演可生成具有磁化率對比的三維地質(zhì)體分布，可有效地用于確定潛在磁性場源體形狀和深度.本研究采用經(jīng)典的3D物性反演算法對研究區(qū)航磁數(shù)據(jù)進(jìn)行了反演成像.反演中將地下劃分為一組棱柱體單元，同時假設(shè)每個單元內(nèi)的磁化率恒定，初始值設(shè)置為磁化率為零的常數(shù)模型，反演過程中，加入了深度約束與物性約束.分析認(rèn)為，異常應(yīng)由地下3 km內(nèi)的場源引起的，所以反演深度范圍為0～3 km.整個測區(qū)被劃分了170×100×30=510000個棱柱體單元，每個棱柱體單元的維度是100 m×100 m×100 m.

圖11為研究區(qū)航磁剩余異常3D反演結(jié)果的透視圖和在AB剖面上的切片.從透視圖11a可以看出，研究區(qū)內(nèi)由西向東存在著三個較大的磁性巖體，磁性異常體及其邊界分別與化極磁異常的高值部分和推斷的斷裂有很好的吻合，進(jìn)一步說明了研究區(qū)磁異常特征很好地反映了磁性巖漿侵入通道和巖體的賦存狀態(tài).從剖面切片圖11b可以看出，在剖面所處位置以下分布著兩個磁性體，其中位于西側(cè)、在兩個鉆孔zk1001和zk2001附近的磁性體體積相對較大，頂部位于地下150 m左右，深部可達(dá)到2800 m，整體呈近似楔形，沿著近北西—南東方向展布，東側(cè)的磁性巖體相對較小較淺，深度范圍為地下150～1600 m，磁性體在上部150～700 m深度時是三個獨(dú)立部分，700 m以下三個獨(dú)立部分相連通構(gòu)成了一個整體.另外，在剖面以外研究區(qū)的東北部存在著一個較大的磁性巖體，從地下350 m一直延伸至2500 m左右，異常體北側(cè)部分沿著近北東向展布，南側(cè)部分沿著近南東向展布.

根據(jù)前人研究，研究區(qū)內(nèi)燕山晚期(127.0±1.8 Ma)侵入的巖漿巖與成礦關(guān)系密切，為成礦提供了含多金屬礦化熱液，對沉積介質(zhì)中的鐵鉛鋅多金屬起到活化、疊加、轉(zhuǎn)移和富集的作用.鉆井巖心顯示，侵入巖巖性為二長花崗巖、花崗閃長巖斑巖、花崗斑巖和輝綠巖等，呈灰色或深灰色，含暗色礦物成分較多，具有相對較高的磁性，是研究區(qū)磁異常的主要場源.因此磁異常三維反演結(jié)果給出了與成礦關(guān)系密切的巖漿巖巖體的空間賦存狀態(tài).

3.2 鉆孔巖心取樣與井旁地震正演模擬

鉆孔zk1001和zk2001位于地震和電法剖面附近(圖12)，其中，zk1001鉆孔距離剖面垂直距離約300 m，井深800 m，zk2001鉆孔距離測線垂直距離約600 m，井深1100 m，兩個鉆孔均在鉆到火成巖后停鉆，并保留有完整的巖心，按照巖心地層序列，我們在這兩口鉆井中進(jìn)行了巖心采樣，在實(shí)驗室內(nèi)采用超聲脈沖法測量了巖石樣品的縱速度，并通過測量標(biāo)本的質(zhì)量和體積，獲得到了巖心的密度，表1是實(shí)驗室獲得的主要巖心的物性參數(shù).

圖12 巖心的巖性及其合成地震記錄(a) zk1001; (b) zk2001.Fig.12 The lithology and the synthetic seismic trace of (a) zk1001, (b) zk2001

根據(jù)鉆井巖心資料以及實(shí)驗室內(nèi)測量的巖心物性參數(shù)，建立了對應(yīng)于兩個井位層位的水平層狀介質(zhì)模型(忽略了厚度小于10m的地層以及溶洞等形態(tài)構(gòu)造)，計算了界面處的反射系數(shù)，采用雷克子波與反射系數(shù)褶積獲得了物性結(jié)構(gòu)的合成地震記錄，如圖13連井地震剖面所示，其中右圖為左圖的放大顯示，合成地震記錄置于地震剖面中，可以看出，剖面的反射同相軸與合成地震記錄具有較好的對應(yīng)關(guān)系.在大理巖出露層位由于速度較低與圍巖差異顯著，形成較強(qiáng)的阻抗，在合成地震記錄上表現(xiàn)為較強(qiáng)的反射.

圖13 連井地震剖面，右圖為左圖的放大顯示Fig.13 Well seismic profiles, the right picture is the enlarged display of the left one

前人研究表明，研究區(qū)內(nèi)礦床的主要含礦巖系為晚石炭系經(jīng)畬組和中二疊系棲霞組的濱?！獪\海相碎屑巖及灰?guī)r.深孔揭露結(jié)果顯示，經(jīng)畬組灰?guī)r大部分變質(zhì)成白色、乳白色大理巖或淺色大理巖化灰?guī)r；棲霞組灰?guī)r大部分也已大理巖化或成為大理巖.

本次鉆井資料分析和巖心物性測量亦發(fā)現(xiàn)，鉛鋅礦和磁鐵礦出現(xiàn)的位置伴隨著大理巖巖層分布，而大理巖或大理巖化灰?guī)r與圍巖縱波速度有較大差異，反映在地震剖面上就是強(qiáng)的反射同相軸.將合成地震記錄并入到地震剖面所在位置，并進(jìn)行能量平衡后，可以看出，井旁合成地震記錄中的強(qiáng)反射與地震剖面中的強(qiáng)反射有很好的對應(yīng)，在磁鐵礦和鉛鋅礦出現(xiàn)的位置，因為大理巖的存在產(chǎn)生了明顯的反射，且振幅較強(qiáng)，雖然磁鐵礦和鉛鋅礦礦體的厚度較薄，頂?shù)追瓷洳荒軌蚍直?，但因為伴生的大理巖巖層的存在，所產(chǎn)生的強(qiáng)反射同相軸為間接識別賦礦層位提供了信息，可據(jù)此配合振幅級別等，推斷橫向主要控礦層位展布和相鄰縱向?qū)游环植?

3.3 基于地震真地表深度偏移成像和CSAMT帶地形二維反演的剖面地質(zhì)解釋與綜合建模

獲得“賦礦層位—控礦構(gòu)造—燕山期巖漿活動”信息是綜合地球物理精細(xì)解釋的目標(biāo)，閩西南大排礦區(qū)屬于推覆控礦構(gòu)造，成礦和控礦的主要地質(zhì)要素之一是推覆斷裂和推覆體內(nèi)部結(jié)構(gòu)，同時，燕山期的巖漿作用是礦區(qū)成礦物質(zhì)的主要來源，也是成礦賦存的重要標(biāo)志.在實(shí)施了剖面精細(xì)探測的基礎(chǔ)上，利用剖面反射地震成像、視電阻率反演和航磁三維反演結(jié)果剖面切片，結(jié)合鉆井巖心和物性資料，對剖面進(jìn)行了聯(lián)合解釋和建模.

地震剖面上因為斷面波和層位產(chǎn)狀的變化明顯，首先進(jìn)行推覆斷層的識別.為了進(jìn)一步判斷所識別斷層的準(zhǔn)確性，將從地震成像結(jié)果中識別出的斷層線拖曳到視電阻率剖面上，可以清晰地看出，大部分?jǐn)鄬泳€的位置位于高阻與低阻的轉(zhuǎn)換位置，以此為依據(jù)，在視電阻率剖面中識別出在地震剖面上特征不夠明顯的斷層，并將其添加到地震剖面上，依此交互識別交互驗證.最后以地表地質(zhì)圖上的已知斷裂向下追蹤對識別出的主要斷層進(jìn)行確定，以保證整體斷層識別的可信度，綜合勾畫出推覆斷層的框架.

由于巖體內(nèi)部的物質(zhì)分布相對均勻，沒有明顯的定向沉積，在地震剖面上呈現(xiàn)隨機(jī)反射的特點(diǎn)，波形雜亂，表現(xiàn)為無明顯反射的空白帶.另外，前述章節(jié)物性統(tǒng)計分析給出，研究區(qū)巖漿巖巖體具有相對較高的電阻率，在視電阻率剖面上表現(xiàn)為高阻，根據(jù)這兩個特征，在地震剖面和視電祖率剖面上交互進(jìn)行了火成巖巖體的識別.首先在地震剖面上圈出明顯無反射的區(qū)域，將其拖曳到視電阻率剖面上，通過視電阻率高阻異常邊界對地震識別的巖體進(jìn)行修改，將修改后的巖體分布拖回到地震剖面，進(jìn)行交互修改交互驗證.進(jìn)一步與磁異常三維反演切片獲得的視磁化率剖面進(jìn)行對比，剖面左側(cè)的巖體在磁化率剖面上顯示高磁性，進(jìn)一步證明了其為巖體的可能，剖面中部的巖體無明顯磁異常，考慮到研究區(qū)的地質(zhì)背景，地震和電法剖面上圈定的這個巖體有可能是無磁性或極弱磁性的酸性巖體，也可能是林地組砂巖巖體，有待于進(jìn)一步驗證，剖面右側(cè)的高磁異常，對應(yīng)了視電阻率剖面上的高阻異常，因此，亦推斷為小型隱伏巖體.

圖14是從地震剖面、電阻率剖面和磁法反演剖面上聯(lián)合識別的推覆斷層、隱伏巖體結(jié)果.結(jié)合圖12和圖13的巖心資料以及井震對比分析，再結(jié)合測區(qū)的地質(zhì)情況以及不同巖性的物理性質(zhì)，在圖15給出了地層、斷層和巖體的綜合解釋結(jié)果.從圖中可以看出，剖面由淺到深，地層依次為：第四系、早二疊紀(jì)童子山組、早二疊紀(jì)文筆山組、早二疊紀(jì)棲霞組、石炭紀(jì)經(jīng)畬組和石炭紀(jì)林地組，地層特征明顯穩(wěn)定并表現(xiàn)出被推覆構(gòu)造錯段向上逆沖的趨勢，其中石炭世經(jīng)畬組到早二疊紀(jì)棲霞組為本區(qū)的主要含礦層位，由于其伴生的大理巖層速度較低，產(chǎn)生了明顯的反射，且振幅較強(qiáng)，通過井資料標(biāo)定后為研究區(qū)控礦地層識別和地層解釋提供了依據(jù).剖面從西到東，發(fā)育了一系列淺層逆沖推覆構(gòu)造，其中，剖面右側(cè)的兩條斷裂向西傾斜，其他斷裂向東傾斜，根據(jù)研究區(qū)地質(zhì)概況，大排外圍北側(cè)的構(gòu)造格局為NNW 傾伏的復(fù)式倒轉(zhuǎn)背斜和南北向斷裂的組合，倒轉(zhuǎn)背斜的軸面傾向北東，走向北北西，沿著背斜兩翼發(fā)育有逆沖斷裂，同時也發(fā)育滑脫斷裂，斷裂以走向近南北為主.由于本研究剖面位于倒轉(zhuǎn)背斜南端末尾，因此，剖面上所看到的這兩種傾斜方向的斷裂反映的是北側(cè)背斜構(gòu)造特征的延續(xù).在背斜構(gòu)造核部即剖面上4 km位置附近——大排鐵鉛鋅多金屬礦區(qū)，推覆構(gòu)造、隱伏巖體和賦礦層位共存，斷裂為花崗巖體的侵入和成礦熱液的運(yùn)移提供了通道，賦礦層位為成礦成藏提供了儲存空間，另外，在12 km處較小的磁性巖體所對應(yīng)的電阻率高阻值異常出現(xiàn)上下倒置現(xiàn)象(位于地表低阻值異常體上方，即較老地層上覆于新地層之上)，推測可能為區(qū)域逆沖推覆體，根據(jù)推覆構(gòu)造控礦特征，推覆體結(jié)構(gòu)的重合是深部找礦的重要依據(jù)，因此推覆體下層是深部找礦需要關(guān)注的區(qū)域.

圖14 地震剖面、電阻率剖面、磁化率剖面綜合解釋斷層和巖體Fig.14 Joint interpretation of faults and magma invasion from seismic profile, resistivity profile and magnetic susceptibility profile

圖15 剖面綜合地質(zhì)解釋Fig.15 Comprehensive geological interpretation of profile AB

4 結(jié)論與認(rèn)識

綜合研究區(qū)成礦理論、找礦模型和綜合地球物理解釋結(jié)果，我們對閩西南大排多金屬礦區(qū)的地球物理響應(yīng)特征及其找礦意義有了如下認(rèn)識：

(1)研究區(qū)存在一系列規(guī)模不等的燕山早中期疊瓦狀逆沖推覆斷裂和層間裂隙，這些斷裂和斷層一方面為巖漿熱液提供了通道，同時，其形成的構(gòu)造斷裂面是礦質(zhì)重新運(yùn)移、沉淀、富集重要部位，是研究區(qū)的主要控礦構(gòu)造.

航磁異常線性構(gòu)造增強(qiáng)平面特征清晰地給出了與隱伏巖體邊界和侵位關(guān)系密切的斷裂構(gòu)造的平面分布信息，據(jù)此推斷出大斷裂6條，次級斷裂16條，其中大斷裂與已有研究結(jié)果基本一致，次級斷裂則多為首次推斷的隱伏斷裂；在反射地震成像剖面上，斷面波形態(tài)以及層位產(chǎn)狀變換等特征明顯，對推覆斷層結(jié)構(gòu)給出了清晰的勾畫，剖面左側(cè)有兩條北西傾向斷層，中部和右側(cè)有多條北東傾向斷層、構(gòu)造面傾角為30°～45°不等，地層被斷層錯斷明顯，表現(xiàn)出向上逆沖的趨勢；在視電祖率剖面上，高阻和低阻的邊界對地震斷層解釋給出了進(jìn)一步的補(bǔ)充和修正，進(jìn)一步通過地表已知斷裂向下追蹤延伸，驗證了剖面斷層解釋的正確性.

(2)研究區(qū)內(nèi)主要巖漿巖為燕山晚期二長花崗斑巖、花崗閃長巖、花崗斑巖和輝綠巖，鉆孔巖心顯示其顏色呈深黑色，說明其所含暗色礦物成分相對較多，物性測量結(jié)果顯示其具有較高的磁化率和電祖率，是研究區(qū)航磁剩余異常和視電祖率高阻異常的主要成因.

磁異常三維反演結(jié)果給出，在約100 km2的磁測區(qū)域內(nèi)有三個較大的隱伏磁性巖體分布；視電阻率剖面上在4 km、9 km、12 km位置處存在三個規(guī)模較大的高阻異常體，視電阻率數(shù)值達(dá)到了8000 Ωm；地震偏移成像剖面上巖體表現(xiàn)為波形雜亂無明顯反射的空白特征，通過地震、電法交互識別并與航磁反演結(jié)果疊合，勾畫出了剖面上的巖體分布特征.其中，平面左側(cè)巖體的北部邊界對應(yīng)于大排礦區(qū)的位置，剖面4 km處的巖體給出了其縱向分布特征，這些隱伏巖體為礦區(qū)成礦作用提供了多金屬礦化熱液，對沉積介質(zhì)中的鐵鉛鋅多金屬礦起到活化、疊加、轉(zhuǎn)移和富集的作用，同時也是開展深部和外圍找礦的重要指示標(biāo)志之一.

(3)研究區(qū)的主要含礦巖系為晚石炭系經(jīng)畬組和中二疊系棲霞組的濱海-淺海相碎屑巖及灰?guī)r.其中，經(jīng)畬組的巖性變化相對較大，棲霞組下部巖性為灰-深灰中厚層狀含燧石灰?guī)r、夾泥灰?guī)r等，其中部分的灰?guī)r已大理巖化或成大理巖；通過鉆井資料分析和巖心物性測量，鉛鋅礦出現(xiàn)的位置往往伴隨著大理巖層分布，而大理巖或大理巖化灰?guī)r表現(xiàn)出了與圍巖縱波速度上的較大差異，這個差異反映在地震剖面上就是強(qiáng)的反射同相軸，為間接判斷賦礦層位提供了依據(jù).

(4)大排礦區(qū)地表起伏大，橫向速度變化強(qiáng)烈，傳統(tǒng)的反射地震基于靜校正的處理以及時間偏移成像方法無法對近地表構(gòu)造進(jìn)行精確成像，本次研究采用真地表成像技術(shù)，通過層析和深層速度體的速度拼合構(gòu)建全深度域速度模型，實(shí)現(xiàn)了從真地表開始的深度偏移成像；同時，采用帶地形的CSAMT數(shù)據(jù)反演，從地表起算的電法成像結(jié)果與真地表地震成像剖面保持了相同的坐標(biāo)位置，為構(gòu)造綜合解釋奠定了很好的基礎(chǔ).

(5)綜合地表地質(zhì)、鉆孔巖心、真地表地震深度偏移成像剖面、帶地形的視電阻率剖面和航磁三維反演開展了剖面地質(zhì)模型建模解釋.地震剖面、視電阻率剖面和航磁反演剖面交互驗證互為約束，勾畫了構(gòu)造形態(tài)、巖體分布，地表地質(zhì)向下延伸、鉆孔巖心橫向標(biāo)定賦予了剖面斷層、巖體和地層的地質(zhì)內(nèi)涵，最終給出了從真地表到3000 m深度的綜合解釋地質(zhì)模型，清晰地給出了推覆構(gòu)造、控礦層位以及巖體分布，為三位一體的找礦思路提供了地球物理響應(yīng)特征.

致謝本研究受國家自然科學(xué)基金重點(diǎn)基金項目“推覆構(gòu)造控礦的地球物理響應(yīng)研究”(41530321)資助.