Robust numerical algorithm to model induction tool responses in planar-layered media

We present a robust, rapid, and computationally efficient forward solver to compute electromagnetic (EM) fields produced by active, induction-regime remote sensing instruments embedded in subsurface, planar-stratified environments containing media of arbitrary anisotropy and loss. The algorithm is based on the robust evaluation of the spectral integral representation for the fields in such complex media, synthesizing the spacedomain EM field as a superposition of plane wave fields. To validate the algorithm´s versatility and efficacy, we present and discuss results involving electromagnetic induction tools used for geophysical prospection and exploration in planar-stratified geophysical formations containing media with isotropic and uniaxial-conductive properties.