3-D linearized scattering of surface waves and a formalism for surface wave holography

Summary. Scattering of surface waves by lateral heterogeneities is analysed in the Born approximation. It is assumed that the background medium is either laterally homogeneous, or smoothly varying in the horizontal direction. A dyadic representation of the Green's function simplifies the theory tremendously. Several examples of the theory are presented. The scattering and mode conversion coefficients are shown for scattering of surface waves by the root of an Alpine-like crustal structure. Furthermore a 'great circle theorem'in a plane geometry is derived. A new proof of Snell's law is given for surface wave scattering by a quarter-space. It is shown how a stationary phase approximation can be used to simplify the Fourier synthesis of the scattered wave in the time domain. Finally a procedure is suggested to do 'surface wave holography'.     

Time-lapse traveltime change of singly scattered acoustic waves

We present a technique based on the single-scattering approximation that relates time-lapse localized changes in the propagation velocity to changes in the traveltime of singly scattered waves. We describe wave propagation in a random medium with homogeneous statistical properties as a single-scattering process where the fluctuations of the velocity with respect to the background velocity are assumed to be weak. This corresponds to one of two end-member regimes of wave propagation in a random medium, the first being single scattering, and the second multiple scattering. We present a formulation that relates the change in the traveltime of the scattered waves to a localized change in the propagation velocity by means of the Born approximation for the scattered wavefield. We validate the methodology with synthetic seismograms calculated with finite differences for 2-D acoustic waves. Potential applications of this technique include non-destructive evaluation of heterogeneous materials and time-lapse monitoring of heterogeneous reservoirs.     

Surface wave tomography: global membrane waves and adjoint methods

We implement the wave equation on a spherical membrane, with a finite-difference algorithm that accounts for finite-frequency effects in the smooth-Earth approximation, and use the resulting 'membrane waves' as an analogue for surface wave propagation in the Earth. In this formulation, we derive fully numerical 2-D sensitivity kernels for phase anomaly measurements, and employ them in a preliminary tomographic application. To speed up the computation of kernels, so that it is practical to formulate the inverse problem also with respect to a laterally heterogeneous starting model, we calculate them via the adjoint method, based on backpropagation, and parallelize our software on a Linux cluster. Our method is a step forward from ray theory, as it surpasses the inherent infinite-frequency approximation. It differs from analytical Born theory in that it does not involve a far-field approximation, and accounts, in principle, for non-linear effects like multiple scattering and wave front healing. It is much cheaper than the more accurate, fully 3-D numerical solution of the Earth's equations of motion, which has not yet been applied to large-scale tomography. Our tomographic results and trade-off analysis are compatible with those found in the ray- and analytical-Born-theory approaches.     

Generalized Born inversion of seismic reflection data

Summary . Born inverse methods give accurate and stable results when the source wavelet is impulsive. However, in many practical applications (reflection seismology) an impulsive source cannot be realized and the inversion needs to be generalized to include an arbitrary source function. In this paper, we present a Born solution to the seismic inverse problem which can accommodate an arbitrary source function and give accurate and stable results. It is shown that the form of the generalized inversion algorithm reduces to a Wiener shaping ***filter, which is solved efficiently using a Levinson recursion algorithm. Numerical examples of synthetic and real field data illustrate the validity of our method.     

Generalized Born scattering of elastic waves in 3-D media

It is well known that when a seismic wave propagates through an elastic medium with gradients in the parameters which describe it (e.g. slowness and density), energy is scattered from the incident wave generating low-frequency partial reflections. Many approximate solutions to the wave equation, e.g. geometrical ray theory (GRT), Maslov theory and Gaussian beams, do not model these signals. The problem of describing partial reflections in 1-D media has been extensively studied in the seismic literature and considerable progress has been made using iterative techniques based on WKBJ, Airy or Langer type ansätze. In this paper we derive a first-order scattering formalism to describe partial reflections in 3-D media. The correction term describing the scattered energy is developed as a volume integral over terms dependent upon the first spatial derivatives (gradients) of the parameters describing the medium and the solution. The relationship we derive could, in principle, be used as the basis for an iterative scheme but the computational expense, particularly for elastic media, will usually prohibit this approach. The result we obtain is closely related to the usual Born approximation, but differs in that the scattering term is not derived from a perturbation to a background model, but rather from the error in an approximate Green's function. We examine analytically the relationship between the results produced by the new formalism and the usual Born approximation for a medium which has no long-wavelength heterogeneities. We show that in such a case the two methods agree approximately as expected, but that in a media with heterogeneities of all wavelengths the new gradient scattering formalism is superior. We establish analytically the connection between the formalism developed here and the iterative approach based on the WKBJ solution which has been used previously in 1-D media. Numerical examples are shown to illustrate the examples discussed.     

Inverse scattering of surface waves: imaging of near-surface heterogeneities

An efficient inverse scattering method is developed for imaging near-surface heterogeneities using scattered surface waves. Three dimensional elastodynamic wave propagation and scattering in a laterally invariant embedding medium is considered. The Born Approximation is used and the scattered wavefield is expressed as a domain type integral representation. The computation time of Green's tensor elements is reduced by considering the radial symmetry of the medium. The method is validated by numerical tests. Ultrasonic laboratory data obtained from a scale model experiment are used for imaging the near-surface inhomogeneities caused by an epoxy-filled hole in the surface of an aluminum block. Both synthetic and the scale model tests show that the location, the actual density contrast and the depth of the inhomogeneities are reasonably well estimated.