地球物理学报  2018, Vol. 61 Issue (1): 1-8   PDF    
短偏移距电磁探测技术概述
何继善1, 薛国强2,3     
1. 中南大学地球科学与信息物理学院, 长沙 410083;
2. 中国科学院矿产资源研究重点实验室, 中国科学院地质与地球物理研究所, 北京 100029;
3. 中国科学院大学, 北京 100049
摘要:由于人工源电磁探测效果与收发距离有一定关系,本文首先分析了收发距离为零的中心回线源瞬变电磁法(Transient Electromagnetic Method,TEM)视电阻率定义和算法特点,回溯了从圆形回线到矩形的修正式中心回线装置的发展过程;概述了从中心回线发展到接地导线源短偏移瞬变电磁(Short-Offset TEM,SOTEM)理论上的可行性;列举了将LOTEM(Long-Offset TEM)从数倍于探测深度的观测区推进到近源区后的探测实例.对于频率域人工源电磁法,介绍了起源于MT(Magneto-Telluric)的CSAMT(Controlled Source Audio-frequency Magneto-Telluric)与电磁频率测深法不同发展道路.电磁频率测深将观测区从远区推进到了中区;广域电磁法通过全区视电阻率新定义不仅提高了解释精度,而且将传统CSAMT远场观测模式推进到中区探测模式;最后指出,如能采用新的点微元假设计算方法,可有望分离自有场和辐射场,实现频域电磁近源探测.
关键词: 近源电磁探测      中心回线TEM      短偏移SOTEM      电磁频率测深      广域电磁法     
Review of the key techniques on short-offset electromagnetic detection
HE JiShan1, XUE GuoQiang2,3     
1. School of Geosciences and Info-Physics of Central South University, Changsha 410083, China;
2. Key Laboratory of Mineral Resources, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;
3. University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: Because the detecting effect of artificial-source electromagnetic methods is related to transceiver distance, taking the analysis expression and apparent resistivity algorithm as starting point, this paper reviews the development process from central-loop to modified central-loop TEM (Transient Electromagnetic Method), and summaries the principle of near-source prospecting from central-loop to grounded-wire source TEM. And the case history of SOTEM (Short-Offset TEM), which expands observation area of LOTEM (Long-Offset TEM) from several times of detection depth to near-source zone, is enumerated. For frequency-domain artificial-source electromagnetic methods, we present CSAMT (Controlled Source Audio-frequency Magnetotelluric) developed from MT (Magnetotelluric) and show the different way of frequency-domain electromagnetic sounding. Frequency sounding method enlarges the observation area from far-to middle-zone while the wide-field electromagnetic method improves interpretation accuracy and realize middle-zone data collection through whole zone apparent resistivity. If self-field and radiation field could be separated from each other based on the point-charge infinitesimal, frequency-domain near-source prospecting would be possible in the future.
Key words: Near-source electromagnetic prospecting    Central-loop TEM    Short-Offset TEM    Frequency-domain sounding    Wide field electromagnetic method    
0 引言

时间域和频率域人工源电磁法,能够可控地改变发射场源的时间序列或频率,使向地下发射的电磁场由浅入深传播,达到探测不同深度地层的目的(Kaufman and Keller, 1983Nabighian,1988).由于较高的信噪比、探测精度和工作效率,人工源电磁法在石油天然气、矿产资源、环境与工程等领域的资源勘查和矿山安全保障方面得到了广泛的应用.与场源来自高空的天然源MT方法不同,由于人工源的加入,观测点与场源的距离(发-收距)成为人工源电磁法勘探中需要考虑的重要因素,甚至是能否获得成功探测的关键.如何处理发-收距问题,不仅是对电磁理论深刻理解的结果(Kaufman and Keller, 1983),也是通过将观测点向近源区推进,减小体积效应、增大信噪比、提高数据观测质量以及提高探测精度的一个重要途径.

