CSAMT的多偶極子源特征與張量測(cè)量

王顯祥，底青云，許誠(chéng)

1 中國(guó)科學(xué)院地質(zhì)與地球物理研究所中國(guó)科學(xué)院工程地質(zhì)力學(xué)院重點(diǎn)實(shí)驗(yàn)室，北京 100029

2 中國(guó)科學(xué)院大學(xué)，北京 100049

1 引言

可控源音頻大地電磁法（Controlled Source audio－frequency Magnetotelluric，CSAMT）是一種從大地電磁法（Magnetotelluric，MT）發(fā)展出的人工源電磁測(cè)深法（何繼善，1990；樸化榮，1990；石昆法，1999；底青云和王若，2008）.該方法由于使用人工場(chǎng)源，具有抗干擾能力強(qiáng)、信噪比高的特點(diǎn).只需改變頻率，即可達(dá)到探測(cè)不同深度的目的，它具有迅速快捷的特點(diǎn).近年來(lái)，CSAMT在金屬礦、地?zé)帷⒌叵滤?、油氣的勘探中發(fā)揮著越來(lái)越重要的作用，已經(jīng)成為一種重要的地球物理勘探方法（底青云等，2001，2002a，2002b；孫英勛，2005；柳建新等，2009；譚儒蛟等，2007；李帝銓，2008）.

然而在CSAMT測(cè)量中，現(xiàn)有的發(fā)射源大都為單電偶極子源，從理論計(jì)算的輻射花樣圖中，可以看到電磁場(chǎng)各分量在360°張角范圍內(nèi)都存在一定的弱區(qū)（何繼善，1990；底青云和王若，2008）.為了保證數(shù)據(jù)質(zhì)量，測(cè)量區(qū)域有著嚴(yán)格的限制，一般要求在偶極源中軸線60°扇形區(qū)域內(nèi)測(cè)量，偶極源長(zhǎng)度在1～3km范圍內(nèi)，除此之外受波區(qū)條件的限制，收發(fā)距一般為10～15km.在崎嶇的山區(qū)，有時(shí)很難滿足這些苛刻條件，從而無(wú)法進(jìn)行正常的測(cè)量工作，即使勉強(qiáng)工作，也很難保證數(shù)據(jù)質(zhì)量.為此，本文發(fā)展了由兩個(gè)偶極子源組成的“L”型裝置源，并在發(fā)射信號(hào)時(shí)采用了一種新的技術(shù)，使得各個(gè)分量在360°張角范圍內(nèi)不再存在弱區(qū)，這也突破了傳統(tǒng)測(cè)量方式中只在偶極源中軸線60°夾角范圍內(nèi)測(cè)量的局限.

CSAMT測(cè)量精度除受發(fā)射源裝置的限制外，還受測(cè)量方式的影響.現(xiàn)有測(cè)量方式大都為標(biāo)量測(cè)量，只觀測(cè)Ex、Hy兩個(gè)分量.嚴(yán)格意義上說(shuō)，這種測(cè)量方式只適用于一維和當(dāng)測(cè)量方向與構(gòu)造方向垂直的二維情況.當(dāng)不滿足以上兩種情況時(shí)，特別是地質(zhì)體為三維或者測(cè)量方向與構(gòu)造方向相差很大時(shí)，直接用標(biāo)量數(shù)據(jù)進(jìn)行反演，必然對(duì)反演結(jié)果帶來(lái)一定的誤差.在以上情況下，有效的測(cè)量方式應(yīng)為張量測(cè)量（Boerner et al.，1993；Bromley，1993；Lu et al.，1997；Wannamaker，1997a，1997b；Li et al.，2000；Garcia et al.，2003），這種測(cè)量方式不僅可以得到任意測(cè)量方向上的視電阻率，還可以提供傾子、二維判別、橢率等一系列評(píng)判地下構(gòu)造復(fù)雜性的參數(shù)，通過(guò)這些參數(shù)可以對(duì)地質(zhì)結(jié)構(gòu)有個(gè)比較系統(tǒng)、詳細(xì)的了解 （陳 樂(lè) 壽 等，1990；Raiche and Hohmann，1974；Homman，1975；Weidelt，1975；Ting et al.，1981；Tripp and Hohmann，1984）.為此，在“L”型源的基礎(chǔ)上，本文也發(fā)展了相關(guān)的張量測(cè)量方法.

