摘要

To eliminate the coupling between air and the earth and to avoid calculation of primary field, we develop a universal algorithm on forward modelling of the frequency-domain electromagnetic (EM) method for different source types and applications. We first present a magnetic vector and a scalar potential with Lorenz gauge based on Maxwell equation, then divide different sources into a lot of short wires. Next, we use staggered grids and the finite volume method to discrete the potential equations and implement the quasi-minimum-residual (QMR) iteration with Jacob to solve the large sparse and symmetrical linear equations system. Finally, we accomplish frequency domain EM modelling for different source types and application areas. Through analyses and comparison of 1D and 3D MT, CSAMT and loop source responses, we demonstrate the efficiency and accuracy of the algorithm proposed in this study. Based on this, we analyze EM responses for different source types. The algorithm and numerical results presented in this paper build a framework for 3D frequency-domain EM forward modelling and inversion.