电磁法的近源勘探,最初由时间域的中心回线TEM装置(蒋邦远,1998牛之琏,2007)实现,然后根据时间域自有场(稳定电荷和电流场在场源附近占优,随源存在和消失)和辐射场(电磁场间的交互感应,场离开场源向外辐射)在时间上的分离性,实现了接地导线源的短偏移距SOTEM勘探(薛国强等,2013).在频率域的CSAMT测量中,如果能用某些方法去掉自有场,也可望有新的近源探测方法出现.由于解析解是电磁法的理论基础、视电阻率是和不同装置相联系的重要参量,定量反演是在正演基础上的程序化过程,因此本文以正演解析解和视电阻率为出发点,针对回线源和接地导线源的方法原理、观测仪器、应用实例等内容,分别概述了时域和频率的近源探测的发展过程和现状,展望进一步的工作.

1 时域电磁法回线装置的近源探测

时间域可控源电磁探测方法的基本原理是:位于地表的发射回线中的双极性阶跃脉冲(图 1)激励,在脉冲关断的间隙观测地下涡流产生的二次场.随着观测时间增加,感应场通过扩散过程穿透到更大深度(Kaufman and Keller, 1983Nabighian,1988),从接收线圈检测到的衰减曲线获得地下不同深度电导率变化信息.

图 1 双极性矩形脉冲激励源 Fig. 1 Bipolar rectangular pulse excitation source

对于典型的共中心回线源瞬变电磁装置,由于采用供电关断后的观测方式,可从时间上将自有场和辐射场分离了开来,实现了频域率中无法实现的同点测量(图 2a, 2b, 2c).为便于求得正演解析公式,将图 2所示方形回线视为圆形回线,圆柱坐标系下,接收线圈中的感生电动势为

图 2 时域TEM回线发射的近源装置 (a)共回线;(b)重叠回线;(c)中心回线;(d)修正中心回线. Fig. 2 Near-source confuguration of time-domain loop-source TEM (a) Common loop; (b) Coincident loop; (c) Central loop; (d) Modified central loop.

(1)

式中,L-1表示Laplace逆变换,t为观测时间,ω为角频率,J1为一阶Bessel函数,a为发射回线半径.此外有:

(2)

其中I为发射电流,Rn为地层总反射系数(石显新,2005).由于实际工程中方形回线铺设方便,发射圆回线半径a与方形回线边长L之间的关系为

(3)

式(1)中的单Bessel函数变为双重Bessel函数可展开为

(4)

式中的bn是接收线圈的半径和匝数.

根据式(4)导出的视电阻率公式为

(5)

式中STSR分别为发射和接收回线的面积.扩散深度为(Spies,1989)

(6)

进一步地,通常认为发射回线中部磁场较为均匀,为了提高实际工作效率,故将发射回线中部1/3作为观测区,如图 2d所示修正中心回线装置.

由于发射回线中心磁场均匀的假设,引起了修正中心回线的边缘效应,当中心点以外的观测点仍然沿用了中心回线视电阻率公式(5),由此产生的解释误差在深部陷落柱、导水小断层等精细勘探中不可忽视.为此,将修正中心回线与大定源回线(图 3)的视电阻率公式统一起来(石显新等,2009),即对于矩形回线内任意测点(图 2d),视电阻率公式为

图 3 大定源回线TEM装置 Fig. 3 Large fixed loop TEM configuration

(7)

式中,u2=Z(a0+a1y+…+a9y9)是与电阻率有关的函数(蒋邦远,1998),LxLy分别是发射回线短边和长边的边长,通常修正中心回线的边长Lx=Ly.

可用于地面回线源瞬变电磁勘探的仪器国外有加拿大Phoenix Geophysics公司的V-8多功能电磁法仪,美国Zonge公司的GDP-32多功能电法工作站,加拿大Geonics公司的PROTEM瞬变电磁系统,澳大利亚EMIT公司的SMARTEM-V;国内有重庆奔腾的WTEM,中国地质大学(武汉)的CUZGTEM,长沙白云的MSD-1,还有西安强源物探研究所的大功率单脉冲EMRS-3等.

水平分层大地上的回线源激发方式,使得地中感应电流仅在水平面中流动、不与地层界面相交,对探测低阻薄层非常有利(Kaufman and Keller, 1983),其中的修正中心回线TEM由于具有极高的工作效率,已经成为煤田水文地质勘探的首选方法.这方面的应用实例,广泛地见于《煤田地质与勘探》、《煤炭学报》、《Applied Geophysics》、《Leading Edge》、《Journal of Environmental & Engineering Geophysics》等刊物.