本文首先推導(dǎo)了“L”型裝置源中各分量的表達(dá)式，驗(yàn)證了公式的正確性，并給出了“L”型源各分量的輻射花樣圖，通過(guò)輻射花樣圖的分析，本文又開發(fā)了一種全新的信號(hào)發(fā)射模式；最后給出一個(gè)張量測(cè)量的實(shí)例，說(shuō)明了張量測(cè)量方式的優(yōu)勢(shì).

2 “L”型裝置源

2.1 模型建立與公式推導(dǎo)

為克服原有單偶極子源布極的缺點(diǎn)，本文設(shè)計(jì)了新的“L”型源，這種源的布極方式如圖1所示，圖中箭頭的方向代表某一時(shí)刻電流的方向，相當(dāng)于原有兩個(gè)相互垂直的單電偶極源組合而成.電偶極源1在 （x0，y0）點(diǎn)產(chǎn)生的電磁場(chǎng)強(qiáng)度為、，電偶極源2在 （x0，y0）產(chǎn)生的電磁場(chǎng)強(qiáng)度為，根據(jù)電磁場(chǎng)的疊加原理，可求出“L”型源在 （x0，y0）產(chǎn)生的電磁場(chǎng)強(qiáng)度為

現(xiàn)把“L”型源簡(jiǎn)化成如圖2所示，假設(shè)源1與源2都可看作電偶極子源，長(zhǎng)度分別為a、b，兩者夾角90°，并共中心點(diǎn)，發(fā)射電流大小都為I.坐標(biāo)方向定義沿源1方向?yàn)閤軸方向，源2方向?yàn)閥軸方向，根據(jù)電偶極子源在均勻半空間的表達(dá)式，它們?cè)冢▁0，y0）處的電場(chǎng)強(qiáng)度分別為（何繼善，1990）

式中上標(biāo)1、2分別代表源1與源2，下標(biāo)x、y分別代表沿x與y方向，I為發(fā)射電流，ρ為均勻半空間電阻率，r為收發(fā)距，k為波數(shù)，θ為r與源1構(gòu)成的夾角.設(shè)α角為：

圖2 “L”型源矢量合成示意圖Fig.2 The sketch map of vector synthesis for“L”type source

根據(jù)（1）式，現(xiàn)計(jì)算源1與源2在 （x0，y0）處，分別沿α和α＋90°方向的電場(chǎng)強(qiáng)度：

電磁場(chǎng)遵從矢量合成原則，即多個(gè)偶極子源可以根據(jù)矢量合成的原則進(jìn)行疊加，疊加后的偶極子源在某點(diǎn)產(chǎn)生的電磁場(chǎng)等于這些偶極子源分別在該點(diǎn)產(chǎn)生的電磁場(chǎng).根據(jù)矢量合成原則，源1與源2合成源長(zhǎng)度應(yīng)為與源1夾角為α，根據(jù)單偶極子源關(guān)系式同樣可得

2.2 模型驗(yàn)證

為了驗(yàn)證公式和程序的正確性，特做了均勻半空間100Ωm的模型做對(duì)比，利用視電阻率公式

圖3 均勻半空間上兩種源得到的視電阻對(duì)比圖Fig.3 Comparison diagram of apparent resistivity for two different sources in homogeneous half space

分別計(jì)算單偶極子源和“L”型裝置源收發(fā)距10km時(shí)，視電阻率從1Hz～104Hz的變化情況，如圖3所示.從圖中可以看出，除在低頻段電阻率值有稍許差別外，兩種裝置源所求得到的視電阻率完全相同，這說(shuō)明本文求“L”型裝置源各分量的公式和程序是正確、可靠的，同時(shí)也表明“L”型裝置源可以很好地應(yīng)用到CSAMT測(cè)量當(dāng)中.

2.3 數(shù)值模擬

