• Proceedings of the KSEEG Conference
• /
• 2001.04a
• /
• pp.53-55
• /
• 2001

• Geophysics and Geophysical Exploration
• /
• v.2 no.2
• /
• pp.86-95
• /
• 1999

This paper presents an efficient three-dimensional (3-D) modeling algorithm using the extended approximation to an electric field integral equation. Numerical evaluations of Green's tensor integral are performed in the spatial wavenumber domain. This approach makes it possible to reduce computing time, to handle smoothly varying conductivity model and to remove singularity problems encountered in the integration of Green's tensor at a source point. The responses obtained by 3-D modeling algorithm developed in this study are compared with those by the full integral equation for a thin-sheet EM scattering. The extensive analyses on the performance of modeling algorithm are made with the conductivity contrasts and source frequencies. These results show that the modeling algorithm are accurate up to the conductivity contrast of 1:16 and the frequency range of 100 Hz-100 kHz. The extended Born approximation, however, may produce inaccurate results for some source and model configurations in which the electric field is discontinuous across the conductivity boundary. We performed the modeling of a composite model of which conductivity varies continuously and this shows the modeling algorithm developed in this study is efficient for 3-D EM modeling. For a cross-hole source-receiver configuration a composite model of which conductivity varies continuously can be successfully simulated using this algorithm.

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

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