P/S separation of OBS data by inversion in a homogeneous medium

Next: Inversion results Up: Theory Previous: Forward and adjoint isotropic

## Inverting for the virtual source model

The inversion is done by least-square fitting of modeled geophone data to the recorded geophone data. The model is a virtual-source array, injected into some locations in a homogeneous medium. There should be some distance between these virtual-sources' locations and the receiver line where the displacement fields will be recorded. This distance is necessary for the recreated wavefield to form.

It is not possible to recreate the entire original wavefield that was recorded by the geophones without accurate knowledge of the acquisition geometry and the medium parameters. However, given a particular acquisition geometry, and some reasonable medium parameters, it is possible to recreate the original wavefield in the vicinity of the geophones. This wavefield will form as a result of the injection of the virtual sources, and will become more similar to the originally recorded wavefield at the receiver locations as it propagates toward them.

The inversion starts from a zero-value initial model of displacement sources . This initial model is injected as a displacement force function, at locations in the medium. The energy is propagated with the forward elastic operator, and recorded at locations in the medium (the receiver locations):

 (18)

The residual that must be minimized is the difference between the calculated data , and the observed data . The objective function is:

 (19)

 (20)

where is the adjoint elastic propagation operator. The data gradient is the forward operator applied to the model gradient:

 (21)

The model gradient and the data residual are updated using an itterative minimization scheme (Claerbout and Fomel, 2011).

The end result of the inversion is the model which, when injected at locations in a homogeneous medium with the parameters used in the inversion, will generate recorded data at locations that are equal to the observed data. Even if the medium parameters that were used to generate the wavefield were not a good approximation to the true parameters, this difference will be minor with regard to the wavefield at a close proximity to the receivers. Therefore, if the data residual is zero, we can assume that the wavefield near the receivers has likewise been reliably recreated.

The displacement field values are calculated for all time steps by running the forward propagation operator using the final virtual-source model as an input:

 (22)

The separation operators in equations 3 and 4 are applied to the displacements fields at the receiver locations, in order to separate P-wave amplitudes from S-wave amplitudes. These operators require calculating spatial derivatives of the displacement fields in the horizontal and vertical directions. We now have omnidirectional wavefield values near the receivers, therefore applying these operators to the recreated wavefields will result in an approximation to the actual P and S-waves that were recorded by the geophones. Furthermore, the displacement separation operators in equations 6-8 can be used to separate P displacements from S displacements.

 P/S separation of OBS data by inversion in a homogeneous medium

Next: Inversion results Up: Theory Previous: Forward and adjoint isotropic

2012-05-10