圖4為電場(chǎng)矢量圖，在圖中其長(zhǎng)度代表電場(chǎng)強(qiáng)度，方向代表了電場(chǎng)方向.圖4a為長(zhǎng)度2km電偶極子源電流矢量圖，從圖中可以看到電場(chǎng)從源的一頭流向另一頭，沿弧形方向流動(dòng)；電場(chǎng)強(qiáng)度在源附近最強(qiáng).圖4b、圖4c為由兩個(gè)長(zhǎng)度2km，夾角90°的電偶極子源組成的“L”型裝置源的電流矢量圖，圖中的紅線代表實(shí)際的電偶極子源，箭頭代表電流方向，根據(jù)矢量疊加原則，可疊加成圖中的粉色源.從圖中可以看到，無(wú)論當(dāng)兩個(gè)電流源方向首尾相接時(shí)，還是兩個(gè)電流源方向相反時(shí)，電流都會(huì)從源的一頭流向另一頭，且中軸線附近集中，軸向稀疏，這些特征都與單偶極子源相符.

圖4 偶極子源與“L”型裝置源電場(chǎng)矢量對(duì)比圖（a）偶極子源電場(chǎng)矢量圖；（b）“L”型裝置源電場(chǎng)矢量圖樣式一；（c）“L”型裝置源電場(chǎng)矢量圖樣式二.圖中的紅色實(shí)線代表真實(shí)存在的偶極子源，粉線代表合成的偶極子源，箭頭代表了偶極子源方向.Fig.4 Comparison diagram of current vector for two different sources in homogeneous half space（a）The diagram of current vector for dipole source；（b）The diagram of one kind of current vector for“L”type source；（c）The diagram of the other kind of current vector for“L”type source.

圖5 偶極子源與“L”型裝置源電場(chǎng)輻射花樣對(duì)比圖（a）lgEx，f＝90Hz；（b）lgEy，f＝90Hz；（c）lgE45°，f＝90Hz；（d）lgE135°，f＝90Hz.圖中的白色實(shí)線代表真實(shí)存在的偶極子源，粉線代表合成的偶極子源，箭頭代表了偶極子源方向.Fig.5 Comparison diagram of radiation patterns for two different sources in homogeneous half spaceThe red solid line represents the true dipole source，the pink line represents the synthesis dipole source，arrowhead represents the direction in the figure.

圖5 為偶極子源與“L”型裝置源電場(chǎng)輻射花樣對(duì)比圖，圖中的白線為實(shí)際的電偶極子源，粉線為應(yīng)合成的源，其中圖5（a、b）為單偶極子源輻射花樣圖，圖5（c、d）為“L”型裝置源輻射花樣圖.通過(guò)比較不難發(fā)現(xiàn)，圖5a與圖5d、圖5b與圖5c除在源附近形態(tài)有些不一致外，整體形態(tài)幾乎完全一致，圖5d和圖5c就像圖5a和圖5b整體往左旋轉(zhuǎn)45°一樣.場(chǎng)源附近不一致的原因，主要是因?yàn)樵?.1節(jié)中推導(dǎo)的公式是在假設(shè)源為偶極子的前提下，當(dāng)離源非常近時(shí)兩個(gè)源不能看作兩個(gè)偶極子，而是每個(gè)偶極子都要被看作由n個(gè)偶極子組成，所以這個(gè)區(qū)域的場(chǎng)是由多個(gè)偶極子合成的，場(chǎng)自然也變得比較復(fù)雜.所幸的時(shí)，CSAMT測(cè)量主要在遠(yuǎn)場(chǎng)區(qū)進(jìn)行，近場(chǎng)區(qū)的形態(tài)對(duì)測(cè)量結(jié)果不會(huì)造成影響.

2.4 “L”型裝置源的應(yīng)用