2 接地导线源短偏移距瞬变电磁法

水平层状大地上接地导线源装置激发的6个场分量,对低阻和高阻地层均有较好的敏感度,且有较大的探测深度.图 4是时变接地导线源除Ez外其余5个分量的场分布(闫述等,2017).与时域回线源不同的是,传统的时域接地导线源TEM继承了频域接地导线源CSAMT的观测方式:在发-收距大于4~6倍探测深度范围内观测,以便像时谐场那样通过距离分离自有场和辐射场,获得时变探测能力.这种接收点远离场源的观测方式称为LOTEM,通常情况下,LOTEM观测ExHz两个场量(Strack et al., 1990严良俊等,1999).众所周知,随收发距离增大场的强度下降,长偏移距在信噪比方面抵消了接地导线源带来的益处.为此,薛国强等(2013)依据时域中心回线TEM的探测原理,选用适当的发射波形,如矩形脉冲激励源关断后(图 1)的观测方式,分离自有场获得近源探测能力,提出短偏移距瞬变电磁法测量方法,并定名为SOTEM,将接地导线源瞬变电磁法的测区由远区扩大到中近区.

图 4 时变接地导线源的场分布(阴影部分为观测区) (a) ExHyHz分量; (b) EyHx分量. Fig. 4 Field distribution of time-varying grounded wire source(stippled areas indicate observation zones) (a) Ex, Hy and Hz components; (b) Ey and Hx components.

当观测点从远离源的区域推进到近源处以后,偶极子假设引起的误差不容忽视(薛国强等,2011; Zhou et al., 2013, 2015b),因此出现了将偶极子叠加求得有限长度接地导线源近似解的方法(李建平等,2007张成范等,2009; Zhou and Xue, 2014).进一步地,为了获得全场区接地导线源的精确解,将时变载流点电荷作为基本微元沿源积分,获得全场区有限长度接地导线源精确解.根据薛国强等(2014)提供的基本解,SOTEM垂直磁场Hz随时间的变化率公式为

(8)

(9)

由式(8)和式(9)所表示的全场区全时期视电阻率参量需用迭代方法求出:在确定的装置参数和时间t情况下,将估算视电阻率值ρ代入式(8)和式(9)后,得到感生电动势ε(t),与实测值ε(t)比较,如果两者间误差达到规定标准,停止迭代,否则迭代继续进行,直至满足精度要求.

SOTEM所用仪器与LOTEM相同,如德国Metrontx公司的DEMS-IV系统,发射机输出电压600 V,最大输出电流200 A;加拿大Phoenix Geophysics公司的V-8、美国Zonge公司的GDP-32等都有LOTEM功能,其中V-8的T-200发射机最大输出电压1600 V,最大输出电流160 A,GDP-32配置的GGT30发射机最大输出电压1000 V,最大输出电流45 A;配TEM67发射机的加拿大Geonics公司PROTEM等.

薛国强等在山西、山东、河南、安徽等地金属矿、煤矿、盐矿采用SOTEM进行了大深度探测(Chen et al., 2015; Zhou et al., 2015a, 2016a, b, 2017).如山东某煤矿的煤系地层顶部和底板含水区圈定,使用仪器为加拿大凤凰公司V-8综合电法仪,发射源AB=1000 m,MN=50 m,收发距r=1000~1500 m,探测深度1500 m,发送基频2.5 Hz,发送电流I=20 A,探头接收面积40000 m2.图 5是第三系底界面以上50 m顺层视电阻率切片图.圈出的低阻区(蓝色区域)与实际资料吻合.

图 5 第三系底界面以上50m顺层视电阻率切片 Fig. 5 Bed-parallel slices of apparent resistivity 50 m above Tertiary base
3 电磁频率测深与广域电磁法 3.1 电磁频率测深

频率域的人工源电磁法主要包括Cagniard视电阻率(Cagniard,1953)定义形式的CSAMT法(何继善,1990)和单分量视电阻率定义形式的频率测深法(陈明生和闫述,1995李毓茂,2012).其中,CSAMT有标量、矢量和张量三种装置形式.最常用的是标量装置,观测ExHy分量(图 4a);矢量装置观测ExEyHxHy分量(图 4b);张量装置是采用两个正交接地导线源的矢量装置.Cagniard电阻率公式为

