萬 偉 王志剛 魯 瑤

(①東華理工大學地球物理與測控技術(shù)學院，江西南昌 330013；②東方地球物理公司綜合物化探處，河北涿州 072751)

0 引言

頻率域可控源電磁法(controlled source electromagnetic，CSEM)通過人工發(fā)射不同頻率的電磁波探測不同深度的地層，相對天然源電磁法具有更強的抗電磁干擾能力和更高分辨率，在油氣與金屬礦產(chǎn)資源勘查、工程地質(zhì)勘察等領(lǐng)域發(fā)揮了重要作用[1-3]。根據(jù)不同的應(yīng)用場景，頻域CSEM分為地面可控源電磁法與海洋可控源電磁法。地面可控源電磁法通常在地表采用固定的單源，在相對較寬的頻段范圍(10-2～104Hz)內(nèi)發(fā)射數(shù)十個不同頻率的電磁波，通過電磁場或視電阻率反演地電結(jié)構(gòu)；海洋可控源電磁法的發(fā)射源通常由船拖拽沿測線移動，在相對較窄頻段范圍(10-1～101Hz)內(nèi)發(fā)射數(shù)個不同量級頻率的電磁波，接收器由船拖拽或置于海底進行數(shù)據(jù)采集[4-5]。海洋可控源電磁法采用移動觀測源，可獲得電場振幅隨偏移距變化的觀測資料，當海底存在相對高阻的油氣儲層時，能夠觀測到明顯的電磁異常，目前已廣泛應(yīng)用于海洋油氣勘探[6-8]。

除海洋可控源電磁法采用移動源觀測方式外，人工源電法如直流電阻率法、激發(fā)極化法等[9-11]，也是通過改變源與測點間的幾何關(guān)系約束地下空間的電性信息，本質(zhì)上也是一種移動源的觀測方式。此外，在長偏移距瞬變電磁法(LOTEM)基礎(chǔ)上發(fā)展起來的多通道瞬變電磁法[12]，采用陣列式布設(shè)，在每個測點上記錄多個偏移距下的時間域響應(yīng)，通過時間域多次疊加和空間域多次覆蓋提高地下結(jié)構(gòu)探測的縱向分辨率[13]，在油氣藏監(jiān)測及礦產(chǎn)勘查中得到了成功應(yīng)用[14-15]。魯瑤等[16]提出了移動源瞬變電磁成像的概念，從空間和時間兩個方面對觀測數(shù)據(jù)進行多次覆蓋和多次疊加，很好地解決了反演非唯一性的問題。瞬變電磁法通過采用不同偏移距的移動源觀測，明顯提高了對地電結(jié)構(gòu)的分辨力。

對于地面頻域可控源電磁法，可控源音頻大地電磁法(controlled-source audio magnetotlleuric, CSAMT)是最早應(yīng)用于實際探測的可控源電磁法[17]。該方法主要基于頻率測深原理，在遠區(qū)觀測相互正交的水平電場與磁場，依據(jù)計算的卡尼亞視電阻率反演地下電性結(jié)構(gòu)。但由于發(fā)射功率的限制，在保證接收點信號強度前提下，即使盡可能增加收發(fā)距，CSAMT觀測資料中仍同時包含近區(qū)場、中間區(qū)場和遠區(qū)場，理論計算和實際觀測均顯示卡尼亞視電阻率在近區(qū)與中間區(qū)會產(chǎn)生畸變，很大程度上降低了CSAMT的探測力。為改善CSAMT探測效果，何繼善[18]提出基于單分量定義全區(qū)視電阻率，并提出廣域電磁法，可同時利用近區(qū)場、中間區(qū)場及遠區(qū)場的觀測資料，相對于傳統(tǒng)CSAMT觀測系統(tǒng)更靈活，顯著拓寬了地面頻域可控源電磁法的應(yīng)用范圍。同時，為提高可控源電磁法的探測精度，王志剛等[19]和何展翔等[20-22]綜合時域與頻域電磁法的優(yōu)點，提出時頻電磁法，該方法廣泛應(yīng)用于陸地油氣勘探及水力壓裂監(jiān)測。趙國澤等[23]和卓賢軍等[24]在高電阻率地區(qū)建立固定大功率發(fā)射臺，利用大功率人工源發(fā)射極低頻電磁波增加地面可控源電磁法的探測深度。此外，為解決因直線源極化存在觀測弱區(qū)的問題，將發(fā)射源由線源變?yōu)椤癓”形源[25]或“十”形源[26]，以拓寬地面可控源電磁法的觀測區(qū)域。

