scholarly journals Reduction of Electromagnetic Reflections in 3D Airborne Transient Electromagnetic Modeling: Application of the CFS-PML in Source-Free Media

2018 ◽  
Vol 2018 ◽  
pp. 1-14 ◽  
Author(s):  
Yanju Ji ◽  
Xuejiao Zhao ◽  
Jiayue Gu ◽  
Dongsheng Li ◽  
Shanshan Guan
Keyword(s):  
Recursive Formula ◽  
Electromagnetic Response ◽  
Complex Frequency ◽  
Electromagnetic Modeling ◽  
Free Medium ◽  
Time Step ◽  
Homogeneous Half Space ◽  
Transient Electromagnetic ◽  
Anomalous Bodies ◽  
Relative Errors

To solve the problem of electromagnetic reflections caused by the termination of finite-difference time-domain (FDTD) grids, we apply the complex frequency-shifted perfectly matched layer (CFS-PML) to airborne transient electromagnetic (ATEM) modeling in a source-free medium. To implement the CFS-PML, two important aspects are improved. First, our method adopts the source-free Maxwell’s equations as the governing equations and introduces the divergence condition, consequently, the discrete form of Maxwell’s third equation is derived with regard to the CFS-PML form. Second, because our method adopts an inhomogeneous time-step, a recursive formula composed of convolution items based on a nonuniform time-step is proposed. The proposed approach is verified via a calculation of the electromagnetic response using homogeneous half-space models with different conductivities. The results show that the CFS-PML can reduce a 60 dB relative errors in late times. Moreover, this approach is also applied to 3D anomalous models; the results indicate that the proposed method can reduce reflections and substantially improve the identification of anomalous bodies. Consequently, the CFS-PML has good implications for ATEM modeling in a source-free medium.

2-D electromagnetic modeling by finite‐element method with a dipole source and topography

Geophysics ◽  
2000 ◽  
Vol 65 (2) ◽  
pp. 465-475 ◽  
Author(s):  
Yuji Mitsuhata
Keyword(s):  
Finite Element ◽  
Delta Function ◽  
Dipole Source ◽  
Electromagnetic Modeling ◽  
Secondary Field ◽  
Controlled Source ◽  
Transient Electromagnetic ◽  
Isoparametric Finite Element ◽  
Anomalous Bodies ◽  
Element Technique

I present a method for calculating frequency‐domain electromagnetic responses caused by a dipole source over a 2-D structure. In modeling controlled‐source electromagnetic data, it is usual to separate the electromagnetic field into a primary (background) and a secondary (scattered) field to avoid a source singularity, and only the secondary field caused by anomalous bodies is computed numerically. However, this conventional scheme is not effective for complex structures lacking a simple background structure. The present modeling method uses a pseudo‐delta function to distribute the dipole source current, and does not need the separation of the primary and the secondary field. In addition, the method employs an isoparametric finite‐element technique to represent realistic topography. Numerical experiments are used to validate the code. Finally, a simulation of a source overprint effect and the response of topography for the long‐offset transient electromagnetic and the controlled‐source magnetotelluric measurements is presented.

A finite‐element solution for the transient electromagnetic response of an arbitrary two‐dimensional resistivity distribution

Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1450-1461 ◽  
Author(s):  
Y. Goldman ◽  
C. Hubans ◽  
S. Nicoletis ◽  
S. Spitz
Keyword(s):  
Finite Element ◽  
Sparse Matrix ◽  
Equation System ◽  
Electromagnetic Response ◽  
Implicit Time ◽  
Two Dimensional ◽  
Time Step ◽  
Transient Electromagnetic ◽  
Exact Boundary ◽  
The Time Domain

We present a numerical method for solving Maxwell’s equations in the case of an arbitrary two‐dimensional resistivity distribution excited by an infinite current line. The electric field is computed directly in the time domain. The computations are carried out in the lower half‐space only because exact boundary conditions are used on the free surface. The algorithm follows the finite‐element approach, which leads (after space discretization) to an equation system with a sparse matrix. Time stepping is done with an implicit time scheme. At each time step, the solution of the equation system is provided by the fast system ICCG(0). The resulting algorithm produces good results even when large resistivity contrasts are involved. We present a test of the algorithm’s performance in the case of a homogeneous earth. With a reasonable grid, the relative error with respect to the analytical solution does not exceed 1 percent, even 2 s after the source is turned off.

Full transmitter waveform transient electromagnetic modeling and inversion for soundings over coal measures