在上文中，從理論和數(shù)值模擬結(jié)果上分別驗(yàn)證了偶極子源滿足矢量合成的原理.這并不是本文的目的，本文的目的主要是利用偶極子源的矢量合成理論改善現(xiàn)有CSAMT測(cè)量中存在的問(wèn)題.“L”型裝置源的第一個(gè)應(yīng)用就是改變現(xiàn)有CSAMT測(cè)量只能在一定張角范圍內(nèi)，而不能在360°張角范圍內(nèi)測(cè)量的現(xiàn)狀；第二個(gè)應(yīng)用主要是用來(lái)減少布設(shè)偶極子源的次數(shù)，矢量測(cè)量要測(cè)量Ex、Ey、Hx、Hy、Hz5個(gè)分量，由偶極子源輻射花樣圖可知，Ex、Hy場(chǎng)強(qiáng)時(shí)，Ey、Hx場(chǎng)弱，為了滿足信噪比的對(duì)等可能要布設(shè)兩個(gè)相互垂直的偶極子源，而張量測(cè)量則必須布設(shè)兩個(gè)方向偶極子源，每個(gè)方向上測(cè)量5個(gè)分量共計(jì)10個(gè)分量，而通過(guò)“L”型裝置源的組合可以輕松減少一次布設(shè)；第三個(gè)應(yīng)用主要在崎嶇山區(qū)，眾所周知，現(xiàn)有CSAMT測(cè)量要滿足很多條件，既要使測(cè)量方向與地質(zhì)體走向垂直，又要選擇合適的收發(fā)距，太小會(huì)落入近區(qū)，太大又使信號(hào)太弱，降低信噪比，在崎嶇山區(qū)，這些苛刻條件有時(shí)是很難滿足的.

圖6 “L”型源一個(gè)周期內(nèi)兩個(gè)不同時(shí)刻的輻射花樣（a，b上半周期；c，d下半周期）（a）lgE45°，f＝90Hz；（b）lgE135°，f＝90Hz；（c）lgE45°，f＝90Hz；（d）lgE135°，f＝90Hz.圖中的白色實(shí)線代表真實(shí)存在的偶極子源，粉線代表合成的偶極子源，箭頭代表了偶極子源方向.Fig.6 Comparison diagram of radiation patterns for“L”type source at various timesThe red solid line represents the true dipole source，the pink line represents the synthesis dipole source，arrowhead represents the direction in the figure.

當(dāng)“L”型源兩個(gè)偶極子源的組合類型如圖1a時(shí)，輻射花樣如圖6（a、b）所示，組合類型如圖1b時(shí)，輻射花樣如圖6（c、d）所示.通過(guò)比較不難發(fā)現(xiàn)，圖6（a、b）和圖6（c、d）有很強(qiáng)的互補(bǔ)性，以45°方向的分量為例進(jìn)行說(shuō)明，在圖6a中電場(chǎng)強(qiáng)的區(qū)域，在圖6c中較弱，反之在圖6c中電場(chǎng)強(qiáng)的區(qū)域在圖6a中較弱.為了充分利用兩種裝置源的互補(bǔ)性，本文重新設(shè)計(jì)了發(fā)射信號(hào)，如圖7a所示，在圖7a的第1個(gè)周期內(nèi)，裝置源相當(dāng)于圖1b源，當(dāng)在第2個(gè)周期時(shí)，源相當(dāng)于圖1a源.現(xiàn)把兩個(gè)周期合并成一個(gè)周期發(fā)射信號(hào)，就可以滿足在一個(gè)周期內(nèi)360°張角范圍內(nèi)，Ex、Ey、Hx、Hy都會(huì)出現(xiàn)最強(qiáng)時(shí)刻，矢量測(cè)量時(shí)只要測(cè)量對(duì)應(yīng)時(shí)刻最強(qiáng)的分量即可，這樣就可以在一個(gè)周期內(nèi)測(cè)量到Ex、Ey、Hx、Hy4個(gè)分量，從而保證了CSAMT矢量測(cè)量信噪比的對(duì)等；張量測(cè)量時(shí)只要把上、下兩個(gè)半周期看作不同的偶極子源，分別測(cè)量5個(gè)分量即可.與傳統(tǒng)要布設(shè)兩次偶極源而實(shí)現(xiàn)CSAMT矢量或者張量測(cè)量相比，本方法具有以下優(yōu)點(diǎn)：因?yàn)榭梢栽谝粋€(gè)周期內(nèi)測(cè)量，而不是分兩次、在不同地方布設(shè)發(fā)射源，可以減少外界隨機(jī)噪音的干擾，還可以減少陰影效應(yīng)、復(fù)印效應(yīng)對(duì)解釋結(jié)果的影響，同時(shí)還減少了人力、物力和時(shí)間的消耗.在改變電流方向的同時(shí)，為了不增加發(fā)射機(jī)的數(shù)目，本文設(shè)計(jì)了如圖7b所示的改變電流方向的裝置.＋號(hào)源選擇性接1、2號(hào)接線柱，－號(hào)源接剩下的另一個(gè)接線柱，就可以達(dá)到改變電流的目的，當(dāng)然要實(shí)現(xiàn)自動(dòng)改變的目標(biāo)還要加一些輔助切換接線柱的設(shè)備.