目前，各種地面頻域可控源電磁法均只采集單個偏移距的電磁場數(shù)據(jù)，若要充分利用電磁數(shù)據(jù)進行頻率測深，需盡可能增加收發(fā)距，而采用過大收發(fā)距易導致接收信號振幅降低，特別是在人文噪聲干擾嚴重的測區(qū)，數(shù)據(jù)信噪比無法滿足處理解釋的要求；若采用小收發(fā)距，觀測資料會過早進入近區(qū)場，此時僅依靠單偏移距的觀測資料無法同時約束不同深度的目標[27]。目前地面可控源電磁法的實際應(yīng)用效果不如海洋可控源電磁法，主要原因是地面可控源電磁觀測更復雜，如地表淺部非目標電性體、不可控人文噪聲及發(fā)射源的布設(shè)等，這些因素的疊加增加了地面可控源電磁法最終探測效果的不確定性[28]。常規(guī)地面可控源電磁法只利用單偏移距的數(shù)據(jù)，在信噪比與探測深度之間較難平衡，也給實際觀測的發(fā)射源位置選擇增加了難度。

針對目前單偏移距地面頻率可控源電磁法存在的上述問題，本文從獲得較高信噪比觀測資料與較大探測深度兩個目的出發(fā)，一方面考慮采用小收發(fā)距獲得高信噪比觀測資料，另一方面采用多偏移距的移動源觀測，通過充分利用幾何測深與頻率測深達到增加探測深度的目的。為檢驗其探測效果，在前人研究基礎(chǔ)上，首先基于有代表性的一維地電模型，對單偏移距與多偏移距移動源可控源電磁法的探測效果進行系統(tǒng)性對比，然后對贛東北某礦區(qū)開展實測數(shù)據(jù)方法實驗，進一步驗證地面移動源電磁法對于解決實際問題的適用性。

1 地面移動源頻率域電磁法及數(shù)據(jù)反演

1.1 地面移動源頻率域電磁法觀測系統(tǒng)

地面移動源頻率域電磁法的觀測系統(tǒng)如圖1所示，以沿x軸方向、離測點相對較近且與測點共線的水平接地電偶極源為發(fā)射信號，觀測時發(fā)射源沿x軸移動，分別在數(shù)個不同收發(fā)距(偏移距)下激發(fā)10-2～103Hz頻段范圍內(nèi)數(shù)十個不同頻率的電磁波，對應(yīng)在測點采集多頻率、多偏移距的水平電場分量Ex。圖1同時展示了傳統(tǒng)可控源電磁法(如CSAMT)的觀測系統(tǒng)，以沿y軸的水平接地電偶極源發(fā)射電磁波，一般情況下源相對較長，如圖中的固定源1，在離發(fā)射源相對較遠的地面同時觀測相互正交的水平電場分量Ey和磁場分量Hx，觀測頻率與前述移動源頻率域電磁法的觀測頻率一致。

圖1 固定源與移動源觀測系統(tǒng)示意圖

基于Ey與Hx可計算阻抗視電阻率

(1)

式中：ω=2πf為角頻率，其中f為頻率；μ0為真空磁導率。式(1)利用Ey與Hx的比值計算阻抗視電阻率，具有一定的壓制噪聲作用，但需注意的是，只有電磁場為平面波場或遠場時，才能用式(1)計算ρa。因此CSAMT數(shù)據(jù)觀測要求發(fā)射源距離測區(qū)足夠遠，使所有頻率的電磁波盡可能為平面波場。

1.2 地面移動源頻率域電磁法探測的物理基礎(chǔ)

在CSAMT實際數(shù)據(jù)觀測中，由于發(fā)射源功率的限制，并不能使所有觀測頻率下的電場均達到平面波場的條件，大部分接收到的低頻電磁場為非平面波場或近區(qū)場，若依據(jù)式(1)計算ρa，計算結(jié)果會在這些頻段失真。以中間層為低阻的三層電阻率模型(圖2)為例，分別計算不同收發(fā)距r時的視電阻率ρa，計算結(jié)果見圖2。水平電偶源在層狀電阻率模型中的電磁場可通過半解析方法求得，再依據(jù)式(1)可計算ρa。

