Three-Dimensional Marine Magnetotelluric Parallel Forward Modeling in Conductive and Magnetic Anisotropic Medium Using Finite-Element Method Based on Secondary Field

The marine magnetotelluric (MMT) method is a significant tool extensively utilized in offshore studies, including the understanding of the Earth’s tectonics and hydrocarbon exploration. Conductive anisotropy and non-zero magnetic susceptibility are common phenomena observed in the Earth’s subsurface, and MMT forward modeling is the basis of practical inversion. However, numerical modeling that incorporates both conductive anisotropy and magnetic susceptibility has received limited attention. Moreover, both accuracy and efficiency are crucial in developing a 3D MMT modeling algorithm. Therefore, we developed a multi-level parallel MMT forward modeling algorithm that is capable of simultaneously modeling conductive and magnetic arbitrary anisotropic models using the vector finite element method based on the secondary field formula. The algorithm’s accuracy was validated through comparisons with previously published results for an arbitrary anisotropic model. The results show that the maximum relative error is below 2%, and the speedup reaches an impressive value of 552.41 when running with 2048 cores. Furthermore, the MMT responses of conductive anisotropy and magnetic susceptibility were comprehensively analyzed by several typical models. Our findings highlight the importance of considering magnetic susceptibility in magnetite-rich regions, particularly as the MMT responses may exhibit opposite responses for anomalies with lower resistivity and higher magnetic susceptibility compared with the surrounding rocks.