安徽銅陵舒家店銅礦三維音頻大地電磁探測(cè)研究

王永清，湯井田，肖曉，周聰，王顯瑩，任政勇

?

安徽銅陵舒家店銅礦三維音頻大地電磁探測(cè)研究

王永清1, 2，湯井田1, 2，肖曉1, 2，周聰1, 2，王顯瑩1, 2，任政勇1, 2

(1. 中南大學(xué)有色金屬成礦預(yù)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室，湖南長(zhǎng)沙，410083;2. 中南大學(xué)地球科學(xué)與信息物理學(xué)院，湖南長(zhǎng)沙，410083)

利用音頻大地電磁法，對(duì)安徽省銅陵礦集區(qū)舒家店斑巖型銅礦床進(jìn)行三維探測(cè)研究。對(duì)經(jīng)過(guò)預(yù)處理的實(shí)測(cè)AMT數(shù)據(jù)利用相位張量進(jìn)行維性分析；通過(guò)挑選適當(dāng)?shù)臏y(cè)點(diǎn)和頻域數(shù)據(jù)，設(shè)計(jì)合理的參數(shù)，利用WSINV3DMT程序進(jìn)行實(shí)測(cè)AMT數(shù)據(jù)3D反演，得到舒家店銅礦床的3D電性模型；利用模型響應(yīng)與實(shí)測(cè)AMT數(shù)據(jù)的擬合對(duì)比對(duì)反演結(jié)果進(jìn)行評(píng)價(jià)；利用不同深度、不同方向的空間切片討論3D反演模型的電性結(jié)構(gòu)特征。研究結(jié)果表明：該區(qū)深部2D/3D構(gòu)造特征明顯；所選反演模型整體精度較高，僅部分局部精度低；深部與代表巖體的高阻體相間共存的低阻體異常極可能是蝕變巖體或直接是礦體的反映。

斑巖型銅礦床；音頻大地電磁法；三維反演；銅陵礦集區(qū)

音頻大地電磁測(cè)深法(audio-frequency magnetotellurics，AMT)發(fā)展于大地電磁測(cè)深法(magnetotellurics，MT)，是以巖礦石的電性差異為基礎(chǔ)，通過(guò)在地表測(cè)量天然音頻電磁場(chǎng)的正交水平分量，并對(duì)其變化規(guī)律進(jìn)行研究，獲取地下電性分布的一種交變電磁勘探方法。該方法探測(cè)深度大，工作效率高，是礦產(chǎn)勘查工作的重要手段之一[1?3]。近年來(lái)，由于方法理論和計(jì)算機(jī)技術(shù)的發(fā)展，三維大地電磁建模與反演取得了長(zhǎng)足進(jìn)步，并逐漸得到應(yīng)用[4?5]。目前三維反演的方法主要非線(xiàn)性共軛梯度法[6?7]、快速松弛反演[8]等，特別是WSINV3DMT代碼[9?10]已得到廣泛應(yīng)用[11]。銅陵是我國(guó)六大有色金屬基地之一，是我國(guó)“矽卡巖型銅礦”的重要礦集區(qū)，所發(fā)現(xiàn)的矽卡巖型銅礦床在安徽省乃至國(guó)內(nèi)都具重要地位[12?14]。銅陵地區(qū)的地表礦產(chǎn)勘查工作程度較高，人們對(duì)區(qū)內(nèi)主要大多數(shù)地質(zhì)、物化探異常區(qū)進(jìn)行了驗(yàn)證或查證，對(duì)地表發(fā)現(xiàn)的礦床、礦點(diǎn)多進(jìn)行了鉆探工程控制；在重要礦床的詳查勘探中，控制深度在1 km以下。在區(qū)內(nèi)以尋找地表礦和淺部隱伏礦的時(shí)代已經(jīng)結(jié)束，20世紀(jì)80年代便進(jìn)入了以尋找深部隱伏礦為主的艱難時(shí)期，自20世紀(jì)90年代對(duì)區(qū)內(nèi)冬瓜山礦床、焦沖金礦詳查結(jié)束后，礦產(chǎn)勘查工作一直未出現(xiàn)明顯進(jìn)展。但進(jìn)入21世紀(jì)后，發(fā)現(xiàn)了姚家?guī)X大型金多金屬礦，新增了可觀的資源儲(chǔ)量，延長(zhǎng)了危機(jī)礦山服務(wù)年限，取得了良好的經(jīng)濟(jì)效益和社會(huì)效益[14?15]。近期，隨著深部找礦的深入開(kāi)展，區(qū)內(nèi)先后發(fā)現(xiàn)了舒家店外圍(深部)斑巖型銅礦、雞冠山銅金礦、桂山鉛鋅多金屬礦、胡村銅鉬礦、蛤蟆嶺金礦等新礦床，證明在銅陵礦集區(qū)深部找礦依然具有較大潛力[16?17]。本文利用音頻大地電磁法在舒家店銅礦床開(kāi)展電性結(jié)構(gòu)研究。在測(cè)區(qū)部署了60個(gè)測(cè)點(diǎn)，利用WSINV3DMT代碼[18?19]進(jìn)行三維反演；分析的侵入巖體的分布形態(tài)，并對(duì)成礦作用和成礦規(guī)律進(jìn)行 討論。

1 地質(zhì)背景

