Three‐dimensional transient electromagnetic responses for a grounded source

Geophysics ◽  
1986 ◽  
Vol 51 (11) ◽  
pp. 2117-2130 ◽  
Author(s):  
Brian M. Gunderson ◽  
Gregory A. Newman ◽  
Gerald W. Hohmann

When the current in a grounded wire is terminated abruptly, currents immediately flow in the Earth to preserve the magnetic field. Initially the current is concentrated near the wire, with a broad zone of return currents below. The electric field maximum broadens and moves downward with time. Currents are channeled into a conductive three‐dimensional body, resulting in anomalous magnetic fields. At early times, when the return currents are channeled into the body, the vertical magnetic field is less than the half‐space field on the far side of the body but is greater than the half‐space field between the source and the body. Later the current in the body reverses; the vertical field is enhanced on the far side of the body and decreased between the source and the body. The horizontal magnetic field has a well‐defined maximum directly over the body at late times, and is a better indicator of the position of the body. The vertical magnetic field and its time derivative change sign with time at receiver locations near the source if a three‐dimensional body is present. These sign reversals present serious problems for one‐dimensional inversion, because decay curves for a layered earth do not change sign. At positions away from the source, the decay curves exhibit no sign reversals—only decreases and enhancements relative to one‐dimensional decay curves. In such cases one‐dimensional inversions may provide useful information, but they are likely to result in fictitious layers and erroneous interpretations.

Geophysics ◽  
1989 ◽  
Vol 54 (5) ◽  
pp. 598-608 ◽  
Author(s):  
Gregory A. Newman ◽  
Walter L. Anderson ◽  
Gerald W. Hohmann

Transient electromagnetic (TEM) borehole responses of 3-D vertical and horizontal tabular bodies in a half‐space are calculated to assess the effect of a conductive host. The transmitter is a large loop at the surface of the earth, and the receiver measures the time derivative of the vertical magnetic field. When the host is conductive (100 Ω ⋅ m), the borehole response is due mainly to current channeled through the body. The observed magnetic‐field response can be visualized as due to galvanic currents that pass through the conductor and return in the half‐space. When the host resistivity is increased, the magnetic field of the conductor is influenced more by vortex currents that flow in closed loops inside the conductor. For a moderately resistive host (1000 Ω ⋅ m), the magnetic field of the body is caused by both vortex and galvanic currents. The galvanic response is observed at early times, followed by the vortex response at later times if the body is well coupled to the transmitter. If the host is very resistive, the galvanic response vanishes; and the response of the conductor is caused only by vortex currents. The shapes of the borehole profiles change considerably with changes in the host resistivity because vortex and galvanic current distributions are very different. When only the vortex response is observed, it is easy to distinguish vertical and horizontal conductors. However, in a conductive host where the galvanic response is dominant, it is difficult to interpret the geometry of the body; only the approximate location of the body can be determined easily. For a horizontal conductor and a single transmitting loop, only the galvanic response enables one to determine whether the conductor is between the transmitter and the borehole or beyond the borehole. A field example shows behavior similar to that of our theoretical results.


Geophysics ◽  
1971 ◽  
Vol 36 (1) ◽  
pp. 101-131 ◽  
Author(s):  
Gerald W. Hohmann

A theoretical solution is developed for the electromagnetic response of a two‐dimensional inhomogeneity in a conductive half‐space, in the field of a line source of current. The solution is in the form of an integral equation, which is reduced to a matrix equation, and solved numerically for the electric field in the body. The electric and magnetic fields at the surface of the half‐space are found by integrating the half‐space Green’s functions over the scattering currents. One advantage of this particular numerical technique is that it is necessary to solve for scattering currents only in the conductor and not throughout the half‐space. The response of a thin, vertical conductor is studied in some detail. Because the only interpretational aids available previously were scale model results for conductors in free space, the results presented here should be useful in interpreting data and in designing new EM systems. As expected, anomalies decay rapidly as depth of burial is increased, due to attenuation in the conductive half‐space. Depth of exploration appears to be greatest for measurements of horizontal magnetic field phase, while vertical field phase is diagnostic of conductivity. Horizontal location and depth of burial are best determined through measurements of vertical or horizontal magnetic field amplitude.


Geophysics ◽  
1985 ◽  
Vol 50 (7) ◽  
pp. 1144-1162 ◽  
Author(s):  
William A. SanFilipo ◽  
Perry A. Eaton ◽  
Gerald W. Hohmann

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.


Geophysics ◽  
1985 ◽  
Vol 50 (12) ◽  
pp. 2849-2861 ◽  
Author(s):  
Jopie I. Adhidjaja ◽  
Gerald W. Hohmann ◽  
Michael L. Oristaglio

