Iterative deblending of simultaneous-source seismic data using seislet-domain shaping regularization |

where and are forward and inverse Fourier transforms, respectively, and is a diagonal block phase-shift operator given by

where denotes the individual phase shift operator for th trace and can be expressed as a diagonal matrix:

Here, denotes the angular frequency, denotes the random dithering time of th trace, and and in equations 10 and 11 denote the number of temporal samples and number of traces, respectively. denotes a diagonal matrix.

Considering that the Fourier operator (with symmetric normalization) and the phase shift operator are both unitary, which means and , it is easy to see that

Thus, we conclude that the dithering operator is also a unitary operator.

Furthermore we notice that

and:

The least-squares solution of equation 3 is therefore:

According to equation 14, an appropriate choice for is simply . This form of the derived backward operator is also referred to as scaled pseudodeblending (Mahdad, 2012).

The dithering operator is unitary only if the time-shift range is small and so the constructive summation of the useful components (events) between different traces can be ignored. Even if this condition isn't fully met, we can use the concept of interference to generalize the meaning of the dithering operator . Although, in this case, the backward operator might not be most appropriately chosen as half of the identity operator, we can still use it as an approximation.

Iterative deblending of simultaneous-source seismic data using seislet-domain shaping regularization |

2014-08-20