Numerical evaluation of Green’s functions for axisymmetric boreholes using leaking modes
By choosing appropriate paths of integration in both the complex frequency (ω) and complex wavenumber (k) planes, exact Green’s functions for elastic wave propagation in axisymmetric boreholes are expressed completely as sums of modes. The integrations with respect to k are performed exactly using Cauchy residue theory. The remaining integrations with respect to ω are then carried out using the fast Fourier transform (FFT). The complete solution, including all possible body waves, is expressed simply as a superposition of modes without any contributions from branch line integrals. There are no spurious arrivals and, provided that the number of points in the FFT can be taken sufficiently large, no restrictions on distance. The method is fast, accurate, and easy to apply.