圖7 “L”型源發(fā)射信號(hào)方式示意圖（a）“L”型源中2個(gè)電偶極源發(fā)射信號(hào)示意圖；（b）改變電流方向裝置示意圖.Fig.7 The sketch map of transmission mode for“L”type source（a）The sketch map of transmission mode of two electrical dipole sources in“L”type source；（b）The sketch map of device to change the current direction.

圖9 偶極源合成示意圖Fig.9 The sketch of composition of two sources

前文中主要討論的是源1與源2之間的夾角為90°時(shí)的情況，實(shí)際上可以把它們之間的夾角推廣到任意角度，仍滿足矢量疊加原則.圖8為兩偶極子源角度呈120°時(shí)的輻射花樣圖，根據(jù)矢量疊加原則合成源應(yīng)沿x方向.從Ex、Ey、Hx、Hy輻射花樣圖上可以看到，其輻射花樣圖與沿x方向電偶極子源輻射花樣圖相比，除在近場(chǎng)區(qū)差別較大外，其他區(qū)域并沒(méi)有太大差別，由此可見其滿足矢量合成原則.實(shí)際上也可以將這種裝置源應(yīng)用到傳統(tǒng)標(biāo)量測(cè)量上，如圖9所示，要在崎嶇山區(qū)布設(shè)源1，由于受地形限制無(wú)法鋪設(shè)，反而有2、3所示兩個(gè)方向的線路可以布設(shè).根據(jù)矢量合成原理，就可以把源分別布設(shè)在2、3線方向上，以達(dá)到鋪設(shè)在方向1上的目的，這樣做雖然損失了部分能量，但卻順利完成了工作.

圖8 “L”型源120°夾角時(shí)輻射花樣圖（a）lgEx，f＝90Hz；（b）lgEy，f＝90Hz；（c）lg Hx，f＝90Hz；（d）lg Hy，f＝90Hz.圖中的白色實(shí)線代表真實(shí)存在的偶極子源，粉線代表合成的偶極子源，箭頭代表了偶極子源方向.Fig.8 Comparison diagram of radiation patterns for"L"type source while two sources at 120°angleThe red solid line represents the true dipole source，the pink line represents the synthesis dipole source，arrowhead represents the direction in the figure.

3 張量測(cè)量

為了說(shuō)明張量測(cè)量的好處，本文設(shè)計(jì)了如圖10所示的模型，源為由兩個(gè)長(zhǎng)度分別為2km的偶極子組成的“L”型裝置源，背景為三層結(jié)構(gòu)，電阻率分別為100、10、1000Ωm，厚度分別為200、100m，在第一層中有一個(gè)邊長(zhǎng)為100m的立方體，電阻率為10Ωm.收發(fā)距為10km，在異常體上方布設(shè)了5條測(cè)線，每條測(cè)線長(zhǎng)400m，線距50m，點(diǎn)距25m.

3.1 構(gòu)造方向的確定

現(xiàn)假設(shè)地質(zhì)結(jié)構(gòu)成二維性，如圖11所示，當(dāng)測(cè)量方向與構(gòu)造主軸方向夾角為任意θ角時(shí)，由電磁場(chǎng)旋轉(zhuǎn)關(guān)系可知：

傳統(tǒng)的標(biāo)量測(cè)量方式中，視電阻為

由公式（13）、（14）可知，當(dāng)測(cè)量方向與主軸方向不同時(shí)，會(huì)帶來(lái)不同的視電阻率偏差，由于只測(cè)量Ex和Hy分量，只能得到視電阻率ρxy，無(wú)法進(jìn)行地質(zhì)體維數(shù)和偏離角θ的判別，如果盲目地進(jìn)行一維或二維的反演，必然會(huì)引起不小的誤差.

