摘要
At present, the main problem in three-dimensional (3D) forward modeling is lower computational efficiency. An accurate, stable and broadband 3D EM forward modeling approach was introduced by Pankratov et al, and has been applied in MT forward modelling successfully. This approach uses Integral Equation (IE) method to solve Maxwell’s equations by modified Neumann series (MNS). It avoids calculating large-scale linear equations. In this paper, for the main problem, GPBi-CG is used to get the solution in modified Neumann series, and the efficiency of iteration is increased. Meanwhile, Green function is divided into two parts and solved in wave number domain. This is better than Fast Hankel Transform method. Finally, two models are tested. By comparing with results of 3D Staggered grid Finite Difference method, this approach is proved correct and effective.