由圖2可見：當收發(fā)距r為無窮大時，頻段10-2～103Hz內(nèi)接收的電磁場均為平面波場，對應(yīng)的阻抗視電阻率ρa(圖2中虛線)與模型基本吻合；收發(fā)距r為12km時，對應(yīng)的阻抗視電阻率見圖2中的黑色實線，可見頻率低于30Hz時，ρa畸變明顯，這是因為低頻段(近場)的觀測電磁場不滿足式(1)的計算條件(即平面波場)。針對這個問題，若直接利用場分量進行反演，則無需計算ρa，可避免近場導致ρa畸變的問題。但場分量相對于ρa更易受噪聲干擾，特別是常規(guī)可控源電磁法采用大收發(fā)距觀測，電磁場振幅相對較小，極易受噪聲干擾，因此縮小收發(fā)距可顯著增大觀測點接收的電磁場振幅，提高信噪比。

圖2 三層電阻率模型及不同收發(fā)距r時計算的視電阻率

以圖2電阻率模型為例，采用常規(guī)可控源電磁法觀測，分別計算不同收發(fā)距下各頻率的Ey振幅(圖3)，可見收發(fā)距從8km減小至2km，整體上Ey振幅提高了約2個數(shù)量級。同時，由于低頻段(近場)電磁場振幅不隨頻率變化，因而低頻段數(shù)據(jù)不再具有頻率測深功能，雖然縮短收發(fā)距顯著提高了電場振幅，但也導致觀測頻段內(nèi)更多的低頻資料進入近場，如圖3中陰影部分：在低頻段(近場)Ey振幅不隨頻率變化，且隨著收發(fā)距減小，低頻段(近場)包含的頻段范圍變寬，意味著更多的低頻數(shù)據(jù)進入近場，同時能夠被用于頻率測深的數(shù)據(jù)也相應(yīng)減少。此時，若僅依靠頻率測深，則有效探測深度將減小，同時無法發(fā)揮近場數(shù)據(jù)的作用。

此外，從圖3還可以看出，近場的電場振幅雖然不隨頻率變化，但對收發(fā)距的變化敏感。基于圖2模型，設(shè)收發(fā)距為2km，分別計算不同中間層電阻率(ρ2=10、100、1000Ω·m)所對應(yīng)的Ey振幅，結(jié)果如圖4所示。由圖可見，近場的Ey振幅對地下電阻率的變化較敏感，即近場的Ey振幅與收發(fā)距和地下電阻率有關(guān)，與頻率無關(guān)。根據(jù)近場的這種特性，本文提出在多頻率觀測基礎(chǔ)上，不斷移動場源，利用改變源與接收點間的幾何位置關(guān)系，獲得不同收發(fā)距、不同頻率的電場響應(yīng)，通過充分利用幾何測深與頻率測深，理論上將有助于提高地面人工源頻率域電磁法的探測效果。

圖3 圖2模型不同收發(fā)距r對應(yīng)的Ey振幅

圖4 圖2模型不同中間層電阻率下的Ey振幅(r=2km)

1.3 合成數(shù)據(jù)及反演方法

為系統(tǒng)對比移動源電磁法與傳統(tǒng)可控源電磁法的探測效果，對典型的一維電阻率模型的合成數(shù)據(jù)進行反演，合成數(shù)據(jù)共5組，來自不同的發(fā)射源、觀測參數(shù)和添加的高斯噪聲水平組合(表1)。第1組數(shù)據(jù)采用CSAMT觀測方式(圖1中固定源1，沿y方向布設(shè))，利用較遠處觀測的Ey和Hx，基于式(1)計算得到ρa，并加載5%的高斯噪聲，得到合成數(shù)據(jù)；第2組數(shù)據(jù)的觀測系統(tǒng)與第1組數(shù)據(jù)相同，區(qū)別是對Ey加載了5%高斯噪聲；第3組數(shù)據(jù)與第1、第2組數(shù)據(jù)的觀測系統(tǒng)相同，區(qū)別是考慮到當收發(fā)距較大、觀測環(huán)境中噪聲較嚴重的情況，對Ey加載了20%的高斯噪聲；第4組數(shù)據(jù)采用圖1中固定源2的布設(shè)方式(沿x軸方向并與測點共線)，采用相對較小收發(fā)距(偏移距)，對Ex加載5%高斯噪聲；第5組數(shù)據(jù)來源于移動源觀測系統(tǒng)(圖1中移動源)，設(shè)計了一系列離測點較近、沿x軸方向并與測點在同一條線上的移動電偶極源，采集Ex，并對Ex加載5%高斯噪聲。這5組數(shù)據(jù)的觀測頻段均為10-2～103Hz，對數(shù)等間隔取40個頻點。后文將對這5組數(shù)據(jù)的反演結(jié)果進行對比分析。