舒家店礦床位于舒家店短軸背斜軸部脊軸轉(zhuǎn)折部位，受巖漿巖與砂質(zhì)頁(yè)巖北西接觸帶巖體內(nèi)側(cè)構(gòu)造控制。舒家店背斜為一軸向北東不對(duì)稱(chēng)的短軸背斜，巖層產(chǎn)狀一般北西翼地層傾角較小，為20°~55°；南東翼地層傾角較大，為70°左右，部分地段有倒轉(zhuǎn)現(xiàn)象；北東外傾轉(zhuǎn)折端因受斷裂影響，產(chǎn)狀較零亂。出露的地層均為志留系上中統(tǒng)墳頭組砂質(zhì)頁(yè)巖、砂巖、粉砂巖及頁(yè)巖[20]。

舒家店背斜受來(lái)自北西與南東2個(gè)方向的力，在礦床附近受北西方向力的影響較大，礦床北東端地層即行向南轉(zhuǎn)折。在南北方向力作用下，礦區(qū)產(chǎn)生了北東向壓性斷裂，巖漿巖即沿此通道上升。巖體軸向約為NE40°，與區(qū)域性的高角度壓性逆斷裂方向基本一致。構(gòu)造破碎帶呈北東方向，也是礦區(qū)主要控礦構(gòu)造。礦區(qū)出露的主要侵入巖可分為3個(gè)巖體6種巖相帶，總體呈巖株?duì)町a(chǎn)于背斜北西翼近軸部的墳頭組砂頁(yè)巖中。結(jié)合區(qū)域巖漿巖活動(dòng)規(guī)律和一些特征微量元素對(duì)比，可判定本區(qū)火成巖屬于同源巖漿分異的產(chǎn)物。由同位素鉀氬法年齡測(cè)定結(jié)果和巖體之間穿插及突變關(guān)系，確認(rèn)3個(gè)巖體形成的先后順序是：最早為花崗閃長(zhǎng)巖體(由中心相花崗閃長(zhǎng)巖和邊緣相花崗閃長(zhǎng)斑巖組成)，年齡為214 Ma；其次為石英閃長(zhǎng)斑巖體(可穿插于花崗閃長(zhǎng)巖體中，亦可被輝石閃長(zhǎng)巖、閃長(zhǎng)巖穿插或包裹)，年齡為178 Ma，屬燕山期；最后為輝石閃長(zhǎng)巖體(由中心相輝石閃長(zhǎng)巖、過(guò)渡相閃長(zhǎng)巖及邊緣相閃長(zhǎng)玢巖組成)，年齡為134~148 Ma，處于燕山中期[19?21]。主要巖體屬于中淺成相，各巖體邊緣相出露較廣，自變質(zhì)作用較強(qiáng)，脈巖較發(fā)育，巖體內(nèi)見(jiàn)圍巖的頂垂體(花崗閃長(zhǎng)巖體的南段較多)，由此可知舒家店巖體屬于淺成剝蝕。

礦床位于繁昌火山盆地邊緣，舒家店短軸背斜軸部脊軸轉(zhuǎn)折部位，為受接觸帶內(nèi)側(cè)巖株體內(nèi)構(gòu)造控制的板狀與脈狀礦體。礦床埋深在150 m以下，剝蝕程度較淺。與區(qū)內(nèi)礦化密切相關(guān)的蝕變主要是鉀長(zhǎng)石化，其次是硅化、綠泥石化，蝕變和礦化程度也一致。黃銅礦與鉀長(zhǎng)石化關(guān)系較密切，鏡鐵礦主要在鈉長(zhǎng)石化、伊利石化等巖石中比較發(fā)育。礦石構(gòu)造以浸染狀為主，細(xì)脈和小細(xì)脈次之，粗脈較少。巖體與圍巖界限不明顯。金屬礦物組分較簡(jiǎn)單，主要是黃銅礦、黃鐵礦，礦石含銅品位較低而穩(wěn)定。礦床賦存的主要巖體為閃長(zhǎng)巖和輝石閃長(zhǎng)巖，無(wú)斑狀結(jié)構(gòu)，雖然閃長(zhǎng)玢巖、石英閃長(zhǎng)斑巖具有斑狀結(jié)構(gòu)，但主要礦體不賦存在此類(lèi)巖石中，礦床成因與輝石閃長(zhǎng)巖、閃長(zhǎng)巖有關(guān)，與閃長(zhǎng)巖關(guān)系更為密切。按成礦作用劃分，舒家店礦床類(lèi)型為斑巖型銅礦床[19, 22]。

2 AMT數(shù)據(jù)處理及分析

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

測(cè)區(qū)簡(jiǎn)化地質(zhì)圖與音頻大地電磁測(cè)點(diǎn)分布如圖1所示，在測(cè)區(qū)內(nèi)共部署音頻大地電磁測(cè)點(diǎn)60個(gè)。數(shù)據(jù)采集使用加拿大鳳凰公司的MTU-5A系統(tǒng)，配備AMTC-30磁傳感器。同時(shí)對(duì)電場(chǎng)和磁場(chǎng)進(jìn)行測(cè)量，共有6套儀器同時(shí)采集，采集頻率范圍為0.35~10 400 Hz。每個(gè)測(cè)點(diǎn)的數(shù)據(jù)采集時(shí)間大于60 min，濾波頻率設(shè)為50 Hz。利用高精度GPS實(shí)時(shí)同步測(cè)量平面坐標(biāo)與高程。利用MTU-5A系統(tǒng)配套的SSMT2000數(shù)據(jù)處理軟件進(jìn)行處理，并對(duì)部分含強(qiáng)干擾的數(shù)據(jù)進(jìn)行時(shí)間域?yàn)V波處理[23?26]和Rhoplus擬合[27]。依據(jù)DZ/T 0173—1997以及SY/T 5820—1999標(biāo)準(zhǔn)對(duì)AMT數(shù)據(jù)進(jìn)行質(zhì)量評(píng)價(jià)。結(jié)果表明數(shù)據(jù)質(zhì)量符合規(guī)范要求，可以進(jìn)行反演解釋。經(jīng)過(guò)處理后的典型數(shù)據(jù)(部分頻率)如圖2中AMT數(shù)據(jù)三維的方形點(diǎn)和圓形點(diǎn)所示。

