3-D Inversion of Z-Axis Tipper Electromagnetic Data Using Finite-Element Method With Unstructured Tetrahedral Grids

Abstract

<inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$Z$ </tex-math></inline-formula> -axis tipper electromagnetic (ZTEM) technique is an airborne electromagnetic method that detects anomalies in the deep earth that are induced by natural sources. Conventional ZTEM forward modeling is generally conducted using structured grids that have limited accuracy and cannot be used to invert complex underground structures and topography. However, unstructured grids can accurately model complex underground structures with fewer cells. Therefore, to effectively model and recover the complex underground structures, we developed a 3-D framework with unstructured tetrahedral grids for ZTEM forward modeling and inversions. The forward problem was formulated using the finite-element method with unstructured tetrahedral grids. To solve the inverse problem, we used a limited-memory quasi-Newton algorithm (L-BFGS) with a parallel direct solver for optimization in order to avoid the explicit calculation of the Hessian matrix and save the memory and computational time. To validate our forward and inversion algorithms, we conducted numerical experiments on three synthetic models and inverted a survey dataset acquired for a tunnel project in Tibet, China. The experimental results demonstrate the effectiveness of our unstructured finite-element and L-BFGS method for modeling and inverting ZTEM data.