表1 一維CSEM合成數(shù)據(jù)參數(shù)表

對表1的數(shù)據(jù)組反演，本文均采用Occam反演方法[29-30]，該反演方法基于正則化思想[31]，在擬合模型響應(yīng)數(shù)據(jù)和實測數(shù)據(jù)的同時，采用最平滑模型保證反演算法的穩(wěn)定性。反演目標函數(shù)為

(2)

式中：m=[m1,m2,…,mM]為模型參數(shù)向量，M表示向量維度；d=[d1,d2,…,dN]為數(shù)據(jù)向量，N表示向量維度；μ為正則化因子；F(m)為模型m的正演響應(yīng)；W表示數(shù)據(jù)協(xié)方差矩陣；?m表示粗糙度矩陣，其定義式有多種方式，這里采用模型的一階粗糙度[30]

(3)

為有效對比不同方法的探測效果，所有觀測數(shù)據(jù)的反演條件保持一致，除了采用平滑模型約束以外，不添加其他約束條件。同時，所有數(shù)據(jù)反演采用的初始模型均為電阻率100Ω·m的均勻半空間。在各數(shù)據(jù)反演迭代過程中，設(shè)定迭代停止條件為擬合誤差rms≤1或前后兩次迭代的rms的差值小于1×10-4。

2 理論模型試驗

對地下電性異常體(層)的識別效果是檢驗電磁法探測效果的重要標準之一。通常情況下，傳統(tǒng)地面可控源電磁法主要基于頻率測深，對低阻目標較敏感，但容易受地表低阻覆蓋層的影響。海洋可控源電磁法主要基于幾何測深，受低阻海水層影響相對小，對高阻目標層的識別效果相對較好[32]。因而，借鑒海洋可控源電磁法的工作原理，地面移動源電磁法是否能通過綜合利用幾何測深與頻率測深有效識別低阻和高阻異常目標層呢？針對這個問題，這里考慮覆蓋層和目標層電阻率不同的情況，建立了表2所示的四個一維地電模型：高阻覆蓋層+高阻目標層(模型1)，高阻覆蓋層+低阻目標層(模型2)，低阻覆蓋層+高阻目標層(模型3)，及低阻覆蓋層+低阻目標層(模型4)。這四個模型的基底電阻率均為200Ω·m，地電模型具體參數(shù)見圖5中黑色曲線。針對這四個模型，分別討論基于表1中5組數(shù)據(jù)對目標層的探測能力。

表2 一維層狀電阻率模型參數(shù)

2.1 傳統(tǒng)地面可控源電磁數(shù)據(jù)反演結(jié)果

針對上述四個模型，采用傳統(tǒng)地面可控源電磁法觀測方式得到一維合成數(shù)據(jù)(表1中數(shù)據(jù)1～數(shù)據(jù)3)，分別對其進行反演，反演結(jié)果見圖5，對應(yīng)的數(shù)據(jù)擬合及誤差rms見圖6。

由圖5和圖6可見：當噪聲干擾較小時，電場Ey(數(shù)據(jù)1)的反演結(jié)果明顯優(yōu)于ρa(數(shù)據(jù)2)，更接近實際模型，且近場畸變導致ρa在低頻段嚴重失真，致使數(shù)據(jù)擬合效果不理想；當Ey噪聲干擾嚴重(數(shù)據(jù)3)時，即使數(shù)據(jù)能夠較好地擬合，最終的反演結(jié)果與實際模型仍存在較大偏差。此外，無論Ey還是ρa的反演結(jié)果均表明，與低阻目標層模型的反演結(jié)果相比，對高阻目標層的識別效果較差。

2.2 小收發(fā)距可控源電磁法數(shù)據(jù)反演結(jié)果

根據(jù)前文分析，相對于ρa，利用Ey進行地下電阻率反演的效果雖然有一定改善，但由于采用較大收發(fā)距，易導致數(shù)據(jù)信噪比低，影響探測效果。若采用小收發(fā)距，雖可提高信噪比，但是否會降低對目標層的識別力？改變收發(fā)距對探測效果有何影響？針對前面四種不同模型，分別開展小收發(fā)距可控源電磁法數(shù)據(jù)反演測試。