(10)

(11)

以上公式仅在远区场条件下成立(杨儒贵,2008).从MT发展起来的CSAMT虽然继承了上述比值电阻率定义,但在勘探工程中,发射机功率、接收机灵敏度和日益严重的电磁干扰,都对大收发距形成了限制,使得远区场的条件不易满足.

电磁频率测深是相对独立发展的.1960年代初,前苏联出现了第一代频率测深仪器和相应的野外工作方法.我国于1967年开始研究电磁频率测深方法,1971年研制出第一台电磁频率测深仪器PC-1型,1983—1994年甘肃煤普查队研制出了ZP-1型瞬变电磁仪器,湖南物探队等单位也研制了电磁频率测深仪器(李毓茂,2012).由于受当时仪器的限制,大多只观测Ex分量,采用全区场公式为

(12)

式中各参数见闫述等(2017).

由于发射电流已知,故可设式(12)中地层参数Gk=Gξ=1并取远区渐进后,有单分量视电阻率为

(13)

实测中,式(13)转为

(14)

式中ABMN分别表示发射导线和接收导线的长度,这是野外工作中的习惯表示.

趋肤深度计算(Spies,1989)公式为

(15)

当进入近区场后,由式(14)计算出的视电阻率是真实电阻率的一半,且电场的假极值依然存在.因此,出现了许多改善视电阻率算法的研究工作(殷长春和朴化荣,1991毛先进和鲍光淑,1996).不过在另一方面,假极值却使得曲线拟合反演有较好的效果,当地层可视为近似水平时,如煤系建造勘探中,将公式(12)代入反演程序,如广义逆矩阵反演(陈明生和闫述,1995),可以取得不错的解释效果.

电磁频率测深以获得较强信号为主要目的,野外施工中当实测曲线进入近区渐近线2~3个频点后结束观测,然后进行曲线反演拟合获得地层结构,最大限度地把观测区从远区扩大到了仍然具备频测能力的中区.为了尽量减小偶极子假设引起的误差(何继善,1990),电磁频率测深采用了赤道偶极装置,即式(12)和(13)中θ=90°,只在发射回线的中垂线上观测,每观测一个测点后,发射和接收同时移动观测下一个测点,工作效率不高.

薛国强等在山西阳曲地区应用中煤科工集团西安研究院研制的DP-1数字频率测深仪进行频率测深勘探, 用于寻找含煤构造.工作中只测电场振幅曲线Ex分量一个参数, 发收距r=2000 m、AB=500 m、MN=200 m,赤道偶极装置,发射电流10 A,观测电场Ex分量.推断解释的煤田储量为10多亿吨,钻探资料证明该结论正确(薛国强和孟令杰,2004).

3.2 广域电磁法

进入2000年之后,随着加拿大PhoenixGeophysics公司的V8、美国Zonge公司的GDP32等仪器占据市场主流地位,国产频率测深仪器不再多用.同时,广域电磁法的优势得到普遍认可,观测范围不再限制于4~6倍探测深度的距离,测点可向近源方移动,像电磁频率测深那样可以充分利用中区勘探能力,进一步提高了观测速度、精度和野外效率(何继善,2010).

广域电磁法视电阻率表达式从均匀大地全场区公式中导出:即令公式(12)中的Gk=Gξ=1,有:

(16)

整理后,得:

(17)

与视电阻率公式(13)不同的是,式(17)右边的分母是ikr的函数,其中仍然含有未知的电阻率ρ.为了将视电阻率提取出来,广域电磁法提出的方法是(何继善,2010余云春,2010袁博,2013):ρ采用计算机迭代的办法.首先取一个可能的电阻率ρ,同时代入发送电流I、导线源的长度l、收发距r、方位角θ和工作频率f到式(17)中.根据得到的Ex与实测到的Ex相差多少.反复修改选取的ρ值,逐次迭代,直到得到的Ex与实测到的Ex符合精度要求为止.并且把最后选取的ρ值作为该装置、该工作条件下大地的视电阻率的最佳值.为了减少迭代次数,尽快地找到最佳值,可以采用数学上的优选法.