圖10 單個(gè)異常體模型示意圖（a）模型斷面圖；（b）觀測(cè)裝置平面圖.Fig.10 The sketch of model with sole conductive anomalous body（a）The cross section of the model；（b）The plane view of the survey line.

圖11 坐標(biāo)系旋轉(zhuǎn)示意圖Fig.11 The sketch map of coordinate system rotation

張量測(cè)量時(shí)：

且有

對(duì)于張量的確定，可以把 “L”型源一個(gè)周期內(nèi)的上、下兩個(gè)半周期看作兩個(gè)不同的電偶極源，這樣（17）、（18）式變?yōu)?/p>

3.2 傾子測(cè)量

根據(jù)電磁場(chǎng)理論，垂直磁場(chǎng)分量可以表示為

其中Tzx、Tzy叫傾子，為復(fù)數(shù)：

設(shè)

在大地電磁中把Tzxrx＾＋Tzyry＾稱為帕金森矢量，具有指向高阻體，或者背離低阻體的指示作用，其長(zhǎng)度代表非均勻性強(qiáng)度.圖10所示地質(zhì)模型測(cè)區(qū)內(nèi)90Hz時(shí)帕金森矢量如圖13a所示，在異常體上方雖表現(xiàn)出一定異常，但整體方向指向電流源的方向，而喪失了指示異常體的作用.本文分析認(rèn)為這主要是受一次場(chǎng)源的影響所致，為了消除其影響，每個(gè)點(diǎn)上的值減去在沒(méi)有異常體時(shí)的值，可得圖13b，此時(shí)恢復(fù)了指示異常體的作用，在異常體上方帕金森矢量背離低阻體，并且在異常體邊緣長(zhǎng)度最長(zhǎng)，表明此處非均勻性最強(qiáng).在實(shí)際野外工作當(dāng)中，一般是很難知道背景場(chǎng)的電阻率的，此時(shí)可以選擇遠(yuǎn)參考點(diǎn)，假設(shè)此點(diǎn)不受異常體影響或受異常體影響很小，本文選擇測(cè)區(qū)（0，200）作為遠(yuǎn)參考點(diǎn)，每點(diǎn)都減去改點(diǎn)上的值可得圖13c，通過(guò)比較圖13b和圖13c可知，兩者相差不大，仍起到了很好地指示異常體的作用.

圖12 張量測(cè)量和標(biāo)量視電阻率對(duì)比圖Fig.12 Comparison curves of apparent resistivity between tensor measurements and scalar measurements

圖13 測(cè)區(qū)內(nèi)帕金森矢量分布圖（a）受電偶極源影響下的帕金森矢量；（b）減去一次場(chǎng)后的帕金森矢量；（c）減去遠(yuǎn)參考點(diǎn)后的帕金森矢量.Fig.13 The distribution map of Parkinson′s vector（a）The distribution map of Parkinson＇s vector under the effect of the electrical dipole source；（b）The distribution map of Parkinson′s vector eliminating the effect of the primary field；（c）The distribution map of Parkinson＇s vector eliminating the value of the far－referential site.

3.3 二維性指數(shù)

根據(jù)電磁理論，二維指數(shù)

具有判別地質(zhì)結(jié)構(gòu)維數(shù)的作用，當(dāng)S越小時(shí)說(shuō)明地下結(jié)構(gòu)越接近二維，愈大愈接近三維，一般當(dāng)S小于0.2～0.4時(shí)可判為二維（陳樂(lè)壽等，1989）.圖10所示地質(zhì)模型測(cè)區(qū)內(nèi)頻率取90Hz時(shí)的二維判別指數(shù)如圖14所示，可以看到在異常體上方呈現(xiàn)三維性，其他地方基本可視為二維結(jié)構(gòu)甚至一維結(jié)構(gòu).當(dāng)不知道地質(zhì)體結(jié)構(gòu)維數(shù)時(shí)，如果直接用異常體上方的測(cè)量數(shù)據(jù)進(jìn)行一維或者二維反演，必然引起很大誤差.其他地方的數(shù)據(jù)由于受異常體的影響微小、甚至可以忽略不計(jì)，所以可以直接進(jìn)行一維或者二維的反演，而不會(huì)引起很大的誤差.