考慮到小收發(fā)距下沿發(fā)射源方向電流密度更高，在觀測時設(shè)定發(fā)射源沿x軸方向，與接收點共線，即對應(yīng)表1中的數(shù)據(jù)4。發(fā)射源長度均為0.1km，分別考慮收發(fā)距為1.2、2.4、3.6、4.8km的情況。對正演模擬Ex加載5%高斯噪聲得到合成數(shù)據(jù)(數(shù)據(jù)4)，模型1～模型4對應(yīng)的數(shù)據(jù)4的反演結(jié)果見圖7。由圖可見：偏移距為1.2km時，所有模型合成數(shù)據(jù)的反演結(jié)果均不能有效反映目標層及圍巖電性特征；偏移距增加到2.4km時，反演結(jié)果有明顯改善；當偏移距增大到3.6km時，所有模型的反演結(jié)果均能反映明顯的電阻率異常，但高阻異常與實際高阻目標層深度不符；隨著偏移距進一步增大到4.8km，各模型反演結(jié)果與偏移距為3.6km時相比沒有明顯變化，并與大收發(fā)距中數(shù)據(jù)2的反演結(jié)果相似。綜合上述結(jié)果分析，對于小收發(fā)距可控源電磁觀測方式，偏移距的變化對探測結(jié)果影響較大，過小的偏移距會嚴重影響探測結(jié)果準確性，恰當?shù)钠凭嗫捎行Ц纳品囱萁Y(jié)果。結(jié)合大收發(fā)距條件下數(shù)據(jù)2的反演結(jié)果(圖5)可知，若采用單發(fā)射源，當偏移距達到一定值后，繼續(xù)增加偏移距對探測結(jié)果改善效果是有限的。需說明的是，此例小收發(fā)距下響應(yīng)數(shù)據(jù)與合成數(shù)據(jù)均擬合較好，對應(yīng)各組數(shù)據(jù)反演的rms誤差值均達到1，因篇幅限制，這里未展示數(shù)據(jù)擬合結(jié)果。

圖5 模型1(a)、模型2(b)、模型3(c)、模型4(d)及對應(yīng)表1中數(shù)據(jù)1～數(shù)據(jù)3的反演結(jié)果

圖6 對應(yīng)圖5中模型1(a)、模型2(b)、模型3(c)、模型4(d)分別對數(shù)據(jù)1(左)、2(中)、3(右)的反演數(shù)據(jù)及擬合誤差rms

圖7 模型1(a)、模型2(b)、模型3(c)、模型4(d)在不同偏移距下對數(shù)據(jù)4的反演結(jié)果

2.3 移動源電磁觀測數(shù)據(jù)反演結(jié)果

針對移動源電磁觀測系統(tǒng)，分析其探測效果，模擬數(shù)據(jù)即表1中的數(shù)據(jù)5。在小收發(fā)距條件下觀測，設(shè)定發(fā)射源長為0.1km，沿x軸方向移動，保持源與測點共線，逐步增大收發(fā)距(偏移距)，可獲得不同偏移距的觀測資料。在移動源觀測系統(tǒng)中，根據(jù)發(fā)射源移動距離、移動次數(shù)的不同，可得到不同偏移距組合的觀測資料，對應(yīng)的探測效果可能存在差異。針對不同偏移距組合情況，這里分別進行討論。

首先，針對模型1～模型4，考慮源移動距離相對較大、源分別移動1次和2次的情況，對應(yīng)的偏移距組合共3組：組合1，源移動1次，偏移距依次為1.2、3.6km；組合2，源移動1次，偏移距依次為2.4、3.6km；組合3，源移動2次，偏移距依次為1.2、2.4、3.6km。以rms達到1為響應(yīng)數(shù)據(jù)與模型數(shù)據(jù)的擬合標準，這些組合的各模型合成數(shù)據(jù)反演結(jié)果見圖8。由圖可見，組合1～組合3的合成數(shù)據(jù)反演結(jié)果均準確地反映了實際模型的電阻率特征。對比圖8與圖5、圖7，可以看出移動源電磁法的探測效果明顯優(yōu)于傳統(tǒng)可控源電磁法和小收發(fā)距可控源電磁法，尤其是對高阻異常的識別效果有明顯提高。此外，對于本文所討論的四個模型，圖8中組合1～組合3的合成數(shù)據(jù)反演結(jié)果相近，說明當源移動距離相對較大時，源只需移動1次(獲得2個不同偏移距的觀測資料)，便可發(fā)揮移動源電磁法準確探測目標的優(yōu)勢。

