E-mail Alert Rss
 

物探与化探, 2022, 46(2): 383-391 doi: 10.11720/wtyht.2022.1541

方法研究·信息处理·仪器研制

地—井瞬变电磁三维响应特征分析与异常体快速定位方法研究

赵友超,1, 张军,1,2, 范涛3, 姚伟华3, 杨洋1,4, 孙怀凤1,4

1.山东大学 岩土与结构工程研究中心,山东 济南 250061

2.山东省交通规划设计院,山东 济南 250031

3.中煤科工集团 西安研究院有限公司,陕西 西安 710077

4.山东省工业技术研究院 先进勘探与透明城市协同创新中心,山东 济南 250061

Analysis of 3D ground-borehole TEM response characteristics and rapid positioning method for anomalous bodies

ZHAO You-Chao,1, ZHANG Jun,1,2, FAN Tao3, YAO Wei-Hua3, YANG Yang1,4, SUN Huai-Feng1,4

1. Geotechnical and Structural Engineering Research Center, Shandong University, Jinan 250061, China

2. Shandong Provincial Communications Planning and Design Institute, Jinan 250031, China

3. Xi'an Research Institute of China Coal Technology & Engineering Group, Xi'an 710077, China

4. Advanced Exploration and Transparent City Innovation Center, Shandong Research Institute of Industrial Technology, Jinan 250061, China

通讯作者: 张军(1995-),男,主要从事隧道工程设计与研究、超前地质预报工作。Email:1035058515@qq.com

责任编辑: 沈效群

收稿日期: 2021-09-27   修回日期: 2021-12-12  

基金资助: 山东省自然科学基金项目(ZR2019MD20)
国家自然科学基金项目(42074145)
国家自然科学基金项目(42004056)

Received: 2021-09-27   Revised: 2021-12-12  

作者简介 About authors

赵友超(1997-),男,主要从事地—井瞬变电磁正反演与应用研究工作。Email: 201914579@mail.sdu.edu.cn

摘要

通过对含异常体的三维地电模型进行系统正演和分析,提出了一种地—井瞬变电磁异常体快速定位方法。在对含异常体的三维地电模型进行正演响应规律分析时发现:纯异常场的X、Y分量曲线零点以及Z分量曲线的极值点与异常体深度对应较好;当异常体的规模、电阻率、埋深发生改变时,X、Y分量曲线形态基本不变;当异常体方位发生变化时,X、Y分量曲线形态发生变化。在此基础上,提出地—井瞬变电磁异常体快速定位方法:首先,依据X、Y曲线的零点确定异常体所在深度;其次,依据X、Y分量曲线形态确定异常体所在象限(90°范围);最后,依据X+Y或者X-Y分量的曲线形态将异常体定位于井周45°范围内。数值试验显示本方法对不同方位异常体模型的定位结果均与设计模型吻合。采用陕北某矿区地—井瞬变电磁实测数据作进一步验证,本文方法推断的结果与实际吻合较好。

关键词： 地—井瞬变电磁; 时域有限差分; 三维正演; 异常体定位

Abstract

Through systematic forward modeling and analysis of a 3D geoelectric model containing anomalous bodies, this study proposed a rapid positioning method of anomalous bodies based on the ground-borehole transient electromagnetic (TEM) method. The analysis of the forward modeling response laws of the 3D geoelectric model containing anomalous bodies shows that the zero points of the X and Y component curves and the extreme points of the Z component curve of a pure anomaly field correspond well to the depths of the anomalous bodies; the morphologies of the X and Y component curves are basically unchanged when the sizes, resistivity, and burial depths of the anomalous bodies change but change when the orientations of anomalous bodies change. On this basis, this study proposed the following method to rapidly position anomalous bodies using the ground-borehole TEM method. First, determine the depths of anomalous bodies according to the zero points of the X and Y curves. Next, determine the quadrants (within 90°) of anomalous bodies according to the morphologies of the X and Y component curves. Finally, position anomalous bodies within 45° of boreholes according to the morphologies of the X+Y or X-Y component curves. Numerical experiments show that the positioning results of models of anomalous bodies with different orientations are consistent with those of the model designed in this study. As further verified using the ground-borehole TEM measured data of a mining area in northern Shaanxi, the inference that there is a water-filled goaf to the northwest of the borehole obtained using the method proposed in this study well agrees with the actual situation.

