Polarization analysis in the wavelet domain based on the adaptive covariance method

We present an improved method for computing polarization attributes of particle motion from multicomponent seismic recordings in the time-frequency domain by using the continuous wavelet transform. This method is based on the analysis of the covariance matrix. We use an approximate analytical formula to compute the elements of the covariance matrix for a time window which is derived from an averaged instantaneous frequency of the multicomponent record. The length of the time-window is automatically and adaptively set to match the dominant period of the analysing wavelet at each time-frequency point. Then the eigenparameters are estimated for each time-frequency point without interpolation. With these key features, our method provides a suitable approach for polarization analysis of dispersive signals or overlapping seismic arrivals in multicomponent seismic data. For polarization analysis in the time domain, we show that the proposed method is consistent with existing polarization analysis methods. We apply the method to real data sets from exploration and earthquake seismology to illustrate some filtering applications and wave type characterizations.