磁性層狀介質(zhì)的CSAMT電場(chǎng)分量響應(yīng)特征

閆述， 薛國(guó)強(qiáng)， 陳明生

1 江蘇大學(xué)計(jì)算機(jī)科學(xué)與通信工程學(xué)院， 江蘇鎮(zhèn)江 212013 2 中國(guó)科學(xué)院地質(zhì)與地球物理研究所, 北京 100029 3 中煤科工集團(tuán)西安研究院有限公司， 西安 710054

?

磁性層狀介質(zhì)的CSAMT電場(chǎng)分量響應(yīng)特征

閆述1， 薛國(guó)強(qiáng)2*， 陳明生3

1 江蘇大學(xué)計(jì)算機(jī)科學(xué)與通信工程學(xué)院， 江蘇鎮(zhèn)江 212013 2 中國(guó)科學(xué)院地質(zhì)與地球物理研究所, 北京 100029 3 中煤科工集團(tuán)西安研究院有限公司， 西安 710054

在傳統(tǒng)的勘探電磁學(xué)理論中，往往把巖層磁導(dǎo)率近似看成空氣中的磁導(dǎo)率，從而得到簡(jiǎn)化后的電磁場(chǎng)響應(yīng)式，主要用于指導(dǎo)非磁性介質(zhì)情況下的電磁探測(cè).本文開(kāi)展了磁性介質(zhì)的電磁響應(yīng)理論、數(shù)值模擬和特征分析等研究，以便實(shí)現(xiàn)對(duì)磁性礦體的電磁精細(xì)探測(cè).首先推導(dǎo)出磁性層情況下的電磁響應(yīng)表達(dá)式；然后，計(jì)算了幾種典型磁性介質(zhì)地質(zhì)斷面的響應(yīng)曲線，并對(duì)其響應(yīng)特征進(jìn)行了分析；研究結(jié)果表明：當(dāng)磁性層位于地表時(shí)，當(dāng)表層相對(duì)磁導(dǎo)率μr>1.2，磁性層對(duì)電磁測(cè)深曲線具體明顯的影響；當(dāng)巖層為中間層時(shí)，只要μr>1.2，對(duì)H型、A型曲線中段會(huì)有影響.對(duì)H型曲線，主要表現(xiàn)在使極小值處的曲線抬高、變寬，以致使曲線形態(tài)發(fā)生變化；當(dāng)磁性層處于最后一層時(shí)，不論巖層的磁導(dǎo)率是多少，對(duì)各類(lèi)二、三層曲線來(lái)講，不會(huì)影響曲線的形態(tài)，只是使?jié)u近線前第一個(gè)極值點(diǎn)的位置右移，極值點(diǎn)至漸近線之間的線段變得較為陡峭.

電磁學(xué)； 磁導(dǎo)率； 勘探； 響應(yīng)； 曲線

1 引言

勘探電磁法是通過(guò)觀測(cè)不同類(lèi)型的電磁信號(hào)以獲取地下目標(biāo)體的電性信息以達(dá)到探測(cè)的目的，傳統(tǒng)的電磁法主要考慮了巖石多種物理特性中的一種，即導(dǎo)電性.目前廣泛應(yīng)用的電磁場(chǎng)理論、勘探技術(shù)及數(shù)據(jù)解釋手段也大部分是在該基礎(chǔ)上形成的.然而，隨著電磁法應(yīng)用領(lǐng)域的越來(lái)越廣，人們對(duì)勘探精度要求的越來(lái)越高，僅僅基于巖石電阻率特性的數(shù)據(jù)解釋在許多情況下已不能滿(mǎn)足特定的勘探需求.這就要求考慮多方面的巖石物理特性，以更加全面地解釋觀測(cè)信號(hào)，提取更加豐富的巖石物性信息.

磁導(dǎo)率作為巖石的另一種重要物性參數(shù)，在通常的理論研究和實(shí)際應(yīng)用中一般不予考慮不同巖石礦物的不同，而將其設(shè)定為與空氣磁導(dǎo)率一樣的常數(shù).這樣的處理方式對(duì)于那些具有順磁性或逆磁性的巖石礦物來(lái)說(shuō)差別微乎其微；但是對(duì)于那些具有鐵磁性或超順磁性的礦物，當(dāng)外加磁場(chǎng)與之作用時(shí)，將會(huì)產(chǎn)生較強(qiáng)的附加磁場(chǎng)，并被接收裝置記錄，若同樣不考慮磁導(dǎo)率帶來(lái)的這些變化，將會(huì)嚴(yán)重影響后期數(shù)據(jù)處理解釋的精度(Ward and Hohamnn，1991).事實(shí)上，相當(dāng)一部分的巖礦石具有較高的相對(duì)磁導(dǎo)率，如磁鐵礦和基性巖發(fā)育區(qū)的礦層和巖層的相對(duì)磁導(dǎo)率往往大于1.

