Elastic wavefield directionality vectors

Next: P-wave and S-wave decomposition Up: Theoretical background Previous: Elastic wavefield modeling

## P-wave and S-wave amplitude separation and displacement decomposition

The Helmholtz amplitude separation is based on the assumption that any vector field can be described as resulting from a combination of a scalar and vector potential fields:

 (4)

Where is the scalar potential field and is the vector potential. is the elastic displacement vector wavefield. The scalar potential generates pressure waves, and the vector potential generates shear waves. Therefore, the Helmholtz method of separating the P-wave amplitude from the S-wave amplitude is to apply a divergence operator and a curl operator to the displacement field:

 (5)

 (6)

The derivation of P-wave and S-wave displacement decomposition by Zhang and McMechan (2010) is done in the wavenumber domain. For an isotropic medium, the linear equation system they arrive at is:

 (7)

 (8)

 (9)

and

 (10)

where is the 3D spatial fourier transform of the displacement field . and are the unknown P and S displacements. is the wavenumber vector that describes the particle displacement direction. , where is the phase velocity in the direction, and is angular frequency.

The solutions for P-wave displacements of these systems in 2D are:

 (11) (12)

The solutions for S-wave displacements are:

 (13) (14)

It is important to note that the in these equations is normalized by the absolute value of the wavenumber . Therefore, if we use to designate the non-normalized wavenumbers, each wavenumber must be divided by . Equations 12 - 14 then take the form:

 (15) (16)

 (17) (18)

As mentioned previously, these are decomposition operators for isotropic media only. I decided to use them initially, in order to evaluate the possibility of using elastic wavefield directionality for the purposes of angle gather creation. The main difference between displacement decomposition and the Helmholtz amplitude separation is that the decomposition is reversible. The sum of the decomposed P and S displacement vectors will result in the original pre-decomposed displacements. However, there is no theoretical way to go from the separated P and S amplitude fields back to the original displacement fields from which they were calculated.

 Elastic wavefield directionality vectors

Next: P-wave and S-wave decomposition Up: Theoretical background Previous: Elastic wavefield modeling

2011-05-24