Efficient magnetotelluric forward modeling in the spherical coordinate system based on a regularization constraint

Guo RongWen,Hou Sheng,Liu JianXin,Li Jian,Wang YongFei,Li RuiQin

DOI: https://doi.org/10.6038/cjg2022Q0651

2023-01-01

Abstract:With the availability of Magnetotelluric (MT) data on transcontinental scale, it is necessary to carry out 3D MT forward modeling and inversion in the spherical coordinate system to avoid the errors caused by the curvature of the Earth. The conductivity term in the curl-curl equation governing the electromagnetic (EM) fields can be neglected, due to the low frequencies used for this kind of forward modeling and inversion problems. And it is difficult to model the charge accumulation caused by the conductivity discontinuity when the numerical methods are used. The convergence of iterative solver would be slowly or even diverged for numerical solution, even preconditioned properly and with the application of divergence correction. To address this problem, the scaled correction term is explicitly added to the curl-curl equation (called the regularization constraint term for convenience) to ensure the divergence of current of each iteration is zero. The method avoids additional solving of the divergence equation and significantly improves the efficiency of 3D MT forward modeling in spherical coordinate. The regularized curl-curl equation is discretized with finite difference methods in the spherical coordinate and a high order spherical harmonic function P-1(0), is used to approximate the MT source for forward modeling. First, a layered model with the analytical solution is designed. Then a simple low resistivity model and a Cascadia model based on an inversion model from the measured data are designed to test the numerical performance of the algorithm. The comparison with analytical solution verifies the correctness of the algorithm. The numerical performance tests show that the algorithm significantly reduces the computational time and the number of iterations compared to the conventional algorithm based on iterative divergence correction. More importantly, the algorithm is stable over all the considered periods.