The traditional three-dimensional (3D) magnetotelluric (MT) forward modeling using Krylov subspace algorithms has the problem of low modeling efficiency. To improve the computational efficiency of 3D MT forward modeling, we present a novel geometric multigrid algorithm for the finite element method. We use the vector finite element to discretize Maxwell’s equations in the frequency domain and apply the Dirichlet boundary conditions to obtain large sparse complex linear equations for the solution of EM responses. To improve the convergence of the solution at low frequencies we use the divergence correction to correct the electric field. Then, we develop a V-cycle geometric multigrid algorithm to solve the linear equations system. To demonstrate the efficiency and effectiveness of our geometric multigrid method, we take three synthetic models (COMMEMI 3D-2 model, Dublin test model 1, modified SEG/EAEG salt dome model) and compare our results with the published ones. Numerical results show that the geometric multigrid algorithm proposed in this paper is much better than the commonly used Krylov subspace algorithms (such as SOR-GMRES, ILU-BICGSTAB, SOR-BICGSTAB) in terms of the iteration number, the solution time, and the stability, and thus is more suitable for large-scale 3D MT forward modeling.

Airborne electromagnetic (AEM) technology is an efficient geophysical exploration tool for investigating subsurface electrical structures. In recent years, 3D inversion of AEM data is being developed rapidly, but it still faces challenges such as low resolution and computational efficiency. To solve these problems, we propose a multiscale shearlet-based regularization inversion algorithm by establishing relationship between spatial resolution and shearlet coefficients in the inversion process. In the initial stage of inversion, the coarse grids and sparse measurement points data are used to recover the main subsurface structure. When the data misfit reaches a certain level, the previous results is used as the coarse scale model in shearlet-domain to recover the model with fine grids and dense measurements. By building this coarse-to-fine inversion scheme, we can well utilize the multiscale information in AEM data and effectively achieve high-resolution inversions. We demonstrate the effectiveness and practicality of our 3D MS inversion algorithm using two synthetic examples and a field dataset from Norway. The numerical experiments show that our inversion method can effectively reduce the computational time and improve inversion resolution.