Printer Friendly

Using the Karhunen-Loeve transform to suppress ground roll in seismic data.


The Sacchi's algorithm (2002) based on the Karhunen-Loeve (K-L) Transform was modified and implemented to suppress Ground Roll without distortion of the reflection signals, it provided better results than conventional techniques for noise removal like f-k, High-Pass and Band Pass Filters. The K-L Transform is well known in other fields as image processing (Levy and Linderbaurn, 2000), face, iris and fingerprint identification. A seismic section is an image of subsurface where the K-L can be useful in seismic processing because spatially uncorrelated signals can be removed providing a clear and coherent image. The algorithm was applied to seismic data generated with hammer, thumper and explosive sources. Conventional processing flows were used, but one replaced filters with K-L Transform, providing stacked sections. The K-L Transform recovers better the reflector amplitudes when compared with others filters, also it removes refractions that cause unreal shallow events and increases the lateral coherence of seismic events showing a more interpretable geology.

Key words: Amplitude, coherent, filter, Ground Roll, noise.


En este estudio se implemento y modifico el Algoritmo de Sacchi (2002), basado en la Transformada Karhunen-Loeve (K-L), como una tecnica para la supresion de Ground Roll, mostrando que aunque la amplitud del Ground Roll es mayor que la amplitud de la senal al usar la Transformada K-L se tiene exito en suprimirlo sin causar distorsion en las senales de reflexion, obteniendo mejores resultados frente a tecnicas convencionales de remocion: f-k, Paso Alto y Pasa Banda. La Transformada K-L ha sido utilizada en el campo de la deteccion y procesamiento de imagenes (Levy and Linderbaurn, 2000), en el reconocimiento de caras, iris y huellas dactilares como factores biometricos de identidad. Una seccion sismica es una imagen del subsuelo por tanto esta transformada puede ser usada en el procesamiento sismico, porque remueve partes correlacionadas espacialmente y proveyendo una imagen clara y coherente.

El algoritmo se aplico a registros generados con Martillo, Thumper y fuentes explosivas. Se hizo un procesamiento convencional hasta el apilado, remplazando los filtros por le uso de la transformada K-L. Los resultados muestran que la Transformada K-L hace una mejor recuperacion de la amplitud de los reflectores, eliminando refracciones que causan eventos someros irreales e incrementando la coherencia lateral de los eventos sismicos, facilitando la interpretacion geologica.

Palabras clave: Amplitud, coherencia, filtro, Onda de Tierra, ruido.


