地球物理学进展  2014, Vol. 29 Issue (1): 418-422   PDF    
环形源形状对时间域电磁响应的影响分析
李耀恩1, 吴云鹏2, 吴小平3    
1. 内蒙古自治区第一水文地质工程地质勘查院, 呼和浩特 010020;
2. 成都信息工程学院计算机学院, 成都 610225;
3. 中南大学地球科学与信息物理学院, 长沙 410083
摘要:时间域电磁测深方法中, 重叠回线、中心回线和大定源等回线源装置, 因其非接地、近场源观测、小装置等优势被广泛应用.但野外观测中, 回线布设时由于野外地形条件的限制, 往往不能布置成规则形状, 回线源的形状对其时间域电磁响应的影响就自然成为一个非常重要因素.作者在文中通过对圆形源、正方形源和正三角形源等回线源的中心点的时间域电磁响应分别按等周长和等面积计算分析, 发现任意形状回线源中心点的时间域电磁响应与回线源的形状无关, 而是与回线源所包围的面积成正比.并由此推断, 任意形状回线源中心点时间域电磁响应观测数据的处理, 可以采取等面积圆形回线源的等效处理与解释.
关键词时间域     电磁测深     回线源     偶极源     电磁响应    
Analysis of the affection for the electromagnetic response of the shape variation of the loop source
LI Yao-en1, WU Yun-peng2, WU Xiao-ping3    
1. Inner Mongolia First Hydrogeology Engineering Geological Prospecting Institute, Hohehot 010020, China;

2. College of Computer Science & Technology of Chengdu University of Information Technology, Cheng du 610225, China;

3. School of Geosciences and Info-physics of Central South University, Changsha 410083, China

Abstract: The Overlap Loop, Central Loop and Fixed Loop are widely used for Time Domain Electromagnetic Sounding because of their non-grounding, near source and small installation superiority. It is very hard to lay a regulation Loop in the field condition because of restriction of field topography. So it is important for analysis of the effects by the variation shape of Loop source. The author found that the central point EM response of Loop source is independent of the shape variation of Loop source and direct ratio to the area of Loop source by the compare of calculating results of circle Loop source, Square Loop source and Triangle Loop source. And get conclusion: it can be processed and interpreted by using Circle Loop source equivalent processing for variation Loop source but the same area.
Key words: time domain     electromagnetic sounding     loop source     dipole source     electromagnetic response    

0 引 言

时间域电磁测深方法以其优越的探测深度和良好的探测效果,越来越得到广泛的应用(考夫曼和凯勒,1983;李建平等,2007滕吉文,2008王华军,2008吴小平等,2010吴小平和吴云鹏,2012).早期的电磁激励方式,主要采用电偶极子或者磁偶极子的激励方式,在远离激励源的远区或者时间域的晚期进行观测处理(考夫曼和凯勒.1983).随着人们对电磁场激励效果研究的深入,越来越多的采用激励效果更佳的回线源方式,如重叠回线、中心回线和大定源等(戚志鹏等,2011石显新等,2009薛国强等,2008).

回线源激励效果固然很好,但由于激励场的时间域响应非常复杂(吴小平和吴云鹏,2012),加上野外观测时回线布设受到野外实际条件的限制,往往实际回线布置并非规则形状,这样给时间域电磁响应的数据处理与解释增加很大困难(吴小平等,2010).有学者曾就任意回线源的电磁响应,采取对任意回线进行等分切割,然后采用回归分析处理的方式,采取正演加反演的联合处理方式来处理任意回线源的时间域观测结果(李建平等,2007).

本文作者通过对圆形回线源、正方形回线源和正三角形回线源的中心点时间域电磁响应的计算(吴小平等,2013)并进行比较分析,发现等周长的不同形状回线源的中心点时间域电磁响应并不一致,而等面积不同形状回线源的中心点时间域电磁响应完全一致.确认不同回线源的中心点的时间域电磁响应与回线源的形状无关,而是与回线源包围面积成正比.并由此推断任意形状回线源中心点时间域电磁响应观测结果可以采用等面积圆形回线源的等效处理与解释方式进行.

1 电偶源与磁偶源的电磁响应

图1所示,为水平电偶极子源和磁偶极子源在均匀半无限空间上方的示意图,当使用阶跃波在偶极子中激励时,位于地面P点的电偶极子源的垂直磁场的时间域电磁响应表达式为(考夫曼和凯勒,1983.)

图 1电偶源与磁偶源示意图 Fig.1The Electrical Dipole and Magnetic Dipole on the Uniform space