其次，當源移動距離相對小時，源移動前、后測點電場振幅的相對變化也較小，若電場振幅的相對變化小到被噪聲掩蓋，此時可視為不同偏移距下觀測的電場振幅一致，理論上與單偏移距觀測的探測效果無區(qū)別，這種情況下開展多偏移距移動源電磁法探測就沒有意義。為詳細地討論源移動距離相對較小情況下的探測效果，針對模型1～模型4，分別設(shè)計了4種偏移距組合：組合4，偏移距依次為2.40、2.41km；組合5，偏移距依次為2.40、2.42km；組合6，偏移距依次為2.40、2.51km；組合7，偏移距依次為2.40、2.52km。經(jīng)過計算，組合4與組合5的源移動前后，測點電場振幅的相對變化范圍為1%～3%；組合6與組合7的源移動前后，測點電場振幅的相對變化均略大于10%。對于含5%高斯噪聲的合成數(shù)據(jù)(數(shù)據(jù)5)，組合4與組合5的源移動前后，測點電場振幅的相對變化被噪聲掩蓋，而組合6與組合7未被噪聲掩蓋。組合4～組合7的各模型合成數(shù)據(jù)反演結(jié)果見圖9。對比圖8與圖9，可見組合6與組合7的合成數(shù)據(jù)反演結(jié)果與組合1～組合3相近，其效果明顯優(yōu)于組合4和組合5。這進一步說明，要發(fā)揮移動源電磁法探測的優(yōu)勢，源移動的距離不能過小，需保證源移動前、后測點電場振幅的相對變化不被噪聲掩蓋。在實際觀測中，環(huán)境噪聲水平通常是未知的，所以在保證接收點信號強度的前提下，建議源移動距離應(yīng)盡可能大一些。

圖8 模型1(a)、模型2(b)、模型3(c)、模型4(d)在偏移距組合1～組合3下對數(shù)據(jù)5的反演結(jié)果

圖9 模型1(a)、模型2(b)、模型3(c)、模型4(d)在偏移距組合4～組合7下對數(shù)據(jù)5的反演結(jié)果

2.4 結(jié)果分析

結(jié)合典型一維模型，比較分析了傳統(tǒng)可控源電磁法及移動源電磁法的探測效果，結(jié)果表明傳統(tǒng)可控源電磁法若采用視電阻率進行反演，受近區(qū)畸變影響嚴重，探測效果較差；若直接采用場分量進行反演，可在一定程度上降低近區(qū)場的影響，但大收發(fā)距會降低接收電場分量的信噪比；通過減小收發(fā)距可提高資料信噪比，但對高阻目標體的識別不夠準確。利用本文提出的移動源電磁法，對各種地電模型均有較好的探測效果，相對于傳統(tǒng)可控源電磁法，能夠更準確地反映目標層電阻率特征，特別是對高阻目標的探測效果提升明顯。

3 實測數(shù)據(jù)分析

為驗證移動源電磁法探測效果，本文結(jié)合實測數(shù)據(jù)展開進一步分析。試驗區(qū)位于贛東北地區(qū)，屬萬載—杭州灣銅金鉛鋅鉭多金屬成礦帶，萬年—德興銅金鉛鋅多金屬成礦亞帶。該地區(qū)也是中國主要的有色金屬和貴金屬成礦區(qū)，已探明的礦種有銅、鉛鋅、金、銀、鎢等多個礦種，這有利于開展人工源電磁法測量進行電性參數(shù)測定，探明深部電性異常特征及空間展布，為鉆探提供物探依據(jù)。

試驗區(qū)內(nèi)巖樣的電性測量結(jié)果為：砂質(zhì)千枚巖、凝灰質(zhì)含碳千枚巖等圍巖的電阻率相對較高，測量的巖樣電阻率值范圍為1550～2160Ω·m；斑巖型銅礦礦體呈低電阻率特征，測量的樣品電阻率值范圍為110～190Ω·m。圍巖與目標礦體(斑巖型銅礦)的電性差異明顯，為電磁勘探奠定了良好的物性基礎(chǔ)。

