• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.184-190
• /
• 2005

We have attempted to give some physical into the factors that control the response of subsurface target to plane wave excitation at VLF (very low frequency) frequencies. Although the VLF technique has some serious disadvantages, such as a lack of penetration depth and limited frequencies, its advantages are also extremely important to have made it by far the most popular electromagnetic technique in current use. In the magnetic-field measurement mode these lightweight, relatively low-priced tools allow us to survey large areas rapidly and inexpensively, to locate and roughly define subsurface electrical conductors. When used in the electric-field mode the technique is, in simple environments, capable of quantitative interpretation and once again the speed and relative cheapness of these devices makes them a natural 'first electromagnetic tool ' to use in reconnaissance mapping.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.17-25
• /
• 1999

A traveltime tomography has been carried out by transforming electromagnetic data in frequency domain to wave-like domain. The transform uniquely relates a field satisfying a diffusion equation to an integral of the corresponding wavefield. But direct transform of frequency domain magnetic fields to wave-field domain is ill-posed problem because the kernel of the integral transform is highly damped. In this study, instead of solving such an unstable problem, it is assumed that wave-fields in transformed domain can be approximated by sum of ray series. And for further simplicity, reflection and refraction energy compared to that of direct wave is weak enough to be neglected. Then first arrival can be approximated by calculating the traveltime of direct wave only. But these assumptions are valid when the conductivity contrast between background medium and the target anomalous body is low enough. So this approach can only be applied to the models with low conductivity contrast. To verify the algorithm, traveltime calculated by this approach was compared to that of direct transform method and exact traveltime, calculated analytically, for homogeneous whole space. The error in first arrival picked by this study was less than that of direct transformation method, especially when the number of frequency samples is less than 10, or when the data are noisy. Layered earth model with varying conductivity contrasts and inclined dyke model have been successfully imaged by applying nonlinear traveltime tomography in 30 iterations within three CPU minutes on a IBM Pentium Pro 200 MHz.

• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.150-161
• /
• 2018

Recently, unmanned aircraft EM (electromagnetic) survey based on ICT (Information and Communication Technology) has been widely utilized because of the efficiency in regional survey. We performed the theoretical study on the unmanned airship EM system developed by KIGAM (Korea Institute of Geoscience and Mineral resources) as part of the practical application of unmanned aircraft EM survey. Since this system has different configurations of transmitting and receiving loops compared to the conventional aircraft EM systems, a new technique is required for the appropriate interpretation of measured responses. Therefore, we proposed a method to calculate the EM field for the arbitrary shaped transmitter and verified its validity through the comparison with analytic solution for circular loop. In addition, to simulate the magnetic responses by three-dimensionally (3D) distributed anomalies, we have adapted our algorithm to 3D frequency-domain EM modeling algorithm based on the edge-FEM (finite element method). Though the analysis on magnetic field responses from a subsurface anomaly, it was found that the response decreases as the depth of the anomaly increases or the flight altitude increases. Also, it was confirmed that the response became smaller as the resistivity of the anomaly increases. However, a nonlinear trend of the out-of-phase component is shown depending on the depth of the anomaly and the used frequency, that makes it difficult to apply simple analysis based on the mapping of the magnitude of the responses and can cause the non-uniqueness problem in calculating the apparent resistivity. Thus, it is a prerequisite to analyze the appropriate frequency band and flight altitude considering the purpose of the survey and the site conditions when conducting a survey using the unmanned aircraft EM system.

• Geophysics and Geophysical Exploration
• /
• v.5 no.2
• /
• pp.84-92
• /
• 2002

In this study, the new modeling scheme has been developed for recently designed and tested electromagnetic survey, which adapts horizontal magnetic dipole with $1\;kHz\~1\;MHz$ frequency range as a source. The 2.5-D secondary field formulation in wavenumber domain was constructed using finite element method and verified through comparing results with layered-earth solutions calculated by integral equations. 2-D conductive- and resistive-block models were constructed for calculating electric field, magnetic field and impedance - the ratio of electric and magnetic fields which are orthogonal each other. This study showed that electric field and impedance are superior in identifying 2-D isolated-body model to magnetic field. In particular, impedance gives more stable results than electric field with similar spatial resolving power, because electric field is divided by magnetic field in impedance. Thus the impedance analysis which uses electric and magnetic fields together would give better result in imaging the shallow anomalies than conventional EM method.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.265-270
• /
• 2006

High-frequency electromagnetic (EM) wave propagation phenomena associated with borehole ground-penetrating radar (GPR) surveys are complex. To improve the understanding of governing physical processes, we present a finite-difference time-domain solution of Maxwell's equations in cylindrical coordinates. This approach allows us to model the full EM wavefield associated with borehole GPR surveys. The algorithm can be easily implemented perfectly matched layers for absorbing boundaries, frequency-dependent media, and finite-length transmitter antenna.

• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.44 no.11
• /
• pp.42-47
• /
• 2007

In this paper, influence of lossy ground and gap variation between lossy ground and UWB antenna on impulse propagation in time domain for impulse ground penetrating radar (GPR) is numerically and experimentally investigated. For this study, a novel planar UWB fat dipole antenna is developed. First, influence of lossy ground and gap variation between lossy ground and UWB antenna is simulated. For verification, a test field of sand and wet clay soil is built and using the developed dipole antenna, transmission behavior is investigated at the test field. With an aid of IDFT (inverse discrete Fourier transform), time domain impulse response for transmission coefficient measured and simulated in frequency domain is obtained. Measurement and simulation show that the frequency of maximum transmission coefficient and transmission coefficient are increased with higher dielectric constant and larger gap distance. In time domain, it is shown that for higher dielectric constant, the amplitude of the received signal in time domain is higher and reflected signals are seriously modified. Also, it is found that variation of gap between antenna and ground surface makes timing of peak value changed.

• Geophysics and Geophysical Exploration
• /
• v.6 no.3
• /
• pp.143-152
• /
• 2003

Among the various geophysical approaches to identify the leakage of leachate with conductivity variation, conventional electrical resistivity survey has been mainly used at waste landfill. We adopted small-loop electromagnetic (EM) survey using multi-frequencies in parallel with electrical resistivity survey to delineate the leakage of leachate through the shallow soil layer at a waste landfill in Jeju Island, and also with self-potential monitoring to detect the streaming potential produced by the movement of leachate. There were no evidences of leakage from waste landfill according to the results of the electrical resistivity survey and SP monitoring, and it was also true from the results of water quality analysis at stream around waste landfill periodically. On the other hand, the results of one-dimensional inversion of spatially-filtered small-loop EM survey data showed the anomalous zone of low resistivity with depth both around and inner waste landfill.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Geophysics and Geophysical Exploration
• /
• v.2 no.2
• /
• pp.96-103
• /
• 1999

Most of saline water inousions have been diagnosed by geophysical or geochemical approach independently. The objective of this study is to provide the effective method to detect the saline water intrusion on the ground water in the vicinity of seashore using these two methods. Schulumberger sounding, frequency domain electromagnetic sounding and geochemical analysis of ground water were carried out to explore saline water intrusion. Schulumberger sounding was implemented in dry surface condition before irrigation water was introduced into the field, while electromagnetic sounding was carried out in wet ground condition after the irrigation. The purpose of duplicated measurements on the equivalent spot at different times was to investigate the variation of anomaly zone depending on the amount of ground water. It was possible to discriminate the anomalous zone due to high water saturation from the low electric resistivity zone by high salt concentrations through this way. For the verification of the geophysical result, the ground water samples in the study area were collected and analysed at the 23 points near the measuring spots. The groundwater at the spot nearest to the sea water intrusion identified by geophysical method indicates higher salinity than the standard limit concentration for agricultural irrigation water (250 mg/1). Isotope analysis of $D({^2}H)$ vs. is ${^18}O$ and PCA analysis were used to discriminate the anthropogenic pollution from those of high salinity from sea water intrusion.

• Geophysics and Geophysical Exploration
• /
• v.1 no.2
• /
• pp.127-135
• /
• 1998

The Born approximation is widely used for solving the complex scattering problems in electromagnetics. Approximating total internal electric field by the background field is reasonable for small material contrasts as long as scatterer is not too large and the frequency is not too high. However in many geophysical applications, moderate and high conductivity contrasts cause both real and imaginary part of internal electric field to differ greatly from background. In the extended Born approximation, which can improve the accuracy of Born approximation dramatically, the total electric field in the integral over the scattering volume is approximated by the background electric field projected to a depolarization tensor. The finite difference and elements methods are usually used in EM scattering problems with a 2D model and a 3D source, due to their capability for simulating complex subsurface conductivity distributions. The price paid for a 3D source is that many wavenumber domain solutions and their inverse Fourier transform must be computed. In these differential equation methods, all the area including homogeneous region should be discretized, which increases the number of nodes and matrix size. Therefore, the differential equation methods need a lot of computing time and large memory. In this study, EM modeling program for a 2D model and a 3D source is developed, which is based on the extended Born approximation. The solution is very fast and stable. Using the program, crosshole EM responses with a vertical magnetic dipole source are obtained and the results are compared with those of 3D integral equation solutions. The agreement between the integral equation solution and extended Born approximation is remarkable within the entire frequency range, but degrades with the increase of conductivity contrast between anomalous body and background medium. The extended Born approximation is accurate in the case conductivity contrast is lower than 1:10. Therefore, the location and conductivity of the anomalous body can be estimated effectively by the extended Born approximation although the quantitative estimate of conductivity is difficult for the case conductivity contrast is too high.