其中PE=Idl为电偶极矩;为概率函数;r为观测点P到电偶极子中心的距离;为自由空间的磁导率;ρ为均匀半空间电阻率.

同样,磁偶极子源的垂直磁场的时间域电磁响应表达式为(考夫曼和凯勒,1983.)

其中PM=IS为磁偶极矩;为概率函数;r为观测点P到磁偶极子中心的距离;为自由空间的磁导率;ρ为均匀半空间电阻率.

从(1)式可知,电偶极源的时间域电磁响应与偶极源中的电流强度和偶极源的长度是成正比例的,要想加大电偶极源的激励效果,就必须有一定的电流强度和偶极长度,但随着偶极长度的增加,观测点就要相应距偶极源更远,否则,偶极源就不具备电偶极子源的特性.然而,如果要加大电流强度,就要使发送设备非常庞大.

从(2)式可以看出,磁偶极子源的时间域电磁响应与偶极子源中的电流强度和偶极子的有效激发面积成正比例,同样,如果想要增大磁偶极子源的激励效果,就必须有一定的电流强度和磁偶极子的有效激发面积,但随着磁偶极子体积的增大,观测点距其的距离就要相应的增加,否则,磁偶极子源将不再具备磁偶极子特性.

随着对时间域近场源激励的研究,人们发现回线源可以充分发挥激励源的激励效果(滕吉文,2008王华军,2008吴小平等,2010吴小平和吴云鹏,2012).重叠回线、中心回线和大定源等激励源方式越来越多的被采用(戚志鹏等,2011石显新等,2009),但近场源的电磁响应非常复杂(李建慧等,2012),目前一般流行采用回线源的特例,即圆形的回线源的理论计算结果来近似其他形状的回线源的电磁响应.

2 环形源的均匀半无限空间上的时间域电磁响应

我们这里主要分析几种典型的回线源,图2所示的是,一个半径为r的圆形回线源、一个边长为s的正方形回线源和一个边长为a的正三角形回线源,我们都分析其中心点的电磁响应.

图 2几种典型环形源示意图 Fig.2Some Typical Loop Sources

对于如图2中的圆形回线源,如果我们任意取一段电偶极电流源,若计算其圆心点的电磁响应就相对简单,由(1)式中,电流方向与径向方向始终垂直,因此,圆形回线源的中心的时间域电磁响应为

其中I为环形源中的电流强度;为概率函数;r为环形源半径;为自由空间的磁导率;ρ为均匀半空间电阻率.

对于方形源和三角源,可以看成是由四边或三边的长线源对中心点电磁响应的叠加总和.可以利用长线源的计算方式(吴小平等,2013),先计算各边的长线源的电磁响应,然后根据电磁场的叠加性原理进行累加.

3 不同形状环形源的电磁响应比较分析

对不同形状环形源的电磁响应的分析,具有十分重要的实际应用价值,因为在实际野外应用中,所有回线源的回线布设都会受到野外的实际条件限制,不可能得到与理论计算假设一样的结果.这就会给实测结果资料的处理与解释带来困难.

图2所示的圆形、正方形和三角形等三种特殊情况,我们先假定他们周长相等,即:

图3所示为半径为100米的圆形源,放置于一个电阻率为100 Ωm的均匀半无限地面上, 线框中电流10 A,在圆形中心点的电磁响应,以及与其等周长的正方形和正三角形的中心点的时间域电磁响应.

图 3等周长不同形状环形源电磁响应 Fig.3EM Response Compare of Variation Loop Source with the same Perimeter

从图中可以看出,其时间域电磁响应曲线形态基本一致,但其响应值不一致,其中圆形电磁响应最强,正三角形响应最弱.

图2中的三种特殊回线源,我们假定其回线面积相等,即:

图4所示为半径为100 m的圆形源,放置于一个电阻率为100 Ωm的均匀半无限地面上,线框中电流10 A,在圆形中心点的电磁响应,以及与其等面积的正方形和正三角形的中心点的时间域电磁响应.
图 4等面积不同形状环形源电磁响应 Fig.4EM Response compare of variation loop source with the same area

从图中可以看出,其时间域电磁响应曲线完全一致.这说明回线源激励时,其中心点的时间域电磁响应与回线源的形状无关,而与回线源所围面积成正比.

由此可以推断,任意形状的回线源,在其中心的时间域电磁响应,可以等效于与其等面积的圆形回线源的时间域电磁响应.这样,对于野外时间域观测中,任意布设的回线源的中心点的电磁响应的数据处理与解释就可以转化为圆形源的数据处理与解释.

