• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.198-200
• /
• 2002

경제적이며 정확한 3차원 전자탐사 모델링을 위해 위해 Habashy et al. (1993)에 의해 제안된 국소 비선형 근사(localized nonlinear approximation)를 이용하여 전자탐사 모델링 알고리듬을 개발하였다. 전자탐사 수치모델링시 많은 계산시간 및 기억용량을 필요로 하는 Green 텐서 적분을 정확하고 빠르게 계산하기 위해, 단일 미소체를 이용한 공간파수 영역에서의 Green 텐서 적분 알고리듬을 제안하였다. 더욱이 Green 텐서의 송수신 방향 및 상반성을 고려하여 각각의 미소체에 의한 전체 미소체에의 Green 텐서 적분을 한 개의 미소체에 의한 전체 미소체에의 Green 텐서 적분 값으로 구하게 하므로 매우 적은 기억용량 만으로 Green 텐서 적분 행렬을 구성할 수 있어, 역산법에 효과적으로 적용할 수 있다. 이 수치 모델링 알고리듬을 기본으로 하여 평활화 제한을 가한 최소자승 역산 알고리듬을 개발하였다. 이 역산 알고리듬을 지표 전자탐사 및 시추공-지표 전자탐사 등에 적용하여 PC에서도 빠르게 3차원 전자탐사 역산이 수행됨을 보였다.

• Proceedings of the KSEEG Conference
• /
• 2001.04a
• /
• pp.123-125
• /
• 2001

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.121-128
• /
• 2017

Interpretation of time-domain electromagnetic (TEM) data has been made mostly based on one-dimensional (1-D) inversion scheme in Korea. A proper interpretation of TEM data should employ 3-D TEM forward and inverse modeling algorithms. This study developed a 3-D TEM modeling algorithm using a finite difference time-domain (FDTD) method with staggered grid. In numerically solving Maxwell equations, fictitious displacement current is included based on an explicit FDTD method using a central difference approximation scheme. The developed modeling algorithm simulated a small-coil source configuration to be verified against analytic solutions for homogeneous half-space models. Further, TEM responses for a 3-D anomaly are modeled and analyzed. We expect that it will contribute greatly to the precise interpretation of TEM data.

• Proceedings of the KSEEG Conference
• /
• 2000.04b
• /
• pp.197-198
• /
• 2000

• Geophysics and Geophysical Exploration
• /
• v.8 no.3
• /
• pp.207-217
• /
• 2005

We present a new approximate formulation of the integral equation (IE) method for models with variable background conductivity. This method overcomes the standard limitation of the conventional If method related to the use of a horizontally layered background only. The new approximate IE method still employs the Green's functions for a horizontally layered 1-D model. However, the new method allows us to use an inhomogeneous background with the IE method. The method was carefully tested for modeling the EM field for complex structures with a known variable background conductivity. It can find wide application in modeling EM data for multiple geological models with some common geoelectrical features, like a known inhomogeneous overburden, or salt dome structures.

• Geophysics and Geophysical Exploration
• /
• v.17 no.3
• /
• pp.163-170
• /
• 2014

Development of a modeling technique for accurately interpreting electromagnetic (EM) data is increasingly required. We introduce finite difference (FD) and finite-element (FE) methods for three-dimensional (3D) frequency-domain EM modeling. In the controlled-source EM methods, formulating the governing equations into a secondary electric field enables us to avoid a singularity problem at the source point. The secondary electric field is discretized using the FD or FE methods for the model region. We represent iterative and direct methods to solve the system of equations resulting from the FD or FE schemes. By applying the static divergence correction in the iterative method, the rate of convergence is dramatically improved, and it is particularly useful to compute a model including surface topography in the FD method. Finally, as an example of an airborne EM survey, we present 3D modeling using the FD method.

• 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.

• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.13-18
• /
• 2000

We describe the concept of the singularity in the integral equation of electromagnetic (EM) modeling in comparison with that in the integral representation of electric fields in EM theory, which would clarify the singular integral problems of the Green's tensor. We have also derived and classified the singular integrals of the Green's tensors in 3-D, 2.5-D and 2-D as well as in the thin sheet integral equations of the EM scattering problem, which have the most important effect on the accuracy of the numerical solution of the problems.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.70-79
• /
• 2001

The integral equation method is a powerful tool for numerical electromagnetic modeling. But the difficulty of this technique is the size of the linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. The conventional Born approximation, where the electric field in the anomalous body is approximated by the background field, is very rapid and easy to compute. However, the technique is inaccurate when the conductivity contrast between the body and the background medium is large. Quasi-linear, quasi-analytical and extended Born approximations are novel approaches to 3-D EM modeling based on the linearization of the integral equations for scattered EM field. These approximation methods are much less time consuming than full integral equation method and more accurate than conventional Born approximation. They we, however, still approximate methods for 3-D EM modeling. Iterative series methods such as modified Born, quasi-linear and quasi-analytical can be used to increase the accuracy of various approximation methods. Comparisons of numerical performance against a full integral equation and various approximation codes show that the iterative series methods are very accurate and almost always converge. Furthermore, they are very fast and easy to implement on a computer. In this study, extended Born series method is developed and it shows more accurate result than that of other series methods. Therefore, Iterative series methods, including extended Born series, open principally new possibilities for fast and accurate 3-D EM modeling and inversion.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.10a
• /
• pp.75-80
• /
• 2009

Despite its flexibility to complex geometry, three-dimensional (3D) electromagnetic(EM) modeling schemes using finite element method (FEM) have been faced to practical limitation due to the resulting large system of equations to be solved. An efficient 3D FEM modeling scheme has been developed, which can adopt either direct or iterative solver depending on the problems. The direct solver PARDISO can reduce the computing time remarkably by incorporating parallel computing on multi-core processor systems, which is appropriate for single frequency multi-source configurations. When limited memory, the iterative solver BiCGSTAB(1) can provide fast and stable convergence. Efficient 3D simulations can be performed by choosing an optimum solver depending on the computing environment and the problems to be solved. This modeling includes various types of controlled-sources and can be exploited as an efficient engine for 3D inversion.