Inverse Q filtering

From SEG Wiki
Jump to: navigation, search
Seismic Data Analysis
Series Investigations in Geophysics
Author Öz Yilmaz
ISBN ISBN 978-1-56080-094-1
Store SEG Online Store

Frequency attenuation caused by the intrinsic attenuation in rocks was discussed in entries dedicated to gain applications and the convolutional model in the time domain. Attenuation causes loss of high frequencies in the propgating waveform with increasing traveltime. This gives rise to a nonstationary behavior in the shape of the wavelets associated with reflection events at different times.

Wave attenuation usually is described by a dimensionless factor Q, which is defined by the ratio of the mean stored energy to the energy loss over a period of time that is equivalent to one cycle of a frequency component of the waveform [1]. Time-variant deconvolution and time-variant spectral whitening discussed in this section are processes that can correct for the time-varying effects of attenuation by spectral flattening. A deterministic alternative to compensate for frequency-dependent attenuation is provided by inverse Q filtering.

The amplitude spectrum of the inverse Q filter is given by (Appendix B.9)


where ω is the angular frequency component associated with the input trace and τ is the time variable associated with the output trace from inverse Q filtering. The phase spectrum of the inverse Q filter usually is assumed to be minimum-phase, which can be computed by taking the Hilbert transform of the amplitude spectrum (Appendix B.4).

Application of the inverse Q filter requires knowledge of the attenuation factor Q, which usually is assumed to be constant. A compilation of laboratory measurements of Q for some rock samples is given by Table 2-28.

Note from Table 2-28 that most measurements have been made at extremely high frequencies compared to the typical frequency band of seismic waves. Nevertheless, by assuming frequency-independent Q factor [1], these measurements can still be considered useful. Also note that the Q factor can vary significantly for limestone, sandstone, and shale rock samples of different origin.

The inverse of the amplitude spectrum defined by equation (48) can be used to obtain a quantitative measure of attenuation. In terms of frequency f, wave velocity v and depth z = , the inverse is


To determine how far in depth the wave has to travel before its amplitude reduces to, say, one-tenth of its amplitude at the surface z = 0, rewrite equation (49) as follows:

Table 2-28. Intrinsic attenuation measurements in rocks adapted from [2].
Rock Type Attenuation Constant, Q Frequency Range (Hz)
Basalt 550 3,000-4,000
Granite 300 20,000-200,000
Quartzite 400 3,000-4,000
Limestone I 200 10,000-15,000
Limestone II 50 2-120
Limestone III 650 4-18,000
Chalk 2 150
Sandstone I 25 550-4,000
Sandstone II 125 20,000
Sandstone III 75 2,500-5,000
Sandstone IV 100 2-40
Shale I 15 75-550
Shale II 75 3,300-12,800

Note that the smaller the Q factor, the lower the velocity and the higher the frequency, the shallower the depth at which the wave amplitude decays to a fraction of the wave amplitude at z = 0. Table 2-29 lists the z values for a frequency of 30 Hz and a velocity of 3000 m/s for a range of Q values.

Table 2-29. Depth at which wave amplitude drops to one-tenth of its original at the surface for a range of Q values (equation 50).
v = 3000 m/s
f = 30 Hz
Q Factor Depth in m
25 1,830
50 3,660
100 7,325
250 18,312
500 36,625

Note that the smaller the Q factor the shallower the depth at which the amplitude drops to the specified value of one-tenth of the original value at the surface z = 0. For very large Q values, that is, for small attenuation, the amplitude reduction to the specified value does not take place until the wave reaches very large depths beyond the exploration objectives.

Unfortunately, there is no reliable way to estimate the attenuation factor Q directly from seismic data. At best, inverse Q filtering can be applied to post- or prestack data (Appendix B.9) using a range of constant Q factors to create a Q panel, much like a filter panel (The 1-D Fourier transform). The factor that yields the flattest spectrum in combination with other signal processing applications — deconvolution and time-variant spectral whitening, is chosen as the optimum Q value.

See also


  1. 1.0 1.1 Kjartansson, 1979, Kjartansson, E., 1979, Constant Q-wave propagation and attenuation: J. Geophys. Res., 84, 4737–4748.
  2. Waters, 1981, Waters, K. H., 1981, Reflection seismology: Second edition, John Wiley & Sons.

External links

find literature about
Inverse Q filtering
SEG button search.png Datapages button.png GeoScienceWorld button.png OnePetro button.png Schlumberger button.png Google button.png AGI button.png