Polarization filters are widely used for denoising seismic data. These filters are applied in the field of seismology, microseismic monitoring, vertical seismic profiling and subsurface imaging. They are primarily useful to suppress ground-roll in seismic reflection data and enhance P and S wave arrivals. Traditional implementations of the polarization filters involved the analysis of the covariance matrix or the SVD decomposition of a three-component seismogram matrix. The linear polarization level, known as rectilinearity, is expressed as a function of the covariance matrix eigenvalues or by the data matrix singular values. Wavefield records that are linearly polarized are amplified while others are attenuated. Besides the described implementations, we introduced a new time domain polarization filter based on the analysis of the seismic data correlation matrix. The principal idea is to extend the notion of the correlation coefficient in 3D space. This new filter avoids the need for diagonalization of the covariance matrix or SVD decomposition of data matrix, which are often time consuming. The new implementation facilitates the choice of the rectilinearity threshold: we demonstrate that linear polarization in 3D can be represented as three classic 2D correlations. A good linear polarization is detected when a high linear correlation between the three seismogram components is mutually observed. The tuning parameters of the new filter are the length of the time window, the filter order, and the rectilinearity threshold. Tests using synthetic seismograms show that optimal results are reached with a filter order that spans between 2 and 4, a rectilinearity threshold between 0.3 and 0.7, and a window length that is equivalent to one to three times the period of the signal wavelet. Compared to covariance-based filters, the new filter can enhance the signal-to-noise ratio by 6 to 20 dB and reduces computational costs by 25%.