圖14 二維性指數(shù)分布圖Fig.14 The distribution map of complex skew

4 結(jié)論

本文建立了一種全的“L”型源布極方式，在此基礎(chǔ)上推導(dǎo)了各個(gè)分量的表達(dá)式.

在公式推導(dǎo)的基礎(chǔ)上，為了充分利用偶極子的矢量疊加原理，設(shè)計(jì)了一種新的“L”型源信號(hào)發(fā)射模式.在新的信號(hào)發(fā)射模式下，360°張角范圍內(nèi)各分量不再存在弱區(qū)，輕松實(shí)現(xiàn)全張角范圍內(nèi)測(cè)量的目的；同時(shí)可以減少矢量測(cè)量或張量測(cè)量時(shí)源的布設(shè)次數(shù)，在節(jié)省時(shí)間的同時(shí)又可提高信噪比.在崎嶇山區(qū)，要滿足測(cè)量的苛刻條件有時(shí)是非常困難的，當(dāng)受地形條件限制無(wú)法找到合適的布極線路時(shí)，可根據(jù)矢量疊加原理選擇多偶極子源合成所需方向，以達(dá)到測(cè)量的目的.

在張量測(cè)量的事例中，計(jì)算了模型的旋轉(zhuǎn)角、傾子、二維判別指數(shù)等，通過(guò)旋轉(zhuǎn)可以得到構(gòu)造主軸上的視電阻率曲線，傾子對(duì)異常體有很好的指示作用，二維判別指數(shù)可以對(duì)地質(zhì)結(jié)構(gòu)的維數(shù)做出很好的判斷，從而選擇合適的反演方法.總之，張量測(cè)量可以提供豐富的信息，模擬的結(jié)果表明其明顯優(yōu)于標(biāo)量測(cè)量.

Boerner D E，Wright J A，Thurlow J G，et al.1993.Tensor CSAMT studies at the Buchans Mine in central Newfoundland.Geophysics，58（1）：12－19.

Bromley C.1993.Tensor CSAMT study of the fault zone between Waikite and Te Kopia geothermal fields.JournalofGeomagnetism andGeoelectricity，45（9）：887－896.

Chen L S，Liu R，Wang T S.1989.The Magnetotelluric Sounding Data Processing and Interpretration Methods （in Chinese）.Beijing：Petroleum Industry.

Chen L S，Wang G E.1990.Magnetotelluric Sounding Method（in Chinese）.Beijing：Geological Publishing House.

Di Q Y，Shi K F，Wang M Y，et al.2001.Water resources exploration with CSAMT and high density electric resistivity Method.ProgressinGeophysics.（in Chinese），16（3）：53－57.

Di Q Y，Wang M Y，Shi K F，et al.2002a.An applied study on prevention of water bursting disaster in mines with the high resolution V6system.ChineseJ.Geophys.（in Chinese），45（5）：744－748.

Di Q Y，Wang M Y，Shi K F，et al.2002b.Mutifunction system—V6and its effectiveness in CSAMT exploration.Progressin Geophysics.（in Chinese），17（4）：663－670.

Di Q Y，Wang R.2008.CSAMT Forward Modeling and Inversion and Its Application（in Chinese）.Beijing：Science Press.

Garcia X，Boerner D，Pedersen L B.2003.Electric and magnetic galvanic distortion decomposition of tensor CSAMT data.Application to data from the Buchans Mine （Newfoundland，Canada）.GeophysicalJournalInternational，154（3）：957－969.

He J S.1990.Control Source Audio Frequency Magnetotelluric（in Chinese）.Changsha：Central South University of Technology Press.

Homman G W.1975.Three－dimensional induced polarization and electromagnetic modeling.Geophysics，40（2）：309－324.

Liu J X，Guo Z W，Guo R W，et al.2009.Application of controlled source audiomagnetotelluric and gravity methods to surveys in the Lionlake hot spring geothermal area.ProgressinGeophys.（in Chinese），24（5）：1868－1873

Li D Q，Di Q Y，Wang G J，et al.2008.Fault detection by CSAMT and its application to new district planning in Beijing.Progress inGeophysics（in Chinese），23（6）：1963－1969.