Geophysics ◽  
1987 ◽  
Vol 52 (3) ◽  
pp. 279-288 ◽  
Author(s):  
Michael W. Asten
Keyword(s):  
Case Depth ◽  
Earth Model ◽  
Electromagnetic Modeling ◽  
Resistive Layer ◽  
Homogeneous Half Space ◽  
Full Waveform ◽  
Transient Electromagnetic ◽  
Electromagnetic Soundings ◽  
And Inversion ◽  
Coal Measures

Modeling of transient electromagnetic soundings requires consideration of the full bipolar transmitter waveform in order to obtain full measurement accuracy from field instruments. The full‐waveform transient may be computed using existing impulse‐response forward algorithms followed by a convolution procedure with truncation and residual estimation after an even number of terms, which results in average increases in run time of about 80 percent. Simple calculations based on approximations for the behavior of transient decay over (a) a conductive layer over a resistive basement, (b) a homogeneous half‐space, and (c) a resistive layer over a conductive basement show that the effect of the full waveform is greatest at sample times near the end of the transmitter off‐time, and for the case of a conductive basement. Failure to consider the full waveform is expected to yield errors in the amplitude of the computed transient in the range of 4 to 100 percent or more, depending upon the sample time and type of earth model. Examples of soundings over moderately conductive coal measures and highly conductive cindered coal show that failure to use the full waveform yields depths and resistivities biased to erroneously high values. The bias is small in the former case, but produces substantial errors in the latter case. Depth to the electrical basement of 300 to 400 m may be estimated to accuracies of approximately 5 percent by computer inversion to three‐layer or four‐layer models.

The harmonic and transient electromagnetic response of a thin dipping dike

Geophysics ◽  
1983 ◽  
Vol 48 (7) ◽  
pp. 934-952 ◽  
Author(s):  
P. Weidelt
Keyword(s):  
Formal Solution ◽  
Harmonic Excitation ◽  
Electromagnetic Response ◽  
Approximate Determination ◽  
Finite Conductivity ◽  
Circular Loop ◽  
Decay Modes ◽  
Transient Electromagnetic ◽  
Free Decay

An exact solution is given for the electromagnetic induction in a dipping dike of finite conductivity, represented as a thin half‐sheet in a nonconducting surrounding. The problem is formulated for arbitrary dipole or circular loop [Formula: see text] configurations. The formal solution obtained by the Wiener‐Hopf technique is cast into a rapidly convergent triple integral suitable for an effective numerical treatment. A good agreement is found between numerical results and analog measurements available for harmonic excitation. The transient response is obtained as a superposition of the half‐sheet free‐decay modes and is illustrated by some numerical examples for coincident loops, including a diagram for the approximate determination of conductance and depth of a vertical dike.

scholarly journals A parallel finite‐difference approach for 3D transient electromagnetic modeling with galvanic sources

Geophysics ◽  
2004 ◽  
Vol 69 (5) ◽  
pp. 1192-1202 ◽  
Author(s):  
Michael Commer ◽  
Gregory Newman
Keyword(s):  
Magnetic Field ◽  
Finite Difference ◽  
Electric Field ◽  
Electric Fields ◽  
Time Derivative ◽  
Stable Solution ◽  
Parallel Computer ◽  
Staggered Grid ◽  
Electromagnetic Modeling ◽  
Transient Electromagnetic

A parallel finite‐difference algorithm for the solution of diffusive, three‐dimensional (3D) transient electromagnetic field simulations is presented. The purpose of the scheme is the simulation of both electric fields and the time derivative of magnetic fields generated by galvanic sources (grounded wires) over arbitrarily complicated distributions of conductivity and magnetic permeability. Using a staggered grid and a modified DuFort‐Frankel method, the scheme steps Maxwell's equations in time. Electric field initialization is done by a conjugate‐gradient solution of a 3D Poisson problem, as is common in 3D resistivity modeling. Instead of calculating the initial magnetic field directly, its time derivative and curl are employed in order to advance the electric field in time. A divergence‐free condition is enforced for both the magnetic‐field time derivative and the total conduction‐current density, providing accurate results at late times. In order to simulate large realistic earth models, the algorithm has been designed to run on parallel computer platforms. The upward continuation boundary condition for a stable solution in the infinitely resistive air layer involves a two‐dimensional parallel fast Fourier transform. Example simulations are compared with analytical, integral‐equation and spectral Lanczos decomposition solutions and demonstrate the accuracy of the scheme.