在工區(qū)開展了圖1中所示的固定源和移動源電磁觀測試驗。采用傳統(tǒng)的CSAMT觀測方式，發(fā)射源沿y方向布設(shè)，長度為1.87km，發(fā)射頻率為0.125～1024Hz，對數(shù)等間隔共40個頻點；為保證信號強度，采用了較大發(fā)射電流10A；測點位于發(fā)射源中垂線上，測點離發(fā)射源較遠，收發(fā)距為10.09km。對兩套數(shù)據(jù)分別進行反演：一是利用Ey與Hx計算得到ρa，對其進行反演；二是直接采用觀測的Ey進行反演。做為對比，開展了移動源觀測試驗，其中發(fā)射源長度為0.2km，沿x軸方向移動且與測點共線，發(fā)射源共移動2次，偏移距依次為2.40、2.51、2.61km，發(fā)射電流保持2A，發(fā)射頻率均為0.02～1024Hz，對數(shù)等間隔取35個頻點。對采集的Ex進行反演。針對移動源電磁法，反演數(shù)據(jù)分別來自兩組偏移距組合方式：組合1為包含3個偏移距(2.40、2.51、2.61km)的Ex觀測資料；組合2為包含2個偏移距(2.40、2.61km)的Ex觀測資料。

傳統(tǒng)可控源電磁法和移動源電磁法共四套實驗數(shù)據(jù)的反演結(jié)果見圖10，對應(yīng)的數(shù)據(jù)擬合結(jié)果見圖11?？梢?，采用CSAMT的視電阻率ρa進行反演時，由于視電阻率在低頻段發(fā)生畸變，反演電阻率在深部出現(xiàn)嚴重失真，模型響應(yīng)與觀測數(shù)據(jù)擬合效果較差；采用CSAMT的Ey直接進行反演，反演電阻率未體現(xiàn)深部(500～1000m)的低阻異常特征。與前兩組反演結(jié)果相比，對基于移動源電磁觀測數(shù)據(jù)進行反演，反演電阻率曲線不僅能反映淺部高阻層的存在，對深部(500～1000m)低阻異常的特征反映得更清楚。

圖10 贛東北地區(qū)實測數(shù)據(jù)電阻率反演結(jié)果

圖11 基于圖10中不同數(shù)據(jù)的反演電阻率模型對應(yīng)的模擬數(shù)據(jù)與觀測數(shù)據(jù)擬合曲線

4 討論與結(jié)論

常規(guī)地面可控源電磁法均采用單發(fā)射源，以頻率測深為主，在觀測數(shù)據(jù)的信噪比與探測深度之間較難平衡。本文提出在地面采用移動源電磁觀測，通過縮小收發(fā)距增加信噪比，同時通過移動發(fā)射源以獲得多個偏移距的觀測數(shù)據(jù)，實現(xiàn)綜合利用幾何測深與頻率測深的優(yōu)勢增加有效探測深度。結(jié)合幾種典型一維地電模型，采用OCCAM反演方法，通過系統(tǒng)對比移動源電磁法與傳統(tǒng)地面可控源電磁法的探測效果發(fā)現(xiàn)：對于大收發(fā)距的地面可控源電磁法，直接采用場分量反演的探測效果明顯優(yōu)于阻抗視電阻率反演，但在大收發(fā)距情況下，只利用場分量反演結(jié)果易受噪聲干擾，并最終影響探測效果，且無論是利用阻抗視電阻率還是場分量反演，均不能解決對高阻異常體不敏感的問題；對于移動源電磁觀測，不同地電模型的反演結(jié)果均優(yōu)于傳統(tǒng)地面可控源電磁法，無論是高阻異常體還是低阻異常體，反演結(jié)果與實際地電模型均吻合較好，且采用移動源電磁法觀測可增加有效探測深度。最后，實際實驗數(shù)據(jù)反演結(jié)果表明，地面移動源電磁法探測效果明顯優(yōu)于常規(guī)地面可控源電磁法。

需要說明的是，本文對地面移動源電磁法的探測效果展開了初步分析，提出采用多偏移距的觀測方式，其發(fā)射源的布設(shè)較靈活，對本文探討的相對簡單的地電模型，采用兩個不同偏移距的發(fā)射源即可明顯改善探測效果，但對于復雜的地電模型或探測對象，如何選擇恰當?shù)膮?shù)以達到最佳的探測效果，還有待進一步的研究。