2.2 維性分析

在進(jìn)行3D反演之前，需對(duì)數(shù)據(jù)進(jìn)行維性分析，以判別3D反演的必要性，并對(duì)觀測(cè)數(shù)據(jù)進(jìn)行定性分析。電性結(jié)構(gòu)維性判別有很多方法，主要有：基于阻抗張量的方法如阻抗極化橢圓分析[28]和基于相位張量的方法[29]?；谙辔粡埩康姆椒梢杂行У叵郎\層電性不均勻體的影響，從而可呈現(xiàn)可靠的深部電性結(jié)構(gòu)的維性。

圖3所示為舒家店礦床60個(gè)測(cè)點(diǎn)在頻率分別為 3 600，900，115和13.7 Hz相位張量橢圓圖。從圖3可見(jiàn)：整體上說(shuō)，舒家店礦床部分測(cè)點(diǎn)在部分頻率下呈現(xiàn)明顯三維構(gòu)造，表明了本文進(jìn)行3D反演的必要性；而在局部，維性分析的結(jié)果呈現(xiàn)了本區(qū)深部地質(zhì)體構(gòu)造分布的定性特征。一維電性特征常對(duì)應(yīng)于橫向均勻分布的地質(zhì)體，如水平地層和巖席、水平巖墻或其他均質(zhì)侵入體；二維電性特征常對(duì)應(yīng)于構(gòu)造走向明確的地質(zhì)體，如褶皺樞紐水平的沉積巖地層、片理方向穩(wěn)定的變質(zhì)巖等；三維電性特征常對(duì)應(yīng)于無(wú)明確構(gòu)造走向的地質(zhì)體，褶皺樞紐傾伏或直立、疊加褶皺、斷裂錯(cuò)距校大的地層或非均質(zhì)地質(zhì)體等，如混合巖化巖體、雜巖體、復(fù)合巖體等。

舒家店斑巖銅礦床位于永村橋—舒家店背斜的軸線(xiàn)上，受控于志留系地層和舒家店閃長(zhǎng)斑巖體[19?20, 22]。平面上，AMT測(cè)線(xiàn)測(cè)點(diǎn)也是分布以上2套地質(zhì)體之上。結(jié)合測(cè)點(diǎn)分布圖(圖1)，從測(cè)點(diǎn)相位張量橢圓分布圖(圖3)可以看出：在高頻段(3 600 Hz)，受地形的影響，淺地表三維特征明顯；在中高頻段(900 Hz)，除在電性梯度帶外，漸顯一維特征，反映物性在不同段趨于單一；在中低頻段(115 Hz)，除個(gè)別頻點(diǎn)受噪聲影響外，一維特征非常明顯，說(shuō)明在此深度范圍內(nèi)物性整體趨于單一；在低頻段(13.7 Hz)，二維特征增強(qiáng)，反映在此深度物性可能明顯區(qū)別于上部地質(zhì)體(115 Hz)。

(a) 測(cè)點(diǎn)S4114視電阻率；(b) 測(cè)點(diǎn)S4115視電阻率；(c) 測(cè)點(diǎn)S4114相位；(d) 測(cè)點(diǎn)S4115相位；(e) 測(cè)點(diǎn)S0412視電阻率；(b) 測(cè)點(diǎn)S4118視電阻率；(c) 測(cè)點(diǎn)S0412相位；(d) 測(cè)點(diǎn)S4118相位

3 AMT數(shù)據(jù)3D反演

3.1 3D反演參數(shù)

3D反演初始模型選擇電阻率為100 Ω·m的均勻半空間。圖4所示為局部坐標(biāo)系中的3D反演網(wǎng)格示意圖。坐標(biāo)原點(diǎn)為60個(gè)AMT測(cè)點(diǎn)的中心點(diǎn)。在局部坐標(biāo)系中，方向垂直向下，方向仍為NS向，方向仍為EW向。圖4(a)中，網(wǎng)格的生成嚴(yán)格滿(mǎn)足測(cè)點(diǎn)位于四邊形單元的中心，并且相鄰測(cè)點(diǎn)必須被2個(gè)四邊形單元隔離。這2種要求保證了計(jì)算的反演算法中正演的精度[4]。沿方向有73個(gè)網(wǎng)格，沿方向有58個(gè)網(wǎng)格，沿方向有20個(gè)網(wǎng)格(＞0)，空氣中網(wǎng)格數(shù)目為7(＜0)。因此，模型網(wǎng)格數(shù)為=73×58×27= 114 318個(gè)。

