9512.net
甜梦文库
当前位置:首页 >> >>

Scale Attributes from Continuous Wavelet Transform


INT 3.1

Scale Attributes from Continuous Wavelet Transform
Satish K. Sinha*, School of Geology and Geophysics, University of Oklahoma, Norman, OK; Partha S. Routh, Dept. of Geosciences, Boise State University, Boise, ID; Phil D. Anno, Seismic Imaging and Prediction, ConocoPhillips., Houston, TX; John P. Castagna, University of Houston, Houston, TX.
Summary Average instantaneous attributes of time-frequency decompositions are useful in revealing the time varying spectral properties of seismic data. In the continuous wavelet transform (CWT), a time signal is decomposed into a time-scale spectrum or a scalogram; unlike a timefrequency spectrum or a spectrogram from the short time Fourier transform (STFT). Although there are various approaches of converting a time-scale spectrum into a timefrequency spectrum we introduce new mathematical formulas to calculate spectral attributes from the scalogram. In this process, we bypass the conversion of a scalogram into a time-frequency spectrum and provide average spectral attributes based on scale. The attributes are: center frequency, dominant frequency, and spectral bandwidth. Since these attributes are based on the CWT, computation of these attributes avoids subjective choice of a window length. Introduction Spectral decomposition of non-stationary signals, like seismic signals, is conventionally achieved by the STFT. Therefore, average instantaneous spectral attributes, such as center frequency, dominant frequency etc., from a spectrogram inherits the properties of the STFT. A relatively new approach of spectral decomposition based on the CWT avoids the subjective choice of a window length and produces a time-scale spectrum also called scalogram. Sinha et al. (2003) converted a scalogram into a timefrequency spectrum called the time-frequency from CWT (TFCWT). One can use the TFCWT to compute instantaneous spectral attributes using the formulas outlined by Barnes (1993). However, we define new formulas to compute the instantaneous spectral attributes directly from a scalogram. In this definition we utilize the fact that the dilating support of the given wavelet, i.e. scale, is inversely proportional to the center frequency of the wavelet (Abry et al., 1993). We use Morlet wavelet for which this mapping is a good approximation. Theory Instantaneous spectral attributes of center frequency, dominant frequency, and spectral bandwidth are defined as various moments of a time-frequency distribution using familiar definitions from probability theory (Barnes, 1993). Mean of an instantaneous power spectrum is called instantaneous center frequency and is defined as the first moment of the time-frequency power spectrum along the frequency axis. Mathematically, it is given by,

f W ( f , t ) df f C (t )
0

2

,

--- (1)

W ( f , t ) df
0

2

where

f C (t ) is instantaneous center frequency. W ( f , t )

is the time-frequency map that can be obtained from STFT, TFCWT or any other spectral decomposition. Square root of the second moment of the time-frequency power spectrum along the frequency axis is called instantaneous dominant frequency.

f 2 W ( f , ) df
2 fD ( ) 0

2

,

--- (2)

W ( f , ) df
0

2

where

f D (t ) is instantaneous dominant frequency.

Standard deviation about the instantaneous center frequency is called instantaneous spectral bandwidth.

(f
2 f BW ( ) 0

f C ( )) 2 W ( f , ) df
, --- (3)

2

W ( f , ) df
0

2

where

f BW (t ) instantaneous spectral bandwidth. All the

above equations with different notations have been adopted from the paper by Barnes (1993) though competing equations for all of them exist in the literature. For a symmetrical wavelet, like Morlet wavelet (Torrence and Compo, 1998), center frequency of the wavelet at each scale can be assumed to be inversely proportional to the scale. Mathematically, it can be written as,

fC

k

,

--- (4)

SEG/Houston 2005 Annual Meeting

779

INT 3.1

Scale Attributes from Continuous Wavelet Transform

where

fC ,

, and

k are center frequency, scale and a

wavelet dependent constant respectively. The center frequency can be considered as a representative frequency for a particular scale. Differentiating equation 4, and replacing f C with f , we obtain

(equations 6, 7, and 8) are computationally efficient at the cost of slight smoother representation of the information. This efficiency translates to saving considerable amount of time depending upon length and number of traces in a large 3D data set. The next logical step is to compare results on a real seismic section and examine if there are any difference in interpretation. Instantaneous dominant frequency from TFCWT using equation 2 is shown in Figure 2. Instantaneous dominant frequency from scalogram using equation 5 is displayed in Figure 3. In both Figures (2 and 3) the frequency range on the colorbar is kept the same. Figure 2 and 3 indicates remarkable similarity and virtually can be considered as identical. Difference in timefrequency resolution is minimal. It is important to note that the computation of Figure 3 is about 5 times faster than Figure 2. Therefore, scale based attributes are more efficient from practical point of view. Figure 4 is produced using the STFT with a “Hanning” window of 200 ms window length.

df

k
2

d .

--- (5)

Changing from frequency to scale in the equations 1, 2, and 3 and assuming the wavelet dependent constant k 1 , we obtain

1 ~ f C (t )
3 0

FW ( , t ) d
, --- (6)

2

1
2 0

FW ( , t ) d FW ( , t ) d
, --- (7)
2

2

1 ~2 f D (t )
and
4 0

1
2 0