Keywords： ground-borehole transient electromagnetic; time-domain finite difference; 3D forward modeling; positioning of anomalous bodies

PDF (7700KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

赵友超, 张军, 范涛, 姚伟华, 杨洋, 孙怀凤. 地—井瞬变电磁三维响应特征分析与异常体快速定位方法研究[J]. 物探与化探, 2022, 46(2): 383-391 doi:10.11720/wtyht.2022.1541

ZHAO You-Chao, ZHANG Jun, FAN Tao, YAO Wei-Hua, YANG Yang, SUN Huai-Feng. Analysis of 3D ground-borehole TEM response characteristics and rapid positioning method for anomalous bodies[J]. Geophysical and Geochemical Exploration, 2022, 46(2): 383-391 doi:10.11720/wtyht.2022.1541

0 引言

瞬变电磁法(TEM)是一种利用电磁感应原理进行勘探的电磁测深方法,通过采集关断间歇产生的纯二次场来分析地下电性分布[1]。瞬变电磁法已被广泛应用于各类勘察工作中,包括矿产资源勘探[2⇓-4],查找地下水[5-6]以及隧道超前地质预报[7⇓⇓-10]等方面。地—井瞬变电磁法是在常规地面瞬变电磁基础上发展起来的地面发射、井中接收的瞬变电磁勘探方法[11]。相比传统地面瞬变电磁方法,地—井瞬变电磁勘探接收探头置于井下,其接收点更加接近目标体,观测到的响应信号会更强,信号反馈的目标体信息也更加真实可靠;同时也在一定程度上规避了地表电磁干扰[12]。

在地—井瞬变电磁正演方面,中外不少学者都进行过深入研究[13⇓⇓⇓-17],而在这一领域,地下异常体定位仍是亟待解决的问题。Dyck等[18]采用物理模型试验的方法并结合数值模拟研究了采用地—井瞬变电磁三分量数据对地下目标体进行定位的可行性,提出同时观测三分量数据有助于解决多解性问题,但并没有给出针对异常体定位的确切方法;J.R.Bishop等[19]将地—井瞬变电磁法用于地下磁黄铁矿探查,并取得了良好的应用效果;孟庆鑫等[20]采用时域有限差分研究了均匀半空间以及低阻覆盖层下均匀半空间中的低阻板状异常的响应规律;徐正玉等[21]采用时域有限差分方法研究了均匀半空间模型和垂直断层面模型的地—井瞬变电磁响应规律。目前,对于地—井瞬变电磁异常体定位的正演研究大多处于定性分析阶段,无法明确定位异常体所处范围。

在上述研究的基础上,针对地—井瞬变电磁正演问题的定性分析做了进一步完善。采用时域有限差分方法建立三维模型进行正演计算,研究地—井瞬变电磁响应特征,并在此基础上提出了一种针对井旁异常体的快速定位方法。本研究在孙怀凤等开发的TEM3DFDTD程序基础上完成。

1 地—井瞬变电磁三维正演与规律分析

基于时域有限差分方法,建立三维地电模型进行正演计算[22-23],分析模型的三分量响应曲线,总结其响应规律,模型如图1所示。采用表1所示模型参数,设异常体为正方体,依据单一变量原则分别改变异常体的电阻率、尺寸、埋深以及水平方位角,对比分析曲线形态特征并总结规律。水平方位角是指异常体中心与坐标原点的连线与X轴正向的夹角。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   模型示意

Fig.1   Schematic diagram of model


表1   三维模型正演参数

Table 1   3-D model forward parameters

模型1模型2模型3模型4
埋深/m505050/100/15050
ρ异常体/(Ω·m)10/50/100101010
水平方位角/(°)45454545/135/225/315
尺寸/m22/4/622
ρ围岩/(Ω·m)1 0001 0001 0001 000
异常体与钻孔
水平间距/m		10101010

新窗口打开| 下载CSV


为研究异常体电阻率的变化对地—井瞬变电磁三分量响应特征的影响,采用表1中的模型1进行正演计算,结果见图2。分析图2,发现X、Y分量响应曲线零点以及Z分量响应曲线极值点对应的深度即异常体所在深度。此外,随着异常体电阻率的增大,三分量的响应曲线形态均无明显变化,X、Y分量始终呈反“S”型,Z分量保持“V”型。但三分量响应曲线的幅值随异常体电阻率增加而逐渐减小,说明地—井瞬变电磁探测对低阻更为敏感,因此低阻模型的感应电动势响应更强烈。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   不同电阻率异常体模型的响应曲线

Fig.2   Response curves of different target resistivity models


为研究异常体尺寸变化对地—井瞬变电磁响应的影响,采用表1中的模型2进行正演计算。对比分析图3,随着异常体尺寸增加,三分量响应曲线的幅值逐渐增加。当异常体边长为2 m和4 m时,三分量响应曲线形态大致保持不变,仍然可以依据X、Y分量零点以及Z分量极值点定位异常体所在深度。但是当异常体边长为6 m时,X、Y分量曲线在异常体中心附近出现多个零点和极值点,Z分量曲线也在异常体中心附近出现多个极值点,推断该现象是由于异常体距离钻孔太近以及网格尺寸太小等因素引起。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   不同规模异常体模型的响应曲线

Fig.3   Response curves of different target size


为研究异常体埋深对三分量响应的影响,采用表1中模型3进行正演计算。分析图4发现,随着埋深增加,X、Y分量的曲线形态始终保持不变,且X、Y分量曲线的零点及Z分量曲线的极值点始终对应异常体所在深度。 随着异常体埋深增加,地—井瞬变电磁三分量响应逐渐减弱。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   不同目标体埋深模型响应曲线

Fig.4   Response curves of different depth of target


为研究异常体水平方位角的变化对三分量响应曲线的影响,采用表1中的模型4进行正演计算。图5显示,当异常体位于不同方位角时,X、Y分量曲线的零点以及Z分量曲线的极值保持不变,均对应异常体中心所在的深度;当异常体位于不同方位时,X、Y分量曲线形态会发生变化。当α=45°时,X、Y分量曲线形态均为反“S”型;当α=135°时,X分量呈“S”型而Y分量呈反“S”型;当α=225°时,X、Y分量均为“S”型;当α=225°时,X分量呈反“S”型而Y分量呈“S”型。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   不同水平方位模型的响应曲线

Fig.5   Response curves of different target azimuth


综合分析以上研究结果,发现如下规律:纯异常场X、Y分量曲线呈“S”型或者反“S”型,Z分量曲线呈“V”型,X、Y分量曲线的零点以及Z分量的极值点对应异常体所在的深度;由于距离钻孔太近以及网格尺寸太小等因素可能会造成响应曲线在异常体所在深度出现多个极值点;当异常体方位角发生改变时,X、Y分量曲线形态会发生变化,其他参数改变时,X、Y分量曲线形态保持不变。

2 异常体定位方法

在地—井瞬变电磁三维正演计算部分,发现只有当异常体方位发生改变时,X、Y分量曲线的形态才会发生改变。基于数值模拟部分得出的结论,不断改变异常体所在方位,进行大量正演计算,绘制其纯异常响应曲线并进行总结发现如下规律:

1) 当异常体位于第一象限时,X、Y分量曲线形态均为反“S”型;当异常体位于第二象限时,X分量呈“S”型而Y分量呈反“S”型;当异常体位于第三象限时,X、Y分量均为“S”型;当异常体位于第四象限时,X分量呈反“S”型而Y分量呈“S”型。