Li X，Oskooi B，Pedersen L B.2000.Inversion of controlled－source tensor magnetotelluric data for a layered earth with azimuthal anisotropy.Geophysics，65（2）：452－464.

Lu X，Unsworth M，Booker J.1997.Two dimensional inversion of tensor CSAMT data／／1997SEG Annual Meeting.

Piao H R.1990.The Principle of Electromagnetic Sounding （in Chinese）.Beijing：Geological Press.

Raiche A P.1974.An integral equation approach to 3Dmodeling.GeophysicalJournalInternational，36（2）：363－376.

Shi K F.1999.Theory and Application of CSAMT Method （in Chinese）.Beijing：Science Press.

Sun Y X.2005.Application of CSAMT electromagnetic image system in survey of deeply－buried long tunnels for expressways.ProgressinGeophysics（in Chinese），20（4）：1184－1189.

Tan R J，Hu R L，Xu W J，et al.2007.CSAMT exploration and geological interpretation of perdue tectonic structures of Longpan deformation slope in Jinsha River.ProgressinGeophysics（in Chinese），22（1）：283－288.

Ting S C，Hohmann G W.1981.Integral equation modeling of three－dimensional magne－totelluric response.Geophysics，46（2）：182－197.

Tripp A C，Hohmann G W.1984.Block diagonalization of the electromagnetic impedance matrix of a symmetric buried body using group theory.IEEETransactionsonGeoscienceand RemoteSensing，22：62－69

Wannamaker P E.1997a.Tensor CSAMT survey over the Sulphur Springs thermal area，Valles Caldera，New Mexico，United States of America，Part I：Implications for structure of the western caldera.Geophysics，62（2）：451－465.

Wannamaker P E.1997b.Tensor CSAMT survey over the Sulphur Springs thermal area，Valles Caldera，New Mexico，United States of America，Part II：Implications for CSAMT methodology.Geophysics，62（2）：466－476.

Weidelt P.1975.Electomagnetic induction in three－dimensional structures.JournalofGeophysics，41：85－109.

附中文參考文獻(xiàn)

陳樂(lè)壽，劉任，王天生.1989.大地電磁測(cè)深資料處理與解釋.石油工業(yè)出版社.

陳樂(lè)壽，王光鍔.1990.大地電磁測(cè)深法.北京：地質(zhì)出版社.

底青云，石昆法，王妙月等.2001.CSAMT法和高密度電法探測(cè)地下水資源.地球物理學(xué)進(jìn)展，16（3）：53－57.

底青云，王妙月，石昆法等.2002a.高分辨率V6系統(tǒng)在礦山頂板涌水隱患中的應(yīng)用研究.地球物理學(xué)報(bào)，45（5）：744－748.

底青云，王妙月，石昆法等.2002b.V6多功能系統(tǒng)及其在CSAMT勘查應(yīng)用中的效果.地球物理學(xué)進(jìn)展，17（4）：663－670.

底青云，王若.2008.可控源音頻大地電磁數(shù)據(jù)正反演及方法應(yīng)用.北京：科學(xué)出版社.

何繼善.1990.可控源音頻大地電磁法.長(zhǎng)沙：中南工業(yè)大學(xué)出版社.

柳建新，郭振威，郭榮文等.2009.CSAMT和重力方法在獅子湖溫泉深部地球物理勘查中的應(yīng)用.地球物理學(xué)進(jìn)展，24（5）：1868－1873.

李帝銓，底青云，王光杰等.2008.CSAMT探測(cè)斷層在北京新區(qū)規(guī)劃中的應(yīng)用.地球物理學(xué)進(jìn)展，23（6）：1963－1969.

樸化榮.1990.電磁測(cè)深法原理.北京：地質(zhì)出版社.

石昆法.1999.可控源音頻大地電磁法理論與應(yīng)用.北京：科學(xué)出版社.

孫英勛.2005.CSAMT法在高速公路長(zhǎng)大深埋隧道勘察中的應(yīng)用研究.地球物理學(xué)進(jìn)展，20（4）：1184－1189.

譚儒蛟，胡瑞林，徐文杰等.2007.金沙江龍?bào)醋冃误w隱伏構(gòu)造CSAMT探測(cè)與解譯.地球物理學(xué)進(jìn)展，22（1）：283－288.