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

• Proceedings of the KSEEG Conference
• /
• 1999.04a
• /
• pp.327-329
• /
• 1999

• Proceedings of the KSEEG Conference
• /
• 1999.04a
• /
• pp.100-102
• /
• 1999

• Geophysics and Geophysical Exploration
• /
• v.1 no.3
• /
• pp.155-160
• /
• 1998

EM tomography technique has been developed. The algorithm used the extended Born approximations for forward modeling and reconstructed a conductivity image by a smoothness constraint least squares inversion method. Observed data, the vertical components of secondary magnetic fields, were simulated with the 3-D integral equation code. The results showed that the location of anomalous body could be imaged very well, but conductivity of the body was lower than real one and the vertical resolution was much higher than the horizontal resolution.

• 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.10 no.4
• /
• pp.308-313
• /
• 2007

The small-loop electromagnetic technique has been used successfully for many geophysical qualitative investigations, particularly for shallow engineering and environmental surveys. Recently, various geophysical imaging methods based on numerical modeling and inversion have been tried in order to get more quantitative subsurface structure. However, conventional 2.5D small loop EM modeling takes a lot of time because responses should be calculated for several wave numbers and transformed into space domain. In this study, we developed a 2.5D HCP small loop EM modeling algorithm using extended Born approximation, which does not require transformation. Also, we checked its validity by comparison with other numerical results.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.316-322
• /
• 2002

A computationally efficient inversion scheme has been developed using the extended Born or localized nonlinear approximation to analyze electromagnetic fields obtained in a single-hole environment. The medium is assumed to be cylindrically symmetric about the borehole, and to maintain the symmetry vertical magnetic dipole source is used throughout. The efficiency and robustness of an inversion scheme is very much dependent on the proper use of Lagrange multiplier, which is often provided manually to achieve desired convergence. In this study, an automatic Lagrange multiplier selection scheme has been developed to enhance the utility of the inversion scheme in handling field data. The inversion scheme has been tested using synthetic data to show its stability and effectiveness.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.