Next: (3) Leastsquares migration/inversion
Up: (A) Noniterative linearized migration/inversion
Previous: (1) Wave equation prestack
(a) Fourier Diffraction Tomography for constant background
Wu and Toksoz (1987) gave the planewave response in the direction from an incident wave :
 
(39) 
where is the 3D Fourier transform of the object function . is some kind of projection.
Comparing this to linear Radon transform, we know that the object function can be accurately restored if the angles of the plane waves continuously change around the object. Fig.3 shows the projection from the real plane wave source and from the virtual plane wave source.
(b) Inverse Generalized Radon Transform for variable background
The scattered wavefield after Bron and WKBJ approximation is of the following form:
 

 (40) 
where and .It is known that the diffractiontime surface in the data space is a counterpart of the isochron surface in the model space. These dual geometric associations naturally give rise to a corresponding pair of projection operators.
Equation (40) can be written as
 
(41) 
The diffraction curve in the data space is a projection of an isochron in the model space. This is a kind of Radon transform Hubral et al. (1996); Miller et al. (1987).
The standard Radon transform and inverse Radon transform in three dimensions are given by
 
(42) 
and
 
(43) 
respectively, where p is the distance from the origin to a plane which cuts through the object body, is the unity direction vector which is normal to the plane, and is a point on the plane.
GRT_fig
Figure 4 The geometry between the incident and scattering rays near the scattering point, or imaging point .

 
Comparing this with the classical Radon transform and its inverse, the final 3D inversion formula can be given as
 
(44) 
In equation (44), the angle variable near the imaging point is used, rather than the measurement configuration at the surface. Fig.4 illustrates this. The angle variable is related to the measurement configuration and reflects the seismic wave illumination aperture. Only if the aperture is large can a high resolution image be obtained. The relative trueamplitude imaging is severely affected by the angle variable. Bleistein and Stockwell (2001); Zhang (2004) gave some similar trueamplitude migration/inversion formulas.
Next: (3) Leastsquares migration/inversion
Up: (A) Noniterative linearized migration/inversion
Previous: (1) Wave equation prestack
Stanford Exploration Project
11/1/2005