Although I have already written a Fortran77 program to solve a system of equations similar to (3) by the conjugate gradient method, Figure 2 shows that this program needs improvement in order to perform the inversion automatically without running into problems of edge effects.

The general procedure for estimating the velocity model will be to iterate until an acceptable velocity model is calculated, as follows:

- The mapping of the seismic velocity data onto the model space gives an initial guess of the model.
- The standard damped least-squares inversion is applied. At iteration
*k*, I- first invert for the model perturbation
(6) - then update the velocity model
(7)

- first invert for the model perturbation

The convergence is achieved when the mistie between the well velocity and the calculated velocities has been eliminated.

11/12/1997