The transient EM step response of a dipping plate in a conductive half‐space

Geophysics ◽  
1992 ◽  
Vol 57 (9) ◽  
pp. 1116-1126 ◽  
Author(s):  
James E. Hanneson

An algorithm for computing the transient electromagnetic (TEM) response of a dipping plate in a conductive half‐space has been developed. For a stationary [Formula: see text] current loop source, calculated profiles simulate the response of the University of Toronto electromagnetic system (UTEM) over a plate in a 1000 Ω ⋅ m half‐space. The objective is to add to knowledge of the galvanic process (causing poloidal plate currents) and the local induction process (causing toroidal currents) by studying host and plate currents with respect to surface profiles. Both processes can occur during TEM surveys. Plates are all [Formula: see text] thick with various depths, dips, and conductances. Calculated host and plate currents provide quantitative examples of several effects. For sufficiently conductive plates, the late time currents are toroidal as for a free‐space host. At earlier times, or at all times for poorly conducting plates, the plate currents are poloidal, and the transitions to toroidal currents, if they occur, are gradual. At very late times, poloidal currents again dominate any toroidal currents but this effect is rarely observed. Stripped, point‐normalized profiles, which reflect secondary fields caused by the anomalous plate currents, illustrate effects such as early time blanking (caused by noninstantaneous diffusion of fields into the target), mid‐time anomaly enhancement (caused by galvanic currents), and late time plate‐in‐free‐space asymptotic behavior.

Geophysics ◽  
1985 ◽  
Vol 50 (6) ◽  
pp. 1002-1019 ◽  
Author(s):  
Peggie R. Gallagher ◽  
Stanley H. Ward ◽  
G. W. Hohmann

The computer program PLATE, developed at the University of Toronto, models the electromagnetic (EM) response of an inductively thin plate in free space. We used PLATE to compute two components of the time derivative of the magnetic field for a range of models for the EM37 fixed‐source transient system ([Formula: see text] loop). Analysis of the response curves produced methods of interpretation for obtaining plate geometry and conductance. The overall width of an anomaly, the distance between peaks and the width of the updip lobes, can provide an estimate of depth. Dip has the dominant effect on the ratio of the peak amplitudes. A rough estimate of plate size and the position in time (early or late) of the currents is essential before proceeding with interpretation. Strike length is not obviously reflected in the shape of the curves, but depth extent is indicated by the rate at which the downdip tail returns to the baseline, except for vertical plates. For vertical plates, curve matching may be the only method of obtaining an estimate of depth extent. Varying conductance for a particular model in free space affects whether a channel represents an early, intermediate, or late time response. The shape of a profile varies with the time of measurement. The estimated time constant can be used to calculate the conductance, provided an estimate of the shortest dimension of the plate is available. Extinction angles appear frequently for plates of small depth extent but do not occur for plates which are of infinite strike and depth extent with respect to the size of the transmitting loop.


Geophysics ◽  
2001 ◽  
Vol 66 (4) ◽  
pp. 1090-1096 ◽  
Author(s):  
Z. Zhang ◽  
J. Xiao

Large‐loop electromagnetic (EM) systems that measure transient EM (TEM) data on the surface or in boreholes have shown increased application in exploration geophysics. Accurate interpretation of borehole TEM data is necessary to discover deep hidden targets that cannot be detected with surface systems. However, the inversion of borehole TEM data has not been fully addressed. In this paper, we study the propagation of the TEM field from a large‐loop EM borehole system inside a layered earth and develop a new inversion algorithm to reconstruct layered conductivity structures from large‐loop TEM data measured with both surface and borehole configurations. The magnetic field and sensitivities are first computed in the frequency domain and are then transformed into the time domain where the inversion is performed. The surface data have a higher S/N ratio at early time channels, while the borehole data have a higher S/N ratio at late time channels. Consequently, the surface data can be inverted to better resolve shallow structures, and the borehole data can be used to better detect deep structures. The merits of joint inversions of borehole and surface data are explored. We test our inversion algorithm using numeric examples.


2019 ◽  
Vol 24 (4) ◽  
pp. 653-663 ◽  
Author(s):  
Xin Wu ◽  
Guangyou Fang ◽  
Guoqiang Xue ◽  
Lihua Liu ◽  
Leisong Liu ◽  
...  

Over the past decade, helicopter-borne transient electromagnetic (HTEM) systems have been rapidly developed. A new HTEM prototype (referred to as a CAS-HTEM) has been developed by the Chinese Academy of Sciences. In terms of hardware, the CAS-HTEM system uses an inflatable structure to carry the transmitting loop, which significantly reduces the weight of the system and makes it easier to transport. A dual gain receiver was innovated to extend the dynamic range of the system. In addition, an observation circuit for transmitting voltage waveform is introduced, so that the derivative waveform of transmitting current with higher SNR could be calculated. In terms of data processing, more reliable early time data could be obtained by band-limited effect removal; a weighted stacking algorithm is introduced to reduce the narrow band noise more effectively and increase the sensitivity of data to the anomaly location; a method based on τ-domain transform is used for late time signal processing. The results of the field test which was carried out in Inner Mongolia were found to be consistent with the drill data, which effectively verified the performance of this HTEM prototype.


Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1462-1471 ◽  
Author(s):  
Brian R. Spies ◽  
Dwight E. Eggers

