Least-squares inversion of a velocity stack

To make the velocity stack usable for other processing steps such as multiple removal, I should try to incorporate the PEF described above into an inversion scheme that compensates for energy removed during the filtering.

Following the ``trial solution'' approach described in (), I change the variables 397#397 and solve equation () for the preconditioned variable 398#398:

   399#399 (158)

Operator 400#400 in equation () is a cascade of two filters described in the previous section. After finding a solution for 398#398,I evaluate 397#397 to obtain the solution of the original problem. To avoid high frequency noise in the model, I introduce a regularization term into the problem and solve the system of equations:

In this paper I used the laplacian as the regularization operator 8#8 in equation ().

Figure  shows the solution to the least-squares problem [equation ()] after 25 iterations of the conjugate-gradient method, the velocity stack after the filtering, and the solution to the preconditioned least-squares problem [equation ()] after 25 iterations of the conjugate-gradient.

 

prec
Figure 5 Left: Velocity stack after 25 iterations without preconditioning. Center: Velocity stack after filtering. Right: Velocity stack after 25 iterations with preconditioning.

Figure  shows the residual for the preconditioned problem [equation ()] and the problem without preconditioning [equation ()].

 

res_ann Figure 6 CG convergence with and without preconditioning

As Figures  and show, although the solution for the preconditioned problem does not have artifacts, it converges much slower than the solution to the problem without preconditioning, which probably makes this method of preconditioning impractical.