• Proceedings of the KSEEG Conference
• /
• 1999.04a
• /
• pp.77-79
• /
• 1999

• Geophysics and Geophysical Exploration
• /
• v.17 no.1
• /
• pp.21-27
• /
• 2014

We analyzed response patterns of test field data acquired with new small loop electromagnetic (EM) system using three-dimensional (3D) electromagnetic modeling code. The size and shape of a conductor was adopted as experimental parameters for EM modeling to understand influencing factors of the response patterns due to a metallic object on the seafloor. Obtaining the responses for four models of difference sizes and shapes through 3D EM modeling, we confirmed that the shape of the object have a more critical factor on the response pattern than size. We also calculated "ppm" values with respect to different altitudes of the sensor and source frequencies. The modeling results show that the consistency of sensor altitude is important and imaginary part of ppm response is more sensitive than real part. We also visualized the contour map of the real and imaginary part of ppm value as a function of frequency and altitude so that we can estimate proper altitude for source frequency band of our survey system. The results of this paper are anticipated to give proper parameters in survey construction for seafloor massive sulfide deposit.

• Journal of the Korean Geophysical Society
• /
• v.5 no.1
• /
• pp.63-71
• /
• 2002

The integral equation method is a powerful tool for electromagnetic numerical modeling. But the difficulty of this technique is the size of their linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. To overcome this limitation, a lot of approximation and series methods, such as conventional Born, modifed Born and extended Born, were developed. But all the methods need volume integration of Green tensor, which is very time consuming. In electromagnetic theory, Green tensor rapidly decreases as the distance between source and field cell increases. Therefore, the source cell which are far away from the field cell does not make an effect on the electric field of the field cell. Consequently, by ignoring the effect of Green tensor due to far away source cells, computing time for electromagnetic numerical modeling can be reduced dramatically. Comparisons of this new method against a full integral equation, extended Born approximation and series code show that the method is accurate enough much less time consuming.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.280-290
• /
• 2002

Three-dimensional (3-D) electromagnetic (EM) modeling algorithm has been developed using finite element method (FEM) to acquire more efficient interpretation techniques of EM data. When FEM based on nodal elements is applied to EM problem, spurious solutions, so called 'vector parasite', are occurred due to the discontinuity of normal electric fields and may lead the completely erroneous results. Among the methods curing the spurious problem, this study adopts vector element of which basis function has the amplitude and direction. To reduce computational cost and required core memory, complex bi-conjugate gradient (CBCG) method is applied to solving complex symmetric matrix of FEM and point Jacobi method is used to accelerate convergence rate. To verify the developed 3-D EM modeling algorithm, its electric and magnetic field for a layered-earth model are compared with those of layered-earth solution. As we expected, the vector based FEM developed in this study does not cause ny vector parasite problem, while conventional nodal based FEM causes lots of errors due to the discontinuity of field variables. For testing the applicability to high frequencies 100 MHz is used as an operating frequency for the layer structure. Modeled fields calculated from developed code are also well matched with the layered-earth ones for a model with dielectric anomaly as well as conductive anomaly. In a vertical electric dipole source case, however, the discontinuity of field variables causes the conventional nodal based FEM to include a lot of errors due to the vector parasite. Even for the case, the vector based FEM gave almost the same results as the layered-earth solution. The magnetic fields induced by a dielectric anomaly at high frequencies show unique behaviors different from those by a conductive anomaly. Since our 3-D EM modeling code can reflect the effect from a dielectric anomaly as well as a conductive anomaly, it may be a groundwork not only to apply high frequency EM method to the field survey but also to analyze the fold data obtained by high frequency EM method.

• Proceedings of the KSEEG Conference
• /
• 2001.04a
• /
• pp.130-133
• /
• 2001

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.361-366
• /
• 2009

Layered-earth Green's functions in electormagnetic (EM) surveys play a key role in modeling the response of exploration targets. They are computed through the Hankel transforms of analytic kernels. Computational precision depends upon the choice of algebraically equivalent forms by which these kemels are expressed. Since three-dimensional (3D) modeling can require a huge number of Green's function evaluations, total computational time can be influenced by computational time for the Hankel transform evaluations. Linear digital filters have proven to be a fast and accurate method of computing these Hankel transforms. In EM modeling for 3D inversion, electric fields are generally evaluated by the secondary field formulation to avoid the singularity problem. In this study, three components of electric fields for five different sources on the surface of homogeneous half-space were derived as primary field solutions. Moreover, reflection coefficients in TE and TM modes were produced to calculate EM responses accurately for a two-layered model having a sea layer. Accurate primary fields should substantially improve accuracy and decrease computation times for Green's function-based problems like MT problems and marine EM surveys.

• Geophysics and Geophysical Exploration
• /
• v.2 no.3
• /
• pp.142-148
• /
• 1999

Three-dimensional electromagnetic modeling algorithm in homogeneous half-space was developed using the extended Born approximation to an electric field integral equation. To examine the performance of the extended Born approximation algorithm, the results were compared with those of the full integral equation results. For a crosshole source-receiver configuration, the agreement between the integral equation and the extended Born approximation was remarkable when the source frequency is lower than 20 kHz and conductivity contrast lower than 1:10. Beyond this conductivity contrast, the simulated results by the extended Born approximation exhibit a difference with respect to those by the integral equation. Therefore, the limit of accuracy lies below contrast of 1:10 in the extended Born approximation. Since for the source frequency range from 20 kHz to 100 kHz, however, the difference is relatively small, the extended Born approximation could be used for a reasonable 3-D EM modeling algorithm.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.66-74
• /
• 2012

This paper presents development of a three-dimensional marine controlled-source electromagnetic (mCSEM) modeling algorithm and its application to a salt and reservoir model to examine detectability of mCSEM for a reservoir under complex subsurface structures. The algorithm is based on the finite difference method, and employs the secondary field formulation for an accurate and fast calculation of modeling responses. The algorithm is verified for a two-layer model by comparing solutions not only with analytic solutions but also with those from other 3D modeling algorithm. We calculate and analyze electric and magnetic fields and their normalized responses for a salt and reservoir model due to three sources located at boundaries between a salt, a reservoir, and background. Numbers and positions of resistive anomalies are informed by normalized responses for three sources, and types of resistive anomalies can be informed when there is a priori information about a salt by seismic exploration.

• Geophysics and Geophysical Exploration
• /
• v.20 no.1
• /
• pp.33-42
• /
• 2017

Recently, the grounded electrical-source airborne transient electromagnetic (GREATEM) system with high power source was introduced to achieve deeper investigation depth and to overcome high noise level. Although the GREATEM is a transient electromagnetic system using a long grounded wire as the transmitter, GREATEM data have been interpreted with 1D earth models because 2D or 3D modeling and inversion of vast airborne data are complicated and expensive to calculate. Generally, 1D inversion is subsequently applied to every survey point and combining 1D images together forms the stitched conductivity-depth image. However, the stitched models often result in abrupt variations in neighboring models. To overcome this problem, laterally constrained inversion (LCI) has been developed in inversion of ATEM data, which can yield layered sections with lateral smooth transitions. In this study, we analysed the GREATEM data through 1D numerical modeling for a curved grounded wire source. Furthermore, we developed a laterally constrained inversion scheme for continuous GREATEM data based on a layered earth model. All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. Applying the developed LCI technique to the GREATEM data, it was confirmed that the laterally constrained inversion can provide laterally smooth model sections that reflect the layering of the survey area effectively.

• Proceedings of the KSEEG Conference
• /
• 1999.04a
• /
• pp.74-76
• /
• 1999