對(duì)于巖礦石磁性對(duì)電磁法帶來(lái)的影響，很早就引起了人們的注意,特別是磁性物質(zhì)對(duì)瞬變電磁響應(yīng)的影響的研究較為深入.Buselli(1982)研究了近地表超順磁物質(zhì)對(duì)瞬變電磁測(cè)深數(shù)據(jù)的影響；Lee(1984a,1984b)研究了超順磁性大地的瞬變電磁響應(yīng)特性；牛之璉(1985)引用Cole-Cole模型研究了瞬變電磁張弛效應(yīng)，指出瞬變電磁張弛效應(yīng)的衰減速度要比導(dǎo)電體瞬變渦流效應(yīng)慢的多，在晚期有可能測(cè)出瞬變電磁張弛效應(yīng)，磁張弛效應(yīng)作為一種附加效應(yīng)，可能成為限制電磁法探測(cè)深度的地質(zhì)噪聲源，但是在有利的地質(zhì)條件下也可以利用磁張弛效應(yīng)評(píng)價(jià)復(fù)雜磁異常；Kozhevnikov和Antonov(2008,2009,2011)分別對(duì)均勻大地、兩層大地及多層大地的磁張弛效應(yīng)對(duì)瞬變電磁數(shù)據(jù)的影響進(jìn)行了研究；殷偉偉(2013)針對(duì)稀疏浸染狀鐵礦的瞬變電磁響應(yīng)特征進(jìn)行了分析，并指出了針對(duì)此類(lèi)鐵礦石的瞬變電磁探測(cè)方法.對(duì)于磁性介質(zhì)對(duì)頻域響應(yīng)的影響，國(guó)內(nèi)外的研究主要集中于頻域航空電磁法中(Huang and Fraser, 2002; 高亮等，2009).習(xí)建軍(2010)和王衛(wèi)平等(2013)都針對(duì)基于磁導(dǎo)率和電阻率的頻率域航空電磁法的正反演進(jìn)行了研究.其他的研究還包括：閆述和陳明生(1997)從全場(chǎng)區(qū)出發(fā)，研究了頻率域測(cè)深中磁性層對(duì)水平電場(chǎng)響應(yīng)及視電阻率帶來(lái)的影響，指出磁性層的存在不可忽視，會(huì)對(duì)視電阻率曲線帶來(lái)較大的影響.

可控源音頻大地電磁法(CSAMT)的波區(qū)理論相對(duì)比較成熟，近年來(lái)，CSAMT單分量觀測(cè)模式的發(fā)展，使得CSAMT的觀測(cè)從單一的遠(yuǎn)區(qū)平面波區(qū)拓展至近區(qū)(Xue et al., 2015).為了獲取更加精確的解釋成果，解決更加復(fù)雜的地質(zhì)勘探任務(wù)，特別是實(shí)現(xiàn)大深度磁性鐵礦的精細(xì)探測(cè)，非常有必要研究巖礦石磁性特征給CSAMT觀測(cè)數(shù)據(jù)及數(shù)據(jù)解釋帶來(lái)的影響.本文首先推導(dǎo)了考慮磁性層的全場(chǎng)區(qū)水平電場(chǎng)頻率域表達(dá)式以及對(duì)應(yīng)的視電阻率表達(dá)式，然后研究了不同相對(duì)磁導(dǎo)率情況下的視電阻率變化特點(diǎn).本文的研究對(duì)于識(shí)別野外實(shí)測(cè)CSAMT數(shù)據(jù)中的磁異常效應(yīng)，提高數(shù)據(jù)解釋精度具有一定的理論和實(shí)際意義.

2 水平磁性層電磁響應(yīng)理論

設(shè)空氣的磁導(dǎo)率為μ0，分層大地各層的磁導(dǎo)率分別為μ1,μ2,…,μn.在距地面高度為h的(0,0,h)點(diǎn)放置一個(gè)諧變水平電偶極子Iax(見(jiàn)圖1)，x表示單位矢量.

