A shape reconstruction method for electromagnetic tomography using adjointfields and level sets

A two-step shape reconstruction method for electromagnetic (EM) tomography is presented which uses adjoint fields and level sets. The inhomogeneous ba ckground permittivity distribution and the values of the permittivities in some penetrable obstacles are assumed to be known, and the number sizes, sh apes, and locations of these obstacles have to be reconstructed given noisy limited-view EM data. The main application we address in the paper is the imaging and monitoring of pollutant plumes in environmental cleanup sites b ased on cross-borehole EM data. The first step of the reconstruction scheme makes use of an inverse scattering solver which recovers equivalent scatte ring sources for a number of experiments, and then calculates from these an approximation for the permittivity distribution in the medium. The second step uses this result as an initial guess for solving the shape reconstruct ion problem. A key point in this second step is the fusion of the 'level se t technique' for representing the shapes of the reconstructed obstacles, an d an 'adjoint field technique' for solving the nonlinear inverse problem. I n each step, a forward and an adjoint Helmholtz problem are solved based on the permittivity distribution which corresponds to the latest best guess f or the representing level set function. A correction for this level set fun ction is then calculated directly by combining the results of these two run s. Numerical experiments are presented which show that the derived method i s able to recover one or more objects with nontrivial shapes given noisy cr oss-borehole EM data.