scholarly journals Finite Difference Algorithm on Non-Uniform Meshes for Modeling 2D Magnetotelluric Responses

Algorithms ◽  
2018 ◽  
Vol 11 (12) ◽  
pp. 203
Author(s):  
Xiaozhong Tong ◽  
Yujun Guo ◽  
Wei Xie
Keyword(s):  
Finite Element ◽  
Finite Difference ◽  
Analytical Solutions ◽  
Good Accuracy ◽  
Validity And Reliability ◽  
Homogeneous Half Space ◽  
1D Model ◽  
Relative Errors ◽  
Future Project ◽  
Finite Difference Approach

A finite-difference approach with non-uniform meshes was presented for simulating magnetotelluric responses in 2D structures. We presented the calculation formula of this scheme from the boundary value problem of electric field and magnetic field, and compared finite-difference solutions with finite-element numerical results and analytical solutions of a 1D model. First, a homogeneous half-space model was tested and the finite-difference approach can provide very good accuracy for 2D magnetotelluric modeling. Then we compared them to the analytical solutions for the two-layered geo-electric model; the relative errors of the apparent resistivity and the impedance phase were both increased when the frequency was increased. To conclude, we compare our finite-difference simulation results with COMMEMI 2D-0 model with the finite-element solutions. Both results are in close agreement to each other. These comparisons can confirm the validity and reliability of our finite-difference algorithm. Moreover, a future project will extend the 2D structures to 3D, where non-uniform meshes should perform especially well.

scholarly journals Diagnosis for Conductor Breaks of Grounding Grids Based on the Wire Loop Method of the Transient Electromagnetic Method

2019 ◽  
Vol 2019 ◽  
pp. 1-11 ◽  
Author(s):  
Shengbao Yu ◽  
Guanliang Dong ◽  
Nannan Liu ◽  
Xiyang Liu ◽  
Chang Xu ◽  
...  
Keyword(s):  
Topological Structure ◽  
Electromagnetic Response ◽  
Height Difference ◽  
Measuring Point ◽  
Wire Loop ◽  
Loop Method ◽  
Transient Electromagnetic ◽  
Grounding Grid ◽  
Grounding Grids ◽  
The Wire

The wire loop method of the transient electromagnetic (TEM) method is used to nondestructively detect conductor breaks of grounding grid. For this purpose, grounding grids serve as an underground wire loop, and the measuring points are arranged on the ground. At each measuring point, a receiving loop is employed to detect the electromagnetic response generated by transmitting the current of the transmitting loop. Conductor breaks can be diagnosed by analyzing the slices of the electromagnetic response. We study the effect of loop size and height difference through the simulation of an intact 2×2 grounding grid, confirming that it is easier to obtain the topological structure using a small transmitting loop and a small height difference. Furthermore, simulations of an intact 4×4 grounding grid and grids with different locations of conductor breaks are also conducted with a small transmitting loop. It is easy to distinguish the topological structure of the grounding grid and the locations of conductor breaks. Finally, the detection method is applied experimentally. The experimental results confirm that the proposed method is an effective technique for conductor break diagnosis.

scholarly journals Numerical Simulation of Transient Electromagnetic Response of Unfavorable Geological Body in Tunnel

2011 ◽  
Vol 90-93 ◽  
pp. 37-40 ◽  
Author(s):  
Lu Bo Meng ◽  
Tian Bin Li ◽  
Zheng Duan
Keyword(s):  
Numerical Simulation ◽  
Finite Element ◽  
Electromagnetic Field ◽  
Electromagnetic Response ◽  
Detection Distance ◽  
Geological Body ◽  
Response Characteristics ◽  
Finite Element Software ◽  
Transient Electromagnetic ◽  
Field Decay

To investigate the transient electromagnetic method of response characteristics in the tunnel geological prediction, the finite element numerical simulation of unfavorable geological body of different location, different resistivity sizes, different shapes, and different volume size were carried out by ANSYS finite element software. The results show that secondary electromagnetic field of different location of unfavorable geological body have same decay rate, when detection distance from 30m to 70m, transient electromagnetic responses are strongest, followed distance from 10m to 30m and from 70m to 90m. The shape, volume and resistivity of unfavorable geological body have strong influence on transient electromagnetic response, unfavorable geological body more sleek, the greater the volume and the smaller the resistivity of unfavorable geological body, the secondary electromagnetic field decay slower.

Sign in / Sign up

Export Citation Format

Share Document