圖1 水平磁性分層介質(zhì)Fig.1 Diagram of layered earth with magnetic permeability

在各層介質(zhì)的分界面上，例如第層i與第i+1層分界面上引入邊界條件：電場(chǎng)強(qiáng)度的切向分量總是連續(xù)的；當(dāng)分界面上無(wú)自由面電流時(shí)磁場(chǎng)強(qiáng)度的切向分量也是連續(xù)的，即

(1)

引入矢量位A求解E和H(張秋光，1983)，公式為

(2)

式中，ω為角頻率，A為矢量位，E代表電場(chǎng)矢量，H代表磁場(chǎng)矢量，下標(biāo)代表不同方向的分量，k為波數(shù)，k2=ω2μ ε+iω μ σ,其中，μ、σ和ε分別表示每一層的磁導(dǎo)率、電導(dǎo)率和介電常數(shù).

根據(jù)(1)式，矢量位A的邊界條件方程可寫(xiě)為

(3)

考慮到邊界條件以及Axi,Axi+1,Azi,Azi+1均滿(mǎn)足齊次亥姆赫茲方程，由此可確定矢量位的解為

(4)

利用(3)可求出(4)式中的4n個(gè)待定系數(shù)，再根據(jù)公式(2)及Ay=0可導(dǎo)出Ex為

(5)

由此確定了個(gè)場(chǎng)強(qiáng)分量的具體表達(dá)式.當(dāng)h=0時(shí)水平分層大地表面上電場(chǎng)強(qiáng)度Ex分量為

(6)

圖2 θ角示意圖Fig.2 Sketch of angle θ

得到水平電場(chǎng)Ex后，就可以求取Ex定義的視電阻率為

(7)

3 含磁性巖層的視電阻率響應(yīng)曲線

圖3 G型斷面視電阻率曲線(a) 第一層磁導(dǎo)率不同； (b) 第二層磁導(dǎo)率不同.Fig.3 Apparent resistivity curves of G-type model(a) Different permeability of the first layer； (b) Different permeability of the second layer.

圖4 D型斷面視電阻率曲線(a) 第一層磁導(dǎo)率不同； (b) 第二層磁導(dǎo)率不同.Fig.4 Apparent resistivity curves of D-type model(a) Different permeability of the first layer； (b) Different permeability of the second layer.

圖5 H型斷面視電阻率曲線(a) 第一層磁導(dǎo)率不同； (b) 第二層磁導(dǎo)率不同； (c) 第三層磁導(dǎo)率不同.Fig.5 Apparent resistivity curves of H-type model(a) Different permeability of the first layer； (b) Different permeability of the second layer； (c) Different permeability of the third layer.

從以上各圖可以看出，當(dāng)最下面一層巖層為磁性層時(shí)，不論是二層曲線還是三層曲線，對(duì)應(yīng)視電阻率曲線的影響都沒(méi)有第一層或中間層具有磁性時(shí)那么明顯了，即隨著基底磁導(dǎo)率的增大，視電阻率的增大幅度較小，而且隨著頻率的降低逐漸與μr3=1(基底無(wú)磁性)時(shí)視電阻率曲線尾部漸近線相重合.

4 含磁性巖層的視電阻率響應(yīng)特征

磁性層的視電阻率曲線響應(yīng)特征基本上在二層G、D和三層H型曲線上表示了出來(lái).當(dāng)某一層的相對(duì)磁導(dǎo)率增大時(shí)，對(duì)應(yīng)的視電阻率值也隨之增大.

(1) 磁性層位于頂層時(shí)

當(dāng)?shù)谝粚訋r層為磁性層時(shí)，對(duì)應(yīng)的視電阻率曲線的首部反映最為明顯.其增大的幅度為第一層巖層真電阻率的μr倍.具體分析如下：

電場(chǎng)強(qiáng)度Ex的公式(6)在遠(yuǎn)區(qū)場(chǎng)情況下，可簡(jiǎn)化為

(8)

而不考慮磁性層時(shí)的遠(yuǎn)區(qū)場(chǎng)公式為

(9)

(8)式與(9)式相比可知，因此視電阻率值是原來(lái)的μr1倍.

