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

Geophysics ◽  
1985 ◽  
Vol 50 (5) ◽  
pp. 798-809 ◽  
Author(s):  
William A. SanFilipo ◽  
Gerald W. Hohmann
Keyword(s):  
Integral Equation ◽  
Frequency Domain ◽  
Half Space ◽  
Time Constant ◽  
Free Space ◽  
Numerical Scheme ◽  
Three Dimensional ◽  
Basis Functions ◽  
The Body ◽  
Spatial Discretization

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.

The effect of a conductive half‐space on the transient electromagnetic response of a three‐dimensional body

Geophysics ◽  
1985 ◽  
Vol 50 (7) ◽  
pp. 1144-1162 ◽  
Author(s):  
William A. SanFilipo ◽  
Perry A. Eaton ◽  
Gerald W. Hohmann
Keyword(s):  
Half Space ◽  
Free Space ◽  
Eddy Currents ◽  
Three Dimensional ◽  
Vertical Component ◽  
The Body ◽  
Electromagnetic Response ◽  
Loop Current ◽  
Transient Electromagnetic ◽  
Space Response

The transient electromagnetic (TEM) response of a three‐dimensional (3-D) prism in a conductive half‐space is not always approximated well by three‐dimensional free‐space or two‐dimensional (2-D) conductive host models. The 3-D conductive host model is characterized by a complex interaction between inductive and current channeling effects. We numerically computed 3-D TEM responses using a time‐domain integral‐equation solution. Models consist of a vertical or horizontal prismatic conductor in conductive half‐space, energized by a rapid linear turn‐off of current in a rectangular loop. Current channeling, characterized by currents that flow through the body, is produced by charges which accumulate on the surface of the 3-D body and results in response profiles that can be much different in amplitude and shape than the corresponding response for the same body in free space, even after subtracting the half‐space response. Responses characterized by inductive (vortex) currents circulating within the body are similar to the response of the body in free space after subtracting the half‐space contribution. The difference between responses dominated by either channeled or vortex currents is subtle for vertical bodies but dramatic for horizontal bodies. Changing the conductivity of the host effects the relative importance of current channeling, the velocity and rate of decay of the primary (half‐space) electric field, and the build‐up of eddy currents in the body. As host conductivity increases, current channeling enhances the amplitude of the response of a vertical body and broadens the anomaly along the profile. For a horizontal body the shape of the anomaly is distorted from the free‐space anomaly by current channeling and is highly sensitive to the resistivity of the host. In the latter case, a 2-D response is similar to the 3-D response only if current channeling effects dominate over inductive effects. For models that are not greatly elongated, TEM responses are more sensitive to the conductivity of the body than galvanic (dc) responses, which saturate at a moderate resistivity contrast. Multicomponent data are preferable to vertical component data because in some cases the presence and location of the target are more easily resolved in the horizontal response and because the horizontal half‐space response decays more quickly than does the corresponding vertical response.

Transient electromagnetic response of a three‐dimensional body in a layered earth

Geophysics ◽  
1986 ◽  
Vol 51 (8) ◽  
pp. 1608-1627 ◽  
Author(s):  
Gregory A. Newman ◽  
Gerald W. Hohmann ◽  
Walter L. Anderson
Keyword(s):  
Integral Equation ◽  
Frequency Domain ◽  
Digital Filters ◽  
Three Dimensional ◽  
The Body ◽  
Late Time ◽  
Near Surface ◽  
Decay Spectrum ◽  
Transient Electromagnetic ◽  
Central Loop

The three‐dimensional (3-D) electromagnetic scattering problem is first formulated in the frequency domain in terms of an electric field volume integral equation. Three‐dimensional responses are then Fourier transformed with sine and cosine digital filters or with the decay spectrum. The digital filter technique is applied to a sparsely sampled frequency sounding, which is replaced by a cubic spline interpolating function prior to convolution with the digital filters. Typically, 20 to 40 frequencies at five to eight points per decade are required for an accurate solution. A calculated transient is usually in error after it has decayed more than six orders in magnitude from early to late time. The decay spectrum usually requires ten frequencies for a satisfactory solution. However, the solution using the decay spectrum appears to be less accurate than the solution using the digital filters, particularly after early times. Checks on the 3-D solution include reciprocity and convergence checks in the frequency domain, and a comparison of Fourier‐transformed responses with results from a direct time‐domain integral equation solution. The galvanic response of a 3-D conductor energized by a large rectangular loop is substantial when host currents are strong near the conductor. The more conductive the host, the longer the galvanic responses will persist. Large galvanic responses occur if a 3-D conductor is in contact with a conductive overburden. For a thin vertical dike embedded within a conductive host, the 3-D response is similar in form but differs in magnitude and duration from the 2-D response generated by two infinite line sources positioned parallel to the strike direction of the 2-D structure. We have used the 3-D solution to study the application of the central‐loop method to structural interpretation. The results suggest variations of thickness of conductive overburden and depth to sedimentary structure beneath volcanics can be mapped with one‐dimensional inversion. Successful 1-D inversions of 3-D transient soundings replace a 3-D conductor by a conducting layer at a similar depth. However, other possibilities include reduced thickness and resistivity of the 1-D host containing the body. Many different 1-D models can be fit to a transient sounding over a 3-D structure. Near‐surface, 3-D geologic noise will not permanently contaminate a central‐loop apparent resistivity sounding. The noise is band‐limited in time and eventually vanishes at late times.

