Electromagnetic divergence correction for 3D anisotropic EM modeling

Until now, the application of divergence correction has been focusing on the electromagnetic (EM) isotropic modeling. However, in regions where the earth demonstrates strong electrical anisotropy due to lamina or faults saturated with water, the traditional divergence corrections may not be working for EM modeling. In this paper, based on the divergence-free property of the current density, we propose a divergence correction technique to speed up the EM modeling process for a 3D arbitrarily anisotropic earth. The volume current density weighting method is adopted for the discretization of divergence correction, so it can be applied to non-uniform grid model. By enforcing the constraints associated with the divergence-free property, the quasi-minimal residual (QMR) iterations and the total computational time for the EM modeling are largely reduced. Furthermore, we discuss how to divide the QMR process into sections and the influence of the number of iterations and the normalized residual on the solution process. We take a helicopter-borne EM system over a 3D earth with arbitrarily electrical anisotropy as an example to examine the effectiveness of our anisotropic divergence correction technique on the EM modelings.