(2) 當(dāng)中間層為磁性層時(shí)，視電阻率曲線的變化比較明顯.對(duì)于三層曲線，當(dāng)相對(duì)磁導(dǎo)率達(dá)到一定程度，如大于4時(shí)，曲線形態(tài)都變得與K型曲線類(lèi)似.

(3) 磁性層位于底層時(shí)

當(dāng)?shù)撞繋r層為磁性層時(shí)，對(duì)應(yīng)的視電阻率曲線的反映不明顯.具體分析如下：

如果低頻時(shí)已進(jìn)入頻率測(cè)深近區(qū)場(chǎng)(k3r≤1)，則公式(6)變?yōu)?/p>

-F0(r,σ1,σ2,…,h1,h2,…)]，

(10)

其中，F(xiàn)n為直流測(cè)深中已知的函數(shù).在三層情況下，其表達(dá)式為

(11)

采用非線性數(shù)字濾波方法計(jì)算上式中的漢克爾變換(Ward and Hohmann, 1991).此時(shí)的近區(qū)場(chǎng)已變得與同裝置的直流電偶極子的場(chǎng)一樣，電場(chǎng)強(qiáng)度已與磁導(dǎo)率無(wú)關(guān).故最后一層巖層為磁性層時(shí)，其視電阻率曲線也總要與非磁性層的曲線尾部漸近線相重合.

可見(jiàn)，為了探測(cè)埋深較深的磁性巖層，應(yīng)該使相應(yīng)的視電阻率曲線段保持在遠(yuǎn)區(qū)或中區(qū)，以達(dá)到能探測(cè)磁性基底的效果.但是，在事先不知巖層參數(shù)的情況下，單純從曲線上仍然不易發(fā)現(xiàn)磁性層的存在，除了依賴(lài)地質(zhì)信息和物性信息外，還可以通過(guò)對(duì)比的方法，掌握磁性層響應(yīng)特性.下面以磁鐵礦為例進(jìn)行說(shuō)明.

圖6 鐵磁性礦層對(duì)H型曲線的影響特征(曲線上的數(shù)值表示第二層的μr2值)Fig.6 Effect of a ferromagnetic layer on apparent resistivity curves of H-type cross section(The μr2 of the second layer is represented by the numerical value of the curve)

設(shè)鐵磁性礦層位于模型的中間層，電阻率為0.001 Ωm，厚度為80 m，上覆地層電阻率為1 Ωm，厚度為10 m，基底電阻率為2 Ωm.圖6給出了中間磁層相對(duì)磁導(dǎo)率不同情況下的正演視電阻率曲線.從圖中可以看出，由于磁鐵礦是具有極低電阻率的礦層，當(dāng)μr=2時(shí)極小值的位置已有很大的移動(dòng)，達(dá)到8時(shí)則使曲線形態(tài)明顯改變了.此時(shí)如果將磁性層作為探測(cè)目標(biāo)時(shí)從上面對(duì)含磁性巖層分層大地頻率測(cè)深視電阻率計(jì)算結(jié)果所繪曲線看出，磁性巖層的存在有著明顯的視電阻率響應(yīng).可以此特征為解釋依據(jù)，將會(huì)獲得好的磁鐵礦探測(cè)效果.

5 結(jié)論

(2) 通過(guò)對(duì)磁性層響應(yīng)分析可知：磁性層埋藏深度越小，磁化特性表現(xiàn)越明顯，磁性層埋藏深度越大，磁化特性表現(xiàn)越不明顯；磁性層的電阻率越大，磁化特性表現(xiàn)越不明顯，磁性層的電阻率越小，磁化特性表現(xiàn)越明顯；遠(yuǎn)場(chǎng)情況下，磁化特性表現(xiàn)比較明顯，近場(chǎng)情況下，磁化特性基本沒(méi)有表現(xiàn).

(3) 具體地，(1)當(dāng)磁性層位于地表時(shí)，當(dāng)表層相對(duì)磁導(dǎo)率μr>1.2，磁性層對(duì)電磁測(cè)深曲線具體明顯的影響；(2)當(dāng)巖層為中間層時(shí)，只要μr>1.2，對(duì)H型、A型曲線中段會(huì)有影響.對(duì)H型曲線，主要表現(xiàn)在使極小值處的曲線抬高、變寬，以致使曲線形態(tài)發(fā)生變化；(3)當(dāng)磁性層處于最后一層時(shí)，不論巖層的磁導(dǎo)率是多少，對(duì)各類(lèi)二、三層曲線來(lái)講，不會(huì)影響曲線的形態(tài)，只是使?jié)u近線前第一個(gè)極值點(diǎn)的位置右移，極值點(diǎn)至漸近線之間的線段變得較為陡峭.此時(shí)磁性層的影響已趨于減弱了，特別對(duì)于μr<2的基性巖、超基性巖已無(wú)明顯影響.此時(shí)如將磁性巖層做為探測(cè)目標(biāo)，則是徒勞的.