反演迭代目標(biāo)均方根誤差(RMS)設(shè)定為1。圖5所示為反演迭代過(guò)程中得到的幾類(lèi)模型參數(shù)。其中，模型范數(shù)表征模型的粗糙程度，均方根誤差表征觀測(cè)數(shù)據(jù)的擬合程度，為L(zhǎng)agrange 乘子。數(shù)據(jù)空間Occam方法的處理步驟[4]是：首先通過(guò)改變使均方根誤差達(dá)到最小值，然后在保持均方根誤差的同時(shí)，通過(guò)改變尋找模型范數(shù)最小的模型。從圖5可見(jiàn)：經(jīng)過(guò)20次迭代后，均方根誤差變化不大。第16次迭代的模型范數(shù)為1 630，均方根誤差為6.136(=1)，該模型均方根誤差接近最小，模型范數(shù)最小。

頻率/Hz：(a)3 600.0；(b) 900.0；(c) 115.0；(d) 13.7

(a) z=0切面；(b) x=0的切面

(a) 模型范數(shù)；(b) 均方根誤差；(c) Lagrange乘子

3.2 3D反演模型

為評(píng)價(jià)3D反演所得模型的精度，可將3D反演模型的響應(yīng)(偏對(duì)角阻抗元素計(jì)算的視電阻率和相位數(shù)據(jù))與實(shí)測(cè)數(shù)據(jù)結(jié)果進(jìn)行對(duì)比，部分結(jié)果如圖2所示?？傮w上，數(shù)據(jù)擬合較好，表明反演結(jié)果相對(duì)可靠，所得模型可被接受，并可作為進(jìn)一步分析解釋的依據(jù)。

圖6所示為舒家店礦床測(cè)區(qū)音頻大地電磁3D反演模型沿著幾條典型測(cè)線(xiàn)的垂直電阻率切面圖。需要說(shuō)明的是，成圖時(shí)，對(duì)3D反演模型進(jìn)行了插值處理，插值方法是空間3D距離倒數(shù)法。從圖6可以看出：縱向上，總體呈現(xiàn)為3層結(jié)構(gòu)，各測(cè)線(xiàn)500 m以?xún)?nèi)覆蓋有低阻層，SJD04線(xiàn)在淺部的低阻背景中含有3個(gè)明顯的高阻異常體；在500~1 700 m處為一相對(duì)高阻層；在1 700~2 000 m處則主要表現(xiàn)為相對(duì)中阻的基底。橫向上，淺部的覆蓋層呈現(xiàn)高、低阻相間特征，但差異較??；而500~1700 m范圍內(nèi)的高、低阻分段特征更為明顯，且相間分布的高、低阻體均展現(xiàn)出較明顯的垂直條帶狀特征；1 700 m以下整體為中低阻背景，其間分布有與上部相貫通的高阻體。

圖7所示為3D反演模型在不同深度的平面切面。從圖6和圖7可知：1) 100 m以淺以低阻覆蓋為主，局部分布有零星的高阻體；2) 100~600 m內(nèi)高低阻分區(qū)域相間分布，整體電阻率隨深度增加逐步升高； 3) 600~1 700 m為相對(duì)中高阻層，高、低阻的分布在橫向上表現(xiàn)出一定的分塊特征；4) 在1 700~2 000 m，電阻率相對(duì)降低，高阻體的范圍縮小。

圖6 舒家店礦床測(cè)區(qū)音頻大地電磁3D反演模型的2D切面

圖7 舒家店礦床測(cè)區(qū)音頻大地電磁3D反演模型的不同深度平面切面

4 討論

在對(duì)地下電性結(jié)構(gòu)所代表的地質(zhì)意義進(jìn)行分析之前，必然要先確定影響電性結(jié)構(gòu)的主控因素。舒家店礦區(qū)地表巖性較簡(jiǎn)單，主要為志留系砂頁(yè)巖和輝石閃長(zhǎng)巖，礦體賦存于巖體之中而巖體的圍巖為志留系地層。由已往舒家店礦區(qū)的物性測(cè)量經(jīng)驗(yàn)可知志留系砂頁(yè)巖電阻率變化較大，但相對(duì)于巖體其電阻率平均值總體偏低，約為2 500 Ω·m，而巖體平均值可達(dá)14 000 Ω·m[18]。值得注意的是，相對(duì)形成于燕山期的巖體，志留系及其下伏地層則經(jīng)歷更長(zhǎng)地質(zhì)歷史時(shí)間的改造而更易形成錯(cuò)綜復(fù)雜的斷裂系統(tǒng)，地質(zhì)流體的下滲充填更易進(jìn)一步降低地層的電阻率。因此，這2種主要地質(zhì)單元之間顯著的物性差異是電磁法探測(cè)的有利基礎(chǔ)。盡管如此，由于礦床在形成過(guò)程中同時(shí)產(chǎn)生了典型的斑巖型銅礦蝕變分帶(蝕變類(lèi)型主要有鈣硅酸鹽化、鉀化、青盤(pán)巖化、長(zhǎng)石分解蝕變，局部分育矽卡巖化)[20]，而巖石的蝕變必然帶來(lái)物性上的劇烈變化，不同的蝕變類(lèi)型和蝕變程度直接影響巖石的電阻率，目前尚無(wú)針對(duì)蝕變巖石的詳細(xì)物性測(cè)定，這給地質(zhì)解釋帶來(lái)一定的困難。對(duì)蝕變影響的正確分析是突破對(duì)地球物理找礦認(rèn)識(shí)限制的必經(jīng)之路。

