3D general-measure inversion of crosswell EM data using a direct solver

Abstract

Abstract We present a three-dimensional (3D) general-measure inversion scheme of crosswell electromagnetic (EM) data in the frequency domain with a direct forward solver. In the forward problem, we discretised the EM Helmholtz equation by the staggered-grid finite difference (SGFD) scheme and solved it using the Intel MKL PARDISO direct solver. By applying a direct solver, we simultaneously solved the multisource forward problems at a given frequency. In the inversion, we integrated a general measure of data misfit and model constraints with linearised least-squares inversion. We reconstructed a model with blocky features by selecting the appropriate measure parameters and model constraints. We used the adjoint equation method to explicitly calculate the Jacobian matrix, which facilitated the determination of an appropriate initial value for the regularisation coefficient in the objective function. We illustrated the inversion scheme using synthetic crosswell EM data with a general measure, the L2 norm, and, specifically, two mixed norms.