(4) 在實(shí)際情況中，制約巖礦石磁導(dǎo)率的因素很多，對(duì)于不同的地質(zhì)環(huán)境，鐵磁巖礦石的磁導(dǎo)率在其平均值的附近有一個(gè)較大的變化范圍.本文計(jì)算時(shí)所采用的巖石物性參數(shù)值只具有一般性，但是計(jì)算結(jié)果不影響對(duì)磁性介質(zhì)響應(yīng)特性的分析.另外，利用磁性層對(duì)頻率測(cè)深視電阻率曲線的響應(yīng)，可探測(cè)近于水平的磁性巖層賦存；基性巖、超基性巖在地中的分布有許多是成水平層狀的，或者是可以看成是水平層狀的.鐵質(zhì)石英巖也可以在一個(gè)比較廣的區(qū)域成層狀分布.因此由水平層狀介質(zhì)計(jì)算出來(lái)的磁性層響應(yīng)特征曲線能很好地反映這類(lèi)磁性巖石的分布.對(duì)于陡產(chǎn)狀及不規(guī)則產(chǎn)狀的磁鐵礦床，盡管他們的相對(duì)磁導(dǎo)率非常大，μr=2.2～25，但不能套用水平情況下得出的結(jié)論，只是有助于定性方面的分析.

Buselli G. 1982. The effect of near-surface superparamagnetic material on electromagnetic measurements.Geophysics, 47(9): 1315-1324.

Gao L, Hu X Y, Wang W Y, et al. 2009. A study of frequency airborne electromagnetic forward under magnetic condition.ChineseJournalofEngineeringGeophysics(in Chinese), 6(4): 399-403.

Huang H P, Fraser D C. 2002. Dielectric permittivity and resistivity mapping using high-frequency, helicopter-borne EM data.Geophysics, 67(3): 727-738.

Kozhevnikov N O, Antonov E Y. 2008. The magnetic relaxation effect on TEM responses of a uniform earth.RussianGeologyandGeophysics, 49(3): 197-205.

Kozhevnikov N O, Antonov E Y. 2009. The magnetic relaxation effect on TEM responses of a two-layer earth.RussianGeologyandGeophysics, 50(10): 895-904.

Kozhevnikov N O, Antonov E Y. 2011. Magnetic relaxation of a horizontal layer: effect on TEM data.RussianGeologyandGeophysics, 52(4): 398-404.

Lee T J. 1984a. The effect of a superparamagnetic layer on the transient electromagnetic response of a ground.GeophysicalProspecting, 32(3): 480-496.

Lee T J. 1984b. The transient electromagnetic response of a magnetic or superparamagnetic ground.Geophysics, 49(7): 854-860.

Niu Z L. 1985. The effect of transient magnetic relaxation and its significance in evaluating the anomaly of complex magnetism.JournalCSIMM(in Chinese), (2): 79-85.

Wang W P, Zeng Z F, Wu C P, et al. 2013. Frequency domain airborne electromagnetic dual-frequency inversion method based on magnetic permeability.GeologicalScienceandTechnologyInformation(in Chinese), 32(2): 181-186.

Ward S H, Hohmann G W. 1991. Electromagnetic theory for geophysical exploration. ∥ Nabighian N ed. Electromagnetic Methods in Applied Geophysics. Tulsa, Oklahoma: Society of Exploration Geophysics, 121-223.

Xi J J. 2010. A study of forward and inversion for frequency-domain airborne electromagnetic method in terms of permeability and resistivity[Master thesis] (in Chinese). Changchun: Jilin University.

Xue G Q, Yan S, Gelius L J, et al. 2015. Discovery of a major coal deposit in china with the use of a modified CSAMT method.JournalofEnvironmental&EngineeringGeophysics, 20(1): 47-56.