分析蝕變規(guī)律發(fā)現(xiàn)：鉀化主要發(fā)育在巖體的深部，淺部不太發(fā)育，其蝕變主要有3種形式，分別為彌漫狀、脈狀及脈體暈，主要蝕變生成礦物為鉀長(zhǎng)石、黑云母等含鉀礦物；而青磐巖化主要分布在巖體的中部，主要生成礦物為綠簾石+綠泥石，其中綠泥石化發(fā)育的范圍較大；長(zhǎng)石分解蝕變(石英—絹云母—綠泥石—黏土化)則疊加在新鮮巖石及早期蝕變組合上，蝕變強(qiáng)度也最大[19]。由此可知，3次主要的蝕變過(guò)程生成的礦物主體為陽(yáng)離子Fe和Na的含水礦物，而這類(lèi)礦物的電阻率相對(duì)蝕變?cè)瓗r中的主要礦物(斜長(zhǎng)石、鉀長(zhǎng)石、石英)要低幾個(gè)數(shù)量級(jí)[30]。而最后一期、強(qiáng)度最大的長(zhǎng)石分解蝕變疊加在早期的蝕變之上更是生成了大量的含水礦物(相對(duì)于白云母含更多水的絹云母、綠泥石、黏土礦物)。在將巖石作為均質(zhì)并且各向同性的多相多晶集合體的Hashin-Shtrikman模型中[31?32]，代入各礦物的電阻率和體積分?jǐn)?shù)可估算出這類(lèi)蝕變會(huì)顯著降低地質(zhì)體的電阻率[33?34]。

與此同時(shí)，由于成礦流體中大量H2O和其他氣體(如CO2)的加入，蝕變通常并非等體積反應(yīng)，蝕變后的地質(zhì)體均存在一定的膨脹系數(shù)，結(jié)合與礦化蝕變有直接成因聯(lián)系的巖體的早期隱爆作用、成礦期后的冷卻收縮至裂作用、礦石礦物的沉淀等過(guò)程，在后期地表水的下滲過(guò)程中均起到降低相應(yīng)地質(zhì)體電阻率的作用。

基于上述討論，結(jié)合舒家店礦床的實(shí)際勘探認(rèn) 識(shí)[35?41]，舒家店銅礦床的三維電性結(jié)構(gòu)具有以下 特征：

1) 測(cè)區(qū)北西和南東部在近地表(?500 m以上)的低阻層代表了志留系地層的分布和延伸。

2) 測(cè)區(qū)中部沿北東—南西走向的近地表高阻體代表了未蝕變的新鮮閃長(zhǎng)類(lèi)巖體；存在低阻體，其原因可能是巖體內(nèi)存在頂垂體或巨型包體、蝕變巖體或礦體。

3) 深部的高阻體基本代表了未蝕變巖體的范圍，但由于志留系中較高電阻率的石英砂巖與新鮮閃長(zhǎng)類(lèi)巖體不存在顯著的物性差別，故不排除高阻體為產(chǎn)狀較陡立的志留系的可能。

4) 深部的低阻體尤其是圍繞高阻體的低阻體異常可能為蝕變巖體或礦體。因?yàn)樵谏畈繋r體演變成孤立的巖株?duì)罡怕瘦^小，這類(lèi)與高阻巖體共存的低阻體便極可能是在巖體蝕變過(guò)程中形成，甚至直接為礦體的反映。

總體而言，3D模型勾畫(huà)了銅陵地區(qū)的空間電性框架。但須注意的是：參與3D反演的實(shí)測(cè)AMT測(cè)點(diǎn)偏少，頻率稀，空間分布不均勻，且數(shù)據(jù)本身存在一定的觀測(cè)誤差，加之反演精度有限，所得到的3D電性模型較粗糙，所得結(jié)果可能存在偏差。

5 結(jié)論

1) 舒家店礦床地下電性主要呈3D分布，進(jìn)行3D反演是必要的。

2) 利用WSINV3DMT得到舒家店測(cè)區(qū)的3D反演電性模型。模型響應(yīng)結(jié)果與實(shí)測(cè)結(jié)果擬合效果較好，表明該反演模型整體合理，可作為解釋的依據(jù)。

3) 舒家店礦床深部與代表巖體的高阻體相間共存的低阻體異常為蝕變巖體或礦體。

4) 所得舒家店礦區(qū)的3D電性模型可為認(rèn)識(shí)舒家店礦床的成礦作用和成礦規(guī)律提供參考。由于參與3D反演實(shí)測(cè)AMT的數(shù)據(jù)測(cè)點(diǎn)少，頻率稀，空間分布不均勻，且數(shù)據(jù)本身存在一定的觀測(cè)誤差，所得3D電性模型較粗糙，難以提供更加精確的地質(zhì)信息，因此，本次得到的3D電性模型仍有待完善，在數(shù)據(jù)密度、反演方法、參數(shù)設(shè)計(jì)、模型約束以及地質(zhì)認(rèn)識(shí)等方面需進(jìn)一步研究。

致謝 安徽省地質(zhì)調(diào)查院對(duì)本文研究提供了幫助，Siripunwaraporn教授提供了WSINV3DMT反演代碼，在此一并致謝！

[1] SMITH R. Electromagnetic induction methods in mining geophysics from 2008 to 2012[J]. Surveys in Geophysics, 2014, 35(1): 123?156.

[2] Chen X B, LV Q T, YAN J Y. 3D electrical structure of porphyry copper deposit:a case study of Shaxi copper deposit[J]. Applied Geophysics, 2012, 9(3): 270?278.

