Fast three-dimensional inversion of magnetotelluric data based on L-BFGS optimization

• Published in: Journal of applied geophysics Vol. 242; p. 105894
• Main Authors: Bai, Ningbo, Zhou, Jiancheng, Hu, Xiangyun, Han, Bo, Zhou, Junjun, Yang, Jian, Zhao, He
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01.11.2025
• Subjects: Hybrid solver; L-BFGS inversion; Magnetotelluric; Solution space dimensionality reduction; Magnetotelluric; Hybrid solver; Solution space dimensionality reduction; L-BFGS inversion
• ISSN: 0926-9851
• DOI: 10.1016/j.jappgeo.2025.105894

Three-dimensional (3D) inversion of magnetotelluric (MT) data is crucial for accurately resolving subsurface conductivity structures and requires robust, computationally efficient inversion techniques. In this study, we propose a novel 3D MT inversion method based on the limited-memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS) optimization algorithm, implemented within an edge-based finite element framework. The L-BFGS algorithm approximates the inverse Hessian using limited memory, thereby avoiding the storage and computation of large, dense matrices. To enhance the computational efficiency of the L-BFGS inversion, we introduce two strategies targeting the most time-consuming components: forward modeling and sensitivity calculations. The first strategy employs solution space dimensionality reduction by mapping the edges of the original forward-modeling grid to a coarser grid with fewer edges but without cells merging. This approach significantly reduces the degrees of freedom needed to solve the forward problem. The second strategy transforms the reduced-dimensional linear system into an equivalent real-valued linear system, which is efficiently solved using a direct–iterative hybrid solver. Numerical experiments demonstrate that solving the reduced-dimensional real linear system with the hybrid solver substantially decreases computational time compared to solving the original complex linear system with the direct solver PARDISO. The validity and robustness of the proposed inversion algorithm were confirmed through applications to both synthetic datasets and field MT data from the Longgang Volcanic Field. •This study realized a fast 3D inversion of MT data based the L-BFGS algorithm.•The SSDR scheme combined with a hybrid solver is used to accelerate the inversion process.•The algorithm is verified by synthetic model and field data.