地孔瞬变电磁法超前探测数值模拟响应特征.pdf

第 44 卷第 10 期煤 炭 学 报Vol. 44 No. 10 2019 年10 月JOURNAL OF CHINA COAL SOCIETYOct. 2019 移动阅读 姚伟华,王鹏,李明星,等. 地孔瞬变电磁法超前探测数值模拟响应特征[J]. 煤炭学报,2019,44103145- 3153. doi10. 13225/ j. cnki. jccs. 2018. 1306 YAO Weihua,WANG Peng,LI Mingxing,et al. Numerical simulation response characteristics of down-hole TEM for advanced detection[J]. Journal of China Coal Society,2019,44103145-3153. doi10. 13225/ j. cnki. jccs. 2018. 1306 地孔瞬变电磁法超前探测数值模拟响应特征 姚伟华,王 鹏,李明星,苏 超,程思远 中煤科工集团西安研究院有限公司,陕西 西安 710077 摘 要为了避免传统井下瞬变电磁超前探测受巷道中电磁干扰信号、锚网、掘进机等影响,提出一 种在地面布设发射回线,在煤矿井下掘进工作面或工作面水平钻孔中测量瞬变电磁场三分量信号 的超前探测方法。 该方法具有径向探测深度大、三分量接收探头离异常体近、异常场占总场的比值 高而更易发现异常体等特点。 基于瞬变电磁三维时域有限差分正演算法,对均匀半空间模型、均匀 半空间含孤立异常体模型以及异常体位于水平钻孔不同方位的地孔瞬变电磁三分量的总场和异常 场响应特征进行了系统研究。 发现三分量总场和异常场之间总体符合烟圈效应,总场和异常场三 分量衰减曲线或者多测道曲线形态随时间的变化特征可以采用不同时刻感应电流环相对于测点或 者测线的位置来进行分析;在不考虑互感现象的情况下,地下接收的不同时刻的三分量响应完全可 以看成是背景电流环和异常电流环的叠加,其中异常场的三分量形态可以采用垂直于一次场方向 相同形状的电流环来拟合;另外当异常体位于水平钻孔不同方位时,异常场 X 分量均表现为 S 型, Y 分量和 Z 分量表现为单峰异常,且不同象限异常场三分量组合形态不同,据此可以判断异常体所 在钻孔的深度和方位。 最后在地面进行了地孔瞬变电磁超前探测方法的相似性试验,验证了理论 模型响应特征的正确性。 关键词地孔瞬变电磁法;三分量响应特征;超前探测 中图分类号P631． 3 文献标志码A 文章编号0253-9993201910-3145-09 收稿日期2018-09-29 修回日期2019-01-16 责任编辑韩晋平 基金项目国家重点研发计划资助项目2017YFC0804105;中国煤炭科工集团有限公司科技创新创业资金专项资助项目2018MS006;中国 煤炭科工集团西安研究院有限公司科技创新创业资金专项资助项目2017XAYZD12 作者简介姚伟华1989,男,陕西西安人,助理研究员,硕士。 E-mailyaoweihua cctegxian． com Numerical simulation response characteristics of down-hole TEM for advanced detection YAO Weihua,WANG Peng,LI Mingxing,SU Chao,CHENG Siyuan Xi’an Research Institute of China Coal Technology three-component response characteristics;advanced detection 我国煤炭开采以地下开采为主,随着煤矿向深部 延伸开采,煤矿开采的水文地质条件变得越发复杂, 水害防治形势愈发严峻。 煤矿井下水害事故主要发 生在巷道掘进期间和工作面回采期间[1]。 针对巷道 掘进前方灾害性地质体的超前物探方法主要有地震 反射波法、瑞雷波法、地质雷达法、矿井瞬变电磁法、 矿井直流电阻率法等[2]。 这些物探方法在以往实践 中能提前发现灾害性地质异常体,在保障矿井安全生 产方面发挥了巨大作用,但是存在探测距离不够长、 假异常多等缺点。 针对这些问题,研究一种在煤矿井 下水平钻孔中应用的地面-水平孔瞬变电磁三分量 超前探测技术。 该技术发射回线位于地表,接收探头 位于地下水平孔中,由于离异常体近,接收的异常响 应强、受电磁干扰小、并且接收三分量数据,能够实现 煤矿井下长距离探查掘进前方隐蔽水源性致灾体。 国外针对地孔瞬变电磁技术研究较早,WOODS D V[3]利用 CRONE 公司的 PEM 系统通过比例模型 实验研究了异常的纵向电导、倾角、几何尺寸和井旁 距离等参数对地井瞬变电磁三分量的影响特征,并总 结了一套地-井瞬变电磁定性解释原则。 EATON P A 和 HOHMANN G W[4],Richard WEST C 和 WARD S H[5]通过对不同导电围岩和低阻覆盖层的地井瞬变 电磁二维、三维模型进行正演计算,分析了导电围岩、 低阻覆盖层对地井瞬变电磁三分量数据的影响特征。 CULL J P[6]论证了瞬变电磁三分量相比于垂直 Z 分 量勘探的优点,最后证明采用地井瞬变电磁三分量勘 探能获得异常体的多个参数,包括倾角、方位角、空间 位置等。 国内、廊坊物化探研究所开展地井瞬变电磁 研究工作较早。 陈锡杰和任怀宗[7]从理论、模型实 验和资料解释 3 个方面对井中等轴状导体的瞬变电 磁响应进行了分析研究,并提出了位于回线磁源下方 的良导球体,对外产生的瞬变磁场可以用一个磁偶极 子来等效,为井中瞬变电磁异常解释提供了理论模 型。 张杰等[8]分别对地井瞬变电磁钻孔位置、薄板 的纵向电导及尺寸、倾角、埋深、回线磁源尺寸等参数 的变化与响应曲线变化之间的关系进行了研究,并提 出半定量的矢量交会技术。 杨毅等[9]采用遗传算法 开展了基于等效涡流的井中瞬变电磁纯异常反演研 究。 徐正玉等[10]对单异常体和多异常体以及在“定 源异井”和“动源定井”两种观测方式下典型地电模 型的响应特征进行了研究。 杨怀杰等[11]对井中不同 导电围岩和含低阻体的电磁场的传播特征进行了研 究。 杨海燕等[12]对均匀半空间和倾斜板状体、含有 低阻覆盖层的倾斜板状体进行了模拟研究。 李术才 等[13]对深部开采矿井中地面电性源发射,井下接收 的瞬变电磁探测方法的响应规律展开研究,分析了不 同收发距和方位角的曲线形态,总结电性源地-井瞬 变电磁响应传播规律。 武军杰等[14]对电性源地井瞬 变电磁法进行了研究,通过数值模拟对三分量曲线形 态进行分析,总结了电性发射源条件下钻孔中地下电 性界面的异常特征。 范涛[15]提出一种巷道孔口旋转 发射、孔内三分量接收的工作方法,结合山西煤矿井 下的工程试验,对该方法超前探测采空区性质、形态 和规模的实用性和有效性进行了检验。 笔者主要从理论模拟方面对地孔瞬变电磁超前 探测方法的三分量总场背景地层加异常体的二次 场响应和异常场响应异常体的二次场响应特征 进行了研究,着力于探讨异常场与总场的响应关系及 传播特征,为该方法后续的精细解释及生产实践提供 理论基础。 6413 第 10 期姚伟华等地孔瞬变电磁法超前探测数值模拟响应特征 1 地孔瞬变电磁超前探测方法 地孔瞬变电磁超前探测方法是指在地面布设发 射回线,在井下掘进工作面或工作面水平钻孔中测量 瞬变电磁场三分量信号,通过对测量的三分量瞬变电 磁信号来探测水平钻孔周围低阻异常体的一种超前 探测方法,工作装置示意如图 1 所示。 相比小回线瞬 变电磁超前探测法,地孔瞬变电磁法由于发射源布设 于地表,一次场能量强,利于激发更强的二次场;在进 行地孔瞬变电磁超前探测时,接收探头位于水平钻孔 中,更接近异常体,测量的信号中,异常场占总场的比 值高;再者由于探头位于钻孔中且远离巷道空间,接 收的数据不易受巷道中电磁干扰信号、锚网、掘进机 等影响,利于发现异常体。 在施工过程中,发射回线 铺设于地表,位于测量钻孔的正上方,尽可能使异常 体最大程度被激发;接收机及接收探头位于井下,测 点点距根据异常体规模而定,一般设置为 1 m。 接收 机与发射机之间可以采用石英钟同步,两者之间无需 通信,发射回线从施工开始至施工结束持续发射正负 交替的梯形波,接收机在井下独立开展测量工作即 可,施工方便。 地孔瞬变电磁超前探测方法其实是地 井瞬变电磁法的另一种形式,区别只是接收点位于掘 进工作面前方的水平孔中。 图 1 地孔瞬变电磁法超前探测示意 Fig． 1 Down-hole TEM advanced detection schematic 2 地孔瞬变电磁三维正演原理 在直角坐标系中将麦克斯韦微分方程组写成各 分量的形式[16]为 ∂Ez ∂y - ∂Ey ∂z - ∂Bx ∂t ∂Ex ∂z - ∂Ez ∂x - ∂By ∂t ∂Ey ∂x - ∂Ex ∂y - ∂Bz ∂t 1 ∂Hz ∂y - ∂Hy ∂z γ ∂Ex ∂t σEx Jsx ∂Hx ∂z - ∂Hz ∂x γ ∂Ey ∂t σEy Jsy ∂Hy ∂x - ∂Hx ∂y γ ∂Ez ∂t σEz 2 式中,γ 与介电常数具有相同的量纲,称为虚拟介电 常数;包含 γ 的项具有电流的量纲,称为虚拟位移电 流;Js代表源电流密度,A/ m2;x,y,z 为 3 个坐标,m; Ex,Ey,Ez为电场强度,V/ m;Hx,Hy,Hz为三分量磁场 强度,A/ m;Bx,By,Bz为三分量磁感应强度,T;S 为电 导率,S/ m。 采用 Yee 晶胞格式对式1,2进行网格离散, 每一个电场磁场分量均由 4 个磁场电场分量包 围,这样的电场、磁场空间分布形式自然符合法拉第 电磁感应定律和安培环路定理的结构形式,同时也满 足 Maxwell 方程组的差分计算要求。 3 地孔瞬变电磁法总场特征 为了得到地孔瞬变电磁法在超前探测中响应的 三分量规律,利用时域有限差分法进行数值试验。 设 置地面回线源尺寸为600 m600 m,电流20 A。 坐标 方向按右手螺旋定律确定,坐标原点位于地面回线中 心。 测线沿 X 增量方向布置在水平煤层中,测点点 距 2 m。 设置背景电阻率为 100 Ωm,在测线上接 收瞬变电磁法二次磁场的 3 个分量,时间窗口为 0． 01 10 ms,时间道数为 100 道。 与地面瞬变电磁法不同,地孔瞬变电磁法接收测 线位于地下,那么测量的三分量响应曲线形态也不 同,为了研究地孔瞬变电磁法的衰减规律,设计半空 间模型,测线位于地下 300 m, 在 Y 轴上的投影 为-10 m,在 X 轴上的投影为[-100,100],计算模型 如图 2 所示。 图 2 均匀半空间模型俯视图 Fig． 2 Uni half-space model top view 图 3 为计算的均匀半空间模型的磁场三分量衰 7413 煤 炭 学 报 2019 年第 44 卷 减曲线与多测道曲线图,其中图3a c计算的衰 减曲线的测点坐标为10,-10,300。 从三分量衰减 曲线可以看出,对水平 X 分量在时间上是由正变负 进行衰减,Y 分量由负变正进行衰减,Z 分量单调递 减。 对于整条测线而言,X 分量成“X”型,Y 分量和 Z 分量呈单峰形态,只是符号相反。 图 3 均匀半空间模型磁场响应曲线 Fig． 3 Magnetic field response curves of a uni half-space model 为便于理解以上描述的现象,可假定地下在断电 后出现一个电流环,该电流环代表了背景空间的整体 涡流,在断电瞬间,它的电流方向和形状与地面发射 框一致,其强度随时间不断减弱,深度不断下降,直径 不断增大,直至能量完全衰减。 在此假定条件下,空 间其他各点的磁场均由该电流环产生,电流环产生的 磁力线按右手螺旋原则确定的方向,对于三分量多测 道曲线,在空间上完全符合单个电流环产生的磁场形 态。 对于衰减曲线,在早期,地下产生的电流环位于 测线上部,其电流方向与发射回线电流方向一致,如 图4 中 t0时刻,此时对于磁场 X 分量与坐标系定义 X 方向相同,所以 X 分量在早期为正,磁场 Z 分量与坐 标系定义 Z 方向一致,所以为正。 随着时间的推移, 电流环向深部转移,直径逐渐增大,强度逐渐减弱。 到晚期时,如图 4 中 t1时刻,电流环位于测线下方, 磁场 X 分量与坐标系定义 X 方向相反,表现为负值, 磁场 Z 分量与坐标系 Z 方向一致,表现为正值,Y 分 量与 X 分量衰减规律一致,在此不再赘述。 4 地孔瞬变电磁法异常场特征 4． 1 异常体感应的异常场特征 为了进一步研究地孔瞬变电磁法三维异常体感 图 4 电流环传播示意 Fig． 4 Schematic diagram of current loop propagation 应涡流磁场随时间的分布规律,利用瞬变电磁三维时 域有限差分程序正演计算异常体三分量响应。 在上 述均匀半空间模型中加入一低阻异常体,异常体的电 阻率为 1 Ωm,规模为 40 m40 m40 m,中心埋深 200 m,位于回线磁源的正下方。 计算了不同时刻 y 0 的 XZ 平面和 Z200 的 XY 平面三分量磁场响应, 将计算的含异常体的三分量响应减去相同背景均匀 全空间的三分量响应,获得异常场三分量响应。 绘制 XZ 平面和 XY 平面异常场三分量矢量分布图,如图 5,6 所示。 由图 5,6 可以看出,一次场消失后,为了维持异 8413 第 10 期姚伟华等地孔瞬变电磁法超前探测数值模拟响应特征 图 5 不同时刻异常场 XZ 平面三分量矢量分布 Fig． 5 Three-component vector distribution of XZ plane of anomalous field at different times 图 6 不同时刻异常场 XY 平面三分量矢量分布 Fig． 6 Three-component vector distribution of XY plane of anomalous field at different times 常体内原来的均匀磁场,异常体内产生感应涡流。 在 XZ 平面早期感应涡流的磁力线围绕两个圆心分布, 整体沿轴 X0 对称。 在 XY 平面早期感应涡流的磁 场由异常体中心指向外,形态呈基本规则的圆环状。 9413 煤 炭 学 报 2019 年第 44 卷 随着时间的推移,XZ 平面感应涡流磁场 X 分量与 Z 分量场强的大小及方向在不断变化,磁力线的分布仍 持续稳定。 XY 平面由中心指向外部的磁场逐渐转 向,到中晚期感应涡流的磁场由外部指向异常体中 心,并基本保持稳定。 除图 5,6 显示的极早期t 0． 038 ms时间以前,低阻异常体感应涡流的磁场在 晚期与水平电流环在 XZ 平面和 XY 平面磁场的分布 特征有明显的相似性。 为了定量的比较异常体涡流场某个时刻与电流 环产生磁场具有等效性,设置测线位于异常体中心正 下方 50 m,正演计算其异常场三分量响应,选取几个 特定时刻比较感应涡流场辐射的磁场与半径为 20 m异常体边长为 40 m水平电流环产生的磁场。 比较时将感应涡流场的磁场与电流环产生磁场的三 分量分别做幅值归一化。 图 7 显示两者三分量的重 合程度,图中用实心圆点代表电流环的磁场,用实线 代表涡流场的磁场,红色表示 X 分量,蓝色表示 Y 分 量,黑色表示 Z 分量。 图 7 不同时刻电流环三分量与感应涡流三分量对比 Fig． 7 Comparison of three components of current loop and induced eddy current at different times 总体来看,图 7 中涡流场的磁场三分量曲线与电 流环产生的磁场曲线重合度较高,也就是说感应涡流 某些时刻产生的磁场可以用电流环产生的磁场来等 效。 NABIGHIAN1979 [17]的烟圈理论指出对于均 匀半空间模型,同一时刻的瞬变电磁二次场可以用一 个电流环的响应等效。 那么在此基础上我们可以做 类似等效,断电瞬间地下产生两类电流环,分别代表 大地背景场和异常场。 背景场电流环的目的是维持 大框内部的磁场,异常场电流环的目的是维持低阻体 内部的初始磁场。 随着时间推移,背景场电流环直径 不断扩大,深度不断下降,电流强度不断减弱。 异常 场电流环的电流强度随时间减弱,但其位置固定,一 直在低阻体内部衰减。 背景电流环只有一个,异常场 电流环可以有多个,每个低阻体内部存在且只存在一 个异常电流环。 如果不考虑互感现象,空间各点的总 磁场为上述背景电流环和异常电流环产生磁场的合 成。 在地孔瞬变电磁勘探中,可以根据此假设来分析 地井瞬变电磁三分量的曲线变化特征。 4． 2 不同象限异常体的响应特征 设置含水异常体分别位于测线截面的 4 个不同 的方位上,如图 8 所示,通过数值模拟获得不同方位 低 阻 异 常 体 的 感 应 二 次 场。 背 景 电 阻 率 为 100 Ωm,异常体电阻率为 1 Ω m, 尺寸为 20 m20 m20 m,中心点在水平钻孔上的投影距离 分别为X-24 m,Y-15 m,Z15 m,水平钻孔位于 地下 300 m 处。 在这 4 组模型中,只通过改变异常体 方位来分析异常体位于钻孔不同方位时异常场的变 化规律。 图 8 不同方位异常体模型示意 Fig． 8 Schematic diagram of different azimuth abnormal body models 图 9 为不同方位异常体正演后获得的感应二次 异常场的 3 个分量,其中图 9a c分别为 1 号异 常体的 X,Y,Z 分量异常测道图,图 9d f分别 为 2 号 异 常 体 的 X, Y, Z 分 量 异 常 测 道 图, 图 9g i分别为 3 号异常体的 X,Y,Z 分量异常 测道图,图 9j l分别为 4 号异常体的 X,Y,Z 分 量异常测道图。 对 4 个异常体产生的异常场而言,异常场的分量 均表现为先正后负或者先负后正的双峰形态,异常场 的 Z 分量均表现为正值或者负值隆起的单峰形态。 0513 第 10 期姚伟华等地孔瞬变电磁法超前探测数值模拟响应特征 图 9 不同方位异常体模型三分量测道图 Fig． 9 Three-component track map of different azimuth abnormal body models X 分量多测道曲线过 0 点位置和 Y,Z 分量极值位置 与异常体中心位置在钻孔上的投影吻合。 不同的是 对于 X 分量,1 号异常和 2 号异常体的表现为先正后 负,3 号和 4 号异常表现为先负后正,对于 Y 分量 1 号异常和 3 号异常表现为正极值异常,2 号和 4 号表 现为负极值异常,对于 Z 分量,4 个模型均表现为正 极值异常。 仔细对比 4 个模型的异常场三分量发现, 异常体位于钻孔不同方位时产生的异常场三分量曲 线组合形态各不相同,据此可以确定异常体位于水平 钻孔的深度及异常体所在钻孔的方位。 对于不同象限异常体的响应特征可以从烟圈理 论来解释,当发射回路断电后,异常体为了维持原来 的一次场,感应出等效涡流即电流环,当异常体位于 钻孔不同方位时,异常体的电流环产生的二次场三分 量响应在测线上的曲线形态组合不同,这说明利用地 孔瞬变电磁的异常场能实现异常体的空间定位。 5 物理试验 为了研究地孔瞬变电磁超前探测的响应特征,在 陕西省蓝田县郊外开展了地孔相似性物理试验,试验 区无高压线,地形平坦,最大落差不超过 5 m。 试验 采用的发射机为 GT30,接收机为 SM24,探头为磁通 门三分量探头。 试验过程中发射回线边长为 60 m 60 m,发射基频为 5 Hz,发射电流为 13 A,电流为逆 时针方向。 异常体为 5 块 7 m1． 2 m 的铁皮拼接而 成,尺寸为 5 m7 m,异常体测线距金属铁板中心水 平距离为5． 3 m,垂直距离为5 m,试验现场照片和测 线布置如图 10,11 所示。 测点点距为 1 m,测线 1 上 有 41 个点,点号编排方式为 1 41。 坐标系规定向 右为 X 正向,向南为 Y 正向,向上为 Z 正向。 图 10 物理试验现场照片 Fig． 10 Physical test site photos 实测的磁场三分量如图 12 所示。 从 12 a c三分量衰减曲线可以看出,对水平 X 分量、Y 1513 煤 炭 学 报 2019 年第 44 卷 图 11 物理试验布置 Fig． 11 Physical test layout 分量和 Z 分量均为单调递减,区别只是 X 分量符号 相反。 对于整条测线的多测道曲线而言,X 分量先负 后正,Y 分量和 Z 分量呈单峰形态。 对于上描述的现象,也可假定地下在断电后在异 常体内和大地中分别出现一个电流环,两个电流环的 感应电流方向与发射回线一致,只是强弱不同,由于 异常体为强低阻,所以异常体内感应的电流环产生的 磁场占主导因素。 测线 1 位于地表,所以一直处于电 流环的上方,所以三分量衰减曲线随时间衰减不会变 号。 在异常体中心的小号点方向,按照右手螺旋定理 推算,X 分量的方向与坐标系 X 分量相反,所以为负; 在异常体中心的大号点方向,X 分量的方向与坐标系 X 分量相同,所以为正。 整条测线上 Y 分量和 Z 分量 方向与坐标系一致均为正。 与第 3 节的理论模型响 应进行对比发现,两者衰减曲线随时间的衰减规律不 一致,原因是理论模型的测线位于地下,随着时间衰 减,感应电流环在早期位于测线上方,在晚期位于测 线下方,水平分量衰减曲线随时间推移出现变号现 象,而物理试验的测线位于地表,感应电流环一致位 于测线下方,所以三分量衰减曲线随时间递推不会出 现变号现象。 综上所述,地孔瞬变电磁物理试验三分 量响应规律与理论模型三分量响应规律一致,两者相 互印证。 图 12 测线 1 磁场响应曲线 Fig． 12 Magnetic field response curve of the Line 1 6 结 论 1通过对地孔瞬变电磁总场和异常场的三分 量响应进行三维数值模拟发现,地孔瞬变电磁三分量 响应随时间的形态变化规律完全遵循“烟圈理论”。 在不考虑异常体与背景地层互感的情况下,地下接收 的不同时刻的瞬变电磁三分量响应可以看成是背景 电流环和异常电流环的叠加,背景电流环只有一个, 异常场电流环可以有多个。 据此可以来分析地孔瞬 变电磁总场和异常场三分量响应特征。 2通过对水平钻孔周围不同方位的异常体进 行正演与物理试验,发现当异常体位于钻孔不同方位 2513 第 10 期姚伟华等地孔瞬变电磁法超前探测数值模拟响应特征 时,其产生的异常场三分量曲线组合形态各不相同。 该现象表明通过在地下水平钻孔中测量瞬变电磁三 分量响应,就能根据异常场形态组合判别异常体相对 于钻孔的方位与深度,实现掘进工作面前方低阻异常 体的超前探测。 3在瞬变电磁三维反演技术还不完全成熟的 情况下,虽然异常场三分量能实现钻孔周围异常体的 精确定位,但是在实际生产中如何根据实测的总场三 分量来提取准确的异常场三分量是该方法后续研究 的一个重要方向。 参考文献References [1] 刘志民,刘希高,张金涛,等. 交流聚焦激电法煤巷超前探测阻 容试验模拟[J]. 煤炭学报,2015,4092144-2151. LIU Zhiming,LIU Xigao,ZHANG Jintao,et al. Experimental simula- tion of resistance-capacitance model for advanced detection in coal roadway based on alternating current focusing induced polariza- tion [J]. Journal of China Coal Society,2015,4092144- 2151. [2] 程久龙,邱浩,叶云涛,等. 矿井瞬变电磁法波场变换与数据处 理方法研究[J]. 煤炭学报,2013,3891646-1650. CHENG Jiulong,QIU Hao,YE Yuntao,et al. Research on wave- field transation and data processing of the mine transient electro- magnetic [J]. Journal of China Coal Society,2013,389 1646-1650. [3] WOODS D V. A model study of the crone borehole pulse electromag- neticPEM system[D]. KingstonQueen’s University,1975. [4] EATON P A,HOHMANN G W. Influence of a conductive host on two-dimensional borehole transient electromagnetic responses [J]. Geophysics,1984,497861-869. [5] WEST R C,WARD S H. Borehole transient electromagnetic response of a three-dimensional fracture zone in a conductive half-space [J]. Geophysics,1988,53111469-1478. [6] CULL J P. Rotation and resolution of three-component DHEM data [J]. Exploration Geophysics,1996,273155-159. [7] 陈锡杰,任怀宗. 钻井中等轴状导体的瞬变电磁响应[J]. 地质 与勘探,19881139-43. CHEN Xijie,REN Haizong. Responses of a buried equiaxial condu- etive body in drill hole transient eleetromagnetie measurements [J]. Geology and Prospecting,19881139-43. [8] 张杰,吕国印,赵敬洗,等. 地-井 TEM 向量交会技术的实现和应 用效果[J]. 物探化探计算技术,2007,29162-166. ZHANG Jie,L Guoying,ZHAO Jingxi,et al. The of surface- borehole TEM vector intersection and its application[J]. Computing Techniques for Geophysical Geochemical Exploration,2007,29 162-166. [9] 杨毅,邓晓红,张杰,等. 一种井中瞬变电磁异常反演方法[J]. 物探与化探,2014,384855-859,864. YANG Yi,DENG Xiaohong,ZHANG Jie,et al. A borehole TEM a- nomaly inversion [J]. Geophysical and Geochemical Explora- tion,2014,384855-859,864. [10] 徐正玉,杨海燕,邓居智,等. 基于异常场的地-井瞬变电磁法 正演研究[J]. 物探与化探,2015,3961176-1182. XU Zhengyu,YANG Haiyan,DENG Juzhi,et al. Research on for- ward simulation of down-hole TEM based on the abnormal field [J]. Geophysical and Geochemical Exploration,2015,39 6 1176-1182. [11] 杨怀杰,潘和平,孟庆鑫,等. 导电围岩对井中三维瞬变电磁响 应的影响规律研究[J]. 石油物探,2016,552288-293,302. YANG Huaijie,PAN Heping,MENG Qingxin,et al. Influence laws of conductive host on borehole 3D transient electromagnetic respon- ses[J]. Geophysical Prospecting for Petroleum,2016,552288- 293,302. [12] 杨海燕,岳建华,徐正玉,等. 覆盖层影响下典型地-井模型瞬 变电磁法正演[J]. 吉林大学学报地球科学版,2016,465 1527-1537. YANG Haiyan,YUE Jianhua,XU Zhengyu,et al. Transient electro- magnetic modeling in ground-borehole model with overbur- den influence[J]. Journal of Jilin UniversityEarth Science Edi- tion,2016,4651527-1537. [13] 李术才,李凯,翟明华,等. 矿井地面-井下电性源瞬变电磁探 测响应规律分析[J]. 煤炭学报,2016,4182024-2032. LI Shucai,LI Kai,ZHAI Minghua,et al. Analysis of grounded tran- sient electromagnetic with surface-tunnel configuration in mining [J]. Journal of China Coal Society,2016,4182024-2032. [14] 武军杰,李貅,智庆全,等. 电性源地-井瞬变电磁异常场响应 特征初步分析[J]. 物探与化探,2017,411129-135. WU Junjie,LI Xiu,ZHI Qingquan,et al. A preliminary analysis of anomalous TEM response characteristics in borehole with electric source transmitter[J]. Geophysical Geochemical Exploration, 2017,411129-135. [15] 范涛. 孔巷瞬变电磁动源定接收方法探测采空区试验[J]. 煤 炭学报,2017,42123229-3238. FAN Tao. Experimental study on the exploration of coal mine goaf by dynamic source and fixed reception Roadway-Borehole TEM de- tection [J]. Journal of China Coal Society,2017,4212 3229-3238. [16] 孙怀凤,李貅,李术才,等. 考虑关断时间的回线源激发 TEM 三 维时域有限差分正演[J]. 地球物理学报,2013,5631049- 1064. SUN Huaifeng, LIXiu, LIShucai, etal. Three-dimensional FDTD modeling of TEM excited by a loop source considering ramp time[J]. CHINESE J. Geophys. ,2013,5631049-1064. [17] NABIGHIAN M N. Quasi-static transient response of a uniaxial conducting half-spaceAn approximate representation[J]. Geo- physics,1979,4410380-414. 3513