THREE‐DIMENSIONAL INDUCED POLARIZATION AND ELECTROMAGNETIC MODELING

Geophysics ◽  
1975 ◽  
Vol 40 (2) ◽  
pp. 309-324 ◽  
Author(s):  
Gerald W. Hohmann
Keyword(s):  
Integral Equation ◽  
Green's Functions ◽  
Green’S Functions ◽  
Three Dimensional ◽  
Induced Polarization ◽  
The Body ◽  
Scale Model ◽  
Two Dimensional ◽  
Electromagnetic Modeling ◽  
The Earth

The induced polarization (IP) and electromagnetic (EM) responses of a three‐dimensional body in the earth can be calculated using an integral equation solution. The problem is formulated by replacing the body by a volume of polarization or scattering current. The integral equation is reduced to a matrix equation, which is solved numerically for the electric field in the body. Then the electric and magnetic fields outside the inhomogeneity can be found by integrating the appropriate dyadic Green’s functions over the scattering current. Because half‐space Green’s functions are used, it is only necessary to solve for scattering currents in the body—not throughout the earth. Numerical results for a number of practical cases show, for example, that for moderate conductivity contrasts the dipole‐dipole IP response of a body five units in strike length approximates that of a two‐dimensional body. Moving an IP line off the center of a body produces an effect similar to that of increasing the depth. IP response varies significantly with conductivity contrast; the peak response occurs at higher contrasts for two‐dimensional bodies than for bodies of limited length. Very conductive bodies can produce negative IP response due to EM induction. An electrically polarizable body produces a small magnetic field, so that it is possible to measure IP with a sensitive magnetometer. Calculations show that horizontal loop EM response is enhanced when the background resistivity in the earth is reduced, thus confirming scale model results.

The Transition Matrix Formalism for the Scattering of an Alluvium on an Elastic Half-Space

2007 ◽  
Author(s):  
Wen-I Liao ◽  
Tsung-Jen Teng ◽  
Shiang-Jung Wang
Keyword(s):  
Half Space ◽  
Transition Matrix ◽  
Plane Waves ◽  
Three Dimensional ◽  
Scattering Matrix ◽  
Elastic Half Space ◽  
Basis Functions ◽  
Matrix Formalism ◽  
T Matrix ◽  
Singular Wave

This paper develops the transition matrix formalism for scattering from an three-dimensional alluvium on an elastic half-space. Betti’s third identity is employed to establish orthogonality conditions among basis functions that are Lamb’s singular wave functions. The total displacements and associated tractions exterior and interior to the surface are expanded in a Rayleigh series. The boundary conditions are applied and the T-matrix is derived. A linear transformation is utilized to construct a set of orthogonal basis functions. The transformed T-matrix is related to the scattering matrix and it is shown that the scattering matrix is symmetric and unitary and that the T-matrix is symmetric. Typical numerical results obtained by incident plane waves for verification are presented.

The Influnce of Velocity for a Moving Load on the Vibration Isolation Using Pile Group

2012 ◽  
Vol 204-208 ◽  
pp. 210-214
Author(s):  
Man Qing Xu ◽  
Bin Xu
Keyword(s):  
Integral Equation ◽  
Fundamental Solution ◽  
Frequency Domain ◽  
Half Space ◽  
Vibration Isolation ◽  
Moving Load ◽  
Free Field ◽  
Isolation Effect ◽  
Moving Loads ◽  
Patch Load

Based on Biot’s theory and integral transform method, the velocity of moving loads impact on the vibration isolation effect of pile rows embedded in a poroelastic half space is investigated in this study. The free field solution for a moving load applied on the surface of a poroelastic half space and the fundamental solution for a harmonic circular patch load applied in the poroelastic half space are derived first. Using Muki’s method and the fundamental solution for the circular patch load as well as the obtained free field solution for the moving load, the second kind of Fredholm integral equation in the frequency domain describing the dynamic interaction between pile rows and the poroelastic half space is developed. Numerical solution of the frequency domain integral equation and numerical inversion of the Fourier transform yield the time domain response of the pile-soil system. Numerical results of this study show that the same pile rows can achieve a better vibration isolation effect for the lower load speed than for the higher speed. Also, the optimal length of piles for higher speed moving loads is shorter than that for lower speed moving loads.

The behavior of the apparent resistivity phase spectrum in the case of two polarizable media

