Electrical Impedance Tomography
SIAM Review
Iterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems
Inverse Problem Theory and Methods for Model Parameter Estimation
Inverse Problem Theory and Methods for Model Parameter Estimation
Magnetotelluric inversion via reverse time migration algorithm of seismic data
Journal of Computational Physics
Resistivity probability tomography imaging at the castle of Zena, Italy
Journal on Image and Video Processing - Special issue on image and video processing for cultural heritage
Hi-index | 31.45 |
We propose an electric resistivity inversion method that is similar to the reverse time migration technique applied to seismic data. For calculating model responses and inversion, we use the mixed finite-element method with the standard P"1-P"0 pair for triangular decompositions, which makes it possible to compute both the electric potential and the electric field vector economically. In order to apply the adjoint state of the Poisson equation in the resistivity inverse problem, we introduce an apparent electric field defined as the dot product between the computed electric field vector and a weighting factor and then defining a virtual source to compute the partial derivative of the electric field vector. We exploit the adjoint state (the symmetry of Green's function) of matrix equations derived from solving the Poisson equation by the mixed finite-element method, for the calculation of the steepest descent direction of our objective function. By computing the steepest descent direction by a dot product of backpropagated residual and virtual source, we can avoid the cumbersome and expensive process of computing the Jacobian matrix directly. We calibrate our algorithm on a synthetic of a buried conductive block and obtain an image that is compatible with the limits of the resistivity method.