r < 6l,即观测点到发射导线中心的距离小于6倍的导线源长度时,视电阻率需沿发射导线积分, 公式为

(18)

式(17)和(18)所定义的视电阻率也称为广域视电阻率,当采用接地导线源观测Ex分量时,称为E-Ex方式的广域电磁法.若式(17)和(18)中的频率为0,则化为直流视电阻率.

显然,如加拿大Phoenix Geophysics公司的V-8、美国Zonge公司的GDP-32等同样适用于广域电磁法勘探.除此以外,湖南继善高科技有限公司自行研发有GY-1型广域电磁仪,最新研制的升级版为JSGY-2广域电磁仪.自2008年以来,随着广域电磁法研究不断深入(余云春,2010袁博,2013),该方法在青海柴达木盆地、大庆徐家围子、内蒙古大杨树盆地、内蒙古查干、新疆吐哈盆地、山西朔州、安徽淮南、山东招远金矿、安徽铜陵冬瓜山铜矿、内蒙古额尔古纳、湖南洞庭湖盆地、河南栾川、重庆城口和贵州六盘水等地开展了大量野外实验与研究工作.在常规油气勘探、深部金属矿勘探、地热调查和煤田采空区勘探等领域都取得了很好的应用效果(王宏宇,2013).

4 展望

在人工源电磁法发展的过程中,从时域回线源实现了频域中无法实现的近源(同点)探测开始,通过选择适当的发射波形和观测时机,将导线源的观测点从远区推进到包括中区和近区在内的全部场区,从LOTEM发展到SOTEM(从SOTEM的全场区观测特性来看,更确切的名称应该为WOTEM,Wide-Offset TEM).WOTEM在现有仪器装备条件下, 为深部资源探测提供了新的技术手段,可望在资源勘查、灾害地质、工程地质、水文地质等方面发挥更大的作用.

在频域电磁探测方面,广域电磁法和所定义的广义视电阻率将观测点从远区扩大到中区,提高了工作效率和解释精度.如果设法分离频域自有场与辐射场,将观测点推向近区不仅可以进一步增大勘探深度、提高信噪比,对人工源电磁场本身的发展也具有相当的价值.相信在国内外同行的共同努力下,短偏移电磁探测将会得到进一步的发展.