2) 在同一象限内,异常体位于不同区域时,X、Y分量响应曲线幅值的相对大小有所不同(这里的幅值是指响应曲线极大值点)。如图6所示,将井周360°范围划分成8个区域,不同区域内X、Y分量幅值相对大小如表2所示。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   区域划分示意

Fig.6   Schematic diagram of regional division


表2   X、Y分量幅值对比

Table 2  Comparison of the amplitude of X and Y components

区域幅值区域幅值
1|dBx/dt|>|dBy/dt|5|dBx/dt|>|dBy/dt|
2|dBx/dt|<|dBy/dt|6|dBx/dt|<|dBy/dt|
3|dBx/dt|<|dBy/dt|7|dBx/dt|<|dBy/dt|
4|dBx/dt|>|dBy/dt|8|dBx/dt|>|dBy/dt|

新窗口打开| 下载CSV


当异常体位于第一、第三象限时,X、Y分量曲线型态一致,可采用X-Y分量来进一步判断异常体所在区域;当异常体位于第二和第四象限时,X、Y分量曲线型态相反,可采用X-(-Y)即X+Y分量来进一步判断异常体所在区域。对上述研究进行归纳总结,提出一种地—井瞬变电磁异常体定位方法,依据该方法可确定异常体所在的深度并将其定位于井周45°范围。该方法对异常体的定位分两步进行:首先依据X、Y分量曲线的零点确定异常体深度;随后,依据X、Y分量曲线形态可以确定异常体所在象限,依据X-Y或者X+Y分量曲线形态可以进一步确定异常体所在的区域,如表3所示。

