Abstract:We propose a full-wave pseudo-analytical numerical electromagnetic (EM) algorithm to model subsurface induction sensors, traversing planar-layered geological formations of arbitrary EM material anisotropy and loss, which are used, for example, in the exploration of hydrocarbon reserves. Unlike past pseudo-analytical planar-layered modeling algorithms that impose parallelism between the formation's bed junctions, our method involves judicious employment of Transformation Optics techniques to address challenges related to modeling relative slope (i.e., tilting) between said junctions (including arbitrary azimuth orientation of each junction). The algorithm exhibits this flexibility, both with respect to loss and anisotropy in the formation layers as well as junction tilting, via employing special planar slabs that coat each “flattened” (i.e., originally tilted) planar interface, locally redirecting the incident wave within the coating slabs to cause wave fronts to interact with the flattened interfaces as if they were still tilted with a specific, user-defined orientation. Moreover, since the coating layers are homogeneous rather than exhibiting continuous material variation, a minimal number of these layers must be inserted and hence reduces added simulation time and computational expense. As said coating layers are not reflectionless however, they do induce artificial field scattering that corrupts legitimate field signatures due to the (effective) interface tilting. Numerical results, for two half-spaces separated by a tilted interface, quantify error trends versus effective interface tilting, material properties, transmitter/receiver spacing, sensor position, coating slab thickness, and transmitter and receiver orientation, helping understand the spurious scattering's effect on reliable (effective) tilting this algorithm can model. Under the effective tilting constraints suggested by the results of said error study, we finally exhibit responses of sensors traversing three-layered media, where we vary the anisotropy, loss, and relative tilting of the formations and explore the sensitivity of the sensor's complex-valued measurements to both the magnitude of effective relative interface tilting (polar rotation) as well as azimuthal orientation of the effectively tilted interfaces.

A damped forward EMI model for a horizontally stratified earth

Modeling the 3D Terrain Effect on MT by the Boundary Element Method

Magnetotelluric Responses of an Anisotropic 1-D Earth with a Layer of Exponentially Varying Conductivity

Strength of the Electric Dipole Source Field in Multilayer Spherical Media

3D Airborne EM Forward Modeling Based on Time-Domain Spectral Element Method

1D forward modeling and inversion algorithm for grounded galvanic source TEM sounding with an arbitrary horizontal wire

3D transient electromagnetic inversion based on explicit finite-difference forward modeling

Rapid surrogate modeling of magnetotelluric in the frequency domain using physics-driven deep neural networks

Magnetotelluric Responses of a Horizontally Nonhomogeneous Model of the Earth with Anisotropic Conductivity Structures

Bayesian Inversion for Layered Spherical Symmetric Earth Conductivity Model from Global Magnetic Data

Using finite dipole lengths in complete earth 3D MT modelling

Understanding Of Near-Fleld Effect Fair Csamt In Electric Azimuthally Anisotropic Half-Space

Analytical Model for Eddy Current Testing of Thickness of Multi-layered Structures and Experimental Verification

Electromagnetic fields excited by arbitrarily orientated dipoles in horizontally layered earth

Forward Calculation of 3D Controlled-Source Electromagnetic Responses Based on Joint Application of Secondary Field and Coupled Potential Formulations

Electromagnetic Field Of A Horizontal In-Finitely Long Magnetic Line Source Over The Earth Coated With A Dielectric Layer

A Generic 1D Forward Modeling and Inversion Algorithm for TEM Sounding with an Arbitrary Horizontal Loop

3D Finite-Element Modeling of Earth Induced Electromagnetic Field and Its Potential Applications for Geomagnetic Satellites

Full-wave algorithm to model effects of bedding slopes on the response of subsurface electromagnetic geophysical sensors near unconformities

Forward modeling of marine DC resistivity method for a layered anisotropic earth

3D finite-difference transient electromagnetic modeling with a whole-space initial field