Geophysics ◽  
1985 ◽  
Vol 50 (5) ◽  
pp. 810-819 ◽  
Author(s):  
Heikki Soininen
Keyword(s):  
Frequency Dependence ◽  
Half Space ◽  
Time Constant ◽  
Apparent Resistivity ◽  
Dispersion Model ◽  
Phase Spectrum ◽  
The Body ◽  
Spectral Parameters ◽  
The Difference ◽  
True Time

I employed numerical modeling to examine the formation of the apparent resistivity phase spectrum first of a polarizable prism situated in a polarizable half‐space, and second of two polarizable prisms joined in an unpolarizable half‐space. The calculations were done using the integral equation technique. The frequency dependence of the resistivity of the polarizable medium is depicted by means of the Cole‐Cole dispersion model. The effect of a weakly polarizable half‐space may be handled by simply adding the phase angle of the half‐space to the apparent phase due to the body. The apparent spectral parameters can be inverted by fitting the sum of two Cole‐Cole dispersion model phase spectra to the apparent phase spectrum. Of the parameters describing the prism, the apparent chargeability is smaller than the chargeability of the original petrophysical spectrum because of geometric attenuation. The apparent frequency dependence, on the other hand, is very close to the value of the original frequency dependence. The apparent time constant is commonly also near the true time constant of the petrophysical spectrum. The values of the apparent spectral parameters of the polarizable half‐space are all close to their petrophysical or true values. The apparent spectrum of two polarizable prisms builds up in a complex fashion. Nevertheless, by measuring the spectra at a number of points along a profile crossing over two formations differing in time constant, the various components can be discriminated from the apparent spectrum even if the difference in time constant is small. As the conductivity contrast decreases, the share of the spectrum of the formation in the apparent spectrum increases. Similarly, the formation with the smaller time constant is in a more advantageous position than the body with the greater time constant.

Numerical Study of Forward-Speed Ship Motion and Added Resistance

2017 ◽  
Author(s):  
D. C. Hong ◽  
T. B. Ha ◽  
K. H. Song
Keyword(s):  
Integral Equation ◽  
Free Surface ◽  
Numerical Study ◽  
Three Dimensional ◽  
The Body ◽  
Experimental Results ◽  
Surface Boundary ◽  
Forward Speed ◽  
Added Resistance ◽  
Following Seas

The added resistance of a ship was calculated using Maruo’s formula [1] involving the three-dimensional Kochin function obtained using the source and normal doublet distribution over the wetted surface of the ship. The density of the doublet distribution was obtained as the solution of the three-dimensional frequency-domain forward-speed Green integral equation containing the exact line integral along the waterline. Numerical results of the Wigley ship models II and III in head seas, obtained by making use of the inner-collocation 9-node second-order boundary element method have been compared with the experimental results reported by Journée [2]. The forward-speed hydrodynamic coefficients of the Wigley models have shown no irregular-frequencylike behavior. The steady disturbance potential due to the constant forward speed of the ship has also been calculated using the Green integral equation associated with the steady forward-speed free-surface Green function since the so-called mj-terms [3] appearing in the body boundary conditions contain the first and second derivatives of the steady potential over the wetted surface of the ship. However, the free-surface boundary condition was kept linear in the present study. The added resistances of the Wigley II and III models in head seas obtained using Maruo’s formula showing acceptable comparison with experimental results, have been presented. The added resistances in following seas obtained using Maruo’s formula have also been presented.

scholarly journals P-FFT Accelerated EFIE with a Novel Diagonal Perturbed ILUT Preconditioner for Electromagnetic Scattering by Conducting Objects in Half Space

2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Lan-Wei Guo ◽  
Yongpin Chen ◽  
Jun Hu ◽  
Joshua Le-Wei Li
Keyword(s):  
Integral Equation ◽  
Electric Field ◽  
Half Space ◽  
Electromagnetic Scattering ◽  
High Performance ◽  
Three Dimensional ◽  
Scattering Problems ◽  
Geometrical Structures ◽  
Incomplete Lu Factorization ◽  
Field Integral

A highly efficient and robust scheme is proposed for analyzing electromagnetic scattering from electrically large arbitrary shaped conductors in a half space. This scheme is based on the electric field integral equation (EFIE) with a half-space Green’s function. The precorrected fast Fourier transform (p-FFT) is first extended to a half space for general three-dimensional scattering problems. A novel enhanced dual threshold incomplete LU factorization (ILUT) is then constructed as an effective preconditioner to improve the convergence of the half-space EFIE. Inspired by the idea of the improved electric field integral operator (IEFIO), the geometrical-optics current/principle value term of the magnetic field integral equation is used as a physical perturbation to stabilize the traditional ILUT perconditioning matrix. The high accuracy of EFIE is maintained, yet good calculating efficiency comparable to the combined field integral equation (CFIE) can be achieved. Furthermore, this approach can be applied to arbitrary geometrical structures including open surfaces and requires no extra types of Sommerfeld integrals needed in the half-space CFIE. Numerical examples are presented to demonstrate the high performance of the proposed solver among several other approaches in typical half-space problems.

Sign in / Sign up

Export Citation Format

Share Document