Full-Wave Algorithm to Model Effects of Bedding Slopes on the Response of Subsurface Electromagnetic Geophysical Sensors near Unconformities

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 however, 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 material and sensor characteristics. 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.

References in corpus (10)

• General Relativity in Electrical Engineering
• Optical design of reflectionless complex media by finite embedded coordinate transformations
• Robust Computation of Dipole Electromagnetic Fields in Arbitrarily-Anisotropic, Planar-Stratified Environments
• Stable pseudoanalytical computation of electromagnetic fields from arbitrarily-oriented dipoles in cylindrically stratified media
• Tensor Green's Function Evaluation in Arbitrarily Anisotropic, Layered Media using Complex-Plane Gauss-Laguerre Quadrature
• Complex-Plane Generalization of Scalar Levin Transforms: A Robust, Rapidly Convergent Method to Compute Potentials and Fields in Multi-Layered Media
• Spectral-Domain-Based Scattering Analysis of Fields Radiated by Distributed Sources in Planar-Stratified Environments with Arbitrarily Anisotropic Layers
• Computation of potentials from current electrodes in cylindrically stratified media: A stable, rescaled semi-analytical formulation
• Interface-Flattening Transform for EM Field Modeling in Tilted, Cylindrically-Stratified Geophysical Media
• Spectral-Domain Computation of Fields Radiated by Sources in Non-Birefringent Anisotropic Media