On the imaging of reflectors in the Earth

Abstract In this paper, I present a modification of the Beylkin inversion operator. This modification accounts for the band-limited nature of the data and makes the role of discontinuities in the sound speed more precise. The inversion presented here partially dispenses with the small-parameter constraint of the Born approximation. This is shown by applying the proposed inversion operator to upward scattered data represented by the Kirchhoff approximation, using the angularly dependent geometrical-optics reflection coefficient. A fully nonlinear estimate of the jump in sound speed may be extracted from the output of this algorithm interpreted in the context of these Kirchhoff-approximate data for the forward problem. The inversion of these data involves integration over the source-receiver surface, the reflecting surface, and frequency. The spatial integrals are computed by the method of stationary phase. The output is asymptotically a scaled singular function of the reflecting surface. The singular function of a surface is a Dirac delta function whose support is on the surface. Thus, knowledge of the singular functions is equivalent to mathematical imaging of the reflector. The scale factor multiplying the singular function is proportional to the geometrical-optics reflection coefficient. In addition to its dependence on the variations in sound speed, this reflection coefficient depends on an opening angle between rays from a source and receiver pair to the reflector. I show how to determine this unknown angle. With the angle determined, the reflection coefficient contains only the sound speed below the reflector as an unknown, and it can be determined. A recursive application of the inversion formalism is possible. That is, starting from the upper surface, each time a major reflector is imaged, the background sound speed is updated to account for the new information and data are processed deeper into the section until a new major reflector is imaged. Hence, the present inversion formalism lends itself to this type of recursive implementation. The inversion proposed here takes the form of a Kirchhoff migration of filtered data traces, with the space-domain amplitude and frequency-domain filter deduced from the inversion theory. Thus, one could view this type of inversion and parameter estimation as a Kirchhoff migration with careful attention to amplitude.