[3] 湯井田, 任政勇, 周聰, 等. 淺部頻率域電磁勘探方法綜述[J]. 地球物理學(xué)報(bào), 2015, 58(8): 2681?2705. TANG Jingtian, REN Zhengyong, ZHOU Cong, et al. Frequency-domain electromagnetic methods for exploration of the shallow subsurface:a review[J]. Chinese Journal of Geophysics, 2015, 58(8): 2681?2705.

[4] SIRIPUNVARAPORN W. Three-dimensional magnetotelluric inversion: an introductory guide for developers and users[J]. Surveys in Geophysics, 2012, 33(1): 5?27.

[5] HAN N, NAM M J, KIM H J, et al. Three-dimensional inversion of magnetotelluric data including sea effects obtained in Pohang,Korea[J]. Journal of Applied Geophysics, 2009, 68(4): 533?545.

[6] NEWMAN G A, ALUMBAUGH D L. Three-dimensional magnetotelluric inversion using non-linear conjugate gradients[J]. Geophysical Journal International, 2000, 140(2): 410?424.

[7] ZHDANOV M S, FANG S, HURSáN G. Electromagnetic inversion using quasi-linear approximation[J]. Geophysics, 2000, 65(5): 1501?1513.

[8] 譚捍東, 余欽范, BOOKER J, 等. 大地電磁法三維快速松弛反演[J]. 地球物理學(xué)報(bào), 2003, 46(6): 850?855. TAN Handong, YU Qinfan, BOOKER J, et al. Three-dimensional rapid relaxation inversion for the magnetotelluric method[J]. Chinese J Geophys, 2003, 46(6): 850?855.

[9] SIRIPUNVARAPORN W, EGBERT G, LENBURY Y, et al. Three-dimensional magnetotelluric inversion: data-space method[J]. Physics of the Earth and Planetary Interiors, 2005, 150(1): 3?14.

[10] SIRIPUNVARAPORN W, EGBERT G. WSINV3DMT: vertical magnetic field transfer function inversion and parallel implementation[J]. Physics of the Earth and Planetary Interiors, 2009, 173(3/4): 317?329.

[11] XIAO Q B, CAI X P, XU X W, et al. Application of the 3D magnetotelluric inversion code in a geologically complex area[J]. Geophysical Prospecting, 2010, 58(6): 1177?1192.

[12] 常印佛, 劉湘培, 吳言昌. 長(zhǎng)江中下游銅鐵成礦帶[M]. 北京: 地質(zhì)出版社, 1991: 349?361.CHANG Yinfo, LIU Xiangpei, WU Yanchang. The copper-iron belt of the lower and middle reaches of the Yangtze River[M]. Beijing: Geological Press, 1991: 349?361.

[13] 翟裕生, 姚書(shū)振, 林新多, 等. 長(zhǎng)江中下游地區(qū)鐵銅(金)成礦規(guī)律[M]. 北京: 地質(zhì)出版社, 1992(1): 1?12.ZHAI Yusheng, YAO Shuzhen, LIN Xinduo, et al. The metallogenic pattern of iron-copper (gold) in the lower Yangtze River belt[M]. Beijing: Geological Publish House, 1992(1): 1?12.

[14] 毛景文, 邵擁軍, 謝桂青, 等. 長(zhǎng)江中下游成礦帶銅陵礦集區(qū)銅多金屬礦床模型[J]. 礦床地質(zhì), 2009, 28(2): 109?119.MAO Jingwen, SHAO Yongjun, XIE Guiqing, et al. Mineral deposit model for porphyry-skarn polymetallic copper deposits in Tongling ore dense district of Middle-Lower Yangtze Valley metallogenic belt[J]. Mineral Deposits, 2009, 28(2): 109?119.

[15] 毛景文, HOLLY S, 杜安道, 等.長(zhǎng)江中下游地區(qū)輝鉬礦Re–Os年齡測(cè)定及其對(duì)成礦作用的指示[J]. 地質(zhì)學(xué)報(bào), 2004, 78(1): 121?131.MAO Jingwen, HOLLY S, DU Andao, et al. Molybdenite Re-Os precise dating for molybdenite from Cu-Au-Mo Deposits in the middle-lower reaches of Yangtze River belt and its implications for mineralization[J]. Acta Geologica Sinica, 2004, 78(1): 121?131.

[16] 儲(chǔ)國(guó)正. 長(zhǎng)江中下游地區(qū)成礦地質(zhì)背景分析[J]. 安徽地質(zhì), 2003, 13(1): 34?43.CHU Guozheng. Analysis of ore-forming geological background in the middle and lower reaches of the Yangtze River[J]. Geology of Anhui, 2003, 13(1): 34?43.

[17] 唐永成, 吳言昌, 儲(chǔ)國(guó)正, 等. 安徽沿江地區(qū)銅金多金屬礦床地質(zhì)[M]. 北京: 地質(zhì)出版社, 1998: 83?98.TANG Yongcheng, WU Yanchang, CHU Guozheng, et al. The polymetal deposit geology of Cu and Au in the region along the river of Anhui[M]. Beijing: Geological Publishing House, 1998: 83?98.

[18] 王世偉, 周濤發(fā), 袁峰, 等.銅陵舒家店斑巖銅礦成礦年代學(xué)研究及其成礦意義[J]. 巖石學(xué)報(bào), 2012, 28(10): 3170?3180.WANG Shiwei, ZHOU Taofa, YUAN Feng, et al. Re-os and 40Ar/39Ar dating of the Shujiadian copper deposit in Tongling, China: implications for regional metallogenesis[J]. Acta Petrologica Sinica, 2012, 28(10): 3170?3180.

