Figure Captions

Methods

A previous Geophysical Corner article (Search and Discovery Article # 40454) showed how low frequency components (specifically that part of the data < 16 Hz) had a higher signal-to-noise ratio. Over the last decade or so, spectral decomposition has become a well-established tool that helps in the analysis of subtle stratigraphic plays and fractured reservoirs.

As the name suggests, spectral decomposition decomposes the seismic data into individual frequency components that fall within the measured seismic bandwidth, so that the same subsurface geology can be seen at different frequencies. Thick beds or features will be tuned and have relatively higher amplitude at lower frequencies, while thin beds will be tuned and have relatively higher amplitude at higher frequencies.

Spectral magnitude highlights features that are tuned, and spectral phase components enhance subtle fault and channel edges that can be used as input to subsequent seismic attribute analysis, such as coherence. Spectral decomposition is done by transforming the seismic data from the time domain into the frequency domain; this can be done simply by using the discrete Fourier transform.

There are other methods that also could be used for the purpose, such as:

Each of these methods has its own applicability and limitations, and the choice of a particular method also could depend on the end objective. For example:

The S-transform method is better than the continuous wavelet transform method, as it yields good temporal and spectral resolution. The matching pursuit method does not need any windowing and so yields both good temporal and spectral resolution. It is, however, computationally more expensive.

There are a number of commercial or proprietary implementations of spectral decomposition that are routinely used in the industry and are based on some variation of the above methods. Using any of the above spectral decomposition methods, the input seismic data volume can be decomposed into amplitude and phase volumes at discrete frequencies within the bandwidth of the data. These discrete frequency volumes are sometimes also referred to as common frequency volumes.

Example

Here, we illustrate the S-transform application of the spectral decomposition method to a case study from western Canada. In Figure 1 we show a comparison of stratal slices through the seismic data (Figure 1a) and 40 Hz spectral decomposition volumes (Figure 1b). The stratal slices were chosen 16 ms below a marker seismic reflector close to 950 ms on seismic data processed with 5-D interpolation used to regularize offsets and azimuths. By design, the original broadband seismic data volume can be reconstructed by a weighted average of the individual component volumes, including the 40 Hz volume shown in Figure 1b.

Note there is greater lateral variation in seismic amplitude in Figure 1b, which in this case is directly related to tuning effects associated with a distributary channel. The comparison of the stratal slices just 4 ms below these slices is shown in Figure 2, where a channel system is seen clearly on the 40 Hz spectral slice (Figure 2b), which is not as clearly seen on the seismic amplitude slice (Figure 2a). Finally, the coherence attribute comparison on the two volumes is shown in Figure 3, where the definition of the individual channels is seen very clearly on the spectral slice.

Conclusion

Spectral decomposition is an effective way of analyzing the seismic response of stratigraphic geologic features. Because of tuning and the variation of the signal-to-noise ratio with frequency, alternative spectral components can provide significant insight into the stratigraphic interpretation.

Return to top.