Yan S, Chen M S. 1997. The responsibility of apparent resistivity curves on frequency sounding for magnetic layers.GeologyandProspecting(in Chinese), 33(3): 38-42.

Yin W W. 2013. Analysis on the response characteristic of the transient electromagnetic method to detect the sparse disseminated

iron ore[Master thesis] (in Chinese). Taiyuan: Taiyuan University of Technology. Zhang Q G. 1983. Field Theory (in Chinese). Beijing: Geology Press.

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

高亮, 胡祥云, 王衛(wèi)一等. 2009. 磁性條件下頻率域航空電磁法正演研究. 工程地球物理學(xué)報(bào), 6(4): 399-403.

牛之璉. 1985. 瞬變磁張弛效應(yīng)及其在評(píng)價(jià)復(fù)雜磁異常中的意義. 中南礦冶學(xué)院學(xué)報(bào), (2): 79-84.

王衛(wèi)平, 曾昭發(fā), 吳成平等. 2013. 基于磁導(dǎo)率的頻率域航空電磁法雙頻反演方法. 地質(zhì)科技情報(bào), 32(2): 181-186.

習(xí)建軍. 2010. 基于磁導(dǎo)率和電阻率的頻率域航空電磁法正反演研究[碩士學(xué)位論文]. 長(zhǎng)春: 吉林大學(xué).

閆述, 陳明生. 1997. 頻率測(cè)深中磁性層的視電阻率響應(yīng). 地質(zhì)與勘探, 33(3): 38-42.

殷偉偉. 2013. 稀疏浸染狀鐵礦瞬變電磁響應(yīng)特征分析[碩士學(xué)位論文]. 太原: 太原理工大學(xué).

張秋光. 1983. 場(chǎng)論. 北京: 地質(zhì)出版社.

(本文編輯 張正峰)

Response of CSAMT electromagnetic component in the magnetized layered earth

YAN Shu1， XUE Guo-Qiang2*， CHEN Ming-Sheng3

1TelecommunicationEngineering，JiangsuUniversity,JiangsuZhenjiang212013,China2InstituteofGeologyandGeophysics，ChineseAcademyofSciences，Beijing100029,China3CCTEMGxianresearchinstitute,Xi′an710054,China

In the traditional exploration electromagnetic theory, the rock permeability is usually regard as the air permeability approximately. As a result, the electromagnetic field response can be derived simply. However, study on electromagnetic response considering rock permeability is required for fine electromagnetic detection of magnetic ore bodies. In this paper, electromagnetic field response expressions of the magnetic layered earth is deduced at first, then the response curves of several kinds of typical magnetic medium geological sections are calculated. The research results show that electromagnetic sounding curves are heavily affected by the rock permeability layer when the earth is covered by a magnetic layer with rock permeabilityμr>1.2. While as H-type and A-type curves are heavily affected by rock permeability (μr>1.2) for the middle layer. On the other hand, when the magnetic layer is buried as the basement, electromagnetic sounding curves are slightly affected by the rock permeability layer.Keywords Electromagnetic research； Permeability； Exploration； Response； Curve

10.6038/cjg20161208.

國(guó)家重大科研裝備研制項(xiàng)目“深部資源探測(cè)核心裝備研發(fā)”(ZDYZ2012-1)-05子項(xiàng)目“多通道大功率電法勘探儀”-04課題“M-TEM資料處理及偏移成像軟件研制”,國(guó)家自然科學(xué)基金(41374129，41474095)聯(lián)合資助.作者簡(jiǎn)介 閆述，女，1953年生，教授，主要從事電磁場(chǎng)與微波技術(shù)研究.E-mail:yanshu@ujs.edu.cn

10.6038/cjg20161208

P631

2015-12-02,2016-10-22收修定稿

閆述， 薛國(guó)強(qiáng)， 陳明生. 2016. 磁性層狀介質(zhì)的CSAMT電場(chǎng)分量響應(yīng)特征. 地球物理學(xué)報(bào)，59(12):4457-4463,

Yan S， Xue G Q， Chen M S. 2016. Response of CSAMT electromagnetic component in the magnetized layered earth.ChineseJ.Geophys. (in Chinese),59(12):4457-4463,doi:10.6038/cjg20161208.

*通訊作者 薛國(guó)強(qiáng)，男，1966年生，研究員，主要從事電磁探測(cè)方法研究.E-mail:ppxueguoqiang@163.com