Seismic waves in stratified anisotropic media

Summary. The response of a structure composed of anisotropic strata can be built up from the reflection and transmission properties of individual interfaces using a slightly modified version of the recursion scheme of Kennett. This scheme is conveniently described in terms of scatterer operators and scatterer products. The effects of a free surface and the introduction of a simple point source at any depth can be accommodated in a manner directly analogous to the treatment for isotropic structures. As in the isotropic case the results so obtained are stable to arbitrary wavenumbers. For isotropic media, synthetic seismograms can be constructed by computing the structure response as a function of frequency and radial wavenumber, then performing the appropriate Fourier and Hankel transforms to obtain the wavefield in time-distance space. Such a scheme is convenient for any system with cylindrical symmétry (including transverse isotropy). Azimuthally anisotropic structures, however, do not display cylindrical symmétry; for these the transverse component of the wavenumber vector will, in general, be non-zero, with the result that phase, group, and energy velocities may all diverge. The problem is then much more conveniently addressed in Cartesian coordinates, with the frequency-wavenumber to time-distance transformation accomplished by 3-D Fourier transform.