![](images/pdficon.gif)
GCSpectral Decomposition for 3-D Geomodeling*
Brad C. Wallet¹
Search and Discovery Article #41408 (2014)
Posted August 11, 2014
*Adapted from the Geophysical Corner column, prepared by the authors, in AAPG Explorer, August, 2014. Editor of Geophysical Corner is Satinder Chopra ([email protected]). Managing Editor of AAPG Explorer is Vern Stefanic.
¹University of Oklahoma, Norman, Oklahoma ([email protected])
Spectral decomposition attributes use a localized time-frequency decomposition within individual traces to characterize seismic data by its spectral content. Since localized spectral content can be linked to physical phenomena of channels such as fill velocity and thickness, such techniques have been popular in the past decade for characterizing and interpreting of fluvial systems including turbidites.
The original work using spectral decomposition focused upon qualitative interpretation upon time slices and/or phantom horizons with a goal of broad scale interpretation and reservoir description. In recent years, however, a number of workflows have been developed to use spectral decomposition for 3-D geomodeling to extract architectural elements such as channels. Most of these methods are proprietary – and, to my knowledge, few have been documented or published.
A typical workflow might involve the use of the instantaneous frequency attribute to determine the dominant frequency associated with a channel. Once this information has been ascertained, an interpreter then might use some spatially associated thresholding (i.e. amplitude clipping) method to extract the channel. These results then would undergo further processing, such as smoothing or denoising, to produce a final interpretation.
♦General statement ♦Figures ♦Method ♦Example ♦Conclusion ♦Acknowledgment
♦General statement ♦Figures ♦Method ♦Example ♦Conclusion ♦Acknowledgment
♦General statement ♦Figures ♦Method ♦Example ♦Conclusion ♦Acknowledgment
♦General statement ♦Figures ♦Method ♦Example ♦Conclusion ♦Acknowledgment
♦General statement ♦Figures ♦Method ♦Example ♦Conclusion ♦Acknowledgment |
Figure 1 illustrates how a workflow like this would be constructed: An examination of these images suggests that the channels might be associated with frequencies in the 40-60 Hz range. This relationship is much stronger for channels "a" and "c," which have better defined tops than channel "b." Note that 50 Hz corresponds to a period of 20 ms, which is the approximate trough to trough thickness associated with all three channels. One possible pitfall is related to the spatial uncertainty associated with spectral attributes. Specifically, since spectral content is calculated using a vertical window, the calculated spectral information might be too imprecisely positioned to allow for good modeling in a vertical sense. This is further compounded by the inclusion of spectral phase, which might allow for the fitting of a rotated wavelet that is not properly centered upon the architectural element of interest. The net result of these problems could be an extracted model that looks geological but that is not properly located vertically. Better understanding of the vertical localization of spectral decomposition attributes is necessary to better understand if this is a true concern. To assess whether a 3-D geomodeling workflow was feasible and valid, I decided to look at the vertical characteristics of some spectral magnitude attributes for my channels of interest. My primary question was, "Are these channels completely and properly isolated vertically using spectral decomposition attributes?" One aspect of the vertical localization will depend upon the algorithm, as there are a number of methods to calculate spectral decomposition attributes. For the purposes of my study, I chose to use a matching pursuit (MP) algorithm, which is regarded as having superior localization. MP spectral decomposition works by decomposing a seismic trace into a set of time-shifted wavelets selected with replacement from a library of precomputed wavelets. The decomposition is done iteratively by fitting component wavelets and calculating residuals until the residuals are sufficiently small. Figure 2 shows some spectral magnitude attributes as calculated using a MP spectral decomposition algorithm calculated using a library of Ricker wavelets. Note that the trough-to-peak distance for all three channels is approximately 30 ms, while the trough-to-trough distance in these channels, as previously noted, is approximately 20 ms. These correspond to frequencies of approximately 30 Hz and 50 Hz, respectively. A low magnitude region matches our intuition for a peak-to-trough region when decomposed using a unimodal wavelet with an internal peak and trough as the internal components will destructively combine when convolved with the two positive sides of the wavelet. Channel "c," while well defined, is not delineated to its proper top. Again, channel "c" is particularly well defined. However, there are some less well defined regions along a portion of the top of channel "a." Either channel b's top or bottom would be somewhat misplaced spatially, depending upon the choice of cutoff for segmenting. Additionally, the inter-fluvial region between channels b and c is sharply defined by a particularly strong low magnitude event. In conclusion, these results show that there is a promise to algorithms using spectral decomposition attributes to do 3-D geomodeling of channels. In the case of the 50 Hz decomposition, the theoretical analysis of the wavelet is based upon the consideration of the bottom of the channel and an internal reflector, while the highlighted region appears to cover the full range of the channel. I believe this has to do with the phase of the fit wavelets. This is a conjecture that suggests future work to better understand the role of phase in modeling channels. Finally, while the results do show vertical localization, all of the images of the attributes are noisy. This emphasizes the importance of good post-processing workflows to provide parsimonious and geologically sound models. Additionally, as some boundaries are misplaced vertically, it is my belief that naïve algorithms will give less desirable results than one that incorporates some heuristical processing to ensure geologically and geophysically consistent models. I would like to thank Victor Aarre and AAPG members Kurt Marfurt and Roger Slatt (Honorary) for many insightful discussions concerning this and related topics. Additionally, I would like to thank Schlumberger for providing the data used in this work. Finally, I would like to thank the sponsors of the Attribute Assisted Seismic Processing and Interpretation (AASPI) Consortium for their ongoing support. |