A synthetic data example is created to test the algorithm to make sure that it works properly. A simple layer-cake earth model is used, shown in Figure . RMS velocities are then created from this model as input data to the algorithm. If the inversion is run on this simple model, the result is almost perfect, as shown in Figure . The inversion is off by a maximum of 3%, which occurs at the bottom-most interface. This error could most likely be reduced further if we decrease the stopping criterion.
Now 1 and 5 percent Gaussian noise is added to the RMS velocities to simulate real data. The inversion of this noisy data with very little smoothness applied is shown in Figures and . The noise introduced to the model shows up as block features. As more noise is added the layers become harder to distinguish from each other.
If we increase the smoothing parameters on the regularization, then much of the noise is smoothed out in the result (Figure and Figure ). If the regularization parameters, are increased further then the result will be even smoother (Figure and Figure ).
As seen in these examples, it is important to correctly choose the regularization parameter to get a good inversion result that is compromise between desired blockiness and introducing spurious elements into the model in the form high spatial frequency events.
It can be seen that not all the noise is smoothed out in either Figure or Figure . This is because if boundaries are sharp then the regularization preserves them. The sharper the boundary, the higher the needs to be smooth them out. Much of the sharp contrast, however, is also smoothed away. From this test it became clear the smoothing along the midpoint is not currently working properly.