Efficient three-dimensional magnetotelluric forward modeling based on a geometric multigrid preconditioner
Wang YongFei,Liu JianXin,Guo RongWen,Liu Rong,Li Jian,Chen Hang,Yang GangQiang
DOI: https://doi.org/10.6038/cjg2022P0110
2022-01-01
Chinese Journal of Geophysics
Abstract:The geometric multigrid (GMG) method approximates the fine grid problems with the large sparse matrix on a coarser grid, which is easier to be solved. However, due to the double curl operator involved in magnetotelluric (MT) forward modeling, the high-frequency error cannot be smoothed effectively, resulting in an algorithm converging slowly or even diverging. For this reason, we introduce a four-color cell-block Gauss-Seidel (GS) method as smoother to improve the convergence efficiency. The algorithm locally satisfies the current divergence-free condition with high parallelism, and no additional divergence correction is needed which improves the efficiency of GMG significantly. However, with the anisotropy of the coefficient matrix (such as sharp change of conductivity, etc.), the efficiency of the GMG method can deteriorate. Krylov subspace methods such as BiCGstab can improve the problem. Therefore, in this paper, an efficient MT forward modeling method based on the staggered grid finite difference method (FDM) is proposed by combining GMG based on the four-color cell-block GS smoother and BiCGstab. In this method, GMG based on the four-color cell-block GS smoother is used as a preconditioner for the BiCGstab solver to significantly improve the forward modeling efficiency. We designed a layered resistivity model to verify the accuracy of the proposed algorithm by comparing with its analytical solution. Then, a double anomaly body resistivity model and a Dublin Model 1 (DTM1) were designed. Based on BiCGstab, the numerical performance of the GMG preconditioner was compared with other traditional preconditioners, such as the symmetric successive over relaxation (SSOR), block incomplete LU (ILU) factorization and the GS. The results show that the algorithm proposed in this paper is much better than the traditional preconditioners in terms of iteration number, compute time and stability. In all examples, the GMG preconditioner achieves convergence within 10 times, and the compute time is reduced by more than 70% compared with the traditional preconditioners, indicating the stability and high efficiency of the proposed method.