FW ( , t ) d
15

2

(a) Seismic Trace

(1 ~2 f BW ( )
0

~ 2 fC )
4

10

Amplitude

5 0 -5 -10 -15

FW ( , t ) d
, --- (8)
2

2

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

(b) Center Frequency
90

2 0

FW ( , t ) d

Frequency (Hz)

1

80 70 60 50 40 30 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

~ ~ ~ where f C , f D , and f BW are instantaneous center
frequency, instantaneous dominant frequency and instantaneous bandwidth respectively from the scalogram FW ( , t ) . Thus equation 6, 7 and 8 allows us to compute these spectral attributes directly from a scale map. Example Comparisons of instantaneous spectral attributes based on frequency and scale for a seismic trace (Figure 1a) are shown in Figures 1-b, c, and d. We note the remarkable similarity between these attributes although the attributes computed from the scalogram seems to be smoother. This is expected because a scale represents a band of frequency and not a single frequency. Therefore, overlap of frequencies from one scale to another produces smoother attributes, but the overall agreement is excellent. The favorable comparison is a consequence of the scale to center frequency mapping for a Morlet wavelet. Computation of attributes based on scalogram is much faster than those based on TFCWT. Therefore, from practical point of view, attribute definitions based on scale

From TFCWT

From Scalogram

90

(c) Dominant Frequency

Ti

( )

Frequency (Hz)

80 70 60 50 40 30 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

From TFCWT

From Scalogram

(d) Bandwidth
40

Bandwidth (Hz)

35 30 25 20 15 10 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

From TFCWT

From Scalogram

0.8

0.9

1

Time (s)

Figure 1: a) A seismic trace of length 1sec sampled at 4ms. It was used to compute frequency and scale based attributes. b) Instantaneous center frequency attribute, c) instantaneous dominant frequency attribute, and d) instantaneous spectral bandwidth computed from TFCWT (shown in black) and from scalogram (shown in red).

SEG/Houston 2005 Annual Meeting

780

INT 3.1

Scale Attributes from Continuous Wavelet Transform

Conclusions Scalogram obtained from continuous wavelet transform can provide useful average measures that are directly interpretable in terms of time-frequency attributes. We obtain new expressions for commonly used attributes such as dominant frequency, center frequency and bandwidth directly from a scalogram. These attributes can be efficiently calculated without converting a scalogram into a time-frequency spectrum. Field example presented in this paper shows a factor of 5 computational efficiency using scale based computation. In addition these attributes avoids the subjective choice of a window length commonly used in STFT. This is because the wavelet acts as an effective window via dilatation process.

Figure 4: Instantaneous dominant frequency from STFT using 200 ms window length.

Acknowledgements The authors acknowledge the Conoco research scientists for their valuable inputs. We would also like to thank Conoco Inc. (now ConocoPhillips) for providing the field dataset. We thank Dr. Dale Cox of Conoco for many useful discussions and for his help in the practical implementation in SU. References
Figure 2: Instantaneous dominant frequency from TFCWT.

Abry, P., Goncalves, P., and Flandrin, P., 1993, Waveletbased spectral analysis of 1/f processes: IEEE Internat. Conf. Acoustic, Speech and Signal Processing, 3, 237-240. Barnes, A. E., 1993, Instantaneous spectral bandwidth and dominant frequency with applications to seismic reflection data: Geophysics, 58, 3, 419-428. Sinha, S., Routh, P., Anno, P. and Castagna, J., 2003, Time-frequency attribute of seismic data using continuous wavelet transform, 73rd Ann. Internat. Mtg.: Soc. of Expl. Geophys., Expanded Abstracts, 1481-1484. Torrence, C., and Compo, G., P., 1998, A practical guide to wavelet analysis: Bulletin of the American Meteorological Society, 79, 1, 61-78.

Figure 3: Instantaneous dominant frequency from scalogram.

SEG/Houston 2005 Annual Meeting

781

EDITED REFERENCES Note: This reference list is a copy-edited version of the reference list submitted by the author. Reference lists for the 2005 SEG Technical Program Expanded Abstracts have been copy edited so that references provided with the online metadata for each paper will achieve a high degree of linking to cited sources that appear on the Web. Scale Attributes from Continuous Wavelet Transform REFERENCES Abry, P., P. Goncalves, and P. Flandrin, 1993, Wavelet-based spectral analysis of 1/f processes: IEEE International Conference on Acoustic, Speech and Signal Processing, 3, 237-240. Barnes, A. E., 1993, Instantaneous spectral bandwidth and dominant frequency with applications to seismic reflection data: Geophysics, 58, 419-428. Sinha, S., P. Routh, P. Anno, and J. Castagna, 2003, Time-frequency attribute of seismic data using continuous wavelet transform: 73rd Annual International Meeting, SEG, Expanded Abstracts, 1481-1484. Torrence, C., and G. P. Compo, 1998, A practical guide to wavelet analysis: Bulletin of the American Meteorological Society, 79, 61-78.


赞助商链接

更多相关文章:
更多相关标签:

All rights reserved Powered by 甜梦文库 9512.net

copyright ©right 2010-2021。
甜梦文库内容来自网络,如有侵犯请联系客服。zhit325@126.com|网站地图