Next: [3] The meaning and Up: Aspects of Linearized migration/inversion Previous: [1] Numerical calculation of

[2] The matrix expression of linearized migration/inversion

The linearized migration/inversion can be formulated from the integral expressed in equation (15). It can be regarded as an inverse generalized Radon transform. Equation (15) can be expressed as a matrix equation. The process for solving the equation set is just the migration/inversion imaging. Following Berkhout (1997), we first give a matrix expression of wave propagation from a source to a scatterer and then to a receiver:
 (58)
where WU is a discretized Green's function for upward wave propagation, WD is a discretized Green's function for downward wave propagation, and is a reflectivity matrix, which is related to the incident angle. If the variation of reflectivity with angle is neglected, becomes a diagonal matrix. The reflectivity in this case is assumed to be the normal reflectivity. In practice, the reflectivity of a reflector varies with the incident angle. This is called an AVO/AVA phenomenon in seismology. The prestack migration/inversion aims at estimating the angle reflectivity to evaluate lithological variations. On the other hand, the residual moveout of the angle reflectivity indicates whether the macro migration/inversion velocity is reasonable or not. The synthetic wave field can be written as follows:
 (59)
In equation (59), each element of the matrix P is a recorded seismic trace in the time domain and a recorded amplitude value for a shot-receiver pair in the frequency domain. Each column is a shot gather, and each row is a common receiver gather. Therefore, equation (59) can be regarded as the matrix expression of equation (15). The classical prestack migration can be formulated as the following:
 (60)
The detailed matrix expression of equation (60) is
 (61)
where is the image of the first layer. In , each row is an angle gather at an imaging point, and each column is a common angle gather. The multiplication of the pth row in the matrix by any column in the matrix corresponds to a detection focusing of a shot gather; the multiplication of the qth column in the matrix by any row in the matrix corresponds to an emission focusing. Then, the image of the second layer can be obtained with
 (62)
Generally, the image of the zith layer is
 (63)
Here, the matrices and are non-recursive. Otherwise, equation (62) and (63) will be of the following forms:
 (64)
or
 (65)
and
 (66)
or
 (67)
Defining the cost function as
 (68)
yields the formula of the linearized migration/inversion:
 (69)
The matrix expression of the migration/inversion in equation (69) is
 (70)
where the denominator term is
 (71)
From equation (61) and (70), the migration/inversion can be locally implemented, because all elements in the matrices WU, , WD and relate only to a given layer. If the matrix is the inverse of the matrix WU, it can be expressed as
 (72)
where the matrix E is an identity matrix. Similarly, if the matrix is the inverse of the matrix WD, we have
 (73)
In practice, and are the conjugates of WUand WD respectively. Therefore, the matrix is a band-width-limited diagonal matrix. The velocity structure and the acquisition geometry affect the inner structure of the matrix. In fact, M is a Hessian which will be discussed in detail later.

Next: [3] The meaning and Up: Aspects of Linearized migration/inversion Previous: [1] Numerical calculation of
Stanford Exploration Project
11/1/2005