参考文献
Cagniard L. 1953. Basic theory of the magneto-telluric method of geophysical prospecting. Geophysics, 18(3): 605-635. DOI:10.1190/1.1437915
Chen M S, Yan S. 1995. Discussion Some Question in the Application of Frequency-Domain Sounding. Beijing: Geological Publishing House.
Chen W Y, Xue G Q, Muhammad Y K, et al. 2015. Application of short-offset TEM (SOTEM) technique in mapping water-enriched zones of coal stratum, an example from East China. Pure and Applied Geophysics, 172(6): 1643-1651. DOI:10.1007/s00024-014-1028-z
He J S. 1990. Controlled Source Audio-Frequency Magnetotelluric Sounding CSAMT. Changsha: Central South University of Technology Press.
He J S. 2010. Wide Electromagnetic Method and Pseudo-Random Signal Electric Method. Beijing: Higher Education Press.
Jiang B Y. 1998. Applied Near Zone Magnetic Source Transient Electromagnetic Exploration. Beijing: Geological Publishing House.
Kaufman A A, Keller G V. 1983. Frequency and Transient Soundings. New York:Elsevier.
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. Progress in Geophysics, 22(6): 1777-1780. DOI:10.3969/j.issn.1004-2903.2007.06.015
Li Y M. 2012. Frequency Sounding Method and the Chart of Frequency Sounding Generated By An Electric Dipole Source. Xuzhou: China University of Mining and Technology Press.
Mao X J, Bao G S. 1996. A direct algorithm for full-wave apparent resistivity from Horizontal electric dipole frequency soundings. Journal of Central South University of Technology, 27(3): 253-256.
Nabighian M N. 1988. Electromagnetic Methods in Applied Geophysics:Vol 2:Application. Tulsa:SEG Books.
Niu Z L. 2007. Theory of Time-Domain Electromagnetic Methods. Changsha: Central South University Press.
Shi X X. 2005. Research on shielding of low-resistivity layer in prospecting by transient electromagnetic method[Ph. D. thesis] (in Chinese). Beijing:China Coal Research Institute CCRI.
Shi X X, Yan S, Fu J M, et al. 2009. Improvement for interpretation of central loop transient electromagnetic method. Chinese Journal of Geophysics, 52(7): 1931-1936. DOI:10.3969/j.issn.0001-5733.2009.07.029
Spies B R. 1989. Depth of investigation in electromagnetic sounding methods. Geophysics, 54(7): 872-888. DOI:10.1190/1.1442716
Strack K M, Lüschen E, Kötz A W. 1990. Long-offset transient electromagnetic (LOTEM) depth soundings applied to crustal studies in the black Forest and Swabian Alb, Federal Republic of Germany. Geophysics, 55(7): 834-842. DOI:10.1190/1.1442897
Wang H Y. 2013. A research of Wide Field electromagnetic method applied in screening oil enrichment area[Master's thesis] (in Chinese). Changsha:Central South University.
Xue G Q, Chen W Y, Zhou N N, et al. 2013. Short-offset TEM technique with a grounded wire source for deep sounding. Chinese Journal of Geophysics, 56(1): 255-261. DOI:10.6038/cjg20130126
Xue G Q, Meng L J. 2004. Application of electromagnetic frequency Sounding method to Yangqu area in Shanxi province. Petroleum Instruments, 18(2): 14-16.
Xue G Q, Wang H Y, Yan S, et al. 2014. Time-domain Green function solution for transient electromagnetic field. Chinese Journal of Geophysics, 57(2): 671-678. DOI:10.6038/cjg20140230
Xue G Q, Yan S, Zhou N N. 2011. Theoretical study on the errors caused by dipole hypothesis of large-loop TEM response. Chinese Journal of Geophysics, 54(9): 2389-2396. DOI:10.3969/j.issn.0001-5733.2011.09.022
Yan L J, Hu W B, Chen Q L, et al. 1999. The estimation and fast inversion of all-time apparent resistivities in long-offset transient electromagnetic sounding. Oil Geophysical Prospecting, 34(5): 532-538.
Yan S, Xue G Q, Qiu W Z, et al. 2017. Interpretation of CSAMT single-component data. Chinese Journal of Geophysics, 60(1): 349-359. DOI:10.6038/cjg20170129
Yang R G. 2008. Advanced Electromagnetic Theory. Beijing: Higher Education Press.
Yin C C, Piao H R. 1991. A study of the definition of apparent resistivity in electromagnetic sounding. Geophysical & Geochemical Exploration, 15(4): 290-299.
Yu Y C. 2010. One-dimensional forward and inversion of Wide field electromagnetic method[Master's thesis] (in Chinese). Changsha:Central South University.
Yuan B. 2013. The numerical calculation of electric field r-direction in wide field electromagnetic method[Master's thesis] (in Chinese). Changsha:Central South University.
Zhang C F, Weng A H, Sun S D, et al. 2009. Computation of whole-time apparent resistivity of large rectangular loop. Journal of Jilin University (Earth Science Edition), 39(4): 755-758.
Zhou N N, Xue G Q. 2014. The ratio apparent resistivity definition of rectangular-loop TEM. Journal of Applied Geophysics, 103: 152-160. DOI:10.1016/j.jappgeo.2014.01.015
Zhou N N, Xue G Q, Chen W Y, et al. 2015a. Large-depth hydrogeological detection in the North China-type coalfield through short-offset grounded-wire TEM. Environmental Earth Sciences, 74(3): 2393-2404. DOI:10.1007/s12665-015-4240-y
Zhou N N, Xue G Q, Chen W Y, et al. 2017. A comparison of TEM data from different near-source systems. Journal of Geophysics and Engineering, 14(3): 487-501. DOI:10.1088/1742-2140/aa5cde
Zhou N N, Xue G Q, Gelius L J, et al. 2015b. Analysis of the near-source error in TEM due to the dipole hypothesis. Journal of Applied Geophysics, 115: 75-83.
Zhou N N, Xue G Q, Hou D Y, et al. 2016a. Short-offset grounded-wire TEM method for efficient detection of mined-out areas in vegetation-covered mountainous coalfields. Exploration Geophysics, 48(4): 374-382.
Zhou N N, Xue G Q, Li H, et al. 2016b. A comparison of different-mode fields generated from grounded-wire source based on the 1D model. Pure and Applied Geophysics, 173(2): 591-606. DOI:10.1007/s00024-015-1088-8
Zhou N N, Xue G Q, Wang H Y. 2013. Comparison of the time-domain electromagnetic field from an infinitesimal point charge and dipole source. Applied Geophysics, 10(3): 349-356. DOI:10.1007/s11770-013-0387-z
陈明生, 阎述. 1995. 论频率测深应用中的几个问题. 北京: 地质出版社.
何继善. 1990. 可控源音频大地电磁测深CSAMT. 长沙: 中南工业大学出版社.
何继善. 2010. 广域电磁法和伪随机信号电法. 北京: 高等教育出版社.
蒋邦远. 1998. 实用近区磁源瞬变电磁法勘探. 北京: 地质出版社.
李建平, 李桐林, 赵雪峰, 等. 2007. 层状介质任意形状回线源瞬变电磁全区视电阻率的研究. 地球物理学进展, 22(6): 1777–1780. DOI:10.3969/j.issn.1004-2903.2007.06.015
李毓茂. 2012. 电磁频率测深方法与电偶源电磁频率测深量板. 徐州: 中国矿业大学出版社.
毛先进, 鲍光淑. 1996. 水平电偶源频率域电磁测深全区视电阻率的直接算法. 中南工业大学学报, 27(3): 253–256.
牛之琏. 2007. 时间域电磁法原理. 长沙: 中南大学出版社.
石显新. 2005. 瞬变电磁法勘探中的低阻层屏蔽问题研究[博士论文]. 北京: 煤炭科学研究总院. http://cdmd.cnki.com.cn/Article/CDMD-83301-2006158772.htm
石显新, 闫述, 傅君眉, 等. 2009. 瞬变电磁法中心回线装置资料解释方法的改进. 地球物理学报, 52(7): 1931–1936. DOI:10.3969/j.issn.0001-5733.2009.07.029
王宏宇. 2013. 广域电磁法探测含油富集区试验研究[硕士论文]. 长沙: 中南大学. http://cdmd.cnki.com.cn/Article/CDMD-10533-1014145101.htm
薛国强, 陈卫营, 周楠楠, 等. 2013. 接地源瞬变电磁短偏移深部探测技术. 地球物理学报, 56(1): 255–261. DOI:10.6038/cjg20130126
薛国强, 孟令杰. 2004. 电磁频率测深法在山西阳曲找煤中的应用. 石油仪器, 18(2): 14–16.
薛国强, 王贺元, 闫述, 等. 2014. 瞬变电磁场时域格林函数解. 地球物理学报, 57(2): 671–678. DOI:10.6038/cjg20140230
薛国强, 闫述, 周楠楠. 2011. 偶极子假设引起的大回线源瞬变电磁响应偏差分析. 地球物理学报, 54(9): 2389–2396. DOI:10.3969/j.issn.0001-5733.2011.09.022
严良俊, 胡文宝, 陈清礼, 等. 1999. 长偏移距瞬变电磁测深的全区视电阻率求取及快速反演方法. 石油地球物理勘探, 34(5): 532–538.
闫述, 薛国强, 邱卫忠, 等. 2017. CSAMT单分量数据解释方法. 地球物理学报, 60(1): 349–359. DOI:10.6038/cjg20170129
杨儒贵. 2008. 高等电磁理论. 北京: 高等教育出版社.
殷长春, 朴化荣. 1991. 电磁测深法视电阻率定义问题的研究. 物探与化探, 15(4): 290–299.
余云春. 2010. 广域电磁法一维正反演[硕士论文]. 长沙: 中南大学.
袁博. 2013. 广域电磁法E-Er数值模拟[硕士论文]. 长沙: 中南大学.
张成范, 翁爱华, 孙世栋, 等. 2009. 计算矩形大定源回线瞬变电磁测深全区视电阻率. 吉林大学学报(地球科学版), 39(4): 755–758.