Problems and misunderstandings arise with the concept of apparent resistivity when the analogy between an apparent resistivity computed from geophysical observations and the true resistivity structure of the subsurface is drawn too tightly. Several definitions of apparent resistivity are available for use in electromagnetic methods; however, those most commonly used do not always exhibit the best behavior. Many of the features of the apparent resistivity curve which have been interpreted as physically significant with one definition disappear when alternative definitions are used. It is misleading to compare the detection or resolution capabilities of different field systems or configurations solely on the basis of the apparent resistivity curve. For the in‐loop transient electromagnetic (TEM) method, apparent resistivity computed from the magnetic field response displays much better behavior than that computed from the induced voltage response. A comparison of “exact” and “asymptotic” formulas for the TEM method reveals that automated schemes for distinguishing early‐time and late‐time branches are at best tenuous, and those schemes are doomed to failure for a certain class of resistivity structures (e.g., the loop size is large compared to the layer thickness). For the magnetotelluric (MT) method, apparent resistivity curves defined from the real part of the impedance exhibit much better behavior than curves based on the conventional definition that uses the magnitude of the impedance. Results of using this new definition have characteristics similar to apparent resistivity obtained from time‐domain processing.


Geophysics ◽  
1984 ◽  
Vol 49 (7) ◽  
pp. 895-901 ◽  
Author(s):  
W. E. Doll ◽  
T. V. Skibicky ◽  
C. S. Clay

Transient electromagnetic (TEM) soundings can be enhanced by additions of matched filtered signals from different coil separations. Results of laboratory metallic model studies are used to demonstrate this technique. A short pulse of current is transmitted from one coil. For each of the signals received in the other coil, the half‐space response (surface wave) is subtracted from the total signal to determine the reflected signal. These reflected signals are passed through matched filters corresponding to the proper separation and a possible layer depth. All of the matched filtered signals corresponding to a particular possible depth are added, and the summed output is stored. The filtering and summing process is repeated for a set of possible interface depths. The sum of matched filtered signals for the correct depth has a peak value of 1 and is symmetric about the peak. We believe that short‐range early‐time values of the transient response can be used to estimate the conductivity of the first layer. This estimate is used for numerical computations of half‐space responses. We have used several magnetic dipoles to represent the finite‐sized coils to show good agreement between theoretical and experimental half‐spaces. Theoretical reflection signals for matched filters can be computed numerically. The stacking of matched filtered signals is an electromagnetic analog of the common‐depth‐point (CDP) seismic technique. It should reduce the effects of nearsurface inhomogeneities and improve the resolution. As an incidental experiment, we show comparisons of the reflected signal and the “image transmission.” The reflected and transmitted signals are essentially the same and support the concept of hybrid‐ray theory for electromagnetic sounding.


Geophysics ◽  
1986 ◽  
Vol 51 (6) ◽  
pp. 1291-1297 ◽  
Author(s):  
Yang Sheng

Early‐time and late‐time apparent resistivity approximations have been widely used for interpretation of long‐offset transient electromagnetic (LOTEM) measurements because it is difficult to find a single apparent resistivity over the whole time range. From a physical point of view, Dr. C. H. Stoyer defined an apparent resistivity for the whole time range. However, there are two problems which hinder its use: one is that there is no explicit formula to calculate the apparent resistivity, and the other is that the apparent resistivity has no single solution. A careful study of the two problems shows that a numerical method can be used to calculate a single apparent resistivity. A formula for the maximum receiver voltage over a uniform earth, when compared with the receiver voltage for a layered earth, leads to the conclusion that, in some cases, a layered earth can produce a larger voltage than any uniform earth can produce. Therefore, our apparent resistivity definition cannot be applied to those cases. In some other cases, the two possible solutions from our definition do not merge, so that neither of them is meaningful for the whole time range.


Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 49-55 ◽  
Author(s):  
Kenneth Duckworth ◽  
Edward S. Krebes

The concept of electromagnetic depth sounding by means of a coincident‐coil frequency‐domain electromagnetic system is developed in theory and demonstrated by means of physical scale modeling. The concept is based on the use of distance from the target as the sounding variable. The theoretical developments are confined to soundings conducted in free‐space with respect to either a homogeneous half‐space or a thin sheet conductor in conditions that approach the resistive limit. The use of distance from the target as the sounding variable becomes practical when the sounding system is a single compact unit of the type that a coincident coil concept inherently provides. In this method of sounding, the distance from the target is determined by taking the ratios of the fields measured at a variety of distances from the target conductor. This permits not only the distance to the target to be determined but also the direction to that target as may be of interest in soundings conducted in mines.


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 ◽  
1984 ◽  
Vol 49 (7) ◽  
pp. 861-869 ◽  
Author(s):  
Perry A. Eaton ◽  
Gerald W. Hohmann

We have computed transient borehole electromagnetic (EM) responses of two‐dimensional (2-D) models using a direct and explicit finite‐difference algorithm. The program computes the secondary electric field which is defined as the difference between the total field and the primary (half‐space) field. The time derivative of the vertical magnetic field in a borehole is computed by numerical differentiation of the total electric field. These models consist of a thin horizontal conductor with a finite width, embedded in a conductive half‐space. Dual line sources energized by a step‐function current lie on the surface of the half‐space and simulate the long sides of a large rectangular loop. Numerical results substantiate several important features of the transient impulse response of such models. The peak response of the target is attenuated as the resistivity of the host decreases. A sign reversal in the secondary electric field occurs later in time as the resistivity of the host decreases. The peak response and the onset of late‐time behavior are delayed in time as well. Secondary responses for models with different host resistivities (10–1000 Ω-m) are approximately the same at late time. If the target is less conductive, the effects of the host, i.e., the attenuation and time delay, are less. It is readily apparent that there exists a time window within which the target’s response is at a maximum relative to the half‐space response. At late time the shape of the borehole anomaly due to a thin conductive 2-D target appears to be independent of the conductivity of the host. The late‐time secondary decay of the target is neither exponential nor power law, and a time constant computed from the slope of a log‐linear decay curve at late time may be much larger than the actual value for the same target in free space.


Sign in / Sign up

Export Citation Format

Share Document