Seismogram synthesis for radially layered media using the generalized reflection/transmission coefficients method: Theory and applications to acoustic logging

A new method based on generalized reflection and transmission coefficients is proposed to calculate the synthetic seismograms in radially multilayered media. This method can be used to efficiently simulate full waveform acoustic logs and crosswell seismic profiles in situations where we need to consider borehole effects. The new formulation is tested by comparing our numerical results with previous available work and shows excellent agreement. Because of the use of the normalized Hankel functions and the normalization factors, this new algorithm for computing seismograms is stable numerically even for high-frequency problems. To show the applicability of this new approach to full waveform sonic logging, we apply it to investigate the effects of complex invaded zones on the geometrical spreading and attenuation estimation for P-waves. We find that a damaged zone (its velocity is slower than the unperturbed formation velocity) exhibits a convergence effect on the P-waves, and a flushed zone (velocity is faster ABSTRACT