The time‐domain electromagnetic (TEM) modeling method of Oristaglio and Hohmann is reformulated here in terms of the secondary field. This finite‐difference method gives a direct, explicit time‐domain solution for a two‐dimensional body in a conductive earth by advancing the field in time with DuFort‐Frankel time‐differencing. As a result, solving for the secondary field, defined as the difference between the total field and field of a half‐space, is not only more efficient but is also simpler and eliminates several problems inherent in the solution for the total field. For example, because the secondary field varies slowly both in space and time, it can be modeled on a coarse grid with large time steps. In addition, for a simple body the field is local; therefore, because the field can be assumed to satisfy a simple boundary condition in the earth computation is greatly simplified. Our tests show that for the same accuracy, the secondary‐field solution is roughly five times faster than the total‐field solution. We compute and analyze the magnetic field impulse response for a suite of models, most of which consist of a thin body embedded in a conductive half‐space—with or without overburden. The results indicate the conductive half‐space will both delay and attenuate the response of the body and even obscure it if the conductivity contrast is small. The results also suggest that the conductive host can alter the decay rate of the response of the body from its free‐space counterpart. Our results for multiple bodies illustrate the importance of early‐time measurements to obtain resolution, particularly for measurements of the horizontal magnetic field. The vertical magnetic field, however, can be used to infer the dip direction of a dipping body by studying the migration of the crossover. The results for models which include overburden show that the effect of a conductive overburden, in addition to the half‐space effect, is to delay the response of the body, because the primary current initially tends to concentrate and slowly diffuse through the overburden, and does not reach the body until later time. This effect also complicates the early‐times profiles, becoming more severe as the conductivity of the overburden is increased.


2021 ◽  
Vol 87 (2) ◽  
Author(s):  
P. Kim ◽  
R. Jorge ◽  
W. Dorland

A simplified analytical form of the on-axis magnetic well and Mercier's criterion for interchange instabilities for arbitrary three-dimensional magnetic field geometries is derived. For this purpose, a near-axis expansion based on a direct coordinate approach is used by expressing the toroidal magnetic flux in terms of powers of the radial distance to the magnetic axis. For the first time, the magnetic well and Mercier's criterion are then written as a one-dimensional integral with respect to the axis arclength. When compared with the original work of Mercier, the derivation here is presented using modern notation and in a more streamlined manner that highlights essential steps. Finally, these expressions are verified numerically using several quasisymmetric and non-quasisymmetric stellarator configurations including Wendelstein 7-X.


2008 ◽  
Vol 8 (3) ◽  
pp. 501-507 ◽  
Author(s):  
G. Prattes ◽  
K. Schwingenschuh ◽  
H. U. Eichelberger ◽  
W. Magnes ◽  
M. Boudjada ◽  
...  

Abstract. We present the results of ground-based Ultra Low Frequency (ULF) magnetic field measurements observed from June to August 2004 during the Bovec earthquake on 12 July 2004. Further we give information about the seismic activity in the local observatory region for an extended time span 2004 and 2005. ULF magnetic field data are provided by the South European Geomagnetic Array (SEGMA) where the experience and heritage from the CHInese MAGnetometer (CHIMAG) fluxgate magnetometer comes to application. The intensities of the horizontal H and vertical Z magnetic field and the polarization ratio R of the vertical and horizontal magnetic field intensity are analyzed taking into consideration three SEGMA observatories located at different close distances and directions from the earthquake epicenter. We observed a significant increase of high polarization ratios during strong seismic activity at the observatory nearest to the Bovec earthquake epicenter. Apart from indirect ionospheric effects electromagnetic noise could be emitted in the lithosphere due to tectonic effects in the earthquake focus region causing anomalies of the vertical magnetic field intensity. Assuming that the measured vertical magnetic field intensities are of lithospheric origin, we roughly estimate the amplitude of electromagnetic noise in the Earths crust considering an average electrical conductivity of <σ>=10−3 S/m and a certain distance of the observatory to the earthquake epicenter.


Geophysics ◽  
2004 ◽  
Vol 69 (5) ◽  
pp. 1192-1202 ◽  
Author(s):  
Michael Commer ◽  
Gregory Newman

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.


2018 ◽  
Vol 32 (32) ◽  
pp. 1850390
Author(s):  
Minos A. Neto ◽  
J. Roberto Viana ◽  
Octavio D. R. Salmon ◽  
E. Bublitz Filho ◽  
José Ricardo de Sousa

The critical frontier of the isotropic antiferromagnetic Heisenberg model in a magnetic field along the z-axis has been studied by mean-field and effective-field renormalization group calculations. These methods, abbreviated as MFRG and EFRG, are based on the comparison of two clusters of different sizes, each of them trying to mimic a specific Bravais lattice. The frontier line in the plane of temperature versus magnetic field was obtained for the simple cubic and the body-centered cubic lattices. Spin clusters with sizes N = 1, 2, 4 were used so as to implement MFRG-12, EFRG-12 and EFRG-24 numerical equations. For the simple cubic lattice, the MFRG frontier exhibits a notorious re-entrant behavior. This problem is improved by the EFRG technique. However, both methods agree at lower fields. For the body-centered cubic lattice, the MFRG method did not work. As in the cubic lattice, all the EFRG results agree at lower fields. Nevertheless, the EFRG-12 approach gave no solution for very low temperatures. Comparisons with other methods have been discussed.


Sign in / Sign up

Export Citation Format

Share Document