Asymptotic pseudounitary stacking operators |

The theoretical least-squares inverse of operator (1) has the well-known form (Tarantola, 1987)

where denotes pseudo-inverse, and the adjoint operator is defined by the dot-product test:

With a specified definition of the dot-product, the generalized inverse minimizes the following quantity, which is the squared norm of the residual:

(15) |

In the case of integral operators, a natural definition of the dot-product is the double integral

(16) |

(17) |

The notion of the adjoint operator completely depends on the
arbitrarily chosen definition of the dot product and norm in the model
and data spaces. A simple way to change those definitions is to find
some positive weights
in the model space and
in
the data space that define the dot products as follows:

To formally define the adjoint of a stacking operator, let us substitute the definition of the stacking operator (1) into the dot product (14), as follows:

Assuming that the function is monotone in

where has the same meaning as in equation (8), and

Comparing equations (21) and (14), we conclude that the adjoint operator is defined by the equality

Thus we have proven that the continuous adjoint of a stacking operator is another stacking operator. The adjoint operator has the same summation path as the asymptotic inverse (8), which guarantees the correct reconstruction of the kinematics of the input wavefield. The amplitude (weighting function) of the adjoint operator is directly proportional to the forward weighting according to equation (22). The coefficient of proportionality is the Jacobian of the transformation of the variables and .

Similar results have been obtained for particular cases of stacking operators: velocity transform (Jedlicka, 1989; Thorson, 1984), Kirchhoff constant-velocity migration (Ji, 1994), and NMO (Crawley, 1995). In the appendix, I exemplify an application of least-squares inversion by reviewing inversion of the Radon operator and showing that it is precisely equivalent to the asymptotic result of the previous section.

Asymptotic pseudounitary stacking operators |

2013-03-03