[19] WANG S, ZHOU T, YUAN F, et al. Geological and geochemical studies of the Shujiadian porphyry Cu deposit, Anhui Province, Eastern China: Implications for ore genesis[J]. Journal of Asian Earth Sciences, 2014, 103: 252?275.

[20] 胡英才, 李桐林, 范翠松, 等.安徽銅陵舒家店銅礦的電磁法試驗(yàn)研究[J]. 地質(zhì)學(xué)報(bào), 2014, 88(4): 612?619.HU Yincai, LI Tonglin, FAN Cuisong, et al. Experiment research of electromagnetic exploration method in Shujiadian Copper Deposit,Tongling,Anhui Province[J]. Acta Geologica Sinica, 2014, 88(4): 612?619.

[21] 徐曉春, 白茹玉, 謝巧勤, 等.安徽銅陵中生代侵入巖地質(zhì)地球化學(xué)特征再認(rèn)識(shí)及成因討論[J]. 巖石學(xué)報(bào), 2012, 28(10): 3139?3169.XU Xiaochun, BAI Ruyu, XIE Qiaoqing, et al. Re-understanding of the geological and geochemical characteristics of the Mesozoic intrusive rocks from Tongling area of Anhui Province, and discussions on their genesis[J]. Acta Petrologica Sinica, 2012, 28(10): 3139?3169.

[22] 賴(lài)小東, 楊曉勇, 孫衛(wèi)東, 等. 銅陵舒家店巖體年代學(xué)、巖石地球化學(xué)特征及成礦意義[J]. 地質(zhì)學(xué)報(bào), 2012, 86(3): 470?485.LAI Xiaodong, YANG Xiaoyong, SUN Weidong, et al. Chronological-geochemical characteristics of the Shujiadian intrusion, Tongling ore cluster field: its significance to metaliogenesis[J]. Acta Geologica Sinica, 2012, 86(3): 470?485.

[23] 湯井田, 李晉, 肖曉, 等. 數(shù)學(xué)形態(tài)濾波與大地電磁噪聲壓制[J]. 地球物理學(xué)報(bào), 2012, 55(5): 1784?1793.TANG Jingtian, LI Jin, XIAO Xiao, et al. Mathematical morphology filtering and noise suppression of magnetotelluric sounding data[J]. Chinese Journal of Geophysics, 2012, 55(5): 1784?1793.

