Multicomponent wave separation using HOSVD/unimodal-ICA subspace method

Multicomponent sensor arrays now are commonly used in seismic acquisition to record polarized waves. In this article, we use a three-mode model polarization mode, distance mode, and temporal mode to take into account the specific structure of signals that are recorded with these arrays, providing a data-structure-preserving processing. With the suggested model, we propose a multilinear decomposition named higher-order singular value decomposition and unimodal independent component analysis HOSVD/unimodal ICA to split the recorded threemode data into two orthogonal subspaces: the signal and noise subspaces.This decomposition allows the separation and identification of polarized waves with infinite apparent horizontal propagation velocity.The HOSVD leads to a definition of a subspace method that is the counterpart of the well-known subspace method for matrices that is driven by singular value decompositionSVD,aclassictoolinmonocomponentarrayprocessing. The proposed three-mode subspace decomposition provides a multicomponentwave-separationalgorithm.Toenhancetheseparation result, when the signal-to-noise ratio is low or when orthogonality constraints are not well adapted to the recorded waves, a unimodal-ICA step is included on the temporal mode. Doingthisreplacestheclassicorthogonalityconstraintsbetween estimated waves with independence constraints that might allow better recovery of recorded seismic waves.Asimulation on realistic two-component 2C geophysical data shows qualitative and quantitative improvements for the wavefield-separation results. The relative-mean-square errors between the original and estimated signal subspaces are, respectively, 52% for SVD applied on each component separately, 27.4% for HOSVD-based technique applied to the whole three-mode dataset, and 7.3% for HOSVD/unimodal-ICA technique. The efficiency of the threemodesubspacedecompositionsalsoisshownonrealthree-component 3C geophysical data. These results emphasize the potentialoftheHOSVD/unimodal-ICAsubspacemethodformulticomponentseismic-waveseparation.

[1]  Pieter M. Kroonenberg,et al.  Three-mode principal component analysis : theory and applications , 1983 .

[2]  Arye Nehorai,et al.  Analysis of a polarized seismic wave model , 1996, IEEE Trans. Signal Process..

[3]  Jean-Louis Lacoume,et al.  Modified singular value decomposition by means of independent component analysis , 2004, Signal Process..

[4]  Joos Vandewalle,et al.  A Multilinear Singular Value Decomposition , 2000, SIAM J. Matrix Anal. Appl..

[5]  Rasmus Bro,et al.  MULTI-WAY ANALYSIS IN THE FOOD INDUSTRY Models, Algorithms & Applications , 1998 .

[6]  A. Laub,et al.  The singular value decomposition: Its computation and some applications , 1980 .

[7]  B. Everitt,et al.  Three-Mode Principal Component Analysis. , 1986 .

[8]  Nikos D. Sidiropoulos,et al.  Parallel factor analysis in sensor array processing , 2000, IEEE Trans. Signal Process..

[9]  Gemma Musacchio,et al.  Polarization filter with singular value decomposition , 2001 .

[10]  Pierre Comon,et al.  Independent component analysis, A new concept? , 1994, Signal Process..

[11]  Arye Nehorai,et al.  Acoustic vector-sensor array processing , 1994, IEEE Trans. Signal Process..

[12]  Arye Nehorai,et al.  Acoustic vector-sensor correlations in ambient noise , 2001 .

[13]  L. Lathauwer,et al.  Signal Processing based on Multilinear Algebra , 1997 .

[14]  G. Giannakis,et al.  A FAST LEAST SQUARES ALGORITHM FOR SEPARATING TRILINEAR MIXTURES , 2004 .

[15]  Tad J. Ulrych,et al.  Application of singular value decomposition to vertical seismic profiling , 1988 .

[16]  L. Tucker,et al.  Some mathematical notes on three-mode factor analysis , 1966, Psychometrika.

[17]  Pierre Comon,et al.  Tensor Decompositions, State of the Art and Applications , 2002 .

[18]  Gene H. Golub,et al.  Matrix computations , 1983 .

[19]  John E. Vidale,et al.  Complex polarization analysis of particle motion , 1986 .

[20]  J. Cardoso,et al.  Blind beamforming for non-gaussian signals , 1993 .

[21]  Arye Nehorai,et al.  Vector-sensor array processing for electromagnetic source localization , 1994, IEEE Trans. Signal Process..

[22]  Nicolas Le Bihan,et al.  Three-mode data set analysis using higher order subspace method: application to sonar and seismo-acoustic signal processing , 2004, Signal Process..

[23]  J. Mari,et al.  4C-Sensor Orientation Using Multicomponent-SVD Method Applied on Shallow Water Data , 2001 .

[24]  G. M. Jackson,et al.  Principal component transforms of triaxial recordings by singular value decomposition , 1991 .

[25]  E. Deprettere SVD and signal processing: algorithms, applications and architectures , 1989 .

[26]  P. Comon Tensor Diagonalization, A useful Tool in Signal Processing , 1994 .

[27]  Jerome Mars,et al.  Advanced signal processing tools for dispersive waves , 2004 .

[28]  K. T. Wong,et al.  Self-initiating MUSIC-based direction finding and polarization estimation in spatio-polarizational beamspace , 2000 .