Integral equation solution for the transient electromagnetic response of a three‐dimensional body in a conductive half‐space

Author:

SanFilipo William A.1,Hohmann Gerald W.2

Affiliation:

1. CRA Exploration Pty. Ltd., 1030 Anderson Road, Karratha, W.A. 6714, Australia

2. Department of Geology and Geophysics, University of Utah, Salt Lake City, Utah 84112

Abstract

The time‐domain integral equation for the three‐dimensional vector electric field is formulated as a convolution of the scattering current with the tensor Green’s function. The convolution integral is divided into a sum of integrals over successive time steps, so that a numerical scheme can be formulated with a time stepping approximation of the convolution of past values of the solution with the system impulse response. This, together with spatial discretization, leads to a matrix equation in which previous solution vectors are multiplied by a series of matrices and fed back into the system by adding to the primary field source vector. The spatial discretization, based on a modification of the usual pulse basis formulation in the frequency domain, includes an additional subset of divergence‐free basis functions generated by integrating the Green’s function around concentric closed rectangular paths. The inductive response of the body is more accurately modeled with these additional basis functions, and a meaningful solution can be obtained for a body in free space. The resulting algorithm produces good results even for large conductivity contrasts. Internal checks, including convergence with respect to spatial and temporal discretization, and reciprocity, demonstrate self‐consistency of the numerical scheme. Independent checks include (a) comparison with results computed for a prism in free space, (b) comparison with results computed for a thin plate, (c) comparison of our conductive half‐space algorithm with an asymptotic solution for a sphere, and (d) comparison with results from inverse Fourier transformation of values computed using a frequency‐domain integral equation algorithm. Qualitative features of the results show that the relative importance of current channeling and confined eddy currents induced in the body depends upon both conductivity contrast and geometry. If the free‐space time constant is less than the time window during which currents in the host have not yet propagated well beyond the body, current channeling dominates the response. In such cases, simple superposition of free‐space results and the background is a poor approximation. In cases where the host currents diffuse beyond the body in a time less than the free‐space time constant of the body, the total response is approximately the sum of the free‐space and background (half‐space) responses.

Publisher

Society of Exploration Geophysicists

Subject

Geochemistry and Petrology,Geophysics

Cited by 72 articles. 订阅此论文施引文献 订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3