Seismic inverse Q filtering is a data processing technology for enhancing the resolution of reflection seismology images. Q is the anelastic attenuation factor or the seismic quality factor, a measure of the energy loss as the seismic wave moves.
Basics Seismic inverse Q-filtering employs a wave propagation reversal procedure that compensates for energy absorption and corrects wavelet distortion due to velocity dispersion. By compensating for amplitude attenuation with a model of the visco-elastic attenuation model type, seismic data can provide true relative-amplitude information for amplitude inversion and subsequent reservoir characterization. By correcting the phase distortion due to velocity dispersion, seismic data with enhanced vertical resolution can yield correct timings for lithological identification. However, Wang's outline of the subject is excellent and to follow his path, inverse Q filtering can be introduced based on the 1-D one-way propagation wave equation. He introduce this equation:.
d U ( r , w ) d r − i k ( w ) U ( r , w ) = 0 ( 1.1 ) {\displaystyle {\frac {dU(r,w)}{dr}}-ik(w)U(r,w)=0\quad (1.1)}
where U(r,w) is the plane wave of radial frequency w at travel distance r, k(w) is the wavenumber and i is the imaginary unit. Reflection seismograms record the reflection wave along the propagation path r from the source to reflector and back to the surface. With this approach Wang assumes that the plane wave U(r,w) has already been attenuated by a Q filter through travel distance r. We must have this in mind when we go to the step of finding a solution of (1.1). It is necessary that the initial U(r,w) either is already created by a forward synthetic Q-filtering process or taken directly from seismic surface data. Wang has introduced this concept in chapter 5 in his book. I think it is necessary to have this in mind also when the inverse theory is developed. Equation (1.1) has an analytical solution given by
U ( r + △ r , w ) = U ( r , w ) exp ( i k ( w ) △ r ) ( 1.2 ) {\displaystyle U(r+\bigtriangleup r,w)=U(r,w)\exp(ik(w)\bigtriangleup r)\quad (1.2)}
Kolsky's attenuation-dispersion model The wavenumber k(w) is an important variable in the solution (1.2). To obtain a solution that can be applied to seismic k(w) must be connected to a function that represent the way U(r,w) propagates in the seismic media. This functions can be regarded as a Q-model. The Kolsky Model is used extensively in seismic inverse Q-filtering. The model assumes the attenuation α(w) to be strictly linear with frequency over the range of measurement:
α = | w | ( 2 c r Q r ) ( 1.3 ) {\displaystyle \alpha ={\frac {|w|}{(2c_{r}Q_{r})}}\quad (1.3)}
And defines the phase velocity as:
1 c ( w ) = 1 c r ( 1 − 1 π Q r l n | w w r | ) ( 1.4 ) {\displaystyle {\frac {1}{c(w)}}={\frac {1}{c_{r}}}(1-{\frac {1}{\pi Q_{r}}}ln|{\frac {w}{w_{r}}}|)\quad (1.4)}
Where cr and Qr are the phase velocity and the Q value at a reference frequency wr. For a large value of Qr >>1 the solution (1.4) can be approximated to
1 c ( w ) = 1 c r | w w r | − γ ( 1.5 ) {\displaystyle {\frac {1}{c(w)}}={\frac {1}{c_{r}}}|{\frac {w}{w_{r}}}|^{-\gamma }\quad (1.5)}
where
γ = ( π Q r ) − 1 {\displaystyle \gamma =(\pi Q_{r})^{-1}}
… excerpt ends here. Continue reading the full article.