The Ground Roll is (Scales and Snider, 1998) generated in shallow layers and is a typical Rayleigh wave. It is characterized by low velocity and frequency and high amplitude, being considered a coherent noise in seismic exploration (D'Agosto et al., 2002). The KL algorithm transforms the seismic data in a reasonable number of independent eigen functions that represents the most important signal characteristics. It has been considered a very advisable tool in the processing and analysis of reflection data (Al-Yahya, 1991).

The conventional filtering techniques: f-k, High-Pass and BandPass are applied in the frequency domain but in the case of Ground Roll amplitude stronger than reflection signals, f-k filters cause serious distortion of the signal (Karsh and Bayrak, 2004). High Pass filter is based on the low frequency of the Ground Roll, with the disadvantage that also eliminates the low-frequency content of reflection signals. Band pass filter depends on the amount of overlaps between the noise and the signal, therefore when the frequency bands concentrated the energy in each band without overlap, the signal/noise ratio will not change (Yilmaz, 1987).

The K-L Transform separates the Ground roll from the reflector signal, permitting to subtract it from seismic data. Its application to real data showed that the Ground Roll suppression increased lateral coherence of seismic events in a stacked section.


The surface waves in the raw register in figure 1 travel with velocities ranging from 100 to 1000 m/s and with frequencies around 10 Hertz and lower than those of reflections and refractions. Note how the near surface effects distort travel time curve at the right flank of reflection "A".


The Karhunen-Loeve Transform

It makes a principal component analysis that is a mathematical way of determining that linear transformation of a sample of points in L-dimensional space which exhibits the properties of the sample most clearly along the coordinate axes. Along the new axes, the sample variance are extremes and uncorrelated, see figure 2. Using a cutoff on the spread along each axis, a sample may be reduced in its dimensionality. This way it can be used to transform independent coordinates into significant and independent ones.


K-L Transform Implementation

The seismic traces [x.sub.i](t) correspond to the rows of the named data matrix [x.sub.n * m], n is the number of traces in the gather and m the number of gathers. The zero-lag covariance matrix [[GAMMA].sub.n * n] is (Jones and Levy, 1987):

[[GAMMA].sub.n * n] [X.sub.n * m] [X.sup.t.sub.m * n] (1)

This expression can be decomposed as

[[GAMMA].sub.n * n] [V.sub.n * p] [A.sub.n * n] [V.sup.t.sub.p * n]. (2)

In equation 2 the columns of the matrix [V.sub.n * p] are the eigenvectors [[GAMMA].sub.n * n] and [A.sub.n * n] is a diagonal matrix with the eigen values on the diagonal. These eigen values are ordered in descending sequence along the principal diagonal of the matrix [A.sub.n * n], [V.sup.t.sub.n * p] is the corresponding transpose matrix. The principal component of the data can be written by:

[M.sub.p * m] [V.sup.t.sub.p * n] [X.sub.n * m] (3)

If we form the matrix [M'.sub.p * m] by selecting the upper m rows of the matrix [M.sub.p * m], m is related with the order of the filter or principal row number, placing zeros in the remaining n-m rows, n is related with the Cut of the Filter, the Ground roll is assumed to be represented by the matrix product between the eigenvectors matrix [V.sub.n * p] and the matrix [M].sub.p * m], that include the principal components in agreement with the following expression:

[X'.sub.n * m] [V.sub.n * p] [M'.sub.p * m] (4)

Finally, the Ground roll [X'.sub.n * m] is subtracted from the data matrix [X.sub.n * m], obtaining the gather without Ground roll

[X.sub.On * m] [X.sub.n * m] [X'.sub.n * m] (5)

Xo is a n * m matrix that represents the filtered seismic gather.

As mentioned before [x.sub.n * m] is the input data matrix and [X'.sub.n * m] is a matrix that represents the Ground roll extracted from the input data.


The capacity of KL Algorithm to suppress Ground Roll at gather level is evaluated using shot gathers with dispersive waves of high amplitude, obtained from shallow and dip seismic surveys using 3 types of source were used: hammer, dynamite and thumper.


The seismic gather (figure 3) is a record corresponding to the UN-01 Seismic Line. Note in the red square like the Ground Roll diminishes before (left) and after (right) KL Algorithm application; the gather is much less noisy in the first 100 ms, with Order of the Filter: 2 and Cut of the Filter: 1. Figure 4 shows the interactive spectral analysis for Seismic Record No 12 of UN-01, before on top and after on bottom the application of the K-L Transform. The K-L Filtering gives good results because it compresses more the amplitude spectrum. Additionally, the maximum amplitude of the raw power spectrum with the K-L is lower into the red circle (10000) compared with one of the same gather without this filter in the black circle (4,5 e + 07), indicating high amplitude concentrations in the spectrum; corresponding to the Ground Roll suppressed with K-L. The record with K-L displays seismic events with coherence and continuity.


Explosive Source

The shot gather in Figure 5, was generated using dynamite, it displays hyperbolic moveout, coherent noise masking the reflection signals. The analyses with f-k, High-Pass and the KL Filters gave the better results for the last one.

The Band-Pass filter is used where the seismic trace contains low-frequency noise, such as Ground Roll and high-frequency ambient noise. The seismic reflection energy usually is confined to a bandwidth between 10 to 70 Hertz, with dominant frequency of 15 Hz.


Comparison between Band-Pass [15-30-55-70] (Figure 6), and KL Filter gave a better result for the KL, the Ground Roll cone has been suppressed and the lateral coherence of deeper seismic events increased.



Figure 7 shows an example, of the shot gather 481 with high content of Ground Roll. The Seismic Line RM-01, was acquired using a thumper source. A substantial increasing in the lateral coherence of the seismic events below the 200 ms is remarked.


Figure 8 shows the amplitude spectrum of the shot gather 481 before and after the KL filter application, a significant reduction of the high amplitude corresponding to Ground Roll is seen.


Shallow Survey

The 2D seismic data was acquired in October-2004 with a split spread configuration in the campus of the National University of Colombia, with a Geode seismic recorder with 24 channels for high resolution acquisition, the source consist in a hammer (8 kg) and geophones with vertical component.

The processing of the UN-01 Seismic Line provided stacked sections shown in Figure 9, on top conventional and on bottom with K-L. The section on bottom presents better results, increases lateral coherence and gives a clearer image. On the other hand, the KL suppresses unreal shallow events (indicated by the arrow), due to refractions that are observed with the conventional processing.


Deep Survey

This 2D data was acquired using a thumper in split spread configuration at Southern of the Middle Magdalena Valley, involving tertiary sediments from the Gualanday Group.

The stacked section in RM-01 in Figure 10 was obtained with conventional processing on bottom and with application of the K-L Transform on top. The top image presents a great rank of amplitude (-3 to 3) compared with the conventional stacked on bottom (-2 to 2) as observed into the red circles; indicating a better amplitude recovering in reflectors with K-L use, although the processing sequence did not include fx Decon. Besides the KL shows coherence and it does not present the whitening effect due to the conventional processing flow (Figure 10. Red square, bottom).


By the other hand, high dips (black arrow), are observed in conventional process image which are not true because the direction of seismic line was N-S, parallels to the beds strikes of the tertiary sediments with very smooth dips; in consequence they must be horizontals in the seismic section, as it is observed in the K-L image in agreement with the subsurface geology.


The Algorithm Karhunen-Loeve is computationally economic and efficient to filter seismic data, providing coherent information of subsurface image. The interactive spectral analysis showed more concentrations of high amplitude due to the Ground Roll, and it verifies that K-L Algorithm suppresses the Ground Roll waves. Also it is verified than the amplitude of the Ground Roll is stronger than the reflection signals. By using K-L Algorithm to extract the Ground Roll it suggests a form to suppress Ground Roll, without distorting the reflection signals.

Additionally, using K-L filter increases the lateral coherence of seismic events, in comparison with the three conventional methods of suppression: Filters f - k, High-Pass and Band-Pass. The algorithm was tested with raw gather containing Ground Roll (high amplitude) obtained from shallow and depth seismic land surveys with different source (hammer, dynamite and thumper. A lateral coherence increase is observed in stacked sections. The K-L use recovers the amplitude of the deep reflectors furnishing a resultant section more approximated to real geology.

(Manuscript received March 2005, Paper accepted October 2005)


* Al-Husseini, M., Glover J. and Barley B., 1981. Dispersion Patterns of the Ground Roll in Eastern Saudi Arabia, Geophysics, Vol 46, p. 121-137.

* Al-Yahya, K.M., 1991. Application of the Partial Karhunen-Loeve Transform to Suppress Random Noise in Seismic Sections: Geophys. Prosp., 39, p. 77-93.

* D'Agosto, C., Marfurt, K. and Steven, J., 2002. Modeling and Removal of Ground Roll from Horizontal Component of C- Waves. SEG.

* Jones, 1. F., and Levy, S., 1987. Signal-to-Noise Enhancement in Multichannel Seismic Data Processing: Geophys. Prosp., 35, p. 12-32.

* Karsh, H. and Bayrak, Y., 2004. Using the Wiener-Levinson Algorithm to Suppress Ground Roll. Journal of Applied Geophysics 55, p 187-197.

* Levy, A. and Linderbaurn, M., 2000. Sequential Karhunen-Loeve Basis Extraction and its Application to Images, IEEE Transactions on Image Processing, Vol 9, No 8, p. 1371-1374.

* Nesvijski, E. G., 2000. On a Possibility of Rayleigh Transformed Sub-Surface Waves Propagation. Center of Technology, Federal University of Santa Maria, Brasil.

* Sacchi, M.D., 2002. Karhunen-Loeve (KL) Filtering of Seismic Data Algorithm. Signal Analysis and Imaging Group (SAIL), Department of Physics, Alberta University.

* Scales, J. and Snider, R., 1998. [??] What is Noise?: Geophysics, Vol 63, p. 1122-1124.

* Yilmaz, O., 1987. Seismic Data Processing: Society of Exploration Geophysicists, Tulsa, p. 9-79.

* Zywicki, D., 1999. Traditional Seismic Surface Wave Tests. Georgia Tech Research Institute (GTRI), p. 11.34.

Ernesto Gomez Londono (1), Luis Castillo Lopez (2) and Thais de Souza Kazmierczak (2)

(1) Ingeominas Bogota, Recursos del Subsuelo, e-mail:

(2) Universidade Federal do Rio Grande do Sul-Brazil. e-mail:, e-mail:
COPYRIGHT 2005 Universidad Nacional de Colombia, Departamento de Geociencias
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2005 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Londono, Ernesto Gomez; Lopez, Luis Castillo; Kazmierczak, Thais de Souza
Publication:Earth Sciences Research Journal
Geographic Code:3COLO
Date:Dec 1, 2005
Previous Article:Wave equation datuming to correct topography effect on foothill seismic data.
Next Article:Application of gravimetry and electric tomography methods to obtain stratigrafic profiles: case study at University of Quindio and Puerto Espejo...

Related Articles
Multiple-valued logic; proceedings.
High-Fidelity Multichannel Audio Coding, 2d ed.
Elastic approach to remove surficial waves from shot gathers.
Evaluation of the T-wave alternans detection methods: a simulation study.
Evaluation of the T-wave alternans detection methods: a simulation study.

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters