Stable evaluation of Green's functions in cylindrically stratified regions with uniaxial anisotropic layers

Abstract

We present a robust algorithm for the computation of electromagnetic fields radiated by point sources (Hertzian dipoles) in cylindrically stratified media where each layer may exhibit material properties (permittivity, permeability, and conductivity) with uniaxial anisotropy. Analytical expressions are obtained based on the spectral representation of the tensor Green's function based on cylindrical Bessel and Hankel eigenfunctions, and extended for layered uniaxial media. Due to the poor scaling of these eigenfunctions for extreme arguments and/or orders, direct numerical evaluation of such expressions can produce numerical instability, i.e., underflow, overflow, and/or round-off errors under finite precision arithmetics. To circumvent these problems, we develop a numerically stable formulation through suitable rescaling of various expressions involved in the computational chain, to yield a robust algorithm for all parameter ranges. Numerical results are presented to illustrate the robustness of the formulation including cases of practical interest to geophysical exploration.

References in corpus (4)

• 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
• Computation of potentials from current electrodes in cylindrically stratified media: A stable, rescaled semi-analytical formulation