表3   目标体区域定位

Table 3  Regional positioning of target

纯异常曲线形态象限X-Y/X+Y曲线形态区域
X:反“S”型,Y:反“S”型第一X-Y:反“S”型1
X:反“S”型,Y:反“S”型第一X-Y:“S”型2
X:“S”型,Y:反“S”型第二X+Y:反“S”型3
X:“S”型,Y:反“S”型第二X+Y:“S”型4
X:“S”型,Y:“S”型第三X-Y:“S”型5
X:“S”型,Y:“S”型第三X-Y:反“S”型6
X:反“S”型,Y:“S”型第四X+Y:“S”型7
X:反“S”型,Y:“S”型第四X+Y:反“S”型8

新窗口打开| 下载CSV


3 数值算例验证

为了验证提出的地—井瞬变电磁异常体快速定位方法的准确性,选择同一个异常体,设置背景电阻率为1 000 Ω·m,异常体电阻率为10 Ω·m,其所在深度为50 m。改变异常体与X轴正向夹角,依次为30°、60°、120°、150°、210°、240°、300°以及330°,使得目标体分别位于划分好的8个区域并进行正演计算,如图7所示。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   异常体方位示意

Fig.7   Schematic diagram of target azimuth


图8(图8-1、图8-2)为8个模型的纯异常场分量响应曲线,依据每组模型响应曲线的零点均可确定异常体深度为50 m,与设计模型异常体所在深度相吻合。根据之前提出的水体定位方法,模型1到模型8异常体区域定位结果如图9所示。对比设计模型,发现所有模型定位结果均与异常体所在区域相吻合,数值试验证明了基于总场的地—井瞬变电磁异常体快速定位方法的准确性和可靠性。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   -1 纯异常分量模型响应曲线

Fig.8   -1 Pure abnormal response curves of 8 models


图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   -2 纯异常分量模型响应曲线

Fig.8   -2 Pure abnormal response curves of 8 models


图9

新窗口打开| 下载原图ZIP| 生成PPT

图9   异常体定位结果

Fig.9   The results of target positioning


4 实测数据验证

为验证提出的异常体快速定位方法的野外适用性,要对其进行现场试验验证。本次现场试验采用的实测数据为文献[24]中的数据,采用本文提出的方法对其进行处理。数据采集所用钻孔位于陕北某矿区,该矿区早期主要采取以掘代采的采煤方法,在地下留下较多的采空区,由于上覆含水层的补给,形成了较多的积水采空区,威胁煤矿生产安全。现已探明位于所选钻孔北偏西15°方向-97 m处存在充水采空区,如图10所示。

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10   采空区所在区域示意

Fig.10   Schematic diagram of the location of the goaf


回线源尺寸设置为360 m×360 m,钻孔位于回线源中间,发射波形为线性关断梯形波,关断时间为0.5 ms,发射电流为20 A。以正东方向为X轴正方向,正北方向为Y方向建立坐标系,采集得到的总场三分量响应曲线如图11所示,提取得到的异常场三分量响应曲线如图12所示。

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11   实测总场三分量响应曲线

Fig.11   Total field response curve of measured data


图12

新窗口打开| 下载原图ZIP| 生成PPT

图12   纯异常场响应曲线

Fig.12   Pure abnormal field response curve of the measured data


图12为纯异常场的三分量响应曲线,图中红色曲线大致描述了实测数据各分量的曲线形态。X、Y分量曲线的零点以及Z分量曲线的零点所在位置为地下97 m附近,据此可判定积水采空区深度为地下97 m附近。此外,X分量曲线大致呈反“S”型,Y分量曲线大致呈“S”型,因此可判定异常体位于第二象限;X+Y分量曲线大致呈“S”型,因此可进一步判定充水采空区位于第三区域,其定位结果如图13所示,与其实际位置相吻合。

图13

新窗口打开| 下载原图ZIP| 生成PPT

图13   充水采空区定位结果