[24] 湯井田, 蔡劍華, 任政勇, 等. Hilbert-Huang變換與大地電磁信號(hào)的時(shí)頻分析[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 35(5): 1399?1405.TANG Jingtian, CAI Jianhua, REN Zhengyong, et al. Hilbert-Huang transform and time-frequency analysis of magnetotelluric signal[J]. Journal of Central South University (Science and Technology), 2009, 35(5): 1399?1405.

[25] 湯井田, 李晉, 肖曉, 等. 基于數(shù)學(xué)形態(tài)濾波的大地電磁強(qiáng)干擾分離方法[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 43(6): 2215?2221.TANG Jingtian, LI Jin, XIAO Xiao, et al. Magnetotelluric sounding data strong interference separation method base on mathematical morphology filtering[J]. Journal of Central South University (Science and Technology), 2012, 43(6): 2215?2221.

[26] 湯井田, 劉子杰, 劉峰屹, 等. 音頻大地電磁法強(qiáng)干擾壓制試驗(yàn)研究[J]. 地球物理學(xué)報(bào), 2015, 58(12): 4636?4647. TANG Jingtian, LIU Zijie, LIU Fengyi, et al. The denoising of the audio magnetotelluric data set with strong interferences[J]. Chinese Journal of Geophysics, 2015, 58(8): 4636?4647.

[27] 周聰, 湯井田, 任政勇, 等. 音頻大地電磁法“死頻帶”畸變數(shù)據(jù)的Rhoplus校正[J]. 地球物理學(xué)報(bào), 2015, 58(12): 4648?4660.ZHOU Cong, TANG Jingtian, REN Zhengyong, et al. Application of the Rhoplus method to audio magnetotelluric dead band distortion data[J]. Chinese Journal of Geophysics, 2015, 58(8): 4648?4660.

[28] BECKEN M, BURKHARDT H. An ellipticity criterion in magnetotelluric tensor analysis[J]. Geophysical Journal International, 2004, 159(1): 69?82.

[29] CALDWELL T G, BIBBY H M, BROWN C. The magnetotelluric phase tensor[J]. Geophysical Journal International, 2004, 158(2): 457?469.

[30] 王彪. 舒家店銅礦床地質(zhì)地球化學(xué)特征及成因分析[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010, 33(6): 906?910.WANG Biao. Geological and geochemical characters of Shujiadian copper deposit and genesis analysis[J]. Journal of Hefei University of Technology, 2010, 33(6): 906?910.

[31] HASHIN Z, SHTRIKMAN S. A variational approach to the theory of the effective magnetic permeability of multiphase materials[J]. Journal of Applied Physics, 1962, 33(10): 3125?3131.

[32] HASHIN Z, SHTRIKMAN S. A variational approach to the theory of elastic behavior of multiphase materials[J]. Journal of the Mechanics and Physics of Solids, 1963, 11(2): 127?140.

[33] JONES A G, EVANS R L, EATON D W.Velocity-conductivity relationships for mantle mineral assemblages in Archean cratonic lithosphere based on a review of laboratory data and Hashin—Shtrikman external bounds[J]. Lithos, 2009, 109(1/2): 131?143.

[34] FULLEA J, MULLER M R, JONES A G. Electrical conductivity of continental lithospheric mantle from integrated geophysical and petrological modeling: application to the Kaapvaal Craton and Rehoboth Terrane, southern Africa[J]. Journal of Geophysical Research, 2011, 116(B10): B10202.

[35] 呂慶田, 楊竹森, 嚴(yán)加永, 等.長(zhǎng)江中下游成礦帶深部成礦潛力、找礦思路與初步嘗試: 以銅陵礦集區(qū)為實(shí)例[J]. 地質(zhì)學(xué)報(bào), 2007, 81(7): 865?881.Lü Qingtian, YANG Zhuseng, YAN Jiayong, et al. The metallogenic potential,prospecting idea and primary attempt in depth of the ore belt of the middle and lower reach of the Yangtze River:a case study of Tongling Ore District[J]. Acta Geologica Sinica, 2007, 81(7): 865?881.

[36] Lü Q T, HOU Z Q, ZHAO J H, et al. Deep seismic reflection profiling reveals complex crustal structure of Tongling ore district[J]. Science in China: Series D, 2004, 47(3): 193?200.

[37] 嚴(yán)加永, 呂慶田, 孟貴祥, 等. 銅陵礦集區(qū)中酸性巖體航磁3D成像及對(duì)深部找礦方向的指示[J]. 礦床地質(zhì), 2009, 28(6): 838?849.YAN Jiayong, Lü Qingtian, MENG Guixiang, et al. Aeromagnetic 3D inversion imaging for intermediate-acid intrusive bodies and its indication significance of deep ore prospecting in Tongling ore concentration district[J]. Mineral Deposits, 2009, 28(6): 838?849.

[38] DENG Jun, WANG Qingfei, XIAO Changhao, et al. Tectonic-magmatic-metallogenic system, Tongling ore cluster region, Anhui Province, China[J]. International Geology Review, 2011, 53(5/6): 449?476.

[39] TANG J T, ZHOU C, WANG X Y, et al. Deep electrical structure and geological significance of Tongling ore district[J]. Tectonophysics, 2013, 606(23): 78?96.

[40] 湯井田, 周聰, 任政勇, 等.安徽銅陵礦集區(qū)大地電磁數(shù)據(jù)三維反演及其構(gòu)造格局[J]. 地質(zhì)學(xué)報(bào), 2014, 88(4): 598?611.TANG Jingtian, ZHOU Cong, REN Zhengyong, et al. Three dimensional magnetotelluric inversion and structural framework of Tongling Ore District, Anhui[J]. Acta Geologica Sinica, 2014, 88(4): 598?611.

[41] 王顯瑩, 湯井田, 張林成, 等. 長(zhǎng)江中下游成礦帶中段巖石圈電性結(jié)構(gòu)研究[J]. 地球物理學(xué)報(bào), 2015, 58(12): 4403?4414.WANG Xianyin, TANG Jingtian, ZHANG Lincheng, et al. Lithospheric electrical structure in the middle and lower reach of Yangtze River metallogenic belt inferred from magnetotelluric sounding[J]. Chinese Journal of Geophysics, 2015, 58(8): 4403?4414.

(編輯 陳燦華)

Three dimensional audio-frequency magnetotellurics prospecting for copper deposit in Shujiadian, Tongling, Anhui Province

WANG Yongqing1, 2, TANG Jingtian1, 2, XIAO Xiao1, 2, ZHOU Cong1, 2, WANG Xianying1, 2, REN Zhengyong1, 2

(1.Key Laboratory of Metallogenic Prediction of Nonferrous Metals, Ministry of Education, Changsha 410083, China;2. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China)

Audio-frequency magnetotelluric method and 3D inversion code were used to study the porphyry copper deposit in Shujiadian, Tongling, Anhui Province. A 3D inversion was demanded. Sites data were selected and frequencies were decided, and the WSINV3DMT code was used with a reasonable parameter design to obtain the 3D conductivity structure. The electrical characteristics and structure framework of the 3D model were discussed by slices of the model with different depths and orientations. The results show that 2D/3D characteristics for the deep structure of the region are obvious. The model is proved to be high in precision by comparison between the observed data and the model response. Combined with the discussion of properties of the main geological body, the deep low resistivity anomaly between the high resistivity rock mass may be the reflection of altered rock or the ore body.

porphyry copper deposit; audio-frequency magnetotellurics; 3D inversion; Tongling ore district

P33

A

1672?7207(2017)04?1018?09

10.11817/j.issn.1672?7207.2017.04.023

2016?06?12；

2016?08?25

國(guó)家深部探測(cè)專(zhuān)項(xiàng)(SinoProbe-3)；國(guó)家自然科學(xué)基金資助項(xiàng)目(41574120, 41174105)；國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目(2014AA06A602)(Project(SinoProbe-3) supported by the Key Program of Deep Exploration in China; Projects(41574120, 41174105 ) supported by the National Natural Science Foundation of China; Project(2014AA06A602) supported by the National High Technology Research and Development Program(863 Program) of China)

肖曉，博士，副教授，從事地球物理信號(hào)處理及反演成像等研究；E-mail：csuxiaox@csu.edu.cn