An efficient algebraic multi-resolution sampling approach to 3D magnetotelluric modeling

Summary Since electromagnetic (EM) fields diffuse more smoothly to greater depth, it physically makes sense to apply fine discretization to model structure at near surface with an increasingly coarser grid both in horizontal and vertical directions as the depth increases for the numerical solution of EM fields. For finite difference magnetotelluric (MT) forward modeling on regular staggered grids, this can lead to difficulties in discretizing the curl-curl equation governing the EM diffusion in the earth at regions, where the grid resolution changes. In this paper, we propose an efficient algebraic multi-resolution sampling (MRS) method for MT forward modeling. In this method, we do not require the generation of physical subgrids and merely subsample the field components. The coefficient matrix for the subsampled components can be obtained by matrix multiplication based on the initial linear system of equations and field interpolation. To guarantee divergence free current during the iterative solution process at low frequencies, which is considered crucial for the development of efficient iterative solvers, our forward modeling is based a regularization equation obtained by augmenting the curl-curl equation with an equivalent scaled grad-div operator for electrical fields (explicitly enforcing the divergence free condition). The correctness of our algebraic MRS algorithm is examined based on a layered model. Its stability and efficiency is exploited by comparing our results with the forward modeling on unsampled staggered grids for the Dublin test model 1 (DTM1) and a complex model with realistic topography, indicating a time reduction of about 42 per cent-82 per cent.