Fig.13   The results of water-filled goaf positioning


5 结论

针对地—井瞬变电磁数据解释以及井旁异常体定位问题,基于时域有限差分方法建立三维模型进行大量正演计算,总结地—井瞬变电磁模型三分量响应规律,在此基础上提出一种井旁异常体快速定位方法。该方法主要依据三分量响应曲线形态以及组合曲线(X-Y或X+Y)形态来确定异常体所在深度并将异常体圈定于井周45°范围,对实际井旁矿体的探测具有指导作用。

地—井瞬变电磁响应特征规律较为复杂,本文提出的方法虽然能够快速圈定井旁异常体,但仍存在一定的局限性且有一定的进步空间。文中仅考虑了含单个异常体的均匀半空间模型的三分量响应规律,而对于背景围岩非均匀以及井旁存在多个异常体等问题仍需进一步深入研究;此外,本文提出的定位方法仅可确定异常体中心深度并将其圈定于井周45°范围,还没有实现异常体的精确定位。下一步工作是进行地—井瞬变电磁反演研究,到达精确定位井旁异常体的需求。

参考文献

朴化荣. 电磁测深法原理[M]. 北京: 地质出版社, 1990.

[本文引用: 1]

Piao H R. Principle of electromagnetic sounding[M]. Beijing: Geological Publishing House, 1990.

[本文引用: 1]

Irvine R J.

Drillhole TEM surveys at Thalanga, Queensland

[J]. Exploration Geophysics, 1987, 18: 285-293.

DOI:10.1071/EG987285      URL     [本文引用: 1]

Xue G Q, Qin K Z, Li X, et al.

Discovery of a Large-scale porphyry molybdenum deposit in tibet through a modified TEM exploration method

[J]. Journal of Environmental and Engineering Geophysics, 2012, 17(1): 19-25.

DOI:10.2113/JEEG17.1.19      URL     [本文引用: 1]

Yang D, Oldenburg D W.

Three-dimensional inversion of airborne time-domain electromagnetic data with applications to a porphyry deposit

[J]. Geophysics, 2012, 77(2): B23-B34.

[本文引用: 1]

Ezersky M.

TEM study of the geoelectrical structure and groundwater salinity of the Nahal Hever sinkhole site, Dead Sea shore, Israel

[J]. Journal of Applied Geophysics, 2011, 75(1): 99-112.

DOI:10.1016/j.jappgeo.2011.06.011      URL     [本文引用: 1]

Keydar S D.

Application of seismic diffraction imaging for detecting near-surface inhomogeneities in the Dead Sea area

[J]. Journal of Applied Geophysics, 2010, 71(2): 47-52.

DOI:10.1016/j.jappgeo.2010.04.001      URL     [本文引用: 1]

李术才, 刘斌, 孙怀凤, 等.

隧道施工超前地质预报研究现状及发展趋势

[J]. 岩石力学与工程学报, 2014, 33(6):1090-1113.

[本文引用: 1]

Li S C, Liu B, Sun H F, et al.

State of ART and trends of advanced geological prediction in tunnel construction

