scholarly journals Electromagnetic modeling of three-dimensional bodies in layered earths using integral equations

1982 ◽  
Author(s):  
P.E. Wannamaker ◽  
G.W. Hohmann
Keyword(s):  
Integral Equations ◽  
Three Dimensional ◽  
Electromagnetic Modeling

Electromagnetic modeling of 3-D structures by the method of system iteration using integral equations

Geophysics ◽  
1992 ◽  
Vol 57 (12) ◽  
pp. 1556-1561 ◽  
Author(s):  
Zonghou Xiong
Keyword(s):  
Integral Equations ◽  
Three Dimensional ◽  
Computation Time ◽  
Iterative Solution ◽  
Matrix Inversion ◽  
Electromagnetic Modeling ◽  
Direct Inversion ◽  
Conductivity Structure ◽  
The Matrix ◽  
Earth Conductivity

A new approach for electromagnetic modeling of three‐dimensional (3-D) earth conductivity structures using integral equations is introduced. A conductivity structure is divided into many substructures and the integral equation governing the scattering currents within a substructure is solved by a direct matrix inversion. The influence of all other substructures are treated as external excitations and the solution for the whole structure is then found iteratively. This is mathematically equivalent to partitioning the scattering matrix into many block submatrices and solving the whole system by a block iterative method. This method reduces computer memory requirements since only one submatrix at a time needs to be stored. The diagonal submatrices that require direct inversion are defined by local scatterers only and thus are generally better conditioned than the matrix for the whole structure. The block iterative solution requires much less computation time than direct matrix inversion or conventional point iterative methods as the convergence depends on the number of the submatrices, not on the total number of unknowns in the solution. As the submatrices are independent of each other, this method is suitable for parallel processing.

scholarly journals Electromagnetic modeling of three‐dimensional bodies in layered earths using integral equations

Geophysics ◽  
1984 ◽  
Vol 49 (1) ◽  
pp. 60-74 ◽  
Author(s):  
Philip E. Wannamaker ◽  
Gerald W. Hohmann ◽  
William A. SanFilipo
Keyword(s):  
Integral Equations ◽  
Green's Functions ◽  
Green’S Functions ◽  
Three Dimensional ◽  
Computation Time ◽  
Vector Current ◽  
Hankel Transform ◽  
Transverse Magnetic ◽  
Electromagnetic Modeling ◽  
Method Of Integral Equations

We have developed an algorithm based on the method of integral equations to simulate the electromagnetic responses of three‐dimensional bodies in layered earths. The inhomogeneities are replaced by an equivalent current distribution which is approximated by pulse basis functions. A matrix equation is constructed using the electric tensor Green’s function appropriate to a layered earth, and it is solved for the vector current in each cell. Subsequently, scattered fields are found by integrating electric and magnetic tensor Green’s functions over the scattering currents. Efficient evaluation of the tensor Green’s functions is a major consideration in reducing computation time. We find that tabulation and interpolation of the six electric and five magnetic Hankel transforms defining the secondary Green’s functions is preferable to any direct Hankel transform calculation using linear filters. A comparison of responses over elongate three‐dimensional (3-D) bodies with responses over two‐dimensional (2-D) bodies of identical cross‐section using plane wave incident fields is the only check available on our solution. Agreement is excellent; however, the length that a 3-D body must have before departures between 2-D transverse electric and corresponding 3-D signatures are insignificant depends strongly on the layering. The 2-D transverse magnetic and corresponding 3-D calculations agree closely regardless of the layered host.

Sign in / Sign up

Export Citation Format

Share Document