The Solution of Transient Electromagnetic Inverse Source Problems Using Time-Domain TLM Method

2012 ◽  
Vol 60 (9) ◽  
pp. 4326-4335
Author(s):  
Yu Zhang ◽  
Mohamed H. Bakr ◽  
Natalia K. Nikolova
Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. E481-E491 ◽  
Author(s):  
Andrei Swidinsky ◽  
Misac Nabighian

Electromagnetic surveys using a vertical transmitter loop are common in land, marine, and airborne geophysical exploration. Most of these horizontal magnetic dipole (HMD) systems operate in the frequency domain, measuring the time derivative of the induced magnetic fields, and therefore a majority of studies have focused on this subset of field measurements. We examine the time-domain electromagnetic response of a HMD including the electric fields and corresponding smoke rings produced in a conductive half-space. Cases of a dipole at the surface and buried within the earth are considered. Results indicate that when the current in the transmitter is rapidly switched off, a single smoke ring is produced within the plane of the vertical transmitter loop, which is then distorted by the air-earth interface. In this situation, the circular smoke ring, which would normally diffuse symmetrically away from the source in a whole space, is approximately transformed into an ellipse, with a vertical major axis at an early time and a horizontal major axis at a late time. As measured from the location of the transmitter, the depth of investigation and lateral footprint of such a system increases with burial depth. It is also observed that the electric field measured in the direction of the magnetic dipole only contains a secondary response related to the charge accumulation on any horizontal conductivity boundaries because the primary field is always absent. This field component can be expressed analytically in terms of a static and time-varying field, the latter term adding spatial complexity to the total horizontal electric field at the earth surface at early times. Applications of this theoretical study include the design of time-domain induction-logging tools, crossborehole electromagnetic surveys, underground mine expansion work, mine rescue procedures, and novel marine electromagnetic experiments.


2021 ◽  
Vol 36 (3) ◽  
pp. 245-251
Author(s):  
Jun Li ◽  
Huaguang Bao ◽  
Dazhi Ding

In order to evaluate scattering from hypersonic vehicles covered with the plasma efficiently, time domain volume shooting and bouncing rays (TDVSBR) is first introduced in this paper. The new method is applied to solve the transient electromagnetic scattering from complex targets, which combines with non-homogeneous dielectric and perfect electric conducting (PEC) bodies. To simplify the problem, objects are discretized into tetrahedrons with different electromagnetic parameters. Then the reflection and transmission coefficients can be obtained by using theory of electromagnetic waves propagation in lossy medium. After that, we simulate the reflection and transmission of rays in different media. At last, the scattered fields or radiation are solved by the last exiting ray from the target. Compared with frequency-domain methods, time-domain methods can obtain the wideband RCS efficiently. Several numerical results are given to demonstrate the high efficiency and accuracy of this proposed scheme.


Geophysics ◽  
2021 ◽  
pp. 1-43
Author(s):  
Qingtao Sun ◽  
Runren Zhang ◽  
Yunyun Hu

To facilitate the modeling of time-domain controlled-source electromagnetic survey, we propose an efficient finite-element method with weighted Laguerre polynomials, which shows a much lower computational complexity than conventional time integration methods. The proposed method allows sampling the field at arbitrary time steps and also its accuracy is determined by the number of polynomials, instead of the time sampling interval. Analysis is given regarding the optimization of the polynomial number to be used and the criterion of selecting the time scale factor. Two numerical examples in marine and land survey environments are included to demonstrate the superiority of the proposed method over the existing backward Euler time integration method. The proposed method is expected to facilitate the modeling of transient electromagnetic surveys in the geophysical regime.


Geophysics ◽  
1992 ◽  
Vol 57 (8) ◽  
pp. 994-1003 ◽  
Author(s):  
Michael Leppin

A numerical method is presented by which the transient electromagnetic response of a two‐dimensional (2-D) conductor, embedded in a conductive host rock and excited by a rectangular current loop, can be modeled. This 2.5-D modeling problem has been formulated in the time domain in terms of a vector diffusion equation for the scattered magnetic induction, which is Fourier transformed into the spatial wavenumber domain in the strike direction of the conductor. To confine the region of solution of the diffusion equation to the conductive earth, boundary values for the components of the magnetic induction on the ground surface have been calculated by means of an integral transform of the vertical component of the magnetic induction at the air‐earth interface. The system of parabolic differential equations for the three magnetic components has been integrated for 9 to 15 discrete spatial wavenumbers ranging from [Formula: see text] to [Formula: see text] using an implicit homogeneous finite‐difference scheme. The discretization of the differential equations on a grid representing a cross‐section of the conductive earth results in a large, sparse system of linear equations, which is solved by the successive overrelaxation method. The three‐dimensional (3-D) response has been computed by an inverse Fourier transformation of the cubic spline interpolated scattered magnetic induction in the wavenumber domain using a digital filtering technique. To test the algorithm, responses have been computed for a two‐layered half‐space and a vertical prism embedded in a conductive host rock. These examples were then compared with results obtained analytically or numerically using frequency‐domain finite‐element and time‐domain integral equation methods. The new numerical procedure gives satisfactory results for a wide range of 2-D conductivity distributions with conductivity ratios exceeding 1:100, provided the grid is sufficiently refined at the corners of the conductivity anomalies.


Sign in / Sign up

Export Citation Format

Share Document