[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(6): 1090-1113.

[本文引用: 1]

薛国强, 李貅.

瞬变电磁隧道超前预报成像技术

[J]. 地球物理学报, 2008, 51(3):894-900.

[本文引用: 1]

Xue G Q, Li X.

The technology of TEM tunnel prediction imaging

[J]. Chinese Journal of Geophysics, 2008, 51(3): 894-900.

[本文引用: 1]

李术才, 孙怀凤, 李貅, 等.

隧道瞬变电磁超前预报平行磁场响应探测方法

[J]. 岩石力学与工程学报, 2014, 33(7):1309-1318.

[本文引用: 1]

Li S C, Sun H F, Li X, et al.

Advanced geology prediction with parallel transient electromagnetic detection in tunneling

[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(7): 1309-1318.

[本文引用: 1]

李貅, 薛国强, 李术才, 等. 瞬变电磁隧道超前预报方法与应用[M]. 北京: 地质出版社, 2013.

[本文引用: 1]

Li X, Xue G Q, Li S C, et al. The method and application of transient electromagnetic in tunnel prediction[M]. Beijing: Geological Publishing House, 2013.

[本文引用: 1]

Dyck A V.

Drill-hole electromagnetic methods

[M]// Electromagnetic Methods in Applied Geophysics. Tulsa: Society of Exploration Geophysicists, 1991.

[本文引用: 1]

牛之链. 时间域电磁法原理[M]. 长沙: 中南大学出版社, 2007.

[本文引用: 1]

Niu Z L. Principle of time domain electromagnetic[M]. Changsha: Central South University Press, 2007

[本文引用: 1]

Woods D V.

A model study of the Crone borehole pulse electromagnetic (PEM) system

[R]. Ontario: Queen’s University, 1975.

[本文引用: 1]

Eaton P A, Hohmann G W.

The influence of a conductive host on two-dimensional borehole transient electromagnetic responses

[J]. Geophysics, 1984, 49(7): 861-869.

DOI:10.1190/1.1441732      URL     [本文引用: 1]

West R C, Ward S H.

The borehole transient electromagnetic response of a three-dimensional fracture zone in a conductive half-space

[J]. Geophysics, 1988, 53(11): 1469-1478.

DOI:10.1190/1.1442427      URL     [本文引用: 1]

Buselli G, Lee S K.

Modelling of drill-hole TEM responses from multiple targets covered by a conductive overburden

[J]. Exploration Geophysics, 1996, 27: 2-3, 141-153.

[本文引用: 1]

李建慧, 刘树才, 焦险峰, 等.

地—井瞬变电磁法三维正演研究

[J]. 石油地球物理勘探, 2015, 50(3):556-564.

[本文引用: 1]

Li J H, Liu S C, Jiao X F, et al.

Three-dimensional forward modeling for surfaceborehole transient electromagnetic method

[J]. Oil Geophysical Prospecting, 2015, 50(3): 556-564.

[本文引用: 1]

Dyck A V, West G F.

The role of simple computer models in interpretations of wide-band, drill-hole electromagnetic surveys in mineral exploration

[J]. Geophysics, 1984, 49(7): 957-980.

DOI:10.1190/1.1441741      URL     [本文引用: 1]

Bishop J R, Lewis J G, Macnae J C.

Down-hole electromagnetic surveys at Renison Bell, Tasmania

[J]. Exploration Geophysics, 1987, 18(3): 265-277.

DOI:10.1071/EG987265      URL     [本文引用: 1]

孟庆鑫, 潘和平.

地—井瞬变电磁响应特征数值模拟分析

[J]. 地球物理学报, 2012, 55(3):1046-1053.

[本文引用: 1]

Meng Q X, Pan H P.

Numerical simulation analysis of surface-hole TEM responses

[J]. Chinese Journal Geophysics, 2012, 55(3): 1046-1053.

[本文引用: 1]

徐正玉, 杨海燕, 邓居智, 等.

回线源三维地—井瞬变电磁法FDTD数值模拟

[J]. 工程地球物理学报, 2015, 12(3):327-332.

[本文引用: 1]

Xu Z Y, Yang H Y, Deng J Z, et al.

Three-dimensions FDTD numerical simulation on the down-hole TEM field with a loop source

[J]. Chinese Journal of Engineering Geophysics, 2015, 12(3): 327-332.

[本文引用: 1]

葛德彪, 闫玉波. 电磁波时域有限差分方法[M].2版. 西安: 西安电子科技大学出版社, 2005.

[本文引用: 1]

Ge D B, Yan Y B. Finite-difference time-domain method for electromagnetic waves[M].2 ed. Xi'an: Xidian University Press, 2005.

[本文引用: 1]

孙怀凤, 李貅, 李术才, 等.

考虑关断时间的回线源激发TEM三维时域有限差分正演

[J]. 地球物理学报, 2013, 56(3):1049-1064.

[本文引用: 1]

Sun H F, Li X, Li S C, et al.

Three-dimensional FDTD modeling of TEM excited by loop source considering ramp time

[J]. Chinese Journal Geophysics, 2013, 56(3): 1049-1064.

[本文引用: 1]

姚伟华, 王鹏, 李明星, 等.

三分量地—孔瞬变电磁法积水采空区探测试验

[J]. 煤田地质与勘探, 2019, 47(5):54-62.

[本文引用: 1]

Yao W H, Wang P, Li M X, et al.

Experimental study of three-component down-hole TEM for detecting water-filled goaf

[J]. Coal Geology & Exploration, 2019, 47(5): 54-62.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址：北京市学院路29号航遥中心 邮编：100083
电话：010-62060192；62060193 E-mail：whtbjb@sina.com