• Journal of IKEEE
• /
• v.23 no.2
• /
• pp.614-621
• /
• 2019

In this paper, we presents the integral equation-fast Fourier transform(IE-FFT) and block matrix preconditioner (BMP) to solve electromagnetic scattering problems of penetrable structures composed of dielectric or magnetic materials. IE-FFT can significantly improve the amount of calculation to solve the matrix equation constructed from the moment method(MoM). Moreover, the iterative method in conjunction with BMP can be significantly reduce the number of iterations required to solve the matrix equations which are constructed from electrically large structures. Numerical results show that IE-FFT and block matrix preconditioner can solve electromagnetic scattering problems for penetrable objects quickly and accurately.

• Journal of electromagnetic engineering and science
• /
• v.19 no.1
• /
• pp.6-12
• /
• 2019

An IE-FFT algorithm is implemented and applied to the electromagnetic (EM) solution of perfect electric conducting (PEC) scattering problems. The solution of the method of moments (MoM), based on the magnetic field integral equation (MFIE), is obtained for PEC objects with closed surfaces. The IE-FFT algorithm uses a uniform Cartesian grid to apply a global fast Fourier transform (FFT), which leads to significantly reduce memory requirement and speed up CPU with an iterative solver. The IE-FFT algorithm utilizes two discretizations, one for the unknown induced surface current on the planar triangular patches of 3D arbitrary geometries and the other on a uniform Cartesian grid for interpolating the free-space Green's function. The uniform interpolation of the Green's functions allows for a global FFT for far-field interaction terms, and the near-field interaction terms should be adequately corrected. A 3D block-Toeplitz structure for the Lagrangian interpolation of the Green's function is proposed. The MFIE formulation with the IE-FFT algorithm, without the help of a preconditioner, is converged in certain iterations with a generalized minimal residual (GMRES) method. The complexity of the IE-FFT is found to be approximately $O(N^{1.5})$and $O(N^{1.5}logN)$ for memory requirements and CPU time, respectively.

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

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.12 no.6
• /
• pp.571-579
• /
• 1987

The integral equations used in electromagnetic fields theory can be used for scattering problems. We can obtain various characteristics of scatterer. Ie, power pattern, scattered field, by finding current distribution on the scatterer. In this paper, current distribution on polygonal cylinder is obtained using integral equations in 2 dimension. For numerical aualysis, the moment method is used with pulse function as a basis function and integral equation is used with extrapolation method, which saves cpu time.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.167-169
• /
• 2002

The extended Born, or localized nonlinear approximation of integral equation (IE) solution has been applied to inverting single-hole electromagnetic (EM) data using a cylindrically symmetric model. The extended Born approximation is less accurate than a full solution but much superior to the simple Born approximation. When applied to the cylindrically symmetric model with a vertical magnetic dipole source, however, the accuracy of the extended Born approximation is greatly improved because the electric field is scalar and continuous everywhere. One of the most important steps in the inversion is the selection of a proper regularization parameter for stability. Occam's inversion (Constable et al., 1987) is an excellent method for obtaining a stable inverse solution. It is extremely slow when combined with a differential equation method because many forward simulations are needed but suitable for the extended Born solution because the Green's functions, the most time consuming part in IE methods, are repeatedly re-usable throughout the inversion. In addition, the If formulation also readily contains a sensitivity matrix, which can be revised at each iteration at little expense. The inversion algorithm developed in this study is quite stable and fast even if the optimum regularization parameter Is sought at each iteration step. Tn this paper we show inversion results using synthetic data obtained from a finite-element method and field data as well.

• Geophysics and Geophysical Exploration
• /
• v.7 no.1
• /
• pp.51-55
• /
• 2004

This paper presents a fast and stable imaging scheme using the localized nonlinear (LN) approximation of integral equation (IE) solutions for inverting electromagnetic data obtained in a crosswell survey. The medium is assumed to be cylindrically symmetric about a source borehole, and to maintain the symmetry a vertical magnetic dipole is used as a source. To find an optimum balance between data fitting and smoothness constraint, we introduce an automatic selection scheme for a Lagrange multiplier, which is sought at each iteration with a least misfit criterion. In this selection scheme, the IE algorithm is quite attractive for saving computing time because Green's functions, whose calculation is a most time-consuming part in IE methods, are repeatedly re-usable throughout the inversion process. The inversion scheme using the LN approximation has been tested to show its stability and efficiency, using both synthetic and field data. The inverted image derived from the field data, collected in a pilot experiment of water-flood monitoring in an oil field, is successfully compared with that derived by a 2.5-dimensional inversion scheme.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.60-67
• /
• 2008

In this paper we consider an application of the method of electromagnetic (EM) migration to the interpretation of a typical marine controlled-source (MCSEM) survey consisting of a set of sea-bottom receivers and a moving electrical bipole transmitter. Three-dimensional interpretation of MCSEM data is a very challenging problem because of the enormous number of computations required in the case of the multi-transmitter and multi-receiver data acquisition systems used in these surveys. At the same time, we demonstrate that the MCSEM surveys with their dense system of transmitters and receivers are extremely well suited for application of the migration method. In order to speed up the computation of the migration field, we apply a fast form of integral equation (IE) solution based on the multigrid quasi-linear (MGQL) approximation which we have developed. The principles of migration imaging formulated in this paper are tested on a typical model of a sea-bottom petroleum reservoir.