A Born-WKBJ inversion method for acoustic reflection data

Density and bulk modulus variations in an acoustic earth are separately recoverable from standard reflection surveys by utilizing the amplitude-versus-offset information present in the observed wave fields. Both earth structure and a variable background velocity can be accounted for by combining the Born and WKBJ approximations, in a "before stack" migration with two output sections, one for density variations and the other for bulk modulus variations. For the inversion, the medium is considered to be composed of a known low-spatial frequency variation (the background) plus an unknown high-spatial frequency variation in bulk modulus and density (the reflectivity). The division between the background and the reflectivity depends upon the frequency content of the source. For constant background parameters, computations are done in the Fourier domain, where the first part of the algorithm includes a frequency shift identical to that in an F-K migration. The modulus and density variations are then determined by observing in a least-squares sense amplitude versus offset wavenumber. For a spatially variable background, WKBJ Green's operators that model the direct wave in a medium with a smoothly varying background are used. A downward continuation with these operators removes the effects of variable velocity from the problem, and, consequently, the remainder of the inversion essentially proceeds as if the background were constant. If the background is strictly depth dependent, the inversion can be expressed in closed form. The method neglects multiples and surface waves and it is restricted to precritical reflections. Density is distinguishable from bulk modulus only if a sufficient range of precritical incident angles is present in the data