The Kolsky basic model and modified model for attenuation and dispersion

The Kolsky basic model and modified model for attenuation and dispersion is the mathematical Q models that is most used in seismic applications. The basic Kolsky model is used for its simplicity in seismic data processing, however it does not rigorously satisfy the Minimum phase criterion and cannot satisfy the Kramers–Kronig relations. But then the Kolsky's modified model comes to our rescue, producing an accurate representation of the velocity dispersion within the seismic frequency band. The basic Kolsky model is presented in Kolsky's book "Stress waves in solids" that are available as a Google book."Stress waves in Solids".(1963 edition) Kolsky's modified model is presented in Wang's book Seismic inverse Q filtering. This book is also available as a Google book:Seismic inverse Q filtering (2008),

Basic

The theoretical background for mathematical Q models can be found in the Wikipedia article: Mathematical Q models. Here we found a function K(w) we can call a propagation constant in line with Futterman.

$K(iw)=k(w)+ia(w)\quad (1.3)$

k(w) can be linked to the phase velocity of the seismic wave with the formula:

$c(w)={\frac {w}{k(w)}}\quad (1.4)$

To obtain a solution that can be applied to seismic k(w) must be connected to a function that represent the way the seismic wave propagates in the seismic media. This functions can be regarded as a Q-model.

In his outline Wang calls the Kolsky-Futterman model the Kolsky model. The model assumes the attenuation α(w) to be strictly linear with frequency over the range of measurement:

$\alpha ={\frac {|w|}{(2c_{r}Q_{r})}}\quad (1.5)$

And defines the phase velocity as:

${\frac {1}{c(w)}}={\frac {1}{c_{r}}}(1-{\frac {1}{\pi Q_{r}}}ln|{\frac {w}{w_{r}}}|)\quad (1.6)$

Where c_r and Q_r are the phase velocity and the Q value at a reference frequency w_r.

For a large value of Qr >>1 the solution (1.6) can be approximated to

${\frac {1}{c(w)}}={\frac {1}{c_{r}}}|{\frac {w}{w_{r}}}|^{-\gamma }\quad (1.7)$

where

$\gamma =(\pi Q_{r})^{-1}$

Kolsky’s model was derived from and fitted well with experimental observations. A requirement in the theory for materials satisfying the linear attenuation assumption is that the reference frequency w_r is a finite (arbitrarily small but nonzero) cut-off on the absorption. According to Kolsky, we are free to choose w_r following the phenomenological criterion that it be small compared with the lowest measured frequency w in the frequency band. Those who want a deeper insight into this concept can go to Futterman (1962)

Computations

Bjørn Ursin and Tommy Toverud published an article where they compared different Q models. They used the Kolsky model as a reference model.

For each of the Q models Ursin B. and Toverud T. presented in their article they computed the attenuation (1.5) and phase velocity (1.6) in the frequency band 0–300 Hz. Fig.1. presents the graph for the Kolsky model - attenuation (left) and phase velocity (right) with c_r=2000 m/s, Q_r=100 and w_r=2π100 Hz.

• 
  Fig.1.Attenuation - dispersion Kolsky model

If we change the value for w_r to a much lower value 2π0.01 Hz, we will get a higher phase velocity for all frequencies:

• 
  Fig.2.Dispersion Kolsky model reference frequency w_r=2π0.01 Hz(green) w_r=2π100 Hz(blue)

Modification to the Kolsky model

The choice of w_r as the lowest frequency in the frequency band will introduce phase errors when we use the Kolsky model as an inverse Q filter. This is very well documented in Wang (2008). So the phase velocity formula in the basic Kolsky model is modified by using the highest frequency w_h as a reference. It could very well be the same as was used by Bjørn Ursin and Tommy Toverud above, w_h=2π100. Hz Then we can get a correct solution with inverse Q filtering with the Kolsky model.

References

• Wang, Yanghua (2008). Seismic inverse Q filtering. Blackwell Pub. ISBN 978-1-4051-8540-0.
• Kolsky, Herbert (1963). Stress Waves in Solids. Courier Dover Publications. ISBN 9780486495347.