4 结 语

时间域回线源电磁测深工作的数据处理与解释工作,往往因野外工作中受野外地形条件的限制,回线源形状呈不规则布置,给资料处理和解释工作带来困难.本文中通过圆形回线源、正方形回线源和正三角形回线源的中心点时间域电磁响应的计算得知:回线源中心点的时间域电磁响应与回线源的形状关系不大,而是与回线源所围面积成正比的,因此,可以得出以下结论:

(1)对于野外任意形状回线源中心电磁响应结果,都可以根据其面积,选择其等面积的圆形回线源的数据处理与解释方式来进行等效处理与解释.

(2)野外回线源布设过程中,可以根据野外实际条件,任意布置回线源的形状,观测结果均可以很方便的等效解释.

参考文献
[1] Combrinck M. 2009. Analysis of numerical differentiation methods applied to time domain electromagnetic (TDEM) geophysical data in the S-layer differential transform[J]. Computer & Geosciences, 35(8): 1563-1573.
[2] Flores C, Peralta-Ortega S A. 2009. Induced polarization with in-loop transient electromagnetic soundings: A case study of mineral discrimination at El Arco porphyry copper, Mexico[J]. Journal of Applied Geophysics, 68(3): 423-436.
[3] Han Z H, Wei W B, Zhang W B. 2008. Study on the exploration depth of transient electromagnetic method in northern China coalfield[J]. Progress in Geophysics (in Chinese), 23(1): 237-241.
[4] Jia D Y, Weng A H, Liu Y H, et al. 2013. Propagation of electromagnetic fields from a horizontal electrical dipole buried in ocean[J]. Progress in Geophys. (in Chinese), 28(1): 507-514, doi: 10. 6038/pg20130518.
[5] Kaufman A A, Keller G A. 1983. Frequency and Transient Sounding [M]. New York: Elsevier Science Publishes.
[6] Lei D. 2010. Studies and applications of 2D CSAMT modeling and inversion with a dipole source and topography[J]. Chinese J. Geophys. (in Chinese), 53(4): 982-993, doi: 10. 3969/j. issn. 0001-5733. 2010. 04.023.
[7] Li J H, Zhu Z Q, Zeng S H, et al. 2012. Progress of forward computation in transient electromagnetic method[J]. Progress in Geophys. (in Chinese), 27(4): 1393-1400, doi: 10. 6038/j. issn. 1004-2903. 2012. 04.013.
[8] Li J P, Li T L, Zhao X F, et al. 2007. Study on the TEM all-time apparent resistivity of arbitrary shape loop source over the layered medium[J]. Progress in Geophysics (in Chinese), 22(6): 1777-1780.
[9] Li X Y, Sun D, Jing Y Q, et al. 2009. Scanning analysis of effective Q using generalized S transform[J]. Progress in Geophys. (in Chinese), 24(5): 1696-1702, doi: 10. 3969/j. issn. 1004-2903. 2009. 05.020.
[10] Lin J Y, Tang J T, Ding M B, et al. 2009. Application of object-oriented technology in geophysical electromagnetic visualization[J]. Progress in Geophys. (in Chinese), 24(6): 2233-2241, doi: 10. 3969/j. issn. 1004-2903. 2009. 06. 043.
[11] Qi Z P, Li X, Zhu H W, et al. 2011. Definition of apparent resistivity for non-center vertical component of Large-loop TEM configuration[J]. Progress in Geophys. (in Chinese), 26(4): 1350-1358, doi: 10. 3969/j. issn. 1004-2903. 2011. 04. 028.
[12] Shi X X, Yan S, Fu J M, et al. 2009. Improvement for interpretation of central loop transient electromagnetic method[J]. Chinese J. Geophys. (in Chinese), 52(7): 1931-1936, doi: 10. 3969/j. issn. 0001-5733. 2009. 07.029.
[13] Tartaras E, Zhdanov M S, Wada K, et al. 2000. Fast imaging of TDEM data based on S-inversion[J]. Journal of Applied Geophysics, 43(1): 15-32.
[14] Teng J W. 2008. The core scientific problems and development direction for the contemporary geophysical research[J]. Progress in Geophysics (in Chinese), 23(3): 637-640.
[15] Wang H J. 2008. Time domain transient electromagnetism all time apparent resistivity translation algorithm[J]. Chinese J. Geophys.(in Chinese), 51(6): 1936-1942.
[16] Wang M, Deng M, Zhang Q S, et al. 2009. The technique of time synchronization operation to control marine electromagnetic emission[J]. Progress in Geophys. (in Chinese), 24(4): 1493-1498, doi: 10. 3969/j. issn. 1004-2903. 2009. 04. 043.
[17] Weng A H, Lu D H, Liu G X. 2003. Definition of whole zone apparent resistivity for transient electromagnetic method of current dipole source [J]. Coal Geology & Exploration, 31(3): 56-59.
[18] Wu X P, He J S, Liu J X. 2006. The effect of the electromagnetic wave’s impedance transform to the results of the ground-probing-radar[J]. Progress in Geophysics (in Chinese), 21(1): 251-255.
[19] Wu X P, Liu J X, Duan W H. 2010a. A utilitarian method for the explanation of time domain electromagnetic sounding[J]. Progress in Geophys. (in Chinese), 25(3): 898-903, doi: 10. 3969/j. issn. 1004-2903. 2010. 03. 024.
[20] Wu X P, Liu J X, Duan W H, et al. 2010b. The effect of topography aberration on the pseudo-section of the dipole-dipole array resistivity sounding[J]. Progress in Geophys. (in Chinese), 25(2): 635-641, doi: 10. 3966/j. issn. 1004-2903. 2010. 02.037.
[21] Wu X P, Wu Y P. 2012a. A numerical calculating method for apparent resistivity of TEM sounding[J]. Progress in Geophys. (in Chinese), 27(6): 2548-2553, doi: 10. 6038/j. issn. 1004-2903. 2012. 06.031.
[22] Wu X P, Wu Y P. 2012b. A comparative analysis for the dipole source and the long line source electromagnetic response[J]. Progress in Geophys. (in Chinese), 27(6): 2602-2608, doi: 10. 6038/j. issn. 1004-2903. 2012. 06.038.
[23] Wu X P, Wu Y P, Liu H L. 2013. A dissected dipole calculation method for the electromagnetic response of long line source[J]. Progress in Geophysics (in Chinese), 28(5): 2667-2672, doi: 10. 6038/pg20130548.
[24] Xie W, Liu J X, Zhang D F. 2009. Development of magnetotelluric data processing software based on VC++[J]. Progress in Geophys. (in Chinese), 24(6): 2242-2249, doi: 10. 3969/j. issn. 1004-2903. 2009. 06. 044.
[25] Xing G L, Wang H J, Yang S D. 2008. The response functions of electromagnetic wave logs in the 2-D axis-symmetric formation[J]. Chinese J. Geophys. (in Chinese), 51(3): 924-932.
[26] Xue G Q, Li X. 2008. Correction of TEM Early-Time sounding date[J]. Progress in Geophysics (in Chinese), 23(1): 220-224.
[27] Xue G Q, Li X, Di Q Y. 2008. Research progress in TEM forward modeling and inversion calculation[J]. Progress in Geophysics (in Chinese), 23(4): 1165-1172.
[28] Yu J H, Yang H Y, Deng J Z. 2012. Whole-space effect on underground transient electromagnetic field in layered media[J]. Progress in Geophys. (in Chinese), 27(4): 1385-1392, doi: 10. 6038/j. issn. 1004-2903. 2012. 04.012.
[29] Zhdanov M S, Pavlov D A. 2001. Analysis and interpretation of anomalous conductivity and magnetic permeability effects in time domain electromagnetic data Part II: Sμ-inversion[J]. Journal of Applied Geophysics, 46(4): 235-248.
[1] 韩自豪, 魏文博, 张文波. 2008. 华北煤田瞬变电磁勘探深度研究[J].   地球物理学进展, 23(1): 237-241.
[2] 贾定宇, 翁爱华, 刘云鹤,等. 2013. 海洋环境中水平电偶极子电磁场特征分析[J].    地球物理学进展, 28(1): 507-514, doi: 10.6038/pg20130518.
[3] 考夫曼, 凯勒. 1987. 频率域和时间域电磁测深[M]. 北京: 地质出版社.
[4] 雷达. 2010. 起伏地形下CSAMT二维正反演研究与应用[J].    地球物理学报, 53(4): 982-993, doi: 10. 3969/j. issn. 0001-5733. 2010. 04.023.
[5] 李建慧, 朱自强, 曾思红,等. 2012. 瞬变电磁法正演计算进展[J].    地球物理学进展, 27(4): 1393-1400, doi: 10. 6038/j. issn. 1004-2903. 2012. 04.013.
[6] 李建平, 李桐林, 赵雪峰,等. 2007. 层状介质任意形状回线源瞬变电磁全区视电阻率的研究[J].   地球物理学进展, 22(6): 1777-1780.
[7] 李雪英, 孙丹, 井涌泉,等. 2009. 利用广义S变换进行等效Q值扫描分析[J].   地球物理学进展, 24(5): 1696-1702, doi: 10. 3969/j. issn. 1004-2903. 2009. 05. 020.
[8] 林家勇, 汤井田, 丁茂斌,等. 2009. 面向对象技术在地球物理电磁成像可视化中的应用[J].   地球物理学进展, 24(6): 2233-2241, doi: 10. 3969/j. issn. 1004-2903. 2009. 06. 043.
[9] 戚志鹏, 李 貅, 朱宏伟,等. 2011. 大定源装置下瞬变电磁法视电阻率定义[J].   地球物理学进展, 26(4): 1350-1358, doi: 10. 3969/ j. issn. 1004-2903. 2011. 04. 028.
[10] 石显新, 闫述, 傅君眉,等. 2009. 瞬变电磁法中心回线装置资料解释方法的改进[J].    地球物理学报, 52(7): 1931-1936, doi: 10. 3969/j. issn. 0001-5733. 2009. 07.029.
[11] 滕吉文. 2008. 当代地球物理学研究的核心科学问题和发展导向[J].   地球物理学进展, 23(3): 637-640.
[12] 王华军. 2008. 时间域瞬变电磁法全区视电阻率的平移算法[J].   地球物理学报, 51(6): 1936-1942.
[13] 王猛, 邓明, 张启升,等. 2009. 控制海底电磁激发脉冲发射的时间同步技术[J].    地球物理学进展, 24(4): 1493-1498, doi: 10. 3969/j. issn. 1004-2903. 2009. 04.043.
[14] 翁爱华, 陆冬华, 刘国兴. 2003. 利用连分式定义瞬变电磁法全区视电阻率研究[J].   煤田地质与勘探, 31(3): 56-59.
[15] 吴小平, 何继善, 柳建新. 2006. 电磁波的阻抗变换作用对地质雷达探测效果的影响[J].   地球物理学进展, 21(1): 251-255.
[16] 吴小平, 柳建新, 段无悔. 2010a. 一种实用的时间域电磁测深解释方法[J].   地球物理学进展, 25(3): 898-903, doi: 10. 3969/j. issn. 1004-2903. 2010. 03. 024.
[17] 吴小平, 柳建新, 段无悔,等. 2010b. 地形畸变对偶极-偶极电阻率测深拟断面图的影响[J].    地球物理学进展, 25(2): 635-641, doi: 10. 3969/j. issn. 1004-2903. 2010. 02.037.
[18] 吴小平, 吴云鹏. 2012a. 时间域电磁法视电阻率的一种数值计算方法[J].   地球物理学进展, 27(6): 2548-2553, doi: 10. 6038/j. issn. 1004-2903. 2012. 06. 031.
[19] 吴小平, 吴云鹏. 2012b. 非零长线源与偶极源的电磁响应比较分析[J].   地球物理学进展, 27(6): 2602-2608, doi: 10. 6038/j. issn. 1004-2903. 2012. 06. 038.
[20] 吴小平, 吴云鹏, 刘洪亮. 2013. 非零长线源的电磁响应的偶极剖分计算方法[J].    地球物理学进展, 28(5): 2667-2672, doi: 10.6038/pg20130548.
[21] 谢维, 柳建新, 张东风. 2009. 基于VC++大地电磁数据处理软件开发[J].   地球物理学进展, 24(6): 2242-2249, doi: 10. 3969/j. issn. 1004-2903. 2009. 06. 044.
[22] 邢光龙, 王宏建, 杨善德. 2008. 二维轴对称介质中电磁波测井的响应函数[J].   地球物理学报, 51(3): 924-932.
[23] 薛国强, 李貅. 2008. 瞬变电磁测深早期数据的修正[J].   地球物理学进展, 23(1): 220-224.
[24] 薛国强, 李貅, 底青云. 2008. 瞬变电磁法正反演问题研究进展[J].   地球物理学进展, 23(4): 1165-1172.
[25] 岳建华, 杨海燕, 邓居智. 2012. 层状介质中地下瞬变电磁场全空间效应[J].    地球物理学进展, 27(4): 1385-1392, doi: 10. 6038/j. issn. 1004-2903. 2012. 04.012.