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.
Electromagnetic modeling of three‐dimensional bodies in layered earths using integral equations