摘要

In this paper, we advance a numerical simulation of electromagnetic responses of three-dimensional (3-D) electrical anomalous body in a layered anisotropic earth using integral equations. First, we decompose the dyadic Green's functions in a horizontal layered antsotropic earth into the singular primary term directly from the source and another non-singular waves reflected and transmitted from the interfaces. Then, iterative method is introduced to reduce computer memory and improve computational efficiency based on Krylov subspace. Finally, we validate the algorithm by comparison with results from former workers and investigate the effects of anisotropy on the electromagnetic responses of 3-D electrical anomalous body in a layered anisotropic earth through numerical models.